The Fitness Effects of Random Mutations in Single-Stranded DNA and RNA Bacteriophages

Mutational fitness effects can be measured with relatively high accuracy in viruses due to their small genome size, which facilitates full-length sequencing and genetic manipulation. Previous work has shown that animal and plant RNA viruses are very sensitive to mutation. Here, we characterize mutational fitness effects in single-stranded (ss) DNA and ssRNA bacterial viruses. First, we performed a mutation-accumulation experiment in which we subjected three ssDNA (ΦX174, G4, F1) and three ssRNA phages (Qβ, MS2, and SP) to plaque-to-plaque transfers and chemical mutagenesis. Genome sequencing and growth assays indicated that the average fitness effect of the accumulated mutations was similar in the two groups. Second, we used site-directed mutagenesis to obtain 45 clones of ΦX174 and 42 clones of Qβ carrying random single-nucleotide substitutions and assayed them for fitness. In ΦX174, 20% of such mutations were lethal, whereas viable ones reduced fitness by 13% on average. In Qβ, these figures were 29% and 10%, respectively. It seems therefore that high mutational sensitivity is a general property of viruses with small genomes, including those infecting animals, plants, and bacteria. Mutational fitness effects are important for understanding processes of fitness decline, but also of neutral evolution and adaptation. As such, these findings can contribute to explain the evolution of ssDNA and ssRNA viruses.     

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.     

The Effect of Mutation and Selection on Codon Adaptation in Escherichia coli Bacteriophage

Studying phage codon adaptation is important not only for understanding the process of translation elongation, but also for reengineering phages for medical and industrial purposes. To evaluate the effect of mutation and selection on phage codon usage, we developed an index to measure selection imposed by host translation machinery, based on the difference in codon usage between all host genes and highly expressed host genes. We developed linear and nonlinear models to estimate the C→T mutation bias in different phage lineages and to evaluate the relative effect of mutation and host selection on phage codon usage. C→T-biased mutations occur more frequently in single-stranded DNA (ssDNA) phages than in double-stranded DNA (dsDNA) phages and affect not only synonymous codon usage, but also nonsynonymous substitutions at second codon positions, especially in ssDNA phages. The host translation machinery affects codon adaptation in both dsDNA and ssDNA phages, with a stronger effect on dsDNA phages than on ssDNA phages. Strand asymmetry with the associated local variation in mutation bias can significantly interfere with codon adaptation in both dsDNA and ssDNA phages.     

Stone age diseases and modern AIDS

Background

Recent reports have indicated that single-stranded DNA (ssDNA) viruses in the taxonomic families Geminiviridae, Parvoviridae and Anellovirus may be evolving at rates of ~10^-4 substitutions per site per year (subs/site/year). These evolution rates are similar to those of RNA viruses and are surprisingly high given that ssDNA virus replication involves host DNA polymerases with fidelities approximately 10 000 times greater than those of error-prone viral RNA polymerases. Although high ssDNA virus evolution rates were first suggested in evolution experiments involving the geminivirus maize streak virus (MSV), the evolution rate of this virus has never been accurately measured. Also, questions regarding both the mechanistic basis and adaptive value of high geminivirus mutation rates remain unanswered.

Results

We determined the short-term evolution rate of MSV using full genome analysis of virus populations initiated from cloned genomes. Three wild type viruses and three defective artificial chimaeric viruses were maintained in planta for up to five years and displayed evolution rates of between 7.4 × 10^-4 and 7.9 × 10^-4 subs/site/year.

Conclusion

These MSV evolution rates are within the ranges observed for other ssDNA viruses and RNA viruses. Although no obvious evidence of positive selection was detected, the uneven distribution of mutations within the defective virus genomes suggests that some of the changes may have been adaptive. We also observed inter-strand nucleotide substitution imbalances that are consistent with a recent proposal that high mutation rates in geminiviruses (and possibly ssDNA viruses in general) may be due to mutagenic processes acting specifically on ssDNA molecules.     

Classification of viral zoonosis through receptor pattern analysis

Background  

Viral zoonosis, the transmission of a virus from its primary vertebrate reservoir species to humans, requires ubiquitous cellular proteins known as receptor proteins. Zoonosis can occur not only through direct transmission from vertebrates to humans, but also through intermediate reservoirs or other environmental factors. Viruses can be categorized according to genotype (ssDNA, dsDNA, ssRNA and dsRNA viruses). Among them, the RNA viruses exhibit particularly high mutation rates and are especially problematic for this reason. Most zoonotic viruses are RNA viruses that change their envelope proteins to facilitate binding to various receptors of host species. In this study, we sought to predict zoonotic propensity through the analysis of receptor characteristics. We hypothesized that the major barrier to interspecies virus transmission is that receptor sequences vary among species--in other words, that the specific amino acid sequence of the receptor determines the ability of the viral envelope protein to attach to the cell.     

Mechanisms for RNA Capture by ssDNA Viruses: Grand Theft RNA

Viruses contain three common types of packaged genomes; double-stranded DNA (dsDNA), RNA (mostly single and occasionally double stranded) and single-stranded DNA (ssDNA). There are relatively straightforward explanations for the prevalence of viruses with dsDNA and RNA genomes, but the evolutionary basis for the apparent success of ssDNA viruses is less clear. The recent discovery of four ssDNA virus genomes that appear to have been formed by recombination between co-infecting RNA and ssDNA viruses, together with the high mutation rate of ssDNA viruses provide possible explanations. RNA–DNA recombination allows ssDNA viruses to access much broader sequence space than through nucleotide substitution and DNA–DNA recombination alone. Multiple non-exclusive mechanisms, all due to the unique replication of ssDNA viruses, are proposed for this unusual RNA capture. RNA capture provides an explanation for the evolutionary success of the ssDNA viruses and may help elucidate the mystery of integrated RNA viruses in viral and cellular DNA genomes.     

Comparison of Five Bacteriophages as Models for Viral Aerosol Studies

Bacteriophages are perceived to be good models for the study of airborne viruses because they are safe to use, some of them display structural features similar to those of human and animal viruses, and they are relatively easy to produce in large quantities. Yet, only a few studies have investigated them as models. It has previously been demonstrated that aerosolization, environmental conditions, and sampling conditions affect viral infectivity, but viral infectivity is virus dependent. Thus, several virus models are likely needed to study their general behavior in aerosols. The aim of this study was to compare the effects of aerosolization and sampling on the infectivity of five tail-less bacteriophages and two pathogenic viruses: MS2 (a single-stranded RNA [ssRNA] phage of the Leviviridae family), Φ6 (a segmented double-stranded RNA [dsRNA] phage of the Cystoviridae family), ΦX174 (a single-stranded DNA [ssDNA] phage of the Microviridae family), PM2 (a double-stranded DNA [dsDNA] phage of the Corticoviridae family), PR772 (a dsDNA phage of the Tectiviridae family), human influenza A virus H1N1 (an ssRNA virus of the Orthomyxoviridae family), and the poultry virus Newcastle disease virus (NDV; an ssRNA virus of the Paramyxoviridae family). Three nebulizers and two nebulization salt buffers (with or without organic fluid) were tested, as were two aerosol sampling devices, a liquid cyclone (SKC BioSampler) and a dry cyclone (National Institute for Occupational Safety and Health two-stage cyclone bioaerosol sampler). The presence of viruses in collected air samples was detected by culture and quantitative PCR (qPCR). Our results showed that these selected five phages behave differently when aerosolized and sampled. RNA phage MS2 and ssDNA phage ΦX174 were the most resistant to aerosolization and sampling. The presence of organic fluid in the nebulization buffer protected phages PR772 and Φ6 throughout the aerosolization and sampling with dry cyclones. In this experimental setup, the behavior of the influenza virus resembled that of phages PR772 and Φ6, while the behavior of NDV was closer to that of phages MS2 and ΦX174. These results provide critical information for the selection of appropriate phage models to mimic the behavior of specific human and animal viruses in aerosols.     

Upper-limit mutation rate estimation for a plant RNA virus

It is generally accepted that mutation rates of RNA viruses are inherently high due to the lack of proofreading mechanisms. However, direct estimates of mutation rate are surprisingly scarce, in particular for plant viruses. Here, based on the analysis of in vivo mutation frequencies in tobacco etch virus, we calculate an upper-bound mutation rate estimation of 3×10⁻⁵ per site and per round of replication; a value which turns out to be undistinguishable from the methodological error. Nonetheless, the value is barely on the lower side of the range accepted for RNA viruses, although in good agreement with the only direct estimate obtained for other plant viruses. These observations suggest that, perhaps, differences in the selective pressures operating during plant virus evolution may have driven their mutation rates towards values lower than those characteristic of other RNA viruses infecting bacteria or animals.     

Single-stranded adenine-rich DNA and RNA retain structural characteristics of their respective double-stranded conformations and show directional differences in stacking pattern

The structural preorganization of isosequential ssDNA and ssRNA hexamers d/r(GAAAAC)(1) [J. Am. Chem. Soc. 2003, 125, 9948] have been investigated by NMR and molecular dynamics simulations. Analysis of the nuclear Overhauser effect spectrometry (NOESY) footprints in the aqueous solution has shown that there is a substantial population of ordered right-handed helical structure in both hexameric single-stranded DNA and RNA, which are reminiscent of their respective right-handed helical duplex form, despite the fact these single-stranded molecules are devoid of any intermolecular hydrogen bonds. The NMR-constrained molecular dynamics (1.5 ns) derived geometries of the adenine-adenine overlaps at each dinucleotide step of the hexameric ssDNA (1a) and ssRNA (1b) show that the relatively electron-rich imidazole stacks above the electron-deficient pyrimidine in 5' to 3' direction in ssDNA (1a) while, in contradistinction, the pyrimidine stacks above the imidazole in the 5' to 3' direction in ssRNA (1b). This also means that the pi-frame of the 5'-pyrimidine can interact with the relatively positively charged imino and amino protons in the 3' direction in ssRNA and in the 5' direction in ssDNA, thereby stabilizing the twist and slide observed in the stacked oligonucleotides. The differently preferred stacking geometries in ssDNA and ssRNA have direct physicochemical implications for self-assembly and pK(a) modulation by the nearest-neighbor interactions, as well as for the dangling-end stabilization effects and imino-proton reactivity.     

Resistance of Aerosolized Bacterial Viruses to Relative Humidity and Temperature

The use of aerosolized bacteriophages as surrogates for hazardous viruses might simplify and accelerate the discovery of links between viral components and their persistence in the airborne state under diverse environmental conditions. In this study, four structurally distinct lytic phages, MS2 (single-stranded RNA [ssRNA]), ϕ6 (double-stranded RNA [dsRNA]), ϕX174 (single-stranded DNA [ssDNA]), and PR772 (double-stranded DNA [dsDNA]), were nebulized into a rotating chamber and exposed to various levels of relative humidity (RH) and temperature as well as to germicidal UV radiation. The aerosolized viral particles were allowed to remain airborne for up to 14 h before being sampled for analysis by plaque assays and quantitative PCRs. Phages ϕ6 and MS2 were the most resistant at low levels of relative humidity, while ϕX174 was more resistant at 80% RH. Phage ϕ6 lost its infectivity immediately after exposure to 30°C and 80% RH. The infectivity of all tested phages rapidly declined as a function of the exposure time to UVC radiation, phage MS2 being the most resistant. Taken altogether, our data indicate that these aerosolized phages behave differently under various environmental conditions and highlight the necessity of carefully selecting viral simulants in bioaerosol studies.     

Viruses and thyroiditis: an update

Background

Despite the demonstration that geminiviruses, like many other single stranded DNA viruses, are evolving at rates similar to those of RNA viruses, a recent study has suggested that grass-infecting species in the genus Mastrevirus may have co-diverged with their hosts over millions of years. This "co-divergence hypothesis" requires that long-term mastrevirus substitution rates be at least 100,000-fold lower than their basal mutation rates and 10,000-fold lower than their observable short-term substitution rates. The credibility of this hypothesis, therefore, hinges on the testable claim that negative selection during mastrevirus evolution is so potent that it effectively purges 99.999% of all mutations that occur.

Results

We have conducted long-term evolution experiments lasting between 6 and 32 years, where we have determined substitution rates of between 2 and 3 × 10^-4 substitutions/site/year for the mastreviruses Maize streak virus (MSV) and Sugarcane streak Réunion virus (SSRV). We further show that mutation biases are similar for different geminivirus genera, suggesting that mutational processes that drive high basal mutation rates are conserved across the family. Rather than displaying signs of extremely severe negative selection as implied by the co-divergence hypothesis, our evolution experiments indicate that MSV and SSRV are predominantly evolving under neutral genetic drift.

Conclusion

The absence of strong negative selection signals within our evolution experiments and the uniformly high geminivirus substitution rates that we and others have reported suggest that mastreviruses cannot have co-diverged with their hosts.     

Quasispecies Theory and the Behavior of RNA Viruses

A large number of medically important viruses, including HIV, hepatitis C virus, and influenza, have RNA genomes. These viruses replicate with extremely high mutation rates and exhibit significant genetic diversity. This diversity allows a viral population to rapidly adapt to dynamic environments and evolve resistance to vaccines and antiviral drugs. For the last 30 years, quasispecies theory has provided a population-based framework for understanding RNA viral evolution. A quasispecies is a cloud of diverse variants that are genetically linked through mutation, interact cooperatively on a functional level, and collectively contribute to the characteristics of the population. Many predictions of quasispecies theory run counter to traditional views of microbial behavior and evolution and have profound implications for our understanding of viral disease. Here, we discuss basic principles of quasispecies theory and describe its relevance for our understanding of viral fitness, virulence, and antiviral therapeutic strategy.     

Patterns of evolution and host gene mimicry in influenza and other RNA viruses

It is well known that the dinucleotide CpG is under-represented in the genomic DNA of many vertebrates. This is commonly thought to be due to the methylation of cytosine residues in this dinucleotide and the corresponding high rate of deamination of 5-methycytosine, which lowers the frequency of this dinucleotide in DNA. Surprisingly, many single-stranded RNA viruses that replicate in these vertebrate hosts also have a very low presence of CpG dinucleotides in their genomes. Viruses are obligate intracellular parasites and the evolution of a virus is inexorably linked to the nature and fate of its host. One therefore expects that virus and host genomes should have common features. In this work, we compare evolutionary patterns in the genomes of ssRNA viruses and their hosts. In particular, we have analyzed dinucleotide patterns and found that the same patterns are pervasively over- or under-represented in many RNA viruses and their hosts suggesting that many RNA viruses evolve by mimicking some of the features of their host's genes (DNA) and likely also their corresponding mRNAs. When a virus crosses a species barrier into a different host, the pressure to replicate, survive and adapt, leaves a footprint in dinucleotide frequencies. For instance, since human genes seem to be under higher pressure to eliminate CpG dinucleotide motifs than avian genes, this pressure might be reflected in the genomes of human viruses (DNA and RNA viruses) when compared to those of the same viruses replicating in avian hosts. To test this idea we have analyzed the evolution of the influenza virus since 1918. We find that the influenza A virus, which originated from an avian reservoir and has been replicating in humans over many generations, evolves in a direction strongly selected to reduce the frequency of CpG dinucleotides in its genome. Consistent with this observation, we find that the influenza B virus, which has spent much more time in the human population, has adapted to its human host and exhibits an extremely low CpG dinucleotide content. We believe that these observations directly show that the evolution of RNA viral genomes can be shaped by pressures observed in the host genome. As a possible explanation, we suggest that the strong selection pressures acting on these RNA viruses are most likely related to the innate immune response and to nucleotide motifs in the host DNA and RNAs.     

Effects of antibiotic resistance alleles on bacterial evolutionary responses to viral parasites

Antibiotic resistance has wide-ranging effects on bacterial phenotypes and evolution. However, the influence of antibiotic resistance on bacterial responses to parasitic viruses remains unclear, despite the ubiquity of such viruses in nature and current interest in therapeutic applications. We experimentally investigated this by exposing various Escherichia coli genotypes, including eight antibiotic-resistant genotypes and a mutator, to different viruses (lytic bacteriophages). Across 960 populations, we measured changes in population density and sensitivity to viruses, and tested whether variation among bacterial genotypes was explained by their relative growth in the absence of parasites, or mutation rate towards phage resistance measured by fluctuation tests for each phage. We found that antibiotic resistance had relatively weak effects on adaptation to phages, although some antibiotic-resistance alleles impeded the evolution of resistance to phages via growth costs. By contrast, a mutator allele, often found in antibiotic-resistant lineages in pathogenic populations, had a relatively large positive effect on phage-resistance evolution and population density under parasitism. This suggests costs of antibiotic resistance may modify the outcome of phage therapy against pathogenic populations previously exposed to antibiotics, but the effects of any co-occurring mutator alleles are likely to be stronger.     

Escape from Prisoner's Dilemma in RNA phage phi6

We previously examined competitive interactions among viruses by allowing the RNA phage phi6 to evolve at high and low multiplicities of infection (ratio of infecting viruses to bacterial cells). Derived high-multiplicity phages were competitively advantaged relative to their ancestors during coinfection, but their fixation caused population fitness to decline. These data conform to the evolution of lowered fitness in a population of defectors, as expected from the Prisoner's Dilemma of game theory. However, the generality of this result is unknown; the evolution of viruses at other multiplicities may alter the fitness payoffs associated with conflicting strategies of cooperation and defection. Here we examine the change in matrix variables by propagating the ancestor under strictly clonal conditions, allowing cooperation the chance to evolve. In competitions involving derived cooperators and their selfish counterparts, our data reveal a new outcome where the two strategies are predicted to coexist in a mixed polymorphism. Thus, we demonstrate that the payoff matrix is not a constant in phi6. Rather, clonal selection allows viruses the opportunity to escape the Prisoner's Dilemma. We discuss mechanisms that may afford selfish genotypes an advantage during intrahost competition and the relevance in our system for alternative ecological interactions among viruses.     

High fidelity of murine hepatitis virus replication is decreased in nsp14 exoribonuclease mutants

Replication fidelity of RNA virus genomes is constrained by the opposing necessities of generating sufficient diversity for adaptation and maintaining genetic stability, but it is unclear how the largest viral RNA genomes have evolved and are maintained under these constraints. A coronavirus (CoV) nonstructural protein, nsp14, contains conserved active-site motifs of cellular exonucleases, including DNA proofreading enzymes, and the severe acute respiratory syndrome CoV (SARS-CoV) nsp14 has 3'-to-5' exoribonuclease (ExoN) activity in vitro. Here, we show that nsp14 ExoN remarkably increases replication fidelity of the CoV murine hepatitis virus (MHV). Replacement of conserved MHV ExoN active-site residues with alanines resulted in viable mutant viruses with growth and RNA synthesis defects that during passage accumulated 15-fold more mutations than wild-type virus without changes in growth fitness. The estimated mutation rate for ExoN mutants was similar to that reported for other RNA viruses, whereas that of wild-type MHV was less than the established rates for RNA viruses in general, suggesting that CoVs with intact ExoN replicate with unusually high fidelity. Our results indicate that nsp14 ExoN plays a critical role in prevention or repair of nucleotide incorporation errors during genome replication. The established mutants are unique tools to test the hypothesis that high replication fidelity is required for the evolution and stability of large RNA genomes.     

The smallest ssDNA phage infecting a marine bacterium

In the marine environment, only a few lytic single-stranded DNA (ssDNA) phages have been isolated and characterized, despite the fact that diverse ssDNA bacteriophages have been discovered via metagenomic studies. In this study, we isolated and characterized a new ssDNA phage, vB_RpoMi-Mini, which infects a marine bacterium Ruegeria pomeroyi DSS-3. With a genome size of 4248 bp and only four putative open reading frames (ORF), vB_RpoMi-Mini becomes the smallest ssDNA phage among the known ssDNA phage isolates and represents the DNA bacteriophage with the least number of ORFs. Genome-wide analysis reveals that bacteriophage Mini is distantly related to the known ssDNA phages and belongs to an unclassified ssDNA phage within the Microviridae family. The presence of peptidase in vB_RpoMi-Mini genome further implies that horizontal gene transfer could be an important driving force in the evolution of ssDNA phages. Bacteriophage Mini seems to have lost the spike protein commonly seen in ssDNA phages, suggesting that ssDNA phage can be more diverse than previously thought. Metagenomic analysis indicates that Mini-like phages are widely distributed in the environments. The discovery of vB_RpoMi-Mini expands our understanding of ssDNA phages in nature, and also indicates our dearth of knowledge regarding of ssDNA phages.