Genome-wide association analyses of post-traumatic stress disorder and its symptom subdomains in the Million Veteran Program.

Manhattan plot for the MVP case-control GWAS (top) and for the MVP/PGC GWAS meta-analysis in EUR samples (bottom).GWAS was performed using logistic regression, co-varying for age, sex, and the first 10 principal components of ancestry. Meta-analysis was conducted with METAL25 using the inverse variance weighting method. Bonferroni correction was used to correct for multiple comparisons; associations with P ≤ 5 × 10−8 (indicated by the horizontal red bar) were considered to be genome-wide significant.

Symptom and polygenic risk scores in veterans of African and European ancestryTop shows density plot of PCL-total scores in veterans of AFR (salmon color) and EUR (teal color) ancestry. Bottom shows density plot of PRS scores (at P-value threshold 0.001) for MVP PCL AFR (salmon color) and MVP PCL EUR (teal color) into PGC PTSD EUR.

Gene Ontology (GO) term and GTEx tissue enrichmenta, Quantile-quantile plots between Gene Ontology (GO) term enrichment (one-sided test for positive relationship between tissue and genetic association) in original PCL-Total and conditioned PCL-Total (blue, autism spectrum disorder; purple, major depression; dark green, anorexia nervosa; light green, anxiety; pink, schizophrenia; light blue, bipolar disorder; orange, attention deficit hyperactivity disorder; red, all eight disorders simultaneously). b, Quantile-quantile relationship between GTEx tissue enrichment (one-sided test for positive relationship between tissue and genetic association) in original PCL-total and conditioned PCL-Total. To avoid over-plotting, enrichment P-values were divided into quantiles. Red diagonal lines indicate a one-to-one relationship between original and conditioned PCL-Total gene set and tissue enrichments. Two-sided tests were used to compare enrichment results.

Genome-wide significant (P < 5 × 10−8) findings, by European (circles) and African (diamonds) ancestry, for PTSD EHR case-control, PCL-Total score and each of the PTSD subcomponents (avoidance, hyperarousal, and re-experiencing).There were no genome-wide significant results for the African case-control and re-experiencing traits. LD independent SNPs for each phenotype and the nearest gene are labeled. The donut chart summarizes the number of hits for each phenotype in the two populations. The genes labelled are significant following regression test for a two-sided P-value and applied Bonferroni-threshold for multiple testing (0.05/k SNPs = P < 5 × 10−8).

Phenotypic (above diagonal) and genetic (below diagonal) correlations between case-control, PCL-Total and subscale scores.Shown are correlation point estimates, 95% CIs, and n (sample size). SNP-heritability (h2) is shown in the left column. For phenotypic correlations, those for case-control are point-biserial correlations, and all others are Pearson correlations.

LDSC genetic correlation (rg) analyses in EUR showing traits from UK Biobank and PGC psychiatric disorders.a, Comparison of rg between PTSD case-control definition and PCL-Total. The grey diagonal indicates perfect linearity between PTSD case-control and PCL-Total genetic correlates. The top ten genetic correlates of PCL-Total are labeled, with purple data points indicating positive rg and green data points indicating negative rg. Data points in a circled in red indicate significant difference in rg magnitude between PTSD and PCL-Total. b, Plot showing a wide range of phenotypes and their rg; the vertical dashed line indicates rg equal to zero.

PCL-Total SNP-heritability (SNP-h2; bars and standard error; right y-axis) and genetic correlation (rg; data points and standard error; left y-axis) relative to original PCL-Total (n = 186,689) after conditioning with each mental health phenotype on the x-axis.

Genomic SEM model with confirmatory factor analysis indicating two correlated factors, the first consisting of PTSD symptoms and the second consisting of anxiety, major depressive disorder, and neuroticism.Comparative fit index (CFI) = 0.999 (typically interpreted as 0.90–0.95 indicating marginal fit).

Association of genetically-regulated transcriptomic changes with PCL-Total (n = 186,689) in the four brain regions identified by the LDSC partitioned heritability analyses.Betas and 95% confidence intervals are reported for each association.

Manhattan plot of MVP AFR case-control GWASHorizontal red line indicates P < 5 × 10−8. P-values are uncorrected. Results are based on logistic regression.

Polygenic risk scores in MVP and PGC-PTSDPolygenic risk score (PRS) from MVP EUR case-control (left) and EUR PCL-total (right) applied to PGC-PTSD13 case-control phenotype with varying P-value thresholds (PT) on the x-axis and explained variance (R2) on the y-axis. The approximate estimate of the explained variance was calculated using a multivariate regression model. P values reported are two sided, and Bonferroni correction accounting for the number of P-value thresholds tested is P = 2.38 × 10−4.