Abstract
Current Centers for Disease Control and Prevention (CDC) guidelines for prescribing pre-exposure prophylaxis (PrEP) to prevent HIV transmission are broad. In order to better characterize groups who may benefit most from PrEP, we reviewed demographics, behaviors, and clinical outcomes for individuals presenting to a publicly-funded sexually transmitted diseases (STD) clinic in Providence, Rhode Island, from 2012 to 2014. Latent class analysis (LCA) was used to identify subgroups of men who have sex with men (MSM) at highest risk for contracting HIV. A total of 1723 individuals presented for testing (75% male; 31% MSM). MSM were more likely to test HIV positive than heterosexual men or women. Among 538 MSM, we identified four latent classes. Class 1 had the highest rates of incarceration (33%), forced sex (24%), but had no HIV infections. Class 2 had <5 anal sex partners in the previous 12 months, the lowest rates of drug/alcohol use during sex and lower HIV prevalence (3%). Class 3 had the highest prevalence of HIV (7%) and other STDs (16%), > 10 anal sex partners in the previous 12 months (69%), anonymous partners (100%), drug/alcohol use during sex (76%), and prior STDs (40%). Class 4 had similar characteristics and HIV prevalence as Class 2. In this population, MSM who may benefit most from PrEP include those who have >10 sexual partners per year, anonymous partners, drug/alcohol use during sex and prior STDs. LCA is a useful tool for identifying clusters of characteristics that may place individuals at higher risk for HIV infection and who may benefit most from PrEP in clinical practice.
Introduction
A
Pre-exposure prophylaxis (PrEP) is a once-daily antiretroviral medication (emtricitabine/tenofovir) that is used to prevent HIV infection. Approved by the US Food and Drug Administration in July 2012, PrEP has been demonstrated to be safe and effective in multiple at-risk populations. 6 Among MSM, PrEP was associated with a greater than 90% efficacy in those with detectable levels of medication in their blood. 7 PrEP is also efficacious in other populations, including high-risk heterosexuals and people who inject drugs (PWID). 7 –10 The success of PrEP in clinical trials now gives way to questions regarding how implementation and dissemination of PrEP will be executed in real-world settings. 11 Given that STD clinics routinely provide care to those at highest-risk for HIV infection, they may be an ideal place to implement PrEP programs. 12
Commonly cited concerns about PrEP implementation include side-effects, cost of medications, increased risk behaviors while taking PrEP, and the potential for drug resistance during seroconversion. 11 Current Centers for Disease Control and Prevention (CDC) guidelines recommend reserving PrEP for individuals at highest risk for acquiring HIV. 13 The guidelines suggest that populations who may benefit most from PrEP include: MSM who have had at least one act of condomless anal sex in the last year; MSM who have been recently diagnosed with an STD; and MSM in a serodiscordant relationship with an HIV-positive partner. These broad criteria may be less than optimal for identifying men who are at highest risk for acquiring HIV. Factors including medication cost, available clinical resources, and potential for long-term side-effects may limit broad PrEP implementation to the general population of MSM. Effective use of resources requires identification of the highest-risk subgroups of MSM to help inform PrEP implementation programs.
In this study, we use latent class analysis (LCA) to characterize risk factors for HIV infection among individuals presenting to an urban STD clinic. We build upon previous studies on individual risk behaviors 14 –17 to identify groups of individuals who are at highest risk for HIV infection and candidates for PrEP. The results of this study will help guide individual clinicians on PrEP use and, importantly, policy-level PrEP planning as programs are implemented in the US.
LCA is a statistical analysis that uncovers subgroups, or “classes” defined by distinct response patterns on multiple HIV risk factors. LCA has been used to evaluate substance use, 18 –23 HIV risk factors, 24 –26 medication adherence, 27 and other HIV behaviors, 28 –30 but has not been widely used to elucidate classes of MSM who may benefit most from PrEP.
Unlike traditional regression approaches that evaluate the average effects of risk factors on HIV risk and can add only a few interaction terms before becoming too difficult to interpret, LCA assesses the differential effects of multiple risk factors simultaneously. In essence, LCA evaluates the interaction and patterns of multiple risk factors to identify individuals at greatest risk for HIV. 31 This methodology allows for the identification of specific subgroups of MSM with distinct patterns of high-risk behaviors, which can then inform targeting of HIV preventions interventions such as PrEP.
The objective of this study was to evaluate the behavioral risk profiles of individuals presenting to an urban STD clinic in order to identify distinct subgroups of individuals at highest risk for acquiring HIV infection, and who may benefit most from PrEP.
Methods
Participants
The Miriam Hospital STD Clinic is the only publicly-funded STD clinic in the state of Rhode Island. All individuals presenting for HIV or STD testing complete a one-page intake form, which includes demographics, behavioral risk factors, and other conventional HIV surveillance information. We collected de-identified patient data for all individuals presenting for testing between January 2012 and May 2014. We reviewed the following surveillance data: age, gender, race, ethnicity, risk behaviors, state of residence, insurance status, prior HIV and STD testing, injection drug use, and number of sex partners in the last 12 months (oral, vaginal, and/or anal). Other behavioral data included sex with anonymous partner(s), alcohol or drug use during sex, sex with someone of unknown HIV status, sex with someone who exchanges sex for drugs or money, history of exchanging drugs or money for sex, having ever been forced to have sex, and prior history of incarceration.
Sexual behavior categories included: MSM (exclusively), men who have sex exclusively with women (MSW), men who have sex with men and women (MSMW), women who have sex exclusively with men (WSM), women who have sex exclusively with women (WSW), and women who have sex with men and women (WSWM). HIV (antibody testing), syphilis (treponemal and non-treponemal), gonorrhea (urine nucleic acid amplification testing, NAAT), and chlamydia (urine NAAT) results were reviewed for all individuals. Rectal and pharyngeal NAAT testing for gonorrhea and chlamydia was performed in a small subset of MSM with symptoms, but otherwise not routinely performed.
Bivariate analyses including the chi-square test were used to compare demographics and behaviors between groups. Significance was defined as alpha less than 0.05 and all reported p-values were two-sided. This study was approved by The Miriam Hospital's Institutional Review Board.
Latent class analysis
LCA was used to identify groups of individuals at highest risk for HIV based on specific variables. The purpose of this analysis was to identify groups of individuals who had a higher prevalence of HIV in order to better understand how to implement PrEP. LCA 32,33 is a statistical approach that generates latent classes based on patterns of data, with the goal of grouping similar individuals. A series of LCA models with between one and five classes was tested using the Mplus software package Version 7.1 (Muthén and Muthén, 1998–2013). To avoid the likelihood of converging on a local maximum, 500 start values were generated for each model.
Indices used to determine the best model fit and thus the optimal LCA solution included the Akaike Information Criterion (AIC), sample size adjusted Bayesian Information Criterion (BIC), which are measures of the relative fit of a statistical model to a set of data compared to other models and used for model selection (lower AIC and BIC values indicate better model fit), and the bootstrapped likelihood ratio test (BLRT) for model fit, 34 which tests the null hypothesis of no improvement in fit for the model under consideration compared to a model with one less class.
Entropy, which measures the extent to which classes are distinct from each other, the average posterior probability of class membership, and interpretability of the classes were also considered. Entropy values above 0.8 indicate better distinction of classes. 35
Variables used in the LCA included: date of last HIV test, number of oral sex partners in the last 12 months, number of anal sex partners in the last 12 months, any anonymous partner(s) in the last 12 months, alcohol or drug use during sex in the last 12 months, sex with a partner of unknown HIV status in the last 12 months, prior STD diagnosis in the last 12 months, exchanging sex for money/drugs or sex with someone who exchanges sex for money/drugs, prior injection drug use, having been ever forced to have sex, and prior history of incarceration.
Demographic covariates included age, race, and ethnicity. Outcomes of primary interest included testing positive for HIV or other STDs (i.e., syphilis, gonorrhea, and/or chlamydia) at the clinic visit. After identifying the LCA model with the optimal number of classes based on these criteria, multinomial logistic regression was used to compare classes on demographic covariates, including age, race, and ethnicity, and on outcomes including testing positive for HIV or another STD.
Results
A total of 1723 individuals presented for HIV and/or STD testing during the study period (Table 1). Among these, 75% were male, 52% were less than 30 years old, 70% were Caucasian/white, 18% African American/black, and 23% Hispanic. Sexual preference included MSM (26%), MSMW (5%), MSW (43%), WSM (22%), WSWM (2%), and WSW (1%). Of the whole cohort, 29% reported a prior STD diagnosis in their lifetime. Seventy-six percent of individuals had been previously tested for HIV. Among those tested, 1.5% were found to be positive for HIV, 4.4% for syphilis, 2.6% for gonorrhea and 8.9% for chlamydia. In total, 11.7% of individuals tested positive for at least one STD.
IDU, injection drug use; MSM, men who have sex with men; MSW, men who have sex exclusively with women (MSW); MSMW, men who have sex with men and women (MSMW); STD, sexually transmitted disease; WSM, women who have sex exclusively with men (FSM); WSW, women who have sex exclusively with women (WSW); WSWM, women who have sex with men and women (WSWM).
Compared to non-MSM, MSM were more likely to test positive for HIV (96%, 22/23 of cases in MSM vs. 4%, 1/23 cases in non-MSM, p < 0.01); syphilis (86%, 59/69 vs. 14%, 10/69, p < 0.01); and gonorrhea (59%, 24/41 vs. 41%, 17/41, p < 0.01, Table 1). Non-MSM were more likely to test positive for chlamydia than MSM (81%, 113/139 of cases in non-MSM vs. 19%, 26/139 cases in MSM, p < 0.01).
Of the 538 MSM presenting for testing, 4.1% were positive for HIV. In the bivariate analyses, MSM who tested positive for HIV were much more likely to have a history of prior IDU (14% vs. 3%, p < 0.01) and of being forced to have sex (27% vs. 5%, p < 0.01, Table 2). There was no difference between MSM who tested positive for HIV compared to those who were HIV negative and age group, race, ethnicity, prior HIV testing, number of sex partners in the preceding 12 months (anal or oral or both), anonymous partners, prior STD diagnosis ever or in the preceding 12 months, exchanging sex for money/drugs, sex with someone who exchanges sex for money/drugs, or history of incarceration. There was no difference between MSM who did and did not test positive for any STD and any of the variables tested.
IDU, injection drug use; PEP, post-exposure prophylaxis; PrEP, pre-exposure prophylaxis; STD, sexually transmitted disease.
In the LCA, a total of 449 MSM were included who had complete demographic, behavioral, and clinical data available. A comparison of model fit indices demonstrated that a four-class solution was preferable (Table 3). Compared to models that contained one to three classes, the four-class solution provided the lowest BIC and AIC values, indicating the best model fit, and good entropy (0.81), suggesting adequate distinction of classes, with a significant improvement in model fit over other models by the BLRT. Although entropy was higher for the five-class model (0.84), the decrease in the AIC and BIC between the four- and five-class model was minimal, and the BLRT suggested no improvement in fit for the five-class model over the four-class model (p = 0.20). The average probability of latent class membership was 0.87 for Class 1, 0.89 for Class 2, 0.90 for Class 3, and 0.92 for class 4 (Table 3).
As shown in Table 4, Class 1 (N = 33) was evenly distributed across age and over one-fifth were Hispanic (21%). This class had the lowest percent of individuals who had never had an HIV test (28%), as well as lowest number of reported oral and anal sex partners (82% and 85% with <1 oral and anal sex partners in the past 12 months, respectively). This class also had the highest proportion of past incarceration (33%), history of having been forced to have sex (24%), and sex with a partner of unknown status (67%). A prior STD in the last 12 months was lowest in this group. Class 1 did not contain any HIV positive individuals and had a low proportion (6%) of members testing positive for any STD (Tables 4 and 5).
STD, sexually transmitted disease.
Testing positive for at least one of the following: HIV, syphilis, gonorrhea, chlamydia.
OR, odds ratio; STD, sexually transmitted disease; REF, reference.
Class 2 (N = 142) had a similar distribution of age and race/ethnicity as Class 1. Approximately half (56%) reported never HIV testing. Class 2 had the highest proportion reporting one or two oral and anal sex partners, and no members reporting five or more sexual partners in the past year. Class 2 had the lowest proportion reporting anonymous partners, alcohol and/or drug use during sex, sex with partners of unknown HIV serostatus, or exchanging sex. This class had the second-lowest HIV prevalence at 3%, and 13% of members of this class were positive for at least one STD (Tables 4 and 5).
Class 3 (N = 55) was evenly distributed across age groups, but had the largest number of Caucasian members (85%). Sixty-eight percent had never had an HIV test and no one in the class had been tested for HIV in the preceding year. Individuals in Class 3 also had the largest number of oral sex partners in the past 12 months; 100% reported >3 partners and 69% >10 partners. These individuals also reported the largest number of anal sex partners in the past 12 months; all had had anal sex and 60% had >10 partners. One-hundred percent of individuals in this class reported having had anonymous sex partners.
This class also had the largest number of individuals who had used alcohol or other drugs during sex (76%), exchanged sex for drugs/money or had sex with someone who exchanges sex for drugs/money (38%), had a STD in the preceding 12 months (40%), or had previously injected drugs (22%). Class 3 had the highest prevalence of HIV (7%) and the highest prevalence of any STD (16%).
Class 4 (N = 219) was evenly distributed across age groups, and had the lowest percent of Hispanic individuals (14%). Similar to Class 2, slightly more than half reported that they had never had an HIV test. Individuals in Class 4 more frequently reported having 5–10 oral and anal sex partners than any other category, the second highest number of partners second to Class 3. A larger percent of Class 4 reported sex with an anonymous partner (74%) and a history of an STD in the last 12 months (19%) than any class except Class 4. No members of Class 4 reported a history of injection drug use or incarceration. Class 4 has the second-highest HIV prevalence at 4% and second-highest STD prevalence at 15% testing positive for any STD (Tables 4 and 5).
Discussion
This LCA of MSM presenting for testing at the only publicly-funded STD clinic in Rhode Island helped identify patients who may benefit most from PrEP by identifying groups of patients with differential HIV prevalence. Overall, 4.1% of MSM presenting to the publicly funded STD clinic in Rhode Island tested positive for HIV. MSM were also more likely to test positive for syphilis and gonorrhea compared to non-MSM.
LCA revealed four distinct classes of MSM demonstrating different characteristics and risk profiles. Among these, Class 1 had no HIV-positive individuals and the lowest rates of STDs (6%). Individuals in Class 1 tended to have been previously tested for HIV and had fewer sex partners compared to other classes. Class 1 also had the highest rate of incarceration (33%) and of having been forced to have sex (24%). Class 3 included MSM with the highest rate of HIV (7%) and other STDs (16%). A significant number (68%) of MSM in Class 3 had never had an HIV test and over 60% had >10 anal sex partners in the preceding year; 100% of this class reported having had anonymous sex partners. Class 3 had the largest number of MSM who were reported using alcohol or drugs during sex (76%), exchanging sex for drugs/money or having sex with someone who exchanges sex for drugs/money (38%), having a prior STD in the preceding 12 months (40%), or using injection drugs (22%). Interestingly, 55% (12/22) of newly diagnosed MSM had heard of PrEP, including 100% of those in the highest-risk class who tested positive.
These analyses helped identify groups of MSM who are at highest risk of HIV and may be ideal candidates for PrEP in the setting of an urban STD clinic. Previous evidence has demonstrated significant HIV risk behavior immediately preceding a new HIV diagnosis. 36 In our STD clinic, when individuals are identified as having behaviors associated with the highest risk class (i.e., a large number of partners, anonymous partners, alcohol and/or drug use during sex, etc.), they receive more intensive counseling and education related to PrEP. In addition, identification of these classes has guided practice by indicating that HIV prevention outreach should focus on areas where individuals meet partners (i.e., venues for meeting anonymous partners) rather than in prisons or jails.
The HIV seropositivity rate at the STD clinic of 1.3% overall and 4.1% among MSM is several magnitudes higher than the general US population of 0.3%, 1 but lower than previous estimates of HIV prevalence among MSM in the US with higher prevalence among MSM of color and significant geographic variation. 37 These trends in HIV prevalence, however, mask more complex latent classes, particularly among MSM. Using latent class analysis, we were able to identify a subgroup of MSM (Class 3) who may most benefit from PrEP. Classes 2 and 4 reflected risk factors and HIV prevalence that placed them in the second highest tier for HIV infection risk. MSM from Class 1 had lower HIV prevalence overall and may be less ideal candidates for PrEP.
Interestingly, incarceration and having been forced to have sex were most common among Class 1, which had no HIV infections in the group. The reasons behind a lack of association between incarceration and forced sex with HIV, which are traditionally considered to be risk factors for HIV, 38 –41 are likely reflective of the current HIV epidemic in Rhode Island. Robust HIV screening and treatment programs in the Rhode Island prison system and community-based needle exchange programs have contributed to a decrease in HIV transmission among incarcerated individuals and PWID in Rhode Island; 42,43 these programs may have had positive spillover effects on patients in Class 1. Individuals in this class reported fewer sexual partners, which may reflect an overall lower HIV risk in this population. LCA identifies clusters of characteristics that may identify individuals who are at higher or lower risk of HIV; in this population, individuals with a history of incarceration and forced sex appear to have lower risk for sexual acquisition of HIV, which may explain this result.
The main characteristics among MSM in the highest risk class included having more than 10 oral and anal sex partners in the past year, history of prior STD(s), having anonymous sex partners, and using drugs or alcohol use before or during sex. These findings are consistent with prior studies focused on individual risk behaviors impact on HIV infection 14,44 –46 and suggests that the combination of these risk factors culminates in elevated HIV risk. Our results suggest that individuals in this latent class may benefit most from PrEP and other HIV prevention interventions.
Previous reports among MSM suggest that prior STDs, drug use, and having more than 10 sexual partners in the last year are among the most robust predictors of HIV acquisition. 17,46 More recent data suggest that social and sexual network phenomena may also play an important role in elevating HIV acquisition risks, particularly among MSM of color. 47,48 Although we did not have network-level data, the LCA results suggest that individuals with similar demographic and behavioral patterns may, based on similar characteristics, be in the same social or sexual networks that facilitate HIV transmission.
Our findings suggest that focusing PrEP interventions on MSM in the highest risk groups is important: MSM who reported concurrent and anonymous partners and who reported drug and alcohol use during sex had the highest prevalence of HIV in our study. HIV-negative individuals in those high-risk groups may benefit most from PrEP. Educational efforts targeting the highest risk networks may be important for expanding PrEP uptake among those populations. Future work should also consider whether identification of the highest risk group via a method such as LCA leads to better uptake of PrEP and engagement in PrEP care compared to identification of higher risk groups via traditional routes, and if this method may ultimately lead to a reduced incidence of HIV.
Our study is subject to several limitations. Some behavioral variables used in the LCA were not routinely collected during the first year of the study but were added later. Therefore, complete demographic and behavioral data were available for most, but not all, patients who presented for testing. There was also a relatively low number of newly diagnosed HIV-positive individuals. This may have limited the power of the study to identify significant risk factors for infection. The results of this study arise from patients presenting for care at a STD clinic, which may differs from individuals receiving care at primary care or other clinics.
While the results of this study may not be generalizable to other populations of MSM, the method for identification of subpopulations that may benefit from PrEP could be replicated in other populations. Nevertheless, our findings suggest that LCA helped identify a group of MSM at highest risk for acquiring HIV infection and most suitable for PrEP. Despite its ability to discern interactions and patterns among multiple HIV risk factors, LCA is largely a descriptive method. Future replication of the analysis in larger, independent samples would be useful to determine whether the same profiles emerge reliably.
Furthermore, it may be important to conduct LCA using additional risk factors, as more informative profiles with greater predictive value may emerge. Future work should also consider using LCA to evaluate incident HIV infection in a prospective setting to determine risk factors in a more rigorous manner.
Although CDC guidelines recommend that PrEP should be considered for all MSM who have had any condomless anal sex in the past 6 months or any STD diagnosed or reported in the past 6 months, 13 results from our analysis suggest that a more nuanced public health approach may be useful for identifying ideal candidates for PrEP. We found that combinations of behaviors in MSM, such as having anonymous sexual partners, having more than 10 sexual partners, and using drugs or alcohol during sex, are among the most useful proxies for identifying those who are at highest risk for contracting HIV. These patients may be the most ideal candidates for PrEP.
The findings from this analysis have informed clinical practice in our own PrEP implementation program; this method and our related findings have important implications for clinical practice and for scaling up PrEP interventions around the country.
Footnotes
Acknowledgments
Philip A. Chan and Amy S. Nunn are supported by Grants 1K23AI096923 and K01AA020228 from the National Institutes of Health. Additional support was provided by Grants R25MH083620 and T32DA013911 from the National Institutes of Health, and the Lifespan/Tufts/Brown Center for AIDS Research funded by the National Institute of Allergy And Infectious Diseases (P30AI042853). Philip Chan, Amy Nunn, and Leandro Mena have received grant support from Gilead Sciences, Inc.
Author Disclosure Statement
Amy Nunn has received consulting fees from Mylan, Inc.
