Abstract
Significant trauma histories and post-traumatic stress disorder (PTSD) are common in persons with substance use disorders (SUD) and often associate with increased SUD severity and poorer response to SUD treatment. As such, this sub-population has been associated with unique risk factors and treatment needs. Understanding the distinct etiological profile of persons with co-occurring SUD and PTSD is therefore crucial for advancing our knowledge of underlying mechanisms and the development of precision treatments. To this end, we employed supervised machine learning algorithms to interrogate the responses of 160 participants with SUD on the multidimensional NIDA Phenotyping Assessment Battery. Significant PTSD symptomatology was correctly predicted in 75% of participants (sensitivity: 80%; specificity: 72.22%) using a classification-based model based on anxiety and depressive symptoms, perseverative thinking styles, and interoceptive awareness. A regression-based machine learning model also utilized similar predictors, but failed to accurately predict severity of PTSD symptoms. These data indicate that even in a population already characterized by elevated negative affect (individuals with SUD), especially severe negative affect was predictive of PTSD symptomatology. In a follow-up analysis of a subset of 102 participants who also completed neurocognitive tasks, comorbidity status was correctly predicted in 86.67% of participants (sensitivity: 91.67%; specificity: 66.67%) based on depressive symptoms and fear-related attentional bias. However, a regression-based analysis did not identify fear-related attentional bias as a splitting factor, but instead split and categorized the sample based on indices of aggression, metacognition, distress tolerance, and interoceptive awareness. These data indicate that within a population of individuals with SUD, aberrations in tolerating and regulating aversive internal experiences may also characterize those with significant trauma histories, akin to findings in persons with anxiety without SUD. The results also highlight the need for further research on PTSD-SUD comorbidity that includes additional comparison groups (i.e., persons with only PTSD), captures additional comorbid diagnoses that may influence the PTSD-SUD relationship, examines additional types of SUDs (e.g., alcohol use disorder), and differentiates between subtypes of PTSD.
Similar content being viewed by others
Introduction
The association between substance use disorders (SUDs) and post-traumatic stress is well-established [1,2,3,4]. Approximately 30–60% of persons seeking treatment for SUDs also meet diagnostic criteria for post-traumatic stress disorder (PTSD) [5, 6], and those with comorbid presentations typically demonstrate a more severe clinical profile [7] with poorer treatment outcomes [8, 9]. Several functional explanations for this association have been put forth [10]. For instance, substance use may serve as a “self-medication” strategy for acutely attenuating post-traumatic stress [11], while ironically worsening emotional health and exacerbating PTSD symptoms [12]. To interrupt this negative cycle, integrative treatments have been developed that focus on mitigating maladaptive avoidance behaviors and improving coping skills [13,14,15]. Unfortunately, while these integrative treatments effectively reduce PTSD symptoms, they show limited efficacy to reduce substance use [16,17,18,19,20].
A clear understanding of the hallmark risk factors and distinct endophenotypic profiles is needed to advance treatment strategies for trauma-exposed SUD populations. Various neurobehavioral risk factors have been linked with comorbid SUD and PTSD, including increased anxiety sensitivity [21], anhedonia [22, 23], impulsivity [24,25,26,27,28], more severe executive functioning deficits [29,30,31], and increased sensitivity to reward and punishment [32]. Yet, while such evidence supports potential etiological pathways toward shared PTSD-SUD vulnerability, the neurobehavioral traits that consistently distinguish comorbid from non-comorbid patients have not yet been identified. Machine learning (ML) is an analytic approach well suited to this task [33,34,35,36,37] and has been utilized by our research group to delineate a multivariate neurobehavioral profile of cocaine dependence [38]. Indeed, ML techniques have considerable advantages over traditional inferential case-control comparisons, including superior accuracy, the ability to handle numerous (and intercorrelating) predictors simultaneously, and the capacity to quantify the relative importance of individual predictors [39].
In the present study, we employed supervised ML to determine neurobehavioral phenotypic features that associated with post-traumatic stress symptoms within a SUD population. Candidate predictors were gathered using the standardized Phenotyping Assessment Battery (PhAB) constructed in collaboration with the National Institute of Drug Abuse (NIDA) [40] based on Research Domain Criteria principles [41, 42] for deployment in clinical trials of SUD pharmacotherapies [43]. The PhAB is a modular compendium of neurofunctional domain-based assessments [40]. Our study sought to identify comorbidity predictors from each of: the negative emotionality domain (e.g., low distress tolerance, anhedonia, anxiety, depression), metacognitive domain (e.g., distorted thought patterns and beliefs), reward domain (e.g., impulsivity, sensitivity to threat), interoceptive domain (e.g., willingness to tolerate aversive body sensations), and sleep domain (e.g., poor sleep quality). In addition, a subset of participants were administered a battery of neurocognitive measures to assess the predictive value of cognitive functioning, including working memory, attention, delay discounting, motoric impulsivity (i.e., response inhibition), and attentional bias towards threat. Consistent with our goal of identifying a neurobehavioral profile of comorbid post-traumatic stress symptoms in SUD, we constructed a series of classification and regression trees predicting post-traumatic stress symptom severity (i.e., score on the Post-traumatic Stress Disorder Checklist for DSM-5 (PCL-5)) as well as binary group identification (i.e., significant comorbid PTSD symptoms vs. no significant comorbid PTSD symptoms).
Materials and methods
The Virginia Commonwealth University Institutional Review Board approved the study and written informed consent was obtained from all participants. Data for the present investigation were collected in the course of a feasibility study for the NIDA PhAB [40].
Subjects
Participants were recruited for the parent study from an established participant registry and through local advertising [40]. Eligibility criteria were intentionally broad to recruit a “real-world” sample of individuals with SUD. To be eligible for participation, individuals were required to be: (a) between 18 and 70 years of age, (b) able to complete forms and interviews in English, and (c) meet DSM-5 diagnostic criteria for a current, primary SUD with the substance of choice being either cocaine, marijuana, or opioids. The MINI International Neuropsychiatric Interview Version 7.0.2 [44] was used to determine DSM-5 SUD diagnoses. The presence of any other condition or illness that, in the opinion of the Principal Investigator or study physician, would preclude safe and/or successful completion of the study was also a cause for exclusion. Participants could not meet DSM-5 criteria for a severe SUD related to other substances. However, due to the ubiquity of polysubstance use in the local population, individuals with mild to moderate comorbid SUD (i.e., alcohol use disorder) were not excluded. Confirmatory urine drug toxicology screens and alcohol breath samples were used in conjunction with self-report to assess recent substance use. Participants were excluded if they met any of the following criteria: (a) current psychosis, mania, suicidal or homicidal ideation, (b) history of seizures (excluding childhood febrile seizures), or (c) loss of consciousness from traumatic brain injury for more than 30 min. A total of 160 participants met inclusion/exclusion criteria and completed all self-report measures, while 102 participants completed additional neurocognitive measures. Participant characteristics are shown in Table 1, and there were no differences between groups with regard to age, sex, or drug of choice.
Assessments
Assessment-derived predictor variables
The NIDA PhAB “core” SUD-relevant domains evaluate cognition, reward, interoception, negative emotionality, metacognition, and sleep [40] (See Supplemental Information for fuller details on questionnaires and cognitive tasks). Nine self-report questionnaires and five neurocognitive tasks were employed here (Table 2). Each measure within the PhAB produced between one and eight scale/subscale scores. The battery of questionnaires/tasks was administered electronically via Redcap and Inquisit Lab software (Millisecond Software LLC, Seattle WA), and the average completion time for the battery was 180 min.
Assessment-derived outcome definitions
Because trauma exposure itself or even subthreshold PTSD symptoms may complicate SUD presentation [45] or treatment [46], the source PhAB protocol administered a dimensional probe of PTSD symptomatology using the PCL-5 [47]. The PCL-5 is a 20-item self-report questionnaire that assesses key DSM-5 criteria for PTSD. Items are rated on a 5-point Likert scale ranging from 0 to 4, and item ratings are summed to produce total scores ranging from 0 to 80 with higher scores reflecting greater symptom severity. Because most previous reports on PTSD-SUD comorbidity utilized clinical diagnostic criteria to obtain a binary (present/absent) diagnosis of PTSD, for harmonization with these studies, our analyses (described below) also featured a proxy dichotomization of (subclinical) vs. clinically-significant PTSD in addition to a continuous PTSD dimensional score. Based on previous findings [48], persons with scores ≥33 tend to have moderate to severe PTSD symptoms and probably meet DSM-5 criteria for PTSD, whereas persons with scores <33 tend to exhibit either a low level of PTSD symptoms or no PTSD symptoms. Thus, to demarcate participants with clinically significant comorbid PTSD, we employed this empirically derived cutoff score of 33. Below a score of 33, the participant was classified with “absent to mild” PTSD symptoms. At or above 33, the participant was classified with “clinically significant” PTSD symptoms. Thus, this binary outcome was employed for supervised ML. The effect sizes of association between each predictor variable and the target variable (no significant PTSD symptoms vs. significant PTSD symptoms) are shown in Supplementary Figs. 1 and 2.
Statistical analysis
We used R version 4.0.3 and Salford Predictive Modeler (SPM) version 8 to perform all statistical analyses, and 0.05 was set as the level of significance for all relevant hypothesis tests. The scores from the PhAB battery described above comprised an initial 72 explanatory (predictor) variables. Due to technical errors and/or participant failures to comply with task instructions, complete and valid neurocognitive task data were available from 102 participants. Our primary analysis centered on the full sample (n = 160) for whom complete questionnaire data were available. Subsequently, we performed follow-up analyses using data from the 102-participant subset (Table 1) wherein the neurocognitive performance variables (Table 2) were added to the model with the questionnaire-based metrics (Table 2). We were, therefore, able to infer whether neurocognitive metrics emerged as appreciable, i.e., top 20 predictors of PTSD symptomotology and whether their inclusion improved the prediction of PTSD symptom status.
The ML algorithms employed in the study were TreeNet and Classification and Regression Trees (CART). TreeNet [49] is an empirical variable selection procedure that can be used to efficiently reduce the number of explanatory variables in a predictive model. TreeNet algorithms [50, 51], also known as stochastic gradient boosting, offer several unique and useful features. These include (a) built-in estimation of prediction accuracy, (b) measures of feature importance, and (c) a measure of similarity between sample inputs. TreeNet improves upon classical decision trees such as CART [52] with retention of the most appealing properties of tree methods. The product of TreeNet is a ranked list of variable importance which is based on predictive models from an ensemble of weak learners in the form of classification trees. This methodology frequently outperforms Random Forest methods in prediction and variable selection [53]. “Boosting” is a method that seeks to convert weak learners into stronger ones. Here, a learner is an algorithm, and in the case of stochastic gradient boosting, that algorithm is CART; thus, boosting aims to improve on CART methods by creating multiple CART trees based on subsets of the data. The result is that TreeNet creates a final model in a gradual, additive, and sequential manner. All default options for SPM were used, and the algorithm was run with 80% of the data randomly selected for model training, and 20% of the data randomly selected for model testing. Based on the TreeNet results, the top 20 most important variables were used as predictor variables for the CART algorithm [54], which was used to predict the target outcome.
CART is a highly useful nonparametric method for building decision trees and predicting novel relationships between phenotypic predictors and biomedical outcomes because there is no requirement to select input variables based on theoretical importance. CART handles structurally complex datasets comprised of both categorical and continuous data with extreme robustness and limited vulnerability to outliers [49]. The three main components of CART involve creating a set of rules for splitting each node in a tree, deciding when a tree is fully grown, and assigning an outcome prediction to each terminal node of the tree [55]. Decision trees produce a clearly interpretable split at each node which is a binary response of some feature in the data set. The basic algorithm for building the decision tree seeks a data feature which maximizes the split between the classes contained in the parent node. CART is a recursive algorithm such that, once an appropriate split resulting in two child nodes is determined, the child nodes then become the new parent nodes, and the process is carried on down the branches of the tree. The CART tree is considered fully grown once a split cannot be identified that reduces model impurity. CART uses cross-validation techniques to determine the accuracy of the decision trees.
We first derived CART models based on a regression approach in which the target variable was the continuous score on the PCL-5. In a regression-based decision tree, the continuous target variable is split using recursive partitioning into a series of bins containing subsets of scores along the continuum. The algorithm does not attempt to predict the exact score on the target variable for each individual, but rather the bin into which that individual’s target score belongs. The average target score within each bin is therefore used as the predicted target score for all individuals who are presumed to belong in the respective bins. The number of bins and boundaries between bins are determined based on quantitative cross-validation aimed at reducing error. For instance, a small number of large bins typically diminishes precision, whereas a large number of small bins increases precision, but can sometimes result in overfitting and insufficient generalizability. We evaluated the accuracy of regression-based models using goodness-of-fit, or R2.
We subsequently derived classification-based CART models, which attempted to classify participants with subclinical levels PTSD symptoms vs. those with clinically-significant levels of PTSD symptoms. As such, the outcome categories were set using the PCL-5 trauma score split into a binary variable based on the empirically derived cutoff at 33. Model performance was examined using confusion matrices, sensitivity and specificity, the F1-score (harmonic mean of precision and recall), as well as area under the receiver operating characteristic (ROC) curve. Based on established benchmarks for the predictive accuracy of psychological tests within the literature [56], we deemed 70% accuracy to reflect acceptable performance, 80% accuracy to reflect good performance, and 90% accuracy to reflect excellent performance.
Results
Full sample analysis
Regression
The goodness-of-fit for the TreeNet training and testing datasets were 0.88 and 0.77, respectively, thus indicating an acceptable degree of explained variance. The 20 most important predictor variables identified by the TreeNet algorithm (Table 3), and the resulting CART algorithm produced a tree containing five primary parent (i.e., splitting) nodes and six child (i.e., terminal) nodes (Fig. 1). Variables used as splitting nodes, in order of importance, were the PROMIS Anxiety, Metacognitions Questionnaire—Uncontrollability/Danger, Buss–Perry Total Score, Buss–Perry Anger, and SUPPS Positive Urgency. Unfortunately, this model failed to demonstrate adequate goodness-of-fit, with R2 values in the training and testing datasets at 0.59 and 0.42, respectively.
White rectangles represent intermediate nodes (i.e., cohorts that could be further split); gray rectangles represent terminal nodes; PROMIS-A PROMIS Anxiety, MCQ-UD Metacognitions Questionnaire—Uncontrollability/Danger, SUPPS Short Urgency-Premeditation-Perseverance-Sensation Seeking-Positive Urgency Impulsive Behavior Scale, BPAS Buss–Perry Aggression Scale.
Classification
Model accuracy for the TreeNet training and testing datasets was 93.88 and 73.88%, respectively, indicating that the variable selection process performed well. Indeed, the examination of the areas under the ROC curves showed good to excellent predictive accuracy, with the value for the training data at 0.97 (sensitivity = 97.44%, specificity = 90.32%) and the value for the testing data at 0.84 (sensitivity = 80.00%, specificity = 66.67%). Confusion matrices and ROC curves for TreeNet analysis are presented in Supplementary Material (Supplementary Tables 1 and 2 and Supplementary Fig. 3). The 20 most important predictor variables identified by the TreeNet algorithm are shown in Table 4. Accuracies for the CART training and testing algorithms were 92.80% (sensitivity = 97.44%, specificity = 87.10%; F1 = 86.36%) and 76.11% (sensitivity = 80.00%, specificity = 72.22%; F1 = 69.57%), respectively (confusion matrices shown in Supplementary Tables 3 and 4). Areas under the ROC curves for both training and testing data (Supplementary Fig. 3) were 0.94 and 0.78, respectively. These results indicate slight model overfitting for the training data. The decision tree derived in CART is displayed in Fig. 2. Variables used as splitting nodes in the CART algorithm consisted of the following (in order of importance): PROMIS Anxiety, Metacognition Questionnaire—Uncontrollability/Danger, Multidimensional Assessment of Interoceptive Awareness—Not Distracting, PROMIS Depression, Metacognition Questionnaire—Total Score, and Multidimensional Assessment of Interoceptive Awareness—Attention Regulation.
White rectangles represent intermediate nodes (i.e., cohorts that could be further split); gray rectangles represent final nodes; Class 0 = PCL-5 < 33; Class 1 = PCL-5 ≥ 33; PROMIS-A PROMIS Anxiety, MCQ-UD Metacognitions Questionnaire—Uncontrollability/Danger, MAIA-ND Multidimensional Assessment of Interoceptive Awareness—Not Distracting, PROMIS-D PROMIS Depression, MAIA-AR Multidimensional Assessment of Interoceptive Awareness—Attention Regulation.
Neurocognitive sub-sample analysis
Regression
The goodness-of-fit for the TreeNet training and testing datasets were 0.69 and 0.64, respectively, thus indicating slightly sub-optimal model performance. The twenty most important predictor variables identified by the TreeNet algorithm are shown in Table 5. Results demonstrated that some neurocognitive task metrics supplanted self-report scales as predictors of trauma symptomatology, namely orienting and conflict effects in the Attention Network Task and happy and fear effects in the Emotional Go-NoGo Task (See Supplemental Material for detailed descriptions of neurocognitive task metrics). The subsequent CART algorithm produced a tree containing six primary splitting nodes and seven terminal nodes (Fig. 3). Variables used as splitting nodes, in order of importance, included PROMIS Anxiety, Metacognitions Questionnaire—Uncontrollability/Danger, Buss-Perry Hostility, Distress Tolerance Scale—Absorption, and Multidimensional Assessment of Interoceptive Awareness—Not Distracting. No neurocognitive assessments emerged as splitting nodes. This model demonstrated potentially adequate goodness-of-fit, with R2 values in the training and testing datasets at 0.72 and 0.64, respectively.
White rectangles represent intermediate nodes (i.e., cohorts that could be further split); gray rectangles represent terminal nodes; PROMIS-A PROMIS Anxiety, MCQ Metacognition Questionnaire, BPAS Buss–Perry Aggression Scale, DTS Distress Tolerance Scale, MAIA Multidimensional Assessment of Interoceptive Awareness.
Classification
Within the 102-participant sub-sample who completed all questionnaires and neurocognitive tasks, supervised ML analysis demonstrated good performance in predicting the presence vs. absence of clinically-significant ratings of PTSD symptoms. Model accuracy for TreeNet training and testing datasets were 86.21 and 73.33%, respectively (confusion matrices are shown in Supplementary Tables 5 and 6). Predictive performance, as evidenced by areas under the receiving operator characteristic curves, were 0.98 for the training data and 0.78 for the testing data (Supplementary Fig. 5). Although results indicated model overfitting for the training data, model accuracy was high enough to warrant further analysis. The 20 most important variables identified by the TreeNet algorithm, which included some neurocognitive assessments, are shown in Table 6. With the top 20 variables (self-report plus cognitive performance) entered as predictors, accuracies for the CART training and testing algorithms were 89.66% (Sensitivity = 100%; Specificity = 79.31%; F1 = 82.86%) and 86.67% (Sensitivity = 66.67%; Specificity = 75.00%; F1 = 41.39%) (confusion matrices shown in Supplementary Tables 7 and 8). Predictive accuracy, as reflected in an area under the curve, was 0.93 for the training data and 0.81 for the testing data (Supplementary Fig. 6). These results indicate good to excellent predictive accuracy. The decision tree derived in this CART analysis is displayed in Fig. 4. The only variables that emerged as splitting nodes in this CART algorithm consisted of the PROMIS Depression and the Fear Effect from the Emotional Go-NoGo Task (in order of importance).
Discussion
The multi-dimensional assessment battery and ML analytics reveal significant insights into neurobehavioral factors associated with comorbid PTSD symptoms among patients with SUD. For participants with only self-report assessment data, the prediction algorithm failed to accurately predict a dimensional PCL-5 score, but did accurately classify individuals as having significant comorbid PTSD symptoms (or not), based on symptoms of anxiety and depression as well as perseverative/intrusive thought patterns, low tolerance of aversive body sensations, and reduced ability to focus on body sensations. In participants who also completed neurocognitive assessments, continuous PTSD symptom severity and significant comorbid PTSD symptoms were predicted with impressive accuracy. Many of the same predictors that contributed to regression algorithms were also included in the classification algorithms, but the neurocognitive subsample classification algorithm uniquely included increased response bias toward by fear-related social information evaluated with the Emotional Go-NoGo Task.
Our results support previous research linking comorbid SUD and PTSD to increased negative affect and cognitive biases. Elevated negative affect is a key feature of PTSD and is strongly associated with SUD, thus our observation that comorbid presentations are associated with increased anxiety and depression is logical. The current study also suggests that PTSD symptomatology in SUD patients is associated with a general deficit in regulating aversive internal states, whether cognitive (e.g., worrisome thoughts), affective (e.g., distress), or interoceptive (e.g., aversive body sensations) in nature. Indeed, this notion is consistent with observations that anxiety sensitivity and sensitivity to punishment are shared vulnerability factors [21, 23, 32]. Anxiety sensitivity reflects the propensity to react negatively or fearfully to anxiety-related sensations, thoughts, emotions, or environmental stimuli and has been linked to PTSD [57], SUD [58], comorbid SUD-PTSD [59], and increased responsivity to the fear-dampening effects of alcohol [60, 61]. Further, there is indeed robust evidence of increased attentional bias toward social threat in anxiety-related disorders [62,63,64] and PTSD [65]. Exaggerated attentional bias toward negative information in anxiety disorders is thought to stem from an imbalance between an exaggerated bottom-up processing of threat or insufficient top-down regulation of threat response by executive control neurocircuitry [66]; this bias is thought to alter subsequent steps of cognitive processing, such as promoting rumination [62]. Consequently, several clinical trials have attempted to reduce attentional bias to threat directly, wherein the patient is trained over multiple trials to disengage from threat stimuli such as by directing gaze to alternative stimuli [67]. Similar bias modification techniques have also targeted attentional capture by substances of abuse [68], including more “gamified” approaches to reduce monotony [69]. These findings collectively lead to the hypothesis that top-down regulation of attention may be bolstered by repeated bias training for one modality (e.g., threat) which may aid in reducing attentional capture by another trigger (e.g., drug-related stimuli).
This study supports existing findings and offers new insights into the relationship between SUD and PTSD. We anticipated that the addition of neurocognitive task scores to predict PTSD symptomatology would improve model performance and that scores from neurocognitive tasks would be utilized as splitting nodes in decision trees. This expectation stems from fact that the constructs measured by neurocognitive assessments, such as response inhibition, tend to be strongly associated with many mental illnesses [70, 71]. However, neurocognitive tasks may also be affected by transient states of fatigue or mental distraction, or conversely prone to uncharacteristically high vigilance summoned under artificial conditions. Self-report questionnaires, with the generally longer time span and real-life related contexts specified in question items, are thought to be more stable than performance on “one-shot” computerized cognitive tasks [72]. As such, it is possible that the presence of PTSD symptoms among SUD patients does not impede cognitive abilities in all contexts, but rather is associated with a subtle shift in cognition that can be functionally important and noticeable over time but are difficult to detect via computerized tests. Another noteworthy feature of our regression-based models was the indication that some constructs may be associated specifically with the extremes of post-traumatic stress symptom severity. For example, aggression-related constructs (from the Buss-Perry Aggression Scale) were only featured in regression algorithms as differentiating very high scores from moderately high scores on the PCL-5. These data may indicate that aggression is only a notable feature in SUD patients with severe PTSD symptoms. Accordingly, perhaps intervention strategies aimed at persons with SUD and mild-to-moderate comorbid PTSD should not be aimed at anger management or social functioning, whereas these may be important components of an intervention for severe comorbid PTSD.
In spite of these insights, the results should be considered in light of some limitations. First, because the study did not include a third group of participants with heightened PTSD symptoms alone (no SUD), we cannot conclude that the traits/risk factors identified truly differentiate comorbid SUD-PTSD from each disorder alone. The predictors highlighted in our results may be reflective of the increased overall disorder burden of comorbid conditions rather than SUD or PTSD alone [73]. Therefore, future studies will benefit from approaches that identify which psychopathological mechanisms differentiate comorbid SUD-PTSD from each disorder on its own. We also did not collect data on comorbid psychiatric diagnoses, which could have influenced study results and could be highly relevant for personalized treatment approaches. Additionally, while the inclusion of persons across multiple types of SUD can enable our predictive variables to serve as potential risk markers that cut across different types of SUDs, our results could presumably vary depending upon the primary substance of abuse. Indeed, individuals with more significant PTSD symptoms related to the hyperarousal symptom cluster were more likely to abuse alcohol, while those with more of avoidance and re-experiencing symptoms were more likely to abuse cocaine [74]. Future studies are therefore needed to compare data across multiple samples of patients who abuse different substances, including substance use disorders not examined in the current study (e.g., alcohol, methamphetamine) and polysubstance use disorders. Relatedly, heterogeneity among symptoms of PTSD may have impacted our findings, such that important subtypes of post-traumatic stress may be associated with unique phenotypic signatures. Multiple subtypes of PTSD have been proposed, including externalizing vs. internalizing subtypes as well as a dissociating subtype [75]. It is therefore possible that participants with comorbid PTSD symptoms were affected by distinct clusters of symptoms, which do not translate well into a linear models of disorder severity, and thus are not easily predicted by our ML approaches. While our predictive algorithms achieved relatively high accuracy, there is a clear need for future studies to deliver greater predictive accuracy by either increasing sample size or including additional biological predictors, particularly neurophysiological variables such as peripheral biomarkers, neuroimaging metrics, and genetic information.
Nevertheless, the performance of our models support a combination of self-report traits and neurocognitive abilities to be involved in the underlying etiology of comorbid SUD and PTSD. Our results from a machine-learning-based approach largely align with previous work, yet also provide a efficient integration of works investigating the role of a single risk factor or etiological feature. Future studies are planned to replicate these results and determine whether targeted treatments aimed at the identified psychopathological processes are associated with improved care for this high-need population.
References
Kessler RC, Nelson CB, McGonagle KA, Edlund MJ, Frank RG, et al. The epidemiology of co-occurring addictive and mental disorders: implications for prevention and service utilization. Am J Orthopsychiatry. 1996;66:17–31.
Swendsen J, Conway KP, Degenhart L, Glantz M, Jin R, Merikangas KR, et al. Mental disorders as risk factors for substance use, abuse and dependence: results from the 10-year follow-up of the National Comorbidity Survey. Addiction. 2010;105:1117–28.
Agosti V, Nunes E, Levin F. Rates of psychiatric comorbidity among U.S. residents with lifetime cannabis dependence. Am J Drug Alcohol Abus. 2002;28:643–52.
Afifi TO, Henriksen CA, Asmundson GJ, Sareen J. Childhood maltreatment and substance use disorders among men and women in a nationally representative sample. Can J Psychiatry. 2012;57:677–86.
Jacobsen LK, Southwick SM, Kosten TR. Substance use disorders in patients with posttraumatic stress disorder: a review of the literature. Am J Psychiatry. 2001;158:1184–90.
Brady KT, Back SE, Coffey SF. Substance abuse and posttraumatic stress disorder. Curr Dir Psychol Sci. 2004;13:206–9.
Schafer I, Najavits LM. Clinical challenges in the treatment of patients with posttraumatic stress disorder and substance abuse. Curr Opin Psychiatry. 2007;20:614–8.
Hien D, Cohen L, Campbell A. Is traumatic stress a vulnerability factor for women with substance use disorders? Clin Psychol Rev. 2005;25:813–23.
Zweben JE, Clark HW, Smith DE. Traumatic experiences and substance abuse: mapping the territory. J Psychoact Drugs. 1994;26:327–44.
Kessler RC. The epidemiology of dual diagnosis. Biol Psychiatry. 2004;56:730–7.
Leeies M, Pagura J, Sareen J, Bolton JM. The use of alcohol and drugs to self-medicate symptoms of posttraumatic stress disorder. Depress Anxiety. 2010;27:731–6.
Khantzian EJ. The self-medication hypothesis of addictive disorders: focus on heroin and cocaine dependence. Am J Psychiatry. 1985;142:1259–64.
McCauley JL, Killeen T, Gros DF, Brady KT, Back SE. Posttraumatic stress disorder and co-occurring substance use disorders: advances in assessment and treatment. Clin Psychol. 2012;19:283.
Coffey SF, Read JP, Norberg MM. Posttraumatic stress disorder and substance use disorder: neuroimaging, neuroendocrine, and psychophysiological findings. In: Stewart SH, Conrod PJ, editors. Anxiety and substance use disorders. New York, NY: Springer; 2008. Ch. 3.
Hien DA, Fitzpatrick S, Saavedra LM, Ebrahimi CT, Norman SB, Tripp J, et al. What’s in a name? A data-driven method to identify optimal psychotherapy classifications to advance treatment research on co-occurring PTSD and substance use disorders. Eur J Psychotraumatol. 2022;13:2001191.
Lenz AS, Henesy R, Callender K. Effectiveness of seeking safety for co‐occurring posttraumatic stress disorder and substance use. J Couns Dev. 2016;94:51–61.
Najavits LM, Hien D. Helping vulnerable populations: A comprehensive review of the treatment outcome literature on substance use disorder and PTSD. J Clin Psychol. 2013;69:433–79.
Simpson TL, Lehavot K, Petrakis IL. No wrong doors: findings from a critical review of behavioral randomized clinical trials for individuals with co-occurring alcohol/drug problems and posttraumatic stress disorder. Alcohol Clin Exp Res. 2017;41:681–702.
van Dam D, Vedel E, Ehring T, Emmelkamp PM. Psychological treatments for concurrent posttraumatic stress disorder and substance use disorder: a systematic review. Clin Psychol Rev. 2012;32:202–14.
Roberts NP, Roberts PA, Jones N, Bisson JI. Psychological interventions for post-traumatic stress disorder and comorbid substance use disorder: a systematic review and meta-analysis. Clin Psychol Rev. 2015;38:25–38.
Vujanovic AA, Farris SG, Bartlett BA, Lyons RC, Haller M, Colvonen PJ, et al. Anxiety sensitivity in the association between posttraumatic stress and substance use disorders: a systematic review. Clin Psychol Rev. 2018;62:37–55.
Maria-Rios CE, Morrow JD. Mechanisms of shared vulnerability to post-traumatic stress disorder and substance use disorders. Front Behav Neurosci. 2020;14:6.
Vujanovic AA, Farris SG, Bartlett BA, Lyons RC, Haller M, Colvonen PJ, et al. Reward functioning in posttraumatic stress and substance use disorders. Curr Opin Psychol. 2017;14:49–55.
Grubbs JB, Chapman H. Predicting gambling situations: the roles of impulsivity, substance use, and post-traumatic stress. Subst Abus. 2019;13:1178221819852641.
James LM, Strom TQ, Leskela J. Risk-taking behaviors and impulsivity among veterans with and without PTSD and mild TBI. Mil Med. 2014;179:357–63.
Falconer E, Bryant R, Felmingham KL, Kemp AH, Gordon E, Peduto A, et al. The neural networks of inhibitory control in posttraumatic stress disorder. J Psychiatry Neurosci. 2008;33:413–22.
Walker J, Bountress KE, Calhoun CD, Metzger IW, Adams Z, Amstadter A, et al. Impulsivity and comorbid PTSD and binge drinking. J Dual Diagn. 2018;14:89–95.
Weiss NH, Tull MT, Anestis MD, Gratz KL. The relative and unique contributions of emotion dysregulation and impulsivity to posttraumatic stress disorder among substance dependent inpatients. Drug Alcohol Depend. 2013;128:45–51.
Jenkins MA, Langlais PJ, Delis D, Cohen RA. Attentional dysfunction associated with posttraumatic stress disorder among rape survivors. Clin Neuropsychol. 2000;14:7–12.
Koenen KC, Driver KL, Oscar-Berman M, Wolfe J, Folsom S, Huang MT, et al. Measures of prefrontal system dysfunction in posttraumatic stress disorder. Brain Cogn. 2001;45:64–78.
Vasterling JJ, Duke LM, Brailey K, Constans JI, Allain AN Jr, Sutker PB. Attention, learning, and memory performances and intellectual resources in Vietnam veterans: PTSD and no disorder comparisons. Neuropsychology. 2002;16:5–14.
Sokhadze E, Singh S, Stewart S, Hollifield M, El-Baz A, Tasman A. Attentional bias to drug- and stress-related pictorial cues in cocaine addiction comorbid with PTSD. J Neurother. 2008;12:205–25.
Lueken U, Straube B, Yang Y, Hahn T, Beesdo-Baum K, Wittchen HU, et al. Separating depressive comorbidity from panic disorder: a combined functional magnetic resonance imaging and machine learning approach. J Affect Disord. 2015;184:182–92.
Jauhar S, Krishnadas MM, Nour D, Cunningham-Owens D, Johnstone EC, Lawrie SM. Is there a symptomatic distinction between the affective psychoses and schizophrenia? A machine learning approach. Schizophr Res. 2018;202:241–7.
Redlich R, Almeida JR, Grotegard D, Opel N, Kugel H, Heindel W, et al. Brain morphometric biomarkers distinguishing unipolar and bipolar depression. A voxel-based morphometry-pattern classification approach. JAMA Psychiatry. 2014;71:1222–30.
Duda M, Ma R, Haber N, Wall DP. Use of machine learning for behavioral distinction of autism and ADHD. Transl Psychiatry. 2016;6:e732.
Duda M, Haber N, Daniels J, Wall DP. Crowdsourced validation of a machine-learning classification system for autism and ADHD. Transl Psychiatry. 2017;7:e1133.
Ahn WY, Ramesh D, Moeller FG, Vassileva J. Utility of machine-learning approaches to identify behavioral markers for substance use disorders: impulsivity dimensions as predictors of current cocaine dependence. Front Psychiatry. 2016;7:34.
Tai AM, Albuquerque A, Carmona NE, Subramanieapullai M, Cha DS, Sheko M, et al. Machine learning and big data: implications for disease modeling and therapeutic discovery in psychiatry. Artif Intell Med. 2019;99:101704.
Keyser-Marcus LA, Ramey T, Bjork J, Adams A, Moeller FG. Development and feasibility study of an addiction-focused phenotyping assessment battery. Am J Addict. 2021;30:398–405.
Cuthbert BN, Insel TR. Toward the future of psychiatric diagnosis: the seven pillars of RDoC. BMC Med. 2013;11:126.
Brooks SJ, Lochner C, Shoptaw S, Stein DJ. Using the research domain criteria (RDoC) to conceptualize impulsivity and compulsivity in relation to addiction. Prog Brain Res. 2017;235:177–218.
Watts AL, Latzmann RD, Boness CL, Kotov R, Keyser-Marcus L, DeYoung CG, et al. New approaches to deep phenotyping in addictions. Psychol Addict Behav. 2023;37:361–75.
Sheehan DV, Lecrubier Y, Sheehan KH, Amorim P, Janavs J, Weiller E, et al. The Mini-International Neuropsychiatric Interview (M.I.N.I.): the development and validation of a structured diagnostic psychiatric interview for DSM-IV and ICD-10. J Clin Psychiatry. 1998;59:22–33.
Tate SR, Norman SB, McQuaid JR, Brown SA. Health problems of substance-dependent veterans with and those without trauma history. J Subst Abus Treat. 2007;33:25–32.
Schneider R, Cronkite R, Timko C. Lifetime physical and sexual abuse and substance use treatment outcomes in men. J Subst Abus Treat. 2008;35:353–61.
Blevins CA, Weathers FW, Davis MT, Witte TK, Domino JL. The posttraumatic stress disorder checklist for DSM‐5 (PCL‐5): development and initial psychometric evaluation. J Trauma Stress. 2015;28:489–98.
Bovin MJ, Marx BP, Weathers FW, Gallagher MW, Rodriguez P, Schnurr PP, et al. Psychometric properties of the PTSD checklist for diagnostic and statistical manual of mental disorders–fifth edition (PCL-5) in veterans. Psychol Assess. 2016;28:1379.
Arbet J, Brokamp C, Meinzen-Derr J, Trinkley KE, Spratt HM. Lessons and tips for designing a machine learning study using EHR data. J Clin Transl Sci. 2020;5:e21.
Friedman JH. Greedy function approximation: a gradient boosting machine. Ann Stat. 2001;29:1189–232.
Friedman JH. Stochastic gradient boosting. Comput Stat Data Anal. 2002;38:367–78.
Breiman L. Classification and regression trees. New York: Routledge, 2017.
Hastie T, Friedman J, Tibshirani R. Boosting and Additive Trees. In: Hastie T, Friedman J, Tibshirani R, editors. The elements of statistical learning: data mining, inference, and prediction. New York: Springer New York; 2009. Ch. 10.
Spratt H, Ju H, Brasier AR. A structured approach to predictive modeling of a two-class problem using multidimensional data sets. Methods. 2013;61:73–85.
Steinberg D, Colla P, CART: tree-structured non-parametric data analysis. San Diego, CA: Salford Systems; 1995.
Youngstrom EA. A primer on receiver operating characteristic analysis and diagnostic efficiency statistics for pediatric psychology: we are ready to ROC. J Pediatr Psychol. 2014;39:204–21.
Taylor S, Koch WJ, McNally RJ. How does anxiety sensitivity vary across the anxiety disorders? J Anxiety Disord. 1992;6:249–59.
Kushner MG, Thuras P, Abrams K, Brekke M, Stritar L. Anxiety mediates the association between anxiety sensitivity and coping-related drinking motives in alcoholism treatment patients. Addict Behav. 2001;26:869–85.
Bonin MF, Norton GR, Asmundson GJ, Dicurzio S, Pidlubney S. Drinking away the hurt: the nature and prevalence of PTSD in substance abuse patients attending a community-based treatment program. J Behav Ther Exp Psychiatry. 2000;31:55–66.
MacDonald AB, Baker JM, Stewart SH, Skinner M. Effects of alcohol on the response to hyperventilation of participants high and low in anxiety sensitivity. Alcohol Clin Exp Res. 2000;24:1656–65.
Stewart SH, Pihl RO. Effects of alcohol administration on psychophysiological and subjective-emotional responses to aversive stimulation in anxiety-sensitive women. Psychol Addict Behav. 1994;8:29.
Mogg K, Bradley BP. Anxiety and threat-related attention: cognitive-motivational framework and treatment. Trends Cogn Sci. 2018;22:225–40.
Armstrong T, Olatunji BO. Eye tracking of attention in the affective disorders: a meta-analytic review and synthesis. Clin Psychol Rev. 2012;32:704–23.
Shechner T, Jarcho JM, Britton JC, Leibenluft E, Pine DS, Nelson EE. Attention bias of anxious youth during extended exposure of emotional face pairs: an eye-tracking study. Depress Anxiety. 2013;30:14–21.
Powers A, Fani N, Murphy L, Briscione M, Bradley B, Tone EB, et al. Attention bias toward threatening faces in women with PTSD: eye tracking correlates by symptom cluster. Eur J Psychotraumatol. 2019;10:1568133.
Schoorl M, Putman P, Van Der Werff S, Van Der Does AW. Attentional bias and attentional control in posttraumatic stress disorder. J Anxiety Disord. 2014;28:203–10.
Hang Y, Zhang G, Wang C, Zhang N, Liu G. Attention bias modification for anxiety disorders in children and adolescents: a systematic review and meta-analysis. Psychiatry Res. 2021;300:113896.
Zhang MW, Ying JB, Song G, Fung DS, Smith HE. Recent advances in attention bias modification for substance addictions. Int J Environ Res Public Health. 2018;15:676.
Zhang MW, Heng S, Song G, Fung DSS, Smith H. Codesigning a mobile gamified attention bias modification intervention: research protocol. BMJ Open. 2019;9:e031897.
Nelson LD, Strickland C, Krieger RF, Arbisi PA, Patrick CJ. Neurobehavioral traits as transdiagnostic predictors of clinical problems. Assessment. 2016;23:75–85.
McTeague LM, Goodkind MS, Etkin A. Transdiagnostic impairment of cognitive control in mental illness. J Psychiatr Res. 2016;83:37–46.
Cyders MA, Coskunpinar A. Measurement of constructs using self-report and behavioral lab tasks: is there overlap in nomothetic span and construct representation for impulsivity? Clin Psychol Rev. 2011;31:965–82.
Caspi A, Moffitt TE. All for one and one for all: mental disorders in one dimension. Am J Psychiatry. 2018;175:831–44.
Saladin ME, Brady KT, Dansky BS, Kilpatrick DG. Understanding comorbidity between PTSD and substance use disorders: two preliminary investigations. Addict Behav. 1995;20:643–55.
Dalenberg CJ, Glaser D, Alhassoon OM. Statistical support for subtypes in posttraumatic stress disorder: the how and why of subtype analysis. Depress Anxiety. 2012;29:671–8.
Acknowledgements
The authors thank Edward Zuniga, Aysha Sami, Mandy Adams, and Erika Lopex for data collection. The authors also thank Dr. Per Carlbring for affording access to the Umea face images, and Dr. Nim Tottenham for advice in task configuration.
Funding
This study was supported by a grant (FGM) from the National Institute of Drug Abuse (U54DA038999) and an Institutional Career Development Award core grant (DCH; 5KL2TR001441). TR was substantially involved in U54DA038999, consistent with her role as Scientific Officer. Any opinion, findings, and conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of the National Institute on Drug Abuse.
Author information
Authors and Affiliations
Contributions
DCH contributed to the analysis, interpretation of results, and writing of the manuscript. All other authors edited the manuscript draft and contributed content, in addition to other unique contributions. HMS analyzed the data. LK-M, TR, and FGM designed the core study recruitment plan and assessment battery. JMB and KAC contributed to the analytic strategy and interpretation of results. GNN contributed to the interpretation of results.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Houghton, D.C., Spratt, H.M., Keyser-Marcus, L. et al. Behavioral and neurocognitive factors distinguishing post-traumatic stress comorbidity in substance use disorders. Transl Psychiatry 13, 296 (2023). https://doi.org/10.1038/s41398-023-02591-3
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41398-023-02591-3
This article is cited by
-
Post-traumatic stress comorbidity in substance use disorder: machine learning analyses of phenotypic drivers
Neuropsychopharmacology (2024)