Mutations in the testis-specific enhancer of SOX9 in the SRY independent sex-determining mechanism in the genus Tokudaia

Identification of a New Enhancer That Promotes Sox9 Expression by a Comparative Analysis of Mouse and Sry-Deficient Amami Spiny Rat

INTRODUCTION

Testis differentiation is initiated by the SRY gene on the Y chromosome in mammalian species. However, the Amami spiny rat, Tokudaia osimensis, lacks both the Y chromosome and the Sry gene and acquired a unique Sox9 regulatory mechanism via a male-specific duplication upstream of Sox9, without Sry. In general mammalian species, the SRY protein binds to a testis-specific enhancer to promote SOX9 gene expression. Several enhancers located upstream of Sox9/SOX9 have been reported in mice and humans. In particular, the binding of SRY to the highly conserved enhancer Enh13 is thought to be a common mechanism underlying testis differentiation and sex determination in mammals.

METHODS

Sequences of T. osimensis homologues of three Sox9 enhancers that were previously reported in mice, Enh8, Enh14, and Enh13, were determined. We performed in vitro assays to confirm enhancer activity involved in Sox9 regulation in T. osimensis.

RESULTS

T. osimensis Enh13 showed enhancer activity when co-transfected with NR5A1 and SOX9. Mouse Enh13 was activated by NR5A1 and SRY; however, T. osimensis Enh13 did not respond to SRY, even though the binding sites of SRY and NR5A1 were conserved. To identify the key sequence that is present in mouse but absent from T. osimensis, we performed reporter gene assays using vectors in which partial sequences of T. osimensis Enh13 were replaced with mouse sequences. For T. osimensis Enh13 in which the second half (approximately 430 bp) was replaced with the corresponding mouse sequence, activity in response to NR5A1 and SRY was recovered. Further, reporter assays revealed that multiple regions in the second half of the mouse Enh13 sequence are required for the response to NR5A1 and SRY. The latter 49 bp was particularly important and contained four binding sites for three transcription factors, POU2F1, HOXA3, and GATA1.

CONCLUSION

We showed that there are unknown sequences responsible for the interaction between NR5A1 and SRY and mEnh13 based on comparative analyses of Sry-dependent and Sry-independent species. Our comparative analyses revealed new molecular mechanisms underlying mammalian sex determination.

Chromosomal-level assembly of Tokudaia osimensis, Tokudaia tokunoshimensis, and Tokudaia muenninki genomes

Herein, we present the first high-quality long-read-based chromosome-level genome assemblies and gene annotations of the genomes of three endangered Tokudaia species: Tokudaia osimensis, Tokudaia tokunoshimensis, and Tokudaia muenninki. These species, which are endemic to different islands of the Ryukyu Islands, Japan, exhibited unique karyotypes and sex chromosomal characteristics. The genome assemblies generated using PacBio, Illumina, and Hi-C sequence data consisted of 13 (corresponded to 12 autosomes and one X chromosome), 23 (corresponded to 22 autosomes and one X chromosome), and 23 (corresponded to 21 autosomes and the neo- and ancestral X regions) chromosome-level scaffolds that contained 2,445, 2,477, and 2,661 Mbp of sequence data, respectively. Annotations of protein-coding genes were performed using RNA-Seq-based, homology-based, and Ab initio methods. BUSCO completeness values for every species exceeded 96% for genomes and 98% for genes. These data can be an important resource for contributing to our understanding of species genomes resulting from allopatric speciation and provide insights into mammalian sex-determination mechanisms and sex chromosome evolution.

The Neo-X Does Not Form a Barr Body but Shows a Slightly Condensed Structure in the Okinawa Spiny Rat (Tokudaia muenninki)

X chromosome inactivation (XCI) is an essential mechanism for gene dosage compensation between male and female cells in mammals. The Okinawa spiny rat (Tokudaia muenninki) is a native rodent in Japan with XX/XY sex chromosomes, like most mammals; however, the X chromosome has acquired a neo-X region (Xp) by fusion with an autosome. We previously reported that dosage compensation has not yet evolved in the neo-X region; however, X-inactive-specific transcript (Xist) RNA (long non-coding RNA required for the initiation of XCI) is partially localized in the region. Here, we show that the neo-X region represents an early chromosomal state in the acquisition of XCI by analyses of heterochromatin and Barr body formation. We found no evidence for heterochromatin formation in the neo-X region by R-banding by acridine orange (RBA) assays and immunostaining of H3K27me3. Double-immunostaining of H3K27me3 and HP1, a component of the Barr body, revealed that the entire ancestral X chromosome region (Xq) showed a bipartite folded structure. By contrast, HP1 was not localized to the neo-X region. However, BAC-FISH revealed that the signals of genes on the neo-X region of the inactive X chromosome were concentrated in a narrow region. These findings indicated that although the neo-X region of the inactive X chromosome does not form a complete Barr body structure (e.g., it lacks HP1), it forms a slightly condensed structure. These findings combined with the previously reported partial binding of Xist RNA suggest that the neo-X region exhibits incomplete inactivation. This may represent an early chromosomal state in the acquisition of the XCI mechanism.

Turnover of mammal sex chromosomes in the Sry-deficient Amami spiny rat is due to male-specific upregulation of Sox9

Mammalian sex chromosomes are highly conserved, and sex is determined by SRY on the Y chromosome. Two exceptional rodent groups in which some species lack a Y chromosome and Sry offer insights into how novel sex genes can arise and replace Sry, leading to sex chromosome turnover. However, intensive study over three decades has failed to reveal the identity of novel sex genes in either of these lineages. We here report our discovery of a male-specific duplication of an enhancer of Sox9 in the Amami spiny rat Tokudaia osimensis, in which males and females have only a single X chromosome (XO/XO) and the Y chromosome and Sry are completely lost. We performed a comprehensive survey to detect sex-specific genomic regions in the spiny rat. Sex-related genomic differences were limited to a male-specific duplication of a 17-kb unit located 430 kb upstream of Sox9 on an autosome. Hi-C analysis using male spiny rat cells showed the duplicated region has potential chromatin interaction with Sox9. The duplicated unit harbored a 1,262-bp element homologous to mouse enhancer 14 (Enh14), a candidate Sox9 enhancer that is functionally redundant in mice. Transgenic reporter mice showed that the spiny rat Enh14 can function as an embryonic testis enhancer in mice. Embryonic gonads of XX mice in which Enh14 was replaced by the duplicated spiny rat Enh14 showed increased Sox9 expression and decreased Foxl2 expression. We propose that male-specific duplication of this Sox9 enhancer substituted for Sry function, defining a novel Y chromosome in the spiny rat.

Decoding sex: Elucidating sex determination and how high-quality genome assemblies are untangling the evolutionary dynamics of sex chromosomes

Sexual reproduction is a diverse and widespread process. In gonochoristic species, the differentiation of sexes occurs through diverse mechanisms, influenced by environmental and genetic factors. In most vertebrates, a master-switch gene is responsible for triggering a sex determination network. However, only a few genes have acquired master-switch functions, and this process is associated with the evolution of sex-chromosomes, which have a significant influence in evolution. Additionally, their highly repetitive regions impose challenges for high-quality sequencing, even using high-throughput, state-of-the-art techniques. Here, we review the mechanisms involved in sex determination and their role in the evolution of species, particularly vertebrates, focusing on sex chromosomes and the challenges involved in sequencing these genomic elements. We also address the improvements provided by the growth of sequencing projects, by generating a massive number of near-gapless, telomere-to-telomere, chromosome-level, phased assemblies, increasing the number and quality of sex-chromosome sequences available for further studies.

Unusual Mammalian Sex Determination Systems: A Cabinet of Curiosities

Therian mammals have among the oldest and most conserved sex-determining systems known to date. Any deviation from the standard XX/XY mammalian sex chromosome constitution usually leads to sterility or poor fertility, due to the high differentiation and specialization of the X and Y chromosomes. Nevertheless, a handful of rodents harbor so-called unusual sex-determining systems. While in some species, fertile XY females are found, some others have completely lost their Y chromosome. These atypical species have fascinated researchers for over 60 years, and constitute unique natural models for the study of fundamental processes involved in sex determination in mammals and vertebrates. In this article, we review current knowledge of these species, discuss their similarities and differences, and attempt to expose how the study of their exceptional sex-determining systems can further our understanding of general processes involved in sex chromosome and sex determination evolution.

Regulation of the Sox3 Gene in an X0/X0 Mammal without Sry, the Amami Spiny Rat, Tokudaia osimensis

Two species of spiny rats, Tokudaia osimensis and Tokudaia tokunoshimensis, show an X0/X0 sex chromosome constitution due to the lack of a Y chromosome. The Sry gene has been completely lost from the genome of these species. We hypothesized that Sox3, which is thought to be originally a homologue of Sry, could function in sex determination in these animals in the absence of Sry. Sox3 was localized in a region of the X chromosome in T. osimensis homologous to mouse. A similar testis- and ovary-specific pattern of expression was observed in mouse and T. osimensis. Although the sequence of the Sox3 gene and its promoter are highly conserved, a 13-bp deletion was specifically found in the promoter region of the 2 spiny rat species. Reporter gene assays were performed to examine the effect of the 13-bp deletion in the promoter region on Sox3 regulation. Although an approximately 60% decrease in activity was observed using the Tokudaia promoters with the 13-bp deletion, the activity was recovered using a mutated promoter in which the deletion was filled with mouse sequence. To evaluate whether SOX3 could regulate Sox9 expression, a reporter gene assay was carried out using testis-specific enhancer of Sox9 core (TESCO). Co-transfection with a combination of mouse SF1 and mouse SOX3 or T. osimensis SOX3 resulted in a greater than 2-fold increase in activity of mouse and T. osimensis TESCO. These results support the idea that the function of SOX3 as a transcription factor, as has been reported in mice and humans, is conserved in T. osimensis. Therefore, we conclude that the Sox3 gene has no function in sex determination in Sry-lacking Tokudaia species.

Spiny rat SRY lacks a long Q-rich domain and is not stable in transgenic mice

BACKGROUND

Although Tokudaia muenninki has multiple extra copies of the Sry gene on the Y chromosome, loss of function of these sequences is indicated. To examine the Sry gene function for sex determining in T. muenninki, we screened a BAC library and identified a clone (SRY26) containing complete SRY coding and promoter sequences.

RESULTS

SRY26 showed high identity to mouse and rat SRY. In an in vitro reporter gene assay, SRY26 was unable to activate testis-specific enhancer of Sox9. Four lines of BAC transgenic mice carrying SRY26 were generated. Although the embryonic gonads of XX transgenic mice displayed sufficient expression levels of SRY26 mRNA, these mice exhibited normal female phenotypes in the external and internal genitalia, and up-regulation of Sox9 was not observed. Expression of the SRY26 protein was confirmed in primate-derived COS7 cells transfected with a SRY26 expression vector. However, the SRY26 protein was not expressed in the gonads of BAC transgenic mice.

CONCLUSIONS

Overall, these results support a previous study demonstrated a long Q-rich domain plays essential roles in protein stabilization in mice. Therefore, the original aim of this study, to examine the function of the Sry gene of this species, was not achieved by creating TG mice.

Chromosomal Evolution in Mole Voles Ellobius (Cricetidae, Rodentia): Bizarre Sex Chromosomes, Variable Autosomes and Meiosis

This study reports on extensive experimental material covering more than 30 years of studying the genetics of mole voles. Sex chromosomes of Ellobius demonstrate an extraordinary case of mammalian sex chromosomes evolution. Five species of mole voles own three types of sex chromosomes; typical for placentals: XY♂/XX♀; and atypical X0♂/X0♀; or XX♂/XX♀. Mechanisms of sex determination in all Ellobius species remain enigmatic. It was supposed that the Y chromosome was lost twice and independently in subgenera Bramus and Ellobius. Previous to the Y being lost, the X chromosome in distinct species obtained some parts of the Y chromosome, with or without Sry, and accumulated one or several copies of the Eif2s3y gene. Along with enormous variations of sex chromosomes, genes of sex determination pathway and autosomes, and five mole vole species demonstrate ability to establish different meiotic mechanisms, which stabilize their genetic systems and make it possible to overcome the evolutionary deadlocks.

Unique XCI evolution in Tokudaia: initial XCI of the neo-X chromosome in Tokudaia muenninki and function loss of XIST in Tokudaia osimensis

X chromosome inactivation (XCI) is an essential mechanism to compensate gene dosage in mammals. Here, we show that XCI has evolved differently in two species of the genus Tokudaia. The Amami spiny rat, Tokudaia osimensis, has a single X chromosome in males and females (XO/XO). By contrast, the Okinawa spiny rat, Tokudaia muenninki, has XX/XY sex chromosomes like most mammals, although the X chromosome has acquired a neo-X region by fusion with an autosome. BAC clones containing the XIST gene, which produces the long non-coding RNA XIST required for XCI, were obtained by screening of T. osimensis and T. muenninki BAC libraries. Each clone was mapped to the homologous region of the X inactivation center in the X chromosome of the two species by BAC-FISH. XIST RNAs were expressed in T. muenninki females, whereas no expression was observed in T. osimensis. The sequence of the XIST RNA was compared with that of mouse, showing that the XIST gene is highly conserved in T. muenninki. XIST RNAs were localized to the ancestral X region (Xq), to the heterochromatic region (pericentromeric region), and partially to the neo-X region (Xp). The hybridization pattern correlated with LINE-1 accumulation in Xq but not in Xp. Dosage of genes located on the neo-X chromosome was not compensated, suggesting that the neo-X region is in an early state of XCI. By contrast, many mutations were observed in the XIST gene of T. osimensis, indicating its loss of function in the XO/XO species.

Reduced Activity of SRY and its Target Enhancer Sox9-TESCO in a Mouse Species with X*Y Sex Reversal

In most eutherian mammals, sex determination is governed by the Y-linked gene Sry, but in African pygmy mice Mus minutoides, Sry action is overridden by a variant X chromosome (X*), yielding X*Y females. We hypothesized that X*Y sex reversal may be underpinned not only by neomorphic X chromosome functionality, but also by a compromised Sry pathway. Here, we show that neither M. minutoides SRY nor its target, the Sox9-TESCO enhancer, had appreciable transcriptional activity in in vitro assays, correlating with sequence degradation compared to Mus musculus counterparts. However, M. minutoides SRY activated its cognate TESCO to a moderate degree, and can clearly engage the male pathway in M. minutoides in the wild, indicating that SRY and TESCO may have co-evolved in M. minutoides to retain function above a threshold level. We suggest that weakening of the SRY/TESCO nexus may have facilitated the rise and spread of a variant X* chromosome carrying female-inducing modifier gene(s).

Molecular mechanism of male differentiation is conserved in the SRY-absent mammal, Tokudaia osimensis

The sex-determining gene SRY induces SOX9 expression in the testes of eutherian mammals via two pathways. SRY binds to testis-specific enhancer of Sox9 (TESCO) with SF1 to activate SOX9 transcription. SRY also up-regulates ER71 expression, and ER71 activates Sox9 transcription. After the initiation of testis differentiation, SOX9 enhances Amh expression by binding to its promoter with SF1. SOX8, SOX9 and SOX10, members of the SOXE gene family, also enhance the activities of the Amh promoter and TESCO. In this study, we investigated the regulation of these sexual differentiation genes in Tokudaia osimensis, which lacks a Y chromosome and the SRY gene. The activity of the AMH promoter was stimulated by SOXE genes and SF1. Mutant AMH promoters, with mutations in its SOX and SF1 binding sites, did not show significant activity by SOX9 and SF1. These results indicate that AMH expression was regulated by the binding of SOX9 and SF1. By contrast, SOXE genes could not enhance TESCO activity. These results indicate that TESCO enhancer activity was lost in this species. Furthermore, the activity of the SOX9 promoter was enhanced by ER71, indicating that ER71 may play an important role in the testis-specific expression of SOX9.

Unique sex chromosome systems in Ellobius: How do male XX chromosomes recombine and undergo pachytene chromatin inactivation?

Most mammalian species have heteromorphic sex chromosomes in males, except for a few enigmatic groups such as the mole voles Ellobius, which do not have the Y chromosome and Sry gene. The Ellobius (XX ♀♂) system of sex chromosomes has no analogues among other animals. The structure and meiotic behaviour of the two X chromosomes were investigated for males of the sibling species Ellobius talpinus and Ellobius tancrei. Their sex chromosomes, despite their identical G-structure, demonstrate short synaptic fragments and crossover-associated MLH1 foci in both telomeric regions only. The chromatin undergoes modifications in the meiotic sex chromosomes. SUMO-1 marks a small nucleolus-like body of the meiotic XX. ATR and ubiH2A are localized in the asynaptic area and the histone γH2AFX covers the entire XX bivalent. The distribution of some markers of chromatin inactivation differentiates sex chromosomes of mole voles from those of other mammals. Sex chromosomes of both studied species have identical recombination and meiotic inactivation patterns. In Ellobius, similar chromosome morphology masks the functional heteromorphism of the male sex chromosomes, which can be seen at meiosis.

Did sex chromosome turnover promote divergence of the major mammal groups?: De novo sex chromosomes and drastic rearrangements may have posed reproductive barriers between monotremes, marsupials and placental mammals

Comparative mapping and sequencing show that turnover of sex determining genes and chromosomes, and sex chromosome rearrangements, accompany speciation in many vertebrates. Here I review the evidence and propose that the evolution of therian mammals was precipitated by evolution of the male‐determining SRY gene, defining a novel XY sex chromosome pair, and interposing a reproductive barrier with the ancestral population of synapsid reptiles 190 million years ago (MYA). Divergence was reinforced by multiple translocations in monotreme sex chromosomes, the first of which supplied a novel sex determining gene. A sex chromosome‐autosome fusion may have separated eutherians (placental mammals) from marsupials 160 MYA. Another burst of sex chromosome change and speciation is occurring in rodents, precipitated by the degradation of the Y. And although primates have a more stable Y chromosome, it may be just a matter of time before the same fate overtakes our own lineage.

Also watch the video abstract.

Porcine SOX9 Gene Expression Is Influenced by an 18 bp Indel in the 5'-Untranslated Region

Sex determining region Y-box 9 (SOX9) is an important regulator of sex and skeletal development and is expressed in a variety of embryonal and adult tissues. Loss or gain of function resulting from mutations within the coding region or chromosomal aberrations of the SOX9 locus lead to a plethora of detrimental phenotypes in humans and animals. One of these phenotypes is the so-called male-to-female or female-to-male sex-reversal which has been observed in several mammals including pig, dog, cat, goat, horse, and deer. In 38,XX sex-reversal French Large White pigs, a genome-wide association study suggested SOX9 as the causal gene, although no functional mutations were identified in affected animals. However, besides others an 18bp indel had been detected in the 5′-untranslated region of the SOX9 gene by comparing affected animals and controls. We have identified the same indel (Δ18) between position +247bp and +266bp downstream the transcription start site of the porcine SOX9 gene in four other pig breeds; i.e., German Large White, Laiwu Black, Bamei, and Erhualian. These animals have been genotyped in an attempt to identify candidate genes for porcine inguinal and/or scrotal hernia. Because the 18bp segment in the wild type 5′-UTR harbours a highly conserved cAMP-response element (CRE) half-site, we analysed its role in SOX9 expression in vitro. Competition and immunodepletion electromobility shift assays demonstrate that the CRE half-site is specifically recognized by CREB. Both binding of CREB to the wild type as well as the absence of the CRE half-site in Δ18 reduced expression efficiency in HEK293T, PK–15, and ATDC5 cells significantly. Transfection experiments of wild type and Δ18 SOX9 promoter luciferase constructs show a significant reduction of RNA and protein levels depending on the presence or absence of the 18bp segment. Hence, the data presented here demonstrate that the 18bp indel in the porcine SOX9 5′-UTR is of functional importance and may therefore indeed be a causative variation in SOX9 associated traits.