Estimating selection pressures on HIV-1 using phylogenetic likelihood models

Human immunodeficiency virus (HIV-1) can rapidly evolve due to selection pressures exerted by HIV-specific immune responses, antiviral agents, and to allow the virus to establish infection in different compartments in the body. Statistical models applied to HIV-1 sequence data can help to elucidate the nature of these selection pressures through comparisons of non-synonymous (or amino acid changing) and synonymous (or amino acid preserving) substitution rates. These models also need to take into account the non-independence of sequences due to their shared evolutionary history. We review how we have developed these methods and have applied them to characterize the evolution of HIV-1 in vivo.To illustrate our methods, we present an analysis of compartment-specific evolution of HIV-1 env in blood and cerebrospinal fluid and of site-to-site variation in the gag gene of subtype C HIV-1

Enrichment of intersubtype HIV-1 recombinants in a dual infection system using HIV-1 strain-specific siRNAs

<p>Abstract</p> <p>Background</p> <p>Intersubtype HIV-1 recombinants in the form of unique or stable circulating recombinants forms (CRFs) are responsible for over 20% of infections in the worldwide epidemic. Mechanisms controlling the generation, selection, and transmission of these intersubtype HIV-1 recombinants still require further investigation. All intersubtype HIV-1 recombinants are generated and evolve from initial dual infections, but are difficult to identify in the human population. In vitro studies provide the most practical system to study mechanisms, but the recombination rates are usually very low in dual infections with primary HIV-1 isolates. This study describes the use of HIV-1 isolate-specific siRNAs to enrich intersubtype HIV-1 recombinants and inhibit the parental HIV-1 isolates from a dual infection.</p> <p>Results</p> <p>Following a dual infection with subtype A and D primary HIV-1 isolates and two rounds of siRNA treatment, nearly 100% of replicative virus was resistant to a siRNA specific for an upstream target sequence in the subtype A envelope (<it>env</it>) gene as well as a siRNA specific for a downstream target sequence in the subtype D <it>env </it>gene. Only 20% (10/50) of the replicating virus had nucleotide substitutions in the siRNA-target sequence whereas the remaining 78% (39/50) harbored a recombination breakpoint that removed both siRNA target sequences, and rendered the intersubtype D/A recombinant virus resistant to the dual siRNA treatment. Since siRNAs target the newly transcribed HIV-1 mRNA, the siRNAs only enrich intersubtype env recombinants and do not influence the recombination process during reverse transcription. Using this system, a strong bias is selected for recombination breakpoints in the C2 region, whereas other HIV-1 env regions, most notably the hypervariable regions, were nearly devoid of intersubtype recombination breakpoints. Sequence conservation plays an important role in selecting for recombination breakpoints, but the lack of breakpoints in many conserved env regions suggest that other mechanisms are at play.</p> <p>Conclusion</p> <p>These findings show that siRNAs can be used as an efficient in vitro tool for enriching recombinants, to facilitate further study on mechanisms of intersubytpe HIV-1 recombination, and to generate replication-competent intersubtype recombinant proteins with a breadth in HIV-1 diversity for future vaccine studies.</p

A novel exploratory method for visual recombination detection

A versatile visual approach for detecting recombination and identifying recombination breakpoints within a sequence alignment is presented. The method is based on two novel diagrams - the highway plot and the occupancy plot - that graphically portray phylogenetic inhomogeneity along an alignment, and can be viewed as a synthesis of two widely used but unrelated methods: bootscanning and quartet-mapping. To illustrate the method, simulated data and HIV-1 and influenza A datasets are investigated

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

Analysis of the Origin and Evolutionary History of HIV-1 CRF28_BF and CRF29_BF Reveals a Decreasing Prevalence in the AIDS Epidemic of Brazil

HIV-1 subtype B and subtype F are prevalent in the AIDS epidemic of Brazil. Recombinations between these subtypes have generated at least four BF circulating recombinant forms (CRFs). CRF28_BF and CRF29_BF are among the first two BF recombinants being identified in Brazil and they contributed significantly to the epidemic. However, the evolution and demographic histories of the CRFs are unclear.A collection of gag and pol sequences sampled within Brazil was screened for CRF28_BF-like and CRF29_BF-like recombination patterns. A Bayesian coalescent framework was employed to delineate the phylogenetic, divergence time and population dynamics of the virus having CRF28_BF-like and CRF29_BF-like genotype. These recombinants were phylogenetically related to each other and formed a well-supported monophyletic clade dated to 1988-1989. The effective number of infections by these recombinants grew exponentially over a five-year period after their emergence, but then decreased toward the present following a logistic model of population growth. The demographic pattern of both recombinants closely resembles those previously reported for CRF31_BC.We revealed that HIV-1 recombinants of the CRF28_BF/CRF29_BF clade are still circulating in the Brazilian population. These recombinants did not exhibit a strong founder effect and showed a decreasing prevalence in the AIDS epidemic of Brazil. Our data suggested that multiple URFs may also play a role in shaping the epidemic of recombinant BF HIV-1 in the region

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

<p>Abstract</p> <p>Background</p> <p>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.</p> <p>Results</p> <p>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.</p> <p>Conclusion</p> <p>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.</p

A jumping profile Hidden Markov Model and applications to recombination sites in HIV and HCV genomes

BACKGROUND: Jumping alignments have recently been proposed as a strategy to search a given multiple sequence alignment A against a database. Instead of comparing a database sequence S to the multiple alignment or profile as a whole, S is compared and aligned to individual sequences from A. Within this alignment, S can jump between different sequences from A, so different parts of S can be aligned to different sequences from the input multiple alignment. This approach is particularly useful for dealing with recombination events. RESULTS: We developed a jumping profile Hidden Markov Model (jpHMM), a probabilistic generalization of the jumping-alignment approach. Given a partition of the aligned input sequence family into known sequence subtypes, our model can jump between states corresponding to these different subtypes, depending on which subtype is locally most similar to a database sequence. Jumps between different subtypes are indicative of intersubtype recombinations. We applied our method to a large set of genome sequences from human immunodeficiency virus (HIV) and hepatitis C virus (HCV) as well as to simulated recombined genome sequences. CONCLUSION: Our results demonstrate that jumps in our jumping profile HMM often correspond to recombination breakpoints; our approach can therefore be used to detect recombinations in genomic sequences. The recombination breakpoints identified by jpHMM were found to be significantly more accurate than breakpoints defined by traditional methods based on comparing single representative sequences

Accurate Detection of Recombinant Breakpoints in Whole-Genome Alignments

We propose a novel method for detecting sites of molecular recombination in multiple alignments. Our approach is a compromise between previous extremes of computationally prohibitive but mathematically rigorous methods and imprecise heuristic methods. Using a combined algorithm for estimating tree structure and hidden Markov model parameters, our program detects changes in phylogenetic tree topology over a multiple sequence alignment. We evaluate our method on benchmark datasets from previous studies on two recombinant pathogens, Neisseria and HIV-1, as well as simulated data. We show that we are not only able to detect recombinant regions of vastly different sizes but also the location of breakpoints with great accuracy. We show that our method does well inferring recombination breakpoints while at the same time maintaining practicality for larger datasets. In all cases, we confirm the breakpoint predictions of previous studies, and in many cases we offer novel predictions