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 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.

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.

No Difference in the Prevalence of HIV-1 gag Cytotoxic T-Lymphocyte-Associated Escape Mutations in Viral Sequences from Early and Late Parts of the HIV-1 Subtype C Pandemic in Botswana.

HIV is known to accumulate escape mutations in the gag gene in response to the immune response from cytotoxic T lymphocytes (CTLs). These mutations can occur within an individual as well as at a population level. The population of Botswana exhibits a high prevalence of HLA*B57 and HLA*B58, which are associated with effective immune control of HIV. In this retrospective cross-sectional investigation, HIV-1 gag gene sequences were analyzed from recently infected participants across two time periods which were 10 years apart: the early time point (ETP) and late time point (LTP). The prevalence of CTL escape mutations was relatively similar between the two time points-ETP (10.6%) and LTP (9.7%). The P17 protein had the most mutations (9.4%) out of the 36 mutations that were identified. Three mutations (A83T, K18R, Y79H) in P17 and T190A in P24 were unique to the ETP sequences at a prevalence of 2.4%, 4.9%, 7.3%, and 5%, respectively. Mutations unique to the LTP sequences were all in the P24 protein, including T190V (3%), E177D (6%), R264K (3%), G248D (1%), and M228L (11%). Mutation K331R was statistically higher in the ETP (10%) compared to the LTP (1%) sequences (p < 0.01), while H219Q was higher in the LTP (21%) compared to the ETP (5%) (p < 0.01). Phylogenetically, the gag sequences clustered dependently on the time points. We observed a slower adaptation of HIV-1C to CTL immune pressure at a population level in Botswana. These insights into the genetic diversity and sequence clustering of HIV-1C can aid in the design of future vaccine strategies.

High concordance in plasma and CSF HIV-1 drug resistance mutations despite high cases of CSF viral escape in individuals with HIV-associated cryptococcal meningitis in Botswana.

OBJECTIVES: We compared the patterns of HIV-1 drug resistance mutations between the CSF and plasma of individuals with HIV-associated cryptococcal meningitis. METHODS: This is a cross-sectional study of archived CSF and plasma samples collected from ART-exposed participants recruited in the Phase 3 AmBisome Therapy Induction Optimisation randomized controlled trial (ISRCTN72509687) conducted in Botswana between 2018 and 2021. HIV-1 RT and protease genes were genotyped using next-generation sequencing and HIV-1 drug resistance mutations were compared between the CSF and plasma compartments stratified by thresholds of ≥20% and <20%. RESULTS: Overall, 66.7% (16/24) of participants had at least one HIV-1 drug resistance mutation in the CSF and/or plasma. A total of 15/22 (68.2%) participants had HIV-1 drug resistance mutations at ≥20% threshold in the plasma and of those, 11 (73.3%) had been on ART longer than 6 months. HIV-1 drug resistance mutations were highly concordant between the CSF and plasma at ≥20% threshold despite a substantial number of individuals experiencing CSF viral escape and with only 54.5% with CSF WBC count ≥20 cells/mm3. Minority HIV-1 drug resistance mutations were detected in 20.8% (5/24) of participants. There were no mutations in the CSF that were not detected in the plasma. CONCLUSIONS: There was high concordance in HIV-1 drug resistance mutations in the CSF and plasma, suggesting intercompartmental mixing and possibly a lack of compartmentalization. Some individuals harboured minority HIV-1 drug resistance mutations, demonstrating the need to employ more sensitive genotyping methods such as next-generation sequencing for the detection of low-abundance mutations.

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-1C env and gag Variation in the Cerebrospinal Fluid and Plasma of Patients with HIV-Associated Cryptococcal Meningitis in Botswana.

HIV-1 compartmentalization in reservoir sites remains a barrier to complete HIV eradication. It is unclear whether there is variation in HIV-1 env and gag between cerebrospinal fluid (CSF) and plasma of individuals with HIV-associated cryptococcal meningitis (CM). We compared HIV-1 env characteristics and the gag cytotoxic T-lymphocyte (CTL) escape mutations from CSF and plasma samples. Employing population-based Sanger sequencing, we sequenced HIV-1 env from CSF of 25 patients and plasma of 26 patients. For gag, 15 CSF and 21 plasma samples were successfully sequenced. Of these, 18 and 9 were paired env and gag CSF/plasma samples, respectively. There was no statistically significant difference in the proportion of CCR5-using strains in the CSF and plasma, (p = 0.50). Discordant CSF/plasma virus co-receptor use was found in 2/18 pairs (11.1%). The polymorphisms in the HIV-1 V3 loop were concordant between the two compartments. From the HIV-1 gag sequences, three pairs had discordant CTL escape mutations in three different epitopes of the nine analyzed. These findings suggest little variation in the HIV-1 env between plasma and CSF and that the CCR5-using strains predominate in both compartments. HIV-1 gag CTL escape mutations also displayed little variation in CSF and plasma suggesting similar CTL selective pressure.

Differences in human immunodeficiency virus-1C viral load and drug resistance mutation between plasma and cerebrospinal fluid in patients with human immunodeficiency virus-associated cryptococcal meningitis in Botswana.

To determine effects of cryptococcal meningitis (CM) on human immunodeficiency virus (HIV)-1C cerebrospinal fluid (CSF) viral escape, CSF/plasma viral discordance, and drug resistance mutation (DRM) discordance between CSF and plasma compartments, we compared CSF and plasma viral load (VL) and DRMs in individuals with HIV-associated CM in Botswana.This cross-sectional study utilized 45 paired CSF/plasma samples from participants in a CM treatment trial (2014-2016). HIV-1 VL was determined and HIV-1 protease and reverse transcriptase genotyping performed. DRMs were determined using the Stanford HIV database. CSF viral escape was defined as HIV-1 ribonucleic acid ≥0.5 log10 higher in CSF than plasma and VL discordance as CSF VL > plasma VL.HIV-1 VL was successfully measured in 39/45 pairs, with insufficient sample volume in 6; 34/39 (87.2%) participants had detectable HIV-1 in plasma and CSF, median 5.1 (interquartile range: 4.7-5.7) and 4.6 (interquartile range:3.7-4.9) log10 copies/mL, respectively (P≤.001). CSF viral escape was present in 1/34 (2.9%) and VL discordance in 6/34 (17.6%). Discordance was not associated with CD4 count, antiretroviral status, fungal burden, CSF lymphocyte percentage nor mental status. Twenty-six of 45 (57.8%) CSF/plasma pairs were successfully sequenced. HIV-1 DRM discordance was found in 3/26 (11.5%); 1 had I84IT and another had M46MI in CSF only. The third had K101E in plasma and V106 M in CSF.Our findings suggest that HIV-1 escape and DRM discordance may occur at lower rates in participants with advanced HIV-disease and CM compared to those with HIV associated neurocognitive impairment.

Genetic diversity of Mycobacterium tuberculosis strains circulating in Botswana.

BACKGROUND: Molecular typing of Mycobacterium tuberculosis (M.tb) isolates can inform Tuberculosis (TB) control programs on the relative proportion of transmission driving the TB epidemic. There is limited data on the M. tb genotypes that are circulating in Botswana. The aim of this study was to generate baseline data on the genetic diversity of M.tb isolates circulating in the country. METHODS: A total of 461 M.tb isolates received at the Botswana National Tuberculosis Reference Laboratory between March 2012 and October 2013 were included in this study. Drug susceptibility testing was conducted using the BD BACTEC MGIT 960 System. M.tb strains were genotyped using spoligotyping and spoligotype patterns were compared with existing patterns in the SITVIT Web database. A subset of drug resistant isolates which formed spoligo clusters (n = 65) was additionally genotyped with 12-loci MIRU. Factors associated with drug resistance and clustering were evaluated using logistic regression. RESULTS: Of the 461 isolates genotyped, 458 showed 108 distinct spoligotype patterns. The predominant M.tb lineages were Lineage 4 (81.9%), Lineage 2 (9%) and Lineage 1 (7.2%). The predominant spoligotype families within Lineage 4 were LAM (33%), S (14%), T (16%), X (16%). Three hundred and ninety-two (86%) isolates could be grouped into 44 clusters (2-46 isolates per cluster); giving a clustering rate of 76%. We identified 173 (37.8%) drug resistant isolates, 48 (10.5%) of these were multi-drug resistant. MIRU typing of the drug resistant isolates allowed grouping of 46 isolates into 14 clusters, giving a clustering rate of 49.2%. There was no association between age, sex, treatment category, region and clustering. CONCLUSIONS: This study highlights the complexity of the TB epidemic in Botswana with multiple strains contributing to disease and provides baseline data on the population structure of M.tb strains in Botswana.