To determine the success of a decade of widespread intervention with highly active antiretroviral therapy (HAART) on controlling the human immunodeficiency virus type 1 (HIV-1) epidemic, we analysed the transmission dynamics of HIV-1 over the past 25 years among men having sex with men (MSM) in the Netherlands.
The first AIDS cases in the Netherlands among MSM were identified in 1981 . HAART, consisting of a combination of drugs from three, and later four, different drug classes, became widely available from 1996 onwards. HAART dramatically reduced plasma and seminal viral load [2,3], resistance  and mortality rates [5,6]. As infectivity is shown to be strongly correlated to viral load [7,8], the widespread use of HAART might thus be expected to have reduced the incidence of HIV infections. Paradoxically, resurgent epidemics have been widely reported post-HAART [9,10]. Increases in risk behaviour [11–13] and syphilis and gonorrhoea diagnoses have also been documented in populations of MSM in several developed countries [13–15]. Earlier mathematical modelling studies have demonstrated that an increase in risk behaviour has the potential to counterbalance the beneficial effect of HAART [16–23].
In the present study, we aim to evaluate the separate impact of risk behaviour, HIV testing behaviour and HAART on the HIV epidemic in Dutch MSM by means of a mathematical model fitted to data recorded within several national databases, which provided extensive information on epidemic trends.
A mathematical model describing HIV transmission and HAART use among MSM in the Netherlands was constructed. The model we used described natural (untreated) disease progression, diagnosis and subsequent use of HAART. The basic structure of the model is illustrated in Fig. 1. The modelling strategy was tailored to the task of analyzing annual HIV and AIDS diagnosis time series, and specifically to tracking changes in per capita transmission rates. The most important factor in this respect is a simultaneous estimation of the prevalence of infectious individuals, weighted by their relative infectiousness, which depends on stage of infection and treatment status, and the incidence of new infections. Mathematical details and analyses of the model, including sensitivity analyses, hypothetical scenarios and predictions, and further data are available on request from the authors.
A method to increase realism in compartmental models is to include a unidirectional flow through several compartments, corresponding to an Erlang survival distribution By fitting such a distribution to data from 130 MSM seroconverters before the HAART era in the Amsterdam Cohort Studies , the maximum likelihood estimate corresponded to five compartments with mean stay in each of 1.89 years. Patients starting their disease progression first spend on average 0.24 year in an extra initial compartment that represents primary infection, and we equated the last stage of infection with AIDS, an approximation that seemed reasonable given the match with the estimated duration of high transmissibility that has been observed prior to death . The relative infectiousness of each stage is calculated from Wawer et al., with primary and last stages of infection being more infectious than the other four stages. As 97% of the MSM population is infected with subtype B, imported infections are all assumed to be in the primary infection stage entering the Netherlands after short holidays, but this assumption was not critical (unpublished analysis available on request).
Highly active antiretroviral therapy
In the HAART era, patients start HAART after being diagnosed, and they will be either on successful treatment (no detectable viral load) or experience therapy failure, with a viral rebound and infectiousness . During successful treatment, HAART is assumed to block both HIV transmission and disease progression [2–8,27,28], and viral blips are not taken into account . Patients experiencing treatment failure are assumed to have periods of apparent successful treatment before failure . After failure, patients go through the unidirectional stages of natural disease progression. We assume that there are three HAART treatment opportunities before patients fail completely and progress to death (representing the diversity of treatments available) . The HAART era started in 1995 with clinical trials and compassionate use, and the mass treatment programme started in 1996 and was fully implemented by 1998. The influence of pre-HAART therapy on HIV viral load and transmission before 1995 is neglected. People start HAART with a rate irrespective of their stage of infection (as multistage disease progression is included in the model, this assumption approximately reproduces the observed pattern of HAART initiation), but at the AIDS stage people are set on HAART immediately. Disease progression is unidirectional. Parameters on HAART use and failure were obtained from the ATHENA national observational cohort .
Transmission rate and risk behaviour
The standardized per infectious capita transmission rate β(t) is a time-varying function that measures the relative rate at which an HIV-positive infectious individual infects new individuals. It is standardized by setting it equal to 1.0 for untreated, undiagnosed individuals in the asymptomatic stage of infection during the first phase of the epidemic (1980–1983), so that all other values are measured relative to this. It is primarily intended as a measure of changes in risk behaviour that can be estimated in our study, and for convenience β(t) will be referred to as risk behaviour rate. β(t) is in fact a compound measure that is affected by changes in the partner change rate, by the rate and nature of risky sex acts within partnerships, by the effect of ‘saturation’ of the susceptible population (when new sexual partners are already previously infected) and by the effect of the changing prevalence of other sexually transmitted infections (STIs) in modulating HIV transmission.
Parameters in the model explicitly adjust for the effect of HAART in reducing infectiousness, for the increased infectiousness during primary and late (AIDS) stages of disease and for the effect of diagnosis in reducing risk behaviour. We assumed that MSM have a 50% reduction in risk behaviour after becoming aware of their seropositive status and implemented this into our model . These assumptions were all encoded as disease-stage-specific scaling parameters of risk behaviour rate β(t).
We define the reproduction number R(t) as the average number of people an infected person at time t would infect over his whole infectious lifespan if conditions remained the same as at time t[31,32]. It incorporates all factors including risk behaviour, effect of diagnosis and the effects of treatment with HAART in preventing infection. If the within-country R(t) is greater than 1, then the epidemic will grow exponentially driven by local transmission, and conversely if this number is less than 1, the epidemic will contract down to a number proportional to the number of imported cases . It is a key aim of public health interventions to avoid a locally driven epidemic and maintain R(t) below one. The state of the epidemic can be characterized by R(t) that can be calculated from the best fit parameters in the model.
We fitted our model simultaneously to the observed time series of annual new diagnoses  and annual new AIDS cases (see below) [14,32,33], which are constrained by the diagnosis rate and the risk behaviour rate β(t), and this made it feasible to estimate both these unknown parameters. Changing these independent parameters has different effects, which differently affect the goodness of fit of the model to the time series. Increasing risk behaviour increases both the number of diagnoses and AIDS cases, whereas increasing the diagnosis rate increases the number of diagnoses in the short term but leads to sustained long-term reductions in the number of diagnoses and AIDS cases.
The analysis was stratified into four distinct historical intervals: 1980–1983, the first AIDS cases were diagnosed ; 1984–1995, serological testing became available, increasing HIV awareness, introduction of first mono-antiretroviral and dual-antiretroviral therapies [3,6]; 1996–1999, early HAART era; and 2000–2004, current HAART era. The diagnosis rate during asymptomatic stages was estimated but was assumed to be zero during the first period (1980–1983). Diagnosis was assumed to be rapid (within 1 month) after AIDS, whereas zero during primary infection. The mean time to diagnosis, defined under conditions at time t, was calculated from the estimated diagnosis rate. The epidemic is assumed to have started with an import of cases in 1980. The model was solved numerically using Runge-Kutta 4 algorithm and was fitted to data by a custom maximum likelihood method. All analyses were performed in Berkeley Madonna, version 8.0.1 (http://www.berkeleymadonna.com).
Annual new AIDS cases
To account for the effect of HAART in preventing progression to AIDS, we used different data sets to simultaneously fit to the following: before 1997, the model is fitted to annual data on AIDS diagnosis among MSM and collected by Statistics Netherlands from the beginning of the HIV epidemic [14,34]; from 1996, the model is fitted to annual data of number of MSM getting HIV diagnosed while having AIDS in the ATHENA national observational cohort .
Annual new diagnoses
From 1984 the model was fitted to data on annual diagnoses per year among MSM in the ATHENA national observational cohort . Since 1998 all HIV patients in the Netherlands have been registered and monitored as part of the ATHENA national observational cohort. The year of first HIV diagnosis is recorded retrospectively at the point of registration into ATHENA. Patients who received HAART and died in the period 1996–1997 were included in the ATHENA database retrospectively. Although there is some uncertainty on the completeness of the retrospective inclusion, it is expected to have only minor bearing on our results. MSM who died before 1996 are not included in the ATHENA database. We explicitly accounted for this data truncation process in our model by implemented chance of surviving until 1996 for the respective stages of infection. In parallel, a prediction is made of the true (not truncated) curve of the number of new diagnoses (Fig. 2a). Data from 2005 are still incomplete, and they are thus not included in the current study.
Source of infection
By the start of 2005, 5516 MSM diagnosed with HIV had been included in the ATHENA observational cohort. Of these registered infections, 8% were reported to have acquired the infection while abroad and 62% from a partner within the Netherlands. Of those born in the Netherlands, 4% were infected abroad and of those born abroad, 41% were infected abroad. We assume that the remaining 30% of infections with an unknown country of infection are split according to these respective ratios respective to their country of birth. Thus, we estimate that overall 14% of diagnosed infections are imported.
For model verification, we compared the model number of prevalent cases with number of living HIV-positive MSM in the ATHENA database. Also, data on the predicted annual number of deaths with documented annual AIDS deaths in the Netherlands  were used for model outcome verification. These data contain information only on sex and no more specific information on risk group. Hence, AIDS deaths among MSM were predicted using the percentage of MSM among male ATHENA participants in 1996 as an estimate, predicting that about 2000 MSM had died of AIDS before 1996  (Fig. 3a). We also compared the estimated proportion of newly diagnosed patients in each disease stage as estimated by our best model fit, with data on CD4 cell count at diagnosis.
Figure 2a and b shows the model curves that fitted best to the observed time series of annual new diagnoses and AIDS diagnoses (data on AIDS at diagnosis not shown). Figure 2d shows the estimated absolute number of new infections per year in the Netherlands. This peaked in 1983 with 802 new infections, and in 2004 with 554 new infections.
Estimates for the risk behaviour rate β(t), the reproduction number R(t) and the mean time to diagnosis are shown in Fig. 2c and d. Over the initial period (1980–1983), the estimate for the reproduction number R(t) is 2.39 [95% CI (confidence interval) 2.17–2.76]. Between 1984 and 1995, the risk behaviour rate declined by 2.3-fold (95% CI 2.03–2.83), indicating large reductions in risk behaviour, and thereby reduced the reproduction number R(t) below one to 0.89 (95% CI 0.85-0.93), that is, just below the epidemic threshold.
After 1995, when HAART was introduced, the reproduction number declined yet further to 0.76 (95% CI 0.7–0.86), but the reduction was not as great as it could have been due to a 18% (95% CI 3–34%) increase of the risk behaviour rate, β(t). The risk behaviour rate is estimated to have increased yet further over the period 2000–2004 and returned to only 29% (95% CI 22–72%) below its value in the initial period 1980–1983. Reductions in the estimated mean time from infection to diagnosis [from 3.71 years (95% CI 3.49–3.97) in 1984–1995 to 2.90 years (95% CI. 2.84–3.03) in 2000–2004] with consecutive reductions in risk behaviour and widespread treatment with HAART resulted in the reproduction number being much lower than in the initial time period 1980-1983. Still, R(t) for the last time period 2000–2004 is estimated to be 1.04 (95% CI 0.98–1.09), near or above the critical epidemic threshold, and thus indicating that HIV may once again be spreading epidemically among MSM in the Netherlands.
From the best fit model, we estimated that 24% of all living HIV-positive MSM were unaware of their HIV-positive status at the start of 2005 and that they account for 90% of new infections. Without both the increase of the risk behaviour rate and the decrease of time to diagnosis, the reproduction number R(t) would have decreased by 24% from 0.89 to 0.68 due to the introduction of HAART. The risk behaviour rate would need to increase by 32% to offset this benefit, with 43% in order to offset the simultaneous benefits of the increase in testing behaviour and with 59% in order to get R(t) equal to one, that is, to revert to epidemic growth. An increase of 66% was measured to have occurred. On the basis of these model estimates, we conclude that HAART has played an important role in limiting transmission but that any gains made have been more than offset by increases in the risk behaviour rate. Had these increases not occurred in the HAART era, the reproduction number R(t) would have declined to 0.6, and the epidemic would have been in convinced decline.
We verified our predictions subjectively for consistency with approximated data on annual number of AIDS deaths in MSM (see Methods), and on the number of currently living diagnosed individuals in the national patient database ATHENA , shown in Fig. 3, and on the number of annual AIDS diagnoses after 1996 (Fig. 2b). We considered the quality of fit acceptable given that the model was not fitted to these data. A qualitative comparison of CD4 cell counts at diagnosis with model predictions in terms of disease stages shows similar trends.
In a sensitivity analysis, the results on the key outcomes β(t) and R(t) appear to be very robust to a wide range of model variants. In particular, model results were consistent when assumptions about the relative infectiousness of disease stages, effect of diagnosis on behaviour, and time from diagnosis to start of therapy were varied. In all model variants, R(t) for 2000–2004 is estimated to be near or above the critical threshold (R = 1), thus implying uncontrolled epidemic spread, with estimates of the current reproduction number ranging between 0.95 and 1.33, depending on the scenario (details available on request).
The joint effect of HAART and risk behaviour on HIV incidence has been previously studied using mathematical models and empirical data [16–21,35]. Although based on different assumptions, all these studies come to the same conclusion regarding the potential for an increase in risk behaviour to offset the benefits of HAART in reducing transmission. Our study provides new evidence that this has actually occurred and quantifies its magnitude and timing within a well studied population of MSM.
A key feature of our study is the existence of several national databases recording diagnoses of HIV infection and AIDS, and deaths, allowing the diagnosis rate to be estimated reliably by simultaneously fitting to these time series within a robust inference framework. We were thus able to confirm that there has indeed been a recent increase in the diagnosis rate, reflecting a more frequent testing as was reported recently, but this was not sufficient to explain the recent increases in the number of people newly diagnosed. Rather, the recent increase in the number of new diagnoses reflects a substantial increase in transmission. Our estimates were corroborated by changing trends in CD4 cell count at diagnosis, where a recent increase in the proportion of newly diagnosed individuals with high CD4 cell counts is apparent.
Testing rates are low in the Netherlands when compared with other developed countries [36,37], and the potential of intervention by frequent testing with the rapid test is not yet fully explored . Our model, however, suggests that the only way to reverse epidemic spread, and get R well below one, is to reduce the risk behaviour rate from current levels. The potential effects of routine use of new diagnostic methods that target primary HIV infection were not explored here and should be explored in future models .
The most likely factor driving changes in the risk behaviour rate parameter β(t) is changing the sexual risk behaviour, both within partnerships and in partner change rates , though related factors such as other STIs acting to enhance transmission, saturation of the susceptible population or even evolution of infectivity could also play a role. Our analysis made it possible to compare the relative changes over time in risk behaviour rate between infectious and negative MSM, the ‘hidden’ information that cannot be measured by survey data, and our results indicate that whatever measures individuals are taking to ‘serosort’  are not proving effective at the population level and have not offset epidemic spread.
The introduction of HAART was accompanied by a decrease in the percentage of resistant strains among new infections [33,41]. However, the recent increase in annual new infections could in turn result in an increasing absolute number of resistant infections .
The widespread use of HAART has led to large reductions in AIDS morbidity and mortality (Figs. 2 and 3). Sustaining these reductions into the future will require either further improvements in treatment efficacy or a response to limit resurgent epidemic spread.
In conclusion, there is an increase in HIV transmission among MSM in the Netherlands, in spite of earlier diagnosis and subsequent effective treatment. The most effective intervention is to bring risk behaviour back to pre-HAART levels.
D.B. was supported by grant 7014 from AIDS Fund Netherlands and by a travel grant from NWO (Netherlands Organisation for Scientific Research). C.F. is funded by the Royal Society. The funders had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript.
1. Prummel MF, ten Berge RJ, Barrowclough H, Cejka V. Kaposi's sarcoma and fatal opportunistic infections in a homosexual man with immunodeficiency. Ned Tijdschr Geneeskd 1983; 127:820–824.
2. Gupta P, Mellors J, Kingsley L, Riddler S, Singh MK, Schreiber S, et al
. High viral load in semen of human immunodeficiency virus type 1-infected men at all stages of disease and its reduction by therapy with protease and nonnucleoside reverse transcriptase inhibitors. J Virol 1997; 71:6271–6275.
3. 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–2301.
4. Vandamme AM, Van Laethem K, De Clercq E. Managing resistance to anti-HIV drugs: an important consideration for effective disease management. Drugs 1999; 57:337–361.
5. Hammer SM, Katzenstein DA, Hughes MD, Gundacker H, Schooley RT, Haubrich RH, et al
. A trial comparing nucleoside monotherapy with combination therapy in HIV-infected adults with CD4 cell counts from 200 to 500 per cubic millimeter. N Engl J Med 1996; 335:1081–1090.
6. Smit C, Geskus R, Uitenbroek D, Mulder D, van den Hoek A, Coutinho RA, et al
. Declining AIDS mortality in Amsterdam: contributions of declining HIV incidence and effective therapy. Epidemiology 2004; 15:536–542.
7. Quinn TC, Wawer MJ, Sewankambo N, Serwadda D, Li CJ, Wabwire-Mangen F, et al
. Viral load and heterosexual transmission of human immunodeficiency virus type 1. N Engl J Med 2000; 342:921–929.
8. Fideli US, Allen SA, Musonda R, Trask S, Hahn BH, Weiss H, et al
. Virologic and immunologic determinants of heterosexual transmission of human immunodeficiency virus type 1 in Africa. AIDS Res Hum Retroviruses 2001; 17:901–910.
9. White PJ, Ward H, Garnett GP. Is HIV out of control in the UK? An example of analysing patterns of HIV spreading using incidence-to-prevalence ratios. AIDS 2006; 20:1898–1901.
10. UNAIDS. Report on the global AIDS epidemic, executive summary UNAIDS/06.29E.
Geneva: Joint United Nations Programme on HIV/AIDS (UNAIDS); 2006. http://www.unaids.org/en/HIV_data/2006GlobalReport/default.asp
11. Stolte IG, Dukers NHTM. Response to ‘High-risk sexual behaviour
increases among London gay men between 1998 and 2001: what is the role of HIV optimism?’. AIDS 2003; 17:2011–2012.
12. Stolte IG, Dukers NHTM, Geskus RB, Coutinho RA, De Wit JBR. Homosexual men
change to risky sex when perceiving less threat of HIV/AIDS since availability of highly active antiretroviral therapy
: a longitudinal study. AIDS 2004; 18:303–309.
13. Truong HHM, Kellogg T, Klausner JD, Katz MH, Dilley J, Knapper K, et al
. Increases in sexually transmitted infections and sexual risk behaviour without a concurrent increase in HIV incidence among men who have sex with men in San Francisco: a suggestion of HIV serosorting? Sex Transm Infect 2006; 82:461–466.
14. van de Laar MJW, de Boer IM, Koedijk FDH, Op de Coul ELM. HIV and Sexually Transmitted Infections in the Netherlands in 2004, RIVM report 441100022/2005
. Bilthoven: National Institute for Public Health and the Environment; 2005. http://www.rivm.nl/Rliotheek/rapporten/rapporten2005.html
15. van der Bij AK, Stolte IG, Coutinho RA, Dukers NH. Increase of sexually transmitted infections, but not HIV, among young homosexual men
in Amsterdam: are STIs still reliable markers for HIV transmission? Sex Transm Infect 2005; 81:34–37.
16. Blower SM, Gershengorn HB, Grant RM. A tale of two futures: HIV and antiretroviral therapy
in San Francisco. Science 2000; 287:650–654.
17. Nagelkerke NJD, Jha P, de Vlas SJ, Korenromp EL, Moses S, Blanchard JF, et al
. Modelling HIV/AIDS epidemics in Botswana and India: impact of interventions to prevent transmission. Bull World Health Organ 2002; 80:89–96.
18. Velasco-Hernandez JX, Gershengorn HB, Blower SM. Could widespread use of combination antiretroviral therapy
eradicate HIV epidemics? Lancet Infect DisV 2 2002:487–493.
19. 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–406.
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–1294.
21. Xiridou M, Geskus R, de Wit J, Coutinho R, Kretzschmar M. The contribution of steady and casual partnerships to the incidence of HIV infection among homosexual men
in Amsterdam. AIDS 2003; 17:1029–1038.
22. Baggaley RF, Ferguson NM, Garnett GP. The epidemiological impact of antiretroviral use predicted by mathematical models
: a review. Emerg Themes Epidemiol 2005; 2:9.
23. Hosseinipour M, Cohen MS, Vernazza PL, Kashuba ADM. Can antiretroviral therapy
be used to prevent sexual transmission of human immunodeficiency virus type 1? Clin Infect Dis 2002; 34:1391–1395.
24. van Griensven GJ, Tielman RA, Goudsmit J, Van der Noordaa J, de Wolf F, de Vroome EM, et al
. Risk factors and prevalence of HIV antibodies in homosexual men
in the Netherlands. Am J Epidemiol 1987; 125:1048–1057.
25. Wawer MJ, Gray RH, Sewankambo NK, Serwadda D, Li XB, Laeyendecker O, et al
. Rates of HIV-1 transmission per coital act, by stage of HIV-1 infection, in Rakai, Uganda. J Infect Dis 2005; 191:1403–1409.
26. Gras L, van Sighem A, van Valkengoed I, Zaheri S, de Wolf F. Monitoring of human immunodeficiency virus (HIV) infection in the Netherlands, Chapter 10, HIV Monitoring Foundation
. Amsterdam: Stichting HIV Monitoring; 2003. http://www.hiv-monitoring.nl
27. van Sighem AI, van de Wiel MA, Ghani AC, Jambroes M, Reiss P, Gyssens IC, et al
. Mortality and progression to AIDS after starting highly active antiretroviral therapy
. AIDS 2003; 17:2227–2236.
28. Barroso PF, Schechter M, Gupta P, Bressan C, Bomfim A, Harrison LH. Adherence to antiretroviral therapy
and persistence of HIV RNA in semen. J Acquir Immune Defic Syndr 2003; 32:435–440.
29. Havlir DV, Bassett R, Levitan D, Gilbert P, Tebas P, Collier AC, et al
. Prevalence and predictive value of intermittent viremia with combination hiv therapy. JAMA 2001; 286:171–179.
30. Marks G, Crepaz N, Senterfitt JW, Janssen RS. Meta-analysis of high-risk sexual behavior in persons aware and unaware they are infected with HIV in the United States: implications for HIV prevention programs. J Acquir Immune Defic Syndr 2005; 39:446–453.
31. Anderson RM, May RM. Infectious diseases of humans: dynamics and control
. Oxford, UK: Oxford University Press; 1991. ISBN 0 19 854040-X.
32. Fraser C. Estimating individual and household reproduction numbers in an emerging epidemic. PLoS ONE 2007; 2:e758.
33. Gras L, van Sighem A, Smit C, Zaheri S, de Wolf F. Monitoring of human immunodeficiency virus (HIV) infection in the Netherlands, Chapter 6, 9&10
. Amsterdam: HIV Monitoring Foundation; 2006. http://www.hiv-monitoring.nl
34. Postma MJ, Jager JC, Dijkgraaf MG, Borleffs JC, Tolley K, Leidl RM. AIDS scenarios for The Netherlands: the economic impact on hospitals. Health Policy 1995; 31:127–150.
35. Katz MH, Schwarcz SK, Kellogg TA, Klausner JD, Dilley JW, Gibson S, et al
. Impact of highly active antiretroviral treatment on HIV seroincidence among men who have sex with men: San Francisco. Am J Public Health 2002; 92:388–394.
36. Dukers NH, Fennema HS, van der Snoek EM, Krol A, Geskus RB, Pospiech M, et al
. HIV incidence and HIV testing behavior in men who have sex with men: using three incidence sources, The Netherlands, 1984–2005. AIDS 2007; 21:491–499.
37. Stolte IG, de Wit JB, Kolader ME, Fennema HS, Coutinho RA, Dukers NH. Low HIV-testing rates among younger high-risk homosexual men
in Amsterdam. Sex Transm Infect 2007; 83:387–391.
38. Hutson S. Is home test for HIV a good idea? New Sci 2005; 188:16.
39. Pilcher CD, Fiscus SA, Nguyen TQ, Foust E, Wolf L, Williams D, et al
. Detection of acute infections during HIV testing in North Carolina. N Engl J Med 2005; 352:1873–1883.
40. van der Bij AK, Kolader ME, de Vries HJ, Prins M, Coutinho RA, Dukers NH. Condom use rather than serosorting explains differences in HIV incidence among men who have sex with men. J Acquir Immune Defic Syndr 2007; 45:574–580.
41. 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–1577.
42. Sanchez MS, Grant RM, Porco TC, Getz WM. HIV drug-resistant strains as epidemiologic sentinels. Emerg Infect Dis 2006; 12:191–197.