The p.R92W variant of NR5A1/Nr5a1 induces testicular development of 46,XX gonads in humans, but not in mice: phenotypic comparison of human patients and mutation-induced mice

Towards a "Testis in a Dish": Generation of Mouse Testicular Organoids that Recapitulate Testis Structure and Expression Profiles

The testis is responsible for sperm production and androgen synthesis. Abnormalities in testis development and function lead to disorders of sex development and male infertility. Currently, no in vitro system exists for modelling the testis. Here, we generated testis organoids from neonatal mouse primary testicular cells using transwell inserts and show that these organoids generate tubule-like structures and cellular organization resembling that of the in vivo testis. Gene expression analysis of organoids demonstrates a profile that recapitulates that observed in in vivo testis. Embryonic testicular cells, but not adult testicular cells are also capable of forming organoids. These organoids can be maintained in culture for 8-9 weeks and shows signs of entry into meiosis. We further developed defined media compositions that promote the immature versus mature Sertoli cell and Leydig cell states, enabling organoid maturation in vitro. These testis organoids are a promising model system for basic research of testes development and function, with translational applications for elucidation and treatment of developmental sex disorders and infertility.

Nuclear Receptor Gene Variants Underlying Disorders/Differences of Sex Development through Abnormal Testicular Development

Gonadal development is the first step in human reproduction. Aberrant gonadal development during the fetal period is a major cause of disorders/differences of sex development (DSD). To date, pathogenic variants of three nuclear receptor genes (NR5A1, NR0B1, and NR2F2) have been reported to cause DSD via atypical testicular development. In this review article, we describe the clinical significance of the NR5A1 variants as the cause of DSD and introduce novel findings from recent studies. NR5A1 variants are associated with 46,XY DSD and 46,XX testicular/ovotesticular DSD. Notably, both 46,XX DSD and 46,XY DSD caused by the NR5A1 variants show remarkable phenotypic variability, to which digenic/oligogenic inheritances potentially contribute. Additionally, we discuss the roles of NR0B1 and NR2F2 in the etiology of DSD. NR0B1 acts as an anti-testicular gene. Duplications containing NR0B1 result in 46,XY DSD, whereas deletions encompassing NR0B1 can underlie 46,XX testicular/ovotesticular DSD. NR2F2 has recently been reported as a causative gene for 46,XX testicular/ovotesticular DSD and possibly for 46,XY DSD, although the role of NR2F2 in gonadal development is unclear. The knowledge about these three nuclear receptors provides novel insights into the molecular networks involved in the gonadal development in human fetuses.

In vitro cellular reprogramming to model gonad development and its disorders

During embryonic development, mutually antagonistic signaling cascades determine gonadal fate toward a testicular or ovarian identity. Errors in this process result in disorders of sex development (DSDs), characterized by discordance between chromosomal, gonadal, and anatomical sex. The absence of an appropriate, accessible in vitro system is a major obstacle in understanding mechanisms of sex-determination/DSDs. Here, we describe protocols for differentiation of mouse and human pluripotent cells toward gonadal progenitors. Transcriptomic analysis reveals that the in vitro-derived murine gonadal cells are equivalent to embryonic day 11.5 in vivo progenitors. Using similar conditions, Sertoli-like cells derived from 46,XY human induced pluripotent stem cells (hiPSCs) exhibit sustained expression of testis-specific genes, secrete anti-Müllerian hormone, migrate, and form tubular structures. Cells derived from 46,XY DSD female hiPSCs, carrying an NR5A1 variant, show aberrant gene expression and absence of tubule formation. CRISPR-Cas9-mediated variant correction rescued the phenotype. This is a robust tool to understand mechanisms of sex determination and model DSDs.

Becoming female: Ovarian differentiation from an evolutionary perspective

Differentiation of the bipotential gonadal primordium into ovaries and testes is a common process among vertebrate species. While vertebrate ovaries eventually share the same functions of producing oocytes and estrogens, ovarian differentiation relies on different morphogenetic, cellular, and molecular cues depending on species. The aim of this review is to highlight the conserved and divergent features of ovarian differentiation through an evolutionary perspective. From teleosts to mammals, each clade or species has a different story to tell. For this purpose, this review focuses on three specific aspects of ovarian differentiation: ovarian morphogenesis, the evolution of the role of estrogens on ovarian differentiation and the molecular pathways involved in granulosa cell determination and maintenance.

A fish with no sex: gonadal and adrenal functions partition between zebrafish NR5A1 co-orthologs

People with NR5A1 mutations experience testicular dysgenesis, ovotestes, or adrenal insufficiency, but we do not completely understand the origin of this phenotypic diversity. NR5A1 is expressed in gonadal soma precursor cells before expression of the sex-determining gene SRY. Many fish have two co-orthologs of NR5A1 that likely partitioned ancestral gene subfunctions between them. To explore ancestral roles of NR5A1, we knocked out nr5a1a and nr5a1b in zebrafish. Single-cell RNA-seq identified nr5a1a-expressing cells that co-expressed genes for steroid biosynthesis and the chemokine receptor Cxcl12a in 1-day postfertilization (dpf) embryos, as does the mammalian adrenal-gonadal (interrenal-gonadal) primordium. In 2dpf embryos, nr5a1a was expressed stronger in the interrenal-gonadal primordium than in the early hypothalamus but nr5a1b showed the reverse. Adult Leydig cells expressed both ohnologs and granulosa cells expressed nr5a1a stronger than nr5a1b. Mutants for nr5a1a lacked the interrenal, formed incompletely differentiated testes, had no Leydig cells, and grew far larger than normal fish. Mutants for nr5a1b formed a disorganized interrenal and their gonads completely disappeared. All homozygous mutant genotypes lacked secondary sex characteristics, including male breeding tubercles and female sex papillae, and had exceedingly low levels of estradiol, 11-ketotestosterone, and cortisol. RNA-seq showed that at 21dpf, some animals were developing as females and others were not, independent of nr5a1 genotype. By 35dpf, all mutant genotypes greatly under-expressed ovary-biased genes. Because adult nr5a1a mutants form gonads but lack an interrenal and conversely, adult nr5a1b mutants lack a gonad but have an interrenal, the adrenal, and gonadal functions of the ancestral nr5a1 gene partitioned between ohnologs after the teleost genome duplication, likely owing to reciprocal loss of ancestral tissue-specific regulatory elements. Identifying such elements could provide hints to otherwise unexplained cases of Differences in Sex Development.

Identification of the first promoter-specific gain-of-function SOX9 missense variant (p.E50K) in a patient with 46,XX ovotesticular disorder of sex development

SOX9, a transcription factor, is expressed in the undifferentiated XX and XY gonads. SRY induces significant upregulation of SOX9 expression in XY gonads. Loss-of-function SOX9 variants cause testicular dysgenesis in 46,XY patients, while duplication of the total gene or the upstream regulatory region results in testicular development in 46,XX patients. However, gain-of-function (GoF) SOX9 variants have not been reported previously. We report the case of a 16-year-old female patient with a 46,XX karyotype who had masculinized external genitalia and unilateral ovotestis. Next-generation sequencing-based genetic screening for disorders of sex development led to the identification of a novel SOX9 variant (p.Glu50Lys), transmitted from the phenotypically normal father. Expression analysis showed that E50K-SOX9 enhanced transactivation of the luciferase reporter containing the testis enhancer sequence core element compared with that containing the wildtype-SOX9. This GoF activity was not observed in the luciferase reporter containing Amh, the gene for anti-Müllerian hormone. We genetically engineered female mice (Sox9^E50K/E50K ), and they showed no abnormalities in the external genitalia or ovaries. In conclusion, a novel SOX9 variant with a promoter-specific GoF activity was identified in vitro; however, the disease phenotype was not recapitulated by the mouse model. At present, the association between the GoF SOX9 variant and the ovotestis phenotype remains unclear. Future studies are needed to verify the possible association.

Male Hypogonadism and Disorders of Sex Development

Disorders of Sex Development (DSD) are congenital anomalies in which there is a discordance between chromosomal, genetic, gonadal, and/or internal/external genital sex. In XY individuals, the process of fetal sex differentiation can be disrupted at the stage of gonadal differentiation, resulting in gonadal dysgenesis, a form of early fetal-onset primary hypogonadism characterized by insufficient androgen and anti-Müllerian hormone (AMH) production, which leads to the development of ambiguous or female genitalia. The process of sex differentiation can also be disrupted at the stage of genital differentiation, due to isolated defects in androgen or AMH secretion, but not both. These are forms of fetal-onset hypogonadism with dissociated gonadal dysfunction. In this review, we present a perspective on impaired testicular endocrine function, i.e., fetal-onset male hypogonadism, resulting in incomplete virilization at birth.

Molecular Characterization of XX Maleness

Androgens and anti-Müllerian hormone (AMH), secreted by the foetal testis, are responsible for the development of male reproductive organs and the regression of female anlagen. Virilization of the reproductive tract in association with the absence of Müllerian derivatives in the XX foetus implies the existence of testicular tissue, which can occur in the presence or absence of SRY. Recent advancement in the knowledge of the opposing gene cascades driving to the differentiation of the gonadal ridge into testes or ovaries during early foetal development has provided insight into the molecular explanation of XX maleness.

A mutation inactivating the distal SF1 binding site on the human anti-Müllerian hormone promoter causes persistent Müllerian duct syndrome

The persistent Müllerian duct syndrome (PMDS) is a 46,XY disorder of sexual development characterized by the persistence of Müllerian duct derivatives, uterus and tubes, in otherwise normally masculinized males. The condition, transmitted as a recessive autosomal trait, is usually due to mutations in either the anti-Müllerian hormone (AMH) gene or its main receptor. Many variants of these genes have been described, all targeting the coding sequences. We report the first case of PMDS due to a regulatory mutation. The AMH promoter contains two binding sites for steroidogenic factor 1 (SF1), one at −102 and the other at −228. Our patient carries a single base deletion at −225, significantly decreasing its capacity for binding SF1, as measured by the electrophoresis mobility shift assay. Furthermore, by linking the AMH promoter to the luciferase gene, we show that the transactivation capacity of the promoter is significantly decreased by the mutation, in contrast to the disruption of the −102 binding site. To explain the difference in impact we hypothesize that SF1 could partially overcome the lack of binding to the −102 binding site by interacting with a GATA4 molecule linked to a nearby response element. We show that disruption of both the −102 SF1 and the −84 GATA response elements significantly decreases the transactivation capacity of the promoter. In conclusion, we suggest that the distance between mutated SF1 sites and potentially rescuing GATA binding motifs might play a role in the development of PMDS.

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.

Loss of Function of the Nuclear Receptor NR2F2, Encoding COUP-TF2, Causes Testis Development and Cardiac Defects in 46,XX Children

Emerging evidence from murine studies suggests that mammalian sex determination is the outcome of an imbalance between mutually antagonistic male and female regulatory networks that canalize development down one pathway while actively repressing the other. However, in contrast to testis formation, the gene regulatory pathways governing mammalian ovary development have remained elusive. We performed exome or Sanger sequencing on 79 46,XX SRY-negative individuals with either unexplained virilization or with testicular/ovotesticular disorders/differences of sex development (TDSD/OTDSD). We identified heterozygous frameshift mutations in NR2F2, encoding COUP-TF2, in three children. One carried a c.103_109delGGCGCCC (p.Gly35Argfs^∗75) mutation, while two others carried a c.97_103delCCGCCCG (p.Pro33Alafs^∗77) mutation. In two of three children the mutation was de novo. All three children presented with congenital heart disease (CHD), one child with congenital diaphragmatic hernia (CDH), and two children with blepharophimosis-ptosis-epicanthus inversus syndrome (BPES). The three children had androgen production, virilization of external genitalia, and biochemical or histological evidence of testicular tissue. We demonstrate a highly significant association between the NR2F2 loss-of-function mutations and this syndromic form of DSD (p = 2.44 × 10⁻⁸). We show that COUP-TF2 is highly abundant in a FOXL2-negative stromal cell population of the fetal human ovary. In contrast to the mouse, these data establish COUP-TF2 as a human “pro-ovary” and “anti-testis” sex-determining factor in female gonads. Furthermore, the data presented here provide additional evidence of the emerging importance of nuclear receptors in establishing human ovarian identity and indicate that nuclear receptors may have divergent functions in mouse and human biology.

Phenotypic Variation in 46,XX Disorders of Sex Development due to the NR5A1 p.R92W Variant: A Sibling Case Report and Literature Review

Recently, a heterozygous missense mutation in NR5A1, p.R92W, was identified as a cause of 46,XX testicular/ovo-testicular disorders of sexual development (DSD). We report a sibling pair with 46,XX DSD due to an NR5A1 mutation with distinct phenotypes, including external and internal genitalia and gonads, for whom different rearing sexes were selected. Thus, the phenotypes of p.R92W vary, even within a family. The father of the patients showed oligozoospermia with the p.R92W mutation, suggesting that in 46,XY individuals, the mutation would cause various gonadal phenotypes. We review and discuss the general role of the R92W mutation in sexual development.

Next generation sequencing and array-based comparative genomic hybridization for molecular diagnosis of pediatric endocrine disorders

Next-generation sequencing (NGS) and array-based comparative genomic hybridization (array CGH) have enabled us to perform high-throughput mutation screening and genome-wide copy number analysis, respectively. These methods can be used for molecular diagnosis of pediatric endocrine disorders. NGS has determined the frequency and phenotypic variation of mutations in several disease-associated genes. Furthermore, whole exome analysis using NGS has successfully identified several novel causative genes for endocrine disorders. Array CGH is currently used as the standard procedure for molecular cytogenetic analysis. Array CGH can detect various submicroscopic genomic rearrangements involving exons or enhancers of disease-associated genes. This review introduces some examples of the use of NGS and array CGH for the molecular diagnosis of pediatric endocrine disorders.

New NR5A1 mutations and phenotypic variations of gonadal dysgenesis

Mutations in NR5A1 have been reported as a frequent cause of 46,XY disorders of sex development (DSD) associated to a broad phenotypic spectrum ranging from infertility, ambiguous genitalia, anorchia to gonadal dygenesis and female genitalia. Here we present the clinical follow up of four 46,XY DSD patients with three novel heterozygous mutations in the NR5A1 gene leading to a p.T40P missense mutation and a p.¹⁸DKVSG²²del nonframeshift deletion in the DNA-binding domain and a familiar p.Y211Tfs*83 frameshift mutation. Functional analysis of the missense and nonframeshift mutation revealed a deleterious character with loss of DNA-binding and transactivation capacity. Both, the mutations in the DNA-binding domain, as well as the familiar frameshift mutation are associated with highly variable endocrine values and phenotypic appearance. Phenotypes vary from males with spontaneous puberty, substantial testosterone production and possible fertility to females with and without Müllerian structures and primary amenorrhea. Exome sequencing of the sibling’s family revealed TBX2 as a possible modifier of gonadal development in patients with NR5A1 mutations.

CRISPR/Cas9-mediated simultaneous knockout of Dmrt1 and Dmrt3 does not recapitulate the 46,XY gonadal dysgenesis observed in 9p24.3 deletion patients

DM domain transcription factors play important roles in sexual development in a wide variety of species from invertebrate to humans. Among seven mammalian family members of DM domain transcription factors, DMRT1 has been studied in mouse and human for its conserved role in male gonadal identity. Chromosomal deletion of 9p24.3, the region in which DMRT1 is located, is associated with 46,XY gonadal dysgenesis. Dmrt1 knockout (KO) mice also showed male-to-female gonadal reprogramming. However, the phenotype of Dmrt1 KO mouse appears only after birth while 46,XY gonadal dysgenesis occurs during the developmental phase, and the cause behind this difference remained unknown. We hypothesized that in human the function of other DMRT genes clustered with DMRT1, namely DMRT3, might also be impaired by the chromosomal deletion, which leads to the gonadal dysgenesis phenotype. Thus, simultaneous loss of multiple DM domain genes in mice could have a more severe impact on gonadal development. To address this issue, we generated double KO mice for Dmrt1 and Dmrt3 via the CRISPR/Cas9 system. Comparing adult and neonatal testes of single and double KO mice, we found that loss of Dmrt1 or Dmrt3, or both, does not have apparent effect on male gonadal formation during embryonic development. Our study demonstrated that the discrepancy between human with 9p24.3 deletion and Dmrt1 KO mouse could not be explained by the simultaneous loss of Dmrt3 gene. CRISPR/Cas9 is a versatile and straightforward approach to elucidate the questions that were otherwise difficult to address with conventional methods.

•

We have generated Dmrt1 KO, Dmrt3 KO and Dmrt1/3 DKO mice via the CRISPR/Cas9 system.

•

.CRISPR/Cas9-mediated gene KO could generate 3 KO lines via single microinjection.

•

Double KO of Dmrt1 and Dmrt3 could not recapitulate human 9p24.3 deletion symptom.