Reactivation of ancestral strains of HIV-1 in the gp120 V3 env region in patients failing antiretroviral therapy and subjected to structured treatment interruption.

We analyzed gp120V3 HIV-1 env region genetic diversity of 27 patients failing antiretrovirals and subjected to 12-week structured treatment interruption (STI). Based on heteroduplex mobility assays, eight patients presented low pre- and post-STI genetic diversity (G1); five presented high pre-STI but low post-STI diversity (G2); five presented low pre-STI and high post-STI diversity (G3); and nine, high pre- and post-STI diversity (G4). One patient from G1, two from G2 and two from G3 were subjected to proviral DNA end-point PCR and sequencing. In three patients, the dramatic disturbance caused by STI resulted in ancestral viral progeny activation, which repopulated the cell reservoir. In two patients presenting highly homogeneous sequences and low immune selective pressure (dN/dS ratio <1), this phenomenon was not observed. The mechanisms involved in viral evolution, in which antiretroviral therapy also applies selective pressure, sometimes affects coreceptor usage of circulating viruses, leading to the suppression of x4 strains.

Fungi evolution revisited: application of the penalized likelihood method to a Bayesian fungal phylogeny provides a new perspective on phylogenetic relationships and divergence dates of Ascomycota groups.

The depiction of evolutionary relationships within phylum Ascomycota is still controversial because of unresolved branching orders in the radiation of major taxa. Here we generated a dataset of 166 small subunit (18S) rDNA sequences, representative of all groups of Fungi and used as input in a Bayesian phylogenetic analysis. This phylogeny suggests that Discomycetes are a basal group of filamentous Ascomycetes and probably maintain ancestor characters since their representatives are intermingled among other filamentous fungi. Also, we show that the evolutionary rate heterogeneity within Ascomycota precludes the assumption of a global molecular clock. Accordingly, we used the penalized likelihood method, and for calibration we included a 400 million-year-old Pyrenomycete fossil considering two distinct scenarios found in the literature, one with an estimated date of 1576 Myr for the plant-animal-fungus split and the other with an estimated date of 965 Myr for the animal-fungus split. Our data show that the current classification of the fossil as a Pyrenomycete is not compatible with the second scenario. Estimates under the first scenario are older than dates proposed in previous studies based on small subunit rDNA sequences but support estimates based on multiprotein analysis, suggesting that the radiation of the major Ascomycota groups occurred into the Proterozoic era.

Analysis of full-length human immunodeficiency virus type 1 genome reveals a variable spectrum of subtypes B and f recombinants in São Paulo, Brazil.

Recombination is one of the major mechanisms contributing to human immunodeficiency virus type 1 (HIV-1) variability. Analysis of pol gene sequences of 215 HIV-1 samples from São Paulo, Brazil classified 189 sequences as subtype B (87.9%), 8 sequences as subtype F (3.7%), and 18 sequences (8.4%) as B/F recombinants. After the analysis of the pol gene, a subset of six recombinant samples composed of sequences with a related recombinant pol structure was selected for full-length genome analysis to identify a possible circulating recombinant form. According to full-length genome analysis, recombination was higher in gag, protease, reverse transcriptase, integrase, and vif. Identification of many distinct recombinant forms and the absence of an identifiable HIV-1 circulating recombinant form suggest that a high frequency of dual infections between HIV-1 subtypes B and F is occurring in São Paulo, Brazil.

Changes in G+C content of a neutrally evolving gene under a non-reversible dynamics measured by computer simulations based on experimental evolution data

To evaluate the effects of non-reversibility on compositional base changes and the distribution of branch lengths along a phylogeny, we extended, by means of computer simulations, our previous sequential PCR in vitro evolution experiment. In that study a 18S rRNA gene evolved neutrally for 280 generations and a homogeneous non-stationary model of base substitution based on a non-reversible dynamics was built from the in vitro evolution data to describe the observed pattern of nucleotide substitutions. Here, the process was extended to 840 generations without selection, using the model parameters calculated from the in vitro evolution experiment. We observed that under a non-reversible model the G+C content of the sequences significantly increases when compared to simulations with a reversible model. The values of mean and variance of the branch lengths are reduced under a non-reversible dynamics although they follow a Poisson distribution. We conclude that the major implication of non-reversibility is the overall decrease of branch lengths, although no transition from a stochastic to an ordered process is observed. According to our model the result of this neutral process will be the increase in the G+C content of the descendant sequences with an overall decrease in the frequency of substitutions.

Experimental phylogeny of neutrally evolving DNA sequences generated by a bifurcate series of nested polymerase chain reactions.

A known phylogeny was generated using a four-step serial bifurcate PCR method. The ancestor sequence (SSU rDNA) evolved in vitro for 280 nested PCR cycles, and the resulting 15 ancestor and 16 terminal sequences (2,238 bp each) were determined. Parsimony, distance, and maximum likelihood analysis of the terminal sequences reconstructed the topology of the real phylogeny and branch lengths accurately. Divergence dates and ancestor sequences were estimated with very small error, particularly at the base of the phylogeny, mostly due to insertion and deletion changes. The substitution patterns along the known phylogeny are not described by reversible models, and accordingly, the probability substitution matrix, based on the observed substitutions from ancestor to terminal nodes along the known phylogeny, was calculated. This approach is an extension of previous studies using bacteriophage serial propagation, because here mutations were allowed to occur neutrally rather than by addition of a mutagenic agent, which produced biased mutational changes. These results provide for the first time biochemical experimental support for phylogenies, divergence date estimates, and an irreversible substitution model based on neutrally evolving DNA sequences. The substitution preferences observed here (A to G and T to C) are consistent with the high G+C content of the Thermus aquaticus genome. This suggests, at least in part, that the method here described, which explores the high Taq DNA polymerase error rate, simulates the evolution of a DNA segment in a thermophilic organism. These organisms include the bacterial rod T. aquaticus and several Archaea, and thus, the method and data set described here may well contribute new insights about the genome evolution of these organisms.