Patterns of recombination in HIV-1M are influenced by selection disfavouring the survival of recombinants with disrupted genomic RNA and protein structures

A Trait-Based Approach to Predicting Viral Host-Range Evolvability

Predicting the evolution of virus host range has proven to be extremely difficult, in part because of the sheer diversity of viruses, each with unique biology and ecological interactions. We have not solved this problem, but to make the problem more tractable, we narrowed our focus to three traits intrinsic to all viruses that may play a role in host-range evolvability: mutation rate, recombination rate, and phenotypic heterogeneity. Although each trait should increase evolvability, they cannot do so unbounded because fitness trade-offs limit the ability of all three traits to maximize evolvability. By examining these constraints, we can begin to identify groups of viruses with suites of traits that make them especially concerning, as well as ecological and environmental conditions that might push evolution toward accelerating host-range expansion.

Elucidation of Early Evolution of HIV-1 Group M in the Congo Basin Using Computational Methods

The Congo Basin region is believed to be the site of the cross-species transmission event that yielded HIV-1 group M (HIV-1M). It is thus likely that the virus has been present and evolving in the region since that cross-species transmission. As HIV-1M was only discovered in the early 1980s, our directly observed record of the epidemic is largely limited to the past four decades. Nevertheless, by exploiting the genetic relatedness of contemporary HIV-1M sequences, phylogenetic methods provide a powerful framework for investigating simultaneously the evolutionary and epidemiologic history of the virus. Such an approach has been taken to find that the currently classified HIV-1 M subtypes and Circulating Recombinant Forms (CRFs) do not give a complete view of HIV-1 diversity. In addition, the currently identified major HIV-1M subtypes were likely genetically predisposed to becoming a major component of the present epidemic, even before the events that resulted in the global epidemic. Further efforts have identified statistically significant hot- and cold-spots of HIV-1M subtypes sequence inheritance in genomic regions of recombinant forms. In this review we provide ours and others recent findings on the emergence and spread of HIV-1M variants in the region, which have provided insights into the early evolution of this virus.

Pervasive and non-random recombination in near full-length HIV genomes from Uganda

Recombination is an important feature of HIV evolution, occurring both within and between the major branches of diversity (subtypes). The Ugandan epidemic is primarily composed of two subtypes, A1 and D, that have been co-circulating for 50 years, frequently recombining in dually infected patients. Here, we investigate the frequency of recombinants in this population and the location of breakpoints along the genome. As part of the PANGEA-HIV consortium, 1,472 consensus genome sequences over 5 kb have been obtained from 1,857 samples collected by the MRC/UVRI & LSHTM Research unit in Uganda, 465 (31.6 per cent) of which were near full-length sequences (>8 kb). Using the subtyping tool SCUEAL, we find that of the near full-length dataset, 233 (50.1 per cent) genomes contained only one subtype, 30.8 per cent A1 (n = 143), 17.6 per cent D (n = 82), and 1.7 per cent C (n = 8), while 49.9 per cent (n = 232) contained more than one subtype (including A1/D (n = 164), A1/C (n = 13), C/D (n = 9); A1/C/D (n = 13), and 33 complex types). K-means clustering of the recombinant A1/D genomes revealed a section of envelope (C2gp120-TMgp41) is often inherited intact, whilst a generalized linear model was used to demonstrate significantly fewer breakpoints in the gag-pol and envelope C2-TM regions compared with accessory gene regions. Despite similar recombination patterns in many recombinants, no clearly supported circulating recombinant form (CRF) was found, there was limited evidence of the transmission of breakpoints, and the vast majority (153/164; 93 per cent) of the A1/D recombinants appear to be unique recombinant forms. Thus, recombination is pervasive with clear biases in breakpoint location, but CRFs are not a significant feature, characteristic of a complex, and diverse epidemic.

Notes on recombination and reassortment in multipartite/segmented viruses

Besides evolving through nucleotide substitution, viruses frequently also evolve by genetic recombination which can occur when related viral variants co-infect the same cells. Viruses with segmented or multipartite genomes can additionally evolve via the reassortment of genomic components. Various computational techniques are now available for identifying and characterizing recombination and reassortment. While these techniques have revealed both that all well studied segmented and multipartite virus species show some capacity for reassortment, and that recombination is common in many multipartite species, they have indicated that recombination is either rare or does not occur in species with segmented genomes. Reassortment and recombination can make it very difficult to study segmented/multipartite viruses using metagenomics-based approaches. Notable challenges include, both the accurate identification and assignment of genomic components to individual genomes, and the differentiation between natural 'real' recombination events and artifactual 'fake' recombination events arising from the inaccurate de novo assembly of genome component sequences determined using short read sequencing.

Patterns of genomic site inheritance in HIV-1M inter-subtype recombinants delineate the most likely genomic sites of subtype-specific adaptation

Recombination between different HIV-1 group M (HIV-1M) subtypes is a major contributor to the ongoing genetic diversification of HIV-1M. However, it remains unclear whether the different genome regions of recombinants are randomly inherited from the different subtypes. To elucidate this, we analysed the distribution within 82 circulating and 201 unique recombinant forms (CRFs/URFs), of genome fragments derived from HIV-1M Subtypes A, B, C, D, F, and G and CRF01_AE. We found that viruses belonging to the analysed HIV-1M subtypes and CRF01_AE contributed certain genome fragments more frequently during recombination than other fragments. Furthermore, we identified statistically significant hot-spots of Subtype A sequence inheritance in genomic regions encoding portions of Gag and Nef, Subtype B in Pol, Tat and Env, Subtype C in Vif, Subtype D in Pol and Env, Subtype F in Gag, Subtype G in Vpu-Env and Nef, and CRF01_AE inheritance in Vpu and Env. The apparent non-randomness in the frequencies with which different subtypes have contributed specific genome regions to known HIV-1M recombinants is consistent with selection strongly impacting the survival of inter-subtype recombinants. We propose that hotspots of genomic region inheritance are likely to demarcate the locations of subtype-specific adaptive genetic variations.

An in vitro Model to Mimic Selection of Replication-Competent HIV-1 Intersubtype Recombination in Dual or Superinfected Patients

The low frequency of HIV-1 recombinants within entire viral populations in both individual patients and culture-based infection models impedes investigation of the underlying factors contributing to either the occurrence of recombinants or the survival of recombinants once they are formed. So far, most of the related studies have no consideration of recombinants' functionality. Here, we established a functional recombinant production (FRP) system to produce pure and functional HIV-1 intersubtype Env recombinants and utilized 454 pyrosequencing to investigate the distribution of over 4000 functional and non-functional recombination breakpoints from either the FRP system or dual infection cultures. The results revealed that most of the breakpoints converged in gp41 (62%) and C1 (25.3%) domains of gp120, which has strong correlation with the similarity between the two recombining sequences. Yet, the breakpoints also appeared in C2 (5.2%) and C5 (4.6%) domains not correlated with the recombining sequence similarity. Interestingly, none of the intersubtype gp120 recombinants recombined between C1 and gp41 regions either from the FRP system or from the dual infection culture, and very few from the HIV epidemic were functional. The present study suggests that the selection of functional Env recombinants is one of the reasons for the predominance of C1 and gp41 Env recombinants in the HIV epidemic, and it provides an in vitro model to mimic the selection of replication-competent HIV-1 intersubtype recombination in dual or superinfected patients.

Paucity of Intact Non-Induced Provirus with Early, Long-Term Antiretroviral Therapy of Perinatal HIV Infection

The latent reservoir is a major barrier to HIV eradication. Reservoir size is emerging as an important biomarker to assess the likelihood of HIV remission in the absence of antiretroviral therapy (ART) and may be reduced by earlier initiation of ART that restricts HIV spread into CD4+ T cells. Reservoir size is traditionally measured with a quantitative viral outgrowth assay (QVOA) that induces replication-competent HIV production through in vitro stimulation of resting CD4+ T cells. However, the recent identification of replication-intact, non-induced proviral genomes (NIPG) suggests the QVOA significantly underestimates (by 62-fold) latent reservoir size in chronically-infected adults. Whether formation and persistence of Intact, NIPG is thwarted by early ART initiation and long-term virologic suppression in perinatal infection is unclear. Here, we show that the latent reservoir in 11 early treated, long-term suppressed perinatally infected children and adolescents was not inducible by QVOA and dominated by defective, NIPG. Single genome analysis of 164 NIPG from 232 million cultured resting CD4+ T cells revealed no replication-intact, near-full length sequences. Forty-three (26%) NIPG contained APOBEC3G-mediated hypermutation, 115 (70%) NIPG contained large internal deletions, one NIPG contained nonsense mutations and indels, and 5 (3%) NIPG were assigned as “Not Evaluable” due to multiple failed sequencing attempts that precluded further classification. The lack of replication competent inducible provirus and intact NIPG in this cohort indicate early, long-term ART of perinatal infection leads to marked diminution of replication-competent HIV-1 reservoirs, creating a favorable state towards interventions aimed at virologic remission.

RDP4: Detection and analysis of recombination patterns in virus genomes

RDP4 is the latest version of recombination detection program (RDP), a Windows computer program that implements an extensive array of methods for detecting and visualising recombination in, and stripping evidence of recombination from, virus genome sequence alignments. RDP4 is capable of analysing twice as many sequences (up to 2,500) that are up to three times longer (up to 10 Mb) than those that could be analysed by older versions of the program. RDP4 is therefore also applicable to the analysis of bacterial full-genome sequence datasets. Other novelties in RDP4 include (1) the capacity to differentiate between recombination and genome segment reassortment, (2) the estimation of recombination breakpoint confidence intervals, (3) a variety of ‘recombination aware’ phylogenetic tree construction and comparison tools, (4) new matrix-based visualisation tools for examining both individual recombination events and the overall phylogenetic impacts of multiple recombination events and (5) new tests to detect the influences of gene arrangements, encoded protein structure, nucleic acid secondary structure, nucleotide composition, and nucleotide diversity on recombination breakpoint patterns. The key feature of RDP4 that differentiates it from other recombination detection tools is its flexibility. It can be run either in fully automated mode from the command line interface or with a graphically rich user interface that enables detailed exploration of both individual recombination events and overall recombination patterns.