ORIGINAL RESEARCH article
Volume 15 - 2021 | https://doi.org/10.3389/fnins.2021.678503
Examining Individual and Synergistic Contributions of PTSD and Genetics to Blood Pressure: A Trans-Ethnic Meta-Analysis
- 1Department of Psychology, University of California, Los Angeles, Los Angeles, CA, United States
- 2Department of Psychiatry, University of California, San Diego, San Diego, CA, United States
- 3Veterans Affairs San Diego Healthcare System, VA Center of Excellence for Stress and Mental Health (CESAMH), San Diego, CA, United States
- 4Department of Psychiatry and Behavioral Sciences, Emory University, Atlanta, GA, United States
- 5Yerkes National Primate Research Center, Atlanta, GA, United States
- 6Department of Psychological Sciences, Case Western Reserve University, Cleveland, OH, United States
- 7National Center on Birth Defects and Developmental Disabilities, Centers for Disease Control and Prevention, Atlanta, GA, United States
- 8NORMENT, Institute of Clinical Medicine, University of Oslo, Oslo, Norway
- 9Division of Mental Health and Addiction, Oslo University Hospital, Oslo, Norway
- 10Duke Molecular Physiology Institute, Duke University, Durham, NC, United States
- 11Durham VA Health Care System, Durham, NC, United States
- 12Department of Psychiatry and Behavioral Sciences, Duke University School of Medicine, Durham, NC, United States
- 13VA Mid-Atlantic Mental Illness Research, Education, and Clinical Center (MIRECC), Genetics Research Laboratory, Durham, NC, United States
- 14Atlanta VA Health Care System, Decatur, GA, United States
- 15Genetic and Developmental Psychiatry Centre, Institute of Psychiatry, Psychology, and Neuroscience, King’s College London, London, United Kingdom
- 16NIHR BRC at the Maudsley, King’s College London, London, United Kingdom
- 17Department of Radiology, University of California, San Diego, San Diego, CA, United States
- 18Department of Neurosciences, University of California, San Diego, San Diego, CA, United States
- 19Department of Psychiatry, Harvard Medical School, Boston, MA, United States
- 20McLean Hospital, Belmont, MA, United States
- 21Department of Psychiatry, Faculty of Medicine and Health Sciences, Stellenbosch University, Cape Town, South Africa
- 22South African Medical Research Council/Stellenbosch University Genomics of Brain Disorders Research Unit, Faculty of Medicine & Health Sciences, Stellenbosch University, Cape Town, South Africa
- 23Department of Psychiatry and Behavioral Neurosciences, School of Medicine, Wayne State University, Detroit, MI, United States
- 24School of Biomedical Sciences, Queensland University of Technology, Kelvin Grove, QLD, Australia
- 25National Center for PTSD, Behavioral Science Division, VA Boston Healthcare System, Boston, MA, United States
- 26Department of Psychiatry, Boston University School of Medicine, Boston, MA, United States
- 27Department of Biostatistics, Boston University School of Public Health, Boston, MA, United States
- 28Biomedical Genetics, Boston University School of Medicine, Boston, MA, United States
- 29Department of Gynecology and Obstetrics, Emory University, Atlanta, GA, United States
- 30Department of Psychological & Brain Sciences, Boston University, Boston, MA, United States
- 31Department of Psychiatry, Case Western Reserve University, Cleveland, OH, United States
- 32GRECC/TRACTS, VA Boston Healthcare System, Boston, MA, United States
- 33Center for Genomics and Personalised Health, Queensland University of Technology, Kelvin Grove, QLD, Australia
- 34Gallipoli Medical Research Foundation, Greenslopes Private Hospital, Brisbane, QLD, Australia
- 35Department of Psychology, University of Washington, Seattle, WA, United States
- 36School of Psychology and Counseling, Queensland University of Technology, Kelvin Grove, QLD, Australia
- 37Department of Psychiatry and Behavioral Sciences, University of Washington, Seattle, WA, United States
Growing research suggests that posttraumatic stress disorder (PTSD) may be a risk factor for poor cardiovascular health, and yet our understanding of who might be at greatest risk of adverse cardiovascular outcomes after trauma is limited. In this study, we conducted the first examination of the individual and synergistic contributions of PTSD symptoms and blood pressure genetics to continuous blood pressure levels. We harnessed the power of the Psychiatric Genomics Consortium-PTSD Physical Health Working Group and investigated these associations across 11 studies of 72,224 trauma-exposed individuals of European (n = 70,870) and African (n = 1,354) ancestry. Genetic contributions to blood pressure were modeled via polygenic scores (PGS) for systolic blood pressure (SBP) and diastolic blood pressure (DBP) that were derived from a prior trans-ethnic blood pressure genome-wide association study (GWAS). Results of trans-ethnic meta-analyses revealed significant main effects of the PGS on blood pressure levels [SBP: β = 2.83, standard error (SE) = 0.06, p < 1E-20; DBP: β = 1.32, SE = 0.04, p < 1E-20]. Significant main effects of PTSD symptoms were also detected for SBP and DBP in trans-ethnic meta-analyses, though there was significant heterogeneity in these results. When including data from the largest contributing study – United Kingdom Biobank – PTSD symptoms were negatively associated with SBP levels (β = −1.46, SE = 0.44, p = 9.8E-4) and positively associated with DBP levels (β = 0.70, SE = 0.26, p = 8.1E-3). However, when excluding the United Kingdom Biobank cohort in trans-ethnic meta-analyses, there was a nominally significant positive association between PTSD symptoms and SBP levels (β = 2.81, SE = 1.13, p = 0.01); no significant association was observed for DBP (β = 0.43, SE = 0.78, p = 0.58). Blood pressure PGS did not significantly moderate the associations between PTSD symptoms and blood pressure levels in meta-analyses. Additional research is needed to better understand the extent to which PTSD is associated with high blood pressure and how genetic as well as contextual factors may play a role in influencing cardiovascular risk.
Cardiovascular disease (CVD) is the leading cause of death and disability worldwide (Virani et al., 2020), and longitudinal research has found that posttraumatic stress disorder (PTSD) – the quintessential trauma-related mental disorder – precedes and predicts the onset of incident CVD (Kubzansky et al., 2007, 2009; Vaccarino et al., 2013; Sumner et al., 2015). Elevated blood pressure may be one mechanism by which PTSD contributes to heightened CVD risk. Meta-analytic evidence suggests that individuals with (versus without) PTSD have significantly greater resting systolic blood pressure (SBP) and diastolic blood pressure (DBP), as well as greater DBP response to trauma-related cues (Buckley and Kaloupek, 2001; Pole, 2007). In addition, PTSD has been associated with elevated blood pressure and hypertension in a variety of trauma-exposed populations, including military veteran and community-based samples (Schnurr et al., 2000; Kibler et al., 2009; Sumner et al., 2016; Burg et al., 2017; Edmondson et al., 2018; Howard et al., 2018).
Although growing evidence suggests that PTSD may be linked to elevated blood pressure – a major risk factor for CVD – there is a lack of understanding of who may be at greatest risk. Harnessing information about genetic liability may help to shed light on who is most vulnerable to developing adverse physical health outcomes after trauma. For example, in a predominantly male sample of younger veterans, Wolf et al. (2017) found that obesity polygenic risk, reflecting additive risk for obesity based on thousands of genotypes across the genome, moderated the relation between PTSD and metabolic syndrome (an indicator of cardiometabolic risk). Specifically, lifetime PTSD symptom severity was more strongly associated with metabolic syndrome severity among individuals with greater polygenic risk for obesity (Wolf et al., 2017). Thus, examining whether genetic liability moderates associations of PTSD with physical health indicators, such as blood pressure (i.e., Gene × Environment interactions) may further our understanding of chronic disease risk and inform more targeted allocation of screening and intervention efforts.
Recent work has documented the genetic underpinnings of blood pressure in both European ancestry and trans-ethnic samples (Franceschini et al., 2013; Warren et al., 2017; Giri et al., 2019). These large-scale genome-wide association studies (GWAS) have identified over 250 genetic loci that are significantly associated with blood pressure outcomes, such as SBP and DBP. In the most recent trans-ethnic meta-analysis in over 750,000 individuals, Giri et al. (2019) replicated 216 previously identified SBP loci and 76 DBP loci, and identified 124 novel loci for SBP and 4 for DBP in analyses of common variants; a smaller number of loci were identified based on analyses of rare variants. Notably, in line with prior research (Franceschini et al., 2013), the direction of effects of these single nucleotide polymorphisms (SNPs) on blood pressure were consistent across race and ethnicity. Furthermore, polygenic scores (PGS) based on the GWAS results were significantly and positively associated with a clinical outcome – hypertension; associations were similar across racial and ethnic groups (Giri et al., 2019).
This study investigated whether blood pressure PGS moderated the associations of PTSD with SBP and DBP. Through a collaborative effort of the Psychiatric Genomics Consortium-PTSD (PGC-PTSD) Physical Health Working Group, we examined main effects of PTSD and blood pressure PGS on continuous blood pressure levels, as well as their interaction, in 72,224 trauma-exposed individuals of European and African ancestry. To estimate associations across studies and examine the role of ancestry, we conducted trans-ethnic, as well as single ancestry, meta-analyses. We hypothesized that both PTSD and blood pressure PGS would be associated with higher SBP and DBP levels and that the association between PTSD and SBP/DBP would be greatest for those with higher blood pressure PGS.
Materials and Methods
Eleven studies participating in the PGC-PTSD consortium and PGC-PTSD Physical Health Working Group contributed data to the current investigation. These studies included the CHOICE study (N = 78), D-cycloserine study (DCS; N = 80), Grady Trauma Project (GTP; N = 737), Gallipoli Medical Research Foundation-Queensland University of Technology (GMRF-QUT; N = 275), Hostility, PTSD, and Physical Health Risk Factors (HOST; N = 121), Marine Resiliency Study (MRS; N = 1,933), Optimizing Treatment for PTSD (OPT; N = 75), Shared Roots Study (SHRS; N = 292), Translational Research Center for TBI and Stress Disorders (TRACTS; N = 312), United Kingdom Biobank (UKBB; N = 67,270), and Vietnam Era Twin Study of Aging (VETSA; N = 1,051). Inclusion criteria were as follows: (1) study administered a measure of PTSD symptoms; (2) blood pressure was measured directly; and (3) genetic data were available. Blood pressure measurement was conducted at the same time or after PTSD assessment. If a study included more than one blood pressure measurement (i.e., measures of blood pressure at more than one time point, in contrast to multiple readings of blood pressure collected at a single time point and averaged to obtain a single blood pressure assessment), investigators selected the earliest blood pressure measure administered to participants in the study that occurred at or after the assessment of PTSD symptoms. In cohorts that received treatment (CHOICE, DCS, and OPT), PTSD and blood pressure measurements used were obtained prior to randomization. Cohort details, including information about measures of PTSD symptoms and blood pressure, are provided in the Supplementary Material.
Each study administered a measure of PTSD symptoms (see Supplementary Material for details). To facilitate comparisons across studies given differences in the measures used, total PTSD symptom severity scores were rescaled to a 0–1 range based on the theoretical ranges of each inventory [(observed score – minimum possible score)/maximum possible score]. This 0–1 scale means that each individual’s score is akin to a percentage of maximum possible symptom severity and thus is directly comparable across cohorts. Examining this rescaled continuous measure of PTSD symptoms, rather than a categorical PTSD indicator, allowed us to model individual differences in severity across studies on a common scale and maximize power (e.g., MacCallum et al., 2002).
We examined SBP and DBP (see Supplementary Material for details) as outcomes in separate analyses. As recommended by Tobin et al. (2005), antihypertensive medication use was accounted for by adding a constant to SBP and DBP values. Specifically, if use of antihypertensive medication was reported for a given participant, 15 mmHg was added to the SBP value and 10 mmHg was added to the DBP value, as in prior research (Tobin et al., 2005; Warren et al., 2017). CHOICE, DCS, and UKBB did not include a measure of antihypertensive medication, and thus blood pressure values were not adjusted in these cohorts.
Genotype Quality Control, Imputation, and Ancestry Determination
Quality control for genotypes was conducted following standard procedures and imputed to the 1,000 Genomes Phase 3 reference panel, as previously described in detail (Bycroft et al., 2018; Nievergelt et al., 2019).
Global ancestry was determined using the program SNPweights (Chen et al., 2013), using a reference panel of 10,000 ancestry informative markers measured in 71 reference populations and 6 continental groups, as previously described (Nievergelt et al., 2019). Pre-quality control genotypes were used for these analyses. The current study focused on participants of European and African ancestry. European ancestry participants were defined as those with >90% European ancestry and <5% ancestry in the remaining reference populations. African American participants were defined as those with >10% African ancestry and <5% ancestry in the other non-European ancestry populations (i.e., placing them on an African-European population gradient), and they were analyzed in the African ancestry group. Participants in SHRS were classified broadly as “South African” and included in the African ancestry group; strong outliers [± 4 standard deviations (SDs) of the mean principal component (PC) for any of the first five ancestry PCs] were removed.
Blood Pressure PGS
Blood pressure PGS were calculated using PRSice v2.3.3 (Choi and O’Reilly, 2019) based on effect sizes from the trans-ethnic GWAS of independent samples from the Million Veteran Program (Giri et al., 2019). SBP and DBP GWAS summary statistics were used to calculate SBP and DBP PGS, respectively1. We filtered the summary statistics to remove SNPs with minor allele frequency (MAF) < 1%, SNPs with effects estimated in less than 90% of the total sample size, or strand ambiguous SNPs. Imputed genotypes for PGC-PTSD datasets were converted from probabilities to hard calls, where the genotype called was the one with the highest imputation probability. If the highest imputation probability was < 0.8, the genotype call was instead set to missing. Only markers with < 5% missing rate and MAF > 1% (calculated within each dataset) were retained. For computation of PGS, variants were linkage disequilibrium (LD) clumped (LD information was calculated within each dataset) based on a 250 kb window where only variants with r2 < 0.1 were retained. PGS were calculated at multiple p-value thresholds (pT), varying by increments of 0.01 and ranging from 5E-8 to 1.
Each study provided information on participant age (in years) and biological sex. In addition, ancestry PCs were calculated. Specifically, within each dataset, related participants were identified using identity by state estimation in PLINK 1.9 (Purcell et al., 2007). For every related pair (pi_hat > 0.2), one individual was removed. PCs were calculated in unrelated participants of the same genetically determined global ancestry using the smartPCA algorithm implemented in Eigenstrat 6.1.4 (Price et al., 2006). For PC calculation, SNPs were excluded if they had MAF < 5%, Hardy Weinberg Equilibrium p > 1E-3, or ambiguous strand, or were located in the major histocompatibility complex region (chr. 6, 25–35 MB) or chromosome 8 inversion (chr. 8, 7–13 MB).
To enhance comparability of the blood pressure PGS used across cohorts, we identified the best fitting PGS pT for the largest sample: UKBB. In this sample, the various PGS were tested for associations with their corresponding blood pressure measure in linear models adjusting for the first five ancestry PCs, age, age-squared, and sex. The best-fitting thresholds in UKBB were then used for all analyses of blood pressure PGS in the remaining datasets. To make effect estimates comparable, PGS were centered to the UKBB cohort sample means and SDs.
Analyses were performed stratified by cohort. As in prior research in the PGC-PTSD (Nievergelt et al., 2019), CHOICE, DCS, and OPT – three treatment studies that contributed pre-treatment assessments of PTSD and blood pressure for the current study – were analyzed together as a single cohort. Main effects of (1) blood pressure PGS and (2) PTSD symptoms were each tested using linear regression models with blood pressure levels as the outcome, adjusting for the first five ancestry PCs, age, age-squared, and sex (single-sex studies did not include a sex covariate). As the VETSA study contained twin pairs, data were analyzed using a linear mixed model accounting for non-independence of observations, as implemented in the lme4 R package (Bates et al., 2015). Interactions of PTSD symptoms and blood pressure PGS were tested by fitting a linear model with main effects of PTSD symptoms, blood pressure PGS, and the interaction term all included in the model. Interactions of all covariates with PTSD symptoms and with blood pressure PGS were included as covariates for these models per the recommendation of Keller (2014), along with covariate main effects2.
Effect estimates from each cohort were meta-analyzed using fixed effects inverse variance weighted meta-analysis, implemented in the metafor R package (Viechtbauer, 2010). We conducted trans-ethnic meta-analyses including all cohorts, in addition to separate meta-analyses of the African ancestry and European ancestry participants. To assess results for heterogeneity, Cochran’s Q and the I2 statistic were calculated (Higgins et al., 2003). Cochran’s Q statistic tests the null hypothesis that the included studies are evaluating the same effect; a p-value is generated by comparing the Q statistic with a χ2 distribution with degrees of freedom equal to the number of studies minus 1. However, the Q statistic is influenced by the number of studies included in a meta-analysis (Higgins et al., 2003). As a result, the I2 statistic, which quantifies the percent of total variability across studies that is due to heterogeneity rather than chance, is frequently also calculated when assessing heterogeneity. The I2 statistic ranges from 0 to 100%, with values of 25, 50, and 75% suggesting low, moderate, and high heterogeneity, respectively (Higgins et al., 2003). Diagnostic plots were produced and assessed for fit and influence of potential outliers. As the UKBB sample alone was substantially larger than the combined remaining cohorts, meta-analyses were performed with and without the UKBB sample, and z-tests were used to compare effect size estimate differences between meta-analyses including and excluding the UKBB sample.
We also conducted several secondary analyses to further explore our results. First, sex-specific PGS and PTSD effects were tested in UKBB (i.e., the largest cohort) by conducting sex-stratified analyses and by fitting Sex × PTSD Symptoms and Sex × Blood Pressure PGS interaction terms in models with PTSD symptoms and blood pressure PGS main effects in the full cohort. Again, per Keller (2014), interactions of all covariates with the variables in interactions were included in models testing interactions. Second, we conducted meta-analyses stratified by cohort type – military or community-based – to investigate whether the nature of the sample influenced results. For these analyses, we excluded the HOST study and combined CHOICE, DCS, and OPT cohort, as they enrolled both military and civilian participants (included cohorts are described in Supplementary Table 2). Third, as some studies (UKBB, CHOICE, and DCS) did not assess antihypertensive medication use and thus contributed unadjusted blood pressure values to meta-analyses, we restricted meta-analyses to those studies with adjusted blood pressure values based on the presence of antihypertensive medication use. With these analyses, we were able to examine the patterns of results that emerged when recommended guidelines for accounting for antihypertensive medication use (Tobin et al., 2005) were followed across all studies.
Demographic information, PTSD symptom severity, and blood pressure levels for each study, stratified by ancestry, are presented in Table 1. There was a mix of military and community-based cohorts, and several of the military cohorts were all-male (e.g., GMRF-QUT, MRS, and VETSA). Although individuals of European and African ancestry were included in these analyses, the overall sample was predominantly of European ancestry due to the large sample size of UKBB. On average, participants were in young and middle adulthood, although GMRF-QUT, UKBB, and VETSA had older average ages. Average PTSD symptoms on the transformed 0–1 scale varied across studies, with the lowest mean symptom levels detected for GMRF-QUT and UKBB and the highest mean symptom levels detected for TRACTS and the three treatment studies: CHOICE, DCS, and OPT. Across studies, mean transformed PTSD symptom levels were higher for African ancestry participants than for European ancestry participants as well. For studies that included both European and African ancestry participants (CHOICE, DCS, HOST, MRS, and OPT), we computed the difference in the mean transformed PTSD symptoms for individuals of African and European ancestry. In a meta-analysis of these difference scores, there was a statistically significant difference such that individuals of African ancestry had higher PTSD symptom levels than individuals of European ancestry (difference = 0.02, SE = 0.00, z = 4.10, p < 1E-4). For the majority of studies, mean SBP and DBP were in the normotensive range. Average blood pressure levels were similar for participants of European and African ancestry. For those studies reporting antihypertensive medication use, prevalence in cohorts ranged from 0% in MRS to 57% in VETSA.
Table 1. Participant characteristics for each contributing study, stratified by ancestry (N = 72,224).
Main Effects of Blood Pressure PGS
Polygenic scores for SBP and DBP were calculated using the Million Veteran Program blood pressure GWAS results (Giri et al., 2019) and were used to predict SBP and DBP, respectively, in the different cohorts. As noted above, the optimal pTs were selected based on results from the largest study: UKBB. The optimal pT for predicting SBP in UKBB was p ≤ 0.02, where PGS explained 2.56% of the variation in SBP (p < 1E-100). Every 1 SD increase of PGS relative to the sample mean was associated with a 3.03-point increase in SBP (SE = 0.07, p < 1E-20). The optimal pT for predicting DBP in UKBB was p ≤ 0.01; PGS explained 1.66% of the variation in DBP (p < 1E-100). For every 1 SD increase of PGS relative to the sample mean, there was a 1.37-point increase in DBP (SE = 0.04, p < 1E-20).
In trans-ethnic meta-analyses, the PGS significantly predicted blood pressure (SBP: β = 2.83, SE = 0.06, p < 1E-20; DBP: β = 1.32, SE = 0.04, p < 1E-20; Figures 1A,B), although there was substantial heterogeneity (Table 2). In the trans-ethnic meta-analyses, effect sizes of the PGS on blood pressure were significantly smaller, albeit still significant and positively related to blood pressure, when UKBB was not included in analyses (SBP: z = −8.17, p < 3.8E-16; DBP: z = −4.27, p < 2.0E-5). Trans-ethnic meta-analyses of studies excluding UKBB indicated a 1.46-point increase in SBP for every 1 SD increase in PGS (SE = 0.18, p = 8.8E-17) and a 0.83-point increase in DBP for every 1 SD increase in PGS (SE = 0.12, p = 2.1E-11; Supplementary Figures 1A,B). No statistically significant heterogeneity was detected in the trans-ethnic meta-analyses of cohorts without UKBB, although there was evidence of low-to-moderate heterogeneity that persisted for DBP. Results of separate meta-analyses of African ancestry and European ancestry participants were similar, overall, to the trans-ethnic meta-analyses without UKBB, although moderate-to-high heterogeneity was observed in the African ancestry meta-analysis for DBP.
Figure 1. Forest plots of trans-ethnic meta-analyses for the main effect of blood pressure polygenic scores on (A) systolic blood pressure and (B) diastolic blood pressure. For each cohort, a square is plotted at the effect estimate value. The size of each plotted square reflects relative precision, where size is inversely proportional to the standard error of a given cohort. Reported effects are beta coefficients and 95% confidence intervals. AA = African ancestry; EA = European ancestry; FE = fixed effects.
Main Effects of PTSD Symptoms
Significant main effects of PTSD symptoms on SBP were detected in three individual cohorts, although the effect directions differed (Table 3). Whereas PTSD symptoms were negatively associated with SBP levels in UKBB, PTSD symptoms were positively associated with SBP levels in European ancestry participants from MRS and in African ancestry participants from the CHOICE/DCS/OPT combined cohort. The trans-ethnic meta-analysis of all cohorts was statistically significant and indicated a negative association between PTSD symptoms and SBP (β = −1.46, SE = 0.44, p = 9.8E-4; Figure 2A). An individual with the highest possible PTSD symptom severity was predicted to have SBP that was 1.46 mmHg lower than an individual with the lowest possible PTSD symptom severity. However, this finding was largely driven by the UKBB results, and there was substantial heterogeneity detected. In the trans-ethnic meta-analysis of all cohorts excluding UKBB, there was a nominally significant positive association between PTSD symptoms and SBP (β = 2.81, SE = 1.13, p = 0.01; Supplementary Figure 2A). The effect sizes for the main effect of PTSD symptoms on SBP levels differed significantly when including vs. excluding UKBB (z = −4.11, p = 4.0E-5). Further, when excluding UKBB in trans-ethnic meta-analyses, heterogeneity was minimal and not statistically significant. A significant positive main effect of PTSD symptoms on SBP was also detected in the European ancestry meta-analysis without UKBB (β = 3.94, SE = 1.39, p = 4.7E-3). No significant main effect of PTSD symptoms was detected for SBP for the African ancestry meta-analysis (β = 0.63, SE = 1.93, p = 0.74).
Figure 2. Forest plots of trans-ethnic meta-analyses for the main effect of posttraumatic stress disorder symptoms on (A) systolic blood pressure and (B) diastolic blood pressure. For each cohort, a square is plotted at the effect estimate value. The size of each plotted square reflects relative precision, where size is inversely proportional to the standard error of a given cohort. Reported effects are beta coefficients and 95% confidence intervals. AA = African ancestry; EA = European ancestry; FE = fixed effects.
For DBP, there were significant main effects of PTSD symptoms in the UKBB and GTP cohorts, but effects were in the opposite direction. In UKBB, PTSD symptoms were positively associated with DBP levels, but PTSD symptoms were negatively associated with DBP levels in GTP. The trans-ethnic meta-analysis was statistically significant and suggested a small positive association between PTSD symptoms and DBP (β = 0.70, SE = 0.26, p = 8.1E-3; Figure 2B). An individual with the highest possible PTSD symptom severity was predicted to have a DBP measure that was 0.70 mmHg higher than an individual with the lowest possible PTSD symptom severity. No statistically significant heterogeneity was detected. In the trans-ethnic meta-analysis excluding UKBB, the association remained positive but was smaller and not statistically significant (β = 0.43, SE = 0.78, p = 0.58; Supplementary Figure 2B). However, the effect sizes for the main effect of PTSD symptoms on DBP levels did not differ significantly when including vs. excluding UKBB (z = 0.37, p = 0.71). A larger positive, albeit non-significant, main effect of PTSD symptoms was observed in the European ancestry meta-analysis without UKBB (β = 1.54, SE = 1.00, p = 0.12), whereas a negative, non-significant main effect of PTSD symptoms was observed in the African ancestry meta-analysis (β = −1.30, SE = 1.25, p = 0.30).
Interactions of PTSD Symptoms and Blood Pressure PGS
Despite nominally significant interaction terms of PTSD symptoms and blood pressure PGS for TRACTS for SBP and for SHRS for DBP, no significant interaction effects emerged in trans-ethnic meta-analyses (Table 4 and Figure 3). Furthermore, there was evidence of low-to-moderate heterogeneity that did not reach the level of statistical significance. For DBP, there was a nominally significant interaction of PTSD symptoms with PGS for the African ancestry meta-analysis (β = −2.44, SE = 1.19, p = 0.04).
Table 4. Interaction effects of PTSD symptoms with polygenic scores on blood pressure, within and across cohorts.
Figure 3. Forest plots of trans-ethnic meta-analyses for the interaction of posttraumatic stress disorder symptoms and blood pressure polygenic scores on (A) systolic blood pressure and (B) diastolic blood pressure. For each cohort, a square is plotted at the effect estimate value. The size of each plotted square reflects relative precision, where size is inversely proportional to the standard error of a given cohort. Reported effects are beta coefficients and 95% confidence intervals. AA = African ancestry; EA = European ancestry; FE = fixed effects.
In the largest study examined – UKBB – we were able to conduct sex-stratified analyses in order to explore potential sex-specific effects (Supplementary Table 1). Women had a stronger positive association of PGS with SBP levels than men, and there was a significant Sex × PGS interaction (β = −0.42, SE = 0.14, p = 2.0E-3). A similar finding was observed for DBP, although the Sex × PGS interaction was only nominally significant (β = −0.20, SE = 0.08, p = 0.01). The negative association between PTSD symptoms and SBP levels was somewhat larger for men than for women, although there was no significant Sex × PTSD Symptoms interaction (β = 1.27, SE = 0.98, p = 0.19). Furthermore, even though PTSD symptoms were significantly positively associated with DBP levels in women only, the Sex × PTSD Symptoms interaction was not significant (β = 0.24, SE = 0.57, p = 0.68). Similar to the overall trans-ethnic meta-analyses, no significant PTSD Symptom × PGS interactions emerged in sex-stratified analyses in UKBB.
The patterns of results from meta-analyses stratified by military or community-based cohort type paralleled those of the trans-ethnic meta-analyses when excluding versus including UKBB, respectively (Supplementary Table 2). Both sets of analyses identified significant positive associations of blood pressure PGS with SBP and DBP, although effects were larger in the community-based cohorts. With respect to the main effect of PTSD symptoms, there was a significant negative effect of PTSD on SBP in the community-based cohorts (β = −2.18, SE = 0.47, p = 3.7E-6) and a significant positive effect of PTSD on SBP in the military cohorts (β = 4.18, SE = 1.41, p = 2.9E-3). For DBP, the positive main effect of PTSD was nominally significant for community-based cohorts (β = 0.64, SE = 0.28, p = 0.02) but not for military cohorts (β = 1.80, SE = 1.02, p = 0.08). No significant PTSD Symptoms × PGS interactions emerged in analyses stratified by military or community-based cohort type.
When restricting analyses to studies that assessed for antihypertensive medication use and were therefore able to adjust blood pressure values for medication use, results of trans-ethnic and single ancestry meta-analyses indicated significant positive main effects of PGS on SBP and DBP (Supplementary Table 3). Additionally, there were significant positive main effects of PTSD on SBP for the trans-ethnic and European ancestry meta-analyses, although no significant effects of PTSD were observed for DBP. No significant PTSD Symptoms × PGS interactions were detected in any of these secondary meta-analyses.
A substantial body of research suggests that PTSD may be a risk factor for poor cardiovascular health (Edmondson and von Känel, 2017; Koenen et al., 2017), and yet our understanding of who might be at greatest risk of adverse cardiovascular outcomes after trauma is limited. In the current study, we conducted the first examination of the individual and synergistic contributions of PTSD symptoms and genetics to continuous blood pressure levels. To explore this research question, we harnessed the power of the collaborative, consortium-based PGC-PTSD Physical Health Working Group and investigated these associations across 11 studies of over 72,000 individuals of European and African ancestry. In trans-ethnic meta-analyses, we identified polygenic main effects based on the results of a trans-ethnic blood pressure GWAS on both SBP and DBP, with higher PGS predicting higher blood pressure levels. Significant main effects of PTSD symptoms were also detected for SBP and DBP in meta-analyses, though there was substantial heterogeneity in the results across cohorts. When including data from the largest contributing study – UKBB – PTSD symptoms were negatively associated with SBP levels and positively associated with DBP levels. However, when excluding UKBB, there was a nominally significant positive association of PTSD symptoms with SBP levels; no significant association was observed for DBP. Furthermore, we did not find evidence that blood pressure PGS significantly moderated the associations between PTSD symptoms and blood pressure levels in meta-analyses.
A notable feature of our study is that we conducted trans-ethnic meta-analyses and used the results from a trans-ethnic blood pressure GWAS from the Million Veteran Program (Giri et al., 2019) to derive our blood pressure PGS. Indeed, compared to other genomics consortia, the PGC-PTSD Group is unique in that participants of diverse ancestries comprise a substantial proportion of the sample (Nievergelt et al., 2019). This aspect of our study thus enhances the generalizability of our findings and is notable given the relative lack of racial and ethnic diversity that is often observed in genetics studies (Popejoy and Fullerton, 2016), which can contribute to health disparities in underrepresented groups. As expected, higher SBP and DBP PGS were each associated with higher SBP and DBP levels, respectively. In their GWAS results, Giri et al. (2019) found that significant sentinel SNPs at independent loci accounted for more variance in SBP (3.56%) than for DBP (1.06%), and we found that effect sizes of PGS on blood pressure levels – although small – were larger for SBP than for DBP in trans-ethnic meta-analyses. However, we found that the PGS based on the best fitting p-value thresholds in UKBB accounted for slightly less variance in SBP and slightly more variance in DBP in this cohort compared to that observed for the sentinel SNPs in the Giri et al. (2019) GWAS. When making these comparisons, it is worth noting though that our approach of generating PGS across thousands of SNPs based on GWAS summary statistics from an independent sample differed from the approach of Giri et al. (i.e., examining variance accounted for by significant sentinel SNPs identified in their own samples). In single ancestry meta-analyses, associations were observed in both European and African ancestry individuals. These findings parallel research on the genetic underpinnings of blood pressure in individuals of diverse backgrounds, which has found similar directions of effects across a number of ancestry groups (Franceschini et al., 2013; Giri et al., 2019). We also found some evidence that PGS for SBP and, to a lesser extent, DBP were more strongly related to the corresponding blood pressure levels in women than in men in the UKBB cohort. Given sex differences in CVD (Virani et al., 2020), such findings are intriguing and merit replication and further investigation.
Although PTSD symptoms were significantly associated with continuous blood pressure levels across studies in trans-ethnic meta-analyses, the results were heterogeneous across studies and for SBP versus DBP. Indeed, directions of the PTSD main effect varied across cohorts, although significant, yet small, main effects were observed in the overall trans-ethnic meta-analyses and were driven by the largest cohort included: UKBB. Interestingly, greater PTSD symptoms were associated with lower SBP and higher DBP levels in trans-ethnic meta-analyses of all cohorts. In particular, the SBP finding was contrary to our hypotheses given the literature documenting associations between PTSD and elevated blood pressure levels and hypertension (Schnurr et al., 2000; Kibler et al., 2009; Sumner et al., 2016; Burg et al., 2017; Edmondson et al., 2018; Howard et al., 2018). However, when excluding the results of UKBB in trans-ethnic and European ancestry meta-analyses, there was a nominally significant association of greater PTSD symptoms with higher SBP levels, and no significant association was observed for DBP. SBP and DBP exhibit different age-related changes over the lifespan; whereas SBP increases continuously with age, DBP has been found to decrease after 50–60 years of age (Franklin et al., 1997). A number of studies suggest that SBP may be a superior predictor of cardiovascular morbidity and mortality compared to DBP (Kannel et al., 1971; Rutan et al., 1989; Borghi et al., 2003). Furthermore, SBP – and not DBP – is incorporated into various global CVD risk algorithms, such as the Framingham Risk Score and Atherosclerotic CVD Risk Calculator (D’Agostino et al., 2008; Muntner et al., 2014).
The disparate results for associations of PTSD symptoms with blood pressure that were obtained when including versus excluding UKBB may be due, in part, to particular aspects of that cohort’s study design. In UKBB, only six symptoms of PTSD were assessed using self-report, which would be expected to limit the variance in PTSD symptoms in this cohort. In addition, no information on antihypertensive medication was available, which would be expected to result in artificially lowered blood pressure readings for those taking such medications. These study characteristics may explain some differences in findings compared to other studies that comprehensively measured symptoms of PTSD and included information on antihypertensive medication use, allowing for adjustment of estimated blood pressure (Tobin et al., 2005). Nevertheless, together these findings suggest that the link between PTSD and high blood pressure may not be robust or may be subject to substantial individual difference factors that we were not able to study as moderators of the PTSD and blood pressure associations. Similarly, in a meta-analysis of associations between PTSD, heart rate, and blood pressure, the effects for blood pressure – although statistically significant – were smaller than those observed for heart rate (Buckley and Kaloupek, 2001). Findings from the epidemiological literature also suggest that the link between PTSD and high blood pressure may be weaker than those with CVD events and diagnoses, such as myocardial infarction, stroke, and heart disease (Sumner et al., 2015, 2016; El-Gabalawy et al., 2018). Some research also suggests that certain PTSD symptom dimensions (e.g., fear-related symptoms) may be more strongly related to high blood pressure than others (Sumner et al., 2020). In future research, a more nuanced approach that includes deep phenotyping of both PTSD and cardiovascular metrics and that examines potential moderators may be needed to better understand how PTSD relates to blood pressure.
Unlike prior research suggesting that PTSD symptom severity may be most strongly related to cardiometabolic risk among those individuals with greater genetic risk for obesity (Wolf et al., 2017), we did not find significant interactions of PTSD and PGS with blood pressure levels in meta-analyses. Thus, the genetic contributions to blood pressure and any PTSD-related contributions to blood pressure were independent of each other, rather than synergistic. Notably, in these models, we accounted for relevant covariates (i.e., age and sex) and interactions of these covariates with our measures of PTSD symptoms and blood pressure PGS. Therefore, it is unlikely that characteristics such as age were obscuring a potential interaction of PTSD symptoms with PGS. However, due to variability in the measures collected across studies, we were unable to examine the influence of other factors such as cumulative trauma burden, time since trauma, or chronicity of PTSD. In the current investigation, we focused on SBP and DBP rather than other metrics like heart rate because these were the most common cardiovascular measures across studies of the PGC-PTSD Physical Health Working Group. However, it is of interest to leverage this model in future work in an effort to better understand risk for adverse physical health consequences among individuals who experience negative psychiatric sequelae of trauma. In particular, this approach may hold promise for informing screening efforts for chronic disease prevention. Given that associations of PTSD with chronic disease outcomes, including CVD, are complex and multifactorial, comprehensive consideration of multiple disease-related factors may be needed to better elucidate risk.
Our findings need to be considered in light of several limitations. First, although we conducted trans-ethnic meta-analyses across 11 different studies, there was variability in the nature of samples (e.g., military vs. community-based samples) and in the measures of PTSD (e.g., interview- vs. questionnaire-based assessments) and blood pressure (e.g., in protocols and monitors). Including a variety of studies increases the external validity but also contributed to the heterogeneity that we observed. Furthermore, this heterogeneity precluded our ability to engage in deep phenotyping and made it difficult to determine if results were due to differences in sample characteristics or methodological factors. Although we considered secondary stratified analyses in military and community-based cohorts to examine whether this study characteristic might explain some of the heterogeneity we observed in overall analyses, it is worth noting that cohort type overlapped with other study factors (e.g., military cohorts were all- or predominantly male and several community-based cohorts were predominantly female). Additionally, it is possible that community-based cohorts included participants who served in the military but did not collect this information and/or that military cohorts included participants who experienced non-military trauma. Second, despite the large sample size, it is possible that we were statistically underpowered to detect PTSD × Blood Pressure PGS interactions. Future research with larger sample sizes is needed. Third, although we focused in particular on individuals of European and African ancestry in the current investigation given the composition of our cohorts, additional work that incorporates individuals from a variety of additional ancestry groups is needed. Fourth, even though we accounted for antihypertensive medication as recommended by Tobin et al. (2005) for most studies in the current analyses, some studies did not assess medication use, and we were unable to determine adherence to these medications even when it was assessed. However, in meta-analyses restricted to those studies that contributed medication-adjusted blood pressure values, we found similar results. Furthermore, previous studies have shown that certain classes of antihypertensive medications (e.g., angiotensin-converting enzyme inhibitors and angiotensin II receptor blockers) are associated with lower PTSD symptoms in trauma-exposed individuals (Khoury et al., 2012; Seligowski et al., 2021), and that the effects of these antihypertensive medications may be dependent on genetic liability (Nylocks et al., 2015). Despite these limitations, our study has several strengths, including generating PGS from a large trans-ethnic GWAS of blood pressure, utilizing continuous measures of PTSD symptom severity, and conducting team science trans-ethnic meta-analyses. Given that much of genetics research has traditionally focused on individuals of European ancestry, which contributes to disparities in understanding of health-relevant factors for individuals from different backgrounds, our study provides an approach for reducing this problem and increasing representation of diverse ancestries.
Posttraumatic stress disorder has been linked to adverse cardiovascular health, and research suggests that shared genetic factors may explain, in part, this association (Sumner et al., 2017; Nievergelt et al., 2019). In the current trans-ethnic meta-analyses of over 72,000 European and African ancestry individuals, PTSD symptoms were associated with SBP and DBP levels but in heterogeneous ways and with small effect sizes. Moreover, although underlying PGS significantly predicted blood pressure levels, this genetic metric did not modify associations of PTSD symptoms with SBP or DBP. Further research is needed to better understand the extent to which PTSD is associated with high blood pressure and how genetic factors might play a role in influencing cardiovascular risk. In addition, future studies merit further investigation of the role of environment and context, along with genetics and PTSD symptoms. In examining determinants of health, social determinants, health behaviors, and access to quality medical care are thought to have a larger impact than heredity, with the largest influence being social factors (Tarlov, 1999; Galea et al., 2011). Thus, it is of interest for future studies to adopt a biopsychosocial approach when considering risk for high blood pressure in individuals who have experienced trauma.
Data Availability Statement
Access to individual-level data for available datasets may be requested through the PGC Data Access Portal at https://www.med.unc.edu/pgc/shared-methods/data-access-portal/. The UKBB data can be obtained by applying here: https://www.ukbiobank.ac.uk/enable-your-research. Other data that contributed to the findings of this study may be available from study Principal Investigators upon request.
The studies involving human participants were reviewed and approved by the Institutional Review Boards of Boston University, Duke University Medical Center, Durham, NC VA, Emory University, Grady Memorial Hospital, Greenslopes Private Hospital, Puget Sound VA Healthcare System, Queensland University of Technology, Stellenbosch University, the Department of Veterans Affairs, University Hospitals, University of California, San Diego, VA Boston Healthcare System, and VA San Diego. The patients/participants provided their written informed consent to participate in this study.
JAS, EW, CN, AM, and MWL conceptualized the study and design. LA, OA, AA-K, DB, JB, BB, GB, JC, AD, MD, NF, CF, MG, CG, GG, MH, SH, TJ, NK, WK, BL, AL, MWL, MJL, JM-K, AM, MRM, ReM, DM, RM, VM, WM, MWM, CM, CN, MP, KR, VR, AR, BR, PR-B, SS, AS, JSS, LV, JV, RY, and LZ assisted with obtaining funding for the included studies, and data collection and preparation. AM conducted the analyses, with assistance from MWL and MG. JAS prepared the first draft of the manuscript. AM, VM, AR, and EW provided edits as members of a writing group. All authors reviewed the manuscript, provided feedback, and approved the manuscript for submission.
This work was supported by the National Heart, Lung, and Blood Institute under award number K01HL130650 to JAS and a Merit Review Award Number I01 CX-001276-01 to EW from the United States Department of Veterans Affairs, Clinical Sciences Research and Development Program. This work was also supported by the Building Interdisciplinary Research Careers in Women’s Health of the National Institutes of Health under Award Number K12HD085850 (JM-K) and UL1TR002378 (Georgia CTSA). Research reported in this publication was supported by the South African Medical Research Council (SAMRC) for the “Shared Roots” Flagship Project, Grant no. MRC-RFA-IFSP-01-2013/SHARED ROOTS through funding received from the South African National Treasury under its Economic Competitiveness and Support Package. The work by LV was supported by funding from the SAMRC through its Division of Research Capacity Development under the SAMRC clinician researcher (MD) scholarship program, funded by the South African National Treasury. VETSA was funded by the National Institute on Aging (award numbers R01AG050595 and R01AG022381). HOST was funded by the National Institute of Mental Health (award number R01MH062482 to JB). The Translational Research Center for TBI and Stress Disorders (TRACTS) was supported by the VA Rehabilitation Research and Development Traumatic Brain Injury Center of Excellence (B9254-C), and the Cooperative Studies Program, Department of Veterans Affairs. Genotyping of the TRACTS cohort was supported by National Institute of Mental Health under award number R21MH102834 to MWM. Analysis of the genetic data was supported by the National Institute of Mental Health under Award Number R01MH106595 to CN. The contents of this research article are the sole responsibility of the authors and do not necessarily represent the official views of the National Institutes of Health, National Heart, Lung, and Blood Institute, National Institute on Aging, National Institute of Mental Health, SAMRC, U.S. Department of Veterans Affairs, U.S. Government, or other funders.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The PGC-PTSD Physical Health Working Group is co-led by JAS and EW. Working group members include Ananda B. Amstadter, DB, JB, Leslie Brick, Karmel W. Choi, Shareefa Dalvie, Terri deRoon-Cassini, Laramie E. Duncan, Bizu Gelaye, Jessica M. Gill, SH, Arash Javanbakht, Amanda J. Khan, NK, Christine L. Larson, ML, Donghao Lu, AM, Samuel A. McLean, DM, VM, Karen S. Mitchell, Peter Naude, CN, Nicole R. Nugent, Sheila A. M. Rauch, AR, SS, Gen Shinozaki, AS, Suraj Thapa, Unnur A. Valdimarsdóttir, Eric Vermetten, and Nagy Youssef. We acknowledge all the participants who contributed to these research studies. Further, for VETSA, the United States Department of Veterans Affairs, Department of Defense; National Personnel Records Center, National Archives and Records Administration; National Opinion Research Center; National Research Council, National Academy of Sciences; and the Institute for Survey Research, Temple University provided invaluable assistance in the creation of the Vietnam Era Twin (VET) Registry. The Cooperative Studies Program of the United States Department of Veterans Affairs provided financial support for development and maintenance of the VET Registry. We would also like to acknowledge the continued cooperation and participation of the members of the VET Registry and their families. This research has been conducted using the United Kingdom Biobank resource under application 41209.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2021.678503/full#supplementary-material
- ^ We examined the genetic correlations of PTSD with SBP and DBP by uploading PTSD GWAS summary statistics from Nievergelt et al. (2019) to LD Hub; GWAS summary statistics for SBP and DBP from UKBB are available on this site. Although statistically significant, the genetic correlation of PTSD with SBP was small (rg = 0.10, p = 4.0E-4); the genetic correlation of PTSD with DBP was not statistically significant (rg = 0.02, p = 0.56). These findings suggest that, by testing PTSD Symptom × Blood Pressure PGS interactions, we were examining Gene × Environment interactions and not environmental factors that were confounded by shared genetics.
- ^ Results of interaction models were similar when including vs. excluding interactions of all covariates with PTSD symptoms and blood pressure PGS. Thus, we present the results of interaction models with these covariate interactions included per the recommendation of Keller (2014).
Borghi, C., Dormi, A., L’italien, G., Lapuerta, P., Franklin, S. S., Collatina, S., et al. (2003). The relationship between systolic blood pressure and cardiovascular risk—results of the Brisighella Heart Study. J. Clin. Hypertens. 5, 47–52. doi: 10.1111/j.1524-6175.2003.01222.x
Buckley, T. C., and Kaloupek, D. G. (2001). A meta-analytic examination of basal cardiovascular activity in posttraumatic stress disorder. Psychosom. Med. 63, 585–594. doi: 10.1097/00006842-200107000-00011
Burg, M. M., Brandt, C., Buta, E., Schwartz, J., Bathulapalli, H., Dziura, J., et al. (2017). Risk for incident hypertension associated with PTSD in military veterans, and the effect of PTSD treatment. Psychosom. Med. 79, 181–188.
Bycroft, C., Freeman, C., Petkova, D., Band, G., Elliott, L. T., Sharp, K., et al. (2018). The UK Biobank resource with deep phenotyping and genomic data. Nature 562, 203–209. doi: 10.1038/s41586-018-0579-z
Chen, C.-Y., Pollack, S., Hunter, D. J., Hirschhorn, J. N., Kraft, P., and Price, A. L. (2013). Improved ancestry inference using weights from external reference panels. Bioinformatics 29, 1399–1406. doi: 10.1093/bioinformatics/btt144
D’Agostino, R. B., Vasan, R. S., Pencina, M. J., Wolf, P. A., Cobain, M., Massaro, J. M., et al. (2008). General cardiovascular risk profile for use in primary care. Circulation 117, 743–753. doi: 10.1161/circulationaha.107.699579
Edmondson, D., Sumner, J. A., Kronish, I. M., Burg, M. M., Oyesiku, L., and Schwartz, J. E. (2018). The association of PTSD with clinic and ambulatory blood pressure in healthy adults. Psychosom. Med. 80, 55–61. doi: 10.1097/psy.0000000000000523
El-Gabalawy, R., Blaney, C., Tsai, J., Sumner, J. A., and Pietrzak, R. H. (2018). Physical health conditions associated with full and subthreshold PTSD in US military veterans: results from the National Health and Resilience in Veterans Study. J. Affect. Disord. 227, 849–853. doi: 10.1016/j.jad.2017.11.058
Franceschini, N., Fox, E., Zhang, Z., Edwards, T. L., Nalls, M. A., Sung, Y. J., et al. (2013). Genome-wide association analysis of blood-pressure traits in African-ancestry individuals reveals common associated genes in African and non-African populations. Am. J. Hum. Genet. 93, 545–554.
Franklin, S. S., Gustin, W. IV, Wong, N. D., Larson, M. G., Weber, M. A., Kannel, W. B., et al. (1997). Hemodynamic patterns of age-related changes in blood pressure: the Framingham Heart Study. Circulation 96, 308–315. doi: 10.1161/01.cir.96.1.308
Galea, S., Tracy, M., Hoggatt, K. J., Dimaggio, C., and Karpati, A. (2011). Estimated deaths attributable to social factors in the United States. Am. J. Public Health 101, 1456–1465. doi: 10.2105/ajph.2010.300086
Giri, A., Hellwege, J. N., Keaton, J. M., Park, J., Qiu, C., Warren, H. R., et al. (2019). Trans-ethnic association study of blood pressure determinants in over 750,000 individuals. Nat. Genet. 51, 51–62.
Howard, J. T., Sosnov, J. A., Janak, J. C., Gundlapalli, A. V., Pettey, W. B., Walker, L. E., et al. (2018). Associations of initial injury severity and posttraumatic stress disorder diagnoses with long-term hypertension risk after combat injury. Hypertension 71, 824–832. doi: 10.1161/hypertensionaha.117.10496
Kannel, W. B., Gordon, T., and Schwartz, M. J. (1971). Systolic versus diastolic blood pressure and risk of coronary heart disease: the Framingham study. Am. J. Cardiol. 27, 335–346. doi: 10.1016/0002-9149(71)90428-0
Keller, M. C. (2014). Gene × environment interaction studies have not properly controlled for potential confounders: the problem and the (simple) solution. Biol. Psychiatry 75, 18–24. doi: 10.1016/j.biopsych.2013.09.006
Khoury, N. M., Marvar, P. J., Gillespie, C. F., Wingo, A., Schwartz, A., Bradley, B., et al. (2012). The renin-angiotensin pathway in posttraumatic stress disorder: angiotensin-converting enzyme inhibitors and angiotensin receptor blockers are associated with fewer traumatic stress symptoms. J. Clin. Psychiatry 73, 849–855. doi: 10.4088/jcp.11m07316
Kibler, J. L., Joshi, K., and Ma, M. (2009). Hypertension in relation to posttraumatic stress disorder and depression in the US National Comorbidity Survey. Behav. Med. 34, 125–132. doi: 10.3200/bmed.34.4.125-132
Koenen, K. C., Sumner, J. A., Gilsanz, P., Glymour, M. M., Ratanatharathorn, A., Rimm, E. B., et al. (2017). Post-traumatic stress disorder and cardiometabolic disease: improving causal inference to inform practice. Psychol. Med. 47, 209–225. doi: 10.1017/s0033291716002294
Kubzansky, L. D., Koenen, K. C., Jones, C., and Eaton, W. W. (2009). A prospective study of posttraumatic stress disorder symptoms and coronary heart disease in women. Health Psychol. 28, 125–130. doi: 10.1037/0278-622.214.171.124
Kubzansky, L. D., Koenen, K. C., Spiro, A., Vokonas, P. S., and Sparrow, D. (2007). Prospective study of posttraumatic stress disorder symptoms and coronary heart disease in the normative aging study. Arch. Gen. Psychiatry 64, 109–116. doi: 10.1001/archpsyc.64.1.109
Muntner, P., Colantonio, L. D., Cushman, M., Goff, D. C., Howard, G., Howard, V. J., et al. (2014). Validation of the atherosclerotic cardiovascular disease Pooled Cohort risk equations. JAMA 311, 1406–1415. doi: 10.1001/jama.2014.2630
Nievergelt, C. M., Maihofer, A. X., Klengel, T., Atkinson, E. G., Chen, C.-Y., Choi, K. W., et al. (2019). International meta-analysis of PTSD genome-wide association studies identifies sex-and ancestry-specific genetic risk loci. Nat. Commun. 10:4558.
Nylocks, K. M., Michopoulos, V., Rothbaum, A. O., Almli, L., Gillespie, C. F., Wingo, A., et al. (2015). An angiotensin-converting enzyme (ACE) polymorphism may mitigate the effects of angiotensin-pathway medications on posttraumatic stress symptoms. Am. J. Med. Genet. B Neuropsychiatr. Genet. 168B, 307–315. doi: 10.1002/ajmg.b.32313
Price, A. L., Patterson, N. J., Plenge, R. M., Weinblatt, M. E., Shadick, N. A., and Reich, D. (2006). Principal components analysis corrects for stratification in genome-wide association studies. Nat. Genet. 38, 904–909. doi: 10.1038/ng1847
Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M. A., Bender, D., et al. (2007). PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81, 559–575. doi: 10.1086/519795
Rutan, G. H., Mcdonald, R. H., and Kuller, L. H. (1989). A historical perspective of elevated systolic vs diastolic blood pressure from an epidemiological and clinical trial viewpoint. J. Clin. Epidemiol. 42, 663–673. doi: 10.1016/0895-4356(89)90010-3
Schnurr, P. P., Spiro, A. III, and Paris, A. H. (2000). Physician-diagnosed medical disorders in relation to PTSD symptoms in older male military veterans. Health Psychol. 19, 91–97. doi: 10.1037/0278-6126.96.36.199
Seligowski, A. V., Duffy, L. A., Merker, J. B., Michopoulos, V., Gillespie, C. F., Marvar, P. J., et al. (2021). The renin–angiotensin system in PTSD: a replication and extension. Neuropsychopharmacology 46, 750–755. doi: 10.1038/s41386-020-00923-1
Sumner, J. A., Duncan, L. E., Wolf, E. J., Amstadter, A. B., Baker, D. G., Beckham, J. C., et al. (2017). Letter to the Editor: posttraumatic stress disorder has genetic overlap with cardiometabolic traits. Psychol. Med. 47, 2036–2039. doi: 10.1017/s0033291717000733
Sumner, J. A., Kubzansky, L. D., Elkind, M. S., Roberts, A. L., Agnew-Blais, J., Chen, Q., et al. (2015). Trauma exposure and posttraumatic stress disorder symptoms predict onset of cardiovascular events in women. Circulation 132, 251–259. doi: 10.1161/circulationaha.114.014492
Sumner, J. A., Kubzansky, L. D., Roberts, A. L., Chen, Q., Rimm, E. B., and Koenen, K. C. (2020). Not all posttraumatic stress disorder symptoms are equal: fear, dysphoria, and risk of developing hypertension in trauma-exposed women. Psychol. Med. 50, 38–47. doi: 10.1017/s0033291718003914
Sumner, J. A., Kubzansky, L. D., Roberts, A. L., Gilsanz, P., Chen, Q., Winning, A., et al. (2016). Posttraumatic stress disorder symptoms and risk of hypertension over 22 years in a large cohort of younger and middle-aged women. Psychol. Med. 46, 3105–3116. doi: 10.1017/s0033291716001914
Tobin, M. D., Sheehan, N. A., Scurrah, K. J., and Burton, P. R. (2005). Adjusting for treatment effects in studies of quantitative traits: antihypertensive therapy and systolic blood pressure. Stat. Med. 24, 2911–2935. doi: 10.1002/sim.2165
Vaccarino, V., Goldberg, J., Rooks, C., Shah, A. J., Veledar, E., Faber, T. L., et al. (2013). Post-traumatic stress disorder and incidence of coronary heart disease: a twin study. J. Am. Coll. Cardiol. 62, 970–978.
Virani, S. S., Alonso, A., Benjamin, E. J., Bittencourt, M. S., Callaway, C. W., Carson, A. P., et al. (2020). Heart disease and stroke statistics—2020 update: a report from the American Heart Association. Circulation 141, e139–e596.
Warren, H. R., Evangelou, E., Cabrera, C. P., Gao, H., Ren, M., Mifsud, B., et al. (2017). Genome-wide association analysis identifies novel blood pressure loci and offers biological insights into cardiovascular risk. Nat. Genet. 49, 403–415.
Wolf, E. J., Miller, D. R., Logue, M. W., Sumner, J., Stoop, T. B., Leritz, E. C., et al. (2017). Contributions of polygenic risk for obesity to PTSD-related metabolic syndrome and cortical thickness. Brain Behav. Immun. 65, 328–336. doi: 10.1016/j.bbi.2017.06.001
Keywords: posttraumatic stress disorder, genetics, blood pressure, trans-ethnic, meta-analysis
Citation: Sumner JA, Maihofer AX, Michopoulos V, Rothbaum AO, Almli LM, Andreassen OA, Ashley-Koch AE, Baker DG, Beckham JC, Bradley B, Breen G, Coleman JRI, Dale AM, Dennis MF, Feeny NC, Franz CE, Garrett ME, Gillespie CF, Guffanti G, Hauser MA, Hemmings SMJ, Jovanovic T, Kimbrel NA, Kremen WS, Lawford BR, Logue MW, Lori A, Lyons MJ, Maples-Keller J, Mavissakalian MR, McGlinchey RE, Mehta D, Mellor R, Milberg W, Miller MW, Morris CP, Panizzon MS, Ressler KJ, Risbrough VB, Rothbaum BO, Roy-Byrne P, Seedat S, Smith AK, Stevens JS, van den Heuvel LL, Voisey J, Young RM, Zoellner LA, PGC-PTSD Physical Health Working Group, Nievergelt CM and Wolf EJ (2021) Examining Individual and Synergistic Contributions of PTSD and Genetics to Blood Pressure: A Trans-Ethnic Meta-Analysis. Front. Neurosci. 15:678503. doi: 10.3389/fnins.2021.678503
Received: 09 March 2021; Accepted: 10 May 2021;
Published: 23 June 2021.
Edited by:Rossen Donev, MicroPharm Ltd., United Kingdom
Reviewed by:Paul J. Marvar, George Washington University, United States
Amy Ferguson, University of Edinburgh, United Kingdom
Copyright © 2021 Sumner, Maihofer, Michopoulos, Rothbaum, Almli, Andreassen, Ashley-Koch, Baker, Beckham, Bradley, Breen, Coleman, Dale, Dennis, Feeny, Franz, Garrett, Gillespie, Guffanti, Hauser, Hemmings, Jovanovic, Kimbrel, Kremen, Lawford, Logue, Lori, Lyons, Maples-Keller, Mavissakalian, McGlinchey, Mehta, Mellor, Milberg, Miller, Morris, Panizzon, Ressler, Risbrough, Rothbaum, Roy-Byrne, Seedat, Smith, Stevens, van den Heuvel, Voisey, Young, Zoellner, PGC-PTSD Physical Health Working Group, Nievergelt and Wolf. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jennifer A. Sumner, email@example.com