Compensatory evolution in response to a novel RNA polymerase: orthologous replacement of a central network gene

A bacteriophage genome was forced to evolve a new system of regulation by replacing its RNA polymerase (RNAP) gene, a central component of the phage developmental pathway, with that of a relative. The experiment used the obligate lytic phage T7 and the RNAP gene of phage T3. T7 RNAP uses 17 phage promoters, which are responsible for all middle and late gene expression, DNA replication, and progeny maturation, but the enzyme has known physical contacts with only 2 other phage proteins. T3 RNAP was supplied in trans by the bacterial host to a T7 genome lacking its own RNAP gene and the phage population was continually propagated on naive bacteria throughout the adaptation. Evolution of the T3 RNAP gene was thereby prevented, and selection was for the evolution of regulatory signals throughout the phage genome. T3 RNAP transcribes from T7 promoters only at low levels, but a single mutation in the promoter confers high expression, providing a ready mechanism for reevolution of gene expression in this system. When selected for rapid growth, fitness of the engineered phage evolved from a low of 5 doublings/h to 33 doublings/h, close to the expected maximum of 37 doublings/h. However, the experiment was terminated before it could be determined accurately that fitness had reached an obvious plateau, and it is not known whether further adaptation could have resulted in complete recovery of fitness. More than 30 mutations were observed in the evolved genome, but changes were found in only 9 of the 16 promoters, and several coding changes occurred in genes with no known contacts with the RNAP. Surprisingly, the T7 genome adapted to T3 RNAP also maintained high fitness when using T7 RNAP, suggesting that the extreme incompatibility of T7 elements with T3 RNAP is not an invariant property of divergence in these expression systems.     

Reversion and T Cell Escape Mutations Compensate the Fitness Loss of a CD8+ T Cell Escape Mutant in Their Cognate Transmitted/Founder Virus

Immune escape mutations that revert back to the consensus sequence frequently occur in newly HIV-1-infected individuals and have been thought to render the viruses more fit. However, their impact on viral fitness and their interaction with other immune escape mutations have not been evaluated in the background of their cognate transmitted/founder (T/F) viral genomes. To precisely determine the role of reversion mutations, we introduced reversion mutations alone or together with CD8⁺ T cell escape mutations in their unmodified cognate T/F viral genome and determined their impact on viral fitness in primary CD4⁺ T cells. Two reversion mutations, V247I and I64T, were identified in Gag and Tat, respectively, but neither had measurable effect on the fitness of their cognate T/F virus. The V247I and G248A mutations that were detected before and concurrently with the potent T cell escape mutation T242N, respectively, were selected by early T cell responses. The V247I or the G248A mutation alone partially restored the fitness loss caused by the T242N mutation. Together they could fully restore the fitness of the T242N mutant to the T/F level. These results demonstrate that the fitness loss caused by a T cell escape mutation could be compensated by preexisting or concurrent reversion and other T cell escape mutations. Our findings indicate that the overall viral fitness is modulated by the complex interplay among T cell escape, compensatory and reversion mutations to maintain the balance between immune escape and viral replication capacity.     

Experimental Evolution Yields Hundreds of Mutations in a Functional Viral Genome

Two lines of the bacteriophage T7 were grown to fix mutations indiscriminately, using a combination of population bottlenecks and mutagenesis. Complete genome sequences revealed 404 and 299 base substitutions in the two lines, the largest number characterized in functional microbial genomes so far. Missense substitutions outnumbered silent substitutions. Silent substitutions occurred at similar rates between essential and nonessential genes, but missense substitutions occurred at a higher rate in nonessential genes than in essential genes, as expected if they were less deleterious in the nonessential genes. Viral fitness declined during this protocol, and subsequent passaging of each mutated line in large population sizes restored some of the lost fitness. Substitution levels during these recoveries were less than 6% of those during the bottleneck phase, and only two changes during recovery were reversions of the original mutations. Exchanges of genomic fragments between the two recovered lines revealed that fitness effects of some substitutions were not additive—that interactions were accumulating which could lead to incompatibility between the diverged genomes. Based on these results, unprecedented high rates of nucleotide and functional divergence in viral genomes should be attainable experimentally by using repeated population bottlenecks at a high mutation rate interspersed with recovery. Present address (M. R. Badgett): Ambion, Inc., 2130 Woodward, Austin, TX 78744, USA Present address (D. Rokyta): Department of Biological Sciences, University of Idaho, Moscow, ID 83844, USA     

Effects of a single amino acid substitution within the p2 region of human immunodeficiency virus type 1 on packaging of spliced viral RNA

Human immunodeficiency virus type 1 encapsidates two copies of viral genomic RNA in the form of a dimer. The dimerization process initiates via a 6-nucleotide palindrome that constitutes the loop of a viral RNA stem-loop structure (i.e., stem loop 1 [SL1], also termed the dimerization initiation site [DIS]) located within the 5' untranslated region of the viral genome. We have now shown that deletion of the entire DIS sequence virtually eliminated viral replication but that this impairment was overcome by four second-site mutations located within the matrix (MA), capsid (CA), p2, and nucleocapsid (NC) regions of Gag. Interestingly, defective viral RNA dimerization caused by the DeltaDIS deletion was not significantly corrected by these compensatory mutations, which did, however, allow the mutated viruses to package wild-type levels of this DIS-deleted viral RNA while excluding spliced viral RNA from encapsidation. Further studies demonstrated that the compensatory mutation T12I located within p2, termed MP2, sufficed to prevent spliced viral RNA from being packaged into the DeltaDIS virus. Consistently, the DeltaDIS-MP2 virus displayed significantly higher levels of infectiousness than did the DeltaDIS virus. The importance of position T12 in p2 was further demonstrated by the identification of four point mutations,T12D, T12E, T12G, and T12P, that resulted in encapsidation of spliced viral RNA at significant levels. Taken together, our data demonstrate that selective packaging of viral genomic RNA is influenced by the MP2 mutation and that this represents a major mechanism for rescue of viruses containing the DeltaDIS deletion.     

Compensatory evolution for a gene deletion is not limited to its immediate functional network

Background  

Genetic disruption of an important phenotype should favor compensatory mutations that restore the phenotype. If the genetic basis of the phenotype is modular, with a network of interacting genes whose functions are specific to that phenotype, compensatory mutations are expected among the genes of the affected network. This perspective was tested in the bacteriophage T3 using a genome deleted of its DNA ligase gene, disrupting DNA metabolism.     

Molecular basis of fitness loss and fitness recovery in vesicular stomatitis virus

Viral populations subjected to repeated genetic bottleneck accumulate deleterious mutations in a process known as Muller's ratchet. Asexual viruses, such as vesicular stomatitis virus (VSV) can recover from Muller's ratchet by replication with large effective population sizes. However, mutants with a history of bottleneck transmissions often show decreased adaptability when compared to non-bottlenecked populations. We have generated a collection of bottlenecked mutants and allowed them to recover by large population passages. We have characterized fitness changes and the complete genomes of these strains. Mutations accumulated during the operation of Muller's ratchet led to the identification of two potential mutational hot spots in the VSV genome. As in other viral systems, transitions were more common than transversions. Both back mutation and compensatory mutations contributed to recovery, although a significant level of fitness increase was observed in nine of the 13 bottlenecked strains with no obvious changes in the consensus sequence. Additional replication of three strains resulted in the fixation of single point mutations. Only two mutations previously found in non-bottlenecked, high-fitness populations that had been adapting to the same environment were identified in the recovered strains.     

HIV-1 Replication Fitness of HLA-B*57/58:01 CTL Escape Variants Is Restored by the Accumulation of Compensatory Mutations in Gag

Expression of HLA-B*57 and the closely related HLA-B*58:01 are associated with prolonged survival after HIV-1 infection. However, large differences in disease course are observed among HLA-B*57/58:01 patients. Escape mutations in CTL epitopes restricted by these HLA alleles come at a fitness cost and particularly the T242N mutation in the TW10 CTL epitope in Gag has been demonstrated to decrease the viral replication capacity. Additional mutations within or flanking this CTL epitope can partially restore replication fitness of CTL escape variants. Five HLA-B*57/58:01 progressors and 5 HLA-B*57/58:01 long-term nonprogressors (LTNPs) were followed longitudinally and we studied which compensatory mutations were involved in the restoration of the viral fitness of variants that escaped from HLA-B*57/58:01-restricted CTL pressure. The Sequence Harmony algorithm was used to detect homology in amino acid composition by comparing longitudinal Gag sequences obtained from HIV-1 patients positive and negative for HLA-B*57/58:01 and from HLA-B*57/58:01 progressors and LTNPs. Although virus isolates from HLA-B*57/58:01 individuals contained multiple CTL escape mutations, these escape mutations were not associated with disease progression. In sequences from HLA-B*57/58:01 progressors, 5 additional mutations in Gag were observed: S126N, L215T, H219Q, M228I and N252H. The combination of these mutations restored the replication fitness of CTL escape HIV-1 variants. Furthermore, we observed a positive correlation between the number of escape and compensatory mutations in Gag and the replication fitness of biological HIV-1 variants isolated from HLA-B*57/58:01 patients, suggesting that the replication fitness of HLA-B*57/58:01 escape variants is restored by accumulation of compensatory mutations.     

Initial Fitness Recovery of HIV-1 Is Associated with Quasispecies Heterogeneity and Can Occur without Modifications in the Consensus Sequence

Background

Fitness recovery of HIV-1 “in vitro” was studied using viral clones that had their fitness decreased as a result of plaque-to-plaque passages.

Principal Findings

After ten large population passages, the viral populations showed an average increase of fitness, although with wide variations among clones. While 5 clones showed significant fitness increases, 3 clones showed increases that were only marginally significant (p<0.1), and 4 clones did not show any change. Fitness recovery was not accompanied by an increase in p24 production, but was associated with an increase in viral titer. Few mutations (an average of 2 mutations per genome) were detected in the consensus nucleotide sequence of the entire genome in all viral populations. Five of the populations did not fix any mutation, and three of them displayed marginally significant fitness increases, illustrating that fitness recovery can occur without detectable alterations of the consensus genomic sequence. The investigation of other possible viral factors associated with the initial steps of fitness recovery, showed that viral quasispecies heterogeneity increased between the initial clones and the passaged populations. A direct statistical correlation between viral heterogeneity and viral fitness was obtained.

Conclusions

Thus, the initial fitness recovery of debilitated HIV-1 clones was mediated by an increase in quasispecies heterogeneity. This observation, together with the invariance of the consensus sequence despite fitness increases demonstrates the relevance of quasispecies heterogeneity in the evolution of HIV-1 in cell culture.     

Progression to AIDS in South Africa is associated with both reverting and compensatory viral mutations

We lack the understanding of why HIV-infected individuals in South Africa progress to AIDS. We hypothesised that in end-stage disease there is a shifting dynamic between T cell imposed immunity and viral immune escape, which, through both compensatory and reverting viral mutations, results in increased viral fitness, elevated plasma viral loads and disease progression. We explored how T cell responses, viral adaptation and viral fitness inter-relate in South African cohorts recruited from Bloemfontein, the Free State (n = 278) and Durban, KwaZulu-Natal (n = 775). Immune responses were measured by γ-interferon ELISPOT assays. HLA-associated viral polymorphisms were determined using phylogenetically corrected techniques, and viral replication capacity (VRC) was measured by comparing the growth rate of gag-protease recombinant viruses against recombinant NL4-3 viruses. We report that in advanced disease (CD4 counts <100 cells/µl), T cell responses narrow, with a relative decline in Gag-directed responses (p<0.0001). This is associated with preserved selection pressure at specific viral amino acids (e.g., the T242N polymorphism within the HLA-B*57/5801 restricted TW10 epitope), but with reversion at other sites (e.g., the T186S polymorphism within the HLA-B*8101 restricted TL9 epitope), most notably in Gag and suggestive of “immune relaxation”. The median VRC from patients with CD4 counts <100 cells/µl was higher than from patients with CD4 counts ≥500 cells/µl (91.15% versus 85.19%, p = 0.0004), potentially explaining the rise in viral load associated with disease progression. Mutations at HIV Gag T186S and T242N reduced VRC, however, in advanced disease only the T242N mutants demonstrated increasing VRC, and were associated with compensatory mutations (p = 0.013). These data provide novel insights into the mechanisms of HIV disease progression in South Africa. Restoration of fitness correlates with loss of viral control in late disease, with evidence for both preserved and relaxed selection pressure across the HIV genome. Interventions that maintain viral fitness costs could potentially slow progression.     

Establishment of COS‐JC cells persistently producing archetype JC polyomavirus

JC polyomavirus (JCPyV) is the causative agent of progressive multifocal leukoencephalopathy (PML), a demyelinating disease of the central nervous system in immunocompromised patients. Archetype JCPyV circulates in the human population. There have been several reports of archetype JCPyV replication in cultured cells, in which propagation was not enough to produce high titers of archetype JCPyV. In this study, we carried out cultivation of the transfected cells with archetype JCPyV DNA MY for more than 2 months to establish COS‐7 cells (designated COS‐JC cells) persistently producing archetype JCPyV. Moreover, JCPyV derived from COS‐JC cells was characterized by analyzing the viral propagation, size of the viral genome, amount of viral DNA, production of viral protein, and structure of the non‐coding control region (NCCR). Southern blotting using a digoxigenin‐labeled JCPyV probe showed two different sizes of the JCPyV genome in COS‐JC cells. For molecular cloning, four of five clones showed a decrease in the size of complete JCPyV genome. Especially, clone No. 10 was generated the large deletion within the Large T antigen. On the other hand, the archetype structure of the NCCR was maintained in COS‐JC cells, although a few point mutations occurred. Quantitative PCR analysis of viral DNA in COS‐JC cells indicated that a high copy number of archetype JCPyV DNA was replicated in COS‐JC cells. These findings suggest that COS‐JC cells could efficiently propagate archetype JCPyV MY and offer a useful tool to study persistent infection of archetype JCPyV in a kidney‐derived system.

     

Identification of a herpesvirus Saimiri cis-acting DNA fragment that permits stable replication of episomes in transformed T cells.

Herpesvirus saimiri is a lymphotropic herpesvirus capable of immortalizing and transforming T cells both in vitro and in vivo. Immortalized and transformed T cells harbor several copies of the viral genome as a persisting genome. The mapping of the cis-acting genetic cis-acting segment (oriP) required for viral episomal maintenance is reported here. Viral DNA fragments that potentially contain oriP were cloned into a plasmid that contains the hygromycin resistance gene. After several round of subcloning followed by transfection, oriP was mapped to a 1.955-kb viral segment. This viral fragment permits stable plasmid replication without deletion or rearrangement as well as episomal maintenance without integration or recombination. The function of oriP depends on a trans-acting factor(s) encoded by the viral genome. The 1.955-kb viral segment includes a dyad symmetry region located between two small nuclear RNA genes and is located upstream of the dihydrofolate reductase gene homolog. Therefore, this oriP contains novel elements distinct from those of other DNA viruses.     

RNA replication kinetics, genetic polymorphism and selection in the case of the hepatitis C virus

We show in a simple theoretical quasispecies model that the replication dynamics of hepatitis C virus and a related model-system, the bovine viral diarrhoea virus, result in an effective reduction of RNA templates in infected cells. Viral fitness does not translate directly into RNA sequence replication efficiency, and hence the abundance of the viral master sequences diminishes over time. Our results suggest that genes not involved in RNA replication accumulate mutations over time because they do not undergo selection during this phase. The selection of viral RNA occurs not only during replication but also during the ensuing stages of the viral life cycle: (i) envelopment of viral RNA and (ii) successful infection of other cells, which also requires functionality of non-replicative genes. In particular, viral fitness requires the ability of the genome to encode structural proteins which do not encounter selective pressure during RNA replication. We conclude by discussing the potential value of antiviral drugs which inhibit selection on parts of the viral genome.     

Evolutionary Robustness of an Optimal Phenotype: Re-evolution of Lysis in a Bacteriophage Deleted for Its Lysin Gene

Optimality models are frequently used to create expectations about phenotypic evolution based on the fittest possible phenotype. However, they often ignore genetic details, which could confound these expectations. We experimentally analyzed the ability of organisms to evolve towards an optimum in an experimentally tractable system, lysis time in bacteriophage T7. T7 lysozyme helps lyse the host cell by degrading its cell wall at the end of infection, allowing viral escape to infect new hosts. Artificial deletion of lysozyme greatly reduced fitness and delayed lysis, but after evolution both phenotypes approached wild-type values. Phage with a lysis-deficient lysozyme evolved similarly. Several mutations were involved in adaptation, but most of the change in lysis timing and fitness increase was mediated by changes in gene 16, an internal virion protein not formerly considered to play a role in lysis. Its muralytic domain, which normally aids genome entry through the cell wall, evolved to cause phage release. Theoretical models suggest there is an optimal lysis time, and lysis more rapid or delayed than this optimum decreases fitness. Artificially constructed lines with very rapid lysis had lower fitness than wild-type T7, in accordance with the model. However, while a slow-lysing line also had lower fitness than wild-type, this low fitness resulted at least partly from genetic details that violated model assumptions.     

Resistance to extinction of low fitness virus subjected to plaque-to-plaque transfers: diversification by mutation clustering.

Plaque-to-plaque transfers of RNA viruses lead to accumulation of mutations and fitness decrease. To test whether continuing plaque-to-plaque transfers would lead to viral extinction, we have subjected several low fitness foot-and-mouth disease virus (FMDV) clones to up to 130 successive plaque transfers, and have analyzed the evolution of plaque titers and genomic nucleotide sequences. No case of viral extinction could be documented. Some low fitness clones that posses an internal poly(A) tract evaded extinction by modifying the length or base composition of the poly(A) tract. The comparison of entire genomic sequences of FMDV clones at increasing plaque transfer number revealed that mutations accumulated at a uniform rate, and that they were distributed unevenly along the genome. Clusters of mutations were identified at different genomic sites in two plaque transfer lineages. Mutation clustering appears to occur stochastically and could not be related to fixation of compensatory mutations. The results document resistance of viral clones to extinction, and suggest that mutation clustering may be a mechanism of genetic diversification of low fitness virus.     

Pyrosequencing reveals restricted patterns of CD8+ T cell escape-associated compensatory mutations in simian immunodeficiency virus

CD8+ T cells play a major role in the containment of human immunodeficiency virus (HIV) and simian immunodeficiency virus (SIV) replication. CD8+ T cell-driven variations in conserved regions under functional constraints result in diminished viral replicative capacity. While compensatory mutations outside an epitope can restore replicative capacity, the kinetics with which they arise remains unknown. Additionally, certain patterns of linked mutations associated with CD8+ T cell epitope escape in these highly conserved regions may lead to variable levels of viral fitness. Here, we used pyrosequencing to investigate the kinetics and patterns of mutations surrounding the Mamu-A1*00101-bound Gag(181-189)CM9 CD8+ T cell epitope. We obtained more than 400 reads for each sequencing time point, allowing us to confidently detect the emergence of viral variants bearing escape mutations with frequencies as low as 1% of the circulating virus. With this level of detail, we demonstrate that compensatory mutations generally arise concomitantly with Gag(181-189)CM9 escape mutations. We observed distinct patterns of linked flanking mutations, most of which were found downstream of Gag(181-189)CM9. Our data indicate that, whereas Gag(181-189)CM9 escape is much more complex that previously appreciated, it occurs in a coordinated fashion, with very specific patterns of flanking mutations required for immune evasion. This is the first detailed report of the ontogeny of compensatory mutations that allow CD8+ T cell epitope escape in infected individuals.     

Sequence and spacing requirements of a retrovirus integration site

Following infection, retroviruses insert a DNA copy of their RNA genome into the host cell genome. This integrative recombination reaction occurs at specific sites on the viral DNA: inverted repeat sequences near the termini of the linear DNA form of the viral genome. We have described elsewhere the generation and analysis of deletion mutations at one of the inverted repeat sequences in Moloney murine leukemia virus. We describe here the effects of insertion mutations made at this locus. Our results show that substantial sequence changes at the site of recombination can be tolerated, and that the spacing between the cleavage sites on the viral DNA can be expanded as well as contracted while still allowing efficient viral integration. After several rounds of virus replication, each of the insertion mutants gave rise to pseudorevertants with new alterations at the integration site.     

Deterministic, compensatory mutational events in the capsid of foot-and-mouth disease virus in response to the introduction of mutations found in viruses from persistent infections

The evolution of foot-and-mouth disease virus (FMDV) (biological clone C-S8c1) in persistently infected cells led to the emergence of a variant (R100) that displayed increased virulence, reduced stability, and other modified phenotypic traits. Some mutations fixed in the R100 genome involved a cluster of highly conserved residues around the capsid pores that participate in interactions with each other and/or between capsid protomers. We have investigated phenotypic and genotypic changes that occurred when these replacements were introduced into the C-S8c1 capsid. The C3007V and M3014L mutations exerted no effect on plaque size or viral yield during lytic infections, or on virion stability, but led to a reduction in biological fitness; the D3009A mutation caused drastic reductions in plaque size and viability. Remarkably, competition of the C3007V mutant with the nonmutated virus invariably resulted in the fixation of the D3009A mutation in the C3007V capsid. In turn, the presence of the D3009A mutation invariably led to the fixation of the M3014L mutation. In both cases, two individually disadvantageous mutations led, together, to an increase in fitness, as the double mutants outcompeted the nonmutated genotype. The higher fitness of C3007V/D3009A was related to a faster multiplication rate. These observations provide evidence for a chain of linked, compensatory mutational events in a defined region of the FMDV capsid. Furthermore, they indicate that the clustering of unique amino acid replacements in viruses from persistent infections may also occur in cytolytic infections in response to changes caused by previous mutations without an involvement of the new mutations in the adaptation to a different environment.     

The rate of compensatory mutation in the DNA bacteriophage phiX174

A compensatory mutation occurs when the fitness loss caused by one mutation is remedied by its epistatic interaction with a second mutation at a different site in the genome. This poorly understood biological phenomenon has important implications, not only for the evolutionary consequences of mutation, but also for the genetic complexity of adaptation. We have carried out the first direct experimental measurement of the average rate of compensatory mutation. An arbitrary selection of 21 missense substitutions with deleterious effects on fitness was introduced by site-directed mutagenesis into the bacteriophage phiX174. For each deleterious mutation, we evolved 8-16 replicate populations to determine the frequency at which a compensatory mutation, instead of the back mutation, was acquired to recover fitness. The overall frequency of compensatory mutation was approximately 70%. Deleterious mutations that were more severe were significantly more likely to be compensated for. Furthermore, experimental reversion of deleterious mutations revealed that compensatory mutations have deleterious effects in a wild-type background. A large diversity of intragenic compensatory mutations was identified from sequencing fitness-recovering genotypes. Subsequent analyses of intragenic mutation diversity revealed a significant degree of clustering around the deleterious mutation in the linear sequence and also within folded protein structures. Moreover, a likelihood analysis of mutation diversity predicts that, on average, a deleterious mutation can be compensated by about nine different intragenic compensatory mutations. We estimate that about half of all compensatory mutations are located extragenically in this organism.     

Independent mechanisms involved in suppression of the Moloney leukemia virus genome during differentiation of murine teratocarcinoma cells

Expression and DNA methylation of the Moloney murine leukemia virus (M-MuLV) genome were investigated in murine teratocarcinoma cells after virus infection. The newly acquired viral genome was devoid of methylation, yet its expression was repressed. The integrated viral genome in undifferentiated teratocarcinoma cells was methylated within 15 days after infection. Although 5-azacytidine decreased the level of DNA methylation, it did not activate M-MuLV in undifferentiated cells. Activation by 5-azacytidine occurred only in differentiated teratocarcinoma cells. Thus two independent mechanisms seem to regulate gene expression during the course of differentiation. The first mechanism operates in undifferentiated cells to block expression of M-MuLV and other exogeneously acquired viral genes, such as SV40 and polyoma virus, and does not depend on DNA methylation. The second mechanism relates only to differentiated cells and represses expression of genes in which DNA is methylated.     

Epistasis and the adaptability of an RNA virus

We have explored the patterns of fitness recovery in the vesicular stomatitis RNA virus. We show that, in our experimental setting, reversions to the wild-type genotype were rare and fitness recovery was at least partially driven by compensatory mutations. We compared compensatory adaptation for genotypes carrying (1) mutations with varying deleterious fitness effects, (2) one or two deleterious mutations, and (3) pairs of mutations showing differences in the strength and sign of epistasis. In all cases, we found that the rate of fitness recovery and the proportion of reversions were positively affected by population size. Additionally, we observed that mutations with large fitness effect were always compensated faster than mutations with small fitness effect. Similarly, compensatory evolution was faster for genotypes carrying a single deleterious mutation than for those carrying pairs of mutations. Finally, for genotypes carrying two deleterious mutations, we found evidence of a negative correlation between the epistastic effect and the rate of compensatory evolution.