Mutations in the tail and rod domains of the neurofilament heavy-chain gene increase the risk of ALS.

OBJECTIVE: Neurofilament heavy-chain gene (NEFH) variants are associated with multiple neurodegenerative diseases, however, their relationship with ALS has not been robustly explored. Still, NEFH is commonly included in genetic screening panels worldwide. We therefore aimed to determine if NEFH variants modify ALS risk. METHODS: Genetic data of 11,130 people with ALS and 7,416 controls from the literature and Project MinE were analysed. We performed meta-analyses of published case-control studies reporting NEFH variants, and variant analysis of NEFH in Project MinE whole-genome sequencing data. RESULTS: Fixed-effects meta-analysis found that rare (MAF <1%) missense variants in the tail domain of NEFH increase ALS risk (OR 4.55, 95% CI 2.13-9.71, p < 0.0001). In Project MinE, ultrarare NEFH variants increased ALS risk (OR 1.37 95% CI 1.14-1.63, p = 0.0007), with rod domain variants (mostly intronic) appearing to drive the association (OR 1.45 95% CI 1.18-1.77, pMadsen-Browning = 0.0007, pSKAT-O = 0.003). While in the tail domain, ultrarare (MAF <0.1%) pathogenic missense variants were also associated with higher risk of ALS (OR 1.94, 95% CI 0.86-4.37, pMadsen-Browning = 0.039), supporting the meta-analysis results. Finally, several tail in-frame deletions were also found to affect disease risk, however, both protective and pathogenic deletions were found in this domain, highlighting an intricate architecture that requires further investigation. INTERPRETATION: We showed that NEFH tail missense and in-frame deletion variants, and intronic rod variants are risk factors for ALS. However, they are not variants of large effect, and their functional impact needs to be clarified in further studies. Therefore, their inclusion in routine genetic screening panels should be reconsidered.

Rare and common genetic determinants of mitochondrial function determine severity but not risk of amyotrophic lateral sclerosis.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease involving selective vulnerability of energy-intensive motor neurons (MNs). It has been unclear whether mitochondrial function is an upstream driver or a downstream modifier of neurotoxicity. We separated upstream genetic determinants of mitochondrial function, including genetic variation within the mitochondrial genome or autosomes; from downstream changeable factors including mitochondrial DNA copy number (mtCN). Across three cohorts including 6,437 ALS patients, we discovered that a set of mitochondrial haplotypes, chosen because they are linked to measurements of mitochondrial function, are a determinant of ALS survival following disease onset, but do not modify ALS risk. One particular haplotype appeared to be neuroprotective and was significantly over-represented in two cohorts of long-surviving ALS patients. Causal inference for mitochondrial function was achievable using mitochondrial haplotypes, but not autosomal SNPs in traditional Mendelian randomization (MR). Furthermore, rare loss-of-function genetic variants within, and reduced MN expression of, ACADM and DNA2 lead to ∼50 % shorter ALS survival; both proteins are implicated in mitochondrial function. Both mtCN and cellular vulnerability are linked to DNA2 function in ALS patient-derived neurons. Finally, MtCN responds dynamically to the onset of ALS independently of mitochondrial haplotype, and is correlated with disease severity. We conclude that, based on the genetic measures we have employed, mitochondrial function is a therapeutic target for amelioration of disease severity but not prevention of ALS.

Co-targeting HSP90 alpha and CDK7 overcomes resistance against HSP90 inhibitors in BCR-ABL1+ leukemia cells.

HSP90 has emerged as an appealing anti-cancer target. However, HSP90 inhibitors (HSP90i) are characterized by limited clinical utility, primarily due to the resistance acquisition via heat shock response (HSR) induction. Understanding the roles of abundantly expressed cytosolic HSP90 isoforms (α and ß) in sustaining malignant cells' growth and the mechanisms of resistance to HSP90i is crucial for exploiting their clinical potential. Utilizing multi-omics approaches, we identified that ablation of the HSP90ß isoform induces the overexpression of HSP90α and extracellular-secreted HSP90α (eHSP90α). Notably, we found that the absence of HSP90α causes downregulation of PTPRC (or CD45) expression and restricts in vivo growth of BCR-ABL1+ leukemia cells. Subsequently, chronic long-term exposure to the clinically advanced HSP90i PU-H71 (Zelavespib) led to copy number gain and mutation (p.S164F) of the HSP90AA1 gene, and HSP90α overexpression. In contrast, acquired resistance toward other tested HSP90i (Tanespimycin and Coumermycin A1) was attained by MDR1 efflux pump overexpression. Remarkably, combined CDK7 and HSP90 inhibition display synergistic activity against therapy-resistant BCR-ABL1+ patient leukemia cells via blocking pro-survival HSR and HSP90α overexpression, providing a novel strategy to avoid the emergence of resistance against treatment with HSP90i alone.

Unsupervised machine learning identifies distinct ALS molecular subtypes in post-mortem motor cortex and blood expression data.

Amyotrophic lateral sclerosis (ALS) displays considerable clinical and genetic heterogeneity. Machine learning approaches have previously been utilised for patient stratification in ALS as they can disentangle complex disease landscapes. However, lack of independent validation in different populations and tissue samples have greatly limited their use in clinical and research settings. We overcame these issues by performing hierarchical clustering on the 5000 most variably expressed autosomal genes from motor cortex expression data of people with sporadic ALS from the KCL BrainBank (N = 112). Three molecular phenotypes linked to ALS pathogenesis were identified: synaptic and neuropeptide signalling, oxidative stress and apoptosis, and neuroinflammation. Cluster validation was achieved by applying linear discriminant analysis models to cases from TargetALS US motor cortex (N = 93), as well as Italian (N = 15) and Dutch (N = 397) blood expression datasets, for which there was a high assignment probability (80-90%) for each molecular subtype. The ALS and motor cortex specificity of the expression signatures were tested by mapping KCL BrainBank controls (N = 59), and occipital cortex (N = 45) and cerebellum (N = 123) samples from TargetALS to each cluster, before constructing case-control and motor cortex-region logistic regression classifiers. We found that the signatures were not only able to distinguish people with ALS from controls (AUC 0.88 ± 0.10), but also reflect the motor cortex-based disease process, as there was perfect discrimination between motor cortex and the other brain regions. Cell types known to be involved in the biological processes of each molecular phenotype were found in higher proportions, reinforcing their biological interpretation. Phenotype analysis revealed distinct cluster-related outcomes in both motor cortex datasets, relating to disease onset and progression-related measures. Our results support the hypothesis that different mechanisms underpin ALS pathogenesis in subgroups of patients and demonstrate potential for the development of personalised treatment approaches. Our method is available for the scientific and clinical community at https://alsgeclustering.er.kcl.ac.uk .

Molecular dynamics analysis of superoxide dismutase 1 mutations suggests decoupling between mechanisms underlying ALS onset and progression.

Mutations in the superoxide dismutase 1 (SOD1) gene are the second most common known cause of ALS. SOD1 variants express high phenotypic variability and over 200 have been reported in people with ALS. It was previously proposed that variants can be broadly classified in two groups, 'wild-type like' (WTL) and 'metal binding region' (MBR) variants, based on their structural location and biophysical properties. MBR variants, but not WTL variants, were associated with a reduction of SOD1 enzymatic activity. In this study we used molecular dynamics and large clinical datasets to characterise the differences in the structural and dynamic behaviour of WTL and MBR variants with respect to the wild-type SOD1, and how such differences influence the ALS clinical phenotype. Our study identified marked structural differences, some of which are observed in both variant groups, while others are group specific. Moreover, collecting clinical data of approximately 500 SOD1 ALS patients carrying variants, we showed that the survival time of patients carrying an MBR variant is generally longer (∼6 years median difference, p < 0.001) with respect to patients with a WTL variant. In conclusion, our study highlighted key differences in the dynamic behaviour between WTL and MBR SOD1 variants, and between variants and wild-type SOD1 at an atomic and molecular level, that could be further investigated to explain the associated phenotypic variability. Our results support the hypothesis of a decoupling between mechanisms of onset and progression of SOD1 ALS, and an involvement of loss-of-function of SOD1 with the disease progression.

VariantSurvival: a tool to identify genotype-treatment response.

Motivation: For a number of neurological diseases, such as Alzheimer's disease, amyotrophic lateral sclerosis, and many others, certain genes are known to be involved in the disease mechanism. A common question is whether a structural variant in any such gene may be related to drug response in clinical trials and how this relationship can contribute to the lifecycle of drug development. Results: To this end, we introduce VariantSurvival, a tool that identifies changes in survival relative to structural variants within target genes. VariantSurvival matches annotated structural variants with genes that are clinically relevant to neurological diseases. A Cox regression model determines the change in survival between the placebo and clinical trial groups with respect to the number of structural variants in the drug target genes. We demonstrate the functionality of our approach with the exemplary case of the SETX gene. VariantSurvival has a user-friendly and lightweight graphical user interface built on the shiny web application package.

SOD1-ALS-Browser: a web-utility for investigating the clinical phenotype in SOD1 amyotrophic lateral sclerosis.

Objective: Variants in the superoxide dismutase (SOD1) gene are among the most common genetic causes of amyotrophic lateral sclerosis. Reflecting the wide spectrum of putatively deleterious variants that have been reported to date, it has become clear that SOD1-linked ALS presents a highly variable age at symptom onset and disease duration.Methods: Here we describe an open access web tool for comparative phenotype analysis in ALS: https://sod1-als-browser.rosalind.kcl.ac.uk/. The tool contains a built-in dataset of clinical information from 1383 people with ALS harboring a SOD1 variant resulting in one of 162 unique amino acid sequence alterations and from a non-SOD1 comparator ALS cohort of 13,469 individuals. We present two examples of analyses possible with this tool, testing how the ALS phenotype relates to SOD1 variants that alter amino acid residue hydrophobicity and to distinct variants at the 94th residue of SOD1, where six are sampled.Results and conclusions: The tool provides immediate access to the datasets and enables bespoke analysis of phenotypic trends associated with different protein variants, including the option for users to upload their own datasets for integration with the server data. The tool can be used to study SOD1-ALS and provides an analytical framework to study the differences between other user-uploaded ALS groups and our large reference database of SOD1 and non-SOD1 ALS. The tool is designed to be useful for clinicians and researchers, including those without programming expertise, and is highly flexible in the analyses that can be conducted.

DNAscan2: a versatile, scalable, and user-friendly analysis pipeline for human next-generation sequencing data.

SUMMARY: The current widespread adoption of next-generation sequencing (NGS) in all branches of basic research and clinical genetics fields means that users with highly variable informatics skills, computing facilities and application purposes need to process, analyse, and interpret NGS data. In this landscape, versatility, scalability, and user-friendliness are key characteristics for an NGS analysis software. We developed DNAscan2, a highly flexible, end-to-end pipeline for the analysis of NGS data, which (i) can be used for the detection of multiple variant types, including SNVs, small indels, transposable elements, short tandem repeats, and other large structural variants; (ii) covers all standard steps of NGS analysis, from quality control of raw data and genome alignment to variant calling, annotation, and generation of reports for the interpretation and prioritization of results; (iii) is highly adaptable as it can be deployed and run via either a graphic user interface for non-bioinformaticians and a command line tool for personal computer usage; (iv) is scalable as it can be executed in parallel as a Snakemake workflow, and; (v) is computationally efficient by minimizing RAM and CPU time requirements. AVAILABILITY AND IMPLEMENTATION: DNAscan2 is implemented in Python3 and is available at https://github.com/KHP-Informatics/DNAscanv2.

Large-scale analyses of CAV1 and CAV2 suggest their expression is higher in post-mortem ALS brain tissue and affects survival.

Introduction: Caveolin-1 and Caveolin-2 (CAV1 and CAV2) are proteins associated with intercellular neurotrophic signalling. There is converging evidence that CAV1 and CAV2 (CAV1/2) genes have a role in amyotrophic lateral sclerosis (ALS). Disease-associated variants have been identified within CAV1/2 enhancers, which reduce gene expression and lead to disruption of membrane lipid rafts. Methods: Using large ALS whole-genome sequencing and post-mortem RNA sequencing datasets (5,987 and 365 tissue samples, respectively), and iPSC-derived motor neurons from 55 individuals, we investigated the role of CAV1/2 expression and enhancer variants in the ALS phenotype. Results: We report a differential expression analysis between ALS cases and controls for CAV1 and CAV2 genes across various post-mortem brain tissues and three independent datasets. CAV1 and CAV2 expression was consistently higher in ALS patients compared to controls, with significant results across the primary motor cortex, lateral motor cortex, and cerebellum. We also identify increased survival among carriers of CAV1/2 enhancer mutations compared to non-carriers within Project MinE and slower progression as measured by the ALSFRS. Carriers showed a median increase in survival of 345 days. Discussion: These results add to an increasing body of evidence linking CAV1 and CAV2 genes to ALS. We propose that carriers of CAV1/2 enhancer mutations may be conceptualised as an ALS subtype who present a less severe ALS phenotype with a longer survival duration and slower progression. Upregulation of CAV1/2 genes in ALS cases may indicate a causal pathway or a compensatory mechanism. Given prior research supporting the beneficial role of CAV1/2 expression in ALS patients, we consider a compensatory mechanism to better fit the available evidence, although further investigation into the biological pathways associated with CAV1/2 is needed to support this conclusion.

Harnessing Transcriptomic Signals for Amyotrophic Lateral Sclerosis to Identify Novel Drugs and Enhance Risk Prediction.

Introduction: Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease. This study integrates the latest ALS genome-wide association study (GWAS) summary statistics with functional genomic annotations with the aim of providing mechanistic insights into ALS risk loci, inferring drug repurposing opportunities, and enhancing prediction of ALS risk and clinical characteristics. Methods: Genes associated with ALS were identified using GWAS summary statistic methodology including SuSiE SNP-based fine-mapping, and transcriptome- and proteome-wide association study (TWAS/PWAS) analyses. Using several approaches, gene associations were integrated with the DrugTargetor drug-gene interaction database to identify drugs that could be repurposed for the treatment of ALS. Furthermore, ALS gene associations from TWAS were combined with observed blood expression in two external ALS case-control datasets to calculate polytranscriptomic scores and evaluate their utility for prediction of ALS risk and clinical characteristics, including site of onset, age at onset, and survival. Results: SNP-based fine-mapping, TWAS and PWAS identified 117 genes associated with ALS, with TWAS and PWAS providing novel mechanistic insights. Drug repurposing analyses identified five drugs significantly enriched for interactions with ALS associated genes, with directional analyses highlighting α-glucosidase inhibitors may exacerbate ALS pathology. Additionally, drug class enrichment analysis showed calcium channel blockers may reduce ALS risk. Across the two observed expression target samples, ALS polytranscriptomic scores significantly predicted ALS risk (R2 = 4%; p-value = 2.1×10-21). Conclusions: Functionally-informed analyses of ALS GWAS summary statistics identified novel mechanistic insights into ALS aetiology, highlighted several therapeutic research avenues, and enabled statistically significant prediction of ALS risk.

Calculating variant penetrance from family history of disease and average family size in population-scale data.

BACKGROUND: Genetic penetrance is the probability of a phenotype when harbouring a particular pathogenic variant. Accurate penetrance estimates are important across biomedical fields including genetic counselling, disease research, and gene therapy. However, existing approaches for penetrance estimation require, for instance, large family pedigrees or availability of large databases of people affected and not affected by a disease. METHODS: We present a method for penetrance estimation in autosomal dominant phenotypes. It examines the distribution of a variant among people affected (cases) and unaffected (controls) by a phenotype within population-scale data and can be operated using cases only by considering family disease history. It is validated through simulation studies and candidate variant-disease case studies. RESULTS: Our method yields penetrance estimates which align with those obtained via existing approaches in the Parkinson's disease LRRK2 gene and pulmonary arterial hypertension BMPR2 gene case studies. In the amyotrophic lateral sclerosis case studies, examining penetrance for variants in the SOD1 and C9orf72 genes, we make novel penetrance estimates which correspond closely to understanding of the disease. CONCLUSIONS: The present approach broadens the spectrum of traits for which reliable penetrance estimates can be obtained. It has substantial utility for facilitating the characterisation of disease risks associated with rare variants with an autosomal dominant inheritance pattern. The yielded estimates avoid any kinship-specific effects and can circumvent ascertainment biases common when sampling rare variants among control populations.

RetroSnake: A modular pipeline to detect human endogenous retroviruses in genome sequencing data.

Human endogenous retroviruses (HERVs) integrated into the human genome as a result of ancient exogenous infections and currently comprise ∼8% of our genome. The members of the most recently acquired HERV family, HERV-Ks, still retain the potential to produce viral molecules and have been linked to a wide range of diseases including cancer and neurodegeneration. Although a range of tools for HERV detection in NGS data exist, most of them lack wet lab validation and they do not cover all steps of the analysis. Here, we describe RetroSnake, an end-to-end, modular, computationally efficient, and customizable pipeline for the discovery of HERVs in short-read NGS data. RetroSnake is based on an extensively wet-lab validated protocol, it covers all steps of the analysis from raw data to the generation of annotated results presented as an interactive html file, and it is easy to use by life scientists without substantial computational training. Availability and implementation: The Pipeline and an extensive documentation are available on GitHub.

The SOD1-mediated ALS phenotype shows a decoupling between age of symptom onset and disease duration.

Superoxide dismutase (SOD1) gene variants may cause amyotrophic lateral sclerosis, some of which are associated with a distinct phenotype. Most studies assess limited variants or sample sizes. In this international, retrospective observational study, we compare phenotypic and demographic characteristics between people with SOD1-ALS and people with ALS and no recorded SOD1 variant. We investigate which variants are associated with age at symptom onset and time from onset to death or censoring using Cox proportional-hazards regression. The SOD1-ALS dataset reports age of onset for 1122 and disease duration for 883 people; the comparator population includes 10,214 and 9010 people respectively. Eight variants are associated with younger age of onset and distinct survival trajectories; a further eight associated with younger onset only and one with distinct survival only. Here we show that onset and survival are decoupled in SOD1-ALS. Future research should characterise rarer variants and molecular mechanisms causing the observed variability.

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.

A polymorphic transcriptional regulatory domain in the amyotrophic lateral sclerosis risk gene CFAP410 correlates with differential isoform expression.

We describe the characterisation of a variable number tandem repeat (VNTR) domain within intron 1 of the amyotrophic lateral sclerosis (ALS) risk gene CFAP410 (Cilia and flagella associated protein 410) (previously known as C21orf2), providing insight into how this domain could support differential gene expression and thus be a modulator of ALS progression or risk. We demonstrated the VNTR was functional in a reporter gene assay in the HEK293 cell line, exhibiting both the properties of an activator domain and a transcriptional start site, and that the differential expression was directed by distinct repeat number in the VNTR. These properties embedded in the VNTR demonstrated the potential for this VNTR to modulate CFAP410 expression. We extrapolated these findings in silico by utilisation of tagging SNPs for the two most common VNTR alleles to establish a correlation with endogenous gene expression. Consistent with in vitro data, CFAP410 isoform expression was found to be variable in the brain. Furthermore, although the number of matched controls was low, there was evidence for one specific isoform being correlated with lower expression in those with ALS. To address if the genotype of the VNTR was associated with ALS risk, we characterised the variation of the CFAP410 VNTR in ALS cases and matched controls by PCR analysis of the VNTR length, defining eight alleles of the VNTR. No significant difference was observed between cases and controls, we noted, however, the cohort was unlikely to contain sufficient power to enable any firm conclusion to be drawn from this analysis. This data demonstrated that the VNTR domain has the potential to modulate CFAP410 expression as a regulatory element that could play a role in its tissue-specific and stimulus-inducible regulation that could impact the mechanism by which CFAP410 is involved in ALS.

Characterisation of retrotransposon insertion polymorphisms in whole genome sequencing data from individuals with amyotrophic lateral sclerosis.

The genetics of an individual is a crucial factor in understanding the risk of developing the neurodegenerative disease amyotrophic lateral sclerosis (ALS). There is still a large proportion of the heritability of ALS, particularly in sporadic cases, to be understood. Among others, active transposable elements drive inter-individual variability, and in humans long interspersed element 1 (LINE1, L1), Alu and SINE-VNTR-Alu (SVA) retrotransposons are a source of polymorphic insertions in the population. We undertook a pilot study to characterise the landscape of non-reference retrotransposon insertion polymorphisms (non-ref RIPs) in 15 control and 15 ALS individuals' whole genomes from Project MinE, an international project to identify potential genetic causes of ALS. The combination of two bioinformatics tools (mobile element locator tool (MELT) and TEBreak) identified on average 1250 Alu, 232 L1 and 77 SVA non-ref RIPs per genome across the 30 analysed. Further PCR validation of individual polymorphic retrotransposon insertions showed a similar level of accuracy for MELT and TEBreak. Our preliminary study did not identify a specific RIP or a significant difference in the total number of non-ref RIPs in ALS compared to control genomes. The use of multiple bioinformatic tools improved the accuracy of non-ref RIP detection and our study highlights the potential importance of studying these elements further in ALS.

Amyotrophic lateral sclerosis and cerebellum.

Amyotrophic lateral sclerosis (ALS) is a devastating, heterogeneous neurodegenerative neuromuscular disease that leads to a fatal outcome within 2-5 years, and yet, a precise nature of the association between its major phenotypes and the cerebellar role in ALS pathology remains unknown. Recently, repeat expansions in several genes in which variants appreciably contribute to cerebellar pathology, including C9orf72, NIPA1, ATXN2 and ATXN1, have been found to confer a significant risk for ALS. To better define this relationship, we performed MAGMA gene-based analysis and tissue enrichment analysis using genome-wide association study summary statistics based on a study of 27,205 people with ALS and 110,881 controls. Our preliminary results imply a striking cerebellar tissue specificity and further support increasing calls for re-evaluation of the cerebellar role in the ALS pathology.

GEOexplorer: a webserver for gene expression analysis and visualisation.

Gene Expression Omnibus (GEO) is a database repository hosting a substantial proportion of publicly available high throughput gene expression data. Gene expression analysis is a powerful tool to gain insight into the mechanisms and processes underlying the biological and phenotypic differences between sample groups. Despite the wide availability of gene expression datasets, their access, analysis, and integration are not trivial and require specific expertise and programming proficiency. We developed the GEOexplorer webserver to allow scientists to access, integrate and analyse gene expression datasets without requiring programming proficiency. Via its user-friendly graphic interface, users can easily apply GEOexplorer to perform interactive and reproducible gene expression analysis of microarray and RNA-seq datasets, while producing a wealth of interactive visualisations to facilitate data exploration and interpretation, and generating a range of publication ready figures. The webserver allows users to search and retrieve datasets from GEO as well as to upload user-generated data and combine and harmonise two datasets to perform joint analyses. GEOexplorer, available at https://geoexplorer.rosalind.kcl.ac.uk, provides a solution for performing interactive and reproducible analyses of microarray and RNA-seq gene expression data, empowering life scientists to perform exploratory data analysis and differential gene expression analysis on-the-fly without informatics proficiency.