Epigenome-wide association of PTSD from heterogeneous cohorts with a common multi-site analysis pipeline

Blood-based DNA methylation and exposure risk scores predict PTSD with high accuracy in military and civilian cohorts

BACKGROUND

Incorporating genomic data into risk prediction has become an increasingly popular approach for rapid identification of individuals most at risk for complex disorders such as PTSD. Our goal was to develop and validate Methylation Risk Scores (MRS) using machine learning to distinguish individuals who have PTSD from those who do not.

METHODS

Elastic Net was used to develop three risk score models using a discovery dataset (n = 1226; 314 cases, 912 controls) comprised of 5 diverse cohorts with available blood-derived DNA methylation (DNAm) measured on the Illumina Epic BeadChip. The first risk score, exposure and methylation risk score (eMRS) used cumulative and childhood trauma exposure and DNAm variables; the second, methylation-only risk score (MoRS) was based solely on DNAm data; the third, methylation-only risk scores with adjusted exposure variables (MoRSAE) utilized DNAm data adjusted for the two exposure variables. The potential of these risk scores to predict future PTSD based on pre-deployment data was also assessed. External validation of risk scores was conducted in four independent cohorts.

RESULTS

The eMRS model showed the highest accuracy (92%), precision (91%), recall (87%), and f1-score (89%) in classifying PTSD using 3730 features. While still highly accurate, the MoRS (accuracy = 89%) using 3728 features and MoRSAE (accuracy = 84%) using 4150 features showed a decline in classification power. eMRS significantly predicted PTSD in one of the four independent cohorts, the BEAR cohort (beta = 0.6839, p=0.006), but not in the remaining three cohorts. Pre-deployment risk scores from all models (eMRS, beta = 1.92; MoRS, beta = 1.99 and MoRSAE, beta = 1.77) displayed a significant (p < 0.001) predictive power for post-deployment PTSD.

CONCLUSION

The inclusion of exposure variables adds to the predictive power of MRS. Classification-based MRS may be useful in predicting risk of future PTSD in populations with anticipated trauma exposure. As more data become available, including additional molecular, environmental, and psychosocial factors in these scores may enhance their accuracy in predicting PTSD and, relatedly, improve their performance in independent cohorts.

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

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.

Methylome-wide and meQTL analysis helps to distinguish treatment response from non-response and pathogenesis markers in schizophrenia

Schizophrenia is a complex condition with entwined genetic and epigenetic risk factors, posing a challenge to disentangle the intermixed pathological and therapeutic epigenetic signatures. To resolve this, we performed 850K methylome-wide and 700K genome-wide studies on the same set of schizophrenia patients by stratifying them into responders, non-responders, and drug-naïve patients. The key genes that signified the response were followed up using real-time gene expression studies to understand the effect of antipsychotics at the gene transcription level. The study primarily implicates hypermethylation in therapeutic response and hypomethylation in the drug-non-responsive state. Several differentially methylated sites and regions colocalized with the schizophrenia genome-wide association study (GWAS) risk genes and variants, supporting the convoluted gene-environment association. Gene ontology and protein-protein interaction (PPI) network analyses revealed distinct patterns that differentiated the treatment response from drug resistance. The study highlights the strong involvement of several processes related to nervous system development, cell adhesion, and signaling in the antipsychotic response. The ability of antipsychotic medications to alter the pathology by modulating gene expression or methylation patterns is evident from the general increase in the gene expression of response markers and histone modifiers and the decrease in class II human leukocyte antigen (HLA) genes following treatment with varying concentrations of medications like clozapine, olanzapine, risperidone, and haloperidol. The study indicates a directional overlap of methylation markers between pathogenesis and therapeutic response, thereby suggesting a careful distinction of methylation markers of pathogenesis from treatment response. In addition, there is a need to understand the trade-off between genetic and epigenetic observations. It is suggested that methylomic changes brought about by drugs need careful evaluation for their positive effects on pathogenesis, course of disease progression, symptom severity, side effects, and refractoriness.

An EWAS of dementia biomarkers and their associations with age, African ancestry, and PTSD

BACKGROUND

Large-scale cohort and epidemiological studies suggest that PTSD confers risk for dementia in later life but the biological mechanisms underlying this association remain unknown. This study examined this question by assessing the influences of PTSD, APOE ε4 genotypes, DNA methylation, and other variables on the age- and dementia-associated biomarkers Aβ40, Aβ42, GFAP, NfL, and pTau-181 measured in plasma. Our primary hypothesis was that PTSD would be associated with elevated levels of these markers.

METHODS

Analyses were based on data from a PTSD-enriched cohort of 849 individuals. We began by performing factor analyses of the biomarkers, the results of which identified a two-factor solution. Drawing from the ATN research framework, we termed the first factor, defined by Aβ40 and Aβ42, "Factor A" and the second factor, defined by GFAP, NfL and pTau-181, "Factor TN." Next, we performed epigenome-wide association analyses (EWAS) of the two-factor scores. Finally, using structural equation modeling (SEM), we evaluated (a) the influence of PTSD, age, APOE ε4 genotype and other covariates on levels of the ATN factors, and (b) tested the mediating influence of the EWAS-significant DNAm loci on these associations.

RESULTS

The Factor A EWAS identified one significant locus, cg13053408, in FANCD2OS. The Factor TN analysis identified 3 EWAS-significant associations: cg26033520 near ASCC1, cg23156469 in FAM20B, and cg15356923 in FAM19A4. The SEM showed age to be related to both factors, more so with Factor TN (β = 0.581, p < 0.001) than Factor A (β = 0.330, p < 0.001). Genotype-determined African ancestry was associated with lower Factor A (β = 0.196, p < 0.001). Contrary to our primary hypothesis, we found a modest negative bivariate correlation between PTSD and the TN factor scores (r = - 0.133, p < 0.001) attributable primarily to reduced levels of GFAP (r = - 0.128, p < 0.001).

CONCLUSIONS

This study identified novel epigenetic associations with ATN biomarkers and demonstrated robust age and ancestral associations that will be essential to consider in future efforts to develop the clinical applications of these tests. The association between PTSD and reduced GFAP, which has been reported previously, warrants further investigation.

Blood-based DNA methylation and exposure risk scores predict PTSD with high accuracy in military and civilian cohorts

Background

Incorporating genomic data into risk prediction has become an increasingly useful approach for rapid identification of individuals most at risk for complex disorders such as PTSD. Our goal was to develop and validate Methylation Risk Scores (MRS) using machine learning to distinguish individuals who have PTSD from those who do not.

Methods

Elastic Net was used to develop three risk score models using a discovery dataset (n = 1226; 314 cases, 912 controls) comprised of 5 diverse cohorts with available blood-derived DNA methylation (DNAm) measured on the Illumina Epic BeadChip. The first risk score, exposure and methylation risk score (eMRS) used cumulative and childhood trauma exposure and DNAm variables; the second, methylation-only risk score (MoRS) was based solely on DNAm data; the third, methylation-only risk scores with adjusted exposure variables (MoRSAE) utilized DNAm data adjusted for the two exposure variables. The potential of these risk scores to predict future PTSD based on pre-deployment data was also assessed. External validation of risk scores was conducted in four independent cohorts.

Results

The eMRS model showed the highest accuracy (92%), precision (91%), recall (87%), and f1-score (89%) in classifying PTSD using 3730 features. While still highly accurate, the MoRS (accuracy = 89%) using 3728 features and MoRSAE (accuracy = 84%) using 4150 features showed a decline in classification power. eMRS significantly predicted PTSD in one of the four independent cohorts, the BEAR cohort (beta = 0.6839, p-0.003), but not in the remaining three cohorts. Pre-deployment risk scores from all models (eMRS, beta = 1.92; MoRS, beta = 1.99 and MoRSAE, beta = 1.77) displayed a significant (p < 0.001) predictive power for post-deployment PTSD.

Conclusion

Results, especially those from the eMRS, reinforce earlier findings that methylation and trauma are interconnected and can be leveraged to increase the correct classification of those with vs. without PTSD. Moreover, our models can potentially be a valuable tool in predicting the future risk of developing PTSD. As more data become available, including additional molecular, environmental, and psychosocial factors in these scores may enhance their accuracy in predicting the condition and, relatedly, improve their performance in independent cohorts.

A novel principal component based method for identifying differentially methylated regions in Illumina Infinium MethylationEPIC BeadChip data

Differentially methylated regions (DMRs) are genomic regions with methylation patterns across multiple CpG sites that are associated with a phenotype. In this study, we proposed a Principal Component (PC) based DMR analysis method for use with data generated using the Illumina Infinium MethylationEPIC BeadChip (EPIC) array. We obtained methylation residuals by regressing the M-values of CpGs within a region on covariates, extracted PCs of the residuals, and then combined association information across PCs to obtain regional significance. Simulation-based genome-wide false positive (GFP) rates and true positive rates were estimated under a variety of conditions before determining the final version of our method, which we have named DMRPC. Then, DMRPC and another DMR method, coMethDMR, were used to perform epigenome-wide analyses of several phenotypes known to have multiple associated methylation loci (age, sex, and smoking) in a discovery and a replication cohort. Among regions that were analysed by both methods, DMRPC identified 50% more genome-wide significant age-associated DMRs than coMethDMR. The replication rate for the loci that were identified by only DMRPC was higher than the rate for those that were identified by only coMethDMR (90% for DMRPC vs. 76% for coMethDMR). Furthermore, DMRPC identified replicable associations in regions of moderate between-CpG correlation which are typically not analysed by coMethDMR. For the analyses of sex and smoking, the advantage of DMRPC was less clear. In conclusion, DMRPC is a new powerful DMR discovery tool that retains power in genomic regions with moderate correlation across CpGs.

The effects of epigenetic age and its acceleration on surface area, cortical thickness, and volume in young adults

DNA methylation age has been used in recent studies as an epigenetic marker of accelerated cellular aging, whose contribution to the brain structural changes was lately acknowledged. We aimed to characterize the association of epigenetic age (i.e. estimated DNA methylation age) and its acceleration with surface area, cortical thickness, and volume in healthy young adults. Using the multi-tissue method (Horvath S. DNA methylation age of human tissues and cell types. 2013. Genome Biol 14), epigenetic age was computed with saliva sample. Epigenetic age acceleration was derived from residuals after adjusting epigenetic age for chronological age. Multiple regression models were computed for 148 brain regions for surface area, cortical thickness, and volume using epigenetic age or accelerated epigenetic age as a predictor and controlling for sex. Epigenetic age was associated with surface area reduction of the left insula. It was also associated with cortical thinning and volume reduction in multiple regions, with prominent changes of cortical thickness in the left temporal regions and of volume in the bilateral orbital gyri. Finally, accelerated epigenetic age was negatively associated with right cuneus gyrus volume. Our findings suggest that understanding the mechanisms of epigenetic age acceleration in young individuals may yield valuable insights into the relationship between epigenetic aging and the cortical change and on the early development of neurocognitive pathology among young adults.

An evaluation of the genome-wide false positive rates of common methods for identifying differentially methylated regions using illumina methylation arrays

Differentially methylated regions (DMRs) are genomic regions with specific methylation patterns across multiple loci that are associated with a phenotype. We examined the genome-wide false positive (GFP) rates of five widely used DMR methods: comb-p, Bumphunter, DMRcate, mCSEA and coMethDMR using both Illumina HumanMethylation450 (450 K) and MethylationEPIC (EPIC) data and simulated continuous and dichotomous null phenotypes (i.e., generated independently of methylation data). coMethDMR provided well-controlled GFP rates (~5%) except when analysing skewed continuous phenotypes. DMRcate generally had well-controlled GFP rates when applied to 450 K data except for the skewed continuous phenotype and EPIC data only for the normally distributed continuous phenotype. GFP rates for mCSEA were at least 0.096 and comb-p yielded GFP rates above 0.34. Bumphunter had high GFP rates of at least 0.35 across conditions, reaching as high as 0.95. Analysis of the performance of these methods in specific regions of the genome found that regions with higher correlation across loci had higher regional false positive rates on average across methods. Based on the false positive rates, coMethDMR is the most recommended analysis method, and DMRcate had acceptable performance when analysing 450 K data. However, as both could display higher levels of FPs for skewed continuous distributions, a normalizing transformation of skewed continuous phenotypes is suggested. This study highlights the importance of genome-wide simulations when evaluating the performance of DMR-analysis methods.

Blood-derived deoxyribonucleic acid methylation clusters associate with adverse social exposures and endophenotypes of stress-related psychiatric illness in a trauma-exposed cohort of women

Adverse social exposures (ASEs) such as low income, low educational attainment, and childhood/adult trauma exposure are associated with variability in brain region measurements of gray matter volume (GMV), surface area (SA), and cortical thickness (CT). These CNS morphometries are associated with stress-related psychiatric illnesses and represent endophenotypes of stress-related psychiatric illness development. Epigenetic mechanisms, such as 5-methyl-cytosine (5mC), may contribute to the biological embedding of the environment but are understudied and not well understood. How 5mC relates to CNS endophenotypes of psychiatric illness is also unclear. In 97 female, African American, trauma-exposed participants from the Grady Trauma Project, we examined the associations of childhood trauma burden (CTQ), adult trauma burden, low income, and low education with blood-derived 5mC clusters and variability in brain region measurements in the amygdala, hippocampus, and frontal cortex subregions. To elucidate whether peripheral 5mC indexes central nervous system (CNS) endophenotypes of psychiatric illness, we tested whether 73 brain/blood correlated 5mC clusters, defined by networks of correlated 5mC probes measured on Illumina's HumanMethylation Epic Beadchip, mediated the relationship between ASEs and brain measurements. CTQ was negatively associated with rostral middle frontal gyrus (RMFG) SA (β =-0.231, p = 0.041). Low income and low education were also associated with SA or CT in a number of brain regions. Seven 5mC clusters were associated with CTQ (pmin = 0.002), two with low education (pmin = 0.010), and three with low income (pmin = 0.007). Two clusters fully mediated the relation between CTQ and RMFG SA, accounting for 47 and 35% of variability, respectively. These clusters were enriched for probes falling in DNA regulatory regions, as well as signal transduction and immune signaling gene ontology functions. Methylome-network analyses showed enrichment of macrophage migration (p = 9 × 10-8), T cell receptor complex (p = 6 × 10-6), and chemokine-mediated signaling (p = 7 × 10-4) pathway enrichment in association with CTQ. Our results support prior work highlighting brain region variability associated with ASEs, while informing a peripheral inflammation-based epigenetic mechanism of biological embedding of such exposures. These findings could also serve to potentiate increased investigation of understudied populations at elevated risk for stress-related psychiatric illness development.

DNA methylation of Nuclear Factor of Activated T Cells 1 mediates the prospective relation between exposure to different traumatic event types and post-traumatic stress disorder

The mechanisms through which exposure to differing trauma types become biologically embedded to shape the risk for post-traumatic stress disorder (PTSD) is unclear. DNA methylation (5-mC), particularly in stress-relevant genes, may play a role in this relationship. Here, we conducted path analysis using generalized structural equation modeling to investigate whether blood-derived 5-mC in Nuclear Factor of Activated T Cells 1 (NFATC1) mediates the prospective association between each of five different trauma types ("assaultive violence", "other injury or shocking experience", "learning of trauma to loved one", "sudden, unexpected death of a close friend or relative", and "other") and lifetime PTSD. All five trauma types were significantly associated with reduced methylation at NFATC1 CpG site, cg17057218. Two of the five trauma types were significantly associated with increased methylation at NFATC1 CpG site, cg22324981. Moreover, methylation at cg17057218 significantly mediated 21-32% of the total effect for four of the five trauma types, while methylation at cg22324981 mediated 27-40% of the total effect for two of the five trauma types. These CpG sites were differentially associated with transcription factor binding sites and chromatin state signatures. NFATC1 5-mC may be a potential mechanism in the relationship between some trauma types and prospective risk for PTSD.

Potential Societal and Cultural Implications of Transgenerational Epigenetic Methylation of Trauma and PTSD: Pathology or Resilience?

Psychological trauma is unique in that it is an environmental event that could induce biological changes and post-traumatic stress disorder (PTSD), depression, or other mood disorders in some patients. On the other hand, there may be no psychopathology (in most cases), or even sometimes post-traumatic growth and resilience. According to the DSM-5, trauma is a prerequisite for PTSD and traumatic stress disorder, but not for depressive episodes or mood disorders, or other psychiatric conditions. This paper brings attention to the preliminary literature on transgenerational inheritance due to trauma exposure and its societal and cultural implications. There is accumulating evidence that exposure to trauma can be passed transgenerationally through epigenetic inheritance leading to changes in gene expression and possible disorders or resilience. The effects of resilience from transgenerational inheritance have not been studied, but should be, for a full understanding not only of the disease risk across generations, but also of its social and cultural implications. The epigenetic pathologic effects across generations also need further studies, as the current research is preliminary; larger replications are needed for definitive and more complete understanding. I present here a glimpse of where we are, a vision of where we should go in terms of future research direction for disease risk transmission, and recommend studies of resilience and post-traumatic growth across generations, as well as other studies related to the societal implications at the population level.

Epigenetic Clock Deceleration and Maternal Reproductive Efforts: Associations With Increasing Gray Matter Volume of the Precuneus

Reproductive efforts, such as pregnancy, delivery, and interaction with children, make maternal brains optimized for child-rearing. However, extensive studies in non-human species revealed a tradeoff between reproductive effort and life expectancy. In humans, large demographic studies have shown that this is the case for the most part; however, molecular marker studies regarding aging remain controversial. There are no studies simultaneously evaluating the relationship between reproductive effort, aging, and brain structures. We therefore examined the associations between reproductive efforts (parity status, number of deliveries, motherhood period, and cumulative motherhood period), DNA methylation age (mAge) acceleration (based on Horvath’s multi-tissue clock and the skin & blood clock), and the regional gray matter volumes (obtained through brain magnetic resonance imaging (MRI) using voxel-based morphometry) in 51 mothers aged 27–46 years of children in early childhood. We found that increasing reproductive efforts were significantly associated with decelerated aging in mothers with one to four children, even after adjusting for the confounding effects in the multiple linear regression models. We also found that the left precuneus gray matter volume was larger as deceleration of aging occurred; increasing left precuneus gray matter volume, on the other hand, mediates the relationship between parity status and mAge deceleration. Our findings suggest that mothers of children in early childhood, who have had less than four children, may benefit from deceleration of aging mediated via structural changes in the precuneus.

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

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.

Epigenome-wide association study of posttraumatic stress disorder identifies novel loci in U.S. military veterans

Posttraumatic stress disorder (PTSD) is a chronic and disabling psychiatric disorder prevalent in military veterans. Epigenetic mechanisms have been implicated in the etiology of PTSD, with DNA methylation being the most studied to identify novel molecular biomarkers associated with this disorder. We performed one of the largest single-sample epigenome-wide association studies (EWAS) of PTSD to date. Our sample included 1135 male European-American U.S. veterans who participated in the National Health and Resilience in Veterans Study (NHRVS). DNA was collected from saliva samples and the Illumina HumanMethylation EPIC BeadChip was used for the methylation analysis. PTSD was assessed using the PTSD Checklist. An EWAS was conducted using linear regression adjusted for age, cell-type proportions, first 10 principal components, and smoking status. After Bonferroni correction, we identified six genome-wide significant (GWS) CpG sites associated with past-month PTSD and three CpGs with lifetime PTSD (p_range = 10^-10-10^-8). These CpG sites map to genes involved in immune function, transcription regulation, axonal guidance, cell signaling, and protein binding. Among these, SENP7, which is involved in transcription regulation and has been linked to risk-taking behavior and alcohol consumption in genome-wide association studies, replicated in an independent veteran cohort and was downregulated in medial orbitofrontal cortex of PTSD postmortem brain tissue. These findings suggest potential epigenetic biomarkers of PTSD that may help inform the pathophysiology of this disorder in veterans and other trauma-affected populations.

Association of Epigenetic Differences Screened in a Few Cases of Monozygotic Twins Discordant for Attention-Deficit Hyperactivity Disorder With Brain Structures

The present study examined the relationship between DNA methylation differences and variations in brain structures involved in the development of attention-deficit hyperactivity disorder (ADHD). First, we used monozygotic (MZ) twins discordant (2 pairs of 4 individuals, 2 boys, mean age 12.5 years) for ADHD to identify candidate DNA methylation sites involved in the development of ADHD. Next, we tried to replicate these candidates in a case-control study (ADHD: N = 18, 15 boys, mean age 10.0 years; Controls: N = 62, 40 boys, mean age 13.9 years). Finally, we examined how methylation rates at those sites relate to the degree of local structural alterations where significant differences were observed between cases and controls. As a result, we identified 61 candidate DNA methylation sites involved in ADHD development in two pairs of discordant MZ twins, among which elevated methylation at a site in the sortilin-related Vps10p domain containing receptor 2 (SorCS2) gene was replicated in the case-control study. We also observed that the ADHD group had significantly reduced gray matter volume (GMV) in the precentral and posterior orbital gyri compared to the control group and that this volume reduction was positively associated with SorCS2 methylation. Furthermore, the reduced GMV regions in children with ADHD are involved in language processing and emotional control, while SorCS2 methylation is also negatively associated with emotional behavioral problems in children. These results indicate that SorCS2 methylation might mediate a reduced GMV in the precentral and posterior orbital gyri and therefore influence the pathology of children with ADHD.

Heterogeneous Mediation Analysis on Epigenomic PTSD and Traumatic Stress in a Predominantly African American Cohort

DNA methylation (DNAm) has been suggested to play a critical role in post-traumatic stress disorder (PTSD), through mediating the relationship between trauma and PTSD. However, this underlying mechanism of PTSD for African Americans still remains unknown. To fill this gap, in this article, we investigate how DNAm mediates the effects of traumatic experiences on PTSD symptoms in the Detroit Neighborhood Health Study (DNHS) (2008-2013) which involves primarily African Americans adults. To achieve this, we develop a new mediation analysis approach for high-dimensional potential DNAm mediators. A key novelty of our method is that we consider heterogeneity in mediation effects across subpopulations. Specifically, mediators in different subpopulations could have opposite effects on the outcome, and thus could be difficult to identify under a traditional homogeneous model framework. In contrast, the proposed method can estimate heterogeneous mediation effects and identifies subpopulations in which individuals share similar effects. Simulation studies demonstrate that the proposed method outperforms existing methods for both homogeneous and heterogeneous data. We also present our mediation analysis results of a dataset with 125 participants and more than 450,000 CpG sites from the DNHS study. The proposed method finds three subgroups of subjects and identifies DNAm mediators corresponding to genes such as HSP90AA1 and NFATC1 which have been linked to PTSD symptoms in literature. Our finding could be useful in future finer-grained investigation of PTSD mechanism and in the development of new treatments for PTSD.

Trauma and posttraumatic stress disorder modulate polygenic predictors of hippocampal and amygdala volume

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.

A Comprehensive Evaluation of Cross-Omics Blood-Based Biomarkers for Neuropsychiatric Disorders

The identification of peripheral multi-omics biomarkers of brain disorders has long been hindered by insufficient sample size and confounder influence. This study aimed to compare biomarker potential for different molecules and diseases. We leveraged summary statistics of five blood quantitative trait loci studies (N = 1980 to 22,609) and genome-wide association studies (N = 9725 to 500,199) from 14 different brain disorders, such as Schizophrenia (SCZ) and Alzheimer's Disease (AD). We applied summary-based and two-sample Mendelian Randomization to estimate the associations between blood molecules and brain disorders. We identified 524 RNA, 807 methylation sites, 29 proteins, seven cytokines, and 22 metabolites having a significant association with at least one of 14 brain disorders. Simulation analyses indicated that a cross-omics combination of biomarkers had better performance for most disorders, and different disorders could associate with different omics. We identified an 11-methylation-site model for SCZ diagnosis (Area Under Curve, AUC = 0.74) by analyzing selected candidate markers in published datasets (total N = 6098). Moreover, we constructed an 18-methylation-sites model that could predict the prognosis of elders with mild cognitive impairment (hazard ratio = 2.32). We provided an association landscape between blood cross-omic biomarkers and 14 brain disorders as well as a suggestion guide for future clinical discovery and application.

A multi-modal MRI analysis of brain structure and function in relation to OXT methylation in maltreated children and adolescents

Child maltreatment dysregulates the brain's oxytocinergic system, resulting in dysfunctional attachment patterns. However, how the oxytocinergic system in children who are maltreated (CM) is epigenetically affected remains unknown. We assessed differences in salivary DNA methylation of the gene encoding oxytocin (OXT) between CM (n = 24) and non-CM (n = 31), alongside its impact on brain structures and functions using multi-modal brain imaging (voxel-based morphometry, diffusion tensor imaging, and task and resting-state functional magnetic resonance imaging). We found that CM showed higher promoter methylation than non-CM, and nine CpG sites were observed to be correlated with each other and grouped into one index (OXTmi). OXTmi was significantly negatively correlated with gray matter volume (GMV) in the left superior parietal lobule (SPL), and with right putamen activation during a rewarding task, but not with white matter structures. Using a random forest regression model, we investigated the sensitive period and type of maltreatment that contributed the most to OXTmi in CM, revealing that they were 5-8 years of age and physical abuse (PA), respectively. However, the presence of PA (PA+) was meant to reflect more severe cases, such as prolonged exposure to multiple types of abuse, than the absence of PA. PA+ was associated with significantly greater functional connectivity between the right putamen set as the seed and the left SPL and the left cerebellum exterior. The results suggest that OXT promoter hypermethylation may lead to the atypical development of reward and visual association structures and functions, thereby potentially worsening clinical aspects raised by traumatic experiences.

Posttraumatic stress disorder: from gene discovery to disease biology

Posttraumatic stress disorder (PTSD) is a complex mental disorder afflicting approximately 7% of the population. The diverse number of traumatic events and the wide array of symptom combinations leading to PTSD diagnosis contribute substantial heterogeneity to studies of the disorder. Genomic and complimentary-omic investigations have rapidly increased our understanding of the heritable risk for PTSD. In this review, we emphasize the contributions of genome-wide association, epigenome-wide association, transcriptomic, and neuroimaging studies to our understanding of PTSD etiology. We also discuss the shared risk between PTSD and other complex traits derived from studies of causal inference, co-expression, and brain morphological similarities. The investigations completed so far converge on stark contrasts in PTSD risk between sexes, partially attributed to sex-specific prevalence of traumatic experiences with high conditional risk of PTSD. To further understand PTSD biology, future studies should focus on detecting risk for PTSD while accounting for substantial cohort-level heterogeneity (e.g. civilian v. combat-exposed PTSD cases or PTSD risk among cases exposed to specific traumas), expanding ancestral diversity among study cohorts, and remaining cognizant of how these data influence social stigma associated with certain traumatic events among underrepresented minorities and/or high-risk populations.

CUE: CpG impUtation ensemble for DNA methylation levels across the human methylation450 (HM450) and EPIC (HM850) BeadChip platforms

DNA methylation at CpG dinucleotides is one of the most extensively studied epigenetic marks. With technological advancements, geneticists can profile DNA methylation with multiple reliable approaches. However, profiling platforms can differ substantially in the CpGs they assess, consequently hindering integrated analysis across platforms. Here, we present CpG impUtation Ensemble (CUE), which leverages multiple classical statistical and modern machine learning methods, to impute from the Illumina HumanMethylation450 (HM450) BeadChip to the Illumina HumanMethylationEPIC (HM850) BeadChip. Data were analysed from two population cohorts with methylation measured both by HM450 and HM850: the Extremely Low Gestational Age Newborns (ELGAN) study (n = 127, placenta) and the VA Boston Posttraumatic Stress Disorder (PTSD) genetics repository (n = 144, whole blood). Cross-validation results show that CUE achieves the lowest predicted root-mean-square error (RMSE) (0.026 in PTSD) and the highest accuracy (99.97% in PTSD) compared with five individual methods tested, including k-nearest-neighbours, logistic regression, penalized functional regression, random forest, and XGBoost. Finally, among all 339,033 HM850-only CpG sites shared between ELGAN and PTSD, CUE successfully imputed 289,604 (85.4%) sites, where success was defined as RMSE < 0.05 and accuracy >95% in PTSD. In summary, CUE is a valuable tool for imputing CpG methylation from the HM450 to HM850 platform.

Gene Expression Analysis in Three Posttraumatic Stress Disorder Cohorts Implicates Inflammation and Innate Immunity Pathways and Uncovers Shared Genetic Risk With Major Depressive Disorder

Posttraumatic stress disorder (PTSD) is a complex psychiatric disorder that can develop following exposure to traumatic events. The Psychiatric Genomics Consortium PTSD group (PGC-PTSD) has collected over 20,000 multi-ethnic PTSD cases and controls and has identified both genetic and epigenetic factors associated with PTSD risk. To further investigate biological correlates of PTSD risk, we examined three PGC-PTSD cohorts comprising 977 subjects to identify differentially expressed genes among PTSD cases and controls. Whole blood gene expression was quantified with the HumanHT-12 v4 Expression BeadChip for 726 OEF/OIF veterans from the Veterans Affairs (VA) Mental Illness Research Education and Clinical Center (MIRECC), 155 samples from the Injury and Traumatic Stress (INTRuST) Clinical Consortium, and 96 Australian Vietnam War veterans. Differential gene expression analysis was performed in each cohort separately followed by meta-analysis. In the largest cohort, we performed co-expression analysis to identify modules of genes that are associated with PTSD and MDD. We then conducted expression quantitative trait loci (eQTL) analysis and assessed the presence of eQTL interactions involving PTSD and major depressive disorder (MDD). Finally, we utilized PTSD and MDD GWAS summary statistics to identify regions that colocalize with eQTLs. Although not surpassing correction for multiple testing, the most differentially expressed genes in meta-analysis were interleukin-1 beta (IL1B), a pro-inflammatory cytokine previously associated with PTSD, and integrin-linked kinase (ILK), which is highly expressed in brain and can rescue dysregulated hippocampal neurogenesis and memory deficits. Pathway analysis revealed enrichment of toll-like receptor (TLR) and interleukin-1 receptor genes, which are integral to cellular innate immune response. Co-expression analysis identified four modules of genes associated with PTSD, two of which are also associated with MDD, demonstrating common biological pathways underlying the two conditions. Lastly, we identified four genes (UBA7, HLA-F, HSPA1B, and RERE) with high probability of a shared causal eQTL variant with PTSD and/or MDD GWAS variants, thereby providing a potential mechanism by which the GWAS variant contributes to disease risk. In summary, we provide additional evidence for genes and pathways previously reported and identified plausible novel candidates for PTSD. These data provide further insight into genetic factors and pathways involved in PTSD, as well as potential regions of pleiotropy between PTSD and MDD.

The Impact of Stress Within and Across Generations: Neuroscientific and Epigenetic Considerations

The impact of stress and trauma on biological systems in humans can be substantial. They can result in epigenetic changes, accelerated brain development and sexual maturation, and predisposition to psychopathology. Such modifications may be accompanied by behavioral, emotional, and cognitive overtones during one's lifetime. Exposure during sensitive periods of neural development may lead to long-lasting effects that may not be affected by subsequent environmental interventions. The cumulative effects of life stressors in an individual may affect offspring's methylome makeup and epigenetic clocks, neurohormonal modulation and stress reactivity, and physiological and reproductive development. While offspring may suffer deleterious effects from parental stress and their own early-life adversity, these factors may also confer traits that prove beneficial and enhance fitness to their own environment. This article synthesizes the data on how stress shapes biological and behavioral dimensions, drawing from preclinical and human models. Advances in this field of knowledge should potentially allow for an improved understanding of how interventions may be increasingly tailored according to individual biomarkers and developmental history.

Associations between childhood family emotional health, fronto-limbic grey matter volume, and saliva 5mC in young adulthood

Background

Poor family emotional health (FEH) during childhood is prevalent and impactful, and likely confers similar neurodevelopmental risks as other adverse social environments. Pointed FEH study efforts are underdeveloped, and the mechanisms by which poor FEH are biologically embedded are unclear. The current exploratory study examined whether variability in 5-methyl-cytosine (5mC) and fronto-limbic grey matter volume may represent pathways through which FEH may become biologically embedded.

Results

In 98 university students aged 18–22 years, retrospective self-reported childhood FEH was associated with right hemisphere hippocampus (b = 10.4, p = 0.005), left hemisphere amygdala (b = 5.3, p = 0.009), and right hemisphere amygdala (b = 5.8, p = 0.016) volumes. After pre-processing and filtering to 5mC probes correlated between saliva and brain, analyses showed that childhood FEH was associated with 49 5mC principal components (module eigengenes; MEs) (p_range = 3 × 10^–6 to 0.047). Saliva-derived 5mC MEs partially mediated the association between FEH and right hippocampal volume (Burlywood ME indirect effect b = − 111, p = 0.014), and fully mediated the FEH and right amygdala volume relationship (Pink4 ME indirect effect b = − 48, p = 0.026). Modules were enriched with probes falling in genes with immune, central nervous system (CNS), cellular development/differentiation, and metabolic functions.

Conclusions

Findings extend work highlighting neurodevelopmental variability associated with adverse social environment exposure during childhood by specifically implicating poor FEH, while informing a mechanism of biological embedding. FEH-associated epigenetic signatures could function as proxies of altered fronto-limbic grey matter volume associated with poor childhood FEH and inform further investigation into primarily affected tissues such as endocrine, immune, and CNS cell types.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-021-01056-y.

A systematic review of childhood maltreatment and DNA methylation: candidate gene and epigenome-wide approaches

Childhood maltreatment is a major risk factor for chronic and severe mental and physical health problems across the lifespan. Increasing evidence supports the hypothesis that maltreatment is associated with epigenetic changes that may subsequently serve as mechanisms of disease. The current review uses a systematic approach to identify and summarize the literature related to childhood maltreatment and alterations in DNA methylation in humans. A total of 100 empirical articles were identified in our systematic review of research published prior to or during March 2020, including studies that focused on candidate genes and studies that leveraged epigenome-wide data in both children and adults. Themes arising from the literature, including consistent and inconsistent patterns of results, are presented. Several directions for future research, including important methodological considerations for future study design, are discussed. Taken together, the literature on childhood maltreatment and DNA methylation underscores the complexity of transactions between the environment and biology across development.

PTSD is associated with increased DNA methylation across regions of HLA-DPB1 and SPATC1L

Posttraumatic stress disorder (PTSD) is characterized by intrusive thoughts, avoidance, negative alterations in cognitions and mood, and arousal symptoms that adversely affect mental and physical health. Recent evidence links changes in DNA methylation of CpG cites to PTSD. Since clusters of proximal CpGs share similar methylation signatures, identification of PTSD-associated differentially methylated regions (DMRs) may elucidate the pathways defining differential risk and resilience of PTSD. Here we aimed to identify epigenetic differences associated with PTSD. DNA methylation data profiled from blood samples using the MethylationEPIC BeadChip were used to perform a DMR analysis in 187 PTSD cases and 367 trauma-exposed controls from the Grady Trauma Project (GTP). DMRs were assessed with R package bumphunter. We identified two regions that associate with PTSD after multiple test correction. These regions were in the gene body of HLA-DPB1 and in the promoter of SPATC1L. The DMR in HLA-DPB1 was associated with PTSD in an independent cohort. Both DMRs included CpGs whose methylation associated with nearby sequence variation (meQTL) and that associated with expression of their respective genes (eQTM). This study supports an emerging literature linking PTSD risk to genetic and epigenetic variation in the HLA region.

RETRACTED: Epigenetic modification of the oxytocin gene is associated with gray matter volume and trait empathy in mothers

This article has been retracted: please see Elsevier Policy on Article Withdrawal (https://www.elsevier.com/about/our-business/policies/article-withdrawal). The Editor has concluded that the acceptance of this article was partly based upon the positive advice of an unreliable reviewer report. The report was provided to the journal by a reviewer suggested by the authors, and there were inappropriate communications between the authors and reviewer during the peer-review process. The Editor has therefore concluded that the review was not appropriate or independent. This manipulation of the peer-review process represents a clear violation of the fundamentals of peer review, our publishing policies, and publishing ethics standards. Apologies are offered to the readers of the journal that this deception was not detected during the submission process.

Epigenome-wide meta-analysis of PTSD across 10 military and civilian cohorts identifies methylation changes in AHRR

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.

PTSD and the klotho longevity gene: Evaluation of longitudinal effects on inflammation via DNA methylation

BACKGROUND

Longevity gene klotho (KL) is associated with age-related phenotypes including lifespan, cardiometabolic disorders, cognition, and brain morphology, in part, by conferring protection against inflammation. We hypothesized that the KL/inflammation association might be altered in the presence of psychiatric stress and operate via epigenetic pathways. We examined KL polymorphisms, and their interaction with posttraumatic stress disorder (PTSD) symptoms, in association with KL DNA methylation in blood. We further examined KL DNA methylation as a predictor of longitudinal changes in a peripheral biomarker of inflammation (C-reactive protein; CRP).

METHODS

The sample comprised 309 white non-Hispanic military veterans (93.5 % male; mean age: 32 years, range: 19-65; 30 % PTSD per structured diagnostic interview); 111 were reassessed approximately two years later.

RESULTS

Analyses revealed a methylation quantitative trait locus at rs9527025 (C370S, previously implicated in numerous studies of aging) in association with a Cytosine-phosphate-Guanine site (cg00129557; B = -.65, p =  1.29 X 10^-20), located within a DNase hypersensitivity site in the body of KL. There was also a rs9527025 x PTSD severity interaction (B = .004, p = .035) on methylation at this locus such that the minor allele was associated with reduced cg00129557 methylation in individuals with few or no PTSD symptoms while this effect was attenuated in those with elevated levels of PTSD. Path models revealed that methylation at cg00129557 was inversely associated with CRP over time (B = -.14, p = .005), controlling for baseline CRP. There was also an indirect effect of rs9527025 X PTSD on subsequent CRP via cg00129557 methylation (indirect B = -.002, p = .033).

CONCLUSIONS

Results contribute to our understanding of the epigenetic correlates of inflammation in PTSD and suggest that KL methylation may be a mechanism by which KL genotype confers risk vs. resilience to accelerated aging in those experiencing traumatic stress.

The relationship between DNA methylation in neurotrophic genes and age as evidenced from three independent cohorts: differences by delirium status

We previously reported the association between DNA methylation (DNAm) of pro-inflammatory cytokine genes and age. In addition, neurotrophic factors are known to be associated with age and neurocognitive disorders. Therefore, we hypothesized that DNAm of neurotrophic genes change with age, especially in delirium patients. DNAm was analyzed using the Illumina HumanMethylation450 or HumanMethylationEPIC BeadChip Kit in 3 independent cohorts: blood from 383 Grady Trauma Project subjects, brain from 21 neurosurgery patients, and blood from 87 inpatients with and without delirium. Both blood and brain samples showed that most of the DNAm of neurotrophic genes were positively correlated with age. Furthermore, DNAm of neurotrophic genes was more positively correlated with age in delirium cases than in non-delirium controls. These findings support our hypothesis that the neurotrophic genes may be epigenetically modulated with age, and this process may be contributing to the pathophysiology of delirium.

An epigenome-wide association study of posttraumatic stress disorder in US veterans implicates several new DNA methylation loci

BACKGROUND

Previous studies using candidate gene and genome-wide approaches have identified epigenetic changes in DNA methylation (DNAm) associated with posttraumatic stress disorder (PTSD).

METHODS

In this study, we performed an EWAS of PTSD in a cohort of Veterans (n = 378 lifetime PTSD cases and 135 controls) from the Translational Research Center for TBI and Stress Disorders (TRACTS) cohort assessed using the Illumina EPIC Methylation BeadChip which assesses DNAm at more than 850,000 sites throughout the genome. Our model included covariates for ancestry, cell heterogeneity, sex, age, and a smoking score based on DNAm at 39 smoking-associated CpGs. We also examined in EPIC-based DNAm data generated from pre-frontal cortex (PFC) tissue from the National PTSD Brain Bank (n = 72).

RESULTS

The analysis of blood samples yielded one genome-wide significant association with PTSD at cg19534438 in the gene G0S2 (p = 1.19 × 10-7, padj = 0.048). This association was replicated in an independent PGC-PTSD-EWAS consortium meta-analysis of military cohorts (p = 0.0024). We also observed association with the smoking-related locus cg05575921 in AHRR despite inclusion of a methylation-based smoking score covariate (p = 9.16 × 10-6), which replicates a previously observed PGC-PTSD-EWAS association (Smith et al. 2019), and yields evidence consistent with a smoking-independent effect. The top 100 EWAS loci were then examined in the PFC data. One of the blood-based PTSD loci, cg04130728 in CHST11, which was in the top 10 loci in blood, but which was not genome-wide significant, was significantly associated with PTSD in brain tissue (in blood p = 1.19 × 10-5, padj = 0.60, in brain, p = 0.00032 with the same direction of effect). Gene set enrichment analysis of the top 500 EWAS loci yielded several significant overlapping GO terms involved in pathogen response, including "Response to lipopolysaccharide" (p = 6.97 × 10-6, padj = 0.042).

CONCLUSIONS

The cross replication observed in independent cohorts is evidence that DNA methylation in peripheral tissue can yield consistent and replicable PTSD associations, and our results also suggest that that some PTSD associations observed in peripheral tissue may mirror associations in the brain.

Epigenetic predictors of all-cause mortality are associated with objective measures of neighborhood disadvantage in an urban population

BACKGROUND

Neighborhood characteristics are robust predictors of overall health and mortality risk for residents. Though there has been some investigation of the role that molecular indicators may play in mediating neighborhood exposures, there has been little effort to incorporate newly developed epigenetic biomarkers into our understanding of neighborhood characteristics and health outcomes.

METHODS

Using 157 participants of the Detroit Neighborhood Health Study with detailed assessments of neighborhood characteristics and genome-wide DNA methylation profiling via the Illumina 450K methylation array, we assessed the relationship between objective neighborhood characteristics and a validated DNA methylation-based epigenetic mortality risk score (eMRS). Associations were adjusted for age, race, sex, ever smoking, ever alcohol usage, education, years spent in neighborhood, and employment. A secondary model additionally adjusted for personal neighborhood perception. We summarized 19 neighborhood quality indicators assessed for participants into 9 principal components which explained over 90% of the variance in the data and served as metrics of objective neighborhood quality exposures.

RESULTS

Of the nine principal components utilized for this study, one was strongly associated with the eMRS (β = 0.15; 95% confidence interval = 0.06-0.24; P = 0.002). This principal component (PC7) was most strongly driven by the presence of abandoned cars, poor streets, and non-art graffiti. Models including both PC7 and individual indicators of neighborhood perception indicated that only PC7 and not neighborhood perception impacted the eMRS. When stratified on neighborhood indicators of greenspace, we observed a potentially protective effect of large mature trees as this feature substantially attenuated the observed association.

CONCLUSION

Objective measures of neighborhood disadvantage are significantly associated with an epigenetic predictor of mortality risk, presenting a potential novel avenue by which neighborhood-level exposures may impact health. Associations were independent of an individual's perception of their neighborhood and attenuated by neighborhood greenspace features. More work should be done to determine molecular risk factors associated with neighborhoods, and potentially protective neighborhood features against adverse molecular effects.

Longitudinal epigenome-wide association studies of three male military cohorts reveal multiple CpG sites associated with post-traumatic stress disorder

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.

Available Software for Meta-analyses of Genome-wide Expression Studies

Advances in transcriptomic methods have led to a large number of published Genome-Wide Expression Studies (GWES), in humans and model organisms. For several years, GWES involved the use of microarray platforms to compare genome-expression data for two or more groups of samples of interest. Meta-analysis of GWES is a powerful approach for the identification of differentially expressed genes in biological topics or diseases of interest, combining information from multiple primary studies. In this article, the main features of available software for carrying out meta-analysis of GWES have been reviewed and seven packages from the Bioconductor platform and five packages from the CRAN platform have been described. In addition, nine previously described programs and four online programs are reviewed. Finally, advantages and disadvantages of these available programs and proposed key points for future developments have been discussed.

A review of epigenetic contributions to post-traumatic stress disorder

Post-traumatic stress disorder (PTSD) is a syndrome which serves as a classic example of psychiatric disorders that result from the intersection of nature and nurture, or gene and environment. By definition, PTSD requires the experience of a traumatic exposure, and yet data suggest that the risk for PTSD in the aftermath of trauma also has a heritable (genetic) component. Thus, PTSD appears to require both a biological (genetic) predisposition that differentially alters how the individual responds to or recovers from trauma exposure. Epigenetics is defined as the study of changes in organisms caused by modification of gene expression rather than alteration of the genetic code itself, and more recently it has come to refer to direct alteration of DNA regulation, but without altering the primary sequence of DNA, or the genetic code. With regards to PTSD, epigenetics provides one way for environmental exposure to be "written" upon the genome, as a direct result of gene and environment (trauma) interactions. This review provides an overview of the main currently understood types of epigenetic regulation, including DNA methylation, histone regulation of chromatin, and noncoding RNA regulation of gene expression. Furthermore, we examine recent literature related to how these methods of epigenetic regulation may be involved in differential risk and resilience for PTSD in the aftermath of trauma.
.

Methylomic profiles reveal sex-specific differences in leukocyte composition associated with post-traumatic stress disorder

Post-traumatic stress disorder (PTSD) is a debilitating mental disorder precipitated by trauma exposure. However, only some persons exposed to trauma develop PTSD. There are sex differences in risk; twice as many women as men develop a lifetime diagnosis of PTSD. Methylomic profiles derived from peripheral blood are well-suited for investigating PTSD because DNA methylation (DNAm) encodes individual response to trauma and may play a key role in the immune dysregulation characteristic of PTSD pathophysiology. In the current study, we leveraged recent methodological advances to investigate sex-specific differences in DNAm-based leukocyte composition that are associated with lifetime PTSD. We estimated leukocyte composition on a combined methylation array dataset (483 participants, ∼450 k CpG sites) consisting of two civilian cohorts, the Detroit Neighborhood Health Study and Grady Trauma Project. Sex-stratified Mann-Whitney U test and two-way ANCOVA revealed that lifetime PTSD was associated with significantly higher monocyte proportions in males, but not in females (Holm-adjusted p-val < 0.05). No difference in monocyte proportions was observed between current and remitted PTSD cases in males, suggesting that this sex-specific difference may reflect a long-standing trait of lifetime history of PTSD, rather than current state of PTSD. Associations with lifetime PTSD or PTSD status were not observed in any other leukocyte subtype and our finding in monocytes was confirmed using cell estimates based on a different deconvolution algorithm, suggesting that our sex-specific findings are robust across cell estimation approaches. Overall, our main finding of elevated monocyte proportions in males, but not in females with lifetime history of PTSD provides evidence for a sex-specific difference in peripheral blood leukocyte composition that is detectable in methylomic profiles and that may reflect long-standing changes associated with PTSD diagnosis.

Translating Molecular and Neuroendocrine Findings in Posttraumatic Stress Disorder and Resilience to Novel Therapies

Many biological systems are altered in association with posttraumatic stress disorder (PTSD) and resilience. However, there are only few approved pharmacological treatments for PTSD, and no approved medications to enhance resilience. This article provides a critical review of select neurobiological findings in PTSD and resilience, and also of pharmacologic approaches that have emerged from this work. The medications summarized involve engagement with targets in the adrenergic, hypothalamic-pituitary-adrenal axis, and neuropeptide Y systems. Other highlighted approaches involve the use of ketamine and 3,4-methylenedioxymethamphetamine-assisted psychotherapy, which recently surfaced as promising strategies for PTSD, though the neurobiological mechanisms underlying their actions, including for promoting resilience, are not yet fully understood. The former approaches fall within the broad concept of "rational pharmacotherapy," in that they attempt to directly target dysregulated systems known to be associated with posttraumatic symptoms. To the extent that use of ketamine and 3,4-methylenedioxymethamphetamine promotes symptom improvement and resilience in PTSD, this provides an opportunity for reverse translation and identification of relevant targets and mechanisms of action through careful study of biological changes resulting from these interventions. Promoting resilience in trauma-exposed individuals may involve more than pharmacologically manipulating dysregulated molecules and pathways associated with developing and sustaining PTSD symptom severity, but also producing a substantial change in mental state that increases the ability to engage with traumatic material in psychotherapy. Neurobiological examination in the context of treatment studies may yield novel targets and promote a greater understanding of mechanisms of recovery from trauma.

Accelerating research on biological aging and mental health: Current challenges and future directions

Aging is associated with complex biological changes that can be accelerated, slowed, or even temporarily reversed by biological and non-biological factors. This article focuses on the link between biological aging, psychological stressors, and mental illness. Rather than comprehensively reviewing this rapidly expanding field, we highlight challenges in this area of research and propose potential strategies to accelerate progress in this field. This effort requires the interaction of scientists across disciplines - including biology, psychiatry, psychology, and epidemiology; and across levels of analysis that emphasize different outcome measures - functional capacity, physiological, cellular, and molecular. Dialogues across disciplines and levels of analysis naturally lead to new opportunities for discovery but also to stimulating challenges. Some important challenges consist of 1) establishing the best objective and predictive biological age indicators or combinations of indicators, 2) identifying the basis for inter-individual differences in the rate of biological aging, and 3) examining to what extent interventions can delay, halt or temporarily reverse aging trajectories. Discovering how psychological states influence biological aging, and vice versa, has the potential to create novel and exciting opportunities for healthcare and possibly yield insights into the fundamental mechanisms that drive human aging.

Correction for multiple testing in candidate-gene methylation studies

Aim: We compared the performance of multiple testing corrections for candidate gene methylation studies, namely Sidak (accurate Bonferroni), false-discovery rate and three adjustments that incorporate the correlation between CpGs: extreme tail theory (ETT), Gao et al. (GEA), and Li and Ji methods. Materials & methods: The experiment-wide type 1 error rate was examined in simulations based on Illumina EPIC and 450K data. Results: For high-correlation genes, Sidak and false-discovery rate corrections were conservative while the Li and Ji method was liberal. The GEA method tended to be conservative unless a threshold parameter was adjusted. The ETT yielded an appropriate type 1 error rate. Conclusion: For genes with substantial correlation across measured CpGs, GEA and ETT can appropriately correct for multiple testing in candidate gene methylation studies.

Investigation of bidirectional longitudinal associations between advanced epigenetic age and peripheral biomarkers of inflammation and metabolic syndrome

Epigenetic age estimations based on DNA methylation (DNAm) can predict human chronological age with a high level of accuracy. These DNAm age algorithms can also be used to index advanced cellular age, when estimated DNAm age exceeds chronological age. Advanced DNAm age has been associated with several diseases and metabolic and inflammatory pathology, but the causal direction of this association is unclear. The goal of this study was to examine potential bidirectional associations between advanced epigenetic age and metabolic and inflammatory markers over time in a longitudinal cohort of 179 veterans with a high prevalence of posttraumatic stress disorder (PTSD) who were assessed over the course of two years. Analyses focused on two commonly investigated metrics of advanced DNAm age derived from the Horvath (developed across multiple tissue types) and Hannum (developed in whole blood) DNAm age algorithms. Results of cross-lagged panel models revealed that advanced Hannum DNAm age at Time 1 (T1) was associated with increased (i.e., accounting for T1 levels) metabolic syndrome (MetS) severity at Time 2 (T2; p = < 0.001). This association was specific to worsening lipid panels and indicators of abdominal obesity (p = 0.001). In contrast, no baseline measures of inflammation or metabolic pathology were associated with changes in advanced epigenetic age over time. No associations emerged between advanced Horvath DNAm age and any of the examined biological parameters. Results suggest that advanced epigenetic age, when measured using an algorithm developed in whole blood, may be a prognostic marker of pathological metabolic processes. This carries implications for understanding pathways linking advanced epigenetic age to morbidity and mortality.

Posttraumatic psychopathology and the pace of the epigenetic clock: a longitudinal investigation

BACKGROUND

Posttraumatic stress disorder (PTSD) and stress/trauma exposure are cross-sectionally associated with advanced DNA methylation age relative to chronological age. However, longitudinal inquiry and examination of associations between advanced DNA methylation age and a broader range of psychiatric disorders is lacking. The aim of this study was to examine if PTSD, depression, generalized anxiety, and alcohol-use disorders predicted acceleration of DNA methylation age over time (i.e. an increasing pace, or rate of advancement, of the epigenetic clock).

METHODS

Genome-wide DNA methylation and a comprehensive set of psychiatric symptoms and diagnoses were assessed in 179 Iraq/Afghanistan war veterans who completed two assessments over the course of approximately 2 years. Two DNA methylation age indices (Horvath and Hannum), each a weighted index of an array of genome-wide DNA methylation probes, were quantified. The pace of the epigenetic clock was operationalized as change in DNA methylation age as a function of time between assessments.

RESULTS

Analyses revealed that alcohol-use disorders (p = 0.001) and PTSD avoidance and numbing symptoms (p = 0.02) at Time 1 were associated with an increasing pace of the epigenetic clock over time, per the Horvath (but not the Hannum) index of cellular aging.

CONCLUSIONS

This is the first study to suggest that posttraumatic psychopathology is longitudinally associated with a quickened pace of the epigenetic clock. Results raise the possibility that accelerated cellular aging is a common biological consequence of stress-related psychopathology, which carries implications for identifying mechanisms of stress-related cellular aging and developing interventions to slow its pace.

Understanding Resilience and Preventing and Treating PTSD

Resilience is defined as the dynamic ability to adapt successfully in the face of adversity, trauma, or significant threat. Some of the key early studies of resilience were observational studies in children. They were followed by research in adults, studies testing interventions to promote resilience in different populations, and a recent upsurge of studies on the underlying genomic and neurobiological mechanisms. Neural and molecular studies in preclinical models of resilience are also increasingly identifying active stress adaptations in resilient animals. Knowledge gained from animal and human studies of resilience can be harnessed to develop new preventive interventions to enhance resilience in at-risk populations. Further, treatment interventions focused on enhancing potentially modifiable protective factors that are consistently linked to psychological resilience can enrich currently available treatment interventions for individuals with posttraumatic stress disorder (PTSD). Translating our expanding knowledge of the neurobiology of resilience additionally promises to yield novel therapeutic strategies for treating this disabling condition. This review summarizes the vast field of resilience research spanning genomic, psychosocial, and neurobiological levels, and discusses how findings have led and can lead to new preventive and treatment interventions for PTSD.

DNA methylation correlates of PTSD: Recent findings and technical challenges

There is increasing evidence that epigenetic factors play a critical role in posttraumatic stress disorder (PTSD), by mediating the impact of environmental exposures to trauma on the regulation of gene expression. DNA methylation is one epigenetic process that has been highly studied in PTSD. This review will begin by providing an overview of DNA methylation (DNAm) methods, and will then highlight two major biological systems that have been identified in the epigenetic regulation in PTSD: (a) the immune system and (b) the stress response system. In addition to candidate gene approaches, we will review novel strategies to study epigenome-wide PTSD-related effects, including epigenome-wide algorithms that distill information from many loci into a single summary score (e.g., measures of "epigenetic age" which have been associated with PTSD). This review will also cover recent epigenome wide association studies (EWAS) of PTSD, and biological pathway models used to identify gene sets enriched in PTSD. Finally, we address technical and methodological advances and challenges to the field, and highlight exciting directions for future research.

Resilience Against Traumatic Stress: Current Developments and Future Directions

Given the high prevalence of stress-related mental disorders, their impact on person, family, and society and the paucity of treatment options for most of these disorders, there is currently a pressing need for innovative approaches to deal with these issues and enhance well-being. One approach which has received increasing attention over the last decade is to shift our scientific and clinical focus from risk factors for psychopathology to factors promoting resilience and mental well-being. In order to summarize and evaluate the current state of scientific affairs on the biological basis of resilience, we provide an overview of the literature on animal and human studies of resilience. Because resilience can only truly be operationalized through longitudinal data collection and analyses, we focus primarily on longitudinal studies. This review shows that the concept of resilience is currently being operationalized, measured and even defined in widely variable manners, both within animal and human studies. We further provide an overview of existing and new strategies that could help promote resilience and which are proposed to be implemented more often in clinical situations. Finally, we summarize the challenges the field is facing and provide recommendations for future research.

Epigenetic meta-analysis across three civilian cohorts identifies NRG1 and HGS as blood-based biomarkers for post-traumatic stress disorder

AIM

Trauma exposure is a necessary, but not deterministic, contributor to post-traumatic stress disorder (PTSD). Epigenetic factors may distinguish between trauma-exposed individuals with versus without PTSD.

MATERIALS & METHODS

We conducted a meta-analysis of PTSD epigenome-wide association studies in trauma-exposed cohorts drawn from civilian contexts. Whole blood-derived DNA methylation levels were analyzed in 545 study participants, drawn from the three civilian cohorts participating in the PTSD working group of the Psychiatric Genomics Consortium.

RESULTS

Two CpG sites significantly associated with current PTSD in NRG1 (cg23637605) and in HGS (cg19577098).

CONCLUSION

PTSD is associated with differential methylation, measured in blood, within HGS and NRG1 across three civilian cohorts.

Epigenetics of Delirium and Aging: Potential Role of DNA Methylation Change on Cytokine Genes in Glia and Blood Along With Aging

Background: Delirium in elderly patients is common and dangerous. Major risk factors include aging and exogenous insults, such as infection or surgery. In animal models, aging enhances pro-inflammatory cytokine release from microglia in response to exogenous insults. The epigenetic mechanism DNA methylation (DNAm) regulates gene expression and changes with age. Older individuals may have methylation changes that influence the increased cytokine upon insult, but the degree to which aging affects DNAm of cytokine genes is not fully understood. Methods: The relationship between DNAm and aging of pro-inflammatory cytokine genes (TNF-alpha, IL1-beta, IL-6) was investigated using methylation array data in two cohorts. Brain and blood samples were collected from a neurosurgery cohort (NSG) of 21 subjects who underwent brain resection. A second cohort, the Grady Trauma Project (GTP), included blood samples from 265 subjects. Results: In the NSG cohort, a significant negative correlation between age and DNAm in brain was found at a CpG in IL-6. With the GTP dataset, significant negative correlations between age and DNAm were seen at most of the CpGs in TNF-alpha. Also, TNF-Alpha expression increases with age. These GTP DNAm correlations were also nominally significant in NSG blood samples. In neuronal negative NSG brain tissue, a similar negative trend was observed. Conclusions: With aging, a decrease in DNAm of cytokines gene CpGs in glia and blood was seen. As this can affect their expression, additional research is needed to fully elucidate the role of DNAm in aging and how it may influence the pathogenesis of delirium.

Robust Findings From 25 Years of PTSD Genetics Research

PURPOSE OF REVIEW

The purpose of this review is to contextualize findings from the first 25 years of PTSD genetics research, focusing on the most robust findings and interpreting results in light of principles that have emerged from modern genetics studies.

RECENT FINDINGS

Genome-wide association studies (GWAS) encompassing tens of thousands of participants enabled the first molecular genetic heritability and genetic correlation estimates for PTSD in 2017. In 2018, highly promising loci for PTSD were reported, including variants in and near the CAMKV, KANSL1, and TCF4 genes. Twin studies from 25 years ago established that PTSD is genetically influenced and foreshadowed the molecular genetic findings of today. Discoveries that were impossible with smaller studies have been achieved via collaborative/team-science efforts. Most promisingly, individual genomic loci offer entirely novel clues about PTSD etiology, providing the raw material for transformative discoveries, and the future of PTSD research is bright.

Traumatic stress and accelerated DNA methylation age: A meta-analysis

BACKGROUND

Recent studies examining the association between posttraumatic stress disorder (PTSD) and accelerated aging, as defined by DNA methylation-based estimates of cellular age that exceed chronological age, have yielded mixed results.

METHODS

We conducted a meta-analysis of trauma exposure and PTSD diagnosis and symptom severity in association with accelerated DNA methylation age using data from 9 cohorts contributing to the Psychiatric Genomics Consortium PTSD Epigenetics Workgroup (combined N = 2186). Associations between demographic and cellular variables and accelerated DNA methylation age were also examined, as was the moderating influence of demographic variables.

RESULTS

Meta-analysis of regression coefficients from contributing cohorts revealed that childhood trauma exposure (when measured with the Childhood Trauma Questionnaire) and lifetime PTSD severity evidenced significant, albeit small, meta-analytic associations with accelerated DNA methylation age (ps = 0.028 and 0.016, respectively). Sex, CD4T cell proportions, and natural killer cell proportions were also significantly associated with accelerated DNA methylation age (all ps < 0.02). PTSD diagnosis and lifetime trauma exposure were not associated with advanced DNA methylation age. There was no evidence of moderation of the trauma or PTSD variables by demographic factors.

CONCLUSIONS

Results suggest that traumatic stress is associated with advanced epigenetic age and raise the possibility that cells integral to immune system maintenance and responsivity play a role in this. This study highlights the need for additional research into the biological mechanisms linking traumatic stress to accelerated DNA methylation age and the importance of furthering our understanding of the neurobiological and health consequences of PTSD.