Two novel cases of biallelic SMPD4 variants with brain structural abnormalities.

Sphingomyelin phosphodiesterase 4 (SMPD4) encodes a member of the Mg2+-dependent, neutral sphingomyelinase family that catalyzes the hydrolysis of the phosphodiester bond of sphingomyelin to form phosphorylcholine and ceramide. Recent studies have revealed that biallelic loss-of-function variants of SMPD4 cause syndromic neurodevelopmental disorders characterized by microcephaly, congenital arthrogryposis, and structural brain anomalies. In this study, three novel loss-of-function SMPD4 variants were identified using exome sequencing (ES) in two independent patients with developmental delays, microcephaly, seizures, and brain structural abnormalities. Patient 1 had a homozygous c.740_741del, p.(Val247Glufs*21) variant and showed profound intellectual disability, hepatomegaly, a simplified gyral pattern, and a thin corpus callosum without congenital dysmorphic features. Patient 2 had a compound heterozygous nonsense c.2124_2125del, p.(Phe709*) variant and splice site c.1188+2dup variant. RNA analysis revealed that the c.1188+2dup variant caused exon 13 skipping, leading to a frameshift (p.Ala406Ser*6). In vitro transcription analysis using minigene system suggested that mRNA transcribed from mutant allele may be degraded by nonsense-mediated mRNA decay system. He exhibited diverse manifestations, including growth defects, muscle hypotonia, respiratory distress, arthrogryposis, insulin-dependent diabetes mellitus, sensorineural hearing loss, facial dysmorphism, and various brain abnormalities, including cerebral atrophy, hypomyelination, and cerebellar hypoplasia. Here, we review previous literatures and discuss the phenotypic diversity of SMPD4-related disorders.

A novel heterozygous TMEM63A variant in a familial case with early onset nystagmus, severe hypomyelination, and a favorable adult prognosis.

Heterozygous transmembrane protein 63A (TMEM63A) variants cause transient infantile hypomyelinating leukodystrophy-19, which features remarkable natural resolution of clinical and imaging findings during childhood. Previous reports have mainly described de novo variants lacking detailed familial cases. Herein, we describe the clinical course of familial cases with a TMEM63A variant. A 5-month-old girl presented with nystagmus, global hypotonia, and difficulty swallowing since birth. Brain magnetic resonance imaging at 1.5 and 5 months revealed diffuse hypomyelination. Her mother, maternal aunt, and grandfather had nystagmus and motor developmental delays in infancy, which resolved spontaneously during childhood. Compared with these cases, the proband's motor developmental delay was profound, and she was the only one with feeding difficulties, necessitating nasogastric tube feeding. Genetic testing revealed a heterozygous TMEM63A variant (NM_014698.3:c.1658G>A, p.(Gly553Asp)) in the proband and her family. This is the first three-generation familial report of a TMEM63A variant that provides insight into its history and heterogeneity.

RNA sequencing and target long-read sequencing reveal an intronic transposon insertion causing aberrant splicing.

More than half of cases with suspected genetic disorders remain unsolved by genetic analysis using short-read sequencing such as exome sequencing (ES) and genome sequencing (GS). RNA sequencing (RNA-seq) and long-read sequencing (LRS) are useful for interpretation of candidate variants and detection of structural variants containing repeat sequences, respectively. Recently, adaptive sampling on nanopore sequencers enables target LRS more easily. Here, we present a Japanese girl with premature chromatid separation (PCS)/mosaic variegated aneuploidy (MVA) syndrome. ES detected a known pathogenic maternal heterozygous variant (c.1402-5A>G) in intron 10 of BUB1B (NM_001211.6), a known responsive gene for PCS/MVA syndrome with autosomal recessive inheritance. Minigene splicing assay revealed that almost all transcripts from the c.1402-5G allele have mis-splicing with 4-bp insertion. GS could not detect another pathogenic variant, while RNA-seq revealed abnormal reads in intron 2. To extensively explore variants in intron 2, we performed adaptive sampling and identified a paternal 3.0 kb insertion. Consensus sequence of 16 reads spanning the insertion showed that the insertion consists of Alu and SVA elements. Realignment of RNA-seq reads to the new reference sequence containing the insertion revealed that 16 reads have 5' splice site within the insertion and 3' splice site at exon 3, demonstrating causal relationship between the insertion and aberrant splicing. In addition, immunoblotting showed severely diminished BUB1B protein level in patient derived cells. These data suggest that detection of transcriptomic abnormalities by RNA-seq can be a clue for identifying pathogenic variants, and determination of insert sequences is one of merits of LRS.

Novel compound heterozygous ATP1A2 variants in a patient with fetal akinesia/hypokinesia sequence.

ATP1A2 encodes a subunit of sodium/potassium-transporting adenosine triphosphatase (Na+ /K+ -ATPase). Heterozygous pathogenic variants of ATP1A2 cause familial hemiplegic migraine, alternating hemiplegia of childhood, and developmental and epileptic encephalopathy. Biallelic loss-of-function variants in ATP1A2 lead to fetal akinesia, respiratory insufficiency, microcephaly, polymicrogyria, and dysmorphic facies, resulting in fetal death. Here, we describe a patient with compound heterozygous ATP1A2 variants consisting of missense and nonsense variants. He survived after birth with brain malformations and the fetal akinesia/hypokinesia sequence. We report a novel type of compound heterozygous variant that might extend the disease spectrum of ATP1A2.

Long-term clinical observation of patients with heterozygous KIF1A variants.

KIF1A-related disorders (KRDs) encompass recessive and dominant variants with wide clinical variability. Recent genetic investigations have expanded the clinical phenotypes of heterozygous KIF1A variants. However, there have been a few long-term observational studies of patients with heterozygous KIF1A variants. A retrospective chart review of consecutive patients diagnosed with spastic paraplegia at Miyagi Children's Hospital from 2016 to 2020 identified six patients with heterozygous KIF1A variants. To understand the long-term changes in clinical symptoms, we examined these patients in terms of their characteristics, clinical symptoms, results of electrophysiological and neuroimaging studies, and genetic testing. The median follow-up period was 30 years (4-44 years). This long-term observational study showed that early developmental delay and equinus gait, or unsteady gait, are the first signs of disease onset, appearing with the commencement of independent walking. In addition, later age-related progression was observed in spastic paraplegia, and the appearance of axonal neuropathy and reduced visual acuity were characteristic features of the late disease phenotype. Brain imaging showed age-related progression of cerebellar atrophy and the appearance of hyperintensity of optic radiation on T2WI and FLAIR imaging. Long-term follow-up revealed a pattern of steady progression and a variety of clinical symptoms, including spastic paraplegia, peripheral neuropathy, reduced visual acuity, and some degree of cerebellar ataxia. Clinical variability between patients was observed to some extent, and therefore, further studies are required to determine the phenotype-genotype correlation.

Sotos syndrome with marked overgrowth in three Japanese patients with heterozygous likely pathogenic NSD1 variants: case reports with review of literature.

We report three Japanese patients with Sotos syndrome accompanied by marked overgrowth, i.e., a 2 8/12-year-old boy with a height of 105.2 cm (+4.4 SD) (patient 1), the mother of patient 1 with a height of 180.8 cm (+4.1 SD) (patient 2), and a 12 10/12-year-old girl with a height of 189.4 cm (+6.3 SD) (patient 3). In addition to the marked overgrowth (tall stature), patients 1-3 exhibited Sotos syndrome-compatible macrocephaly and characteristic features, whereas intellectual and developmental disabilities remained at a borderline level in patient 1 and were apparently absent from patients 2 and 3. Thus, whole exome sequencing was performed to confirm the diagnosis, revealing a likely pathogenic c.6356A>G:p.(Asp2119Gly) variant in NSD1 of patients 1 and 2, and a likely pathogenic c.6599dupT:p.(Ser2201Valfs*4) variant in NSD1 of patient 3 (NM_022455.5). The results, in conjunction with the previously reported data in nine patients with marked overgrowth (≥4.0 SD), imply that several patients with Sotos syndrome have extreme tall stature even in adulthood. Thus, it is recommended to examine NSD1 in patients with marked overgrowth as the salient feature.

Identification of pathogenic deep intronic variant and exonic LINE-1 insertion in a patient with Meckel syndrome.

Biallelic CC2D2A variants are associated with a wide range of neurodevelopmental disorders including Meckel syndrome. Here we report a Japanese girl with Meckel syndrome harboring a pathogenic deep intronic variant (NM_001378615.1:c.1149+3569A>G) and an exonic LINE-1 insertion, which was predicted to cause aberrant splicing by SpliceAI and was detected by TEMP2 program, respectively. RNA analysis using urine-derived cells (UDCs) showed retention of 149-bp intronic sequences, leading to frameshift. Immunoblotting showed marked reduction of CC2D2A protein in the patient. Our report demonstrated that utilization of transposon detection tool and functional analysis using UDCs will increase diagnostic yield of genome sequencing.

De Novo Truncating Variants in the Last Exon of SEMA6B Cause Progressive Myoclonic Epilepsy.

De novo variants (DNVs) cause many genetic diseases. When DNVs are examined in the whole coding regions of genes in next-generation sequencing analyses, pathogenic DNVs often cluster in a specific region. One such region is the last exon and the last 50 bp of the penultimate exon, where truncating DNVs cause escape from nonsense-mediated mRNA decay [NMD(-) region]. Such variants can have dominant-negative or gain-of-function effects. Here, we first developed a resource of rates of truncating DNVs in NMD(-) regions under the null model of DNVs. Utilizing this resource, we performed enrichment analysis of truncating DNVs in NMD(-) regions in 346 developmental and epileptic encephalopathy (DEE) trios. We observed statistically significant enrichment of truncating DNVs in semaphorin 6B (SEMA6B) (p value: 2.8 × 10-8; exome-wide threshold: 2.5 × 10-6). The initial analysis of the 346 individuals and additional screening of 1,406 and 4,293 independent individuals affected by DEE and developmental disorders collectively identified four truncating DNVs in the SEMA6B NMD(-) region in five individuals who came from unrelated families (p value: 1.9 × 10-13) and consistently showed progressive myoclonic epilepsy. RNA analysis of lymphoblastoid cells established from an affected individual showed that the mutant allele escaped NMD, indicating stable production of the truncated protein. Importantly, heterozygous truncating variants in the NMD(+) region of SEMA6B are observed in general populations, and SEMA6B is most likely loss-of-function tolerant. Zebrafish expressing truncating variants in the NMD(-) region of SEMA6B orthologs displayed defective development of brain neurons and enhanced pentylenetetrazole-induced seizure behavior. In summary, we show that truncating DNVs in the final exon of SEMA6B cause progressive myoclonic epilepsy.

De novo CLCN3 variants affecting Gly327 cause severe neurodevelopmental syndrome with brain structural abnormalities.

A recent study revealed that monoallelic missense or biallelic loss-of-function variants in the chloride voltage-gated channel 3 (CLCN3) cause neurodevelopmental disorders resulting in brain abnormalities. Functional studies suggested that some missense variants had varying gain-of-function effects on channel activity. Meanwhile, two patients with homozygous frameshift variants showed severe neuropsychiatric disorders and a range of brain structural abnormalities. Here we describe two patients with de novo CLCN3 variants affecting the same amino acid, Gly327 (p.(Gly327Ser) and p.(Gly327Asp)). They showed severe neurological phenotypes including global developmental delay, intellectual disability, hypotonia, failure to thrive, and various brain abnormalities. They also presented with characteristic brain and ophthalmological abnormalities, hippocampal and retinal degradation, which were observed in patients harboring homozygous loss-of-function variants. These findings were also observed in CLCN3-deficient mice, indicating that the monoallelic missense variant may also have a dominant negative effect. This study will expand the phenotypic spectrum of CLCN3-related disorders.

A novel de novo KCNB1 variant altering channel characteristics in a patient with periventricular heterotopia, abnormal corpus callosum, and mild seizure outcome.

KCNB1 encodes the α-subunit of Kv2.1, the main contributor to neuronal delayed rectifier potassium currents. The subunit consists of six transmembrane α helices (S1-S6), comprising the voltage-sensing domain (S1-S4) and the pore domain (S5-P-S6). Heterozygous KCNB1 pathogenic variants are associated with developmental and epileptic encephalopathy. Here we report an individual who shows the milder phenotype compared to the previously reported cases, including delayed language development, mild intellectual disability, attention deficit hyperactivity disorder, late-onset epilepsy responsive to an antiepileptic drug, elevation of serum creatine kinase, and peripheral axonal neuropathy. On the other hand, his brain MRI showed characteristic findings including periventricular heterotopia, polymicrogyria, and abnormal corpus callosum. Exome sequencing identified a novel de novo KCNB1 variant c.574G>A, p.(Ala192Thr) located in the S1 segment of the voltage-sensing domain. Functional analysis using the whole-cell patch-clamp technique in Neuro2a cells showed that the Ala192Thr mutant reduces both activation and inactivation of the channel at membrane voltages in the range of -50 to -30 mV. Our case could expand the phenotypic spectrum of patients with KCNB1 variants, and suggested that variants located in the S1 segment might be associated with a milder outcome of seizures.

A deep intronic TCTN2 variant activating a cryptic exon predicted by SpliceRover in a patient with Joubert syndrome.

The recent introduction of genome sequencing in genetic analysis has led to the identification of pathogenic variants located in deep introns. Recently, several new tools have emerged to predict the impact of variants on splicing. Here, we present a Japanese boy of Joubert syndrome with biallelic TCTN2 variants. Exome sequencing identified only a heterozygous maternal nonsense TCTN2 variant (NM_024809.5:c.916C >T, p.(Gln306Ter)). Subsequent genome sequencing identified a deep intronic variant (c.1033+423G>A) inherited from his father. The machine learning algorithms SpliceAI, Squirls, and Pangolin were unable to predict alterations in splicing by the c.1033+423G>A variant. SpliceRover, a tool for splice site prediction using FASTA sequence, was able to detect a cryptic exon which was 85-bp away from the variant and within the inverted Alu sequence while SpliceRover scores for these splice sites showed slight increase (donor) or decrease (acceptor) between the reference and mutant sequences. RNA sequencing and RT-PCR using urinary cells confirmed inclusion of the cryptic exon. The patient showed major symptoms of TCTN2-related disorders such as developmental delay, dysmorphic facial features and polydactyly. He also showed uncommon features such as retinal dystrophy, exotropia, abnormal pattern of respiration, and periventricular heterotopia, confirming these as one of features of TCTN2-related disorders. Our study highlights usefulness of genome sequencing and RNA sequencing using urinary cells for molecular diagnosis of genetic disorders and suggests that database of cryptic splice sites predicted in introns by SpliceRover using the reference sequences can be helpful in extracting candidate variants from large numbers of intronic variants in genome sequencing.

The Structural Abnormalities Are Deeply Involved in the Cause of RPGRIP1-Related Retinal Dystrophy in Japanese Patients.

Leber congenital amaurosis (LCA) is the most severe form of inherited retinal dystrophy. RPGRIP1-related LCA accounts for 5-6% of LCA. We performed whole-exome sequencing and whole-genome sequencing (WGS) on 29 patients with clinically suspected LCA and examined ophthalmic findings in patients with biallelic pathogenic variants of RPGRIP1. In addition to five previously reported cases, we identified five cases from four families with compound heterozygous RPGRIP1 variants using WGS. Five patients had null variants comprising frameshift variants, an Alu insertion, and microdeletions. A previously reported 1339 bp deletion involving exon 18 was found in four cases, and the deletion was relatively prevalent in the Japanese population (allele frequency: 0.002). Microdeletions involving exon 1 were detected in four cases. In patients with RPGRIP1 variants, visual acuity remained low, ranging from light perception to 0.2, and showed no correlation with age. In optical coherence tomography images, the ellipsoid zone (EZ) length decreased with age in all but one case of unimpaired EZ. The retinal structure was relatively preserved in all cases; however, there were cases with great differences in visual function compared to their siblings and a 56-year-old patient who still had a faint EZ line. Structural abnormalities may be important genetic causes of RPGRIP1-related retinal dystrophy in Japanese patients, and WGS was useful for detecting them.

An intronic GNAO1 variant leading to in-frame insertion cause movement disorder controlled by deep brain stimulation.

GNAO1 variants are associated with a wide range of neurodevelopmental disorders including epileptic encephalopathies and movement disorders. It has been reported that some GNAO1 variants are associated with movement disorders, and the 207-246 amino acid region was proposed as a mutational hotspot. Here, we report an intronic variant (NM_020988.3:c.724-8G>A) in GNAO1 in a Japanese girl who showed mild developmental delay and movement disorders including dystonia and myoclonus. Her movement disorders were improved by deep brain stimulation treatment as previously reported. This variant has been recurrently reported in four patients and was transmitted from her mother who possessed the variant as low-prevalent mosaicism. Using RNA extracted from lymphoblastoid cells derived from the patient, we demonstrated that the variant caused abnormal splicing of in-frame 6-bp intronic retention, leading to 2 amino acid insertion (p.Thr241_Asn242insProGln). Immunoblotting and immunostaining using WT and mutant GNAO1 vectors showed no significant differences in protein expression levels, but the cellular localization pattern of this mutant was partially shifted to the cytoplasm whereas WT was exclusively localized in the cellular membrane. Our report first clarified abnormal splicing and resulting mutant protein caused by the c.724-8G>A variant.

Elucidation of pathological mechanism caused by human disease mutation in CaMKIIß.

Recently, we have identified CaMKIIα and CaMKIIß mutations in patients with neurodevelopmental disorders by whole exome sequencing study. Most CaMKII mutants have increased phosphorylation of Thr286/287, which induces autonomous activity of CaMKII, using cell culture experiments. In this study, we explored the pathological mechanism of motor dysfunction observed exclusively in a patient with Pro213Leu mutation in CaMKIIß using a mouse model of the human disease. The homozygous CaMKIIß Pro213Leu knockin mice showed age-dependent motor dysfunction and growth failure from 2 weeks after birth. In the cerebellum, the mutation did not alter the mRNA transcript level, but the CaMKIIß protein level was dramatically decreased. Furthermore, in contrast to previous result from cell culture, Thr287 phosphorylation of CaMKIIß was also reduced. CaMKIIß Pro213Leu knockin mice showed similar motor dysfunction as CaMKIIß knockout mice, newly providing evidence for a loss of function rather than a gain of function. Our disease model mouse showed similar phenotypes of the patient, except for epileptic seizures. We clearly demonstrated that the pathological mechanism is a reduction of mutant CaMKIIß in the brain, and the physiological aspects of mutation were greatly different between in vivo and cell culture.

ACAN biallelic variants in a girl with severe idiopathic short stature.

Although ACAN heterozygous loss-of-function variants often cause idiopathic short stature (ISS) phenotype, there is no report describing ISS phenotype caused by ACAN biallelic loss-of-function variants. We encountered a 4 1/12-year-old Japanese girl with a height of 80.4 cm (-5.2 SD), a weight of 11.4 kg (-1.9 SD), a head circumference of 48.7 cm (-0.6 SD), and an arm span/height ratio of 1.0 (+1.1 SD). Endocrine studies and bone survey showed no abnormal findings. Whole exome sequencing revealed biallelic rare variants in ACAN, i.e., NM_013227.4:c.4214delC:p.(Pro1405Leufs*3) derived from her father and paternal grandfather with short stature (-2.9 and -2.0 SD, respectively) and NM_013227.4:c.7124 A>G:p.(Gln2375Arg) inherited from her mother and maternal grandmother with short stature (-2.1 and -3.0 SD, respectively). The frameshift variant underwent nonsense mediated mRNA decay, and the missense variant was assessed to have high pathogenicity. The results imply for the first time that ACAN biallelic loss-of-function variants can cause severe ISS phenotype.

Genome sequencing and RNA sequencing of urinary cells reveal an intronic FBN1 variant causing aberrant splicing.

Exome sequencing and panel testing have improved diagnostic yield in genetic analysis by comprehensively detecting pathogenic variants in exonic regions. However, it is important to identify non-exonic pathogenic variants to further improve diagnostic yield. Here, we present a female proband and her father who is diagnosed with Marfan syndrome, a systemic connective tissue disorder caused by pathogenic variants in FBN1. There are also two affected individuals in the siblings of the father, indicating the genetic basis in this family. However, panel testing performed by two institutions reported no causal variants. To further explore the genetic basis of the family, we performed genome sequencing of the proband and RNA sequencing of urinary cells derived from urine samples of the proband and her father because FBN1 is strongly expressed in urinary cells though it is poorly expressed in peripheral blood mononuclear cells. Genome sequencing identified a rare intronic variant (c.5789-15G>A) in intron 47 of FBN1 (NM_000138.4), which was transmitted from her father. RNA sequencing revealed allelic imbalance (monoallelic expression) of FBN1, retention of intron 47, and fewer aberrant transcripts utilizing new acceptor sites within exon 48, which were confirmed by RT-PCR. These results highlighted urinary cells as clinically accessible tissues for RNA sequencing if disease-causing genes are not sufficiently expressed in the blood, and the usefulness of multi-omics analysis for molecular diagnosis of genetic disorders.

Retrotransposition disrupting EBP in a girl and her mother with X-linked dominant chondrodysplasia punctata.

X-linked dominant chondrodysplasia punctata (CDPX2) is a rare congenital disorder caused by pathogenic variants in EBP on Xp11.23. We encountered a girl and her mother with CDPX2-compatible phenotypes including punctiform calcification in the neonatal period of the girl, and asymmetric limb shortening and ichthyosis following the Blaschko lines in both subjects. Although Sanger direct sequencing failed to reveal a disease-causing variant in EBP, whole genome sequencing (WGS) followed by Manta analysis identified a ~ 4.5 kb insertion at EBP exon 2 of both subjects. The insertion was associated with the hallmarks of retrotransposition such as an antisense poly(A) tail, a target site duplication, and a consensus endonuclease cleavage site, and the inserted sequence harbored full-length SVA_F1 element with 5'- and 3'-transductions containing the Alu sequence. The results imply the relevance of retrotransposition to the human genetic diseases and the usefulness of WGS in the identification of retrotransposition.

Generation of Flag/DYKDDDDK Epitope Tag Knock-In Mice Using i-GONAD Enables Detection of Endogenous CaMKIIα and ß Proteins.

Specific antibodies are necessary for cellular and tissue expression, biochemical, and functional analyses of protein complexes. However, generating a specific antibody is often time-consuming and effort-intensive. The epitope tagging of an endogenous protein at an appropriate position can overcome this problem. Here, we investigated epitope tag position using AlphaFold2 protein structure prediction and developed Flag/DYKDDDDK tag knock-in CaMKIIα and CaMKIIß mice by combining CRISPR-Cas9 genome editing with electroporation (i-GONAD). With i-GONAD, it is possible to insert a small fragment of up to 200 bp into the genome of the target gene, enabling efficient and convenient tagging of a small epitope. Experiments with commercially available anti-Flag antibodies could readily detect endogenous CaMKIIα and ß proteins by Western blotting, immunoprecipitation, and immunohistochemistry. Our data demonstrated that the generation of Flag/DYKDDDDK tag knock-in mice by i-GONAD is a useful and convenient choice, especially if specific antibodies are unavailable.

De novo variants in CELF2 that disrupt the nuclear localization signal cause developmental and epileptic encephalopathy.

We report heterozygous CELF2 (NM_006561.3) variants in five unrelated individuals: Individuals 1-4 exhibited developmental and epileptic encephalopathy (DEE) and Individual 5 had intellectual disability and autistic features. CELF2 encodes a nucleocytoplasmic shuttling RNA-binding protein that has multiple roles in RNA processing and is involved in the embryonic development of the central nervous system and heart. Whole-exome sequencing identified the following CELF2 variants: two missense variants [c.1558C>T:p.(Pro520Ser) in unrelated Individuals 1 and 2, and c.1516C>G:p.(Arg506Gly) in Individual 3], one frameshift variant in Individual 4 that removed the last amino acid of CELF2 c.1562dup:p.(Tyr521Ter), possibly resulting in escape from nonsense-mediated mRNA decay (NMD), and one canonical splice site variant, c.272-1G>C in Individual 5, also probably leading to NMD. The identified variants in Individuals 1, 2, 4, and 5 were de novo, while the variant in Individual 3 was inherited from her mosaic mother. Notably, all identified variants, except for c.272-1G>C, were clustered within 20 amino acid residues of the C-terminus, which might be a nuclear localization signal. We demonstrated the extranuclear mislocalization of mutant CELF2 protein in cells transfected with mutant CELF2 complementary DNA plasmids. Our findings indicate that CELF2 variants that disrupt its nuclear localization are associated with DEE.

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

Disorders of sex development (DSDs) are defined as congenital conditions in which chromosomal, gonadal or anatomical sex is atypical. In many DSD cases, genetic causes remain to be elucidated. Here, we performed a case-control exome sequencing study comparing gene-based burdens of rare damaging variants between 26 DSD cases and 2625 controls. We found exome-wide significant enrichment of rare heterozygous truncating variants in the MYRF gene encoding myelin regulatory factor, a transcription factor essential for oligodendrocyte development. All three variants occurred de novo. We identified an additional 46,XY DSD case of a de novo damaging missense variant in an independent cohort. The clinical symptoms included hypoplasia of Müllerian derivatives and ovaries in 46,XX DSD patients, defective development of Sertoli and Leydig cells in 46,XY DSD patients and congenital diaphragmatic hernia in one 46,XY DSD patient. As all of these cells and tissues are or partly consist of coelomic epithelium (CE)-derived cells (CEDC) and CEDC developed from CE via proliferaiton and migration, MYRF might be related to these processes. Consistent with this hypothesis, single-cell RNA sequencing of foetal gonads revealed high expression of MYRF in CE and CEDC. Reanalysis of public chromatin immunoprecipitation sequencing data for rat Myrf showed that genes regulating proliferation and migration were enriched among putative target genes of Myrf. These results suggested that MYRF is a novel causative gene of 46,XY and 46,XX DSD and MYRF is a transcription factor regulating CD and/or CEDC proliferation and migration, which is essential for development of multiple organs.