Readily-accessible and standardised capture of genotypic variation has revolutionised our understanding of the genetic contribution to disease. Unfortunately, the corresponding systematic capture of patient phenotypic variation needed to fully interpret the impact of genetic variation has lagged far behind. Exploiting deep and systematic phenotyping of a cohort of 197 patients presenting with heterogeneous developmental disorders and whose genomes harbour de novo CNVs, we systematically applied a range of commonly-used functional genomics approaches to identify the underlying molecular perturbations and their phenotypic impact. Grouping patients into 408 non-exclusive patient-phenotype groups, we identified a functional association amongst the genes disrupted in 209 (51%) groups. We find evidence for a significant number of molecular interactions amongst the association-contributing genes, including a single highly-interconnected network disrupted in 20% of patients with intellectual disability, and show using microcephaly how these molecular networks can be used as baits to identify additional members whose genes are variant in other patients with the same phenotype. Exploiting the systematic phenotyping of this cohort, we observe phenotypic concordance amongst patients whose variant genes contribute to the same functional association but note that (i) this relationship shows significant variation across the different approaches used to infer a commonly perturbed molecular pathway, and (ii) that the phenotypic similarities detected amongst patients who share the same inferred pathway perturbation result from these patients sharing many distinct phenotypes, rather than sharing a more specific phenotype, inferring that these pathways are best characterized by their pleiotropic effects.

Mutations in more than a hundred genes have been reported to cause X-linked recessive intellectual disability (ID) mainly in males. In contrast, the number of identified X-linked genes in which de novo mutations specifically cause ID in females is limited. Here, we report 17 females with de novo loss-of-function mutations in USP9X, encoding a highly conserved deubiquitinating enzyme. The females in our study have a specific phenotype that includes ID/developmental delay (DD), characteristic facial features, short stature, and distinct congenital malformations comprising choanal atresia, anal abnormalities, post-axial polydactyly, heart defects, hypomastia, cleft palate/bifid uvula, progressive scoliosis, and structural brain abnormalities. Four females from our cohort were identified by targeted genetic testing because their phenotype was suggestive for USP9X mutations. In several females, pigment changes along Blaschko lines and body asymmetry were observed, which is probably related to differential (escape from) X-inactivation between tissues. Expression studies on both mRNA and protein level in affected-female-derived fibroblasts showed significant reduction of USP9X level, confirming the loss-of-function effect of the identified mutations. Given that some features of affected females are also reported in known ciliopathy syndromes, we examined the role of USP9X in the primary cilium and found that endogenous USP9X localizes along the length of the ciliary axoneme, indicating that its loss of function could indeed disrupt cilium-regulated processes. Absence of dysregulated ciliary parameters in affected female-derived fibroblasts, however, points toward spatiotemporal specificity of ciliary USP9X (dys-)function.

Intellectual disability (ID) affects approximately 1%-3% of humans with a gender bias toward males. Previous studies have identified mutations in more than 100 genes on the X chromosome in males with ID, but there is less evidence for de novo mutations on the X chromosome causing ID in females. In this study we present 35 unique deleterious de novo mutations in DDX3X identified by whole exome sequencing in 38 females with ID and various other features including hypotonia, movement disorders, behavior problems, corpus callosum hypoplasia, and epilepsy. Based on our findings, mutations in DDX3X are one of the more common causes of ID, accounting for 1%-3% of unexplained ID in females. Although no de novo DDX3X mutations were identified in males, we present three families with segregating missense mutations in DDX3X, suggestive of an X-linked recessive inheritance pattern. In these families, all males with the DDX3X variant had ID, whereas carrier females were unaffected. To explore the pathogenic mechanisms accounting for the differences in disease transmission and phenotype between affected females and affected males with DDX3X missense variants, we used canonical Wnt defects in zebrafish as a surrogate measure of DDX3X function in vivo. We demonstrate a consistent loss-of-function effect of all tested de novo mutations on the Wnt pathway, and we further show a differential effect by gender. The differential activity possibly reflects a dose-dependent effect of DDX3X expression in the context of functional mosaic females versus one-copy males, which reflects the complex biological nature of DDX3X mutations.

Revertant mosaicism is an infrequently observed phenomenon caused by spontaneous correction of a pathogenic allele. We have observed such reversions caused by mitotic recombination of mutant TERC (telomerase RNA component) alleles in six patients from four families affected by dyskeratosis congenita (DC). DC is a multisystem disorder characterized by mucocutaneous abnormalities, dystrophic nails, bone-marrow failure, lung fibrosis, liver cirrhosis, and cancer. We identified a 4 nt deletion in TERC in a family with an autosomal-dominant form of DC. In two affected brothers without bone-marrow failure, sequence analysis revealed pronounced overrepresentation of the wild-type allele in blood cells, whereas no such skewing was observed in the other tissues tested. These observations suggest that this mosaic pattern might have resulted from somatic reversion of the mutated allele to the normal allele in blood-forming cells. SNP-microarray analysis on blood DNA from the two brothers indeed showed independent events of acquired segmental isodisomy of chromosome 3q, including TERC, indicating that the reversions must have resulted from mitotic recombination events. Subsequently, after developing a highly sensitive method of detecting mosaic homozygosity, we have found four additional cases with a mosaic-reversion pattern in blood cells; these four cases are part of a cohort of 17 individuals with germline TERC mutations. This shows that revertant mosaicism is a recurrent event in DC. This finding has important implications for improving diagnostic testing and understanding the variable phenotype of DC.

Aberrant activation or inhibition of potassium (K+) currents across the plasma membrane of cells has been causally linked to altered neurotransmission, cardiac arrhythmias, endocrine dysfunction, and (more rarely) perturbed developmental processes. The K+ channel subfamily K member 4 (KCNK4), also known as TRAAK (TWIK-related arachidonic acid-stimulated K+ channel), belongs to the mechano-gated ion channels of the TRAAK/TREK subfamily of two-pore-domain (K2P) K+ channels. While K2P channels are well known to contribute to the resting membrane potential and cellular excitability, their involvement in pathophysiological processes remains largely uncharacterized. We report that de novo missense mutations in KCNK4 cause a recognizable syndrome with a distinctive facial gestalt, for which we propose the acronym FHEIG (facial dysmorphism, hypertrichosis, epilepsy, intellectual disability/developmental delay, and gingival overgrowth). Patch-clamp analyses documented a significant gain of function of the identified KCNK4 channel mutants basally and impaired sensitivity to mechanical stimulation and arachidonic acid. Co-expression experiments indicated a dominant behavior of the disease-causing mutations. Molecular dynamics simulations consistently indicated that mutations favor sealing of the lateral intramembrane fenestration that has been proposed to negatively control K+ flow by allowing lipid access to the central cavity of the channel. Overall, our findings illustrate the pleiotropic effect of dysregulated KCNK4 function and provide support to the hypothesis of a gating mechanism based on the lateral fenestrations of K2P channels.

Chromatin remodeling is of crucial importance during brain development. Pathogenic alterations of several chromatin remodeling ATPases have been implicated in neurodevelopmental disorders. We describe an index case with a de novo missense mutation in CHD3, identified during whole genome sequencing of a cohort of children with rare speech disorders. To gain a comprehensive view of features associated with disruption of this gene, we use a genotype-driven approach, collecting and characterizing 35 individuals with de novo CHD3 mutations and overlapping phenotypes. Most mutations cluster within the ATPase/helicase domain of the encoded protein. Modeling their impact on the three-dimensional structure demonstrates disturbance of critical binding and interaction motifs. Experimental assays with six of the identified mutations show that a subset directly affects ATPase activity, and all but one yield alterations in chromatin remodeling. We implicate de novo CHD3 mutations in a syndrome characterized by intellectual disability, macrocephaly, and impaired speech and language.

Histones mediate dynamic packaging of nuclear DNA in chromatin, a process that is precisely controlled to guarantee efficient compaction of the genome and proper chromosomal segregation during cell division and to accomplish DNA replication, transcription, and repair. Due to the important structural and regulatory roles played by histones, it is not surprising that histone functional dysregulation or aberrant levels of histones can have severe consequences for multiple cellular processes and ultimately might affect development or contribute to cell transformation. Recently, germline frameshift mutations involving the C-terminal tail of HIST1H1E, which is a widely expressed member of the linker histone family and facilitates higher-order chromatin folding, have been causally linked to an as-yet poorly defined syndrome that includes intellectual disability. We report that these mutations result in stable proteins that reside in the nucleus, bind to chromatin, disrupt proper compaction of DNA, and are associated with a specific methylation pattern. Cells expressing these mutant proteins have a dramatically reduced proliferation rate and competence, hardly enter into the S phase, and undergo accelerated senescence. Remarkably, clinical assessment of a relatively large cohort of subjects sharing these mutations revealed a premature aging phenotype as a previously unrecognized feature of the disorder. Our findings identify a direct link between aberrant chromatin remodeling, cellular senescence, and accelerated aging.

The human RNA helicase DDX6 is an essential component of membrane-less organelles called processing bodies (PBs). PBs are involved in mRNA metabolic processes including translational repression via coordinated storage of mRNAs. Previous studies in human cell lines have implicated altered DDX6 in molecular and cellular dysfunction, but clinical consequences and pathogenesis in humans have yet to be described. Here, we report the identification of five rare de novo missense variants in DDX6 in probands presenting with intellectual disability, developmental delay, and similar dysmorphic features including telecanthus, epicanthus, arched eyebrows, and low-set ears. All five missense variants (p.His372Arg, p.Arg373Gln, p.Cys390Arg, p.Thr391Ile, and p.Thr391Pro) are located in two conserved motifs of the RecA-2 domain of DDX6 involved in RNA binding, helicase activity, and protein-partner binding. We use functional studies to demonstrate that the first variants identified (p.Arg373Gln and p.Cys390Arg) cause significant defects in PB assembly in primary fibroblast and model human cell lines. These variants' interactions with several protein partners were also disrupted in immunoprecipitation assays. Further investigation via complementation assays included the additional variants p.Thr391Ile and p.Thr391Pro, both of which, similarly to p.Arg373Gln and p.Cys390Arg, demonstrated significant defects in P-body assembly. Complementing these molecular findings, modeling of the variants on solved protein structures showed distinct spatial clustering near known protein binding regions. Collectively, our clinical and molecular data describe a neurodevelopmental syndrome associated with pathogenic missense variants in DDX6. Additionally, we suggest DDX6 join the DExD/H-box genes DDX3X and DHX30 in an emerging class of neurodevelopmental disorders involving RNA helicases.

We report 15 individuals with de novo pathogenic variants in WDR26. Eleven of the individuals carry loss-of-function mutations, and four harbor missense substitutions. These 15 individuals comprise ten females and five males, and all have intellectual disability with delayed speech, a history of febrile and/or non-febrile seizures, and a wide-based, spastic, and/or stiff-legged gait. These subjects share a set of common facial features that include a prominent maxilla and upper lip that readily reveal the upper gingiva, widely spaced teeth, and a broad nasal tip. Together, these features comprise a recognizable facial phenotype. We compared these features with those of chromosome 1q41q42 microdeletion syndrome, which typically contains WDR26, and noted that clinical features are consistent between the two subsets, suggesting that haploinsufficiency of WDR26 contributes to the pathology of 1q41q42 microdeletion syndrome. Consistent with this, WDR26 loss-of-function single-nucleotide mutations identified in these subjects lead to nonsense-mediated decay with subsequent reduction of RNA expression and protein levels. We derived a structural model of WDR26 and note that missense variants identified in these individuals localize to highly conserved residues of this WD-40-repeat-containing protein. Given that WDR26 mutations have been identified in ∼1 in 2,000 of subjects in our clinical cohorts and that WDR26 might be poorly annotated in exome variant-interpretation pipelines, we would anticipate that this disorder could be more common than currently appreciated.

KCNQ5 is a highly conserved gene encoding an important channel for neuronal function; it is widely expressed in the brain and generates M-type current. Exome sequencing identified de novo heterozygous missense mutations in four probands with intellectual disability, abnormal neurological findings, and treatment-resistant epilepsy (in two of four). Comprehensive analysis of this potassium channel for the four variants expressed in frog oocytes revealed shifts in the voltage dependence of activation, including altered activation and deactivation kinetics. Specifically, both loss-of-function and gain-of-function KCNQ5 mutations, associated with increased excitability and decreased repolarization reserve, lead to pathophysiology.

Recurrent somatic variants in SPOP are cancer specific; endometrial and prostate cancers result from gain-of-function and dominant-negative effects toward BET proteins, respectively. By using clinical exome sequencing, we identified six de novo pathogenic missense variants in SPOP in seven individuals with developmental delay and/or intellectual disability, facial dysmorphisms, and congenital anomalies. Two individuals shared craniofacial dysmorphisms, including congenital microcephaly, that were strikingly different from those of the other five individuals, who had (relative) macrocephaly and hypertelorism. We measured the effect of SPOP variants on BET protein amounts in human Ishikawa endometrial cancer cells and patient-derived cell lines because we hypothesized that variants would lead to functional divergent effects on BET proteins. The de novo variants c.362G>A (p.Arg121Gln) and c. 430G>A (p.Asp144Asn), identified in the first two individuals, resulted in a gain of function, and conversely, the c.73A>G (p.Thr25Ala), c.248A>G (p.Tyr83Cys), c.395G>T (p.Gly132Val), and c.412C>T (p.Arg138Cys) variants resulted in a dominant-negative effect. Our findings suggest that these opposite functional effects caused by the variants in SPOP result in two distinct and clinically recognizable syndromic forms of intellectual disability with contrasting craniofacial dysmorphisms.

Next-generation sequencing is a powerful tool for the discovery of genes related to neurodevelopmental disorders (NDDs). Here, we report the identification of a distinct syndrome due to de novo or inherited heterozygous mutations in Tousled-like kinase 2 (TLK2) in 38 unrelated individuals and two affected mothers, using whole-exome and whole-genome sequencing technologies, matchmaker databases, and international collaborations. Affected individuals had a consistent phenotype, characterized by mild-borderline neurodevelopmental delay (86%), behavioral disorders (68%), severe gastro-intestinal problems (63%), and facial dysmorphism including blepharophimosis (82%), telecanthus (74%), prominent nasal bridge (68%), broad nasal tip (66%), thin vermilion of the upper lip (62%), and upslanting palpebral fissures (55%). Analysis of cell lines from three affected individuals showed that mutations act through a loss-of-function mechanism in at least two case subjects. Genotype-phenotype analysis and comparison of computationally modeled faces showed that phenotypes of these and other individuals with loss-of-function variants significantly overlapped with phenotypes of individuals with other variant types (missense and C-terminal truncating). This suggests that haploinsufficiency of TLK2 is the most likely underlying disease mechanism, leading to a consistent neurodevelopmental phenotype. This work illustrates the power of international data sharing, by the identification of 40 individuals from 26 different centers in 7 different countries, allowing the identification, clinical delineation, and genotype-phenotype evaluation of a distinct NDD caused by mutations in TLK2.

We describe six persons from three families with three homozygous protein truncating variants in PUS7: c.89_90del (p.Thr30Lysfs∗20), c.1348C>T (p.Arg450∗), and a deletion of the penultimate exon 15. All these individuals have intellectual disability with speech delay, short stature, microcephaly, and aggressive behavior. PUS7 encodes the RNA-independent pseudouridylate synthase 7. Pseudouridylation is the most abundant post-transcriptional modification in RNA, which is primarily thought to stabilize secondary structures of RNA. We show that the disease-related variants lead to abolishment of PUS7 activity on both tRNA and mRNA substrates. Moreover, pus7 knockout in Drosophila melanogaster results in a number of behavioral defects, including increased activity, disorientation, and aggressiveness supporting that neurological defects are caused by PUS7 variants. Our findings demonstrate that RNA pseudouridylation by PUS7 is essential for proper neuronal development and function.

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.

Calcium (Ca2+) is a universal second messenger involved in synaptogenesis and cell survival; consequently, its regulation is important for neurons. ATPase plasma membrane Ca2+ transporting 1 (ATP2B1) belongs to the family of ATP-driven calmodulin-dependent Ca2+ pumps that participate in the regulation of intracellular free Ca2+. Here, we clinically describe a cohort of 12 unrelated individuals with variants in ATP2B1 and an overlapping phenotype of mild to moderate global development delay. Additional common symptoms include autism, seizures, and distal limb abnormalities. Nine probands harbor missense variants, seven of which were in specific functional domains, and three individuals have nonsense variants. 3D structural protein modeling suggested that the variants have a destabilizing effect on the protein. We performed Ca2+ imaging after introducing all nine missense variants in transfected HEK293 cells and showed that all variants lead to a significant decrease in Ca2+ export capacity compared with the wild-type construct, thus proving their pathogenicity. Furthermore, we observed for the same variant set an incorrect intracellular localization of ATP2B1. The genetic findings and the overlapping phenotype of the probands as well as the functional analyses imply that de novo variants in ATP2B1 lead to a monogenic form of neurodevelopmental disorder.

Neurodevelopmental disorders are highly heterogenous conditions resulting from abnormalities of brain architecture and/or function. FBXW7 (F-box and WD-repeat-domain-containing 7), a recognized developmental regulator and tumor suppressor, has been shown to regulate cell-cycle progression and cell growth and survival by targeting substrates including CYCLIN E1/2 and NOTCH for degradation via the ubiquitin proteasome system. We used a genotype-first approach and global data-sharing platforms to identify 35 individuals harboring de novo and inherited FBXW7 germline monoallelic chromosomal deletions and nonsense, frameshift, splice-site, and missense variants associated with a neurodevelopmental syndrome. The FBXW7 neurodevelopmental syndrome is distinguished by global developmental delay, borderline to severe intellectual disability, hypotonia, and gastrointestinal issues. Brain imaging detailed variable underlying structural abnormalities affecting the cerebellum, corpus collosum, and white matter. A crystal-structure model of FBXW7 predicted that missense variants were clustered at the substrate-binding surface of the WD40 domain and that these might reduce FBXW7 substrate binding affinity. Expression of recombinant FBXW7 missense variants in cultured cells demonstrated impaired CYCLIN E1 and CYCLIN E2 turnover. Pan-neuronal knockdown of the Drosophila ortholog, archipelago, impaired learning and neuronal function. Collectively, the data presented herein provide compelling evidence of an F-Box protein-related, phenotypically variable neurodevelopmental disorder associated with monoallelic variants in FBXW7.

Missense and truncating variants in the X-chromosome-linked CLCN4 gene, resulting in reduced or complete loss-of-function (LOF) of the encoded chloride/proton exchanger ClC-4, were recently demonstrated to cause a neurocognitive phenotype in both males and females. Through international clinical matchmaking and interrogation of public variant databases we assembled a database of 90 rare CLCN4 missense variants in 90 families: 41 unique and 18 recurrent variants in 49 families. For 43 families, including 22 males and 33 females, we collated detailed clinical and segregation data. To confirm causality of variants and to obtain insight into disease mechanisms, we investigated the effect on electrophysiological properties of 59 of the variants in Xenopus oocytes using extended voltage and pH ranges. Detailed analyses revealed new pathophysiological mechanisms: 25% (15/59) of variants demonstrated LOF, characterized by a "shift" of the voltage-dependent activation to more positive voltages, and nine variants resulted in a toxic gain-of-function, associated with a disrupted gate allowing inward transport at negative voltages. Functional results were not always in line with in silico pathogenicity scores, highlighting the complexity of pathogenicity assessment for accurate genetic counselling. The complex neurocognitive and psychiatric manifestations of this condition, and hitherto under-recognized impacts on growth, gastrointestinal function, and motor control are discussed. Including published cases, we summarize features in 122 individuals from 67 families with CLCN4-related neurodevelopmental condition and suggest future research directions with the aim of improving the integrated care for individuals with this diagnosis.

THOC6 variants are the genetic basis of autosomal recessive THOC6 Intellectual Disability Syndrome (TIDS). THOC6 is critical for mammalian Transcription Export complex (TREX) tetramer formation, which is composed of four six-subunit THO monomers. The TREX tetramer facilitates mammalian RNA processing, in addition to the nuclear mRNA export functions of the TREX dimer conserved through yeast. Human and mouse TIDS model systems revealed novel THOC6-dependent, species-specific TREX tetramer functions. Germline biallelic Thoc6 loss-of-function (LOF) variants result in mouse embryonic lethality. Biallelic THOC6 LOF variants reduce the binding affinity of ALYREF to THOC5 without affecting the protein expression of TREX members, implicating impaired TREX tetramer formation. Defects in RNA nuclear export functions were not detected in biallelic THOC6 LOF human neural cells. Instead, mis-splicing was detected in human and mouse neural tissue, revealing novel THOC6-mediated TREX coordination of mRNA processing. We demonstrate that THOC6 is required for key signaling pathways known to regulate the transition from proliferative to neurogenic divisions during human corticogenesis. Together, these findings implicate altered RNA processing in the developmental biology of TIDS neuropathology.

Clusters of functionally related genes can be disrupted by a single copy number variant (CNV). We demonstrate that the simultaneous disruption of multiple functionally related genes is a frequent and significant characteristic of de novo CNVs in patients with developmental disorders (P = 1 × 10(-3)). Using three different functional networks, we identified unexpectedly large numbers of functionally related genes within de novo CNVs from two large independent cohorts of individuals with developmental disorders. The presence of multiple functionally related genes was a significant predictor of a CNV's pathogenicity when compared to CNVs from apparently healthy individuals and a better predictor than the presence of known disease or haploinsufficient genes for larger CNVs. The functionally related genes found in the de novo CNVs belonged to 70% of all clusters of functionally related genes found across the genome. De novo CNVs were more likely to affect functional clusters and affect them to a greater extent than benign CNVs (P = 6 × 10(-4)). Furthermore, such clusters of functionally related genes are phenotypically informative: Different patients possessing CNVs that affect the same cluster of functionally related genes exhibit more similar phenotypes than expected (P < 0.05). The spanning of multiple functionally similar genes by single CNVs contributes substantially to how these variants exert their pathogenic effects.

Intellectual disability (ID) and autism spectrum disorders (ASD) are genetically heterogeneous, and a significant number of genes have been associated with both conditions. A few mutations in POGZ have been reported in recent exome studies; however, these studies do not provide detailed clinical information. We collected the clinical and molecular data of 25 individuals with disruptive mutations in POGZ by diagnostic whole-exome, whole-genome, or targeted sequencing of 5,223 individuals with neurodevelopmental disorders (ID primarily) or by targeted resequencing of this locus in 12,041 individuals with ASD and/or ID. The rarity of disruptive mutations among unaffected individuals (2/49,401) highlights the significance (p = 4.19 × 10(-13); odds ratio = 35.8) and penetrance (65.9%) of this genetic subtype with respect to ASD and ID. By studying the entire cohort, we defined common phenotypic features of POGZ individuals, including variable levels of developmental delay (DD) and more severe speech and language delay in comparison to the severity of motor delay and coordination issues. We also identified significant associations with vision problems, microcephaly, hyperactivity, a tendency to obesity, and feeding difficulties. Some features might be explained by the high expression of POGZ, particularly in the cerebellum and pituitary, early in fetal brain development. We conducted parallel studies in Drosophila by inducing conditional knockdown of the POGZ ortholog row, further confirming that dosage of POGZ, specifically in neurons, is essential for normal learning in a habituation paradigm. Combined, the data underscore the pathogenicity of loss-of-function mutations in POGZ and define a POGZ-related phenotype enriched in specific features.