Coxsackievirus B3 mutator strains are attenuated in vivo

Molecular Evolutionary Dynamics of Coxsackievirus A6 Causing Hand, Foot, and Mouth Disease From 2021 to 2023 in China: Genomic Epidemiology Study

Background

Hand, foot, and mouth disease (HFMD) is a global public health concern, notably within the Asia-Pacific region. Recently, the primary pathogen causing HFMD outbreaks across numerous countries, including China, is coxsackievirus (CV) A6, one of the most prevalent enteroviruses in the world. It is a new variant that has undergone genetic recombination and evolution, which might not only induce modifications in the clinical manifestations of HFMD but also heighten its pathogenicity because of nucleotide mutation accumulation.

Objective

The study assessed the epidemiological characteristics of HFMD in China and characterized the molecular epidemiology of the major pathogen (CV-A6) causing HFMD. We attempted to establish the association between disease progression and viral genetic evolution through a molecular epidemiological study.

Methods

Surveillance data from the Chinese Center for Disease Control and Prevention from 2021 to 2023 were used to analyze the epidemiological seasons and peaks of HFMD in Henan, China, and capture the results of HFMD pathogen typing. We analyzed the evolutionary characteristics of all full-length CV-A6 sequences in the NCBI database and the isolated sequences in Henan. To characterize the molecular evolution of CV-A6, time-scaled tree and historical population dynamics regarding CV-A6 sequences were estimated. Additionally, we analyzed the isolated strains for mutated or missing amino acid sites compared to the prototype CV-A6 strain.

Results

The 2021-2023 epidemic seasons for HFMD in Henan usually lasted from June to August, with peaks around June and July. The monthly case reporting rate during the peak period ranged from 20.7% (4854/23,440) to 35% (12,135/34,706) of the total annual number of cases. Analysis of the pathogen composition of 2850 laboratory-confirmed cases identified 8 enterovirus serotypes, among which CV-A6 accounted for the highest proportion (652/2850, 22.88%). CV-A6 emerged as the major pathogen for HFMD in 2022 (203/732, 27.73%) and 2023 (262/708, 37.01%). We analyzed all CV-A6 full-length sequences in the NCBI database and the evolutionary features of viruses isolated in Henan. In China, the D3 subtype gradually appeared from 2011, and by 2019, all CV-A6 virus strains belonged to the D3 subtype. The VP1 sequences analyzed in Henan showed that its subtypes were consistent with the national subtypes. Furthermore, we analyzed the molecular evolutionary features of CV-A6 using Bayesian phylogeny and found that the most recent common ancestor of CV-A6 D3 dates back to 2006 in China, earlier than the 2011 HFMD outbreak. Moreover, the strains isolated in 2023 had mutations at several amino acid sites compared to the original strain.

Conclusions

The CV-A6 virus may have been introduced and circulating covertly within China prior to the large-scale HFMD outbreak. Our laboratory testing data confirmed the fluctuation and periodic patterns of CV-A6 prevalence. Our study provides valuable insights into understanding the evolutionary dynamics of CV-A6.

Fidelity Characterization of Highly Pathogenic Porcine Reproductive and Respiratory Syndrome Virus and NADC30-like Strain

The porcine reproductive and respiratory syndrome virus (PRRSV) has significantly impacted the global pork industry for over three decades. Its high mutation rates and frequent recombination greatly intensifies its epidemic and threat. To explore the fidelity characterization of Chinese highly pathogenic PRRSV JXwn06 and the NADC30-like strain CHsx1401, self-recombination and mutation in PAMs, MARC-145 cells, and pigs were assessed. In vitro, CHsx1401 displayed a higher frequency of recombination junctions and a greater diversity of junction types than JXwn06. In vivo, CHsx1401 exhibited fewer junction types yet maintained a higher junction frequency. Notably, JXwn06 showed more accumulation of mutations. To pinpoint the genomic regions influencing their fidelity, chimeric viruses were constructed, with the exchanged nsp9-10 regions between JXwn06 and CHsx1401. The SJn9n10 strain, which incorporates JXwn06's nsp9-10 into the CHsx1401 genome, demonstrated reduced sensitivity to nucleotide analogs compared to CHsx1401. Conversely, compared with JXwn06, the JSn9n10 strain showed increased sensitivity to these inhibitors. The swapped nsp9-10 also influences the junction frequency and accumulated mutations as their donor strains. The results indicate a propensity for different types of genetic variations between these two strains and further highlight the nsp9-10 region as a critical determinant of their fidelity.

Distinct chikungunya virus polymerase palm subdomains contribute to virus replication and virion assembly

Alphaviruses encode an error-prone RNA-dependent RNA polymerase (RdRp), nsP4, required for genome synthesis, yet how the RdRp functions in the complete alphavirus life cycle is not well-defined. Previous work using chikungunya virus (CHIKV) has established the importance of the nsP4 residue cysteine 483 in maintaining viral genetic fidelity. Given the location of residue C483 in the nsP4 palm domain, we hypothesized that other residues within this domain and surrounding subdomains would also contribute to polymerase function. To test this hypothesis, we designed a panel of nsP4 variants via homology modeling based on the Coxsackievirus B3 3 polymerase. We rescued each variant in both mammalian and mosquito cells and discovered that the palm domain and ring finger subdomain contribute to polymerase host-specific replication and genetic stability. Surprisingly, in mosquito cells, these variants in the ring finger and palm domain were replication competent and produced viral structural proteins, but they were unable to produce infectious progeny, indicating a yet uncharacterized role for the polymerase in viral assembly. Finally, we have identified additional residues in the nsP4 palm domain that influence the genetic diversity of the viral progeny, potentially via an alteration in NTP binding and/or discrimination by the polymerase. Taken together, these studies highlight that distinct nsP4 subdomains regulate multiple processes of the alphavirus life cycle, placing nsP4 in a central role during the switch from RNA synthesis to packaging and assembly.

Inter and intra-host diversity of RSV in hematopoietic stem cell transplant adults with normal and delayed viral clearance

Respiratory syncytial virus (RSV) infection in immunocompromised individuals often leads to prolonged illness, progression to severe lower respiratory tract infection, and even death. How the host immune environment of the hematopoietic stem cell transplant (HCT) adults can affect viral genetic variation during an acute infection is not understood well. In the present study, we performed whole genome sequencing of RSV/A or RSV/B from samples collected longitudinally from HCT adults with normal (<14 days) and delayed (≥14 days) RSV clearance who were enrolled in a ribavirin trial. We determined the inter-host and intra-host genetic variation of RSV and the effect of mutations on putative glycosylation sites. The inter-host variation of RSV is centered in the attachment (G) and fusion (F) glycoprotein genes followed by polymerase (L) and matrix (M) genes. Interestingly, the overall genetic variation was constant between normal and delayed clearance groups for both RSV/A and RSV/B. Intra-host variation primarily occurred in the G gene followed by non-structural protein (NS1) and L genes; however, gain or loss of stop codons and frameshift mutations appeared only in the G gene and only in the delayed viral clearance group. Potential gain or loss of O-linked glycosylation sites in the G gene occurred both in RSV/A and RSV/B isolates. For RSV F gene, loss of N-linked glycosylation site occurred in three RSV/B isolates within an antigenic epitope. Both oral and aerosolized ribavirin did not cause any mutations in the L gene. In summary, prolonged viral shedding and immune deficiency resulted in RSV variation, especially in structural mutations in the G gene, possibly associated with immune evasion. Therefore, sequencing and monitoring of RSV isolates from immunocompromised patients are crucial as they can create escape mutants that can impact the effectiveness of upcoming vaccines and treatments.

Persistent Coxsackievirus B3 Infection in Pancreatic Ductal Cells In Vitro Downregulates Cellular Polyamine Metabolism

Picornaviruses infect a wide variety of cell types in vitro, with rapid replication kinetics and pronounced cytopathic effect. Coxsackievirus B3 (CVB3) can also establish a persistent infection in vivo that can lead to pathology, including dilated cardiomyopathy and myocarditis. One model system to study persistent infection is the pancreatic ductal cell line PANC-1, which CVB3 infects and is maintained indefinitely. We have characterized this model for CVB3 infection to study persistent infection for over 6 months. We find that CVB3 rapidly replicates within PANC-1 cells without robust cytopathic effect, and after 1 month in culture, titers stabilize. We find that infection does not significantly affect cellular viability. Persistent virus reverts to lytic infection when transferred to Huh7 or Vero cells. We find that persistent CVB3 adapts to PANC-1 cells via mutation of its capsid proteins and, curiously, the viral polymerase (3Dpol) to generate a high-fidelity polymerase. Persistent infection is associated with reduced cleavage of eIF4G, reduced plaque size, and decreasing particle infectivity. We further find that polyamine metabolism is altered in persistently infected cells, with the rate-limiting enzyme ornithine decarboxylase (ODC1) reduced in translation. We further find that targeting polyamine synthesis reduces persistent infection without affecting the viability of the PANC-1 cells. Finally, we find that viral fidelity is essential to maintaining CVB3 infection, and targeting viral fidelity reduces persistent virus infection. Together, these data highlight a novel role for polyamines and fidelity in persistent CVB3 infection and suggest avenues for therapeutic development to target persistent infection. IMPORTANCE Enteroviruses are significant human pathogens that can cause severe disease, including cardiomyopathies. Viruses like coxsackievirus B3 (CVB3) can cause tissue damage by lytically infecting cells; however, CVB3 can also persistently infect, which has been associated with several pathologies. Studying persistent infection in vitro is challenging, as CVB3 lytically infects most cellular model systems. Here, we show that CVB3 establishes persistent infection in pancreatic ductal cells in vitro, similar to prior studies on other coxsackieviruses. We also show that this infection results in adaptation of the virus to these cells, as well as changes to cellular metabolism of polyamines.

The E2 glycoprotein holds key residues for Mayaro virus adaptation to the urban Aedes aegypti mosquito

Adaptation to mosquito vectors suited for transmission in urban settings is a major driver in the emergence of arboviruses. To better anticipate future emergence events, it is crucial to assess their potential to adapt to new vector hosts. In this work, we used two different experimental evolution approaches to study the adaptation process of an emerging alphavirus, Mayaro virus (MAYV), to Ae. aegypti, an urban mosquito vector of many other arboviruses. We identified E2-T179N as a key mutation increasing MAYV replication in insect cells and enhancing transmission after escaping the midgut of live Ae. aegypti. In contrast, this mutation decreased viral replication and binding in human fibroblasts, a primary cellular target of MAYV in humans. We also showed that MAYV E2-T179N generates reduced viremia and displays less severe tissue pathology in vivo in a mouse model. We found evidence in mouse fibroblasts that MAYV E2-T179N is less dependent on the Mxra8 receptor for replication than WT MAYV. Similarly, exogenous expression of human apolipoprotein receptor 2 and Mxra8 enhanced WT MAYV replication compared to MAYV E2-T179N. When this mutation was introduced in the closely related chikungunya virus, which has caused major outbreaks globally in the past two decades, we observed increased replication in both human and insect cells, suggesting E2 position 179 is an important determinant of alphavirus host-adaptation, although in a virus-specific manner. Collectively, these results indicate that adaptation at the T179 residue in MAYV E2 may result in increased vector competence-but coming at the cost of optimal replication in humans-and may represent a first step towards a future emergence event.

Genetic characterization and molecular epidemiology of Coxsackievirus A12 from mainland China during 2010-2019

Coxsackievirus A12 (CVA12) is an enterovirus that has been isolated in many countries in recent years. However, studies on CVA12 are limited, and its effective population size, evolutionary dynamics and recombination patterns have not been clarified now. In this study, we described the phylogenetic characteristics of 16 CVA12 strains isolated from pediatric HFMD patients in mainland China from 2010 to 2019. Comparison of the nucleotide sequences and amino acid sequences with the CVA12 prototype strain revealed that the 16 CVA12 strains are identical in 78.8-79% and 94-94.2%, respectively. A phylodynamic analysis based on the 16 full-length VP1 sequences from this study and 21 sequences obtained from GenBank revealed a mean substitution rate of 6.61 × 10^-3 substitutions/site/year (95% HPD: 5.16-8.20 × 10^-3), dating the time to most recent common ancestor (tMRCA) of CVA12 back to 1946 (95% HPD: 1942-1947). The Bayesian skyline plot showed that the effective population size has experienced twice dynamic fluctuations since 2007. Phylogeographic analysis identified two significant migration pathways, indicating the existence of cross-provincial transmission of CVA12 in mainland China. Recombination analysis revealed two recombination patterns between 16 CVA12 strains and other EV-A, suggesting that there may be extensive genetic exchange between CVA12 and other enteroviruses. In summary, a total of 16 full-length CVA12 strains were reported in this study, providing valuable references for further studies of CVA12 worldwide.

RNA-Dependent RNA Polymerase of the Second Human Pegivirus Exhibits a High-Fidelity Feature

The virus-encoded RNA-dependent RNA polymerase (RdRp) is responsible for viral replication, and its fidelity is closely related to viral diversity, pathogenesis, virulence, and fitness. Hepatitis C virus (HCV) and the second human pegivirus (HPgV-2) belong to the family Flaviviridae and share some features, including similar viral genome structure. Unlike HCV, HPgV-2 preserves a highly conserved genome sequence and low intrahost variation. However, the underlying mechanism remains to be elucidated. In this study, we evaluated the fidelity of HPgV-2 and HCV RdRp in an in vitro RNA polymerase reaction system. The results showed higher fidelity of HPgV-2 RdRp than HCV NS5B with respect to the misincorporation rate due to their difference in recognizing nucleoside triphosphate (NTP) substrates. Furthermore, HPgV-2 RdRp showed lower sensitivity than HCV to sofosbuvir, a nucleotide inhibitor against HCV RdRp, which explained the insusceptibility of HPgV-2 to direct-acting antiviral (DAA) therapy against HCV infection. Our results indicate that HPgV-2 could be an excellent model for studying the mechanisms involved in viral polymerase fidelity as well as RNA virus diversity and evolution. IMPORTANCE RNA viruses represent the most important pathogens for humans and animals and exhibit rapid evolution and high adaptive capacity, which is due to the high mutation rates for using the error-prone RNA-dependent RNA polymerase (RdRp) during replication. The fidelity of RdRp is closely associated with viral diversity, fitness, and pathogenesis. Previous studies have shown that the second human pegivirus (HPgV-2) exhibits a highly conserved genome sequence and low intrahost variation, which might be due to the fidelity of HPgV-2 RdRp. In this work, we used a series of in vitro RNA polymerase assays to evaluate the in vitro fidelity of HPgV-2 RdRp and compared it with that of HCV RdRp. The results indicated that HPgV-2 RdRp preserves significantly higher fidelity than HCV RdRp, which might contribute to the conservation of the HPgV-2 genome. The unique feature of HPgV-2 RdRp fidelity provides a new model for investigation of viral RdRp fidelity.

The SARS-CoV nsp12 Polymerase Active Site Is Tuned for Large-Genome Replication

Positive-strand RNA viruses replicate their genomes using virally encoded RNA-dependent RNA polymerases (RdRP) with a common active-site structure and closure mechanism upon which replication speed and fidelity can evolve to optimize virus fitness. Coronaviruses (CoV) form large multicomponent RNA replication-transcription complexes containing a core RNA synthesis machine made of the nsp12 RdRP protein with one nsp7 and two nsp8 proteins as essential subunits required for activity. We show that assembly of this complex can be accelerated 5-fold by preincubation of nsp12 with nsp8 and further optimized with the use of a novel nsp8L7 heterodimer fusion protein construct. Using rapid kinetics methods, we measure elongation rates of up to 260 nucleotides (nt)/s for the core replicase, a rate that is unusually fast for a viral polymerase. To address the origin of this fast rate, we examined the roles of two CoV-specific residues in the RdRP active site: Ala547, which replaces a conserved glutamate above the bound NTP, and Ser759, which mutates the palm domain GDD sequence to SDD. Our data show that Ala547 allows for a doubling of replication rate, but this comes at a fidelity cost that is mitigated by using a SDD sequence in the palm domain. Our biochemical data suggest that fixation of mutations in polymerase motifs F and C played a key role in nidovirus evolution by tuning replication rate and fidelity to accommodate their large genomes. IMPORTANCE Replicating large genomes represents a challenge for RNA viruses because fast RNA synthesis is needed to escape innate immunity defenses, but faster polymerases are inherently low-fidelity enzymes. Nonetheless, the coronaviruses replicate their ≈30-kb genomes using the core polymerase structure and mechanism common to all positive-strand RNA viruses. The classic explanation for their success is that the large-genome nidoviruses have acquired an exonuclease-based repair system that compensates for the high polymerase mutation rate. In this work, we establish that the nidoviral polymerases themselves also play a key role in maintaining genome integrity via mutations at two key active-site residues that enable very fast replication rates while maintaining typical mutation rates. Our findings further demonstrate the evolutionary plasticity of the core polymerase platform by showing how it has adapted during the expansion from short-genome picornaviruses to long-genome nidoviruses.

The Second Human Pegivirus, a Non-Pathogenic RNA Virus with Low Prevalence and Minimal Genetic Diversity

The second human pegivirus (HPgV-2) is a virus discovered in the plasma of a hepatitis C virus (HCV)-infected patient in 2015 belonging to the pegiviruses of the family Flaviviridae. HPgV-2 has been proved to be epidemiologically associated with and structurally similar to HCV but unrelated to HCV disease and non-pathogenic, but its natural history and tissue tropism remain unclear. HPgV-2 is a unique RNA virus sharing the features of HCV and the first human pegivirus (HPgV-1 or GBV-C). Moreover, distinct from most RNA viruses such as HCV, HPgV-1 and human immunodeficiency virus (HIV), HPgV-2 exhibits much lower genomic diversity, with a high global sequence identity ranging from 93.5 to 97.5% and significantly lower intra-host variation than HCV. The mechanisms underlying the conservation of the HPgV-2 genome are not clear but may include efficient innate immune responses, low immune selection pressure and, possibly, the unique features of the viral RNA-dependent RNA polymerase (RdRP). In this review, we summarize the prevalence, pathogenicity and genetic diversity of HPgV-2 and discuss the possible reasons for the uniformity of its genome sequence, which should elucidate the implications of RNA virus fidelity for attenuated viral vaccines.

Intra-Host Diversity of Dengue Virus in Mosquito Vectors

Dengue virus (DENV) is the most common arbovirus, causing a significant burden on both the economy and global healthcare systems. The virus is transmitted by Aedes species of mosquitoes as a swarm of closely related virus genomes, collectively referred to as a quasispecies. The level of genomic diversity within this quasispecies varies as DENV moves through various ecological niches within its transmission cycle. Here, the factors that influence the level of DENV quasispecies diversity during the course of infection in the mosquito vectors are reviewed.

The Role of the Flavivirus Replicase in Viral Diversity and Adaptation

Flaviviruses include several emerging and re-emerging arboviruses which cause millions of infections each year. Although relatively well-studied, much remains unknown regarding the mechanisms and means by which these viruses readily alternate and adapt to different hosts and environments. Here, we review a subset of the different aspects of flaviviral biology which impact host switching and viral fitness. These include the mechanism of replication and structural biology of the NS3 and NS5 proteins, which reproduce the viral genome; rates of mutation resulting from this replication and the role of mutational frequency in viral fitness; and the theory of quasispecies evolution and how it contributes to our understanding of genetic and phenotypic plasticity.

Influence of Ribavirin on Mumps Virus Population Diversity

Frequent mumps outbreaks in vaccinated populations and the occurrence of neurological complications (e.g., aseptic meningitis or encephalitis) in patients with mumps indicate the need for the development of more efficient vaccines as well as specific antiviral therapies. RNA viruses are genetically highly heterogeneous populations that exist on the edge of an error threshold, such that additional increases in mutational burden can lead to extinction of the virus population. Deliberate modulation of their natural mutation rate is being exploited as an antiviral strategy and a possibility for rational vaccine design. The aim of this study was to examine the ability of ribavirin, a broad-spectrum antiviral agent, to introduce mutations in the mumps virus (MuV) genome and to investigate if resistance develops during long-term in vitro exposure to ribavirin. An increase in MuV population heterogeneity in the presence of ribavirin has been observed after one passage in cell culture, as well as a bias toward C-to-U and G-to-A transitions, which have previously been defined as ribavirin-related. At higher ribavirin concentration, MuV loses its infectivity during serial passaging and does not recover. At low ribavirin concentration, serial passaging leads to a more significant increase in population diversity and a stronger bias towards ribavirin-related transitions, independently of viral strain or cell culture. In these conditions, the virus retains its initial growth capacity, without development of resistance at a whole-virus population level.

Application Route and Immune Status of the Host Determine Safety and Oncolytic Activity of Oncolytic Coxsackievirus B3 Variant PD-H

The coxsackievirus B3 strain PD-0 has been proposed as a new oncolytic virus for the treatment of colorectal carcinoma. Here, we generated a cDNA clone of PD-0 and analyzed the virus PD-H, newly generated from this cDNA, in xenografted and syngenic models of colorectal cancer. Replication and cytotoxic assays revealed that PD-H replicated and lysed colorectal carcinoma cell lines in vitro as well as PD-0. Intratumoral injection of PD-H into subcutaneous DLD-1 tumors in nude mice resulted in strong inhibition of tumor growth and significantly prolonged the survival of the animals, but virus-induced systemic infection was observed in one of the six animals. In a syngenic mouse model of subcutaneously growing Colon-26 tumors, intratumoral administration of PD-H led to a significant reduction of tumor growth, the prolongation of animal survival, the prevention of tumor-induced cachexia, and the elevation of CD3+ and dendritic cells in the tumor microenvironment. No virus-induced side effects were observed. After intraperitoneal application, PD-H induced weak pancreatitis and myocarditis in immunocompetent mice. By equipping the virus with target sites of miR-375, which is specifically expressed in the pancreas, organ infections were prevented. Moreover, employment of this virus in a syngenic mouse model of CT-26 peritoneal carcinomatosis resulted in a significant reduction in tumor growth and an increase in animal survival. The results demonstrate that the immune status of the host, the route of virus application, and the engineering of the virus with target sites of suitable microRNAs are crucial for the use of PD-H as an oncolytic virus.

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.

Imprecise recombinant viruses evolve via a fitness-driven, iterative process of polymerase template-switching events

Recombination is a common feature of many positive-strand RNA viruses, playing an important role in virus evolution. However, to date, there is limited understanding of the mechanisms behind the process. Utilising in vitro assays, we have previously shown that the template-switching event of recombination is a random and ubiquitous process that often leads to recombinant viruses with imprecise genomes containing sequence duplications. Subsequently, a process termed resolution, that has yet to be mechanistically studied, removes these duplicated sequences resulting in a virus population of wild type length genomes. Using defined imprecise recombinant viruses together with Oxford Nanopore and Illumina high throughput next generation sequencing technologies we have investigated the process of resolution. We show that genome resolution involves subsequent rounds of template-switching recombination with viral fitness resulting in the survival of a small subset of recombinant genomes. This alters our previously held understanding that recombination and resolution are independent steps of the process, and instead demonstrates that viruses undergo frequent and continuous recombination events over a prolonged period until the fittest viruses, predominantly those with wild type length genomes, dominate the population.

Viruses with positive-sense RNA genomes, such as poliovirus, have several mechanisms by which they evolve. One of these is the process of recombination involving the large-scale exchange of genetic information. Recombination occurs during replication when the viral polymerase, bound to the nascent RNA chain, switches from copying one genome to another. However, the polymerase does not always accurately switch between the two, resulting in sequence duplications or deletions, and genomes that are referred to as imprecise. Over multiple rounds of replication sequence duplications are lost and genomes are resolved to wild type length, but it is unclear how this occurs. Here we used synthetic polioviruses containing defined sequence duplications to determine that the genome population undergoes repeated rounds of recombination until sequence duplications are lost and viruses with precise, wild type length genomes are selected for. This selection is based on the overall fitness of the virus population, with less fit imprecise viruses evolving more quickly. Our study suggests that recombination is a continual process where virus fitness drives the selection of a small subset of recombinant variants. These data are important for understanding how novel viruses evolve via recombination and how this process can be blocked to prevent novel and dangerous pathogens from arising.

Tradeoffs for a viral mutant with enhanced replication speed

RNA viruses exist as genetically heterogeneous populations due to high mutation rates, and many of these mutations reduce fitness and/or replication speed. However, it is unknown whether mutations can increase replication speed of a virus already well adapted to replication in cultured cells. By sequentially passaging coxsackievirus B3 in cultured cells and collecting the very earliest progeny, we selected for increased replication speed. We found that a single mutation in a viral capsid protein, VP1-F106L, was sufficient for the fast-replication phenotype. Characterization of this mutant revealed quicker genome release during entry compared to wild-type virus, highlighting a previously unappreciated infection barrier. However, this mutation also reduced capsid stability in vitro and reduced replication and pathogenesis in mice. These results reveal a tradeoff between overall replication speed and fitness. Importantly, this approach-selecting for the earliest viral progeny-could be applied to a variety of viral systems and has the potential to reveal unanticipated inefficiencies in viral replication cycles.

Allosteric and dynamic control of RNA-dependent RNA polymerase function and fidelity

All RNA viruses encode an RNA-dependent RNA polymerase (RdRp) responsible for genome replication. It is now recognized that enzymes in general, and RdRps specifically, are dynamic macromolecular machines such that their moving parts, including active site loops, play direct functional roles. While X-ray crystallography has provided deep insight into structural elements important for RdRp function, this methodology generally provides only static snapshots, and so is limited in its ability to report on dynamic fluctuations away from the lowest energy conformation. Nuclear magnetic resonance (NMR), molecular dynamics (MD) simulations and other biophysical techniques have brought new insight into RdRp function by their ability to characterize the trajectories, kinetics and thermodynamics of conformational motions. In particular, these methodologies have identified coordinated motions among conserved structural motifs necessary for nucleotide selection and incorporation. Disruption of these motions through amino acid substitutions or inhibitor binding impairs RdRp function. Understanding and re-engineering these motions thus provides exciting new avenues for anti-viral strategies. This chapter outlines the basics of these methodologies, summarizes the dynamic motions observed in different RdRps important for nucleotide selection and incorporation, and illustrates how this information can be leveraged towards rational vaccine strain development and anti-viral drug design.

Genome Evolution of Two Genetically Homogeneous Infectious Bursal Disease Virus Strains During Passages in vitro and ex vivo in the Presence of a Mutagenic Nucleoside Analog

The avibirnavirus infectious bursal disease virus (IBDV) is responsible for a highly contagious and sometimes lethal disease of chickens (Gallus gallus). IBDV genetic variation is well-described for both field and live-attenuated vaccine strains, however, the dynamics and selection pressures behind this genetic evolution remain poorly documented. Here, genetically homogeneous virus stocks were generated using reverse genetics for a very virulent strain, rvv, and a vaccine-related strain, rCu-1. These viruses were serially passaged at controlled multiplicities of infection in several biological systems, including primary chickens B cells, the main cell type targeted by IBDV in vivo. Passages were also performed in the absence or presence of a strong selective pressure using the antiviral nucleoside analog 7-deaza-2'-C-methyladenosine (7DMA). Next Generation Sequencing (NGS) of viral genomes after the last passage in each biological system revealed that (i) a higher viral diversity was generated in segment A than in segment B, regardless 7DMA treatment and viral strain, (ii) diversity in segment B was increased by 7DMA treatment in both viruses, (iii) passaging of IBDV in primary chicken B cells, regardless of 7DMA treatment, did not select cell-culture adapted variants of rvv, preserving its capsid protein (VP2) properties, (iv) mutations in coding and non-coding regions of rCu-1 segment A could potentially associate to higher viral fitness, and (v) a specific selection, upon 7DMA addition, of a Thr329Ala substitution occurred in the viral polymerase VP1. The latter change, together with Ala270Thr change in VP2, proved to be associated with viral attenuation in vivo. These results identify genome sequences that are important for IBDV evolution in response to selection pressures. Such information will help tailor better strategies for controlling IBDV infection in chickens.

An MHV-68 Mutator Phenotype Mutant Virus, Confirmed by CRISPR/Cas9-Mediated Gene Editing of the Viral DNA Polymerase Gene, Shows Reduced Viral Fitness

Drug resistance studies on human γ-herpesviruses are hampered by the absence of an in vitro system that allows efficient lytic viral replication. Therefore, we employed murine γ-herpesvirus-68 (MHV-68) that efficiently replicates in vitro as a model to study the antiviral resistance of γ-herpesviruses. In this study, we investigated the mechanism of resistance to nucleoside (ganciclovir (GCV)), nucleotide (cidofovir (CDV), HPMP-5azaC, HPMPO-DAPy) and pyrophosphate (foscarnet (PFA)) analogues and the impact of these drug resistance mutations on viral fitness. Viral fitness was determined by dual infection competition assays, where MHV-68 drug-resistant viral clones competed with the wild-type virus in the absence and presence of antivirals. Using next-generation sequencing, the composition of the viral populations was determined at the time of infection and after 5 days of growth. Antiviral drug resistance selection resulted in clones harboring mutations in the viral DNA polymerase (DP), denoted Y383S^GCV, Q827R^HPMP-5azaC, G302W^PFA, K442T^PFA, G302W+K442T^PFA, C297W^HPMPO-DAPy and C981Y^CDV. Without antiviral pressure, viral clones Q827R^HPMP-5azaC, G302W^PFA, K442T^PFA and G302W+K442T^PFA grew equal to the wild-type virus. However, in the presence of antivirals, these mutants had a growth advantage over the wild-type virus that was moderately to very strongly correlated with antiviral resistance. The Y383S^GCV mutant was more fit than the wild-type virus with and without antivirals, except in the presence of brivudin. The C297W and C981Y changes were associated with a mutator phenotype and had a severely impaired viral fitness in the absence and presence of antivirals. The mutator phenotype caused by C297W in MHV-68 DP was validated by using a CRISPR/Cas9 genome editing approach.

Development of a Genetically Stable Live Attenuated Influenza Vaccine Strain Using an Engineered High-Fidelity Viral Polymerase

RNA viruses demonstrate a vast range of variants, called quasispecies, due to error-prone replication by viral RNA-dependent RNA polymerase. Although live attenuated vaccines are effective in preventing RNA virus infection, there is a risk of reversal to virulence after their administration. To test the hypothesis that high-fidelity viral polymerase reduces the diversity of influenza virus quasispecies, resulting in inhibition of reversal of the attenuated phenotype, we first screened for a high-fidelity viral polymerase using serial virus passages under selection with a guanosine analog ribavirin. Consequently, we identified a Leu66-to-Val single amino acid mutation in polymerase basic protein 1 (PB1). The high-fidelity phenotype of PB1-L66V was confirmed using next-generation sequencing analysis and biochemical assays with the purified influenza viral polymerase. As expected, PB1-L66V showed at least two-times-lower mutation rates and decreased misincorporation rates, compared to the wild type (WT). Therefore, we next generated an attenuated PB1-L66V virus with a temperature-sensitive (ts) phenotype based on FluMist, a live attenuated influenza vaccine (LAIV) that can restrict virus propagation by ts mutations, and examined the genetic stability of the attenuated PB1-L66V virus using serial virus passages. The PB1-L66V mutation prevented reversion of the ts phenotype to the WT phenotype, suggesting that the high-fidelity viral polymerase could contribute to generating an LAIV with high genetic stability, which would not revert to the pathogenic virus.IMPORTANCE The LAIV currently in use is prescribed for actively immunizing individuals aged 2 to 49 years. However, it is not approved for infants and elderly individuals, who actually need it the most, because it might prolong virus propagation and cause an apparent infection in these individuals, due to their weak immune systems. Recently, reversion of the ts phenotype of the LAIV strain currently in use to a pathogenic virus was demonstrated in cultured cells. Thus, the generation of mutations associated with enhanced virulence in LAIV should be considered. In this study, we isolated a novel influenza virus strain with a Leu66-to-Val single amino acid mutation in PB1 that displayed a significantly higher fidelity than the WT. We generated a novel LAIV candidate strain harboring this mutation. This strain showed higher genetic stability and no ts phenotype reversion. Thus, our high-fidelity strain might be useful for the development of a safer LAIV.

Attenuated dengue viruses are genetically more diverse than their respective wild-type parents

Dengue poses a significant burden of individual health, health systems and the economy in dengue endemic regions. As such, dengue vaccine development has been an active area of research. Previous studies selected attenuated vaccine candidates based on plaque size. However, these candidates led to mixed safety outcome in clinical trials, suggesting it is insufficiently informative as an indicator of dengue virus (DENV) attenuation. In this study, we examined the genome diversity of wild-type DENVs and their attenuated derivatives developed by Mahidol University and tested in phase 1 clinical trials. We found that the attenuated DENVs, in particular the strain under clinical development by Takeda Vaccines, DENV2 PDK53, showed significantly higher genome diversity than its wild-type parent, DENV2 16681. The determinant of genomic diversity was intrinsic to the PDK53 genome as infectious clone of PDK53 showed greater genomic diversity after a single in vitro passage compared to 16681 infectious clone. Similar trends were observed with attenuated DENV1 and DENV4, both of which were shown to be attenuated clinically, but not DENV3 that was not adequately attenuated clinically. Taken together, evidence presented here suggests that genome diversity could be developed into a marker of DENV attenuation.

Increased RNA virus population diversity improves adaptability

The replication machinery of most RNA viruses lacks proofreading mechanisms. As a result, RNA virus populations harbor a large amount of genetic diversity that confers them the ability to rapidly adapt to changes in their environment. In this work, we investigate whether further increasing the initial population diversity of a model RNA virus can improve adaptation to a single selection pressure, thermal inactivation. For this, we experimentally increased the diversity of coxsackievirus B3 (CVB3) populations across the capsid region. We then compared the ability of these high diversity CVB3 populations to achieve resistance to thermal inactivation relative to standard CVB3 populations in an experimental evolution setting. We find that viral populations with high diversity are better able to achieve resistance to thermal inactivation at both the temperature employed during experimental evolution as well as at a more extreme temperature. Moreover, we identify mutations in the CVB3 capsid that confer resistance to thermal inactivation, finding significant mutational epistasis. Our results indicate that even naturally diverse RNA virus populations can benefit from experimental augmentation of population diversity for optimal adaptation and support the use of such viral populations in directed evolution efforts that aim to select viruses with desired characteristics.

Polymerase Fidelity Contributes to Foot-and-Mouth Disease Virus Pathogenicity and Transmissibility In Vivo

The low fidelity of foot-and-mouth disease virus (FMDV) RNA-dependent RNA polymerase allows FMDV to exhibit high genetic diversity. Previously, we showed that the genetic diversity of FMDV plays an important role in virulence in suckling mice. Here, we mutated the amino acid residue Phe257, located in the finger domain of FMDV polymerase and conserved across FMDV serotypes, to a cysteine (F257C) to study the relationship between viral genetic diversity, virulence, and transmissibility in natural hosts. The single amino acid substitution in FMDV polymerase resulted in a high-fidelity virus variant, rF257C, with growth kinetics indistinguishable from those of wild-type (WT) virus in cell culture, but it displayed smaller plaques and impaired fitness in direct competition assays. Furthermore, we found that rF257C was attenuated in vivo in both suckling mice and pigs (one of its natural hosts). Importantly, contact exposure experiments showed that the rF257C virus exhibited reduced transmissibility compared to that of wild-type FMDV in the porcine model. This study provides evidence that FMDV genetic diversity is important for viral virulence and transmissibility in susceptible animals. Given that type O FMDV exhibits the highest genetic diversity among all seven serotypes of FMDV, we propose that the lower polymerase fidelity of the type O FMDV could contribute to its dominance worldwide.IMPORTANCE Among the seven serotypes of FMDV, serotype O FMDV have the broadest distribution worldwide, which could be due to their high virulence and transmissibility induced by high genetic diversity. In this paper, we generated a single amino acid substitution FMDV variant with a high-fidelity polymerase associated with viral fitness, virulence, and transmissibility in a natural host. The results highlight that maintenance of viral population diversity is essential for interhost viral spread. This study provides evidence that higher genetic diversity of type O FMDV could increase both virulence and transmissibility, thus leading to their dominance in the global epidemic.

Genetic Adaptation of Coxsackievirus B1 during Persistent Infection in Pancreatic Cells

Coxsackie B (CVB) viruses have been associated with type 1 diabetes. We have recently observed that CVB1 was linked to the initiation of the autoimmune process leading to type 1 diabetes in Finnish children. Viral persistency in the pancreas is currently considered as one possible mechanism. In the current study persistent infection was established in pancreatic ductal and beta cell lines (PANC-1 and 1.1B4) using four different CVB1 strains, including the prototype strain and three clinical isolates. We sequenced 5′ untranslated region (UTR) and regions coding for structural and non-structural proteins and the second single open reading frame (ORF) protein of all persisting CVB1 strains using next generation sequencing to identify mutations that are common for all of these strains. One mutation, K257R in VP1, was found from all persisting CVB1 strains. The mutations were mainly accumulated in viral structural proteins, especially at BC, DE, EF loops and C-terminus of viral capsid protein 1 (VP1), the puff region of VP2, the knob region of VP3 and infection-enhancing epitope of VP4. This showed that the capsid region of the viruses sustains various changes during persistency some of which could be hallmark(s) of persistency.

Residues within the Foot-and-Mouth Disease Virus 3Dpol Nuclear Localization Signal Affect Polymerase Fidelity

Many RNA viruses encode a proof-reading deficient, low-fidelity RNA-dependent polymerase (RdRp), which generates genetically diverse populations that can adapt to changing environments and thwart antiviral therapies. 3D^pol, the RdRp of the foot-and-mouth disease virus (FMDV), is responsible for replication of viral genomes. The 3D^pol N terminus encodes a nuclear localization signal (NLS) sequence,MRKTKLAPT, important for import of the protein to host nucleus. Previous studies showed that substitutions at residues 18 and 20 of the NLS are defective in proper incorporation of nucleotides and RNA binding. Here, we use a systematic alanine scanning mutagenesis approach to understand the role of individual residues of the NLS in nuclear localization and nucleotide incorporation activities of 3D^pol We identify two residues of 3D^pol NLS, T19 and L21, that are important for the maintenance of enzyme fidelity. The 3D^pol NLS alanine substitutions of T19 and L21 results in aberrant incorporation of nucleoside analogs, conferring a low fidelity phenotype of the enzyme. A molecular dynamics simulation of RNA- and mutagen (RTP)-bound 3D^pol revealed that the T19 residue participates in a hydrogen bond network, including D165 in motif F and R416 at the C terminus of the FMDV 3D^pol and RNA template-primer. Based on these findings and previous studies, we conclude that at least the first six residues of theMRKTKLAPT sequence motif play a vital role in the maintenance of faithful RNA synthesis activity (fidelity) of FMDV 3D^pol, suggesting that the role of the NLS motif in similar viral polymerases needs to be revisited.IMPORTANCE In this study, we employed genetic and molecular dynamics approaches to analyze the role of individual amino acids of the FMDV 3D^pol nuclear localization signal (NLS). The NLS residues were mutated to alanine using a type A full-genome cDNA clone, and the virus progeny was analyzed for defects in growth and in competition with the parental virus. We identified two mutants in 3D^pol, T19A and L21A, that exhibited high rate of mutation, were sensitive to nucleotide analogs, and displayed reduced replicative fitness compared to the parental virus. Using molecular dynamics simulation, we demonstrated that residues T19 and L21 played a role in the structural configuration of the interaction network at the 3D^pol palm subdomain. Cumulatively, our data suggest that the T19 and L21 3D^pol amino acids are important for maintaining the fidelity of the FMDV polymerase and ensuring faithful replication of the FMDV genome.

Picornaviral polymerase domain exchanges reveal a modular basis for distinct biochemical activities of viral RNA-dependent RNA polymerases

Picornaviral RNA-dependent RNA polymerases (RdRPs) have low replication fidelity that is essential for viral fitness and evolution. Their global fold consists of the classical "cupped right hand" structure with palm, fingers, and thumb domains, and these RdRPs also possess a unique contact between the fingers and thumb domains. This interaction restricts movements of the fingers, and RdRPs use a subtle conformational change within the palm domain to close their active sites for catalysis. We have previously shown that this core RdRP structure and mechanism provide a platform for polymerases to fine-tune replication rates and fidelity to optimize virus fitness. Here, we further elucidated the structural basis for differences in replication rates and fidelity among different viruses by generating chimeric RdRPs from poliovirus and coxsackievirus B3. We designed these chimeric polymerases by exchanging the fingers, pinky finger, or thumb domains. The results of biochemical, rapid-quench, and stopped-flow assays revealed that differences in biochemical activity map to individual modular domains of this polymerase. We found that the pinky finger subdomain is a major regulator of initiation and that the palm domain is the major determinant of catalytic rate and nucleotide discrimination. We further noted that thumb domain interactions with product RNA regulate translocation and that the palm and thumb domains coordinately control elongation complex stability. Several RdRP chimeras supported the growth of infectious poliovirus, providing insights into enterovirus species-specific protein-protein interactions required for virus replication.

A Low Fidelity Virus Shows Increased Recombination during the Removal of an Alphavirus Reporter Gene

Reporter genes for RNA viruses are well-known to be unstable due to putative RNA recombination events that excise inserted nucleic acids. RNA recombination has been demonstrated to be co-regulated with replication fidelity in alphaviruses, but it is unknown how recombination events at the minority variant level act, which is important for vaccine and trans-gene delivery design. Therefore, we sought to characterize the removal of a reporter gene by a low-fidelity alphavirus mutant over multiple replication cycles. To examine this, GFP was inserted into TC-83, a live-attenuated vaccine for the alphavirus Venezuelan equine encephalitis virus, as well as a low-fidelity variant of TC-83, and passaged until fluorescence was no longer observed. Short-read RNA sequencing using ClickSeq was performed to determine which regions of the viral genome underwent recombination and how this changed over multiple replication cycles. A rapid removal of the GFP gene was observed, where minority variants in the virus population accumulated small deletions that increased in size over the course of passaging. Eventually, these small deletions merged to fully remove the GFP gene. The removal was significantly enhanced during the passaging of low-fidelity TC-83, suggesting that increased levels of recombination are a defining characteristic of this mutant.

Within-Host Viral Diversity: A Window into Viral Evolution

The evolutionary dynamics of a virus can differ within hosts and across populations. Studies of within-host evolution provide an important link between experimental studies of virus evolution and large-scale phylodynamic analyses. They can determine the extent to which global processes are recapitulated on local scales and how accurately experimental infections model natural ones. They may also inform epidemiologic models of disease spread and reveal how host-level dynamics contribute to a virus's evolution at a larger scale. Over the last decade, advances in viral sequencing have enabled detailed studies of viral genetic diversity within hosts. I review how within-host diversity is sampled, measured, and expressed, and how comparative studies of viral diversity can be leveraged to elucidate a virus's evolutionary dynamics. These concepts are illustrated with detailed reviews of recent research on the within-host evolution of influenza virus, dengue virus, and cytomegalovirus.

Measuring Alphavirus Fidelity Using Non-Infectious Virus Particles

Mutations are incorporated into the genomes of RNA viruses at an optimal frequency and altering this precise frequency has been proposed as a strategy to create live-attenuated vaccines. However, determining the effect of specific mutations that alter fidelity has been difficult because of the rapid selection of the virus population during replication. By deleting residues of the structural polyprotein PE2 cleavage site, E3Δ56-59, in Venezuelan equine encephalitis virus (VEEV) TC-83 vaccine strain, non-infectious virus particles were used to assess the effect of single mutations on mutation frequency without the interference of selection that results from multiple replication cycles. Next-generation sequencing analysis revealed a significantly lower frequency of transversion mutations and overall mutation frequency for the fidelity mutants compared to VEEV TC-83 E3Δ56-59. We demonstrate that deletion of the PE2 cleavage site halts virus infection while making the virus particles available for downstream sequencing. The conservation of the site will allow the evaluation of suspected fidelity mutants across alphaviruses of medical importance.

A nucleobase-binding pocket in a viral RNA-dependent RNA polymerase contributes to elongation complex stability

The enterovirus 71 (EV71) 3D^pol is an RNA-dependent RNA polymerase (RdRP) that plays the central role in the viral genome replication, and is an important target in antiviral studies. Here, we report a crystal structure of EV71 3D^pol elongation complex (EC) at 1.8 Å resolution. The structure reveals that the 5′-end guanosine of the downstream RNA template interacts with a fingers domain pocket, with the base sandwiched by H44 and R277 side chains through hydrophobic stacking interactions, and these interactions are still maintained after one in-crystal translocation event induced by nucleotide incorporation, implying that the pocket could regulate the functional properties of the polymerase by interacting with RNA. When mutated, residue R277 showed an impact on virus proliferation in virological studies with residue H44 having a synergistic effect. In vitro biochemical data further suggest that mutations at these two sites affect RNA binding, EC stability, but not polymerase catalytic rate (k_cat) and apparent NTP affinity (K_M,NTP). We propose that, although rarely captured by crystallography, similar surface pocket interaction with nucleobase may commonly exist in nucleic acid motor enzymes to facilitate their processivity. Potential applications in antiviral drug and vaccine development are also discussed.

Quasispecies dynamics in disease prevention and control

Medical interventions to prevent and treat viral disease constitute evolutionary forces that may modify the genetic composition of viral populations that replicate in an infected host and influence the genomic composition of those viruses that are transmitted and progress at the epidemiological level. Given the adaptive potential of viruses in general and the RNA viruses in particular, the selection of viral mutants that display some degree of resistance to inhibitors or vaccines is a tangible challenge. Mutant selection may jeopardize control of the viral disease. Strategies intended to minimize vaccination and treatment failures are proposed and justified based on fundamental features of viral dynamics explained in the preceding chapters. The recommended use of complex, multiepitopic vaccines, and combination therapies as early as possible after initiation of infection falls under the general concept that complexity cannot be combated with simplicity. It also follows that sociopolitical action to interrupt virus replication and spread as soon as possible is as important as scientifically sound treatment designs to control viral disease on a global scale.

Molecular basis of genetic variation of viruses

Genetic variation is a necessity of all biological systems. Viruses use all known mechanisms of variation; mutation, several forms of recombination, and segment reassortment in the case of viruses with a segmented genome. These processes are intimately connected with the replicative machineries of viruses, as well as with fundamental physical-chemical properties of nucleotides when acting as template or substrate residues. Recombination has been viewed as a means to rescue viable genomes from unfit parents or to produce large modifications for the exploration of phenotypic novelty. All types of genetic variation can act conjointly as blind processes to provide the raw materials for adaptation to the changing environments in which viruses must replicate. A distinction is made between mechanistically unavoidable and evolutionarily relevant mutation and recombination.

RNA-Dependent RNA Polymerase Speed and Fidelity are not the Only Determinants of the Mechanism or Efficiency of Recombination

Using the RNA-dependent RNA polymerase (RdRp) from poliovirus (PV) as our model system, we have shown that Lys-359 in motif-D functions as a general acid in the mechanism of nucleotidyl transfer. A K359H (KH) RdRp derivative is slow and faithful relative to wild-type enzyme. In the context of the KH virus, RdRp-coding sequence evolves, selecting for the following substitutions: I331F (IF, motif-C) and P356S (PS, motif-D). We have evaluated IF-KH, PS-KH, and IF-PS-KH viruses and enzymes. The speed and fidelity of each double mutant are equivalent. Each exhibits a unique recombination phenotype, with IF-KH being competent for copy-choice recombination and PS-KH being competent for forced-copy-choice recombination. Although the IF-PS-KH RdRp exhibits biochemical properties within twofold of wild type, the virus is impaired substantially for recombination in cells. We conclude that there are biochemical properties of the RdRp in addition to speed and fidelity that determine the mechanism and efficiency of recombination. The interwoven nature of speed, fidelity, the undefined property suggested here, and recombination makes it impossible to attribute a single property of the RdRp to fitness. However, the derivatives described here may permit elucidation of the importance of recombination on the fitness of the viral population in a background of constant polymerase speed and fidelity.

2'-C-methylated nucleotides terminate virus RNA synthesis by preventing active site closure of the viral RNA-dependent RNA polymerase

The 2'-C-methyl ribonucleosides are nucleoside analogs representing an important class of antiviral agents, especially against positive-strand RNA viruses. Their value is highlighted by the highly successful anti-hepatitis C drug sofosbuvir. When appropriately phosphorylated, these nucleotides are successfully incorporated into RNA by the virally encoded RNA-dependent RNA polymerase (RdRp). This activity prevents further RNA extension, but the mechanism is poorly characterized. Previously, we had identified NMR signatures characteristic of formation of RdRp-RNA binary and RdRp-RNA-NTP ternary complexes for the poliovirus RdRp, including an open-to-closed conformational change necessary to prepare the active site for catalysis of phosphoryl transfer. Here we used these observations as a framework for interpreting the effects of 2'-C-methyl adenosine analogs on RNA chain extension in solution-state NMR spectroscopy experiments, enabling us to gain additional mechanistic insights into 2'-C-methyl ribonucleoside-mediated RNA chain termination. Contrary to what has been proposed previously, poliovirus RdRp that was bound to RNA with an incorporated 2'-C-methyl nucleotide could still bind to the next incoming NTP. Our results also indicated that incorporation of the 2'-C-methyl nucleotide does not disrupt RdRp-RNA interactions and does not prevent translocation. Instead, incorporation of the 2'-C-methyl nucleotide blocked closure of the RdRp active site upon binding of the next correct incoming NTP, which prevented further nucleotide addition. We propose that other nucleotide analogs that act as nonobligate chain terminators may operate through a similar mechanism.

Tyr82 Amino Acid Mutation in PB1 Polymerase Induces an Influenza Virus Mutator Phenotype

In various positive-sense single-stranded RNA viruses, a low-fidelity viral RNA-dependent RNA polymerase (RdRp) confers attenuated phenotypes by increasing the mutation frequency. We report a negative-sense single-stranded RNA virus RdRp mutant strain with a mutator phenotype. Based on structural data of RdRp, rational targeting of key residues, and screening of fidelity variants, we isolated a novel low-fidelity mutator strain of influenza virus that harbors a Tyr82-to-Cys (Y82C) single-amino-acid substitution in the PB1 polymerase subunit. The purified PB1-Y82C polymerase indeed showed an increased frequency of misincorporation compared with the wild-type PB1 in an in vitro biochemical assay. To further investigate the effects of position 82 on PB1 polymerase fidelity, we substituted various amino acids at this position. As a result, we isolated various novel mutators other than PB1-Y82C with higher mutation frequencies. The structural model of influenza virus polymerase complex suggested that the Tyr82 residue, which is located at the nucleoside triphosphate entrance tunnel, may influence a fidelity checkpoint. Interestingly, although the PB1-Y82C variant replicated with wild-type PB1-like kinetics in tissue culture, the 50% lethal dose of the PB1-Y82C mutant was 10 times lower than that of wild-type PB1 in embryonated chicken eggs. In conclusion, our data indicate that the Tyr82 residue of PB1 has a crucial role in regulating polymerase fidelity of influenza virus and is closely related to attenuated pathogenic phenotypes in vivo IMPORTANCE Influenza A virus rapidly acquires antigenic changes and antiviral drug resistance, which limit the effectiveness of vaccines and drug treatments, primarily owing to its high rate of evolution. Virus populations formed by quasispecies can contain resistance mutations even before a selective pressure is applied. To study the effects of the viral mutation spectrum and quasispecies, high- and low-fidelity variants have been isolated for several RNA viruses. Here, we report the discovery of a low-fidelity RdRp variant of influenza A virus that contains a substitution at Tyr82 in PB1. Viruses containing the PB1-Y82C substitution showed growth kinetics and viral RNA synthesis levels similar to those of the wild-type virus in cell culture; however, they had significantly attenuated phenotypes in a chicken egg infection experiment. These data demonstrated that decreased RdRp fidelity attenuates influenza A virus in vivo, which is a desirable feature for the development of safer live attenuated vaccine candidates.

Recombination in Enteroviruses, a Multi-Step Modular Evolutionary Process

RNA recombination is a major driving force in the evolution and genetic architecture shaping of enteroviruses. In particular, intertypic recombination is implicated in the emergence of most pathogenic circulating vaccine-derived polioviruses, which have caused numerous outbreaks of paralytic poliomyelitis worldwide. Recent experimental studies that relied on recombination cellular systems mimicking natural genetic exchanges between enteroviruses provided new insights into the molecular mechanisms of enterovirus recombination and enabled to define a new model of genetic plasticity for enteroviruses. Homologous intertypic recombinant enteroviruses that were observed in nature would be the final products of a multi-step process, during which precursor nonhomologous recombinant genomes are generated through an initial inter-genomic RNA recombination event and can then evolve into a diversity of fitter homologous recombinant genomes over subsequent intra-genomic rearrangements. Moreover, these experimental studies demonstrated that the enterovirus genome could be defined as a combination of genomic modules that can be preferentially exchanged through recombination, and enabled defining the boundaries of these recombination modules. These results provided the first experimental evidence supporting the theoretical model of enterovirus modular evolution previously elaborated from phylogenetic studies of circulating enterovirus strains. This review summarizes our current knowledge regarding the mechanisms of recombination in enteroviruses and presents a new evolutionary process that may apply to other RNA viruses.

Rational Control of Poliovirus RNA-Dependent RNA Polymerase Fidelity by Modulating Motif-D Loop Conformational Dynamics

The conserved structural motif D is an important determinant of the speed and fidelity of viral RNA-dependent RNA polymerases (RdRps). Structural and computational studies have suggested that conformational changes in the motif-D loop that help to reposition the catalytic lysine represent critical steps in nucleotide selection and incorporation. Conformations of the motif-D loop in the poliovirus RdRp are likely controlled in part by noncovalent interactions involving the motif-D residue Glu364. This residue swivels between making interactions with Lys228 and Asn370 to stabilize the open and closed loop conformations, respectively. We show here that we can rationally control the motif-D loop conformation by breaking these interactions. The K228A variant favors a more active closed conformation, leading to increased nucleotide incorporation rates and decreased nucleotide selectivity, and the N370A variant favors a less active open conformation, leading to decreased nucleotide incorporation rates and increased nucleotide selectivity. Similar competing interactions likely control nucleotide incorporation rates and fidelity in other viral RdRps. Rational engineering of these interactions may be important in the generation of live, attenuated vaccine strains, considering the established relationships between RdRp function and viral pathogenesis.

Evolutionary and Structural Overview of Human Picornavirus Capsid Antibody Evasion

Picornaviruses constitute one of the most relevant viral groups according to their impact on human and animal health. Etiologic agents of a broad spectrum of illnesses with a clinical presentation that ranges from asymptomatic to fatal disease, they have been the cause of uncountable epidemics throughout history. Picornaviruses are small naked RNA-positive single-stranded viruses that include some of the most important pillars in the development of virology, comprising poliovirus, rhinovirus, and hepatitis A virus. Picornavirus infectious particles use the fecal-oral or respiratory routes as primary modes of transmission. In this regard, successful viral spread relies on the capability of viral capsids to (i) shelter the viral genome, (ii) display molecular determinants for cell receptor recognition, (iii) facilitate efficient genome delivery, and (iv) escape from the immune system. Importantly, picornaviruses display a substantial amount of genetic variability driven by both mutation and recombination. Therefore, the outcome of their replication results in the emergence of a genetically diverse cloud of individuals presenting phenotypic variance. The host humoral response against the capsid protein represents the most active immune pressure and primary weapon to control the infection. Since the preservation of the capsid function is deeply rooted in the virus evolutionary dynamics, here we review the current structural evidence focused on capsid antibody evasion mechanisms from that perspective.

Naturally occurring mutations in PB1 affect influenza A virus replication fidelity, virulence, and adaptability

Background

Mutations in the PB1 subunit of RNA-dependent RNA polymerase (RdRp) of influenza A virus can affect replication fidelity. Before the influenza A/H1N1 pandemic in 2009, most human influenza A/H1N1 viruses contained the avian-associated residue, serine, at position 216 in PB1. However, near the onset of the 2009 pandemic, human viruses began to acquire the mammalian-associated residue, glycine, at PB1–216, and PB1–216G became predominant in human viruses thereafter.

Methods

Using entropy-based analysis algorithm, we have previously identified several host-specific amino-acid signatures that separated avian and swine viruses from human influenza viruses. The presence of these host-specific signatures in human influenza A/H1N1 viruses suggested that these mutations were the result of adaptive genetic evolution that enabled these influenza viruses to circumvent host barriers, which resulted in cross-species transmission. We investigated the biological impact of this natural avian-to-mammalian signature substitution at PB1–216 in human influenza A/H1N1 viruses.

Results

We found that PB1–216G viruses had greater mutation potential, and were more sensitive to ribavirin than PB1–216S viruses. In oseltamivir-treated HEK293 cells, PB1–216G viruses generated mutations in viral neuraminidase at a higher rate than PB1–216S viruses. By contrast, PB1–216S viruses were more virulent in mice than PB1–216G viruses. These results suggest that the PB1-S216G substitution enhances viral epidemiological fitness by increasing the frequency of adaptive mutations in human influenza A/H1N1 viruses.

Conclusions

Our results thus suggest that the increased adaptability and epidemiological fitness of naturally arising human PB1–216G viruses, which have a canonical low-fidelity replicase, were the biological mechanisms underlying the replacement of PB1–216S viruses with a high-fidelity replicase following the emergence of pdmH1N1. We think that continued surveillance of such naturally occurring PB1–216 variants among others is warranted to assess the potential impact of changes in RdRp fidelity on the adaptability and epidemiological fitness of human A/H1N1 influenza viruses.

Electronic supplementary material

The online version of this article (10.1186/s12929-019-0547-4) contains supplementary material, which is available to authorized users.

Rhinovirus RNA Polymerase

Human rhinovirus is responsible for causing 50% of common cold infections in infants and adults. It belongs to the picornavirus family of nonenveloped positive-strand RNA viruses. The RNA synthesis of rhinovirus is carried out by RNA-dependent RNA polymerase, also known as 3D^Pol. It catalyzes the synthesis of negative-strand RNA using a positive-strand template. The structure of the enzyme consists of three domains: palm, fingers, and thumb domains and Mg²⁺ in the active site. These conserved structural features of the enzyme help in catalyzing phosphodiester bond formation between the two consecutive nucleotide units complimentary to the template RNA using a VPg primer. Owing to the presence of over 100 serotypes of the enzyme, designing specific inhibitors targeting the polymerase is a challenging task and until now no clinically approved antirhino viral drug is reported. In this review, we have given detailed information about the structure and function of the enzyme and also discussed some of the inhibitors and their in vivo activity against 3D^Pol.

A unique intra-molecular fidelity-modulating mechanism identified in a viral RNA-dependent RNA polymerase

Typically not assisted by proofreading, the RNA-dependent RNA polymerases (RdRPs) encoded by the RNA viruses may need to independently control its fidelity to fulfill virus viability and fitness. However, the precise mechanism by which the RdRP maintains its optimal fidelity level remains largely elusive. By solving 2.1-2.5 Å resolution crystal structures of the classical swine fever virus (CSFV) NS5B, an RdRP with a unique naturally fused N-terminal domain (NTD), we identified high-resolution intra-molecular interactions between the NTD and the RdRP palm domain. In order to dissect possible regulatory functions of NTD, we designed mutations at residues Y471 and E472 to perturb key interactions at the NTD-RdRP interface. When crystallized, some of these NS5B interface mutants maintained the interface, while the others adopted an 'open' conformation that no longer retained the intra-molecular interactions. Data from multiple in vitro RdRP assays indicated that the perturbation of the NTD-RdRP interactions clearly reduced the fidelity level of the RNA synthesis, while the processivity of the NS5B elongation complex was not affected. Collectively, our work demonstrates an explicit and unique mode of polymerase fidelity modulation and provides a vivid example of co-evolution in multi-domain enzymes.

RNA Virus Fidelity Mutants: A Useful Tool for Evolutionary Biology or a Complex Challenge?

RNA viruses replicate with low fidelity due to the error-prone nature of the RNA-dependent RNA polymerase, which generates approximately one mutation per round of genome replication. Due to the large population sizes produced by RNA viruses during replication, this results in a cloud of closely related virus variants during host infection, of which small increases or decreases in replication fidelity have been shown to result in virus attenuation in vivo, but not typically in vitro. Since the discovery of the first RNA virus fidelity mutants during the mid-aughts, the field has exploded with the identification of over 50 virus fidelity mutants distributed amongst 7 RNA virus families. This review summarizes the current RNA virus fidelity mutant literature, with a focus upon the definition of a fidelity mutant as well as methods to confirm any mutational changes associated with the fidelity mutant. Due to the complexity of such a definition, in addition to reports of unstable virus fidelity phenotypes, the future translational utility of these mutants and applications for basic science are examined.

Contribution of a Multifunctional Polymerase Region of Foot-and-Mouth Disease Virus to Lethal Mutagenesis

Viral RNA-dependent RNA polymerases (RdRps) are major determinants of high mutation rates and generation of mutant spectra that mediate RNA virus adaptability. The RdRp of the picornavirus foot-and-mouth disease virus (FMDV), termed 3D, is a multifunctional protein that includes a nuclear localization signal (NLS) in its N-terminal region. Previous studies documented that some amino acid substitutions within the NLS altered nucleotide recognition and enhanced the incorporation of the mutagenic purine analogue ribavirin in viral RNA, but the mutants tested were not viable and their response to lethal mutagenesis could not be studied. Here we demonstrate that NLS amino acid substitution M16A of FMDV serotype C does not affect infectious virus production but accelerates ribavirin-mediated virus extinction. The mutant 3D displays polymerase activity, RNA binding, and copying processivity that are similar to those of the wild-type enzyme but shows increased ribavirin-triphosphate incorporation. Crystal structures of the mutant 3D in the apo and RNA-bound forms reveal an expansion of the template entry channel due to the replacement of the bulky Met by Ala. This is a major difference with other 3D mutants with altered nucleotide analogue recognition. Remarkably, two distinct loop β9-α11 conformations distinguish 3Ds that exhibit higher or lower ribavirin incorporation than the wild-type enzyme. This difference identifies a specific molecular determinant of ribavirin sensitivity of FMDV. Comparison of several polymerase mutants indicates that different domains of the molecule can modify nucleotide recognition and response to lethal mutagenesis. The connection of this observation with current views on quasispecies adaptability is discussed.IMPORTANCE The nuclear localization signal (NLS) of the foot-and-mouth disease virus (FMDV) polymerase includes residues that modulate the sensitivity to mutagenic agents. Here we have described a viable NLS mutant with an amino acid replacement that facilitates virus extinction by ribavirin. The corresponding polymerase shows increased incorporation of ribavirin triphosphate and local structural modifications that implicate the template entry channel. Specifically, comparison of the structures of ribavirin-sensitive and ribavirin-resistant FMDV polymerases has identified loop β9-α11 conformation as a determinant of sensitivity to ribavirin mutagenesis.

Extinction of Zika Virus and Usutu Virus by Lethal Mutagenesis Reveals Different Patterns of Sensitivity to Three Mutagenic Drugs

Flaviviruses constitute an increasing source of public health concern, with growing numbers of pathogens causing disease and geographic spread to temperate climates. Despite a large body of evidence supporting mutagenesis as a conceivable antiviral strategy, there are currently no data on the sensitivity to increased mutagenesis for Zika virus (ZIKV) and Usutu virus (USUV), two emerging flaviviral threats.

Flaviviruses constitute an increasing source of public health concern, with growing numbers of pathogens causing disease and geographic spread to temperate climates. Despite a large body of evidence supporting mutagenesis as a conceivable antiviral strategy, there are currently no data on the sensitivity to increased mutagenesis for Zika virus (ZIKV) and Usutu virus (USUV), two emerging flaviviral threats. In this study, we demonstrate that both viruses are sensitive to three ribonucleosides, favipiravir, ribavirin, and 5-fluorouracil, that have shown mutagenic activity against other RNA viruses while remaining unaffected by a mutagenic deoxyribonucleoside. Serial cell culture passages of ZIKV in the presence of these compounds resulted in the rapid extinction of infectivity, suggesting elevated sensitivity to mutagenesis. USUV extinction was achieved when a 10-fold dilution was applied between every passage, but not in experiments involving undiluted virus, indicating an overall lower susceptibility than ZIKV. Although the two viruses are inhibited by the same three drugs, ZIKV is relatively more susceptive to serial passage in the presence of purine analogues (favipiravir and ribavirin), while USUV replication is suppressed more efficiently by 5-fluorouracil. These differences in sensitivity typically correlate with the increases in the mutation frequencies observed in each nucleoside treatment. These results are relevant to the development of efficient therapies based on lethal mutagenesis and support the rational selection of different mutagenic nucleosides for each pathogen. We will discuss the implications of these results to the fidelity of flavivirus replication and the design of antiviral therapies based on lethal mutagenesis.

Mechanisms and Concepts in RNA Virus Population Dynamics and Evolution

RNA viruses are unique in their evolutionary capacity, exhibiting high mutation rates and frequent recombination. They rapidly adapt to environmental changes, such as shifts in immune pressure or pharmacological challenge. The evolution of RNA viruses has been brought into new focus with the recent developments of genetic and experimental tools to explore and manipulate the evolutionary dynamics of viral populations. These studies have uncovered new mechanisms that enable viruses to overcome evolutionary challenges in the environment and have emphasized the intimate relationship of viral populations with evolution. Here, we review some of the emerging viral and host mechanisms that underlie the evolution of RNA viruses. We also discuss new studies that demonstrate that the relationship between evolutionary dynamics and virus biology spans many spatial and temporal scales, affecting transmission dynamics within and between hosts as well as pathogenesis.

A speed-fidelity trade-off determines the mutation rate and virulence of an RNA virus

Mutation rates can evolve through genetic drift, indirect selection due to genetic hitchhiking, or direct selection on the physicochemical cost of high fidelity. However, for many systems, it has been difficult to disentangle the relative impact of these forces empirically. In RNA viruses, an observed correlation between mutation rate and virulence has led many to argue that their extremely high mutation rates are advantageous because they may allow for increased adaptability. This argument has profound implications because it suggests that pathogenesis in many viral infections depends on rare or de novo mutations. Here, we present data for an alternative model whereby RNA viruses evolve high mutation rates as a byproduct of selection for increased replicative speed. We find that a poliovirus antimutator, 3D^G64S, has a significant replication defect and that wild-type (WT) and 3D^G64S populations have similar adaptability in 2 distinct cellular environments. Experimental evolution of 3D^G64S under selection for replicative speed led to reversion and compensation of the fidelity phenotype. Mice infected with 3D^G64S exhibited delayed morbidity at doses well above the lethal level, consistent with attenuation by slower growth as opposed to reduced mutational supply. Furthermore, compensation of the 3D^G64S growth defect restored virulence, while compensation of the fidelity phenotype did not. Our data are consistent with the kinetic proofreading model for biosynthetic reactions and suggest that speed is more important than accuracy. In contrast with what has been suggested for many RNA viruses, we find that within-host spread is associated with viral replicative speed and not standing genetic diversity.

Why organisms have different mutation rates is a longstanding question in evolutionary biology. The polymerases of RNA viruses generally lack proofreading activity and exhibit extremely high mutation rates. Because most mutations are deleterious and mutation rates are typically tuned by natural selection, we asked why RNA viruses haven’t evolved a polymerase with a lower mutation rate. We used experimental evolution and a murine infection model to show that RNA virus mutation rates may actually be too high and are not necessarily adaptive. Rather, our data indicate that viral mutation rates have evolved to be higher as a result of selection for viruses with faster replication kinetics. We suggest that viruses have high mutation rates, not because they facilitate adaptation but because it is hard to be both fast and accurate and these viruses have prioritized speed over fidelity.

Emergency Services of Viral RNAs: Repair and Remodeling

Reproduction of RNA viruses is typically error-prone due to the infidelity of their replicative machinery and the usual lack of proofreading mechanisms. The error rates may be close to those that kill the virus. Consequently, populations of RNA viruses are represented by heterogeneous sets of genomes with various levels of fitness. This is especially consequential when viruses encounter various bottlenecks and new infections are initiated by a single or few deviating genomes. Nevertheless, RNA viruses are able to maintain their identity by conservation of major functional elements. This conservatism stems from genetic robustness or mutational tolerance, which is largely due to the functional degeneracy of many protein and RNA elements as well as to negative selection. Another relevant mechanism is the capacity to restore fitness after genetic damages, also based on replicative infidelity. Conversely, error-prone replication is a major tool that ensures viral evolvability. The potential for changes in debilitated genomes is much higher in small populations, because in the absence of stronger competitors low-fit genomes have a choice of various trajectories to wander along fitness landscapes. Thus, low-fit populations are inherently unstable, and it may be said that to run ahead it is useful to stumble. In this report, focusing on picornaviruses and also considering data from other RNA viruses, we review the biological relevance and mechanisms of various alterations of viral RNA genomes as well as pathways and mechanisms of rehabilitation after loss of fitness. The relationships among mutational robustness, resilience, and evolvability of viral RNA genomes are discussed.

Low-fidelity Venezuelan equine encephalitis virus polymerase mutants to improve live-attenuated vaccine safety and efficacy

During RNA virus replication, there is the potential to incorporate mutations that affect virulence or pathogenesis. For live-attenuated vaccines, this has implications for stability, as replication may result in mutations that either restore the wild-type phenotype via reversion or compensate for the attenuating mutations by increasing virulence (pseudoreversion). Recent studies have demonstrated that altering the mutation rate of an RNA virus is an effective attenuation tool. To validate the safety of low-fidelity mutations to increase vaccine attenuation, several mutations in the RNA-dependent RNA-polymerase (RdRp) were tested in the live-attenuated Venezuelan equine encephalitis virus vaccine strain, TC-83. Next generation sequencing after passage in the presence of mutagens revealed a mutant containing three mutations in the RdRp, TC-83 3x, to have decreased replication fidelity, while a second mutant, TC-83 4x displayed no change in fidelity, but shared many phenotypic characteristics with TC-83 3x. Both mutants exhibited increased, albeit inconsistent attenuation in an infant mouse model, as well as increased immunogenicity and complete protection against lethal challenge of an adult murine model compared with the parent TC-83. During serial passaging in a highly permissive model, the mutants increased in virulence but remained less virulent than the parent TC-83. These results suggest that the incorporation of low-fidelity mutations into the RdRp of live-attenuated vaccines for RNA viruses can confer increased immunogenicity whilst showing some evidence of increased attenuation. However, while in theory such constructs may result in more effective vaccines, the instability of the vaccine phenotype decreases the likelihood of this being an effective vaccine strategy.

Foot-and-mouth disease virus type O specific mutations determine RNA-dependent RNA polymerase fidelity and virus attenuation

Previous studies have shown that the FMDV Asia1/YS/CHA/05 high-fidelity mutagen-resistant variants are attenuated (Zeng et al., 2014). Here, we introduced the same single or multiple-amino-acid substitutions responsible for increased 3D^pol fidelity of type Asia1 FMDV into the type O FMDV O/YS/CHA/05 infectious clone. The rescued viruses O-DA and O-DAMM are lower replication fidelity mutants and showed an attenuated phenotype. These results demonstrated that the same amino acid substitution of 3D^pol in different serotypes of FMDV strains had different effects on viral fidelity. In addition, nucleoside analogues were used to select high-fidelity mutagen-resistant type O FMDV variants. The rescued mutagen-resistant type O FMDV high-fidelity variants exhibited significantly attenuated fitness and a reduced virulence phenotype. These results have important implications for understanding the molecular mechanism of FMDV evolution and pathogenicity, especially in developing a safer modified live-attenuated vaccine against FMDV.