Introduction

Major depressive disorder (MDD) severely limits psychosocial functioning, diminishes quality of life, and is a leading cause of disability worldwide1. The 12-month prevalence of MDD is ~6%2 and similar when comparing high-income countries with low-income and middle-income countries, indicating that MDD is neither a simple consequence of modern day lifestyle in developed countries, nor poverty3,4. Furthermore, although social and cultural factors such as socioeconomic status can have a role in major depression, genomic and other underlying biological factors ultimately drive the occurrence of this condition5. Twin studies have provided heritability estimates of the MDD of ~30–40%6. One of the molecular mechanisms through which environmental and genetic factors can modulate a disease outcome is epigenetics, with DNA methylation being one of the most studied modifications of the genome.

Recent epigenome-wide association studies (EWAS) showed an association of whole blood DNA methylation levels with depressive symptoms7,8 as well as MDD9,10,11, but little is known about brain epigenetic markers of depression or MDD. The EWAS of depressive symptoms were both conducted in late middle-aged and elderly people from the general population (mean age 70 years8 and 65 years7), an age group with an increased risk of developing dementia12. However, both studies could not determine whether their findings were confounded by dementia, which is known to be highly associated with late-life depression13. On the other hand, recent EWAS of MDD were performed in younger participants (mean age 42 years)9,10,11 and it is unclear if their findings can be generalized across age groups. Furthermore, most previous EWAS on depression were limited owing to measuring DNA methylation changes in blood7,8,9,10. Two recent studies conducted EWAS of MDD in 206 post mortem brain samples, but the MDD diagnosis was based on information obtained from a family member and there was no information on dementia9,11. Thus, there is need to understand the epigenetic changes in the human brain that are associated with late-life MDD and to determine whether these changes are independent of dementia.

In this study, we investigated associations between both brain tissue-based individual CpGs as well as regions of differential methylation and late-life MDD in 608 participants from the Religious Order Study and Rush Memory and Aging Project (ROS/MAP) cohorts. To reduce the risk of confounding by cognitive status, we excluded participants with a diagnosis of dementia at the time of MDD assessment and adjusted for cognitive status at the last follow-up visit (closest to methylation assessment) in our analyses. Furthermore, we performed a stratified analysis for men and women to investigate the sex-specific methylation patterns of MDD.

Methods

Study design and study population

The study population included deceased subjects from two large, prospectively followed cohorts recruited by investigators at Rush Alzheimer’s Disease Center in Chicago, IL: The Religious Orders Study (ROS) and the Rush Memory and Aging Project (MAP)14,15. Participants provided informed consent, an Anatomic Gift Act for organ donation, and a repository consent to allow their data to be repurposed. Both studies were approved by an Institutional Review Board of Rush University Medical Center. To be included in the present study, participants must have been assessed for MDD and have available genotype data and methylation profiles derived from the dorsolateral prefrontal cortex. Furthermore, we excluded participants with a diagnosis of dementia at the time of MDD assessment (at baseline evaluation). As in previous publications, the ROS and MAP data were analyzed jointly since much of the phenotypic data collected are identical at the item level in both studies and collected by the same investigative team14,16.

DNA methylation

DNA methylation was measured from the dorsolateral prefrontal cortex (dPFC; Broadman area 46) as previously described in 737 ROS/MAP participant samples14. DNA was extracted from cortically dissected sections of dPFC and DNA methylation was measured using the Illumina HumanMethylation450 Beadchip array. Initial data processing, including color channel normalization, and background removal, was performed using the Illumina GenomeStudio software. The raw IDAT files were obtained from Synapse (www.synapse.org; Synapse ID: syn7357283) and the following probes were removed: (1) probes with a detection p value > 0.01 in any sample, (2) probes annotated to the X and Y chromosomes by Illumina, (3) probes that cross-hybridize with other probes owing to sequence similarity, (3) non-CpG site probes, and (4) probes that overlap with common SNPs. After this filtering, the remaining CpG sites were normalized using the BMIQ algorithm in Watermelon R package17, and the ComBat function from the sva R package was used to adjust for batch effects18. After quality control, 408,689 discrete CpG dinucleotides in 608 subjects were used for analysis.

Genotype data

Genotyping data were generated using two microarrays, Affymetrix GeneChip 6.0 (Affymetrix, Inc, Santa Clara, CA, USA), and Illumina HumanOmniExpress (Illumina, Inc, San Diego, CA, USA) as described previously19. Genotyping was imputed to the 1000 Genome Project Phase 3 using the Michigan Imputation Server20, and the following filtering criteria were applied minor allele frequency (MAF) > 5%, Hardy–Weinberg p value > 10−5 and genotype imputation R2 > 0.3. Principal components were estimated using common (MAF > 0.05) unlinked (R2 < 0.1) autosomal markers by EIGENSTRAT21.

Diagnosis of MDD

A clinical diagnosis of current MDD was rendered by an examining clinician. The diagnosis was based on clinical interview using the criteria of the Diagnostic and Statistical Manual of Mental Disorders, 3rd Edition, Revised (DSM-III-R)22. The MDD diagnosis included present versus not present. In this study, we focused on the diagnosis of MDD at the baseline assessment to reduce the risk that our findings are confounded by dementia.

Clinical diagnosis of cognitive status

A clinical diagnosis of dementia status was rendered based on a three-stage process including computer scoring of cognitive tests, clinical judgment by a neuropsychologist, and diagnostic classification by a clinician. All participants undergo a uniform, structured, clinical evaluation including a battery of 21 cognitive tests of which 19 are in common. These tests were scored by computer using a decision tree designed to mimic clinical judgment and a rating of severity of impairment was given for five cognitive domains. A neuropsychologist, blinded to participant demographics, reviews the impairment ratings and other clinical information and renders a clinical judgment regarding the presence of impairment and dementia. A clinician (neurologist, geriatrician, neuropsychologist, or geriatric nurse practitioner) then reviews all available data and examines the participant and renders a final diagnostic classification. Clinical diagnosis of dementia and clinical Alzheimer’s disease (AD) are based on criteria of the joint working group of the National Institute of Neurological and Communicative Disorders and Stroke and the Alzheimer’s Disease and Related Disorders Association (NINCDS/ADRDA). The diagnosis of AD requires evidence of a meaningful decline in cognitive function relative to a previous level of performance with impairment in memory and at least one other area of cognition. Diagnosis of mild cognitive impairment (MCI) is rendered for persons who are judged to have cognitive impairment by the neuropsychologist but are judged to not meet criteria for dementia by the clinician. Persons without dementia or MCI are categorized as having no cognitive impairment.

Statistical analysis

For the brain EWAS of MDD, we ran a multivariate robust linear regression model with empirical Bayes from the R package limma (version 3.40.6)23 using clinical diagnosis of MDD at baseline as the independent variable and each CpG methylation as a dependent variable, adjusting for age at death, sex, post mortem interval (PMI), proportion of neurons, and the first three genetic principal components. The effect estimates from the adjusted models (Δ beta) refer to the difference in mean DNA methylation beta values between groups (with and without MDD). We applied a Bonferroni threshold to correct for multiple testing based on the number of tested CpG sites (threshold: 0.05/408,689 = 1.22 × 10−07). Fine mapping of our epigenome-wide associations was done with coMET24, which is a visualization tool of EWAS results with functional genomic annotations and estimation of co-methylation patterns. We conducted the following sensitivity analyses: (1) We included the cognitive status at the last follow-up visit (closest to methylation assessment) as a covariate to investigate if our findings were confounded by dementia, (2) We confirmed our associations using linear regression with p values obtained from normal theory (lm() function in R) as well as from a permutation test, (3) we corrected the p values for inflation and bias using a Bayesian method for estimation of the empirical null distribution as implemented in the R/Bioconductor package bacon25, and (4) We adjusted our association models for a polygenic risk score (PRS) for MDD (calculated with PRSice26 and UK Biobank summary statistics from27 with a p value < 0.05) to test if our EWAS findings were independent of genetic risk for MDD. Furthermore, we investigated the overlap between CpG sites associated with MDD and those associated with the PRS for MDD.

CpG sites that reach epigenome-wide significance were replicated using the summary statistics from a cell type-specific EWAS of MDD, which is based on methylation enrichment-based sequencing data from three collections of human post mortem brain (n = 206)11. This replication was used to validate our findings and to provide mechanistic insights about the most relevant cell types for our associations.

Differentially methylated regions (DMRs) in MDD were identified using DMRcate, that identifies DMRs from tunable kernel smoothing process of association signals28. Input files were our single-CpG EWAS results on MDD including regression coefficients, standard deviations, and uncorrected p values. DMRs were defined based on the following criteria: (a) a DMR should contain more than one probe; (b) regional information can be combined from probes within 1000 bp; (c) the region showed FDR corrected p value < 0.05.

To identify plausible pathways associated with MDD, we performed an over-representation analysis based on the 1000 CpGs with the lowest p values for the association with MDD. We used the R Bioconductor package missMethyl (version 1.18.0 gometh function), which performs one-sided hypergeometric tests taking into account and correcting for any bias derived from the use of differing numbers of probes per gene interrogated by the array29.

Results

Description of study participants

There were 608 ROS/MAP participants included in this study with an average age at baseline visit of 81 years and an average age of death of 86 years (Table 1). Sixty-four percent of the participants were female. At baseline, 5% of the participants were diagnosed with MDD, which is consistent with the 12-month prevalence rate of MDD in the general population2. Women showed a slightly higher prevalence of MDD than men (5.4% versus 4.1%, difference not significant).

Table 1 Study characteristics.

Differentially methylated CpG sites in brain tissue are associated with late-life depression

Differentially methylated CpG sites in the YOD1 (cg25594636, p value = 2.55 × 10−11; cg03899372, p value = 3.12 × 10−09; cg12796440, p value = 1.51 × 10−08, cg23982678, p value = 7.94 × 10−08), UGT8 (cg18921206, p value = 1.75 × 10−08), FNDC3B (cg20367479, p value = 4.97 × 10−08) and SLIT2 (cg10946669, p value =  8.01 × 10−08) loci were associated with MDD at the epigenome-wide significance level (Bonferroni-adjustment) after adjusting for sex, PMI, proportion of neurons, first three genetic principal components, and age at death (Table 2, Fig. 1a, Table S1). These associations were robust to additional adjustment for dementia diagnosis assessed at the last follow-up visit (Table 2). Overall, four CpG sites in YOD1 were significantly associated with late-life MDD (Table 2) and these were all located in the same CpG island, but only moderately correlated with each other (Fig. 2). This CpG island is located in an exon of YOD1 and in an intron of PFKFB2. The distribution of the DNA methylation beta values of the four most significant CpG sites in the YOD1 locus stratified by MDD diagnosis is shown in Fig. S1 and the distribution of the DNA methylation beta values of the other three significant CpG sites (cg18921206, cg20367479, and cg10946669) is shown in Fig. S2. The significant associations were confirmed in sensitivity analyses using linear regression models, permutation tests (Table S2) as well as correcting p values for potential inflation and bias (Fig. S3). Furthermore, associations were robust to additional adjustment for a PRS for MDD, which shows that our EWAS findings were independent of a genetic risk for MDD (Table S3).

Table 2 Significant associations between DNA methylation and MDD.
Fig. 1: Manhattan and QQ-plots EWAS on clinical diagnosis of MDD.
figure 1

Adjusted for age at death, sex, PMI, neuron proportion and the first three principal components from the genotype data. Bonferroni threshold: 1.22 × 10−07.

Fig. 2: Fine mapping of the association between DNA methylation in YOD1 and MDD.
figure 2

EWAS results of the association between CpG sites and MDD adjusted for age at death, sex, PMI, neuron proportions and the first three principal components from the genotype data. The most significant CpG site (cg25594636) is marked in purple. The three CpG sites marked in yellow belong to a DMR (p value = 5.06 × 1010, Table S3). The y axis indicates the strength of association in terms of negative logarithm of the association P value. Each circle represents a CpG site. Red dashed line within the graph indicates the genome-wide significance threshold (Bonferroni threshold: 1.22 × 10−07). The regulatory information and correlation matrix of other CpG sites in the region with the top hit are shown below the x axis. Color intensity marks the strength of the correlation and color indicates the direction of the correlation.

Only two of the seven CpG sites that were significantly associated with MDD in ROS/MAP were included in the cell type-specific EWAS published in ref. 11 (cg18921206 and cg20367479, Table S4). Of these, only cg20367479 was nominally significant for bulk brain (p value = 0.022) and the effects were not robust across different cell types in the replication cohort11 and not in the same direction as in ROS/MAP.

Analyses of differentially methylated regions

We identified one significant DMR from our EWAS results on late-life MDD that is located in the YOD1/PFKFB2 locus (Fig. 2, Table S5, minimum FDR p value for the region = 5.06 × 10−10), which is not surprising given the differential CpG site analysis. This DMR includes three CpG sites that are located downstream of the most significant CpG site from our EWAS on late-life MDD (cg25594636, Table 2).

Associations are stronger in men than in women

Interestingly, we found more methylation sites associated with MDD in men than in women (Fig. 1, Table 3, Table S1, Fig. S4), although the sample size was much smaller in men (N = 220 men vs N = 388 women). Differentially methylated CpG sites in YOD1 were more strongly associated with late-life MDD in men than in women (e.g., for cg03899372, men: beta = 0.041, p value = 8.80 × 10−09; women: beta = 0.010, p value = 0.0024; p value sex interaction = 4.51 × 10−06; Table 3). Methylation in PRICKLE4 (p value sex interaction = 1.26 × 10−09), GFAP (p value sex interaction = 6.88 × 10−05), RP11-1E3.1 (p value sex interaction = 1.11 × 10−07) and UBB (p value sex interaction = 1.54 × 10−11) was only associated with MDD in men, but not in women or in both men and women (Table 3).

Table 3 Significant associations between DNA methylation and MDD in male.

Pathway analysis

No significantly enriched pathway was found among the 1000 most significant CpG sites from the EWAS of late-life MDD (Table S6). The smallest p value (6 × 10−5) was reached for calmodulin-dependent protein phosphatase activity (GO:0033192). The genes that belong to this pathway are PPM1A (chr12), PPM1F (chr16), PPP3CA (chr3), PPP3CB (chr14), and PPP3CC (chr14) (Table S7).

To identify biological pathways associated with MDD across both sexes, we conducted a pathway analysis of the CpG sites that reached at least nominal significance in both sex-specific EWAS (p value < 0.05). In total, 807 CpG sites reached at least nominal significance in both sex-specific EWAS (Fig. S5), but no significantly enriched pathway was found among these (Table S8).

Differential DNA methylation in association with polygenic burden for MDD

DNA methylation was not associated with the PRS for MDD at the epigenome-wide significance level (Bonferroni-adjustment) after adjusting for sex, PMI, proportion of neurons, first three genetic principal components, age at death and additional unmeasured confounding using the R package “cate” (Fig. S6, Table S9). Furthermore, the PRS for MDD was not associated with late-life MDD in our study (Fig. S7) and there was only minimal overlap between MDD-associated differentially methylated positions and those associated with the PRS for MDD (Fig. S8).

Discussion

In this study, we found epigenome-wide significant associations between brain tissue-based DNA methylation and late-life MDD in >600 participants from the ROS/MAP cohorts. The most significant and robust association was found with altered methylation levels in the YOD1/PFKFB2 loci. This association was not confounded by dementia or a genetic risk for MDD and significant in both the single site and region-based analysis. Interestingly, brain tissue-based methylation levels were stronger associated with late-life MDD in men than in women.

The most significant CpG sites were found in a region covering an exon of YOD1 and an intron of PFKFB2. YOD1 is a highly conserved deubiquitinase similar to yeast OTU130 that is associated with regulation of the endoplasmic reticulum-associated degradation to maintain the proper folded state of proteins31. In addition, YOD1 is a negative regulator of TRAF6/p62-triggered IL-1 signaling32 and IL-1 has an important role in the regulation of inflammatory responses as well as in depression33,34,35. Together, these suggest that YOD1 is associated with depression perhaps via influencing the inflammatory responses. At last, previous studies suggest that YOD1 contributes to pathogenesis of neurodegenerative disease like Huntington disease and Parkinson’s disease30. PFKFB2 has been studied in the context of brain tumors36,37, but there is no evidence for an association with neuropsychological disease. Therefore, we hypothesize that the CpG sites we found to be associated with late-life MDD are most likely linked to YOD1 regulation.

Further associations with MDD were found for CpG sites in the UGT8, FNDC3B, and SLIT2 loci. UGT8 is a known blood biomarker gene for low mood with evidence of differential expression in human post mortem brains from mood disorder subjects38. In addition, lower expression of UGT8 have been shown in brain tissue from subjects with MDD compared with normal controls39. Therefore, our study extends the current literature by highlighting that not only gene expression, but also brain tissue-based methylation in UGT8 is linked to MDD. FNDC3B and SLIT2 have been discussed in association with brain tumors40, but there is no evidence for an association with neuropsychological disease.

We found stronger associations between brain tissue-based methylation levels and late-life MDD in men than in women. This finding is in line with previous studies showing sex-specific differences in serum biomarkers, mRNA expression, and brain activity of MDD cases, demonstrating that sex has an important role in the molecular heterogeneity of MDD41,42,43. Our findings expand the existing literature by adding DNA methylation from brain tissue to the list of biological patterns that differ between women and men with MDD, which may have important implications for diagnosis as well as treatment strategies.

This is to our knowledge, the first brain tissue-based epigenome-wide study of late-life MDD in a community-based study, and the first EWAS of MDD, which incorporates cognitive status at time of MDD diagnosis as well as at time of death. Two previous EWAS investigated the association with depressive symptoms in middle-aged and elderly people using methylation levels from whole blood7,8. Beside the difference in phenotype definition, the biggest difference between these studies and ours is the tissue in which methylation was measured. In line with a previous study comparing signals from blood and brain tissue44, we could not replicate the whole blood methylation signals from ref. 7,8 in our brain tissue-based EWAS (Tables S10 and S11). In two recent brain tissue-based EWAS of MDD, differential methylation was measured in 206 post mortem brain samples by enrichment-based sequencing9,11. However, as in ROS/MAP methylation was measured with the Illumina HumanMethylation450 Beadchip array, loci overlapping between blood and brain in9 (chr2: 208,230,169; chr9: 101,119,679; chr4: 71,632,888) were not available in ROS/MAP; therefore, we could not use the data set from Aberg et al. for replication purposes. Using the same samples, Chan et al.11 conducted a cell type-specific EWAS, in which none of the CpGs reached epigenome-wide significance for neurons and glia or for bulk brain. Owing to the different assessment of methylation (array-based versus sequencing), only two of our seven significant CpG sites were available in ref. 11 and both of them were not successfully replicated.

Strengths of this study include the ROS/MAP cohort itself, which is notable for its longitudinal nature with very high follow-up rates, prospective collection of data, a community-based cohort design, and high autopsy rates. Furthermore, the 12-month prevalence of MDD in our study population matches that in the general population2, which makes our findings generalizable beyond our study population. Another strength of our study is the analysis of methylation levels from brain tissue, which is the most relevant tissue for the pathophysiology of depression. In addition, we reduced the risk of confounding by cognitive status by excluding participants with a diagnosis of dementia at time of MDD assessment and by adjusting our analyses for cognitive status at the last follow-up visit (closest to methylation assessment).

The study is potentially limited by the use of bulk tissue analysis, which might obscure signals from different cell populations. This problem was mitigated in our analysis by adjusting for cell type composition. Future studies should investigate the role of YOD1, UGT8, FNDC3B, and SLIT2 in specific cell types from brain and investigate whether there is a causal relationship between gene dysregulation and MDD in animal models. Up to now, there is only one study analyzing cell type-specific associations between DNA methylation and MDD. However, the authors did not find any significant associations with MDD and owing to the different approach of assessing differential DNA methylation (array-based versus sequencing), our most significant CpG sites were not available in their data11. Therefore, there is an urgent need for more large-scale brain tissue-based EWAS of MDD to validate our and the previous9,11 findings and to better understand the consequences of MDD on the human brain. Another limitation of our EWAS was the small number of MDD cases in our study population. However, to reduce the risk of false positive findings owing to the imbalanced study design, we validated our findings by using different modelling approaches (limma, linear regression, permutation tests, DMR analysis).

In conclusion, we have presented evidence for brain-based DNA methylation in association with late-life MDD. We identified methylation in YOD1, UGT8, FNDC3B, and SLIT2 as new epigenetic factors associated with late-life MDD, which are not confounded by cognitive status or a genetic risk for MDD and stronger associated with MDD in male than in female.