Next Article in Journal
Use of Saliva for Diagnosis and Monitoring the SARS-CoV-2: A General Perspective
Previous Article in Journal
Combined Use of Febuxostat and Colchicine Does Not Increase Acute Hepatotoxicity in Patients with Gout: A Retrospective Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Summary-Based Methylome-Wide Association Analyses Suggest Potential Genetically Driven Epigenetic Heterogeneity of Alzheimer’s Disease

by
Alireza Nazarian
*,
Anatoliy I. Yashin
and
Alexander M. Kulminski
*
Biodemography of Aging Research Unit, Social Science Research Institute, Duke University, Durham, NC 27705, USA
*
Authors to whom correspondence should be addressed.
J. Clin. Med. 2020, 9(5), 1489; https://doi.org/10.3390/jcm9051489
Submission received: 14 February 2020 / Revised: 30 April 2020 / Accepted: 13 May 2020 / Published: 15 May 2020
(This article belongs to the Section Mental Health)

Abstract

:
Alzheimer’s disease (AD) is a progressive neurodegenerative disorder with no curative treatment available. Exploring the genetic and non-genetic contributors to AD pathogenesis is essential to better understand its underlying biological mechanisms, and to develop novel preventive and therapeutic strategies. We investigated potential genetically driven epigenetic heterogeneity of AD through summary data-based Mendelian randomization (SMR), which combined results from our previous genome-wide association analyses with those from two publicly available methylation quantitative trait loci studies of blood and brain tissue samples. We found that 152 probes corresponding to 113 genes were epigenetically associated with AD at a Bonferroni-adjusted significance level of 5.49E-07. Of these, 10 genes had significant probes in both brain-specific and blood-based analyses. Comparing males vs. females and hypertensive vs. non-hypertensive subjects, we found that 22 and 79 probes had group-specific associations with AD, respectively, suggesting a potential role for such epigenetic modifications in the heterogeneous nature of AD. Our analyses provided stronger evidence for possible roles of four genes (i.e., AIM2, C16orf80, DGUOK, and ST14) in AD pathogenesis as they were also transcriptionally associated with AD. The identified associations suggest a list of prioritized genes for follow-up functional studies and advance our understanding of AD pathogenesis.

1. Introduction

Alzheimer’s disease (AD) is the major cause of dementia and is projected to affect more than 13 million people in the United States by 2050, thus imposing huge health and economic burdens [1,2]. Late onset AD is believed to be a multifactorial disease caused by complex interactions between various genetic and non-genetic factors [3]. Many genetic variants mapped to several chromosomal regions and genes have thus far been associated with AD by genome-wide association studies (GWAS) [4,5]; although, the vast majority of AD cases cannot be etiologically attributed to these variants [2,6]. Also, none of non-genetic AD-associated factors (e.g., age, cardiovascular risk factors, head trauma, depression, and educational attainment) has been proven to have a strong causal relationship with AD [7,8].
Epigenetic modifications of gene expression in interaction with non-genetic factors are hypothesized to contribute to AD development [6,9], particularly in light of the heterogeneous clinical manifestations of AD observed among patients with similar or identical genetic backgrounds [10]. The potential role of epigenetic mechanisms in AD pathogenesis has been widely investigated in cell lines, mouse models, post-mortem brain tissue, and blood cells [6,10,11,12,13]. Several studies have explored the global DNA methylation in AD cases compared with controls, although their findings have been inconclusive, with some reporting global hypomethylation in AD, some suggesting global hypermethylation in AD, and the others reporting no significant differences between cases and controls [12]. Previous studies have also provided many lines of evidence of associations between AD and gene-specific epigenetic modifications. They mainly investigated the DNA methylation and histone modification differences between AD cases and unaffected controls using candidate gene or genome-wide analysis approaches (e.g., pyrosequencing and array hybridization) which revealed AD-associated epigenetic modifications in some well-known AD genes, such as amyloid-β precursor protein (APP), Microtubule Associated Protein Tau (MAPT) [14], and Apolipoprotein E (APOE) [15], as well as in other genes [12]. For instance, Iwata et al. discovered CpG hypermethylation in APP and MAPT in post-mortem brain samples from AD patients, which were suggested to contribute to neural dysfunction and AD development [14]. Foraker et al. found that AD patients had a lower mean methylation level in 76 CpG sites across APOE gene compared with age-matched controls when hippocampus and frontal lobe samples were analyzed. However, APOE methylation was not statistically different between cases and controls in samples obtained from their cerebellum [15].
In most cases, epigenetically dysregulated genes were uniquely found in a single study [6,10,12,13], although AD-associated epigenetic modifications of some genes have been replicated in independent studies. For instance, several studies have reported CpG hypermethylation in the ANK1 gene in different brain regions, such as entorhinal and prefrontal cortices, superior temporal gyrus, and/or hippocampus in AD patients [16,17,18]. Hypermethylated regions overlapping DUSP22 gene were previously detected in entorhinal and dorsolateral prefrontal cortices and/or hippocampus of AD affected individuals [18,19], and CpG hypermethylation of SORBS3 was detected in the cerebral cortex of AD patients and transgenic AD mouse models [11,20]. Moreover, differentially methylated regions overlapping CDH23, RHBDF2, and RPL13 genes were reported in previous studies [16,17,21]. The mRNA expressions of these genes were also found to be altered in AD patients [16]. In addition, several genes whose associations with AD were replicated by independent GWAS [2], such as ABCA7, BIN1, CLU, HLA-DRB5, SLC24A4, and SORL1, are epigenetically implicated in AD as well [16,22,23]. The case-control studies and cell/animal models may not, however, reflect genetic contributions to AD-associated epigenetic modifications as they are more likely to identify the environmentally induced epigenetic alterations [6,9]. In addition to the studies using individual-level data, several epigenetically AD-associated genes, such as BIN1, APOC1, HLA-DRB1, HLA-DRB5, and TOMM40, have been reported by summary data-based analyses [24,25] which reflect genetically driven (i.e., through cis acting variants) epigenetic alterations [26].
In this study, we performed methylome-wide association (MWA) analyses of AD using the summary data-based Mendelian randomization (SMR) method [26] to investigate genetically driven epigenetic contributors to AD pathogenesis. Instead of analyzing individual-level data, the SMR method integrates the summary results from previous GWAS [27,28] and methylation quantitative trait loci (mQTLs) studies using blood samples [29] and brain tissue [30] in order to identify associations between AD and methylation alterations that may mediate the genetic associations examined by GWAS. Central to our study was to investigate potential genetically driven epigenetic heterogeneity of AD. Therefore, summary results from our previous GWAS which aimed to analyze genetic heterogeneity of AD in contrasting groups of subjects stratified based on their sex and history of hypertension (HTN) were used for our MWA analyses. Sex has been identified as a risk factor for AD and there are many reports highlighting sex disparities in epidemiological and clinical features of AD [31,32,33,34,35,36,37]. HTN is also a major cardiovascular risk factor for AD that may be involved in initiation and progression of the disease by causing structural and functional damages to cerebral microvasculature and promoting amyloid plaques formation [8,38,39]. By detecting several group-specific AD-associated single-nucleotide polymorphisms (SNPs) at P < 5E-06, our GWAS suggested that differences in the genetic architecture of AD between these contrasting groups may differentially contribute to AD pathogenesis [27,28]. Thus, the current study using summary results from these two GWAS may provide novel insights into potential genetically driven epigenetic heterogeneity of AD. To further validate significant findings, we compared our MWA results with those from our previous transcriptome-wide association (TWA) analyses of AD [27,28] that implemented the SMR method using the same GWAS summary results along with data from blood-based [40] and brain-specific [30,41] expression quantitative trait loci (eQTLs) studies.

2. Methods

2.1. GWAS Data

This study makes use of the results of our previous genome-wide association meta-analyses [27,28]. Briefly, these meta-analyses were performed using genotype and phenotype data from four independent datasets: (1) Cardiovascular Health Study (CHS) [42]; (2) Framingham Heart Study (FHS) [43,44]; (3) Late-Onset Alzheimer’s Disease Family Study (LOADFS) from the National Institute on Aging [45], available to the research community through the dbGaP repository (https://www.ncbi.nlm.nih.gov/gap); and (4) Health and Retirement Study (HRS) [46], which can be accessed through dbGaP and the University of Michigan restricted access webpage (http://hrsonline.isr.umich.edu/index.php?p=data). These meta-analyses were performed under five analysis plans in which the genetic basis of AD was investigated among: (1) all subjects in each dataset, (2) only males, (3) only females [27], (4) only subjects with a history of HTN, or (5) only subjects with no history of HTN [28]. AD patients were mainly diagnosed clinically based on neurologic findings (e.g., using National Institute of Neurological and Communicative Disorders and Stroke-Alzheimer’s Disease and Related Disorders Association (NINCDS-ADRDA) criteria [47]) and were either identified directly (LOADFS and FHS datasets) or reported indirectly (CHS and HRS datasets) through the International Classification of Disease codes, Ninth revision (i.e., ICD-9:331.0 code). The numbers of AD cases were 2741, 952, 1789, 1262, and 796 under plans 1–5, respectively; and the numbers of unaffected controls were 14739, 6337, 8402, 9608, and 4010, respectively. The studied subjects were all of Caucasian ancestry to make samples more homogeneous.
For each analysis plan, the additive genetic associations of ~2 million SNPs with AD were investigated by fitting logistic regression (CHS and HRS cohorts with population-based design) [48] or generalized mixed logistic regression (LOADFS and FHS cohorts with family-based design) [49] models. The top five principal components of genotype data, birth year, and sex (except plans 2 and 3) of subjects were considered as fixed-effects covariates. In the case of LOADFS and FHS cohorts, family identifier was also included as a random-effects covariate in the fitted models to adjust for potential confounding from family structure. Individual GWAS results from the four datasets were then combined by inverse-variance meta-analysis [50]. Under plans 2–5 that aimed to investigate the genetic heterogeneity of AD through stratified analyses of datasets under consideration, group-specific SNPs effects were identified by a Wald chi-square test (df = 1) [51] which was performed for any SNPs with significant association signals in only one of the contrasting groups in order to determine whether the SNPs odds ratios were significantly different between males and females (plans 2 and 3) [27] and between hypertensive and non-hypertensive subjects (plans 4 and 5) [28].
χ 2 = ( b 1 b 2 ) 2 s e 1 2 + s e 2 2
where b1 (se1) and b2 (se2) are the beta coefficients (and their standard errors) of a SNP in each of the two contrasting groups.

2.2. mQTLs Data

The summary results from two previous mQTLs studies using blood samples (n = 1980) [29] and human brain tissue (n = 1160 from a meta-analysis of three independent brain-specific mQTLs data of mostly dorsolateral prefrontal cortex and fetal brain samples) [30] were also used for our analyses. The mQTLs studies provided genome-wide CpG methylation data using the Illumina Human Methylation 450 K array. The mQTLs data in the format compatible for MWA analyses can be downloaded at: https://cnsgenomics.com/software/smr/#DataResource. The annotation of probes was in accordance with the Illumina support files for Human Methylation 450K array. Probes which were located in the inter-genic regions (IGRs) (i.e., not located within any gene or within 1.5 kb of the transcription start site of any gene [52]) were annotated to their closest genes.

2.3. MWA Analysis

Under each of the five analysis plans, two sets of MWA analyses (i.e., blood-based and brain-specific) were performed by combining the results from our GWAS with publicly available summary results from the two mQTLs studies. MWA analyses were performed by the SMR package (v 0.710) [26] to identify SNPs that might be pleiotropically associated with AD and DNA methylation changes. The SMR package was run using default input arguments. Probes that had at least one significant mQTL (i.e., a SNP with PmQTL < 5E-08) that was also among the SNPs in our GWAS were included. This resulted in the inclusion of sets of up to 90,357 and 90,848 probes with significant cis-mQTLs from blood-based and brain-specific mQTLs studies under the five analysis plans.
Associations of any probes with AD were first sought through a SMR test, and significant associations were determined at a Bonferroni-adjusted significance level of 5.49E-07 (i.e., 0.05/91000) to account for multiple comparisons. Probes with significant PSMR were then selected for heterogeneity in dependent instruments (HEIDI) testing to identify associations that were likely to arise from the pleiotropic effects of a single locus on both methylation changes and AD status (i.e., probes with PHEIDI ≥ 0.05) and not from the linkage between adjacent variants that affected AD susceptibility and methylation patterns separately (i.e., probes with PHEIDI < 0.05) [26]. Here, HRS was used as the reference panel for estimating pair-wise linkage disequilibrium and SNP clumping.
To examine the consistency of probe effects in blood-based and brain-specific analyses, the bSMR of any probes were compared between these analyses using the chi-square test mentioned above in the GWAS data section. In addition, probes that were detected in either males or females and in either hypertensive or non-hypertensive groups were subject to the chi-square test to find out whether their bSMR were significantly different between the two contrasting groups (i.e., they had group-specific effects).
Finally, lists of AD-associated genes from MWA analyses were compared to those from our previous blood-based and brain-specific TWA analyses [27,28] to identify any overlaps between epigenetically and transcriptionally AD-associated genes.

2.4. Pathway Enrichment Analysis

Pathway enrichment analysis was performed to correlate nominally AD-associated genes in our MWA results with biological processes that might contribute to AD pathogenesis. Pathway-based analyses were performed by the GSA-SNP2 (i.e., gene set analysis-single nucleotide polymorphism2) package [53] using 1329 canonical pathways provided by the Broad Institute gene set enrichment analysis (GSEA) website [54] based on information from several pathway databases such as Kyoto Encyclopedia of Genes and Genomes (KEGG) [55], REACTOME pathway knowledgebase [56], Pathway Interaction Database (PID) [57], and Matrisome Project [58]. Significant AD-associated pathways were determined using plan-specific false discovery rates (FDR) [59] at which the numbers of possible false-positively detected pathways were smaller than 1 (i.e., FDR levels between 0.05 and 0.25).

2.5. Ethics Approval

This study focuses on secondary analysis of data obtained from dbGaP and the University of Michigan [42,43,44,45,46] (please see the Supporting Acknowledgment in Additional File 1) and does not involve gathering data from human subjects directly. The study was performed according to the Duke University Institutional Review Board (IRB) guidelines.

3. Results

3.1. Blood-Based MWA Analyses

We found that 8, 31, 9, 6 and 84 probes passed both SMR at a Bonferroni-adjusted level of 5.49E-07 (PSMR between 8.73E-20 and 5.26E-07) and HEIDI (PHEIDI ≥ 0.05) tests under analysis plans 1–5, respectively (Additional File 1: Table S1). These probes were mapped to 5, 21, 9, 5, and 66 genes (71 chromosomal regions, i.e., cytogenetic bands, in total), respectively. Seventeen genes had more than one significant probe (2–9 probes per gene that were 51–61,765 base pairs apart and, in most cases, had the same top mQTLs). Top mQTLs corresponding to these probes were nominally significant (6.45E-06 ≤ PGWAS) in our genome-wide meta-analyses [27,28], except for the cg06750524 probe corresponding to the APOE gene whose top mQTLs had 2.15E-83 ≤ PGWAS ≤ 8.19E-30 under the five analysis plans of interest.

3.2. Brain-Specific MWA Analyses

There were 2, 6, 4, 4, and 27 probes that passed both SMR at a Bonferroni-adjusted threshold of 5.49E-07 (PSMR between 1.52E-12 and 5.17E-07) and HEIDI (PHEIDI ≥ 0.05) tests under plans 1–5, respectively (Additional File 1: Table S2). These probes were mapped to 2, 5, 3, 4, and 24 genes (located in 26 chromosomal regions), respectively. Six genes had more than one significant probe (2–4 probes per gene that were 5–740 base pairs apart and were mostly influenced by the same genetic signal). Again, the top mQTLs corresponding to these probes were nominally significant (4.88E-05 ≤ PGWAS) in our GWAS except for the one corresponding to cg02613937 probe, which had 3.73E-63 ≤ PGWAS ≤ 9.84E-24. This probe was mapped to the TOMM40 gene, which is near the APOE gene.

3.3. Comparison of Blood-Based and Brain-Specific MWA Results

The consistency of blood-based and brain-specific results was examined by comparing the probes effect sizes and directions (i.e., the magnitudes and signs of bSMR) between the two analyses. The directions of effects were the same for ~77% of probes in both analyses and across five plans of interest. When the blood-based and brain-specific bSMR were compared using a Wald chi-square test, less than 1% of probes (i.e., 0.006–0.073% across the five study plans) had significantly different effects at the Bonferroni-adjusted significance level. Probes corresponding to the following 10 genes were significantly associated with AD in both blood-based and brain-specific analyses (Table 1 and Table 2): NANOS2 (plan 2), HLA-DQB2 (plan 3), FAM193B (plan 4), SLC6A7, BPGM, PSTK, KRTAP5-11, LECT1, ZNF598, and C16orf80 (plan 5). All but BPGM and KRTAP5-11 had common probes in the two analyses, with directions of effects being the same and not significantly different at Bonferroni-adjusted level. The top mQTLs in blood-based and brain-specific analyses were the same for probes corresponding to NANOS2, HLA-DQB2, FAM193B, SLC6A7, KRTAP5-11, and ZNF598.

3.4. Group-Specific Findings

No probes/genes outside the APOE cluster genes region (i.e., chromosome 19q13.32) were significant in both males and females (i.e., plans 2 and 3). LOC154449 (chromosome 6q27 region) was the only gene outside the APOE cluster genes region that had AD-associated probes in blood-based MWA analyses of both hypertensive and non-hypertensive subjects (i.e., plans 4 and 5).
When the bSMR of probes were compared using a Wald chi-square test, we found that 16 of 38 blood-based probes and six of eight brain-specific probes that were detected either in males or females had sex-specific effects at Bonferroni-adjusted significance levels of 0.00132 and 0.00625, respectively (Additional File 1: Tables S3 and S4). Among 88 and 29 blood-based and brain-specific probes that were detected in either hypertensive or non-hypertensive subjects, 58 and 21 probes had significantly different effects in the two groups at Bonferroni-adjusted significance levels of 0.00057 and 0.00172, respectively (Additional File 1: Tables S5 and S6).

3.5. Comparison of MWA and GWAS Results

To investigate the novelty of our findings with respect to their potential implication in AD pathogenesis, we determined whether there were AD-associated SNPs with significant PGWAS at genome-wide (PGWAS < 5E-08) or suggestive (5E-08 ≤ PGWAS < 5E-06) significance levels within ±1 Mb regions and/or chromosomal regions of the detected probes in our genome-wide meta-analyses or in other studies reported by GRASP [4] and NHGRI-EBI GWAS [5] catalogs.
We identified AD-associated SNPs with PGWAS < 5E-08 within ±1 Mb of probes corresponding to APOE, TOMM40, and NANOS2 genes (all within the chromosome 19q13.32 region) in our genome-wide meta-analyses and previous GWAS [4,5]. No SNPs with PGWAS < 5E-08 were found within ±1 Mb flanking regions of any other probes in our meta-analyses. However, AD-associated SNPs with PGWAS < 5E-08 were previously reported by other studies within ±1 Mb of several other probes [4,5]. These probes were mapped to 22 genes (all outside the chromosome 19q13.32 region): CLIC1, BRD2, HLA-DPB1, ITIH2, PHLDA1 (plan 2), HLA-DQA2, HLA-DQB2, LECT1 (plan 3), and SLC25A2, PPT2-EGFL8, EGFL8, COL11A2, TREM1, NDUFA4, ZNF394, CHRNA2, ITIH2, LECT1, CMIP, NGFR, LOC100288866, MUM1, SIGLEC12, and EBF4 (plan 5).
In addition, the ±1 Mb flanking regions of several other probes attained 5E-08 ≤ PGWAS < 5E-06 in our or previous GWAS. Detailed information about these probes/genes can be found in Additional File 1: Tables S1 and S2. For instance, there were AD-associated SNPs at suggestive significance levels within ±1 Mb of probes corresponding to AP2A2, ADCY8, HLA-DQA2, HLA-DQB2, and SLC35C1 (all outside the chromosome 19q13.32 region) in our GWA meta-analyses.

3.6. Comparison of MWA and TWA Results

Analysis of overlaps between MWA and our previous TWA results [27,28] revealed that, among the potential epigenetically AD-associated genes, four genes also had significant AD-associated probes in TWA analyses (Table 3). These four genes, AIM2, DGUOK, ST14, and C16orf80, had significant probes in subjects with no history of HTN (i.e., plan 5). Of these genes, C16orf80 had significant probes in both blood-based and brain-specific MWA analyses; DGUOK and ST14 had AD-associated probes in blood-based analyses; and AIM2 had significant probes in brain-specific analyses. With respect to the TWA analyses, C16orf80 had significant probes in brain-specific analyses; AIM2, and DGUOK had significant probes in blood-based analyses; and ST14 had AD-associated probes in TWA analyses of both blood samples and brain tissue.

3.7. Pathway Enrichment Analyses

Pathway-based analyses (Table 4 and Table 5) revealed that AD-associated probes/genes from blood-based MWA analyses were enriched in 16 pathways (i.e., 7, 4, 4, and 3 pathways under plans 1, 2, 3, and 5, respectively). Of these, two pathways (i.e., GABA-B receptor activation (plans 1 and 3) and GABA receptor activation (plans 2 and 3)) were significant in more than one plan. We also found that nine pathways (i.e., 1, 2, 2, 3, and 3 significant pathways in plans 1–5, respectively) were associated with AD when brain-specific MWA results were enriched. Of these, two pathways (MHC class II antigen presentation (plans 1 and 3) and type II diabetes mellitus (plans 2 and 4)) were significant in more than one plan and were also enriched in both brain-specific and blood-based analyses.

4. Discussion

Despite the detection of many genetic variants and identification of several non-genetic factors that may play roles in AD susceptibility, the definitive underlying mechanisms in most AD cases is unclear. Thus, epigenetic mechanisms may be key contributors to the heterogeneous nature of AD [9,10,13,23]. The epigenetic architecture of AD has been widely investigated in case-control studies and cell/animal models [12]. The AD-associated epigenetic modifications found in these studies can be environmentally induced or genetically driven (i.e., through cis acting variants).
We combined the results from our previous GWAS [27,28] with data from two publicly available mQTLs studies of brain tissue [30] and blood samples [29] to identify genes that might be epigenetically associated with AD. In contrast to studies using individual-level data, epigenetic associations detected by summary data-based analyses are all genetically driven [26]. A major focus of our study was to explore potential genetically driven epigenetic heterogeneity of AD based on its two main risk factors (i.e., sex [31,32,33,34,35,36,37] and HTN [8,38,39]). Therefore, in order to investigate sex-specific and HTN-specific epigenetic changes, our MWA analyses were performed under five alternative plans in which summary results from GWAS on either all subjects, only males, only females [27], only subjects with a history of HTN, or only subjects with no history of HTN [28] were included in analyses.
Our analyses demonstrated that 152 probes corresponding to 113 genes were epigenetically associated with AD. The top mQTLs corresponding to these probes were mostly nominally significant in our genome-wide meta-analyses. This might be in part due to suboptimal statistical power of our analyses which can be improved by analyzing larger datasets or more importantly due to the genetic heterogeneity of AD within and between the analyzed cohorts (i.e., LOADFS, CHS, FHS, and HRS). The ±1 Mb flanking regions of ~18% and ~34% of detected probes had attained PGWAS < 5E-08 and 5E-08 ≤ PGWAS < 5E-06, respectively, in our genome-wide meta-analyses or other studies reported by GWAS databases [4,5]. Comparing our findings with those detected in other SMR-based analyses of AD [24,25] revealed that TOMM40, which had significant probes in brain-specific analyses under all five plans of our study, was epigenetically associated with AD in a previous study [24].
Investigating group-specific epigenetic alterations, we found that probes corresponding to APOE and TOMM40 genes (i.e., inside the chromosome 19q13.32 region) were significant in blood-based and brain-specific analyses, respectively, of both males and females (i.e., plans 2 and 3) and both hypertensive and non-hypertensive groups (i.e., plans 4 and 5). However, several probes (all outside the chromosome 19q13.32 region, except cg05206559 corresponding to NANOS2 gene in males) were group-specifically associated with AD, indicating potential genetically driven epigenetic heterogeneity of AD based on the two studied risk factors. For instance, we found that among 38 and eight probes that were detected in blood-based and brain-specific analyses, respectively, in either males or females, 22 probes had sex-specific effects when their bSMR were compared between the two sexes using a Wald chi-square test (Additional File 1: Tables S3 and S4). Comparing results from hypertensive and non-hypertensive groups, we found that there were 88 (blood-based analyses) and 29 (brain-specific analyses) significant probes outside the APOE region which were not in common between these two groups. Of these, 79 probes had group-specific effects when their bSMR were compared between hypertensive and non-hypertensive groups (Additional File 1: Tables S5 and S6). Addressing genetic and epigenetic heterogeneities of AD is essential for understanding its pathogenesis and developing more efficient and personalized medical interventions tailored to the genetic and epigenetic profiles of individuals.
Our MWA analyses were performed using both brain-specific and blood-based mQTLs data which provided the opportunity to assess the consistency of potential AD-associated epigenetic changes detected in these analyses. Although the pattern of DNA methylation can be tissue- or cell-specific [6,60], previous studies have demonstrated the utility of blood samples for investigating AD-associated epigenetic modifications by reporting global or gene-specific methylation changes in AD subjects compared with matched healthy controls [61,62,63,64,65]. This might be due to the systemic sequelae of AD, as AD may extensively impact cellular and molecular processes in peripheral tissues and nonneural cells including red blood cells, leukocytes, and platelets [66,67,68,69,70,71]. In addition, blood-based analyses may provide more statistical power than brain-specific studies, which generally have smaller sample sizes due to difficulties in obtaining brain samples from living subjects. Consistent with previous reports, our findings supported the feasibility of using data from blood samples to investigate epigenetic changes involved in AD. The direction of blood-based and brain-specific effects were the same for ~77% of probes and the effects of less than 1% of probes were significantly different between the two analyses across the five analysis plans of interest. We also found that probes corresponding to 10 genes were associated with AD in both blood-based and brain-speficic MWA analyses (Table 1 and Table 2). Most of these genes were previously implicated in AD at genome-wide or suggestive significance levels by GWAS [4,5], except SLC6A7, PSTK, and KRTAP5-11. AD-associated SNPs at PGWAS < 5E-08 were found within ±1 Mb of probes mapped to NANOS2, HLA-DQB2, and LECT1 in our meta-analyses and/or previous GWAS. SNPs with 5E-08 ≤ PGWAS < 5E-06 were found within ±1 Mb flanking regions of probes corresponding to FAM193B, BPGM, ZNF598, and C16orf80. Moreover, empirical evidence links some of these genes to AD in humans and animal models (e.g., SLC6A7 [72] and BPGM [71]).
It should be stressed that the identified AD-associated genes in summary-based analyses do not prove any definitive causal relationships. Instead, they suggest a list of prioritized genes whose potential roles in AD pathogenesis need to be validated by further functional studies [26]. In a recent study, Hannon et al. detected overlapping mQTL and eQTL signals with functional implications for several complex diseases/traits, such as Crohn’s disease, ulcerative colitis, blood lipids, height, and schizophrenia by comparing their SMR-based analyses [73]. Therefore, to further pinpoint potential targets, we compared the list of epigenetically AD-associated genes identified from MWA analyses with transcriptionally AD-associated genes identified from our previous TWA analyses [27,28].
Our comparisons identified a short list of four potentially AD-associated genes that had significant probes in both MWA and TWA analyses (i.e., AIM2, DGUOK, ST14, and C16orf80 in non-hypertensive subjects with PSMR between 4.62E-07 and 1.35E-10 in MWA analyses and between 2.18E-05 and 7.78E-07 in TWA analyses [28]). Probes corresponding to all genes but AIM2 had group-specific effects when their bSMR were compared between hypertensive and non-hypertensive groups using a Wald chi-square test (Additional File 1: Tables S5 and S6). AD-associated SNPs with PGWAS < 5E-08 were not found within ±1 Mb flanking regions of these probes in our meta-analyses or other studies in GWAS databases [4,5], although several SNPs with 5E-08 ≤ PGWAS < 5E-06 were previously reported within ±1 Mb of probes corresponding to AIM2 [74] and C16orf80 [75,76]. In addition, chromosomal regions corresponding to ST14 [77] (i.e., 11q24.3 region) contained previously reported AD-associated SNPs at P < 5E-08.
A review of the literature provided additional insights, strengthening the potential roles of these four genes in AD. For instance, AIM2 encodes a protein involved in regulating cell proliferation and innate immunity [78]. SNPs mapped to this gene were previously associated with white blood cells count at PGWAS < 5E-08 [79]. AIM2, along with several other proteins, were suggested to initiate inflammasome formation in response to stimuli such as viruses, bacteria, and damaged cells. Inflammasomes mediate the release of pro-inflammatory cytokines, such as IL-1β and IL-18, that are believed to be involved in AD development [80,81,82]. IL-1β may increase in the blood, cerebrospinal fluid, and brain of AD patients and blood level of IL-18 may increase in early stages of AD. IL-1β can activate astrocytes and microglia cells and stimulate the release of APP and amyloid-β () from neurons. Also, IL-18, which is overexpressed in astrocytes, microglia, and neurons around plaques, may promote formation and mediate tau protein hyper-phosphorylation [82]. It was reported that methylene blue (MB), an inhibitor of inflammasome proteins such as AIM2, NLRP3, and NLRC4 [80], can decelerate the production of plaques and neurofibrillary tangles. Thus, MB-based medications were suggested as potential treatments for AD [83]. Moreover, Wu et al. reported that AIM2 knock-out mice exhibited behavioral changes and impaired auditory fear memory [84].
DGUOK encodes a mitochondrial enzyme involved in the purine metabolism pathway [78]. Mutations in this gene were linked to some mitochondrial disorders with Mendelian inheritance, such as mitochondrial depletion syndrome [85]. Mitochondrial dysfunction has also been reported as an important finding in neurons of AD patients [86,87]. Ansoleaga et al. showed that DGUOK was downregulated in the precuneus and entorhinal cortex of patients in AD stages III-IV and V-VI (Braak and Braak staging system [88]), respectively, compared with matched healthy controls [89]. In addition, SNPs mapped to DGUOK were associated with systemic lupus erythematosus at PGWAS < 5E-08 [90]. The risks of developing AD and vascular dementia slightly increases among patients with autoimmune disorders, such as lupus erythematosus [91].
ST14 encodes a membrane serine protease with tumor suppressor activity [78] that was not associated with AD or its risk factors at PGWAS < 5E-06 by previous GWAS [4,5]. However, Wirz et al. found that the ortholog of ST14 is overexpressed (i.e., 5.39-fold change with p < 0.008) in the frontal cortex of APPswe/PS1dE9 transgenic mice harboring mutant forms of APP and PSEN1 in response to plaque development [92]. Yin et al. reported that the mouse ortholog of ST14 was upregulated in plaque-associated microglia cells in 5XFAD transgenic mice harboring mutant forms of APP and PSEN1 genes compared with aged-matched control mice [93].
C16orf80 (also known as BUG22 and CFAP20) encodes a highly conserved protein involved in the post-translational modification of Tubulin subunits of microtubules. Such modifications might be essential for microtubule function and stability in ciliated cells, such as sperm, and in neurons [94]. Microtubules are major component of neuronal transport machinery, in which defects can lead to neurodegenerative diseases (e.g., the role of microtubule-associated proteins, such as tau protein, in AD) [95,96]. In a previous study, Mendes Maia et al. reported that Drosophila melanogaster carrying mutant copies of the ortholog of C16orf80 had a short lifespan and defects in body morphology, climbing activity, and locomotion, which were mostly reversed when gene expression was restored in the nervous system [94]. However, C16orf80 was not previously associated with AD or its risk factors at PGWAS < 5E-06 [4,5].
Our pathway enrichment analyses of the brain-specific and blood-based MWA results revealed that nine and 16 pathways were associated with AD, respectively. These pathways were mostly involved in biological processes such as immune system responses (e.g., MHC class II antigen presentation), mitochondrial function (e.g., TCA cycle and respiratory electron transport), neurogenesis, synaptic function, and neurotransmitter signaling (e.g., L1CAM interactions, GABA receptor activation, neurotransmitter receptors and postsynaptic signal transmission, and transmission across chemical synapses pathways) that have been implicated in AD pathogenesis [87,97,98,99,100,101,102,103]. Two enriched pathways (i.e., MHC class II antigen presentation and type II diabetes mellitus) were common between the brain-speficic and blood-based MWA analyses, highlighting potential links between AD and immune system responses [102,103] and type II diabetes mellitus as an important vascular risk factor for AD [104].
Despite its rigor, we acknowledge that this study has limitations that could be addressed by future research using different methodologies and data. Using summary results from GWAS with larger sample sizes is likely to increase the statistical power of analyses. However, it should be noted that increasing sample sizes may not necessarily result in considerably increased power of GWAS due to the genetic heterogeneity underlying complex diseases. As mentioned above, the summary-based methylome-/transcriptome-wide approaches cannot draw definitive causal relationships between the disease of interest and detected genes [26]. Such analyses can only help generate hypotheses regarding the possible involvement of a short list of genes in the pathogenesis of the studied disorder, which need to be validated empirically. Analyzing individual-level data which provide gene expressions and epigenetic profiles for the same case and control subjects would help obtaining a more definitive view of the underlying biological processes of AD and, in addition, may allow investigating the roles of non-genetic factors (e.g., smoking, medications that interfere with DNA methylation, exposure to metals, nutritional ingredients) in the observed transcriptome and epigenome changes. This is particularly important because epigenetic alterations can be environmentally induced [6,9]. It would also be interesting to investigate whether detected epigenome changes are associated with AD progression. This requires data from different AD stages [88] with sufficient sample sizes. The CHS, FHS, HRS, and LOADFS datasets analyzed in our study do not provide disease staging information for AD subjects. Finally, investigating cell-specific (i.e., neurons and different glial cells) epigenetic alterations may provide valuable additional insights into the epigenetic architecture of AD, although small sample sizes and insufficient statistical power can be a major problem for such studies.

5. Conclusions

Our MWA analyses revealed associations between AD and probes corresponding to 113 genes. Most of these genes were not associated with AD in previous GWAS and the ±1 Mb flanking regions of ~45% of detected probes did not attain PGWAS < 5E-06 previously. The top mQTLs corresponding to these probes were mostly nominally significant in our GWAS which might be due to suboptimal sample sizes and statistical power of our analyses and/or the genetic heterogeneity of AD within and between the analyzed cohorts. Performing MWA analyses under five plans provided the opportunity to explore potential genetically driven epigenetic heterogeneity of AD in contrasting groups of subjects based on their sex and history of HTN. Comparing the MWA results from plans 2 and 3 (i.e., males vs. females) and from plans 4 and 5 (i.e., hypertensive vs. non-hypertensive subjects), we found that 22 and 79 probes were group-specifically associated with AD, respectively. Thus, this study suggests a role for genetically driven epigenetic modifications as contributing factors to the heterogeneous nature of AD, addressing of which may have translational impacts for implementing more efficient and personalized medical interventions (e.g., developing sex-specific therapeutic targets). The potential AD-genes associations detected here do not imply casualty and should only be used as a short list to prioritize candidate genes for future studies. The comparison of MWA and TWA results together with additional information from empirical studies strengthened the possible roles of four genes (i.e., AIM2, C16orf80, DGUOK, and ST14) in AD pathogenesis and helped further prioritize the list of potentially AD-associated genes for follow-up studies. Consistent with previous reports, our findings demonstrated the applicability of blood-based mQTLs data for the study of epigenetics mechanisms of AD as several genes and pathways were associated with AD in both brain-specific and blood-based MWA analyses and the probe effects detected in these analyses did not show significant differences.

Supplementary Materials

The following are available online at https://www.mdpi.com/2077-0383/9/5/1489/s1, Additional File 1 containing Supporting Acknowledgment, Table S1: Blood-based methylome-wide association results; Table S2: Brain-specific methylome-wide association results; Table S3: Wald chi-square test to compare probes effects between males and females for probes that were significant in blood-based analyses of only one of the two groups; Table S4: Wald chi-square test to compare probes effects between males and females for probes that were significant in brain-specific analyses of only one of the two groups; Table S5: Wald chi-square test to compare probes effects between hypertensive and non-hypertensive subjects for probes that were significant in blood-based analyses of only one of the two groups; Table S6: Wald chi-square test to compare probes effects between hypertensive and non-hypertensive subjects for probes that were significant in brain-specific analyses of only one of the two groups.

Author Contributions

The authors’ responsibilities were as follows: A.N. and A.M.K. designed the study, A.N. analyzed data, A.M.K. and A.I.Y. provided critical feedback, A.N., A.M.K. and A.I.Y. wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by grants from the National Institute on Aging (P01AG043352, R01AG047310, and R01AG065477). The funders had no role in study design, data collection and analysis, decision to publish, or manuscript preparation. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

Acknowledgments

Please see the Supporting Acknowledgment in Additional File 1 regarding the CHS, FHS, HRS, and LOADFS datasets used for genome-wide association meta-analyses.

Conflicts of Interest

The authors declare no competing interests.

Abbreviations

ABCA7ATP Binding Cassette Subfamily A Member 7
AD Alzheimer’s Disease
ADCY8Adenylate Cyclase 8
AIM2Absent in Melanoma 2
ANK1Ankyrin 1
AP2A2Adaptor Related Protein Complex 2 Subunit Alpha 2
APOC1Apolipoprotein C1
APOEApolipoprotein E
APPAmyloid Beta Precursor Protein
Amyloid-β
BIN1Bridging Integrator 1
BPGMBisphosphoglycerate Mutase
BRD2Bromodomain Containing 2
BUG22Basal Body Upregulated Gene 22
C10orf54Chromosome 10 Open Reading Frame 54
C16orf80Chromosome 16 Open Reading Frame 80
CDH23Cadherin Related 23
CFAP20Cilia and Flagella Associated Protein 20
CHRNA2Cholinergic Receptor Nicotinic Alpha 2 Subunit
CHS Cardiovascular Health Study
CLIC1Chloride Intracellular Channel 1
CLUClusterin
CMIPC-Maf Inducing Protein
COL11A2Collagen Type XI Alpha 2 Chain
dbGaP The Database of Genotypes and Phenotypes
DGUOKDeoxyguanosine Kinase
DUSP22Dual Specificity Phosphatase 22
EBF4Early B Cell Factor Family Member 4
EGFL8Epidermal Growth Factor-Like Like Domain Multiple 8
eQTL Expression Quantitative trait Locus
FAM193BFamily with Sequence Similarity 193 Member B
FDR False Discovery Rate
FHS Framingham Heart Study
GABAGamma-Aminobutyric Acid
GRASP Genome-Wide Repository of Associations Between SNPs and Phenotypes
GSA Gene Set Analysis
GSA-SNP2 Gene Set Analysis-Single-Nucleotide-Polymorphism-2
GSEA Gene Set Enrichment Analysis
GWAS Genome-Wide Association Study
HEIDI Heterogeneity in Dependent Instruments
HLA-DPB1Human Leukocyte Antigen Class II, DP Beta 1
HLA-DQA2Human Leukocyte Antigen Class II, DQ Alpha 2
HLA-DQB2Human Leukocyte Antigen Class II, DQ Beta 2
HLA-DRB1Human Leukocyte Antigen Class II, DR Beta 1
HLA-DRB5Human Leukocyte Antigen Class II, DR Beta 5
HRS Health and Retirement Study
HTN Hypertension
ICD-9 International Classification of Disease codes, Ninth revision
IGR Inter-Genic Region
IL-18Interleukin 18
IL-1βInterleukin 1 Beta
IRB Institutional Review Board
ITIH2Inter-Alpha-Trypsin Inhibitor Heavy Chain 2
KEGG Kyoto Encyclopedia of Genes and Genomes
KRTAP5-11Keratin Associated Protein 5-11
L1CAML1 Cell Adhesion Molecule
LECT1Leukocyte Cell Derived Chemotaxin 1
LOADFS Late-Onset Alzheimer's Disease Family Study
LOC100288866Uncharacterized LOC100288866
LOC154449Uncharacterized LOC154449
MAPTMicrotubule Associated Protein Tau
MB Methylene Blue
MHCMajor Histocompatibility Complex
mQTL Methylation Quantitative trait Locus
MUM1Melanoma Ubiquitous Mutated Protein 1
MWA Methylome-Wide Association
NABA Matrisome Project
NANOS2Nanos C2HC-Type Zinc Finger 2
NDUFA4NDUFA4 Mitochondrial Complex Associated
NGFRNerve Growth Factor Receptor
NHGRI-EBI GWAS National Human Genome Research Institute-European Bioinformatics Institute Genome-Wide Association Studies Catalog
NINCDS-ADRDA National Institute of Neurological and Communicative Disorders and Stroke of the United States-the Alzheimer’s Disease and Related Disorders Association
NLRC4Nucleotide-Binding Oligomerization Domain, Leucine Rich Repeat and Caspase Recruitment Domain Containing 4
NLRP3Nucleotide-Binding Oligomerization Domain, Leucine Rich Repeat and Pyrin Domain Containing 3
PHLDA1Pleckstrin Homology Like Domain Family A Member 1
PID Pathway Interaction Database
PPT2-EGFL8Palmitoyl-Protein Thioesterase 2-Epidermal Growth Factor-Like Like Domain Multiple 8 Readthrough
PSEN1Presenilin 1
PSTKPhosphoseryl-TRNA Kinase
RHBDF2Rhomboid 5 Homolog 2
RPL13Ribosomal Protein L13
SIGLEC12Sialic Acid Binding Immunoglobulin Like Lectin 12
SLC24A4Solute Carrier Family 24 Member 4
SLC25A2Solute Carrier Family 25 Member 2
SLC35C1Solute Carrier Family 35 Member C1
SLC6A7Solute Carrier Family 6 Member 7
SMR Summary Data-Based Mendelian Randomization
SNP Single-Nucleotide Polymorphism
SORBS3Sorbin And SH3 Domain Containing 3
SORL1Sortilin Related Receptor 1
ST14Suppression of Tumorigenicity 14
TCA Tricarboxylic Acid
TOMM40Translocase of Outer Mitochondrial Membrane 40
TREM1Triggering Receptor Expressed on Myeloid Cells 1
TWA Transcriptome-Wide Association
ZNF394Zinc Finger Protein 394
ZNF598Zinc Finger Protein 598

References

  1. Alzheimer’s Association. 2016 Alzheimer’s disease facts and figures. Alzheimers Dement. 2016, 12, 459–509. [Google Scholar] [CrossRef] [PubMed]
  2. Ridge, P.G.; Hoyt, K.B.; Boehme, K.; Mukherjee, S.; Crane, P.K.; Haines, J.L.; Mayeux, R.; Farrer, L.A.; Pericak-Vance, M.A.; Schellenberg, G.D.; et al. Assessment of the genetic variance of late-onset Alzheimer’s disease. Neurobiol. Aging 2016, 41, 200.e13–200.e20. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Raghavan, N.; Tosto, G. Genetics of Alzheimer’s disease: The importance of polygenic and epistatic components. Curr. Neurol. Neurosci. Rep. 2017, 17, 78. [Google Scholar] [CrossRef] [PubMed]
  4. Leslie, R.; O’Donnell, C.J.; Johnson, A.D. GRASP: Analysis of genotype-phenotype results from 1390 genome-wide association studies and corresponding open access database. Bioinformatics 2014, 30, i185–i194. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. MacArthur, J.; Bowler, E.; Cerezo, M.; Gil, L.; Hall, P.; Hastings, E.; Junkins, H.; McMahon, A.; Milano, A.; Morales, J.; et al. The new NHGRI-EBI Catalog of published genome-wide association studies (GWAS Catalog). Nucleic Acids Res. 2017, 45, D896–D901. [Google Scholar] [CrossRef]
  6. Sanchez-Mut, J.V.; Gräff, J. Epigenetic alterations in Alzheimer’s disease. Front. Behav. Neurosci. 2015, 9, 347. [Google Scholar] [CrossRef]
  7. Daviglus, M.L.; Bell, C.C.; Berrettini, W.; Bowen, P.E.; Connolly, E.S.; Cox, N.J.; Dunbar-Jacob, J.M.; Granieri, E.C.; Hunt, G.; McGarry, K.; et al. NIH state-of-the-science conference statement: Preventing Alzheimer’s disease and cognitive decline. NIH Consens. State Sci. Statements 2010, 27, 1–30. [Google Scholar]
  8. Power, M.C.; Weuve, J.; Gagne, J.J.; McQueen, M.B.; Viswanathan, A.; Blacker, D. The association between blood pressure and incident Alzheimer disease: A systematic review and meta-analysis. Epidemiology 2011, 22, 646–659. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Lahiri, D.K.; Zawia, N.H.; Greig, N.H.; Sambamurti, K.; Maloney, B. Early-life events may trigger biochemical pathways for Alzheimer’s disease: The “LEARn” model. Biogerontology 2008, 9, 375–379. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Yokoyama, A.S.; Rutledge, J.C.; Medici, V. DNA methylation alterations in Alzheimer’s disease. Environ. Epigenet. 2017, 3, dvx008. [Google Scholar] [CrossRef]
  11. Sanchez-Mut, J.V.; Aso, E.; Panayotis, N.; Lott, I.; Dierssen, M.; Rabano, A.; Urdinguio, R.G.; Fernandez, A.F.; Astudillo, A.; Martin-Subero, J.I.; et al. DNA methylation map of mouse and human brain identifies target genes in Alzheimer’s disease. Brain 2013, 136, 3018–3027. [Google Scholar] [CrossRef] [PubMed]
  12. Wen, K.-X.; Miliç, J.; El-Khodor, B.; Dhana, K.; Nano, J.; Pulido, T.; Kraja, B.; Zaciragic, A.; Bramer, W.M.; Troup, J.; et al. The role of DNA methylation and histone modifications in neurodegenerative diseases: A systematic review. PLoS ONE 2016, 11, e0167201. [Google Scholar] [CrossRef]
  13. Liu, X.; Jiao, B.; Shen, L. The epigenetics of Alzheimer’s disease: Factors and therapeutic implications. Front. Genet. 2018, 9, 579. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Iwata, A.; Nagata, K.; Hatsuta, H.; Takuma, H.; Bundo, M.; Iwamoto, K.; Tamaoka, A.; Murayama, S.; Saido, T.; Tsuji, S. Altered CpG methylation in sporadic Alzheimer’s disease is associated with APP and MAPT dysregulation. Hum. Mol. Genet. 2014, 23, 648–656. [Google Scholar] [CrossRef] [PubMed]
  15. Foraker, J.; Millard, S.P.; Leong, L.; Thomson, Z.; Chen, S.; Keene, C.D.; Bekris, L.M.; Yu, C.-E. The APOE gene is differentially methylated in Alzheimer’s disease. J. Alzheimers Dis. 2015, 48, 745–755. [Google Scholar] [CrossRef] [PubMed]
  16. De Jager, P.L.; Srivastava, G.; Lunnon, K.; Burgess, J.; Schalkwyk, L.C.; Yu, L.; Eaton, M.L.; Keenan, B.T.; Ernst, J.; McCabe, C.; et al. Alzheimer’s disease: Early alterations in brain DNA methylation at ANK1, BIN1, RHBDF2 and other loci. Nat. Neurosci. 2014, 17, 1156–1163. [Google Scholar] [CrossRef] [PubMed]
  17. Lunnon, K.; Smith, R.; Hannon, E.; De Jager, P.L.; Srivastava, G.; Volta, M.; Troakes, C.; Al-Sarraj, S.; Burrage, J.; Macdonald, R.; et al. Methylomic profiling implicates cortical deregulation of ANK1 in Alzheimer’s disease. Nat. Neurosci. 2014, 17, 1164–1170. [Google Scholar] [CrossRef] [Green Version]
  18. Semick, S.A.; Bharadwaj, R.A.; Collado-Torres, L.; Tao, R.; Shin, J.H.; Deep-Soboslay, A.; Weiss, J.R.; Weinberger, D.R.; Hyde, T.M.; Kleinman, J.E.; et al. Integrated DNA methylation and gene expression profiling across multiple brain regions implicate novel genes in Alzheimer’s disease. Acta Neuropathol. 2019, 137, 557–569. [Google Scholar] [CrossRef]
  19. Sanchez-Mut, J.V.; Aso, E.; Heyn, H.; Matsuda, T.; Bock, C.; Ferrer, I.; Esteller, M. Promoter hypermethylation of the phosphatase DUSP22 mediates PKA-dependent TAU phosphorylation and CREB activation in Alzheimer’s disease. Hippocampus 2014, 24, 363–368. [Google Scholar] [CrossRef] [Green Version]
  20. Siegmund, K.D.; Connor, C.M.; Campan, M.; Long, T.I.; Weisenberger, D.J.; Biniszkiewicz, D.; Jaenisch, R.; Laird, P.W.; Akbarian, S. DNA methylation in the human cerebral cortex is dynamically regulated throughout the life span and involves differentiated neurons. PLoS ONE 2007, 2, e895. [Google Scholar] [CrossRef]
  21. Lord, J.; Cruchaga, C. The epigenetic landscape of Alzheimer’s disease. Nat. Neurosci. 2014, 17, 1138–1140. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Yu, L.; Chibnik, L.B.; Srivastava, G.P.; Pochet, N.; Yang, J.; Xu, J.; Kozubek, J.; Obholzer, N.; Leurgans, S.E.; Schneider, J.A.; et al. Association of Brain DNA methylation in SORL1, ABCA7, HLA-DRB5, SLC24A4, and BIN1 with pathological diagnosis of Alzheimer disease. JAMA Neurol. 2015, 72, 15–24. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Fetahu, I.S.; Ma, D.; Rabidou, K.; Argueta, C.; Smith, M.; Liu, H.; Wu, F.; Shi, Y.G. Epigenetic signatures of methylated DNA cytosine in Alzheimer’s disease. Sci. Adv. 2019, 5, eaaw2880. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Marioni, R.E.; Harris, S.E.; Zhang, Q.; McRae, A.F.; Hagenaars, S.P.; Hill, W.D.; Davies, G.; Ritchie, C.W.; Gale, C.R.; Starr, J.M.; et al. GWAS on family history of Alzheimer’s disease. Transl. Psychiatry 2018, 8, 99. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Zhao, T.; Hu, Y.; Zang, T.; Wang, Y. Integrate GWAS, eQTL, and mQTL data to identify Alzheimer’s disease-related genes. Front. Genet. 2019, 10, 1021. [Google Scholar] [CrossRef] [Green Version]
  26. Zhu, Z.; Zhang, F.; Hu, H.; Bakshi, A.; Robinson, M.R.; Powell, J.E.; Montgomery, G.W.; Goddard, M.E.; Wray, N.R.; Visscher, P.M.; et al. Integration of summary data from GWAS and eQTL studies predicts complex trait gene targets. Nat. Genet. 2016, 48, 481–487. [Google Scholar] [CrossRef]
  27. Nazarian, A.; Yashin, A.I.; Kulminski, A.M. Genome-wide analysis of genetic predisposition to Alzheimer’s disease and related sex disparities. Alzheimers Res. Ther. 2019, 11, 5. [Google Scholar] [CrossRef] [Green Version]
  28. Nazarian, A.; Arbeev, K.G.; Yashkin, A.P.; Kulminski, A.M. Genetic heterogeneity of Alzheimer’s disease in subjects with and without hypertension. GeroScience 2019, 41, 137–154. [Google Scholar] [CrossRef]
  29. McRae, A.F.; Marioni, R.E.; Shah, S.; Yang, J.; Powell, J.E.; Harris, S.E.; Gibson, J.; Henders, A.K.; Bowdler, L.; Painter, J.N.; et al. Identification of 55,000 Replicated DNA Methylation QTL. Sci. Rep. 2018, 8, 17605. [Google Scholar] [CrossRef] [Green Version]
  30. Qi, T.; Wu, Y.; Zeng, J.; Zhang, F.; Xue, A.; Jiang, L.; Zhu, Z.; Kemper, K.; Yengo, L.; Zheng, Z.; et al. Identifying gene targets for brain-related traits using transcriptomic and methylomic data from blood. Nat. Commun. 2018, 9, 2282. [Google Scholar] [CrossRef] [Green Version]
  31. Genin, E.; Hannequin, D.; Wallon, D.; Sleegers, K.; Hiltunen, M.; Combarros, O.; Bullido, M.J.; Engelborghs, S.; De Deyn, P.; Berr, C.; et al. APOE and Alzheimer disease: A major gene with semi-dominant inheritance. Mol. Psychiatry 2011, 16, 903–907. [Google Scholar] [CrossRef] [PubMed]
  32. Andersen, K.; Launer, L.J.; Dewey, M.E.; Letenneur, L.; Ott, A.; Copeland, J.R.; Dartigues, J.F.; Kragh-Sorensen, P.; Baldereschi, M.; Brayne, C.; et al. Gender differences in the incidence of AD and vascular dementia: The EURODEM Studies. EURODEM incidence research group. Neurology 1999, 53, 1992–1997. [Google Scholar] [CrossRef] [PubMed]
  33. Carter, C.L.; Resnick, E.M.; Mallampalli, M.; Kalbarczyk, A. Sex and gender differences in Alzheimer’s disease: Recommendations for future research. J. Womens Health (Larchmt.) 2012, 21, 1018–1023. [Google Scholar] [CrossRef] [PubMed]
  34. Mayeux, R. Epidemiology of neurodegeneration. Annu. Rev. Neurosci. 2003, 26, 81–104. [Google Scholar] [CrossRef]
  35. Mielke, M.M.; Vemuri, P.; Rocca, W.A. Clinical epidemiology of Alzheimer’s disease: Assessing sex and gender differences. Clin. Epidemiol. 2014, 6, 37–48. [Google Scholar] [CrossRef] [Green Version]
  36. Henderson, V.W.; Buckwalter, J.G. Cognitive deficits of men and women with Alzheimer’s disease. Neurology 1994, 44, 90–96. [Google Scholar] [CrossRef]
  37. Barnes, L.L.; Wilson, R.S.; Bienias, J.L.; Schneider, J.A.; Evans, D.A.; Bennett, D.A. Sex differences in the clinical manifestations of Alzheimer disease pathology. Arch. Gen. Psychiatry 2005, 62, 685–691. [Google Scholar] [CrossRef]
  38. Faraco, G.; Iadecola, C. Hypertension: A harbinger of stroke and dementia. Hypertension 2013, 62, 810–817. [Google Scholar] [CrossRef] [Green Version]
  39. Csiszar, A.; Tarantini, S.; Fülöp, G.A.; Kiss, T.; Valcarcel-Ares, M.N.; Galvan, V.; Ungvari, Z.; Yabluchanskiy, A. Hypertension impairs neurovascular coupling and promotes microvascular injury: Role in exacerbation of Alzheimer’s disease. GeroScience 2017, 39, 359–372. [Google Scholar] [CrossRef]
  40. Lloyd-Jones, L.R.; Holloway, A.; McRae, A.; Yang, J.; Small, K.; Zhao, J.; Zeng, B.; Bakshi, A.; Metspalu, A.; Dermitzakis, M.; et al. The genetic architecture of gene expression in peripheral blood. Am. J. Hum. Genet. 2017, 100, 228–237. [Google Scholar] [CrossRef] [Green Version]
  41. GTEx Consortium Genetic effects on gene expression across human tissues. Nature 2017, 550, 204–213. [CrossRef] [PubMed]
  42. Fried, L.P.; Borhani, N.O.; Enright, P.; Furberg, C.D.; Gardin, J.M.; Kronmal, R.A.; Kuller, L.H.; Manolio, T.A.; Mittelmark, M.B.; Newman, A. The cardiovascular health study: Design and rationale. Ann. Epidemiol. 1991, 1, 263–276. [Google Scholar] [CrossRef]
  43. Dawber, T.R.; Meadors, G.F.; Moore, F.E. Epidemiological approaches to heart disease: The Framingham study. Am. J. Public Health Nations Health 1951, 41, 279–286. [Google Scholar] [CrossRef] [PubMed]
  44. Feinleib, M.; Kannel, W.B.; Garrison, R.J.; McNamara, P.M.; Castelli, W.P. The Framingham offspring study: Design and preliminary data. Prev. Med. 1975, 4, 518–525. [Google Scholar] [CrossRef]
  45. Lee, J.H.; Cheng, R.; Graff-Radford, N.; Foroud, T.; Mayeux, R. Analyses of the national institute on aging late-onset Alzheimer’s disease family study: Implication of additional loci. Arch. Neurol. 2008, 65, 1518–1526. [Google Scholar] [CrossRef] [Green Version]
  46. Sonnega, A.; Faul, J.D.; Ofstedal, M.B.; Langa, K.M.; Phillips, J.W.; Weir, D.R. Cohort profile: The health and retirement study (HRS). Int. J. Epidemiol. 2014, 43, 576–585. [Google Scholar] [CrossRef] [Green Version]
  47. McKhann, G.; Drachman, D.; Folstein, M.; Katzman, R.; Price, D.; Stadlan, E.M. Clinical diagnosis of Alzheimer’s disease: Report of the NINCDS-ADRDA Work Group under the auspices of Department of Health and Human Services Task Force on Alzheimer’s Disease. Neurology 1984, 34, 939–944. [Google Scholar] [CrossRef] [Green Version]
  48. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.R.; Bender, D.; Maller, J.; Sklar, P.; de Bakker, P.I.W.; Daly, M.J.; et al. PLINK: A tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 2007, 81, 559–575. [Google Scholar] [CrossRef] [Green Version]
  49. Bates, D.; Mächler, M.; Bolker, B.; Walker, S. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 2015, 67, 1–48. [Google Scholar] [CrossRef]
  50. Mägi, R.; Morris, A.P. GWAMA: Software for genome-wide association meta-analysis. BMC Bioinform. 2010, 11, 288. [Google Scholar] [CrossRef] [Green Version]
  51. Allison, P.D. Comparing logit and probit coefficients across groups. Sociol. Methods Res. 1999, 28, 186–208. [Google Scholar] [CrossRef]
  52. Hannon, E.; Gorrie-Stone, T.J.; Smart, M.C.; Burrage, J.; Hughes, A.; Bao, Y.; Kumari, M.; Schalkwyk, L.C.; Mill, J. Leveraging DNA-methylation quantitative-trait loci to characterize the relationship between methylomic variation, gene expression, and complex traits. Am. J. Hum. Genet. 2018, 103, 654–665. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Yoon, S.; Nguyen, H.C.T.; Yoo, Y.J.; Kim, J.; Baik, B.; Kim, S.; Kim, J.; Kim, S.; Nam, D. Efficient pathway enrichment and network analysis of GWAS summary data using GSA-SNP2. Nucleic Acids Res. 2018, 46, e60. [Google Scholar] [CrossRef] [Green Version]
  54. Subramanian, A.; Tamayo, P.; Mootha, V.K.; Mukherjee, S.; Ebert, B.L.; Gillette, M.A.; Paulovich, A.; Pomeroy, S.L.; Golub, T.R.; Lander, E.S.; et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. USA 2005, 102, 15545–15550. [Google Scholar] [CrossRef] [Green Version]
  55. Kanehisa, M.; Goto, S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef] [PubMed]
  56. Fabregat, A.; Jupe, S.; Matthews, L.; Sidiropoulos, K.; Gillespie, M.; Garapati, P.; Haw, R.; Jassal, B.; Korninger, F.; May, B.; et al. The reactome pathway knowledgebase. Nucleic Acids Res. 2018, 46, D649–D655. [Google Scholar] [CrossRef] [PubMed]
  57. Schaefer, C.F.; Anthony, K.; Krupa, S.; Buchoff, J.; Day, M.; Hannay, T.; Buetow, K.H. PID: The pathway interaction database. Nucleic Acids Res. 2009, 37, D674–D679. [Google Scholar] [CrossRef]
  58. Naba, A.; Clauser, K.R.; Hoersch, S.; Liu, H.; Carr, S.A.; Hynes, R.O. The matrisome: In silico definition and in vivo characterization by proteomics of normal and tumor extracellular matrices. Mol. Cell. Proteom. 2012, 11, M111.014647. [Google Scholar] [CrossRef] [Green Version]
  59. Benjamini, Y.; Hochberg, Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B Methodol. 1995, 57, 289–300. [Google Scholar] [CrossRef]
  60. Smith, A.K.; Kilaru, V.; Kocak, M.; Almli, L.M.; Mercer, K.B.; Ressler, K.J.; Tylavsky, F.A.; Conneely, K.N. Methylation quantitative trait loci (meQTLs) are consistently detected across ancestry, developmental stage, and tissue type. BMC Genom. 2014, 15, 145. [Google Scholar] [CrossRef] [Green Version]
  61. Bollati, V.; Galimberti, D.; Pergoli, L.; Dalla Valle, E.; Barretta, F.; Cortini, F.; Scarpini, E.; Bertazzi, P.A.; Baccarelli, A. DNA methylation in repetitive elements and Alzheimer disease. Brain Behav. Immun. 2011, 25, 1078–1083. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Chang, L.; Wang, Y.; Ji, H.; Dai, D.; Xu, X.; Jiang, D.; Hong, Q.; Ye, H.; Zhang, X.; Zhou, X.; et al. Elevation of peripheral BDNF promoter methylation links to the risk of Alzheimer’s disease. PLoS ONE 2014, 9, e110773. [Google Scholar] [CrossRef] [PubMed]
  63. Di Francesco, A.; Arosio, B.; Falconi, A.; Micioni Di Bonaventura, M.V.; Karimi, M.; Mari, D.; Casati, M.; Maccarrone, M.; D’Addario, C. Global changes in DNA methylation in Alzheimer’s disease peripheral blood mononuclear cells. Brain Behav. Immun. 2015, 45, 139–144. [Google Scholar] [CrossRef] [PubMed]
  64. Nagata, T.; Kobayashi, N.; Ishii, J.; Shinagawa, S.; Nakayama, R.; Shibata, N.; Kuerban, B.; Ohnuma, T.; Kondo, K.; Arai, H.; et al. Association between DNA Methylation of the BDNF Promoter Region and Clinical Presentation in Alzheimer’s Disease. Dement. Geriatr. Cogn. Dis. Extra 2015, 5, 64–73. [Google Scholar] [CrossRef] [PubMed]
  65. Ji, H.; Wang, Y.; Liu, G.; Xu, X.; Dai, D.; Chen, Z.; Zhou, D.; Zhou, X.; Han, L.; Li, Y.; et al. OPRK1 promoter hypermethylation increases the risk of Alzheimer’s disease. Neurosci. Lett. 2015, 606, 24–29. [Google Scholar] [CrossRef]
  66. Blass, J.P.; Hanin, I.; Barclay, L.; Kopp, U.; Reding, M.J. Red blood cell abnormalities in Alzheimer disease. J. Am. Geriatr. Soc. 1985, 33, 401–405. [Google Scholar] [CrossRef]
  67. Sevush, S.; Jy, W.; Horstman, L.L.; Mao, W.W.; Kolodny, L.; Ahn, Y.S. Platelet activation in Alzheimer disease. Arch. Neurol. 1998, 55, 530–536. [Google Scholar] [CrossRef]
  68. Etcheberrigaray, R.; Ibarreta, D. Ionic channels and second messenger alterations in Alzheimer’s disease. Relevance of studies in nonneuronal cells. Rev. Neurol. 2001, 33, 740–749. [Google Scholar] [CrossRef] [PubMed]
  69. Gibson, G.E.; Huang, H.-M. Oxidative processes in the brain and non-neuronal tissues as biomarkers of Alzheimer’s disease. Front. Biosci. 2002, 7, d1007–d1015. [Google Scholar] [CrossRef] [Green Version]
  70. Catricala, S.; Torti, M.; Ricevuti, G. Alzheimer disease and platelets: How’s that relevant. Immun. Ageing 2012, 9, 20. [Google Scholar] [CrossRef] [Green Version]
  71. Kaminsky, Y.G.; Reddy, V.P.; Ashraf, G.M.; Ahmad, A.; Benberin, V.V.; Kosenko, E.A.; Aliev, G. Age-related defects in erythrocyte 2,3-diphosphoglycerate metabolism in dementia. Aging Dis. 2013, 4, 244–255. [Google Scholar] [CrossRef] [PubMed]
  72. Hokama, M.; Oka, S.; Leon, J.; Ninomiya, T.; Honda, H.; Sasaki, K.; Iwaki, T.; Ohara, T.; Sasaki, T.; LaFerla, F.M.; et al. Altered expression of diabetes-related genes in Alzheimer’s disease brains: The Hisayama study. Cereb. Cortex 2014, 24, 2476–2488. [Google Scholar] [CrossRef] [PubMed]
  73. Hannon, E.; Weedon, M.; Bray, N.; O’Donovan, M.; Mill, J. Pleiotropic effects of trait-associated genetic variation on DNA methylation: Utility for refining GWAS loci. Am. J. Hum. Genet. 2017, 100, 954–959. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Li, H.; Wetten, S.; Li, L.; St Jean, P.L.; Upmanyu, R.; Surh, L.; Hosford, D.; Barnes, M.R.; Briley, J.D.; Borrie, M.; et al. Candidate single-nucleotide polymorphisms from a genomewide association study of Alzheimer disease. Arch. Neurol. 2008, 65, 45–53. [Google Scholar] [CrossRef] [Green Version]
  75. Han, M.-R.; Schellenberg, G.D.; Wang, L.-S. Alzheimer’s Disease Neuroimaging Initiative Genome-wide association reveals genetic effects on human Aβ42 and τ protein levels in cerebrospinal fluids: A case control study. BMC Neurol. 2010, 10, 90. [Google Scholar] [CrossRef] [Green Version]
  76. Beecham, G.W.; Hamilton, K.; Naj, A.C.; Martin, E.R.; Huentelman, M.; Myers, A.J.; Corneveaux, J.J.; Hardy, J.; Vonsattel, J.-P.; Younkin, S.G.; et al. Genome-wide association meta-analysis of neuropathologic features of Alzheimer’s disease and related dementias. PLoS Genet. 2014, 10, e1004606. [Google Scholar] [CrossRef] [Green Version]
  77. Sherva, R.; Tripodis, Y.; Bennett, D.A.; Chibnik, L.B.; Crane, P.K.; de Jager, P.L.; Farrer, L.A.; Saykin, A.J.; Shulman, J.M.; Naj, A.; et al. Genome-wide association study of the rate of cognitive decline in Alzheimer’s disease. Alzheimers Dement. 2014, 10, 45–52. [Google Scholar] [CrossRef] [Green Version]
  78. Stelzer, G.; Rosen, N.; Plaschkes, I.; Zimmerman, S.; Twik, M.; Fishilevich, S.; Stein, T.I.; Nudel, R.; Lieder, I.; Mazor, Y.; et al. The GeneCards suite: From gene data mining to disease genome sequence analyses. Curr. Protoc. Bioinform. 2016, 54, 1.30.1–1.30.33. [Google Scholar] [CrossRef]
  79. Keller, M.F.; Reiner, A.P.; Okada, Y.; van Rooij, F.J.A.; Johnson, A.D.; Chen, M.-H.; Smith, A.V.; Morris, A.P.; Tanaka, T.; Ferrucci, L.; et al. Trans-ethnic meta-analysis of white blood cell phenotypes. Hum. Mol. Genet. 2014, 23, 6944–6960. [Google Scholar] [CrossRef]
  80. Ahn, H.; Kang, S.G.; Yoon, S.; Ko, H.-J.; Kim, P.-H.; Hong, E.-J.; An, B.-S.; Lee, E.; Lee, G.-S. Methylene blue inhibits NLRP3, NLRC4, AIM2, and non-canonical inflammasome activation. Sci. Rep. 2017, 7, 12409. [Google Scholar] [CrossRef]
  81. Freeman, L.C.; Ting, J.P.-Y. The pathogenic role of the inflammasome in neurodegenerative diseases. J. Neurochem. 2016, 136 (Suppl. 1), 29–38. [Google Scholar] [CrossRef] [PubMed]
  82. Liu, L.; Chan, C. The role of inflammasome in Alzheimer’s disease. Ageing Res. Rev. 2014, 15, 6–15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Oz, M.; Lorke, D.E.; Petroianu, G.A. Methylene blue and Alzheimer’s disease. Biochem. Pharmacol. 2009, 78, 927–932. [Google Scholar] [CrossRef] [PubMed]
  84. Wu, P.-J.; Liu, H.-Y.; Huang, T.-N.; Hsueh, Y.-P. AIM2 inflammasomes regulate neuronal morphology and influence anxiety and memory in mice. Sci. Rep. 2016, 6, 32405. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Tadiboyina, V.T.; Rupar, A.; Atkison, P.; Feigenbaum, A.; Kronick, J.; Wang, J.; Hegele, R.A. Novel mutation in DGUOK in hepatocerebral mitochondrial DNA depletion syndrome associated with cystathioninuria. Am. J. Med. Genet. A 2005, 135, 289–291. [Google Scholar] [CrossRef]
  86. Lunnon, K.; Keohane, A.; Pidsley, R.; Newhouse, S.; Riddoch-Contreras, J.; Thubron, E.B.; Devall, M.; Soininen, H.; Kłoszewska, I.; Mecocci, P.; et al. Mitochondrial genes are altered in blood early in Alzheimer’s disease. Neurobiol. Aging 2017, 53, 36–47. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Querfurth, H.W.; LaFerla, F.M. Alzheimer’s Disease. N. Engl. J. Med. 2010, 362, 329–344. [Google Scholar] [CrossRef] [Green Version]
  88. Braak, H.; Braak, E. Staging of Alzheimer’s disease-related neurofibrillary changes. Neurobiol. Aging 1995, 16, 271–278; discussion 278–284. [Google Scholar] [CrossRef]
  89. Ansoleaga, B.; Jové, M.; Schlüter, A.; Garcia-Esparcia, P.; Moreno, J.; Pujol, A.; Pamplona, R.; Portero-Otín, M.; Ferrer, I. Deregulation of purine metabolism in Alzheimer’s disease. Neurobiol. Aging 2015, 36, 68–80. [Google Scholar] [CrossRef]
  90. Yang, W.; Tang, H.; Zhang, Y.; Tang, X.; Zhang, J.; Sun, L.; Yang, J.; Cui, Y.; Zhang, L.; Hirankarn, N.; et al. Meta-analysis followed by replication identifies loci in or near CDKN1B, TET3, CD80, DRAM1, and ARID5B as associated with systemic lupus erythematosus in Asians. Am. J. Hum. Genet. 2013, 92, 41–51. [Google Scholar] [CrossRef] [Green Version]
  91. Wotton, C.J.; Goldacre, M.J. Associations between specific autoimmune diseases and subsequent dementia: Retrospective record-linkage cohort study, UK. J. Epidemiol. Community Health 2017, 71, 576–583. [Google Scholar] [CrossRef] [PubMed]
  92. Wirz, K.T.S.; Bossers, K.; Stargardt, A.; Kamphuis, W.; Swaab, D.F.; Hol, E.M.; Verhaagen, J. Cortical beta amyloid protein triggers an immune response, but no synaptic changes in the APPswe/PS1dE9 Alzheimer’s disease mouse model. Neurobiol. Aging 2013, 34, 1328–1342. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  93. Yin, Z.; Raj, D.; Saiepour, N.; Van Dam, D.; Brouwer, N.; Holtman, I.R.; Eggen, B.J.L.; Möller, T.; Tamm, J.A.; Abdourahman, A.; et al. Immune hyperreactivity of Aβ plaque-associated microglia in Alzheimer’s disease. Neurobiol. Aging 2017, 55, 115–122. [Google Scholar] [CrossRef] [PubMed]
  94. Mendes Maia, T.; Gogendeau, D.; Pennetier, C.; Janke, C.; Basto, R. Bug22 influences cilium morphology and the post-translational modification of ciliary microtubules. Biol. Open 2014, 3, 138–151. [Google Scholar] [CrossRef] [Green Version]
  95. Baird, F.J.; Bennett, C.L. Microtubule defects & neurodegeneration. J. Genet. Syndr. Gene 2013, 4, 203. [Google Scholar] [CrossRef] [Green Version]
  96. Brunden, K.R.; Lee, V.M.-Y.; Smith, A.B.; Trojanowski, J.Q.; Ballatore, C. Altered microtubule dynamics in neurodegenerative disease: Therapeutic potential of microtubule-stabilizing drugs. Neurobiol. Dis. 2017, 105, 328–335. [Google Scholar] [CrossRef]
  97. Atamna, H.; Frey, W.H. Mechanisms of mitochondrial dysfunction and energy deficiency in Alzheimer’s disease. Mitochondrion 2007, 7, 297–310. [Google Scholar] [CrossRef]
  98. Rodríguez, J.J.; Verkhratsky, A. Neurogenesis in Alzheimer’s disease. J. Anat. 2011, 219, 78–89. [Google Scholar] [CrossRef]
  99. Hroudová, J.; Singh, N.; Fišar, Z. Mitochondrial dysfunctions in neurodegenerative diseases: Relevance to Alzheimer’s disease. Biomed. Res. Int. 2014, 2014, 175062. [Google Scholar] [CrossRef]
  100. Li, Y.; Sun, H.; Chen, Z.; Xu, H.; Bu, G.; Zheng, H. Implications of GABAergic Neurotransmission in Alzheimer’s Disease. Front. Aging Neurosci. 2016, 8, 31. [Google Scholar] [CrossRef] [Green Version]
  101. Tucsek, Z.; Noa Valcarcel-Ares, M.; Tarantini, S.; Yabluchanskiy, A.; Fülöp, G.; Gautam, T.; Orock, A.; Csiszar, A.; Deak, F.; Ungvari, Z. Hypertension-induced synapse loss and impairment in synaptic plasticity in the mouse hippocampus mimics the aging phenotype: Implications for the pathogenesis of vascular cognitive impairment. GeroScience 2017, 39, 385–406. [Google Scholar] [CrossRef] [PubMed]
  102. Schetters, S.T.T.; Gomez-Nicola, D.; Garcia-Vallejo, J.J.; Van Kooyk, Y. Neuroinflammation: Microglia and T Cells get ready to tango. Front. Immunol. 2018, 8, 1905. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  103. Cao, W.; Zheng, H. Peripheral immune system in aging and Alzheimer’s disease. Mol. Neurodegener. 2018, 13, 51. [Google Scholar] [CrossRef] [PubMed]
  104. Chatterjee, S.; Mudher, A. Alzheimer’s disease and type 2 diabetes: A critical assessment of the shared pathological traits. Front. Neurosci. 2018, 12, 383. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Table 1. Blood-based methylome-wide association results for genes that had significant probes in both brain-specific and blood-based analyses.
Table 1. Blood-based methylome-wide association results for genes that had significant probes in both brain-specific and blood-based analyses.
ProbeIDChrProbePosGeneSNPPosA1FreqPGWASPmQTLbSMRSESMRPSMRPHEIDINHEIDICurrent?Previous?Region?
Plan 2: Only Males
cg0520655919q13.3245913997NANOS2rs6652968745914171A0.1331.83E-042.84E-410.7230.1257.67E-093.28E-0120GGG
cg2567358419q13.3245914293NANOS2rs6652968745914171A0.1331.83E-044.40E-300.8490.1522.45E-081.23E-0120GGG
cg1419229919q13.3245914381NANOS2rs6652968745914171A0.1331.83E-046.71E-420.7180.1247.30E-091.08E-0120GGG
cg1970280219q13.3245914471NANOS2rs6652968745914171A0.1331.83E-043.22E-390.7430.1299.10E-091.03E-0120GGG
Plan 3: Only Females
cg102185466p21.3232762046HLA-DQB2rs776853832762044C0.4266.15E-051.30E-126−0.3040.0603.27E-076.43E-0220SGG
Plan 4: Hypertensive Subjects
cg233957495q35.3177557245FAM193Brs1001530177558514G0.0463.36E-042.55E-26−0.4840.0883.08E-088.77E-025NSS
Plan 5: Non-hypertensive Subjects
cg086313575q32150209647SLC6A7rs10076748150209303A0.1071.77E-031.54E-1930.2880.0563.18E-072.02E-0120NNG
cg238910497q33134679117BPGMrs73441994134679118A0.0214.26E-021.18E-229−0.1560.0301.70E-076.07E-014NSS
cg2463573610q26.13122979534PSTKrs2421140123027854A0.0298.09E-032.67E-77−0.3460.0606.12E-097.16E-018NNN
cg0536084711q13.471576873KRTAP5-11rs1182720871578103T0.0201.70E-033.47E-13−0.9420.1593.50E-092.02E-014NNS
cg1763229913q14.352738831LECT1rs488594752735009C0.0371.23E-037.51E-540.5920.0852.67E-121.34E-0120NGG
cg0955731313q14.352739039LECT1rs488594752735009C0.0371.23E-031.46E-400.6750.1001.37E-111.02E-0120NGG
cg0939729316p13.32005032ZNF598rs727666392005819A0.1741.69E-045.78E-510.6880.1163.06E-092.85E-0120NSG
cg2680489116p13.32005241ZNF598rs112489051999727T0.1814.88E-053.56E-980.5390.0801.62E-117.60E-0220NSG
cg0857618516p13.32005683ZNF598rs727666392005819A0.1741.69E-044.06E-440.7400.1264.76E-093.59E-0120NSG
cg1047020816p13.32008700ZNF598rs10584741998795T0.1816.82E-056.56E-191.1120.2091.02E-077.58E-0214NSG
cg0699836116q2158110599C16orf80rs1044502658109349G0.0695.00E-045.61E-97−0.4420.0691.35E-102.53E-0120NSS
Genomic coordinates are based on Human Genome version 38 (hg38). Chr: chromosomal region (i.e., cytogenetic band); ProbePos: probe position; Gene: the gene or closest gene corresponding to the probe; SNP: top methylation quantitative trait locus (mQTL); Pos: SNP position; A1/Freq: SNP’s effect allele and its frequency; PGWAS: p-value of the SNP in genome-wide association meta-analysis; PmQTL: p-value of the SNP in mQTLs analysis; bSMR, SESMR, and PSMR: beta coefficient, its standard error, and p-value of the probe in summary data-based Mendelian randomization (SMR) test; PHEIDI: p-value of the heterogeneity in dependent instruments (HEIDI) test; NHEIDI: number of single-nucleotide polymorphisms used for HEIDI test; Current?: whether there is any AD-associated SNP within ±1 Mb of the probe in the current genome-wide meta-analysis (N: None, G: SNP with PGWAS < 5E-08, and S: SNP with 5E-08 ≤ PGWAS < 5E-06); Previous?: whether there is any AD-associated SNP within ±1 Mb of the probe in previous GWAS (N: None, G: SNP with PGWAS < 5E-08, and S: SNP with 5E-08 ≤ PGWAS < 5E-06); Region?: whether there is any AD-associated SNP within the chromosomal region (i.e., cytogenetic band) corresponding to the probe (N: None, G: SNP with PGWAS < 5E-08, and S: SNP with 5E-08 ≤ PGWAS < 5E-06).
Table 2. Brain-specific methylome-wide association results for genes that had significant probes in both brain-specific and blood-based analyses.
Table 2. Brain-specific methylome-wide association results for genes that had significant probes in both brain-specific and blood-based analyses.
ProbeIDChrProbePosGeneSNPPosA1FreqPGWASPmQTLbSMRSESMRPSMRPHEIDINHEIDICurrent?Previous?Region?
Plan 2: Only Males
cg0520655919q13.3245913997NANOS2rs6652968745914171G0.8671.83E-045.86E-2980.2720.0432.96E-108.50E-0119GGG
Plan 3: Only Females
cg043221116p21.3232761987HLA-DQB2rs776853832762044A0.5746.15E-050−0.2010.0392.21E-078.61E-0220SGG
cg102185466p21.3232762046HLA-DQB2rs776853832762044A0.5746.15E-050−0.1980.0382.18E-078.32E-0220SGG
Plan 4: Hypertensive Subjects
cg233957495q35.3177557245FAM193Brs1001530177558514A0.9543.36E-042.34E-15−0.7910.1575.17E-071.01E-015NSS
Plan 5: Non-hypertensive Subjects
cg086313575q32150209647SLC6A7rs10076748150209303C0.8931.77E-032.82E-2950.2300.0452.76E-072.24E-0118NNG
cg103086297q33134670051BPGMrs73439998134663724C0.9793.01E-029.28E-48−0.5200.1012.88E-072.57E-013NSS
cg2463573610q26.13122979534PSTKrs13328826122992107A0.9706.26E-032.48E-20−0.3740.0721.68E-078.24E-013NNN
cg1556736011q13.471611653KRTAP5-11rs1182720871578103C0.9801.70E-039.66E-10−0.6790.1301.67E-073.71E-013NNS
cg0955731313q14.352739039LECT1rs488596152755200C0.9604.63E-036.93E-310.5470.1031.06E-075.67E-017NGG
cg0701131816p13.32004943ZNF598rs727666392005819G0.8261.69E-0400.2910.0461.96E-101.12E-0117NSG
cg0939729316p13.32005032ZNF598rs727666392005819G0.8261.69E-0400.2820.0441.86E-101.13E-0118NSG
cg0521118916p13.32005402ZNF598rs115423021986934T0.8197.26E-0500.2830.0437.47E-111.01E-0118NSG
cg0857618516p13.32005683ZNF598rs727666392005819G0.8261.69E-0400.2950.0462.00E-109.02E-0216NSG
cg0699836116q2158110599C16orf80rs7401979058107923T0.9315.00E-044.77E-20−0.5910.1095.49E-086.81E-0111NSS
Genomic coordinates are based on Human Genome version 38 (hg38). Chr: chromosomal region (i.e., cytogenetic band); ProbePos: probe position; Gene: the gene or closest gene corresponding to the probe; SNP: top methylation quantitative trait locus (mQTL); Pos: SNP position; A1/Freq: SNP’s effect allele and its frequency; PGWAS: p-value of the SNP in genome-wide association meta-analysis; PmQTL: p-value of the SNP in mQTLs analysis; bSMR, SESMR, and PSMR: beta coefficient, its standard error, and p-value of the probe in summary data-based Mendelian randomization (SMR) test; PHEIDI: p-value of the heterogeneity in dependent instruments (HEIDI) test; NHEIDI: number of single-nucleotide polymorphisms used for HEIDI test; Current?: whether there is any AD-associated SNP within ±1 Mb of the probe in the current genome-wide meta-analysis (N: None, G: SNP with PGWAS < 5E-08, and S: SNP with 5E-08 ≤ PGWAS < 5E-06); Previous?: whether there is any AD-associated SNP within ±1 Mb of the probe in previous GWAS (N: None, G: SNP with PGWAS < 5E-08, and S: SNP with 5E-08 ≤ PGWAS < 5E-06); Region?: whether there is any AD-associated SNP within the chromosomal region (i.e., cytogenetic band) corresponding to the probe (N: None, G: SNP with PGWAS < 5E-08, and S: SNP with 5E-08 ≤ PGWAS < 5E-06).
Table 3. Methylome-wide association results for the four genes that had epigenetically and transcriptionally AD-associated probes.
Table 3. Methylome-wide association results for the four genes that had epigenetically and transcriptionally AD-associated probes.
ProbeIDChrProbePosGeneSNPPosA1FreqPGWASPmQTLbSMRSESMRPSMRPHEIDINHEIDICurrent?Previous?Region?
Blood-based Analyses
cg030635112p13.173930386DGUOKrs673715673932607C0.0365.62E-032.71E-227−0.2470.0412.74E-091.09E-0111NNN
cg0285071511q24.3130159317ST14rs34008994130165703T0.0961.55E-041.21E-26−0.8120.1384.14E-097.87E-0120NNG
cg2102976911q24.3130159620ST14rs34008994130165703T0.0961.55E-044.09E-18−1.0060.1844.58E-089.16E-0120NNG
cg0699836116q2158110599C16orf80rs1044502658109349G0.0695.00E-045.61E-97−0.4420.0691.35E-102.53E-0120NSS
Brain-specific Analyses
cg110031331q23.1159076601AIM2rs16841642159077008G0.9525.30E-036.30E-82−0.3120.0624.62E-073.40E-0118NSN
cg0699836116q2158110599C16orf80rs7401979058107923T0.9315.00E-044.77E-20−0.5910.1095.49E-086.81E-0111NSS
Genomic coordinates are based on Human Genome version 38 (hg38). Chr: chromosomal region (i.e., cytogenetic band); ProbePos: probe position; Gene: the gene or closest gene corresponding to the probe; SNP: top methylation quantitative trait locus (mQTL); Pos: SNP position; A1/Freq: SNP’s effect allele and its frequency; PGWAS: p-value of the SNP in genome-wide association meta-analysis; PmQTL: p-value of the SNP in mQTLs analysis; bSMR, SESMR, and PSMR: beta coefficient, its standard error, and p-value of the probe in summary data-based Mendelian randomization (SMR) test; PHEIDI: p-value of the heterogeneity in dependent instruments (HEIDI) test; NHEIDI: number of single-nucleotide polymorphisms used for HEIDI test; Current?: whether there is any AD-associated SNP within ±1 Mb of the probe in the current genome-wide meta-analysis (N: None, G: SNP with PGWAS < 5E-08, and S: SNP with 5E-08 ≤ PGWAS < 5E-06); Previous?: whether there is any AD-associated SNP within ±1 Mb of the probe in previous GWAS (N: None, G: SNP with PGWAS < 5E-08, and S: SNP with 5E-08 ≤ PGWAS < 5E-06); Region?: whether there is any AD-associated SNP within the chromosomal region (i.e., cytogenetic band) corresponding to the probe (N: None, G: SNP with PGWAS < 5E-08, and S: SNP with 5E-08 ≤ PGWAS < 5E-06).
Table 4. Pathway-enrichment of blood-based methylome-wide association results.
Table 4. Pathway-enrichment of blood-based methylome-wide association results.
PathwayPathway SourceGSEA IDSizeCountZ-Scorep-Valueq-Value
Plan 1: All Subjects
Type II diabetes mellitusKEGGM1970847144.0172.95E-057.35E-03
MHC class II antigen presentationREACTOMEM70591163.5571.87E-042.33E-02
Host Interactions of HIV factorsREACTOMEM5283132113.2026.81E-045.65E-02
LysosomeKEGGM11266121113.1119.31E-045.80E-02
GABA-B receptor activation REACTOMEM95438103.0081.31E-036.54E-02
L1CAM interactionsREACTOMEM87286172.9871.41E-036.54E-02
Vascular smooth muscle contractionKEGGM9387115222.8522.17E-037.73E-02
Plan 2: Only Males
Neurotransmitter receptors and postsynaptic signal transmissionREACTOMEM752137253.3693.77E-041.02E-01
Transmission across chemical synapses REACTOMEM15514186343.2875.06E-041.02E-01
GABA receptor activationREACTOMEM97652113.0411.18E-031.06E-01
Phospholipase C-mediated cascadeREACTOMEM85654122.7542.94E-031.98E-01
Plan 3: Only Females
GABA-B receptor activation REACTOMEM95438113.6981.09E-042.66E-02
O-linked glycosylation of mucinsREACTOMEM54659103.4183.15E-043.86E-02
GABA receptor activationREACTOMEM97652133.3643.84E-043.86E-02
extracellular matrix (ECM) regulatorsNABAM3468238413.3613.88E-043.86E-02
Plan 5: Non-hypertensive Subjects
Retinoblastoma 1 pathwayPIDM27965103.711.04E-043.20E-02
Circadian clockREACTOMEM93853123.5082.26E-043.48E-02
Alzheimer’s diseaseKEGGM16024169243.0111.30E-031.34E-01
GSEA: Gene Set Enrichment Analysis; Size: number of genes in the pathway; Count: number of enriched genes in the pathway; KEGG: Kyoto Encyclopedia of Genes and Genomes; REACTOME: REACTOME pathway knowledgebase; PID: Pathway Interaction Database; NABA: Matrisome Project. The false discovery rate thresholds were 0.1, 0.2, 0.05, and 0.15 for plans 1, 2, 3, and 5, respectively.
Table 5. Pathway-enrichment of brain-specific methylome-wide association results.
Table 5. Pathway-enrichment of brain-specific methylome-wide association results.
PathwayPathway SourceGSEA IDSizeCountZ-Scorep-Valueq-Value
Plan 1: All Subjects
MHC class II antigen presentationREACTOMEM70591143.34.84E-041.07E-01
Plan 2: Only Males
Ubiquitin mediated proteolysisKEGGM15247138143.1986.91E-041.54E-01
Type II diabetes mellitusKEGGM1970847172.891.93E-032.15E-01
Plan 3: Only Females
MHC class II antigen presentationREACTOMEM70591183.1388.50E-041.56E-01
Transport of inorganic cations/anions and amino acids/oligopeptidesREACTOMEM82394112.8492.19E-032.02E-01
Plan 4: Hypertensive Subjects
DNA repairREACTOMEM15434112103.875.44E-051.26E-02
Type II diabetes mellitusKEGGM1970847103.6221.46E-041.69E-02
Extracellular matrix (ECM) affiliated proteinsNABAM5880171223.0191.27E-039.77E-02
Plan 5: Non-hypertensive Subjects
Respiratory electron transport, ATP synthesis by chemiosmotic coupling, and heat production by uncoupling proteinsREACTOMEM102598103.8515.89E-051.66E-02
Hematopoietic cell lineageKEGGM685688133.0031.33E-031.88E-01
The citric acid (TCA) cycle and respiratory electron transport REACTOMEM516141142.9331.68E-031.88E-01
GSEA: Gene Set Enrichment Analysis; Size: number of genes in the pathway; Count: number of enriched genes in the pathway; KEGG: Kyoto Encyclopedia of Genes and Genomes; REACTOME: REACTOME pathway knowledgebase; PID: Pathway Interaction Database; NABA: Matrisome Project. The false discovery rate thresholds were 0.2, 0.25, 0.25, 0.1, and 0.2 for plans 1–5, respectively.

Share and Cite

MDPI and ACS Style

Nazarian, A.; Yashin, A.I.; Kulminski, A.M. Summary-Based Methylome-Wide Association Analyses Suggest Potential Genetically Driven Epigenetic Heterogeneity of Alzheimer’s Disease. J. Clin. Med. 2020, 9, 1489. https://doi.org/10.3390/jcm9051489

AMA Style

Nazarian A, Yashin AI, Kulminski AM. Summary-Based Methylome-Wide Association Analyses Suggest Potential Genetically Driven Epigenetic Heterogeneity of Alzheimer’s Disease. Journal of Clinical Medicine. 2020; 9(5):1489. https://doi.org/10.3390/jcm9051489

Chicago/Turabian Style

Nazarian, Alireza, Anatoliy I. Yashin, and Alexander M. Kulminski. 2020. "Summary-Based Methylome-Wide Association Analyses Suggest Potential Genetically Driven Epigenetic Heterogeneity of Alzheimer’s Disease" Journal of Clinical Medicine 9, no. 5: 1489. https://doi.org/10.3390/jcm9051489

APA Style

Nazarian, A., Yashin, A. I., & Kulminski, A. M. (2020). Summary-Based Methylome-Wide Association Analyses Suggest Potential Genetically Driven Epigenetic Heterogeneity of Alzheimer’s Disease. Journal of Clinical Medicine, 9(5), 1489. https://doi.org/10.3390/jcm9051489

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop