Yong Heng Wang1,2†
Pan Pan Luo1†
Ao Yi Geng1Xinwei Li3
Tai-Hang Liu1,2Yi Jie He1Lin Huang1
Ya Qin Tang1*- 1Department of Bioinformatics, School of Basic Medical Sciences, Chongqing Medical University, Chongqing, China
- 2Joint International Research Laboratory of Reproduction and Development, Chongqing Medical University, Chongqing, China
- 3School of Microelectronics and Communication Engineering, Chongqing University, Chongqing, China
Numerous genetic variants associated with Alzheimer’s disease (AD) have been identified through genome-wide association studies (GWAS), but their interpretation is hindered by the strong linkage disequilibrium (LD) among the variants, making it difficult to identify the causal variants directly. To address this issue, the transcriptome-wide association study (TWAS) was employed to infer the association between gene expression and a trait at the genetic level using expression quantitative trait locus (eQTL) cohorts. In this study, we applied the TWAS theory and utilized the improved Joint-Tissue Imputation (JTI) approach and Mendelian Randomization (MR) framework (MR-JTI) to identify potential AD-associated genes. By integrating LD score, GTEx eQTL data, and GWAS summary statistic data from a large cohort using MR-JTI, a total of 415 AD-associated genes were identified. Then, 2873 differentially expressed genes from 11 AD-related datasets were used for the Fisher test of these AD-associated genes. We finally obtained 36 highly reliable AD-associated genes, including APOC1, CR1, ERBB2, and RIN3. Moreover, the GO and KEGG enrichment analysis revealed that these genes are primarily involved in antigen processing and presentation, amyloid-beta formation, tau protein binding, and response to oxidative stress. The identification of these potential AD-associated genes not only provides insights into the pathogenesis of AD but also offers biomarkers for early diagnosis of the disease.
Introduction
Alzheimer’s disease (AD) is the most common form of dementia, accounting for approximately 50–60% of dementia cases. As a progressive neurodegenerative disorder, AD is characterized by a gradual loss of memory and cognition, with amyloid plaques and neurofibrillary tangles being the primary pathological features (Hansson et al., 2018). The recent data suggest that there are around 55 million people worldwide living with dementia, and the prevalence of dementia is expected to triple by 2050 (Khan et al., 2020). This increase will impose a significant burden on both families and society, making it crucial to analyze the pathogenic mechanism and identify potential risk factors for AD diagnosis and treatment.
Genome-wide Association Studies (GWAS) have emerged as a powerful tool for investigating complex diseases, leading to the discovery of over 40 AD-associated risk alleles through large cohort studies on AD (Jansen et al., 2019; Scheltens et al., 2021). However, due to the strong linkage disequilibrium (LD) among variants, the loci identified by GWAS cannot be interpreted directly, which further obscures the causality between variants and phenotypes. Therefore, GWAS data alone is insufficient for determining the causal genes and underlying regulatory mechanisms. To fill this gap, transcriptome-wide association studies (TWAS) are valuable developed for studying potential gene regulatory mechanisms associated with variable traits by integrating transcriptomic and genetic data (Tang et al., 2021). Traditional TWAS analysis typically involves three stages: (i) Construction of association between single nucleotide polymorphisms (SNPs) and gene expression using weighted calculations based on reference panels such as the GTEx database or other cohort data containing both genotyping and expression (known as the expression quantitative trait locus or eQTL model). (ii) Filling the lacked gene expression data in the large-scale GWAS cohort using the trained eQTL model. (iii) Utilizing the filled gene expression data to infer associations between gene expression and traits (Luningham et al., 2020). To achieve more accurate results, several improved methods have been developed based on the traditional TWAS approach, including PrediXcan (Gamazon et al., 2015), S-PrediXcan (Barbeira et al., 2018), UTMOST (Rodriguez-Fontenla and Carracedo, 2021), CoMM (Yang et al., 2019), PMR-Egger (Yuan et al., 2020), moPMR-Egger (Liu et al., 2021), VC-TWAS (Tang et al., 2021), TIGAR (Nagpal et al., 2019). TWAS has proven successful in integrating transcriptomic and genetic data to study various complex human diseases, such as schizophrenia, breast cancer, prostate cancer, and Crohn’s disease (Tang et al., 2021).
In 2020, Zhou et al. (2020) presented a new method to optimize TWAS, namely the Joint-Tissue Imputation (JTI) approach and a Mendelian Randomization (MR) framework for causal inference (MR-JTI). Traditional TWAS methods often fail to fully capitalize on the shared biological characteristics across multiple tissues in the GTEx dataset during the training of the prediction model. Consequently, the prediction accuracy decreases. In contrast, JTI effectively leverages the common regulatory architecture of gene expression across multiple tissues. When the transcriptional regulation of the target gene in simple tissue is specific, it will automatically restore the model to a single tissue prediction model, PrediXcan, thus improving tissue specificity in prediction. These advancements make JTI a superior prediction method compared to PrediXcan, BSLMM, TIGAR, and UTMOST models (Zhou et al., 2013). Furthermore, JTI incorporates Mendelian Randomization (MR) into its framework, significantly enhancing the evidence level and credibility of association analysis. This integration helps address false positive issues caused by horizontal gene pleiotropy and potential confounding factors. In this study, we conducted tissue-specific TWAS for AD across 13 brain regions and blood. We utilized GWAS data and eQTL cohorts from GTEx (version 8) and employed the MR-JTI method. The experimental flow chart detailing the process is depicted in Figure 1. Our findings provide more reliable potential biomarkers and targets for investigating the pathogenesis of AD.
Figure 1. The work flow chart of the experiment.
Materials and methods
Public GWAS statistics data of AD
The GWAS summary statistics data for AD were obtained from the published study in 2022, which included 111,326 AD cases and 677,663 controls (Bellenguez et al., 2022). The data were downloaded from the European Bioinformatics Institute GWAS Catalog1 with accession number GCST90027158. For detailed information about the sample collection, analysis methods, and results, please refer to the original publication (Bellenguez et al., 2022).
GTEx quantitative trait loci (eQTL) data
The GTEx project is the continuously updated public database of human genetic resources. Its latest version V8, contains massive sequencing data of 54 non-diseased tissue sites from 948 donors (Consortium et al., 2017). Using the gene expression and genotype data, the GTEx project performed the eQTL analysis of each tissue site, and the eQTL models can be obtained from its website.2 Zhou et al. (2020) developed the joint-tissue imputation (JTI) approach by considering the shared genetic regulation effects across different tissues and the unique genetic regulation in the target tissues. They trained the GTEx eQTL data with JTI to enhance the prediction performance of instrumental variables for genes. The JTI pre-training models was downloadable, including the eQTL summary statistics and the SNP-SNP covariance matrices3 (Zhou et al., 2020). Here, we obtained the trained pre-trained prediction models for 13 brain regions and blood, including brain amygdala, brain anterior cingulate cortex, brain caudate basal ganglia, brain cerebellar hemisphere, brain cerebellum, brain cortex, brain frontal cortex, brain hippocampus, brain hypothalamus, brain nucleus accumbens basal ganglia, brain putamen basal ganglia, brain spinal cord cervical, and brain substantia nigra.
Joint-tissue imputation (JTI) and Mendelian randomization (MR) integrative analysis (MR-JTI)
The TWAS analysis was performed by JTI using the pre-training eQTL models and GWAS summary statistic data. The P-values of genes in JTI analysis were adjusted using the p.adjust function in R (version 4.1.3) with the false discovery rate (FDR) method. Genes with FDR < 0.05 were identified as AD-associated genes. While JTI established the relationship between gene expression and AD, it remained unclear whether the differential expression of these genes was the cause or the result of AD. Consequently, the identified AD-associated genes were employed for subsequent MR analysis to elucidate the causal relationship between candidate genes and AD. To mitigate confounding effects, LD scores were calculated. These LD scores were obtained using the GCTA software, utilizing data from the 1000 Genomes Project,4 which created a catalog of common human genetic variation by using openly consented samples from people who declared themselves to be healthy. The MR-JTI analysis combined LD scores, GTEx eQTL data, and GWAS summary statistics to obtain potential AD-associated genes.
Analysis of differentially expressed genes in AD
Gene Expression Omnibus (GEO)5 database is a public genomic database containing the entire gene expression data, chips, and microarrays. The 11 published datasets of AD-associated were retrieved and downloaded from the GEO database (Table 1). Differentially expressed genes (DEGs) were identified using a P < 0.05 and | log fold-change (FC)| > 1. The gene symbols of datasets were annotated using DAVID online software.6 Finally, the AD-associated genes identified by MR-JTI and the AD-associated DEGs were tested by the Fisher test to confirm the notable intersection genes.
Table 1. Details of AD related GEO datasets.
Pathway and functional enrichment analysis of AD-associated genes
The Kyoto encyclopedia of genes and genomes (KEGG)7 and Gene Ontology (GO)8 were used for the pathway and functional enrichment analysis of AD-associated genes by R package ClusterProfiler. The GO enrichment analysis includes biological pathways (biological process, BP), cellular components (CC), and molecular function (MF). Moreover, the R package Circlize was used to visualize the high reliable genes and the significantly enriched pathways (adjusted P < 0.05) in which they were involved. Whether the high reliable genes and these pathways have been confirmed to be associated with immunity and AD were evaluated by manual literature search. The high reliable genes were further compared with the immune genes in innateDB9 (Breuer et al., 2013) and ImmPort10 (Bhattacharya et al., 2018) database to confirm their relationship with immunity.
Results
Identification of 415 AD-associated genes by MR-JTI
To identify risk factors associated with AD, we performed the TWAS analysis using JTI pre-training models and GWAS summary statistic data from AD in 2022 (Zhou et al., 2020; Bellenguez et al., 2022). JTI integrates shared genetic regulation effects across multiple tissues and tissue-specific genetic regulation effects, providing prediction models for each tissue. While AD symptoms are primarily linked to hippocampal and frontal cortex lesions, brain amyloid plaques and atrophy can occur throughout the brain in AD patients (Khan et al., 2020). Furthermore, there has been increased attention on the use of blood biomarkers for AD diagnosis and early screening (Blennow and Zetterberg, 2018). Consequently, we downloaded the corresponding JTI pre-training model and performed TWAS analysis to identify AD risk genes in 13 brain regions and blood, respectively. After removing the duplicate genes, a total of 670 potential genes (FDR < 0.05, Supplementary Table 1) were obtained from 13 brain regions and blood. Since JTI has finally established a relationship between gene expression and AD, it is speculated that the expression pattern of these risk genes is associated with the genetic risk of AD.
However, we are uncertain whether the expression changes of these AD-associated risk factors are the cause or outcome of the disease. Therefore, the MR-JTI analysis was conducted on the genes identified by JTI, integrating LD score, eQTL data, and GWAS summary statistic data. After removing the duplicate genes, 415 AD potential causal risk genes (FDR < 0.05, Supplementary Table 2) were identified from 13 brain regions and blood (Supplementary Table 1). The Manhattan plot illustrates the AD-associated risk factors in 13 brain regions and blood screened by MR-JTI (Figure 2 and Supplementary Figure 1). Among the top five genes in the brain amygdala region, RAB8B and HLA-DOB have been confirmed to be involved in the pathogenesis of AD (Figure 2A; Patel et al., 2021; Martinez et al., 2023). Strong evidences indicates that abnormalities of CR1, APOC1, APOC2, LACTB, and ABCA7 were closely related to AD (Figures 2B, C; Karch and Goate, 2015; Shao et al., 2018; De Roeck et al., 2019; Kulminski et al., 2022; Yu et al., 2022). It is worth noting that the top five genes in the blood have been confirmed to be associated with AD, including HAL-DQB2, PLCALM, APH1B, CNN2, and CEACAM19 (Figure 2D; Logue et al., 2018; Seripa et al., 2018; Gockley et al., 2021; Park et al., 2021; Zhang et al., 2022b).
Figure 2. Manhattan plots illustrating MR-JTI results in different brain regions and blood. The vertical axis is the corresponding “-log (P-value)” of each gene in the JTI result; the higher the corresponding value of the gene point, the higher the association between the gene and AD. Brain amygdala (A); brain cortex (B), brain nucleus accumbens basal ganglia (C), and blood (D).
Enrichment analysis of AD-associated genes
To further explore the relationship between the 415 genes and AD, the pathway and functional enrichment analysis was performed by GO and KEGG. The GO results revealed that the 415 genes were enriched in immune-related pathways (MHC class II protein complex assembly, leukocyte mediated immunity, and T cell activation), ERK1 and ERK2 cascade, Tau protein binding, and Ubiquitin binding, especially enriched in amyloid-beta formation (including genes CLU, APH1B, ABCA7, BIN1, and PICALM) (Figure 3A and Supplementary Table 3). The results of KEGG functional enrichment analysis showed that the AD-associated genes enriched in antigen processing and presentation, lysosome, and Th1 and Th2 cell differentiation (Figure 3B and Supplementary Table 4).
Figure 3. Enrichment analysis of AD-associated genes. (A) The GO enrichment analysis of 415 AD-associated genes, including molecular biological process (BP), cellular components (CC), and function (MF). (B) The KEGG enrichment analysis of 415 AD-associated genes.
Fisher test for AD-associated genes and DEGs
To enhance the reliability of screening AD-associated risk factors, we conducted a Fisher test on 415 AD-associated genes using differentially expressed genes (DEGs) in AD. Firstly, we searched the GEO database and analyzed the DEGs in 11 AD-related datasets, comprising eight different brain region datasets and three peripheral blood datasets (Table 1). A total of 2737 unique DEGs (P < 0.05 and | log2FC| > 1) were identified from eight AD-related GEO datasets of brain regions (Supplementary Table 5), and 147 unique DEGs (P < 0.05 and | log2FC| > 1) were screened from three peripheral blood GEO datasets (Supplementary Table 5). Venn analysis of 378 AD-associated genes in the brain regions and 2737 DEGs in the brain regions showed that there were 34 common highly reliable genes associated with AD (Figure 4A and Supplementary Table 6). These genes contain APOC1, CR1, CISD1, and others. Additionally, there were two common highly reliable genes associated with AD that were shared between the 72 AD-associated genes in the blood and the 147 DEGs in the blood (Figure 4B), including LAT2 and NDUFS2. Their roles in AD pathology are listed in Table 2.
Figure 4. Overlapping differentially expressed genes (DEGs) with causal genes in brain regions (A) and blood (B).
Table 2. The details of highly reliable common genes in AD.
To provide insights for further research on the function of these highly reliable genes in AD, the relationship between the pathways involved in the highly reliable genes and AD was evaluated through manual literature search (Figure 5). The results showed that most of these pathways have either been previously reported or are potentially associated with the pathology of AD, including MHC protein complex assembly, and antigen processing and presentation. Furthermore, it was found that the majority of these pathways were related to immunity. To confirm whether these pathways are indeed immune-related, an additional literature search was conducted.
Figure 5. Correspondence between highly reliable genes and enriched pathways. From the outer circle to inner circle, the first circle represents an index containing 20 genes and 95 pathways (enriched by the 415 potential AD-associated genes with adjusted P < 0.05); the second circle denotes the gene or pathway type; the third circle is whether the gene/pathway has been reported to be related to AD; the fourth circle indicates whether the gene/pathway is immune-related. Gray lines indicate the correspondence between genes and pathways.
Discussion
With the increase in human life expectancy and the intensification of the aging society, the prevalence of AD is also on the rise (Alzheimer’s Association, 2021). However, due to the complex etiology and unknown pathogenesis of AD, effective intervention measures in the clinic are lacking. Additionally, there is a lack of effective targets, creating obstacles for drug discovery. Therefore, we performed TWAS using JTI to identify more potential AD risk factors. A total of 670 AD potential risk genes were identified from 13 brain regions and blood, and further causal inference with MR-JTI identified 415 AD-associated risk genes. Currently, the mainstream view is that one of the causes of AD is the deposition of β-amyloid protein (Aβ) and Tau protein resulting in the death of massive neurons (Masters and Selkoe, 2012). These 415 risk genes also enriched in related pathways, including amyloid-beta formation and Tau protein binding. These pathways involve CLU, APH1B, ABCA7, BIN1, PICALM, and MARK4, all of which have been widely proven to be closely related to the pathogenesis of AD by GWAS and numbers experiments (Jansen et al., 2019; Nott et al., 2019; Li et al., 2021; Schwartzentruber et al., 2021). For example, ABCA7 may stimulate cholesterol efflux from cells into lipoprotein particles and further inhibiting/clearing Aβ aggregates and influencing the risk of AD (Chan et al., 2008; Wildsmith et al., 2013). Phosphorylated MARK4 was highly expressed in AD brain tissue, and it could phosphorylate tau at Ser262/356, contributing to tau accumulation, toxicity, and subsequent neurodegeneration (Oba et al., 2020; Waseem et al., 2021). Notably, many genes enriched in immune-related pathways, mainly involving the members of the histocompatibility Complex, Class II, such as HLA-DQA2, HLA-DOB and HLA-DRB1 (Supplementary Tables 3, 4). Some of them have been proven to be contributed to the onset of AD, such as HLA-DRB1, HLA-DQA1, and HLA-DQB1 (McGeer et al., 1988; Payton et al., 2006; Mansouri et al., 2015; Zhang et al., 2022b).
To obtain more reliable AD risk factors, the DEGs in AD were used to conduct the Fisher test on 415 AD-associated risk genes. The 34 common genes were obtained in the brain regions, and two common genes in the blood. The deposition of β-amyloid protein in the brain plays a crucial role in the pathogenesis of AD. Generally, under normal physiological conditions, the production and clearance of Aβ maintain a dynamic balance. However, under pathological conditions, Aβ production increases or clearance decreases, disrupting the balance and leading to excessive deposition of Aβ in the brain. This, in turn, triggers a series of pathological processes, such as mitochondrial dysfunction, oxidative stress, and neurofibrillary tangles (Selkoe, 1993). Here, several risk factors that affect Aβ deposition and formation were also identified, including CNTNAP1, CR1, CTSH, and IL34 (Table 2). CR1, encoding a type-I transmembrane glycoprotein, is one of the most important risk genes for late-onset AD, playing multiple roles in the onset of AD, such as Aβ clearance, neuroinflammation, and tauopathy (Zhu et al., 2015). Additionally, autophagy disorders were proved to be associated with AD, and our study also identified autophagy-related genes, including ERRB2 and RAB8B, which have been preliminarily confirmed to be associated with the pathogenesis of AD (Wang et al., 2017; Martinez et al., 2023). Growing evidence from clinical and pathological studies indicates the important relationships between the ongoing deterioration of brain cholesterol metabolic disturbance and AD pathophysiology (Vestergaard et al., 2010). APOC1 and APOC2, both belonging to the apolipoprotein family, and their gene polymorphisms have been reported to be associating with the onset of AD by many studies (Cervantes et al., 2011; Kulminski et al., 2022). In summary, 21 out of 36 highly reliable genes can participate in the pathology of AD through various mechanisms.
More importantly, we have identified many new AD risk factors (Table 2), and their specific involvement in AD has not been clearly reported. By analyzing the relationship between highly reliable genes and enrichment pathways, clues were provided for future research on the functions of these newly identified genes in AD (Figure 5). Most of these pathways have been either been reported or potentially associated with the pathology of AD, including multiple innate and adaptive immune pathways, cell death regulation, DNA repair, cell adhesion, lipoprotein metabolism, protein endocytosis and exocytosis. Therefore, these pathways deserve special attention in future AD research. The newly identified AD risk genes may participate in the occurrence and development of AD through these pathways. For instance, several studies have reported abnormal structure and function of mitochondria in the AD brain, leading to abnormal energy metabolism (Kerr et al., 2017). This abnormal energy metabolism can impact the synaptic plasticity of neurons, thereby affecting memory and learning. CISD1 and DNA2 were reported to be related to mitochondrial function (Geldenhuys et al., 2014; Ding and Liu, 2015), but further detailed and in-depth research is needed to understand their involvement in AD. This implies that significant efforts are still required to explore and investigate the risk factors of AD.
Although many AD-associated targets have been identified by GWAS, high-throughput sequencing, molecular epidemiology, and other methods in the past few decades, it is still a drop in the bucket to improve the awareness of the pathogenesis of AD and the diagnosis and treatment of AD. In this study, we utilized the improved TWAS method, MR-JTI, to integrate and analyze LD score, GTEx eQTL data, and GWAS summary statistic data. A total of 415 AD-associated genes were identified, and the 36 more reliable AD risk was further confirmed by using 11 AD-associated datasets and the Fisher test. The identification of these genes is not only the verification of reported AD-associated genes but also provides new potential AD biomarkers for follow-up research.
Data availability statement
The raw data used in this study come from public databases, and the link and accession numbers are listed in the text. The resulting data presented in this study are included in the article/Supplementary material.
Author contributions
YW and YT: conceptualization. YW, PL, T-HL, and YT: methodology. PL, XL, YH, and LH: data collection. YW, PL, YH, and YT: data analysis. YW, PL, XL, and YT: data curation. YW, PL, XL, and LH: visualization. YW, PL, AG, T-HL, and YT: drafted the manuscript. All authors contributed to the article and approved the submitted version.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnagi.2023.1183119/full#supplementary-material
Footnotes
- ^ https://www.ebi.ac.uk/gwas/
- ^ https://www.gtexportal.org/home/
- ^ https://zenodo.org/record/3842289#.Y9iX_XZByyx
- ^ https://ctg.cncr.nl/software/MAGMA/ref_data/
- ^ http://www.ncbi.nlm.nih.gov/geo
- ^ https://david.ncifcrf.gov/tools.jsp
- ^ https://www.genome.jp/kegg/
- ^ http://geneontology.org
- ^ https://www.innatedb.com/
- ^ https://www.immport.org/home
References
Alzheimer’s Association (2021). 2021 Alzheimer’s disease facts and figures. Alzheimers Dement. 17, 327–406. doi: 10.1002/alz.12328
Antonell, A., Llado, A., Altirriba, J., Botta-Orfila, T., Balasa, M., Fernandez, M., et al. (2013). A preliminary study of the whole-genome expression profile of sporadic and monogenic early-onset Alzheimer’s disease. Neurobiol. Aging 34, 1772–1778. doi: 10.1016/j.neurobiolaging.2012.12.026
Barbeira, A. N., Dickinson, S. P., Bonazzola, R., Zheng, J., Wheeler, H. E., Torres, J. M., et al. (2018). Exploring the phenotypic consequences of tissue specific gene expression variation inferred from GWAS summary statistics. Nat. Commun. 9:1825.
Bellenguez, C., Kucukali, F., Jansen, I. E., Kleineidam, L., Moreno-Grau, S., Amin, N., et al. (2022). New insights into the genetic etiology of Alzheimer’s disease and related dementias. Nat. Genet. 54, 412–436. doi: 10.1038/s41588-022-01024-z
Bhattacharya, S., Dunn, P., Thomas, C. G., Smith, B., Schaefer, H., Chen, J., et al. (2018). ImmPort, toward repurposing of open access immunological assay data for translational and clinical research. Sci. Data 5:180015. doi: 10.1038/sdata.2018.15
Blalock, E. M., Geddes, J. W., Chen, K. C., Porter, N. M., Markesbery, W. R., and Landfield, P. W. (2004). Incipient Alzheimer’s disease: microarray correlation analyses reveal major transcriptional and tumor suppressor responses. Proc. Natl. Acad. Sci. U S A. 101, 2173–2178. doi: 10.1073/pnas.0308512100
Blennow, K., and Zetterberg, H. (2018). Biomarkers for Alzheimer’s disease: current status and prospects for the future. J. Intern. Med. 284, 643–663. doi: 10.1111/joim.12816
Breuer, K., Foroushani, A. K., Laird, M. R., Chen, C., Sribnaia, A., Lo, R., et al. (2013). InnateDB: systems biology of innate immunity and beyond–recent updates and continuing curation. Nucleic Acids Res. 41, D1228–D1233. doi: 10.1093/nar/gks1147
Cervantes, S., Samaranch, L., Vidal-Taboada, J. M., Lamet, I., Bullido, M. J., Frank-Garcia, A., et al. (2011). Genetic variation in APOE cluster region and Alzheimer’s disease risk. Neurobiol. Aging 32, 2107, e7-e17. doi: 10.1016/j.neurobiolaging.2011.05.023
Chan, S. L., Kim, W. S., Kwok, J. B., Hill, A. F., Cappai, R., Rye, K. A., et al. (2008). ATP-binding cassette transporter A7 regulates processing of amyloid precursor protein in vitro. J. Neurochem. 106, 793–804. doi: 10.1111/j.1471-4159.2008.05433.x
Consortium, Laboratory, Data Analysis and Coordinating Center (LDACC)—Analysis Working Group, Statistical Methods groups—Analysis Working Group, Enhancing GTEx (eGTEx) groups, NIH Common Fund, and NIH/NCI. (2017). Genetic effects on gene expression across human tissues. Nature 550, 204–213. doi: 10.1038/nature24277
De Roeck, A., Van Broeckhoven, C., and Sleegers, K. (2019). The role of ABCA7 in Alzheimer’s disease: evidence from genomics, transcriptomics and methylomics. Acta Neuropathol. 138, 201–220. doi: 10.1007/s00401-019-01994-1
Ding, L., and Liu, Y. (2015). Borrowing nuclear DNA helicases to protect mitochondrial DNA. Int. J. Mol. Sci. 16, 10870–10887. doi: 10.3390/ijms160510870
Dunckley, T., Beach, T. G., Ramsey, K. E., Grover, A., Mastroeni, D., Walker, D. G., et al. (2006). Gene expression correlates of neurofibrillary tangles in Alzheimer’s disease. Neurobiol. Aging 27, 1359–1371. doi: 10.1016/j.neurobiolaging.2005.08.013
Gamazon, E. R., Wheeler, H. E., Shah, K. P., Mozaffari, S. V., Aquino-Michaels, K., Carroll, R. J., et al. (2015). A gene-based association method for mapping traits using reference transcriptome data. Nat. Genet. 47, 1091–1098. doi: 10.1038/ng.3367
Geldenhuys, W. J., Leeper, T. C., and Carroll, R. T. (2014). mitoNEET as a novel drug target for mitochondrial dysfunction. Drug Discov. Today 19, 1601–1606. doi: 10.1016/j.drudis.2014.05.001
Gockley, J., Montgomery, K. S., Poehlman, W. L., Wiley, J. C., Liu, Y., Gerasimov, E., et al. (2021). Multi-tissue neocortical transcriptome-wide association study implicates 8 genes across 6 genomic loci in Alzheimer’s disease. Genome Med. 13:76. doi: 10.1186/s13073-021-00890-2
Hansson, O., Seibyl, J., Stomrud, E., Zetterberg, H., Trojanowski, J. Q., Bittner, T., et al. (2018). CSF biomarkers of Alzheimer’s disease concord with amyloid-beta PET and predict clinical progression: a study of fully automated immunoassays in BioFINDER and ADNI cohorts. Alzheimers Dement. 14, 1470–1481. doi: 10.1016/j.jalz.2018.01.010
Hokama, M., Oka, S., Leon, J., Ninomiya, T., Honda, H., Sasaki, K., et al. (2014). Altered expression of diabetes-related genes in Alzheimer’s disease brains: the Hisayama study. Cereb. Cortex 24, 2476–2488. doi: 10.1093/cercor/bht101
Hur, J. Y., Teranishi, Y., Kihara, T., Yamamoto, N. G., Inoue, M., Hosia, W., et al. (2012). Identification of novel gamma-secretase-associated proteins in detergent-resistant membranes from brain. J. Biol. Chem. 287, 11991–12005. doi: 10.1074/jbc.M111.246074
Jansen, I. E., Savage, J. E., Watanabe, K., Bryois, J., Williams, D. M., Steinberg, S., et al. (2019). Genome-wide meta-analysis identifies new loci and functional pathways influencing Alzheimer’s disease risk. Nat. Genet. 51, 404–413. doi: 10.1038/s41588-018-0311-9
Jin, S., Wang, X., Xiang, X., Wu, Y., Hu, J., Li, Y., et al. (2021). Inhibition of GPR17 with cangrelor improves cognitive impairment and synaptic deficits induced by Abeta(1-42) through Nrf2/HO-1 and NF-kappaB signaling pathway in mice. Int. Immunopharmacol. 101(Pt B):108335. doi: 10.1016/j.intimp.2021.108335
Kant, S., Stopa, E. G., Johanson, C. E., Baird, A., and Silverberg, G. D. (2018). Choroid plexus genes for CSF production and brain homeostasis are altered in Alzheimer’s disease. Fluids Barriers CNS 15:34. doi: 10.1186/s12987-018-0120-7
Karch, C. M., and Goate, A. M. (2015). Alzheimer’s disease risk genes and mechanisms of disease pathogenesis. Biol. Psychiatry 77, 43–51. doi: 10.1016/j.biopsych.2014.05.006
Kerr, J. S., Adriaanse, B. A., Greig, N. H., Mattson, M. P., Cader, M. Z., Bohr, V. A., et al. (2017). Mitophagy and Alzheimer’s disease: cellular and molecular mechanisms. Trends Neurosci. 40, 151–166. doi: 10.1016/j.tins.2017.01.002
Khan, S., Barve, K. H., and Kumar, M. S. (2020). Recent advancements in pathogenesis, diagnostics and treatment of Alzheimer’s disease. Curr. Neuropharmacol. 18, 1106–1125. doi: 10.2174/1570159X18666200528142429
Kibinge, N. K., Relton, C. L., Gaunt, T. R., and Richardson, T. G. (2020). Characterizing the causal pathway for genetic variants associated with neurological phenotypes using human brain-derived proteome data. Am. J. Hum. Genet. 106, 885–892. doi: 10.1016/j.ajhg.2020.04.007
Kulminski, A. M., Jain-Washburn, E., Loiko, E., Loika, Y., Feng, F., Culminskaya, I., et al. (2022). Associations of the APOE epsilon2 and epsilon4 alleles and polygenic profiles comprising APOE-TOMM40-APOC1 variants with Alzheimer’s disease biomarkers. Aging 14, 9782–9804. doi: 10.18632/aging.204384
Li, Y., Laws, S. M., Miles, L. A., Wiley, J. S., Huang, X., Masters, C. L., et al. (2021). Genomics of Alzheimer’s disease implicates the innate and adaptive immune systems. Cell Mol. Life Sci. 78, 7397–7426. doi: 10.1007/s00018-021-03986-5
Li, Y., Shi, H., Chen, T., Xue, J., Wang, C., Peng, M., et al. (2022). Establishing a competing endogenous RNA (ceRNA)-immunoregulatory network associated with the progression of Alzheimer’s disease. Ann. Transl. Med. 10:65. doi: 10.21037/atm-21-6762
Li, Y., Xu, M., Xiang, B. L., Li, X., Zhang, D. F., Zhao, H., et al. (2023). Functional genomics identify causal variant underlying the protective CTSH locus for Alzheimer’s disease. Neuropsychopharmacology Online ahead of print. doi: 10.1038/s41386-023-01542-2
Liu, L., Zeng, P., Xue, F., Yuan, Z., and Zhou, X. (2021). Multi-trait transcriptome-wide association studies with probabilistic Mendelian randomization. Am. J. Hum. Genet. 108, 240–256. doi: 10.1016/j.ajhg.2020.12.006
Logue, M. W., Lancour, D., Farrell, J., Simkina, I., Fallin, M. D., Lunetta, K. L., et al. (2018). Targeted sequencing of Alzheimer disease genes in african americans implicates novel risk variants. Front. Neurosci. 12:592. doi: 10.3389/fnins.2018.00592
Lu, R. C., Yang, W., Tan, L., Sun, F. R., Tan, M. S., Zhang, W., et al. (2017). Association of HLA-DRB1 polymorphism with Alzheimer’s disease: a replication and meta-analysis. Oncotarget 8, 93219–93226. doi: 10.18632/oncotarget.21479
Luningham, J. M., Chen, J., Tang, S., De Jager, P. L., Bennett, D. A., Buchman, A. S., et al. (2020). Bayesian genome-wide TWAS method to leverage both cis- and trans-eQTL information through summary statistics. Am. J. Hum. Genet. 107, 714–726. doi: 10.1016/j.ajhg.2020.08.022
Lyons, C. E., Zhou, X., Razzoli, M., Chen, M., Xia, W., Ashe, K., et al. (2022). Lifelong chronic psychosocial stress induces a proteomic signature of Alzheimer’s disease in wildtype mice. Eur. J. Neurosci. 55, 2971–2985. doi: 10.1111/ejn.15329
Maes, O. C., Schipper, H. M., Chertkow, H. M., and Wang, E. (2009). Methodology for discovery of Alzheimer’s disease blood-based biomarkers. J. Gerontol. A Biol. Sci. Med. Sci. 64, 636–645. doi: 10.1093/gerona/glp045
Maes, O. C., Schipper, H. M., Chong, G., Chertkow, H. M., and Wang, E. (2010). A GSTM3 polymorphism associated with an etiopathogenetic mechanism in Alzheimer disease. Neurobiol. Aging 31, 34–45. doi: 10.1016/j.neurobiolaging.2008.03.007
Maes, O. C., Xu, S., Yu, B., Chertkow, H. M., Wang, E., and Schipper, H. M. (2007). Transcriptional profiling of Alzheimer blood mononuclear cells by microarray. Neurobiol. Aging 28, 1795–1809. doi: 10.1016/j.neurobiolaging.2006.08.004
Mansouri, L., Messalmani, M., Klai, S., Bedoui, I., Derbali, H., Gritli, N., et al. (2015). Association of HLA-DR/DQ polymorphism with Alzheimer’s disease. Am. J. Med. Sci. 349, 334–337. doi: 10.1097/MAJ.0000000000000416
Martinez, R. A. S., Pinky, P. D., Harlan, B. A., and Brewer, G. J. (2023). GTP energy dependence of endocytosis and autophagy in the aging brain and Alzheimer’s disease. Geroscience 45, 757–780. doi: 10.1007/s11357-022-00717-x
Masters, C. L., and Selkoe, D. J. (2012). Biochemistry of amyloid beta-protein and amyloid deposits in Alzheimer disease. Cold Spring Harb. Perspect. Med. 2:a006262. doi: 10.1101/cshperspect.a006262
McGeer, P. L., Itagaki, S., Boyes, B. E., and McGeer, E. G. (1988). Reactive microglia are positive for HLA-DR in the substantia nigra of Parkinson’s and Alzheimer’s disease brains. Neurology 38, 1285–1291. doi: 10.1212/WNL.38.8.1285
Nagpal, S., Meng, X., Epstein, M. P., Tsoi, L. C., Patrick, M., Gibson, G., et al. (2019). TIGAR: an improved bayesian tool for transcriptomic data imputation enhances gene mapping of complex traits. Am. J. Hum. Genet. 105, 258–266. doi: 10.1016/j.ajhg.2019.05.018
Nott, A., Holtman, I. R., Coufal, N. G., Schlachetzki, J. C. M., Yu, M., Hu, R., et al. (2019). Brain cell type-specific enhancer-promoter interactome maps and disease-risk association. Science 366, 1134–1139. doi: 10.1126/science.aay0793
Nunez-Iglesias, J., Liu, C. C., Morgan, T. E., Finch, C. E., and Zhou, X. J. (2010). Joint genome-wide profiling of miRNA and mRNA expression in Alzheimer’s disease cortex reveals altered miRNA regulation. PLoS One 5:e8898. doi: 10.1371/journal.pone.0008898
Oba, T., Saito, T., Asada, A., Shimizu, S., Iijima, K. M., and Ando, K. (2020). Microtubule affinity-regulating kinase 4 with an Alzheimer’s disease-related mutation promotes tau accumulation and exacerbates neurodegeneration. J. Biol. Chem. 295, 17138–17147. doi: 10.1074/jbc.RA120.014420
Park, Y. H., Pyun, J. M., Hodges, A., Jang, J. W., Bice, P. J., Kim, S., et al. (2021). Dysregulated expression levels of APH1B in peripheral blood are associated with brain atrophy and amyloid-beta deposition in Alzheimer’s disease. Alzheimers Res. Ther. 13:183. doi: 10.1186/s13195-021-00919-z
Patel, D., Zhang, X., Farrell, J. J., Chung, J., Stein, T. D., Lunetta, K. L., et al. (2021). Cell-type-specific expression quantitative trait loci associated with Alzheimer disease in blood and brain tissue. Transl. Psychiatry 11:250. doi: 10.1038/s41398-021-01373-z
Patel, H., Hodges, A. K., Curtis, C., Lee, S. H., Troakes, C., Dobson, R. J. B., et al. (2019). Transcriptomic analysis of probable asymptomatic and symptomatic alzheimer brains. Brain Behav. Immun. 80, 644–656. doi: 10.1016/j.bbi.2019.05.009
Payton, A., van den Boogerd, E., Davidson, Y., Gibbons, L., Ollier, W., Rabbitt, P., et al. (2006). Influence and interactions of cathepsin D, HLA-DRB1 and APOE on cognitive abilities in an older non-demented population. Genes Brain Behav. 5(Suppl. 1), 23–31. doi: 10.1111/j.1601-183X.2006.00191.x
Prikas, E., Paric, E., Asih, P. R., Stefanoska, K., Stefen, H., Fath, T., et al. (2022). Tau target identification reveals NSF-dependent effects on AMPA receptor trafficking and memory formation. EMBO J. 41:e10242. doi: 10.15252/embj.2021110242
Rodriguez-Fontenla, C., and Carracedo, A. (2021). UTMOST, a single and cross-tissue TWAS (Transcriptome wide association study), reveals new ASD (autism spectrum disorder) associated genes. Transl. Psychiatry 11:256. doi: 10.1038/s41398-021-01378-8
Scheltens, P., De Strooper, B., Kivipelto, M., Holstege, H., Chetelat, G., Teunissen, C. E., et al. (2021). Alzheimer’s disease. Lancet 397, 1577–1590. doi: 10.1016/S0140-6736(20)32205-4
Schwartzentruber, J., Cooper, S., Liu, J. Z., Barrio-Hernandez, I., Bello, E., Kumasaka, N., et al. (2021). Genome-wide meta-analysis, fine-mapping and integrative prioritization implicate new Alzheimer’s disease risk genes. Nat. Genet. 53, 392–402. doi: 10.1038/s41588-020-00776-w
Selkoe, D. J. (1993). Physiological production of the beta-amyloid protein and the mechanism of Alzheimer’s disease. Trends Neurosci. 16, 403–409. doi: 10.1016/0166-2236(93)90008-A
Seripa, D., Panza, F., Paroni, G., D’Onofrio, G., Bisceglia, P., Gravina, C., et al. (2018). Role of CLU, PICALM, and TNK1 genotypes in aging with and without Alzheimer’s disease. Mol. Neurobiol. 55, 4333–4344. doi: 10.1007/s12035-017-0547-x
Shao, Y., Shaw, M., Todd, K., Khrestian, M., D’Aleo, G., Barnard, P. J., et al. (2018). DNA methylation of TOMM40-APOE-APOC2 in Alzheimer’s disease. J. Hum. Genet. 63, 459–471. doi: 10.1038/s10038-017-0393-8
Shen, L., Liao, L., Chen, C., Guo, Y., Song, D., Wang, Y., et al. (2017). Proteomics Analysis of Blood Serums from Alzheimer’s Disease Patients Using iTRAQ Labeling Technology. J. Alzheimers Dis. 56, 361–378. doi: 10.3233/JAD-160913
Shen, R., Zhao, X., He, L., Ding, Y., Xu, W., Lin, S., et al. (2020). Upregulation of RIN3 induces endosomal dysfunction in Alzheimer’s disease. Transl. Neurodegener. 9:26. doi: 10.1186/s40035-020-00206-1
Stopa, E. G., Tanis, K. Q., Miller, M. C., Nikonova, E. V., Podtelezhnikov, A. A., Finney, E. M., et al. (2018). Comparative transcriptomics of choroid plexus in Alzheimer’s disease, frontotemporal dementia and Huntington’s disease: implications for CSF homeostasis. Fluids Barriers CNS 15:18. doi: 10.1186/s12987-018-0102-9
Tan, M. G., Chua, W. T., Esiri, M. M., Smith, A. D., Vinters, H. V., and Lai, M. K. (2010). Genome wide profiling of altered gene expression in the neocortex of Alzheimer’s disease. J. Neurosci. Res. 88, 1157–1169. doi: 10.1002/jnr.22290
Tang, S., Buchman, A. S., De Jager, P. L., Bennett, D. A., Epstein, M. P., and Yang, J. (2021). Novel variance-component TWAS method for studying complex human diseases with applications to Alzheimer’s dementia. PLoS Genet. 17:e1009482. doi: 10.1371/journal.pgen.1009482
Vestergaard, M., Hamada, T., Morita, M., and Takagi, M. (2010). Cholesterol, lipids, amyloid Beta, and Alzheimer’s. Curr. Alzheimer Res. 7, 262–270. doi: 10.2174/156720510791050821
Wang, B. J., Her, G. M., Hu, M. K., Chen, Y. W., Tung, Y. T., Wu, P. Y., et al. (2017). ErbB2 regulates autophagic flux to modulate the proteostasis of APP-CTFs in Alzheimer’s disease. Proc. Natl. Acad. Sci. U S A. 114, E3129–E3138. doi: 10.1073/pnas.1618804114
Wang, H., Devadoss, D., Nair, M., Chand, H. S., and Lakshmana, M. K. (2022). Novel Alzheimer risk factor IQ motif containing protein K is abundantly expressed in the brain and is markedly increased in patients with Alzheimer’s disease. Front. Cell Neurosci. 16:954071. doi: 10.3389/fncel.2022.954071
Waseem, R., Anwar, S., Khan, S., Shamsi, A., Hassan, M. I., Anjum, F., et al. (2021). MAP/Microtubule affinity regulating kinase 4 inhibitory potential of irisin: a new therapeutic strategy to combat cancer and Alzheimer’s disease. Int. J. Mol. Sci. 22:10986. doi: 10.3390/ijms222010986
Wildsmith, K. R., Holley, M., Savage, J. C., Skerrett, R., and Landreth, G. E. (2013). Evidence for impaired amyloid beta clearance in Alzheimer’s disease. Alzheimers Res. Ther. 5:33. doi: 10.1186/alzrt187
Yang, C., Wan, X., Lin, X., Chen, M., Zhou, X., and Liu, J. (2019). CoMM: a collaborative mixed model to dissecting genetic contributions to complex traits by leveraging regulatory information. Bioinformatics 35, 1644–1652. doi: 10.1093/bioinformatics/bty865
Yokoyama, J. S., Wang, Y., Schork, A. J., Thompson, W. K., Karch, C. M., Cruchaga, C., et al. (2016). Association between genetic traits for immune-mediated diseases and Alzheimer disease. JAMA Neurol. 73, 691–697. doi: 10.1001/jamaneurol.2016.0150
Yu, L., Boyle, P. A., Wingo, A. P., Yang, J., Wang, T., Buchman, A. S., et al. (2022). Neuropathologic correlates of human cortical proteins in Alzheimer disease and related dementias. Neurology 98, e1031–e1039. doi: 10.1212/WNL.0000000000013252
Yuan, Z., Zhu, H., Zeng, P., Yang, S., Sun, S., Yang, C., et al. (2020). Testing and controlling for horizontal pleiotropy with probabilistic Mendelian randomization in transcriptome-wide association studies. Nat. Commun. 11:3861. doi: 10.1038/s41467-020-17668-6
Zhang, X., Yu, W., Cao, X., Wang, Y., Zhu, C., and Guan, J. (2022a). Identification of serum biomarkers in patients with Alzheimer’s disease by 2D-DIGE proteomics. Gerontology 68, 686–698. doi: 10.1159/000520961
Zhang, X., Zou, M., Wu, Y., Jiang, D., Wu, T., Zhao, Y., et al. (2022b). Regulation of the late onset Alzheimer’s disease associated HLA-DQA1/DRB1 expression. Am. J. Alzheimers Dis. Other Demen. 37:15333175221085066. doi: 10.1177/15333175221085066
Zhou, D., Jiang, Y., Zhong, X., Cox, N. J., Liu, C., and Gamazon, E. R. (2020). A unified framework for joint-tissue transcriptome-wide association and Mendelian randomization analysis. Nat. Genet. 52, 1239–1246. doi: 10.1038/s41588-020-0706-2
Zhou, S., Ma, G., Luo, H., Shan, S., Xiong, J., and Cheng, G. (2022). Identification of five potential predictive biomarkers for Alzheimer’s disease by integrating the unified test for molecular signatures and weighted gene co-expression network analysis. J. Gerontol. A Biol. Sci. Med. Sci. 78, 653–658. doi: 10.1093/gerona/glac179
Zhou, X., Carbonetto, P., and Stephens, M. (2013). Polygenic modeling with bayesian sparse linear mixed models. PLoS Genet. 9:e1003264. doi: 10.1371/journal.pgen.1003264
Zhu, X. C., Yu, J. T., Jiang, T., Wang, P., Cao, L., and Tan, L. (2015). CR1 in Alzheimer’s disease. Mol. Neurobiol. 51, 753–765. doi: 10.1007/s12035-014-8723-8
Keywords: GWAS, TWAS, Mendelian Randomization, eQTL, Alzheimer’s disease
Citation: Wang YH, Luo PP, Geng AY, Li X, Liu T-H, He YJ, Huang L and Tang YQ (2023) Identification of highly reliable risk genes for Alzheimer’s disease through joint-tissue integrative analysis. Front. Aging Neurosci. 15:1183119. doi: 10.3389/fnagi.2023.1183119
Received: 11 March 2023; Accepted: 30 May 2023;
Published: 21 June 2023.
Edited by:
Emilia Vitale, National Research Council (CNR), ItalyReviewed by:
Chunyu Li, Sichuan University, ChinaEva Bagyinszky, Gachon University, Republic of Korea
Copyright © 2023 Wang, Luo, Geng, Li, Liu, He, Huang and Tang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ya Qin Tang, yqtang@cqmu.edu.cn
†These authors have contributed equally to this work