Neuroinflammation impacts on the progression of amyotrophic lateral sclerosis (ALS), a fatal neurodegenerative disorder. Specialized pro-resolving mediators trigger the resolution of inflammation. We investigate the specialized pro-resolving mediator blood profile and their receptors' expression in peripheral blood mononuclear cells in relation to survival in ALS. People living with ALS (pwALS) were stratified based on bulbar versus limb onset and on key progression metrics using a latent class model, to separate faster progressing from slower progressing ALS. Specialized pro-resolving mediator blood concentrations were measured at baseline and in one additional visit in 20 pwALS and 10 non-neurological controls (Cohort 1). Flow cytometry was used to study the GPR32 and GPR18 resolvin receptors' expression in peripheral blood mononuclear cells from 40 pwALS and 20 non-neurological controls (Cohort 2) at baseline and in two additional visits in 17 pwALS. Survival analysis was performed using Cox proportional hazards models, including known clinical predictors and GPR32 and GPR18 mononuclear cell expression. Differential expression and linear discriminant analyses showed that plasma resolvins were able to distinguish phenotypic variants of ALS from non-neurological controls. RvE3 was elevated in blood from pwALS, whilst RvD1, RvE3, RvT4 and RvD1(n-3 DPA) were upregulated in A-S and RvD2 in A-F. Compared to non-neurological controls, GPR32 was upregulated in monocytes expressing the active inflammation-suppressing CD11b(+) integrin from fast-progressing pwALS, including those with bulbar onset disease (P < 0.0024), whilst GPR32 and GPR18 were downregulated in most B and T cell subtypes. Only GPR18 was upregulated in na & iuml;ve double positive Tregs, memory cytotoxic Tregs, senescent late memory B cells and late senescent CD8(+) T cells from pwALS compared to non-neurological controls (P < 0.0431). Higher GPR32 and GPR18 median expression in blood mononuclear cells was associated with longer survival, with GPR32 expression in classical monocytes (hazard ratio: 0.11, P = 0.003) and unswitched memory B cells (hazard ratio: 0.44, P = 0.008) showing the most significant association, along with known clinical predictors. Low levels of resolvins and downregulation of their membrane receptors in blood mononuclear cells are linked to a faster progression of ALS. Higher mononuclear cell expression of resolvin receptors is a predictor of longer survival. These findings suggest a lipid-mediated neuroprotective response that could be harnessed to develop novel therapeutic strategies and biomarkers for ALS.
Continued methodological advances have enabled numerous statistical approaches for the analysis of summary statistics from genome-wide association studies. Genetic correlation analysis within specific regions enables a new strategy for identifying pleiotropy. Genomic regions with significant ‘local’ genetic correlations can be investigated further using state-of-the-art methodologies for statistical fine-mapping and variant colocalisation. We explored the utility of a genome-wide local genetic correlation analysis approach for identifying genetic overlaps between the candidate neuropsychiatric disorders, Alzheimer’s disease (AD), amyotrophic lateral sclerosis (ALS), frontotemporal dementia, Parkinson’s disease, and schizophrenia. The correlation analysis identified several associations between traits, the majority of which were loci in the human leukocyte antigen region. Colocalisation analysis suggested that disease-implicated variants in these loci often differ between traits and, in one locus, indicated a shared causal variant between ALS and AD. Our study identified candidate loci that might play a role in multiple neuropsychiatric diseases and suggested the role of distinct mechanisms across diseases despite shared loci. The fine-mapping and colocalisation analysis protocol designed for this study has been implemented in a flexible analysis pipeline that produces HTML reports and is available at: https://github.com/ThomasPSpargo/COLOC-reporter.
Background Amyotrophic lateral sclerosis (ALS) is a fatal heterogeneous neurodegenerative disease that typically leads to death from respiratory failure within two to five years. Despite the identification of several genetic risk factors, the biological processes involved in ALS pathogenesis remain poorly understood. The motor cortex is an ideal region to study dysregulated pathological processes in ALS as it is affected from the earliest stages of the disease. In this study, we investigated motor-cortex gene expression of cases and controls to gain new insight into the molecular footprint of ALS. Methods We performed a large case-control differential expression analysis of two independent post-mortem motor cortex bulk RNA-sequencing (RNAseq) datasets from the King’s College London BrainBank (N = 171) and TargetALS (N = 132). Differentially expressed genes from both datasets were subjected to gene and pathway enrichment analysis. Genes common to both datasets were also reviewed for their involvement with known mechanisms of ALS pathogenesis to identify potential candidate genes. Finally, we performed a correlation analysis of genes implicated in pathways enriched in both datasets with clinical outcomes such as the age of onset and survival. Results Differential expression analysis identified 2,290 and 402 differentially expressed genes in KCL BrainBank and TargetALS cases, respectively. Enrichment analysis revealed significant synapse-related processes in the KCL BrainBank dataset, while the TargetALS dataset carried an immune system-related signature. There were 44 differentially expressed genes which were common to both datasets, which represented previously recognised mechanisms of ALS pathogenesis, such as lipid metabolism, mitochondrial energy homeostasis and neurovascular unit dysfunction. Differentially expressed genes in both datasets were significantly enriched for the neuropeptide signalling pathway. By looking at the relationship between the expression of neuropeptides and their receptors with clinical measures, we found that in both datasets NPBWR1, TAC3 and SSTR1 correlated with age of onset, and GNRH1, TACR1 with survival. We provide access to gene-level expression results to the broader research community through a publicly available web application ( https://alsgeexplorer.er.kcl.ac.uk ). Conclusion This study identified motor-cortex specific pathways altered in ALS patients, potential molecular targets for therapeutic disease intervention and a set of neuropeptides and receptors for investigation as potential biomarkers.
Abstract 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 .
ABSTRACT Background Amyotrophic lateral sclerosis (ALS) displays considerable clinical, genetic and molecular heterogeneity. Machine learning approaches have shown potential to disentangle complex disease landscapes and they have been utilised for patient stratification in ALS. However, lack of independent validation in different populations and in pre-mortem tissue samples have greatly limited their use in clinical and research settings. We overcame such issues by performing a large-scale study of over 600 post-mortem brain and blood samples of people with ALS from four independent datasets from the UK, Italy, the Netherlands and the US. Methods Hierarchical clustering was performed on the 5000 most variably expressed autosomal genes identified from post-mortem motor cortex expression data of people with sporadic ALS from the KCL BrainBank (N=112). The molecular architectures of each cluster were investigated with gene enrichment, network and cell composition analysis. Methylation and genetic data were also used to assess if other omics measures differed between individuals. Validation of these clusters was achieved by applying linear discriminant analysis models based on the KCL BrainBank to the TargetALS US motor cortex (N=93), as well as Italian (N=15) and Dutch (N=397) blood expression datasets. Phenotype analysis was also performed to assess cluster-specific differences in clinical outcomes. Results We identified three molecular phenotypes, which reflect the proposed major mechanisms of ALS pathogenesis: synaptic and neuropeptide signalling, excitotoxicity and oxidative stress, and neuroinflammation. Known ALS risk genes were identified among the informative genes of each cluster, suggesting potential for genetic profiling of the molecular phenotypes. Cell types which are known to be associated with specific molecular phenotypes were found in higher proportions in those clusters. These molecular phenotypes were validated in independent motor cortex and blood datasets. Phenotype analysis identified distinct cluster-related outcomes associated with progression, survival and age of death. We developed a public webserver ( https://alsgeclustering.er.kcl.ac.uk ) that allows users to stratify samples with our model by uploading their expression data. Conclusions We have identified three molecular phenotypes, driven by different cell types, which reflect the proposed major mechanisms of ALS pathogenesis. Our results support the hypothesis of biological heterogeneity in ALS where different mechanisms underly ALS pathogenesis in a subgroup of patients that can be identified by a specific expression signature. These molecular phenotypes show potential for stratification of clinical trials, the development of biomarkers and personalised treatment approaches.
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.
Abstract Background: Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease characterised by a highly variable clinical presentation and multifaceted genetic and biological bases that translate into great patient heterogeneity. The identification of homogeneous subgroups of patients in terms of both clinical presentation and biological causes, could favour the development of effective treatments, healthcare, and clinical trials. We aimed to identify and characterise homogenous clinical subgroups of ALS, examining whether they represent underlying biological trends. Methods: Latent class clustering analysis, an unsupervised machine-learning method, was used to identify homogenous subpopulations in 6,523 people with ALS from Project MinE, using widely collected ALS-related clinical variables. The clusters were validated using 7,829 independent patients from STRENGTH. We tested whether the identified subgroups were associated with biological trends in genetic variation across genes previously linked to ALS, polygenic risk scores of ALS and related neuropsychiatric traits, and in gene expression data from post-mortem motor cortex samples. Results: We identified five ALS subgroups based on patterns in clinical data which were general across international datasets. Distinct genetic trends were observed for rare variants in the SOD1 and C9orf72 genes, and across genes implicated in biological processes relevant to ALS. Polygenic risk scores of ALS, schizophrenia and Parkinson's disease were also higher in distinct clusters with respect to controls. Gene expression analysis identified different altered biological processes across clusters reflecting the genetic differences. We developed a machine learning classifier based on our model to assign subgroup membership using clinical data available at first visit, and made it available on a public webserver at http://latentclusterals.er.kcl.ac.uk. Conclusion: ALS subgroups characterised by highly distinct clinical presentations were discovered and validated in two large independent international datasets. Such groups were also characterised by different underlying genetic architectures and biology. Our results showed that data-driven patient stratification into more clinically and biologically homogeneous subtypes of ALS is possible and could help develop more effective and targeted approaches to the biomedical and clinical study of ALS.
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.
Background Over the past three decades there have been numerous molecular biology developments that have led to an explosion in the number of gene expression studies being performed. Many of these gene expression studies publish their data to the public database GEO, making them freely available. By analysing gene expression datasets, researchers can identify genes that are differentially expressed between two groups. This can provide insights that lead to the development of new tests and treatments for diseases. Despite the wide availability of gene expression datasets, analysing them is difficult for several reasons. These reasons include the fact that most methods for performing gene expression analysis require programming proficiency. Results We developed the GEOexplorer software package to overcome several of the difficulties in performing gene expression analysis. GEOexplorer was therefore developed as a web application, that can perform interactive and reproducible microarray gene expression analysis, while producing a wealth of interactive visualisations to facilitate result exploration. GEOexplorer is implemented in R using the Shiny framework and is fully integrated with the existing core structures of the Bioconductor project. Users can perform the essential steps of exploratory data analysis and differential gene expression analysis intuitively and generate a broad spectrum of publication ready outputs. Conclusion GEOexplorer is distributed as an R package in the Bioconductor project ( http://bioconductor.org/packages/GEOexplorer/ ). GEOexplorer provides a solution for performing interactive and reproducible analyses of microarray gene expression data, empowering life scientists to perform exploratory data analysis and differential gene expression analysis on GEO microarray datasets.