Novel Genomic Variants, Atypical Phenotypes and Evidence of a Digenic/Oligogenic Contribution to Disorders/Differences of Sex Development in a Large North African Cohort

NR5A1/SF-1 Collaborates with Inhibin α and the Androgen Receptor

Steroidogenic factor 1 (SF-1) is a nuclear receptor that regulates steroidogenesis and reproductive development. NR5A1/SF-1 variants are associated with a broad spectrum of phenotypes across individuals with disorders of sex development (DSDs). Oligogenic inheritance has been suggested as an explanation. SF-1 interacts with numerous partners. Here, we investigated a constellation of gene variants identified in a 46,XY severely undervirilized individual carrying an ACMG-categorized 'pathogenic' NR5A1/SF-1 variant in comparison to the healthy carrier father. Candidate genes were revealed by whole exome sequencing, and pathogenicity was predicted by different in silico tools. We found variants in NR1H2 and INHA associated with steroidogenesis, sex development, and reproduction. The identified variants were tested in cell models. Novel SF-1 and NR1H2 binding sites in the AR and INHA gene promoters were found. Transactivation studies showed that wild-type NR5A1/SF-1 regulates INHA and AR gene expression, while the NR5A1/SF-1 variant had decreased transcriptional activity. NR1H2 was found to regulate AR gene transcription; however, the NR1H2 variant showed normal activity. This study expands the NR5A1/SF-1 network of interacting partners, while not solving the exact interplay of different variants that might be involved in revealing the observed DSD phenotype. It also illustrates that understanding complex genetics in DSDs is challenging.

Deleterious ZNRF3 germline variants cause neurodevelopmental disorders with mirror brain phenotypes via domain-specific effects on Wnt/β-catenin signaling

Zinc and RING finger 3 (ZNRF3) is a negative-feedback regulator of Wnt/β-catenin signaling, which plays an important role in human brain development. Although somatically frequently mutated in cancer, germline variants in ZNRF3 have not been established as causative for neurodevelopmental disorders (NDDs). We identified 12 individuals with ZNRF3 variants and various phenotypes via GeneMatcher/Decipher and evaluated genotype-phenotype correlation. We performed structural modeling and representative deleterious and control variants were assessed using in vitro transcriptional reporter assays with and without Wnt-ligand Wnt3a and/or Wnt-potentiator R-spondin (RSPO). Eight individuals harbored de novo missense variants and presented with NDD. We found missense variants associated with macrocephalic NDD to cluster in the RING ligase domain. Structural modeling predicted disruption of the ubiquitin ligase function likely compromising Wnt receptor turnover. Accordingly, the functional assays showed enhanced Wnt/β-catenin signaling for these variants in a dominant negative manner. Contrarily, an individual with microcephalic NDD harbored a missense variant in the RSPO-binding domain predicted to disrupt binding affinity to RSPO and showed attenuated Wnt/β-catenin signaling in the same assays. Additionally, four individuals harbored de novo truncating or de novo or inherited large in-frame deletion variants with non-NDD phenotypes, including heart, adrenal, or nephrotic problems. In contrast to NDD-associated missense variants, the effects on Wnt/β-catenin signaling were comparable between the truncating variant and the empty vector and between benign variants and the wild type. In summary, we provide evidence for mirror brain size phenotypes caused by distinct pathomechanisms in Wnt/β-catenin signaling through protein domain-specific deleterious ZNRF3 germline missense variants.

Evidence for NR2F2/COUP-TFII involvement in human testis development

NR2F2 encodes COUP-TFII, an orphan nuclear receptor required for the development of the steroidogenic lineages of the murine fetal testes and ovaries. Pathogenic variants in human NR2F2 are associated with testis formation in 46,XX individuals, however, the function of COUP-TFII in the human testis is unknown. We report a de novo heterozygous variant in NR2F2 (c.737G > A, p.Arg246His) in a 46,XY under-masculinized boy with primary hypogonadism. The variant, located within the ligand-binding domain, is predicted to be highly damaging. In vitro studies indicated that the mutation does not impact the stability or subcellular localization of the protein. NR5A1, a related nuclear receptor that is a key factor in gonad formation and function, is known to physically interact with COUP-TFII to regulate gene expression. The mutant protein did not affect the physical interaction with NR5A1. However, in-vitro assays demonstrated that the mutant protein significantly loses the inhibitory effect on NR5A1-mediated activation of both the LHB and INSL3 promoters. The data support a role for COUP-TFII in human testis formation. Although mutually antagonistic sets of genes are known to regulate testis and ovarian pathways, we extend the list of genes, that together with NR5A1 and WT1, are associated with both 46,XX and 46,XY DSD.

Mini-Puberty, Physiological and Disordered: Consequences, and Potential for Therapeutic Replacement

There are 3 physiological waves of central hypothalamic-pituitary-gonadal (HPG) axis activity over the lifetime. The first occurs during fetal life, the second-termed "mini-puberty"-in the first months after birth, and the third at puberty. After adolescence, the axis remains active all through adulthood. Congenital hypogonadotropic hypogonadism (CHH) is a rare genetic disorder characterized by a deficiency in hypothalamic gonadotropin-releasing hormone (GnRH) secretion or action. In cases of severe CHH, all 3 waves of GnRH pulsatility are absent. The absence of fetal HPG axis activation manifests in around 50% of male newborns with micropenis and/or undescended testes (cryptorchidism). In these boys, the lack of the mini-puberty phase accentuates testicular immaturity. This is characterized by a low number of Sertoli cells, which are important for future reproductive capacity. Thus, absent mini-puberty will have detrimental effects on later fertility in these males. The diagnosis of CHH is often missed in infants, and even if recognized, there is no consensus on optimal therapeutic management. Here we review physiological mini-puberty and consequences of central HPG axis disorders; provide a diagnostic approach to allow for early identification of these conditions; and review current treatment options for replacement of mini-puberty in male infants with CHH. There is evidence from small case series that replacement with gonadotropins to mimic "mini-puberty" in males could have beneficial outcomes not only regarding testis descent, but also normalization of testis and penile sizes. Moreover, such therapeutic replacement regimens in disordered mini-puberty could address both reproductive and nonreproductive implications.

Worldwide cohort study of 46, XY differences/disorders of sex development genetic diagnoses: geographic and ethnic differences in variants

Differences/disorders of sex development (DSDs) in individuals with a 46, XY karyotype are a group of congenital disorders that manifest as male gonadal hypoplasia or abnormalities of the external genitalia. Approximately 50% of patients with 46, XY DSDs cannot obtain a molecular diagnosis. The aims of this paper were to review the most common causative genes and rare genes in patients with 46, XY DSDs, analyze global molecular diagnostic cohorts for the prevalence and geographic distribution of causative genes, and identify the factors affecting cohort detection results. Although the spectrum of genetic variants varies across regions and the severity of the clinical phenotype varies across patients, next-generation sequencing (NGS), the most commonly used detection method, can still reveal genetic variants and aid in diagnosis. A comparison of the detection rates of various sequencing modalities revealed that whole-exome sequencing (WES) facilitates a greater rate of molecular diagnosis of the disease than panel sequencing. Whole-genome sequencing (WGS), third-generation sequencing, and algorithm advancements will contribute to the improvement of detection efficiency. The most commonly mutated genes associated with androgen synthesis and action are AR, SR5A2, and HSD17B3, and the most commonly mutated genes involved in gonadal formation are NR5A1 and MAP3K1. Detection results are affected by differences in enrollment criteria and sequencing technologies.

PTPN11 and FLNA variants in a boy with ambiguous genitalia, short stature, and non-specific dysmorphic features

Noonan syndrome is a congenital disorder characterized by distinctive facial appearance, congenital heart defects, short stature, and skeletal dysplasia. Although boys with Noonan syndrome frequently exhibit cryptorchidism, a mild form of 46,XY disorders of sex development (DSD), they barely manifest more severe genital abnormalities. Here, we report a boy with ambiguous genitalia, short stature, and non-specific dysmorphic features. He had no cardiac abnormalities or skeletal dysplasia. His score in the Noonan syndrome diagnostic criteria (36 of 157 points, 23%) was lower than the cutoff for diagnosis (50%). Whole-exome sequencing identified a de novo heterozygous variant (c.922A>G: p.Asn308Asp) in PTPN11 and a maternally inherited hemizygous variant (c.1439C>T: p.Pro480Leu) in FLNA. The PTPN11 variant was a known causative mutation for Noonan syndrome. FLNA is a causative gene for neurodevelopmental and skeletal abnormalities and has also been implicated in 46,XY DSD. The p.Pro480Leu variant of FLNA was assessed as deleterious by in silico analyses. These results provide evidence that whole-exome sequencing is a powerful tool for diagnosing patients with atypical disease manifestations. Furthermore, our data suggest a possible role of digenic mutations as phenotypic modifiers of Noonan syndrome.

Clinical and genetic characteristics of a large international cohort of individuals with rare NR5A1/SF-1 variants of sex development

BACKGROUND

Steroidogenic factor 1 (SF-1/NR5A1) is essential for human sex development. Heterozygous NR5A1/SF-1 variants manifest with a broad range of phenotypes of differences of sex development (DSD), which remain unexplained.

METHODS

We conducted a retrospective analysis on the so far largest international cohort of individuals with NR5A1/SF-1 variants, identified through the I-DSD registry and a research network.

FINDINGS

Among 197 individuals with NR5A1/SF-1 variants, we confirmed diverse phenotypes. Over 70% of 46, XY individuals had a severe DSD phenotype, while 90% of 46, XX individuals had female-typical sex development. Close to 100 different novel and known NR5A1/SF-1 variants were identified, without specific hot spots. Additionally, likely disease-associated variants in other genes were reported in 32 individuals out of 128 tested (25%), particularly in those with severe or opposite sex DSD phenotypes. Interestingly, 48% of these variants were found in known DSD or SF-1 interacting genes, but no frequent gene-clusters were identified. Sex registration at birth varied, with <10% undergoing reassignment. Gonadectomy was performed in 30% and genital surgery in 58%. Associated organ anomalies were observed in 27% of individuals with a DSD, mainly concerning the spleen. Intrafamilial phenotypes also varied considerably.

INTERPRETATION

The observed phenotypic variability in individuals and families with NR5A1/SF-1 variants is large and remains unpredictable. It may often not be solely explained by the monogenic pathogenicity of the NR5A1/SF-1 variants but is likely influenced by additional genetic variants and as-yet-unknown factors.

FUNDING

Swiss National Science Foundation (320030-197725) and Boveri Foundation Zürich, Switzerland.

Clinical characteristics and genetic expansion of 46,XY disorders of sex development children in a Chinese prospective study

Diagnosis and management strategy of disorders of sex development (DSD) are difficult and various due to heterogeneous phenotype and genotype. Under widespread use of genomic sequencing technologies, multiple genes and mechanisms have been identified and proposed as genetic causes of 46,XY DSD. In this study, 178 46,XY DSD patients were enrolled and underwent gene sequencing (either whole-exome sequencing or targeted panel gene sequencing). Detailed clinical phenotype and genotype information were summarized which showed that the most common clinical manifestations were micropenis (56.74%, 101/178), cryptorchidism (34.27%, 61/178), and hypospadias (17.42%, 31/178). Androgen synthesis/action disorders and idiopathic hypogonadotropic hypogonadism were the most frequent clinical diagnoses, accounting, respectively, for 40.90 and 21.59%. From all next-generation sequencing results, 103 candidate variants distributed across 32 genes were identified in 88 patients. The overall molecular detection rate was 49.44% (88/178), including 35.96% (64/178) pathogenic/likely pathogenic variants and 13.48% (24/178) variants of uncertain significance. Of all, 19.42% (20/103) variants were first reported in 46,XY DSD patients. Mutation c.680G>A (p.R227Q) on SRD5A2 (steroid 5-alpha-reductase 2) (36.67%, 11/30) was a hotspot mutation in the Chinese population. Novel candidate genes related to DSD (GHR (growth hormone receptor) and PHIP (pleckstrin homology domain-interacting protein)) were identified. Overall, this was a large cohort of 46,XY DSD patients with a common clinical classification and phenotype spectrum of Chinese patients. Targeted gene panel sequencing covered most of the genes contributing to DSD, whereas whole-exome sequencing detected more candidate genes.

Heterozygous rare variants in NR2F2 cause a recognizable multiple congenital anomaly syndrome with developmental delays

Nuclear receptor subfamily 2 group F member 2 (NR2F2 or COUP-TF2) encodes a transcription factor which is expressed at high levels during mammalian development. Rare heterozygous Mendelian variants in NR2F2 were initially identified in individuals with congenital heart disease (CHD), then subsequently in cohorts of congenital diaphragmatic hernia (CDH) and 46,XX ovotesticular disorders/differences of sexual development (DSD); however, the phenotypic spectrum associated with pathogenic variants in NR2F2 remains poorly characterized. Currently, less than 40 individuals with heterozygous pathogenic variants in NR2F2 have been reported. Here, we review the clinical and molecular details of 17 previously unreported individuals with rare heterozygous NR2F2 variants, the majority of which were de novo. Clinical features were variable, including intrauterine growth restriction (IUGR), CHD, CDH, genital anomalies, DSD, developmental delays, hypotonia, feeding difficulties, failure to thrive, congenital and acquired microcephaly, dysmorphic facial features, renal failure, hearing loss, strabismus, asplenia, and vascular malformations, thus expanding the phenotypic spectrum associated with NR2F2 variants. The variants seen were predicted loss of function, including a nonsense variant inherited from a mildly affected mosaic mother, missense and a large deletion including the NR2F2 gene. Our study presents evidence for rare, heterozygous NR2F2 variants causing a highly variable syndrome of congenital anomalies, commonly associated with heart defects, developmental delays/intellectual disability, dysmorphic features, feeding difficulties, hypotonia, and genital anomalies. Based on the new and previous cases, we provide clinical recommendations for evaluating individuals diagnosed with an NR2F2-associated disorder.

Lessons Learned from 17 Years of Multidisciplinary Care for Differences of Sex Development Patients at a Single Indonesian Center

BACKGROUND

Our multidisciplinary team (MDT) is a large specialized team based in Semarang, Indonesia, that cares for a wide variety of pediatric and adult individuals with differences of sex development (DSD) from across Indonesia. Here, we describe our work over the last 17 years.

METHODS

We analyzed phenotypic, hormonal, and genetic findings from clinical records for all patients referred to our MDT during the period 2004-2020.

RESULTS

Among 1,184 DSD patients, 10% had sex chromosome DSD, 67% had 46,XY DSD, and 23% had 46,XX DSD. The most common sex chromosome anomaly was Turner syndrome (45,X) (55 cases). For patients with 46,XY DSD under-masculinization was the most common diagnosis (311 cases), and for 46,XX DSD, a defect of Müllerian development was most common (131 cases) followed by congenital adrenal hyperplasia (CAH) (116 cases). Sanger sequencing, MLPA, and targeted gene sequencing of 257 patients with 46,XY DSD found likely causative variants in 21% (55 cases), with 13 diagnostic genes implicated. The most affected gene codes for the androgen receptor. Molecular analysis identified a diagnosis for 69 of 116 patients with CAH, with 62 carrying variants in CYP21A2 including four novel variants, and 7 patients carrying variants in CYP11B1. In many cases, these genetic diagnoses influenced the clinical management of patients and their families.

CONCLUSIONS

Our work has highlighted the occurrence of different DSDs in Indonesia. By applying sequencing technologies as part of our clinical care, we have delivered a number of genetic diagnoses and identified novel pathogenic variants in some genes, which may be clinically specific to Indonesia. Genetics can inform many aspects of DSD clinical management, and while many of our patients remain undiagnosed, we hope that future testing may provide answers for even more.

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

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

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.

Additional evidence for the role of chromosomal imbalances and SOX8, ZNRF3 and HHAT gene variants in early human testis development

BACKGROUND

Forty-six ,XY Differences/Disorders of Sex Development (DSD) are characterized by a broad phenotypic spectrum ranging from typical female to male with undervirilized external genitalia, or more rarely testicular regression with a typical male phenotype. Despite progress in the genetic diagnosis of DSD, most 46,XY DSD cases remain idiopathic.

METHODS

To determine the genetic causes of 46,XY DSD, we studied 165 patients of Tunisian ancestry, who presented a wide range of DSD phenotypes. Karyotyping, candidate gene sequencing, and whole-exome sequencing (WES) were performed.

RESULTS

Cytogenetic abnormalities, including a high frequency of sex chromosomal anomalies (85.4%), explained the phenotype in 30.9% (51/165) of the cohort. Sanger sequencing of candidate genes identified a novel pathogenic variant in the SRY gene in a patient with 46,XY gonadal dysgenesis. An exome screen of a sub-group of 44 patients with 46,XY DSD revealed pathogenic or likely pathogenic variants in 38.6% (17/44) of patients.

CONCLUSION

Rare or novel pathogenic variants were identified in the AR, SRD5A2, ZNRF3, SOX8, SOX9 and HHAT genes. Overall our data indicate a genetic diagnosis rate of 41.2% (68/165) in the group of 46,XY DSD.