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.

Exploring unsolved cases of lissencephaly spectrum: integrating exome and genome sequencing for higher diagnostic yield.

Lissencephaly is a rare brain malformation characterized by abnormal neuronal migration during cortical development. In this study, we performed a comprehensive genetic analysis using next-generation sequencing in 12 unsolved Japanese lissencephaly patients, in whom PAFAH1B1, DCX, TUBA1A, and ARX variants were excluded using the Sanger method. Exome sequencing (ES) was conducted on these 12 patients, identifying pathogenic variants in CEP85L, DYNC1H1, LAMC3, and DCX in four patients. Next, we performed genome sequencing (GS) on eight unsolved patients, and structural variants in PAFAH1B1, including an inversion and microdeletions involving several exons, were detected in three patients. Notably, these microdeletions in PAFAH1B1 could not to be detected by copy number variation (CNV) detection tools based on the depth of coverage methods using ES data. The density of repeat sequences, including Alu sequences or segmental duplications, which increase the susceptibility to structural variations, is very high in some lissencephaly spectrum genes (PAFAH1B1, TUBA1A, DYNC1H1). These missing CNVs were due to the limitations of detecting repeat sequences in ES-based CNV detection tools. Our study suggests that a combined approach integrating ES with GS can contribute to a higher diagnostic yield and a better understanding of the genetic landscape of the lissencephaly spectrum.

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.

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.

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.

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.

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.

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.

Missense and truncating variants in CHD5 in a dominant neurodevelopmental disorder with intellectual disability, behavioral disturbances, and epilepsy.

Located in the critical 1p36 microdeletion region, the chromodomain helicase DNA-binding protein 5 (CHD5) gene encodes a subunit of the nucleosome remodeling and deacetylation (NuRD) complex required for neuronal development. Pathogenic variants in six of nine chromodomain (CHD) genes cause autosomal dominant neurodevelopmental disorders, while CHD5-related disorders are still unknown. Thanks to GeneMatcher and international collaborations, we assembled a cohort of 16 unrelated individuals harboring heterozygous CHD5 variants, all identified by exome sequencing. Twelve patients had de novo CHD5 variants, including ten missense and two splice site variants. Three familial cases had nonsense or missense variants segregating with speech delay, learning disabilities, and/or craniosynostosis. One patient carried a frameshift variant of unknown inheritance due to unavailability of the father. The most common clinical features included language deficits (81%), behavioral symptoms (69%), intellectual disability (64%), epilepsy (62%), and motor delay (56%). Epilepsy types were variable, with West syndrome observed in three patients, generalized tonic-clonic seizures in two, and other subtypes observed in one individual each. Our findings suggest that, in line with other CHD-related disorders, heterozygous CHD5 variants are associated with a variable neurodevelopmental syndrome that includes intellectual disability with speech delay, epilepsy, and behavioral problems as main features.

Biallelic Variants in CNPY3, Encoding an Endoplasmic Reticulum Chaperone, Cause Early-Onset Epileptic Encephalopathy.

Early-onset epileptic encephalopathies, including West syndrome (WS), are a group of neurological disorders characterized by developmental impairments and intractable seizures from early infancy. We have now identified biallelic CNPY3 variants in three individuals with WS; these include compound-heterozygous missense and frameshift variants in a family with two affected siblings (individuals 1 and 2) and a homozygous splicing variant in a consanguineous family (individual 3). All three individuals showed hippocampal malrotation. In individuals 1 and 2, electroencephalography (EEG) revealed characteristic fast waves and diffuse sharp- and slow-wave complexes. The fast waves were clinically associated with seizures. CNPY3 encodes a co-chaperone in the endoplasmic reticulum and regulates the subcellular distribution and responses of multiple Toll-like receptors. The amount of CNPY3 in lymphoblastoid cells derived from individuals 1 and 2 was severely lower than that in control cells. Cnpy3-knockout mice exhibited spastic or dystonic features under resting conditions and hyperactivity and anxiolytic behavior during the open field test. Also, their resting EEG showed enhanced activity in the fast beta frequency band (20-35 Hz), which could mimic the fast waves in individuals 1 and 2. These data suggest that CNPY3 and Cnpy3 perform essential roles in brain function in addition to known Toll-like receptor-dependent immune responses.

A boy with biallelic frameshift variants in TTC5 and brain malformation resembling tubulinopathies.

Brain malformations have heterogeneous genetic backgrounds. Tubulinopathies are a wide range of brain malformations caused by variants in tubulin and microtubules-associated genes. Recently biallelic variants in TTC5, also known as stress responsive activator of p300, have been reported in 11 patients from seven families with developmental delay, intellectual disability, and brain malformations. Here, we report compound heterozygous frameshift variants in TTC5 in a Japanese boy who showed severe psychomotor developmental delay and pseudobulbar palsy with growth failure. Brain magnetic resonance imaging showed a simplified gyral pattern and undetectable anterior limb of the internal capsule, suggesting tubulinopathies. Immunoblotting using lymphoblastoid cells derived from the patient showed undetectable TTC5 protein. Ttc5 silencing by RNA interference in Neuro2a cells reduced Tubulin ß3 protein level and caused abnormal cell cycle. Our report suggests a possible link between TTC5-related brain malformation and tubulinopathies.

Identification of two novel de novo TUBB variants in cases with brain malformations: case reports and literature review.

Heterozygous variants in TUBB encoding one of ß-tubulin isotypes are known to cause two overlapping developmental brain disorders, complex cortical dysplasia with other brain malformations (CDCBM) and congenital symmetric circumferential skin creases (CSCSC). To date, six cases of CSCSC and eight cases of CDCBM caused by nine heterozygous variants have been reported. Here we report two cases with novel de novo missense TUBB variants (NM_178014.4:c.863A>G, p.(Glu288Gly) and c.869C>T, p.(Thr290Ile)). Case 1 presented brain malformations consistent with tubulinopathies including abnormalities in cortex, basal ganglia, corpus callosum, brain stem, and cerebellum along with other systemic features such as coloboma, facial dysmorphisms, vesicoureteral reflux, hypoplastic kidney, and cutis laxa-like mild skin loosening. Another case presented abnormalities of the corpus callosum, brain stem, and cerebellum along with facial dysmorphisms. We reviewed previous literature and suggest the diversity of clinical findings of TUBB-related disorders.

Global developmental delay, systemic dysmorphism and epilepsy in a patient with a de novo U2AF2 variant.

U2 small nuclear RNA auxiliary factor 2 (U2AF2) is an essential pre-mRNA splicing factor in an early step of splicing. Alternative splicing plays an important role in neuronal development, and disorders of RNA processing steps are implicated in neurological disorders. Recently, the large trio whole-exome sequencing study reported U2AF2 as a novel gene significantly associated with developmental disorders: however, the clinical details of patients with U2AF2 variants were not available. Here, we report an individual with a de novo U2AF2 variant (c.445C>T, p.(Arg149Trp)) using trio-based whole-exome sequencing. This residue was positioned in the RNA recognition motif 1 which recognizes a polypyrimidine-tract splice site signal. The patient showed global developmental delay, intellectual disability, epilepsy, short stature, microcephaly, facial dysmorphism, intermittent exotropia, bilateral ptosis, muscle hypotonia and thin corpus callosum, indicating that U2AF2-related disorder could include systemic dysmorphisms, epilepsy and brain malformation along with global developmental delay.

Comprehensive genetic analysis confers high diagnostic yield in 16 Japanese patients with corpus callosum anomalies.

Corpus callosum anomalies (CCA) is a common congenital brain anomaly with various etiologies. Although one of the most important etiologies is genetic factors, the genetic background of CCA is heterogenous and diverse types of variants are likely to be causative. In this study, we analyzed 16 Japanese patients with corpus callosum anomalies to delineate clinical features and the genetic background of CCAs. We observed the common phenotypes accompanied by CCAs: intellectual disability (100%), motor developmental delay (93.8%), seizures (60%), and facial dysmorphisms (50%). Brain magnetic resonance imaging showed colpocephaly (enlarged posterior horn of the lateral ventricles, 84.6%) and enlarged supracerebellar cistern (41.7%). Whole exome sequencing revealed genetic alterations in 9 of the 16 patients (56.3%), including 8 de novo alterations (2 copy number variants and variants in ARID1B, CDK8, HIVEP2, and TCF4) and a recessive variant of TBCK. De novo ARID1B variants were identified in three unrelated individuals, suggesting that ARID1B variants are major genetic causes of CCAs. A de novo TCF4 variant and somatic mosaic deletion at 18q21.31-qter encompassing TCF4 suggest an association of TCF4 abnormalities with CCAs. This study, which analyzes CCA patients usung whole exome sequencing, demonstrates that comprehensive genetic analysis would be useful for investigating various causal variants of CCAs.

Genetic and phenotypic analysis of 101 patients with developmental delay or intellectual disability using whole-exome sequencing.

Whole-exome sequencing (WES) enables identification of pathogenic variants, including copy number variants (CNVs). In this study, we performed WES in 101 Japanese patients with unexplained developmental delay (DD) or intellectual disability (ID) (63 males and 38 females), 98 of them with trio-WES. Pathogenic variants were identified in 54 cases (53.5%), including four cases with pathogenic CNVs. In one case, a pathogenic variant was identified by reanalysis of exome data; and in two cases, two molecular diagnoses were identified. Among 58 pathogenic variants, 49 variants occurred de novo in 48 patients, including two somatic variants. The accompanying autism spectrum disorder and external ear anomalies were associated with detection of pathogenic variants with odds ratios of 11.88 (95% confidence interval [CI] 2.52-56.00) and 3.46 (95% CI 1.23-9.73), respectively. These findings revealed the importance of reanalysis of WES data and detection of CNVs and somatic variants in increasing the diagnostic yield for unexplained DD/ID. In addition, genetic testing is recommended when patients suffer from the autism spectrum disorder or external ear anomalies, which potentially suggests the involvement of genetic factors associated with gene expression regulation.