High prevalence of pre-treatment and acquired HIV-1 drug resistance mutations among non-citizens living with HIV in Botswana.

Background: Approximately 30,000 non-citizens are living with HIV in Botswana, all of whom as of 2020 are eligible to receive free antiretroviral treatment (ART) within the country. We assessed the prevalence of HIV-1 mutational profiles [pre-treatment drug resistance (PDR) and acquired drug resistance (ADR)] among treatment-experienced (TE) and treatment-naïve (TN) non-citizens living with HIV in Botswana. Methods: A total of 152 non-citizens living with HIV were enrolled from a migrant HIV clinic at Independence Surgery, a private practice in Botswana from 2019-2021. Viral RNA isolated from plasma samples were genotyped for HIV drug resistance (HIVDR) using Sanger sequencing. Major known HIV drug resistance mutations (DRMs) in the pol region were determined using the Stanford HIV Drug Resistance Database. The proportions of HIV DRMs amongst TE and TN non-citizens were estimated with 95% confidence intervals (95% CI) and compared between the two groups. Results: A total of 60/152 (39.5%) participants had a detectable viral load (VL) >40 copies/mL and these were included in the subsequent analyses. The median age at enrollment was 43 years (Q1, Q3: 38-48). Among individuals with VL > 40 copies/mL, 60% (36/60) were treatment-experienced with 53% (19/36) of them on Atripla. Genotyping had a 62% (37/60) success rate - 24 were TE, and 13 were TN. A total of 29 participants (78.4, 95% CI: 0.12-0.35) had major HIV DRMs, including at least one non-nucleoside reverse transcriptase inhibitor (NNRTI) associated DRM. In TE individuals, ADR to any antiretroviral drug was 83.3% (20/24), while for PDR was 69.2% (9/13). The most frequent DRMs were nucleoside reverse transcriptase inhibitors (NRTIs) M184V (62.1%, 18/29), NNRTIs V106M (41.4%, 12/29), and K103N (34.4%, 10/29). No integrase strand transfer inhibitor-associated DRMs were reported. Conclusion: We report high rates of PDR and ADR in ART-experienced and ART-naïve non-citizens, respectively, in Botswana. Given the uncertainty of time of HIV acquisition and treatment adherence levels in this population, routine HIV-1C VL monitoring coupled with HIVDR genotyping is crucial for long-term ART success.

Rapid dynamic changes of FL.2 variant: A case report of COVID-19 breakthrough infection.

We investigated intra-host genetic evolution using two SARS-CoV-2 isolates from a fully vaccinated (primary schedule x2 doses of AstraZeneca plus a booster of Pfizer), >70-year-old woman with a history of lymphoma and hypertension who presented a SARS-CoV-2 infection for 3 weeks prior to death due to COVID-19. Two full genome sequences were determined from samples taken 13 days apart with both belonging to Pango lineage FL.2: the first detection of this Omicron sub-variant in Botswana. FL.2 is a sub-lineage of XBB.1.9.1. The repertoire of mutations and minority variants in the Spike protein differed between the two time points. Notably, we also observed deletions within the ORF1a and Membrane proteins; both regions are associated with high T-cell epitope density. The internal milieu of immune-suppressed individuals may accelerate SARS-CoV-2 evolution; hence, close monitoring is warranted.

Predicted resistance to broadly neutralizing antibodies (bnAbs) and associated HIV-1 envelope characteristics among seroconverting adults in Botswana.

We used HIV-1C sequences to predict (in silico) resistance to 33 known broadly neutralizing antibodies (bnAbs) and evaluate the different HIV-1 Env characteristics that may affect virus neutralization. We analyzed proviral sequences from adults with documented HIV-1 seroconversion (N = 140) in Botswana (2013-2018). HIV-1 env sequences were used to predict bnAb resistance using bNAb-ReP, to determine the number of potential N-linked glycosylation sites (PNGS) and evaluate Env variable region characteristics (VC). We also assessed the presence of signature mutations that may affect bnAb sensitivity in vitro. We observe varied results for predicted bnAb resistance among our cohort. 3BNC117 showed high predicted resistance (72%) compared to intermediate levels of resistance to VRC01 (57%). We predict low resistance to PGDM100 and 10-1074 and no resistance to 4E10. No difference was observed in the frequency of PNGS by bNAb susceptibility patterns except for higher number of PNGs in V3 bnAb resistant strains. Associations of VC were observed for V1, V4 and V5 loop length and net charge. We also observed few mutations that have been reported to confer bnAb resistance in vitro. Our results support use of sequence data and machine learning tools to predict the best bnAbs to use within populations.

Fostemsavir resistance-associated polymorphisms in HIV-1 subtype C in a large cohort of treatment-naïve and treatment-experienced individuals in Botswana.

IMPORTANCE: Fostemsavir (FTR) is a newly licensed antiretroviral drug that has been shown to have activity against HIV-1. The mechanism of action of FTR is different from all currently available antiretrovirals (ARVs), and as such, it offers hope for HIV-1 suppression in those people with HIV (PWH) who harbor HIV-1 variants with drug resistance mutations to currently used ARVs. Using 6,030 HIV-1 sequences covering the HIV-1 envelope from PWH in Botswana who are antiretroviral therapy (ART) naïve as well as those who are failing ART, we explored the sequences for FTR resistance-associated polymorphisms. We found the prevalence of FTR polymorphisms to be similar in both ART-naïve and ART-experienced individuals with VF in this setting, with no prior FTR exposure. Further studies on the phenotypic impact of these polymorphisms are warranted to guide how to monitor for FTR resistance.

Predicted broadly neutralizing antibody (bnAb) resistance and associated envelope characteristics of adults with HIV-1 seroconversion in Botswana.

We used HIV-1C sequences to predict (in silico) resistance to 33 known broadly neutralizing antibodies (bNAbs) and evaluate the different HIV-1 env characteristics that may affect virus neutralization. We analyzed proviral sequences from adults with documented HIV-1 seroconversion (N=140) in Botswana (2013-2018). HIV-1 env sequences were used to predict bnAb resistance using bNAb-ReP, to determine the number of potential N-linked glycosylation sites (PNGS) and evaluate env variable region characteristics (VC). We also assessed the presence of signature mutations that may affect bnAb sensitivity in vitro. We observe varied results for predicted bnAb resistance among our cohort. 3BNC117 showed high predicted resistance (72%) compared to intermediate levels of resistance to VRC01 (57%). We predict low resistance to PGDM100 and 10-1074 and no resistance to 4E10. No difference was observed in the frequency of PNGS by bNAb susceptibility patterns except for higher number of PNGs in V3 bnAb resistant strains. Associations of VC were observed for V1, V4 and V5 loop length and net charge. We also observed few mutations that have been reported to confer bnAb resistance in vitro. Our results support use of sequence data and machine learning tools to predict the best bnAbs to use within populations.

Archived rilpivirine-associated resistance mutations among ART-naive and virologically suppressed people living with HIV-1 subtype C in Botswana: implications for cabotegravir/rilpivirine use.

OBJECTIVES: Pre-existing rilpivirine resistance-associated mutations (RVP-RAMs) have been found to predict HIV-1 virological failure in those switching to long-acting injectable cabotegravir/rilpivirine. We here evaluated the prevalence of archived RPV-RAMs in a cohort of people living with HIV (PWH). METHODS: We analysed near full-length HIV-1 pol sequences from proviral DNA for the presence of RPV-RAMs, which were defined according to the 2022 IAS-USA drug resistance mutation list and Stanford HIV drug resistance database. RESULTS: RPV-RAMs were identified in 757/5805 sequences, giving a prevalence of 13.0% (95% CI 12%-13.9%). Amongst the ART-naive group, 137/1281 (10.7%, 95% CI 9.1%-12.5%) had at least one RPV-RAM. Of the 4524 PWH with viral suppression on ART (VL <400 copies/mL), 620 (13.7%, 95% CI 12.7%-14.7%) had at least one RPV-RAM. E138A was the most prevalent RPV-RAM in the ART-naive group (7.9%) and the ART-suppressed group (9.3%). The rest of the mutations observed (L100I, K101E, E138G, E138K, E138Q, Y181C, H221Y, M230L, A98G, V179D, G190A, G190E and M230I) were below a prevalence of 1%. CONCLUSIONS: RPV-RAMs were present in 10.7% of ART-naive and 13.7% of ART-suppressed PWH in Botswana. The most common RPV-RAM in both groups was E138A. Since individuals with the E138A mutation may be more likely to fail cabotegravir/rilpivirine, monitoring RPV-RAMs will be crucial for effective cabotegravir/rilpivirine implementation in this setting.

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

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.

Doravirine-associated resistance mutations in antiretroviral therapy naïve and experienced adults with HIV-1 subtype C infection in Botswana.

OBJECTIVES: There are limited data on the prevalence of doravirine (DOR)-associated drug resistance mutations in people with HIV (PWH) in Botswana. This cross-sectional, retrospective study aimed to explore the prevalence of DOR-associated resistance mutations among ART-naïve and -experienced PWH in Botswana enrolled in the population-based Botswana Combination Prevention Project (BCPP). METHODS: A total of 6078 HIV-1C pol sequences were analysed for DOR-associated resistance mutations using the Stanford HIV drug resistance database, and their levels were predicted according to the Stanford DRM penalty scores and resistance interpretation. Virologic failure was defined as HIV-1 RNA load (VL) >400 copies/mL. RESULTS: Among 6078 PWH, 5999 (99%) had known ART status, and 4529/5999 (79%) were on ART at time of sampling. The suppression rate among ART-experienced was 4517/4729 (96%). The overall prevalence of any DOR-associated resistance mutations was 181/1473 (12.3% [95% confidence interval {CI}: 10.7-14.1]); by ART status: 42/212 (19.8% [95% CI: 14.7-25.4]) among ART-failing individuals (VL ≥400 copies/mL) and 139/1261 (11.0% [95% CI: 9.3-12.9]) among ART-naïve individuals (P < 0.01). Intermediate DOR-associated resistance mutations were observed in 106/1261 (7.8% [95% CI: 6.9-10.1]) in ART-naïve individuals and 29/212 (13.7% [95% CI: 9.4-8.5]) among ART-experienced participants (P < 0.01). High-level DOR-associated resistance mutations were observed in 33/1261 (2.6% [95% CI: 1.8-3.7]) among ART-naïve and 13/212 (6.1% [95% CI: 3.6-10.8]) among ART-failing PWH (P < 0.01). PWH failing ART with at least one EFV/NVP-associated resistance mutation had high prevalence 13/67 (19.4%) of high-level DOR-associated resistance mutations. CONCLUSION: DOR-associated mutations were rare (11.0%) among ART-naive PWH but present in 62.7% of Botswana individuals who failed NNRTI-based ART with at least one EFV/NVP-associated resistance mutation. Testing for HIV drug resistance should underpin the use of DOR in PWH who have taken first-generation NNRTIs.

Reversal of CSF HIV-1 Escape during Treatment of HIV-Associated Cryptococcal Meningitis in Botswana.

Cerebrospinal fluid (CSF) viral escape has been poorly described among people with HIV-associated cryptococcal meningitis. We determined the prevalence of CSF viral escape and HIV-1 viral load (VL) trajectories in individuals treated for HIV-associated cryptococcal meningitis. A retrospective longitudinal study was performed using paired CSF and plasma collected prior to and during the antifungal treatment of 83 participants recruited at the Botswana site of the phase-3 AMBITION-cm trial (2018−2021). HIV-1 RNA levels were quantified then CSF viral escape (CSF HIV-1 RNA ≥ 0.5 log10 higher than plasma) and HIV-1 VL trajectories were assessed. CSF viral escape occurred in 20/62 (32.3%; 95% confidence interval [CI]: 21.9−44.6%), 13/52 (25.0%; 95% CI: 15.2−38.2%) and 1/33 (3.0%; 95% CI: 0.16−15.3%) participants at days 1, 7 and 14 respectively. CSF viral escape was significantly lower on day 14 compared to days 1 and 7, p = 0.003 and p = 0.02, respectively. HIV-1 VL decreased significantly from day 1 to day 14 post antifungal therapy in the CSF but not in the plasma (ß = −0.47; 95% CI: −0.69 to −0.25; p < 0.001). CSF viral escape is high among individuals presenting with HIV-associated cryptococcal meningitis; however, antifungal therapy may reverse this, highlighting the importance of rapid initiation of antifungal therapy in these patients.

HIV-1 drug resistance mutations among individuals with low-level viraemia while taking combination ART in Botswana.

OBJECTIVES: To assess whether a single instance of low-level viraemia (LLV) is associated with the presence of drug resistance mutations (DRMs) and predicts subsequent virological failure (VF) in adults receiving ART in 30 communities participating in the Botswana Combination Prevention Project. METHODS: A total of 6078 HIV-1 C pol sequences were generated and analysed using the Stanford HIV drug resistance database. LLV was defined as plasma VL = 51-999 copies/mL and VF was defined as plasma VL ≥ 1000 copies/mL. RESULTS: Among 6078 people with HIV (PWH), 4443 (73%) were on ART for at least 6 months. Of the 332 persons on ART with VL > 50 copies/mL, 175 (4%) had VL ≥ 1000 copies/mL and 157 (4%) had LLV at baseline. The prevalence of any DRM was 57 (36%) and 78 (45%) in persons with LLV and VL ≥ 1000 copies/mL, respectively. Major DRMs were found in 31 (20%) with LLV and 53 (30%) with VL ≥ 1000 copies/mL (P = 0.04). Among the 135 PWH with at least one DRM, 17% had NRTI-, 35% NNRTI-, 6% PI- and 3% INSTI-associated mutations. Among the 3596 participants who were followed up, 1709 (48%) were on ART for ≥6 months at entry and had at least one subsequent VL measurement (median 29 months), 43 (3%) of whom had LLV. The OR of experiencing VF in persons with LLV at entry was 36-fold higher than in the virally suppressed group. CONCLUSIONS: A single LLV measurement while on ART strongly predicted the risk of future VF, suggesting the use of VL > 50 copies/mL as an indication for more intensive adherence support with more frequent VL monitoring.

Genetic diversity in L1 ORF of human papillomavirus in women with cervical cancer with and without human immunodeficiency virus in Botswana and Kenya.

BACKGROUND: The variation of human papillomavirus (HPV) genotypes shapes the risks of cervical cancer and these variations are not well defined in Africa. Nucleotide changes within the L1 gene, nucleotide variability, and phylogeny were explored in relation to HIV in samples from Botswana and Kenya. METHODS: A total of 98 HPV-positive cervical samples were sequenced to identify different HPV variants. Phylogenetic inferences were used to determine HPV genotypes and investigate the clustering of sequences between women living with HIV (WLWHIV) and -women not living with HIV (WNLWHIV). RESULTS: Out of 98 generated sequences, 83.7% (82/98) participants had high-risk (HR) HPV genotypes while 16.3% (16/98) had low-risk (LR) HPV genotypes. Among participants with HR-HPV genotypes, 47.6% (39/82) were coinfected with HIV. The prevalence of HR-HPV genotypes was statistically higher in the Botswana population compared to Kenya (p-value < 0.001). Multiple amino acid mutations were identified in both countries. Genetic diversity differed considerably among WLWHIV and WNLWHIV. The mean pairwise distances between HPV-16 between HIV and HIV/HPV as well as for HPV-18 were statistically significant. Six (6) new deleterious mutations were identified in the HPV genotypes based on the sequencing of the L1 region, HPV-16 (L441P, S343P), HPV-18 (S424P), HPV-45 (Q366H, Y365F), and HPV-84 (F458L). The majority of the patients with these mutations were co-infected with HIV. CONCLUSIONS: Genomic diversity and different genomic variants of HPV sequences were demonstrated. Candidate novel mutations within the L1 gene were identified in both countries which can be further investigated using functional assays.

Patterns of pretreatment drug resistance mutations of very early diagnosed and treated infants in Botswana.

OBJECTIVE: To describe the occurrence of HIV drug resistance mutations (DRMs) in both intact and defective HIV-1 cell-associated DNA (HIV-1 CAD) among early-treated infants. DESIGN: The Botswana EIT Study (ClinicalTrials.gov NCT02369406) initiated antiretroviral therapy (ART) in the first week of life and evaluated HIV-1 in plasma and peripheral blood mononuclear cells (PBMCs). METHODOLOGY: We analyzed 257 near-HIV-1 full-length sequences (nFLS) obtained by Illumina next-generation sequencing from infants near birth. Sanger sequencing of pol was performed for mothers at delivery and children with clinical failure through 96 weeks. DRMs were identified using the Stanford HIV Drug Resistance Database. RESULTS: In 27 infants, median PBMC HIV-1 proviral load was 492 copies/ml [IQR: 78-1246 copies/ml] at a median of 2 days (range 1-32); 18 (66.7%) had no DRMs detected; six (22.2%) had DRMs detected in defective DNA only, and three (11.1%) had DRMs in both defective and intact DNA (P = 0.09). A total of 60 of 151 (37.7%) defective sequences had at least one DRM: 31.8% NNRTI, 15.2% NRTI, 5.3% protease inhibitor, and 15.5% INSTI-associated mutations. In intact sequences, 33 of 106 (31.1%) had at least 1 DRM: 29.2% NNRTI, 7.5% NRTI, 0.9% protease inhibitor, and no INSTI-associated mutations. For all three infants with intact sequence DRMs, corresponding DRMs occurred in maternal plasma at delivery. Archived DRMs were detectable at a later clinical rebound on only one occasion. CONCLUSION: Defective HIV-1 cell-associated DNA sequences may overestimate the prevalence of drug resistance among early-treated children. The impact of DRMs from intact proviruses on long-term treatment outcomes warrants further investigation.

Low rates of nucleoside reverse transcriptase inhibitor and nonnucleoside reverse transcriptase inhibitor drug resistance in Botswana.

BACKGROUND: Scale-up of antiretroviral therapy (ART) and introduction of treat-all strategy necessitates population-level monitoring of acquired HIV drug resistance (ADR) and pretreatment drug resistance (PDR) mutations. METHODS: Blood samples were collected from 4973 HIV-positive individuals residing in 30 communities across Botswana who participated in the Botswana Combination Prevention Project (BCPP) in 2013-2018. HIV sequences were obtained by long-range HIV genotyping. Major drug-resistance mutations (DRMs) and surveillance drug resistance mutations (SDRMs) associated with nucleoside reverse transcriptase inhibitors (NRTI) and nonnucleoside reverse transcriptase inhibitors (NNRTI) were analyzed according to the Stanford University HIV Drug Resistance Database. Viral sequences were screened for G-to-A hypermutations. A threshold of 2% was used for hypermutation adjustment. Viral suppression was considered at HIV-1 RNA load ≤400 copies/ml. RESULTS: Among 4973 participants with HIV-1C sequences, ART data were available for 4927 (99%) including 3858 (78%) on ART. Among those on ART, 3435 had viral load data and 3297 (96%) were virologically suppressed. Among 1069 (22%) HIV-infected individuals not on ART, we found NRTI-associated and NNRTI-associated SDRMs were found in 1.5% (95% confidence interval [CI] 1.0-2.5%) and 2.9% (95% CI 2.0-4.2%), respectively. Of the 138 (4%) of individuals who had detectable HIV-1 RNA, we found NRTI-associated and NNRTI-associated drug resistance mutations in 16% (95% CI 10-25%) and 33% (95% CI 25-42%), respectively. CONCLUSION: We found a low prevalence of NRTI-associated and NNRTI-associated PDR-resistance mutations among residents of rural and peri-urban communities across Botswana. However, individuals on ART with detectable virus had ADR NRTI and NNRTI mutations above 15%.