Recombination rate and selection strength in HIV intra-patient evolution

The evolutionary dynamics of HIV during the chronic phase of infection is driven by the host immune response and by selective pressures exerted through drug treatment. To understand and model the evolution of HIV quantitatively, the parameters governing genetic diversification and the strength of selection need to be known. While mutation rates can be measured in single replication cycles, the relevant effective recombination rate depends on the probability of coinfection of a cell with more than one virus and can only be inferred from population data. However, most population genetic estimators for recombination rates assume absence of selection and are hence of limited applicability to HIV, since positive and purifying selection are important in HIV evolution. Here, we estimate the rate of recombination and the distribution of selection coefficients from time-resolved sequence data tracking the evolution of HIV within single patients. By examining temporal changes in the genetic composition of the population, we estimate the effective recombination to be r=1.4e-5 recombinations per site and generation. Furthermore, we provide evidence that selection coefficients of at least 15% of the observed non-synonymous polymorphisms exceed 0.8% per generation. These results provide a basis for a more detailed understanding of the evolution of HIV. A particularly interesting case is evolution in response to drug treatment, where recombination can facilitate the rapid acquisition of multiple resistance mutations. With the methods developed here, more precise and more detailed studies will be possible, as soon as data with higher time resolution and greater sample sizes is available.Comment: to appear in PLoS Computational Biolog

The Evolutionary Analysis of Emerging Low Frequency HIV-1 CXCR4 Using Variants through Time—An Ultra-Deep Approach

Large-scale parallel pyrosequencing produces unprecedented quantities of sequence data. However, when generated from viral populations current mapping software is inadequate for dealing with the high levels of variation present, resulting in the potential for biased data loss. In order to apply the 454 Life Sciences' pyrosequencing system to the study of viral populations, we have developed software for the processing of highly variable sequence data. Here we demonstrate our software by analyzing two temporally sampled HIV-1 intra-patient datasets from a clinical study of maraviroc. This drug binds the CCR5 coreceptor, thus preventing HIV-1 infection of the cell. The objective is to determine viral tropism (CCR5 versus CXCR4 usage) and track the evolution of minority CXCR4-using variants that may limit the response to a maraviroc-containing treatment regimen. Five time points (two prior to treatment) were available from each patient. We first quantify the effects of divergence on initial read k-mer mapping and demonstrate the importance of utilizing population-specific template sequences in relation to the analysis of next-generation sequence data. Then, in conjunction with coreceptor prediction algorithms that infer HIV tropism, our software was used to quantify the viral population structure pre- and post-treatment. In both cases, low frequency CXCR4-using variants (2.5–15%) were detected prior to treatment. Following phylogenetic inference, these variants were observed to exist as distinct lineages that were maintained through time. Our analysis, thus confirms the role of pre-existing CXCR4-using virus in the emergence of maraviroc-insensitive HIV. The software will have utility for the study of intra-host viral diversity and evolution of other fast evolving viruses, and is available from http://www.bioinf.manchester.ac.uk/segminator/

Sequence determinants of breakpoint location during HIV-1 intersubtype recombination

Retroviral recombination results from strand switching, during reverse transcription, between the two copies of genomic RNA present in the virus. We analysed recombination in part of the envelope gene, between HIV-1 subtype A and D strains. After a single infection cycle, breakpoints clustered in regions corresponding to the constant portions of Env. With some exceptions, a similar distribution was observed after multiple infection cycles, and among recombinant sequences in the HIV Sequence Database. We compared the experimental data with computer simulations made using a program that only allows recombination to occur whenever an identical base is present in the aligned parental RNAs. Experimental recombination was more frequent than expected on the basis of simulated recombination when, in a region spanning 40 nt from the 5′ border of a breakpoint, no more than two discordant bases between the parental RNAs were present. When these requirements were not fulfilled, breakpoints were distributed randomly along the RNA, closer to the distribution predicted by computer simulation. A significant preference for recombination was also observed for regions containing homopolymeric stretches. These results define, for the first time, local sequence determinants for recombination between divergent HIV-1 isolates

Identifying recombination hotspots in the HIV-1 genome

HIV-1 infection is characterised by the rapid generation of genetic diversity that facilitates viral escape from immune selection and antiretroviral therapy. Despite recombination's crucial role in viral diversity and evolution, little is known about the genomic factors that influence recombination between highly similar genomes. In this study, we use a minimally modified full length HIV-1 genome and high throughput sequence analysis to study recombination in gag and pol in T cells. We find that recombination is favoured at a number of recombination hotspots, where recombination occurs six times more frequently than at corresponding coldspots. Interestingly, these hotspots occur near important features of the HIV-1 genome, but do not occur at sites immediately around protease inhibitor or reverse transcriptase inhibitor drug resistance mutations. We show that the recombination hot and cold spots are consistent across five blood donors and are independent of co-receptor mediated entry. Finally, we check common experimental confounders and find that these are not driving the location of recombination hotspots. This is the first study to identify the location of recombination hotspots, between two similar viral genomes with great statistical power and under conditions that closely reflect natural recombination events amongst HIV-1 quasispecies

Inhibition of Akt sensitises neuroblastoma cells to gold(III) porphyrin 1a, a novel antitumour drug induced apoptosis and growth inhibition

Background:Gold(III) porphyrin 1a is a new class of anticancer drug, which inhibits cell proliferation of wide range of human cancer cell lines and induces apoptosis in human nasopharyngeal carcinoma cells. However, the underlying signalling mechanism by which gold(III) porphyrin 1a modifies the intracellular apoptosis pathways in tumour cells has not been explained in detail in neuroblastoma cells.Methods:Cell proliferation and apoptosis were determined by measuring 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) and Annexin V binding, respectively. Western blot assay was used to detect proteins involved in apoptotic and Akt pathways. In vivo tumour growth was assessed by inoculating tumour cells to nude mice subcutaneously, and gold(III) porphyrin 1a was administrated intravenously.Results:This study assessed the antitumour effect and mechanism of gold(III) porphyrin 1a on neuroblastoma in vitro and in vivo. Gold(III) porphyrin 1a displayed a growth inhibition and induction of apoptosis in neuroblastoma cells effectively in vitro, which was accompanied with release of cytochrome c and Smac/DIABLO and caspases activation. Further studies indicated that gold(III) porphyrin 1a inhibited X-linked inhibitor of apoptosis (XIAP). However, we found that gold(III) porphyrin 1a can induce a survival signal, Akt activation within minutes and could last for at least 24 h. To further confirm association between activation of Akt and the effectiveness of gold(III) porphyrin 1a, neuroblastoma cells were treated with API-2, an Akt-specific inhibitor. API-2 sensitised cells to gold(III) porphyrin 1a-induced apoptosis and growth inhibition.Conclusion:These results suggested that Akt may be considered as a molecular brake that neuroblastoma cells rely on to slow down gold(III) porphyrin 1a-induced apoptosis and antiproliferation. Gold(III) porphyrin 1a is a mitochondrial apoptotic stimulus but also activates Akt, suggesting an involvement of Akt in mediating the effectiveness to growth inhibition and apoptosis by gold(III) porphyrin 1a and that inhibition of Akt can enhance the anticancer activity of gold(III) porphyrin 1a in neuroblastoma. © 2009 Cancer Research UK.published_or_final_versio

Accurately Measuring Recombination between Closely Related HIV-1 Genomes

Retroviral recombination is thought to play an important role in the generation of immune escape and multiple drug resistance by shuffling pre-existing mutations in the viral population. Current estimates of HIV-1 recombination rates are derived from measurements within reporter gene sequences or genetically divergent HIV sequences. These measurements do not mimic the recombination occurring in vivo, between closely related genomes. Additionally, the methods used to measure recombination make a variety of assumptions about the underlying process, and often fail to account adequately for issues such as co-infection of cells or the possibility of multiple template switches between recombination sites. We have developed a HIV-1 marker system by making a small number of codon modifications in gag which allow recombination to be measured over various lengths between closely related viral genomes. We have developed statistical tools to measure recombination rates that can compensate for the possibility of multiple template switches. Our results show that when multiple template switches are ignored the error is substantial, particularly when recombination rates are high, or the genomic distance is large. We demonstrate that this system is applicable to other studies to accurately measure the recombination rate and show that recombination does not occur randomly within the HIV genome

The Role of Recombination for the Coevolutionary Dynamics of HIV and the Immune Response

The evolutionary implications of recombination in HIV remain not fully understood. A plausible effect could be an enhancement of immune escape from cytotoxic T lymphocytes (CTLs). In order to test this hypothesis, we constructed a population dynamic model of immune escape in HIV and examined the viral-immune dynamics with and without recombination. Our model shows that recombination (i) increases the genetic diversity of the viral population, (ii) accelerates the emergence of escape mutations with and without compensatory mutations, and (iii) accelerates the acquisition of immune escape mutations in the early stage of viral infection. We see a particularly strong impact of recombination in systems with broad, non-immunodominant CTL responses. Overall, our study argues for the importance of recombination in HIV in allowing the virus to adapt to changing selective pressures as imposed by the immune system and shows that the effect of recombination depends on the immunodominance pattern of effector T cell responses

In silico modeling indicates the development of HIV-1 resistance to multiple shRNA gene therapy differs to standard antiretroviral therapy

<p>Abstract</p> <p>Background</p> <p>Gene therapy has the potential to counter problems that still hamper standard HIV antiretroviral therapy, such as toxicity, patient adherence and the development of resistance. RNA interference can suppress HIV replication as a gene therapeutic via expressed short hairpin RNAs (shRNAs). It is now clear that multiple shRNAs will likely be required to suppress infection and prevent the emergence of resistant virus.</p> <p>Results</p> <p>We have developed the first biologically relevant stochastic model in which multiple shRNAs are introduced into CD34+ hematopoietic stem cells. This model has been used to track the production of gene-containing CD4+ T cells, the degree of HIV infection, and the development of HIV resistance in lymphoid tissue for 13 years. In this model, we found that at least four active shRNAs were required to suppress HIV infection/replication effectively and prevent the development of resistance. The inhibition of incoming virus was shown to be critical for effective treatment. The low potential for resistance development that we found is largely due to a pool of replicating wild-type HIV that is maintained in non-gene containing CD4+ T cells. This wild-type HIV effectively out-competes emerging viral strains, maintaining the viral <it>status quo</it>.</p> <p>Conclusions</p> <p>The presence of a group of cells that lack the gene therapeutic and is available for infection by wild-type virus appears to mitigate the development of resistance observed with systemic antiretroviral therapy.</p