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

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.

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.

Differential Conservation and Loss of Chicken Repeat 1 (CR1) Retrotransposons in Squamates Reveal Lineage-Specific Genome Dynamics Across Reptiles

Transposable elements (TEs) are repetitive DNA sequences which create mutations and generate genetic diversity across the tree of life. In amniote vertebrates, TEs have been mainly studied in mammals and birds, whose genomes generally display low TE diversity. Squamates (Order Squamata; including ∼11,000 extant species of lizards and snakes) show as much variation in TE abundance and activity as they do in species and phenotypes. Despite this high TE activity, squamate genomes are remarkably uniform in size. We hypothesize that novel, lineage-specific genome dynamics have evolved over the course of squamate evolution. To understand the interplay between TEs and host genomes, we analyzed the evolutionary history of the chicken repeat 1 (CR1) retrotransposon, a TE family found in most tetrapod genomes which is the dominant TE in most reptiles. We compared 113 squamate genomes to the genomes of turtles, crocodilians, and birds and used ancestral state reconstruction to identify shifts in the rate of CR1 copy number evolution across reptiles. We analyzed the repeat landscapes of CR1 in squamate genomes and determined that shifts in the rate of CR1 copy number evolution are associated with lineage-specific variation in CR1 activity. We then used phylogenetic reconstruction of CR1 subfamilies across amniotes to reveal both recent and ancient CR1 subclades across the squamate tree of life. The patterns of CR1 evolution in squamates contrast other amniotes, suggesting key differences in how TEs interact with different host genomes and at different points across evolutionary history.

Macro- and micro-anatomical investigation of the oropharyngeal roof of landform greek tortoise (Testudo graeca graeca) and semi-aquatic red-eared slider turtle (Trachemys scripta elegans)

The present investigation examined the oropharyngeal roof of two turtles having different feeding behaviors: the landform Greek tortoise (Testudo graeca graeca) primarily herbivores and the semi-aquatic red-eared slider turtle (Trachemys scripta elegans) lives in freshwater that opportunistic omnivorous grossly and by scanning and light microscopes. Grossly, the Greek tortoise had a V-shaped roof consisting of the upper rhamphotheca, peri-palatine region, upper alveolar ridge, peripheral palatine ridge, median palatine ridge, vomer, choanae, caudal palatine part, and pharynx. At the same time, the red-eared slider had a semilunar roof consisting of upper rhamphotheca, two peripheral palatine ridges, core of palatine ridges, upper alveolar band, vomer, choanae, caudal palatine part, and pharynx. SEM revealed that the red-eared slider roof appeared more straightforward. The upper rhamphotheca is sharp, with a median premaxillary notch in the red-eared slider that gives a powerful bite for cutting to compensate absence of the teeth. Additionally, the red-eared slider's upper alveolar band is interrupted by a single upper alveolar ridge that appears spiky, pointed, and longer as it needs powerful chewing of prey and there are two types of teeth-like projections at its peri-palatine area for food-crushing and chewing. The Greek tortoise palatine region had numerous ridges and folds to provide roughness for food processing. Greek tortoises had small-sized choanae with two choanal folds to minimize choanal openings when eating dusty grasses. Histologically, Greek tortoise palate was rostrally thicker and more keratinized than caudally, and the caudal palatine region was characterized by a single pair of circumvallate-like papilla with multiple mucous openings and secretions, while red-eared slider palate was slightly keratinized at the peri-choanal region, and the rest of the palate was non-keratinized with few mucous openings. The current investigation found various structural oropharyngeal roof adaptations to feeding behavior in the omnivore red-eared slide compared to the herbivorous Greek turtle.

Differential Conservation and Loss of CR1 Retrotransposons in Squamates Reveals Lineage-Specific Genome Dynamics across Reptiles

Transposable elements (TEs) are repetitive DNA sequences which create mutations and generate genetic diversity across the tree of life. In amniotic vertebrates, TEs have been mainly studied in mammals and birds, whose genomes generally display low TE diversity. Squamates (Order Squamata; ~11,000 extant species of lizards and snakes) show as much variation in TE abundance and activity as they do in species and phenotypes. Despite this high TE activity, squamate genomes are remarkably uniform in size. We hypothesize that novel, lineage-specific dynamics have evolved over the course of squamate evolution to constrain genome size across the order. Thus, squamates may represent a prime model for investigations into TE diversity and evolution. To understand the interplay between TEs and host genomes, we analyzed the evolutionary history of the CR1 retrotransposon, a TE family found in most tetrapod genomes. We compared 113 squamate genomes to the genomes of turtles, crocodilians, and birds, and used ancestral state reconstruction to identify shifts in the rate of CR1 copy number evolution across reptiles. We analyzed the repeat landscapes of CR1 in squamate genomes and determined that shifts in the rate of CR1 copy number evolution are associated with lineage-specific variation in CR1 activity. We then used phylogenetic reconstruction of CR1 subfamilies across amniotes to reveal both recent and ancient CR1 subclades across the squamate tree of life. The patterns of CR1 evolution in squamates contrast other amniotes, suggesting key differences in how TEs interact with different host genomes and at different points across evolutionary history.

The tongue of the red-eared slider (Trachemys scripta elegans): morphological characterization through gross, light, scanning electron, and immunofluorescence microscopic examination

The red-eared slider (Trachemys scripta elegans) is renowned for its remarkable adaptations, yet much of its complex biology remains unknown. In this pioneering study, we utilized a combination of gross anatomy, scanning electron microscopy (SEM), light microscopy, and immunofluorescence techniques to examine the tongue's omnivorous adaptation in this species. This research bridges a critical knowledge gap, enhancing our understanding of this intriguing reptile. Gross examination revealed a unique arrowhead-shaped tongue with a median lingual fissure and puzzle-piece-shaped tongue papillae. SEM unveiled rectangular filiform, conical, and fungiform papillae, with taste pores predominantly on the dorsal surface and mucous cells on the lateral surface of the papillae. Histologically, the tongue's apex featured short rectangular filiform and fungiform papillae, while the body exhibited varying filiform shapes and multiple taste buds on fungiform papillae. The tongue's root contained lymphatic tissue with numerous lymphocytes surrounding the central crypt, alongside lingual skeletal musculature, blood and lymph vessels, and Raffin corpuscles in the submucosa. The lingual striated muscle bundles had different orientations, and the lingual hyaline cartilage displayed a bluish coloration of the ground substance, along with a characteristic isogenous group of chondrocytes. Our research represents the first comprehensive application of immunofluorescence techniques to investigate the cellular intricacies of the red-eared slider's tongue by employing seven distinct antibodies, revealing a wide array of compelling and significant findings. Vimentin revealed the presence of taste bud cells, while synaptophysin provided insights into taste bud and nerve bundle characteristics. CD34 and PDGFRα illuminated lingual stromal cells, and SOX9 and PDGFRα shed light on chondrocytes within the tongue's cartilage. CD20 mapped B-cell lymphocyte distribution in the lingual tonsil, while alpha smooth actin (α-SMA) exposed the intricate myofibroblast and smooth muscle network surrounding the lingual blood vessels and salivary glands. In conclusion, our comprehensive study advances our knowledge of the red-eared slider's tongue anatomy and physiology, addressing a significant research gap. These findings not only contribute to the field of turtle biology but also deepen our appreciation for the species' remarkable adaptations in their specific ecological niches.

Matrilineal phylogeny and habitat suitability of the endangered spotted pond turtle (Geoclemys hamiltonii; Testudines: Geoemydidae): a two-dimensional approach to forecasting future conservation consequences

The spotted pond turtle (Geoclemys hamiltonii) is a threatened and less explored species endemic to Bangladesh, India, Nepal, and Pakistan. To infer structural variation and matrilineal phylogenetic interpretation, the present research decoded the mitogenome of G. hamiltonii (16,509 bp) using next-generation sequencing technology. The mitogenome comprises 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and one AT-rich control region (CR) with similar strand symmetry in vertebrates. The ATG was identified as a start codon in most of the PCGs except Cytochrome oxidase subunit 1 (cox1), which started with the GTG codon. The non-coding CR of G. hamiltonii was determined to have a unique structure and variation in different domains and stem-loop secondary structure as compared with other Batagurinae species. The PCGs-based Bayesian phylogeny inferred strong monophyletic support for all Batagurinae species and confirmed the sister relationship of G. hamiltonii with Pangshura and Batagur taxa. We recommend generating more mitogenomic data for other Batagurinae species to confirm their population structure and evolutionary relationships. In addition, the present study aims to infer the habitat suitability and habitat quality of G. hamiltonii in its global distribution, both in the present and future climatic scenarios. We identify that only 58,542 km2 (7.16%) of the total range extent (817,341 km2) is suitable for this species, along with the fragmented habitats in both the eastern and western ranges. Comparative habitat quality assessment suggests the level of patch shape in the western range is higher (71.3%) compared to the eastern range. Our results suggest a massive decline of approximately 65.73% to 70.31% and 70.53% to 75.30% under ssp245 and ssp585 future scenarios, respectively, for the years between 2021-2040 and 2061-2080 compared with the current distribution. The present study indicates that proper conservation management requires greater attention to the causes and solutions to the fragmented distribution and safeguarding of this endangered species in the Indus, Ganges, and Brahmaputra (IGB) river basins.

A novel herpesvirus from a wild-caught Madagascar spider tortoise shows evidence of host-viral coevolution with a duplication event in Durocryptodira

Herpesviruses can be significant reptile pathogens. Herpesviral infection in a wild-caught, male spider tortoise (Pyxis arachnoides) under human care was detected during a routine wellness examination prior to transition between zoologic organizations. The tortoise had no clinical signs of illness. Oral swabs obtained during a physical examination as part of pre-shipment risk mitigation for infectious disease were submitted for consensus herpesvirus PCR assay and sequencing. Based on comparative sequence analysis, the novel herpesvirus identified is a member of the subfamily Alphaherpesvirinae. Studies of herpesviral phylogeny in chelonian species support branching patterns of turtle herpesviruses that closely mirror those of their hosts. The symmetry of these patterns is suggestive of close codivergence of turtle herpesviruses with their host species. The distribution of these viruses in both tortoises and emydids suggests a phylogenetic duplication event in the herpesviruses after host divergence of the Pleurodira and basal to the divergence of Americhelydia. Herpesviral infections have been documented to cause higher morbidity when introduced to aberrant host species, and significant consideration must be given to the presence of herpesviruses in the management of tortoise collections, particularly collections that include various species of testudines.

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.

Genome Evolution and the Future of Phylogenomics of Non-Avian Reptiles

Non-avian reptiles comprise a large proportion of amniote vertebrate diversity, with squamate reptiles-lizards and snakes-recently overtaking birds as the most species-rich tetrapod radiation. Despite displaying an extraordinary diversity of phenotypic and genomic traits, genomic resources in non-avian reptiles have accumulated more slowly than they have in mammals and birds, the remaining amniotes. Here we review the remarkable natural history of non-avian reptiles, with a focus on the physical traits, genomic characteristics, and sequence compositional patterns that comprise key axes of variation across amniotes. We argue that the high evolutionary diversity of non-avian reptiles can fuel a new generation of whole-genome phylogenomic analyses. A survey of phylogenetic investigations in non-avian reptiles shows that sequence capture-based approaches are the most commonly used, with studies of markers known as ultraconserved elements (UCEs) especially well represented. However, many other types of markers exist and are increasingly being mined from genome assemblies in silico, including some with greater information potential than UCEs for certain investigations. We discuss the importance of high-quality genomic resources and methods for bioinformatically extracting a range of marker sets from genome assemblies. Finally, we encourage herpetologists working in genomics, genetics, evolutionary biology, and other fields to work collectively towards building genomic resources for non-avian reptiles, especially squamates, that rival those already in place for mammals and birds. Overall, the development of this cross-amniote phylogenomic tree of life will contribute to illuminate interesting dimensions of biodiversity across non-avian reptiles and broader amniotes.

Cerebellar Inputs in the American Alligator (Alligator mississippiensis)

Crocodilians (alligators, crocodiles, and gharials) are the closet living relatives to birds and, as such, represent a key clade to understand the evolution of the avian brain. However, many aspects of crocodilian neurobiology remain unknown. In this paper, we address an important knowledge gap as there are no published studies of cerebellar connections in any crocodilian species. We used injections of retrograde tracers into the cerebellum of the American alligator (Alligator mississippiensis) to describe for the first time the origin of climbing and mossy fiber inputs. We found that inputs to the cerebellum in the American alligator are similar to those of other nonavian reptiles and birds. Retrograde labeled cells were found in the spinal cord, inferior olive, reticular formation, vestibular and cerebellar nuclei, as well as in nucleus ruber and surrounding tegmentum. Additionally, we found no retrogradely labeled cells in the anterior rhombencephalon which suggest that, like other nonavian reptiles, crocodilians may lack pontine nuclei. Similar to birds and other nonavian reptiles, we found inputs to the cerebellum from the pretectal nucleus lentiformis mesencephali. Additionally, we found retrogradely labeled neurons in two nuclei in the pretectum: the nucleus circularis and the interstitial nucleus of the posterior commissure. These pretectal projections have not been described in any other nonavian reptile to date, but they do resemble projections from the nucleus spiriformis medialis of birds. Our results show that many inputs to the cerebellum are highly conserved among sauropsids and that extensive pretectal inputs to the cerebellum are not exclusive to the avian brain. Finally, we suggest that the pontine nuclei of birds are an evolutionary novelty that may have evolved after the last common ancestor between birds and crocodilians, and may represent an intriguing case of convergent evolution with mammals.

Reference Genome of the Northwestern Pond Turtle, Actinemys marmorata

The northwestern pond turtle, Actinemys marmorata, and its recently recognized sister species, the southwestern pond turtle, A. pallida, are the sole aquatic testudines occurring over most of western North America and the only living representatives of the genus Actinemys. Although it historically ranged from Washington state through central California, USA, populations of the northwestern pond turtle have been in decline for decades and the species is afforded state-level protection across its range; it is currently being considered for protection under the US Endangered Species Act. Here, we report a new, chromosome-level assembly of A. marmorata as part of the California Conservation Genomics Project (CCGP). Consistent with the reference genome strategy of the CCGP, we used Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology to produce a de novo assembled genome. The assembly comprises 198 scaffolds spanning 2,319,339,408 base pairs, has a contig N50 of 75 Mb, a scaffold N50 of 146Mb, and BUSCO complete score of 96.7%, making it the most complete testudine assembly of the 24 species from 13 families that are currently available. In combination with the A. pallida reference genome that is currently under construction through the CCGP, the A. marmorata genome will be a powerful tool for documenting landscape genomic diversity, the basis of adaptations to salt tolerance and thermal capacity, and hybridization dynamics between these recently diverged species.

Genetic variation among sea turtle life stages and species suggests connectivity among ocean basins

Regional genetic differentiation of mitochondrial lineages occurs in migratory species with natal philopatry such as sea turtles. However, early juvenile dispersal represents a key opportunity for gene flow and colonization of new regions through founder events, making it an important yet under-studied life stage. To assess connectivity among sea turtle life stages and ocean basins, we sequenced mitochondrial DNA (mtDNA) fragments from 35 juveniles sampled in the Gulf of Mexico from the rarely observed dispersal stage across three species: green turtles (Chelonia mydas; n = 30), hawksbills (Eretmochelys imbricata; n = 3), and loggerheads (Caretta caretta; n = 2). We estimated green turtle rookery contributions using a many-to-many Bayesian mixed stock analysis that incorporated dispersal probabilities based on rookery size and transport via ocean currents. We assembled a gene tree including 709 distinct mtDNA control region haplotypes from the literature for all seven extant sea turtle species to assess gaps in life-stage data across ocean basins, as well as contextualize the lineages we sampled from dispersing juveniles. Our results indicate a high likelihood that green turtles sampled in the Gulf of Mexico originated from rookeries along the coast of Mexico, with smaller contributions from Costa Rica and Suriname. The gene tree analysis yielded species-level relationships consistent with those presented previously, while intra-species relationships between lineages and ocean basins differed, particularly within loggerhead and green turtle clades. Our results highlight the lack of genetic data from juvenile sea turtles, especially the early dispersal stage, and the potential for these data to answer broader questions of connectivity and diversification across species and lineages.

Chromosome-level genome assembly for the Aldabra giant tortoise enables insights into the genetic health of a threatened population

BACKGROUND

The Aldabra giant tortoise (Aldabrachelys gigantea) is one of only two giant tortoise species left in the world. The species is endemic to Aldabra Atoll in Seychelles and is listed as Vulnerable on the International Union for Conservation of Nature Red List (v2.3) due to its limited distribution and threats posed by climate change. Genomic resources for A. gigantea are lacking, hampering conservation efforts for both wild and ex situpopulations. A high-quality genome would also open avenues to investigate the genetic basis of the species' exceptionally long life span.

FINDINGS

We produced the first chromosome-level de novo genome assembly of A. gigantea using PacBio High-Fidelity sequencing and high-throughput chromosome conformation capture. We produced a 2.37-Gbp assembly with a scaffold N50 of 148.6 Mbp and a resolution into 26 chromosomes. RNA sequencing-assisted gene model prediction identified 23,953 protein-coding genes and 1.1 Gbp of repetitive sequences. Synteny analyses among turtle genomes revealed high levels of chromosomal collinearity even among distantly related taxa. To assess the utility of the high-quality assembly for species conservation, we performed a low-coverage resequencing of 30 individuals from wild populations and two zoo individuals. Our genome-wide population structure analyses detected genetic population structure in the wild and identified the most likely origin of the zoo-housed individuals. We further identified putatively deleterious mutations to be monitored.

CONCLUSIONS

We establish a high-quality chromosome-level reference genome for A. gigantea and one of the most complete turtle genomes available. We show that low-coverage whole-genome resequencing, for which alignment to the reference genome is a necessity, is a powerful tool to assess the population structure of the wild population and reveal the geographic origins of ex situ individuals relevant for genetic diversity management and rewilding efforts.

A systematic compendium of turtle mandibular anatomy using digital dissections of soft tissue and osteology

Turtles are a charismatic reptile group with a peculiar body plan, which most notably includes the shell. Anatomists have often focused descriptive efforts on the shell and other strongly derived body parts, such as the akinetic skull, or the cervical vertebrae. Other parts of turtle osteology, like the girdles, limbs, and mandibles, are documented with less rigor and detail. The mandible is the primary skeletal element involved in food acquisition and initial food processing of turtles, and its features are thus likely linked to feeding ecology. In addition, the mandible of turtles is composed of up to seven bones (sometimes fused to as little as three) and has thus anatomical complexity that may be insightful for systematic purposes and phylogenetic research. Despite apparent complexity and diversity to the mandible of turtles, this anatomical system has not been systematically studied, not even in search of characters that might improve phylogenetic resolution. Here, we describe the mandibular osteology for all major subclades of extant turtles with the help of digitally dissected 3D models derived from high-resolution computed tomography (μCT) scans of 70 extant species. We provide 31 fully segmented mandibles, as well as 3D models of the mandibular musculature, innervation, and arterial circulation of the cryptodire Dermatemys mawii. We synthesize observed variation into 51 morphological characters, which we optimize onto a molecular phylogeny. This analysis shows some mandibular characters to have high systematic value, whereas others are highly homoplastic and may underlie ecological influences or other factors invoking variation.

Phylogenomics and Diversification of the Schistosomatidae Based on Targeted Sequence Capture of Ultra-Conserved Elements

Schistosomatidae Stiles and Hassall 1898 is a medically significant family of digenetic trematodes (Trematoda: Digenea), members of which infect mammals or birds as definitive hosts and aquatic or amphibious gastropods as intermediate hosts. Currently, there are 17 named genera, for many of which evolutionary interrelationships remain unresolved. The lack of a resolved phylogeny has encumbered our understanding of schistosomatid evolution, specifically patterns of host-use and the role of host-switching in diversification. Here, we used targeted sequence capture of ultra-conserved elements (UCEs) from representatives of 13 of the 17 named genera and 11 undescribed lineages that are presumed to represent either novel genera or species to generate a phylogenomic dataset for the estimation of schistosomatid interrelationships. This study represents the largest phylogenetic effort within the Schistosomatidae in both the number of loci and breadth of taxon sampling. We present a near-comprehensive family-level phylogeny providing resolution to several clades of long-standing uncertainty within Schistosomatidae, including resolution for the placement of the North American mammalian schistosomes, implying a second separate capture of mammalian hosts. Additionally, we present evidence for the placement of Macrobilharzia at the base of the Schistosoma + Bivitellobilharzia radiation. Patterns of definitive and intermediate host use and a strong role for intermediate host-switching are discussed relative to schistosomatid diversification.

Scratching the Itch: Updated Perspectives on the Schistosomes Responsible for Swimmer's Itch around the World

Although most studies of digenetic trematodes of the family Schistosomatidae dwell on representatives causing human schistosomiasis, the majority of the 130 identified species of schistosomes infect birds or non-human mammals. The cercariae of many of these species can cause swimmer's itch when they penetrate human skin. Recent years have witnessed a dramatic increase in our understanding of schistosome diversity, now encompassing 17 genera with eight more lineages awaiting description. Collectively, schistosomes exploit 16 families of caenogastropod or heterobranch gastropod intermediate hosts. Basal lineages today are found in marine gastropods and birds, but subsequent diversification has largely taken place in freshwater, with some reversions to marine habitats. It seems increasingly likely that schistosomes have on two separate occasions colonized mammals. Swimmer's itch is a complex zoonotic disease manifested through several different routes of transmission involving a diversity of different host species. Swimmer's itch also exemplifies the value of adopting the One Health perspective in understanding disease transmission and abundance because the schistosomes involved have complex life cycles that interface with numerous species and abiotic components of their aquatic environments. Given the progress made in revealing their diversity and biology, and the wealth of questions posed by itch-causing schistosomes, they provide excellent models for implementation of long-term interdisciplinary studies focused on issues pertinent to disease ecology, the One Health paradigm, and the impacts of climate change, biological invasions and other environmental perturbations.

Matamatas Chelus spp. (Testudines, Chelidae) have a remarkable evolutionary history of sex chromosomes with a long-term stable XY microchromosome system

The genus Chelus, commonly known as Matamata is one of the most emblematic and remarkable species among the Neotropical chelids. It is an Amazonian species with an extensive distribution throughout Negro/Orinoco and Amazonas River basins. Currently, two species are formally recognized: Chelus orinocensis and Chelus fimbriata and although it is still classified as "Least Concern" in the IUCN, the Matamatas are very appreciated and illegally sold in the international pet trade. Regardless, little is known regarding many aspects of its natural history. Chromosomal features for Chelus, for instance, are meagre and practically restricted to the description of the diploid number (2n = 50) for Chelus fimbriata, and its sex determining strategies are yet to be fully investigated. Here, we examined the karyotype of Chelus fimbriata and the newly described Chelus orinocensis, applying an extensive conventional and molecular cytogenetic approach. This allowed us to identify a genetic sex determining mechanism with a micro XY sex chromosome system in both species, a system that was likely present in their most common recent ancestor Chelus colombiana. Furthermore, the XY system found in Chelus orinocensis and Chelus fimbriata, as seen in other chelid species, recruited several repeat motifs, possibly prior to the split of South America and Australasian lineages, indicating that such system indeed dates back to the earliest lineages of Chelid species.

A novel herpesvirus detected in 3 species of chelonians

Herpesviruses are found in free-living and captive chelonian populations, often in association with morbidity and mortality. To date, all known chelonian herpesviruses fall within the subfamily Alphaherpesvirinae. We detected a novel herpesvirus in 3 species of chelonians: a captive leopard tortoise (Stigmochelys pardalis) in western TX, USA; a steppe tortoise (Testudo [Agrionemys] horsfieldii) found near Fort Irwin, CA, USA; and 2 free-living, three-toed box turtles (Terrapene mexicana triunguis) found in Forest Park, St. Louis, MO. The leopard tortoise was coinfected with the tortoise intranuclear coccidian and had clinical signs of upper respiratory tract disease. The steppe tortoise had mucopurulent nasal discharge and lethargy. One of the three-toed box turtles had no clinical signs; the other was found dead with signs of trauma after being observed with blepharedema, tympanic membrane swelling, cervical edema, and other clinical signs several weeks prior to death. Generally, the branching order of the turtle herpesviruses mirrors the divergence patterns of their hosts, consistent with codivergence. Based on phylogenetic analysis, this novel herpesvirus clusters with a clade of viruses that infect emydid hosts and is likely of box turtle origin. Therefore, we suggest the name terrapene alphaherpesvirus 3 (TerAHV3) for the novel virus. This virus also has the ability to host-jump to tortoises, and previously documented herpesviral morbidity tends to be more common in aberrant hosts. The relationship between clinical signs and infection with TerAHV3 in these animals is unclear, and further investigation is merited.

A Genomic Perspective on the Evolutionary Diversification of Turtles

To examine phylogenetic heterogeneity in turtle evolution, we collected thousands of high-confidence single-copy orthologs from 19 genome assemblies representative of extant turtle diversity and estimated a phylogeny with multispecies coalescent and concatenated partitioned methods. We also collected next-generation sequences from 26 turtle species and assembled millions of biallelic markers to reconstruct phylogenies based on annotated regions from the western painted turtle (Chrysemys picta bellii) genome (coding regions, introns, untranslated regions, intergenic, and others). We then measured gene tree-species tree discordance, as well as gene and site heterogeneity at each node in the inferred trees, and tested for temporal patterns in phylogenomic conflict across turtle evolution. We found strong and consistent support for all bifurcations in the inferred turtle species phylogenies. However, a number of genes, sites, and genomic features supported alternate relationships between turtle taxa. Our results suggest that gene tree-species tree discordance in these datasets is likely driven by population-level processes such as incomplete lineage sorting. We found very little effect of substitutional saturation on species tree topologies, and no clear phylogenetic patterns in codon usage bias and compositional heterogeneity. There was no correlation between gene and site concordance, node age, and DNA substitution rate across most annotated genomic regions. Our study demonstrates that heterogeneity is to be expected even in well resolved clades such as turtles, and that future phylogenomic studies should aim to sample as much of the genome as possible in order to obtain accurate phylogenies for assessing conservation priorities in turtles.

Divergence and hybridization in sea turtles: Inferences from genome data show evidence of ancient gene flow between species

Reconstructing past events of hybridization and population size changes are required to understand speciation mechanisms and current patterns of genetic diversity, and ultimately contribute to species' conservation. Sea turtles are ancient species currently facing anthropogenic threats including climate change, fisheries, and illegal hunting. Five of the seven extant sea turtle species are known to currently hybridize, especially along the Brazilian coast where some populations can have ~32%-42% of hybrids. Although frequently observed today, it is not clear what role hybridization plays in the evolutionary diversification of this group of reptiles. In this study, we generated whole genome resequencing data of the five globally distributed sea turtle species to estimate a calibrated phylogeny and the population size dynamics, and to understand the role of hybridization in shaping the genomes of these ancient species. Our results reveal discordant species divergence dates between mitochondrial and nuclear genomes, with a high frequency of conflicting trees throughout the nuclear genome suggesting that some sea turtle species frequently hybridized in the past. The reconstruction of the species' demography showed a general decline in effective population sizes with no signs of recovery, except for the leatherback sea turtle. Furthermore, we discuss the influence of reference bias in our estimates. We show long-lasting ancestral gene flow events within Chelonioidea that continued for millions of years after initial divergence. Speciation with gene flow is a common pattern in marine species, and it raises questions whether current hybridization events should be considered as a part of these species' evolutionary history or a conservation issue.

Sex Chromosomes and Master Sex-Determining Genes in Turtles and Other Reptiles

Among tetrapods, the well differentiated heteromorphic sex chromosomes of birds and mammals have been highly investigated and their master sex-determining (MSD) gene, Dmrt1 and SRY, respectively, have been identified. The homomorphic sex chromosomes of reptiles have been the least studied, but the gap with birds and mammals has begun to fill. This review describes our current knowledge of reptilian sex chromosomes at the cytogenetic and molecular level. Most of it arose recently from various studies comparing male to female gene content. This includes restriction site-associated DNA sequencing (RAD-Seq) experiments in several male and female samples, RNA sequencing and identification of Z- or X-linked genes by male/female comparative transcriptome coverage, and male/female transcriptomic or transcriptome/genome substraction approaches allowing the identification of Y- or W-linked transcripts. A few putative master sex-determining (MSD) genes have been proposed, but none has been demonstrated yet. Lastly, future directions in the field of reptilian sex chromosomes and their MSD gene studies are considered.

A brief review of vertebrate sex evolution with a pledge for integrative research: towards 'sexomics'

Triggers and biological processes controlling male or female gonadal differentiation vary in vertebrates, with sex determination (SD) governed by environmental factors or simple to complex genetic mechanisms that evolved repeatedly and independently in various groups. Here, we review sex evolution across major clades of vertebrates with information on SD, sexual development and reproductive modes. We offer an up-to-date review of divergence times, species diversity, genomic resources, genome size, occurrence and nature of polyploids, SD systems, sex chromosomes, SD genes, dosage compensation and sex-biased gene expression. Advances in sequencing technologies now enable us to study the evolution of SD at broader evolutionary scales, and we now hope to pursue a sexomics integrative research initiative across vertebrates. The vertebrate sexome comprises interdisciplinary and integrated information on sexual differentiation, development and reproduction at all biological levels, from genomes, transcriptomes and proteomes, to the organs involved in sexual and sex-specific processes, including gonads, secondary sex organs and those with transcriptional sex-bias. The sexome also includes ontogenetic and behavioural aspects of sexual differentiation, including malfunction and impairment of SD, sexual differentiation and fertility. Starting from data generated by high-throughput approaches, we encourage others to contribute expertise to building understanding of the sexomes of many key vertebrate species. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)'.

Turtles and Tortoises Are in Trouble

Turtles and tortoises (chelonians) have been integral components of global ecosystems for about 220 million years and have played important roles in human culture for at least 400,000 years. The chelonian shell is a remarkable evolutionary adaptation, facilitating success in terrestrial, freshwater and marine ecosystems. Today, more than half of the 360 living species and 482 total taxa (species and subspecies combined) are threatened with extinction. This places chelonians among the groups with the highest extinction risk of any sizeable vertebrate group. Turtle populations are declining rapidly due to habitat loss, consumption by humans for food and traditional medicines and collection for the international pet trade. Many taxa could become extinct in this century. Here, we examine survival threats to turtles and tortoises and discuss the interventions that will be needed to prevent widespread extinction in this group in coming decades.

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

Living turtles are characterized by extraordinarily low species diversity given their age. The clade's extensive fossil record indicates that climate and biogeography may have played important roles in determining their diversity. We investigated this hypothesis by collecting a molecular dataset for 591 individual turtles that, together, represent 80% of all turtle species, including representatives of all families and 98% of genera, and used it to jointly estimate phylogeny and divergence times. We found that the turtle tree is characterized by relatively constant diversification (speciation minus extinction) punctuated by a single threefold increase. We also found that this shift is temporally and geographically associated with newly emerged continental margins that appeared during the Eocene-Oligocene transition about 30 million years before present. In apparent contrast, the fossil record from this time period contains evidence for a major, but regional, extinction event. These seemingly discordant findings appear to be driven by a common global process: global cooling and drying at the time of the Eocene-Oligocene transition. This climatic shift led to aridification that drove extinctions in important fossil-bearing areas, while simultaneously exposing new continental margin habitat that subsequently allowed for a burst of speciation associated with these newly exploitable ecological opportunities.

Evolutionary history of mental glands in turtles reveals a single origin in an aquatic ancestor and recurrent losses independent of macrohabitat

Despite the relevance of chemical communication in vertebrates, comparative examinations of macroevolutionary trends in chemical signaling systems are scarce. Many turtle and tortoise species are reliant on chemical signals to communicate in aquatic and terrestrial macrohabitats, and many of these species possess specialized integumentary organs, termed mental glands (MGs), involved in the production of chemosignals. We inferred the evolutionary history of MGs and tested the impact of macrohabitat on their evolution. Inference of ancestral states along a time-calibrated phylogeny revealed a single origin in the ancestor of the subclade Testudinoidea. Thus, MGs represent homologous structures in all descending lineages. We also inferred multiple independent losses of MGs in both terrestrial and aquatic clades. Although MGs first appeared in an aquatic turtle (the testudinoid ancestor), macrohabitat seems to have had little effect on MG presence or absence in descendants. Instead, we find clade-specific evolutionary trends, with some clades showing increased gland size and morphological complexity, whereas others exhibiting reduction or MG loss. In sister clades inhabiting similar ecological niches, contrasting patterns (loss vs. maintenance) may occur. We conclude that the multiple losses of MGs in turtle clades have not been influenced by macrohabitat and that other factors have affected MG evolution.

Draft Genome of the Common Snapping Turtle, Chelydra serpentina, a Model for Phenotypic Plasticity in Reptiles

Turtles are iconic reptiles that inhabit a range of ecosystems from oceans to deserts and climates from the tropics to northern temperate regions. Yet, we have little understanding of the genetic adaptations that allow turtles to survive and reproduce in such diverse environments. Common snapping turtles, Chelydra serpentina, are an ideal model species for studying adaptation to climate because they are widely distributed from tropical to northern temperate zones in North America. They are also easy to maintain and breed in captivity and produce large clutch sizes, which makes them amenable to quantitative genetic and molecular genetic studies of traits like temperature-dependent sex determination. We therefore established a captive breeding colony and sequenced DNA from one female using both short and long reads. After trimming and filtering, we had 209.51Gb of Illumina reads, 25.72Gb of PacBio reads, and 21.72 Gb of Nanopore reads. The assembled genome was 2.258 Gb in size and had 13,224 scaffolds with an N50 of 5.59Mb. The longest scaffold was 27.24Mb. BUSCO analysis revealed 97.4% of core vertebrate genes in the genome. We identified 3.27 million SNPs in the reference turtle, which indicates a relatively high level of individual heterozygosity. We assembled the transcriptome using RNA-Seq data and used gene prediction software to produce 22,812 models of protein coding genes. The quality and contiguity of the snapping turtle genome is similar to or better than most published reptile genomes. The genome and genetic variants identified here provide a foundation for future studies of adaptation to climate.

The role of selection in the evolution of marine turtles mitogenomes

Sea turtles are the only extant chelonian representatives that inhabit the marine environment. One key to successful colonization of this habitat is the adaptation to different energetic demands. Such energetic requirement is intrinsically related to the mitochondrial ability to generate energy through oxidative phosphorylation (OXPHOS) process. Here, we estimated Testudines phylogenetic relationships from 90 complete chelonian mitochondrial genomes and tested the adaptive evolution of 13 mitochondrial protein-coding genes of sea turtles to determine how natural selection shaped mitochondrial genes of the Chelonioidea clade. Complete mitogenomes showed strong support and resolution, differing at the position of the Chelonioidea clade in comparison to the turtle phylogeny based on nuclear genomic data. Codon models retrieved a relatively increased dN/dS (ω) on three OXPHOS genes for sea turtle lineages. Also, we found evidence of positive selection on at least three codon positions, encoded by NADH dehydrogenase genes (ND4 and ND5). The accelerated evolutionary rates found for sea turtles on COX2, ND1 and CYTB and the molecular footprints of positive selection found on ND4 and ND5 genes may be related to mitochondrial molecular adaptation to stress likely resulted from a more active lifestyle in sea turtles. Our study provides insight into the adaptive evolution of the mtDNA genome in sea turtles and its implications for the molecular mechanism of oxidative phosphorylation.

The Amazonian Red Side-Necked Turtle Rhinemys rufipes (Spix, 1824) (Testudines, Chelidae) Has a GSD Sex-Determining Mechanism with an Ancient XY Sex Microchromosome System

The Amazonian red side-necked turtle Rhynemis rufipes is an endemic Amazonian Chelidae species that occurs in small streams throughout Colombia and Brazil river basins. Little is known about various biological aspects of this species, including its sex determination strategies. Among chelids, the greatest karyotype diversity is found in the Neotropical species, with several 2n configurations, including cases of triploidy. Here, we investigate the karyotype of Rhinemys rufipes by applying combined conventional and molecular cytogenetic procedures. This allowed us to discover a genetic sex-determining mechanism that shares an ancestral micro XY sex chromosome system. This ancient micro XY system recruited distinct repeat motifs before it diverged from several South America and Australasian species. We propose that such a system dates back to the earliest lineages of the chelid species before the split of South America and Australasian lineages.

Ancient tropical extinctions at high latitudes contributed to the latitudinal diversity gradient

Global biodiversity currently peaks at the equator and decreases toward the poles. Growing fossil evidence suggest this hump-shaped latitudinal diversity gradient (LDG) has not been persistent through time, with similar diversity across latitudes flattening out the LDG during past greenhouse periods. However, when and how diversity declined at high latitudes to generate the modern LDG remains an open question. Although diversity-loss scenarios have been proposed, they remain mostly undemonstrated. We outline the "asymmetric gradient of extinction and dispersal" framework that contextualizes previous ideas behind the LDG under a time-variable scenario. Using phylogenies and fossils of Testudines, Crocodilia, and Lepidosauria, we find that the hump-shaped LDG could be explained by (1) disproportionate extinctions of high-latitude tropical-adapted clades when climate transitioned from greenhouse to icehouse, and (2) equator-ward biotic dispersals tracking their climatic preferences when tropical biomes became restricted to the equator. Conversely, equivalent diversification rates across latitudes can account for the formation of an ancient flat LDG. The inclusion of fossils in macroevolutionary studies allows revealing time-dependent extinction rates hardly detectable from phylogenies only. This study underscores that the prevailing evolutionary processes generating the LDG during greenhouses differed from those operating during icehouses.

Distribution of GFAP in Squamata: Extended Immunonegative Areas, Astrocytes, High Diversity, and Their Bearing on Evolution

Squamata is one of the richest and most diverse extant groups. The present study investigates the glial fibrillary acidic protein (GFAP)-immunopositive elements of five lizard and three snake species; each represents a different family. The study continues our former studies on bird, turtle, and caiman brains. Although several studies have been published on lizards, they usually only investigated one species. Almost no data are available on snakes. The animals were transcardially perfused. Immunoperoxidase reactions were performed with a mouse monoclonal anti-GFAP (Novocastra). The original radial ependymoglia is enmeshed by secondary, non-radial processes almost beyond recognition in several brain areas like in other reptiles. Astrocytes occur but only as complementary elements like in caiman but unlike in turtles, where astrocytes are absent. In most species, extended areas are free of GFAP—a meaningful difference from other reptiles. The predominance of astrocytes and the presence of areas free of GFAP immunopositivity are characteristic of birds and mammals; therefore, they must be apomorphic features of Squamata, which appeared independently from the evolution of avian glia. However, these features show a high diversity; in some lizards, they are even absent. There was no principal difference between the glial structures of snakes and lizards. In conclusion, the glial structure of Squamata seems to be the most apomorphic one among reptiles. The high diversity suggests that its evolution is still intense. The comparison of identical brain areas with different GFAP contents in different species may promote understanding the role of GFAP in neuronal networks. Our findings are in accordance with the supposal based on our previous studies that the GFAP-free areas expand during evolution.

Transcriptome-based target-enrichment baits for stony corals (Cnidaria: Anthozoa: Scleractinia)

Despite the ecological and economic significance of stony corals (Scleractinia), a robust understanding of their phylogeny remains elusive due to patchy taxonomic and genetic sampling, as well as the limited availability of informative markers. To increase the number of genetic loci available for phylogenomic analyses in Scleractinia, we designed 15,919 DNA enrichment baits targeting 605 orthogroups (mean 565 ± SD 366 bp) over 1,139 exon regions. A further 236 and 62 barcoding baits were designed for COI and histone H3 genes respectively for quality and contamination checks. Hybrid capture using these baits was performed on 18 coral species spanning the presently understood scleractinian phylogeny, with two corallimorpharians as outgroup. On average, 74% of all loci targeted were successfully captured for each species. Barcoding baits were matched unambiguously to their respective samples and revealed low levels of cross-contamination in accordance with expectation. We put the data through a series of stringent filtering steps to ensure only scleractinian and phylogenetically informative loci were retained, and the final probe set comprised 13,479 baits, targeting 452 loci (mean 531 ± SD 307 bp) across 865 exon regions. Maximum likelihood, Bayesian and species tree analyses recovered maximally supported, topologically congruent trees consistent with previous phylogenomic reconstructions. The phylogenomic method presented here allows for consistent capture of orthologous loci among divergent coral taxa, facilitating the pooling of data from different studies and increasing the phylogenetic sampling of scleractinians in the future.

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

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

Genomic analyses reveal two species of the matamata (Testudines: Chelidae: Chelus spp.) and clarify their phylogeography

The matamata is one of the most charismatic turtles on earth, widely distributed in northern South America. Debates have occurred over whether or not there should be two subspecies or species recognized due to its geographic variation in morphology. Even though the matamata is universally known, its natural history, conservation status and biogeography are largely unexplored. In this study we examined the phylogeographic differentiation of the matamata based on three mitochondrial DNA fragments (2168 bp of the control region, cytochrome oxidase subunit I, and the cytochrome b gene), one nuclear genomic DNA fragment (1068 bp of the R35 intron) and 1661 Single Nucleotide Polymorphisms (SNPs). Our molecular and morphological analyses revealed the existence of two distinct, genetically deeply divergent evolutionary lineages of matamatas that separated in the late Miocene (approximately 12.7 million years ago), corresponding well to the time when the Orinoco Basin was established. As a result of our analyses, we describe the genetically and morphologically highly distinct matamata from the Orinoco and Río Negro Basins and the Essequibo drainage as a species new to science (Chelus orinocensis sp. nov.). Chelus fimbriata sensu stricto is distributed in the Amazon Basin and the Mahury drainage. Additionally, the analyses revealed that each species displays phylogeographic differentiation. For C. orinocensis, there is moderate mitochondrial differentiation between the Orinoco and the Río Negro. For C. fimbriata, there is more pronounced differentiation matching different river systems. One mitochondrial clade was identified from the Amazon, Ucayali, and Mahury Rivers, and another one from the Madeira and Jaci Paraná Rivers. The C. orinocensis in the Essequibo and Branco Rivers have haplotypes that constitute a third clade clustering with C. fimbriata. Phylogenetic analyses of the R35 intron and SNP data link the matamatas from the Essequibo and Branco with the new species, suggesting past gene flow and old mitochondrial introgression. Chelus orinocensis is collected for the pet trade in Colombia and Venezuela. However, neither the extent of the harvest nor its impact are known. Hence, it is crucial to gather more information and to assess its exploitation throughout its distribution range to obtain a better understanding of its conservation status and to design appropriate conservation and management procedures. RESUMEN: La matamata es una de las tortugas más carismáticas del mundo, ampliamente distribuida en el norte de Sudamérica. Debido a su variación morfológica geográfica, se debate sobre el reconocimiento de dos subespecies o especies. A pesar de que la matamata es universalmente conocida, su historia natural, estado de conservación y biogeografía han sido muy poco estudiados. En este estudio examinamos la diferenciación filogeográfica de las matamatas en base ​​a tres fragmentos de ADN mitocondrial (2168 pb de la región de control, la subunidad I del citocromo oxidasa y el gen del citocromo b), un fragmento de ADN genómico nuclear (1068 pb del intrón R35) y 1661 polimorfismos de nucleótido único (SNPs). Nuestros análisis moleculares y morfológicos revelaron la existencia de dos linajes evolutivos distintos de matamatas, genéticamente divergentes que se separaron en el Mioceno tardio (hace aproximadamente 12.7 millones de años), correspondiendo al tiempo en que se estableció la cuenca del Orinoco. Como resultado de nuestros análisis, describimos las genéticamente y morfológicamente distintas matamatas de las cuencas del Orinoco, Río Negro y Essequibo como una especie nueva para la ciencia (Chelus orinocensis sp. nov.). Chelus fimbriata sensu stricto se distribuye en la cuenca del Amazonas y en el drenaje del Mahury. Adicionalmente, los análisis revelaron que cada especie muestra diferenciación filogeográfica. Para C. orinocensis, hay una moderada diferenciación mitocondrial entre el Orinoco y el Río Negro. Para C. fimbriata, hay una diferenciación más pronunciada, concordando con los diferentes sistemas fluviales. Se identificó un clado de los ríos Amazonas, Ucayali y Mahury y otro de los ríos Madeira y Jaci Paraná. Las C. orinocensis de los ríos Essequibo y Branco tienen haplotipos que constituyen un tercer clado que se agrupa con C. fimbriata. Los análisis filogenéticos del intrón R35 y los datos de SNP asocian las matamatas de Essequibo y Branco con la nueva especie, sugiriendo flujo de genes pasado ​​e introgresión mitocondrial antigua. Chelus orinocensis se colecta para el comercio de mascotas en Colombia y Venezuela. Sin embargo, ni se conoce el alcance de las colectas ni su impacto. Por lo tanto, es crucial recopilar más información y evaluar su explotación en todo su rango de distribución, comprender mejor su estado de conservación y para diseñar acciones apropiadas de conservación y manejo.

Zebrin Expression in the Cerebellum of Two Crocodilian Species

While in birds and mammals the cerebellum is a highly convoluted structure that consists of numerous transverse lobules, in most amphibians and reptiles it consists of only a single unfolded sheet. Orthogonal to the lobules, the cerebellum is comprised of sagittal zones that are revealed in the pattern of afferent inputs, the projection patterns of Purkinje cells, and Purkinje cell response properties, among other features. The expression of several molecular markers, such as aldolase C, is also parasagittally organized. Aldolase C, also known as zebrin II (ZII), is a glycolytic enzyme expressed in the cerebellar Purkinje cells of the vertebrate cerebellum. In birds, mammals, and some lizards (Ctenophoresspp.), ZII is expressed in a heterogenous fashion of alternating sagittal bands of high (ZII+) and low (ZII-) expression Purkinje cells. In contrast, turtles and snakes express ZII homogenously (ZII+) in their cerebella, but the pattern in crocodilians is unknown. Here, we examined the expression of ZII in two crocodilian species (Crocodylus niloticus and Alligator mississippiensis) to help determine the evolutionary origin of striped ZII expression in vertebrates. We expected crocodilians to express ZII in a striped (ZII+/ZII-) manner because of their close phylogenetic relationship to birds and their larger and more folded cerebellum compared to that of snakes and turtles. Contrary to our prediction, all Purkinje cells in the crocodilian cerebellum had a generally homogenous expression of ZII (ZII+) rather than clear ZII+/- stripes. Our results suggest that either ZII stripes were lost in three groups (snakes, turtles, and crocodilians) or ZII stripes evolved independently three times (lizards, birds, and mammals).

Using the Fossil Record to Evaluate Timetree Timescales

The fossil and geologic records provide the primary data used to established absolute timescales for timetrees. For the paleontological evaluation of proposed timetree timescales, and for node-based methods for constructing timetrees, the fossil record is used to bracket divergence times. Minimum brackets (minimum ages) can be established robustly using well-dated fossils that can be reliably assigned to lineages based on positive morphological evidence. Maximum brackets are much harder to establish, largely because it is difficult to establish definitive evidence that the absence of a taxon in the fossil record is real and not just due to the incompleteness of the fossil and rock records. Five primary methods have been developed to estimate maximum age brackets, each of which is discussed. The fact that the fossilization potential of a group typically decreases the closer one approaches its time of origin increases the challenge of estimating maximum age brackets. Additional complications arise: 1) because fossil data actually bracket the time of origin of the first relevant fossilizable morphology (apomorphy), not the divergence time itself; 2) due to the phylogenetic uncertainty in the placement of fossils; 3) because of idiosyncratic temporal and geographic gaps in the rock and fossil records; and 4) if the preservation potential of a group changed significantly during its history. In contrast, uncertainties in the absolute ages of fossils are typically relatively unimportant, even though the vast majority of fossil cannot be dated directly. These issues and relevant quantitative methods are reviewed, and their relative magnitudes assessed, which typically correlate with the age of the group, its geographic range, and species richness.

Smooth Muscle in Cardiac Chambers is Common in Turtles and Extensive in the Emydid Turtle, Trachemys scripta

A prominent layer of smooth muscle lining the luminal side of the atria of freshwater turtles (Emydidae) was described more than a century ago. We recently demonstrated that this smooth muscle provides a previously unrecognized mechanism to change cardiac output in the emydid red-eared slider (Trachemys scripta) that possibly contributes to their tremendous diving capacity. The purpose of the present immunohistochemical study was firstly to screen major groups of vertebrates for the presence of cardiac smooth muscle. Secondly, we investigated the phylogenetic distribution of cardiac smooth muscle within the turtle order (Testudines), including terrestrial and aquatic species. Atrial smooth muscle was not detected in a range of vertebrates, including Xenopus laevis, Alligator mississippiensis, and Caiman crocodilus, all of which have pronounced diving capacities. However, we confirmed earlier reports that traces of smooth muscle are found in human atrial tissue. Only within the turtles (eight species) was there substantial amounts of nonvascular smooth muscle in the heart. This amount was greatest in the atria, while the amount in proportion to cardiac muscle was greater in the sinus venosus than in other chambers. T. scripta had more smooth muscle in the sinus venosus and atria than the other turtles. In some specimens, there was some smooth muscle in the ventricle and the pulmonary vein. Our study demonstrates that cardiac smooth muscle likely appeared early in turtle evolution and has become extensive within the Emydidae family, possibly in association with diving. Across other tetrapod clades, cardiac smooth muscle might not associate with diving. Anat Rec, 303:1327-1336, 2020. © 2019 The Authors. The Anatomical Record published by Wiley Periodicals, Inc. on behalf of American Association for Anatomy.

Different mechanisms of serum complement activation in the plasma of common (Chelydra serpentina) and alligator (Macrochelys temminckii) snapping turtles

Reptiles are declining worldwide yet our understanding of their immune function lags far behind other taxa. The innate immune system is the primary mode of defense in reptiles, and the serum complement cascade is its major component. We assessed serum complement activity of plasma in two closely related aquatic turtle species, the common snapping turtle (CST; Chelydra serpentina) and alligator snapping turtle (AST; Macrochelys temminckii). We used a sheep red blood cell (SRBC) hemolysis assay to assess serum complement activity. Although the antibacterial activities of the plasma of these turtle species are similar, the hemolytic activity was much stronger in CST than AST. Treatment with inhibitors of the serum complement cascade indicated differences in the mechanisms of complement activation between the turtle species. We subjected plasma from both turtle species to mannan affinity chromatography and analyzed the eluate with SDS-PAGE, which revealed that plasma from the CSTs contained only small amounts of one C-type lectin protein while the AST plasma contained high concentrations of two C-type lectins (31.0 and 35.9 kDa). Edman degradation analyses confirmed that the two AST proteins contained identical N-terminal sequences. Thus, the CST appears to rely more heavily on the alternative mechanism of serum complement activation, while the AST appears to rely more on the lectin-mediated pathway, which is a pattern recognition response to prokaryotes not activated by the SRBCs. These results are unique in that the use of serum complement pathways are generally assumed to be conserved within clades.

Trabecular architecture in the humeral metaphyses of non-avian reptiles (Crocodylia, Squamata and Testudines): Lifestyle, allometry and phylogeny

The lifestyle of extinct tetrapods is often difficult to assess when clear morphological adaptations such as swimming paddles are absent. According to the hypothesis of bone functional adaptation, the architecture of trabecular bone adapts sensitively to physiological loadings. Previous studies have already shown a clear relation between trabecular architecture and locomotor behavior, mainly in mammals and birds. However, a link between trabecular architecture and lifestyle has rarely been examined. Here, we analyzed trabecular architecture of different clades of reptiles characterized by a wide range of lifestyles (aquatic, amphibious, generalist terrestrial, fossorial, and climbing). Humeri of squamates, turtles, and crocodylians have been scanned with microcomputed tomography. We selected spherical volumes of interest centered in the proximal metaphyses and measured trabecular spacing, thickness and number, degree of anisotropy, average branch length, bone volume fraction, bone surface density, and connectivity density. Only bone volume fraction showed a significant phylogenetic signal and its significant difference between squamates and other reptiles could be linked to their physiologies. We found negative allometric relationships for trabecular thickness and spacing, positive allometries for connectivity density and trabecular number and no dependence with size for degree of anisotropy and bone volume fraction. The different lifestyles are well separated in the morphological space using linear discriminant analyses, but a cross-validation procedure indicated a limited predictive ability of the model. The trabecular bone anisotropy has shown a gradient in turtles and in squamates: higher values in amphibious than terrestrial taxa. These allometric scalings, previously emphasized in mammals and birds, seem to be valid for all amniotes. Discriminant analysis has offered, to some extent, a distinction of lifestyles, which however remains difficult to strictly discriminate. Trabecular architecture seems to be a promising tool to infer lifestyle of extinct tetrapods, especially those involved in the terrestrialization.

Mitogenomics of historical type specimens of Australasian turtles: clarification of taxonomic confusion and old mitochondrial introgression

Diagnosability is central to taxonomy as are type specimens which define taxa. New advances in technologies and the discovery of new informative traits must be matched with previous taxonomic decisions based on name-bearing type specimens. Consequently, the challenge of sequencing highly degraded DNA from historical types becomes an inevitability to resolve the very many taxonomic issues arising from, by modern standards, poor historical species descriptions leading to difficulties to assign names to genetic clusters identified from fresh material. Here we apply high-throughput parallel sequencing and sequence baiting to reconstruct the mitogenomes from 18 type specimens of Australasian side-necked turtles (Chelidae). We resolve a number of important issues that have confused the taxonomy of this family, and analyse the mitogenomes of the types and those of fresh material to improve our understanding of the phylogenetic relationships of this morphologically conservative group. Together with previously published nuclear genomic data, our study provides evidence for multiple old mitochondrial introgressions.

Performance in three shell functions predicts the phenotypic distribution of hard-shelled turtles

Adaptive landscapes have served as fruitful guides to evolutionary research for nearly a century. Current methods guided by landscape frameworks mostly utilize evolutionary modeling (e.g., fitting data to Ornstein-Uhlenbeck models) to make inferences about adaptive peaks. Recent alternative methods utilize known relationships between phenotypes and functional performance to derive information about adaptive landscapes; this information can then help explain the distribution of species in phenotypic space and help infer the relative importance of various functions for guiding diversification. Here, data on performance for three turtle shell functions-strength, hydrodynamic efficiency, and self-righting ability-are used to develop a set of predicted performance optima in shell shape space. The distribution of performance optima shows significant similarity to the distribution of existing turtle species and helps explain the absence of shells in otherwise anomalously empty regions of morphospace. The method outperforms a modeling-based approach in inferring the location of reasonable adaptive peaks and in explaining the shape of the phenotypic distributions of turtle shells. Performance surface-based methods allow researchers to more directly connect functional performance with macroevolutionary diversification, and to explain the distribution of species (including presences and absences) across phenotypic space.

Variation Across Mitochondrial Gene Trees Provides Evidence for Systematic Error: How Much Gene Tree Variation Is Biological?

The use of large genomic data sets in phylogenetics has highlighted extensive topological variation across genes. Much of this discordance is assumed to result from biological processes. However, variation among gene trees can also be a consequence of systematic error driven by poor model fit, and the relative importance of biological vs. methodological factors in explaining gene tree variation is a major unresolved question. Using mitochondrial genomes to control for biological causes of gene tree variation, we estimate the extent of gene tree discordance driven by systematic error and employ posterior prediction to highlight the role of model fit in producing this discordance. We find that the amount of discordance among mitochondrial gene trees is similar to the amount of discordance found in other studies that assume only biological causes of variation. This similarity suggests that the role of systematic error in generating gene tree variation is underappreciated and critical evaluation of fit between assumed models and the data used for inference is important for the resolution of unresolved phylogenetic questions.

Antibacterial activities of plasma from the common (Chelydra serpentina) and alligator snapping turtle (Macrochelys temminckii)

Innate immunity provides a fast-acting and nonspecific defense against microbial infection, and appears to have particular importance in the immune response of ectothermic vertebrates. Chelonians are a globally distributed and diverse group, yet little is known about their basic immune function. The chelonian family Chelydridae is made up of two genera (Chelydra and Macrochelys), represented in our study by the widespread common snapping turtle ( Chelydra serpentina; CST) and the southeast USA endemic alligator snapping turtle ( Macrochelys temminckii; AST). Our goal was to quantify the innate immune response of the family Chelydridae, using the antibacterial activity of plasma as a measure of immune function. Our results show that the plasma of both species has strong antibacterial properties, but CST plasma kills a higher percentage of bacteria than AST plasma. In addition, while both species showed the highest antibacterial activity at 25 to 30°C, CST plasma retained its antibacterial properties at lower and higher temperatures than AST plasma. Our results indicate that, like many ectotherms, Chelydridae have a relatively strong innate immune response. The stronger, more robust immune response of CSTs compared with ASTs is likely correlated to the differences in geographic ranges but may also have implications for each species' tolerance to anthropogenic habitat degradation and global climate change.

Phylogeny, biogeography and diversification patterns of side-necked turtles (Testudines: Pleurodira)

Pleurodires or side-necked turtles are today restricted to freshwater environments of South America, Africa-Madagascar and Australia, but in the past they were distributed much more broadly, being found also on Eurasia, India and North America, and marine environments. Two hypotheses were proposed to explain this distribution; in the first, vicariance would have shaped the current geographical distribution and, in the second, extinctions constrained a previously widespread distribution. Here, we aim to reconstruct pleurodiran biogeographic history and diversification patterns based on a new phylogenetic hypothesis recovered from the analysis of the largest morphological dataset yet compiled for the lineage, testing which biogeographical process prevailed during its evolutionary history. The resulting topology generally agrees with previous hypotheses of the group and shows that most diversification shifts were related to the exploration of new niches, e.g. littoral or marine radiations. In addition, as other turtles, pleurodires do not seem to have been much affected by either the Cretaceous-Palaeogene or the Eocene-Oligocene mass extinctions. The biogeographic analyses highlight the predominance of both anagenetic and cladogenetic dispersal events and support the importance of transoceanic dispersals as a more common driver of area changes than previously thought, agreeing with previous studies with other non-turtle lineages.