Background. Antiretroviral therapy (ART) reduces human immunodeficiency virus (HIV) RNA load and the probability of transmitting HIV to an HIV-uninfected partner. However, the potential reduction in secondary transmission associated with ART may be offset by the longer duration of infectiousness.
Methods. To estimate the effects of ART on the secondary transmission of HIV among men who have sex with men, we used a previously published state-transition model of HIV disease to simulate the clinical and virologic course of HIV infection among 2 cohorts of men who have sex with men: (1) a cohort of individuals who were not receiving ART and (2) a cohort of individuals treated with US guideline—concordant ART. The model tracked the number of acts of unprotected insertive anal intercourse, transmission risk per act as determined by HIV RNA level, and the number of secondary cases generated in each cohort.
Results. The estimated mean number of secondary transmissions from an HIV-infected individual after 10, 20, and 30 years of infection were 1.9, 2.5, and 2.5, respectively, in the untreated cohort, compared with 1.4, 1.8, and 2.3, respectively, in the treated cohort. The total number of transmissions for the treated cohort began to exceed the total number of transmissions for the untreated cohort 33 years after infection; over the entire course of infection, treatment with ART led to a 23% increase in secondary infections. All estimates of the impact of ART on secondary transmission were sensitive to changes in risk behaviors.
Conclusions. These results suggest that ART must be accompanied by effective HIV-related risk reduction interventions. Programs that target prevention to decrease further HIV transmission are crucial to epidemic control.
It is uncertain whether treatment of HIV infection with antiretroviral therapy (ART) will slow the US HIV epidemic. Because ART is associated with a reduction in HIV RNA levels in both blood and semen [1–3], and because lower levels of HIV RNA are correlated with lower HIV transmission rates [4, 5], ART might curb transmission from individuals receiving therapy. However, because ART does not eradicate HIV RNA [1, 3, 6], the potential reduction in secondary transmission by ART may be offset by the longer duration of infectiousness [7], which may be compounded if HIV-infected individuals do not fully adhere to therapy or reduce their HIV risk behavior.
HIV transmission may occur during every stage of disease, from primary infection to late-stage AIDS. Levels of HIV RNA are highest during primary HIV infection [8, 9]. However, the duration of time spent in this high-infectivity stage is low. One recent report [10] suggests that late-stage disease may also play an important role in transmission. These offsetting dynamics must be taken into account in evaluating the overall impact of ART on HIV transmission.
Several mathematical models have been developed to determine the effect of ART on HIV transmission by estimating the basic reproductive number (R0), which represents the number of individuals that a single infectious person will infect when introduced into a completely susceptible population [11–14]. These models are limited by their simplification of the biological, clinical, and behavioral components of HIV transmission. To address these issues, we modified an existing state-transition model to simulate HIV transmission from HIV-infected individuals to uninfected individuals, accounting for the relationships between ART, HIV RNA level, disease progression, behavioral changes, and transmission risk.
Cost-effectiveness of preventing AIDS complications model (CEPAC). We employed a previously published [15] state-transition simulation model of HIV disease, CEPAC, to estimate the effects of ART on HIV transmission. Primary infection was analyzed separately from chronic and late-stage infection. For chronic infection, men enter the model with a mean CD4+ cell count of 744 cells/mm3 and an HIV RNA level close to set-point, as determined by the Multicenter AIDS Cohort Study [16]. Men were defined as having late-stage infection when they had a CD4+ cell count <50 cells/mm3 and were not receiving ART or when they had a CD4+ cell count <50 cells/mm3 and had experienced failure of their most recent ART regimen.
In CEPAC, HIV RNA level determines the rate of CD4+ cell count decrease, and CD4+ cell count determines the monthly risk of opportunistic infection (OI) and AIDS-related mortality. Simulations within CEPAC may include various treatment regimens, including both OI prophylaxis and ART. When ART reduces the HIV RNA level of an individual, the CD4+ cell count increases at a rate reported in clinical trials [17–19]. When ART fails (defined by virologic rebound, decrease in CD4+ cell count, or a new OI), HIV RNA level slowly reverts to set point, until a new ART regimen is initiated [20].
For each simulated individual, the model tallies the number of months spent in each health state, the number of clinical events, the time spent in each HIV RNA level stratum, and the time of death. After the entire simulated cohort has cycled through the model, life expectancy and mean time spent in each HIV RNA level stratum are estimated. Data on disease progression, ART initiation, and ART efficacy are based on previously published data [21]. We included sequential lines of ART, each with decreasing efficacy of HIV RNA level suppression and CD4+ cell count reconstitution (table 1 and table 2).
The basic reproductive number (R0) by year in the untreated (thick blue line) and treated (red line) cohorts are shown to compare the time required to reach a given R0. The proportion of the original cohort that is alive in the untreated (thin blue line) and treated (pink line) cohort and the proportion of the original cohort receiving ART in the treated cohort (green line) are also shown. This graph is truncated at 60 years after infection.
Summary of input parameters for the model: HIV RNA distribution and CD4+ cell count decrease by HIV RNA level.
Estimation of R. We used the notation R to denote the cumulative number of secondary transmissions from a typical infectious individual by t years after infection, if all contacts were with susceptible individuals. In this notation, R would be equivalent to the traditional basic reproductive number; we consider the cumulative transmission by a given time after infection to focus attention on the nearer (more predictable) term. To estimate R of HIV under different treatment strategies in a population of HIV-infected men who have sex with men (MSM), we added a transmission-tracking component to CEPAC. Monthly acts of unprotected insertive anal intercourse (UIAI) were determined by CD4+ cell count and summed across the entire course of infection. The model then multiplied the number of acts by the probability of HIV transmission per act, which is driven by HIV RNA level at the time of the act. This formula, (number of acts) × (transmission per act) = transmitted cases, generates the number of secondary cases per person per month. R0 was calculated as the mean number of secondary cases up to time t after infection for a cohort of 1 million individuals. Because of the uncertainty of future ART regimens, we focused the primary analysis on the first 10 years after infection, but we also present results for 20 years after infection, 30 years after infection, and the lifetime of the infected cohort.
To estimate the portion of the basic reproductive number attributable to chronic infection, we ran the model under 2 conditions: with no ART and with ART for individuals with CD4+ cell counts <350 cells/mm3. We then added the number of secondary infections attributable to primary infection.
Transmission probability parameters. Because the study population consisted of MSM, we estimated transmission probabilities for unprotected anal intercourse, the primary means of HIV transmission in this population [22], by HIV RNA level. Using data [23] collected from discordant heterosexual partners, we categorized per-contact heterosexual transmission probabilities by HIV RNA level (table 3). To estimate the probability of transmission during primary and late-stage infection, we used estimates from Wawer and colleagues [24, 25]. We compared these heterosexual transmission probabilities to the estimates of per-contact infectivity for UIAI obtained from Vittinghoff et al. [26] and DeGruttola et al. [27] and estimated the probability of HIV transmission during UIAI in MSMs to be 8-fold higher than the probability of transmission observed for male-to-female vaginal intercourse (table 3).
Transmission model parameters: transmission probability per act of unprotected insertive anal intercourse (UIAI).
Sexual behavior parameter estimates.We determined, using the Multicenter AIDS Cohort Study public use dataset [16, 28], that CD4+ cell count, not HIV RNA level, was related to the probability of being sexually active and engaging in insertive anal intercourse. To estimate the number of unprotected sex acts among MSMs stratified by CD4+ cell count, we used data from the Boston Partners Study [22, 29] (table 4).
Transmission model parameters: number of acts of unprotected insertive anal intercourse (UIAI) per month.
Because sexual behavior may differ during primary infection, defined as the first 3 months after seroconversion [30, 31] with HIV RNA level>100,000 copies/mL [30, 32], we obtained a separate estimate of the mean monthly number of acts of UIAI during this period. We estimated this number from data for the 6 months prior to the first positive HIV test result [33]; the estimate was similar to the number of acts in most of the CD4+ cell count categories (table 4). Because sensitivity analyses showed that smoothing out the relationship between CD4+ cell count and the number of acts per month had very little impact on the overall results, we used actual data gathered from the Boston Partners Study categorized by CD4+ cell count. Because we were particularly interested in estimating the effect of primary infection on secondary transmission, we kept primary infection as a separate category.
For the analyses, we assumed that each act of UIAI occurred with a different partner. Although it is unlikely that all individuals change partners with each sexual encounter, the number of transmission events “double-counted” in our analysis should be small because of the rapid rate of partner turnover in the population of interest. In the untreated population, the annual probability of transmission was ∼0.25; if transmission events were randomly distributed, this would amount to ∼4 years between transmissions. This is ∼16 times the mean duration of a partnership (3 months) in the Boston Partners Study [22, 29].
To evaluate the relative impact of the assumptions in the model, we conducted a number of sensitivity analyses. First, we used previously published values for the efficacy of therapy [34–37] and performed sensitivity analyses to evaluate the effect of ART failure rates on the number of secondary transmissions. Second, we explored the effect of varying the 8 : 1 ratio of MSM to heterosexual per-contact transmission probability for all disease stages, including primary and late-stage HIV infection. Third, we considered departures from the primary analysis assumption that men did not change their sexual risk behavior after initiation of therapy. Thus, we evaluated the impact of a potential increase or decrease in the frequency of UIAI both after ART initiation for an individual and after the availability of ART to the population. Fourth, a secondary case that occurs early in the course of infection would contribute more to the epidemic than a case that occurs later [38]. Therefore, we explored the effect of decreasing the impact of secondary cases exponentially with time, so that a case occurring earlier in the epidemic is weighted more than a case occurring later in the epidemic. Finally, because in some studies risky sexual behavior has been found to be negatively associated with age [39–41], we examined the effect of exponentially decreasing risk behavior with age on the lifetime estimates—an analysis that was technically equivalent to that which reduced the impact of secondary cases exponentially with time. We examined these parameters both individually and jointly to understand their impact.
Baseline model results. On average, the untreated cohort of MSM had a much shorter survival (mean survival, 145.4 months), compared with the cohort treated with ART (mean survival, 366.7 months), and engaged in fewer total lifetime acts of UIAI (5). For both the untreated and treated cohorts, an average of 0.36 secondary cases occurred per person during primary infection. Within the first 10 years of infection, primary infection accounted for 19% of the secondary cases transmitted by the untreated cohort, compared with 26% of the secondary cases transmitted by the treated cohort. With ART still widely used 10 years after infection, <0.01% of secondary cases occurred during the last 6 months of life, compared with 5% in untreated men (table 5).
Ten-year and lifetime projected outcomes for untreated and treated cohorts per 1 million population.
Secondary transmissions from an untreated cohort and a treated cohort. We estimated an untreated R of 1.9 and a treated R of 1.4, suggesting that ART reduces transmission in the first 10 years after infection. ART has a strong, beneficial impact on secondary transmissions over the first 20 years after infection, with an untreated R of 2.5, compared with a treated R of 1.8. However, at longer time horizons after infection, the benefit of treatment decreased, because viral load reductions were offset by longer duration of infectiousness (untreated R = 2.5 and treated R = 2.3 at 30 years), assuming that no new HIV treatments become available. The treated and untreated R were equivalent at 33 years, because of the eventual failure of ART and the longer survival in the treated cohort. Over the entire lifetime of the 2 cohorts, treatment led to a 23% increase in R0 = R (figure 1).
The graph of R (figure 1) shows that, in the untreated cohort, R increases steadily but begins to level off as individuals die, reaching a plateau by 26 years. In contrast, the R for the treated cohort shows initial, steady growth similar to that of the R for the untreated cohort; growth then slows after ∼10 years, as individuals begin treatment and their HIV RNA level decreases. The slope of R begins to increase again after ∼20 years of infection, as a large proportion of individuals exhaust available treatment options.
Sensitivity analyses. Because most of the assumptions that we tested were similar across both treated and untreated populations, they had no effect on the relative magnitude (ratio) of R for untreated and treated cohorts. Changing the rate of ART failure had little effect on the number of secondary transmissions at 10 years after infection (table 6). Population-wide changes in sexual activity levels increased or decreased R0 by the same factor regardless of ART use, as did changes in the assumption that per-contact infectivity of UIAI was 8-fold higher than vaginal intercourse. The probability of transmission during primary, chronic, and late-stage infection had a multiplicative effect on the number of secondary cases during each of these disease stages (table 6). However, changing the overall probability of transmission during primary infection and late-stage infection had little effect on the number of secondary cases, because together they account for ∼25% of all secondary cases that occurred within the first 10 years of infection. Comparatively, there was a greater impact on transmission when we increased or decreased the probability of transmission during the chronic stage of infection, because the majority of secondary cases arise during this long disease stage.
Changes in risk behavior because of the availability of ART or because of starting ART had a substantial effect on R in the treated group (table 6). For example, if all of the men in the treated cohort decreased risky sexual behavior by 50%, compared with the untreated cohort, simply because ART was available and regardless of whether they had yet begun treatment, R would decrease to 0.7, indicating a 63% decrease in secondary cases, compared with no treatment. If these men doubled their risky sexual behavior, there would be a mean of 2.8 secondary cases per person after 10 years (i.e., 1.5 times the number of cases found without treatment). However, the effects of changing sexual behavior were not as substantial over a 10-year period if we assumed that only the men receiving therapy changed their behavior (table 6).
When we assumed that individuals decrease their sexual activity as they age, there was a disproportionate reduction in secondary cases among patients in the treated cohort, who live longer and generate relatively more secondary cases at older ages. The overall difference in the number of secondary cases in the treated cohort, compared with the untreated cohort, decreased with greater decreases in sexual activity with age, until—at a 3% annual decrease in sexual activity—there were more secondary cases in the untreated cohort than in the treated cohort (table 6). These same analyses could assess the impact of changes in secondary transmission in the context of a growing epidemic [8, 38] (table 6).
Using a detailed natural history and treatment model of HIV disease, combined with monthly probabilities of sexual activity and per-act transmission risk, we estimated that the use of ART resulted in a reduction in the number of secondary HIV transmissions from 1.9 to 1.4 transmissions per person during the initial 10 years after infection, assuming no increase in risk behavior and no changes in available therapy. These results support more-aggressive and more-active case-finding approaches, such as routine voluntary HIV screening and appropriate ART. However, much of this benefit may be reduced over the lifetime of HIV-infected patients because of the combined effect of increased survival and the eventual failure of ART, increasing R0 from 2.5 to 3.1. However, this increase is not inevitable and may be attenuated by the identification of new and effective ART regimens and other effective treatment options, as well as by decreased sexual activity resulting either from the aging of the ART cohort or from specific risk-reduction interventions.
An alternative to estimating R0 would be to estimate the effects of ART on incidence. Although these estimates are closely related, the effects of treatment on incidence depend on a number of additional factors, such as the current phase of the epidemic, the types of sexual mixing in the population, and other aspects of the epidemic in a particular population. We have chosen to limit our analysis to effects on R0 to make use of the detailed natural history model CEPAC to calculate a single measure of infectiousness.
The values of R0 obtained here (and in other existing estimates of treatment effects) represent averages across the population. It is well known that behavioral heterogeneity increases the effective value of R0 above its mean value in the population, because highly sexually active individuals are more likely to become infected and to transmit infection [42]. As with many of the other uncertainties in model inputs, this issue affects the magnitude of R0 in both treated and untreated populations, but it should not affect the relationship between them.
A published mathematical model has concluded that widespread ART use could eradicate the HIV epidemic in the MSM population in San Francisco [12, 13] if the incidence of unsafe sexual behavior remained stable. Velasco-Hernandez et al. [13] assumed a 2-fold to 100-fold decrease in infectivity for individuals treated with ART and estimated an R0 of 1.0 with the assumption of no change in sexual risk behavior. Our model provided a more detailed method for estimating time-dependent reductions in infectiousness and the increased duration of infectiousness achieved by ART by explicitly incorporating biological, clinical, and behavioral components. Our results, although still predicting that treatment would reduce infectiousness over the first 10 years, are less optimistic than earlier efforts regarding the likelihood of HIV eradication [43].
Gray et al. [14] estimated a 10-year R0 of 1.44 for an untreated population, which decreased to 1.09 when HIV-infected individuals with HIV RNA levels>55,000 copies/mL were treated with ART. Our 10-year estimates of R0 differed from those reported by Gray et al. [14], who simulated a heterosexual population and modeled ART use as a function of HIV RNA level rather than of CD4+ cell count.
Several studies have focused on the impact of primary [8, 9] and late-stage infection [10] on transmission. We found that the number of secondary cases that occurred during primary infection was a small percentage of the total number of secondary cases that occurred during the first 10 years of infection. Late-stage infection contributes even fewer secondary cases. Therefore, the extended duration of chronic infection contributes the largest proportion of secondary cases and provides more opportunities for preventive interventions.
We assumed that all sexual acts were with a seronegative partner consistent with the definition of R0. To calculate the expected number of secondary cases in a population in which only a percentage of partners were seronegative, R in both untreated and treated individuals would be multiplied by the same percentage if we assume the probability of encountering a seronegative partner is equal in the treated and untreated cohorts. Thus, the ratio of secondary transmissions for an untreated versus treated cohort remains the same over the cohort's lifetime.
Another important assumption in the model was that the number of acts of UIAI each month does not decrease with age. Variations in this assumption had the largest effect on lifetime estimates of R0. Because treated individuals live longer, a decrease in sexual activity with increasing age has a more powerful impact on the treated population. Sexual activity would have to decrease by ∼3% per year before the R0 of the treated cohort is less than that of the untreated cohort.
Because a secondary case early in the course of infection contributes more to future cases than does a case that occurs later [8, 38], the same sensitivity analysis used for age could be used to consider the impact of time-weighted transmission. In the context of an epidemic that is growing by ≥3% per year (but in the absence of age-related behavior change), the contribution of a treated person to epidemic growth would be less than that of an untreated person. Put another way, the benefits of treatment are greater in a growing epidemic, because the delay in transmission by treated persons reduces their contribution to the epidemic, even in situations in which the total number of secondary cases resulting from transmission from a treated person might be greater.
An implicit assumption of this model is that the relationship between plasma HIV RNA level and per-contact infectiousness was the same for treated and untreated individuals. Among individuals receiving ART, HIV RNA levels have been shown to be lower in semen than in plasma [44]; if this phenomenon can be confirmed, the effect of treatment in reducing secondary transmission might be greater than was observed in our model, which was driven by reductions in plasma HIV RNA level. A priority for future research would be to measure the probability of transmission given seminal HIV RNA levels.
In conclusion, the relationship between ART and lifetime infectiousness is complex, with an initial benefit over the short term but an eventual increase in total infections transmitted over the lifetime of an HIV-infected cohort. Our best estimate of R0 in patients receiving ART always exceeds 1.0, suggesting that ART alone will not eradicate the HIV epidemic. Thus, it will be important to implement complementary programs that target reduction in secondary transmission, in addition to ART, to further decrease HIV transmission.
Data in this manuscript were collected by the Multicenter AIDS Cohort Study. Principal Investigators included Joseph B. Margolick and Lisa Jacobson (The Johns Hopkins University Bloomberg School of Public Health; Baltimore, MD), John Phair (Howard Brown Health Center and Northwestern University Medical School; Chicago, IL), Roger Detels (University of California, Los Angeles; Los Angeles, CA), and Charles Rinaldo (University of Pittsburgh; Pittsburgh, PA).
Financial support. National Institute of Allergy and Infectious Diseases (T32-AI-07433, R01-AI-058736, R37-AI-42006, K24-AI-062476, K23-AI-01794, K25-AI-50436), the National Institute of Mental Health (R01-MH-65869-03), and the Doris Duke Charitable Foundation (Clinical Scientist Development Award). The Multicenter AIDS Cohort Study is funded by the National Institute of Allergy and Infectious Diseases, with additional supplemental funding from the National Cancer Institute and the National Heart, Lung, and Blood Institute (UO1-AI-35042, 5-MO1-RR-00722 [GCRC], UO1-AI-35043, UO1-AI-37984, UO1-AI-35039, UO1-AI-35040, UO1-AI-37613, and UO1-AI-35041).
Potential conflicts of interest. All authors: no conflicts.
IDSA Members: For your free access to this journal, log in via the IDSA members area.
Open access options for authors visit Oxford Open
This journal enables compliance with the NIH Public Access Policy