Shark tales: A molecular species-level phylogeny of sharks (Selachimorpha, Chondrichthyes)

Assessing the relationships between phylogenetic and functional singularities in sharks (Chondrichthyes)

The relationships between diversity and ecosystem functioning have become a major focus of science. A crucial issue is to estimate functional diversity, as it is intended to impact ecosystem dynamics and stability. However, depending on the ecosystem, it may be challenging or even impossible to directly measure ecological functions and thus functional diversity. Phylogenetic diversity was recently under consideration as a proxy for functional diversity. Phylogenetic diversity is indeed supposed to match functional diversity if functions are conservative traits along evolution. However, in case of adaptive radiation and/or evolutive convergence, a mismatch may appear between species phylogenetic and functional singularities. Using highly threatened taxa, sharks, this study aimed to explore the relationships between phylogenetic and functional diversities and singularities. Different statistical computations were used in order to test both methodological issue (phylogenetic reconstruction) and overall a theoretical questioning: the predictive power of phylogeny for function diversity. Despite these several methodological approaches, a mismatch between phylogeny and function was highlighted. This mismatch revealed that (i) functions are apparently nonconservative in shark species, and (ii) phylogenetic singularity is not a proxy for functional singularity. Functions appeared to be not conservative along the evolution of sharks, raising the conservational challenge to identify and protect both phylogenetic and functional singular species. Facing the current rate of species loss, it is indeed of major importance to target phylogenetically singular species to protect genetic diversity and also functionally singular species in order to maintain particular functions within ecosystem.

Draft sequencing and assembly of the genome of the world's largest fish, the whale shark: Rhincodon typus Smith 1828

BACKGROUND

The whale shark (Rhincodon typus) has by far the largest body size of any elasmobranch (shark or ray) species. Therefore, it is also the largest extant species of the paraphyletic assemblage commonly referred to as fishes. As both a phenotypic extreme and a member of the group Chondrichthyes - the sister group to the remaining gnathostomes, which includes all tetrapods and therefore also humans - its genome is of substantial comparative interest. Whale sharks are also listed as an endangered species on the International Union for Conservation of Nature's Red List of threatened species and are of growing popularity as both a target of ecotourism and as a charismatic conservation ambassador for the pelagic ecosystem. A genome map for this species would aid in defining effective conservation units and understanding global population structure.

RESULTS

We characterised the nuclear genome of the whale shark using next generation sequencing (454, Illumina) and de novo assembly and annotation methods, based on material collected from the Georgia Aquarium. The data set consisted of 878,654,233 reads, which yielded a draft assembly of 1,213,200 contigs and 997,976 scaffolds. The estimated genome size was 3.44Gb. As expected, the proteome of the whale shark was most closely related to the only other complete genome of a cartilaginous fish, the holocephalan elephant shark. The whale shark contained a novel Toll-like-receptor (TLR) protein with sequence similarity to both the TLR4 and TLR13 proteins of mammals and TLR21 of teleosts. The data are publicly available on GenBank, FigShare, and from the NCBI Short Read Archive under accession number SRP044374.

CONCLUSIONS

This represents the first shotgun elasmobranch genome and will aid studies of molecular systematics, biogeography, genetic differentiation, and conservation genetics in this and other shark species, as well as providing comparative data for studies of evolutionary biology and immunology across the jawed vertebrate lineages.

Functional analysis of the musculo-skeletal system of the gill apparatus in Heptranchias perlo (Chondrichthyes: Hexanchidae)

Musculo-skeletal morphology is an indispensable source for understanding functional adaptations. Analysis of morphology of the branchial apparatus of Hexanchiform sharks can provide insight into aspects of their respiration that are difficult to observe directly. In this study, I compare the structure of the musculo-skeletal system of the gill apparatus of Heptranchias perlo and Squalus acanthias in respect to their adaptation for one of two respiratory mechanisms known in sharks, namely, the active two-pump (oropharyngeal and parabranchial) ventilation and the ram-jet ventilation. In both species, the oropharyngeal pump possesses two sets of muscles, one for compression and the other for expansion. The parabranchial pump only has constrictors. Expansion of this pump occurs only due to passive elastic recoil of the extrabranchial cartilages. In Squalus acanthias the parabranchial chambers are large and equipped by powerful superficial constrictors. These muscles and the outer walls of the parabranchial chambers are much reduced in Heptranchias perlo, and thus it likely cannot use this pump. However, this reduction allows for vertical elongation of outer gill slits which, along with greater number of gill pouches, likely decreases branchial resistance and, at the same time, increases the gill surface area, and can be regarded as an adaptation for ram ventilation at lower speeds.

Fin modules: an evolutionary perspective on appendage disparity in basal vertebrates

BACKGROUND

Fishes are extremely speciose and also highly disparate in their fin configurations, more specifically in the number of fins present as well as their structure, shape, and size. How they achieved this remarkable disparity is difficult to explain in the absence of any comprehensive overview of the evolutionary history of fish appendages. Fin modularity could provide an explanation for both the observed disparity in fin configurations and the sequential appearance of new fins. Modularity is considered as an important prerequisite for the evolvability of living systems, enabling individual modules to be optimized without interfering with others. Similarities in developmental patterns between some of the fins already suggest that they form developmental modules during ontogeny. At a macroevolutionary scale, these developmental modules could act as evolutionary units of change and contribute to the disparity in fin configurations. This study addresses fin disparity in a phylogenetic perspective, while focusing on the presence/absence and number of each of the median and paired fins.

RESULTS

Patterns of fin morphological disparity were assessed by mapping fin characters on a new phylogenetic supertree of fish orders. Among agnathans, disparity in fin configurations results from the sequential appearance of novel fins forming various combinations. Both median and paired fins would have appeared first as elongated ribbon-like structures, which were the precursors for more constricted appendages. Among chondrichthyans, disparity in fin configurations relates mostly to median fin losses. Among actinopterygians, fin disparity involves fin losses, the addition of novel fins (e.g., the adipose fin), and coordinated duplications of the dorsal and anal fins. Furthermore, some pairs of fins, notably the dorsal/anal and pectoral/pelvic fins, show non-independence in their character distribution, supporting expectations based on developmental and morphological evidence that these fin pairs form evolutionary modules.

CONCLUSIONS

Our results suggest that the pectoral/pelvic fins and the dorsal/anal fins form two distinct evolutionary modules, and that the latter is nested within a more inclusive median fins module. Because the modularity hypotheses that we are testing are also supported by developmental and variational data, this constitutes a striking example linking developmental, variational, and evolutionary modules.

Squamation and ecology of thelodonts

Thelodonts are an enigmatic group of Paleozoic jawless vertebrates that have been well studied from taxonomical, biostratigraphic and paleogeographic points of view, although our knowledge of their ecology and mode of life is still scant. Their bodies were covered by micrometric scales whose morphology, histology and the developmental process are extremely similar to those of extant sharks. Based on these similarities and on the well-recognized relationship between squamation and ecology in sharks, here we explore the ecological diversity and lifestyles of thelodonts. For this we use classic morphometrics and discriminant analysis to characterize the squamation patterns of a significant number of extant shark species whose ecology is well known. Multivariate analyses have defined a characteristic squamation pattern for each ecological group, thus establishing a comparative framework for inferring lifestyles in thelodonts. We then use this information to study the squamation of the currently described 147 species of thelodonts, known from both articulated and disarticulated remains. Discriminant analysis has allowed recognizing squamation patterns comparable to those of sharks and links them to specific ecological groups. Our results suggest a remarkable ecological diversity in thelodonts. A large number of them were probably demersal species inhabiting hard substrates, within caves and crevices in rocky environments or reefs, taking advantage of the flexibility provided by their micromeric squamations. Contrary to classical interpretations, only few thelodonts were placed among demersal species inhabiting sandy and muddy substrates. Schooling species with defensive scales against ectoparasites could be also abundant suggesting that social interactions and pressure of ectoparasites were present in vertebrates as early the Silurian. The presence of species showing scales suggestive of low to moderate speed and a lifestyle presumably associated with open water environments indicates adaptation of thelodonts to deep water habitats. Scale morphology suggests that some other thelodonts were strong-swimming pelagic species, most of them radiating during the Early Devonian in association with the Nekton Revolution.

Cutting blade dentitions in squaliform sharks form by modification of inherited alternate tooth ordering patterns

The squaliform sharks represent one of the most speciose shark clades. Many adult squaliforms have blade-like teeth, either on both jaws or restricted to the lower jaw, forming a continuous, serrated blade along the jaw margin. These teeth are replaced as a single unit and successor teeth lack the alternate arrangement present in other elasmobranchs. Micro-CT scans of embryos of squaliforms and a related outgroup (Pristiophoridae) revealed that the squaliform dentition pattern represents a highly modified version of tooth replacement seen in other clades. Teeth of Squalus embryos are arranged in an alternate pattern, with successive tooth rows containing additional teeth added proximally. Asynchronous timing of tooth production along the jaw and tooth loss prior to birth cause teeth to align in oblique sets containing teeth from subsequent rows; these become parallel to the jaw margin during ontogeny, so that adult Squalus has functional tooth rows comprising obliquely stacked teeth of consecutive developmental rows. In more strongly heterodont squaliforms, initial embryonic lower teeth develop into the oblique functional sets seen in adult Squalus, with no requirement to form, and subsequently lose, teeth arranged in an initial alternate pattern.

Substitutions in the Glycogenin-1 Gene Are Associated with the Evolution of Endothermy in Sharks and Tunas

Despite 400–450 million years of independent evolution, a strong phenotypic convergence has occurred between two groups of fish: tunas and lamnid sharks. This convergence is characterized by centralization of red muscle, a distinctive swimming style (stiffened body powered through tail movements) and elevated body temperature (endothermy). Furthermore, both groups demonstrate elevated white muscle metabolic capacities. All these traits are unusual in fish and more likely evolved to support their fast-swimming, pelagic, predatory behavior. Here, we tested the hypothesis that their convergent evolution was driven by selection on a set of metabolic genes. We sequenced white muscle transcriptomes of six tuna, one mackerel, and three shark species, and supplemented this data set with previously published RNA-seq data. Using 26 species in total (including 7,032 tuna genes plus 1,719 shark genes), we constructed phylogenetic trees and carried out maximum-likelihood analyses of gene selection. We inferred several genes relating to metabolism to be under selection. We also found that the same one gene, glycogenin-1, evolved under positive selection independently in tunas and lamnid sharks, providing evidence of convergent selective pressures at gene level possibly underlying shared physiology.

Complete mitochondrial genome of the porbeagle shark, Lamna nasus (Chondrichthyes, Lamnidae)

We report for the first time, the complete mitochondrial genome sequence of the porbeagle shark, Lamna nasus, from a specimen collected from offshore waters of New England, USA in the western North Atlantic Ocean. The genome structure of this species is similar to the other reported shark mitogenomes. The genome sequence has a total length of 16,697 bases; similar in size to the mtDNA genomes reported for other lamnid species. A Bayesian phylogenetic tree was reconstructed for the Lamnidae family using mitogenome sequences available in the Genbank.

Do saline taxa evolve faster? Comparing relative rates of molecular evolution between freshwater and marine eukaryotes

The major branches of life diversified in the marine realm, and numerous taxa have since transitioned between marine and freshwaters. Previous studies have demonstrated higher rates of molecular evolution in crustaceans inhabiting continental saline habitats as compared with freshwaters, but it is unclear whether this trend is pervasive or whether it applies to the marine environment. We employ the phylogenetic comparative method to investigate relative molecular evolutionary rates between 148 pairs of marine or continental saline versus freshwater lineages representing disparate eukaryote groups, including bony fish, elasmobranchs, cetaceans, crustaceans, mollusks, annelids, algae, and other eukaryotes, using available protein-coding and noncoding genes. Overall, we observed no consistent pattern in nucleotide substitution rates linked to habitat across all genes and taxa. However, we observed some trends of higher evolutionary rates within protein-coding genes in freshwater taxa-the comparisons mainly involving bony fish-compared with their marine relatives. The results suggest no systematic differences in substitution rate between marine and freshwater organisms.

Repeated and Widespread Evolution of Bioluminescence in Marine Fishes

Bioluminescence is primarily a marine phenomenon with 80% of metazoan bioluminescent genera occurring in the world's oceans. Here we show that bioluminescence has evolved repeatedly and is phylogenetically widespread across ray-finned fishes. We recover 27 independent evolutionary events of bioluminescence, all among marine fish lineages. This finding indicates that bioluminescence has evolved many more times than previously hypothesized across fishes and the tree of life. Our exploration of the macroevolutionary patterns of bioluminescent lineages indicates that the present day diversity of some inshore and deep-sea bioluminescent fish lineages that use bioluminescence for communication, feeding, and reproduction exhibit exceptional species richness given clade age. We show that exceptional species richness occurs particularly in deep-sea fishes with intrinsic bioluminescent systems and both shallow water and deep-sea lineages with luminescent systems used for communication.

Biofluorescence in Catsharks (Scyliorhinidae): Fundamental Description and Relevance for Elasmobranch Visual Ecology

Biofluorescence has recently been found to be widespread in marine fishes, including sharks. Catsharks, such as the Swell Shark (Cephaloscyllium ventriosum) from the eastern Pacific and the Chain Catshark (Scyliorhinus retifer) from the western Atlantic, are known to exhibit bright green fluorescence. We examined the spectral sensitivity and visual characteristics of these reclusive sharks, while also considering the fluorescent properties of their skin. Spectral absorbance of the photoreceptor cells in these sharks revealed the presence of a single visual pigment in each species. Cephaloscyllium ventriosum exhibited a maximum absorbance of 484 ± 3 nm and an absorbance range at half maximum (λ1/2max) of 440-540 nm, whereas for S. retifer maximum absorbance was 488 ± 3 nm with the same absorbance range. Using the photoreceptor properties derived here, a "shark eye" camera was designed and developed that yielded contrast information on areas where fluorescence is anatomically distributed on the shark, as seen from other sharks' eyes of these two species. Phylogenetic investigations indicate that biofluorescence has evolved at least three times in cartilaginous fishes. The repeated evolution of biofluorescence in elasmobranchs, coupled with a visual adaptation to detect it; and evidence that biofluorescence creates greater luminosity contrast with the surrounding background, highlights the potential importance of biofluorescence in elasmobranch behavior and biology.

Complete mitochondrial genome and the phylogenetic position of the graceful catshark Proscyllium habereri (Carcharhiniformes: Proscylliidae)

In this study, the complete mitogenome of Proscyllium habereri (Carcharhiniformes: Proscylliidae) is first determined. It is 16,708 bp in length, containing 37 genes with typical order to that of most other vertebrates. Its overall base composition of the H-strand is A 30.9%; C 23.7%; G 14.2%; T 31.2%. Two start codons (ATG and GTG) and two stop codons (TAG and TAA/T) are found in the protein-coding genes. The 22 tRNA genes range from 67 bp (tRNA-Cys, tRNA-Ser2) to 75 bp (tRNA-Leu1). The phylogenetic result showed that P. habereri was clustered to Pseudotriakis microdon.

Complete mitochondrial DNA genome of bonnethead shark, Sphyrna tiburo, and phylogenetic relationships among main superorders of modern elasmobranchs

Elasmobranchs are one of the most diverse groups in the marine realm represented by 18 orders, 55 families and about 1200 species reported, but also one of the most vulnerable to exploitation and to climate change. Phylogenetic relationships among main orders have been controversial since the emergence of the Hypnosqualean hypothesis by Shirai (1992) that considered batoids as a sister group of sharks. The use of the complete mitochondrial DNA (mtDNA) may shed light to further validate this hypothesis by increasing the number of informative characters. We report the mtDNA genome of the bonnethead shark Sphyrna tiburo, and compare it with mitogenomes of other 48 species to assess phylogenetic relationships. The mtDNA genome of S. tiburo, is quite similar in size to that of congeneric species but also similar to the reported mtDNA genome of other Carcharhinidae species. Like most vertebrate mitochondrial genomes, it contained 13 protein coding genes, two rRNA genes and 22 tRNA genes and the control region of 1086 bp (D-loop). The Bayesian analysis of the 49 mitogenomes supported the view that sharks and batoids are separate groups.

An exceptionally preserved Eocene shark and the rise of modern predator-prey interactions in the coral reef food web

BACKGROUND

Following extreme climatic warming events, Eocene Lagerstätten document aquatic and terrestrial vertebrate faunas surprisingly similar to modern counterparts. This transition in marine systems is best documented in the earliest teleost-dominated coral reef assemblage of Pesciara di Bolca, northern Italy, from near the end of the Eocene Climatic Optimum. Its rich fauna shows similarities with that of the modern Great Barrier Reef in niche exploitation by and morphological disparity among teleost primary consumers. However, such paleoecological understanding has not transcended trophic levels above primary consumers, particularly in carcharhiniform sharks.

RESULTS

We report an exceptionally preserved fossil school shark (Galeorhinus cuvieri) from Pesciara di Bolca. In addition to the spectacular preservation of soft tissues, including brain, muscles, and claspers, this male juvenile shark has stomach contents clearly identifiable as a sphyraenid acanthomorph (barracuda). This association provides evidence that a predator-prey relationship between Galeorhinus and Sphyraena in the modern coral reefs has roots in the Eocene. A growth curve of the living species of Galeorhinus fitted to G. cuvieri suggests that all specimens of G. cuvieri from the lagoonal deposits of Bolca represent sexually and somatically immature juveniles.

CONCLUSION

The modern trophic association between higher-degree consumers (Galeorhinus and Sphyraena) has a counterpart in the Eocene Bolca, just as Bolca and the Great Barrier Reef show parallels among teleost primary consumers. Given the age of Bolca, trophic networks among consumers observed in modern coral reefs arose by the exit from the Climatic Optimum. The biased representation of juveniles suggests that the Bolca Lagerstätte served as a nursery habitat for G. cuvieri. Ultraviolet photography may be useful in probing for exceptional soft tissue preservation before common acid preparation methods.

The phylogenomic position of the grey nurse shark Carcharias taurus Rafinesque, 1810 (Lamniformes, Odontaspididae) inferred from the mitochondrial genome

The complete mitochondrial genome of the grey nurse shark Carcharias taurus is described from 25 963 828 sequences obtained using Illumina NGS technology. Total length of the mitogenome is 16 715 bp, consisting of 2 rRNAs, 13 protein-coding regions, 22 tRNA and 2 non-coding regions thus updating the previously published mitogenome for this species. The phylogenomic reconstruction inferred from the mitogenome of 15 species of Lamniform and Carcharhiniform sharks supports the inclusion of C. taurus in a clade with the Lamnidae and Cetorhinidae. This complete mitogenome contributes to ongoing investigation into the monophyly of the Family Odontaspididae.

Predicting loss of evolutionary history: Where are we?

The Earth's evolutionary history is threatened by species loss in the current sixth mass extinction event in Earth's history. Such extinction events not only eliminate species but also their unique evolutionary histories. Here we review the expected loss of Earth's evolutionary history quantified by phylogenetic diversity (PD) and evolutionary distinctiveness (ED) at risk. Due to the general paucity of data, global evolutionary history losses have been predicted for only a few groups, such as mammals, birds, amphibians, plants, corals and fishes. Among these groups, there is now empirical support that extinction threats are clustered on the phylogeny; however this is not always a sufficient condition to cause higher loss of phylogenetic diversity in comparison to a scenario of random extinctions. Extinctions of the most evolutionarily distinct species and the shape of phylogenetic trees are additional factors that can elevate losses of evolutionary history. Consequently, impacts of species extinctions differ among groups and regions, and even if global losses are low within large groups, losses can be high among subgroups or within some regions. Further, we show that PD and ED are poorly protected by current conservation practices. While evolutionary history can be indirectly protected by current conservation schemes, optimizing its preservation requires integrating phylogenetic indices with those that capture rarity and extinction risk. Measures based on PD and ED could bring solutions to conservation issues, however they are still rarely used in practice, probably because the reasons to protect evolutionary history are not clear for practitioners or due to a lack of data. However, important advances have been made in the availability of phylogenetic trees and methods for their construction, as well as assessments of extinction risk. Some challenges remain, and looking forward, research should prioritize the assessment of expected PD and ED loss for more taxonomic groups and test the assumption that preserving ED and PD also protects rare species and ecosystem services. Such research will be useful to inform and guide the conservation of Earth's biodiversity and the services it provides.

Mitochondrial genome and phylogenetic position of the sliteye shark Loxodon macrorhinus

The sliteye shark Loxodon macrorhinus is the only member of the genus Loxodon in the family Carcharhinidae. In this study, we first present the complete mitochondrial genome of L. macrorhinus and determine its phylogenetic position within Carcharhinidae based on relative mitogenomes. The mitochondrial genome was 16 702 bp in length with the typical gene order in vertebrates. The overall base composition of the H-strand was 31.7% A, 25.8% C, 13.1% G, and 29.4% T. Two start codons (ATG and GTG) and three stop codons (TAG, AGG, and TAA/T) were found in the protein-coding genes. The tRNA genes ranged from 67 bp to 75 bp. Loxodon macrorhinus was placed as sister to the genus Scoliodon in the Bayesian tree.

Global pattern of phylogenetic species composition of shark and its conservation priority

The diversity of marine communities is in striking contrast with the diversity of terrestrial communities. In all oceans, species richness is low in tropical areas and high at latitudes between 20 and 40°. While species richness is a primary metric used in conservation and management strategies, it is important to take into account the complex phylogenetic patterns of species compositions within communities. We measured the phylogenetic skew and diversity of shark communities throughout the world. We found that shark communities in tropical seas were highly phylogenetically skewed, whereas temperate sea communities had phylogenetically diversified species compositions. Interestingly, although geographically distant from one another, tropical sea communities were all highly skewed toward requiem sharks (Carcharhinidae), hammerhead sharks (Sphyrnidae), and whale sharks (Rhincodon typus). Worldwide, the greatest phylogenetic evenness in terms of clades was found in the North Sea and coastal regions of countries in temperate zones, such as the United Kingdom, Ireland, southern Australia, and Chile. This study is the first to examine patterns of phylogenetic diversity of shark communities on a global scale. Our findings suggest that when establishing conservation activities, it is important to take full account of phylogenetic patterns of species composition and not solely use species richness as a target. Protecting areas of high phylogenetic diversity in sharks, which were identified in this study, could form a broader strategy for protecting other threatened marine species.

Molecular phylogeny of Squaliformes and first occurrence of bioluminescence in sharks

Background

Squaliform sharks represent approximately 27 % of extant shark diversity, comprising more than 130 species with a predominantly deep-dwelling lifestyle. Many Squaliform species are highly specialized, including some that are bioluminescent, a character that is reported exclusively from Squaliform sharks within Chondrichthyes. The interfamiliar relationships within the order are still not satisfactorily resolved. Herein we estimate the phylogenetic interrelationships of a generic level sampling of “squaloid” sharks and closely related taxa using aligned sequences derived from a targeted gene capture approach. The resulting phylogenetic estimate is further used to evaluate the age of first occurrence of bioluminescence in Squaliformes.

Results

Our dataset comprised 172 putative ortholog exon sequences. Phylogenetic estimates result in a fully resolved tree supporting a monophyletic lineage of Squaliformes excluding Echinorhinus. Non-luminous Squalidae are inferred to be the sister to a clade comprising all remaining Squaliform families. Our results suggest that the origin of photophores is coincident with an elevated diversification rate and the splitting of families Dalatiidae, Etmopteridae, Oxynotidae and Somniosidae at the transition of the Lower to the Upper Cretaceous. The presence of luminous organs was confirmed for the Sleeper shark genus Zameus. These results indicate that bioluminescence in sharks is not restricted solely to the families Etmopteridae and Dalatiidae as previously believed.

Conclusions

The sister-clade to non-luminous Squalidae comprises five families. The presence of photophores is reported for extant members of three out of these five families based on results of this study, i.e. Lantern sharks (Etmopteridae), Kitefin sharks (Dalatiidae) and Sleeper sharks (Somniosidae). Our results suggest that the origin of luminous organs arose during the rapid diversification event that gave rise to the extant Squaliform families. These inferences are consistent with the idea of diversification of Squaliform sharks being associated with the emergence of new deep-sea habitats in the Lower Cretaceous, which may have been facilitated by the evolution of bioluminescence.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0446-6) contains supplementary material, which is available to authorized users.

Complete mitochondrial genome and the phylogenetic position of the Blotchy swell shark Cephaloscyllium umbratile

In this study, the complete mitochondrial genome of the Blotchy swell shark Cephaloscyllium umbratile was determined. It was a circle molecular (16 698 bp), contained 37 genes with typical order to that of most other vertebrates. The nucleotide composition was 31.0% A, 24.0% C, 14.0% G, and 31.3% T. There were 26 bp short intergenic spaces located in 11 gene junctions and 28 bp overlaps located in 7 gene junctions in the whole mitogenome. Two start codons (GTG and ATG) and two stop codons (TAG and TAA/T) were used in the protein-coding genes. The phylogenetic result showed that C. umbratile was clustered with Scyliorhinus canicula and formed the Scyliorhinidae clade, which was the most basal clade within Carcharhiniformes, and Carcharhinidae is not monophyletic.

Digestive enzyme activities are higher in the shortfin mako shark, Isurus oxyrinchus, than in ectothermic sharks as a result of visceral endothermy

Lamnid sharks are regionally endothermic fishes that maintain visceral temperatures elevated above the ambient water temperature. Visceral endothermy is thought to increase rates of digestion and food processing and allow thermal niche expansion. We tested the hypothesis that, at in vivo temperatures, the endothermic shortfin mako shark, Isurus oxyrinchus, has higher specific activities of three digestive enzymes-gastric pepsin and pancreatic trypsin and lipase-than the thresher shark, Alopias vulpinus, and the blue shark, Prionace glauca, neither of which can maintain elevated visceral temperatures. Homogenized stomach or pancreas tissue obtained from sharks collected by pelagic longline was incubated at both 15 and 25 °C, at saturating substrate concentrations, to quantify tissue enzymatic activity. The mako had significantly higher enzyme activities at 25 °C than did the thresher and blue sharks at 15 °C. This difference was not a simple temperature effect, because at 25 °C the mako had higher trypsin activity than the blue shark and higher activities for all enzymes than the thresher shark. We also hypothesized that the thermal coefficient, or Q 10 value, would be higher for the mako shark than for the thresher and blue sharks because of its more stable visceral temperature. However, the mako and thresher sharks had similar Q 10 values for all enzymes, perhaps because of their closer phylogenetic relationship. The higher in vivo digestive enzyme activities in the mako shark should result in higher rates of food processing and may represent a selective advantage of regional visceral endothermy.

An unusual mitochondrial genome structure of the tonguefish, Cynoglossus trigrammus: Control region translocation and a long additional non-coding region inversion

Flatfishes (Pleuronectiformes) exhibit different types of large-scale gene rearrangements. In the present study, the mitochondrial (mt) genome (18,369bp) of a tonguefish, Cynoglossus trigrammus, was determined using de novo mitochondrion genome sequencing. Compared with other flatfishes, the mt genome of C. trigrammus revealed distinct mitogenome architectures that primarily included two striking findings: 1) insertion of an additional long non-coding region (1647bp) making it the second largest genome length among Pleuronectiformes and 2) the translocation of the control region. The reconstructed phylogenetic tree based on 13 mt protein-coding gene sequences recovered the monophyletic suborder Pleuronectoidei and the family Cynoglossidae. These data provide useful information for a better understanding of the mitogenomic diversities and evolution in fish as well as novel genetic markers for studying population genetics and species identification.

'Fish' (Actinopterygii and Elasmobranchii) diversification patterns through deep time

Actinopterygii (ray-finned fishes) and Elasmobranchii (sharks, skates and rays) represent more than half of today's vertebrate taxic diversity (approximately 33000 species) and form the largest component of vertebrate diversity in extant aquatic ecosystems. Yet, patterns of 'fish' evolutionary history remain insufficiently understood and previous studies generally treated each group independently mainly because of their contrasting fossil record composition and corresponding sampling strategies. Because direct reading of palaeodiversity curves is affected by several biases affecting the fossil record, analytical approaches are needed to correct for these biases. In this review, we propose a comprehensive analysis based on comparison of large data sets related to competing phylogenies (including all Recent and fossil taxa) and the fossil record for both groups during the Mesozoic-Cainozoic interval. This approach provides information on the 'fish' fossil record quality and on the corrected 'fish' deep-time phylogenetic palaeodiversity signals, with special emphasis on diversification events. Because taxonomic information is preserved after analytical treatment, identified palaeodiversity events are considered both quantitatively and qualitatively and put within corresponding palaeoenvironmental and biological settings. Results indicate a better fossil record quality for elasmobranchs due to their microfossil-like fossil distribution and their very low diversity in freshwater systems, whereas freshwater actinopterygians are diverse in this realm with lower preservation potential. Several important diversification events are identified at familial and generic levels for elasmobranchs, and marine and freshwater actinopterygians, namely in the Early-Middle Jurassic (elasmobranchs), Late Jurassic (actinopterygians), Early Cretaceous (elasmobranchs, freshwater actinopterygians), Cenomanian (all groups) and the Paleocene-Eocene interval (all groups), the latter two representing the two most exceptional radiations among vertebrates. For each of these events along with the Cretaceous-Paleogene extinction, we provide an in-depth review of the taxa involved and factors that may have influenced the diversity patterns observed. Among these, palaeotemperatures, sea-levels, ocean circulation and productivity as well as continent fragmentation and environment heterogeneity (reef environments) are parameters that largely impacted on 'fish' evolutionary history, along with other biotic constraints.

Comparative analyses of animal-tracking data reveal ecological significance of endothermy in fishes

Despite long evolutionary separations, several sharks and tunas share the ability to maintain slow-twitch, aerobic red muscle (RM) warmer than ambient water. Proximate causes of RM endothermy are well understood, but ultimate causes are unclear. Two advantages often proposed are thermal niche expansion and elevated cruising speeds. The thermal niche hypothesis is generally supported, because fishes with RM endothermy often exhibit greater tolerance to broad temperature ranges. In contrast, whether fishes with RM endothermy cruise faster, and achieve any ecological benefits from doing so, remains unclear. Here, we compiled data recorded by modern animal-tracking tools for a variety of free-swimming marine vertebrates. Using phylogenetically informed allometry, we show that both cruising speeds and maximum annual migration ranges of fishes with RM endothermy are 2-3 times greater than fishes without it, and comparable to nonfish endotherms (i.e., penguins and marine mammals). The estimated cost of transport of fishes with RM endothermy is twice that of fishes without it. We suggest that the high energetic cost of RM endothermy in fishes is offset by the benefit of elevated cruising speeds, which not only increase prey encounter rates, but also enable larger-scale annual migrations and potentially greater access to seasonally available resources.

Development and evolution of dentition pattern and tooth order in the skates and rays (batoidea; chondrichthyes)

Shark and ray (elasmobranch) dentitions are well known for their multiple generations of teeth, with isolated teeth being common in the fossil record. However, how the diverse dentitions characteristic of elasmobranchs form is still poorly understood. Data on the development and maintenance of the dental patterning in this major vertebrate group will allow comparisons to other morphologically diverse taxa, including the bony fishes, in order to identify shared pattern characters for the vertebrate dentition as a whole. Data is especially lacking from the Batoidea (skates and rays), hence our objective is to compile data on embryonic and adult batoid tooth development contributing to ordering of the dentition, from cleared and stained specimens and micro-CT scans, with 3D rendered models. We selected species (adult and embryonic) spanning phylogenetically significant batoid clades, such that our observations may raise questions about relationships within the batoids, particularly with respect to current molecular-based analyses. We include developmental data from embryos of recent model organisms Leucoraja erinacea and Raja clavata to evaluate the earliest establishment of the dentition. Characters of the batoid dentition investigated include alternate addition of teeth as offset successional tooth rows (versus single separate files), presence of a symphyseal initiator region (symphyseal tooth present, or absent, but with two parasymphyseal teeth) and a restriction to tooth addition along each jaw reducing the number of tooth families, relative to addition of successor teeth within each family. Our ultimate aim is to understand the shared characters of the batoids, and whether or not these dental characters are shared more broadly within elasmobranchs, by comparing these to dentitions in shark outgroups. These developmental morphological analyses will provide a solid basis to better understand dental evolution in these important vertebrate groups as well as the general plesiomorphic vertebrate dental condition.

Feeding performance of king Mackerel, Scomberomorus cavalla

Feeding performance is an organism's ability to capture and handle prey. Although bite force is a commonly used metric of feeding performance, other factors such as bite pressure and strike speed are also likely to affect prey capture. Therefore, this study investigated static bite force, dynamic speeds, and predator and prey forces resulting from ram strikes, as well as bite pressure of the king mackerel, Scomberomorus cavalla, in order to examine their relative contributions to overall feeding performance. Theoretical posterior bite force ranged from 14.0-318.7 N. Ram speed, recorded with a rod and reel incorporated with a line counter and video camera, ranged from 3.3-15.8B L/s. Impact forces on the prey ranged from 0.1-1.9 N. Bite pressure, estimated using theoretical bite forces at three gape angles and tooth cross-sectional areas, ranged from 1.7-56.9 MPa. Mass-specific bite force for king mackerel is relatively low in comparison with other bony fishes and sharks, with relatively little impact force applied to the prey during the strike. This suggests that king mackerel rely on high velocity chases and high bite pressure generated via sharp, laterally compressed teeth to maximize feeding performance.

Constructional morphology within the head of hammerhead sharks (sphyrnidae)

The study of functional trade-offs is important if a structure, such as the cranium, serves multiple biological roles, and is, therefore, shaped by multiple selective pressures. The sphyrnid cephalofoil presents an excellent model for investigating potential trade-offs among sensory, neural, and feeding structures. In this study, hammerhead shark species were chosen to represent differences in head form through phylogeny. A combination of surface-based geometric morphometrics, computed tomography (CT) volumetric analysis, and phylogenetic analyses were utilized to investigate potential trade-offs within the head. Hammerhead sharks display a diversity of cranial morphologies where the position of the eyes and nares vary among species, with only minor changes in shape, position, and volume of the feeding apparatus through phylogeny. The basal winghead shark, Eusphyra blochii, has small anteriorly positioned eyes. Through phylogeny, the relative size and position of the eyes change, such that derived species have larger, more medially positioned eyes. The lateral position of the external nares is highly variable, showing no phylogenetic trend. Mouth size and position are conserved, remaining relatively unchanged. Volumetric CT analyses reveal no trade-offs between the feeding apparatus and the remaining cranial structures. The few trade-offs were isolated to the nasal capsule volume's inverse correlation with braincase, chondrocranial, and total cephalofoil volume. Eye volume also decreased as cephalofoil width increased. These data indicate that despite considerable changes in head shape, much of the head is morphologically conserved through sphyrnid phylogeny, particularly the jaw cartilages and their associated feeding muscles, with shape change and morphological trade-offs being primarily confined to the lateral wings of the cephalofoil and their associated sensory structures.

Characterization of the Macropodus opercularis complete mitochondrial genome and family Channidae taxonomy using Illumina-based de novo transcriptome sequencing

In this study, the complete mitochondrial genome of Macropodus opercularis was sequenced using Illumina-based de novo transcriptome technology and annotated using bioinformatic tools. The circular mitochondrial genome was 16,496bp in length and contained two ribosomal RNAs, 13 protein-coding genes, 22 transfer RNA genes, and the control region. The gene composition and order were similar to suborder Anabantoidei. Phylogenetic analyses using concatenated amino acid and nucleotide sequences of the 13 protein-coding genes with two different methods (Neighbor-joining and Bayesian analysis) both highly supported the close relationship of M. opercularis to M. ocellatus, consistent with previous classifications based on morphological and molecular studies. Furthermore, family Channidae and Parachanna insignis were clustered in the same clade. Our results supported the inclusion of family Channidae in suborder Channoidei. The complete mitochondrial genome of M. opercularis will provide genetic markers for better understanding species identification, population genetics and phylogeographics of freshwater fishes.

Molecular phylogeny of the bamboo sharks (Chiloscyllium spp.)

Chiloscyllium, commonly called bamboo shark, can be found inhabiting the waters of the Indo-West Pacific around East Asian countries such as Malaysia, Myanmar, Thailand, Singapore, and Indonesia. The International Union for Conservation of Nature (IUCN) Red List has categorized them as nearly threatened sharks out of their declining population status due to overexploitation. A molecular study was carried out to portray the systematic relationships within Chiloscyllium species using 12S rRNA and cytochrome b gene sequences. Maximum parsimony and Bayesian were used to reconstruct their phylogeny trees. A total of 381 bp sequences' lengths were successfully aligned in the 12S rRNA region, with 41 bp sites being parsimony-informative. In the cytochrome b region, a total of 1120 bp sites were aligned, with 352 parsimony-informative characters. All analyses yield phylogeny trees on which C. indicum has close relationships with C. plagiosum. C. punctatum is sister taxon to both C. indicum and C. plagiosum while C. griseum and C. hasseltii formed their own clade as sister taxa. These Chiloscyllium classifications can be supported by some morphological characters (lateral dermal ridges on the body, coloring patterns, and appearance of hypobranchials and basibranchial plate) that can clearly be used to differentiate each species.

Dinosaur physiology. Evidence for mesothermy in dinosaurs

Were dinosaurs ectotherms or fast-metabolizing endotherms whose activities were unconstrained by temperature? To date, some of the strongest evidence for endothermy comes from the rapid growth rates derived from the analysis of fossil bones. However, these studies are constrained by a lack of comparative data and an appropriate energetic framework. Here we compile data on ontogenetic growth for extant and fossil vertebrates, including all major dinosaur clades. Using a metabolic scaling approach, we find that growth and metabolic rates follow theoretical predictions across clades, although some groups deviate. Moreover, when the effects of size and temperature are considered, dinosaur metabolic rates were intermediate to those of endotherms and ectotherms and closest to those of extant mesotherms. Our results suggest that the modern dichotomy of endothermic versus ectothermic is overly simplistic.

DNA barcoding of shark meats identify species composition and CITES-listed species from the markets in Taiwan

BACKGROUND

An increasing awareness of the vulnerability of sharks to exploitation by shark finning has contributed to a growing concern about an unsustainable shark fishery. Taiwan's fleet has the 4th largest shark catch in the world, accounting for almost 6% of the global figures. Revealing the diversity of sharks consumed by Taiwanese is important in designing conservation plans. However, fins make up less than 5% of the total body weight of a shark, and their bodies are sold as filets in the market, making it difficult or impossible to identify species using morphological traits.

METHODS

In the present study, we adopted a DNA barcoding technique using a 391-bp fragment of the mitochondrial cytochrome oxidase I (COI) gene to examine the diversity of shark filets and fins collected from markets and restaurants island-wide in Taiwan.

RESULTS

Amongst the 548 tissue samples collected and sequenced, 20 major clusters were apparent by phylogenetic analyses, each of them containing individuals belonging to the same species (most with more than 95% bootstrap values), corresponding to 20 species of sharks. Additionally, Alopias pelagicus, Carcharhinus falciformis, Isurus oxyrinchus, and Prionace glauca consisted of 80% of the samples we collected, indicating that these species might be heavily consumed in Taiwan. Approximately 5% of the tissue samples used in this study were identified as species listed in CITES Appendix II, including two species of Sphyrna, C. longimanus and Carcharodon carcharias.

CONCLUSION

DNA barcoding provides an alternative method for understanding shark species composition when species-specific data is unavailable. Considering the global population decline, stock assessments of Appendix II species and highly consumed species are needed to accomplish the ultimate goal of shark conservation.

Evolutionary relations of Hexanchiformes deep-sea sharks elucidated by whole mitochondrial genome sequences

Hexanchiformes is regarded as a monophyletic taxon, but the morphological and genetic relationships between the five extant species within the order are still uncertain. In this study, we determined the whole mitochondrial DNA (mtDNA) sequences of seven sharks including representatives of the five Hexanchiformes, one squaliform, and one carcharhiniform and inferred the phylogenetic relationships among those species and 12 other Chondrichthyes (cartilaginous fishes) species for which the complete mitogenome is available. The monophyly of Hexanchiformes and its close relation with all other Squaliformes sharks were strongly supported by likelihood and Bayesian phylogenetic analysis of 13,749 aligned nucleotides of 13 protein coding genes and two rRNA genes that were derived from the whole mDNA sequences of the 19 species. The phylogeny suggested that Hexanchiformes is in the superorder Squalomorphi, Chlamydoselachus anguineus (frilled shark) is the sister species to all other Hexanchiformes, and the relations within Hexanchiformes are well resolved as Chlamydoselachus, (Notorynchus, (Heptranchias, (Hexanchus griseus, H. nakamurai))). Based on our phylogeny, we discussed evolutionary scenarios of the jaw suspension mechanism and gill slit numbers that are significant features in the sharks.

Body wall development in lamprey and a new perspective on the origin of vertebrate paired fins

Classical hypotheses regarding the evolutionary origin of paired appendages propose transformation of precursor structures (gill arches and lateral fin folds) into paired fins. During development, gnathostome paired appendages form as outgrowths of body wall somatopleure, a tissue composed of somatic lateral plate mesoderm (LPM) and overlying ectoderm. In amniotes, LPM contributes connective tissue to abaxial musculature and forms ventrolateral dermis of the interlimb body wall. The phylogenetic distribution of this character is uncertain because lineage analyses of LPM have not been generated in anamniotes. We focus on the evolutionary history of the somatopleure to gain insight into the tissue context in which paired fins first appeared. Lampreys diverged from other vertebrates before the acquisition of paired fins and provide a model for investigating the preappendicular condition. We present vital dye fate maps that suggest the somatopleure is eliminated in lamprey as the LPM is separated from the ectoderm and sequestered to the coelomic linings during myotome extension. We also examine the distribution of postcranial mesoderm in catshark and axolotl. In contrast to lamprey, our findings support an LPM contribution to the trunk body wall of these taxa, which is similar to published data for amniotes. Collectively, these data lead us to hypothesize that a persistent somatopleure in the lateral body wall is a gnathostome synapomorphy, and the redistribution of LPM was a key step in generating the novel developmental module that ultimately produced paired fins. These embryological criteria can refocus arguments on paired fin origins and generate hypotheses testable by comparative studies on the source, sequence, and extent of genetic redeployment.

Revealing less derived nature of cartilaginous fish genomes with their evolutionary time scale inferred with nuclear genes

Cartilaginous fishes, divided into Holocephali (chimaeras) and Elasmoblanchii (sharks, rays and skates), occupy a key phylogenetic position among extant vertebrates in reconstructing their evolutionary processes. Their accurate evolutionary time scale is indispensable for better understanding of the relationship between phenotypic and molecular evolution of cartilaginous fishes. However, our current knowledge on the time scale of cartilaginous fish evolution largely relies on estimates using mitochondrial DNA sequences. In this study, making the best use of the still partial, but large-scale sequencing data of cartilaginous fish species, we estimate the divergence times between the major cartilaginous fish lineages employing nuclear genes. By rigorous orthology assessment based on available genomic and transcriptomic sequence resources for cartilaginous fishes, we selected 20 protein-coding genes in the nuclear genome, spanning 2973 amino acid residues. Our analysis based on the Bayesian inference resulted in the mean divergence time of 421 Ma, the late Silurian, for the Holocephali-Elasmobranchii split, and 306 Ma, the late Carboniferous, for the split between sharks and rays/skates. By applying these results and other documented divergence times, we measured the relative evolutionary rate of the Hox A cluster sequences in the cartilaginous fish lineages, which resulted in a lower substitution rate with a factor of at least 2.4 in comparison to tetrapod lineages. The obtained time scale enables mapping phenotypic and molecular changes in a quantitative framework. It is of great interest to corroborate the less derived nature of cartilaginous fish at the molecular level as a genome-wide phenomenon.

Estimating the rate of irreversibility in protein evolution

Whether or not evolutionary change is inherently irreversible remains a controversial topic. Some examples of evolutionary irreversibility are known; however, this question has not been comprehensively addressed at the molecular level. Here, we use data from 221 human genes with known pathogenic mutations to estimate the rate of irreversibility in protein evolution. For these genes, we reconstruct ancestral amino acid sequences along the mammalian phylogeny and identify ancestral amino acid states that match known pathogenic mutations. Such cases represent inherent evolutionary irreversibility because, at the present moment, reversals to these ancestral amino acid states are impossible for the human lineage. We estimate that approximately 10% of all amino acid substitutions along the mammalian phylogeny are irreversible, such that a return to the ancestral amino acid state would lead to a pathogenic phenotype. For a subset of 51 genes with high rates of irreversibility, as much as 40% of all amino acid evolution was estimated to be irreversible. Because pathogenic phenotypes do not resemble ancestral phenotypes, the molecular nature of the high rate of irreversibility in proteins is best explained by evolution with a high prevalence of compensatory, epistatic interactions between amino acid sites. Under such mode of protein evolution, once an amino acid substitution is fixed, the probability of its reversal declines as the protein sequence accumulates changes that affect the phenotypic manifestation of the ancestral state. The prevalence of epistasis in evolution indicates that the observed high rate of irreversibility in protein evolution is an inherent property of protein structure and function.

Sharing and re-use of phylogenetic trees (and associated data) to facilitate synthesis

BACKGROUND

Recently, various evolution-related journals adopted policies to encourage or require archiving of phylogenetic trees and associated data. Such attention to practices that promote sharing of data reflects rapidly improving information technology, and rapidly expanding potential to use this technology to aggregate and link data from previously published research. Nevertheless, little is known about current practices, or best practices, for publishing trees and associated data so as to promote re-use.

FINDINGS

Here we summarize results of an ongoing analysis of current practices for archiving phylogenetic trees and associated data, current practices of re-use, and current barriers to re-use. We find that the technical infrastructure is available to support rudimentary archiving, but the frequency of archiving is low. Currently, most phylogenetic knowledge is not easily re-used due to a lack of archiving, lack of awareness of best practices, and lack of community-wide standards for formatting data, naming entities, and annotating data. Most attempts at data re-use seem to end in disappointment. Nevertheless, we find many positive examples of data re-use, particularly those that involve customized species trees generated by grafting to, and pruning from, a much larger tree.

CONCLUSIONS

The technologies and practices that facilitate data re-use can catalyze synthetic and integrative research. However, success will require engagement from various stakeholders including individual scientists who produce or consume shareable data, publishers, policy-makers, technology developers and resource-providers. The critical challenges for facilitating re-use of phylogenetic trees and associated data, we suggest, include: a broader commitment to public archiving; more extensive use of globally meaningful identifiers; development of user-friendly technology for annotating, submitting, searching, and retrieving data and their metadata; and development of a minimum reporting standard (MIAPA) indicating which kinds of data and metadata are most important for a re-useable phylogenetic record.

An analytical approach for estimating fossil record and diversification events in sharks, skates and rays

Background

Modern selachians and their supposed sister group (hybodont sharks) have a long and successful evolutionary history. Yet, although selachian remains are considered relatively common in the fossil record in comparison with other marine vertebrates, little is known about the quality of their fossil record. Similarly, only a few works based on specific time intervals have attempted to identify major events that marked the evolutionary history of this group.

Methodology/Principal Findings

Phylogenetic hypotheses concerning modern selachians’ interrelationships are numerous but differ significantly and no consensus has been found. The aim of the present study is to take advantage of the range of recent phylogenetic hypotheses in order to assess the fit of the selachian fossil record to phylogenies, according to two different branching methods. Compilation of these data allowed the inference of an estimated range of diversity through time and evolutionary events that marked this group over the past 300 Ma are identified. Results indicate that with the exception of high taxonomic ranks (orders), the selachian fossil record is by far imperfect, particularly for generic and post-Triassic data. Timing and amplitude of the various identified events that marked the selachian evolutionary history are discussed.

Conclusion/Significance

Some identified diversity events were mentioned in previous works using alternative methods (Early Jurassic, mid-Cretaceous, K/T boundary and late Paleogene diversity drops), thus reinforcing the efficiency of the methodology presented here in inferring evolutionary events. Other events (Permian/Triassic, Early and Late Cretaceous diversifications; Triassic/Jurassic extinction) are newly identified. Relationships between these events and paleoenvironmental characteristics and other groups’ evolutionary history are proposed.