Turtle origins: insights from phylogenetic retrofitting and molecular scaffolds

Adding new taxa to morphological phylogenetic analyses without substantially revising the set of included characters is a common practice, with drawbacks (undersampling of relevant characters) and potential benefits (character selection is not biased by preconceptions over the affinities of the ‘retrofitted’ taxon). Retrofitting turtles (Testudines) and other taxa to recent reptile phylogenies consistently places turtles with anapsid‐grade parareptiles (especially Eunotosaurus and/or pareiasauromorphs), under both Bayesian and parsimony analyses. This morphological evidence for turtle–parareptile affinities appears to contradict the robust genomic evidence that extant (living) turtles are nested within diapsids as sister to extant archosaurs (birds and crocodilians). However, the morphological data are almost equally consistent with a turtle–archosaur clade: enforcing this molecular scaffold onto the morphological data does not greatly increase tree length (parsimony) or reduce likelihood (Bayesian inference). Moreover, under certain analytic conditions, Eunotosaurus groups with turtles and thus also falls within the turtle–archosaur clade. This result raises the possibility that turtles could simultaneously be most closely related to a taxon traditionally considered a parareptile (Eunotosaurus) and still have archosaurs as their closest extant sister group.     

Reptile phylogeny and the interrelationships of turtles

A comprehensive analysis of amniote interrelationships is presented in an attempt to test turtle interrelationships. The results refute earlier hypotheses that turtles are related to parareptiles, i.e. to procolophonids or pareiasaurs. Instead, turtles are shown to be the sister-group of Sauropterygia, the two clades being nested within Sauria as sister-group of Lepidosauriformes. This scenario is also supported by several developmental and soft tissue characters which are shown to be congruent with the current phylogeny. The analysis strongly supports a monophyletic Parareptilia, sister-group of a monophylctic Eurcptilia. The Diapsida, however, is paraphyletic unless it includes turtles and sauropterygians. Additionally, the position of turtles within Diapsida has major implications for the evolutionary history and/or significance of many characters, i.e. temporal fenestration.     

Congruence,non‐homology,and the phylogeny of basal turtles

Joyce, W.G. and Sterli J. 2010. Congruence, non‐homology, and the phylogeny of basal turtles.–Acta Zoologica (Stockholm) Modern cladistic analysis is characterized by the assembly of increasingly larger data sets coupled with the use of congruence as the final test of homology. Some critics of this development have recently called for a return to more detailed primary homology analysis while questioning the utility of congruence. This discussion appears to be central to the debate regarding the phylogenetic relationships of basal turtles, as the large data sets developed by us have been criticized recently for utilizing poorly constructed characters and including too many homoplasy‐prone characters. Our analysis of this critique reveals that (1) new information regarding poorly understood taxa has a greater impact on the outcome of turtle phylogenies than the characters under dispute; (2) most current turtle phylogenies differ in taxon sampling, not character sampling, and so it appears illogical to condemn a particular analysis for its character sampling; (3) even evolutionary taxonomists should agree that key characters utilized to resolve basal turtle relationships cannot be thought to be ‘infallible’; (4) whereas various criteria provide positive evidence for homology, only congruence provides positive evidence for non‐homology; and (5) a stalemate between conflicting camps within a congruence frame work is preferable to the ad hoc dismissal of data sets, because authoritative statements are untestable.     

Molecular evidence for a clade of turtles.

Although turtles have been generally grouped with the most primitive reptile species, the origin and phylogenetic relationships of turtles have remained unresolved to date. To confirm the phylogenetic position of turtles in amniotes, we have cloned and determined the cDNA sequences encoding for skink lactate dehydrogenase (LDH)-A and LDH-B, snake LDH-A, and African clawed frog LDH-A; four alpha-enolase cDNA sequences from turtle, alligator, skink, and snake were also cloned and determined. All of these eight cDNA sequences, as well as the previously published LDH-A, LDH-B, and alpha-enolase of mammals, birds, reptiles, and African clawed frog, were analyzed by the phylogenetic tree reconstruction methods of neighbor-joining, maximum parsimony, and maximum likelihood. In the phylogenetic analyses, the turtle was found to be closely related to the alligator. Also, we found that the turtle had diverged after the divergence of squamates and birds. This departs from previous hypotheses of turtle evolution and further suggests that turtles are the latest of divergent reptiles, having been derived from an ancestor of crocodilian lineage within the last 200 million years.     

Turtle phylogeny: insights from a novel nuclear intron

Introns have gained considerable popularity as markers for molecular phylogenetics. However, no primers exist for a nuclear intron that amplifies across all turtles. Available data from morphology and mitochondrial DNA have not unambiguously resolved relationships within the superfamily Trionychoidea and the family Chelidae, which together form a large portion of extant turtle diversity. We tested the phylogenetic utility of a novel intron from the RNA fingerprint protein 35 (R35) as applied to these two areas of turtle systematics. We found the intron to be a single-copy locus that provides excellent resolving power for lineages among turtles, though problems with alignment made it impossible to infer deeper amniote relationships. Maximum parsimony and maximum likelihood both demonstrated the polyphyly of Trionychoidea and the reciprocal monophyly of Australian/New Guinea and South American chelid turtles. This is the first study to resolve such relationships with strong statistical support, and we suggest that R35 holds great promise for resolving additional persistent problems in the phylogeny of living turtles.     

Jaw musculature during the dawn of turtle evolution

Using a new approach to study muscle anatomy in vertebrates, the fully differentiated jaw musculature of 42 turtle species was studied and character mappings were performed. Soft tissue arrangements were correlated to the temporal openings (emarginations) of the skull and the trochlearis system of the jaw apparatus among turtle taxa. When compared to the cranial anatomy of stem Testudines, most characters detected as diagnostic of particular extant groups have to be considered as being evolved first within Testudines. Hence, jaw muscle anatomy of extant turtles is difficult to compare to that of other reptilian taxa. Moreover, the high number of apomorphic character changes speaks for a divergating turtle and saurian morphotype of jaw musculature, which could indicate either a position of turtles outside of Sauria or a highly derived, undetectable origin within that group. In general, a low direct correlation of most soft and hard tissue characters was detected. This finding could imply that both character complexes are more integrated to each other driven by functional morphology; i.e., the composition of muscle fibre types. That condition highlights the difficulty in using gross anatomy of jaw muscle characters to interpret temporal bone arrangements among amniotes in general.     

Derived appendicular morphology in Early Cretaceous aquatic turtles evidenced by the track record

A well‐preserved isolated manus impression evidences the presence of derived aquatic adaptations in the forelimb morphology of Berriasian turtles. Size and the abundant co‐occurring turtle fauna indicate that the track was left by a large‐sized, bottom‐walking basal pancryptodiran. The footprint shows an interesting mosaic of primitive and derived features. The basic turtle pattern is reflected by a short and broad autopodium but thin phalanges, claw reduction and extensive webbing enclosing the digits and lateral and medial margins of the distal metapodium are clearly specialized features related to the aquatic environment. Some characters including the proportional elongation of the first digit, which reached as far distally as the other digits and was probably highly mobile at the metacarpal‐phalangeal joint, indicate an even higher degree of specialization. The specimen gives evidence on a hitherto unknown range of appendicular modification in Early Cretaceous turtles, and redundantly proves the aquatic habitat and locomotion of large‐sized turtles from the Berriasian of northwestern Germany.     

Comparison of functional aspects of the coagulation cascade in human and sea turtle plasmas

Functional hemostatic pathways are critical for the survival of all vertebrates and have been evolving for more than 400 million years. The overwhelming majority of studies of hemostasis in vertebrates have focused on mammals with very sparse attention paid to reptiles. There have been virtually no studies of the coagulation pathway in sea turtles whose ancestors date back to the Jurassic period. Sea turtles are often exposed to rapidly altered environmental conditions during diving periods. This may reduce their blood pH during prolonged hypoxic dives. This report demonstrates that five species of turtles possess only one branch of the mammalian coagulation pathway, the extrinsic pathway. Mixing studies of turtle plasmas with human factor-deficient plasmas indicate that the intrinsic pathway factors VIII and IX are present in turtle plasma. These two factors may play a significant role in supporting the extrinsic pathway by feedback loops. The intrinsic factors, XI and XII are not detected which would account for the inability of reagents to induce coagulation via the intrinsic pathway in vitro. The analysis of two turtle factors, factor II (prothrombin) and factor X, demonstrates that they are antigenically/functionally similar to the corresponding human factors. The turtle coagulation pathway responds differentially to both pH and temperature relative to each turtle species and relative to human samples. The coagulation time (prothrombin time) increases as the temperature decreases between 37 and 15 °C. The increased time follows a linear relationship, with similar slopes for loggerhead, Kemps ridley and hawksbill turtles as well as for human samples. Leatherback turtle samples show a dramatic nonlinear increased time below 23 °C, and green turtle sample responses were similar but less dramatic. All samples also showed increased prothrombin times as the pH decreased from 7.8 to 6.4, except for three turtle species. The prothrombin times decreased, to varying extents, in a linear fashion relative to reduced pH with the rate of change greatest in leatherbacks>greenloggerhead turtles. All studies were conducted with reagents developed for human samples which would impact on the quantitative results with the turtle samples, but are not likely to alter the qualitative results. These comparative studies of the coagulation pathway in sea turtles and humans could enhance our knowledge of structure/function relationships and evolution of coagulation factors.     

Development and characterization of novel microsatellite markers from the olive ridley sea turtle (Lepidochelys olivacea)

Olive ridley turtles, although widely distributed globally and in Indian coastal waters, have undergone declines in recent years due to anthropogenic factors, particularly fishery‐related mortality. Assessment of genetic variability in existing populations is critical to the development of effective conservation strategies. Here we describe the development of six highly polymorphic microsatellite loci from a simple sequence repeat‐enriched genomic DNA library of olive ridley turtle. Characterization of five of these loci using 83 individual olive ridley turtles revealed eight to 24 alleles per locus, high observed and expected heterozygosity values and broad cross‐species amplifications. The sixth microsatellite was found to be monomorphic in the olive ridley samples but was polymorphic in two related marine turtle species. These microsatellites thus provide efficient genetic markers to understand the population structure, phylogeography and species relationships of olive ridley and other marine turtle species.     

A new, nearly complete stem turtle from the Jurassic of South America with implications for turtle evolution

Turtles have been known since the Upper Triassic (210Myr old); however, fossils recording the first steps of turtle evolution are scarce and often fragmentary. As a consequence, one of the main questions is whether living turtles (Testudines) originated during the Late Triassic (210Myr old) or during the Middle to Late Jurassic (ca 160Myr old). The discovery of the new fossil turtle, Condorchelys antiqua gen. et sp. nov. from the Middle to Upper Jurassic (ca 160-146Myr old) of South America (Patagonia, Argentina), presented here sheds new light on early turtle evolution. An updated cladistic analysis of turtles shows that C. antiqua and other fossil turtles are not crown turtles, but stem turtles. This cladistic analysis also shows that stem turtles were more diverse than previously thought, and that until the Middle to Upper Jurassic there were turtles without the modern jaw closure mechanism.     

A novel group of negative-sense RNA viruses associated with epizootics in managed and free-ranging freshwater turtles in Florida,USA

Few aquatic animal negative-sense RNA viruses have been characterized, and their role in disease is poorly understood. Here, we describe a virus isolated from diseased freshwater turtles from a Florida farm in 2007 and from an ongoing epizootic among free-ranging populations of Florida softshell turtles (Apalone ferox), Florida red-bellied cooters (Pseudemys nelsoni), and peninsula cooters (Pseudemys peninsularis). Affected turtles presented with similar neurological signs, oral and genital ulceration, and secondary microbial infections. Microscopic lesions were most severe in the softshell turtles and included heterophilic/histiocytic meningoencephalitis, multi-organ vasculitis, and cytologic observation of leukocytic intracytoplasmic inclusions. The virus was isolated using Terrapene heart (TH-1) cells. Ultrastructurally, viral particles were round to pleomorphic and acquired an envelope with prominent surface projections by budding from the cell membrane. Viral genomes were sequenced from cDNA libraries of two nearly identical isolates and determined to be bi-segmented, with an ambisense coding arrangement. The larger segment encodes a predicted RNA-directed RNA polymerase (RdRP) and a putative zinc-binding matrix protein. The smaller segment encodes a putative nucleoprotein and an envelope glycoprotein precursor (GPC). Thus, the genome organization of this turtle virus resembles that of arenaviruses. Phylogenetic analysis shows that the RdRP of the turtle virus is highly diverged from the RdRPs of all known negative-sense RNA viruses and forms a deep branch within the phylum Negarnaviricota, that is not affiliated with any known group of viruses, even at the class level. In contrast, the GPC protein of the turtle virus is confidently affiliated with homologs from a distinct group of fish hantaviruses. Thus, the turtle virus is expected to become the founder of a new taxon of negative-sense RNA viruses, at least with a family rank, but likely, an order or even a class. These viruses probably evolved either by reassortment or by intrasegment recombination between a virus from a distinct branch of negarnaviruses distant from all known groups and a hanta-like aquatic virus. We suggest the provisional name Tosoviridae for the putative new family, with Turtle fraservirus 1 (TFV1) as the type species within the genus Fraservirus. A conventional RT-PCR assay, targeting the TFV1 RdRP, confirmed the presence of viral RNA in multiple tissues and exudates from diseased turtles. The systemic nature of the TFV1 infection was further supported by labeling of cells within lesions using in situ hybridization targeting the RNA of the TFV1 RdRP.     

Diversity and structure of turtle assemblages: associations with wetland characters across a floodplain landscape

We studied 24 wetlands along a 296‐km portion of the lower Missouri River floodplain, USA during 1996–1998. Our goal was to examine relationships between wetland characters and freshwater turtle diversity. We collected data on a total of 2201 individual captures of six freshwater turtle species. Ten species variables (richness, evenness, proportion of lentic species combined, proportion of lotic species combined, and proportions of six individual species) described the turtle assemblage structure and were associated with six broad wetland categories (remnant oxbow, remnant slough, connected scour, isolated scour, temporary farmed, and temporary not farmed) as well as with six abiotic (water depth, distance from the river, duration dry, duration flooded, size, and turbidity) and two biotic (primary and secondary productivity) wetland characters. Repeated‐measures MANOVA and ANOVA revealed highest diversity in remnant sloughs and oxbows as well as in newly created scour wetlands and lowest diversity in temporary wetlands. Multiple regressions and correlations indicated that the single most important wetland characteristic for high turtle diversity was a low annual duration of drying. Oxbows supported the highest proportions of lentic turtles; while connected scours supported the highest proportions of lentic turtles. Lentic species proportions were highest in wetlands that were frequently flooded, relatively distant from the Missouri River, and insect rich, while lotic species proportions were highest in wetlands that were turbid, close to the river, and relatively insect poor. The cohabitation patterns of several turtle species suggest that direct or indirect interactions between species may affect assemblage structure. Based on our broad array of wetlands, we identified crucial wetland types and characteristics that were associated with maximum diversity for freshwater turtles. Empirical studies similar to ours across broad landscapes are needed for a variety of wetland floral and faunal groups before sound conservation and management recommendations can be formulated.     

Using Genes as Characters and a Parsimony Analysis to Explore the Phylogenetic Position of Turtles

The phylogenetic position of turtles within the vertebrate tree of life remains controversial. Conflicting conclusions from different studies are likely a consequence of systematic error in the tree construction process, rather than random error from small amounts of data. Using genomic data, we evaluate the phylogenetic position of turtles with both conventional concatenated data analysis and a “genes as characters” approach. Two datasets were constructed, one with seven species (human, opossum, zebra finch, chicken, green anole, Chinese pond turtle, and western clawed frog) and 4584 orthologous genes, and the second with four additional species (soft-shelled turtle, Nile crocodile, royal python, and tuatara) but only 1638 genes. Our concatenated data analysis strongly supported turtle as the sister-group to archosaurs (the archosaur hypothesis), similar to several recent genomic data based studies using similar methods. When using genes as characters and gene trees as character-state trees with equal weighting for each gene, however, our parsimony analysis suggested that turtles are possibly sister-group to diapsids, archosaurs, or lepidosaurs. None of these resolutions were strongly supported by bootstraps. Furthermore, our incongruence analysis clearly demonstrated that there is a large amount of inconsistency among genes and most of the conflict relates to the placement of turtles. We conclude that the uncertain placement of turtles is a reflection of the true state of nature. Concatenated data analysis of large and heterogeneous datasets likely suffers from systematic error and over-estimates of confidence as a consequence of a large number of characters. Using genes as characters offers an alternative for phylogenomic analysis. It has potential to reduce systematic error, such as data heterogeneity and long-branch attraction, and it can also avoid problems associated with computation time and model selection. Finally, treating genes as characters provides a convenient method for examining gene and genome evolution.     

Integration of morphological data sets for phylogenetic analysis of Amniota: the importance of integumentary characters and increased taxonomic sampling

Several mutually exclusive hypotheses have been advanced to explain the phylogenetic position of turtles among amniotes. Traditional morphology-based analyses place turtles among extinct anapsids (reptiles with a solid skull roof), whereas more recent studies of both morphological and molecular data support an origin of turtles from within Diapsida (reptiles with a doubly fenestrated skull roof). Evaluation of these conflicting hypotheses has been hampered by nonoverlapping taxonomic samples and the exclusion of significant taxa from published analyses. Furthermore, although data from soft tissues and anatomical systems such as the integument may be particularly relevant to this problem, they are often excluded from large-scale analyses of morphological systematics. Here, conflicting hypotheses of turtle relationships are tested by (1) combining published data into a supermatrix of morphological characters to address issues of character conflict and missing data; (2) increasing taxonomic sampling by more than doubling the number of operational taxonomic units to test internal relationships within suprageneric ingroup taxa; and (3) increasing character sampling by approximately 25% by adding new data on the osteology and histology of the integument, an anatomical system that has been historically underrepresented in morphological systematics. The morphological data set assembled here represents the largest yet compiled for Amniota. Reevaluation of character data from prior studies of amniote phylogeny favors the hypothesis that turtles indeed have diapsid affinities. Addition of new ingroup taxa alone leads to a decrease in overall phylogenetic resolution, indicating that existing characters used for amniote phylogeny are insufficient to explain the evolution of more highly nested taxa. Incorporation of new data from the soft and osseous components of the integument, however, helps resolve relationships among both basal and highly nested amniote taxa. Analysis of a data set compiled from published sources and data original to this study supports monophyly of Amniota, Synapsida, Reptilia, Parareptilia, Eureptilia, Eosuchia, Diapsida, Neodiapsida, Sauria, Lepidosauria, and Archosauriformes, as well as several more highly nested divisions within the latter two clades. Turtles are here resolved as the sister taxon to a monophyletic Lepidosauria (squamates + Sphenodon), a novel phylogenetic position that nevertheless is consistent with recent molecular and morphological studies that have hypothesized diapsid affinities for this clade.     

Pareiasaur phylogeny and the origin of turtles

The evolutionary relationship of all the valid species (and thus genera) of pareiasaurs are assessed through a phylogcnctic analysis of these taxa together with turtles, Owenetta, Barasaurus, Sclerosaurus, procolophonids, lanthanosuchids, nyctiphruretids, and nycterolctcrids. 128 os-teological characters were used, and almost all relevant taxa were examined. The results confirm that among these taxa, pareiasaurs and turtles form a robust clade, to the exclusion of all other taxa including procolophonids. However, pareiasaurs might not be the mono-phyletic sister group of turtles, as previously suggested. Rather, there is some evidence that pareiasaurs are paraphyletic with respect to (i.e. ‘ancestral to’) turtles. Among pareiasaurs, the early, large, heavily ossified forms such as Brady.saurus are most distantly related to turtles. These forms are characterized by rather smooth skulls, and dermal armour restricted to the dorsal midline. More closely related to turtles are forms such as Scutosaurus, Pareiasuchus, and Elginia. These taxa form a distinct clade of pareiasaurs, characterized by a very ‘mammallike’ pelvis, elaborate cranial ornamentation and a loose covering of osteoderms over the entire dorsum. The late, dwarf pareiasaurs Nanoparia, Anthodon, and Pumiliopareia are the nearest relatives of turtles. These forms exhibit otherwise uniquely turtle features such as a rigid covering of dermal armour over the entire dorsal region, expanded flattened ribs, cylindrical scapula blade, great reduction of humeral torsion (to 25^o), greatly developed trochanter major, offset femoral head, and reduced cnemial crest of the tibia. Thus, many features thought to be restricted to turtles (and thus to have evolved simultaneously with the turtle shell) actually arose earlier, at various points along the pareiasaurian stem lineage. The identification of the nature and sequence of anatomical changes leading to the origin of turtles, and the possibility that turtles are derived from dwarf pareiasaurs, should have important implications for speculations on the evolutionary biology of turtle origins.     

Triassic turtle tracks and the origin of turtles

Turtle (Testudines) tracks, Chelonipus torquatus, reported from the early Middle Triassic (Anisian) of Germany, and Chelonipus isp. from the late Early Triassic (Spathian) of Wyoming and Utah, are the oldest fossil evidence of turtles, but have been omitted in recent discussions of turtle origins. These tracks provide significant clues as to how early the turtle Bauplan originated. Turtle trackways are quite distinctive: the manus and pes form tracks nearly parallel to the midline and indicate an unusually wide gait in which the trackway width is nearly equal to the stride length. These tracks do not fit what would be expected to be made by Triassic Pappochelys or Odontochelys, a supposed prototurtle and an early turtle, respectively. In contrast, these tracks are consistent with what would be expected from the Triassic turtles Proganochelys and Palaeochersis. The features inferred to be present in Triassic turtle tracks support the notion that Odontochelys is a derived aquatic branch of the turtle stem lineage rather than the ancestral state of all turtles. Chelonipus also resembles the Permian track Pachypes dolomiticus, generally assigned to a pareiasaur trackmaker. These revelations highlight the need to consider all available evidence regarding turtle origins, rather than just the body fossils.     

Integrating fossil preservation biases in the selection of calibrations for molecular divergence time estimation

The selection of fossil data to use as calibration age priors in molecular divergence time estimates inherently links neontological methods with paleontological theory. However, few neontological studies have taken into account the possibility of a taphonomic bias in the fossil record when developing approaches to fossil calibration selection. The Sppil-Rongis effect may bias the first appearance of a lineage toward the recent causing most objective calibration selection approaches to erroneously exclude appropriate calibrations or to incorporate multiple calibrations that are too young to accurately represent the divergence times of target lineages. Using turtles as a case study, we develop a Bayesian extension to the fossil selection approach developed by Marshall (2008. A simple method for bracketing absolute divergence times on molecular phylogenies using multiple fossil calibrations points. Am. Nat. 171:726-742) that takes into account this taphonomic bias. Our method has the advantage of identifying calibrations that may bias age estimates to be too recent while incorporating uncertainty in phylogenetic parameter estimates such as tree topology and branch lengths. Additionally, this method is easily adapted to assess the consistency of potential calibrations to any one calibration in the candidate pool.     

Relationship between tropical cyclones and the distribution of sea turtle nesting grounds

Aim This study examines the relationship between the distribution of existing sea turtle nesting sites and historical patterns of tropical cyclone events to investigate whether cyclones influence the current distribution of sea turtle nesting sites. The results, together with information on predicted cyclone activity and other key environmental variables, will help in the identification and prediction of future nesting sites for sea turtles as changes to the coastal environment continue. Location Queensland, Australia. Methods We used data on the nesting distribution of seven populations of four species of sea turtles [green (Chelonia mydas), flatback (Natator depressus), hawksbill (Eretmochelys imbricata) and loggerhead (Caretta caretta)] from the eastern Queensland coast, and tropical cyclone track data from 1969 to 2007 to explore the relationship between (1) sea turtle nesting phenology and cyclone season, and (2) sea turtle nesting sites and cyclone distribution. Furthermore, using two green turtle populations as a case study, we investigated the relationship between cyclone disturbance and sea turtle reproductive output, nesting site and season. Bootstrapping was used to explore if current sea turtle nesting sites are located in areas with lower or higher cyclone frequency than areas where turtles are currently not nesting. Results All populations of sea turtles studied here were disturbed by cyclone activity during the study period. The exposure (frequency) of tropical cyclones that crossed each nesting site varied greatly among and within the various sea turtle populations. This was mainly a result of the spatial distribution of each population’s nesting sites. Bootstrapping indicated that nesting sites generally have experienced lower cyclone activity than other areas that are available for nesting. Main conclusions Tropical cyclones might have been sufficiently detrimental to sea turtle hatching success on the eastern Queensland coast that through a natural selection process turtles in this region are now nesting in areas with lower cyclone activity. Therefore, it is important that future studies that predict climate or range shifts for sea turtle nesting distributions consider future cyclone activity as one of the variables in their model.