The distribution of HIV-1 recombination breakpoints

Identification of the Near Full-length Genome of a Novel HIV-1 CRF01_AE/CRF07_BC Recombinant with a Complex Genomic Structure Isolated in Hebei Province, China

BACKGROUND

During HIV genotypic drug resistance testing of patient samples in Baoding, Hebei Province, China, in 2022, a recombinant fragment was detected in the pol region of an HIV-1 strain.

OBJECTIVE

The objective of the study was to analyze the near full-length genome of a novel HIV-1 CRF01_AE/CRF07_BC recombinant with a complex genomic structure.

METHODS

Viral RNA was extracted from the blood of the infected individual and reverse transcribed to cDNA. Two overlapping segments of the HIV-1 genome were amplified using a nearendpoint dilution method and sequenced. Recombinant breakpoints were determined using RIP, jpHMM, and SimPlot 3.5.1 software. MEGA 6.0 software was used to construct a neighbor-joining phylogenetic tree.

RESULTS

We obtained the near full-length genome sequence (8680 bp) of a novel HIV-1 CRF01_AE/CRF07_BC recombinant. Recombination analysis showed that the genome comprised at least 12 overlapping segments, including six CRF07_BC and six CRF01_AE segments, with CRF07_BC as the backbone. The emergence of CRF01_AE/CRF07_BC recombinant strains indicated that HIV-1 co-infection is common. However, the increasing genetic complexity of the HIV-1 epidemic in China warrants continued investigation.

CONCLUSION

The increase in CRF01_AE/CRF07_BC recombinant viruses suggests that HIV-1 has a high genetic mutation rate in Hebei, China. This highlights the need for close monitoring of HIV-1 molecular epidemiologic changes to provide accurate, up-to-date information for effective disease control.

WGS- versus ORF5-Based Typing of PRRSV: A Belgian Case Study

Porcine reproductive and respiratory syndrome virus (PRRSV) is the causative agent of one of the most widespread and economically devastating diseases in the swine industry. Typing circulating PRRSV strains by means of sequencing is crucial for developing adequate control strategies. Most genetic studies only target the highly variable open reading frame (ORF) 5, for which an extensive database is available. In this study, we performed whole-genome sequencing (WGS) on a collection of 124 PRRSV-1 positive serum samples that were collected over a 5-year period (2015–2019) in Belgium. Our results show that (nearly) complete PRRSV genomes can be obtained directly from serum samples with a high success rate. Analysis of the coding regions confirmed the exceptionally high genetic diversity, even among Belgian PRRSV-1 strains. To gain more insight into the added value of WGS, we performed phylogenetic cluster analyses on separate ORF datasets as well as on a single, concatenated dataset (CDS) containing all ORFs. A comparison between the CDS and ORF clustering schemes revealed numerous discrepancies. To explain these differences, we performed a large-scale recombination analysis, which allowed us to identify a large number of potential recombination events that were scattered across the genome. As PRRSV does not contain typical recombination hot-spots, typing PRRSV strains based on a single ORF is not recommended. Although the typing accuracy can be improved by including multiple regions, our results show that the full genetic diversity among PRRSV strains can only be captured by analysing (nearly) complete genomes. Finally, we also identified several vaccine-derived recombinant strains, which once more raises the question of the safety of these vaccines.

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.

Characterization of a Novel HIV-1 CRF01_AE/CRF07_BC Recombinant Strain Among Men Who Have Sex with Men in Liaoning, China

CRF01_AE and CRF07_BC are two widespread human immunodeficiency virus type 1 (HIV-1) strains among different high-risk populations, including men who have sex with men (MSM), in China. This co-epidemic of various HIV strains enables the production of second-generation recombinants. In this study, we detected a novel HIV-1 CRF01_AE/CRF07_BC recombinant from LN321945, an MSM lived in Liaoning province, northeast China. The phylogenetic and recombination analyses indicated the near full-length genome (NFLG) sequence of LN321945 had six recombination breakpoints, with three CRF07_BC fragments inserted into a CRF01_AE backbone. Further subregion trees analysis revealed that both CRF01_AE and CRF07_BC fragments were derived from two predominant HIV-1 strains among MSM. In addition, the NFLG of LN321945 was revealed to be clustered closely to another CRF01_AE/CRF07_BC recombinant previously identified in Shaanxi province, northwest China, but these two recombinants had distinct recombination structure and origin of CRF01_AE fragments. Hence, this study identified a second-generation recombinant between the main strains circulating among MSM, indicating more complicated trend of HIV-1 epidemic in China.

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.

Near full genome characterization of HIV-1 unique recombinant forms in Cameroon reveals dominant CRF02_AG and F2 recombination patterns

INTRODUCTION

In Cameroon, a manifold diversity of HIV strains exists with CRF02_AG and unique recombinant forms (URFs) being the predominant strains. In recent years, a steady increase in URFs and clade F2 viruses has been monitored through partial genome sequencing. There is an information gap in the characterization of emerging URFs along the full genome, which is needed to address the challenges URFs pose towards diagnosis, treatment and HIV-1 vaccine design.

METHOD

Eighteen Cameroonian URFs from samples collected between the years 2000 and 2015 were studied using a newly developed near full genome sequencing (NFGS) protocol based on variable nested RT-PCRs with a versatile primer set. Near full genomes were characterized for recombination patterns and sequence signatures with possible impact on antiretroviral treatment or Env-directed immune responses. Third-generation sequencing (3GS) of near full or half genomes (HGs) gave insight into intra-patient URF diversity.

RESULTS

The characterized URFs were composed of a broad variety of subtypes and recombinants including A, F, G, CRF01_AE, CRF02_AG and CRF22_01A1. Phylogenetic analysis unveiled dominant CRF02_AG and F2 recombination patterns. 3GS indicated a high intra-patient URF diversity with up to four distinct viral sub-populations present in plasma at the same time. URF pol genomic analysis revealed a number of accessory drug resistance mutations (DRMs) in the ART-naïve participants. Genotypic env analysis suggests CCR5 usage in 14/18 samples and identified deviations at residues, critical for gp120/gp41 interphase and CD4 binding site broadly neutralizing antibodies in more than half of the studied URFs. V1V2 sites of immune pressure in the human RV144 vaccine study varied in more than a third of URFs.

CONCLUSIONS

This study identified novel mosaic patterns in URFs in Cameroon. In line with the regional predominance of CRF_02AG and the increased prevalence of clade F2, prominent CRF_02AG and F2 background patterns were observed underlying the URFs. In the context of the novel mosaic genomes, the impact of the identified accessory DRMs and Env epitope variations on treatment and immune control remains elusive. The evolving diversity of HIV-1 URFs in Cameroon requires continuous monitoring to respond to the increasing challenges for diagnosis, antiretroviral treatment and prevention.

Development of a Versatile, Near Full Genome Amplification and Sequencing Approach for a Broad Variety of HIV-1 Group M Variants

Near full genome sequencing (NFGS) of HIV-1 is required to assess the genetic composition of HIV-1 strains comprehensively. Population-wide, it enables a determination of the heterogeneity of HIV-1 and the emergence of novel/recombinant strains, while for each individual it constitutes a diagnostic instrument to assist targeted therapeutic measures against viral components. There is still a lack of robust and adaptable techniques for efficient NFGS from miscellaneous HIV-1 subtypes. Using rational primer design, a broad primer set was developed for the amplification and sequencing of diverse HIV-1 group M variants from plasma. Using pure subtypes as well as diverse, unique recombinant forms (URF), variable amplicon approaches were developed for NFGS comprising all functional genes. Twenty-three different genomes composed of subtypes A (A1), B, F (F2), G, CRF01_AE, CRF02_AG, and CRF22_01A1 were successfully determined. The NFGS approach was robust irrespective of viral loads (≥306 copies/mL) and amplification method. Third-generation sequencing (TGS), single genome amplification (SGA), cloning, and bulk sequencing yielded similar outcomes concerning subtype composition and recombinant breakpoint patterns. The introduction of a simple and versatile near full genome amplification, sequencing, and cloning method enables broad application in phylogenetic studies of diverse HIV-1 subtypes and can contribute to personalized HIV therapy and diagnosis.

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.

The diversification of begomovirus populations is predominantly driven by mutational dynamics

Begomoviruses (single-stranded DNA plant viruses) are responsible for serious agricultural threats. Begomovirus populations exhibit a high degree of within-host genetic variation and evolve as quickly as RNA viruses. Although the recombination-prone nature of begomoviruses has been extensively demonstrated, the relative contribution of recombination and mutation to the genetic variation of begomovirus populations has not been assessed. We estimated the genetic variability of begomovirus datasets from around the world. An uneven distribution of genetic variation across the length of the cp and rep genes due to recombination was evident from our analyses. To estimate the relative contributions of recombination and mutation to the genetic variability of begomoviruses, we mapped all substitutions over maximum likelihood trees and counted the number of substitutions on branches which were associated with recombination (η_r) and mutation (η_μ). In addition, we also estimated the per generation relative rates of both evolutionary mechanisms (r/μ) to express how frequently begomovirus genomes are affected by recombination relative to mutation. We observed that the composition of genetic variation in all begomovirus datasets was dominated by mutation. Additionally, the low correlation between the estimates indicated that the relative contributions of recombination and mutation are not necessarily a function of their relative rates. Our results show that, although a considerable fraction of the genetic variation levels could be assigned to recombination, it was always lower than that due to mutation, indicating that the diversification of begomovirus populations is predominantly driven by mutational dynamics.

Viral reverse transcriptases

Reverse transcriptases (RTs) play a major role in the replication of Retroviridae, Metaviridae, Pseudoviridae, Hepadnaviridae and Caulimoviridae. RTs are enzymes that are able to synthesize DNA using RNA or DNA as templates (DNA polymerase activity), and degrade RNA when forming RNA/DNA hybrids (ribonuclease H activity). In retroviruses and LTR retrotransposons (Metaviridae and Pseudoviridae), the coordinated action of both enzymatic activities converts single-stranded RNA into a double-stranded DNA that is flanked by identical sequences known as long terminal repeats (LTRs). RTs of retroviruses and LTR retrotransposons are active as monomers (e.g. murine leukemia virus RT), homodimers (e.g. Ty3 RT) or heterodimers (e.g. human immunodeficiency virus type 1 (HIV-1) RT). RTs lack proofreading activity and display high intrinsic error rates. Besides, high recombination rates observed in retroviruses are promoted by poor processivity that causes template switching, a hallmark of reverse transcription. HIV-1 RT inhibitors acting on its polymerase activity constitute the backbone of current antiretroviral therapies, although novel drugs, including ribonuclease H inhibitors, are still necessary to fight HIV infections. In Hepadnaviridae and Caulimoviridae, reverse transcription leads to the formation of nicked circular DNAs that will be converted into episomal DNA in the host cell nucleus. Structural and biochemical information on their polymerases is limited, although several drugs inhibiting HIV-1 RT are known to be effective against the human hepatitis B virus polymerase. In this review, we summarize current knowledge on reverse transcription in the five virus families and discuss available biochemical and structural information on RTs, including their biosynthesis, enzymatic activities, and potential inhibition.

Analysis of HIV-1 intersubtype recombination breakpoints suggests region with high pairing probability may be a more fundamental factor than sequence similarity affecting HIV-1 recombination

Background

With increasing data on HIV-1, a more relevant molecular model describing mechanism details of HIV-1 genetic recombination usually requires upgrades. Currently an incomplete structural understanding of the copy choice mechanism along with several other issues in the field that lack elucidation led us to perform an analysis of the correlation between breakpoint distributions and (1) the probability of base pairing, and (2) intersubtype genetic similarity to further explore structural mechanisms.

Methods

Near full length sequences of URFs from Asia, Europe, and Africa (one sequence/patient), and representative sequences of worldwide CRFs were retrieved from the Los Alamos HIV database. Their recombination patterns were analyzed by jpHMM in detail. Then the relationships between breakpoint distributions and (1) the probability of base pairing, and (2) intersubtype genetic similarities were investigated.

Results

Pearson correlation test showed that all URF groups and the CRF group exhibit the same breakpoint distribution pattern. Additionally, the Wilcoxon two-sample test indicated a significant and inexplicable limitation of recombination in regions with high pairing probability. These regions have been found to be strongly conserved across distinct biological states (i.e., strong intersubtype similarity), and genetic similarity has been determined to be a very important factor promoting recombination. Thus, the results revealed an unexpected disagreement between intersubtype similarity and breakpoint distribution, which were further confirmed by genetic similarity analysis. Our analysis reveals a critical conflict between results from natural HIV-1 isolates and those from HIV-1-based assay vectors in which genetic similarity has been shown to be a very critical factor promoting recombination.

Conclusions

These results indicate the region with high-pairing probabilities may be a more fundamental factor affecting HIV-1 recombination than sequence similarity in natural HIV-1 infections. Our findings will be relevant in furthering the understanding of HIV-1 recombination mechanisms.

Electronic supplementary material

The online version of this article (doi:10.1186/s12985-016-0616-1) contains supplementary material, which is available to authorized users.

Coalescent Inference Using Serially Sampled, High-Throughput Sequencing Data from Intrahost HIV Infection

Human immunodeficiency virus (HIV) is a rapidly evolving pathogen that causes chronic infections, so genetic diversity within a single infection can be very high. High-throughput "deep" sequencing can now measure this diversity in unprecedented detail, particularly since it can be performed at different time points during an infection, and this offers a potentially powerful way to infer the evolutionary dynamics of the intrahost viral population. However, population genomic inference from HIV sequence data is challenging because of high rates of mutation and recombination, rapid demographic changes, and ongoing selective pressures. In this article we develop a new method for inference using HIV deep sequencing data, using an approach based on importance sampling of ancestral recombination graphs under a multilocus coalescent model. The approach further extends recent progress in the approximation of so-called conditional sampling distributions, a quantity of key interest when approximating coalescent likelihoods. The chief novelties of our method are that it is able to infer rates of recombination and mutation, as well as the effective population size, while handling sampling over different time points and missing data without extra computational difficulty. We apply our method to a data set of HIV-1, in which several hundred sequences were obtained from an infected individual at seven time points over 2 years. We find mutation rate and effective population size estimates to be comparable to those produced by the software BEAST. Additionally, our method is able to produce local recombination rate estimates. The software underlying our method, Coalescenator, is freely available.

A model-based information sharing protocol for profile Hidden Markov Models used for HIV-1 recombination detection

BACKGROUND

In many applications, a family of nucleotide or protein sequences classified into several subfamilies has to be modeled. Profile Hidden Markov Models (pHMMs) are widely used for this task, modeling each subfamily separately by one pHMM. However, a major drawback of this approach is the difficulty of dealing with subfamilies composed of very few sequences. One of the most crucial bioinformatical tasks affected by the problem of small-size subfamilies is the subtyping of human immunodeficiency virus type 1 (HIV-1) sequences, i.e., HIV-1 subtypes for which only a small number of sequences is known.

RESULTS

To deal with small samples for particular subfamilies of HIV-1, we introduce a novel model-based information sharing protocol. It estimates the emission probabilities of the pHMM modeling a particular subfamily not only based on the nucleotide frequencies of the respective subfamily but also incorporating the nucleotide frequencies of all available subfamilies. To this end, the underlying probabilistic model mimics the pattern of commonality and variation between the subtypes with regards to the biological characteristics of HI viruses. In order to implement the proposed protocol, we make use of an existing HMM architecture and its associated inference engine.

CONCLUSIONS

We apply the modified algorithm to classify HIV-1 sequence data in the form of partial HIV-1 sequences and semi-artificial recombinants. Thereby, we demonstrate that the performance of pHMMs can be significantly improved by the proposed technique. Moreover, we show that our algorithm performs significantly better than Simplot and Bootscanning.

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

Genetic recombination is a major contributor to the ongoing diversification of HIV. It is clearly apparent that across the HIV-genome there are defined recombination hot and cold spots which tend to co-localise both with genomic secondary structures and with either inter-gene boundaries or intra-gene domain boundaries. There is also good evidence that most recombination breakpoints that are detectable within the genes of natural HIV recombinants are likely to be minimally disruptive of intra-protein amino acid contacts and that these breakpoints should therefore have little impact on protein folding. Here we further investigate the impact on patterns of genetic recombination in HIV of selection favouring the maintenance of functional RNA and protein structures. We confirm that chimaeric Gag p24, reverse transcriptase, integrase, gp120 and Nef proteins that are expressed by natural HIV-1 recombinants have significantly lower degrees of predicted folding disruption than randomly generated recombinants. Similarly, we use a novel single-stranded RNA folding disruption test to show that there is significant, albeit weak, evidence that natural HIV recombinants tend to have genomic secondary structures that more closely resemble parental structures than do randomly generated recombinants. These results are consistent with the hypothesis that natural selection has acted both in the short term to purge recombinants with disrupted RNA and protein folds, and in the longer term to modify the genome architecture of HIV to ensure that recombination prone sites correspond with those where recombination will be minimally deleterious.

U3 region in the HIV-1 genome adopts a G-quadruplex structure in its RNA and DNA sequence

Genomic regions rich in G residues are prone to adopt G-quadruplex structure. Multiple Sp1-binding motifs arranged in tandem have been suggested to form this structure in promoters of cancer-related genes. Here, we demonstrate that the G-rich proviral DNA sequence of the HIV-1 U3 region, which serves as a promoter of viral transcription, adopts a G-quadruplex structure. The sequence contains three binding elements for transcription factor Sp1, which is involved in the regulation of HIV-1 latency, reactivation, and high-level virus expression. We show that the three Sp1 binding motifs can adopt different forms of G-quadruplex structure and that the Sp1 protein can recognize and bind to its site folded into a G-quadruplex. In addition, a c-kit2 specific antibody, designated hf2, binds to two different G-quadruplexes formed in Sp1 sites. Since U3 is encoded at both viral genomic ends, the G-rich sequence is also present in the RNA genome. We demonstrate that the RNA sequence of U3 forms dimers with characteristics known for intermolecular G-quadruplexes. Together with previous reports showing G-quadruplex dimers in the gag and cPPT regions, these results suggest that integrity of the two viral genomes is maintained through numerous intermolecular G-quadruplexes formed in different RNA genome locations. Reconstituted reverse transcription shows that the potassium-dependent structure formed in U3 RNA facilitates RT template switching, suggesting that the G-quadruplex contributes to recombination in U3.

Mechanism of HIV-1 RNA dimerization in the central region of the genome and significance for viral evolution

The genome of HIV-1 consists of two identical or nearly identical RNA molecules. The RNA genomes are held in the same, parallel orientation by interactions at the dimer initiation site (DIS). Previous studies showed that in addition to interactions at DIS, sequences located 100 nucleotides downstream from the 5' splice site can dimerize in vitro through an intermolecular G-quartet structure. Here we report that the highly conserved G-rich sequence in the middle portion of the HIV-1 genome near the central polypurine tract (cPPT) dimerizes spontaneously under high ionic strength in the absence of protein. The antisense RNA does not dimerize, strongly indicating that RNA dimerization does not exclusively involve A:U and G:C base pairing. The cation-dependent reverse transcriptase pausing profile, CD spectra profile, and cation-dependent association and thermal dissociation characteristics indicate G-quartet structures. Different forms of G-quartets are formed including monomers and, significantly, intermolecular dimers. Our results indicate that RNA genome dimerization and parallel alignment initiated through interactions at DIS may be greatly expanded and stabilized by formation of an intermolecular G-quartet at a distant site near the cPPT. It is likely that formation of G-quartet structure near the cPPT in vivo keeps the RNA genomes in proximity over a long range, promoting genetic recombination in numerous hot spots.

Molecular evolution of the HIV-1 Thai epidemic between the time of RV144 immunogen selection to the execution of the vaccine efficacy trial

The RV144 HIV-1 vaccine trial (Thailand, 2003 to 2009), using immunogens genetically matched to the regional epidemic, demonstrated the first evidence of efficacy for an HIV-1 vaccine. Here we studied the molecular evolution of the HIV-1 epidemic from the time of immunogen selection to the execution of the efficacy trial. We studied HIV-1 genetic diversity among 390 volunteers who were deferred from enrollment in RV144 due to preexisting HIV-1 infection using a multiregion hybridization assay, full-genome sequencing, and phylogenetic analyses. The subtype distribution was 91.7% CRF01_AE, 3.5% subtype B, 4.3% B/CRF01_AE recombinants, and 0.5% dual infections. CRF01_AE strains were 31% more diverse than the ones from the 1990s Thai epidemic. Sixty-nine percent of subtype B strains clustered with the cosmopolitan Western B strains. Ninety-three percent of B/CRF01_AE recombinants were unique; recombination breakpoint analysis showed that these strains were highly embedded within the larger network that integrates recombinants from East/Southeast Asia. Compared to Thai sequences from the early 1990s, the distance to the RV144 immunogens increased 52% to 68% for CRF01_AE Env immunogens and 12% to 29% for subtype B immunogens. Forty-three percent to 48% of CRF01_AE sequences differed from the sequence of the vaccine insert in Env variable region 2 positions 169 and 181, which were implicated in vaccine sieve effects in RV144. In conclusion, compared to the molecular picture at the early stages of vaccine development, our results show an overall increase in the genetic complexity of viruses in the Thai epidemic and in the distance to vaccine immunogens, which should be considered at the time of the analysis of the trial results.

Stochastic simulations suggest that HIV-1 survives close to its error threshold

The use of mutagenic drugs to drive HIV-1 past its error threshold presents a novel intervention strategy, as suggested by the quasispecies theory, that may be less susceptible to failure via viral mutation-induced emergence of drug resistance than current strategies. The error threshold of HIV-1, , however, is not known. Application of the quasispecies theory to determine poses significant challenges: Whereas the quasispecies theory considers the asexual reproduction of an infinitely large population of haploid individuals, HIV-1 is diploid, undergoes recombination, and is estimated to have a small effective population size in vivo. We performed population genetics-based stochastic simulations of the within-host evolution of HIV-1 and estimated the structure of the HIV-1 quasispecies and . We found that with small mutation rates, the quasispecies was dominated by genomes with few mutations. Upon increasing the mutation rate, a sharp error catastrophe occurred where the quasispecies became delocalized in sequence space. Using parameter values that quantitatively captured data of viral diversification in HIV-1 patients, we estimated to be substitutions/site/replication, ∼2–6 fold higher than the natural mutation rate of HIV-1, suggesting that HIV-1 survives close to its error threshold and may be readily susceptible to mutagenic drugs. The latter estimate was weakly dependent on the within-host effective population size of HIV-1. With large population sizes and in the absence of recombination, our simulations converged to the quasispecies theory, bridging the gap between quasispecies theory and population genetics-based approaches to describing HIV-1 evolution. Further, increased with the recombination rate, rendering HIV-1 less susceptible to error catastrophe, thus elucidating an added benefit of recombination to HIV-1. Our estimate of may serve as a quantitative guideline for the use of mutagenic drugs against HIV-1.

Currently available antiretroviral drugs curtail HIV infection but fail to eradicate the virus. A strategy of intervention radically different from that employed by current drugs has been proposed by the molecular quasispecies theory. The theory predicts that increasing the viral mutation rate beyond a critical value, called the error threshold, would cause a severe loss of genetic information, potentially leading to viral clearance. Several chemical mutagens are now being developed that can increase the mutation rate of HIV-1. Their success depends on reliable estimates of the error threshold of HIV-1, which are currently lacking. The quasispecies theory cannot be applied directly to HIV-1: the theory considers an infinitely large population of asexually reproducing haploid individuals, whereas HIV-1 is diploid, undergoes recombination, and is estimated to have a small effective population size in vivo. We employed detailed stochastic simulations that overcome the limitations of the quasispecies theory and accurately mimic HIV-1 evolution in vivo. With these simulations, we estimated the error threshold of HIV-1 to be ∼2–6-fold higher than its natural mutation rate, suggesting that HIV-1 survives close to its error threshold and may be readily susceptible to mutagenic drugs.

Genetic architecture of HIV-1 genes circulating in north India & their functional implications

This review presents data on genetic and functional analysis of some of the HIV-1 genes derived from HIV-1 infected individuals from north India (Delhi, Punjab and Chandigarh). We found evidence of novel B/C recombinants in HIV-1 LTR region showing relatedness to China/Myanmar with 3 copies of Nfκb sites; B/C/D mosaic genomes for HIV-1 Vpr and novel B/C Tat. We reported appearance of a complex recombinant form CRF_02AG of HIV-1 envelope sequences which is predominantly found in Central/Western Africa. Also one Indian HIV-1 envelope subtype C sequence suggested exclusive CXCR4 co-receptor usage. This extensive recombination, which is observed in about 10 per cent HIV-1 infected individuals in the Vpr genes, resulted in remarkably altered functions when compared with prototype subtype B Vpr. The Vpu C was found to be more potent in causing apoptosis when compared with Vpu B when analyzed for subG1 DNA content. The functional implications of these changes as well as in other genes of HIV-1 are discussed in detail with possible implications for subtype-specific pathogenesis highlighted.

Molecular characterization of HIV type 1 strains from newly diagnosed patients in Cyprus (2007-2009) recovers multiple clades including unique recombinant strains and lack of transmitted drug resistance

HIV-1 evolution generates substantial genetic diversity among isolates, the majority of which are represented in areas where multiple strains cocirculate. A heterogeneous genetic HIV-1 pool has been found in Cyprus, directing us to determine the dynamics of the local HIV-1 infection by characterizing strains isolated from 74 subjects during 2007-2009, representing 88% of the known-living HIV-1-infected population, of whom 53 are newly diagnosed therapy-naive patients and 21 are chronic patients, according to the European HIV Resistance guidelines. Near full-length genome sequences were amplified by RT-nested PCR using diluted RNA from all HIV-1 seropositives and sequenced using a newly designed assay. Resistant mutations were not found among the population of the newly diagnosed therapy-naive patients either to protease, reverse transcriptase, or integrase inhibitors. Phylogenetic analyses indicated subtype B as the main subtype (48.6%), followed by subtype A (18.9%), subtype C (10.8%), CRF02_AG (8.1%), CRF11_cpx (2.7%), and (sub)subtype F1 and CRF37_cpx (1.4% each). Six HIV-1 isolates (8.1%) were not classified in any pure (sub)subtype or circulating recombinant form (CRF). Complete phylogenetic and bootscanning analyses revealed that each isolate had a new, unique recombinant pattern and is distinct from all other CRFs or unique recombinant forms (URFs) reported so far. Two of the six isolates have the same mosaic pattern. Analogous to results of the earlier epidemiological studies, this study expands on the HIV-1 sequence database and reveals the high degree of diversity of HIV-1 infection in Cyprus.

Mechanisms and factors that influence high frequency retroviral recombination

With constantly changing environmental selection pressures, retroviruses rely upon recombination to reassort polymorphisms in their genomes and increase genetic diversity, which improves the chances for the survival of their population. Recombination occurs during DNA synthesis, whereby reverse transcriptase undergoes template switching events between the two copackaged RNAs, resulting in a viral recombinant with portions of the genetic information from each parental RNA. This review summarizes our current understanding of the factors and mechanisms influencing retroviral recombination, fidelity of the recombination process, and evaluates the subsequent viral diversity and fitness of the progeny recombinant. Specifically, the high mutation rates and high recombination frequencies of HIV-1 will be analyzed for their roles in influencing HIV-1 global diversity, as well as HIV-1 diagnosis, drug treatment, and vaccine development.

RecDraw: a software package for the representation of HIV-1 recombinant structures

The crucial role of recombination in HIV-1 biology is being increasingly recognized. In vitro studies have shown that up to 30 strand-transfer events may occur per viral replication cycle. Thus, recombination may surpass mutation as a major mechanism driving HIV-1 evolution. Currently, recombinant strains comprise 37% of the full-genome HIV-1 sequence database, including sequences representing 47 Circulating Recombinant Forms (CRFs) and more than 250 different Unique Recombinant Forms (URFs). Mapping of recombination breakpoints helps establish relationships among strains that are related by descent, such as CRF07_BC and CRF08_BC in China, and sheds light on their origin and epidemic spread. Additionally, unrelated recombinants sharing common breakpoints may reflect recombination hotspots within the viral genome. Here we present a software tool, RecDraw, for the graphical representation and efficient comparison of recombinant HIV-1 structures and breakpoints. RecDraw is a platform-flexible, Java stand-alone application available through http://www.hivresearch.org/research.php?ServiceID = 5&SubServiceID = 6 .

Experimental evidence of recombination in murine noroviruses

Based on sequencing data, norovirus (NoV) recombinants have been described, but no experimental evidence of recombination in NoVs has been documented. Using the murine norovirus (MNV) model, we investigated the occurrence of genetic recombination between two co-infecting wild-type MNV isolates in RAW cells. The design of a PCR-based genotyping tool allowed accurate discrimination between the parental genomes and the detection of a viable recombinant MNV (Rec MNV) in the progeny viruses. Genetic analysis of Rec MNV identified a homologous-recombination event located at the ORF1-ORF2 overlap. Rec MNV exhibited distinct growth curves and produced smaller plaques than the wild-type MNV in RAW cells. Here, we demonstrate experimentally that MNV undergoes homologous recombination at the previously described recombination hot spot for NoVs, suggesting that the MNV model might be suitable for in vitro studies of NoV recombination. Moreover, the results show that exchange of genetic material between NoVs can generate viruses with distinct biological properties from the parental viruses.

Longitudinal quasispecies analysis of viral variants in HIV type 1 dually infected individuals highlights the importance of sequence identity in viral recombination

Little is known regarding the likelihood of recombination between any given pair of nonidentical HIV-1 viruses in vivo. The present study analyzes the HIV-1 quasispecies in the C1C2 region of env, the vif-vpr-vpu accessory gene region, and the reverse transcriptase region of pol. These sequences were amplified from samples obtained sequentially over a 12- to 33-month period from five dually HIV-1-infected subjects. Analysis of an average of 248 clones amplified from each subject revealed no recombinants within the three loci studied of the subtype-discordant infecting strains, whose genetic diversity was >11% in env. In contrast, two subjects who were initially coinfected by two subtype-concordant variants with genetic diversity of 7.4% in env were found to harbor 10 unique recombinants of these strains, as exhibited by analysis of the env gene. The frequent recombination observed among the subtype-concordant strains studied herein correlates with prior sequence analyses that have commonly found higher rates of recombination at loci bearing the most conserved sequences, demonstrating an important role for sequence identity in HIV-1 recombination. Viral load analysis revealed that the samples studied contained an average of 8125 virus copies/ml (range, 882-31,626 copies/ml), signifying that the amount of viral RNA in the samples was not limiting for studying virus diversity. These data reveal that recombination between genetically distant strains may not be an immediate or common outcome to dual infection in vivo and suggest critical roles for viral and host factors such as viral fitness, virus diversity, and host immune responses that may contribute to limiting the frequency of intersubtype recombination during in vivo dual infection.

The remarkable frequency of human immunodeficiency virus type 1 genetic recombination

The genetic diversity of human immunodeficiency virus type 1 (HIV-1) results from a combination of point mutations and genetic recombination, and rates of both processes are unusually high. This review focuses on the mechanisms and outcomes of HIV-1 genetic recombination and on the parameters that make recombination so remarkably frequent. Experimental work has demonstrated that the process that leads to recombination--a copy choice mechanism involving the migration of reverse transcriptase between viral RNA templates--occurs several times on average during every round of HIV-1 DNA synthesis. Key biological factors that lead to high recombination rates for all retroviruses are the recombination-prone nature of their reverse transcription machinery and their pseudodiploid RNA genomes. However, HIV-1 genes recombine even more frequently than do those of many other retroviruses. This reflects the way in which HIV-1 selects genomic RNAs for coencapsidation as well as cell-to-cell transmission properties that lead to unusually frequent associations between distinct viral genotypes. HIV-1 faces strong and changeable selective conditions during replication within patients. The mode of HIV-1 persistence as integrated proviruses and strong selection for defective proviruses in vivo provide conditions for archiving alleles, which can be resuscitated years after initial provirus establishment. Recombination can facilitate drug resistance and may allow superinfecting HIV-1 strains to evade preexisting immune responses, thus adding to challenges in vaccine development. These properties converge to provide HIV-1 with the means, motive, and opportunity to recombine its genetic material at an unprecedented high rate and to allow genetic recombination to serve as one of the highest barriers to HIV-1 eradication.

Near full-length genetic analysis of HIV sequences derived from Cyprus: evidence of a highly polyphyletic and evolving infection

The molecular epidemiology of HIV-1 infection was previously studied in Cyprus but the degree of HIV-1 diversity has remained indefinable. The main objective of the present study is to examine HIV-1 strains isolated from 77 HIV-1-infected individuals representing 38% of the known infected population in Cyprus in the period 1986 to 2006. DNA of the near full-length genome encoding gag, pol, vif, vpr, vpu, tat, rev, env, and 5'-end of nef was amplified by nested PCR/RT-PCR from all HIV-1 seropositives and sequenced using a newly designed assay. Detailed phylogenetic and bootscanning analyses were performed to determine phylogenetic associations and subtype assignments. Phylogenetic analyses of the obtained viral sequences indicated that subtype B was the dominant subtype (61%), followed by subtype A (23.3%), subtype C (5.2%), CRF02_AG (3.9%), and subtype D, CRF01_AE, and CRF04_cpx (1.3% each). Two HIV-1 isolates (2.6%), originating from the Democratic Republic of Congo (DRC), were not classified in any pure (sub)subtype or circulating recombinant form (CRF). Complete phylogenetic and bootscanning analyses revealed that one of these isolates had a new, unique recombinant pattern, comprising segments of subtypes D and G, and is distinct from any other CRFs or URFs reported so far. Detailed analyses of the sequence of the second isolate, which could not be classified, reveal that it is close to subtype K reference sequences but clusters near the root of the clade. At least two epidemiologically unrelated HIV-1 seropositives with an HIV-1 variant similar to this isolate are required to designate this variant as a novel HIV-1 subtype or subsubtype of subtype K. Analogous to results of the earlier epidemiological studies, these data exhibit the extensive heterogeneity of HIV-1 infection in Cyprus, which is being fueled by a continuous entry of new strains from other countries, creating an evolving and polyphyletic infection.

Molecular mechanisms of recombination restriction in the envelope gene of the human immunodeficiency virus

The ability of pathogens to escape the host's immune response is crucial for the establishment of persistent infections and can influence virulence. Recombination has been observed to contribute to this process by generating novel genetic variants. Although distinctive recombination patterns have been described in many viral pathogens, little is known about the influence of biases in the recombination process itself relative to selective forces acting on newly formed recombinants. Understanding these influences is important for determining how recombination contributes to pathogen genome and proteome evolution. Most previous research on recombination-driven protein evolution has focused on relatively simple proteins, usually in the context of directed evolution experiments. Here, we study recombination in the envelope gene of HIV-1 between primary isolates belonging to subtypes that recombine naturally in the HIV/AIDS pandemic. By characterizing the early steps in the generation of recombinants, we provide novel insights into the evolutionary forces that shape recombination patterns within viral populations. Specifically, we show that the combined effects of mechanistic processes that determine the locations of recombination breakpoints across the HIV-1 envelope gene, and purifying selection acting against dysfunctional recombinants, can explain almost the entire distribution of breakpoints found within this gene in nature. These constraints account for the surprising paucity of recombination breakpoints found in infected individuals within this highly variable gene. Thus, the apparent randomness of HIV evolution via recombination may in fact be relatively more predictable than anticipated. In addition, the dominance of purifying selection in localized areas of the HIV genome defines regions where functional constraints on recombinants appear particularly strong, pointing to vulnerable aspects of HIV biology.

Recombination allows mixing portions of genomes of different origins, generating chimeric genes and genomes. With respect to the random generation of new mutations, it can lead to the simultaneous insertion of several substitutions, introducing more drastic changes in the genome. Furthermore, recombination is expected to yield a higher proportion of functional products since it combines variants that already exist in the population and that are therefore compatible with the survival of the organism. However, when recombination involves genetically distant strains, it can be constrained by the necessity to retain the functionality of the resulting products. In pathogens, which are subjected to strong selective pressures, recombination is particularly important, and several viruses, such as the human immunodeficiency virus (HIV), readily recombine. Here, we demonstrate the existence of preferential regions for recombination in the HIV-1 envelope gene when crossing sequences representative of strains observed to recombine in vivo. Furthermore, some recombinants give a decreased proportion of functional products. When considering these factors, one can retrace the history of most natural HIV recombinants. Recombination in HIV appears not so unpredictable, therefore, and the existence of recombinants that frequently generate nonfunctional products highlights previously unappreciated limits of the genetic flexibility of HIV.

Widely conserved recombination patterns among single-stranded DNA viruses

The combinatorial nature of genetic recombination can potentially provide organisms with immediate access to many more positions in sequence space than can be reached by mutation alone. Recombination features particularly prominently in the evolution of a diverse range of viruses. Despite rapid progress having been made in the characterization of discrete recombination events for many species, little is currently known about either gross patterns of recombination across related virus families or the underlying processes that determine genome-wide recombination breakpoint distributions observable in nature. It has been hypothesized that the networks of coevolved molecular interactions that define the epistatic architectures of virus genomes might be damaged by recombination and therefore that selection strongly influences observable recombination patterns. For recombinants to thrive in nature, it is probably important that the portions of their genomes that they have inherited from different parents work well together. Here we describe a comparative analysis of recombination breakpoint distributions within the genomes of diverse single-stranded DNA (ssDNA) virus families. We show that whereas nonrandom breakpoint distributions in ssDNA virus genomes are partially attributable to mechanistic aspects of the recombination process, there is also a significant tendency for recombination breakpoints to fall either outside or on the peripheries of genes. In particular, we found significantly fewer recombination breakpoints within structural protein genes than within other gene types. Collectively, these results imply that natural selection acting against viruses expressing recombinant proteins is a major determinant of nonrandom recombination breakpoint distributions observable in most ssDNA virus families.

Conserved positive selection signals in gp41 across multiple subtypes and difference in selection signals detectable in gp41 sequences sampled during acute and chronic HIV-1 subtype C infection

BACKGROUND

The high diversity of HIV variants driving the global AIDS epidemic has caused many to doubt whether an effective vaccine against the virus is possible. However, by identifying the selective forces that are driving the ongoing diversification of HIV and characterising their genetic consequences, it may be possible to design vaccines that pre-empt some of the virus' more common evasion tactics. One component of such vaccines might be the envelope protein, gp41. Besides being targeted by both the humoral and cellular arms of the immune system this protein mediates fusion between viral and target cell membranes and is likely to be a primary determinant of HIV transmissibility.

RESULTS

Using recombination aware analysis tools we compared site specific signals of selection in gp41 sequences from different HIV-1 M subtypes and circulating recombinant forms and identified twelve sites evolving under positive selection across multiple major HIV-1 lineages. To identify evidence of selection operating during transmission our analysis included two matched datasets sampled from patients with acute or chronic subtype C infections. We identified six gp41 sites apparently evolving under different selection pressures during acute and chronic HIV-1 infections. These sites mostly fell within functional gp41 domains, with one site located within the epitope recognised by the broadly neutralizing antibody, 4E10.

CONCLUSION

Whereas these six sites are potentially determinants of fitness and are therefore good candidate targets for subtype-C specific vaccines, the twelve sites evolving under diversifying selection across multiple subtypes might make good candidate targets for broadly protective vaccines.

Identifying the important HIV-1 recombination breakpoints

Recombinant HIV-1 genomes contribute significantly to the diversity of variants within the HIV/AIDS pandemic. It is assumed that some of these mosaic genomes may have novel properties that have led to their prevalence, particularly in the case of the circulating recombinant forms (CRFs). In regions of the HIV-1 genome where recombination has a tendency to convey a selective advantage to the virus, we predict that the distribution of breakpoints—the identifiable boundaries that delimit the mosaic structure—will deviate from the underlying null distribution. To test this hypothesis, we generate a probabilistic model of HIV-1 copy-choice recombination and compare the predicted breakpoint distribution to the distribution from the HIV/AIDS pandemic. Across much of the HIV-1 genome, we find that the observed frequencies of inter-subtype recombination are predicted accurately by our model. This observation strongly indicates that in these regions a probabilistic model, dependent on local sequence identity, is sufficient to explain breakpoint locations. In regions where there is a significant over- (either side of the env gene) or under- (short regions within gag, pol, and most of env) representation of breakpoints, we infer natural selection to be influencing the recombination pattern. The paucity of recombination breakpoints within most of the envelope gene indicates that recombinants generated in this region are less likely to be successful. The breakpoints at a higher frequency than predicted by our model are approximately at either side of env, indicating increased selection for these recombinants as a consequence of this region, or at least part of it, having a tendency to be recombined as an entire unit. Our findings thus provide the first clear indication of the existence of a specific portion of the genome that deviates from a probabilistic null model for recombination. This suggests that, despite the wide diversity of recombinant forms seen in the viral population, only a minority of recombination events appear to be of significance to the evolution of HIV-1.

Multiple variants of HIV can infect the same cell, and because each viral particle contains two copies of the viral genomic RNA, RNAs from different viruses can occasionally be incorporated together within a viral particle. When this virus subsequently infects another cell, genetic exchange (recombination) may occur between these two divergent copies of genomic RNA as a result of a switch between the RNA molecules while they are copied into DNA. This process is very important to understand as it contributes to the generation of new HIV variants. In this study, we have analysed a set of recombinant HIV genomes generated in the laboratory to construct a probabilistic model of the propensity for the switch to take place in specific regions of the genome, dependent on the local similarity of the parental viral sequences. This model allows us to predict the locations where recombination should occur more frequently. By comparing these predictions to the patterns of recombination observed in the HIV-1 pandemic, we identify the genomic regions in which recombination has been more important, in that it has provided an evolutionary advantage to the virus.

Implications of recombination for HIV diversity

The human immunodeficiency virus (HIV) population is characterised by extensive genetic variability that results from high error and recombination rates of the reverse transcription process, and from the fast turnover of virions in HIV-infected individuals. Among the viral variants encountered at the global scale, recombinant forms are extremely abundant. Some of these recombinants (known as circulating recombinant forms) become fixed and undergo rapid expansion in the population. The reasons underlying their epidemiological success remain at present poorly understood and constitute a fascinating area for future research to improve our understanding of immune escape, pathogenicity and transmission. Recombinant viruses are generated during reverse transcription as a consequence of template switching between the two genetically different genomic RNAs present in a heterozygous virus. Recombination can thereby generate shortcuts in evolution by producing mosaic reverse transcription products of parental genomes. Therefore, in a single infectious cycle multiple mutations that are positively selected can be combined or, conversely, negatively selected mutations can be removed. Recombination is therefore involved in different aspects of HIV evolution, adaptation to its host, and escape from antiviral treatments.