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.