Hostname: page-component-74656f9884-kv2lp Total loading time: 0 Render date: 2026-11-05T20:00:46.696Z Has data issue: false hasContentIssue false
  • English
  • Français

Epigenetic aging and PTSD outcomes in the immediate aftermath of trauma

Published online by Cambridge University Press:  23 March 2023

Anthony S. Zannas Open the ORCID record for Anthony S. Zannas [Opens in a new window] ,
Sarah D. Linnstaedt ,
Xinming An ,
Jennifer S. Stevens ,
Nathaniel G. Harnett ,
Alyssa R. Roeckner ,
Katelyn I. Oliver ,
David R. Rubinow ,
Elisabeth B. Binder  and
Karestan C. Koenen
...Show all authors
Anthony S. Zannas*
Affiliation:
Department of Psychiatry, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA Carolina Stress Initiative, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA Institute for Trauma Recovery, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA
Sarah D. Linnstaedt
Affiliation:
Institute for Trauma Recovery, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA Department of Anesthesiology, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA
Xinming An
Affiliation:
Institute for Trauma Recovery, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA Department of Anesthesiology, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA
Jennifer S. Stevens
Affiliation:
Department of Psychiatry and Behavioral Sciences, Emory University School of Medicine, Atlanta, GA, USA
Nathaniel G. Harnett
Affiliation:
Department of Psychiatry, McLean Hospital, Harvard Medical School, Boston, MA, USA
Alyssa R. Roeckner
Affiliation:
Department of Psychiatry and Behavioral Sciences, Emory University School of Medicine, Atlanta, GA, USA
Katelyn I. Oliver
Affiliation:
Department of Psychiatry and Behavioral Sciences, Emory University School of Medicine, Atlanta, GA, USA
David R. Rubinow
Affiliation:
Department of Psychiatry, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA
Elisabeth B. Binder
Affiliation:
Department of Translational Research in Psychiatry, Max Planck Institute of Psychiatry, Munich, Germany
Karestan C. Koenen
Affiliation:
Department of Epidemiology, Harvard School of Public Health, Harvard University, Boston, MA, USA Department of Social and Behavioral Sciences, Harvard School of Public Health, Harvard University, Boston, MA, USA
Kerry J. Ressler
Affiliation:
Department of Psychiatry, McLean Hospital, Harvard Medical School, Boston, MA, USA
Samuel A. McLean
Affiliation:
Department of Psychiatry, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA Institute for Trauma Recovery, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA Department of Anesthesiology, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA
*
Author for correspondence: Anthony S. Zannas, E-mail: anthony_zannas@med.unc.edu
Rights & Permissions [Opens in a new window]

Abstract

Background

Psychological trauma exposure and posttraumatic stress disorder (PTSD) have been associated with advanced epigenetic age. However, whether epigenetic aging measured at the time of trauma predicts the subsequent development of PTSD outcomes is unknown. Moreover, the neural substrates underlying posttraumatic outcomes associated with epigenetic aging are unclear.

Methods

We examined a multi-ancestry cohort of women and men (n = 289) who presented to the emergency department (ED) after trauma. Blood DNA was collected at ED presentation, and EPIC DNA methylation arrays were used to assess four widely used metrics of epigenetic aging (HorvathAge, HannumAge, PhenoAge, and GrimAge). PTSD symptoms were evaluated longitudinally at the time of ED presentation and over the ensuing 6 months. Structural and functional neuroimaging was performed 2 weeks after trauma.

Results

After covariate adjustment and correction for multiple comparisons, advanced ED GrimAge predicted increased risk for 6-month probable PTSD diagnosis. Secondary analyses suggested that the prediction of PTSD by GrimAge was driven by worse trajectories for intrusive memories and nightmares. Advanced ED GrimAge was also associated with reduced volume of the whole amygdala and specific amygdala subregions, including the cortico-amygdaloid transition and the cortical and accessory basal nuclei.

Conclusions

Our findings shed new light on the relation between biological aging and trauma-related phenotypes, suggesting that GrimAge measured at the time of trauma predicts PTSD trajectories and is associated with relevant brain alterations. Furthering these findings has the potential to enhance early prevention and treatment of posttraumatic psychiatric sequelae.

Keywords

AgingDNA methylationepigeneticsPTSDstress

Information

Type
Original Article
Information
Psychological Medicine , Volume 53 , Issue 15 , November 2023 , pp. 7170 - 7179
Check for updates
Copyright
Copyright © The Author(s), 2023. Published by Cambridge University Press

Introduction

Psychological trauma and stress-related phenotypes have long been linked with accelerated aging. This link has repeatedly captured the imagination of literary writers (Hugo, Reference Hugo2012; Zannas, Reference Zannas2019b), has been observed in clinical settings (Bersani, Mellon, Reus, & Wolkowitz, Reference Bersani, Mellon, Reus and Wolkowitz2019), and has been supported by epidemiological studies (Felitti et al., Reference Felitti, Anda, Nordenberg, Williamson, Spitz, Edwards and Marks1998; Vaccarino et al., Reference Vaccarino, Goldberg, Rooks, Shah, Veledar, Faber and Bremner2013). Dissecting the underlying mechanisms is important and timely, given the high prevalence of trauma exposure and global aging of the human population (Breslau et al., Reference Breslau, Kessler, Chilcoat, Schultz, Davis and Andreski1998; U.S. Department of Health and Human Services, 2013). Among plausible mechanisms, epigenetics – the chemical changes that regulate genomic function without altering the genetic code – has emerged as a key link between stress exposure and health outcomes and as a molecular hallmark of the aging process (Cavalli & Heard, Reference Cavalli and Heard2019; Gassen, Chrousos, Binder, & Zannas, Reference Gassen, Chrousos, Binder and Zannas2016; Lopez-Otin, Blasco, Partridge, Serrano, & Kroemer, Reference Lopez-Otin, Blasco, Partridge, Serrano and Kroemer2013). In particular, a critical epigenetic modification is DNA methylation (DNAm) in the cytosine-guanine (CpG) context, which through array technology has become widely studied in humans (Yong, Hsu, & Chen, Reference Yong, Hsu and Chen2016; Zannas, Reference Zannas2019a).

DNAm patterns have been shown to change extensively with age (Fraga et al., Reference Fraga, Ballestar, Paz, Ropero, Setien, Ballestar and Esteller2005; Horvath & Raj, Reference Horvath and Raj2018), and composite (multi-CpG) methylomic markers (so-called ‘epigenetic aging’) that combine the DNAm status of multiple age-regulated CpG sites can predict not only chronological age (Hannum et al., Reference Hannum, Guinney, Zhao, Zhang, Hughes, Sadda and Zhang2013; Horvath, Reference Horvath2013) but also diverse health outcomes (Hillary et al., Reference Hillary, Stevenson, McCartney, Campbell, Walker, Howard and Marioni2020; Joyce et al., Reference Joyce, Gao, Zheng, Ma, Hwang, Liu and Lloyd-Jones2021; Levine et al., Reference Levine, Lu, Quach, Chen, Assimes, Bandinelli and Horvath2018; Lu et al., Reference Lu, Quach, Wilson, Reiner, Aviv, Raj and Horvath2019; McCrory et al., Reference McCrory, Fiorito, Hernandez, Polidoro, O'Halloran, Hever and Kenny2021; Zheng et al., Reference Zheng, Joyce, Colicino, Liu, Zhang, Dai and Hou2016) [reviewed in (Horvath & Raj, Reference Horvath and Raj2018; Palma-Gudiel, Fañanás, Horvath, & Zannas, Reference Palma-Gudiel, Fañanás, Horvath and Zannas2020)]. The early (first-generation) metrics by Horvath and Hannum et al. were derived by regression models aiming to predict chronological age, and the difference between DNAm-predicted and chronological age was proposed as a measure of an individual's biological aging (Hannum et al., Reference Hannum, Guinney, Zhao, Zhang, Hughes, Sadda and Zhang2013; Horvath, Reference Horvath2013). Subsequently developed (second-generation) metrics, such as the widely used PhenoAge and GrimAge, further aimed to predict healthspan and lifespan by including in their regression model clinical biomarkers and mortality endpoints (Levine et al., Reference Levine, Lu, Quach, Chen, Assimes, Bandinelli and Horvath2018; Lu et al., Reference Lu, Quach, Wilson, Reiner, Aviv, Raj and Horvath2019).

Leveraging these markers, we and others previously linked various types of stress and trauma exposure with advanced epigenetic age (Belsky et al., Reference Belsky, Caspi, Corcoran, Sugden, Poulton, Arseneault and Moffitt2022; Boks et al., Reference Boks, van Mierlo, Rutten, Radstake, De Witte, Geuze and Vermetten2015; Brody, Yu, Chen, Beach, & Miller, Reference Brody, Yu, Chen, Beach and Miller2016; Copeland, Shanahan, McGinnis, Aberg, & van den Oord, Reference Copeland, Shanahan, McGinnis, Aberg and van den Oord2022; Harvanek, Fogelman, Xu, & Sinha, Reference Harvanek, Fogelman, Xu and Sinha2021; Katrinli et al., Reference Katrinli, Stevens, Wani, Lori, Kilaru, van Rooij and Smith2020; Lim, Nzegwu, & Wright, Reference Lim, Nzegwu and Wright2022; Wolf et al., Reference Wolf, Maniates, Nugent, Maihofer, Armstrong, Ratanatharathorn and Logue2018; Zannas et al., Reference Zannas, Arloth, Carrillo-Roa, Iurato, Roh, Ressler and Mehta2015a). Moreover, published work to date has associated advanced epigenetic age with posttraumatic stress disorder (PTSD) (Katrinli et al., Reference Katrinli, Stevens, Wani, Lori, Kilaru, van Rooij and Smith2020; Kuan et al., Reference Kuan, Ren, Clouston, Yang, Jonas, Kotov and Luft2021; Mehta et al., Reference Mehta, Bruenig, Pierce, Sathyanarayanan, Stringfellow, Miller and Shakespeare-Finch2022; Na et al., Reference Na, Montalvo-Ortiz, Nagamatsu, Southwick, Krystal, Gelernter and Pietrzak2022; Wang et al., Reference Wang, Hui, Goldberg, Smith, Kaseer, Murrah and Sun2022; Wolf et al., Reference Wolf, Maniates, Nugent, Maihofer, Armstrong, Ratanatharathorn and Logue2018), though a lack of and even an opposite direction of association have been reported (Boks et al., Reference Boks, van Mierlo, Rutten, Radstake, De Witte, Geuze and Vermetten2015; Mehta et al., Reference Mehta, Bruenig, Lawford, Harvey, Carrillo-Roa, Morris and Voisey2018; Verhoeven et al., Reference Verhoeven, Yang, Wolkowitz, Bersani, Lindqvist, Mellon and Hammamieh2018). Among studies reporting positive associations, findings further vary depending on the timing of PTSD diagnosis: several cohorts have associated advanced epigenetic age with either lifetime or current PTSD (Katrinli et al., Reference Katrinli, Stevens, Wani, Lori, Kilaru, van Rooij and Smith2020; Kuan et al., Reference Kuan, Ren, Clouston, Yang, Jonas, Kotov and Luft2021; Mehta et al., Reference Mehta, Bruenig, Pierce, Sathyanarayanan, Stringfellow, Miller and Shakespeare-Finch2022; Na et al., Reference Na, Montalvo-Ortiz, Nagamatsu, Southwick, Krystal, Gelernter and Pietrzak2022; Wang et al., Reference Wang, Hui, Goldberg, Smith, Kaseer, Murrah and Sun2022), but a meta-analysis found this association to be significant for lifetime PTSD only (Wolf et al., Reference Wolf, Maniates, Nugent, Maihofer, Armstrong, Ratanatharathorn and Logue2018). Such variable findings suggest that epigenetic aging and PTSD risk are linked through a complex relationship, the direction of which remains unclear. Importantly, no studies to date have examined whether epigenetic aging measured at the time of trauma predicts the subsequent development of PTSD outcomes. This hypothesis is plausible, given that several hallmarks of aging, including chronic inflammation, metabolic dysregulation, stem cell dysfunction, and epigenetic alterations, are thought to play key roles in PTSD pathogenesis (Kao et al., Reference Kao, He, Zannas, Hahn, Kühne, Reichel and Turck2016; López-Otín, Blasco, Partridge, Serrano, & Kroemer, Reference López-Otín, Blasco, Partridge, Serrano and Kroemer2022; Mellon, Gautam, Hammamieh, Jett, & Wolkowitz, Reference Mellon, Gautam, Hammamieh, Jett and Wolkowitz2018; Seo et al., Reference Seo, Park, Park, Kim, Kim and Kim2019; Zannas, Provençal, & Binder, Reference Zannas, Provençal and Binder2015b). Moreover, there is a paucity of studies integrating epigenomic and phenotypic assessments with neuroimaging measures to uncover the neural substrates underlying posttraumatic outcomes associated with epigenetic aging.

To address these knowledge gaps, the present study leverages the AURORA (Advancing Understanding of RecOvery afteR traumA) cohort (McLean et al., Reference McLean, Ressler, Koenen, Neylan, Germine, Jovanovic and Kessler2020), a multi-ancestry cohort of women and men who presented to the emergency department (ED) after trauma. Participant assessments included blood collection at ED presentation, longitudinal PTSD symptoms during the 6 months following trauma exposure, and structural and functional neuroimaging 2 weeks after trauma. To capture the potentially different aspects and unique contributions of epigenetic aging markers, we here examine all four aforementioned, widely used epigenetic aging markers (HorvathAge, HannumAge, PhenoAge, and GrimAge) (Hannum et al., Reference Hannum, Guinney, Zhao, Zhang, Hughes, Sadda and Zhang2013; Horvath, Reference Horvath2013; Levine et al., Reference Levine, Lu, Quach, Chen, Assimes, Bandinelli and Horvath2018; Lu et al., Reference Lu, Quach, Wilson, Reiner, Aviv, Raj and Horvath2019). Given that prevention and treatment of psychiatric sequelae would greatly benefit from biomarkers available early after trauma exposure, we first examine if epigenetic aging at ED presentation predicts the development of PTSD outcomes during follow-up. We then assess structural and functional neural correlates of epigenetic aging that may be relevant for PTSD outcomes. In particular, our structural MRI analyses focus on the amygdala and the hippocampus, two brain regions with established roles in stress and trauma-related phenotypes (Del Casale et al., Reference Del Casale, Ferracuti, Barbetti, Bargagna, Zega, Iannuccelli and Pompili2022; Morey et al., Reference Morey, Gold, LaBar, Beall, Brown, Haswell and McCarthy2012), whereas functional MRI analyses explore alterations in network connectivity, which have been linked with PTSD outcomes (Korgaonkar et al., Reference Korgaonkar, Chakouch, Breukelaar, Erlinger, Felmingham, Forbes and Bryant2020; Sheynin et al., Reference Sheynin, Duval, King, Angstadt, Phan, Simon and Liberzon2020).

Methods

Study participants

All data for the present report are obtained from the AURORA (McLean et al., Reference McLean, Ressler, Koenen, Neylan, Germine, Jovanovic and Kessler2020), a large multi-ancestry cohort study (total n > 3000) that involves women and men presenting to the ED within 72 h after exposure to psychological trauma. Inclusion and exclusion criteria for AURORA participants were as follows. Patients aged 18–75 years who presented to the ED within 72 h of trauma exposure at participating ED sites were screened for study eligibility. Trauma exposures automatically qualifying for study enrollment were motor vehicle collision, physical assault, sexual assault, fall greater than 10 feet, or mass casualty incidents. Other trauma exposures also qualified if: (1) the individual responded to a screener question that they experienced the exposure as involving actual or threatened serious injury, sexual violence, or death, either by direct exposure, witnessing, or learning about it; and (2) the research assistant agreed that the exposure was a plausible qualifying event. Exclusion criteria included administration of general anesthesia, long bone fractures, laceration with significant hemorrhage, solid organ injury > American Association for the Surgery of Trauma Grade 1, not alert and oriented at the time of enrollment, not fluent in written or spoken English, visual or auditory impairment precluding completion of web-based neurocognitive evaluations and/or telephone follow-ups, self-inflicted or occupational injury, prisoners, individuals pregnant or breastfeeding, individuals reporting ongoing domestic violence, and individuals taking >20 mg morphine or equivalent per day. To be eligible for the study, patients also needed to have an iOS or Android-compatible smartphone with internet access and an email address that they check regularly.

The present study is focused on a subset of AURORA participants (n = 289) in whom epigenetic assessments of ED blood samples were performed. A smaller subset of these individuals (n = 63) also underwent neuroimaging assessments 2 weeks after trauma. Clinicodemographic characteristics of included participants are presented in Table 1.

Table 1.

Participant demographic and clinical variables

MRI, magnetic resonance imaging; n, number; r, Pearson correlation coefficient; s.d., standard deviation.

PTSD at 6 months was defined using the PTSD Checklist for DSM-5 (PCL-5) and a previously established score threshold of ⩾31 (Blevins et al., Reference Blevins, Weathers, Davis, Witte and Domino2015; Bovin et al., Reference Bovin, Marx, Weathers, Gallagher, Rodriguez, Schnurr and Keane2016; Kessler et al., Reference Kessler, Ressler, House, Beaudoin, An, Stevens and McLean2021).

Phenotypic measures

Probable PTSD diagnosis at 6 months was defined using the PTSD Checklist for DSM-5 (PCL-5) – a 20-item self-report scale that uses a 0–4 response format asking how much the participant was ‘bothered by’ each PTSD symptom (0–4 scale) in the past 30 days – and a previously established PCL-5 score threshold of ⩾31 (Blevins, Weathers, Davis, Witte, & Domino, Reference Blevins, Weathers, Davis, Witte and Domino2015; Bovin et al., Reference Bovin, Marx, Weathers, Gallagher, Rodriguez, Schnurr and Keane2016; Kessler et al., Reference Kessler, Ressler, House, Beaudoin, An, Stevens and McLean2021). Distinct PTSD-related symptom trajectories (intrusive memories, hyperarousal, avoidance, nightmares, and sleep disturbance) were characterized using data from symptom assessments collected from AURORA participants at 10 timepoints during the 6 months after trauma exposure (McLean et al., Reference McLean, Ressler, Koenen, Neylan, Germine, Jovanovic and Kessler2020). Each of these symptoms is assessed through smartphone-based surveys using 2–3 items asking participants to rate the symptom severity or frequency (0–4 scale) experienced over the preceding days. For each symptom, latent trajectory classes were defined using growth mixture models. Given the smaller cohort subset with epigenetic data, statistical power was increased here by comparing the combined trajectories with low or moderate recovery symptoms v. those with moderate or high persistent symptoms over the follow-up duration. Further details are provided in online Supplementary Methods and in Beaudoin et al. (Reference Beaudoin, An, Basu, Ji, Liu, Kessler and McLean2023). PTSD symptoms at ED presentation and over the lifetime were assessed using the abbreviated (six-item) civilian version of PCL-5 (PCL-C), and the presence of significant PTSD symptoms at ED presentation was defined using a previously established score threshold of ⩾14 (Lang & Stein, Reference Lang and Stein2005). Childhood trauma history was assessed using a modified, 11-item survey derived from the short version of the childhood trauma questionnaire (CTQ) (Bernstein et al., Reference Bernstein, Stein, Newcomb, Walker, Pogge, Ahluvalia and Zule2003). The survey included two items each from the physical neglect, emotional neglect, emotional abuse, and physical abuse subtype and three items from the sexual abuse subtype, with each item asking the frequency (0–4 scale) of traumatic experience during the participant's childhood. Lifetime trauma burden was defined as the sum of all questionnaire items comprising the Life Events Checklist (LEC-5) (Gray, Litz, Hsu, & Lombardo, Reference Gray, Litz, Hsu and Lombardo2004; Weathers et al., Reference Weathers, Blake, Schnurr, Kaloupek, Marx and Keane2013). General mental and physical health at the time of ED presentation were derived using normative scores based on questions from the 12-item Short Form Health Survey (SF-12) (Ware, Kosinski, & Keller, Reference Ware, Kosinski and Keller1996). Quantity of daily tobacco and alcohol use at ED presentation was assessed using the PhenX toolkit (Hamilton et al., Reference Hamilton, Strader, Pratt, Maiese, Hendershot, Kwok and Haines2011).

DNA methylation

Blood samples were collected in the ED using DNA PAXgene tubes (Qiagen, Germantown MD, USA), frozen at −20°C at each collection study site, and then batch-shipped on dry ice to the NIMH repository (Piscataway, NJ, USA). DNA was isolated upon arrival to the repository using chemagen magnetic bead technology via Chemagic 360 instrumentation (PerkinElmer, Waltham, MA, USA). DNA concentration and purity were determined using UV/Vis on a Lunatic reader (Unchained Labs, Pleasanton, CA, USA). Bisulfite conversion of the isolated DNA was performed at the University of Minnesota Genomics Center, St. Paul, MN using EZ-96 DNA Methylation Kits (Zymo Research, Irvine, CA, USA) according to the manufacturer's instructions. DNAm was quantified using the Infinium Human MethylationEPIC BeadChip (Illumina Inc., San Diego, CA, USA). All steps were performed manually except for hybridization and staining steps, which were performed by a liquid handling robot (Tecan, Männedorf, Switzerland). To account for potential technical batch effects, DNA samples from different outcomes were randomized across beadchips. Quality control (QC) was performed using the CHAMP package in RStudio (Tian et al., Reference Tian, Morris, Webster, Yang, Beck, Feber and Teschendorff2017). ChAMP is an integrated analysis pipeline that filters low-quality probes and samples, adjusts for Infinium I and Infinium II probe design, and corrects for batch effects. Methylation data were cleaned by removing: (i) probes with low detection p > 0.1 or with beadcount <3 in at least 5% of samples; (ii) previously identified cross-reactive and polymorphic probes; (iii) probes containing SNPs that overlap with a CpG site, at single base extension sites, or when the CpG probe was located near short insertions or deletions; and (iv) probes located on the X and Y chromosomes. Data were then visually inspected by singular value decomposition and remaining batch effects were removed using ComBat (Johnson, Li, & Rabinovic, Reference Johnson, Li and Rabinovic2007; Leek, Johnson, Parker, Jaffe, & Storey, Reference Leek, Johnson, Parker, Jaffe and Storey2012). QC-processed DNAm data (β values) were used to calculate epigenetic aging with the online Horvath calculator, and downstream analyses focused on the four widely used first- and second-generation epigenetic aging markers: HorvathAge (Horvath, Reference Horvath2013), HannumAge (Hannum et al., Reference Hannum, Guinney, Zhao, Zhang, Hughes, Sadda and Zhang2013), PhenoAge (Levine et al., Reference Levine, Lu, Quach, Chen, Assimes, Bandinelli and Horvath2018), and GrimAge (Lu et al., Reference Lu, Quach, Wilson, Reiner, Aviv, Raj and Horvath2019). Since different immune cell types have distinct epigenetic profiles and their blood distribution can be influenced by stress (Adalsteinsson et al., Reference Adalsteinsson, Gudnason, Aspelund, Harris, Launer, Eiriksdottir and Gudnason2012; Beis et al., Reference Beis, von Känel, Heimgartner, Zuccarella-Hackl, Bürkle, Ehlert and Wirtz2018), array DNAm data and standard procedures were used to estimate blood cell proportions (CD8 + T cells, CD4 + T cells, B cells, natural killer cells, granulocytes, monocytes) (Houseman et al., Reference Houseman, Accomando, Koestler, Christensen, Marsit, Nelson and Kelsey2012) that were included as covariates to adjust for potential confounding in all regression models.

Genetic ancestry principal components

DNA samples were genotyped using the Infinium Global Screening Array-24 v1.0 (Illumina Inc.) at the Stanley Center/Broad Institute. Data were quality-controlled and principal component analysis (PCA) was implemented using the Plink1.9 program. We performed PCA jointly on our samples with samples from the 1000 Genomes Project (1000G) and extracted the top 10 principal components (PCs) (plink.eigenvec and plink.eigenval files were generated) of the variance-standardized relationship matrix. With the 1000G PCs, we trained a decision tree and used the rpart (Recursive Partitioning and Regression Trees) package in RStudio (Version 2022.02.3) to test AURORA PCs to predicate the genetic ancestry and the possible genetic ancestry probabilities.

Neuroimaging measures

T1-weighted structural MRI and resting-state functional MRI (rs-fMRI) data were collected across five sites with harmonized acquisition parameters at the 2-week follow-up timepoint (Harnett et al., Reference Harnett, van Rooij, Ely, Lebois, Murty, Jovanovic and Stevens2021). Data were preprocessed using FMRIPREP v1.2.2 (Esteban et al., Reference Esteban, Markiewicz, Blair, Moodie, Isik, Erramuzpe and Gorgolewski2019). Brain surfaces were reconstructed, and subcortical volumes were extracted using FreeSurfer v6.0.1 (Dale, Fischl, & Sereno, Reference Dale, Fischl and Sereno1999), in order to generate volumetric data for the hippocampus and the amygdala (left and right). Each hippocampus was subdivided into 21 subregions (Iglesias et al., Reference Iglesias, Augustinack, Nguyen, Player, Player, Wright and Van Leemput2015) and each amygdala into nine nuclei (Saygin et al., Reference Saygin, Kliemann, Iglesias, van der Kouwe, Boyd, Reuter and Augustinack2017). For each of these metrics and to minimize the number of comparisons, analyses examined the sum volume of both sides (left and right) as the variable of interest. The rs-fMRI data (TR = 2.36 s, 230 volumes, 9:05 min scan time) were processed using ICA-AROMA as part of the FMRIPREP pipeline, which has been shown to handle motion artifacts in a robust, data-driven fashion that performs equal to or better than standard scrubbing or censoring procedures (Pruim et al., Reference Pruim, Mennes, van Rooij, Llera, Buitelaar and Beckmann2015b; Pruim, Mennes, Buitelaar, & Beckmann, Reference Pruim, Mennes, Buitelaar and Beckmann2015a). The rs-fMRI data were further processed within the Analysis for Functional NeuroImages program 3dTproject to perform linear detrending, censoring of non-steady state volumes identified by FMRIPREP, bandpass filtering (0.01–0.1 Hz), and regression of white matter, corticospinal fluid, and global signal to account for potential physiological noise. Network connectivity was estimated by correlating the mean fMRI time-course from regions of interest (ROIs) in the Yeo 7-Network atlas (Yeo et al., Reference Yeo, Krienen, Sepulcre, Sabuncu, Lashkari, Hollinshead and Buckner2011). Independent Pearson correlation coefficients were calculated for each pair of ROIs to represent the strength of network-to-network functional connectivity. Pearson correlations were z-transformed prior to statistical analyses.

Statistical analysis

All statistical analyses were performed in R version 4.2.0. Logistic regression models tested ED epigenetic aging markers as the primary (four total) independent variables and 6-month probable PTSD diagnosis as the primary dependent variable of interest. These primary analyses were corrected for multiple testing with the conservative Bonferroni method (Bland & Altman, Reference Bland and Altman1995), leading to an adjusted α (p value threshold) of 0.05/4 = 0.0125. Secondary analyses evaluated associations between epigenetic aging and distinct PTSD-related symptom trajectories and between individual GrimAge CpG sites and PTSD outcomes. Other secondary analyses evaluated associations between epigenetic aging and amygdala and hippocampus volumes as well as resting-state functional connectivity. These analyses were corrected for multiple testing using the false discovery rate (FDR) method (Benjamini & Hochberg, Reference Benjamini and Hochberg1995). All models were adjusted for age, sex, self-reported race/ethnicity, educational level, marital status, and DNAm-estimated blood cell proportions. We further controlled all models involving MRI data for imaging site and models involving structural MRI for total intracranial volume. Additional sensitivity analyses adjusted for CTQ and lifetime trauma scores, PTSD symptoms at the time of ED presentation and over the lifetime, SF-12 mental and physical health scores, quantity of tobacco and alcohol use, income, body mass index (BMI), and the first three genetic ancestry (SNP-based) PCs.

Results

Cohort overview and clinicodemographic characteristics

The characteristics of participants (n = 289) are summarized in Table 1. Most participants were non-Hispanic Black women less than 40 years of age. As expected, all four epigenetic aging markers (HorvathAge, HannumAge, PhenoAge, and GrimAge) were strongly correlated with chronological age (all pairwise correlation p values < 2.2 × 10−16). Approximately one in four participants had probable PTSD 6 months after trauma. There were no significant differences in chronological age, sex, or race/ethnicity between participants with and without 6-month PTSD (all p values > 0.46).

Advanced GrimAge at the time of trauma predicts PTSD outcomes in the ensuing 6 months

Primary analyses tested whether any of the four measures of epigenetic aging (HorvathAge, HannumAge, PhenoAge, and GrimAge) predicted 6-month PTSD outcome. After adjusting for age, sex, self-reported race/ethnicity, educational level, marital status, and DNAm-estimated blood cell proportions, and correcting for multiple comparisons (Bonferroni-adjusted α < 0.0125), only ED GrimAge significantly predicted 6-month probable PTSD diagnosis (n = 244, β = 0.11, s.e. = 0.04, z = 2.5, p = 0.0114). This effect remained significant after further stepwise adjustment for childhood and lifetime trauma burden, PTSD symptoms at the time of ED presentation and over the lifetime, general mental and physical health, quantity of tobacco and alcohol use, income, BMI, and genetic ancestry PCs (p values between 0.0092 and 0.0414). Moreover, a dose–response relationship between epigenetic aging and PTSD was observed, with individuals in the highest GrimAge tertile having 17% greater risk than those in the medium and 44% greater risk than individuals in the lowest tertile (Fig. 1a, b). Secondary analyses assessed whether ED GrimAge predicts distinct PTSD-related symptom trajectories: intrusive memories, hyperarousal, avoidance, nightmares, and sleep disturbance. As expected, pairwise positive correlations were observed for scores obtained for all five symptom categories at each of the 10 timepoints (all p values < 4.8 × 10−7). In adjusted analyses, advanced GrimAge significantly predicted worse trajectories of intrusive memories (n = 289, β = 0.10, s.e. = 0.04, z = 2.5, p = 0.0125) and nightmares (n = 289, β = 0.09, s.e. = 0.04, z = 2.1, p = 0.0319) during the 6 months after trauma, and these predictions did not change after further controlling for PTSD symptoms at ED presentation (p = 0.0177 for intrusive memories and 0.0405 for nightmares). In contrast, no significant findings were observed for hyperarousal, avoidance, and sleep (p values > 0.24). Taken together, these findings suggest that advanced epigenetic age is associated with heightened PTSD risk and, in particular, increased intrusive memories and nightmares in the aftermath of trauma.

Fig. 1.

Advanced GrimAge at emergency department (ED) presentation predicts increased risk for posttraumatic stress disorder (PTSD) 6 months later. (a) Proportion of participants with 6-month probable PTSD stratified by ED GrimAge tertile. (b) Graph depicting how ED GrimAge predicts the probability for 6-month PTSD diagnosis. Predicted probabilities have been adjusted for covariates (details in Methods section). Depicted GrimAge residuals are after adjusting for chronological age. Additional statistics are provided in Results.

GrimAge CpG sites and PTSD outcomes

The prediction of PTSD trajectories by the composite (multi-CpG) GrimAge marker also prompted us to examine whether it may be driven by specific CpG sites. To address this question, we separately tested if 6-month probable PTSD was associated with ED methylation levels at each of the 1030 CpGs that comprise GrimAge (Lu et al., Reference Lu, Quach, Wilson, Reiner, Aviv, Raj and Horvath2019). After adjusting for age, sex, self-reported race/ethnicity, educational level, marital status, and DNAm-estimated blood cell proportions, a total of 39 CpGs sites were found to predict PTSD at the nominal but not the FDR-adjusted level of statistical significance (online Supplementary Table S4). Among these 39 CpG sites, ED methylation at only two sites was also significantly associated with the GrimAge-predicted intrusive memories and nightmares symptom trajectories (cg06722193, cg02716826; online Supplementary Table S4). To rule out the possibility that associations merely reflect correlation strength of these CpGs with GrimAge in our data, we performed pairwise correlations of each GrimAge CpG with the composite GrimAge metric. When ranking correlations from smaller to larger p values, cg06722193 and cg02716826 were respectively ranked #286 and #53 among all GrimAge CpGs. The association with cg06722193 (shown in Fig. 2) is particularly interesting, given that this CpG is located within IRX6, a gene that is involved in neuronal development and is epigenetically regulated by stress exposure (Del Corvo et al., Reference Del Corvo, Bongiorni, Stefanon, Sgorlon, Valentini, Ajmone Marsan and Chillemi2020; Leung et al., Reference Leung, Kosyk, Welter, Dietrich, Archer and Zannas2022; Star et al., Reference Star, Zhu, Shi, Liu, Pashmforoush, Sauve and Chow2012). These findings provide limited evidence for genomic site-specific prediction of PTSD outcomes and rather support this prediction as an overall (i.e. emergent) property of the composite GrimAge marker.

Fig. 2.

Example of a GrimAge CpG, the methylation levels of which predict 6-month posttraumatic stress disorder (PTSD) and related symptom trajectories at the nominal level of statistical significance. Predicted probabilities have been adjusted for covariates (as described in Methods). Detailed statistics are provided in online Supplementary Table S4.

Advanced GrimAge is associated with brain alterations relevant for PTSD outcomes

To identify posttraumatic brain alterations associated with advanced epigenetic age (assessed via ED GrimAge), we leveraged structural and functional neuroimaging (MRI) data available in a smaller subset of individuals at the 2-week follow-up timepoint (n = 63; 40 women, 23 men; age mean, 35.1; age range = 18–67). Structural MRI analyses focused on the amygdala and the hippocampus, two brain regions with established roles in stress and trauma-related phenotypes (Del Casale et al., Reference Del Casale, Ferracuti, Barbetti, Bargagna, Zega, Iannuccelli and Pompili2022; Morey et al., Reference Morey, Gold, LaBar, Beall, Brown, Haswell and McCarthy2012). After adjusting for age, sex, self-reported race/ethnicity, educational level, marital status, DNAm-estimated blood cell proportions, imaging site, and total intracranial volume, advanced GrimAge was significantly associated with reduced whole amygdala volume (n = 62, β = −22.7, s.e. = 10.8, t = −2.1, p = 0.0413) but not hippocampal volume (p = 0.21). Secondary analyses examined if this association was driven by selected amygdala subregions (nine measures total). After covariate adjustment and FDR correction for multiple testing, advanced GrimAge was significantly associated with smaller volumes of the cortico-amygdaloid transition and the cortical and accessory basal nuclei (FDR-adjusted p values between 0.0013 and 0.0169; example shown in Fig. 3). We next explored associations between advanced GrimAge and resting-state functional network connectivity (21 measures total). After covariate adjustment and FDR correction, GrimAge was not significantly associated with any network connectivity measure (n = 63; all FDR-adjusted p values > 0.07).

Fig. 3.

Advanced GrimAge is associated with reduced amygdala subregion volumes. Example is shown for volume of the accessory basal nucleus of the amygdala. The brain image on the right depicts the accessory basal nucleus from a single subject. Depicted GrimAge residuals are after adjusting for chronological age and DNA methylation-estimated blood cell proportions. Additional statistics are provided in Results.

Discussion

Prior research has linked stress exposure and PTSD with advanced epigenetic age (Belsky et al., Reference Belsky, Caspi, Corcoran, Sugden, Poulton, Arseneault and Moffitt2022; Boks et al., Reference Boks, van Mierlo, Rutten, Radstake, De Witte, Geuze and Vermetten2015; Brody et al., Reference Brody, Yu, Chen, Beach and Miller2016; Copeland et al., Reference Copeland, Shanahan, McGinnis, Aberg and van den Oord2022; Harvanek et al., Reference Harvanek, Fogelman, Xu and Sinha2021; Katrinli et al., Reference Katrinli, Stevens, Wani, Lori, Kilaru, van Rooij and Smith2020; Kuan et al., Reference Kuan, Ren, Clouston, Yang, Jonas, Kotov and Luft2021; Lim et al., Reference Lim, Nzegwu and Wright2022; Mehta et al., Reference Mehta, Bruenig, Pierce, Sathyanarayanan, Stringfellow, Miller and Shakespeare-Finch2022; Na et al., Reference Na, Montalvo-Ortiz, Nagamatsu, Southwick, Krystal, Gelernter and Pietrzak2022; Wang et al., Reference Wang, Hui, Goldberg, Smith, Kaseer, Murrah and Sun2022; Wolf et al., Reference Wolf, Maniates, Nugent, Maihofer, Armstrong, Ratanatharathorn and Logue2018; Zannas et al., Reference Zannas, Arloth, Carrillo-Roa, Iurato, Roh, Ressler and Mehta2015a), but no studies have examined whether epigenetic aging measured at the time of trauma predicts subsequent development of PTSD outcomes. Leveraging the AURORA cohort (McLean et al., Reference McLean, Ressler, Koenen, Neylan, Germine, Jovanovic and Kessler2020), the present study showed that advanced epigenetic age (measured with GrimAge) at ED presentation predicts increased PTSD risk and worse intrusive memories and nightmares symptom trajectories in the ensuing 6 months. In a cohort subset with neuroimaging data 2 weeks after trauma, advanced ED epigenetic age was further associated with reduced whole amygdala and amygdala subregion volumes.

Recent work in the AURORA study has identified promising clinical predictors (Kessler et al., Reference Kessler, Ressler, House, Beaudoin, An, Stevens and McLean2021; Ziobrowski et al., Reference Ziobrowski, Kennedy, Ustun, House, Beaudoin, An and van Rooij2021), but personalized interventions would also benefit from molecular markers of risk for developing distinct posttraumatic outcomes (Howie, Rijal, & Ressler, Reference Howie, Rijal and Ressler2019; Linnstaedt, Zannas, McLean, Koenen, & Ressler, Reference Linnstaedt, Zannas, McLean, Koenen and Ressler2020; Smith et al., Reference Smith, Ratanatharathorn, Maihofer, Naviaux, Aiello, Amstadter and Nievergelt2020; Zannas et al., Reference Zannas, Provençal and Binder2015b). Epigenetic signatures have been proposed as prime candidate markers of posttraumatic vulnerability (Howie et al., Reference Howie, Rijal and Ressler2019; Linnstaedt et al., Reference Linnstaedt, Zannas, McLean, Koenen and Ressler2020; Smith et al., Reference Smith, Ratanatharathorn, Maihofer, Naviaux, Aiello, Amstadter and Nievergelt2020; Zannas et al., Reference Zannas, Provençal and Binder2015b), given the epigenome's role as a molecular interface between environment and health (Cavalli & Heard, Reference Cavalli and Heard2019; Gassen et al., Reference Gassen, Chrousos, Binder and Zannas2016). Our finding that GrimAge is a predictive marker of PTSD extends previous studies observing advanced GrimAge in individuals with current or lifetime PTSD (Katrinli et al., Reference Katrinli, Stevens, Wani, Lori, Kilaru, van Rooij and Smith2020; Na et al., Reference Na, Montalvo-Ortiz, Nagamatsu, Southwick, Krystal, Gelernter and Pietrzak2022; Wolf et al., Reference Wolf, Maniates, Nugent, Maihofer, Armstrong, Ratanatharathorn and Logue2018; Yang et al., Reference Yang, Wu, Verhoeven, Gautam, Reus, Kang and Wolkowitz2021). Moreover, our prospective observations build on previous longitudinal studies indicating that both trauma exposure and increased PTSD symptoms are associated with accelerated epigenetic aging over time (Belsky et al., Reference Belsky, Caspi, Corcoran, Sugden, Poulton, Arseneault and Moffitt2022; Mehta et al., Reference Mehta, Bruenig, Pierce, Sathyanarayanan, Stringfellow, Miller and Shakespeare-Finch2022; Sumner, Colich, Uddin, Armstrong, & McLaughlin, Reference Sumner, Colich, Uddin, Armstrong and McLaughlin2019; Wolf et al., Reference Wolf, Logue, Morrison, Wilcox, Stone, Schichman and Miller2019; Yang et al., Reference Yang, Wu, Verhoeven, Gautam, Reus, Kang and Wolkowitz2021). The unique ability of GrimAge to capture vulnerability in the immediate aftermath of trauma may stem from its development as a predictor of healthspan and lifespan, which likely makes it more amenable to environmental stressors as life advances. In contrast, we speculate that DNAm-based predictors of chronological age may undergo more tightly programmed epigenomic changes as a result of advancing age. Intriguingly, our longitudinal data spanning distinct symptom categories further show that GrimAge predicts worse trajectories only for intrusive memories and nightmares, thereby suggesting that advanced epigenetic age contributes to psychiatric risk through select PTSD symptoms that develop and persist longitudinally after trauma exposure.

Leveraging neuroimaging data available in a subset of study participants, we also found that GrimAge is associated with reduced volume of the whole amygdala and specific amygdala subregions, including the cortico-amygdaloid transition and the cortical and accessory basal nuclei. Reduced volume in the whole amygdala and in select amygdala subregions has been previously observed in PTSD (Morey et al., Reference Morey, Gold, LaBar, Beall, Brown, Haswell and McCarthy2012, Reference Morey, Clarke, Haswell, Phillips, Clausen, Mufford and LaBar2020). Work in animal models further shows that amygdala subregions can shrink in stress-exposed mice and further predispose to exacerbated behavioral sequelae after stress exposure (Golub et al., Reference Golub, Kaltwasser, Mauch, Herrmann, Schmidt, Holsboer and Wotjak2011; Yang et al., Reference Yang, Mozhui, Karlsson, Cameron, Williams and Holmes2008). Our findings thus suggest that advanced epigenetic age is associated with structural alterations in the amygdala and related increased vulnerability for PTSD development and persistence. This possibility is congruent with a prior study linking GrimAge with brain region-specific cortical atrophy (Katrinli et al., Reference Katrinli, Stevens, Wani, Lori, Kilaru, van Rooij and Smith2020). It is important to note that all neuroimaging measures in the present study were obtained 2 weeks after ED presentation. While this timepoint was in part selected due to challenges inherent to performing MRI at the time of trauma exposure, it also lies temporally between ED presentation and the 6-month follow-up. Thus, an intriguing hypothesis is that early structural brain alterations associated with advanced GrimAge could predispose to worse PTSD trajectories in the aftermath of trauma. However, the current study design limited our ability to test this hypothesis, given the lack of neuroimaging measures before trauma exposure that would be necessary to temporally disentangle the observed associations.

Additional limitations should be considered when interpreting the findings reported herein. Although our analyses adjusted for several potential confounders, the study's unique design and observational setting did not allow us to include a control (non-trauma) group that would disentangle the extent to which advanced GrimAge was a direct consequence of or already present before the traumatic event. However, given that blood samples were collected within hours of trauma exposure and neuroimaging was performed at the 2-week follow-up, it is likely that epigenetic patterns and brain alterations were already present prior to trauma. Moreover, the study design precluded us from testing whether PTSD symptoms accelerate epigenetic aging, which is a more commonly studied direction of association. Analyses were adjusted for key potential confounders, including childhood and lifetime trauma burden, ED and lifetime PTSD symptoms, and general mental and physical health, but the possibility that other undocumented confounders could in part account for the observed associations cannot be ruled out. As expected, all symptom scores showed significant positive pairwise correlations and thus do not represent independent signals; however, our findings suggest that advanced GrimAge is specifically associated with worse intrusive memories and nightmares symptom trajectories. Our sample size was modest, especially for analyses involving neuroimaging data. This limited our power for conclusively testing genomic site-specific predictions and precluded us from testing if brain alterations statistically mediate the prediction of PTSD outcomes by GrimAge. Epigenetic assessments were conducted in whole blood, and while our analyses adjusted for blood cell composition, the findings' mechanistic relevance remains to be dissected in brain tissues with direct phenotypic relevance, such as the amygdala. The presented findings will thus benefit by replication and further dissection in larger independent cohorts and postmortem datasets.

In sum, the findings presented here shed new light on the relation between biological aging and trauma-related phenotypes, suggesting that GrimAge measured immediately after trauma predicts subsequent PTSD trajectories and is associated with relevant brain alterations. Furthering these findings has the potential to enhance early prevention and treatment of posttraumatic psychiatric sequelae.

Supplementary material

The supplementary material for this article can be found at https://doi.org/10.1017/S0033291723000636

Acknowledgements

The investigators wish to thank the trauma survivors participating in the AURORA Study. Their time and effort during a challenging period of their lives make our efforts to improve recovery for future trauma survivors possible. This project was supported by NIMH under U01MH110925, the US Army MRMC, One Mind, and The Mayday Fund. This work was also supported by a NARSAD Young Investigator grant (#24135) and a Foundation of Hope for Research and Treatment of Mental Illness grant awarded to ASZ. Verily Life Sciences and Mindstrong Health provided some of the hardware and software used to perform study assessments. Bio-samples for this publication were processed by and obtained from the NIMH Repository & Genomics Resource, supported by cooperative agreement U24MH068457. Data and/or research tools used in the preparation of this manuscript were obtained from the National Institute of Mental Health (NIMH) Data Archive (NDA) under DOI: 10.15154/1528075. NDA is a collaborative informatics system created by the National Institutes of Health to provide a national resource to support and accelerate research in mental health. This manuscript reflects the views of the authors and may not reflect the opinions or views of any of the funders or of the Submitters submitting original data to NDA. The authors would like to thank Dr Ake T. Lu and Dr. Steve Horvath for providing the list of CpG sites that comprise the GrimAge marker.

Author contributions

Conceptualization, analyses, and manuscript draft: ASZ. Genomic data: SDL. Phenotypic data and statistical supervision: XA. Neuroimaging data: JSS, NGH, ARR, KIO. Critical editing and supervision: DRR, EBB, KCK, KJR, SAM. Funding: ASZ, KCK, KJR, SAM. Manuscript editing and approval: all authors.

Financial support

This project was supported by NIMH under U01MH110925, the US Army MRMC, One Mind, and The Mayday Fund. This work was also supported by a NARSAD Young Investigator grant (#24135) and a Foundation of Hope for Research and Treatment of Mental Illness grant awarded to ASZ. Bio-samples for this publication were processed by and obtained from the NIMH Repository & Genomics Resource, supported by cooperative agreement U24MH068457.

Conflict of interest

None.

Ethical standards

The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008.

References

Adalsteinsson, B. T., Gudnason, H., Aspelund, T., Harris, T. B., Launer, L. J., Eiriksdottir, G., … Gudnason, V. (2012). Heterogeneity in white blood cells has potential to confound DNA methylation measurements. PLoS ONE, 7(10), e46705. doi: 10.1371/journal.pone.0046705CrossRefGoogle ScholarPubMed
Beaudoin, F. L., An, X., Basu, A., Ji, Y., Liu, M., Kessler, R. C., … McLean, S. A. (2023). Use of serial smartphone-based assessments to characterize diverse neuropsychiatric symptom trajectories in a large trauma survivor cohort. Translational Psychiatry, 13(1), 4. doi: 10.1038/s41398-022-02289-yCrossRefGoogle Scholar
Beis, D., von Känel, R., Heimgartner, N., Zuccarella-Hackl, C., Bürkle, A., Ehlert, U., & Wirtz, P. H. (2018). The role of norepinephrine and α-adrenergic receptors in acute stress-induced changes in granulocytes and monocytes. Psychosomatic Medicine, 80(7), 649–658.CrossRefGoogle ScholarPubMed
Belsky, D. W., Caspi, A., Corcoran, D. L., Sugden, K., Poulton, R., Arseneault, L., … Moffitt, T. E. (2022). DunedinPACE, a DNA methylation biomarker of the pace of aging. Elife, 11, e73420. doi: 10.7554/eLife.73420CrossRefGoogle ScholarPubMed
Benjamini, Y., & Hochberg, Y. (1995). Controlling the false discovery rate: A practical and powerful approach to multiple testing. Journal of the Royal statistical society: series B (Methodological), 57(1), 289–300.Google Scholar
Bernstein, D. P., Stein, J. A., Newcomb, M. D., Walker, E., Pogge, D., Ahluvalia, T., … Zule, W. (2003). Development and validation of a brief screening version of the Childhood Trauma Questionnaire. Child Abuse and Neglect, 27(2), 169–190. doi: 10.1016/s0145-2134(02)00541-0CrossRefGoogle ScholarPubMed
Bersani, F. S., Mellon, S. H., Reus, V. I., & Wolkowitz, O. M. (2019). Accelerated aging in serious mental disorders. Current Opinion in Psychiatry, 32(5), 381–387. doi: 10.1097/yco.0000000000000525CrossRefGoogle ScholarPubMed
Bland, J. M., & Altman, D. G. (1995). Multiple significance tests: The Bonferroni method. British Medical Journal, 310(6973), 170. doi: 10.1136/bmj.310.6973.170CrossRefGoogle ScholarPubMed
Blevins, C. A., Weathers, F. W., Davis, M. T., Witte, T. K., & Domino, J. L. (2015). The posttraumatic stress disorder checklist for DSM-5 (PCL-5): Development and initial psychometric evaluation. Journal of Traumatic Stress, 28(6), 489–498. doi: 10.1002/jts.22059CrossRefGoogle ScholarPubMed
Boks, M. P., van Mierlo, H. C., Rutten, B. P., Radstake, T. R., De Witte, L., Geuze, E., … Vermetten, E. (2015). Longitudinal changes of telomere length and epigenetic age related to traumatic stress and post-traumatic stress disorder. Psychoneuroendocrinology, 51, 506–512. doi: 10.1016/j.psyneuen.2014.07.011CrossRefGoogle ScholarPubMed
Bovin, M. J., Marx, B. P., Weathers, F. W., Gallagher, M. W., Rodriguez, P., Schnurr, P. P., & Keane, T. M. (2016). Psychometric properties of the PTSD checklist for diagnostic and statistical manual of mental disorders-fifth edition (PCL-5) in veterans. Psychological Assessment, 28(11), 1379–1391. doi: 10.1037/pas0000254CrossRefGoogle ScholarPubMed
Breslau, N., Kessler, R. C., Chilcoat, H. D., Schultz, L. R., Davis, G. C., & Andreski, P. (1998). Trauma and posttraumatic stress disorder in the community: The 1996 Detroit area survey of trauma. Archives of General Psychiatry, 55(7), 626–632. doi: 10.1001/archpsyc.55.7.626CrossRefGoogle ScholarPubMed
Brody, G. H., Yu, T., Chen, E., Beach, S. R., & Miller, G. E. (2016). Family-centered prevention ameliorates the longitudinal association between risky family processes and epigenetic aging. Journal of Child Psychology and Psychiatry and Allied Disciplines, 57(5), 566–574. doi: 10.1111/jcpp.12495CrossRefGoogle ScholarPubMed
Cavalli, G., & Heard, E. (2019). Advances in epigenetics link genetics to the environment and disease. Nature, 571(7766), 489–499. doi: 10.1038/s41586-019-1411-0CrossRefGoogle Scholar
Copeland, W. E., Shanahan, L., McGinnis, E. W., Aberg, K. A., & van den Oord, E. (2022). Early adversities accelerate epigenetic aging into adulthood: A 10-year, within-subject analysis. Journal of Child Psychology and Psychiatry and Allied Disciplines, 63(11), 1308–1315. doi: 10.1111/jcpp.13575CrossRefGoogle ScholarPubMed
Dale, A. M., Fischl, B., & Sereno, M. I. (1999). Cortical surface-based analysis. I. Segmentation and surface reconstruction. Neuroimage, 9(2), 179–194. doi: 10.1006/nimg.1998.0395CrossRefGoogle ScholarPubMed
Del Casale, A., Ferracuti, S., Barbetti, A. S., Bargagna, P., Zega, P., Iannuccelli, A., … Pompili, M. (2022). Grey matter volume reductions of the left hippocampus and amygdala in PTSD: A coordinate-based meta-analysis of magnetic resonance imaging studies. Neuropsychobiology, 81(4), 257–264. doi: 10.1159/000522003CrossRefGoogle ScholarPubMed
Del Corvo, M., Bongiorni, S., Stefanon, B., Sgorlon, S., Valentini, A., Ajmone Marsan, P., … Chillemi, G. (2020). Genome-wide DNA methylation and gene expression profiles in cows subjected to different stress level as assessed by cortisol in milk. Genes (Basel), 11(8), 850. doi: 10.3390/genes11080850CrossRefGoogle ScholarPubMed
Esteban, O., Markiewicz, C. J., Blair, R. W., Moodie, C. A., Isik, A. I., Erramuzpe, A., … Gorgolewski, K. J. (2019). fMRIPrep: A robust preprocessing pipeline for functional MRI. Nature Methods, 16(1), 111–116. doi: 10.1038/s41592-018-0235-4CrossRefGoogle ScholarPubMed
Felitti, V. J., Anda, R. F., Nordenberg, D., Williamson, D. F., Spitz, A. M., Edwards, V., … Marks, J. S. (1998). Relationship of childhood abuse and household dysfunction to many of the leading causes of death in adults. The adverse childhood experiences (ACE) study. American Journal of Preventive Medicine, 14(4), 245–258.CrossRefGoogle ScholarPubMed
Fraga, M. F., Ballestar, E., Paz, M. F., Ropero, S., Setien, F., Ballestar, M. L., … Esteller, M. (2005). Epigenetic differences arise during the lifetime of monozygotic twins. Proceedings of the National Academy of Sciences of the USA, 102(30), 10604–10609. doi: 10.1073/pnas.0500398102CrossRefGoogle ScholarPubMed
Gassen, N. C., Chrousos, G. P., Binder, E. B., & Zannas, A. S. (2016). Life stress, glucocorticoid signaling, and the aging epigenome: Implications for aging-related diseases. Neuroscience and Biobehavioral Reviews, 74(Pt B), 356–365. doi: 10.1016/j.neubiorev.2016.06.003CrossRefGoogle ScholarPubMed
Golub, Y., Kaltwasser, S. F., Mauch, C. P., Herrmann, L., Schmidt, U., Holsboer, F., … Wotjak, C. T. (2011). Reduced hippocampus volume in the mouse model of posttraumatic stress disorder. Journal of Psychiatric Research, 45(5), 650–659. doi: 10.1016/j.jpsychires.2010.10.014CrossRefGoogle ScholarPubMed
Gray, M. J., Litz, B. T., Hsu, J. L., & Lombardo, T. W. (2004). Psychometric properties of the life events checklist. Assessment, 11(4), 330–341. doi: 10.1177/1073191104269954CrossRefGoogle ScholarPubMed
Hamilton, C. M., Strader, L. C., Pratt, J. G., Maiese, D., Hendershot, T., Kwok, R. K., … Haines, J. (2011). The PhenX Toolkit: Get the most from your measures. American Journal of Epidemiology, 174(3), 253–260. doi: 10.1093/aje/kwr193CrossRefGoogle ScholarPubMed
Hannum, G., Guinney, J., Zhao, L., Zhang, L., Hughes, G., Sadda, S., … Zhang, K. (2013). Genome-wide methylation profiles reveal quantitative views of human aging rates. Molecular Cell, 49(2), 359–367. doi: 10.1016/j.molcel.2012.10.016CrossRefGoogle ScholarPubMed
Harnett, N. G., van Rooij, S. J. H., Ely, T. D., Lebois, L. A. M., Murty, V. P., Jovanovic, T., … Stevens, J. S. (2021). Prognostic neuroimaging biomarkers of trauma-related psychopathology: Resting-state fMRI shortly after trauma predicts future PTSD and depression symptoms in the AURORA study. Neuropsychopharmacology, 46(7), 1263–1271. doi: 10.1038/s41386-020-00946-8CrossRefGoogle ScholarPubMed
Harvanek, Z. M., Fogelman, N., Xu, K., & Sinha, R. (2021). Psychological and biological resilience modulates the effects of stress on epigenetic aging. Translational Psychiatry, 11(1), 601. doi: 10.1038/s41398-021-01735-7CrossRefGoogle ScholarPubMed
Hillary, R. F., Stevenson, A. J., McCartney, D. L., Campbell, A., Walker, R. M., Howard, D. M., … Marioni, R. E. (2020). Epigenetic measures of ageing predict the prevalence and incidence of leading causes of death and disease burden. Clinical Epigenetics, 12(1), 115. doi: 10.1186/s13148-020-00905-6CrossRefGoogle ScholarPubMed
Horvath, S. (2013). DNA methylation age of human tissues and cell types. Genome Biology, 14(10), R115. doi: 10.1186/gb-2013-14-10-r115CrossRefGoogle ScholarPubMed
Horvath, S., & Raj, K. (2018). DNA methylation-based biomarkers and the epigenetic clock theory of ageing. Nature Reviews Genetics, 19(6), 371–384. doi: 10.1038/s41576-018-0004-3CrossRefGoogle ScholarPubMed
Houseman, E. A., Accomando, W. P., Koestler, D. C., Christensen, B. C., Marsit, C. J., Nelson, H. H., … Kelsey, K. T. (2012). DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinformatics, 13, 86. doi: 10.1186/1471-2105-13-86CrossRefGoogle ScholarPubMed
Howie, H., Rijal, C. M., & Ressler, K. J. (2019). A review of epigenetic contributions to post-traumatic stress disorder. Dialogues in Clinical Neuroscience, 21(4), 417–428. doi: 10.31887/DCNS.2019.21.4/kresslerCrossRefGoogle ScholarPubMed
Hugo, V. (2012). Les miserables. San Diego, CA: Canterbury Classics.Google Scholar
Iglesias, J. E., Augustinack, J. C., Nguyen, K., Player, C. M., Player, A., Wright, M., … Van Leemput, K. (2015). A computational atlas of the hippocampal formation using ex vivo, ultra-high resolution MRI: Application to adaptive segmentation of in vivo MRI. Neuroimage, 115, 117–137. doi: 10.1016/j.neuroimage.2015.04.042CrossRefGoogle ScholarPubMed
Johnson, W. E., Li, C., & Rabinovic, A. (2007). Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics, 8(1), 118–127. doi: 10.1093/biostatistics/kxj037CrossRefGoogle ScholarPubMed
Joyce, B., Gao, T., Zheng, Y., Ma, J., Hwang, S. J., Liu, L., … Lloyd-Jones, D. (2021). Epigenetic age acceleration reflects long-term cardiovascular health. Circulation Research, 129(8), 770–781. doi: 10.1161/circresaha.121.318965CrossRefGoogle ScholarPubMed
Kao, C. Y., He, Z., Zannas, A. S., Hahn, O., Kühne, C., Reichel, J. M., … Turck, C. W. (2016). Fluoxetine treatment prevents the inflammatory response in a mouse model of posttraumatic stress disorder. Journal of Psychiatric Research, 76, 74–83. doi: 10.1016/j.jpsychires.2016.02.003CrossRefGoogle Scholar
Katrinli, S., Stevens, J., Wani, A. H., Lori, A., Kilaru, V., van Rooij, S. J. H., … Smith, A. K. (2020). Evaluating the impact of trauma and PTSD on epigenetic prediction of lifespan and neural integrity. Neuropsychopharmacology, 45(10), 1609–1616. doi: 10.1038/s41386-020-0700-5CrossRefGoogle ScholarPubMed
Kessler, R. C., Ressler, K. J., House, S. L., Beaudoin, F. L., An, X., Stevens, J. S., … McLean, S. A. (2021). Socio-demographic and trauma-related predictors of PTSD within 8 weeks of a motor vehicle collision in the AURORA study. Molecular Psychiatry, 26(7), 3108–3121. doi: 10.1038/s41380-020-00911-3CrossRefGoogle ScholarPubMed
Korgaonkar, M. S., Chakouch, C., Breukelaar, I. A., Erlinger, M., Felmingham, K. L., Forbes, D., … Bryant, R. A. (2020). Intrinsic connectomes underlying response to trauma-focused psychotherapy in post-traumatic stress disorder. Translational Psychiatry, 10(1), 270. doi: 10.1038/s41398-020-00938-8CrossRefGoogle ScholarPubMed
Kuan, P. F., Ren, X., Clouston, S., Yang, X., Jonas, K., Kotov, R., … Luft, B. J. (2021). PTSD is associated with accelerated transcriptional aging in World Trade Center responders. Translational Psychiatry, 11(1), 311. doi: 10.1038/s41398-021-01437-0CrossRefGoogle ScholarPubMed
Lang, A. J., & Stein, M. B. (2005). An abbreviated PTSD checklist for use as a screening instrument in primary care. Behaviour Research and Therapy, 43(5), 585–594. doi: 10.1016/j.brat.2004.04.005CrossRefGoogle ScholarPubMed
Leek, J. T., Johnson, W. E., Parker, H. S., Jaffe, A. E., & Storey, J. D. (2012). The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics, 28(6), 882–883. doi: 10.1093/bioinformatics/bts034CrossRefGoogle ScholarPubMed
Leung, C. S., Kosyk, O., Welter, E. M., Dietrich, N., Archer, T. K., & Zannas, A. S. (2022). Chronic stress-driven glucocorticoid receptor activation programs key cell phenotypes and functional epigenomic patterns in human fibroblasts. iScience, 25(9), 104960. doi: 10.1016/j.isci.2022.104960CrossRefGoogle ScholarPubMed
Levine, M. E., Lu, A. T., Quach, A., Chen, B. H., Assimes, T. L., Bandinelli, S., … Horvath, S. (2018). An epigenetic biomarker of aging for lifespan and healthspan. Aging, 10(4), 573–591. doi: 10.18632/aging.101414CrossRefGoogle ScholarPubMed
Lim, S., Nzegwu, D., & Wright, M. L. (2022). The impact of psychosocial stress from life trauma and racial discrimination on epigenetic aging-A systematic review. Biological Research for Nursing, 24(2), 202–215. doi: 10.1177/10998004211060561CrossRefGoogle ScholarPubMed
Linnstaedt, S. D., Zannas, A. S., McLean, S. A., Koenen, K. C., & Ressler, K. J. (2020). Literature review and methodological considerations for understanding circulating risk biomarkers following trauma exposure. Molecular Psychiatry, 25(9), 1986–1999. doi: 10.1038/s41380-019-0636-5CrossRefGoogle ScholarPubMed
Lopez-Otin, C., Blasco, M. A., Partridge, L., Serrano, M., & Kroemer, G. (2013). The hallmarks of aging. Cell, 153(6), 1194–1217. doi: 10.1016/j.cell.2013.05.039CrossRefGoogle ScholarPubMed
López-Otín, C., Blasco, M. A., Partridge, L., Serrano, M., & Kroemer, G. (2022). Hallmarks of aging: An expanding universe. Cell, 186(2), 243–278. doi: 10.1016/j.cell.2022.11.001CrossRefGoogle Scholar
Lu, A. T., Quach, A., Wilson, J. G., Reiner, A. P., Aviv, A., Raj, K., … Horvath, S. (2019). DNA methylation GrimAge strongly predicts lifespan and healthspan. Aging, 11(2), 303–327. doi: 10.18632/aging.101684CrossRefGoogle ScholarPubMed
McCrory, C., Fiorito, G., Hernandez, B., Polidoro, S., O'Halloran, A. M., Hever, A., … Kenny, R. A. (2021). GrimAge outperforms other epigenetic clocks in the prediction of age-related clinical phenotypes and all-cause mortality. Journals of Gerontology. Series A: Biological Sciences and Medical Sciences, 76(5), 741–749. doi: 10.1093/gerona/glaa286CrossRefGoogle ScholarPubMed
McLean, S. A., Ressler, K., Koenen, K. C., Neylan, T., Germine, L., Jovanovic, T., … Kessler, R. (2020). The AURORA study: A longitudinal, multimodal library of brain biology and function after traumatic stress exposure. Molecular Psychiatry, 25(2), 283–296. doi: 10.1038/s41380-019-0581-3CrossRefGoogle ScholarPubMed
Mehta, D., Bruenig, D., Lawford, B., Harvey, W., Carrillo-Roa, T., Morris, C. P., … Voisey, J. (2018). Accelerated DNA methylation aging and increased resilience in veterans: The biological cost for soldiering on. Neurobiology of Stress, 8, 112–119. doi: 10.1016/j.ynstr.2018.04.001CrossRefGoogle Scholar
Mehta, D., Bruenig, D., Pierce, J., Sathyanarayanan, A., Stringfellow, R., Miller, O., … Shakespeare-Finch, J. (2022). Recalibrating the epigenetic clock after exposure to trauma: The role of risk and protective psychosocial factors. Journal of Psychiatric Research, 149, 374–381. doi: 10.1016/j.jpsychires.2021.11.026CrossRefGoogle ScholarPubMed
Mellon, S. H., Gautam, A., Hammamieh, R., Jett, M., & Wolkowitz, O. M. (2018). Metabolism, metabolomics, and inflammation in posttraumatic stress disorder. Biological Psychiatry, 83(10), 866–875. doi: 10.1016/j.biopsych.2018.02.007CrossRefGoogle ScholarPubMed
Morey, R. A., Clarke, E. K., Haswell, C. C., Phillips, R. D., Clausen, A. N., Mufford, M. S., … LaBar, K. S. (2020). Amygdala nuclei volume and shape in military veterans with posttraumatic stress disorder. Biological Psychiatry: Cognitive Neuroscience and Neuroimaging, 5(3), 281–290. doi: 10.1016/j.bpsc.2019.11.016Google ScholarPubMed
Morey, R. A., Gold, A. L., LaBar, K. S., Beall, S. K., Brown, V. M., Haswell, C. C., … McCarthy, G. (2012). Amygdala volume changes in posttraumatic stress disorder in a large case-controlled veterans group. Archives of General Psychiatry, 69(11), 1169–1178. doi: 10.1001/archgenpsychiatry.2012.50CrossRefGoogle Scholar
Na, P. J., Montalvo-Ortiz, J. L., Nagamatsu, S. T., Southwick, S. M., Krystal, J. H., Gelernter, J., … Pietrzak, R. H. (2022). Association of symptoms of posttraumatic stress disorder and GrimAge, an epigenetic marker of mortality risk, in US military veterans. Journal of Clinical Psychiatry, 83(4), 21br14309. doi: 10.4088/JCP.21br14309Google ScholarPubMed
Palma-Gudiel, H., Fañanás, L., Horvath, S., & Zannas, A. S. (2020). Psychosocial stress and epigenetic aging. International Review of Neurobiology, 150, 107–128. doi: 10.1016/bs.irn.2019.10.020CrossRefGoogle ScholarPubMed
Pruim, R. H. R., Mennes, M., Buitelaar, J. K., & Beckmann, C. F. (2015a). Evaluation of ICA-AROMA and alternative strategies for motion artifact removal in resting state fMRI. Neuroimage, 112, 278–287. doi: 10.1016/j.neuroimage.2015.02.063CrossRefGoogle ScholarPubMed
Pruim, R. H. R., Mennes, M., van Rooij, D., Llera, A., Buitelaar, J. K., & Beckmann, C. F. (2015b). ICA-AROMA: A robust ICA-based strategy for removing motion artifacts from fMRI data. Neuroimage, 112, 267–277. doi: 10.1016/j.neuroimage.2015.02.064CrossRefGoogle ScholarPubMed
Saygin, Z. M., Kliemann, D., Iglesias, J. E., van der Kouwe, A. J. W., Boyd, E., Reuter, M., … Augustinack, J. C. (2017). High-resolution magnetic resonance imaging reveals nuclei of the human amygdala: Manual segmentation to automatic atlas. Neuroimage, 155, 370–382. doi: 10.1016/j.neuroimage.2017.04.046CrossRefGoogle ScholarPubMed
Seo, J. H., Park, H. S., Park, S. S., Kim, C. J., Kim, D. H., & Kim, T. W. (2019). Physical exercise ameliorates psychiatric disorders and cognitive dysfunctions by hippocampal mitochondrial function and neuroplasticity in post-traumatic stress disorder. Experimental Neurology, 322, 113043. doi: 10.1016/j.expneurol.2019.113043CrossRefGoogle ScholarPubMed
Sheynin, J., Duval, E. R., King, A. P., Angstadt, M., Phan, K. L., Simon, N. M., … Liberzon, I. (2020). Associations between resting-state functional connectivity and treatment response in a randomized clinical trial for posttraumatic stress disorder. Depression and Anxiety, 37(10), 1037–1046. doi: 10.1002/da.23075CrossRefGoogle Scholar
Smith, A. K., Ratanatharathorn, A., Maihofer, A. X., Naviaux, R. K., Aiello, A. E., Amstadter, A. B., … Nievergelt, C. M. (2020). Epigenome-wide meta-analysis of PTSD across 10 military and civilian cohorts identifies methylation changes in AHRR. Nature Communications, 11(1), 5965. doi: 10.1038/s41467-020-19615-xCrossRefGoogle ScholarPubMed
Star, E. N., Zhu, M., Shi, Z., Liu, H., Pashmforoush, M., Sauve, Y., … Chow, R. L. (2012). Regulation of retinal interneuron subtype identity by the Iroquois homeobox gene Irx6. Development, 139(24), 4644–4655. doi: 10.1242/dev.081729CrossRefGoogle ScholarPubMed
Sumner, J. A., Colich, N. L., Uddin, M., Armstrong, D., & McLaughlin, K. A. (2019). Early experiences of threat, but not deprivation, are associated with accelerated biological aging in children and adolescents. Biological Psychiatry, 85(3), 268–278. doi: 10.1016/j.biopsych.2018.09.008CrossRefGoogle Scholar
Tian, Y., Morris, T. J., Webster, A. P., Yang, Z., Beck, S., Feber, A., & Teschendorff, A. E. (2017). ChAMP: Updated methylation analysis pipeline for Illumina BeadChips. Bioinformatics, 33(24), 3982–3984.CrossRefGoogle ScholarPubMed
U.S. Department of Health and Human Services (2013). The state of aging and health in America 2013. Atlanta, GA: Centers for Disease Control and Prevention.Google Scholar
Vaccarino, V., Goldberg, J., Rooks, C., Shah, A. J., Veledar, E., Faber, T. L., … Bremner, J. D. (2013). Post-traumatic stress disorder and incidence of coronary heart disease: A twin study. Journal of the American College of Cardiology, 62(11), 970–978. doi: 10.1016/j.jacc.2013.04.085CrossRefGoogle ScholarPubMed
Verhoeven, J. E., Yang, R., Wolkowitz, O. M., Bersani, F. S., Lindqvist, D., Mellon, S. H., … Hammamieh, R. (2018). Epigenetic age in male combat-exposed war veterans: Associations with posttraumatic stress disorder status. Molecular Neuropsychiatry, 4(2), 90–99. doi: 10.1159/000491431Google Scholar
Wang, Z., Hui, Q., Goldberg, J., Smith, N., Kaseer, B., Murrah, N., … Sun, Y. V. (2022). Association between posttraumatic stress disorder and epigenetic age acceleration in a sample of twins. Psychosomatic Medicine, 84(2), 151–158. doi: 10.1097/psy.0000000000001028CrossRefGoogle Scholar
Ware, J. Jr., Kosinski, M., & Keller, S. D. (1996). A 12-item short-form health survey: Construction of scales and preliminary tests of reliability and validity. Medical Care, 34(3), 220–233. doi: 10.1097/00005650-199603000-00003CrossRefGoogle ScholarPubMed
Weathers, F. W., Blake, D. D., Schnurr, P. P., Kaloupek, D. G., Marx, B. P., & Keane, T. M. (2013). The life events checklist for DSM-5 (LEC-5) – Extended [measurement instrument]. Retrieved from https://www.ptsd.va.gov/professional/assessment/te-measures/life_events_checklist.aspGoogle Scholar
Wolf, E. J., Logue, M. W., Morrison, F. G., Wilcox, E. S., Stone, A., Schichman, S. A., … Miller, M. W. (2019). Posttraumatic psychopathology and the pace of the epigenetic clock: A longitudinal investigation. Psychological Medicine, 49(5), 791–800. doi: 10.1017/s0033291718001411CrossRefGoogle Scholar
Wolf, E. J., Maniates, H., Nugent, N., Maihofer, A. X., Armstrong, D., Ratanatharathorn, A., … Logue, M. W. (2018). Traumatic stress and accelerated DNA methylation age: A meta-analysis. Psychoneuroendocrinology, 92, 123–134. doi: 10.1016/j.psyneuen.2017.12.007CrossRefGoogle ScholarPubMed
Yang, R., Wu, G. W. Y., Verhoeven, J. E., Gautam, A., Reus, V. I., Kang, J. I., … Wolkowitz, O. M. (2021). A DNA methylation clock associated with age-related illnesses and mortality is accelerated in men with combat PTSD. Molecular Psychiatry, 26(9), 4999–5009. doi: 10.1038/s41380-020-0755-zCrossRefGoogle ScholarPubMed
Yang, R. J., Mozhui, K., Karlsson, R. M., Cameron, H. A., Williams, R. W., & Holmes, A. (2008). Variation in mouse basolateral amygdala volume is associated with differences in stress reactivity and fear learning. Neuropsychopharmacology, 33(11), 2595–2604. doi: 10.1038/sj.npp.1301665CrossRefGoogle ScholarPubMed
Yeo, B. T., Krienen, F. M., Sepulcre, J., Sabuncu, M. R., Lashkari, D., Hollinshead, M., … Buckner, R. L. (2011). The organization of the human cerebral cortex estimated by intrinsic functional connectivity. Journal of Neurophysiology, 106(3), 1125–1165. doi: 10.1152/jn.00338.2011Google ScholarPubMed
Yong, W. S., Hsu, F. M., & Chen, P. Y. (2016). Profiling genome-wide DNA methylation. Epigenetics & Chromatin, 9, 26. doi: 10.1186/s13072-016-0075-3CrossRefGoogle ScholarPubMed
Zannas, A. S. (2019a). Decoding the life story of our epigenome. Epigenomics, 11(11), 1233–1236. doi: 10.2217/epi-2019-0155CrossRefGoogle ScholarPubMed
Zannas, A. S. (2019b). Epigenetics as a key link between psychosocial stress and aging: Concepts, evidence, mechanisms. Dialogues in Clinical Neuroscience, 21(4), 389–396. doi: 10.31887/DCNS.2019.21.4/azannasCrossRefGoogle ScholarPubMed
Zannas, A. S., Arloth, J., Carrillo-Roa, T., Iurato, S., Roh, S., Ressler, K. J., … Mehta, D. (2015a). Lifetime stress accelerates epigenetic aging in an urban, African American cohort: Relevance of glucocorticoid signaling. Genome Biology, 16, 266. doi: 10.1186/s13059-015-0828-5CrossRefGoogle Scholar
Zannas, A. S., Provençal, N., & Binder, E. B. (2015b). Epigenetics of posttraumatic stress disorder: Current evidence, challenges, and future directions. Biological Psychiatry, 78(5), 327–335. doi: 10.1016/j.biopsych.2015.04.003CrossRefGoogle ScholarPubMed
Zheng, Y., Joyce, B. T., Colicino, E., Liu, L., Zhang, W., Dai, Q., … Hou, L. (2016). Blood epigenetic age may predict cancer incidence and mortality. EBioMedicine, 5, 68–73. doi: 10.1016/j.ebiom.2016.02.008CrossRefGoogle ScholarPubMed
Ziobrowski, H. N., Kennedy, C. J., Ustun, B., House, S. L., Beaudoin, F. L., An, X., … van Rooij, S. J. H. (2021). Development and validation of a model to predict posttraumatic stress disorder and major depression after a motor vehicle collision. JAMA Psychiatry, 78(11), 1228–1237. doi: 10.1001/jamapsychiatry.2021.2427CrossRefGoogle Scholar
Figure 0

Table 1. Participant demographic and clinical variables

Figure 1

Fig. 1. Advanced GrimAge at emergency department (ED) presentation predicts increased risk for posttraumatic stress disorder (PTSD) 6 months later. (a) Proportion of participants with 6-month probable PTSD stratified by ED GrimAge tertile. (b) Graph depicting how ED GrimAge predicts the probability for 6-month PTSD diagnosis. Predicted probabilities have been adjusted for covariates (details in Methods section). Depicted GrimAge residuals are after adjusting for chronological age. Additional statistics are provided in Results.

Figure 2

Fig. 2. Example of a GrimAge CpG, the methylation levels of which predict 6-month posttraumatic stress disorder (PTSD) and related symptom trajectories at the nominal level of statistical significance. Predicted probabilities have been adjusted for covariates (as described in Methods). Detailed statistics are provided in online Supplementary Table S4.

Figure 3

Fig. 3. Advanced GrimAge is associated with reduced amygdala subregion volumes. Example is shown for volume of the accessory basal nucleus of the amygdala. The brain image on the right depicts the accessory basal nucleus from a single subject. Depicted GrimAge residuals are after adjusting for chronological age and DNA methylation-estimated blood cell proportions. Additional statistics are provided in Results.

Supplementary material: File

Zannas et al. supplementary material 1
Download undefined(File)
File 45.2 KB
Supplementary material: File

Zannas et al. supplementary material 2
Download undefined(File)
File 34.3 KB