Post-traumatic stress and future substance use outcomes: leveraging antecedent factors to stratify risk
1 Department of Social and Behavioral Sciences, Harvard T.H. Chan School of Public Health, Boston, MA, United States
2 Department of Psychiatry and Behavioral Sciences, State University of New York Downstate Medical Center, New York City, NY, United States
3 Department of Anesthesiology, Institute for Trauma Recovery, University of North Carolina at Chapel Hill, Chapel Hill, NC, United States
4 Department of Emergency Medicine, Washington University School of Medicine, St. Louis, MO, United States
5 Department of Epidemiology, Brown University, Providence, RI, United States
6 Department of Emergency Medicine, Brown University, Providence, RI, United States
7 Department of Biostatistics, Gillings School of Global Public Health, University of North Carolina, Chapel Hill, NC, United States
8 Departments of Psychiatry and Neurology, University of California San Francisco, San Francisco, CA, United States
9 Department of Biomedical Informatics, Emory University School of Medicine, Atlanta, GA, United States
10 Department of Biomedical Engineering, Georgia Institute of Technology and Emory University, Atlanta, GA, United States
11 Department of Psychiatry and Behavioral Neurosciences, Wayne State University, Detroit, MI, United States
12 Institute for Technology in Psychiatry, McLean Hospital, Belmont, MA, United States
13 The Many Brains Project, Belmont, MA, United States
14 Department of Psychiatry, Harvard Medical School, Boston, MA, United States
15 Department of Psychology and Neuroscience, University of North Carolina at Chapel Hill, Chapel Hill, NC, United States
16 Department of Sociology, University of North Carolina at Chapel Hill, Chapel Hill, NC, United States
17 Department of Psychiatry, McLean Hospital, Belmont, MA, United States
18 Department of Emergency Medicine, University of Massachusetts Chan Medical School, Worcester, MA, United States
19 Department of Emergency Medicine, Vanderbilt University Medical Center, Nashville, TN, United States
20 Department of Emergency Medicine, Henry Ford Health System, Detroit, MI, United States
21 Department of Emergency Medicine, Indiana University School of Medicine, Indianapolis, IN, United States
22 Department of Emergency Medicine, University of Florida College of Medicine -Jacksonville, Jacksonville, FL, United States
23 Department of Emergency Medicine, Cooper Medical School of Rowan University, Camden, NJ, United States
24 Department of Emergency Medicine, Ohio State University College of Medicine, Columbus, OH, United States
25 Department of Emergency Medicine, Oakland University William Beaumont School of Medicine, Rochester, MI, United States
26 Department of Emergency Medicine, Lewis Katz School of Medicine, Temple University, Philadelphia, PA, United States
27 Department of Emergency Medicine, Emory University School of Medicine, Atlanta, GA, United States
28 Department of Surgery, University of Pennsylvania, Philadelphia, PA, United States
29 Department of Neurosurgery, University of Pennsylvania, Philadelphia, PA, United States
30 Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, United States
31 Department of Surgery, Division of Traumatology, Surgical Critical Care and Emergency Surgery, University of Pennsylvania, Philadelphia, PA, United States
32 Department of Emergency Medicine, Einstein Medical Center, Philadelphia, PA, United States
33 Department of Emergency Medicine, Wayne State University, Ascension St. John Hospital, Detroit, MI, United States
34 Department of Emergency Medicine, Massachusetts General Hospital, Boston, MA, United States
35 Department of Emergency Medicine, Trinity Health-Ann Arbor, Ypsilanti, MI, United States
36 Department of Emergency Medicine, University of Massachusetts Medical School-Baystate, Springfield, MA, United States
37 Department of Emergency Medicine, Wayne State University, Detroit Receiving Hospital, Detroit, MI, United States
38 Department of Emergency Medicine, McGovern Medical School at The University of Texas Health Science Center, Houston, TX, United States
39 Department of Emergency Medicine, Harvard Medical School, Boston, MA, United States
40 Department of Psychological Sciences, University of Missouri - St. Louis, St. Louis, MO, United States
41 Department of Psychology, Yale University, New Haven, CT, United States
42 Department of Anesthesiology, University of Michigan Medical School, Ann Arbor, MI, United States
43 Department of Internal Medicine-Rheumatology, University of Michigan Medical School, Ann Arbor, MI, United States
44 Department of Emergency Medicine, University of North Carolina at Chapel Hill, Chapel Hill, NC, United States
45 Department of Psychiatry, Institute for Trauma Recovery, University of North Carolina at Chapel Hill, Chapel Hill, NC, United States
46 Department of Epidemiology, Harvard T.H. Chan School of Public Health, Harvard University, Boston, MA, United States
*Correspondence: Henri M. Garrison-Desany, hgarrisondesany@hsph.harvard.eduAbstract
Background
Post-traumatic stress disorder (PTSD) and substance use (tobacco, alcohol, and cannabis) are highly comorbid. Many factors affect this relationship, including sociodemographic and psychosocial characteristics, other prior traumas, and physical health. However, few prior studies have investigated this prospectively, examining new substance use and the extent to which a wide range of factors may modify the relationship to PTSD.
Methods
The Advancing Understanding of RecOvery afteR traumA (AURORA) study is a prospective cohort of adults presenting at emergency departments (N = 2,943). Participants self-reported PTSD symptoms and the frequency and quantity of tobacco, alcohol, and cannabis use at six total timepoints. We assessed the associations of PTSD and future substance use, lagged by one timepoint, using the Poisson generalized estimating equations. We also stratified by incident and prevalent substance use and generated causal forests to identify the most important effect modifiers of this relationship out of 128 potential variables.
Results
At baseline, 37.3% (N = 1,099) of participants reported likely PTSD. PTSD was associated with tobacco frequency (incidence rate ratio (IRR): 1.003, 95% CI: 1.00, 1.01, p = 0.02) and quantity (IRR: 1.01, 95% CI: 1.001, 1.01, p = 0.01), and alcohol frequency (IRR: 1.002, 95% CI: 1.00, 1.004, p = 0.03) and quantity (IRR: 1.003, 95% CI: 1.001, 1.01, p = 0.001), but not with cannabis use. There were slight differences in incident compared to prevalent tobacco frequency and quantity of use; prevalent tobacco frequency and quantity were associated with PTSD symptoms, while incident tobacco frequency and quantity were not. Using causal forests, lifetime worst use of cigarettes, overall self-rated physical health, and prior childhood trauma were major moderators of the relationship between PTSD symptoms and the three substances investigated.
Conclusion
PTSD symptoms were highly associated with tobacco and alcohol use, while the association with prospective cannabis use is not clear. Findings suggest that understanding the different risk stratification that occurs can aid in tailoring interventions to populations at greatest risk to best mitigate the comorbidity between PTSD symptoms and future substance use outcomes. We demonstrate that this is particularly salient for tobacco use and, to some extent, alcohol use, while cannabis is less likely to be impacted by PTSD symptoms across the strata.
1Introduction
Traumatic events affect nearly 70% of people in a given year (1). It is estimated that 4.6% of trauma-exposed people report subthreshold post-traumatic stress (PTS) symptoms (2), and 5.6% will experience post-traumatic stress disorder (PTSD) (3). PTS is associated with numerous behavioral health outcomes, including substance use behaviors. Substance use has a higher prevalence rate among those with traumatic event exposures compared to those without (4), but the direction of effects remains unclear. Two primary hypotheses have emerged: individuals who are exposed to a traumatic event may use substances as a way to cope with the stress (5, 6), or individuals who use substances already are more likely to develop PTSD symptoms later (7), potentially through an increased risk of future traumatic events and re-traumatization, such as motor vehicle collision due to intoxication (8). To disentangle this relationship, it is important to use prospective cohort data to follow trends through time and identify differences between incident patterns, defined as new cases of substance use only when previously there was no use, and prevalent patterns, defined as both new and existing behaviors.
The relationship between substance use and PTSD, however, can be influenced by a wide range of factors. Not only is there evidence that sociodemographic factors, such as gender and age, modify this relationship (9–11), but measures of personal resiliency (12) and social support (13) have also been shown to modify this risk. However, it is not clear which of these factors, or others that have previously not been considered, are most important for moderating the risk of substance use with increased PTS symptoms. Understanding this risk stratification is critical for addressing substance use among communities exposed to traumatic events and will aid in developing treatment protocols and guide policy decisions at the population level by helping direct resources to interventions focused on groups at higher risk of substance use after traumatic events.
Prior studies, however, have also lacked prospectively collected data following a traumatic event to stratify this risk longitudinally to fully address this question. Rather, many observational studies examining the comorbidity over time have i) sampled from populations in substance use disorder or PTSD treatment programs (14, 15) and may demonstrate the most severe forms of exposure and outcome; ii) been focused primarily in veteran cohorts, which may not be generalizable to a civilian population (16–20); or iii) been analyzed cross-sectional data when PTSD and substance use were ascertained and therefore lacked a temporal relationship (21–23).
Therefore, to address these gaps in knowledge, we used a national sample of individuals presenting within 72 hours of a traumatic event to 29 emergency departments throughout the country with prospective follow-up, which allowed us to investigate incident and prevalent associations after trauma between substance use and PTSD symptoms. We hypothesized that not only sociodemographic factors but also variables related to psychological resiliency, overall health and sleep, and others may statistically modify the relationship between PTSD and substance use. We hypothesized that these effects will be the most pronounced for new cases of substance use, compared to ongoing substance use.
2Methods
2.1Study sample
The AURORA cohort has been previously described in detail (24). In brief, AURORA is a prospective cohort comprised of nearly 3,000 individuals who present at one of 29 participating emergency departments (EDs) within 72 hours of experiencing a traumatic event. Adults aged 18 to 75 years were excluded if they were administered general anesthesia at the time; experienced long bone fractures, significant hemorrhage, or solid organ injury; or were not alert or oriented at the time of enrollment. Individuals were observed for 1 year, with follow-up at week 2, week 8, month 3, month 6, and month 12 following enrollment in the ED. AURORA included a self-reported questionnaire and biospecimen collection of blood and saliva. Participants also consented to have their medical records included in the data collection.
2.2Measures
Exposure was defined as the raw score of the self-reported PTSD checklist for DSM-5 (PCL-5) (25), reported for the 30 days prior to the event at enrollment ED visit, in the 30 days prior to week 2, week 8, month 3, month 6, and month 12. The PCL-5 was primarily used as a raw score of 0 to 80; however, we also investigated the association of baseline PCL-5 symptoms above the validated threshold of 33 (26). Using the continuous score, we maintained granularity in describing the effects across a spectrum, including subclinical presentations that still affect daily life.
Substance use was defined across three primary outcomes: tobacco, alcohol, and cannabis. These were ascertained as a self-reported count of the frequency of use and as a count of quantity of use in the last 30 days, assessed at all six timepoints using the PhenX Toolkit for evaluation (27). The quantity of cannabis was not collected as part of the PhenX Substance Abuse and Addiction Core Tier 1 questions used in the parent study protocol; therefore, our analysis focused on frequency for comparability across substances.
The following potential confounders were examined: participant age, participant gender identity (defined as cisgender men, cisgender women, and transgender/non-binary people), participant marital status (defined as never married, married, separated/divorced, and widowed/other), participant education (defined as not attending high school, attending only high school, attending college, and attending graduate school), and household income status (according to the following categories: ≤$35,000/year, >$35,000 and ≤$75,000/year, >$75,000/year, and “did not know”). Participant race/ethnicity was defined as Hispanic, non-Hispanic White, non-Hispanic Black/African-American, non-Hispanic Asian, and Native American/American Indian; other races not listed were grouped as a single category.
We investigated 128 variables in total as potential modifiers of the relationship between PTSD and substance use ( Supplementary Table 3 ). The main domains of interest that the study questionnaire covered included the Connor–Davidson Resilience Scale 10 (CD-RISC-10) (28), the Five Facet Mindfulness Questionnaire (FFMQ) as a three-item numeric scale (29), the PROMIS item bank for depression and anxiety (30, 31), the Short Form Survey (SF)-12 to measure overall patient wellbeing (32), the Childhood Trauma Questionnaire Short-Form (CTSQF) (33, 34), the emotional support section of the Perceived Social Support scale (35), the Pittsburgh Sleep Quality Index (36), the Area Deprivation Index 2019 (37) linked via census tract based on a self-reported address at baseline and analyzed as national percentiles reflecting neighborhood deprivation, and the previously reported sociodemographic factors.
2.3Statistical methods
Descriptive statistics were generated for exposure, outcomes, and covariates of interest. Bivariate associations between sociodemographic variables and likely PTSD diagnosis at baseline were tested using an appropriate Student’s t-test for continuous variables and a chi-squared test for categorical variables.
A sensitivity analysis, reported in the Supplementary Material , examined additional substances in the AURORA cohort, including opiates, cocaine, hallucinogens, stimulant drugs, barbiturates, and sedatives, in addition to tobacco, alcohol, and cannabis, using principal component analysis ( Supplementary Figure 1 ). Results indicated a primary factor indexing tobacco, alcohol, and cannabis and a secondary factor indexing all other substances (opiates, cocaine, hallucinogens, stimulant drugs, barbiturates, and sedatives). These additional substances had a low frequency of use overall, with small sample sizes (N’s ranging from 27 to 105). We also assessed correlations between pairwise combinations of substances using a correlogram and found similarly that tobacco, alcohol, and cannabis use frequencies were correlated with one another but not correlated with any other individual substances ( Supplementary Figure 2 ). Therefore, we proceeded with the analysis focused on the factor 1 substances.
We used generalized estimating equations (GEEs) to account for the longitudinal nature of our data, which included multiple observations for each participant. We determined autocorrelation as exchangeable via graphs and estimated the autocorrelation coefficient. We estimated GEE models, controlling for major sociodemographic covariates (gender, age, race/ethnicity, income, and marital status) and two versions of a cross-lagged model.
The first cross-lagged model used PTS symptoms from timepoints week 2, week 8, month 3, and month 6 to predict substance use at timepoints week 8, month 3, month 6, and month 12. Antecedent risk factors were defined at baseline ED visit and week 2. Therefore, at week 2, there was an overlap in risk factor responses and PTS symptom ascertainment, and the model used four timepoints. This model used the maximum number of timepoints while still maintaining the correct temporal order, although there was an overlap between antecedent factors and PTS symptoms in week 2.
The second model used antecedent risk factors at ED baseline and week 2. We used PTS symptoms at week 8, month 3, and month 6. We predicted substance use at month 3, month 6, and month 12. This model had no overlap in PTS symptoms and earlier risk factors and used three timepoints as a more stringent temporal order. We fit both models, given that the former is closer temporally to the index trauma event but has overlapping ascertainment of the risk stratification variables and the exposure, while the latter model may show less association due to increased temporal distance from the index event but maintains greater separation between antecedent and concurrent timepoints. We used the Poisson GEEs with robust standard errors given that counts of frequency and quantity of use were our primary outcomes of interest. We considered using the negative binomial GEE model as well, given that this relaxes dispersion assumptions; however, it could not converge for the data, including when amending the optimizer and attempting different model specifications.
We also compared incident and prevalent substance use and associations with PTS symptoms. “Incident” substance use was defined as having reported no past-month substance use at baseline but with later use at one or more subsequent timepoints. Only past 30-day use was considered due to the very low number of participants reporting no substance use in their lifetime. “Prevalent” substance use was defined as endorsing use at baseline and at least one subsequent timepoint.
2.4Causal forest models
We conducted an honest causal forest analysis to identify the most important factors that stratify the risk of increased substance use due to PTSD symptoms. Causal forests are similar to random forests in that they aggregate a number of causal trees that iteratively maximize the heterogeneity of the average treatment effect across the strata (38). We defined the treatment as PTSD and the outcome as substance use, and we first tested a number of parameters to identify the best method for the causal forests given the data. This included testing a minimum node size of 5, 30, 50, and 100; testing the number of considered variables as default , 10% of p, and 30% of p considered for a given split, with p defined as the number of variables in the full covariate matrix (p = 128); and testing the number of trees as 1,000 trees, 2,000 trees, and 5,000 trees per forest. We found no major differences in fit when using the default , with 2,000 trees and five variables in the node size and proceeded with these parameters. We reported the variable importance for tobacco, alcohol, and cannabis use frequency outcomes. We generated variable importance using the grf package, based on the weighted counts of the proportion of splits on the variable of interest to a depth of 4 (39). We also used the causal forests similarly to generate doubly robust scores to stratify the risk as high versus low and to compare the conditional average treatment effect (CATE) between these strata.
In this way, the CATE can be used similarly to a propensity score. The most notable difference is that it was developed from the causal forest, which seeks to maximize the difference in the relationship between PTSD and substance use between the strata (rather than associated with only the exposure or outcome). Therefore, the CATE as a score represents a stratification tool that similarly maximizes the difference in the effect of PTSD and substance use. We used the median as the cutoff to stratify, theoretically, that the two strata are as maximally apart as the data indicate. We provide a sensitivity analysis using additional cutoff (25th and 75th, and 10th and 90th) percentiles of the CATE as a score in Supplementary Table 4 .
We conducted an omnibus evaluation of the calibration by regressing the scores on the treatment effect, whereby a mean forest prediction nearing 1 indicated good calibration in the mean model and tested the statistical significance of a main treatment effect and a differential forest prediction nearing 1 indicated good calibration of the CATE model for both high/low score strata and tested the statistical significance of a differential treatment effect (39, 40).
2.5Missingness
We assessed missingness and determined that it was likely to be missing at random (MAR) using graphical methods and t-tests for whether a variable was missing (yes/no) and our primary measures of the PCL-5 and frequency of substance use. We conducted multiple imputation by chained equations for 20 datasets across 30 iterations each. We pooled all model estimates using Rubin’s rules (41).
3Results
In our sample, we had 2,943 participants overall, with 1,844 without PCL-5 symptoms indicative of likely PTSD diagnosis and 1,099 with PCL-5 symptom severity indicative of likely diagnosis at baseline, defined as a score of 33 or greater ( Table 1 ). Participants without likely PTSD tended to be younger (mean = 34.8 years) compared to those with likely PTSD (mean = 37.7 years, t-test p-value<0.001). Overall, our sample was primarily non-Hispanic Black/African-American, with 49.5% (N = 1,458), followed by non-Hispanic White/European American (34.7%, n = 1,020) and Hispanic (11.6%, n = 342). The majority of our participants were female assigned at birth (61.8%, n = 1,818) and identified as cisgender female (61.7%, n = 1,815). Being cisgender female was overrepresented in our likely PTSD sample (70.2%, n = 772, p-value <0.001).
| Unlikely PTSD 1
N = 1,844 | Likely PTSD N = 1,099 | Overall N = 2,943 | p-Value 2 | |
|---|---|---|---|---|
| Age, mean (SD) | 34.8 (13.1) | 37.7 (13.5) | 34.9 (11.3) | <0.001 |
| Median (min, max) | 31.5 (18.0, 74.0) | 34.0 (18.0, 73.0) | 32 (18, 74) | |
| Race/ethnicity | ||||
| Hispanic | 207 (11.2%) | 135 (12.3%) | 342 (11.6%) | <0.001 |
| Non-Hispanic Black | 969 (52.5%) | 489 (44.5%) | 1,458 (49.5%) | |
| Non-Hispanic other | 67 (3.6%) | 44 (4.0%) | 111 (3.8%) | |
| Non-Hispanic White | 593 (32.2%) | 427 (38.9%) | 1,020 (34.7%) | |
| Missing | 8 (0.4%) | 4 (0.4%) | 12 (0.4%) | |
| Sex assigned at birth | <0.001 | |||
| Male | 796 (43.2%) | 328 (29.8%) | 1,124 (38.2%) | |
| Female | 1,047 (56.8%) | 771 (70.2%) | 1,818 (61.8%) | |
| Missing | 1 (0.1%) | 0 | 1 (0.0%) | |
| Gender identity | <0.001 | |||
| Male | 795 (43.1%) | 325 (29.6%) | 1,120 (38.1%) | |
| Female | 1,043 (56.6%) | 772 (70.2%) | 1,815 (61.7%) | |
| Transgender | 2 (0.1%) | 2 (0.2%) | 4 (0.1%) | |
| None | 3 (0.2%) | 0 (0%) | 3 (0.1%) | |
| Missing | 1 (0.1%) | 0 (0%) | 1 (0.0%) | |
| Marital status | 0.002 | |||
| Divorced | 252 (13.7%) | 198 (18.0%) | 450 (15.3%) | |
| Married | 375 (20.3%) | 233 (21.2%) | 608 (20.7%) | |
| Never married | 1,175 (63.7%) | 636 (57.9%) | 1,811 (61.5%) | |
| Widowed | 30 (1.6%) | 27 (2.5%) | 57 (1.9%) | |
| Missing | 12 (0.7%) | 5 (0.5%) | 17 (0.6%) | |
| Education | <0.001 | |||
| No HS | 11 (0.6%) | 4 (0.4%) | 15 (0.5%) | |
| Some/finished HS | 1,243 (67.4%) | 661 (60.1%) | 1,904 (64.7%) | |
| Some/finished college | 476 (25.8%) | 331 (30.1%) | 807 (27.4%) | |
| Graduate school | 108 (5.9%) | 100 (9.1%) | 208 (7.1%) | |
| Missing | 6 (0.3%) | 3 (0.3%) | 9 (0.3%) | |
| Income | <0.001 | |||
| ≤$35,000/year | 986 (53.5%) | 658 (59.9%) | 1,644 (55.9%) | |
| >$35,000 and ≤$75,000/year | 349 (18.9%) | 220 (20.0%) | 569 (19.3%) | |
| >$75,000/year | 224 (12.1%) | 144 (13.1%) | 368 (12.5%) | |
| Don’t know | 285 (15.5%) | 77 (7.0%) | 362 (12.3%) | |
| Traumatic events | 0.84 | |||
| Assault | 174 (9.5%) | 113 (10.3%) | 288 (9.8%) | |
| Collision | 1,414 (76.7%) | 833 (75.8%) | 2,247 (76.4%) | |
| Fall | 136 (7.4%) | 77 (7.0%) | 213 (7.2%) | |
| Other | 118 (6.4%) | 75 (6.8%) | 193 (6.6%) | |
| Missing | 1 (0.1%) | 1 (0.1%) | 2(0.1%) | |
| Tobacco | 0.06 | |||
| None | 1,205 (65.3%) | 679 (61.8%) | 1,884 (64.0%) | |
| Any | 632 (34.3%) | 414 (37.7%) | 1,046 (35.5%) | |
| Missing | 7 (0.4%) | 6 (0.5%) | 13 (0.4%) | |
| Alcohol | 0.10 | |||
| None | 721 (39.1%) | 395 (35.9%) | 1,116 (37.9%) | |
| Any | 1,116 (60.5%) | 699 (63.6%) | 1,815 (61.7%) | |
| Missing | 7 (0.4%) | 5 (0.5%) | 12 (0.4%) | |
| Cannabis | 0.90 | |||
| None | 1,290 (70.0%) | 774 (70.4%) | 2,064 (70.1%) | |
| Any | 539 (29.2%) | 319 (29.0%) | 858 (29.2%) | |
| Missing | 15 (0.8%) | 6 (0.5%) | 21 (0.7%) | |
Most of our participants had experienced a collision of some kind (primarily motor vehicle), representing 76.4% (n = 2,247) of our sample ( Supplementary Table 1 ). There was no statistical evidence of a difference in index traumatic events between participants with likely PTSD and those without likely PTSD (p = 0.84). Notably, most participants reported alcohol use in the past 30 days prior to the ED visit (61.7%, n = 1,815) and were considered prevalent alcohol use cases, although this did not differ by PTSD symptoms (p = 0.10). Many also used cannabis in the prior 30 days (29.2%, n = 858) and/or smoked tobacco (35.5%, n = 1,046). These were defined as prevalent cannabis and prevalent tobacco groups, respectively. For incident use (defined as no use at baseline but later use at one or more subsequent timepoints), 162 (10.0% of those not smoking at baseline) participants were incident tobacco users, 296 (30.4% of those not drinking at baseline) were incident alcohol users, and 141 (7.9% of those not using cannabis at baseline) were incident cannabis users.
When examining substance use outcomes at week 8, month 3, month 6, and month 12 in the four-timepoint lagged model ( Table 2 ), we found significant associations between PTSD symptoms and future tobacco use (incidence rate ratio (IRR): 1.003, 95% CI: 1.000, 1.005, p = 0.02) and alcohol use (IRR: 1.002, 95% CI: 1.000, 1.004, p = 0.03). There was no significant association with future cannabis use (IRR: 1.002, 95% CI: 0.999, 1.005, p = 0.13). PTSD symptoms were also associated with future quantity of tobacco smoked (IRR: 1.005, 95% CI: 1.001, 1.01) and quantity of alcohol consumed (IRR: 1.003, 95% CI: 1.001, 1.005). Similar patterns in the three-timepoint lagged model were demonstrated when we used only month 3, month 6, and month 12 time periods ( Supplementary Table 2 ).
| Tobacco frequency 2 | Tobacco quantity 3 | Alcohol frequency | Alcohol quantity 3 | Cannabis frequency | ||||||
|---|---|---|---|---|---|---|---|---|---|---|
| Incidence rate (95% CI) | p-Value | Incidence Rate (95% CI) | p-Value | Incidence rate (95% CI) | p-Value | Incidence rate (95% CI) | p-Value | Incidence rate (95% CI) | p-Value | |
| PTSD symptoms 1 | 1.003 (1, 1.005) | 0.02 | 1.005 (1.001, 1.01) | 0.01 | 1.002 (1.000, 1.004) | 0.03 | 1.003 (1.001, 1.01) | 0.001 | 1.002 (0.999, 1.005) | 0.13 |
| Time | 0.966 (0.938, 0.996) | 0.03 | 0.98 (0.93, 1.02) | 0.29 | 0.99 (0.96, 1.012) | 0.28 | 0.97 (0.94, 0.997) | 0.03 | 0.988 (0.949, 1.029) | 0.55 |
| Marital status | ||||||||||
| Married (Ref) | ||||||||||
| Never married | 1.09 (0.87, 1.36) | 0.46 | 1.09 (0.87, 1.36) | 0.46 | 1.15 (1.01, 1.3) | 0.04 | 1.14 (1.01, 1.28) | 0.03 | 1.42 (1.14, 1.78) | 0.002 |
| Divorced | 1.46 (1.14, 1.88) | 0.003 | 1.46 (1.14, 1.88) | 0.003 | 0.997 (0.85, 1.17) | 0.97 | 1.12 (0.98, 1.29) | 0.11 | 1.35 (1.04, 1.75) | 0.02 |
| Widowed/other | 1.002 (0.62, 1.61) | 0.99 | 1.002 (0.62, 1.61) | 0.99 | 1.2 (0.88, 1.64) | 0.24 | 1.17 (0.88, 1.55) | 0.29 | 1.48 (0.91, 2.38) | 0.11 |
| Gender | ||||||||||
| Cisgender female (Ref) | ||||||||||
| Cisgender male | 1.48 (1.27, 1.72) | <0.001 | 1.48 (1.27, 1.72) | <0.001 | 1.18 (1.08, 1.3) | <0.001 | 1.25 (1.15, 1.36) | <0.001 | 1.35 (1.19, 1.53) | <0.001 |
| Transgender | 2.99 (0.67, 13.31) | 0.15 | 2.99 (0.67, 13.31) | 0.15 | 0.87 (0.43, 1.73) | 0.68 | 1.31 (0.69, 2.46) | 0.41 | 2.05 (0.77, 5.48) | 0.15 |
| Race/ethnicity | ||||||||||
| Hispanic (Ref) | ||||||||||
| Non-Hispanic Black | 1.16 (0.97, 1.38) | 0.10 | 0.7 (0.24, 2.03) | 0.51 | 0.78 (0.37, 1.66) | 0.52 | 0.75 (0.43, 1.32) | 0.32 | 0.97 (0.79, 1.17) | 0.72 |
| Non-Hispanic other | 1.02 (0.72, 1.45) | 0.91 | 0.55 (0.18, 1.75) | 0.31 | 0.68 (0.31, 1.51) | 0.35 | 0.63 (0.35, 1.14) | 0.13 | 0.8 (0.55, 1.16) | 0.23 |
| Non-Hispanic White | 1.44 (1.2, 1.73) | <0.001 | 1.04 (0.35, 3.05) | 0.94 | 0.86 (0.4, 1.83) | 0.69 | 0.76 (0.43, 1.33) | 0.33 | 1.04 (0.85, 1.28) | 0.71 |
| Age | 0.998 (0.993, 1.003) | 0.39 | 1.004 (0.997, 1.01) | 0.25 | 0.999 (0.995, 1.003) | 0.69 | 0.99 (0.99, 0.996) | <0.001 | 0.98 (0.98, 0.99) | <0.001 |
| Income | ||||||||||
| >$35k and ≤$75k (Ref) | ||||||||||
| ≤$35k | 1.38 (1.18, 1.61) | <0.001 | 1.45 (1.17, 1.8) | <0.001 | 0.91 (0.8, 1.02) | 0.12 | 0.95 (0.86, 1.06) | 0.37 | 1.17 (0.97, 1.4) | 0.09 |
| >$75k | 0.66 (0.51, 0.85) | 0.002 | 0.46 (0.32, 0.67) | <0.001 | 1.34 (1.15, 1.56) | <0.001 | 1.11 (0.98, 1.26) | 0.11 | 0.8 (0.6, 1.07) | 0.13 |
| Did not report | 1.2 (0.99, 1.46) | 0.07 | 1.24 (0.93, 1.64) | 0.14 | 1.06 (0.91, 1.23) | 0.45 | 0.97 (0.85, 1.11) | 0.71 | 1.13 (0.9, 1.42) | 0.30 |
There were some differences between the effect of prevalent and incident substance use among individuals who did not use the given substance in the past 30 days prior to the emergency department baseline visit when examining substance use outcomes at month 3, month 6, and month 12 ( Figure 1 ). Increased PTSD symptoms were associated with tobacco frequency (IRR: 1.004, 95% CI: 1.001, 1.01) and quantity (IRR: 1.01, 95% CI: 1.002, 1.01) among prevalent tobacco smokers but were not significantly associated among incident smokers who did not smoke at ED baseline (95% CI for frequency: 0.998, 1.01; 95% CI for quantity: 0.996, 1.01). Alcohol frequency did not show significant associations with PTSD, but alcohol quantity was associated with PTSD among both prevalent and incident drinkers. Incident alcohol consumption had a slightly higher point estimate association with increased PTSD symptoms compared to prevalent consumption (IRR: 1.004 vs. IRR: 1.003), although the respective confidence intervals overlapped. Cannabis showed a slight association among prevalent cases (IRR: 1.003, 95% CI: 1.0001, 1.01), although this was not significant for incident cannabis users.
Across causal forests for tobacco, alcohol, and cannabis frequency using the three-timepoint lagging, we found that lifetime worst cigarette use was a major predictor of all causal forests ( Figure 2 ). The SF-12, which measures overall health effects on quality of life, identified at baseline and week 2 was also important for all substances. For alcohol, lifetime physical abuse, the total score of the childhood trauma questionnaire, and lifetime years of substance use were the next most important factors. For tobacco, resiliency as measured by the CD-RISC-10 was the third most important predictor, and education was the fifth most important predictor, and they were not identified as important for either alcohol or cannabis. For cannabis, race and ethnicity were important stratification factors, which were not identified as important above the mean importance factor for the other substances.
Participants were stratified as low- and high-risk strata using the causal forests similar to propensity scores using the forest prediction probabilities for likely vs. unlikely substance use ( Table 3 ). There was a significant difference in the conditional average treatment effect for tobacco frequency when stratifying low- and high-risk participants (difference in CATE: 0.03, 95% CI: 0.01, 0.06); however, this was not significant for alcohol (difference in CATE: 0.02, 95% CI: −0.001, 0.04) or cannabis (difference in CATE: 0.004, 95% CI: −0.01, 0.02). Bias estimates calculated as per Athey et al. (42) are reported in Supplementary Figure 3 , which were all centered around the null.
| Mean forest effect 1 | p-Value | Differential forest effect 2 | p-Value | ATE among high strata | ATE among low strata | Difference in PTSD effect between strata | |
|---|---|---|---|---|---|---|---|
| Tobacco | 0.28 | 0.38 | 0.29 | 0.001 | 0.03 | −0.01 | 0.03 (0.01, 0.06) |
| Alcohol | 0.94 | 0.04 | 0.74 | 0.07 | 0.001 | −0.02 | 0.02 (−0.0005, 0.04) |
| Cannabis | 1.07 | 0.12 | 0.17 | 0.06 | 0.006 | 0.002 | 0.004 (−0.01, 0.02) |
4Discussion
Our study used a novel cohort of trauma-exposed adults to understand the impact of post-traumatic symptoms on substance use behaviors. We found that PTS symptom scores were highly associated with future tobacco and alcohol use after controlling for sociodemographic factors, but there was little evidence of a relationship with future cannabis use. We also found that prior substance use behaviors were the most important stratification factor of the risk of future substance use. Even for alcohol or cannabis frequency, lifetime reported worst cigarette use was the most important predictor, followed by self-rated overall health at ED and week 2 visits. Tobacco frequency demonstrated the greatest difference in the CATE between the strata determined by the causal forest-predicted probability. Our study has important implications for clinical providers to identify areas of intervention for patients who may have comorbid substance use and PTSD symptoms, as well as avenues for future research into these differential effects of PTSD on behavioral health.
Tobacco showed differences in prevalent use compared to incident use when modeling both frequency and quantity, suggesting that prevalent users are at the highest risk of changes in tobacco consumption with PTSD symptoms. Tobacco also demonstrated the largest CATE based on the causal forest stratification, although the calibration demonstrated that the fit may be insufficient in this causal forest, as demonstrated by the mean forest predictions and differential forest predictions both being less than 1. This fit was not improved when toggling parameters (such as number of trees and node size). Lifetime worst use of cigarettes (e.g., a higher number of cigarettes used when describing one’s most prolific lifetime period of smoking) was the most important variable for tobacco as well as alcohol and cannabis causal forests. Prior tobacco consumption behaviors and initiation make it more likely for participants to use substances again. Intervening on tobacco use may benefit the overall amelioration of substance use, although this should be investigated directly through experimental studies. As a potential area of intervention, while some treatment protocols for comorbid PTSD and tobacco dependence have described addressing both tobacco use and PTSD symptoms concurrently, these have not necessarily been tailored to stress response or the underlying causes of this association (43). This dual-treatment protocol has been demonstrated as effective in military populations (44); however, it has not been clear to what extent this is generalizable to civilian populations.
Overall health as measured by the 12-item Short Form Health Survey (32) and resiliency measured by the Connor–Davidson Resilience Scale-10 (28) were some of the highest importance variables for the relationship between PTSD symptoms and tobacco use. Given the physical effects of tobacco smoking over time, the importance of physical health may be related to the length of smoking. For example, individuals smoking for decades may experience worse health and also report higher past 30-day tobacco use compared to those who do not, and individuals experiencing higher rates of distress related to PTSD symptoms have also reported worse physical health in longitudinal studies. In one longitudinal study, there was evidence that PTSD symptoms mediate the relationship between traumatic event exposures (such as the number of traumatic events experienced) and poor physical health, as well as the number of traumatic events and substance use (45).
For alcohol quantity used, we saw consistent associations with the PCL-5 symptoms. In contrast, alcohol use frequency did not demonstrate consistent associations. The ATE was significant in the causal forests, suggesting that there was an overall effect of PTSD on alcohol quantity, but that there was no major differential effect in our sample, given that the CATE difference was not significant. This may be related to the rate of alcohol use in our sample, whereby it is not possible to differentiate between groups because the risk remains high overall. This may be related to selection bias in our sample, as many individuals presenting at the ED did so after a motor vehicle accident, potentially related to intoxicated driving. We did not directly ask participants if they were intoxicated during these accidents in order to reduce social desirability bias in baseline alcohol consumption responses, but future studies with more diverse traumatic index events in their sample may be better poised to identify relationships with alcohol. Notably, therefore, alcohol may precipitate the event that leads to PTSD symptoms, and those PTSD symptoms may precipitate future alcohol use. Multiple studies have identified that this relationship is likely to be bidirectional (46). It is hypothesized that different risk pathways may be involved in which develops first (47). For chronic pain, for instance, this is considered the “mutual maintenance” hypothesis, whereby distress and disability increase both chronic pain and PTSD, and a host of internal psychological factors further this relationship (48). No study to date has directly tested the mutual maintenance of alcohol use and PTSD symptoms; however, it stands to reason that a similar phenomenon can occur.
It was notable that physical abuse and childhood trauma variables were in the top 5 causal modifiers. Childhood trauma especially has been previously investigated as an exacerbating factor for alcohol dependence later in life (49, 50), and it is notable that this was more important than more proximate adult variables. This suggests that while there is evidence that facets such as mindfulness (51–53) or resiliency (12, 54) are related to alcohol use, these are unlikely to mitigate this relationship in our sample. However, in a high-risk sample, these moderating factors may have only slight effects overall. Future studies should investigate prior trauma in the context of other proximate psychological modifiers in samples with more diverse substance use behaviors.
Our study found no evidence that PTSD symptoms affected future cannabis use. This comports with a recent study from the AURORA cohort that there was little difference in the relationship between PTSD symptoms and groups already using cannabis (55). There is a growing interest in the therapeutic uses of cannabis for a range of mental health disorders, including PTSD, as evidenced by the increased legalization of cannabis in the United States. Notably, PTSD diagnosis is a qualifying condition for citizens to obtain a medical marijuana license within 29 states as of 2022 (56). Under this dynamic, we would expect that prior PTSD symptoms would increase future cannabis use if participants were treating their symptoms with cannabis. To this hypothesis, a systematic review of 10 papers found suggestive but overall low-quality evidence of cannabis’ effectiveness in reducing PTSD symptoms (57).
For instance, in a randomized controlled trial of different cannabis formulations, there was no difference in PTSD symptoms between cannabis treatment and placebo groups (58). However, in an observational study that observed participants for a year, they did report reduced PTSD symptoms associated with self-reported cannabis use (59). One of the largest studies of 2,276 veterans in intensive PTSD treatment programs found an increase in PTSD symptomology among those who started or continued to use cannabis after discharge from the treatment program (60). Our findings, however, represent a large civilian cohort not in intensive treatment programs, and therefore, more research should investigate the generalizability between our study and others and how individuals may relate to cannabis for therapeutic purposes.
Notably, the timing of cannabis use vs. PTSD symptoms is also critical in interpreting these findings. We lagged PTSD symptoms behind cannabis exposure in order to directly test a potential self-medication hypothesis, where one would expect increased prior PTSD symptoms associated with future cannabis use. If, instead, current cannabis use decreases PTSD symptoms, we would see a negative effect; in our sensitivity analysis using cross-sectional data (e.g., PTSD symptoms and cannabis use ascertained at the same timepoint), we continued to not see an association, however. This suggests that the association between PTSD symptoms and cannabis use may be subject to confounding or mediation by other factors in our sample that we have controlled for while other studies have not, or other studies may have over-adjusted models and opened epidemiological collider stratification biases not subject in our study. In post-hoc analyses, we had the power to detect reasonable effect sizes, which suggests that it is unlikely our findings represent a false negative alone.
Our findings demonstrate that different substances may have different relationships with PTSD symptoms. This comports with and expands upon prior findings (61–63). For example, a 2018 study found that different symptom types for PTSD had specific associations, alcohol use was associated with avoidance symptoms, and hyperarousal symptoms were elevated among participants with drug use disorders (64). It further stands to theoretically reason that, due to different subjective effects of substances, those with hyperarousal symptoms primarily may be drawn to a particularly depressant substance that results in greater feelings of numbing or hypoarousal, while those with avoidant symptoms may seek another subjective response. While our study aggregated across symptom dimensions using a validated questionnaire, future studies may select specific dimensions of PTSD symptoms and their associations with substance use behaviors. This would provide information that could personalize intervention recommendations for patients even further.
In examining psychosocial factors, which are often the basis for interventions (e.g., mindfulness and meditation training techniques, and recommendations for building supportive relationships to mitigate symptoms), we found that they did not demonstrate the most important associations. Rather, areas of general health, childhood traumatic events, and even sleep quality showed some of the highest associations before areas such as resiliency and mindfulness. This suggests that simply building cognitive or psychological resiliency may be better supported with concurrent physical health considerations and sleep protocols. Physical wellbeing has been previously identified as an important resilience factor (65–67), as well as larger sociodemographic factors including income and education (51, 54, 68, 69). Our findings highlight the holistic nature of trauma and substance use, and their effects on the body and behavior cannot be ignored in favor of a purely mental understanding.
Our use of causal forests also demonstrates a novel approach to investigating effect modification. With a large cohort, we have a sufficient sample size to fit the generalized random forest algorithm and estimate conditional average treatment effects. To our knowledge, only one study to date has previously applied this method related to post-traumatic stress, which focused on a prospective cohort of older Japanese adults after the 2011 Great East Japan Earthquake (70). This paper may not be generalizable to other populations who have not experienced a single mass-casualty incident, and they used 51 predictors of the CATE, compared to our 128. Using a “bottom-up” approach may identify potential effects that would not be investigated from a “top-down” theory-driven approach alone based on established literature. Notably, the use of forests is similar to allowing for multi-order interactions (such as three-way and four-way) based on the splitting algorithm; however, by aggregating over thousands of forests, it reduces issues of overfitting that may otherwise arise (39). However, it is not without its limits: as a non-parametric estimation, it may have issues with confidence interval construction based on undersmoothing versus a bias-correction method and, similar to more standard methods, experience edge effects in slope estimates.
Our study had a number of limitations, including its reliance primarily on self-report. Toxicological testing was not available for participants at the time of their recruitment, and we relied on self-reported questionnaires, which may be subject to social desirability bias (71). The majority of our participants were recruited following a motor vehicle accident in which they were the driver, which may reduce the number of substance use endorsements in our sample, which would bias our findings toward null hypotheses. Future research may consider toxicological testing at recruitment and follow-up via urine or hair samples to ascertain outcomes. Despite this, we still reported high levels of substance use, which suggests that many did feel comfortable reporting their past 30-day use.
While our sample included many different index traumatic events for eligibility, the majority of our sample represented motor vehicle accident cases, which may not be generalizable to samples with other traumatic events, such as mass casualty events or assaults. This was likely related to our recruitment strategy from EDs, and future studies may consider sampling strategies that prioritize other trauma-exposed adults to understand whether these findings replicate in these additional samples. Our sample is unlikely to reflect individuals who do not experience a traumatic event resulting in physical injury, although our PTSD assessment did not require that symptoms be solely related to the index traumatic event. Therefore, while our PTSD metric may capture some effect of other traumatic events that participants may have experienced, this remains a limitation of the sample and our findings.
There also may be measurement bias in our definition of PTSD and other psychosocial metrics. Given known limitations in the DSM-5 (72) and other measures (73), which may not capture all dimensions of PTSD, depression, anxiety, or other constructs that we considered in our causal forest analysis, our findings may be biased. While we examined many factors, there may be additional unmeasured confounding. Finally, we were limited by the sample size of other substances and were unable to investigate whether these findings held also for opioids, stimulants, etc.
Our study also had numerous strengths; primarily, the cohort represents a national civilian population with a range of psychosocial data and prospective follow-ups. Prior longitudinal cohorts investigating both PTSD and SU to date have been made up of participants primarily in treatment for PTSD or SU (74–76) or veteran samples (77, 78). This allowed us to understand the effect of these various constructs over time and track changes in the frequency of use with greater granularity than had we been limited to binarized variables. We also considered the three most popular substances in the United States in the same sample, allowing for comparisons to be made between substances. Finally, we were able to include the largest reported number of psychosocial factors as potential modifiers using a novel statistical method, which allowed for a comprehensive understanding of our participants compared to smaller cohorts.
The public health implications of this work may benefit directing relevant interventions to patients with the highest risk of PTSD symptoms and future substance use behaviors. We identified prior substance use behaviors, overall health, and childhood traumatic events as highly important. Future studies should consider integrating these screening factors and examine if they can effectively stratify risk in a clinical setting to benefit patient outcomes. Identifying people at risk early may direct interventions for primary and secondary prevention against future or worsening substance use after a traumatic event. While informing patients they should be careful of any substance behaviors, providing insight into specific substances may better address their needs beyond generalized advice.
5Conclusion
We demonstrated differential effects of PTSD on future tobacco use using causal forest modeling, with prior cigarette use being the most important factor. We also found that future alcohol use increases with increased PTSD symptoms. Finally, we found that future cannabis frequency of use did not demonstrate a consistent relationship with PTSD symptoms, suggesting that there may not be a clear association when controlling for potential confounders or when accounting for multiple moderators. Taken together, these findings show that substance use following a traumatic event can vary based on personal substance use histories and depending on how the substance use itself is measured. This has important implications both for future epidemiological research that may consider different metrics of substance use and for clinicians interested in identifying exacerbating factors that may moderate this relationship in their patient populations.
In the future, we recommend similar investigations with sample sizes spanning a wider range of traumatic events that would allow for sufficient subgroup analyses. In our sample, the focus on motor vehicle incidents may limit generalizability or obscure event-specific relationships with substance use behaviors. We also recommend studies that examine differences between substance use disorders and subclinical behaviors and subscale-specific PTSD evaluations; this would essentially create a matrix of potential facets of PTSD that then may be associated with subclinical vs. clinical disordered substance use, providing even more tailored insights. Future work should build upon findings of the most important variables to develop risk stratification methods that may eventually aid in personalized recommendations to patients following a traumatic event to mitigate their individual substance use risks.
Data availability statement
The datasets analyzed for this study can be found in the NIMH data archive [https://nda.nih.gov/edit_collection.html?id=2526].
Ethics statement
The studies involving humans were approved by Harvard Longwood Campus Institutional Review Board. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.
Conflict of interest
KK’s research has been supported by the Robert Wood Johnson Foundation, the Kaiser Family Foundation, the Harvard Center on the Developing Child, Stanley Center for Psychiatric Research at the Broad Institute of MIT and Harvard, the National Institutes of Health, One Mind, the Anonymous Foundation, and Cohen Veterans Bioscience. She has been a paid consultant for Baker Hostetler, Discovery Vitality, and the Department of Justice. She has been a paid external reviewer for the Chan Zuckerberg Foundation, the University of Cape Town, and Capita Ireland. She has had paid speaking engagements in the last three years with the American Psychological Association, European Central Bank. Sigmund Freud University – Milan, Cambridge Health Alliance, and Coverys. She receives royalties from Guilford Press and Oxford University Press. TN has received research support from NIH, VA, and Rainwater Charitable Foundation, and consulting income from Jazz Pharmaceuticals. In the last three years, GDC has received research funding from the NSF, NIH and LifeBell AI, and unrestricted donations from AliveCor Inc, Amazon Research, the Center for Discovery, the Gates Foundation, Google, the Gordon and Betty Moore Foundation, MathWorks, Microsoft Research, Nextsense Inc, One Mind Foundation, the Rett Research Foundation, and Samsung Research. GDC has financial interest in AliveCor Inc and Nextsense Inc. He also is the CTO of MindChild Medical and CSO of LifeBell AI and has ownership in both companies. These relationships are unconnected to the current work. SR reports grants from NIH during the conduct of the study; personal fees from SOBP Society of Biological Psychiatry paid role as secretary, other from Oxford University Press royalties, other from APP American Psychiatric Publishing Inc. royalties, other from VA Veterans Administration per diem for oversight committee, and other from Community Psychiatry/Mindpath Health paid board service, including equity outside the submitted work; other from National Association of Behavioral Healthcare for paid Board service; other from Springer Publishing royalties; and Leadership roles on Board or Council for SOBP, ADAA Anxiety and Depression Association of America, and NNDC National Network of Depression Centers. SS has received funding from the Florida Medical Malpractice Joint Underwriter’s Association Dr. Alvin E. Smith Safety of Healthcare Services Grant; Allergan Foundation; the NIH/NIA-funded Jacksonville Aging Studies Center JAX-ASCENT; R33AG05654; and the Substance Abuse and Mental Health Services Administration 1H79TI083101-01; and the Florida Blue Foundation. CJ has no competing interests related to this work, though he has been an investigator on studies funded by AstraZeneca, Vapotherm, Abbott, and Ophirex. JJ receives consulting payments from Janssen Pharmaceuticals. SEH has no competing interests related to this work, though in the last three years he has received research funding from Aptinyx and Arbor Medical Innovations, and consulting payments from Aptinyx, Heron Therapeutics, and Eli Lilly. SM served as a consultant for Walter Reed Army Institute for Research and for Arbor Medical Innovations.
The remaining 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 handling editor Y-lT declared a shared affiliation with the author(s) GDC.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpsyt.2024.1249382/full#supplementary-material