MYRF haploinsufficiency causes 46,XY and 46,XX disorders of sex development: bioinformatics consideration

MYRF: A unique transmembrane transcription factor- from proteolytic self-processing to its multifaceted roles in animal development

The Myelin Regulator Factor (MYRF) is a master regulator governing myelin formation and maintenance in the central nervous system. The conservation of MYRF across metazoans and its broad tissue expression suggest it has functions extending beyond the well-established role in myelination. Loss of MYRF results in developmental lethality in both invertebrates and vertebrates, and MYRF haploinsufficiency in humans causes MYRF-related Cardiac Urogenital Syndrome, underscoring its importance in animal development; however, these mechanisms are largely unexplored. MYRF, an unconventional transcription factor, begins embedded in the membrane and undergoes intramolecular chaperone mediated trimerization, which triggers self-cleavage, allowing its N-terminal segment with an Ig-fold DNA-binding domain to enter the nucleus for transcriptional regulation. Recent research suggests developmental regulation of cleavage, yet the mechanisms remain enigmatic. While some parts of MYRF's structure have been elucidated, others remain obscure, leaving questions about how these motifs are linked to its intricate processing and function.

The genetic spectrum of a Chinese series of patients with 46, XY disorders of the sex development

PURPOSE

The etiology of 46, XY disorders of sex development (46, XY DSD) is complex, and studies have shown that different series of patients with 46, XY DSD has different genetic spectrum. In this study, we aimed to investigate the underlying genetic etiology in a Chinese series of patients with 46, XY DSD by whole exome sequencing (WES).

METHODS

Seventy patients with 46, XY DSD were enrolled from the Peking Union Medical College Hospital (Beijing, China). The detailed clinical characteristics were evaluated, and peripheral blood was collected for WES to find the patients' rare variants (RVs) of genes related to 46, XY DSD. The clinical significance of the RVs was annotated according to American College of Medical Genetics and Genomics (ACMG) guidelines.

RESULTS

A total of 57 RVs from nine genes were identified in 56 patients with 46, XY DSD, which include 21 novel RVs and 36 recurrent RVs. Based on the American ACMG guidelines, 43 variants were classified as pathogenic(P) or likely pathogenic (LP) variants and 14 variants were defined as variants of uncertain significance (VUS). P or LP variants were identified in 64.3% (45/70) patients of the series. Thirty-nine, 14, and 4 RVs were involved in the process of androgen synthesis and action, testicular determination and developmental process, and syndromic 46, XY DSD, respectively. The top three genes most frequently affected to cause 46, XY DSD were AR, SRD5A2, and NR5A1. Seven patients were found harboring RVs of the 46, XY DSD pathogenic genes identified in recent years, namely DHX37 in four patients, MYRF in two patients, and PPP2R3C in one patient.

CONCLUSION

We identified 21 novel RVs of nine genes, which extended the genetic spectrum of 46, XY DSD pathogenic variants. Our study showed that 60% of the patients were caused by AR, SRD5A2 or NR5A1 P/LP variants. Therefore, polymerase chain reaction (PCR) amplification and Sanger sequencing of these three genes could be performed first to identify the pathogeny of the patients. For those patients whose pathogenic variants had not been found, whole-exome sequencing could be helpful in determining the etiology.

CRISPR-based knockout and base editing confirm the role of MYRF in heart development and congenital heart disease

High-throughput DNA sequencing studies increasingly associate DNA variants with congenital heart disease (CHD). However, functional modeling is a crucial prerequisite for translating genomic data into clinical care. We used CRISPR-Cas9-mediated targeting of 12 candidate genes in the vertebrate model medaka (Oryzias latipes), five of which displayed a novel cardiovascular phenotype spectrum in F0 (crispants): mapre2, smg7, cdc42bpab, ankrd11 and myrf, encoding a transcription factor recently linked to cardiac-urogenital syndrome. Our myrf mutant line showed particularly prominent embryonic cardiac defects recapitulating phenotypes of pediatric patients, including hypoplastic ventricle. Mimicking human mutations, we edited three sites to generate specific myrf single-nucleotide variants via cytosine and adenine base editors. The Glu749Lys missense mutation in the conserved intramolecular chaperon autocleavage domain fully recapitulated the characteristic myrf mutant phenotype with high penetrance, underlining the crucial function of this protein domain. The efficiency and scalability of base editing to model specific point mutations accelerate gene validation studies and the generation of human-relevant disease models.

Identification of a Novel Variant in Myelin Regulatory Growth Factor by Next-Generation Sequencing Led to the Detection of a Clinically Inapparent Congenital Heart Defect in a Patient with a 46,XY Disorder of Sex Development

In patients with 46,XY disorders of sex development (DSDs), next-generation sequencing (NGS) has high diagnostic efficiency. One contribution to this diagnostic approach is the possibility of applying reverse phenotyping when a variant in a gene associated with multiple organ hits is found. Our aim is to report a case of a patient with 46,XY DSDs in whom the identification of a novel variant in MYRF led to the detection of a clinically inapparent congenital heart defect. A full-term newborn presented with ambiguous genitalia, as follows: a 2 cm phallus, penoscrotal hypospadias, partially fused labioscrotal folds, an anogenital distance of 1.2 cm, and non-palpable gonads. The karyotype was 46,XY, serum testosterone and AMH were low, whereas LH and FSH were high, leading to the diagnosis of dysgenetic DSD. Whole exome sequencing identified a novel, heterozygous, nonsense variant in MYRF, classified as pathogenic according to the ACMG criteria. MYRF encodes a membrane-bound transcriptional factor expressed in several tissues associated with OCUGS syndrome (ophthalmic, cardiac, and urogenital anomalies). In the patient, oriented clinical assessment ruled out ophthalmic defects, but ultrasonography confirmed meso/dextrocardia. We report a novel MYRF variant in a patient with 46,XY DSDs, allowing us to identify a clinically inapparent congenital heart defect by reverse phenotyping.

Bioinformatics analysis and verification of hub genes in 46,XY, disorders of sexual development

CONTEXT

46,XY, disorders of sexual development (46,XY, DSD) is a congenital genetic disease whose pathogenesis is complex and clinical manifestations are diverse. The existing molecular research has often focused on single-centre sequencing data, instead of prediction based on big data.

AIMS

This work aimed to fully understand the pathogenesis of 46,XY, DSD, and summarise the key pathogenic genes.

METHODS

Firstly, the potential pathogenic genes were identified from public data. Secondly, bioinformatics was used to predict pathogenic genes, including hub gene analysis, protein-protein interaction (PPI) and function enrichment analysis. Lastly, the genomic DNA from two unrelated families were recruited, next-generation sequencing and Sanger sequencing were performed to verify the hub genes.

KEY RESULTS

A total of 161 potential pathogenic genes were selected from MGI and PubMed gene sets. The PPI network was built which included 144 nodes and 194 edges. MCODE 4 was selected from PPI which scored the most significant P -value. The top 15 hub genes were ranked and identified by Cytoscape. Furthermore, three variants were found on SRD5A2 gene by genome sequencing, which belonged to the prediction hub genes.

CONCLUSIONS

Our results indicate that occurrence of 46,XY, DSD is attributed to a variety of genes. Bioinformatics analysis can help us predict the hub genes and find the most core network MCODE model.

IMPLICATIONS

Bioinformatic predictions may provide a novel perspective on better understanding the pathogenesis of 46,XY, DSD.

Single-cell transcriptomic profiling redefines the origin and specification of early adrenogonadal progenitors

Adrenal cortex and gonads represent the two major steroidogenic organs in mammals. Both tissues are considered to share a common developmental origin characterized by the expression of Nr5a1/Sf1. The precise origin of adrenogonadal progenitors and the processes driving differentiation toward the adrenal or gonadal fate remain, however, elusive. Here, we provide a comprehensive single-cell transcriptomic atlas of early mouse adrenogonadal development including 52 cell types belonging to twelve major cell lineages. Trajectory reconstruction reveals that adrenogonadal cells emerge from the lateral plate rather than the intermediate mesoderm. Surprisingly, we find that gonadal and adrenal fates have already diverged prior to Nr5a1 expression. Finally, lineage separation into gonadal and adrenal fates involves canonical versus non-canonical Wnt signaling and differential expression of Hox patterning genes. Thus, our study provides important insights into the molecular programs of adrenal and gonadal fate choice and will be a valuable resource for further research into adrenogonadal ontogenesis.

Truncation mutations in MYRF underlie primary angle closure glaucoma

Mutations in myelin regulatory factor (MYRF), a gene mapped to 11q12-q13.3, are responsible for autosomal dominant high hyperopia and seem to be associated with angle closure glaucoma, which is one of the leading causes of irreversible blindness worldwide. Whether there is a causal link from the MYRF mutations to the pathogenesis of primary angle-closure glaucoma (PACG) remains unclear at this time. Six truncation mutations, including five novel and one previously reported, in MYRF are identified in seven new probands with hyperopia, of whom all six adults have glaucoma, further confirming the association of MYRF mutations with PACG. Immunofluorescence microscopy demonstrates enriched expression of MYRF in the ciliary body and ganglion cell layer in humans and mice. Myrf^mut/+ mice have elevated IOP and fewer ganglion cells along with thinner retinal nerve fiber layer with ganglion cell layer than wild-type. Transcriptome sequencing of Myrf^mut/+ retinas shows downregulation of Dnmt3a, a gene previously associated with PACG. Co-immunoprecipitation demonstrates a physical association of DNMT3A with MYRF. DNA methylation sequencing identifies several glaucoma-related cell events in Myrf^mut/+ retinas. The interaction between MYRF and DNMT3A underlies MYRF-associated PACG and provides clues for pursuing further investigation into the pathogenesis of PACG and therapeutic target.

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.

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.

Case Report: De novo variant in myelin regulatory factor in a Chinese child with 46,XY disorder/difference of sex development, cardiac and urogenital anomalies, and short stature

The myelin regulatory factor (MYRF; MIM# 608329) gene was first identified as a critical transcription factor involved in oligodendrocyte differentiation and central nervous system myelination. With the recent development of exome sequencing, pathogenic variants of MYRF had been considered as the cause of cardiac-urogenital syndrome (CUGS), 46,XY and 46,XX disorders/differences of sex development (DSDs), and nanophthalmos. Herein, we described a 4-year-7-month-old "girl" with ventricular septal defect, atrial septal defect, patent ductus arteriosus, severe pulmonary hypertension, moderate-to-severe tricuspid regurgitation, enlarged coronary sinus, left superior vena cava, and right lung hypoplasia at birth. Later, the patient developed short stature and amblyopia. Further examination revealed a karyotype 46,XY and visible uterus, whereas the presence of gonads were not explored. Laparoscopy revealed dysplasia of testicular gonad. Whole-exome sequencing (WES) was performed and a de novo heterozygous mutation in MYRF was identified, known as c.2817G > A/p. W939* (NM_001127392.3). Therefore, this case report presented multiple clinical manifestations with syndromic symptoms of CUGS, 46,XY DSD, and ocular symptoms. These new data expanded the phenotype of the MYRF variant and may benefit to characterize the phenotypes caused by the variants of this gene.

CNS Hypomyelination Disrupts Axonal Conduction and Behavior in Larval Zebrafish

Myelination is essential for central nervous system (CNS) formation, health and function. As a model organism, larval zebrafish have been extensively employed to investigate the molecular and cellular basis of CNS myelination, because of their genetic tractability and suitability for non-invasive live cell imaging. However, it has not been assessed to what extent CNS myelination affects neural circuit function in zebrafish larvae, prohibiting the integration of molecular and cellular analyses of myelination with concomitant network maturation. To test whether larval zebrafish might serve as a suitable platform with which to study the effects of CNS myelination and its dysregulation on circuit function, we generated zebrafish myelin regulatory factor (myrf) mutants with CNS-specific hypomyelination and investigated how this affected their axonal conduction properties and behavior. We found that myrf mutant larvae exhibited increased latency to perform startle responses following defined acoustic stimuli. Furthermore, we found that hypomyelinated animals often selected an impaired response to acoustic stimuli, exhibiting a bias toward reorientation behavior instead of the stimulus-appropriate startle response. To begin to study how myelination affected the underlying circuitry, we established electrophysiological protocols to assess various conduction properties along single axons. We found that the hypomyelinated myrf mutants exhibited reduced action potential conduction velocity and an impaired ability to sustain high-frequency action potential firing. This study indicates that larval zebrafish can be used to bridge molecular and cellular investigation of CNS myelination with multiscale assessment of neural circuit function.SIGNIFICANCE STATEMENT Myelination of CNS axons is essential for their health and function, and it is now clear that myelination is a dynamic life-long process subject to modulation by neuronal activity. However, it remains unclear precisely how changes to myelination affects animal behavior and underlying action potential conduction along axons in intact neural circuits. In recent years, zebrafish have been employed to study cellular and molecular mechanisms of myelination, because of their relatively simple, optically transparent, experimentally tractable vertebrate nervous system. Here we find that changes to myelination alter the behavior of young zebrafish and action potential conduction along individual axons, providing a platform to integrate molecular, cellular, and circuit level analyses of myelination using this model.

Rare and de novo variants in 827 congenital diaphragmatic hernia probands implicate LONP1 as candidate risk gene

Congenital diaphragmatic hernia (CDH) is a severe congenital anomaly that is often accompanied by other anomalies. Although the role of genetics in the pathogenesis of CDH has been established, only a small number of disease-associated genes have been identified. To further investigate the genetics of CDH, we analyzed de novo coding variants in 827 proband-parent trios and confirmed an overall significant enrichment of damaging de novo variants, especially in constrained genes. We identified LONP1 (lon peptidase 1, mitochondrial) and ALYREF (Aly/REF export factor) as candidate CDH-associated genes on the basis of de novo variants at a false discovery rate below 0.05. We also performed ultra-rare variant association analyses in 748 affected individuals and 11,220 ancestry-matched population control individuals and identified LONP1 as a risk gene contributing to CDH through both de novo and ultra-rare inherited largely heterozygous variants clustered in the core of the domains and segregating with CDH in affected familial individuals. Approximately 3% of our CDH cohort who are heterozygous with ultra-rare predicted damaging variants in LONP1 have a range of clinical phenotypes, including other anomalies in some individuals and higher mortality and requirement for extracorporeal membrane oxygenation. Mice with lung epithelium-specific deletion of Lonp1 die immediately after birth, most likely because of the observed severe reduction of lung growth, a known contributor to the high mortality in humans. Our findings of both de novo and inherited rare variants in the same gene may have implications in the design and analysis for other genetic studies of congenital anomalies.

Crystal structure of the MyRF ICA domain with its upstream β-helical stalk reveals the molecular mechanisms underlying its trimerization and self-cleavage

Myelin gene regulatory factor (MyRF), a novel membrane transcription factor expressed on the endoplasmic reticulum membrane, functions as a trimer. The trimerization of MyRF is associated with a fragment between the DNA binding domain and transmembrane domain that shares homology with the triple-β-helix and intramolecular chaperone autocleavage (ICA) domain of phage tailspike proteins. The molecular details of these domains in eukaryotes have not been elucidated. Here, we present the crystal structure of the MyRF ICA domain with its upstream β-helical stalk, determined at 2.4Å resolution. The structure showed that its upstream β-helical stalk is different from the triple β-helix reported before. This is the first structure of the mammalian protein with a triple β-helix. Structure analysis demonstrated that the triple α-helical coiled-coil formed at the MyRF ICA domain C-terminal was the main driving force for the trimerization. Additionally, our findings showed that MyRF was cleaved via a highly conserved serine-lysine catalytic dyad mechanism and that cleavage would be activated only if the ICA domains were organized as trimers. In contrast to the viral ICA domain, almost no interaction was found between the MyRF ICA domain and its upstream neighboring β-helix of the stalk; thus, activation of self-cleavage may not be triggered by the upstream region of the ICA domain, contrary to the observations made in phages. These findings provided an important insight into the molecular mechanisms of MyRF trimerization and self-cleavage.

The embryonic ontogeny of the gonadal somatic cells in mice and monkeys

In the early fetal stage, the gonads are bipotent and only later become the ovary or testis, depending on the genetic sex. Despite many studies examining how sex determination occurs from biopotential gonads, the spatial and temporal organization of bipotential gonads and their progenitors is poorly understood. Here, using lineage tracing in mice, we find that the gonads originate from a T⁺ primitive streak through WT1⁺ posterior intermediate mesoderm and appear to share origins anteriorly with the adrenal glands and posteriorly with the metanephric mesenchyme. Comparative single-cell transcriptomic analyses in mouse and cynomolgus monkey embryos reveal the convergence of the lineage trajectory and genetic programs accompanying the specification of biopotential gonadal progenitor cells. This process involves sustained expression of epithelial genes and upregulation of mesenchymal genes, thereby conferring an epithelial-mesenchymal hybrid state. Our study provides key resources for understanding early gonadogenesis in mice and primates.

Functional mechanisms of MYRF DNA-binding domain mutations implicated in birth defects

Myrf is a pleiotropic membrane-bound transcription factor that plays critical roles in diverse organisms, including in oligodendrocyte differentiation, embryonic development, molting, and synaptic plasticity. Upon autolytic cleavage, the Myrf N-terminal fragment enters the nucleus as a homo-trimer and functions as a transcription factor. Homo-trimerization is essential for this function because it imparts DNA-binding specificity and affinity. Recent exome sequencing studies have implicated four de novo MYRF DNA-binding domain (DBD) mutations (F387S, Q403H, G435R, and L479V) in novel syndromic birth defects involving the diaphragm, heart, and the urogenital tract. It remains unknown whether and how these four mutations alter the transcription factor function of MYRF. Here, we studied them by introducing homologous mutations to the mouse Myrf protein. We found that the four DBD mutations abolish the transcriptional activity of the Myrf N-terminal fragment by interfering with its homo-trimerization ability by perturbing the DBD structure. Since the Myrf N-terminal fragment strictly functions as a homo-trimer, any loss-of-function mutation has the potential to act as a dominant negative. We observed that one copy of Myrf-F387S, Myrf-Q403H, or Myrf-L479V, but not Myrf-G435R, was tolerated by the Myrf N-terminal homo-trimer for structural and functional integrity. These data suggest that F387S, Q403H, and L479V cause birth defects by haploinsufficiency, while G435R does so via dominant negative functionality.

Functional mechanism and pathogenic potential of MYRF ICA domain mutations implicated in birth defects

Myrf is a membrane-bound transcription factor that plays a key role in various biological processes. The Intramolecular Chaperone Auto-processing (ICA) domain of Myrf forms a homo-trimer, which carries out the auto-cleavage of Myrf. The ICA homo-trimer-mediated auto-cleavage of Myrf is a prerequisite for its transcription factor function in the nucleus. Recent exome sequencing studies have implicated two MYRF ICA domain mutations (V679A and R695H) in a novel syndromic form of birth defects. It remains unknown whether and how the two mutations impact the transcription factor function of Myrf and, more importantly, how they are pathogenic for congenital anomalies. Here, we show that V679A and R695H cripple the ICA domain, blocking the auto-cleavage of Myrf. Consequently, Myrf-V679A and Myrf-R695H do not exhibit any transcriptional activity. Molecular modeling suggests that V679A and R695H abrogate the auto-cleavage function of the ICA homo-trimer by destabilizing its homo-trimeric assembly. We also found that the ICA homo-trimer can tolerate one copy of Myrf-V679A or Myrf-R695H for its auto-cleavage function, indicating that V679A and R695H are not dominant negatives. Thus, if V679A and R695H in a heterozygous state caused birth defects, it would be via haploinsufficiency of MYRF.