Dating cryptodiran nodes: origin and diversification of the turtle superfamily Testudinoidea

Turtle body size evolution is determined by lineage-specific specializations rather than global trends

Organisms display a considerable variety of body sizes and shapes, and macroevolutionary investigations help to understand the evolutionary dynamics behind such variations. Turtles (Testudinata) show great body size disparity, especially when their rich fossil record is accounted for. We explored body size evolution in turtles, testing which factors might influence the observed patterns and evaluating the existence of long-term directional trends. We constructed the most comprehensive body size dataset for the group to date, tested for correlation with paleotemperature, estimated ancestral body sizes, and performed macroevolutionary model-fitting analyses. We found no evidence for directional body size evolution, even when using very flexible models, thereby rejecting the occurrence of Cope's rule. We also found no significant effect of paleotemperature on overall through-time body size patterns. In contrast, we found a significant influence of habitat preference on turtle body size. Freshwater turtles display a rather homogeneous body size distribution through time. In contrast, terrestrial and marine turtles show more pronounced variation, with terrestrial forms being restricted to larger body sizes, up to the origin of testudinids in the Cenozoic, and marine turtles undergoing a reduction in body size disparity after the extinctions of many groups in the mid-Cenozoic. Our results, therefore, suggest that long-term, generalized patterns are probably explained by factors specific to certain groups and related at least partly to habitat use.

Phylogenomics reconciles molecular data with the rich fossil record on the origin of living turtles

Although a consensus exists that all living turtles fall within either Pleurodira or Cryptodira clades, estimating when these lineages split is still under debate. Most molecular studies date the split in the Triassic Period, whereas a Jurassic age is unanimous among morphological studies. Each hypothesis implies different paleobiogeographical scenarios to explain early turtle evolution. Here we explored the rich turtle fossil record with the Fossilized Birth-Death (FBD) and the traditional node dating (ND) methods using complete mitochondrial genomes (147 taxa) and a set of nuclear orthologs with over 10 million bp (25 taxa) to date the major splits in Testudines. Our results support an Early Jurassic split (191-182 Ma) for the crown Testudines with great consistency across different dating methods and datasets, with a narrow confidence interval. This result is independently supported by the oldest fossils of Testudines that postdate the Middle Jurassic (174 Ma), which were not used for calibration in this study. This age coincides with the Pangaea fragmentation and the formation of saltwater barriers such as the Atlantic Ocean and the Turgai Strait, supporting that diversification in Testudines was triggered by vicariance. Our ages of the splits in Pleurodira coincide with the geologic events of the Late Jurassic and Early Cretaceous. Conversely, the early Cryptodira radiation remained in Laurasia, and its diversification ensued as all its major lineages expanded their distribution into every continent during the Cenozoic. We provide the first detailed hypothesis of the evolution of Cryptodira in the Southern Hemisphere, in which our time estimates are correlated with each contact between landmasses derived from Gondwana and Laurasia. Although most South American Cryptodira arrived through the Great American Biotic Interchange, our results indicate that the Chelonoidis ancestor probably arrived from Africa through the chain islands of the South Atlantic during the Paleogene. Together, the presence of ancient turtle diversity and the vital role that turtles occupy in marine and terrestrial ecosystems underline South America as a chief area for conservation.

Geometric morphometrics and anatomical network analyses reveal ecospace partitioning among geoemydid turtles from the Uinta Formation, Utah

We present new fossil records of the geoemydid turtle Bridgeremys pusilla from the Uinta Formation of Utah. Turtles are abundant throughout the unit, and known taxa are similar to those from the older strata in the Upper Green River Basin in Wyoming from the Bridger and Washakie Formations. B. pusilla is known from Bridgerian deposits but was not previously known from after the Turtle Bluff Member of the Bridger Formation. The taxon was coveal with two species of the geoemydid Echmatemys (E. callopyge and E. wyomingensis), a common genus of extinct pond turtles known primarily from lacustrine and fluvial deposits in western North America, including the Uinta Basin. In addition to previously documented morphological differences, our geometric morphometric analyses revealed significant differences in epiplastral morphology between B. pusilla and the two coeval Echmatemys species. Bridgeremys pusilla shared several morphological characters with Testudinidae. However, our anatomical network analysis suggests that the carapace of B. pusilla distributed stress forces in a manner more similar to emydids (basal and derived) than to derived testudinoids (Testudinidae and Emydidae), including Echmatemys species. This finding changes our understanding of the ecology of the species and sheds light onto how geoemydid turtles of the Uinta Formation may have partitioned the available ecospace. These new Uintan records extend the geographic range of B. pusilla into the Uinta Basin and stratigraphically through the top of the Uinta Formation, extending the temporal range of the taxon by more than 4 million years through the Uintan North American Land Mammal Age to the base of the Duchesne River Formation.

The phylogenetic relationships of geoemydid turtles from the Eocene Messel Pit Quarry: a first assessment using methods for continuous and discrete characters

The geoemydid turtles of the Eocoene Messel Pit Quarry of Hesse, Germany, are part of a rich Western European fossil record of testudinoids. Originally referred to as "Ocadia" kehreri and "Ocadia" messeliana, their systematic relationships remain unclear. A previous study proposed that a majority of the Western European geoemydids, including the Messel geoemydids, are closely related to the Recent European representatives of the clade Mauremys. Another study hypothesised that the Western European geoemydid fauna is more phylogenetically diverse, and that the Messel geoemydids are closely related to the East Asian turtles Orlitia and Malayemys. Here we present the first quantitative analyses to date that investigate this question. We use continuous characters in the form of ratios to estimate the placement of the Messel geoemydids in a reference tree that was estimated from molecular data. We explore the placement error obtained from that data with maximum likelihood and Bayesian methods, as well as linear parsimony in combination with discrete characters. We find good overall performance with Bayesian and parsimony analyses. Parsimony performs even better when we also incorporated discrete characters. Yet, we cannot pin down the position of the Messel geoemydids with high confidence. Depending on how intraspecific variation of the ratio characters is treated, parsimony favours a placement of the Messel fossils sister to Orlitia borneensis or sister to Geoemyda spengleri, with weak bootstrap support. The latter placement is suspect because G. spengleri is a phylogenetically problematic species with molecular and morphological data. There is even less support for placements within the Mauremys clade.

Redescription of the skull of the Australian flatback sea turtle, Natator depressus, provides new morphological evidence for phylogenetic relationships among sea turtles (Chelonioidea)

Chelonioidea (sea turtles) are a group where available morphological evidence for crown-group relationships are incongruent with those established using molecular data. However, morphological surveys of crown-group taxa tend to focus on a recurring subset of the extant species. The Australian flatback sea turtle, Natator depressus, is often excluded from comparisons and it is the most poorly known of the seven extant species of Chelonioidea. Previous descriptions of its skull morphology are limited and conflict. Here we describe three skulls of adult N. depressus and re-examine the phylogenetic relationships according to morphological character data. Using X-ray micro Computed Tomography we describe internal structures of the braincase and identify new phylogenetically informative characters not previously reported. Phylogenetic analysis using a Bayesian approach strongly supports a sister-group relationship between Chelonia mydas and N. depressus, a topology that was not supported by previous analyses of morphological data but one that matches the topology supported by analysis of molecular data. Our results highlight the general need to sample the morphological anatomy of crown-group taxa more thoroughly before concluding that morphological and molecular evidence are incongruous.

Climatic and topographic changes since the Miocene influenced the diversification and biogeography of the tent tortoise (Psammobates tentorius) species complex in Southern Africa

BACKGROUND

Climatic and topographic changes function as key drivers in shaping genetic structure and cladogenic radiation in many organisms. Southern Africa has an exceptionally diverse tortoise fauna, harbouring one-third of the world's tortoise genera. The distribution of Psammobates tentorius (Kuhl, 1820) covers two of the 25 biodiversity hotspots in the world, the Succulent Karoo and Cape Floristic Region. The highly diverged P. tentorius represents an excellent model species for exploring biogeographic and radiation patterns of reptiles in Southern Africa.

RESULTS

We investigated genetic structure and radiation patterns against temporal and spatial dimensions since the Miocene in the Psammobates tentorius species complex, using multiple types of DNA markers and niche modelling analyses. Cladogenesis in P. tentorius started in the late Miocene (11.63-5.33 Ma) when populations dispersed from north to south to form two geographically isolated groups. The northern group diverged into a clade north of the Orange River (OR), followed by the splitting of the group south of the OR into a western and an interior clade. The latter divergence corresponded to the intensification of the cold Benguela current, which caused western aridification and rainfall seasonality. In the south, tectonic uplift and subsequent exhumation, together with climatic fluctuations seemed responsible for radiations among the four southern clades since the late Miocene. We found that each clade occurred in a habitat shaped by different climatic parameters, and that the niches differed substantially among the clades of the northern group but were similar among clades of the southern group.

CONCLUSION

Climatic shifts, and biome and geographic changes were possibly the three major driving forces shaping cladogenesis and genetic structure in Southern African tortoise species. Our results revealed that the cladogenesis of the P. tentorius species complex was probably shaped by environmental cooling, biome shifts and topographic uplift in Southern Africa since the late Miocene. The Last Glacial Maximum (LGM) may have impacted the distribution of P. tentorius substantially. We found the taxonomic diversify of the P. tentorius species complex to be highest in the Greater Cape Floristic Region. All seven clades discovered warrant conservation attention, particularly Ptt-B-Ptr, Ptt-A and Pv-A.

The complete mitochondrial genome of the endangered Assam Roofed Turtle, Pangshura sylhetensis (Testudines: Geoemydidae): Genomic features and phylogeny

The Assam Roofed Turtle, Pangshura sylhetensis is an endangered and least studied species endemic to India and Bangladesh. The present study decodes the first complete mitochondrial genome of P. sylhetensis (16,568 bp) by using next-generation sequencing. The assembly encodes 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and one control region (CR). Most of the genes were encoded on the majority strand, except NADH dehydrogenase subunit 6 (nad6) and eight tRNAs. All PCGs start with an ATG initiation codon, except for Cytochrome oxidase subunit 1 (cox1) and NADH dehydrogenase subunit 5 (nad5), which both start with GTG codon. The study also found the typical cloverleaf secondary structures in most of the predicted tRNA structures, except for serine (trnS1) which lacks of conventional DHU arm and loop. Both Bayesian and maximum-likelihood phylogenetic inference using 13 concatenated PCGs demonstrated strong support for the monophyly of all 52 Testudines species within their respective families and revealed Batagur trivittata as the nearest neighbor of P. sylhetensis. The mitogenomic phylogeny with other amniotes is congruent with previous research, supporting the sister relationship of Testudines and Archosaurians (birds and crocodilians). Additionally, the mitochondrial Gene Order (GO) analysis indicated plesiomorphy with the typical vertebrate GO in most of the Testudines species.

The earliest known crown-Testudo tortoise from the late Miocene (Vallesian, 9 Ma) of Greece

We here report on fossil remains of the earliest known crown-Testudo, an extant clade of Mediterranean testudinid tortoises from the late Miocene (Vallesian, MN 10) from the hominoid locality Ravin de la Pluie (RPl) in Greece. The material studied is a small, nearly complete carapace with a clearly distinct hypo-xiphiplastral hinge. This supports the sensu stricto generic assignment. This new terrestrial testudinid specimen is characterized by a possible tectiform, narrow, elongated shell with a pentagonal pygal and a long, posteriorly elevated, lenticular and rounded dorsal epiplastral lip. These unique features differ from those of other known Mediterranean hinged forms and allow the erection of the new species Testudo hellenica sp. nov. This taxon is phylogenetically close to two Greek species, the extant T. marginata and the fossil T. marmorum (Turolian, around 7.3 Ma). This record provides evidence for the first appearance of the genus Testudo sensu stricto at a minimum age of 9 Ma.

Sex is determined by XX/XY sex chromosomes in Australasian side-necked turtles (Testudines: Chelidae)

Turtles demonstrate variability in sex determination and, hence, constitute an excellent model for the evolution of sex chromosomes. Notably, the sex determination of the freshwater turtles from the family Chelidae, a species-rich group with wide geographical distribution in the southern hemisphere, is still poorly explored. Here we documented the presence of an XX/XY sex determination system in seven species of the Australasian chelid genera Chelodina, Emydura, and Elseya by conventional (karyogram reconstruction, C-banding) and molecular cytogenetic methods (comparative genome hybridization, in situ hybridization with probes specific for GATA microsatellite motif, the rDNA loci, and the telomeric repeats). The sex chromosomes are microchromosomes in all examined species of the genus Chelodina. In contrast, the sex chromosomes are the 4^th largest pair of macrochromosomes in the genera Emydura and Elseya. Their X chromosomes are submetacentric, while their Y chromosomes are metacentric. The chelid Y chromosomes contain a substantial male-specific genomic region with an accumulation of the GATA microsatellite motif, and occasionally, of the rDNA loci and telomeric repeats. Despite morphological differences between sex chromosomes, we conclude that male heterogamety was likely already present in the common ancestor of Chelodina, Emydura and Elseya in the Mesozoic period.

Comparative analysis of the shape and size of the middle ear cavity of turtles reveals no correlation with habitat ecology

The middle ear of turtles differs from other reptiles in being separated into two distinct compartments. Several ideas have been proposed as to why the middle ear is compartmentalized in turtles, most suggesting a relationship with underwater hearing. Extant turtle species span fully marine to strictly terrestrial habitats, and ecomorphological hypotheses of turtle hearing predict that this should correlate with variation in the structure of the middle ear due to differences in the fluid properties of water and air. We investigate the shape and size of the air-filled middle ear cavity of 56 extant turtles using 3D data and phylogenetic comparative analysis to test for correlations between habitat preferences and the shape and size of the middle ear cavity. Only weak correlations are found between middle ear cavity size and ecology, with aquatic taxa having proportionally smaller cavity volumes. The middle ear cavity of turtles exhibits high shape diversity among species, but we found no relationship between this shape variation and ecology. Surprisingly, the estimated acoustic transformer ratio, a key functional parameter of impedance-matching ears in vertebrates, also shows no relation to habitat preferences (aquatic/terrestrial) in turtles. We suggest that middle ear cavity shape may be controlled by factors unrelated to hearing, such as the spatial demands of surrounding cranial structures. A review of the fossil record suggests that the modern turtle ear evolved during the Early to Middle Jurassic in stem turtles broadly adapted to freshwater and terrestrial settings. This, combined with our finding that evolutionary transitions between habitats caused only weak evolutionary changes in middle ear structure, suggests that tympanic hearing in turtles evolved as a compromise between subaerial and underwater hearing.

The Complete Mitochondrial Genome of Platysternon megacephalum peguense and Molecular Phylogenetic Analysis

Platysternon megacephalum is the only living representative species of Platysternidae and only three subspecies remain: P. m. megalorcephalum, P. m. shiui, and P. m. peguense. However, previous reports implied that P. m. peguense has distinct morphological and molecular features. The characterization of the mitogenome has been accepted as an efficient means of phylogenetic and evolutionary analysis. Hence, this study first determined the complete mitogenome of P. m. peguense with the aim to identify the structure and variability of the P. m. peguense mitogenome through comparative analysis. Furthermore, the phylogenetic relationship of the three subspecies was tested. Based on different tRNA gene loss and degeneration of these three subspecies, their rearrangement pathways have been inferred. Phylogenetic analysis showed that P. m. peguense is a sister group to (P. m. megalorcephalum and P. m. shiui). Furthermore, the divergence time estimation of these three subspecies coincided with the uplift of the Tibetan Plateau. This study shows that the genetic distances between P. m. peguense and the other two subspecies are comparable to interspecific genetic distances, for example within Mauremys. In general, this study provides new and meaningful insights into the evolution of the three Platysternidae subspecies.

A new testudinoid turtle from the middle to late Eocene of Vietnam

BACKGROUND

Testudinoidea is a major clade of turtles that has colonized different ecological environments across the globe throughout the Tertiary. Aquatic testudinoids have a particularly rich fossil record in the Tertiary of the northern hemisphere, but little is known about the evolutionary history of the group, as the phylogenetic relationships of most fossils have not been established with confidence, in part due to high levels of homoplasy and polymorphism.

METHODS

We here focus on describing a sample of 30 testudinoid shells, belonging to a single population that was collected from lake sediments from the middle to late Eocene (35-39 Ma) Na Duong Formation in Vietnam. The phylogenetic placement of this new material is investigated by integrating it and 11 other species of putative geoemydids from the Eocene and Oligocene to a recently published matrix of geoemydid turtles, that embraces the use of polymorphic characters, and then running a total-evidence analysis.

RESULTS

The new material is highly polymorphic, but can be inferred with confidence to be a new taxon, Banhxeochelys trani gen. et sp. nov. It shares morphological similarities with other southeastern Asian testudinoids, Isometremys lacuna and Guangdongemys pingi, but is placed phylogenetically at the base of Pan-Testuguria when fossils are included in the analysis, or as a stem geoemydid when other fossils are deactivated from the matrix. The vast majority of other putative fossil geoemydids are placed at the base of Pan-Testuguria as well.

DISCUSSION

The phylogenetic placement of fossil testudinoids used in the analysis is discussed individually and each species compared to Banhxeochelys trani gen. et sp. nov. The high levels of polymorphism observed in the new taxon is discussed in terms of ontogenetic and random variability. This is the first time that a large sample of fossil testudinoids has its morphological variation described in detail.

Total evidence analysis and body size evolution of extant and extinct tortoises (Testudines: Cryptodira: Pan-Testudinidae)

Testudinidae (tortoises) is an extant clade of terrestrial turtles of worldwide distribution and with a rich fossil record that provides an exceptional context for studying their evolutionary history. Because of the lack of global phylogenetic analyses integrating extinct taxa, our current knowledge of the relationships of the total clade of Testudinidae is rather poor. To resolve this issue, we performed the first total evidence analysis of Pan-Testudinidae. The total evidence trees are congruent with the molecular topology and agree on the dichotomy of derived Testudinidae (=Testudininae; Converted Clade Name) into two previously recognized major clades, Testudona and Geochelona (New Clade Name). The integration of extinct taxa into the analysis allowed the stratigraphic fit of the total evidence trees, indicating that crown Testudininae, Testudona and Geochelona all originated by the Late Eocene, in agreement with recent molecular estimates. Ghost lineage analysis indicates high diversification in the Late Eocene and in the Miocene. The age of crown Testudo is Late Miocene, again in accordance with some molecular dates. Phylogenetic placement of fossils demonstrates that giant body size independently evolved in multiple continental mainland taxa and confirms recent results deduced from living taxa-giantism in Testudinidae is not linked to the insular effect. An unexpected outcome is the recovery of miniaturization in Testudona (<30 cm carapace length) that emerged sometime between the Oligocene and Early Miocene. No clear correlation between body size evolution and climate is apparent, but increased taxon sampling may nevertheless demonstrate the role of cooling and warming as one of many influential variables.

A new slider turtle (Testudines: Emydidae: Deirochelyinae: Trachemys) from the late Hemphillian (late Miocene/early Pliocene) of eastern Tennessee and the evolution of the deirochelyines

Trachemys (Testudines: Emydidae) represents one of the most well-known turtle genera today. The evolution of Trachemys, while being heavily documented with fossil representatives, is not well understood. Numerous fossils from the late Hemphillian Gray Fossil Site (GFS) in northeastern Tennessee help to elucidate its evolution. The fossil Trachemys at the GFS represent a new species. The new taxon, Trachemys haugrudi, is described, and currently represents the most thoroughly described fossil emydid species known. A phylogenetic analysis, including 31 species, focusing on the subfamily Deirochelyinae is performed that includes the new fossil species, along with numerous other modern and fossil deirochelyine species, representing the first phylogenetic analysis published that includes several fossil deirochelyines. The phylogenetic analysis, utilizing morphological evidence, provides monophyletic clades of all modern deirochelyines, including Chrysemys, Deirochelys, Pseudemys, Malaclemys, Graptemys, and Trachemys. A strict consensus tree finds the recently described fossil species Graptemys kerneri to be part of a clade of Graptemys + Malaclemys. Three fossil taxa, including one previously referred to Pseudemys (Pseudemys caelata) and two to Deirochelys (Deirochelys carri and Deirochelys floridana) are found to form a clade with modern Deirochelys reticularia reticularia, with D. floridana sister to the other members of the clade. Chrysemys is found to be part of a basal polytomy with Deirochelys in relation to other deirochelyine taxa. Two fossil taxa previously referred to Chrysemys (Chrysemys timida and Chrysemys williamsi) form a paraphyly with the modern Chrysemys picta picta and Deirochelys, and may be referable to distinct genera. Additionally, fossil taxa previously attributed to Trachemys (Trachemys hillii, Trachemys idahoensis, Trachemys inflata, and Trachemys platymarginata) and T. haugrudi are found to form a clade separate from clades of northern and southern Trachemys species, potentially suggesting a distinct lineage of Trachemys with no modern survivors. Hypotheses of phylogenetic relationships mostly agree between the present study and previous ones, although the inclusion of fossil taxa provides further clues to the evolution of parts of the Deirochelyinae. The inclusion of more fossil taxa and characters may help resolve the placement of some taxa, and further elucidate the evolution of these New World turtles.

Molecular phylogeny and divergence of the map turtles (Emydidae: Graptemys)

The map turtles (genus Graptemys) comprise a morphologically diverse clade that forms a major component of the southeastern US hotspot of chelonian diversity. Map turtles have experienced both recent and rapid diversification resulting in long-standing uncertainty regarding species boundaries and phylogenetic relationships within the genus as well as timing of their divergence. We present a phylogeny for the group that includes geographically representative sampling for all described species and subspecies. We make use of an empirical prior on rates of molecular evolution to estimate divergence times with a molecular clock under a coalescent framework. Together, the phylogeny and divergence time estimates suggest that diversification has been both more recent and more rapid than has so far been suspected. We provide a well-supported evolutionary framework for Graptemys that is necessary for understanding map turtle diversity, biogeography, and for conservation of this threatened clade of turtles.

Phylotranscriptomic consolidation of the jawed vertebrate timetree

Phylogenomics is extremely powerful but introduces new challenges as no agreement exists on "standards" for data selection, curation and tree inference. We use jawed vertebrates (Gnathostomata) as model to address these issues. Despite considerable efforts in resolving their evolutionary history and macroevolution, few studies have included a full phylogenetic diversity of gnathostomes and some relationships remain controversial. We tested a novel bioinformatic pipeline to assemble large and accurate phylogenomic datasets from RNA sequencing and find this phylotranscriptomic approach successful and highly cost-effective. Increased sequencing effort up to ca. 10Gbp allows recovering more genes, but shallower sequencing (1.5Gbp) is sufficient to obtain thousands of full-length orthologous transcripts. We reconstruct a robust and strongly supported timetree of jawed vertebrates using 7,189 nuclear genes from 100 taxa, including 23 new transcriptomes from previously unsampled key species. Gene jackknifing of genomic data corroborates the robustness of our tree and allows calculating genome-wide divergence times by overcoming gene sampling bias. Mitochondrial genomes prove insufficient to resolve the deepest relationships because of limited signal and among-lineage rate heterogeneity. Our analyses emphasize the importance of large curated nuclear datasets to increase the accuracy of phylogenomics and provide a reference framework for the evolutionary history of jawed vertebrates.

Phylogenomic analyses of 539 highly informative loci dates a fully resolved time tree for the major clades of living turtles (Testudines)

Accurate time-calibrated phylogenies are the centerpiece of many macroevolutionary studies, and the relationship between the size and scale of molecular data sets and the density and accuracy of fossil calibrations is a key element of time tree studies. Here, we develop a target capture array specifically for living turtles, compare its efficiency to an ultraconserved element (UCE) dataset, and present a time-calibrated molecular phylogeny based on 539 nuclear loci sequenced from 26 species representing the breadth of living turtle diversity plus outgroups. Our gene array, based on three fully sequenced turtle genomes, is 2.4 times more variable across turtles than a recently published UCE data set for an identical subset of 13 species, confirming that taxon-specific arrays return more informative data per sequencing effort than UCEs. We used our genomic data to estimate the ages of living turtle clades including a mid-late Triassic origin for crown turtles and a mid-Carboniferous split of turtles from their sister group, Archosauria. By specifically excluding several of the earliest potential crown turtle fossils and limiting the age of fossil calibration points to the unambiguous crown lineage Caribemys oxfordiensis from the Late Jurassic (Oxfordian, 163.5-157.3Ma) we corroborate a relatively ancient age for living turtles. We also provide novel age estimates for five of the ten testudine families containing more than a single species, as well as several intrafamilial clades. Most of the diversity of crown turtles appears to date to the Paleogene, well after the Cretaceous-Paleogene mass extinction 66mya.

Phylogenetic relationships and divergence dates of softshell turtles (Testudines: Trionychidae) inferred from complete mitochondrial genomes

The softshell turtles (Trionychidae) are one of the most widely distributed reptile groups in the world, and fossils have been found on all continents except Antarctica. The phylogenetic relationships among members of this group have been previously studied; however, disagreements regarding its taxonomy, its phylogeography and divergence times are still poorly understood as well. Here, we present a comprehensive mitogenomic study of softshell turtles. We sequenced the complete mitochondrial genomes of 10 softshell turtles, in addition to the GenBank sequence of Dogania subplana, Lissemys punctata, Trionyx triunguis, which cover all extant genera within Trionychidae except for Cyclanorbis and Cycloderma. These data were combined with other mitogenomes of turtles for phylogenetic analyses. Divergence time calibration and ancestral reconstruction were calculated using BEAST and RASP software, respectively. Our phylogenetic analyses indicate that Trionychidae is the sister taxon of Carettochelyidae, and support the monophyly of Trionychinae and Cyclanorbinae, which is consistent with morphological data and molecular analysis. Our phylogenetic analyses have established a sister taxon relationship between the Asian Rafetus and the Asian Palea + Pelodiscus + Dogania + Nilssonia + Amyda, whereas a previous study grouped the Asian Rafetus with the American Apalone. The results of divergence time estimates and area ancestral reconstruction show that extant Trionychidae originated in Asia at around 108 million years ago (MA), and radiations mainly occurred during two warm periods, namely Late Cretaceous-Early Eocene and Oligocene. By combining the estimated divergence time and the reconstructed ancestral area of softshell turtles, we determined that the dispersal of softshell turtles out of Asia may have taken three routes. Furthermore, the times of dispersal seem to be in agreement with the time of the India-Asia collision and opening of the Bering Strait, which provide evidence for the accuracy of our estimation of divergence time. Overall, the mitogenomes of this group were used to explore the origin and dispersal route of Trionychidae and have provided new insights on the evolution of this group.

Stable Cretaceous sex chromosomes enable molecular sexing in softshell turtles (Testudines: Trionychidae)

Turtles demonstrate variability in sex determination ranging from environmental sex determination (ESD) to highly differentiated sex chromosomes. However, the evolutionary dynamics of sex determining systems in this group is not well known. Differentiated ZZ/ZW sex chromosomes were identified in two species of the softshell turtles (Trionychidae) from the subfamily Trionychinae and Z-specific genes were identified in a single species. We tested Z-specificity of a subset of these genes by quantitative PCR comparing copy gene numbers in male and female genomes in 10 species covering the phylogenetic diversity of trionychids. We demonstrated that differentiated ZZ/ZW sex chromosomes are conserved across the whole family and that they were already present in the common ancestor of the extant trionychids. As the sister lineage, Carettochelys insculpta, possess ESD, we can date the origin of the sex chromosomes in trionychids between 200 Mya (split of Trionychidae and Carettochelyidae) and 120 Mya (basal splitting of the recent trionychids). The results support the evolutionary stability of differentiated sex chromosomes in some lineages of ectothermic vertebrates. Moreover, our approach determining sex-linkage of protein coding genes can be used as a reliable technique of molecular sexing across trionychids useful for effective breeding strategy in conservation projects of endangered species.

The Historical Speciation of Mauremys Sensu Lato: Ancestral Area Reconstruction and Interspecific Gene Flow Level Assessment Provide New Insights

Mauremys sensu lato was divided into Mauremys, Chinemys, Ocadia, and Annamemys based on earlier research on morphology. Phylogenetic research on this group has been controversial because of disagreements regarding taxonomy, and the historical speciation is still poorly understood. In this study, 32 individuals of eight species that are widely distributed in Eurasia were collected. The complete mitochondrial (mt) sequences of 14 individuals of eight species were sequenced. Phylogenetic relationships, interspecific divergence times, and ancestral area reconstructions were explored using mt genome data (10,854 bp). Subsequent interspecific gene flow level assessment was performed using five unlinked polymorphic microsatellite loci. The Bayesian and maximum likelihood analyses revealed a paraphyletic relationship among four old genera (Mauremys, Annamemys, Chinemys, and Ocadia) and suggested the four old genera should be merged into the genus (Mauremys). Ancestral area reconstruction and divergence time estimation suggested Southeast Asia may be the area of origin for the common ancestral species of this genus and genetic drift may have played a decisive role in species divergence due to the isolated event of a glacial age. However, M. japonica may have been speciated due to the creation of the island of Japan. The detection of extensive gene flow suggested no vicariance occurred between Asia and Southeast Asia. Inconsistent results between gene flow assessment and phylogenetic analysis revealed the hybrid origin of M. mutica (Southeast Asian). Here ancestral area reconstruction and interspecific gene flow level assessment were first used to explore species origins and evolution of Mauremys sensu lato, which provided new insights on this genus.

Feeding behaviour in a 'basal' tortoise provides insights on the transitional feeding mode at the dawn of modern land turtle evolution

Almost all extant testudinids are highly associated with terrestrial habitats and the few tortoises with high affinity to aquatic environments are found within the genus Manouria. Manouria belongs to a clade which forms a sister taxon to all remaining tortoises and is suitable as a model for studying evolutionary transitions within modern turtles. We analysed the feeding behaviour of Manouria emys and due to its phylogenetic position, we hypothesise that the species might have retained some ancestral features associated with an aquatic lifestyle. We tested whether M. emys is able to feed both in aquatic and terrestrial environments. In fact, M. emys repetitively tried to reach submerged food items in water, but always failed to grasp them—no suction feeding mechanism was applied. When feeding on land, M. emys showed another peculiar behaviour; it grasped food items by its jaws—a behaviour typical for aquatic or semiaquatic turtles—and not by the tongue as generally accepted as the typical feeding mode in all tortoises studied so far. In M. emys, the hyolingual complex remained retracted during all food uptake sequences, but the food transport was entirely lingual based. The kinematical profiles significantly differed from those described for other tortoises and from those proposed from the general models on the function of the feeding systems in lower tetrapods. We conclude that the feeding behaviour of M. emys might reflect a remnant of the primordial condition expected in the aquatic ancestor of the tortoises.

Climate-mediated diversification of turtles in the Cretaceous

Chelonians are ectothermic, with an extensive fossil record preserved in diverse palaeoenvironmental settings: consequently, they represent excellent models for investigating organismal response to long-term environmental change. We present the first Mesozoic chelonian taxic richness curve, subsampled to remove geological/collection biases, and demonstrate that their palaeolatitudinal distributions were climate mediated. At the Jurassic/Cretaceous transition, marine taxa exhibit minimal diversity change, whereas non-marine diversity increases. A Late Cretaceous peak in 'global' non-marine subsampled richness coincides with high palaeolatitude occurrences and the Cretaceous thermal maximum (CTM): however, this peak also records increased geographic sampling and is not recovered in continental-scale diversity patterns. Nevertheless, a model-detrended richness series (insensitive to geographic sampling) also recovers a Late Cretaceous peak, suggesting genuine geographic range expansion among non-marine turtles during the CTM. Increased Late Cretaceous diversity derives from intensive North American sampling, but subsampling indicates that Early Cretaceous European/Asian diversity may have exceeded that of Late Cretaceous North America.

Population genomics of the endangered giant Galápagos tortoise

BACKGROUND

The giant Galápagos tortoise, Chelonoidis nigra, is a large-sized terrestrial chelonian of high patrimonial interest. The species recently colonized a small continental archipelago, the Galápagos Islands, where it has been facing novel environmental conditions and limited resource availability. To explore the genomic consequences of this ecological shift, we analyze the transcriptomic variability of five individuals of C. nigra, and compare it to similar data obtained from several continental species of turtles.

RESULTS

Having clarified the timing of divergence in the Chelonoidis genus, we report in C. nigra a very low level of genetic polymorphism, signatures of a weakened efficacy of purifying selection, and an elevated mutation load in coding and regulatory sequences. These results are consistent with the hypothesis of an extremely low long-term effective population size in this insular species. Functional evolutionary analyses reveal a reduced diversity of immunity genes in C. nigra, in line with the hypothesis of attenuated pathogen diversity in islands, and an increased selective pressure on genes involved in response to stress, potentially related to the climatic instability of its environment and its elongated lifespan. Finally, we detect no population structure or homozygosity excess in our five-individual sample.

CONCLUSIONS

These results enlighten the molecular evolution of an endangered taxon in a stressful environment and point to island endemic species as a promising model for the study of the deleterious effects on genome evolution of a reduced long-term population size.

The girdles of the oldest fossil turtle, Proterochersis robusta, and the age of the turtle crown

Background

Proterochersis robusta from the Late Triassic (Middle Norian) of Germany is the oldest known fossil turtle (i.e. amniote with a fully formed turtle shell), but little is known about its anatomy. A newly prepared, historic specimen provides novel insights into the morphology of the girdles and vertebral column of this taxon and the opportunity to reassess its phylogenetic position.

Results

The anatomy of the pectoral girdle of P. robusta is similar to that of other primitive turtles, including the Late Triassic (Carnian) Proganochelys quenstedti, in having a vertically oriented scapula, a large coracoid foramen, a short acromion process, and bony ridges that connect the acromion process with the dorsal process, glenoid, and coracoid, and by being able to rotate along a vertical axis. The pelvic elements are expanded distally and suturally attached to the shell, but in contrast to modern pleurodiran turtles the pelvis is associated with the sacral ribs.

Conclusions

The primary homology of the character “sutured pelvis” is unproblematic between P. robusta and extant pleurodires. However, integration of all new observations into the most complete phylogenetic analysis that support the pleurodiran nature of P. robusta reveals that this taxon is more parsimoniously placed along the phylogenetic stem of crown Testudines. All current phylogenetic hypotheses therefore support the basal placement of this taxon, imply that the sutured pelvis of this taxon developed independently from that of pleurodires, and conclude that the age of the turtle crown is Middle Jurassic.

A new xinjiangchelyid turtle from the Middle Jurassic of Xinjiang, China and the evolution of the basipterygoid process in Mesozoic turtles

BACKGROUND

Most turtles from the Middle and Late Jurassic of Asia are referred to the newly defined clade Xinjiangchelyidae, a group of mostly shell-based, generalized, small to mid-sized aquatic froms that are widely considered to represent the stem lineage of Cryptodira. Xinjiangchelyids provide us with great insights into the plesiomorphic anatomy of crown-cryptodires, the most diverse group of living turtles, and they are particularly relevant for understanding the origin and early divergence of the primary clades of extant turtles.

RESULTS

Exceptionally complete new xinjiangchelyid material from the ?Qigu Formation of the Turpan Basin (Xinjiang Autonomous Province, China) provides new insights into the anatomy of this group and is assigned to Xinjiangchelys wusu n. sp. A phylogenetic analysis places Xinjiangchelys wusu n. sp. in a monophyletic polytomy with other xinjiangchelyids, including Xinjiangchelys junggarensis, X. radiplicatoides, X. levensis and X. latiens. However, the analysis supports the unorthodox, though tentative placement of xinjiangchelyids and sinemydids outside of crown-group Testudines. A particularly interesting new observation is that the skull of this xinjiangchelyid retains such primitive features as a reduced interpterygoid vacuity and basipterygoid processes.

CONCLUSIONS

The homology of basipterygoid processes is confidently demonstrated based on a comprehensive review of the basicranial anatomy of Mesozoic turtles and a new nomenclatural system is introduced for the carotid canal system of turtles. The loss of the basipterygoid process and the bony enclosure of the carotid circulation system occurred a number of times independently during turtle evolution suggesting that the reinforcement of the basicranial region was essential for developing a rigid skull, thus paralleling the evolution of other amniote groups with massive skulls.

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.

The determinants of the molecular substitution process in turtles

Neutral rates of molecular evolution vary across species, and this variation has been shown to be related to biological traits. One of the first patterns to be observed in vertebrates has been an inverse relationship between body mass (BM) and substitution rates. The effects of three major life-history traits (LHT) that covary with BM - metabolic rate, generation time and longevity (LON) - have been invoked to explain this relationship. However, most of the theoretical and empirical evidence supporting this relationship comes from endothermic vertebrates, that is, mammals and birds, in which the environmental conditions, especially temperature, do not have a direct impact on cellular and molecular biology. We analysed the variations in mitochondrial and nuclear rates of synonymous substitution across 224 turtle species and examined their correlation with two LHT (LON and BM) and two environmental variables [latitude (LAT) and habitat]. Our analyses indicate that in turtles, neutral rates of molecular evolution are hardly correlated with LON or BM. Rather, both the mitochondrial and nuclear substitution rates are significantly correlated with LAT - faster evolution in the tropics - and especially so for aquatic species. These results question the generality of the relationships reported in mammals and birds and suggest that environmental factors might be the strongest determinants of the mutation rate in ectotherms.

Phylogenomic analyses support the position of turtles as the sister group of birds and crocodiles (Archosauria)

BACKGROUND

The morphological peculiarities of turtles have, for a long time, impeded their accurate placement in the phylogeny of amniotes. Molecular data used to address this major evolutionary question have so far been limited to a handful of markers and/or taxa. These studies have supported conflicting topologies, positioning turtles as either the sister group to all other reptiles, to lepidosaurs (tuatara, lizards and snakes), to archosaurs (birds and crocodiles), or to crocodilians. Genome-scale data have been shown to be useful in resolving other debated phylogenies, but no such adequate dataset is yet available for amniotes.

RESULTS

In this study, we used next-generation sequencing to obtain seven new transcriptomes from the blood, liver, or jaws of four turtles, a caiman, a lizard, and a lungfish. We used a phylogenomic dataset based on 248 nuclear genes (187,026 nucleotide sites) for 16 vertebrate taxa to resolve the origins of turtles. Maximum likelihood and Bayesian concatenation analyses and species tree approaches performed under the most realistic models of the nucleotide and amino acid substitution processes unambiguously support turtles as a sister group to birds and crocodiles. The use of more simplistic models of nucleotide substitution for both concatenation and species tree reconstruction methods leads to the artefactual grouping of turtles and crocodiles, most likely because of substitution saturation at third codon positions. Relaxed molecular clock methods estimate the divergence between turtles and archosaurs around 255 million years ago. The most recent common ancestor of living turtles, corresponding to the split between Pleurodira and Cryptodira, is estimated to have occurred around 157 million years ago, in the Upper Jurassic period. This is a more recent estimate than previously reported, and questions the interpretation of controversial Lower Jurassic fossils as being part of the extant turtles radiation.

CONCLUSIONS

These results provide a phylogenetic framework and timescale with which to interpret the evolution of the peculiar morphological, developmental, and molecular features of turtles within the amniotes.