A new implication of quasispecies dynamics: Broad virus diversification in absence of external perturbations

Genetic variations and gene expression profiles of Rice Black-streaked dwarf virus (RBSDV) in different host plants and insect vectors: insights from RNA-Seq analysis

Rice black-streaked dwarf virus (RBSDV) is an etiological agent of a destructive disease infecting some economically important crops from the Gramineae family in Asia. While RBSDV causes high yield losses, genetic characteristics of replicative viral populations have not been investigated within different host plants and insect vectors. Herein, eleven publicly available RNA-Seq datasets from Chinese RBSDV-infected rice, maize, and viruliferous planthopper (Laodelphax striatellus) were obtained from the NCBI database. The patterns of SNP and RNA expression profiles of expected RBSDV populations were analyzed by CLC Workbench 20 and Geneious Prime software. These analyses discovered 2,646 mutations with codon changes in RBSDV whole transcriptome and forty-seven co-mutated hotspots with high variant frequency within the crucial regions of S5-1, S5-2, S6, S7-1, S7-2, S9, and S10 open reading frames (ORFs) which are responsible for some virulence and host range functions. Moreover, three joint mutations are located on the three-dimensional protein of P9-1. The infected RBSDV-susceptible rice cultivar KTWYJ3 and indigenous planthopper datasets showed more co-mutated hotspot numbers than others. Our analyses showed the expression patterns of viral genomic fragments varied depending on the host type. Unlike planthopper, S5-1, S2, S6, and S9-1 ORFs, respectively had the greatest read numbers in host plants; and S5-2, S9-2, and S7-2 were expressed in the lowest level. These findings underscore virus/host complexes are effective in the genetic variations and gene expression profiles of plant viruses. Our analysis revealed no evidence of recombination events. Interestingly, the negative selection was observed at 12 RBSDV ORFs, except for position 1015 in the P1 protein, where a positive selection was detected. The research highlights the potential of SRA datasets for analysis of the virus cycle and enhances our understanding of RBSDV's genetic diversity and host specificity.

mRNA Technology and Mucosal Immunization

Current mRNA vaccines are mainly administered via intramuscular injection, which induces good systemic immunity but limited mucosal immunity. Achieving mucosal immunity through mRNA vaccination could diminish pathogen replication at the entry site and reduce interhuman transmission. However, delivering mRNA vaccines to mucosae faces challenges like mRNA degradation, poor entry into cells, and reactogenicity. Encapsulating mRNA in extracellular vesicles may protect the mRNA and reduce reactogenicity, making mucosal mRNA vaccines possible. Plant-derived extracellular vesicles from edible fruits have been investigated as mRNA carriers. Studies in animals show that mRNA vehiculated in orange-derived extracellular vesicles can elicit both systemic and mucosal immune responses when administered by the oral, nasal, or intramuscular routes. Once lyophilized, these products show remarkable stability. The optimization of mRNA to improve translation efficiency, immunogenicity, reactogenicity, and stability can be obtained through adjustments of the 5'cap region, poly-A tail, codons selection, and the use of nucleoside analogues. Recent studies have also proposed self-amplifying RNA vaccines containing an RNA polymerase as well as circular mRNA constructs. Data from parenterally primed animals demonstrate the efficacy of nasal immunization with non-adjuvanted protein, and studies in humans indicate that the combination of a parenteral vaccine with the natural exposure of mucosae to the same antigen provides protection and reduces transmission. Hence, mucosal mRNA vaccination would be beneficial at least in organisms pre-treated with parenteral vaccines. This practice could have wide applications for the treatment of infectious diseases.

Viruses as archaeological tools for uncovering ancient molecular relationships

The entry of a virus into the host cell always implies the alteration of certain intracellular molecular relationships, some of which may involve the recovery of ancient cellular activities. In this sense, viruses are archaeological tools for identifying unexpressed activities in noninfected cells. Among these, activities that hinder virus propagation may represent cellular defense mechanisms, for example, activities that mutagenize the viral genome such as ADAR-1 or APOBEC activities. Instead, those that facilitate virus propagation can be interpreted as the result of viral adaptation to-or mimicking-cellular structures, enabling the virus to perform anthropomorphic activities, including hijacking, manipulating, and reorganizing cellular factors for their own benefit. The alternative we consider here is that some of these second set of cellular activities were already in the uninfected cell but silenced, under the negative control of the cell or lineage, and that they represent a necessary precondition for viral infection. For example, specifically loading an amino acid at the 3'-end of the mRNA of some plant viruses by aminoacyl-tRNA synthetases has proved essential for virus infection despite this reaction not occurring with cellular mRNAs. Other activities of this type are discussed here, together with the biological context in which they acquire a coherent meaning, that is, genetic latency and molecular conflict.

Epidemiological Dynamics of Foot-and-Mouth Disease in the Horn of Africa: The Role of Virus Diversity and Animal Movement

The Horn of Africa is a large area of arid and semi-arid land, holding about 10% of the global and 40% of the entire African livestock population. The region's livestock production system is mainly extensive and pastoralist. It faces countless problems, such as a shortage of pastures and watering points, poor access to veterinary services, and multiple endemic diseases like foot-and-mouth disease (FMD). Foot-and-mouth disease is one of the most economically important livestock diseases worldwide and is endemic in most developing countries. Within Africa, five of the seven serotypes of the FMD virus (FMDV) are described, but serotype C is not circulating anymore, a burden unseen anywhere in the world. The enormous genetic diversity of FMDV is favored by an error-prone RNA-dependent RNA polymerase, intra-typic and inter-typic recombination, as well as the quasi-species nature of the virus. This paper describes the epidemiological dynamics of foot-and-mouth disease in the Horn of Africa with regard to the serotypes and topotypes distribution of FMDV, the livestock production systems practiced, animal movement, the role of wildlife, and the epidemiological complexity of FMD. Within this review, outbreak investigation data and serological studies confirm the endemicity of the disease in the Horn of Africa. Multiple topotypes of FMDV are described in the literature as circulating in the region, with further evolution of virus diversity predicted. A large susceptible livestock population and the presence of wild ungulates are described as complicating the epidemiology of the disease. Further, the husbandry practices and legal and illegal trading of livestock and their products, coupled with poor biosecurity practices, are also reported to impact the spread of FMDV within and between countries in the region. The porosity of borders for pastoralist herders fuels the unregulated transboundary livestock trade. There are no systematic control strategies in the region except for sporadic vaccination with locally produced vaccines, while literature indicates that effective control measures should also consider virus diversity, livestock movements/biosecurity, transboundary trade, and the reduction of contact with wild, susceptible ungulates.

Hepatitis C virus fitness can influence the extent of infection-mediated epigenetic modifications in the host cells

Introduction

Cellular epigenetic modifications occur in the course of viral infections. We previously documented that hepatitis C virus (HCV) infection of human hepatoma Huh-7.5 cells results in a core protein-mediated decrease of Aurora kinase B (AURKB) activity and phosphorylation of Serine 10 in histone H3 (H3Ser10ph) levels, with an affectation of inflammatory pathways. The possible role of HCV fitness in infection-derived cellular epigenetic modifications is not known.

Methods

Here we approach this question using HCV populations that display a 2.3-fold increase in general fitness (infectious progeny production), and up to 45-fold increase of the exponential phase of intracellular viral growth rate, relative to the parental HCV population.

Results

We show that infection resulted in a HCV fitness-dependent, average decrease of the levels of H3Ser10ph, AURKB, and histone H4 tri-methylated at Lysine 20 (H4K20m3) in the infected cell population. Remarkably, the decrease of H4K20m3, which is a hallmark of cellular transformation, was significant upon infection with high fitness HCV but not upon infection with basal fitness virus.

Discussion

Here we propose two mechanisms ─which are not mutually exclusive─ to explain the effect of high viral fitness: an early advance in the number of infected cells, or larger number of replicating RNA molecules per cell. The implications of introducing HCV fitness as an influence in virus-host interactions, and for the course of liver disease, are warranted. Emphasis is made in the possibility that HCV-mediated hepatocellular carcinoma may be favoured by prolonged HCV infection of a human liver, a situation in which viral fitness is likely to increase.

Viral Fitness, Population Complexity, Host Interactions, and Resistance to Antiviral Agents

Fitness of viruses has become a standard parameter to quantify their adaptation to a biological environment. Fitness determinations for RNA viruses (and some highly variable DNA viruses) meet with several uncertainties. Of particular interest are those that arise from mutant spectrum complexity, absence of population equilibrium, and internal interactions among components of a mutant spectrum. Here, concepts, fitness measurements, limitations, and current views on experimental viral fitness landscapes are discussed. The effect of viral fitness on resistance to antiviral agents is covered in some detail since it constitutes a widespread problem in antiviral pharmacology, and a challenge for the design of effective antiviral treatments. Recent evidence with hepatitis C virus suggests the operation of mechanisms of antiviral resistance additional to the standard selection of drug-escape mutants. The possibility that high replicative fitness may be the driver of such alternative mechanisms is considered. New broad-spectrum antiviral designs that target viral fitness may curtail the impact of drug-escape mutants in treatment failures. We consider to what extent fitness-related concepts apply to coronaviruses and how they may affect strategies for COVID-19 prevention and treatment.

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.

SARS-CoV-2 Mutant Spectra at Different Depth Levels Reveal an Overwhelming Abundance of Low Frequency Mutations

Populations of RNA viruses are composed of complex and dynamic mixtures of variant genomes that are termed mutant spectra or mutant clouds. This applies also to SARS-CoV-2, and mutations that are detected at low frequency in an infected individual can be dominant (represented in the consensus sequence) in subsequent variants of interest or variants of concern. Here we briefly review the main conclusions of our work on mutant spectrum characterization of hepatitis C virus (HCV) and SARS-CoV-2 at the nucleotide and amino acid levels and address the following two new questions derived from previous results: (i) how is the SARS-CoV-2 mutant and deletion spectrum composition in diagnostic samples, when examined at progressively lower cut-off mutant frequency values in ultra-deep sequencing; (ii) how the frequency distribution of minority amino acid substitutions in SARS-CoV-2 compares with that of HCV sampled also from infected patients. The main conclusions are the following: (i) the number of different mutations found at low frequency in SARS-CoV-2 mutant spectra increases dramatically (50- to 100-fold) as the cut-off frequency for mutation detection is lowered from 0.5% to 0.1%, and (ii) that, contrary to HCV, SARS-CoV-2 mutant spectra exhibit a deficit of intermediate frequency amino acid substitutions. The possible origin and implications of mutant spectrum differences among RNA viruses are discussed.

Drug Sensitivity of Vaccine-Derived Rubella Viruses and Quasispecies Evolution in Granulomatous Lesions of Two Ataxia-Telangiectasia Patients Treated with Nitazoxanide

A strong association between rubella virus (RuV) and chronic granulomas, in individuals with inborn errors of immunity, has been recently established. Both the RA27/3 vaccine and wild-type RuV strains were highly sensitive to a broad-spectrum antiviral drug, nitazoxanide (NTZ), in vitro. However, NTZ treatment, used as a salvage therapy, resulted in little or no improvements of RuV-associated cutaneous granulomas in patients. Here, we report investigations of possible causes of treatment failures in two ataxia-telangiectasia patients. Although a reduction in RuV RNA in skin lesions was detected by real-time RT-PCR, live immunodeficiency-related vaccine-derived rubella viruses (iVDRV) were recovered from granulomas, before and after the treatments. Tizoxanide, an active NTZ metabolite, inhibited replications of all iVDRVs in cultured A549 cells, but the 50% and 90% inhibitory concentrations were 10–40 times higher than those for the RA27/3 strain. There were no substantial differences in iVDRV sensitivities, neither before nor after treatments. Analysis of quasispecies in the E1 gene, a suspected NTZ target, showed no effect of NTZ treatments on quasispecies’ complexity in lesions. Thus, failures of NTZ therapies were likely due to low sensitivities of iVDRVs to the drug, and not related to the emergence of resistance, following long-term NTZ treatments.

A Two-Level, Intramutant Spectrum Haplotype Profile of Hepatitis C Virus Revealed by Self-Organized Maps

RNA viruses replicate as complex mutant spectra termed viral quasispecies. The frequency of each individual genome in a mutant spectrum depends on its rate of generation and its relative fitness in the replicating population ensemble. The advent of deep sequencing methodologies allows for the first-time quantification of haplotype abundances within mutant spectra. There is no information on the haplotype profile of the resident genomes and how the landscape evolves when a virus replicates in a controlled cell culture environment. Here, we report the construction of intramutant spectrum haplotype landscapes of three amplicons of the NS5A-NS5B coding region of hepatitis C virus (HCV). Two-dimensional (2D) neural networks were constructed for 44 related HCV populations derived from a common clonal ancestor that was passaged up to 210 times in human hepatoma Huh-7.5 cells in the absence of external selective pressures. The haplotype profiles consisted of an extended dense basal platform, from which a lower number of protruding higher peaks emerged. As HCV increased its adaptation to the cells, the number of haplotype peaks within each mutant spectrum expanded, and their distribution shifted in the 2D network. The results show that extensive HCV replication in a monotonous cell culture environment does not limit HCV exploration of sequence space through haplotype peak movements. The landscapes reflect dynamic variation in the intramutant spectrum haplotype profile and may serve as a reference to interpret the modifications produced by external selective pressures or to compare with the landscapes of mutant spectra in complex in vivo environments. IMPORTANCE The study provides for the first time the haplotype profile and its variation in the course of virus adaptation to a cell culture environment in the absence of external selective constraints. The deep sequencing-based self-organized maps document a two-layer haplotype distribution with an ample basal platform and a lower number of protruding peaks. The results suggest an inferred intramutant spectrum fitness landscape structure that offers potential benefits for virus resilience to mutational inputs.

Historical Perspective on the Discovery of the Quasispecies Concept

Viral quasispecies are dynamic distributions of nonidentical but closely related mutant and recombinant viral genomes subjected to a continuous process of genetic variation, competition, and selection that may act as a unit of selection. The quasispecies concept owes its theoretical origins to a model for the origin of life as a collection of mutant RNA replicators. Independently, experimental evidence for the quasispecies concept was obtained from sampling of bacteriophage clones, which revealed that the viral populations consisted of many mutant genomes whose frequency varied with time of replication. Similar findings were made in animal and plant RNA viruses. Quasispecies became a theoretical framework to understand viral population dynamics and adaptability. The evidence came at a time when mutations were considered rare events in genetics, a perception that was to change dramatically in subsequent decades. Indeed, viral quasispecies was the conceptual forefront of a remarkable degree of biological diversity, now evident for cell populations and organisms, not only for viruses. Quasispecies dynamics unveiled complexities in the behavior of viral populations,with consequences for disease mechanisms and control strategies. This review addresses the origin of the quasispecies concept, its major implications on both viral evolution and antiviral strategies, and current and future prospects.

Mutation Rates, Mutation Frequencies, and Proofreading-Repair Activities in RNA Virus Genetics

The error rate displayed during template copying to produce viral RNA progeny is a biologically relevant parameter of the replication complexes of viruses. It has consequences for virus-host interactions, and it represents the first step in the diversification of viruses in nature. Measurements during infections and with purified viral polymerases indicate that mutation rates for RNA viruses are in the range of 10-3 to 10-6 copying errors per nucleotide incorporated into the nascent RNA product. Although viruses are thought to exploit high error rates for adaptation to changing environments, some of them possess misincorporation correcting activities. One of them is a proofreading-repair 3' to 5' exonuclease present in coronaviruses that may decrease the error rate during replication. Here we review experimental evidence and models of information maintenance that explain why elevated mutation rates have been preserved during the evolution of RNA (and some DNA) viruses. The models also offer an interpretation of why error correction mechanisms have evolved to maintain the stability of genetic information carried out by large viral RNA genomes such as the coronaviruses.

Population Disequilibrium as Promoter of Adaptive Explorations in Hepatitis C Virus

Replication of RNA viruses is characterized by exploration of sequence space which facilitates their adaptation to changing environments. It is generally accepted that such exploration takes place mainly in response to positive selection, and that further diversification is boosted by modifications of virus population size, particularly bottleneck events. Our recent results with hepatitis C virus (HCV) have shown that the expansion in sequence space of a viral clone continues despite prolonged replication in a stable cell culture environment. Diagnosis of the expansion was based on the quantification of diversity indices, the occurrence of intra-population mutational waves (variations in mutant frequencies), and greater individual residue variations in mutant spectra than those anticipated from sequence alignments in data banks. In the present report, we review our previous results, and show additionally that mutational waves in amplicons from the NS5A-NS5B-coding region are equally prominent during HCV passage in the absence or presence of the mutagenic nucleotide analogues favipiravir or ribavirin. In addition, by extending our previous analysis to amplicons of the NS3- and NS5A-coding region, we provide further evidence of the incongruence between amino acid conservation scores in mutant spectra from infected patients and in the Los Alamos National Laboratory HCV data banks. We hypothesize that these observations have as a common origin a permanent state of HCV population disequilibrium even upon extensive viral replication in the absence of external selective constraints or changes in population size. Such a persistent disequilibrium—revealed by the changing composition of the mutant spectrum—may facilitate finding alternative mutational pathways for HCV antiviral resistance. The possible significance of our model for other genetically variable viruses is discussed.

Intrapatient Evolutionary Dynamics in an Individual Infected with HIV-1 CRF01_AE Who Experienced Periods of Treatment Failure

Although previous studies have analyzed cross-level CRF01_AE viral genomic data in populations, less is known about intrapatient viral evolutionary dynamics during antiretroviral therapy (ART) failure. We longitudinally sampled plasma and peripheral blood mononuclear cells (PBMC) at different time points from one human immunodeficiency virus type 1 infected patient. The evolution of viral quasispecies was inferred from viral phylogenies. Before treatment, no drug-resistant mutations were found in this patient's plasma, and all viruses had C-C chemokine receptor type 5 (CCR5) tropism. Two months after treatment, the majority of the virus population in plasma and PBMC were drug resistant and X4-tropic. By 5 months after treatment, the viral load increased significantly, and viruses reversed tropism from X4 to R5 in plasma and PBMC. During treatment failure, the effective population of the pol DNA reservoir in PBMC remained stable, whereas the env DNA reservoir increased. The effective population of the R5 tropism virus increased more rapidly than that of the X4 tropism virus. The ratio of non-synonymous to synonymous substitutions in the env gene of R5 tropism virus (0.43) was lower than X4 tropism (0.52). However, four env positive selection sites were identified in R5 tropism viruses (HXB2: 364, 398, 399, and 400) but none were identified in X4 tropism viruses. Our data demonstrated the different intrapatient evolutionary dynamics patterns of env and pol genes in an individual who experienced periods of ART failure. Our findings also suggest the importance of the R5 tropism virus in the DNA reservoir during ART failure.

Amino Acid Substitutions Associated with Treatment Failure for Hepatitis C Virus Infection

Despite the high virological response rates achieved with current directly acting antiviral agents (DAAs) against hepatitis C virus (HCV), around 2% to 5% of treated patients do not achieve a sustained viral response. The identification of amino acid substitutions associated with treatment failure requires analytical designs, such as subtype-specific ultradeep sequencing (UDS) methods, for HCV characterization and patient management. Using this procedure, we have identified six highly represented amino acid substitutions (HRSs) in NS5A and NS5B of HCV, which are not bona fide resistance-associated substitutions (RAS), from 220 patients who failed therapy. They were present frequently in basal and posttreatment virus of patients who failed different DAA-based therapies. Contrary to several RAS, HRSs belong to the acceptable subset of substitutions according to the PAM250 replacement matrix. Their mutant frequency, measured by the number of deep sequencing reads within the HCV quasispecies that encode the relevant substitutions, ranged between 90% and 100% in most cases. They also have limited predicted disruptive effects on the three-dimensional structures of the proteins harboring them. Possible mechanisms of HRS origin and dominance, as well as their potential predictive value for treatment response, are discussed.

Dissimilar Conservation Pattern in Hepatitis C Virus Mutant Spectra, Consensus Sequences, and Data Banks

The influence of quasispecies dynamics on long-term virus diversification in nature is a largely unexplored question. Specifically, whether intra-host nucleotide and amino acid variation in quasispecies fit the variation observed in consensus sequences or data bank alignments is unknown. Genome conservation and dynamics simulations are used for the computational design of universal vaccines, therapeutic antibodies and pan-genomic antiviral agents. The expectation is that selection of escape mutants will be limited when mutations at conserved residues are required. This strategy assumes long-term (epidemiologically relevant) conservation but, critically, does not consider short-term (quasispecies-dictated) residue conservation. We calculated mutant frequencies of individual loci from mutant spectra of hepatitis C virus (HCV) populations passaged in cell culture and from infected patients. Nucleotide or amino acid conservation in consensus sequences of the same populations, or in the Los Alamos HCV data bank did not match residue conservation in mutant spectra. The results relativize the concept of sequence conservation in viral genetics and suggest that residue invariance in data banks is an insufficient basis for the design of universal viral ligands for clinical purposes. Our calculations suggest relaxed mutational restrictions during quasispecies dynamics, which may contribute to higher calculated short-term than long-term viral evolutionary rates.

LUCS: a high-resolution nucleic acid sequencing tool for accurate long-read analysis of individual DNA molecules

Nucleic acid sequence analyses are fundamental to all aspects of biological research, spanning aging, mitochondrial DNA (mtDNA) and cancer, as well as microbial and viral evolution. Over the past several years, significant improvements in DNA sequencing, including consensus sequence analysis, have proven invaluable for high-throughput studies. However, all current DNA sequencing platforms have limited utility for studies of complex mixtures or of individual long molecules, the latter of which is crucial to understanding evolution and consequences of single nucleotide variants and their combinations. Here we report a new technology termed LUCS (Long-molecule UMI-driven Consensus Sequencing), in which reads from third-generation sequencing are aggregated by unique molecular identifiers (UMIs) specific for each individual DNA molecule. This enables in-silico reconstruction of highly accurate consensus reads of each DNA molecule independent of other molecules in the sample. Additionally, use of two UMIs enables detection of artificial recombinants (chimeras). As proof of concept, we show that application of LUCS to assessment of mitochondrial genomes in complex mixtures from single cells was associated with an error rate of 1X10^-4 errors/nucleotide. Thus, LUCS represents a major step forward in DNA sequencing that offers high-throughput capacity and high-accuracy reads in studies of long DNA templates and nucleotide variants in heterogenous samples.