Search

We assessed the performance and clinical relevance of Illumina MiSeq next-generation sequencing (NGS) for HIV-1 genotyping compared with Sanger sequencing (SS). We analyzed 167 participants, 45 with virologic failure (VL ≥ 1000 copies/mL), i.e., cases, and 122 time-matched participants with virologic suppression (VL < 1000 copies/mL), i.e., controls, 12 months post-ART initiation. Major surveillance drug resistance mutations (SDRMs) detected by SS were all detectable by NGS. Among cases at 12 months, SS identified SDRMs in 32/45 (71.1%) while NGS identified SDRMs among 35/45 (77.8%), increasing the number of cases with SDRMs by 3/45 (6.7%). Participants identified with, and proportions of major SDRMs increased when NGS was used. NGS vs. SS at endpoint revealed for NNRTIs: 36/45 vs. 33/45; Y181C: 26/45 vs. 24/45; K103N: 9/45 vs. 6/45 participants with SDRMs, respectively. At baseline, NGS revealed major SDRMs in 9/45 (20%) cases without SDRMs by SS. Participant MBL/043, among the nine, the following major SDRMs existed: L90M to PIs, K65R and M184V to NRTIs, and Y181C and K103N to NNRTIs. The SDRMs among the nine increased SDRMs to NRTIs, NNRTIs, and PIs. Only 43/122 (25.7%) of participants had pre-treatment minority SDRMs. Also, 24.4% of the cases vs. 26.2 of controls had minority SDRMs (p = 0.802); minority SDRMs were not associated with virologic failure. NGS agreed with SS in HIV-1 genotyping but detected additional major SDRMs and identified more participants harboring major SDRMs, expanding the HIV DRM profile of this cohort. NGS could improve HIV genotyping to guide treatment decisions for enhancing ART efficacy, a cardinal pre-requisite in the pursuit of the UNAIDS 95-95-95 targets.

The Sustainable East Africa Research in Community Health (SEARCH) trial was a universal test-and-treat (UTT) trial in rural Uganda and Kenya, aiming to lower regional HIV-1 incidence. Here, we quantify breakthrough HIV-1 transmissions occurring during the trial from population-based, dried blood spot samples. Between 2013 and 2017, we obtained 549 gag and 488 pol HIV-1 consensus sequences from 745 participants: 469 participants infected prior to trial commencement and 276 SEARCH-incident infections. Putative transmission clusters, with a 1.5% pairwise genetic distance threshold, were inferred from maximum likelihood phylogenies; clusters arising after the start of SEARCH were identified with Bayesian time-calibrated phylogenies. Our phylodynamic approach identified nine clusters arising after the SEARCH start date: eight pairs and one triplet, representing mostly opposite-gender linked (6/9), within-community transmissions (7/9). Two clusters contained individuals with non-nucleoside reverse transcriptase inhibitor (NNRTI) resistance, both linked to intervention communities. The identification of SEARCH-incident, within-community transmissions reveals the role of unsuppressed individuals in sustaining the epidemic in both arms of a UTT trial setting. The presence of transmitted NNRTI resistance, implying treatment failure to the efavirenz-based antiretroviral therapy (ART) used during SEARCH, highlights the need to improve delivery and adherence to up-to-date ART recommendations, to halt HIV-1 transmission.

Detailed characterization of transmitted HIV-1 variants in Uganda is fundamentally important to inform vaccine design, yet studies on the transmitted full-length strains of subtype D viruses are limited. Here, we amplified single genomes and characterized viruses, some of which were previously classified as subtype D by sub-genomic pol sequencing that were transmitted in Uganda between December 2006 to June 2011. Analysis of 5′ and 3′ half genome sequences showed 73% (19/26) of infections involved single virus transmissions, whereas 27% (7/26) of infections involved multiple variant transmissions based on predictions of a model of random virus evolution. Subtype analysis of inferred transmitted/founder viruses showed a high transmission rate of inter-subtype recombinants (69%, 20/29) involving mainly A1/D, while pure subtype D variants accounted for one-third of infections (31%, 9/29). Recombination patterns included a predominance of subtype D in the gag/pol region and a highly recombinogenic envelope gene. The signal peptide-C1 region and gp41 transmembrane domain (Tat2/Rev2 flanking region) were hotspots for A1/D recombination events. Analysis of a panel of 14 transmitted/founder molecular clones showed no difference in replication capacity between subtype D viruses (n = 3) and inter-subtype mosaic recombinants (n = 11). However, individuals infected with high replication capacity viruses had a faster CD4 T cell loss. The high transmission rate of unique inter-subtype recombinants is striking and emphasizes the extraordinary challenge for vaccine design and, in particular, for the highly variable and recombinogenic envelope gene, which is targeted by rational designs aimed to elicit broadly neutralizing antibodies.

Phylogenetic inference is useful in characterising HIV transmission networks and assessing where prevention is likely to have the greatest impact. However, estimating parameters that influence the network structure is still scarce, but important in evaluating determinants of HIV spread. We analyzed 2017 HIV pol sequences (728 Lake Victoria fisherfolk communities (FFCs), 592 female sex workers (FSWs) and 697 general population (GP)) to identify transmission networks on Maximum Likelihood (ML) phylogenetic trees and refined them using time-resolved phylogenies. Network generative models were fitted to the observed degree distributions and network parameters, and corrected Akaike Information Criteria and Bayesian Information Criteria values were estimated. 347 (17.2%) HIV sequences were linked on ML trees (maximum genetic distance ≤4.5%, ≥95% bootstrap support) and, of these, 303 (86.7%) that consisted of pure A1 (n = 168) and D (n = 135) subtypes were analyzed in BEAST v1.8.4. The majority of networks (at least 40%) were found at a time depth of ≤5 years. The waring and yule models fitted best networks of FFCs and FSWs respectively while the negative binomial model fitted best networks in the GP. The network structure in the HIV-hyperendemic FFCs is likely to be scale-free and shaped by preferential attachment, in contrast to the GP. The findings support the targeting of interventions for FFCs in a timely manner for effective epidemic control. Interventions ought to be tailored according to the dynamics of the HIV epidemic in the target population and understanding the network structure is critical in ensuring the success of HIV prevention programs.

The ability to efficiently establish a new infection is a critical property for human immunodeficiency virus type 1 (HIV-1). Although the envelope protein of the virus plays an essential role in receptor binding and internalization of the infecting virus, the structural proteins, the polymerase and the assembly of new virions may also play a role in establishing and spreading viral infection in a new host. We examined Ugandan viruses from newly infected patients and focused on the contribution of the Gag-Pol genes to replication capacity. A panel of Gag-Pol sequences generated using single genome amplification from incident HIV-1 infections were cloned into a common HIV-1 NL4.3 pol/env backbone and the influence of Gag-Pol changes on replication capacity was monitored. Using a novel protein domain approach, we then documented diversity in the functional protein domains across the Gag-Pol region and identified differences in the Gag-p6 domain that were frequently associated with higher in vitro replication.

The General Population Cohort (GPC) in south-western Uganda has a low HIV-1 incidence rate (pol sequences (GPC: n = 1418, non-GPC sites: n = 1223, Central Uganda: n = 1010 and Eastern Uganda: n = 145) generated between 2003–2015 were analysed using phylogenetic methods with demographic data to understand HIV-1 transmission in this cohort and inform the epidemic response. HIV-1 subtype A1 was the most prevalent strain in the GPC area (GPC and non-GPC sites) (39.8%), central (45.9%) and eastern (52.4%) Uganda. However, in the GPC alone, subtype D was the predominant subtype (39.1%). Of the 524 transmission clusters identified by Cluster Picker, all large clusters (≥5 individuals, n = 8) involved individuals from the GPC. In a multivariate analysis, clustering was strongly associated with being female (adjusted Odds Ratio, aOR = 1.28; 95% CI, 1.06–1.54), being >25 years (aOR = 1.52; 95% CI, 1.16–2.0) and being a resident in the GPC (aOR = 6.90; 95% CI, 5.22–9.21). Phylogeographic analysis showed significant viral dissemination (Bayes Factor test, BF > 3) from the GPC without significant viral introductions (BF < 3) into the GPC. The findings suggest localized HIV-1 transmission in the GPC. Intensifying geographically focused combination interventions in the GPC would contribute towards controlling HIV-1 infections.

Across sub-Saharan Africa, key populations with elevated HIV-1 incidence and/or prevalence have been identified, but their contribution to disease spread remains unclear. We performed viral deep-sequence phylogenetic analyses to quantify transmission dynamics between the general population (GP), fisherfolk communities (FF), and women at high risk of infection and their clients (WHR) in central and southwestern Uganda. Between August 2014 and August 2017, 6185 HIV-1 positive individuals were enrolled in 3 GP and 10 FF communities, 3 WHR enrollment sites. A total of 2531 antiretroviral therapy (ART) naïve participants with plasma viral load >1000 copies/mL were deep-sequenced. One hundred and twenty-three transmission networks were reconstructed, including 105 phylogenetically highly supported source–recipient pairs. Only one pair involved a WHR and male participant, suggesting that improved population sampling is needed to assess empirically the role of WHR to the transmission dynamics. More transmissions were observed from the GP communities to FF communities than vice versa, with an estimated flow ratio of 1.56 (95% CrI 0.68–3.72), indicating that fishing communities on Lake Victoria are not a net source of transmission flow to neighboring communities further inland. Men contributed disproportionally to HIV-1 transmission flow regardless of age, suggesting that prevention efforts need to better aid men to engage with and stay in care.