Molecular and morphological differentiation of testes and ovaries in relation to the thermosensitive period of gonad development in the snapping turtle, Chelydra serpentina

Gene expression dynamics during temperature-dependent sex determination in a sea turtle

Fifty years ago, researchers discovered a link between ambient temperature and the sex of turtle embryos. More recently, significant progress has been made in understanding the influence of temperature on freshwater turtles. However, our understanding of the key genetic factors in other turtle groups, such as sea turtles, remains limited. To address this gap, we conducted RNA-seq analyses on embryonic tissues from the sea olive ridley turtle during the thermosensitive period (stages 21-26) at temperatures known to produce males (26 °C) and females (33 °C). Our findings revealed that incubation temperatures primarily influence genes with broad expression across tissues due to differential cell division rates and later have an effect regulating gonad-specific transcripts. This effect is mostly related to gene activation rather than transcription repression. We performed transcriptome analyses following shifts in incubation temperatures of bi-potential gonads. This approach allowed us to identify genes that respond rapidly and may be closer to the beginning of the temperature-sensing pathway. Notably, we observed swift adaptations in the expression levels of chromatin modifiers JARID2 and KDM6B, as well as the splicing factor SRSF5, and transcription regulators THOC2, DDX3X and CBX3, but little impact in the overall gonad-specific pathways, indicating that temperature-sensing genes may change rapidly but the rewiring of the gonad's developmental fate is complex and resilient. AUTHOR SUMMARY: Sea turtles, one of the most iconic creatures of our oceans, confront a troubling reality of endangerment, a peril magnified by the looming specter of climate change. This climatic shift is gradually increasing the temperature of the nesting beaches thus causing dramatic male/female population biases. Conservation efforts will need genetic and molecular information to reverse the negative effects of climate change on the populations. In this study, we conducted the first transcriptomic analysis of embryonic tissues, including gonads, brain, liver, and mesonephros, in the olive ridley sea turtle during the critical thermosensitive period spanning stages 21-26. We examined both male-producing (26 °C) and female-producing (33 °C) temperatures and found that incubation temperatures influence temperature-sensitive genes that are either expressed globally or specifically associated with the gonads. These findings indicate that incubation temperatures predominantly sway genes with broad expression patterns due to differential cell division rates. This natural process was opted in the gonads to drive sex determination. We also identified genes that are rapidly capable of sensing temperature changes and that could play a role in the activation of the sex determination pathway. Overall, our study sheds light on the intricate interplay between temperature and gene expression during sea turtle development, revealing dynamic changes in the transcriptome and highlighting the involvement of key genetic players in sex determination.

Histomorphic analysis and expression of mRNA and miRNA in embryonic gonadal differentiation in Chinese soft-shelled turtle (Pelodiscus sinensis)

The Chinese soft-shelled turtle (Pelodiscus sinensis) is extensively cultured in Asia for its nutritional and medical value. Gonadal differentiation is fantastic in turtles, whereas morphologic, mRNA, and miRNA expressions were insufficient in the turtle. In this study, ovaries and testes histomorphology analysis of 14-23 stage embryos were performed, and mRNA and miRNA expression profiles were analyzed. Histomorphology analysis revealed that gonads were undifferentiated at embryonic stage 14. Ovarian morphological differentiation became evident from stage 15, which was characterized by the development of the cortical region and degeneration of the medullary region. Concurrently, testicular morphological differentiation was apparent from stage 15, marked by the development of the medullary region and degeneration of the cortical region. qRT-PCR results showed that Cyp19a1 and Foxl2 exhibited female-specific expression at stage 15 and the expression increased throughout most of the embryonic development. Dmrt1, Amh, and Sox9 displayed male-specific expression at stage 15 and tended to increase substantially at later developmental stages. The expression of miR-8356 and miR-3299 in ZZ gonads were significantly higher than that in ZW gonads at stage 15, 17 and 19, and they had the highest expression at stage 15. While the expression of miR-8085 and miR-7982 had the highest expression at stage 19. Furthermore, chromatin remodeler genes showed differential expression in female and male P. sinensis gonads. These results of master sex-differentiation genes and morphological characteristics would provide a reference for the research of sex differentiation and sex reversal in turtles. Additionally, the expression of chromatin remodeler genes indicated they might be involved in gonadal differentiation of P. sinensis.

The -KTS splice variant of WT1 is essential for ovarian determination in mice

Sex determination in mammals depends on the differentiation of the supporting lineage of the gonads into Sertoli or pregranulosa cells that govern testis and ovary development, respectively. Although the Y-linked testis-determining gene Sry has been identified, the ovarian-determining factor remains unknown. In this study, we identified -KTS, a major, alternatively spliced isoform of the Wilms tumor suppressor WT1, as a key determinant of female sex determination. Loss of -KTS variants blocked gonadal differentiation in mice, whereas increased expression, as found in Frasier syndrome, induced precocious differentiation of ovaries independently of their genetic sex. In XY embryos, this antagonized Sry expression, resulting in male-to-female sex reversal. Our results identify -KTS as an ovarian-determining factor and demonstrate that its time of activation is critical in gonadal sex differentiation.

Testicular heat stress, a historical perspective and two postulates for why male germ cells are heat sensitive

Herein, we compare the different experimental regimes used to induce testicular heat stress and summarise their impact on sperm production and male fertility. Irrespective of the protocol used, scrotal heat stress causes loss of sperm production. This is first seen 1-2 weeks post heat stress, peaking 4-5 weeks thereafter. The higher the temperature, or the longer the duration of heat, the more pronounced germ cell loss becomes, within extreme cases this leads to azoospermia. The second, and often underappreciated impact of testicular hyperthermia is the production of poor-quality spermatozoa. Typically, those cells that survive hyperthermia develop into morphologically abnormal and poorly motile spermatozoa. While both apoptotic and non-apoptotic pathways are known to contribute to hyperthermic germ cell loss, the mechanisms leading to formation of poor-quality sperm remain unclear. Mechanistically, it is unlikely that testicular hyperthermia affects messenger RNA (mRNA) abundance, as a comparison of four different mammalian studies shows no consistent single gene changes. Using available evidence, we propose two novel models to explain how testicular hyperthermia impairs sperm formation. Our first model suggests aberrant alternative splicing, while the second model proposes a loss of RNA repression. Importantly, neither model requires consistent changes in RNA species.

Sex Determination and Ovarian Development in Reptiles and Amphibians: From Genetic Pathways to Environmental Influences

BACKGROUND

Reptiles and amphibians provide untapped potential for discovering how a diversity of genetic pathways and environmental conditions are incorporated into developmental processes that can lead to similar functional outcomes. These groups display a multitude of reproductive strategies, and whereas many attributes are conserved within groups and even across vertebrates, several aspects of sexual development show considerable variation.

SUMMARY

In this review, we focus our attention on the development of the reptilian and amphibian ovary. First, we review and describe the events leading to ovarian development, including sex determination and ovarian maturation, through a comparative lens. We then describe how these events are influenced by environmental factors, focusing on temperature and exposure to anthropogenic chemicals. Lastly, we identify critical knowledge gaps and future research directions that will be crucial to moving forward in our understanding of ovarian development and the influences of the environment in reptiles and amphibians.

KEY MESSAGES

Reptiles and amphibians provide excellent models for understanding the diversity of sex determination strategies and reproductive development. However, a greater understanding of the basic biology of these systems is necessary for deciphering the adaptive and potentially disruptive implications of embryo-by-environment interactions in a rapidly changing world.

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.

Role of Alternative Splicing in Sex Determination in Vertebrates

During the process of sex determination, a germ-cell-containing undifferentiated gonad is converted into either a male or a female reproductive organ. Both the composition of sex chromosomes and the environment determine sex in vertebrates. It is assumed that transcription level regulation drives this cascade of mechanisms; however, transcription factors can alter gene expression beyond transcription initiation by controlling pre-mRNA splicing and thereby mRNA isoform production. Using the key time window in sex determination and gonad development in mice, it has been reported that new non-transcriptional events, such as alternative splicing, could play a key role in sex determination in mammals. We know the role of key regulatory factors, like WT1(+/-KTS) or FGFR2(b/c) in pre-mRNA splicing and sex determination, indicating that important steps in the vertebrate sex determination process probably operate at a post-transcriptional level. Here, we discuss the role of pre-mRNA splicing regulators in sex determination in vertebrates, focusing on the new RNA-seq data reported from mice fetal gonadal transcriptome.

Evolutionary Turnover in Wnt Gene Expression but Conservation of Wnt Signaling during Ovary Determination in a TSD Reptile

Temperature-dependent sex determination (TSD) is a well-known characteristic of many reptilian species. However, the molecular processes linking ambient temperature to determination of gonad fate remain hazy. Here, we test the hypothesis that Wnt expression and signaling differ between female- and male-producing temperatures in the snapping turtle Chelydra serpentina. Canonical Wnt signaling involves secretion of glycoproteins called WNTs, which bind to and activate membrane bound receptors that trigger β-catenin stabilization and translocation to the nucleus where β-catenin interacts with TCF/LEF transcription factors to regulate expression of Wnt targets. Non-canonical Wnt signaling occurs via 2 pathways that are independent of β-catenin: one involves intracellular calcium release (the Wnt/Ca2+ pathway), while the other involves activation of RAC1, JNK, and RHOA (the Wnt/planar cell polarity pathway). We screened 20 Wnt genes for differential expression between female- and male-producing temperatures during sex determination in the snapping turtle. Exposure of embryos to the female-producing temperature decreased expression of 7 Wnt genes but increased expression of 2 Wnt genes and Rspo1 relative to embryos at the male-producing temperature. Temperature also regulated expression of putative Wnt target genes in vivo and a canonical Wnt reporter (6x TCF/LEF sites drive H2B-GFP expression) in embryonic gonadal cells in vitro. Results indicate that Wnt signaling was higher at the female- than at the male-producing temperature. Evolutionary analyses of all 20 Wnt genes revealed that thermosensitive Wnts, as opposed to insensitive Wnts, were less likely to show evidence of positive selection and experienced stronger purifying selection within TSD species.

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.

Temperature-Dependent Alternative Splicing of Precursor mRNAs and Its Biological Significance: A Review Focused on Post-Transcriptional Regulation of a Cold Shock Protein Gene in Hibernating Mammals

Multiple mRNA isoforms are often generated during processing such as alternative splicing of precursor mRNAs (pre-mRNA), resulting in a diversity of generated proteins. Alternative splicing is an essential mechanism for the functional complexity of eukaryotes. Temperature, which is involved in all life activities at various levels, is one of regulatory factors for controlling patterns of alternative splicing. Temperature-dependent alternative splicing is associated with various phenotypes such as flowering and circadian clock in plants and sex determination in poikilothermic animals. In some specific situations, temperature-dependent alternative splicing can be evoked even in homothermal animals. For example, the splicing pattern of mRNA for a cold shock protein, cold-inducible RNA-binding protein (CIRP or CIRBP), is changed in response to a marked drop in body temperature during hibernation of hamsters. In this review, we describe the current knowledge about mechanisms and functions of temperature-dependent alternative splicing in plants and animals. Then we discuss the physiological significance of hypothermia-induced alternative splicing of a cold shock protein gene in hibernating and non-hibernating animals.

Post-Transcriptional Mechanisms Respond Rapidly to Ecologically Relevant Thermal Fluctuations During Temperature-Dependent Sex Determination

An organism's ability to integrate transient environmental cues experienced during development into molecular and physiological responses forms the basis for adaptive shifts in phenotypic trajectories. During temperature-dependent sex determination (TSD), thermal cues during discrete periods in development coordinate molecular changes that ultimately dictate sexual fate and contribute to patterns of inter- and intra-sexual variation. How these mechanisms interface with dynamic thermal environments in nature remain largely unknown. By deploying thermal loggers in wild nests of the American alligator (Alligator mississippiensis) over two consecutive breeding seasons, we observed that 80% of nests exhibit both male- and female-promoting thermal cues during the thermosensitive period, and of these nests, all exhibited both male- and female-promoting temperatures within the span of a single day. These observations raise a critical question-how are opposing environmental cues integrated into sexually dimorphic transcriptional programs across short temporal scales? To address this question, alligator embryos were exposed to fluctuating temperatures based on nest thermal profiles and sampled over the course of a daily thermal fluctuation. We examined the expression dynamics of upstream genes in the temperature-sensing pathway and find that post-transcriptional alternative splicing and transcript abundance of epigenetic modifier genes JARID2 and KDM6B respond rapidly to thermal fluctuations while transcriptional changes of downstream effector genes, SOX9 and DMRT1, occur on a delayed timescale. Our findings reveal how the basic mechanisms of TSD operate in an ecologically relevant context. We present a hypothetical hierarchical model based on our findings as well as previous studies, in which temperature-sensitive alternative splicing incrementally influences the epigenetic landscape to affect the transcriptional activity of key sex-determining genes.

Thermosensitive period for sex determination of the tropical freshwater turtle Malayemys macrocephala

Many egg-laying reptiles possess temperature-dependent sex determination (TSD) in which outcome of gonadogenesis is determined by incubation temperature during a temperature-sensitive period of development. Prior studies on Malayemys macrocephala showed that incubation temperatures influence gonadal development and suggested that M. macrocephala exhibits TSD. However, information on the temperature-sensitivity period in this species was unknown until the current study. Turtle eggs were collected from rice fields in central Thailand from December 2016 to February 2017. In the laboratory, eggs were incubated at male-biased temperature (26 °C) and shifted to female-biased temperature (32 °C), or vice versa. Single shift experiments were performed systematically during embryonic stages 13-21. After hatching, sex of individual turtles was determined by histological analysis. We found that the sex determination of M. macrocephala is affected by temperature up to stage 16 of embryonic development.

Embryonic Temperature Programs Phenotype in Reptiles

Reptiles are critically affected by temperature throughout their lifespan, but especially so during early development. Temperature-induced changes in phenotype are a specific example of a broader phenomenon called phenotypic plasticity in which a single individual is able to develop different phenotypes when exposed to different environments. With climate change occurring at an unprecedented rate, it is important to study temperature effects on reptiles. For example, the potential impact of global warming is especially pronounced in species with temperature-dependent sex determination (TSD) because temperature has a direct effect on a key phenotypic (sex) and demographic (population sex ratios) trait. Reptiles with TSD also serve as models for studying temperature effects on the development of other traits that display continuous variation. Temperature directly influences metabolic and developmental rate of embryos and can have permanent effects on phenotype that last beyond the embryonic period. For instance, incubation temperature programs post-hatching hormone production and growth physiology, which can profoundly influence fitness. Here, we review current knowledge of temperature effects on phenotypic and developmental plasticity in reptiles. First, we examine the direct effect of temperature on biophysical processes, the concept of thermal performance curves, and the process of thermal acclimation. After discussing these reversible temperature effects, we focus the bulk of the review on developmental programming of phenotype by temperature during embryogenesis (i.e., permanent developmental effects). We focus on oviparous species because eggs are especially susceptible to changes in ambient temperature. We then discuss recent work probing the role of epigenetic mechanisms in mediating temperature effects on phenotype. Based on phenotypic effects of temperature, we return to the potential impact of global warming on reptiles. Finally, we highlight key areas for future research, including the identification of temperature sensors and assessment of genetic variation for thermosensitivity.

Role for androgens in determination of ovarian fate in the common snapping turtle, Chelydra serpentina

Sex steroids are involved in sex determination in almost all vertebrates, including species with temperature-dependent sex determination (TSD). It is well established that aromatase and estrogens are involved in ovary determination in TSD species. In contrast, the role of non-aromatizable androgens in TSD is less clear. In this study, we used dihydrotestosterone (DHT) and an antagonist of the mammalian androgen receptor (flutamide) to examine the impact of androgens on sex determination in the snapping turtle. We incubated eggs at a male-producing temperature and treated embryos with drug delivery vehicle (5 L ethanol), DHT in vehicle, or flutamide in vehicle during the sex-determining period. We then measured expression of markers for ovarian and testicular development and genes involved in steroidogenesis. A subset of embryos and hatchlings were collected for histological analysis of gonad differentiation and sex determination. DHT and flutamide both induced ovarian development: 100% of vehicle-treated hatchlings had testes, while 60% of DHT-treated and 32% flutamide-treated hatchlings had ovaries. DHT and flutamide treatments also had feminizing effects on gene expression patterns and the structure of embryonic gonads. DHT treatment increased expression of FoxL2, androgen receptor, aromatase and several steroidogenic genes. Flutamide produced a similar, but weaker, pattern of gene expression. Genes involved in testis development (Sox9 and Amh) were influenced by flutamide treatment. Our findings support the hypothesis that androgens and the androgen receptor are involved in ovary determination in the common snapping turtle.

Thyroid hormone modulates offspring sex ratio in a turtle with temperature-dependent sex determination

The adaptive significance of temperature-dependent sex determination (TSD) has attracted a great deal of research, but the underlying mechanisms by which temperature determines the sex of a developing embryo remain poorly understood. Here, we manipulated the level of a thyroid hormone (TH), triiodothyronine (T₃), during embryonic development (by adding excess T₃ to the eggs of the red-eared slider turtle Trachemys scripta, a reptile with TSD), to test two competing hypotheses on the proximate basis for TSD: the developmental rate hypothesis versus the hormone hypothesis Exogenous TH accelerated embryonic heart rate (and hence metabolic rate), developmental rate, and rates of early post-hatching growth. More importantly, hyperthyroid conditions depressed expression of Cyp19a1 (the gene encoding for aromatase) and levels of oestradiol, and induced more male offspring. This result is contrary to the direction of sex-ratio shift predicted by the developmental rate hypothesis, but consistent with that predicted by the hormone hypothesis Our results suggest an important role for THs in regulating sex steroid hormones, and therefore, in affecting gonadal sex differentiation in TSD reptiles. Our study has implications for the conservation of TSD reptiles in the context of global change because environmental contaminants may disrupt the activity of THs, and thereby affect offspring sex in TSD reptiles.

Dimorphic DNA methylation during temperature-dependent sex determination in the sea turtle Lepidochelys olivacea

Sex determination in vertebrates depends on the expression of a conserved network of genes. Sea turtles such as Lepidochelys olivacea have temperature-dependent sex determination. The present work analyses some of the epigenetic processes involved in this. We describe sexual dimorphism in global DNA methylation patterns between ovaries and testes of L. olivacea and show that the differences may arise from a combination of DNA methylation and demethylation events that occur during sex determination. Irrespective of incubation temperature, 5-hydroxymethylcytosine was abundant in the bipotential gonad; however, following sex determination, this modification was no longer found in pre-Sertoli cells in the testes. These changes correlate with the establishment of the sexually dimorphic DNA methylation patterns, down regulation of Sox9 gene expression in ovaries and irreversible gonadal commitment towards a male or female differentiation pathway. Thus, DNA methylation changes may be necessary for the stabilization of the gene expression networks that drive the differentiation of the bipotential gonad to form either an ovary or a testis in L. olivacea and probably among other species that manifest temperature-dependent sex determination.

A Novel Candidate Gene for Temperature-Dependent Sex Determination in the Common Snapping Turtle

Temperature-dependent sex determination (TSD) was described nearly 50 years ago. Researchers have since identified many genes that display differential expression at male- vs. female-producing temperatures. Yet, it is unclear whether these genes (1) are involved in sex determination per se, (2) are downstream effectors involved in differentiation of ovaries and testes, or (3) are thermo-sensitive but unrelated to gonad development. Here we present multiple lines of evidence linking CIRBP to sex determination in the snapping turtle, Chelydra serpentina We demonstrate significant associations between a single nucleotide polymorphism (SNP) (c63A > C) in CIRBP, transcript levels in embryonic gonads during specification of gonad fate, and sex in hatchlings from a thermal regime that produces mixed sex ratios. The A allele was induced in embryos exposed to a female-producing temperature, while expression of the C allele did not differ between female- and male-producing temperatures. In accord with this pattern of temperature-dependent, allele-specific expression, AA homozygotes were more likely to develop ovaries than AC heterozygotes, which, in turn, were more likely to develop ovaries than CC homozygotes. Multiple regression using SNPs in CIRBP and adjacent loci suggests that c63A > C may be the causal variant or closely linked to it. Differences in CIRBP allele frequencies among turtles from northern Minnesota, southern Minnesota, and Texas reflect small and large-scale latitudinal differences in TSD pattern. Finally, analysis of CIRBP protein localization reveals that CIRBP is in a position to mediate temperature effects on the developing gonads. Together, these studies strongly suggest that CIRBP is involved in determining the fate of the bipotential gonad.