Human immunodeficiency virus type 1 genetic recombination is more frequent than that of Moloney murine leukemia virus despite similar template switching rates

Retroviral recombinants result from template switching between copackaged viral genomes. Here, marker reassortment between coexpressed vectors was measured during single replication cycles, and human immunodeficiency virus type 1 (HIV-1) recombination was observed six- to sevenfold more frequently than murine leukemia virus (MLV) recombination. Template switching was also assayed by using transduction-type vectors in which donor and acceptor template regions were joined covalently. In this situation, where RNA copackaging could not vary, MLV and HIV-1 template switching rates were indistinguishable. These findings argue that MLV's lower intermolecular recombination frequency does not reflect enzymological differences. Instead, these data suggest that recombination rates differ because coexpressed MLV RNAs are less accessible to the recombination machinery than are coexpressed HIV RNAs. This hypothesis provides a plausible explanation for why most gammaretrovirus recombinants, although relatively rare, display evidence of multiple nonselected crossovers. By implying that recombinogenic template switching occurs roughly four times on average during the synthesis of every MLV or HIV-1 DNA, these results suggest that virtually all products of retroviral replication are biochemical recombinants.     

Testing the hypothesis of a recombinant origin of human immunodeficiency virus type 1 subtype E

The human immunodeficiency virus type 1 (HIV-1) epidemic in Southeast Asia has been largely due to the emergence of clade E (HIV-1E). It has been suggested that HIV-1E is derived from a recombinant lineage of subtype A (HIV-1A) and subtype E, with multiple breakpoints along the E genome. We obtained complete genome sequences of clade E viruses from Thailand (93TH057 and 93TH065) and from the Central African Republic (90CF11697 and 90CF4071), increasing the total number of HIV-1E complete genome sequences available to seven. Phylogenetic analysis of complete genomes showed that subtypes A and E are themselves monophyletic, although together they also form a larger monophyletic group. The apparent phylogenetic incongruence at different regions of the genome that was previously taken as evidence of recombination is shown to be not statistically significant. Furthermore, simulations indicate that bootscanning and pairwise distance results, previously used as evidence for recombination, can be misleading, particularly when there are differences in substitution or evolutionary rates across the genomes of different subtypes. Taken jointly, our analyses suggest that there is inadequate support for the hypothesis that subtype E variants are derived from a recombinant lineage. In contrast, many other HIV strains claimed to have a recombinant origin, including viruses for which only a single parental strain was employed for analysis, do indeed satisfy the statistical criteria we propose. Thus, while intersubtype recombinant HIV strains are indeed circulating, the criteria for assigning a recombinant origin to viral structures should include statistical testing of alternative hypotheses to avoid inappropriate assignments that would obscure the true evolutionary properties of these viruses.     

Genetic recombination between human immunodeficiency virus type 1 (HIV-1) and HIV-2, two distinct human lentiviruses

Human immunodeficiency virus type 1 (HIV-1) and HIV-2 are genetically distinct viruses that each can cause AIDS. Approximately 1 million people are infected with both HIV-1 and HIV-2. Additionally, these two viruses use the same receptor and coreceptors and can therefore infect the same target cell populations. To explore potential genetic interactions, we first examined whether RNAs from HIV-1 and HIV-2 can be copackaged into the same virion. We used modified near-full-length viruses that each contained a green fluorescent protein gene (gfp) with a different inactivating mutation. Thus, a functional gfp could be reconstituted via recombination, which was used to detect the copackaging of HIV-1 and HIV-2 RNAs. The GFP-positive (GFP⁺) phenotype was detected in approximately 0.2% of the infection events, which was 35-fold lower than the intrasubtype HIV-1 rates. We isolated and characterized 54 GFP⁺ single-cell clones and determined that all of them contained proviruses with reconstituted gfp. We then mapped the general structures of the recombinant viruses and characterized the recombination junctions by DNA sequencing. We observed several different recombination patterns, including those that had crossovers only in gfp. The most common hybrid genomes had heterologous long terminal repeats. Although infrequent, crossovers in the viral sequences were also identified. Taken together, our study demonstrates that HIV-1 and HIV-2 can recombine, albeit at low frequencies. These observations indicate that multiple factors are likely to restrict the generation of viable hybrid HIV-1 and HIV-2 viruses. However, considering the large coinfected human population and the high viral load in patients, these rare events could provide the basis for the generation of novel human immunodeficiency viruses.     

An exploratory algorithm to identify intra-host recombinant viral sequences

Since recombination leads to the generation of mosaic genomes that violate the assumption of traditional phylogenetic methods that sequence evolution can be accurately described by a single tree, results and conclusions based on phylogenetic analysis of data sets including recombinant sequences can be severely misleading. Many methods are able to adequately detect recombination between diverse sequences, for example between different HIV-1 subtypes. More problematic is the identification of recombinants among closely related sequences such as a viral population within a host. We describe a simple algorithmic procedure that enables detection of intra-host recombinants based on split-decomposition networks and a robust statistical test for recombination. By applying this algorithm to several published HIV-1 data sets we conclude that intra-host recombination was significantly underestimated in previous studies and that up to one-third of the env sequences longitudinally sampled from a given subject can be of recombinant origin. The results show that our procedure can be a valuable exploratory tool for detection of recombinant sequences before phylogenetic analysis, and also suggest that HIV-1 recombination in vivo is far more frequent and significant than previously thought.     

Recombination in AIDS viruses

Recombination contributes to the generation of genetic diversity in human immunodeficiency viruses (HIV) but can only occur between viruses replicating within the same cell. Since individuals have not been found to be simultaneously coinfected with multiple divergent strains of HIV-1 or HIV-2, recombination events have been thought to be restricted to the rather closely related members of the quasispecies that evolves during the course of HIV infection. Here we describe examples of both HIV-1 and HIV-2 genomes that appear to be hybrids of genetically quite divergent viruses. Phylogenetic analyses were used to examine the evolutionary relationships among multiple HIV strains. Evolutionary trees derived from different genomic regions were consistent with respect to most of the viruses investigated. However, some strains of HIV-1 and HIV-2 exhibited significantly discordant branching orders indicative of genetic exchanges during their evolutionary histories. The crossover points of these putative recombination events were mapped by examining the distribution of phylogenetically informative sites supporting alternative tree topologies. A similar example of a recombinant simian immunodeficiency virus identified in West African green monkeys has also been described recently. These results indicate that coinfection with highly divergent viral strains can occur in HIV-infected humans and SIV-infected primates and could lead to the generation of hybrid genomes with significantly altered biological properties. Thus, future characterization of primate lentiviruses should include careful phylogenetic investigation of possible genomic mosaicism.Correspondence to: P.M. Sharp     

Frequent Intra-Subtype Recombination among HIV-1 Circulating in Tanzania

The study estimated the prevalence of HIV-1 intra-subtype recombinant variants among female bar and hotel workers in Tanzania. While intra-subtype recombination occurs in HIV-1, it is generally underestimated. HIV-1 env gp120 V1-C5 quasispecies from 45 subjects were generated by single-genome amplification and sequencing (median (IQR) of 38 (28–50) sequences per subject). Recombination analysis was performed using seven methods implemented within the recombination detection program version 3, RDP3. HIV-1 sequences were considered recombinant if recombination signals were detected by at least three methods with p-values of ≤0.05 after Bonferroni correction for multiple comparisons. HIV-1 in 38 (84%) subjects showed evidence for intra-subtype recombination including 22 with HIV-1 subtype A1, 13 with HIV-1 subtype C, and 3 with HIV-1 subtype D. The distribution of intra-patient recombination breakpoints suggested ongoing recombination and showed selective enrichment of recombinant variants in 23 (60%) subjects. The number of subjects with evidence of intra-subtype recombination increased from 29 (69%) to 36 (82%) over one year of follow-up, although the increase did not reach statistical significance. Adjustment for intra-subtype recombination is important for the analysis of multiplicity of HIV infection. This is the first report of high prevalence of intra-subtype recombination in the HIV/AIDS epidemic in Tanzania, a region where multiple HIV-1 subtypes co-circulate. HIV-1 intra-subtype recombination increases viral diversity and presents additional challenges for HIV-1 vaccine design.     

Mortality among Human Immunodeficiency Virus Type 2-Positive Villagers in Rural Guinea-Bissau Is Correlated with Viral Genotype

We present the results of a 6-year study of 131 human immunodeficiency virus (HIV) type 2 (HIV-2)-infected individuals from a rural population in Guinea-Bissau. Proviral DNA sequences 1.3 kb in length were obtained from each individual and, together with clinical data, including proviral load and CD4 and CD8 levels, were used to assess whether viral genotype influences clinical outcome. With a phylogenetic model, a correlation was found between viral genotype and mortality; this correlation was not due to confounding factors, such as age-specific viral strains or cohabitation of patients. The data provide strong evidence for the involvement of viral genetic factors in determining HIV disease progression in vivo. The pattern of association found suggests that virulence factors are multiple and scattered throughout the HIV-2 genome and can be rapidly gained or lost by the virus through a combination of mutation and recombination. These findings may lead to the identification of viral determinants of HIV disease progression.     

Global Patterns of Recombination across Human Viruses

Viral recombination is a major evolutionary mechanism driving adaptation processes, such as the ability of host-switching. Understanding global patterns of recombination could help to identify underlying mechanisms and to evaluate the potential risks of rapid adaptation. Conventional approaches (e.g., those based on linkage disequilibrium) are computationally demanding or even intractable when sequence alignments include hundreds of sequences, common in viral data sets. We present a comprehensive analysis of recombination across 30 genomic alignments from viruses infecting humans. In order to scale the analysis and avoid the computational limitations of conventional approaches, we apply newly developed topological data analysis methods able to infer recombination rates for large data sets. We show that viruses, such as ZEBOV and MARV, consistently displayed low levels of recombination, whereas high levels of recombination were observed in Sarbecoviruses, HBV, HEV, Rhinovirus A, and HIV. We observe that recombination is more common in positive single-stranded RNA viruses than in negatively single-stranded RNA ones. Interestingly, the comparison across multiple viruses suggests an inverse correlation between genome length and recombination rate. Positional analyses of recombination breakpoints along viral genomes, combined with our approach, detected at least 39 nonuniform patterns of recombination (i.e., cold or hotspots) in 18 viral groups. Among these, noteworthy hotspots are found in MERS-CoV and Sarbecoviruses (at spike, Nucleocapsid and ORF8). In summary, we have developed a fast pipeline to measure recombination that, combined with other approaches, has allowed us to find both common and lineage-specific patterns of recombination among viruses with potential relevance in viral adaptation.     

Human immunodeficiency virus type 1 genomic RNA sequences in the female genital tract and blood: compartmentalization and intrapatient recombination

Investigation of human immunodeficiency virus type 1 (HIV-1) in the genital tract of women is crucial to the development of vaccines and therapies. Previous analyses of HIV-1 in various anatomic sites have documented compartmentalization, with viral sequences from each location that were distinct yet phylogenetically related. Full-length RNA genomes derived from different compartments in the same individual, however, have not yet been studied. Furthermore, although there is evidence that intrapatient recombination may occur frequently, recombinants comprising viruses from different sites within one individual have rarely been documented. We compared full-length HIV-1 RNA sequences in the plasma and female genital tract, focusing on a woman with high HIV-1 RNA loads in each compartment who had been infected heterosexually and then transmitted HIV-1 by the same route. We cloned and sequenced 10 full-length HIV-1 RNA genomes from her genital tract and 10 from her plasma. We also compared viral genomes from the genital tract and plasma of four additional heterosexually infected women, sequencing 164 env and gag clones obtained from the two sites. Four of five women, including the one whose complete viral sequences were determined, displayed compartmentalized HIV-1 genomes. Analyses of full-length, compartmentalized sequences made it possible to document complex intrapatient HIV-1 recombinants that were composed of alternating viral sequences characteristic of each site. These findings demonstrate that the genital tract and blood harbor genetically distinct populations of replicating HIV-1 and provide evidence that recombination between strains from the two compartments contributes to rapid evolution of viral sequence variation in infected individuals.     

Deciphering human immunodeficiency virus type 1 transmission and early envelope diversification by single-genome amplification and sequencing

Accurate identification of the transmitted virus and sequences evolving from it could be instrumental in elucidating the transmission of human immunodeficiency virus type 1 (HIV-1) and in developing vaccines, drugs, or microbicides to prevent infection. Here we describe an experimental approach to analyze HIV-1 env genes as intact genetic units amplified from plasma virion RNA by single-genome amplification (SGA), followed by direct sequencing of uncloned DNA amplicons. We show that this strategy precludes in vitro artifacts caused by Taq-induced nucleotide substitutions and template switching, provides an accurate representation of the env quasispecies in vivo, and has an overall error rate (including nucleotide misincorporation, insertion, and deletion) of less than 8 x 10(-5). Applying this method to the analysis of virus in plasma from 12 Zambian subjects from whom samples were obtained within 3 months of seroconversion, we show that transmitted or early founder viruses can be identified and that molecular pathways and rates of early env diversification can be defined. Specifically, we show that 8 of the 12 subjects were each infected by a single virus, while 4 others acquired more than one virus; that the rate of virus evolution in one subject during an 80-day period spanning seroconversion was 1.7 x 10(-5) substitutions per site per day; and that evidence of strong immunologic selection can be seen in Env and overlapping Rev sequences based on nonrandom accumulation of nonsynonymous mutations. We also compared the results of the SGA approach with those of more-conventional bulk PCR amplification methods performed on the same patient samples and found that the latter is associated with excessive rates of Taq-induced recombination, nucleotide misincorporation, template resampling, and cloning bias. These findings indicate that HIV-1 env genes, other viral genes, and even full-length viral genomes responsible for productive clinical infection can be identified by SGA analysis of plasma virus sampled at intervals typical in large-scale vaccine trials and that pathways of viral diversification and immune escape can be determined accurately.     

Covariation of viral recombination with single nucleotide variants during virus evolution revealed by CoVaMa

Adaptation of viruses to their environments occurs through the acquisition of both novel single-nucleotide variants (SNV) and recombination events including insertions, deletions, and duplications. The co-occurrence of SNVs in individual viral genomes during their evolution has been well-described. However, unlike covariation of SNVs, studying the correlation between recombination events with each other or with SNVs has been hampered by their inherent genetic complexity and a lack of bioinformatic tools. Here, we expanded our previously reported CoVaMa pipeline (v0.1) to measure linkage disequilibrium between recombination events and SNVs within both short-read and long-read sequencing datasets. We demonstrate this approach using long-read nanopore sequencing data acquired from Flock House virus (FHV) serially passaged in vitro. We found SNVs that were either correlated or anti-correlated with large genomic deletions generated by nonhomologous recombination that give rise to Defective-RNAs. We also analyzed NGS data from longitudinal HIV samples derived from a patient undergoing antiretroviral therapy who proceeded to virological failure. We found correlations between insertions in the p6^Gag and mutations in Gag cleavage sites. This report confirms previous findings and provides insights on novel associations between SNVs and specific recombination events within the viral genome and their role in viral evolution.     

Multiple barriers to recombination between divergent HIV-1 variants revealed by a dual-marker recombination assay

Recombination is a major force for generating human immunodeficiency virus type 1 (HIV-1) diversity and produces numerous recombinants circulating in the human population. We previously established a cell-based system using green fluorescent protein gene (gfp) as a reporter to study the mechanisms of HIV-1 recombination. We now report an improved system capable of detecting recombination using authentic viral sequences. Frameshift mutations were introduced into the gag gene so that parental viruses do not express full-length Gag; however, recombination can generate a progeny virus that expresses a functional Gag. We demonstrate that this Gag reconstitution assay can be used to detect recombination between two group M HIV-1 variants of the same or of different subtypes. Using both gfp and gag assays, we found that, similar to group M viruses, group O viruses also recombine frequently. When recombination between a group M virus and a group O virus was examined, we found three distinct barriers for intergroup recombination. First, similar to recombination within group M viruses, intergroup recombination is affected by the identity of the dimerization initiation signal (DIS); variants with the same DIS recombined at a higher rate than those with different DIS. Second, using the gfp recombination assay, we showed that intergroup recombination occurs much less frequently than intragroup recombination, even though the gfp target sequence is identical in all viruses. Finally, Gag reconstitution between variants from different groups is further reduced compared with green fluorescent protein, indicating that sequence divergence interferes with recombination efficiency in the gag gene. Compared with identical sequences, we estimate that recombination rates are reduced by 3-fold and by 10- to 13-fold when the target regions in gag contain 91% and 72-73% sequence identities, respectively. These results show that there are at least three distinct mechanisms preventing exchange of genetic information between divergent HIV-1 variants from different groups.     

Comparison of Major and Minor Viral SNPs Identified through Single Template Sequencing and Pyrosequencing in Acute HIV-1 Infection

Massively parallel sequencing (MPS) technologies, such as 454-pyrosequencing, allow for the identification of variants in sequence populations at lower levels than consensus sequencing and most single-template Sanger sequencing experiments. We sought to determine if the greater depth of population sampling attainable using MPS technology would allow detection of minor variants in HIV founder virus populations very early in infection in instances where Sanger sequencing detects only a single variant. We compared single nucleotide polymorphisms (SNPs) during acute HIV-1 infection from 32 subjects using both single template Sanger and 454-pyrosequencing. Pyrosequences from a median of 2400 viral templates per subject and encompassing 40% of the HIV-1 genome, were compared to a median of five individually amplified near full-length viral genomes sequenced using Sanger technology. There was no difference in the consensus nucleotide sequences over the 3.6kb compared in 84% of the subjects infected with single founders and 33% of subjects infected with multiple founder variants: among the subjects with disagreements, mismatches were found in less than 1% of the sites evaluated (of a total of nearly 117,000 sites across all subjects). The majority of the SNPs observed only in pyrosequences were present at less than 2% of the subject’s viral sequence population. These results demonstrate the utility of the Sanger approach for study of early HIV infection and provide guidance regarding the design, utility and limitations of population sequencing from variable template sources, and emphasize parameters for improving the interpretation of massively parallel sequencing data to address important questions regarding target sequence evolution.     

Determining the frequency and mechanisms of HIV-1 and HIV-2 RNA copackaging by single-virion analysis

HIV-1 and HIV-2 are derived from two distinct primate viruses and share only limited sequence identity. Despite this, HIV-1 and HIV-2 Gag polyproteins can coassemble into the same particle and their genomes can undergo recombination, albeit at an extremely low frequency, implying that HIV-1 and HIV-2 RNA can be copackaged into the same particle. To determine the frequency of HIV-1 and HIV-2 RNA copackaging and to dissect the mechanisms that allow the heterologous RNA copackaging, we directly visualized the RNA content of each particle by using RNA-binding proteins tagged with fluorescent proteins to label the viral genomes. We found that when HIV-1 and HIV-2 RNA are present in viral particles at similar ratios, ～10% of the viral particles encapsidate both HIV-1 and HIV-2 RNAs. Furthermore, heterologous RNA copackaging can be promoted by mutating the 6-nucleotide (6-nt) dimer initiation signal (DIS) to discourage RNA homodimerization or to encourage RNA heterodimerization, indicating that HIV-1 and HIV-2 RNA can heterodimerize prior to packaging using the DIS sequences. We also observed that the coassembly of HIV-1 and HIV-2 Gag proteins is not required for the heterologous RNA copackaging; HIV-1 Gag proteins are capable of mediating HIV-1 and HIV-2 RNA copackaging. These results define the cis- and trans-acting elements required for and affecting the heterologous RNA copackaging, a prerequisite for the generation of chimeric viruses by recombination, and also shed light on the mechanisms of RNA-Gag recognition essential for RNA encapsidation.