- Open Access
Genome-wide DNA methylation profiling in the superior temporal gyrus reveals epigenetic signatures associated with Alzheimer’s disease
© Watson et al. 2016
- Received: 19 May 2015
- Accepted: 29 December 2015
- Published: 19 January 2016
Alzheimer’s disease affects ~13 % of people in the United States 65 years and older, making it the most common neurodegenerative disorder. Recent work has identified roles for environmental, genetic, and epigenetic factors in Alzheimer’s disease risk.
We performed a genome-wide screen of DNA methylation using the Illumina Infinium HumanMethylation450 platform on bulk tissue samples from the superior temporal gyrus of patients with Alzheimer’s disease and non-demented controls. We paired a sliding window approach with multivariate linear regression to characterize Alzheimer’s disease-associated differentially methylated regions (DMRs).
We identified 479 DMRs exhibiting a strong bias for hypermethylated changes, a subset of which were independently associated with aging. DMR intervals overlapped 475 RefSeq genes enriched for gene ontology categories with relevant roles in neuron function and development, as well as cellular metabolism, and included genes reported in Alzheimer’s disease genome-wide and epigenome-wide association studies. DMRs were enriched for brain-specific histone signatures and for binding motifs of transcription factors with roles in the brain and Alzheimer’s disease pathology. Notably, hypermethylated DMRs preferentially overlapped poised promoter regions, marked by H3K27me3 and H3K4me3, previously shown to co-localize with aging-associated hypermethylation. Finally, the integration of DMR-associated single nucleotide polymorphisms with Alzheimer’s disease genome-wide association study risk loci and brain expression quantitative trait loci highlights multiple potential DMRs of interest for further functional analysis.
We have characterized changes in DNA methylation in the superior temporal gyrus of patients with Alzheimer’s disease, highlighting novel loci that facilitate better characterization of pathways and mechanisms underlying Alzheimer’s disease pathogenesis, and improve our understanding of epigenetic signatures that may contribute to the development of disease.
- Superior Temporal Gyrus
- GWAS SNPs
- Bivalent Domain
- Illumina Infinium HumanMethylation450
- False Discovery Rate Cutoff
Alzheimer’s disease (AD) is the most common neurodegenerative disorder and the leading cause of dementia in the elderly . Diagnosis of AD is based on the presence of neurofibrillary tangles and amyloid plaques , and symptoms typically include memory loss and impaired cognitive ability. Although the pathological hallmarks associated with dementia-related symptoms in AD appear largely similar between both the early-onset and late-onset forms of the disease, their underlying etiologies contrast . Whereas early-onset AD is a familial autosomal dominant disorder caused by rare, highly penetrant mutations in one of a small set of genes (APP, PSEN1, and PSEN2), the more common late-onset form of the disease (accounting for 90–95 % of cases) occurs sporadically, and risk is determined by complex underlying mechanisms [3–6]. Estimates based on twin concordance rates suggest heritability of late-onset AD is as high as 70 %, implicating major roles for genetic as well as non-genetic factors . Indeed, through candidate gene studies, as well as more recent genome-wide association studies (GWASs) and whole-exome sequencing, both common and rare variants associated with the late-onset form of AD have been identified [7–11]. Collectively, however, common GWAS variants account for only a modest proportion (~30 %) of the underlying variance in disease susceptibility . Several environmental factors are also thought to play a role [5, 6], yet exactly how these contribute to risk, onset, and progression remains poorly defined.
Recently, there has been increasing interest in the role of epigenetic mechanisms in the interaction between the genome and environment in human diseases [13, 14], including AD [15, 16]. Epigenetic alterations can be defined as modifications to DNA that impact gene expression and phenotype without a change in the nucleotide sequence. These changes can arise within cells of an individual and be maintained through mitosis , as well as passed from parent to offspring meiotically . One of the best-studied epigenetic modifications involves changes in DNA methylation at CpG dinucleotides. The establishment of DNA methylation is essential for normal cell development and differentiation , and impacts many key cellular processes, including gene regulation , X chromosome inactivation [21, 22], and genomic imprinting . The application of genome-wide methylation profiling techniques has led to a rapid increase in the characterization of methylation patterns across the genome, and has facilitated the identification of differentially methylated regions (DMRs) associated with evolutionary processes , human aging [25–29], cancer [30, 31], and complex disease [32–35].
Several lines of evidence point to the influence of DNA methylation in AD pathogenesis , including direct connections between AD and DNA methylation that have been observed both globally and at specific loci. For example, differences in tissue-wide methylation patterns in disease-relevant brain regions have been reported in patients with AD compared to controls, as well as in monozygotic twins discordant for AD [36, 37]. Locus-specific examples from targeted gene studies also demonstrate a role for DNA methylation changes in AD and include observed disease-associated differences at ribosomal RNA gene promoters , genome-wide LINE-1 elements , and known AD susceptibility genes . More recently, three epigenome-wide studies (EWAS) have been conducted in AD, collectively generating DNA methylation profiles from three different brain regions of patients with AD [41–43], observing both cross-tissue and tissue-specific effects. Importantly, each of these studies discovered differentially methylated CpGs outside of well-established AD genetic risk loci, highlighting the potential utility of EWAS in the characterization of novel genes and pathways underlying disease processes.
In the present study, we have used the Illumina Infinium HumanMethylation450 array platform to conduct a genome-wide screen of DNA methylation in the superior temporal gyrus (STG) of 34 patients with AD and 34 controls, a brain region recently demonstrated to be a site of significant AD-associated gene dysregulation . Building on previous EWAS in AD, which have primarily focused on single CpG analysis, we applied our recently developed pipeline that aims to identify DMRs harboring multiple statistically significant CpGs exhibiting concordant disease-associated changes in methylation . Using this approach we have identified novel and robust DMRs associated with >400 coding transcripts, many of which have known roles in brain function and AD pathology. Additionally, we show that identified DMRs co-localize with other functional epigenetic signatures in brain tissues, overlap with risk loci identified in AD GWASs and previous EWASs, and harbor expression quantitative trait loci (eQTLs) associated with changes in brain gene expression.
Study subjects and sample preparation
Tissue samples from the STG of 34 patients with confirmed late-onset AD and 34 non-demented controls matched by age of death (AOD), race, and gender were obtained from the Mount Sinai Brain Bank (www.mssm.edu/research/labs/neuropathology-and-brain-banking). Each donor had previously undergone a battery of pathological evaluations, and diagnosis of AD was based on both clinical and neuropathological criteria  (Additional file 1: Table S1). Patients with AD had a mean clinical dementia rating of 3.3, a Braak Stage score average of 5.7 (see  for staging classification), a mean cortical plaque density of 19.8 based on measurements from five cortical brain regions , and a mean AOD of 79.1 years (range, 66–92 years). Controls were determined to have negligible cortical plaque densities (mean, 0.64), and either no evidence of or only mild clinical symptoms of dementia (mean clinical dementia rating, 0.87; mean Braak Stage score, 1.4), with a mean AOD of 80.5 years (range, 66–95 years). Diagnostic and dementia assessment consent procedures were approved by the institutional review boards of Mount Sinai Medical Center, Jewish Home and Hospital, and the JJ Peters VA Medical Center. Consents for brain donation were obtained in writing from the legal next of kin of all donors.
Tissue dissections and sample preparations were carried out following previously published protocols (see [44–47]). Briefly, following dissections, samples were subjected to proteinase K digestion and treatment with RNAse A. Genomic DNA was then isolated using standard phenol/chloroform extraction and ethanol precipitation methods.
Infinium HumanMethylation450 BeadChip processing
One microgram of DNA from each sample was sodium bisulfate-treated using the EZ DNA Methylation Kit (Zymo Research, Irvine, CA, USA) and processed for analysis on the Illumina Infinium HumanMethylation450 (Illumina, San Diego, CA, USA) array platform at the Mount Sinai Icahn School of Medicine genomics core facility (New York, NY, USA). Subjects were distributed across six BeadChips (12 samples/array) taking into account AOD, gender, race, and case–control status to mitigate anomalies resulting from potential batch effects. The GenomeStudio Methylation Module Package (version 1.9, Illumina) was used for initial data processing, allowing for the calculation of methylation values (expressed as β-values, ranging from 0 to 1) and detection P-values for 482,421 individual probes spanning the 22 autosomes and sex chromosomes. Owing to differences in sex chromosome number between males and females, and the fact that our cohort was of mixed gender, only autosomal loci were considered here.
Before proceeding to statistical analysis, data were processed further following the pipeline developed by Huynh et al. . Probes meeting the following criteria known to impact array performance were excluded: (1) those mapping to more than one position in the human reference genome (build NCBI36; hg18) using BSMAP , allowing a maximum of two mismatches and three gaps; and (2) those probes for which a 1000 Genomes Project [49, 50] single nucleotide polymorphism (SNP; minor allele frequency ≥ 0.05) mapped to within 5 base pairs (bp) of the probe-targeted CpG. In addition, on a per sample basis, individual β-values for a given CpG were not considered if their detection P-value was > 0.01. Data for the 461,272 remaining autosomal CpGs passing our exclusion criteria in the 68 individuals were color and background adjusted, and quantile normalized using lumi and methylumi, implemented in R [51, 52] (www.R-project.org). The Beta Mixture Quantile Method, as implemented in BMIQ version 1.3 , was also applied to the data to correct Infinium I/II probe type bias. Plots from principal component analysis using autosomal methylation profiles from the 68 samples did not reveal major batch effects or anomalous samples. Raw and processed data for all samples have been deposited in GEO under accession GSE76105.
Identification of DMRs associated with Alzheimer’s disease in the superior temporal gyrus
We first employed linear regression to delineate disease-specific effects on methylation between cases and controls at each of the 461,272 autosomal CpGs. Linear models were developed to account for various independent variables in addition to disease status, including AOD, gender, race, array/batch, and neuronal/glial cell composition. Post-mortem interval was not included in the test model because there was no significant difference observed between cases and controls (t-test, P = 0.881). Cell proportions were estimated from our bulk tissue samples using the CETS R package, developed from 450 K array profiles of sorted neuronal and glial cell subsets from 59 adult individuals . When evaluating results of linear regression, neuronal proportion (NP) was considered only for 154,874 CpGs previously reported to be differentially methylated between the two cell types . For these CpGs, we used the following model: β = AOD + gender + race + array + NP + disease status. For all remaining CpGs, NP was not included in the model. CpGs exhibiting average increases in methylation within the AD group as compared to controls (based on regression coefficients associated with disease status) were defined as hypermethylated (hyper) and those exhibiting decreases as hypomethylated (hypo).
Significant CpGs were clustered into DMRs using a 1 kilobase (kb) sliding window, modified from the genome-tiling method described by Bock , and previously developed and implemented by Huynh et al. ; a window size of 1 kb was previously shown to be optimal for Illumina Infinium HumanMethylation450 array data based on correlations between methylation levels of closely neighboring CpGs . Fisher’s method was used to combine one-sided regression P-values across neighboring CpGs within a given 1 kb window, while taking into account methylation state (i.e., hyper vs. hypo). The positions of CpGs on either end of a significant window demarcated the coordinates of each DMR. Combined P-values for each window/DMR were corrected using the Benjamini–Hochberg method for false discovery rate (FDR) ; DMRs meeting a 1 % FDR cutoff were used for downstream analyses.
CpG/DMR annotation and overlap with genomic features
Probes/CpGs were annotated based on their overlap with specific genomic features using BedTools version 2.1 . Features included functional RefSeq (hg18) genes and promoters (defined as ± 2 kb from the transcriptional start site, TSS), where CpGs were considered to be intergenic if they overlapped neither gene bodies nor promoters; CpG islands (CpGi; UCSC/hg18 annotation), shores (±2 kb from CpGi), shelves (±2 kb from CpG shores), and sea (not within islands, shores, or shelves); and DNaseI hypersensitivity sites and histone marks (H3K9ac, H3K27ac, H3K27me3, H3K4me1, and H3K4me3) in various human brain datasets generated as part of the ENCODE and REMC projects and curated as described previously [58, 59].
DMR enrichments in specific gene-related and CpGi-related features were tested using the χ-square test, by comparing the proportions of DMR-CpGs within the genomic features to the overlap of these features with a background list of 461,272 autosomal CpGs from the 450 K array. For instances in which DMRs overlapped both a promoter and gene body, the promoter annotation was given precedent. The enrichment of DMRs overlapping ENCODE and REMC datasets was tested using INRICH ; DMRs (hyper and hypo were considered separately) were used as test regions, the curated histone marks were used as target regions, and again the background set of 450 K CpGs was used as the map file. INRICH estimates if the DMRs overlap curated histone marks more than expected by chance; DMRs were permuted within the genome but matched to the associated DMRs in terms of the number of DMR sites and the number of overlapping histone marks. Empirical P-values were estimated based on 10,000 permutations.
RefSeq gene promoters (±2 kb TSSs) overlapped by DMRs were assessed for potential enrichments of defined transcription factor binding sites (TFBS) characterized in human lymphoblastoid cell lines (LCLs; n = 282 TFBS motifs), and a human medulloblastoma cell line (n = 258 TFBS motifs). Specifically, we used LCL TFBS reported by Pique-Regi et al.  characterized using the CENTIPEDE algorithm, after removing binding motifs lacking specifically assigned TFs; medulloblastoma TFBSs were defined as sites where evolutionarily conserved binding motifs of human/mouse/rat TFs overlapped regions of open chromatin based on DNaseI/formaldehyde-assisted isolation of regulatory elements/chromatin immunoprecipitation synthesis [58, 62]. The UCSC tracks used to compile the medulloblastoma dataset were found at http://genome.ucsc.edu/cgi-bin/hgFileUi?db=hg19&g=wgEncodeOpenChromSynth and http://hgdownload.cse.ucsc.edu/goldenPath/hg19/database/tfbsConsSites. Enrichments were calculated by comparing the counts of TFBSs within our set of DMR-associated promoters to the number of counts occurring in a background list of all RefSeq promoters overlapped by all sampled 450 K probes. The significance of enrichments was assessed using Fisher’s exact test (P-values were Bonferonni corrected, with a threshold set to P < 0.01), including only motifs found in at least 5 % of the 276 tested DMR-associated gene promoters. Gene ontology (GO) enrichments for DMR-associated genes were assessed using GOrilla .
Assessing DMRs in the context of GWAS SNPs and brain eQTLs
The GWAS SNPs used were downloaded from the NHGRI GWAS Catalog , using entries under “Alzheimer’s Disease” (access date: December 2014), including SNPs recently reported in a large meta-analysis ; only SNPs with P-values < 10−6 were considered. To generate brain eQTLs, we used the gene expression and genotyping datasets Braincloud  [GEO accession number: GSE30272], NIA/NIH  [GEO accession number: GSE15745], Harvard Brain Tissue Resource Center  [GEO accession number: GSE44772], and UK Brain Expression Consortium  [GEO accession number: GSE46706]. Brain eQTLs were determined using methods published previously . To detect the overlap among the index GWAS AD SNPs and eQTLs we used the regulatory trait concordance (RTC) approach [70, 71]. The RTC method detects the overlap of disease-associated variants with functional SNPs, accounting for the correlation structure in the genome (i.e., linkage disequilibrium, LD). RTC scores range from 0 to 1, with values ≥0.9 indicating likely causal regulatory effects, as demonstrated previously [70, 71]. For downstream analysis described here we consider pairs of AD GWAS SNP–eQTLs with RTC ≥0.9. The enrichment analysis of DMRs with AD-associated eQTLs was conducted using the GoShifter package (https://www.broadinstitute.org/mpg/goshifter/). GoShifter estimates the significance of overlap between trait-associated variants (AD-associated eQTLs) and epigenome annotations (DMRs), by generating null distributions of randomly shifting annotations locally within a tested region. For this analysis we used 10,000 permutations.
Technical validation of CpG methylation
Primers for locus-specific Sequenom MassARRAY EpiTYPER (Sequenom, San Diego, CA, USA) assays were designed using the EpiDesigner primer design software (http://www.epidesigner.com/start3.html). The same bisulfite-converted DNA samples used for array analysis were used for EpiTYPER PCR amplification following manufacturer’s specifications, and post-PCR sample processing and imaging were carried out at the Einstein College of Medicine Genomics Core (New York, NY, USA). β-values generated by the two technologies were compared using the Pearson’s correlation coefficient (r).
Identification of DMRs in the superior temporal gyrus of patients with Alzheimer’s disease
Although many brain regions are affected throughout the progression of AD, an extensive study of gene expression changes associated with late-onset AD severity across 15 brain regions recently found the STG to be a site of significant gene dysregulation , motivating our focus on this specific region in the present study. We conducted genome-wide profiling of DNA methylation in STG bulk tissue samples from 34 patients with AD and 34 non-demented controls. Quality control processing (see “Methods”) ultimately resulted in high-quality methylation data for 461,272 autosomal CpGs in each of the 68 individuals for differential methylation analysis. Before proceeding to tests for differential methylation, we tested the reproducibility of our 450 K array data by assessing the extent of technical variation at CpGs within five genomic regions using independent locus-specific EpiTYPER PCR-based assays in 30–55 individuals from our cohort. Methylation estimates for the six CpGs tested on both platforms were significantly correlated between the two technologies (P < 0.005; Additional file 2: Table S2).
Recent studies of genome-wide DNA methylation have revealed considerable effects of age [25–29], gender , ethnicity , and cellular composition [19, 54, 74, 75]. Although our disease and control samples were relatively well matched for AOD, gender, and ethnicity, our estimations of neuronal versus glial cell proportions revealed a smaller proportion of neuronal cells in our AD samples (AD mean = 0.247; Control mean = 0.303; t-test, P = 0.00099; Additional file 3: Figure S1). This is consistent with histological studies reporting neuronal loss in the brains of patients with AD . Thus, we used multivariate linear regression to delineate significant AD-associated effects on DNA methylation while accounting for potential effects of these variables. It is worth noting that when we repeated our analysis without considering neuronal proportions using a t-test (data not shown), we found that a large fraction of the CpGs that were significantly associated with AD (~62 %) were also found to have significant differences in methylation (>5 %) between neuronal and glial cells , demonstrating the importance of incorporating cell composition information into methylation studies in DNA extracted from bulk tissue.
We first compared results from our linear regression analysis to the top 100 CpGs recently reported from an EWAS conducted in STG tissue in a separate cohort of patients with AD and controls . Of the 96 CpGs also screened in our study, 22 (22.9 %) were found to be differentially methylated with the same directional change in patients with AD in our cohort (P < 0.05, one-tailed). We also assessed the degree of methylation differences associated with disease in our samples compared to the Lunnon et al.  cohort for these same 96 CpGs, and observed a significant correlation between the two datasets (r = 0.34; P = 0.00067; Additional file 4: Figure S2). The extent of replication observed between our two cohorts is comparable to that initially reported by Lunnon et al.  between their samples and other independent cohorts. The incomplete overlap across studies is to be expected given the smaller sample sizes studied to date. Analogous to what has been observed in genetic studies such as GWAS, with increases in cohort sizes, we should expect to see stronger and broader reproducibility of EWAS results.
Globally, population-wide (AD and controls) β-value averages across all 461,272 autosomal CpGs showed a bimodal distribution (Fig. 1b), with the majority of values falling either below 0.2 or above 0.8. The DMRs identified in our study were strongly biased toward hypermethylated changes (increased in AD; hyper-DMRs = 321, hypo-DMRs = 158; Fig. 1c). Given that AD is linked to aging, it is interesting that CpG DNA methylation has also been shown to increase with age in multiple studies of the human brain [25, 28]. Significant AD-associated CpG methylation was also recently reported to independently correlate with age . We further investigated potential links between AD-DMR CpGs and aging in our dataset by assessing the effects of sample AOD on methylation at the top significant CpG within each DMR (hyper, n = 321; hypo, n = 158) in control samples (n = 34; ages = 66–95), again using multivariate linear regression to account for effects of gender, ethnicity, array/batch, and neuronal/glial proportions. Of the hypermethylated DMR-CpGs, ~21.8 % were significantly associated with AOD (P < 0.05), compared to only ~12 % of hypomethylated DMR-CpGs (Additional file 7: Figure S4A). The degree of significance (−log10 P-value) and absolute estimated regression coefficients were also higher on average for hypermethylated DMR-CpGs (Additional file 7: Figure S4B).
Similar to recent reports in AD and other complex diseases [32–35, 41–43], excluding cancer, the average effect of disease state on CpG methylation was modest (Fig. 1c), with an average absolute β-value change of 0.021 at significant CpGs within DMRs. When only the top CpG per DMR with respect to β-value change was considered, this mean difference increased slightly to 0.03. Importantly, however, even modest differences in methylation have been shown to associate with significant alterations in gene expression [32, 34, 41–43].
Top 25 differentially methylated regions associated with Alzheimer’s disease in the superior temporal gyrus
Number of CpGs
Largest β difference
Most significant P valuea
intergenic (nearest gene, TRIM40)
AGPAT1, RNF5, RNF5P1
LOC100507547, PPT2, PPT2-EGFL8, PRRT1
Analysis of DMRs in the context of gene ontology and functional genomic datasets in the human brain
The top DMR in our dataset (Table 1) overlapped the RNA helicase gene MOV10L1. Although little is known about the specific function of this gene in the brain, many other genes associated with top DMRs have reported roles in brain function, such as RNF39, KLK7, DUSP6, NAV2, and NRG2. In the context of AD pathology, the protein KLK7 (Fig. 2b; Table 1) was recently shown to cleave and degrade β-amyloid (Aβ) and mitigate Aβ-mediated toxicity in vitro , possibly consistent with the observed negative correlation between KLK7 expression and AD disease severity . DUSP6 was recently shown to be a target of the AD-associated microRNA miR-125b, exhibiting decreased expression in AD brains; notably, knockdown of DUSP6 in primary hippocampal neurons lead to a significant increase in tau protein hyperphosphorylation, a key hallmark of AD .
It is also interesting that several of the top DMR-associated genes are involved in adiposity, fat distribution, and the synthesis and metabolism of cholesterol and lipids (PRDM16, TBX15, ELOVL1, and AGPAT1) [81, 82]. GO enrichments in categories related to cholesterol/lipid metabolism are among the highest observed for AD risk genes identified by GWAS ; alterations in the expression of related genes have been observed previously in the same cohort studied here [84–86]. Experimental data also highlight the potential importance of genes involved in these processes in AD pathology, with alterations in lipid and cholesterol levels having been observed in the blood, cerebrospinal fluid, and brains of patients with AD [87, 88], and associated with cognitive performance . Furthermore, at the molecular level, both cholesterol and lipids have important roles in modulating the production and aggregation of Aβ via interactions with well-known mediators of AD pathogenesis and risk, such as APP, APOE, PSEN1, and BACE1 [90, 91].
We next assessed the association of our DMRs with functional genome-wide datasets generated in non-diseased brain tissue/cell lines. During development and cellular differentiation, DNA methylation is known to act in concert with chromatin alterations such as histone methylation and acetylation to modify gene expression programs . Furthermore, the occurrence of histone modifications can illuminate genomic regions with functional properties in the context of disease . Thus, to explore whether our DMRs overlapped functional regions relevant to the human brain, we tested for enrichment of AD DMRs in regions with repressive and permissive histone modification profiles (H3K9ac, H3K27ac, H3K27me3, H3K4me1, and H3K4me3) generated from fetal and adult bulk brain tissue, as well as iPS-derived neurons (Fig. 4b; Additional file 9: Table S5). We found significant enrichments for hyper DMRs in poised promoters (corrected P = 0.001), also referred to as bivalent domains, characterized by the occupancy of H3K27me3 and H3K4me3. Bivalent domains are generally thought to take on repressed states, while remaining “poised” for activation, and these regions are known to have important roles in cell development and pluripotency . It is interesting to note that, although CpGs within poised promoters are typically characterized by hypomethylation, increased DNA hypermethylation associated with human aging has been shown to occur preferentially in bivalent domains in various tissues, including the brain [26, 29, 94]; such changes have also been noted in cancer and cell culture, and may suggest that hypermethylation of bivalent domains results in a reduction of cell pluripotency [95, 96]. Given the strong connection between AD and aging, the overlap observed here between DMRs and poised promoters could have important implications for understanding molecular mechanisms underlying disease onset and progression.
DNA methylation is also known to play complex roles in TF binding, in some cases either hindering or facilitating interactions between DNA motifs and proteins . Thus, we tested for enrichment in AD-DMR-associated RefSeq gene promoters (n = 276) of TFBS from two independent datasets, one curated from LCLs, and a second brain dataset consisting of evolutionarily conserved TFBSs residing within regions of open chromatin (active) in a medulloblastoma cell line (see “Methods”). After applying a multiple-testing correction and stringent filter (Bonferroni P < 0.01), we noted 29 and 28 significant TF motif enrichments for the medulloblastoma and LCL datasets, respectively (Fig. 4c). The strongest enrichments were found for motifs of NFAT in the medulloblastoma dataset (fold-enrichment = 3.45, Bonferroni = 0.003), and MAZR in the LCL dataset (fold-enrichment = 2.97; Bonferroni = 2.18 × 10−10); the NFAT transcription factor family, in particular, has demonstrated roles in AD pathology [98, 99]. Several other motifs/TFs known to regulate pathways involved in brain function were also identified, such as PPARG, PPARA, and SP1. Specifically in the context of previous findings in AD, SP1 has been shown to regulate enriched gene sets that exhibit expression changes associated with memory impairment in patients with AD .
Analysis of DMRs in the context of Alzheimer’s disease GWAS SNPs and brain eQTLs
Overlap between loci of differential methylation/expression, methylation QTLs and eQTLs, and GWAS regions has been observed previously in complex disease. Such findings demonstrate that in addition to the ability to identify novel epigenetic signatures underlying risk or disease progression, EWAS data also provide an opportunity to potentially inform the assignment of putative function to genetic variants associated with disease risk, and may help guide functional analyses of GWAS loci . Of the 479 DMRs identified, 15 fell within ±250 kb of a previously reported GWAS SNP, including those associated with CLU, DIP2C, FRMD4A, HLA-DRB1, HLA-DQB1, CTNNA2, and KLK7 (Additional file 10: Table S6); DMRs overlapping promoters of CLU and FRMD4A fell within ±2 kb of a GWAS SNP.
To further investigate potential functional links with AD GWAS, we integrated DMRs and brain eQTLs with AD GWAS regions, using the regulatory trait approach [70, 71]. We identified 129 risk AD loci that were associated with gene expression of at least one transcript at RTC ≥ 0.9. We examined the enrichment of these AD-associated eQTLs with DMRs using GoShifter. There was no significant enrichment with DMRs (hypomethylated or hypermethylated; empirical P > 0.3). However, we found that three AD-associated eQTLs (and SNPs with r 2 > 0.8) fell inside of AD-DMRs (Additional file 11: Table S7). This included eQTLs associated with the expression of AGPAT1, TAP2, and CLU. Additional file 12: Figure S5A–C shows spatial relationships between DMRs, GWAS SNPs, and eQTLs. Two of these DMRs encompass SNPs in LD with a single GWAS SNP (rs111418223) but two distinct eQTL signals, both of which impact the expression of genes in the vicinity of the Human Leukocyte Antigen (HLA) gene region. TAP2 haplotypes have previously been shown to contribute to AD risk via interactions with APOE4 polymorphisms, with speculated involvement of TAP2 in connections between herpes simplex virus-1 infection and AD . Both AGPAT1 and CLU have roles in lipid/cholesterol metabolism. Specifically, variants in AGPAT1 are associated with variation in levels of circulating sphingolipids and phospholipids in human plasma , and impairments of both Agpat1 expression and cholesterol metabolism have also been observed in a rat model of Huntington’s disease . A role for CLU in AD pathology has long been suspected, solidified by the identification of variants contributing to AD risk by multiple independent GWAS . Several of these risk variants have been linked to alterations of CLU expression and alternative splicing in AD [103–105].
We have conducted an epigenome-wide screen in the STG of patients with AD to characterize clusters of CpGs exhibiting concordant disease-associated changes in DNA methylation. After accounting for effects of sample age, gender, ethnicity, and neuronal/glial cell proportions, we identified 479 autosomal DMRs, the majority of which were defined by hypermethylation in AD cases compared to controls. Although the degree of average disease-associated methylation differences was relatively modest within DMRs, this observation is consistent with previous findings in AD and other complex disease [41–43], including those shown to influence gene expression at both the transcript and protein levels [32, 41–43].
An important consideration is that, although we accounted for potential effects of cellular composition differences between AD cases and controls in our characterization of differential methylation, additional investigation will be required to assess whether the methylation signatures observed are cell type-specific; this is in fact a limitation of all AD EWAS studies conducted to date using bulk tissue . Despite such limitations, our systems-level analyses of DMRs and associated genes provided evidence for likely roles of these regions in AD pathology. This included observed enrichments of DMR-associated genes for GO terms related to the development and function of neurons, as well as cellular metabolism, both relevant to known molecular and neurological impairments in AD, and consistent with findings from AD transcriptome studies and GWAS [83–86, 106].
We also found that DMRs were distributed non-randomly in the genome, with biases in their co-localization within gene and CpG island annotations, as well as preferential overlap with specific brain histone methylation signatures and gene promoters harboring brain-relevant TFs. Most notably, we found significant enrichments specifically for hypermethylated DMRs in poised promoters of the adult brain; these loci, characterized by the presence of both H3K27me3 and H3K4me3, mark regulatory regions associated with developmental genes that have key roles in cellular differentiation and pluripotency [93, 95]. Poised promoters have also been shown to preferentially overlap CpGs that become hypermethylated with age [26, 94]. Interestingly, we found a strong bias for CpG hypermethylation among AD-DMRs, and that these CpGs were enriched for age-associated methylation changes in our control samples when compared to CpGs within hypomethylated DMRs. Taken together, these points highlight a potential interplay between disease-associated epigenetic alterations and aging in AD pathology, and suggest more targeted research in this area may be warranted.
Finally, our results indicate that the study of epigenetic signatures can aid in the characterization of novel genomic regions associated with disease, particularly those overlooked by alternative approaches. Future challenges in the field include the development of effective strategies for integrating epigenetic and transcriptomic profiles with genetic datasets, as a means to better understand the roles of different forms of variation in AD .
This work was supported by a New Investigator Research Grant from the Alzheimer’s Association (2012ALZNIRG69983) to AJS; NIH grants DA033660, HG006696, HD073731, and MH097018 to AJS; research grant 6-FY13-92 from the March of Dimes to AJS; NIH grant R01AG050986 to PR; a Brain Behavior Research Foundation grant 20540 to PR; an Alzheimer's Association grant NIRG-340998 to PR; a Veterans Affairs Merit grant BX002395 to PR; and NIH grants AG05138 and HHSN271291300031 to VH. The authors thank Dr. Zachary Kaminsky for useful discussion regarding the estimation of neuronal/glial cell proportions from bulk tissue methylation profiles.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Alzheimer’s Association. Alzheimer’s disease facts and figures. Alzheimers Dement. 2012;2012(8):131–68.Google Scholar
- Braak H, Braak E. Diagnostic criteria for neuropathologic assessment of Alzheimer’s disease. Neurobiol Aging. 1997;18:S85–8.PubMedView ArticleGoogle Scholar
- Bird TD. Genetic aspects of Alzheimer disease. Genet Med. 2008;10:231–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Bertram L, Lill CM, Tanzi RE. The genetics of Alzheimer disease: back to the future. Neuron. 2010;68:270–81.PubMedView ArticleGoogle Scholar
- Migliore L, Coppedè F. Genetics, environmental factors and the emerging role of epigenetics in neurodegenerative diseases. Mutat Res. 2009;667:82–97.PubMedView ArticleGoogle Scholar
- Gatz M, Reynolds CA, Fratiglioni L, Johansson B, Mortimer JA, Berg S, et al. Role of genes and environments for explaining Alzheimer disease. Arch Gen Psychiatry. 2006;63:168–74.PubMedView ArticleGoogle Scholar
- Harold D, Abraham R, Hollingworth P, Sims R, Gerrish A, Hamshere ML, et al. Genome-wide association study identifies variants at CLU and PICALM associated with Alzheimer’s disease. Nat Genet. 2009;41:1088–93.PubMed CentralView ArticlePubMedGoogle Scholar
- Bertram L, Tanzi RE. Genome-wide association studies in Alzheimer’s disease. Hum Mol Genet. 2009;18:R137–45.PubMed CentralPubMedView ArticleGoogle Scholar
- Chapman J, Rees E, Harold D, Ivanov D, Gerrish A, Sims R, et al. A genome-wide study shows a limited contribution of rare copy number variants to Alzheimer’s disease risk. Hum Mol Genet. 2013;22:816–24.PubMed CentralPubMedView ArticleGoogle Scholar
- Lambert J-C, Grenier-Boley B, Harold D, Zelenika D, Chouraki V, Kamatani Y, et al. Genome-wide haplotype association study identifies the FRMD4A gene as a risk locus for Alzheimer’s disease. Mol Psychiatry. 2013;18:461–70.PubMed CentralPubMedView ArticleGoogle Scholar
- Lambert JC, Ibrahim-Verbaas CA, Harold D, Naj AC, Sims R, Bellenguez C, et al. Meta-analysis of 74,046 individuals identifies 11 new susceptibility loci for Alzheimer’s disease. Nat Genet. 2013;45:1452–8.PubMed CentralPubMedView ArticleGoogle Scholar
- Lee SH, Harold D, Nyholt DR, Goddard ME, Zondervan KT, Williams J, et al. Estimation and partitioning of polygenic variation captured by common SNPs for Alzheimer’s disease, multiple sclerosis and endometriosis. Hum Mol Genet. 2013;22:832–41.PubMed CentralPubMedView ArticleGoogle Scholar
- Portela A, Esteller M. Epigenetic modifications and human disease. Nat Biotechnol. 2010;28:1057–68.View ArticlePubMedGoogle Scholar
- Handel AE, Ebers GC, Ramagopalan SV. Epigenetics: molecular mechanisms and implications for disease. Trends Mol Med. 2010;16:7–16.PubMedView ArticleGoogle Scholar
- Mastroeni D, Grover A, Delvaux E, Whiteside C, Coleman PD, Rogers J. Epigenetic mechanisms in Alzheimer’s disease. Neurobiol Aging. 2011;32:1161–80.PubMed CentralPubMedView ArticleGoogle Scholar
- Lunnon K, Mill J. Epigenetic studies in Alzheimer’s disease: current findings, caveats, and considerations for future studies. Am J Med Genet B Neuropsychiatr Genet. 2013;162B(8):789–99.PubMedView ArticleGoogle Scholar
- Probst AV, Dunleavy E, Almouzni G. Epigenetic inheritance during the cell cycle. Nat Rev Mol Cell Biol. 2009;10:192–206.PubMedView ArticleGoogle Scholar
- Guilmatre A, Sharp AJ. Parent of origin effects. Clin Genet. 2012;81:201–9.PubMedView ArticleGoogle Scholar
- Ziller MJ, Gu H, Müller F, Donaghey J, Tsai LT-Y, Kohlbacher O, et al. Charting a dynamic DNA methylation landscape of the human genome. Nature. 2013;500:477–81.PubMed CentralPubMedView ArticleGoogle Scholar
- Bell JT, Pai AA, Pickrell JK, Gaffney DJ, Pique-Regi R, Degner JF, et al. DNA methylation patterns associate with genetic and gene expression variation in HapMap cell lines. Genome Biol. 2011;12:R10.PubMed CentralPubMedView ArticleGoogle Scholar
- Sharp AJ, Stathaki E, Migliavacca E, Brahmachary M, Montgomery SB, Dupre Y, et al. DNA methylation profiles of human active and inactive X chromosomes. Genome Res. 2011;21:1592–600.PubMed CentralView ArticlePubMedGoogle Scholar
- Bala Tannan N, Brahmachary M, Garg P, Borel C, Alnefaie R, Watson CT, et al. DNA methylation profiling in X;autosome translocations supports a role for L1 repeats in the spread of X chromosome inactivation. Hum Mol Genet. 2014;23:1224–36.PubMed CentralPubMedView ArticleGoogle Scholar
- Edwards CA, Ferguson-Smith AC. Mechanisms regulating imprinted genes in clusters. Curr Opin Cell Biol. 2007;19:281–9.PubMedView ArticleGoogle Scholar
- Hernando-Herraez I, Prado-Martinez J, Garg P, Fernandez-Callejo M, Heyn H, Hvilsom C, et al. Dynamics of DNA methylation in recent human and great ape evolution. PLoS Genet. 2013;9:e1003763.PubMed CentralView ArticlePubMedGoogle Scholar
- Numata S, Ye T, Hyde TM, Guitart-Navarro X, Tao R, Wininger M, et al. DNA methylation signatures in development and aging of the human prefrontal cortex. Am J Hum Genet. 2012;90:260–72.PubMed CentralPubMedView ArticleGoogle Scholar
- Rakyan VK, Down TA, Maslau S, Andrew T, Yang T-P, Beyan H, et al. Human aging-associated DNA hypermethylation occurs preferentially at bivalent chromatin domains. Genome Res. 2010;20:434–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Hannum G, Guinney J, Zhao L, Zhang L, Hughes G, Sadda S, et al. resource genome-wide methylation profiles reveal quantitative views of human aging rates. Mol Cell. 2013;49:1–9.View ArticleGoogle Scholar
- Hernandez DG, Nalls MA, Gibbs JR, Arepalli S, van der Brug M, Chong S, et al. Distinct DNA methylation changes highly correlated with chronological age in the human brain. Hum Mol Genet. 2011;20:1164–72.PubMed CentralView ArticlePubMedGoogle Scholar
- Teschendorff AE, Menon U, Gentry-Maharaj A, Ramus SJ, Weisenberger DJ, Shen H, et al. Age-dependent DNA methylation of genes that are suppressed in stem cells is a hallmark of cancer. Genome Res. 2010;20:440–6.PubMed CentralView ArticlePubMedGoogle Scholar
- Heyn H, Sayols S, Moutinho C, Vidal E, Sanchez-Mut JV, Stefansson OA, et al. Linkage of DNA methylation quantitative trait loci to human cancer risk. Cell Rep. 2014;7:331–8.PubMedView ArticleGoogle Scholar
- Stefansson OA, Moran S, Gomez A, Sayols S, Arribas-Jorba C, Sandoval J, et al. A DNA methylation-based definition of biologically distinct breast cancer subtypes. Mol Oncol. 2014;9:1–14.Google Scholar
- Huynh JL, Garg P, Thin TH, Yoo S, Dutta R, Trapp BD, et al. Epigenome-wide differences in pathology-free regions of multiple sclerosis-affected brains. Nat Neurosci. 2014;17:121–30.PubMed CentralPubMedView ArticleGoogle Scholar
- Pidsley R, Viana J, Hannon E, Spiers H, Troakes C, Al-Saraj S, et al. Methylomic profiling of human brain tissue supports a neurodevelopmental origin for schizophrenia. Genome Biol. 2014;15:483.PubMed CentralPubMedView ArticleGoogle Scholar
- Van Eijk KR, de Jong S, Strengman E, Buizer-Voskamp JE, Kahn RS, Boks MP, et al. Identification of schizophrenia-associated loci by combining DNA methylation and gene expression data from whole blood. Eur J Hum Genet. 2014;2013:1–5.Google Scholar
- Liu Y, Aryee MJ, Padyukov L, Fallin MD, Hesselberg E, Runarsson A, et al. Epigenome-wide association data implicate DNA methylation as an intermediary of genetic risk in rheumatoid arthritis. Nat Biotechnol. 2013;31:142–7.PubMed CentralPubMedView ArticleGoogle Scholar
- Mastroeni D, McKee A, Grover A, Rogers J, Coleman PD. Epigenetic differences in cortical neurons from a pair of monozygotic twins discordant for Alzheimer’s disease. PLoS One. 2009;4:e6617.PubMed CentralPubMedView ArticleGoogle Scholar
- Chouliaras L, Mastroeni D, Delvaux E, Grover A, Kenis G, Hof PR, et al. Consistent decrease in global DNA methylation and hydroxymethylation in the hippocampus of Alzheimer’s disease patients. Neurobiol Aging. 2013;34:2091–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Pietrzak M, Rempala G, Nelson PT, Zheng J-J, Hetman M. Epigenetic silencing of nucleolar rRNA genes in Alzheimer’s disease. PLoS One. 2011;6:e22585.PubMed CentralPubMedView ArticleGoogle Scholar
- Bollati V, Galimberti D, Pergoli L, Dalla Valle E, Barretta F, Cortini F, et al. DNA methylation in repetitive elements and Alzheimer disease. Brain Behav Immun. 2011;25:1078–83.PubMed CentralPubMedView ArticleGoogle Scholar
- Yu L, Chibnik LB, Srivastava GP, Pochet N, Yang J, Xu J, et al. Association of brain DNA methylation in SORL1, ABCA7, HLA-DRB5, SLC24A4, and BIN1 with pathological diagnosis of Alzheimer disease. JAMA Neurol. 2014;60612:15–24.Google Scholar
- Bakulski KM, Dolinoy DC, Sartor MA, Paulson HL, Konen JR, Lieberman AP, et al. Genome-wide DNA methylation differences between late-onset Alzheimer’s disease and cognitively normal controls in human frontal cortex. J Alzheimers Dis. 2012;29:571–88.PubMedGoogle Scholar
- Lunnon K, Smith R, Hannon E, De Jager PL, Srivastava G, Volta M, et al. Methylomic profiling implicates cortical deregulation of ANK1 in Alzheimer’s disease. Nat Neurosci. 2014;17:1164–70.PubMed CentralPubMedView ArticleGoogle Scholar
- De Jager PL, Srivastava G, Lunnon K, Burgess J, Schalkwyk LC, Yu L, et al. Alzheimer’s disease: early alterations in brain DNA methylation at ANK1, BIN1, RHBDF2 and other loci. Nat Neurosci. 2014;17:1156–63.PubMed CentralPubMedView ArticleGoogle Scholar
- Haroutunian V, Katsel P, Schmeidler J. Transcriptional vulnerability of brain regions in Alzheimer’s disease and dementia. Neurobiol Aging. 2009;30:561–73.PubMed CentralPubMedView ArticleGoogle Scholar
- Haroutunian V, Perl D. Regional distribution of neuritic plaques in the nondemented elderly and subjects with very mild Alzheimer disease. Arch Neurol. 1998;55:1185–91.PubMedView ArticleGoogle Scholar
- Haroutunian V, Purohit D. Neurofibrillary tangles in nondemented elderly subjects and mild Alzheimer disease. Arch Neurol. 1999;56(June 1999):713–8.PubMedView ArticleGoogle Scholar
- Haroutunian V, Davies P, Vianna C, Buxbaum JD, Purohit DP. Tau protein abnormalities associated with the progression of Alzheimer disease type dementia. Neurobiol Aging. 2007;28:1–7.PubMedView ArticleGoogle Scholar
- Xi Y, Li W. BSMAP: whole genome bisulfite sequence MAPping program. BMC Bioinformatics. 2009;10:232.PubMed CentralPubMedView ArticleGoogle Scholar
- Abecasis GR, Altshuler D, Auton A, Brooks LD, Durbin RM, Gibbs RA, et al. A map of human genome variation from population-scale sequencing. Nature. 2010;467:1061–73.PubMedView ArticleGoogle Scholar
- Abecasis GR, Auton A, Brooks LD, DePristo MA, Durbin RM, Handsaker RE, et al. An integrated map of genetic variation from 1,092 human genomes. Nature. 2012;491:56–65.PubMedView ArticleGoogle Scholar
- Du P, Kibbe WA, Lin SM. lumi: a pipeline for processing Illumina microarray. Bioinformatics. 2008;24:1547–8.PubMedView ArticleGoogle Scholar
- Du P, Zhang X, Huang C-C, Jafari N, Kibbe WA, Hou L, et al. Comparison of Beta-value and M-value methods for quantifying methylation levels by microarray analysis. BMC Bioinformatics. 2010;11:587.PubMed CentralPubMedView ArticleGoogle Scholar
- Teschendorff AE, Marabita F, Lechner M, Bartlett T, Tegner J, Gomez-Cabrero D, et al. A beta-mixture quantile normalization method for correcting probe design bias in Illumina Infinium 450 k DNA methylation data. Bioinformatics. 2013;29:189–96.PubMed CentralPubMedView ArticleGoogle Scholar
- Guintivano J, Aryee MJ, Kaminsky ZA. A cell epigenotype specific model for the correction of brain cellular heterogeneity bias and its application to age, brain region and major depression. Epigenetics. 2013;8:290–302.PubMed CentralView ArticlePubMedGoogle Scholar
- Bock C. Analysing and interpreting DNA methylation data. Nat Rev Genet. 2012;13:705–19.PubMedView ArticleGoogle Scholar
- Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B. 1995;57:289–300.Google Scholar
- Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26:841–2.PubMed CentralPubMedView ArticleGoogle Scholar
- Maurano MT, Humbert R, Rynes E, Thurman RE, Haugen E, Wang H, et al. Systematic localization of common disease-associated variation in regulatory DNA. Science. 2012;337:1190–6.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhu J, Adli M, Zou JY, Verstappen G, Coyne M, Zhang X, et al. Genome-wide chromatin state transitions associated with developmental and environmental cues. Cell. 2013;152:642–54.PubMed CentralPubMedView ArticleGoogle Scholar
- Lee PH, O’Dushlaine C, Thomas B, Purcell SM. INRICH: interval-based enrichment analysis for genome-wide association studies. Bioinformatics. 2012;28:1797–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Pique-Regi R, Degner JF, Pai AA, Gaffney DJ, Gilad Y, Pritchard JK. Accurate inference of transcription factor binding from DNA sequence and chromatin accessibility data. Genome Res. 2011;21:447–55.PubMed CentralPubMedView ArticleGoogle Scholar
- Good PJ, Guyer MS, Kamholz S, Liefer L, Wetterstrand K, Kampa D, et al. The ENCODE ( ENCyclopedia Of DNA Elements ) Project. Science. 2004;306(5696):636–40.View ArticleGoogle Scholar
- Eden E, Navon R, Steinfeld I, Lipson D, Yakhini Z. GOrilla: a tool for discovery and visualization of enriched GO terms in ranked gene lists. BMC Bioinformatics. 2009;10:48.PubMed CentralPubMedView ArticleGoogle Scholar
- Welter D, MacArthur J, Morales J, Burdett T, Hall P, Junkins H, et al. The NHGRI GWAS Catalog, a curated resource of SNP-trait associations. Nucleic Acids Res. 2014;42(Database issue):D1001–6.PubMed CentralPubMedView ArticleGoogle Scholar
- Colantuoni C, Lipska BK, Ye T, Hyde TM, Tao R, Leek JT, et al. Temporal dynamics and genetic control of transcription in the human prefrontal cortex. Nature. 2011;478:519–23.PubMed CentralPubMedView ArticleGoogle Scholar
- Gibbs JR, van der Brug MP, Hernandez DG, Traynor BJ, Nalls MA, Lai S-L, et al. Abundant quantitative trait loci exist for DNA methylation and gene expression in human brain. PLoS Genet. 2010;6:e1000952.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhang B, Gaiteri C, Bodea L-G, Wang Z, McElwee J, Podtelezhnikov AA, et al. Integrated systems approach identifies genetic nodes and networks in late-onset Alzheimer’s disease. Cell. 2013;153:707–20.PubMed CentralView ArticlePubMedGoogle Scholar
- Ramasamy A, Trabzuni D, Guelfi S, Varghese V, Smith C, Walker R, et al. Genetic variability in the regulation of gene expression in ten regions of the human brain. Nat Neurosci. 2014;17:1418–28.PubMed CentralPubMedView ArticleGoogle Scholar
- Roussos P, Mitchell AC, Voloudakis G, Fullard JF, Pothula VM, Tsang J, et al. A role for noncoding variation in schizophrenia. Cell Rep. 2014;9:1417–29.PubMed CentralPubMedView ArticleGoogle Scholar
- Nica AC, Montgomery SB, Dimas AS, Stranger BE, Beazley C, Barroso I, et al. Candidate causal regulatory effects by integration of expression QTLs with complex trait genetic associations. PLoS Genet. 2010;6:e1000895.PubMed CentralPubMedView ArticleGoogle Scholar
- Grundberg E, Small KS, Hedman ÅK, Nica AC, Buil A, Keildson S, et al. Mapping cis- and trans-regulatory effects across multiple tissues in twins. Nat Genet. 2012;44:1084–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Boks MP, Derks EM, Weisenberger DJ, Strengman E, Janson E, Sommer IE, et al. The relationship of DNA methylation with age, gender and genotype in twins and healthy controls. PLoS One. 2009;4:e6767.PubMed CentralPubMedView ArticleGoogle Scholar
- Heyn H, Moran S, Hernando-Herraez I, Sayols S, Gomez A, Sandoval J, et al. DNA methylation contributes to natural human variation. Genome Res. 2013;23:1363–72.PubMed CentralView ArticlePubMedGoogle Scholar
- Koestler DC, Christensen B, Karagas MR, Marsit CJ, Langevin SM, Kelsey KT, et al. Blood-based profiles of DNA methylation predict the underlying distribution of cell types: a validation analysis. Epigenetics. 2013;8:816–26.PubMed CentralPubMedView ArticleGoogle Scholar
- Kozlenkov A, Roussos P, Timashpolsky A, Barbu M, Rudchenko S, Bibikova M, et al. Differences in DNA methylation between human neuronal and glial cells are concentrated in enhancers and non-CpG sites. Nucleic Acids Res. 2014;42:109–27.PubMed CentralPubMedView ArticleGoogle Scholar
- Gómez-Isla T, Hollister R, West H, Mui S, Growdon JH, Petersen RC, et al. Neuronal loss correlates with but exceeds neurofibrillary tangles in Alzheimer’s disease. Ann Neurol. 1997;41:17–24.PubMedView ArticleGoogle Scholar
- Robinson MD, Kahraman A, Law CW, Lindsay H, Nowicka M, Weber LM, et al. Statistical methods for detecting differentially methylated loci and regions. Front Genet. 2014;5:324.PubMed CentralPubMedView ArticleGoogle Scholar
- Shropshire TD, Reifert J, Rajagopalan S, Baker D, Feinstein SC, Daugherty PS. Amyloid β peptide cleavage by kallikrein 7 attenuates fibril growth and rescues neurons from Aβ-mediated toxicity in vitro. Biol Chem. 2014;395:109–18.PubMedView ArticleGoogle Scholar
- Bossers K, Wirz KTS, Meerhoff GF, Essing AHW, van Dongen JW, Houba P, et al. Concerted changes in transcripts in the prefrontal cortex precede neuropathology in Alzheimer’s disease. Brain. 2010;133(Pt 12):3699–723.PubMedView ArticleGoogle Scholar
- Banzhaf-Strathmann J, Benito E, May S, Arzberger T, Tahirovic S, Kretzschmar H, et al. MicroRNA-125b induces tau hyperphosphorylation and cognitive deficits in Alzheimer’s disease. EMBO J. 2014;33:1667–80.PubMed CentralPubMedView ArticleGoogle Scholar
- Ohno Y, Suto S, Yamanaka M, Mizutani Y, Mitsutake S, Igarashi Y, et al. ELOVL1 production of C24 acyl-CoAs is linked to C24 sphingolipid synthesis. Proc Natl Acad Sci USA. 2010;107:18439–44.PubMed CentralPubMedView ArticleGoogle Scholar
- Demirkan A, van Duijn CM, Ugocsai P, Isaacs A, Pramstaller PP, Liebisch G, et al. Genome-wide association study identifies novel loci associated with circulating phospho- and sphingolipid concentrations. PLoS Genet. 2012;8:e1002490.PubMed CentralPubMedView ArticleGoogle Scholar
- Jones L, Holmans PA, Hamshere ML, Harold D, Moskvina V, Ivanov D, et al. Genetic evidence implicates the immune system and cholesterol metabolism in the aetiology of Alzheimer’s disease. PLoS One. 2010;5:e13950.PubMed CentralPubMedView ArticleGoogle Scholar
- Katsel P, Li C, Haroutunian V. Gene Expression alterations in the sphingolipid metabolism pathways during progression of dementia and Alzheimer’s Disease: a shift toward ceramide accumulation at the earliest recognizable stages of Alzheimer’s disease? Neurochem Res. 2007;32:845–56.PubMedView ArticleGoogle Scholar
- Akram A, Schmeidler J, Katsel P, Hof PR, Haroutunian V. Increased expression of RXR a in dementia: an early harbinger for the cholesterol dyshomeostasis. Mol Neurodegener. 2010;5:1–14.View ArticleGoogle Scholar
- Akram A, Schmeidler J, Katsel P, Hof PR, Haroutunian V. Association of ApoE and LRP mRNA levels with dementia and AD neuropathology. Neurobiol Aging. 2012;33:628.e1–628.e14.View ArticleGoogle Scholar
- Cutler RG, Kelly J, Storie K, Pedersen WA, Tammara A, Hatanpaa K, et al. Involvement of oxidative stress-induced abnormalities in ceramide and cholesterol metabolism in brain aging and Alzheimer’s disease. Proc Natl Acad Sci USA. 2004;101(Track I):1–6.Google Scholar
- He X, Huang Y, Li B, Gong C-X, Schuchman EH. Deregulation of sphingolipid metabolism in Alzheimer’s disease. Neurobiol Aging. 2010;31:398–408.PubMed CentralPubMedView ArticleGoogle Scholar
- Mielke MM, Haughey NJ, Bandaru VVR, Zetterberg H, Blennow K, Andreasson U, et al. Cerebrospinal fluid sphingolipids, β-amyloid, and tau in adults at risk for Alzheimer’s disease. Neurobiol Aging. 2014;35:2486–94.PubMed CentralPubMedView ArticleGoogle Scholar
- Foley P. Lipids in Alzheimer’s disease: a century-old story. Biochim Biophys Acta. 2010;1801:750–3.PubMedView ArticleGoogle Scholar
- Di Paolo G, Kim T-W. Linking lipids to Alzheimer’s disease: cholesterol and beyond. Nat Rev Neurosci. 2011;12:284–96.PubMed CentralPubMedView ArticleGoogle Scholar
- Cedar H, Bergman Y. Linking DNA methylation and histone modification: patterns and paradigms. Nat Rev Genet. 2009;10:295–304.PubMedView ArticleGoogle Scholar
- Bernstein BE, Mikkelsen TS, Xie X, Kamal M, Huebert DJ, Cuff J, et al. A bivalent chromatin structure marks key developmental genes in embryonic stem cells. Cell. 2006;125:315–26.View ArticlePubMedGoogle Scholar
- Watson CT, Disanto G, Sandve GK, Breden F, Giovannoni G, Ramagopalan SV. Age-associated hyper-methylated regions in the human brain overlap with bivalent chromatin domains. PLoS One. 2012;7:e43840.PubMed CentralPubMedView ArticleGoogle Scholar
- Meissner A, Mikkelsen TS, Gu H, Wernig M, Hanna J, Sivachenko A, et al. Genome-scale DNA methylation maps of pluripotent and differentiated cells. Nature. 2008;454:766–70.PubMed CentralPubMedGoogle Scholar
- Ohm JE, Baylin SB. Stem cell chromatin patterns: an instructive mechanism for DNA hypermethylation? Cell Cycle. 2014;6:1040–3.View ArticleGoogle Scholar
- Hu S, Wan J, Su Y, Song Q, Zeng Y, Nguyen HN, et al. DNA methylation presents distinct binding sites for human transcription factors. Elife. 2013;2:e00726.PubMed CentralPubMedGoogle Scholar
- Hudry E, Wu HY, Arbel-Ornath M, Hashimoto T, Matsouaka R, Fan Z, et al. Inhibition of the NFAT pathway alleviates amyloid β neurotoxicity in a mouse model of Alzheimer's disease. J Neurosci. 2012;32:3176–92.PubMed CentralPubMedView ArticleGoogle Scholar
- Abdul HM, Sama MA, Furman JL, Mathis DM, Beckett TL, Weidner AM, et al. Cognitive decline in Alzheimer’s disease is associated with selective changes in calcineurin/NFAT signaling. J Neurosci. 2012;29:12957–69.View ArticleGoogle Scholar
- Ramanan VK, Kim S, Holohan K, Shen L, Nho K, Risacher SL, et al. Genome-wide pathway analysis of memory impairment in the Alzheimer’s Disease Neuroimaging Initiative (ADNI) cohort implicates gene candidates, canonical pathways, and networks. Brain Imaging Behav. 2012;6:634–48.PubMed CentralPubMedView ArticleGoogle Scholar
- Bullido MJ, Martínez-García A, Artiga MJ, Aldudo J, Sastre I, Gil P, et al. A TAP2 genotype associated with Alzheimer’s disease in APOE4 carriers. Neurobiol Aging. 2007;28:519–23.PubMedView ArticleGoogle Scholar
- Cong W, Cai H, Wang R, Daimon CM, Maudsley S, Raber K, et al. Altered hypothalamic protein expression in a rat model of Huntington’s disease. PLoS One. 2012;7:e47240.PubMed CentralPubMedView ArticleGoogle Scholar
- Yu J-T, Tan L. The Role of clusterin in Alzheimer’s disease: pathways, pathogenesis, and therapy. Mol Neurobiol. 2012;45:314–26.PubMedView ArticleGoogle Scholar
- Szymanski M, Wang R, Bassett SS, Avramopoulos D. Alzheimer’s risk variants in the clusterin gene are associated with alternative splicing. Transl Psychiatry. 2011;1:e18–7.PubMed CentralPubMedView ArticleGoogle Scholar
- Schürmann B, Wiese B, Bickel H, Weyerer S, Riedel-Heller SG, Pentzek M, et al. Association of the Alzheimer’s disease clusterin risk allele with plasma clusterin concentration. J Alzheimers Dis. 2011;25:421–4.PubMedGoogle Scholar
- Miller JA, Woltjer RL, Goodenbour JM, Horvath S, Geschwind DH. Genes and pathways underlying regional and cell type changes in Alzheimer’s disease. Genome Med. 2013;5:48.PubMed CentralPubMedView ArticleGoogle Scholar