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.

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.

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.

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.

Safety of AZD1222 COVID-19 vaccine and low Incidence of SARS-CoV-2 infection in Botswana following ChAdOx1(AZD1222) vaccination: A single-arm open-label interventional study - final study results.

Objectives: We report the final analysis of the single-arm open-label study evaluating the safety and COVID-19 incidence after AZD1222 vaccination in Botswana conducted between September 2021 and August 2022. Methods: The study included three groups of adults (>18 years), homologous AZD1222 primary series and booster (AZ2), heterologous primary series with one dose AZD1222, and AZD1222 booster (HPS), and primary series other than AZD1222 and AZD1222 booster (OPS). We compared the incidence of AEs in participants with and without prior COVID-19 infection using an exact test for rate ratios. Results: Among 10,894 participants, 9192 (84.4%) were enrolled at first vaccine dose, 521 (4.8%) at second vaccine, and 1181 (10.8%) at the booster vaccine. Of 10,855 included in the full analysis set, 1700 received one dose of AZD1222; 5377 received two doses; 98 received a heterologous series including one AZD1222 and a booster; 30 in the HPS group; 1058 in the OPS group; and 2592 in the AZ2 group. No laboratory-confirmed COVID-19 hospitalizations or deaths were reported. The incidence of laboratory-confirmed symptomatic COVID infection for the AZ2 group was 6.22 (95% confidence interval: 2.51-12.78) per 1000 participant-years (1000-PY) and 3.5 (95% confidence interval: 0.42-12.57) per 1000-PY for AZ2+booster group. Most adverse events were mild, with higher incidence in participants with prior COVID-19 infection. Individuals with prior COVID-19 exposure exhibited higher binding antibody responses. No differences in outcomes were observed by HIV status. Conclusion: AZD1222 is safe, effective, and immunogenic for people living with and without HIV.

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.

Atypical Hepatitis B Virus Serology Profile-Hepatitis B Surface Antigen-Positive/Hepatitis B Core Antibody-Negative-In Hepatitis B Virus/HIV Coinfected Individuals in Botswana.

(1) Background: Hepatitis B core antibodies (anti-HBc) are a marker of hepatitis B virus (HBV) exposure; hence, a normal HBV serology profile is characterized by HBV surface antigen (HBsAg) and anti-HBc positivity. However, atypical HBV serologies occur, and we aimed to determine the prevalence of an atypical profile (HBsAg+/anti-HBc-) in a cohort of people with HIV-1 (PWH) in Botswana. (2) Methods: Plasma samples from an HIV-1 cohort in Botswana (2013-2018) were used. The samples were screened for HBsAg and anti-HBc. Next-generation sequencing was performed using the GridION platform. The Wilcoxon rank-sum test and Chi-squared tests were used for the comparison of continuous and categorical variables, respectively. (3) Results: HBsAg+/anti-HBc- prevalence was 13.7% (95% CI 10.1-18.4) (36/263). HBsAg+/anti-HBc- participants were significantly younger (p < 0.001), female (p = 0.02) and ART-naïve (p = 0.04) and had a detectable HIV viral load (p = 0.02). There was no statistically significant difference in the number of mutations observed in participants with HBsAg+/anti-HBc- vs. those with HBsAg+/anti-HBc+ serology. (4) Conclusions: We report a high HBsAg+/anti-HBc- atypical serology profile prevalence among PWH in Botswana. We caution against HBV-testing algorithms that consider only anti-HBc+ samples for HBsAg testing, as they are likely to underestimate HBV prevalence. Studies to elucidate the mechanisms and implications of this profile are warranted.

Near-complete genome of SARS-CoV-2 Delta variant of concern identified in a symptomatic dog (Canis lupus familiaris) in Botswana.

We sought to investigate whether SARS-CoV-2 was present, and to perform full-length genomic sequencing, in a 5-year-old male crossbreed dog from Gaborone, Botswana that presented overt clinical signs (flu-like symptoms, dry hacking cough and mild dyspnoea). It was only sampled a posteriori, because three adult owners were diagnosed with SARS-CoV-2 infection. Next-generation sequencing based on Oxford Nanopore Technology (ONT) was performed on amplicons that were generated using a reverse transcriptase real-time polymerase chain reaction (RT-qPCR) of confirmed positive SARS-CoV-2 nasopharyngeal and buccal swabs, as well as a bronchoalveolar lavage with mean real cycle threshold (qCt) value of 36 based on the Nucleocapsid (N) gene. Descriptive comparisons to known sequences in Botswana and internationally were made using mutation profiling analysis and phylogenetic inferences. Human samples were not available. A near-full length SARS-CoV-2 genome (∼90% coverage) was successfully genotyped and classified under clade 20 O and Pango-Lineage AY.43 (Pango v.4.0.6 PLEARN-v1.3; 2022-04-21), which is a sublineage of the Delta variant of concern (VOC) (formerly called B.1.617.2, first detected in India). We did not identify novel mutations that may be used to distinguish SARS-CoV-2 isolates from the dog and humans. In addition to Spike (S) region mutation profiling, we performed phylogenetic analysis including 30 Delta sequences publicly available reference also isolated from dogs. In addition, we performed another exploratory analysis to investigate the phylogenetic relatedness of sequence isolated from dog with those from humans in Botswana (n = 1303) as of 31 March 2022 and of same sublineage. Expectedly, the sequence formed a cluster with Delta sublineages - AY.43, AY.116 and B.1.617.2 - circulating in same time frame. This is the first documented report of human-associated SARS-CoV-2 infection in a dog in Botswana. Although the direction of transmission remains unknown, this study further affirms the need for monitoring pets during different COVID-19 waves for possible clinically relevant SARS-CoV-2 transmissions between species.

High Prevalence of Hepatitis B Virus Infection Among People With HIV in Rural and Periurban Communities in Botswana.

Background: We aimed to determine the prevalence of hepatitis B virus (HBV) infection among people with human immunodeficiency virus (PWH) in rural and periurban communities in Botswana. Methods: PWH from a previous population-based study, the Botswana Prevention Combination Project, which enrolled adults in 30 communities across Botswana (2013-2018), were screened for HBV surface antigen (HBsAg) and HBV core antibody (anti-HBc). HBsAg-positive (HBsAg+) samples were further screened for HBV core immunoglobulin M antibodies (anti-HBc immunoglobulin M [IgM]) and HBV e antigen (HBeAg). We quantified HBV viral load on participants who tested positive (n = 148) and negative for HBsAg (n = 381). Results: Of 3304 participants tested, 271 (8% [95% confidence interval {CI}, 7%-9%]) were HBsAg+ while 1788 (56% [95% CI, 54%-57%]) of 3218 PWH whom we tested had positive anti-HBc. Approximately 88% of HBsAg+ participants were on antiretroviral therapy (ART), 40% and 56% of whom were receiving lamivudine- and tenofovir-containing ART, respectively. Male sex (relative risk ratio [RRR], 1.8 [95% CI, 1.2-2.7]) and the northern geographic region (RRR, 2.5 [95% CI, 1.4-4.7]) were independent predictors of HBV infection (HBsAg+). Of 381 persons with negative HBsAg who were tested for occult HBV, 126 (33% [95% CI, 29%-38%]) had positive HBV DNA. Eleven participants were highly viremic with high HBV viral load while on a lamivudine- or tenofovir-containing regimen. Ten (91%) of these participants also had positive HBeAg serology, while 4 (36%) had positive anti-HBc IgM serology. Conclusions: The prevalence of HBV was high among PWH in Botswana while on ART regimens with activity against HBV.

Whole genome sequencing reveals population diversity and variation in HIV-1 specific host genes.

HIV infection continues to be a major global public health issue. The population heterogeneity in susceptibility or resistance to HIV-1 and progression upon infection is attributable to, among other factors, host genetic variation. Therefore, identifying population-specific variation and genetic modifiers of HIV infectivity can catapult the invention of effective strategies against HIV-1 in African populations. Here, we investigated whole genome sequences of 390 unrelated HIV-positive and -negative individuals from Botswana. We report 27.7 million single nucleotide variations (SNVs) in the complete genomes of Botswana nationals, of which 2.8 million were missing in public databases. Our population structure analysis revealed a largely homogenous structure in the Botswana population. Admixture analysis showed elevated components shared between the Botswana population and the Niger-Congo (65.9%), Khoe-San (32.9%), and Europeans (1.1%) ancestries in the population of Botswana. Statistical significance of the mutational burden of deleterious and loss-of-function variants per gene against a null model was estimated. The most deleterious variants were enriched in five genes: ACTRT2 (the Actin Related Protein T2), HOXD12 (homeobox D12), ABCB5 (ATP binding cassette subfamily B member 5), ATP8B4 (ATPase phospholipid transporting 8B4) and ABCC12 (ATP Binding Cassette Subfamily C Member 12). These genes are enriched in the glycolysis and gluconeogenesis (p < 2.84e-6) pathways and therefore, may contribute to the emerging field of immunometabolism in which therapy against HIV-1 infection is being evaluated. Published transcriptomic evidence supports the role of the glycolysis/gluconeogenesis pathways in the regulation of susceptibility to HIV, and that cumulative effects of genetic modifiers in glycolysis/gluconeogenesis pathways may potentially have effects on the expression and clinical variability of HIV-1. Identified genes and pathways provide novel avenues for other interventions, with the potential for informing the design of new therapeutics.