Secondary Logo

Journal Logo

Prevalence and clinical impacts of HIV-1 intersubtype recombinants in Uganda revealed by near-full-genome population and deep sequencing approaches

Lee, Guinevere Q.a,*; Bangsberg, David R.b,c,d,e,f; Mo, Theresaa; Lachowski, Chrisa; Brumme, Chanson J.a; Zhang, Wendya; Lima, Viviane D.a; Boum, Yap IIb; Mwebesa, Bosco Bwanab; Muzoora, Conradb; Andia, Irenb; Mbalibulha, Yonab; Kembabazi, Annetb; Carroll, Ryanb,c,e; Siedner, Mark J.b,c,e; Haberer, Jessica E.e; Mocello, A. Raing; Kigozi, Simone H.b; Hunt, Peter W.g; Martin, Jeffrey N.g; Harrigan, P. Richarda,*

doi: 10.1097/QAD.0000000000001619
CLINICAL SCIENCE
Free
SDC

Objectives: HIV-1 subtypes A1 and D cocirculate in a rural community in Mbarara, Uganda. This study examines HIV-1 intersubtype recombination in this community under a full-genome sequencing context. We aim to estimate prevalence, examine time trends, and test for clinical correlates and outcomes associated with intersubtype recombinants.

Methods: Near-full-genome HIV-1 Sanger sequence data were collected from plasma samples of 504 treatment-naïve individuals, who then received protease inhibitor or nonnucleoside reverse transcriptase inhibitor-containing regimens and were monitored for up to 7.5 years. Subtypes were inferred by Los Alamos Recombinant Identification Program (RIP) 3.0 and compared with Sanger/REGA and MiSeq/RIP. ‘Nonrecombinants’ and ‘recombinants’ infections were compared in terms of pretherapy viral load, CD4+ cell count, posttherapy time to virologic suppression, virologic rebound, first CD4+ rise above baseline and sustained CD4+ recovery.

Results: Prevalence of intersubtype recombinants varied depending on the genomic region examined: gag (15%), prrt (11%), int (8%), vif (10%), vpr (2%), vpu (9%), GP120 (8%), GP41 (18%), and nef (4%). Of the 200 patients with near-full-genome data, prevalence of intersubtype recombination was 46%; the most frequently observed recombinant was A1-D (25%). Sanger/REGA and MiSeq/RIP yielded generally consistent results. Phylogenetic tree revealed most recombinants did not share common ancestors. No temporal trend was observed (all P > 0.1). Subsequent subtype switches were detected in 27 of 143 (19%) study participants with follow-up sequences. Nonrecombinant versus recombinants infections were not significantly different in any pre nor posttherapy clinical correlates examined (all P > 0.2).

Conclusion: Intersubtype recombination was highly prevalent (46%) in Uganda if the entire HIV genome was considered, but was neither associated with clinical correlates nor therapy outcomes.

aBC Centre for Excellence in HIV/AIDS, Vancouver, British Columbia, Canada

bMbarara University of Science and Technology, Mbarara, Uganda

cHarvard Medical School

dHarvard School of Public Health

eMassachusetts General Hospital, Boston, Massachusetts

fOregon Health Sciences University, Portland, Oregon

gUniversity of California, San Francisco, California, USA.

*P. Richard Harrigan and Guinevere Q. Lee equally contributed to this article.

Correspondence to P. Richard Harrigan, PhD, BC Centre for Excellence in HIV/AIDS, 603-1081 Burrard Street, Vancouver, BC V6Z 1Y6, Canada. Tel: +604 806 8281; fax: +604 806 8464; e-mail: prharrigan@cfenet.ubc.ca

Received 8 March, 2017

Revised 28 June, 2017

Accepted 29 June, 2017

Supplemental digital content is available for this article. Direct URL citations appear in the printed text and are provided in the HTML and PDF versions of this article on the journal's Website (http://www.AIDSonline.com).

Back to Top | Article Outline

Introduction

HIV-1 Group M, which currently dominates the global epidemic, is classified into subtypes (or clades) A1, A2, A3, A4, B, C, D, F1, F2, G, H, J, K and various circulating recombinant forms such as AB, AE, and AG [1]. Because the virus’ life cycle involves packing two full-length RNA genomes during viral particles assembly, it provides chances for an event called ‘template switching’ during reverse transcription, leading to the generation of recombinant daughter genomes that contains portions of the parental templates [2].

It has been estimated that HIV-1 circulating recombinant forms and unique recombinant forms are currently responsible for 18–20% of the infections worldwide, and are especially prominent in African, Asian, and South American countries where multiple subtypes cocirculate [3]. However, most studies that report the prevalence of intersubtype recombinant viruses examine only a specific part of the HIV-1 genome, mainly from pol because of its usage in drug resistance testing and its availability in public databases.

Relevant to this study, prevalence of A1-D recombinants in rural Uganda was estimated to be 6–19% based on pol and/or GP41 sequences conglomerated from multiple cohorts and sequence database in multiple studies [4–6]. A small-scale study examined near-full-genome data from 46 patients in Rakai, Uganda and reported 30% recombinants [7]. Another study reported a 30% prevalence of AD recombinants by pol and its enrichment in severely septic patients [8]. Little is known about pre and posttreatment virologic and immunologic impacts associated with infections by these recombinant viruses.

The Uganda AIDS Rural Treatment Outcomes cohort consists of over 500 HIV-infected patients in Mbarara, Uganda, where HIV-1 subtypes A1 and D cocirculate [9–11]. Plasma samples were available before and after treatment initiation, and virologic and immunologic outcome data were collected for over 7 years posttherapy. As such, this cohort provides an excellent opportunity to observe the natural prevalence and clinical impact of HIV recombinants.

The objective of this study is to estimate prevalence, examine time trends, and test for clinical correlates and outcomes associated with infections by HIV-1 intersubtype recombinants. We hypothesize that prevalence of recombinants is higher than previously reported if examined under a near-full-genome-sequencing context, and that infection with HIV-1 recombinants, compared with nonrecombinants, is associated with negative pretherapy clinical correlates and inferior posttherapy virologic and immunologic responses.

Back to Top | Article Outline

Methods

Ethics statement

The study was approved by the Mbarara University of Science and Technology Human Subjects Committee and Partners Healthcare Human Subjects Committee, the Uganda Council of Science and Technology, the University of British Columbia/ Providence Healthcare Research Ethics Board (H11–01642), and the University of California Human Research Subjects Committee. All participants provided written informed consent.

Back to Top | Article Outline

Cohort description

The Uganda AIDS Rural Treatment Outcomes [12–14] is a cross-sectional cohort of 504 initially treatment-naïve HIV-1 infected study participants. They were followed primarily at the Immune Suppression Syndrome Clinic in Mbarara, Uganda, a rural community 4.5 h by automobile from the capital city of Kampala. Study participants were enrolled just before the start of antiretroviral regimen from 27 June 2005 to 8 April 2010, and were longitudinally followed every 3–6 months to receive viral load and CD4+ cell count monitoring for up to 7.5 years until 11 January 2013 or until lost to follow-up. Among the 504 participants, 296 posttherapy follow-up samples were available for HIV-1 RNA sequencing from 143 unique individuals.

Back to Top | Article Outline

HIV-1 near-full-genome PCR amplification

Total nucleic acid was extracted from 500 μl of plasma samples using NucliSENS easyMag (bioMérieux, Cambridge, Massachusetts, USA). All PCR and sequencing primer sequences and thermocycler methods are listed in Supplementary Tables 1–3, http://links.lww.com/QAD/B149. Briefly, reverse transcription and nested-PCR reactions were performed using a three-reaction five-amplicon approach with near-full-genome coverage. The five overlapping amplicons covered gag to protease (HXB2 coordinate 680–2724), protease to reverse transcriptase (2011–3798), reverse transcriptase to vpu (3626–5980), vpr to GP120 (5549–7760), and GP41 to nef (7652–9610).

Back to Top | Article Outline

Sanger sequencing

Bulk sequencing was performed on ABI 3730 DNA Sequencer using BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Foster City, California, USA). Chromatograms were aligned against HXB2 reference sequences by our in-house automated alignment and base-calling program RECall [15], resulting in seven sequence databases for each of gag, prrt, int, vif, vpr, vpu, partial GP120, GP41, and nef, which were concatenated to create near-full-genome sequences. For more details on alignment, concatenation and quality control, refer to Supplementary Material and Methods, http://links.lww.com/QAD/B149. Los Alamos PhyML 3.0 was used to construct phylogenetic trees from these concatenated near-full-genome sequences along with relevant Los Alamos 2010 HIV-1 subtype references (Fig. 1).

Fig. 1

Fig. 1

Back to Top | Article Outline

Subtyping and recombinant inferences

Los Alamos Recombinant Identification Program (RIP) 3.0 was used in our primary analysis.In our gene-by-gene examination, RIP window size 400 with confidence interval (CI) 95% was arbitrary selected to infer subtype for sequences with HXB2 reference sequence more than 600 nucleotides in length (gag, prrt, int, GP41, nef), otherwise, RIP window size 100 with CI 90% was used (vif, vpr, vpu, partial GP120). In our near-full-genome examination, RIP window size 400 with CI 95% was used. A sample was called a ‘nonrecombinant’ when RIP returned a single subtype inference (e.g. A1) and a ‘recombinant’ when RIP detected multiple subtypes within a fragment (e.g. A1-D). We defined a ‘subtype switch’ event as having any discordant subtyping results in a series of longitudinal samples from a single patient (e.g. switching from A1-D recombinant to D). For comparison, we repeated subtyping and recombinants inferences for prrt using another algorithm REGA 2.0 (BIOAFRICA) with default settings.

Back to Top | Article Outline

Definitions of therapy outcomes

Virologic suppression was defined as less than 400 copies HIV RNA/ml to reflect the viral load detection limit during the initial years of the follow-up period. Four therapy outcomes were examined: time to virologic suppression (first of two consecutive viral load <400 copies HIV RNA/ml), time to posttherapy virologic rebound (first of two consecutive viral load ≥400 copies HIV RNA/ml postsuppression, defined as the number of days since the first of two consecutive virologic suppression event), time to first CD4+ rise (any posttherapy CD4+ cell count above baseline), and time to sustained CD4+ recovery (first of two consecutive posttherapy CD4+ cell count increase of more than 200 cells/μl from pretherapy count, or first of two consecutive posttherapy CD4+ cell count >350 cells/μl). ‘Lost to follow-up’ was defined as the lack of an event until study cutoff and having the last clinic visit more than 18 months (548 days) before study cutoff. A sample R script on the definitions and data extraction for virologic outcome is available in Gist (https://gist.github.com/guineverelee/public). Transmitted drug resistance data for this cohort was defined by the WHO list [16] and was previously published by our group [9]. Posttreatment drug resistance was defined using Stanford HIVdb [17].

Back to Top | Article Outline

Statistical methods

All statistical analyses were performed using R, SAS (SAS Institute Inc., Cary, North Carolina, USA) and/or GraphPad Prism 5.0 (GraphPad Software, La Jolla, California, USA). Two-tailed Mann–Whitney tests were used to compare age, pretherapy viral load, pretherapy CD4+ cell count, and follow-up duration. Log-rank tests were used to compare Kaplan–Meier curves of time to virologic suppression, virologic rebound, CD4+ rise and CD4+ recovery. Two-tailed Fisher's exact test was used to compare sex distribution and the extent of lost to follow-up. Multivariate Cox Proportional Hazard confounder models were used for time to virologic suppression and time to postsuppression virologic rebound. The initial list of variables included in the full Cox models were ‘recombinants or nonrecombinants,’ ‘age at enrollment,’ ‘ sex,’ ‘baseline viral load,’ ‘baseline CD4+ cell count,’ ‘number of follow-up visits,’ ‘mean visit interval in days (per 10 increase),’ ‘year of therapy start (dichotomized into >2007 versus ≤2007),’ and ‘type of first regimen (nevirapine-based versus others).’ Then, variables were dropped one-at-a-time using the lowest relative change in the coefficient for the variable related to the outcome as a criterion, until the maximum change from the full model exceeded 5% [18]. Statistical significance was defined as P < 0.05 in all analyses. Note, Bonferroni correction for multiple comparisons could arguably be used instead (P = 0.05/26 = 0.002); this value is provided here for benchmarking purpose.

Back to Top | Article Outline

MiSeq (Illumina) deep sequencing

For each sample, the five nested second-round PCR amplicons from three reactions were pooled (5 μl each) and purified using AMPure Beads (Agilent). Library was prepared with Nextera XT kits (Illumina) to manufacturer's protocol and sequenced with MiSeq Reagent Kit V2 (500 cycles) (Illumina, San Diego, California, USA) with a target coverage depth at 8000. FASTQ outputs were processed by our in-house bioinformatics pipeline (version 6.8). Quality cutoff q15 was chosen. Each genomic region (gag, protease, reverse transcriptase, int, vif, vpr, vpu, env, and nef) was aligned using bowtie separately. Briefly, shortgun reads were initially aligned to its corresponding HXB2 reference sequences followed by a reiterative process to obtain sample-specific reference sequences. Then, one consensus sequence/sample/genomic region was created at a 20% nucleotide mixtures cutoff. Sequences with an average coverage depth less than 10 were excluded from subsequent analyses. ‘Concordance between Sanger and MiSeq’ was defined as ‘either having completely identical subtype inference results (e.g. Sanger A1 versus MiSeq A1) or having partially concordant results (e.g. Sanger A1 versus MiSeq A1-D recombinant).’ The following sequences were excluded from our analyses for quality control: sequences with an average coverage depth less than 10 across the target genomic region, sequences that were shorter than 400 or 100 nucleotides in length for RIP window size 400 and 100 analyses, respectively, and sequences which RIP failed to yield a subtype inference and returned as ‘none significant.’ A secondary analysis specifically examined samples with coverage depth between 10 and 100.

Back to Top | Article Outline

Results

Pretherapy baseline characteristics (n = 504)

Cross-sectional pretherapy baseline HIV+ plasma samples were collected from 2005 to 2010 from 504 HIV-infected patients immediately before initiation of antiretroviral therapy. At baseline, 69% were women, median age was 35 (Q1–Q3 29–39), median baseline viral load and CD4+ cell count were 1 × 105 copies HIV RNA/ml (Q1–Q3 4 × 104–4 × 105) and 132 cells/μl (Q1–Q3 75–200). Initial regimens were primarily nevirapine (NVP) (86%) and efavirenz (EFV) based (12%) in combination with lamivudine (3TC) and zidovudine (ZDV).

Back to Top | Article Outline

Prevalence of HIV-1 intersubtype recombination

When each genomic region was individually examined, intersubtype recombinants were detected at the following frequencies in pretreatment samples: gag (15%), prrt (11%), int (8%), vif (10%), vpr (2%), vpu (9%), GP120 (8%), GP41 (18%), and nef (4%) as shown in Table 1 column gag to nef. Of the 200 patients who had sequence data available across all genomic regions, nucleotides from all genomic regions were concatenated to produce near-full-genome HIV data to be reanalyzed by RIP. Prevalence of intersubtype recombination detected anywhere along the genome was 46% (Table 1, column ‘near-full-genome’). The most frequently detected recombinant was A1-D (25%). Stratification by year revealed no temporal trend in the prevalence of recombinants (all P > 0.1, Supplementary Table 4, http://links.lww.com/QAD/B149). Phylogenetic analysis by maximum-likelihood tree of these 200 near-full-genome sequences showed that most A1-D and D-A1 recombinants did not cluster into monophyletic groups and did not share common ancestor(s), suggesting multiple recombination events (Fig. 1). Recombinant breakpoints scattered across the genome (Fig. 2). To assess the consistency of RIP's subtype inferences compared with other subtyping algorithms, we submitted all prrt sequences to REGA 2.0 (BIOAFRICA) for a representative comparison. Overall, REGA predicted a higher prevalence of recombinants in prrt (20%) than RIP (11%) and their results were 89% concordant.

Table 1

Table 1

Fig. 2

Fig. 2

Back to Top | Article Outline

Longitudinal switches in HIV-1 subtypes

Longitudinal HIV sequence data was available from 143 individuals, but none of them were among the 200 individuals who had near-complete full-genome data. Therefore, we proceeded to estimate prevalence of subtype switches in individual genomic regions. Longitudinal sequences were available from 77 study participants (gag), 90 (prrt), 67 (int), 63 (vif), 59 (vpr), 77 (vpu), 31 (GP120), 80 (GP41), and 86 (nef). Longitudinal subtype switches were observed in 5% of the study participants (gag), 9% (prrt), 6% (int), 2% (vif), 3% (vpr), 3% (vpu), 13% (GP120), 8% (GP41) and 7% (nef). Details of these 37 cases of subtype switches in 27 of 143 study participants (19%) are listed in Supplementary Tables 5a–i, http://links.lww.com/QAD/B149. To explore potential reasons for switches, we examined phylogenetic trees, and found that one discordant sample from patient MBA1100 was an immediate neighbor of MBA1101 in a phylogenetic tree of all available vif sequences, suggesting potential labeling error (Supplementary Table 5b, http://links.lww.com/QAD/B149), whereas in another case nef sequences of patient MBA1435 switched from A1 to G, but all sequences formed a tight monophyletic group on a phylogenetic tree of the entire cohort's nef sequences, suggesting potential RIP artifacts (Supplementary Table 5i, http://links.lww.com/QAD/B149). No potential explanation for other switch cases was found.

Back to Top | Article Outline

Comparison with Illumina (MiSeq) near-full-genome deep sequencing data

To examine whether another sequencing method would yield consistent subtyping results, we performed near-full-genome Nextera XT MiSeq deep sequencing, consensus nucleotide sequence generation, and RIP subtyping on 23 randomly selected samples. Of these, 20 were pretherapy baseline samples from 20 unique individuals, and three were posttherapy samples from three unique individuals. Paired Sanger and MiSeq data were successfully obtained for 22 (gag), 22 (prrt), 11 (int), 19 (vif), 19 (vpr), 22 (vpu), 19 (partial GP120), 16 (GP41), 20 (nef), and 19 (near-full-genome). RIP subtype inferences were always concordant between Sanger sequences and MiSeq-derived consensus in all of gag, prrt, int, vpr, vpu, GP120, GP41, nef, and full-genome data; the two discordant cases were observed in vif (first case Sanger ‘D’ but MiSeq ‘C’; second case Sanger ‘D’ but MiSeq ‘B’). Among the 19 pretherapy samples with near-full-genome MiSeq data, 11 (58%) were inferred by RIP as recombinants, closely resembling the value estimate by Sanger near-full-genome data (46%, Table 1).

Furthermore, among all the above paired results, we observed 30 cases in which the average MiSeq coverage depths were extremely low between 10 and 100, mainly in int (9/30, 30%) and vif (12/30, 40%). Interestingly, in 29/30 cases, RIP subtype inferences at this low coverage depth remained completely concordant between paired Sanger and MiSeq-derived sequences. The remaining 1/30 case was a GP120 sequence; RIP failed to associate both the Sanger and MiSeq-derived sequences with any known subtype or recombinants. Also of note was that of the 19 paired GP120 Sanger and MiSeq sequences that covered the first 160 nucleotides of GP120, MiSeq successfully yielded full-length GP120 consensus sequences that included all the variable loops in 16/19 (84%) cases. In contrast, we were unable to obtain clean Sanger sequence data beyond approximately the 160th nucleotide in GP120 because of the high sequence diversity in env resulting in base-calling ambiguities.

Back to Top | Article Outline

Baseline clinical correlates and therapy outcomes (n = 200)

We dichotomized the 200 patients with near-full-genome data into ‘nonrecombinants’ (54%) and ‘recombinants’ (46%) HIV-1 infections, and compared their pre and posttherapy-initiation clinical correlates. At pretherapy baseline, the two groups were not significantly different in sex (67 versus 70% women, P = 0.8 Fisher two tailed), baseline viral load (Fig. 3a, P = 0.7 Mann–Whitney), and baseline CD4+ cell count (Fig. 3b, P = 0.2 Mann–Whitney). Study participants infected with recombinants were slightly younger (median age 35 versus 34, P = 0.04 Mann–Whitney). At posttherapy, univariate tests showed marginally significant differences in ‘time to virologic suppression’ (Fig. 3c nonadjusted, P = 0.03 log rank), but were not significantly different in ‘time to postsuppression virologic rebound’ (Fig. 3d nonadjusted, P = 0.1 log rank), ‘time to first CD4+ rise above baseline’ (Fig. 3e, P = 0.3 log rank), and ‘time to sustained CD4+ recovery’ (Fig. 3f, P = 0.6 log rank). Neither the ‘proportion of study participants lost to follow-up’ nor ‘duration of follow-up’ were significantly different between groups (all P > 0.2). Next, we further explored the marginal differences observed in virologic outcomes using multivariate Cox Proportional Hazard confounder models: after adjustment, both ‘time to virologic suppression’ and ‘time to posttherapy virologic rebound’ were not significantly different between groups (Fig. 3c adjusted P = 0.3, hazard ratio recombinants/nonrecombinants 0.8, 95% CI 0.6–1.1, controlling for visit interval, year of therapy start, and type of first regimen; Fig. 3d adjusted P = 0.4, hazard ratio 1.6, 95% CI 0.6–4.5, controlling for age, sex, baseline viral load, and visit interval). Finally, we compared the prevalence of transmitted drug resistance and the prevalence of any major drug resistance mutations in recombinants versus nonrecombinant groups and did not find significant difference between groups (P = 0.7 and 1.0 Fisher exact two tail).

Fig. 3

Fig. 3

Back to Top | Article Outline

Discussion

In summary, our near-full-genome sequencing approach revealed a high prevalence of infections by intersubtype HIV-1 recombinants (46%, without time trend) in a rural African community where multiple HIV subtypes cocirculate. We also provided evidence that most of these recombinants arose from independent recombination events, and found evidence of longitudinal subtype switches. Importantly, our study provided evidence that infections with recombinant HIV was neither associated with any negative pre nor posttherapy virologic nor immunologic clinical correlates.

Our reported prevalence of Uganda's HIV-1 recombinants in prrt (11%) fell within the range previously reported (6–19%) [4–6]. Importantly, we showed that the prevalence of intersubtype recombinants is a lot higher under a full-genome context (Table 1, column gag to nef). This finding has several implications. First, it points out that studies which used only part of the genome to estimate recombination prevalence may result in underestimation. Second, it suggests that HIV subtype studies on disease progression and clinical impacts could be biased depending on the genomic region used for subtyping. For instance, numerous studies reported that subtype D infections (by pol) as being more aggressive with a faster disease progression [19–21]. However, in our current report, multiple patients had subtype D in pol, but A1 in GP41. This may also have an implication on subtype-specific vaccine development strategies. Third, our observation that these recombinants did not tend to share common ancestors in a phylogenetic tree suggests that recombination is very frequent in HIV natural biology, and may even contribute to the change of subtype over time that we have observed (on top of other explanations such as undetected sample mixed up and/or superinfections). This knowledge of recombination frequency may help further our understanding about the natural ‘template switch’ frequency during the reverse transcription step in HIV's life cycle.

Our conclusion on the prevalence of recombinants is limited by our choice of subtyping algorithm and its settings. We have arbitrary chosen RIP with a window size 400 and CI of 95% for genes that are more than 600 nucleotides in length; otherwise window size 100 and CI of 90% is chosen. It should be noted that a lower ‘input sequence length to RIP window size’ ratio corresponds to a decreased sensitivity for recombinants detection. As these ratios differ in each genomic region examined in our gene-by-gene analyses (Table 1, column gag to nef), comparison of percentage prevalence of recombinants between genomic regions would not be appropriate. A fair estimation of the prevalence of recombinants (46%) was obtained in our near-full-genome concatenation approach: All sequences were 8625 nucleotides in length, and a constant RIP window size of 400 with 95% CI was used. In addition, we also showed that a different subtyping algorithm, REGA, called more recombinants in prrt than RIP did, which highlighted the variability introduced by different algorithms and settings.

Another technical limitation relates to primer mismatch bias, which could potentially lead to a preferential amplification of one of the amplicons over another in our duplex A or duplex B reactions, and/or preferential amplification of one subtype over another. To address this concern, we reviewed our dataset to compare the number of successful amplifications and the distribution of subtypes across amplicons (Table 1). First, we observed that the number of successful amplifications of each region were comparable, ranging from 387 to 486 successes (excluding GP120, of which Sanger sequence data was uncallable over the variant stretch of the genome, likely because of excessive indels). As all 504 samples were subjected to the same reactions A1/A2, B1/B2 and C, this comparable success rate implies that there was at least a low degree of bias in amplification efficiency between amplicons. Next, Table 1 shows that the percentage distribution of subtypes between genomic regions were also comparable across amplicons (duplex A1 contained gag; A2 vpu; duplex B1 prrt; B2 GP41, nef; reaction C: int, vif, vpr). Across these amplicons, %subtype-A1 ranged from 43.0 to 53.4%, %C ranged from 3.6 to 10.1%, and %D ranged from 30.3 to 39.4%. In other words, we did not observe one subtype being over represented in one amplicon compared with another. Although we cannot rule out the absence/presence of primer amplification bias because of primer mismatches, these observations suggest comparable amplification efficiency across genome and across subtypes.

Our conclusion on the clinical impacts of recombinant HIV-1 is limited by sample size and sampling intervals. We only had 200 patients who had near-full-genome sequencing data available and were thus included in the clinical outcome analyses. A much larger study with more patients will be needed to increase statistical power. Furthermore, these patients were scheduled to receive virologic and immunologic monitoring every 3–6 months, which resulted in the ‘staircase-like’ shapes observed in the Kaplan–Meier curves (instead of smooth curves) in Fig. 3c–f. This factor, compounded with the relatively low number of study participants in both groups (n ≤ 100), might have compromised the statistical comparisons.

Finally, our conclusion that Sanger and MiSeq produced highly concordant subtyping results is limited by the lower PCR/sequencing success rate in int and vif. This was potentially because of our ‘equal volume pooling’ approach (instead of an ‘equal mole pooling’ approach) during the MiSeq library preparation step, which reflects different PCR amplification efficiency between amplicons for this particular set of samples. However, our results also suggest that an average MiSeq coverage depth of 10–100 can still very accurately predict Sanger subtyping results across all genomic regions in these Ugandan nonsubtype B HIV samples. This agrees with previous observations that MiSeq can accurately predict Sanger-derived reverse transcriptase sequences at coverage depth below 100 [22]. Furthermore, we showed that MiSeq was more successful in yielding full-length GP120 sequences compared with Sanger sequencing, supporting the move to deep sequencing in HIV env genetic studies.

In conclusion, this study revealed a high prevalence of HIV-1 infections by intersubtype recombinants in a rural African community where subtype A1 and D cocirculate, but showed that infections by recombinants did not impact pre or posttherapy virologic or immunologic clinical correlates and therapy outcomes. Future studies should continue to monitor intersubtype recombinants in this and other similar communities under a full-genome sequencing context to keep check of their spread, evolution, and clinical impacts.

Back to Top | Article Outline

Acknowledgements

We thank all patients participating in the UARTO cohort whom have made this study possible; we also thank Mr Conan Woods for sequence database management and quality control assistance, and Ms Benita Yip for advice on data analysis.

The work presented here was carried out in collaboration between all authors. The study was conceptualized and designed by G.Q.L. and P.R.H. Plasma samples, baseline and follow-up data were collected and managed by A.R.M., D.R.B., J.N.M., P.W.H., S.H.K., and Y.B. HIV-1 genotyping laboratory work was done by C.L., G.Q.L., and T.M. Results were analyzed by C.J.B., G.Q.L., W.Z., and V.D.L. G.Q.L. wrote the manuscript; Presented in part at: the Annual Canadian Conference on HIV/AIDS Research (CAHR 2015; abstract BS20) and the 8th IAS Conference on HIV Pathogenesis, Treatment, and Prevention, Vancouver, British Columbia, Canada (IAS 2015; abstract MOPEA040).

P.R.H. has received consulting fees from ViiV/Pfizer and Quest; holds stock in Merck, Gilead and Illumina. This work was funded by the Canadian Institutes of Health Research (CIHR), National Institutes of Health (NIH) Centers for AIDS Research (CFAR) Program, National Institute of Mental Health, and National Institute on Alcohol Abuse and Alcoholism (R01 MH054907, P30 AI027763, U01 CA066529, UM1 CA181255, R21 AA014784).

Back to Top | Article Outline

Conflicts of interest

There are no conflicts of interest.

Back to Top | Article Outline

References

1. Robertson DL, Anderson JP, Bradac JA, Carr JK, Foley B, Funkhouser RK, et al. HIV-1 nomenclature proposal. Science 2000; 288:55–56.
2. Kuzembayeva M, Dilley K, Sardo L, Hu WS. Life of psi: how full-length HIV-1 RNAs become packaged genomes in the viral particles. Virology 2014; 454–455:362–370.
3. Lau KA, Wong JJ. Current trends of HIV recombination worldwide. Infect Dis Rep 2013; 5:e4.
4. Yebra G, Ragonnet-Cronin M, Ssemwanga D, Parry CM, Logue CH, Cane PA, et al. Analysis of the history and spread of HIV-1 in Uganda using phylodynamics. J Gen Virol 2015; 96:1890–1898.
5. Ssemwanga D, Ndembi N, Lyagoba F, Bukenya J, Seeley J, Vandepitte J, et al. HIV type 1 subtype distribution, multiple infections, sexual networks, and partnership histories in female sex workers in Kampala, Uganda. AIDS Res Hum Retroviruses 2012; 28:357–365.
6. Gale CV, Yirrell DL, Campbell E, Van der Paal L, Grosskurth H, Kaleebu P. Genotypic variation in the pol gene of HIV type 1 in an antiretroviral treatment-naive population in rural southwestern Uganda. AIDS Res Hum Retroviruses 2006; 22:985–992.
7. Harris ME, Serwadda D, Sewankambo N, Kim B, Kigozi G, Kiwanuka N, et al. Among 46 near full length HIV type 1 genome sequences from Rakai District, Uganda, subtype D and AD recombinants predominate. AIDS Res Hum Retroviruses 2002; 18:1281–1290.
8. Doka NI, Jacob ST, Banura P, Moore CC, Meya D, Mayanja-Kizza H, et al. Enrichment of HIV-1 subtype AD recombinants in a Ugandan cohort of severely septic patients. PLoS One 2012; 7:e48356.
9. Lee GQ, Bangsberg DR, Muzoora C, Boum Y, Oyugi JH, Emenyonu N, et al. Prevalence and virologic consequences of transmitted HIV-1 drug resistance in Uganda. AIDS Res Hum Retroviruses 2014; 30:896–906.
10. Lee GQ. Standard and enhanced genotypic assessment of HIV drug resistance and tropism in subtype-B and nonsubtype-B HIV-1 Infections. 2014
11. Lee GQ, Lachowski C, Cai E, et al. Non-R5-tropic HIV-1 in subtype A and D infections were associated with lower pretherapy CD4 count but not with therapy outcomes. AIDS 2016; 30:1781–1788.
12. Kaida A, Matthews LT, Kanters S, Kabakyenga J, Muzoora C, Mocello AR, et al. Incidence and predictors of pregnancy among a cohort of HIV-positive women initiating antiretroviral therapy in Mbarara, Uganda. PLoS One 2013; 8:e63411.
13. Hunt PW, Cao HL, Muzoora C, Ssewanyana I, Bennett J, Emenyonu N, et al. Impact of CD8+ T-cell activation on CD4+ T-cell recovery and mortality in HIV-infected Ugandans initiating antiretroviral therapy. AIDS 2011; 25:2123–2131.
14. Byakwaga H, Boum Y, Huang Y, Muzoora C, Kembabazi A, Weiser SD, et al. The kynurenine pathway of tryptophan catabolism, CD4+ T-cell recovery, and mortality among HIV-infected ugandans initiating antiretroviral therapy. J Infect Dis 2014; 210:383–391.
15. Woods CK, Brumme CJ, Liu TF, Chui CKS, Chu AL, Wynhoven B, et al. Automating HIV drug resistance genotyping with RECall, a freely accessible sequence analysis tool. J Clin Microbiol 2012; 50:1936–1942.
16. Bennett DE, Camacho RJ, Otelea D, Kuritzkes DR, Fleury H, Kiuchi M, et al. Drug resistance mutations for surveillance of transmitted HIV-1 drug-resistance: 2009 update. PLoS One 2009; 4:e4724.
17. Tang MW, Liu TF, Shafer RW. The HIVdb system for HIV-1 genotypic resistance interpretation. Intervirology 2012; 55:98–101.
18. Lima VD, Geller J, Bangsberg DR, Patterson TL, Daniel M, Kerr T, et al. The effect of adherence on the association between depressive symptoms and mortality among HIV-infected individuals first initiating HAART. AIDS 2007; 21:1175–1183.
19. Eller MA, Opollo MS, Liu M, Redd AD, Eller LA, Kityo C, et al. HIV type 1 disease progression to AIDS and death in a rural Ugandan cohort is primarily dependent on viral load despite variable subtype and T-cell immune activation levels. J Infect Dis 2015; 211:1574–1584.
20. Santoro MM, Perno CF. HIV-1 genetic variability and clinical implications. ISRN Microbiol 2013; 2013:481314.
21. Pant Pai N, Shivkumar S, Cajas JM. Does genetic diversity of HIV-1 non-B subtypes differentially impact disease progression in treatment-naive HIV-1-infected individuals? A systematic review of evidence: 1996–2010. J Acquir Immune Defic Syndr 2012; 59:382–388.
22. Lapointe H, Dong W, Lee GQ, et al. HIV drug resistance testing by high-multiplex ‘wide’ sequencing on the Illumina MiSeq. Antimicrob Agents Chemother 2015; 9:6824–6833.
Keywords:

Africa; clinical outcomes; consequence; deep sequencing; full-genome sequencing; HIV-1; non-B subtypes; recombinants; Uganda; virologic outcomes

Supplemental Digital Content

Back to Top | Article Outline
Copyright © 2017 Wolters Kluwer Health, Inc.