Deogratius Ssemwanga

Geographic and Temporal Trends in the Molecular Epidemiology and Genetic Mechanisms of Transmitted HIV-1 Drug Resistance: An Individual-Patient- and Sequence-Level Meta-Analysis

Background Regional and subtype-specific mutational patterns of HIV-1 transmitted drug resistance (TDR) are essential for informing first-line antiretroviral (ARV) therapy guidelines and designing diagnostic assays for use in regions where standard genotypic resistance testing is not affordable. We sought to understand the molecular epidemiology of TDR and to identify the HIV-1 drug-resistance mutations responsible for TDR in different regions and virus subtypes. Methods and Findings We reviewed all GenBank submissions of HIV-1 reverse transcriptase sequences with or without protease and identified 287 studies published between March 1, 2000, and December 31, 2013, with more than 25 recently or chronically infected ARV-naïve individuals. These studies comprised 50,870 individuals from 111 countries. Each set of study sequences was analyzed for phylogenetic clustering and the presence of 93 surveillance drug-resistance mutations (SDRMs). The median overall TDR prevalence in sub-Saharan Africa (SSA), south/southeast Asia (SSEA), upper-income Asian countries, Latin America/Caribbean, Europe, and North America was 2.8%, 2.9%, 5.6%, 7.6%, 9.4%, and 11.5%, respectively. In SSA, there was a yearly 1.09-fold (95% CI: 1.05–1.14) increase in odds of TDR since national ARV scale-up attributable to an increase in non-nucleoside reverse transcriptase inhibitor (NNRTI) resistance. The odds of NNRTI-associated TDR also increased in Latin America/Caribbean (odds ratio [OR] = 1.16; 95% CI: 1.06–1.25), North America (OR = 1.19; 95% CI: 1.12–1.26), Europe (OR = 1.07; 95% CI: 1.01–1.13), and upper-income Asian countries (OR = 1.33; 95% CI: 1.12–1.55). In SSEA, there was no significant change in the odds of TDR since national ARV scale-up (OR = 0.97; 95% CI: 0.92–1.02). An analysis limited to sequences with mixtures at less than 0.5% of their nucleotide positions—a proxy for recent infection—yielded trends comparable to those obtained using the complete dataset. Four NNRTI SDRMs—K101E, K103N, Y181C, and G190A—accounted for >80% of NNRTI-associated TDR in all regions and subtypes. Sixteen nucleoside reverse transcriptase inhibitor (NRTI) SDRMs accounted for >69% of NRTI-associated TDR in all regions and subtypes. In SSA and SSEA, 89% of NNRTI SDRMs were associated with high-level resistance to nevirapine or efavirenz, whereas only 27% of NRTI SDRMs were associated with high-level resistance to zidovudine, lamivudine, tenofovir, or abacavir. Of 763 viruses with TDR in SSA and SSEA, 725 (95%) were genetically dissimilar; 38 (5%) formed 19 sequence pairs. Inherent limitations of this study are that some cohorts may not represent the broader regional population and that studies were heterogeneous with respect to duration of infection prior to sampling. Conclusions Most TDR strains in SSA and SSEA arose independently, suggesting that ARV regimens with a high genetic barrier to resistance combined with improved patient adherence may mitigate TDR increases by reducing the generation of new ARV-resistant strains. A small number of NNRTI-resistance mutations were responsible for most cases of high-level resistance, suggesting that inexpensive point-mutation assays to detect these mutations may be useful for pre-therapy screening in regions with high levels of TDR. In the context of a public health approach to ARV therapy, a reliable point-of-care genotypic resistance test could identify which patients should receive standard first-line therapy and which should receive a protease-inhibitor-containing regimen.

Molecular identification, cloning and characterization of transmitted/founder HIV-1 subtype A, D and A/D infectious molecular clones

We report the molecular identification, cloning and initial biological characterization of 12 full-length HIV-1 subtype A, D and A/D recombinant transmitted/founder (T/F) genomes. T/F genomes contained intact canonical open reading frames and all T/F viruses were replication competent in primary human T-cells, although subtype D virus replication was more efficient (p<0.05). All 12 viruses utilized CCR5 but not CXCR4 as a co-receptor for entry and exhibited a neutralization profile typical of tier 2 primary virus strains, with significant differences observed between subtype A and D viruses with respect to sensitivity to monoclonal antibodies VRC01, PG9 and PG16 and polyclonal subtype C anti-HIV IgG (p<0.05 for each). The present report doubles the number of T/F HIV-1 clones available for pathogenesis and vaccine research and extends their representation to include subtypes A, B, C and D.

Multiple HIV-1 infections with evidence of recombination in heterosexual partnerships in a low risk rural clinical cohort in Uganda.

We report on the frequency of multiple infections, generation of recombinants and consequences on disease progression in 35 HIV-1 infected individuals from 7 monogamous and 6 polygamous partnerships within a Rural Clinical Cohort in Uganda. The env-C2V3, gag-p24 and pol-IN genes were sequenced. Single genome amplified half genome sequences were used to map recombination breakpoints. Three participants were dually infected with subtypes A and D, one case with subtype A and A/D recombinant and the fifth with 2 phylogenetically distinct A/D recombinants. Occurrence of A/D recombination was observed in two multiple infected individuals. Rate of late stage WHO events using Cox regression was 3 times greater amongst multiple infected compared to singly infected individuals (hazard ratio 3.35; 95% CI 1.09, 10.3; p = 0.049). We have shown that polygamous relationships involving subtype discordant partnerships was a major contributor of multiple infections with generation of inter subtype recombinants in our cohort.

Rates of HIV-1 superinfection and primary HIV-1 infection are similar in female sex workers in Uganda

Objective:To determine and compare the rates of HIV superinfection and primary HIV infection in high-risk female sex workers (FSWs) in Kampala, Uganda. Design:A retrospective analysis of individuals who participated in a clinical cohort study among high-risk FSWs in Kampala, Uganda. Methods:Plasma samples from HIV-infected FSWs in Kampala, Uganda were examined with next-generation sequencing of the p24 and gp41HIV genomic regions for the occurrence of superinfection. Primary HIV incidence was determined from initially HIV-uninfected FSWs from the same cohort, and incidence rate ratios were compared. Results:The rate of superinfection in these women (7/85; 3.4/100 person-years) was not significantly different from the rate of primary infection in the same population (3.7/100 person-years; incidence rate ratio = 0.91, P = 0.42). Seven women also entered the study dual-infected (16.5% either dual or superinfected). The women with any presence of dual infection were more likely to report sex work as their only source of income (P = 0.05), and trended to be older and more likely to be widowed (P = 0.07). Conclusions:In this cohort of FSWs, HIV superinfection occurred at a high rate and was similar to that of primary HIV infection. These results differ from a similar study of high-risk female bar workers in Kenya that found the rate of superinfection to be significantly lower than the rate of primary HIV infection.

Analysis of the history and spread of HIV-1 in Uganda using phylodynamics

HIV prevalence has decreased in Uganda since the 1990s, but remains substantial within high-risk groups. Here, we reconstruct the history and spread of HIV subtypes A1 and D in Uganda and explore the transmission dynamics in high-risk populations. We analysed HIV pol sequences from female sex workers in Kampala (n = 42), Lake Victoria fisher-folk (n = 46) and a rural clinical cohort (n = 74), together with publicly available sequences from adjacent regions in Uganda (n = 412) and newly generated sequences from samples taken in Kampala in 1986 (n = 12). Of the sequences from the three Ugandan populations, 60 (37.1 %) were classified as subtype D, 54 (33.3 %) as subtype A1, 31 (19.1 %) as A1/D recombinants, six (3.7 %) as subtype C, one (0.6 %) as subtype G and 10 (6.2 %) as other recombinants. Among the A1/D recombinants we identified a new candidate circulating recombinant form. Phylodynamic and phylogeographic analyses using BEAST indicated that the Ugandan epidemics originated in 1960 (1950–1968) for subtype A1 and 1973 (1970–1977) for D, in rural south-western Uganda with subsequent spread to Kampala. They also showed extensive interconnection with adjacent countries. The sequence analysis shows both epidemics grew exponentially during the 1970s–1980s and decreased from 1992, which agrees with HIV prevalence reports in Uganda. Inclusion of sequences from the 1980s indicated the origin of both epidemics was more recent than expected and substantially narrowed the confidence intervals in comparison to previous estimates. We identified three transmission clusters and ten pairs, none of them including patients from different populations, suggesting active transmission within a structured transmission network.

Effect of HIV-1 Subtypes on Disease Progression in Rural Uganda: A Prospective Clinical Cohort Study

Objective We examined the association of HIV-1 subtypes with disease progression based on three viral gene regions. Design A prospective HIV-1 clinical cohort study in rural Uganda. Methods Partial gag, env and pol genes were sequenced. Cox proportional hazard regression modelling was used to estimate adjusted hazard ratios (aHRs) of progression to: CD4≤250, AIDS onset and death, adjusted for sex, age and CD4 count at enrolment. Results Between 1990 and 2010, 292 incident cases were subtyped: 25% had subtype A, 45% had D, 26% had A/D recombinants, 1% had C and 4% were other recombinant forms. Of the 278 incident cases included in the disease progression analysis, 62% progressed to CD4≤250, 32% to AIDS, and 34% died with a higher proportion being among subtype D cases. The proportions of individuals progressing to the three endpoints were significantly higher among individuals infected with subtype D. Throughout the study period, individuals infected with subtype D progressed faster to CD4≤250, adjusted HR (aHR), (95% CI) = 1.72 (1.16–2.54), but this was mainly due to events in the period before antiretroviral therapy (ART) introduction, when individuals infected with subtype D significantly progressed faster to CD4≤250 than subtype A cases; aHR (95% CI) = 1.78 (1.01–3.14). Conclusions In this population, HIV-1 subtype D was the most prevalent and was associated with faster HIV-1 disease progression than subtype A. Further studies are needed to examine the effect of HIV-1 subtypes on disease progression in the ART period and their effect on the virological and immunological ART outcomes.

Low Drug Resistance Levels Among Drug-Naive Individuals with Recent HIV Type 1 Infection in a Rural Clinical Cohort in Southwestern Uganda

To investigate the prevalence of transmitted drug resistance (TDR) among individuals with recent HIV-1 infection between February 2004 and January 2010 in a rural clinical cohort, samples from 72 participants were analyzed. Results from the 72 participants showed no protease inhibitor and nucleoside reverse transcriptase inhibitor-associated mutations. One participant (1.4%, 95% CI: 0.04-7.5%) had two nonnucleoside reverse transcriptase inhibitor mutations (G190E and P225H). HIV-1 subtype frequencies were A 22 (30.6%), D 38 (52.8%), and C 1 (1.4%); 11 (15.3%) were A/D unique recombinant forms. Seven years after the scale up of antiretroviral therapy (ART) in a rural clinical cohort in Uganda, the prevalence of TDR among recently HIV-1-infected individuals was low at 1.4%. Since our findings from an HIV study cohort may not be generalizable to the general population, routine TDR surveys in specific populations may be necessary to inform policy on the magnitude and prevention strategies of TDR.

Transmitted Antiretroviral Drug Resistance Among Drug-Naive Female Sex Workers With Recent Infection in Kampala, Uganda

During 2006-2007, transmitted human immunodeficiency virus (HIV) drug resistance (TDR) among drug-naive women with newly diagnosed HIV infection and likely to be recently infected when attending antenatal clinics in Entebbe was found to be <5% with use of the World Health Organization (WHO) survey method. Using the same method, we attempted to classify TDR among women who seroconverted during 2008-2010 and who were identified from a cohort of recently infected sex workers in Kampala, Uganda. TDR mutations were identified using the 2009 WHO TDR mutations list. The WHO survey method could not be used to classify TDR because the necessary sample size was not reached during the survey period. However, a point prevalence estimate of 2.6% (95% confidence interval, 0.07%-13.8%) nonnucleoside reverse-transcriptase inhibitor TDR was determined.

Using nearly full-genome HIV sequence data improves phylogeny reconstruction in a simulated epidemic

HIV molecular epidemiology studies analyse viral pol gene sequences due to their availability, but whole genome sequencing allows to use other genes. We aimed to determine what gene(s) provide(s) the best approximation to the real phylogeny by analysing a simulated epidemic (created as part of the PANGEA_HIV project) with a known transmission tree. We sub-sampled a simulated dataset of 4662 sequences into different combinations of genes (gag-pol-env, gag-pol, gag, pol, env and partial pol) and sampling depths (100%, 60%, 20% and 5%), generating 100 replicates for each case. We built maximum-likelihood trees for each combination using RAxML (GTR + Γ), and compared their topologies to the corresponding true tree’s using CompareTree. The accuracy of the trees was significantly proportional to the length of the sequences used, with the gag-pol-env datasets showing the best performance and gag and partial pol sequences showing the worst. The lowest sampling depths (20% and 5%) greatly reduced the accuracy of tree reconstruction and showed high variability among replicates, especially when using the shortest gene datasets. In conclusion, using longer sequences derived from nearly whole genomes will improve the reliability of phylogenetic reconstruction. With low sample coverage, results can be highly variable, particularly when based on short sequences.

Deciphering the Complex Distribution of Human Immunodeficiency Virus Type 1 Subtypes among Different Cohorts in Northern Tanzania

Background Increased understanding of the genetic diversity of HIV-1 is challenging but important in the development of an effective vaccine. We aimed to describe the distribution of HIV-1 subtypes in northern Tanzania among women enrolled in studies preparing for HIV-1 prevention trials (hospitality facility-worker cohorts), and among men and women in an open cohort demographic surveillance system (Kisesa cohort). Methods The polymerase encompassing partial reverse transcriptase was sequenced and phylogenetic analysis performed and subtype determined. Questionnaires documented demographic data. We examined factors associated with subtype using multinomial logistic regression, adjusted for study, age, and sex. Results Among 140 individuals (125 women and 15 men), subtype A1 predominated (54, 39%), followed by C (46, 33%), D (25, 18%) and unique recombinant forms (URFs) (15, 11%). There was weak evidence to suggest different subtype frequencies by study (for example, 18% URFs in the Kisesa cohort versus 5–9% in the hospitality facility-worker cohorts; adjusted relative-risk ratio (aRR) = 2.35 [95% CI 0.59,9.32]; global p = 0.09). Compared to men, women were less likely to have subtype D versus A (aRR = 0.12 [95% CI 0.02,0.76]; global p = 0.05). There was a trend to suggest lower relative risk of subtype D compared to A with older age (aRR = 0.44 [95% CI 0.23,0.85] per 10 years; global p = 0.05). Conclusions We observed multiple subtypes, confirming the complex genetic diversity of HIV-1 strains circulating in northern Tanzania, and found some differences between cohorts and by age and sex. This has important implications for vaccine design and development, providing opportunity to determine vaccine efficacy in diverse HIV-1 strains.