Prediction of Coreceptor Tropism in HIV-1 Subtype C in Botswana.

It remains unknown whether the C-C motif chemokine receptor type 5 (CCR5) coreceptor is still the predominant coreceptor used by Human Immunodeficiency Virus-1 (HIV-1) in Botswana, where the HIV-1 subtype C predominates. We sought to determine HIV-1C tropism in Botswana using genotypic tools, taking into account the effect of antiretroviral treatment (ART) and virologic suppression. HIV-1 gp120 V3 loop sequences from 5602 participants were analyzed for viral tropism using three coreceptor use predicting algorithms/tools: Geno2pheno, HIV-1C Web Position-Specific Score Matrices (WebPSSM) and the 11/25 charge rule. We then compared the demographic and clinical characteristics of people living with HIV (PLWH) harboring R5- versus X4-tropic viruses using χ2 and Wilcoxon rank sum tests for categorical and continuous data analysis, respectively. The three tools congruently predicted 64% of viruses as either R5-tropic or X4-tropic. Geno2pheno and the 11/25 charge rule had the highest concordance at 89%. We observed a significant difference in ART status between participants harboring X4- versus R5-tropic viruses. X4-tropic viruses were more frequent among PLWH receiving ART (χ2 test, p = 0.03). CCR5 is the predominant coreceptor used by HIV-1C strains circulating in Botswana, underlining the strong potential for CCR5 inhibitor use, even in PLWH with drug resistance. We suggest that the tools for coreceptor prediction should be used in combination.

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.

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.

Epidemiological and viral characteristics of undiagnosed HIV infections in Botswana.

BACKGROUND: HIV-1 is endemic in Botswana. The country's primary challenge is identifying people living with HIV who are unaware of their status. We evaluated factors associated with undiagnosed HIV infection using HIV-1 phylogenetic, behavioural, and demographic data. METHODS: As part of the Botswana Combination Prevention Project, 20% of households in 30 villages were tested for HIV and followed from 2013 to 2018. A total of 12,610 participants were enrolled, 3596 tested HIV-positive at enrolment, and 147 participants acquired HIV during the trial. Extensive socio-demographic and behavioural data were collected from participants and next-generation sequences were generated for HIV-positive cases. We compared three groups of participants: (1) those previously known to be HIV-positive at enrolment (n = 2995); (2) those newly diagnosed at enrolment (n = 601) and (3) those who tested HIV-negative at enrolment but tested HIV-positive during follow-up (n = 147). We searched for differences in demographic and behavioural factors between known and newly diagnosed group using logistic regression. We also compared the topology of each group in HIV-1 phylogenies and used a genetic diversity-based algorithm to classify infections as recent (< 1 year) or chronic (≥ 1 year). RESULTS: Being male (aOR = 2.23) and younger than 35 years old (aOR = 8.08) was associated with undiagnosed HIV infection (p < 0.001), as was inconsistent condom use (aOR = 1.76). Women were more likely to have undiagnosed infections if they were married, educated, and tested frequently. For men, being divorced increased their risk. The genetic diversity-based algorithm classified most incident infections as recent (75.0%), but almost none of known infections (2.0%). The estimated proportion of recent infections among new diagnoses was 37.0% (p < 0.001). CONCLUSION: Our results indicate that those with undiagnosed infections are likely to be young men and women who do not use condoms consistently. Among women, several factors were predictive: being married, educated, and testing frequently increased risk. Men at risk were more difficult to delineate. A sizeable proportion of undiagnosed infections were recent based on a genetic diversity-based classifier. In the era of "test and treat all", pre-exposure prophylaxis may be prioritized towards individuals who self-identify or who can be identified using these predictors in order to halt onward transmission in time.

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.

The role of CYP2B6 516G>T polymorphism on efavirenz/nevirapine toxicity. Implications on treatment outcomes: Lessons from Botswana.

ABSTRACT: The two non-nucleoside reverse transcriptase inhibitors (NNRTIs), efavirenz (EFV) and nevirapine (NVP), are currently the core antiretroviral drugs for treatment of HIV in sub-Saharan Africa including Botswana. The drugs are metabolized by Cytochrome P450 2B6 (CYP2B6) liver enzyme. The CYP2B6 gene that encodes for metabolism of these drugs is known to be highly polymorphic. One of the polymorphism in the CYP2B6 gene, 516G>T, particularly the 516T allele, is known to confer poor metabolism of EFV and NVP. This may lead to high levels of plasma drug concentrations and development of treatment toxicities, like central nervous system toxicities, and cutaneous and hepatic toxicities, for EFV and NVP, respectively. The CYP2B6 516G allele on the other hand is associated with an extensive metabolism of the two NNRTIs drugs. We sought to establish association between possible developments of NNRTIs toxicities with CYP2B6 516G>T variation in Botswana.A total of 316 peripheral blood mononuclear cells samples were used in a retrospective view. All the samples were from participants on EFV/NVP-containing regimen with known toxicity output. TaqMan Real-Time PCR approach was applied for assessing CYP2B6 516 allele variation in cases with treatment toxicity and those without. Analysis was performed by chi-square statistics and logistic regression analysis.The rate of poor metabolizers among participants with toxicity and those without toxicity was 18.4% and 15.1%, respectively. The CYP2B6 516 genotype distribution comparisons between the participants with toxicity and those without were not statistically different (chi-square = .326; P = .568).CYP2B6 516 variation was not associated with NNRTI toxicity. No other factors were associated with toxicity when considering age, baseline body mass index, baseline CD4, baseline HIV viral load and adherence. The results were discussed in the context of all the studies done in Botswana to date.

Health Impact and Cost-Effectiveness of HIV Testing, Linkage, and Early Antiretroviral Treatment in the Botswana Combination Prevention Project.

BACKGROUND: The Botswana Combination Prevention Project tested the impact of combination prevention (CP) on HIV incidence in a community-randomized trial. Each trial arm had ∼55,000 people, 26% HIV prevalence, and 72% baseline ART coverage. Results showed intensive testing and linkage campaigns, expanded antiretroviral treatment (ART), and voluntary male medical circumcision referrals increased coverage and decreased incidence over ∼29 months of follow-up. We projected lifetime clinical impact and cost-effectiveness of CP in this population. SETTING: Rural and periurban communities in Botswana. METHODS: We used the Cost-Effectiveness of Preventing AIDS Complications model to estimate lifetime health impact and cost of (1) earlier ART initiation and (2) averting an HIV infection, which we applied to incremental ART initiations and averted infections calculated from trial data. We determined the incremental cost-effectiveness ratio [US$/quality-adjusted life-years (QALY)] for CP vs. standard of care. RESULTS: In CP, 1418 additional people with HIV initiated ART and an additional 304 infections were averted. For each additional person started on ART, life expectancy increased 0.90 QALYs and care costs increased by $869. For each infection averted, life expectancy increased 2.43 QALYs with $9200 in care costs saved. With CP, an additional $1.7 million were spent on prevention and $1.2 million on earlier treatment. These costs were mostly offset by decreased care costs from averted infections, resulting in an incremental cost-effectiveness ratio of $79 per QALY. CONCLUSIONS: Enhanced HIV testing, linkage, and early ART initiation improve life expectancy, reduce transmission, and can be cost-effective or cost-saving in settings like Botswana.

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.

Human Immunodeficiency Virus (HIV) Genetic Diversity Informs Stage of HIV-1 Infection Among Patients Receiving Antiretroviral Therapy in Botswana.

BACKGROUND: Human immunodeficiency virus (HIV)-1 genetic diversity increases during infection and can help infer the time elapsed since infection. However, the effect of antiretroviral treatment (ART) on the inference remains unknown. METHODS: Participants with estimated duration of HIV-1 infection based on repeated testing were sourced from cohorts in Botswana (n = 1944). Full-length HIV genome sequencing was performed from proviral deoxyribonucleic acid. We optimized a machine learning model to classify infections as < or >1 year based on viral genetic diversity, demographic, and clinical data. RESULTS: The best predictive model included variables for genetic diversity of HIV-1 gag, pol, and env, viral load, age, sex, and ART status. Most participants were on ART. Balanced accuracy was 90.6% (95% confidence interval, 86.7%-94.1%). We tested the algorithm among newly diagnosed participants with or without documented negative HIV tests. Among those without records, those who self-reported a negative HIV test within <1 year were more frequently classified as recent than those who reported a test >1 year previously. There was no difference in classification between those self-reporting a negative HIV test <1 year, whether or not they had a record. CONCLUSIONS: These results indicate that recency of HIV-1 infection can be inferred from viral sequence diversity even among patients on suppressive ART.

Association of CYP2B6 Genetic Variation with Efavirenz and Nevirapine Drug Resistance in HIV-1 Patients from Botswana.

PURPOSE: CYP2B6 liver enzyme metabolizes the two non-nucleoside reverse transcriptase inhibitors Efavirenz (EFV) and Nevirapine (NVP) used in the antiretroviral therapy (ART) regimens for HIV-infected individuals. Polymorphisms of the CYP2B6 gene influence drug levels in plasma and possibly virological outcomes. The aim of this study was to explore the potential impact of CYP2B6 genotype and haplotype variation on the risk of developing EFV/NVP drug resistance mutations (DRMs) in HIV-1 patients receiving EFV-/NVP-containing regimens in Botswana. PATIENTS AND METHODS: Participants were a sub-sample of a larger study (Tshepo study) conducted in Gaborone, Botswana, among HIV-infected individuals taking EFV/NVP containing ART. Study samples were retrieved and assigned to cases (with DRMs) and controls (without DRMs). Four single-nucleotide polymorphisms (SNPs) in the CYP2B6 gene (-82T>C; 516G>T; 785A>G; 983T>C) were genotyped, the haplotypes reconstructed, and the metabolic score assigned. The possible association between drug resistance and several independent factors (baseline characteristics and CYP2B6 genotypes) was assessed by Binary Logistic Regression (BLR) analysis. EFV/NVP resistance status and CYP2B6 haplotypes were also analyzed using Z-test, chi-square and Fisher's exact test statistics. RESULTS: Two hundred and twenty-seven samples were analysed (40 with DRMs, 187 without DRMs). BLR analysis showed an association between EFV/NVP resistance and CYP2B6 516G allele (OR: 2.26; 95% CI: 1.27-4.01; P=0.005). Moreover, haplotype analysis revealed that the proportion of EFV/NVP-resistant infections was higher among CYP2B6 fast than extensive/slow metabolizers (30.8% vs 16.8%; P=0.035), with the 516G allele more represented in the haplotypes of fast than extensive/slow metabolizers (100.0% vs 53.8%; P<0.001). CONCLUSION: We demonstrated that the CYP2B6 516G allele, and even more when combined in fast metabolic haplotypes, is associated with the presence of EFV/NVP resistance, strengthening the need to assess the CYP2B6 genetic profiles in HIV-infected patients in order to improve the virologic outcomes of NNRTI containing ART.

Incidence of hepatitis B virus infection among human immunodeficiency virus-infected treatment naïve adults in Botswana.

Hepatitis B virus (HBV) and human immunodeficiency virus (HIV) coinfection is highest in sub-Saharan Africa and results in accelerated clinical outcomes compared with HBV or HIV mono-infection. HBV clearance rates are higher in healthy adults; however, in sub-Saharan Africa, there are limited data on clearance of incident HBV in HIV-infected adults. Therefore, we sought to estimate HBV incidence and HBV surface antigen (HBsAg) clearance in HIV-infected adults in Botswana.This was a retrospective longitudinal study of 442 HIV-1C infected treatment naïve patients enrolled in a previous Botswana Harvard AIDS Institute Partnership study. Archived plasma samples from 435 HIV-infected treatment naïve participants were screened for HBsAg and HBV core antibody (anti-HBc). HBsAg was evaluated annually over a 4-year period, and HBV deoxyribonucleic acid (DNA) levels of HBsAg-positive chronic and incident patients were quantified.Baseline median CD4+ T-cell count was 458 cells/µL [Q1, Q3: 373, 593], and median HIV viral load was 4.15 copies/mL [Q1, Q3: 3.46, 4.64]. Twenty two HBV incident cases occurred, representing an incidence of 3.6/100 person-years [95% CI: 2.2-5.6]. All incident HBV cases with a follow-up sample available for screening (13/22) cleared HBsAg. Detectable HBV viral loads among chronic and incident cases ranged between 5.15 × 10 to 1.4 × 10 IU/L and 1.80 × 10 to 1.7 × 10 IU/mL, respectively.We report high HBV incidence associated with elevated HBV DNA levels despite high CD4+ T-cell counts in HIV-infected patients in Botswana. These incidence cases represent a potential source of HBV transmission in the population. Scaling-up of HIV treatment strategies utilizing antiretroviral therapy regimens with anti-HBV activity coupled with screening for HBV infections in households of the HBsAg-positive cases is recommended.

Use of a mutation-specific genotyping method to assess for HIV-1 drug resistance in antiretroviral-naïve HIV-1 Subtype C-infected patients in Botswana.

Background: HIV-1 drug resistance poses a major threat to the success of antiretroviral therapy. The high costs of available HIV drug resistance assays prohibit their routine usage in resource-limited settings. Pan-degenerate amplification and adaptation (PANDAA), a focused genotyping approach based on quantitative PCR (qPCR), promises a fast and cost-effective way to detect HIV drug resistance mutations (HIVDRMs).  Given the high cost of current genotyping methods, we sought to use PANDAA for screening key HIVDRMs in antiretroviral-naïve individuals at codons 103, 106 and 184 of the HIV-1 reverse transcriptase gene. Mutations selected at these positions have been shown to be the most common driver mutations in treatment failure.  Methods: A total of 103 samples from antiretroviral-naïve individuals previously genotyped by Sanger population sequencing were used to assess and verify the performance of PANDAA. PANDAA samples were run on the ABI 7500 Sequence Detection System to genotype the K103N, V106M and M184V HIVDRMs. In addition, the cost per sample and reaction times were compared. Results: Sanger population sequencing and PANDAA detected K103N mutation in three (2.9%) out of 103 participants.  There was no evidence of baseline V106M and M184V mutations observed in our study. To genotype the six HIVDRMs it costs approximately 40 USD using PANDAA, while the reagents cost per test for Sanger population sequencing is approximately 100 USD per sample. PANDAA was performed quicker compared to Sanger sequencing, 2 hours for PANDAA versus 15 hours for Sanger sequencing. Conclusion: The performance of PANDAA and Sanger population sequencing demonstrated complete concordance. PANDAA could improve patient management by providing quick and relatively cheap access to drug-resistance information.

High incidence of tuberculosis in the first year of antiretroviral therapy in the Botswana National antiretroviral therapy programme between 2011 and 2015.

OBJECTIVE: Tuberculosis (TB) remains one of the leading causes of mortality and morbidity among people living with HIV. We sought to estimate the incidence of TB in a national database of HIV-infected patients receiving antiretroviral therapy (ART) in Botswana. DESIGN: A retrospective analysis of HIV-infected adult patients (≥18years) who initiated ART between 2011 and 2015 in the Botswana ART program. METHODS: Multivariable analysis using Cox regression included sex, age, viral load and CD4 T-cell counts. RESULTS: Of 45 729 patients, with a median follow-up of 1.7 years Q1 : Q3, 0.5, 3.1), 1791 patients developed TB over a median of 1.5 years (Q1 : Q3, 0.3, 3.1) of follow-up (incidence rate 1.9 per 100 person-years; 95% CI 1.8-2.0). At baseline, the median CD4 T-cell count was 272 cells/µl (Q1, Q3 146, 403). The risk of TB was greatest within the first year of ART (incidence rate 2.9 per 100 person-years; 95% CI 2.7-3.1) and in patients with CD4 T-cell counts below 50 cells/µl (incidence rate 8.3/100 person-years; 95% CI 7.1-9.7). Patients with viral loads above 10 000 copies/ml at 3 months post-ART initiation had two times higher risk of TB, hazard ratio 2.5 (95% CI 1.8-2.3). CONCLUSION: We report a high incidence of TB within the first year of ART and in patients with advanced immunodeficiency. Improved screening strategies and virologic monitoring during this early period on ART, coupled with TB preventive treatment, will reduce the burden of TB.

Molecular characterization of hepatitis C virus in liver disease patients in Botswana: a retrospective cross-sectional study.

BACKGROUND: Hepatitis C virus (HCV) infection is a major cause of chronic liver disease globally. Direct acting antivirals (DAAs) have proven effective in curing HCV. However, the current standard of care (SOC) in Botswana remains PEGylated interferon-α (IFN-α) with ribavirin. Several mutations have been reported to confer resistance to interferon-based treatments. Therefore, there is a need to determine HCV genotypes in Botswana, as these data will guide new treatment guidelines and understanding of HCV epidemiology in Botswana. METHODS: This was a retrospective cross-sectional pilot study utilizing plasma obtained from 55 participants from Princess Marina Hospital in Gaborone, Botswana. The partial core region of HCV was amplified, and genotypes were determined using phylogenetic analysis. RESULTS: Four genotype 5a and two genotype 4v sequences were identified. Two significant mutations - K10Q and R70Q - were observed in genotype 5a sequences and have been associated with increased risk of hepatocellular carcinoma (HCC), while R70Q confers resistance to interferon-based treatments. CONCLUSION: Genotypes 5a and 4v are circulating in Botswana. The presence of mutations in genotype 5 suggests that some patients may not respond to IFN-based regimens. The information obtained in this study, in addition to the World health organization (WHO) recommendations, can be utilized by policy makers to implement DAAs as the new SOC for HCV treatment in Botswana.

Self-reported risky sexual practices among adolescents and young adults in Botswana.

BACKGROUND: Adolescents and young adults account for more than one-third of incident Human Immunodeficiency Virus (HIV) infections globally. Understanding sexual practices of this high-risk group is critical in designing HIV targeted prevention programming. OBJECTIVES: To describe self-reported risky sexual practices of adolescents and young adults aged 16-24 years from 30 Botswana communities. METHODS: Cross-sectional, self-reported age at sexual debut; number of sexual partners; condom and alcohol use during sex; intergenerational sex; and transactional sex data were collected. Modified Poisson estimating equations were used to obtain univariate and multivariate-adjusted prevalence ratios (PR) and 95% confidence intervals (CI) comparing engagement in different sexual practices according to gender, accounting for the clustered design of the study. RESULTS: Among the 3380 participants, 2311 reported being sexually active with more females reporting being sexually active compared to males (65% vs. 35%, respectively; p < 0.0001). In univariate analyses, female participants were more likely to report inconsistent condom use (PR 1.61; 95% CI 1.44-1.80), intergenerational sex (PR 9.00; 95% CI 5.84-13.88) and transactional sex (PR 3.46; 95% CI 2.07-5.77) than males, yet less likely to report engaging in sex before age 15 years (PR 0.59; 95% CI: 0.41-0.85), using alcohol around the time of intercourse (PR: 0.59; 95% CI 0.45-0.76) or having ≥ two partners in the last 12 months (PR 0.65; 95% CI 0.57-0.74). CONCLUSIONS: Self-reported risky sexual practices of adolescents and young adults in Botswana differed significantly between males and females. Gender-specific risky sexual practices highlight the importance of developing tailored HIV prevention programming.

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

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

Molecular characterization of hepatitis B virus in blood donors in Botswana.

Hepatitis B virus (HBV) poses a significant threat to blood transfusion safety in sub-Saharan Africa (SSA) where allogeneic blood donations are screened serologically, and more sensitive nucleic acid tests (NATs) are utilized infrequently. HBV strains circulating among blood donors in Botswana are not yet characterized. We designed a cross-sectional study to determine the HBV sub-genotypes and prevalence of hepatitis B surface antigen (HBsAg) among blood donors between November 2014 and October 2015. A total of 12,575 blood donations were screened for HBsAg and 50 consecutive plasma samples were selected for genotyping from confirmed HBsAg+ donations. Overlapping Pol and complete S (Pol/S) open reading frames (ORFs) were sequenced from extracted HBV DNA. To identify any signature amino acids, mutations were compared to sequences from a cohort of chronic HBV patients co-infected with HIV and were treatment naïve. The prevalence of HBsAg+ blood donors was 1.02% (95% CI 0.9-1.2%), and the circulating sub-genotypes were A1 serotype adw2 (36.1%), D2 serotype ayw2 (2.9%), and D3 serotypes ayw 1/2 (58.3%). Prevalence of escape mutations was 14% from HBV isolates of blood donors and 15% from isolates of HBV/HIV co-infected patients (p = 0.6926). The escape mutations sP120L, sG130R, sY134H, and sD144A were identified predominantly among HBV isolates from blood donors. These escape mutations have been associated with accelerated HBV sequelae [e.g., liver cirrhosis (LC) and hepatocellular carcinoma (HCC)], failure to detect HBsAg, inability to respond to immunoglobulin (Ig) therapy, and HBV vaccine escape. Characterizing the HBV burden, circulating sub-genotypes, and clinically relevant mutations among blood donors in Botswana is important to elucidate the efficacy of currently available vaccines, predicting HBV-transmission patterns, understanding the cohort's risk to HBV-related complications, and to developing prevention strategies and effective genotype-based antiretroviral therapies.

Cross-sectional estimates revealed high HIV incidence in Botswana rural communities in the era of successful ART scale-up in 2013-2015.

BACKGROUND: Botswana is close to reaching the UNAIDS "90-90-90" HIV testing, antiretroviral treatment (ART), and viral suppression goals. We sought to determine HIV incidence in this setting with both high HIV prevalence and high ART coverage. METHODS: We used a cross-sectional approach to assessing HIV incidence. A random, population-based sample of adults age 16-64 years was enrolled in 30 rural and peri-urban communities as part of the Botswana Combination Prevention Project (BCPP), from October 2013 -November 2015. Data and samples from the baseline household survey were used to estimate cross-sectional HIV incidence, following an algorithm that combined Limiting-Antigen Avidity Assay (LAg-Avidity EIA), ART status (documented or by testing ARV drugs in plasma) and HIV-1 RNA load. The LAg-Avidity EIA cut-off normalized optical density (ODn) was set at 1.5. The HIV-1 RNA cut-off was set at 400 copies/mL. For estimation purposes, the Mean Duration of Recent Infection was 130 days and the False Recent Rate (FRR) was evaluated at values of either 0 or 0.39%. RESULTS: Among 12,610 individuals participating in the baseline household survey, HIV status was available for 12,570 participants and 3,596 of them were HIV positive. LAg-Avidity EIA data was generated for 3,581 (99.6%) of HIV-positive participants. Of 326 participants with ODn ≤1.5, 278 individuals were receiving ART verified through documentation and were considered to represent longstanding HIV infections. Among the remaining 48 participants who reported no use of ART, 14 had an HIV-1 RNA load ≤400 copies/mL (including 3 participants with ARVs in plasma) and were excluded, as potential elite/viremic controllers or undisclosed ART. Thus, 34 LAg-Avidity-EIA-recent, ARV-naïve individuals with detectable HIV-1 RNA (>400 copies/mL) were classified as individuals with recent HIV infections. The annualized HIV incidence among 16-64 year old adults was estimated at 1.06% (95% CI 0.68-1.45%) with zero FRR, and at 0.64% (95% CI 0.24-1.04%) using a previously defined FRR of 0.39%. Within a subset of younger individuals 16-49 years old, the annualized HIV incidence was estimated at 1.29% (95% CI 0.82-1.77%) with zero FRR, and at 0.90% (95% CI 0.42-1.38%) with FRR set to 0.39%. CONCLUSIONS: Using a cross-sectional estimate of HIV incidence from 2013-2015, we found that at the time of near achievement of the UNAIDS 90-90-90 targets, ~1% of adults (age 16-64 years) in Botswana's rural and peri-urban communities became HIV infected annually.

Molecular Characterization of Near Full-Length Genomes of Hepatitis B Virus Isolated from Predominantly HIV Infected Individuals in Botswana.

The World Health Organization plans to eliminate hepatitis B and C Infections by 2030. Therefore, there is a need to study and understand hepatitis B virus (HBV) epidemiology and viral evolution further, including evaluating occult (HBsAg-negative) HBV infection (OBI), given that such infections are frequently undiagnosed and rarely treated. We aimed to molecularly characterize HBV genomes from 108 individuals co-infected with human immunodeficiency virus (HIV) and chronic hepatitis B (CHB) or OBI identified from previous HIV studies conducted in Botswana from 2009 to 2012. Full-length (3.2 kb) and nearly full-length (~3 kb) genomes were amplified by nested polymerase chain reaction (PCR). Sequences from OBI participants were compared to sequences from CHB participants and GenBank references to identify OBI-unique mutations. HBV genomes from 50 (25 CHB and 25 OBI) individuals were successfully genotyped. Among OBI participants, subgenotype A1 was identified in 12 (48%), D3 in 12 (48%), and E in 1 (4%). A similar genotype distribution was observed in CHB participants. Whole HBV genome sequences from Botswana, representing OBI and CHB, were compared for the first time. There were 43 OBI-unique mutations, of which 26 were novel. Future studies using larger sample sizes and functional analysis of OBI-unique mutations are warranted.