Next Article in Journal
Illicit Substance Use and the COVID-19 Pandemic in the United States: A Scoping Review and Characterization of Research Evidence in Unprecedented Times
Previous Article in Journal
Nicotine Exerts Cytotoxic Effects in a Panel of Healthy Cell Lines and Strong Irritating Potential on Blood Vessels
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Identification of Subphenotypes and Associations with Health Outcomes in Patients with Opioid-Related Emergency Department Encounters Using Latent Class Analysis

1
Division of Medical Toxicology, Department of Emergency Medicine, Cook County Health, Chicago, IL 60612, USA
2
Department of Emergency Medicine, Rush Medical College, Rush University, Chicago, IL 60612, USA
3
Addiction Data Science Laboratory, Department of Psychiatry & Behavioral Science, Rush University Medical Center, Chicago, IL 60612, USA
4
Department of Psychology, Olivet Nazarene University, Bourbonnais, IL 60914, USA
5
Doctor of Medicine Program, Rush Medical College, Rush University, Chicago, IL 60612, USA
6
School of Medicine and Public Health, University of Wisconsin, Madison, WI 53715, USA
7
Department of Medicine, University of Wisconsin-Madison, Madison, WI 53715, USA
8
Institute for Juvenile Research, Department of Psychiatry, University of Illinois Chicago, Chicago, IL 60612, USA
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2022, 19(14), 8882; https://doi.org/10.3390/ijerph19148882
Submission received: 17 June 2022 / Revised: 17 July 2022 / Accepted: 19 July 2022 / Published: 21 July 2022

Abstract

:
The emergency department (ED) is a critical setting for the treatment of patients with opioid misuse. Detecting relevant clinical profiles allows for tailored treatment approaches. We sought to identify and characterize subphenotypes of ED patients with opioid-related encounters. A latent class analysis was conducted using 14,057,302 opioid-related encounters from 2016 through 2017 using the National Emergency Department Sample (NEDS), the largest all-payer ED database in the United States. The optimal model was determined by face validity and information criteria-based metrics. A three-step approach assessed class structure, assigned individuals to classes, and examined characteristics between classes. Class associations were determined for hospitalization, in-hospital death, and ED charges. The final five-class model consisted of the following subphenotypes: Chronic pain (class 1); Alcohol use (class 2); Depression and pain (class 3); Psychosis, liver disease, and polysubstance use (class 4); and Pregnancy (class 5). Using class 1 as the reference, the greatest odds for hospitalization occurred in classes 3 and 4 (Ors 5.24 and 5.33, p < 0.001) and for in-hospital death in class 4 (OR 3.44, p < 0.001). Median ED charges ranged from USD 2177 (class 1) to USD 2881 (class 4). These subphenotypes provide a basis for examining patient-tailored approaches for this patient population.

1. Introduction

Drug overdose is a leading cause of death in the United States, with the majority involving an opioid [1]. There were approximately 50,000 opioid overdose deaths reported in 2019, and the recent COVID-19 pandemic has seen increases in opioid-related deaths [2,3]. Patients with opioid misuse disproportionately utilize emergency medical services, with almost 1.5 million opioid-related emergency department (ED) encounters annually [4,5,6]. ED patients have been found to have more severe substance use patterns compared with patients seen in the primary care setting [7]. The one-year mortality rate for patients following an opioid-related ED visit may be greater than 5% [8,9]. Therefore, the ED represents a critical point of contact between patients with opioid misuse and the health care system as ED encounters are opportunities for interventions to prevent future morbidity and mortality [10]. Several interventions initiated in the ED setting have been effective in improving outpatient treatment engagement and preventing overdose, including buprenorphine initiation, take-home naloxone, team-based care coordination, and peer recovery support services [11,12,13,14,15,16]. While the utilization and expansion of these services have improved the care of ED patients with opioid misuse, none have been shown to be uniformly effective.
Opioid misuse is characterized by a heterogeneous pattern influenced by several factors including type of opioid used, severity of illness, exposure route, comorbid conditions, and social determinants of health [17]. In short, opioid misuse is not a one-size-fits-all condition. It is therefore unsurprising that ED-based interventions have not been uniformly effective for all patients. Current evidence suggests that treatment responses and clinical outcomes differ between specific subphenotypes of patients with opioid misuse [18,19]. Such subphenotypes are defined by varying risks for poor outcome or shared underlying biologic factors. Ideally, management approaches for patients with opioid-related encounters would be tailored to the needs of individual patients. To personalize treatment approaches within the ED setting, clinically relevant subphenotypes of patients with opioid misuse must first be discovered and described. Distinct clusters of patients with opioid misuse, defined by shared characteristics, have been identified in the general population and clinical settings outside of the ED [20,21,22]. There is a need to determine if such subphenotypes exist within the ED setting, and if these subphenotypes differ in patient characteristics and outcomes. These subphenotypes would inform treatment considerations, resource allocation, public policy, and future research.
To identify and describe subphenotypes of patients with opioid-related ED encounters, we applied latent class analysis (LCA) to a nationally representative sample of ED encounters. LCA utilizes multivariate categorical data to surface informative patterns in a dataset to cluster subsets of patient encounters based on shared characteristics [23,24]. LCA has been used to detect clinically informative subphenotypes of patients for health conditions including alcohol misuse, e-cigarette use, and respiratory failure, among others [25,26,27]. We hypothesize that by using LCA we will detect distinct and clinically relevant subphenotypes of ED patients with opioid-related encounters; these subphenotypes signal key differences in patient characteristics with implications for patient outcomes as well as potential treatment and intervention pathways.

2. Materials and Methods

2.1. Data Source

Clinical data were obtained from the National Emergency Department Sample (NEDS), the largest all-payer emergency department database in the United States. NEDS is a nationally representative sample developed and maintained by the Healthcare Cost and Utilization Project (HCUP) [28]. NEDS contains systematically stratified discharge data from approximately 1000 contributing hospitals from all regions of the United States, using 20% sampling across strata (geographic region, urban or rural location, teaching status, ownership, and trauma-level designation) from all emergency department visits nationwide. The database contains a combination of structured billing and demographic variables.

2.2. Study Population

The sample population included patients over 12 years of age with an opioid-related International Classification of Diseases, Tenth Revision (ICD-10) diagnosis code (Table S1) as a primary or secondary diagnosis excluding codes negating opioid use, indicating remote use, or indicating disease in remission. Opioid-related ICD-10 codes were adopted from the Healthcare Cost and Utilization Project and represent a commonly used definition for health services research [29]. ICD-10 codes reflect final billing diagnosis codes used for claims by payers. For the purposes of the study, the sample was not stratified by encounter reason, such as opioid overdose or withdrawal, as the goal of the investigation was to identify subphenotypes by patient characteristics rather than encounter-specific characteristics, which are more likely to change over time. Adolescent encounters were included as substance misuse because such use often begins in the adolescent period, and many characteristics of substance use do not greatly differ between adolescents and adults [30]. NEDS does not utilize patient identifiers, so multiple encounters by the same patient were included as independent observations. We merged NEDS files from 2016 and 2017 for analysis. This time frame was chosen due to the increasing contribution of synthetic opioids to the opioid epidemic in recent years and for data consistency as 2016 represents the first entire year that included ICD-10 diagnosis codes, covariates, and sample weights. The year 2017 is the most recent year with available data at the time of analysis.

2.3. Variables and Covariates

Candidate variables within NEDS included those with previously established associations with opioid misuse, sufficient interpretability for providers, and conceptual association with overdose risk [20,21]. These class-defining variables for the LCA models consisted of groups of diagnosis codes for the following disease processes: (1) chronic pain; (2) alcohol use disorders; (3) psychoses; (4) depression; (5) liver disease; (6) pregnancy; (7) cocaine use; and (8) amphetamine use (Table S2). These variables were chosen a priori based on their previous use in effectively clustering patients with substance misuse using LCA in settings outside the ED [20,21,31]. Demographic variables utilized as covariates included age, sex, payer, hospital urban–rural designation, and median household income for the patient’s ZIP code (stratified by quartile); these variables were selected because of their known associations with substance use [32,33]. The six-category urban–rural classification scheme was developed by the National Center for Health Statistics to differentiate between central and fringe metropolitan areas with smaller areas subdivided by population. Race, ethnicity, and level of education are not variables collected by NEDS and were therefore not utilized as covariates.

2.4. Statistical Analyses

A series of nested LCA models with K + 1 classes were iteratively fit up to a 10-class model utilizing the dichotomous variables derived from ICD-10 codes and weights accounting for the complex sampling design of NEDS. Models were fit by maximum likelihood estimation using the expectation–maximization algorithm. The optimal number of latent classes was determined using clinical interpretability in addition to the exploration of information criteria-based metrics and model entropy. The goal was to identify the minimum number of classes that best described and fit the data. Among models with similar performance, the most parsimonious model was chosen. The metrics for model selection included the adjusted Bayesian information criterion (aBIC), Akaike information criterion (AIC), and chi-squared goodness of fit. Initial LCA was explored on all patients across three years (2015–2017), although because covariate information and sample weights were not available for 2015, weighted analyses were conducted only with the 2016 and 2017 data after verifying the robustness of class structure within this reduced sample.
Class sizes were determined by class membership using the highest posterior probability for each patient encounter to identify class sizes that were either a too-large or too-small proportion of the population to be clinically informative. Previous simulations have shown that classes accounting for less than 5% of the total cohort are less likely to be informative [34]. Face validity of the subphenotypes was determined by clinical interpretability on manual review by researchers with expertise in emergency care, substance use research, and/or LCA methodology (NC, MA, CM, and NK).
Covariates were explored using the three-step approach due to recent concerns regarding utilizing covariates in the initial exploration of class structure increasing the potential for bias [35,36]. The three-step approach initially assesses class structure independent of covariates, assigns individuals to classes based on posterior probability of class membership, and then explores the covariates’ relationships with class membership. We also explored the outcomes of inpatient hospitalization, in-hospital death, and ED charges (USD) across the latent classes. To identify differences in demographics and outcomes between subphenotypes, descriptive statistics and odds ratios were calculated for the demographic variables and the patient outcomes. Odds ratios were determined by multinomial logistic regression predicting class membership by covariate while correcting for classification error. They represent raw associations between covariate and class and were not adjusted for other covariates used in the latent class analysis. Due to the large sample size, traditional significance tests and p-values were uninformative (all p-values < 0.001) and standard errors were generally < 0.001. Emphasis was thus placed on the descriptive statistics and odds ratios for class membership or class makeup based on covariates. For interpretability, the outcome of ED charges was determined as median raw charges for each class, and for between-class comparisons, ED charges were transformed to a natural logarithmic scale to account for the strong positive skew and non-normal distribution. Analyses were performed using Stata Statistical Software: Release 17 (College Station, TX, USA: StataCorp LLC. StataCorp. 2019) and Mplus (Los Angeles, CA: Muthén & Muthén). Missingness occurred mostly in the ED charge field, a limitation noted by HCUP for the 2016 dataset, and was handled by listwise deletion [28]. The protocol was reviewed and deemed exempt by the institutional review board of Rush University Medical Center. This study conforms, where applicable, to the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement guidelines (Table S3) [37].

3. Results

3.1. Model Selection

A total of 14,057,302 ED encounters met the inclusion criteria for analysis. The information criteria and smallest class sizes from LCA models with classes from 1 to 10 are shown in Table 1. All models with six or more classes contained multiple classes representing less than 5% of the study population. Information criteria generally indicated preference for larger numbers of classes (Figure S1); however, the five-class model was the most parsimonious with optimal clinical interpretability and class sizes. After consideration of model fit indices, relative class sizes, and interpretability, the five-class model was chosen as the final model (Figure 1). Table 2 displays the class estimate probabilities for the inclusion of each class-defining variable for each class in the final model.
While the seven-class LCA model had more favorable metrics, including model entropy, these gains were minimal compared with the increased complexity of the model. Manual review of this model for face validity revealed almost-identical classes as the five-class LCA model with the addition of two small classes (4.1% and 1.9% of the study population), psychoses and polysubstance use, which were subsumed by other classes in the final five-class model. The mean posterior latent class probabilities for class membership exceeded 0.94 for all classes in the five-class model. The five classes are summarized as the following patient subphenotypes: Chronic pain (class 1); Alcohol use (class 2); Depression with chronic pain (class 3, hereafter “depression with pain”); Psychosis, liver disease, and polysubstance use (class 4); and Pregnancy (class 5).

3.2. Subphenotype Characteristics

Descriptive statistics of the patient characteristics for each subphenotype are shown in Table 3. Females constituted the majority of patients in classes 1 (chronic pain), 3 (depression and pain), and 5 (pregnancy). Class 5 (pregnancy) contained the highest proportion of patients located in central metropolitan areas, while class 3 (depression and pain) had the highest proportion located in areas outside of central metropolitan locations. Class 2 (alcohol use) contained the highest proportion of patients designated as self-pay. For between-group comparisons, odds ratios are provided in the multimedia appendices (Table S4).

3.3. Outcomes and Cost Analysis

The results for patient outcomes of hospitalization, in-hospital death, and ED charge are presented in Table 4. Class 1 (chronic pain) was chosen as the reference category as it represented the largest subphenotype. Compared with the reference class, the odds for hospitalization, in-hospital death, and higher ED charges were higher for all other classes. The only exception to this general trend was the pregnancy class (class 5) with regard to the outcome of death, as a very small proportion of patients in the pregnancy class died (OR = 0.00, 95% CI 0.00-0.00). Classes 2 (alcohol use), 3 (depression and pain), and 4 (psychosis, liver disease, and polysubstance use) had the greatest odds of hospital admission, while class 4 carried the greatest odds of in-hospital death. The median ED charge per encounter ranged from USD 2177 to USD 2881, with the lowest and highest charges belonging to class 1 and class 4, respectively.

4. Discussion

The ED is a critical setting for the initiation of treatment and harm-reduction interventions for patients with opioid misuse. Identifying subphenotypes of patients with opioid-related encounters is a necessary step towards personalized interventions in this setting. Our analysis, utilizing a nationally representative population of ED patients, provides analytic support for five clinically relevant subphenotypes of opioid-related encounters. These subphenotypes best represented the population of ED patients with opioid-related encounters based on a combination of information criteria-based metrics, class sizes, and clinical interpretability. The high mean posterior latent class probabilities for membership in all classes of the model indicated a strong fit of these classes with the dataset. The subphenotypes detected provide a critical framework for conceptualizing the heterogeneity among ED patients with opioid misuse and carry implications for the development of personalized treatment pathways.
Each subphenotype demonstrated distinct demographic profiles and variations in clinical outcomes. The largest subphenotype consisted of those with an opioid-related diagnosis and chronic pain. This subphenotype had a majority female membership and among the highest proportion of all classes for private insurance (excluding subphenotype 5, consisting of pregnant patients). The outcomes of hospital admission, in-hospital death, and total ED charges were among the lowest for this subphenotype. Addressing opioid misuse within this subphenotype can be difficult as patients may rely on opioids for the relief of chronic pain, and opioids may be the cause of both benefit and harm within individual patients. Additionally, chronic pain is a broad diagnosis affecting patients across demographic groups with multiple, often overlapping causes, which further complicates treatment decisions. Previous estimates of opioid misuse in patients with chronic pain have ranged from 21% to 29% [38]. These estimates, along with the current findings, suggest a role for opioid misuse screening in patients with chronic pain syndromes. The recognition of chronic pain diagnoses as a class-defining variable in the largest subphenotype of opioid-related ED encounters necessitates further exploration of how chronic pain, opioid use, and emergency department utilization are related.
Two opioid-related subphenotypes were defined primarily by the co-occurring diagnosis of psychiatric disorders (depression and psychoses). The associations of opioid misuse with both depression and psychosis have been noted previously, and both conditions are associated with an increased risk for opioid overdose [39]. These classes demonstrated the highest proportion of public insurance and the greatest odds for hospital admission. These classes, in aggregate, account for over a quarter of all opioid-related ED encounters, and their presence underscores the potential opportunities in addressing uncontrolled psychiatric disorders in addition to the risks from opioid misuse during an ED encounter.
The subphenotype defined by psychosis also had the highest incidence of stimulant-related disorders (cocaine and amphetamines) and liver disease. As substance use can induce psychosis, it is unclear whether this finding represents primary psychosis, substance-induced psychosis, hepatic encephalopathy, or a combination of these, further complicating treatment. The risks of opioid misuse and misuse of other substances are increasingly seen as multiplicative. It is therefore unsurprising that this class had the greatest odds of in-hospital death and among the greatest odds of hospitalization. The presence of liver disease as a class-defining variable for this class, as opposed to the class defined by alcohol use, likely speaks to more severe substance use patterns including injection drug use within this class and possibly co-infection with the hepatitis C virus (HCV). Patients with HCV have a disproportionately high prevalence of substance use disorders [40]. The presence of liver disease in this population, along with the emergence of improved treatments for HCV, provides an additional rationale for expanded screening for HCV infection in patients with substance use disorders [41]. While a high prevalence of HCV has been noted in urban Eds and ED-based HCV screening programs have been effective, the presence of HCV screening in this setting is highly variable. The high mortality noted within this class provides an argument for expanding HCV screening programs and programs addressing co-occurring substance misuse within the ED setting. The associations between opioid misuse, psychiatric diagnoses, and polysubstance use are complex and often related to social determinants of health such as housing insecurity, poverty, and race [33,42]. Given this complexity, patients within this subphenotype may particularly benefit from ED-based behavioral approaches in addition to pharmaceutical ones. An understanding of the interplay between these factors is critical for the development of effective public health interventions, which should involve not only those that are opioid-specific but also those which address social and psychiatric needs.
The smallest latent class consisted of opioid-related encounters associated with an alcohol use diagnosis. This class had the highest proportions of male membership, lack of insurance, and membership within the lowest income quartile. Alcohol co-involvement in opioid-involved overdose deaths is common and associated with binge drinking [43]. Given the relatively high odds for hospitalization and in-hospital death within this subphenotype, clinicians should routinely evaluate for problematic alcohol use behaviors in patients with opioid misuse. Effective outpatient interventions for problematic alcohol use should be considered and evaluated within the ED setting. Similarly, policy initiatives in opioid overdose prevention should attempt to address alcohol misuse with accessible intervention programs including resources for those with no insurance or from disadvantaged financial backgrounds.
The identification of the pregnancy subphenotype of opioid-related ED encounters, which represented nearly 12.0% of all opioid-related ED encounters, is consistent with the increased health care engagement and prevalence of opioid use disorder previously described during pregnancy [44]. The incidence of hospital births complicated by opioid use disorder increased more than fourfold from 1999 to 2014 [45]. ED encounters related to pregnancy and opioid use provide opportunities for discussions regarding the treatment options available to prevent both fetal and maternal opioid-related complications during pregnancy. Some of these options have been shown to decrease the severity of neonatal abstinence syndrome and decrease emergency health care utilization in the postpartum period [44,46]. Eds should foster partnerships with clinicians caring for pregnant patients to ensure that appropriate follow-up for both opioid misuse and pregnancy can be assured following discharge from the emergency setting. The detection of opioid misuse in pregnancy represents an opportunity for both the initiation of effective treatments and cost savings for the health care system.
The five-class model has similarities but also important differences when compared with those from other clinical settings. In the hospitalized setting, for example, Liu et al. detected five subphenotypes in their analysis of opioid-related hospitalizations in Pennsylvania [21]. Some subphenotypes were markedly similar to ours, specifically those defined by pregnancy and by polysubstance use with co-occurring psychiatric disorders. However, our model noted additional subphenotypes defined by depression, alcohol use, and chronic pain that were not previously described. These differences are likely due to differences in the variables chosen for analysis and differences between the ED patient population and hospitalized patients. In our analysis, the largest subphenotype of opioid-related ED encounters was defined by chronic pain. This subphenotype likely has overlap with the largest subphenotype identified in hospitalized patients, which was defined by opioid use disorder without co-occurring psychiatric disorders. The recognition of chronic pain diagnoses as a class-defining variable for the largest opioid-misuse subphenotype provides critical context for interventions and barriers to treatment for patients within this subphenotype.
Although not all ED patients with opioid misuse will fall neatly into one of the five subphenotypes identified in this analysis, the study results provide an initial framework through which to understand and tailor care for this patient population. This model may be particularly useful in assisting clinicians in evaluating comorbidities and identifying appropriate pharmaceutical and non-pharmaceutical interventions. This five-class model represents a first step towards the development of personalized treatment approaches for ED patients with opioid misuse and invites further research into the relationship between subphenotypes, specific interventions, and outcomes.
There are several limitations to the current study. The diagnosis codes used as class-defining variables are frequently insensitive and may miss cases of opioid misuse not captured by billing codes. The entry of diagnosis codes by providers requires the detection of a disease process and the entry of the diagnosis into the electronic health record, which contributes to the low sensitivity. Additionally, diagnosis codes do not reliably differentiate between types of opioids or reasons for use. The selection of class-defining variables was based on previous research utilizing LCA in other settings [20,21]. Other potentially informative class-defining variables were possibly excluded from model development. The selected class-defining variables, including chronic pain, are themselves simplifications of complex disease processes and may add some degree of heterogeneity to the subphenotypes described. Individual observations were assigned to a class by their highest posterior probability to facilitate data presentation and interpretation. In practice, however, subphenotypes are not truly mutually exclusive, and patients may have characteristics associated with multiple classes. The models developed were based on nationally representative data. Model validation is needed to ensure that the detected subphenotypes and associated outcomes persist within smaller geographic subdivisions. Last, latent class analysis relies on a combination of the analytic interpretation and evaluation of clinical interpretability; therefore, agreement among subject matter experts may vary.

5. Conclusions

This latent class analysis of ED encounters in the United States demonstrated five distinct and clinically relevant subphenotypes of ED patients with opioid-related encounters: Chronic Pain (class 1); Alcohol use (class 2); Depression and pain (class 3); Psychosis, liver disease, and polysubstance use (class 4); and Pregnancy (class 5). It is important to recognize that the characteristics and needs of ED patients with opioid-related encounters may differ greatly. The detection of these subphenotypes should inform treatment considerations, resource allocation, and future research on the efficacy of patient-tailored ED-initiated interventions for opioid misuse.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijerph19148882/s1, Figure S1: Information criteria-based metrics by number of classes in model; Table S1: ICD-10-CM diagnosis codes defining opioid-related encounters; Table S2: ICD-10 codes; Table S3: STROBE Statement—checklist; Table S4: Patient characteristics odds ratios by latent class.

Author Contributions

Conceptualization, N.C., M.A. and N.S.K.; methodology, N.C., D.L.S., J.A., B.S. and M.A.; software, D.L.S.; validation, N.C. and D.L.S.; formal analysis, D.L.S.; investigation, N.C., D.L.S., C.M.M. and J.A.; resources, D.L.S., H.M.T., M.A. and N.S.K.; data curation, N.S.K. and D.L.S.; writing—original draft preparation, N.C. and C.M.M.; writing—review and editing, all authors; visualization, N.S.K. and D.L.S.; supervision, B.S., H.M.T., M.A. and N.S.K.; project administration, N.S.K.; funding acquisition, H.M.T., M.A. and N.S.K. All authors have read and agreed to the published version of the manuscript.

Funding

This project is supported by research grants from the National Institute on Drug Abuse (UG1-DA049467, R01-DA051464, K23-DA055061), National Center for Advancing Translational Sciences (UL1-TR002398, KL2-TR002387), and National Institute on Alcohol Abuse and Alcoholism (K23-AA024503). Additional support was provided by Rush University Medical Center.

Institutional Review Board Statement

Ethical review and approval were waived for this study because there were no human subjects.

Informed Consent Statement

Not applicable.

Data Availability Statement

Primary data supporting reported results can be obtained from the Healthcare Cost and Utilization Project: Available online: https://www.hcup-us.ahrq.gov/nedsoverview.jsp (accessed on 23 October 2020).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Centers for Disease Control and Prevention, National Center for Health Statistics. Multiple Cause of Death 1999–2017 on CDC WONDER Online Database, Released January 2020. Available online: https://wonder.cdc.gov/wonder/help/mcd.html (accessed on 3 March 2021).
  2. CDC Health Alert Network Increase in Fatal Drug Overdoses Across the United States Driven by Synthetic Opioids Before and During the COVID-19 Pandemic 2020. Available online: https://emergency.cdc.gov/han/2020/han00438.asp (accessed on 3 March 2021).
  3. National Institute on Drug Abuse. Overdose Death Rates. 2021. Available online: https://nida.nih.gov/research-topics/trends-statistics/overdose-death-rates (accessed on 3 March 2021).
  4. Langabeer, J.R.; Stotts, A.L.; Bobrow, B.J.; Wang, H.E.; Chambers, K.A.; Yatsco, A.J.; Cardenas-Turanzas, M.; Champagne-Langabeer, T. Prevalence and Charges of Opioid-Related Visits to U.S. Emergency Departments. Drug Alcohol Depend. 2021, 221, 108568. [Google Scholar] [CrossRef] [PubMed]
  5. Abbasi, A.B.; Salisbury-Afshar, E.; Jovanov, D.; Berberet, C.; Arunkumar, P.; Aks, S.E.; Layden, J.E.; Pho, M.T. Health Care Utilization of Opioid Overdose Decedents with No Opioid Analgesic Prescription History. J. Urban Health 2019, 96, 38–48. [Google Scholar] [CrossRef] [PubMed]
  6. Balio, C.P.; Wiley, K.K.; Greene, M.S.; Vest, J.R. Opioid-Related Emergency Department Encounters: Patient, Encounter, and Community Characteristics Associated With Repeated Encounters. Ann. Emerg. Med. 2020, 75, 568–575. [Google Scholar] [CrossRef] [PubMed]
  7. O’Grady, M.A.; Kapoor, S.; Kwon, N.; Morley, J.; Auerbach, M.; Neighbors, C.J.; Conigliaro, J.; Morgenstern, J. Substance Use Screening and Brief Intervention: Evaluation of Patient and Implementation Differences between Primary Care and Emergency Department Settings. J. Eval. Clin. Pract. 2019, 25, 441–447. [Google Scholar] [CrossRef] [PubMed]
  8. Weiner, S.G.; Baker, O.; Bernson, D.; Schuur, J.D. One-Year Mortality of Patients After Emergency Department Treatment for Nonfatal Opioid Overdose. Ann. Emerg. Med. 2020, 75, 13–17. [Google Scholar] [CrossRef] [Green Version]
  9. Leece, P.; Chen, C.; Manson, H.; Orkin, A.M.; Schwartz, B.; Juurlink, D.N.; Gomes, T. One-Year Mortality After Emergency Department Visit for Nonfatal Opioid Poisoning: A Population-Based Analysis. Ann. Emerg. Med. 2020, 75, 20–28. [Google Scholar] [CrossRef] [Green Version]
  10. D’Onofrio, G.; McCormack, R.P.; Hawk, K. Emergency Departments-A 24/7/365 Option for Combating the Opioid Crisis. N. Engl. J. Med. 2018, 379, 2487–2490. [Google Scholar] [CrossRef] [Green Version]
  11. D’Onofrio, G.; Chawarski, M.C.; O’Connor, P.G.; Pantalon, M.V.; Busch, S.H.; Owens, P.H.; Hawk, K.; Bernstein, S.L.; Fiellin, D.A. Emergency Department-Initiated Buprenorphine for Opioid Dependence with Continuation in Primary Care: Outcomes During and After Intervention. J. Gen. Intern. Med. 2017, 32, 660–666. [Google Scholar] [CrossRef] [Green Version]
  12. Kim, H.S.; Aks, S.E. Take-Home Naloxone and the Need for a Publicly Funded Naloxone Supply. J. Addict. Med. 2021, 16, 1–3. [Google Scholar] [CrossRef]
  13. Waye, K.M.; Goyer, J.; Dettor, D.; Mahoney, L.; Samuels, E.A.; Yedinak, J.L.; Marshall, B.D.L. Implementing Peer Recovery Services for Overdose Prevention in Rhode Island: An Examination of Two Outreach-Based Approaches. Addict. Behav. 2019, 89, 85–91. [Google Scholar] [CrossRef]
  14. Kaczorowski, J.; Bilodeau, J.; Orkin, A.M.; Dong, K.; Daoust, R.; Kestler, A. Emergency Department-Initiated Interventions for Patients with Opioid Use Disorder: A Systematic Review. Acad. Emerg. Med. 2020, 27, 1173–1182. [Google Scholar] [CrossRef] [PubMed]
  15. Gunn, A.H.; Smothers, Z.P.W.; Schramm-Sapyta, N.; Freiermuth, C.E.; MacEachern, M.; Muzyk, A.J. The Emergency Department as an Opportunity for Naloxone Distribution. West. J. Emerg. Med. 2018, 19, 1036–1042. [Google Scholar] [CrossRef] [PubMed]
  16. Kaucher, K.A.; Caruso, E.H.; Sungar, G.; Gawenus, L.; Hurlbut, K.; Sanchez, D.C.; Broderick, K. Evaluation of an Emergency Department Buprenorphine Induction and Medication-Assisted Treatment Referral Program. Am. J. Emerg. Med. 2020, 38, 300–304. [Google Scholar] [CrossRef] [PubMed]
  17. Chan, G.; Gelernter, J.; Oslin, D.; Farrer, L.; Kranzler, H.R. Empirically Derived Subtypes of Opioid Use and Related Behaviors. Addiction 2011, 106, 1146–1154. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Choi, S.; Biello, K.B.; Bazzi, A.R.; Drainoni, M.-L. Age Differences in Emergency Department Utilization and Repeat Visits among Patients with Opioid Use Disorder at an Urban Safety-Net Hospital: A Focus on Young Adults. Drug Alcohol Depend. 2019, 200, 14–18. [Google Scholar] [CrossRef]
  19. Weiss, R.D.; Rao, V. The Prescription Opioid Addiction Treatment Study: What Have We Learned. Drug Alcohol Depend. 2017, 173 (Suppl. S1), S48–S54. [Google Scholar] [CrossRef]
  20. Afshar, M.; Joyce, C.; Dligach, D.; Sharma, B.; Kania, R.; Xie, M.; Swope, K.; Salisbury-Afshar, E.; Karnik, N.S. Subtypes in Patients with Opioid Misuse: A Prognostic Enrichment Strategy Using Electronic Health Record Data in Hospitalized Patients. PLoS ONE 2019, 14, e0219717. [Google Scholar] [CrossRef]
  21. Liu, S.J.; Mair, C.; Songer, T.J.; Krans, E.E.; Wahed, A.; Talbott, E. Opioid-Related Hospitalizations in Pennsylvania: A Latent Class Analysis. Drug Alcohol Depend. 2019, 202, 185–190. [Google Scholar] [CrossRef]
  22. Barton, A.W.; Reinhart, C.A.; Campbell, C.C.; Smith, D.C.; Albarracin, D. Opioid Use at the Transition to Emerging Adulthood: A Latent Class Analysis of Non-Medical Use of Prescription Opioids and Heroin Use. Addict. Behav. 2021, 114, 106757. [Google Scholar] [CrossRef]
  23. Mori, M.; Krumholz, H.M.; Allore, H.G. Using Latent Class Analysis to Identify Hidden Clinical Phenotypes. JAMA 2020, 324, 700. [Google Scholar] [CrossRef]
  24. Grant, R.W.; McCloskey, J.; Hatfield, M.; Uratsu, C.; Ralston, J.D.; Bayliss, E.; Kennedy, C.J. Use of Latent Class Analysis and K-Means Clustering to Identify Complex Patient Profiles. JAMA Netw. Open 2020, 3, e2029068. [Google Scholar] [CrossRef] [PubMed]
  25. Lee, Y.O.; Morgan-Lopez, A.A.; Nonnemaker, J.M.; Pepper, J.K.; Hensel, E.C.; Robinson, R.J. Latent Class Analysis of E-Cigarette Use Sessions in Their Natural Environments. Nicotine Tob. Res. 2019, 21, 1408–1413. [Google Scholar] [CrossRef] [PubMed]
  26. Calfee, C.S.; Delucchi, K.; Parsons, P.E.; Thompson, B.T.; Ware, L.B.; Matthay, M.A. NHLBI ARDS Network Subphenotypes in Acute Respiratory Distress Syndrome: Latent Class Analysis of Data from Two Randomised Controlled Trials. Lancet Respir. Med. 2014, 2, 611–620. [Google Scholar] [CrossRef] [Green Version]
  27. Skogen, J.C.; Bøe, T.; Thørrisen, M.M.; Riper, H.; Aas, R.W. Sociodemographic Characteristics Associated with Alcohol Consumption and Alcohol-Related Consequences, a Latent Class Analysis of The Norwegian WIRUS Screening Study. BMC Public Health 2019, 19, 1364. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. The Hcup Nationwide Emergency Department Sample (NEDS). 2016. Available online: https://www.hcup-us.ahrq.gov/db/nation/neds/NEDS_Introduction_2016.jsp (accessed on 29 November 2021).
  29. Weiss, A.J.; McDermott, K.W.; Heslin, K.C. Table 1, ICD-10-CM Diagnosis Codes Defining Different Opioid-Related Conditions. Available online: http://www.ncbi.nlm.nih.gov/books/NBK538344/table/sb247.tab1/ (accessed on 18 January 2021).
  30. Bonar, E.E.; Coughlin, L.; Roche, J.S.; Philyaw-Kotov, M.L.; Bixler, E.A.; Sinelnikov, S.; Kolosh, A.; Cihak, M.J.; Cunningham, R.M.; Walton, M.A. Prescription Opioid Misuse among Adolescents and Emerging Adults in the United States: A Scoping Review. Prev. Med. 2020, 132, 105972. [Google Scholar] [CrossRef]
  31. Ridley, K.; Coleman, M. The Epidemiology of Amphetamine Type Stimulant-Related Admissions in Albany, Western Australia: 2008–2013. Australas. Psychiatry 2015, 23, 241–244. [Google Scholar] [CrossRef]
  32. Altekruse, S.F.; Cosgrove, C.M.; Altekruse, W.C.; Jenkins, R.A.; Blanco, C. Socioeconomic Risk Factors for Fatal Opioid Overdoses in the United States: Findings from the Mortality Disparities in American Communities Study (MDAC). PLoS ONE 2020, 15, e0227966. [Google Scholar] [CrossRef] [Green Version]
  33. van Draanen, J.; Tsang, C.; Mitra, S.; Karamouzian, M.; Richardson, L. Socioeconomic Marginalization and Opioid-Related Overdose: A Systematic Review. Drug Alcohol Depend. 2020, 214, 108127. [Google Scholar] [CrossRef]
  34. Nylund, K.L.; Asparouhov, T.; Muthén, B.O. Deciding on the Number of Classes in Latent Class Analysis and Growth Mixture Modeling: A Monte Carlo Simulation Study. Struct. Equ. Modeling A Multidiscip. J. 2007, 14, 535–569. [Google Scholar] [CrossRef]
  35. Wurpts, I.C.; Geiser, C. Is Adding More Indicators to a Latent Class Analysis Beneficial or Detrimental? Results of a Monte-Carlo Study. Front. Psychol. 2014, 5, 920. [Google Scholar] [CrossRef] [Green Version]
  36. Stamovlasis, D.; Papageorgiou, G.; Tsitsipis, G.; Tsikalas, T.; Vaiopoulou, J. Illustration of Step-Wise Latent Class Modeling With Covariates and Taxometric Analysis in Research Probing Children’s Mental Models in Learning Sciences. Front. Psychol. 2018, 9, 532. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. von Elm, E.; Altman, D.G.; Egger, M.; Pocock, S.J.; Gøtzsche, P.C.; Vandenbroucke, J.P. STROBE Initiative The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) Statement: Guidelines for Reporting Observational Studies. Lancet 2007, 370, 1453–1457. [Google Scholar] [CrossRef]
  38. Vowles, K.E.; McEntee, M.L.; Julnes, P.S.; Frohe, T.; Ney, J.P.; van der Goes, D.N. Rates of Opioid Misuse, Abuse, and Addiction in Chronic Pain: A Systematic Review and Data Synthesis. Pain 2015, 156, 569–576. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Fendrich, M.; Becker, J.; Hernandez-Meier, J. Psychiatric Symptoms and Recent Overdose among People Who Use Heroin or Other Opioids: Results from a Secondary Analysis of an Intervention Study. Addict. Behav. Rep. 2019, 10, 100212. [Google Scholar] [CrossRef]
  40. Cucciare, M.A.; Cheung, R.C.; Rongey, C. Treating Substance Use Disorders in Patients with Hepatitis C. Addiction 2015, 110, 1057–1059. [Google Scholar] [CrossRef] [Green Version]
  41. Hsieh, Y.-H.; Rothman, R.E.; Laeyendecker, O.B.; Kelen, G.D.; Avornu, A.; Patel, E.U.; Kim, J.; Irvin, R.; Thomas, D.L.; Quinn, T.C. Evaluation of the Centers for Disease Control and Prevention Recommendations for Hepatitis C Virus Testing in an Urban Emergency Department. Clin. Infect. Dis. 2016, 62, 1059–1065. [Google Scholar] [CrossRef]
  42. Ali, M.M.; Sutherland, H.; Rosenoff, E. Comorbid Health Conditions and Treatment Utilization among Individuals with Opioid Use Disorder Experiencing Homelessness. Subst. Use Misuse 2021, 56, 571–574. [Google Scholar] [CrossRef]
  43. Tori, M.E.; Larochelle, M.R.; Naimi, T.S. Alcohol or Benzodiazepine Co-Involvement With Opioid Overdose Deaths in the United States, 1999–2017. JAMA Netw. Open 2020, 3, e202361. [Google Scholar] [CrossRef] [Green Version]
  44. Ahrens, K.A.; McBride, C.A.; O’Connor, A.; Meyer, M.C. Medication for Addiction Treatment and Postpartum Health Care Utilization Among Pregnant Persons with Opioid Use Disorder. J. Addict. Med. 2021, 16, 56–64. [Google Scholar] [CrossRef]
  45. Haight, S.C.; Ko, J.Y.; Tong, V.T.; Bohm, M.K.; Callaghan, W.M. Opioid Use Disorder Documented at Delivery Hospitalization-United States, 1999–2014. MMWR Morb. Mortal. Wkly. Rep. 2018, 67, 845–849. [Google Scholar] [CrossRef]
  46. Jones, H.E.; Arria, A.M.; Baewert, A.; Heil, S.H.; Kaltenbach, K.; Martin, P.R.; Coyle, M.G.; Selby, P.; Stine, S.M.; Fischer, G. Buprenorphine Treatment of Opioid-Dependent Pregnant Women: A Comprehensive Review. Addiction 2012, 107, 5–27. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Radar plots for five subphenotypes of ED patients with opioid-related encounters.
Figure 1. Radar plots for five subphenotypes of ED patients with opioid-related encounters.
Ijerph 19 08882 g001
Table 1. Information criterion for all models.
Table 1. Information criterion for all models.
ModelAICaBICχ2EntropySmallest Class Size (%)
1 class83,278,91283,279,01516,759,968 (501)NANA
2 class75,387,09675,387,3129,191,275 (492)0.90846.6%
3 class71,463,07971,462,4095,424,310 (481)0.93411.7%
4 class69,347,56968,348,0133,396,531 (470)0.94819.0%
5 class68,231,30768,231,8652,326,259 (459)0.97810.6%
6 class67,675,77867,676,4501,796,618 (450)0.9634.5%
7 class67,162,51067,163,2961,086,526 (442)0.9901.9%
8 class66,415,22766,416,127538,386 (432)0.9520.6%
9 class66,124,28966,125,302307,504 (415)0.9681.0%
10 class64,841,15565,842,28232,548 (412)0.9800.4%
AIC, Akaike information criterion; aBIC, adjusted Bayesian information criterion.
Table 2. Heat map of 5-class LCA model by class estimate of candidate class-defining variables.
Table 2. Heat map of 5-class LCA model by class estimate of candidate class-defining variables.
Class 1
Chronic Pain
Class 2
Alcohol Use
Class 3
Depression & Pain
Class 4
Psychosis, Liver
Disease &
Polysubstance Use
Class 5
Pregnancy
Class size48.9%11.0%17.1%11.1%11.9%
Class-defining variable
Chronic pain1.0000.1040.2210.0390.038
Alcohol use 0.0001.0000.0000.0120.000
Psychoses0.0060.0310.0340.3220.001
Depression0.0000.1581.0000.0220.009
Liver disease0.0100.1080.0250.3290.000
Pregnancy0.0020.0010.0030.0031.000
Cocaine use0.0020.0480.0110.1030.001
Amphetamine use0.0010.0150.0070.0850.001
LCA, latent class analysis.
Table 3. Covariate analysis of patient characteristics by latent class.
Table 3. Covariate analysis of patient characteristics by latent class.
Class 1
Chronic Pain
Class 2
Alcohol Use
Class 3
Depression & Pain
Class 4
Psychosis, Liver Disease &
Polysubstance Use
Class 5
Pregnancy
n6,477,2231,377,5262,234,7011,288,1141,565,534
Age (median, IQR)48 (32–62)48 (35–57)52 (35–67)47 (32–60)27 (22–31)
Sex
Female59%30%67%43%100%
Male41%70%33%57%0%
Payer
Medicare29%19%41%32%1%
Medicaid25%33%24%32%55%
Private28%23%24%18%31%
Self-pay13%21%7%14%9%
No charge0%1%0%1%0%
Other5%4%3%3%3%
Median income
Top quartile39%36%34%40%42%
2nd quartile27%25%28%26%27%
3rd quartile20%21%22%20%19%
4th quartile14%19%17%14%12%
Urbanicity
Central metropolitan28%34%24%34%39%
Fringe metropolitan20%22%22%21%20%
250–999 K22%21%24%21%20%
50–250 K11%10%11%9%9%
Micropolitan11%8%12%9%8%
Non-core8%5%7%6%5%
Note: Discharge weights were used in generating descriptive statistics.
Table 4. Patient outcomes by latent class membership.
Table 4. Patient outcomes by latent class membership.
Latent ClassDescriptorHospital AdmissionIn-Hospital DeathED Charges
%OR (95% CI)Count
(per 1000)
OR (95% CI)Median (USD)OR a (95% CI)
Class 1Chronic pain10.4%ref.1.7ref.$2177ref.
Class 2Alcohol use32.9%4.38
(4.36–4.40)
6.51.98
(1.95–2.00)
$28171.26
(1.26–1.27)
Class 3Depression & pain37.0%5.24
(5.22–5.27)
6.72.01
(1.99–2.04)
$26451.22
(1.22–1.22)
Class 4Psychosis, liver disease & polysubstance use37.1%5.33
(5.31–5.36)
18.93.4
(3.39–3.48)
$28811.29
(1.29–1.30)
Class 5Pregnancy12.5%1.24
(1.24–1.25)
<0.10.00
(0.00–0.00)
$26051.07
(1.07–1.08)
All odds ratios unadjusted; ED, emergency department; OR, odds ratio based on three -step procedure; USD, United States dollar; ref, reference category. a Natural log transformed.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chhabra, N.; Smith, D.L.; Maloney, C.M.; Archer, J.; Sharma, B.; Thompson, H.M.; Afshar, M.; Karnik, N.S. The Identification of Subphenotypes and Associations with Health Outcomes in Patients with Opioid-Related Emergency Department Encounters Using Latent Class Analysis. Int. J. Environ. Res. Public Health 2022, 19, 8882. https://doi.org/10.3390/ijerph19148882

AMA Style

Chhabra N, Smith DL, Maloney CM, Archer J, Sharma B, Thompson HM, Afshar M, Karnik NS. The Identification of Subphenotypes and Associations with Health Outcomes in Patients with Opioid-Related Emergency Department Encounters Using Latent Class Analysis. International Journal of Environmental Research and Public Health. 2022; 19(14):8882. https://doi.org/10.3390/ijerph19148882

Chicago/Turabian Style

Chhabra, Neeraj, Dale L. Smith, Caitlin M. Maloney, Joseph Archer, Brihat Sharma, Hale M. Thompson, Majid Afshar, and Niranjan S. Karnik. 2022. "The Identification of Subphenotypes and Associations with Health Outcomes in Patients with Opioid-Related Emergency Department Encounters Using Latent Class Analysis" International Journal of Environmental Research and Public Health 19, no. 14: 8882. https://doi.org/10.3390/ijerph19148882

APA Style

Chhabra, N., Smith, D. L., Maloney, C. M., Archer, J., Sharma, B., Thompson, H. M., Afshar, M., & Karnik, N. S. (2022). The Identification of Subphenotypes and Associations with Health Outcomes in Patients with Opioid-Related Emergency Department Encounters Using Latent Class Analysis. International Journal of Environmental Research and Public Health, 19(14), 8882. https://doi.org/10.3390/ijerph19148882

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop