Etiological involvement of KCND1 variants in an X-linked neurodevelopmental disorder with variable expressivity.

Utilizing trio whole-exome sequencing and a gene matching approach, we identified a cohort of 18 male individuals from 17 families with hemizygous variants in KCND1, including two de novo missense variants, three maternally inherited protein-truncating variants, and 12 maternally inherited missense variants. Affected subjects present with a neurodevelopmental disorder characterized by diverse neurological abnormalities, mostly delays in different developmental domains, but also distinct neuropsychiatric signs and epilepsy. Heterozygous carrier mothers are clinically unaffected. KCND1 encodes the α-subunit of Kv4.1 voltage-gated potassium channels. All variant-associated amino acid substitutions affect either the cytoplasmic N- or C-terminus of the channel protein except for two occurring in transmembrane segments 1 and 4. Kv4.1 channels were functionally characterized in the absence and presence of auxiliary ß subunits. Variant-specific alterations of biophysical channel properties were diverse and varied in magnitude. Genetic data analysis in combination with our functional assessment shows that Kv4.1 channel dysfunction is involved in the pathogenesis of an X-linked neurodevelopmental disorder frequently associated with a variable neuropsychiatric clinical phenotype.

The detection of a strong episignature for Chung-Jansen syndrome, partially overlapping with Börjeson-Forssman-Lehmann and White-Kernohan syndromes.

Chung-Jansen syndrome is a neurodevelopmental disorder characterized by intellectual disability, behavioral problems, obesity and dysmorphic features. It is caused by pathogenic variants in the PHIP gene that encodes for the Pleckstrin homology domain-interacting protein, which is part of an epigenetic modifier protein complex. Therefore, we hypothesized that PHIP haploinsufficiency may impact genome-wide DNA methylation (DNAm). We assessed the DNAm profiles of affected individuals with pathogenic and likely pathogenic PHIP variants with Infinium Methylation EPIC arrays and report a specific and sensitive DNAm episignature biomarker for Chung-Jansen syndrome. In addition, we observed similarities between the methylation profile of Chung-Jansen syndrome and that of functionally related and clinically partially overlapping genetic disorders, White-Kernohan syndrome (caused by variants in DDB1 gene) and Börjeson-Forssman-Lehmann syndrome (caused by variants in PHF6 gene). Based on these observations we also proceeded to develop a common episignature biomarker for these disorders. These newly defined episignatures can be used as part of a multiclass episignature classifier for screening of affected individuals with rare disorders and interpretation of genetic variants of unknown clinical significance, and provide further insights into the common molecular pathophysiology of the clinically-related Chung-Jansen, Börjeson-Forssman-Lehmann and White-Kernohan syndromes.

De novo variants in ATXN7L3 lead to developmental delay, hypotonia and distinctive facial features.

Deubiquitination is critical for the proper functioning of numerous biological pathways such as DNA repair, cell cycle progression, transcription, signal transduction, and autophagy. Accordingly, pathogenic variants in deubiquitinating enzymes (DUBs) have been implicated in neurodevelopmental disorders (ND) and congenital abnormalities. ATXN7L3 is a component of the DUB module of the SAGA complex, and two other related DUB modules, and serves as an obligate adaptor protein of 3 ubiquitin-specific proteases (USP22, USP27X or USP51). Through exome sequencing and GeneMatching, we identified nine individuals with heterozygous variants in ATXN7L3. The core phenotype included global motor and language developmental delay, hypotonia, and distinctive facial characteristics including hypertelorism, epicanthal folds, blepharoptosis, a small nose and mouth, and low-set posteriorly rotated ears. In order to assess pathogenicity, we investigated the effects of a recurrent nonsense variant [c.340C>T; p.(Arg114Ter)] in fibroblasts of an affected individual. ATXN7L3 protein levels were reduced, and deubiquitylation was impaired, as indicated by an increase in histone H2Bub1 levels. This is consistent with the previous observation of increased H2Bub1 levels in Atxn7l3-null mouse embryos, which have developmental delay and embryonic lethality. In conclusion, we present clinical information and biochemical characterization supporting ATXN7L3 variants in the pathogenesis of a rare syndromic ND.

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.

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.

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.

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.

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.

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.