Epigenome-wide association studies identify novel DNA methylation sites associated with PTSD: A meta-analysis of 23 military and civilian cohorts
=================================================================================================================================================

* Seyma Katrinli
* Agaz H Wani
* Adam X Maihofer
* Andrew Ratanatharathorn
* Nikolaos P Daskalakis
* Janitza Montalvo-Ortiz
* Diana L Núñez-Ríos
* Anthony S Zannas
* Xiang Zhao
* Allison E Aiello
* Allison E Ashley-Koch
* Diana Avetyan
* Dewleen G Baker
* Jean C Beckham
* Marco P Boks
* Leslie A Brick
* Evelyn Bromet
* Frances A Champagne
* Chia-Yen Chen
* Shareefa Dalvie
* Michelle F Dennis
* Segun Fatumo
* Catherine Fortier
* Sandro Galea
* Melanie E Garrett
* Elbert Geuze
* Gerald Grant
* Michael A Hauser
* Jasmeet P Hayes
* Sian MJ Hemmings
* Bertrand Russel Huber
* Aarti Jajoo
* Stefan Jansen
* Ronald C Kessler
* Nathan A Kimbrel
* Anthony P King
* Joel E Kleinman
* Nastassja Koen
* Karestan C Koenen
* Pei-Fen Kuan
* Israel Liberzon
* Sarah D Linnstaedt
* Adriana Lori
* Benjamin J Luft
* Jurjen J Luykx
* Christine E Marx
* Samuel A McLean
* Divya Mehta
* William Milberg
* Mark W Miller
* Mary S Mufford
* Clarisse Musanabaganwa
* Jean Mutabaruka
* Leon Mutesa
* Charles B Nemeroff
* Nicole R Nugent
* Holly K Orcutt
* Xue-Jun Qin
* Sheila A M Rauch
* Kerry J Ressler
* Victoria B Risbrough
* Eugène Rutembesa
* Bart P F Rutten
* Soraya Seedat
* Dan J Stein
* Murray B Stein
* Sylvanus Toikumo
* Robert J Ursano
* Annette Uwineza
* Mieke H Verfaellie
* Eric Vermetten
* Christiaan H Vinkers
* Erin B Ware
* Derek E Wildman
* Erika J Wolf
* Ross McD Young
* Ying Zhao
* Leigh L van den Heuvel
* PGC-PTSD Epigenetics Workgroup
* PsychENCODE PTSD Brainomics Project
* Traumatic Stress Brain Research Group
* Monica Uddin
* Caroline M Nievergelt
* Alicia K Smith
* Mark W Logue

## Abstract

**Background** The occurrence of post-traumatic stress disorder (PTSD) following a traumatic event is associated with biological differences that can represent the susceptibility to PTSD, the impact of trauma, or the sequelae of PTSD itself. These effects include differences in DNA methylation (DNAm), an important form of epigenetic gene regulation, at multiple CpG loci across the genome. Moreover, these effects can be shared or specific to both central and peripheral tissues. Here, we aim to identify blood DNAm differences associated with PTSD and characterize the underlying biological mechanisms by examining the extent to which they mirror associations across multiple brain regions.

**Methods** As the Psychiatric Genomics Consortium (PGC) PTSD Epigenetics Workgroup, we conducted the largest cross-sectional meta-analysis of epigenome-wide association studies (EWASs) of PTSD to date, involving 5077 participants (2156 PTSD cases and 2921 trauma-exposed controls) from 23 civilian and military studies. PTSD diagnosis assessments were harmonized following the standardized guidelines established by the PGC-PTSD Workgroup. DNAm was assayed from blood using either Illumina HumanMethylation450 or MethylationEPIC (850K) BeadChips. A common QC pipeline was applied. Within each cohort, DNA methylation was regressed on PTSD, sex (if applicable), age, blood cell proportions, and ancestry. An inverse variance-weighted meta-analysis was performed. We conducted replication analyses in tissue from multiple brain regions, neuronal nuclei, and a cellular model of prolonged stress.

**Results** We identified 11 CpG sites associated with PTSD in the overall meta-analysis (1.44e-09 < *p* < 5.30e-08), as well as 14 associated in analyses of specific strata (military vs civilian cohort, sex, and ancestry), including CpGs in *AHRR* and *CDC42BPB*. Many of these loci exhibit blood-brain correlation in methylation levels and cross-tissue associations with PTSD in multiple brain regions. Methylation at most CpGs correlated with their annotated gene expression levels.

**Conclusions** This study identifies 11 PTSD-associated CpGs, also leverages data from postmortem brain samples, GWAS, and genome-wide expression data to interpret the biology underlying these associations and prioritize genes whose regulation differs in those with PTSD.

Keywords
*   PTSD
*   trauma
*   DNA methylation
*   postmortem brain
*   GWAS
*   gene expression

## Introduction

Posttraumatic stress disorder (PTSD) is a serious psychiatric disorder characterized by intrusive memories of the traumatic event(s), avoidance of or numbing to situations that trigger those memories, and hyperarousal symptoms that can disturb mental and physical health (1).

These symptoms are associated with lower levels of self-care, lower compliance with medical treatment, and higher rates of substance use (2, 3). Thus, it is not surprising that PTSD increases the risk for chronic medical conditions, such as cardiovascular disorders, independent of lifestyle factors (e.g., substance use and sleep quality)(4, 5). Although most individuals experience at least one traumatic event, only a small fraction develop PTSD(6). Genetic and environmental factors contribute to this differential susceptibility in PTSD development upon trauma exposure (7, 8).

Genome-wide association studies (GWAS) of PTSD demonstrated remarkable success at identifying relevant genes, many of which are involved in the stress response or immune function (see reviews(9, 10)). The recent Psychiatric Genomics Consortium PTSD Workgroup (PGC-PTSD) Freeze 3 GWAS identified 95 genomic loci associated with PTSD, implicating genes involved in stress, immune, fear, and threat-related processes (11). Nonetheless, genetic differences do not fully account for an individual’s susceptibility to PTSD. Trauma exposure has been shown to alter epigenetic patterns in both animal and human studies, prompting the need to conduct epigenetic studies of PTSD in addition to genetic studies (12, 13). Epigenetic mechanisms are chemical modifications that can dictate the timing and magnitude of gene expression without altering the DNA sequence (14). The most widely studied epigenetic mechanism is DNA methylation (DNAm), which is defined as the addition of a methyl group to cytosine bases, particularly at cytosine-guanine dinucleotides (CpG sites). DNAm patterns respond to changes in the environment, are potentially reversible, and can be targeted for disease therapies (15, 16). Environmental influences on DNAm are apparent across the life span and may provide insight into the biological response to trauma(17).

Which specific DNAm sites differ across individuals and how they correlate with exposures and gene expression can vary across tissues (18). DNAm in human brain tissue, which is most relevant to the study of psychiatric disorders, is not easily accessible in living patients and hence is not a viable PTSD biomarker for clinical use. However, correlation has been observed between peripheral tissues (e.g., blood) and brain DNAm levels at specific genomic loci, and hence blood DNAm can potentially serve as a robust biomarker for implementing early intervention and developing improved preventative or therapeutic strategies for PTSD (19, 20).

Moreover, PTSD symptoms have been linked to the components of the peripheral immune system (21, 22) that can be readily assessed in blood DNAm. Multiple peripheral epigenome-wide association studies (EWASs) of PTSD identified CpGs in genes related to the immune system and neurotransmission(23–28). While prior EWASs of PTSD have reported promising results, the small sample sizes and variability of analysis methods across studies make it difficult to combine and interpret the findings effectively. Recent meta-analyses led by the PGC-PTSD Epigenetics Workgroup minimized these issues by increasing sample size, increasing sample diversity, and using a common quality control and analysis pipeline(29–33). These meta-analyses identified multiple new loci associated with PTSD, including *NRG1, AHRR, MAD1L1*, and *TBXAS1*, implicating immune dysregulation in those with PTSD (30–33).

Building on the prior work by Smith et al. (31), which conducted an EWAS meta-analysis in 1896 participants from 10 cohorts, this study includes 13 additional cohorts with a denser and more comprehensive DNAm array, bringing the sample up to 5077 participants from 23 civilian and military cohorts. Our current investigation replicated the findings of the initial PGC-PTSD epigenome-wide meta-analysis, reporting lower *AHRR* methylation in those with PTSD, and identified 8 new (11 total) PTSD-associated loci. We also leveraged data from postmortem brain samples, a cellular model of prolonged stress, GWAS, and genome-wide gene expression studies to interpret the biology underlying these associations and prioritize genes whose regulation differs in those with PTSD.

## Methods

### Cohorts and post-traumatic stress disorder assessments

The study includes 2156 current PTSD cases and 2921 trauma-exposed controls from 9 civilian cohorts: BEAR, DNHS, DCHS, GTP, NIU, Shared Roots, AURORA, H3A_Rwanda, WTC; and 9 military cohorts: GMRFQUT, MRS, PRISMO, Army STARRS, PROGrESS, NCPTSD/TRACTS, INTRuST, and VA cohorts (VA-M-AA & VA-M-EA). For DNHS, GTP, MRS, PRISMO, and Army STARRS, two different datasets were available based on the DNAm array. Two different datasets for these five cohorts did not have any overlapping samples and were treated as independent studies. Sample characteristics for the 23 studies that participated in the meta-analysis are summarized in Table 1.

View this table:
[Table 1.](http://medrxiv.org/content/early/2024/07/15/2024.07.15.24310422/T1)

Table 1. Overview of the studies

The sample is heterogeneous in terms of sex (41% female), ancestry (46% European, 43% African, and 11% of other ancestries), and cohort type (56% military cohort). Civilian cohorts skew towards being more female (67%) and African (70%), whereas military cohorts are predominantly male (80%) and European (64%).

All participants were exposed to a traumatic event and 42% met the criteria for current PTSD. The current PTSD diagnosis was assessed by each study following the standardized guidelines established by the PGC-PTSD Workgroup (7). Briefly, current PTSD diagnosis was determined based on the specific criteria set by the principal investigator of each study. All control subjects had experienced trauma. Subjects without a current PTSD but with a prior history of PTSD (i.e., remitted PTSD), were excluded. Detailed descriptions of cohorts and PTSD assessments are provided in the Supplement. All participants in these studies gave informed consent. The institutional review boards of each respective institution approved these studies.

### DNA methylation

Whole blood DNAm was measured using the Illumina MethylationEPIC BeadChip in 14 studies, and the HumanMethylation450 BeadChip in 9 studies (Table 1). All studies used a standardized consortium-developed quality control (QC) pipeline that differed somewhat depending on which chip was used. The 450K array pipeline (29, 31) is described in eMethods in the Supplement.

The EPIC pipeline (available at [https://github.com/PGC-PTSD-EWAS/EPIC_QC](https://github.com/PGC-PTSD-EWAS/EPIC_QC)) was similar to the 450K array pipeline. Samples with probe detection call rates lower than 90% and average intensity values that were either less than 50% of the overall sample mean or below 2000 arbitrary units (AU) were excluded. Probes with detection *p*-values > 0.01 were considered low quality and treated as missing. Probes that were missing in > 10% of the samples within the studies and were cross-hybridizing were removed (34). Data was normalized using single-sample Noob (ssNoob) implemented in R package *minfi* (*35*). *ComBat* was used to account for batch effects of chip and position while preserving PTSD, age, and sex effects (if applicable) (36). Blood-cell composition (i.e., the proportion of CD8+T, CD4+T, natural killer (NK), B cells, monocytes, and neutrophils) was estimated using the Robust Partial Correlation (RPC) method in *Epidish* (37) with a reference data specific to EPIC array (38). For studies without genome-wide genotype data (DNHS, NIU, Shared Roots, AURORA, H3A_Rwanda, GMRFQUT, PROGrESS), we estimated ancestry principal components (PCs) from DNAm data, using the method developed by Barfield et al. (39), as previously described (33). PCs 2 and 3, which were the components that correlate most with self-reported ancestry, were included as covariates (33, 39). In cohorts with available genome-wide genotype data, PCs 1-3 from GWAS were used to adjust for ancestry. We used R package *bacon* to control inflation, only if doing so results in the lambda being closer to 1 (40). To predict smoking status, a DNAm-based smoking score was calculated, as previously described (27) for cohorts with EPIC array data. A detailed description of DNAm-based ancestry PCs and smoking score calculation is provided in eMethods in the Supplement.

### Epigenome-wide association analysis

The association between PTSD and DNAm was tested using multivariable linear regression models within cohorts with balanced plate designs. For studies in which plate layouts were not balanced (Shared Roots, H3A_Rwanda, GMRFQUT), we conducted mixed-effect regression models, including chip as a random effect term. R package *CpGassoc* was used to fit the models (41). The models were adjusted for age, sex (if applicable), blood cell composition (i.e., CD8T, CD4T, NK, B cell, and monocyte cell proportions), and ancestry PCs. A post-hoc sensitivity analysis was performed by including a covariate for smoking: DNAm-based smoking score in studies with EPIC data and current smoking status for studies with 450K data. Furthermore, we conducted stratified analyses for both sexes, ancestry (European and African ancestry), and cohort type (civilian and military cohorts).

To combine results across studies, we performed inverse-variance weighting (IVW) meta-analysis in *meta* (42). Meta-analysis tested 411,786 CpGs common to 450K and EPIC arrays (23 studies), and 404,794 EPIC array specific CpGs (14 studies). Epigenome-wide significance threshold recommended for the MethylationEPIC BeadChip (*p* < 9.0E-08) was used to determine statistical significance (43). Gene Ontology (GO) enrichment analyses were conducted using the top 1000 CpGs in *missMethyl* (44). An FDR threshold of 5% was used to identify significant GO terms.

### Cross-tissue association analyses

#### Blood – Brain Correlations

The Blood Brain DNA Methylation Comparison Tool (19) was used to assess the correlations between methylation in blood and prefrontal cortex (PFC), entorhinal cortex (EC), superior temporal gyrus (STG), and cerebellum.

#### Postmortem brain DNAm

DNAm measured from post-mortem brains was obtained from two studies, each of which examined a unique but not necessarily distinct set of brain regions and cohorts: the National Center for PTSD Brain Bank cohort (NCPTSD-BB (45)) and the PsychENCODE Consortium for PTSD (PEC-PTSD) Brainomics cohort (46) (see eMethods in the Supplement for details), both of which were sourced from the Lieber Institute for Brain Development.

Methylation at PTSD-associated CpGs from the EWAS was tested for association with PTSD in DNA extracted from postmortem dorsolateral prefrontal cortex (dlPFC, BA9/46), ventromedial prefrontal cortex (vmPFC, BA12/32), amygdala, and dentate gyrus (DG). DNAm in the post-mortem tissue was measured using the EPIC array. We examined the associations with PTSD in dlPFC and vmPFC of 42 PTSD cases and 30 controls from the NCPTSD-BB. The associations between DNAm and PTSD in amygdala and DG were tested in 77 PTSD cases and 77 controls from the PEC-PTSD.

#### Neuronal nuclei

We examined cross-tissue association from neuronal nuclei isolation from the orbitofrontal cortex (OFC) of 25 PTSD cases and 13 healthy controls collected at the VA’s NCPTSD-BB (47). Fluorescence-Activated Nuclei Sorting (FANS) protocol was employed to isolate NeuN+ cells and the nuclei underwent reduced representation oxidative bisulfite-sequencing (RRoxBS), as previously described (48). We examined whether there was differential methylation and hydroxymethylation within 500 bp of CpGs from the epigenome-wide association analyses. Eight CpG sites match between the EPIC array and RRoxBS, two evaluating the same position (cg21566642 and cg25691167), and six within 500 bp (cg05575921, cg15148933, cg19558029, cg21161138, cg23576855, and cg26599989).

#### Cellular model of prolonged stress

We explored the associations between each of the PTSD-associated CpGs resulting from epigenome-wide association analyses and a cellular model of prolonged stress in which fibroblasts were subjected to physiological stress hormone (cortisol) levels for a prolonged period (51 days) as previously described (49, 50).

### Gene Regulation

Correlations between PTSD-associated CpG sites’ DNAm levels and expression levels of the corresponding gene (as determined by the EPIC array annotation) were tested in whole-blood RNA-sequencing (RNA-seq) data from participants in the BEAR (n=127), AURORA (n=173), NCPTSD Merit (n=204), and MRS cohorts (n=128 with multiple visits totaling 357 samples). The results were meta-analyzed using the IVW method using a Bonferroni correction for the number of CpGs examined. Detailed information about cohort-specific RNA-seq data generation is described in eMethods in the Supplement.

### Genetic effects

To evaluate the effect of nearby (<1 MB) polymorphisms on DNAm levels of CpGs associated with PTSD, we used cis-methylation quantitative trait locus (cis-meQTL) data from GoDMC(51) and meQTL EPIC(52) databases. For both databases, their default multiple testing adjustment was utilized: an FDR threshold of 5% in meQTL EPIC and *p* < 1e-08 in GoDMC. We checked the associations between the identified cis-SNPs and PTSD in the recent Freeze 3 GWAS from PGC-PTSD(11). Finally, we evaluated genetic associations between the genes with PTSD-associated DNAm changes and PTSD, using the gene-based test results from the recent PGC-PTSD Freeze 3 GWAS(11) using a Bonferroni correction for the 11 genes examined.

## Results

### Epigenome-wide association meta-analysis

We identified 11 PTSD-associated CpGs that passed the epigenome-wide significance threshold (*p* < 9e-08, Table 2, Figure 1, eFigure 1 in the Supplement). All CpG sites, except one site (cg21161138) near *AHRR*, remained nominally significant (1.52e-03 < *p* < 2.10e-07) with the same direction of association in the sensitivity analysis adjusted for smoking score (eTable 1 and eFigure 1 in the Supplement).

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/15/2024.07.15.24310422/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2024/07/15/2024.07.15.24310422/F1)

Figure 1. Manhattan Plot of the epigenome-wide association meta-analyses.
The x-axis depicts chromosomes and the location of each CpG site across the genome. The y-axis depicts the -log10 of the unadjusted *p*-value for the association with current PTSD. Each dot represents a CpG site. The solid blue line indicates the epigenome-wide statistical significance at *p* < 9.0e-8.

View this table:
[Table 2.](http://medrxiv.org/content/early/2024/07/15/2024.07.15.24310422/T2)

Table 2. CpG sites associated with current PTSD in the primary meta-analysis

### Stratified analyses for sex, ancestry, and cohort type

The directions of associations of the 11 significant CpGs were consistent at nominal significance (p < 0.05) when stratified by both sex, cohort type (civilian and military), and European ancestry. However, in the African ancestry stratified analysis, only 8 CpGs remained associated (p < 0.05) with PTSD (eTable 2 in the Supplement). We identified 1 epigenome-wide significant CpG site associated with PTSD in the female-stratified analysis, 1 in the male-stratified analysis, 5 in European ancestry-stratified analysis, 1 in the African ancestry-stratified analysis, 1 in the analysis for civilian-cohort analysis, and 5 in the military-cohort analysis (Table 3, eFigures 2-8 in the Supplement). Of note, some CpGs identified in the stratified analyses were specific to their stratum and were not associated with PTSD (*p* > 0.05) in the other stratum. For instance, cg25691167 in *FERD3L* was associated with PTSD in females (*p* = 4.24E-08), but not in males (*p* = 0.56) or the primarily male military cohorts (*p* = 0.13).

![Fig 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/15/2024.07.15.24310422/F2.medium.gif)

[Fig 2.](http://medrxiv.org/content/early/2024/07/15/2024.07.15.24310422/F2)

Fig 2. Summary of all analyses and findings.
Figure combines CpGs from the main analysis (gold) and stratified analyses for sex, ancestry, and trauma type (light gold); and summarizes the results of blood and brain correlations (rose); cross-tissue associations for multiple brain regions (light blue), neuronal nuclei (blue), and a fibroblast model of prolonged stress (aqua), gene expression (purple); and genetic effects, including methylation quantitative trait loci (meQTL) analyses (light green) and genetic associations from the recent PGC-PTSD GWAS (dark green). Positive findings (*p* < 0.05) are indicated with the specific color of the respective category. Asterisk (*) indicates epigenome-wide significance (*p* < 9e-8). Gray represents the CpGs or genes were not present in the respective datasets. PFC: prefrontal cortex. EC: entorhinal cortex. STG: superior temporal gyrus. CER: cerebellum. dlPFC: dorsolateral prefrontal cortex. vmPFC: ventromedial prefrontal cortex. DG: dentate gyrus. 5mC: 5-Methylcytosine. 5hmC: 5-Hydroxymethylcytosine. GC: glucocorticoid.

View this table:
[Table 3.](http://medrxiv.org/content/early/2024/07/15/2024.07.15.24310422/T3)

Table 3. CpG sites associated with current PTSD in the stratified analyses

### Gene enrichment analysis

We identified 3 GO term enrichments in the female-stratified analysis (FDR < 0.05), including nervous system development (eTable 3 in the Supplement). We did not identify any significant enrichments in the main analysis or other strata.

### Cross-tissue associations

Of the 19 significant CpGs from the primary and stratified analyses, 3 were not available in the Blood Brain DNA Methylation Comparison Tool database. Methylation at 9 of the remaining 16 CpGs was correlated between blood and at least one brain region (Figure 2, eTable4 in the Supplement). The strongest correlation was observed between the blood and PFC (r = 0.91) for cg23576855 (*AHRR*).

We next evaluated whether the PTSD-associated CpGs from blood were also associated with PTSD in the brain regions of dlPFC, vmPFC, amygdala, and DG, as well as neuronal nuclei. Out of 19 CpGs from the primary and stratified analyses in bulk tissue, 9 were nominally associated with PTSD in at least one brain region or the neuronal nuclei (*p* < 0.05, Figure 2, eTables 5 and 6 in the Supplement). For instance, PTSD cases had higher cg04987734 (*CDC42BPB*) methylation both in the blood (*p* = 3.26e-8) and the dlPFC (*p* = 3.9e-3). In addition, 7 CpGs exhibited nominally significant methylation changes (5 CpGs with the same direction of association) in fibroblasts when subjected to physiological stress hormone levels in a cellular model of prolonged stress (0.024 < ΔDNAm *<* 0.225*, p* < 0.05, Figure 2, eTable 7 in the Supplement)(49). Notably, the CpG sites (cg21161138 [*AHRR*], cg25320328 [*GFI*], cg19719391 [*intergenic*], cg25691167 [*FERD3L*]) that were consistently associated with PTSD in blood and postmortem brain also demonstrated a significant change in this cell model.

### Gene Regulation

Of the 11 PTSD-associated CpGs annotated to a gene, 10 correlated with their gene expression in at least one of the 4 RNA-seq cohorts (Figure 2, eTable 8 in the Supplement). Meta-analysis across cohorts identified 4 CpGs whose methylation levels were negatively correlated with their annotated gene expression after multiple test corrections: cg05575921, cg21161138, and cg23576855 with *AHRR*, and cg27541344 with *BCL11B* (p < 4.5e-03, -16.69 ≤ z ≤ -4.86 eTable 8 in the Supplement).

### Genetic regulation of PTSD-associated CpGs

Of the 19 CpGs from the primary and stratified analyses, 12 CpGs were associated with at least one nearby SNP within 1 MB according to GoDMC(51) and meQTL EPIC(53) databases (eTable 9 in the Supplement). In total, we identified 41 lead meQTLs, of which 8 were nominally associated with PTSD (*p* < 0.05, albeit in the opposite direction) in the PGC-PTSD Freeze 3 GWAS(11) (eTable 9 in the Supplement). For instance, lower methylation at the intergenic cg14753356 site is associated with PTSD (z = -5.72, *p* = 1.09E-08, Table 2) and higher methylation at cg14753356 is associated with rs28986310 T allele (beta = 0.31, *p* < 5e-324), which increases the risk of PTSD (z = 6.73, *p* = 1.68e-11, eTable 9 in the Supplement).

In addition, of the 11 genes harboring PTSD-associated CpGs, *AHRR*, *CDC42BPB*, and *BCL11B* were implicated in the recent PGC-PTSD gene-based analysis (3.48e-03 < *p* < 5.40e-05, eTable 10 in the Supplement). All genes are still significant after Bonferroni correction for 11 tests (p < 4.5e-03).

## Discussion

In the present epigenome-wide meta-analysis of blood DNAm levels, we identified 11 CpG sites associated with PTSD. Many of these are also associated with PTSD in multiple brain regions and a cellular model of prolonged stress exposure (49). Many of the PTSD-associated CpGs were associated with the expression of their respective genes. An examination of the most recent PGC-PTSD GWAS (11) indicated that many SNPs near significant CpGs were associated with both DNAm levels and PTSD diagnosis. Moreover, many of the EWAS-implicated genes were significant in gene-based tests from the GWAS. Thus, the GWAS results represent further evidence supporting the role of these genes in PTSD pathogenesis. However, the lower significance of the SNPs in association with PTSD relative to the corresponding CpGs (eTable 10 in the Supplement) indicates that the association between PTSD and the CpGs is likely not simply a byproduct of PTSD-associated SNPs.

Our *CDC42BPB* (CDC42 binding protein kinase beta) findings are of particular interest. *CDC42BPB* is involved in the regulation of cytoskeletal rearrangement, cell migration, and neurodevelopment(54). In the current study, increased *CDC42BPB* methylation at cg04987734 was associated with PTSD in blood and the dlPFC. Methylation at this site was also positively correlated with *CDC42BPB* expression in blood. Notably, higher methylation at cg04987734 has been associated with depressive symptoms (55) and increased C-reactive protein (CRP) levels (56, 57), which is perhaps not surprising given the bi-directional genetic association between PTSD and CRP levels (58). Multiple studies reported increased CRP levels and other inflammatory markers in those with PTSD, suggesting inflammation as an important component of PTSD (59–61). Future studies are warranted to investigate *CDC42BPB* methylation as a potential mediator of the relationship between PTSD and inflammation.

We observed 3 CpGs (cg05575921, cg21161138, and cg23576855) in *AHRR* (aryl-hydrocarbon receptor repressor) that were identified in the initial PGC-PTSD EWAS (31) and an independent study of US veterans (45). The aryl hydrocarbon receptor (AhR) plays a role in immunomodulation, including the regulation of T lymphocytes, B cell maturation, and the activity of macrophages, dendritic cells, and neutrophils (62), supporting the link between the immune system and PTSD. Although methylation at the *AHRR* CpGs is known to be influenced by smoking (63), cg05575921 and cg23576855 remained nominally significant (*p* = 1.52E-03 and p = 7.59E-04, respectively) in our sensitivity analysis adjusting for DNAm-based smoking scores. Notably, cg21161138 DNAm was also associated with PTSD in postmortem dlPFC and neuronal nuclei and demonstrated a significant change in the cellular model of prolonged stress (49), supporting the notion that there is an association between methylation at this locus and PTSD which is independent of the effects of smoking.

The stratified analyses identified DNAm-PTSD associations specific to sex, ancestry, and cohort type. The PTSD-associated site cg25691167 (*FERD3L*) in females (*p* = 4.24e-08) was not associated with PTSD in males (*p* = 0.57), suggesting that DNAm changes in cg25691167 might be sex-specific. Similarly, cg27541344 (*BCL11B*) was associated with PTSD in the civilian (*p* = 7.21e-08), but not the military cohorts (*p* = 0.38), whereas 3 out of 5 PTSD-associated CpGs in the military cohorts were not significant in the civilian cohorts (p > 0.05). We speculate that some of these DNAm changes specific to sex, ancestry, or cohort type might be due to unique characteristics of those strata, such as hormonal factors and environmental exposures.

## Strengths and Limitations

To our knowledge, this is the largest EWAS of PTSD to date. Our sample is diverse in terms of sex, ancestry, and cohort type. We leveraged data from postmortem brain samples, a cellular model of prolonged stress, GWAS, and genome-wide expression data to support our findings.

However, the study is not without limitations. First, methylation BeadChips only assess a subset of CpG sites in the genome; therefore, we may not capture all PTSD-associated CpG sites. Second, this is a cross-sectional study of participants with prior exposure to a traumatic event; thus, we were not able to assess whether the differences in DNAm between individuals with and without PTSD are a cause or consequence of PTSD or both. Third, our primary meta-analysis was performed using measures of blood DNAm. While this strategy provided valuable insights for future research on biomarkers of PTSD, it might not accurately represent the DNAm patterns within other tissues that are likely the most relevant to PTSD. However, the majority of the PTSD-associated CpG sites’ methylation levels were correlated between blood and at least one brain region. In addition, most PTSD-associated CpGs in blood were also associated with PTSD in one or more brain regions. Fourth, we used bulk tissue and adjusted for cellular heterogeneity, which might have obscured some signals, given the alterations in cell composition in those with PTSD(64, 65). Additionally, the brain regions examined varied between the online databases used to examine the blood-brain correlation of methylation values and the gene expression differences associated with PTSD, which can complicate interpretation. Finally, most cohorts that participated in the meta-analysis did not have detailed physical or psychiatric information on participants, including detailed information on trauma type and timing, PTSD symptom course, and treatment, making it challenging to evaluate and adjust for potential confounders, including substance use, comorbidities, or medication use.

## Conclusions

Taken together, this study replicates our previous findings and identifies novel PTSD-associated CpGs. Supporting data from multiple sources suggest that epigenetic mechanisms, particularly methylation in *AHRR* and *CDC42BPB*, may contribute to the complex relationship between the immune system and PTSD.

## Supporting information

Supplemental Files [[supplements/310422_file02.pdf]](pending:yes)

## Data Availability

The main summary statistics data that support the findings of this study will be available within Supplementary Data upon publication. Individual-level data from the cohorts or cohort-level summary statistics will be made available to researchers following an approved analysis proposal through the PGC-PTSD Epigenetics Workgroup with agreement of the cohort PIs. The raw data for the GTP cohort is available in the Gene Expression Omnibus database with the accession code GSE132203. Owing to military cohort data sharing restrictions, data from the VA, VA MIRECC, MRS, Army STARRS, PRISMO, PROGrESS, and NCPTSD/TRACTS cannot be publicly posted. However, such data can be provided in de-identified from a data repository through a data use agreement following applicable guidelines on data sharing and privacy protection. For additional information on access to these data, including PI contact information for the contributing cohorts, please contact the corresponding author.

## List of abbreviations

CpG
:   Cytosine-guanine dinucleotides
DG
:   Dentate gyrus
dlPFC
:   Dorsolateral prefrontal cortex
DNAm
:   DNA methylation
EWAS
:   Epigenome-wide association study
FANS
:   Fluorescence-Activated Nuclei Sorting
GWAS
:   Genome-wide association study
meQTL
:   Methylation quantitative trait locus
NK
:   Natural killer
OFC
:   Orbitofrontal cortex
PC
:   Principal component
PGC-PTSD
:   Psychiatric Genomics Consortium PTSD Workgroup
PTSD
:   Posttraumatic stress disorder
RPC
:   Robust Partial Correlation
RRoxBS
:   Reduced representation oxidative bisulfite-sequencing
vmPFC
:   Ventromedial prefrontal cortex

## Declarations

### Ethics approval and consent to participate

All participants in these studies gave informed consent. The institutional review boards of each respective institution approved these studies.

### Consent for publication

Not applicable.

### Availability of data and materials

The main summary statistics data that support the findings of this study will be available within Supplementary Data upon publication. Individual-level data from the cohorts or cohort-level summary statistics will be made available to researchers following an approved analysis proposal through the PGC-PTSD Epigenetics Workgroup with agreement of the cohort PIs. The raw data for the GTP cohort is available in the Gene Expression Omnibus database with the accession code GSE132203. Owing to military cohort data sharing restrictions, data from the VA, VA MIRECC, MRS, Army STARRS, PRISMO, PROGrESS, and NCPTSD/TRACTS cannot be publicly posted. However, such data can be provided in de-identified from a data repository through a data use agreement following applicable guidelines on data sharing and privacy protection. For additional information on access to these data, including PI contact information for the contributing cohorts, please contact the corresponding author.

### Competing interests

CYC is an employee of Biogen. NPD has served on scientific advisory boards for BioVie Pharma, Circular Genomics and Sentio Solutions for unrelated work. NRN serves as an unpaid member of the Ilumivu advisory board. SAMR receives support from the Wounded Warrior Project (WWP), Department of Veterans Affairs (VA), National Institute of Health (NIH), McCormick Foundation, Tonix Pharmaceuticals, Woodruff Foundation, and Department of Defense (DOD). Dr. Rauch receives royalties from Oxford University Press and American. KJR serves as a consultant for Acer, Bionomics, and Jazz Pharma; SABs for Sage, Boehringer Ingelheim, and Senseye. DJS has received consultancy honoraria from Discovery Vitality, Johnson & Johnson, Kanna, L’Oreal, Lundbeck, Orion, Sanofi, Servier, Takeda and Vistagen. MBS has in the past 3 years received consulting income from Acadia Pharmaceuticals, Aptinyx, atai Life Sciences, BigHealth, Biogen, Bionomics, BioXcel Therapeutics, Boehringer Ingelheim, Clexio, Delix Therapeutics, Eisai, EmpowerPharm, Engrail Therapeutics, Janssen, Jazz Pharmaceuticals, NeuroTrauma Sciences, PureTech Health, Sage Therapeutics, Sumitomo Pharma, and Roche/Genentech. MBS has stock options in Oxeia Biopharmaceuticals and EpiVario. MBS has been paid for his editorial work on Depression and Anxiety (Editor-in-Chief), Biological Psychiatry (Deputy Editor), and UpToDate (Co-Editor-in-Chief for Psychiatry). MBS has also received research support from NIH, Department of Veterans Affairs, and the Department of Defense. MBS is on the scientific advisory board for the Brain and Behavior Research Foundation and the Anxiety and Depression Association of America.

### Funding

This work was supported by the National Institute of Mental Health (NIMH; R01MH108826 and R01MH106595). This work was also supported by I01BX003477, a Department of Veterans Affairs BLR&D grant to MWL; 1R03AG051877, 1R21AG061367-01, RF1AG068121, and 1I01CX001276-01A2 to EJW; R21MH102834, 5I01CX000431, 5R01MH079806 to MWM; R01MD011728 to MU; R01MH105379 to NRN; R01MH117291, R01MH117292, R01MH117293 to FAC; R01MH108826 to NPD; U01MH115485 to LM, SF, SJ, JM, AU, and DEW; R01MH117291, R01MH117292, R01MH117293 to JEK; CDC/NIOSH U01OH012466 to PFK; R01MH117291, R01MH117292, R01MH117293 to CBN; R01MH093500 to CMN; 1R15MH099521-01, 5R21MH085436-02, 1R15HD049907-01A1 to HKO; R01MH108826 to AKS; I01BX002577, IK2CX000525, and lK6BX003777 to NAK; U01MH110925 to SAM; R21DA050160 and 1DP1DA058737 to JMO; Department of Veterans Affairs B9254-C to WM, B3001-C to CF; Department of Defense W81XWH-11-1-0073 and the National Center for Advancing Translational Sciences of the NIH UL1TR000433PEC-PTSD to SAMR; BX006186; BX005872 to VBR; VIDI award (91718336) from the Netherlands Scientific Organization to BPFR; Dutch Research Council (NWO) VIDI grant (09150171910042) to CHV. Brainomics work was supported by R01MH117291, R01MH117292, and R01MH117293. DJS and NK was supported by the South African Medical Research Council and the Bill and Melinda Gates Foundation (OPP 1017641). SS is supported by the South African Medical Research Council. SK is supported by the Building Interdisciplinary Research Careers in Women’s Health of the National Institutes of Health under Award Number K12HD085850.

### Author Contributions

PGC-PTSD writing group: S.K., M.W.L., C.M.N., A.K.S., and M.U.

Study PI or co-PI: A.E.A., A.E.A.-K., D.G.B., J.C.B., E.B., C.F., S.G., E.G., G.G., M.A.H., S.J., R.C.K., N.A.K., J.E.K., K.C.K., P.-F.K., M.W.L., B.L., C.E.M., S.A.M., W.M., M.W.M., L.M., C.B.N., C.M.N., N.R.N., H.K.O., S.A.M.R., K.J.R., V.R., S.S., A.K.S., D.J.S., M.B.S., M.U., R.J.U., E.V., D.E.W., E.J.W., and R.M.Y

Obtained funding for studies: J.C.B., M.P.B., S.F., C.F., E.G., M.A.H., R.C.K., K.C.K., M.W.L., J.J.L., S.A.M., M.W.M., C.M.N., N.R.N., H.K.O., S.A.M.R., K.J.R., B.P.F.R., A.K.S., M.U., R.J.U., and E.V.

Clinical: D.G.B., J.C.B., M.F.D., S.F., C.F., E.G., J.P.H., N.A.K., N.K., C.M., J.M., S.A.M.R., K.J.R., E.R., A.U., M.H.V., E.V., E.J.W., and L.L.V.D.H.

Contributed data: A.E.A.-K., D.G.B., J.C.B., F.A.C., N.P.D., M.F.D., C.F., M.A.H., J.P.H., S.M.J.H., B.R.H., S.J., S.K., N.A.K., A.P.K., K.C.K., P.-F.K., I.L., A.L, B.L., J.J.L., D.M., W.M., M.W.M., J.M.-O., C.M., L.M., C.M.N., N.R.N., S.A.M.R., V.R., S.S., A.K.S., M.B.S., S.T., M.U., A.U., M.H.V., D.E.W., E.J.W., A.S.Z., and L.L.V.D.H.

Statistical analysis: A.E.A.-K., D.A., M.P.B., L.B., C.-Y.C., S.D., N.P.D., S.F., M.E.G., A.J., S.K., A.P.K., I.L., M.W.L., A.X.M., D.M., M.S.M., C.M.N., D.L.N.-R., X-J.Q., A.R., B.P.F.R., A.K.S., C.H.V., A.H.W., E.B.W., A.S.Z., X.Z., and Y.Z.

Bioinformatics: M.P.B., L.B., C.-Y.C., N.P.D., M.E.G., M.A.H., A.J., S.K., S.D.L., M.W.L., A.X.M., A.R., B.P.F.R., A.H.W., X.Z., and Y.Z.

Genomics: M.P.B., M.A.H., S.D.L., A.X.M., B.P.F.R., C.H.V., A.S.Z., and X.Z. PI of the EWAS group: M.W.L., C.M.N., A.K.S., and M.U.

## Acknowledgments

PGC-PTSD Epigenetics Workgroup: Reid S Alisch, Ananda B Amstadter, Don Armstrong, Archana Basu, Jean C Beckham, Nicole L Bjorklund, Barbara H Chaiyachati, Judith B M Ensink, Segun Fatumo, Leland L Fleming, Sandro Galea, Joel Gelernter, Ryan J Herringa, Sonia Jain, Diana L Juvinao-Quintero, Seyma Katrinli, Elizabeth Ketema, José J Martínez-Magaña, Burook Misganaw, Shiela Tiemi Nagamatsu, Danny M Nispeling, John Pfeiffer, Christian Schmahl, Gen Shinozaki, Clara Snijders, Jennifer A Sumner, Patricia C Swart, Audrey Tyrka, Robert J Ursano, Mirjam van Zuiden, Eric Vermetten, Jaqueline S Womersley, Nagy A Youssef, Yuanchao Zheng, Yiwen Zhu, Lea Zillich PsychENCODE PTSD Brainomics Project: Dhivya Arasappan, Sabina Berretta, Rahul A. Bharadwaj, Frances A. Champagne, Leonardo Collado-Torres, Christos Chatzinakos, Nikolaos P. Daskalakis, Chris P. DiPietro, Duc M. Duong, Amy Deep-Soboslay, Nick Eagles, Louise Huuki, Thomas Hyde, Artemis Iatrou, Aarti Jajoo, Joel E. Kleinman, Charles B. Nemeroff, Geo Pertea, Deanna Ross, Nicholas T. Seyfried, Joo Heon Shin, Kerry J.Ressler, Clara Snijders, Ran Tao, Daniel R. Weinberger, Stefan Wuchty, Dennis Wylie Traumatic Stress Brain Research Group: Victor E. Alvarez, David Benedek, Alicia Che, Dianne A. Cruz, David A. Davis, Matthew J. Girgenti, Ellen Hoffman, Paul E. Holtzheimer, Bertrand R. Huber, Alfred Kaye, John H. Krystal, Adam T. Labadorf, Terence M. Keane, Ann McKee, Brian Marx, Crystal Noller, Meghan Pierce, William K. Scott, Paula Schnurr, Krista DiSano, Thor Stein, Douglas E. Williamson, Keith A. Young

*   Received July 15, 2024.
*   Revision received July 15, 2024.
*   Accepted July 15, 2024.


*   © 2024, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## References

1.  1.American Psychiatric Association. Diagnostic and Statistical Manual of Mental Disorders. 2013.
    
    

2.  2.Zen AL, Whooley MA, Zhao S, Cohen BE. Post-traumatic stress disorder is associated with poor health behaviors: findings from the heart and soul study. Health Psychol. 2012;31(2):194–201.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1037/a0025989&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22023435&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F15%2F2024.07.15.24310422.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000301392800009&link_type=ISI) 

3.  3.Jacobsen LK, Southwick SM, Kosten TR. Substance use disorders in patients with posttraumatic stress disorder: a review of the literature. Am J Psychiatry. 2001;158(8):1184–90.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1176/appi.ajp.158.8.1184&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11481147&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F15%2F2024.07.15.24310422.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000170246800002&link_type=ISI) 

4.  4.O’Donovan A, Cohen BE, Seal KH, Bertenthal D, Margaretten M, Nishimi K, et al. Elevated risk for autoimmune disorders in iraq and afghanistan veterans with posttraumatic stress disorder. Biol Psychiatry. 2015;77(4):365–74.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2014.06.015&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25104173&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F15%2F2024.07.15.24310422.atom) 

5.  5.Ryder AL, Azcarate PM, Cohen BE. PTSD and Physical Health. Curr Psychiatry Rep. 2018;20(12):116.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s11920-018-0977-9&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30367276&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F15%2F2024.07.15.24310422.atom) 

6.  6.Benjet C, Bromet E, Karam EG, Kessler RC, McLaughlin KA, Ruscio AM, et al. The epidemiology of traumatic event exposure worldwide: results from the World Mental Health Survey Consortium. Psychol Med. 2016;46(2):327–43.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0033291715001981&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F15%2F2024.07.15.24310422.atom) 

7.  7.Logue MW, Amstadter AB, Baker DG, Duncan L, Koenen KC, Liberzon I, et al. The Psychiatric Genomics Consortium Posttraumatic Stress Disorder Workgroup: Posttraumatic Stress Disorder Enters the Age of Large-Scale Genomic Collaboration. Neuropsychopharmacology. 2015;40(10):2287–97.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/npp.2015.118&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25904361&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F15%2F2024.07.15.24310422.atom) 

8.  8.Nievergelt CM, Ashley-Koch AE, Dalvie S, Hauser MA, Morey RA, Smith AK, et al. Genomic Approaches to Posttraumatic Stress Disorder: The Psychiatric Genomic Consortium Initiative. Biol Psychiatry. 2018;83(10):831–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2018.01.020&link_type=DOI) 

9.  9.Polimanti R, Wendt FR. Posttraumatic stress disorder: from gene discovery to disease biology. Psychol Med. 2021;51(13):2178–88.
    
    

10. 10.Duncan LE, Cooper BN, Shen H. Robust Findings From 25 Years of PTSD Genetics Research. Curr Psychiatry Rep. 2018;20(12):115.
    
    

11. 11.Nievergelt CM, Maihofer AX, Atkinson EG, Chen C-Y, Choi KW, Coleman JR, et al. Discovery of 95 PTSD loci provides insight into genetic architecture and neurobiology of trauma and stress-related disorders. medRxiv. 2023:2023.08.31.23294915.
    
    

12. 12.Daskalakis NP, Rijal CM, King C, Huckins LM, Ressler KJ. Recent Genetics and Epigenetics Approaches to PTSD. Curr Psychiatry Rep. 2018;20(5):30.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F15%2F2024.07.15.24310422.atom) 

13. 13.Zannas AS, Provencal N, Binder EB. Epigenetics of Posttraumatic Stress Disorder: Current Evidence, Challenges, and Future Directions. Biol Psychiatry. 2015;78(5):327–35.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2015.04.003&link_type=DOI) 

14. 14.Bonasio R, Tu S, Reinberg D. Molecular signals of epigenetic states. Science. 2010;330(6004):612-6.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzMzAvNjAwNC82MTIiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wNy8xNS8yMDI0LjA3LjE1LjI0MzEwNDIyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

15. 15.Baker-Andresen D, Ratnu VS, Bredy TW. Dynamic DNA methylation: a prime candidate for genomic metaplasticity and behavioral adaptation. Trends Neurosci. 2013;36(1):3–13.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tins.2012.09.003&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23041052&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F15%2F2024.07.15.24310422.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000314258500002&link_type=ISI) 

16. 16.Rasmi Y, Shokati A, Hassan A, Aziz SG, Bastani S, Jalali L, et al. The role of DNA methylation in progression of neurological disorders and neurodegenerative diseases as well as the prospect of using DNA methylation inhibitors as therapeutic agents for such disorders. IBRO Neurosci Rep. 2023;14:28–37.
    
    

17. 17.Murgatroyd C, Wu Y, Bockmuhl Y, Spengler D. Genes learn from stress: how infantile trauma programs us for depression. Epigenetics. 2010;5(3):194–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4161/epi.5.3.11375&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20339319&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F15%2F2024.07.15.24310422.atom) 

18. 18.Lokk K, Modhukur V, Rajashekar B, Martens K, Magi R, Kolde R, et al. DNA methylome profiling of human tissues identifies global and tissue-specific methylation patterns. Genome Biol. 2014;15(4):r54.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/gb-2014-15-4-r54&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24690455&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F15%2F2024.07.15.24310422.atom) 

19. 19.Hannon E, Lunnon K, Schalkwyk L, Mill J. Interindividual methylomic variation across blood, cortex, and cerebellum: implications for epigenetic studies of neurological and neuropsychiatric phenotypes. Epigenetics. 2015;10(11):1024–32.
    
    

20. 20.Daskalakis NP, Cohen H, Nievergelt CM, Baker DG, Buxbaum JD, Russo SJ, et al. New translational perspectives for blood-based biomarkers of PTSD: From glucocorticoid to immune mediators of stress susceptibility. Exp Neurol. 2016;284(Pt B):133–40.
    
    

21. 21.Morrison FG, Miller MW, Logue MW, Assef M, Wolf EJ. DNA methylation correlates of PTSD: Recent findings and technical challenges. Prog Neuropsychopharmacol Biol Psychiatry. 2019;90:223–34.
    
    

22. 22.Katrinli S, Oliveira NCS, Felger JC, Michopoulos V, Smith AK. The role of the immune system in posttraumatic stress disorder. Transl Psychiatry. 2022;12(1):313.
    
    

23. 23.Uddin M, Aiello AE, Wildman DE, Koenen KC, Pawelec G, de Los Santos R, et al. Epigenetic and immune function profiles associated with posttraumatic stress disorder. Proc Natl Acad Sci U S A. 2010;107(20):9470–5.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTA3LzIwLzk0NzAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wNy8xNS8yMDI0LjA3LjE1LjI0MzEwNDIyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

24. 24.Smith AK, Conneely KN, Kilaru V, Mercer KB, Weiss TE, Bradley B, et al. Differential immune system DNA methylation and cytokine regulation in post-traumatic stress disorder. Am J Med Genet B Neuropsychiatr Genet. 2011;156B(6):700–8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ajmg.b.31212&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21714072&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F15%2F2024.07.15.24310422.atom) 

25. 25.Kuan PF, Waszczuk MA, Kotov R, Marsit CJ, Guffanti G, Gonzalez A, et al. An epigenome-wide DNA methylation study of PTSD and depression in World Trade Center responders. Transl Psychiatry. 2017;7(6):e1158.
    
    

26. 26.Rutten BPF, Vermetten E, Vinkers CH, Ursini G, Daskalakis NP, Pishva E, et al. Longitudinal analyses of the DNA methylome in deployed military servicemen identify susceptibility loci for post-traumatic stress disorder. Mol Psychiatry. 2018;23(5):1145–56.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/mp.2017.120&link_type=DOI) 

27. 27.Logue MW, Miller MW, Wolf EJ, Huber BR, Morrison FG, Zhou Z, et al. An epigenome-wide association study of posttraumatic stress disorder in US veterans implicates several new DNA methylation loci. Clin Epigenetics. 2020;12(1):46.
    
    

28. 28.Hammamieh R, Chakraborty N, Gautam A, Muhie S, Yang R, Donohue D, et al. Whole-genome DNA methylation status associated with clinical PTSD measures of OIF/OEF veterans. Transl Psychiatry. 2017;7(7):e1169.
    
    

29. 29.Ratanatharathorn A, Boks MP, Maihofer AX, Aiello AE, Amstadter AB, Ashley-Koch AE, et al. Epigenome-wide association of PTSD from heterogeneous cohorts with a common multi-site analysis pipeline. Am J Med Genet B Neuropsychiatr Genet. 2017;174(6):619–30.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ajmg.b.32568&link_type=DOI) 

30. 30.Uddin M, Ratanatharathorn A, Armstrong D, Kuan PF, Aiello AE, Bromet EJ, et al. Epigenetic meta-analysis across three civilian cohorts identifies NRG1 and HGS as blood-based biomarkers for post-traumatic stress disorder. Epigenomics. 2018;10(12):1585–601.
    
    

31. 31.Smith AK, Ratanatharathorn A, Maihofer AX, Naviaux RK, Aiello AE, Amstadter AB, et al. Epigenome-wide meta-analysis of PTSD across 10 military and civilian cohorts identifies methylation changes in AHRR. Nat Commun. 2020;11(1):5965.
    
    

32. 32.Snijders C, Maihofer AX, Ratanatharathorn A, Baker DG, Boks MP, Geuze E, et al. Longitudinal epigenome-wide association studies of three male military cohorts reveal multiple CpG sites associated with post-traumatic stress disorder. Clin Epigenetics. 2020;12(1):11.
    
    

33. 33.Katrinli S, Maihofer AX, Wani AH, Pfeiffer JR, Ketema E, Ratanatharathorn A, et al. Epigenome-wide meta-analysis of PTSD symptom severity in three military cohorts implicates DNA methylation changes in genes involved in immune system and oxidative stress. Mol Psychiatry. 2022;27(3):1720–8.
    
    

34. 34.McCartney DL, Walker RM, Morris SW, McIntosh AM, Porteous DJ, Evans KL. Identification of polymorphic and off-target probe binding sites on the Illumina Infinium MethylationEPIC BeadChip. Genom Data. 2016;9:22–4.
    
    

35. 35.Fortin JP, Triche TJ, Jr., Hansen KD. Preprocessing, normalization and integration of the Illumina HumanMethylationEPIC array with minfi. Bioinformatics. 2017;33(4):558–60.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btw691&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F15%2F2024.07.15.24310422.atom) 

36. 36.Leek JT, Johnson WE, Parker HS, Jaffe AE, Storey JD. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics. 2012;28(6):882–3.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/bts034&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22257669&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F15%2F2024.07.15.24310422.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000301972900020&link_type=ISI) 

37. 37.Teschendorff AE, Breeze CE, Zheng SC, Beck S. A comparison of reference-based algorithms for correcting cell-type heterogeneity in Epigenome-Wide Association Studies. BMC Bioinformatics. 2017;18(1):105.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12859-017-1511-5&link_type=DOI) 

38. 38.Salas LA, Koestler DC, Butler RA, Hansen HM, Wiencke JK, Kelsey KT, et al. An optimized library for reference-based deconvolution of whole-blood biospecimens assayed using the Illumina HumanMethylationEPIC BeadArray. Genome Biol. 2018;19(1):64.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13059-018-1448-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F15%2F2024.07.15.24310422.atom) 

39. 39.Barfield RT, Almli LM, Kilaru V, Smith AK, Mercer KB, Duncan R, et al. Accounting for population stratification in DNA methylation studies. Genet Epidemiol. 2014;38(3):231–41.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/gepi.21789&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24478250&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F15%2F2024.07.15.24310422.atom) 

40. 40.van Iterson M, van Zwet EW, Consortium B, Heijmans BT. Controlling bias and inflation in epigenome- and transcriptome-wide association studies using the empirical null distribution. Genome Biol. 2017;18(1):19.
    
    

41. 41.Barfield RT, Kilaru V, Smith AK, Conneely KN. CpGassoc: an R function for analysis of DNA methylation microarray data. Bioinformatics. 2012;28(9):1280–1.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/bts124&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22451269&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F15%2F2024.07.15.24310422.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000303338200016&link_type=ISI) 

42. 42.Balduzzi S, Rucker G, Schwarzer G. How to perform a meta-analysis with R: a practical tutorial. Evid Based Ment Health. 2019;22(4):153–60.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F15%2F2024.07.15.24310422.atom) 

43. 43.Mansell G, Gorrie-Stone TJ, Bao Y, Kumari M, Schalkwyk LS, Mill J, et al. Guidance for DNA methylation studies: statistical insights from the Illumina EPIC array. BMC Genomics. 2019;20(1):366.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12864-019-5761-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31088362&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F15%2F2024.07.15.24310422.atom) 

44. 44.Phipson B, Maksimovic J, Oshlack A. missMethyl: an R package for analyzing data from Illumina’s HumanMethylation450 platform. Bioinformatics. 2016;32(2):286–8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btv560&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26424855&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F15%2F2024.07.15.24310422.atom) 

45. 45.Logue MW, Miller MW, Wolf EJ, Huber BR, Morrison FG, Zhou Z, et al. An epigenome-wide association study of posttraumatic stress disorder in US veterans implicates several new DNA methylation loci. Clinical Epigenetics. 2020; IN PRESS.
    
    

46. 46.Daskalakis NP, Iatrou A, Chatzinakos C, Jajoo A, Snijders C, DiPietro CP, et al. Systems Biology Dissection of PTSD and MDD across brain regions, cell-types and blood. Science. In Revision.
    
    

47. 47.Friedman MJ, Huber BR, Brady CB, Ursano RJ, Benedek DM, Kowall NW, et al. VA’s National PTSD Brain Bank: a National Resource for Research. Curr Psychiatry Rep. 2017;19(10):73.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s11920-017-0822-6&link_type=DOI) 

48. 48.Rompala G, Nagamatsu ST, Martinez-Magana JJ, Nunez-Rios DL, Wang J, Girgenti MJ, et al. Profiling neuronal methylome and hydroxymethylome of opioid use disorder in the human orbitofrontal cortex. Nat Commun. 2023;14(1):4544.
    
    

49. 49.Zannas AS. Naturalistic Stress Hormone Levels Drive Cumulative Epigenomic Changes along the Cellular Lifespan. Int J Mol Sci. 2021;22(16).
    
    

50. 50.Leung CS, Kosyk O, Welter EM, Dietrich N, Archer TK, Zannas AS. Chronic stress-driven glucocorticoid receptor activation programs key cell phenotypes and functional epigenomic patterns in human fibroblasts. iScience. 2022;25(9):104960.
    
    

51. 51.Min JL, Hemani G, Hannon E, Dekkers KF, Castillo-Fernandez J, Luijk R, et al. Genomic and phenotypic insights from an atlas of genetic effects on DNA methylation. Nat Genet. 2021;53(9):1311–21.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-021-00923-x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F15%2F2024.07.15.24310422.atom) 

52. 52.Villicana S, Castillo-Fernandez J, Hannon E, Christiansen C, Tsai PC, Maddock J, et al. Genetic impacts on DNA methylation help elucidate regulatory genomic processes. Genome Biol. 2023;24(1):176.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13059-023-03011-x&link_type=DOI) 

53. 53.Villicaña S, Castillo-Fernandez J, Hannon E, Christiansen C, Tsai P-C, Maddock J, et al. Genetic impacts on DNA methylation help elucidate regulatory genomic processes. bioRxiv. 2023:2023.03.31.535045.
    
    

54. 54.Chen XQ, Tan I, Leung T, Lim L. The myotonic dystrophy kinase-related Cdc42-binding kinase is involved in the regulation of neurite outgrowth in PC12 cells. J Biol Chem. 1999;274(28):19901–5.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjEyOiIyNzQvMjgvMTk5MDEiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wNy8xNS8yMDI0LjA3LjE1LjI0MzEwNDIyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

55. 55.Story Jovanova O, Nedeljkovic I, Spieler D, Walker RM, Liu C, Luciano M, et al. DNA Methylation Signatures of Depressive Symptoms in Middle-aged and Elderly Persons: Meta-analysis of Multiethnic Epigenome-wide Studies. JAMA Psychiatry. 2018;75(9):949–59.
    
    

56. 56.Ligthart S, Marzi C, Aslibekyan S, Mendelson MM, Conneely KN, Tanaka T, et al. DNA methylation signatures of chronic low-grade inflammation are associated with complex diseases. Genome Biol. 2016;17(1):255.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13059-016-1119-5&link_type=DOI) 

57. 57.Wielscher M, Mandaviya PR, Kuehnel B, Joehanes R, Mustafa R, Robinson O, et al. DNA methylation signature of chronic low-grade inflammation and its role in cardio-respiratory diseases. Nat Commun. 2022;13(1):2408.
    
    

58. 58.Muniz Carvalho C, Wendt FR, Maihofer AX, Stein DJ, Stein MB, Sumner JA, et al. Dissecting the genetic association of C-reactive protein with PTSD, traumatic events, and social support. Neuropsychopharmacology. 2021;46(6):1071–7.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41386-020-0655-6&link_type=DOI) 

59. 59.Eraly SA, Nievergelt CM, Maihofer AX, Barkauskas DA, Biswas N, Agorastos A, et al. Assessment of plasma C-reactive protein as a biomarker of posttraumatic stress disorder risk. JAMA Psychiatry. 2014;71(4):423–31.
    
    

60. 60.Yang JJ, Jiang W. Immune biomarkers alterations in post-traumatic stress disorder: A systematic review and meta-analysis. J Affect Disord. 2020;268:39–46.
    
    

61. 61.Passos IC, Vasconcelos-Moreno MP, Costa LG, Kunz M, Brietzke E, Quevedo J, et al. Inflammatory markers in post-traumatic stress disorder: a systematic review, meta-analysis, and meta-regression. Lancet Psychiatry. 2015;2(11):1002–12.
    
    

62. 62.Hao N, Whitelaw ML. The emerging roles of AhR in physiology and immunity. Biochem Pharmacol. 2013;86(5):561–70.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bcp.2013.07.004&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23856287&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F15%2F2024.07.15.24310422.atom) 

63. 63.Wilson R, Wahl S, Pfeiffer L, Ward-Caviness CK, Kunze S, Kretschmer A, et al. The dynamics of smoking-related disturbed methylation: a two time-point study of methylation change in smokers, non-smokers and former smokers. BMC Genomics. 2017;18(1):805.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12864-017-4198-0&link_type=DOI) 

64. 64.Katrinli S, Smith AK. Immune system regulation and role of the human leukocyte antigen in posttraumatic stress disorder. Neurobiol Stress. 2021;15:100366.
    
    

65. 65.Kim GS, Smith AK, Xue F, Michopoulos V, Lori A, Armstrong DL, et al. Methylomic profiles reveal sex-specific differences in leukocyte composition associated with post-traumatic stress disorder. Brain Behav Immun. 2019;81:280–91.