De novo variants in the non-coding spliceosomal snRNA gene RNU4-2 are a frequent cause of syndromic neurodevelopmental disorders.

Around 60% of individuals with neurodevelopmental disorders (NDD) remain undiagnosed after comprehensive genetic testing, primarily of protein-coding genes1. Increasingly, large genome-sequenced cohorts are improving our ability to discover new diagnoses in the non-coding genome. Here, we identify the non-coding RNA RNU4-2 as a novel syndromic NDD gene. RNU4-2 encodes the U4 small nuclear RNA (snRNA), which is a critical component of the U4/U6.U5 tri-snRNP complex of the major spliceosome2. We identify an 18 bp region of RNU4-2 mapping to two structural elements in the U4/U6 snRNA duplex (the T-loop and Stem III) that is severely depleted of variation in the general population, but in which we identify heterozygous variants in 119 individuals with NDD. The vast majority of individuals (77.3%) have the same highly recurrent single base-pair insertion (n.64_65insT). We estimate that variants in this region explain 0.41% of individuals with NDD. We demonstrate that RNU4-2 is highly expressed in the developing human brain, in contrast to its contiguous counterpart RNU4-1 and other U4 homologs, supporting RNU4-2's role as the primary U4 transcript in the brain. Overall, this work underscores the importance of non-coding genes in rare disorders. It will provide a diagnosis to thousands of individuals with NDD worldwide and pave the way for the development of effective treatments for these individuals.

Narrowing the diagnostic gap: Genomes, episignatures, long-read sequencing, and health economic analyses in an exome-negative intellectual disability cohort.

PURPOSE: Genome sequencing (GS)-specific diagnostic rates in prospective tightly ascertained exome sequencing (ES)-negative intellectual disability (ID) cohorts have not been reported extensively. METHODS: ES, GS, epigenetic signatures, and long-read sequencing diagnoses were assessed in 74 trios with at least moderate ID. RESULTS: The ES diagnostic yield was 42 of 74 (57%). GS diagnoses were made in 9 of 32 (28%) ES-unresolved families. Repeated ES with a contemporary pipeline on the GS-diagnosed families identified 8 of 9 single-nucleotide variations/copy-number variations undetected in older ES, confirming a GS-unique diagnostic rate of 1 in 32 (3%). Episignatures contributed diagnostic information in 9% with GS corroboration in 1 of 32 (3%) and diagnostic clues in 2 of 32 (6%). A genetic etiology for ID was detected in 51 of 74 (69%) families. Twelve candidate disease genes were identified. Contemporary ES followed by GS cost US$4976 (95% CI: $3704; $6969) per diagnosis and first-line GS at a cost of $7062 (95% CI: $6210; $8475) per diagnosis. CONCLUSION: Performing GS only in ID trios would be cost equivalent to ES if GS were available at $2435, about a 60% reduction from current prices. This study demonstrates that first-line GS achieves higher diagnostic rate than contemporary ES but at a higher cost.

Integration of EpiSign, facial phenotyping, and likelihood ratio interpretation of clinical abnormalities in the re-classification of an ARID1B missense variant.

Heterozygous ARID1B variants result in Coffin-Siris syndrome. Features may include hypoplastic nails, slow growth, characteristic facial features, hypotonia, hypertrichosis, and sparse scalp hair. Most reported cases are due to ARID1B loss of function variants. We report a boy with developmental delay, feeding difficulties, aspiration, recurrent respiratory infections, slow growth, and hypotonia without a clinical diagnosis, where a previously unreported ARID1B missense variant was classified as a variant of uncertain significance. The pathogenicity of this variant was refined through combined methodologies including genome-wide methylation signature analysis (EpiSign), Machine Learning (ML) facial phenotyping, and LIRICAL. Trio exome sequencing and EpiSign were performed. ML facial phenotyping compared facial images using FaceMatch and GestaltMatcher to syndrome-specific libraries to prioritize the trio exome bioinformatic pipeline gene list output. Phenotype-driven variant prioritization was performed with LIRICAL. A de novo heterozygous missense variant, ARID1B p.(Tyr1268His), was reported as a variant of uncertain significance. The ACMG classification was refined to likely pathogenic by a supportive methylation signature, ML facial phenotyping, and prioritization through LIRICAL. The ARID1B genotype-phenotype has been expanded through an extended analysis of missense variation through genome-wide methylation signatures, ML facial phenotyping, and likelihood-ratio gene prioritization.

De novo MCM6 variants in neurodevelopmental disorders: a recognizable phenotype related to zinc binding residues.

The minichromosome maintenance (MCM) complex acts as a DNA helicase during DNA replication, and thereby regulates cell cycle progression and proliferation. In addition, MCM-complex components localize to centrosomes and play an independent role in ciliogenesis. Pathogenic variants in genes coding for MCM components and other DNA replication factors have been linked to growth and developmental disorders as Meier-Gorlin syndrome and Seckel syndrome. Trio exome/genome sequencing identified the same de novo MCM6 missense variant p.(Cys158Tyr) in two unrelated individuals that presented with overlapping phenotypes consisting of intra-uterine growth retardation, short stature, congenital microcephaly, endocrine features, developmental delay and urogenital anomalies. The identified variant affects a zinc binding cysteine in the MCM6 zinc finger signature. This domain, and specifically cysteine residues, are essential for MCM-complex dimerization and the induction of helicase activity, suggesting a deleterious effect of this variant on DNA replication. Fibroblasts derived from the two affected individuals showed defects both in ciliogenesis and cell proliferation. We additionally traced three unrelated individuals with de novo MCM6 variants in the oligonucleotide binding (OB)-fold domain, presenting with variable (neuro)developmental features including autism spectrum disorder, developmental delay, and epilepsy. Taken together, our findings implicate de novo MCM6 variants in neurodevelopmental disorders. The clinical features and functional defects related to the zinc binding residue resemble those observed in syndromes related to other MCM components and DNA replication factors, while de novo OB-fold domain missense variants may be associated with more variable neurodevelopmental phenotypes. These data encourage consideration of MCM6 variants in the diagnostic arsenal of NDD.

Dominant-negative variants in CBX1 cause a neurodevelopmental disorder.

PURPOSE: This study aimed to establish variants in CBX1, encoding heterochromatin protein 1ß (HP1ß), as a cause of a novel syndromic neurodevelopmental disorder. METHODS: Patients with CBX1 variants were identified, and clinician researchers were connected using GeneMatcher and physician referrals. Clinical histories were collected from each patient. To investigate the pathogenicity of identified variants, we performed in vitro cellular assays and neurobehavioral and cytological analyses of neuronal cells obtained from newly generated Cbx1 mutant mouse lines. RESULTS: In 3 unrelated individuals with developmental delay, hypotonia, and autistic features, we identified heterozygous de novo variants in CBX1. The identified variants were in the chromodomain, the functional domain of HP1ß, which mediates interactions with chromatin. Cbx1 chromodomain mutant mice displayed increased latency-to-peak response, suggesting the possibility of synaptic delay or myelination deficits. Cytological and chromatin immunoprecipitation experiments confirmed the reduction of mutant HP1ß binding to heterochromatin, whereas HP1ß interactome analysis demonstrated that the majority of HP1ß-interacting proteins remained unchanged between the wild-type and mutant HP1ß. CONCLUSION: These collective findings confirm the role of CBX1 in developmental disabilities through the disruption of HP1ß chromatin binding during neurocognitive development. Because HP1ß forms homodimers and heterodimers, mutant HP1ß likely sequesters wild-type HP1ß and other HP1 proteins, exerting dominant-negative effects.

De Novo ZMYND8 variants result in an autosomal dominant neurodevelopmental disorder with cardiac malformations.

PURPOSE: ZMYND8 encodes a multidomain protein that serves as a central interactive hub for coordinating critical roles in transcription regulation, chromatin remodeling, regulation of super-enhancers, DNA damage response and tumor suppression. We delineate a novel neurocognitive disorder caused by variants in the ZMYND8 gene. METHODS: An international collaboration, exome sequencing, molecular modeling, yeast two-hybrid assays, analysis of available transcriptomic data and a knockdown Drosophila model were used to characterize the ZMYND8 variants. RESULTS: ZMYND8 variants were identified in 11 unrelated individuals; 10 occurred de novo and one suspected de novo; 2 were truncating, 9 were missense, of which one was recurrent. The disorder is characterized by intellectual disability with variable cardiovascular, ophthalmologic and minor skeletal anomalies. Missense variants in the PWWP domain of ZMYND8 abolish the interaction with Drebrin and missense variants in the MYND domain disrupt the interaction with GATAD2A. ZMYND8 is broadly expressed across cell types in all brain regions and shows highest expression in the early stages of brain development. Neuronal knockdown of the DrosophilaZMYND8 ortholog results in decreased habituation learning, consistent with a role in cognitive function. CONCLUSION: We present genomic and functional evidence for disruption of ZMYND8 as a novel etiology of syndromic intellectual disability.

Germline variants in tumor suppressor FBXW7 lead to impaired ubiquitination and a neurodevelopmental syndrome.

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.

Clinically Responsive Genomic Analysis Pipelines: Elements to Improve Detection Rate and Efficiency.

Massively parallel sequencing has markedly improved mendelian diagnostic rates. This study assessed the effects of custom alterations to a diagnostic genomic bioinformatic pipeline in response to clinical need and derived practice recommendations relative to diagnostic rates and efficiency. The Genomic Annotation and Interpretation Application (GAIA) bioinformatics pipeline was designed to detect panel, exome, and genome sample integrity and prioritize gene variants in mendelian disorders. Reanalysis of selected negative cases was performed after improvements to the pipeline. GAIA improvements and their effect on sensitivity are described, including addition of a PubMed search for gene-disease associations not in the Online Mendelian Inheritance of Man database, inclusion of a process for calling low-quality variants (known as QPatch), and gene symbol nomenclature consistency checking. The new pipeline increased the diagnostic rate and reduced staff costs, resulting in a saving of US$844.34 per additional diagnosis. Recommendations for genomic analysis pipeline requirements are summarized. Clinically responsive bioinformatics pipeline improvements increase diagnostic sensitivity and increase cost-effectiveness.

Germline AGO2 mutations impair RNA interference and human neurological development.

ARGONAUTE-2 and associated miRNAs form the RNA-induced silencing complex (RISC), which targets mRNAs for translational silencing and degradation as part of the RNA interference pathway. Despite the essential nature of this process for cellular function, there is little information on the role of RISC components in human development and organ function. We identify 13 heterozygous mutations in AGO2 in 21 patients affected by disturbances in neurological development. Each of the identified single amino acid mutations result in impaired shRNA-mediated silencing. We observe either impaired RISC formation or increased binding of AGO2 to mRNA targets as mutation specific functional consequences. The latter is supported by decreased phosphorylation of a C-terminal serine cluster involved in mRNA target release, increased formation of dendritic P-bodies in neurons and global transcriptome alterations in patient-derived primary fibroblasts. Our data emphasize the importance of gene expression regulation through the dynamic AGO2-RNA association for human neuronal development.

A Mild PUM1 Mutation Is Associated with Adult-Onset Ataxia, whereas Haploinsufficiency Causes Developmental Delay and Seizures.

Certain mutations can cause proteins to accumulate in neurons, leading to neurodegeneration. We recently showed, however, that upregulation of a wild-type protein, Ataxin1, caused by haploinsufficiency of its repressor, the RNA-binding protein Pumilio1 (PUM1), also causes neurodegeneration in mice. We therefore searched for human patients with PUM1 mutations. We identified eleven individuals with either PUM1 deletions or de novo missense variants who suffer a developmental syndrome (Pumilio1-associated developmental disability, ataxia, and seizure; PADDAS). We also identified a milder missense mutation in a family with adult-onset ataxia with incomplete penetrance (Pumilio1-related cerebellar ataxia, PRCA). Studies in patient-derived cells revealed that the missense mutations reduced PUM1 protein levels by ∼25% in the adult-onset cases and by ∼50% in the infantile-onset cases; levels of known PUM1 targets increased accordingly. Changes in protein levels thus track with phenotypic severity, and identifying posttranscriptional modulators of protein expression should identify new candidate disease genes.

Integrating exome sequencing into a diagnostic pathway for epileptic encephalopathy: Evidence of clinical utility and cost effectiveness.

BACKGROUND: Epileptic encephalopathies are a devastating group of neurological conditions in which etiological diagnosis can alter management and clinical outcome. Exome sequencing and gene panel testing can improve diagnostic yield but there is no cost-effectiveness analysis of their use or consensus on how to best integrate these tests into clinical diagnostic pathways. METHODS: We conducted a retrospective cost-effectiveness study comparing trio exome sequencing with a standard diagnostic approach, for a well-phenotyped cohort of 32 patients with epileptic encephalopathy, who remained undiagnosed after "first-tier" testing. Sensitivity analysis was included with a range of commercial exome and multigene panels. RESULTS: The diagnostic yield was higher for the exome sequencing (16/32; 50%) than the standard arm (2/32; 6.2%). The trio exome sequencing pathway was cost-effective compared to the standard diagnostic pathway with a cost saving of AU$5,236 (95% confidence intervals $2,482; $9,784) per additional diagnosis; the standard pathway cost approximately 10 times more per diagnosis. Sensitivity analysis demonstrated that the majority of commercial exome sequencing and multigene panels studied were also cost-effective. The clinical utility of all diagnoses was reported. CONCLUSION: Our study supports the integration of exome sequencing and gene panel testing into the diagnostic pathway for epileptic encephalopathy, both in terms of cost effectiveness and clinical utility. We propose a diagnostic pathway that integrates initial rapid screening for treatable causes and comprehensive genomic screening. This study has important implications for health policy and public funding for epileptic encephalopathy and other neurological conditions.

A De Novo Mutation in the Sodium-Activated Potassium Channel KCNT2 Alters Ion Selectivity and Causes Epileptic Encephalopathy.

Early infantile epileptic encephalopathies (EOEE) are a debilitating spectrum of disorders associated with cognitive impairments. We present a clinical report of a KCNT2 mutation in an EOEE patient. The de novo heterozygous variant Phe240Leu SLICK was identified by exome sequencing and confirmed by Sanger sequencing. Phe240Leu rSlick and hSLICK channels were electrophysiologically, heterologously characterized to reveal three significant alterations to channel function. First, [Cl-]i sensitivity was reversed in Phe240Leu channels. Second, predominantly K+-selective WT channels were made to favor Na+ over K+ by Phe240Leu. Third, and consequent to altered ion selectivity, Phe240Leu channels had larger inward conductance. Further, rSlick channels induced membrane hyperexcitability when expressed in primary neurons, resembling the cellular seizure phenotype. Taken together, our results confirm that Phe240Leu is a "change-of-function" KCNT2 mutation, demonstrating unusual altered selectivity in KNa channels. These findings establish pathogenicity of the Phe240Leu KCNT2 mutation in the reported EOEE patient.

Comparison of whole-exome sequencing of matched fresh and formalin fixed paraffin embedded melanoma tumours: implications for clinical decision making.

The identification of recurrent driver mutations by whole-exome sequencing (WES) of fresh-frozen human cancers and the subsequent development of novel targeted therapies have recently transformed the treatment of many cancers including melanoma. In routine clinical practice, fresh-frozen tissue is rarely available and mutation testing usually needs to be carried out on archival formalin fixed, paraffin embedded (FFPE) tissue, from which DNA is typically fragmented, cross-linked and of lower quality. In this study we aimed to determine whether WES data generated from genomic DNA (gDNA) extracted from FFPE tissues can be produced reliably and of clinically-actionable standard. In this study of ten melanoma patients, we compared WES data produced from analysis of gDNA isolated from FFPE tumour tissue with that isolated from fresh-frozen tumour tissue from the same specimen. FFPE samples were sequenced using both Illumina's Nextera and NimbleGen SeqCap exome capture kits. To examine mutations between the two tissue sources and platforms, somatic mutations in the FFPE exomes were called using the matched fresh tissue sequence as a reference. Of the 10 FFPE DNA samples, seven Nextera and four SeqCap samples passed library preparation. On average, there were 5341 and 2246 variants lost in FFPE compared to matched fresh tissue utilising Nextera and SeqCap kits, respectively. In order to explore the feasibility of future clinical implementation of WES, FFPE variants in 27 genes of important clinical relevance in melanoma were assessed. The average concordance rate was 43.2% over a total of 1299 calls for the chosen genes in the FFPE DNA. For the current clinically most important melanoma mutations, 0/3 BRAF and 6/8 (75%) NRAS FFPE calls were concordant with the fresh tissue result, which was confirmed using a Sequenom OncoCarta Panel. The poor performance of FFPE WES indicates that specialised library construction to account for low quality DNA and further refinements will be necessary before this approach could be used for routine clinical decision making over currently preferred techniques.

Asparagine Synthetase Deficiency causes reduced proliferation of cells under conditions of limited asparagine.

Asparagine Synthetase Deficiency is a recently described cause of profound intellectual disability, marked progressive cerebral atrophy and variable seizure disorder. To date there has been limited functional data explaining the underlying pathophysiology. We report a new case with compound heterozygous mutations in the ASNS gene (NM_183356.3:c. [866G>C]; [1010C>T]). Both variants alter evolutionarily conserved amino acids and were predicted to be pathogenic based on in silico protein modelling that suggests disruption of the critical ATP binding site of the ASNS enzyme. In patient fibroblasts, ASNS expression as well as protein and mRNA stability are not affected by these variants. However, there is markedly reduced proliferation of patient fibroblasts when cultured in asparagine-limited growth medium, compared to parental and wild type fibroblasts. Restricting asparagine replicates the physiology within the blood-brain-barrier, with limited transfer of dietary derived asparagine, resulting in reliance of neuronal cells on intracellular asparagine synthesis by the ASNS enzyme. These functional studies offer insight into the underlying pathophysiology of the dramatic progressive cerebral atrophy associated with Asparagine Synthetase Deficiency.

A simple and novel method for RNA-seq library preparation of single cell cDNA analysis by hyperactive Tn5 transposase.

BACKGROUND: Deep sequencing of single cell-derived cDNAs offers novel insights into oncogenesis and embryogenesis. However, traditional library preparation for RNA-seq analysis requires multiple steps with consequent sample loss and stochastic variation at each step significantly affecting output. Thus, a simpler and better protocol is desirable. The recently developed hyperactive Tn5-mediated library preparation, which brings high quality libraries, is likely one of the solutions. RESULTS AND CONCLUSIONS: Here, we tested the applicability of hyperactive Tn5-mediated library preparation to deep sequencing of single cell cDNA, optimized the protocol, and compared it with the conventional method based on sonication. This new technique does not require any expensive or special equipment, which secures wider availability. A library was constructed from only 100 ng of cDNA, which enables the saving of precious specimens. Only a few steps of robust enzymatic reaction resulted in saved time, enabling more specimens to be prepared at once, and with a more reproducible size distribution among the different specimens. The obtained RNA-seq results were comparable to the conventional method. Thus, this Tn5-mediated preparation is applicable for anyone who aims to carry out deep sequencing for single cell cDNAs.

Dense genotyping identifies and localizes multiple common and rare variant association signals in celiac disease.

Using variants from the 1000 Genomes Project pilot European CEU dataset and data from additional resequencing studies, we densely genotyped 183 non-HLA risk loci previously associated with immune-mediated diseases in 12,041 individuals with celiac disease (cases) and 12,228 controls. We identified 13 new celiac disease risk loci reaching genome-wide significance, bringing the number of known loci (including the HLA locus) to 40. We found multiple independent association signals at over one-third of these loci, a finding that is attributable to a combination of common, low-frequency and rare genetic variants. Compared to previously available data such as those from HapMap3, our dense genotyping in a large sample collection provided a higher resolution of the pattern of linkage disequilibrium and suggested localization of many signals to finer scale regions. In particular, 29 of the 54 fine-mapped signals seemed to be localized to single genes and, in some instances, to gene regulatory elements. Altogether, we define the complex genetic architecture of the risk regions of and refine the risk signals for celiac disease, providing the next step toward uncovering the causal mechanisms of the disease.

Identification of type 1 diabetes-associated DNA methylation variable positions that precede disease diagnosis.

Monozygotic (MZ) twin pair discordance for childhood-onset Type 1 Diabetes (T1D) is ∼50%, implicating roles for genetic and non-genetic factors in the aetiology of this complex autoimmune disease. Although significant progress has been made in elucidating the genetics of T1D in recent years, the non-genetic component has remained poorly defined. We hypothesized that epigenetic variation could underlie some of the non-genetic component of T1D aetiology and, thus, performed an epigenome-wide association study (EWAS) for this disease. We generated genome-wide DNA methylation profiles of purified CD14+ monocytes (an immune effector cell type relevant to T1D pathogenesis) from 15 T1D-discordant MZ twin pairs. This identified 132 different CpG sites at which the direction of the intra-MZ pair DNA methylation difference significantly correlated with the diabetic state, i.e. T1D-associated methylation variable positions (T1D-MVPs). We confirmed these T1D-MVPs display statistically significant intra-MZ pair DNA methylation differences in the expected direction in an independent set of T1D-discordant MZ pairs (P = 0.035). Then, to establish the temporal origins of the T1D-MVPs, we generated two further genome-wide datasets and established that, when compared with controls, T1D-MVPs are enriched in singletons both before (P = 0.001) and at (P = 0.015) disease diagnosis, and also in singletons positive for diabetes-associated autoantibodies but disease-free even after 12 years follow-up (P = 0.0023). Combined, these results suggest that T1D-MVPs arise very early in the etiological process that leads to overt T1D. Our EWAS of T1D represents an important contribution toward understanding the etiological role of epigenetic variation in type 1 diabetes, and it is also the first systematic analysis of the temporal origins of disease-associated epigenetic variation for any human complex disease.