Rapid viral quasispecies evolution: implications for vaccine and drug strategies

Host Cell Targets for Unconventional Antivirals against RNA Viruses

The recent COVID-19 crisis has highlighted the importance of RNA-based viruses. The most prominent members of this group are SARS-CoV-2 (coronavirus), HIV (human immunodeficiency virus), EBOV (Ebola virus), DENV (dengue virus), HCV (hepatitis C virus), ZIKV (Zika virus), CHIKV (chikungunya virus), and influenza A virus. With the exception of retroviruses which produce reverse transcriptase, the majority of RNA viruses encode RNA-dependent RNA polymerases which do not include molecular proofreading tools, underlying the high mutation capacity of these viruses as they multiply in the host cells. Together with their ability to manipulate the immune system of the host in different ways, their high mutation frequency poses a challenge to develop effective and durable vaccination and/or treatments. Consequently, the use of antiviral targeting agents, while an important part of the therapeutic strategy against infection, may lead to the selection of drug-resistant variants. The crucial role of the host cell replicative and processing machinery is essential for the replicative cycle of the viruses and has driven attention to the potential use of drugs directed to the host machinery as therapeutic alternatives to treat viral infections. In this review, we discuss small molecules with antiviral effects that target cellular factors in different steps of the infectious cycle of many RNA viruses. We emphasize the repurposing of FDA-approved drugs with broad-spectrum antiviral activity. Finally, we postulate that the ferruginol analog (18-(phthalimide-2-yl) ferruginol) is a potential host-targeted antiviral.

Exploring the application of immunotherapy against HIV infection in the setting of malignancy: A detailed review article

According to the Joint United Nations Programme on HIV/AIDS (UNAIDS), as of 2019, approximately 42.2 million people have died from acquired immunodeficiency syndrome (AIDS)-related illnesses since the start of the epidemic. Antiretroviral therapy (ART) has significantly reduced mortality, morbidity, and incidence of the human immunodeficiency virus (HIV)/AIDS-defining cancers, taming once-dreaded disease into a benign chronic infection. Although the treatment has prolonged the patients' survival, general HIV prevalence has increased and this increase has dovetailed with an increasing incidence of Non-AIDS-defining cancers (NADCs) among people living with HIV (PLWH). This is happening when new promising approaches in both oncology and HIV infection are being developed. This review focuses on recent progress witnessed in immunotherapy approaches against HIV-related, Non-AIDS-defining cancers (NADCs), and HIV infection.

VERSO: A comprehensive framework for the inference of robust phylogenies and the quantification of intra-host genomic diversity of viral samples

We introduce VERSO, a two-step framework for the characterization of viral evolution from sequencing data of viral genomes, which is an improvement on phylogenomic approaches for consensus sequences. VERSO exploits an efficient algorithmic strategy to return robust phylogenies from clonal variant profiles, also in conditions of sampling limitations. It then leverages variant frequency patterns to characterize the intra-host genomic diversity of samples, revealing undetected infection chains and pinpointing variants likely involved in homoplasies. On simulations, VERSO outperforms state-of-the-art tools for phylogenetic inference. Notably, the application to 6,726 amplicon and RNA sequencing samples refines the estimation of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) evolution, while co-occurrence patterns of minor variants unveil undetected infection paths, which are validated with contact tracing data. Finally, the analysis of SARS-CoV-2 mutational landscape uncovers a temporal increase of overall genomic diversity and highlights variants transiting from minor to clonal state and homoplastic variants, some of which fall on the spike gene. Available at: https://github.com/BIMIB-DISCo/VERSO.

Mutational signatures and heterogeneous host response revealed via large-scale characterization of SARS-CoV-2 genomic diversity

To dissect the mechanisms underlying the inflation of variants in the Severe Acute Respiratory Syndrome CoronaVirus 2 (SARS-CoV-2) genome, we present a large-scale analysis of intra-host genomic diversity, which reveals that most samples exhibit heterogeneous genomic architectures, due to the interplay between host-related mutational processes and transmission dynamics. The decomposition of minor variants profiles unveils three non-overlapping mutational signatures related to nucleotide substitutions and likely ruled by APOlipoprotein B Editing Complex (APOBEC), Reactive Oxygen Species (ROS), and Adenosine Deaminase Acting on RNA (ADAR), highlighting heterogeneous host responses to SARS-CoV-2 infections. A corrected-for-signatures dN/dS analysis demonstrates that such mutational processes are affected by purifying selection, with important exceptions. In fact, several mutations appear to transit toward clonality, defining new clonal genotypes that increase the overall genomic diversity. Furthermore, the phylogenomic analysis shows the presence of homoplasies and supports the hypothesis of transmission of minor variants. This study paves the way for the integrated analysis of intra-host genomic diversity and clinical outcomes of SARS-CoV-2 infections.

Infection of calves with in-vivo passaged bovine parainfluenza-3 virus, alone or in combination with bovine respiratory syncytial virus and bovine coronavirus

Bovine respiratory disease complex is etiologically complex and usually involves co-infection by several agents, including bovine parainfluenza virus-3 (BPIV-3), bovine respiratory syncytial virus (BRSV), and bovine coronavirus (BCoV). Traditionally, vaccines have been tested in seronegative calves infected with a single in vitro-passaged agent, often with little disease, resulting in unvaccinated subjects. To overcome the potential problem of attenuation coincident with in vitro culture of the viruses, cocktails of field isolates of BPIV-3s and BCoVs were passaged in the lungs of neonatal colostrum-deprived calves. Lung lavage fluids were used as inocula, alone and in combination with in-vivo passaged BRSV, and aerosolized into a trailer containing conventionally reared 9-week-old weaned Holstein calves with decayed, but still measurable, maternal antibodies. Calves developed acute respiratory disease of variable severity. Upon necropsy, there were characteristic gross and histologic lesions in the respiratory tract, associated immunohistochemically with BPIV-3, BRSV, and BCoV. In-vivo passage of viruses is an alternative to in vitro culture to produce inocula to better study the pathogenesis of infection and more rigorously and relevantly assess vaccine efficacy.

Comparative mutational analyses of influenza A viruses

The error-prone RNA-dependent RNA polymerase (RdRP) and external selective pressures are the driving forces for RNA viral diversity. When confounded by selective pressures, it is difficult to assess if influenza A viruses (IAV) that have a wide host range possess comparable or distinct spontaneous mutational frequency in their RdRPs. We used in-depth bioinformatics analyses to assess the spontaneous mutational frequencies of two RdRPs derived from human seasonal (A/Wuhan/359/95; Wuhan) and H5N1 (A/Vietnam/1203/04; VN1203) viruses using the mini-genome system with a common firefly luciferase reporter serving as the template. High-fidelity reverse transcriptase was applied to generate high-quality mutational spectra which allowed us to assess and compare the mutational frequencies and mutable motifs along a target sequence of the two RdRPs of two different subtypes. We observed correlated mutational spectra (τ correlation P < 0.0001), comparable mutational frequencies (H3N2:5.8 ± 0.9; H5N1:6.0 ± 0.5), and discovered a highly mutable motif "(A)AAG" for both Wuhan and VN1203 RdRPs. Results were then confirmed with two recombinant A/Puerto Rico/8/34 (PR8) viruses that possess RdRP derived from Wuhan or VN1203 (RG-PR8×Wuhan(PB2, PB1, PA, NP) and RG-PR8×VN1203(PB2, PB1, PA, NP)). Applying novel bioinformatics analysis on influenza mutational spectra, we provide a platform for a comprehensive analysis of the spontaneous mutation spectra for an RNA virus.

Viral quasispecies evolution

Evolution of RNA viruses occurs through disequilibria of collections of closely related mutant spectra or mutant clouds termed viral quasispecies. Here we review the origin of the quasispecies concept and some biological implications of quasispecies dynamics. Two main aspects are addressed: (i) mutant clouds as reservoirs of phenotypic variants for virus adaptability and (ii) the internal interactions that are established within mutant spectra that render a virus ensemble the unit of selection. The understanding of viruses as quasispecies has led to new antiviral designs, such as lethal mutagenesis, whose aim is to drive viruses toward low fitness values with limited chances of fitness recovery. The impact of quasispecies for three salient human pathogens, human immunodeficiency virus and the hepatitis B and C viruses, is reviewed, with emphasis on antiviral treatment strategies. Finally, extensions of quasispecies to nonviral systems are briefly mentioned to emphasize the broad applicability of quasispecies theory.

Applications of next-generation sequencing technologies to diagnostic virology

Novel DNA sequencing techniques, referred to as “next-generation” sequencing (NGS), provide high speed and throughput that can produce an enormous volume of sequences with many possible applications in research and diagnostic settings. In this article, we provide an overview of the many applications of NGS in diagnostic virology. NGS techniques have been used for high-throughput whole viral genome sequencing, such as sequencing of new influenza viruses, for detection of viral genome variability and evolution within the host, such as investigation of human immunodeficiency virus and human hepatitis C virus quasispecies, and monitoring of low-abundance antiviral drug-resistance mutations. NGS techniques have been applied to metagenomics-based strategies for the detection of unexpected disease-associated viruses and for the discovery of novel human viruses, including cancer-related viruses. Finally, the human virome in healthy and disease conditions has been described by NGS-based metagenomics.

Models of RNA virus evolution and their roles in vaccine design

Viruses are fast evolving pathogens that continuously adapt to the highly variable environments they live and reproduce in. Strategies devoted to inhibit virus replication and to control their spread among hosts need to cope with these extremely heterogeneous populations and with their potential to avoid medical interventions. Computational techniques such as phylogenetic methods have broadened our picture of viral evolution both in time and space, and mathematical modeling has contributed substantially to our progress in unraveling the dynamics of virus replication, fitness, and virulence. Integration of multiple computational and mathematical approaches with experimental data can help to predict the behavior of viral pathogens and to anticipate their escape dynamics. This piece of information plays a critical role in some aspects of vaccine development, such as viral strain selection for vaccinations or rational attenuation of viruses. Here we review several aspects of viral evolution that can be addressed quantitatively, and we discuss computational methods that have the potential to improve vaccine design.

Mutant spectra in virus behavior

RNA viruses replicate as complex mutant spectra, also termed ‘mutant clouds’, known as viral quasispecies. While this is a widely observed viral population structure, it is less known that a number of biologically relevant features of this important group of viral pathogens depend on (or are strongly influenced by) the complexity and composition of mutant spectra. Among them, fitness increase or decrease depending on intrapopulation complementation or interference, selection triggered by memory genomes, pathogenic potential of viruses, disease evolution and the response to antiviral treatments. Quasispecies represent the recognition of complex behavior in viruses, and it is an oversimplification to equate such a population structure with the classic polymorphism of population biology. Darwinian principles acting on genome collectivities that replicate with high error rates provide a unique population structure prone to flexible and largely unpredictable behavior.

Massively parallel sequencing for monitoring genetic consistency and quality control of live viral vaccines

Intrinsic genetic instability of RNA viruses may lead to the accumulation of revertants during manufacture of live viral vaccines, requiring rigorous quality control to ensure vaccine safety. Each lot of oral poliovirus vaccine (OPV) is tested for neurovirulence in animals and also for the presence of neurovirulent revertants. Mutant analysis by PCR and restriction enzyme cleavage (MAPREC) is used to measure the frequency of neurovirulent mutations at the 5' untranslated region (UTR) of the viral genome that correlate with the level of neurovirulence determined by the monkey neurovirulence test. However, MAPREC can only monitor mutations at a few genomic loci and miss mutations at other sites that could adversely affect vaccine quality. Here we propose to use massively parallel sequencing (MPS) for sensitive detection and quantification of all mutations in the entire genome of attenuated viruses. Analysis of vaccine samples and reference preparations demonstrated a perfect agreement with MAPREC results. Quantitative MPS analysis of validated reference preparations tested by MAPREC produced identical results, suggesting that the method could take advantage of the existing reference materials and be used as a replacement for the MAPREC procedure in lot release of OPV. Patterns of mutations present at a low level in vaccine preparations were characteristic of seed viruses used for their manufacture and could be used for identification of individual batches. This approach may represent the ultimate tool for monitoring genetic consistency of live viral vaccines.

Information catastrophe in RNA viruses through replication thresholds

RNA viruses are known to replicate at very high mutation rates. These rates are actually known to be close to their so-called error threshold. This threshold is in fact a critical point beyond which genetic information is lost through a so-called error catastrophe. However, the transition from a stable quasispecies to genetic drift and loss of information can also occur by crossing replication thresholds, below some replication rates, the viral population is suddenly unable to survive. Available data from hepatitis C virus population analysis [Mas, A., Ulloa, E., Bruguera, M., Furcić, I., Garriga, D., Fábregas, S., Andreu, D., Saiz, J.C., Díez, J., 2004. Hepatitis C virus population analysis of a single-source nosocomial outbreak reveals an inverse correlation between viral load and quasispecies complexity. J. Gen. Virol. 85, 3619-3626] can be interpreted through this theoretical view, providing evidence for such a replication threshold. Here a simple model is used in order to provide evidence for such a phenomenon, consistent with available data.

Strategies and mechanisms for host and pathogen survival in acute and persistent viral infections

Persistent viral infections causing serious diseases derive, primarily, from altered function of the immune system. Knowledge of the very complex composition and function of the innate and adaptive branches of the immune system is essential to understanding persistent infection. The best solution to the problem of persistent infection is by prevention using prophylactic vaccines. Hit and run viruses evade immune destruction by infecting new hosts and rarely persist. Hit and stay viruses evade immune control by sequestration, blockade of antigen presentation, cytokine escape, evasion of natural killer cell activities, escape from apoptosis, and antigenic change. Twelve prophylactic vaccines against hit and run agents exist, and there are only three vaccines against hit and stay viruses, all of which are of DNA composition. Several new vaccines against hit and stay viruses are feasible, but protective vaccines against RNA HIV and hepatitis C agents are highly unlikely, short of a major breakthrough. Therapeutic vaccines are very improbable without a magnitude of favorable new discoveries. In the meantime, antiviral chemotherapy, chemotherapy/prophylactic vaccination, and short interfering RNA silencing are worthy of intense investigation.

Evolution of foot-and-mouth disease virus

Foot-and-mouth disease virus evolution is strongly influenced by high mutation rates and a quasispecies dynamics. Mutant swarms are subjected to positive selection, negative selection and random drift of genomes. Adaptation is the result of selective amplification of subpopulations of genomes. The extent of adaptation to a given environment is quantified by a relative fitness value. Fitness values depend on the virus and its physical and biological environment. Generally, infections involving large population passages result in fitness gain and population bottlenecks lead to fitness loss. Very different types of mutations tend to accumulate in the foot-and-mouth disease virus (FMDV) genome depending on the virus population size during replication. Quasispecies dynamics predict higher probability of success of antiviral strategies based on multivalent vaccines and combination therapy, and this has been supported by clinical and veterinary practice. Quasispecies suggest also new antiviral strategies based on virus entry into error catastrophe, and such procedures are under investigation. Studies with FMDV have contributed to the understanding of quasispecies dynamics and some of its biological implications.

Virus population dynamics, fitness variations and the control of viral disease: an update

Viral quasispecies dynamics and variations of viral fitness are reviewed in connection with viral disease control. Emphasis is put on resistance of human immunodeficiency virus and some human DNA viruses to antiviral inhibitors. Future trends in multiple target antiviral therapy and new approaches based on virus entry into error catastrophe (extinction mutagenesis) are discussed.

Quasispecies development by high frequency RNA recombination during MHV persistence

Recent studies suggest that infectious viruses and particularly persisting viral RNAs often exist as diverse populations or "quasispecies". We have developed an approach to characterize populations of the murine coronavirus mouse hepatitis virus (MHV) generated during persistent infection which has allowed us to begin to address the role of the viral quasispecies in MHV pathogenesis. We analyzed the population of persisting viral RNAs using reverse-transcription polymerase chain reaction amplification (RT-PCR) of the S1 "hypervariable" region of the spike gene followed by differential colony hybridization to identify spike deletion variants (SDVs) from acute and persistently infected mice. Sequence analysis revealed that mice with the most severe chronic paralysis harbored the most complex quasispecies. Mapping of the SDVs to the predicted RNA secondary structure of the spike RNA revealed that an isolated stem loop structure is frequently deleted. Overall, these results are consistent with high frequency recombination at sites of RNA secondary structure contributing to expansion of the viral quasispecies and persisting viral pathogenesis.

Quasispecies and the implications for virus persistence and escape

BACKGROUND

In the 1970s Manfred Eigen and colleagues proposed a new model of molecular evolution to explain adaptability and rapid evolution of simple replicons, as those that probably populated the earth at the onset of life. This model of evolution placed emphasis on mutant generation, to the point of invalidating the concept of wild-type genomes as a defined sequence of nucleotides. In striking similarity with the proposals for such early replicons, present-day RNA viruses consist of complex distributions of nonidentical but closely related genomes termed quasispecies.

OBJECTIVES

To discuss indeterminations inherent to a quasispecies structure and to the analytical procedures to define it, biological implications of quasispecies, and the need to take into account this type of population structure, in order to design effective strategies to prevent and control diseases caused by highly variable viruses.

RESULTS

Quasispecies have many biological implications, extending from viral pathogenesis to the emergence of new pathogens, rapid antigenic variation, and alterations in cell tropism, virulence, host range and viral gene expression.

CONCLUSIONS

Diseases caused by highly variable RNA viruses prove very difficult to control and vaccine development against such viruses are largely unsuccessful. It is important to understand quasispecies composition and dynamics, as quasispecies are an important step in the natural history of RNA viruses.

Quasispecies structure and persistence of RNA viruses

Viral quasispecies are closely related (but nonidentical) mutant and recombinant viral genomes subjected to continuous genetic variation, competition, and selection. Quasispecies structure and dynamics of replicating RNA enable virus populations to persist in their hosts and cause disease. We review mechanisms of viral persistence in cells, organisms, and populations of organisms and suggest that the critical interplay between host and viral influences (including in some cases the quasispecies organization) is the main driving force for long-term survival of viruses in nature.

RNA virus mutations and fitness for survival

RNA viruses exploit all known mechanisms of genetic variation to ensure their survival. Distinctive features of RNA virus replication include high mutation rates, high yields, and short replication times. As a consequence, RNA viruses replicate as complex and dynamic mutant swarms, called viral quasispecies. Mutation rates at defined genomic sites are affected by the nucleotide sequence context on the template molecule as well as by environmental factors. In vitro hypermutation reactions offer a means to explore the functional sequence space of nucleic acids and proteins. The evolution of a viral quasispecies is extremely dependent on the population size of the virus that is involved in the infections. Repeated bottleneck events lead to average fitness losses, with viruses that harbor unusual, deleterious mutations. In contrast, large population passages result in rapid fitness gains, much larger than those so far scored for cellular organisms. Fitness gains in one environment often lead to fitness losses in an alternative environment. An important challenge in RNA virus evolution research is the assignment of phenotypic traits to specific mutations. Different constellations of mutations may be associated with a similar biological behavior. In addition, recent evidence suggests the existence of critical thresholds for the expression of phenotypic traits. Epidemiological as well as functional and structural studies suggest that RNA viruses can tolerate restricted types and numbers of mutations during any specific time point during their evolution. Viruses occupy only a tiny portion of their potential sequence space. Such limited tolerance to mutations may open new avenues for combating viral infections.

Generation of coronavirus spike deletion variants by high-frequency recombination at regions of predicted RNA secondary structure

Coronavirus RNA evolves in the central nervous systems (CNS) of mice during persistent infection. This evolution can be monitored by detection of a viral quasispecies of spike deletion variants (SDVs) (C. L. Rowe, S. C. Baker, M. J. Nathan, and J. O. Fleming, J. Virol. 71:2959-2969, 1997). We and others have found that the deletions cluster in the region from 1,200 to 1,800 nucleotides from the 5' end of the spike gene sequence, termed the "hypervariable" region. To address how SDVs might arise, we generated the predicted folding structures of the positive- and negative-strand senses of the entire 4,139-nt spike RNA sequence. We found that a prominent, isolated stem-loop structure is coincident with the hypervariable region in each structure. To determine if this predicted stem-loop is a "hot spot" for RNA recombination, we assessed whether this region of the spike is more frequently deleted than three other selected regions of the spike sequence in a population of viral sequences isolated from the CNS of acutely and persistently infected mice. Using differential colony hybridization of cloned spike reverse transcription-PCR products, we detected SDVs in which the hot spot was deleted but did not detect SDVs in which other regions of the spike sequence were exclusively deleted. Furthermore, sequence analysis and mapping of the crossover sites of 25 distinct patterns of SDVs showed that the majority of crossover sites clustered to two regions at the base of the isolated stem-loop, which we designated as high-frequency recombination sites 1 and 2. Interestingly, the majority of the left and right crossover sites of the SDVs were directly across from or proximal to one another, suggesting that these SDVs are likely generated by intramolecular recombination. Overall, our results are consistent with there being an important role for the spike RNA secondary structure as a contributing factor in the generation of SDVs during persistent infection.

RNA virus fitness

RNA viruses constitute the most abundant group of pathogens of man, animals and plants. They share high mutation rates which are in the range 10(-3) to 10(-5) misincorporations per nucleotide site and round of copying. This is due to the absence or low efficiency of proofreading-repair or postreplicative repair activities associated with replicating RNA. Populations of RNA viruses are extremely heterogeneous and form dynamic mutant swarms termed viral quasispecies. This genetic organisation implies that any individual mutant has only a fleeting existence; that is, RNA viral genomes are statistically defined but individually indeterminate. RNA viruses are able to accommodate their average nucleotide sequences to changes in environment. A parameter used to quantitate adaptation is fitness, or the relative ability of a virus to produce infectious progeny. Repeated transfers of one or a few particles (bottleneck events) generally lead to fitness losses. In contrast, large population passages allow competitive optimisation of mutant genomes and fitness gains. Of relevance to medical practice is the ability of viral quasispecies to overcome selective pressures imposed by vaccines and antiviral agents. Particularly dramatic have been the systematic isolations of HIV-1 mutants resistant to antiretroviral inhibitors in treated individuals. In addition to the ability of HIV-1 quasispecies to generate many mutant genomes in short times, calculations of mutation frequencies in the pol gene of HIV-1 populations have documented that mutations related to resistance to antiretroviral inhibitors preexist in the mutant swarms of HIV-1 quasispecies. It is not possible at present to anticipate whether a suitable drug cocktail may be capable of sustained inhibition of HIV-1 replication without selection of mutants resistant to the combination of antiviral agents. Copyright 1997 John Wiley & Sons, Ltd.

Rapid evolution of viral RNA genomes

Mutation rates during RNA virus replication are several orders of magnitude larger than those operating during replication of cellular DNA. This results in the continuous generation of mutant genomes and in their rating in competition with other variants present and arising in the population. The dynamic mutant distributions that constitute RNA virus populations are termed quasispecies. This concept has facilitated links between population genetics and virology and has a number of important implications for viral pathogenesis and the control of viral disease. One of them is that the mutant spectra in RNA viruses constitute large reservoirs of genetic and phenotypic variants with potentially altered biological properties. Individual mutants kept in a low proportion under a set of environmental conditions may become dominant following an environmental change. Relevant to this review are possible links between the alteration of quasispecies distributions and nutritional deficiencies and oxidative stress in cells. In addition to being a possible mechanism of viral pathogenesis, oxidative stress, and other environmental modifications resulting from nutritional imbalances, may promote population disequilibrium in replicating viruses. In particular, the increased mutagenesis mediated by oxidative DNA damage could also affect replicating RNA and integrated provirus, extending the mutant repertoire of viruses. Also, the impairment of humoral and cellular immune functions may delay or prevent viral clearance, leading to an expanded representation of viral mutants in the infected organism. Thus, nutritional deficiencies are a potential source of viral mutants with altered biological properties.