Back to Journals » Infection and Drug Resistance » Volume 15

HIV-1C in-House RNA-Based Genotyping Assay for Detection of Drug Resistance Mutations in Samples with Low-Level Viral Loads

Authors Bareng OT, Choga WT, Maphorisa ST, Seselamarumo S, Seatla KK , Mokgethi PT, Maruapula D , Mogwele ML, Ditshwanelo D, Moraka NO, Gobe I, Motswaledi MS, Makhema JM, Musonda R, Shapiro R, Essex M, Novitsky V, Moyo S , Gaseitsiwe S

Received 19 September 2022

Accepted for publication 8 November 2022

Published 22 December 2022 Volume 2022:15 Pages 7565—7576

DOI https://doi.org/10.2147/IDR.S388816

Checked for plagiarism Yes

Review by Single anonymous peer review

Peer reviewer comments 2

Editor who approved publication: Prof. Dr. Héctor Mora-Montes



Ontlametse T Bareng,1,2 Wonderful T Choga,1,2 Segomotso T Maphorisa,3 Sekgabo Seselamarumo,1 Kaelo K Seatla,1 Patrick T Mokgethi,1,4 Dorcas Maruapula,1,4 Mompati L Mogwele,1 Doreen Ditshwanelo,1,5 Natasha O Moraka,1 Irene Gobe,2 Modisa S Motswaledi,2 Joseph M Makhema,1,6 Rosemary Musonda,1 Roger Shapiro,1,6 Max Essex,1,6 Vlad Novitsky,1 Sikhulile Moyo,1,2,6 Simani Gaseitsiwe1,6

1Botswana Harvard AIDS Institute Partnership, Gaborone, Botswana; 2School of Allied Health Professions, Faculty of Health Sciences, University of Botswana, Gaborone, Botswana; 3Ministry of Health and Wellness, Republic of Botswana, Gaborone, Botswana; 4Department of Biological Sciences, Faculty of Science, University of Botswana, Gaborone, Botswana; 5Department of Biological Science and Biotechnology, Botswana International University of Science and Technology, Palapye, Botswana; 6Department of Immunology and Infectious Diseases, Harvard T.H. Chan School of Public Health, Boston, MA, USA

Correspondence: Simani Gaseitsiwe, Botswana Harvard AIDS Institute Partnership, Private Bag BO320, Bontleng, Gaborone, Botswana, Tel +267 390 2671, Fax +267 390 1284, Email [email protected]

Purpose: Monitoring HIV-1 drug resistance mutations (DRM) in treated patients on combination antiretroviral therapy (cART) with a detectable HIV-1 viral load (VL) is important for the selection of appropriate cART. Currently, there is limited data on HIV DRM at low-level viremia (LLV) (VL 401– 999 copies/mL) due to the use of a threshold of VL ≥ 1000 copies/mL for HIV DRM testing. We here assess the performance of an in-house HIV drug resistance genotyping assay using plasma for the detection of DRM at LLV.
Methods: We used a total of 96 HIV plasma samples from the population-based Botswana Combination Prevention Project (BCPP). The samples were stratified by VL groups: 50 samples had LLV, defined as 401– 999 copies/mL, and 46 had ≥ 1000 copies/mL. HIV pol (PR and RT) region was amplified and sequenced using an in-house genotyping assay with BigDye sequencing chemistry. Known HIV DRMs were identified using the Stanford HIV Drug Resistance Database. Genotyping success rate between the two groups was estimated and compared using the comparison of proportions test.
Results: The overall genotyping success rate was 79% (76/96). For VL groups, the genotyping success was 72% (36/50) at LLV and 87% (40/46) at VL ≥ 1000 copies/mL. Among generated sequences, the overall prevalence of individuals with at least 1 major or intermediate-associated DRM was 24% (18/76). The proportions of NNRTI-, NRTI- and PI-associated resistance mutations were 28%, 24%, and 0%, respectively. The most predominant mutations detected were K103N (18%) and M184V (12%) in NNRTI- and NRTI-associated mutations, respectively. The prevalence of DRM was 17% (6/36) at LLV and 30% (12/40) at VL ≥ 1000 copies/mL.
Conclusion: The in-house HIV genotyping assay successfully genotyped 72% of LLV samples and was able to detect 17% of DRM amongst them. Our results highlight the possibility and clinical significance of genotyping HIV among individuals with LLV.

Keywords: in-house genotyping, low-level viremia, samples, HIV-1C drug resistance testing

Introduction

A subset of people with human immunodeficiency virus (PWH) on potent antiretroviral therapy (ART) have low-level viremia (LLV), a detectable HIV viral load (VL) below 1000 copies/mL.1,2 Lack of HIV suppression may be attributed to many factors such as the presence of drug resistance mutations (DRM), metabolic complications affecting the pharmacokinetics of the drugs, and lack of medication adherence.2–4 Recent studies have found the presence of DRMs in LLV to be a strong predictor of subsequent virologic failure (VF).5–8 Despite these recent findings, no clear conclusion can be drawn on the impact of DRM on individuals with LLV and its long-term clinical outcomes because of the limited availability of HIV drug resistance genotyping assays validated for samples with LLV.

Most commercial RNA-based genotyping assays require HIV plasma viral loads > 1000 copies/mL for increased amplification success rates and accurate results.9,10 The other reason for limited data on DRM at LLV is linked to the threshold of 1000 copies/mL to define virological failure and threshold for genotyping testing, especially in low and middle-income settings. HIV proviral DNA is very useful in LLV or when plasma sequencing is not successful.11–13 However, proviral DNA-based assays are not routinely utilized for drug resistance testing for clinical monitoring purposes, which highlights the need for optimization of RNA-based genotyping assays for LLV. There is the possibility that genotyping samples with low VL using RNA-based approaches are limited by small copy numbers, which may render inaccurate genotyping results. However, there are optimized HIV genotyping assays for LLV with high success rates which utilize RNA-based approaches.10,14–17 The ability to genotype LLV could contribute to improved management of PWH experiencing low-level viremia, especially for the selection of appropriate ART, preserving future ART options. Therefore, our study aimed to optimize and assess the performance of an in-house RNA-based HIV genotyping assay for LLV in HIV-1 subtype C which predominates in Botswana and is the most prevalent HIV-1 subtype globally. Botswana is one of the countries with a high HIV prevalence of 18.6% among 15–49-year-olds with 7200 individuals newly infected in 2021.18 However, Botswana has had major successes in the country-wide implementation of antiretroviral therapy (ART) and high viral load suppression rates.19 The country has reached UNAIDS 95–95-95 UNAIDS targets where 95.1% of the population is aware of their HIV status, 98.0% are on ART and 97.9% are virally suppressed.20 Despite major advances in the development of antiretroviral (ARV) drugs and ARV therapy (ART) treatment guidelines,21–23 Botswana like other middle-income countries continues to face challenges including development, transmission and spread of HIV drug resistance mutations.24 This highlights the need for HIV drug resistance monitoring at any detectable VL with LLV included.

Methods

Study Population

Stored plasma samples from whole blood collected from PWH aged 16–64 years who were enrolled in the Botswana Combination Prevention Project (BCPP) from 2013 to 2018 and residing in 30 communities across central, northern, and southern parts of Botswana were utilized. BCPP was a community-randomized trial evaluating whether a package of standard HIV prevention trials would lessen the frequency of HIV cases over time among ART-experienced and -naïve individuals in Botswana, as previously described.19,25

Selection of Study Samples

Abbott m2000sp/rt (Abbott Laboratories, Abbott Park, IL, USA) was used for the quantification of HIV-1 VL with the lowest limit of detection at 40 copies/mL and the highest limit at 10,000,000 copies/mL. A total of 96 samples with VL > 400 copies/mL and enough volume ≥ 200μL were conveniently sampled from remnants of the BCPP cohort for the assay optimization, with 52% at LLV and 48% at VL ≥ 1000 copies/mL. Samples with VL 401–999 copies/mL were categorized as LLV, while those with VL ≥ 1000 copies/mL were considered to be in virologic failure (VF),26 per WHO guidelines.

Nucleic Acid Purification

HIV RNA was manually extracted from 200μL plasma samples stored at −80°C using the Zymo Research Quick-RNA Viral kit (Zymo Research, Pretoria, South Africa) as per the manufacturer’s instructions. The isolated RNA was used as a first-round PCR template and stored at −80°C for future use.

Amplification of HIV Pol

The HIV Pol (PR and RT) region of 1.1 kb, spanning the entire protease and the first 800 bps of RT, was amplified with a transcriptor one-step RT-PCR kit (Roche Applied Science, Penzberg, Germany). A one-step RT-PCR reaction mixture was prepared to consist of 0.5µL Transcriptase Enzyme Roche One Step, 7µL of RNase-free water, 5µL of transcriptor one-step RT-PCR 5X Buffer, 2.5µL mixture of primers CWF1-LNA2(5’-GAA GGA CCA AAT GAA AGA YTG-3’ (2 μM) and CWR1-LNA3 (5’-GCA TAC TTY CCG TTT TCA G-3’ (2 μM), and 10µL of HIV RNA, making a total of 25µL. The PCR parameters for the initial one-step PCR were reverse transcribed at 50°C for 30 minutes, initial denaturation at 94°C for 7 mins, 10 cycles consisting of a denaturation stage at 94°C for 10 seconds, annealing stage of 55°C for 30 seconds and an extension step of 68°C for 2 minutes and then 35 cycles with denaturation at 94°C for 10 seconds, annealing stage of 55.5°C for 30 seconds and final extension step of 68°C for 2 minutes, increasing each cycle by 10 seconds. The last stage was the final elongation at 68°C for 5 minutes with a hold stage at 4°C for a maximum of 18 hours. The nested PCR was performed using 1µL aliquot of first-round PCR product mixed with 9.0µL of RNase-free water, 12.5µL of Phusion High-Fidelity PCR Master Mix with HF Buffer, and 2.5µL of primers CWF1-LNA2(5’-GAA GGA CCA AAT GAA AGA YTG-3’ (2 μM) and RT-20C (5’CTG CCA ATT TCT AAC TGC CTT C-3 ‘(2 μM), making 25µL of reaction volume. The PCR parameters for the nested PCR were initial denaturation at 98°C for 30 seconds, 35 cycles made of denaturation at 98°C for 10 seconds, an annealing stage of 62°C for 30 seconds, an extension step of 72°C for 20 seconds, and the final elongation at 72°C for 10 minutes with a hold stage at 4°C for a maximum of 18 hours. Amplification was confirmed by electrophoresis in a 1% agarose gel in Tris-Borate-EDTA (TBE) buffer and ran at 90 volts for 45 minutes stained with 5μL ethidium bromide (0.5 mg/mL) and visualized under a UV source (260 nm). Failed samples were re-extracted and PCR amplification was repeated once using a set of rescue primers (Table 1) with the same PCR parameters and conditions.

Table 1 Rescue Primers Used When Samples Failed to Amplify

Sequencing

Successfully amplified HIV pol was purified using Exo-CIP™ Rapid PCR Clean-up Kit (New England Biolabs, Ipswich, USA) according to the manufacturer’s instructions. The purified amplicons were sequenced using BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Foster City, USA) on a 3031xl genetic analyzer using primers CWF1-LNA2(5’-GAA GGA CCA AAT GAA AGA YTG-3’), CWCS2(5’-AGAACTCAAGACTTTTGGG-3’), CWCS3(5’-TGCTGGGTGCGGTATTC-3’), CWCS5(5’TGGTAAATTTGATGTCCAT-3’), Seq2.1(5’GGCCAGGAATTTTCTTCAGAGC-3’), Seq6(5’-CCATCCCTGTGGAAGCACATTA-3’) and RT-20C (5’CTG CCA ATT TCT AAC TGC CTT C-3’). The cycle sequencing reaction mix contained 4.8μL of RNase-free water, 3μL of Big Dye 5X sequencing buffer, 1μL BigDye terminator, 0.2 μL of 2 μM of each sequencing primer, and 1μL of the purified amplicon, to make a total reaction volume of 10μL. The cycle sequencing reaction conditions were as follows: 25 cycles at 96°C for 10 seconds, 50°C for 5 seconds, 60°C for 4 minutes and hold at 4°C. The BigDye XTerminator purification kit (Applied Biosystems, Foster City, USA) was used to purify the sequencing reaction by adding 10μL of the BigDye XTerminator and 45μL SAM solution to cycle sequencing products. The reaction plate was vortexed at 1800 rpm for 45 minutes and then centrifuged at 3000g for 3 minutes at room temperature. Next, 65μL of purified cycle sequencing products in a reaction plate were analyzed using an ABI 3031xl genetic analyzer (Applied Biosystems, Foster City, USA)

Bioinformatics Analysis

The quality and read length of sequences obtained from the ABI 3130XL Genetic Analyzer were assessed using Sequencer Version 5.0 by manually trimming the beginning and end of each sequence to remove ambiguous nucleotides. Only sequences with quality scores of ≥73% were assembled into a contiguous sequence (contig). Mixed bases or ambiguous nucleotides were adopted or confirmed with the sequences covering the same position in the contig. The same software was used to assemble multiple reads of each sequence into a single contig(consensus sequence). The targeted sequence length was 1.1 kb, covering the partial Pol region. Sequence imputation and quality control assessment was assessed using AliView v1.26. Reference-based multiple sequence alignment (MSA) was constructed using muscle v3.8.31 implemented in AliView v1.26 and HIV-1 reference strain sequence (HXB2). MSA was further used for downstream analyses including phylogenetic inferences based on maximum-likelihood (ML), clade assessment, and DRM analysis.

All generated sequences that passed the QC were included in the phylogenetic tree. A maximum-likelihood tree topology was inferred from the reference-based MSA in IQ-TREE2 using General Time Reverse plus F plus R4 (GTR+F+R4) as a best-fit model of nucleotide substitution. A total of 1000 bootstrap replicates were employed to infer support for branches in the resulting tree topology using Booster. Tree visuals and annotations were performed in FigTree v1.4.4. Posterior probabilities 0.90 and above were noted as statistically significant.

HIV Drug Resistance Analysis

The MSA was imported into the Stanford HIV Drug Resistance Database to determine the HIV drug resistance profiles of each of the consensus sequences. HIV pol region was analysed for known DRM associated with nucleoside reverse transcriptase inhibitors (NRTIs), non-nucleoside reverse transcriptase inhibitors (NNRTIs) and protease inhibitors (PIs) according to the Stanford University HIV Drug Resistance Database.27 In sequences generated from ART naïve individuals, Calibrated Population Resistance in the Stanford University HIV database (https://hivdb.stanford.edu/cpr/form/PRRT/) was used to identify known surveillance drug resistance mutations.

Statistical Analysis

Patients’ general characteristics among the two VL groups were compared using Wilcoxon rank-sum test for continuous variables and Fisher’s exact test for categorical variables. The success rate of genotyping was calculated using the number of successfully genotyped samples in each VL category over the total number in that category multiplied by 100%. Genotyping success rates at LLV and VL ≥1000 copies/mL were compared using a comparison of proportions test to determine any statistical difference in the performance of the assay between these two groups. The prevalence of mutations was estimated with 95% confidence intervals using the binomial exact method for each group and further compared between the two VL groups. All of the analyses were done using STATA version 15 and p-values < 0.05 were considered statistically significant.

Results

Baseline and Clinical Characteristics

Among 96 participants included, 72% were females and the median age for all participants was 33.5, IQR: of 26.4–41 years. About 67% were antiretroviral therapy-experienced and the median HIV VL was 982, IQR: 669–12,422 copies/mL. The study participants were stratified into two VL groups: LLV (VL: 401–999 copies/mL) and VL ≥ 1000 copies/mL. Between these 2 groups, there was a statistically significant difference in ART status, whereby 64% of the LLV group were ART naïve and all participants with VL ≥ 1000 copies/mL were on ART (Table 2).

Table 2 Baseline Demographics and Clinical Characteristics for Study Participants

Amplification and Sequencing Success Rate

The lowest VL that the assay was able to successfully genotype was 423 copies/mL. The overall genotyping success rate was 79% (Figure 1), 69% (66/96) of samples were successfully genotyped by the original protocol while 10% (10/96) were genotyped using the rescue protocol described in Table 1. Stratifying by the VL group, a total of 72% (36/50) of samples with VL 401–999 copies/mL were successfully genotyped, compared to 87% (40/46) at VL ≥ 1000 copies/mL (p=0.07). The median log VL of successfully genotyped samples was 3.1 VL log10 copies/mL, IQR: 2.9–4.2, while the median for the failed samples was 2.9 VL log10 copies/mL, IQR: 2.7–3.1 (Figure 2). The genotyping outcomes of the in-house genotyping assay were not associated with the viral load (p=0.2).

Figure 1 Genotyping success of in-house RNA-based HIV genotyping assay, overall and stratified by VL group.

Abbreviations: LLV, low-level viremia (VL 401–999 copies/mL); VL, viral load.

Figure 2 Association of genotyping outcomes with viral load measurement.

Abbreviations: Q1, first quartile; Q3, third quartile; VL, viral load.

Phylogenetic Analysis

Sequence quality control and clustering were performed using phylogenetic analysis. Using pairwise distance analysis, the 76 sequences generated by the assay were unique. The maximum likelihood phylogenetic tree was constructed from HIV-1C RT/PR sequences generated using our in-house genotyping assay (Figure 3). There was no clustering of sequences by their VL group.

Figure 3 Maximum likelihood phylogenetic tree constructed using 1000 bootstrap values showing the uniqueness of the sequences and that there is no clustering of sequences by their VL group. The light-blue dot represents standard HIV reference strain.

Abbreviation: HXB2, standard HIV reference strain.

Drug Resistance Patterns

Among the generated sequences, the overall prevalence of individuals with at least 1 major-associated DRM was 24% (18/76) (Table 3). By VL group, individuals with LLV who were both ART naïve and experienced had a 17% prevalence of DRM compared to 30% among those with VL ≥ 1000 copies/mL who were all on ART (p-value=0.18). In the LLV group, 69% (25/36) of successfully genotyped individuals were ART naïve of whom 2 had DRM (1 with only NNRTI -associated mutations K101E and K103N while another one had only NRTI-associated mutation M41L). The prevalence of DRMs were NNRTI (28%), NRTI (24%) and PI (0%)-associated resistance mutations. Only 22% (4/18) had NRTI-, 39% (7/18) had NNRTI-, and 78% (14/18) had both NNRTI- and NRTI-associated resistance mutations. The most predominant mutations detected were K103N (18%) and M184V (12%) in NNRTI- and NRTI-associated mutations, respectively (Figure 4). Both NRTI- and NNRTI-associated resistance mutations were stratified by the VL group (Figures 5 and 6).

Table 3 Prevalence of HIV DRM, Overall and by VL Groups

Figure 4 HIV DRM classified by NRTI and NNRTI resistance classes (N=76) using an in-house RNA-based HIV genotyping assay with detectable plasma HIV-1 RNA > 400 copies/mL.

Abbreviations: DRM, drug resistance mutations; NRTI, nucleotide reverse transcriptase inhibitors; NNRTI, non-nucleotide reverse transcriptase inhibitors.

Figure 5 NRTI-associated resistance mutations stratified by VL groups.

Abbreviations: DRM, drug resistance mutations; LLV, low-level viremia; NRTI, nucleoside reverse transcriptase inhibitors; VL, viral load.

Figure 6 NNRTI-associated resistance mutations are stratified by VL groups.

Abbreviations: DRM, Drug resistance mutations; LLV, low level viremia; NNRTI, non-nucleoside reverse transcriptase inhibitors; VL, viral load.

Discussion

Our in-house RNA-based HIV genotyping assay showed that routine HIV genotyping of low-level viremia can be performed with a success rate of 72% and with a lower limit of genotyping of 423 copies/mL. The assay was able to detect 17% of LLV participants with at least one drug resistance mutation.

The performance of the in-house genotyping assay was higher (87%) among samples with VL ≥ 1000 copies/mL, compared to the LLV group (72%), however this difference was not statistically significant (p=0.07). Our study findings highlight the potential of the assay to genotype samples with LLV and a similar trend was reported in other studies.10,16,28–30 It is difficult to compare our study with others on how a lower LLV range determines genotyping success due to the use of LLV 401–999 copies/mL when other studies went as low as 51 copies/mL. The same threshold of VL > 400 copies/mL was utilized in a South African cohort where the overall genotyping success rate was 83%,31 which was not statistically different from the 79% reported in our study (p=0.13).

Although our study had a modest sample size, the genotyping success was similar to 72.4% among 1072 samples28 and 67% among 756 samples, in two studies done in France.32 However, our genotyping success rate was relatively lower compared to other assays utilized in LLV settings with success ≥ 95%.10,15,16,30,33,34 A higher plasma input volume of ≥ 1mL is usually recommended for LLV samples, to increase HIV copies for amplification;35–38 however, our assay was able to successfully genotype about 72% of LLV samples using a plasma input volume of 200 µL for extraction. Nevertheless, this is the first LLV genotyping assay to have a success rate of 72% with non-concentrated virus copies using lower sample input volume.

Current HIV treatment guidelines on genotyping of detectable VL below 1000 copies/mL testing are unclear because it might be unsuccessful.39 Furthermore, these guidelines recommend strongly against testing for patients with VL 200–1000 copies/mL using RNA-based approaches, but rather recommend using proviral DNA. Additionally, resistance assay kits are only validated to test samples with VL ≥ 100040,41 or VL ≥ 200042,43 copies/mL. However, our results indicate that resistance testing of samples with viral loads below 1000 copies/mL provides clinically relevant information.

The association between VL and the overall prevalence of DRM was statistically similar when compared among individuals with LLV who were both ART naïve and experienced (17%) against those with VL ≥ 1000 copies/mL (30%) who were all on ART. However, a statistically higher prevalence of NNRTI-associated resistance was reported among individuals with VL ≥ 1000 copies/mL compared to LLV. The overall prevalence of DRM among individuals with LLV included 69% (25/36) of successfully genotyped individuals who were ART naïve, which might have lowered the NNRTI-associated resistance in this group. The main objective was to assess the practicability of genotyping LLV samples. Therefore, both ART-naïve and -experienced individuals were included in the study. We report no statistically significant difference in the prevalence of NNRTI-associated resistance mutations between 2 VL groups when ART-experienced LLV was compared to VL ≥ 1000 copies/mL. These findings were concordant with our previous study that reported no association between VL and DRM in the same cohort, although the DRMs were detected in proviral DNA.44 Most studies have reported concordance in DRM detected in both viral RNA and proviral DNA;45,46 therefore both can be useful in LLV for the detection of DRM.

The assay was not assessed on non-HIV-1C viruses because the HIV-1C epidemic is predominant in sub-Saharan Africa, which may make the assay applicable in the region. One of the study limitations is that the performance of the optimized assay was not compared to any commercial HIV genotyping assay, however, most of the commercial genotyping assays are validated for VL measurements of ≥1000 copies/mL which makes the comparison difficult in LLV samples. One of the study limitations is that most of the ART-experienced individuals have missing ART regimen information. In BCPP cohort, 60% of PWH were on efavirenz (EFV), 26.6% on nevirapine (NVP) and 7% on dolutegravir (DTG)-based regimens while 5% were on second-line therapy and 1% on salvage therapy among individuals who had available ARV regimen information. Those who were on EFV and NVP, were also on 2 NRTI regimens such as lamivudine (3TC) and zidovudine (ZDV) or tenofovir disoproxil fumarate (TDF) and emtricitabine (FTC). In this study, 64 participants were on ART, where ART regimen was available for 17 individuals (8 – EFV, 4-NVP, 3-lopinavir (LPV) and 1-DTG). Despite the majority of individuals missing ART regimen information, the higher prevalence of NNRTI and NRTI-associated DRM reported in our study supports the exposure of NNRTI and NRTI-containing regimens among these individuals. However, our findings confirmed that individuals with LLV in the settings of HIV-1C can be genotyped, and harbor DRM as reported in other studies. Taken together, these results highlight the concept that genotyping resistance testing may be useful in managing individuals with LLV as resistance data can be used for selecting appropriate ART, especially for those with persisting LLV.

In conclusion, an in-house RNA-based HIV genotyping assay was able to successfully genotype LLV samples (VL 401–999 copies/mL) with the lowest viral load of 423 copies/mL and predicted 17% DRM amongst them. However, the assay can be further optimized to increase the genotyping success rate and for lower VL, and further evaluated for reproducibility and specificity and compared with a commercial genotyping assay.

Data Sharing Statement

All relevant data are within the paper, including the figures and tables. The sequences obtained in our study were submitted to National Center for Biotechnology Information (NCBI) GenBank and the accession numbers are OP922678 – OP922753.

Ethics Approval

The BCPP study was approved by the Institutional Review Boards (IRBs) at the US Centers for Disease Control and Prevention and the Botswana Health Research and Development Committee and is registered at ClinicalTrials.gov (NCT01965470). The study was conducted according to the principle expressed in the declaration of Helsinki. All BCPP participants provided written informed consent for their samples to be used in future studies.

Acknowledgments

We thank all the study participants. We thank the BCPP study team for their contribution to this study. We thank the Ministry of Health and Wellness, Botswana Harvard AIDS Institute Partnership, CDC Botswana, and US CDC for their excellent support and contributions to the BCPP study. We also thank Katlego Selelo, Pearl Kaumba and Charlene Raphaka for their assistance with the laboratory work and Lendsey Melton for his valuable edits.

Author Contributions

All authors made a significant contribution to the work reported, whether that is in the conception, study design, execution, acquisition of data, analysis and interpretation, or in all these areas; took part in drafting, revising or critically reviewing the article; gave final approval of the version to be published; have agreed on the journal to which the article has been submitted; and agree to be accountable for all aspects of the work.

Funding

The BCPP Impact Evaluation was funded by the President’s Emergency Plan for AIDS Relief (PEPFAR) through the Centers for Disease Control and Prevention (CDC, cooperative agreements U01 GH000447 and U2G GH001911). OTB was supported by the Fogarty International Center at the US National Institutes of Health (D43 TW009610). SM is partially supported by the Fogarty International Center at the US National Institutes of Health (1 K43 TW 012350-01). SG and WTC were partially supported by H3ABioNet. H3ABioNet is supported by the National Institutes of Health Common Fund (U41 HG006941). H3ABioNet is an initiative of the Human Health and Heredity in Africa Consortium (H3Africa) program of the African Academy of Science (AAS). SM, RM, DD & OTB were supported by the Trials of Excellence in Southern Africa (TESA III), which is part of the EDCTP2 program supported by the European Union (CSA2020NoE-3104 TESAIII CSA2020NoE). SG, SM, WTC and NOM are partly supported through the Sub-Saharan African Network for TB/HIV Research Excellence (SANTHE 2.0) from the Bill and Melinda Gates Foundation (INV-033558). PANGEA-HIV is funded primarily by the Bill and Melinda Gates Foundation. The views expressed in this publication are those of the authors and not necessarily those of the funding agencies. The funders had no role in the study design, data collection, decision to publish, or preparation of this manuscript.

Disclosure

The authors report no conflicts of interest in relation to this work.

References

1. Sklar PA, Ward DJ, Baker RK., et al. Prevalence and clinical correlates of HIV viremia (‘blips’) in patients with previous suppression below the limits of quantification. Aids. 2002;16(15):2035–2041. doi:10.1097/00002030-200210180-00008

2. Tobin NH, Learn GH, Holte SE, et al. Evidence that low-level viremias during effective highly active antiretroviral therapy result from two processes: expression of archival virus and replication of virus. J Virol. 2005;79(15):9625–9634. doi:10.1128/jvi.79.15.9625-9634.2005

3. Bessong PO, Matume ND, Tebit DM. Potential challenges to sustained viral load suppression in the HIV treatment programme in South Africa: a narrative overview. AIDS Res Ther. 2021;18(1):1. doi:10.1186/s12981-020-00324-w

4. Thaker HK. HIV viral suppression in the era of antiretroviral therapy. Postgrad Med J. 2003;79(927):36. doi:10.1136/pmj.79.927.36

5. Jordan MR, Winsett J, Tiro A, et al. HIV Drug Resistance Profiles and Clinical Outcomes in Patients with Viremia Maintained at Very Low Levels. World Journal of AIDS. 2013;3(2):71–78. doi:10.4236/wja.2013.32010

6. Li JZ, Gallien S, Do TD, et al. Prevalence and significance of HIV-1 drug resistance mutations among patients on antiretroviral therapy with detectable low-level viremia. Antimicrob Agents Chemother. 2012;56(11):5998–6000. doi:10.1128/aac.01217-12

7. Swenson LC, Min JE, Woods CK, et al. HIV drug resistance detected during low-level viraemia is associated with subsequent virologic failure. AIDS. 2014;28(8):1125–1134. doi:10.1097/qad.0000000000000203

8. Inzaule SC, Bertagnolio S, Kityo CM, et al. The relative contributions of HIV drug resistance, nonadherence and low-level viremia to viremic episodes on antiretroviral therapy in sub-Saharan Africa. AIDS. 2020;34(10):1559–1566. doi:10.1097/QAD.0000000000002588

9. Rosemary A, Chika O, Jonathan O, et al. Genotyping performance evaluation of commercially available HIV-1 drug resistance test. PLoS One. 2018;13(6):e0198246. doi:10.1371/journal.pone.0198246

10. Santoro MM, Fabeni L, Armenia D, et al. Reliability and clinical relevance of the HIV-1 drug resistance test in patients with low viremia levels. Clin Infect Dis. 2014;58(8):1156–1164. doi:10.1093/cid/ciu020

11. Villalobos C, Ceballos ME, Ferres M, Palma C. Drug resistance mutations in proviral DNA of HIV-infected patients with low level of viremia. J Clin Virol. 2020;132:104657. doi:10.1016/j.jcv.2020.104657

12. Panel on Antiretroviral Guidelines for Adults and Adolescents. Guidelines for the Use of Antiretroviral Agents in Adults and Adolescents with HIV. US Department of Human and Health Services; 2022.

13. Novitsky V, Zahralban-Steele M, McLane MF, et al. Long-Range HIV Genotyping Using Viral RNA and Proviral DNA for Analysis of HIV Drug Resistance and HIV Clustering. J Clin Microbiol. 2016;4:1175. doi:10.1128/JCM.00190-16

14. Zash R, Holmes L, Diseko M, et al. Neural-Tube Defects and Antiretroviral Treatment Regimens in Botswana. N Engl J Med. 2019;381(9):827–840. doi:10.1056/NEJMoa1905230

15. Yang Z, Morrison R, Oates C, et al. HIV-1 Genotypic Resistance Testing on Low Viral Load Specimens Using the Abbott ViroSeq HIV-1 Genotyping System. Lab Med. 2008;39(11):671–673. doi:10.1309/lmvovu1xrb9o3jzv

16. Gonzalez-Serna A, Min JE, Woods C, et al. Performance of HIV-1 drug resistance testing at low-level viremia and its ability to predict future virologic outcomes and viral evolution in treatment-naive individuals. Clin Infect Dis. 2014;58(8):1165–1173. doi:10.1093/cid/ciu019

17. Richman DD. Editorial commentary: extending HIV drug resistance testing to low levels of plasma viremia. Clin Infect Dis. 2014;58(8):1174–1175. doi:10.1093/cid/ciu025

18. AIDS JUNPoH. IN DANGER: UNAIDS Global AIDS Update 2022; 2022.

19. Gaolathe T, Wirth KE, Holme MP, et al. Botswana’s progress toward achieving the 2020 UNAIDS 90- 90-90 antiretroviral therapy and virological suppression goals: a population-based survey. Lancet HIV. 2016;3(5):e221–30. doi:10.1016/s2352-3018(16

20. Stafford MM, Laws K, Marima RL. PELBC01. Botswana achieved the Joint United Nations Programme on HIV/AIDS (UNAIDS) 95‐95‐95 targets: results from the Fifth Botswana HIV/AIDS Impact Survey (BAIS V). J Int AIDS Soc. 2022;2022:231–232.

21. Society EAC. EACS Guidelines; 2018.

22. services DoHaH. Guidelines for the Use of Antiretroviral Agents in Adults and Adolescents Living with HIV; 2019.

23. WHO. Guidelines for Managing Advanced HIV Disease and Rapid Initiation of Antiretroviral Therapy; 2017.

24. Rowley CF, MacLeod IJ, Maruapula D, et al. Sharp increase in rates of HIV transmitted drug resistance at antenatal clinics in Botswana demonstrates the need for routine surveillance. J Antimicrob Chemother. 2016;71(5):1361–1366. doi:10.1093/jac/dkv500

25. Makhema J, Wirth KE, Pretorius Holme M, et al. Universal Testing, Expanded Treatment, and Incidence of HIV Infection in Botswana. N Engl J Med. 2019;381(3):230–242. doi:10.1056/NEJMoa1812281

26. World Health Organization. Consolidated Guidelines on the Use of Antiretroviral Drugs for Treating and Preventing HIV Infection: Recommendations for a Public Health Approach. 2nd. World Health Organization; 2016.

27. Rhee SY. Human immunodeficiency virus reverse transcriptase and protease sequence database. Nucleic Acids Res. 2003;31(1):298–303. doi:10.1093/nar/gkg100

28. Assoumou L, Charpentier C, Recordon-Pinson P, et al. Prevalence of HIV-1 drug resistance in treated patients with viral load >50 copies/mL: a 2014 French nationwide study. J Antimicrob Chemother. 2017;72(6):1769–1773. doi:10.1093/jac/dkx042

29. Parczewski M, Leszczyszyn-Pynka M, Witak-Jędra M, Maciejewska K, Urbańska A. Efficacy of genotypic drug resistance testing in patients with low-level plasma HIV-1 viremia. HIV AIDS Rev. 2015;14(3):80–83. doi:10.1016/j.hivar.2015.02.001

30. Mackie N, Dustan S, McClure MO, Weber JN, Clarke JR. Detection of HIV-1 antiretroviral resistance from patients with persistently low but detectable viraemia. J Virol Methods. 2004;119(2):73–78. doi:10.1016/j.jviromet.2004.02.015

31. Bangalee A, Hans L, Steegen K. Feasibility and clinical relevance of HIV-1 drug resistance testing in patients with low-level viraemia in South Africa. J Antimicrob Chemother. 2021;76(10):2659–2665. doi:10.1093/jac/dkab220

32. Assoumou L, Descamps D, Yerly S, et al. Prevalence of HIV-1 drug resistance in treated patients with viral load >50 copies/mL in 2009: a French nationwide study. J Antimicrob Chemother. 2013;68(6):1400–1405. doi:10.1093/jac/dkt033

33. Waters L, Mandalia S, Asboe D. Successful use of genotypic resistance testing in HIV-1-infected individuals with detectable viraemia between 50 and 1000 copies/mL. AIDS. 2006;20(5):778–779. doi:10.1097/01.aids.0000216381.37679.69

34. Cane PA, Kaye S, Smit E, et al. Genotypic antiretroviral drug resistance testing at low viral loads in the UK. HIV Med. 2008;9(8):673–676. doi:10.1111/j.1468-1293.2008.00607.x

35. Mens H, Kearney M, Wiegand A, et al. Amplifying and quantifying HIV-1 RNA in HIV infected individuals with viral loads below the limit of detection by standard clinical assays. J Vis Exp. 2011;(55):2960. doi:10.3791/2960

36. Hatano H, Delwart EL, Norris PJ, et al. Evidence of persistent low-level viremia in long-term HAART-suppressed, HIV-infected individuals. AIDS. 2010;24(16):2535–2539. doi:10.1097/QAD.0b013e32833dba03

37. Vancoillie L, Mortier V, Demecheleer E, et al. Drug resistance is rarely the cause or consequence of long-term persistent low-level viraemia in HIV-1-infected patients on ART. Antivir Ther. 2015;20(8):789–794. doi:10.3851/imp2966

38. Villahermosa ML, Thomson M, Vázquez de Parga E, et al. Improved conditions for extraction and amplification of human immunodeficiency virus type 1 RNA from plasma samples with low viral load. J Hum Virol. 2000;3(1):27–34.

39. Günthard HF, Calvez V, Paredes R, et al. Human Immunodeficiency Virus Drug Resistance: 2018 Recommendations of the International Antiviral Society–USA Panel. Clin Infect Dis. 2019;68(2):177–187. doi:10.1093/cid/ciy463

40. Dorward J, Sookrajh Y, Ngobese H, et al. Protocol for a randomised feasibility study of Point-Of-care HIV viral load testing to Enhance Re-suppression in South Africa: the POwER study. BMJ Open. 2021;11(2):e045373. doi:10.1136/bmjopen-2020-045373

41. Hirsch MS, Günthard HF, Schapiro JM, et al. Antiretroviral Drug Resistance Testing in Adult HIV-1 Infection: 2008 Recommendations of an International AIDS Society–USA Panel. Clin Infect Dis. 2008;47(2):266–285. doi:10.1086/589297

42. Cunningham S, Ank B, Lewis D, et al. Performance of the applied biosystems ViroSeq human immunodeficiency virus type 1 (HIV-1) genotyping system for sequence-based analysis of HIV-1 in pediatric plasma samples. J Clin Microbiol. 2001;39(4):1254–1257. doi:10.1128/jcm.39.4.1254-1257.2001

43. Mukaide M, Sugiura W, Matuda M, et al. Evaluation of Viroseq-HIV version 2 for HIV drug resistance. Jpn J Infect Dis. 2000;53(5):203–205.

44. Bareng OT, Moyo S, Zahralban-Steele M, et al. HIV-1 drug resistance mutations among individuals with low-level viraemia while taking combination ART in Botswana. J Antimicrob Chemother. 2022;77(5):1385–1395. doi:10.1093/jac/dkac056

45. Moyo S, Gaseitsiwe S, Zahralban-Steele M, et al. Low rates of nucleoside reverse transcriptase inhibitor and nonnucleoside reverse transcriptase inhibitor drug resistance in Botswana. Aids. 2019;33(6):1073–1082. doi:10.1097/qad.0000000000002166

46. Bon I, Alessandrini F, Borderi M, Gorini R, Re MC. Analysis of HIV-1 drug-resistant variants in plasma and peripheral blood mononuclear cells from untreated individuals: implications for clinical management. New Microbiol. 2007;30(3):313–317.

Creative Commons License © 2022 The Author(s). This work is published by Dove Medical Press Limited, and licensed under a Creative Commons Attribution License. The full terms of the License are available at http://creativecommons.org/licenses/by/4.0/. The license permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.