Trypanorhynch Assemblages Indicate Ecological and Phylogenetical Attributes of Their Elasmobranch Final Hosts

This study explores non-metric multidimensional scaling (nMDS) as a tool for investigating parasites as indicators of the elasmobranch biology. An attractive feature of nMDS is its ability to allow assemblage-level parasite data to be simultaneously applied to questions of host biology. This method was examined using the tapeworm order Trypanorhyncha Diesing, 1863, which is known to be transmitted among their hosts through the marine food web (via predation), can unambiguously be identified in the intermediate and final hosts, and has the potential as an indicator of the host feeding biology. Our analyses focused on trypanorhynch assemblages in elasmobranchs as definitive hosts. The relationships between trypanorhynch assemblages and the depth, feeding ecology, habitat, and phylogeny for all sharks were complex, but we found that depth distribution, diet composition and habitat type were the major influencing factors. Several species of sharks showed different characters than known from their descriptions that could be attributed to the change of shark behavior or the trypanorhynch host path. The relationship between the trypanorhynch assemblage and factors for carcharhiniform species alone was more robust than for all sharks. In the carcharhiniform analysis, the relationship between habitat type and trypanorhynch assemblage was most remarkable. Overlapping host ecology was evident even in phylogenetically-distant related hosts.


Introduction
Elasmobranchs are difficult to observe in nature and information on rare or less-frequently caught species is scarce.This especially concerns aspects of their ecology, main habitat, depth range, and most important prey items.Their prey is rapidly digested and affiliation, e.g., to mainly fish or invertebrate-feeding activities, is difficult.Fish parasites have been widely used as biological indicators for the host ecology, including their zoogeographical distribution, migration, and feeding behavior [1].However, the methods applied for the elasmobranchs are, so far, limited due to often restricted and unpredictable catches and less availability of specimens to study.
The cestode order Trypanorhyncha Diesing, 1863, is considered to be a basal tapeworm group of elasmobranchs.Its members can be characterized by a scolex with a tentacular apparatus that can be used for classification [2][3][4][5].The most recent revision of the trypanorhynchs subdivided the order into five different superfamilies and 15 families [6].The first large-scale molecular phylogenetic analyses of this group included trypanorhynch species belonging to all superfamilies and 12 of the 15 recognized families [7,8].The analysis of complete simple sequence repeats DNA (ssrDNA) and partial long sequence repeat DNA (lsrDNA) (D1-D3) revealed high resolution and the presence of a basal split, representing the shark-and ray-inhabiting lineages as new suborders.These results gave, in part, support for the trypanorhynch classification by Palm [6] as well as more insight into the intrafamily phylogenetic relationships within the Trypanorhyncha.The host-parasite-checklist of these cestodes indicated that some of the superfamilies and families seem to have host preferences [6], while the patterns of others can be explained by host switches, even between rajiform and selachiform elasmobranchs [7].
The ~300 described species of the order Trypanorhyncha so far exhibit the following general life-cycle pattern [6,[9][10][11].The first intermediate host becomes infected when it consumes an oncosphere or free-swimming coracidium larva.The first parasitic larval stage is a procercoid; the limited available data suggest that this stage parasitizes copepods.Second intermediate hosts include a wide array of marine animals, but invertebrates and teleosts figure particularly prominent at this level.Unlike most groups of tapeworms, the trypanorhynch final larval stage, though having a different morphology (plerocercoid, plerocercus, merocercoid), is generally identifiable to species because it exhibits the adult hook pattern, a pattern on which much of the taxonomy of the group is based [4][5][6].Some trypanorhynchs use a paratenic host following the final intermediate host; this host serves to bridge food-web gaps.As adults, trypanorhynchs parasitize the spiral intestine of sharks and rays (elasmobranchs).Fortunately, it is common for individual elasmobranch species to host multiple trypanorhynch species and, thus, assemblage-level data are readily available for many elasmobranchs that serve as definitive hosts.
Our study uses non-metric multidimensional scaling (nMDS) as a tool for investigating parasites as indicators for selected ecological attributes and the phylogeny of their elasmobranch final hosts.A particularly novel and attractive feature of nMDS is the possibility to allow assemblage-level parasite data to be simultaneously applied to questions of host biology.The method was examined using trypanorhynch tapeworm assemblages, which have obvious potential as indicators of host feeding biology in the elasmobranch.The complex life-cycles of these parasites, each of which involves a minimum of three different hosts and three distinct life-cycle stages that vary in their degree of host specificity [11] also suggest trypanorhynchs as appropriate indicators of other aspects of their host biology.
Among elasmobranchs, sharks (Selachii) were specifically chosen because they vary substantially in diet and habitat, which allows a wide spectrum of aspects of host biology to be explored.Within the Selachii, more detailed analyses were conducted on sharks of the order Carcharhiniformes (ground sharks).Comprehensive trypanorhynch data were available for this shark order, and provided a unique opportunity to perform nMDS analyses on a more precise taxonomic level.Factors concerning the shark biology included (1) taxonomy (order for the Selachii analyses, or family for the Carcharhiniform analyses); (2) habitat (benthic vs. pelagic); (3) depth distribution (above or below 200 m depth); and (4) diet (primarily invertebrates vs. primarily teleosts).Trypanorhynch assemblage dependence on these four factors was assessed by superimposing the different variables into the two-dimensional ordination plots resulting from nMDS analysis of the hierarchical agglomerative classification of shark species generated from their trypanorhynch assemblages.Trypanorhynch assemblages were considered to be appropriate as indicators of the biological factors that were found to be significantly correlated with the agglomerative classification.

All Sharks
The result of the hierarchical, agglomerative classification (CLUSTER) of trypanorhynch assemblages of 50 species of Selachii revealed five major clusters of shark species at an arbitrary dissimilarity level of 78% (Figure 1).The two-dimensional ordination plots resulting from nMDS of this agglomerative classification, with sharks coded for taxonomic order, diet, habitat, and depth distribution are shown in Figure 2. The analysis of similarities (ANOSIM) revealed significant differences in trypanorhynch assemblage composition between sharks that primarily consume invertebrates versus those that primarily consume teleosts (global R = 0.279; p < 0.01), between those from benthic and pelagic habitats (global R = 0.275; p < 0.01), and between sharks occurring in deep versus shallow water (global R = 0.327; p < 0.01), thus supporting the significance of the major clusters resulting from the classification and ordination methods.With respect to taxonomy, the ANOSIM showed significant differences among shark orders (global R = 0.381; p < 0.01), suggesting a relationship between trypanorhynch genera and shark orders.In ANOSIM pairwise comparison, the trypanorhych assemblages of the Carcharhiniformes were significantly different from those of the Squaliformes (global R = 0.546; p < 0.01), the Orectolobiformes (global R = 0.466; p < 0.1), the Lamniformes (global R = 0.312; p < 0.01), and the Hexanchiformes (global R = 0.275; p < 0.05).In addition, the trypanorhynch assemblages of the Lamniformes were significantly different from those of the Orectolobiformes (global R = 0.432; p < 0.01) and Squaliformes (global R = 0.355; p < 0.05), and the assemblages of the Orectolobiforms were significantly different from those of the Squaliformes (global R = 0.620; p < 0.01).No significant differences were found between the trypanorhynch assemblages of any of the other 16 pairs of shark orders.

Carcharhiniformes
The result of the hierarchical, agglomerative classification (CLUSTER) of trypanorhynch assemblages of Carcharhiniformes, the order with the most comprehensive trypanorhynch dataset in Palm [6], is shown in Figure 3. Two major clusters of sharks at an arbitrary dissimilarity level of 82% were recovered (I and II).The nMDS plots of the classification, with sharks coded for family, diet, and habitat are shown in Figure 4. Analysis of similarities revealed significant differences in trypanorhynch assemblage composition between sharks that primarily consume invertebrates versus those that primarily consume teleosts (global R = 0.683; p < 0.01), and also between those from benthic and pelagic habitats (global R = 0.590; p < 0.01).Regarding taxonomy, however, ANOSIM pairwise comparisons, detected significant differences on trypanorhynch assemblages only between Carcharhinidae and Triakidae (global R = 0.678; p < 0.01) and between Sphyrnidae and Triakidae (global R = 0.748; p < 0.01).No significant differences were found between the trypanorhynch assemblages of any of the other eight pairs of carcharhiniform families.

All Sharks
Our results demonstrate that depth distribution, diet composition and habitat type (benthic or pelagic) were the major influencing factors responsible for the composition of trypanorhynch assemblages and the groupings (cluster) of the sharks.Depth distribution affected the formation of all shark clusters (Figure 2): Cluster I included a single shallow/deep water species, Cluster II was dominated by deep water species, and Cluster III by shallow water species.Clusters IV and V were generated within Cluster III, and were dominated by invertebrate or teleost feeders, respectively.Cluster IV and Cluster V were also influenced by habitat type (benthic or pelagic), respectively.Moreover, Cluster II was dominated by piscivorous sharks.In general, no obvious pattern of the shark orders was found in the analysis.For instance, the most represented order Carcharhiniformes occurred in Clusters II and III (IV and V), and Lamniformes were found in Clusters I, II, and III (V).
Host-parasite associations are essential to gain information on the ecology of the hosts.The four shallow water shark species, Squatina squatina, Pristiophorus cirratus, Carcharias taurus, and Scyliorhinus stellaris, became an exception in Cluster II, because they were included in the deep water species group.That was probably due to the occurrence of Nybelinia, Grillotia, and Gilquinia that generally characterized the deep water species group, but which were also found in these four exception species.Thus, we can conclude that these four species which had been coded as shallow water types [12,13] might be collated to shallow/deep water shark species.Additionally, two shark species (Deania hystricosa and Centrophorus squamosus) that were part of Cluster II had been coded as unknown diet species.Grillotia that was found in both species characterized them as belonging into the fish feeder shark group.In Cluster III, Carcharhinus amblyrhynchoides was recognized as an unknown depth species.The trypanorhynchs Heteronybelinia, Callitetrarhynchus and Dasyrhynchus were found in this shark.Based on SIMPER analysis and Mann-Whitney U test, Callitetrarhynchus and Dasyrhynchus are characteristic trypanorhynchs of shallow water sharks.Hence, we consider C. amblyrhynchoides as a shallow water species.Three exception species were found in Cluster IV that was dominated by benthic and invertebratefeeding sharks, Galeorhinus galeus, Triakis semifasciata and Sutorectus tentaculatus.This exception leads us to suspect that the habitat type of G. galeus needs additional scrutiny; although Compagno [13] indicates this to be a pelagic species.Its trypanorhynch fauna suggests that it may be more benthic in its habits than currently recognized.From further analysis of the trypanorhynch assemblages we suggest that G. galeus, T. semifasciata, and S. tentaculatus prefer invertebrates over fish, although that had been coded as fish/invertebrate feeders.
In the last cluster V that was dominated by pelagic and piscivorous sharks, we found seven exception species.These were Hymenolepis microstoma, Mustelus manazo, Notorynchus cepedianus, Stegostoma fasciatum, Nebrius ferrugineus, Alopias superciliosus, and Carcharhinus melanopterus.C. melanopterus that had been coded as fish/invertebrate feeder seems to prefer more fish over invertebrates based on its trypanorhynch fauna.In A. superciliosus that was pelagic/benthic realm, we found Nybelinia, Mixonybelinia and Sphyriocephalus.Nybelinia, and Mixonybelinia characterized both groups, pelagic and benthic environment, while Sphyriocephalus characterized the pelagic environment only.Hence, we conclude that A. superciliosus, where adult species of Sphyriocephalus can be found, prefers the pelagic habitat.We also noted that adult Hepatoxylon was found in No. cepedianus [6] which suggests this shark to live in a pelagic/benthic habitat, although coded as benthic [12].Hepatoxylon is a sphyriocephalid trypanorhynch occurring in oceanic but also deep water hosts.Although Sphyriocephalus characterises pelagic species, the species number of Sphyriocephalus was not significant different between pelagic and benthic sharks.We found parasite genera that characterised the pelagic and benthic group in M. manazo, S. fasciatum, and N. ferrugineus.The parasite genera found in these three sharks also characterized fish and invertebrate feeder.We suspect the three shark species were changing from invertebrate feeder and benthic habitat to be fish/invertebrate feeder and occurring in the pelagic/benthic habitat.
The last exception species in cluster V was H. microstoma.Based on its trypanorhynch fauna, H. microstoma was a benthic and invertebrate feeder.It became part of the group that was dominated by pelagic and fish feeder species.A plausible reason for this was the genus Nybelinia that was found in H. microstoma.Nybelinia was found in pelagic and benthic habitat species, as well as fish and invertebrate feeder species.The species numbers of Nybelinia in pelagic species were significantly higher than the numbers in benthic species.However, this was not significantly different between fish and invertebrate feeders.Another genus found in H. microstoma was Trigonolobium (one species).The species numbers of Trigonolobium in invertebrate feeders was significantly higher than the numbers in fish feeders, but it was not significantly different between pelagic and benthic sharks.We conclude that based on H. microstoma, habitat was a more dominant factor than diet in generating the Clusters IV and V within Cluster III.
Several shark species, based on the trypanorhynch fauna, were differently characterized than known from the literature [12][13][14][15].An obvious reason for this is the poor knowledge on many shark species.Another explanation can be that the shark species can changed their behavior due to environmental changes (e.g., availability of fish, temperature increase, less production).Sharks might change their diet with regard to their size and prey diversity [16,17].It has been stated that sharks are one of the most affected species concerning habitat loss or habitat change caused by climate change [18], which could explain their habitat changes.Another factor causing difficult affiliation of the sharks is that the knowledge on the trypanorhynch fauna is also far from complete, with new host records and new species descriptions added year by year.In general, trypanorhynchs have low (second) intermediate host and high final host specificity [11].Although many trypanorhynchs have higher final host specificity, when they change the intermediate host (e.g., from invertebrate to fish or from fish to other fish that has a different habitat) they have a high possibility to reach further final host.This leads to regular host switches that have been recognized in trypanorhynch phylogeny.When this situation happens, then the trypanorhynch fauna makes it more difficult to characterize the shark species and results in conflicting evidence.

Carcharhiniformes
Among the factors investigated in the carcharhiniform analysis, the relationship between habitat and trypanorhynch assemblage was most remarkable.In fact, with the exception of G. galeus, the two major clusters of sharks recovered from the CLUSTER analysis of trypanorhynch assemblage data exactly corresponded to habitat occupied by the sharks.However, G. galeus coded as a pelagic shark was identified as a benthic shark (see Discussion above).
nMDS suggests that the trypanorhynch fauna has substantial potential to indicate the habit type of their elasmobranch final hosts.However, it is interesting to explore this in more depth by considering the individual genera comprising those assemblages.The trypanorhynch fauna of pelagic carcharhiniform sharks is primarily characterised by species of the genera Otobothrium, Dasyrhynchus, Heteronybelinia, Tentacularia, and Floriceps.In contrast, the trypanorhynch fauna of benthic carcharhiniform sharks is primarily characterized by species of Dollfusiella, Lacistorhynchus, Diesingium, and Eutetrarhynchus.All nine of these genera can be considered to be useful discriminatory taxa of environmental tendencies as the former five genera are abundant in pelagic sharks and scarce in benthic sharks while the latter four genera are scarce in pelagic sharks and abundant in benthic species.It should be noted that Nybelinia is abundant in sharks from both pelagic and benthic environments, as it occurs likewise in sharks and rays and, thus, it must be carefully used in combination with other taxa as a discriminatory taxon assigning the elasmobranch to a habit type.Thus a combination of tentaculariids, lacistorhynchids and otobothriids characterizes the Charcharhiniformes.According to Palm [6] these trypanorhynch taxa genera represent at least two different life-cycles.The intermediate hosts of tentaculariids, such as Tentacularia and Nybelinia, are pelagic or oceanic teleosts as well as squids as paratenic hosts.These trypanorhynchs develop robust plerocercoids in the intermediate host, and are characterised by low host specificity [11].Lacistorhynchoids and otobothroids develop plerocerci with blastocysts, and the intermediate hosts are often pelagic fish.For example, the otobothriid Otobothrium penetrans mainly infests belonids as second intermediate hosts that are most common in the upper layer (10 m) of the water column [19,20].Eutetrarhynchoids and gymnorhynchoids are generally missing within the pelagic carcharhiniforms.The former have life cycles that involve invertebrates as intermediate hosts, and are most commonly found on rajiform elasmobranchs, often in the benthic environment.The latter develop merocercoids and live in pelagic hosts.
The relationship between diet and trypanorhynch assemblage was similarly notable.Of the 28 species of carcharhiniform sharks included in the analysis, the diets of all but four species were coincident with the two major clusters of sharks recovered from the CLUSTER analysis of trypanorhynch assemblage data (Figure 4).Remarkably, these exceptions were the three species (C.melanopterus, T. semifasciata, and G. galeus) coded as consuming both fish and invertebrates based on the diet data available in Compagno [13].These results suggest that fish play a larger role in the diet of C. melanopterus, than do invertebrates.They further suggest that the diets of T. semifasciatea and G. galeus emphasize invertebrates over fish.
At the individual genus level, Otobothrium, Dasyrhynchus, Heteronybelinia, Callitetrarhynchus, Tentacularia, and Floriceps appear to be excellent indicators of a preference for fish as they are particularly abundant in piscivorous sharks and rare to non-existent in sharks that consume invertebrates.In contrast, Dollfusiella, Lacistorhynchus, Diesingium, and Trigonolobium appear to be particularly significant in invertebrate feeders, and rare in piscivorous shark species.According to food preferences, the SIMPER analysis showed a similar picture.This is in line with the typical life cycle of the trypanorhynchs (see above).According to the trypanorhynch fauna, the studied carcharhiniforms of group 2, Carcharhinidae, and Sphyrnidae, obviously, are mainly pelagic fish feeders.
Carcharhiniforms that were separately clustered from the pelagic fish feeders were the Triakidae, Scyliorhinidae and Hemigaleidae.The benthic carcharhiniforms of group 1 are characterized by the trypanorhynch genera Dollfusiella, Lacistorhynchus, Nybelinia and Diesingium, while the invertebrate feeders were characterized by Dollfusiella, Lacistorhynchus, Diesingium, and Trigonolobium.The life cycle of the eutetrarhynchoid Dollfusiella is known to involve benthic invertebrates such as shrimps and crabs as second intermediate host, while the second intermediate host of Trigonolobium is yet unknown [6].However, our results give a first indication that Trigonolobium might also involve marine invertebrates within its life cycle.The importance of Lacistorhynchus and Diesingium to distinguish this group of carcharhiniforms gives conflicting evidence.Both genera are typical trypanorhynchs from triakid sharks, with the plerocercus stage of the former often infecting a range of, mainly demersal or bottom-dwelling fish species.Thus, Lacistorhynchus and also the lacistorhynchid genus Diesingium might indicate the benthic environment of triakid carcharhiniform sharks, however, they cannot be utilized to analyze the feeding ecology of their hosts.
The two major clusters of sharks recovered from CLUSTER analysis were consistent with the carcharhiniform family.In no instance were species of the same family found in both clusters.In addition, the ANOSIM test showed that the trypanorhynch assemblage of the Carcharhinidae significantly differed from that of the Triakidae.The lack of significance between the Carcharhinidae and the Scyliorhinidae and the Hemigaleidae is likely to be the result of the lack of power of this test caused by the analyzed dataset, as the analysis included only a single species in each of the latter families.However, an additional level of resolution emerges if the phylogenetic relationships among carcharhiniform families will be considered.Based on the relatively robust, albeit preliminary molecular study Vélez-Zuazo and Agnarsson [21] concluded that hammerhead sharks (Sphyrnidae) were grouped with whaler sharks (Carcharhinidae) in phylogenetic relationship.The largest cluster (Cluster II) resulting from the agglomerative classification consists solely of members of this clade of sharks.In addition, ANOSIM failed to detect a significant difference between the assemblages of the Sphyrnidae and the Carcharhinidae.This test would be interesting to explore further with additional sphyrnid species because it included 16 species of carcharhinid sharks, but only three of the nine species of Sphyrnidae.

Materials and Methods
To study the host-parasite relationships between the Trypanorhyncha and their elasmobranch hosts, trypanorhynch data were analyzed from ca. 4000 host/parasite records compiled by Palm [6].Only those shark species for which at least two adult trypanorhynch records had been available, were included in the analyses (See Table S1).Based on this minimum criterion, the trypanorhynch data consisted of 307 host/parasite records for a total of 50 shark species representing seven of the eight selachian orders (Carcharhiniformes, Lamniformes, Hexanchiformes, Orectilobiformes, Pristiophoriformes, Squaliformes, and Squatiniformes).These data, however, conspicuously emphasized carcharhiniform sharks, with 207 of the host/parasite records belonging to 28 species of this most speciose shark order.
A matrix was generated by scoring each shark species according to the total number of species in each trypanorhynch genus reported from any of the shark species.The absence of a trypanorhynch genus from a shark species was scored 0 in the matrix.A resemblance matrix among shark species was generated from this matrix of trypanorhynch assemblages choosing the Bray-Curtis measure [22,23].All further multivariate data analyses were performed using the PRIMER software (release 6, Primer-E Ltd. 6.1.11,Meadow View, Plymouth, UK; Clarke and Warwick 2001) [24].Sharks were classified by hierarchical agglomerative cluster analysis of this resemblance matrix using the group-average linking method CLUSTER [23].The results of these analyses were then ordinated using nMDS two-dimensional ordination plots of shark species which were generated from trypanorhynch assemblage data for each of the factors under consideration.These plots were produced from both the full matrix of data for the selachian analyses, and also from the more restricted matrix consisting exclusively of carcharhiniform shark data.In each case, the conditions of the factor under consideration were superimposed on the ordination plot for that factor by coding each shark species with a symbol representing the condition it exhibited.For each factor, ANOSIM [25,26] was performed to test significant differences among conditions of the factor as reflected by trypanorhynch assemblage.In instances in which significant differences were seen, SIMPER analyses [27] were performed to explore which trypanorhynch genera contributed most to these differences.In these analyses trypanorhynch genera were ordered according to decreasing contribution (%) to the total similarity observed.Genera accounting for high levels of similarity were determined to be good candidates as individual indicators.
For the analyses, each shark species was categorized for habitat (benthic or pelagic) and diet (primarily invertebrates and primarily teleosts).For the analysis of the Selachii, sharks were also coded for depth of distribution (below 200 m as deepwater, 200 m or less as shallow water species).In the carcharhiniform analyses taxonomy was assessed at the level of family; in the Selachii analyses taxonomy was assessed at the level of order?Sharks were categorized for these factors based on information published in the work of Compagno [12][13][14][15], because this provides detailed data for all of the target species.Shark taxonomic position at order and family levels was also based on Compagno [14].Sharks were categorised for diet, habitat, and depth distribution for the Hexanchiformes, Squaliformes, and Squatiniformes [11], for the Carcharhiniformes [12], and for the Lamniformes, Orectilobiformes, and Pristiophoriformes [14].Categorization by habitat type was fairly straightforward as most shark species were either benthic or pelagic.Coding for depth was similarly unambiguous with the recorded depths of most sharks being either above, or below 200 m.Categorisation of sharks according to diet, however, was somewhat more ambiguous.While many sharks had a clear preference for either invertebrate or teleost prey, other species consumed both types of prey with some regularity.In such instances, sharks were coded as exhibiting both conditions of this factor.For the purpose of this study, two mutually exclusive conditions for this factor were chosen, based on the final intermediate hosts utilized by trypanorhynchs (i.e., invertebrates vs. teleosts).
Kruskal-Wallis test and Mann-Whitney U test were performed to compare the number of parasite species from each genus in each shark group [28,29].The Kruskal-Wallis test was employed for more than two categories of sharks (e.g., selachian order; Carcharhiniformes, Lamniformes, Hexanchiformes, Orectilobiformes, Pristiophoriformes, Squaliformes, and Squatiniformes), and Mann-Whitney U test was employed for two categories of sharks (e.g., habitat; benthic and pelagic).

Conclusions
This study demonstrates that nMDS of the elasmobranch parasitic Trypanorhyncha Diesing, 1863, is a useful tool for investigating parasites as indicators of the host biology in marine ecosystems.Transmission through the marine food web (via predation) and an unambiguous identification in the final (sharks and rays) and intermediate (teleosts, other marine invertebrates) hosts allow conclusions to be made on the feeding biology of the host.The relationship between trypanorhynch assemblage and investigated factors for all sharks was complex, but we found that depth distribution, diet, and habitat type were the major influencing factors.The relationships within carcharhiniform species were more robust, especially if the factor habitat is considered.Overlapping host ecology was evident even in phylogenetically-distant hosts, demonstrating host switches and low host specificity in some of the trypanorhynch taxa.Considering that the underlying dataset is based on a host-parasite list from 2004 [6] and more species descriptions and host records are added continuously nMDS proves to be a powerful tool to attribute aspects of shark biology (phylogeny, habitat, depth distribution, diet) to rare and difficult-to-observe elasmobranch species in future.

Supplementary Materials:
The following are available online at www.mdpi.com/2410-3888/2/2/8/s1,Table S1: List of shark species analysed in the recent study; Table S2: SIMPER analysis results for similarities in each habitat of all sharks; Table S3: SIMPER analysis results for similarities in each diet of all sharks; Table S4: SIMPER analysis results for similarities in each depth of all sharks; Table S5: SIMPER analysis results for similarities in each order of all sharks; Table S6: SIMPER analysis results for similarities in each habitat of carcharhiniform; Table S7: SIMPER analysis results for similarities in diet of carcharhiniform; Table S8: SIMPER analysis results for similarities in each family of carcharhiniform.

Figure 1 .
Figure 1.Cluster analysis of all sharks according to Trypanorhyncha assemblages.Roman numerals I-V indicate the five major clusters.

Figure 1 .
Figure 1.Cluster analysis of all sharks according to Trypanorhyncha assemblages.Roman numerals I-V indicate the five major clusters.

Figure 2 .
Figure 2. Multidimensional scaling plot for all shark species based on Bray-Curtis similarity and grouped by depth (A), habitat (B), diet (C), and order (D).Clusters I-III marked in blue, cluster IV in red, and cluster V in green.

Figure 2 .
Figure 2. Multidimensional scaling plot for all shark species based on Bray-Curtis similarity and grouped by depth (A), habitat (B), diet (C), and order (D).Clusters I-III marked in blue, cluster IV in red, and cluster V in green.

Figure 3 .Figure 3 .
Figure 3. Cluster analysis of carcharhiniform sharks according to Trypanorhyncha assemblages.Roman letters I-II indicate the two major clusters.

Figure 4 .
Figure 4. Multidimensional scaling plot for carcharhiniform species based on Bray-Curtis similarity and grouped by habitat (A), diet (B), and family (C).Clusters I and II marked in blue.Figure 4. Multidimensional scaling plot for carcharhiniform species based on Bray-Curtis similarity and grouped by habitat (A), diet (B), and family (C).Clusters I and II marked in blue.

Figure 4 .
Figure 4. Multidimensional scaling plot for carcharhiniform species based on Bray-Curtis similarity and grouped by habitat (A), diet (B), and family (C).Clusters I and II marked in blue.Figure 4. Multidimensional scaling plot for carcharhiniform species based on Bray-Curtis similarity and grouped by habitat (A), diet (B), and family (C).Clusters I and II marked in blue.