Multiple-Tissue Integrative Transcriptome-Wide Association Studies Discovered New Genes Associated With Amyotrophic Lateral Sclerosis

TBK1, a prioritized drug repurposing target for amyotrophic lateral sclerosis: evidence from druggable genome Mendelian randomization and pharmacological verification in vitro

BACKGROUND

There is a lack of effective therapeutic strategies for amyotrophic lateral sclerosis (ALS); therefore, drug repurposing might provide a rapid approach to meet the urgent need for treatment.

METHODS

To identify therapeutic targets associated with ALS, we conducted Mendelian randomization (MR) analysis and colocalization analysis using cis-eQTL of druggable gene and ALS GWAS data collections to determine annotated druggable gene targets that exhibited significant associations with ALS. By subsequent repurposing drug discovery coupled with inclusion criteria selection, we identified several drug candidates corresponding to their druggable gene targets that have been genetically validated. The pharmacological assays were then conducted to further assess the efficacy of genetics-supported repurposed drugs for potential ALS therapy in various cellular models.

RESULTS

Through MR analysis, we identified potential ALS druggable genes in the blood, including TBK1 [OR 1.30, 95%CI (1.19, 1.42)], TNFSF12 [OR 1.36, 95%CI (1.19, 1.56)], GPX3 [OR 1.28, 95%CI (1.15, 1.43)], TNFSF13 [OR 0.45, 95%CI (0.32, 0.64)], and CD68 [OR 0.38, 95%CI (0.24, 0.58)]. Additionally, we identified potential ALS druggable genes in the brain, including RESP18 [OR 1.11, 95%CI (1.07, 1.16)], GPX3 [OR 0.57, 95%CI (0.48, 0.68)], GDF9 [OR 0.77, 95%CI (0.67, 0.88)], and PTPRN [OR 0.17, 95%CI (0.08, 0.34)]. Among them, TBK1, TNFSF12, RESP18, and GPX3 were confirmed in further colocalization analysis. We identified five drugs with repurposing opportunities targeting TBK1, TNFSF12, and GPX3, namely fostamatinib (R788), amlexanox (AMX), BIIB-023, RG-7212, and glutathione as potential repurposing drugs. R788 and AMX were prioritized due to their genetic supports, safety profiles, and cost-effectiveness evaluation. Further pharmacological analysis revealed that R788 and AMX mitigated neuroinflammation in ALS cell models characterized by overly active cGAS/STING signaling that was induced by MSA-2 or ALS-related toxic proteins (TDP-43 and SOD1), through the inhibition of TBK1 phosphorylation.

CONCLUSIONS

Our MR analyses provided genetic evidence supporting TBK1, TNFSF12, RESP18, and GPX3 as druggable genes for ALS treatment. Among the drug candidates targeting the above genes with repurposing opportunities, FDA-approved drug-R788 and AMX served as effective TBK1 inhibitors. The subsequent pharmacological studies validated the potential of R788 and AMX for treating specific ALS subtypes through the inhibition of TBK1 phosphorylation.

Advances in sequencing technologies for amyotrophic lateral sclerosis research

Amyotrophic lateral sclerosis (ALS) is caused by upper and lower motor neuron loss and has a fairly rapid disease progression, leading to fatality in an average of 2-5 years after symptom onset. Numerous genes have been implicated in this disease; however, many cases remain unexplained. Several technologies are being used to identify regions of interest and investigate candidate genes. Initial approaches to detect ALS genes include, among others, linkage analysis, Sanger sequencing, and genome-wide association studies. More recently, next-generation sequencing methods, such as whole-exome and whole-genome sequencing, have been introduced. While those methods have been particularly useful in discovering new ALS-linked genes, methodological advances are becoming increasingly important, especially given the complex genetics of ALS. Novel sequencing technologies, like long-read sequencing, are beginning to be used to uncover the contribution of repeat expansions and other types of structural variation, which may help explain missing heritability in ALS. In this review, we discuss how popular and/or upcoming methods are being used to discover ALS genes, highlighting emerging long-read sequencing platforms and their role in aiding our understanding of this challenging disease.

Identifying Candidate Genes Associated with Sporadic Amyotrophic Lateral Sclerosis via Integrative Analysis of Transcriptome-Wide Association Study and Messenger RNA Expression Profile

Amyotrophic lateral sclerosis, a fatal neurodegeneration disease affecting motor neurons in the brain and spinal cord, is difficult to diagnose and treat. The objective of this study is to identify novel candidate genes related to ALS. Transcriptome-wide association study of ALS was conducted by integrating the genome-wide association study summary data (including 1234 ALS patients and 2850 controls) and pre-computed gene expression weights of different tissues. The ALS-associated genes identified by TWAS were further compared with the differentially expressed genes detected by the mRNA expression profiles of the sporadic ALS. Functional enrichment and annotation analysis of identified genes were performed by an R package and the functional mapping and annotation software. TWAS identified 761 significant genes (P_TWAS < 0.05), 627 Gene ontology terms, and 8 Kyoto Encyclopedia of Genes and Genomes pathways for ALS, such as C9orf72, with three expression quantitative trait loci were found significantly: rs2453554 (P_{TWAS CBRS} = 4.68 × 10^-10, P_{TWAS CBRS} = 2.54 × 10^-9), rs10967976 (P_{TWAS CBRS} = 7.85 × 10^-10, P_{TWAS CBRS} = 8.91 × 10^-9, P_{TWAS CBRS} = 1.49 × 10^-7, P_{TWAS CBRS} = 5.59 × 10^-7), rs3849946 (P_{TWAS CBRS} = 7.69 × 10^-4, P_{TWAS YBL} = 4.02 × 10^-2), Mitochondrion (P_adj = 4.22 × 10^-16), and Cell cycle (P_adj = 2.04 × 10^-3). Moreover, 107 common genes, 4 KEGG pathways and 41 GO terms were detected by integrating mRNA expression profiles of sALS, such as CPVL (FC = 2.06, P_mRNA = 6.99 × 10^-6, P_{TWAS CBR} = 2.88 × 10^-2, P_{TWAS CBR} = 4.37 × 10^-2), Pyrimidine Metabolism (P_adj = 2.43 × 10^-2), and Cell Activation (P_adj = 5.54 × 10^-3). Multiple candidate genes and pathways were detected for ALS. Our findings may provide novel clues for understanding the genetic mechanism of ALS.

A comprehensive comparison of multilocus association methods with summary statistics in genome-wide association studies

BACKGROUND

Multilocus analysis on a set of single nucleotide polymorphisms (SNPs) pre-assigned within a gene constitutes a valuable complement to single-marker analysis by aggregating data on complex traits in a biologically meaningful way. However, despite the existence of a wide variety of SNP-set methods, few comprehensive comparison studies have been previously performed to evaluate the effectiveness of these methods.

RESULTS

We herein sought to fill this knowledge gap by conducting a comprehensive empirical comparison for 22 commonly-used summary-statistics based SNP-set methods. We showed that only seven methods could effectively control the type I error, and that these well-calibrated approaches had varying power performance under the simulation scenarios. Overall, we confirmed that the burden test was generally underpowered and score-based variance component tests (e.g., sequence kernel association test) were much powerful under the polygenic genetic architecture in both common and rare variant association analyses. We further revealed that two linkage-disequilibrium-free P value combination methods (e.g., harmonic mean P value method and aggregated Cauchy association test) behaved very well under the sparse genetic architecture in simulations and real-data applications to common and rare variant association analyses as well as in expression quantitative trait loci weighted integrative analysis. We also assessed the scalability of these approaches by recording computational time and found that all these methods can be scalable to biobank-scale data although some might be relatively slow.

CONCLUSION

In conclusion, we hope that our findings can offer an important guidance on how to choose appropriate multilocus association analysis methods in post-GWAS era. All the SNP-set methods are implemented in the R package called MCA, which is freely available at https://github.com/biostatpzeng/ .

Causal Inference of Genetic Variants and Genes in Amyotrophic Lateral Sclerosis

Amyotrophic lateral sclerosis (ALS) is a fatal progressive multisystem disorder with limited therapeutic options. Although genome-wide association studies (GWASs) have revealed multiple ALS susceptibility loci, the exact identities of causal variants, genes, cell types, tissues, and their functional roles in the development of ALS remain largely unknown. Here, we reported a comprehensive post-GWAS analysis of the recent large ALS GWAS (n = 80,610), including functional mapping and annotation (FUMA), transcriptome-wide association study (TWAS), colocalization (COLOC), and summary data-based Mendelian randomization analyses (SMR) in extensive multi-omics datasets. Gene property analysis highlighted inhibitory neuron 6, oligodendrocytes, and GABAergic neurons (Gad1/Gad2) as functional cell types of ALS and confirmed cerebellum and cerebellar hemisphere as functional tissues of ALS. Functional annotation detected the presence of multiple deleterious variants at three loci (9p21.2, 12q13.3, and 12q14.2) and highlighted a list of SNPs that are potentially functional. TWAS, COLOC, and SMR identified 43 genes at 24 loci, including 23 novel genes and 10 novel loci, showing significant evidence of causality. Integrating multiple lines of evidence, we further proposed that rs2453555 at 9p21.2 and rs229243 at 14q12 functionally contribute to the development of ALS by regulating the expression of C9orf72 in pituitary and SCFD1 in skeletal muscle, respectively. Together, these results advance our understanding of the biological etiology of ALS, feed into new therapies, and provide a guide for subsequent functional experiments.

Integrative eQTL-weighted hierarchical Cox models for SNP-set based time-to-event association studies

BACKGROUND

Integrating functional annotations into SNP-set association studies has been proven a powerful analysis strategy. Statistical methods for such integration have been developed for continuous and binary phenotypes; however, the SNP-set integrative approaches for time-to-event or survival outcomes are lacking.

METHODS

We here propose IEHC, an integrative eQTL (expression quantitative trait loci) hierarchical Cox regression, for SNP-set based survival association analysis by modeling effect sizes of genetic variants as a function of eQTL via a hierarchical manner. Three p-values combination tests are developed to examine the joint effects of eQTL and genetic variants after a novel decorrelated modification of statistics for the two components. An omnibus test (IEHC-ACAT) is further adapted to aggregate the strengths of all available tests.

RESULTS

Simulations demonstrated that the IEHC joint tests were more powerful if both eQTL and genetic variants contributed to association signal, while IEHC-ACAT was robust and often outperformed other approaches across various simulation scenarios. When applying IEHC to ten TCGA cancers by incorporating eQTL from relevant tissues of GTEx, we revealed that substantial correlations existed between the two types of effect sizes of genetic variants from TCGA and GTEx, and identified 21 (9 unique) cancer-associated genes which would otherwise be missed by approaches not incorporating eQTL.

CONCLUSION

IEHC represents a flexible, robust, and powerful approach to integrate functional omics information to enhance the power of identifying association signals for the survival risk of complex human cancers.

Leveraging Methylation Alterations to Discover Potential Causal Genes Associated With the Survival Risk of Cervical Cancer in TCGA Through a Two-Stage Inference Approach

BACKGROUND

Multiple genes were previously identified to be associated with cervical cancer; however, the genetic architecture of cervical cancer remains unknown and many potential causal genes are yet to be discovered.

METHODS

To explore potential causal genes related to cervical cancer, a two-stage causal inference approach was proposed within the framework of Mendelian randomization, where the gene expression was treated as exposure, with methylations located within the promoter regions of genes serving as instrumental variables. Five prediction models were first utilized to characterize the relationship between the expression and methylations for each gene; then, the methylation-regulated gene expression (MReX) was obtained and the association was evaluated via Cox mixed-effect model based on MReX. We further implemented the aggregated Cauchy association test (ACAT) combination to take advantage of respective strengths of these prediction models while accounting for dependency among the p-values.

RESULTS

A total of 14 potential causal genes were discovered to be associated with the survival risk of cervical cancer in TCGA when the five prediction models were separately employed. The total number of potential causal genes was brought to 23 when conducting ACAT. Some of the newly discovered genes may be novel (e.g., YJEFN3, SPATA5L1, IMMP1L, C5orf55, PPIP5K2, ZNF330, CRYZL1, PPM1A, ESCO2, ZNF605, ZNF225, ZNF266, FICD, and OSTC). Functional analyses showed that these genes were enriched in tumor-associated pathways. Additionally, four genes (i.e., COL6A1, SYDE1, ESCO2, and GIPC1) were differentially expressed between tumor and normal tissues.

CONCLUSION

Our study discovered promising candidate genes that were causally associated with the survival risk of cervical cancer and thus provided new insights into the genetic etiology of cervical cancer.

An Integrative Transcriptome-Wide Analysis of Amyotrophic Lateral Sclerosis for the Identification of Potential Genetic Markers and Drug Candidates

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative neuromuscular disease. Although genome-wide association studies (GWAS) have successfully identified many variants significantly associated with ALS, it is still difficult to characterize the underlying biological mechanisms inducing ALS. In this study, we performed a transcriptome-wide association study (TWAS) to identify disease-specific genes in ALS. Using the largest ALS GWAS summary statistic (n = 80,610), we identified seven novel genes using 19 tissue reference panels. We conducted a conditional analysis to verify the genes’ independence and to confirm that they are driven by genetically regulated expressions. Furthermore, we performed a TWAS-based enrichment analysis to highlight the association of important biological pathways, one in each of the four tissue reference panels. Finally, utilizing a connectivity map, a database of human cell expression profiles cultured with bioactive small molecules, we discovered functional associations between genes and drugs to identify 15 bioactive small molecules as potential drug candidates for ALS. We believe that, by integrating the largest ALS GWAS summary statistic with gene expression to identify new risk loci and causal genes, our study provides strong candidates for molecular basis experiments in ALS.

Detection of Genetic Overlap Between Rheumatoid Arthritis and Systemic Lupus Erythematosus Using GWAS Summary Statistics

Background

Clinical and epidemiological studies have suggested systemic lupus erythematosus (SLE) and rheumatoid arthritis (RA) are comorbidities and common genetic etiologies can partly explain such coexistence. However, shared genetic determinations underlying the two diseases remain largely unknown.

Methods

Our analysis relied on summary statistics available from genome-wide association studies of SLE (N = 23,210) and RA (N = 58,284). We first evaluated the genetic correlation between RA and SLE through the linkage disequilibrium score regression (LDSC). Then, we performed a multiple-tissue eQTL (expression quantitative trait loci) weighted integrative analysis for each of the two diseases and aggregated association evidence across these tissues via the recently proposed harmonic mean P-value (HMP) combination strategy, which can produce a single well-calibrated P-value for correlated test statistics. Afterwards, we conducted the pleiotropy-informed association using conjunction conditional FDR (ccFDR) to identify potential pleiotropic genes associated with both RA and SLE.

Results

We found there existed a significant positive genetic correlation (r_g = 0.404, P = 6.01E-10) via LDSC between RA and SLE. Based on the multiple-tissue eQTL weighted integrative analysis and the HMP combination across various tissues, we discovered 14 potential pleiotropic genes by ccFDR, among which four were likely newly novel genes (i.e., INPP5B, OR5K2, RP11-2C24.5, and CTD-3105H18.4). The SNP effect sizes of these pleiotropic genes were typically positively dependent, with an average correlation of 0.579. Functionally, these genes were implicated in multiple auto-immune relevant pathways such as inositol phosphate metabolic process, membrane and glucagon signaling pathway.

Conclusion

This study reveals common genetic components between RA and SLE and provides candidate associated loci for understanding of molecular mechanism underlying the comorbidity of the two diseases.

Aggregating multiple expression prediction models improves the power of transcriptome-wide association studies

Transcriptome-wide association study (TWAS) is an important integrative method for identifying genes that are causally associated with phenotypes. A key step of TWAS involves the construction of expression prediction models for every gene in turn using its cis-SNPs as predictors. Different TWAS methods rely on different models for gene expression prediction, and each such model makes a distinct modeling assumption that is often suitable for a particular genetic architecture underlying expression. However, the genetic architectures underlying gene expression vary across genes throughout the transcriptome. Consequently, different TWAS methods may be beneficial in detecting genes with distinct genetic architectures. Here, we develop a new method, HMAT, which aggregates TWAS association evidence obtained across multiple gene expression prediction models by leveraging the harmonic mean P-value combination strategy. Because each expression prediction model is suited to capture a particular genetic architecture, aggregating TWAS associations across prediction models as in HMAT improves accurate expression prediction and enables subsequent powerful TWAS analysis across the transcriptome. A key feature of HMAT is its ability to accommodate the correlations among different TWAS test statistics and produce calibrated P-values after aggregation. Through numerical simulations, we illustrated the advantage of HMAT over commonly used TWAS methods as well as ad hoc P-value combination rules such as Fisher's method. We also applied HMAT to analyze summary statistics of nine common diseases. In the real data applications, HMAT was on average 30.6% more powerful compared to the next best method, detecting many new disease-associated genes that were otherwise not identified by existing TWAS approaches. In conclusion, HMAT represents a flexible and powerful TWAS method that enjoys robust performance across a range of genetic architectures underlying gene expression.