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.

Multi-ancestry genome-wide association study of major depression aids locus discovery, fine mapping, gene prioritization and causal inference.

Most genome-wide association studies (GWAS) of major depression (MD) have been conducted in samples of European ancestry. Here we report a multi-ancestry GWAS of MD, adding data from 21 cohorts with 88,316 MD cases and 902,757 controls to previously reported data. This analysis used a range of measures to define MD and included samples of African (36% of effective sample size), East Asian (26%) and South Asian (6%) ancestry and Hispanic/Latin American participants (32%). The multi-ancestry GWAS identified 53 significantly associated novel loci. For loci from GWAS in European ancestry samples, fewer than expected were transferable to other ancestry groups. Fine mapping benefited from additional sample diversity. A transcriptome-wide association study identified 205 significantly associated novel genes. These findings suggest that, for MD, increasing ancestral and global diversity in genetic studies may be particularly important to ensure discovery of core genes and inform about transferability of findings.

A Mendelian randomization study of genetic liability to post-traumatic stress disorder and risk of ischemic stroke.

Observational studies have shown an association between post-traumatic stress disorder (PTSD) and ischemic stroke (IS) but given the susceptibility to confounding it is unclear if these associations represent causal effects. Mendelian randomization (MR) facilitates causal inference that is robust to the influence of confounding. Using two sample MR, we investigated the causal effect of genetic liability to PTSD on IS risk. Ancestry-specific genetic instruments of PTSD and four quantitative sub-phenotypes of PTSD, including hyperarousal, avoidance, re-experiencing, and total symptom severity score (PCL-Total) were obtained from the Million Veteran Programme (MVP) using a threshold P value (P) of <5 × 10-7, clumping distance of 1000 kilobase (Mb) and r2 < 0.01. Genetic association estimates for IS were obtained from the MEGASTROKE consortium (Ncases = 34,217, Ncontrols = 406,111) for European ancestry individuals and from the Consortium of Minority Population Genome-Wide Association Studies of Stroke (COMPASS) (Ncases = 3734, Ncontrols = 18,317) for African ancestry individuals. We used the inverse-variance weighted (IVW) approach as the main analysis and performed MR-Egger and the weighted median methods as pleiotropy-robust sensitivity analyses. In European ancestry individuals, we found evidence of an association between genetic liability to PTSD avoidance, and PCL-Total and increased IS risk (odds ratio (OR)1.04, 95% Confidence Interval (CI) 1.007-1.077, P = 0.017 for avoidance and (OR 1.02, 95% CI 1.010-1.040, P = 7.6 × 10-4 for PCL total). In African ancestry individuals, we found evidence of an association between genetically liability to PCL-Total and reduced IS risk (OR 0.95 (95% CI 0.923-0.991, P = 0.01) and hyperarousal (OR 0.83 (95% CI 0.691-0.991, P = 0.039) but no association was observed for PTSD case-control, avoidance, or re-experiencing. Similar estimates were obtained with MR sensitivity analyses. Our findings suggest that specific sub-phenotypes of PTSD, such as hyperarousal, avoidance, PCL total, may have a causal effect on people of European and African ancestry's risk of IS. This shows that the molecular mechanisms behind the relationship between IS and PTSD may be connected to symptoms of hyperarousal and avoidance. To clarify the precise biological mechanisms involved and how they may vary between populations, more research is required.

Differentially Expressed Heterogeneous Overdispersion Genes Testing for Count Data.

The mRNA-seq data analysis is a powerful technology for inferring information from biological systems of interest. Specifically, the sequenced RNA fragments are aligned with genomic reference sequences, and we count the number of sequence fragments corresponding to each gene for each condition. A gene is identified as differentially expressed (DE) if the difference in its count numbers between conditions is statistically significant. Several statistical analysis methods have been developed to detect DE genes based on RNA-seq data. However, the existing methods could suffer decreasing power to identify DE genes arising from overdispersion and limited sample size. We propose a new differential expression analysis procedure: heterogeneous overdispersion genes testing (DEHOGT) based on heterogeneous overdispersion modeling and a post-hoc inference procedure. DEHOGT integrates sample information from all conditions and provides a more flexible and adaptive overdispersion modeling for the RNA-seq read count. DEHOGT adopts a gene-wise estimation scheme to enhance the detection power of differentially expressed genes. DEHOGT is tested on the synthetic RNA-seq read count data and outperforms two popular existing methods, DESeq and EdgeR, in detecting DE genes. We apply the proposed method to a test dataset using RNAseq data from microglial cells. DEHOGT tends to detect more differently expressed genes potentially related to microglial cells under different stress hormones treatments.

Serum lead, mercury, manganese, and copper and DNA methylation age among adults in Detroit, Michigan.

Although the effects of lead, mercury, manganese, and copper on individual disease processes are well understood, estimating the health effects of long-term exposure to these metals at the low concentrations often observed in the general population is difficult. In addition, the health effects of joint exposure to multiple metals are difficult to estimate. Biological aging refers to the integrative progression of multiple physiologic and molecular changes that make individuals more at risk of disease. Biomarkers of biological aging may be useful to estimate the population-level effects of metal exposure prior to the development of disease in the population. We used data from 290 participants in the Detroit Neighborhood Health Study to estimate the effect of serum lead, mercury, manganese, and copper on three DNA methylation-based biomarkers of biological aging (Horvath Age, PhenoAge, and GrimAge). We used mixed models and Bayesian kernel machine regression and controlled for participant sex, race, ethnicity, cigarette use, income, educational attainment, and block group poverty. We observed consistently positive estimates of the effects between lead and GrimAge acceleration and mercury and PhenoAge acceleration. In contrast, we observed consistently negative associations between manganese and PhenoAge acceleration and mercury and Horvath Age acceleration. We also observed curvilinear relationships between copper and both PhenoAge and GrimAge acceleration. Increasing total exposure to the observed mixture of metals was associated with increased PhenoAge and GrimAge acceleration and decreased Horvath Age acceleration. These findings indicate that an increase in serum lead or mercury from the 25th to 75th percentile is associated with a ∼0.25-year increase in two epigenetic markers of all-cause mortality in a population of adults in Detroit, Michigan. While few of the findings were statistically significant, their consistency and novelty warrant interest.

Community engagement in epigenomic and neurocognitive research on post-traumatic stress disorder in Rwandans exposed to the 1994 genocide against the Tutsi: lessons learned.

Epigenomic and neurocognitive studies have provided new perspectives on post-traumatic stress disorder and its intergenerational transmission. This article outlines the lessons learned from community engagement (CE) in such research on Rwandan genocide survivors. A strong trauma-related response was observed within the research project-targeted community (genocide survivors) during explanation of the project. CE also revealed privacy concerns, as community members worried that any leakage of genetic/(epi)genomic data could affect not only themselves but also their close relatives. Adopting a culture of CE in the process of research implementation enables the prioritization of targeted community needs and interests. Furthermore, CE has stimulated the development of mental healthcare interventions, which married couples can apply to protect their offspring and thus truly break the cycle of inherited vulnerability.

Studies of how human genes are affected by the environment (epigenomic studies) have provided new perspectives on post-traumatic stress disorder and its intergenerational transmission. This article describes the lessons learned from community engagement (CE) in this type of research in a Rwandan genocide-exposed population. A strong trauma-related response was observed within the community while explaining the project. CE also revealed the participants' privacy concerns related to leakage of genetic/(epi)genomic data that could also affect their close relatives. Adopting a culture of CE in the process of research implementation enables the prioritization of community needs and interests. CE has furthermore stimulated the development of preventive interventions for married couples to protect their offspring and thus truly break the cycle of inherited vulnerability.

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.

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.

Leukocyte methylomic imprints of exposure to the genocide against the Tutsi in Rwanda: a pilot epigenome-wide analysis.

Aim & methods: We conducted a pilot epigenome-wide association study of women from Tutsi ethnicity exposed to the genocide while pregnant and their resulting offspring, and a comparison group of women who were pregnant at the time of the genocide but living outside of Rwanda.Results: Fifty-nine leukocyte-derived DNA samples survived quality control: 33 mothers (20 exposed, 13 unexposed) and 26 offspring (16 exposed, 10 unexposed). Twenty-four significant differentially methylated regions (DMRs) were identified in mothers and 16 in children. Conclusions:In utero genocide exposure was associated with CpGs in three of the 24 DMRs: BCOR, PRDM8 and VWDE, with higher DNA methylation in exposed versus unexposed offspring. Of note, BCOR and VWDE show significant correlation between brain and blood DNA methylation within individuals, suggesting these peripherally derived signals of genocide exposure may have relevance to the brain.

Lay abstract The 1994 Rwandan genocide against ethnic Tutsi has been associated with adverse mental health outcomes in survivors decades later, but the molecular mechanisms that contribute to this association remain poorly characterized. Epigenetic mechanisms such as DNA methylation regulate gene function and change in response to life experiences. We identified differentially methylated regions (DMRs) in genocide-exposed versus unexposed mothers and children. In utero genocide exposure was linked with methylation differences in three maternal DMRs, with higher methylation in exposed offspring. Two of three DMRs show correlation between brain and blood methylation within individuals, suggesting that peripherally derived signals of genocide exposure may be relevant to the brain.

The impact of psychopathology, social adversity and stress-relevant DNA methylation on prospective risk for post-traumatic stress: A machine learning approach.

BACKGROUND: A range of factors have been identified that contribute to greater incidence, severity, and prolonged course of post-traumatic stress disorder (PTSD), including: comorbid and/or prior psychopathology; social adversity such as low socioeconomic position, perceived discrimination, and isolation; and biological factors such as genomic variation at glucocorticoid receptor regulatory network (GRRN) genes. This complex etiology and clinical course make identification of people at higher risk of PTSD challenging. Here we leverage machine learning (ML) approaches to identify a core set of factors that may together predispose persons to PTSD. METHODS: We used multiple ML approaches to assess the relationship among DNA methylation (DNAm) at GRRN genes, prior psychopathology, social adversity, and prospective risk for PTS severity (PTSS). RESULTS: ML models predicted prospective risk of PTSS with high accuracy. The Gradient Boost approach was the top-performing model with mean absolute error of 0.135, mean square error of 0.047, root mean square error of 0.217, and R2 of 95.29%. Prior PTSS ranked highest in predicting the prospective risk of PTSS, accounting for >88% of the prediction. The top ranked GRRN CpG site was cg05616442, in AKT1, and the top ranked social adversity feature was loneliness. CONCLUSION: Multiple factors including prior PTSS, social adversity, and DNAm play a role in predicting prospective risk of PTSS. ML models identified factors accounting for increased PTSS risk with high accuracy, which may help to target risk factors that reduce the likelihood or course of PTSD, potentially pointing to approaches that can lead to early intervention. LIMITATION: One of the limitations of this study is small sample size.

Evaluating the impact of trauma and PTSD on epigenetic prediction of lifespan and neural integrity.

Post-traumatic stress disorder (PTSD) is a debilitating disorder that develops in some people following trauma exposure. Trauma and PTSD have been associated with accelerated cellular aging. This study evaluated the effect of trauma and PTSD on accelerated GrimAge, an epigenetic predictor of lifespan, in traumatized civilians. This study included 218 individuals with current PTSD, 427 trauma-exposed controls without any history of PTSD and 209 subjects with lifetime PTSD history who are not categorized as current PTSD cases. The Traumatic Events Inventory (TEI) and Clinician-Administered PTSD Scale (CAPS) were used to measure lifetime trauma burden and PTSD, respectively. DNA from whole blood was interrogated using the MethylationEPIC or HumanMethylation450 BeadChips. GrimAge estimates were calculated using the methylation age calculator. Cortical thickness of 69 female subjects was assessed by using T1-weighted structural MRI images. Associations between trauma exposure, PTSD, cortical thickness, and GrimAge acceleration were tested with multiple regression models. Lifetime trauma burden (p = 0.03), current PTSD (p = 0.02) and lifetime PTSD (p = 0.005) were associated with GrimAge acceleration, indicative of a shorter predicted lifespan. The association with lifetime PTSD was replicated in an independent cohort (p = 0.04). In the MRI sub sample, GrimAge acceleration also associated with cortical atrophy in the right lateral orbitofrontal cortex (padj = 0.03) and right posterior cingulate (padj = 0.04), brain areas associated with emotion-regulation and threat-regulation. Our findings suggest that lifetime trauma and PTSD may contribute to a higher epigenetic-based mortality risk. We also demonstrate a relationship between cortical atrophy in PTSD-relevant brain regions and shorter predicted lifespan.

RANDOMIZE: A Web Server for Data Randomization.

The microarray-based Illumina Infinium MethylationEpic BeadChip (Epic 850k) has become a useful and standard tool for epigenome wide deoxyribonucleic acid (DNA) methylation profiling. Data from this technology may suffer from batch effects due to improper handling of the samples during the plating process. Batch effects are a significant issue and can give rise to spurious and inaccurate results and reduction in power to detect real biological differences. Careful study design, such as randomizing the samples to uniformly distribute the samples across the factors responsible for batch effects, is crucial to address batch effects and other technical artifacts. Randomization helps to reduce the likelihood of bias and impact of difference among groups. This process of randomizing the samples can be a tedious, error-prone, and time-consuming task without a user-friendly and efficient tool. We present RANDOMIZE, a web-based application designed to perform randomization of relevant metadata to evenly distribute samples across the factors typically responsible for batch effects in DNA methylation microarrays, such as rows, chips and plates. We demonstrate that the tool is efficient, fast and easy to use. The tool is freely available online at https://coph-usf.shinyapps.io/RANDOMIZE/ and can be accessed using any web browser. Sample data and tutorial is also available with the tool.