Human sex reversal is caused by duplication or deletion of core enhancers upstream of SOX9

Two redundant transcription factor binding sites in a single enhancer are essential for mammalian sex determination

Male development in mammals depends on the activity of the two SOX gene: Sry and Sox9, in the embryonic testis. As deletion of Enhancer 13 (Enh13) of the Sox9 gene results in XY male-to-female sex reversal, we explored the critical elements necessary for its function and hence, for testis and male development. Here, we demonstrate that while microdeletions of individual transcription factor binding sites (TFBS) in Enh13 lead to normal testicular development, combined microdeletions of just two SRY/SOX binding motifs can alone fully abolish Enh13 activity leading to XY male-to-female sex reversal. This suggests that for proper male development to occur, these few nucleotides of non-coding DNA must be intact. Interestingly, we show that depending on the nature of these TFBS mutations, dramatically different phenotypic outcomes can occur, providing a molecular explanation for the distinct clinical outcomes observed in patients harboring different variants in the same enhancer.

Decoding the epigenetic mechanism of mammalian sex determination

Sex determination embodies a dynamic and intricate developmental process wielding significant influence over the destiny of bipotential gonads, steering them towards male or female gonads. Gonadal differentiation and the postnatal manifestation of the gonadal phenotype involve a sophisticated interplay of transcription factors such as SOX9 and FOXL2. Central to this interplay are chromatin modifiers regulating the mutual antagonism during this interplay. In this review, the key findings and knowledge gaps in DNA methylation, histone modification, and non-coding RNA-mediated control throughout mammalian gonadal development are covered. Furthermore, it explores the role of the developing brain in playing a pivotal role in the initiation of gonadogenesis and the subsequent involvement of gonadal hormone/hormone receptor in fine-tuning sexual differentiation. Based on promising facts, the role of the developing brain through the hypothalamic pituitary gonadal axis is explained and suggested as a novel hypothesis. The article also discusses the potential impact of ecological factors on the human epigenome in relation to sex determination and trans-generational epigenetics in uncovering novel genes and mechanisms involved in sex determination and gonadal differentiation. We have subtly emphasized the disruptions in epigenetic regulations contributing to sexual disorders, which further allows us to raise certain questions, decipher approaches for handling these questions and setting up the direction of future research.

Diagnosis and management of non-CAH 46,XX disorders/differences in sex development

Prenatal-onset androgen excess leads to abnormal sexual development in 46,XX individuals. This androgen excess can be caused endogenously by the adrenals or gonads or by exposure to exogenous androgens. The most common cause of 46,XX disorders/differences in sex development (DSD) is congenital adrenal hyperplasia (CAH) due to 21-hydroxylase deficiency, comprising >90% of 46,XX DSD cases. Deficiencies of 11β-hydroxylase, 3β-hydroxysteroid dehydrogenase, and P450-oxidoreductase (POR) are rare types of CAH, resulting in 46,XX DSD. In all CAH forms, patients have normal ovarian development. The molecular genetic causes of 46,XX DSD, besides CAH, are uncommon. These etiologies include primary glucocorticoid resistance (PGCR) and aromatase deficiency with normal ovarian development. Additionally, 46,XX gonads can differentiate into testes, causing 46,XX testicular (T) DSD or a coexistence of ovarian and testicular tissue, defined as 46,XX ovotesticular (OT)-DSD. PGCR is caused by inactivating variants in NR3C1, resulting in glucocorticoid insensitivity and the signs of mineralocorticoid and androgen excess. Pathogenic variants in the CYP19A1 gene lead to aromatase deficiency, causing androgen excess. Many genes are involved in the mechanisms of gonadal development, and genes associated with 46,XX T/OT-DSD include translocations of the SRY; copy number variants in NR2F2, NR0B1, SOX3, SOX9, SOX10, and FGF9, and sequence variants in NR5A1, NR2F2, RSPO1, SOX9, WNT2B, WNT4, and WT1. Progress in cytogenetic and molecular genetic techniques has significantly improved our understanding of the etiology of non-CAH 46,XX DSD. Nonetheless, uncertainties about gonadal function and gender outcomes may make the management of these conditions challenging. This review explores the intricate landscape of diagnosing and managing these conditions, shedding light on the unique aspects that distinguish them from other types of DSD.

The single-cell chromatin landscape in gonadal cell lineage specification

Gonad development includes sex determination and divergent maturation of the testes and ovaries. Recent advances in measuring gene expression in single cells are providing new insights into this complex process. However, the underlying epigenetic regulatory mechanisms remain unclear. Here, we profiled chromatin accessibility in mouse gonadal cells of both sexes from embryonic day 11.5 to 14.5 using single-cell assay for transposase accessible chromatin by sequencing (scATAC-seq). Our results showed that individual cell types can be inferred by the chromatin landscape, and that cells can be temporally ordered along developmental trajectories. Integrative analysis of transcriptomic and chromatin-accessibility maps identified multiple putative regulatory elements proximal to key gonadal genes Nr5a1, Sox9 and Wt1. We also uncover cell type-specific regulatory factors underlying cell type specification. Overall, our results provide a better understanding of the epigenetic landscape associated with the progressive restriction of cell fates in the gonad.

A conserved NR5A1-responsive enhancer regulates SRY in testis-determination

The Y-linked SRY gene initiates mammalian testis-determination. However, how the expression of SRY is regulated remains elusive. Here, we demonstrate that a conserved steroidogenic factor-1 (SF-1)/NR5A1 binding enhancer is required for appropriate SRY expression to initiate testis-determination in humans. Comparative sequence analysis of SRY 5' regions in mammals identified an evolutionary conserved SF-1/NR5A1-binding motif within a 250 bp region of open chromatin located 5 kilobases upstream of the SRY transcription start site. Genomic analysis of 46,XY individuals with disrupted testis-determination, including a large multigenerational family, identified unique single-base substitutions of highly conserved residues within the SF-1/NR5A1-binding element. In silico modelling and in vitro assays demonstrate the enhancer properties of the NR5A1 motif. Deletion of this hemizygous element by genome-editing, in a novel in vitro cellular model recapitulating human Sertoli cell formation, resulted in a significant reduction in expression of SRY. Therefore, human NR5A1 acts as a regulatory switch between testis and ovary development by upregulating SRY expression, a role that may predate the eutherian radiation. We show that disruption of an enhancer can phenocopy variants in the coding regions of SRY that cause human testis dysgenesis. Since disease causing variants in enhancers are currently rare, the regulation of gene expression in testis-determination offers a paradigm to define enhancer activity in a key developmental process.

FOXL2 interaction with different binding partners regulates the dynamics of ovarian development

The transcription factor FOXL2 is required in ovarian somatic cells for female fertility. Differential timing of Foxl2 deletion, in embryonic versus adult mouse ovary, leads to distinctive outcomes, suggesting different roles across development. Here, we comprehensively investigated FOXL2's role through a multi-omics approach to characterize gene expression dynamics and chromatin accessibility changes, coupled with genome-wide identification of FOXL2 targets and on-chromatin interacting partners in somatic cells across ovarian development. We found that FOXL2 regulates more targets postnatally, through interaction with factors regulating primordial follicle formation and steroidogenesis. Deletion of one interactor, ubiquitin-specific protease 7 (Usp7), results in impairment of somatic cell differentiation, germ cell nest breakdown, and ovarian development, leading to sterility. Our datasets constitute a comprehensive resource for exploration of the molecular mechanisms of ovarian development and causes of female infertility.

Reprograming skin fibroblasts into Sertoli cells: a patient-specific tool to understand effects of genetic variants on gonadal development

BACKGROUND

Disorders/differences of sex development (DSD) are congenital conditions in which the development of chromosomal, gonadal, or anatomical sex is atypical. With overlapping phenotypes and multiple genes involved, poor diagnostic yields are achieved for many of these conditions. The current DSD diagnostic regimen can be augmented by investigating transcriptome/proteome in vivo, but it is hampered by the unavailability of affected gonadal tissue at the relevant developmental stage. We try to mitigate this limitation by reprogramming readily available skin tissue-derived dermal fibroblasts into Sertoli cells (SC), which could then be deployed for different diagnostic strategies. SCs form the target cell type of choice because they act like an organizing center of embryonic gonadal development and many DSD arise when these developmental processes go awry.

METHODS

We employed a computational predictive algorithm for cell conversions called Mogrify to predict the transcription factors (TFs) required for direct reprogramming of human dermal fibroblasts into SCs. We established trans-differentiation culture conditions where stable transgenic expression of these TFs was achieved in 46, XY adult dermal fibroblasts using lentiviral vectors. The resulting Sertoli like cells (SLCs) were validated for SC phenotype using several approaches.

RESULTS

SLCs exhibited Sertoli-like morphological and cellular properties as revealed by morphometry and xCelligence cell behavior assays. They also showed Sertoli-specific expression of molecular markers such as SOX9, PTGDS, BMP4, or DMRT1 as revealed by IF imaging, RNAseq and qPCR. The SLC transcriptome shared about two thirds of its differentially expressed genes with a human adult SC transcriptome and expressed markers typical of embryonic SCs. Notably, SLCs lacked expression of most markers of other gonadal cell types such as Leydig, germ, peritubular myoid or granulosa cells.

CONCLUSIONS

The trans-differentiation method was applied to a variety of commercially available 46, XY fibroblasts derived from patients with DSD and to a 46, XX cell line. The DSD SLCs displayed altered levels of trans-differentiation in comparison to normal 46, XY-derived SLCs, thus showcasing the robustness of this new trans-differentiation model. Future applications could include using the SLCs to improve definitive diagnosis of DSD in patients with variants of unknown significance.

SOX on tumors, a comfort or a constraint?

The sex-determining region Y (SRY)-related high-mobility group (HMG) box (SOX) family, composed of 20 transcription factors, is a conserved family with a highly homologous HMG domain. Due to their crucial role in determining cell fate, the dysregulation of SOX family members is closely associated with tumorigenesis, including tumor invasion, metastasis, proliferation, apoptosis, epithelial-mesenchymal transition, stemness and drug resistance. Despite considerable research to investigate the mechanisms and functions of the SOX family, confusion remains regarding aspects such as the role of the SOX family in tumor immune microenvironment (TIME) and contradictory impacts the SOX family exerts on tumors. This review summarizes the physiological function of the SOX family and their multiple roles in tumors, with a focus on the relationship between the SOX family and TIME, aiming to propose their potential role in cancer and promising methods for treatment.

Clinical utility of early rapid genome sequencing in the evaluation of patients with differences of sex development

Establishing an early and accurate genetic diagnosis among patients with differences of sex development (DSD) is crucial in guiding the complex medical and psychosocial care they require. Genetic testing routinely utilized in clinical practice for this population is predicated upon physical exam findings and biochemical and endocrine profiling. This approach, however, is inefficient and unstandardized. Many patients with DSD, particularly those with 46,XY DSD, never receive a molecular genetic diagnosis. Rapid genome sequencing (rGS) is gaining momentum as a first-tier diagnostic instrument in the evaluation of patients with DSD given its ability to provide greater diagnostic yield and timely results. We present the case of a patient with nonbinary genitalia and systemic findings for whom rGS identified a novel variant of the WT1 gene and resulted in a molecular diagnosis within two weeks of life. This timeframe of diagnosis for syndromic DSD is largely unprecedented at our institution. Rapid GS expedited mobilization of a multidisciplinary medical team; enabled early understanding of clinical trajectory; informed planning of medical and surgical interventions; and guided individualized psychosocial support provided to the family. This case highlights the potential of early rGS in transforming the evaluation and care of patients with DSD.

Identification of a New Enhancer That Promotes Sox9 Expression by a Comparative Analysis of Mouse and Sry-Deficient Amami Spiny Rat

INTRODUCTION

Testis differentiation is initiated by the SRY gene on the Y chromosome in mammalian species. However, the Amami spiny rat, Tokudaia osimensis, lacks both the Y chromosome and the Sry gene and acquired a unique Sox9 regulatory mechanism via a male-specific duplication upstream of Sox9, without Sry. In general mammalian species, the SRY protein binds to a testis-specific enhancer to promote SOX9 gene expression. Several enhancers located upstream of Sox9/SOX9 have been reported in mice and humans. In particular, the binding of SRY to the highly conserved enhancer Enh13 is thought to be a common mechanism underlying testis differentiation and sex determination in mammals.

METHODS

Sequences of T. osimensis homologues of three Sox9 enhancers that were previously reported in mice, Enh8, Enh14, and Enh13, were determined. We performed in vitro assays to confirm enhancer activity involved in Sox9 regulation in T. osimensis.

RESULTS

T. osimensis Enh13 showed enhancer activity when co-transfected with NR5A1 and SOX9. Mouse Enh13 was activated by NR5A1 and SRY; however, T. osimensis Enh13 did not respond to SRY, even though the binding sites of SRY and NR5A1 were conserved. To identify the key sequence that is present in mouse but absent from T. osimensis, we performed reporter gene assays using vectors in which partial sequences of T. osimensis Enh13 were replaced with mouse sequences. For T. osimensis Enh13 in which the second half (approximately 430 bp) was replaced with the corresponding mouse sequence, activity in response to NR5A1 and SRY was recovered. Further, reporter assays revealed that multiple regions in the second half of the mouse Enh13 sequence are required for the response to NR5A1 and SRY. The latter 49 bp was particularly important and contained four binding sites for three transcription factors, POU2F1, HOXA3, and GATA1.

CONCLUSION

We showed that there are unknown sequences responsible for the interaction between NR5A1 and SRY and mEnh13 based on comparative analyses of Sry-dependent and Sry-independent species. Our comparative analyses revealed new molecular mechanisms underlying mammalian sex determination.

The smallest likely pathogenic duplication of a SOX9 enhancer identified to date in a family with 46,XX testicular differences of sex development

Copy number variants that duplicate distal upstream enhancer elements of the SOX9 gene cause 46,XX testicular differences of sex development (DSD) which is characterized by a 46,XX karyotype in an individual presenting with either ambiguous genitalia or genitalia with varying degrees of virilization, including those resembling typical male genitalia. Reported duplications in this region range in size from 24 to 780 kilobases (kb). Here we report a family with two affected individuals, the proband and his maternal uncle, harboring a 3.7 kb duplication of a SOX9 enhancer identified by clinical genome sequencing. Prior fluorescence in situ hybridization (FISH) for SRY and a multi-gene panel for ambiguous genitalia were non-diagnostic. The unaffected mother also carries this duplication, consistent with previously described incomplete penetrance. To our knowledge, this is the smallest duplication identified to-date, most of which resides in a 5.2 kb region that has been previously shown to possess enhancer activity that promotes the expression of SOX9. The duplication was confirmed by quantitative-PCR and shown to be in tandem by bidirectional Sanger sequencing breakpoint analysis. This finding highlights the importance of non-coding variant interrogation in suspected genetic disorders.

ERK1/2-SOX9/FOXL2 axis regulates ovarian steroidogenesis and favors the follicular-luteal transition

Estradiol and progesterone are the primary sex steroids produced by the ovary. Upon luteinizing hormone surge, estradiol-producing granulosa cells convert into progesterone-producing cells and eventually become large luteal cells of the corpus luteum. Signaling pathways and transcription factors involved in the cessation of estradiol and simultaneous stimulation of progesterone production in granulosa cells are not clearly understood. Here, we decipher that phosphorylated ERK1/2 regulates granulosa cell steroidogenesis by inhibiting estradiol and inducing progesterone production. Down-regulation of transcription factor FOXL2 and up-regulation of SOX9 by ERK underpin its differential steroidogenic function. Interestingly, the incidence of SOX9 is largely uncovered in ovarian cells and is found to regulate FOXL2 along with CYP19A1 and STAR genes, encoding rate-limiting enzymes of steroidogenesis, in cultured granulosa cells. We propose that the novel ERK1/2-SOX9/FOXL2 axis in granulosa cells is a critical regulator of ovarian steroidogenesis and may be considered when addressing pathophysiologies associated with inappropriate steroid production and infertility in humans and animals.

Three-dimensional genome structure and function

Linear DNA undergoes a series of compression and folding events, forming various three-dimensional (3D) structural units in mammalian cells, including chromosomal territory, compartment, topologically associating domain, and chromatin loop. These structures play crucial roles in regulating gene expression, cell differentiation, and disease progression. Deciphering the principles underlying 3D genome folding and the molecular mechanisms governing cell fate determination remains a challenge. With advancements in high-throughput sequencing and imaging techniques, the hierarchical organization and functional roles of higher-order chromatin structures have been gradually illuminated. This review systematically discussed the structural hierarchy of the 3D genome, the effects and mechanisms of cis-regulatory elements interaction in the 3D genome for regulating spatiotemporally specific gene expression, the roles and mechanisms of dynamic changes in 3D chromatin conformation during embryonic development, and the pathological mechanisms of diseases such as congenital developmental abnormalities and cancer, which are attributed to alterations in 3D genome organization and aberrations in key structural proteins. Finally, prospects were made for the research about 3D genome structure, function, and genetic intervention, and the roles in disease development, prevention, and treatment, which may offer some clues for precise diagnosis and treatment of related diseases.

Structural elements promote architectural stripe formation and facilitate ultra-long-range gene regulation at a human disease locus

Enhancer clusters overlapping disease-associated mutations in Pierre Robin sequence (PRS) patients regulate SOX9 expression at genomic distances over 1.25 Mb. We applied optical reconstruction of chromatin architecture (ORCA) imaging to trace 3D locus topology during PRS-enhancer activation. We observed pronounced changes in locus topology between cell types. Subsequent analysis of single-chromatin fiber traces revealed that these ensemble-average differences arise through changes in the frequency of commonly sampled topologies. We further identified two CTCF-bound elements, internal to the SOX9 topologically associating domain, which promote stripe formation, are positioned near the domain's 3D geometric center, and bridge enhancer-promoter contacts in a series of chromatin loops. Ablation of these elements results in diminished SOX9 expression and altered domain-wide contacts. Polymer models with uniform loading across the domain and frequent cohesin collisions recapitulate this multi-loop, centrally clustered geometry. Together, we provide mechanistic insights into architectural stripe formation and gene regulation over ultra-long genomic ranges.

Genetic control of typical and atypical sex development

Sex development relies on the sex-specific action of gene networks to differentiate the bipotential gonads of the growing fetus into testis or ovaries, followed by the differentiation of internal and external genitalia depending on the presence or absence of hormones. Differences in sex development (DSD) arise from congenital alterations during any of these processes, and are classified depending on sex chromosomal constitution as sex chromosome DSD, 46,XY DSD or 46,XX DSD. Understanding the genetics and embryology of typical and atypical sex development is essential for diagnosing, treating and managing DSD. Advances have been made in understanding the genetic causes of DSD over the past 10 years, especially for 46,XY DSD. Additional information is required to better understand ovarian and female development and to identify further genetic causes of 46,XX DSD, besides congenital adrenal hyperplasia. Ongoing research is focused on the discovery of further genes related to typical and atypical sex development and, therefore, on improving diagnosis of DSD.

FGF9 variant in 46,XY DSD patient suggests a role for dimerization in sex determination

46,XY gonadal dysgenesis (GD) is a Disorder/Difference of Sex Development (DSD) that can present with phenotypes ranging from ambiguous genitalia to complete male-to-female sex reversal. Around 50% of 46,XY DSD cases receive a molecular diagnosis. In mice, Fibroblast growth factor 9 (FGF9) is an important component of the male sex-determining pathway. Two FGF9 variants reported to date disrupt testis development in mice, but not in humans. Here, we describe a female patient with 46,XY GD harbouring the rare FGF9 variant (missense mutation), NM_002010.2:c.583G > A;p.(Asp195Asn) (D195N). By biochemical and cell-based approaches, the D195N variant disrupts FGF9 protein homodimerisation and FGF9-heparin-binding, and reduces both Sertoli cell proliferation and Wnt4 repression. XY Fgf9D195N/D195N foetal mice show a transient disruption of testicular cord development, while XY Fgf9D195N/- foetal mice show partial male-to-female gonadal sex reversal. In the general population, the D195N variant occurs at an allele frequency of 2.4 × 10-5 , suggesting an oligogenic basis for the patient's DSD. Exome analysis of the patient reveals several known and novel variants in genes expressed in human foetal Sertoli cells at the time of sex determination. Taken together, our results indicate that disruption of FGF9 homodimerization impairs testis determination in mice and, potentially, also in humans in combination with other variants.

Genomic structural variation: A complex but important driver of human evolution

Structural variants (SVs)-including duplications, deletions, and inversions of DNA-can have significant genomic and functional impacts but are technically difficult to identify and assay compared with single-nucleotide variants. With the aid of new genomic technologies, it has become clear that SVs account for significant differences across and within species. This phenomenon is particularly well-documented for humans and other primates due to the wealth of sequence data available. In great apes, SVs affect a larger number of nucleotides than single-nucleotide variants, with many identified SVs exhibiting population and species specificity. In this review, we highlight the importance of SVs in human evolution by (1) how they have shaped great ape genomes resulting in sensitized regions associated with traits and diseases, (2) their impact on gene functions and regulation, which subsequently has played a role in natural selection, and (3) the role of gene duplications in human brain evolution. We further discuss how to incorporate SVs in research, including the strengths and limitations of various genomic approaches. Finally, we propose future considerations in integrating existing data and biospecimens with the ever-expanding SV compendium propelled by biotechnology advancements.

SOX3 duplication in a boy with 46,XX ovotesticular disorder of sex development and his 46,XX sister with atypical genitalia: Probable germline mosaicism

Ovotesticular disorders of sex development (OT-DSD) are characterized by ovarian follicles and seminiferous tubules in the same individual, with a wide range of atypical genitalia. We report on two sibs with atypical genitalia and SRY-negative 46,XX DSD, OT-DSD was confirmed only in the boy, while the girl had bilateral ovaries. Chromosome microarray analysis (CMA) showed a 737-kb duplication at Xq27.1 including the entire SOX3 gene in both sibs, which was confirmed by quantitative real time PCR. Also, X chromosome inactivation assay showed random inactivation in both sibs. Whole exome sequencing revealed no pathogenic or likely pathogenic variant. CMA of the parents showed normal results for both, suggesting that germline mosaicism could be the reason of recurrence of this duplication in the siblings. Our results support a pathogenic role of SOX3 overexpression in 46,XX subjects leading to variable DSD phenotypes.

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.

SOX9 and SRY binding sites on mouse mXYSRa/Enh13 enhancer redundantly regulate Sox9 expression to varying degrees

Sox9 plays an essential role in mammalian testis formation. It has been reported that gene expression in the testes is regulated by enhancers. Among them, mXYSRa/Enh13-which is located at far upstream of the transcription start site-plays a critical role, wherein its deletion causes complete male-to-female sex reversal in mice. It has been proposed that the binding sites (BSs) of SOX9 and SRY, the latter of which is the sex determining gene on the Y chromosome, are associated with mXYSRa/Enh13. They function as an enhancer, whereby the sequences are evolutionarily conserved and in vivo binding of SOX9 and SRY to mXYSRa/Enh13 has been demonstrated previously. However, their precise in vivo functions have not been examined to date. To this end, this study generated mice with substitutions on the SOX9 and SRY BSs to reveal their in vivo functions. Homozygous mutants of SOX9 and SRY BS were indistinguishable from XY males, whereas double mutants had small testes, suggesting that these functions are redundant and that there is another functional sequence on mXYSRa/Enh13, since mXYSRa/Enh13 deletion mice are XY females. In addition, the majority of hemizygous mice with substitutions in SOX9 BS and SRY BS were female and male, respectively, suggesting that SOX9 BS contributes more to SRY BS for mXYSRa/Enh13 to function. The additive effect of SOX9 and SRY via these BSs was verified using an in vitro assay. In conclusion, SOX9 BS and SRY BS function redundantly in vivo, and at least one more functional sequence should exist in mXYSRa/Enh13.

Turnover of mammal sex chromosomes in the Sry-deficient Amami spiny rat is due to male-specific upregulation of Sox9

Mammalian sex chromosomes are highly conserved, and sex is determined by SRY on the Y chromosome. Two exceptional rodent groups in which some species lack a Y chromosome and Sry offer insights into how novel sex genes can arise and replace Sry, leading to sex chromosome turnover. However, intensive study over three decades has failed to reveal the identity of novel sex genes in either of these lineages. We here report our discovery of a male-specific duplication of an enhancer of Sox9 in the Amami spiny rat Tokudaia osimensis, in which males and females have only a single X chromosome (XO/XO) and the Y chromosome and Sry are completely lost. We performed a comprehensive survey to detect sex-specific genomic regions in the spiny rat. Sex-related genomic differences were limited to a male-specific duplication of a 17-kb unit located 430 kb upstream of Sox9 on an autosome. Hi-C analysis using male spiny rat cells showed the duplicated region has potential chromatin interaction with Sox9. The duplicated unit harbored a 1,262-bp element homologous to mouse enhancer 14 (Enh14), a candidate Sox9 enhancer that is functionally redundant in mice. Transgenic reporter mice showed that the spiny rat Enh14 can function as an embryonic testis enhancer in mice. Embryonic gonads of XX mice in which Enh14 was replaced by the duplicated spiny rat Enh14 showed increased Sox9 expression and decreased Foxl2 expression. We propose that male-specific duplication of this Sox9 enhancer substituted for Sry function, defining a novel Y chromosome in the spiny rat.

Modeling Human Gonad Development in Organoids

BACKGROUND

Our learning about human reproductive development is greatly hampered due to the absence of an adequate model. Animal studies cannot truthfully recapitulate human developmental processes, and studies of human fetal tissues are limited by their availability and ethical restrictions. Innovative three-dimensional (3D) organoid technology utilizing human pluripotent stem cells (hPSCs) offered a new approach to study tissue and organ development in vitro. However, a system for modeling human gonad development has not been established, thus, limiting our ability to study causes of infertility.

METHODS

In our study we utilized the 3D hPSC organoid culture in mini-spin bioreactors. Relying on intrinsic self-organizing and differentiation capabilities of stem cells, we explored whether organoids could mimic the development of human embryonic and fetal gonad.

RESULTS

We have developed a simple, bioreactor-based organoid system for modeling early human gonad development. Male hPSC-derived organoids follow the embryonic gonad developmental trajectory and differentiate into multipotent progenitors, which further specialize into testicular supporting and interstitial cells. We demonstrated functional activity of the generated cell types by analyzing the expression of cell type-specific markers. Furthermore, the specification of gonadal progenitors in organoid culture was accompanied by the characteristic architectural tissue organization.

CONCLUSION

This organoid system opens the opportunity for detailed studies of human gonad and germ cell development that can advance our understanding of sex development disorders. Implementation of human gonad organoid technology could be extended to modeling causes of infertility and regenerative medicine applications.

Fadrozole-mediated sex reversal in the embryonic chicken gonad involves a PAX2 positive undifferentiated supporting cell state

Gonadal sex differentiation among vertebrates involves divergent fates of a common group of progenitor cells present in both presumptive ovaries and testes. The first cell type to differentiate gives rise to pre-Sertoli cells in the testis, and pre-follicular cells in the ovary. These cells derive from a common lineage of so-called “supporting cells”. In birds and other egg-laying vertebrates, locally synthesised estrogen has a central role in ovarian development and influences the fate of these supporting cells. Manipulation of estrogen levels during embryonic development induces gonadal sex reversal, providing an experimental setting to evaluate the process of gonadal sex differentiation. Recently, we identified PAX2 as a novel marker of the undifferentiated supporting cell lineage in the chicken embryo, expressed in both sexes prior to overt gonadal sex differentiation. PAX2 expression is downregulated at the onset of gonadal sex differentiation in both males and females. The analysis of this undifferentiated supporting cell marker, together with Sertoli (male) and pre-granulosa (female) will enhance our understanding of supporting cell differentiation. Here we characterized the supporting cells differentiation process and identified undifferentiated supporting cells in estrogen-mediated sex reversal experiments. Female embryos treated with the aromatase inhibitor fadrozole developed into ovotestis, containing pre-granulosa cells, Sertoli cells and PAX2 positive undifferentiated supporting cells. In contrast, male embryos treated with 17β-estradiol showed no PAX2+ undifferentiated gonadal supporting cells. Fadrozole time-course as well as multiple dose analysis suggests that supporting cell transdifferentiation involves a dedifferentiation event into a PAX2+ undifferentiated supporting cell state, followed by a redifferentiation towards the opposite sex lineage.

SOX9 in organogenesis: shared and unique transcriptional functions

The transcription factor SOX9 is essential for the development of multiple organs including bone, testis, heart, lung, pancreas, intestine and nervous system. Mutations in the human SOX9 gene led to campomelic dysplasia, a haploinsufficiency disorder with several skeletal malformations frequently accompanied by 46, XY sex reversal. The mechanisms underlying the diverse SOX9 functions during organ development including its post-translational modifications, the availability of binding partners, and tissue-specific accessibility to target gene chromatin. Here we summarize the expression, activities, and downstream target genes of SOX9 in molecular genetic pathways essential for organ development, maintenance, and function. We also provide an insight into understanding the mechanisms that regulate the versatile roles of SOX9 in different organs.

MYRF: A New Regulator of Cardiac and Early Gonadal Development—Insights from Single Cell RNA Sequencing Analysis

De novo variants in the myelin regulatory factor (MYRF), a transcription factor involved in the differentiation of oligodendrocytes, have been linked recently to the cardiac and urogenital syndrome, while familiar variants are associated with nanophthalmos. Here, we report for the first time on a patient with a de novo stop-gain variant in MYRF (p.Q838*) associated with Scimitar syndrome, 46,XY partial gonadal dysgenesis (GD) and severe hyperopia. Since variants in MYRF have been described in both 46,XX and 46,XY GD, we assumed a role of MYRF in the early development of the bipotential gonad. We used publicly available single cell sequencing data of human testis and ovary from different developmental stages and analysed them for MYRF expression. We identified MYRF expression in the subset of coelomic epithelial cells at stages of gonadal ridge development in 46,XX and 46,XY individuals. Differential gene expression analysis revealed significantly upregulated genes. Within these, we identified CITED2 as a gene containing a MYRF binding site. It has been shown that Cited2^−/− mice have gonadal defects in both testis and ovary differentiation, as well as defects in heart development and establishment of the left–right axis. This makes MYRF a potential candidate as an early regulator of gonadal and heart development via upregulation of the transcriptional cofactor CITED2.

A Ctnnb1 enhancer regulates neocortical neurogenesis by controlling the abundance of intermediate progenitors

β-catenin-dependent canonical Wnt signaling plays a plethora of roles in neocortex (Ncx) development, but its function in regulating the abundance of intermediate progenitors (IPs) is elusive. Here we identified neCtnnb1, an evolutionarily conserved cis-regulatory element with typical enhancer features in developing Ncx. neCtnnb1 locates 55 kilobase upstream of and spatially close to the promoter of Ctnnb1, the gene encoding β-catenin. CRISPR/Cas9-mediated activation or interference of the neCtnnb1 locus enhanced or inhibited transcription of Ctnnb1. neCtnnb1 drove transcription predominantly in the subventricular zone of developing Ncx. Knock-out of neCtnnb1 in mice resulted in compromised expression of Ctnnb1 and the Wnt reporter in developing Ncx. Importantly, knock-out of neCtnnb1 lead to reduced production and transit-amplification of IPs, which subsequently generated fewer upper-layer Ncx projection neurons (PNs). In contrast, enhancing the canonical Wnt signaling by stabilizing β-catenin in neCtnnb1-active cells promoted the production of IPs and upper-layer Ncx PNs. ASH2L was identified as the key trans-acting factor that associates with neCtnnb1 and Ctnnb1’s promoter to maintain Ctnnb1’s transcription in both mouse and human Ncx progenitors. These findings advance understanding of transcriptional regulation of Ctnnb1, and provide insights into mechanisms underlying Ncx expansion during development.

DHX37 and 46,XY DSD: A New Ribosomopathy?

Recently, a series of recurrent missense variants in the RNA-helicase DHX37 have been reported associated with either 46,XY gonadal dysgenesis, 46,XY testicular regression syndrome (TRS), or anorchia. All affected children have non-syndromic forms of disorders/differences of sex development (DSD). These variants, which involve highly conserved amino acids within known functional domains of the protein, are predicted by in silico tools to have a deleterious effect on helicase function. DHX37 is required for ribosome biogenesis in eukaryotes, and how these variants cause DSD is unclear. The relationship between DHX37 and human congenital disorders is complex as compound heterozygous as well as de novo heterozygous missense variants in DHX37 are also associated with a complex congenital developmental syndrome (NEDBAVC, neurodevelopmental disorder with brain anomalies and with or without vertebral or cardiac anomalies; OMIM 618731), consisting of microcephaly, global developmental delay, seizures, facial dysmorphia, and kidney and cardiac anomalies. Here, we will give a brief overview of ribosome biogenesis and the role of DHX37 in this process. We will discuss variants in DHX37, their contribution to human disease in the general context of human ribosomopathies, and the possible disease mechanisms that may be involved.

Genome-wide identification of functional enhancers and their potential roles in pig breeding

Background

The pig is an economically important livestock species and is a widely applied large animal model in medical research. Enhancers are critical regulatory elements that have fundamental functions in evolution, development and disease. Genome-wide quantification of functional enhancers in the pig is needed.

Results

We performed self-transcribing active regulatory region sequencing (STARR-seq) in the porcine kidney epithelial PK15 and testicular ST cell lines, and reliably identified 2576 functional enhancers. Most of these enhancers were located in repetitive sequences and were enriched within silent and lowly expressed genes. Enhancers poorly overlapped with chromatin accessibility regions and were highly enriched in chromatin with the repressive histone modification H3K9me3, which is different from predicted pig enhancers detected using ChIP-seq for H3K27ac or/and H3K4me1 modified histones. This suggests that most pig enhancers identified with STARR-seq are endogenously repressed at the chromatin level and may function during cell type-specific development or at specific developmental stages. Additionally, the PPP3CA gene is associated with the loin muscle area trait and the QKI gene is associated with alkaline phosphatase activity that may be regulated by distal functional enhancers.

Conclusions

In summary, we generated the first functional enhancer map in PK15 and ST cells for the pig genome and highlight its potential roles in pig breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-022-00726-y.

SOX Genes and Their Role in Disorders of Sex Development

SOX genesare master regulatory genes controlling development and are fundamental to the establishment of sex determination in a multitude of organisms. The discovery of the master sex-determining gene SRY in 1990 was pivotal for the understanding of how testis development is initiated in mammals. With this discovery, an entire family of SOX factors were uncovered that play crucial roles in cell fate decisions during development. The importance of SOX genes in human reproductive development is evident from the various disorders of sex development (DSD) upon loss or overexpression of SOX gene function. Here, we review the roles that SOX genes play in gonad development and their involvement in DSD. We start with an overview of sex determination and differentiation, DSDs, and the SOX gene family and function. We then provide detailed information and discussion on SOX genes that have been implicated in DSDs, both at the gene and regulatory level. These include SRY, SOX9, SOX3, SOX8, and SOX10. This review provides insights on the crucial balance of SOX gene expression levels needed for gonad development and maintenance and how changes in these levels can lead to DSDs.

Genetic variation in CSF2 (5q31.1) is associated with longitudinal susceptibility to pediatric malaria, severe malarial anemia, and all-cause mortality in a high-burden malaria and HIV region of Kenya

Plasmodium falciparum infections remain among the leading causes of morbidity and mortality in holoendemic transmission areas. Located within region 5q31.1, the colony-stimulating factor 2 gene (CSF2) encodes granulocyte–macrophage colony-stimulating factor (GM-CSF), a hematopoietic growth factor that mediates host immune responses. Since the effect of CSF2 variation on malaria pathogenesis remains unreported, we investigated the impact of two genetic variants in the 5q31.1 gene region flanking CSF2:g-7032 G > A (rs168681:G > A) and CSF2:g.64544T > C (rs246835:T > C) on the rate and timing of malaria and severe malarial anemia (SMA, Hb < 5.0 g/dL) episodes over 36 months of follow-up. Children (n = 1654, aged 2–70 months) were recruited from a holoendemic P. falciparum transmission area of western Kenya. Decreased incidence rate ratio (IRR) for malaria was conferred by inheritance of the CSF2:g.64544 TC genotype (P = 0.0277) and CSF2 AC/GC diplotype (P = 0.0015). Increased IRR for malaria was observed in carriers of the CSF2 AT/GC diplotype (P = 0.0237), while the inheritance of the CSF2 AT haplotype increased the IRR for SMA (P = 0.0166). A model estimating the longitudinal risk of malaria showed decreased hazard rates among CSF2 AC haplotype carriers (P = 0.0045). Investigation of all-cause mortality revealed that inheritance of the GA genotype at CSF2:g-7032 increased the risk of mortality (P = 0.0315). Higher risk of SMA and all-cause mortality were observed in younger children (P < 0.0001 and P = 0.0015), HIV-1(+) individuals (P < 0.0001 and P < 0.0001), and carriers of HbSS (P = 0.0342 and P = 0.0019). Results from this holoendemic P. falciparum area show that variation in gene region 5q31.1 influences susceptibility to malaria, SMA, and mortality, as does age, HIV-1 status, and inheritance of HbSS.

Construction of Copy Number Variation Map Identifies Small Regions of Overlap and Candidate Genes for Atypical Female Genitalia Development

Copy number variations (CNVs) have been implicated in various conditions of differences of sexual development (DSD). Generally, larger genomic aberrations are more often considered disease-causing or clinically relevant, but over time, smaller CNVs have been associated with various forms of DSD. The main objective of this study is to identify small CNVs and the smallest regions of overlap (SROs) in patients with atypical female genitalia (AFG) and build a CNV map of AFG. We queried the DECIPHER database for recurrent duplications and/or deletions detected across the genome of AFG individuals. From these data, we constructed a chromosome map consisting of SROs and investigated such regions for genes that may be associated with the development of atypical female genitalia. Our study identified 180 unique SROs (7.95 kb to 45.34 Mb) distributed among 22 chromosomes. The most SROs were found in chromosomes X, 17, 11, and 22. None were found in chromosome 3. From these SROs, we identified 22 genes as potential candidates. Although none of these genes are currently associated with AFG, a literature review indicated that almost half were potentially involved in the development and/or function of the reproductive system, and only one gene was associated with a disorder that reported an individual patient with ambiguous genitalia. Our data regarding novel SROs requires further functional investigation to determine the role of the identified candidate genes in the development of atypical female genitalia, and this paper should serve as a catalyst for downstream molecular studies that may eventually affect the genetic counseling, diagnosis, and management of these DSD patients.

Deciphering Sex-Specific Differentiation of Human Fetal Gonads: Insight From Experimental Models

Sex-specific gonadal differentiation is initiated by the expression of SRY in male foetuses. This promotes a signalling pathway directing testicular development, while in female foetuses the absence of SRY and expression of pro-ovarian factors promote ovarian development. Importantly, in addition to the initiation of a sex-specific signalling cascade the opposite pathway is simultaneously inhibited. The somatic cell populations within the gonads dictates this differentiation as well as the development of secondary sex characteristics via secretion of endocrine factors and steroid hormones. Opposing pathways SOX9/FGF9 (testis) and WNT4/RSPO1 (ovary) controls the development and differentiation of the bipotential mouse gonad and even though sex-specific gonadal differentiation is largely considered to be conserved between mice and humans, recent studies have identified several differences. Hence, the signalling pathways promoting early mouse gonad differentiation cannot be directly transferred to human development thus highlighting the importance of also examining this signalling in human fetal gonads. This review focus on the current understanding of regulatory mechanisms governing human gonadal sex differentiation by combining knowledge of these processes from studies in mice, information from patients with differences of sex development and insight from manipulation of selected signalling pathways in ex vivo culture models of human fetal gonads.

Sequence-based modeling of three-dimensional genome architecture from kilobase to chromosome scale

To learn how genomic sequence influences multiscale three-dimensional (3D) genome architecture, this manuscript presents a sequence-based deep-learning approach, Orca, that predicts directly from sequence the 3D genome architecture from kilobase to whole-chromosome scale. Orca captures the sequence dependencies of structures including chromatin compartments and topologically associating domains, as well as diverse types of interactions from CTCF-mediated to enhancer-promoter interactions and Polycomb-mediated interactions with cell-type specificity. Orca enables various applications including predicting structural variant effects on multiscale genome organization and it recapitulated effects of experimentally studied variants at varying sizes (300 bp to 90 Mb). Moreover, Orca enables in silico virtual screens to probe the sequence basis of 3D genome organization at different scales. At the submegabase scale, it predicted specific transcription factor motifs underlying cell-type-specific genome interactions. At the compartment scale, virtual screens of sequence activities suggest a model for the sequence basis of chromatin compartments with a prominent role of transcription start sites.

Identification of Small Regions of Overlap from Copy Number Variable Regions in Patients with Hypospadias

Hypospadias is a common form of congenital atypical sex development that is often associated with other congenital comorbidities. Many genes have been associated with the condition, most commonly single sequence variations. Further investigations of recurrent and overlapping copy number variations (CNVs) have resulted in the identification of genes and chromosome regions associated with various conditions, including differences of sex development (DSD). In this retrospective study, we investigated the DECIPHER database, as well as an internal institutional database, to identify small recurrent CNVs among individuals with isolated and syndromic hypospadias. We further investigated these overlapping recurrent CNVs to identify 75 smallest regions of overlap (SROs) on 18 chromosomes. Some of the genes within these SROs may be considered potential candidate genes for the etiology of hypospadias and, occasionally, additional comorbid phenotypes. This study also investigates for the first time additional common phenotypes among individuals with hypospadias and overlapping CNVs. This study provides data that may aid genetic counseling and management of individuals with hypospadias, as well as improve understanding of its underlying genetic etiology and human genital development overall.

Expression and Potential Prognostic Value of SOX9, MCL-1 and SPOCK1 in Gastric Adenocarcinoma

Gastric cancer is a common malignancy and remains one of the leading causes of cancer-related deaths, though its incidence is in decline in most developed countries. One of the major challenges of treating gastric cancer is tumor heterogeneity, which portends a high degree of prognostic variance and the necessity for different treatment modalities. Tumor heterogeneity is at least in part due to divergent differentiation of tumor cells to clones harboring different molecular alterations. Here we studied the expression of emerging prognostic markers SOX9, MCL-1, and SPOCK1 (Testican-1) in a cohort of gastric cancer by immunohistochemistry and investigated how individual biomarkers and their combinations predict disease prognosis. We found frequent expression of SPOCK1 (in both nuclei and cytoplasm), MCL-1 and SOX9 in gastric cancer. In univariate analysis, nuclear SPOCK1 expression and pathologic TNM stage were negative prognostic markers in this cohort. In multivariate analysis, SOX9 expression stood out as a predictor of poor prognosis. Further subgroup analysis suggested prognostic value of SOX9 expression in poorly differentiated gastric adenocarcinoma. MCL-1 showed no prognostic role in this cohort.

H3K27ac chromatin acetylation and gene expression analysis reveal sex- and situs-related differences in developing chicken gonads

BACKGROUND

Birds exhibit a unique asymmetry in terms of gonad development. The female left gonad generates a functional ovary, whereas the right gonad regresses. In males, both left and right gonads would develop into testes. How is this left/right asymmetry established only in females but not in males remains unknown. The epigenetic regulation of gonadal developmental genes may contribute to this sex disparity. The modification of histone tails such as H3K27ac is tightly coupled to chromatin activation and gene expression. To explore whether H3K27ac marked chromatin activation is involved in the asymmetric development of avian gonads, we probed genome-wide H3K27ac occupancy in left and right gonads from both sexes and related chromatin activity profile to the expression of gonadal genes. Furthermore, we validated the effect of chromatin activity on asymmetric gonadal development by manipulating the chromatin histone acetylation levels.

METHODS

The undifferentiated gonads from both sides of each sex were collected and subjected to RNA-Seq and H3K27ac ChIP-Seq experiments. Integrated analysis of gene expression and active chromatin regions were performed to identify the sex- and situs-specific regulation and expression of gonadal genes. The histone deacetylase inhibitor trichostatin A (TSA) was applied to the undifferentiated female right gonads to assess the effect of chromatin activation on gonadal gene expression and cell proliferation.

RESULTS

Even before sex differentiation, the gonads already show divergent gene expression between different sexes and between left/right sides in females. The sex-specific H3K27ac chromatin distributions coincide with the higher expression of male/female specification genes in each sex. Unexpectedly, the H3K27ac marked chromatin activation show a dramatic difference between left and right gonads in both sexes, although the left/right asymmetric gonadal development was observed only in females but not in males. In females, the side-specific H3K27ac occupancy instructs the differential expression of developmental genes between the pair of gonads and contributes to the development of left but not right gonad. However, in males, the left/right discrepancy of H3K27ac chromatin distribution does not drive the side-biased gene expression or gonad development. The TSA-induced retention of chromatin acetylation causes up-regulation of ovarian developmental genes and increases cell proliferation in the female right gonad.

CONCLUSIONS

We revealed that left/right asymmetry in H3K27ac marked chromatin activation exists in both sexes, but this discrepancy gives rise to asymmetric gonadal development only in females. Other mechanisms overriding the chromatin activation would control the symmetric development of male gonads in chicken.

Developmental disorders caused by haploinsufficiency of transcriptional regulators: a perspective based on cell fate determination

Many human birth defects and neurodevelopmental disorders are caused by loss-of-function mutations in a single copy of transcription factor (TF) and chromatin regulator genes. Although this dosage sensitivity has long been known, how and why haploinsufficiency (HI) of transcriptional regulators leads to developmental disorders (DDs) is unclear. Here I propose the hypothesis that such DDs result from defects in cell fate determination that are based on disrupted bistability in the underlying gene regulatory network (GRN). Bistability, a crucial systems biology concept to model binary choices such as cell fate decisions, requires both positive feedback and ultrasensitivity, the latter often achieved through TF cooperativity. The hypothesis explains why dosage sensitivity of transcriptional regulators is an inherent property of fate decisions, and why disruption of either positive feedback or cooperativity in the underlying GRN is sufficient to cause disease. I present empirical and theoretical evidence in support of this hypothesis and discuss several issues for which it increases our understanding of disease, such as incomplete penetrance. The proposed framework provides a mechanistic, systems-level explanation of HI of transcriptional regulators, thus unifying existing theories, and offers new insights into outstanding issues of human disease. This article has an associated Future Leader to Watch interview with the author of the paper.

OUP accepted manuscript

The pathways of gametogenesis encompass elaborate cellular specialization accompanied by precise partitioning of the genome content in order to produce fully matured spermatozoa and oocytes. Transcription factors are an important class of molecules that function in gametogenesis to regulate intrinsic gene expression programs, play essential roles in specifying (or determining) germ cell fate and assist in guiding full maturation of germ cells and maintenance of their populations. Moreover, in order to reinforce or redirect cell fate in vitro, it is transcription factors that are most frequently induced, over-expressed or activated. Many reviews have focused on the molecular development and genetics of gametogenesis, in vivo and in vitro, in model organisms and in humans, including several recent comprehensive reviews: here, we focus specifically on the role of transcription factors. Recent advances in stem cell biology and multi-omic studies have enabled deeper investigation into the unique transcriptional mechanisms of human reproductive development. Moreover, as methods continually improve, in vitro differentiation of germ cells can provide the platform for robust gain- and loss-of-function genetic analyses. These analyses are delineating unique and shared human germ cell transcriptional network components that, together with somatic lineage specifiers and pluripotency transcription factors, function in transitions from pluripotent stem cells to gametes. This grand theme review offers additional insight into human infertility and reproductive disorders that are linked predominantly to defects in the transcription factor networks and thus may potentially contribute to the development of novel treatments for infertility.

Tenuous Transcriptional Threshold of Human Sex Determination. I. SRY and Swyer Syndrome at the Edge of Ambiguity

Male sex determination in mammals is initiated by SRY, a Y-encoded transcription factor. The protein contains a high-mobility-group (HMG) box mediating sequence-specific DNA bending. Mutations causing XY gonadal dysgenesis (Swyer syndrome) cluster in the box and ordinarily arise de novo. Rare inherited variants lead to male development in one genetic background (the father) but not another (his sterile XY daughter). De novo and inherited mutations occur at an invariant Tyr adjoining the motif's basic tail (box position 72; Y127 in SRY). In SRY-responsive cell lines CH34 and LNCaP, de novo mutations Y127H and Y127C reduced SRY activity (as assessed by transcriptional activation of principal target gene Sox9) by 5- and 8-fold, respectively. Whereas Y127H impaired testis-specific enhancer assembly, Y127C caused accelerated proteasomal proteolysis; activity was in part rescued by proteasome inhibition. Inherited variant Y127F was better tolerated: its expression was unperturbed, and activity was reduced by only twofold, a threshold similar to other inherited variants. Biochemical studies of wild-type (WT) and variant HMG boxes demonstrated similar specific DNA affinities (within a twofold range), with only subtle differences in sharp DNA bending as probed by permutation gel electrophoresis and fluorescence resonance-energy transfer (FRET); thermodynamic stabilities of the free boxes were essentially identical. Such modest perturbations are within the range of species variation. Whereas our cell-based findings rationalize the de novo genotype-phenotype relationships, a molecular understanding of inherited mutation Y127F remains elusive. Our companion study uncovers cryptic biophysical perturbations suggesting that the para-OH group of Y127 anchors a novel water-mediated DNA clamp.

Genetics of 46,XY gonadal dysgenesis

In 46,XY men, testis is determined by a genetic network(s) that both promotes testis formation and represses ovarian development. Disruption of this process results in a lack of testis-determination and affected individuals present with 46,XY gonadal dysgenesis (GD), a part of the spectrum of Disorders/Differences of Sex Development/Determination (DSD). A minority of all cases of GD are associated with pathogenic variants in key players of testis-determination, SRY, SOX9, MAP3K1 and NR5A1. However, most of the cases remain unexplained. Recently, unbiased exome sequencing approaches have revealed new genes and loci that may cause 46,XY GD. We critically evaluate the evidence to support causality of these factors and describe how functional studies are continuing to improve our understanding of genotype-phenotype relationships in genes that are established causes of GD. As genomic data continues to be generated from DSD cohorts, we propose several recommendations to help interpret the data and establish causality.

Tenuous transcriptional threshold of human sex determination. II. SRY exploits water-mediated clamp at the edge of ambiguity

Y-encoded transcription factor SRY initiates male differentiation in therian mammals. This factor contains a high-mobility-group (HMG) box, which mediates sequence-specific DNA binding with sharp DNA bending. A companion article in this issue described sex-reversal mutations at box position 72 (residue 127 in human SRY), invariant as Tyr among mammalian orthologs. Although not contacting DNA, the aromatic ring seals the domain's minor wing at a solvent-exposed junction with a basic tail. A seeming paradox was posed by the native-like biochemical properties of inherited Swyer variant Y72F: its near-native gene-regulatory activity is consistent with the father's male development, but at odds with the daughter's XY female somatic phenotype. Surprisingly, aromatic rings (Y72, F72 or W72) confer higher transcriptional activity than do basic or polar side chains generally observed at solvated DNA interfaces (Arg, Lys, His or Gln). Whereas biophysical studies (time-resolved fluorescence resonance energy transfer and heteronuclear NMR spectroscopy) uncovered only subtle perturbations, dissociation of the Y72F complex was markedly accelerated relative to wild-type. Studies of protein-DNA solvation by molecular-dynamics (MD) simulations of an homologous high-resolution crystal structure (SOX18) suggest that Y72 para-OH anchors a network of water molecules at the tail-DNA interface, perturbed in the variant in association with nonlocal conformational fluctuations. Loss of the Y72 anchor among SRY variants presumably "unclamps" its basic tail, leading to (a) rapid DNA dissociation despite native affinity and (b) attenuated transcriptional activity at the edge of sexual ambiguity. Conservation of Y72 suggests that this water-mediated clamp operates generally among SRY and metazoan SOX domains.

Generation and mutational analysis of a transgenic mouse model of human SRY

SRY is the Y-chromosomal gene that determines male sex development in humans and most other mammals. After three decades of study, we still lack a detailed understanding of which domains of the SRY protein are required to engage the pathway of gene activity leading to testis development. Some insight has been gained from the study of genetic variations underlying differences/disorders of sex determination (DSD), but the lack of a system of experimentally generating SRY mutations and studying their consequences in vivo has limited progress in the field. To address this issue, we generated a mouse model carrying a human SRY transgene able to drive testis determination in XX mice. Using CRISPR-Cas9 gene editing, we generated novel genetic modifications in each of SRY's three domains (N-terminal, HMG box, and C-terminal) and performed a detailed analysis of their molecular and cellular effects on embryonic testis development. Our results provide new functional insights unique to human SRY and present a versatile and powerful system in which to functionally analyze variations of SRY including known and novel pathogenic variants found in DSD.

Labile sex chromosomes in the Australian freshwater fish family Percichthyidae

Sex-specific ecology has management implications, but rapid sex-chromosome turnover in fishes hinders sex-marker development for monomorphic species. We used annotated genomes and reduced-representation sequencing data for two Australian percichthyids, Macquarie perch Macquaria australasica and golden perch M. ambigua, and whole genome resequencing for 50 Macquarie perch of each sex, to identify sex-linked loci and develop an affordable sexing assay. In silico pool-seq tests of 1,492,004 Macquarie perch SNPs revealed that a 275-kb scaffold was enriched for gametologous loci. Within this scaffold, 22 loci were sex-linked in a predominantly XY system, with females being homozygous for the X-linked allele at all 22, and males having the Y-linked allele at >7. Seven XY-gametologous loci (all males, but no females, are heterozygous or homozygous for the male-specific allele) were within a 146-bp region. A PCR-RFLP sexing assay targeting one Y-linked SNP, tested in 66 known-sex Macquarie perch and two of each sex of three confamilial species, plus amplicon sequencing of 400 bp encompassing the 146-bp region, revealed that the few sex-linked positions differ between species and between Macquarie perch populations. This indicates sex-chromosome lability in Percichthyidae, supported by nonhomologous scaffolds containing sex-linked loci for Macquarie- and golden perches. The present resources facilitate genomic research in Percichthyidae, including formulation of hypotheses about candidate genes of interest such as transcription factor SOX1b that occurs in the 275-kb scaffold ~38 kb downstream of the 146-bp region containing seven XY-gametologous loci. Sex-linked markers will be useful for determining genetic sex in some populations and studying sex chromosome turnover.

Generation of an induced pluripotent stem cell line GZHMCi008-A derived from a patient with SRY-positive 46,XX testicular disorder of sex development

The testicular disorder of sex development (TDSD) is a rare condition characterized by a male appearance with a female karyotype. The most frequent cause of TDSD is misplacement of the sex determining region Y (SRY) gene on the X chromosome. Here, we report the generation of an induced pluripotent stem cell (iPSC) line from peripheral blood mononuclear cells of a patient with SRY-positive 46,XX TDSD. This cell line offers an unprecedented cellular model to investigate the profound manifestations like infertility of the male sex reversal patients, and serves as a useful tool to develop therapies for the disease.

A Novel WT1 Mutation Identified in a 46,XX Testicular/Ovotesticular DSD Patient Results in the Retention of Intron 9

Simple Summary

Disorders/differences of sexual development are very diverse. Among them is a condition characterized by the presence of testicular tissue in people with female chromosomes, which is typically manifested by male or ambiguous genitalia. While genetic counseling is beneficial for these people and their families, the genetic causes of these cases are only partially understood. We describe a new mutation in the WT1 gene that results in the presence of testicular tissue in a child with a female karyotype. We propose molecular mechanisms disrupted by this mutation. This finding widened our understanding of processes that govern sexual development and can be used to develop diagnostic tests for disorders/differences of sexual development.

Abstract

The 46,XX testicular DSD (disorder/difference of sexual development) and 46,XX ovotesticular DSD (46,XX TDSD and 46,XX OTDSD) phenotypes are caused by genetic rearrangements or point mutations resulting in imbalance between components of the two antagonistic, pro-testicular and pro-ovarian pathways; however, the genetic causes of 46,XX TDSD/OTDSD are not fully understood, and molecular diagnosis for many patients with the conditions is unavailable. Only recently few mutations in the WT1 (WT1 transcription factor; 11p13) gene were described in a group of 46,XX TDSD and 46,XX OTDSD individuals. The WT1 protein contains a DNA/RNA binding domain consisting of four zinc fingers (ZnF) and a three-amino acid (KTS) motif that is present or absent, as a result of alternative splicing, between ZnF3 and ZnF4 (±KTS isoforms). Here, we present a patient with 46,XX TDSD/OTDSD in whom whole exome sequencing revealed a heterozygous de novo WT1 c.1437A>G mutation within an alternative donor splice site which is used for −KTS WT1 isoform formation. So far, no mutation in this splice site has been identified in any patient group. We demonstrated that the mutation results in the retention of intron 9 in the mature mRNA of the 46,XX TDSD/OTDSD patient. In cases when the erroneous mRNA is translated, exclusively the expression of a truncated WT1 +KTS protein lacking ZnF4 and no −KTS protein occurs from the mutated allele of the patient. We discuss potential mechanisms and pathways which can be disturbed upon two conditions: Absence of Zn4F and altered +KTS/−KTS ratio.

Long-Range Regulation of Key Sex Determination Genes

The development of sexually dimorphic gonads is a unique process that starts with the specification of the bipotential genital ridges and culminates with the development of fully differentiated ovaries and testes in females and males, respectively. Research on sex determination has been mostly focused on the identification of sex determination genes, the majority of which encode for proteins and specifically transcription factors such as SOX9 in the testes and FOXL2 in the ovaries. Our understanding of which factors may be critical for sex determination have benefited from the study of human disorders of sex development (DSD) and animal models, such as the mouse and the goat, as these often replicate the same phenotypes observed in humans when mutations or chromosomic rearrangements arise in protein-coding genes. Despite the advances made so far in explaining the role of key factors such as SRY, SOX9, and FOXL2 and the genes they control, what may regulate these factors upstream is not entirely understood, often resulting in the inability to correctly diagnose DSD patients. The role of non-coding DNA, which represents 98% of the human genome, in sex determination has only recently begun to be fully appreciated. In this review, we summarize the current knowledge on the long-range regulation of 2 important sex determination genes, SOX9 and FOXL2, and discuss the challenges that lie ahead and the many avenues of research yet to be explored in the sex determination field.

Cis-Regulatory Control of Mammalian Sex Determination

Sex determination is the process by which an initial bipotential gonad adopts either a testicular or ovarian cell fate. The inability to properly complete this process leads to a group of developmental disorders classified as disorders of sex development (DSD). To date, dozens of genes were shown to play roles in mammalian sex determination, and mutations in these genes can cause DSD in humans or gonadal sex reversal/dysfunction in mice. However, exome sequencing currently provides genetic diagnosis for only less than half of DSD patients. This points towards a major role for the non-coding genome during sex determination. In this review, we highlight recent advances in our understanding of non-coding, cis-acting gene regulatory elements and discuss how they may control transcriptional programmes that underpin sex determination in the context of the 3-dimensional folding of chromatin. As a paradigm, we focus on the Sox9 gene, a prominent pro-male factor and one of the most extensively studied genes in gonadal cell fate determination.

Targeting the Non-Coding Genome for the Diagnosis of Disorders of Sex Development

Disorders of sex development (DSD) are a complex group of conditions with highly variable clinical phenotypes, most often caused by failure of gonadal development. DSD are estimated to occur in around 1.7% of all live births. Whilst the understanding of genes involved in gonad development has increased exponentially, approximately 50% of patients with a DSD remain without a genetic diagnosis, possibly implicating non-coding genomic regions instead. Here, we review how variants in the non-coding genome of DSD patients can be identified using techniques such as array comparative genomic hybridization (CGH) to detect copy number variants (CNVs), and more recently, whole genome sequencing (WGS). Once a CNV in a patient's non-coding genome is identified, putative regulatory elements such as enhancers need to be determined within these vast genomic regions. We will review the available online tools and databases that can be used to refine regions with potential enhancer activity based on chromosomal accessibility, histone modifications, transcription factor binding site analysis, chromatin conformation, and disease association. We will also review the current in vitro and in vivo techniques available to demonstrate the functionality of the identified enhancers. The review concludes with a clinical update on the enhancers linked to DSD.

Genetic Regulation of Avian Testis Development

As in other vertebrates, avian testes are the site of spermatogenesis and androgen production. The paired testes of birds differentiate during embryogenesis, first marked by the development of pre-Sertoli cells in the gonadal primordium and their condensation into seminiferous cords. Germ cells become enclosed in these cords and enter mitotic arrest, while steroidogenic Leydig cells subsequently differentiate around the cords. This review describes our current understanding of avian testis development at the cell biology and genetic levels. Most of this knowledge has come from studies on the chicken embryo, though other species are increasingly being examined. In chicken, testis development is governed by the Z-chromosome-linked DMRT1 gene, which directly or indirectly activates the male factors, HEMGN, SOX9 and AMH. Recent single cell RNA-seq has defined cell lineage specification during chicken testis development, while comparative studies point to deep conservation of avian testis formation. Lastly, we identify areas of future research on the genetics of avian testis development.