A phenome-wide association study of methylated GC-rich repeats identifies a GCC repeat expansion in AFF3 as a significant cause of intellectual disability.

GC-rich tandem repeat expansions (TREs) are often associated with DNA methylation, gene silencing and folate-sensitive fragile sites and underlie several congenital and late-onset disorders. Through a combination of DNA methylation profiling and tandem repeat genotyping, we identified 24 methylated TREs and investigated their effects on human traits using PheWAS in 168,641 individuals from the UK Biobank, identifying 156 significant TRE:trait associations involving 17 different TREs. Of these, a GCC expansion in the promoter of AFF3 was linked with a 2.4-fold reduced probability of completing secondary education, an effect size comparable to several recurrent pathogenic microdeletions. In a cohort of 6,371 probands with neurodevelopmental problems of suspected genetic etiology, we observed a significant enrichment of AFF3 expansions compared to controls. With a population prevalence that is at least 5-fold higher than the TRE that causes fragile X syndrome, AFF3 expansions represent a significant cause of neurodevelopmental delay.

Genetic variability in sporadic amyotrophic lateral sclerosis.

With the advent of gene therapies for amyotrophic lateral sclerosis (ALS), there is a surge in gene testing for this disease. Although there is ample experience with gene testing for C9orf72, SOD1, FUS and TARDBP in familial ALS, large studies exploring genetic variation in all ALS-associated genes in sporadic ALS (sALS) are still scarce. Gene testing in a diagnostic setting is challenging, given the complex genetic architecture of sALS, for which there are genetic variants with large and small effect sizes. Guidelines for the interpretation of genetic variants in gene panels and for counselling of patients are lacking. We aimed to provide a thorough characterization of genetic variability in ALS genes by applying the American College of Medical Genetics and Genomics (ACMG) criteria on whole genome sequencing data from a large cohort of 6013 sporadic ALS patients and 2411 matched controls from Project MinE. We studied genetic variation in 90 ALS-associated genes and applied customized ACMG-criteria to identify pathogenic and likely pathogenic variants. Variants of unknown significance were collected as well. In addition, we determined the length of repeat expansions in C9orf72, ATXN1, ATXN2 and NIPA1 using the ExpansionHunter tool. We found C9orf72 repeat expansions in 5.21% of sALS patients. In 50 ALS-associated genes, we did not identify any pathogenic or likely pathogenic variants. In 5.89%, a pathogenic or likely pathogenic variant was found, most commonly in SOD1, TARDBP, FUS, NEK1, OPTN or TBK1. Significantly more cases carried at least one pathogenic or likely pathogenic variant compared to controls (odds ratio 1.75; P-value 1.64 × 10-5). Isolated risk factors in ATXN1, ATXN2, NIPA1 and/or UNC13A were detected in 17.33% of cases. In 71.83%, we did not find any genetic clues. A combination of variants was found in 2.88%. This study provides an inventory of pathogenic and likely pathogenic genetic variation in a large cohort of sALS patients. Overall, we identified pathogenic and likely pathogenic variants in 11.13% of ALS patients in 38 known ALS genes. In line with the oligogenic hypothesis, we found significantly more combinations of variants in cases compared to controls. Many variants of unknown significance may contribute to ALS risk, but diagnostic algorithms to reliably identify and weigh them are lacking. This work can serve as a resource for counselling and for the assembly of gene panels for ALS. Further characterization of the genetic architecture of sALS is necessary given the growing interest in gene testing in ALS.

Unexpected frequency of the pathogenic AR CAG repeat expansion in the general population.

CAG repeat expansions in exon 1 of the AR gene on the X chromosome cause spinal and bulbar muscular atrophy, a male-specific progressive neuromuscular disorder associated with a variety of extra-neurological symptoms. The disease has a reported male prevalence of approximately 1:30 000 or less, but the AR repeat expansion frequency is unknown. We established a pipeline, which combines the use of the ExpansionHunter tool and visual validation, to detect AR CAG expansion on whole-genome sequencing data, benchmarked it to fragment PCR sizing, and applied it to 74 277 unrelated individuals from four large cohorts. Our pipeline showed sensitivity of 100% [95% confidence interval (CI) 90.8-100%], specificity of 99% (95% CI 94.2-99.7%), and a positive predictive value of 97.4% (95% CI 84.4-99.6%). We found the mutation frequency to be 1:3182 (95% CI 1:2309-1:4386, n = 117 734) X chromosomes-10 times more frequent than the reported disease prevalence. Modelling using the novel mutation frequency led to estimate disease prevalence of 1:6887 males, more than four times more frequent than the reported disease prevalence. This discrepancy is possibly due to underdiagnosis of this neuromuscular condition, reduced penetrance, and/or pleomorphic clinical manifestations.

Whole genome sequencing analysis reveals post-zygotic mutation variability in monozygotic twins discordant for amyotrophic lateral sclerosis.

Amyotrophic lateral sclerosis is a heterogeneous, fatal neurodegenerative disease, characterized by motor neuron loss and in 50% of cases also by cognitive and/or behavioral changes. Mendelian forms of ALS comprise approximately 10-15% of cases. The majority is however considered sporadic, but also with a high contribution of genetic risk factors. To explore the contribution of somatic mutations and/or epigenetic changes to disease risk, we performed whole genome sequencing and methylation analyses using samples from multiple tissues on a cohort of 26 monozygotic twins discordant for ALS, followed by in-depth validation and replication experiments. The results of these analyses implicate several mechanisms in ALS pathophysiology, which include a role for de novo mutations, defects in DNA damage repair and accelerated aging.

The impact of age on genetic testing decisions in amyotrophic lateral sclerosis.

Amyotrophic lateral sclerosis (ALS) is a heterogeneous neurodegenerative syndrome. In up to 20% of cases, a family history is observed. Although Mendelian disease gene variants are found in apparently sporadic ALS, genetic testing is usually restricted to those with a family history or younger patients with sporadic disease. With the advent of therapies targeting genetic ALS, it is important that everyone treatable is identified. We therefore sought to determine the probability of a clinically actionable ALS genetic test result by age of onset, globally, but using the UK as an exemplar. Blood-derived DNA was sequenced for ALS genes, and the probability of a clinically actionable genetic test result estimated. For a UK subset, age- and sex-specific population incidence rates were used to determine the number of such results missed by restricting testing by age of onset according to UK's National Genomic Test Directory criteria. There were 6274 people with sporadic ALS, 1551 from the UK. The proportion with a clinically actionable genetic test result ranged between 0.21 [95% confidence interval (CI) 0.18-0.25] in the youngest age group to 0.15 (95% CI 0.13-0.17) in the oldest age group for a full gene panel. For the UK, the equivalent proportions were 0.23 (95% CI 0.13-0.33) in the youngest age group to 0.17 (95% CI 0.13-0.21) in the oldest age group. By limiting testing in those without a family history to people with onset below 40 years, 115 of 117 (98% of all, 95% CI 96%-101%) clinically actionable test results were missed. There is a significant probability of a clinically actionable genetic test result in people with apparently sporadic ALS at all ages. Although some countries limit testing by age, doing so results in a significant number of missed pathogenic test results. Age of onset and family history should not be a barrier to genetic testing in ALS.

REViewer: haplotype-resolved visualization of read alignments in and around tandem repeats.

BACKGROUND: Expansions of short tandem repeats are the cause of many neurogenetic disorders including familial amyotrophic lateral sclerosis, Huntington disease, and many others. Multiple methods have been recently developed that can identify repeat expansions in whole genome or exome sequencing data. Despite the widely recognized need for visual assessment of variant calls in clinical settings, current computational tools lack the ability to produce such visualizations for repeat expansions. Expanded repeats are difficult to visualize because they correspond to large insertions relative to the reference genome and involve many misaligning and ambiguously aligning reads. RESULTS: We implemented REViewer, a computational method for visualization of sequencing data in genomic regions containing long repeat expansions and FlipBook, a companion image viewer designed for manual curation of large collections of REViewer images. To generate a read pileup, REViewer reconstructs local haplotype sequences and distributes reads to these haplotypes in a way that is most consistent with the fragment lengths and evenness of read coverage. To create appropriate training materials for onboarding new users, we performed a concordance study involving 12 scientists involved in short tandem repeat research. We used the results of this study to create a user guide that describes the basic principles of using REViewer as well as a guide to the typical features of read pileups that correspond to low confidence repeat genotype calls. Additionally, we demonstrated that REViewer can be used to annotate clinically relevant repeat interruptions by comparing visual assessment results of 44 FMR1 repeat alleles with the results of triplet repeat primed PCR. For 38 of these alleles, the results of visual assessment were consistent with triplet repeat primed PCR. CONCLUSIONS: Read pileup plots generated by REViewer offer an intuitive way to visualize sequencing data in regions containing long repeat expansions. Laboratories can use REViewer and FlipBook to assess the quality of repeat genotype calls as well as to visually detect interruptions or other imperfections in the repeat sequence and the surrounding flanking regions. REViewer and FlipBook are available under open-source licenses at https://github.com/illumina/REViewer and https://github.com/broadinstitute/flipbook respectively.

Genome-wide study of DNA methylation shows alterations in metabolic, inflammatory, and cholesterol pathways in ALS.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease with an estimated heritability between 40 and 50%. DNA methylation patterns can serve as proxies of (past) exposures and disease progression, as well as providing a potential mechanism that mediates genetic or environmental risk. Here, we present a blood-based epigenome-wide association study meta-analysis in 9706 samples passing stringent quality control (6763 patients, 2943 controls). We identified a total of 45 differentially methylated positions (DMPs) annotated to 42 genes, which are enriched for pathways and traits related to metabolism, cholesterol biosynthesis, and immunity. We then tested 39 DNA methylation-based proxies of putative ALS risk factors and found that high-density lipoprotein cholesterol, body mass index, white blood cell proportions, and alcohol intake were independently associated with ALS. Integration of these results with our latest genome-wide association study showed that cholesterol biosynthesis was potentially causally related to ALS. Last, DNA methylation at several DMPs and blood cell proportion estimates derived from DNA methylation data were associated with survival rate in patients, suggesting that they might represent indicators of underlying disease processes potentially amenable to therapeutic interventions.

Structural variation analysis of 6,500 whole genome sequences in amyotrophic lateral sclerosis.

There is a strong genetic contribution to Amyotrophic lateral sclerosis (ALS) risk, with heritability estimates of up to 60%. Both Mendelian and small effect variants have been identified, but in common with other conditions, such variants only explain a little of the heritability. Genomic structural variation might account for some of this otherwise unexplained heritability. We therefore investigated association between structural variation in a set of 25 ALS genes, and ALS risk and phenotype. As expected, the repeat expansion in the C9orf72 gene was identified as associated with ALS. Two other ALS-associated structural variants were identified: inversion in the VCP gene and insertion in the ERBB4 gene. All three variants were associated both with increased risk of ALS and specific phenotypic patterns of disease expression. More than 70% of people with respiratory onset ALS harboured ERBB4 insertion compared with 25% of the general population, suggesting respiratory onset ALS may be a distinct genetic subtype.

Telomere length analysis in amyotrophic lateral sclerosis using large-scale whole genome sequence data.

Background: Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease characterized by the loss of upper and lower motor neurons, leading to progressive weakness of voluntary muscles, with death following from neuromuscular respiratory failure, typically within 3 to 5 years. There is a strong genetic contribution to ALS risk. In 10% or more, a family history of ALS or frontotemporal dementia is obtained, and the Mendelian genes responsible for ALS in such families have now been identified in about 50% of cases. Only about 14% of apparently sporadic ALS is explained by known genetic variation, suggesting that other forms of genetic variation are important. Telomeres maintain DNA integrity during cellular replication, differ between sexes, and shorten naturally with age. Sex and age are risk factors for ALS and we therefore investigated telomere length in ALS. Methods: Samples were from Project MinE, an international ALS whole genome sequencing consortium that includes phenotype data. For validation we used donated brain samples from motor cortex from people with ALS and controls. Ancestry and relatedness were evaluated by principal components analysis and relationship matrices of DNA microarray data. Whole genome sequence data were from Illumina HiSeq platforms and aligned using the Isaac pipeline. TelSeq was used to quantify telomere length using whole genome sequence data. We tested the association of telomere length with ALS and ALS survival using Cox regression. Results: There were 6,580 whole genome sequences, reducing to 6,195 samples (4,315 from people with ALS and 1,880 controls) after quality control, and 159 brain samples (106 ALS, 53 controls). Accounting for age and sex, there was a 20% (95% CI 14%, 25%) increase of telomere length in people with ALS compared to controls (p = 1.1 × 10-12), validated in the brain samples (p = 0.03). Those with shorter telomeres had a 10% increase in median survival (p = 5.0×10-7). Although there was no difference in telomere length between sporadic ALS and familial ALS (p=0.64), telomere length in 334 people with ALS due to expanded C9orf72 repeats was shorter than in those without expanded C9orf72 repeats (p = 5.0×10-4). Discussion: Although telomeres shorten with age, longer telomeres are a risk factor for ALS and worsen prognosis. Longer telomeres are associated with ALS.

Common and rare variant association analyses in amyotrophic lateral sclerosis identify 15 risk loci with distinct genetic architectures and neuron-specific biology.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease with a lifetime risk of one in 350 people and an unmet need for disease-modifying therapies. We conducted a cross-ancestry genome-wide association study (GWAS) including 29,612 patients with ALS and 122,656 controls, which identified 15 risk loci. When combined with 8,953 individuals with whole-genome sequencing (6,538 patients, 2,415 controls) and a large cortex-derived expression quantitative trait locus (eQTL) dataset (MetaBrain), analyses revealed locus-specific genetic architectures in which we prioritized genes either through rare variants, short tandem repeats or regulatory effects. ALS-associated risk loci were shared with multiple traits within the neurodegenerative spectrum but with distinct enrichment patterns across brain regions and cell types. Of the environmental and lifestyle risk factors obtained from the literature, Mendelian randomization analyses indicated a causal role for high cholesterol levels. The combination of all ALS-associated signals reveals a role for perturbations in vesicle-mediated transport and autophagy and provides evidence for cell-autonomous disease initiation in glutamatergic neurons.

Chromosomal scale assembly of parasitic wasp genome reveals symbiotic virus colonization.

Endogenous viruses form an important proportion of eukaryote genomes and a source of novel functions. How large DNA viruses integrated into a genome evolve when they confer a benefit to their host, however, remains unknown. Bracoviruses are essential for the parasitism success of parasitoid wasps, into whose genomes they integrated ~103 million years ago. Here we show, from the assembly of a parasitoid wasp genome at a chromosomal scale, that bracovirus genes colonized all ten chromosomes of Cotesia congregata. Most form clusters of genes involved in particle production or parasitism success. Genomic comparison with another wasp, Microplitis demolitor, revealed that these clusters were already established ~53 mya and thus belong to remarkably stable genomic structures, the architectures of which are evolutionary constrained. Transcriptomic analyses highlight temporal synchronization of viral gene expression without resulting in immune gene induction, suggesting that no conflicts remain between ancient symbiotic partners when benefits to them converge.

Genetic analysis of ALS cases in the isolated island population of Malta.

Genetic isolates are compelling tools for mapping genes of inherited disorders. The archipelago of Malta, a sovereign microstate in the south of Europe is home to a geographically and culturally isolated population. Here, we investigate the epidemiology and genetic profile of Maltese patients with amyotrophic lateral sclerosis (ALS), identified throughout a 2-year window. Cases were largely male (66.7%) with a predominant spinal onset of symptoms (70.8%). Disease onset occurred around mid-age (median age: 64 years, men; 59.5 years, female); 12.5% had familial ALS (fALS). Annual incidence rate was 2.48 (95% CI 1.59-3.68) per 100,000 person-years. Male-to-female incidence ratio was 1.93:1. Prevalence was 3.44 (95% CI 2.01-5.52) cases per 100,000 inhabitants on 31st December 2018. Whole-genome sequencing allowed us to determine rare DNA variants that change the protein-coding sequence of ALS-associated genes. Interestingly, the Maltese ALS patient cohort was found to be negative for deleterious variants in C9orf72, SOD1, TARDBP or FUS genes, which are the most commonly mutated ALS genes globally. Nonetheless, ALS-associated repeat expansions were identified in ATXN2 and NIPA1. Variants predicted to be damaging were also detected in ALS2, DAO, DCTN1, ERBB4, SETX, SCFD1 and SPG11. A total of 40% of patients with sporadic ALS had a rare and deleterious variant or repeat expansion in an ALS-associated gene, whilst the genetic cause of two thirds of fALS cases could not be pinpointed to known ALS genes or risk loci. This warrants further studies to elucidate novel genes that cause ALS in this unique population isolate.

The Effect of SMN Gene Dosage on ALS Risk and Disease Severity.

OBJECTIVE: The role of the survival of motor neuron (SMN) gene in amyotrophic lateral sclerosis (ALS) is unclear, with several conflicting reports. A decisive result on this topic is needed, given that treatment options are available now for SMN deficiency. METHODS: In this largest multicenter case control study to evaluate the effect of SMN1 and SMN2 copy numbers in ALS, we used whole genome sequencing data from Project MinE data freeze 2. SMN copy numbers of 6,375 patients with ALS and 2,412 controls were called from whole genome sequencing data, and the reliability of the calls was tested with multiplex ligation-dependent probe amplification data. RESULTS: The copy number distribution of SMN1 and SMN2 between cases and controls did not show any statistical differences (binomial multivariate logistic regression SMN1 p = 0.54 and SMN2 p = 0.49). In addition, the copy number of SMN did not associate with patient survival (Royston-Parmar; SMN1 p = 0.78 and SMN2 p = 0.23) or age at onset (Royston-Parmar; SMN1 p = 0.75 and SMN2 p = 0.63). INTERPRETATION: In our well-powered study, there was no association of SMN1 or SMN2 copy numbers with the risk of ALS or ALS disease severity. This suggests that changing SMN protein levels in the physiological range may not modify ALS disease course. This is an important finding in the light of emerging therapies targeted at SMN deficiencies. ANN NEUROL 2021;89:686-697.

ATXN1 repeat expansions confer risk for amyotrophic lateral sclerosis and contribute to TDP-43 mislocalization.

Increasingly, repeat expansions are being identified as part of the complex genetic architecture of amyotrophic lateral sclerosis. To date, several repeat expansions have been genetically associated with the disease: intronic repeat expansions in C9orf72, polyglutamine expansions in ATXN2 and polyalanine expansions in NIPA1. Together with previously published data, the identification of an amyotrophic lateral sclerosis patient with a family history of spinocerebellar ataxia type 1, caused by polyglutamine expansions in ATXN1, suggested a similar disease association for the repeat expansion in ATXN1. We, therefore, performed a large-scale international study in 11 700 individuals, in which we showed a significant association between intermediate ATXN1 repeat expansions and amyotrophic lateral sclerosis (P = 3.33 × 10-7). Subsequent functional experiments have shown that ATXN1 reduces the nucleocytoplasmic ratio of TDP-43 and enhances amyotrophic lateral sclerosis phenotypes in Drosophila, further emphasizing the role of polyglutamine repeat expansions in the pathophysiology of amyotrophic lateral sclerosis.

ExpansionHunter Denovo: a computational method for locating known and novel repeat expansions in short-read sequencing data.

Repeat expansions are responsible for over 40 monogenic disorders, and undoubtedly more pathogenic repeat expansions remain to be discovered. Existing methods for detecting repeat expansions in short-read sequencing data require predefined repeat catalogs. Recent discoveries emphasize the need for methods that do not require pre-specified candidate repeats. To address this need, we introduce ExpansionHunter Denovo, an efficient catalog-free method for genome-wide repeat expansion detection. Analysis of real and simulated data shows that our method can identify large expansions of 41 out of 44 pathogenic repeats, including nine recently reported non-reference repeat expansions not discoverable via existing methods.

Pharmacogenetic interactions in amyotrophic lateral sclerosis: a step closer to a cure?

Genetic mutations related to amyotrophic lateral sclerosis (ALS) act through distinct pathophysiological pathways, which may lead to varying treatment responses. Here we assess the genetic interaction between C9orf72, UNC13A, and MOBP with creatine and valproic acid treatment in two clinical trials. Genotypic data was available for 309 of the 338 participants (91.4%). The UNC13A genotype affected mortality (p = 0.012), whereas C9orf72 repeat-expansion carriers exhibited a faster rate of decline in overall (p = 0.051) and bulbar functioning (p = 0.005). A dose-response pharmacogenetic interaction was identified between creatine and the A allele of the MOBP genotype (p = 0.027), suggesting a qualitative interaction in a recessive model (HR 3.96, p = 0.015). Not taking genetic information into account may mask evidence of response to treatment or be an unrecognized source of bias. Incorporating genetic data could help investigators to identify critical treatment clues in patients with ALS.

ExpansionHunter: a sequence-graph-based tool to analyze variation in short tandem repeat regions.

SUMMARY: We describe a novel computational method for genotyping repeats using sequence graphs. This method addresses the long-standing need to accurately genotype medically important loci containing repeats adjacent to other variants or imperfect DNA repeats such as polyalanine repeats. Here we introduce a new version of our repeat genotyping software, ExpansionHunter, that uses this method to perform targeted genotyping of a broad class of such loci. AVAILABILITY AND IMPLEMENTATION: ExpansionHunter is implemented in C++ and is available under the Apache License Version 2.0. The source code, documentation, and Linux/macOS binaries are available at https://github.com/Illumina/ExpansionHunter/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.