Correlation between mutation rate and genome size in riboviruses: mutation rate of bacteriophage Qβ

Fidelity Characterization of Highly Pathogenic Porcine Reproductive and Respiratory Syndrome Virus and NADC30-like Strain

The porcine reproductive and respiratory syndrome virus (PRRSV) has significantly impacted the global pork industry for over three decades. Its high mutation rates and frequent recombination greatly intensifies its epidemic and threat. To explore the fidelity characterization of Chinese highly pathogenic PRRSV JXwn06 and the NADC30-like strain CHsx1401, self-recombination and mutation in PAMs, MARC-145 cells, and pigs were assessed. In vitro, CHsx1401 displayed a higher frequency of recombination junctions and a greater diversity of junction types than JXwn06. In vivo, CHsx1401 exhibited fewer junction types yet maintained a higher junction frequency. Notably, JXwn06 showed more accumulation of mutations. To pinpoint the genomic regions influencing their fidelity, chimeric viruses were constructed, with the exchanged nsp9-10 regions between JXwn06 and CHsx1401. The SJn9n10 strain, which incorporates JXwn06's nsp9-10 into the CHsx1401 genome, demonstrated reduced sensitivity to nucleotide analogs compared to CHsx1401. Conversely, compared with JXwn06, the JSn9n10 strain showed increased sensitivity to these inhibitors. The swapped nsp9-10 also influences the junction frequency and accumulated mutations as their donor strains. The results indicate a propensity for different types of genetic variations between these two strains and further highlight the nsp9-10 region as a critical determinant of their fidelity.

Evolutionary Adaptation of an RNA Bacteriophage to Repeated Freezing and Thawing Cycles

Bacteriophage fitness is determined by factors influencing both their replication within bacteria and their ability to maintain infectivity between infections. The latter becomes particularly crucial under adverse environmental conditions or when host density is low. In such scenarios, the damage experienced by viral particles could lead to the loss of infectivity, which might be mitigated if the virus undergoes evolutionary optimization through replication. In this study, we conducted an evolution experiment involving bacteriophage Qβ, wherein it underwent 30 serial transfers, each involving a cycle of freezing and thawing followed by replication of the surviving viruses. Our findings show that Qβ was capable of enhancing its resistance to this selective pressure through various adaptive pathways that did not impair the virus replicative capacity. Notably, these adaptations predominantly involved mutations located within genes encoding capsid proteins. The adapted populations exhibited higher resistance levels than individual viruses isolated from them, and the latter surpassed those observed in single mutants generated via site-directed mutagenesis. This suggests potential interactions among mutants and mutations. In conclusion, our study highlights the significant role of extracellular selective pressures in driving the evolution of phages, influencing both the genetic composition of their populations and their phenotypic properties.

Mutation rate, selection, and epistasis inferred from RNA virus haplotypes via neural posterior estimation

RNA viruses are particularly notorious for their high levels of genetic diversity, which is generated through the forces of mutation and natural selection. However, disentangling these two forces is a considerable challenge, and this may lead to widely divergent estimates of viral mutation rates, as well as difficulties in inferring the fitness effects of mutations. Here, we develop, test, and apply an approach aimed at inferring the mutation rate and key parameters that govern natural selection, from haplotype sequences covering full-length genomes of an evolving virus population. Our approach employs neural posterior estimation, a computational technique that applies simulation-based inference with neural networks to jointly infer multiple model parameters. We first tested our approach on synthetic data simulated using different mutation rates and selection parameters while accounting for sequencing errors. Reassuringly, the inferred parameter estimates were accurate and unbiased. We then applied our approach to haplotype sequencing data from a serial passaging experiment with the MS2 bacteriophage, a virus that parasites Escherichia coli. We estimated that the mutation rate of this phage is around 0.2 mutations per genome per replication cycle (95% highest density interval: 0.051-0.56). We validated this finding with two different approaches based on single-locus models that gave similar estimates but with much broader posterior distributions. Furthermore, we found evidence for reciprocal sign epistasis between four strongly beneficial mutations that all reside in an RNA stem loop that controls the expression of the viral lysis protein, responsible for lysing host cells and viral egress. We surmise that there is a fine balance between over- and underexpression of lysis that leads to this pattern of epistasis. To recap, we have developed an approach for joint inference of the mutation rate and selection parameters from full haplotype data with sequencing errors and used it to reveal features governing MS2 evolution.

Quasispecies Nature of RNA Viruses: Lessons from the Past

Viral quasispecies are distinct but closely related mutants formed by the disparity in viral genomes due to recombination, mutations, competition, and selection pressure. Theoretical derivation for the origin of a quasispecies is owed to the error-prone replication by polymerase and mutants of RNA replicators. Here, we briefly addressed the theoretical and mathematical origin of quasispecies and their dynamics. The impact of quasispecies for major salient human pathogens is reviewed. In the current global scenario, rapid changes in geographical landscapes favor the origin and selection of mutants. It comes as no surprise that a cauldron of mutants poses a significant risk to public health, capable of causing pandemics. Mutation rates in RNA viruses are magnitudes higher than in DNA organisms, explaining their enhanced virulence and evolvability. RNA viruses cause the most devastating pandemics; for example, members of the Orthomyxoviridae family caused the great influenza pandemic (1918 flu or Spanish flu), the SARS (severe acute respiratory syndrome) and MERS (Middle East respiratory syndrome) outbreak, and the human immunodeficiency viruses (HIV), lentiviruses of the Retroviridae family, caused worldwide devastation. Rapidly evolving RNA virus populations are a daunting challenge for the designing of effective control measures like vaccines. Developing awareness of the evolutionary dispositions of RNA viral mutant spectra and what influences their adaptation and virulence will help curtail outbreaks of past and future pathogens.

Exhaustive mutational analysis of severe acute respiratory syndrome coronavirus 2 ORF3a: An essential component in the pathogen's infectivity cycle

Detailed knowledge of a protein's key residues may assist in understanding its function and designing inhibitors against it. Consequently, such knowledge of one of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)'s proteins is advantageous since the virus is the etiological agent behind one of the biggest health crises of recent times. To that end, we constructed an exhaustive library of bacteria differing from each other by the mutated version of the virus's ORF3a viroporin they harbor. Since the protein is harmful to bacterial growth due to its channel activity, genetic selection followed by deep sequencing could readily identify mutations that abolish the protein's function. Our results have yielded numerous mutations dispersed throughout the sequence that counteract ORF3a's ability to slow bacterial growth. Comparing these data with the conservation pattern of ORF3a within the coronavirinae provided interesting insights: Deleterious mutations obtained in our study corresponded to conserved residues in the protein. However, despite the comprehensive nature of our mutagenesis coverage (108 average mutations per site), we could not reveal all of the protein's conserved residues. Therefore, it is tempting to speculate that our study unearthed positions in the protein pertinent to channel activity, while other conserved residues may correspond to different functionalities of ORF3a. In conclusion, our study provides important information on a key component of SARS-CoV-2 and establishes a procedure to analyze other viroporins comprehensively.

Plant Virus Adaptation to New Hosts: A Multi-scale Approach

Viruses are studied at each level of biological complexity: from within-cells to ecosystems. The same basic evolutionary forces and principles operate at each level: mutation and recombination, selection, genetic drift, migration, and adaptive trade-offs. Great efforts have been put into understanding each level in great detail, hoping to predict the dynamics of viral population, prevent virus emergence, and manage their spread and virulence. Unfortunately, we are still far from this. To achieve these ambitious goals, we advocate for an integrative perspective of virus evolution. Focusing in plant viruses, we illustrate the pervasiveness of the above-mentioned principles. Beginning at the within-cell level, we describe replication modes, infection bottlenecks, and cellular contagion rates. Next, we move up to the colonization of distal tissues, discussing the fundamental role of random events. Then, we jump beyond the individual host and discuss the link between transmission mode and virulence. Finally, at the community level, we discuss properties of virus-plant infection networks. To close this review we propose the multilayer network theory, in which elements at different layers are connected and submit to their own dynamics that feed across layers, resulting in new emerging properties, as a way to integrate information from the different levels.

Standing Genetic Diversity and Transmission Bottleneck Size Drive Adaptation in Bacteriophage Qβ

A critical issue to understanding how populations adapt to new selective pressures is the relative contribution of the initial standing genetic diversity versus that generated de novo. RNA viruses are an excellent model to study this question, as they form highly heterogeneous populations whose genetic diversity can be modulated by factors such as the number of generations, the size of population bottlenecks, or exposure to new environment conditions. In this work, we propagated at nonoptimal temperature (43 °C) two bacteriophage Qβ populations differing in their degree of heterogeneity. Deep sequencing analysis showed that, prior to the temperature change, the most heterogeneous population contained some low-frequency mutations that had previously been detected in the consensus sequences of other Qβ populations adapted to 43 °C. Evolved populations with origin in this ancestor reached similar growth rates, but the adaptive pathways depended on the frequency of these standing mutations and the transmission bottleneck size. In contrast, the growth rate achieved by populations with origin in the less heterogeneous ancestor did depend on the transmission bottleneck size. The conclusion is that viral diversification in a particular environment may lead to the emergence of mutants capable of accelerating adaptation when the environment changes.

Inferring the distribution of fitness effects in patient-sampled and experimental virus populations: two case studies

We here propose an analysis pipeline for inferring the distribution of fitness effects (DFE) from either patient-sampled or experimentally-evolved viral populations, that explicitly accounts for non-Wright-Fisher and non-equilibrium population dynamics inherent to pathogens. We examine the performance of this approach via extensive power and performance analyses, and highlight two illustrative applications - one from an experimentally-passaged RNA virus, and the other from a clinically-sampled DNA virus. Finally, we discuss how such DFE inference may shed light on major research questions in virus evolution, ranging from a quantification of the population genetic processes governing genome size, to the role of Hill-Robertson interference in dictating adaptive outcomes, to the potential design of novel therapeutic approaches to eradicate within-patient viral populations via induced mutational meltdown.

Distinct Retrotransposon Evolution Profile in the Genome of Rabbit (Oryctolagus cuniculus)

Although the rabbit genome has already been annotated, it is mobilome remains largely unknown. Here, multiple pipelines were used to de novo mine and annotate the mobilome in rabbit. Four families and 19 subfamilies of LINE1s, two families and nine subfamilies of SINEs, and 12 ERV families were defined in rabbit based on sequence identity, structural organization, and phylogenetic tree. The analysis of insertion age and polymerase chain reaction suggests that a number of families are very young and may remain active, such as L1B, L1D, OcuSINEA, and OcuERV1. RepeatMasker annotation revealed a distinct transposable element landscape within the genome, with approximately two million copies of SINEs, representing the greatest proportion of the genome (19.61%), followed by LINEs (15.44%), and LTRs (4.11%), respectively, considerably different from most other mammal mobilomes except hedgehog and tree shrew, in which LINEs have the highest proportion. Furthermore, a very high rate of insertion polymorphisms (>85%) for the youngest subfamily (OcuSINEA1) was identified by polymerase chain reaction. The majority of retrotransposon insertions overlapped with protein-coding regions (>80%) and lncRNA (90%) genes. Genomic distribution bias was observed for retrotransposons, with those immediately upstream (-1 kb) and downstream (1 kb) of genes significantly depleted. Local GC content in 50-kb widows had significantly negative correlations with LINE (rs=-0.996) and LTR (rs=-0.829) insertions. The current study revealed a distinct mobilome landscape in rabbit, which will assist in the elucidation of the evolution of the genome of lagomorphs, and even other mammals.

Intra-Population Competition during Adaptation to Increased Temperature in an RNA Bacteriophage

Evolution of RNA bacteriophages of the family Leviviridae is governed by the high error rates of their RNA-dependent RNA polymerases. This fact, together with their large population sizes, leads to the generation of highly heterogeneous populations that adapt rapidly to most changes in the environment. Throughout adaptation, the different mutants that make up a viral population compete with each other in a non-trivial process in which their selective values change over time due to the generation of new mutations. In this work we have characterised the intra-population dynamics of a well-studied levivirus, Qβ, when it is propagated at a higher-than-optimal temperature. Our results show that adapting populations experienced rapid changes that involved the ascent of particular genotypes and the loss of some beneficial mutations of early generation. Artificially reconstructed populations, containing a fraction of the diversity present in actual populations, fixed mutations more rapidly, illustrating how population bottlenecks may guide the adaptive pathways. The conclusion is that, when the availability of beneficial mutations under a particular selective condition is elevated, the final outcome of adaptation depends more on the occasional occurrence of population bottlenecks and how mutations combine in genomes than on the selective value of particular mutations.

The enzymes for genome size increase and maintenance of large (+)RNA viruses

With sizes <50 kb, viral RNA genomes are at the crossroads of genetic, biophysical, and biochemical stability in their host cell. Here, we analyze the enzymatic assets accompanying large RNA genome viruses, mostly based on recent scientific advances in Coronaviridae. We argue that, in addition to the presence of an RNA exonuclease (ExoN), two markers for the large size of viral RNA genomes are (i) the presence of one or more RNA methyltransferases (MTases) and (ii) a specific architecture of the RNA-dependent RNA polymerase active site. We propose that RNA genome expansion and maintenance are driven by an evolutionary ménage-à-trois made of fast and processive RNA polymerases, RNA repair ExoNs, and RNA MTases that relates to the transition between RNA- to DNA-based life.

SARS-CoV-2 biology and variants: anticipation of viral evolution and what needs to be done

The global propagation of SARS-CoV-2 and the detection of a large number of variants, some of which have replaced the original clade to become dominant, underscores the fact that the virus is actively exploring its evolutionary space. The longer high levels of viral multiplication occur - permitted by high levels of transmission -, the more the virus can adapt to the human host and find ways to success. The third wave of the COVID-19 pandemic is starting in different parts of the world, emphasizing that transmission containment measures that are being imposed are not adequate. Part of the consideration in determining containment measures is the rationale that vaccination will soon stop transmission and allow a return to normality. However, vaccines themselves represent a selection pressure for evolution of vaccine-resistant variants, so the coupling of a policy of permitting high levels of transmission/virus multiplication during vaccine roll-out with the expectation that vaccines will deal with the pandemic, is unrealistic. In the absence of effective antivirals, it is not improbable that SARS-CoV-2 infection prophylaxis will involve an annual vaccination campaign against 'dominant' viral variants, similar to influenza prophylaxis. Living with COVID-19 will be an issue of SARS-CoV-2 variants and evolution. It is therefore crucial to understand how SARS-CoV-2 evolves and what constrains its evolution, in order to anticipate the variants that will emerge. Thus far, the focus has been on the receptor-binding spike protein, but the virus is complex, encoding 26 proteins which interact with a large number of host factors, so the possibilities for evolution are manifold and not predictable a priori. However, if we are to mount the best defence against COVID-19, we must mount it against the variants, and to do this, we must have knowledge about the evolutionary possibilities of the virus. In addition to the generic cellular interactions of the virus, there are extensive polymorphisms in humans (e.g. Lewis, HLA, etc.), some distributed within most or all populations, some restricted to specific ethnic populations and these variations pose additional opportunities for/constraints on viral evolution. We now have the wherewithal - viral genome sequencing, protein structure determination/modelling, protein interaction analysis - to functionally characterize viral variants, but access to comprehensive genome data is extremely uneven. Yet, to develop an understanding of the impacts of such evolution on transmission and disease, we must link it to transmission (viral epidemiology) and disease data (patient clinical data), and the population granularities of these. In this editorial, we explore key facets of viral biology and the influence of relevant aspects of human polymorphisms, human behaviour, geography and climate and, based on this, derive a series of recommendations to monitor viral evolution and predict the types of variants that are likely to arise.

cDNA-Derived RNA Phage Assembly Reveals Critical Residues in the Maturation Protein of the Pseudomonas aeruginosa Leviphage PP7

PP7 is a leviphage, with a single-stranded RNA genome, that infects Pseudomonas aeruginosa PAO1. A reverse genetic system for PP7 was previously created by using reverse-transcribed cDNA (PP7_O) from a virion-derived RNA genome. Here, we have found that the PP7_O cDNA contained 20 nucleotide differences from the PP7 genome sequence deposited in the database. We created another reverse genetic system exploiting chemically synthesized cDNA (PP7_S) based on the database sequence. Unlike PP7_O, which yielded infectious PP7 virions, PP7_S-derived particles were incapable of plaque formation on PAO1 cells, which was restored in the PAO1 cells expressing the maturation protein (MP) from PP7_O Using this reverse genetic system, we revealed two amino acid residues involved in the known roles of MP (i.e., adsorption and genome replication), fortuitously providing a lesson that the viral RNA genome sequencing needs functional verification, possibly by a reverse genetic system.IMPORTANCE The biological significance of RNA phages has been largely ignored, ironically, because few studies have focused on RNA phages. As an initial attempt to properly represent RNA phages in the phageome, we previously created, by using reverse-transcribed cDNA, a reverse genetic system for the small RNA phage PP7, which infects the opportunistic human pathogen Pseudomonas aeruginosa We report another system by using chemically synthesized cDNA based on the database genome that has 20 nucleotide differences from the previous cDNA. Investigation of those cDNA-derived phage virions revealed that two amino acids of the maturation protein are crucial for the normal phage lifecycle at different steps. Our study provides insight into the molecular basis for the RNA phage lifecycle and a lesson that the RNA genome sequencing needs to be carefully validated by cDNA-based phage assembly systems.

New Perspectives on Antimicrobial Agents: Remdesivir Treatment for COVID-19

Remdesivir was recently approved by the Food and Drug Administration for the treatment of hospitalized patients with coronavirus disease 2019 (COVID-19). Remdesivir is the prodrug of an adenosine analogue that inhibits viral replication of several RNA virus families, including Coronaviridae. Preclinical data in animal models of coronavirus diseases, including COVID-19, have demonstrated that early treatment with remdesivir leads to improved survival, decreased lung injury, and decreased levels of viral RNA.

Remdesivir was recently approved by the Food and Drug Administration for the treatment of hospitalized patients with coronavirus disease 2019 (COVID-19). Remdesivir is the prodrug of an adenosine analogue that inhibits viral replication of several RNA virus families, including Coronaviridae. Preclinical data in animal models of coronavirus diseases, including COVID-19, have demonstrated that early treatment with remdesivir leads to improved survival, decreased lung injury, and decreased levels of viral RNA. Recent clinical data have demonstrated the clinical activity of remdesivir in terms of faster time to recovery in patients with severe COVID-19 and higher odds of improved clinical status in patients with moderate COVID-19. Here, clinical trials published to date are presented and appraised. Remdesivir’s potential benefits and its favorable adverse-event profile make it an option for the treatment of COVID-19. This article examines the available literature describing remdesivir’s pharmacology, pharmacokinetics, and preclinical and clinical data.

The Role of Temperature in Shaping Mosquito-Borne Viruses Transmission

Mosquito-borne diseases having the greatest impact on human health are typically prevalent in the tropical belt of the world. However, these diseases are conquering temperate regions, raising the question of the role of temperature on their dynamics and expansion. Temperature is one of the most significant abiotic factors affecting, in many ways, insect vectors and the pathogens they transmit. Here, we debate the veracity of this claim by synthesizing current knowledge on the effects of temperature on arboviruses and their vectors, as well as the outcome of their interactions.

Network analysis, sequence and structure dynamics of key proteins of coronavirus and human host, and molecular docking of selected phytochemicals of nine medicinal plants

The novel coronavirus of 2019 (nCoV-19) has become a pandemic, affecting over 205 nations with over 7,410,000 confirmed cases which has resulted to over 418,000 deaths worldwide. This study aimed to identify potential therapeutic compounds and phytochemicals of medicinal plants that have potential to modulate the expression network of genes that are involve in SARS-CoV-2 pathology in human host and to understand the dynamics key proteins involved in the virus-host interactions. The method used include gene network analysis, molecular docking, and sequence and structure dynamics simulations. The results identified DNA-dependent protein kinase (DNA-PK) and Protein kinase CK2 as key players in SARS-CoV-2 lifecycle. Among the predicted drugs compounds, clemizole, monorden, spironolactone and tanespimycin showed high binding energies; among the studied repurposing compounds, remdesivir, simeprevir and valinomycin showed high binding energies; among the predicted acidic compounds, acetylursolic acid and hardwickiic acid gave high binding energies; while among the studied anthraquinones and glycosides compounds, ellagitannin and friedelanone showed high binding energies against 3-Chymotrypsin-like protease (3CL^pro), Papain-like protease (PL^pro), helicase (nsp13), RNA-dependent RNA polymerase (nsp12), 2'-O-ribose methyltransferase (nsp16) of SARS-CoV-2 and DNA-PK and CK2alpha in human. The order of affinity for CoV proteins is 5Y3E > 6NUS > 6JYT > 2XYR > 3VB6. Finally, medicinal plants with phytochemicals such as caffeine, ellagic acid, quercetin and their derivatives could possibly remediate COVID-19.Communicated by Ramaswamy H. Sarma.

Compounds with Therapeutic Potential against Novel Respiratory 2019 Coronavirus

Currently, the expansion of the novel human respiratory coronavirus (known as SARS-CoV-2 [severe acute respiratory syndrome coronavirus 2], COVID-2019 [coronavirus disease 2019], or 2019-nCoV [2019 novel coronavirus]) has stressed the need for therapeutic alternatives to alleviate and stop this new epidemic. The previous epidemics of infections by high-morbidity human coronaviruses, such as SARS-CoV in 2003 and the Middle East respiratory syndrome coronavirus (MERS-CoV) in 2012, prompted the characterization of compounds that could be potentially active against the currently emerging novel coronavirus, SARS-CoV-2. The most promising compound is remdesivir (GS-5734), a nucleotide analog prodrug currently in clinical trials for treating Ebola virus infections. Remdesivir inhibited the replication of SARS-CoV and MERS-CoV in tissue cultures, and it displayed efficacy in nonhuman animal models. In addition, a combination of the human immunodeficiency virus type 1 (HIV-1) protease inhibitors lopinavir/ritonavir and interferon beta (LPV/RTV-IFN-β) was shown to be effective in patients infected with SARS-CoV. LPV/RTV-IFN-β also improved clinical parameters in marmosets and mice infected with MERS-CoV. Remarkably, the therapeutic efficacy of remdesivir appeared to be superior to that of LPV/RTV-IFN-β against MERS-CoV in a transgenic humanized mouse model. The relatively high mortality rates associated with these three novel human coronavirus infections, SARS-CoV, MERS-CoV, and SARS-CoV-2, have suggested that proinflammatory responses might play a role in the pathogenesis. It remains unknown whether the generated inflammatory state should be targeted. Therapeutics that target the coronavirus alone might not be able to reverse highly pathogenic infections. This minireview aims to provide a summary of therapeutic compounds that have shown potential in fighting SARS-CoV-2 infections.

Compounds with Therapeutic Potential against Novel Respiratory 2019 Coronavirus

Currently, the expansion of the novel human respiratory coronavirus (known as SARS-CoV-2 [severe acute respiratory syndrome coronavirus 2], COVID-2019 [coronavirus disease 2019], or 2019-nCoV [2019 novel coronavirus]) has stressed the need for therapeutic alternatives to alleviate and stop this new epidemic. The previous epidemics of infections by high-morbidity human coronaviruses, such as SARS-CoV in 2003 and the Middle East respiratory syndrome coronavirus (MERS-CoV) in 2012, prompted the characterization of compounds that could be potentially active against the currently emerging novel coronavirus, SARS-CoV-2. The most promising compound is remdesivir (GS-5734), a nucleotide analog prodrug currently in clinical trials for treating Ebola virus infections. Remdesivir inhibited the replication of SARS-CoV and MERS-CoV in tissue cultures, and it displayed efficacy in nonhuman animal models. In addition, a combination of the human immunodeficiency virus type 1 (HIV-1) protease inhibitors lopinavir/ritonavir and interferon beta (LPV/RTV-IFN-β) was shown to be effective in patients infected with SARS-CoV. LPV/RTV-IFN-β also improved clinical parameters in marmosets and mice infected with MERS-CoV. Remarkably, the therapeutic efficacy of remdesivir appeared to be superior to that of LPV/RTV-IFN-β against MERS-CoV in a transgenic humanized mouse model. The relatively high mortality rates associated with these three novel human coronavirus infections, SARS-CoV, MERS-CoV, and SARS-CoV-2, have suggested that proinflammatory responses might play a role in the pathogenesis. It remains unknown whether the generated inflammatory state should be targeted. Therapeutics that target the coronavirus alone might not be able to reverse highly pathogenic infections. This minireview aims to provide a summary of therapeutic compounds that have shown potential in fighting SARS-CoV-2 infections.

Compounds with Therapeutic Potential against Novel Respiratory 2019 Coronavirus

Currently, the expansion of the novel human respiratory coronavirus (known as SARS-CoV-2 [severe acute respiratory syndrome coronavirus 2], COVID-2019 [coronavirus disease 2019], or 2019-nCoV [2019 novel coronavirus]) has stressed the need for therapeutic alternatives to alleviate and stop this new epidemic. The previous epidemics of infections by high-morbidity human coronaviruses, such as SARS-CoV in 2003 and the Middle East respiratory syndrome coronavirus (MERS-CoV) in 2012, prompted the characterization of compounds that could be potentially active against the currently emerging novel coronavirus, SARS-CoV-2. The most promising compound is remdesivir (GS-5734), a nucleotide analog prodrug currently in clinical trials for treating Ebola virus infections. Remdesivir inhibited the replication of SARS-CoV and MERS-CoV in tissue cultures, and it displayed efficacy in nonhuman animal models. In addition, a combination of the human immunodeficiency virus type 1 (HIV-1) protease inhibitors lopinavir/ritonavir and interferon beta (LPV/RTV-IFN-β) was shown to be effective in patients infected with SARS-CoV. LPV/RTV-IFN-β also improved clinical parameters in marmosets and mice infected with MERS-CoV. Remarkably, the therapeutic efficacy of remdesivir appeared to be superior to that of LPV/RTV-IFN-β against MERS-CoV in a transgenic humanized mouse model. The relatively high mortality rates associated with these three novel human coronavirus infections, SARS-CoV, MERS-CoV, and SARS-CoV-2, have suggested that proinflammatory responses might play a role in the pathogenesis. It remains unknown whether the generated inflammatory state should be targeted. Therapeutics that target the coronavirus alone might not be able to reverse highly pathogenic infections. This minireview aims to provide a summary of therapeutic compounds that have shown potential in fighting SARS-CoV-2 infections.

Molecular basis of genetic variation of viruses

Genetic variation is a necessity of all biological systems. Viruses use all known mechanisms of variation; mutation, several forms of recombination, and segment reassortment in the case of viruses with a segmented genome. These processes are intimately connected with the replicative machineries of viruses, as well as with fundamental physical-chemical properties of nucleotides when acting as template or substrate residues. Recombination has been viewed as a means to rescue viable genomes from unfit parents or to produce large modifications for the exploration of phenotypic novelty. All types of genetic variation can act conjointly as blind processes to provide the raw materials for adaptation to the changing environments in which viruses must replicate. A distinction is made between mechanistically unavoidable and evolutionarily relevant mutation and recombination.

Viral quasispecies

Viral quasispecies refers to a population structure that consists of extremely large numbers of variant genomes, termed mutant spectra, mutant swarms or mutant clouds. Fueled by high mutation rates, mutants arise continually, and they change in relative frequency as viral replication proceeds. The term quasispecies was adopted from a theory of the origin of life in which primitive replicons) consisted of mutant distributions, as found experimentally with present day RNA viruses. The theory provided a new definition of wild type, and a conceptual framework for the interpretation of the adaptive potential of RNA viruses that contrasted with classical studies based on consensus sequences. Standard clonal analyses and deep sequencing methodologies have confirmed the presence of myriads of mutant genomes in viral populations, and their participation in adaptive processes. The quasispecies concept applies to any biological entity, but its impact is more evident when the genome size is limited and the mutation rate is high. This is the case of the RNA viruses, ubiquitous in our biosphere, and that comprise many important pathogens. In virology, quasispecies are defined as complex distributions of closely related variant genomes subjected to genetic variation, competition and selection, and that may act as a unit of selection. Despite being an integral part of their replication, high mutation rates have an upper limit compatible with inheritable information. Crossing such a limit leads to RNA virus extinction, a transition that is the basis of an antiviral design termed lethal mutagenesis.

Harnessing host-virus evolution in antiviral therapy and immunotherapy

Pathogen resistance and development costs are major challenges in current approaches to antiviral therapy. The high error rate of RNA synthesis and reverse‐transcription confers genome plasticity, enabling the remarkable adaptability of RNA viruses to antiviral intervention. However, this property is coupled to fundamental constraints including limits on the size of information available to manipulate complex hosts into supporting viral replication. Accordingly, RNA viruses employ various means to extract maximum utility from their informationally limited genomes that, correspondingly, may be leveraged for effective host‐oriented therapies. Host‐oriented approaches are becoming increasingly feasible because of increased availability of bioactive compounds and recent advances in immunotherapy and precision medicine, particularly genome editing, targeted delivery methods and RNAi. In turn, one driving force behind these innovations is the increasingly detailed understanding of evolutionarily diverse host–virus interactions, which is the key concern of an emerging field, neo‐virology. This review examines biotechnological solutions to disease and other sustainability issues of our time that leverage the properties of RNA and DNA viruses as developed through co‐evolution with their hosts.

Evolutionary Dynamics in the RNA Bacteriophage Qβ Depends on the Pattern of Change in Selective Pressures

The rate of change in selective pressures is one of the main factors that determines the likelihood that populations can adapt to stress conditions. Generally, the reduction in the population size that accompanies abrupt environmental changes makes it difficult to generate and select adaptive mutations. However, in systems with high genetic diversity, as happens in RNA viruses, mutations with beneficial effects under new conditions can already be present in the population, facilitating adaptation. In this work, we have propagated an RNA bacteriophage (Qβ) at temperatures higher than the optimum, following different patterns of change. We have determined the fitness values and the consensus sequences of all lineages throughout the evolutionary process in order to establish correspondences between fitness variations and adaptive pathways. Our results show that populations subjected to a sudden temperature change gain fitness and fix mutations faster than those subjected to gradual changes, differing also in the particular selected mutations. The life-history of populations prior to the environmental change has great importance in the dynamics of adaptation. The conclusion is that in the bacteriophage Qβ, the standing genetic diversity together with the rate of temperature change determine both the rapidity of adaptation and the followed evolutionary pathways.

STABILITY OF A STOCHASTICALLY PERTURBED MODEL OF INTRACELLULAR SINGLE-STRANDED RNA VIRUS REPLICATION

Compared to the replication of double-stranded RNA and DNA viruses, the replication of single-stranded viruses requires the production of a number of intermediate strands that serve as templates for the synthesis of genomic-sense strands. Two theoretical extreme mechanisms for replication for such single-stranded viruses have been proposed; one extreme being represented by the so-called linear stamping machine and the opposite extreme by the exponential growth. Of course, real systems are more complex and examples have been described in which a combination of such extreme mechanisms can also occur: a fraction of the produced progeny resulting from a stamping-machine type of replication that uses the parental genome as template, whereas other fraction of the progeny results from the replication of other progeny genomes. Martínez et al. 1 , Sardanyés et al. 2 and Fornés et al. 3 suggested and analyzed a deterministic model of single-stranded RNA (ssRNA) virus intracellular replication that incorporated variability in the replication mechanisms.

To explore how stochasticity can affect this mixed-model principal properties, in this paper, we consider the stability of a stochastically perturbed model of ssRNA virus replication within a cell. Using the direct Lyapunov method, we found sufficient conditions for the stability in probability of equilibrium states for this model. This result confirms that this heterogeneous model of single-stranded RNA virus replication is stable with respect to stochastic perturbations of the environment.

Stochastic Modeling and Simulation of Viral Evolution

RNA viruses comprise vast populations of closely related, but highly genetically diverse, entities known as quasispecies. Understanding the mechanisms by which this extreme diversity is generated and maintained is fundamental when approaching viral persistence and pathobiology in infected hosts. In this paper, we access quasispecies theory through a mathematical model based on the theory of multitype branching processes, to better understand the roles of mechanisms resulting in viral diversity, persistence and extinction. We accomplish this understanding by a combination of computational simulations and the theoretical analysis of the model. In order to perform the simulations, we have implemented the mathematical model into a computational platform capable of running simulations and presenting the results in a graphical format in real time. Among other things, we show that the establishment of virus populations may display four distinct regimes from its introduction into new hosts until achieving equilibrium or undergoing extinction. Also, we were able to simulate different fitness distributions representing distinct environments within a host which could either be favorable or hostile to the viral success. We addressed the most used mechanisms for explaining the extinction of RNA virus populations called lethal mutagenesis and mutational meltdown. We were able to demonstrate a correspondence between these two mechanisms implying the existence of a unifying principle leading to the extinction of RNA viruses.

Evolutionary adaptation of an RNA bacteriophage to the simultaneous increase in the within-host and extracellular temperatures

Bacteriophages are the most numerous biological entities on Earth. They are on the basis of most ecosystems, regulating the diversity and abundance of bacterial populations and contributing to the nutrient and energy cycles. Bacteriophages have two well differentiated phases in their life cycle, one extracellular, in which they behave as inert particles, and other one inside their hosts, where they replicate to give rise to a progeny. In both phases they are exposed to environmental conditions that often act as selective pressures that limit both their survival in the environment and their ability to replicate, two fitness traits that frequently cannot be optimised simultaneously. In this study we have analysed the evolutionary ability of an RNA bacteriophage, the bacteriophage Qβ, when it is confronted with a temperature increase that affects both the extracellular and the intracellular media. Our results show that Qβ can optimise its survivability when exposed to short-term high temperature extracellular heat shocks, as well as its replicative ability at higher-than-optimal temperature. Mutations responsible for simultaneous adaptation were the same as those selected when adaptation to each condition proceeded separately, showing the absence of important trade-offs between survival and reproduction in this virus.

Structural and Functional Basis of the Fidelity of Nucleotide Selection by Flavivirus RNA-Dependent RNA Polymerases

Viral RNA-dependent RNA polymerases (RdRps) play a central role not only in viral replication, but also in the genetic evolution of viral RNAs. After binding to an RNA template and selecting 5'-triphosphate ribonucleosides, viral RdRps synthesize an RNA copy according to Watson-Crick base-pairing rules. The copy process sometimes deviates from both the base-pairing rules specified by the template and the natural ribose selectivity and, thus, the process is error-prone due to the intrinsic (in)fidelity of viral RdRps. These enzymes share a number of conserved amino-acid sequence strings, called motifs A-G, which can be defined from a structural and functional point-of-view. A co-relation is gradually emerging between mutations in these motifs and viral genome evolution or observed mutation rates. Here, we review our current knowledge on these motifs and their role on the structural and mechanistic basis of the fidelity of nucleotide selection and RNA synthesis by Flavivirus RdRps.

Proofreading-Deficient Coronaviruses Adapt for Increased Fitness over Long-Term Passage without Reversion of Exoribonuclease-Inactivating Mutations

The coronavirus (CoV) RNA genome is the largest among the single-stranded positive-sense RNA viruses. CoVs encode a proofreading 3'-to-5' exoribonuclease within nonstructural protein 14 (nsp14-ExoN) that is responsible for CoV high-fidelity replication. Alanine substitution of ExoN catalytic residues [ExoN(-)] in severe acute respiratory syndrome-associated coronavirus (SARS-CoV) and murine hepatitis virus (MHV) disrupts ExoN activity, yielding viable mutant viruses with defective replication, up to 20-fold-decreased fidelity, and increased susceptibility to nucleoside analogues. To test the stability of the ExoN(-) genotype and phenotype, we passaged MHV-ExoN(-) 250 times in cultured cells (P250), in parallel with wild-type MHV (WT-MHV). Compared to MHV-ExoN(-) P3, MHV-ExoN(-) P250 demonstrated enhanced replication and increased competitive fitness without reversion at the ExoN(-) active site. Furthermore, MHV-ExoN(-) P250 was less susceptible than MHV-ExoN(-) P3 to multiple nucleoside analogues, suggesting that MHV-ExoN(-) was under selection for increased replication fidelity. We subsequently identified novel amino acid changes within the RNA-dependent RNA polymerase and nsp14 of MHV-ExoN(-) P250 that partially accounted for the reduced susceptibility to nucleoside analogues. Our results suggest that increased replication fidelity is selected in ExoN(-) CoVs and that there may be a significant barrier to ExoN(-) reversion. These results also support the hypothesis that high-fidelity replication is linked to CoV fitness and indicate that multiple replicase proteins could compensate for ExoN functions during replication.IMPORTANCE Uniquely among RNA viruses, CoVs encode a proofreading exoribonuclease (ExoN) in nsp14 that mediates high-fidelity RNA genome replication. Proofreading-deficient CoVs with disrupted ExoN activity [ExoN(-)] either are nonviable or have significant defects in replication, RNA synthesis, fidelity, fitness, and virulence. In this study, we showed that ExoN(-) murine hepatitis virus can adapt during long-term passage for increased replication and fitness without reverting the ExoN-inactivating mutations. Passage-adapted ExoN(-) mutants also demonstrate increasing resistance to nucleoside analogues that is explained only partially by secondary mutations in nsp12 and nsp14. These data suggest that enhanced resistance to nucleoside analogues is mediated by the interplay of multiple replicase proteins and support the proposed link between CoV fidelity and fitness.

Different rates of spontaneous mutation of chloroplastic and nuclear viroids as determined by high-fidelity ultra-deep sequencing

Mutation rates vary by orders of magnitude across biological systems, being higher for simpler genomes. The simplest known genomes correspond to viroids, subviral plant replicons constituted by circular non-coding RNAs of few hundred bases. Previous work has revealed an extremely high mutation rate for chrysanthemum chlorotic mottle viroid, a chloroplast-replicating viroid. However, whether this is a general feature of viroids remains unclear. Here, we have used high-fidelity ultra-deep sequencing to determine the mutation rate in a common host (eggplant) of two viroids, each representative of one family: the chloroplastic eggplant latent viroid (ELVd, Avsunviroidae) and the nuclear potato spindle tuber viroid (PSTVd, Pospiviroidae). This revealed higher mutation frequencies in ELVd than in PSTVd, as well as marked differences in the types of mutations produced. Rates of spontaneous mutation, quantified in vivo using the lethal mutation method, ranged from 1/1000 to 1/800 for ELVd and from 1/7000 to 1/3800 for PSTVd depending on sequencing run. These results suggest that extremely high mutability is a common feature of chloroplastic viroids, whereas the mutation rates of PSTVd and potentially other nuclear viroids appear significantly lower and closer to those of some RNA viruses.

Spontaneous mutations are the ultimate source of genetic variation and their characterization provides fundamental information about evolutionary processes. The highest mutation rate so far described corresponds to a hammerhead viroid infecting plant chloroplasts. Viroids are plant-exclusive parasites constituted by 250–400 nt-long, non-protein-coding RNAs, and are divided into two families with distinct mechanisms of replication and localization: chloroplastic (Avsunviroidae), and nuclear (Pospiviroidae). Here, we have used high-fidelity ultra-deep sequencing to compare side by side the mutation rates of one representative member of each viroid family in the same host. We found that the mutation rate of the nuclear viroid was several fold lower than that of the chloroplastic viroid.

Veterinary pharmaceuticals in aqueous systems and associated effects: an update

Environmental studies have shown that pharmaceuticals can contaminate aqueous matrices, such as groundwater, surface water, sediment as well as aquatic flora and fauna. Effluents from sewage and wastewater treatment plants, pharmaceutical industries and hospitals have been implicated in such contamination. Recent studies have however revealed significant concentrations of pharmaceuticals in wastewater from animal facilities in proximal aquatic habitats. Furthermore, epidemiological studies have shown a consistent positive correlation between exposure to some drugs of veterinary importance and increased adverse effects in aquatic biota largely due to induction of endocrine disruption, antibiotic resistance, neurotoxicity, genotoxicity and oxidative stress. The aquatic habitats and associated biota are important in the maintenance of global ecosystem and food chain. For this reason, anything that compromises the integrity and functions of the aquatic environment may lead to major upset in the world's ecosystems. Therefore, knowledge about this route of exposure cannot be neglected and monitoring of their occurrence in the environment is required. This review focuses on scientific evidence that link the presence of pharmaceuticals in aqueous matrices to animal production facilities and presents means to reduce the occurrence of veterinary pharmaceutical residues in the aquatic habitats.

Mechanisms of viral mutation

The remarkable capacity of some viruses to adapt to new hosts and environments is highly dependent on their ability to generate de novo diversity in a short period of time. Rates of spontaneous mutation vary amply among viruses. RNA viruses mutate faster than DNA viruses, single-stranded viruses mutate faster than double-strand virus, and genome size appears to correlate negatively with mutation rate. Viral mutation rates are modulated at different levels, including polymerase fidelity, sequence context, template secondary structure, cellular microenvironment, replication mechanisms, proofreading, and access to post-replicative repair. Additionally, massive numbers of mutations can be introduced by some virus-encoded diversity-generating elements, as well as by host-encoded cytidine/adenine deaminases. Our current knowledge of viral mutation rates indicates that viral genetic diversity is determined by multiple virus- and host-dependent processes, and that viral mutation rates can evolve in response to specific selective pressures.

Effect of mismatch repair on the mutation rate of bacteriophage ϕX174

Viral mutation rates vary widely in nature, yet the mechanistic and evolutionary determinants of this variability remain unclear. Small DNA viruses mutate orders of magnitude faster than their hosts despite using host-encoded polymerases for replication, which suggests these viruses may avoid post-replicative repair. Supporting this, the genome of bacteriophage ϕX174 is completely devoid of GATC sequence motifs, which are required for methyl-directed mismatch repair in Escherichia coli. Here, we show that restoration of the randomly expected number of GATC sites leads to an eightfold reduction in the rate of spontaneous mutation of the phage, without severely impairing its replicative capacity over the short term. However, the efficacy of mismatch repair in the presence of GATC sites is limited by inefficient methylation of the viral DNA. Therefore, both GATC avoidance and DNA under-methylation elevate the mutation rate of the phage relative to that of the host. We also found that the effects of GATC sites on the phage mutation rate vary extensively depending on their specific location within the phage genome. Finally, the mutation rate reduction afforded by GATC sites is fully reverted under stress conditions, which up-regulate repair pathways and expression of error-prone host polymerases such as heat and treatment with the base analog 5-fluorouracil, suggesting that access to repair renders the phage sensitive to stress-induced mutagenesis.

Oligomeric viral proteins: small in size, large in presence

Viruses are obligate parasites that rely heavily on host cellular processes for replication. The small number of proteins typically encoded by a virus is faced with selection pressures that lead to the evolution of distinctive structural properties, allowing each protein to maintain its function under constraints such as small genome size, high mutation rate, and rapidly changing fitness conditions. One common strategy for this evolution is to utilize small building blocks to generate protein oligomers that assemble in multiple ways, thereby diversifying protein function and regulation. In this review, we discuss specific cases that illustrate how oligomerization is used to generate a single defined functional state, to modulate activity via different oligomeric states, or to generate multiple functional forms via different oligomeric states.

Impact of increased mutagenesis on adaptation to high temperature in bacteriophage Qβ

RNA viruses replicate with very high error rates, which makes them more sensitive to additional increases in this parameter. This fact has inspired an antiviral strategy named lethal mutagenesis, which is based on the artificial increase of the error rate above a threshold incompatible with virus infectivity. A relevant issue concerning lethal mutagenesis is whether incomplete treatments might enhance the adaptive possibilities of viruses. We have addressed this question by subjecting an RNA virus, the bacteriophage Qβ, to different transmission regimes in the presence or the absence of sublethal concentrations of the mutagenic nucleoside analogue 5-azacytidine (AZC). Populations obtained were subsequently exposed to a non-optimal temperature and analyzed to determine their consensus sequences. Our results show that previously mutagenized populations rapidly fixed a specific set of mutations upon propagation at the new temperature, suggesting that the expansion of the mutant spectrum caused by AZC has an influence on later evolutionary behavior.

Single-Cell Analysis of RNA Virus Infection Identifies Multiple Genetically Diverse Viral Genomes within Single Infectious Units

Genetic diversity enables a virus to colonize novel hosts, evade immunity, and evolve drug resistance. However, viral diversity is typically assessed at the population level. Given the existence of cell-to-cell variation, it is critical to understand viral genetic structure at the single-cell level. By combining single-cell isolation with ultra-deep sequencing, we characterized the genetic structure and diversity of a RNA virus shortly after single-cell bottlenecks. Full-length sequences from 881 viral plaques derived from 90 individual cells reveal that sequence variants pre-existing in different viral genomes can be co-transmitted within the same infectious unit to individual cells. Further, the rate of spontaneous virus mutation varies across individual cells, and early production of diversity depends on the viral yield of the very first infected cell. These results unravel genetic and structural features of a virus at the single-cell level, with implications for viral diversity and evolution.

•

Genetic diversity is found within viral units infecting an individual cell

•

Viral particles within the same infectious unit interact genetically and functionally

•

The production of viral genetic diversity is marked by widespread cell-to-cell variation

•

Single-cell virus yields determine the early production viral genetic diversity

Genetic Instability of RNA Viruses

Despite having very limited coding capacity, RNA viruses are able to withstand challenge of antiviral drugs, cause epidemics in previously exposed human populations, and, in some cases, infect multiple host species. They are able to achieve this by virtue of their ability to multiply very rapidly, coupled with their extraordinary degree of genetic heterogeneity. RNA viruses exist not as single genotypes, but as a swarm of related variants, and this genomic diversity is an essential feature of their biology. RNA viruses have a variety of mechanisms that act in combination to determine their genetic heterogeneity. These include polymerase fidelity, error-mitigation mechanisms, genomic recombination, and different modes of genome replication. RNA viruses can vary in their ability to tolerate mutations, or “genetic robustness,” and several factors contribute to this. Finally, there is evidence that some RNA viruses exist close to a threshold where polymerase error rate has evolved to maximize the possible sequence space available, while avoiding the accumulation of a lethal load of deleterious mutations. We speculate that different viruses have evolved different error rates to complement the different “life-styles” they possess.

How Do Species Interactions Affect Evolutionary Dynamics Across Whole Communities?

Theories of how species evolve in changing environments mostly consider single species in isolation or pairs of interacting species. Yet all organisms live in diverse communities containing many hundreds of species. This review discusses how species interactions influence the evolution of constituent species across whole communities. When species interactions are weak or inconsistent, evolutionary dynamics should be predictable by factors identified by single-species theory. Stronger species interactions, however, can alter evolutionary outcomes and either dampen or promote evolution of constituent species depending on the number of species and the distribution of interaction strengths across the interaction network. Genetic interactions, such as horizontal gene transfer, might also affect evolutionary outcomes. These evolutionary mechanisms in turn affect whole-community properties, such as the level of ecosystem functioning. Successful management of both ecosystems and focal species requires new understanding of evolutionary interactions across whole communities.

Extremely High Mutation Rate of HIV-1 In Vivo

Rates of spontaneous mutation critically determine the genetic diversity and evolution of RNA viruses. Although these rates have been characterized in vitro and in cell culture models, they have seldom been determined in vivo for human viruses. Here, we use the intrapatient frequency of premature stop codons to quantify the HIV-1 genome-wide rate of spontaneous mutation in DNA sequences from peripheral blood mononuclear cells. This reveals an extremely high mutation rate of (4.1 ± 1.7) × 10⁻³ per base per cell, the highest reported for any biological entity. Sequencing of plasma-derived sequences yielded a mutation frequency 44 times lower, indicating that a large fraction of viral genomes are lethally mutated and fail to reach plasma. We show that the HIV-1 reverse transcriptase contributes only 2% of mutations, whereas 98% result from editing by host cytidine deaminases of the A3 family. Hypermutated viral sequences are less abundant in patients showing rapid disease progression compared to normal progressors, highlighting the antiviral role of A3 proteins. However, the amount of A3-mediated editing varies broadly, and we find that low-edited sequences are more abundant among rapid progressors, suggesting that suboptimal A3 activity might enhance HIV-1 genetic diversity and pathogenesis.

The rate of spontaneous mutation of the HIV-1 genome within its human host is exceptionally high, is mostly driven by host cytidine deaminases, and probably plays a role in disease progression.

The high levels of genetic diversity of the HIV-1 virus grant it the ability to escape the immune system, to rapidly evolve drug resistance, and to circumvent vaccination strategies. However, our knowledge of HIV-1 mutation rates has been largely restricted to in vitro and cell culture studies because of the inherent complexity of measuring these rates in vivo. Here, by analyzing the frequency of premature stop codons, we show that the HIV-1 mutation rate in vivo is two orders of magnitude higher than that predicted by in vitro studies, making it the highest reported mutation rate for any biological system. A large component of this rate is from host cellular cytidine deaminases, which induce mutations in the viral DNA as a defense mechanism. While the HIV-1 genome is hypermutated in blood cells, only a very small fraction of these mutations reach the plasma, indicating that many viruses are defective as a result of the extremely high mutation load. In addition, we find that the HIV-1 mutation rate tends to be higher in patients showing normal disease progression than in those undergoing rapid progression, emphasizing the negative impact on viral fitness of hypermutation by host cytidine deaminases. However, we also observe subpopulations of weakly-mutated viral genomes whose sequence diversity may influence viral pathogenesis. Our work highlights the fine balance for HIV-1 between enough mutation to evade host responses and too much mutation that can inactivate the virus.

Changes in protein domains outside the catalytic site of the bacteriophage Qβ replicase reduce the mutagenic effect of 5-azacytidine

The high genetic heterogeneity and great adaptability of RNA viruses are ultimately caused by the low replication fidelity of their polymerases. However, single amino acid substitutions that modify replication fidelity can evolve in response to mutagenic treatments with nucleoside analogues. Here, we investigated how two independent mutants of the bacteriophage Qβ replicase (Thr210Ala and Tyr410His) reduce sensitivity to the nucleoside analogue 5-azacytidine (AZC). Despite being located outside the catalytic site, both mutants reduced the mutation frequency in the presence of the drug. However, they did not modify the type of AZC-induced substitutions, which was mediated mainly by ambiguous base pairing of the analogue with purines. Furthermore, the Thr210Ala and Tyr410His substitutions had little or no effect on replication fidelity in untreated viruses. Also, both substitutions were costly in the absence of AZC or when the action of the drug was suppressed by adding an excess of natural pyrimidines (uridine or cytosine). Overall, the phenotypic properties of these two mutants were highly convergent, despite the mutations being located in different domains of the Qβ replicase. This suggests that treatment with a given nucleoside analogue tends to select for a unique functional response in the viral replicase.

IMPORTANCE

In the last years, artificial increase of the replication error rate has been proposed as an antiviral therapy. In this study, we investigated the mechanisms by which two substitutions in the Qβ replicase confer partial resistance to the mutagenic nucleoside analogue AZC. As opposed to previous work with animal viruses, where different mutations selected sequentially conferred nucleoside analogue resistance through different mechanisms, our results suggest that there are few or no alternative AZC resistance phenotypes in Qβ. Also, despite resistance mutations being highly costly in the absence of the drug, there was no sequential fixation of secondary mutations. Bacteriophage Qβ is the virus with the highest reported mutation rate, which should make it particularly sensitive to nucleoside analogue treatments, probably favoring resistance mutations even if they incur high costs. The results are also relevant for understanding the possible pathways by which fidelity of the replication machinery can be modified.

Contribution of silent mutations to thermal adaptation of RNA bacteriophage Qβ

Changes in protein function and other biological properties, such as RNA structure, are crucial for adaptation of organisms to novel or inhibitory environments. To investigate how mutations that do not alter amino acid sequence may be positively selected, we performed a thermal adaptation experiment using the single-stranded RNA bacteriophage Qβ in which the culture temperature was increased from 37.2°C to 41.2°C and finally to an inhibitory temperature of 43.6°C in a stepwise manner in three independent lines. Whole-genome analysis revealed 31 mutations, including 14 mutations that did not result in amino acid sequence alterations, in this thermal adaptation. Eight of the 31 mutations were observed in all three lines. Reconstruction and fitness analyses of Qβ strains containing only mutations observed in all three lines indicated that five mutations that did not result in amino acid sequence changes but increased the amplification ratio appeared in the course of adaptation to growth at 41.2°C. Moreover, these mutations provided a suitable genetic background for subsequent mutations, altering the fitness contribution from deleterious to beneficial. These results clearly showed that mutations that do not alter the amino acid sequence play important roles in adaptation of this single-stranded RNA virus to elevated temperature.

IMPORTANCE

Recent studies using whole-genome analysis technology suggested the importance of mutations that do not alter the amino acid sequence for adaptation of organisms to novel environmental conditions. It is necessary to investigate how these mutations may be positively selected and to determine to what degree such mutations that do not alter amino acid sequences contribute to adaptive evolution. Here, we report the roles of these silent mutations in thermal adaptation of RNA bacteriophage Qβ based on experimental evolution during which Qβ showed adaptation to growth at an inhibitory temperature. Intriguingly, four synonymous mutations and one mutation in the untranslated region that spread widely in the Qβ population during the adaptation process at moderately high temperature provided a suitable genetic background to alter the fitness contribution of subsequent mutations from deleterious to beneficial at a higher temperature.

RNA virus evolution at variable error rate

ABSTRACT: 

RNA viruses replicate their genomes with very high error rates, which leads to the generation of a large genetic diversity that makes them highly adaptable to most environmental pressures, including antiviral drugs and immune responses. However, since most mutations are deleterious, an excess of errors can be very negative for RNA viruses, entailing that error rates must be finely regulated. Currently, the manipulation of the error rate is emerging as a promising antiviral therapy that could minimize the problem of virus adaptation to classical treatments. This review provides a detailed analysis of the different outcomes that can result from the variation of the error rate in RNA viruses, on the basis of the more relevant findings obtained in experimental studies.

Thinking Outside the Triangle: Replication Fidelity of the Largest RNA Viruses

When judged by ubiquity, adaptation, and emergence of new diseases, RNA viruses are arguably the most successful biological organisms. This success has been attributed to a defect of sorts: high mutation rates (low fidelity) resulting in mutant swarms that allow rapid selection for fitness in new environments. Studies of viruses with small RNA genomes have identified fidelity determinants in viral RNA-dependent RNA polymerases and have shown that RNA viruses likely replicate within a limited fidelity range to maintain fitness. In this review we compare the fidelity of small RNA viruses with that of the largest RNA viruses, the coronaviruses. Coronaviruses encode the first known viral RNA proofreading exoribonuclease, a function that likely allowed expansion of the coronavirus genome and that dramatically increases replication fidelity and the range of tolerated variation. We propose models for regulation of coronavirus fidelity and discuss the implications of altered fidelity for RNA virus replication, pathogenesis, and evolution.

Variability in the mutation rates of RNA viruses

ABSTRACT: 

It is well established that RNA viruses show extremely high mutation rates, but less attention has been paid to the fact that their mutation rates also vary strongly, from 10-6 to 10-4 substitutions per nucleotide per cell infection. The causes explaining this variability are still poorly understood, but candidate factors are the viral genome size and polarity, host-specific gene expression patterns, or the intracellular environment. Differences between animal and plant viruses, or between arthropod-borne and directly transmitted viruses have also been postulated. Finally, RNA viruses may be able to regulate the rate at which new mutations spread in the population by modifying features of the viral infection cycle, such as lysis time.

Delayed lysis confers resistance to the nucleoside analogue 5-fluorouracil and alleviates mutation accumulation in the single-stranded DNA bacteriophage ϕX174

Rates of spontaneous mutation determine viral fitness and adaptability. In RNA viruses, treatment with mutagenic nucleoside analogues selects for polymerase variants with increased fidelity, showing that viral mutation rates can be adjusted in response to imposed selective pressures. However, this type of resistance is not possible in viruses that do not encode their own polymerases, such as single-stranded DNA viruses. We previously showed that serial passaging of bacteriophage ϕX174 in the presence of the nucleoside analogue 5-fluorouracil (5-FU) favored substitutions in the lysis protein E (P. Domingo-Calap, M. Pereira-Gomez, and R. Sanjuán, J. Virol. 86:: 9640-9646, 2012, doi:10.1128/JVI.00613-12). Here, we found that approximately half (6/12) of the amino acid replacements in the N-terminal region of this protein led to delayed lysis, and two of these changes (V2A and D8A) also conferred partial resistance to 5-FU. By delaying lysis, the V2A and D8A substitutions allowed the virus to increase the burst size per cell in the presence of 5-FU. Furthermore, these substitutions tended to alleviate drug-induced mutagenesis by reducing the number of rounds of copying required for population growth, revealing a new mechanism of resistance. This form of mutation rate regulation may also be utilized by other viruses whose replication mode is similar to that of bacteriophage ϕX174.

IMPORTANCE

Many viruses display high rates of spontaneous mutations due to defects in proofreading or postreplicative repair, allowing them to rapidly adapt to changing environments. Viral mutation rates may have been optimized to achieve high adaptability without incurring an excessive genetic load. Supporting this, RNA viruses subjected to chemical mutagenesis treatments have been shown to evolve higher-fidelity polymerases. However, many viruses cannot modulate replication fidelity because they do not encode their own polymerase. Here, we show a new mechanism for regulating viral mutation rates. We found that, under mutagenic conditions, the single-stranded bacteriophage ϕX174 evolved delayed lysis, and that this allowed the virus to increase the amount of progeny produced per cell. As a result, the viral population was amplified in fewer infection cycles, reducing the chances for mutation appearance.

Cell tropism predicts long-term nucleotide substitution rates of mammalian RNA viruses

The high rates of RNA virus evolution are generally attributed to replication with error-prone RNA-dependent RNA polymerases. However, these long-term nucleotide substitution rates span three orders of magnitude and do not correlate well with mutation rates or selection pressures. This substitution rate variation may be explained by differences in virus ecology or intrinsic genomic properties. We generated nucleotide substitution rate estimates for mammalian RNA viruses and compiled comparable published rates, yielding a dataset of 118 substitution rates of structural genes from 51 different species, as well as 40 rates of non-structural genes from 28 species. Through ANCOVA analyses, we evaluated the relationships between these rates and four ecological factors: target cell, transmission route, host range, infection duration; and three genomic properties: genome length, genome sense, genome segmentation. Of these seven factors, we found target cells to be the only significant predictors of viral substitution rates, with tropisms for epithelial cells or neurons (P<0.0001) as the most significant predictors. Further, one-tailed t-tests showed that viruses primarily infecting epithelial cells evolve significantly faster than neurotropic viruses (P<0.0001 and P<0.001 for the structural genes and non-structural genes, respectively). These results provide strong evidence that the fastest evolving mammalian RNA viruses infect cells with the highest turnover rates: the highly proliferative epithelial cells. Estimated viral generation times suggest that epithelial-infecting viruses replicate more quickly than viruses with different cell tropisms. Our results indicate that cell tropism is a key factor in viral evolvability.

RNA viruses are the fastest evolving human pathogens, making their treatment and control difficult. Compared to DNA viruses, RNA viruses replicate with much lower fidelity, which can explain why RNA viruses evolve significantly faster than most DNA viruses. However, there is tremendous variation among the evolutionary rates of different RNA viruses, which is not explained by variation in mutation rates. Here we present a survey of mammalian RNA virus rates of evolution, and a comprehensive comparison of these rates to different properties of virus genomic architecture and ecology. We found that cell tropism is the most significant predictor of long-term rates of mammalian RNA virus evolution. For instance, viruses targeting epithelial cells evolve significantly faster than viruses that target neurons. Our results provide mechanistic insight into why viruses that infect respiratory and gastrointestinal epithelia have been difficult to control.