A Solve-RD ClinVar-based reanalysis of 1522 index cases from ERN-ITHACA reveals common pitfalls and misinterpretations in exome sequencing

PURPOSE: Within the Solve-RD project (https://solve-rd.eu/), the European Reference Network for Intellectual disability, TeleHealth, Autism and Congenital Anomalies aimed to investigate whether a reanalysis of exomes from unsolved cases based on ClinVar annotations could establish additional diagnoses. We present the results of the “ClinVar low-hanging fruit” reanalysis, reasons for the failure of previous analyses, and lessons learned. METHODS: Data from the first 3576 exomes (1522 probands and 2054 relatives) collected from European Reference Network for Intellectual disability, TeleHealth, Autism and Congenital Anomalies was reanalyzed by the Solve-RD consortium by evaluating for the presence of single-nucleotide variant, and small insertions and deletions already reported as (likely) pathogenic in ClinVar. Variants were filtered according to frequency, genotype, and mode of inheritance and reinterpreted. RESULTS: We identified causal variants in 59 cases (3.9%), 50 of them also raised by other approaches and 9 leading to new diagnoses, highlighting interpretation challenges: variants in genes not known to be involved in human disease at the time of the first analysis, misleading genotypes, or variants undetected by local pipelines (variants in off-target regions, low quality filters, low allelic balance, or high frequency). CONCLUSION: The “ClinVar low-hanging fruit” analysis represents an effective, fast, and easy approach to recover causal variants from exome sequencing data, herewith contributing to the reduction of the diagnostic deadlock

BOD1 Is Required for Cognitive Function in Humans and <i>Drosophila</i>

Here we report a stop-mutation in the BOD1 (Biorientation Defective 1) gene, which co-segregates with intellectual disability in a large consanguineous family, where individuals that are homozygous for the mutation have no detectable BOD1 mRNA or protein. The BOD1 protein is required for proper chromosome segregation, regulating phosphorylation of PLK1 substrates by modulating Protein Phosphatase 2A (PP2A) activity during mitosis. We report that fibroblast cell lines derived from homozygous BOD1 mutation carriers show aberrant localisation of the cell cycle kinase PLK1 and its phosphatase PP2A at mitotic kinetochores. However, in contrast to the mitotic arrest observed in BOD1-siRNA treated HeLa cells, patient-derived cells progressed through mitosis with no apparent segregation defects but at an accelerated rate compared to controls. The relatively normal cell cycle progression observed in cultured cells is in line with the absence of gross structural brain abnormalities in the affected individuals. Moreover, we found that in normal adult brain tissues BOD1 expression is maintained at considerable levels, in contrast to PLK1 expression, and provide evidence for synaptic localization of Bod1 in murine neurons. These observations suggest that BOD1 plays a cell cycle-independent role in the nervous system. To address this possibility, we established two Drosophila models, where neuron-specific knockdown of BOD1 caused pronounced learning deficits and significant abnormalities in synapse morphology. Together our results reveal novel postmitotic functions of BOD1 as well as pathogenic mechanisms that strongly support a causative role of BOD1 deficiency in the aetiology of intellectual disability. Moreover, by demonstrating its requirement for cognitive function in humans and Drosophila we provide evidence for a conserved role of BOD1 in the development and maintenance of cognitive features

Structural and non-coding variants increase the diagnostic yield of clinical whole genome sequencing for rare diseases

BACKGROUND: Whole genome sequencing is increasingly being used for the diagnosis of patients with rare diseases. However, the diagnostic yields of many studies, particularly those conducted in a healthcare setting, are often disappointingly low, at 25–30%. This is in part because although entire genomes are sequenced, analysis is often confined to in silico gene panels or coding regions of the genome. METHODS: We undertook WGS on a cohort of 122 unrelated rare disease patients and their relatives (300 genomes) who had been pre-screened by gene panels or arrays. Patients were recruited from a broad spectrum of clinical specialties. We applied a bioinformatics pipeline that would allow comprehensive analysis of all variant types. We combined established bioinformatics tools for phenotypic and genomic analysis with our novel algorithms (SVRare, ALTSPLICE and GREEN-DB) to detect and annotate structural, splice site and non-coding variants. RESULTS: Our diagnostic yield was 43/122 cases (35%), although 47/122 cases (39%) were considered solved when considering novel candidate genes with supporting functional data into account. Structural, splice site and deep intronic variants contributed to 20/47 (43%) of our solved cases. Five genes that are novel, or were novel at the time of discovery, were identified, whilst a further three genes are putative novel disease genes with evidence of causality. We identified variants of uncertain significance in a further fourteen candidate genes. The phenotypic spectrum associated with RMND1 was expanded to include polymicrogyria. Two patients with secondary findings in FBN1 and KCNQ1 were confirmed to have previously unidentified Marfan and long QT syndromes, respectively, and were referred for further clinical interventions. Clinical diagnoses were changed in six patients and treatment adjustments made for eight individuals, which for five patients was considered life-saving. CONCLUSIONS: Genome sequencing is increasingly being considered as a first-line genetic test in routine clinical settings and can make a substantial contribution to rapidly identifying a causal aetiology for many patients, shortening their diagnostic odyssey. We have demonstrated that structural, splice site and intronic variants make a significant contribution to diagnostic yield and that comprehensive analysis of the entire genome is essential to maximise the value of clinical genome sequencing

Novel diagnostic DNA methylation episignatures expand and refine the epigenetic landscapes of Mendelian disorders

Overlapping clinical phenotypes and an expanding breadth and complexity of genomic associations are a growing challenge in the diagnosis and clinical management of Mendelian disorders. The functional consequences and clinical impacts of genomic variation may involve unique, disorder-specific, genomic DNA methylation episignatures. In this study, we describe 19 novel episignature disorders and compare the findings alongside 38 previously established episignatures for a total of 57 episignatures associated with 65 genetic syndromes. We demonstrate increasing resolution and specificity ranging from protein complex, gene, sub-gene, protein domain, and even single nucleotide-level Mendelian episignatures. We show the power of multiclass modeling to develop highly accurate and disease-specific diagnostic classifiers. This study significantly expands the number and spectrum of disorders with detectable DNA methylation episignatures, improves the clinical diagnostic capabilities through the resolution of unsolved cases and the reclassification of variants of unknown clinical significance, and provides further insight into the molecular etiology of Mendelian conditions

Gain and loss of function variants in EZH1 disrupt neurogenesis and cause dominant and recessive neurodevelopmental disorders.

Genetic variants in chromatin regulators are frequently found in neurodevelopmental disorders, but their effect in disease etiology is rarely determined. Here, we uncover and functionally define pathogenic variants in the chromatin modifier EZH1 as the cause of dominant and recessive neurodevelopmental disorders in 19 individuals. EZH1 encodes one of the two alternative histone H3 lysine 27 methyltransferases of the PRC2 complex. Unlike the other PRC2 subunits, which are involved in cancers and developmental syndromes, the implication of EZH1 in human development and disease is largely unknown. Using cellular and biochemical studies, we demonstrate that recessive variants impair EZH1 expression causing loss of function effects, while dominant variants are missense mutations that affect evolutionarily conserved aminoacids, likely impacting EZH1 structure or function. Accordingly, we found increased methyltransferase activity leading to gain of function of two EZH1 missense variants. Furthermore, we show that EZH1 is necessary and sufficient for differentiation of neural progenitor cells in the developing chick embryo neural tube. Finally, using human pluripotent stem cell-derived neural cultures and forebrain organoids, we demonstrate that EZH1 variants perturb cortical neuron differentiation. Overall, our work reveals a critical role of EZH1 in neurogenesis regulation and provides molecular diagnosis for previously undefined neurodevelopmental disorders

100,000 Genomes Pilot on Rare-Disease Diagnosis in Health Care — Preliminary Report

BACKGROUND: The U.K. 100,000 Genomes Project is in the process of investigating the role of genome sequencing in patients with undiagnosed rare diseases after usual care and the alignment of this research with health care implementation in the U.K. National Health Service. Other parts of this project focus on patients with cancer and infection. METHODS: We conducted a pilot study involving 4660 participants from 2183 families, among whom 161 disorders covering a broad spectrum of rare diseases were present. We collected data on clinical features with the use of Human Phenotype Ontology terms, undertook genome sequencing, applied automated variant prioritization on the basis of applied virtual gene panels and phenotypes, and identified novel pathogenic variants through research analysis. RESULTS: Diagnostic yields varied among family structures and were highest in family trios (both parents and a proband) and families with larger pedigrees. Diagnostic yields were much higher for disorders likely to have a monogenic cause (35%) than for disorders likely to have a complex cause (11%). Diagnostic yields for intellectual disability, hearing disorders, and vision disorders ranged from 40 to 55%. We made genetic diagnoses in 25% of the probands. A total of 14% of the diagnoses were made by means of the combination of research and automated approaches, which was critical for cases in which we found etiologic noncoding, structural, and mitochondrial genome variants and coding variants poorly covered by exome sequencing. Cohortwide burden testing across 57,000 genomes enabled the discovery of three new disease genes and 19 new associations. Of the genetic diagnoses that we made, 25% had immediate ramifications for clinical decision making for the patients or their relatives. CONCLUSIONS: Our pilot study of genome sequencing in a national health care system showed an increase in diagnostic yield across a range of rare diseases. (Funded by the National Institute for Health Research and others.)

AMFR dysfunction causes autosomal recessive spastic paraplegia in human that is amenable to statin treatment in a preclinical model

Hereditary spastic paraplegias (HSP) are rare, inherited neurodegenerative or neurodevelopmental disorders that mainly present with lower limb spasticity and muscle weakness due to motor neuron dysfunction. Whole genome sequencing identified bi-allelic truncating variants in AMFR, encoding a RING-H2 finger E3 ubiquitin ligase anchored at the membrane of the endoplasmic reticulum (ER), in two previously genetically unexplained HSP-affected siblings. Subsequently, international collaboration recognized additional HSP-affected individuals with similar bi-allelic truncating AMFR variants, resulting in a cohort of 20 individuals from 8 unrelated, consanguineous families. Variants segregated with a phenotype of mainly pure but also complex HSP consisting of global developmental delay, mild intellectual disability, motor dysfunction, and progressive spasticity. Patient-derived fibroblasts, neural stem cells (NSCs), and in vivo zebrafish modeling were used to investigate pathomechanisms, including initial preclinical therapy assessment. The absence of AMFR disturbs lipid homeostasis, causing lipid droplet accumulation in NSCs and patient-derived fibroblasts which is rescued upon AMFR re-expression. Electron microscopy indicates ER morphology alterations in the absence of AMFR. Similar findings are seen in amfra-/- zebrafish larvae, in addition to altered touch-evoked escape response and defects in motor neuron branching, phenocopying the HSP observed in patients. Interestingly, administration of FDA-approved statins improves touch-evoked escape response and motor neuron branching defects in amfra-/- zebrafish larvae, suggesting potential therapeutic implications. Our genetic and functional studies identify bi-allelic truncating variants in AMFR as a cause of a novel autosomal recessive HSP by altering lipid metabolism, which may potentially be therapeutically modulated using precision medicine with statins

A new microdeletion syndrome involving TBC1D24, ATP6V0C, and PDPK1 causes epilepsy, microcephaly, and developmental delay

PURPOSE: Contiguous gene deletions are known to cause several neurodevelopmental syndromes, many of which are caused by recurrent events on chromosome 16. However, chromosomal microarray studies (CMA) still yield copy-number variants (CNVs) of unknown clinical significance. We sought to characterize eight individuals with overlapping 205-kb to 504-kb 16p13.3 microdeletions that are distinct from previously published deletion syndromes. METHODS: Clinical information on the patients and bioinformatic scores for the deleted genes were analyzed. RESULTS: All individuals in our cohort displayed developmental delay, intellectual disability, and various forms of seizures. Six individuals were microcephalic and two had strabismus. The deletion was absent in all 13 parents who were available for testing. The area of overlap encompasses seven genes including TBC1D24, ATP6V0C, and PDPK1 (also known as PDK1). Bi-allelic TBC1D24 pathogenic variants are known to cause nonsyndromic deafness, epileptic disorders, or DOORS syndrome (deafness, onychodystrophy, osteodystrophy, mental retardation, seizures). Sanger sequencing of the nondeleted TBC1D24 allele did not yield any additional pathogenic variants. CONCLUSIONS: We propose that 16p13.3 microdeletions resulting in simultaneous haploinsufficiencies of TBC1D24, ATP6V0C, and PDPK1 cause a novel rare contiguous gene deletion syndrome of microcephaly, developmental delay, intellectual disability, and epilepsy

Widespread Contribution of Gdf7 Lineage to Cerebellar Cell Types and Implications for Hedgehog-Driven Medulloblastoma Formation

The roof plate is a specialized embryonic midline tissue of the central nervous system that functions as a signaling center regulating dorsal neural patterning. In the developing hindbrain, roof plate cells express Gdf7 and previous genetic fate mapping studies showed that these cells contribute mostly to non-neural choroid plexus epithelium. We demonstrate here that constitutive activation of the Sonic hedgehog signaling pathway in the Gdf7 lineage invariably leads to medulloblastoma. Lineage tracing analysis reveals that Gdf7-lineage cells not only are a source of choroid plexus epithelial cells, but are also present in the cerebellar rhombic lip and contribute to a subset of cerebellar granule neuron precursors, the presumed cell-of-origin for Sonic hedgehog-driven medulloblastoma. We further show that Gdf7-lineage cells also contribute to multiple neuronal and glial cell types in the cerebellum, including glutamatergic granule neurons, unipolar brush cells, Purkinje neurons, GABAergic interneurons, Bergmann glial cells, and white matter astrocytes. These findings establish hindbrain roof plate as a novel source of diverse neural cell types in the cerebellum that is also susceptible to oncogenic transformation by deregulated Sonic hedgehog signaling

Mediator Kinase Inhibition Further Activates Super-Enhancer Associated Genes in AML

Super-enhancers (SEs), which are composed of large clusters of enhancers densely loaded with the Mediator complex, transcription factors (TFs), and chromatin regulators, drive high expression of genes implicated in cell identity and disease, such as lineage-controlling TFs and oncogenes 1, 2. BRD4 and CDK7 are positive regulators of SE-mediated transcription3,4,5. In contrast, negative regulators of SE-associated genes have not been well described. Here we report that Mediator-associated kinases cyclin-dependent kinase 8 (CDK8) and CDK19 restrain increased activation of key SE-associated genes in acute myeloid leukaemia (AML) cells. We determined that the natural product cortistatin A (CA) selectively inhibited Mediator kinases, had antileukaemic activity in vitro and in vivo, and disproportionately induced upregulation of SE-associated genes in CA-sensitive AML cell lines but not in CA-insensitive cell lines. In AML cells, CA upregulated SE-associated genes with tumour suppressor and lineage-controlling functions, including the TFs CEBPA, IRF8, IRF1 and ETV6 6, 7, 8. The BRD4 inhibitor I-BET151 downregulated these SE-associated genes, yet also has antileukaemic activity. Individually increasing or decreasing expression of these TFs suppressed AML cell growth, providing evidence that leukaemia cells are sensitive to dosage of SE-associated genes. Our results demonstrate that Mediator kinases can negatively regulate SE-associated gene expression in specific cell types and can be pharmacologically targeted as a therapeutic approach to AML