|Home | About | Journals | Submit | Contact Us | Français|
This study characterized the extent and patterns self-reported drug use among aging adults with and without HIV, assessed differences in patterns by HIV status, and examined pattern correlates. Data derived from 6351 HIV infected and uninfected adults enrolled in an eight-site matched cohort, the Veterans Aging Cohort Study (VACS). Using clinical variables from electronic medical records and sociodemographics, drug use consequences, and frequency of drug use from baseline surveys, we performed latent class analyses (LCA) stratified by HIV status and adjusted for clinical and socio-demographic covariates. Participants were, on average, age 50 (range 22–86), primarily male (95%) and African-American (64%). Five distinct patterns emerged: non-users, past primarily marijuana users, past multidrug users, current high consequence multidrug users, and current low consequence primarily marijuana users. HIV status strongly influenced class membership. Non -users were most p revalent among HIV uninfected (36.4%) and current high consequence multidrug users (25.5%) were most prevalent among HIV infected. While problems of obesity marked those not currently u sing drugs, current users experienced higher prevalences of medical or mental health disorders. Multimorbidity was highest among past and current multidrug users. HIV-infected participants were more likely than HIV-uninfected participants to be current low consequence primarily marijuana users. In this sample, active drug use and abuse were common. HIV infected and uninfected Veterans differed on extent and patterns of drug use and on important characteristics within identified classes. Findings have the potential to inform screening and intervention efforts in aging drug users with and without HIV.
Drug use among older adults in the US is occurring at unprecedented high levels--and growing. The baby boom generation (i.e., people born 1946–1964) has higher prevalences of lifetime drug use compared to all previous generations and appears to be continuing their drug use later into life (Han et al., 2009a). In 2007, all adults aged 50–59 were baby boomers; recent national data reveal that in the same year the prevalence of self-reported past year illicit drug use in this age group nearly doubled from 2002 figures (Han et al., 2009b). Updated projections forecast that by 2020, the number of adults 50 or older with substance use disorder will more than double from 2.8 million in 2002–2006 to 5.7 million (Han et al., 2009b), a projection 30% higher than previously estimated (Colliver et al., 2006; Gfroerer et al., 2003). The enormity of this problem is compounded by under -recognition of substance use in older people (Levin and Kruger, 2002), limited data on the types and patterns of their drug use (Levy, 1998; Lieberman, 2000), the lack of substance abuse treatment services for older adults (Office of Applied Studies, 2007), underutilization of screening tools for primary care practitioners to identify substance abuse problems (Friedmann et al., 2001; Kim et al., 2007; Lewis, 1997), and the medical and psychological comorbidities that often accompany and complicate the health status of older drug users (Geriatric Mental Health Foundation,; Reid and Anderson, 1997; Reid et al., 2002; Simoni-Wastila and Yang, 2006).
At the same time, HIV is increasingly affecting older populations because HIV infected persons are living longer on more effective treatments while HIV incidence in this age group is on the rise. Approximately 10% of new HIV infections in the US occur among people over 50 years of age (Centers for Disease Control and Prevention, 2008). Drug use, especially by injection, is one of the primary modes of HIV and hepatitis transmission, and thus one might expect HIV infected older adults to have higher prevalences of injection drug use (IDU) compared to their uninfected peers. On the other hand, HIV incidence among IDUs has declined by 50% or greater over the past 20 years (Des Jarlais et al., 2000; Lee et al., 2003), and, in many p laces, HIV transmission among IDUs has been associated with sexual transmission rather than unsafe injection (Centers for Disease Control and Prevention, 2009; Kral et al., 2001; Strathdee et al., 2001). Furthermore, most illicit drug use in the US is non-injection (Substance Abuse and Mental Health Services Administration, 2009). Regardless of HIV status and transmission category, patterns of current drug use are poorly understood in older adults (Rabkin et al., 2004; Rosenberg, 1995; Schlaerth, 2007).
The nexus of these three issues — drug use, HIV, and aging — is an understudied yet looming problem. Evidence suggests that HIV may speed aging and exacerbate problems of older age. Active drug use contributes to poor health outcomes for people with and without HIV alike (AETC, 2006; Bruce et al., 2008; Kapadia et al., 2005), and may be particularly harmful if addictive disorders persist later in life. Older adults with HIV may thus be d isproportionately burdened by the negative consequences of ongoing drug abuse (Zanjani et al., 2007). For example, both older age and HIV infection are risk factors for unintentional fatal drug overdose (Darke and Zador, 1996; Tardiff et al., 1997; Wang et al., 2005). Biological mechanisms of this elevated risk may be due to changes in metabolism and body water content in older adults that can lead to higher serum drug concentrations over prolonged periods, thereby increasing risks of unintentional injuries such as drug overdose.
Few studies have been conducted with sufficient power to contrast drug use by HIV status and age; even fewer provide generalizable data (Levy, 1998; Lieberman, 2000) or detailed information regarding medical and psychiatric comorbidities related to substance use problems. One study, limited by sample size, found that lifetime drug dependence was higher for the HIV infected than uninfected individuals, and that expected age-associated declines in diagnoses of dependence were not detected among the HIV infected (Rabkin et al., 2004). Analysis of data from a larger, more diverse sample is needed to assist clinicians and policy makers in developing evidence-based responses to changing patterns of drug use in aging populations.
In this paper, we used data from the Veterans Aging Cohort Study, (VACS)(Justice et al., 2006a), a longitudinal, prospectively consented, observational cohort study conducted among HIV infected Veterans and age-, race- and site-matched HIV uninfected controls receiving medical care at one of eight US Veterans Administration (VA) facilities to: 1) characterize patterns of self-reported drug use in a large and diverse cohort of aging HIV infected and uninfected US Veterans and analyze their construct validity; 2) assess pattern differences by HIV status; and 3) describe demographic, medical and psychiatric pattern correlates within this population. To accomplish these aims, we employed latent class analysis (LCA) (Agrawal et al., 2007; Lubke and Muthen, 2005; Muthen, 2006) to identify patterns of behavior, using all available baseline VACS information. We chose this empirical approach to minimize the potential for measurement error, bias, or oversimplification (Kapadia et al., 2005) inherent when a priori definitions are used to aggregate and categorize complex behaviors.
VACS has been described extensively elsewhere (Conigliaro et al., 2004; Justice, 2006; Justice et al., 2006a; Justice et al., 2001; Justice et al., 2006b); design and population characteristics pertinent to this analysis are reported here. Initiated in 2002, the eight-site VACS aimed to explore comorbidities and behaviors affecting medical outcomes in the broader context of aging and HIV(Justice et al., 2006a; Justice et al., 2001). HIV infected participants were recruited from VA infectious disease clinics; uninfected participants were recruited from VA general medicine clinics. Participants are largely male, representative of the VA population in care for HIV(Justice et al., 2006a). VACS is characterized by a network of administrative, clinical, and research resources that allows linkage of annual self-reported survey data to electronic medical records (EMR), vital statistics, pharmacy prescriptions, and healthcare use. Annual loss to follow up is <5%.
The VACS baseline survey included items on frequency of use (Centers for Disease Control and Prevention, 2002) of marijuana, heroin/opioids, cocaine/crack, and stimulants (amphetamines, speed, crystal methamphetamine). Responses for marijuana and stimulant/cocaine use were: never, no use in the past year, some use in the past year (monthly or less often), and heavy use in the past year (weekly to daily use). For heroin/opioids the response options were never, not in the past year, any past year use. Due to sparseness of data and lack of clinically meaningful differences in their use, stimulant use was grouped with cocaine/crack use. The degree of consequence (high, low) experienced due to drug use was derived from an item within the Michigan Composite International Diagnostic Interview (CIDI-SF) (Kessler et al., 1998): “In the past 12 months, did your use of drugs ever interfere with you r work at school, or a job, or at home?”. Responses were never drug use, use more than one year ago (problems may or may not have been experienced), no problems experienced in the past year, or problems experienced in the past year. Together, these items comprised the indicators for the LCA.
Variables known to be associated with drug use that might differentiate the patterns included socio-demographics, comorbidities (hepatitis C (HCV), medical and psychiatric conditions) (Kilbourne et al., 2001), history of and past year injection drug use, and past year receipt of VA substance abuse and mental health treatment. The selection of candidate covariates of the LCA was a theoretically driven one that was tested statistically in the LCA. The theoretical model informing our analyses was Gelberg’s Behavioral Model for Vulnerable Population s, which amends Andersen’s Behavioral Model to explicitly include domains especially relevant to vulnerable populations (e.g., the homeless) to explain and predict personal health practices (Gelberg et al., 2000). In this model, health outcomes and behaviors are functions of predisposing factors, enabling resources, and need. Predisposing domains traditionally include demographics like age, gender, and ethnicity, but were expanded to include vulnerable domains such as lifestyle instability. In this analysis, we considered HIV positive Veterans as a vulnerable population. Substance abuse, race/ethnicity, and gender are conceptualized as predisposing factors; enabling resources include having a regular care provider, HAART, and income; need encompasses health status, presence of depressive symptoms, and comorbidities; health behaviors act as intermediate outcomes and include adherence to medication; and outcomes are disease progression and death. The focus of the p resent analysis is on the interplay of predisposing, enabling, and need factors as covariates of the LCA; the remaining domains are the subject of future analyses. Socio-demographics included age, race/ethnicity, education level, income, marital status, smoking, body mass index ((BMI) calculated based on self-reported height and weight) and having ever been homeless or marginally housed. Current problem drinking was assessed by the AUDIT-C (Bush et al., 1998). Other covariates included lifestyle instability (Gordon et al., 2006; Kim et al., 2006; McFarland et al., 2003; Shapiro et al., 1999; Song et al., 2000), approximated by past month homelessness and concern about food availability, and quality of life (mental, physical composite scores) which was captured by the Short Form-12 (SF-12) (Ware et al., 1996). Comorbidities were identified u sing International Classification of Diseases, 9th revision (ICD-9-CM) codes from the participant’s EMR, grouped in to clusters of medical an d psychiatric disease diagnoses, as previously reported (Goulet et al., 2005; Goulet et al., 2007). For participants with HIV, concurrent AIDS-associated illness occurring one year before and six months after baseline was defined according to the Centers for Disease Control and Prevention’s 1993 case definition (Centers for Disease Control and Prevention, 1993).
We used LCA, a statistical method for creating meaningful subgroups, or latent classes, based on similarities in responses to a set of observed indicators. It is assumed that the observed variables are independent from one another once conditioned on the latent variable, an assumption known as local independence (Hagenaars and McCutcheon, 2002). LCA produces: (a) the prevalence of each latent class (i.e., prior probability that a randomly chosen person will be in each class) and (b) the probability of response to an indicator, conditional on the latent class. Class categorizations are based on the prevalence of the latent classes and are mutually exclusive and exhaustive. While an individual has a probability of membership in all of the latent classes, theoretically, they belong to one true class. Based on the profile of the response indicators, a putative name for the class can be derived. LCA has a number of advantages. It can incorporate covariates that may influence class membership (Muthen, 2004). LCA has been used previously in evaluations of subjects with substance use disorders (e.g., (Monga et al., 2007)). Also, LCA employs a full information maximum likelihood approach that allows for analysis of all available data, assuming data are missing at random. Therefore we were able to model u sing all available baseline cases instead of deleting those with missing drug use data.
First, we constructed a LCA model u sing the entire sample (regardless of HIV status). We systematically built from simple unadjusted to more complex adjusted models. Unadjusted models considered one to seven classes. To identify the optimal number of latent classes and the best fitting model, we sought to minimize the Bayesian Information Criteria (BIC), maximize entropy, a summary measure that conveys classification quality (Muthen and Muthen, 2007), and yield interpretable latent classes of >1% prevalence. Once a final model was selected, we tested for differences by HIV status and then incorporated covariates in an adjusted LCA model. Results describe the empirically-derived classes and their construct validity. Mplus version 5.1 software (Muthen and Muthen, 2007) was used to fit LCA models.
Adjusted latent classes were compared to other classifications of drug abuse and dependence obtainable through medical records and VACS surveys. These included, respectively, ICD-9-CM diagnoses and CIDI-SF score ≥4, indicating high probability of fulfilling DSM-IV drug dependence criteria (Kessler et al., 1998; Robins et al., 1988). A search of the participant’s EMR on e year prior and six months after the survey visit date was conducted of VA administrative records at the Austin Automation Center (Austin, TX) for ICD-9-CM diagnosis codes 291.00 to 292.99 and 303.00 to 305.99 (Druss and Rosenheck, 2000). Patients were coded dichotomously according to whether or not their EMR noted abuse and/or dependence of alcohol or drugs. Using χ2 and Kappa statistics, we assessed how these related indicators mapped onto the derived latent classes (Orlando et al., 2004, 2005; Tucker et al., 2006) as “validating indicators”.
By including HIV status as a covariate that could have class-varying probabilities and comparing the fit of nested models (Muthen and Muthen, 2007) using χ2 tests we assessed for differences by HIV status. Direct effects (i.e., HIV effects drug use indicators directly versus through the latent class) were tested in similar fashion. Based on these results, subsequent LCA analyses were stratified by HIV status, to characterize disease-specific class differences and correlates.
To examine correlates of class membership, we conducted latent class regression. We adjusted for covariates by employing a stepped modeling approach (Delucchi et al., 2004). First, taking the modal predicted most likely class from the unadjusted LCA as the outcome, class-specific differences were tested for all covariates using bivariate analyses (Wald χ2 tests for categorical variables, ANOVAs or Kruskall-Wallace test for continuous variables), stratified by HIV status. Due to the large number of comparisons made, the overall α level was Bonferroni-corrected to p < .0013 (α = .05/39). Variables found to be significant were then entered into a multinomial logistic regression of the modal predicted most likely class to determine factors independently associated with membership (Ghandour et al., 2008). The resulting significant covariates were considered candidates for the adjusted latent class regression model. Then, HIV-specific LCAs were re-estimated incorporating the candidate variables and checking that covariates excluded earlier in the analysis d id not significantly improve fit, reduce variance, or measurably (≥10%) alter parameter estimates. Last, by testing nested models, any candidate variables that did not help differentiate classes were dropped from the final adjusted latent class regression models. Results p resent the class- and HIV-specific socio-demographic and health status variables. This study was approved by institutional review boards at the Yale School of Medicine and the VA Connecticut Healthcare System, West Haven.
Participants were, on average, age 50 (range 22–86), primarily male and African-American (Table 1). HIV-infected participants were significantly more likely to be African-American, HCV seropositive, not married or living with a partner, and to have past and recent histories of homelessness. Alcohol use was similarly high across HIV status. Drug use was reported by HIV infected and uninfected participants, with a significantly higher prevalence (p<.0001 for all indicators) among the former. Self-reported drug use was higher than that captured by ICD-9 diagnoses (Figure 1) and other drug abuse/dependence criteria (Table 1).
For the full sample (N=6351), a five-class unadjusted model provided the best fit to the observed data (entropy 5-class=0.94, BIC 5-class= 39472.3 vs. BIC 4-class= 40126.5) and defined clinically interpretable and relevant subgroups labeled for convenience based on their profiles of response to the latent class indicators as non -users, past primarily marijuana users, past multidrug users, current high consequence multidrug users, an d current low consequence primarily marijuana users. Local independence was checked by looking at bivariate residuals conditional on class membership associated with each pair of drug use indicators. P-values associated with these tests remained below critical values for the stratified HIV positive model and 2 violations were found with critical values flagging a significant dependence in the HIV negative model. To account for these few large residuals, direct effects were tested between the indicators (cocaine/stimulate use, heroin/opiate use and drug use consequences). Inclusion of the direct effects did not further improve model fit or significantly change the estimated model parameters or probabilities so they were not retained in the final models. HIV status strongly influenced class membership, evidenced by significantly improved model fit (χ2difference=279.56, degrees of freedom =4, p<.0001). HIV status significantly (p<.001) differentiated 18 of the 20 possible class membership combinations (5 comparisons, k-1 classes). Compared to non-users (reference group), HIV infected subjects were more likely to be classified as past multidrug users (OR 2.3), current multidrug users (OR 2.4) or, most likely, as current primarily marijuana users (OR 3.3, all p<.0001).
Based on these strong HIV effects, we conducted unadjusted LCAs stratified by HIV status, and found that the five class model again fit best (HIV infected: BIC 5- class=21215.3 vs. BIC 4-class=21463.2; HIV uninfected: BIC 5-class=21215.3 vs. BIC 4- class=21463.2). When the resultant classes were treated as outcomes in bivariate and multivariable analyses, 19 candidate covariates were identified for consideration in adjusted LCA models. In fitting the multinomial regressions, we detected nonlinearity in effects for age, anticipated from known relationships between aging and drug use and possible generational effects (Rosenberg, 1995; Schlaerth, 2007). At age 60, probability of membership in a current drug use class wa s 0.19 for the HIV infected and 0.087 for the HIV infected; at age 70 these probabilities were 0.066 and 0.032, respectively.
In the final latent class regression modeling step, adjusting the LCA models retained the underlying structures of the unadjusted models, improved model fit, and maintained high entropy (HIV infected: entropy=.92, BIC 5 class=17035.4; HIV uninfected: entropy=0.93, BIC 5 class=15258.5). While the item-response probability profiles were similar by HIV status for the non-users and past multidrug user classes, they differed by HIV status for past primarily marijuana and all current user classes (Table 2). Compared to the HIV uninfected, the probability of reporting high consequence drug use was lower for the HIV infected current multidrug users (0.64 HIV+ versus 0.74 HIV−). For the low consequence primarily marijuana use class, HIV infected participants were more likely to have a history of cocaine/stimulant use (0.43 HIV+ versus 0.32 HIV−) whereas the HIV uninfected were more likely non-users of cocaine/stimulants (0.47 HIV− versus 0.35 HIV+).
For both HIV infected and uninfected, Kappa statistics were moderately high between the high consequence multidrug users and drug dependence according to CIDI (K=0.55 HIV+, K=0.68 HIV−) and ICD-9 (K=0.45 HIV+, K=0.59 HIV−). Lower bounds of the Kappa statistics for all measures of agreement lay above the suggested .40 minimum (Landis and Koch, 1977; Sim and Wright, 2005). The remaining classes had small and negative Kappas with the two drug dependence criteria. Regardless of HIV status, the non-user class had the most negative Kappas (range: −0.32 to −0.18) with the drug dependence criteria. Negative Kappas mean agreement occurred less often than predicted by chance, suggesting genuine d isagreement (Juurlink and Detsky, 2005), which is expected in the non-use and low drug use categories. Most members of the current multidrug users class exhibited abuse/dependence criteria, regardless of HIV status, though this was consistently lower among the HIV-infected. In sum, results suggest construct validity of the latent classes was achieved.
Table 3 displays socio-demographic and health status characteristics of the adjusted drug classes by HIV status. Hispanic/other ethnicity, sex (HIV infected) and education level (HIV uninfected) were not different across classes. Health status variables not statistically significant across class were BMI categories for overweight and underweight, renal disease, pulmonary disease, and, for the HIV uninfected, the medical disease cluster. For HIV infected participants, vascular disease prevalence did not differ by class. All other variables reported were statistically significant across classes by Wald tests using a Bonferroni-adjusted p -value of 0.0013. Class descriptions below summarize findings reported in Tables 2–4 but focus on drug use and health status results only, in the interest of space.
Measures of mental quality of life were generally high for non-users, while for the HIV uninfected, the median physical composite score suggested a lower quality of life. Obesity in this class was the second highest at 15.7% of the HIV infected and 42.0% of the HIV uninfected. Past year alcohol use was common but AUDIT-C scores were low and, as expected, non-users had a low prevalence of substance abuse or mental health treatment. Non-users had high prevalences of chronic medical disease (i.e., diabetes, hypertension and vascular disease) and the lowest prevalence of psychiatric diseases of all classes. HCV was relatively low for HIV infected and uninfected non -users. Compared to all other classes, HIV infected subjects in this class experienced the fewest AIDS-associated illnesses.
The prevalence of this class was similar by HIV status. Past primarily marijuana users shared similar characteristics with the non-user class, regardless of HIV status. Differences were primarily socio-demographic (more non-minorities, higher educational status among HIV+; highest income and employment). Additionally, past primarily marijuana users more often reported more past year alcohol use although their AUDIT-C scores were similar to those of non-users. Medical but not psychiatric diseases were less prevalent than in non-users. The HIV infected had a higher prevalence of AIDS-associated illnesses than in non-users but this prevalence was low relative to other classes with drug use histories.
Overall prevalence of this class was similar by HIV status, but the populations reported different health and drug use experiences. A history of injection drug use was common among the past multidrug users, more so among the HIV infected. Past multidrug users had the lowest reported alcohol use of all classes but there were higher AUDIT-C scores among those continuing to d rink. This class also reported more recent treatment for substance abuse and mental health problems than classes 1, 2 and 5. Past multidrug users suffered from significant health problems, reflected in low SF-12 physical and mental scores. The highest prevalences of obesity and medical diseases were found in this class. Among the HIV infected, chronic medical diseases (i.e., diabetes, hypertension, vascular disease) and HCV peaked in this class.
History of and current drug injection distinguished this class, with current injection more common among the HIV infected. Evidence suggested that some members of this class transitioned from injecting to non-injecting drug use (42.3% ever to 19.8% past year HIV−, 57.9% ever to 27.1% past year HIV+), though available data cannot discern when such changes took p lace. Continued drug use was primarily through smoking marijuana and non-injection use of heroin/opioids and cocaine (including crack). Class members also reported current problem drinking. History of recent substance abuse treatment and methadone treatment was greater in this than for all other classes and they reported poor quality of life scores. Prevalence of obesity was low as were medical diseases associated with obesity, especially diabetes, hypertension, and vascular disease (HIV− only). On the other hand, prevalence of HCV and psychiatric diagnoses-in particular, major depression and bipolar - was high. Co-occurrence of medical and psychiatric disease (i.e., multimorbidity) and AIDS-associated illness was highest in this class.
The distinct differences by disease status in this class were notable. Among the HIV infected in this class, socio-demographics (e.g., non-minority, educational status, income, and stability) and low AUDIT-C scores were similar to prevalences reported for class 1 and 2. For the HIV uninfected, socio-demographics (e.g., employment, history of homelessness, educational status) and high AUDIT-C scores were more akin to those of multidrug users (class 3, 4). Past year alcohol use was reported by approximately 80% of this class and, generally, members were less likely to have been treated recently for substance abuse.
Regarding health status, the HIV-infected in class 5 had the lowest prevalence of obesity and a (non significant) high prevalence of being underweight. In contrast, prevalence of obesity was similar among the HIV uninfected in this class to those not currently u sing drugs. The HIV uninfected had lower physical than mental quality of life composite scores. Regardless of HIV status, class 5 exhibited poor quality of life on the SF-12 mental composite akin to current multidrug users yet unlike them receipt of mental health treatment was low. Medical and psychiatric disease prevalence in the HIV uninfected fell between those of classes 1, 2 and of 3, 4. A lower prevalence of AIDS-associated illness and medical and psychiatric disease was associated with this class’s HIV infected, in line with or lower than classes of users with out a history of multidrug use.
Taken together, the Table 2 and and33 data suggested that, more than a year ago, some of class 5 members transitioned away from hard drug use (i.e., HIV−, HIV+: 0.16, 0.21 probability of past heroin/opioid, 0.32, 0.43 past cocaine/stimulant) and injecting (i.e., 1.2%, 0.8% current injecting, 10.8%, 21.0% history of injecting) to primarily marijuana use. Other responses (e.g., never use) suggested initiation or continued use of primarily marijuana.
This study used LCA to discover subpopulations of drug users and la tent class regression to explicitly incorporate the influence of HIV status on their composition. We observed considerable drug use among an older adult population, one which is often not recognized as HIV infected or drug using (Bhavan et al., 2008; Effros et al., 2008; Goulet et al., 2007; Hartel et al., 2006). A table in the online Appendix1 summarizes the findings including similarities among classes, distinguishing characteristics, and HIV-specific results.
The continued drug use and specific patterns detected in our sample extend findings reported from other studies of aging, drug use, and HIV (Hartel et al., 2006; Hser et al., 2001; Lieberman, 2000; Pitts et al., 2005; Rosenberg, 1995; Schlaerth, 2007; Schlaerth et al., 2004). Similarities in drug use of HIV infected and uninfected populations exist, suggesting that the general medicine clinic controls match well to the HIV infected. However, key differences by HIV status emerged from this analysis, with clinical and public health implications for HIV and aging populations.
We found that drug use and abuse were lower but not uncommon among those of older ages, with illness, or even with HIV. These findings reinforce the need for age-appropriate treatment and harm reduction services and for physicians and caregivers to continue screening for substance abuse and dependence, regard less of their patient’s age or disease status. Health maintenance guidelines for geriatric care (e.g., http://www.americangeriatrics.org) and for people aging with HIV make limited mention of substance abuse screening (Bhavan et al., 2008).
Within the HIV infected and uninfected in the VACS, psychiatric comorbidities, medical comorbidities, and their co-occurrence were more prevalent in people with current and past histories of multidrug use but were far less prevalent among current primarily marijuana users. The “multimorbidity” among Veteran s with HIV an d substance use disorder has been described previously (Goulet et al., 2007). Our results expand upon this work by highlighting the heterogeneity and range of multimorbidity by pattern of drug use. This study d id not capture severity of comorbidities, however, for similar prevalences of medical and psychiatric comorbidity diagnoses, members of the current multidrug users class reported the lowest physical and mental quality of life scores. In addition, AIDS-associated illness was highest in the current multi-drug users class. Results suggest that effectively treating the older adult HIV infected patient in the VA, as elsewhere, w ill entail addressing multimorbidities as well as employing special efforts such as interventions to improve medication adherence and reduce depressive symptoms (Altice et al., 2007; Altice et al., 2004; Safren et al., 2009), to keep AIDS-associated illness at bay in older adults who continue multidrug use.
The low consequence, primarily marijuana users appear to be a unique drug using group. The low prevalence of drug dependence and the absence of negative consequences despite chronic use runs counter to research on cannabis abuse/dependence and their subtypes (Grant et al., 2006). However, other LCA studies have detected the existence of a marijuana-only drug use class amongst nationally representative samples of lifetime drug users (Agrawal et al., 2007). Motivations for primarily marijuana use may be a mixture of transitions away from use of other drugs, (re)initiation, self medication, or evidence of diffusion of medical marijuana as innovative treatment (e.g., appetite stimulation, management of medication side effects (de Jong et al., 2005; Haney et al., 2007; Kuo et al., 2004)), which might be better known to the HIV community and their physicians. Whether this class represents a later stage of drug use, as suggested by similarities to the multidrug user classes for the HIV uninfected, a therapeutic disease-specific class, as suggested by similarities to the non - user and past marijuana users for the HIV infected, or if differences indicate two separate pathways into this class, is open for inquiry.
The discovered classes appear to be valid, as indicated by the association between the clinical diagnostic criteria in the expected directions and by the covariates predictive of class membership. For example, one would have expected, as was found, that recent injecting drug use would be strongly related to the high consequence, multidrug users (class 4) but not to any other class. Adjusting for covariates improved our models’ fit, aligning results closer to the true class membership and highlighting the interplay of predisposing (e.g., race, gender, homelessness), enabling (e.g., income, education), an d need (e.g., comorbidities) factors expected under Gelberg’s Behavioral Model for Vulnerable Populations (Gelberg et al., 2000). Other reported LCAs of drug abuse (e.g., (Agrawal et al., 2007; Monga et al., 2007) have not adjusted for covariates.
Some of our findings concur with previously reported LCAs, which detected classes of marijuana only and multidrug users (Agrawal et al., 2007), while there are also some classes we could not detect, such as prescription drug abusers. Differences may be explained by choice of study sample (age, illness status), exposure assessment (lifetime vs. past year use, types of drugs assessed), choice of model (dichotomous, ordinal, etc.), and inclusion of covariates, among others.
Future studies should explore longitudinal transitions between classes and test hypotheses about class membership. Little is known about how frequently older HIV infected and uninfected drug users move in and out of treatment, initiate and terminate drug use, and switch to other substances (Hser et al., 2007). It has been hypothesized that the increasing isolation experienced as one ages may promote transitions from “harder” drugs (e.g., heroin) which may be used alone to “softer” drugs (e.g., marijuana) used more socially (Fogarty et al., 2007; Osborne and Fogel, 2008; Schlaerth et al., 2004). Future endeavors should also test the utility and validity of our derived drug use classes in predicting key health outcomes.
This study has several limitations. The crucial limitation is its cross-sectional nature. One cannot conclude from our data, for example, that drug use decreases with increasing age but rather that, at the ages represented in VACS, probability of certain class membership is lower at older ages. In cross-sectional designs, it is also not possible to determine whether the classes of current and past drug users are stable or dynamic, whether they represent switching of substances or harm minimizing use, intensification or diminution of particular substances used over time, or whether members are incident, chronic p revalent or relapse case members. Longitudinal analyses are underway to address such questions. In the meantime, this study p resents hypothesis-generating, initial descriptions and comparisons by HIV status of an aging, considerably active drug using population. Also, VACS focuses on an aging population, thus our study may not be representative of young HIV infected Veterans and women, who could have different drug use patterns. The largely male older Veteran population of the VACS may also not be representative of aging adults more generally. However, drug use in older adults occurs predominantly among males (Armstrong, 2007; Han et al., 2009a; Zanjani et al., 2007) and the high degree of drug use found in the VACS HIV uninfected population shares similarities with use among older people in the general population described by others (Han et al., 2009b). Prevalence of club drug and benzodiazepine use was not assessed, hence we may have underestimated the true prevalence of drug use. Prescription opioid abuse was not measured separately from heroin so any potential transition to these substances would not be detected in item-response probabilities. Finally, the extent to which our findings generalize to non-Veterans is unknown. While Veterans in care generally differ substantially from the US population receiving health care, Veterans in care with HIV infection have much in common with others receiving treatment for HIV infection (Conigliaro et al., 2006; Kohli et al., 2005; Lohse et al., 2007; Palella et al., 2006). As evidence of this fact, VACS data was directly compared to other large HIV cohorts in Europe and North America and deemed consistent enough with these other samples to be combined in two high profile cross cohort analyses (Kitahata et al., 2009; Sterne et al., 2009).
The prevalence of any current drug use (classes 4 and 5) in the VACS was 44.5% among HIV infected and 30% among HIV uninfected. Relying on clinical reports of drug dependence would instead p lace the estimate at 25% (ICD-9 code) or 14% (CIDI-SF), with no statistical differences by HIV status. Moreover, there was substantial divergence between EMR and self-reported data (CIDI-SF or LCA) for the past multidrug user class, implying misclassification of past users as current users in the EMR. An LCA approach to exposure measurement serves as a descriptive and analytic tool to find consensus in patterns of behavior reported by participants in a way that incorporates rather than punishes for measurement error and missing data typical in observational study. As the application of LCA continues to expand (Muthen, 1992, 2008), our study suggests that LCA can be used to detect statistically and clinically valid subpopulations of drug users without having to administer lengthy interviews or having to rely solely on medical record-generated proxies of drug abuse/dependence.
Our findings contribute to the increasing literature on drug use, aging, and HIV (Bhavan et al., 2008; Effros et al., 2008; Ettenhofer et al., 2009; Han et al., 2009a; Han et al., 2009b) and are directly relevant to the VA, which provides treatment and care to more than 19,000 people with HIV and 213,000 people diagnosed with Substance Abuse Disorders as well as to Veterans with undetected diagnoses. Specifically, intervention should be prioritized for the problematic multidrug user class, who showed considerable need, poor health, and low quality of life. Results lend support for efforts to integrate mental health, chronic disease, and HIV care; implementation of screening for improved detection of drug abuse and potential referral for specialized care; and expansion of age-appropriate treatment and harm reduction options for drug users.
Disclaimer: The views expressed in this article are those of the authors and do not necessarily reflect the position or policy of the Department of Veterans Affairs.
*A summary table of the results is available with the online version of this article. Please see Appendix A at doi:xxx/j.drugalcdep.xxx …
1Appendix - A summary of results is available as a table with the online version of this article at doi:xxx/j.drugalcdep.xxx …
Publisher's Disclaimer: This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.