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

Mutation of HIV-1 Genomes in a Clinical Population Treated with the Mutagenic Nucleoside KP1461

The deoxycytidine analog KP1212, and its prodrug KP1461, are prototypes of a new class of antiretroviral drugs designed to increase viral mutation rates, with the goal of eventually causing the collapse of the viral population. Here we present an extensive analysis of viral sequences from HIV-1 infected volunteers from the first “mechanism validation” phase II clinical trial of a mutagenic base analog in which individuals previously treated with antiviral drugs received 1600 mg of KP1461 twice per day for 124 days. Plasma viral loads were not reduced, and overall levels of viral mutation were not increased during this short-term study, however, the mutation spectrum of HIV was altered. A large number (N = 105 per sample) of sequences were analyzed, each derived from individual HIV-1 RNA templates, after 0, 56 and 124 days of therapy from 10 treated and 10 untreated control individuals (>7.1 million base pairs of unique viral templates were sequenced). We found that private mutations, those not found in more than one viral sequence and likely to have occurred in the most recent rounds of replication, increased in treated individuals relative to controls after 56 (p = 0.038) and 124 (p = 0.002) days of drug treatment. The spectrum of mutations observed in the treated group showed an excess of A to G and G to A mutations (p = 0.01), and to a lesser extent T to C and C to T mutations (p = 0.09), as predicted by the mechanism of action of the drug. These results validate the proposed mechanism of action in humans and should spur development of this novel antiretroviral approach.Koronis Pharmaceutical

GENOMEPOP: A program to simulate genomes in populations

<p>Abstract</p> <p>Background</p> <p>There are several situations in population biology research where simulating DNA sequences is useful. Simulation of biological populations under different evolutionary genetic models can be undertaken using backward or forward strategies. Backward simulations, also called coalescent-based simulations, are computationally efficient. The reason is that they are based on the history of lineages with surviving offspring in the current population. On the contrary, forward simulations are less efficient because the entire population is simulated from past to present. However, the coalescent framework imposes some limitations that forward simulation does not. Hence, there is an increasing interest in forward population genetic simulation and efficient new tools have been developed recently. Software tools that allow efficient simulation of large DNA fragments under complex evolutionary models will be very helpful when trying to better understand the trace left on the DNA by the different interacting evolutionary forces. Here I will introduce GenomePop, a forward simulation program that fulfills the above requirements. The use of the program is demonstrated by studying the impact of intracodon recombination on global and site-specific <it>dN/dS </it>estimation.</p> <p>Results</p> <p>I have developed algorithms and written software to efficiently simulate, forward in time, different Markovian nucleotide or codon models of DNA mutation. Such models can be combined with recombination, at inter and intra codon levels, fitness-based selection and complex demographic scenarios.</p> <p>Conclusion</p> <p>GenomePop has many interesting characteristics for simulating SNPs or DNA sequences under complex evolutionary and demographic models. These features make it unique with respect to other simulation tools. Namely, the possibility of forward simulation under General Time Reversible (GTR) mutation or GTR×MG94 codon models with intra-codon recombination, arbitrary, user-defined, migration patterns, diploid or haploid models, constant or variable population sizes, etc. It also allows simulation of fitness-based selection under different distributions of mutational effects. Under the 2-allele model it allows the simulation of recombination hot-spots, the definition of different frequencies in different populations, etc. GenomePop can also manage large DNA fragments. In addition, it has a scaling option to save computation time when simulating large sequences and population sizes under complex demographic and evolutionary situations. These and many other features are detailed in its web page <abbrgrp><abbr bid="B1">1</abbr></abbrgrp>.</p

The Evolutionary Dynamics of a Rapidly Mutating Virus within and between Hosts: The Case of Hepatitis C Virus

Many pathogens associated with chronic infections evolve so rapidly that strains found late in an infection have little in common with the initial strain. This raises questions at different levels of analysis because rapid within-host evolution affects the course of an infection, but it can also affect the possibility for natural selection to act at the between-host level. We present a nested approach that incorporates within-host evolutionary dynamics of a rapidly mutating virus (hepatitis C virus) targeted by a cellular cross-reactive immune response, into an epidemiological perspective. The viral trait we follow is the replication rate of the strain initiating the infection. We find that, even for rapidly evolving viruses, the replication rate of the initial strain has a strong effect on the fitness of an infection. Moreover, infections caused by slowly replicating viruses have the highest infection fitness (i.e., lead to more secondary infections), but strains with higher replication rates tend to dominate within a host in the long-term. We also study the effect of cross-reactive immunity and viral mutation rate on infection life history traits. For instance, because of the stochastic nature of our approach, we can identify factors affecting the outcome of the infection (acute or chronic infections). Finally, we show that anti-viral treatments modify the value of the optimal initial replication rate and that the timing of the treatment administration can have public health consequences due to within-host evolution. Our results support the idea that natural selection can act on the replication rate of rapidly evolving viruses at the between-host level. It also provides a mechanistic description of within-host constraints, such as cross-reactive immunity, and shows how these constraints affect the infection fitness. This model raises questions that can be tested experimentally and underlines the necessity to consider the evolution of quantitative traits to understand the outcome and the fitness of an infection

Evolutionary Genomics of a Temperate Bacteriophage in an Obligate Intracellular Bacteria (Wolbachia)

Genome evolution of bacteria is usually influenced by ecology, such that bacteria with a free-living stage have large genomes and high rates of horizontal gene transfer, while obligate intracellular bacteria have small genomes with typically low amounts of gene exchange. However, recent studies indicate that obligate intracellular species that host-switch frequently harbor agents of horizontal transfer such as mobile elements. For example, the temperate double-stranded DNA bacteriophage WO in Wolbachia persistently transfers between bacterial coinfections in the same host. Here we show that despite the phage's rampant mobility between coinfections, the prophage's genome displays features of constraint related to its intracellular niche. First, there is always at least one intact prophage WO and usually several degenerate, independently-acquired WO prophages in each Wolbachia genome. Second, while the prophage genomes are modular in composition with genes of similar function grouping together, the modules are generally not interchangeable with other unrelated phages and thus do not evolve by the Modular Theory. Third, there is an unusual core genome that strictly consists of head and baseplate genes; other gene modules are frequently deleted. Fourth, the prophage recombinases are diverse and there is no conserved integration sequence. Finally, the molecular evolutionary forces acting on prophage WO are point mutation, intragenic recombination, deletion, and purifying selection. Taken together, these analyses indicate that while lateral transfer of phage WO is pervasive between Wolbachia with occasional new gene uptake, constraints of the intracellular niche obstruct extensive mixture between WO and the global phage population. Although the Modular Theory has long been considered the paradigm of temperate bacteriophage evolution in free-living bacteria, it appears irrelevant in phages of obligate intracellular bacteria

Surface-modified polyacrylonitrile nanofibers as supports

Polyacrylonitrile nanofibers (PAN-nfs) are one of the most studied nanofibres because of their excellent characteristics, such as good mechanical strength, chemical resistance, and good thermal stability. Due to the easy dissolution in polar organic solvents, PAN-nfs are mostly produced via electrospinning technique. The electrospun PAN-nfs surfaces are relatively in-active and hydrophobic, and, therefore, hinder some potential applications; however, chemical surface modification reactions, such as amination, reduction, hydrolysis, and amidoximation, have been carried out on them. These reactions bring about functional groups, such as amine, hydroxyl, carboxylic, imine etc, to the surface PAN-nfs and invariably make their surfaces active and hydrophilic. The surface-modified PAN-nfs have been used as supports for organic compounds, enzymes, and antibodies in biological studies. They have also been used for immobilization of various organic ligands for adsorption of metal ions in water. Furthermore, because of their ability to complex metal ions, several surface-modified PAN-nfs have also been used as supports for transition metal catalysts in Fenton’s chemistry.IS

Computational analyses of an evolutionary arms race between mammalian immunity mediated by immunoglobulin A and its subversion by bacterial pathogens

IgA is the predominant immunoglobulin isotype in mucosal tissues and external secretions, playing important roles both in defense against pathogens and in maintenance of commensal microbiota. Considering the complexity of its interactions with the surrounding environment, IgA is a likely target for diversifying or positive selection. To investigate this possibility, the action of natural selection on IgA was examined in depth with six different methods: CODEML from the PAML package and the SLAC, FEL, REL, MEME and FUBAR methods implemented in the Datamonkey webserver. In considering just primate IgA, these analyses show that diversifying selection targeted five positions of the Cα1 and Cα2 domains of IgA. Extending the analysis to include other mammals identified 18 positively selected sites: ten in Cα1, five in Cα2 and three in Cα3. All but one of these positions display variation in polarity and charge. Their structural locations suggest they indirectly influence the conformation of sites on IgA that are critical for interaction with host IgA receptors and also with proteins produced by mucosal pathogens that prevent their elimination by IgA-mediated effector mechanisms. Demonstrating the plasticity of IgA in the evolution of different groups of mammals, only two of the eighteen selected positions in all mammals are included in the five selected positions in primates. That IgA residues subject to positive selection impact sites targeted both by host receptors and subversive pathogen ligands highlights the evolutionary arms race playing out between mammals and pathogens, and further emphasizes the importance of IgA in protection against mucosal pathogens.The Portuguese Foundation for Science and Technology (www.fct.pt) supported the doctoral fellowship of Ana Pinheiro (SFRH/BD/71252/2010) and the post-doctoral fellowship of Pedro J. Esteves (SPRH/BPD/27021/2006). This work was also supported by a project of the Portuguese Foundation for Science and Technology (PTDC/BIA-BEC/103158/2008).Peer Reviewe