- Open Access
Longitudinal epigenome-wide association studies of three male military cohorts reveal multiple CpG sites associated with post-traumatic stress disorder
Clinical Epigenetics volume 12, Article number: 11 (2020)
Epigenetic mechanisms have been suggested to play a role in the development of post-traumatic stress disorder (PTSD). Here, blood-derived DNA methylation data (HumanMethylation450 BeadChip) collected prior to and following combat exposure in three cohorts of male military members were analyzed to assess whether DNA methylation profiles are associated with the development of PTSD. A total of 123 PTSD cases and 143 trauma-exposed controls were included in the analyses. The Psychiatric Genomics Consortium (PGC) PTSD EWAS QC pipeline was used on all cohorts, and results were combined using a sample size weighted meta-analysis in a two-stage design. In stage one, we jointly analyzed data of two new cohorts (N = 126 and 78) for gene discovery, and sought to replicate significant findings in a third, previously published cohort (N = 62) to assess the robustness of our results. In stage 2, we aimed at maximizing power for gene discovery by combining all three cohorts in a meta-analysis.
Stage 1 analyses identified four CpG sites in which, conditional on pre-deployment DNA methylation, post-deployment DNA methylation was significantly associated with PTSD status after epigenome-wide adjustment for multiple comparisons. The most significant (intergenic) CpG cg05656210 (p = 1.0 × 10−08) was located on 5q31 and significantly replicated in the third cohort. In addition, 19 differentially methylated regions (DMRs) were identified, but failed replication. Stage 2 analyses identified three epigenome-wide significant CpGs, the intergenic CpG cg05656210 and two additional CpGs located in MAD1L1 (cg12169700) and HEXDC (cg20756026). Interestingly, cg12169700 had an underlying single nucleotide polymorphism (SNP) which was located within the same LD block as a recently identified PTSD-associated SNP in MAD1L1. Stage 2 analyses further identified 12 significant differential methylated regions (DMRs), 1 of which was located in MAD1L1 and 4 were situated in the human leukocyte antigen (HLA) region.
This study suggests that the development of combat-related PTSD is associated with distinct methylation patterns in several genomic positions and regions. Our most prominent findings suggest the involvement of the immune system through the HLA region and HEXDC, and MAD1L1 which was previously associated with PTSD.
Post-traumatic stress disorder (PTSD) is a debilitating psychiatric disorder that can develop following direct or indirect exposure to a potentially life-threatening traumatic incident. Symptoms include persistent re-experiencing of the trauma, avoidance behavior, hyperarousal, and negative mood . Although most individuals have the potential to withstand negative effects of trauma exposure on long-term mental health and to recover promptly, some are more vulnerable and at increased risk of developing PTSD. Understanding the molecular and neurobiological underpinnings of this differential susceptibility is currently receiving considerable attention, and epigenetic mediation of environmental influences has been proposed as a potential key mechanism [2,3,4].
Several epigenome-wide association studies (EWAS) have aimed to identify differentially methylated CpGs in PTSD [5,6,7,8]. However, most of these studies are based on association analyses where methylation was assessed at a single time point (cross-sectional), with limited ability to adjust for confounding variables. Only one PTSD study to date reported longitudinal changes in methylation profiles across a period of combat exposure in order to capture changes in DNA methylation over time in relation to phenotypic changes . We made use of the Prospective Research In Stress-related Military Operations (PRISMO) study, which in the present study was used as a replication cohort.
Here, we followed a previously published two-stage design  where we first meta-analyzed two longitudinal, USA-based military cohorts in order to identify associations between changes in methylation levels from pre-deployment to post-deployment and the development of PTSD. We then sought replication of our significant findings in PRISMO. In the second stage, we combined all three cohorts in a meta-analysis. This two-stage approach allows us to investigate the robustness of our findings through replication in stage 1, while increasing power for gene discovery by combining all three studies in stage 2. For all three cohorts, DNA methylation data and phenotypic data were collected prior to and following a 4–7-month deployment to an active ware zone in Iraq or Afghanistan. All studies selected PTSD cases and controls at post-deployment and only included subjects without PTSD at pre-deployment. Of the significant CpGs in the second analysis stage, we assessed associations with nearby single nucleotide polymorphisms (SNPs) and gene expression data, and examined correlations between blood and brain methylation status. To the best of our knowledge, this is the largest study aimed at detecting methylation changes associated with the development of PTSD. This prospective and longitudinal analysis permits us to more accurately capture dynamic changes in DNA methylation in relation to PTSD development while minimizing confounding due to intra-individual variability.
Three military cohorts were included in this study, i.e., the US Marine Resiliency Study (MRS), the US Army Study to Assess Risk and Resilience in Servicemembers (Army STARRS), and the Dutch PRISMO study. Demographic and clinical characteristics of subjects from all three cohorts (total N subjects = 266) can be found in Table 1. All subjects were male, and the majority were of European ancestry (N = 211, 79%). Within each cohort, cases and controls did not differ significantly in terms of age. Pre-deployment PTSD symptoms were significantly different between cases and controls from MRS only, with cases scoring slightly higher on the Clinician Administered PTSD Scale (CAPS) as compared to controls (p = .002; Table 1). In MRS and Army STARRS, cases were exposed to more traumatic events as compared to controls (p < .001 for both cohorts).
Stage 1: meta-analysis of MRS and Army STARRS
Data from MRS and Army STARRS were combined in order to identify CpG sites in which, conditional on baseline DNA methylation, post-deployment methylation was associated with PTSD status. Four genome-wide significant CpG sites (i.e., differentially methylated positions, DMPs) were identified using a conservative Bonferroni threshold of p = 1.13 × 10−07 for the 450K EWAS array (Table 2). These sites were located near SPRY4, in SDK1, CTRC, and CDH15. The direction of DNA methylation profiles associated with PTSD development was different for each site (Additional file 1: Figures S1-4). Additionally, after the Bonferroni correction for ~ 26,000 predefined regions, 19 differentially methylated regions (DMRs) were identified in which, conditional on baseline DNA methylation, post-deployment methylation was significantly associated with PTSD status (Table 3).
Replication in PRISMO
After the Bonferroni correction for the four significant DMPs and when using a one-sided test, the association of one CpG site, the intergenic site cg05656210, was replicated in PRISMO (p = 2.0 × 10−02; Table 2). Both the discovery meta-analysis and replication analysis show decreased DNA methylation in association with PTSD status. None of the 19 significant DMRs were replicated in PRISMO (Table 3).
Stage 2: meta-analysis of MRS, Army STARRS, and PRISMO
When combining MRS, Army STARRS, and PRISMO, the DNA methylation profile of three CpG sites was significantly associated with post-deployment PTSD status (Table 2, Fig. 1). The intergenic CpG that replicated in PRISMO remained the most significant (Z = − 6.14, p = 8.1 × 10−10). The other sites were located in the gene body regions of MAD1L1 and HEXDC (Additional file 1: Figures S1, S5, S6). Sensitivity analyses for the potentially confounding effects of changes in smoking and alcohol use did not substantially affect these results (Additional file 1: Table S1). Furthermore, 12 DMRs were identified (Additional file 1: Figures S7-18, Fig. 1), 7 of which were also significant in stage 1 and 4 were located in the human leukocyte antigen (HLA) region (Table 3).
Genetic effects and gene expression
Using MRS data, genetic effects on DNA methylation levels of the significant DMPs were assessed by testing for associations with SNPs within 500 kb of the DMPs. All DMPs had significantly associated SNPs which explained approximately 80% of the variation in methylation (p < 2 × 10−16) and were located within 1 bp of their respective CpG sites (Additional file 1: Table S2). However, adjusting for genotypes in the main model to assess the impact of SNPs on the association between DNA methylation and PTSD did not substantially affect the observed findings (Additional file 1: Table S3). We further assessed the association between baseline or post-deployment methylation signatures of these DMPs and corresponding blood-derived gene expression data which was available for MRS . At baseline, methylation levels of the CpGs located in HEXDC and MAD1L1 were significantly correlated with gene expression data (Table 4), with an inverse correlation between DNA methylation and expression in HEXDC and a positive correlation between methylation and expression of MAD1L1. Post-deployment DNA methylation values correlated to gene expression of HEXDC only. Methylation of the intergenic site cg05656210 was not significantly associated with expression of the closest gene at any time point. We additionally assessed the association between changes in expression and changes in methylation for these DMPs and found no significant associations (Table 4).
Blood-brain correlations of PTSD-associated CpGs
Blood-brain correlations of methylation levels of the significant stage 2 DMPs were examined using a publicly available database . For all three DMPs, blood DNA methylation levels correlated strongly with those in the prefrontal cortex, entorhinal cortex, superior temporal gyrus, and cerebellum (r ≥ 0.93 for all sites; p values ranging between 1.48 × 10−32 and 5.32 × 10−72; Additional file 1: Table S4, Additional file 1: Figure S19 for cg05656210).
Exposure to trauma is a prerequisite for the development of PTSD, yet not all individuals develop PTSD following trauma . The underlying biological mechanisms of this differential susceptibility have not yet been fully identified, and even the largest genome-wide association studies (GWAS) to date explain only a small proportion of the disease liability [13, 14]. Epigenetic changes have been studied as one potential mechanism, but most association studies have used cross-sectional designs which render it impossible to establish causality. Here, we use a more powerful longitudinal design to investigate associations of DNA methylation with post-deployment PTSD status across very similar military cohorts deployed to combat in Iraq and Afghanistan. We started with combining the USA-based MRS and Army STARRS cohorts and sought replication using the previously published Dutch PRISMO study . To maximize power for new discoveries, we also performed a meta-analysis across all three cohorts. The first analysis stage revealed 4 genome-wide significant DMPs and 19 DMRs which were linked to post-deployment PTSD status. One of these DMPs replicated in PRISMO. In the second stage, a meta-analysis of all 3 studies revealed that the replicating DMP and 7 DMRs remained significant, and 2 additional DMPs and 12 DMRs were identified.
The replicating DMP cg05656210 remained the top-ranked significant marker in the second analysis stage. This CpG site is an intergenic site annotated near SPRY4. SPRY4 is a member of the Sprouty proteins which are mainly involved in inhibiting receptor tyrosine kinase (RTK) signaling . Upon activation by growth factor ligands, RTK signaling has a wide variety of downstream effects ranging from the regulation of cell proliferation and differentiation to the modulation of cellular metabolism . One particular receptor involved in RTK signaling, i.e., receptor tyrosine kinase B (TrkB), and its main ligand, brain-derived neurotrophic factor (BDNF), have repeatedly been shown to be affected in stress-related disorders such as depression . Consistently, two independent studies reported decreased mRNA levels of BDNF and TrkB in the prefrontal cortex and hippocampus of individuals who committed suicide as compared to healthy control subjects [18, 19]. Moreover, SPRY4-IT1, a long non-coding RNA derived from the second intron of SPRY4 , has been shown to interact with SKA2 , a gene that was suggested to be a promising biomarker for suicidal behavior [22, 23], stress susceptibility, and stress-related disorders such as PTSD [23,24,25]. Finally, SPRY4 was previously found differentially methylated in the blood of patients diagnosed with schizophrenia . Together, these results suggest that alterations within SPRY4 could contribute to psychiatric disorders such as depression, PTSD, and schizophrenia and potentially play a role in suicidal behavior. The question as to whether and how the identified DMP influences the expression of SPRY4 and is involved in these phenotypes still remains to be answered.
The second top significant probe, cg12169700, and one DMR are located within MAD1L1. MAD1L1 is part of the mitotic spindle-assembly checkpoint (SAC) which monitors the proper attachment of chromosomes to the microtubule spindle apparatus, delays the start of anaphase until all chromosomes are properly attached, and in doing so, ensures correct chromosome separation [27, 28]. Malfunctions of the MAD1L1 protein could therefore contribute to aneuploidy and chromosomal instability. Specific SNPs within this gene have previously been associated with bipolar disorder [29, 30], schizophrenia [30,31,32], and depression . Interestingly, MAD1L1 was recently identified in a PTSD GWAS of the Million Veteran Program (MVP) . The SNP that underlies cg12169700, rs11761270, is located in the same large linkage disequilibrium (LD) block as the MVP finding. In the MVP, carriers of the minor allele of rs11761270 showed decreased levels of methylation and were at increased risk of having PTSD. This corresponds to our own findings in which PTSD cases show a reduction in methylation from pre- to post-deployment. Moreover, using expression data from MRS, we found that methylation at this site was positively associated with gene expression of MAD1L1. This also aligns with previous findings that showed that blood levels of MAD1L1 were decreased in highly stress-susceptible individuals . Together, these findings suggest that specific methylation profiles within MAD1L1 may be regarded as a risk factor for PTSD in addition to several other psychiatric disorders . However, the underlying mechanisms through which such disturbances could contribute to psychiatric disorders warrant further research.
The third CpG site is located in HEXDC which to date has no known implications in any psychiatric disease. Although the biological functions of its product, hexosaminidase D, remain largely unknown, it is believed to be a glycosidase and previous studies found associations with rheumatoid arthritis [37, 38]. Interestingly, recent studies found significant pleiotropy between rheumatoid arthritis and PTSD , and rheumatoid arthritis and schizophrenia . Other studies found associations between PTSD symptoms and rheumatoid arthritis in a twin population  and in an epidemiological study of military veterans . This link between mental disorders and immune-related processes is discussed more in detail below. The DMP of HEXDC was located directly adjacent to rs4789774, a known expression quantitative trait locus (eQTL) that regulates the expression of HEXDC in the human brain cortex and of NARF and NARF-IT1 in a number of tissue types including blood (http://genome.ucsc.edu/). Moreover, a modest negative correlation was found between methylation of this site and gene expression of HEXDC, both at baseline and at post-deployment.
Twelve significant DMRs were found in the second phase of the analysis. Our strongest finding was in the HLA region, a gene-dense region which contains over 200 genes that encode human leukocyte antigen complex proteins in charge of presenting peptide antigens to trigger immune reactions, among other (non-)immune functions . Their non-immunological roles include processes such as neurodevelopment, synaptic plasticity, learning, memory, and stress reactivity [44, 45]. It is therefore not entirely surprising that several epigenetic modifications and genetic variants within this region have repeatedly been found implicated in neuropsychiatric disorders such as schizophrenia and bipolar disorder (recently reviewed in ). Along with our HEXDC finding, these observations further enhance the existing notion that immune factors play an important role and should continue to be studied in relation to mental disorders such as PTSD . Although it is now clear that immune imbalances are present in PTSD, questions related to causality and further implications for prevention strategies and treatment options largely remain to be answered.
Follow-up analyses were done using the significant DMPs from the stage 2 meta-analysis only. The discovery that methylation levels at the top three PTSD-associated CpGs were highly associated with the genotype of the nearby SNPs led us to question whether the associations between methylation and PTSD status were mainly driven by genotype. However, direct adjustment for genotype in a sensitivity analysis did not attenuate the associations between DNA methylation and PTSD status and our current sample size limits our ability to conduct analyses specific to genotype strata to further investigate interaction effects between SNPs and methylation.
Since our methylation data were based on DNA from peripheral blood, we further examined correlations between blood and several brain regions, i.e., the prefrontal cortex, entorhinal cortex, superior temporal gyrus, and cerebellum. The results indicate that blood-brain correlations of all top DMPs were strong for all four brain regions suggesting that these findings could potentially also be relevant for tissues other than blood. Assessing these correlations is relevant when dealing with disorders such as PTSD which are characterized by functional and structural alterations within the brain but for which the accessibility to human brain tissue is limited. However, these and similar findings will need to be confirmed using postmortem brain tissue and their precise role in PTSD development will need to be investigated further.
The main limitation of the present study is its small sample size which likely captures only a fraction of all implicated CpGs and renders additional analyses such as pathway and network analyses underpowered. It further needs to be emphasized that this study used data generated with Illumina’s 450K arrays which only assess a subset of all CpG sites. Next, although examining blood-derived DNA methylation is informative when seeking relatively easily accessible biomarkers, follow-up studies are needed in order to assess these methylation patterns within the tissue of interest, i.e., the brain. Another important limitation of the study is the fact that subjects that developed PTSD were exposed to significantly more traumatic events than control subjects in MRS and Army STARRS. One way to address these differences would be to include trauma exposure as a covariate. However, given the high correlation between trauma exposure and PTSD, this would likely diminish the association of methylation changes and PTSD. Therefore, we acknowledge that our analyses may include associations with trauma exposure in addition to PTSD. Next, the limited replication of findings from the US cohorts in the Dutch PRISMO study, and vice versa , may point to type I errors, or be partially due to small sample sizes and/or heterogeneity in study designs, study environments, and potential confounders such as immune status or medication use. For example, findings reported from PRISMO  point towards the involvement of DMPs and DMRs located within ZFP57, RNF39, and HIST1H2APS2, which were not implicated in our analyses. However, these genes are located within the HLA region, a region which is indeed implicated in our analyses. Furthermore, whereas the PRISMO study is entirely based on subjects of European ancestry, the US MRS and Army STARRS studies include more ancestral diversity. Genes in the HLA region tend to be highly polymorphic, particularly with respect to ancestral background . Thus, it is possible that this heterogeneity may contribute to the lack of replication. Further analyses should be ancestry-specific, once sample sizes are adequate. Furthermore, at this stage, it is unclear whether the identified differential methylation patterns in PTSD cases have any functional consequences. Although they may influence gene expression, the current dataset has limited power to establish causality. In order to make claims regarding causation, performing functional studies in vitro and animal studies will be needed in order to unravel precise biological mechanisms. Finally, to maximize power for discovery, the present cohorts were chosen to be highly similar in regard to demographics, type of trauma, and time since trauma exposure. Thus, the degree to which these findings on active duty, predominantly European-ancestry military men, may generalize to females, civilians, or other ancestries is unclear.
In summary, this is the largest study on methylation changes associated with the development of combat-related PTSD to date. Our observations point towards the implication of biologically interesting genes such as HLA region which is involved in immune-related processes, HEXDC which also has been suggested to play a role in immunity, and MAD1L1, a PTSD-related gene recently identified in the large MVP. These findings strengthen the notion that specific DNA methylation profiles are involved in the development of combat-related PTSD. Larger longitudinal studies and integrative efforts are now needed to build upon these preliminary findings in order to understand their functional consequences and integrate them more broadly into our current understanding of the (epi)genomic basis of PTSD.
Marine Resiliency Study
The Marine Resiliency Study (MRS)  is a prospective PTSD study of Marines and Navy personnel deployed to Iraq or Afghanistan. PTSD symptoms were assessed approximately 1 month before deployment, and 3 and/or 6 months post-deployment using the CAPS and the PTSD Checklist (PCL) for DSM-IV. Biological samples including whole blood were collected at all time points. Information on smoking and alcohol use was collected on a self-report basis. Combat exposure was assessed approximately 1 week post-deployment using the Deployment Risk and Resilience Inventory (DRRI). A subset of 63 PTSD cases and 63 controls was selected for the methylation assays and inclusion in the present study. All subjects were free of a PTSD diagnosis at pre-deployment and had CAPS scores ≤ 25. After return from a ~ 7-month deployment period, blood samples from PTSD cases (following the DSM-IV full or partially stringent diagnosis [50, 51]) were selected either at the 3- or the 6-month follow-up visit, based on when these subjects had their highest recorded CAPS scores. Subsequently, controls were frequency matched to the selected cases for age, ancestry, and time of post-deployment visit. The study was approved by the institutional review boards of the University of California San Diego, VA San Diego Research Service, and Naval Health Research Center. All subjects provided informed consent.
The Army Study to Assess Risk and Resilience in Servicemembers (Army STARRS) is a prospective study among US Army personnel gathering information on risk and resilience factors for suicidality and psychopathology . All subjects completed a computerized version of the Composite International Diagnostic Interview screening scales (CIDI-SC) and the PCL6 screener for DSM-IV approximately 6 weeks before deployment to Afghanistan, and the PCL-C at 1, 2, and 6 months post-deployment. PTSD diagnosis was assigned using multiple imputation methods that relied on PCL and CIDI-SC data , and information on trauma exposure was gathered from self-administered questions on childhood-, adult-, and military-related events. Information on smoking and alcohol use was collected on a self-report basis. Biological samples including whole blood were collected approximately 6 weeks before deployment and 1 month post-deployment. A subset of 31 cases and 47 controls were selected for the methylation assays and inclusion in this analysis. All subjects were free of a PTSD diagnosis at pre-deployment. PTSD cases were selected based on their PTSD diagnosis at 6 months post-deployment. Controls were PTSD-free subjects matched on age, deployment stress, and childhood adversity. The study procedures were approved by the Institutional Review Boards of all collaborating organizations. All subjects provided informed consent.
Replication dataset: PRISMO
Replication data was obtained from the Prospective Research In Stress-related Military Operations (PRISMO) study, a prospective study of Dutch military soldiers deployed to Afghanistan [54, 55]. The severity of current PTSD symptoms was assessed using the Self-Report Inventory for PTSD (SRIP), and blood samples were collected approximately 1 month before and 1 and 6 months after deployment. Traumatic stress exposure during deployment to Afghanistan was assessed with a deployment experiences checklist. Information on smoking and alcohol use was collected on a self-report basis. A subset of 29 cases and 33 controls was selected for the methylation assays and inclusion in this analysis (see  for selection criteria). The study was approved by the ethical committee of the University Medical Center Utrecht, and was conducted in accordance with the Declaration of Helsinki. All subjects provided informed consent.
In all cohorts, longitudinal whole blood DNA methylation levels were measured using the Illumina HumanMethylation450K BeadChip. The Psychiatric Genomics Consortium (PGC)-EWAS quality control pipeline was used on all three cohorts . Briefly, samples were excluded when having a probe detection call rate < 90% and an average intensity value < 50% of the overall sample mean or < 2,000 arbitrary units (AU). Individual probes with detection p values > 0.001 or those based on less than three beads were set to missing. Remaining probes were excluded when cross-reactivity occurred between autosomal and sex chromosomes. CpG sites with missing data for > 10% of samples within cohorts were excluded (Additional file 1: Table S5). After filtering, the β-values reflecting methylation levels of individual cytosine residues were normalized to correct for differences between type I and type II probes using Beta Mixture Quantile Normalization (BMIQ) . ComBat  was used to correct for remaining issues such as batch and plate effects. To account for differences in cell type composition between samples, proportions of CD8, CD4, NK, B cells, monocytes, and granulocytes were estimated for each individual using their unique DNA methylation profiles. This was estimated using the estimatecellcounts function in minfi .
The normalized β-values were logit transformed to M values which were used for linear regression analysis. Post-deployment DNA methylation was modeled as a function of post-deployment PTSD status while adjusting for pre-deployment DNA methylation; age; changes in CD4T, CD8T, NK, B cell, and monocyte cell proportions; and principal components (PCs) for ancestry. For MRS and Army STARRS, the PCs were derived from available GWAS and PCs 1–3 were included. For PRISMO, the method described by Barfield and colleagues  was used to derive PCs from the EWAS data and PCs 2–4 (see ) were included. HC3 standard errors were calculated using the sandwich R library . Analyses were performed on each cohort independently, and the obtained p values were combined using a sample size weighted meta-analysis. Significance was declared at p < 1.13 × 10−7 after a stringent Bonferroni correction for 439,897 probes. Possible confounding effects of changes in smoking and alcohol use were assessed as a sensitivity analysis.
DMR analysis was performed on a set of 26,000 pre-defined gene regions within gene bodies, promoter regions, and CpG islands using the mCSEA version 1.2 package for R . Regions were included when annotated to having at least five CpGs. For each study, EWAS p values, methylation level values, and a phenotype and covariate data matrix were supplied as program inputs. p values were derived using 100,000 permutations. A sample size weighted meta-analysis of DMRs was performed based on z-score transformations of permutation p values. Significances of DMRs (p < 1.92 × 10−6) were derived based on a Bonferroni correction for the 26,000 tests performed. All positions and regions were in reference to GRCh37/hg19.
We considered replication as significant when the effect directions matched between the discovery and replication samples and the p values held up to the Bonferroni correction for the number of replications attempted (i.e., 4 for the DMPs, 19 for the DMRs) using a one-sided test.
Detecting genetic effects and links with gene expression
Associations between baseline levels of methylation of each significant CpG from the second analysis stage and nearby SNPs (within 500 kilobases; kb) were assessed in the MRS dataset using PLINK  to detect the potential influence of genetic effects on DNA methylation. For a given CpG site, the SNP with the lowest p value was carried forward as an additional covariate in the regression models as a sensitivity analysis.
For CpGs annotated to genes, we estimated the correlations between baseline or post-deployment CpG methylation levels and corresponding blood gene expression in MRS data. We additionally estimated the correlations between changes in expression and changes in methylation for these CpGs. Details of mRNA expression measurement in MRS can be found elsewhere . Briefly, from 124 MRS Marine participants, blood was drawn at ~ 1 month pre-deployment and mRNA levels were analyzed through sequencing on a Illumina Hi-Seq 2000.
We used the UCSC genome browser tool (http://genome.ucsc.edu/) to identify if SNPs associated with our CpGs influenced expression in other tissue types based on combined expression eQTL data from 44 tissues from GTEx v6 .
The Blood Brain DNA Methylation Comparison Tool (http://epigenetics.iop.kcl.ac.uk/bloodbrain/) was used to assess correlations between the methylation status of the top hits of the combined meta-analysis in blood and brain . Specifically, this tool yields Pearson’s correlation coefficients (r) and associated p values for the association of the methylation status of individual CpG sites in blood and the prefrontal cortex, entorhinal cortex, superior temporal gyrus, and cerebellum.
Availability of data and materials
Since the present study makes use of military datasets which are typically not shared publicly due to confidentiality reasons, the datasets used and/or analyzed during the current study are available from the corresponding author upon reasonable request only.
- Army STARRS:
Army Study to Assess Risk and Resilience in Servicemembers
Beta Mixture Quantile Normalization
Clinician Administered PTSD Scale
Differentially methylated position
Differentially methylated region
Deployment Risk and Resilience Inventory
Expression quantitative trait locus
Epigenome-wide association study
Eenome-wide association studies
Human leukocyte antigen
Marine Resiliency Study
Psychiatric Genomics Consortium
Prospective Research In Stress-related Military Operations
Post-traumatic stress disorder
Single nucleotide polymorphism
Self-Report Inventory for PTSD
Association AP. Diagnostic and statistical manual of mental disorders. 5th ed. Arlington: American Psychiatric Publishing; 2013.
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.
Daskalakis NP, Rijal CM, King C, Huckins LM, Ressler KJ. Recent genetics and epigenetics approaches to PTSD. Curr Psychiatry Rep. 2018;20(5):30.
Sheerin CM, Lind MJ, Bountress K, Nugent NR, Amstadter AB. The genetics and epigenetics of ptsd: overview, recent advances, and future directions. Curr Opin Psychol. 2017;14:5–11.
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.
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.
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.
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.
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.
Nievergelt CM, Maihofer AX, Mustapic M, Yurgil KA, Schork NJ, Miller MW, et al. Genomic predictors of combat stress vulnerability and resilience in U.S. Marines: a genome-wide association study across multiple ancestries implicates PRTFDC1 as a potential PTSD gene. (1873-3360 (Electronic)). 2015;51:459-71.
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.
Yehuda R, Hoge CW, McFarlane AC, Vermetten E, Lanius RA, Nievergelt CM, et al. Post-traumatic stress disorder. Nat Rev Dis Primers. 2015;1:15057.
Gelernter J, Sun N, Polimanti R, Pietrzak R, Levey DF, Bryois J, et al. Genome-wide association study of post-traumatic stress disorder reexperiencing symptoms in >165,000 US veterans. Nature Neuroscience. 2019;22(9):1394–401.
Nievergelt CM, Maihofer AX, Klengel T, Atkinson EG, Chen C-Y, Choi KW, et al. International meta-analysis of PTSD genome-wide association studies identifies sex- and ancestry-specific genetic risk loci. Nature Communications. 2019;10(1):4558.
Mason JM, Morrison DJ, Basson MA, Licht JD. Sprouty proteins: multifaceted negative-feedback regulators of receptor tyrosine kinase signaling. Trends Cell Biol. 2006;16(1):45–54.
Lemmon MA, Schlessinger J. Cell signaling by receptor tyrosine kinases. Cell. 2010;141(7):1117–34.
Zhang JC, Yao W, Hashimoto K. Brain-derived neurotrophic factor (BDNF)-TrkB signaling in inflammation-related depression and potential therapeutic targets. Curr Neuropharmacol. 2016;14(7):721–31.
Dwivedi Y, Rizavi HS, Conley RR, Roberts RC, Tamminga CA, Pandey GN. Altered gene expression of brain-derived neurotrophic factor and receptor tyrosine kinase B in postmortem brain of suicide subjects. Arch Gen Psychiatry. 2003;60(8):804–15.
Pandey GN, Ren X, Rizavi HS, Conley RR, Roberts RC, Dwivedi Y. Brain-derived neurotrophic factor and tyrosine kinase B receptor signalling in post-mortem brain of teenage suicide victims. Int J Neuropsychopharmacol. 2008;11(8):1047–61.
Khaitan D, Dinger ME, Mazar J, Crawford J, Smith MA, Mattick JS, et al. The melanoma-upregulated long noncoding RNA SPRY4-IT1 modulates apoptosis and invasion. Cancer Res. 2011;71(11):3852–62.
He XJ, Bian EB, Ma CC, Wang C, Wang HL, Zhao B. Long non-coding RNA SPRY4-IT1 promotes the proliferation and invasion of U251 cells through upregulation of SKA2. Oncol Lett. 2018;15(3):3977–84.
Guintivano J, Brown T, Newcomer A, Jones M, Cox O, Maher BS, et al. Identification and replication of a combined epigenetic and genetic biomarker predicting suicide and suicidal behaviors. Am J Psychiatry. 2014;171(12):1287–96.
Kaminsky Z, Wilcox HC, Eaton WW, Van Eck K, Kilaru V, Jovanovic T, et al. Epigenetic and genetic variation at SKA2 predict suicidal behavior and post-traumatic stress disorder. Transl Psychiatry. 2015;5:e627.
Boks MP, Rutten BP, Geuze E, Houtepen LC, Vermetten E, Kaminsky Z, et al. SKA2 methylation is involved in cortisol stress reactivity and predicts the development of post-traumatic stress disorder (PTSD) After Military Deployment. Neuropsychopharmacology. 2016;41(5):1350–6.
Sadeh N, Spielberg JM, Logue MW, Wolf EJ, Smith AK, Lusk J, et al. SKA2 methylation is associated with decreased prefrontal cortical thickness and greater PTSD severity among trauma-exposed veterans. Mol Psychiatry. 2016;21(3):357–63.
Zaharieva I, Georgieva L, Nikolov I, Kirov G, Owen MJ, O'Donovan MC, et al. Association study in the 5q31-32 linkage region for schizophrenia using pooled DNA genotyping. BMC Psychiatry. 2008;8:11.
Musacchio A, Salmon ED. The spindle-assembly checkpoint in space and time. Nat Rev Mol Cell Biol. 2007;8(5):379–93.
Bolanos-Garcia VM. Formation of multiprotein assemblies in the nucleus: the spindle assembly checkpoint. Int Rev Cell Mol Biol. 2014;307:151–74.
Trost S, Diekhof EK, Mohr H, Vieker H, Kramer B, Wolf C, et al. Investigating the impact of a genome-wide supported bipolar risk variant of MAD1L1 on the human reward system. Neuropsychopharmacology. 2016;41(11):2679–87.
Zhao L, Chang H, Zhou DS, Cai J, Fan W, Tang W, et al. Replicated associations of FADS1, MAD1L1, and a rare variant at 10q26.13 with bipolar disorder in Chinese population. Transl Psychiatry. 2018;8(1):270.
Ripke S, O'Dushlaine C, Chambert K, Moran JL, Kahler AK, Akterin S, et al. Genome-wide association analysis identifies 13 new risk loci for schizophrenia. Nat Genet. 2013;45(10):1150–9.
Su L, Shen T, Huang G, Long J, Fan J, Ling W, et al. Genetic association of GWAS-supported MAD1L1 gene polymorphism rs12666575 with schizophrenia susceptibility in a Chinese population. Neurosci Lett. 2016;610:98–103.
Wray NR, Ripke S, Mattheisen M, Trzaskowski M, Byrne EM, Abdellaoui A, et al. Genome-wide association analyses identify 44 risk variants and refine the genetic architecture of major depression. Nat Genet. 2018;50(5):668–81.
Gelernter J. ea. Genome-wide association study of posttraumatic stress disorder (PTSD) re-experiencing symptoms in >165,000 US veterans. Nature Neuroscience. 2019;In press.
Le-Niculescu H, Roseberry K, Levey DF, Rogers J, Kosary K, Prabha S, et al. Towards precision medicine for stress disorders: diagnostic biomarkers and targeted drugs. Molecular Psychiatry. 2019.
Levey DF, Gelernter J, Polimanti R, Zhou H, Cheng Z, Aslan M, et al. Reproducible risk loci and psychiatric comorbidities in anxiety: results from ~200,000 million veteran program participants. bioRxiv. 2019;540245.
Pasztoi M, Sodar B, Misjak P, Paloczi K, Kittel A, Toth K, et al. The recently identified hexosaminidase D enzyme substantially contributes to the elevated hexosaminidase activity in rheumatoid arthritis. Immunol Lett. 2013;149(1-2):71–6.
Alteen MG, Oehler V, Nemcovicova I, Wilson IB, Vocadlo DJ, Gloster TM. Mechanism of human nucleocytoplasmic hexosaminidase D. Biochemistry. 2016;55(19):2735–47.
Stein MB, Chen C-Y, Ursano RJ, Cai T, Gelernter J, Heeringa SG, et al. Genome-wide association studies of posttraumatic stress disorder in 2 cohorts of US army soldiers. JAMA Psychiatry. 2016;73(7):695–704.
Wang Q, Yang C, Gelernter J, Zhao H. Pervasive pleiotropy between psychiatric disorders and immune disorders revealed by integrative analysis of multiple GWAS. Hum Genet. 2015;134(11-12):1195–209.
Boscarino JA, Forsberg CW, Goldberg J. A twin study of the association between ptsd symptoms and rheumatoid arthritis. Psychosom Med. 2010;72(5):481–6.
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. Biological Psychiatry. 2015;77(4):365–74.
Beck S, Trowsdale J. The human major histocompatibility complex: lessons from the DNA sequence. Annu Rev Genomics Hum Genet. 2000;1(1):117–37.
Huh GS, Boulanger LM, Du H, Riquelme PA, Brotz TM, Shatz CJ. Functional requirement for class I MHC in CNS development and plasticity. Science (New York, NY). 2000;290(5499):2155–9.
Sankar A, MacKenzie RN, Foster JA. Loss of class I MHC function alters behavior and stress reactivity. J Neuroimmunol. 2012;244(1):8–15.
Debnath M, Berk M, Leboyer M, Tamouza R. The MHC/HLA gene complex in major psychiatric disorders: emerging roles and implications. Curr Behav Neurosci Rep. 2018;5(2):179–88.
Wang Z, Caughron B, Young MRI. Posttraumatic stress disorder: an immunological disorder? Front Psychiatry. 2017;8:222.
Katrinli S, Lori A, Kilaru V, Carter S, Powers A, Gillespie CF, et al. Association of HLA locus alleles with posttraumatic stress disorder. Brain Behav Immun. 2019;81:655–8.
Baker DG, Nash WP, Litz BT, Geyer MA, Risbrough VB, Nievergelt CM, et al. Predictors of risk and resilience for posttraumatic stress disorder among ground combat Marines: methods of the Marine Resiliency Study. Prev Chronic Dis. 2012;9:E97.
Sher L. Recognizing post-traumatic stress disorder. QJM. 2004;97(1):1–5.
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 novel methylation loci. bioRxiv. 2019;585109.
Ursano RJ, Colpe LJ, Heeringa SG, Kessler RC, Schoenbaum M, Stein MB, et al. The Army study to assess risk and resilience in servicemembers (Army STARRS). Psychiatry. 2014;77(2):107–19.
Kessler RC, Santiago PN, Colpe LJ, Dempsey CL, First MB, Heeringa SG, et al. Clinical reappraisal of the Composite International Diagnostic Interview Screening Scales (CIDI-SC) in the Army Study to Assess Risk and Resilience in Servicemembers (Army STARRS). Int J Methods Psychiatr Res. 2013;22(4):303–21.
Reijnen A, Rademaker AR, Vermetten E, Geuze E. Prevalence of mental health symptoms in Dutch military personnel returning from deployment to Afghanistan: a 2-year longitudinal analysis. Eur Psychiatry. 2015;30(2):341–6.
Eekhout I, Reijnen A, Vermetten E, Geuze E. Post-traumatic stress symptoms 5 years after military deployment to Afghanistan: an observational cohort study. Lancet Psychiatry. 2016;3(1):58–64.
Teschendorff AE, Marabita F, Lechner M, Bartlett T, Tegner J, Gomez-Cabrero D, et al. A beta-mixture quantile normalization method for correcting probe design bias in Illumina Infinium 450 k DNA methylation data. Bioinformatics. 2013;29(2):189–96.
Johnson WE, Li C, Rabinovic A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2007;8(1):118–27.
Houseman EA, Accomando WP, Koestler DC, Christensen BC, Marsit CJ, Nelson HH, et al. DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinformatics. 2012;13:86.
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.
Zeileis A. Econometric Computing with HC and HAC covariance matrix estimators. J Stat Softw. 2004;11:10.
Martorell-Marugan J, Gonzalez-Rumayor V, Carmona-Saez P. mCSEA: Detecting subtle differentially methylated regions; 2018. p. 293381.
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. American journal of human genetics. 2007;81(3):559–75.
Breen MS, Maihofer AX, Glatt SJ, Tylee DS, Chandler SD, Tsuang MT, et al. Gene networks specific for innate immunity define post-traumatic stress disorder. Molecular psychiatry. 2015;20(12):1538–45.
Consortium GT. The genotype-tissue expression (GTEx) project. Nature Genet. 2013;45(6):580–5.
Acknowledged are Mark A. Geyer (UCSD), Daniel T. O’Connor (UCSD), and all MRS investigators, as well as the MRS administrative core and data collection staff. The authors also thank the Marine and Navy Corpsmen volunteers for military service and participation in MRS.
This work was supported by the US Army Medical Research and Materiel Command as well as the National Institute of Mental Health (NIMH; R01MH108826; R01MH106595).
Funding for MRS was provided by the Marine Corps, Navy Bureau of Medicine and Surgery (BUMED), VA Health Research and Development (HSR&D), Veterans Affairs San Diego Healthcare System, Center of Excellence for Stress and Mental Health, and NIH R01MH093500. Data collection of PRISMO was funded by the Dutch Ministry of Defense, and DNA methylation analyses were funded by the VENI Award fellowship from the Netherlands Organization for Scientific Research (NWO, grant number 916.11.086). Army STARRS was sponsored by the Department of the Army and funded under cooperative agreement number U01MH087981 (2009–2015) with the National Institutes of Health, National Institute of Mental Health (NIH/NIMH).
Ethics approval and consent to participate
The MRS study was approved by the institutional review boards of the University of California San Diego, VA San Diego Research Service, and Naval Health Research Center. For Army STARRS, the study procedures were approved by the Institutional Review Boards of all collaborating organizations. The PRISMO study was approved by the ethical committee of the University Medical Center Utrecht, and was conducted in accordance with the Declaration of Helsinki. All subjects provided informed consent.
Consent for publication
RCK received support for his epidemiological studies from Sanofi Aventis; was a consultant for Johnson & Johnson Wellness and Prevention, Sage Pharmaceuticals, Shire, Takeda; and served on an advisory board for the Johnson & Johnson Services Inc. Lake Nona Life Project. Kessler is a co-owner of DataStat, Inc., a market research firm that carries out healthcare research. MBS has in the past 3 years been a consultant for Actelion, Aptinyx, Bionomics, Dart Neuroscience, Healthcare Management Technologies, Janssen, Neurocrine Biosciences, Oxeia Biopharmaceuticals, Pfizer, and Resilience Therapeutics. The other authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1. Differentially methylated positions (DMPs) in MRS, Army STARRS and PRISMO with and without corrections for smoking status and alcohol use. Table S2. SNPs within 500 bps upstream or downstream of the significant DMPs. Table S3. Differentially methylated positions (DMPs) in MRS with and without correction for main associated SNPs. Table S4. Correlations between blood and brain methylation levels for the top CpG sites based on external data. Table S5. Number of probes and samples removed at each stage of the quality control pipeline. Figure S1. Methylation values (B values) at cg05656210 for each cohort separately. Figure S2. Methylation values (B values) at cg05901543 for each cohort separately. Figure S3. Methylation values (B values) at cg16956686 for each cohort separately. Figure S4. Methylation values (B values) at cg18917957 for each cohort separately. Figure S5. Methylation values (B values) at cg12169700 for each cohort separately. Figure S6. Methylation values (B values) at cg20756026 for each cohort separately. Figure S7. Differentially methylated region at HLA-DBP1: 33048416-33048814. Figure S8. Differentially methylated region at HLA-DPB1: 33043976-33054001. Figure S9. Differentially methylated region at HLA-DRB1: 32547019-32557404. Figure S10. Differentially methylated region at HLA-DRB1: 32551851-32552331. Figure S11. Differentially methylated region at HOXA4: 27169572-27170638. Figure S12. Differentially methylated region at HOXA4: 27169740-27171528. Figure S13. Differentially methylated region at KCNE1: 35827824-35884508. Figure S14. Differentially methylated region at KCNE1: 35831697-35832365. Figure S15. Differentially methylated region at MAD1L1: 1885033-1885402. Figure S16. Differentially methylated region at NTRK1: 156814881-156815792. Figure S17. Differentially methylated region at SLC17A3: 25882327-25882560. Figure S18. Differentially methylated region at TRMT12: 125461772-125464547. Figure S19. Example of blood-brain correlations of methylation levels in cg05656210. PFC: prefrontal cortex, EC: entorhinal cortex, STG: superior temporal gyrus, CER: cerebellum.
About this article
Cite this article
Snijders, C., Maihofer, A.X., Ratanatharathorn, A. et al. Longitudinal epigenome-wide association studies of three male military cohorts reveal multiple CpG sites associated with post-traumatic stress disorder. Clin Epigenet 12, 11 (2020) doi:10.1186/s13148-019-0798-7
- DNA methylation