Epistatic interactions between neuraminidase mutations facilitated the emergence of the oseltamivir-resistant H1N1 influenza viruses

Ancestral Reconstruction and the Evolution of Protein Energy Landscapes

A protein's sequence determines its conformational energy landscape. This, in turn, determines the protein's function. Understanding the evolution of new protein functions therefore requires understanding how mutations alter the protein energy landscape. Ancestral sequence reconstruction (ASR) has proven a valuable tool for tackling this problem. In ASR, one phylogenetically infers the sequences of ancient proteins, allowing characterization of their properties. When coupled to biophysical, biochemical, and functional characterization, ASR can reveal how historical mutations altered the energy landscape of ancient proteins, allowing the evolution of enzyme activity, altered conformations, binding specificity, oligomerization, and many other protein features. In this article, we review how ASR studies have been used to dissect the evolution of energy landscapes. We also discuss ASR studies that reveal how energy landscapes have shaped protein evolution. Finally, we propose that thinking about evolution from the perspective of an energy landscape can improve how we approach and interpret ASR studies.

Recombinant A(H6N1)-H274Y avian influenza virus with dual drug resistance does not require permissive mutations to retain the replicative fitness in vitro and in ovo

The possible emergence of drug-resistant avian flu raises concerns over the limited effectiveness of currently approved antivirals (neuraminidase inhibitors - NAIs) in the hypothetical event of a zoonotic spillover. Our study demonstrated that the recombinant avian A(H6N1) viruses showed reduced inhibition (RI) by multiple NAI drugs following the introduction of point mutations found predominantly in the neuraminidase gene (NA) of NAI-resistant human influenza strains (E119V, R292K and H274Y; N2 numbering). Moreover, A(H6N1)-H274Y showed increased replication efficiency in vitro, and a fitness advantage over wild-type (WT) when co-inoculated into embryonated hen's eggs. The results presented in our study together with the zoonotic potential of the A(H6N1) virus as evidenced by the human infection from 2013, highlight the need for enhanced monitoring of NAI resistance-associated signatures in circulating LPAI (low pathogenic avian influenza) globally.

Molecular modeling and phylogenetic analyses highlight the role of amino acid 347 of the N1 subtype neuraminidase in influenza virus host range and interspecies adaptation

The N1 neuraminidases (NAs) of avian and pandemic human influenza viruses contain tyrosine and asparagine, respectively, at position 347 on the rim of the catalytic site; the biological significance of this difference is not clear. Here, we used molecular dynamics simulation to model the effects of amino acid 347 on N1 NA interactions with sialyllacto-N-tetraoses 6'SLN-LC and 3'SLN-LC, which represent NA substrates in humans and birds, respectively. Our analysis predicted that Y347 plays an important role in the NA preference for the avian-type substrates. The Y347N substitution facilitates hydrolysis of human-type substrates by resolving steric conflicts of the Neu5Ac2-6Gal moiety with the bulky side chain of Y347, decreasing the free energy of substrate binding, and increasing the solvation of the Neu5Ac2-6Gal bond. Y347 was conserved in all N1 NA sequences of avian influenza viruses in the GISAID EpiFlu database with two exceptions. First, the Y347F substitution was present in the NA of a specific H6N1 poultry virus lineage and was associated with the substitutions G228S and/or E190V/L in the receptor-binding site (RBS) of the hemagglutinin (HA). Second, the highly pathogenic avian H5N1 viruses of the Gs/Gd lineage contained sporadic variants with the NA substitutions Y347H/D, which were frequently associated with substitutions in the HA RBS. The Y347N substitution occurred following the introductions of avian precursors into humans and pigs with N/D347 conserved during virus circulation in these hosts. Comparative evolutionary analysis of site 347 revealed episodic positive selection across the entire tree and negative selection within most host-specific groups of viruses, suggesting that substitutions at NA position 347 occurred during host switches and remained under pervasive purifying selection thereafter. Our results elucidate the role of amino acid 347 in NA recognition of sialoglycan substrates and emphasize the significance of substitutions at position 347 as a marker of host range and adaptive evolution of influenza viruses.

Co-evolution of immunity and seasonal influenza viruses

Seasonal influenza viruses cause recurring global epidemics by continually evolving to escape host immunity. The viral constraints and host immune responses that limit and drive the evolution of these viruses are increasingly well understood. However, it remains unclear how most of these advances improve the capacity to reduce the impact of seasonal influenza viruses on human health. In this Review, we synthesize recent progress made in understanding the interplay between the evolution of immunity induced by previous infections or vaccination and the evolution of seasonal influenza viruses driven by the heterogeneous accumulation of antibody-mediated immunity in humans. We discuss the functional constraints that limit the evolution of the viruses, the within-host evolutionary processes that drive the emergence of new virus variants, as well as current and prospective options for influenza virus control, including the viral and immunological barriers that must be overcome to improve the effectiveness of vaccines and antiviral drugs.

Prevalence and mechanisms of evolutionary contingency in human influenza H3N2 neuraminidase

Neuraminidase (NA) of human influenza H3N2 virus has evolved rapidly and been accumulating mutations for more than half-century. However, biophysical constraints that govern the evolutionary trajectories of NA remain largely elusive. Here, we show that among 70 natural mutations that are present in the NA of a recent human H3N2 strain, >10% are deleterious for an ancestral strain. By mapping the permissive mutations using combinatorial mutagenesis and next-generation sequencing, an extensive epistatic network is revealed. Biophysical and structural analyses further demonstrate that certain epistatic interactions can be explained by non-additive stability effect, which in turn modulates membrane trafficking and enzymatic activity of NA. Additionally, our results suggest that other biophysical mechanisms also contribute to epistasis in NA evolution. Overall, these findings not only provide mechanistic insights into the evolution of human influenza NA and elucidate its sequence-structure-function relationship, but also have important implications for the development of next-generation influenza vaccines.

Predicting Permissive Mutations That Improve the Fitness of A(H1N1)pdm09 Viruses Bearing the H275Y Neuraminidase Substitution

Oseltamivir-resistant influenza viruses arise due to amino acid mutations in key residues of the viral neuraminidase (NA). These changes often come at a fitness cost; however, it is known that permissive mutations in the viral NA can overcome this cost. This result was observed in former seasonal A(H1N1) viruses in 2007 which expressed the H275Y substitution (N1 numbering) with no apparent fitness cost and lead to widespread oseltamivir resistance. Therefore, this study aims to predict permissive mutations that may similarly enable fit H275Y variants to arise in currently circulating A(H1N1)pdm09 viruses. The first approach in this study utilized in silico analyses to predict potentially permissive mutations. The second approach involved the generation of a virus library which encompassed all possible NA mutations while keeping H275Y fixed. Fit variants were then selected by serially passaging the virus library either through ferrets by transmission or passaging once in vitro. The fitness impact of selected substitutions was further evaluated experimentally. The computational approach predicted three candidate permissive NA mutations which, in combination with each other, restored the replicative fitness of an H275Y variant. The second approach identified a stringent bottleneck during transmission between ferrets; however, three further substitutions were identified which may improve transmissibility. A comparison of fit H275Y variants in vitro and in experimentally infected animals showed a statistically significant correlation in the variants that were positively selected. Overall, this study provides valuable tools and insights into potential permissive mutations that may facilitate the emergence of a fit H275Y A(H1N1)pdm09 variant.

IMPORTANCE Oseltamivir (Tamiflu) is the most widely used antiviral for the treatment of influenza infections. Therefore, resistance to oseltamivir is a public health concern. This study is important as it explores the different evolutionary pathways available to current circulating influenza viruses that may lead to widespread oseltamivir resistance. Specifically, this study develops valuable experimental and computational tools to evaluate the fitness landscape of circulating A(H1N1)pmd09 influenza viruses bearing the H275Y mutation. The H275Y substitution is most commonly reported to confer oseltamivir resistance but also leads to loss of virus replication and transmission fitness, which limits its spread. However, it is known from previous influenza seasons that influenza viruses can evolve to overcome this loss of fitness. Therefore, this study aims to prospectively predict how contemporary A(H1N1)pmd09 influenza viruses may evolve to overcome the fitness cost of bearing the H275Y NA substitution, which could result in widespread oseltamivir resistance.

Structure-based design of stabilized recombinant influenza neuraminidase tetramers

Influenza virus neuraminidase (NA) is a major antiviral drug target and has recently reemerged as a key target of antibody-mediated protective immunity. Here we show that recombinant NAs across non-bat subtypes adopt various tetrameric conformations, including an "open" state that may help explain poorly understood variations in NA stability across viral strains and subtypes. We use homology-directed protein design to uncover the structural principles underlying these distinct tetrameric conformations and stabilize multiple recombinant NAs in the "closed" state, yielding two near-atomic resolution structures of NA by cryo-EM. In addition to enhancing thermal stability, conformational stabilization improves affinity to protective antibodies elicited by viral infection, including antibodies targeting a quaternary epitope and the broadly conserved catalytic site. Stabilized NAs can also be integrated into viruses without affecting fitness. Our findings provide a deeper understanding of NA structure, stability, and antigenicity, and establish design strategies for reinforcing the conformational integrity of recombinant NA proteins.

Identification of a permissive secondary mutation that restores the enzymatic activity of oseltamivir resistance mutation H275Y

Many oseltamivir resistance mutations exhibit fitness defects in the absence of drug pressure that hinders their propagation in hosts. Secondary permissive mutations can rescue fitness defects and facilitate the segregation of resistance mutations in viral populations. Previous studies have identified a panel of permissive or compensatory mutations in neuraminidase (NA) that restore the growth defect of the predominant oseltamivir resistance mutation (H275Y) in H1N1 influenza A. In prior work, we identified a hyperactive mutation (Y276F) that increased NA activity by approximately 70%. While Y276F had not been previously identified as a permissive mutation, we hypothesized that Y276F may counteract the defects caused by H275Y by buffering its reduced NA expression and enzyme activity. In this study we measured the relative fitness, NA activity, and surface expression, as well as sensitivity to oseltamivir, for several oseltamivir resistance mutations including H275Y in the wildtype or Y276F genetic background. Our results demonstrate that Y276F selectively rescues the fitness defect of H275Y by restoring its NA surface expression and enzymatic activity, elucidating the local compensatory structural impacts of Y276F on the adjacent H275Y. Importance The potential for influenza A virus (IAV) to cause pandemics makes understanding evolutionary mechanisms that impact drug resistance critical for developing surveillance and treatment strategies. Oseltamivir is the most widely used therapeutic strategy to treat IAV infections, but mutations in IAV can lead to drug resistance. The main oseltamivir resistance mutation, H275Y, occurs in the neuraminidase (NA) protein of IAV and reduces drug binding as well as NA function. Here, we identify a new helper mutation, Y276F that can rescue the functional defects of H275Y and contribute to the evolution of drug resistance in IAV.

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.

Second sialic acid-binding site of influenza A virus neuraminidase: binding receptors for efficient release

Influenza A viruses (IAVs) are a major cause of human respiratory tract infections and cause significant disease and mortality. Human IAVs originate from animal viruses that breached the host species barrier. IAV particles contain sialoglycan receptor-binding hemagglutinin (HA) and receptor-destroying neuraminidase (NA) in their envelope. When IAV crosses the species barrier, the functional balance between HA and NA needs to be adjusted to the sialoglycan repertoire of the novel host species. Relatively little is known about the role of NA in host adaptation in contrast to the extensively studied HA. NA prevents virion aggregation and facilitates release of (newly assembled) virions from cell surfaces and from decoy receptors abundantly present in mucus and cell glycocalyx. In addition to a highly conserved catalytic site, NA carries a second sialic acid-binding site (2SBS). The 2SBS preferentially binds α2,3-linked sialic acids and enhances activity of the neighboring catalytic site by bringing/keeping multivalent substrates in close contact with this site. In this way, the 2SBS contributes to the HA-NA balance of virus particles and affects virus replication. The 2SBS is highly conserved in all NA subtypes of avian IAVs, with some notable exceptions associated with changes in the receptor-binding specificity of HA and host tropism. Conservation of the 2SBS is invariably lost in human (pandemic) viruses and in several other viruses adapted to mammalian host species. Preservation or loss of the 2SBS is likely to be an important factor of the viral host range.

The mechanism of action of T-705 as a unique delayed chain terminator on influenza viral polymerase transcription

Favipiravir (T-705) has been developed as a potent anti-influenza drug and exhibited a strong inhibition effect against a broad spectrum of RNA viruses. Its active form, ribofuranosyl-triphosphate (T-705-RTP), functions as a competitive substrate for the RNA-dependent RNA polymerase (RdRp) of the influenza A virus (IAV). However, the exact inhibitory mechanisms of T-705 remain elusive and subject to a long-standing debate. Although T-705 has been proposed to inhibit transcription by acting as a chain terminator, it is also paradoxically suggested to be a mutagen towards IAV RdRp by inducing mutations due to its ambiguous base pairing of C and U. Here, we combined biochemical assay with molecular dynamics (MD) simulations to elucidate the molecular mechanism underlying the inhibitory functions exerted by T-705 in IAV RdRp. Our in vitro transcription assay illustrated that IAV RdRp could recognize T-705 as a purine analogue and incorporate it into the nascent RNA strand. Incorporating a single T-705 is incapable of inhibiting transcription as extra natural nucleotides can be progressively added. However, when two consecutive T-705 are incorporated, viral transcription is completely terminated. MD simulations reveal that the sequential appearance of two T-705 in the nascent strand destabilizes the active site and disrupts the base stacking of the nascent RNA. Altogether, our results provide a plausible explanation for the inhibitory roles of T-705 targeting IAV RdRp by integrating the computational and experimental methods. Our study also offers a comprehensive platform to investigate the inhibition effect of antivirals and a novel explanation for the designing of anti-flu drugs.

Interferon Lambda Delays the Emergence of Influenza Virus Resistance to Oseltamivir

Influenza viruses are a leading cause of morbidity and mortality worldwide. These air-borne pathogens are able to cross the species barrier, leading to regular seasonal epidemics and sporadic pandemics. Influenza viruses also possess a high genetic variability, which allows for the acquisition of resistance mutations to antivirals. Combination therapies with two or more drugs targeting different mechanisms of viral replication have been considered an advantageous option to not only enhance the effectiveness of the individual treatments, but also reduce the likelihood of resistance emergence. Using an in vitro infection model, we assessed the barrier to viral resistance of a combination therapy with the neuraminidase inhibitor oseltamivir and human interferon lambda against the pandemic H1N1 A/Netherlands/602/2009 (H1N1pdm09) virus. We serially passaged the virus in a cell line derived from human bronchial epithelial cells in the presence or absence of increasing concentrations of oseltamivir alone or oseltamivir plus interferon lambda. While the treatment with oseltamivir alone quickly induced the emergence of antiviral resistance through a single mutation in the neuraminidase gene, the co-administration of interferon lambda delayed the emergence of drug-resistant influenza virus variants. Our results suggest a possible clinical application of interferon lambda in combination with oseltamivir to treat influenza.

Secondary substitutions in the hemagglutinin and neuraminidase genes associated with neuraminidase inhibitor resistance are rare in the Influenza Resistance Information Study (IRIS)

Amino acid substitutions in influenza virus neuraminidase (NA) that cause resistance to neuraminidase inhibitors (NAI) generally result in virus attenuation. However, influenza viruses may acquire secondary substitutions in the NA and hemagglutinin (HA) proteins that can restore viral fitness. To assess to which extent this happens, the emergence of NAI resistance substitutions and secondary - potentially compensatory - substitutions was quantified in influenza viruses of immunocompetent individuals included in the Influenza Resistance Information Study (IRIS; NCT00884117). Known resistance substitutions were detected by mutation specific RT-PCR in viruses of 57 of 1803 (3.2%) oseltamivir-treated individuals, including 39 individuals infected with A/H1N1pdm09 [H275Y] virus and 18 with A/H3N2 [R292K] virus. A total of fifteen and ten other amino acid substitutions were acquired in HA and NA respectively, of A/H1N1pdm09, A/H3N2 and influenza B viruses upon treatment with oseltamivir but none of these was associated with resistance to oseltamivir. All cultured viruses with the known resistance substitutions H275Y or R292K showed reduced susceptibility to oseltamivir in the NA-star assay. Upon next-generation sequencing, the vast majority of NAI resistant A/H1N1pdm09 and A/H3N2 viruses had no resistance-associated secondary substitutions at high frequency. Only in two A/H1N1pdm09 [H275Y] viruses, the potentially compensatory secondary substitutions HA-D52N and NA-R152K were detected. We conclude that the emergence of secondary substitutions that may restore viral fitness upon the emergence of known influenza virus NAI resistance substitutions was a rare event in this immunocompetent population.

Functional neuraminidase inhibitor resistance motifs in avian influenza A(H5Nx) viruses

Neuraminidase inhibitors (NAIs) are antiviral agents recommended worldwide to treat or prevent influenza virus infections in humans. Past influenza virus pandemics seeded by zoonotic infection by avian influenza viruses (AIV) as well as the increasing number of human infections with AIV have shown the importance of having information about resistance to NAIs by avian NAs that could cross the species barrier. In this study we introduced four NAI resistance-associated mutations (N2 numbering) previously found in human infections into the NA of three current AIV subtypes of the H5Nx genotype that threaten the poultry industry and human health: highly pathogenic H5N8, H5N6 and H5N2. Using the established MUNANA assay we showed that a R292K substitution in H5N6 and H5N2 viruses significantly reduced susceptibility to three licenced NAIs: oseltamivir, zanamivir and peramivir. In contrast the mutations E119V, H274Y and N294S had more variable effects with NAI susceptibility being drug- and strain-specific. We measured the replicative fitness of NAI resistant H5N6 viruses and found that they replicated to comparable or significantly higher titres in primary chicken cells and in embryonated hens' eggs as compared to wild type - despite the NA activity of the viral neuraminidase proteins being reduced. The R292K and N294S drug resistant H5N6 viruses had single amino acid substitutions in their haemagglutinin (HA): Y98F and A189T, respectively (H3 numbering) which reduced receptor binding properties possibly balancing the reduced NA activity seen. Our results demonstrate that the H5Nx viruses can support drug resistance mutations that confer reduced susceptibility to licenced NAIs and that these H5N6 viruses did not show diminished replicative fitness in avian cell cultures. Our results support the requirement for on-going surveillance of these strains in bird populations to include motifs associated with human drug resistance.

Differential Viral-Host Immune Interactions Associated with Oseltamivir-Resistant H275Y and Wild-Type H1N1 A(pdm09) Influenza Virus Pathogenicity

Oseltamivir is a common therapy against influenza A virus (IAV) infections. The acquisition of oseltamivir resistance (OR) mutations, such as H275Y, hampers viral fitness. However, OR H1N1 viruses have demonstrated the ability to spread throughout different populations. The objective of this work was to compare the fitness of two strains of OR (R6 and R7) containing the H275Y mutation, and a wild-type (F) pandemic influenza A (H1N1) 2009 (pdm09) virus both in vitro and in vivo in mice and to select one OR strain for a comparison with F in ferrets. R6 showed faster replication and pathogenicity than R7 in vitro and in mice. Subsequently, R6 was selected for the fitness comparison with the F strain in ferrets. Ferrets infected with the F virus showed more severe clinical signs, histopathological lung lesions, and viral quantification when compared to OR R6-infected animals. More importantly, differential viral kinetics correlated with differential pro-inflammatory host immune responses in the lungs of infected ferrets, where OR-infected animals developed a protective higher expression of type I IFN and Retinoid acid Inducible Gene I (RIG-I) genes early after infection, resulting in the development of milder disease. These results suggest the presence of early specific viral-host immune interactions relevant in the development of influenza-associated lung pathology.

Genetic variations on 31 and 450 residues of influenza A nucleoprotein affect viral replication and translation

Background

Influenza A viruses cause epidemics/severe pandemics that pose a great global health threat. Among eight viral RNA segments, the multiple functions of nucleoprotein (NP) play important roles in viral replication and transcription.

Methods

To understand how NP contributes to the virus evolution, we analyzed the NP gene of H3N2 viruses in Taiwan and 14,220 NP sequences collected from Influenza Research Database. The identified genetic variations were further analyzed by mini-genome assay, virus growth assay, viral RNA and protein expression as well as ferret model to analyze their impacts on viral replication properties.

Results

The NP genetic analysis by Taiwan and global sequences showed similar evolution pattern that the NP backbones changed through time accompanied with specific residue substitutions from 1999 to 2018. Other than the conserved residues, fifteen sporadic substitutions were observed in which the 31R, 377G and 450S showed higher frequency. We found 31R and 450S decreased polymerase activity while the dominant residues (31 K and 450G) had higher activity. The 31 K and 450G showed better viral translation and replication in vitro and in vivo.

Conclusions

These findings indicated variations identified in evolution have roles in modulating viral replication in vitro and in vivo. This study demonstrates that the interaction between variations of NP during virus evolution deserves future attention.

Mutations in Influenza A Virus Neuraminidase and Hemagglutinin Confer Resistance against a Broadly Neutralizing Hemagglutinin Stem Antibody

Influenza A virus (IAV), a major cause of human morbidity and mortality, continuously evolves in response to selective pressures. Stem-directed, broadly neutralizing antibodies (sBnAbs) targeting the influenza virus hemagglutinin (HA) are a promising therapeutic strategy, but neutralization escape mutants can develop. We used an integrated approach combining viral passaging, deep sequencing, and protein structural analyses to define escape mutations and mechanisms of neutralization escape in vitro for the F10 sBnAb. IAV was propagated with escalating concentrations of F10 over serial passages in cultured cells to select for escape mutations. Viral sequence analysis revealed three mutations in HA and one in neuraminidase (NA). Introduction of these specific mutations into IAV through reverse genetics confirmed their roles in resistance to F10. Structural analyses revealed that the selected HA mutations (S123G, N460S, and N203V) are away from the F10 epitope but may indirectly impact influenza virus receptor binding, endosomal fusion, or budding. The NA mutation E329K, which was previously identified to be associated with antibody escape, affects the active site of NA, highlighting the importance of the balance between HA and NA function for viral survival. Thus, whole-genome population sequencing enables the identification of viral resistance mutations responding to antibody-induced selective pressure.IMPORTANCE Influenza A virus is a public health threat for which currently available vaccines are not always effective. Broadly neutralizing antibodies that bind to the highly conserved stem region of the influenza virus hemagglutinin (HA) can neutralize many influenza virus strains. To understand how influenza virus can become resistant or escape such antibodies, we propagated influenza A virus in vitro with escalating concentrations of antibody and analyzed viral populations by whole-genome sequencing. We identified HA mutations near and distal to the antibody binding epitope that conferred resistance to antibody neutralization. Additionally, we identified a neuraminidase (NA) mutation that allowed the virus to grow in the presence of high concentrations of the antibody. Virus carrying dual mutations in HA and NA also grew under high antibody concentrations. We show that NA mutations mediate the escape of neutralization by antibodies against HA, highlighting the importance of a balance between HA and NA for optimal virus function.

Rapid and simple detection of Tamiflu-resistant influenza virus: Development of oseltamivir derivative-based lateral flow biosensor for point-of-care (POC) diagnostics

We have developed a novel oseltamivir derivative (oseltamivir hexylthiol; OHT) that exhibits a higher binding affinity for Tamiflu-resistant virus (Tamiflu resistance) than for the wild-type virus (Tamiflu-susceptible virus; WT) as an antibody. First, OHT-modified gold nanoparticles (OHT-GNPs) are used in a simple colorimetric assay as nanoprobes for the Tamiflu-resistant virus. In the presence of Tamiflu-resistant virus, they show a colorimetric change from deep red to purple because of the OHT-GNP aggregation driven by strong interactions between OHT and neuraminidase (NA) on the surface of the Tamiflu-resistance. Moreover, the color gradually turns purple as the concentration of the Tamiflu-resistant virus increases, allowing the determination of the presence of the virus with the naked eye. Furthermore, an OHT-based lateral flow assay (LFA) has been developed as a rapid and easy detection device for Tamiflu resistance. It shows detection specificity for various virus concentrations of Tamiflu-resistant virus even for the mixture of WT and Tamiflu-resistant viruses, where the limit of detection (LOD) is 5 × 10² ~ 10³ PFU per test (=1 × 10⁴ PFU/mL). It has been confirmed that this platform can provide accurate information on whether a virus exhibits Tamiflu resistance, thus supporting the selection of appropriate treatments using point-of-care (POC) diagnostics.

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.

Predictable Phenotypes of Antibiotic Resistance Mutations

Antibiotic-resistant bacteria represent a major threat to our ability to treat bacterial infections. Two factors that determine the evolutionary success of antibiotic resistance mutations are their impact on resistance level and the fitness cost. Recent studies suggest that resistance mutations commonly show epistatic interactions, which would complicate predictions of their stability in bacterial populations. We analyzed 13 different chromosomal resistance mutations and 10 host strains of Salmonella enterica and Escherichia coli to address two main questions. (i) Are there epistatic interactions between different chromosomal resistance mutations? (ii) How does the strain background and genetic distance influence the effect of chromosomal resistance mutations on resistance and fitness? Our results show that the effects of combined resistance mutations on resistance and fitness are largely predictable and that epistasis remains rare even when up to four mutations were combined. Furthermore, a majority of the mutations, especially target alteration mutations, demonstrate strain-independent phenotypes across different species. This study extends our understanding of epistasis among resistance mutations and shows that interactions between different resistance mutations are often predictable from the characteristics of the individual mutations.

The spread of antibiotic-resistant bacteria imposes an urgent threat to public health. The ability to forecast the evolutionary success of resistant mutants would help to combat dissemination of antibiotic resistance. Previous studies have shown that the phenotypic effects (fitness and resistance level) of resistance mutations can vary substantially depending on the genetic context in which they occur. We conducted a broad screen using many different resistance mutations and host strains to identify potential epistatic interactions between various types of resistance mutations and to determine the effect of strain background on resistance phenotypes. Combinations of several different mutations showed a large amount of phenotypic predictability, and the majority of the mutations displayed strain-independent phenotypes. However, we also identified a few outliers from these patterns, illustrating that the choice of host organism can be critically important when studying antibiotic resistance mutations.

Comparative studies on variability, phylogenesis, and correlated mutations of neuraminidases from influenza virus type A

Neuraminidase (NA) is an important protein for the replication cycle of influenza A viruses. NA is an enzyme that cleaves the sialic acid receptors; this process plays a significant role in viral life cycle. Blocking NA with a specific inhibitor is an effective way to treat the flu. However, some strains show resistance to current drugs. Therefore, NA is the focus for the intense research for new antiviral drugs and also for the explanation of the functions of new mutations. This research focuses on determining the profile of variability and phylogenetic analysis and finding the correlated mutations within a set of 149 sequences of NA belonging to various strains of influenza A virus. In this study, we have used the original programs (Corm, Consensus Constructor, and SSSSg) and also other bioinformatics software. NA proteins are characterized by various levels of variability in different regions, which was presented in detail with the aid of ConSurf. The use of four independent methods to create the phylogenetic trees gave some new data on the evolutionary relationship within the NA family proteins. The search for correlated mutations shows several potentially important correlated positions that were not reported previously to be significant. The use of such an approach can be potentially important and gives new information regarding NA proteins of influenza A virus.

Influenza Virus: Small Molecule Therapeutics and Mechanisms of Antiviral Resistance

BACKGROUND

Influenza viruses cause severe upper respiratory illness in children and the elderly during seasonal epidemics. Influenza viruses from zoonotic reservoirs can also cause pandemics with significant loss of life in all age groups. Although vaccination is one of the most effective methods to protect against seasonal epidemics, seasonal vaccines vary in efficacy, can be ineffective in the elderly population, and do not provide protection against novel strains. Small molecule therapeutics are a critical part of our antiviral strategies to control influenza virus epidemics and pandemics as well as to ameliorate disease in elderly and immunocompromised individuals.

OBJECTIVE

This review aims to summarize the existing antiviral strategies for combating influenza viruses, the mechanisms of antiviral resistance for available drugs, and novel therapeutics currently in development.

METHODS

We systematically evaluated and synthesized the published scientific literature for mechanistic detail into therapeutic strategies against influenza viruses.

RESULTS

Current IAV strains have developed resistance to neuraminidase inhibitors and nearly complete resistance to M2 ion channel inhibitors, exacerbated by sub-therapeutic dosing used for treatment and chemoprophylaxis. New tactics include novel therapeutics targeting host components and combination therapy, which show potential for fighting influenza virus disease while minimizing viral resistance.

CONCLUSION

Antiviral drugs are crucial for controlling influenza virus disease burden, but their efficacy is limited by human misuse and the capacity of influenza viruses to circumvent antiviral barriers. To relieve the public health hardship of influenza virus, emerging therapies must be selected for their capacity to impede not only influenza virus disease, but also the development of antiviral resistance.

In vivo mallard experiments indicate that zanamivir has less potential for environmental influenza A virus resistance development than oseltamivir

Neuraminidase inhibitors are a cornerstone of influenza pandemic preparedness before vaccines can be mass-produced and thus a neuraminidase inhibitor-resistant pandemic is a serious threat to public health. Earlier work has demonstrated the potential for development and persistence of oseltamivir resistance in influenza A viruses exposed to environmentally relevant water concentrations of the drug when infecting mallards, the natural influenza reservoir that serves as the genetic base for human pandemics. As zanamivir is the major second-line neuraminidase inhibitor treatment, this study aimed to assess the potential for development and persistence of zanamivir resistance in an in vivo mallard model; especially important as zanamivir will probably be increasingly used. Our results indicate less potential for development and persistence of resistance due to zanamivir than oseltamivir in an environmental setting. This conclusion is based on: (1) the lower increase in zanamivir IC₅₀ conferred by the mutations caused by zanamivir exposure (2-17-fold); (2) the higher zanamivir water concentration needed to induce resistance (at least 10 µg l^-1); (3) the lack of zanamivir resistance persistence without drug pressure; and (4) the multiple resistance-related substitutions seen during zanamivir exposure (V116A, A138V, R152K, T157I and D199G) suggesting lack of one straight-forward evolutionary path to resistance. Our study also adds further evidence regarding the stability of the oseltamivir-induced substitution H275Y without drug pressure, and demonstrates the ability of a H275Y-carrying virus to acquire secondary mutations, further boosting oseltamivir resistance when exposed to zanamivir. Similar studies using influenza A viruses of the N2-phylogenetic group of neuraminidases are recommended.

Identification of neuraminidase inhibitors against dual H274Y/I222R mutant strains

Influenza is an annual seasonal epidemic that has continually drawn public attentions, due to the potential death toll and drug resistance. Neuraminidase, which is essential for the spread of influenza virus, has been regarded as a valid target for the treatment of influenza infection. Although neuraminidase drugs have been developed, they are susceptible to drug-resistant mutations in the sialic-binding site. In this study, we established computational models (site-moiety maps) of H1N1 and H5N1 to determine properties of the 150-cavity, which is adjacent to the drug-binding site. The models reveal that hydrogen-bonding interactions with residues R118, D151, and R156 and van der Waals interactions with residues Q136, D151, and T439 are important for identifying 150-cavitiy inhibitors. Based on the models, we discovered three new inhibitors with IC₅₀ values <10 μM that occupies both the 150-cavity and sialic sites. The experimental results identified inhibitors with similar activities against both wild-type and dual H274Y/I222R mutant neuraminidases and showed little cytotoxic effects. Furthermore, we identified three new inhibitors situated at the sialic-binding site with inhibitory effects for normal neuraminidase, but lowered effects for mutant strains. The results suggest that the new inhibitors can be used as a starting point to combat drug-resistant strains.

Drug resistance in influenza A virus: the epidemiology and management

Influenza A virus (IAV) is the sole cause of the unpredictable influenza pandemics and deadly zoonotic outbreaks and constitutes at least half of the cause of regular annual influenza epidemics in humans. Two classes of anti-IAV drugs, adamantanes and neuraminidase (NA) inhibitors (NAIs) targeting the viral components M2 ion channel and NA, respectively, have been approved to treat IAV infections. However, IAV rapidly acquired resistance against both classes of drugs by mutating these viral components. The adamantane-resistant IAV has established itself in nature, and a majority of the IAV subtypes, especially the most common H1N1 and H3N2, circulating globally are resistant to adamantanes. Consequently, adamantanes have become practically obsolete as anti-IAV drugs. Similarly, up to 100% of the globally circulating IAV H1N1 subtypes were resistant to oseltamivir, the most commonly used NAI, until 2009. However, the 2009 pandemic IAV H1N1 subtype, which was sensitive to NAIs and has now become one of the dominant seasonal influenza virus strains, has replaced the pre-2009 oseltamivir-resistant H1N1 variants. This review traces the epidemiology of both adamantane- and NAI-resistant IAV subtypes since the approval of these drugs and highlights the susceptibility status of currently circulating IAV subtypes to NAIs. Further, it provides an overview of currently and soon to be available control measures to manage current and emerging drug-resistant IAV. Finally, this review outlines the research directions that should be undertaken to manage the circulation of IAV in intermediate hosts and develop effective and alternative anti-IAV therapies.

Evolution of the neuraminidase gene of seasonal influenza A and B viruses in Thailand between 2010 and 2015

The neuraminidase inhibitors (NAIs) oseltamivir and zanamivir are commonly used for the treatment and control of influenza A and B virus infection. However, the emergence of new influenza virus strains with reduced susceptibility to NAIs may appear with the use of these antivirals or even naturally. We therefore screened the neuraminidase (NA) sequences of seasonal influenza virus A(H1N1), A(H1N1)pdm09, A(H3N2), and influenza B virus strains identified in Thailand for the presence of substitutions previously reported to reduce susceptibility to NAIs. We initially examined oseltamivir resistance (characterized by the H275Y mutation in the NA gene) in 485 A(H1N1)pdm09 strains circulating in Thailand and found that 0.82% (4/485) had this substitution. To further evaluate the evolution of the NA gene, we also randomly selected 98 A(H1N1)pdm09, 158 A(H3N2), and 69 influenza B virus strains for NA gene amplification and sequencing, which revealed various amino acid mutations in the active site of the NA protein previously shown to be associated with reduced susceptibility to NAIs. Phylogenetic analysis of the influenza virus strains from this study and elsewhere around the world, together with the estimations of nucleotide substitution rates and selection pressure, and the predictions of B-cell epitopes and N-linked glycosylation sites all provided evidence for the ongoing evolution of NA. The overall rates of NA evolution for influenza A viruses were higher than for influenza B virus at the nucleotide level, although influenza B virus possessed more genealogical diversity than that of influenza A viruses. The continual surveillance of the antigenic changes associated with the NA protein will not only contribute to the influenza virus database but may also provide a better understanding of selection pressure exerted by antiviral use.

The inhibition of H1N1 influenza virus-induced apoptosis by silver nanoparticles functionalized with zanamivir

As one of the most effective drugs for influenza virus infection, clinical application of zanamivir is restricted with the emergence of resistant influenza virus.

Widespread Historical Contingency in Influenza Viruses

In systems biology and genomics, epistasis characterizes the impact that a substitution at a particular location in a genome can have on a substitution at another location. This phenomenon is often implicated in the evolution of drug resistance or to explain why particular "disease-causing" mutations do not have the same outcome in all individuals. Hence, uncovering these mutations and their locations in a genome is a central question in biology. However, epistasis is notoriously difficult to uncover, especially in fast-evolving organisms. Here, we present a novel statistical approach that replies on a model developed in ecology and that we adapt to analyze genetic data in fast-evolving systems such as the influenza A virus. We validate the approach using a two-pronged strategy: extensive simulations demonstrate a low-to-moderate sensitivity with excellent specificity and precision, while analyses of experimentally validated data recover known interactions, including in a eukaryotic system. We further evaluate the ability of our approach to detect correlated evolution during antigenic shifts or at the emergence of drug resistance. We show that in all cases, correlated evolution is prevalent in influenza A viruses, involving many pairs of sites linked together in chains; a hallmark of historical contingency. Strikingly, interacting sites are separated by large physical distances, which entails either long-range conformational changes or functional tradeoffs, for which we find support with the emergence of drug resistance. Our work paves a new way for the unbiased detection of epistasis in a wide range of organisms by performing whole-genome scans.

Risk of resistant avian influenza A virus in wild waterfowl as a result of environmental release of oseltamivir

Oseltamivir is the best available anti-influenza drug and has therefore been stockpiled worldwide in large quantities as part of influenza pandemic preparedness planning. The active metabolite oseltamivir carboxylate (OC) is stable and is not removed by conventional sewage treatment. Active OC has been detected in river water at concentrations up to 0.86 µg/L. Although the natural reservoir hosts of influenza A virus (IAV) are wild waterfowl that reside in aquatic environments, the ecologic risks associated with environmental OC release and its potential to generate resistant viral variants among wild birds has largely been unknown. However, in recent years a number of in vivo mallard (Anas platyrhynchos) studies have been conducted regarding the potential of avian IAVs to become resistant to OC in natural reservoir birds if these are drug exposed. Development of resistance to OC was observed both in Group 1 (N1) and Group 2 (N2, N9) neuraminidase subtypes, when infected ducks were exposed to OC at concentrations between 0.95 and 12 µg/L in their water. All resistant variants maintained replication and transmission between ducks during drug exposure. In an A(H1N1)/H274Y virus, the OC resistance mutation persisted without selective drug pressure, demonstrating the potential of an IAV with a permissive genetic background to acquire and maintain OC resistance, potentially allowing circulation of the resistant variant among wild birds. The experimental studies have improved the appreciation of the risks associated with the environmental release of OC related to resistance development of avian IAVs among wild birds. Combined with knowledge of efficient methods for improved sewage treatment, the observations warrant implementation of novel efficient wastewater treatment methods, rational use of anti-influenza drugs, and improved surveillance of IAV resistance in wild birds.

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.

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.

Oseltamivir Resistance in Influenza A(H6N2) Caused by an R292K Substitution in Neuraminidase Is Not Maintained in Mallards without Drug Pressure

Background

Wild waterfowl is the natural reservoir of influenza A virus (IAV); hosted viruses are very variable and provide a source for genetic segments which can reassort with poultry or mammalian adapted IAVs to generate novel species crossing viruses. Additionally, wild waterfowl act as a reservoir for highly pathogenic IAVs. Exposure of wild birds to the antiviral drug oseltamivir may occur in the environment as its active metabolite can be released from sewage treatment plants to river water. Resistance to oseltamivir, or to other neuraminidase inhibitors (NAIs), in IAVs of wild waterfowl has not been extensively studied.

Aim and Methods

In a previous in vivo Mallard experiment, an influenza A(H6N2) virus developed oseltamivir resistance by the R292K substitution in the neuraminidase (NA), when the birds were exposed to oseltamivir. In this study we tested if the resistance could be maintained in Mallards without drug exposure. Three variants of resistant H6N2/R292K virus were each propagated during 17 days in five successive pairs of naïve Mallards, while oseltamivir exposure was decreased and removed. Daily fecal samples were analyzed for viral presence, genotype and phenotype.

Results and Conclusion

Within three days without drug exposure no resistant viruses could be detected by NA sequencing, which was confirmed by functional NAI sensitivity testing. We conclude that this resistant N2 virus could not compete in fitness with wild type subpopulations without oseltamivir drug pressure, and thus has no potential to circulate among wild birds. The results of this study contrast to previous observations of drug induced resistance in an avian H1N1 virus, which was maintained also without drug exposure in Mallards. Experimental observations on persistence of NAI resistance in avian IAVs resemble NAI resistance seen in human IAVs, in which resistant N2 subtypes do not circulate, while N1 subtypes with permissive mutations can circulate without drug pressure. We speculate that the phylogenetic group N1 NAs may easier compensate for NAI resistance than group N2 NAs, though further studies are needed to confirm such conclusions.

Clinical Implications of Antiviral Resistance in Influenza

Influenza is a major cause of severe respiratory infections leading to excessive hospitalizations and deaths globally; annual epidemics, pandemics, and sporadic/endemic avian virus infections occur as a result of rapid, continuous evolution of influenza viruses. Emergence of antiviral resistance is of great clinical and public health concern. Currently available antiviral treatments include four neuraminidase inhibitors (oseltamivir, zanamivir, peramivir, laninamivir), M2-inibitors (amantadine, rimantadine), and a polymerase inhibitor (favipiravir). In this review, we focus on resistance issues related to the use of neuraminidase inhibitors (NAIs). Data on primary resistance, as well as secondary resistance related to NAI exposure will be presented. Their clinical implications, detection, and novel therapeutic options undergoing clinical trials are discussed.

Coordinated Evolution of Influenza A Surface Proteins

The surface proteins hemagglutinin (HA) and neuraminidase (NA) of human influenza A virus evolve under selection pressures to escape adaptive immune responses and antiviral drug treatments. In addition to these external selection pressures, some mutations in HA are known to affect the adaptive landscape of NA, and vice versa, because these two proteins are physiologically interlinked. However, the extent to which evolution of one protein affects the evolution of the other one is unknown. Here we develop a novel phylogenetic method for detecting the signatures of such genetic interactions between mutations in different genes – that is, inter-gene epistasis. Using this method, we show that influenza surface proteins evolve in a coordinated way, with mutations in HA affecting subsequent spread of mutations in NA and vice versa, at many sites. Of particular interest is our finding that the oseltamivir-resistance mutations in NA in subtype H1N1 were likely facilitated by prior mutations in HA. Our results illustrate that the adaptive landscape of a viral protein is remarkably sensitive to its genomic context and, more generally, that the evolution of any single protein must be understood within the context of the entire evolving genome.

The fitness of an organism depends on the coordinated function of many genes. Thus, how a mutation in one gene affects fitness often depends on what mutations are present in other genes. This dependence is called “genetic interaction” or “epistasis”. The prevalence and type of such interactions are not well understood. Epistasis can be inferred from time-series sequencing data when a mutation in one gene is observed to facilitate the spread of a mutation in another gene. However, the situation is much more complicated when new combinations of genes are formed by processes such as recombination or reassortment. In such cases, deducing the time and order of genetic changes is difficult. Here, we devise a method to infer pairs of mutations in different genes which closely follow one another in the presence of reassortment. We apply it to evolution of two surface proteins of influenza A virus, hemagglutinin and neuraminidase, which are important targets for the human immune system and drugs. We show that mutations in one of these proteins are often facilitated by prior mutations, or compensated by subsequent mutations, in the other protein. In particular, drug-resistance mutations in neuraminidase were likely made possible by prior mutation in hemagglutinin. Knowledge of such interactions is necessary to fully understand and predict evolution.

Expected Effect of Deleterious Mutations on Within-Host Adaptation of Pathogens

Adaptation is a common theme in both pathogen emergence, for example, in zoonotic cross-species transmission, and pathogen control, where adaptation might limit the effect of the immune response and antiviral treatment. When such evolution requires deleterious intermediate mutations, fitness ridges and valleys arise in the pathogen's fitness landscape. The effect of deleterious intermediate mutations on within-host pathogen adaptation is examined with deterministic calculations, appropriate for pathogens replicating in large populations with high error rates. The effect of deleterious intermediate mutations on pathogen adaptation is smaller than their name might suggest: when two mutations are required and each individual single mutation is fully deleterious, the pathogen can jump across the fitness valley by obtaining two mutations at once, leading to a proportion of adapted mutants that is 20-fold lower than that in the situation where the fitness of all mutants is neutral. The negative effects of deleterious intermediates are typically substantially smaller and outweighed by the fitness advantages of the adapted mutant. Moreover, requiring a specific mutation order has a substantially smaller effect on pathogen adaptation than the effect of all intermediates being deleterious. These results can be rationalized when the number of routes of mutation available to the pathogen is calculated, providing a simple approach to estimate the effect of deleterious mutations. The calculations discussed here are applicable when the effect of deleterious mutations on the within-host adaptation of pathogens is assessed, for example, in the context of zoonotic emergence, antigenic escape, and drug resistance.

IMPORTANCE

Adaptation is critical for pathogens after zoonotic transmission into a new host species or to achieve antigenic immune escape and drug resistance. Using a deterministic approach, the effects of deleterious intermediate mutations on pathogen adaptation were calculated while avoiding commonly made simplifications that do not apply to large pathogen populations replicating with high mutation rates. Perhaps unexpectedly, pathogen adaptation does not halt when the intermediate mutations are fully deleterious. The negative effects of deleterious mutations are substantially outweighed by the fitness gains of adaptation. To gain an understanding of the effect of deleterious mutations on pathogen adaptation, a simple approach that counts the number of routes available to the pathogen with and without deleterious intermediate mutations is introduced. This methodology enables a straightforward calculation of the proportion of the pathogen population that will cross a fitness valley or traverse a fitness ridge, without reverting to more complicated mathematical models.

Influenza A(H7N9) virus acquires resistance-related neuraminidase I222T substitution when infected mallards are exposed to low levels of oseltamivir in water

Influenza A virus (IAV) has its natural reservoir in wild waterfowl, and new human IAVs often contain gene segments originating from avian IAVs. Treatment options for severe human influenza are principally restricted to neuraminidase inhibitors (NAIs), among which oseltamivir is stockpiled in preparedness for influenza pandemics. There is evolutionary pressure in the environment for resistance development to oseltamivir in avian IAVs, as the active metabolite oseltamivir carboxylate (OC) passes largely undegraded through sewage treatment to river water where waterfowl reside. In an in vivo mallard (Anas platyrhynchos) model, we tested if low-pathogenic avian influenza A(H7N9) virus might become resistant if the host was exposed to low levels of OC. Ducks were experimentally infected, and OC was added to their water, after which infection and transmission were maintained by successive introductions of uninfected birds. Daily fecal samples were tested for IAV excretion, genotype, and phenotype. Following mallard exposure to 2.5 μg/liter OC, the resistance-related neuraminidase (NA) I222T substitution, was detected within 2 days during the first passage and was found in all viruses sequenced from subsequently introduced ducks. The substitution generated 8-fold and 2.4-fold increases in the 50% inhibitory concentration (IC50) for OC (P < 0.001) and zanamivir (P = 0.016), respectively. We conclude that OC exposure of IAV hosts, in the same concentration magnitude as found in the environment, may result in amino acid substitutions, leading to changed antiviral sensitivity in an IAV subtype that can be highly pathogenic to humans. Prudent use of oseltamivir and resistance surveillance of IAVs in wild birds are warranted.

Influenza A viruses of swine circulating in the United States during 2009-2014 are susceptible to neuraminidase inhibitors but show lineage-dependent resistance to adamantanes

Antiviral drug susceptibility is one of the evaluation criteria of pandemic potential posed by an influenza virus. Influenza A viruses of swine (IAV-S) can play an important role in generating novel variants, yet limited information is available on the drug resistance profiles of IAV-S circulating in the U.S. Phenotypic analysis of the IAV-S isolated in the U.S. (2009-2011) (n=105) revealed normal inhibition by the neuraminidase (NA) inhibitors (NAIs) oseltamivir, zanamivir, and peramivir. Screening NA sequences from IAV-S collected in the U.S. (1930-2014) showed 0.03% (1/3396) sequences with clinically relevant H274Y-NA substitution. Phenotypic analysis of IAV-S isolated in the U.S. (2009-2011) confirmed amantadine resistance caused by the S31N-M2 and revealed an intermediate level of resistance caused by the I27T-M2. The majority (96.7%, 589/609) of IAV-S with the I27T-M2 in the influenza database were isolated from pigs in the U.S. The frequency of amantadine-resistant markers among IAV-S in the U.S. was high (71%), and their distribution was M-lineage dependent. All IAV-S of the Eurasian avian M lineage were amantadine-resistant and possessed either a single S31N-M2 substitution (78%, 585/747) or its combination with the V27A-M2 (22%, 162/747). The I27T-M2 substitution accounted for 43% (429/993) of amantadine resistance in classic swine M lineage. Phylogenetic analysis showed that both S31N-M2 and I27T-M2 emerged stochastically but appeared to be fixed in the U.S. IAV-S population. This study defines a drug-susceptibility profile, identifies the frequency of drug-resistant markers, and establishes a phylogenetic approach for continued antiviral-susceptibility monitoring of IAV-S in the U.S.

Oseltamivir-resistant influenza A (H1N1) virus strain with an H274Y mutation in neuraminidase persists without drug pressure in infected mallards

Influenza A virus (IAV) has its natural reservoir in wild waterfowl, and emerging human IAVs often contain gene segments from avian viruses. The active drug metabolite of oseltamivir (oseltamivir carboxylate [OC]), stockpiled as Tamiflu for influenza pandemic preparedness, is not removed by conventional sewage treatment and has been detected in river water. There, it may exert evolutionary pressure on avian IAV in waterfowl, resulting in the development of resistant viral variants. A resistant avian IAV can circulate among wild birds only if resistance does not restrict viral fitness and if the resistant virus can persist without continuous drug pressure. In this in vivo mallard (Anas platyrhynchos) study, we tested whether an OC-resistant avian IAV (H1N1) strain with an H274Y mutation in the neuraminidase (NA-H274Y) could retain resistance while drug pressure was gradually removed. Successively infected mallards were exposed to decreasing levels of OC, and fecal samples were analyzed for the neuraminidase sequence and phenotypic resistance. No reversion to wild-type virus was observed during the experiment, which included 17 days of viral transmission among 10 ducks exposed to OC concentrations below resistance induction levels. We conclude that resistance in avian IAV that is induced by exposure of the natural host to OC can persist in the absence of the drug. Thus, there is a risk that human-pathogenic IAVs that evolve from IAVs circulating among wild birds may contain resistance mutations. An oseltamivir-resistant pandemic IAV would pose a substantial public health threat. Therefore, our observations underscore the need for prudent oseltamivir use, upgraded sewage treatment, and surveillance for resistant IAVs in wild birds.