Posttraumatic stress disorder (PTSD) is a debilitating disorder that has been associated with brain abnormalities, including white matter alterations. However, little is known about the effect of treatment on these brain alterations. To investigate the course of white matter alterations in PTSD, we used a longitudinal design investigating treatment effects on white matter integrity using diffusion tensor imaging (DTI). Diffusion tensor and magnetization transfer images were obtained pre- and posttreatment from veterans with (n=39) and without PTSD (n=22). After treatment, 16 PTSD patients were remitted, and 23 had persistent PTSD based on PTSD diagnosis. The dorsal and hippocampal cingulum bundle, stria terminalis, and fornix were investigated as regions of interest. Exploratory whole-brain analyses were also performed. Groups were compared with repeated-measures ANOVA for fractional anisotropy (FA), and magnetization transfer ratio. Persistently symptomatic PTSD patients had increasing FA of the dorsal cingulum over time, and at reassessment these FA values were higher than both combat controls and the remitted PTSD group. Group-by-time interactions for FA were found in the hippocampal cingulum, fornix, and stria terminalis, posterior corona radiata, and superior longitudinal fasciculus. Our results indicate that higher FA of the dorsal cingulum bundle may be an acquired feature of persistent PTSD that develops over time. Furthermore, treatment might have differential effects on the hippocampal cingulum, fornix, stria terminalis, posterior corona radiata, and superior longitudinal fasciculus in remitted vs persistent PTSD patients. This study contributes to a better understanding of the neural underpinnings of PTSD treatment outcome.

Background

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.

Results

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.

Conclusions

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.

Epigenetic clocks are widely used aging biomarkers calculated from DNA methylation data, but this data can be surprisingly unreliable. Here we show technical noise produces deviations up to 9 years between replicates for six prominent epigenetic clocks, limiting their utility. We present a computational solution to bolster reliability, calculating principal components from CpG-level data as input for biological age prediction. Our retrained principal-component versions of six clocks show agreement between most replicates within 1.5 years, improved detection of clock associations and intervention effects, and reliable longitudinal trajectories in vivo and in vitro. This method entails only one additional step compared to traditional clocks, requires no replicates or prior knowledge of CpG reliabilities for training, and can be applied to any existing or future epigenetic biomarker. The high reliability of principal component-based clocks is critical for applications to personalized medicine, longitudinal tracking, in vitro studies, and clinical trials of aging interventions.

Epigenetic factors modify the effects of environmental factors on biological outcomes. Identification of epigenetic changes that associate with PTSD is therefore a crucial step in deciphering mechanisms of risk and resilience. In this study, our goal is to identify epigenetic signatures associated with PTSD symptom severity (PTSS) and changes in PTSS over time, using whole blood DNA methylation (DNAm) data (MethylationEPIC BeadChip) of military personnel prior to and following combat deployment. A total of 429 subjects (858 samples across 2 time points) from three male military cohorts were included in the analyses. We conducted two different meta-analyses to answer two different scientific questions: one to identify a DNAm profile of PTSS using a random effects model including both time points for each subject, and the other to identify a DNAm profile of change in PTSS conditioned on pre-deployment DNAm. Four CpGs near four genes (F2R, CNPY2, BAIAP2L1, and TBXAS1) and 88 differentially methylated regions (DMRs) were associated with PTSS. Change in PTSS after deployment was associated with 15 DMRs, of those 2 DMRs near OTUD5 and ELF4 were also associated with PTSS. Notably, three PTSS-associated CpGs near F2R, BAIAP2L1 and TBXAS1 also showed nominal evidence of association with change in PTSS. This study, which identifies PTSD-associated changes in genes involved in oxidative stress and immune system, provides novel evidence that epigenetic differences are associated with PTSS.

Mild Traumatic brain injury (mTBI) is a signature wound in military personnel, and repetitive mTBI has been linked to age-related neurogenerative disorders that affect white matter (WM) in the brain. However, findings of injury to specific WM tracts have been variable and inconsistent. This may be due to the heterogeneity of mechanisms, etiology, and comorbid disorders related to mTBI. Non-negative matrix factorization (NMF) is a data-driven approach that detects covarying patterns (components) within high-dimensional data. We applied NMF to diffusion imaging data from military Veterans with and without a self-reported TBI history. NMF identified 12 independent components derived from fractional anisotropy (FA) in a large dataset (n = 1,475) gathered through the ENIGMA (Enhancing Neuroimaging Genetics through Meta-Analysis) Military Brain Injury working group. Regressions were used to examine TBI- and mTBI-related associations in NMF-derived components while adjusting for age, sex, post-traumatic stress disorder, depression, and data acquisition site/scanner. We found significantly stronger age-dependent effects of lower FA in Veterans with TBI than Veterans without in four components (q < 0.05), which are spatially unconstrained by traditionally defined WM tracts. One component, occupying the most peripheral location, exhibited significantly stronger age-dependent differences in Veterans with mTBI. We found NMF to be powerful and effective in detecting covarying patterns of FA associated with mTBI by applying standard parametric regression modeling. Our results highlight patterns of WM alteration that are differentially affected by TBI and mTBI in younger compared to older military Veterans.

Compelling evidence suggests that epigenetic mechanisms such as DNA methylation play a role in stress regulation and in the etiologic basis of stress related disorders such as Post traumatic Stress Disorder (PTSD). Here we describe the purpose and methods of an international consortium that was developed to study the role of epigenetics in PTSD. Inspired by the approach used in the Psychiatric Genomics Consortium, we brought together investigators representing seven cohorts with a collective sample size of N = 1147 that included detailed information on trauma exposure, PTSD symptoms, and genome-wide DNA methylation data. The objective of this consortium is to increase the analytical sample size by pooling data and combining expertise so that DNA methylation patterns associated with PTSD can be identified. Several quality control and analytical pipelines were evaluated for their control of genomic inflation and technical artifacts with a joint analysis procedure established to derive comparable data over the cohorts for meta-analysis. We propose methods to deal with ancestry population stratification and type I error inflation and discuss the advantages and disadvantages of applying robust error estimates. To evaluate our pipeline, we report results from an epigenome-wide association study (EWAS) of age, which is a well-characterized phenotype with known epigenetic associations. Overall, while EWAS are highly complex and subject to similar challenges as genome-wide association studies (GWAS), we demonstrate that an epigenetic meta-analysis with a relatively modest sample size can be well-powered to identify epigenetic associations. Our pipeline can be used as a framework for consortium efforts for EWAS.

Traumatic brain injury (TBI) is common among military personnel and the civilian population and is often followed by a heterogeneous array of clinical, cognitive, behavioral, mood, and neuroimaging changes. Unlike many neurological disorders that have a characteristic abnormal central neurologic area(s) of abnormality pathognomonic to the disorder, a sufficient head impact may cause focal, multifocal, diffuse or combination of injury to the brain. This inconsistent presentation makes it difficult to establish or validate biological and imaging markers that could help improve diagnostic and prognostic accuracy in this patient population. The purpose of this manuscript is to describe both the challenges and opportunities when conducting military-relevant TBI research and introduce the Enhancing NeuroImaging Genetics through Meta-Analysis (ENIGMA) Military Brain Injury working group. ENIGMA is a worldwide consortium focused on improving replicability and analytical power through data sharing and collaboration. In this paper, we discuss challenges affecting efforts to aggregate data in this patient group. In addition, we highlight how "big data" approaches might be used to understand better the role that each of these variables might play in the imaging and functional phenotypes of TBI in Service member and Veteran populations, and how data may be used to examine important military specific issues such as return to duty, the late effects of combat-related injury, and alteration of the natural aging processes.

Traumatic brain injury (TBI) in military populations can cause disruptions in brain structure and function, along with cognitive and psychological dysfunction. Diffusion magnetic resonance imaging (dMRI) can detect alterations in white matter (WM) microstructure, but few studies have examined brain asymmetry. Examining asymmetry in large samples may increase sensitivity to detect heterogeneous areas of WM alteration in mild TBI. Through the Enhancing Neuroimaging Genetics Through Meta-Analysis Military-Relevant Brain Injury working group, we conducted a mega-analysis of neuroimaging and clinical data from 16 cohorts of Active Duty Service Members and Veterans (n = 2598). dMRI data were processed together along with harmonized demographic, injury, psychiatric, and cognitive measures. Fractional anisotropy in the cingulum showed greater asymmetry in individuals with deployment-related TBI, driven by greater left lateralization in TBI. Results remained significant after accounting for potentially confounding variables including posttraumatic stress disorder, depression, and handedness, and were driven primarily by individuals whose worst TBI occurred before age 40. Alterations in the cingulum were also associated with slower processing speed and poorer set shifting. The results indicate an enhancement of the natural left laterality of the cingulum, possibly due to vulnerability of the nondominant hemisphere or compensatory mechanisms in the dominant hemisphere. The cingulum is one of the last WM tracts to mature, reaching peak FA around 42 years old. This effect was primarily detected in individuals whose worst injury occurred before age 40, suggesting that the protracted development of the cingulum may lead to increased vulnerability to insults, such as TBI.

Epigenetic differences may help to distinguish between PTSD cases and trauma-exposed controls. Here, we describe the results of the largest DNA methylation meta-analysis of PTSD to date. Ten cohorts, military and civilian, contribute blood-derived DNA methylation data from 1,896 PTSD cases and trauma-exposed controls. Four CpG sites within the aryl-hydrocarbon receptor repressor (AHRR) associate with PTSD after adjustment for multiple comparisons, with lower DNA methylation in PTSD cases relative to controls. Although AHRR methylation is known to associate with smoking, the AHRR association with PTSD is most pronounced in non-smokers, suggesting the result was independent of smoking status. Evaluation of metabolomics data reveals that AHRR methylation associated with kynurenine levels, which are lower among subjects with PTSD. This study supports epigenetic differences in those with PTSD and suggests a role for decreased kynurenine as a contributor to immune dysregulation in PTSD.

The volume of subcortical structures represents a reliable, quantitative, and objective phenotype that captures genetic effects, environmental effects such as trauma, and disease effects such as posttraumatic stress disorder (PTSD). Trauma and PTSD represent potent exposures that may interact with genetic markers to influence brain structure and function. Genetic variants, associated with subcortical volumes in two large normative discovery samples, were used to compute polygenic scores (PGS) for the volume of seven subcortical structures. These were applied to a target sample enriched for childhood trauma and PTSD. Subcortical volume PGS from the discovery sample were strongly associated in our trauma/PTSD enriched sample (n = 7580) with respective subcortical volumes of the hippocampus (p = 1.10 × 10^-20), thalamus (p = 7.46 × 10^-10), caudate (p = 1.97 × 10^-18), putamen (p = 1.7 × 10^-12), and nucleus accumbens (p = 1.99 × 10^-7). We found a significant association between the hippocampal volume PGS and hippocampal volume in control subjects from our sample, but was absent in individuals with PTSD (GxE; (beta = -0.10, p = 0.027)). This significant GxE (PGS × PTSD) relationship persisted (p < 1 × 10^-19) in four out of five threshold peaks (0.024, 0.133, 0.487, 0.730, and 0.889) used to calculate hippocampal volume PGSs. We detected similar GxE (G × ChildTrauma) relationships in the amygdala for exposure to childhood trauma (rs4702973; p = 2.16 × 10^-7) or PTSD (rs10861272; p = 1.78 × 10^-6) in the CHST11 gene. The hippocampus and amygdala are pivotal brain structures in mediating PTSD symptomatology. Trauma exposure and PTSD modulate the effect of polygenic markers on hippocampal volume (GxE) and the amygdala volume PGS is associated with PTSD risk, which supports the role of amygdala volume as a risk factor for PTSD.