Positive selection of synonymous mutations in vesicular stomatitis virus

Prevailing evolutionary forces are typically deduced from the pattern of differences in synonymous and non-synonymous mutations, under the assumption of neutrality in the absence of amino acid change. We determined the complete sequence of ten vesicular stomatitis virus populations evolving under positive selection. A significant number of the mutations occurred independently in two or more strains, a process known as parallel evolution, and a substantial fraction of the parallel mutations were silent. Parallel evolution was also identified in non-coding regions. These results indicate that silent mutations can significantly contribute to adaptation in RNA viruses, and relative frequencies of synonymous and non-synonymous substitutions may not be useful to resolve their evolutionary history.     

The human immunodeficiency virus type 1 encapsidation site is a multipartite RNA element composed of functional hairpin structures.

We analyzed the leader region of human immunodeficiency virus type 1 (HIV-1) RNA to decipher the nature of the cis-acting E/psi element required for encapsidation of viral RNA into virus particles. Our data indicate that, for RNA encapsidation, there are at least two functional subregions in the leader region. One subregion is located at a position immediately proximal to the major splice donor, and the second is located between the splice donor and the beginning of the gag gene. This suggests that at least two discrete cis-acting elements are recognition signals for encapsidation. To determine whether specific putative RNA secondary structures serve as the signal(s) for encapsidation, we constructed primary base substitution mutations that would be expected to destabilize these potential structures and second-site compensatory mutations that would restore secondary structure. Analysis of these mutants allowed the identification of two discrete hairpins that facilitate RNA encapsidation in vivo. Thus, the HIV-1 E/psi region is a multipartite element composed of specific and functional RNA secondary structures. Compensation of the primary mutations by the second-site mutations could not be attained in trans. This indicates that interstrand base pairing between these two stem regions within the hairpins does not appear to be the basis for HIV-1 RNA dimer formation. Comparison of the hypothetical RNA secondary structures from 10 replication-competent HIV-1 strains suggests that a subset of the hydrogen-bonded base pairs within the stems of the hairpins is likely to be required for function in cis.     

Identification of three feline immunodeficiency virus (FIV) env gene subtypes and comparison of the FIV and human immunodeficiency virus type 1 evolutionary patterns.

Feline immunodeficiency virus (FIV) is a lentivirus associated with AIDS-like illnesses in cats. As such, FIV appears to be a feline analog of human immunodeficiency virus (HIV). A hallmark of HIV infection is the large degree of viral genetic diversity that can develop within an infected individual and the even greater and continually increasing level of diversity among virus isolates from different individuals. Our goal in this study was to determine patterns of FIV genetic diversity by focusing on a 684-nucleotide region encompassing variable regions V3, V4, and V5 of the FIV env gene in order to establish parallels and distinctions between FIV and HIV type 1 (HIV-1). Our data demonstrate that, like HIV-1, FIV can be separated into distinct envelope sequence subtypes (three are described here). Similar to that found for HIV-1, the pairwise sequence divergence within an FIV subtype ranged from 2.5 to 15.0%, whereas that between subtypes ranged from 17.8 to 26.2%. However, the high number of synonymous nucleotide changes among FIV V3 to V5 env sequences may also include a significant number of back mutations and suggests that the evolutionary distances among FIV subtypes are underestimated. Although only a few subtype B viruses were available for examination, the pattern of diversity between the FIV A and B subtypes was found to be significantly distinct; subtype B sequences had proportionally fewer mutations that changed amino acids, compared with silent changes, suggesting a more advanced state of adaptation to the host. No similar distinction was evident for HIV-1 subtypes. The diversity of FIV genomes within individual infected cats was found to be as high as 3.7% yet twofold lower than that within HIV-1-infected people over a comparable region of the env gene. Despite these differences, significant parallels between patterns of FIV evolution and HIV-1 evolution exist, indicating that a wide array of potentially divergent virus challenges need to be considered in FIV vaccine and pathogenesis studies.     

Founder virus population related to route of virus transmission: a determinant of intrahost human immunodeficiency virus type 1 evolution?

We and others have shown that in individual human immunodeficiency virus type 1 (HIV-1) infection, the adaptive evolution of HIV-1 is influenced by host immune competence. In this study, we tested the hypothesis that in addition to selective forces operating within the host, transmission bottlenecks have an impact on HIV-1 intrahost evolution. Therefore, we studied the intrahost evolution of the V3 region of the external glycoprotein gp120 of HIV-1 during the 3- and 5-year periods following seroconversion after parenteral versus sexual (male-to-male) transmission in 41 participants of the Amsterdam prospective cohorts of homosexual men (n = 31) and intravenous drug users (IVDUs; n = 10) who were AIDS free and had comparable numbers of CD4+ cells. We observed that HIV-1 strains in homosexual men accumulated over 5 years more nonsynonymous substitutions within the V3 loop than HIV-1 strains in IVDUs as a result of lower rates of nonsynonymous evolution in both the initial 3-year period from seroconversion and the following 2-year period as well as a larger proportion of nonsynonymous back substitutions in IVDUs. The mean numbers of synonymous substitutions did not differ between the two risk groups. Since HIV-1 strains in IVDUs could be distinguished from the viruses of homosexual men based on several nucleotide substitutions of which the most conserved is a synonymous substitution at the tip of the V3 loop (GGC pattern), we studied whether the founder virus population itself has an impact on the intrahost evolution of HIV-1. The mean number of nonsynonymous substitutions accumulated over 5 years within the V3 loop was lower in 10 IVDUs infected by the HIV-1 strains with the GGC signature than in 4 IVDUs infected by HIV-1 strains lacking this pattern, while the mean numbers of synonymous substitutions were similar in the two groups.     

A phylogenetic method for detecting positive epistasis in gene sequences and its application to RNA virus evolution

RNA virus genomes are compact, often containing multiple overlapping reading frames and functional secondary structure. Consequently, it is thought that evolutionary interactions between nucleotide sites are commonplace in the genomes of these infectious agents. However, the role of epistasis in natural populations of RNA viruses remains unclear. To investigate the pervasiveness of epistasis in RNA viruses, we used a parsimony-based computational method to identify pairs of co-occurring mutations along phylogenies of 177 RNA virus genes. This analysis revealed widespread evidence for positive epistatic interactions at both synonymous and nonsynonymous nucleotide sites and in both clonal and recombining viruses, with the majority of these interactions spanning very short sequence regions. These findings have important implications for understanding the key aspects of RNA virus evolution, including the dynamics of adaptation. Additionally, many comparative analyses that utilize the phylogenetic relationships among gene sequences assume that mutations represent independent, uncorrelated events. Our results show that this assumption may often be invalid.     

A Short Sequence Motif in the 5′ Leader of the HIV-1 Genome Modulates Extended RNA Dimer Formation and Virus Replication

The 5′ leader of the HIV-1 RNA genome encodes signals that control various steps in the replication cycle, including the dimerization initiation signal (DIS) that triggers RNA dimerization. The DIS folds a hairpin structure with a palindromic sequence in the loop that allows RNA dimerization via intermolecular kissing loop (KL) base pairing. The KL dimer can be stabilized by including the DIS stem nucleotides in the intermolecular base pairing, forming an extended dimer (ED). The role of the ED RNA dimer in HIV-1 replication has hardly been addressed because of technical challenges. We analyzed a set of leader mutants with a stabilized DIS hairpin for in vitro RNA dimerization and virus replication in T cells. In agreement with previous observations, DIS hairpin stability modulated KL and ED dimerization. An unexpected previous finding was that mutation of three nucleotides immediately upstream of the DIS hairpin significantly reduced in vitro ED formation. In this study, we tested such mutants in vivo for the importance of the ED in HIV-1 biology. Mutants with a stabilized DIS hairpin replicated less efficiently than WT HIV-1. This defect was most severe when the upstream sequence motif was altered. Virus evolution experiments with the defective mutants yielded fast replicating HIV-1 variants with second site mutations that (partially) restored the WT hairpin stability. Characterization of the mutant and revertant RNA molecules and the corresponding viruses confirmed the correlation between in vitro ED RNA dimer formation and efficient virus replication, thus indicating that the ED structure is important for HIV-1 replication.     

Codon conservation in the influenza A virus genome defines RNA packaging signals

Genome segmentation facilitates reassortment and rapid evolution of influenza A virus. However, segmentation complicates particle assembly as virions must contain all eight vRNA species to be infectious. Specific packaging signals exist that extend into the coding regions of most if not all segments, but these RNA motifs are poorly defined. We measured codon variability in a large dataset of sequences to identify areas of low nucleotide sequence variation independent of amino acid conservation in each segment. Most clusters of codons showing very little synonymous variation were located at segment termini, consistent with previous experimental data mapping packaging signals. Certain internal regions of conservation, most notably in the PA gene, may however signify previously unidentified functions in the virus genome. To experimentally test the bioinformatics analysis, we introduced synonymous mutations into conserved codons within known packaging signals and measured incorporation of the mutant segment into virus particles. Surprisingly, in most cases, single nucleotide changes dramatically reduced segment packaging. Thus our analysis identifies cis-acting sequences in the influenza virus genome at the nucleotide level. Furthermore, we propose that strain-specific differences exist in certain packaging signals, most notably the haemagglutinin gene; this finding has major implications for the evolution of pandemic viruses.     

Dimer initiation signal of human immunodeficiency virus type 1: its role in partner selection during RNA copackaging and its effects on recombination

Frequent human immunodeficiency virus type 1 (HIV-1) recombination occurs during DNA synthesis when portions of the two copackaged RNAs are used as templates to generate a hybrid DNA copy. Therefore, the frequency of copackaging of genomic RNAs from two different viruses (heterozygous virion formation) affects the generation of genotypically different recombinants. We hypothesized that the selection of copackaged RNA partners is largely determined by Watson-Crick pairing at the dimer initiation signal (DIS), a 6-nucleotide palindromic sequence at the terminal loop of stem-loop 1 (SL1). To test our hypothesis, we examined whether heterozygous virion formation could be encouraged by manipulation of the DIS. Three pairs of viruses were generated with compensatory DIS mutations, designed so that perfect DIS base pairing could only occur between RNAs derived from different viruses, not between RNAs from the same virus. We observed that vector pairs with compensatory DIS mutations had an almost twofold increase in recombination rates compared with wild-type viruses. These data suggest that heterozygous virion formation was enhanced in viruses with compensatory DIS mutations (from 50% to more than 90% in some viral pairings). The role of the SL1 stem in heterozygous virion formation was also tested; our results indicated that the intermolecular base pairing of the stem sequences does not affect RNA partner selection. In summary, our results demonstrate that the Watson-Crick pairing of the DIS is a major determinant in the selection of the copackaged RNA partner, and altering the base pairing of the DIS can change the proportion of heterozygous viruses in a viral population. These results also strongly support the hypothesis that HIV-1 RNA dimers are formed prior to encapsidation.     

Full-length RNA structure prediction of the HIV-1 genome reveals a conserved core domain

A distance constrained secondary structural model of the ≈10 kb RNA genome of the HIV-1 has been predicted but higher-order structures, involving long distance interactions, are currently unknown. We present the first global RNA secondary structure model for the HIV-1 genome, which integrates both comparative structure analysis and information from experimental data in a full-length prediction without distance constraints. Besides recovering known structural elements, we predict several novel structural elements that are conserved in HIV-1 evolution. Our results also indicate that the structure of the HIV-1 genome is highly variable in most regions, with a limited number of stable and conserved RNA secondary structures. Most interesting, a set of long distance interactions form a core organizing structure (COS) that organize the genome into three major structural domains. Despite overlapping protein-coding regions the COS is supported by a particular high frequency of compensatory base changes, suggesting functional importance for this element. This new structural element potentially organizes the whole genome into three major domains protruding from a conserved core structure with potential roles in replication and evolution for the virus.     

Antisense-mediated inhibition of human immunodeficiency virus (HIV) replication by use of an HIV type 1-based vector results in severely attenuated mutants incapable of developing resistance

We have constructed a human immunodeficiency virus type 1 (HIV-1)-based lentiviral vector expressing a 937-base antisense sequence against the HIV-1 envelope gene. Transduction of CD4(+) T lymphocytes with this vector results in expression of the therapeutic antisense sequence and subsequent inhibition of productive HIV-1 replication. In this report, we examined the effect of antisense-mediated suppression on the potential development of virus escape mutants using a permissive T-cell line cultured under conditions that over serial passages specifically allowed for generation and amplification of mutants selected for by antisense pressure. In the resulting virus clones, we found a significant increase in the number of deletions at the envelope target region (91% compared to 27.5% in wild-type HIV). Deletions were most often greater than 1 kb in length. These data demonstrate for the first time that during antisense-mediated suppression of HIV, mutants develop as a direct result of selective pressure on the HIV genomic RNA. Interestingly, in clones where deletions were not observed, there was a high rate of A-G transitions in mutants at the antisense target region but not outside this region, which is consistent with those mutations that are predicted as a result of antisense-mediated modification of double-stranded RNA by the enzyme double-stranded RNA-specific adenosine deaminase. These clones were not found to be escape mutants, as their replicative ability was severely attenuated, and they did not replicate in the presence of vector.     

Long-term experimental evolution of HIV-1 reveals effects of environment and mutational history

An often-returning question for not only HIV-1, but also other organisms, is how predictable evolutionary paths are. The environment, mutational history, and random processes can all impact the exact evolutionary paths, but to which extent these factors contribute to the evolutionary dynamics of a particular system is an open question. Especially in a virus like HIV-1, with a large mutation rate and large population sizes, evolution is expected to be highly predictable if the impact of environment and history is low, and evolution is not neutral. We investigated the effect of environment and mutational history by analyzing sequences from a long-term evolution experiment, in which HIV-1 was passaged on 2 different cell types in 8 independent evolutionary lines and 8 derived lines, 4 of which involved a switch of the environment. The experiments lasted for 240–300 passages, corresponding to approximately 400–600 generations or almost 3 years. The sequences show signs of extensive parallel evolution—the majority of mutations that are shared between independent lines appear in both cell types, but we also find that both environment and mutational history significantly impact the evolutionary paths. We conclude that HIV-1 evolution is robust to small changes in the environment, similar to a transmission event in the absence of an immune response or drug pressure. We also find that the fitness landscape of HIV-1 is largely smooth, although we find some evidence for both positive and negative epistatic interactions between mutations.

Analysis of the longest evolutionary experiment with HIV-1 to-date reveals continuous viral adaptation over several years. The authors quantify the environment-specific mutations that arise and determine the fraction of mutations that co-occur with significantly different frequencies than expected by chance.     

Human immunodeficiency virus type 1 genetic evolution in children with different rates of development of disease.

The rate of development of disease varies considerably among human immunodeficiency virus type 1 (HIV-1)-infected children. The reasons for these observed differences are not clearly understood but most probably depend on the dynamic interplay between the HIV-1 quasispecies virus population and the immune constraints imposed by the host. To study the relationship between disease progression and genetic diversity, we analyzed the evolution of viral sequences within six perinatally infected children by examining proviral sequences spanning the C2 through V5 regions of the viral envelope gene by PCR of blood samples obtained at sequential visits. PCR product DNAs from four sample time points per child were cloned, and 10 to 13 clones from each sample were sequenced. Greater genetic distances relative to the time of infection were found for children with low virion-associated RNA burdens and slow progression to disease relative to those found for children with high virion-associated RNA burdens and rapid progression to disease. The greater branch lengths observed in the phylogenetic reconstructions correlated with a higher accumulation rate of nonsynonymous base substitutions per potential nonsynonymous site, consistent with positive selection for change rather than a difference in replication kinetics. Viral sequences from children with slow progression to disease also showed a tendency to form clusters that associated with different sampling times. These progressive shifts in the viral population were not found in viral sequences from children with rapid progression to disease. Therefore, despite the HIV-1 quasispecies being a diverse, rapidly evolving, and competing population of genetic variants, different rates of genetic evolution could be found under different selective constraints. These data suggest that the evolutionary dynamics exhibited by the HIV-1 quasispecies virus populations are compatible with a Darwinian system evolving under the constraints of natural selection.