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

Deep mutational scanning reveals the functional constraints and evolutionary potential of the influenza A virus PB1 protein

IMPORTANCE

The influenza virus polymerase is important for adaptation to new hosts and, as a determinant of mutation rate, for the process of adaptation itself. We performed a deep mutational scan of the polymerase basic 1 (PB1) protein to gain insights into the structural and functional constraints on the influenza RNA-dependent RNA polymerase. We find that PB1 is highly constrained at specific sites that are only moderately predicted by the global structure or larger domain. We identified a number of beneficial mutations, many of which have been shown to be functionally important or observed in influenza virus' natural evolution. Overall, our atlas of PB1 mutations and their fitness impacts serves as an important resource for future studies of influenza replication and evolution.

Mutual information networks reveal evolutionary relationships within the influenza A virus polymerase

The influenza A virus (IAV) RNA polymerase is an essential driver of IAV evolution. Mutations that the polymerase introduces into viral genome segments during replication are the ultimate source of genetic variation, including within the three subunits of the IAV polymerase (polymerase basic protein 2, polymerase basic protein 1, and polymerase acidic protein). Evolutionary analysis of the IAV polymerase is complicated, because changes in mutation rate, replication speed, and drug resistance involve epistatic interactions among its subunits. In order to study the evolution of the human seasonal H3N2 polymerase since the 1968 pandemic, we identified pairwise evolutionary relationships among ∼7000 H3N2 polymerase sequences using mutual information (MI), which measures the information gained about the identity of one residue when a second residue is known. To account for uneven sampling of viral sequences over time, we developed a weighted MI (wMI) metric and demonstrate that wMI outperforms raw MI through simulations using a well-sampled severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) dataset. We then constructed wMI networks of the H3N2 polymerase to extend the inherently pairwise wMI statistic to encompass relationships among larger groups of residues. We included hemagglutinin (HA) in the wMI network to distinguish between functional wMI relationships within the polymerase and those potentially due to hitch-hiking on antigenic changes in HA. The wMI networks reveal coevolutionary relationships among residues with roles in replication and encapsidation. Inclusion of HA highlighted polymerase-only subgraphs containing residues with roles in the enzymatic functions of the polymerase and host adaptability. This work provides insight into the factors that drive and constrain the rapid evolution of influenza viruses.

Evolution of the SARS-CoV-2 Mutational Spectrum

SARS-CoV-2 evolves rapidly in part because of its high mutation rate. Here, we examine whether this mutational process itself has changed during viral evolution. To do this, we quantify the relative rates of different types of single-nucleotide mutations at 4-fold degenerate sites in the viral genome across millions of human SARS-CoV-2 sequences. We find clear shifts in the relative rates of several types of mutations during SARS-CoV-2 evolution. The most striking trend is a roughly 2-fold decrease in the relative rate of G→T mutations in Omicron versus early clades, as was recently noted by Ruis et al. (2022. Mutational spectra distinguish SARS-CoV-2 replication niches. bioRxiv, doi:10.1101/2022.09.27.509649). There is also a decrease in the relative rate of C→T mutations in Delta, and other subtle changes in the mutation spectrum along the phylogeny. We speculate that these changes in the mutation spectrum could arise from viral mutations that affect genome replication, packaging, and antagonization of host innate-immune factors, although environmental factors could also play a role. Interestingly, the mutation spectrum of Omicron is more similar than that of earlier SARS-CoV-2 clades to the spectrum that shaped the long-term evolution of sarbecoviruses. Overall, our work shows that the mutation process is itself a dynamic variable during SARS-CoV-2 evolution and suggests that human SARS-CoV-2 may be trending toward a mutation spectrum more similar to that of other animal sarbecoviruses.

Evolution of the SARS-CoV-2 mutational spectrum

SARS-CoV-2 evolves rapidly in part because of its high mutation rate. Here we examine whether this mutational process itself has changed during viral evolution. To do this, we quantify the relative rates of different types of single nucleotide mutations at four-fold degenerate sites in the viral genome across millions of human SARS-CoV-2 sequences. We find clear shifts in the relative rates of several types of mutations during SARS-CoV-2 evolution. The most striking trend is a roughly two-fold decrease in the relative rate of G→T mutations in Omicron versus early clades, as was recently noted by Ruis et al (2022). There is also a decrease in the relative rate of C→T mutations in Delta, and other subtle changes in the mutation spectrum along the phylogeny. We speculate that these changes in the mutation spectrum could arise from viral mutations that affect genome replication, packaging, and antagonization of host innate-immune factors-although environmental factors could also play a role. Interestingly, the mutation spectrum of Omicron is more similar than that of earlier SARS-CoV-2 clades to the spectrum that shaped the long-term evolution of sarbecoviruses. Overall, our work shows that the mutation process is itself a dynamic variable during SARS-CoV-2 evolution, and suggests that human SARS-CoV-2 may be trending towards a mutation spectrum more similar to that of other animal sarbecoviruses.

Next-generation T cell-activating vaccination increases influenza virus mutation prevalence

To determine the potential for viral adaptation to T cell responses, we probed the full influenza virus genome by next-generation sequencing directly ex vivo from infected mice, in the context of an experimental T cell-based vaccine, an H5N1-based viral vectored vaccinia vaccine Wyeth/IL-15/5Flu, versus the current standard-of-care, seasonal inactivated influenza vaccine (IIV) and unvaccinated conditions. Wyeth/IL-15/5Flu vaccination was coincident with increased mutation incidence and frequency across the influenza genome; however, mutations were not enriched within T cell epitope regions, but high allele frequency mutations within conserved hemagglutinin stem regions and PB2 mammalian adaptive mutations arose. Depletion of CD4⁺ and CD8⁺ T cell subsets led to reduced frequency of mutants in vaccinated mice; therefore, vaccine-mediated T cell responses were important drivers of virus diversification. Our findings suggest that Wyeth/IL-15/5Flu does not generate T cell escape mutants but increases stochastic events for virus adaptation by stringent bottlenecks.

A paucigranulocytic asthma host environment promotes the emergence of virulent influenza viral variants

Influenza virus has a high mutation rate, such that within one host different viral variants can emerge. Evidence suggests that influenza virus variants are more prevalent in pregnant and/or obese individuals due to their impaired interferon response. We have recently shown that the non-allergic, paucigranulocytic subtype of asthma is associated with impaired type I interferon production. Here, we seek to address if this is associated with an increased emergence of influenza virus variants. Compared to controls, mice with paucigranulocytic asthma had increased disease severity and an increased emergence of influenza virus variants. Specifically, PB1 mutations exclusively detected in asthmatic mice were associated with increased polymerase activity. Furthermore, asthmatic host-derived virus led to increased disease severity in wild-type mice. Taken together, these data suggest that at least a subset of patients with asthma may be more susceptible to severe influenza and may be a possible source of new influenza virus variants.

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.

Influenza Virus Polymerase Mutation Stabilizes a Foreign Gene Inserted into the Virus Genome by Enhancing the Transcription/Replication Efficiency of the Modified Segment

We previously attempted to establish a reporter influenza virus by inserting the gene for the Venus fluorescent protein into the NS segment of influenza A/Puerto Rico/8/34 (PR8, H1N1) virus to yield WT-Venus-PR8. Although the inserted Venus gene was deleted during serial passages of WT-Venus-PR8, we discovered that the PB2-E712D mutation stabilizes the Venus gene. Here, we explored the mechanisms by which Venus gene deletion occurs and how the polymerase mutation stabilizes the Venus gene. Deep sequencing analysis revealed that PB2-E712D does not cause an appreciable change in the mutation rate, suggesting that the stability of the Venus gene is not affected by polymerase fidelity. We found by using quantitative real-time PCR that WT-Venus-PR8 induces high-level interferon beta (IFN-β) expression. The induction of IFN-β expression seemed to result from the reduced transcription/replication efficiency of the modified NS segment in WT-Venus-PR8. In contrast, the transcription/replication efficiency of the modified NS segment was enhanced by the PB2-E712D mutation. Loss of the Venus gene in WT-Venus-PR8 appeared to be caused by internal deletions in the NS segment. Moreover, to further our understanding of the Venus stabilization mechanisms, we identified additional amino acid mutations in the virus polymerase complex that stabilize the Venus gene. We found that some of these amino acids are located near the template exit or the product exit of the viral polymerase, suggesting that these amino acids contribute to the stability of the Venus gene by affecting the binding affinity between the polymerase complex and the RNA template and product.IMPORTANCE The reverse genetics method of influenza virus generation has enabled us to generate recombinant viruses bearing modified viral proteins. Recombinant influenza viruses expressing foreign genes have become useful tools in basic research, and such viruses can be utilized as efficient virus vectors or multivalent vaccines. However, the insertion of a foreign gene into the influenza virus genome often impairs virus replication, and the inserted genes are unstable. Elucidation of the mechanisms of foreign gene stabilization will help us to establish useful recombinant influenza viruses.

Viral RNA-dependent RNA polymerase mutants display an altered mutation spectrum resulting in attenuation in both mosquito and vertebrate hosts

The presence of bottlenecks in the transmission cycle of many RNA viruses leads to a severe reduction of number of virus particles and this occurs multiple times throughout the viral transmission cycle. Viral replication is then necessary for regeneration of a diverse mutant swarm. It is now understood that any perturbation of the mutation frequency either by increasing or decreasing the accumulation of mutations in an RNA virus results in attenuation of the virus. To determine if altering the rate at which a virus accumulates mutations decreases the probability of a successful virus infection due to issues traversing host bottlenecks, a series of mutations in the RNA-dependent RNA polymerase of Venezuelan equine encephalitis virus (VEEV), strain 68U201, were tested for mutation rate changes. All RdRp mutants were attenuated in both the mosquito and vertebrate hosts, while showing no attenuation during in vitro infections. The rescued viruses containing these mutations showed some evidence of change in fidelity, but the phenotype was not sustained following passaging. However, these mutants did exhibit changes in the frequency of specific types of mutations. Using a model of mutation production, these changes were shown to decrease the number of stop codons generated during virus replication. This suggests that the observed mutant attenuation in vivo may be due to an increase in the number of unfit genomes, which may be normally selected against by the accumulation of stop codons. Lastly, the ability of these attenuated viruses to transition through a bottleneck in vivo was measured using marked virus clones. The attenuated viruses showed an overall reduction in the number of marked clones for both the mosquito and vertebrate hosts, as well as a reduced ability to overcome the known bottlenecks in the mosquito. This study demonstrates that any perturbation of the optimal mutation frequency whether through changes in fidelity or by alterations in the mutation frequency of specific nucleotides, has significant deleterious effects on the virus, especially in the presence of host bottlenecks.

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.

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.

Mutation and Epistasis in Influenza Virus Evolution

Influenza remains a persistent public health challenge, because the rapid evolution of influenza viruses has led to marginal vaccine efficacy, antiviral resistance, and the annual emergence of novel strains. This evolvability is driven, in part, by the virus's capacity to generate diversity through mutation and reassortment. Because many new traits require multiple mutations and mutations are frequently combined by reassortment, epistatic interactions between mutations play an important role in influenza virus evolution. While mutation and epistasis are fundamental to the adaptability of influenza viruses, they also constrain the evolutionary process in important ways. Here, we review recent work on mutational effects and epistasis in influenza viruses.

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.