Molecular epidemiology for human pathogens has been successfully used as a tool to identify clusters in order to assess the modalities of the spread of infection and to monitor the impact of therapy on the transmission of drug resistance [1–3]. These, in turn, may translate into prevention measures focused on populations at risk.
Transmission of HIV-1 infection is largely influenced by viraemic levels . Individuals with primary HIV-1 infection (PHI) harbour high viral load in both blood and genital secretions [5–8] and it has been postulated that spread of HIV infection is driven in large part by these individuals [9,10]. However, because of the difficulty in identifying individuals with acute infection, there are only anecdotal reports on secondary transmission in this clinical setting [11,12].
Suboptimal therapy only partially contains viral replication and selects for drug-resistant variants. Non-nucleoside reverse transcriptase inhibitors and protease inhibitors (PI) are more potent than previously used drugs, which might have an impact on transmission of drug-resistant variants. Another factor that may influence the rate of transmission of drug-resistant variants is the prevalence of infection with non-B subtypes, since these non-B subtypes originate from countries where antiretroviral drugs are not readily available [13–15]. In Western Europe and in the United States, the prevalence of drug-resistant variants in patients with recent HIV infection has been approximately 10% in recent years [16–22]. This relatively high transmission rate of drug-resistant variants is a major concern for the efficacy of antiretroviral therapy in the future and has led to the recommendation of resistance testing in recently infected individuals [23,24].
In order to evaluate the impact of PHI in the spread of HIV and to assess the rate of transmission of drug-resistant variants, we have performed phylogenetic analyses based on pol sequences and collected epidemiological data in 197 individuals with documented PHI in Switzerland and neighbouring France over four years (1996–1999).
The study included all individuals with documented PHI identified in six AIDS centres of University Hospitals in Switzerland (Basel, Geneva, Lausanne, Lugano, St Gallen and Zurich) and two AIDS Centres of Hospitals close to Geneva (Annemasse and Annecy, France). For each patient, the clinician in charge completed an epidemiological questionnaire labelled with a code number. In the total of 197 individuals infected between January 1996 and January 2000, PHI was documented by evolving HIV antibody response and/or symptoms consistent with an acute retroviral syndrome  within 3 months in 70% of individuals and by seroconversion within 12 months of presentation in 30%. Sequence analyses were performed on plasma samples collected at the first visit, before initiation of antiretroviral treatment.
The procedure of viral isolation, reverse transcription, amplification and sequencing were described earlier . Briefly, viral RNA was isolated from 100 μl plasma using Amplicor HCV Cobas reagents (Roche, Basel, Switzerland). Viral RNA was transcribed into complementary DNA (cDNA), which was subjected to nested polymerase chain reaction (PCR) for reverse transcriptase (RT) and protease genes. Direct double-stranded sequencing was performed on an automatic sequencer with Big Dye terminator kit (Applied Biosystems, Foster City, California, USA). Alignment of sequences was performed using CLUSTAL W version 1.7 . Phylogenetic analyses were carried out using the PHYLIP package . Neighbour-joining trees with 100 replicates were used in the bootstrap analysis , and values above 85 (85%) were considered as supporting the grouping. These clusters were confirmed using Fitch–Margoliash and least-squares analysis with 100 bootstrap replicates . HIV-1 subtypes were identified on combined RT plus protease sequences using SimPlot software, which calculates and plots the percentage identity of the query sequence to a panel of reference sequences .
Quantification of HIV-1 RNA levels were performed using Amplicor HIV Monitor (Roche) according to the manufacturer's instruction. For the analysis of viraemia over time for each patient in the Swiss HIV Cohort Study, the first viraemia value for each year was used.
CD3, CD4 and CD8 lymphocyte cell counts were determined by flow cytometry (Coulter EPICS IV, Basel, Switzerland) using fluorescein-labelled DAKO-T3, DAKO-T8, and R-phycoerythrin DAKO-CD4 (Dako, Glostrup, Denmark).
A total of 197 individuals with documented PHI were included. The mean age was 34 years (range, 18–71); 80% were male and 94% were Caucasians. The risk factors for HIV-1 transmission were homosexual contacts in 42%, heterosexual contacts in 42%, intravenous drug use (IDU) in 13% and unknown in 3%. Breakdown of these data per year of HIV-1 infection is reported in Table 1. At the time of sequence analysis, the mean HIV-1 RNA was 5.08 log10 copies/ml (range, 3.0–7.45) and the mean CD4 cell count was 593 × 106 cells/l (range, 100–1817).
Sequences of the RT and protease genes were recovered for 191 and 176 individuals, respectively. All sequences showed a full-length open reading frame. HIV-1 B subtypes were identified using SimPlot analysis in 143/193 (74%) individuals. Non-B subtypes were detected in 23% of individuals in 1996, 22% in 1997, 21% in 1998 and 35% in 1999 (Fig. 1). The percentages of non-B subtypes were 1.6% A, 3.6% C, 2% D, 1% F, 0.5% G, 5.8% J, 4.8% AE and 6.6% AG. Non-B subtypes were more frequently detected in female than in male patients (41% versus 21%). Among Caucasians, 22% were infected with non-B subtypes compared with 60% of non-Caucasians. Non-B subtypes were detected in 7.5% of individuals infected by homosexual contacts, 40% by heterosexual contacts and 26% by IDU. An increase in HIV-1 subtype diversity over time was observed in the heterosexual risk group. The J subtype, already observed in 1996 in IDU, was first detected in 1999 in the heterosexual contact risk group, suggesting transmission of this subtype from IDU to heterosexuals.
Phylogenetic analyses and epidemiological linkages
A phylogenetic tree constructed by using neighbour-joining method for the 193 available RT sequences revealed significant clustering (bootstrap value ≥ 85/100) for 58 (30%) individuals (data not shown). Sequences from these individuals were re-analyzed using the Fitch–Margoliash method. The resulting phylogenetic tree demonstrates significant clustering (bootstrap value ≥ 90/100) for 56 (29%) individuals (Fig. 2). Similar results were obtained using phylogenetic analyses of the protease gene (data not shown). Eighteen clusters of genetically related viruses, each including 2 to 11 PHI individuals were observed. Clusters were more frequently observed in IDU than in homosexuals and heterosexuals (48%, 30% and 26%, respectively).
Epidemiological linkages were formally documented for three individuals included in a cluster of six individuals infected in the St Gallen (SG) area. The index case was a female infected through rape. She transmitted the infection to her regular sexual partner and to at least one other individual through needle sharing at the time she presented symptoms of acute HIV infection. The six individuals shared drugs or had common friends. Another large cluster (HIV subtype J) of 11 individuals was observed in the western part of Switzerland. This cluster was subdivided into three subclusters (Fig. 2). Evidence of secondary infection through needle sharing was documented by contact tracing in two individuals.
The 16 other clusters included two to four individuals. These individuals lived in the same geographical areas (Annecy and Lausanne are about 60 km from Geneva), and had the same risk factor for HIV infection, except two (homosexual contact and IDU). The delay between the estimated dates of infection was less than 12 months for individuals included in the same cluster. Epidemiological linkages were documented for six couples (five heterosexual and one homosexual couple). Four couples (two non-Caucasian and two Caucasian couples) presented a similar history. In each case, the husband had a short stay (less than 2 months) in Africa, had sexual intercourse with commercial sex workers, developed acute retroviral syndrome shortly after returning in Switzerland and infected his wife within 10 days. All the female partners developed an acute retroviral syndrome within 3 to 6 weeks. In one case, the wife started highly active antiretroviral therapy (HAART) within 10 days after the return of her husband, before seroconversion or other evidence of HIV infection developed. Treatment was stopped after 1 month and acute HIV infection developed in the wife. Among the 56 PHI individuals with significant sequence clustering, contact tracing was documented in 17 individuals, demonstrating secondary HIV infection at the time of PHI. There were no epidemiological linkages between individuals with unrelated sequences.
Transmission of drug-resistant variants
Mutations previously described as associated with antiretroviral drug resistance  were observed in 17/193 (8.8%) individuals. The prevalence of genotypic drug-resistant variants among patients with PHI was 8.6% (PI resistant, 3.0%) in 1996, 14.6% (PI resistant, 7.7%) in 1997, 8.8% (PI resistant, 2.0%) in 1998 and 5.0% (PI resistant, 1.9%) in 1999. A similar trend was observed when considering only the 143 PHI individuals infected with HIV-1 B subtype (8.0% in 1996, 15.6% in 1997, 11.6% in 1998 and 7.5% in 1999). Eleven (5.8%) individuals harboured variants with mutations associated with resistance to zidovudine and or stavudine (M41L, D67N, K70R, L210W and T215Y/F). The lamivudine-resistant mutation M184V was detected in three (1.6%) individuals and mutations associated with resistance to non-nucleoside reverse transcriptase inhibitors (G190A and/or Y181C) were detected in two (1.1%) individuals. Major mutations associated with PI resistance were detected in 6/176 (3.4%) individuals: the V82A/F mutation was detected in four of these and the L90M mutation in two. In these six individuals, two to five of the minor mutations associated with PI resistance (L10I, K20R, M36I, I54V, L63P, A71T/V and I84V) were also detected. Amino acid insertions in the protease (35ED, 37NN) were observed in three individuals; epidemiological linkage was documented for two of them. Drug-resistant variants were identified more frequently in male than female (10.4% versus 2.6%), and more frequently in IDU (13%) and homosexuals (11.3%) than in heterosexuals (6.1%).
There was a decrease in the prevalence of the transmission of drug-resistant variants at the later stages of the study. The proportion of patients with undetectable viraemia in the Swiss HIV Cohort Study increased from less than 10% to more than 50% between 1996 and 1999 (Fig. 3).
It has long been suspected that individuals with PHI are more likely to transmit infection, owing to their high viraemia and their unawareness of their infection status [9,10]. In this study, we found that approximately one-third of individuals with PHI harbour viral variants genotypically related to other individuals with recent HIV transmission and this was confirmed by contact tracing in one-third of them. This is likely a minimal estimate since during the 4 years of the study period there have been approximately 3000 newly reported HIV-1 infections in Switzerland  and we have identified only 197 individuals at the time of PHI.
In this large and systematic study, a peak in transmission of antiretroviral drug-resistant variants, reaching nearly 15%, was observed in 1997, followed by decreased prevalence in 1998 (8.8%) and 1999 (5.0%). As in previous studies, genotypic resistance to RT inhibitors was more frequent than to PI [16–22]. The prevalence of transmission of drug-resistant variants was higher in IDU (13%) and in the homosexual group (11%) than in the heterosexual group (6%). The high transmission rate of drug-resistant variants in the first two risk groups is likely because of the initial spreading of HIV-1 infection and, therefore, earlier antiretroviral treatment within these two risk groups. The lower transmission rate observed in the heterosexual risk group is also partially explained by the high percentage of women infected with HIV-1 non-B subtypes (41%).
The observed recent decrease in the prevalence of the transmission of drug-resistant variants is unexpected. Several factors might be contributive. Viral load is the most important factor determining the probability of HIV transmission among untreated serodiscordant individuals  as well as for vertical transmission [33,34]. It is, therefore, likely that any decrease in viraemia induced by HAART decreases infectiousness. Transmission of drug-resistant variants probably originates in a population of individuals who continue to be viraemic despite HAART. Conversely, transmission of drug-resistant variants is unlikely if HAART is effective, or if new cases of HIV infection derive from infected individuals who have not been exposed to antiretroviral drugs, such as African immigrants. We provide evidence that both factors might be involved in Switzerland. The proportion of patients with undetectable viraemia in the Swiss HIV Cohort Study increased from less than 10% to more than 50% between 1996 and 1999 (Fig. 3). Moreover in the period to 1997, a large proportion of patients were on mono- or bitherapy with nucleoside analogue inhibitors of RT. In this period, 563 (72%) of these patients had virological failure, with viraemia > 400 copies/ml 6 months after treatment initiation. In 67% of these patients, a switch of nucleoside analogue inhibitor and addition of a PI and/or non-nucleoside inhibitor of RT was associated with a decrease in viraemia to < 400 copies/ml in 1997–1999. Overall, this resulted in a decrease of the population of infected patients susceptible to transmit drug-resistant variants. An additional factor contributing to the decrease in the transmission of drug-resistant variants is the steady increase of individuals infected with HIV-1 non-B subtypes. These non-B subtypes were observed at high prevalence, but not exclusively, in African immigrants and in women infected through heterosexual contacts.
New guidelines favour a less aggressive strategy for initiation of treatment, based more on CD4 cell counts than on viraemic levels, and treatment interruptions within or outside clinical trials are increasingly popular because of drug adverse events. This obviously leads to an increased number of infectious individuals. Measures to prevent HIV transmission should, therefore, be reinforced. Finally, if mutations conferring multiple drug resistance result in reduced viral fitness , the limitations in the replication capacity might also reduce the risk of transmission.
The data presented here underline two main points: first, it is not inevitability that transmission of antiretroviral drug resistance will increase, since preventive programmes and treatment availability may decrease it; second, individuals with acute HIV infection are frequently involved in the spread of HIV in a geographical area where newly diagnosed infections decreased over recent years. Hence, the identification of newly infected individuals and contact tracing should be considered as an important public health measure to reduce HIV transmission further.
We thank C. Gaille and W. Caveng for excellent technical help.
1. Godfrey-Faussett P, Sonnenberg P, Shearer SC. et al
. Tuberculosis control and molecular epidemiology in a South African gold-mining community. Lancet 2000, 356: 1066–1071.
2. Katzenstein TL, Jorgensen LB, Permin H. et al
. Nosocomial HIV-transmission in an outpatient clinic detected by epidemiological and phylogenetic analyses. J AIDS 1999, 13: 1737–1744.
3. Belec L, Si Mohamed A, Muller-Trutwin MC. et al
. Genetically related human immunodeficiency virus type 1 in three adults of a family with no identified risk factor for intrafamilial transmission. J Virol 1998, 72: 5831–5839.
4. Quinn TC, Wawer MJ, Sewankambo N. et al
. Viral load and heterosexual transmission of human immunodeficiency virus type 1.Rakai Project Study Group.
N Engl J Med 2000, 342: 921–929.
5. Daar ES, Moudgil T, Meyer RD, Ho DD. Transient high levels of viremia in patients with primary human immunodeficiency virus type 1 infection. N Engl J Med 1991, 324: 961–964.
6. Clark SJ, Saag MS, Decker WD. et al
. High titers of cytopathic virus in plasma of patients with symptomatic primary HIV-1 infection. N Engl J Med 1991, 324: 954–960.
7. Tindall B, Evans L, Cunningham P. et al
. Identification of HIV-1 in semen following primary HIV-1 infection. AIDS 1992, 6: 949–952.
8. Vernazza P, Perrin L, Cooper D et al.Quest Study. High genital shedding of HIV during primary HIV infection supports the nee for early recognition and intervention.Seventh European Conference on Clinical Aspects and Treatment of HIV Infection
. Lisbon, October 1999 [abstract 401].
9. Koopman JS, Jacquez JA, Welch GW. et al
. The role of early HIV infection in the spread of HIV through populations. J Acquir Immune Defic Syndr Hum Retrovirol 1997, 14: 249–258.
10. Jacquez JA, Koopman JS, Simon CP, Longini IM Jr. Role of the primary infection in epidemics of HIV infection in gay cohorts. J Acquir Immune Defic Syndr 1994, 7: 1169–1184.
11. Clumeck N, Taelman H, Hermans P, Piot P, Schoumacher M, De Wit S. A cluster of HIV infection among heterosexual people without apparent risk factors. N Engl J Med 1989, 321: 1460–1462.
12. Schockmel GA, Balavoine JF, Yerly S, Perrin L. Use of highly sensitive assays for the evaluation of post-exposure HIV prophylaxis. AIDS 1998, 12: 1726–1727.
13. Boni J, Pyra H, Gebhardt M. et al
. High frequency of non-B subtypes in newly diagnosed HIV-1 infections in Switzerland. J Acquir Immune Defic Syndr 1999, 22: 174–179.
14. Barin F, Courouce AM, Pillonel J, Buzelay L. Increasing diversity of HIV-1M serotypes in French blood donors over a 10-year period (1985–1995).Retrovirus Study Group of the French Society of Blood Transfusion.
AIDS 1997, 11: 1503–1508.
15. Couturier E, Damond F, Roques P. et al
. HIV-1 diversity in France, 1996–1998.The AC 11 laboratory network.
AIDS 2000, 14: 289–296.
16. Yerly S, Kaiser L, Race E, Bru JP, Clavel F, Perrin L. Transmission of antiretroviral-drug-resistant HIV-1 variants. Lancet 1999, 354: 729–733.
17. Little SJ, Daar ES, D'Aquila RT. et al
. Reduced antiretroviral drug susceptibility among patients with primary HIV infection. JAMA 1999, 282: 1142–1149.
18. Boden D, Hurley A, Zhang L. et al
. HIV-1 drug resistance in newly infected individuals. JAMA 1999, 282: 1135–1141.
19. Wegner SA, Brodine SK, Mascola JR. et al
. Prevalence of genotypic and phenotypic resistance to anti-retroviral drugs in a cohort of therapy-naive HIV-1 infected US military personnel. AIDS 2000, 14: 1009–1015.
20. Tamalet C, Pasquier C, Yahi N. et al
. Prevalence of drug resistant mutants and virological response to combination therapy in patients with primary HIV-1 infection. J Med Virol 2000, 61: 181–186.
21. Salomon H, Wainberg MA, Brenner B. et al
. Prevalence of HIV-1 resistant to antiretroviral drugs in 81 individuals newly infected by sexual contact or injecting drug use.Investigators of the Quebec Primary Infection Study.
AIDS 2000, 14: F17–F23.
22. Brodine SK, Shaffer RA, Starkey MJ. et al
. Drug resistance patterns, genetic subtypes, clinical features, and risk factors in military personnel with HIV-1 seroconversion. Ann Intern Med 1999, 131: 502–506.
23. Hirsch MS, Brun-Vezinet F, D'Aquila RT. et al
. Antiretroviral drug resistance testing in adult HIV-1 infection: recommendations of an International AIDS Society-USA Panel. JAMA 2000, 283: 2417–2426.
24. Therapeutic management of HIV-infected subjects. Report of the expert group under the direction of Prof. J. F. Delfraissy, 1999.
Presse Med 2000, 29: 25–28.
25. Kinloch-de Loës S, de Saussure P. et al
. Symptomatic primary infection due to human immunodeficiency virus type 1: Review of 31 cases. Clin Infect Dis 1993, 17: 59–65.
26. Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, positions-specific gap penalties and weight matrix choice. Nucl Acids Res 1994, 22: 4673–4680.
27. Felsenstein J. PHYLIP Inference Package Version 3.5.
Seattle: Department of Genetics, University of Washington; 1993.
28. Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol 1987, 4: 406–425.
29. Fitch WM, Margoliash E. Construction of phylogenetic trees. Science 1967, 155: 279–284.
30. Lole KS, Bollinger RC, Paranjape RS. et al
. Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol 1999, 73: 152–160.
31. Schinazi RF, Larder BA, Mellors J. Resistance table: mutations in retroviral genes associated with drug resistance: 2000–2001 update. Int Antiviral News 2000, 8: 65–91.
32. OFSP. Sida-information. Bull OFSP 2000, 5: 108.108.
33. Connor EM, Sperling RS, Gelber R. et al
. Reduction of maternal-infant transmission of human immunodeficiency virus type 1 with zidovudine treatment.Pediatric AIDS Clinical Trials Group Protocol 076 Study Group.
N Engl J Med 1994, 331: 1173–1180.
34. The European Collaborative Study. Maternal viral load and vertical transmission of HIV-1: an important factor but not the only one. AIDS 1999, 13: 1377–1385.
35. Kaufmann D, Munoz M, Bleiber G. et al
. Virological and immunological characteristics of HIV treatment failure. AIDS 2000, 14: 1767–1774.
The members of the Swiss HIV Cohort Study are R. Amiet, M. Battegay (Chairman of the Scientific Board), E. Bernasconi, H. Bucher, Ph. Bürgisser, M. Egger, P. Erb, W. Fierz, M. Flepp (Chairman of the Clinical and Laboratory Committee), P. Francioli (President of the SHCS, Centre Hospitalier Universitaire Vaudois, CH-1011-Lausanne), H. J. Furrer, M. Gorgievski, H. Günthard, P. Grob, B. Hirschel, C. Kind, Th. Klimkait, B. Ledergerber, U. Lauper, M. Opravil, F. Paccaud, G. Pantaleo, L. Perrin, W. Pichler, J.-C. Piffaretti, M. Rickenbach (Head of Data Centre), C. Rudin (Chairman of the Mother & Child Substudy), P. Sudre, V. Schiffer, J. Schupbach, A. Telenti, P. Vernazza and R. Weber.