Prevalence of comorbidities in individuals with neurodevelopmental disorders from the aggregated phenomics data of 51,227 pediatric individuals.

The prevalence of comorbidities in individuals with neurodevelopmental disorders (NDDs) is not well understood, yet these are important for accurate diagnosis and prognosis in routine care and for characterizing the clinical spectrum of NDD syndromes. We thus developed PhenomAD-NDD, an aggregated database containing the comorbid phenotypic data of 51,227 individuals with NDD, all harmonized into Human Phenotype Ontology (HPO), with in total 3,054 unique HPO terms. We demonstrate that almost all congenital anomalies are more prevalent in the NDD population than in the general population, and the NDD baseline prevalence allows for an approximation of the enrichment of symptoms. For example, such analyses of 33 genetic NDDs show that 32% of enriched phenotypes are currently not reported in the clinical synopsis in the Online Mendelian Inheritance in Man (OMIM). PhenomAD-NDD is open to all via a visualization online tool and allows us to determine the enrichment of symptoms in NDD.

Expanding the PRAAS spectrum: De novo mutations of immunoproteasome subunit ß-type 10 in six infants with SCID-Omenn syndrome.

Mutations in proteasome ß-subunits or their chaperone and regulatory proteins are associated with proteasome-associated autoinflammatory disorders (PRAAS). We studied six unrelated infants with three de novo heterozygous missense variants in PSMB10, encoding the proteasome ß2i-subunit. Individuals presented with T-B-NK± severe combined immunodeficiency (SCID) and clinical features suggestive of Omenn syndrome, including diarrhea, alopecia, and desquamating erythematous rash. Remaining T cells had limited T cell receptor repertoires, a skewed memory phenotype, and an elevated CD4/CD8 ratio. Bone marrow examination indicated severely impaired B cell maturation with limited V(D)J recombination. All infants received an allogeneic stem cell transplant and exhibited a variety of severe inflammatory complications thereafter, with 2 peri-transplant and 2 delayed deaths. The single long-term transplant survivor showed evidence for genetic rescue through revertant mosaicism overlapping the affected PSMB10 locus. The identified variants (c.166G>C [p.Asp56His] and c.601G>A/c.601G>C [p.Gly201Arg]) were predicted in silico to profoundly disrupt 20S immunoproteasome structure through impaired ß-ring/ß-ring interaction. Our identification of PSMB10 mutations as a cause of SCID-Omenn syndrome reinforces the connection between PRAAS-related diseases and SCID.

De novo variants predicting haploinsufficiency for DIP2C are associated with expressive speech delay.

The disconnected (disco)-interacting protein 2 (DIP2) gene was first identified in D. melanogaster and contains a DNA methyltransferase-associated protein 1 (DMAP1) binding domain, Acyl-CoA synthetase domain and AMP-binding sites. DIP2 regulates axonal bifurcation of the mushroom body neurons in D. melanogaster and is required for axonal regeneration in the neurons of C. elegans. The DIP2 homologues in vertebrates, Disco-interacting protein 2 homolog A (DIP2A), Disco-interacting protein 2 homolog B (DIP2B), and Disco-interacting protein 2 homolog C (DIP2C), are highly conserved and expressed widely in the central nervous system. Although there is evidence that DIP2C plays a role in cognition, reports of pathogenic variants in these genes are rare and their significance is uncertain. We present 23 individuals with heterozygous DIP2C variants, all manifesting developmental delays that primarily affect expressive language and speech articulation. Eight patients had de novo variants predicting loss-of-function in the DIP2C gene, two patients had de novo missense variants, three had paternally inherited loss of function variants and six had maternally inherited loss-of-function variants, while inheritance was unknown for four variants. Four patients had cardiac defects (hypertrophic cardiomyopathy, atrial septal defects, and bicuspid aortic valve). Minor facial anomalies were inconsistent but included a high anterior hairline with a long forehead, broad nasal tip, and ear anomalies. Brainspan analysis showed elevated DIP2C expression in the human neocortex at 10-24 weeks after conception. With the cases presented herein, we provide phenotypic and genotypic data supporting the association between loss-of-function variants in DIP2C with a neurocognitive phenotype.

TREX tetramer disruption alters RNA processing necessary for corticogenesis in THOC6 Intellectual Disability Syndrome.

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.

Refining the 9q34.3 microduplication syndrome reveals mild neurodevelopmental features associated with a distinct global DNA methylation profile.

Precise regulation of gene expression is important for correct neurodevelopment. 9q34.3 deletions affecting the EHMT1 gene result in a syndromic neurodevelopmental disorder named Kleefstra syndrome. In contrast, duplications of the 9q34.3 locus encompassing EHMT1 have been suggested to cause developmental disorders, but only limited information has been available. We have identified 15 individuals from 10 unrelated families, with 9q34.3 duplications <1.5 Mb in size, encompassing EHMT1 entirely. Clinical features included mild developmental delay, mild intellectual disability or learning problems, autism spectrum disorder, and behavior problems. The individuals did not consistently display dysmorphic features, congenital anomalies, or growth abnormalities. DNA methylation analysis revealed a weak DNAm profile for the cases with 9q34.3 duplication encompassing EHMT1, which could segregate the majority of the affected cases from controls. This study shows that individuals with 9q34.3 duplications including EHMT1 gene present with mild non-syndromic neurodevelopmental disorders and DNA methylation changes different from Kleefstra syndrome.

USP27X variants underlying X-linked intellectual disability disrupt protein function via distinct mechanisms.

Neurodevelopmental disorders with intellectual disability (ND/ID) are a heterogeneous group of diseases driving lifelong deficits in cognition and behavior with no definitive cure. X-linked intellectual disability disorder 105 (XLID105, #300984; OMIM) is a ND/ID driven by hemizygous variants in the USP27X gene encoding a protein deubiquitylase with a role in cell proliferation and neural development. Currently, only four genetically diagnosed individuals from two unrelated families have been described with limited clinical data. Furthermore, the mechanisms underlying the disorder are unknown. Here, we report 10 new XLID105 individuals from nine families and determine the impact of gene variants on USP27X protein function. Using a combination of clinical genetics, bioinformatics, biochemical, and cell biology approaches, we determined that XLID105 variants alter USP27X protein biology via distinct mechanisms including changes in developmentally relevant protein-protein interactions and deubiquitylating activity. Our data better define the phenotypic spectrum of XLID105 and suggest that XLID105 is driven by USP27X functional disruption. Understanding the pathogenic mechanisms of XLID105 variants will provide molecular insight into USP27X biology and may create the potential for therapy development.

Variant-specific pathophysiological mechanisms of AFF3 differently influence transcriptome profiles.

Background: We previously described the KINSSHIP syndrome, an autosomal dominant disorder associated with intellectual disability (ID), mesomelic dysplasia and horseshoe kidney,caused by de novo variants in the degron of AFF3. Mouse knock-ins and overexpression in zebrafish provided evidence for a dominant-negative (DN) mode-of-action, wherein an increased level of AFF3 resulted in pathological effects. Methods: Evolutionary constraints suggest that other mode-of-inheritance could be at play. We challenged this hypothesis by screening ID cohorts for individuals with predicted-to-be deleterious variants in AFF3. We used both animal and cellular models to assess the deleteriousness of the identified variants. Results: We identified an individual with a KINSSHIP-like phenotype carrying a de novo partial duplication of AFF3 further strengthening the hypothesis that an increased level of AFF3 is pathological. We also detected seventeen individuals displaying a milder syndrome with either heterozygous LoF or biallelic missense variants in AFF3. Consistent with semi-dominance, we discovered three patients with homozygous LoF and one compound heterozygote for a LoF and a missense variant, who presented more severe phenotypes than their heterozygous parents. Matching zebrafish knockdowns exhibit neurological defects that could be rescued by expressing human AFF3 mRNA, confirming their association with the ablation of aff3. Conversely, some of the human AFF3 mRNAs carrying missense variants identified in affected individuals did not complement. Overexpression of mutated AFF3 mRNAs in zebrafish embryos produced a significant increase of abnormal larvae compared to wild-type overexpression further demonstrating deleteriousness. To further assess the effect of AFF3 variation, we profiled the transcriptome of fibroblasts from affected individuals and engineered isogenic cells harboring +/+, DN/DN, LoF/+, LoF/LoF or DN/LoF AFF3 genotypes. The expression of more than a third of the AFF3 bound loci is modified in either the DN/DN or the LoF/LoF lines. While the same pathways are affected, only about one-third of the differentially expressed genes are common to these homozygote datasets, indicating that AFF3 LoF and DN variants largely modulate transcriptomes differently, e.g. the DNA repair pathway displayed opposite modulation. Conclusions: Our results and the high pleiotropy shown by variation at this locus suggest that minute changes in AFF3 function are deleterious.

Further delineation of the rare GDACCF (global developmental delay, absent or hypoplastic corpus callosum, dysmorphic facies syndrome): genotype and phenotype of 22 patients with ZNF148 mutations.

BACKGROUND: Pathogenic variants in the zinc finger protein coding genes are rare causes of intellectual disability and congenital malformations. Mutations in the ZNF148 gene causing GDACCF syndrome (global developmental delay, absent or hypoplastic corpus callosum, dysmorphic facies; MIM #617260) have been reported in five individuals so far. METHODS: As a result of an international collaboration using GeneMatcher Phenome Central Repository and personal communications, here we describe the clinical and molecular genetic characteristics of 22 previously unreported individuals. RESULTS: The core clinical phenotype is characterised by developmental delay particularly in the domain of speech development, postnatal growth retardation, microcephaly and facial dysmorphism. Corpus callosum abnormalities appear less frequently than suggested by previous observations. The identified mutations concerned nonsense or frameshift variants that were mainly located in the last exon of the ZNF148 gene. Heterozygous deletion including the entire ZNF148 gene was found in only one case. Most mutations occurred de novo, but were inherited from an affected parent in two families. CONCLUSION: The GDACCF syndrome is clinically diverse, and a genotype-first approach, that is, exome sequencing is recommended for establishing a genetic diagnosis rather than a phenotype-first approach. However, the syndrome may be suspected based on some recurrent, recognisable features. Corpus callosum anomalies were not as constant as previously suggested, we therefore recommend to replace the term 'GDACCF syndrome' with 'ZNF148-related neurodevelopmental disorder'.

Mobile element insertions in rare diseases: a comparative benchmark and reanalysis of 60,000 exome samples.

Mobile element insertions (MEIs) are a known cause of genetic disease but have been underexplored due to technical limitations of genetic testing methods. Various bioinformatic tools have been developed to identify MEIs in Next Generation Sequencing data. However, most tools have been developed specifically for genome sequencing (GS) data rather than exome sequencing (ES) data, which remains more widely used for routine diagnostic testing. In this study, we benchmarked six MEI detection tools (ERVcaller, MELT, Mobster, SCRAMble, TEMP2 and xTea) on ES data and on GS data from publicly available genomic samples (HG002, NA12878). For all the tools we evaluated sensitivity and precision of different filtering strategies. Results show that there were substantial differences in tool performance between ES and GS data. MELT performed best with ES data and its combination with SCRAMble increased substantially the detection rate of MEIs. By applying both tools to 10,890 ES samples from Solve-RD and 52,624 samples from Radboudumc we were able to diagnose 10 patients who had remained undiagnosed by conventional ES analysis until now. Our study shows that MELT and SCRAMble can be used reliably to identify clinically relevant MEIs in ES data. This may lead to an additional diagnosis for 1 in 3000 to 4000 patients in routine clinical ES.

A recurrent de novo MAX p.Arg60Gln variant causes a syndromic overgrowth disorder through differential expression of c-Myc target genes.

Cyclin D2 (CCND2) stabilization underpins a range of macrocephaly-associated disorders through mutation of CCND2 or activating mutations in upstream genes encoding PI3K-AKT pathway components. Here, we describe three individuals with overlapping macrocephaly-associated phenotypes who carry the same recurrent de novo c.179G>A (p.Arg60Gln) variant in Myc-associated factor X (MAX). The mutation, located in the b-HLH-LZ domain, causes increased intracellular CCND2 through increased transcription but it does not cause stabilization of CCND2. We show that the purified b-HLH-LZ domain of MAXArg60Gln (Max∗Arg60Gln) binds its target E-box sequence with a lower apparent affinity. This leads to a more efficient heterodimerization with c-Myc resulting in an increase in transcriptional activity of c-Myc in individuals carrying this mutation. The recent development of Omomyc-CPP, a cell-penetrating b-HLH-LZ-domain c-Myc inhibitor, provides a possible therapeutic option for MAXArg60Gln individuals, and others carrying similar germline mutations resulting in dysregulated transcriptional c-Myc activity.

Systematic analysis of paralogous regions in 41,755 exomes uncovers clinically relevant variation.

The short lengths of short-read sequencing reads challenge the analysis of paralogous genomic regions in exome and genome sequencing data. Most genetic variants within these homologous regions therefore remain unidentified in standard analyses. Here, we present a method (Chameleolyser) that accurately identifies single nucleotide variants and small insertions/deletions (SNVs/Indels), copy number variants and ectopic gene conversion events in duplicated genomic regions using whole-exome sequencing data. Application to a cohort of 41,755 exome samples yields 20,432 rare homozygous deletions and 2,529,791 rare SNVs/Indels, of which we show that 338,084 are due to gene conversion events. None of the SNVs/Indels are detectable using regular analysis techniques. Validation by high-fidelity long-read sequencing in 20 samples confirms >88% of called variants. Focusing on variation in known disease genes leads to a direct molecular diagnosis in 25 previously undiagnosed patients. Our method can readily be applied to existing exome data.

Biallelic MAD2L1BP (p31comet) mutation is associated with mosaic aneuploidy and juvenile granulosa cell tumors.

MAD2L1BP-encoded p31comet mediates Trip13-dependent disassembly of Mad2- and Rev7-containing complexes and, through this antagonism, promotes timely spindle assembly checkpoint (SAC) silencing, faithful chromosome segregation, insulin signaling, and homology-directed repair (HDR) of DNA double-strand breaks. We identified a homozygous MAD2L1BP nonsense variant, R253*, in 2 siblings with microcephaly, epileptic encephalopathy, and juvenile granulosa cell tumors of ovary and testis. Patient-derived cells exhibited high-grade mosaic variegated aneuploidy, slowed-down proliferation, and instability of truncated p31comet mRNA and protein. Corresponding recombinant p31comet was defective in Trip13, Mad2, and Rev7 binding and unable to support SAC silencing or HDR. Furthermore, C-terminal truncation abrogated an identified interaction of p31comet with tp53. Another homozygous truncation, R227*, detected in an early-deceased patient with low-level aneuploidy, severe epileptic encephalopathy, and frequent blood glucose elevations, likely corresponds to complete loss of function, as in Mad2l1bp-/- mice. Thus, human mutations of p31comet are linked to aneuploidy and tumor predisposition.

Mechanisms of mRNA processing defects in inherited THOC6 intellectual disability syndrome.

THOC6 is the genetic basis of autosomal recessive THOC6 Intellectual Disability Syndrome (TIDS). THOC6 facilitates the formation of the Transcription Export complex (TREX) tetramer, composed of four THO monomers. The TREX tetramer supports mammalian mRNA processing that is distinct from yeast TREX dimer functions. Human and mouse TIDS model systems allow novel THOC6-dependent TREX tetramer functions to be investigated. Biallelic loss-of-functon(LOF) THOC6 variants do not influence the expression and localization of TREX members in human cells, but our data suggests reduced binding affinity of ALYREF. Impairment of TREX nuclear export functions were not detected in cells with biallelic THOC6 LOF. Instead, mRNA mis-splicing was observed in human and mouse neural tissue, revealing novel insights into THOC6-mediated TREX coordination of mRNA processing. We demonstrate that THOC6 is required for regulation of key signaling pathways in human corticogenesis that dictate the transition from proliferative to neurogenic divisions that may inform TIDS neuropathology.

PhenoScore quantifies phenotypic variation for rare genetic diseases by combining facial analysis with other clinical features using a machine-learning framework.

Several molecular and phenotypic algorithms exist that establish genotype-phenotype correlations, including facial recognition tools. However, no unified framework that investigates both facial data and other phenotypic data directly from individuals exists. We developed PhenoScore: an open-source, artificial intelligence-based phenomics framework, combining facial recognition technology with Human Phenotype Ontology data analysis to quantify phenotypic similarity. Here we show PhenoScore's ability to recognize distinct phenotypic entities by establishing recognizable phenotypes for 37 of 40 investigated syndromes against clinical features observed in individuals with other neurodevelopmental disorders and show it is an improvement on existing approaches. PhenoScore provides predictions for individuals with variants of unknown significance and enables sophisticated genotype-phenotype studies by testing hypotheses on possible phenotypic (sub)groups. PhenoScore confirmed previously known phenotypic subgroups caused by variants in the same gene for SATB1, SETBP1 and DEAF1 and provides objective clinical evidence for two distinct ADNP-related phenotypes, already established functionally.

De novo PHF5A variants are associated with craniofacial abnormalities, developmental delay, and hypospadias.

PURPOSE: The SF3B splicing complex is composed of SF3B1-6 and PHF5A. We report a developmental disorder caused by de novo variants in PHF5A. METHODS: Clinical, genomic, and functional studies using subject-derived fibroblasts and a heterologous cellular system were performed. RESULTS: We studied 9 subjects with congenital malformations, including preauricular tags and hypospadias, growth abnormalities, and developmental delay who had de novo heterozygous PHF5A variants, including 4 loss-of-function (LOF), 3 missense, 1 splice, and 1 start-loss variant. In subject-derived fibroblasts with PHF5A LOF variants, wild-type and variant PHF5A mRNAs had a 1:1 ratio, and PHF5A mRNA levels were normal. Transcriptome sequencing revealed alternative promoter use and downregulated genes involved in cell-cycle regulation. Subject and control fibroblasts had similar amounts of PHF5A with the predicted wild-type molecular weight and of SF3B1-3 and SF3B6. SF3B complex formation was unaffected in 2 subject cell lines. CONCLUSION: Our data suggest the existence of feedback mechanisms in fibroblasts with PHF5A LOF variants to maintain normal levels of SF3B components. These compensatory mechanisms in subject fibroblasts with PHF5A or SF3B4 LOF variants suggest disturbed autoregulation of mutated splicing factor genes in specific cell types, that is, neural crest cells, during embryonic development rather than haploinsufficiency as pathomechanism.

ANK2 loss-of-function variants are associated with epilepsy, and lead to impaired axon initial segment plasticity and hyperactive network activity in hiPSC-derived neuronal networks.

PURPOSE: To characterize a novel neurodevelopmental syndrome due to loss-of-function (LoF) variants in Ankyrin 2 (ANK2), and to explore the effects on neuronal network dynamics and homeostatic plasticity in human-induced pluripotent stem cell-derived neurons. METHODS: We collected clinical and molecular data of 12 individuals with heterozygous de novo LoF variants in ANK2. We generated a heterozygous LoF allele of ANK2 using CRISPR/Cas9 in human-induced pluripotent stem cells (hiPSCs). HiPSCs were differentiated into excitatory neurons, and we measured their spontaneous electrophysiological responses using micro-electrode arrays (MEAs). We also characterized their somatodendritic morphology and axon initial segment (AIS) structure and plasticity. RESULTS: We found a broad neurodevelopmental disorder (NDD), comprising intellectual disability, autism spectrum disorders and early onset epilepsy. Using MEAs, we found that hiPSC-derived neurons with heterozygous LoF of ANK2 show a hyperactive and desynchronized neuronal network. ANK2-deficient neurons also showed increased somatodendritic structures and altered AIS structure of which its plasticity is impaired upon activity-dependent modulation. CONCLUSIONS: Phenotypic characterization of patients with de novo ANK2 LoF variants defines a novel NDD with early onset epilepsy. Our functional in vitro data of ANK2-deficient human neurons show a specific neuronal phenotype in which reduced ANKB expression leads to hyperactive and desynchronized neuronal network activity, increased somatodendritic complexity and AIS structure and impaired activity-dependent plasticity of the AIS.

The clinical and molecular spectrum of the KDM6B-related neurodevelopmental disorder.

De novo variants are a leading cause of neurodevelopmental disorders (NDDs), but because every monogenic NDD is different and usually extremely rare, it remains a major challenge to understand the complete phenotype and genotype spectrum of any morbid gene. According to OMIM, heterozygous variants in KDM6B cause "neurodevelopmental disorder with coarse facies and mild distal skeletal abnormalities." Here, by examining the molecular and clinical spectrum of 85 reported individuals with mostly de novo (likely) pathogenic KDM6B variants, we demonstrate that this description is inaccurate and potentially misleading. Cognitive deficits are seen consistently in all individuals, but the overall phenotype is highly variable. Notably, coarse facies and distal skeletal anomalies, as defined by OMIM, are rare in this expanded cohort while other features are unexpectedly common (e.g., hypotonia, psychosis, etc.). Using 3D protein structure analysis and an innovative dual Drosophila gain-of-function assay, we demonstrated a disruptive effect of 11 missense/in-frame indels located in or near the enzymatic JmJC or Zn-containing domain of KDM6B. Consistent with the role of KDM6B in human cognition, we demonstrated a role for the Drosophila KDM6B ortholog in memory and behavior. Taken together, we accurately define the broad clinical spectrum of the KDM6B-related NDD, introduce an innovative functional testing paradigm for the assessment of KDM6B variants, and demonstrate a conserved role for KDM6B in cognition and behavior. Our study demonstrates the critical importance of international collaboration, sharing of clinical data, and rigorous functional analysis of genetic variants to ensure correct disease diagnosis for rare disorders.