Posttraumatic widespread pain (PTWP) and posttraumatic stress symptoms (PTSS) are frequent comorbid sequelae of trauma that occur at different rates in women and men. We sought to identify microRNA (miRNA) that may contribute to sex-dependent differences in vulnerability to these outcomes. Monte Carlo simulations (x10,000) identified miRNA in which predicted targeting of PTWP or PTSS genes was most enriched. Expression of the leading candidate miRNA to target PTWP/PTSS-related genes, miR-19b, has been shown to be influenced by estrogen and stress exposure. We evaluated whether peritraumatic miR-19b blood expression levels predicted PTWP and PTSS development in women and men experiencing trauma of motor vehicle collision (n = 179) and in women experiencing sexual assault trauma (n = 74). A sex-dependent relationship was observed between miR-19b expression levels and both PTWP (β = -2.41, P = 0.034) and PTSS (β = -3.01, P = 0.008) development 6 months after motor vehicle collision. The relationship between miR-19b and PTSS (but not PTWP) was validated in sexual assault survivors (β = -0.91, P = 0.013). Sex-dependent expression of miR-19b was also observed in blood and nervous tissue from 2 relevant animal models. Furthermore, in support of increasing evidence indicating a role for the circadian rhythm (CR) in PTWP and PTSS pathogenesis, miR-19b targets were enriched in CR gene transcripts. Human cohort and in vitro analyses assessing miR-19b regulation of key CR transcripts, CLOCK and RORA, supported the potential importance of miR-19b to regulating the CR pathway. Together, these results highlight the potential role that sex-dependent expression of miR-19b might play in PTWP and PTSS development after trauma/stress exposure.

Posttraumatic stress disorder (PTSD) is a heritable (h² = 24-71%) psychiatric illness. Copy number variation (CNV) is a form of rare genetic variation that has been implicated in the etiology of psychiatric disorders, but no large-scale investigation of CNV in PTSD has been performed. We present an association study of CNV burden and PTSD symptoms in a sample of 114,383 participants (13,036 cases and 101,347 controls) of European ancestry. CNVs were called using two calling algorithms and intersected to a consensus set. Quality control was performed to remove strong outlier samples. CNVs were examined for association with PTSD within each cohort using linear or logistic regression analysis adjusted for population structure and CNV quality metrics, then inverse variance weighted meta-analyzed across cohorts. We examined the genome-wide total span of CNVs, enrichment of CNVs within specified gene-sets, and CNVs overlapping individual genes and implicated neurodevelopmental regions. The total distance covered by deletions crossing over known neurodevelopmental CNV regions was significant (beta = 0.029, SE = 0.005, P = 6.3 × 10^-8). The genome-wide neurodevelopmental CNV burden identified explains 0.034% of the variation in PTSD symptoms. The 15q11.2 BP1-BP2 microdeletion region was significantly associated with PTSD (beta = 0.0206, SE = 0.0056, P = 0.0002). No individual significant genes interrupted by CNV were identified. 22 gene pathways related to the function of the nervous system and brain were significant in pathway analysis (FDR q < 0.05), but these associations were not significant once NDD regions were removed. A larger sample size, better detection methods, and annotated resources of CNV are needed to explore this relationship further.

Background

Posttraumatic stress disorder (PTSD) and major depressive disorder (MDD) are commonly reported co-occurring mental health consequences of psychological trauma exposure. The disorders have high genetic overlap. Trauma is a complex phenotype but research suggests that trauma sensitivity has a heritable basis. We investigated whether sensitivity to trauma in those with MDD reflects a similar genetic component in those with PTSD.

Methods

Genetic correlations between PTSD and MDD in individuals reporting trauma and MDD in individuals not reporting trauma were estimated, as well as with recurrent MDD and single-episode MDD, using genome-wide association study (GWAS) summary statistics. Genetic correlations were replicated using PTSD data from the Psychiatric Genomics Consortium and the Million Veteran Program. Polygenic risk scores were generated in UK Biobank participants who met the criteria for lifetime MDD (N = 29 471). We investigated whether genetic loading for PTSD was associated with reporting trauma in these individuals.

Results

Genetic loading for PTSD was significantly associated with reporting trauma in individuals with MDD [OR 1.04 (95% CI 1.01-1.07), Empirical-p = 0.02]. PTSD was significantly more genetically correlated with recurrent MDD than with MDD in individuals not reporting trauma (rg differences = ~0.2, p < 0.008). Participants who had experienced recurrent MDD reported significantly higher rates of trauma than participants who had experienced single-episode MDD (χ2 > 166, p < 0.001).

Conclusions

Our findings point towards the existence of genetic variants associated with trauma sensitivity that might be shared between PTSD and MDD, although replication with better powered GWAS is needed. Our findings corroborate previous research highlighting trauma exposure as a key risk factor for recurrent MDD.

The risk of posttraumatic stress disorder (PTSD) following trauma is heritable, but robust common variants have yet to be identified. In a multi-ethnic cohort including over 30,000 PTSD cases and 170,000 controls we conduct a genome-wide association study of PTSD. We demonstrate SNP-based heritability estimates of 5-20%, varying by sex. Three genome-wide significant loci are identified, 2 in European and 1 in African-ancestry analyses. Analyses stratified by sex implicate 3 additional loci in men. Along with other novel genes and non-coding RNAs, a Parkinson's disease gene involved in dopamine regulation, PARK2, is associated with PTSD. Finally, we demonstrate that polygenic risk for PTSD is significantly predictive of re-experiencing symptoms in the Million Veteran Program dataset, although specific loci did not replicate. These results demonstrate the role of genetic variation in the biology of risk for PTSD and highlight the necessity of conducting sex-stratified analyses and expanding GWAS beyond European ancestry populations.

Background

Posttraumatic stress disorder (PTSD) is heritable and a potential consequence of exposure to traumatic stress. Evidence suggests that a quantitative approach to PTSD phenotype measurement and incorporation of lifetime trauma exposure (LTE) information could enhance the discovery power of PTSD genome-wide association studies (GWASs).

Methods

A GWAS on PTSD symptoms was performed in 51 cohorts followed by a fixed-effects meta-analysis (N = 182,199 European ancestry participants). A GWAS of LTE burden was performed in the UK Biobank cohort (N = 132,988). Genetic correlations were evaluated with linkage disequilibrium score regression. Multivariate analysis was performed using Multi-Trait Analysis of GWAS. Functional mapping and annotation of leading loci was performed with FUMA. Replication was evaluated using the Million Veteran Program GWAS of PTSD total symptoms.

Results

GWASs of PTSD symptoms and LTE burden identified 5 and 6 independent genome-wide significant loci, respectively. There was a 72% genetic correlation between PTSD and LTE. PTSD and LTE showed largely similar patterns of genetic correlation with other traits, albeit with some distinctions. Adjusting PTSD for LTE reduced PTSD heritability by 31%. Multivariate analysis of PTSD and LTE increased the effective sample size of the PTSD GWAS by 20% and identified 4 additional loci. Four of these 9 PTSD loci were independently replicated in the Million Veteran Program.

Conclusions

Through using a quantitative trait measure of PTSD, we identified novel risk loci not previously identified using prior case-control analyses. PTSD and LTE have a high genetic overlap that can be leveraged to increase discovery power through multivariate methods.