Access to combination antiretroviral therapy (ART) for HIV-1-infected persons and antiretroviral (ARV) prophylaxis for prevention of mother-to-child transmission of HIV-1 (PMTCT) in sub-Saharan Africa has greatly expanded during the past 5 years.1 In resource-limited settings where access to routine HIV-1 viral load monitoring is lacking and where the unregulated use of ARV drugs may be common, the selection of drug-resistant HIV-1 variants2 and their subsequent transmission to newly infected individuals3-6 is of particular concern, especially since second-line treatment options are limited. Few studies have assessed the mutational patterns associated with HIV-1 drug resistance among pre-treatment populations in sub-Saharan Africa in which drug pressure is increasing after ARV rollout.7-12
The government of Zambia, a country in southern Africa, which is among the countries worst affected by the HIV-1 pandemic, initiated a comprehensive HIV-1 care and treatment program with support from international agencies.13 By the end of 2007, nation-wide ART coverage was 46% of those in need.1 Standard first-line ART regimens combine a dual nucleoside/nucleotide reverse transcriptase inhibitor (NRTI) backbone with a nonnucleoside reverse transcriptase inhibitor (NNRTI).14
The objective of this study was to evaluate the mutational patterns and factors associated with baseline drug resistance in HIV-1-infected individuals present at time of initiating first-line ART in the geographical setting of Lusaka, Zambia, where ART first became available in the country.
Study Population and Design
The PharmAccess African Studies to Evaluate Resistance Monitoring Study is a multicenter prospective observational cohort of HIV-1-infected patients who receive ART in routine circumstances at 13 clinical sites in 6 African countries.15 We conducted a cross-sectional analysis including 3 clinical sites in Lusaka, Zambia (Fig. 1) as follows: Lusaka Trust Hospital, a private general hospital (Woodlands area); KARA Clinic, a free nongovernment sector clinic (city center); and, Coptic Hospital, a free faith-based general hospital (Manda Hill area). The 3 sites have provided HIV-1 care and treatment since 1997, 2004, and 2006, respectively. The Academic Medical Center Institutional Review Board and the University of Zambia Research Ethics Committee approved all study procedures. Confirmed HIV-1 seropositive individuals aged ≥18 years who were eligible to initiate first-line ART as defined by national guidelines (ie, advanced immunodeficiency as defined by CD4 count <200 cells/μL or advanced disease according to the World Health Organization (WHO) clinical stages)14 were consecutively enrolled. All participants provided written informed consent for use of routinely collected demographic, clinical, and laboratory data and additional phlebotomy for assessment of HIV-1 RNA and genotypic resistance. Exclusion criteria were pregnancy at study screening and re-initiation on first-line ART less than 30 days after stopping previous first-line ART. Re-initiation on first-line ART more than 30 days after stopping previous first-line ART and/or any previous use of ARV prophylaxis or nonsuppressive mono/dual therapy were not exclusion criteria (ARV-experienced group).
HIV-1 RNA determination was performed on EDTA-anticoagulated plasma using the NucliSens EasyQ real-time assay version 1.2 (bioMérieux, Lyon, France). Population-based sequencing was performed on all plasma specimens which had HIV-1 RNA >1000 copies per milliliter using an in-house method.16 Briefly, HIV-1 RNA was extracted from 200 μL of plasma using the automated Roche MagNa Pure LC analyzer and the MagNA Pure LC Total Nucleic Acid Isolation Kit (Roche, Germany). Genotyping encompassed protease and codons 1-230 of reverse-transcriptase, using an in-house sequencing method with an ABI 3730 Genetic Analyzer (Applied Biosystems, Foster City, CA). Sequences were assembled and manually edited using Sequencher version 4.5 software (Genecodes, Ann Arbor, MI). GenBank accession numbers: HM119603-HM120150.
Genotypic Resistance Analysis and Subtyping
Baseline resistance was defined as the presence of ≥1 Drug resistance-associated mutation (DRM) according to the WHO 2009 Surveillance Drug Resistance Mutation list,17 using the Stanford Calibrated Population Resistance analysis tool (version 4.1 beta, available at http://hivdb.stanford.edu/). HIV-1 subtypes were determined using the REGA HIV-1 subtyping algorithm (version 2.0, available at http://www.bioafrica.net/subtypetool/html/)18 and additional phylogenetic analysis using neighbor-joining method if required.
Univariate and multivariate analyses were used to determine factors associated with drug resistance using logistic regression and expressed as odds ratios (ORs) (95% confidence interval, CI) and P values (P < 0.05 statistically significant). Prevalence values were calculated with a CI based on the binomial distribution. Categorical data were compared using χ2 test. Continuous data were investigated using Student t test. All analyses were performed using Stata version 10 (StataCorp LP, TX).
Between March 2007 and September 2008, a total of 839 adult men and nonpregnant women were recorded to initiate first-line ART at the 3 sites. Screening efforts resulted in the enrolment in the study of 584 individuals who met eligibility criteria and provided consent (ie, 70% recruitment rate). A valid baseline HIV-1 RNA result was available for 576 (98.6%) participants. HIV-1 RNA was >1000 copies per milliliter in 556 (96.5%); of those, sequence analysis was successful in 548 (98.6%). Table 1 summarizes the demographic characteristics of all 584 participants. Females comprised 54.8% (n = 320) and were younger than males [36.1 (SD: 8.8) versus 40.0 (SD 8.8) years, P < 0.0001]. All participants were native Zambians. 321 (55.0%) participants had advanced stage disease (ie, WHO stage III or IV). Median CD4 count was 132 cells per microliter (interquartile range, IQR: 69-203). Mean HIV-1 RNA was 5.0 (SD: 0.9) log10 copies/mL. Twenty-seven (4.6%) participants had previous ARV experience, either as (highly active) ART (n = 14), single-dose nevirapine for PMTCT (n = 5), combination therapy for PMTCT (n = 1), or unspecified (n = 7). Patient characteristics neither differed between ARV-naive versus ARV-experienced participants (Table 1) nor between sites (data not shown).
Frequencies of Subtypes and Drug Resistance-Associated Mutations
Subtype C was identified in 98.0% of participants (537 of 548). Other subtypes and circulating recombinant forms (CRFs) were subtype A1 (0.5%, 3 of 548), CRF02_AG (0.5%, 3 of 548), G (0.4%, 2 of 548), CRF09_cpx (0.4%, 2 of 548), and D (0.2%, 1 of 548). The overall prevalence of resistance was 5.7% (31 of 548 participants; CI: 3.9% to 7.9%); the prevalence of DRMs associated with NRTIs, NNRTIs, and protease inhibitors (PIs) was 1.1% (6 of 548), 4.0% (22 of 548), and 1.1% (6 of 548), respectively (Fig. 1). Among ARV-naive participants, the prevalence of resistance was 5.2% (27 of 523; CI: 3.4% to 7.4%); the prevalence of DRMs associated with NRTIs, NNRTIs, and PIs was 1.0% (5 of 523), 3.6% (19 of 523), and 1.1% (6 of 523), respectively. Among ARV-experienced participants, the prevalence of resistance was 16.0% (4 of 25 participants; CI: 4.5% to 36.1%); the prevalence of DRMs associated with NRTIs, NNRTIs, and PIs was 4.0% (1 of 25), 12.0% (3 of 25), and 0.0% (0 of 25), respectively. Detected DRMs were Y181C (1.6%, 9 of 548), K103N (1.3%, 7 of 548), K103S (0.7%, 4 of 548), G190A (0.7%, 4 of 548), L100I (0.4%, 2 of 548), K101E (0.4%, 2 of 548), M184V (0.4%, 2 of 548), V106M (0.2%, 1 of 548), Y188C (0.2%, 1 of 548), G190S (0.2%, 1 of 548), K65R (0.2%, 1 of 548), T69D (0.2%, 1 of 548), K70R (0.2%, 1 of 548), K70E (0.2%, 1 of 548), L74I (0.2%, 1 of 548), V75S (0.2%, 1 of 548), V75T (0.2%, 1 of 548), and K219E (0.2%, 1 of 548) in reverse transcriptase and L90M (0.7%, 4 of 548), I85V (0.2%, 1 of 548), and I50L (0.2%, 1 of 548) in protease. DRM frequencies did not differ significantly across sites (data not shown). Table 2 provides an overview of demographic and virologic characteristics of the 31 participants who harbored ≥1 DRM. Dual-class resistance to NRTIs and NNRTIs was detected in 0.6% (3 of 523) of ARV-naive participants, and no triple-class resistance was observed. One ARV-naive participant harbored M184V plus 2 thymidine analogue mutations (TAMs) by the TAMII pathway (K70R and K219E).
Factors Associated With Baseline Drug Resistance
In univariate analysis, participants with versus without resistance did not differ for age, sex, WHO clinical stage, median serum hemoglobin level, median CD4 count, mean HIV-1 RNA, and subtype, whereas previous ARV experience was more frequent among participants with drug resistance compared with those without resistance (12.9% versus 4.1%, P = 0.031) (Table 3). Multiple logistic regression analysis of patient characteristics demonstrated that previous ARV use (versus ARV-naive status; OR: 4.32, CI: 1.34 to 14.0, P = 0.015) and hemoglobin level (highest tertile versus lowest tertile; OR: 2.74, CI: 1.09 to 6.89, P = 0.033) were independently associated with the presence of baseline resistance (Table 3).
In this study, we investigated patterns of drug-resistant HIV-1 present at time of initiating first-line ART among 584 predominantly HIV-1 subtype C-infected patients in Lusaka, Zambia, enrolled in the The PharmAccess African Studies to Evaluate Resistance Monitoring study during 2007-2008. Six (CI 3.9 to 7.9) percent of participants were found to harbor ≥1 DRM, based on population genotyping. Baseline resistance may reflect the combined effect of drug-resistant strains transmitted during infection and acquired during previous ARV exposure. The majority (95.4%) of participants were reported to be ARV naive at baseline; compared with ARV-naive status, resistance was more frequent after prior use of ART or PMTCT (16.0% versus 5.2%, P = 0.022). NNRTI-associated DRMs were observed in the highest frequency (4.0%), whereas dual-class or triple-class resistance was rarely observed. An important strength of the study was its large multisite patient sample, representative of a variety of clinic populations within the geographic setting of Lusaka.
Primary HIV-1 drug resistance in Zambia had only been assessed in 1 small study before ARV drugs became widely available, reporting no major drug resistance-associated mutations (DRMs) among 28 ARV-naive persons.19 Population surveys7-12 and mathematical models20 to date have reported low levels of drug-resistant infections in African populations with increasing selective ARV drug pressure, but recent preliminary reports have suggested transmission of resistant strains directly after national ARV rollout programs.21,22 The predominance of NNRTI-associated DRMs (mostly Y181C and K103N/S) in this baseline study probably results from the widespread use of NNRTIs as part of standard first-line ART and PMTCT regimens. Most solitary NNRTI-associated DRMs cause a complete loss of activity of efavirenz and nevirapine,23 which may compromise the initial response to the standard first-line therapy.22 Moreover, NNRTI-associated DRMs have been shown to have only a modest impact on viral replicative fitness,23,24 allowing them to persist in the absence of the drug and to establish infection in a new host. Sex was not found to be associated with NNRTI resistance, suggesting either a limited effect on baseline resistance from PMTCT in females or secondary transmission of PMTCT-acquired NNRTI-resistant strains.
We found baseline NRTI-associated resistance to be limited. This contrasts not only with industrialized countries where transmission of NRTI-associated resistance is predominant3-6 but also with several recent African studies in subtype C-infected patients experiencing treatment failure which reported considerable rates of NRTI-associated DRMs, including TAMs and K65R.25-28 Several mechanisms could be expected to play a role. First, ARV rollout in sub-Saharan Africa is based on (highly active) ART, and widespread access has been established only recently, whereas in the industrialized world, ARV drugs have been widely used for many years including nonsuppressive mono and dual therapies with thymidine analogues in the past. Second, the reduced replicative fitness of variants harboring multiple TAMs and/or K65R might reduce transmission efficiency.29 Third, underestimation of NRTI-associated DRMs is possible due to reversion to wild type and/or outgrowth of minority wild-type species over time in absence of drug-selective pressure resulting in a reduction of the mutant strains with poor replicative fitness to minor variants below the limit of detection of genotypic analysis.30-33
Because of the infrequent use of PIs, we observed very few significant DRMs in protease. A few participants harbored clades D, G, CRF02_AG, and CRF09_cpx, which to our knowledge have not been described in Zambia before. Twenty specimens (3.5%) had unexpectedly low (<1000 copies/mL) HIV-1 RNA levels, which is most likely due to interindividual variations in viral replication rates and immunologic control34; other reasons could include partial viral suppression because of any undisclosed recent ARV exposure and varied assay performance between viral subtypes.35-37 As an additional observation, the presence of drug-resistant virus at baseline was found to be associated with a nonreduced serum hemoglobin level. A noncausal relation seems plausible, that is, advancing HIV infection leading to anemia as a result of bone marrow suppression,38 in parallel to the diminution of poorly replicating drug-resistant minor variants, as described above.30-33
The study has several limitations. First, it cannot be completely ruled out that reportedly ARV-naive participants had unknown previous exposure to therapy and/or prophylaxis. Second, the potential for selection bias exists, although the lack of heterogeneity in resistance pattern between the established private versus the 2 more recently introduced free ART programs argues against this. Third, data on route, country and duration of HIV-1 infection, and the source's ARV history, as possible associated factors of resistance, were not available.
In conclusion, this study on baseline HIV-1 resistance in routine ART programs in Lusaka, Zambia, adds important information regarding the predicted population-level response to standard first-line ART. Patients with previous ARV-experience are particularly at risk for a compromised initial response to standard NNRTI-based regimens. If baseline NNRTI resistance levels further increase, reassessment of first-line guidelines may be warranted to maintain individual and population-level benefits of ART. It is mandatory to monitor worldwide for the presence and spread of drug-resistant HIV-1.
The authors thank the study participants, the site staff, and the support staff at PharmAccess Foundation and Contract Laboratory Services. Special thanks to Ferdinand Wit (Academic Medical Center) and to David van de Vijver (Erasmus University Medical Center) for assistance in data analysis. The PharmAccess African Studies to Evaluate Resistance is part of the LAASER program (Linking African and Asian Societies for an Enhanced Response to HIV/AIDS), a partnership of Stichting Aids Fonds, The Foundation for AIDS Research (amfAR)-TREAT Asia, PharmAccess Foundation, and International Civil Society Support.
2. Gupta RK, Hill A, Sawyer AW, et al. Virological monitoring and resistance to first-line highly active antiretroviral therapy
in adults infected with HIV-1 treated under WHO guidelines: a systematic review and meta-analysis. Lancet Infect Dis
3. Grant RM, Hecht FM, Warmerdam M, et al. Time trends in primary HIV-1 drug resistance
among recently infected persons. JAMA
4. Little SJ, Holte S, Routy JP, et al. Antiretroviral-drug resistance among patients recently infected with HIV. N Engl J Med
5. Wensing AM, van dV, Angarano G, et al. Prevalence of drug-resistant HIV-1 variants in untreated individuals in Europe: implications for clinical management. J. Infect. Dis
6. Vercauteren J, Wensing AM, van de Vijver DA, et al. Transmission of drug-resistant HIV-1 is stabilizing in Europe. J Infect Dis
7. Abegaz WE, Grossman Z, Wolday D, et al. Threshold survey evaluating transmitted HIV drug resistance among public antenatal clinic clients in Addis Ababa, Ethiopia. Antivir Ther
. 2008;13(Suppl 2):89-94.
8. Kamoto K, Aberle-Grasse J. Surveillance
of transmitted HIV drug resistance with the World Health Organization threshold survey method in Lilongwe, Malawi. Antivir Ther
. 2008;13(Suppl 2):83-87.
9. Maphalala G, Okello V, Mndzebele S, et al. Surveillance
of transmitted HIV drug resistance in the Manzini-Mbabane corridor, Swaziland, in 2006. Antivir Ther
. 2008;13(Suppl 2):95-100.
10. Ndembi N, Lyagoba F, Nanteza B, et al. Transmitted antiretroviral drug resistance surveillance
among newly HIV type 1-diagnosed women attending an antenatal clinic in Entebbe, Uganda. AIDS Res Hum Retroviruses
11. Pillay V, Ledwaba J, Hunt G, et al. Antiretroviral drug resistance surveillance
among drug-naive HIV-1-infected individuals in Gauteng Province, South Africa in 2002 and 2004. Antivir Ther
. 2008;13(Suppl 2):101-107.
12. Somi GR, Kibuka T, Diallo K, et al. Surveillance
of transmitted HIV drug resistance among women attending antenatal clinics in Dar es Salaam, Tanzania. Antivir Ther
. 2008;13(Suppl 2):77-82.
13. Stringer JS, Zulu I, Levy J, et al. Rapid scale-up of antiretroviral therapy
at primary care sites in Zambia
: feasibility and early outcomes. JAMA
14. Ministry of Health Zambia
. Antiretroviral Therapy for Chronic HIV Infection in Adults and Adolescents. New ART Protocols
. Lusaka, Zambia
: Ministry of Health Zambia
15. Hamers RL, Kityo C, Siwale M, et al. Cohort profile: the PharmAccess African Studies to Evaluate Resistance Monitoring Study (PASER-M): HIV-1 Drug Resistance
in Sub-Saharan Africa. Presented at: XVIII International AIDS Conference; July 18-23, 2010; Vienna, Austria. Abstract no. THPE0236.
16. Wallis CL, Papathanasopoulos MA, Lakhi S, et al. Affordable in-house antiretroviral drug resistance assay with good performance in non-subtype B HIV-1. J Virol Methods
17. Bennett DE, Camacho RJ, Otelea D, et al. Drug resistance mutations for surveillance
of transmitted HIV-1 drug-resistance: 2009 update. PLoS One
18. de Oliveira T, Deforche K, Cassol S, et al. An automated genotyping system for analysis of HIV-1 and other microbial sequences. Bioinformatics
19. Handema R, Terunuma H, Kasolo F, et al. Prevalence of drug-resistance-associated mutations in antiretroviral drug-naive Zambians infected with subtype C HIV-1. AIDS Res Hum Retroviruses
20. Blower S, Bodine E, Kahn J, et al. The antiretroviral rollout and drug-resistant HIV in Africa: insights from empirical data and theoretical models. AIDS
21. Apisarnthanarak A, Jirayasethpong T, Sa-nguansilp C, et al. Antiretroviral drug resistance among antiretroviral-naive persons with recent HIV infection in Thailand. HIV Med
22. Barth RE, Wensing AM, Tempelman HA, et al. Rapid accumulation of nonnucleoside reverse transcriptase inhibitor-associated resistance: evidence of transmitted resistance in rural South Africa. AIDS
23. Richman DD, Havlir D, Corbeil J, et al. Nevirapine resistance mutations of human immunodeficiency virus type 1 selected during therapy. J. Virol
24. Nijhuis M, Deeks S, Boucher C. Implications of antiretroviral resistance on viral fitness. Curr. Opin Infect Dis
25. Hosseinipour MC, van Oosterhout JJ, Weigel R, et al. The public health approach to identify antiretroviral therapy
failure: high-level nucleoside reverse transcriptase inhibitor resistance among Malawians failing first-line antiretroviral therapy
26. Marconi VC, Sunpath H, Lu Z, et al. Prevalence of HIV-1 drug resistance
after failure of a first highly active antiretroviral therapy
regimen in KwaZulu Natal, South Africa. Clin Infect Dis
27. Orrell C, Walensky RP, Losina E, et al. HIV type-1 clade C resistance genotypes in treatment-naive patients and after first virological failure in a large community antiretroviral therapy
programme. Antivir Ther
28. Wallis C, Sanne I, Venter F, et al.Varied patterns of HIV-1 drug resistance
on failing first-line antiretroviral therapy
in South Africa. J Acquir Immune Defic Syndr
29. Yerly S, Jost S, Telenti A, et al. Infrequent transmission of HIV-1 drug-resistant variants. Antivir Ther
30. Weber J, Chakraborty B, Weberova J, et al. Diminished replicative fitness of primary human immunodeficiency virus type 1 isolates harboring the K65R mutation. J Clin Microbiol
31. White KL, Margot NA, Wrin T, et al. Molecular mechanisms of resistance to human immunodeficiency virus type 1 with reverse transcriptase mutations K65R and K65R+M184V and their effects on enzyme function and viral replication capacity. Antimicrob Agents Chemother
32. Johnson JA, Li JF, Wei X, et al. Minority HIV-1 drug resistance
mutations are present in antiretroviral treatment-naive populations and associate with reduced treatment efficacy. PLoS Med
33. Schuurman R, Brambilla D, de Groot T, et al. Underestimation of HIV type 1 drug resistance mutations: results from the ENVA-2 genotyping proficiency program. AIDS Res Hum Retroviruses
34. Novitsky V, Woldegabriel E, Kebaabetswe L, et al. Viral load and CD4+
T-cell dynamics in primary HIV-1 subtype C infection. J Acquir Immune Defic Syndr
35. Gottesman BS, Grosman Z, Lorber M, et al. Measurement of HIV RNA in patients infected by subtype C by assays optimized for subtype B results in an underestimation of the viral load. J Med Virol
36. Holguin A, Lopez M, Molinero M, et al. Performance of three commercial viral load assays, Versant human immunodeficiency virus type 1 (HIV-1) RNA bDNA v3.0, Cobas AmpliPrep/Cobas TaqMan HIV-1, and NucliSens HIV-1 EasyQ v1.2, testing HIV-1 non-B subtypes and recombinant variants. J Clin Microbiol
37. Rouet F, Foulongne V, Viljoen J, et al. Comparison of the Generic HIV Viral Load assay with the Amplicor HIV-1 monitor v1.5 and Nuclisens HIV-1 EasyQ v1.2 techniques for plasma HIV-1 RNA quantitation of non-B subtypes: the Kesho Bora preparatory study. J Virol Methods
38. Claster S. Biology of anemia, differential diagnosis, and treatment options in human immunodeficiency virus infection. J Infect Dis
. 2002;185(Suppl 2):S105-S109.