Effective lethal mutagenesis of influenza virus by three nucleoside analogs

Ribavirin inhibits the replication of infectious bursal disease virus predominantly through depletion of cellular guanosine pool

Introduction

The antiviral activity of different mutagens against single-stranded RNA viruses is well documented; however, their activity on the replication of double-stranded RNA viruses remains unexplored. This study aims to investigate the effect of different antivirals on the replication of a chicken embryo fibroblast-adapted Infectious Bursal Disease virus, FVSKG2. This study further explores the antiviral mechanism utilized by the most effective anti-IBDV agent.

Methods

The cytotoxicity and anti-FVSKG2 activity of different antiviral agents (ribavirin, 5-fluorouracil, 5-azacytidine, and amiloride) were evaluated. The virus was serially passaged in chicken embryo fibroblasts 11 times at sub-cytotoxic concentrations of ribavirin, 5-fluorouracil or amiloride. Further, the possible mutagenic and non-mutagenic mechanisms utilized by the most effective anti-FVSKG2 agent were explored.

Results and Discussion

Ribavirin was the least cytotoxic on chicken embryo fibroblasts, followed by 5-fluorouracil, amiloride and 5-azacytidine. Ribavirin inhibited the replication of FVSKG2 in chicken embryo fibroblasts significantly at concentrations as low as 0.05 mM. The extinction of FVSKG2 was achieved during serial passage of the virus in chicken embryo fibroblasts at ≥0.05 mM ribavirin; however, the emergence of a mutagen-resistant virus was not observed until the eleventh passage. Further, no mutation was observed in 1,898 nucleotides of the FVSKG2 following its five passages in chicken embryo fibroblasts in the presence of 0.025 mM ribavirin. Ribavarin inhibited the FVSKG2 replication in chicken embryo fibroblasts primarily through IMPDH-mediated depletion of the Guanosine Triphosphate pool of cells. However, other mechanisms like ribavirin-mediated cytokine induction or possible inhibition of viral RNA-dependent RNA polymerase through its interaction with the enzyme's active sites enhance the anti-IBDV effect. Ribavirin inhibits ds- RNA viruses, likely through IMPDH inhibition and not mutagenesis. The inhibitory effect may, however, be augmented by other non-mutagenic mechanisms, like induction of antiviral cytokines in chicken embryo fibroblasts or interaction of ribavirin with the active sites of RNA-dependent RNA polymerase of the virus.

The combined effect of oseltamivir and favipiravir on influenza a virus evolution in patients hospitalized with severe influenza

Our previous study shows favipiravir and oseltamivir combination therapy may accelerate clinical recovery compared to oseltamivir monotherapy in severe influenza, but its effect on virological evolution and resistance mutation against oseltamivir is still unknown. In this study, we collected longitudinal respiratory samples from influenza patients who underwent combination therapy and applied them to next generation sequencing of the whole genome of the influenza A virus (IAV). We also included a cohort untreated with any antivirals to serve as the control. In total, 62 samples from 19 patients treated with combination therapy and 20 samples from 20 patients untreated were successfully sequenced. The nucleotide diversity in the whole genome of IAV in the combination group showed no difference compared to that in the control group (P > 0.05). Moreover, we observed 174 kinds of nonsynonymous nucleotide substitutions in patients with combination therapy, mostly in NA (n = 44) and HA (n = 43). Of them, the G→A transition was the dominant variant type (27%) and 46/174 (26%) was reported to have biological effects, such as increased pathogenicity and polymerase activity. Among the 29 mutations conferring reduction in oseltamivir sensitivity we investigated, H275Y was the only mutation detected in the 4 samples from 1 of 19 patients and demonstrated increasing frequency during the treatment. Mutations conferring favipiravir resistance were not observed. Our studies showed combination therapy of favipiravir and oseltamivir has little effect on virus nucleotide diversity, nor prevents the increase of oseltamivir-resistant variants.

Heterogeneity in viral populations increases the rate of deleterious mutation accumulation

RNA viruses have high mutation rates, with the majority of mutations being deleterious. We examine patterns of deleterious mutation accumulation over multiple rounds of viral replication, with a focus on how cellular coinfection and heterogeneity in viral output affect these patterns. Specifically, using agent-based intercellular simulations we find, in agreement with previous studies, that coinfection of cells by viruses relaxes the strength of purifying selection, and thereby increases the rate of deleterious mutation accumulation. We further find that cellular heterogeneity in viral output exacerbates the rate of deleterious mutation accumulation, regardless of whether this heterogeneity in viral output is stochastic or is due to variation in cellular multiplicity of infection. These results highlight the need to consider the unique life histories of viruses and their population structure to better understand observed patterns of viral evolution.

Screening of SARS-CoV-2 antivirals through a cell-based RNA-dependent RNA polymerase (RdRp) reporter assay

COVID-19 (Coronavirus Disease 2019) caused by SARS-CoV-2 (Severe Acute Respiratory Syndrome CoronaVirus-2) continues to pose an international public health threat and thus far, has resulted in greater than 6.4 million deaths worldwide. Vaccines are critical tools to limit COVID-19 spread, but antiviral drug development is an ongoing global priority due to fast-spreading COVID-19 variants that may elude vaccine efficacies. The RNA-dependent RNA polymerase (RdRp) of SARS-CoV-2 is an essential enzyme of viral replication and transcription machinery complex. Therefore, the RdRp is an attractive target for the development of effective anti-COVID-19 therapeutics. In this study, we developed a cell-based assay to determine the enzymatic activity of SARS-CoV-2 RdRp through a luciferase reporter system. The SARS-CoV-2 RdRp reporter assay was validated using known inhibitors of RdRp polymerase, remdesivir along with other anti-virals including ribavirin, penciclovir, rhoifolin, 5′CT, and dasabuvir. Dasabuvir (an FDA-approved drug) exhibited promising RdRp inhibitory activity among these inhibitors. Anti-viral activity of dasabuvir was also tested on the replication of SARS-CoV-2 through infection of Vero E6 cells. Dasabuvir inhibited the replication of SARS-CoV-2, USA-WA1/2020 as well as B.1.617.2 (delta variant) in Vero E6 cells in a dose-dependent manner with EC₅₀ values 9.47 μM and 10.48 μM, for USA-WA1/2020 and B.1.617.2 variants, respectively. Our results suggest that dasabuvir can be further evaluated as a therapeutic drug for COVID-19. Importantly, this system provides a robust, target-specific, and high-throughput screening compatible (z- and z’-factors of >0.5) platforms that will be a valuable tool for screening SARS-CoV-2 RdRp inhibitors.

Potential and action mechanism of favipiravir as an antiviral against Junin virus

Favipiravir is a nucleoside analogue that inhibits the replication and transcription of a broad spectrum of RNA viruses, including pathogenic arenaviruses. In this study, we isolated a favipiravir-resistant mutant of Junin virus (JUNV), which is the causative agent of Argentine hemorrhagic fever, and analyzed the antiviral mechanism of favipiravir against JUNV. Two amino acid substitutions, N462D in the RNA-dependent RNA polymerase (RdRp) and A168T in the glycoprotein precursor GPC, were identified in the mutant. GPC-A168T substitution enhanced the efficiency of JUNV internalization, which explains the robust replication kinetics of the mutant in the virus growth analysis. Although RdRp-N462D substitution did not affect polymerase activity levels in a minigenome system, comparisons of RdRp error frequencies showed that the virus with RdRp-D462 possessed a significantly higher fidelity. Our next generation sequence (NGS) analysis showed a gradual accumulation of both mutations as we passaged the virus in presence of favipiravir. We also provided experimental evidence for the first time that favipiravir inhibited JUNV through the accumulation of transition mutations, confirming its role as a purine analogue against arenaviruses. Moreover, we showed that treatment with a combination of favipiravir and either ribavirin or remdesivir inhibited JUNV replication in a synergistic manner, blocking the generation of the drug-resistant mutant. Our findings provide new insights for the clinical management and treatment of Argentine hemorrhagic fever.

Development of antivirals requires cautious and extensive assessment of action mechanism as well as potential for emergence of resistant phenotype of the virus. In recent years, favipiravir has been put forward as a promising candidate for the treatment of Argentine hemorrhagic fever (AHF) caused by Junin virus (JUNV). We, therefore, aimed to provide experimental evidence on action mechanism of favipiravir to help guide its clinical use. Here we show that favipiravir causes lethal mutation that impairs virus infectivity. More importantly, we demonstrate that the virus has the capability to escape favipiravir selective pressure by acquiring two amino acid substitutions on glycoprotein precursor and polymerase proteins. This observation raises concern over the use of only favipiravir in therapeutic regimens. To overcome this risk, we show that combination of favipiravir with other nucleoside analogues demonstrates a synergistic effect and suppresses the ability of JUNV to escape drug pressure. Favipiravir, ribavirin, and remdesivir have a broad spectrum of antiviral activity. Therefore, combination therapies of these drugs would be expected to have potential therapeutic effects for not only AHF but also the diseases caused by a variety of viruses, including emerging RNA viruses.

A Virus Is a Community: Diversity within Negative-Sense RNA Virus Populations

Negative-sense RNA virus populations are composed of diverse viral components that interact to form a community and shape the outcome of virus infections. At the genomic level, RNA virus populations consist not only of a homogeneous population of standard viral genomes but also of an extremely large number of genome variants, termed viral quasispecies, and nonstandard viral genomes, which include copy-back viral genomes, deletion viral genomes, mini viral RNAs, and hypermutated RNAs. At the particle level, RNA virus populations are composed of pleomorphic particles, particles missing or having additional genomes, and single particles or particle aggregates. As we continue discovering more about the components of negative-sense RNA virus populations and their crucial functions during virus infection, it will become more important to study RNA virus populations as a whole rather than their individual parts. In this review, we will discuss what is known about the components of negative-sense RNA virus communities, speculate how the components of the virus community interact, and summarize what vaccines and antiviral therapies are being currently developed to target or harness these components.

Lethal Mutagenesis of RNA Viruses and Approved Drugs with Antiviral Mutagenic Activity

In RNA viruses, a small increase in their mutation rates can be sufficient to exceed their threshold of viability. Lethal mutagenesis is a therapeutic strategy based on the use of mutagens, driving viral populations to extinction. Extinction catastrophe can be experimentally induced by promutagenic nucleosides in cell culture models. The loss of HIV infectivity has been observed after passage in 5-hydroxydeoxycytidine or 5,6-dihydro-5-aza-2′-deoxycytidine while producing a two-fold increase in the viral mutation frequency. Among approved nucleoside analogs, experiments with polioviruses and other RNA viruses suggested that ribavirin can be mutagenic, although its mechanism of action is not clear. Favipiravir and molnupiravir exert an antiviral effect through lethal mutagenesis. Both drugs are broad-spectrum antiviral agents active against RNA viruses. Favipiravir incorporates into viral RNA, affecting the G→A and C→U transition rates. Molnupiravir (a prodrug of β-d-N⁴-hydroxycytidine) has been recently approved for the treatment of SARS-CoV-2 infection. Its triphosphate derivative can be incorporated into viral RNA and extended by the coronavirus RNA polymerase. Incorrect base pairing and inefficient extension by the polymerase promote mutagenesis by increasing the G→A and C→U transition frequencies. Despite having remarkable antiviral action and resilience to drug resistance, carcinogenic risks and genotoxicity are important concerns limiting their extended use in antiviral therapy.

Low Pathogenicity H7N3 Avian Influenza Viruses Have Higher Within-Host Genetic Diversity Than a Closely Related High Pathogenicity H7N3 Virus in Infected Turkeys and Chickens

Within-host viral diversity offers a view into the early stages of viral evolution occurring after a virus infects a host. In recent years, advances in deep sequencing have allowed for routine identification of low-frequency variants, which are important sources of viral genetic diversity and can potentially emerge as a major virus population under certain conditions. We examined within-host viral diversity in turkeys and chickens experimentally infected with closely related H7N3 avian influenza viruses (AIVs), specifically one high pathogenicity AIV (HPAIV) and two low pathogenicity AIV (LPAIVs) with different neuraminidase protein stalk lengths. Consistent with the high mutation rates of AIVs, an abundance of intra-host single nucleotide variants (iSNVs) at low frequencies of 2–10% was observed in all samples collected. Furthermore, a small number of common iSNVs were observed between turkeys and chickens, and between directly inoculated and contact-exposed birds. Notably, the LPAIVs have significantly higher iSNV diversities and frequencies of nonsynonymous changes than the HPAIV in both turkeys and chickens. These findings highlight the dynamics of AIV populations within hosts and the potential impact of genetic changes, including mutations in the hemagglutinin gene that confers the high pathogenicity pathotype, on AIV virus populations and evolution.

Favipiravir Inhibits Hepatitis A Virus Infection in Human Hepatocytes

Hepatitis A virus (HAV) is a causative agent of acute hepatitis and can occasionally induce acute liver failure. However, specific potent anti-HAV drug is not available on the market currently. Thus, we investigated several novel therapeutic drugs through a drug repositioning approach, targeting ribonucleic acid (RNA)-dependent RNA polymerase and RNA-dependent deoxyribonucleic acid polymerase. In the present study, we examined the anti-HAV activity of 18 drugs by measuring the HAV subgenomic replicon and HAV HA11-1299 genotype IIIA replication in human hepatoma cell lines, using a reporter assay and real-time reverse transcription polymerase chain reaction, respectively. Mutagenesis of the HAV 5’ untranslated region was also examined by next-generation sequencing. These specific parameters were explored because lethal mutagenesis has emerged as a novel potential therapeutic approach to treat RNA virus infections. Favipiravir inhibited HAV replication in both Huh7 and PLC/PRF/5 cells, although ribavirin inhibited HAV replication in only Huh7 cells. Next-generation sequencing demonstrated that favipiravir could introduce nucleotide mutations into the HAV genome more than ribavirin. In conclusion, favipiravir could introduce nucleotide mutations into the HAV genome and work as an antiviral against HAV infection. Provided that further in vivo experiments confirm its efficacy, favipiravir would be useful for the treatment of severe HAV infection.

Peering into Avian Influenza A(H5N8) for a Framework towards Pandemic Preparedness

2014 marked the first emergence of avian influenza A(H5N8) in Jeonbuk Province, South Korea, which then quickly spread worldwide. In the midst of the 2020-2021 H5N8 outbreak, it spread to domestic poultry and wild waterfowl shorebirds, leading to the first human infection in Astrakhan Oblast, Russia. Despite being clinically asymptomatic and without direct human-to-human transmission, the World Health Organization stressed the need for continued risk assessment given the nature of Influenza to reassort and generate novel strains. Given its promiscuity and easy cross to humans, the urgency to understand the mechanisms of possible species jumping to avert disastrous pandemics is increasing. Addressing the epidemiology of H5N8, its mechanisms of species jumping and its implications, mutational and reassortment libraries can potentially be built, allowing them to be tested on various models complemented with deep-sequencing and automation. With knowledge on mutational patterns, cellular pathways, drug resistance mechanisms and effects of host proteins, we can be better prepared against H5N8 and other influenza A viruses.

Tradition to Pathogenesis: A Novel Hypothesis for Elucidating the Pathogenesis of Diseases Based on the Traditional Use of Medicinal Plants

Traditional medicines embody knowledge on medicinal plants that has been accumulated through cultural evolution over millennia. In the latter half of the 20th century, two approaches to medicinal plant research have been established: the "Bench to Bedside" and the "Bedside to Bench" approaches which serve primarily for the development of more efficient therapeutics. Here, we propose a third, novel approach: from "Tradition to Pathogenesis" which aims to understand the pathogenesis of diseases based on the cultural evolution of their respective empirical treatments. We analyse multiple examples of diseases where the acting mechanism of traditional treatments across multiple cultures points to the pathogenesis of the respective disease. E.g., many cultures traditionally treat rheumatism with anti-bacterial botanical drugs, which is at odds with our current understanding that rheumatism is an aseptic inflammation. Furthermore, gastric ailments have traditionally been treated with anti-infectious botanical drugs indicating local infections, as demonstrated by the discovery of Helicobacter pylori as a common cause of gastric ulcer. Understanding traditional treatments can thus help to elucidate the pathogenesis of the disease.

Quaternary ammonium salts based on (-)-borneol as effective inhibitors of influenza virus

A series of compounds containing a 1,7,7-trimethylbicyclo[2.2.1]heptane fragment were evaluated for their antiviral activity against influenza A virus strain A/Puerto Rico/8/34 (H1N1) in vitro. The most potent antiviral compound proved to be a quaternary ammonium salt based on (-)-borneol, 10a. In in vitro experiments, compound 10a inhibited influenza A viruses (H1, H1pdm09, and H3 subtypes), with an IC₅₀ value of 2.4-16.8 µM (depending on the virus), and demonstrated low toxicity (CC₅₀ = 1311 µM). Mechanism-of-action studies for compound 10a revealed it to be most effective when added at the early stages of the viral life cycle. In direct haemolysis inhibition tests, compound 10a was shown to decrease the membrane-disrupting activity of influenza A virus strain A/Puerto Rico/8/34. According to molecular modelling results, the lead compound 10a can bind to different sites in the stem region of the viral hemagglutinin.

SARS-CoV-2 biology and variants: anticipation of viral evolution and what needs to be done

The global propagation of SARS-CoV-2 and the detection of a large number of variants, some of which have replaced the original clade to become dominant, underscores the fact that the virus is actively exploring its evolutionary space. The longer high levels of viral multiplication occur - permitted by high levels of transmission -, the more the virus can adapt to the human host and find ways to success. The third wave of the COVID-19 pandemic is starting in different parts of the world, emphasizing that transmission containment measures that are being imposed are not adequate. Part of the consideration in determining containment measures is the rationale that vaccination will soon stop transmission and allow a return to normality. However, vaccines themselves represent a selection pressure for evolution of vaccine-resistant variants, so the coupling of a policy of permitting high levels of transmission/virus multiplication during vaccine roll-out with the expectation that vaccines will deal with the pandemic, is unrealistic. In the absence of effective antivirals, it is not improbable that SARS-CoV-2 infection prophylaxis will involve an annual vaccination campaign against 'dominant' viral variants, similar to influenza prophylaxis. Living with COVID-19 will be an issue of SARS-CoV-2 variants and evolution. It is therefore crucial to understand how SARS-CoV-2 evolves and what constrains its evolution, in order to anticipate the variants that will emerge. Thus far, the focus has been on the receptor-binding spike protein, but the virus is complex, encoding 26 proteins which interact with a large number of host factors, so the possibilities for evolution are manifold and not predictable a priori. However, if we are to mount the best defence against COVID-19, we must mount it against the variants, and to do this, we must have knowledge about the evolutionary possibilities of the virus. In addition to the generic cellular interactions of the virus, there are extensive polymorphisms in humans (e.g. Lewis, HLA, etc.), some distributed within most or all populations, some restricted to specific ethnic populations and these variations pose additional opportunities for/constraints on viral evolution. We now have the wherewithal - viral genome sequencing, protein structure determination/modelling, protein interaction analysis - to functionally characterize viral variants, but access to comprehensive genome data is extremely uneven. Yet, to develop an understanding of the impacts of such evolution on transmission and disease, we must link it to transmission (viral epidemiology) and disease data (patient clinical data), and the population granularities of these. In this editorial, we explore key facets of viral biology and the influence of relevant aspects of human polymorphisms, human behaviour, geography and climate and, based on this, derive a series of recommendations to monitor viral evolution and predict the types of variants that are likely to arise.

The Roles of Host 5-Methylcytosine RNA Methyltransferases during Viral Infections

Eukaryotic 5-methylcytosine RNA methyltransferases catalyze the transfer of a methyl group to the fifth carbon of a cytosine base in RNA sequences to produce 5-methylcytosine (m5C). m5C RNA methyltransferases play a crucial role in the maintenance of functionality and stability of RNA. Viruses have developed a number of strategies to suppress host innate immunity and ensure efficient transcription and translation for the replication of new virions. One such viral strategy is to use host m5C RNA methyltransferases to modify viral RNA and thus to affect antiviral host responses. Here, we summarize the latest findings concerning the roles of m5C RNA methyltransferases, namely, NOL1/NOP2/SUN domain (NSUN) proteins and DNA methyltransferase 2/tRNA methyltransferase 1 (DNMT2/TRDMT1) during viral infections. Moreover, the use of m5C RNA methyltransferase inhibitors as an antiviral therapy is discussed.

Synthesis of 1,2,3-triazolyl nucleoside analogues and their antiviral activity

Abstract

Based on the fact that a search for influenza antivirals among nucleoside analogues has drawn very little attention of chemists, the present study reports the synthesis of a series of 1,2,3-triazolyl nucleoside analogues in which a pyrimidine fragment is attached to the ribofuranosyl-1,2,3-triazol-4-yl moiety by a polymethylene linker of variable length. Target compounds were prepared by the Cu alkyne-azide cycloaddition (CuAAC) reaction. Derivatives of uracil, 6-methyluracil, 3,6-dimethyluracil, thymine and quinazolin-2,4-dione with ω-alkyne substituent at the N1 (or N5) atom and azido 2,3,5-tri-O-acetyl-D-β-ribofuranoside were used as components of the CuAAC reaction. All compounds synthesized were evaluated for antiviral activity against influenza virus A/PR/8/34/(H1N1) and coxsackievirus B3. The best values of IC₅₀ (inhibiting concentration) and SI (selectivity index) were demonstrated by the lead compound 4i in which the 1,2,3-triazolylribofuranosyl fragment is attached to the N1 atom of the quinazoline-2,4-dione moiety via a butylene linker (IC₅₀ = 30 μM, SI = 24) and compound 8n in which the 1,2,3-triazolylribofuranosyl fragment is attached directly to the N5 atom of the 6-methyluracil moiety (IC₅₀ = 15 μM, SI = 5). According to theoretical calculations, the antiviral activity of the 1,2,3-triazolyl nucleoside analogues 4i and 8n against H1N1 (A/PR/8/34) influenza virus can be explained by their influence on the functioning of the polymerase acidic protein (PA) of RNA-dependent RNA polymerase (RdRP).

Graphic abstract

Electronic supplementary material

The online version of this article (10.1007/s11030-020-10141-y) contains supplementary material, which is available to authorized users.

Selected nucleos(t)ide-based prescribed drugs and their multi-target activity

Nucleos(t)ide analogues play pivotal roles as antiviral, cytotoxic or immunosuppressive agents. Here, we review recent reports of nucleoside analogues that exhibit broad-spectrum activity towards multiple life-threatening RNA and DNA viruses. We also present a discussion about nucleoside antimetabolites-approved antineoplastic agents-that have recently been shown to have antiviral and/or antibacterial activity. The approved drugs and drug combinations, as well as recently identified candidates for investigation and/or experimentation, are discussed. Several examples of repurposed drugs that have already been approved for use are presented. This strategy can be crucial for the first-line treatment of acute infections or coinfections and for the management of drug-resistant strains.

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.

Harnessing host-virus evolution in antiviral therapy and immunotherapy

Pathogen resistance and development costs are major challenges in current approaches to antiviral therapy. The high error rate of RNA synthesis and reverse‐transcription confers genome plasticity, enabling the remarkable adaptability of RNA viruses to antiviral intervention. However, this property is coupled to fundamental constraints including limits on the size of information available to manipulate complex hosts into supporting viral replication. Accordingly, RNA viruses employ various means to extract maximum utility from their informationally limited genomes that, correspondingly, may be leveraged for effective host‐oriented therapies. Host‐oriented approaches are becoming increasingly feasible because of increased availability of bioactive compounds and recent advances in immunotherapy and precision medicine, particularly genome editing, targeted delivery methods and RNAi. In turn, one driving force behind these innovations is the increasingly detailed understanding of evolutionarily diverse host–virus interactions, which is the key concern of an emerging field, neo‐virology. This review examines biotechnological solutions to disease and other sustainability issues of our time that leverage the properties of RNA and DNA viruses as developed through co‐evolution with their hosts.

Existing Host Range Mutations Constrain Further Emergence of RNA Viruses

RNA viruses are capable of rapid host shifting, typically due to a point mutation that confers expanded host range. As additional point mutations are necessary for further expansions, epistasis among host range mutations can potentially affect the mutational neighborhood and frequency of niche expansion. We mapped the mutational neighborhood of host range expansion using three genotypes of the double-stranded RNA (dsRNA) bacteriophage φ6 (wild type and two isogenic host range mutants) on the novel host Pseudomonas syringae pv. atrofaciens. Both Sanger sequencing of 50 P. syringae pv. atrofaciens mutant clones for each genotype and population Illumina sequencing revealed the same high-frequency mutations allowing infection of P. syringae pv. atrofaciens. Wild-type φ6 had at least nine different ways of mutating to enter the novel host, eight of which are in p3 (host attachment protein gene), and 13/50 clones had unchanged p3 genes. However, the two isogenic mutants had dramatically restricted neighborhoods: only one or two mutations, all in p3. Deep sequencing revealed that wild-type clones without mutations in p3 likely had changes in p12 (morphogenic protein), a region that was not polymorphic for the two isogenic host range mutants. Sanger sequencing confirmed that 10/13 of the wild-type φ6 clones had nonsynonymous mutations in p12, and 2 others had point mutations in p9 and p5. None of these genes had previously been associated with host range expansion in φ6. We demonstrate, for the first time, epistatic constraint in an RNA virus due to host range mutations themselves, which has implications for models of serial host range expansion.IMPORTANCE RNA viruses mutate rapidly and frequently expand their host ranges to infect novel hosts, leading to serial host shifts. Using an RNA bacteriophage model system (Pseudomonas phage φ6), we studied the impact of preexisting host range mutations on another host range expansion. Results from both clonal Sanger and Illumina sequencing show that extant host range mutations dramatically narrow the neighborhood of potential host range mutations compared to that of wild-type φ6. This research suggests that serial host-shifting viruses may follow a small number of molecular paths to enter additional novel hosts. We also identified new genes involved in φ6 host range expansion, expanding our knowledge of this important model system in experimental evolution.

Determining the Mutation Bias of Favipiravir in Influenza Virus Using Next-Generation Sequencing

New antiviral drugs are needed as a first line of defense in the event of a novel influenza pandemic. Favipiravir is a broad-spectrum antiviral which is effective against influenza. The exact mechanism of how favipiravir works to inhibit influenza is still unclear. We used next-generation sequencing (NGS) to demonstrate that favipiravir causes mutations in influenza RNA. The greater depth of NGS sequence information over traditional sequencing methods allowed us to precisely determine the bias of particular mutations caused by favipiravir. NGS can also be used in a standard diagnostic pipeline to show that favipiravir is acting on the virus by revealing the mutation bias pattern typical to the drug. Our work will aid in testing whether viruses are resistant to favipiravir and may help demonstrate the effect of favipiravir on viruses in a clinical setting. This will be important if favipiravir is used during a future influenza pandemic.

Favipiravir is a broad-spectrum antiviral drug that may be used to treat influenza. Previous research has identified that favipiravir likely acts as a mutagen, but the precise mutation bias that favipiravir induces in influenza virus RNAs has not been described. Here, we use next-generation sequencing (NGS) with barcoding of individual RNA molecules to accurately and quantitatively detect favipiravir-induced mutations and to sample orders of magnitude more mutations than would be possible through Sanger sequencing. We demonstrate that favipiravir causes mutations and show that favipiravir primarily acts as a guanine analogue and secondarily as an adenine analogue resulting in the accumulation of transition mutations. We also use a standard NGS pipeline to show that the mutagenic effect of favipiravir can be measured by whole-genome sequencing of virus.

IMPORTANCE New antiviral drugs are needed as a first line of defense in the event of a novel influenza pandemic. Favipiravir is a broad-spectrum antiviral which is effective against influenza. The exact mechanism of how favipiravir works to inhibit influenza is still unclear. We used next-generation sequencing (NGS) to demonstrate that favipiravir causes mutations in influenza RNA. The greater depth of NGS sequence information over traditional sequencing methods allowed us to precisely determine the bias of particular mutations caused by favipiravir. NGS can also be used in a standard diagnostic pipeline to show that favipiravir is acting on the virus by revealing the mutation bias pattern typical to the drug. Our work will aid in testing whether viruses are resistant to favipiravir and may help demonstrate the effect of favipiravir on viruses in a clinical setting. This will be important if favipiravir is used during a future influenza pandemic.

A Defect in Influenza A Virus Particle Assembly Specific to Primary Human Macrophages

Identification of host cell determinants promoting or suppressing replication of viruses has been aided by analyses of host cells that impose inherent blocks on viral replication. In this study, we show that primary human MDM, which are not permissive to IAV replication, fail to support virus particle formation. This defect is specific to primary human macrophages, since a human monocytic cell line differentiated to macrophage-like cells supports IAV particle formation. We further identified association between two viral transmembrane proteins, HA and M2, on the cell surface as a discrete assembly step, which is defective in MDM. Defective HA-M2 association and particle budding, but not virus release, in MDM are rescued by disruption of actin cytoskeleton, revealing a previously unknown, negative role for actin, which specifically targets an early step in the multistep IAV production. Overall, our study uncovered a host-mediated restriction of association between viral transmembrane components during IAV assembly.

Influenza A virus (IAV) propagates efficiently in epithelial cells, its primary target in the respiratory tract. In contrast, productive infection of most IAV strains is either blocked or highly inefficient in macrophages. The exact nature of the defect in IAV replication in human macrophages remains unknown. In this study, we showed that even compared to a monocytic cell line differentiated to macrophage-like cells, primary human monocyte-derived macrophages (MDM) are inefficient in IAV production, despite comparable levels of expression of viral glycoproteins at the plasma membrane. Correlative fluorescence scanning electron microscopy revealed that formation of budding structures at the cell surface is inefficient in MDM even though clustering of a viral glycoprotein, hemagglutinin (HA), is observed, suggesting that a step in IAV particle assembly is blocked in MDM. Using an in situ proximity ligation assay, we further determined that HA associates with neuraminidase (NA) but fails to associate with another viral transmembrane protein, M2, at the MDM plasma membrane. Notably, the defects in HA-M2 association and particle assembly in MDM were reversed upon cytochalasin D treatment that inhibits actin polymerization. These results suggest that HA-M2 association on the plasma membrane is a discrete step in IAV production, which is susceptible to suppression by actin cytoskeleton in MDM. Virus release remained inefficient in MDM upon cytochalasin D treatment, suggesting the presence of an additional defect(s) in virus release in this cell type. Overall, our study revealed the presence of multiple cell-type-specific mechanisms negatively regulating IAV production at the plasma membrane in MDM.

The mechanism of resistance to favipiravir in influenza

Favipiravir is a broad-spectrum antiviral that has shown promise in treatment of influenza virus infections, in particular due to the apparent lack of emergence of resistance mutations against the drug in cell culture or animal studies. We demonstrate here that a mutation in a conserved region of the viral RNA polymerase confers resistance to favipiravir in vitro and in cell culture. The resistance mutation has a cost to viral fitness, but this can be restored by a compensatory mutation in the polymerase. Our findings support the development of favipiravir-resistance diagnostic and surveillance testing strategies and reinforce the importance of considering combinations of therapies to treat influenza infections.

Favipiravir is a broad-spectrum antiviral that has shown promise in treatment of influenza virus infections. While emergence of resistance has been observed for many antiinfluenza drugs, to date, clinical trials and laboratory studies of favipiravir have not yielded resistant viruses. Here we show evolution of resistance to favipiravir in the pandemic H1N1 influenza A virus in a laboratory setting. We found that two mutations were required for robust resistance to favipiravir. We demonstrate that a K229R mutation in motif F of the PB1 subunit of the influenza virus RNA-dependent RNA polymerase (RdRP) confers resistance to favipiravir in vitro and in cell culture. This mutation has a cost to viral fitness, but fitness can be restored by a P653L mutation in the PA subunit of the polymerase. K229R also conferred favipiravir resistance to RNA polymerases of other influenza A virus strains, and its location within a highly conserved structural feature of the RdRP suggests that other RNA viruses might also acquire resistance through mutations in motif F. The mutations identified here could be used to screen influenza virus-infected patients treated with favipiravir for the emergence of resistance.

Viral replication modes in single-peak fitness landscapes: A dynamical systems analysis

Positive-sense, single-stranded RNA viruses are important pathogens infecting almost all types of organisms. Experimental evidence from distributions of mutations and from viral RNA amplification suggest that these pathogens may follow different RNA replication modes, ranging from the stamping machine replication (SMR) to the geometric replication (GR) mode. Although previous theoretical work has focused on the evolutionary dynamics of RNA viruses amplifying their genomes with different strategies, little is known in terms of the bifurcations and transitions involving the so-called error threshold (mutation-induced dominance of mutants) and lethal mutagenesis (extinction of all sequences due to mutation accumulation and demographic stochasticity). Here we analyze a dynamical system describing the intracellular amplification of viral RNA genomes evolving on a single-peak fitness landscape focusing on three cases considering neutral, deleterious, and lethal mutants. We analytically derive the critical mutation rates causing lethal mutagenesis and error threshold, governed by transcritical bifurcations that depend on parameters α (parameter introducing the mode of replication), replicative fitness of mutants (k₁), and on the spontaneous degradation rates of the sequences (ϵ). Our results relate the error catastrophe with lethal mutagenesis in a model with continuous populations of viral genomes. The former case involves dominance of the mutant sequences, while the latter, a deterministic extinction of the viral RNAs during replication due to increased mutation. For the lethal case the critical mutation rate involving lethal mutagenesis is μ_c=1-ɛ/α. Here, the SMR involves lower critical mutation rates, being the system more robust to lethal mutagenesis replicating closer to the GR mode. This result is also found for the neutral and deleterious cases, but for these later cases lethal mutagenesis can shift to the error threshold once the replication mode surpasses a threshold given by α=ϵ/k₁.

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.

Orally Efficacious Broad-Spectrum Ribonucleoside Analog Inhibitor of Influenza and Respiratory Syncytial Viruses

Morbidity and mortality resulting from influenza-like disease are a threat, especially for older adults. To improve case management, next-generation broad-spectrum antiviral therapeutics that are efficacious against major drivers of influenza-like disease, including influenza viruses and respiratory syncytial virus (RSV), are urgently needed. Using a dual-pathogen high-throughput screening protocol for influenza A virus (IAV) and RSV inhibitors, we have identified N⁴-hydroxycytidine (NHC) as a potent inhibitor of RSV, influenza B viruses, and IAVs of human, avian, and swine origins. Biochemical in vitro polymerase assays and viral RNA sequencing revealed that the ribonucleotide analog is incorporated into nascent viral RNAs in place of cytidine, increasing the frequency of viral mutagenesis. Viral passaging in cell culture in the presence of an inhibitor did not induce robust resistance. Pharmacokinetic profiling demonstrated dose-dependent oral bioavailability of 36 to 56%, sustained levels of the active 5'-triphosphate anabolite in primary human airway cells and mouse lung tissue, and good tolerability after extended dosing at 800 mg/kg of body weight/day. The compound was orally efficacious against RSV and both seasonal and highly pathogenic avian IAVs in mouse models, reducing lung virus loads and alleviating disease biomarkers. Oral dosing reduced IAV burdens in a guinea pig transmission model and suppressed virus spread to uninfected contact animals through direct transmission. Based on its broad-spectrum efficacy and pharmacokinetic properties, NHC is a promising candidate for future clinical development as a treatment option for influenza-like diseases.

Evidence for the Selective Basis of Transition-to-Transversion Substitution Bias in Two RNA Viruses

The substitution rates of transitions are higher than expected by chance relative to those of transversions. Many have argued that selection disfavors transversions, as nonsynonymous transversions are less likely to conserve biochemical properties of the original amino acid. Only recently has it become feasible to directly test this selective hypothesis by comparing the fitness effects of a large number of transition and transversion mutations. For example, a recent study of six viruses and one beta-lactamase gene did not find evidence supporting the selective hypothesis. Here, we analyze the relative fitness effects of transition and transversion mutations from our recently published genome-wide study of mutational fitness effects in influenza virus. In contrast to prior work, we find that transversions are significantly more detrimental than transitions. Using what we believe to be an improved statistical framework, we also identify a similar trend in two HIV data sets. We further demonstrate a fitness difference in transition and transversion mutations using four deep mutational scanning data sets of influenza virus and HIV, which provided adequate statistical power. We find that three of the most commonly cited radical/conservative amino acid categories are predictive of fitness, supporting their utility in studies of positive selection and codon usage bias. We conclude that selection is a major contributor to the transition:transversion substitution bias in viruses and that this effect is only partially explained by the greater likelihood of transversion mutations to cause radical as opposed to conservative amino acid changes.

An Integrative Approach to Virus-Host Protein-Protein Interactions

Since cell regulation and protein expression can be dramatically altered upon infection by viruses, studying the mechanisms by which viruses infect cells and the regulatory networks they disrupt is essential to understanding viral pathogenicity. This line of study can also lead to discoveries about the workings of host cells themselves. Computational methods are rapidly being developed to investigate viral-host interactions, and here we highlight recent methods and the insights that they have revealed so far, with a particular focus on methods that integrate different types of data. We also review the challenges of working with viruses compared with traditional cellular biology, and the limitations of current experimental and informatics methods.

Limits of variation, specific infectivity, and genome packaging of massively recoded poliovirus genomes

Computer design and chemical synthesis generated viable variants of poliovirus type 1 (PV1), whose ORF (6,189 nucleotides) carried up to 1,297 "Max" mutations (excess of overrepresented synonymous codon pairs) or up to 2,104 "SD" mutations (randomly scrambled synonymous codons). "Min" variants (excess of underrepresented synonymous codon pairs) are nonviable except for P2Min, a variant temperature-sensitive at 33 and 39.5 °C. Compared with WT PV1, P2Min displayed a vastly reduced specific infectivity (si) (WT, 1 PFU/118 particles vs. P2Min, 1 PFU/35,000 particles), a phenotype that will be discussed broadly. Si of haploid PV presents cellular infectivity of a single genotype. We performed a comprehensive analysis of sequence and structures of the PV genome to determine if evolutionary conserved cis-acting packaging signal(s) were preserved after recoding. We showed that conserved synonymous sites and/or local secondary structures that might play a role in determining packaging specificity do not survive codon pair recoding. This makes it unlikely that numerous "cryptic, sequence-degenerate, dispersed RNA packaging signals mapping along the entire viral genome" [Patel N, et al. (2017) Nat Microbiol 2:17098] play the critical role in poliovirus packaging specificity. Considering all available evidence, we propose a two-step assembly strategy for +ssRNA viruses: step I, acquisition of packaging specificity, either (a) by specific recognition between capsid protein(s) and replication proteins (poliovirus), or (b) by the high affinity interaction of a single RNA packaging signal (PS) with capsid protein(s) (most +ssRNA viruses so far studied); step II, cocondensation of genome/capsid precursors in which an array of hairpin structures plays a role in virion formation.

Epistatic Interactions within the Influenza A Virus Polymerase Complex Mediate Mutagen Resistance and Replication Fidelity

RNA viruses exist as genetically diverse populations. This standing genetic diversity gives them the potential to adapt rapidly, evolve resistance to antiviral therapeutics, and evade immune responses. Viral mutants with altered mutation rates or mutational tolerance have provided insights into how genetic diversity arises and how it affects the behavior of RNA viruses. To this end, we identified variants within the polymerase complex of influenza virus that are able tolerate drug-mediated increases in viral mutation rates. We find that drug resistance is highly dependent on interactions among mutations in the polymerase complex. In contrast to other viruses, influenza virus counters the effect of higher mutation rates primarily by maintaining high levels of genome replication. These findings suggest the importance of maintaining large population sizes for viruses with high mutation rates and show that multiple proteins can affect both mutation rate and genome synthesis.

Lethal mutagenesis is a broad-spectrum antiviral strategy that employs mutagenic nucleoside analogs to exploit the high mutation rate and low mutational tolerance of many RNA viruses. Studies of mutagen-resistant viruses have identified determinants of replicative fidelity and the importance of mutation rate to viral population dynamics. We have previously demonstrated the effective lethal mutagenesis of influenza A virus using three nucleoside analogs as well as the virus’s high genetic barrier to mutagen resistance. Here, we investigate the mutagen-resistant phenotypes of mutations that were enriched in drug-treated populations. We find that PB1 T123A has higher replicative fitness than the wild type, PR8, and maintains its level of genome production during 5-fluorouracil (2,4-dihydroxy-5-fluoropyrimidine) treatment. Surprisingly, this mutagen-resistant variant also has an increased baseline rate of C-to-U and G-to-A mutations. A second drug-selected mutation, PA T97I, interacts epistatically with PB1 T123A to mediate high-level mutagen resistance, predominantly by limiting the inhibitory effect of nucleosides on polymerase activity. Consistent with the importance of epistatic interactions in the influenza virus polymerase, our data suggest that nucleoside analog resistance and replication fidelity are strain dependent. Two previously identified ribavirin {1-[(2R,3R,4S,5R)-3,4-dihydroxy-5-(hydroxymethyl)oxolan-2-yl]-1H-1,2,4-triazole-3-carboxamide} resistance mutations, PB1 V43I and PB1 D27N, do not confer drug resistance in the PR8 background, and the PR8-PB1 V43I polymerase exhibits a normal baseline mutation rate. Our results highlight the genetic complexity of the influenza A virus polymerase and demonstrate that increased replicative capacity is a mechanism by which an RNA virus can counter the negative effects of elevated mutation rates.

IMPORTANCE RNA viruses exist as genetically diverse populations. This standing genetic diversity gives them the potential to adapt rapidly, evolve resistance to antiviral therapeutics, and evade immune responses. Viral mutants with altered mutation rates or mutational tolerance have provided insights into how genetic diversity arises and how it affects the behavior of RNA viruses. To this end, we identified variants within the polymerase complex of influenza virus that are able to tolerate drug-mediated increases in viral mutation rates. We find that drug resistance is highly dependent on interactions among mutations in the polymerase complex. In contrast to other viruses, influenza virus counters the effect of higher mutation rates primarily by maintaining high levels of genome replication. These findings suggest the importance of maintaining large population sizes for viruses with high mutation rates and show that multiple proteins can affect both mutation rate and genome synthesis.

A novel twelve class fluctuation test reveals higher than expected mutation rates for influenza A viruses

Influenza virus’ low replicative fidelity contributes to its capacity for rapid evolution. Clonal sequencing and fluctuation tests have suggested that the influenza virus mutation rate is 2.7 × 10^–6 - 3.0 × 10^–5 substitutions per nucleotide per strand copied (s/n/r). However, sequencing assays are biased toward mutations with minimal fitness impacts and fluctuation tests typically investigate only a subset of all possible single nucleotide mutations. We developed a fluctuation test based on reversion to fluorescence in a set of virally encoded mutant green fluorescent proteins, which allowed us to measure the rates of selectively neutral mutations representative of the twelve different mutation types. We measured an overall mutation rate of 1.8 × 10^–4 s/n/r for PR8 (H1N1) and 2.5 × 10^–4 s/n/r for Hong Kong 2014 (H3N2) and a transitional bias of 2.7–3.6. Our data suggest that each replicated genome will have an average of 2–3 mutations and highlight the importance of mutational load in influenza virus evolution.

DOI:http://dx.doi.org/10.7554/eLife.26437.001

Progress of small molecular inhibitors in the development of anti-influenza virus agents

The influenza pandemic is a major threat to human health, and highly aggressive strains such as H1N1, H5N1 and H7N9 have emphasized the need for therapeutic strategies to combat these pathogens. Influenza anti-viral agents, especially active small molecular inhibitors play important roles in controlling pandemics while vaccines are developed. Currently, only a few drugs, which function as influenza neuraminidase (NA) inhibitors and M2 ion channel protein inhibitors, are approved in clinical. However, the acquired resistance against current anti-influenza drugs and the emerging mutations of influenza virus itself remain the major challenging unmet medical needs for influenza treatment. It is highly desirable to identify novel anti-influenza agents. This paper reviews the progress of small molecular inhibitors act as antiviral agents, which include hemagglutinin (HA) inhibitors, RNA-dependent RNA polymerase (RdRp) inhibitors, NA inhibitors and M2 ion channel protein inhibitors etc. Moreover, we also summarize new, recently reported potential targets and discuss strategies for the development of new anti-influenza virus drugs.

Distinct Effects of T-705 (Favipiravir) and Ribavirin on Influenza Virus Replication and Viral RNA Synthesis

T-705 (favipiravir) is a new antiviral agent in advanced clinical development for influenza therapy. It is supposed to act as an alternative substrate for the viral polymerase, causing inhibition of viral RNA synthesis or virus mutagenesis. These mechanisms were also proposed for ribavirin, an established and broad antiviral drug that shares structural similarity with T-705. We here performed a comparative analysis of the effects of T-705 and ribavirin on influenza virus and host cell functions. Influenza virus-infected cell cultures were exposed to T-705 or ribavirin during single or serial virus passaging. The effects on viral RNA synthesis and infectious virus yield were determined and mutations appearing in the viral genome were detected by whole-genome virus sequencing. In addition, the cellular nucleotide pools as well as direct inhibition of the viral polymerase enzyme were quantified. We demonstrate that the anti-influenza virus effect of ribavirin is based on IMP dehydrogenase inhibition, which results in fast and profound GTP depletion and an imbalance in the nucleotide pools. In contrast, T-705 acts as a potent and GTP-competitive inhibitor of the viral polymerase. In infected cells, viral RNA synthesis is completely inhibited by T-705 or ribavirin at ≥50 μM, whereas exposure to lower drug concentrations induces formation of noninfectious particles and accumulation of random point mutations in the viral genome. This mutagenic effect is 2-fold higher for T-705 than for ribavirin. Hence, T-705 and ribavirin both act as purine pseudobases but profoundly differ with regard to the mechanism behind their antiviral and mutagenic effects on influenza virus.

A Novel Endonuclease Inhibitor Exhibits Broad-Spectrum Anti-Influenza Virus Activity In Vitro

Antiviral drugs are important in preventing and controlling influenza, particularly when vaccines are ineffective or unavailable. A single class of antiviral drugs, the neuraminidase inhibitors (NAIs), is recommended for treating influenza. The limited therapeutic options and the potential risk of antiviral resistance are driving the search for additional small-molecule inhibitors that act on influenza virus proteins. The acid polymerase (PA) of influenza viruses is a promising target for new antivirals because of its essential role in initiating virus transcription. Here, we characterized a novel compound, RO-7, identified as a putative PA endonuclease inhibitor. RO-7 was effective when added before the cessation of genome replication, reduced polymerase activity in cell-free systems, and decreased relative amounts of viral mRNA and genomic RNA during influenza virus infection. RO-7 specifically inhibited the ability of the PA endonuclease domain to cleave a nucleic acid substrate. RO-7 also inhibited influenza A viruses (seasonal and 2009 pandemic H1N1 and seasonal H3N2) and B viruses (Yamagata and Victoria lineages), zoonotic viruses (H5N1, H7N9, and H9N2), and NAI-resistant variants in plaque reduction, yield reduction, and cell viability assays in Madin-Darby canine kidney (MDCK) cells with nanomolar to submicromolar 50% effective concentrations (EC50s), low toxicity, and favorable selective indices. RO-7 also inhibited influenza virus replication in primary normal human bronchial epithelial cells. Overall, RO-7 exhibits broad-spectrum activity against influenza A and B viruses in multiple in vitro assays, supporting its further characterization and development as a potential antiviral agent for treating influenza.

The Mutational Robustness of Influenza A Virus

A virus’ mutational robustness is described in terms of the strength and distribution of the mutational fitness effects, or MFE. The distribution of MFE is central to many questions in evolutionary theory and is a key parameter in models of molecular evolution. Here we define the mutational fitness effects in influenza A virus by generating 128 viruses, each with a single nucleotide mutation. In contrast to mutational scanning approaches, this strategy allowed us to unambiguously assign fitness values to individual mutations. The presence of each desired mutation and the absence of additional mutations were verified by next generation sequencing of each stock. A mutation was considered lethal only after we failed to rescue virus in three independent transfections. We measured the fitness of each viable mutant relative to the wild type by quantitative RT-PCR following direct competition on A549 cells. We found that 31.6% of the mutations in the genome-wide dataset were lethal and that the lethal fraction did not differ appreciably between the HA- and NA-encoding segments and the rest of the genome. Of the viable mutants, the fitness mean and standard deviation were 0.80 and 0.22 in the genome-wide dataset and best modeled as a beta distribution. The fitness impact of mutation was marginally lower in the segments coding for HA and NA (0.88 ± 0.16) than in the other 6 segments (0.78 ± 0.24), and their respective beta distributions had slightly different shape parameters. The results for influenza A virus are remarkably similar to our own analysis of CirSeq-derived fitness values from poliovirus and previously published data from other small, single stranded DNA and RNA viruses. These data suggest that genome size, and not nucleic acid type or mode of replication, is the main determinant of viral mutational fitness effects.

Like other RNA viruses, influenza virus has a very high mutation rate. While high mutation rates may increase the rate at which influenza virus will adapt to a new host, acquire a new route of transmission, or escape from host immune surveillance, data from model systems suggest that most new viral mutations are either lethal or highly detrimental. Mutational robustness refers to the ability of a virus to tolerate, or buffer, these mutations. The mutational robustness of a virus will determine which mutations are maintained in a population and may have a greater impact on viral evolution than mutation rate. We defined the mutational robustness of influenza A virus by measuring the fitness of a large number of viruses, each with a single point mutation. We found that the overall robustness of influenza was similar to that of poliovirus and other viruses of similar size. Interestingly, mutations appeared to be more easily accommodated in hemagglutinin and neuraminidase than elsewhere in the genome. This work will inform models of influenza evolution at the global and molecular scale.

The Influenza Virus Polymerase Complex: An Update on Its Structure, Functions, and Significance for Antiviral Drug Design

Influenza viruses cause seasonal epidemics and pandemic outbreaks associated with significant morbidity and mortality, and a huge cost. Since resistance to the existing anti‐influenza drugs is rising, innovative inhibitors with a different mode of action are urgently needed. The influenza polymerase complex is widely recognized as a key drug target, given its critical role in virus replication and high degree of conservation among influenza A (of human or zoonotic origin) and B viruses. We here review the major progress that has been made in recent years in unravelling the structure and functions of this protein complex, enabling structure‐aided drug design toward the core regions of the PA endonuclease, PB1 polymerase, or cap‐binding PB2 subunit. Alternatively, inhibitors may target a protein–protein interaction site, a cellular factor involved in viral RNA synthesis, the viral RNA itself, or the nucleoprotein component of the viral ribonucleoprotein. The latest advances made for these diverse pharmacological targets have yielded agents in advanced (i.e., favipiravir and VX‐787) or early clinical testing, besides several experimental inhibitors in various stages of development, which are all covered here.

Homology-Based Identification of a Mutation in the Coronavirus RNA-Dependent RNA Polymerase That Confers Resistance to Multiple Mutagens

Positive-sense RNA viruses encode RNA-dependent RNA polymerases (RdRps) essential for genomic replication. With the exception of the large nidoviruses, such as coronaviruses (CoVs), RNA viruses lack proofreading and thus are dependent on RdRps to control nucleotide selectivity and fidelity. CoVs encode a proofreading exonuclease in nonstructural protein 14 (nsp14-ExoN), which confers a greater-than-10-fold increase in fidelity compared to other RNA viruses. It is unknown to what extent the CoV polymerase (nsp12-RdRp) participates in replication fidelity. We sought to determine whether homology modeling could identify putative determinants of nucleotide selectivity and fidelity in CoV RdRps. We modeled the CoV murine hepatitis virus (MHV) nsp12-RdRp structure and superimposed it on solved picornaviral RdRp structures. Fidelity-altering mutations previously identified in coxsackie virus B3 (CVB3) were mapped onto the nsp12-RdRp model structure and then engineered into the MHV genome with [nsp14-ExoN(+)] or without [nsp14-ExoN(-)] ExoN activity. Using this method, we identified two mutations conferring resistance to the mutagen 5-fluorouracil (5-FU): nsp12-M611F and nsp12-V553I. For nsp12-V553I, we also demonstrate resistance to the mutagen 5-azacytidine (5-AZC) and decreased accumulation of mutations. Resistance to 5-FU, and a decreased number of genomic mutations, was effectively masked by nsp14-ExoN proofreading activity. These results indicate that nsp12-RdRp likely functions in fidelity regulation and that, despite low sequence conservation, some determinants of RdRp nucleotide selectivity are conserved across RNA viruses. The results also indicate that, with regard to nucleotide selectivity, nsp14-ExoN is epistatic to nsp12-RdRp, consistent with its proposed role in a multiprotein replicase-proofreading complex.

IMPORTANCE

RNA viruses have evolutionarily fine-tuned replication fidelity to balance requirements for genetic stability and diversity. Responsibility for replication fidelity in RNA viruses has been attributed to the RNA-dependent RNA polymerases, with mutations in RdRps for multiple RNA viruses shown to alter fidelity and attenuate virus replication and virulence. Coronaviruses (CoVs) are the only known RNA viruses to encode a proofreading exonuclease (nsp14-ExoN), as well as other replicase proteins involved in regulation of fidelity. This report shows that the CoV RdRp (nsp12) likely functions in replication fidelity; that residue determinants of CoV RdRp nucleotide selectivity map to similar structural regions of other, unrelated RNA viral polymerases; and that for CoVs, the proofreading activity of the nsp14-ExoN is epistatic to the function of the RdRp in fidelity.

Impact of increased mutagenesis on adaptation to high temperature in bacteriophage Qβ

RNA viruses replicate with very high error rates, which makes them more sensitive to additional increases in this parameter. This fact has inspired an antiviral strategy named lethal mutagenesis, which is based on the artificial increase of the error rate above a threshold incompatible with virus infectivity. A relevant issue concerning lethal mutagenesis is whether incomplete treatments might enhance the adaptive possibilities of viruses. We have addressed this question by subjecting an RNA virus, the bacteriophage Qβ, to different transmission regimes in the presence or the absence of sublethal concentrations of the mutagenic nucleoside analogue 5-azacytidine (AZC). Populations obtained were subsequently exposed to a non-optimal temperature and analyzed to determine their consensus sequences. Our results show that previously mutagenized populations rapidly fixed a specific set of mutations upon propagation at the new temperature, suggesting that the expansion of the mutant spectrum caused by AZC has an influence on later evolutionary behavior.

Synergistic antiviral activity of ribavirin and interferon-α against parrot bornaviruses in avian cells

Avian bornaviruses are the causative agents of proventricular dilatation disease (PDD), a widely distributed and often fatal disease in captive psittacines. Because neither specific prevention measures nor therapies against PDD and bornavirus infections are currently available, new antiviral strategies are required to improve animal health. We show here that the nucleoside analogue ribavirin inhibited bornavirus activity in a polymerase reconstitution assay and reduced viral load in avian cell lines infected with two different parrot bornaviruses. Furthermore, we observed that ribavirin enhanced type I IFN signalling in avian cells. Combined treatment of avian bornavirus-infected cells with ribavirin and recombinant IFN-α strongly enhanced the antiviral efficiency compared to either drug alone. The combined use of ribavirin and type I IFN might represent a promising new strategy for therapeutic treatment of captive parrots persistently infected with avian bornaviruses.

Measurements of Intrahost Viral Diversity Are Extremely Sensitive to Systematic Errors in Variant Calling

With next-generation sequencing technologies, it is now feasible to efficiently sequence patient-derived virus populations at a depth of coverage sufficient to detect rare variants. However, each sequencing platform has characteristic error profiles, and sample collection, target amplification, and library preparation are additional processes whereby errors are introduced and propagated. Many studies account for these errors by using ad hoc quality thresholds and/or previously published statistical algorithms. Despite common usage, the majority of these approaches have not been validated under conditions that characterize many studies of intrahost diversity. Here, we use defined populations of influenza virus to mimic the diversity and titer typically found in patient-derived samples. We identified single-nucleotide variants using two commonly employed variant callers, DeepSNV and LoFreq. We found that the accuracy of these variant callers was lower than expected and exquisitely sensitive to the input titer. Small reductions in specificity had a significant impact on the number of minority variants identified and subsequent measures of diversity. We were able to increase the specificity of DeepSNV to >99.95% by applying an empirically validated set of quality thresholds. When applied to a set of influenza virus samples from a household-based cohort study, these changes resulted in a 10-fold reduction in measurements of viral diversity. We have made our sequence data and analysis code available so that others may improve on our work and use our data set to benchmark their own bioinformatics pipelines. Our work demonstrates that inadequate quality control and validation can lead to significant overestimation of intrahost diversity.

IMPORTANCE

Advances in sequencing technology have made it feasible to sequence patient-derived viral samples at a level sufficient for detection of rare mutations. These high-throughput, cost-effective methods are revolutionizing the study of within-host viral diversity. However, the techniques are error prone, and the methods commonly used to control for these errors have not been validated under the conditions that characterize patient-derived samples. Here, we show that these conditions affect measurements of viral diversity. We found that the accuracy of previously benchmarked analysis pipelines was greatly reduced under patient-derived conditions. By carefully validating our sequencing analysis using known control samples, we were able to identify biases in our method and to improve our accuracy to acceptable levels. Application of our modified pipeline to a set of influenza virus samples from a cohort study provided a realistic picture of intrahost diversity and suggested the need for rigorous quality control in such studies.

Antiviral drug resistance as an adaptive process

Antiviral drug resistance is a matter of great clinical importance that, historically, has been investigated mostly from a virological perspective. Although the proximate mechanisms of resistance can be readily uncovered using these methods, larger evolutionary trends often remain elusive. Recent interest by population geneticists in studies of antiviral resistance has spurred new metrics for evaluating mutation and recombination rates, demographic histories of transmission and compartmentalization, and selective forces incurred during viral adaptation to antiviral drug treatment. We present up-to-date summaries on antiviral resistance for a range of drugs and viral types, and review recent advances for studying their evolutionary histories. We conclude that information imparted by demographic and selective histories, as revealed through population genomic inference, is integral to assessing the evolution of antiviral resistance as it pertains to human health.

New strategies against drug resistance to herpes simplex virus

Herpes simplex virus (HSV), a member of the Herpesviridae family, is a significant human pathogen that results in mucocutaneous lesions in the oral cavity or genital infections. Acyclovir (ACV) and related nucleoside analogues can successfully treat HSV infections, but the emergence of drug resistance to ACV has created a barrier for the treatment of HSV infections, especially in immunocompromised patients. There is an urgent need to explore new and effective tactics to circumvent drug resistance to HSV. This review summarises the current strategies in the development of new targets (the DNA helicase/primase (H/P) complex), new types of molecules (nature products) and new antiviral mechanisms (lethal mutagenesis of Janus-type nucleosides) to fight the drug resistance of HSV.

Trends in Antiviral Strategies

Viral populations are true moving targets regarding the genomic sequences to be targeted in antiviral designs. Experts from different fields have expressed the need of new paradigms for antiviral interventions and viral disease control. This chapter reviews several strategies that aim at counteracting the adaptive capacity of viral quasispecies. The proposed designs are based on combinations of different antiviral drugs and immune modulators, or in the administration of virus-specific mutagenic agents, in an approach termed lethal mutagenesis of viruses. It consists of decreasing viral fitness by an excess of mutations that render viral proteins sub-optimal or non-functional. Viral extinction by lethal mutagenesis involves several sequential, overlapping steps that recapitulate the major concepts of intra-population interactions and genetic information stability discussed in preceding chapters. Despite the magnitude of the challenge, the chapter closes with some optimistic prospects for an effective control of viruses displaying error-prone replication, based on the combined targeting of replication fidelity and the induction of the innate immune response.

Mutations in coronavirus nonstructural protein 10 decrease virus replication fidelity

Coronaviruses (CoVs) are unique in encoding a 3'→5' exoribonuclease within nonstructural protein 14 (nsp14-ExoN) that is required for high-fidelity replication, likely via proofreading. nsp14 associates with the CoV RNA-dependent RNA polymerase (nsp12-RdRp), and nsp14-ExoN activity is enhanced by binding nsp10, a small nonenzymatic protein. However, it is not known whether nsp10 functions in the regulation of CoV replication fidelity. To test this, we engineered single and double alanine substitution mutations into the genome of murine hepatitis virus (MHV-A59) containing ExoN activity [ExoN(+)] at positions within nsp10 known to disrupt the nsp10-nsp14 interaction in vitro. We show that an nsp10 mutant, R80A/E82A-ExoN(+), was five to ten times more sensitive to treatment with the RNA mutagen 5-fluorouracil (5-FU) than wild-type (WT)-ExoN(+), suggestive of decreased replication fidelity. This decreased-fidelity phenotype was confirmed using two additional nucleoside analogs, 5-azacytidine and ribavirin. R80A/E82A-ExoN(+) reached a peak titer similar to and demonstrated RNA synthesis kinetics comparable to those seen with WT-ExoN(+). No change in 5-FU sensitivity was observed for R80A/E82A-ExoN(-) relative to MHV-ExoN(-), indicating that the decreased-fidelity phenotype of R80A/E82A-ExoN(-) is linked to the presence of ExoN activity. Our results demonstrate that nsp10 is important for CoV replication fidelity and support the hypothesis that nsp10 functions to regulate nsp14-ExoN activity during virus replication.

IMPORTANCE

The adaptive capacity of CoVs, as well as all other RNA viruses, is partially attributed to the presence of extensive population genetic diversity. However, decreased fidelity is detrimental to CoV replication and virulence; mutant CoVs with decreased replication fidelity are attenuated and more sensitive to inhibition by RNA mutagens. Thus, identifying the viral protein determinants of CoV fidelity is important for understanding CoV replication, pathogenesis, and virulence. In this report, we show that nsp10, a small, nonenzymatic viral protein, contributes to CoV replication fidelity. Our data support the hypothesis that CoVs have evolved multiple proteins, in addition to nsp14-ExoN, that are responsible for maintaining the integrity of the largest known RNA genomes.