Participant Characteristics
This study utilized samples from 404 participants from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) study with DNA methylation and structural neuroimaging data available. At baseline, 121 participants were considered CN, 236 were diagnosed with MCI, and 47 were diagnosed with AD. ADNI is a longitudinal study of subjects across the US that includes multimodal neuroimaging, blood biomarkers, and clinical markers of AD and has been described previously10,11. As part of the study, patients undergo a rigorous clinical exam, which includes neurologic examination, neuropsychiatric evaluation, cognitive testing, blood sampling, and brain MRI12. Participants were diagnosed as CN or with MCI or AD based on a structured protocol that integrated clinical data, cognitive testing, and brain MRI results, which has been described previously13. Clinical characteristics are detailed in Table 1. Clinical symptoms severity was measured using the Clinical Dementia Rating Scale Sum of Boxes (CDR-SB) score14. Cognitive changes over time were measured using the Montreal Cognitive Assessment (MoCA)15. Written and informed consent was obtained from study participants in compliance with local IRB and ADNI protocols as has been previously described16.
Image Processing
Structural MRI images from all participants were segmented using two pipelines to evaluate cortical structure and white matter disease. Cortical segmentation and structural analysis was performed using FreeSurfer version 5.1 using T1-weighted images as previously described17,18. Briefly, all images were segmented using FreeSurfer’s automated pipeline and then manually checked for segmentation accuracy, with segmentation errors manually corrected19. Cortical thickness measurements from 68 cortical regions of interest (ROI; 34 for each hemisphere) in the Desikan-Killiany atlas18 were retained for analysis. White matter disease was estimated using a pipeline developed at the University of California, Davis, which quantifies WMH volumes via Bayesian modeling based on participants’ 3D T1 and FLAIR images20.
DNA Methylation Assessment
DNA methylation was profiled from blood samples of ADNI participants using the Illumina Infinium Human MethylationEPIC V1 BeadChip Array, which covers ~ 866,000 CpGs (illumina.com). Briefly, data was normalized using the dasen method in the wateRmelon R package21, and quality control procedures were performed, including removal of samples with abnormal CpG detection p-values > 0.05, checking the ratio of X/Y chromosome probe intensities for sex concordance, and comparison of targeted SNP genotypes to genotype microarray data as previously described22. Estimates of epigenetic age for the Levine 513 CpG site DNAmPhenoAge clock8 and DNA methylation-based mortality risk assessment (DNAmGrimAge)7 were calculated using R scripts. Mean imputation was utilized for missing values. The first time point for each peripheral blood sample was used for epigenetic age estimation. For a subset of patients, technical replicates were present, and in these situations, we averaged the mean epigenetic age across replicates for all downstream analyses.
Statistical Analyses
All statistical analyses were completed using R 4.1.223. Longitudinal analyses of the effect of epigenetic age on disease progression were performed using Cox proportional hazard modeling with the R package ‘survival’24, controlling for APOE ε4 allele dosage, sex, years of education, and CDR-SB. Conversion to disease was analyzed in CN participants and was defined as a change in clinical diagnosis from CN to MCI or CN to AD. The proportional hazards assumption was tested for each model using the ‘cox.zph’ function and was not statistically significant (p > 0.05 level for all analyses shown). Mixed-effects linear regression analyses were used to assess the relationship between baseline epigenetic age and longitudinal MoCA scores, controlling for baseline and time interactions of baseline MoCA score, sex, years of education, and APOE 𝜀4 dosage. We used the following linear mixed-effect model:
ΔMoCA = 𝛽0 + 𝛽1Δ𝑡 + 𝛽2MoCAbaseline*Δ𝑡 + 𝛽3Sex*Δ𝑡 + 𝛽4Education*Δ𝑡 + 𝛽5𝐴𝑃𝑂𝐸𝜀4*Δ𝑡 + 𝑒
Associations between epigenetic age and cortical thickness were conducted using multiple regression after covarying for APOE ε4 dosage, sex, years of education, CDR-SB, and total intracranial volume. Volume-rendered images of multiple regression analysis results were created using the R package ‘fsbrain’25.
As has been done in prior work26 and in line with our hypotheses, we first performed all analyses using epigenetic age alone to predict disease progression and neuroimaging biomarkers. For outcomes which demonstrated a significant relationship at p < 0.05, we specifically tested whether epigenetic age provided independent information above and beyond chronologic age by covarying for chronologic age in addition to all previously listed covariates.
Where applicable, correction for multiple comparisons was completed using the false discovery rate (FDR) technique27.
Approval and Consent
All ADNI research activities were approved by Institutional Review Boards at each patient recruitment site. Written informed consent was obtained from all patients.