Elucidating the Transcriptional States of Spermatogenesis-Joint Analysis of Germline and Supporting Cell, Mice and Human, Normal and Perturbed, Bulk and Single-Cell RNA-Seq

In studying the molecular underpinning of spermatogenesis, we expect to understand the fundamental biological processes better and potentially identify genes that may lead to novel diagnostic and therapeutic strategies toward precision medicine in male infertility. In this review, we emphasized our perspective that the path forward necessitates integrative studies that rely on complementary approaches and types of data. To comprehensively analyze spermatogenesis, this review proposes four axes of integration. First, spanning the analysis of spermatogenesis in the healthy state alongside pathologies. Second, the experimental analysis of model systems (in which we can deploy treatments and perturbations) alongside human data. Third, the phenotype is measured alongside its underlying molecular profiles using known markers augmented with unbiased profiles. Finally, the testicular cells are studied as ecosystems, analyzing the germ cells alongside the states observed in the supporting somatic cells. Recently, the study of spermatogenesis has been advancing using single-cell RNA sequencing, where scientists have uncovered the unique stages of germ cell development in mice, revealing new regulators of spermatogenesis and previously unknown cell subtypes in the testis. An in-depth analysis of meiotic and postmeiotic stages led to the discovery of marker genes for spermatogonia, Sertoli and Leydig cells and further elucidated all the other germline and somatic cells in the testis microenvironment in normal and pathogenic conditions. The outcome of an integrative analysis of spermatogenesis using advanced molecular profiling technologies such as scRNA-seq has already propelled our biological understanding, with additional studies expected to have clinical implications for the study of male fertility. By uncovering new genes and pathways involved in abnormal spermatogenesis, we may gain insights into subfertility or sterility.

Differential gene expression and immune cell infiltration in maedi-visna virus-infected lung tissues

BACKGROUND

Maedi-visna virus (MVV) is a lentivirus that infects monocyte/macrophage lineage cells in sheep, goats, and wild ruminants and causes pneumonia, mastitis, arthritis, and encephalitis. The immune response to MVV infection is complex, and a complete understanding of its infection and pathogenesis is lacking. This study investigated the in vivo transcriptomic patterns of lung tissues in sheep exposed to MVV using the RNA sequencing technology.

RESULT

The results indicated that 2,739 genes were significantly differentially expressed, with 1,643 downregulated genes and 1,096 upregulated genes. Many variables that could be unique to MVV infections were discovered. Gene Ontology analysis revealed that a significant proportion of genes was enriched in terms directly related to the immune system and biological responses to viral infections. Kyoto Encyclopedia of Genes and Genomes analysis revealed that the most enriched pathways were related to virus-host cell interactions and inflammatory responses. Numerous immune-related genes, including those encoding several cytokines and interferon regulatory factors, were identified in the protein-protein interaction network of differentially expressed genes (DEGs). The expression of DEGs was evaluated using real-time polymerase chain reaction and western blot analysis. CXCL13, CXCL6, CXCL11, CCR1, CXCL8, CXCL9, CXCL10, TNFSF8, TNFRSF8, IL7R, IFN-γ, CCL2, and MMP9 were upregulated. Immunohistochemical analysis was performed to identify the types of immune cells that infiltrated MVV-infected tissues. B cells, CD4+ and CD8+ T cells, and macrophages were the most prevalent immune cells correlated with MVV infection in the lungs.

CONCLUSION

Overall, the findings of this study provide a comprehensive understanding of the in vivo host response to MVV infection and offer new perspectives on the gene regulatory networks that underlie pathogenesis in natural hosts.

Direct male development in chromosomally ZZ zebrafish

The genetics of sex determination varies across taxa, sometimes even within a species. Major domesticated strains of zebrafish (Danio rerio), including AB and TU, lack a strong genetic sex determining locus, but strains more recently derived from nature, like Nadia (NA), possess a ZZ male/ZW female chromosomal sex-determination system. AB fish pass through a juvenile ovary stage, forming oocytes that survive in fish that become females but die in fish that become males. To understand mechanisms of gonad development in NA zebrafish, we studied histology and single cell transcriptomics in developing ZZ and ZW fish. ZW fish developed oocytes by 22 days post-fertilization (dpf) but ZZ fish directly formed testes, avoiding a juvenile ovary phase. Gonads of some ZW and WW fish, however, developed oocytes that died as the gonad became a testis, mimicking AB fish, suggesting that the gynogenetically derived AB strain is chromosomally WW. Single-cell RNA-seq of 19dpf gonads showed similar cell types in ZZ and ZW fish, including germ cells, precursors of gonadal support cells, steroidogenic cells, interstitial/stromal cells, and immune cells, consistent with a bipotential juvenile gonad. In contrast, scRNA-seq of 30dpf gonads revealed that cells in ZZ gonads had transcriptomes characteristic of testicular Sertoli, Leydig, and germ cells while ZW gonads had granulosa cells, theca cells, and developing oocytes. Hematopoietic and vascular cells were similar in both sex genotypes. These results show that juvenile NA zebrafish initially develop a bipotential gonad; that a factor on the NA W chromosome, or fewer than two Z chromosomes, is essential to initiate oocyte development; and without the W factor, or with two Z doses, NA gonads develop directly into testes without passing through the juvenile ovary stage. Sex determination in AB and TU strains mimics NA ZW and WW zebrafish, suggesting loss of the Z chromosome during domestication. Genetic analysis of the NA strain will facilitate our understanding of the evolution of sex determination mechanisms.

Fetal Leydig cells: What we know and what we don't

During male fetal development, testosterone plays an essential role in the differentiation and maturation of the male reproductive system. Deficient fetal testosterone production can result in variations of sex differentiation that may cause infertility and even increased tumor incidence later in life. Fetal Leydig cells in the fetal testis are the major androgen source in mammals. Although fetal and adult Leydig cells are similar in their functions, they are two distinct cell types, and therefore, the knowledge of adult Leydig cells cannot be directly applied to understanding fetal Leydig cells. This review summarizes our current knowledge of fetal Leydig cells regarding their cell biology, developmental biology, and androgen production regulation in rodents and human. Fetal Leydig cells are present in basement membrane-enclosed clusters in between testis cords. They originate from the mesonephros mesenchyme and the coelomic epithelium and start to differentiate upon receiving a Desert Hedgehog signal from Sertoli cells or being released from a NOTCH signal from endothelial cells. Mature fetal Leydig cells produce androgens. Human fetal Leydig cell steroidogenesis is LHCGR (Luteinizing Hormone Chronic Gonadotropin Receptor) dependent, while rodents are not, although other Gαs -protein coupled receptors might be involved in rodent steroidogenesis regulation. Fetal steroidogenesis ceases after sex differentiation is completed, and some fetal Leydig cells dedifferentiate to serve as stem cells for adult testicular cell types. Significant gaps are acknowledged: (1) Why are adult and fetal Leydig cells different? (2) What are bona fide progenitor and fetal Leydig cell markers? (3) Which signaling pathways and transcription factors regulate fetal Leydig cell steroidogenesis? It is critical to discover answers to these questions so that we can understand vulnerable targets in fetal Leydig cells and the mechanisms for androgen production that when disrupted, leads to variations in sex differentiation that range from subtle to complete sex reversal.

Direct Male Development in Chromosomally ZZ Zebrafish

The genetics of sex determination varies across taxa, sometimes even within a species. Major domesticated strains of zebrafish ( Danio rerio ), including AB and TU, lack a strong genetic sex determining locus, but strains more recently derived from nature, like Nadia (NA), possess a ZZ male/ZW female chromosomal sex-determination system. AB strain fish pass through a juvenile ovary stage, forming oocytes that survive in fish that become females but die in fish that become males. To understand mechanisms of gonad development in NA zebrafish, we studied histology and single cell transcriptomics in developing ZZ and ZW fish. ZW fish developed oocytes by 22 days post-fertilization (dpf) but ZZ fish directly formed testes, avoiding a juvenile ovary phase. Gonads of some ZW and WW fish, however, developed oocytes that died as the gonad became a testis, mimicking AB fish, suggesting that the gynogenetically derived AB strain is chromosomally WW. Single-cell RNA-seq of 19dpf gonads showed similar cell types in ZZ and ZW fish, including germ cells, precursors of gonadal support cells, steroidogenic cells, interstitial/stromal cells, and immune cells, consistent with a bipotential juvenile gonad. In contrast, scRNA-seq of 30dpf gonads revealed that cells in ZZ gonads had transcriptomes characteristic of testicular Sertoli, Leydig, and germ cells while ZW gonads had granulosa cells, theca cells, and developing oocytes. Hematopoietic and vascular cells were similar in both sex genotypes. These results show that juvenile NA zebrafish initially develop a bipotential gonad; that a factor on the NA W chromosome or fewer than two Z chromosomes is essential to initiate oocyte development; and without the W factor or with two Z doses, NA gonads develop directly into testes without passing through the juvenile ovary stage. Sex determination in AB and TU strains mimics NA ZW and WW zebrafish, suggesting loss of the Z chromosome during domestication. Genetic analysis of the NA strain will facilitate our understanding of the evolution of sex determination mechanisms.

Crucial Roles of the Mesenchymal Androgen Receptor in Wolffian Duct Development

Wolffian duct (WD) maintenance and differentiation is predominantly driven by the androgen action, which is mediated by the androgen receptor (AR). It is well established that the mesenchyme indicates the fate and differentiation of epithelial cells. However, in vivo developmental requirement of mesenchymal AR in WD development is still undefined. By designing a mesenchyme-specific Ar knockout (ARcKO), we discovered that the loss of mesenchymal Ar led to the bilateral or unilateral degeneration of caudal WDs and cystic formation at the cranial WDs. Ex vivo culture of ARcKO WDs invariably resulted in bilateral defects, suggesting that some factor(s) originating from surrounding tissues in vivo might promote WD survival and growth even in the absence of mesenchymal Ar. Mechanistically, we found cell proliferation was significantly reduced in both epithelial and mesenchymal compartments; but cell apoptosis was not affected. Transcriptomic analysis by RNA sequencing of E14.5 mesonephroi revealed 131 differentially expressed genes. Multiple downregulated genes (Top2a, Wnt9b, Lama2, and Lamc2) were associated with morphological and cellular changes in ARcKO male embryos (ie, reduced cell proliferation and decreased number of epithelial cells). Mesenchymal differentiation into smooth muscle cells that are critical for morphogenesis was also impaired in ARcKO male embryos. Taken together, our results demonstrate the crucial roles of the mesenchymal AR in WD maintenance and morphogenesis in mice.

The construction of a testis transcriptional cell atlas from embryo to adult reveals various somatic cells and their molecular roles

BACKGROUND

The testis is a complex organ that undergoes extensive developmental changes from the embryonic stage to adulthood. The development of germ cells, which give rise to spermatozoa, is tightly regulated by the surrounding somatic cells.

METHODS

To better understand the dynamics of these changes, we constructed a transcriptional cell atlas of the testis, integrating single-cell RNA sequencing data from over 26,000 cells across five developmental stages: fetal germ cells, infants, childhood, peri-puberty, and adults. We employed various analytical techniques, including clustering, cell type assignments, identification of differentially expressed genes, pseudotime analysis, weighted gene co-expression network analysis, and evaluation of paracrine cell-cell communication, to comprehensively analyze this transcriptional cell atlas of the testis.

RESULTS

Our analysis revealed remarkable heterogeneity in both somatic and germ cell populations, with the highest diversity observed in Sertoli and Myoid somatic cells, as well as in spermatogonia, spermatocyte, and spermatid germ cells. We also identified key somatic cell genes, including RPL39, RPL10, RPL13A, FTH1, RPS2, and RPL18A, which were highly influential in the weighted gene co-expression network of the testis transcriptional cell atlas and have been previously implicated in male infertility. Additionally, our analysis of paracrine cell-cell communication supported specific ligand-receptor interactions involved in neuroactive, cAMP, and estrogen signaling pathways, which support the crucial role of somatic cells in regulating germ cell development.

CONCLUSIONS

Overall, our transcriptional atlas provides a comprehensive view of the cell-to-cell heterogeneity in the testis and identifies key somatic cell genes and pathways that play a central role in male fertility across developmental stages.

Emerging concepts on Leydig cell development in fetal and adult testis

Leydig cells (Lc) reside in the interstitial compartment of the testis and are the target of Luteinising hormone (LH) for Testosterone (T) production, thus critically regulates male fertility. Classical histological studies have identified two morphologically different populations of Lc during testicular development [fetal (FLc) and adult (ALc)]. Recent progress in ex vivo cell/organ culture, genome-wide analysis, genetically manipulated mouse models, lineage tracing, and single-cell RNA-seq experiments have revealed the diverse cellular origins with differential transcriptomic and distinct steroidogenic outputs of these populations. FLc originates from both coelomic epithelium and notch-active Nestin-positive perivascular cells located at the gonad-mesonephros borders, and get specified as Nr5a1 (previously known as Ad4BP/SF-1) expressing cells by embryonic age (E) 12.5 days in fetal mouse testes. These cells produce androstenedione (precursor of T, due to lack of HSD17β3 enzyme) and play critical a role in initial virilization and patterning of the male external genitalia. However, in neonatal testis, FLc undergoes massive regression/dedifferentiation and gradually gets replaced by T-producing ALc. Very recent studies suggest a small fraction (5-20%) of FLc still persists in adult testis. Both Nestin-positive perivascular cells and FLc are considered to be the progenitor populations for ALc. This minireview article summarizes the current understanding of Lc development in fetal and adult testes highlighting their common or diverse cellular (progenitor/stem) origins with respective functional significance in both rodents and primates. (227 words).

Cell adhesion function was altered during the seasonal regression of the seminiferous epithelium in the mink species Neovison vison

Minks are seasonal breeders whose seminiferous epithelium undergoes regression through massive germ cell death, leaving only Sertoli cells and spermatogonial cells in the tubules. However, the molecular mechanisms that control this biological process remain largely unknown. This study describes a transcriptomic analysis of mink testes at various reproductive stages (active, regressing, and inactive). A comparison of seminiferous epithelium at different stages of reproduction shows that cell adhesion is altered during regression. In addition, genes and proteins involved in forming the blood-testis barrier (BTB) were examined in sexually active and inactive minks. The seminiferous epithelium in the testes of sexually inactive minks expressed occludin, but this expression was not discernibly observed in the testes of sexually active minks. There was no discernible expression of CX43 in the seminiferous epithelium in the testes of sexually inactive minks, but CX43 was expressed in the testes of sexually active minks. During the regression process, we observed a remarkable increase in the expression levels of Claudin-11, which is associated with Sertoli-germ cell junctions. In conclusion, these findings suggest a loss of Sertoli-germ cell adhesion, which may regulate postmeiotic cell shedding during testicular regression in mink.

Decoding the Spermatogenesis Program: New Insights from Transcriptomic Analyses

Spermatogenesis is a complex differentiation process coordinated spatiotemporally across and along seminiferous tubules. Cellular heterogeneity has made it challenging to obtain stage-specific molecular profiles of germ and somatic cells using bulk transcriptomic analyses. This has limited our ability to understand regulation of spermatogenesis and to integrate knowledge from model organisms to humans. The recent advancement of single-cell RNA-sequencing (scRNA-seq) technologies provides insights into the cell type diversity and molecular signatures in the testis. Fine-grained cell atlases of the testis contain both known and novel cell types and define the functional states along the germ cell developmental trajectory in many species. These atlases provide a reference system for integrated interspecies comparisons to discover mechanistic parallels and to enable future studies. Despite recent advances, we currently lack high-resolution data to probe germ cell-somatic cell interactions in the tissue environment, but the use of highly multiplexed spatial analysis technologies has begun to resolve this problem. Taken together, recent single-cell studies provide an improvedunderstanding of gametogenesis to examine underlying causes of infertility and enable the development of new therapeutic interventions.

Targeted Disruption of Lats1 and Lats2 in Mice Impairs Testis Development and Alters Somatic Cell Fate

Hippo signaling plays an essential role in the development of numerous tissues. Although it was previously shown that the transcriptional effectors of Hippo signaling Yes-associated protein (YAP) and transcriptional coactivator with PDZ-binding motif (TAZ) can fine-tune the regulation of sex differentiation genes in the testes, the role of Hippo signaling in testis development remains largely unknown. To further explore the role of Hippo signaling in the testes, we conditionally deleted the key Hippo kinases large tumor suppressor homolog kinases 1 and -2 (Lats1 and Lats2, two kinases that antagonize YAP and TAZ transcriptional co-regulatory activity) in the somatic cells of the testes using an Nr5a1-cre strain (Lats1flox/flox;Lats2flox/flox;Nr5a1-cre). We report here that early stages of testis somatic cell differentiation were not affected in this model but progressive testis cord dysgenesis was observed starting at gestational day e14.5. Testis cord dysgenesis was further associated with the loss of polarity of the Sertoli cells and the loss of SOX9 expression but not WT1. In parallel with testis cord dysgenesis, a loss of steroidogenic gene expression associated with the appearance of myofibroblast-like cells in the interstitial space was also observed in mutant animals. Furthermore, the loss of YAP phosphorylation, the accumulation of nuclear TAZ (and YAP) in both the Sertoli and interstitial cell populations, and an increase in their transcriptional co-regulatory activity in the testes suggest that the observed phenotype could be attributed at least in part to YAP and TAZ. Taken together, our results suggest that Hippo signaling is required to maintain proper differentiation of testis somatic cells.

Identification of sex determination locus in sea cucumber Apostichopus japonicus using genome-wide association study

Background

Sex determination mechanisms are complicated and diverse across taxonomic categories. Sea cucumber Apostichopus japonicus is a benthic echinoderm, which is the closest group of invertebrates to chordate, and important economic and ecologically aquaculture species in China. A. japonicus is dioecious, and no phenotypic differences between males and females can be detected before sexual maturation. Identification of sex determination locus will broaden knowledge about sex-determination mechanism in echinoderms, which allows for the identification of sex-linked markers and increases the efficiency of sea cucumber breeding industry.

Results

Here, we integrated assembly of a novel chromosome-level genome and resequencing of female and male populations to investigate the sex determination mechanisms of A. japonicus. We built a chromosome-level genome assembly AJH1.0 using Hi-C technology. The assembly AJH1.0 consists of 23 chromosomes ranging from 22.4 to 60.4 Mb. To identify the sex-determination locus of A. japonicus, we conducted genome-wide association study (GWAS) and analyses of distribution characteristics of sex-specific SNPs and fixation index F_ST. The GWAS analysis showed that multiple sex-associated loci were located on several chromosomes, including chromosome 4 (24.8%), followed by chromosome 9 (10.7%), chromosome 17 (10.4%), and chromosome 18 (14.1%). Furthermore, analyzing the homozygous and heterozygous genotypes of plenty of sex-specific SNPs in females and males confirmed that A. japonicus might have a XX/XY sex determination system. As a physical region of 10 Mb on chromosome 4 included the highest number of sex-specific SNPs and higher F_ST values, this region was considered as the candidate sex determination region (SDR) in A. japonicus.

Conclusions

In the present study, we integrated genome-wide association study and analyses of sex-specific variations to investigate sex determination mechanisms. This will bring novel insights into gene regulation during primitive gonadogenesis and differentiation and identification of master sex determination gene in sea cucumber. In the sea cucumber industry, investigation of molecular mechanisms of sex determination will be helpful for artificial fertilization and precise breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08632-3.

Distinctive functioning of STARD1 in the fetal Leydig cells compared to adult Leydig and adrenal cells. Impact of Hedgehog signaling via the primary cilium

STARD1 stimulates cholesterol transfer to mitochondrial CYP11A1 for conversion to pregnenolone. A cholesterol-binding START domain is guided by an N-terminal domain in a cell selective manner. Fetal and adult Leydig cells (FLC, ALC) show distinct Stard1 regulation. sm- FISH microscopy, which resolves individual molecules of Stard1 mRNA, shows uniformly high basal expression in each FLC. In ALC, in vivo, and cultured MA-10 cells, basal Stard1 expression is minimal. PKA activates loci asynchronously, with delayed splicing/export of 3.5 kb mRNA to mitochondria. After 60 min, ALC transition to an integrated mRNA delivery to mitochondria that is seen in FLC. Sertoli cells cooperate in Stard1 stimulation in FLC by delivering DHH to the primary cilium. There PTCH, SMO and cholesterol cooperate to release GLI3 to activate the Stard1 locus, probably by directing histone changes. ALC lack cilia. PKA then primes locus activation. FLC and ALC share similar SIK/CRTC/CREB regulation characterized for adrenal cells.

Cholesterol Contributes to Male Sex Differentiation Through Its Developmental Role in Androgen Synthesis and Hedgehog Signaling

Two specialized functions of cholesterol during fetal development include serving as a precursor to androgen synthesis and supporting hedgehog (HH) signaling activity. Androgens are produced by the testes to facilitate masculinization of the fetus. Recent evidence shows that intricate interactions between the HH and androgen signaling pathways are required for optimal male sex differentiation and defects of either can cause birth anomalies indicative of 46,XY male variations of sex development (VSD). Further, perturbations in cholesterol synthesis can cause developmental defects, including VSD, that phenocopy those caused by disrupted androgen or HH signaling, highlighting the functional role of cholesterol in promoting male sex differentiation. In this review, we focus on the role of cholesterol in systemic androgen and local HH signaling events during fetal masculinization and their collective contributions to pediatric VSD.

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.

Testicular somatic cell-like cells derived from embryonic stem cells induce differentiation of epiblasts into germ cells

Regeneration of the testis from pluripotent stem cells is a real challenge, reflecting the complexity of the interaction of germ cells and somatic cells. Here we report the generation of testicular somatic cell-like cells (TesLCs) including Sertoli cell-like cells (SCLCs) from mouse embryonic stem cells (ESCs) in xeno-free culture. We find that Nr5a1/SF1 is critical for interaction between SCLCs and PGCLCs. Intriguingly, co-culture of TesLCs with epiblast-like cells (EpiLCs), rather than PGCLCs, results in self-organised aggregates, or testicular organoids. In the organoid, EpiLCs differentiate into PGCLCs or gonocyte-like cells that are enclosed within a seminiferous tubule-like structure composed of SCLCs. Furthermore, conditioned medium prepared from TesLCs has a robust inducible activity to differentiate EpiLCs into PGCLCs. Our results demonstrate conditions for in vitro reconstitution of a testicular environment from ESCs and provide further insights into the generation of sperm entirely in xeno-free culture.

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.

STARD1 Functions in Mitochondrial Cholesterol Metabolism and Nascent HDL Formation. Gene Expression and Molecular mRNA Imaging Show Novel Splicing and a 1:1 Mitochondrial Association

STARD1 moves cholesterol (CHOL) from the outer mitochondrial membrane (OMM) to the inner membrane (IMM) in steroidogenic cells. This activity is integrated into CHOL trafficking and synthesis homeostasis, involving uptake through SR-B1 and LDL receptors and distribution through endosomes, ER, and lipid droplets. In adrenal cells, STARD1 is imported into the mitochondrial matrix accompanied by delivery of several hundred CHOL molecules. This transfer limits CYP11A1-mediated generation of pregnenolone. CHOL transfer is coupled to translation of STARD1 mRNA at the OMM. In testis cells, slower CHOL trafficking seems to be limiting. STARD1 also functions in a slower process through ER OMM contacts. The START domain of STARD1 is utilized by a family of genes, which includes additional STARD (forms 3-6) and GRAMD1B proteins that transfer CHOL. STARD forms 2 and 7 deliver phosphatidylcholine. STARD1 and STARD7 target their respective activities to mitochondria, via N-terminal domains (NTD) of over 50 amino acids. The NTD is not essential for steroidogenesis but exerts tissue-selective enhancement (testis>>adrenal). Three conserved sites for cleavage by the mitochondrial processing protease (MPP) generate three forms, each potentially with specific functions, as demonstrated in STARD7. STARD1 is expressed in macrophage and cardiac repair fibroblasts. Additional functions include CHOL metabolism by CYP27A1 that directs activation of LXR and CHOL export processes. STARD1 generates 3.5- and 1.6-kb mRNA from alternative polyadenylation. The 3.5-kb form exclusively binds the PKA-induced regulator, TIS11b, which binds at conserved sites in the extended 3'UTR to control mRNA translation and turnover. STARD1 expression also exhibits a novel, slow splicing that delayed splicing delivery of mRNA to mitochondria. Stimulation of transcription by PKA is directed by suppression of SIK forms that activate a CRTC/CREB/CBP promoter complex. This process is critical to pulsatile hormonal activation in vivo. sm-FISH RNA imaging shows a flow of single STARD1 mRNA particles from asymmetric accumulations of primary transcripts at gene loci to 1:1 complex of 3.5-kb mRNA with peri-nuclear mitochondria. Adrenal cells are similar but distinguished from testis cells by appreciable basal expression prior to hormonal activation. This difference is conserved in culture and in vivo.

Patterns, Profiles, and Parsimony: Dissecting Transcriptional Signatures From Minimal Single-Cell RNA-Seq Output With SALSA

Single-cell RNA sequencing (scRNA-seq) technologies have precipitated the development of bioinformatic tools to reconstruct cell lineage specification and differentiation processes with single-cell precision. However, current start-up costs and recommended data volumes for statistical analysis remain prohibitively expensive, preventing scRNA-seq technologies from becoming mainstream. Here, we introduce single-cell amalgamation by latent semantic analysis (SALSA), a versatile workflow that combines measurement reliability metrics with latent variable extraction to infer robust expression profiles from ultra-sparse sc-RNAseq data. SALSA uses a matrix focusing approach that starts by identifying facultative genes with expression levels greater than experimental measurement precision and ends with cell clustering based on a minimal set of Profiler genes, each one a putative biomarker of cluster-specific expression profiles. To benchmark how SALSA performs in experimental settings, we used the publicly available 10X Genomics PBMC 3K dataset, a pre-curated silver standard from human frozen peripheral blood comprising 2,700 single-cell barcodes, and identified 7 major cell groups matching transcriptional profiles of peripheral blood cell types and driven agnostically by < 500 Profiler genes. Finally, we demonstrate successful implementation of SALSA in a replicative scRNA-seq scenario by using previously published DropSeq data from a multi-batch mouse retina experimental design, thereby identifying 10 transcriptionally distinct cell types from > 64,000 single cells across 7 independent biological replicates based on < 630 Profiler genes. With these results, SALSA demonstrates that robust pattern detection from scRNA-seq expression matrices only requires a fraction of the accrued data, suggesting that single-cell sequencing technologies can become affordable and widespread if meant as hypothesis-generation tools to extract large-scale differential expression effects.

Specific expression and alternative splicing of mouse genes during spermatogenesis

Considering the high abundance of spliced RNAs in testis compared to other tissues, it is needed to construct the landscape of alternative splicing during spermatogenesis. However, there is still a lack of the systematic analysis of alternative RNA splicing in spermatogenesis. Here, we constructed a landscape of alternative RNA splicing during mouse spermatogenesis based on integrated RNA-seq data sets. Our results presented several novel alternatively spliced genes (Eif2s3y, Erdr1 Uty and Zfy1) in the Y chromosome with a specific expression pattern. Remarkably, the alternative splicing genes were grouped into co-expression networks involved in the microtubule cytoskeleton organization and post-transcriptional regulation of the gene expression, indicating the potential pathway to germ cell generation. Furthermore, based on the co-expression networks, we identified Atxn2l as a potential key gene in spermatogenesis, which presented dynamic expression patterns in different alternative splicing types. Ultimately, we proposed splicing regulatory networks for understanding novel and innovative alternative splicing regulation mechanisms during spermatogenesis. In summary, our research provides a systematic analysis of alternative RNA splicing and some novel spliced genes related to spermatogenesis.

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.

Thyroid hormone action in the developing testis: intergenerational epigenetics

Male fertility involves the successful transmission of the genetic code to the next generation. It requires appropriately timed cellular processes during testis development, adequate support of spermatogenesis by hormonal cues from the reproductive axis and cellular cross-talk between germ and somatic cells. In addition to being the vessel of the father’s genome, increasing evidence shows that the mature sperm carries valuable epigenetic information – the epigenome – that, after fecundation, influences the development of the next generation, affecting biological traits and disease susceptibility. The epigenome of the germ line is susceptible to environmental factors, including exogenous chemicals and diet, but it is also affected by endogenous molecules and pathophysiological conditions. Factors affecting testis development and the epigenetic information of the germ line are critical for fertility and of relevance to the non-genetic but heritable component in the etiology of complex conditions. Thyroid hormones are one of those factors and their action, when untimely, produces profound effects on the developing testis, affecting spermatogenesis, steroidogenesis, testis size, reproductive hormones and fertility. Altered thyroid hormone states can also change the epigenetic information of the male germ line, with phenotypic consequences for future generations. In the context of past literature concerning the consequences of altered thyroid hormone action for testis development, here we review recent findings about the pathophysiological roles of the principal determinants of testicular thyroid hormone action. We also discuss limited work on the effects of thyroid hormone on the male germ line epigenome and the implications for the intergenerational transmission of phenotypes via epigenetic mechanisms.

Spermatogonial Type 3 Deiodinase Regulates Thyroid Hormone Target Genes in Developing Testicular Somatic Cells

Premature overexposure to thyroid hormone causes profound effects on testis growth, spermatogenesis, and male fertility. We used genetic mouse models of type 3 deiodinase (DIO3) deficiency to determine the genetic programs affected by premature thyroid hormone action and to define the role of DIO3 in regulating thyroid hormone economy in testicular cells. Gene expression profiling in the neonatal testis of DIO3-deficient mice identified 5699 differentially expressed genes. Upregulated and downregulated genes were, respectively, involved according to DAVID analysis with cell differentiation and proliferation. They included anti-Müllerian hormone and genes involved in the formation of the blood-testis barrier, which are specific to Sertoli cells (SCs). They also included steroidogenic genes, which are specific to Leydig cells. Comparison with published data sets of genes enriched in SCs and spermatogonia, and responsive to retinoic acid (RA), identified a subset of genes that were regulated similarly by RA and thyroid hormone. This subset of genes showed an expression bias, as they were downregulated when enriched in spermatogonia and upregulated when enriched in SCs. Furthermore, using a genetic approach, we found that DIO3 is not expressed in SCs, but spermatogonia-specific inactivation of DIO3 led to impaired testis growth, reduced SC number, decreased cell proliferation and, especially during neonatal development, altered gene expression specific to somatic cells. These findings indicate that spermatogonial DIO3 protects testicular cells from untimely thyroid hormone signaling and demonstrate a mechanism of cross-talk between somatic and germ cells in the neonatal testis that involves the regulation of thyroid hormone availability and action.

Comprehensive Transcriptomic Analysis of Mouse Gonadal Development Involving Sexual Differentiation, Meiosis and Gametogenesis

Background

Mammalian gonadal development is crucial for fertility. Sexual differentiation, meiosis and gametogenesis are critical events in the process of gonadal development. Abnormalities in any of these events may cause infertility. However, owing to the complexity of these developmental events, the underlying molecular mechanisms are not fully understood and require further research.

Results

In this study, we employed RNA sequencing to examine transcriptome profiles of murine female and male gonads at crucial stages of these developmental events. By bioinformatics analysis, we identified a group of candidate genes that may participate in sexual differentiation, including Erbb3, Erbb4, and Prkg2. One hundred and two and 134 candidate genes that may be important for female and male gonadal development, respectively, were screened by analyzing the global gene expression patterns of developing female and male gonads. Weighted gene co-expression network analysis was performed on developing female gonads, and we identified a gene co-expression module related to meiosis. By alternative splicing analysis, we found that cassette-type exon and alternative start sites were the main forms of alternative splicing in developing gonads. A considerable portion of differentially expressed and alternatively spliced genes were involved in meiosis.

Conclusion

Taken together, our findings have enriched the gonadal transcriptome database and provided novel candidate genes and avenues to research the molecular mechanisms of sexual differentiation, meiosis, and gametogenesis.

Supplementary information

Supplementary information accompanies this paper at 10.1186/s12575-019-0108-y.

Development of fetal and adult Leydig cells

BACKGROUND

In mammals, two distinct Leydig cell populations, fetal Leydig cells (FLCs) and adult Leydig cells (ALCs), appear in the prenatal and postnatal testis, respectively. Although the functional differences between these cell types have been well described, the developmental relationship between FLCs and ALCs has not been fully understood. In this review, I focus on the cellular origins of FLCs and ALCs as well as the developmental and functional links between them.

METHODS

I surveyed previous reports about FLC and/or ALC development and summarized the findings.

MAIN FINDINGS

Fetal Leydig cells and ALCs were identified to have separate origins in the fetal and neonatal testis, respectively. However, several studies suggested that FLCs and ALCs share a common progenitor pool. Moreover, perturbation of FLC development at the fetal stage induces ALC dysfunction in adults, suggesting a functional link between FLCs and ALCs. Although the lineage relationship between FLCs and ALCs remains controversial, a recent study suggested that some FLCs dedifferentiate at the fetal stage, and that these cells serve as ALC stem cells.

CONCLUSION

Findings obtained from animal studies might provide clues to the causative mechanisms of male reproductive dysfunctions such as testicular dysgenesis syndrome in humans.

Nr5a1 suppression during the murine fetal period optimizes ovarian development by fine-tuning Notch signaling

The nuclear receptor NR5A1 is equally expressed and required for development of the gonadal primordia of both sexes, but, after sex determination, it is upregulated in XY testes and downregulated in XX ovaries. We have recently demonstrated, in mice, that this downregulation is mediated by forkhead box L2 (FOXL2) and hypothesized that adequate suppression of Nr5a1 is essential for normal ovarian development. Further, analysis of human patients with disorders/differences of sex development suggests that overexpression of NR5A1 can result in XX (ovo)testicular development. Here, we tested the role of Nr5a1 by overexpression in fetal gonads using a Wt1-BAC (bacterial artificial chromosome) transgene system. Enforced Nr5a1 expression compromised ovarian development in 46,XX mice, resulting in late-onset infertility, but did not induce (ovo)testis differentiation. The phenotype was similar to that of XX mice lacking Notch signaling. The expression level of Notch2 was significantly reduced in Nr5a1 transgenic mice, and the ovarian phenotype was almost completely rescued by in utero treatment with a NOTCH2 agonist. We conclude that suppression of Nr5a1 during the fetal period optimizes ovarian development by fine-tuning Notch signaling.

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.

At the Crossroads of Fate-Somatic Cell Lineage Specification in the Fetal Gonad

The reproductive endocrine systems are vastly different between males and females. This sexual dimorphism of the endocrine milieu originates from sex-specific differentiation of the somatic cells in the gonads during fetal life. Most gonadal somatic cells arise from the adrenogonadal primordium. After separation of the adrenal and gonadal primordia, the gonadal somatic cells initiate sex-specific differentiation during gonadal sex determination with the specification of the supporting cell lineages: Sertoli cells in the testis vs granulosa cells in the ovary. The supporting cell lineages then facilitate the differentiation of the steroidogenic cell lineages, Leydig cells in the testis and theca cells in the ovary. Proper differentiation of these cell types defines the somatic cell environment that is essential for germ cell development, hormone production, and establishment of the reproductive tracts. Impairment of lineage specification and function of gonadal somatic cells can lead to disorders of sexual development (DSDs) in humans. Human DSDs and processes for gonadal development have been successfully modeled using genetically modified mouse models. In this review, we focus on the fate decision processes from the initial stage of formation of the adrenogonadal primordium in the embryo to the maintenance of the somatic cell identities in the gonads when they become fully differentiated in adulthood.

Recent progress in understanding the mechanisms of Leydig cell differentiation

Leydig cells in fetal and adult testes play pivotal roles in eliciting male characteristics by producing androgen. Although numerous studies of Leydig cells have been performed, the mechanisms for differentiation of the two cell types (fetal Leydig and adult Leydig cells), their developmental and functional relationship, and their differential characteristics remain largely unclear. Based on recent technical progress in genome-wide analysis and in vitro investigation, novel and fascinating observations concerning the issues above have been obtained. Focusing on fetal and adult Leydig cells, this review summarizes the recent progress that has advanced our understanding of the cells.

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.

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.

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.

Development, function and fate of fetal Leydig cells

During fetal testis development, fetal Leydig cells (FLCs) are found to be originated from multiple progenitor cells. FLC specification and function are under tight regulation of specific genes and signaling proteins. Furthermore, Sertoli cells play a crucial role to regulate FLC differentiation during fetal testis development. FLC progenitor- and FLC-produced biomolecules are also involved in the differentiation and activity of rodent FLCs. The main function of FLCs is to produce androgens to masculinize XY embryos. However, FLCs are capable of producing androstenedione but not testosterone due to the lack of 17β-HSD (17β-hydroxysteroid dehydrogenase), but fetal Sertoli cells express 17β-HSD which thus transforms androstenedione to testosterone in the fetal testis. On the other hand, FLCs produce activin A to regulate Sertoli cell proliferation, and Sertoli cells in turn modulate testis cord expansion. It is now generally accepted that adult Leydig cells (ALCs) gradually replace FLCs during postnatal development to produce testosterone to support spermatogenesis as FLCs undergo degeneration in neonatal and pre-pubertal testes. However, based on studies using genetic tracing mouse models, FLCs are found to persist in adult testes, making up ∼20% of total Leydig cells. In this review, we evaluate the latest findings regarding the development, function and fate of FLCs during fetal and adult testis development.

Mapping lineage progression of somatic progenitor cells in the mouse fetal testis

Testis morphogenesis is a highly orchestrated process involving lineage determination of male germ cells and somatic cell types. Although the origin and differentiation of germ cells are known, the developmental course specific for each somatic cell lineage has not been clearly defined. Here, we construct a comprehensive map of somatic cell lineage progression in the mouse testis. Both supporting and interstitial cell lineages arise from WT1⁺ somatic progenitor pools in the gonadal primordium. A subpopulation of WT1⁺ progenitor cells acquire SOX9 expression and become Sertoli cells that form testis cords, whereas the remaining WT1⁺ cells contribute to progenitor cells in the testis interstitium. Interstitial progenitor cells diversify through the acquisition of HES1, an indication of Notch activation, at the onset of sex determination. HES1⁺ interstitial progenitors, through the action of Sertoli cell-derived Hedgehog signals, become positive for GLI1. The GLI1⁺ interstitial cells eventually develop into two cell lineages: steroid-producing fetal Leydig cells and non-steroidogenic cells. The fetal Leydig cell population is restricted by Notch2 signaling from the neighboring somatic cells. The non-steroidogenic progenitor cells retain their undifferentiated state during fetal stage and become adult Leydig cells in post-pubertal testis. These results provide the first lineage progression map that illustrates the sequential establishment of somatic cell populations during testis morphogenesis.

Isolation and Characterization of Fetal Leydig Progenitor Cells of Male Mice

Fetal and adult Leydig cells develop in mammalian prenatal and postnatal testes, respectively. In mice, fetal Leydig cells (FLCs) emerge in the interstitial space of the testis at embryonic day 12.5 and thereafter increase in number, possibly through differentiation from progenitor cells. However, the progenitor cells have not yet been identified. Previously, we established transgenic mice in which FLCs are labeled strongly with enhanced green fluorescent protein (EGFP). Interestingly, fluorescence-activated cell sorting provided us with weakly EGFP-labeled cells as well as strongly EGFP-labeled FLCs. In vitro reconstruction of fetal testes demonstrated that weakly EGFP-labeled cells contain FLC progenitors. Transcriptome from the 2 cell populations revealed, as expected, marked differences in the expression of genes required for growth factor/receptor signaling and steroidogenesis. In addition, genes for energy metabolisms such as glycolytic pathways and the citrate cycle were activated in strongly EGFP-labeled cells, suggesting that metabolism is activated during FLC differentiation.