CoChat
Sign Up
Cited by 253Epigenetics and DNA MethylationAlzheimer's disease research and treatmentsDementia and Cognitive Impairment Research

The DOI matches its Crossref record; 80 of the 80 references it cites were found in a registry or an index.

Epigenetic dysregulation of enhancers in neurons is associated with Alzheimer’s disease pathology and cognitive symptoms

Open access
Li, Peipei, Marshall, Lee, Oh, Gabriel, Jakubowski, Jennifer L., Groot, Daniel, He, Yu, Wang, Ting, Petronis, Arturas, Labrie, Viviane·2019·Nature Communications
DOI 10.1038/s41467-019-10101-7PMC PMC6529540
Open this paper in CoChat — free accountOpen PDF

Abstract

Epigenetic control of enhancers alters neuronal functions and may be involved in Alzheimer’s disease (AD). Here, we identify enhancers in neurons contributing to AD by comprehensive fine-mapping of DNA methylation at enhancers, genome-wide. We examine 1.2 million CpG and CpH sites in enhancers in prefrontal cortex neurons of individuals with no/mild, moderate, and severe AD pathology ( n = 101). We identify 1224 differentially methylated enhancer regions; most of which are hypomethylated at CpH sites in AD neurons. CpH methylation losses occur in normal aging neurons, but are accelerated in AD. Integration of epigenetic and transcriptomic data demonstrates a pro-apoptotic reactivation of the cell cycle in post-mitotic AD neurons. Furthermore, AD neurons have a large cluster of significantly hypomethylated enhancers in the DSCAML1 gene that targets BACE1 . Hypomethylation of these enhancers in AD is associated with an upregulation of BACE1 transcripts and an increase in amyloid plaques, neurofibrillary tangles, and cognitive decline.

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1038/s41467-019-10101-7

Cites

80 of 80 references verified0 citation errors · 0 not checkable

Statistics

Figures

Fig. 1 Fine-mapping of DNA methylation changes at enhancers, genome-wide, in AD neurons. In prefrontal cortex neurons of individuals with no/mild, moderate, and severe AD neurofibrillary tangle pathology (Braak stage 1–6; n = 101 individuals), DNA methylation was profiled at enhancer and promoter regions across the genome using bisulfite padlock probe sequencing. a Manhattan plot showing differentially methylated regions associated with the severity of neurofibrillary tangle pathology, after controlling for sex, age, postmortem interval, and neuronal subtype proportion. −log10(P) refers to the significance of differentially methylated cytosine enrichment at enhancers, with the sign corresponding to the direction of methylation change (hyper or hypomethylated). Threshold for genome-wide significance (dashed line) is q < 0.05. Enhancers in DSCAML1 highlighted in red. b Image of ETS motif, the most strongly overrepresented TF motif within the differentially methylated enhancers in AD (p < 10−18, binomial distribution model in oPOSSUM). Motif sequence logo is provided by MotifMap80. c Contribution of CpG and CpH sites to DNA methylation differences at enhancers in AD. Motif analysis depicts a predominant contribution of CpH sites, especially CpA. d DNA methylation changes at enhancer regions in the DSCAML1 gene in neurons that correspond to the severity of AD pathology. Dashed red line is genome-wide significance threshold (q < 0.05). Red dots denote enhancers that were also differentially methylated before the arrival of neurofibrillary tangles in the prefrontal cortex (Braak 1–4). Tracks for neuronal (NeuN + ) H3K27ac and H3K4me3, from PsychENCODE (pink, n = 9 individuals) are shown. e Averaged % DNA methylation across Braak stage for the most significantly disrupted enhancer region in DSCAML1 (chr11: 117,504,514–117,506,898) in AD neurons. F = female; M = male. *q < 10−14, robust linear regression model followed by hypergeometric test. The boxplot center line is the median, the lower and upper limits are the first and third quartiles (25th and 75th percentiles), and the whiskers are 1.5× the interquartile range
Fig. 1 Fine-mapping of DNA methylation changes at enhancers, genome-wide, in AD neurons. In prefrontal cortex neurons of individuals with no/mild, moderate, and severe AD neurofibrillary tangle pathology (Braak stage 1–6; n = 101 individuals), DNA methylation was profiled at enhancer and promoter regions across the genome using bisulfite padlock probe sequencing. a Manhattan plot showing differentially methylated regions associated with the severity of neurofibrillary tangle pathology, after controlling for sex, age, postmortem interval, and neuronal subtype proportion. −log10(P) refers to the significance of differentially methylated cytosine enrichment at enhancers, with the sign corresponding to the direction of methylation change (hyper or hypomethylated). Threshold for genome-wide significance (dashed line) is q < 0.05. Enhancers in DSCAML1 highlighted in red. b Image of ETS motif, the most strongly overrepresented TF motif within the differentially methylated enhancers in AD (p < 10−18, binomial distribution model in oPOSSUM). Motif sequence logo is provided by MotifMap80. c Contribution of CpG and CpH sites to DNA methylation differences at enhancers in AD. Motif analysis depicts a predominant contribution of CpH sites, especially CpA. d DNA methylation changes at enhancer regions in the DSCAML1 gene in neurons that correspond to the severity of AD pathology. Dashed red line is genome-wide significance threshold (q < 0.05). Red dots denote enhancers that were also differentially methylated before the arrival of neurofibrillary tangles in the prefrontal cortex (Braak 1–4). Tracks for neuronal (NeuN + ) H3K27ac and H3K4me3, from PsychENCODE (pink, n = 9 individuals) are shown. e Averaged % DNA methylation across Braak stage for the most significantly disrupted enhancer region in DSCAML1 (chr11: 117,504,514–117,506,898) in AD neurons. F = female; M = male. *q < 10−14, robust linear regression model followed by hypergeometric test. The boxplot center line is the median, the lower and upper limits are the first and third quartiles (25th and 75th percentiles), and the whiskers are 1.5× the interquartile rangeCC BY · PubMed Central
Fig. 2 Chromatin conformation analysis in cortical neurons identifying the gene targets of the differentially methylated enhancers in AD. Chromatin interactions, genome-wide, were determined using Hi-C data from human prefrontal cortex. a Circos plot showing significant interactions between differentially methylated enhancers in AD neurons and gene promoters (TSS ± 2 kb). AD enhancers acted in cis for all their target gene promoters (average 587 kb ± 37 kb enhancer–target promoter distance; n = 1224 enhancers). b, c Pathways affected by genes with differentially methylated enhancers in AD, as determined by MetaCore. Top ten GO processes and disease pathways show that in our unbiased, genome-wide analysis, the enhancers dysregulated in AD neurons affected genes involved in neurogenesis and amyloid neuropathies (q < 0.05, hypergeometric distribution; red dashed line). d Chromatin interactions of the enhancers in DSCAML1 that were differentially methylated in AD neurons (all interactions ± 400 kb of DSCAML1 shown). Enhancers in DSCAML1 targeted the promoters of BACE1 (shown in red) as well as other genes
Fig. 2 Chromatin conformation analysis in cortical neurons identifying the gene targets of the differentially methylated enhancers in AD. Chromatin interactions, genome-wide, were determined using Hi-C data from human prefrontal cortex. a Circos plot showing significant interactions between differentially methylated enhancers in AD neurons and gene promoters (TSS ± 2 kb). AD enhancers acted in cis for all their target gene promoters (average 587 kb ± 37 kb enhancer–target promoter distance; n = 1224 enhancers). b, c Pathways affected by genes with differentially methylated enhancers in AD, as determined by MetaCore. Top ten GO processes and disease pathways show that in our unbiased, genome-wide analysis, the enhancers dysregulated in AD neurons affected genes involved in neurogenesis and amyloid neuropathies (q < 0.05, hypergeometric distribution; red dashed line). d Chromatin interactions of the enhancers in DSCAML1 that were differentially methylated in AD neurons (all interactions ± 400 kb of DSCAML1 shown). Enhancers in DSCAML1 targeted the promoters of BACE1 (shown in red) as well as other genesCC BY · PubMed Central
Fig. 3 Omics analysis discovered high-confidence networks involved in AD. Using the FOREST OmicsIntegrator77, we merged our DNA methylation sequencing dataset profiling enhancers in AD neurons and transcriptomic sequencing data generated from the no/mild, moderate, and severe AD prefrontal cortex (Braak stage 1–6). Molecular pathways affected in AD neurons, as determined by network analysis are shown, and highlight two major networks involved in pathological progression in AD neurons. Pathway analysis by MetaCore of these two networks show their primary involvement in cell cycle reentry and amyloid neuropathies, respectively (q < 0.001, hypergeometric distribution). Hub genes: UBC and CUL3; subnetwork genes connected to CUL3: APP, CTNNB1, PIK3R1, ESR1, HSP90AB1, and PPP2CA are shown in black
Fig. 3 Omics analysis discovered high-confidence networks involved in AD. Using the FOREST OmicsIntegrator77, we merged our DNA methylation sequencing dataset profiling enhancers in AD neurons and transcriptomic sequencing data generated from the no/mild, moderate, and severe AD prefrontal cortex (Braak stage 1–6). Molecular pathways affected in AD neurons, as determined by network analysis are shown, and highlight two major networks involved in pathological progression in AD neurons. Pathway analysis by MetaCore of these two networks show their primary involvement in cell cycle reentry and amyloid neuropathies, respectively (q < 0.001, hypergeometric distribution). Hub genes: UBC and CUL3; subnetwork genes connected to CUL3: APP, CTNNB1, PIK3R1, ESR1, HSP90AB1, and PPP2CA are shown in blackCC BY · PubMed Central
Fig. 4 DNA methylation at DSCAML1 intron 3 is linked to BACE1 gene expression and to the pathology and cognitive symptoms of AD. a Correlation of DNA methylation at a DSCAML1 CpG site (cg07533617) with BACE1 mRNA expression in the AD prefrontal cortex. DNA methylation status at cg07533617 (situated in DSCAML1 intron 3 at an enhancer interacting with the BACE1 gene promoter) was inversely correlated with BACE1 mRNA expression during the early stages of AD (Braak stage ≤ 4, n = 101 AD patients; p < 0.005, robust linear regression) but not in late stage (Braak stage ≥ 5, n = 76 AD patients; p = 0.21, robust linear regression). Analysis performed using the ROSMAP dataset, which contained both transcriptomic data generated by RNA-seq and genome-wide DNA methylation data generated by 450K Illumina DNA methylation arrays from the same AD patients. Correlation between residuals of DNA methylation (beta values of probe cg07533617, x-axis) and BACE1 mRNA levels (FPKM, y-axis) shown. b Significance of correlation between DNA methylation at cg07533617 in DSCAML1 intron 3 and AD pathology/cognitive symptoms; determined using ROSMAP DNA methylation and pathological/clinical data (n = 465 AD and controls; 251 AD and 214 controls). DNA hypomethylation at cg07533617 in DSCAML1 intron 3 was significantly correlated with increased amyloid pathology and neurofibrillary tangle density, as well as a decline in episodic memory, perceptual speed, and global cognitive function (*p < 0.05, linear mixed model with annual cognitive measures as the longitudinal outcomes and DNA methylation as the predictor). a, b All analyses of ROSMAP data adjusted for sex, age, postmortem interval, years of education, neuronal cell proportion, as well as in a RIN
Fig. 4 DNA methylation at DSCAML1 intron 3 is linked to BACE1 gene expression and to the pathology and cognitive symptoms of AD. a Correlation of DNA methylation at a DSCAML1 CpG site (cg07533617) with BACE1 mRNA expression in the AD prefrontal cortex. DNA methylation status at cg07533617 (situated in DSCAML1 intron 3 at an enhancer interacting with the BACE1 gene promoter) was inversely correlated with BACE1 mRNA expression during the early stages of AD (Braak stage ≤ 4, n = 101 AD patients; p < 0.005, robust linear regression) but not in late stage (Braak stage ≥ 5, n = 76 AD patients; p = 0.21, robust linear regression). Analysis performed using the ROSMAP dataset, which contained both transcriptomic data generated by RNA-seq and genome-wide DNA methylation data generated by 450K Illumina DNA methylation arrays from the same AD patients. Correlation between residuals of DNA methylation (beta values of probe cg07533617, x-axis) and BACE1 mRNA levels (FPKM, y-axis) shown. b Significance of correlation between DNA methylation at cg07533617 in DSCAML1 intron 3 and AD pathology/cognitive symptoms; determined using ROSMAP DNA methylation and pathological/clinical data (n = 465 AD and controls; 251 AD and 214 controls). DNA hypomethylation at cg07533617 in DSCAML1 intron 3 was significantly correlated with increased amyloid pathology and neurofibrillary tangle density, as well as a decline in episodic memory, perceptual speed, and global cognitive function (*p < 0.05, linear mixed model with annual cognitive measures as the longitudinal outcomes and DNA methylation as the predictor). a, b All analyses of ROSMAP data adjusted for sex, age, postmortem interval, years of education, neuronal cell proportion, as well as in a RINCC BY · PubMed Central
Fig. 5 Genetic variation within DSCAML1 intron 3 enhancers affects BACE1 mRNA levels in the prefrontal cortex. ROSMAP data containing AD patients and controls (n = 278 individuals) that had both genotype data (using genome-wide SNP arrays) and transcriptome analysis (by RNA-seq). The extended genomic area around BACE1 (±500 kb) is shown. Haplotypes (n = 53) determined by Haploview. Analysis of haplotype association with BACE1 expression that controlled for age, sex, postmortem interval, Braak stage, years of education, RIN, and neuronal cell proportion was performed. Pink track shows the significance of the haplotype association with BACE1 expression (dashed red line is q < 0.05). BACE1 expression in the prefrontal cortex was influenced by two haplotypes in DSCAML1 (q < 10−4, robust linear regression model). These two haplotypes overlapped four enhancers determined to be epigenetically dysregulated in AD neurons, supporting their involvement in BACE1 regulation. Three haplotypes near the BACE1 TSS were also associated with BACE1 expression. Haplotype information (including SNP ID) and their population frequencies are shown for the five haplotypes significantly associated with BACE1 expression. D’ values, measure of linkage disequilibrium, are shown in the boxes, with red boxes indicative of a strong linkage disequilibrium between SNPs
Fig. 5 Genetic variation within DSCAML1 intron 3 enhancers affects BACE1 mRNA levels in the prefrontal cortex. ROSMAP data containing AD patients and controls (n = 278 individuals) that had both genotype data (using genome-wide SNP arrays) and transcriptome analysis (by RNA-seq). The extended genomic area around BACE1 (±500 kb) is shown. Haplotypes (n = 53) determined by Haploview. Analysis of haplotype association with BACE1 expression that controlled for age, sex, postmortem interval, Braak stage, years of education, RIN, and neuronal cell proportion was performed. Pink track shows the significance of the haplotype association with BACE1 expression (dashed red line is q < 0.05). BACE1 expression in the prefrontal cortex was influenced by two haplotypes in DSCAML1 (q < 10−4, robust linear regression model). These two haplotypes overlapped four enhancers determined to be epigenetically dysregulated in AD neurons, supporting their involvement in BACE1 regulation. Three haplotypes near the BACE1 TSS were also associated with BACE1 expression. Haplotype information (including SNP ID) and their population frequencies are shown for the five haplotypes significantly associated with BACE1 expression. D’ values, measure of linkage disequilibrium, are shown in the boxes, with red boxes indicative of a strong linkage disequilibrium between SNPsCC BY · PubMed Central
Fig. 6 Accelerated CpH methylation changes with aging in AD neurons. Age-dependent changes in CpH methylation in neurons from individuals with no/mild, moderate, and severe AD pathology (Braak stage 1–2, 3–4, 5–6, and n = 38 individuals, 32, and 31 individuals, respectively). Age analysis was performed on CpH sites (n = 3661 CpHs) in enhancers relevant to AD (n = 848 enhancers showing both epigenetic and associated gene transcript differences in AD, as identified by Omics integration analysis). a Aging changes in % CpH methylation at enhancers in neurons. CpH methylation decreased with age in individuals with no/mild AD pathology (p < 0.005, robust linear regression model), but not in moderate or severe AD. Analysis of % CpH methylation changes with age adjusted for sex, postmortem interval, and neuronal proportion. b Box plot showing the difference between the epigenetic age and chronological age in neurons of moderate and severe AD cases. The epigenetic aging calculator17 used CpH sites to determine the epigenetic ages of the moderate and severe AD groups. CpH methylation aging was accelerated in severe AD neurons (p < 0.05, paired t-test). The boxplot center line is the median, the lower and upper limits are the first and third quartiles (25th and 75th percentiles), and the whiskers are 1.5× the interquartile range. c CpH methylation age relative to chronological age in neurons of severe AD cases (Braak stage 5–6). Each dot represents a severe AD case, and the line demarks where chronological age matches epigenetic age. Purple dots are samples in which epigenetic age is greater than chronological age. Severe AD cases typically show an epigenetic age of 80 years or older regardless of chronological age
Fig. 6 Accelerated CpH methylation changes with aging in AD neurons. Age-dependent changes in CpH methylation in neurons from individuals with no/mild, moderate, and severe AD pathology (Braak stage 1–2, 3–4, 5–6, and n = 38 individuals, 32, and 31 individuals, respectively). Age analysis was performed on CpH sites (n = 3661 CpHs) in enhancers relevant to AD (n = 848 enhancers showing both epigenetic and associated gene transcript differences in AD, as identified by Omics integration analysis). a Aging changes in % CpH methylation at enhancers in neurons. CpH methylation decreased with age in individuals with no/mild AD pathology (p < 0.005, robust linear regression model), but not in moderate or severe AD. Analysis of % CpH methylation changes with age adjusted for sex, postmortem interval, and neuronal proportion. b Box plot showing the difference between the epigenetic age and chronological age in neurons of moderate and severe AD cases. The epigenetic aging calculator17 used CpH sites to determine the epigenetic ages of the moderate and severe AD groups. CpH methylation aging was accelerated in severe AD neurons (p < 0.05, paired t-test). The boxplot center line is the median, the lower and upper limits are the first and third quartiles (25th and 75th percentiles), and the whiskers are 1.5× the interquartile range. c CpH methylation age relative to chronological age in neurons of severe AD cases (Braak stage 5–6). Each dot represents a severe AD case, and the line demarks where chronological age matches epigenetic age. Purple dots are samples in which epigenetic age is greater than chronological age. Severe AD cases typically show an epigenetic age of 80 years or older regardless of chronological ageCC BY · PubMed Central
Fig. 7 Proposed model displaying the contribution of enhancer dysregulation to the development and progression of AD. In aging, neurons gradually lose CpH methylation at enhancers, especially at enhancers regulating genes involved in the early cell cycle. Aging and other AD risk factors that trigger improper attempts of neurons to renter the cell cycle promotes the formation of neurofibrillary tangles in the AD brain and is neurotoxic. In AD neurons there is also a hypomethylation of enhancers in DSCAML1 that occurs early in disease. Enhancers in DSCAML1 activate BACE1 expression, which leads to the increased production of the Aβ peptides. Aβ peptides promotes an accelerated loss of DNA methylation at enhancers in AD neurons, thereby facilitating aberrant cell cycle reactivation. Aβ peptides also forms plaques, which fuels the spread of neurofibrillary tangle pathology across the brain. Consequently, epigenetic changes at enhancers in AD neurons contributes to the pathology, neurodegeneration, and cognitive symptoms of AD
Fig. 7 Proposed model displaying the contribution of enhancer dysregulation to the development and progression of AD. In aging, neurons gradually lose CpH methylation at enhancers, especially at enhancers regulating genes involved in the early cell cycle. Aging and other AD risk factors that trigger improper attempts of neurons to renter the cell cycle promotes the formation of neurofibrillary tangles in the AD brain and is neurotoxic. In AD neurons there is also a hypomethylation of enhancers in DSCAML1 that occurs early in disease. Enhancers in DSCAML1 activate BACE1 expression, which leads to the increased production of the Aβ peptides. Aβ peptides promotes an accelerated loss of DNA methylation at enhancers in AD neurons, thereby facilitating aberrant cell cycle reactivation. Aβ peptides also forms plaques, which fuels the spread of neurofibrillary tangle pathology across the brain. Consequently, epigenetic changes at enhancers in AD neurons contributes to the pathology, neurodegeneration, and cognitive symptoms of ADCC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

Cited by

Sorted by citations

See all 253

Related

See all

Metadata

Authors
Li, P., Marshall, L., Oh, G., Jakubowski, J. L., Groot, D., He, Y., Wang, T., Petronis, A., & Labrie, V.
Affiliations
  • Van Andel Institute (US)
  • Centre for Addiction and Mental Health (CA)
  • Washington University in St. Louis (US)
Corresponding author
Labrie, Viviane
ORCID
Venue
Nature Communications
Publisher
Springer Science and Business Media LLC
Published
2019-05-21
Volume
10
Issue
1
Pages
2246-2246
ISSN
2041-1723 (ISSN-L)
DOI
10.1038/s41467-019-10101-7
PMC
PMC6529540
OpenAlex
W2945063373
Cited by
234
References
80
Citations by year
2026: 15 · 2025: 22 · 2024: 31 · 2023: 30 · 2022: 41 · 2021: 51 · 2020: 34 · 2019: 10
FWCI
12.73
Type
journal article
Language
en
Open access
Licence
Locations
  • Nature Communications (publishedVersion, cc-by)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • Digital Commons@Becker (Washington University School of Medicine) (submittedVersion, cc-by)
  • RePEc: Research Papers in Economics (submittedVersion)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion, cc-by-sa)
Funders
  • Scottish Rite Charitable Foundation of Canada (15110)
  • U.S. Department of Defense (PD170089)
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Epigenetics and DNA Methylation
  • Life Sciences › Neuroscience › Cellular and Molecular Neuroscience › Alzheimer's disease research and treatments
  • Health Sciences › Medicine › Psychiatry and Mental health › Dementia and Cognitive Impairment Research
Keywords
Alzheimer's disease, prefrontal cortex neurons, BACE1, amyloid plaques, neurofibrillary tangles, cognitive decline, cell cycle re-entry, apoptosis, fine mapping
MeSH terms
Aged, Aged, 80 and over, Alzheimer Disease, Enhancer Elements, Genetic, Female, Humans, Male, Middle Aged, Neurons, Cohort Studies, Cell Adhesion Molecules, Up-Regulation, Aspartic Acid Endopeptidases, Prefrontal Cortex, CpG Islands, DNA Methylation, Epigenesis, Genetic, Amyloid Precursor Protein Secretases, Cognitive Dysfunction

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review