Duration and fitness dependence of quasispecies memory.

The duration and fitness dependence of memory in viral quasispecies evolving in cell culture have been investigated using two genetic markers of foot-and-mouth disease virus (FMDV). In lineages of antigenic variant FMDV RED, which reverted to FMDV RGD, memory FMDV RED genomes were detected after 50 infectious cycles, and memory level was fitness dependent. In growth-competition experiments between a reference FMDV RGD and two different FMDV RED populations, a 7.6-fold higher fitness of the initial FMDV RED population resulted in 30 to 100-fold higher memory level. In lineages of low-fitness clones containing an elongated internal polyadenylate tract, revertants lacking excess adenylate residues became dominant by passage 20. However, genomes including a larger number of adenylate residues were detected as memory genomes after at least 150 infectious cycles. Thus, quasispecies memory can be durable and is fitness dependent, as predicted from the growth competition of two mutant forms of a genome. An understanding of factors influencing quasispecies memory levels and duration may have implications for the extended diagnosis of viruses based on the quantification of minority genomes.     

Memory in viral quasispecies

Biological adaptive systems share some common features: variation among their constituent elements and continuity of core information. Some of them, such as the immune system, are endowed with memory of past events. In this study we provide direct evidence that evolving viral quasispecies possess a molecular memory in the form of minority components that populate their mutant spectra. The experiments have involved foot-and-mouth disease virus populations with known evolutionary histories. The composition and behavior of the viral population in response to a selective constraint were influenced by past evolutionary history in a way that could not be predicted from examination of consensus nucleotide sequences of the viral populations. The molecular memory of the viral quasispecies influenced both the nature and the intensity of the response of the virus to a selective constraint.     

Fitness increase of memory genomes in a viral quasispecies

Viral quasispecies may contain a subset of minority genomes that reflect those genomic sequences that were dominant at an early phase of quasispecies evolution. Such minority genomes are referred to as memory in viral quasispecies. A memory marker previously characterized in foot-and-mouth disease virus (FMDV) is an internal oligoadenylate tract of variable length that became dominant upon serial plaque-to-plaque transfers of FMDV clones. During large population passages, genomes with internal oligoadenylate were outcompeted by wild-type revertants but remained in the mutant spectra as memory genomes. Here, we report a quantification of relative fitness of several FMDV clones, harboring internal oligoadenylate tracts of different length, and that were retrieved at early or late times (passage number) after implementation of memory. The results show that for any given length range of the oligoadenylate, maintenance in memory resulted in an increase in relative fitness, comparable to the increase undergone by the entire population. The fitness increase is in agreement with the Red Queen hypothesis, and implies a replicative memory mechanism. Thus, permanence of memory genomes may be a source of high fitness variants despite their initial low fitness, and despite having remained hidden in mutant spectra. This reinforces the interest of diagnosing minority genomes during chronic human and animal viral infections.     

The quasispecies (extremely heterogeneous) nature of viral RNA genome populations: biological relevance--a review

We review evidence that cloned (or uncloned) populations of most RNA viruses do not consist of a single genome species of defined sequence, but rather of heterogeneous mixtures of related genomes (quasispecies). Due to very high mutation rates, genomes of a quasispecies virus population share a consensus sequence but differ from each other and from the consensus sequence by one, several, or many mutations. Viral genome analyses by sequencing, fingerprinting, cDNA cloning etc. indicate that most viral RNA populations (quasispecies) contain all possible single and double genomic site mutations and varying proportions of triple, quadruple, etc. site mutations. This quasispecies structure of RNA virus populations has many important theoretical and practical implications because mutations at only one or a few sites may alter the phenotype of an RNA virus.     

Memory in retroviral quasispecies: experimental evidence and theoretical model for human immunodeficiency virus

Viral quasispecies may possess a molecular memory of their past evolutionary history, imprinted on minority components of the mutant spectrum. Here we report experimental evidence and a theoretical model for memory in retroviral quasispecies in vivo. Apart from replicative memory associated with quasispecies dynamics, retroviruses may harbour a "cellular" or "anatomical" memory derived from their integrative cycle and the presence of viral reservoirs in body compartments. Three independent sets of data exemplify the two kinds of memory in human immunodeficiency virus type 1 (HIV-1). The data provide evidence of re-emergence of sequences that were hidden in cellular or anatomical compartments for extended periods of infection, and recovery of a quasispecies from pre-existing genomes. We develop a three-component model that incorporates the essential features of the quasispecies dynamics of retroviruses exposed to selective pressures. Significantly, a numerical study based on this model is in agreement with the experimental data, further supporting the existence of both replicative and reservoir memory in retroviral quasispecies.     

Influence of the mutant spectrum in viral evolution: focused selection of antigenic variants in a reconstructed viral quasispecies

RNA viruses replicate as complex mutant distributions termed viral quasispecies. Despite this, studies on virus populations subjected to positive selection have generally been performed and analyzed as if the viral population consisted of a defined genomic nucleotide sequence; such a simplification may not reflect accurately the molecular events underlying the selection process. In the present study, we have reconstructed a foot-and-mouth disease virus quasispecies with multiple, low-frequency, genetically distinguishable mutants that can escape neutralization by a monoclonal antibody. Some of the mutants included an amino acid substitution that affected an integrin recognition motif that overlaps with the antibody-binding site, whereas other mutants included an amino acid substitution that affected antibody binding but not integrin recognition. We have monitored consensus and clonal nucleotide sequences of populations passaged either in the absence or the presence of the neutralizing antibody. In both cases, the populations focused toward a specific mutant that was surrounded by a cloud of mutants with different antigenic and cell recognition specificities. In the absence of antibody selection, an antigenic variant that maintained integrin recognition became dominant, but the mutant cloud included as one of its minority components a variant with altered integrin recognition. Conversely, in the presence of antibody selection, a variant with altered integrin recognition motif became dominant, but it was surrounded by a cloud of antigenic variants that maintained integrin recognition. The results have documented that a mutant spectrum can exert an influence on a viral population subjected to a sustained positive selection pressure and have unveiled a mechanism of antigenic flexibility in viral populations, consisting in the presence in the selected quasispecies of mutants with different antigenic and cell recognition specificities.     

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, μ_c, and the critical interaction rate between master genomes and replicases, γ_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 (μ,γ) in which the master sequence survives is reduced when DIPs are introduced into the system.     

Combinatorial analysis and algorithms for quasispecies reconstruction using next-generation sequencing

Background  

Next-generation sequencing (NGS) offers a unique opportunity for high-throughput genomics and has potential to replace Sanger sequencing in many fields, including de-novo sequencing, re-sequencing, meta-genomics, and characterisation of infectious pathogens, such as viral quasispecies. Although methodologies and software for whole genome assembly and genome variation analysis have been developed and refined for NGS data, reconstructing a viral quasispecies using NGS data remains a challenge. This application would be useful for analysing intra-host evolutionary pathways in relation to immune responses and antiretroviral therapy exposures. Here we introduce a set of formulae for the combinatorial analysis of a quasispecies, given a NGS re-sequencing experiment and an algorithm for quasispecies reconstruction. We require that sequenced fragments are aligned against a reference genome, and that the reference genome is partitioned into a set of sliding windows (amplicons). The reconstruction algorithm is based on combinations of multinomial distributions and is designed to minimise the reconstruction of false variants, called in-silico recombinants.     

A trade-off between neutrality and adaptability limits the optimization of viral quasispecies

Theoretical studies of quasispecies usually focus on two properties of those populations at the mutation-selection equilibrium, namely asymptotic growth rate and population diversity. It has been postulated that, as a consequence of the high error rate of quasispecies replication, an increase of neutrality facilitates population optimization by reducing the amount of mutations with a deleterious effect on fitness. In this study we analyse how the optimization of equilibrium properties is affected when a quasispecies evolves in an environment perturbed through frequent bottleneck events. By means of a simple model we demonstrate that high neutrality may be detrimental when the population has to overcome repeated reductions in the population size, and that the property to be optimized in this situation is the time required to regenerate the quasispecies, i.e. its adaptability. In the scenario described, neutrality and adaptability cannot be simultaneously optimized. When fitness is equated with long-term survivability, high neutrality is the appropriate strategy in constant environments, while populations evolving in fluctuating environments are fitter when their neutrality is low, such that they can respond faster to perturbations. Our results might be relevant to better comprehend how a minority virus could displace the circulating quasispecies, a fact observed in natural infections and essential in viral evolution.     

Dynamics of mutation and recombination in a replicating population of complementing, defective viral genomes

In a previous study, we documented that serial passage of a biological clone of foot-and-mouth disease virus (FMDV) at high multiplicity of infection (moi) in cell culture resulted in viral populations dominated by defective genomes that included internal in-frame deletions, affecting the L and capsid-coding regions, and were infectious by complementation. In the present study, analyses of the defective genomes present in individual viral plaques, and of consensus nucleotide sequences determined for the entire genomes of sequential samples, have revealed a continuous dynamics of mutation and recombination. At some points of high genetic instability, multiple minority genomes with different internal deletions co-existed in the population. At later passages, a new defective RNA arose and displaced a related, previously dominant RNA. Nucleotide sequences of the different genomic forms found in sequential isolates have revealed an accumulation of mutations at an average rate of 0.12 substitutions per genome per passage. At the regions around the deletion sites, substantial, minor or no nucleotide sequence identity is found, suggesting relaxed sequence requirements for the occurrence of internal deletions. Competition experiments indicate a selective advantage of late phase defective genomes over their precursor forms. The defective genome-based FMDV retained an expansion of host cell tropism, undergone by the standard virus at a previous stage of the same evolutionary lineage. Thus, despite a complex dynamics of mutation and recombination, and phases of high genetic instability, a biologically relevant phenotypic trait was stably maintained after the evolutionary transition towards a primitive genome segmentation. The results extend the concept of a complex spectrum of mutant genomes to a complex spectrum of defective genomes in some evolutionary transitions of RNA viruses.     

Probability of fixation of an advantageous mutant in a viral quasispecies

The probability that an advantageous mutant rises to fixation in a viral quasispecies is investigated in the framework of multitype branching processes. Whether fixation is possible depends on the overall growth rate of the quasispecies that will form if invasion is successful rather than on the individual fitness of the invading mutant. The exact fixation probability can be calculated only if the fitnesses of all potential members of the invading quasispecies are known. Quasispecies fixation has two important characteristics: First, a sequence with negative selection coefficient has a positive fixation probability as long as it has the potential to grow into a quasispecies with an overall growth rate that exceeds that of the established quasispecies. Second, the fixation probabilities of sequences with identical fitnesses can nevertheless vary over many orders of magnitudes. Two approximations for the probability of fixation are introduced. Both approximations require only partial knowledge about the potential members of the invading quasispecies. The performance of these two approximations is compared to the exact fixation probability on a network of RNA sequences with identical secondary structure.     

Antigenic stimulation specifically reactivates the replication of archived simian immunodeficiency virus genomes in chronically infected macaques

Human immunodeficiency virus/simian immunodeficiency virus (SIV) diversification is a direct consequence of viral replication and occurs principally in secondary lymphoid organs where CD4(+) T cells are activated and proliferate. However, the evolution of viral quasispecies may also be driven by various nonexclusive mechanisms, including adaptation to specific immune responses and modification of viral fitness. Analysis of viral quasispecies in SIV-infected macaques subjected to repeated antigenic stimulations allowed us to demonstrate transient expansions of SIV populations that were highly dependent upon activation of antigen-specific T cells. T-cell clones expanded in response to a particular antigen were infected by a specific viral population and persisted for prolonged periods. Upon a second stimulation by encounter with the same antigen, these specific genomes were at the origin of a new burst of replication, leading to rapid but transient replacement of the viral quasispecies in blood. Finally, longitudinal analysis of SIV sequence variation during and between antigenic stimulations revealed that viral evolution is mostly constrained to periods of strong immunological activity.     

Replicative homeostasis II: Influence of polymerase fidelity on RNA virus quasispecies biology: Implications for immune recognition, viral autoimmunity and other "virus receptor" diseases

Much of the worlds' population is in active or imminent danger from established infectious pathogens, while sporadic and pandemic infections by these and emerging agents threaten everyone. RNA polymerases (RNA_pol) generate enormous genetic and consequent antigenic heterogeneity permitting both viruses and cellular pathogens to evade host defences. Thus, RNA_pol causes more morbidity and premature mortality than any other molecule. The extraordinary genetic heterogeneity defining viral quasispecies results from RNA_pol infidelity causing rapid cumulative genomic RNA mutation a process that, if uncontrolled, would cause catastrophic loss of sequence integrity and inexorable quasispecies extinction. Selective replication and replicative homeostasis, an epicyclical regulatory mechanism dynamically linking RNApol fidelity and processivity with quasispecies phenotypic diversity, modulating polymerase fidelity and, hence, controlling quasispecies behaviour, prevents this happening and also mediates immune escape. Perhaps more importantly, ineluctable generation of broad phenotypic diversity after viral RNA is translated to protein quasispecies suggests a mechanism of disease that specifically targets, and functionally disrupts, the host cell surface molecules – including hormone, lipid, cell signalling or neurotransmitter receptors – that viruses co-opt for cell entry. This mechanism – "Viral Receptor Disease (VRD)" – may explain so-called "viral autoimmunity", some classical autoimmune disorders and other diseases, including type II diabetes mellitus, and some forms of obesity. Viral receptor disease is a unifying hypothesis that may also explain some diseases with well-established, but multi-factorial and apparently unrelated aetiologies – like coronary artery and other vascular diseases – in addition to diseases like schizophrenia that are poorly understood and lack plausible, coherent, pathogenic explanations.     

Semiconservative quasispecies equations for polysomic genomes: the haploid case

This paper develops the semiconservative quasispecies equations for genomes consisting of an arbitrary number of chromosomes. We assume that the chromosomes are distinguishable, so that we are effectively considering haploid genomes. We derive the quasispecies equations under the assumption of arbitrary lesion repair efficiency, and consider the cases of both random and immortal strand chromosome segregation. We solve the model in the limit of infinite sequence length for the case of the static single fitness peak landscape, where the master genome has a first-order growth rate constant of k>1, and all other genomes have a first-order growth rate constant of 1. If we assume that each chromosome can tolerate an arbitrary number of lesions, so that only one master copy of the strands is necessary for a functional chromosome, then for random chromosome segregation we obtain an equilibrium mean fitness of [equation in text] below the error catastrophe, while for immortal strand co-segregation we obtain kappa (t=infinity)=k[e(-mu(1-lambda/2))+e(-mulambda/2)-1] (N denotes the number of chromosomes, lambda denotes the lesion repair efficiency, and mu is identical with epsilonL, where epsilon is the per base-pair mismatch probability, and L is the total genome length). It follows that immortal strand co-segregation leads to significantly better preservation of the master genome than random segregation when lesion repair is imperfect. Based on this result, we conjecture that certain classes of tumor cells exhibit immortal strand co-segregation.     

Fitness alteration of foot-and-mouth disease virus mutants: measurement of adaptability of viral quasispecies.

We document the rapid alteration of fitness of two foot-and-mouth disease virus (FMDV) mutants resistant to a neutralizing monoclonal antibody. Both mutants showed a selective disadvantage in BHK-21 cells when passaged in competition with their parental FMDV. Upon repeated replication of the mutants alone, they acquired a selective advantage over the parental FMDV and fixed additional genomic substitutions without reversion of the monoclonal antibody-resistant phenotype. Thus, variants that were previously kept at low frequency in the mutant spectrum of a viral quasispecies rapidly became the master sequence of a new genomic distribution and dominated the viral population.     

QSdpR: Viral quasispecies reconstruction via correlation clustering

RNA viruses are characterized by high mutation rates that give rise to populations of closely related genomes, known as viral quasispecies. Underlying heterogeneity enables the quasispecies to adapt to changing conditions and proliferate over the course of an infection. Determining genetic diversity of a virus (i.e., inferring haplotypes and their proportions in the population) is essential for understanding its mutation patterns, and for effective drug developments. Here, we present QSdpR, a method and software for the reconstruction of quasispecies from short sequencing reads. The reconstruction is achieved by solving a correlation clustering problem on a read-similarity graph and the results of the clustering are used to estimate frequencies of sub-species; the number of sub-species is determined using pseudo F index. Extensive tests on both synthetic datasets and experimental HIV-1 and Zika virus data demonstrate that QSdpR compares favorably to existing methods in terms of various performance metrics.