Mapping of HIV-1C Transmission Networks Reveals Extensive Spread of Viral Lineages Across Villages in Botswana Treatment-as-Prevention Trial.

BACKGROUND: Phylogenetic mapping of HIV-1 lineages circulating across defined geographical locations is promising for better understanding HIV transmission networks to design optimal prevention interventions. METHODS: We obtained near full-length HIV-1 genome sequences from people living with HIV (PLWH), including participants on antiretroviral treatment in the Botswana Combination Prevention Project, conducted in 30 Botswana communities in 2013-2018. Phylogenetic relationships among viral sequences were estimated by maximum likelihood. RESULTS: We obtained 6078 near full-length HIV-1C genome sequences from 6075 PLWH. We identified 984 phylogenetically distinct HIV-1 lineages (molecular HIV clusters) circulating in Botswana by mid-2018, with 2-27 members per cluster. Of these, dyads accounted for 62%, approximately 32% (n = 316) were found in single communities, and 68% (n = 668) were spread across multiple communities. Men in clusters were approximately 3 years older than women (median age 42 years, vs 39 years; P < .0001). In 65% of clusters, men were older than women, while in 35% of clusters women were older than men. The majority of identified viral lineages were spread across multiple communities. CONCLUSIONS: A large number of circulating phylogenetically distinct HIV-1C lineages (molecular HIV clusters) suggests highly diversified HIV transmission networks across Botswana communities by 2018.

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.

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