Original Research ARTICLE
Shared Biological Pathways Between Alzheimer’s Disease and Ischemic Stroke
- 1Department of Neurology, Tianjin Neurological Institute, Tianjin Medical University General Hospital, Tianjin, China
- 2Key Laboratory of Post-neurotrauma Neuro-repair and Regeneration in Central Nervous System, Tianjin Neurological Institute, Ministry of Education and Tianjin City, Tianjin, China
Alzheimer’s disease (AD) and ischemic stroke (IS) are an immense socioeconomic burden worldwide. There is a possibility that shared genetic factors lead to their links at epidemiological and pathophysiological levels. Although recent genome-wide association studies (GWAS) have provided profound insights into the genetics of AD and IS, no shared genetic variants have been identified to date. This prompted us to initiate this study, which sought to identify shared pathways linking AD and IS. We took advantage of large-scale GWAS summary data of AD (17,008 AD cases and 37,154 controls) and IS (10,307 cases and 19,326 controls) to conduct pathway analyses using genetic pathways from multiple well-studied databases, including GO, KEGG, PANTHER, Reactome, and Wikipathways. Collectively, we discovered that AD and IS shared 179 GO categories (56 biological processes, 95 cellular components, and 28 molecular functions); and the following pathways: six KEGG pathways; two PANTHER pathways; four Reactome pathways; and one in Wikipathways pathway. The more fine-grained GO terms were mainly summarized into different functional categories: transcriptional and post-transcriptional regulation, synapse, endocytic membrane traffic through the endosomal system, signaling transduction, immune process, multi-organism process, protein catabolic metabolism, and cell adhesion. The shared pathways were roughly classified into three categories: immune system; cancer (NSCLC and glioma); and signal transduction pathways involving the cadherin signaling pathway, Wnt signaling pathway, G-protein signaling and downstream signaling mediated by phosphoinositides (PIPs). The majority of these common pathways linked to both AD and IS were supported by convincing evidence from the literature. In conclusion, our findings contribute to a better understanding of common biological mechanisms underlying AD and IS and serve as a guide to direct future research.
Alzheimer’s disease (AD) is the most common cause of dementia in the elderly. It is characterized pathologically by extracellular deposits of amyloid-β peptide (Aβ) and intracellular neurofibrillary tangles (NFTs) containing tau protein (Cuyvers and Sleegers, 2016). Ischemic stroke (IS) is a major health-threatening cerebrovascular disease with severe complications, such as post-stroke infection, physical disability, and cognitive deficits (Donnan et al., 2008). These two neurological diseases have immense socioeconomic impact worldwide (Norrving and Kissela, 2013; Alzheimer’sAssociation, 2016).
Mounting evidence suggests that there is a potential correspondence between AD and IS. Firstly, epidemiological studies have revealed that AD is a contributing factor to the development of IS (Chi et al., 2013; Tolppanen et al., 2013), and vice versa (Gamaldo et al., 2006). Secondly, AD and IS have common risk factors (e.g., hypertension, diabetes, obesity, and hyperlipidemia) (de Bruijn and Ikram, 2014; Boehme et al., 2017). Thirdly, recent evidence has indicated that brain ischemia can promote the development of AD by inducing β-/γ-secretase-mediated Aβ accumulation and tau protein gene alterations (Salminen et al., 2017; Pluta et al., 2018). Fourth, neuroinflammation elicited by the immune system is thought to play essential roles in the development and progression of both AD and IS (Liu Y.H. et al., 2013; Anrather and Iadecola, 2016). Lastly, research shows that abnormal tau protein also plays crucial roles in IS (Tuo et al., 2017). Together, these findings support the hypothesis that shared genetic risk factors link AD and IS at the epidemiological and pathophysiological levels.
Recent genome-wide association studies (GWASs) have provided profound insights into the complex genetic architecture of AD (Shen et al., 2015; Li et al., 2016; Liu et al., 2016; Zhang et al., 2016; Jun et al., 2017; Sims et al., 2017) and IS (Traylor et al., 2012; Williams et al., 2013; Kilarski et al., 2014; Malik et al., 2016; Network et al., 2016; Crawford et al., 2018). Most recently, 29 risk loci for AD (Jansen et al., 2018) and 28 loci associated with IS and its subtypes [i.e., large vessel disease (LVD), cardioembolic stroke (CE), and small vessel disease (SVD)] (Malik et al., 2018) have been identified. Despite these findings, no significant genetic variants shared by AD and IS or its subtypes were discovered. From the perspective of pathways, however, four gene ontology (GO) categories jointly associated with AD and SVD were identified (Traylor et al., 2016). Moreover, single-disease pathway analysis demonstrates that AD and IS share common pathways that involve natural killer (NK) cells, i.e., NK cell mediated cytotoxicity in AD (Lambert et al., 2010; Liu et al., 2014; Chen et al., 2016; Jiang et al., 2017) and NK cell signaling in IS (Malik et al., 2016). These findings imply that the link between AD and IS may have to do with shared genetic signals at the pathway level.
Complex diseases like AD and IS are mostly driven by the joint interactions of associated genes affected by large proportions of SNPs well below genome-wide significance, and the crosstalk of regulatory pathways (Furlong, 2013). Based on this concept, pathway-based analysis has been an effective strategy to investigate the potential mechanisms of complex disease (Luo et al., 2010; Jin et al., 2014), thus is widely used to unravel either single disease etiology or pleiotropism of clinically distinct diseases (Lesnick et al., 2007; Liu G. et al., 2012; Liu Y. et al., 2013; Liu et al., 2014; Bao et al., 2015; Xiang et al., 2015). Taken together, these considerations prompted us to perform pathway analyses using AD and IS GWAS data to further dissect their common molecular mechanisms.
Materials and Methods
We obtained AD GWAS summary data from the International Genomics of Alzheimer’s Project (IGAP). IGAP performed a large, two-stage joint association analysis of AD on individuals of European descent (Lambert et al., 2013). In stage 1 (discovery stage), 17,008 AD cases and 37,154 controls from four previously published GWAS samples were included. These were the European Alzheimer’s disease Initiative (EADI), the Alzheimer Disease Genetics Consortium (ADGC), the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium (CHARGE), and the Genetic and Environmental Risk in AD Consortium (GERAD). After quality control procedures, 7,055,881 SNPs remained for further analysis. In stage 2 (replication stage), the top SNPs with a P-value less than 1E-3 in stage 1 were selected for replication in an independent sample (8,572 cases and 11,312 controls). Finally, the results of the two stages were combined in a meta-analysis in order to identify loci reaching genome-wide significance (P < 5.00E-08). In the present study, we used the GWAS summary data from stage 1.
We selected the IS GWAS summary results from the METASTROKE collaboration published by Malik et al (Malik et al., 2016). In the discovery stage, METASTROKE conducted a meta-analysis of 12 case-control GWAS comprising 10,307 Caucasian IS cases and 19,326 Caucasian controls, producing quality-controlled 8.3 million SNPs for further analysis. In the replication stage, the top SNPs (P < 1.00E-05) in the discovery meta-analysis were selected for replication in three independent samples consisting of Caucasian (13,435 cases and 29,269 controls) and South Asian (2,385 cases and 5,193 controls) samples. Ultimately, a transethnic meta-analysis combining the two stages was conducted. Here, we used the IS GWAS summary statistics from the discovery phase.
Firstly, we computed gene-level association statistics through gene-based association tests. Next, gene-based pathway analysis was performed to identify pathways shared between AD and IS.
Gene-Based Testing for AD and IS GWAS Datasets
We performed gene-level analysis implemented in the versatile gene-based association study software (VEGAS). By incorporating the effects of a full set of markers within a gene and correcting for linkage disequilibrium (LD) between SNPs, we were able to calculate gene-based P-values. Compared with other gene-based testing approaches like PLINK, a set-based test based on permutations (Purcell et al., 2007), VEGAS is much more computationally efficient by using simulations from multivariate normal distributions (Liu et al., 2010). In addition, because VEGAS only requires SNP rs IDs and P-values, it is particularly appropriate for GWAS summary results, in which individual genotyped information is not available. We used VEGAS2, a recently updated version of VEGAS that uses LD estimates from the 1000 Genomes Project phase I panel (Mishra and Macgregor, 2015). We chose the following defaults to define the gene boundaries for locating SNPs in VEGAS2: “SNPs within a gene plus SNPs outside of the gene with r2 > 0.8 with SNPs in the gene” (Mishra and Macgregor, 2015). By using this default setting, we were able to capture the effects of regulatory SNPs and minimize non-specificity caused by large boundaries (i.e., ±50 kb) (Mishra and Macgregor, 2015).
Recently, Boyle et al. proposed the hypothesis that most complex traits are largely driven by the vast majority of peripheral genes with smaller effects by propagating through regulatory networks to core genes (Boyle et al., 2017). As such, genes in AD and IS with P < 0.05, rather than with stringent threshold like gene-wide statistical significance, were selected for subsequent pathway analysis.
Pathway-Based Analysis for Shared Disease Pathways
It is likely that genes of AD and IS are involved in one pathway, but in different proportions, such as in the upstream and downstream of the pathway, respectively. Hence, we performed pathway analysis for AD and IS separately by integrating genes for each disease from VEGAS2 (P < 0.05) into known functional annotations using WebGestalt (Wang et al., 2017). Next, their respective enrichment results were overlapped to identify shared pathways jointly associated with AD and IS.
The latest version of WebGestalt 2017 included genetic pathways from GO annotations (Ashburner et al., 2000), KEGG (Kanehisa et al., 2004), PANTHER (Mi et al., 2016), Reactome (Fabregat et al., 2016), and WikiPathways (Pico et al., 2008). We chose over-representation analysis (ORA), a hypergeometric test, to identify the enrichment of the disease-related genes among all the genes in a given pathway. The P-value for observing at least m disease-related genes in a given pathway can be calculated as follows:
where M is the total number of genes of interest that are associated with a given disease, N is the number of reference genes, and n is the number of genes in the pathway. The smaller the P-value is, the more likely the list of disease-related genes of interest will be overrepresented in this pathway. To avoid testing overly narrow or broad pathways, which could result in some false positives, we chose pathways comprising at least 20 and at most 300 genes for analysis (Jiang et al., 2017). The threshold of statistical significance for a pathway associated with AD and IS was P < 0.05 in both diseases.
Meta-Analysis Using Fisher’s Method
Here, Fisher’s method was used to combine the P-values of each shared pathway in AD and IS into one test statistic. For a given pathway, the Fisher equation for the statistic is:
where Pi is the P-value of the pathway in the ith study and k is the number of total studies. x2 follows chi-square distribution with 2k degrees of freedom (Begum et al., 2012).
Gene-Based Testing for AD and IS GWAS Datasets and Validation of VEGAS2
We identified 1915 AD genes and 1288 IS genes with P < 0.05. More detailed results are provided in Supplementary Table S1. To verify the reliability of the VEGAS2 gene-based testing approach, we compared our gene-based association results of AD with previously established AD risk loci identified through conventional GWAS approaches. We replicated 15 genes (APOE, TOMM40, APOC1, EXOC3L2, PVRL2, CR1, HLA-DRB6, EPHA1, CLU, PICALM, ABCA7, CD33, MS4A2, MS4A6A, BIN1) with gene-wide significance (P < 2.35E-6, Bonferroni corrected for 21,244 genes). These genes belong to 12 established risk loci for AD (i.e., APOE/TOMM40/APOC1 locus; EXOC3L2/BLOC1S3/MARK4 locus; PVRL2; CR1; HLA-DRB5-DRB1 region; EPHA1; CLU; PICALM; ABCA7; CD33; MS4A locus; BIN1). In addition, we replicated 11 genes (SORL1, CD2AP, MS4A6E, MS4A4A, MS4A3, HLA-DQB1, HLA-DRB1, HLA-DQA1, PTK2B, CELF1, and SLC24A4) with suggestive association (P < 1.00E-04); these belong to seven known significant loci (SORL1, CD2AP, MS4A locus, HLA-DRB5-DRB1 region, PTK2B, CELF1, and SLC24A4). Thus, the replication of these established findings shows that VEGAS2 can be trusted as a validated tool for our subsequent pathway analysis.
Pathway-Based Analysis for Shared Pathways Between AD and IS
Using the 1915 AD genes and 1288 IS genes having a P < 0.05, pathway analysis for AD and IS were separately conducted followed by overlapping their respective significant enrichment results.
In GO, we identified 179 common GO categories associated with both AD and IS: 56 biological processes, 95 cellular components, and 28 molecular functions (Supplementary Table S2). Here, we paid more attention to the more fine-grained GO terms and weeded out their ancestors to solve the problem of redundancy. Taken as a whole, the remaining GO terms of biological process, cellular component and molecular function mainly fall into categories listed below: transcriptional and post-transcriptional regulation, synapse, endocytic membrane traffic through the endosomal system, signaling transduction, immune process, multi-organism process, protein catabolic metabolism, cell adhesion, and others (Figure 1 and Supplementary Table S2). The detailed information about complete GO annotation results and summarized GO annotation results as well as their relevant GO terms is displayed in Supplementary Table S2.
FIGURE 1. Summarized GO annotation results shared by AD and IS. The complete GO annotation results shared by AD and IS are summarized and reorganized here, with only more fine-grained terms labeled and classified into different functional categories: transcriptional and post-transcriptional regulation, synapse, endocytic membrane traffic, signaling transduction, immune process, multi-organism process, protein catabolic metabolism, cell adhesion, and others, which are shown in nodes. The relevant GO terms of biological process, cellular component, and molecular function are shown in targets. The threshold for statistical significance for a GO term is P < 0.05 in both AD and IS. The edges connecting the nodes and targets are attributed by the value of –log10 P.
In KEGG, we identified six significant pathways shared by AD and IS: NK cell mediated cytotoxicity [hsa04650], Toll-like receptor signaling pathway [hsa04620], non-small cell lung cancer [hsa05223], glioma [hsa05214], phospholipase D signaling pathway [hsa04072], and hepatitis B [hsa05161]. Besides these, we also identified 2 shared significant PANTHER pathways between AD and IS: cadherin signaling pathway [P00012] and Wnt signaling pathway [P00057]. Four significant Reactome pathways shared by AD and IS were immunoregulatory interactions between a lymphoid and a non-lymphoid cell [R-HSA-198933], synthesis of PIPs at the plasma membrane [R-HSA-1660499], PI metabolism [R-HSA-1483255] and cooperation of PDCL (PhLP1) and TRiC/CCT in G-protein beta folding [R-HSA-6814122]. In the WikiPathways dataset, one significant pathway was identified: signaling pathways in glioblastoma [WP2261]. The detailed results are described in Table 1.
Large-scale GWAS of AD and IS have identified a set of risk loi with genome-wide significance and provided deep insights into the genetics of AD and IS (Cuyvers and Sleegers, 2016; Boehme et al., 2017). Despite these discoveries, a large proportion of heritability remains underappreciated (Zhang et al., 2013; Markus and Bevan, 2014). To date, no shared genetic determinants jointly associated with AD and IS have been revealed. However, AD and IS have common biological pathways that may account for their similarities at the epidemiological, neuropathological, and molecular levels. Indeed, our pathway-based association tests using large-scale GWAS summary datasets for AD and IS prove the genetic links between these two diseases at the pathway level.
In order to fully understand the biological features shared by AD and IS, we included in our analyses different functional annotation databases (i.e., GO annotation, KEGG, PANHTER, Reactome, and Wikipathways). GO annotation sets represent the most comprehensive spectrum of functional categories but are mainly computationally inferred, with less than 1% being experimentally confirmed (Ashburner et al., 2000). In contrast, KEGG, PANTHER, and Reactome pathway databases are manually curated from the literature based on experimentally validated evidence (Kanehisa et al., 2004; Fabregat et al., 2016; Mi et al., 2016). Compared to the centrally curated databases above, open community-based Wikipathways allows for broader exploration into the entire biological spectrum (Pico et al., 2008). Altogether, we uncovered 179 GO annotation terms, 6 shared KEGG pathways, 2 PANTHER pathways, 4 Reactome pathways, and one pathway in Wikipathways jointly associated with AD and IS.
Shared GO Annotations
The GO annotation provides a knowledge base for understanding the roles of genes and proteins in cells. Complete GO schema generates a huge number of annotation terms with parent–child relationships, which are non-independently arranged as a directed acyclic graph (DAG). To assist in interpretation, we tried to focus on more fine-grained terms, allowing more specific biological mechanisms to be dissected. These descendants turned out to be relevant to the following functional categories: transcriptional and post-transcriptional regulation, synapse, endocytic membrane traffic through the endosomal system, signaling transduction, immune process, multi-organism process, protein catabolic metabolism, cell adhesion and others (Figure 1). In particular, dysfunction of endocytic membrane dynamics plays a key role in AD pathogenesis by enhancing the processing of amyloid precursor protein (APP) into Aβ (Wu and Yao, 2009; Rajendran and Annaert, 2012). However, the relationship between the endosomal system and IS remains unknown and warrants further extensive investigation.
Shared KEGG Pathways
Together with one Reactome pathway [R-HSA-198933 immunoregulatory interactions between a lymphoid and a non-lymphoid cell], two KEGG pathways are related to the immune system (i.e., Toll-like receptor signaling pathway [hsa04620] and NK cell mediated cytotoxicity [hsa04650]). Toll-like receptor signaling pathway has also been reported to be associated with AD (Liu et al., 2014). It is now widely accepted that Toll-like receptors (TLRs) play essential roles in Aβ-induced microglial inflammatory activation in AD (Liu S. et al., 2012). The persistent inflammatory stimuli fueled by activated microglia in turn leads to a disturbance in microglial clearance of Aβ, eventually causing neuronal degeneration (Heneka et al., 2015). Following IS, damage-associated molecular patterns (DAMPs) released from the intracellular compartment of dying and dead cells activate microglia and perivascular macrophages by pattern recognition receptors (PRRs), especially TLRs. The triggering cascades of innate immunity further exacerbate stroke lesions in the acute phase of IS (Iadecola and Anrather, 2011).
Natural killer cell mediated cytotoxicity has been reported by previous GWAS and brain expression datasets of AD (Lambert et al., 2010; Liu et al., 2014; Chen et al., 2016; Jiang et al., 2017). The association for IS and NK cell signaling was also observed (Malik et al., 2016), and supported by evidence that infiltrating NK cells exacerbate ischemia lesions and increase neuronal death either in experimental IS models or in human IS (Gan et al., 2014). The inflammatory responses in both AD and IS are initiated locally within the central nervous system (CNS) under aberrant local conditions, and presumably elicited by the innate immune system. Therefore, the similarity of their pathophysiological processes may provide a reasonable explanation for our findings that AD and IS share common biological pathways involved in the immune system.
Another two KEGG pathways are related to cancer (i.e., glioma [hsa05214] and non-small cell lung cancer [hsa05223]). Convincing evidence suggests an inverse relationship and significant genetic overlap between AD and cancer (Roe et al., 2005; Driver et al., 2012; Musicco et al., 2013; Ou et al., 2013). For instance, the cellular behavior of proliferative cells in cancer is essentially the opposite of the processes leading to degenerating neurons in AD (Behrens et al., 2009; Driver and Lu, 2010; Driver, 2014). Bioinformatics analyses has unraveled shared genes and signaling pathways connecting AD and glioblastoma (Liu T. et al., 2013). The one significant pathway we found in the Wikipathways dataset (signaling pathways in glioblastoma [WP2261]) also supports this linkage. Though some illustrative cases suggested that glioblastoma should be considered in the etiology of acute IS (Pina et al., 2014; Lasocki and Gaillard, 2016), little is known about the relationship between IS and glioma or other prevalent cancers.
The biological significance of the phospholipase D (PLD) signaling pathway [hsa04072] has been suggested for both AD and IS (Oliveira and Di Paolo, 2010; Stegner et al., 2013; Frohman, 2015). Hepatitis B pathway [hsa05161] has already been revealed in hippocampal CA1 region of AD patients (Mastroeni et al., 2018). Epidemiology studies suggested that hepatitis B virus (HBV) infection was associated with decreased risk of IS, although the underlying mechanism was unclear (Sung et al., 2007; Kuo et al., 2017). Part of the GO terms we deciphered, i.e., positive regulation of viral transcription [GO:0050434], positive regulation of viral life cycle [GO:1903902] and defense response to virus [GO:0051607] appeared to support this linkage.
Shared PANTHER Pathways
In the present study, we uncovered 2 significant PANTHER pathways shared by AD and IS (cadherin signaling pathway [P00012] and Wnt signaling pathway [P00057]). Similar to cadherin signaling pathway, cadherin binding involved in cell-cell adhesion [GO:0098641] was implied in GO. Both enriched PANTHER pathways primarily arose from the Pcdhα gene cluster (Pcdhα) shared by AD and IS in our VEGAS2 gene-based analysis. Protocadherins (Pcdhs) constitute the largest group of the cadherin superfamily of cell-adhesion molecules. Pcdhs are predominantly expressed in the CNS (i.e., neurons, astrocytes, pericytes, choroid plexus epithelial cells, and brain microvascular endothelial cells), and are crucial for neuronal survival, neural circuit assembly, and maintenance of the blood-brain barrier (BBB) (Takeichi, 2007; Dilling et al., 2017). Moreover, the interactions between Pcdhs and Wnt signaling pathway have recently been uncovered (Mah and Weiner, 2017). Wnt signaling pathway plays critical roles in neurogenesis and synaptic transmission and plasticity, whereby its deregulation shows significant linkage to the development of AD (Caruso et al., 2006; De Ferrari et al., 2007; Inestrosa and Varela-Nallar, 2014). For IS, experimental evidence has demonstrated that Wnt signaling activation may exert effective neuroprotection by providing an appropriate compensatory microenvironment for neurogenesis or neuronal survival in focal ischemic injury (Shruster et al., 2012; Wu et al., 2015).
Shared Reactome Pathways
Two shared Reactome pathways are PI (phosphatidylinositol) metabolism [R-HSA-1483255] and synthesis of PIPs (phosphoinositides) at the plasma membrane [R-HSA-1660499], which are connected by a parent–child relationship. Similarly, phosphatidylinositol binding [GO:0035091] was also identified in GO. PIPs consist of seven diverse phosphorylated forms of phosphatidylinositol. Among these, PI(4,5)P2 and PI(3,4,5)P3 occur in relatively high abundance at the plasma membrane (Di Paolo and De Camilli, 2006; Riehle et al., 2013). PI(4,5)P2 primarily serves as a pool for the formation of intracellular second messengers—inositol 1,4,5-trisphosphate (IP3), diacylglycerol (DAG), and PI(3,4,5)P3— to mediate downstream signal cascades [e.g., Ca2+ response, PKC activation and PI(3,4,5)P3/Akt/mTOR pathway] from plasma membrane receptors like GPCRs (Riehle et al., 2013). PIPs play important roles in the nervous system by regulating lipid signaling, receptor signaling and membrane dynamics (e.g., endocytic membrane traffic). Thus, deregulation of these pathways has been implicated in neurological diseases including AD and IS (Waugh, 2015).
Another significant Reactome pathway shared by AD and IS is cooperation of PDCL (PhLP1) and TRiC/CCT in G-protein beta folding [R-HSA-6814122]. Recent studies reveal the interplay of the protein-folding chaperone CCT (also called TRiC) and its co-chaperone PhLP1 in G-protein βγ (Gβγ) dimer formation, which is a particularly critical step in G-protein signaling (Plimpton et al., 2015). Comparably, heterotrimeric G-protein complex [GO:0005834] was also identified in GO. In AD, G-protein signal transduction pathways mediated by G protein-coupled receptors (GPCRs) are strongly related to α-/β-/γ-secretase-mediated APP processing through the endosomal system (Teng et al., 2010; Thathiah and De Strooper, 2011; Heese, 2013).
In summary, we discovered several shared significant pathways jointly associated with AD and IS, which can roughly be classified into three categories: immune system, cancer (NSCLC and glioma), and signal transduction pathways. The latter includes cadherin signaling pathway, Wnt signaling pathway, G-protein signaling, and downstream signaling mediated by PIPs. Convincing evidence in the literature corroborates the majority of these shared pathways.
Despite these interesting and consequential findings, our study had some limitations. Firstly, we were able to analyze only the overall IS summary data. However, IS is genetically heterogeneous, as nearly all its risk loci are specific to individual subtypes. Secondly, the input lists of significant genes for pathway analysis were selected by setting an arbitrary threshold (P < 0.05), which might lead to information loss of some marginally less significant genes (e.g., P = 0.051). Lastly, the statistical model we used for pathway analysis was ORA, a hypergeometric test, which treats each gene equally and does not take into account interdependencies between genes in a pathway. However, the significance of each gene (P-value) can be informative in assigning different weights to input genes in the pathway analysis. In the future, we aim to conduct more robust statistical analysis using multiple pathway analysis methods and multiple datasets to replicate and expand our findings.
Collectively, in our gene-based tests and comprehensive pathway analyses of AD and IS summary GWAS datasets, we discovered several shared novel functional pathways linking AD and IS. By providing an in-depth investigation here of the shared biological mechanism associated with these two neurological disorders, our findings may advance the current understanding of the biology of AD and IS and guide future research on AD or IS in new directions.
PC accomplished by the research. JH applied for the GWAS data of AD and IS. PC completed the data analysis and drafted the manuscript. XM, HL, and WL were responsible for manuscript revision. The final version was approved for submission by all listed authors.
This work was supported by the National Natural Science Foundation of China (Grant Nos. 81571600, 81322018, 81273287, and 81100887 to JH); the Youth Topnotch Talent Support Program (to JH); the National Key Clinical Specialty Construction Project of China (to JH); and the Natural Science Foundation of Tianjin (Grant No. 17JCZDJC35500 to JH).
Conflict of Interest Statement
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.
We thank the International Genomics of Alzheimer’s Project (IGAP) and METASTROKE collaboration for providing summary results data for these analyses.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2018.00605/full#supplementary-material
TABLE S1 | Gene-based testing results for the GWAS datasets of AD and IS.
TABLE S2 | GO annotation results shared by AD and IS.
Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., et al. (2000). Gene ontology: tool for the unification of biology. The gene ontology consortium. Nat. Genet. 25, 25–29. doi: 10.1038/75556
Bao, X., Liu, G., Jiang, Y., Jiang, Q., Liao, M., Feng, R., et al. (2015). Cell adhesion molecule pathway genes are regulated by cis-regulatory SNPs and show significantly altered expression in Alzheimer’s disease brains. Neurobiol. Aging 36, 2904.e1–2904.e7. doi: 10.1016/j.neurobiolaging.2015.06.006
Begum, F., Ghosh, D., Tseng, G. C., and Feingold, E. (2012). Comprehensive literature review and statistical considerations for GWAS meta-analysis. Nucleic Acids Res. 40, 3777–3784. doi: 10.1093/nar/gkr1255
Caruso, A., Motolese, M., Iacovelli, L., Caraci, F., Copani, A., Nicoletti, F., et al. (2006). Inhibition of the canonical Wnt signaling pathway by apolipoprotein E4 in PC12 cells. J. Neurochem. 98, 364–371. doi: 10.1111/j.1471-4159.2006.03867.x
Chen, J., Xie, C., Zhao, Y., Li, Z., Xu, P., and Yao, L. (2016). Gene expression analysis reveals the dysregulation of immune and metabolic pathways in Alzheimer’s disease. Oncotarget 7, 72469–72474. doi: 10.18632/oncotarget.12505
Chi, N. F., Chien, L. N., Ku, H. L., Hu, C. J., and Chiou, H. Y. (2013). Alzheimer disease and risk of stroke: a population-based cohort study. Neurology 80, 705–711. doi: 10.1212/WNL.0b013e31828250af
Crawford, K. M., Gallego-Fabrega, C., Kourkoulis, C., Miyares, L., Marini, S., Flannick, J., et al. (2018). Cerebrovascular disease knowledge portal: an open-access data resource to accelerate genomic discoveries in stroke. Stroke 49, 470–475. doi: 10.1161/STROKEAHA.117.018922
Cuyvers, E., and Sleegers, K. (2016). Genetic variations underlying Alzheimer’s disease: evidence from genome-wide association studies and beyond. Lancet Neurol. 15, 857–868. doi: 10.1016/S1474-4422(16)00127-7
De Ferrari, G. V., Papassotiropoulos, A., Biechele, T., Wavrant De-Vrieze, F., Avila, M. E., Major, M. B., et al. (2007). Common genetic variation within the low-density lipoprotein receptor-related protein 6 and late-onset Alzheimer’s disease. Proc. Natl. Acad. Sci. U.S.A. 104, 9434–9439. doi: 10.1073/pnas.0603523104
Dilling, C., Roewer, N., Förster, C., and Burek, M. (2017). Multiple protocadherins are expressed in brain microvascular endothelial cells and might play a role in tight junction protein regulation. J. Cereb. Blood Flow Metab. 37, 3391–3400. doi: 10.1177/0271678X16688706
Driver, J. A., Beiser, A., Au, R., Kreger, B. E., Splansky, G. L., Kurth, T., et al. (2012). Inverse association between cancer and Alzheimer’s disease: results from the Framingham Heart Study. BMJ 344:e1442. doi: 10.1136/bmj.e1442
Gamaldo, A., Moghekar, A., Kilada, S., Resnick, S. M., Zonderman, A. B., and O’brien, R. (2006). Effect of a clinical stroke on the risk of dementia in a prospective cohort. Neurology 67, 1363–1369. doi: 10.1212/01.wnl.0000240285.89067.3f
Gan, Y., Liu, Q., Wu, W., Yin, J. X., Bai, X. F., Shen, R., et al. (2014). Ischemic neurons recruit natural killer cells that accelerate brain infarction. Proc. Natl. Acad. Sci. U.S.A. 111, 2704–2709. doi: 10.1073/pnas.1315943111
Jansen, I., Savage, J., Watanabe, K., Bryois, J., Williams, D., Steinberg, S., et al. (2018). Genetic meta-analysis identifies 9 novel loci and functional pathways for Alzheimers disease risk. bioRxiv [Preprint]. doi: 10.1101/258533
Jiang, Q., Jin, S., Jiang, Y., Liao, M., Feng, R., Zhang, L., et al. (2017). Alzheimer’s disease variants with the genome-wide significance are significantly enriched in immune pathways and active in immune cells. Mol. Neurobiol. 54, 594–600. doi: 10.1007/s12035-015-9670-8
Jin, L., Zuo, X. Y., Su, W. Y., Zhao, X. L., Yuan, M. Q., Han, L. Z., et al. (2014). Pathway-based analysis tools for complex diseases: a review. Genomics Proteomics Bioinformatics 12, 210–220. doi: 10.1016/j.gpb.2014.10.002
Jun, G. R., Chung, J., Mez, J., Barber, R., Beecham, G. W., Bennett, D. A., et al. (2017). Transethnic genome-wide scan identifies novel Alzheimer’s disease loci. Alzheimers Dement. 13, 727–738. doi: 10.1016/j.jalz.2016.12.012
Kilarski, L. L., Achterberg, S., Devan, W. J., Traylor, M., Malik, R., Lindgren, A., et al. (2014). Meta-analysis in more than 17,900 cases of ischemic stroke reveals a novel association at 12q24.12. Neurology 83, 678–685. doi: 10.1212/WNL.0000000000000707
Kuo, C. S., Chen, Y. T., Hsu, C. Y., Chang, C. C., Chou, R. H., Li, S. Y., et al. (2017). The impact of chronic hepatitis B infection on major adverse cardiovascular events and all-cause mortality in patients with diabetes: a nationwide population-based study from Taiwan. BMJ Open 7:e016179. doi: 10.1136/bmjopen-2017-016179
Lambert, J. C., Grenier-Boley, B., Chouraki, V., Heath, S., Zelenika, D., Fievet, N., et al. (2010). Implication of the immune system in Alzheimer’s disease: evidence from genome-wide pathway analysis. J. Alzheimers Dis. 20, 1107–1118. doi: 10.3233/JAD-2010-100018
Lambert, J. C., Ibrahim-Verbaas, C. A., Harold, D., Naj, A. C., Sims, R., Bellenguez, C., et al. (2013). Meta-analysis of 74,046 individuals identifies 11 new susceptibility loci for Alzheimer’s disease. Nat. Genet. 45, 1452–1458. doi: 10.1038/ng.2802
Lesnick, T. G., Papapetropoulos, S., Mash, D. C., Ffrench-Mullen, J., Shehadeh, L., De Andrade, M., et al. (2007). A genomic pathway approach to a complex disease: axon guidance and Parkinson disease. PLoS Genet. 3:e98. doi: 10.1371/journal.pgen.0030098
Li, Y., Song, D., Jiang, Y., Wang, J., Feng, R., Zhang, L., et al. (2016). CR1 rs3818361 polymorphism contributes to Alzheimer’s disease susceptibility in Chinese population. Mol. Neurobiol. 53, 4054–4059. doi: 10.1007/s12035-015-9343-7
Liu, G., Jiang, Y., Wang, P., Feng, R., Jiang, N., Chen, X., et al. (2012). Cell adhesion molecules contribute to Alzheimer’s disease: multiple pathway analyses of two genome-wide association studies. J. Neurochem. 120, 190–198. doi: 10.1111/j.1471-4159.2011.07547.x
Liu, S., Liu, Y., Hao, W., Wolf, L., Kiliaan, A. J., Penke, B., et al. (2012). TLR2 is a primary receptor for Alzheimer’s amyloid beta peptide to trigger neuroinflammatory activation. J. Immunol. 188, 1098–1107. doi: 10.4049/jimmunol.1101121
Liu, G., Xu, Y., Jiang, Y., Zhang, L., Feng, R., and Jiang, Q. (2016). PICALM rs3851179 variant confers susceptibility to Alzheimer’s disease in Chinese population. Mol. Neurobiol. 54, 3131–3136. doi: 10.1007/s12035-016-9886-2
Liu, G., Yao, L., Liu, J., Jiang, Y., Ma, G., Chen, Z., et al. (2014). Cardiovascular disease contributes to Alzheimer’s disease: evidence from large-scale genome-wide association studies. Neurobiol. Aging 35, 786–792. doi: 10.1016/j.neurobiolaging.2013.10.084
Liu, J. Z., Mcrae, A. F., Nyholt, D. R., Medland, S. E., Wray, N. R., Brown, K. M., et al. (2010). A versatile gene-based test for genome-wide association studies. Am. J. Hum. Genet. 87, 139–145. doi: 10.1016/j.ajhg.2010.06.009
Liu, T., Ren, D., Zhu, X., Yin, Z., Jin, G., Zhao, Z., et al. (2013). Transcriptional signaling pathways inversely regulated in Alzheimer’s disease and glioblastoma multiform. Sci. Rep. 3:3467. doi: 10.1038/srep03467
Liu, Y., Li, Z., Zhang, M., Deng, Y., Yi, Z., and Shi, T. (2013). Exploring the pathogenetic association between schizophrenia and type 2 diabetes mellitus diseases based on pathway analysis. BMC Med. Genomics 6(Suppl. 1):S17. doi: 10.1186/1755-8794-6-S1-S17
Liu, Y. H., Zeng, F., Wang, Y. R., Zhou, H. D., Giunta, B., Tan, J., et al. (2013). Immunity and Alzheimer’s disease: immunological perspectives on the development of novel therapies. Drug Discov. Today 18, 1212–1220. doi: 10.1016/j.drudis.2013.07.020
Malik, R., Chauhan, G., Traylor, M., Sargurupremraj, M., Okada, Y., Mishra, A., et al. (2018). Multiancestry genome-wide association study of 520,000 subjects identifies 32 loci associated with stroke and stroke subtypes. Nat. Genet. 50, 524–537. doi: 10.1038/s41588-018-0058-3
Malik, R., Traylor, M., Pulit, S. L., Bevan, S., Hopewell, J. C., Holliday, E. G., et al. (2016). Low-frequency and common genetic variation in ischemic stroke: the METASTROKE collaboration. Neurology 86, 1217–1226. doi: 10.1212/WNL.0000000000002528
Mastroeni, D., Nolz, J., Sekar, S., Delvaux, E., Serrano, G., Cuyugan, L., et al. (2018). Laser-captured microglia in the Alzheimer’s and Parkinson’s brain reveal unique regional expression profiles and suggest a potential role for hepatitis B in the Alzheimer’s brain. Neurobiol. Aging 63, 12–21. doi: 10.1016/j.neurobiolaging.2017.10.019
Mi, H., Poudel, S., Muruganujan, A., Casagrande, J. T., and Thomas, P. D. (2016). PANTHER version 10: expanded protein families and functions, and analysis tools. Nucleic Acids Res. 44, D336–D342. doi: 10.1093/nar/gkv1194
Musicco, M., Adorni, F., Di Santo, S., Prinelli, F., Pettenati, C., Caltagirone, C., et al. (2013). Inverse occurrence of cancer and Alzheimer disease: a population-based incidence study. Neurology 81, 322–328. doi: 10.1212/WNL.0b013e31829c5ec1
Network, N. S. G., International Stroke Genetics Consortium, Pulit, S. L., Mcardle, P. F., Wong, Q., Malik, R., et al. (2016). The NINDS stroke genetics network: a genome-wide association study of ischemic stroke and its subtypes. Lancet Neurol. 15, 174–184. doi: 10.1016/S1474-4422(15)00338-5
Ou, S. M., Lee, Y. J., Hu, Y. W., Liu, C. J., Chen, T. J., Fuh, J. L., et al. (2013). Does Alzheimer’s disease protect against cancers? A nationwide population-based study. Neuroepidemiology 40, 42–49. doi: 10.1159/000341411
Pina, S., Carneiro, A., Rodrigues, T., Samoes, R., Taipa, R., Melo-Pires, M., et al. (2014). Acute ischemic stroke secondary to glioblastoma. A case report. Neuroradiol. J. 27, 85–90. doi: 10.15274/NRJ-2014-10009
Plimpton, R. L., Cuellar, J., Lai, C. W., Aoba, T., Makaju, A., Franklin, S., et al. (2015). Structures of the Gbeta-CCT and PhLP1-Gbeta-CCT complexes reveal a mechanism for G-protein beta-subunit folding and Gbetagamma dimer assembly. Proc. Natl. Acad. Sci. U.S.A. 112, 2413–2418. doi: 10.1073/pnas.1419595112
Pluta, R., Bogucka-Kocka, A., Ułamek-Kozioł, M., Bogucki, J., Januszewski, S., Kocki, J., et al. (2018). Ischemic tau protein gene induction as an additional key factor driving development of Alzheimer’s phenotype changes in CA1 area of hippocampus in an ischemic model of Alzheimer’s disease. Pharmacol. Rep. 70, 881–884. doi: 10.1016/j.pharep.2018.03.004
Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M. A., Bender, D., et al. (2007). PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81, 559–575. doi: 10.1086/519795
Salminen, A., Kauppinen, A., and Kaarniranta, K. (2017). Hypoxia/ischemia activate processing of Amyloid Precursor Protein: impact of vascular dysfunction in the pathogenesis of Alzheimer’s disease. J. Neurochem. 140, 536–549. doi: 10.1111/jnc.13932
Shen, N., Chen, B., Jiang, Y., Feng, R., Liao, M., Zhang, L., et al. (2015). An updated analysis with 85,939 samples confirms the association between CR1 rs6656401 polymorphism and Alzheimer’s disease. Mol. Neurobiol. 51, 1017–1023. doi: 10.1007/s12035-014-8761-2
Shruster, A., Ben-Zur, T., Melamed, E., and Offen, D. (2012). Wnt signaling enhances neurogenesis and improves neurological function after focal ischemic injury. PLoS One 7:e40843. doi: 10.1371/journal.pone.0040843
Sims, R., Van Der Lee, S. J., Naj, A. C., Bellenguez, C., Badarinarayanm, N., Jakobsdottir, J., et al. (2017). Rare coding variants in PLCG2, ABI3, and TREM2 implicate microglial-mediated innate immunity in Alzheimer’s disease. Nat. Genet. 49, 1373–1384. doi: 10.1038/ng.3916
Stegner, D., Thielmann, I., Kraft, P., Frohman, M. A., Stoll, G., and Nieswandt, B. (2013). Pharmacological inhibition of phospholipase D protects mice from occlusive thrombus formation and ischemic stroke–brief report. Arterioscler. Thromb. Vasc. Biol. 33, 2212–2217. doi: 10.1161/ATVBAHA.113.302030
Sung, J., Song, Y. M., Choi, Y. H., Ebrahim, S., and Davey Smith, G. (2007). Hepatitis B virus seropositivity and the risk of stroke and myocardial infarction. Stroke 38, 1436–1441. doi: 10.1161/STROKEAHA.106.466268
Teng, L., Zhao, J., Wang, F., Ma, L., and Pei, G. (2010). A GPCR/secretase complex regulates beta- and gamma-secretase specificity for Abeta production and contributes to AD pathogenesis. Cell Res. 20, 138–153. doi: 10.1038/cr.2010.3
Tolppanen, A. M., Lavikainen, P., Solomon, A., Kivipelto, M., Soininen, H., and Hartikainen, S. (2013). Incidence of stroke in people with Alzheimer disease: a national register-based approach. Neurology 80, 353–358. doi: 10.1212/WNL.0b013e31827f08c5
Traylor, M., Adib-Samii, P., Harold, D., Alzheimer’s Disease Neuroimaging Initiative, International Stroke Genetics Consortium (ISGC), and Uk Young Lacunar Source Dna resource (2016). Shared genetic contribution to Ischaemic Stroke and Alzheimer’s Disease. Ann. Neurol. doi: 10.1002/ana.24621 [Epub ahead of print].
Traylor, M., Farrall, M., Holliday, E. G., Sudlow, C., Hopewell, J. C., Cheng, Y. C., et al. (2012). Genetic risk factors for ischaemic stroke and its subtypes (the METASTROKE collaboration): a meta-analysis of genome-wide association studies. Lancet Neurol. 11, 951–962. doi: 10.1016/S1474-4422(12)70234-X
Tuo, Q. Z., Lei, P., Jackman, K. A., Li, X. L., Xiong, H., Li, X. L., et al. (2017). Tau-mediated iron export prevents ferroptotic damage after ischemic stroke. Mol. Psychiatry 22, 1520–1530. doi: 10.1038/mp.2017.171
Wang, J., Vasaikar, S., Shi, Z., Greer, M., and Zhang, B. (2017). WebGestalt 2017: a more comprehensive, powerful, flexible and interactive gene set enrichment analysis toolkit. Nucleic Acids Res. 45, W130–W137. doi: 10.1093/nar/gkx356
Williams, F. M., Carter, A. M., Hysi, P. G., Surdulescu, G., Hodgkiss, D., Soranzo, N., et al. (2013). Ischemic stroke is associated with the ABO locus: the EuroCLOT study. Ann. Neurol. 73, 16–31. doi: 10.1002/ana.23838
Wu, C., Chen, J., Chen, C., Wang, W., Wen, L., Gao, K., et al. (2015). Wnt/β-catenin coupled with HIF-1α/VEGF signaling pathways involved in galangin neurovascular unit protection from focal cerebral ischemia. Sci. Rep. 5:16151. doi: 10.1038/srep16151
Xiang, Z., Xu, M., Liao, M., Jiang, Y., Jiang, Q., Feng, R., et al. (2015). Integrating genome-wide association study and brain expression data highlights cell adhesion molecules and purine metabolism in Alzheimer’s disease. Mol. Neurobiol. 52, 514–521. doi: 10.1007/s12035-014-8884-5
Zhang, B., Gaiteri, C., Bodea, L.-G., Wang, Z., Mcelwee, J., Podtelezhnikov, A. A., et al. (2013). Integrated systems approach identifies genetic nodes and networks in late-onset Alzheimer’s disease. Cell 153, 707–720. doi: 10.1016/j.cell.2013.03.030
Zhang, S., Li, X., Ma, G., Jiang, Y., Liao, M., Feng, R., et al. (2016). CLU rs9331888 polymorphism contributes to Alzheimer’s disease susceptibility in caucasian but not east Asian populations. Mol. Neurobiol. 53, 1446–1451. doi: 10.1007/s12035-015-9098-1
Keywords: Alzheimer’s disease, ischemic stroke, genome-wide association studies, gene-based analysis, pathway-based analysis
Citation: Cui P, Ma X, Li H, Lang W and Hao J (2018) Shared Biological Pathways Between Alzheimer’s Disease and Ischemic Stroke. Front. Neurosci. 12:605. doi: 10.3389/fnins.2018.00605
Received: 22 June 2018; Accepted: 10 August 2018;
Published: 07 September 2018.
Edited by:Yan Huang, Harvard Medical School, United States
Reviewed by:Dapeng Hao, Baylor College of Medicine, United States
Jianzhong Su, Wenzhou Medical University, China
Copyright © 2018 Cui, Ma, Li, Lang and Hao. 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: Junwei Hao, firstname.lastname@example.org