HIV Drug-resistant Strains as Epidemiologic Sentinels
http://www.100md.com
《传染病的形成》
University of California, Berkeley, California, USA
Gladstone Institute of Virology and Immunology, San Francisco, California, USA
California Department of Health Services, Berkeley, California, USA
Observed declines in drug resistance to nucleoside reverse transcriptase inhibitors among persons recently infected with HIV-1 in monitored subpopulations can be interpreted as a positive sign and lead public health officials to decrease efforts towards HIV prevention. By means of a mathematical model, we identified 3 processes that can account for the observed decline: increase in high-risk behavior, decrease in proportion of acutely infected persons whose conditions are treated, and change in treatment efficacy. These processes, singly or in combination, can lead to increases or decreases in disease and drug-resistance prevalence in the general population. We discuss the most appropriate public health response under each scenario and emphasize how further data collection and analyses are required to more reliably evaluate the observed time trends and the relative importance of forces shaping the epidemic. Our study highlights how drug resistance markers can be used as epidemiologic sentinels to devise public health solutions.
In recent years investigators have begun monitoring the HIV epidemic by reporting changes in the proportion of newly infected persons who are carrying an HIV-1 drug-resistant strain, i.e., the primary or acute resistant fraction. Several studies report decreases in this primary resistant fraction (1–7), including transient decreases (8–14). Unfortunately, none of these studies included precise longitudinal data on the exact number and type of infected persons or of the fraction of the total population that is screened for acute infection or resistance, and as we show here, making direct interpretations from data collected from a subset of the population can lead to erroneous conclusions. Given the potentially serious clinical implications of drug resistance for HIV-infected persons, public health officials and other authorities need to know whether the decline in drug resistance among acutely infected persons in the monitored subpopulations corresponds to a real decline in drug resistance in the general population and whether this effect is sustainable over time. The decline might be interpreted as a positive sign caused, for example, by less high-risk activity by HIV-positive persons infected with a drug-resistant variant. This explanation could lead public health officials to decrease their support for HIV surveillance and prevention programs targeted at impeding the spread of drug-resistant HIV strains such as drug-resistance testing or adherence counseling. The amount of resources that should be dedicated to drug resistance monitoring and reporting is a controversial issue in light of the recent isolation of a highly virulent multidrug-resistant strain in New York City (15). Here we show how drug-resistance data can offer not only clinical information regarding appropriate treatment regimens for individual patients, but also critical insights into an epidemic's course.
Treatment History of HIV-1 and Its Impact on High-Risk Behavior
Different types of drugs have been developed to fight HIV. Zidovudine (AZT), a nucleoside reverse transcriptase inhibitor (NRTI), was first administered in 1987, and until 1995, monotherapy or dual therapy with NRTIs were the only treatments available. The first protease inhibitor (PI), saquinavir, was approved for treatment in 1995, followed closely in 1996 by a nonnucleoside reverse transcriptase inhibitor (NNRTI), nevirapine. These new drugs generated a major change in the treatment strategy against HIV—highly active antiretroviral therapy (HAART)—that coincided with the start of the monitoring periods in several of the studies mentioned above (1995–1996). With HAART, at least 3 drugs are administered at the same time, which substantially reduces viral load and, compared to results of earlier regimens, increases the life expectancy of patients. These advantages follow because the mutations necessary to confer resistance to HAART are generated at a slower rate and are lost more rapidly than those conferring resistance to monotherapy or dual therapy. Moreover, viral strains resistant to HAART are not as efficient at completing their own life cycle (e.g., their replication rates are lower), they may generate less illness and lower proportion of deaths among infected persons, and the viral strains are less likely to be transmitted to other persons.
The primary resistance time trends observed for NRTIs do not match those observed for the other 2 drug types. For example, in North America some researchers report a decrease in the proportion of persons recently infected with a drug-sensitive HIV strain resistant to NRTIs; this decrease is followed by an increase and subsequent decrease (11–13). The pattern is particularly noticeable in the study by Grant et al. (12), in which NRTI genotypic resistance decreased from ≈30% in 1997 to 5% in 1999, rose to 20% in 2000, and fell to 15% in 2001. Little et al. also found this trend in primary resistance to NNRTIs (11). These 2 studies (11,12) also documented a steady increase in the proportion of persons newly infected with a virus strain resistant to PIs. Here we focused on identifying the likely forces responsible for the time trends exhibited by viral strains resistant to NRTIs. We did so because their time trends are expected to provide better insight into the long-term dynamics of the epidemic than strains resistant to PIs or NNRTIs, given that NRTIs have been administered to more HIV-infected persons and for a longer period than the other 2 types of drugs (16,17).
Treatment optimism after the initial successes of HAART likely affected the subsequent dynamics of HIV because these favorable treatment outcomes led some persons to increase their high-risk behavior. Later it became apparent that HAART does not completely eliminate HIV from an infected person or impede its transmission. Moreover, when HAART first became implemented, the best strategy was believed to be "hit hard, hit early," because the medical community was trying to limit the expansion of HIV within an infected person's body and ameliorate the gradual deterioration of the patient's immune system. However, HAART can have considerable negative side effects, which affect the functioning of the gastrointestinal system, renal system, pancreas, and liver and produce changes in blood count, allergies, lactic acidosis, and other problems. As a result, treatment began to be delayed to balance the following factors: 1) containing the viral load, 2) minimizing the risk of drug-resistant mutants developing by limiting the amount of treatment time, and 3) reducing negative side effects.
Modeling Drug Resistance in HIV
Empirical studies have shown that antiretroviral treatment (ARV) produces substantial changes in the viral dynamics at the within-host level that translate into substantial changes at the between-host level (8,18,19). Mathematical and computational models permit us to create simplified versions of complex realities that we can manipulate to further our understanding of their dynamic behavior. Consequently, numerous theoretical studies have investigated the impact of drug therapy on HIV dynamics at both levels (e.g., 20). Initial HIV treatment models (e.g., 21–23) addressed how ARVs might affect the infectiousness of treated persons, and the spread of HIV and its disease-induced deaths. The magnitude of the public health threat created by drug-resistant HIV strains was only recognized later. As a result, Zaric et al. presented a novel model that showed that adhering to treatment regimens would discourage the emergence of multidrug-resistant HIV strains in heterogeneous populations (24). Blower et al. developed a relatively simple but revealing deterministic compartmental framework (25) that has served as the reference point for most of the modeling studies subsequently done to investigate the effect of ARVs on disease incidence and prevalence, drug-resistance transmission and prevalence, AIDS death rate, and the potential to eradicate the HIV epidemic (26–31). Dangerfield et al. built a detailed HAART treatment model that accounts for persons in all 4 HIV stages, and the last is partitioned in early- and late-stage AIDS (32). They investigated the effects of HAART on HIV incidence and prevalence, assuming different average efficacious periods and assimilation times for HAART, different infectivity probabilities when receiving HAART, and different increases in the mean number of sexual partners.
Blower and Volberding reviewed mathematical studies used to understand the dynamics of a drug-resistant HIV epidemic, predict the incidence and prevalence of drug-resistant HIV strains, evaluate cost-benefit strategies, and assess the impact of public health policies (33). The general approach to these studies had been to construct a descriptive simplification of the epidemic by identifying critical categories and processes and to use this structure to make predictions, given a set of assumptions regarding the parameter values. In this regard, several studies have characterized the epidemic's trends as monotonic, including the fraction of new HIV infections that are drug resistant (e.g., 16,25,28,34). Goudsmit et al. conducted an analysis in which including changes in treatment rates explained the nonmonotonic trends of zidovudine resistance observed in a cohort of newly infected homosexual men enrolled in the Amsterdam Cohort Study (1).
In our previous study (35), we extended the basic modeling framework detailed in (25) to incorporate additional complexity, including 2–3 separate categories of acutely infected persons, depending on whether a person was infected with a drug-sensitive HIV strain, a strain resistant to monotherapy, or a strain resistant to triple-drug therapy. In doing so, we were able to distinguish among acutely infected persons, who are clinically and epidemiologically distinct from uninfected and chronically infected persons (e.g., we can consider them to engage more frequently in high-risk behavior (2) than chronically infected persons). We also counted these categories separately and tracked their temporal trends and better channeled the different categories through the model, according to the different processes acting on them (such as a decrease in the proportion of persons receiving treatment among those recently infected, an effect that did not occur in persons in the chronic phase). We also created 10 subcategories in the chronically infected stage to more accurately represent the progression of persons from the acute stage of infection to AIDS (35). Moreover, because our intent was to explain observed trends rather than to make future predictions, we adopted a specific approach that consisted in altering individually the value of each parameter during a given simulation (as was done in [1] with treatment rates of all persons), rather than running simulations with a set of fixed parameters and comparing outcomes across runs.
Primary Resistant Fraction
The most direct explanation for the decrease in the observed proportion of newly infected persons infected with NRTI-resistant HIV-1 strains is that resistance to NRTIs in the recently HIV-infected population is declining. Unfortunately, the data need to be further evaluated because the HIV infection status of every person in the general population has not been monitored. We therefore do not have absolute numbers for these time trends, but only the relative numbers obtained from monitored subpopulations, which consist of consenting persons enrolled in research programs at specific locations. To be eligible to participate in these programs, patients had to display symptoms typical of an acute HIV seroconversion syndrome or have recently engaged in risky activities that could have placed them at risk of contracting HIV. Accordingly, what has decreased is the fraction of drug-resistant carriers among the pool of recently infected HIV patients who are willing to participate in particular research programs and attend clinics involved in these studies. The time trends exhibited by the variable representing the actual counts of all the newly infected persons who are carriers for a resistant strain in the general population may or may not be a direct match to those of the monitored subpopulations.
These points are best illustrated by considering the fraction of recently infected persons who are carriers of a drug-resistant strain (primary resistant fraction, FR), defined mathematically as
,
where S is the number of persons initially infected with a drug-sensitive strain, and R is the number of persons initially infected with a drug-resistant strain. This fraction may decrease because R decreases (fewer newly infected persons have a strain that is drug resistant), or because S increases (more newly infected persons carry drug-sensitive strains). If both S and R increase or decrease by the same proportion, FR remains unchanged. As explained below, the benefit of using this variable's time trends to further our understanding of the past, present, and future of the HIV epidemic is that underlying alterations in the relative values of drug-sensitive and drug-resistant strains may arise from a variety of mechanisms with critically different epidemiologic outcomes.
To determine which processes could have caused the observed decrease in FR, we built a mathematical model of HIV transmission ; a more mathematically detailed explanation of our analysis can be found in our previous study (35). We then simulated the epidemic using this model and varied each of the parameters shown in Figure 1 (e.g., the average number of high-risk contacts in 1 year, the likelihood of transmitting HIV given a high-risk contact, the fraction of persons with acute or chronic HIV infection that are placed on treatment each year, the likelihood of generating or losing drug resistance in 1 year). Once we determined which processes can cause a decrease in the acute fraction infected with a drug-resistant strain, we evaluated whether the process had occurred in industrialized countries in recent years. If so, we could consider the process as a potential contributor to the observed trends. As a result, we identified 3 independent processes that caused a decrease in FR and were consistent with the history of the HIV-1 epidemic in industrialized countries from 1995 to 2001: 1) overall increase in risky behavior, 2) decrease in the fraction of individuals in the acute phase who are placed on treatment, and 3) increase in the efficacy of treatment. Goudsmit et al. also found that discontinuation of monotherapy with zidovudine in 1996 explained the observed drop in zidovudine resistance in patients newly infected with HIV in the Amsterdam Cohort Study (1).
Course of the HIV Epidemic
How do these 3 processes (increase in high-risk behavior, treatment delay, and greater treatment efficacy) impact the course of the HIV epidemic Are these effects comparable or are they sufficiently different such that the policy implications will vary according to which one we interpret to be the leading cause for the observed decline To address these questions, we determined the short- and long-term changes induced by these processes on 2 variables of critical public health importance for which we do not have reliable measurements: prevalence of disease (i.e., overall fraction of persons infected with HIV in the general population) and prevalence of drug resistance (i.e., overall fraction of persons infected with a drug-resistant HIV strain in the general population). By doing so, we use time trend changes in the relative prevalence values of 2 viral strains (i.e., the fraction FR defined above) to make inferences about changes in the absolute values of disease and resistance prevalence in the population (i.e., actual counts of infected persons and carriers of drug-resistant infections in the general population). Figures 2B and 2C show that the 3 factors we identified as causing a decrease in the primary resistant fraction are predicted to force different and permanent, long-term changes in disease and drug-resistance prevalence. A synopsis of our findings is provided in the Table.
Policy Implications and Comparison with the Data
Our results demonstrate how a decrease in the fraction of persons recently infected with a drug-resistant HIV strain can occur not only when the epidemiologic conditions improve (i.e., disease and drug-resistance prevalence in the population decrease), but also when the epidemic worsens (i.e., disease and drug-resistance prevalence increase). The 3 processes that can generate the decrease in primary resistance are not mutually exclusive, and probably all have contributed to the observed time course of FR. The challenge now is to identify which one has had the greatest effect on the recent trends of disease and drug-resistance prevalence in the HIV-1 epidemic.
If an increase in high-risk behavior has dominated HIV-1 epidemiology since the onset of HAART, then the decrease in primary resistance, counter to intuition, signals a worsening of the epidemic: a greater number of persons may have become infected, and a greater number of persons may be infected with a viral strain resistant to drug therapy. Other studies have obtained similar conclusions (e.g., 28,29,31,32). If this is the case, the public health response to the decrease in drug-resistance levels among the acutely infected should be to expand programs aimed at reducing high-risk behavior.
Determining the most appropriate public health response is difficult if a decrease in the fraction of acutely ill persons receiving treatment is the main driving force of the HIV epidemic. Under this scenario, the indications for treating acutely infected persons may need to be modified by taking into account the potential balance between an increased number of infected persons as opposed to a decreased number of carriers of drug-resistant infections. Cost-benefit analyses of this nature are an intrinsic part of public health policy (33). In any case, we do not expect this process to be the main driving force responsible for recent trends in disease and drug resistance in the HIV-1 epidemic because the number of persons in the acute phase of HIV infection is much smaller than that in the chronic phase (1), and a large proportion of chronically infected patients received HAART at the beginning of the study period (36). Moreover, new treatment regimens, such as structured treatment interruptions and drug holidays, may have affected recent drug-resistance trends (36).
The most favorable outcome occurs if the increase in treatment efficacy brought about with HAART is the most important process determining recent HIV-1 trends. Now both the prevalence of the disease and of drug resistance in the population are decreasing, and therefore the decline in drug-resistance prevalence among the acutely infected is a positive sign (25,26,32). Results under this scenario underscore the importance of public health interventions directed toward increasing the number of persons receiving treatment (we must keep in mind that these results assume that the fraction of persons treated remains constant).
However, the uncertainty in the parameters does not allow us to readily distinguish between the 3 likely scenarios. Moreover, we cannot be sure that the trends observed correspond to a real decrease in FR, or are simply fluctuations due to stochastic or sampling phenomena around a monotonically increasing time trend (11,18,25,34,37–39). Little et al. and Grant et al. report trends consistent with the first scenario: the decrease in the fraction of persons infected with an NRTI-resistant strain is followed by an increase and subsequent decrease (11–13). Several authors report biphasic patterns of alternating trends (1,2,5,9,10,14), which are not necessarily correlated with an increase in non-B subtypes (8). Other studies report overall increasing (38), stable (40), and decreasing (3,4,6,7) trends in the proportion of persons recently infected with an NRTI drug-resistant HIV-1 strain. These studies, together with our results, highlight why disease surveillance must be increased, with additional data collection and analyses, to fully understand the present and future course of the HIV epidemic. In this regard, mathematical modeling can provide a crucial tool for the correct interpretation of epidemiologic data by identifying the processes responsible for generating observed time trends and characterizing their potential implications for public health programs.
Conclusion
Our mathematical analysis shows that the observed time trends of measurable quantities from particular subgroups of infected persons (such as primary drug resistance in monitored subpopulations) can correspond to different and unexpected time trends of variables of critical public health interest that are not measured directly in the general population. On the other hand, with the appropriate analyses, information on drug resistant strains can be used not only to guide treatment in individual patients, but also as epidemiologic sentinels to help devise public health solutions. Because changes in the relative value of 2 strains that vary in any of their life history traits (such as their ability to be transmitted, to be suppressed when in the presence of drug therapies, or to lose mutations that confer drug resistance) can show information on an epidemic's trends, the reasoning and methods we used in this study can be applied equally well to understand the epidemiology of any genetically variable microbe.
Acknowledgments
Special thanks go to Jason Barbour and 2 anonymous reviewers for their careful reading of the manuscript and their insightful suggestions. We also acknowledge Joyce Troiano, Holly Ganz, Eran Karmon, James O. Lloyd-Smith, and other members of Wayne M. Getz's laboratory for their help throughout this study.
This research was supported by NIH-NIDA grant no. POHC01000726 (M.S.S.), NIH-NIDA R01-D10135 grant (T.C.P. and W.M.G.) and the James S. McDonnell Foundation (W.M.G.).
Dr Sánchez is a postdoctoral fellow in the Department of Environmental Science, Policy and Management, University of California, Berkeley. Her main research interests include the theoretical study of infectious disease epidemics and applications to public health intervention, as well as the evolutionary implications of host-parasite associations.
References
Goudsmit J, Weverling GJ, van der Hoek L, de Ronde A, Miedema F, Coutinho RA, et al. Carrier rate of zidovudine-resistant HIV-1: the impact of failing therapy on transmission of resistant strains. AIDS. 2001;15:2293–301.
Yerly S, Vora S, Rizzardi P, Chave JP, Vernazza PL, Flepp M, et al. Acute HIV infection: impact on the spread of HIV and transmission of drug resistance. AIDS. 2001;15:2287–92.
de Mendoza C, del Romero J, Rodríguez C, Corral A,Soriano V. Decline in the rate of genotypic resistance to antiretroviral drugs in recent HIV seroconverters in Madrid. AIDS. 2002;16:1830–2.
Simon V, Vanderhoeven J, Hurley A, Ramratnam B, Louie M, Dawson K, et al. Evolving patterns of HIV-1 resistance to antiretroviral agents in newly infected individuals. AIDS. 2002;16:1511–9.
Ammaranond P, Cunningham P, Oelrichs R, Suzuki K, Harris C, Leas L, et al. No increase in protease resistance and a decrease in reverse transcriptase resistance mutations in primary HIV-1 infection: 1992–2001. AIDS. 2003;17:264–7.
Bezemer D, Jurriaans S, Prins M, van der Hoek L, Prins JM, de Wolf F, et al. Declining trend in transmission of drug-resistant HIV-1 in Amsterdam. AIDS. 2004;18:1571–7.
Routy JP, Machouf N, Edwardes MD, Brenner BG, Thomas R, Trottier B, et al. Factors associated with a decrease in the prevalence of drug resistance in newly HIV-1 infected individuals in Montreal. AIDS. 2004;18:2305–12.
de Mendoza C, Rodríguez C, Eiros JM, Colomina J, García F, Leiva P, et al. Antiretroviral recommendations may influence the rate of transmission of drug-resistant HIV type 1. Clin Infect Dis. 2005;41:227–32.
Boden D, Hurley A, Zhang L, Cao Y, Guo Y, Jones E, et al. HIV-1 drug resistance in newly infected individuals. JAMA. 1999;282:1135–41.
Harzic M, Pellegrin I, Deveau C, Chaix ML, Dubeaux B, Garrigue I, et al. Genotypic drug resistance during HIV-1 primary infection in France (1996–1999): Frequency and response to treatment. AIDS. 2002;16:793–6.
Little SJ, Holte S, Routy JP, Daar ES, Markowitz M, Collier AC, et al. Antiretroviral-drug resistance among patients recently infected with HIV. N Engl J Med. 2002;347:385–94.
Grant RM, Hecht FM, Warmerdam M, Liu L, Liegler T, Petropoulos CJ, et al. Time trends in primary HIV-1 drug resistance among recently infected persons. JAMA. 2002;288:181–8.
Grant RM, Liegler T, Spotts G, Hecht FM. Declining nucleoside reverse transcriptase inhibitor primary resistance in San Francisco, 2000–2003. In: XII International HIV Drug Resistance Workshop: Basic Principles and Clinical Implications, Los Cabos, Mexico, 2003.
Violin M, Velleca R, Cozzi-Lepri A, Riva C, Grossi PA, Carnevale G, et al. Prevalence of HIV-1 primary drug resistance in seroconverters of the ICONA cohort over the period 1996–2001. J Acquir Immune Defic Syndr. 2004;36:761–4.
Basu P. Aggressive HIV strain sets off dubious public health measure. Nat Med. 2005;11:360.
Briones C, Perez-Olmeda M, Rodriguez C, del Romero J, Hertogs K, Soriano V. Primary genotypic and phenotypic HIV-1 drug resistance in recent seroconverters in Madrid. J Acquir Immune Defic Syndr. 2001;26:145–50.
Holodniy M, Charlebois ED, Bangsberg DR, Zolopa AR, Schulte M, Moss AR. Prevalence of antiretroviral drug resistance in the HIV-1-infected urban indigent population in San Francisco: a representative study. Int J STD AIDS. 2004;15:543–51.
Leigh Brown AJ, Frost SD, Mathews WC, Dawson K, Hellmann NS, Daar ES, et al. Transmission fitness of drug-resistant human immunodeficiency virus and the prevalence of resistance in the antiretroviral-treated population. J Infect Dis. 2003;187:683–6.
de Mendoza C, Rodríguez C, Corral A, del Romero J, Gallego O, Soriano V. Evidence for differences in the sexual transmission efficiency of HIV strains with distinct drug resistance genotypes. Clin Infect Dis. 2004;39:1231–8.
Nowak MA, May RM. Virus dynamics: mathematical principles of immunology and virology. Oxford, UK: Oxford University Press; 2000.
Anderson RM, Gupta S, May RM. Potential of community-wide chemotherapy or immunotherapy to control the spread of HIV-1. Nature. 1991;350:356–9.
Garnett GP, Anderson RM. Antiviral therapy and the transmission dynamics of HIV-1. J Antimicrob Chemother. 1996;37:135–50.
Levin BR, Bull JJ, Stewart FM. The intrinsic rate of increase of HIV AIDS: Epidemiological and evolutionary implications. Math Biosci. 1996;132:69–96.
Zaric GS, Brandeau ML, Bayoumi AM, Owens DK. The effects of protease inhibitors on the spread of HIV and the development of drug-resistant HIV strains: a simulation study. Simulation. 1998;71:262–75.
Blower SM, Gershengorn HB, Grant RM. A tale of two futures: HIV and antiretroviral therapy in San Francisco. Science. 2000;287:650–4.
Blower S, Ma L, Farmer P, Koenig S. Predicting the impact of antiretrovirals in resource-poor settings: preventing HIV infections whilst controlling drug resistance. Curr Drug Targets Infect Disord. 2003;3:345–53.
Velasco-Hernández JX, Gershengorn HB, Blower SM. Could widespread use of combination antiretroviral therapy eradicate HIV epidemics Lancet Infect Dis. 2002;2:487–93.
Blower SM, Aschenbach AN, Gershengorn HB, Kahn JO. Predicting the unpredictable: transmission of drug-resistant HIV. Nat Med. 2001;7:1016–20.
Law MG, Prestage G, Grulich A, van de Ven P, Kippax S. Modelling the effect of combination antiretroviral treatments on HIV incidence. AIDS. 2001;15:1287–94.
Law MG, Prestage G, Grulich A, van de Ven P, Kippax S. Modelling HIV incidence in gay men: increased treatment, unsafe sex and sexually transmissible infections. AIDS. 2002;16:499–501.
Clements MS, Prestage G, Grulich A, van de Ven P, Kippax S, Law MG. Modeling trends in HIV incidence among homosexual men in Australia 1995–2006. J Acquir Immune Defic Syndr. 2004;35:401–6.
Dangerfield BC, Fang YX, Roberts CA. Model-based scenarios for the epidemiology of HIV/AIDS: the consequences of highly active antiretroviral therapy. System Dynamics Review. 2001;17:119–50.
Blower S, Volberding P. What can modeling tell us about the threat of antiviral drug resistance Curr Opin Infect Dis. 2002;15:609–14.
Blower SM, Aschenbach AN, Kahn JO. Predicting the transmission of drug-resistant HIV: comparing theory with data. Lancet Infect Dis. 2003;3:10–1.
Sánchez MS, Grant RM, Porco TC, Gross KL, Getz WM. A decrease in drug resistance levels of the HIV epidemic can be bad news. Bull Math Biol. 2005;67:761–82.
de Mendoza C, Martín-Carbonero L, Gallego O, Corral A, González-Lahoz J, Soriano V. Relationship between drug resistance mutations, plasma viremia, and CD4(+) T-cell counts in patients with chronic HIV infection. J Med Virol. 2005;76:1–6.
Duwe S, Brunn M, Altmann D, Hamouda O, Schmidt B, Walter H, et al. Frequency of genotypic and phenotypic drug-resistant HIV-1 among therapy-naive patients of the German seroconverter study. J Acquir Immune Defic Syndr. 2001;26:266–73.
UK Collaborative Group on Monitoring the Transmission of HIV Drug Resistance. Analysis of prevalence of HIV-1 drug resistance in primary infections in the United Kingdom. BMJ. 2001;322:1087–8.
Weinstock HS, Zaidi I, Heneine W, Bennett D, Garcia-Lerma JG, Douglas JM, et al. The epidemiology of antiretroviral drug resistance among drug-naive HIV-1-infected persons in 10 US cities. J Infect Dis. 2004;189:2174–80.
Chaix ML, Descamps D, Harzic M, Schneider W, Deveau C, Tamalet C, et al. Stable prevalence of genotypic drug resistance mutations but increase in non-B virus among patients with primary HIV-1 infection in France. AIDS. 2003;17:2635–43.(María S. Sánchez, Robert )
Gladstone Institute of Virology and Immunology, San Francisco, California, USA
California Department of Health Services, Berkeley, California, USA
Observed declines in drug resistance to nucleoside reverse transcriptase inhibitors among persons recently infected with HIV-1 in monitored subpopulations can be interpreted as a positive sign and lead public health officials to decrease efforts towards HIV prevention. By means of a mathematical model, we identified 3 processes that can account for the observed decline: increase in high-risk behavior, decrease in proportion of acutely infected persons whose conditions are treated, and change in treatment efficacy. These processes, singly or in combination, can lead to increases or decreases in disease and drug-resistance prevalence in the general population. We discuss the most appropriate public health response under each scenario and emphasize how further data collection and analyses are required to more reliably evaluate the observed time trends and the relative importance of forces shaping the epidemic. Our study highlights how drug resistance markers can be used as epidemiologic sentinels to devise public health solutions.
In recent years investigators have begun monitoring the HIV epidemic by reporting changes in the proportion of newly infected persons who are carrying an HIV-1 drug-resistant strain, i.e., the primary or acute resistant fraction. Several studies report decreases in this primary resistant fraction (1–7), including transient decreases (8–14). Unfortunately, none of these studies included precise longitudinal data on the exact number and type of infected persons or of the fraction of the total population that is screened for acute infection or resistance, and as we show here, making direct interpretations from data collected from a subset of the population can lead to erroneous conclusions. Given the potentially serious clinical implications of drug resistance for HIV-infected persons, public health officials and other authorities need to know whether the decline in drug resistance among acutely infected persons in the monitored subpopulations corresponds to a real decline in drug resistance in the general population and whether this effect is sustainable over time. The decline might be interpreted as a positive sign caused, for example, by less high-risk activity by HIV-positive persons infected with a drug-resistant variant. This explanation could lead public health officials to decrease their support for HIV surveillance and prevention programs targeted at impeding the spread of drug-resistant HIV strains such as drug-resistance testing or adherence counseling. The amount of resources that should be dedicated to drug resistance monitoring and reporting is a controversial issue in light of the recent isolation of a highly virulent multidrug-resistant strain in New York City (15). Here we show how drug-resistance data can offer not only clinical information regarding appropriate treatment regimens for individual patients, but also critical insights into an epidemic's course.
Treatment History of HIV-1 and Its Impact on High-Risk Behavior
Different types of drugs have been developed to fight HIV. Zidovudine (AZT), a nucleoside reverse transcriptase inhibitor (NRTI), was first administered in 1987, and until 1995, monotherapy or dual therapy with NRTIs were the only treatments available. The first protease inhibitor (PI), saquinavir, was approved for treatment in 1995, followed closely in 1996 by a nonnucleoside reverse transcriptase inhibitor (NNRTI), nevirapine. These new drugs generated a major change in the treatment strategy against HIV—highly active antiretroviral therapy (HAART)—that coincided with the start of the monitoring periods in several of the studies mentioned above (1995–1996). With HAART, at least 3 drugs are administered at the same time, which substantially reduces viral load and, compared to results of earlier regimens, increases the life expectancy of patients. These advantages follow because the mutations necessary to confer resistance to HAART are generated at a slower rate and are lost more rapidly than those conferring resistance to monotherapy or dual therapy. Moreover, viral strains resistant to HAART are not as efficient at completing their own life cycle (e.g., their replication rates are lower), they may generate less illness and lower proportion of deaths among infected persons, and the viral strains are less likely to be transmitted to other persons.
The primary resistance time trends observed for NRTIs do not match those observed for the other 2 drug types. For example, in North America some researchers report a decrease in the proportion of persons recently infected with a drug-sensitive HIV strain resistant to NRTIs; this decrease is followed by an increase and subsequent decrease (11–13). The pattern is particularly noticeable in the study by Grant et al. (12), in which NRTI genotypic resistance decreased from ≈30% in 1997 to 5% in 1999, rose to 20% in 2000, and fell to 15% in 2001. Little et al. also found this trend in primary resistance to NNRTIs (11). These 2 studies (11,12) also documented a steady increase in the proportion of persons newly infected with a virus strain resistant to PIs. Here we focused on identifying the likely forces responsible for the time trends exhibited by viral strains resistant to NRTIs. We did so because their time trends are expected to provide better insight into the long-term dynamics of the epidemic than strains resistant to PIs or NNRTIs, given that NRTIs have been administered to more HIV-infected persons and for a longer period than the other 2 types of drugs (16,17).
Treatment optimism after the initial successes of HAART likely affected the subsequent dynamics of HIV because these favorable treatment outcomes led some persons to increase their high-risk behavior. Later it became apparent that HAART does not completely eliminate HIV from an infected person or impede its transmission. Moreover, when HAART first became implemented, the best strategy was believed to be "hit hard, hit early," because the medical community was trying to limit the expansion of HIV within an infected person's body and ameliorate the gradual deterioration of the patient's immune system. However, HAART can have considerable negative side effects, which affect the functioning of the gastrointestinal system, renal system, pancreas, and liver and produce changes in blood count, allergies, lactic acidosis, and other problems. As a result, treatment began to be delayed to balance the following factors: 1) containing the viral load, 2) minimizing the risk of drug-resistant mutants developing by limiting the amount of treatment time, and 3) reducing negative side effects.
Modeling Drug Resistance in HIV
Empirical studies have shown that antiretroviral treatment (ARV) produces substantial changes in the viral dynamics at the within-host level that translate into substantial changes at the between-host level (8,18,19). Mathematical and computational models permit us to create simplified versions of complex realities that we can manipulate to further our understanding of their dynamic behavior. Consequently, numerous theoretical studies have investigated the impact of drug therapy on HIV dynamics at both levels (e.g., 20). Initial HIV treatment models (e.g., 21–23) addressed how ARVs might affect the infectiousness of treated persons, and the spread of HIV and its disease-induced deaths. The magnitude of the public health threat created by drug-resistant HIV strains was only recognized later. As a result, Zaric et al. presented a novel model that showed that adhering to treatment regimens would discourage the emergence of multidrug-resistant HIV strains in heterogeneous populations (24). Blower et al. developed a relatively simple but revealing deterministic compartmental framework (25) that has served as the reference point for most of the modeling studies subsequently done to investigate the effect of ARVs on disease incidence and prevalence, drug-resistance transmission and prevalence, AIDS death rate, and the potential to eradicate the HIV epidemic (26–31). Dangerfield et al. built a detailed HAART treatment model that accounts for persons in all 4 HIV stages, and the last is partitioned in early- and late-stage AIDS (32). They investigated the effects of HAART on HIV incidence and prevalence, assuming different average efficacious periods and assimilation times for HAART, different infectivity probabilities when receiving HAART, and different increases in the mean number of sexual partners.
Blower and Volberding reviewed mathematical studies used to understand the dynamics of a drug-resistant HIV epidemic, predict the incidence and prevalence of drug-resistant HIV strains, evaluate cost-benefit strategies, and assess the impact of public health policies (33). The general approach to these studies had been to construct a descriptive simplification of the epidemic by identifying critical categories and processes and to use this structure to make predictions, given a set of assumptions regarding the parameter values. In this regard, several studies have characterized the epidemic's trends as monotonic, including the fraction of new HIV infections that are drug resistant (e.g., 16,25,28,34). Goudsmit et al. conducted an analysis in which including changes in treatment rates explained the nonmonotonic trends of zidovudine resistance observed in a cohort of newly infected homosexual men enrolled in the Amsterdam Cohort Study (1).
In our previous study (35), we extended the basic modeling framework detailed in (25) to incorporate additional complexity, including 2–3 separate categories of acutely infected persons, depending on whether a person was infected with a drug-sensitive HIV strain, a strain resistant to monotherapy, or a strain resistant to triple-drug therapy. In doing so, we were able to distinguish among acutely infected persons, who are clinically and epidemiologically distinct from uninfected and chronically infected persons (e.g., we can consider them to engage more frequently in high-risk behavior (2) than chronically infected persons). We also counted these categories separately and tracked their temporal trends and better channeled the different categories through the model, according to the different processes acting on them (such as a decrease in the proportion of persons receiving treatment among those recently infected, an effect that did not occur in persons in the chronic phase). We also created 10 subcategories in the chronically infected stage to more accurately represent the progression of persons from the acute stage of infection to AIDS (35). Moreover, because our intent was to explain observed trends rather than to make future predictions, we adopted a specific approach that consisted in altering individually the value of each parameter during a given simulation (as was done in [1] with treatment rates of all persons), rather than running simulations with a set of fixed parameters and comparing outcomes across runs.
Primary Resistant Fraction
The most direct explanation for the decrease in the observed proportion of newly infected persons infected with NRTI-resistant HIV-1 strains is that resistance to NRTIs in the recently HIV-infected population is declining. Unfortunately, the data need to be further evaluated because the HIV infection status of every person in the general population has not been monitored. We therefore do not have absolute numbers for these time trends, but only the relative numbers obtained from monitored subpopulations, which consist of consenting persons enrolled in research programs at specific locations. To be eligible to participate in these programs, patients had to display symptoms typical of an acute HIV seroconversion syndrome or have recently engaged in risky activities that could have placed them at risk of contracting HIV. Accordingly, what has decreased is the fraction of drug-resistant carriers among the pool of recently infected HIV patients who are willing to participate in particular research programs and attend clinics involved in these studies. The time trends exhibited by the variable representing the actual counts of all the newly infected persons who are carriers for a resistant strain in the general population may or may not be a direct match to those of the monitored subpopulations.
These points are best illustrated by considering the fraction of recently infected persons who are carriers of a drug-resistant strain (primary resistant fraction, FR), defined mathematically as
,
where S is the number of persons initially infected with a drug-sensitive strain, and R is the number of persons initially infected with a drug-resistant strain. This fraction may decrease because R decreases (fewer newly infected persons have a strain that is drug resistant), or because S increases (more newly infected persons carry drug-sensitive strains). If both S and R increase or decrease by the same proportion, FR remains unchanged. As explained below, the benefit of using this variable's time trends to further our understanding of the past, present, and future of the HIV epidemic is that underlying alterations in the relative values of drug-sensitive and drug-resistant strains may arise from a variety of mechanisms with critically different epidemiologic outcomes.
To determine which processes could have caused the observed decrease in FR, we built a mathematical model of HIV transmission ; a more mathematically detailed explanation of our analysis can be found in our previous study (35). We then simulated the epidemic using this model and varied each of the parameters shown in Figure 1 (e.g., the average number of high-risk contacts in 1 year, the likelihood of transmitting HIV given a high-risk contact, the fraction of persons with acute or chronic HIV infection that are placed on treatment each year, the likelihood of generating or losing drug resistance in 1 year). Once we determined which processes can cause a decrease in the acute fraction infected with a drug-resistant strain, we evaluated whether the process had occurred in industrialized countries in recent years. If so, we could consider the process as a potential contributor to the observed trends. As a result, we identified 3 independent processes that caused a decrease in FR and were consistent with the history of the HIV-1 epidemic in industrialized countries from 1995 to 2001: 1) overall increase in risky behavior, 2) decrease in the fraction of individuals in the acute phase who are placed on treatment, and 3) increase in the efficacy of treatment. Goudsmit et al. also found that discontinuation of monotherapy with zidovudine in 1996 explained the observed drop in zidovudine resistance in patients newly infected with HIV in the Amsterdam Cohort Study (1).
Course of the HIV Epidemic
How do these 3 processes (increase in high-risk behavior, treatment delay, and greater treatment efficacy) impact the course of the HIV epidemic Are these effects comparable or are they sufficiently different such that the policy implications will vary according to which one we interpret to be the leading cause for the observed decline To address these questions, we determined the short- and long-term changes induced by these processes on 2 variables of critical public health importance for which we do not have reliable measurements: prevalence of disease (i.e., overall fraction of persons infected with HIV in the general population) and prevalence of drug resistance (i.e., overall fraction of persons infected with a drug-resistant HIV strain in the general population). By doing so, we use time trend changes in the relative prevalence values of 2 viral strains (i.e., the fraction FR defined above) to make inferences about changes in the absolute values of disease and resistance prevalence in the population (i.e., actual counts of infected persons and carriers of drug-resistant infections in the general population). Figures 2B and 2C show that the 3 factors we identified as causing a decrease in the primary resistant fraction are predicted to force different and permanent, long-term changes in disease and drug-resistance prevalence. A synopsis of our findings is provided in the Table.
Policy Implications and Comparison with the Data
Our results demonstrate how a decrease in the fraction of persons recently infected with a drug-resistant HIV strain can occur not only when the epidemiologic conditions improve (i.e., disease and drug-resistance prevalence in the population decrease), but also when the epidemic worsens (i.e., disease and drug-resistance prevalence increase). The 3 processes that can generate the decrease in primary resistance are not mutually exclusive, and probably all have contributed to the observed time course of FR. The challenge now is to identify which one has had the greatest effect on the recent trends of disease and drug-resistance prevalence in the HIV-1 epidemic.
If an increase in high-risk behavior has dominated HIV-1 epidemiology since the onset of HAART, then the decrease in primary resistance, counter to intuition, signals a worsening of the epidemic: a greater number of persons may have become infected, and a greater number of persons may be infected with a viral strain resistant to drug therapy. Other studies have obtained similar conclusions (e.g., 28,29,31,32). If this is the case, the public health response to the decrease in drug-resistance levels among the acutely infected should be to expand programs aimed at reducing high-risk behavior.
Determining the most appropriate public health response is difficult if a decrease in the fraction of acutely ill persons receiving treatment is the main driving force of the HIV epidemic. Under this scenario, the indications for treating acutely infected persons may need to be modified by taking into account the potential balance between an increased number of infected persons as opposed to a decreased number of carriers of drug-resistant infections. Cost-benefit analyses of this nature are an intrinsic part of public health policy (33). In any case, we do not expect this process to be the main driving force responsible for recent trends in disease and drug resistance in the HIV-1 epidemic because the number of persons in the acute phase of HIV infection is much smaller than that in the chronic phase (1), and a large proportion of chronically infected patients received HAART at the beginning of the study period (36). Moreover, new treatment regimens, such as structured treatment interruptions and drug holidays, may have affected recent drug-resistance trends (36).
The most favorable outcome occurs if the increase in treatment efficacy brought about with HAART is the most important process determining recent HIV-1 trends. Now both the prevalence of the disease and of drug resistance in the population are decreasing, and therefore the decline in drug-resistance prevalence among the acutely infected is a positive sign (25,26,32). Results under this scenario underscore the importance of public health interventions directed toward increasing the number of persons receiving treatment (we must keep in mind that these results assume that the fraction of persons treated remains constant).
However, the uncertainty in the parameters does not allow us to readily distinguish between the 3 likely scenarios. Moreover, we cannot be sure that the trends observed correspond to a real decrease in FR, or are simply fluctuations due to stochastic or sampling phenomena around a monotonically increasing time trend (11,18,25,34,37–39). Little et al. and Grant et al. report trends consistent with the first scenario: the decrease in the fraction of persons infected with an NRTI-resistant strain is followed by an increase and subsequent decrease (11–13). Several authors report biphasic patterns of alternating trends (1,2,5,9,10,14), which are not necessarily correlated with an increase in non-B subtypes (8). Other studies report overall increasing (38), stable (40), and decreasing (3,4,6,7) trends in the proportion of persons recently infected with an NRTI drug-resistant HIV-1 strain. These studies, together with our results, highlight why disease surveillance must be increased, with additional data collection and analyses, to fully understand the present and future course of the HIV epidemic. In this regard, mathematical modeling can provide a crucial tool for the correct interpretation of epidemiologic data by identifying the processes responsible for generating observed time trends and characterizing their potential implications for public health programs.
Conclusion
Our mathematical analysis shows that the observed time trends of measurable quantities from particular subgroups of infected persons (such as primary drug resistance in monitored subpopulations) can correspond to different and unexpected time trends of variables of critical public health interest that are not measured directly in the general population. On the other hand, with the appropriate analyses, information on drug resistant strains can be used not only to guide treatment in individual patients, but also as epidemiologic sentinels to help devise public health solutions. Because changes in the relative value of 2 strains that vary in any of their life history traits (such as their ability to be transmitted, to be suppressed when in the presence of drug therapies, or to lose mutations that confer drug resistance) can show information on an epidemic's trends, the reasoning and methods we used in this study can be applied equally well to understand the epidemiology of any genetically variable microbe.
Acknowledgments
Special thanks go to Jason Barbour and 2 anonymous reviewers for their careful reading of the manuscript and their insightful suggestions. We also acknowledge Joyce Troiano, Holly Ganz, Eran Karmon, James O. Lloyd-Smith, and other members of Wayne M. Getz's laboratory for their help throughout this study.
This research was supported by NIH-NIDA grant no. POHC01000726 (M.S.S.), NIH-NIDA R01-D10135 grant (T.C.P. and W.M.G.) and the James S. McDonnell Foundation (W.M.G.).
Dr Sánchez is a postdoctoral fellow in the Department of Environmental Science, Policy and Management, University of California, Berkeley. Her main research interests include the theoretical study of infectious disease epidemics and applications to public health intervention, as well as the evolutionary implications of host-parasite associations.
References
Goudsmit J, Weverling GJ, van der Hoek L, de Ronde A, Miedema F, Coutinho RA, et al. Carrier rate of zidovudine-resistant HIV-1: the impact of failing therapy on transmission of resistant strains. AIDS. 2001;15:2293–301.
Yerly S, Vora S, Rizzardi P, Chave JP, Vernazza PL, Flepp M, et al. Acute HIV infection: impact on the spread of HIV and transmission of drug resistance. AIDS. 2001;15:2287–92.
de Mendoza C, del Romero J, Rodríguez C, Corral A,Soriano V. Decline in the rate of genotypic resistance to antiretroviral drugs in recent HIV seroconverters in Madrid. AIDS. 2002;16:1830–2.
Simon V, Vanderhoeven J, Hurley A, Ramratnam B, Louie M, Dawson K, et al. Evolving patterns of HIV-1 resistance to antiretroviral agents in newly infected individuals. AIDS. 2002;16:1511–9.
Ammaranond P, Cunningham P, Oelrichs R, Suzuki K, Harris C, Leas L, et al. No increase in protease resistance and a decrease in reverse transcriptase resistance mutations in primary HIV-1 infection: 1992–2001. AIDS. 2003;17:264–7.
Bezemer D, Jurriaans S, Prins M, van der Hoek L, Prins JM, de Wolf F, et al. Declining trend in transmission of drug-resistant HIV-1 in Amsterdam. AIDS. 2004;18:1571–7.
Routy JP, Machouf N, Edwardes MD, Brenner BG, Thomas R, Trottier B, et al. Factors associated with a decrease in the prevalence of drug resistance in newly HIV-1 infected individuals in Montreal. AIDS. 2004;18:2305–12.
de Mendoza C, Rodríguez C, Eiros JM, Colomina J, García F, Leiva P, et al. Antiretroviral recommendations may influence the rate of transmission of drug-resistant HIV type 1. Clin Infect Dis. 2005;41:227–32.
Boden D, Hurley A, Zhang L, Cao Y, Guo Y, Jones E, et al. HIV-1 drug resistance in newly infected individuals. JAMA. 1999;282:1135–41.
Harzic M, Pellegrin I, Deveau C, Chaix ML, Dubeaux B, Garrigue I, et al. Genotypic drug resistance during HIV-1 primary infection in France (1996–1999): Frequency and response to treatment. AIDS. 2002;16:793–6.
Little SJ, Holte S, Routy JP, Daar ES, Markowitz M, Collier AC, et al. Antiretroviral-drug resistance among patients recently infected with HIV. N Engl J Med. 2002;347:385–94.
Grant RM, Hecht FM, Warmerdam M, Liu L, Liegler T, Petropoulos CJ, et al. Time trends in primary HIV-1 drug resistance among recently infected persons. JAMA. 2002;288:181–8.
Grant RM, Liegler T, Spotts G, Hecht FM. Declining nucleoside reverse transcriptase inhibitor primary resistance in San Francisco, 2000–2003. In: XII International HIV Drug Resistance Workshop: Basic Principles and Clinical Implications, Los Cabos, Mexico, 2003.
Violin M, Velleca R, Cozzi-Lepri A, Riva C, Grossi PA, Carnevale G, et al. Prevalence of HIV-1 primary drug resistance in seroconverters of the ICONA cohort over the period 1996–2001. J Acquir Immune Defic Syndr. 2004;36:761–4.
Basu P. Aggressive HIV strain sets off dubious public health measure. Nat Med. 2005;11:360.
Briones C, Perez-Olmeda M, Rodriguez C, del Romero J, Hertogs K, Soriano V. Primary genotypic and phenotypic HIV-1 drug resistance in recent seroconverters in Madrid. J Acquir Immune Defic Syndr. 2001;26:145–50.
Holodniy M, Charlebois ED, Bangsberg DR, Zolopa AR, Schulte M, Moss AR. Prevalence of antiretroviral drug resistance in the HIV-1-infected urban indigent population in San Francisco: a representative study. Int J STD AIDS. 2004;15:543–51.
Leigh Brown AJ, Frost SD, Mathews WC, Dawson K, Hellmann NS, Daar ES, et al. Transmission fitness of drug-resistant human immunodeficiency virus and the prevalence of resistance in the antiretroviral-treated population. J Infect Dis. 2003;187:683–6.
de Mendoza C, Rodríguez C, Corral A, del Romero J, Gallego O, Soriano V. Evidence for differences in the sexual transmission efficiency of HIV strains with distinct drug resistance genotypes. Clin Infect Dis. 2004;39:1231–8.
Nowak MA, May RM. Virus dynamics: mathematical principles of immunology and virology. Oxford, UK: Oxford University Press; 2000.
Anderson RM, Gupta S, May RM. Potential of community-wide chemotherapy or immunotherapy to control the spread of HIV-1. Nature. 1991;350:356–9.
Garnett GP, Anderson RM. Antiviral therapy and the transmission dynamics of HIV-1. J Antimicrob Chemother. 1996;37:135–50.
Levin BR, Bull JJ, Stewart FM. The intrinsic rate of increase of HIV AIDS: Epidemiological and evolutionary implications. Math Biosci. 1996;132:69–96.
Zaric GS, Brandeau ML, Bayoumi AM, Owens DK. The effects of protease inhibitors on the spread of HIV and the development of drug-resistant HIV strains: a simulation study. Simulation. 1998;71:262–75.
Blower SM, Gershengorn HB, Grant RM. A tale of two futures: HIV and antiretroviral therapy in San Francisco. Science. 2000;287:650–4.
Blower S, Ma L, Farmer P, Koenig S. Predicting the impact of antiretrovirals in resource-poor settings: preventing HIV infections whilst controlling drug resistance. Curr Drug Targets Infect Disord. 2003;3:345–53.
Velasco-Hernández JX, Gershengorn HB, Blower SM. Could widespread use of combination antiretroviral therapy eradicate HIV epidemics Lancet Infect Dis. 2002;2:487–93.
Blower SM, Aschenbach AN, Gershengorn HB, Kahn JO. Predicting the unpredictable: transmission of drug-resistant HIV. Nat Med. 2001;7:1016–20.
Law MG, Prestage G, Grulich A, van de Ven P, Kippax S. Modelling the effect of combination antiretroviral treatments on HIV incidence. AIDS. 2001;15:1287–94.
Law MG, Prestage G, Grulich A, van de Ven P, Kippax S. Modelling HIV incidence in gay men: increased treatment, unsafe sex and sexually transmissible infections. AIDS. 2002;16:499–501.
Clements MS, Prestage G, Grulich A, van de Ven P, Kippax S, Law MG. Modeling trends in HIV incidence among homosexual men in Australia 1995–2006. J Acquir Immune Defic Syndr. 2004;35:401–6.
Dangerfield BC, Fang YX, Roberts CA. Model-based scenarios for the epidemiology of HIV/AIDS: the consequences of highly active antiretroviral therapy. System Dynamics Review. 2001;17:119–50.
Blower S, Volberding P. What can modeling tell us about the threat of antiviral drug resistance Curr Opin Infect Dis. 2002;15:609–14.
Blower SM, Aschenbach AN, Kahn JO. Predicting the transmission of drug-resistant HIV: comparing theory with data. Lancet Infect Dis. 2003;3:10–1.
Sánchez MS, Grant RM, Porco TC, Gross KL, Getz WM. A decrease in drug resistance levels of the HIV epidemic can be bad news. Bull Math Biol. 2005;67:761–82.
de Mendoza C, Martín-Carbonero L, Gallego O, Corral A, González-Lahoz J, Soriano V. Relationship between drug resistance mutations, plasma viremia, and CD4(+) T-cell counts in patients with chronic HIV infection. J Med Virol. 2005;76:1–6.
Duwe S, Brunn M, Altmann D, Hamouda O, Schmidt B, Walter H, et al. Frequency of genotypic and phenotypic drug-resistant HIV-1 among therapy-naive patients of the German seroconverter study. J Acquir Immune Defic Syndr. 2001;26:266–73.
UK Collaborative Group on Monitoring the Transmission of HIV Drug Resistance. Analysis of prevalence of HIV-1 drug resistance in primary infections in the United Kingdom. BMJ. 2001;322:1087–8.
Weinstock HS, Zaidi I, Heneine W, Bennett D, Garcia-Lerma JG, Douglas JM, et al. The epidemiology of antiretroviral drug resistance among drug-naive HIV-1-infected persons in 10 US cities. J Infect Dis. 2004;189:2174–80.
Chaix ML, Descamps D, Harzic M, Schneider W, Deveau C, Tamalet C, et al. Stable prevalence of genotypic drug resistance mutations but increase in non-B virus among patients with primary HIV-1 infection in France. AIDS. 2003;17:2635–43.(María S. Sánchez, Robert )