Genomics of sexual cell fate transdifferentiation in the mouse gonad

Diminished DNA binding affinity of DMRT1 caused by heterozygous DM domain mutations is a cause of male infertility

The most severe form of male infertility is idiopathic non-obstructive azoospermia (NOA), a complete sperm absence in the ejaculate. We performed exome sequencing in the Croatian infertile brothers with NOA and found a variant in DMRT1 (Doublesex and mab-3 related transcription factor 1) gene that was further assessed by the EMSA assay and molecular dynamic simulations. We additionally screened for DMRT1 mutations in 1940 infertile men diagnosed with spermatogenic failure, 644 normozoospermic controls, and 105 females with primary ovarian insufficiency (POI) recruited to the GEnetics of Male INfertility Initiative (GEMINI) or Estonian Andrology (ESTAND) cohorts. DMRT1 p.Pro74Leu (chr9:g.842059C > T) variant was detected in infertile brothers in the highly conserved position within the DNA binding DM domain of the protein. EMSA assay showed reduced DNA binding of DMRT1P74L and molecular dynamic simulations showed differences in structural and dynamical properties between the wild type protein and DMRT1P74L. Plausible disease-causing DMRT1 variants were only identified in infertile men (13/1940; 0.67%), and none in 639 fertile controls. Burden testing showed an excess of rare deleterious DM domain mutations in the infertility cohort compared to gnomAD v.4.0 population-based controls (Fisher's exact test, p = 1.44 x 10-5). Three rare deleterious variants in DMRT1 were found in 104 cases of POI. The findings of this study strengthen the evidence of DMRT1 variants being a causal factor for male infertility and provide the distribution of likely pathogenic variants across the gene. This is also the first study to suggest that DMRT1 variants may also be linked to POI.

Postnatal Ovarian Transdifferentiation in the Absence of Estrogen Receptor Signaling Is Dependent on Genetic Background

Normal ovarian function requires the expression of estrogen receptors α (ESR1) and β (ESR2) in distinct cell types within the ovary. The double estrogen receptor knockout (αβERKO) ovary had the appearance of seminiferous tubule-like structures that expressed SOX9; this phenotype was lost when the animals were repeatedly backcrossed to the C57BL/6J genetic background. A new line of ERKO mice, Ex3αβERKO, was developed for targeted disruption on a mixed genetic background. Histological examination of the ovaries in the Ex3αβERKO showed the appearance of seminiferous tubule-like structures in mice aged 6 to 12 months. These dismorphogenic regions have cells that no longer express granulosa cell-specific FOXL2, while other cells express Sertoli cell-specific SOX9 as examined by immunohistochemistry. Whole ovarian gene expression analysis in Ex3αERKO, Ex3βRKO, and Ex3αβERKO found many genes differentially expressed compared to controls with one Esr1 and Esr2 allele. The genes specific to the Ex3αβERKO ovary were compared to other models of postnatal ovarian transdifferentiation, identifying 21 candidate genes. To examine the genetic background contributions, DNA was isolated from αβERKO mice that did not show ovarian transdifferentiation and compared to DNA from Ex3αβERKO using Mouse Diversity Array. A genomic region putatively associated with transdifferentiation was identified on Chr18 (5-15 M) and genes in this region were compared to the genes differentially expressed in models of ovarian transdifferentiation. This work demonstrates the importance of ESRs in maintaining granulosa cell differentiation within the ovary, identifies several potential gene candidates, and suggests that genetic background can be a confounding factor.

Loss of NR5A1 in mouse Sertoli cells after sex determination changes cellular identity and induces cell death by anoikis

To investigate the role of the nuclear receptor NR5A1 in the testis after sex determination, we analyzed mice lacking NR5A1 in Sertoli cells (SCs) from embryonic day (E) 13.5 onwards. Ablation of Nr5a1 impaired the expression of genes characteristic of SC identity (e.g. Sox9 and Amh), caused SC death from E14.5 onwards through a Trp53-independent mechanism related to anoikis, and induced disorganization of the testis cords. Together, these effects caused germ cells to enter meiosis and die. Single-cell RNA-sequencing experiments revealed that NR5A1-deficient SCs changed their molecular identity: some acquired a 'pre-granulosa-like' cell identity, whereas other reverted to a 'supporting progenitor-like' cell identity, most of them being 'intersex' because they expressed both testicular and ovarian genes. Fetal Leydig cells (LCs) did not display significant changes, indicating that SCs are not required beyond E14.5 for their emergence or maintenance. In contrast, adult LCs were absent from postnatal testes. In addition, adult mutant males displayed persistence of Müllerian duct derivatives, decreased anogenital distance and reduced penis length, which could be explained by the loss of AMH and testosterone synthesis due to SC failure.

A comparative analysis of genes differentially expressed between rete testis cells and Sertoli cells of the mouse testis

The rete testis (RT) is a region of the mammalian testis that plays an important role in testicular physiology. The RT epithelium consists of cells sharing some well-known gene markers with supporting Sertoli cells (SCs). However, little is known about the differences in gene expression between these two cell populations. Here, we used fluorescence-activated cell sorting (FACS) to obtain pure cultures of neonatal RT cells and SCs and identified differentially expressed genes (DEGs) between these cell types using RNA sequencing (RNA-seq). We then compared our data with the RNA-seq data of other studies that examined RT cells and SCs of mice of different ages and generated a list of DEGs permanently upregulated in RT cells throughout testis development and in culture, which included 86 genes, and a list of 79 DEGs permanently upregulated in SCs. The analysis of studies on DMRT1 function revealed that nearly half of the permanent DEGs could be regulated by this SC upregulated transcription factor. We suggest that useful cell lineage markers and candidate genes for the specification of both RT cells and SCs may be present among these permanent DEGs.