Understanding sex determination in the mouse: genetics, epigenetics and the story of mutual antagonisms

Epigenetic disruptions in the offspring hypothalamus in response to maternal infection

DNA methylation is an epigenetic modification that can alter gene expression, and the incidence can vary across developmental stages, inflammatory conditions, and sexes. The effects of viral maternal viral infection and sex on the DNA methylation patterns were studied in the hypothalamus of a pig model of immune activation during development. DNA methylation at single-base resolution in regions of high CpG density was measured on 24 individual hypothalamus samples using reduced representation bisulfite sequencing. Differential over- and under-methylated sites were identified and annotated to proximal genes and corresponding biological processes. A total of 120 sites were differentially methylated (FDR-adjusted p-value < 0.05) between maternal infection or sex groups. Among the 66 sites differentially methylated between groups exposed to inflammatory signals and control, most sites were over-methylated in the challenged group and included sites in the promoter regions of genes SIRT3 and NRBP1. Among the 54 differentially methylated sites between females and males, most sites were over-methylated in females and included sites in the promoter region of genes TNC and EIF4G1. The analysis of the genes proximal to the differentially methylated sites suggested that biological processes potentially impacted include immune response, neuron migration and ensheathment, peptide signaling, adaptive thermogenesis, and tissue development. These results suggest that translational studies should consider that the prolonged effect of maternal infection during gestation may be enacted through epigenetic regulatory mechanisms that may differ between sexes.

Decoding the epigenetic mechanism of mammalian sex determination

Sex determination embodies a dynamic and intricate developmental process wielding significant influence over the destiny of bipotential gonads, steering them towards male or female gonads. Gonadal differentiation and the postnatal manifestation of the gonadal phenotype involve a sophisticated interplay of transcription factors such as SOX9 and FOXL2. Central to this interplay are chromatin modifiers regulating the mutual antagonism during this interplay. In this review, the key findings and knowledge gaps in DNA methylation, histone modification, and non-coding RNA-mediated control throughout mammalian gonadal development are covered. Furthermore, it explores the role of the developing brain in playing a pivotal role in the initiation of gonadogenesis and the subsequent involvement of gonadal hormone/hormone receptor in fine-tuning sexual differentiation. Based on promising facts, the role of the developing brain through the hypothalamic pituitary gonadal axis is explained and suggested as a novel hypothesis. The article also discusses the potential impact of ecological factors on the human epigenome in relation to sex determination and trans-generational epigenetics in uncovering novel genes and mechanisms involved in sex determination and gonadal differentiation. We have subtly emphasized the disruptions in epigenetic regulations contributing to sexual disorders, which further allows us to raise certain questions, decipher approaches for handling these questions and setting up the direction of future research.

Deciphering Sex-Specific Differentiation of Human Fetal Gonads: Insight From Experimental Models

Sex-specific gonadal differentiation is initiated by the expression of SRY in male foetuses. This promotes a signalling pathway directing testicular development, while in female foetuses the absence of SRY and expression of pro-ovarian factors promote ovarian development. Importantly, in addition to the initiation of a sex-specific signalling cascade the opposite pathway is simultaneously inhibited. The somatic cell populations within the gonads dictates this differentiation as well as the development of secondary sex characteristics via secretion of endocrine factors and steroid hormones. Opposing pathways SOX9/FGF9 (testis) and WNT4/RSPO1 (ovary) controls the development and differentiation of the bipotential mouse gonad and even though sex-specific gonadal differentiation is largely considered to be conserved between mice and humans, recent studies have identified several differences. Hence, the signalling pathways promoting early mouse gonad differentiation cannot be directly transferred to human development thus highlighting the importance of also examining this signalling in human fetal gonads. This review focus on the current understanding of regulatory mechanisms governing human gonadal sex differentiation by combining knowledge of these processes from studies in mice, information from patients with differences of sex development and insight from manipulation of selected signalling pathways in ex vivo culture models of human fetal gonads.

Broad-spectrum XX and XY gonadal dysgenesis in patients with a homozygous L193S variant in PPP2R3C

CONTEXT

Homozygous and heterozygous variants in PPP2R3C are associated with syndromic 46,XY complete gonadal dysgenesis (Myo-Ectodermo-Gonadal Dysgenesis (MEGD) syndrome), and impaired spermatogenesis, respectively. This study expands the role of PPP2R3C in the aetiology of gonadal dysgenesis (GD).

METHOD

We sequenced the PPP2R3C gene in four new patients from three unrelated families. The clinical, laboratory, and molecular characteristics were investigated. We have also determined the requirement for Ppp2r3c in mice (C57BL6/N) using CRISPR/Cas9 genome editing.

RESULTS

A homozygous c.578T>C (p.L193S) PPP2R3C variant was identified in one 46,XX girl with primary gonadal insufficiency, two girls with 46,XY complete GD, and one undervirilised boy with 46,XY partial GD. The patients with complete GD had low gonadal and adrenal androgens, low anti-Müllerian hormone, and high follicle-stimulating hormone and luteinizing hormone concentrations. All patients manifested characteristic features of MEGD syndrome. Heterozygous Ppp2r3c knockout mice appeared overtly normal and fertile. Inspection of homozygous embryos at 14.5, 9.5, and 8.5 days post coitum(dpc) revealed evidence of dead embryos. We conclude that loss of function of Ppp2r3c is not compatible with viability in mice and results in embryonic death from 7.5 dpc or earlier.

CONCLUSION

Our data indicate the essential roles for PPP2R3C in mouse and human development. Germline homozygous variants in human PPP2R3C are associated with distinctive syndromic GD of varying severity in both 46,XY and 46,XX individuals.

Sex-Biased Gene Expression and Isoform Profile of Brine Shrimp Artemia franciscana by Transcriptome Analysis

Simple Summary

The brine shrimp Artemia is a promising model organism for ZW sex determination system, but the genes related to sex determination and differentiation of Artemia have not yet been examined in detail. In this study, the first isoform-level transcriptome sequencing was performed on female and male Artemia franciscana. By using PacBio Iso-Seq and Illumina RNA-Seq technologies, we found 39 candidate sex determination genes that showed sex-biased gene expression. The male-biased expressed genes included DMRT1 and Sad genes, which had three and seven isoforms, respectively. Among these, the Sad gene is an ecdysteroid biosynthetic pathway gene associated with arthropod molting and metamorphosis. We propose the importance and the necessity of further research on genes involved in ecdysteroid biosynthesis. These results will contribute to understand sex determination and differentiation of Artemia and other crustaceans having ZW systems.

Abstract

The brine shrimp Artemia has a ZW sex determination system with ZW chromosomes in females and ZZ chromosomes in males. Artemia has been considered a promising model organism for ZW sex-determining systems, but the genes involved in sex determination and differentiation of Artemia have not yet been identified. Here, we conducted transcriptome sequencing of female and male A. franciscana using PacBio Iso-Seq and Illumina RNA-Seq techniques to identify candidate sex determination genes. Among the 42,566 transcripts obtained from Iso-Seq, 23,514 were analyzed. Of these, 2065 (8.8%) were female specific, 2513 (10.7%) were male specific, and 18,936 (80.5%) were co-expressed in females and males. Based on GO enrichment analysis and expression values, we found 10 female-biased and 29 male-biased expressed genes, including DMRT1 and Sad genes showing male-biased expression. Our results showed that DMRT1 has three isoforms with five exons, while Sad has seven isoforms with 2–11 exons. The Sad gene is involved in ecdysteroid signaling related to molting and metamorphosis in arthropods. Further studies on ecdysteroid biosynthetic genes are needed to improve our understanding of Artemia sex determination. This study will provide a valuable resource for sex determination and differentiation studies on Artemia and other crustaceans with ZW systems.

The effect of sex on the mouse lens transcriptome

The transcriptome of mammalian tissues differs between males and females, and these differences can change across the lifespan, likely regulating known sexual dimorphisms in disease prevalence and severity. Cataract, the most prevalent disease of the ocular lens, occurs at similar rates in young individuals, but its incidence is elevated in older women compared to men of the same age. However, the influence of sex on the lens transcriptome was unknown. RNAseq based transcriptomic profiling of young adult C57BL/6J mouse lens epithelial and fiber cells revealed that few genes are differentially expressed between the sexes. In contrast, lens cells from aged (24 month old) male and female C57BL/6J mice differentially expressed many genes, including several whose expression is lens preferred. Like cataracts, posterior capsular opacification (PCO), a major sequela of cataract surgery, may also be more prevalent in women. Lens epithelial cells isolated from mouse eyes 24 h after lens fiber cell removal exhibited numerous transcriptomic differences between the sexes, including genes implicated in complement cascades and extracellular matrix regulation, and these differences are much more pronounced in aged mice than in young mice. These results provide an unbiased basis for future studies on how sex affects the lens response to aging, cataract development, and cataract surgery.

Dynamics of the transcriptional landscape during human fetal testis and ovary development

STUDY QUESTION

Which transcriptional program triggers sex differentiation in bipotential gonads and downstream cellular events governing fetal testis and ovary development in humans?

SUMMARY ANSWER

The characterization of a dynamically regulated protein-coding and non-coding transcriptional landscape in developing human gonads of both sexes highlights a large number of potential key regulators that show an early sexually dimorphic expression pattern.

WHAT IS KNOWN ALREADY

Gonadal sex differentiation is orchestrated by a sexually dimorphic gene expression program in XX and XY developing fetal gonads. A comprehensive characterization of its non-coding counterpart offers promising perspectives for deciphering the molecular events underpinning gonad development and for a complete understanding of the etiology of disorders of sex development in humans.

STUDY DESIGN, SIZE, DURATION

To further investigate the protein-coding and non-coding transcriptional landscape during gonad differentiation, we used RNA-sequencing (RNA-seq) and characterized the RNA content of human fetal testis (N = 24) and ovaries (N = 24) from 6 to 17 postconceptional week (PCW), a key period in sex determination and gonad development.

PARTICIPANTS/MATERIALS, SETTING, METHODS

First trimester fetuses (6-12 PCW) and second trimester fetuses (13-14 and 17 PCW) were obtained from legally induced normally progressing terminations of pregnancy. Total RNA was extracted from whole human fetal gonads and sequenced as paired-end 2 × 50 base reads. Resulting sequences were mapped to the human genome, allowing for the assembly and quantification of corresponding transcripts.

MAIN RESULTS AND THE ROLE OF CHANCE

This RNA-seq analysis of human fetal testes and ovaries at seven key developmental stages led to the reconstruction of 22 080 transcripts differentially expressed during testicular and/or ovarian development. In addition to 8935 transcripts displaying sex-independent differential expression during gonad development, the comparison of testes and ovaries enabled the discrimination of 13 145 transcripts that show a sexually dimorphic expression profile. The latter include 1479 transcripts differentially expressed as early as 6 PCW, including 39 transcription factors, 40 long non-coding RNAs and 20 novel genes. Despite the use of stringent filtration criteria (expression cut-off of at least 1 fragment per kilobase of exon model per million reads mapped, fold change of at least 2 and false discovery rate adjusted P values of less than <1%), the possibility of assembly artifacts and of false-positive differentially expressed transcripts cannot be fully ruled out.

LARGE-SCALE DATA

Raw data files (fastq) and a searchable table (.xlss) containing information on genomic features and expression data for all refined transcripts have been submitted to the NCBI GEO under accession number GSE116278.

LIMITATIONS, REASONS FOR CAUTION

The intrinsic nature of this bulk analysis, i.e. the sequencing of transcripts from whole gonads, does not allow direct identification of the cellular origin(s) of the transcripts characterized. Potential cellular dilution effects (e.g. as a result of distinct proliferation rates in XX and XY gonads) may account for a few of the expression profiles identified as being sexually dimorphic. Finally, transcriptome alterations that would result from exposure to pre-abortive drugs cannot be completely excluded. Although we demonstrated the high quality of the sorted cell populations used for experimental validations using quantitative RT-PCR, it cannot be totally excluded that some germline expression may correspond to cell contamination by, for example, macrophages.

WIDER IMPLICATIONS OF THE FINDINGS

For the first time, this study has led to the identification of 1000 protein-coding and non-coding candidate genes showing an early, sexually dimorphic, expression pattern that have not previously been associated with sex differentiation. Collectively, these results increase our understanding of gonad development in humans, and contribute significantly to the identification of new candidate genes involved in fetal gonad differentiation. The results also provide a unique resource that may improve our understanding of the fetal origin of testicular and ovarian dysgenesis syndromes, including cryptorchidism and testicular cancers.

STUDY FUNDING/COMPETING INTEREST(S)

This work was supported by the French National Institute of Health and Medical Research (Inserm), the University of Rennes 1, the French School of Public Health (EHESP), the Swiss National Science Foundation [SNF n° CRS115_171007 to B.J.], the French National Research Agency [ANR n° 16-CE14-0017-02 and n° 18-CE14-0038-02 to F.C.], the Medical Research Council [MR/L010011/1 to P.A.F.] and the European Community's Seventh Framework Programme (FP7/2007-2013) [under grant agreement no 212885 to P.A.F.] and from the European Union's Horizon 2020 Research and Innovation Programme [under grant agreement no 825100 to P.A.F. and S.M.G.]. There are no competing interests related to this study.

Sex determination, gonadal sex differentiation, and plasticity in vertebrate species

A diverse array of sex determination (SD) mechanisms, encompassing environmental to genetic, have been found to exist among vertebrates, covering a spectrum from fixed SD mechanisms (mammals) to functional sex change in fishes (sequential hermaphroditic fishes). A major landmark in vertebrate SD was the discovery of the SRY gene in 1990. Since that time, many attempts to clone an SRY ortholog from nonmammalian vertebrates remained unsuccessful, until 2002, when DMY/dmrt1by was discovered as the SD gene of a small fish, medaka. Surprisingly, however, DMY/dmrt1by was found in only 2 species among more than 20 species of medaka, suggesting a large diversity of SD genes among vertebrates. Considerable progress has been made over the last 3 decades, such that it is now possible to formulate reasonable paradigms of how SD and gonadal sex differentiation may work in some model vertebrate species. This review outlines our current understanding of vertebrate SD and gonadal sex differentiation, with a focus on the molecular and cellular mechanisms involved. An impressive number of genes and factors have been discovered that play important roles in testicular and ovarian differentiation. An antagonism between the male and female pathway genes exists in gonads during both sex differentiation and, surprisingly, even as adults, suggesting that, in addition to sex-changing fishes, gonochoristic vertebrates including mice maintain some degree of gonadal sexual plasticity into adulthood. Importantly, a review of various SD mechanisms among vertebrates suggests that this is the ideal biological event that can make us understand the evolutionary conundrums underlying speciation and species diversity.

Transcriptomic analysis of mRNA expression and alternative splicing during mouse sex determination

Mammalian sex determination hinges on sexually dimorphic transcriptional programs in developing fetal gonads. A comprehensive view of these programs is crucial for understanding the normal development of fetal testes and ovaries and the etiology of human disorders of sex development (DSDs), many of which remain unexplained. Using strand-specific RNA-sequencing, we characterized the mouse fetal gonadal transcriptome from 10.5 to 13.5 days post coitum, a key time window in sex determination and gonad development. Our dataset benefits from a greater sensitivity, accuracy and dynamic range compared to microarray studies, allows global dynamics and sex-specificity of gene expression to be assessed, and provides a window to non-transcriptional events such as alternative splicing. Spliceomic analysis uncovered female-specific regulation of Lef1 splicing, which may contribute to the enhanced WNT signaling activity in XX gonads. We provide a user-friendly visualization tool for the complete transcriptomic and spliceomic dataset as a resource for the field.

The molecular pathways underlying early gonadal development

The body of knowledge surrounding reproductive development spans the fields of genetics, anatomy, physiology and biomedicine, to build a comprehensive understanding of the later stages of reproductive development in humans and animal models. Despite this, there remains much to learn about the bi-potential progenitor structure that the ovary and testis arise from, known as the genital ridge (GR). This tissue forms relatively late in embryonic development and has the potential to form either the ovary or testis, which in turn produce hormones required for development of the rest of the reproductive tract. It is imperative that we understand the genetic networks underpinning GR development if we are to begin to understand abnormalities in the adult. This is particularly relevant in the contexts of disorders of sex development (DSDs) and infertility, two conditions that many individuals struggle with worldwide, with often no answers as to their aetiology. Here, we review what is known about the genetics of GR development. Investigating the genetic networks required for GR formation will not only contribute to our understanding of the genetic regulation of reproductive development, it may in turn open new avenues of investigation into reproductive abnormalities and later fertility issues in the adult.

Single cell transcriptome sequencing: A new approach for the study of mammalian sex determination

Mammalian sex determination is a highly complex developmental process that is particularly difficult to study due to the limited number of gonadal cells present at the bipotential stage, the large cellular heterogeneity in both testis and ovaries and the rapid sex-dependent differentiation processes. Single-cell RNA-sequencing (scRNA-seq) circumvents the averaging artifacts associated with methods traditionally used to profile bulk populations of cells. It is a powerful tool that allows the identification and classification of cell populations in a comprehensive and unbiased manner. In particular, scRNA-seq enables the tracing of cells along developmental trajectories and characterization of the transcriptional dynamics controlling their differentiation. In this review, we describe the current state-of-the-art experimental methods used for scRNA-seq and discuss their strengths and limitations. Additionally, we summarize the multiple key insights that scRNA-seq has provided to the understanding of mammalian sex determination. Finally, we briefly discuss the future of this technology, as well as complementary applications in single cell -omics in the context of mammalian sex determination.

SOX4 regulates gonad morphogenesis and promotes male germ cell differentiation in mice

The group C SOX transcription factors SOX4, -11 and -12 play important and mutually overlapping roles in development of a number of organs. Here, we examined the role of SoxC genes during gonadal development in mice. All three genes were expressed in developing gonads of both sexes, predominantly in somatic cells, with Sox4 being most strongly expressed. Sox4 deficiency resulted in elongation of both ovaries and testes, and an increased number of testis cords. While female germ cells entered meiosis normally, male germ cells showed reduced levels of differentiation markers Nanos2 and Dnmt3l and increased levels of pluripotency genes Cripto and Nanog, suggesting that SOX4 may normally act to restrict the pluripotency period of male germ cells and ensure their proper differentiation. Finally, our data reveal that SOX4 (and, to a lesser extent, SOX11 and -12) repressed transcription of the sex-determining gene Sox9 via an upstream testis-specific enhancer core (TESCO) element in fetal gonads, raising the possibility that SOXC proteins may function as transcriptional repressors in a context-dependent manner.