RNA Virus Populations Research Articles

Extensive Genome-Wide Variability of Human Cytomegalovirus in Congenitally Infected Infants

Research has shown that RNA virus populations are highly variable, most likely due to low fidelity replication of RNA genomes. It is generally assumed that populations of DNA viruses will be less complex and show reduced variability when compared to RNA viruses. Here, we describe the use of high throughput sequencing for a genome wide study of viral populations from urine samples of neonates with congenital human cytomegalovirus (HCMV) infections. We show that HCMV intrahost genomic variability, both at the nucleotide and amino acid level, is comparable to many RNA viruses, including HIV. Within intrahost populations, we find evidence of selective sweeps that may have resulted from immune-mediated mechanisms. Similarly, genome wide, population genetic analyses suggest that positive selection has contributed to the divergence of the HCMV species from its most recent ancestor. These data provide evidence that HCMV, a virus with a large dsDNA genome, exists as a complex mixture of genome types in humans and offer insights into the evolution of the virus.

Molecular Dynamics Simulations of Viral RNA Polymerases Link Conserved and Correlated Motions of Functional Elements to Fidelity

The viral RNA-dependent RNA polymerase (RdRp) is essential for multiplication of all RNA viruses. The sequence diversity of an RNA virus population contributes to its ability to infect the host. This diversity emanates from errors made by the RdRp during RNA synthesis. The physical basis for RdRp fidelity is unclear but is linked to conformational changes occurring during the nucleotide-addition cycle. To understand RdRp dynamics that might influence RdRp function, we have analyzed all-atom molecular dynamics simulations on the nanosecond timescale of four RdRps from the picornavirus family that exhibit 30–74% sequence identity. Principal component analysis showed that the major motions observed during the simulations derived from conserved structural motifs and regions of known function. The dynamics of residues participating in the same biochemical property, for example, RNA binding, nucleotide binding or catalysis, were correlated even when spatially distant on the RdRp structure. The conserved and correlated dynamics of functional structural elements suggest coevolution of dynamics with structure and function of the RdRp. Crystal structures of all picornavirus RdRps exhibit a template–nascent RNA duplex channel too small to fully accommodate duplex RNA. Simulations revealed opening and closing motions of the RNA and nucleoside triphosphate channels, which might be relevant to nucleoside triphosphate entry, inorganic pyrophosphate exit and translocation. A role for nanosecond timescale dynamics in RdRp fidelity is supported by the altered dynamics of the high-fidelity G64S derivative of PV RdRp relative to wild-type enzyme.

Population Diversity and its Relationship with Infectious and Tumor Diseases

Evolution is based on mutation processes and the subsequent selection of the best-adapted individuals, i.e., “survival of the fittest”. Thus, diversity is a common feature in living things. Pathogens and their hosts, as well as tumoral processes are no exception. Changes in the genetic material of microorganisms allow them to colonize new ecological niches while hosts try to adapt to the new residents by either eliminating them or by establishing mutualistic interactions. The permanent coexistence of microorganisms and their hosts leads to continuous evolvement necessary for their survival, following a strategy that was initially proposed for coevolution models in predator-prey interactions known as the “Red Queen Hypothesis”. Similar conclusions can be established for tumoral processes. Thus, in general, tumors, microorganisms, and hosts can be considered at the population level as complex distributions of related although not identical individuals. Initially, these complex systems were described in a theoretical model that tried to explain the functioning of populations of primitive replicons during the origin of life. This theory was later experimentally confirmed by analyzing populations of positive-strand RNA viruses or RNA(+). In both cases, for replicons and RNA(+) viruses, the populations are described as quasispecies. Here, we will review the most recent discoveries regarding interactions between population diversity and pathogenicity, both in microorganisms and tumors. Keywords: Quasispecies, virus, cancer, Tumor, Pathogens, microorganisms, positive-strand RNA, nucleic acid, nucleotide, Red Queen Hypothesis, Muller's ratchet, Q bacteriophage, foot-and-mouth disease, prokaryotes, eukaryotes, human immunodeficiency virus, hepatitis C virus, HIV, reverse transcriptase, lymphocytic choriomeningitis virus, highly active antiretroviral therapy, HAART, metastasis, human oncogenes, Cellular DNA, mutation, heterogeneity, HAACT, Darwinian model, Pseudomonas aeruginosa, Pf4 bacteriophage

Metagenomic Analysis of RNA Viruses in a Fresh Water Lake

Freshwater lakes and ponds present an ecological interface between humans and a variety of host organisms. They are a habitat for the larval stage of many insects and may serve as a medium for intraspecies and interspecies transmission of viruses such as avian influenza A virus. Furthermore, freshwater bodies are already known repositories for disease-causing viruses such as Norwalk Virus, Coxsackievirus, Echovirus, and Adenovirus. While RNA virus populations have been studied in marine environments, to this date there has been very limited analysis of the viral community in freshwater. Here we present a survey of RNA viruses in Lake Needwood, a freshwater lake in Maryland, USA. Our results indicate that just as in studies of other aquatic environments, the majority of nucleic acid sequences recovered did not show any significant similarity to known sequences. The remaining sequences are mainly from viral types with significant similarity to approximately 30 viral families. We speculate that these novel viruses may infect a variety of hosts including plants, insects, fish, domestic animals and humans. Among these viruses we have discovered a previously unknown dsRNA virus closely related to Banna Virus which is responsible for a febrile illness and is endemic to Southeast Asia. Moreover we found multiple viral sequences distantly related to Israeli Acute Paralysis virus which has been implicated in honeybee colony collapse disorder. Our data suggests that due to their direct contact with humans, domestic and wild animals, freshwater ecosystems might serve as repositories of a wide range of viruses (both pathogenic and non-pathogenic) and possibly be involved in the spread of emerging and pandemic diseases.

Detecting natural selection in RNA virus populations using sequence summary statistics

At present, most analyses that aim to detect the action of natural selection upon viral gene sequences use phylogenetic estimates of the ratio of silent to replacement mutations. Such methods, however, are impractical to compute on large data sets comprising hundreds of complete viral genomes, which are becoming increasingly common due to advances in genome sequencing technology. Here we investigate the statistical performance of computationally efficient tests that are based on sequence summary statistics, and explore their applicability to RNA virus data sets in two ways. Firstly, we perform extensive simulations in order to measure the type I error of two well-known summary statistic methods – Tajima's D and the McDonald–Kreitman test – under a range of virus-like mutational and demographic scenarios. Secondly, we apply these methods to a compilation of ∼100 RNA virus alignments that represent natural RNA virus populations. In addition, we develop and introduce a new implementation of the McDonald–Kreitman test and show that it greatly improves the test's statistical reliability on typical viral data sets. Our results suggest that variants of the McDonald–Kreitman test could prove useful in the analysis of very large sets of highly diverse viral genetic data.

Complementation and Epistasis in Viral Coinfection Dynamics

Coinfection in RNA virus populations results in two important phenomena, complementation and recombination. Of the two, complementation has a strong effect on selection against deleterious mutations, as has been confirmed in earlier studies. As complementation delays the purging of less-fit mutations, coinfection may be detrimental to the evolution of a virus population. Here we employ both deterministic modeling and stochastic simulation to explore the mechanisms underlying the interactions between complementation and other evolutionary factors, namely, mutation, selection, and epistasis. We find that strong complementation reduces slightly the overall fitness of a virus population but substantially enhances its diversity and robustness, especially when interacting with selection and epistasis.

The 30th anniversary of quasispecies

The meeting on ‘Quasispecies: past, present and future’ took place between 17 and 18 November 2008, in Barcelona, Spain, and was organized by J. Gomez, C. Lopez‐Galindez, M.A. Martinez & A. Mas. ![][1] A meeting was held in Barcelona, Spain, in November 2008 to celebrate the 30th anniversary of the publication of the article that described the extensive genetic heterogeneity of bacteriophage Qβ (Domingo et al , 1978), which is considered to mark the beginning of experimental studies on viral quasispecies. This meeting was held at the impressive fifteenth century building of the ancient Hospital de la Santa Creu in the old town of Barcelona, which is now the headquarters of the Institut d'Estudis Catalans, and was attended by C. Weissmann (Jupiter, FL, USA), M. Billeter (Zurich, Switzerland) and E. Domingo (Madrid, Spain), who were three early protagonists of the phage Qβ work at the University of Zurich in the 1970s. Several speakers presented their results on the theoretical aspects of the population dynamics of cells and viruses, the clinical implications of quasispecies, and extensions of the quasispecies concept to cellular genes and prions. The meeting was introduced by A. Mas (Albacete, Spain), who reflected on the increasing impact that quasispecies have had in the scientific literature over the past three decades, and quoted some of the key references on viral quasispecies (Martell et al , 1992; Meyerhans et al , 1989; Najera et al , 1995; Vignuzzi et al , 2006; for a historical review of the impact of quasispecies in virology, see Holland, 2006). The scientific presentations were opened by Weissmann and Domingo, who were the last and first authors of the 1978 paper, respectively. Their talks conveyed the scientific atmosphere of the 1970s—when molecular biology was carried out with few recombinant‐DNA techniques—to a young audience. Nucleic‐acid sequencing was in … [1]: /embed/graphic-1.gif

Rapid Sequence Change and Geographical Spread of Human Parvovirus B19: Comparison of B19 Virus Evolution in Acute and Persistent Infections

Parvovirus B19 is a common human pathogen maintained by horizontal transmission between acutely infected individuals. However, B19 virus can also be detected in tissues throughout the life of the host, although little is understood about the nature of such persistence. In the current study, we created large VP1/2 sequence data sets of plasma- and tissue (autopsy)-derived variants of B19 virus with known sample dates to compare the rates of sequence change in exogenous virus populations with those in persistently infected individuals. By using linear regression and likelihood-based methods (such as the BEAST program), we found that plasma-derived B19 virus showed a substitution rate of 4 x 10(-4) and an unconstrained (synonymous)-substitution rate of 18 x 10(-4) per site per year, several times higher than previously estimated and within the range of values for mammalian RNA viruses. The underlying high mutation frequency implied by these substitution rates may enable rapid adaptive changes that are more commonly ascribed to RNA virus populations. These revised estimates predict that the last common ancestor for currently circulating genotype 1 variants of B19 virus existed around 1956 to 1959, fitting well with previous analyses of the B19 virus "bioportfolio" that support a complete cessation of genotype 2 infections and their replacement by genotype 1 infections in the 1960s. In contrast, the evolution of B19 virus amplified from tissue samples was best modeled by using estimated dates of primary infection rather than sample dates, consistent with slow or absent sequence change during persistence. Determining what epidemiological or biological factors led to such a complete and geographically extensive population replacement over this short period is central to further understanding the nature of parvovirus evolution.

Effects of Spatial Competition on the Diversity of a Quasispecies

The diversity harbored by populations of RNA viruses results from high mutation rates, as well as from the characteristics of the environment where they evolve. By means of a simple model for structured quasispecies, we quantify how competition for space among phenotypic types shapes their distribution at the mutation-selection equilibrium. We introduce a general framework to treat this problem and relate mutation rate and competition strength to the quasispecies composition. For diffusion limited competition, diversity typically increases and the asymptotic growth rate of the population diminishes as diffusion decreases. Limited mobility confers a relative advantage to worse competitors. The stationary state is characterized by an over-production of viral particles. Empirical data allow an estimation of mutation rates compatible with the diversity observed in viral populations infecting cellular monolayers.

Evolutionary Genomics of Host Adaptation in Vesicular Stomatitis Virus

Populations experiencing similar selection pressures can sometimes diverge in the genetic architectures underlying evolved complex traits. We used RNA virus populations of large size and high mutation rate to study the impact of historical environment on genome evolution, thus increasing our ability to detect repeatable patterns in the evolution of genetic architecture. Experimental vesicular stomatitis virus populations were evolved on HeLa cells, on MDCK cells, or on alternating hosts. Turner and Elena (2000. Cost of host radiation in an RNA virus. Genetics. 156:1465-1470.) previously showed that virus populations evolved in single-host environments achieved high fitness on their selected hosts but failed to increase in fitness relative to their ancestor on the unselected host and that alternating-host-evolved populations had high fitness on both hosts. Here we determined the complete consensus sequence for each evolved population after 95 generations to gauge whether the parallel phenotypic changes were associated with parallel genomic changes. We also analyzed the patterns of allele substitutions to discern whether differences in fitness across hosts arose through true pleiotropy or the presence of not only a mutation that is beneficial in both hosts but also 1 or more mutations at other loci that are costly in the unselected environment (mutation accumulation [MA]). We found that ecological history may influence to what extent pleiotropy and MA contribute to fitness asymmetries across environments. We discuss the degree to which current genetic architecture is expected to constrain future evolution of complex traits, such as host use by RNA viruses.

Fitness Declines in Tobacco Etch Virus upon Serial Bottleneck Transfers

It has been well established that populations of RNA viruses transmitted throughout serial bottlenecks suffer from significant fitness declines as a consequence of the accumulation of deleterious mutations by the onset of Muller's ratchet. Bottlenecks are unavoidably linked to different steps of the infectious cycle of most plant RNA viruses, such as vector-mediated transmissions and systemic colonization of new leaves. Here we report evidence for fitness declines by the accumulation of deleterious mutations in the potyvirus Tobacco etch virus (TEV). TEV was inoculated into the nonsystemic host Chenopodium quinoa, and local lesions were isolated and used to initiate 20 independent mutation accumulation lineages. Weekly, a random lesion from each lineage was isolated and used to inoculate the next set of plants. At each transfer, the Malthusian growth rate was estimated. After 11 consecutive transfers, all lineages suffered significant fitness losses, and one even became extinct. The average rate of fitness decline was 5% per day. The average pattern of fitness decline was consistent with antagonistic epistasis between deleterious mutations, as postulated for antiredundant genomes. Temporal fitness fluctuations were not explained by random noise but reflected more complex underlying processes related to emergence and self-organization phenomena.

Phylogenetic Evidence for Deleterious Mutation Load in RNA Viruses and Its Contribution to Viral Evolution

Populations of RNA viruses are often characterized by abundant genetic variation. However, the relative fitness of these mutations is largely unknown, although this information is central to our understanding of viral emergence, immune evasion, and drug resistance. Here we develop a phylogenetic method, based on the distribution of nonsynonymous and synonymous changes, to assess the relative fitness of polymorphisms in the structural genes of 143 RNA viruses. This reveals that a substantial proportion of the amino acid variation observed in natural populations of RNA viruses comprises transient deleterious mutations that are later purged by purifying selection, potentially limiting virus adaptability. We also demonstrate, for the first time, the existence of a relationship between amino acid variability and the phylogenetic distribution of polymorphisms. From this relationship, we propose an empirical threshold for the maximum viable deleterious mutation load in RNA viruses.

Emerging epidemic diseases of frogs in Britain are dependent on the source of ranavirus agent and the route of exposure

A series of transmission studies was conducted to investigate the aetiology, or aetiologies, of emerging fatal epidemic disease syndromes affecting the common frog (Rana temporaria) in Britain. The syndromes, characterized by skin ulceration or systemic haemorrhages, were induced upon exposure to lesion homogenates or cultured ranavirus. The re-isolation of ranavirus from experimentally affected frogs fulfilled Koch's postulates. Aeromonas hydrophila, previously associated with similar lesions, was not significant to disease development. Unexpectedly, disease outcomes were influenced by both the source of agent and the route of exposure, indicating that different ranaviruses with different tissue tropisms and pathogeneses (possibly similar to quasi-species in RNA virus populations) are circulating in the British common frog population. Our findings confirm that ranavirus disease has emerged as an important cause of amphibian mortality in Britain.

A Phylogenetic Method for Detecting Positive Epistasis in Gene Sequences and Its Application to RNA Virus Evolution

RNA virus genomes are compact, often containing multiple overlapping reading frames and functional secondary structure. Consequently, it is thought that evolutionary interactions between nucleotide sites are commonplace in the genomes of these infectious agents. However, the role of epistasis in natural populations of RNA viruses remains unclear. To investigate the pervasiveness of epistasis in RNA viruses, we used a parsimony-based computational method to identify pairs of co-occurring mutations along phylogenies of 177 RNA virus genes. This analysis revealed widespread evidence for positive epistatic interactions at both synonymous and nonsynonymous nucleotide sites and in both clonal and recombining viruses, with the majority of these interactions spanning very short sequence regions. These findings have important implications for understanding the key aspects of RNA virus evolution, including the dynamics of adaptation. Additionally, many comparative analyses that utilize the phylogenetic relationships among gene sequences assume that mutations represent independent, uncorrelated events. Our results show that this assumption may often be invalid.

Quasispecies diversity determines pathogenesis through cooperative interactions in a viral population

An RNA virus population does not consist of a single genotype; rather, it is an ensemble of related sequences, termed quasispecies. Quasispecies arise from rapid genomic evolution powered by the high mutation rate of RNA viral replication. Although a high mutation rate is dangerous for a virus because it results in nonviable individuals, it has been hypothesized that high mutation rates create a 'cloud' of potentially beneficial mutations at the population level, which afford the viral quasispecies a greater probability to evolve and adapt to new environments and challenges during infection. Mathematical models predict that viral quasispecies are not simply a collection of diverse mutants but a group of interactive variants, which together contribute to the characteristics of the population. According to this view, viral populations, rather than individual variants, are the target of evolutionary selection. Here we test this hypothesis by examining the consequences of limiting genomic diversity on viral populations. We find that poliovirus carrying a high-fidelity polymerase replicates at wild-type levels but generates less genomic diversity and is unable to adapt to adverse growth conditions. In infected animals, the reduced viral diversity leads to loss of neurotropism and an attenuated pathogenic phenotype. Notably, using chemical mutagenesis to expand quasispecies diversity of the high-fidelity virus before infection restores neurotropism and pathogenesis. Analysis of viruses isolated from brain provides direct evidence for complementation between members in the quasispecies, indicating that selection indeed occurs at the population level rather than on individual variants. Our study provides direct evidence for a fundamental prediction of the quasispecies theory and establishes a link between mutation rate, population dynamics and pathogenesis.

Remote Site Control of an Active Site Fidelity Checkpoint in a ViralRNA-dependent RNAPolymerase

The kinetic, thermodynamic, and structural basis for fidelity of nucleic acid polymerases remains controversial. An understanding of viral RNA-dependent RNA polymerase (RdRp) fidelity has become a topic of considerable interest as a result of recent experiments that show that a 2-fold increase in fidelity attenuates viral pathogenesis and a 2-fold decrease in fidelity reduces viral fitness. Here we show that a conformational change step preceding phosphoryl transfer is a key fidelity checkpoint for the poliovirus RdRp (3Dpol). We provide evidence that this conformational change step is orientation of the triphosphate into a conformation suitable for catalysis, suggesting a kinetic and structural model for RdRp fidelity that can be extrapolated to other classes of nucleic acid polymerases. Finally, we show that a site remote from the catalytic center can control this checkpoint, which occurs at the active site. Importantly, similar connections between a remote site and the active site exist in a wide variety of viral RdRps. The capacity for sites remote from the catalytic center to alter fidelity suggests new possibilities for targeting the viral RdRp for antiviral drug development.

Resistance is futile

A new study of dominant negative functions in cells infected by a positive-strand RNA virus detects an array of locus- and allele-specific effects. Exploiting subunit defects in multi-component complex assemblies provides a new approach to identifying targets for antiviral therapies that may inhibit the emergence of drug-resistant RNA virus populations.

Following the very initial growth of biological RNA viral clones

Due to their extremely high genetic diversity, which is a direct consequence of high mutation rates, RNA viruses are often described as molecular quasispecies. According to this theory, RNA virus populations cannot be understood in terms of individual viral clones, as they are clouds of interconnected mutants, but this prediction has not yet been demonstrated experimentally. The goal of this study was to determine the fitness of individual clones sampled from a given RNA virus population, a necessary previous step to test the above prediction. To do so, limiting dilutions of a vesicular stomatitis virus population were employed to isolate single viral clones and their initial growth dynamics were followed, corresponding to the release of the first few hundred viral particles. This technique is useful for estimating basic fitness parameters, such as intracellular growth rate, viral yield per cell, rate at which cells are infected and time spent in cell-to-cell transmission. A combination of these parameters allows estimation of the fitness of individual clones, which seems to be determined mainly by their ability to complete infection cycles more quickly. Interestingly, fitness was systematically higher for initial clones than for their derived populations. In addition to environmental changes, such as cellular defence mechanisms, these differences are attributable to high RNA virus mutation rates.

Ribavirin and lethal mutagenesis of poliovirus: molecular mechanisms, resistance and biological implications

Positive strand RNA virus populations are a collection of similar but genetically different viruses. They exist as viral quasispecies due to the high mutation rates of the low fidelity viral RNA-dependent RNA polymerase (RdRp). It is thought that this genomic heterogeneity is advantageous to the population, allowing for adaptation to rapidly changing environments that present varying types and degrees of selective pressure. However, one consequence of this extensive diversity is the susceptibility to mutagens that further increase sequence variation. Since RNA viruses live at the edge of maximal variability, an increase in the mutation rate is likely to force the virus beyond the tolerable mutation frequency into ‘error catastrophe’. One such mutagen, ribavirin, is an antiviral nucleoside analog that is mutagenic to several RNA viruses. Ribavirin is incorporated into the viral genome causing lethal mutagenesis and a subsequent decrease in the specific infectivity. Even so, passaging poliovirus in the presence of low to intermediate concentrations of the drug leads to the emergence of a viral population resistant to the effects of ribavirin. These viruses have a point mutation in the RdRp that increases the overall polymerase fidelity. Interestingly, as predicted by the quasispecies theory, ribavirin resistant viruses are less adaptable, as they are more susceptible to other non-mutagenic antiviral drugs and are highly attenuated in vivo. Here, we review the mechanism of action of ribavirin on poliovirus and other RNA viruses, the possibility for escape via increased fidelity of the viral polymerase, the consequences of this response on viral population dynamics, and the biological implications for the therapeutic use of mutagenic antiviral agents.

Do Taq-generated RT-PCR products from RNA viruses accurately reflect viral genetic heterogeneity?

Since the first report of genetically heterogeneous, or quasispecies, populations of RNA viruses, the genetic heterogeneity of the RNA genomes of major viral pathogens has been extensively studied. These studies aim to provide insights into the evolutionary pressures that act upon viruses, in order to define windows where anti-viral therapies will be most effective, to take prognostic values from viral genetic distributions at a given time, and to aid the development of novel therapeutic compounds that may tilt viral replication towards information loss. Many methodologies are employed to analyse genetic distributions of a virus in a given sample, but all involve the generation, and subsequent analysis, of the sequence information contained in a reverse-transcription-polymerase chain reaction (RT-PCR) product. Despite the fact that the aim of these RT-PCRs is to obtain sequence information from viral genomes, their application to this task is approached without adequate consideration of this end-goal. The establishment of an RT-PCR for a specific viral target genome generally proceeds in the same fashion as one would apply to establishing a PCR to determine the presence or absence of a specific target sequence in a given sample. However, it is becoming increasingly apparent that RT-PCR products generated by amplification with the ubiquitous thermostable DNA polymerase Taq, coupled with standard cloning and sequencing methodologies, has the potential to yield inaccurate and misleading data as pertains to the information content of populations of RNA viral genomes. This review discusses varying approaches employed to analyse heterogeneous populations of hepatitis C virus RNA genomes.