A global phylogeny of turtles reveals a burst of climate-associated diversification on continental margins

A longitudinal analysis of function annotations of the human proteome reveals consistently high biases

The resources required to study gene function are limited, especially when considering the number of genes in the human genome and the complexity of their function. Therefore, genes are prioritized for experimental studies based on many different considerations, including, but not limited to, perceived biomedical importance, such as disease-associated genes, or the understanding of biological processes, such as cell signalling pathways. At the same time, most genes are not studied or are under-characterized, which hampers our understanding of their function and potential effects on human health and wellness. Understanding function annotation disparity is a necessary first step toward understanding how much functional knowledge is gained from the human genome, and toward guidelines for better targeting future studies of the genes in the human genome effectively. Here, we present a comprehensive longitudinal analysis of the human proteome utilizing data analysis tools from economics and information theory. Specifically, we view the human proteome as a population of proteins within a knowledge economy: we treat the quantified knowledge of the protein's function as the analogue of wealth and examine the distribution of information in a population of proteins in the proteome in the same manner distribution of wealth is studied in societies. Our results show a highly skewed distribution of information about human proteins over the last decade, in which the inequality in the annotations given to the proteins remains high. Additionally, we examine the correlation between the knowledge about protein function as captured in databases and the interest in proteins as reflected by mentions in the scientific literature. We show a large gap between knowledge and interest and dissect the factors leading to this gap. In conclusion, our study shows that research efforts should be redirected to less studied proteins to mitigate the disparity among human proteins both in databases and literature.

Using the axial skeleton as armor: mechanical behavior of sea turtle carapaces throughout ontogeny

The shells of turtles serve as protection, yet shell shape and natural history widely vary among turtles. Here, we identify the mechanical behavior that provides marine turtles, species characterized with fusiform shells, with biomechanical strength and resilience. The multi-layered carapacial bone structure seemingly serves a protective role for the muscles, nerves and viscera it houses. What are the shell's material properties that provide protection? Most previous work has focused on non-marine turtles, which differ in natural history and shell morphology from marine species. We measured carapacial mechanical behavior of green turtle (Chelonia mydas), loggerhead (Caretta caretta) and Kemp's ridleys (Lepidochelys kempii) across a range of body sizes in juveniles, subadults and adults. Carapace samples were tested using quasi-static compression to quantify stiffness (Young's modulus), yield strength and toughness. The mechanical characteristics of marine turtle shells are grossly akin to those of other turtles and driven by the bone's sandwich structure. Yet, the material properties indicate that marine turtle shells are less stiff and strong than those of their freshwater and terrestrial counterparts. We hypothesize that increased flexibility of the shell may reflect tradeoffs for life that include experiencing pressure from diving somewhat deeply in marine environments. Shell material properties also differ among species and ontogenetically. Green turtles have the stiffest, strongest and toughest shells while loggerhead carapaces are the most compliant. Stiffness and yield strength show positive relationships with body size which are most pronounced in green turtles and Kemp's ridleys. Phylogenetic histories and ecological differences likely drive this interspecific variation.

Testing Relationships between Multiple Regional Features and Biogeographic Processes of Speciation, Extinction, and Dispersal

The spatial and environmental features of regions where clades are evolving are expected to impact biogeographic processes such as speciation, extinction, and dispersal. Any number of regional features (such as elevation, distance, area, etc.) may be directly or indirectly related to these processes. For example, it may be that distances or differences in elevation or both may limit dispersal rates. However, it is difficult to disentangle which features are most strongly related to rates of different processes. Here, we present an extensible Multifeature Feature-Informed GeoSSE (MultiFIG) model that allows for the simultaneous investigation of any number of regional features. MultiFIG provides a conceptual framework for incorporating large numbers of features of different types, including categorical, quantitative, within-region, and between-region features, along with a mathematical framework for translating those features into biogeographic rates for statistical hypothesis testing. Using traditional Bayesian parameter estimation and reversible-jump Markov chain Monte Carlo, MultiFIG allows for the exploration of models with different numbers and combinations of feature-effect parameters and generates estimates for the strengths of relationships between each regional feature and core process. We validate this model with a simulation study covering a range of scenarios with different numbers of regions, tree sizes, and feature values. We also demonstrate the application of MultiFIG with an empirical case study of the South American lizard genus Liolaemus, investigating 16 regional features related to area, distance, and elevation. Our results show two important feature-process relationships: a negative distance/dispersal relationship and a negative area/extinction relationship. Interestingly, although speciation rates were found to be higher in Andean versus non-Andean regions, the model did not assign significance to Andean- or elevation-related parameters. These results highlight the need to consider multiple regional features in biogeographic hypothesis testing.

A Unified Method for Detecting Phylogenetic Signals in Continuous, Discrete, and Multiple Trait Combinations

Phylogenetic signals are widely used in ecological and evolutionary studies. Trait data used to detect phylogenetic signals can be continuous or discrete, but existing indices are designed for either type, not both. Moreover, most existing methods can only perform phylogenetic detection of individual traits, despite the fact that biological functions are often the result of interactions among multiple traits. Some attempts to detect phylogenetic signals across multiple trait combinations have employed alternative indicators, which may not align perfectly with the rigorous criteria for defining phylogenetic signals. In this study, we developed a new index (the M statistic) to detect phylogenetic signals for continuous traits, discrete traits, and multiple trait combinations. This capability is inherited from Gower's distance, which is used in the calculation of the M statistic to convert various types of traits into distances. The M statistic strictly adheres to the definition of phylogenetic signals and detects them by comparing these distances from phylogenies and traits. Using simulated data, we compared the performance of our new approach with that of existing commonly used indices. The results show that our method is not inferior to the existing methods. It performs well in handling continuous variables, discrete variables, and multiple trait combinations. We used trait data of turtles (Testudines) to demonstrate the utility of our new method. We suggest this new index as an original method for the detection of phylogenetic signals across various variable types. We provide an R package called "phylosignalDB" to facilitate all calculations.

Amphibians and Reptiles Exhibit Different Ecological and Evolutionary Spatial Patterns in the Amazon Basin

Understanding spatial variability in ecological and evolutionary patterns is key to Amazonian biodiversity conservation. This study examined taxonomic, phylogenetic, and functional diversity across amphibians and reptiles, assessing the influence of elevation, interrelationships among metrics, and distribution across five Amazon Basin ecoregions, exploring the "cradle" (speciation) and "museum" (lineage preservation) hypotheses. We analyzed 1011 amphibian species from three lineages and 828 reptile species from four lineages. Integrating distribution maps, phylogenies, and trait data, we calculated phylogenetic (PD), functional (FD), and taxonomic (TD) diversity, including mean phylogenetic (PDmntd) and functional (FDmntd) distance to the nearest taxon. We examined spatial regressions between diversity metrics and elevation, assessed correlations among metrics, and compared diversity metrics across ecoregions for each lineage. Diversity metrics across amphibian and reptile lineages reveal distinct geographical patterns related to elevation. Anurans exhibit higher PD, FD, and TD in the western Amazon, while squamates show hotspots at low altitudes. Testudines are linked to major rivers, and crocodilians display high PD near the equator. Anurans and squamates show elevated PDmntd and FDmntd in the Andes, whereas testudines are found in cratonic regions. Significant correlations and notable differences among ecoregions were found, especially in the Andes and low regions of the Amazon Basin. This study highlights the diverse eco-evolutionary patterns among amphibian and reptile lineages in the Amazon Basin, each exhibiting distinct hotspots distributed across ecoregions. The findings align with the cradle-museum hypothesis, suggesting that some regions serve as centers of ongoing diversification, others preserve ancient lineages, or serve as both. The cradle-museum hypothesis should be carefully analyzed, as each taxon presents a distinct pattern. This research underscores the necessity for targeted conservation strategies tailored to distinct ecological and evolutionary dynamics across ecoregions.

Comparative analysis of the mitochondrial genomes of the soft-shelled turtles Palea steindachneri and Pelodiscus axenaria and phylogenetic implications for Trionychia

Soft-shelled turtles, or Trionychia, are an enigmatic and fascinating group due to their specific morphological features and ecological adaptations. Based on mitochondrial DNA (mtDNA) and/or nuclear markers, previous studies showed the incongruent phylogenetic topologies within Trionychia (e.g., the Palea and its closely related species). In order to resolve the equivocal relationships and obtain some "genome-level" common evolutionary characters of soft-shelled turtles, in this study, we assembled and annotated the complete mitochondrial genomes of Palea steindachneri and Pelodiscus axenaria, both naturally distributed in Asia. The sizes of the two mitochondrial genomes were 16,811 bp and 17,143 bp, respectively. Typical vertebrate animal mtDNA features were observed, such as the usual gene components and arrangements (37 genes with a non-coding control region) and the A + T biased nucleotide compositions on the light strand (61.5% and 62.7%, respectively). All conserved blocks common to the vertebrates control region except for the extended terminal associated sequences (ETAS2) were found in the two soft-shelled turtles. The ω ratio averaged over all sites of each protein-coding gene (PCG) was below 1, which indicated purifying selection at the gene-wide level. However, a positive selection site at the 350-codon position in the cytb gene was detected, as estimated by Bayes empirical Bayes (BEB) analysis. Compared with the gene subsets, the mitogenomes provided the most robust phylogenetic resolution. The monophyly of the clades Amydona, Gigantaesuarochelys, and Apalonia was well supported. Topology discrepancies were observed among different datasets (e.g., the positions of Lissemys and Palea), reflecting the heterogeneous phylogenetic signals in the soft-shelled turtle mitogenomes. Precise date estimation based on Bayesian relaxed clock analyses indicated that the crown group age of extant Trionychia was approximately 115.84 Ma (95% HPD: 91.33-142.18 Ma). Paleoclimate changes, especially the Eocene - Oligocene transition, could be responsible for the speciation in these groups. Our results reiterated the necessity and effectiveness of incorporating entire mitochondrial genomes to delineate phylogenetic relationships in chelonian phylogeny studies.

Genomic insights into marine environment adaptation and conservation of the threatened olive ridley turtle (Lepidochelys olivacea)

Sea turtles are marine flagship species and most of them are currently in a threatened state. Long-term surviving in the ocean has driven significant morphological and physiological changes for this group, which makes them an ideal model for studying adaptive evolution of marine environments. Herein, we present a chromosome-scale genome of Lepidochelys olivacea with a genome size of 2.22 Gb and a contig N50 of 97.3 Mb. Comparative genomic analyses uncovered a suite of adaptive changes in genes related to olfaction, vision, virus defense, and longevity, which may help explain the genetic underpinnings of its marine environment adaptation. We also observed that the genome-wide heterozygosity of L. olivacea was low (6.45e-4), consistent with its prolonged population decline. Overall, our study provides valuable genetic resources for understanding evolutionary adaptations to aquatic environment and for the conservation of this threatened species.

Complex Models of Sequence Evolution Improve Fit, But Not Gene Tree Discordance, for Tetrapod Mitogenomes

Variation in gene tree estimates is widely observed in empirical phylogenomic data and is often assumed to be the result of biological processes. However, a recent study using tetrapod mitochondrial genomes to control for biological sources of variation due to their haploid, uniparentally inherited, and non-recombining nature found that levels of discordance among mitochondrial gene trees were comparable to those found in studies that assume only biological sources of variation. Additionally, they found that several of the models of sequence evolution chosen to infer gene trees were doing an inadequate job of fitting the sequence data. These results indicated that significant amounts of gene tree discordance in empirical data may be due to poor fit of sequence evolution models and that more complex and biologically realistic models may be needed. To test how the fit of sequence evolution models relates to gene tree discordance, we analyzed the same mitochondrial data sets as the previous study using 2 additional, more complex models of sequence evolution that each include a different biologically realistic aspect of the evolutionary process: A covarion model to incorporate site-specific rate variation across lineages (heterotachy), and a partitioned model to incorporate variable evolutionary patterns by codon position. Our results show that both additional models fit the data better than the models used in the previous study, with the covarion being consistently and strongly preferred as tree size increases. However, even these more preferred models still inferred highly discordant mitochondrial gene trees, thus deepening the mystery around what we label the "Mito-Phylo Paradox" and leading us to ask whether the observed variation could, in fact, be biological in nature after all.

Conservation of OFD1 Protein Motifs: Implications for Discovery of Novel Interactors and the OFD1 Function

OFD1 is a protein involved in many cellular processes, including cilia biogenesis, mitotic spindle assembly, translation, autophagy and the repair of double-strand DNA breaks. Despite many potential interactors identified in high-throughput studies, only a few have been directly confirmed with their binding sites identified. We performed an analysis of the evolutionary conservation of the OFD1 sequence in three clades: 80 Tetrapoda, 144 Vertebrata or 26 Animalia species, and identified 59 protein-binding motifs localized in the OFD1 regions conserved in various clades. Our results indicate that OFD1 contains 14 potential post-translational modification (PTM) sites targeted by at least eight protein kinases, seven motifs bound by proteins recognizing phosphorylated aa residues and a binding site for phosphatase 2A. Moreover, OFD1 harbors both a motif that enables its phosphorylation by mitogen-activated protein kinases (MAPKs) and a specific docking site for these proteins. Generally, our results suggest that OFD1 forms a scaffold for interaction with many proteins and is tightly regulated by PTMs and ligands. Future research on OFD1 should focus on the regulation of OFD1 function and localization.

A review of age estimation methods in non-avian reptiles by growth marks in hard tissues

Age and growth-related data are basic biological parameters, essential in population ecology, evolution, and conservation biology. There is a growing body of published information on reptile demography derived from sclerochronology, a technique based on counting the growth layers deposited in bones (skeletochronology) and other hard body structures. Since the data are not always easily available, we compiled the existing published data, described the current status of knowledge, synthetized the conclusions of disparate studies, and identified patterns of research and information gaps, prioritizing the needs for future research. Our database includes the results of 468 published studies covering 236 reptile species from 41 families. These represent less than 2% of the total number of known extant species. Turtles and crocodiles are proportionally better studied, while snakes are the least examined group. The distribution of the research does not reflect conservation needs; we found an important geographic bias, with an overrepresentation of Northern temperate species. Only 23% of the studies checked the assumption of periodicity of growth marks deposition, and the method was found to be reliable or adequate in 79% of the cases. Overall, the data obtained through sclerochronology can be considered robust, especially if validation methods are employed, since the general goal is to characterize population parameters, trends, and dynamics, rather than determining the exact age of any specimen in particular.

Discovery of the first sea turtle adenovirus and turtle associated circoviruses

Turtles are an evolutionarily unique and morphologically distinctive order of reptiles, and many species are globally endangered. Although a high diversity of adenoviruses in scaled reptiles is well-documented, turtle adenoviruses remain largely understudied. To investigate their molecular diversity, we focused on the identification and characterisation of adenoviruses in turtle-derived organ, swab and egg samples. Since reptile circoviruses have been scarcely reported and no turtle circoviruses have been documented to date, we also screened our samples for circoviruses. Host-virus coevolution is a common feature of these viral families, so we aimed to investigate possible signs of this as well. Two screening projects were conducted: one on Brazilian samples collected from animals in their natural habitat, and the other on Hungarian pet shop samples. Nested PCR systems were used for the detection of adeno- and circoviruses and purified PCR products were Sanger sequenced. Phylogenetic trees for the viruses were reconstructed based on the adenoviral DNA polymerase and hexon genes, circoviral Rep genes, and for the turtle hosts based on mitochondrial cytochrome b amino acid sequences. During the screening, testadeno-, siadeno-, and circovirus strains were detected. The circovirus strains were classified into the genus Circovirus, exhibiting significant evolutionary divergence but forming a monophyletic clade within a group of fish circoviruses. The phylogenetic tree of turtles reflected their taxonomic relationships, showing a deep bifurcation between suborders and distinct monophyletic clades corresponding to families. A similar clustering pattern was observed among the testadenovirus strains in their phylogenetic tree. As a result, this screening of turtle samples revealed at least three new testadenoviruses, including the first sea turtle adenovirus, evidence of coevolution between testadenoviruses and their hosts, and the first turtle associated circoviruses. These findings underscore the need for further research on viruses in turtles, and more broadly in reptiles, to better understand their viral diversity and the evolutionary processes shaping host-virus interactions.

Functional and Character Disparity Are Decoupled in Turtle Mandibles

Turtles have high shape variation of their mandibles, likely reflecting adaptations to a broad variety of food items and ingestion strategies. Here, we compare functional disparity measured by biomechanical proxies and character disparity measured by discrete morphological characters. Functional and character disparities vary between clades and ecological groups and are thus decoupled. Comparisons with cranial disparity also indicate decoupled patterns within the turtle skull. Exploration of mandibular patterns reveals that several biomechanical configurations or character state combinations can lead to the same feeding type (i.e., convergence) or that high functional disparity can be achieved at a low exhaustion of character state combinations (e.g., cryptodires). Dietary specialists show larger functional disparity than generalists, but the phylogenetically widespread generalist ecology leads to high character disparity signals in the ecotype. Whereas character disparity generally shows high phylogenetic signal, functional disparity patterns correspond to dietary specializations, which may occur convergently across different groups. Despite this, individual functional measurements have overlapping ranges across ecogroups and do not always conform to biomechanical expectations. Jaw opening and closing biomechanical advantages model trade-offs between force transmission and opening/closing speeds, and turtles show a variety of combinations of values that we try to synthesize into several "jaw types". Closing mechanical advantage shows that turtles retain high levels of force transmission at the anterior jaw end compared with other groups (e.g., pseudosuchians). This can possibly be explained as an evolutionary adaptation to retain high bite forces at small head sizes.

A systematic comparative description of extant turtle humeri, with comments on humerus disparity and evolution based on fossil comparisons

The humerus is central for locomotion in turtles as quadrupedal animals. Osteological variation across testudine clades remains poorly documented. Here, we systematically describe the humerus anatomy for all major extant turtle clades based on 38 species representing the phylogenetic and ecological diversity of crown turtles. Three Late Triassic species of shelled stem turtles (Testudindata) are included to establish the plesiomorphic humerus morphology. Our work is based on 3D models, establishing a publicly available digital database. Previously defined terms for anatomical sides of the humerus (e.g., dorsal, ventral) are often not aligned with the respective body sides in turtles and other quadrupedal animals with sprawling gait. We propose alternative anatomical directional terms to simplify communication: radial and ulnar (the sides articulating with the radius/ulna), capitular (the side bearing the humeral head), and intertubercular (opposite to capitular surface). Turtle humeri show low morphological variation with exceptions concentrated in locomotory specialists. We propose 15 discrete characters to summarize osteological variation for future phylogenetic studies. Disparity analyses comparing non-shelled and shelled turtles indicate that the presence of the shell constrains humerus variation. Flippered aquatic turtles are released from this constraint and significantly increase overall disparity. Ontogenetic changes of turtle humeri are related to increased ossification and pronunciation of the proximal processes, the distal articulation areas, and the closure of the ectepicondylar groove to a foramen. Some turtle species retain juvenile features into adulthood and provide evidence for paedomorphic evolution. We review major changes of turtle humerus morphology throughout the evolution of its stem group.

De novo genome assemblies of two cryptodiran turtles with ZZ/ZW and XX/XY sex chromosomes provide insights into patterns of genome reshuffling and uncover novel 3D genome folding in amniotes

Understanding the evolution of chromatin conformation among species is fundamental to elucidate the architecture and plasticity of genomes. Nonrandom interactions of linearly distant loci regulate gene function in species-specific patterns, affecting genome function, evolution, and, ultimately, speciation. Yet, data from nonmodel organisms are scarce. To capture the macroevolutionary diversity of vertebrate chromatin conformation, here we generate de novo genome assemblies for two cryptodiran (hidden-neck) turtles via Illumina sequencing, chromosome conformation capture, and RNA-seq: Apalone spinifera (ZZ/ZW, 2n = 66) and Staurotypus triporcatus (XX/XY, 2n = 54). We detected differences in the three-dimensional (3D) chromatin structure in turtles compared to other amniotes beyond the fusion/fission events detected in the linear genomes. Namely, whole-genome comparisons revealed distinct trends of chromosome rearrangements in turtles: (1) a low rate of genome reshuffling in Apalone (Trionychidae) whose karyotype is highly conserved when compared to chicken (likely ancestral for turtles), and (2) a moderate rate of fusions/fissions in Staurotypus (Kinosternidae) and Trachemys scripta (Emydidae). Furthermore, we identified a chromosome folding pattern that enables "centromere-telomere interactions" previously undetected in turtles. The combined turtle pattern of "centromere-telomere interactions" (discovered here) plus "centromere clustering" (previously reported in sauropsids) is novel for amniotes and it counters previous hypotheses about amniote 3D chromatin structure. We hypothesize that the divergent pattern found in turtles originated from an amniote ancestral state defined by a nuclear configuration with extensive associations among microchromosomes that were preserved upon the reshuffling of the linear genome.

Chromosome-level genome of diamondback terrapin provides insight into the genetic basis of salinity adaptation

Diamondback terrapins (Malaclemys terrapin centrata) exhibit strong environmental adaptability and live in both freshwater and saltwater. However, the genetic basis of this adaptability has not been the focus of research. In this study, we successfully constructed a ∼2.21-Gb chromosome-level genome assembly for M. t. centrata using high-coverage and high-depth genomic sequencing data generated on multiple platforms. The M. t. centrata genome contains 25 chromosomes and the scaffold N50 of ∼143.75 Mb, demonstrating high continuity and accuracy. In total, 53.82% of the genome assembly was composed of repetitive sequences, and 22 435 protein-coding genes were predicted. Our phylogenetic analysis indicated that M. t. centrata was closely related to the red-eared slider turtle (Trachemys scripta elegans), with divergence approximately ∼23.6 million years ago (Mya) during the early Neogene period of the Cenozoic era. The population size of M. t. centrata decreased significantly over the past ∼14 Mya during the Cenozoic era. Comparative genomic analysis indicated that 36 gene families related to ion transport were expanded and several genes (AQP3, solute carrier subfamily, and potassium channel genes) underwent specific amino acid site mutations in the M. t. centrata genome. Changes to these ion transport-related genes may have contributed to the remarkable salinity adaptability of diamondback terrapin. The results of this study not only provide a high-quality reference genome for M. t. centrata but also elucidate the possible genetic basis for salinity adaptation in this species.

Living with a tumor: A case of osteosarcoma involving the medullary region in Phrynops cf. P. geoffroanus (Testudines: Chelidae)

The individual Geoffroy's side-necked turtle, Phrynops cf. P. geoffroanus, was diagnosed postmortem with osteosarcoma associated with the forelimb through morphological and histological analysis. Osteosarcoma stands as the most prevalent primary malignant bone tumor in tetrapods. The tumor presents itself as a large mass in the distal epiphysis, characterized by spicular outgrowths and a rugose external texture. Histologically, the afflicted humerus displayed a high degree of vascularity and exhibited an extensive bone resorption process involving the medullary and endosteal regions. Notably, a clear transition between the bone marrow and cortical bone was absent, indicative of a remodeling process featuring Haversian bone system apposition. Additionally, the diaphyseal region displayed the progression of neoplastic bone tissue along the bone. For comparative purposes, we describe a humeral thin section from a healthy specimen revealing compact primary bone interrupted by cyclical growth marks which differs from the continuous growth observed in the neoplastic humerus. To assess the neoplastic bone growth rate at the mid-diaphysis level, phylogenetic eigenvector maps (PEM) were employed, utilizing osteocyte density and vascular density as explanatory variables. The findings indicated that the osteosarcoma exhibited a slow-growing nature, suggesting that the turtle had to live with this condition for years. As the neoplasia continued to expand, it likely led to disadvantages for the pathological Phrynops individual due to humeral deformity. Furthermore, malignancy was associated with angiogenesis and the invasion of the medullary region by neoplastic bone tissue, raising the likelihood of metastasis as an additional factor contributing to the individual's sickness. The presence of numerous vascular canals in the diaphyseal thin section suggested a low-grade central osteosarcoma. It is worth noting that osseous neoplasms are rarely documented in Testudines, making this case of osteosarcoma in a South American freshwater chelid specimen a unique and rare occurrence.

Skull osteology, neuroanatomy, and jaw-related myology of the pig-nosed turtle Carettochelys insculpta (Cryptodira, Trionychia)

The osteology, neuroanatomy, and musculature are known for most primary clades of turtles (i.e., "families"), but knowledge is still lacking for one particular clade, the Carettochelyidae. Carettochelyids are represented by only one living taxon, the pig-nosed turtle Carettochelys insculpta. Here, we use micro-computed tomography of osteological and contrast-enhanced stained specimens to describe the cranial osteology, neuroanatomy, circulatory system, and jaw musculature of Carettochelys insculpta. The jaw-related myology is described in detail for the first time for this taxon, including m. zygomaticomandibularis, a muscular unit only found in trionychians. We also document a unique arterial pattern for the internal carotid artery and its subordinate branches and provide an extensive list of osteological ontogenetic differences. The present work provides new insights into the craniomandibular anatomy of turtles and will allow a better understanding of the evolutionary history of the circulatory system of trionychians and intraspecific variation among turtles.

Growth and life history of freshwater chelydrid turtles (Testudines: Cryptodira): A bone histological approach

The current study examines the growth pattern and lifestyle habits of the freshwater snapping turtles Chelydra and Macrochelys based on limb bone histology. Femora, humeri, and tibiae of 25 individuals selected from a range of ontogenetic stages were assessed to determine inter-element and intraskeletal histological variation. Osteohistological assessment of multiple elements is consistent with overall moderate growth rates as revealed by the dominance of parallel-fibered bone. However, the growth was cyclical as shown by deposition of multiple lines of arrested growths in the compacta. It appears that the bone tissue of C. serpentina is more variable through ontogeny with intermittent higher growth rates. M. temminckii appears to grow more slowly than C. serpentina possessing compact and thick cortices in accordance with their larger size. Overall, vascularization decreases through ontogeny with humeri and femora being well-vascularized in both species. Contrarily, epipodials are poorly vascularized, though simple longitudinal and radial canals are present, suggesting differences in growth patterns when compared with associated diaphyseal sections. The tibiae were found to be the least remodeled of the limb bones and therefore better suited for skeletochronology for snapping turtles. Intra-elementally, femora and humeri preserved higher cortical vascularity ventrally, suggestive of faster relative growth. We hypothesize that the differential growth pattern in limb bones of snapping turtles may relate to differential functional constraints, where forelimbs are operational in swimming while the hindlimbs provide stability.

Diversity and biogeography of South American mud turtles elucidated by multilocus DNA sequencing (Testudines: Kinosternidae)

Kinosternon is the most speciose genus of extant turtles, with 22 currently recognized species, distributed across large parts of the Americas. Most species have small distributions, but K. leucostomum and K. scorpioides range from Mexico to South America. Previous studies have found discordance between mitochondrial and nuclear phylogenies in some kinosternid groups, with the current taxonomy following the nuclear-based results. Herein, based on extended molecular, geographic, and taxonomic sampling, we explore the phylogeographic structure and taxonomic limits for K. leucostomum and the K. scorpioides group and present a fossil-calibrated nuclear time tree for Kinosternon. Our results reveal contrasting differentiation patterns for the K. scorpioides group and K. leucostomum, despite overlapping distributions. Kinosternon leucostomum shows only shallow geographic divergence, whereas the K. scorpioides group is polyphyletic with up to 10 distinct taxa, some of them undescribed. We support the elevation of K. s. albogulare and K. s. cruentatum to species level. Given the deep divergence within the genus Kinosternon, we propose the recognition of three subgenera, Kinosternon, Cryptochelys and Thyrosternum, and the abandonment of the group-based classification, at least for the K. leucostomum and K. scorpioides groups. Our results show an initial split in Kinosternon that gave rise to two main radiations, one Nearctic and one mainly Neotropical. Most speciation events in Kinosternon occurred during the Quaternary and we hypothesize that they were mediated by both climatic and geological events. Additionally, our data imply that at least three South American colonizations occurred, two in the K. leucostomum group, and one in the K. scorpioides group. Additionally, we hypothesize that discordance between mitochondrial and nuclear phylogenetic signal is due to mitochondrial capture from an extinct kinosternine lineage.

The topological organization of the turtle cranium is constrained and conserved over long evolutionary timescales

The cranium of turtles (Testudines) is characterized by the secondary reduction of temporal fenestrae and loss of cranial joints (i.e., characteristics of anapsid, akinetic skulls). Evolution and ontogeny of the turtle cranium are associated with shape changes. Cranial shape variation among Testudines can partially be explained by dietary and functional adaptations (neck retraction), but it is unclear if cranial topology shows similar ecomorphological signal, or if it is decoupled from shape evolution. We assess the topological arrangement of cranial bones (i.e., number, relative positioning, connections), using anatomical network analysis. Non-shelled stem turtles have similar cranial arrangements to archosauromorph outgroups. Shelled turtles (Testudinata) evolve a unique cranial organization that is associated with bone losses (e.g., supratemporal, lacrimal, ectopterygoid) and an increase in complexity (i.e., densely and highly interconnected skulls with low path lengths between bones), resulting from the closure of skull openings and establishment of unusual connections such as a parietal-pterygoid contact in the secondary braincase. Topological changes evolutionarily predate many shape changes. Topological variation and taxonomic morphospace discrimination among crown turtles are low, indicating that cranial topology may be constrained. Observed variation results from repeated losses of nonintegral bones (i.e., premaxilla, nasal, epipterygoid, quadratojugal), and changes in temporal emarginations and palate construction. We observe only minor ontogenetic changes. Topology is not influenced by diet and habitat, contrasting cranial shape. Our results indicate that turtles have a unique cranial topology among reptiles that is conserved after its initial establishment, and shows that cranial topology and shape have different evolutionary histories.

The meaning of mass extinctions and what the fossil record tells us about angiosperm survival at K-Pg: a reply to Hagen (2024)

Last year, we published research using phylogenetic comparative methods (PCMs) to reveal no phylogenetic evidence for elevated lineage-level extinction rates in angiosperms across K-Pg (Thompson JB, Ramírez-Barahona S. 2023 No phylogenetic evidence for angiosperm mass extinction at the Cretaceous-Palaeogene (K-Pg) boundary. Biol. Lett. 19, 20230314. (doi:10.1098/rsbl.2023.0314)), results that are in step with the global angiosperm fossil record. In a critique of our paper (Hagen ER. 2024 A critique of Thompson and Ramírez-Barahona (2023) or: how I learned to stop worrying and love the fossil record. EcoEvoRxiv. (doi:10.32942/X2631W)), simulation work is presented to argue we erred in our methodological choices and interpretations, and that we should have deferred to fossil evidence. In our opinion, underlying this critique are poor methodological choices on simulations and philosophical problems surrounding the definition of a mass extinction event, which leads to incorrect interpretations of both the fossil record and PCMs. We further argue that deferring to one source of evidence in favour of the other shuts the door to important evolutionary and philosophical questions.

Convergent evolution toward a slow pace of life predisposes insular endotherms to anthropogenic extinctions

Island vertebrates have evolved a number of morphological, physiological, and life history characteristics that set them apart from their mainland relatives. However, to date, the evolution of metabolism and its impact on the vulnerability to extinction of insular vertebrates remains poorly understood. This study used metabolic data from 2813 species of tetrapod vertebrates, including 695 ectothermic and 2118 endothermic species, to reveal that island mammals and birds evolved convergent metabolic strategies toward a slow pace of life. Insularity was associated with shifts toward slower metabolic rates and greater generation lengths in endotherms, while insularity just drove the evolution of longer generation lengths in ectotherms. Notably, a slow pace of life has exacerbated the extinction of insular endemic species in the face of anthropogenic threats. These findings have important implications for understanding physiological adaptations associated with the island syndrome and formulating conservation strategies across taxonomic groups with different metabolic modes.

Geographic variation in incubation temperatures promoting viable offspring production in broadly co-distributed turtles

Organisms whose early life stages are environmentally sensitive produce offspring within a relatively narrow range of suitable abiotic conditions. In reptiles, development rate and survival are often maximized if incubation temperatures remain under 31°C, though this upper bound may vary within and among species. We addressed this expectation by comparing responses to egg incubation at 30°C versus 33°C in congeneric turtle species pairs with broad syntopic geographic distributions. In the two softshell turtles (Apalone spp.), the greatest changes in development rate and phenotypic variance were observed in the northernmost population, which had a low survival rate (40%) at 33°C. The presumably suboptimal temperature (33°C) for northern populations otherwise yielded 76%-93% survival rates and fast swimming speeds in more southern populations. Still, in one species, northern hatchlings incubated at 33°C matched the elevated speeds of their southern counterparts, revealing a countergradient response. In northern populations of the two map turtles (Graptemys spp.), survival was also reduced (28%-60%) at 33°C and the development rate (relative to 30°C) increased by up to 75%. Our experiments on divergent taxa with similar nesting ecologies substantiate that the optimal thermal range for offspring production is variable. These findings encourage further work on how population- and species-level differences relate to local adaptation in widely distributed oviparous species.

Two Major Extinction Events in the Evolutionary History of Turtles: One Caused by an Asteroid, the Other by Hominins

AbstractWe live in a time of accelerated biological extinctions that has the potential to mirror past mass extinction events. However, the rarity of mass extinctions and the restructuring of diversity they cause complicate direct comparisons between the current extinction crisis and earlier events. Among animals, turtles (Testudinata) are one of few groups that have both a rich fossil record and sufficiently stable ecological and functional roles to enable meaningful comparisons between the end-Cretaceous mass extinction (∼66 Ma) and the ongoing wave of extinctions. Here we analyze the fossil record of the entire turtle clade and identify two peaks in extinction rates over their evolutionary history. The first coincides with the Cretaceous-Paleogene transition, reflecting patterns previously reported for other taxa. The second major extinction event started in the Pliocene and continues until now. This peak is detectable only for terrestrial turtles and started much earlier in Africa and Eurasia than elsewhere. On the basis of the timing, geography, and functional group of this extinction event, we postulate a link to co-occurring hominins rather than climate change as the cause. These results lend further support to the view that negative biodiversity impacts were already incurred by our ancestors and related lineages and demonstrate the severity of this continued impact through human activities.

The cranial and postcranial morphology of Hutchemys rememdium and its impact on the phylogenetic relationships of Plastomenidae (Testudinata, Trionychidae)

Hutchemys rememdium is a poorly understood softshell turtle (Trionychidae) from the mid Paleocene of the Williston Basin of North America previously known only from postcranial remains. A particularly rich collection of previously undescribed material from the Tiffanian 4 North American Land Mammal Age (NALMA) of North Dakota is here presented consisting of numerous shells that document new variation, some non-shell postcrania, and cranial remains, which are described based on 3D models extracted from micro-CT data. Although the observed shell variation weakens previously noted differences with the younger species Hutchemys arctochelys from the Clarkforkian NALMA, the two taxa are still recognized as distinct. Parsimony and Bayesian phylogenetic analyses reaffirm the previously challenged placement of Hutchemys rememdium within the clade Plastomenidae, mostly based on novel observations of cranial characters made possible by the new material and the micro-CT data. The new topology supports the notion that the well-ossified plastron of plastomenids originated twice in parallel near the Cretaceous/Paleogene boundary, once in the Hutchemys lineage and once in the Gilmoremys/Plastomenus lineage. Hutchemys rememdium is notable for being the only documented species of trionychid in the mid Paleocene of the Williston Basin. The presence of multiple individuals in a carbonaceous claystone indicates this taxon lived in swamps and lakes and its expanded triturating surface suggests it had a durophagous diet.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13358-024-00315-8.

Morphological Diversity of Turtle Hyoid Apparatus is Linked to Feeding Behavior

The hyoid apparatus of tetrapods is highly diverse in its morphology. It plays an important role in feeding, breathing, sound production, and various other behaviors. Among turtles, the diversity of the hyoid apparatus has been recurrently linked to their habitat. The ossification of the hyoid corpus is often the main trait used in correlations with "niche" occupancy, an ossified corpus being associated with aquatic environments and a cartilaginous corpus with terrestrial life. Most studies conducted so far have focused on species belonging to Testudinoidea, the clade that occupies the biggest diversity of habitats (i.e., terrestrial, semi-terrestrial, and aquatic animals), while other turtle lineages have been largely understudied. We assessed the adult anatomy of the hyoid apparatus of 92 turtle species from all "families", together with ossification sequences from embryological series of 11 species, some described for the first time here. Using nearly 40 different discrete anatomical characters, we discuss the evolutionary patterns and the biological significance of morphological transformations in the turtle hyoid elements. Morphological changes are strongly associated to feeding modes, with several instances of convergent evolution within and outside the Testudines clade, and are not as strongly connected to habitat as previously thought. Some of the hyoid character states we describe are diagnostic of specific turtle clades, thus providing phylogenetically relevant information.

Multiple horizontal transfer events of a DNA transposon into turtles, fishes, and a frog

Horizontal transfer of transposable elements (HTT) has been reported across many species and the impact of such events on genome structure and function has been well described. However, few studies have focused on reptilian genomes, especially HTT events in Testudines (turtles). Here, as a consequence of investigating the repetitive content of Malaclemys terrapin terrapin (Diamondback turtle) we found a high similarity DNA transposon, annotated in RepBase as hAT-6_XT, shared between other turtle species, ray-finned fishes, and a frog. hAT-6_XT was notably absent in reptilian taxa closely related to turtles, such as crocodiles and birds. Successful invasion of DNA transposons into new genomes requires the conservation of specific residues in the encoded transposase, and through structural analysis, these residues were identified indicating some retention of functional transposition activity. We document six recent independent HTT events of a DNA transposon in turtles, which are known to have a low genomic evolutionary rate and ancient repeats.

Digital skull anatomy of the Oligocene North American tortoise Stylemys nebrascensis with taxonomic comments on the species and comparisons with extant testudinids of the Gopherus-Manouria clade

The anatomy of North American tortoises is poorly understood, despite a rich fossil record from the Eocene and younger strata. Stylemys nebrascensis is a particularly noteworthy turtle in this regard, as hundreds of specimens are known from Oligocene deposits, and as this species is one of the earliest fossil turtles to have been described in the scientific literature. Since its initial description based on a shell, many specimens with more complete material have been referred to Stylemys nebrascensis. Here, we review and confirm the referral of an important historic specimen to Stylemys nebrascensis, which includes shell, non-shell postcranial, and skull material. This allows us to document unique skull features of Stylemys nebrascensis (e.g., an unusual 'poststapedial canal' that connects the posterior skull surface with the cavum acustico-jugulare) and to refer another well-preserved skull to the species. Based on computed-tomography scanning of these two skulls, we provide a detailed description of the cranial and mandibular osteology of Stylemys nebrascensis. Stylemys nebrascensis has a combination of plesiomorphic skull characteristics (e.g., retention of a medial jugal process) and derived traits shared with extant gopher tortoises (e.g., median premaxillary ridge) that suggest it may be a stem-representative of the gopher tortoise lineage. This supports the hypothesis that extant and fossil tortoises from North America form a geographically restricted clade that split from Asian relatives during the Paleogene.

The latest freshwater giants: a new Peltocephalus (Pleurodira: Podocnemididae) turtle from the Late Pleistocene of the Brazilian Amazon

Overkill of large mammals is recognized as a key driver of Pleistocene megafaunal extinctions in the Americas and Australia. While this phenomenon primarily affected mega-mammals, its impact on large Quaternary reptiles has been debated. Freshwater turtles, due to the scarcity of giant forms in the Quaternary record, have been largely neglected in such discussions. Here we present a new giant podocnemidid turtle, Peltocephalus maturin sp. nov., from the Late Pleistocene Rio Madeira Formation in the Brazilian Amazon, that challenges this assumption. Morphological and phylogenetic analyses of the holotype, a massive partial lower jaw, reveal close affinities to extant Amazonian species and suggest an omnivorous diet. Body size regressions indicate Pe. maturin possibly reached about 180 cm in carapace length and is among the largest freshwater turtles ever found. This finding presents the latest known occurrence of giant freshwater turtles, hinting at coexistence with early human inhabitants in the Amazon.

ReptTraits: a comprehensive dataset of ecological traits in reptiles

Trait datasets are increasingly being used in studies investigating eco-evolutionary theory and global conservation initiatives. Reptiles are emerging as a key group for studying these questions because their traits are crucial for understanding the ability of animals to cope with environmental changes and their contributions to ecosystem processes. We collected data from earlier databases, and the primary literature to create an up-to-date dataset of reptilian traits, encompassing 40 traits from 12060 species of reptiles (Archelosauria: Crocodylia and Testudines, Rhynchocephalia, and Squamata: Amphisbaenia, Sauria, and Serpentes). The data were gathered from 1288 sources published between 1820 and 2023. The dataset includes morphological, physiological, behavioral, and life history traits, as well as information on the availability of genetic data, IUCN Red List assessments, and population trends.

Diversity and evolution of the vertebrate chemoreceptor gene repertoire

Chemoreception - the ability to smell and taste - is an essential sensory modality of most animals. The number and type of chemical stimuli that animals can perceive depends primarily on the diversity of chemoreceptors they possess and express. In vertebrates, six families of G protein-coupled receptors form the core of their chemosensory system, the olfactory/pheromone receptor gene families OR, TAAR, V1R and V2R, and the taste receptors T1R and T2R. Here, we study the vertebrate chemoreceptor gene repertoire and its evolutionary history. Through the examination of 1,527 vertebrate genomes, we uncover substantial differences in the number and composition of chemoreceptors across vertebrates. We show that the chemoreceptor gene families are co-evolving, highly dynamic, and characterized by lineage-specific expansions (for example, OR in tetrapods; TAAR, T1R in teleosts; V1R in mammals; V2R, T2R in amphibians) and losses. Overall, amphibians, followed by mammals, are the vertebrate clades with the largest chemoreceptor repertoires. While marine tetrapods feature a convergent reduction of chemoreceptor numbers, the number of OR genes correlates with habitat in mammals and birds and with migratory behavior in birds, and the taste receptor repertoire correlates with diet in mammals and with aquatic environment in fish.

Decoupling speciation and extinction reveals both abiotic and biotic drivers shaped 250 million years of diversity in crocodile-line archosaurs

Whereas living representatives of Pseudosuchia, crocodylians, number fewer than 30 species, more than 700 pseudosuchian species are known from their 250-million-year fossil record, displaying far greater ecomorphological diversity than their extant counterparts. With a new time-calibrated tree of >500 species, we use a phylogenetic framework to reveal that pseudosuchian evolutionary history and diversification dynamics were directly shaped by the interplay of abiotic and biotic processes over hundreds of millions of years, supported by information theory analyses. Speciation, but not extinction, is correlated with higher temperatures in terrestrial and marine lineages, with high sea level associated with heightened extinction in non-marine taxa. Low lineage diversity and increased speciation in non-marine species is consistent with opportunities for niche-filling, whereas increased competition may have led to elevated extinction rates. In marine lineages, competition via increased lineage diversity appears to have driven both speciation and extinction. Decoupling speciation and extinction, in combination with ecological partitioning, reveals a more complex picture of pseudosuchian evolution than previously understood. As the number of species threatened with extinction by anthropogenic climate change continues to rise, the fossil record provides a unique window into the drivers that led to clade success and those that may ultimately lead to extinction.

Turtle skull development unveils a molecular basis for amniote cranial diversity

Amniote skulls display diverse architectural patterns including remarkable variations in the number of temporal arches surrounding the upper and lower temporal fenestrae. However, the cellular and molecular basis underlying this diversification remains elusive. Turtles are a useful model to understand skull diversity due to the presence of secondarily closed temporal fenestrae and different extents of temporal emarginations (marginal reduction of dermal bones). Here, we analyzed embryos of three turtle species with varying degrees of temporal emargination and identified shared widespread coexpression of upstream osteogenic genes Msx2 and Runx2 and species-specific expression of more downstream osteogenic genes Sp7 and Sparc in the head. Further analysis of representative amniote embryos revealed differential expression patterns of osteogenic genes in the temporal region, suggesting that the spatiotemporal regulation of Msx2, Runx2, and Sp7 distinguishes the temporal skull morphology among amniotes. Moreover, the presence of Msx2- and/or Runx2-positive temporal mesenchyme with osteogenic potential may have contributed to their extremely diverse cranial morphology in reptiles.

Turtle Shell Kinesis Underscores Constraints and Opportunities in the Evolution of the Vertebrate Musculoskeletal System

Species groups that feature traits with a low number of potentially variable (evolvable) character states are more likely to repeatedly evolve similar phenotypes, that is, convergence. To evaluate this phenomenon, this present paper addresses anatomical alterations in turtles that convergently evolved shell kinesis, for example, the movement of shell bones to better shield the head and extremities. Kinesis constitutes a major departure from the evolutionarily conserved shell of modern turtles, yet it has arisen independently at least 8 times. The hallmark signature of kinesis is the presence of shell bone articulations or "hinges," which arise via similar skeletal remodeling processes in species that do not share a recent common ancestor. Still, the internal biomechanical components that power kinesis may differ in such distantly related species. Complex diarthrodial joints and modified muscle connections expand the functional boundaries of the limb girdles and neck in a lineage-specific manner. Some lineages even exhibit mobility of thoracic and sacral vertebrae to facilitate shell closure. Depending on historical contingency and structural correlation, a myriad of anatomical alterations has yielded similar functional outcomes, that is, many-to-one mapping, during the convergent evolution of shell kinesis. The various iterations of this intricate phenotype illustrate the potential for the vertebrate musculoskeletal system to undergo evolutionary change, even when constraints are imposed by the development and structural complexity of a shelled body plan. Based on observations in turtles and comparisons to other vertebrates, a hypothetical framework that implicates functional interactions in the origination of novel musculoskeletal traits is presented.

Evolution of posture in amniotes-Diving into the trabecular architecture of the femoral head

Extant amniotes show remarkable postural diversity. Broadly speaking, limbs with erect (strongly adducted, more vertically oriented) posture are found in mammals that are particularly heavy (graviportal) or show good running skills (cursorial), while crouched (highly flexed) limbs are found in taxa with more generalized locomotion. In Reptilia, crocodylians have a "semi-erect" (somewhat adducted) posture, birds have more crouched limbs and lepidosaurs have sprawling (well-abducted) limbs. Both synapsids and reptiles underwent a postural transition from sprawling to more erect limbs during the Mesozoic Era. In Reptilia, this postural change is prominent among archosauriforms in the Triassic Period. However, limb posture in many key Triassic taxa remains poorly known. In Synapsida, the chronology of this transition is less clear, and competing hypotheses exist. On land, the limb bones are subject to various stresses related to body support that partly shape their external and internal morphology. Indeed, bone trabeculae (lattice-like bony struts that form the spongy bone tissue) tend to orient themselves along lines of force. Here, we study the link between femoral posture and the femoral trabecular architecture using phylogenetic generalized least squares. We show that microanatomical parameters measured on bone cubes extracted from the femoral head of a sample of amniote femora depend strongly on body mass, but not on femoral posture or lifestyle. We reconstruct ancestral states of femoral posture and various microanatomical parameters to study the "sprawling-to-erect" transition in reptiles and synapsids, and obtain conflicting results. We tentatively infer femoral posture in several hypothetical ancestors using phylogenetic flexible discriminant analysis from maximum likelihood estimates of the microanatomical parameters. In general, the trabecular network of the femoral head is not a good indicator of femoral posture. However, ancestral state reconstruction methods hold great promise for advancing our understanding of the evolution of posture in amniotes.

A late-surviving phytosaur from the northern Atlantic rift reveals climate constraints on Triassic reptile biogeography

BACKGROUND

The origins of all major living reptile clades, including the one leading to birds, lie in the Triassic. Following the largest mass extinction in Earth's history at the end of the Permian, the earliest definite members of the three major living reptile clades, the turtles (Testudines), crocodylians and birds (Archosauria), and lizards, snakes, amphisbaenians, and Tuatara (Lepidosauria) appeared. Recent analyses of the Triassic reptile fossil record suggest that the earliest diversifications in all three of these clades were tightly controlled by abrupt paleoclimate fluctuations and concordant environmental changes. Yet, this has only been preliminarily tested using information from evolutionary trees. Phytosauria consists of superficially crocodylian-like archosaurs that either form the sister to the crown or are the earliest divergence on the crocodylian stem and are present throughout the Triassic, making this clade an excellent test case for examining this biogeographic hypothesis.

RESULTS

Here, I describe a new phytosaur, Jupijkam paleofluvialis gen. et sp. nov., from the Late Triassic of Nova Scotia, Canada, which at that time sat in northern Pangaea near the northern terminus of the great central Pangean rift. As one of the northernmost occurrences of Phytosauria, J. paleofluvialis provides critical new biogeographic data that enables revised estimations of phytosaur historical biogeography along phylogenies of this clade built under multiple methodologies. Reconstructions of phytosaur historical biogeography based on different phylogenies and biogeographic models suggest that phytosaurs originated in northern Pangaea, spread southward, and then dispersed back northward at least once more during the Late Triassic.

CONCLUSIONS

The results presented in this study link phytosaur biogeography to major changes to Triassic global climate and aridity. Together with the earliest dinosaurs and several other reptile lineages, phytosaur diversification and migration appear to have been restricted by the formation and loss of arid belts across the Pangean supercontinent.

Somatic Growth and Maturity for Four Species of River Cooter Including Pseudemys concinna suwanniensis, P. nelsoni, P. peninsularis, and P. texana

Pseudemys is a genus of commonly occurring freshwater turtles with limited growth information across their long lifespans. We used 11,361 mark-recapture events to estimate the somatic growth rates of P. nelsoni, P. peninsularis, P. concinna suwanniensis, and P. texana from freshwater springs and developed a Bayesian growth model to estimate the species-specific, site-specific, and individual effects on growth. We corroborated evidence for fast juvenile growth and slower adult growth in Pseudemys but found uncommonly fast growth rates, with turtles doubling or tripling in size in the first year. P. texana males had the smallest average maximum size (L∞, 243 mm), while P. c. suwanniensis females had the largest (423 mm). Environmental conditions at springs had significant effects on k, the growth coefficient, but not L∞. We derived, using a ratio of length at maturity to L∞ (71.7% and 87%, males and females), that females matured 1.15-1.57 times older than males except for P. c. suwanniensis, which matured three times older. Given the local abundance declines in many Pseudemys from anthropogenic impacts, this study provides important baseline life history information for Pseudemys species for use in ongoing conservation efforts and presents a novel hierarchical modeling approach using a long-term mark-recapture dataset.

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.

Low coverage sequencing provides insights into the key features of the nuclear and mitochondrial genomes of the Alligator Snapping Turtle Macrochelys temminckii

The alligator snapping turtle Macrochelys temminckii is a culturally, ecologically, and evolutionary relevant species of conservation concern. In this study, we conducted a genome survey of M. temminckii. Using a low-coverage short read sequencing strategy, this study estimated the genome size, repetitive genome content, annotated and quantified repetitive elements, assembled the 45S rRNA DNA operon, and characterized in detail the mitochondrial genome of M. temminckii. Using a k-mer strategy, the estimated haploid genome size varied between 3.77 and 3.19 Gbp, which is within the range previously reported for other representatives of the family Chelydridae. Repetitive genome content estimates using different k-mers (21 to 51) indicated that more than 75% of the genome of M. temminckii comprised repetitive elements. Taking into account only annotated repetitive elements, the most common repetitive elements were classified as Class I - Long Interspersed Nuclear Element (LINE) which were more abundant than Class I - Penelope and Class I - Long Terminal Repeat (LTR) Ty3-gypsy mobile elements. Less abundant repeat element families in the nuclear genome of M. temminckii included Class I - DIRS mobile elements and Satellite DNA. The nuclear ribosomal operon was partially assembled into two contigs, one encoding the complete ssrDNA and a second comprising the full lsrDNA. The AT-rich complete mitochondrial genome was 16,570 bp long. These new genomic resources are of utmost importance to aid in the development of conservation plans for this freshwater turtle.

Genomics reveals broad hybridization in deeply divergent Palearctic grass and water snakes (Natrix spp.)

Understanding speciation is one of the cornerstones of biological diversity research. Currently, speciation is often understood as a continuous process of divergence that continues until genetic or other incompatibilities minimize or prevent interbreeding. The Palearctic snake genus Natrix is an ideal group to study speciation, as it comprises taxa representing distinct stages of the speciation process, ranging from widely interbreeding parapatric taxa through parapatric species with very limited gene flow in narrow hybrid zones to widely sympatric species. To understand the evolution of reproductive isolation through time, we have sequenced the genomes of all five species within this genus and two additional subspecies. We used both long-read and short-read methods to sequence and de-novo-assemble two high-quality genomes (Natrix h. helvetica, Natrix n. natrix) to their 1.7 Gb length with a contig N50 of 4.6 Mbp and 1.5 Mbp, respectively, and used these as references to assemble the remaining short-read-based genomes. Our phylogenomic analyses yielded a well-supported dated phylogeny and evidence for a surprisingly complex history of interspecific gene flow, including between widely sympatric species. Furthermore, evidence for gene flow was also found for currently allopatric species pairs. Genetic exchange among these well-defined, distinct, and several million-year-old reptile species emphasizes that speciation and maintenance of species distinctness can occur despite continued genetic exchange.

Cranial and mandibular anatomy of Plastomenus thomasii and a new time-tree of trionychid evolution

Trionychid (softshell) turtles have a peculiar bauplan, which includes shell reductions and cranial elongation. Despite a rich fossil record dating back to the Early Cretaceous, the evolutionary origin of the trionychid bauplan is poorly understood, as even old fossils show great anatomical similarities to extant species. Documenting structural detail of fossil trionychids may help resolve the evolutionary history of the group. Here, we study the cranial and mandibular anatomy of Plastomenus thomasii using µCT scanning. Plastomenus thomasii belongs to the Plastomenidae, a long-lived (Santonian-Eocene) clade with uncertain affinities among trionychid subclades. The skulls of known plastomenids are characterized by unusual features otherwise not known among trionychids, such as extremely elongated, spatulate mandibular symphyses. We use anatomical observations for updated phylogenetic analyses using both parsimony and Bayesian methods. There is strong support across methods for stem-cyclanorbine affinities for plastomenids. The inclusion of stratigraphic data in our Bayesian analysis indicates that a range of Cretaceous Asian fossils including Perochelys lamadongensis may be stem-trionychids, suggesting that many features of trionychid anatomy evolved prior to the appearance of the crown group. Divergence time estimates from Bayesian tip-dating for the origin of crown Trionychia (134.0 Ma) and Pan-Trionychidae (123.8 Ma) constrain the evolutionary time span during which the trionychid bauplan has evolved to a range of < 11 million years. Bayesian rate estimation implies high morphological rates during early softshell turtle evolution. If correct, plastomenids partially fill the stratigraphic gap which results from shallow divergence times of crown cyclanorbines during the late Eocene.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13358-023-00267-5.

A quantum leap in the evolution of platyhelminths: host-switching from turtles to hippopotamuses illustrated from a phylogenetic meta-analysis of polystomes (Monogenea, Polystomatidae)

While monogenean worms are mainly parasites of the gills and skin of fish, and to a lesser extent parasites of the oral cavity, urinary bladder, and/or conjunctival sacs of amphibians and freshwater turtles, Oculotrema hippopotami Stunkard, 1924 is the single monogenean polystome reported from a mammal, the common hippopotamus (Hippopotamus amphibius Linnaeus). Several hypotheses have been suggested in the last decade to explain the origin of this enigmatic parasite which infects the conjunctival sacs of H. amphibius. Based on a molecular phylogeny inferred from nuclear (28S and 18S) and mitochondrial (12S and COI) sequences of O. hippopotami and chelonian polystomes, we found a sister group relationship between O. hippopotami and Apaloneotrema moleri (Du Preez & Morrison, 2012). This result suggests lateral parasite transfer between freshwater turtles and hippopotamuses, thus likely reflecting one of the most exceptional known examples of host-switching in the course of vertebrate evolution. It also demonstrates that the proximity in the ecological habitat of parasites within host species is an important feature for their speciation and diversification. Because A. moleri and its host, the Florida softshell turtle (Apalone ferox (Schneider)), are restricted to the USA, we suggest that an ancestral stock of parasites may have been isolated on primitive African trionychids after they diverged from their American relatives, and then switched to hippopotamuses or anthracotheres in Africa.

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.

Locomotor and postural diversity among reptiles viewed through the prism of femoral microanatomy: Palaeobiological implications for some Permian and Mesozoic taxa

The water-to-land transition by the first tetrapod vertebrates represents a key stage in their evolution. Selection pressures exerted by this new environment on animals led to the emergence of new locomotor and postural strategies that favoured access to different ecological niches and contributed to their evolutionary success. Today, amniotes show great locomotor and postural diversity, particularly among Reptilia, whose extant representatives include parasagittally locomoting erect and crouched bipeds (birds), sub-parasagittal 'semi-erect' quadrupeds (crocodylians) and sprawling quadrupeds (squamates and turtles). But the different steps leading to such diversity remain enigmatic and the type of locomotion adopted by many extinct species raises questions. This is notably the case of certain Triassic taxa such as Euparkeria and Marasuchus. The exploration of the bone microanatomy in reptiles could help to overcome these uncertainties. Indeed, this locomotor and postural diversity is accompanied by great microanatomical disparity. On land, the bones of the appendicular skeleton support the weight of the body and are subject to multiple constraints that partly shape their external and internal morphology. Here we show how microanatomical parameters measured in cross-section, such as bone compactness or the position of the medullocortical transition, can be related to locomotion. We hypothesised that this could be due to variations in cortical thickness. Using statistical methods that take phylogeny into account (phylogenetic flexible discriminant analyses), we develop different models of locomotion from a sample of femur cross-sections from 51 reptile species. We use these models to infer locomotion and posture in 7 extinct reptile taxa for which they remain debated or not fully clear. Our models produced reliable inferences for taxa that preceded and followed the quadruped/biped and sprawling/erect transitions, notably within the Captorhinidae and Dinosauria. For taxa contemporary with these transitions, such as Terrestrisuchus and Marasuchus, the inferences are more questionable. We use linear models to investigate the effect of body mass and functional ecology on our inference models. We show that body mass seems to significantly impact our model predictions in most cases, unlike the functional ecology. Finally, we illustrate how taphonomic processes can impact certain microanatomical parameters, especially the eccentricity of the section, while addressing some other potential limitations of our methods. Our study provides insight into the evolution of enigmatic locomotion in various early reptiles. Our models and methods could be used by palaeontologists to infer the locomotion and posture in other extinct reptile taxa, especially when considered in combination with other lines of evidence.

Divergent sensory and immune gene evolution in sea turtles with contrasting demographic and life histories

Sea turtles represent an ancient lineage of marine vertebrates that evolved from terrestrial ancestors over 100 Mya. The genomic basis of the unique physiological and ecological traits enabling these species to thrive in diverse marine habitats remains largely unknown. Additionally, many populations have drastically declined due to anthropogenic activities over the past two centuries, and their recovery is a high global conservation priority. We generated and analyzed high-quality reference genomes for the leatherback (Dermochelys coriacea) and green (Chelonia mydas) turtles, representing the two extant sea turtle families. These genomes are highly syntenic and homologous, but localized regions of noncollinearity were associated with higher copy numbers of immune, zinc-finger, and olfactory receptor (OR) genes in green turtles, with ORs related to waterborne odorants greatly expanded in green turtles. Our findings suggest that divergent evolution of these key gene families may underlie immunological and sensory adaptations assisting navigation, occupancy of neritic versus pelagic environments, and diet specialization. Reduced collinearity was especially prevalent in microchromosomes, with greater gene content, heterozygosity, and genetic distances between species, supporting their critical role in vertebrate evolutionary adaptation. Finally, diversity and demographic histories starkly contrasted between species, indicating that leatherback turtles have had a low yet stable effective population size, exhibit extremely low diversity compared with other reptiles, and harbor a higher genetic load compared with green turtles, reinforcing concern over their persistence under future climate scenarios. These genomes provide invaluable resources for advancing our understanding of evolution and conservation best practices in an imperiled vertebrate lineage.

Ancient DNA elucidates the lost world of western Indian Ocean giant tortoises and reveals a new extinct species from Madagascar

Before humans arrived, giant tortoises occurred on many western Indian Ocean islands. We combined ancient DNA, phylogenetic, ancestral range, and molecular clock analyses with radiocarbon and paleogeographic evidence to decipher their diversity and biogeography. Using a mitogenomic time tree, we propose that the ancestor of the extinct Mascarene tortoises spread from Africa in the Eocene to now-sunken islands northeast of Madagascar. From these islands, the Mascarenes were repeatedly colonized. Another out-of-Africa dispersal (latest Eocene/Oligocene) produced on Madagascar giant, large, and small tortoise species. Two giant and one large species disappeared c. 1000 to 600 years ago, the latter described here as new to science using nuclear and mitochondrial DNA. From Madagascar, the Granitic Seychelles were colonized (Early Pliocene) and from there, repeatedly Aldabra (Late Pleistocene). The Granitic Seychelles populations were eradicated and later reintroduced from Aldabra. Our results underline that integrating ancient DNA data into a multi-evidence framework substantially enhances the knowledge of the past diversity of island faunas.

100 million years of turtle paleoniche dynamics enable the prediction of latitudinal range shifts in a warming world

Past responses to environmental change provide vital baseline data for estimating the potential resilience of extant taxa to future change. Here, we investigate the latitudinal range contraction that terrestrial and freshwater turtles (Testudinata) experienced from the Late Cretaceous to the Paleogene (100.5-23.03 mya) in response to major climatic changes. We apply ecological niche modeling (ENM) to reconstruct turtle niches, using ancient and modern distribution data, paleogeographic reconstructions, and the HadCM3L climate model to quantify their range shifts in the Cretaceous and late Eocene. We then use the insights provided by these models to infer their probable ecological responses to future climate scenarios at different representative concentration pathways (RCPs 4.5 and 8.5 for 2100), which project globally increased temperatures and spreading arid biomes at lower to mid-latitudes. We show that turtle ranges are predicted to expand poleward in the Northern Hemisphere, with decreased habitat suitability at lower latitudes, inverting a trend of latitudinal range contraction that has been prevalent since the Eocene. Trionychids and freshwater turtles can more easily track their niches than Testudinidae and other terrestrial groups. However, habitat destruction and fragmentation at higher latitudes will probably reduce the capability of turtles and tortoises to cope with future climate changes.

A Solemys Skull from the Late Cretaceous of Southern France

The first skull of Solemys (Testudines: Helochlydridae) is reported from the Late Cretaceous (Rognacian) Bastide Neuve locality, Fox Amphoux, Var, France. It is assigned to Solemys gaudryi (Matheron, 1869) on the basis of associated shell elements. Our study provides new insights regarding the skull morphology of the family Helochelydridae and suggests that Helochelydra from England and Naomichelys from North America appear to be closer to each other than to Solemys.

Mitochondrial DNA and Distribution Modelling Evidenced the Lost Genetic Diversity and Wild-Residence of Star Tortoise, Geochelone elegans (Testudines: Testudinidae) in India

The Indian star tortoise (Geochelone elegans) is a massively traded animal in South Asia. To mitigate this risk, the conservation agencies recommended guidelines to safeguard this charismatic species in nature. We adopted mitochondrial DNA-based investigation and performed species distribution modeling of G. elegans throughout its distribution range in the Indian subcontinent. The genetic analyses revealed weak genetic landscape shape interpolations, low intraspecific distances (0% to 1.5%) with mixed haplotype diversity, and a single molecular operational taxonomic unit (MOTU) in the cytochrome b gene dataset. The star tortoise, G. elegans, and its sister species Geochelone platynota showed a monophyletic clustering in the Bayesian (BA) phylogeny. We also attempt to understand the habitat suitability and quality of G. elegans in its distribution range. Our results suggest that, out of the extant area, only 56,495 km² (9.90%) is suitable for this species, with regions of highest suitability in Sri Lanka. Comparative habitat quality estimation suggests the patch shape complexity and habitat fragmentation are greater in the western and southern ranges of India, which have been greatly influenced by an increased level of urbanization and agriculture practices. We have also provided a retrospect on the potential threat to G. elegans related to the wildlife trade on the regional and international spectrum. Our results detected multiple trading hubs and junctions overlying within the suitable ranges which need special attention in the vicinity. The present study calls for a proper conservation strategy to combat the fragmented distribution and explicitly recommends intensive genetic screening of founder individuals or isolated adult colonies, implementing scientific breeding, and subsequent wild release to restore the lost genetic diversity of star tortoises.