Phylodynamic Structure in the Botswana HIV Epidemic.

Background Studying viral sequences can provide insights into the structure of host contact networks through which the virus is transmitted. Uncovering the population structure of the HIV-1 epidemic in Botswana will help optimise public health interventions and may identify hidden sub-epidemics. We sought to determine the phylodynamic structure of the Botswana HIV-1 epidemic from viral sequence genetic data. Methods The Botswana Combination Prevention Project (BCPP) randomly sampled 20% of households in 30 villages in Botswana between 2013-2018 and tested for HIV-1. Extensive demographic data were collected from all participants and next-generation full-genome HIV-1 sequences were generated from HIV-1 positive participants (n = 4,164), 78% of whom were on antiretroviral treatment (ART). We inferred the stage of infection (< or > 1 year) among HIV-1 cases based on nucleotide diversity and clinical data using a previously trained machine learning model. We then reconstructed time-resolved gag and pol phylogenies from sequences, other Botswana cohorts and publicly available sequences that were genetically close to those from Botswana. We statistically explored phylogenies for partitions with diverging patterns of coalescence, indicating sub-epidemics, and estimated viral effective population size through time, a measure of viral incidence, for each partition. Finally, we compared the demographic makeup, clinical and geographic characteristics across partitions using χ2, ANOVA tests and Tukey analysis. Results We identified three partitions of time-resolved gag and pol phylogenies, revealing divergent patterns of coalescence and HIV-1 transmission. In both gag and pol phylogenies, partitions with persistent growth and transmission were characterised by lower treatment coverage and more recent infections when compared to other partitions. The Southern and South East regions of Botswana were over-represented in the fast-growing partitions. Conclusion Our findings suggest that transmission is slowing in segments of the population that have high ART coverage. However, recent infections are over-represented in ongoing sub-epidemics. The phylodynamic structure suggests that there are districts with higher growth and prioritising these in the deployment of public health interventions might curb new infections. Nonetheless the high mobility of Botswana residents should be taken into consideration in implementing effective interventions to combat HIV-1.

Hepatitis Delta Virus Clade 8 Is the Predominant Clade Circulating in Botswana amongst People Living with HIV.

Hepatitis delta virus (HDV) co-infections more often result in severe hepatitis compared to hepatitis B virus (HBV) infections alone. Despite a high HDV prevalence (7.1%), information regarding circulating HDV clades is very limited in Botswana. We extracted total nucleic acid from confirmed HDV-positive samples and quantified their viral load. We then sequenced the large hepatitis delta antigen (L-HDAg) using Oxford Nanopore Technology (ONT). Genotyping was performed using the HDV Database, and HDV mutation profiling was performed on AliView. All participants with HBV genotypic information belonged to sub-genotype A1, and 80% (4/5) of them had a higher HDV viral load and a lower HBV viral load. We sequenced 75% (9/12) of the HDV-positive samples, which belonged to HDV clade 8. A total of 54 mutations were discovered, with the most prevalent being Q148R (16%), D149P (16%) and G151D (16%). Known mutations such as S117A, K131R, R139K and G151D were detected, while the other mutations were novel. Our results reveal that HDV clade 8 is the predominant clade in Botswana. The significance of all mutations remains unclear. Future studies with a larger sample size to detect other HDV clades that might be circulating in Botswana and functionally characterize the detected mutations are warranted.

Impact of Hepatitis Delta Virus Infection on the Selection of Hepatitis B Surface Antigen Mutations.

The interaction of multiple viruses in one host is thought to enhance the development of mutations. However, the impact of hepatitis D virus (HDV) positivity on the development of unique hepatitis B virus (HBV) mutations among people living with human immunodeficiency virus (HIV) (PLWH) remains poorly understood in African countries, including Botswana. We used HBV sequences generated from the Botswana Combination Prevention Project (BCPP), which is the largest pair-matched cluster-randomized HIV trial in Botswana. Only participants with available HBV sequences (n = 55) were included in our study ([HIV/HBV-positive (n = 50) and HIV/HBV/HDV-positive (n = 5)]. Geno2pheno was used to determine HBV genotypes, and HBV surface region sequences (all subgenotype A1) were aligned in AliView for mutational analysis, while the impact of mutations was assessed using Phyre2. Our results identified 182 common mutations between the two groups. In the HIV/HBV/HDV cohort, only three mutations (L95W, W156Q, C221Y) were classified as deleterious, with only L95W being the most frequent. In the HIV/HBV cohort, four mutations (W199R, C221A, C221S, W223G) were also classified as deleterious. Our results demonstrate the presence of unique HBV mutations among the HIV/HBV/HDV-positive cohort. Functional characterization of these mutations is recommended to determine their effect on HDV.

Low Prevalence of Nirmatrelvir-Ritonavir Resistance-Associated Mutations in SARS-CoV-2 Lineages From Botswana.

Background: We evaluated naturally occurring nirmatrelvir-ritonavir (NTV/r) resistance-associated mutations (RAMs) among severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) strains from Botswana, a country with no NTV/r use to date, in order to recommend the usage of the agent for high-risk patients with coronavirus disease 2019 (COVID-19). Methods: We conducted a retrospective analysis using 5254 complete SARS-CoV-2 sequences from Botswana (September 2020-September 2023). We evaluated the mutational landscape of SARS-CoV-2 3-Chymotrypsin-like protease (3CLpro) relative to the highlighted list of RAMs granted Food and Drug Administration Emergency Use Authorization in 2023. Results: The sequenced 5254 samples included Beta variants of concerns (VOCs; n = 323), Delta VOCs (n = 1314), and Omicron VOCs (n = 3354). Overall, 77.8% of the sequences exhibited at least 1 polymorphism within 76/306 amino acid positions in the nsp5 gene. NTV/rRAMs were identified in 34/5254 (0.65%; 95% CI, 0.43%-0.87%) and occurred at 5 distinct positions. Among the NTV/r RAMS detected, A191V was the most prevalent (24/34; 70.6%). Notably, T21I mutation had a prevalence of 20.6% (7/34) and coexisted with either K90R (n = 3) polymorphism in Beta sequences with RAMs or P132H (n = 3) polymorphism for Omicron sequences with RAMs. Other NTV/r RAMs detected included P108S, with a prevalence of 5.88% (2/34), and L50F, with a prevalence of 2.94% (1/34). NTV/r RAMs were significantly higher (P < .001) in Delta (24/35) compared with Beta (4/34) and Omicron (6/34) sequences. Conclusions: The frequency of NTV/r RAMs in Botswana was low. Higher rates were observed in Delta VOCs compared to Omicron and Beta VOCs. As NTV/r use expands globally, continuous surveillance for drug-resistant variants is essential, given the RAMs identified in our study.

Molecular Characterization of Hepatitis B Virus in People Living with HIV in Rural and Peri-Urban Communities in Botswana.

(1) Background: Hepatitis B virus (HBV) sequencing data are important for monitoring HBV evolution. We aimed to molecularly characterize HBV sequences from participants with HBV surface antigen-positive (HBsAg+) serology and occult hepatitis B infection (OBI+). (2) Methods: We utilized archived plasma samples from people living with human immunodeficiency virus (PLWH) in Botswana. HBV DNA was sequenced, genotyped and analyzed for mutations. We compared mutations from study sequences to those from previously generated HBV sequences in Botswana. The impact of OBI-associated mutations on protein function was assessed using the Protein Variation Effect Analyzer. (3) Results: Sequencing success was higher in HBsAg+ than in OBI+ samples [86/128 (67.2%) vs. 21/71 (29.2%)]. Overall, 93.5% (100/107) of sequences were genotype A1, 2.8% (3/107) were D3 and 3.7% (4/107) were E. We identified 13 escape mutations in 18/90 (20%) sequences with HBsAg coverage, with K122R having the highest frequency. The mutational profile of current sequences differed from previous Botswana HBV sequences, suggesting possible mutational changes over time. Mutations deemed to have an impact on protein function were tpQ6H, surfaceV194A and preCW28L. (4) Conclusions: We characterized HBV sequences from PLWH in Botswana. Escape mutations were prevalent and were not associated with OBI. Longitudinal HBV studies are needed to investigate HBV natural evolution.

OPTIMIZATION OF DRIED BLOOD SPOT ASSAYS FOR HEPATITIS B VIRUS SURFACE ANTIBODY QUANTIFICATION.

Dried blood spot (DBS) cards can be used as an alternative sample collection method to plasma, however, there is no optimized elution protocol for DBS cards specifically for hepatitis B surface antibody (anti-HBs) testing. The study aimed to develop a DBS elution protocol for anti-HBs quantification. Our study sought to determine the ideal phosphate-buffered saline (PBS) buffer volume to use by comparing three PBS volumes (300uL, 450uL, and 500uL), and the optimal time to agitate DBS discs on a plate shaker (1hr, 2hrs, 3hrs, and 4hrs) to yield DBS anti-HBs concentrations that are comparable to corresponding plasma anti-HBs concentrations. The optimal DBS storage temperature (25°C, -20°C, and -80°C) was investigated to determine the ideal long-term storage temperature of the cards. Residual samples were used for optimization (2019-2021). A total of 50 DBS-plasma pairs was used throughout the study, with plasma anti-HBs concentrations being used as the golden standard to compare. The analysis of results was carried out by determining the p-values of the Wilcoxon sign rank. A two-way analysis of variance (ANOVA) was also performed to determine the impact of PBS elution volumes, elution time, and storage temperature on the anti-HBs concentration of DBS samples on STATA Version 15.0. No statistically significant difference between the DBS-plasma anti-HBs pairs was observed when using 450 or 500uL of PBS buffer and when samples were agitated for 3 hours (p=0.594, p=0.499 respectively). The optimal storage temperature for DBS cards was 25°C because the results showed no statistically significant difference between DBS-plasma anti-HBs titers (p=0.594). The two-way ANOVA analysis showed that elution volumes and time had no statistically significant impact on the DBS anti-HBs concentrations, p=0.948 and p=0.381 respectively. Storage temperature had a statistically significant impact on the DBS anti-HBs concentrations, p=0.002. The optimized DBS elution protocol for anti-HBs quantification will help monitor vaccine efficacy in infants due to the low sample volumes required compared to plasma and also can be used for anti-HBs testing in resource-limited areas around the country.

Hepatitis B Virus Prevalence among HIV-Uninfected People Living in Rural and Peri-Urban Areas in Botswana.

(1) Background: we determined the prevalence of the hepatitis B virus (HBV) amongst people without human immunodeficiency virus (HIV) in rural and peri-urban areas in Botswana. (2) Methods: We screened for the hepatitis B surface antigen (HBsAg) from archived plasma samples of people without HIV (n = 2135) randomly selected from the Botswana Combination Prevention Program (BCPP) (2013-2018). We sequenced 415 bp of the surface region using BigDye sequencing chemistry. (3) Results: The median age of participants was 31 (IQR: 24-46) and 64% (1360/2135) were female. HBV prevalence was 4.0% (86/2135) [95% CI: 3.3-4.9]) and ranged between 0-9.2%. Older participants (>35 years) had increased odds of HBV positivity (OR: 1.94; 95% CI: [1.32-2.86]; p = 0.001). Thirteen samples were sequenced and seven (53.8%) were genotype A, three (23.1%) were genotype D and genotype E each. Clinically significant mutations were identified in the surface region, but no classic drug resistance mutations were identified. (4) Conclusions: We report an HBV prevalence of 4.0% (95% CI 3.3-4.9) among people without HIV in rural and peri-urban communities in Botswana with varying rates in different communities. A comprehensive national HBV program is required in Botswana to guide HBV prevention, testing and management.

Persistence and risk factors of occult hepatitis B virus infections among antiretroviral therapy-naïve people living with HIV in Botswana.

Aim: This study aimed to determine the kinetics of occult hepatitis B virus infections (OBI) among people with HIV (PWH). Methods: The study used archived plasma samples from longitudinal HIV natural history studies. We identified new OBI cases and assessed risk factors for OBI using Cox proportional hazards regression analysis. Results: At baseline, 8 of 382 [(2.1%) (95% CI: 1.06-4.1)] samples tested positive for hepatitis B surface antigen (HBsAg+). Of the 374 HBsAg-negative samples, 76 had sufficient sample volume for HBV DNA screening. OBI positivity (OBI+) at baseline was reported in 11 of 76 [14.7 95% CI (8.3-24.1)] HBsAg-negative (HBsAg-) participants. Baseline HBsAg-negative samples with sufficient follow-up samples (n = 90) were used for analysis of newly identified OBI cases. Participants contributed 129.74 person-years to the study and were followed for a median of 1.02 years (IQR: 1.00-2.00). Cumulatively, there were 34 newly identified OBI cases from the 90 participants, at the rate of 26.2/100 person-years (95% CI: 18.7-36.7). Newly identified OBI cases were more common among men than women (61.1% vs. 31.9%) and among participants with CD4+ T-cell counts ≤450 cells/mL (p-value = 0.02). Most of the newly identified OBI cases [55.9% (19/34)] were possible reactivations as they were previously HBV core antibody positive. Conclusion: There was a high rate of newly identified OBI among young PWH in Botswana, especially in men and in participants with lower CD4+ T-cell counts. OBI screening in PWH should be considered because of the risk of transmission, possible reactivation, and risk factors for the development of chronic liver disease, including hepatocellular carcinoma.

Low-Level Viremia among Adults Living with HIV on Dolutegravir-Based First-Line Antiretroviral Therapy Is a Predictor of Virological Failure in Botswana.

We evaluated subsequent virologic outcomes in individuals experiencing low-level virem ia (LLV) on dolutegravir (DTG)-based first-line antiretroviral therapy (ART) in Botswana. We used a national dataset from 50,742 adults who initiated on DTG-based first-line ART from June 2016-December 2022. Individuals with at least two viral load (VL) measurements post three months on DTG-based first-line ART were evaluated for first and subsequent episodes of LLV (VL:51-999 copies/mL). LLV was sub-categorized as low-LLV (51-200 copies/mL), medium-LLV (201-400 copies/mL) and high-LLV (401-999 copies/mL). The study outcome was virologic failure (VF) (VL ≥ 1000 copies/mL): virologic non-suppression defined as single-VF and confirmed-VF defined as two-consecutive VF measurements after an initial VL < 1000 copies/mL. Cox regression analysis identified predictive factors of subsequent VF. The prevalence of LLV was only statistically different at timepoints >6-12 (2.8%) and >12-24 (3.9%) (p-value < 0.01). LLV was strongly associated with both virologic non-suppression (adjusted hazards ratio [aHR] = 2.6; 95% CI: 2.2-3.3, p-value ≤ 0.001) and confirmed VF (aHR = 2.5; 95% CI: 2.4-2.7, p-value ≤ 0.001) compared to initially virally suppressed PLWH. High-LLV (HR = 3.3; 95% CI: 2.9-3.6) and persistent-LLV (HR = 6.6; 95% CI: 4.9-8.9) were associated with an increased hazard for virologic non-suppression than low-LLV and a single-LLV episode, respectively. In a national cohort of PLWH on DTG-based first-line ART, LLV > 400 copies/mL and persistent-LLV had a stronger association with VF. Frequent VL testing and adherence support are warranted for individuals with VL > 50 copies/mL.

High Prevalence of Hepatitis B Virus Drug Resistance Mutations to Lamivudine among People with HIV/HBV Coinfection in Rural and Peri-Urban Communities in Botswana.

(1) Background: We aimed to determine the prevalence of hepatitis B virus (HBV) resistance-associated mutations (RAMs) in people with HBV and human immunodeficiency virus (HBV/HIV) in Botswana. (2) Methods: We sequenced HBV deoxyribonucleic acid (DNA) from participants with HBV/HIV from the Botswana Combination Prevention Project study (2013-2018) using the Oxford Nanopore GridION platform. Consensus sequences were analyzed for genotypic and mutational profiles. (3) Results: Overall, 98 HBV sequences had evaluable reverse transcriptase region coverage. The median participant age was 43 years (IQR: 37, 49) and 66/98 (67.4%) were female. Most participants, i.e., 86/98 (87.8%) had suppressed HIV viral load (VL). HBV RAMs were identified in 61/98 (62.2%) participants. Most RAMs were in positions 204 (60.3%), 180 (50.5%), and 173 (33.3%), mostly associated with lamivudine resistance. The triple mutations rtM204V/L180M/V173L were the most predominant (17/61 [27.9%]). Most participants (96.7%) with RAMs were on antiretroviral therapy for a median duration of 7.5 years (IQR: 4.8, 10.5). Approximately 27.9% (17/61) of participants with RAMs had undetectable HBV VL, 50.8% (31/61) had VL < 2000 IU/mL, and 13/61 (21.3%) had VL ≥ 2000 IU/mL. (4) Conclusions: The high prevalence of lamivudine RAMs discourages the use of ART regimens with 3TC as the only HBV-active drug in people with HIV/HBV.

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.

Hepatitis B surface antigen loss in individuals with chronic hepatitis B virus and HIV-1 infections in Botswana.

OBJECTIVES: We sought to determine hepatitis B surface antigen (HBsAg) loss and its predictors among people with chronic hepatitis B (CHB) infections and HIV (PWH) in Botswana. METHODS: Archived plasma samples from a cohort of PWH in Botswana (2013-2018) with 3 yearly time-points were used. Samples were screened for HBsAg, immunoglobulin M HBV core antibodies (anti-HBc IgM) and HBV e-antigen (HBeAg) at all time points. HBV deoxyribonucleic acid (DNA) quantification was done at baseline. The Wilcoxon rank-sum was used to compare continuous variables while the chi-squared test and Fishers exact test were used for categorical data wherever appropriate. Logistic regression was used to assess predictors of seroclearance. RESULTS: Of 141 participants with HBsAg-positive serology (HBsAg+) at baseline, 92.2% (131/141) [95% confidence interval (CI) 87.4-96.1] were persistently HBsAg+ at year 1. We report a HBsAg loss of 7.1% (10/141) (95% CI 3.9-12.6) among participants with negative HBeAg and negative IgM serologies. HBsAg loss was 6.3% (7/111) among antiretroviral therapy (ART)-experienced participants and 10.7% (3/28) (95% CI 0.4-5.0) in ART-naive participants. Most participants who had positive anti-HBc IgM serology and did not lose HBsAg were on either lamivudine (3TC)-based therapy or non-tenofovir disoproxil fumarate (TDF)-based therapy, except for one participant. The participants also had varying HBeAg status. HBsAg loss was independent of HIV viral load, CD4 + cell count, age, and sex. CONCLUSION: We report a HBsAg loss of 6.3% over a 3-year period among ART-experienced CHB participants. Future studies that focus on HBsAg loss in mono-infected patients and the possible correlation between HBeAg status and HBsAg loss are warranted.

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.

Correction: Youth and healthcare workers' perspectives on the feasibility and acceptability of self-testing for HIV, Hepatitis and Syphilis among young people: Qualitative findings from a pilot study in Gaborone, Botswana.

[This corrects the article DOI: 10.1371/journal.pone.0288971.].

Developing a diversity, equity and inclusion compass to guide scientific capacity strengthening efforts in Africa.

Diversity, equity and inclusion (DEI) in science is vital to improve the scientific process and ensure societal uptake and application of scientific results. DEI challenges include a full spectrum of issues from the lack of, and promotion of, women in science, to the numerous barriers in place that limit representation of African scientists in global scientific efforts. DEI principles in African science remain relatively underdeveloped, with limited engagement and discussion among all stakeholders to ensure that initiatives are relevant to local environments. The Sub-Saharan African Network for TB/HIV research Excellence (SANTHE) is a network of African-led research in HIV, tuberculosis (TB), associated co-morbidities, and emerging pathogens, now based in eight African countries. Our aim, as a scientific capacity strengthening network, was to collaboratively produce a set of DEI guidelines and to represent them visually as a DEI compass. We implemented a consortium-wide survey, focus group discussions and a workshop where we were able to identify the key DEI challenges as viewed by scientists and support staff within the SANTHE network. Three thematic areas were identified: 1. Conquering Biases, 2. Respecting the Needs of a Diverse Workforce (including mental health challenges, physical disability, career stability issues, demands of parenthood, and female-specific challenges), and 3. Promotion of African Science. From this we constructed a compass that included proposed steps to start addressing these issues. The use of the compass metaphor allows 're-adjustment/re-positioning' making this a dynamic output. The compass can become a tool to establish an institution's DEI priorities and then to progress towards them.

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.

High prevalence of hepatitis delta virus among people with hepatitis B virus and HIV coinfection in Botswana.

BACKGROUND: Approximately 15-20 million people worldwide are infected with hepatitis delta virus (HDV), which is approximately 5 % of people with chronic hepatitis B virus (HBV). Sub-Saharan Africa has high HDV prevalence, leading to worse clinical outcomes among people who are HIV/HBV/HDV tri-infected. There are limited data on HDV prevalence among people with HIV (PWH) who are HBV-infected and uninfected in Botswana. We, therefore, determined HDV prevalence among PWH in Botswana. METHODS: This was a retrospective cross-sectional study utilizing archived plasma samples from PWH with results for HBV markers such as hepatitis B surface antigen (HBsAg), hepatitis B core antibody (anti-HBc), immunoglobulin M antibody to hepatitis B core antigen (IgM anti-HBc) and hepatitis B e antigen (HBeAg). Samples were categorized according to their HBsAg status and screened for anti-HDV antibodies. Total nucleic acid was extracted from samples with a single positive anti-HDV result, and HDV ribonucleic acid (RNA) load was quantified using the Altona Diagnostic RealStar® HDV RT-PCR kit. Statistical analysis was performed using STATA version 14.0 where p-values < 0.05 were considered statistically significant. RESULTS: The study cohort (n = 478) included both HBsAg positive (44 %) and negative (56 %) participants, with a median age of 42 [IQR; 41-43]. Anti-HDV prevalence of (15/211) [7.1 %, 95 % CI: 4.4 - 11.4] was recorded among HBsAg positive participants, all of whom were IgM anti-HBc negative, while 5/6 participants were HBeAg negative. HDV RNA load was detected in 11/12 (92 %) anti-HDV-positive participants. No HDV prevalence was recorded among participants who were HBsAg negative, therefore, the overall HDV prevalence was (15/478) [3.1 %, 95 % CI: 1.9 - 5.1]. HIV viral load suppression was statistically insignificant, irrespective of HDV status. CONCLUSIONS: We report high HDV prevalence among HBsAg-positive PWH in Botswana. Most HDV-positive participants had active HDV infection, therefore, we recommend HDV screening in this cohort to guide their clinical care.

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.