International meta-analysis of PTSD genome-wide association studies identifies sex- and ancestry-specific genetic risk loci.

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

Assessment of Brain Age in Posttraumatic Stress Disorder: Findings from the ENIGMA PTSD and Brain Age Working Groups

Background Posttraumatic stress disorder (PTSD) is associated with markers of accelerated aging. Estimates of brain age, compared to chronological age, may clarify the effects of PTSD on the brain and may inform treatment approaches targeting the neurobiology of aging in the context of PTSD. Method Adult subjects (N = 2229; 56.2% male) aged 18–69 years (mean = 35.6, SD = 11.0) from 21 ENIGMA-PGC PTSD sites underwent T1-weighted brain structural magnetic resonance imaging, and PTSD assessment (PTSD+, n = 884). Previously trained voxel-wise (brainageR) and region-of-interest (BARACUS and PHOTON) machine learning pipelines were compared in a subset of control subjects (n = 386). Linear mixed effects models were conducted in the full sample (those with and without PTSD) to examine the effect of PTSD on brain predicted age difference (brain PAD; brain age − chronological age) controlling for chronological age, sex, and scan site. Results BrainageR most accurately predicted brain age in a subset (n = 386) of controls (brainageR: ICC = 0.71, R = 0.72, MAE = 5.68; PHOTON: ICC = 0.61, R = 0.62, MAE = 6.37; BARACUS: ICC = 0.47, R = 0.64, MAE = 8.80). Using brainageR, a three-way interaction revealed that young males with PTSD exhibited higher brain PAD relative to male controls in young and old age groups; old males with PTSD exhibited lower brain PAD compared to male controls of all ages. Discussion Differential impact of PTSD on brain PAD in younger versus older males may indicate a critical window when PTSD impacts brain aging, followed by age-related brain changes that are consonant with individuals without PTSD. Future longitudinal research is warranted to understand how PTSD impacts brain aging across the lifespan

Smaller total and subregional cerebellar volumes in posttraumatic stress disorder:a mega-analysis by the ENIGMA-PGC PTSD workgroup

Although the cerebellum contributes to higher-order cognitive and emotional functions relevant to posttraumatic stress disorder (PTSD), prior research on cerebellar volume in PTSD is scant, particularly when considering subregions that differentially map on to motor, cognitive, and affective functions. In a sample of 4215 adults (PTSD n = 1642; Control n = 2573) across 40 sites from the ENIGMA-PGC PTSD working group, we employed a new state-of-the-art deep-learning based approach for automatic cerebellar parcellation to obtain volumetric estimates for the total cerebellum and 28 subregions. Linear mixed effects models controlling for age, gender, intracranial volume, and site were used to compare cerebellum volumes in PTSD compared to healthy controls (88% trauma-exposed). PTSD was associated with significant grey and white matter reductions of the cerebellum. Compared to controls, people with PTSD demonstrated smaller total cerebellum volume, as well as reduced volume in subregions primarily within the posterior lobe (lobule VIIB, crus II), vermis (VI, VIII), flocculonodular lobe (lobule X), and corpus medullare (all p -FDR &lt; 0.05). Effects of PTSD on volume were consistent, and generally more robust, when examining symptom severity rather than diagnostic status. These findings implicate regionally specific cerebellar volumetric differences in the pathophysiology of PTSD. The cerebellum appears to play an important role in higher-order cognitive and emotional processes, far beyond its historical association with vestibulomotor function. Further examination of the cerebellum in trauma-related psychopathology will help to clarify how cerebellar structure and function may disrupt cognitive and affective processes at the center of translational models for PTSD.</p

Smaller total and subregional cerebellar volumes in posttraumatic stress disorder:a mega-analysis by the ENIGMA-PGC PTSD workgroup

Although the cerebellum contributes to higher-order cognitive and emotional functions relevant to posttraumatic stress disorder (PTSD), prior research on cerebellar volume in PTSD is scant, particularly when considering subregions that differentially map on to motor, cognitive, and affective functions. In a sample of 4215 adults (PTSD n = 1642; Control n = 2573) across 40 sites from the ENIGMA-PGC PTSD working group, we employed a new state-of-the-art deep-learning based approach for automatic cerebellar parcellation to obtain volumetric estimates for the total cerebellum and 28 subregions. Linear mixed effects models controlling for age, gender, intracranial volume, and site were used to compare cerebellum volumes in PTSD compared to healthy controls (88% trauma-exposed). PTSD was associated with significant grey and white matter reductions of the cerebellum. Compared to controls, people with PTSD demonstrated smaller total cerebellum volume, as well as reduced volume in subregions primarily within the posterior lobe (lobule VIIB, crus II), vermis (VI, VIII), flocculonodular lobe (lobule X), and corpus medullare (all p -FDR &lt; 0.05). Effects of PTSD on volume were consistent, and generally more robust, when examining symptom severity rather than diagnostic status. These findings implicate regionally specific cerebellar volumetric differences in the pathophysiology of PTSD. The cerebellum appears to play an important role in higher-order cognitive and emotional processes, far beyond its historical association with vestibulomotor function. Further examination of the cerebellum in trauma-related psychopathology will help to clarify how cerebellar structure and function may disrupt cognitive and affective processes at the center of translational models for PTSD.</p

A Comparison of Methods to Harmonize Cortical Thickness Measurements Across Scanners and Sites

Results of neuroimaging datasets aggregated from multiple sites may be biased by site-specific profiles in participants’ demographic and clinical characteristics, as well as MRI acquisition protocols and scanning platforms. We compared the impact of four different harmonization methods on results obtained from analyses of cortical thickness data: (1) linear mixed-effects model (LME) that models site-specific random intercepts (LME INT), (2) LME that models both site-specific random intercepts and age-related random slopes (LME INT+SLP), (3) ComBat, and (4) ComBat with a generalized additive model (ComBat-GAM). Our test case for comparing harmonization methods was cortical thickness data aggregated from 29 sites, which included 1,340 cases with posttraumatic stress disorder (PTSD) (6.2–81.8 years old) and 2,057 trauma-exposed controls without PTSD (6.3–85.2 years old). We found that, compared to the other data harmonization methods, data processed with ComBat-GAM was more sensitive to the detection of significant case-control differences (Χ 2(3) = 63.704, p < 0.001) as well as case-control differences in age-related cortical thinning (Χ 2(3) = 12.082, p = 0.007). Both ComBat and ComBat-GAM outperformed LME methods in detecting sex differences (Χ 2(3) = 9.114, p = 0.028) in regional cortical thickness. ComBat-GAM also led to stronger estimates of age-related declines in cortical thickness (corrected p-values < 0.001), stronger estimates of case-related cortical thickness reduction (corrected p-values < 0.001), weaker estimates of age-related declines in cortical thickness in cases than controls (corrected p-values < 0.001), stronger estimates of cortical thickness reduction in females than males (corrected p-values < 0.001), and stronger estimates of cortical thickness reduction in females relative to males in cases than controls (corrected p-values < 0.001). Our results support the use of ComBat-GAM to minimize confounds and increase statistical power when harmonizing data with non-linear effects, and the use of either ComBat or ComBat-GAM for harmonizing data with linear effects

Remodeling of the Cortical Structural Connectome in Posttraumatic Stress Disorder:Results from the ENIGMA-PGC PTSD Consortium

BACKGROUND: Posttraumatic stress disorder (PTSD) is accompanied by disrupted cortical neuroanatomy. We investigated alteration in covariance of structural networks associated with PTSD in regions that demonstrate the case-control differences in cortical thickness (CT) and surface area (SA). METHODS: Neuroimaging and clinical data were aggregated from 29 research sites in >1,300 PTSD cases and >2,000 trauma-exposed controls (age 6.2-85.2 years) by the ENIGMA-PGC PTSD working group. Cortical regions in the network were rank-ordered by effect size of PTSD-related cortical differences in CT and SA. The top-n (n = 2 to 148) regions with the largest effect size for PTSD > non-PTSD formed hypertrophic networks, the largest effect size for PTSD < non-PTSD formed atrophic networks, and the smallest effect size of between-group differences formed stable networks. The mean structural covariance (SC) of a given n-region network was the average of all positive pairwise correlations and was compared to the mean SC of 5,000 randomly generated n-region networks. RESULTS: Patients with PTSD, relative to non-PTSD controls, exhibited lower mean SC in CT-based and SA-based atrophic networks. Comorbid depression, sex and age modulated covariance differences of PTSD-related structural networks. CONCLUSIONS: Covariance of structural networks based on CT and cortical SA are affected by PTSD and further modulated by comorbid depression, sex, and age. The structural covariance networks that are perturbed in PTSD comport with converging evidence from resting state functional connectivity networks and networks impacted by inflammatory processes, and stress hormones in PTSD

Potential causal association between gut microbiome and posttraumatic stress disorder

Background: The causal effects of gut microbiome and the development of posttraumatic stress disorder (PTSD) are still unknown. This study aimed to clarify their potential causal association using mendelian randomization (MR). Methods: The summary-level statistics for gut microbiome were retrieved from a genome-wide association study (GWAS) of the MiBioGen consortium. As to PTSD, the Freeze 2 datasets were originated from the Psychiatric Genomics Consortium Posttraumatic Stress Disorder Working Group (PGC-PTSD), and the replicated datasets were obtained from FinnGen consortium. Single nucleotide polymorphisms meeting MR assumptions were selected as instrumental variables. The inverse variance weighting (IVW) method was employed as the main approach, supplemented by sensitivity analyses to evaluate potential pleiotropy and heterogeneity and ensure the robustness of the MR results. We also performed reverse MR analyses to explore PTSD’s causal effects on the relative abundances of specific features of the gut microbiome. Results: In Freeze 2 datasets from PGC-PTSD, eight bacterial traits revealed a potential causal association between gut microbiome and PTSD (IVW, all P < 0.05). In addition, Genus.Dorea and genus.Sellimonas were replicated in FinnGen datasets, in which eight bacterial traits revealed a potential causal association between gut microbiome and the occurrence of PTSD. The heterogeneity and pleiotropy analyses further supported the robustness of the IVW findings, providing additional evidence for their reliability. Conclusion: Our study provides the potential causal impact of gut microbiomes on the development of PTSD, shedding new light on the understanding of the dysfunctional gut-brain axis in this disorder. Our findings present novel evidence and call for investigations to confirm the association between their links, as well as to illuminate the underlying mechanisms

Enhancing Discovery of Genetic Variants for Posttraumatic Stress Disorder Through Integration of Quantitative Phenotypes and Trauma Exposure Information

Funding Information: This work was supported by the National Institute of Mental Health / U.S. Army Medical Research and Development Command (Grant No. R01MH106595 [to CMN, IL, MBS, KJRe, and KCK], National Institutes of Health (Grant No. 5U01MH109539 to the Psychiatric Genomics Consortium ), and Brain & Behavior Research Foundation (Young Investigator Grant [to KWC]). Genotyping of samples was provided in part through the Stanley Center for Psychiatric Genetics at the Broad Institute supported by Cohen Veterans Bioscience . Statistical analyses were carried out on the LISA/Genetic Cluster Computer ( https://userinfo.surfsara.nl/systems/lisa ) hosted by SURFsara. This research has been conducted using the UK Biobank resource (Application No. 41209). This work would have not been possible without the financial support provided by Cohen Veterans Bioscience, the Stanley Center for Psychiatric Genetics at the Broad Institute, and One Mind. Funding Information: MBS has in the past 3 years received consulting income from Actelion, Acadia Pharmaceuticals, Aptinyx, Bionomics, BioXcel Therapeutics, Clexio, EmpowerPharm, GW Pharmaceuticals, Janssen, Jazz Pharmaceuticals, and Roche/Genentech and has stock options in Oxeia Biopharmaceuticals and Epivario. In the past 3 years, NPD has held a part-time paid position at Cohen Veterans Bioscience, has been a consultant for Sunovion Pharmaceuticals, and is on the scientific advisory board for Sentio Solutions for unrelated work. In the past 3 years, KJRe has been a consultant for Datastat, Inc., RallyPoint Networks, Inc., Sage Pharmaceuticals, and Takeda. JLM-K has received funding and a speaking fee from COMPASS Pathways. MU has been a consultant for System Analytic. HRK is a member of the Dicerna scientific advisory board and a member of the American Society of Clinical Psychopharmacology Alcohol Clinical Trials Initiative, which during the past 3 years was supported by Alkermes, Amygdala Neurosciences, Arbor Pharmaceuticals, Dicerna, Ethypharm, Indivior, Lundbeck, Mitsubishi, and Otsuka. HRK and JG are named as inventors on Patent Cooperative Treaty patent application number 15/878,640, entitled “Genotype-guided dosing of opioid agonists,” filed January 24, 2018. RP and JG are paid for their editorial work on the journal Complex Psychiatry. OAA is a consultant to HealthLytix. All other authors report no biomedical financial interests or potential conflicts of interest. Funding Information: This work was supported by the National Institute of Mental Health/ U.S. Army Medical Research and Development Command (Grant No. R01MH106595 [to CMN, IL, MBS, KJRe, and KCK], National Institutes of Health (Grant No. 5U01MH109539 to the Psychiatric Genomics Consortium), and Brain & Behavior Research Foundation (Young Investigator Grant [to KWC]). Genotyping of samples was provided in part through the Stanley Center for Psychiatric Genetics at the Broad Institute supported by Cohen Veterans Bioscience. Statistical analyses were carried out on the LISA/Genetic Cluster Computer (https://userinfo.surfsara.nl/systems/lisa) hosted by SURFsara. This research has been conducted using the UK Biobank resource (Application No. 41209). This work would have not been possible without the financial support provided by Cohen Veterans Bioscience, the Stanley Center for Psychiatric Genetics at the Broad Institute, and One Mind. This material has been reviewed by the Walter Reed Army Institute of Research. There is no objection to its presentation and/or publication. The opinions or assertions contained herein are the private views of the authors and are not to be construed as official or as reflecting true views of the U.S. Department of the Army or the Department of Defense. We thank the investigators who comprise the PGC-PTSD working group and especially the more than 206,000 research participants worldwide who shared their life experiences and biological samples with PGC-PTSD investigators. We thank Mark Zervas for his critical input. Full acknowledgments are in Supplement 1. MBS has in the past 3 years received consulting income from Actelion, Acadia Pharmaceuticals, Aptinyx, Bionomics, BioXcel Therapeutics, Clexio, EmpowerPharm, GW Pharmaceuticals, Janssen, Jazz Pharmaceuticals, and Roche/Genentech and has stock options in Oxeia Biopharmaceuticals and Epivario. In the past 3 years, NPD has held a part-time paid position at Cohen Veterans Bioscience, has been a consultant for Sunovion Pharmaceuticals, and is on the scientific advisory board for Sentio Solutions for unrelated work. In the past 3 years, KJRe has been a consultant for Datastat, Inc. RallyPoint Networks, Inc. Sage Pharmaceuticals, and Takeda. JLM-K has received funding and a speaking fee from COMPASS Pathways. MU has been a consultant for System Analytic. HRK is a member of the Dicerna scientific advisory board and a member of the American Society of Clinical Psychopharmacology Alcohol Clinical Trials Initiative, which during the past 3 years was supported by Alkermes, Amygdala Neurosciences, Arbor Pharmaceuticals, Dicerna, Ethypharm, Indivior, Lundbeck, Mitsubishi, and Otsuka. HRK and JG are named as inventors on Patent Cooperative Treaty patent application number 15/878,640, entitled ?Genotype-guided dosing of opioid agonists,? filed January 24, 2018. RP and JG are paid for their editorial work on the journal Complex Psychiatry. OAA is a consultant to HealthLytix. All other authors report no biomedical financial interests or potential conflicts of interest. Publisher Copyright: © 2021 Society of Biological PsychiatryBackground: 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.publishersversionpublishe

Rare copy number variation in posttraumatic stress disorder

Posttraumatic stress disorder (PTSD) is a heritable (h2 = 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