Abstract
Background
We assessed the prevalence of acquired HIV drug resistance (HIVDR) and associated factors among patients receiving first-line antiretroviral therapy (ART) in Rwanda.
Methods
This cross-sectional study included 702 patients receiving first-line ART for at least 6 months with last viral load (VL) results ≥1000 copies/mL. Blood plasma samples were subjected to VL testing; specimens with unsuppressed VL were genotyped to identify HIVDR-associated mutations. Data were analysed using STATA/SE.
Results
Median time on ART was 86.4 months (interquartile range [IQR], 44.8–130.2 months), and median CD4 count at ART initiation was 311 cells/mm3 (IQR, 197–484 cells/mm3). Of 414 (68.2%) samples with unsuppressed VL, 378 (88.3%) were genotyped. HIVDR included 347 (90.4%) non-nucleoside reverse transcriptase inhibitor- (NNRTI), 291 (75.5%) nucleoside reverse transcriptase inhibitor- (NRTI) and 13 (3.5%) protease inhibitor (PI) resistance-associated mutations. The most common HIVDR mutations were K65R (22.7%), M184V (15.4%) and D67N (9.8%) for NRTIs and K103N (34.4%) and Y181C/I/V/YC (7%) for NNRTIs. Independent predictors of acquired HIVDR included current ART regimen of zidovudine + lamivudine + nevirapine (adjusted odds ratio [aOR], 3.333 [95% confidence interval (CI): 1.022–10.870]; p = 0.046) for NRTI resistance and current ART regimen of tenofovir + emtricitabine + nevirapine (aOR, 0.148 [95% CI: 0.028–0.779]; p = 0.025), zidovudine + lamivudine + efavirenz (aOR, 0.105 [95% CI: 0.016–0.693]; p = 0.020) and zidovudine + lamivudine + nevirapine (aOR, 0.259 [95% CI: 0.084–0.793]; p = 0.019) for NNRTI resistance. History of ever switching ART regimen was associated with NRTI resistance (aOR, 2.53 [95% CI: 1.198–5.356]; p = 0.016) and NNRTI resistance (aOR, 3.23 [95% CI: 1.435–7.278], p = 0.005).
Conclusion
The prevalence of acquired HIV drug resistance (HIVDR) was high among patient failing to re-suppress VL and was associated with current ART regimen and ever switching ART regimen. The findings of this study support the current WHO guidelines recommending that patients on an NNRTI-based regimen should be switched based on a single viral load test and suggests that national HIV VL monitoring of patients receiving ART has prevented long-term treatment failure that would result in the accumulation of TAMs and potential loss of efficacy of all NRTI used in second-line ART as the backbone in combination with either dolutegravir or boosted PIs.
Introduction
HIV prevention and monitoring HIV drug resistance (HIVDR) maximize the long-term effectiveness of first-line antiretroviral therapy (ART), optimize patient outcomes, minimize transmission of drug-resistant HIV and ensure sustainability of ART programs.1,2 However, resource-limited settings often have limited ART regimens and limited access to routine viral load (VL) testing to monitor treatment outcomes. 1 Even in settings with optimal ART program management, HIVDR is expected to emerge in populations receiving ART, which could affect the response to second-line ART and contribute to transmission of drug-resistant HIV. 3
In July 2016, the Rwanda ART program initiated the Treat All policy, in which all people living with HIV (PLHIV) initiate ART within the same week (preferably on the same day) as diagnosis and undergo VL testing at 6 months following ART initiation. Subsequently, patients receiving ART undergo VL testing every 12 months. Patients with potential virologic failure receive 3 months of enhanced adherence counselling and close monitoring to address adherence-related factors and undergo VL testing after potential adherence issues are addressed. 5
The standard first-line ART regimens for adults in Rwanda were tenofovir/lamivudine in combination with one non-nucleoside reverse transcriptase inhibitor (NNRTI; efavirenz as the first option or nevirapine) or abacavir/lamivudine plus one NNRTI (efavirenz as the first option or nevirapine). In July 2018, dolutegravir was introduced as first-line ART for ART-naive adults and adolescents who weigh >35 kg. The standard second-line ART regimen is zidovudine, lamivudine and atazanavir/ritonavir or lopinavir/ritonavir according to national ART guidelines adopted in 2016 or a dolutegravir-based regimen after HIV genotyping. 5
At the end of 2018 in Rwanda, of the >200,000 PLHIV, 83.1% were receiving ART in 539 health facilities. 4 The extent of ART coverage in Rwanda increases the likelihood of HIVDR among PLHIV receiving first-line regimens. Previous studies in Rwanda have characterized transmitted drug resistance 6 and acquired drug resistance to first-line 7 and second-line regimens. 8 However, previous studies on HIVDR or ART failure in Rwanda are outdated, had small sample sizes or included only patients with a short treatment duration7,9,10; no nationally representative estimate is available for HIVDR in patients for whom first-line ART has failed. Although the World Health Organization (WHO) recommends that HIV treatment scale-up be accompanied by a robust assessment of drug resistance emergence and transmission, in Rwanda, these factors are monitored only when second-line ART has failed and patients need to switch to a third-line ART regimen. The most recent study on acquired HIVDR for patients receiving first-line regimens in Rwanda was conducted in 2011. We conducted a survey to estimate national prevalence of HIVDR and VL re-suppression among adults living with HIV who had received first-line ART for at least 6 months and who had failed to re-suppress their VL (≥1000 copies/mL).
Methods
Study setting, design and sampling
We conducted a cross-sectional national study and collected blood specimens and minimal information, including sociodemographic and clinical characteristics of patients. Study design followed a two-stage sampling survey design without neither clustering nor stratification (referred as the standard complex survey design). 11 In the first stage, 66 health facilities were selected from the total national sites dispensing ART in the country using probability proportional to health facility size. In the second stage, eligible patients were selected randomly based on health facility’s allocated minimum sample size. The sample size was calculated using formula assumptions and conditions adapted from the WHO generic protocol, which was adjusted to the specific context of Rwanda 11 : The WHO assumed prevalence of viral suppression of 91% which was replaced with a viral re-suppression of 70%,12–14. drug resistance prevalence of 50%, 15 assumed percentage of patients on first-line therapy at time of survey of 95.8%, assumed polymerase chain reaction amplification rate of 80%, desired precision of ±5% reduced to ±3.75, design effect for sampling of 1.5 and response rate percent of 75%13,14 calculated using the following formula.
Assumed chain reaction amplification rate × (1-% of LTFU) × (1-% of deaths) × (1-% of net transfers) using lost to follow-up (LTFU) per year: 2.56%, death per year: 0.91%, transfer in per year; 4.7%, net transfers per year: 1.1%, which is equal to (% of transfer outs − % of transfer ins).
The survey was conducted in 67 health facilities, selected from more than 539 health facilities providing ART in Rwanda. Based on these sample size estimation assumptions and an anticipated NNRTI drug resistance prevalence of 50%, a sample of 750 individuals will provide an effective sample size of 203 individuals across 67 health facilities with a precision that corresponds with an alpha of ±7% and a 95% confidence interval (CI) level. An effective sample size of 377 individuals used in this study provided an estimate with a precision that corresponds with an alpha of ±5%.
The survey included all PLHIV aged ≥18 years who had received HIV first-line ART for at least 6 months, had VL ≥1000 copies/mL from the most recent test at the time of the survey inception and whose latest VL test results were obtained not more than 3 months before the date of sample collection. The study did not include individuals aged <18 years, patients reinitiating ART because of treatment failure, patients receiving second-line or third-line regimens and patients who did not provide written informed consent for participation in the survey.
Data collection
The list of patients receiving first-line ART who met the eligibility criteria in each health facility was obtained from the ART patient register. All eligible patients were contacted and invited to participate in the study during an appointment at the health facilities. We used a structured questionnaire programmed before in Personal Digital Assistant devices to collect patients’ demographic and clinical information from patient charts/registers. We followed WHO recommendations to collect blood samples from study participants for HIVDR testing. 16
Patient flow
Eligible patients with latest VL results ≥1000 copies/mL (recorded in patient file) were listed and contacted for blood draw at the health facilities, and samples were stored at the national reference laboratory. Collected samples were subjected to a second VL test, after which all samples with VL≥1000 copies/mL were genotyped to identify potential HIVDR mutations (Figure 1). Cascade of patient enrolment based on eligibility criteria, blood samples with unsuppressed VL (VL ≥1000 copies/mL) for genotyping and number of patients with drug-resistant mutation types (based on weighted %). * The presented % are weighted prevalence estimates.
HIVDR genotyping
Whole-blood samples were processed at the Rwanda National Reference Laboratory for VL testing and HIVDR genotyping. VL testing was conducted using the Roche COBAS AmpliPrep/COBAS TaqMan HIV-1 for plasma samples. Specimens with unsuppressed VL were subjected to Thermo Fisher Scientific HIVDR genotyping test to detect ART-resistant mutations via sequencing reverse transcriptase and protease genes. Additionally, 65.6% (248) of the HIVDR sequences’ raw data and Fasta files were shared with the WHO-designated Specialized Drug Resistance Laboratory (SDRL) at the U.S. Centers for Disease Control and Prevention (CDC), Atlanta, Georgia, USA for quality control and to ensure correct interpretation of individual patients’ HIVDR-genotyping results used for this study. The analysis of the HIVDR genotyping results and interpretation for the drug resistance profiles was done using Stanford University HIV Drug Resistance Database. Based on the Stanford HIVdb algorithm the following scoring criteria were used: There are 5 drug resistance levels: 1 indicates susceptible, 2 indicates potential low-level resistance, 3 indicates low-level resistance, 4 indicates intermediate resistance, 5 indicates high-level resistance. The scores are the sum of each mutation penalty score for a drug. Scores less than 10 indicate susceptible; scores between 10 and 14 indicate potential low-level resistance; scores between 15 and 29 indicate low-level resistance; scores between 30 and 59 indicate intermediate resistance. Scores of 60 or greater indicate high-level resistance.
Data analysis
Univariate descriptive statistics, specifically frequency counts and corresponding percentages, were used to describe sample characteristics, social demographics and clinical characteristics, acquired HIVDR testing cascade and various types of mutations present in the sample. Additionally, bivariate analysis with chi-squared tests was used to identify factors associated with failure to re-suppress VL and each acquired drug-resistant mutation class (nucleoside reverse transcriptase inhibitors [NRTIs], NNRTIs and protease inhibitors [PIs]), with 0.05 threshold of alpha level associated with 95% CI level (p-value). Multivariate logistic regression models were fit for all significant predictors, with p-value <0.05 considered statistically significant for each independent outcome (failure to re-suppress and each acquired drug-resistant mutation class). Stepwise backward logistic regression modelling was used to produce adjusted odds ratios and corresponding 95% CI. Per WHO guidelines, 11 the analyses accounted for the standard complex survey design. The standard complex survey design was employed because neither clustering nor stratification was done during sampling. 11 STATA/SE software (version 16, Stata Corp LLC) was used for data analysis.
Ethical consideration
This survey was reviewed and approved by the Rwanda National Ethics Committee and was reviewed in accordance with the US Centers for Disease Control and Prevention (CDC) human research protection procedures and was determined to be research, but CDC investigators did not interact with human subjects or have access to identifiable data or specimens for research purposes. All participants signed the informed consent form before enrolling in the study.
Results
Of the 702 patients included in this survey, 469 (66.8%) were women and most (334 [47.6%]) were aged 35–49 years. Median CD4 count at ART initiation was 311 cells/mm3 (interquartile range [IQR], 197–484 cells/mm3), and median time on ART was 86.4 months (IQR, 44.8–130.2 months). Of the participants, 68.3% had received treatment for ≥5 years and 37.3% of those receiving ART had a history of switching ART regimens at least once. The most common regimens at the time of specimen collection were tenofovir + lamivudine + efavirenz (39.3%), tenofovir + lamivudine + nevirapine (17.2%) and zidovudine + lamivudine + nevirapine (16.8%).
HIVDR-testing cascade and geographical distribution of HIVDR mutations
Nucleotide sequencing in the reverse transcriptase gene identified NNRTIs and NRTIs. The number and weighted prevalence of any drug-resistant mutations was 355 [92.7%, 95% CI: 86.2–96.3], of NRTIs was 291 [75.5%, 95% CI: 66.7–82.5], of NNRTIs was 347 [90.4%, 95% CI: 81.5–95.2), of PIs was 13 [3.5%, 95% CI: 1.4–8.7] (Figure 1). Figure 2 presents the geographical distribution by proportion of NRTI-resistant and NNRTI-resistant mutations with corresponding patient number (N) enrolled in study health facilities. Both NRTI-resistant and NNRTI-resistant mutations were concentrated in the City of Kigali and in the Western province. Maps of Rwanda showing the geographic distribution and prevalence of NNRTI-resistant (panel A) and NRTI-resistant (panel B) mutations in samples tested per site with VL >1000 copies/mL.
Quality Control of HIVDR genotyping
The quality control (QC) and quality assurance process used in this study examined the sequence data at the batch and survey level according to recommended WHO HIVResNet HIVDR laboratory operational framework. 17 At the batch level, sequences were analysed by the checking the raw chromatogram quality, sequence length, single-stranded coverage, presence of stop codons, out-of-frame insertions or deletions, missing or gaps in the sequences, highly unusual mutations, excessive ambiguity, excessive number of mixed bases and APOBEC mutations. However, at survey level, all sequences were examined to identify expected and unexpected sequence similarity between pairs of sequences that were analysed in different batches to generate a sequence identity matrix with confirmation done by using the online drug resistance quality control tool developed by WHO and British Columbia Centre for Excellence in HIV/AIDS. The sequences that passed QC were submitted to Stanford's HIV drug resistance database (HIVdb) where the same parameters were checked again in addition to obtaining the drug resistance profile for this study. A total of 1572 raw data (ab1 files) and 262 assembled sequences (Fasta files) were submitted from NRL and examined by the CDC International Laboratory Branch through workflow process described above for the final drug resistance profiles, of which 158 sequences (60%) were 100% identical, 80 sequences (31%) with 1–3 nucleotides and 24 sequences (9%) with up to 15 nucleotides differed mostly because of nucleotide mixture calling. The overall greater proportion of identical sequences and some with a few differing nucleotides (>90%) indicates that NRL staff are well trained and competent in the genotyping process to generate highly reliable HIVDR results reported in this study.
Drug-resistant mutations by drug class and current ART regimen
NRTI-resistant mutations associated with first-line ART regimens among adults living with HIV (n = 286) in Rwanda.
Abbreviations: NNRTI, non-nucleoside reverse transcriptase inhibitors; ART, antiretroviral therapy; AZT, zidovudine; 3TC, lamivudine; EFV, efavirenz; NVP, nevirapine; TDF, tenofovir disoproxil fumarate; FTC, emtricitabine; LPV/r, lopinavir/ritonavir; ABC, abacavir; DTG, dolutegravir.
NNRTI-resistant mutations associated with first-line ART regimens among adults living with HIV (n = 342) in Rwanda.
Abbreviations: NNRTI, non-nucleoside reverse transcriptase inhibitors; ART, antiretroviral therapy; AZT, zidovudine; 3TC, lamivudine; EFV, efavirenz; NVP, nevirapine; TDF, tenofovir disoproxil fumarate; FTC, emtricitabine; LPV/r, lopinavir/ritonavir; ABC, abacavir; DTG, dolutegravir.
Multivariable logistic regression to determine factors associated with failure to re-suppress the VL among adults living with HIV in Rwanda.
Predictors of failure to re-suppress VL
In the univariate analysis, latest CD4 categories and current ART regimen were significantly associated with failure to re-suppress (VL ≥1000 copies/mL). Those with latest CD4 ranging between 200 and 349 cells/
Predictors of HIVDR mutations
Univariate and multivariable analyses were conducted to evaluate variables associated with HIVDR mutations.
Multivariable logistic regression to determine factors associated with any NRTI-resistant mutations among adults living with HIV (n = 378) in Rwanda.
Abbreviations: NRTI, nucleoside reverse transcriptase inhibitors; aOR, adjusted odds ratio; CI, confidence interval; ART, antiretroviral therapy, AZT, zidovudine; 3TC, lamivudine; EFV, efavirenz; NVP, nevirapine; TDF, tenofovir disoproxil fumarate; FTC, emtricitabine; LPV/r, lopinavir/ritonavir; ABC, abacavir.
Multivariable logistic regression to determine factors associated with any NNRTI-resistance mutations among adults living with HIV (n = 378) in Rwanda.
Abbreviations: NNRTI, non-nucleoside reverse transcriptase inhibitors; aOR, adjusted odds ratio; CI, confidence interval; ART, antiretroviral therapy; AZT, zidovudine; 3TC, lamivudine; EFV, efavirenz; NVP, nevirapine; TDF, tenofovir disoproxil fumarate; FTC, emtricitabine; LPV/r, lopinavir/ritonavir; ABC, abacavir.
Discussion
This is the first national cross-sectional study examining the prevalence and distribution of HIVDR mutations associated with first-line regimen failure in the Rwanda national HIV prevention, care and treatment program. Of all patients with unsuppressed VL at study inception using VL collected 3 months prior to the study, 31.8% had suppressed VL at the second VL measurement (VL test done using samples collected in the study), indicating that most patients enrolled had failed to re-suppress their VL. These patients once on an NNRTI-based regimen should be switched to a different ART regimen based on a single viral load test and switched based on a second VL test after enhanced adherence counselling once on a non-NNRTI based regimen as recommended by the 2021 WHO guideline. 18
In this study, 92.7% of patients had at least one drug-resistant mutation against NRTIs, NNRTIs or PIs used in the national program. The resistant mutations identified against different drug classes were 90.4% for NNRTIs, 75.5% for NRTIs and 3.5% for PIs. Other studies conducted elsewhere in adults aged ≥18 years have found varying results for prevalence of NRTI-resistant and NNRTI-resistant mutations. Studies have found 93.11% of patients with NNRTI-resistant mutations, 74.40% of patients with NRTI-resistant mutations in Hunan South China and 93% and 96% of patients with NRTI-resistant and NNRTI-resistant mutations, respectively, in resource-limited settings (i.e. Thailand, South Africa, India, Malawi and Tanzania).19,20 Another study found substantially lower drug-resistant mutation rates: 11% of patients had NRTI-associated drug-resistant mutations and 17% of patients had no known mutations, depending on the VL monitoring practice in place. 21
The most common HIVDR mutations identified in our study were K65R, M184V and D67N for NRTIs and K103N, A98G and Y181C/I/V for NNRTIs. It is likely that some of the patients on TDF-based first-line therapy previously acquired the D67N mutation while receiving zidovudine. Overall, the NRTI-resistant mutations observed in our study allow effective NRTIs to remain as backbone options for second-line ART. Tenofovir disoproxil fumarate genotypic resistance, defined by K65R/N and/or K70E/Q/G, occurs in 20%–60% of individuals with unsuppressed VL receiving WHO-recommended tenofovir-containing first-line regimens. 22
In our study, 22% of patients had the K65R/N mutation and 8% had K70E/G/Q. This suggests that one of three NRTI-resistant mutations developed against tenofovir. This high prevalence could be influenced by the fact that most patients in our study were receiving tenofovir-based ART. K65R/N and K70E/G/Q are rarely combined because of the reduced virus replication capacity when the two mutations occur together. The second most prevalent NRTI-resistant mutation was M184V/I (15%), which is selected by lamivudine or emtricitabine and is reported as the most common NRTI-resistant mutation after first-line ART failure.23,24 Although M184V causes high-level resistance to lamivudine and emtricitabine, the M184V/I mutation also is associated with reduced viral fitness and increased virus susceptibility to zidovudine and tenofovir. For this reason, lamivudine/emtricitabine is usually maintained in the backbone of second-line combination ART. 25 Our finding that K65R was more common than M184V is unique compared to other studies which were conducted in prior years and merits further exploration and validation as part of drug resistance surveillance in other countries. 26
In addition to common NRTI-resistant mutations, our study identified A62AV mutations (2.1%), which is also a tenofovir accessory resistant mutation that often occurs in combination with the multi-NRTI-resistant mutations K65R or Q151M. A62V is widespread in subtype A. However, we also found D67D/N and K219K/Q/R (0.3%), which are thymidine analogues mutations (TAMs) and not tenofovir induced. M41l, E44D, K70R, M184V, K219Q (0.3%), M41Ml (0.3.1%) and V75M (0.7%) are TAMs and are not selected by tenofovir. These findings suggest that first-line ART failed or that there might have been some patients transferred to different health facilities and initiated tenofovir-based ART after a zidovudine/stavudine-based regimen failed, but the patients did not disclose their ART history.
Furthermore, TAMS/multiple TAMS in the study were rare in our study, though some of the patients had received zidovudine for a long period. This suggests that national HIV VL monitoring of patients receiving ART has prevented long-term treatment failure that would result in the accumulation of TAMs and potential loss of efficacy of all NRTI. In contrast to our NRTI-resistant mutations findings, our NNRTIs results would exclude any further use of first-generation NNRTIs (efavirenz and nevirapine) because of the cross-resistance of NNRTI mutations. The emergence of Y181C/I/V, mostly selected by nevirapine, would also limit the use of etravirine, which is reserved as a rescue drug for patients for whom second-line regimens have failed. K101E decreases HIV susceptibility to all NNRTIs (efavirenz, nevirapine and etravirine). This suggests that national HIV VL monitoring of patients receiving ART has prevented long-term treatment failure that would result in the accumulation of TAMs and potential loss of efficacy of all NRTI used in second-line ART as the backbone in combination with either dolutegravir or boosted PIs.
We found that switching the initial ART regimen was significantly associated with developing NRTI-resistant or NNRTI-resistant mutations. These results could be because of the fact that patients switched their regimens either because of side effects or ART toxicity. The low barrier to resistance of NNRTIs used in first-line ART might lead to resistant mutations from the time of clinical side effects/toxicity to the time of switching the initial ART regimen. In addition, we found that patients receiving tenofovir, lamivudine and efavirenz were more likely to develop NNRTI-resistant mutations than patients on other regimens. In 2009, the tenofovir-based regimen was introduced as the preferred first-line ART in Rwanda. During our study, patients received zidovudine-based ART if they had contraindications to use tenofovir or nevirapine. It is possible that patients stopped their initial ART before the treatment was changed and that NNRTI-resistant mutations developed during that period.
Our findings are subject to several limitations. Some of the assumptions used for sample size determination were posteriorly determined, and therefore, conservative assumption was used to maximize sample size. The sample size estimation was therefore based on more stringent assumptions derived from credible and relevant references in order not to under-estimate the required sample size. We included only patients with unsuppressed VL, for whom first-line ART had failed. This VL threshold cannot capture drug-resistant mutations that occur at lower load viraemia and that are associated with subsequent virologic failure. Because our study was conducted at the time dolutegravir was first introduced in Rwanda, it was not possible to characterize potential dolutegravir-resistant mutations. Our logistic regression analysis only included variables that were available. Because the data were collected from patient charts (files), which are associated with missing variable bias, we might have missed or not captured variables that might have changed our results.
NRTI-resistant mutations preserve effective second-line ART with selected NRTIs as the backbone in combination with either dolutegravir or boosted PI. 22 The prevalence of acquired HIV drug resistance (HIVDR) was high among patient failing to re-suppress VL and was associated with current ART regimen and ever switching ART regimen. The findings of this study support the current WHO guidelines recommending that patients on an NNRTI-based regimen should be switched based on a single viral load test and suggests that national HIV VL monitoring of patients receiving ART has prevented long-term treatment failure that would result in the accumulation of TAMs and potential loss of efficacy of all NRTI used in second-line ART as the backbone in combination with either dolutegravir or boosted PIs. This study further supports the bi-annual revisions and updates to the national HIV treatment guidelines.
Footnotes
Acknowledgements
We would like to thank the Government of Rwanda, particularly the Rwanda Biomedical Center (RBC) survey team, for their dedication to this study and the participating health facilities and patients for their collaboration. In addition, we would also like to acknowledge technical support received from our partners at CDC and WHO.
Authors’ contributions
GM contributed to study conception, study design, data collection and analysis and led manuscript writing. MR, EK, SSM, RCNM, JMS, ABS and ER contributed to study conception and study design, drafted and edited the manuscript and interpreted the results. ET, SSM, RCNM, AK, GM, GNR and ER were responsible for conducting statistical analysis, interpreting results and drafting and editing the manuscript. GM, MRH, JD, CB, MR, SDH and ER contributed to collecting data, managing field activities, manuscript editing and interpreting results. MR, GNR and ER supervised the study and protocol implementation and contributed to manuscript development and interpretation of results. All authors read and approved the final submitted version of the manuscript.
Declaration of conflicting interests
The author(s) declared the following potential conflicts of interest with respect to the research, authorship, and/or publication of this article: We declare a conflict of interest. Josh DeVos is a co-inventor in U.S. patent US10053741B2.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study is supported by the President’s Emergency Plan for AIDS Relief (PEPFAR) through the Centers for Disease Control and Prevention (CDC) under the terms of cooperative agreement #GH001612-04. The findings and conclusions in this manuscript are those of the authors and do not necessarily represent the official position of the funding agencies.
