Comprehensive mapping of adaptation of the avian influenza polymerase protein PB2 to humans

Highly pathogenic avian influenza H5N1 virus infections in pinnipeds and seabirds in Uruguay: Implications for bird-mammal transmission in South America

The highly pathogenic avian influenza viruses of clade 2.3.4.4b have caused unprecedented deaths in South American wild birds, poultry, and marine mammals. In September 2023, pinnipeds and seabirds appeared dead on the Uruguayan Atlantic coast. Sixteen influenza virus strains were characterized by real-time reverse transcription PCR and genome sequencing in samples from sea lions (Otaria flavescens), fur seals (Arctocephalus australis), and terns (Sterna hirundinacea). Phylogenetic and ancestral reconstruction analysis showed that these strains have pinnipeds most likely as the ancestral host, representing a recent introduction of clade 2.3.4.4b in Uruguay. The Uruguayan and closely related strains from Peru (sea lions) and Chile (sea lions and a human case) carry mammalian adaptative residues 591K and 701N in the viral polymerase basic protein 2 (PB2). Our findings suggest that clade 2.3.4.4b strains in South America may have spread from mammals to mammals and seabirds, revealing a new transmission route.

popDMS infers mutation effects from deep mutational scanning data

Deep mutational scanning (DMS) experiments provide a powerful method to measure the functional effects of genetic mutations at massive scales. However, the data generated from these experiments can be difficult to analyze, with significant variation between experimental replicates. To overcome this challenge, we developed popDMS, a computational method based on population genetics theory, to infer the functional effects of mutations from DMS data. Through extensive tests, we found that the functional effects of single mutations and epistasis inferred by popDMS are highly consistent across replicates, comparing favorably with existing methods. Our approach is flexible and can be widely applied to DMS data that includes multiple time points, multiple replicates, and different experimental conditions.

Advancing Antibody Engineering through Synthetic Evolution and Machine Learning

Abs are versatile molecules with the potential to achieve exceptional binding to target Ags, while also possessing biophysical properties suitable for therapeutic drug development. Protein display and directed evolution systems have transformed synthetic Ab discovery, engineering, and optimization, vastly expanding the number of Ab clones able to be experimentally screened for binding. Moreover, the burgeoning integration of high-throughput screening, deep sequencing, and machine learning has further augmented in vitro Ab optimization, promising to accelerate the design process and massively expand the Ab sequence space interrogated. In this Brief Review, we discuss the experimental and computational tools employed in synthetic Ab engineering and optimization. We also explore the therapeutic challenges posed by developing Abs for infectious diseases, and the prospects for leveraging machine learning-guided protein engineering to prospectively design Abs resistant to viral escape.

Computational multigene interactions in virus growth and infection spread

Viruses persist in nature owing to their extreme genetic heterogeneity and large population sizes, which enable them to evade host immune defenses, escape antiviral drugs, and adapt to new hosts. The persistence of viruses is challenging to study because mutations affect multiple virus genes, interactions among genes in their impacts on virus growth are seldom known, and measures of viral fitness are yet to be standardized. To address these challenges, we employed a data-driven computational model of cell infection by a virus. The infection model accounted for the kinetics of viral gene expression, functional gene-gene interactions, genome replication, and allocation of host cellular resources to produce progeny of vesicular stomatitis virus, a prototype RNA virus. We used this model to computationally probe how interactions among genes carrying up to eleven deleterious mutations affect different measures of virus fitness: single-cycle growth yields and multicycle rates of infection spread. Individual mutations were implemented by perturbing biophysical parameters associated with individual gene functions of the wild-type model. Our analysis revealed synergistic epistasis among deleterious mutations in their effects on virus yield; so adverse effects of single deleterious mutations were amplified by interaction. For the same mutations, multicycle infection spread indicated weak or negligible epistasis, where single mutations act alone in their effects on infection spread. These results were robust to simulation in high- and low-host resource environments. Our work highlights how different types and magnitudes of epistasis can arise for genetically identical virus variants, depending on the fitness measure. More broadly, gene-gene interactions can differently affect how viruses grow and spread.

ProteinGym: Large-Scale Benchmarks for Protein Design and Fitness Prediction

Predicting the effects of mutations in proteins is critical to many applications, from understanding genetic disease to designing novel proteins that can address our most pressing challenges in climate, agriculture and healthcare. Despite a surge in machine learning-based protein models to tackle these questions, an assessment of their respective benefits is challenging due to the use of distinct, often contrived, experimental datasets, and the variable performance of models across different protein families. Addressing these challenges requires scale. To that end we introduce ProteinGym, a large-scale and holistic set of benchmarks specifically designed for protein fitness prediction and design. It encompasses both a broad collection of over 250 standardized deep mutational scanning assays, spanning millions of mutated sequences, as well as curated clinical datasets providing high-quality expert annotations about mutation effects. We devise a robust evaluation framework that combines metrics for both fitness prediction and design, factors in known limitations of the underlying experimental methods, and covers both zero-shot and supervised settings. We report the performance of a diverse set of over 70 high-performing models from various subfields (eg., alignment-based, inverse folding) into a unified benchmark suite. We open source the corresponding codebase, datasets, MSAs, structures, model predictions and develop a user-friendly website that facilitates data access and analysis.

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.

Genetic Diversity of Type A Influenza Viruses Found in Swine Herds in Northwestern Poland from 2017 to 2019: The One Health Perspective

Influenza A viruses (IAV) are still a cause of concern for public health and veterinary services worldwide. With (-) RNA-segmented genome architecture, influenza viruses are prone to reassortment and can generate a great variety of strains, some capable of crossing interspecies barriers. Seasonal IAV strains continuously spread from humans to pigs, leading to multiple reassortation events with strains endemic to swine. Due to its high adaptability to humans, a reassortant strain based on "human-like" genes could potentially be a carrier of avian origin segments responsible for high virulence, and hence become the next pandemic strain with unseen pathogenicity. The rapid evolution of sequencing methods has provided a fast and cost-efficient way to assess the genetic diversity of IAV. In this study, we investigated the genetic diversity of swine influenza viruses (swIAVs) collected from Polish farms. A total of 376 samples were collected from 11 farms. The infection was confirmed in 112 cases. The isolates were subjected to next-generation sequencing (NGS), resulting in 93 full genome sequences. Phylogenetic analysis classified 59 isolates as genotype T (H1avN2g) and 34 isolates as genotype P (H1pdmN1pdm), all of which had an internal gene cassette (IGC) derived from the H1N1pdm09-like strain. These data are consistent with evolutionary trends in European swIAVs. The applied methodology proved to be useful in monitoring the genetic diversity of IAV at the human-animal interface.

Jointly modeling deep mutational scans identifies shifted mutational effects among SARS-CoV-2 spike homologs

Deep mutational scanning (DMS) is a high-throughput experimental technique that measures the effects of thousands of mutations to a protein. These experiments can be performed on multiple homologs of a protein or on the same protein selected under multiple conditions. It is often of biological interest to identify mutations with shifted effects across homologs or conditions. However, it is challenging to determine if observed shifts arise from biological signal or experimental noise. Here, we describe a method for jointly inferring mutational effects across multiple DMS experiments while also identifying mutations that have shifted in their effects among experiments. A key aspect of our method is to regularize the inferred shifts, so that they are nonzero only when strongly supported by the data. We apply this method to DMS experiments that measure how mutations to spike proteins from SARS-CoV-2 variants (Delta, Omicron BA.1, and Omicron BA.2) affect cell entry. Most mutational effects are conserved between these spike homologs, but a fraction have markedly shifted. We experimentally validate a subset of the mutations inferred to have shifted effects, and confirm differences of > 1,000-fold in the impact of the same mutation on spike-mediated viral infection across spikes from different SARS-CoV-2 variants. Overall, our work establishes a general approach for comparing sets of DMS experiments to identify biologically important shifts in mutational effects.

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.

Avian Influenza Virus Tropism in Humans

An influenza pandemic happens when a novel influenza A virus is able to infect and transmit efficiently to a new, distinct host species. Although the exact timing of pandemics is uncertain, it is known that both viral and host factors play a role in their emergence. Species-specific interactions between the virus and the host cell determine the virus tropism, including binding and entering cells, replicating the viral RNA genome within the host cell nucleus, assembling, maturing and releasing the virus to neighboring cells, tissues or organs before transmitting it between individuals. The influenza A virus has a vast and antigenically varied reservoir. In wild aquatic birds, the infection is typically asymptomatic. Avian influenza virus (AIV) can cross into new species, and occasionally it can acquire the ability to transmit from human to human. A pandemic might occur if a new influenza virus acquires enough adaptive mutations to maintain transmission between people. This review highlights the key determinants AIV must achieve to initiate a human pandemic and describes how AIV mutates to establish tropism and stable human adaptation. Understanding the tropism of AIV may be crucial in preventing virus transmission in humans and may help the design of vaccines, antivirals and therapeutic agents against the virus.

Joint visualization of seasonal influenza serology and phylogeny to inform vaccine composition

Seasonal influenza vaccines must be updated regularly to account for mutations that allow influenza viruses to escape our existing immunity. A successful vaccine should represent the genetic diversity of recently circulating viruses and induce antibodies that effectively prevent infection by those recent viruses. Thus, linking the genetic composition of circulating viruses and the serological experimental results measuring antibody efficacy is crucial to the vaccine design decision. Historically, genetic and serological data have been presented separately in the form of static visualizations of phylogenetic trees and tabular serological results to identify vaccine candidates. To simplify this decision-making process, we have created an interactive tool for visualizing serological data that has been integrated into Nextstrain's real-time phylogenetic visualization framework, Auspice. We show how the combined interactive visualizations may be used by decision makers to explore the relationships between complex data sets for both prospective vaccine virus selection and retrospectively exploring the performance of vaccine viruses.

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

The influenza A (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 (PB2, PB1, and PA). 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 metric (wMI) and demonstrate that wMI outperforms raw MI through simulations using a well-sampled 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 HA in the wMI network to distinguish between functional wMI relationships within the polymerase and those potentially due to hitchhiking 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.

Mutational fitness landscape of human influenza H3N2 neuraminidase

Influenza neuraminidase (NA) has received increasing attention as an effective vaccine target. However, its mutational tolerance is not well characterized. Here, the fitness effects of >6,000 mutations in human H3N2 NA are probed using deep mutational scanning. Our result shows that while its antigenic regions have high mutational tolerance, there are solvent-exposed regions with low mutational tolerance. We also find that protein stability is a major determinant of NA mutational fitness. The deep mutational scanning result correlates well with mutational fitness inferred from natural sequences using a protein language model, substantiating the relevance of our findings to the natural evolution of circulating strains. Additional analysis further suggests that human H3N2 NA is far from running out of mutations despite already evolving for >50 years. Overall, this study advances our understanding of the evolutionary potential of NA and the underlying biophysical constraints, which in turn provide insights into NA-based vaccine design.

Generation of an Attenuated Chimeric Bat Influenza A Virus Live-Vaccine Prototype

Recurring epizootic influenza A virus (IAV) infections in domestic livestock such as swine and poultry are associated with a substantial economic burden and pose a constant threat to human health. Therefore, universally applicable and safe animal vaccines are urgently needed. We recently demonstrated that a reassortment-incompatible chimeric bat H17N10 virus harboring the A/swan/Germany/R65/2006 (H5N1) surface glycoproteins hemagglutinin (HA) and neuraminidase (NA) can be efficiently used as a modified live influenza vaccine (MLIV). To ensure vaccine safety and, thus, improve the applicability of this novel MLIV for mammalian usage, we performed consecutive passaging in eggs and chickens. Following passaging, we identified mutations in the viral polymerase subunits PB2 (I382S), PB1 (Q694H and I695K), and PA (E141K). Strikingly, recombinant chimeric viruses encoding these mutations showed no growth deficiencies in avian cells but displayed impaired growth in human cells and mice. Homologous prime-boost immunization of mice with one of these avian-adapted chimeric viruses, designated rR65_mono/H17N10_EP18, elicited a strong neutralizing antibody response and conferred full protection against lethal highly pathogenic avian influenza virus (HPAIV) H5N1 challenge infection. Importantly, the insertion of the avian-adaptive mutations into the conventional avian-like A/SC35M/1980 (H7N7) and prototypic human A/PR/8/34 (H1N1) viruses led to attenuated viral growth in human cells and mice. Collectively, our data show that the polymerase mutations identified here can be utilized to further improve the safety of our novel H17N10-based MLIV candidates for future mammalian applications. IMPORTANCE Recurring influenza A virus outbreaks in livestock, particularly in swine and chickens, pose a constant threat to humans. Live attenuated influenza vaccines (LAIVs) might be a potent tool to prevent epizootic outbreaks and the resulting human IAV infections; however, LAIVs have several disadvantages, especially in terms of reassortment with circulating IAVs. Notably, the newly identified bat influenza A viruses H17N10 and H18N11 cannot reassort with conventional IAVs. Chimeric bat influenza A viruses encoding surface glycoproteins of conventional IAV subtypes might thus function as safe and applicable modified live influenza vaccines (MLIVs).

Host adaptive mutations in the 2009 H1N1 pandemic influenza A virus PA gene regulate translation efficiency of viral mRNAs via GRSF1

Avian species are the major natural reservoir from which pandemic influenza A viruses can be introduced to humans. Avian influenza A virus genes, including the three viral polymerase genes, PA, PB1 and PB2, require host-adaptive mutations to allow for viral replication and transmission in humans. Previously, PA from the 2009 pH1N1 viral polymerase was found to harbor host-adaptive mutations leading to enhanced viral polymerase activity. By quantifying translation and mRNA transcription, we found that the 2009 pH1N1 PA, and the associated host-adaptive mutations, led to greater translation efficiency. This was due to enhanced cytosolic accumulation of viral mRNA, which was dependent on the host RNA binding protein GRSF1. Mutations to the GRSF1 binding site in viral mRNA, as well as GRSF1 knockdown, reduced cytosolic accumulation and translation efficiency of viral mRNAs. This study identifies a previously unrecognized mechanism by which host-adaptive mutations in PA regulate viral replication and host adaptation. Importantly, these results provide greater insight into the host adaptation process of IAVs and reveal the importance of GRSF1 in the lifecycle of IAV.

Effect of the Interaction between Viral PB2 and Host SphK1 on H9N2 AIV Replication in Mammals

The H9N2 avian influenza virus (AIV) is currently widespread worldwide, posing a severe threat to the poultry industry and public health. Reassortment is an important way for influenza viruses to adapt to a new host. In 2007, the PB2 gene of H9N2 AIV in China was reassorted, and the DK1-like lineage replaced the F/98-like lineage, forming a dominant genotype of G57. This genotype and its reassortants (such as H7N9, H10N8 and H5N6) showed higher mammalian adaptation, and caused increased human infections. However, the adaptive mechanisms of the DK1-like lineage PB2 gene remain unclear. Here, we confirmed that the PB2 lineage of the H9N2 AIV currently prevalent in China still belongs to the DK1-like lineage and, compared with the previously predominant F/98-like lineage, the DK1-like lineage PB2 gene significantly enhances H9N2 AIV to mammalian adaptation. Through transcriptomic analysis and qRT–PCR and western blot experiments, we identified a host factor, sphingosine kinase 1 (SphK1), that is closely related to viral replication. SphK1 inhibits the replication of DK1-like PB2 gene H9N2 AIV, but the ability of SphK1 protein to bind DK1-like PB2 protein is weaker than that of F/98-like PB2 protein, which may contribute to H9N2 AIV containing the DK1-like PB2 gene to escape the inhibitory effect of host factor SphK1 for efficient infection. This study broadens our understanding of the adaptive evolution of H9N2 AIV and highlights the necessity to pay close attention to the AIV that contains the adaptive PB2 protein in animals and humans.

Evolutionary Landscapes of Host-Virus Arms Races

Vertebrate immune systems suppress viral infection using both innate restriction factors and adaptive immunity. Viruses mutate to escape these defenses, driving hosts to counterevolve to regain fitness. This cycle recurs repeatedly, resulting in an evolutionary arms race whose outcome depends on the pace and likelihood of adaptation by host and viral genes. Although viruses evolve faster than their vertebrate hosts, their proteins are subject to numerous functional constraints that impact the probability of adaptation. These constraints are globally defined by evolutionary landscapes, which describe the fitness and adaptive potential of all possible mutations. We review deep mutational scanning experiments mapping the evolutionary landscapes of both host and viral proteins engaged in arms races. For restriction factors and some broadly neutralizing antibodies, landscapes favor the host, which may help to level the evolutionary playing field against rapidly evolving viruses. We discuss the biophysical underpinnings of these landscapes and their therapeutic implications.

Swine H1N1 Influenza Virus Variants with Enhanced Polymerase Activity and HA Stability Promote Airborne Transmission in Ferrets

Diverse IAVs circulate in animals, yet few acquire the viral traits needed to start a human pandemic. A stabilized HA and mammalian-adapted polymerase have been shown to promote the adaptation of IAVs to humans and ferrets (the gold-standard model for IAV replication, pathogenicity, and transmissibility).

Antigenic evolution of human influenza H3N2 neuraminidase is constrained by charge balancing

As one of the main influenza antigens, neuraminidase (NA) in H3N2 virus has evolved extensively for more than 50 years due to continuous immune pressure. While NA has recently emerged as an effective vaccine target, biophysical constraints on the antigenic evolution of NA remain largely elusive. Here, we apply combinatorial mutagenesis and next-generation sequencing to characterize the local fitness landscape in an antigenic region of NA in six different human H3N2 strains that were isolated around 10 years apart. The local fitness landscape correlates well among strains and the pairwise epistasis is highly conserved. Our analysis further demonstrates that local net charge governs the pairwise epistasis in this antigenic region. In addition, we show that residue coevolution in this antigenic region is correlated with the pairwise epistasis between charge states. Overall, this study demonstrates the importance of quantifying epistasis and the underlying biophysical constraint for building a model of influenza evolution.

Influenza A Virus-Host Specificity: An Ongoing Cross-Talk Between Viral and Host Factors

One big threat from influenza A viruses (IAVs) is that novel viruses emerge from mutation alongside reassortment. Some of them have gained the capability to transmit into human from the avian reservoir. Understanding the molecular events and the involved factors in breaking the cross-species barrier holds important implication for the surveillance and prevention of potential influenza outbreaks. In this review, we summarize recent progresses, including several ground-breaking findings, in how the interaction between host and viral factors, exemplified by the PB2 subunit of the influenza virus RNA polymerase co-opting host ANP32 protein to facilitate transcription and replication of the viral genome, shapes the evolution of IAVs from host specificity to cross-species infection.

Amino acid sites related to the PB2 subunits of IDV affect polymerase activity

BACKGROUND

In 2011, a new influenza virus, named Influenza D Virus (IDV), was isolated from pigs, and then cattle, presenting influenza-like symptoms. IDV is one of the causative agents of Bovine Respiratory Disease (BRD), which causes high morbidity and mortality in feedlot cattle worldwide. To date, the molecular mechanisms of IDV pathogenicity are unknown. Recent IDV outbreaks in cattle, along with serological and genetic evidence of IDV infection in humans, have raised concerns regarding the zoonotic potential of this virus. Influenza virus polymerase is a determining factor of viral pathogenicity to mammals.

METHODS

Here we take a prospective approach to this question by creating a random mutation library about PB2 subunit of the IDV viral polymerase to test which amino acid point mutations will increase viral polymerase activity, leading to increased pathogenicity of the virus.

RESULTS

Our work shows some exact sites that could affect polymerase activities in influenza D viruses. For example, two single-site mutations, PB2-D533S and PB2-G603Y, can independently increase polymerase activity. The PB2-D533S mutation alone can increase the polymerase activity by 9.92 times, while the PB2-G603Y mutation increments the activity by 8.22 times.

CONCLUSION

Taken together, our findings provide important insight into IDV replication fitness mediated by the PB2 protein, increasing our understanding of IDV replication and pathogenicity and facilitating future studies.

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.

Characterization of changes in the hemagglutinin that accompanied the emergence of H3N2/1968 pandemic influenza viruses

The hemagglutinin (HA) of A/H3N2 pandemic influenza viruses (IAVs) of 1968 differed from its inferred avian precursor by eight amino acid substitutions. To determine their phenotypic effects, we studied recombinant variants of A/Hong Kong/1/1968 virus containing either human-type or avian-type amino acids in the corresponding positions of HA. The precursor HA displayed receptor binding profile and high conformational stability typical for duck IAVs. Substitutions Q226L and G228S, in addition to their known effects on receptor specificity and replication, marginally decreased HA stability. Substitutions R62I, D63N, D81N and N193S reduced HA binding avidity. Substitutions R62I, D81N and A144G promoted viral replication in human airway epithelial cultures. Analysis of HA sequences revealed that substitutions D63N and D81N accompanied by the addition of N-glycans represent common markers of avian H3 HA adaptation to mammals. Our results advance understanding of genotypic and phenotypic changes in IAV HA required for avian-to-human adaptation and pandemic emergence.

Multiple polymerase acidic (PA) I38X substitutions in influenza A(H1N1)pdm09 virus permit polymerase activity and cause reduced baloxavir inhibition

BACKGROUND

Baloxavir marboxil is an antiviral drug that targets the endonuclease activity of the influenza virus polymerase acidic (PA) protein. PA I38T/M/F substitutions reduce its antiviral efficacy.

OBJECTIVES

To understand the effects of the 19 possible amino acid (AA) substitutions at PA 38 on influenza A(H1N1)pdm09 polymerase activity and inhibition by baloxavir acid, the active metabolite of baloxavir marboxil.

METHODS

Influenza A(H1N1)pdm09 viral polymerase complexes containing all 19 I38X AA substitutions were reconstituted in HEK293T cells in a mini-replicon assay. Polymerase complex activity and baloxavir inhibitory activity were measured in the presence or absence of 50 nM baloxavir acid.

RESULTS

Only three substitutions (R, K, P) reduced polymerase activity to <79% of I38-WT. When compared with the prototypical baloxavir marboxil resistance marker T38, 5 substitutions conferred 10%-35% reductions in baloxavir acid inhibitory activity (M, L, F, Y, C) and 11 substitutions conferred >50% reductions (R, K, S, N, G, W, A, Q, E, D, H), while two substitutions (V, P) maintained baloxavir acid inhibitory activity.

CONCLUSIONS

Most PA 38 substitutions permit a functional replication complex retaining some drug resistance in the mini-replicon assay. This study provides a targeted approach for virus rescue and analysis of novel baloxavir marboxil reduced-susceptibility markers, supports the consideration of a broader range of these markers during antiviral surveillance and adds to the growing knowledge of baloxavir marboxil resistance profiles.

Comprehensive Profiling of Mutations to Influenza Virus PB2 That Confer Resistance to the Cap-Binding Inhibitor Pimodivir

Antivirals are used not only in the current treatment of influenza but are also stockpiled as a first line of defense against novel influenza strains for which vaccines have yet to be developed. Identifying drug resistance mutations can guide the clinical deployment of the antiviral and can additionally define the mechanisms of drug action and drug resistance. Pimodivir is a first-in-class inhibitor of the polymerase basic protein 2 (PB2) subunit of the influenza A virus polymerase complex. A number of resistance mutations have previously been identified in treated patients or cell culture. Here, we generate a complete map of the effect of all single-amino-acid mutations to an avian PB2 on resistance to pimodivir. We identified both known and novel resistance mutations not only in the previously implicated cap-binding and mid-link domains, but also in the N-terminal domain. Our complete map of pimodivir resistance thus enables the evaluation of whether new viral strains contain mutations that will confer pimodivir resistance.

Inhibition of Influenza Virus Polymerase by Interfering with Its Protein-Protein Interactions

Influenza (flu) virus is a serious threat to global health with the potential to generate devastating pandemics. The availability of broad spectrum antiviral drugs is an unequaled weapon during pandemic events, especially when a vaccine is still not available. One of the most promising targets for the development of new antiflu therapeutics is the viral RNA-dependent RNA polymerase (RdRP). The assembly of the flu RdRP heterotrimeric complex from the individual polymerase acidic protein (PA), polymerase basic protein 1 (PB1), and polymerase basic protein 2 (PB2) subunits is a prerequisite for RdRP functions, such as mRNA synthesis and genome replication. In this Review, we report the known protein-protein interactions (PPIs) occurring by RdRP that could be disrupted by small molecules and analyze their benefits and drawbacks as drug targets. An overview of small molecules able to interfere with flu RdRP functions exploiting the PPI inhibition approach is described. In particular, an update on the most recent inhibitors targeting the well-consolidated RdRP PA-PB1 subunit heterodimerization is mainly reported, together with pioneer inhibitors targeting other virus-virus or virus-host interactions involving RdRP subunits. As demonstrated by the PA-PB1 interaction inhibitors discussed herein, the inhibition of flu RdRP functions by PPI disrupters clearly represents a valid means to identify compounds endowed with a broad spectrum of action and a reduced propensity to develop drug resistance, which are the main issues of antiviral drugs.

Mammalian cells use the autophagy process to restrict avian influenza virus replication

Host adaptive mutations in the influenza A virus (IAV) PB2 protein are critical for human infection, but their molecular action is not well understood. We observe that when IAV containing avian PB2 infects mammalian cells, viral ribonucleoprotein (vRNP) aggregates that localize to the microtubule-organizing center (MTOC) are formed. These vRNP aggregates resemble LC3B-associated autophagosome structures, with aggresome-like properties, in that they cause the re-distribution of vimentin. However, electron microscopy reveals that these aggregates represent an accumulation of autophagic vacuoles. Compared to mammalian-PB2 virus, avian-PB2 virus induces higher autophagic flux in infected cells, indicating an increased rate of autophagosomes containing avian vRNPs fusing with lysosomes. We found that p62 is essential for the formation of vRNP aggregates and that the Raptor-interacting region of p62 is required for interaction with vRNPs through the PB2 polymerase subunit. Selective autophagic sequestration during late-stage virus replication is thus an additional strategy for host restriction of avian-PB2 IAV.

Applications of Deep Mutational Scanning in Virology

Several recently developed high-throughput techniques have changed the field of molecular virology. For example, proteomics studies reveal complete interactomes of a viral protein, genome-wide CRISPR knockout and activation screens probe the importance of every single human gene in aiding or fighting a virus, and ChIP-seq experiments reveal genome-wide epigenetic changes in response to infection. Deep mutational scanning is a relatively novel form of protein science which allows the in-depth functional analysis of every nucleotide within a viral gene or genome, revealing regions of importance, flexibility, and mutational potential. In this review, we discuss the application of this technique to RNA viruses including members of the Flaviviridae family, Influenza A Virus and Severe Acute Respiratory Syndrome Coronavirus 2. We also briefly discuss the reverse genetics systems which allow for analysis of viral replication cycles, next-generation sequencing technologies and the bioinformatics tools that facilitate this research.

A single N342D substitution in Influenza B Virus NA protein determines viral pathogenicity in mice

Influenza B virus (IBV) is one of the most important human respiratory viruses: it causes approximately one-third of the global influenza-related disease burden each year. However, compared with the several pathogenicity-related molecular markers that have been identified for influenza A virus (IAV), little is known about potential IBV pathogenicity-related markers. Here, although the IBV strain B/Anhui-Tunxi/1528/2014 (AH1528/14) exhibited a more efficient replication ability in vitro and higher pathogenicity in vivo compared with IBV strain B/Anhui-Baohe/127/2015 (AH127/15), only three amino acids differences (HA_A390E, NA_N342D and PB1_V212I) were observed among their full genomes. The contributions of each amino acid difference to the virus pathogenicity were further investigated. Compared with the wild type IBV virus rAH127, the recombinant virus harbouring a single substitution of HA_A390E had a similar phenotype, whereas the recombinant virus harbouring PB1_V212I replicated to a moderately higher titre in both MDCK cells and in mice. Notably, the virus harbouring NA_N342D showed significantly better growth properties in MDCK cells and higher fatality rates in mice. In addition, the presence of NA_N342D dramatically enhanced the viral neuraminidase activity. In conclusion, our study identified a novel IBV molecular marker, NA_N342D, that could significantly increase the virulence of IBV in mice.

Deep Mutational Scanning of Viral Glycoproteins and Their Host Receptors

Deep mutational scanning or deep mutagenesis is a powerful tool for understanding the sequence diversity available to viruses for adaptation in a laboratory setting. It generally involves tracking an in vitro selection of protein sequence variants with deep sequencing to map mutational effects based on changes in sequence abundance. Coupled with any of a number of selection strategies, deep mutagenesis can explore the mutational diversity available to viral glycoproteins, which mediate critical roles in cell entry and are exposed to the humoral arm of the host immune response. Mutational landscapes of viral glycoproteins for host cell attachment and membrane fusion reveal extensive epistasis and potential escape mutations to neutralizing antibodies or other therapeutics, as well as aiding in the design of optimized immunogens for eliciting broadly protective immunity. While less explored, deep mutational scans of host receptors further assist in understanding virus-host protein interactions. Critical residues on the host receptors for engaging with viral spikes are readily identified and may help with structural modeling. Furthermore, mutations may be found for engineering soluble decoy receptors as neutralizing agents that specifically bind viral targets with tight affinity and limited potential for viral escape. By untangling the complexities of how sequence contributes to viral glycoprotein and host receptor interactions, deep mutational scanning is impacting ideas and strategies at multiple levels for combatting circulating and emergent virus strains.

Experimental Approaches to Identify Host Factors Important for Influenza Virus

An ever-expanding toolkit of experimental methods provides the means to discover and characterize host factors important for influenza virus. Here, we describe common methods for investigating genetic relationships and physical interactions between virus and host. A comprehensive knowledge of host:virus interactions is key to understanding how influenza virus exploits the host cell and to potentially identify vulnerabilities that may be manipulated to prevent or treat disease.

Molecular basis of host-adaptation interactions between influenza virus polymerase PB2 subunit and ANP32A

Avian influenza polymerase undergoes host adaptation in order to efficiently replicate in human cells. Adaptive mutants are localised on the C-terminal (627-NLS) domains of the PB2 subunit. In particular, mutation of PB2 residue 627 from E to K rescues polymerase activity in mammalian cells. A host transcription regulator ANP32A, comprising a long C-terminal intrinsically disordered domain (IDD), is responsible for this adaptation. Human ANP32A IDD lacks a 33 residue insertion compared to avian ANP32A, and this deletion restricts avian influenza polymerase activity. We used NMR to determine conformational ensembles of E627 and K627 forms of 627-NLS of PB2 in complex with avian and human ANP32A. Human ANP32A IDD transiently binds to the 627 domain, exploiting multivalency to maximise affinity. E627 interrupts the polyvalency of the interaction, an effect compensated by an avian-unique motif in the IDD. The observed binding mode is maintained in the context of heterotrimeric influenza polymerase, placing ANP32A in the immediate vicinity of known host-adaptive PB2 mutants.

Avian influenza polymerase undergoes host adaptation in order to efficiently replicate in human cells. Here, the authors use NMR spectroscopy and quantitative ensemble modelling to describe the highly dynamic assemblies formed by the human-adapted or avian-adapted C-terminal domains with the respective ANP32A host proteins.

Quantifying within-host diversity of H5N1 influenza viruses in humans and poultry in Cambodia

Avian influenza viruses (AIVs) periodically cross species barriers and infect humans. The likelihood that an AIV will evolve mammalian transmissibility depends on acquiring and selecting mutations during spillover, but data from natural infection is limited. We analyze deep sequencing data from infected humans and domestic ducks in Cambodia to examine how H5N1 viruses evolve during spillover. Overall, viral populations in both species are predominated by low-frequency (<10%) variation shaped by purifying selection and genetic drift, and half of the variants detected within-host are never detected on the H5N1 virus phylogeny. However, we do detect a subset of mutations linked to human receptor binding and replication (PB2 E627K, HA A150V, and HA Q238L) that arose in multiple, independent humans. PB2 E627K and HA A150V were also enriched along phylogenetic branches leading to human infections, suggesting that they are likely human-adaptive. Our data show that H5N1 viruses generate putative human-adapting mutations during natural spillover infection, many of which are detected at >5% frequency within-host. However, short infection times, genetic drift, and purifying selection likely restrict their ability to evolve extensively during a single infection. Applying evolutionary methods to sequence data, we reveal a detailed view of H5N1 virus adaptive potential, and develop a foundation for studying host-adaptation in other zoonotic viruses.

H5N1 avian influenza viruses can cross species barriers and cause severe disease in humans. H5N1 viruses currently cannot replicate and transmit efficiently among humans, but animal infection studies and modeling experiments have suggested that human adaptation may require only a few mutations. However, data from natural spillover infection has been limited, posing a challenge for risk assessment. Here, we analyze a unique dataset of deep sequence data from H5N1 virus-infected humans and domestic ducks in Cambodia. We find that well-known markers of human receptor binding and replication arise in multiple, independent humans. We also find that 3 mutations detected within-host are enriched along phylogenetic branches leading to human infections, suggesting that they are likely human-adapting. However, we also show that within-host evolution in both humans and ducks are shaped heavily by purifying selection and genetic drift, and that a large fraction of within-host variation is never detected on the H5N1 phylogeny. Taken together, our data show that H5N1 viruses do generate human-adapting mutations during natural infection. However, short infection times, purifying selection, and genetic drift may severely limit how much H5N1 viruses can evolve during the course of a single infection.

Deep Mutational Scanning Comprehensively Maps How Zika Envelope Protein Mutations Affect Viral Growth and Antibody Escape

Functional constraints on viral proteins are often assessed by examining sequence conservation among natural strains, but this approach is relatively ineffective for Zika virus because all known sequences are highly similar. Here, we take an alternative approach to map functional constraints on Zika virus's envelope (E) protein by using deep mutational scanning to measure how all amino acid mutations to the E protein affect viral growth in cell culture. The resulting sequence-function map is consistent with existing knowledge about E protein structure and function but also provides insight into mutation-level constraints in many regions of the protein that have not been well characterized in prior functional work. In addition, we extend our approach to completely map how mutations affect viral neutralization by two monoclonal antibodies, thereby precisely defining their functional epitopes. Overall, our study provides a valuable resource for understanding the effects of mutations to this important viral protein and also offers a roadmap for future work to map functional and antigenic selection to Zika virus at high resolution.IMPORTANCE Zika virus has recently been shown to be associated with severe birth defects. The virus's E protein mediates its ability to infect cells and is also the primary target of the antibodies that are elicited by natural infection and vaccines that are being developed against the virus. Therefore, determining the effects of mutations to this protein is important for understanding its function, its susceptibility to vaccine-mediated immunity, and its potential for future evolution. We completely mapped how amino acid mutations to the E protein affected the virus's ability to grow in cells in the laboratory and escape from several antibodies. The resulting maps relate changes in the E protein's sequence to changes in viral function and therefore provide a valuable complement to existing maps of the physical structure of the protein.

Species-Specific Host-Virus Interactions: Implications for Viral Host Range and Virulence

A growing number of studies indicate that host species-specific and virus strain-specific interactions of viral molecules with the host innate immune system play a pivotal role in determining virus host range and virulence. Because interacting proteins are likely constrained in their evolution, mutations that are selected to improve virus replication in one species may, by chance, alter the ability of a viral antagonist to inhibit immune responses in hosts the virus has not yet encountered. Based on recent findings of host-species interactions of poxvirus, herpesvirus, and influenza virus proteins, we propose a model for viral fitness and host range which considers the full interactome between a specific host species and a virus, resulting from the combination of all interactions, positive and negative, that influence whether a virus can productively infect a cell and cause disease in different hosts.

Profiling host ANP32A splicing landscapes to predict influenza A virus polymerase adaptation

Species’ differences in cellular factors limit avian influenza A virus (IAV) zoonoses and human pandemics. The IAV polymerase, vPol, harbors evolutionary sites to overcome restriction and determines virulence. Here, we establish host ANP32A as a critical driver of selection, and identify host-specific ANP32A splicing landscapes that predict viral evolution. We find that avian species differentially express three ANP32A isoforms diverging in a vPol-promoting insert. ANP32As with shorter inserts interact poorly with vPol, are compromised in supporting avian-like IAV replication, and drive selection of mammalian-adaptive vPol sequences with distinct kinetics. By integrating selection data with multi-species ANP32A splice variant profiling, we develop a mathematical model to predict avian species potentially driving (swallow, magpie) or maintaining (goose, swan) mammalian-adaptive vPol signatures. Supporting these predictions, surveillance data confirm enrichment of several mammalian-adaptive vPol substitutions in magpie IAVs. Profiling host ANP32A splicing could enhance surveillance and eradication efforts against IAVs with pandemic potential.

Polymorphisms in the avian influenza A virus (IAV) polymerase restrict its host range during transmission from birds to mammals. Here, the authors investigate differences in the host chromatin regulator ANP32A regarding IAV polymerase adaptation, and profile ANP32A splicing to predict avian species associated with pre-adaptive human-signatures in the virus.

PB2 mutations arising during H9N2 influenza evolution in the Middle East confer enhanced replication and growth in mammals

Avian influenza virus H9N2 has been endemic in birds in the Middle East, in particular in Egypt with multiple cases of human infections since 1998. Despite concerns about the pandemic threat posed by H9N2, little is known about the biological properties of H9N2 in this epicentre of infection. Here, we investigated the evolutionary dynamics of H9N2 in the Middle East and identified phylogeny-associated PB2 mutations that acted cooperatively to increase H9N2 replication/transcription in human cells. The accumulation of PB2 mutations also correlated with an increase in H9N2 virus growth in the upper and lower airways of mice and in virulence. These mutations clustered on a solvent-exposed region in the PB2-627 domain in proximity to potential interfaces with host factors. These PB2 mutations have been found at high prevalence during evolution of H9N2 in the field, indicating that they have provided a selective advantage for viral adaptation to infect poultry. Therefore, continuous prevalence of H9N2 virus in the Middle East has generated a far more fit or optimized replication phenotype, leading to an expanded viral host range, including to mammals, which may pose public health risks beyond the current outbreaks.

The G1-like clade of H9N2 influenza viruses can undergo genetic reassortment with other influenza virus subtypes to produce novel zoonotic viruses, such as the Gs/GD lineage H5N1, H7N9, H10N8, and H5N8 viruses. Since 1998, the G1-like subclade of H9N2 influenza virus has been widely circulating in birds in Central Asia and the Middle East and a number of human cases have been reported. However, little is known about the biological properties of H9N2 viruses in this epicentre of infection. Our data showed that, during about two decades of evolution in nature, G1-like subclade strains evolved to produce strains with appreciably higher replication phenotypes in Central Asia and the Middle East, which led to their expanded host range, including to humans. Therefore, G1-like subclade strains in these areas may accumulate mutations to produce novel viruses and the large gene pool in these areas would enable reassortment with other influenza viruses. This study indicated the need for studies of H9N2 viruses in such areas to monitor their evolutionary dynamics and possible genetic changes.