Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2025 Feb 24;16(1):1908.
doi: 10.1038/s41467-024-55356-x.

Insights from a methylome-wide association study of antidepressant exposure

Affiliations

Insights from a methylome-wide association study of antidepressant exposure

E Davyson et al. Nat Commun. .

Abstract

This study tests the association of whole-blood DNA methylation and antidepressant exposure in 16,531 individuals from Generation Scotland (GS), using self-report and prescription-derived measures. We identify 8 associations and a high concordance of results between self-report and prescription-derived measures. Sex-stratified analyses observe nominally significant increased effect estimates in females for four CpGs. There is observed enrichment for genes expressed in the Amygdala and annotated to synaptic vesicle membrane ontology. Two CpGs (cg15071067; DGUOK-AS1 and cg26277237; KANK1) show correlation between DNA methylation with the time in treatment. There is a significant overlap in the top 1% of CpGs with another independent methylome-wide association study of antidepressant exposure. Finally, a methylation profile score trained on this sample shows a significant association with antidepressant exposure in a meta-analysis of eight independent external datasets. In this large investigation of antidepressant exposure and DNA methylation, we demonstrate robust associations which warrant further investigation to inform on the design of more effective and tolerated treatments for depression.

PubMed Disclaimer

Conflict of interest statement

Competing interests: R.E.M. is an advisor to the Epigenetic Clock Development Foundation and Optima Partners. D.L.M. was employed by Optima Partners Ltd in a part-time capacity. H.J.G. has received travel grants and speaker honoraria from Fresenius Medical Care, Neuraxpharm, Servier, Indorsia and Janssen Cilag, as well as research funding from Fresenius Medical Care. H.J.G. had personal contracts approved by the university administration for speaker honoraria and one IIT with Fresenius Medical Care. T.K. received unrestricted educational grants from Servier, Janssen, Recordati, Aristo, Otsuka, neuraxpharm. B.W.J.H.P. has received research funding (not related to the current paper) from Boehringer Ingelheim and Jansen Research. All other authors report no biomedical financial interests or potential conflicts of interest.

Figures

Fig. 1
Fig. 1. Study design of investigating antidepressant exposure and DNA methylation.
Data: Participants in Generation Scotland (GS) provided blood samples from which DNA methylation was measured. Their antidepressant exposure status was measured using both self-report questionnaires and prescription-derived measures. Prescription-derived measures: Repeated regular prescriptions over time (X axis) for antidepressants (purple bars) are merged to form active antidepressant treatment periods (blue bars). Individuals in an active treatment period at the time of blood sample (black cross) are classed as antidepressant exposed. Methylome-wide association studies: A methylome-wide association study (MWAS), subsequent regional analysis and functional annotation was performed for both measures of antidepressant exposure. Additionally, an enrichment analysis was done using MBD-Sequencing data in the Netherlands Depression and Anxiety (NESDA) cohort for the self-report antidepressant exposure. Methylation Profile Score: Weights for a methylation profile score (MPS) of self-reported antidepressant exposure was calculated in GS using a least absolute shrinkage and selection operator (LASSO) model. Eight independent datasets then tested the association of this MPS with self-reported antidepressant exposure. Created in BioRender. Davyson, E. (2024) https://BioRender.com/s43g100.
Fig. 2
Fig. 2. Methylome-wide association studies (MWAS) of self-reported and prescription-derived antidepressant exposure.
Manhattan plots of the MWAS of self-report (A) and prescription-derived (B) antidepressant exposure using a Mixed-linear-model Omics-based Analysis (MOA) model. Significance was assessed using the p-value threshold 9.42 × 10−8, as recommended for case-control MWAS analyses. C The MOA MWAS standardised effect sizes and 95% confidence-intervals (effect estimate +/- 1.96*Standard Error) for associated CpGs (p < 9.42 × 108) for the full sample (Nself-report = 16,531, Nprescription-derived = 7951) and MDD-subgroup sample (Nself-report = 2268, Nprescription-derived = 792). Data are presented as the MWAS effect estimates +/- the 95% confidence intervals. Effect sizes represent a per-1 standard-deviation increase in CpG methylation M-values.
Fig. 3
Fig. 3. An antidepressant exposure methylation profile score (MPS) and antidepressant exposure in external cohorts.
A The number of participants exposed and unexposed to antidepressants in each cohort (NNTR = 3087, NSTRADL = 658, NFTC = 1678, NSHIP-TREND = 495, NALSPAC = 801, NLBC1936 = 889, NMARS-UniDep = 312, NERISK = 1658, NFOR2107 = 658). B Nagelkerke’s pseudo R2, representing an estimate of how much variance in the antidepressant exposure outcome that is explained by the MPS in each cohort. C The effect estimates and 95% confidence intervals (effect estimate +/− 1.96*Standard Error) of MPS ~ antidepressant exposure association in each cohort, using either a generalised linear model (FOR2107 and ALSPAC), a generalised linear mixed model (SHIP-Trend, LBC1936, MARS-UniDep, STRADL and E-Risk) or a generalised estimation equation model (FTC and NTR). The square size = study weight in the random-effects meta-analysis. The pooled effect estimate interval, calculated from a random-effects meta-analysis, is represented by the blue diamond (Npooled=10,236). Created in BioRender. Davyson, E. (2024) https://BioRender.com/q70l792.

Update of

References

    1. GBD 2017 Dis Injury Incidence Pr, James, S. L. G. et al. Global, regional, and national incidence, prevalence, and years lived with disability for 354 diseases and injuries for 195 countries and territories, 1990-2017: a systematic analysis for the Global Burden of Disease Study 2017. The Lancet. 2018;392:1789-1858. 10.1016/S0140-6736(18)32279-7. - PMC - PubMed
    1. Braund, T. A. et al. Antidepressant side effects and their impact on treatment outcome in people with major depressive disorder: an iSPOT-D report. Transl. Psychiatry11, 1–8 (2021). - PMC - PubMed
    1. Cipriani, A. et al. Comparative efficacy and acceptability of 21 antidepressant drugs for the acute treatment of adults with major depressive disorder: a systematic review and network meta-analysis. Lancet391, 1357–1366 (2018). - PMC - PubMed
    1. Souery, D. et al. Switching antidepressant class does not improve response or remission in treatment-resistant depression. J. Clin. Psychopharmacol.31, 512–516 (2011). - PubMed
    1. McIntyre, R. S. et al. Treatment-resistant depression: definition, prevalence, detection, management, and investigational interventions. World Psychiatry22, 394–412 (2023). - PMC - PubMed

Substances

Grants and funding

LinkOut - more resources