Exome copy number variant detection, analysis, and classification in a large cohort of families with undiagnosed rare genetic disease.

Copy number variants (CNVs) are significant contributors to the pathogenicity of rare genetic diseases and, with new innovative methods, can now reliably be identified from exome sequencing. Challenges still remain in accurate classification of CNV pathogenicity. CNV calling using GATK-gCNV was performed on exomes from a cohort of 6,633 families (15,759 individuals) with heterogeneous phenotypes and variable prior genetic testing collected at the Broad Institute Center for Mendelian Genomics of the Genomics Research to Elucidate the Genetics of Rare Diseases consortium and analyzed using the seqr platform. The addition of CNV detection to exome analysis identified causal CNVs for 171 families (2.6%). The estimated sizes of CNVs ranged from 293 bp to 80 Mb. The causal CNVs consisted of 140 deletions, 15 duplications, 3 suspected complex structural variants (SVs), 3 insertions, and 10 complex SVs, the latter two groups being identified by orthogonal confirmation methods. To classify CNV variant pathogenicity, we used the 2020 American College of Medical Genetics and Genomics/ClinGen CNV interpretation standards and developed additional criteria to evaluate allelic and functional data as well as variants on the X chromosome to further advance the framework. We interpreted 151 CNVs as likely pathogenic/pathogenic and 20 CNVs as high-interest variants of uncertain significance. Calling CNVs from existing exome data increases the diagnostic yield for individuals undiagnosed after standard testing approaches, providing a higher-resolution alternative to arrays at a fraction of the cost of genome sequencing. Our improvements to the classification approach advances the systematic framework to assess the pathogenicity of CNVs.

Neurodevelopmental disorders associated variants in ADAT3 disrupt the activity of the ADAT2/ADAT3 tRNA deaminase complex and impair neuronal migration.

The ADAT2/ADAT3 complex catalyzes the adenosine to inosine modification at the wobble position of eukaryotic tRNAs. Mutations in ADAT3 , the catalytically inactive subunit of the ADAT2/ADAT3 complex, have been identified in patients presenting with severe neurodevelopmental disorders (NDDs). Yet, the physiological function of ADAT2/ADAT3 complex during brain development remains totally unknown. Here we showed that maintaining a proper level of ADAT2/ADAT3 catalytic activity is required for correct radial migration of projection neurons in the developing mouse cortex. In addition, we not only reported 7 new NDD patients carrying biallelic variants in ADAT3 but also deeply characterize the impact of those variants on ADAT2/ADAT3 structure, biochemical properties, enzymatic activity and tRNAs editing and abundance. We demonstrated that all the identified variants alter both the expression and the activity of the complex leading to a significant decrease of I 34 with direct consequence on their steady-state. Using in vivo complementation assays, we correlated the severity of the migration phenotype with the degree of the loss of function caused by the variants. Altogether, our results indicate a critical role of ADAT2/ADAT3 during cortical development and provide cellular and molecular insights into the pathogenicity of ADAT3-related neurodevelopmental disorder.

Somatic Mosaicism in Amyotrophic Lateral Sclerosis and Frontotemporal Dementia Reveals Widespread Degeneration from Focal Mutations.

Although mutations in dozens of genes have been implicated in familial forms of amyotrophic lateral sclerosis (fALS) and frontotemporal degeneration (fFTD), most cases of these conditions are sporadic (sALS and sFTD), with no family history, and their etiology remains obscure. We tested the hypothesis that somatic mosaic mutations, present in some but not all cells, might contribute in these cases, by performing ultra-deep, targeted sequencing of 88 genes associated with neurodegenerative diseases in postmortem brain and spinal cord samples from 404 individuals with sALS or sFTD and 144 controls. Known pathogenic germline mutations were found in 20.6% of ALS, and 26.5% of FTD cases. Predicted pathogenic somatic mutations in ALS/FTD genes were observed in 2.7% of sALS and sFTD cases that did not carry known pathogenic or novel germline mutations. Somatic mutations showed low variant allele fraction (typically <2%) and were often restricted to the region of initial discovery, preventing detection through genetic screening in peripheral tissues. Damaging somatic mutations were preferentially enriched in primary motor cortex of sALS and prefrontal cortex of sFTD, mirroring regions most severely affected in each disease. Somatic mutation analysis of bulk RNA-seq data from brain and spinal cord from an additional 143 sALS cases and 23 controls confirmed an overall enrichment of somatic mutations in sALS. Two adult sALS cases were identified bearing pathogenic somatic mutations in DYNC1H1 and LMNA, two genes associated with pediatric motor neuron degeneration. Our study suggests that somatic mutations in fALS/fFTD genes, and in genes associated with more severe diseases in the germline state, contribute to sALS and sFTD, and that mosaic mutations in a small fraction of cells in focal regions of the nervous system can ultimately result in widespread degeneration.

Exome copy number variant detection, analysis and classification in a large cohort of families with undiagnosed rare genetic disease.

Copy number variants (CNVs) are significant contributors to the pathogenicity of rare genetic diseases and with new innovative methods can now reliably be identified from exome sequencing. Challenges still remain in accurate classification of CNV pathogenicity. CNV calling using GATK-gCNV was performed on exomes from a cohort of 6,633 families (15,759 individuals) with heterogeneous phenotypes and variable prior genetic testing collected at the Broad Institute Center for Mendelian Genomics of the GREGoR consortium. Each family's CNV data was analyzed using the seqr platform and candidate CNVs classified using the 2020 ACMG/ClinGen CNV interpretation standards. We developed additional evidence criteria to address situations not covered by the current standards. The addition of CNV calling to exome analysis identified causal CNVs for 173 families (2.6%). The estimated sizes of CNVs ranged from 293 bp to 80 Mb with estimates that 44% would not have been detected by standard chromosomal microarrays. The causal CNVs consisted of 141 deletions, 15 duplications, 4 suspected complex structural variants (SVs), 3 insertions and 10 complex SVs, the latter two groups being identified by orthogonal validation methods. We interpreted 153 CNVs as likely pathogenic/pathogenic and 20 CNVs as high interest variants of uncertain significance. Calling CNVs from existing exome data increases the diagnostic yield for individuals undiagnosed after standard testing approaches, providing a higher resolution alternative to arrays at a fraction of the cost of genome sequencing. Our improvements to the classification approach advances the systematic framework to assess the pathogenicity of CNVs.

Exome Sequencing and the Identification of New Genes and Shared Mechanisms in Polymicrogyria.

Importance: Polymicrogyria is the most commonly diagnosed cortical malformation and is associated with neurodevelopmental sequelae including epilepsy, motor abnormalities, and cognitive deficits. Polymicrogyria frequently co-occurs with other brain malformations or as part of syndromic diseases. Past studies of polymicrogyria have defined heterogeneous genetic and nongenetic causes but have explained only a small fraction of cases. Objective: To survey germline genetic causes of polymicrogyria in a large cohort and to consider novel polymicrogyria gene associations. Design, Setting, and Participants: This genetic association study analyzed panel sequencing and exome sequencing of accrued DNA samples from a retrospective cohort of families with members with polymicrogyria. Samples were accrued over more than 20 years (1994 to 2020), and sequencing occurred in 2 stages: panel sequencing (June 2015 to January 2016) and whole-exome sequencing (September 2019 to March 2020). Individuals seen at multiple clinical sites for neurological complaints found to have polymicrogyria on neuroimaging, then referred to the research team by evaluating clinicians, were included in the study. Targeted next-generation sequencing and/or exome sequencing were performed on probands (and available parents and siblings) from 284 families with individuals who had isolated polymicrogyria or polymicrogyria as part of a clinical syndrome and no genetic diagnosis at time of referral from clinic, with sequencing from 275 families passing quality control. Main Outcomes and Measures: The number of families in whom genetic sequencing yielded a molecular diagnosis that explained the polymicrogyria in the family. Secondarily, the relative frequency of different genetic causes of polymicrogyria and whether specific genetic causes were associated with co-occurring head size changes were also analyzed. Results: In 32.7% (90 of 275) of polymicrogyria-affected families, genetic variants were identified that provided satisfactory molecular explanations. Known genes most frequently implicated by polymicrogyria-associated variants in this cohort were PIK3R2, TUBB2B, COL4A1, and SCN3A. Six candidate novel polymicrogyria genes were identified or confirmed: de novo missense variants in PANX1, QRICH1, and SCN2A and compound heterozygous variants in TMEM161B, KIF26A, and MAN2C1, each with consistent genotype-phenotype relationships in multiple families. Conclusions and Relevance: This study's findings reveal a higher than previously recognized rate of identifiable genetic causes, specifically of channelopathies, in individuals with polymicrogyria and support the utility of exome sequencing for families affected with polymicrogyria.

Bi-allelic variants in INTS11 are associated with a complex neurological disorder.

The Integrator complex is a multi-subunit protein complex that regulates the processing of nascent RNAs transcribed by RNA polymerase II (RNAPII), including small nuclear RNAs, enhancer RNAs, telomeric RNAs, viral RNAs, and protein-coding mRNAs. Integrator subunit 11 (INTS11) is the catalytic subunit that cleaves nascent RNAs, but, to date, mutations in this subunit have not been linked to human disease. Here, we describe 15 individuals from 10 unrelated families with bi-allelic variants in INTS11 who present with global developmental and language delay, intellectual disability, impaired motor development, and brain atrophy. Consistent with human observations, we find that the fly ortholog of INTS11, dIntS11, is essential and expressed in the central nervous systems in a subset of neurons and most glia in larval and adult stages. Using Drosophila as a model, we investigated the effect of seven variants. We found that two (p.Arg17Leu and p.His414Tyr) fail to rescue the lethality of null mutants, indicating that they are strong loss-of-function variants. Furthermore, we found that five variants (p.Gly55Ser, p.Leu138Phe, p.Lys396Glu, p.Val517Met, and p.Ile553Glu) rescue lethality but cause a shortened lifespan and bang sensitivity and affect locomotor activity, indicating that they are partial loss-of-function variants. Altogether, our results provide compelling evidence that integrity of the Integrator RNA endonuclease is critical for brain development.

TMEM161B regulates cerebral cortical gyration, Sonic Hedgehog signaling, and ciliary structure in the developing central nervous system.

Sonic hedgehog signaling regulates processes of embryonic development across multiple tissues, yet factors regulating context-specific Shh signaling remain poorly understood. Exome sequencing of families with polymicrogyria (disordered cortical folding) revealed multiple individuals with biallelic deleterious variants in TMEM161B, which encodes a multi-pass transmembrane protein of unknown function. Tmem161b null mice demonstrated holoprosencephaly, craniofacial midline defects, eye defects, and spinal cord patterning changes consistent with impaired Shh signaling, but were without limb defects, suggesting a CNS-specific role of Tmem161b. Tmem161b depletion impaired the response to Smoothened activation in vitro and disrupted cortical histogenesis in vivo in both mouse and ferret models, including leading to abnormal gyration in the ferret model. Tmem161b localizes non-exclusively to the primary cilium, and scanning electron microscopy revealed shortened, dysmorphic, and ballooned ventricular zone cilia in the Tmem161b null mouse, suggesting that the Shh-related phenotypes may reflect ciliary dysfunction. Our data identify TMEM161B as a regulator of cerebral cortical gyration, as involved in primary ciliary structure, as a regulator of Shh signaling, and further implicate Shh signaling in human gyral development.

Missense MED12 variants in 22 males with intellectual disability: From nonspecific symptoms to complete syndromes.

We describe the phenotype of 22 male patients (20 probands) carrying a hemizygous missense variant in MED12. The phenotypic spectrum is very broad ranging from nonspecific intellectual disability (ID) to the three well-known syndromes: Opitz-Kaveggia syndrome, Lujan-Fryns syndrome, or Ohdo syndrome. The identified variants were randomly distributed throughout the gene (p = 0.993, χ2 test), but mostly outside the functional domains (p = 0.004; χ2 test). Statistical analyses did not show a correlation between the MED12-related phenotypes and the locations of the variants (p = 0.295; Pearson correlation), nor the protein domain involved (p = 0.422; Pearson correlation). In conclusion, establishing a genotype-phenotype correlation in MED12-related diseases remains challenging. Therefore, we think that patients with a causative MED12 variant are currently underdiagnosed due to the broad patients' clinical presentations.

ANKLE2-related microcephaly: A variable microcephaly syndrome resembling Zika infection.

OBJECTIVE: This study delineates the clinical and molecular spectrum of ANKLE2-related microcephaly (MIC), as well as highlights shared pathological mechanisms between ANKLE2 and the Zika virus. METHODS: We identified 12 individuals with MIC and variants in ANKLE2 with a broad range of features. Probands underwent thorough phenotypic evaluations, developmental assessments, and anthropometric measurements. Brain imaging studies were systematically reviewed for developmental abnormalities. We functionally interrogated a subset of identified ANKLE2 variants in Drosophila melanogaster. RESULTS: All individuals had MIC (z-score ≤ -3), including nine with congenital MIC. We identified a broad range of brain abnormalities including simplified cortical gyral pattern, full or partial callosal agenesis, increased extra-axial spaces, hypomyelination, cerebellar vermis hypoplasia, and enlarged cisterna magna. All probands had developmental delays in at least one domain, with speech and language delays being the most common. Six probands had skin findings characteristic of ANKLE2 including hyper- and hypopigmented macules. Only one individual had scalp rugae. Functional characterization in Drosophila recapitulated the human MIC phenotype. Of the four variants tested, p.Val229Gly, p.Arg236*, and p.Arg536Cys acted as partial-loss-of-function variants, whereas the c.1421-1G>C splicing variant demonstrated a strong loss-of-function effect. INTERPRETATION: Deleterious variants in the ANKLE2 gene cause a unique MIC syndrome characterized by congenital or postnatal MIC, a broad range of structural brain abnormalities, and skin pigmentary changes. Thorough functional characterization has identified shared pathogenic mechanisms between ANKLE2-related MIC and congenital Zika virus infection. This study further highlights the importance of a thorough diagnostic evaluation including molecular diagnostic testing in individuals with MIC.

Somatic genomic changes in single Alzheimer's disease neurons.

Dementia in Alzheimer's disease progresses alongside neurodegeneration1-4, but the specific events that cause neuronal dysfunction and death remain poorly understood. During normal ageing, neurons progressively accumulate somatic mutations5 at rates similar to those of dividing cells6,7 which suggests that genetic factors, environmental exposures or disease states might influence this accumulation5. Here we analysed single-cell whole-genome sequencing data from 319 neurons from the prefrontal cortex and hippocampus of individuals with Alzheimer's disease and neurotypical control individuals. We found that somatic DNA alterations increase in individuals with Alzheimer's disease, with distinct molecular patterns. Normal neurons accumulate mutations primarily in an age-related pattern (signature A), which closely resembles 'clock-like' mutational signatures that have been previously described in healthy and cancerous cells6-10. In neurons affected by Alzheimer's disease, additional DNA alterations are driven by distinct processes (signature C) that highlight C>A and other specific nucleotide changes. These changes potentially implicate nucleotide oxidation4,11, which we show is increased in Alzheimer's-disease-affected neurons in situ. Expressed genes exhibit signature-specific damage, and mutations show a transcriptional strand bias, which suggests that transcription-coupled nucleotide excision repair has a role in the generation of mutations. The alterations in Alzheimer's disease affect coding exons and are predicted to create dysfunctional genetic knockout cells and proteostatic stress. Our results suggest that known pathogenic mechanisms in Alzheimer's disease may lead to genomic damage to neurons that can progressively impair function. The aberrant accumulation of DNA alterations in neurodegeneration provides insight into the cascade of molecular and cellular events that occurs in the development of Alzheimer's disease.

Variants in ADD1 cause intellectual disability, corpus callosum dysgenesis, and ventriculomegaly in humans.

PURPOSE: Adducins interconnect spectrin and actin filaments to form polygonal scaffolds beneath the cell membranes and form ring-like structures in neuronal axons. Adducins regulate mouse neural development, but their function in the human brain is unknown. METHODS: We used exome sequencing to uncover ADD1 variants associated with intellectual disability (ID) and brain malformations. We studied ADD1 splice isoforms in mouse and human neocortex development with RNA sequencing, super resolution imaging, and immunoblotting. We investigated 4 variant ADD1 proteins and heterozygous ADD1 cells for protein expression and ADD1-ADD2 dimerization. We studied Add1 functions in vivo using Add1 knockout mice. RESULTS: We uncovered loss-of-function ADD1 variants in 4 unrelated individuals affected by ID and/or structural brain defects. Three additional de novo copy number variations covering the ADD1 locus were associated with ID and brain malformations. ADD1 is highly expressed in the neocortex and the corpus callosum, whereas ADD1 splice isoforms are dynamically expressed between cortical progenitors and postmitotic neurons. Human variants impair ADD1 protein expression and/or dimerization with ADD2. Add1 knockout mice recapitulate corpus callosum dysgenesis and ventriculomegaly phenotypes. CONCLUSION: Our human and mouse genetics results indicate that pathogenic ADD1 variants cause corpus callosum dysgenesis, ventriculomegaly, and/or ID.

Biallelic loss-of-function variants in WDR11 are associated with microcephaly and intellectual disability.

Heterozygous missense variants in the WD repeat domain 11 (WDR11) gene are associated with hypogonadotropic hypogonadism in humans. In contrast, knockout of both alleles of Wdr11 in mice results in a more severe phenotype with growth and developmental delay, features of holoprosencephaly, heart defects and reproductive disorders. Similar developmental defects known to be associated with aberrant hedgehog signaling and ciliogenesis have been found in zebrafish after Wdr11 knockdown. We here report biallelic loss-of-function variants in the WDR11 gene in six patients from three independent families with intellectual disability, microcephaly and short stature. The findings suggest that biallelic WDR11 variants in humans result in an overlapping but milder phenotype compared to Wdr11-deficient animals. However, the observed human phenotype differs significantly from dominantly inherited variants leading to hypogonadotropic hypogonadism, suggesting that recessive WDR11 variants result in a clinically distinct entity.

Early role for a Na+,K+-ATPase (ATP1A3) in brain development.

Osmotic equilibrium and membrane potential in animal cells depend on concentration gradients of sodium (Na+) and potassium (K+) ions across the plasma membrane, a function catalyzed by the Na+,K+-ATPase α-subunit. Here, we describe ATP1A3 variants encoding dysfunctional α3-subunits in children affected by polymicrogyria, a developmental malformation of the cerebral cortex characterized by abnormal folding and laminar organization. To gain cell-biological insights into the spatiotemporal dynamics of prenatal ATP1A3 expression, we built an ATP1A3 transcriptional atlas of fetal cortical development using mRNA in situ hybridization and transcriptomic profiling of ∼125,000 individual cells with single-cell RNA sequencing (Drop-seq) from 11 areas of the midgestational human neocortex. We found that fetal expression of ATP1A3 is most abundant to a subset of excitatory neurons carrying transcriptional signatures of the developing subplate, yet also maintains expression in nonneuronal cell populations. Moving forward a year in human development, we profiled ∼52,000 nuclei from four areas of an infant neocortex and show that ATP1A3 expression persists throughout early postnatal development, most predominantly in inhibitory neurons, including parvalbumin interneurons in the frontal cortex. Finally, we discovered the heteromeric Na+,K+-ATPase pump complex may form nonredundant cell-type-specific α-ß isoform combinations, including α3-ß1 in excitatory neurons and α3-ß2 in inhibitory neurons. Together, the developmental malformation phenotype of affected individuals and single-cell ATP1A3 expression patterns point to a key role for α3 in human cortex development, as well as a cell-type basis for pre- and postnatal ATP1A3-associated diseases.

A recurrent, homozygous EMC10 frameshift variant is associated with a syndrome of developmental delay with variable seizures and dysmorphic features.

PURPOSE: The endoplasmic reticulum membrane complex (EMC) is a highly conserved, multifunctional 10-protein complex related to membrane protein biology. In seven families, we identified 13 individuals with highly overlapping phenotypes who harbor a single identical homozygous frameshift variant in EMC10. METHODS: Using exome, genome, and Sanger sequencing, a recurrent frameshift EMC10 variant was identified in affected individuals in an international cohort of consanguineous families. Multiple families were independently identified and connected via Matchmaker Exchange and internal databases. We assessed the effect of the frameshift variant on EMC10 RNA and protein expression and evaluated EMC10 expression in normal human brain tissue using immunohistochemistry. RESULTS: A homozygous variant EMC10 c.287delG (Refseq NM_206538.3, p.Gly96Alafs*9) segregated with affected individuals in each family, who exhibited a phenotypic spectrum of intellectual disability (ID) and global developmental delay (GDD), variable seizures and variable dysmorphic features (elongated face, curly hair, cubitus valgus, and arachnodactyly). The variant arose on two founder haplotypes and results in significantly reduced EMC10 RNA expression and an unstable truncated EMC10 protein. CONCLUSION: We propose that a homozygous loss-of-function variant in EMC10 causes a novel syndromic neurodevelopmental phenotype. Remarkably, the recurrent variant is likely the result of a hypermutable site and arose on distinct founder haplotypes.