Evolutionary reconstruction of a hairpin deleted from the genome of an RNA virus

Genome-regulated Assembly of a ssRNA Virus May Also Prepare It for Infection

Many single-stranded, positive-sense RNA viruses regulate assembly of their infectious virions by forming multiple, cognate coat protein (CP)-genome contacts at sites termed Packaging Signals (PSs). We have determined the secondary structures of the bacteriophage MS2 ssRNA genome (gRNA) frozen in defined states using constraints from X-ray synchrotron footprinting (XRF). Comparison of the footprints from phage and transcript confirms the presence of multiple PSs in contact with CP dimers in the former. This is also true for a virus-like particle (VLP) assembled around the gRNA in vitro in the absence of the single-copy Maturation Protein (MP) found in phage. Since PS folds are present at many sites across gRNA transcripts, it appears that this genome has evolved to facilitate this mechanism of assembly regulation. There are striking differences between the gRNA-CP contacts seen in phage and the VLP, suggesting that the latter are inappropriate surrogates for aspects of phage structure/function. Roughly 50% of potential PS sites in the gRNA are not in contact with the protein shell of phage. However, many of these sit adjacent to, albeit not in contact with, PS-binding sites on CP dimers. We hypothesize that these act as PSs transiently during assembly but subsequently dissociate. Combining the XRF data with PS locations from an asymmetric cryo-EM reconstruction suggests that the genome positions of such dissociations are non-random and may facilitate infection. The loss of many PS-CP interactions towards the 3' end of the gRNA would allow this part of the genome to transit more easily through the narrow basal body of the pilus extruding machinery. This is the known first step in phage infection. In addition, each PS-CP dissociation event leaves the protein partner trapped in a non-lowest free-energy conformation. This destabilizes the protein shell which must disassemble during infection, further facilitating this stage of the life-cycle.

Computational based design and tracking of synthetic variants of Porcine circovirus reveal relations between silent genomic information and viral fitness

Viral genomes not only code the protein content, but also include silent, overlapping codes which are important to the regulation of the viral life cycle and affect its evolution. Due to the high density of these codes, their non-modular nature and the complex intracellular processes they encode, the ability of current approaches to decipher them is very limited. We describe the first computational-experimental pipeline for studying the effects of viral silent and non-silent information on its fitness. The pipeline was implemented to study the Porcine Circovirus type 2 (PCV2), the shortest known eukaryotic virus, and includes the following steps: (1) Based on the analyses of 2100 variants of PCV, suspected silent codes were inferred. (2) Five hundred variants of the PCV2 were designed to include various ‘smart’ silent mutations. (3) Using state of the art synthetic biology approaches, the genomes of these five hundred variants were generated. (4) Competition experiments between the variants were performed in Porcine kidney-15 (PK15) cell-lines. (5) The variant titers were analyzed based on novel next-generation sequencing (NGS) experiments. (6) The features related to the titer of the variants were inferred and their analyses enabled detection of various novel silent functional sequence and structural motifs. Furthermore, we demonstrate that 50 of the silent variants exhibit higher fitness than the wildtype in the analyzed conditions.

Emergency Services of Viral RNAs: Repair and Remodeling

Reproduction of RNA viruses is typically error-prone due to the infidelity of their replicative machinery and the usual lack of proofreading mechanisms. The error rates may be close to those that kill the virus. Consequently, populations of RNA viruses are represented by heterogeneous sets of genomes with various levels of fitness. This is especially consequential when viruses encounter various bottlenecks and new infections are initiated by a single or few deviating genomes. Nevertheless, RNA viruses are able to maintain their identity by conservation of major functional elements. This conservatism stems from genetic robustness or mutational tolerance, which is largely due to the functional degeneracy of many protein and RNA elements as well as to negative selection. Another relevant mechanism is the capacity to restore fitness after genetic damages, also based on replicative infidelity. Conversely, error-prone replication is a major tool that ensures viral evolvability. The potential for changes in debilitated genomes is much higher in small populations, because in the absence of stronger competitors low-fit genomes have a choice of various trajectories to wander along fitness landscapes. Thus, low-fit populations are inherently unstable, and it may be said that to run ahead it is useful to stumble. In this report, focusing on picornaviruses and also considering data from other RNA viruses, we review the biological relevance and mechanisms of various alterations of viral RNA genomes as well as pathways and mechanisms of rehabilitation after loss of fitness. The relationships among mutational robustness, resilience, and evolvability of viral RNA genomes are discussed.

Viral attenuation by engineered protein fragmentation

A possible but untested method of viral attenuation is protein fragmentation, engineering wild-type proteins as two or more peptides that self-assemble after translation. Here, the bacteriophage T7 was engineered to encode its essential RNA polymerase as two peptides. Initial fitness was profoundly suppressed. Subjecting the engineered virus to over 100 generations of adaptation by serial transfer resulted in a large fitness increase, still remaining below that of evolved wild-type. The fitness increase was accompanied by three substitutions in the fragmented peptides as well as six mutations in other parts of the genome, but the fragmentation was retained. This study thereby demonstrates the feasibility of using gene fragmentation as a possibly permanent method of attenuation, but the initial fitness of the engineered genome may be a poor measure of its fitness on extended adaptation.

Dynamics of molecular evolution in RNA virus populations depend on sudden versus gradual environmental change

Understanding the dynamics of molecular adaptation is a fundamental goal of evolutionary biology. While adaptation to constant environments has been well characterized, the effects of environmental complexity remain seldom studied. One simple but understudied factor is the rate of environmental change. Here we used experimental evolution with RNA viruses to investigate whether evolutionary dynamics varied based on the rate of environmental turnover. We used whole-genome next-generation sequencing to characterize evolutionary dynamics in virus populations adapting to a sudden versus gradual shift onto a novel host cell type. In support of theoretical models, we found that when populations evolved in response to a sudden environmental change, mutations of large beneficial effect tended to fix early, followed by mutations of smaller beneficial effect; as predicted, this pattern broke down in response to a gradual environmental change. Early mutational steps were highly parallel across replicate populations in both treatments. The fixation of single mutations was less common than sweeps of associated "cohorts" of mutations, and this pattern intensified when the environment changed gradually. Additionally, clonal interference appeared stronger in response to a gradual change. Our results suggest that the rate of environmental change is an important determinant of evolutionary dynamics in asexual populations.

Synonymous mutations reduce genome compactness in icosahedral ssRNA viruses

Recent studies have shown that single-stranded (ss) viral RNAs fold into more compact structures than random RNA sequences with similar chemical composition and identical length. Based on this comparison, it has been suggested that wild-type viral RNA may have evolved to be atypically compact so as to aid its encapsidation and assist the viral assembly process. To further explore the compactness selection hypothesis, we systematically compare the predicted sizes of >100 wild-type viral sequences with those of their mutants, which are evolved in silico and subject to a number of known evolutionary constraints. In particular, we enforce mutation synonynimity, preserve the codon-bias, and leave untranslated regions intact. It is found that progressive accumulation of these restricted mutations still suffices to completely erase the characteristic compactness imprint of the viral RNA genomes, making them in this respect physically indistinguishable from randomly shuffled RNAs. This shows that maintaining the physical compactness of the genome is indeed a primary factor among ssRNA viruses' evolutionary constraints, contributing also to the evidence that synonymous mutations in viral ssRNA genomes are not strictly neutral.

Rapid evolution of in vivo-selected sequences and structures replacing 20% of a subviral RNA

The 356 nt noncoding satellite RNA C (satC) of Turnip crinkle virus (TCV) is composed of 5' sequences from a second TCV satRNA (satD) and 3' sequences derived from TCV. SHAPE structure mapping revealed that 76 nt in the poorly-characterized satD-derived region form an extended hairpin (H2). Pools of satC in which H2 was replaced with 76, 38, or 19 random nt were co-inoculated with TCV helper virus onto plants and satC fitness assessed using in vivo functional selection (SELEX). The most functional progeny satCs, including one as fit as wild-type, contained a 38-39 nt H2 region that adopted a hairpin structure and exhibited an increased ratio of dimeric to monomeric molecules. Some progeny of satC with H2 deleted featured a duplication of 38 nt, partially rebuilding the deletion. Therefore, H2 can be replaced by a 38-39 nt hairpin, sufficient for overall structural stability of the 5' end of satC.

A recombinant RNA bacteriophage system to identify functionally important nucleotides in a self-cleaving ribozyme

Background

RNA bacteriophages like Qbeta and MS2 are well known for their high mutation rate, short infection cycle and strong selection against foreign inserts. The hammerhead ribozyme (HHRz) is a small self-cleaving RNA molecule whose active residues have previously been identified by mutational analysis of each individual base. Here the functionally important bases of HHRz were determined in a single screening experiment by inserting the HHRz into the genome of MS2.

Findings

The minimal HHRz of satellite Tobacco ringspot virus was cloned into the genome of RNA bacteriophage MS2. Sequence analysis of the surviving phages revealed that the majority had acquired single base-substitutions that apparently inactivated the HHRz. The positions of these substitutions exactly matched that of the previously determined core residues of the HHRz.

Conclusions

Natural selection against a ribozyme in the genome of MS2 can be used to quickly identify nucleotides required for self-cleavage.

Influence of genome-scale RNA structure disruption on the replication of murine norovirus--similar replication kinetics in cell culture but attenuation of viral fitness in vivo

Mechanisms by which certain RNA viruses, such as hepatitis C virus, establish persistent infections and cause chronic disease are of fundamental importance in viral pathogenesis. Mammalian positive-stranded RNA viruses establishing persistence typically possess genome-scale ordered RNA secondary structure (GORS) in their genomes. Murine norovirus (MNV) persists in immunocompetent mice and provides an experimental model to functionally characterize GORS. Substitution mutants were constructed with coding sequences in NS3/4- and NS6/7-coding regions replaced with sequences with identical coding and (di-)nucleotide composition but disrupted RNA secondary structure (F1, F2, F1/F2 mutants). Mutants replicated with similar kinetics to wild-type (WT) MNV3 in RAW264.7 cells and primary macrophages, exhibited similar (highly restricted) induction and susceptibility to interferon-coupled cellular responses and equal replication fitness by serial passaging of co-cultures. In vivo, both WT and F1/F2 mutant viruses persistently infected mice, although F1, F2 and F1/F2 mutant viruses were rapidly eliminated 1-7 days post-inoculation in competition experiments with WT. F1/F2 mutants recovered from tissues at 9 months showed higher synonymous substitution rates than WT and nucleotide substitutions that potentially restored of RNA secondary structure. GORS plays no role in basic replication of MNV but potentially contributes to viral fitness and persistence in vivo.

Inactivation and tailing during UV254 disinfection of viruses: contributions of viral aggregation, light shielding within viral aggregates, and recombination

UV disinfection of viruses frequently leads to tailing after an initial exponential decay. Aggregation, light shielding, recombination, or resistant virus subpopulations have been proposed as explanations; however, none of these options has been conclusively demonstrated. This study investigates how aggregation affects virus inactivation by UV(254) in general, and the tailing phenomenon in particular. Bacteriophage MS2 was aggregated by lowering the solution pH before UV(254) disinfection. Aggregates were redispersed prior to enumeration to obtain the remaining fraction of individual infectious viruses. Results showed that initial inactivation kinetics were similar for viruses incorporated in aggregates (up to 1000 nm in radius) and dispersed viruses; however, aggregated viruses started to tail more readily than dispersed ones. Neither light shielding, nor the presence of resistant subpopulations could account for the tailing. Instead, tailing was consistent with recombination arising from the simultaneous infection of the host by several impaired viruses. We argue that UV(254) treatment of aggregates permanently fused a fraction of viruses, which increased the likelihood of multiple infection of a host cell and ultimately enabled the production of infective viruses via recombination.

Genomics of bacterial and archaeal viruses: dynamics within the prokaryotic virosphere

Prokaryotes, bacteria and archaea, are the most abundant cellular organisms among those sharing the planet Earth with human beings (among others). However, numerous ecological studies have revealed that it is actually prokaryotic viruses that predominate on our planet and outnumber their hosts by at least an order of magnitude. An understanding of how this viral domain is organized and what are the mechanisms governing its evolution is therefore of great interest and importance. The vast majority of characterized prokaryotic viruses belong to the order Caudovirales, double-stranded DNA (dsDNA) bacteriophages with tails. Consequently, these viruses have been studied (and reviewed) extensively from both genomic and functional perspectives. However, albeit numerous, tailed phages represent only a minor fraction of the prokaryotic virus diversity. Therefore, the knowledge which has been generated for this viral system does not offer a comprehensive view of the prokaryotic virosphere. In this review, we discuss all families of bacterial and archaeal viruses that contain more than one characterized member and for which evolutionary conclusions can be attempted by use of comparative genomic analysis. We focus on the molecular mechanisms of their genome evolution as well as on the relationships between different viral groups and plasmids. It becomes clear that evolutionary mechanisms shaping the genomes of prokaryotic viruses vary between different families and depend on the type of the nucleic acid, characteristics of the virion structure, as well as the mode of the life cycle. We also point out that horizontal gene transfer is not equally prevalent in different virus families and is not uniformly unrestricted for diverse viral functions.

The fitness effects of synonymous mutations in DNA and RNA viruses

Despite being silent with respect to protein sequence, synonymous nucleotide substitutions can be targeted by natural selection directly at the DNA or RNA level. However, there has been no systematic assessment of how frequent this type of selection is. Here, we have constructed 53 single random synonymous substitution mutants of the bacteriophages Qβ and ΦX174 by site-directed mutagenesis and assayed their fitness. Analysis of this mutant collection and of previous studies undertaken with a variety of single-stranded (ss) viruses demonstrates that selection at synonymous sites is stronger in RNA viruses than in DNA viruses. We estimate that this type of selection contributes approximately 18% of the overall mutational fitness effects in ssRNA viruses under our assay conditions and that random synonymous substitutions have a 5% chance of being lethal to the virus, whereas in ssDNA viruses, these figures drop to 1.4% and 0%, respectively. In contrast, the effects of nonsynonymous substitutions appear to be similar in ssRNA and ssDNA viruses.

Experimental evolution, loss-of-function mutations, and "the first rule of adaptive evolution"

Adaptive evolution can cause a species to gain, lose, or modify a function; therefore, it is of basic interest to determine whether any of these modes dominates the evolutionary process under particular circumstances. Because mutation occurs at the molecular level, it is necessary to examine the molecular changes produced by the underlying mutation in order to assess whether a given adaptation is best considered as a gain, loss, or modification of function. Although that was once impossible, the advance of molecular biology in the past half century has made it feasible. In this paper, I review molecular changes underlying some adaptations, with a particular emphasis on evolutionary experiments with microbes conducted over the past four decades. I show that by far the most common adaptive changes seen in those examples are due to the loss or modification of a pre-existing molecular function, and I discuss the possible reasons for the prominence of such mutations.

HIV-1 evolution: frustrating therapies, but disclosing molecular mechanisms

Replication of HIV-1 under selective pressure frequently results in the evolution of virus variants that replicate more efficiently under the applied conditions. For example, in patients on antiretroviral therapy, such evolution can result in variants that are resistant to the HIV-1 inhibitors, thus frustrating the therapy. On the other hand, virus evolution can help us to understand the molecular mechanisms that underlie HIV-1 replication. For example, evolution of a defective virus mutant can result in variants that overcome the introduced defect by restoration of the original sequence or by the introduction of additional mutations in the viral genome. Analysis of the evolution pathway can reveal the requirements of the element under study and help to understand its function. Analysis of the escape routes may generate new insight in the viral life cycle and result in the identification of unexpected biological mechanisms. We have developed in vitro HIV-1 evolution into a systematic research tool that allows the study of different aspects of the viral replication cycle. We will briefly review this method of forced virus evolution and provide several examples that illustrate the power of this approach.

How RNA viruses maintain their genome integrity

RNA genomes are vulnerable to corruption by a range of activities, including inaccurate replication by the error-prone replicase, damage from environmental factors, and attack by nucleases and other RNA-modifying enzymes that comprise the cellular intrinsic or innate immune response. Damage to coding regions and loss of critical cis-acting signals inevitably impair genome fitness; as a consequence, RNA viruses have evolved a variety of mechanisms to protect their genome integrity. These include mechanisms to promote replicase fidelity, recombination activities that allow exchange of sequences between different RNA templates, and mechanisms to repair the genome termini. In this article, we review examples of these processes from a range of RNA viruses to showcase the diverse approaches that viruses have evolved to maintain their genome sequence integrity, focusing first on mechanisms that viruses use to protect their entire genome, and then concentrating on mechanisms that allow protection of the genome termini, which are especially vulnerable. In addition, we discuss examples in which it might be beneficial for a virus to 'lose' its genomic termini and reduce its replication efficiency.

Virus evolution as a tool to study HIV-1 biology

Mutational analysis of the viral genome is frequently used to study the role of sequence or structural elements in HIV-1 replication. Many laboratories that use this approach have occasionally come across revertant viruses that overcome an introduced defect either by restoration of the original sequence or by the introduction of additional mutations in the viral genome. Similarly, replication of a wild type virus under selective pressure, due to the presence of inhibitors or due to specific culture settings, may result in the appearance of evolved variants that replicate more efficiently under the applied conditions. We have developed in vitro HIV-1 evolution from an anecdotal event to a systematic research tool to study different aspects of the viral replication cycle. In this manuscript, we will briefly review the method of forced virus evolution to study HIV-1 biology and provide several examples that illustrate the power of this method, as it frequently yielded interesting and unexpected information about the mechanism of virus replication.

Second site escape of a T20-dependent HIV-1 variant by a single amino acid change in the CD4 binding region of the envelope glycoprotein

Background

We previously described the selection of a T20-dependent human immunodeficiency virus type-1 (HIV-1) variant in a patient on T20 therapy. The fusion inhibitor T20 targets the viral envelope (Env) protein by blocking a conformational switch that is critical for viral entry into the host cell. T20-dependent viral entry is the result of 2 mutations in Env (GIA-SKY), creating a protein that undergoes a premature conformational switch, and the presence of T20 prevents this premature switch and rescues viral entry. In the present study, we performed 6 independent evolution experiments with the T20-dependent HIV-1 variant in the absence of T20, with the aim to identify second site compensatory changes, which may provide new mechanistic insights into Env function and the T20-dependence mechanism.

Results

Escape variants with improved replication capacity appeared within 42 days in 5 evolution cultures. Strikingly, 3 cultures revealed the same single amino acid change in the CD4 binding region of Env (glycine at position 431 substituted for arginine: G431R). This mutation was sufficient to abolish the T20-dependence phenotype and restore viral replication in the absence of T20. The GIA-SKY-G431R escape variant produces an Env protein that exhibits reduced syncytia formation and reduced cell-cell fusion activity. The escape variant was more sensitive to an antibody acting on an early gp41 intermediate, suggesting that the G431R mutation helps preserve a pre-fusion Env conformation, similar to T20 action. The escape variant was also less sensitive to soluble CD4, suggesting a reduced CD4 receptor affinity.

Conclusion

The forced evolution experiments indicate that the premature conformational switch of the T20-dependent HIV-1 Env variant (GIA-SKY) can be corrected by a second site mutation in Env (GIA-SKY-G431R) that affects the interaction with the CD4 receptor.

Structural constraints and mutational bias in the evolutionary restoration of a severe deletion in RNA phage MS2

A 4-nucleotide (nt) deletion was made in the 36-nt-long intercistronic region separating the coat and replicase genes of the single-stranded RNA phage MS2. This region is the focus of several RNA structures conferring high fitness. One such element is the operator hairpin, which, in the course of infection, will bind a coat-protein dimer, thereby precluding further replicase synthesis and initiating encapsidation. Another structure is a long-distance base pairing (MJ) controlling replicase expression. The 4-nt deletion does not directly affect the operator hairpin but it disrupts the MJ pairing. Its main effect, however, is a frame shift in the overlapping lysis gene. This gene starts in the upstream coat gene, runs through the 36-nt-long intercistronic region, and ends in the downstream replicase cistron. Here we report and interpret the spectrum of solutions that emerges when the crippled phage is evolved. Four different solutions were obtained by sequencing 40 plaques. Three had cured the frame shift in the lysis gene by inserting one nt in the loop of the operator hairpin causing its inactivation. Yet these low-fitness revertants could further improve themselves when evolved. The inactivated operator was replaced by a substitute and thereafter these revertants found several ways to restore control over the replicase gene. To allow for the evolutionary enrichment of low-probability but high-fitness revertants, we passaged lysate samples before plating. Revertants obtained in this way also restored the frame shift, but not at the expense of the operator. By taking larger and larger lysates samples for such bulk evolution, ever higher-fitness and lower-frequency revertants surfaced. Only one made it back to wild type. As a rule, however, revertants moved further and further away from the wild-type sequence because restorative mutations are, in the majority of cases, selected for their capacity to improve the phenotype by optimizing one of several potential alternative RNA foldings that emerge as a result of the initial deletion. This illustrates the role of structural constraints which limit the path of subsequent restorative mutations.

What is a quasispecies?

The concept of the quasispecies as a society formed from a clone of an asexually reproducing organism is reviewed. A broad spectrum of mutants is generated that compete one with another. Eventually a steady state is formed where each mutant type is represented according to its fitness and its formation by mutation. This quasispecies has a defined wild type sequence, which is the weighted average of all genotypes present. The quasispecies concept has been shown to affect the pathway of evolution and has been studied on RNA viruses which have a particularly high mutation rate. They (and possibly the majority of other species) operate close to the error threshold that allows maximum exploration of sequence space while conserving the information content of the genotype. The consequences of the quasispecies concept for the new 'evolutionary technology' are discussed.

The error threshold

Basic principles underlying the population dynamics of bacteria and viruses are presented, with emphasis on RNA viruses. Concepts reviewed here include fitness, mutant generation, competition, selection, sequence space and the theoretical origins of quasispecies. A "wild-type" virus is no longer viewed as the fittest type, but as the center of gravity of a population landscape. A consequence of quasispecies is the existence of an error threshold for selective competence. The error threshold has a justification in quasispecies theory and lends itself to exploration of strategies to drive viral populations over the error threshold, the central theme of this volume.

Adaptation to the deleterious effects of antimicrobial drug resistance mutations by compensatory evolution

Compensatory mutations, due to their ability to mask the deleterious effects of another mutation, are important for the adaptation and evolution of most organisms. Resistance to antibiotics, antivirals, antifungals, herbicides and insecticides is usually associated with a fitness cost. As a result of compensatory evolution, the initial fitness costs conferred by resistance mutations (or other deleterious mutations) can often be rapidly and efficiently reduced. Such compensatory evolution is potentially of importance for (i) the long-term persistence of drug resistance, (ii) reducing the rate of fitness loss associated with the accumulation of deleterious mutations in small asexual populations, and (iii) the evolution of complexity of cellular processes.

Compensatory adaptation to the deleterious effect of antibiotic resistance in Salmonella typhimurium

Most chromosomal mutations that cause antibiotic resistance impose fitness costs on the bacteria. This biological cost can often be reduced by compensatory mutations. In Salmonella typhimurium, the nucleotide substitution AAA42 --> AAC in the rpsL gene confers resistance to streptomycin. The resulting amino acid substitution (K42N) in ribosomal protein S12 causes an increased rate of ribosomal proofreading and, as a result, the rate of protein synthesis, bacterial growth and virulence are decreased. Eighty-one independent lineages of the low-fitness, K42N mutant were evolved in the absence of antibiotic to ameliorate the costs. From the rate of fixation of compensated mutants and their fitness, the rate of compensatory mutations was estimated to be > or = 10-7 per cell per generation. The size of the population bottleneck during evolution affected fitness of the adapted mutants: a larger bottleneck resulted in higher average fitness. Only four of the evolved lineages contained streptomycin-sensitive revertants. The remaining 77 lineages contained mutants that were still fully streptomycin resistant, had retained the original resistance mutation and also acquired compensatory mutations. Most of the compensatory mutations, resulting in at least 35 different amino acid substitutions, were novel single-nucleotide substitutions in the rpsD, rpsE, rpsL or rplS genes encoding the ribosomal proteins S4, S5, S12 and L19 respectively. Our results show that the deleterious effects of a resistance mutation can be compensated by an unexpected variety of mutations.

Nucleotide sequence of a ssRNA phage from Acinetobacter: kinship to coliphages

The complete nucleotide sequence of ssRNA phage AP205 propagating in Acinetobacter species is reported. The RNA has three large ORFs, which code for the following homologues of the RNA coliphage proteins: the maturation, coat and replicase proteins. Their gene order is the same as that in coliphages. RNA coliphages or Leviviridae fall into two genera: the alloleviviruses, like Q(beta), which have a coat read-through protein, and the leviviruses, like MS2, which do not have this coat protein extension. AP205 has no read-through protein and may therefore be classified as a levivirus. A major digression from the known leviviruses is the apparent absence of a lysis gene in AP205 at the usual position, overlapping the coat and replicase proteins. Instead, two small ORFs are present at the 5' terminus, preceding the maturation gene. One of these might encode a lysis protein. The other is of unknown function. Other new features concern the 3'-terminal sequence. In all ssRNA coliphages, there are always three cytosine residues at the 3' end, but in AP205, there is only a single terminal cytosine. Distantly related viruses, like AP205 and the coliphages, do not have significant sequence identity; yet, important secondary structural features of the RNA are conserved. This is shown here for the 3' UTR and the replicase-operator hairpin. Interestingly, although AP205 has the genetic map of a levivirus, its 3' UTR has the length and RNA secondary structure of an allolevivirus. Sharing features with both MS2 and Q(beta) suggests that, in an evolutionary sense, AP205 should be placed between Q(beta) and MS2. A phylogenetic tree for the ssRNA phages is presented.

Experimental genomic evolution: extensive compensation for loss of DNA ligase activity in a virus

Deletion of the viral ligase gene drastically reduced the fitness of bacteriophage T7 on a ligase-deficient host. Viral evolution recovered much of this fitness during long-term passage, but the final fitness remained below that of the intact virus. Compensatory changes occurred chiefly in genes involved in DNA metabolism: the viral endonuclease, helicase, and DNA polymerase. Two other compensatory changes of unknown function also occurred. Using a method to distinguish compensatory mutations from other beneficial mutations, five additional substitutions from the recovery were shown to enhance adaptation to culture conditions and were not compensatory for the deletion. In contrast to the few previous studies of viral recovery from deletions, the compensatory changes in T7 did not restore the deletion or duplicate major regions of the genome. The ability of this deleted genome to recover much of the lost fitness via mutations in its remaining genes reveals a considerable evolutionary potential to modify the interactions of its elements in maintaining an essential set of functions.

Forced retroevolution of an RNA bacteriophage

The operator hairpin ahead of the replicase gene in RNA bacteriophage MS2 contains overlapping signals for binding the coat protein and ribosomes. Coat protein binding inhibits further translation of the gene and forms the first step in capsid formation. The hairpin sequence was partially randomized to assess the importance of this structure element for the bacteriophage and to monitor alternative solutions that would evolve on the passaging of mutant phages. The evolutionary reconstruction of the operator failed in the majority of mutants. Instead, a poor imitation developed containing only some of the recognition signals for the coat protein. Three mutants were of particular interest in that they contained double nonsense codons in the lysis reading frame that runs through the operator hairpin. The simultaneous reversion of two stop codons into sense codons has a very low probability of occurring. Therefore the phage solved the problem by deleting the nonsense signals and, in fact, the complete operator, except for the initiation codon of the replicase gene. Several revertants were isolated with activities ranging from 1% to 20% of wild type. The operator, long thought to be a critical regulator, now appears to be a dispensable element. In addition, the results indicate how RNA viruses can be forced to step back to an attenuated form.

A second-site mutation that restores replication of a Tat-defective human immunodeficiency virus

We previously constructed a large set of mutants of the human immunodeficiency virus type 1 (HIV-1) regulatory protein Tat with conservative amino acid substitutions in the activation domain. These Tat variants were analyzed in the context of the infectious virus, and several mutants were found to be defective for replication. In an attempt to obtain second-site suppressor mutations that could provide information on the Tat protein structure, some of the replication-impaired viruses were used as a parent for the isolation of revertant viruses with improved replication capacity. Sequence analysis of revertant viruses frequently revealed changes within the tat gene, most often first-site reversions either to the wild-type amino acid or to related amino acids that restore, at least partially, the Tat function and virus replication. Of 30 revertant cultures, we identified only one second-site suppressor mutation. The inactive Y26A mutant yielded the second-site suppressor mutation Y47N that partially restored trans-activation activity and virus replication. Surprisingly, when the suppressor mutation was introduced in the wild-type Tat background, it also improved the trans-activation function of this protein about twofold. We conclude that the gain of function measured for the Y47N change is not specific for the Y26A mutant, arguing against a direct interaction of Tat amino acids 26 and 47 in the three-dimensional fold of this protein. Other revertant viruses did not contain any additional Tat changes, and some viruses revealed putative second-site Tat mutations that did not significantly improve Tat function and virus replication. We reason that these mutations were introduced by chance through founder effects or by linkage to suppressor mutations elsewhere in the virus genome. In conclusion, the forced evolution of mutant HIV-1 genomes, which is an efficient approach for the analysis of RNA regulatory motifs, seems less suited for the analysis of the structure of this small transcription factor, although protein variants with interesting properties can be generated.

Dynamic competition between alternative structures in viroid RNAs simulated by an RNA folding algorithm

The folding pathways of viroid RNAs were studied using computer simulations by the genetic algorithm for RNA folding. The folding simulations were performed for PSTVd RNAs of both polarities, using the wild-type sequence and some previously known mutants with suggested changes in the stable or metastable structures. It is shown that metastable multihairpin foldings in the minus strand replicative intermediates are established due to the specific folding pathway that ensures the absence of the most stable rod-like structure. Simulations of the PSTVd minus strand folding during transcription reveal a metastable hairpin, formed in the left terminal domain region of the PSTVd. Despite high sequence variability, this hairpin is conserved in all known large viroids of both subgroups of PSTVd type, and is presumably necessary to guide the folding of the HPII hairpin which is functional in the minus strand. The folding simulations are able to demonstrate the changes in the balance between metastable and stable structures in mutant PSTVd RNAs. The stable rod-like structure of the circular viroid (+) RNA is also folded via a dynamic folding pathway. Furthermore, the simulations show that intermediate steps in the forced evolution of a shortened PSTVd replicon may be reconstructed by a mechanistic model of different folding pathway requirements in plus- and minus-strand RNAs. Thus the formation of viroid RNA structure strongly depends on dynamics of competition between alternative RNA structures. This also suggests that the replication efficiency of viroid sequences may be estimated by a simulation of the folding process.

Rescue of the RNA phage genome from RNase III cleavage

The secondary structure of the RNA from the single-stranded RNA bacteriophages, like MS2 and Qb, has evolved to serve a variety of functions such as controlling gene expression, exposing binding sites for the replicase and capsid proteins, allowing strand separation and so forth. On the other hand, all of these foldings have to perform in bacterial cells in which various RNA splitting enzymes are present. We therefore examined whether phage RNA structure is under selective pressure by host RNases. Here we show this to be true for RNase III. A fully double-stranded hairpin of 17 bp, which is an RNase III target, was inserted into a non-coding region of the MS2 RNA genome. In an RNase III-host these phages survived but in wild-type bacteria they did not. Here the stem underwent Darwinian evolution to a structure that was no longer a substrate for RNase III. This was achieved in three different ways: (i) the perfect stem was maintained but shortened by removing all or most of the insert; (ii) the stem acquired suppressor mutations that replaced Watson-Crick base pairs by mismatches; (iii) the stem acquired small deletions or insertions that created bulges. These insertions consist of short stretches of non-templated A or U residues. Their origin is ascribed to polyadenylation at the site of the RNase III cut (in the + or - strand) either by Escherichia coli poly(A) polymerase or by idling MS2 replicase.

Altered 3'-terminal RNA structure in phage Qbeta adapted to host factor-less Escherichia coli

The RNA phage Qbeta requires for the replication of its genome an RNA binding protein called Qbeta host factor or Hfq protein. Our previous results suggested that this protein mediates the access of replicase to the 3'-end of the Qbeta plus strand RNA. Here we report the results of an evolutionary experiment in which phage Qbeta was adapted to an Escherichia coli Q13 host strain with an inactivated host factor (hfq) gene. This strain initially produced phage at a titer approximately 10,000-fold lower than the wild-type strain and with minute plaque morphology, but after 12 growth cycles, phage titer and plaque size had evolved to levels near those of the wild-type host. RNAs isolated from adapted Qbeta mutants were efficient templates for replicase without host factor in vitro. Electron microscopy showed that mutant RNAs, in contrast to wild-type RNA, efficiently interacted with replicase at the 3'-end in the absence of host factor. The same set of four mutations in the 3'-terminal third of the genome was found in several independently evolved phage clones. One mutation disrupts the base pairing of the 3'-terminal CCCOH sequence, suggesting that the host factor stimulates activity of the wild-type RNA template by melting out its 3'-end.