Expression profiling of purified mouse gonadal somatic cells during the critical time window of sex determination reveals novel candidate genes for human sexual dysgenesis syndromes

Genetic reanalysis of patients with a difference of sex development carrying the NR5A1/SF-1 variant p.Gly146Ala has discovered other likely disease-causing variations

NR5A1/SF-1 (Steroidogenic factor-1) variants may cause mild to severe differences of sex development (DSD) or may be found in healthy carriers. The NR5A1/SF-1 c.437G>C/p.Gly146Ala variant is common in individuals with a DSD and has been suggested to act as a susceptibility factor for adrenal disease or cryptorchidism. Since the allele frequency is high in the general population, and the functional testing of the p.Gly146Ala variant revealed inconclusive results, the disease-causing effect of this variant has been questioned. However, a role as a disease modifier is still possible given that oligogenic inheritance has been described in patients with NR5A1/SF-1 variants. Therefore, we performed next generation sequencing (NGS) in 13 DSD individuals harboring the NR5A1/SF-1 p.Gly146Ala variant to search for other DSD-causing variants and clarify the function of this variant for the phenotype of the carriers. Panel and whole-exome sequencing was performed, and data were analyzed with a filtering algorithm for detecting variants in NR5A1- and DSD-related genes. The phenotype of the studied individuals ranged from scrotal hypospadias and ambiguous genitalia in 46,XY DSD to opposite sex in both 46,XY and 46,XX. In nine subjects we identified either a clearly pathogenic DSD gene variant (e.g. in AR) or one to four potentially deleterious variants that likely explain the observed phenotype alone (e.g. in FGFR3, CHD7). Our study shows that most individuals carrying the NR5A1/SF-1 p.Gly146Ala variant, harbor at least one other deleterious gene variant which can explain the DSD phenotype. This finding confirms that the NR5A1/SF-1 p.Gly146Ala variant may not contribute to the pathogenesis of DSD and qualifies as a benign polymorphism. Thus, individuals, in whom the NR5A1/SF-1 p.Gly146Ala gene variant has been identified as the underlying genetic cause for their DSD in the past, should be re-evaluated with a NGS method to reveal the real genetic diagnosis.

Functional Analysis of Mmd2 and Related PAQR Genes During Sex Determination in Mice

INTRODUCTION

Sex determination in eutherian mammals is controlled by the Y-linked gene Sry, which drives the formation of testes in male embryos. Despite extensive study, the genetic steps linking Sry action and male sex determination remain largely unknown. Here, we focused on Mmd2, a gene that encodes a member of the progestin and adipoQ receptor (PAQR) family. Mmd2 is expressed during the sex-determining period in XY but not XX gonads, suggesting a specific role in testis development.

METHODS

We used CRISPR to generate mouse strains deficient in Mmd2 and its 2 closely related PAQR family members, Mmd and Paqr8, which are also expressed during testis development. Following characterization of Mmd2 expression in the developing testis, we studied sex determination in embryos from single knockout as well as Mmd2;Mmd and Mmd2;Paqr8 double knockout lines using quantitative RT-PCR and immunofluorescence.

RESULTS

Analysis of knockout mice deficient in Sox9 and Nr5a1 revealed that Mmd2 operates downstream of these known sex-determining genes. However, fetal testis development progressed normally in Mmd2-null embryos. To determine if other genes might have compensated for the loss of Mmd2, we analyzed Paqr8 and Mmd-null embryos and confirmed that in both knockout lines, sex determination occurred normally. Finally, we generated Mmd2;Mmd and Mmd2;Paqr8 double-null embryos and again observed normal testis development.

DISCUSSION

These results may reflect functional redundancy among PAQR factors, or their dispensability in gonadal development. Our findings highlight the difficulties involved in identifying genes with a functional role in sex determination and gonadal development through expression screening and loss-of-function analyses of individual candidate genes and may help to explain the paucity of genes in which variations have been found to cause human disorders/differences of sex development.

Gonadal Sex Differentiation and Ovarian Organogenesis along the Cortical-Medullary Axis in Mammals

In most mammals, the sex of the gonads is based on the fate of the supporting cell lineages, which arises from the proliferation of coelomic epithelium (CE) that surfaces on the bipotential genital ridge in both XY and XX embryos. Recent genetic studies and single-cell transcriptome analyses in mice have revealed the cellular and molecular events in the two-wave proliferation of the CE that produce the supporting cells. This proliferation contributes to the formation of the primary sex cords in the medullary region of both the testis and the ovary at the early phase of gonadal sex differentiation, as well as to that of the secondary sex cords in the cortical region of the ovary at the perinatal stage. To support gametogenesis, the testis forms seminiferous tubules in the medullary region, whereas the ovary forms follicles mainly in the cortical region. The medullary region in the ovary exhibits morphological and functional diversity among mammalian species that ranges from ovary-like to testis-like characteristics. This review focuses on the mechanism of gonadal sex differentiation along the cortical-medullary axis and compares the features of the cortical and medullary regions of the ovary in mammalian species.

Candidate genes for infertility: an in-silico study based on cytogenetic analysis

Background

The cause of infertility remains unclear in a significant proportion of reproductive-age couples who fail to conceive naturally. Chromosomal aberrations have been identified as one of the main genetic causes of male and female infertility. Structural chromosomal aberrations may disrupt the functioning of various genes, some of which may be important for fertility. The present study aims to identify candidate genes and putative functional interaction networks involved in male and female infertility using cytogenetic data from cultured peripheral blood lymphocytes of infertile patients.

Methods

Karyotypic analyses was done in 201 infertile patients (100 males and 101 females) and 201 age and gender matched healthy controls (100 males and 101 females) after 72 h peripheral lymphocyte culturing and GTG banding, followed by bioinformatic analysis using Cytoscape v3.8.2 and Metascape.

Results

Several chromosomal regions with a significantly higher frequency of structural aberrations were identified in the infertile males (5q2, 10q2, and 17q2) and females (6q2, 16q2, and Xq2). Segregation of the patients based on type of infertility (primary v/s secondary infertility) led to the identification of chromosomal regions with a significantly higher frequency of structural aberrations exclusively within the infertile males (5q2, 17q2) and females (16q2) with primary infertility. Cytoscape identified two networks specific to these regions: a male specific network with 99 genes and a female specific network with 109 genes. The top enriched GO terms within the male and female infertility networks were “skeletal system morphogenesis” and “mRNA transport” respectively. PSME3, PSMD3, and CDC27 were the top 3 hub genes identified within the male infertility network. Similarly, UPF3B, IRF8, and PSMB1 were the top 3 hub genes identified with the female infertility network. Among the hub genes identified in the male- and female-specific networks, PSMB1, PSMD3, and PSME3 are functional components of the proteasome complex. These hub genes have a limited number of reports related to their respective roles in maintenance of fertility in mice model and humans and require validation in further studies.

Conclusion

The candidate genes predicted in the present study can serve as targets for future research on infertility.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01320-x.

Somatic regulation of female germ cell regeneration and development in planarians

Female germ cells develop into oocytes, with the capacity for totipotency. In most animals, these remarkable cells are specified during development and cannot be regenerated. By contrast, planarians, known for their regenerative prowess, can regenerate germ cells. To uncover mechanisms required for female germ cell development and regeneration, we generated gonad-specific transcriptomes and identified genes whose expression defines progressive stages of female germ cell development. Strikingly, early female germ cells share molecular signatures with the pluripotent stem cells driving planarian regeneration. We observe spatial heterogeneity within somatic ovarian cells and find that a regionally enriched foxL homolog is required for oocyte differentiation, but not specification, suggestive of functionally distinct somatic compartments. Unexpectedly, a neurotransmitter-biosynthetic enzyme, aromatic L-amino acid decarboxylase (AADC), is also expressed in somatic gonadal cells, and plays opposing roles in female and male germ cell development. Thus, somatic gonadal cells deploy conserved factors to regulate germ cell development and regeneration in planarians.

Unlike most animals, planarians can regenerate germ cells. Here, Khan and Newmark characterize gene expression in the planarian ovary, identifying genes expressed at progressive stages of female germ cell development and in somatic ovarian cells. Functional characterization revealed somatically expressed genes required for specification or differentiation of female germ cells.

Sertoli Cell-Germ Cell Interactions Within the Niche: Paracrine and Juxtacrine Molecular Communications

Male germ cell development depends on multiple biological events that combine epigenetic reprogramming, cell cycle regulation, and cell migration in a spatio-temporal manner. Sertoli cells are a crucial component of the spermatogonial stem cell niche and provide essential growth factors and chemokines to developing germ cells. This review focuses mainly on the activation of master regulators of the niche in Sertoli cells and their targets, as well as on novel molecular mechanisms underlying the regulation of growth and differentiation factors such as GDNF and retinoic acid by NOTCH signaling and other pathways.

Sex Maintenance in Mammals

The crucial event in mammalian sexual differentiation occurs at the embryonic stage of sex determination, when the bipotential gonads differentiate as either testes or ovaries, according to the sex chromosome constitution of the embryo, XY or XX, respectively. Once differentiated, testes produce sexual hormones that induce the subsequent differentiation of the male reproductive tract. On the other hand, the lack of masculinizing hormones in XX embryos permits the formation of the female reproductive tract. It was long assumed that once the gonad is differentiated, this developmental decision is irreversible. However, several findings in the last decade have shown that this is not the case and that a continuous sex maintenance is needed. Deletion of Foxl2 in the adult ovary lead to ovary-to-testis transdifferentiation and deletion of either Dmrt1 or Sox9/Sox8 in the adult testis induces the opposite process. In both cases, mutant gonads were genetically reprogrammed, showing that both the male program in ovaries and the female program in testes must be actively repressed throughout the individual's life. In addition to these transcription factors, other genes and molecular pathways have also been shown to be involved in this antagonism. The aim of this review is to provide an overview of the genetic basis of sex maintenance once the gonad is already differentiated.

Specific Transcriptomic Signatures and Dual Regulation of Steroidogenesis Between Fetal and Adult Mouse Leydig Cells

Leydig cells (LC) are the main testicular androgen-producing cells. In eutherian mammals, two types of LCs emerge successively during testicular development, fetal Leydig cells (FLCs) and adult Leydig cells (ALCs). Both display significant differences in androgen production and regulation. Using bulk RNA sequencing, we compared the transcriptomes of both LC populations to characterize their specific transcriptional and functional features. Despite similar transcriptomic profiles, a quarter of the genes show significant variations in expression between FLCs and ALCs. Non-transcriptional events, such as alternative splicing was also observed, including a high rate of intron retention in FLCs compared to ALCs. The use of single-cell RNA sequencing data also allowed the identification of nine FLC-specific genes and 50 ALC-specific genes. Expression of the corticotropin-releasing hormone 1 (Crhr1) receptor and the ACTH receptor melanocortin type 2 receptor (Mc2r) specifically in FLCs suggests a dual regulation of steroidogenesis. The androstenedione synthesis by FLCs is stimulated by luteinizing hormone (LH), corticotrophin-releasing hormone (CRH), and adrenocorticotropic hormone (ACTH) whereas the testosterone synthesis by ALCs is dependent exclusively on LH. Overall, our study provides a useful database to explore LC development and functions.

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.

Applying Single-Cell Analysis to Gonadogenesis and DSDs (Disorders/Differences of Sex Development)

The gonads are unique among the body's organs in having a developmental choice: testis or ovary formation. Gonadal sex differentiation involves common progenitor cells that form either Sertoli and Leydig cells in the testis or granulosa and thecal cells in the ovary. Single-cell analysis is now shedding new light on how these cell lineages are specified and how they interact with the germline. Such studies are also providing new information on gonadal maturation, ageing and the somatic-germ cell niche. Furthermore, they have the potential to improve our understanding and diagnosis of Disorders/Differences of Sex Development (DSDs). DSDs occur when chromosomal, gonadal or anatomical sex are atypical. Despite major advances in recent years, most cases of DSD still cannot be explained at the molecular level. This presents a major pediatric concern. The emergence of single-cell genomics and transcriptomics now presents a novel avenue for DSD analysis, for both diagnosis and for understanding the molecular genetic etiology. Such -omics datasets have the potential to enhance our understanding of the cellular origins and pathogenesis of DSDs, as well as infertility and gonadal diseases such as cancer.

Regulation of GDNF expression in Sertoli cells

Sertoli cells regulate male germ cell proliferation and differentiation and are a critical component of the spermatogonial stem cell (SSC) niche, where homeostasis is maintained by the interplay of several signaling pathways and growth factors. These factors are secreted by Sertoli cells located within the seminiferous epithelium, and by interstitial cells residing between the seminiferous tubules. Sertoli cells and peritubular myoid cells produce glial cell line-derived neurotrophic factor (GDNF), which binds to the RET/GFRA1 receptor complex at the surface of undifferentiated spermatogonia. GDNF is known for its ability to drive SSC self-renewal and proliferation of their direct cell progeny. Even though the effects of GDNF are well studied, our understanding of the regulation its expression is still limited. The purpose of this review is to discuss how GDNF expression in Sertoli cells is modulated within the niche, and how these mechanisms impact germ cell homeostasis.

Dissecting Cell Lineage Specification and Sex Fate Determination in Gonadal Somatic Cells Using Single-Cell Transcriptomics

Sex determination is a unique process that allows the study of multipotent progenitors and their acquisition of sex-specific fates during differentiation of the gonad into a testis or an ovary. Using time series single-cell RNA sequencing (scRNA-seq) on ovarian Nr5a1-GFP⁺ somatic cells during sex determination, we identified a single population of early progenitors giving rise to both pre-granulosa cells and potential steroidogenic precursor cells. By comparing time series single-cell RNA sequencing of XX and XY somatic cells, we provide evidence that gonadal supporting cells are specified from these early progenitors by a non-sex-specific transcriptomic program before pre-granulosa and Sertoli cells acquire their sex-specific identity. In XX and XY steroidogenic precursors, similar transcriptomic profiles underlie the acquisition of cell fate but with XX cells exhibiting a relative delay. Our data provide an important resource, at single-cell resolution, for further interrogation of the molecular and cellular basis of mammalian sex determination.

Broad Phenotypes of Disorders/Differences of Sex Development in MAMLD1 Patients Through Oligogenic Disease

Disorders/differences of sex development (DSD) are the result of a discordance between chromosomal, gonadal, and genital sex. DSD may be due to mutations in any of the genes involved in sex determination and development in general, as well as gonadal and/or genital development specifically. MAMLD1 is one of the recognized DSD genes. However, its role is controversial as some MAMLD1 variants are present in normal individuals, several MAMLD1 mutations have wild-type activity in functional studies, and the Mamld1-knockout male mouse presents with normal genitalia and reproduction. We previously tested nine MAMLD1 variants detected in nine 46,XY DSD patients with broad phenotypes for their functional activity, but none of the mutants, except truncated L210X, had diminished transcriptional activity on known target promoters CYP17A1 and HES3. In addition, protein expression of MAMLD1 variants was similar to wild-type, except for the truncated L210X. We hypothesized that MAMLD1 variants may not be sufficient to explain the phenotype in 46,XY DSD individuals, and that further genetic studies should be performed to search for additional hits explaining the broad phenotypes. We therefore performed whole exome sequencing (WES) in seven of these 46,XY patients with DSD and in one 46,XX patient with ovarian insufficiency, who all carried MAMLD1 variants. WES data were filtered by an algorithm including disease-tailored lists of MAMLD1-related and DSD-related genes. Fifty-five potentially deleterious variants in 41 genes were identified; 16/55 variants were reported in genes in association with hypospadias, 8/55 with cryptorchidism, 5/55 with micropenis, and 13/55 were described in relation with female sex development. Patients carried 1-16 variants in 1-16 genes together with their MAMLD1 variation. Network analysis of the identified genes revealed that 23 genes presented gene/protein interactions with MAMLD1. Thus, our study shows that the broad phenotypes of individual DSD might involve multiple genetic variations contributing towards the complex network of sexual development.

CBX2 is required to stabilize the testis pathway by repressing Wnt signaling

XX and XY fetal gonads are initially bipotential, poised between the ovary and testis fate. Multiple lines of evidence suggest that commitment to testis fate requires the repression of genes associated with ovary fate. It was previously shown that loss of CBX2, the subunit of the Polycomb Repressive Complex 1 (PRC1) that binds H3K27me3 and mediates silencing, leads to ovary development in XY mice and humans. While it had been proposed that CBX2 is an activator of the testis-determining gene Sry, we investigated the alternative possibility that CBX2 has a direct role as a repressor of the antagonistic ovary-promoting pathway. To investigate this possibility, we developed a quantitative genome-wide profile of the repressive histone mark H3K27me3 and its active counterpart H3K4me3 in isolated XY and XX gonadal supporting cells before and after sex determination. We show that testis and ovary sex-determining (SD) genes are bivalent before sex determination, providing insight into how the bipotential state of the gonad is established at the epigenetic level. After sex determination, many SD genes of the alternate pathway remain bivalent, possibly contributing to the ability of these cells to transdifferentiate even in adults. The finding that many genes in the Wnt signaling pathway were targeted for H3K27me3-mediated repression in Sertoli cells led us to test whether deletion of Wnt4 could rescue testis development in Cbx2 mutants. We show that Sry expression and testis development were rescued in XY Cbx2^-/-;Wnt4^-/- mice. Furthermore, we show that CBX2 directly binds the downstream Wnt signaler Lef1, an ovary-promoting gene that remains bivalent in Sertoli cells. Our results suggest that stabilization of the testis fate requires CBX2-mediated repression of bivalent ovary-determining genes, which would otherwise block testis development.

During development, the bipotential fetal gonad can commit to the testis fate or to the ovary fate. Mutation of the epigenetic regulator CBX2 leads to ovary development in XY embryos, suggesting a critical role for chromatin remodeling during sex determination. However, the epigenetic mechanisms that regulate the testis vs. ovary cell-fate decision in the mammalian bipotential gonad are poorly understood. In this study, we developed a genome-wide profile of two histone modifications that play critical roles during development: H3K27me3 (repressive) and H3K4me3 (active). We find that sex-determining genes that are initially co-expressed in XX and XY bipotential gonads are bivalent (marked by both H3K4me3 and H3K27me3) prior to sex determination, poised to engage either the testis or ovary fate. Remarkably, after sex determination, repressed genes that promote the alternate fate remain bivalent. We show that stabilization of the testis fate requires CBX2-mediated repression of bivalent ovary-determining genes, which would otherwise block testis development. Our study provides insight into how the bipotential state of the gonad is established at the epigenetic level, and how the testis fate is stabilized by repression of the ovary fate during sex determination.

Transcriptome profiling reveals male- and female-specific gene expression pattern and novel gene candidates for the control of sex determination and gonad development in Xenopus laevis

Xenopus laevis is an amphibian (frog) species widely used in developmental biology and genetics. To unravel the molecular machinery regulating sex differentiation of Xenopus gonads, we analyzed for the first time the transcriptome of developing amphibian gonads covering sex determination period. We applied microarray at four developmental stages: (i) NF50 (undifferentiated gonad during sex determination), (ii) NF53 (the onset of sexual differentiation of the gonads), (iii) NF56 (sexual differentiation of the gonads), and (iv) NF62 (developmental progression of differentiated gonads). Our analysis showed that during the NF50, the genetic female (ZW) gonads expressed more sex-specific genes than genetic male (ZZ) gonads, which suggests that a robust genetic program is realized during female sex determination in Xenopus. However, a contrasting expression pattern was observed at later stages (NF56 and NF62), when the ZW gonads expressed less sex-specific genes than ZZ gonads, i.e., more genes may be involved in further development of the male gonads (ZZ). We identified sexual dimorphism in the expression of several functional groups of genes, including signaling factors, proteases, protease inhibitors, transcription factors, extracellular matrix components, extracellular matrix enzymes, cell adhesion molecules, and epithelium-specific intermediate filaments. In addition, our analysis detected a sexually dimorphic expression of many uncharacterized genes of unknown function, which should be studied further to reveal their identity and if/how they regulate gonad development, sex determination, and sexual differentiation. Comparison between genes sex-specifically expressed in developing gonads of Xenopus and available transcriptome data from zebrafish, two reptile species, chicken, and mouse revealed significant differences in the genetic control of sex determination and gonad development. This shows that the genetic control of gonad development is evolutionarily malleable.

Integrated mRNA and miRNA expression profile analyses reveal the potential roles of sex-biased miRNA-mRNA pairs in gonad tissues of the Chinese concave-eared torrent frog (Odorrana tormota)

The Chinese concave-eared torrent frog (Odorrana tormota) is typically sexually dimorphic. Females are significantly less common than males in the wild. Until now, the molecular mechanisms of reproduction and sex differentiation of frogs remain unclear. Here, we integrated mRNA and microRNA (miRNA) expression profiles to reveal the molecular mechanisms of reproduction and sex differentiation in O. tormota. We identified 234 differentially expressed miRNAs (DEMs) and 18,551 differentially expressed transcripts. Of these, 12,053 mRNAs and 64 miRNAs were upregulated in testes, and 6,498 mRNAs and 170 miRNAs were upregulated in ovaries. Integrated analysis of the miRNA and mRNA expression profiles predicted 75,602 potential miRNA-mRNA interaction sites, with 42,065 negative miRNA-mRNA interactions. We found 36 differentially expressed genes (DEGs) related to reproduction and sex differentiation, of which 15 DEGs formed 92 negative miRNA-mRNA interactions with 34 known DEMs. Thus, miRNAs may play other important roles in O. tormota. Furthermore, Gene Ontology enrichment and Kyoto Encyclopedia of Genes and Genomes pathway analyses showed reproductive-related processes, such as the gonadotropinreleasing hormone signaling pathway and ovarian steroidogenesis. Based on functional annotation and the literature, the retinoic acid signaling pathway, the SOX9-AMH pathway, and the process of spermatogenesis may be involved in the molecular mechanisms of reproduction and sex differentiation in O. tormota, and may be regulated by miRNAs. The miRNA-mRNA pairs described may provide further understanding of the regulatory mechanisms associated with reproduction and sex differentiation, and the molecular mechanism of reproduction in O. tormota.

Genetic Control of Gonadal Sex Determination and Development

Sex determination is the process by which the bipotential gonads develop as either testes or ovaries. With two distinct potential outcomes, the gonadal primordium offers a unique model for the study of cell fate specification and how distinct cell populations diverge from multipotent progenitors. This review focuses on recent advances in our understanding of the genetic programs and epigenetic mechanisms that regulate gonadal sex determination and the regulation of cell fate commitment in the bipotential gonads. We rely primarily on mouse data to illuminate the complex and dynamic genetic programs controlling cell fate decision and sex-specific cell differentiation during gonadal formation and gonadal sex determination.

Gonadal supporting cells acquire sex-specific chromatin landscapes during mammalian sex determination

Cis-regulatory elements are critical for the precise spatiotemporal regulation of genes during development. However, identifying functional regulatory sites that drive cell differentiation in vivo has been complicated by the high numbers of cells required for whole-genome epigenetic assays. Here, we identified putative regulatory elements during sex determination by performing ATAC-seq and ChIP-seq for H3K27ac in purified XX and XY gonadal supporting cells before and after sex determination in mice. We show that XX and XY supporting cells initiate sex determination with similar chromatin landscapes and acquire sex-specific regulatory elements as they commit to the male or female fate. To validate our approach, we identified a functional gonad-specific enhancer downstream of Bmp2, an ovary-promoting gene. This work increases our understanding of the complex regulatory network underlying mammalian sex determination and provides a powerful resource for identifying non-coding regulatory elements that could harbor mutations that lead to Disorders of Sexual Development.

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.

Sex reversal following deletion of a single distal enhancer of Sox9

Cell fate decisions require appropriate regulation of key genes. Sox9, a direct target of SRY, is pivotal in mammalian sex determination. In vivo high-throughput chromatin accessibility techniques, transgenic assays, and genome editing revealed several novel gonadal regulatory elements in the 2-megabase gene desert upstream of Sox9 Although others are redundant, enhancer 13 (Enh13), a 557-base pair element located 565 kilobases 5' from the transcriptional start site, is essential to initiate mouse testis development; its deletion results in XY females with Sox9 transcript levels equivalent to those in XX gonads. Our data are consistent with the time-sensitive activity of SRY and indicate a strict order of enhancer usage. Enh13 is conserved and embedded within a 32.5-kilobase region whose deletion in humans is associated with XY sex reversal, suggesting that it is also critical in humans.

Loss of Glis3 causes dysregulation of retrotransposon silencing and germ cell demise in fetal mouse testis

Fetal germ cell development is regulated by an elaborate combination of cell-extrinsic and cell-intrinsic signals. Here we identify a novel role for the Krüppel-like transcription factor Gli-Similar 3 (Glis3) in male germ cell development in the mouse embryos. Glis3 is expressed in male germ cells during the brief window of time prior to initiation of piRNA-dependent retrotransposon surveillance. Disruption of Glis3 function led to a widespread reduction in retrotransposon silencing factors, aberrant retrotransposon expression and pronounced germ cell loss. Experimental induction of precocious Glis3 expression in vivo before its normal expression resulted in premature expression of several piRNA pathway members, suggesting that GLIS3 is necessary for the activation of the retrotransposon silencing programs. Our findings reveal an unexpected role for GLIS3 in the development of male germ cells and point to a central role for GLIS3 in the control of retrotransposon silencing in the fetal germline.

New transcriptomic tools to understand testis development and functions

The testis plays a central role in the male reproductive system - secreting several hormones including male steroids and producing male gametes. A complex and coordinated molecular program is required for the proper differentiation of testicular cell types and maintenance of their functions in adulthood. The testicular transcriptome displays the highest levels of complexity and specificity across all tissues in a wide range of species. Many studies have used high-throughput sequencing technologies to define the molecular dynamics and regulatory networks in the testis as well as to identify novel genes or gene isoforms expressed in this organ. This review intends to highlight the complementarity of these transcriptomic studies and to show how the use of different sequencing protocols contribute to improve our global understanding of testicular biology.

ZNRF3 functions in mammalian sex determination by inhibiting canonical WNT signaling

Sex determination involves antagonistic interactions between the testis-determining (SRY-SOX9-FGF9) and ovary-promoting (RSPO1-WNT/β-catenin-FOXL2) pathways, but the underlying molecular mechanisms remain unclear. We show that ZNRF3, an E3 ubiquitin ligase that inhibits WNT signaling and is a direct target of RSPO1-mediated membrane clearance, is testis-determining in mice. Testis determination defects in the absence of ZNRF3 arise due to ectopic canonical WNT signaling in XY gonads at the sex-determining stage. We identify human ZNRF3 sequence variants in cases of 46,XY disorders of sex development with XY female presentation. In vitro functional assays show that these variants disrupt ZNRF3 function. Our data reveal a sex-determining role for ZNRF3 and indicate that interactions between ZNRF3 and RSPO1 regulate mammalian sex determination.

Mammalian sex determination is controlled by the antagonistic interactions of two genetic pathways: The SRY-SOX9-FGF9 network promotes testis determination partly by opposing proovarian pathways, while RSPO1/WNT-β-catenin/FOXL2 signals control ovary development by inhibiting SRY-SOX9-FGF9. The molecular basis of this mutual antagonism is unclear. Here we show that ZNRF3, a WNT signaling antagonist and direct target of RSPO1-mediated inhibition, is required for sex determination in mice. XY mice lacking ZNRF3 exhibit complete or partial gonadal sex reversal, or related defects. These abnormalities are associated with ectopic WNT/β-catenin activity and reduced Sox9 expression during fetal sex determination. Using exome sequencing of individuals with 46,XY disorders of sex development, we identified three human ZNRF3 variants in very rare cases of XY female presentation. We tested two missense variants and show that these disrupt ZNRF3 activity in both human cell lines and zebrafish embryo assays. Our data identify a testis-determining function for ZNRF3 and indicate a mechanism of direct molecular interaction between two mutually antagonistic organogenetic pathways.

Cell adhesion molecules expression pattern indicates that somatic cells arbitrate gonadal sex of differentiating bipotential fetal mouse gonad

Unlike other organ anlagens, the primordial gonad is sexually bipotential in all animals. In mouse, the bipotential gonad differentiates into testis or ovary depending on the genetic sex (XY or XX) of the fetus. During gonad development cells segregate, depending on genetic sex, into distinct compartments: testis cords and interstitium form in XY gonad, and germ cell cysts and stroma in XX gonad. However, our knowledge of mechanisms governing gonadal sex differentiation remains very vague. Because it is known that adhesion molecules (CAMs) play a key role in organogenesis, we suspected that diversified expression of CAMs should also play a crucial role in gonad development. Using microarray analysis we identified 129 CAMs and factors regulating cell adhesion during sexual differentiation of mouse gonad. To identify genes expressed differentially in three cell lines in XY and XX gonads: i) supporting (Sertoli or follicular cells), ii) interstitial or stromal cells, and iii) germ cells, we used transgenic mice expressing EGFP reporter gene and FACS cell sorting. Although a large number of CAMs expressed ubiquitously, expression of certain genes was cell line- and genetic sex-specific. The sets of CAMs differentially expressed in supporting versus interstitial/stromal cells may be responsible for segregation of these two cell lines during gonadal development. There was also a significant difference in CAMs expression pattern between XY supporting (Sertoli) and XX supporting (follicular) cells but not between XY and XX germ cells. This indicates that differential CAMs expression pattern in the somatic cells but not in the germ line arbitrates structural organization of gonadal anlagen into testis or ovary.

Ex vivo cultures combined with vivo-morpholino induced gene knockdown provide a system to assess the role of WT1 and GATA4 during gonad differentiation

Gonad morphogenesis relies on the correct spatiotemporal expression of a number of genes that together fulfill the differentiation of the bipotential gonad into testes or ovaries. As such, the transcription factors WT1 and GATA4 are pivotal for proper gonadal development. Here we address the contributions of GATA4 and WT1 to the sex differentiation phase in testes and ovaries. We applied an ex vivo technique for cultivating gonads in hanging droplets of media that were supplemented with vivo-morpholinos to knockdown WT1 and GATA4 either alone or in combination at the same developmental stage. We show that WT1 is equally important for both, the initial establishment and the maintenance of the sex-specific gene expression signature in testes and ovaries. We further identified Foxl2 as a novel putative downstream target gene of WT1. Moreover, knockdown of WT1 reduced mRNA levels of several molecular components of the hedgehog signaling pathway in XY gonads, whereas Gata4 vivo-morpholino treatment increased transcripts of Dhh and Ptch1 in embryonic testes. The data suggest that for its proper function, WT1 relies on the correct expression of the GATA4 protein. Furthermore, GATA4 down-regulates several ovarian promoting genes in testes, such as Ctnnb1, Fst, and Bmp2, suggesting that this repression is required for maintaining the male phenotype. In conclusion, this study provides novel insights into the role of WT1 and GATA4 during the sex differentiation phase and represents an approach that can be applied to assess other proteins with as yet unknown functions during gonadal development.

A genomic atlas of human adrenal and gonad development

BACKGROUND

In humans, the adrenal glands and gonads undergo distinct biological events between 6-10 weeks post conception (wpc), such as testis determination, the onset of steroidogenesis and primordial germ cell development. However, relatively little is currently known about the genetic mechanisms underlying these processes. We therefore aimed to generate a detailed genomic atlas of adrenal and gonad development across these critical stages of human embryonic and fetal development.

METHODS

RNA was extracted from 53 tissue samples between 6-10 wpc (adrenal, testis, ovary and control). Affymetrix array analysis was performed and differential gene expression was analysed using Bioconductor. A mathematical model was constructed to investigate time-series changes across the dataset. Pathway analysis was performed using ClueGo and cellular localisation of novel factors confirmed using immunohistochemistry.

RESULTS

Using this approach, we have identified novel components of adrenal development (e.g. ASB4, NPR3) and confirmed the role of SRY as the main human testis-determining gene. By mathematical modelling time-series data we have found new genes up-regulated with SOX9 in the testis (e.g. CITED1), which may represent components of the testis development pathway. We have shown that testicular steroidogenesis has a distinct onset at around 8 wpc and identified potential novel components in adrenal and testicular steroidogenesis (e.g. MGARP, FOXO4, MAP3K15, GRAMD1B, RMND2), as well as testis biomarkers (e.g. SCUBE1). We have also shown that the developing human ovary expresses distinct subsets of genes (e.g. OR10G9, OR4D5), but enrichment for established biological pathways is limited.

CONCLUSION

This genomic atlas is revealing important novel aspects of human development and new candidate genes for adrenal and reproductive disorders.

A genomic atlas of human adrenal and gonad development

Background: In humans, the adrenal glands and gonads undergo distinct biological events between 6-10 weeks post conception (wpc), such as testis determination, the onset of steroidogenesis and primordial germ cell development. However, relatively little is currently known about the genetic mechanisms underlying these processes. We therefore aimed to generate a detailed genomic atlas of adrenal and gonad development across these critical stages of human embryonic and fetal development.

Methods: RNA was extracted from 53 tissue samples between 6-10 wpc (adrenal, testis, ovary and control). Affymetrix array analysis was performed and differential gene expression was analysed using Bioconductor. A mathematical model was constructed to investigate time-series changes across the dataset. Pathway analysis was performed using ClueGo and cellular localisation of novel factors confirmed using immunohistochemistry.

Results: Using this approach, we have identified novel components of adrenal development (e.g. ASB4, NPR3) and confirmed the role of SRY as the main human testis-determining gene. By mathematical modelling time-series data we have found new genes up-regulated with SOX9 in the testis (e.g. CITED1), which may represent components of the testis development pathway. We have shown that testicular steroidogenesis has a distinct onset at around 8 wpc and identified potential novel components in adrenal and testicular steroidogenesis (e.g. MGARP, FOXO4, MAP3K15, GRAMD1B, RMND2), as well as testis biomarkers (e.g. SCUBE1). We have also shown that the developing human ovary expresses distinct subsets of genes (e.g. OR10G9, OR4D5), but enrichment for established biological pathways is limited.

Conclusion: This genomic atlas is revealing important novel aspects of human development and new candidate genes for adrenal and reproductive disorders.

Leveraging Online Resources to Prioritize Candidate Genes for Functional Analyses: Using the Fetal Testis as a Test Case

With each new microarray or RNA-seq experiment, massive quantities of transcriptomic information are generated with the purpose to produce a list of candidate genes for functional analyses. Yet an effective strategy remains elusive to prioritize the genes on these candidate lists. In this review, we outline a prioritizing strategy by taking a step back from the bench and leveraging the rich range of public databases. This in silico approach provides an economical, less biased, and more effective solution. We discuss the publicly available online resources that can be used to answer a range of questions about a gene. Is the gene of interest expressed in the system of interest (using expression databases)? Where else is this gene expressed (using added-value transcriptomic resources)? What pathways and processes is the gene involved in (using enriched gene pathway analysis and mouse knockout databases)? Is this gene correlated with human diseases (using human disease variant databases)? Using mouse fetal testis as an example, our strategies identified 298 genes annotated as expressed in the fetal testis. We cross-referenced these genes to existing microarray data and narrowed the list down to cell-type-specific candidates (35 for Sertoli cells, 11 for Leydig cells, and 25 for germ cells). Our strategies can be customized so that they allow researchers to effectively and confidently prioritize genes for functional analysis.

The Battle of the Sexes: Human Sex Development and Its Disorders

The process of sexual differentiation is central for reproduction of almost all metazoan and therefore for maintenance of practically all multicellular organisms. In sex development we can distinguish two different processes: First, sex determination is the developmental decision that directs the undifferentiated embryo into a sexually dimorphic individual. In mammals, sex determination equals gonadal development. The second process known as sex differentiation takes place once the sex determination decision has been made through factors produced by the gonads that determine the development of the phenotypic sex. Most of the knowledge on the factors involved in sexual development came from animal models and from studies of cases in whom the genetic or the gonadal sex does not match the phenotypical sex, i.e., patients affected by disorders of sex development (DSD). Generally speaking, factors influencing sex determination are transcriptional regulators, whereas factors important for sex differentiation are secreted hormones and their receptors. This review focuses on the factors involved in gonadal determination, and whenever possible, references on the "prismatic" clinical cases are given.

The Sertoli cell: one hundred fifty years of beauty and plasticity

It has been one and a half centuries since Enrico Sertoli published the seminal discovery of the testicular 'nurse cell', not only a key cell in the testis, but indeed one of the most amazing cells in the vertebrate body. In this review, we begin by examining the three phases of morphological research that have occurred in the study of Sertoli cells, because microscopic anatomy was essentially the only scientific discipline available for about the first 75 years after the discovery. Biochemistry and molecular biology then changed all of biological sciences, including our understanding of the functions of Sertoli cells. Immunology and stem cell biology were not even topics of science in 1865, but they have now become major issues in our appreciation of Sertoli cell's role in spermatogenesis. We end with the universal importance and plasticity of function by comparing Sertoli cells in fish, amphibians, and mammals. In these various classes of vertebrates, Sertoli cells have quite different modes of proliferation and epithelial maintenance, cystic vs. tubular formation, yet accomplish essentially the same function but in strikingly different ways.

Methods for the Study of Gonadal Development

Current knowledge on gonadal development and sex determination is the product of many decades of research involving a variety of scientific methods from different biological disciplines such as histology, genetics, biochemistry, and molecular biology. The earliest embryological investigations, followed by the invention of microscopy and staining methods, were based on histological examinations. The most robust development of histological staining techniques occurred in the second half of the nineteenth century and resulted in structural descriptions of gonadogenesis. These first studies on gonadal development were conducted on domesticated animals; however, currently the mouse is the most extensively studied species. The next key point in the study of gonadogenesis was the advancement of methods allowing for the in vitro culture of fetal gonads. For instance, this led to the description of the origin of cell lines forming the gonads. Protein detection using antibodies and immunolabeling methods and the use of reporter genes were also invaluable for developmental studies, enabling the visualization of the formation of gonadal structure. Recently, genetic and molecular biology techniques, especially gene expression analysis, have revolutionized studies on gonadogenesis and have provided insight into the molecular mechanisms that govern this process. The successive invention of new methods is reflected in the progress of research on gonadal development.

A Boolean network model of human gonadal sex determination

BACKGROUND

Gonadal sex determination (GSD) in humans is a complex biological process that takes place in early stages of embryonic development when the bipotential gonadal primordium (BGP) differentiates towards testes or ovaries. This decision is directed by one of two distinct pathways embedded in a GSD network activated in a population of coelomic epithelial cells, the Sertoli progenitor cells (SPC) and the granulosa progenitor cells (GPC). In males, the pathway is activated when the Sex-Determining Region Y (SRY) gene starts to be expressed, whereas in females the WNT4/ β-catenin pathway promotes the differentiation of the GPCs towards ovaries. The interactions and dynamics of the elements that constitute the GSD network are poorly understood, thus our group is interested in inferring the general architecture of this network as well as modeling the dynamic behavior of a set of genes associated to this process under wild-type and mutant conditions.

METHODS

We reconstructed the regulatory network of GSD with a set of genes directly associated with the process of differentiation from SPC and GPC towards Sertoli and granulosa cells, respectively. These genes are experimentally well-characterized and the effects of their deficiency have been clinically reported. We modeled this GSD network as a synchronous Boolean network model (BNM) and characterized its attractors under wild-type and mutant conditions.

RESULTS

Three attractors with a clear biological meaning were found; one of them corresponding to the currently known gene expression pattern of Sertoli cells, the second correlating to the granulosa cells and, the third resembling a disgenetic gonad.

CONCLUSIONS

The BNM of GSD that we present summarizes the experimental data on the pathways for Sertoli and granulosa establishment and sheds light on the overall behavior of a population of cells that differentiate within the developing gonad. With this model we propose a set of regulatory interactions needed to activate either the SRY or the WNT4/ β-catenin pathway as well as their downstream targets, which are critical for further sex differentiation. In addition, we observed a pattern of altered regulatory interactions and their dynamics that lead to some disorders of sex development (DSD).

Identification of candidate gonadal sex differentiation genes in the chicken embryo using RNA-seq

BACKGROUND

Despite some advances in recent years, the genetic control of gonadal sex differentiation during embryogenesis is still not completely understood. To identify new candidate genes involved in ovary and testis development, RNA-seq was used to define the transcriptome of embryonic chicken gonads at the onset of sexual differentiation (day 6.0/stage 29).

RESULTS

RNA-seq revealed more than 1000 genes that were transcribed in a sex-biased manner at this early stage of gonadal differentiation. Comparison with undifferentiated gonads revealed that sex biased expression was derived primarily from autosomal rather than sex-linked genes. Gene ontology and pathway analysis indicated that many of these genes encoded proteins involved in extracellular matrix function and cytoskeletal remodelling, as well as tubulogenesis. Several of these genes are novel candidate regulators of gonadal sex differentiation, based on sex-biased expression profiles that are altered following experimental sex reversal. We further characterised three female-biased (ovarian) genes; calpain-5 (CAPN5), G-protein coupled receptor 56 (GPR56), and FGFR3 (fibroblast growth factor receptor 3). Protein expression of these candidates in the developing ovaries suggests that they play an important role in this tissue.

CONCLUSIONS

This study provides insight into the earliest steps of vertebrate gonad sex differentiation, and identifies novel candidate genes for ovarian and testicular development.

Rapid screening of gene function by systemic delivery of morpholino oligonucleotides to live mouse embryos

Traditional gene targeting methods in mice are complex and time consuming, especially when conditional deletion methods are required. Here, we describe a novel technique for assessing gene function by injection of modified antisense morpholino oligonucleotides (MOs) into the heart of mid-gestation mouse embryos. After allowing MOs to circulate through the embryonic vasculature, target tissues were explanted, cultured and analysed for expression of key markers. We established proof-of-principle by partially phenocopying known gene knockout phenotypes in the fetal gonads (Stra8, Sox9) and pancreas (Sox9). We also generated a novel double knockdown of Gli1 and Gli2, revealing defects in Leydig cell differentiation in the fetal testis. Finally, we gained insight into the roles of Adamts19 and Ctrb1, genes of unknown function in sex determination and gonadal development. These studies reveal the utility of this method as a means of first-pass analysis of gene function during organogenesis before committing to detailed genetic analysis.

Genetic regulation of mammalian gonad development

Sex-specific gonadal development starts with formation of the bipotential gonad, which then differentiates into either a mature testis or an ovary. This process is dependent on activation of either the testis-specific or the ovary-specific pathway while the opposite pathway is continuously repressed. A network of transcription factors tightly regulates initiation and maintenance of these distinct pathways; disruption of these networks can lead to disorders of sex development in humans and male-to-female or female-to-male sex reversal in mice. Sry is the Y-linked master switch that is both required and sufficient to drive the testis-determining pathway. Another key component of the testis pathway is Sox9, which acts immediately downstream of Sry. In contrast to the testis pathway, no single sex-determining factor has been identified in the ovary pathway; however, multiple genes, such as Foxl2, Rspo1, Ctnnb1, and Wnt4, seem to work synergistically and in parallel to ensure proper ovary development. Our understanding of the regulatory networks that underpin testis and ovary development has grown substantially over the past two decades.

Long-term effects of early-life exposure to environmental oestrogens on ovarian function: role of epigenetics

Oestrogens play an important role in development and function of the brain and reproductive tract. Accordingly, it is considered that developmental exposure to environmental oestrogens can disrupt neural and reproductive tract development, potentially resulting in long-term alterations in neurobehaviour and reproductive function. Many chemicals have been shown to have oestrogenic activity, whereas others affect oestrogen production and turnover, resulting in the disruption of oestrogen signalling pathways. However, these mechanisms and the concentrations required to induce these effects cannot account for the myriad adverse effects of environmental toxicants on oestrogen-sensitive target tissues. Hence, alternative mechanisms are assumed to underlie the adverse effects documented in experimental animal models and thus could be important to human health. In this review, the epigenetic regulation of gene expression is explored as a potential target of environmental toxicants including oestrogenic chemicals. We suggest that toxicant-induced changes in epigenetic signatures are important mechanisms underlying the disruption of ovarian follicular development. In addition, we discuss how exposure to environmental oestrogens during early life can alter gene expression through effects on epigenetic control potentially leading to permanent changes in ovarian physiology.

Genome-wide identification of AR-regulated genes translated in Sertoli cells in vivo using the RiboTag approach

An understanding of the molecular mechanisms by which androgens drive spermatogenesis has been thwarted by the fact that few consistent androgen receptor (AR) target genes have been identified. Here, we addressed this issue using next-generation sequencing coupled with the RiboTag approach, which purifies translated mRNAs expressed in cells that express cyclic recombinase (CRE). Using RiboTag mice expressing CRE in Sertoli cells (SCs), we identified genes expressed specifically in SCs in both prepubertal and adult mice. Unexpectedly, this analysis revealed that the SC-specific gene program is already largely defined at the initiation of spermatogenesis despite the subsequent dramatic maturational changes known to occur in SCs. To identify AR-regulated genes, we generated triple-mutant mice in which the SCs express the RiboTag but lack ARs. RNA sequencing analysis revealed hundreds of SC-expressed AR-regulated genes that had previously gone unnoticed, including suppressed genes involved in ovarian development. Comparison of the SC-enriched dataset with that from the whole testes allowed us to classify genes in terms of their degree of expression in SCs. This revealed that a greater fraction of AR-up-regulated genes than AR-down-regulated genes were expressed predominantly in SCs. Our results also revealed that AR signaling in SCs causes a large number of genes not detectably expressed in SCs to undergo altered expression, thereby providing genome-wide evidence for wide-scale communication between SCs and other cells. Taken together, our results identified novel classes of genes expressed in a hormone-dependent manner in different testicular cell subsets and highlight a new approach to analyze cell type-specific gene regulation.

FOXL2 transcriptionally represses Sf1 expression by antagonizing WT1 during ovarian development in mice

Steroidogenic factor 1 (SF1; Ad4BP/NR5A1) plays key roles in gonadal development. Initially, the Sf1 gene is expressed in mouse fetal gonads of both sexes, but later is up-regulated in testes and down-regulated in ovaries. While Sf1 expression is activated and maintained by Wilms tumor 1 (WT1) and LIM homeobox 9 (LHX9), the mechanism of sex-specific regulation remains unclear. We hypothesized that Sf1 is repressed by the transcription factor Forkhead box L2 (FOXL2) during ovarian development. In an in vitro system (TM3 cells), up-regulation of Sf1 by the WT1 splice variant WT1-KTS was antagonized by FOXL2, as determined by quantitative RT-PCR. Using reporter assays, we localized the Sf1 proximal promoter region involved in this antagonism to a 674-bp interval. A conserved FOXL2 binding site was identified in this interval by in vitro chromatin immunoprecipitation. Introducing mutations into this site abolished negative regulation by FOXL2 in reporter assays. Finally, in Foxl2-null mice, Sf1 expression was increased 2-fold relative to wild-type XX fetal gonads. Our results support the hypothesis that FOXL2 negatively regulates Sf1 expression by antagonizing WT1-KTS during early ovarian development in mice.

SOX9 regulates microRNA miR-202-5p/3p expression during mouse testis differentiation

MicroRNAs are important regulators of developmental gene expression, but their contribution to fetal gonad development is not well understood. We have identified the evolutionarily conserved gonadal microRNAs miR-202-5p and miR-202-3p as having a potential role in regulating mouse embryonic gonad differentiation. These microRNAs are expressed in a sexually dimorphic pattern as the primordial XY gonad differentiates into a testis, with strong expression in Sertoli cells. In vivo, ectopic expression of pri-miR-202 in XX gonads did not result in molecular changes to the ovarian determination pathway. Expression of the primary transcript of miR-202-5p/3p remained low in XY gonads in a conditional Sox9-null mouse model, suggesting that pri-miR-202 transcription is downstream of SOX9, a transcription factor that is both necessary and sufficient for male sex determination. We identified the pri-miR-202 promoter that is sufficient to drive expression in XY but not XX fetal gonads ex vivo. Mutation of SOX9 and SF1 binding sites reduced ex vivo transactivation of the pri-miR-202 promoter, demonstrating that pri-miR-202 may be a direct transcriptional target of SOX9/SF1 during testis differentiation. Our findings indicate that expression of the conserved gonad microRNA, miR-202-5p/3p, is downstream of the testis-determining factor SOX9, suggesting an early role in testis development.

Translational genetics for diagnosis of human disorders of sex development

Disorders of sex development (DSDs) are congenital conditions with discrepancies between the chromosomal, gonadal, and phenotypic sex of the individual. Such disorders have historically been difficult to diagnose and cause great stress to patients and their families. Genetic analysis of human samples has been instrumental in elucidating the molecules and pathways involved in the development of the bipotential gonad into a functioning testis or ovary. However, many DSD patients still do not receive a genetic diagnosis. New genetic and genomic technologies are expanding our knowledge of the underlying mechanism of DSDs and opening new avenues for clinical diagnosis. We review the genetic technologies that have elucidated the genes that are well established in sex determination in humans, discuss findings from more recent genomic technologies, and propose a new paradigm for clinical diagnosis of DSDs.

Fine time course expression analysis identifies cascades of activation and repression and maps a putative regulator of mammalian sex determination

In vertebrates, primary sex determination refers to the decision within a bipotential organ precursor to differentiate as a testis or ovary. Bifurcation of organ fate begins between embryonic day (E) 11.0-E12.0 in mice and likely involves a dynamic transcription network that is poorly understood. To elucidate the first steps of sexual fate specification, we profiled the XX and XY gonad transcriptomes at fine granularity during this period and resolved cascades of gene activation and repression. C57BL/6J (B6) XY gonads showed a consistent ~5-hour delay in the activation of most male pathway genes and repression of female pathway genes relative to 129S1/SvImJ, which likely explains the sensitivity of the B6 strain to male-to-female sex reversal. Using this fine time course data, we predicted novel regulatory genes underlying expression QTLs (eQTLs) mapped in a previous study. To test predictions, we developed an in vitro gonad primary cell assay and optimized a lentivirus-based shRNA delivery method to silence candidate genes and quantify effects on putative targets. We provide strong evidence that Lmo4 (Lim-domain only 4) is a novel regulator of sex determination upstream of SF1 (Nr5a1), Sox9, Fgf9, and Col9a3. This approach can be readily applied to identify regulatory interactions in other systems.

NOTCH signaling in Sertoli cells regulates gonocyte fate

Gonocytes (or prospermatogonia) are the precursors to spermatogonial stem cells (SSCs), which provide the foundation for spermatogenesis through their ability to both self-renew and generate daughter cells. Despite their relative importance, the regulatory mechanisms that govern gonocyte maintenance and transition to SSCs are poorly understood. Recently, we reported that constitutive activation of NOTCH1 signaling in Sertoli cells causes gonocyte exit from quiescence--the first suggestion of the potential role of this signaling pathway in the testis. This Extra View will review what is known about NOTCH signaling, particularly in Sertoli cells and germ cells in the testes, by providing a background on germ cell biology and a summary of our recently published data on NOTCH1 signaling in Sertoli cells. We also describe additional data showing that aberrant proliferation and differentiation of gonocytes in response to constitutive activation of NOTCH1 signaling in Sertoli cells involves de novo expression of cell cycle proteins and a marked upregulation of the KIT receptor. These data further suggest that NOTCH signaling orchestrates a dynamic balance between maintenance and differentiation of gonocytes in the perinatal testis.

Gadd45g is essential for primary sex determination, male fertility and testis development

In humans and most mammals, differentiation of the embryonic gonad into ovaries or testes is controlled by the Y-linked gene SRY. Here we show a role for the Gadd45g protein in this primary sex differentiation. We characterized mice deficient in Gadd45a, Gadd45b and Gadd45g, as well as double-knockout mice for Gadd45ab, Gadd45ag and Gadd45bg, and found a specific role for Gadd45g in male fertility and testis development. Gadd45g-deficient XY mice on a mixed 129/C57BL/6 background showed varying degrees of disorders of sexual development (DSD), ranging from male infertility to an intersex phenotype or complete gonadal dysgenesis (CGD). On a pure C57BL/6 (B6) background, all Gadd45g(-/-) XY mice were born as completely sex-reversed XY-females, whereas lack of Gadd45a and/or Gadd45b did not affect primary sex determination or testis development. Gadd45g expression was similar in female and male embryonic gonads, and peaked around the time of sex differentiation at 11.5 days post-coitum (dpc). The molecular cause of the sex reversal was the failure of Gadd45g(-/-) XY gonads to achieve the SRY expression threshold necessary for testes differentiation, resulting in ovary and Müllerian duct development. These results identify Gadd45g as a candidate gene for male infertility and 46,XY sex reversal in humans.

Disorders of sex development: new genes, new concepts

Formerly known as 'intersex' conditions, disorders of sex development (DSDs) are congenital conditions in which chromosomal, gonadal or anatomical sex is atypical. A complete revision of the nomenclature and classification of DSDs has been undertaken, which emphasizes the genetic aetiology of these disorders and discards pejorative terms. Uptake of the new terminology is widespread. DSDs affecting gonadal development are perhaps the least well understood. Unravelling the molecular mechanisms underlying gonadal development has revealed new causes of DSDs, although a specific molecular diagnosis is made in only ∼20% of patients. Conversely, identification of the molecular causes of DSDs has provided insight into the mechanisms of gonadal development. Studies of N-ethyl-N-nitrosourea mutagenesis in the mouse, and multigene diagnostic screening and genome-wide approaches, such as array-comparative genomic hybridization and next-generation sequencing, in patients with DSDs are accelerating the discovery of genes involved in gonadal development and DSDs. Furthermore, long-range gene regulatory mutations and multiple gene mutations are emerging as new causes of DSDs. Patients with DSDs, their parents and medical staff are confronted with challenging decisions regarding gender assignment, genital surgery and lifelong care. These advances are refining prognostic prediction and systematically improving the diagnosis and long-term management of children with DSDs.

A complete gonadal dysgenesis case with mental retardation, congenital hip dislocation, severe vertebra rotoscoliosis, pectus excavatus, and spina bifida occulta

BACKGROUND

46,XY, or Swyer syndrome, is a complete gonadal dysgenesis. Patients usually presents with primary amenorrhea with underdeveloped secondary sex characteristics. Phenotypes of these patients are female. In this report, a Swyer syndrome case is reported with novel clinical features that are classified as connective tissue disorders. This case and the 2 other previously reported Swyer syndrome cases with ascendant aortic aneurysm and diaphragmatic hernia are suggest that the Y chromosome has an important role in the structure of connective tissue.

CASE

Here we report a case of a 17-year-old with clinical features of 46,XY complete gonadal dysgenesis including external female genitalia, hypoplastic uterus, hypergonadotrophic hypogonadism, incomplete secondary sex characterics, primary amenorrhea, and normal male karyotype. In addition, she had mild mental retardation, severe rotoscoliosis, pectus excavatus, spina bifida occulta, hip dislocation, and long, slender extremities. She had a rudimentary uterus and streak gonads; after giving her cyclic estrogen and progesterone pills, she was able to menstruate.

SUMMARY AND CONCLUSION

In this report, a Swyer syndrome case was discussed regarding clinical features, especially those are not characteristic for Swyer syndrome after a review of the literature.

Insulin and IGF1 receptors are essential for XX and XY gonadal differentiation and adrenal development in mice

Mouse sex determination provides an attractive model to study how regulatory genetic networks and signaling pathways control cell specification and cell fate decisions. This study characterizes in detail the essential role played by the insulin receptor (INSR) and the IGF type I receptor (IGF1R) in adrenogenital development and primary sex determination. Constitutive ablation of insulin/IGF signaling pathway led to reduced proliferation rate of somatic progenitor cells in both XX and XY gonads prior to sex determination together with the downregulation of hundreds of genes associated with the adrenal, testicular, and ovarian genetic programs. These findings indicate that prior to sex determination somatic progenitors in Insr;Igf1r mutant gonads are not lineage primed and thus incapable of upregulating/repressing the male and female genetic programs required for cell fate restriction. In consequence, embryos lacking functional insulin/IGF signaling exhibit (i) complete agenesis of the adrenal cortex, (ii) embryonic XY gonadal sex reversal, with a delay of Sry upregulation and the subsequent failure of the testicular genetic program, and (iii) a delay in ovarian differentiation so that Insr;Igf1r mutant gonads, irrespective of genetic sex, remained in an extended undifferentiated state, before the ovarian differentiation program ultimately is initiated at around E16.5.