Crocodilian phylogeny inferred from twelve mitochondrial protein-coding genes, with new complete mitochondrial genomic sequences for Crocodylus acutus and Crocodylus novaeguineae

Phylogenetic analysis of a new morphological dataset elucidates the evolutionary history of Crocodylia and resolves the long-standing gharial problem

First appearing in the latest Cretaceous, Crocodylia is a clade of semi-aquatic, predatory reptiles, defined by the last common ancestor of extant alligators, caimans, crocodiles, and gharials. Despite large strides in resolving crocodylian interrelationships over the last three decades, several outstanding problems persist in crocodylian systematics. Most notably, there has been persistent discordance between morphological and molecular datasets surrounding the affinities of the extant gharials, Gavialis gangeticus and Tomistoma schlegelii. Whereas molecular data consistently support a sister taxon relationship, in which they are more closely related to crocodylids than to alligatorids, morphological data indicate that Gavialis is the sister taxon to all other extant crocodylians. Here we present a new morphological dataset for Crocodylia based on a critical reappraisal of published crocodylian character data matrices and extensive firsthand observations of a global sample of crocodylians. This comprises the most taxonomically comprehensive crocodylian dataset to date (144 OTUs scored for 330 characters) and includes a new, illustrated character list with modifications to the construction and scoring of characters, and 46 novel characters. Under a maximum parsimony framework, our analyses robustly recover Gavialis as more closely related to Tomistoma than to other extant crocodylians for the first time based on morphology alone. This result is recovered regardless of the weighting strategy and treatment of quantitative characters. However, analyses using continuous characters and extended implied weighting (with high k-values) produced the most resolved, well-supported, and stratigraphically congruent topologies overall. Resolution of the gharial problem reveals that: (1) several gavialoids lack plesiomorphic features that formerly drew them towards the stem of Crocodylia; and (2) more widespread similarities occur between species traditionally divided into tomistomines and gavialoids, with these interpreted here as homology rather than homoplasy. There remains significant temporal incongruence regarding the inferred divergence timing of the extant gharials, indicating that several putative gavialids ('thoracosaurs') are incorrectly placed and require future re-appraisal. New alligatoroid interrelationships include: (1) support for a North American origin of Caimaninae in the latest Cretaceous; (2) the recovery of the early Paleogene South American taxon Eocaiman as a 'basal' alligatoroid; and (3) the paraphyly of the Cenozoic European taxon Diplocynodon. Among crocodyloids, notable results include modifications to the taxonomic content of Mekosuchinae, including biogeographic affinities of this clade with latest Cretaceous-early Paleogene Asian crocodyloids. In light of our new results, we provide a comprehensive review of the evolutionary and biogeographic history of Crocodylia, which included multiple instances of transoceanic and continental dispersal.

Near-complete phylogeny of extant Crocodylia (Reptilia) using mitogenome-based data

Species of the order Crocodylia are mostly large, predatory and semi-aquatic reptiles. Crocodylia, the closest living relatives of birds, first appeared in the Late Cretaceous period. In the present study, the complete mitochondrial (mt) genomes of 19 Crocodylia species, including two species (Melanosuchus niger and Caiman yacare) that have not been previously sequenced for mitogenomes, were processed through Illumina sequencing to offer genetic resources and compare with the mitogenomes of Crocodylia species reported previously. In addition, a high-resolution phylogenetic tree of nearly all current recognized species of Crocodylia is constructed based on mitogenomic data. Phylogenetic analyses support monophyly of three families: Alligatoridae (four genera: Alligator, Caiman, Melanosuchus and Paleosuchus), Crocodylidae (three genera: Crocodylus, Mecistops and Osteolaemus) and Gavialidae (two genera: Gavialis and Tomistoma). The tree topology is generally similar to previous studies. Molecular dating suggests that the first split within Crocodylia date back to the Upper Cretaceous (approx. 86.75 Mya). The estimated time to the most recent common ancestor (TMRCA) of Alligatoridae is 53.33 Mya and that of Crocodylidae and Gavialidae is 50.13 Mya, which might be closely linked to climate changes during the Late Palaeocene and Early Eocene. Additionally, this study proves that the diversification rate within Crocodylia began to increase from the Late Eocene (about 36 Mya) and two diversification peak periods of Crocodylia (0–10 Mya and 10–20 Mya) are disclosed, which is roughly consistent with the estimated crocodylian species richness through time. Combining all these clues, we can suggest that climate fluctuation may have played a decisive role in the speciation of Crocodylia.

Reptile enamel matrix proteins: Selection, divergence, and functional constraint

The three major enamel matrix proteins (EMPs): amelogenin (AMEL), ameloblastin (AMBN), and enamelin (ENAM), are intrinsically linked to tooth development in tetrapods. However, reptiles and mammals exhibit significant differences in dental patterning and development, potentially affecting how EMPs evolve in each group. In most reptiles, teeth are replaced continuously throughout life, while mammals have reduced replacement to only one or two generations. Reptiles also form structurally simple, aprismatic enamel while mammalian enamel is composed of highly organized hydroxyapatite prisms. These differences, combined with reported low sequence homology in reptiles, led us to predict that reptiles may experience lower selection pressure on their EMPs as compared with mammals. However, we found that like mammals, reptile EMPs are under moderate purifying selection, with some differences evident between AMEL, AMBN, and ENAM. We also demonstrate that sequence homology in reptile EMPs is closely associated with divergence times, with more recently diverged lineages exhibiting high homology, along with strong phylogenetic signal. Lastly, despite sequence divergence, none of the reptile species in our study exhibited mutations consistent with diseases that cause degeneration of enamel (e.g. amelogenesis imperfecta). Despite short tooth retention time and simplicity in enamel structure, reptile EMPs still exhibit purifying selection required to form durable enamel.

Genetic evidence supports a distinct lineage of American crocodile (Crocodylus acutus) in the Greater Antilles

Four species of true crocodile (genus Crocodylus) have been described from the Americas. Three of these crocodile species exhibit non-overlapping distributions—Crocodylus intermedius in South America, C. moreletii along the Caribbean coast of Mesoamerica, and C. rhombifer confined to Cuba. The fourth, C. acutus, is narrowly sympatric with each of the other three species. In this study, we sampled 113 crocodiles across Crocodylus populations in Cuba, as well as exemplar populations in Belize and Florida (USA), and sequenced three regions of the mitochondrial genome (D-loop, cytochrome b, cytochrome oxidase I; 3,626 base pair long dataset) that overlapped with published data previously collected from Colombia, Jamaica, and the Cayman Islands. Phylogenetic analyses of these data revealed two, paraphyletic lineages of C. acutus. One lineage, found in the continental Americas, is the sister taxon to C. intermedius, while the Greater Antillean lineage is most closely related to C. rhombifer. In addition to the paraphyly of the two C. acutus lineages, we recovered a 5.4% estimate of Tamura-Nei genetic divergence between the Antillean and continental clades. The reconstructed paraphyly, distinct phylogenetic affinities and high genetic divergence between Antillean and continental C. acutus populations are consistent with interspecific differentiation within the genus and suggest that the current taxon recognized as C. acutus is more likely a complex of cryptic species warranting a reassessment of current taxonomy. Moreover, the inclusion, for the first time, of samples from the western population of the American crocodile in Cuba revealed evidence for continental mtDNA haplotypes in the Antilles, suggesting this area may constitute a transition zone between distinct lineages of C. acutus. Further study using nuclear character data is warranted to more fully characterize this cryptic diversity, resolve taxonomic uncertainty, and inform conservation planning in this system.

Vertebral morphometrics and lung structure in non-avian dinosaurs

The lung-air sac system of modern birds is unique among vertebrates. However, debate surrounds whether an avian-style lung is restricted to birds or first appeared in their dinosaurian ancestors, as common osteological correlates for the respiratory system offer limited information on the lungs themselves. Here, we shed light on these issues by using axial morphology as a direct osteological correlate of lung structure, and quantifying vertebral shape using geometric morphometrics in birds, crocodilians and a wide range of dinosaurian taxa. Although fully avian lungs were a rather late innovation, we quantitatively show that non-avian dinosaurs and basal dinosauriforms possessed bird-like costovertebral joints and a furrowed thoracic ceiling. This would have immobilized the lung's dorsal surface, a structural prerequisite for a thinned blood-gas barrier and increased gas exchange potential. This could have permitted high levels of aerobic and metabolic activity in dinosaurs, even in the hypoxic conditions of the Mesozoic, contributing to their successful radiation.

The complete mitochondrial genomes of Tarsiger cyanurus and Phoenicurus auroreus: a phylogenetic analysis of Passeriformes

Passeriformes is the largest group within aves and the phylogenetic relationships between Passeriformes have caused major disagreement in ornithology. Particularly, the phylogenetic relationships between muscicapoidea and sylvioidea are complex, and their taxonomic boundaries have not been clearly defined. Our aim was to study the status of two bird species: Tarsiger cyanurus and Phoenicurus auroreus. Furthermore, we analyzed the phylogenetic relationships of Passeriformes. Complete mitochondrial DNA (mtDNA) sequences of both species were determined and the lengths were 16,803 (T. cyanurus) and 16,772 bp (P. auroreus), respectively. Thirteen protein-coding genes, 22 tRNA genes, two rRNA genes, and one control region were identified in these mtDNAs. The contents of A and T at the base compositions was significantly higher than the content of G and C, and this AT skew was positive, while the GC skew was negative. The monophyly of Passeriformes is divided into four major clades: Corvoidea, Sylvioidea, Passeroidea, and Musicicapoidea. Paridae should be separated from the superfamily Sylvioidea and placed within the superfamily Muscicapoidea. The family Muscicapidae and Corvida were paraphyly, while Carduelis and Emberiza were grouped as a sister taxon. The relationships between some species of the order passeriformes may remain difficult to resolve despite an effort to collect additional characters for phylogenetic analysis. Current research of avian phylogeny should focus on adding characters and taxa and use both effectively to obtain a better resolution for deeper and shallow nodes.

Role of Chromosome Changes in Crocodylus Evolution and Diversity

The karyotypes of most species of crocodilians were studied using conventional and molecular cytogenetics. These provided an important contribution of chromosomal rearrangements for the evolutionary processes of Crocodylia and Sauropsida (birds and reptiles). The karyotypic features of crocodilians contain small diploid chromosome numbers (30~42), with little interspecific variation of the chromosome arm number (fundamental number) among crocodiles (56~60). This suggested that centric fusion and/or fission events occurred in the lineage, leading to crocodilian evolution and diversity. The chromosome numbers of Alligator, Caiman, Melanosuchus, Paleosuchus, Gavialis, Tomistoma, Mecistops, and Osteolaemus were stable within each genus, whereas those of Crocodylus (crocodylians) varied within the taxa. This agreed with molecular phylogeny that suggested a highly recent radiation of Crocodylus species. Karyotype analysis also suggests the direction of molecular phylogenetic placement among Crocodylus species and their migration from the Indo-Pacific to Africa and The New World. Crocodylus species originated from an ancestor in the Indo-Pacific around 9~16 million years ago (MYA) in the mid-Miocene, with a rapid radiation and dispersion into Africa 8~12 MYA. This was followed by a trans-Atlantic dispersion to the New World between 4~8 MYA in the Pliocene. The chromosomes provided a better understanding of crocodilian evolution and diversity, which will be useful for further study of the genome evolution in Crocodylia.

Biased gene conversion and GC-content evolution in the coding sequences of reptiles and vertebrates

Mammalian and avian genomes are characterized by a substantial spatial heterogeneity of GC-content, which is often interpreted as reflecting the effect of local GC-biased gene conversion (gBGC), a meiotic repair bias that favors G and C over A and T alleles in high-recombining genomic regions. Surprisingly, the first fully sequenced nonavian sauropsid (i.e., reptile), the green anole Anolis carolinensis, revealed a highly homogeneous genomic GC-content landscape, suggesting the possibility that gBGC might not be at work in this lineage. Here, we analyze GC-content evolution at third-codon positions (GC3) in 44 vertebrates species, including eight newly sequenced transcriptomes, with a specific focus on nonavian sauropsids. We report that reptiles, including the green anole, have a genome-wide distribution of GC3 similar to that of mammals and birds, and we infer a strong GC3-heterogeneity to be already present in the tetrapod ancestor. We further show that the dynamic of coding sequence GC-content is largely governed by karyotypic features in vertebrates, notably in the green anole, in agreement with the gBGC hypothesis. The discrepancy between third-codon positions and noncoding DNA regarding GC-content dynamics in the green anole could not be explained by the activity of transposable elements or selection on codon usage. This analysis highlights the unique value of third-codon positions as an insertion/deletion-free marker of nucleotide substitution biases that ultimately affect the evolution of proteins.

An expanded combined evidence approach to the Gavialis problem using geometric morphometric data from crocodylian braincases and Eustachian systems

The phylogenetic position of the Indian gharial (Gavialis gangeticus) is disputed--morphological characters place Gavialis as the sister to all other extant crocodylians, whereas molecular and combined analyses find Gavialis and the false gharial (Tomistoma schlegelii) to be sister taxa. Geometric morphometric techniques have only begun to be applied to this issue, but most of these studies have focused on the exterior of the skull. The braincase has provided useful phylogenetic information for basal crurotarsans, but has not been explored for the crown group. The Eustachian system is thought to vary phylogenetically in Crocodylia, but has not been analytically tested. To determine if gross morphology of the crocodylian braincase proves informative to the relationships of Gavialis and Tomistoma, we used two- and three-dimensional geometric morphometric approaches. Internal braincase images were obtained using high-resolution computerized tomography scans. A principal components analysis identified that the first component axis was primarily associated with size and did not show groupings that divide the specimens by phylogenetic affinity. Sliding semi-landmarks and a relative warp analysis indicate that a unique Eustachian morphology separates Gavialis from other extant members of Crocodylia. Ontogenetic expansion of the braincase results in a more dorsoventrally elongate median Eustachian canal. Changes in the shape of the Eustachian system do provide phylogenetic distinctions between major crocodylian clades. Each morphometric dataset, consisting of continuous morphological characters, was added independently to a combined cladistic analysis of discrete morphological and molecular characters. The braincase data alone produced a clade that included crocodylids and Gavialis, whereas the Eustachian data resulted in Gavialis being considered a basally divergent lineage. When each morphometric dataset was used in a combined analysis with discrete morphological and molecular characters, it generated a tree that matched the topology of the molecular phylogeny of Crocodylia.

The complete mitochondrial genome of Amolops ricketti (Amphidia, Anura, Ranidae)

The complete mitochondrial genome from the South China torrent frog Amolops ricketti was determined. This mitogenome was 17,771 bp in length, containing 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes and a control region (CR). All the protein-coding genes in A. ricketti were distributed on the H-strand, except for the ND6 subunit gene and eight tRNA genes which were encoded on the L-strand.

Complete mitogenome of Asiatic lion resolves phylogenetic status within Panthera

BACKGROUND

The origin, evolution and speciation of the lion, has been subject of interest, debate and study. The present surviving lions of the genus Panthera comprise of eight sub-species inclusive of Asiatic lion Panthera leo persica of India's Gir forest. Except for the Asiatic lion, the other seven subspecies are found in different parts of Africa. There have been different opinions regarding the phylogenetic status of Panthera leo, as well as classifying lions of different geographic regions into subspecies and races. In the present study, mitogenome sequence of P. leo persica deduced, using Ion Torrent PGM to assess phylogeny and evolution which may play an increasingly important role in conservation biology.

RESULTS

The mtDNA sequence of P. leo persica is 17,057 bp in length with 40.8% GC content. Annotation of mitogenome revealed total 37 genes, including 13 protein coding, 2 rRNA and 22 tRNA. Phylogenetic analysis based on whole mitogenome, suggests Panthera pardus as a neighbouring species to P. leo with species divergence at ~2.96 mya.

CONCLUSION

This work presents first report on complete mitogenome of Panthera leo persica. It sheds light on the phylogenetic and evolutionary status within and across Felidae members. The result compared and evaluated with earlier reports of Felidae shows alteration of phylogenetic status and species evolution. This study may provide information on genetic diversity and population stability.

Global taxonomic diversity of living reptiles

Reptiles are one of the most ecologically and evolutionarily remarkable groups of living organisms, having successfully colonized most of the planet, including the oceans and some of the harshest and more environmentally unstable ecosystems on earth. Here, based on a complete dataset of all the world's diversity of living reptiles, we analyse lineage taxonomic richness both within and among clades, at different levels of the phylogenetic hierarchy. We also analyse the historical tendencies in the descriptions of new reptile species from Linnaeus to March 2012. Although (non-avian) reptiles are the second most species-rich group of amniotes after birds, most of their diversity (96.3%) is concentrated in squamates (59% lizards, 35% snakes, and 2% amphisbaenians). In strong contrast, turtles (3.4%), crocodilians (0.3%), and tuataras (0.01%) are far less diverse. In terms of species discoveries, most turtles and crocodilians were described early, while descriptions of lizards, snakes and amphisbaenians are multimodal with respect to time. Lizard descriptions, in particular, have reached unprecedented levels during the last decade. Finally, despite such remarkably asymmetric distributions of reptile taxonomic diversity among groups, we found that the distributions of lineage richness are consistently right-skewed, with most clades (monophyletic families and genera) containing few lineages (monophyletic genera and species, respectively), while only a few have radiated greatly (notably the families Colubridae and Scincidae, and the lizard genera Anolis and Liolaemus). Therefore, such consistency in the frequency distribution of richness among clades and among phylogenetic levels suggests that the nature of reptile biodiversity is fundamentally fractal (i.e., it is scale invariant). We then compared current reptile diversity with the global reptile diversity and taxonomy known in 1980. Despite substantial differences in the taxonomies (relative to 2012), the patterns of lineage richness remain qualitatively identical, hence reinforcing our conclusions about the fractal nature of reptile biodiversity.

Perspectives on the convergent evolution of tetrapod salt glands

Since their discovery in 1958, the function of specialized salt-secreting glands in tetrapods has been studied in great detail, and such studies continue to contribute to a general understanding of transport mechanisms of epithelial water and ions. Interestingly, during that same time period, there have been only few attempts to understand the convergent evolution of this tissue, likely as a result of the paucity of taxonomic, embryological, and molecular data available. In this review, we synthesize the available data regarding the distribution of salt glands across extant and extinct tetrapod lineages and the anatomical position of the salt gland in each taxon. Further, we use these data to develop hypotheses about the various factors that have influenced the convergent evolution of salt glands across taxa with special focus on the variation in the anatomical position of the glands and on the molecular mechanisms that may have facilitated the development of a salt gland by co-option of a nonsalt-secreting ancestral gland. It is our hope that this review will stimulate renewed interest in the topic of the convergent evolution of salt glands and inspire future empirical studies aimed at evaluating the hypotheses we lay out herein.

Insights into the ecology and evolutionary success of crocodilians revealed through bite-force and tooth-pressure experimentation

Background

Crocodilians have dominated predatory niches at the water-land interface for over 85 million years. Like their ancestors, living species show substantial variation in their jaw proportions, dental form and body size. These differences are often assumed to reflect anatomical specialization related to feeding and niche occupation, but quantified data are scant. How these factors relate to biomechanical performance during feeding and their relevance to crocodilian evolutionary success are not known.

Methodology/Principal Findings

We measured adult bite forces and tooth pressures in all 23 extant crocodilian species and analyzed the results in ecological and phylogenetic contexts. We demonstrate that these reptiles generate the highest bite forces and tooth pressures known for any living animals. Bite forces strongly correlate with body size, and size changes are a major mechanism of feeding evolution in this group. Jaw shape demonstrates surprisingly little correlation to bite force and pressures. Bite forces can now be predicted in fossil crocodilians using the regression equations generated in this research.

Conclusions/Significance

Critical to crocodilian long-term success was the evolution of a high bite-force generating musculo-skeletal architecture. Once achieved, the relative force capacities of this system went essentially unmodified throughout subsequent diversification. Rampant changes in body size and concurrent changes in bite force served as a mechanism to allow access to differing prey types and sizes. Further access to the diversity of near-shore prey was gained primarily through changes in tooth pressure via the evolution of dental form and distributions of the teeth within the jaws. Rostral proportions changed substantially throughout crocodilian evolution, but not in correspondence with bite forces. The biomechanical and ecological ramifications of such changes need further examination.

Complete mitochondrial genome sequence from an endangered Indian snake, Python molurus molurus (Serpentes, Pythonidae)

This paper reports the complete mitochondrial genome sequence of an endangered Indian snake, Python molurus molurus (Indian Rock Python). A typical snake mitochondrial (mt) genome of 17258 bp length comprising of 37 genes including the 13 protein coding genes, 22 tRNA genes, and 2 ribosomal RNA genes along with duplicate control regions is described herein. The P. molurus molurus mt. genome is relatively similar to other snake mt. genomes with respect to gene arrangement, composition, tRNA structures and skews of AT/GC bases. The nucleotide composition of the genome shows that there are more A-C % than T-G% on the positive strand as revealed by positive AT and CG skews. Comparison of individual protein coding genes, with other snake genomes suggests that ATP8 and NADH3 genes have high divergence rates. Codon usage analysis reveals a preference of NNC codons over NNG codons in the mt. genome of P. molurus. Also, the synonymous and non-synonymous substitution rates (ka/ks) suggest that most of the protein coding genes are under purifying selection pressure. The phylogenetic analyses involving the concatenated 13 protein coding genes of P. molurus molurus conformed to the previously established snake phylogeny.