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

Using a Handful of Transcriptomes to Detect Sex-Linked Markers and Develop Molecular Sexing Assays in a Species with Homomorphic Sex Chromosomes

To understand the biology of a species, it is often crucial to be able to differentiate males and females. However, many species lack easily identifiable sexually dimorphic traits. In those that possess sex chromosomes, molecular sexing offers a good alternative, and molecular sexing assays can be developed through the comparison of male and female genomic sequences. However, in many nonmodel species, sex chromosomes are poorly differentiated, and identifying sex-linked sequences and developing sexing assays can be challenging. In this study, we highlight a simple transcriptome-based procedure for the detection of sex-linked markers suitable for the development of sexing assays that circumvents limitations of more commonly used approaches. We apply it to the spotted snow skink Carinascincus ocellatus, a viviparous lizard with homomorphic XY chromosomes that has environmentally induced sex reversal. With transcriptomes from three males and three females alone, we identify thousands of putative Y-linked sequences. We confirm linkage through alignment of assembled transcripts to a distantly related lizard genome and readily design multiple single locus polymerase chain reaction primers to sex C. ocellatus and related species. Our approach also facilitates valuable comparisons of sex determining systems on a broad taxonomic scale.

Lack of Dosage Balance and Incomplete Dosage Compensation in the ZZ/ZW Gila Monster (Heloderma suspectum) Revealed by De Novo Genome Assembly

Reptiles exhibit a variety of modes of sex determination, including both temperature-dependent and genetic mechanisms. Among those species with genetic sex determination, sex chromosomes of varying heterogamety (XX/XY and ZZ/ZW) have been observed with different degrees of differentiation. Karyotype studies have demonstrated that Gila monsters (Heloderma suspectum) have ZZ/ZW sex determination and this system is likely homologous to the ZZ/ZW system in the Komodo dragon (Varanus komodoensis), but little else is known about their sex chromosomes. Here, we report the assembly and analysis of the Gila monster genome. We generated a de novo draft genome assembly for a male using 10X Genomics technology. We further generated and analyzed short-read whole genome sequencing and whole transcriptome sequencing data for three males and three females. By comparing female and male genomic data, we identified four putative Z chromosome scaffolds. These putative Z chromosome scaffolds are homologous to Z-linked scaffolds identified in the Komodo dragon. Further, by analyzing RNAseq data, we observed evidence of incomplete dosage compensation between the Gila monster Z chromosome and autosomes and a lack of balance in Z-linked expression between the sexes. In particular, we observe lower expression of the Z in females (ZW) than males (ZZ) on a global basis, though we find evidence suggesting local gene-by-gene compensation. This pattern has been observed in most other ZZ/ZW systems studied to date and may represent a general pattern for female heterogamety in vertebrates.

Multiple Transitions between Y Chromosome and Autosome in Tago's Brown Frog Species Complex

Sex chromosome turnover is the transition between sex chromosomes and autosomes. Although many cases have been reported in poikilothermic vertebrates, their evolutionary causes and genetic mechanisms remain unclear. In this study, we report multiple transitions between the Y chromosome and autosome in the Japanese Tago's brown frog complex. Using chromosome banding and molecular analyses (sex-linked and autosomal single nucleotide polymorphisms, SNPs, from the nuclear genome), we investigated the frogs of geographic populations ranging from northern to southern Japan of two species, Rana tagoi and Rana sakuraii (2n = 26). Particularly, the Chiba populations of East Japan and Akita populations of North Japan in R. tagoi have been, for the first time, investigated here. As a result, we identified three different sex chromosomes, namely chromosomes 3, 7, and 13, in the populations of the two species. Furthermore, we found that the transition between the Y chromosome (chromosome 7) and autosome was repeated through hybridization between two or three different populations belonging to the two species, followed by restricted chromosome introgression. These dynamic sex chromosome turnovers represent the first such findings in vertebrates and imply that speciation associated with inter- or intraspecific hybridization plays an important role in sex chromosome turnover in frogs.

Sex-Linked Loci on the W Chromosome in the Multi-Ocellated Racerunner (Eremias multiocellata) Confirm Genetic Sex-Determination Stability in Lacertid Lizards

The multi-ocellated racerunner, Eremias multiocellata, was considered to have temperature-dependent sex determination (TSD), as its sex ratio can be influenced at different temperatures. However, such an observation contrasts with recent findings that suggest TSD is less common than previously thought. Here, a genotyping-by-sequencing (GBS) approach was employed to identify sex-linked markers in the E. multiocellata, for which the mechanism choice of TSD or GSD is still controversial. We preliminarily identified 119 sex-linked markers based on sex-associated sex-specific sequences, 97% of which indicated female heterogamety. After eliminating the false positives, 38 sex-linked markers were recognized, all of which showed the ZW/ZZ system. Then, eight of the novel markers were verified by PCR amplification from 15 populations of E. multiocellata, which support the GSD in E. multiocellata without geographic variation. To test the conservation of sex chromosome in Eremias, the eight markers were further cross-tested by PCR amplification in 10 individuals of the Mongolian racerunner (Eremias argus), two of which exhibited cross-utility. The novel sex-linked markers could be mapped on the W chromosome of the sand lizard (Lacerta agilis). Our finding that the sex-linked markers are shared in closely related species, along with a conserved synteny of the W chromosome, further supports the homology and conservation of sex chromosomes in the lacertid lizards.

Context-dependent thermolability of sex determination in a lacertid lizard with heteromorphic sex chromosomes

Developmental conditions can profoundly impact key life history traits of the individual. In cases where offspring sex is driven by developmental reaction norms, permanent changes to the phenotype can fundamentally alter life history trajectories. Sex determination mechanisms in reptiles are remarkably diverse, including well-characterised genetic and temperature-dependent sex determination. In rarer, but increasingly more commonly documented cases, sex can also be determined by a combination of the two, with temperature overriding the genetically determined sex. Thus, sex-by-temperature interactions is a mechanism that can be contextually labile, where reaction norms of sex against developmental environment might only be observable under certain conditions. We examine the effects of incubation temperature on hatchling sex in an oviparous lizard with clearly defined heteromorphic sex chromosomes presumed to determine sex solely on a genetic basis. We also test the repeatability of our results by replicating incubation experiments across 3 years. We show that warmer temperatures may override chromosomal sex and cause an overproduction of daughters. However, this effect was inconsistent among years, with high temperature only resulting in a daughter-significant bias in one year. Warm-incubated daughters were more efficient at converting yolk into tissue, which would allow for greater resource allocation to other fitness-related processes, such as growth. This suggests that thermolabile sex determination could be a trait under selection. More energy-efficient embryos also produced faster-growing offspring, suggesting that energy utilization patterns of the embryo were maintained into the juvenile stage, which could have important implications for the ontogenetic development and evolution of life histories.

Incomplete dosage balance and dosage compensation in the ZZ/ZW Gila monster (Heloderma suspectum) revealed by de novo genome assembly

Reptiles exhibit a variety of modes of sex determination, including both temperature-dependent and genetic mechanisms. Among those species with genetic sex determination, sex chromosomes of varying heterogamety (XX/XY and ZZ/ZW) have been observed with different degrees of differentiation. Karyotype studies have demonstrated that Gila monsters (Heloderma suspectum) have ZZ/ZW sex determination and this system is likely homologous to the ZZ/ZW system in the Komodo dragon (Varanus komodoensis), but little else is known about their sex chromosomes. Here, we report the assembly and analysis of the Gila monster genome. We generated a de novo draft genome assembly for a male using 10X Genomics technology. We further generated and analyzed short-read whole genome sequencing and whole transcriptome sequencing data for three males and three females. By comparing female and male genomic data, we identified four putative Z-chromosome scaffolds. These putative Z-chromosome scaffolds are homologous to Z-linked scaffolds identified in the Komodo dragon. Further, by analyzing RNAseq data, we observed evidence of incomplete dosage compensation between the Gila monster Z chromosome and autosomes and a lack of balance in Z-linked expression between the sexes. In particular, we observe lower expression of the Z in females (ZW) than males (ZZ) on a global basis, though we find evidence suggesting local gene-by-gene compensation. This pattern has been observed in most other ZZ/ZW systems studied to date and may represent a general pattern for female heterogamety in vertebrates.

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.

Sex-linked markers in an Australian frog Platyplectrum ornatum (Limnodynastidae) with a small genome and homomorphic sex chromosomes

Amphibians have highly diverse sex-determining modes leading to a notable interest in vertebrate sex determination and sex chromosome evolution. The identification of sex-determining systems in amphibians, however, is often difficult as a vast majority consist of homomorphic sex chromosomes making them hard to distinguish. In this study, we used Diversity Array Technology sequencing (DArTseq) to identify the sex-determining system in the ornate burrowing frog from Australia, Platyplectrum ornatum. We applied DArTseq to 44 individuals, 19 males and 25 females, collected from two locations to develop sex-linked markers. Unexpectedly, these 44 individuals were classified into two distinct population clusters based on our SNP analyses, 36 individuals in cluster 1, and 8 individuals in cluster 2. We then performed sex-linkage analyses separately in each cluster. We identified 35 sex-linked markers from cluster 1, which were all associated with maleness. Therefore, P. ornatum cluster 1 is utilising a male heterogametic (XX/XY) sex-determining system. On the other hand, we identified 210 sex-linked markers from cluster 2, of which 89 were male specific, i.e., identifying XX/XY sex determining system and 111 were female specific, i.e., identifying ZZ/ZW sex determining system, suggesting existence of either male or female heterogametic sex determining system in cluster 2. We also performed cytogenetic analyses in 1 male and 1 female from cluster 1; however, we did not detect any visible differentiation between the X and Y sex chromosomes. We also mapped sex-linked markers from the two clusters against the P. ornatum genome and our comparative analysis indicated that the sex chromosomes in both clusters shared homologies to chromosome 10 (autosome) of Rana temporaria and ZWY sex chromosome of Xenopus tropicalis. Our preliminary data suggest that it is plausible that the cluster 2 has a potential to be either male or female heterogamety in sex determination, requiring further investigation.

Highly rapid and diverse sex chromosome evolution in the Odorrana frog species complex

Sex chromosomes in poikilothermal vertebrates are characterized by rapid and diverse evolution at the species or population level. Our previous study revealed that the Taiwanese frog Odorrana swinhoana (2n = 26) has a unique system of multiple sex chromosomes created by three sequential translocations among chromosomes 1, 3, and 7. To reveal the evolutionary history of sex chromosomes in the Odorrana species complex, we first identified the original, homomorphic sex chromosomes, prior to the occurrence of translocations, in the ancestral-type population of O. swinhoana. Then, we extended the investigation to a closely related Japanese species, Odorrana utsunomiyaorum, which is distributed on two small islands. We used a high-throughput nuclear genomic approach to analyze single-nucleotide polymorphisms and identify the sex-linked markers. Those isolated from the O. swinhoana ancestral-type population were found to be aligned to chromosome 1 and showed male heterogamety. In contrast, almost all the sex-linked markers isolated from O. utsunomiyaorum were heterozygous in females and homozygous in males and were aligned to chromosome 9. Morphologically, we confirmed chromosome 9 to be heteromorphic in females, showing a ZZ-ZW sex determination system, in which the W chromosomes were heterochromatinized in a stripe pattern along the chromosome axis. These results indicated that after divergence of the two species, the ancestral homomorphic sex chromosome 1 underwent highly rapid and diverse evolution, i.e., sequential translocations with two autosomes in O. swinhoana, and turnover to chromosome 9 in O. utsunomiyaorum, with a transition from XY to ZW heterogamety and change to heteromorphy.

Genetic determination and JARID2 over-expression in a thermal incubation experiment in Casque-Headed Lizard

Non-avian reptiles, unlike mammals and birds, have undergone numerous sex determination changes. Casque-Headed Lizards have replaced the ancestral XY system shared across pleurodonts with a new pair of XY chromosomes. However, the evolutionary forces that triggered this transition have remained unclear. An interesting hypothesis suggests that species with intermediate states, with sex chromosomes but also thermal-induced sex reversal at specific incubation temperatures, could be more susceptible to sex determination turnovers. We contrasted genotypic data (presence/absence of the Y chromosome) against the histology of gonads of embryos from stages 35–37 incubated at various temperatures, including typical male-producing (26°C) and female-producing (32°C) temperatures. Our work apparently reports for the first time the histology of gonads, including morphological changes, from stages 35–37 of development in the family Corytophanidae. We also observed that all embryos developed hemipenes, suggesting sex-linked developmental heterochrony. We observed perfect concordance between genotype and phenotype at all temperatures. However, analysis of transcriptomic data from embryos incubated at 26°C and 32°C identified transcript variants of the chromatin modifiers JARID2 and KDM6B that have been linked to temperature-dependent sex determination in other reptiles. Our work tested the validity of a mixed sex determination system in the family Corytophanidae. We found that XY chromosomes are dominant; however, our work supports the hypothesis of a conserved transcriptional response to incubation temperatures across non-avian reptiles that could be a reminiscence of an ancestral sex determination system.

Sex reversal explains some, but not all, climate-mediated sex ratio variation within a viviparous reptile

Evolutionary transitions in sex-determining systems have occurred frequently yet understanding how they occur remains a major challenge. In reptiles, transitions from genetic to temperature-dependent sex determination can occur if the gene products that determine sex evolve thermal sensitivity, resulting in sex-reversed individuals. However, evidence of sex reversal is limited to oviparous reptiles. Here we used thermal experiments to test whether sex reversal is responsible for differences in sex determination in a viviparous reptile, Carinascincus ocellatus, a species with XY sex chromosomes and population-specific sex ratio response to temperature. We show that sex reversal is occurring and that its frequency is related to temperature. Sex reversal was unidirectional (phenotypic males with XX genotype) and observed in both high- and low-elevation populations. We propose that XX-biased genotypic sex ratios could produce either male- or female-biased phenotypic sex ratios as observed in low-elevation C. ocellatus under variable rates of XX sex reversal. We discuss reasons why sex reversal may not influence sex ratios at high elevation. Our results suggest that the mechanism responsible for evolutionary transitions from genotypic to temperature-dependent sex determination is more complex than can be explained by a single process such as sex reversal.

Labile sex chromosomes in the Australian freshwater fish family Percichthyidae

Sex-specific ecology has management implications, but rapid sex-chromosome turnover in fishes hinders sex-marker development for monomorphic species. We used annotated genomes and reduced-representation sequencing data for two Australian percichthyids, Macquarie perch Macquaria australasica and golden perch M. ambigua, and whole genome resequencing for 50 Macquarie perch of each sex, to identify sex-linked loci and develop an affordable sexing assay. In silico pool-seq tests of 1,492,004 Macquarie perch SNPs revealed that a 275-kb scaffold was enriched for gametologous loci. Within this scaffold, 22 loci were sex-linked in a predominantly XY system, with females being homozygous for the X-linked allele at all 22, and males having the Y-linked allele at >7. Seven XY-gametologous loci (all males, but no females, are heterozygous or homozygous for the male-specific allele) were within a 146-bp region. A PCR-RFLP sexing assay targeting one Y-linked SNP, tested in 66 known-sex Macquarie perch and two of each sex of three confamilial species, plus amplicon sequencing of 400 bp encompassing the 146-bp region, revealed that the few sex-linked positions differ between species and between Macquarie perch populations. This indicates sex-chromosome lability in Percichthyidae, supported by nonhomologous scaffolds containing sex-linked loci for Macquarie- and golden perches. The present resources facilitate genomic research in Percichthyidae, including formulation of hypotheses about candidate genes of interest such as transcription factor SOX1b that occurs in the 275-kb scaffold ~38 kb downstream of the 146-bp region containing seven XY-gametologous loci. Sex-linked markers will be useful for determining genetic sex in some populations and studying sex chromosome turnover.

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.

Viviparous Reptile Regarded to Have Temperature-Dependent Sex Determination Has Old XY Chromosomes

The water skinks Eulamprus tympanum and Eulamprus heatwolei show thermally induced sex determination where elevated temperatures give rise to male offspring. Paradoxically, Eulamprus species reproduce in temperatures of 12–15 °C making them outliers when compared with reptiles that use temperature as a cue for sex determination. Moreover, these two species are among the very few viviparous reptiles reported to have thermally induced sex determination. Thus, we tested whether these skinks possess undetected sex chromosomes with thermal override. We produced transcriptome and genome data for E. heatwolei. We found that E. heatwolei presents XY chromosomes that include 14 gametologs with regulatory functions. The Y chromosomal region is 79–116 Myr old and shared between water and spotted skinks. Our work provides clear evidence that climate could be useful to predict the type of sex determination systems in reptiles and it also indicates that viviparity is strictly associated with sex chromosomes.

Are Geckos Special in Sex Determination? Independently Evolved Differentiated ZZ/ZW Sex Chromosomes in Carphodactylid Geckos

Amniotes possess astonishing variability in sex determination ranging from environmental sex determination (ESD) to genotypic sex determination (GSD) with highly differentiated sex chromosomes. Geckos are one of the few amniote groups with substantial variability in sex determination. What makes them special in this respect? We hypothesized that the extraordinary variability of sex determination in geckos can be explained by two alternatives: 1) unusual lability of sex determination, predicting that the current GSD systems were recently formed and are prone to turnovers; and 2) independent transitions from the ancestral ESD to later stable GSD, which assumes that geckos possessed ancestrally ESD, but once sex chromosomes emerged, they remain stable in the long term. Here, based on genomic data, we document that the differentiated ZZ/ZW sex chromosomes evolved within carphodactylid geckos independently from other gekkotan lineages and remained stable in the genera Nephrurus, Underwoodisaurus, and Saltuarius for at least 15 Myr and potentially up to 45 Myr. These results together with evidence for the stability of sex chromosomes in other gekkotan lineages support more our second hypothesis suggesting that geckos do not dramatically differ from the evolutionary transitions in sex determination observed in the majority of the amniote lineages.

Temperature-Induced Sex Reversal in Reptiles: Prevalence, Discovery, and Evolutionary Implications

Sex reversal is the process by which an individual develops a phenotypic sex that is discordant with its chromosomal or genotypic sex. It occurs in many lineages of ectothermic vertebrates, such as fish, amphibians, and at least one agamid and one scincid reptile species. Sex reversal is usually triggered by an environmental cue that alters the genetically determined process of sexual differentiation, but it can also be caused by exposure to exogenous chemicals, hormones, or pollutants. Despite the occurrence of both temperature-dependent sex determination (TSD) and genetic sex determination (GSD) broadly among reptiles, only 2 species of squamates have thus far been demonstrated to possess sex reversal in nature (GSD with overriding thermal influence). The lack of species with unambiguously identified sex reversal is not necessarily a reflection of a low incidence of this trait among reptiles. Indeed, sex reversal may be relatively common in reptiles, but little is known of its prevalence, the mechanisms by which it occurs, or the consequences of sex reversal for species in the wild under a changing climate. In this review, we present a roadmap to the discovery of sex reversal in reptiles, outlining the various techniques that allow new occurrences of sex reversal to be identified, the molecular mechanisms that may be involved in sex reversal and how to identify them, and approaches for assessing the impacts of sex reversal in wild populations. We discuss the evolutionary implications of sex reversal and use the central bearded dragon (Pogona vitticeps) and the eastern three-lined skink (Bassiana duperreyi) as examples of how species with opposing patterns of sex reversal may be impacted differently by our rapidly changing climate. Ultimately, this review serves to highlight the importance of understanding sex reversal both in the laboratory and in wild populations and proposes practical solutions to foster future research.

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.

Pleistocene divergence in the absence of gene flow among populations of a viviparous reptile with intraspecific variation in sex determination

Polymorphisms can lead to genetic isolation if there is differential mating success among conspecifics divergent for a trait. Polymorphism for sex-determining system may fall into this category, given strong selection for the production of viable males and females and the low success of heterogametic hybrids when sex chromosomes differ (Haldane's rule). Here we investigated whether populations exhibiting polymorphism for sex determination are genetically isolated, using the viviparous snow skink Carinascincus ocellatus. While a comparatively high elevation population has genotypic sex determination, in a lower elevation population there is an additional temperature component to sex determination. Based on 11,107 SNP markers, these populations appear genetically isolated. "Isolation with Migration" analysis also suggests these populations diverged in the absence of gene flow, across a period encompassing multiple Pleistocene glaciations and likely greater geographic proximity of populations. However, further experiments are required to establish whether genetic isolation may be a cause or consequence of differences in sex determination. Given the influence of temperature on sex in one lineage, we also discuss the implications for the persistence of this polymorphism under climate change.

High elevation increases the risk of Y chromosome loss in Alpine skink populations with sex reversal

The view that has genotypic sex determination and environmental sex determination as mutually exclusive states in fishes and reptiles has been contradicted by the discovery that chromosomal sex and environmental influences can co-exist within the same species, hinting at a continuum of intermediate states. Systems where genes and the environment interact to determine sex present the opportunity for sex reversal to occur, where the phenotypic sex is the opposite of that predicted by their sex chromosome complement. The skink Bassiana duperreyi has XX/XY sex chromosomes with sex reversal of the XX genotype to a male phenotype, in laboratory experiments, and in field nests, in response to exposure to cold incubation temperatures. Here we studied the frequency of sex reversal in adult populations of B. duperreyi in response to climatic variation, using elevation as a surrogate for environmental temperatures. We demonstrate sex reversal in the wild for the first time in adults of a reptile species with XX/XY sex determination. The highest frequency of sex reversal occurred at the highest coolest elevation location, Mt Ginini (18.46%) and decreased in frequency to zero with decreasing elevation. We model the impact of this under Fisher's frequency-dependent selection to show that, at the highest elevations, populations risk the loss of the Y chromosome and a transition to temperature-dependent sex determination. This study contributes to our understanding of the risks of extinction from climate change in species subject to sex reversal by temperature, and will provide focus for future research to test on-the-ground management strategies to mitigate the effects of climate in local populations.

Do the Offspring of Sex Reversals Have Higher Sensitivity to Environmental Perturbations?

Sex determination systems in vertebrates vary along a continuum from genetic (GSD) to environmental sex determination (ESD). Individuals that show a sexual phenotype opposite to their genotypic sex are called sex reversals. Aside from genetic elements, temperature, sex steroids, and exogenous chemicals are common factors triggering sex reversal, a phenomenon that may occur even in strict GSD species. In this paper, we review the literature on instances of sex reversal in fish, amphibians, reptiles, birds, and mammals. We focus on the offspring of sex-reversed parents in the instances that they can be produced, and show that in all cases studied the offspring of these sex-reversed parents exhibit a higher sensitivity to environmental perturbations than the offspring of non-sex-reversed parents. We suggest that the inheritance of this sensitivity, aside from possible genetic factors, is likely to be mediated by epigenetic mechanisms such as DNA methylation, since these mechanisms are responsive to environmental cues, and epigenetic modifications can be transmitted to the subsequent generations. Species with a chromosomal GSD system with environmental sensitivity and availability of genetic sex markers should be employed to further test whether offspring of sex-reversed parents have greater sensitivity to environmental perturbations. Future studies could also benefit from detailed whole-genome data in order to elucidate the underlying molecular mechanisms. Finally, we discuss the consequences of such higher sensitivity in the context of global climate change.

With or Without W? Molecular and Cytogenetic Markers are Not Sufficient for Identification of Environmentally-Induced Sex Reversal in the Bearded Dragon

Transitions from environmental sex determination (ESD) to genotypic sex determination (GSD) require an intermediate step of sex reversal, i.e., the production of individuals with a mismatch between the ancestral genotypic and the phenotypic sex. Among amniotes, the sole well-documented transition in this direction was shown in the laboratory in the central bearded dragon, Pogona vitticeps, where very high incubation temperatures led to the production of females with the male-typical (ZZ) genotype. These sex-reversed females then produced offspring whose sex depended on the incubation temperature. Sex-reversed animals identified by molecular and cytogenetic markers were also reported in the field, and their increasing incidence was speculated as a climate warming-driven transition in sex determination. We show that the molecular and cytogenetic markers normally sex-linked in P. vitticeps are also sex-linked in P. henrylawsoni and P. minor, which points to quite ancient sex chromosomes in this lineage. Nevertheless, we demonstrate, based on a crossing experiment with a male bearded dragon who possesses a mismatch between phenotypic sex and genotype, that the used cytogenetic and molecular markers might not be reliable for the identification of sex reversal. Sex reversal should not be considered as the only mechanism causing a mismatch between genetic sex-linked markers and phenotypic sex, which can emerge also by other processes, here most likely by a rare recombination between regions of sex chromosomes which are normally sex-linked. We warn that sex-linked, even apparently for a long evolutionary time, and sex-specific molecular and cytogenetic markers are not a reliable tool for the identification of sex-reversed individuals in a population and that sex reversal has to be verified by other approaches, particularly by observation of the sex ratio of the progeny.

Differences in Homomorphic Sex Chromosomes Are Associated with Population Divergence in Sex Determination in Carinascincus ocellatus (Scincidae: Lygosominae)

Sex determination directs development as male or female in sexually reproducing organisms. Evolutionary transitions in sex determination have occurred frequently, suggesting simple mechanisms behind the transitions, yet their detail remains elusive. Here we explore the links between mechanisms of transitions in sex determination and sex chromosome evolution at both recent and deeper temporal scales (<1 Myr; ~79 Myr). We studied a rare example of a species with intraspecific variation in sex determination, Carinascincus ocellatus, and a relative, Liopholis whitii, using c-banding and mapping of repeat motifs and a custom Y chromosome probe set to identify the sex chromosomes. We identified both unique and conserved regions of the Y chromosome among C. ocellatus populations differing in sex determination. There was no evidence for homology of sex chromosomes between C. ocellatus and L. whitii, suggesting independent evolutionary origins. We discuss sex chromosome homology between members of the subfamily Lygosominae and propose links between sex chromosome evolution, sex determination transitions, and karyotype evolution.

Poorly differentiated XX/XY sex chromosomes are widely shared across skink radiation

Differentiated sex chromosomes are believed to be evolutionarily stable, while poorly differentiated sex chromosomes are considered to be prone to turnovers. With around 1700 currently known species forming ca 15% of reptile species diversity, skinks (family Scincidae) are a very diverse group of squamates known for their large ecological and morphological variability. Skinks generally have poorly differentiated and cytogenetically indistinguishable sex chromosomes, and their sex determination was suggested to be highly variable. Here, we determined X-linked genes in the common sandfish (Scincus scincus) and demonstrate that skinks have shared the same homologous XX/XY sex chromosomes across their wide phylogenetic spectrum for at least 85 million years, approaching the age of the highly differentiated ZZ/ZW sex chromosomes of birds and advanced snakes. Skinks thus demonstrate that even poorly differentiated sex chromosomes can be evolutionarily stable. The conservation of sex chromosomes across skinks allows us to introduce the first molecular sexing method widely applicable in this group.

An Investigation of ZZ/ZW and XX/XY Sex Determination Systems in North African Catfish (Clarias gariepinus, )

An investigation of sex-specific loci may provide important insights into fish sex determination strategies. This may be useful for biotechnological purposes, for example, to produce all-male or all-female fish for commercial breeding. The North African catfish species, Clarias gariepinus, has been widely adopted for aquaculture because its superior growth and disease resistance render the species suitable for hybridization with other catfish to improve the productivity and quality of fish meat. This species has either a ZZ/ZW or XX/XY sex determination system. Here, we investigate and characterize these systems using high-throughput genome complexity reduction sequencing as Diversity Arrays Technology. This approach was effective in identifying moderately sex-linked loci with both single-nucleotide polymorphisms (SNPs) and restriction fragment presence/absence (PA) markers in 30 perfectly sexed individuals of C. gariepinus. However, SNPs based markers were not found in this study. In total, 41 loci met the criteria for being moderately male-linked (with male vs. female ratios 80:20 and 70:30), while 25 loci were found to be moderately linked to female sex. No strictly male- or female-linked loci were detected. Seven moderately male-linked loci were partially homologous to some classes of transposable elements and three moderately male-linked loci were partially homologous to functional genes. Our data showed that the male heterogametic XX/XY sex determination system should co-exist with the ZZ/ZW system in C. gariepinus. Our finding of the co-existence of XX/XY and ZZ/ZW systems can be applied to benefit commercial breeding of this species in Thailand. This approach using moderately sex-linked loci provides a solid baseline for revealing sex determination mechanisms and identify potential sex determination regions in catfish, allowing further investigation of genetic improvements in breeding programs.

Genome-wide SNP analysis of Siamese cobra (Naja kaouthia) reveals the molecular basis of transitions between Z and W sex chromosomes and supports the presence of an ancestral super-sex chromosome in amniotes

Elucidation of the process of sex chromosome differentiation is necessary to understand the dynamics of evolutionary mechanisms in organisms. The W sex chromosome of the Siamese cobra (Naja kaouthia) contains a large number of repeats and shares amniote sex chromosomal linkages. Diversity Arrays Technology provides an effective approach to identify sex-specific loci that are epoch-making, to understand the dynamics of molecular transitions between the Z and W sex chromosomes in a snake lineage. From a total of 543 sex-specific loci, 90 showed partial homology with sex chromosomes of several amniotes and 89 loci were homologous to transposable elements. Two loci were confirmed as W-specific nucleotides after PCR amplification. These loci might result from a sex chromosome differentiation process and involve putative sex-determination regions in the Siamese cobra. Sex-specific loci shared linkage homologies among amniote sex chromosomes, supporting an ancestral super-sex chromosome.

Degrees of change: between and within population variation in thermal reaction norms of phenology in a viviparous lizard

As the earth warms, populations will be faced with novel environments to which they may not be adapted. In the short term, populations can be buffered against the negative effects, or maximize the beneficial effects, of such environmental change via phenotypic plasticity and, in the longer term, via adaptive evolution. However, the extent and direction of these population-level responses will be dependent on the degree to which responses vary among the individuals within them (i.e., within population variation in plasticity), which is, itself, likely to vary among populations. Despite this, we have estimates of among-individual variation in plastic responses across multiple populations for only a few systems. This lack of data limits our ability to predict the consequences of environmental change for population and species persistence accurately. Here, we utilized a 16-yr data set from climatically distinct populations of the viviparous skink Niveoscincus ocellatus tracking over 1,200 litters from more than 600 females from each population to examine inter- and intrapopulation variability in the response of parturition date to environmental temperature. We found that these populations share a common population-mean reaction norm but differ in the degree to which reaction norms vary among individuals. These results suggest that even where populations share a common mean-level response, we cannot assume that they will be affected similarly by altered environmental conditions. If we are to assess how changing climates will impact species and populations accurately, we require estimates of how plastic responses vary both among and within populations.

Sex determination systems in reptiles are related to ambient temperature but not to the level of climatic fluctuation

Background

Vertebrates exhibit diverse sex determination systems and reptiles stand out by having highly variable sex determinations that include temperature-dependent and genotypic sex determination (TSD and GSD, respectively). Theory predicts that populations living in either highly variable or cold climatic conditions should evolve genotypic sex determination to buffer the populations from extreme sex ratios, yet these fundamental predictions have not been tested across a wide range of taxa.

Results

Here, we use phylogenetic analyses of 213 reptile species representing 38 families (TSD = 101 species, GSD = 112 species) and climatic data to compare breeding environments between reptiles with GSD versus TSD. We show that GSD and TSD are confronted with the same level of climatic fluctuation during breeding seasons. However, TSD reptiles are significantly associated with warmer climates. We found a strong selection on the breeding season length that minimises exposure to cold and fluctuating climate. Phylogenetic path analyses comparing competing evolutionary hypotheses support that transitions in sex determination systems influenced the ambient temperature at which the species reproduces and nests. In turn, this interaction affects other variables such as the duration of the breeding season and life-history traits.

Conclusions

Taken together, our results challenge long-standing hypotheses about the association between sex determination and climate variability. We also show that ambient temperature is important during breeding seasons and it helps explain the effects of sex determination systems on the geographic distribution of extant reptile species.

Reproductive phenotype predicts adult bite-force performance in sex-reversed dragons (Pogona vitticeps)

Sex-related differences in morphology and behavior are well documented, but the relative contributions of genes and environment to these traits are less well understood. Species that undergo sex reversal, such as the central bearded dragon (Pogona vitticeps), offer an opportunity to better understand sexually dimorphic traits because sexual phenotypes can exist on different chromosomal backgrounds. Reproductively female dragons with a discordant sex chromosome complement (sex reversed), at least as juveniles, exhibit traits in common with males (e.g., longer tails and greater boldness). However, the impact of sex reversal on sexually dimorphic traits in adult dragons is unknown. Here, we investigate the effect of sex reversal on bite-force performance, which may be important in resource acquisition (e.g., mates and/or food). We measured body size, head size, and bite force of the three sexual phenotypes in a colony of captive animals. Among adults, we found that males (ZZm) bite more forcefully than either chromosomally concordant females (ZWf) or sex-reversed females (ZZf), and this difference is associated with having relatively larger head dimensions. Therefore, adult sex-reversed females, despite apparently exhibiting male traits as juveniles, do not develop the larger head and enhanced bite force of adult male bearded dragons. This pattern is further illustrated in the full sample by a lack of positive allometry of bite force in sex-reversed females that is observed in males. The results reveal a close association between reproductive phenotype and bite force performance, regardless of sex chromosome complement.

Australian lizards are outstanding models for reproductive biology research

Australian lizards are a diverse group distributed across the continent and inhabiting a wide range of environments. Together, they exhibit a remarkable diversity of reproductive morphologies, physiologies, and behaviours that is broadly representative of vertebrates in general. Many reproductive traits exhibited by Australian lizards have evolved independently in multiple lizard lineages, including sociality, complex signalling and mating systems, viviparity, and temperature-dependent sex determination. Australian lizards are thus outstanding model organisms for testing hypotheses about how reproductive traits function and evolve, and they provide an important basis of comparison with other animals that exhibit similar traits. We review how research on Australian lizard reproduction has contributed to answering broader evolutionary and ecological questions that apply to animals in general. We focus on reproductive traits, processes, and strategies that are important areas of current research, including behaviours and signalling involved in courtship; mechanisms involved in mating, egg production, and sperm competition; nesting and gestation; sex determination; and finally, birth in viviparous species. We use our review to identify important questions that emerge from an understanding of this body of research when considered holistically. Finally, we identify additional research questions within each topic that Australian lizards are well suited for reproductive biologists to address.

Identifying sex-linked markers in Litoria aurea: a novel approach to understanding sex chromosome evolution in an amphibian

Few taxa exhibit the variability of sex-determining modes as amphibians. However, due to the presence of homomorphic sex chromosomes in many species, this phenomenon has been difficult to study. The Australian frog, Litoria aurea, has been relatively well studied over the past 20 years due to widespread declines largely attributable to chytrid fungus. However, it has been subject to few molecular studies and its mode of sex determination remained unknown. We applied DArTseq™ to develop sex-linked single nucleotide polymorphisms (SNPs) and restriction fragment presence/absence (PA) markers in 44 phenotypically sexed L. aurea individuals from the Molonglo River in NSW, Australia. We conclusively identified a male heterogametic (XX-XY) sex determination mode in this species, identifying 11 perfectly sex-linked SNP and six strongly sex-linked PA markers. We identified a further 47 moderately sex-linked SNP loci, likely serving as evidence indicative of XY recombination. Furthermore, within these 47 loci, a group of nine males were found to have a feminised Y chromosome that significantly differed to all other males. We postulate ancestral sex-reversal as a means for the evolution of this now pseudoautosomal region on the Y chromosome. Our findings present new evidence for the ‘fountain of youth’ hypothesis for the retention of homomorphic sex chromosomes in amphibians and describe a novel approach for the study of sex chromosome evolution in amphibia.

Temperature and telomeres: thermal treatment influences telomere dynamics through a complex interplay of cellular processes in a cold-climate skink

Telomere dynamics vary fundamentally between endothermic populations and species as a result of differences in life history, yet we know little about these patterns in ectotherms. In ectotherms, the relationships between climate, metabolism and life history suggest that telomere attrition should be higher at relatively high environmental temperatures compared to relatively low environmental temperatures, but these effects may vary between populations due to local adaptation. To address this hypothesis, we sampled reactive oxygen species (ROS) and telomere length of lizards from warm lowland and cool highland populations of a climatically widespread lizard species that we exposed to hot or cold basking treatments. The hot treatment increased relative telomere length compared to the cold treatment independent of climatic origin or ROS levels. Lizards from the cool highland region had lower ROS levels than those from the warm lowland region. Within the highland lizards, ROS increased more in the cold basking treatment than the hot basking treatment. These results are in the opposite direction to those predicted, suggesting that the relationships between temperature, metabolism, ROS and telomere dynamics are not straightforward. Future work incorporating detailed understanding of the thermal reaction norms of these and other linked traits is needed to fully understand these processes.

Little evidence for switches to environmental sex determination and turnover of sex chromosomes in lacertid lizards

Amniotes possess variability in sex determination, from environmental sex determination (ESD), where no sex chromosomes are present, to genotypic sex determination (GSD) with highly differentiated sex chromosomes. Some evolutionary scenarios postulate high stability of differentiated sex chromosomes and rare transitions from GSD to ESD. However, sex chromosome turnovers and two independent transitions from highly differentiated ZZ/ZW sex chromosomes to ESD were previously reported in the lacertid lizards. Here, we examined the homology of sex chromosomes in the wide phylogenetic spectrum of lacertids and their outgroups by comparing gene copy numbers between sexes in genes previously found to be Z-specific in some lacertids. Our current sampling covers 45 species from 26 genera including lineages supposed to possess a derived sex determining systems. We found that all tested lacertids share homologous differentiated ZZ/ZW sex chromosomes, which were present already in their common ancestor living around 85 million years ago. These differentiated sex chromosomes are not present in amphisbaenians and teiid lizards, the close relatives of lacertids. Our study demonstrates how inaccuracies in data can influence the outcome of phylogenetic reconstructions of evolution of sex determination, in this case they overestimated the number of shifts from GSD to ESD and the rate in turnovers of sex chromosomes.

Developmental asynchrony and antagonism of sex determination pathways in a lizard with temperature-induced sex reversal

Vertebrate sex differentiation follows a conserved suite of developmental events: the bipotential gonads differentiate and shortly thereafter sex specific traits become dimorphic. However, this may not apply to squamates, a diverse vertebrate lineage comprising of many species with thermosensitive sexual development. Of the three species with data on the relative timing of gonad differentiation and genital dimorphism, the females of two (Niveoscincus ocellatus and Barisia imbricata) exhibit a phase of temporary pseudohermaphroditism or TPH (gonads have differentiated well before genital dimorphism). We report a third example of TPH in Pogona vitticeps, an agamid with temperature-induced male to female sex reversal. These findings suggest that for female squamates, genital and gonad development may not be closely synchronised, so that TPH may be common. We further observed a high frequency of ovotestes, a usually rare gonadal phenotype characterised by a mix of male and female structures, exclusively associated with temperature-induced sex reversal. We propose that ovotestes are evidence of a period of antagonism between male and female sex-determining pathways during sex reversal. Female sexual development in squamates is considerably more complex than has been appreciated, providing numerous avenues for future exploration of the genetic and hormonal cues that govern sexual development.

Plastic rates of development and the effect of thermal extremes on offspring fitness in a cold-climate viviparous lizard

Populations at the climatic margins of a species' distribution can be exposed to conditions that cause developmental stress, resulting in developmental abnormalities. Even within the thermal range of normal development, phenotypes often vary with developmental temperature (i.e., thermal phenotypic plasticity). These effects can have significant consequences for organismal fitness and, thus, population persistence. Reptiles, as ectotherms, are particularly vulnerable to thermal effects on development and are, therefore, considered to be at comparatively high risk from changing climates. Understanding the extent and direction of thermal effects on phenotypes and their fitness consequences is crucial if we are to make meaningful predictions of how populations and species will respond as climates warm. Here, we experimentally manipulated the thermal conditions experienced by females from a high-altitude, cold-adapted population of the viviparous skink, Niveoscincus ocellatus, to examine the consequences of thermal conditions at the margins of this population's normal temperature range. We found strong effects of thermal conditions on the development of key phenotypic traits that have implications for fitness. Specifically, we found that offspring born earlier as a result of high temperatures during gestation had increased growth over the first winter of life, but there was no effect on offspring survival, nor was there an effect of developmental temperature on the incidence of developmental abnormalities. Combined, our results suggest that advancing birth dates that result from warming climates may have positive effects in this population via increased growth.