Establishment and characterization of turtle liver organoids provides a potential model to decode their unique adaptations

Painted turtles are remarkable for their freeze tolerance and supercooling ability along with their associated resilience to hypoxia/anoxia and oxidative stress, rendering them an ideal biomedical model for hypoxia-induced injuries (including strokes), tissue cooling during surgeries, and organ cryopreservation. Yet, such research is hindered by their seasonal reproduction and slow maturation. Here we developed and characterized adult stem cell-derived turtle liver organoids (3D self-assembled in vitro structures) from painted, snapping, and spiny softshell turtles spanning ~175My of evolution, with a subset cryopreserved. This development is, to the best of our knowledge, a first for this vertebrate Order, and complements the only other non-avian reptile organoids from snake venom glands. Preliminary characterization, including morphological, transcriptomic, and proteomic analyses, revealed organoids enriched in cholangiocytes. Deriving organoids from distant turtles and life stages demonstrates that our techniques are broadly applicable to chelonians, permitting the development of functional genomic tools currently lacking in herpetological research. Such platform could potentially support studies including genome-to-phenome mapping, gene function, genome architecture, and adaptive responses to climate change, with implications for ecological, evolutionary, and biomedical research.

Transcriptomic thermal plasticity underlying gonadal development in a turtle with ZZ/ZW sex chromosomes despite canalized genotypic sex determination

Understanding genome-wide responses to environmental conditions during embryogenesis is essential for discerning the evolution of developmental plasticity and canalization, two processes generating phenotypic variation targeted by natural selection. Here, we present the first comparative trajectory analysis of matched transcriptomic developmental time series from two reptiles incubated under identical conditions, a turtle with a ZZ/ZW system of genotypic sex determination (GSD), Apalone spinifera, and a turtle with temperature-dependent sex determination (TSD), Chrysemys picta. Results from our genome-wide, hypervariate gene expression analysis of sexed embryos across five developmental stages revealed that substantial transcriptional plasticity in the developing gonads can persist for >145 Myr, long after the canalization of sex determination via the evolution of sex chromosomes, while some gene-specific thermal sensitivity drifts or evolves anew. Such standing thermosensitivity represents an underappreciated evolutionary potential harbored by GSD species that may be co-opted during future adaptive shifts in developmental programing, such as a GSD to TSD reversal, if favored by ecological conditions. Additionally, we identified novel candidate regulators of vertebrate sexual development in GSD reptiles, including sex-determining candidate genes in a ZZ/ZW turtle.

A Cautionary Tale of Sexing by Methylation: Hybrid Bisulfite-Conversion Sequencing of Immunoprecipitated Methylated DNA in Chrysemys picta Turtles with Temperature-Dependent Sex Determination Reveals Contrasting Patterns of Somatic and Gonadal Methylation, but No Unobtrusive Sex Diagnostic

Background: The gonads of Chrysemys picta, a turtle with temperature-dependent sex determination (TSD), exhibit differential DNA methylation between males and females, but whether the same is true in somatic tissues remains unknown. Such differential DNA methylation in the soma would provide a non-lethal sex diagnostic for TSD turtle hatchings who lack visually detectable sexual dimorphism when young. Methods: Here, we tested multiple approaches to study DNA methylation in tail clips of Chrysemys picta hatchlings, to identify differentially methylated candidate regions/sites that could serve as molecular sex markers To detect global differential methylation in the tails we used methylation-sensitive ELISA, and to test for differential local methylation we developed a novel hybrid method by sequencing immunoprecipitated and bisulfite converted DNA (MeDIP-BS-seq) followed by PCR validation of candidate regions/sites after digestion with a methylation-sensitive restriction enzyme. Results: We detected no global differences in methylation between males and females via ELISA. While we detected inter-individual variation in DNA methylation in the tails, this variation was not sexually dimorphic, in contrast with hatchling gonads. Conclusions: Results highlight that differential DNA methylation is tissue-specific and plays a key role in gonadal formation (primary sexual development) and maintenance post-hatching, but not in the somatic tail tissue.

Garcias M, Kratochvíl L, Valenzuela N, Georges A, Waters PD, Ruiz-Herrera A. Meiotic chromosome dynamics and double strand break formation in reptiles

During meiotic prophase I, tightly regulated processes take place, from pairing and synapsis of homologous chromosomes to recombination, which are essential for the generation of genetically variable haploid gametes. These processes have canonical meiotic features conserved across different phylogenetic groups. However, the dynamics of meiotic prophase I in non-mammalian vertebrates are poorly known. Here, we compare four species from Sauropsida to understand the regulation of meiotic prophase I in reptiles: the Australian central bearded dragon (Pogona vitticeps), two geckos (Paroedura picta and Coleonyx variegatus) and the painted turtle (Chrysemys picta). We first performed a histological characterization of the spermatogenesis process in both the bearded dragon and the painted turtle. We then analyzed prophase I dynamics, including chromosome pairing, synapsis and the formation of double strand breaks (DSBs). We show that meiosis progression is highly conserved in reptiles with telomeres clustering forming the bouquet, which we propose promotes homologous pairing and synapsis, along with facilitating the early pairing of micro-chromosomes during prophase I (i.e., early zygotene). Moreover, we detected low levels of meiotic DSB formation in all taxa. Our results provide new insights into reptile meiosis.

Cadmium Ecotoxic Effects on Embryonic Dmrt1 and Aromatase Expression in Chrysemys picta Turtles May Implicate Changes in DNA Methylation

Temperature-dependent sex determination (TSD) decides the sex fate of an individual based on incubation temperature. However, other environmental factors, such as pollutants, could derail TSD sexual development. Cadmium is one such contaminant of soils and water bodies known to affect DNA methylation, an epigenetic DNA modification with a key role in sexual development of TSD vertebrate embryos. Yet, whether cadmium alters DNA methylation of genes underlying gonadal formation in turtles remains unknown. Here, we investigated the effects of cadmium on the expression of two gene regulators of TSD in the painted turtle, Chrysemys picta, incubated at male-producing and female-producing temperatures using qPCR. Results revealed that cadmium alters transcription of Dmrt1 and aromatase, overriding the normal thermal effects during embryogenesis, which could potentially disrupt the sexual development of TSD turtles. Results from a preliminary DNA methylation-sensitive PCR assay implicate changes in DNA methylation of Dmrt1 as a potential cause that requires further testing (aromatase methylation assays were precluded).

Turtle Nest-Site Choice, Anthropogenic Challenges, and Evolutionary Potential for Adaptation

Oviparous animals, such as turtles, lay eggs whose success or demise depends on environmental conditions that influence offspring phenotype (morphology, physiology, and in many reptiles, also sex determination), growth, and survival, while in the nest and post-hatching. Consequently, because turtles display little parental care, maternal provisioning of the eggs and female nesting behavior are under strong selection. But the consequences of when and where nests are laid are affected by anthropogenic habitat disturbances that alter suitable nesting areas, expose eggs to contaminants in the wild, and modify the thermal and hydric environment experienced by developing embryos, thus impacting hatchling survival and the sexual fate of taxa with temperature-dependent sex determination (TSD) and genotypic sex determination (GSD). Indeed, global and local environmental change influences air, water, and soil temperature and moisture, which impact basking behavior, egg development, and conditions within the nest, potentially rendering current nesting strategies maladaptive as offspring mortality increases and TSD sex ratios become drastically skewed. Endocrine disruptors can sex reverse TSD and GSD embryos alike. Adapting to these challenges depends on genetic variation, and little to no heritability has been detected for nest-site behavior. However, modest heritability in threshold temperature (above and below which females or males develop in TSD taxa, respectively) exists in the wild, as well as interpopulation differences in the reaction norm of sex ratio to temperature, and potentially also in the expression of gene regulators of sexual development. If this variation reflects additive genetic components, some adaptation might be expected, provided that the pace of environmental change does not exceed the rate of evolution. Research remains urgently needed to fill current gaps in our understanding of the ecology and evolution of nest-site choice and its adaptive potential, integrating across multiple levels of organization.

Thermosensitive sex chromosome dosage compensation in ZZ/ZW softshell turtles, Apalone spinifera

Sex chromosome dosage compensation (SCDC) overcomes gene-dose imbalances that disturb transcriptional networks, as when ZW females or XY males are hemizygous for Z/X genes. Mounting data from non-model organisms reveal diverse SCDC mechanisms, yet their evolution remains obscure, because most informative lineages with variable sex chromosomes are unstudied. Here, we discovered SCDC in turtles and an unprecedented thermosensitive SCDC in eukaryotes. We contrasted RNA-seq expression of Z-genes, their autosomal orthologues, and control autosomal genes in Apalone spinifera (ZZ/ZW) and Chrysemys picta turtles with temperature-dependent sex determination (TSD) (proxy for ancestral expression). This approach disentangled chromosomal context effects on Z-linked and autosomal expression, from lineage effects owing to selection or drift. Embryonic Apalone SCDC is tissue- and age-dependent, regulated gene-by-gene, complete in females via Z-upregulation in both sexes (Type IV) but partial and environmentally plastic via Z-downregulation in males (accentuated at colder temperature), present in female hatchlings and a weakly suggestive in adult liver (Type I). Results indicate that embryonic SCDC evolved with/after sex chromosomes in Apalone's family Tryonichidae, while co-opting Z-gene upregulation present in the TSD ancestor. Notably, Apalone's SCDC resembles pygmy snake's, and differs from the full-SCDC of Anolis lizards who share homologous sex chromosomes (XY), advancing our understanding of how XX/XY and ZZ/ZW systems compensate gene-dose imbalance. 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 II)'.

Podocnemis expansa Turtles Hint to a Unifying Explanation for the Evolution of Temperature-Dependent Sex Determination in Long-Lived and Short-Lived Vertebrates

The adaptive significance of temperature-dependent sex determination (TSD) remains elusive for many long-lived reptiles. Various hypotheses proposed potential ecological drivers of TSD. The Charnov-Bull'77 model remains the most robust and explains the maintenance of TSD in short-lived vertebrates, where sex ratios correlate with seasonal temperatures within years that confer sex-specific fitness (colder springs produce females who grow larger and gain in fecundity, whereas warmer summers produce males who mature at smaller size). Yet, evidence of fitness differentials correlated with incubation temperature is scarce for long-lived taxa. Here, it is proposed that the Charnov-Bull'77 model applies similarly to long-lived taxa, but at a longer temporal scale, by revisiting ecological and genetic data from the long-lived turtle Podocnemis expansa. After ruling out multiple alternatives, it is hypothesized that warmer-drier years overproduce females and correlate with optimal resource availability in the flood plains, benefitting daughters more than sons, whereas resources are scarcer (due to reduced flowering/fruiting) during colder-rainier years that overproduce males, whose fitness is less impacted by slower growth rates. New technical advances and collaborative interdisciplinary efforts are delineated that should facilitate testing this hypothesis directly, illuminating the understanding of TSD evolution in P. expansa and other long-lived TSD reptiles.

Heterochromatin and microsatellites detection in karyotypes of four sea turtle species: Interspecific chromosomal differences

The wide variation in size and content of eukaryotic genomes is mainly attributed to the accumulation of repetitive DNA sequences, like microsatellites, which are tandemly repeated DNA sequences. Sea turtles share a diploid number (2n) of 56, however recent molecular cytogenetic data have shown that karyotype conservatism is not a rule in the group. In this study, the heterochromatin distribution and the chromosomal location of microsatellites (CA)n, (GA)n, (CAG)n, (GATA)n, (GAA)n, (CGC)n and (GACA)n in Chelonia mydas, Caretta caretta, Eretmochelys imbricata and Lepidochelys olivacea were comparatively investigated. The obtained data showed that just the (CA)n, (GA)n, (CAG)n and (GATA)n microsatellites were located on sea turtle chromosomes, preferentially in heterochromatic regions of the microchromosomes (mc). Variations in the location of heterochromatin and microsatellites sites, especially in some pericentromeric regions of macrochromosomes, corroborate to proposal of centromere repositioning occurrence in Cheloniidae species. Furthermore, the results obtained with the location of microsatellites corroborate with the temperature sex determination mechanism proposal and the absence of heteromorphic sex chromosomes in sea turtles. The findings are useful for understanding part of the karyotypic diversification observed in sea turtles, especially those that explain the diversification of Carettini from Chelonini species.

Effects of agricultural pesticides on the reproductive system of aquatic wildlife species, with crocodilians as sentinel species

Agricultural pesticides represent a significant class of endocrine-disrupting chemicals (EDCs) to which non-target organisms around the world are constantly exposed. Laboratory studies have found strong evidence showing the endocrine-disruptive potential of these pesticides at environmentally relevant exposure levels. Since the field of endocrine disruption continues to grow in richness and complexity, this review aims to provide an update on the effects of two agricultural pesticides that act as EDCs: atrazine and endosulfan. We will focus mainly on the effects on crocodilians due to their worldwide occurrence in tropical and sub-tropical wetland ecosystems and their ecological and physiological features, which render them vulnerable to exposure to pesticides with endocrine-disrupting action at all life stages. The results here reviewed provide important insights into the effects of hormonally active agricultural pesticides at cellular, tissue, and organ levels in the reproductive system of crocodiles. A better understanding of the effects of exposure to environmentally relevant doses of EDCs on the reproductive system of crocodilians will contribute to protect and improve the health of both wildlife species and humans.

Karyotypic Evolution of Sauropsid Vertebrates Illuminated by Optical and Physical Mapping of the Painted Turtle and Slider Turtle Genomes

Recent sequencing and software enhancements have advanced our understanding of the evolution of genomic structure and function, especially addressing novel evolutionary biology questions. Yet fragmentary turtle genome assemblies remain a challenge to fully decipher the genetic architecture of adaptive evolution. Here, we use optical mapping to improve the contiguity of the painted turtle (Chrysemys picta) genome assembly and use de novo fluorescent in situ hybridization (FISH) of bacterial artificial chromosome (BAC) clones, BAC-FISH, to physically map the genomes of the painted and slider turtles (Trachemys scripta elegans). Optical mapping increased C. picta's N50 by ~242% compared to the previous assembly. Physical mapping permitted anchoring ~45% of the genome assembly, spanning 5544 genes (including 20 genes related to the sex determination network of turtles and vertebrates). BAC-FISH data revealed assembly errors in C. picta and T. s. elegans assemblies, highlighting the importance of molecular cytogenetic data to complement bioinformatic approaches. We also compared C. picta's anchored scaffolds to the genomes of other chelonians, chicken, lizards, and snake. Results revealed a mostly one-to-one correspondence between chromosomes of painted and slider turtles, and high homology among large syntenic blocks shared with other turtles and sauropsids. Yet, numerous chromosomal rearrangements were also evident across chelonians, between turtles and squamates, and between avian and non-avian reptiles.

Molecular Cytogenetic Characterization of the Sicilian Endemic Pond Turtle Emys trinacris and the Yellow-Bellied Slider Trachemys scripta scripta (Testudines, Emydidae)

Turtles, a speciose group consisting of more than 300 species, demonstrate karyotypes with diploid chromosome numbers ranging from 2n = 26 to 2n = 68. However, cytogenetic analyses have been conducted only to 1/3rd of the turtle species, often limited to conventional staining methods. In order to expand our knowledge of the karyotype evolution in turtles, we examined the topology of the (TTAGGG)_n telomeric repeats and the rDNA loci by fluorescence in situ hybridization (FISH) on the karyotypes of two emydids: the Sicilian pond turtle, Emys trinacris, and the yellow-bellied slider, Trachemys scripta scripta (family Emydidae). Furthermore, AT-rich and GC-rich chromosome regions were detected by DAPI and CMA₃ stains, respectively. The cytogenetic analysis revealed that telomeric sequences are restricted to the terminal ends of all chromosomes and the rDNA loci are localized in one pair of microchromosomes in both species. The karyotype of the Sicilian endemic E. trinacris with diploid number 2n = 50, consisting of 13 pairs of macrochromosomes and 12 pairs of microchromosomes, is presented here for first time. Our comparative examination revealed similar cytogenetic features in Emys trinacris and the closely related E. orbicularis, as well as to other previously studied emydid species, demonstrating a low rate of karyotype evolution, as chromosomal rearrangements are rather infrequent in this group of turtles.

Turtle Insights into the Evolution of the Reptilian Karyotype and the Genomic Architecture of Sex Determination

Sex chromosome evolution remains an evolutionary puzzle despite its importance in understanding sexual development and genome evolution. The seemingly random distribution of sex-determining systems in reptiles offers a unique opportunity to study sex chromosome evolution not afforded by mammals or birds. These reptilian systems derive from multiple transitions in sex determination, some independent, some convergent, that lead to the birth and death of sex chromosomes in various lineages. Here we focus on turtles, an emerging model group with growing genomic resources. We review karyotypic changes that accompanied the evolution of chromosomal systems of genotypic sex determination (GSD) in chelonians from systems under the control of environmental temperature (TSD). These transitions gave rise to 31 GSD species identified thus far (out of 101 turtles with known sex determination), 27 with a characterized sex chromosome system (13 of those karyotypically). These sex chromosomes are varied in terms of the ancestral autosome they co-opted and thus in their homology, as well as in their size (some are macro-, some are micro-chromosomes), heterogamety (some are XX/XY, some ZZ/ZW), dimorphism (some are virtually homomorphic, some heteromorphic with larger-X, larger W, or smaller-Y), age (the oldest system could be ~195 My old and the youngest < 25 My old). Combined, all data indicate that turtles follow some tenets of classic theoretical models of sex chromosome evolution while countering others. Finally, although the study of dosage compensation and molecular divergence of turtle sex chromosomes has lagged behind research on other aspects of their evolution, this gap is rapidly decreasing with the acceleration of ongoing research and growing genomic resources in this group.

Alternative splicing and thermosensitive expression of Dmrt1 during urogenital development in the painted turtle, Chrysemys picta

BACKGROUND

The doublesex and mab-3 related transcription factor 1 (Dmrt1) is a highly conserved gene across numerous vertebrates and invertebrates in sequence and function. Small aminoacid changes in Dmrt1 are associated with turnovers in sex determination in reptiles. Dmrt1 is upregulated in males during gonadal development in many species, including the painted turtle, Chrysemys picta, a reptile with temperature-dependent sex determination (TSD). Dmrt1 is reported to play different roles during sex determination and differentiation, yet whether these functions are controlled by distinct Dmrt1 spliceoforms remains unclear. While Dmrt1 isoforms have been characterized in various vertebrates, no study has investigated their existence in any turtle.

METHODS

We examine the painted turtle to identify novel Dmrt1 isoforms that may be present during urogenital development using PCR, profile their expression by RNA-seq across five embryonic stages at male- and female-producing temperatures, and validate their expression pattern via qPCR with transcript-specific fluorescent probes.

RESULTS

A novel Dmrt1 spliceoform was discovered for the first time in chelonians, lacking exons 2 and 3 (Dmrt1 ΔEx2Ex3). Dmrt1 canonical and ΔEx2Ex3 transcripts were differentialy expressed by temperature at stages 19 and 22 in developing gonads of painted turtles, after the onset of sex determination, and displayed a significant male-biased expression pattern. This transcriptional pattern differs from studies in other turtles and vertebrates that reported Dmrt1 differential expression before or at the onset of sex determination. This study provides the first insight into Dmrt1 transcriptional diversity in turtles and opens the door for future functional studies of the alternative Dmrt1 transcript uncovered here.

CONCLUSIONS

The discovery of an isoform in turtles indicate that alternative splicing may be a common feature of Dmrt1 across vertebrates, as isoforms are also found in crocodilians, birds, mammals and fish, and this variation remains unexplained. The relatively late-onset of Dmrt1 expression observed here contrasts with other turtles, indicating that Dmrt1 is not the topmost male sex -determining factor in C. picta. When placed in a phylogenetic context, this discrepancy underscores the divergent regulation of Dmrt1, and of sexual development more generally, across vertebrates.

Comparative Cytogenetics of Four Sea Turtle Species (Cheloniidae): G-Banding Pattern and in situ Localization of Repetitive DNA Units

Sea turtles are considered flagship species for marine biodiversity conservation and are considered to be at varying risk of extinction globally. Cases of hybridism have been reported in sea turtles, but chromosomal analyses are limited to classical karyotype descriptions and a few molecular cytogenetic studies. In order to compare karyotypes and understand evolutive mechanisms related to chromosome dif­ferentiation in this group, <i>Chelonia mydas</i>, <i>Caretta caretta</i>, <i>Eretmochelys imbricata</i>, and <i>Lepidochelys olivacea</i> were cytogenetically characterized in the present study. When the obtained cytogenetic data were compared with the putative ancestral Cryptodira karyotype, the studied species showed the same diploid number (2n) of 56 chromosomes, with some variations in chromosomal morphology (karyotypic formula) and minor changes in longitudinal band locations. In situ localization using a 18S ribosomal DNA probe indicated a homeologous microchromosome pair bearing a 45S ribosomal DNA locus and size heteromorphism in all 4 species. Interstitial telomeric sites were identified in a microchromosome pair in <i>C. mydas</i> and <i>C. caretta</i>. The data showed that interspecific variations occurred in chromosomal sets among the Cheloniidae species, in addition to other Cryptodira karyotypes. These variations generated lineage-specific karyotypic diversification in sea turtles, which will have considerable implications for hybrid recognition and for the study, the biology, ecology, and evolutionary history of regional and global populations. Furthermore, we demonstrated that some chromosome rearrangements occurred in sea turtle species, which is in conflict with the hypothesis of conserved karyotypes in this group.

Cytogenetic Analysis Did Not Reveal Differentiated Sex Chromosomes in Ten Species of Boas and Pythons (Reptilia: Serpentes)

Homologous and differentiated ZZ/ZW sex chromosomes (or derived multiple neo-sex chromosomes) were often described in caenophidian snakes, but sex chromosomes were unknown until recently in non-caenophidian snakes. Previous studies revealed that two species of boas (Boa imperator, B. constrictor) and one species of python (Python bivittatus) independently evolved XX/XY sex chromosomes. In addition, heteromorphic ZZ/ZW sex chromosomes were recently revealed in the Madagascar boa (Acrantophis sp. cf. dumerili) and putatively also in the blind snake Myriopholis macrorhyncha. Since the evolution of sex chromosomes in non-caenophidian snakes seems to be more complex than previously thought, we examined ten species of pythons and boas representing the families Boidae, Calabariidae, Candoiidae, Charinidae, Pythonidae, and Sanziniidae by conventional and molecular cytogenetic methods, aiming to reveal their sex chromosomes. Our results show that all examined species do not possess sex-specific differences in their genomes detectable by the applied cytogenetic methods, indicating the presence of poorly differentiated sex chromosomes or even the absence of sex chromosomes. Interestingly, fluorescence in situ hybridization with telomeric repeats revealed extensive distribution of interstitial telomeric repeats in eight species, which are likely a consequence of intra-chromosomal rearrangements.

Discovery of Putative XX/XY Male Heterogamety in Emydura subglobosa Turtles Exposes a Novel Trajectory of Sex Chromosome Evolution in Emydura

The discovery of sex chromosome systems in non-model organisms has elicited growing recognition that sex chromosomes evolved via diverse paths that are not fully elucidated. Lineages with labile sex determination, such as turtles, hold critical cues, yet data are skewed toward hide-neck turtles (suborder Cryptodira) and scant for side-neck turtles (suborder Pleurodira). Here, we used classic and molecular cytogenetics to investigate Emydura subglobosa (ESU), an unstudied side-neck turtle with genotypic sex determination from the family Chelidae, where extensive morphological divergence exists among XX/XY systems. Our data represent the first cytogenetic description for ESU. Similarities were found between ESU and E. macquarii (EMA), such as identical chromosome number (2n = 50), a single and dimorphic nucleolus organizer region (NOR) localized in a microchromosome pair (ESU14) of both sexes (detected via FISH of 18S rDNA). Only the larger NOR is active (detected by silver staining). As in EMA, comparative genome hybridization revealed putative macro XX/XY chromosomes in ESU (the 4th largest pair). Our comparative analyses and revaluation of previous data strongly support the hypothesis that Emydura's XX/XY system evolved via fusion of an ancestral micro-Y (retained by Chelodina longicollis) onto a macro-autosome. This evolutionary trajectory differs from the purported independent evolution of XX/XY from separate ancestral autosomes in Chelodina and Emydura that was previously reported. Our data permit dating this Y-autosome fusion to at least the split of Emydura around 45 Mya and add critical information about the evolution of the remarkable diversity of sex-determining mechanisms in turtles, reptiles, and vertebrates.

Deciphering the Origin and Evolution of the X1X2Y System in Two Closely-Related Oplegnathus Species (Oplegnathidae and Centrarchiformes)

Oplegnathus fasciatus and O. punctatus (Teleostei: Centrarchiformes: Oplegnathidae), are commercially important rocky reef fishes, endemic to East Asia. Both species present an X1X2Y sex chromosome system. Here, we investigated the evolutionary forces behind the origin and differentiation of these sex chromosomes, with the aim to elucidate whether they had a single or convergent origin. To achieve this, conventional and molecular cytogenetic protocols, involving the mapping of repetitive DNA markers, comparative genomic hybridization (CGH), and whole chromosome painting (WCP) were applied. Both species presented similar 2n, karyotype structure and hybridization patterns of repetitive DNA classes. 5S rDNA loci, besides being placed on the autosomal pair 22, resided in the terminal region of the long arms of both X1 chromosomes in females, and on the X1 and Y chromosomes in males. Furthermore, WCP experiments with a probe derived from the Y chromosome of O. fasciatus (OFAS-Y) entirely painted the X1 and X2 chromosomes in females and the X1, X2, and Y chromosomes in males of both species. CGH failed to reveal any sign of sequence differentiation on the Y chromosome in both species, thereby suggesting the shared early stage of neo-Y chromosome differentiation. Altogether, the present findings confirmed the origin of the X1X2Y sex chromosomes via Y-autosome centric fusion and strongly suggested their common origin.

Extreme thermal fluctuations from climate change unexpectedly accelerate demographic collapse of vertebrates with temperature-dependent sex determination

Global climate is warming rapidly, threatening vertebrates with temperature-dependent sex determination (TSD) by disrupting sex ratios and other traits. Less understood are the effects of increased thermal fluctuations predicted to accompany climate change. Greater fluctuations could accelerate feminization of species that produce females under warmer conditions (further endangering TSD animals), or counter it (reducing extinction risk). Here we use novel experiments exposing eggs of Painted Turtles (Chrysemys picta) to replicated profiles recorded in field nests plus mathematically-modified profiles of similar shape but wider oscillations, and develop a new mathematical model for analysis. We show that broadening fluctuations around naturally male-producing (cooler) profiles feminizes developing embryos, whereas embryos from warmer profiles remain female or die. This occurs presumably because wider oscillations around cooler profiles expose embryos to very low temperatures that inhibit development, and to feminizing temperatures where most embryogenesis accrues. Likewise, embryos incubated under broader fluctuations around warmer profiles experience mostly feminizing temperatures, some dangerously high (which increase mortality), and fewer colder values that are insufficient to induce male development. Therefore, as thermal fluctuations escalate with global warming, the feminization of TSD turtle populations could accelerate, facilitating extinction by demographic collapse. Aggressive global CO₂ mitigation scenarios (RCP2.6) could prevent these risks, while intermediate actions (RCP4.5 and RCP6.0 scenarios) yield moderate feminization, highlighting the peril that insufficient reductions of greenhouse gas emissions pose for TSD taxa. If our findings are generalizable, TSD squamates, tuatara, and crocodilians that produce males at warmer temperatures could suffer accelerated masculinization, underscoring the broad taxonomic threats of climate change.

Molecular evidence for sex reversal in wild populations of green frogs (Rana clamitans)

In vertebrates, sex determination occurs along a continuum from strictly genotypic (GSD), where sex is entirely guided by genes, to strictly environmental (ESD), where rearing conditions, like temperature, determine phenotypic sex. Along this continuum are taxa which have combined genetic and environmental contributions to sex determination (GSD + EE), where some individuals experience environmental effects which cause them to sex reverse and develop their phenotypic sex opposite their genotypic sex. Amphibians are often assumed to be strictly GSD with sex reversal typically considered abnormal. Despite calls to understand the relative natural and anthropogenic causes of amphibian sex reversal, sex reversal has not been closely studied across populations of any wild amphibian, particularly in contrasting environmental conditions. Here, we use sex-linked molecular markers to discover sex reversal in wild populations of green frogs (Rana clamitans) inhabiting ponds in either undeveloped, forested landscapes or in suburban neighborhoods. Our work here begins to suggest that sex reversal may be common within and across green frog populations, occurring in 12 of 16 populations and with frequencies of 2–16% of individuals sampled within populations. Additionally, our results also suggest that intersex phenotypic males and sex reversal are not correlated with each other and are also not correlated with suburban land use. While sex reversal and intersex are often considered aberrant responses to human activities and associated pollution, we found no such associations here. Our data perhaps begin to suggest that, relative to what is often suggested, sex reversal may be a relatively natural process in amphibians. Future research should focus on assessing interactions between genes and the environment to understand the molecular and exogenous basis of sex determination in green frogs and in other amphibians.

Insights into Emydid Turtle Cytogenetics: The European Pond Turtle as a Model Species

Our knowledge of Testudines evolution is limited by the lack of modern cytogenetic data. Compared to other reptiles, there is little information even on chromosome banding, let alone molecular cytogenetic data. Here, we provide detailed information on the karyotype of the European pond turtle Emys orbicularis, a model Emydidae, employing both chromosome banding and molecular cytogenetics. We provide a high-resolution G-banded karyotype and a map of rDNA genes and telomeric sequences using fluorescence in situ hybridization. We test hypotheses of sex-determining mechanisms in Emys by comparative genomic hybridization to determine if Emys has a cryptic sex-specific region. Our results provide valuable data to guide future efforts on genome sequencing and anchoring in Emydidae and for understanding karyotype evolution in Testudines.

Thermal Response of Epigenetic Genes Informs Turtle Sex Determination with and without Sex Chromosomes

Vertebrate sexual fate can be established by environmental cues (e.g., temperature-dependent sex determination, TSD) or by genetic content (genotypic sex determination, GSD). While methylation is implicated in TSD, the influence of broader epigenetic processes in sexual development remains obscure. Here, we investigated for the first time the embryonic gonadal expression of the genome-wide epigenetic machinery in turtles, including genes and noncoding RNAs (ncRNAs) involved in DNA/histone acetylation, methylation, ubiquitination, phosphorylation, and RNAi. This machinery was active and differentially thermosensitive in TSD versus GSD (ZZ/ZW) turtles. Methylation and histone acetylation genes responded the strongest. The results suggest these working hypotheses: (i) TSD might be mediated by epigenetically controlled hormonal pathways (via acetylation, methylation, and ncRNAs), or by (ii) hormonally controlled epigenetic processes, and (iii) key epigenetic events prior to the canonical thermosensitive period may explain differences between TSD and GSD. Novel epigenetic candidate regulators other than methylation were identified, including previously unknown ncRNAs that could potentially mediate gonadogenesis. These findings illuminate the molecular ecology of reptilian sex determination and permitted hypothesis building to help guide future functional studies on the epigenetic transduction of external cues in TSD versus GSD systems.

Conservation of Sex-Linked Markers among Conspecific Populations of a Viviparous Skink, Niveoscincus ocellatus, Exhibiting Genetic and Temperature-Dependent Sex Determination

Sex determination systems are exceptionally diverse and have undergone multiple and independent evolutionary transitions among species, particularly reptiles. However, the mechanisms underlying these transitions have not been established. Here, we tested for differences in sex-linked markers in the only known reptile that is polymorphic for sex determination system, the spotted snow skink, Niveoscincus ocellatus, to quantify the genomic differences that have accompanied this transition. In a highland population, sex is determined genetically, whereas in a lowland population, offspring sex ratio is influenced by temperature. We found a similar number of sex-linked loci in each population, including shared loci, with genotypes consistent with male heterogamety (XY). However, population-specific linkage disequilibrium suggests greater differentiation of sex chromosomes in the highland population. Our results suggest that transitions between sex determination systems can be facilitated by subtle genetic differences.

Chromosomal Context Affects the Molecular Evolution of Sex-linked Genes and Their Autosomal Counterparts in Turtles and Other Vertebrates

Sex chromosomes evolve differently from autosomes because natural selection acts distinctly on them given their reduced recombination and smaller population size. Various studies of sex-linked genes compared with different autosomal genes within species support these predictions. Here, we take a novel alternative approach by comparing the rate of evolution between subsets of genes that are sex-linked in selected reptiles/vertebrates and the same genes located in autosomes in other amniotes. We report for the first time the faster evolution of Z-linked genes in a turtle (the Chinese softshell turtle Pelodiscus sinensis) relative to autosomal orthologs in other taxa, including turtles with temperature-dependent sex determination (TSD). This faster rate was absent in its close relative, the spiny softshell turtle (Apalone spinifera), thus revealing important lineage effects, and was only surpassed by mammalian-X linked genes. In contrast, we found slower evolution of X-linked genes in the musk turtle Staurotypus triporcatus (XX/XY) and homologous Z-linked chicken genes. TSD lineages displayed overall faster sequence evolution than taxa with genotypic sex determination (GSD), ruling out global effects of GSD on molecular evolution beyond those by sex-linkage. Notably, results revealed a putative selective sweep around two turtle genes involved in vertebrate gonadogenesis (Pelodiscus-Z-linked Nf2 and Chrysemys-autosomal Tspan7). Our observations reveal important evolutionary changes at the gene level mediated by chromosomal context in turtles despite their low overall evolutionary rate and illuminate sex chromosome evolution by empirically testing expectations from theoretical models. Genome-wide analyses are warranted to test the generality and prevalence of the observed patterns.

Temperature- vs. estrogen-induced sex determination in Caiman latirostris embryos: Both females, but with different expression patterns of key molecules involved in ovarian development

Caiman latirostris is a species with temperature dependent sex determination (TSD), which implies that the incubation temperature of the eggs is the main factor that determines the sex during a thermo-sensitive period (TSP). However, estrogens play a critical role in this process. The administration of 17β-estradiol (E₂) previous to TSP overrides the effects of male incubation temperature, producing phenotypic females. This effect has been defined as sex reversal or estrogen-induced sex determination (E₂SD). The aim of the present study is to describe similarities and differences in the effects of TSD and E₂SD treatment conditions on ovary development. Our results show that the two treatment conditions studied are able to produce different ovaries. Treatment with E₂ modified the expression pattern of estrogen receptor alpha and progesterone receptor, and expression of the enzyme aromatase. Moreover, in E₂SD females, the proliferation/apoptosis dynamic was also altered and high expression of TAp63 was observed suggesting the presence of greater DNA damage in germ cells. To the best of our knowledge, this is the first report that describes the morphology of the female gonad of C. latirostris in three stages of embryonic development and shows the expression of TAp63 during the gonad development of a reptile. It is important to emphasize that the changes demonstrated in E₂SD female gonads of embryos show that environmental compounds with proven estrogenic activity alter the follicular dynamics of C. latirostris in neonatal as much as in juvenile animals, endangering their reproductive health and possibly bringing consequences to ecology and evolution.

MeDIP-seq and nCpG analyses illuminate sexually dimorphic methylation of gonadal development genes with high historic methylation in turtle hatchlings with temperature-dependent sex determination

Background

DNA methylation alters gene expression but not DNA sequence and mediates some cases of phenotypic plasticity. Temperature-dependent sex determination (TSD) epitomizes phenotypic plasticity where environmental temperature drives embryonic sexual fate, as occurs commonly in turtles. Importantly, the temperature-specific transcription of two genes underlying gonadal differentiation is known to be induced by differential methylation in TSD fish, turtle and alligator. Yet, how extensive is the link between DNA methylation and TSD remains unclear. Here we test for broad differences in genome-wide DNA methylation between male and female hatchling gonads of the TSD painted turtle Chrysemys picta using methyl DNA immunoprecipitation sequencing, to identify differentially methylated candidates for future study. We also examine the genome-wide nCpG distribution (which affects DNA methylation) in painted turtles and test for historic methylation in genes regulating vertebrate gonadogenesis.

Results

Turtle global methylation was consistent with other vertebrates (57% of the genome, 78% of all CpG dinucleotides). Numerous genes predicted to regulate turtle gonadogenesis exhibited sex-specific methylation and were proximal to methylated repeats. nCpG distribution predicted actual turtle DNA methylation and was bimodal in gene promoters (as other vertebrates) and introns (unlike other vertebrates). Differentially methylated genes, including regulators of sexual development, had lower nCpG content indicative of higher historic methylation.

Conclusions

Ours is the first evidence suggesting that sexually dimorphic DNA methylation is pervasive in turtle gonads (perhaps mediated by repeat methylation) and that it targets numerous regulators of gonadal development, consistent with the hypothesis that it may regulate thermosensitive transcription in TSD vertebrates. However, further research during embryogenesis will help test this hypothesis and the alternative that instead, most differential methylation observed in hatchlings is the by-product of sexual differentiation and not its cause.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-017-0136-2) contains supplementary material, which is available to authorized users.

Transcriptomic responses to environmental temperature by turtles with temperature-dependent and genotypic sex determination assessed by RNAseq inform the genetic architecture of embryonic gonadal development

Vertebrate sexual fate is decided primarily by the individual's genotype (GSD), by the environmental temperature during development (TSD), or both. Turtles exhibit TSD and GSD, making them ideal to study the evolution of sex determination. Here we analyze temperature-specific gonadal transcriptomes (RNA-sequencing validated by qPCR) of painted turtles (Chrysemys picta TSD) before and during the thermosensitive period, and at equivalent stages in soft-shell turtles (Apalone spinifera-GSD), to test whether TSD's and GSD's transcriptional circuitry is identical but deployed differently between mechanisms. Our data show that most elements of the mammalian urogenital network are active during turtle gonadogenesis, but their transcription is generally more thermoresponsive in TSD than GSD, and concordant with their sex-specific function in mammals [e.g., upregulation of Amh, Ar, Esr1, Fog2, Gata4, Igf1r, Insr, and Lhx9 at male-producing temperature, and of β-catenin, Foxl2, Aromatase (Cyp19a1), Fst, Nf-kb, Crabp2 at female-producing temperature in Chrysemys]. Notably, antagonistic elements in gonadogenesis (e.g., β-catenin and Insr) were thermosensitive only in TSD early-embryos. Cirbp showed warm-temperature upregulation in both turtles disputing its purported key TSD role. Genes that may convert thermal inputs into sex-specific development (e.g., signaling and hormonal pathways, RNA-binding and heat-shock) were differentially regulated. Jak-Stat, Nf-κB, retinoic-acid, Wnt, and Mapk-signaling (not Akt and Ras-signaling) potentially mediate TSD thermosensitivity. Numerous species-specific ncRNAs (including Xist) were differentially-expressed, mostly upregulated at colder temperatures, as were unannotated loci that constitute novel TSD candidates. Cirbp showed warm-temperature upregulation in both turtles. Consistent transcription between turtles and alligator revealed putatively-critical reptilian TSD elements for male (Sf1, Amh, Amhr2) and female (Crabp2 and Hspb1) gonadogenesis. In conclusion, while preliminary, our data helps illuminate the regulation and evolution of vertebrate sex determination, and contribute genomic resources to guide further research into this fundamental biological process.

Cytogenetic Insights into the Evolution of Chromosomes and Sex Determination Reveal Striking Homology of Turtle Sex Chromosomes to Amphibian Autosomes

Turtle karyotypes are highly conserved compared to other vertebrates; yet, variation in diploid number (2n = 26-68) reflects profound genomic reorganization, which correlates with evolutionary turnovers in sex determination. We evaluate the published literature and newly collected comparative cytogenetic data (G- and C-banding, 18S-NOR, and telomere-FISH mapping) from 13 species spanning 2n = 28-68 to revisit turtle genome evolution and sex determination. Interstitial telomeric sites were detected in multiple lineages that underwent diploid number and sex determination turnovers, suggesting chromosomal rearrangements. C-banding revealed potential interspecific variation in centromere composition and interstitial heterochromatin at secondary constrictions. 18S-NORs were detected in secondary constrictions in a single chromosomal pair per species, refuting previous reports of multiple NORs in turtles. 18S-NORs are linked to ZW chromosomes in Apalone and Pelodiscus and to X (not Y) in Staurotypus. Notably, comparative genomics across amniotes revealed that the sex chromosomes of several turtles, as well as mammals and some lizards, are homologous to components of Xenopus tropicalis XTR1 (carrying Dmrt1). Other turtle sex chromosomes are homologous to XTR4 (carrying Wt1). Interestingly, all known turtle sex chromosomes, except in Trionychidae, evolved via inversions around Dmrt1 or Wt1. Thus, XTR1 appears to represent an amniote proto-sex chromosome (perhaps linked ancestrally to XTR4) that gave rise to turtle and other amniote sex chromosomes.

Sex determination, longevity, and the birth and death of reptilian species

Vertebrate sex-determining mechanisms (SDMs) are triggered by the genotype (GSD), by temperature (TSD), or occasionally, by both. The causes and consequences of SDM diversity remain enigmatic. Theory predicts SDM effects on species diversification, and life-span effects on SDM evolutionary turnover. Yet, evidence is conflicting in clades with labile SDMs, such as reptiles. Here, we investigate whether SDM is associated with diversification in turtles and lizards, and whether alterative factors, such as lifespan's effect on transition rates, could explain the relative prevalence of SDMs in turtles and lizards (including and excluding snakes). We assembled a comprehensive dataset of SDM states for squamates and turtles and leveraged large phylogenies for these two groups. We found no evidence that SDMs affect turtle, squamate, or lizard diversification. However, SDM transition rates differ between groups. In lizards TSD-to-GSD surpass GSD-to-TSD transitions, explaining the predominance of GSD lizards in nature. SDM transitions are fewer in turtles and the rates are similar to each other (TSD-to-GSD equals GSD-to-TSD), which, coupled with TSD ancestry, could explain TSD's predominance in turtles. These contrasting patterns can be explained by differences in life history. Namely, our data support the notion that in general, shorter lizard lifespan renders TSD detrimental favoring GSD evolution in squamates, whereas turtle longevity permits TSD retention. Thus, based on the macro-evolutionary evidence we uncovered, we hypothesize that turtles and lizards followed different evolutionary trajectories with respect to SDM, likely mediated by differences in lifespan. Combined, our findings revealed a complex evolutionary interplay between SDMs and life histories that warrants further research that should make use of expanded datasets on unexamined taxa to enable more conclusive analyses.

Physical Mapping and Refinement of the Painted Turtle Genome (Chrysemys picta) Inform Amniote Genome Evolution and Challenge Turtle-Bird Chromosomal Conservation

Comparative genomics continues illuminating amniote genome evolution, but for many lineages our understanding remains incomplete. Here, we refine the assembly (CPI 3.0.3 NCBI AHGY00000000.2) and develop a cytogenetic map of the painted turtle (Chrysemys picta-CPI) genome, the first in turtles and in vertebrates with temperature-dependent sex determination. A comparison of turtle genomes with those of chicken, selected nonavian reptiles, and human revealed shared and novel genomic features, such as numerous chromosomal rearrangements. The largest conserved syntenic blocks between birds and turtles exist in four macrochromosomes, whereas rearrangements were evident in these and other chromosomes, disproving that turtles and birds retain fully conserved macrochromosomes for greater than 300 Myr. C-banding revealed large heterochromatic blocks in the centromeric region of only few chromosomes. The nucleolar-organizing region (NOR) mapped to a single CPI microchromosome, whereas in some turtles and lizards the NOR maps to nonhomologous sex-chromosomes, thus revealing independent translocations of the NOR in various reptilian lineages. There was no evidence for recent chromosomal fusions as interstitial telomeric-DNA was absent. Some repeat elements (CR1-like, Gypsy) were enriched in the centromeres of five chromosomes, whereas others were widespread in the CPI genome. Bacterial artificial chromosome (BAC) clones were hybridized to 18 of the 25 CPI chromosomes and anchored to a G-banded ideogram. Several CPI sex-determining genes mapped to five chromosomes, and homology was detected between yet other CPI autosomes and the globally nonhomologous sex chromosomes of chicken, other turtles, and squamates, underscoring the independent evolution of vertebrate sex-determining mechanisms.