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.

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.

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.

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.

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.

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.

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.

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.

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.

Decreased hepatitis B virus vaccine response among HIV-positive infants compared with HIV-negative infants in Botswana.

OBJECTIVES: We sought to determine vaccine antibody titres and the prevalence of hepatitis B surface antigen (HBsAg) in both HIV-positive and HIV-negative infants born to HIV-positive mothers in Botswana. DESIGN: This was a retrospective cross-sectional study using 449 archived dried blood spot samples from both HIV-positive and HIV-negative infants collected between 2016 and 2018. METHODS: We screened dried blood spot samples for HBsAg and determined hepatitis B surface antibody titres. We determined hepatitis B virus (HBV) genotypes by amplifying 415 base-pairs of the surface region. RESULTS: HIV-positive infants mounted a significantly lower immune response to the HBV vaccine (P  < 0.001). Furthermore, a lower proportion of HIV-positive infants had protective hepatitis B surface antibody titres (74.5%) than HIV-negative infants (89.2%) (P < 0.001). HIV-positive infants were older and 50.9% of them had completed vaccination (P = 0.018). Of the 449 infant samples tested, three (0.67%) were positive for HBsAg. Of the three HBsAg-positive infants, two had protective titres (>10 mIU/ml). Two of the three HBV-positive infants were infected with genotype D3 and had no drug-resistance or escape mutations. CONCLUSION: Vaccine response was lower among HIV-positive infants compared with HIV-negative infants. HBV infections were observed in both HIV-positive and HIV-negative infants in Botswana. Studies to investigate additional preventive strategies to reduce HBV mother-to-child transmission are recommended.

Rapid epidemic expansion of the SARS-CoV-2 Omicron variant in southern Africa.

The SARS-CoV-2 epidemic in southern Africa has been characterized by three distinct waves. The first was associated with a mix of SARS-CoV-2 lineages, while the second and third waves were driven by the Beta (B.1.351) and Delta (B.1.617.2) variants, respectively1-3. In November 2021, genomic surveillance teams in South Africa and Botswana detected a new SARS-CoV-2 variant associated with a rapid resurgence of infections in Gauteng province, South Africa. Within three days of the first genome being uploaded, it was designated a variant of concern (Omicron, B.1.1.529) by the World Health Organization and, within three weeks, had been identified in 87 countries. The Omicron variant is exceptional for carrying over 30 mutations in the spike glycoprotein, which are predicted to influence antibody neutralization and spike function4. Here we describe the genomic profile and early transmission dynamics of Omicron, highlighting the rapid spread in regions with high levels of population immunity.

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.