Josep Sardanyés

Simple genomes, complex interactions: Epistasis in RNA virus

Owed to their reduced size and low number of proteins encoded, RNA viruses and other subviral pathogens are often considered as being genetically too simple. However, this structural simplicity also creates the necessity for viral RNA sequences to encode for more than one protein and for proteins to carry out multiple functions, all together resulting in complex patterns of genetic interactions. In this work we will first review the experimental studies revealing that the architecture of viral genomes is dominated by antagonistic interactions among loci. Second, we will also review mathematical models and provide a description of computational tools for the study of RNA virus dynamics and evolution. As an application of these tools, we will finish this review article by analyzing a stochastic bit-string model of in silico virus replication. This model analyzes the interplay between epistasis and the mode of replication on determining the population load of deleterious mutations. The model suggests that, for a given mutation rate, the deleterious mutational load is always larger when epistasis is predominantly antagonistic than when synergism is the rule. However, the magnitude of this effect is larger if replication occurs geometrically than if it proceeds linearly.

Tempo and Mode of Plant RNA Virus Escape from RNA Interference-Mediated Resistance

ABSTRACT A biotechnological application of artificial microRNAs (amiRs) is the generation of plants that are resistant to virus infection. This resistance has proven to be highly effective and sequence specific. However, before these transgenic plants can be deployed in the field, it is important to evaluate the likelihood of the emergence of resistance-breaking mutants. Two issues are of particular interest: (i) whether such mutants can arise in nontransgenic plants that may act as reservoirs and (ii) whether a suboptimal expression level of the transgene, resulting in subinhibitory concentrations of the amiR, would favor the emergence of escape mutants. To address the first issue, we experimentally evolved independent lineages of Turnip mosaic virus (TuMV) (family Potyviridae) in fully susceptible wild-type Arabidopsis thaliana plants and then simulated the spillover of the evolving virus to fully resistant A. thaliana transgenic plants. To address the second issue, the evolution phase took place with transgenic plants that expressed the amiR at subinhibitory concentrations. Our results show that TuMV populations replicating in susceptible hosts accumulated resistance-breaking alleles that resulted in the overcoming of the resistance of fully resistant plants. The rate at which resistance was broken was 7 times higher for TuMV populations that experienced subinhibitory concentrations of the antiviral amiR. A molecular characterization of escape alleles showed that they all contained at least one nucleotide substitution in the target sequence, generally a transition of the G-to-A and C-to-U types, with many instances of convergent molecular evolution. To better understand the viral population dynamics taking place within each host, as well as to evaluate relevant population genetic parameters, we performed in silico simulations of the experiments. Together, our results contribute to the rational management of amiR-based antiviral resistance in plants.

Replication Mode and Landscape Topology Differentially Affect RNA Virus Mutational Load and Robustness

ABSTRACT Regardless of genome polarity, intermediaries of complementary sense must be synthesized and used as templates for the production of new genomic strands. Depending on whether these new genomic strands become themselves templates for producing extra antigenomic ones, thus giving rise to geometric growth, or only the firstly synthesized antigenomic strands can be used to this end, thus following Lurias stamping machine model, the abundances and distributions of mutant genomes will be different. Here we propose mathematical and bit string models that allow distinguishing between stamping machine and geometric replication. We have observed that, regardless the topology of the fitness landscape, the critical mutation rate at which the master sequence disappears increases as the mechanism of replication switches from purely geometric to stamping machine. We also found that, for a wide range of mutation rates, large-effect mutations do not accumulate regardless the scheme of replication. However, mild mutations accumulate more in the geometric model. Furthermore, at high mutation rates, geometric growth leads to a population collapse for intermediate values of mutational effects at which the stamping machine still produces master genomes. We observed that the critical mutation rate was weakly dependent on the strength of antagonistic epistasis but strongly dependent on synergistic epistasis. In conclusion, we have shown that RNA viruses may increase their robustness against the accumulation of deleterious mutations by replicating as stamping machines and that the magnitude of this benefit depends on the topology of the fitness landscape assumed.

Dynamics of a Plant RNA Virus Intracellular Accumulation: Stamping Machine vs. Geometric Replication

The tremendous evolutionary potential of RNA viruses allows them to thrive despite host defense mechanisms and endows them with properties such as emergence, host switching, and virulence. The frequency of mutant viruses after an infectious process results from the interplay between the error rate of the viral replicase, from purifying mechanisms acting after transcription on aberrant RNAs, and from the amplification dynamics of virus RNA positive (+) and negative (–) strands. Two extreme scenarios describing viral RNA amplification are the geometric growth, in which each RNA strand serves as template for the synthesis of complementary strands with the same efficiency, and the stamping machine, where a strand is reiteratively used as template to synthesize multiple copies of the complementary. The resulting mutation frequencies are completely different, being geometric growth largely more mutagenic than stamping machine. In this work we evaluate the contribution of geometric growth and stamping machine to the overall genome amplification of the plant (+)-strand RNA virus turnip mosaic virus (TuMV). By means of transfection experiments of Nicotiana benthamiana protoplasts with a TuMV cDNA infectious clone and by using strand-specific quantitative real-time PCR, we determined the amplification dynamics of viral (+) and (–) RNA during a single-cell infectious process. A mathematical model describing the amplification of each viral strand was fitted to the data. Analyses of the model parameters showed that TuMV (+) and (–) RNA amplification occurs through a mixed strategy with ∼93% of genomes produced via stamping machine and only ∼7% resulting from geometric growth.

Quasispecies Spatial Models for RNA Viruses with Different Replication Modes and Infection Strategies

Empirical observations and theoretical studies suggest that viruses may use different replication strategies to amplify their genomes, which impact the dynamics of mutation accumulation in viral populations and therefore, their fitness and virulence. Similarly, during natural infections, viruses replicate and infect cells that are rarely in suspension but spatially organized. Surprisingly, most quasispecies models of virus replication have ignored these two phenomena. In order to study these two viral characteristics, we have developed stochastic cellular automata models that simulate two different modes of replication (geometric vs stamping machine) for quasispecies replicating and spreading on a two-dimensional space. Furthermore, we explored these two replication models considering epistatic fitness landscapes (antagonistic vs synergistic) and different scenarios for cell-to-cell spread, one with free superinfection and another with superinfection inhibition. We found that the master sequences for populations replicating geometrically and with antagonistic fitness effects vanished at low critical mutation rates. By contrast, the highest critical mutation rate was observed for populations replicating geometrically but with a synergistic fitness landscape. Our simulations also showed that for stamping machine replication and antagonistic epistasis, a combination that appears to be common among plant viruses, populations further increased their robustness by inhibiting superinfection. We have also shown that the mode of replication strongly influenced the linkage between viral loci, which rapidly reached linkage equilibrium at increasing mutations for geometric replication. We also found that the strategy that minimized the time required to spread over the whole space was the stamping machine with antagonistic epistasis among mutations. Finally, our simulations revealed that the multiplicity of infection fluctuated but generically increased along time.

Error threshold in rna quasispecies models with complementation

A general assumption of quasispecies models of replicons dynamics is that the fitness of a genotype is entirely determined by its sequence. However, a more biologically plausible situation is that fitness depends on the proteins that catalyze metabolic reactions, including replication. In a stirred population of replicons, such as viruses replicating and accumulating within the same cell, the association between a given genome and the proteins it encodes is not tight as it can be replicated by proteins translated from other genomes. We have investigated how this complementation phenomenon affects the error threshold in simple quasispecies mean field models. We first studied a model in which the master and the mutant genomes code for wild-type and mutant replicases, respectively. We assume that the mutant replicase has a reduced activity and that the wild-type replicase does not have increased affinity for the master genome. The whole pool of replicases can bind and replicate both genomes. We then analyze a different model considering a more extreme case of mutant genomes, the defective interfering particles (DIPs) described in many cases of viral infection. DIPs, with a higher replication rate owed to their shorter genomes, do not code for replicase, but they are able of using the replicase translated from the master genome. Our models allow to study how the probability of interaction between the genomes and the whole pool of replicases affects the error threshold. In both systems we characterize the scenario of coexistence between master and mutant genomes, providing the critical values of mutation rate, mu(c), and the critical interaction rate between master genomes and replicases, gamma(c), at which the quasispecies enters into error catastrophe, a situation in which the mutant genomes dominate the population. In both cases, we showed that the error-threshold transition is given by transcritical-like bifurcations, suggesting a continuous phase transition. We have also found that the region in the parameter space (mu,gamma) in which the master sequence survives is reduced when DIPs are introduced into the system.

GHOSTS IN THE ORIGINS OF LIFE

The so-called bottleneck or ghost can appear after a saddle-node bifurcation, leaving a region in phase space by which the flow is attracted although no fixed points are present. Such ghosts, displayed by some dynamical systems, actually cause a delay of the flow. In this paper, we analyze a saddle-node ghost found in a biological model for the two-member hypercycle dynamics. The model predicts a scaling law of the dynamic delay caused by the ghost near the threshold: τ ~ ϕ-1/2, consistent with previous results in physical systems. Possible biological meanings for such a dynamical phenomenon are outlined.

General scaling law in the saddle–node bifurcation: a complex phase space study

Saddle–node bifurcations have been described in a multitude of nonlinear dynamical systems modeling physical, chemical, as well as biological systems. Typically, this type of bifurcation involves the transition of a given set of fixed points from the real to the complex phase space. After the bifurcation, a saddle remnant can continue influencing the flows and generically, for non-degenerate saddle–node bifurcations, the time the flows spend in the bottleneck region of the ghost follows the inverse square root scaling law. Here we analytically derive this scaling law for a general one-dimensional, analytical, autonomous dynamical system undergoing a not necessarily non-degenerate saddle–node bifurcation, in terms of the degree of degeneracy by using complex variable techniques. We then compare the analytic calculations with a one-dimensional equation modeling the dynamics of an autocatalytic replicator. The numerical results are in agreement with the analytical solution.

Can a minimal replicating construct be identified as the embodiment of cancer

Genomic instability is a hallmark of cancer. Cancer cells that exhibit abnormal chromosomes are characteristic of most advanced tumours, despite the potential threat represented by accumulated genetic damage. Carcinogenesis involves a loss of key components of the genetic and signalling molecular networks; hence some authors have suggested that this is part of a trend of cancer cells to behave as simple, minimal replicators. In this study, we explore this conjecture and suggest that, in the case of cancer, genomic instability has an upper limit that is associated with a minimal cancer cell network. Such a network would include (for a given microenvironment) the basic molecular components that allow cells to replicate and respond to selective pressures. However, it would also exhibit internal fragilities that could be exploited by appropriate therapies targeting the DNA repair machinery. The implications of this hypothesis are discussed.

Differences in Accumulation and Virulence Determine the Outcome of Competition during Tobacco etch virus Coinfection

Understanding the evolution of virulence for RNA viruses is essential for developing appropriate control strategies. Although it has been usually assumed that virulence is a consequence of within-host replication of the parasite, viral strains may be highly virulent without experiencing large accumulation as a consequence of immunopathological host responses. Using two strains of Tobacco etch potyvirus (TEV) that show a negative relationship between virulence and accumulation rate, we first explored the evolution of virulence and fitness traits during simple and mixed infections. Short-term evolution experiments initiated with each strain independently confirmed the genetic and evolutionary stability of virulence and viral load, although infectivity significantly increased for both strains. Second, competition experiments between hypo- and hypervirulent TEV strains have shown that the outcome of competition is driven by differences in replication rate. A simple mathematical model has been developed to analyze the dynamics of these two strains during coinfection. The model qualitatively reproduced the experimental results using biologically meaningful parameters. Further analyses of the model also revealed a wide parametric region in which a low-fitness but hypovirulent virus can still outcompete a high-fitness but hypervirulent one. These results provide additional support to the observation that virulence and within-host replication may not necessarily be strongly tied in plant RNA viruses.