Characterization of Influenza B Virus Variants with Reduced Neuraminidase Inhibitor Susceptibility

Clinical characteristics of outpatients with influenza-B-associated pneumonia and molecular evolution of influenza B virus in Beijing, China, during the 2021-2022 influenza season

We analyzed the clinical characteristics of outpatients with influenza-B-associated pneumonia during the 2021-2022 influenza season and analyzed the molecular epidemiology and evolution of influenza B virus. The presence of influenza B virus was confirmed by reverse transcription polymerase chain reaction (RT-PCR). Electronic medical records were used to collect and analyze data of outpatients. The HA and NA genes were phylogenetically analyzed using ClustalW 2.10 and MEGA 11.0. Out of 1569 outpatients who tested positive for influenza B virus, 11.7% (184/1569) developed pneumonia, and of these, 19.0% (35/184) had underlying diseases. Fever, cough, and sore throat were the most common symptoms. Among the complications, acute respiratory distress syndrome (ARDS), acute kidney injury (AKI), and shock accounted for 2.7% (5/184), 4.9% (9/184), and 1.6% (3/184), respectively. Of the outpatients, 2.7% (5/184) were admitted to the hospital, and 0.5% (1/184) of them died. All of the strains from Beijing were identified as belonging to the B/Victoria lineage. The HA and NA gene sequences of 41 influenza B viruses showed high similarity to each other, and all of them belonged to clade 1A.3. Compared with the vaccine strain B/Washington/02/2019, all of the isolates contained N150K, G181E, and S194D mutations. S194D, E195K, and K200R mutations were detected in the 190 helix of the receptor binding region of HA. Co-mutations of H122Q, A127T, P144L, N150K, G181E, S194D, and K200R in HA and D53N, N59S, and G233E in NA were detected in 78.0% (32/41) of the isolates, and 56.3% (18/32) of these were from outpatients with influenza-B-associated pneumonia. Influenza outpatients with underlying diseases were more likely to develop pneumonia. No significant differences were observed in clinical symptoms or laboratory results between outpatients with and without pneumonia, so testing for influenza virus seems to be a good choice. The observed amino acid variations suggest that current vaccines might not provide effective protection.

In vitro and in vivo characterization of SARS-CoV-2 strains resistant to nirmatrelvir

Nirmatrelvir, an oral antiviral agent that targets a SARS-CoV-2 main protease (3CLpro), is clinically useful against infection with SARS-CoV-2 including its omicron variants. Since most omicron subvariants have reduced sensitivity to many monoclonal antibody therapies, potential SARS-CoV-2 resistance to nirmatrelvir is a major public health concern. Several amino acid substitutions have been identified as being responsible for reduced susceptibility to nirmatrelvir. Among them, we selected L50F/E166V and L50F/E166A/L167F in the 3CLpro because these combinations of substitutions are unlikely to affect virus fitness. We prepared and characterized delta variants possessing Nsp5-L50F/E166V and Nsp5-L50F/E166A/L167F. Both mutant viruses showed decreased susceptibility to nirmatrelvir and their growth in VeroE6/TMPRSS2 cells was delayed. Both mutant viruses showed attenuated phenotypes in a male hamster infection model, maintained airborne transmissibility, and were outcompeted by wild-type virus in co-infection experiments in the absence of nirmatrelvir, but less so in the presence of the drug. These results suggest that viruses possessing Nsp5-L50F/E166V and Nsp5-L50F/E166A/L167F do not become dominant in nature. However, it is important to closely monitor the emergence of nirmatrelvir-resistant SARS-CoV-2 variants because resistant viruses with additional compensatory mutations could emerge, outcompete the wild-type virus, and become dominant.

Genetic characterization and whole-genome sequencing-based genetic analysis of influenza virus in Jining City during 2021-2022

Background

The influenza virus poses a significant threat to global public health due to its high mutation rate. Continuous surveillance, development of new vaccines, and public health measures are crucial in managing and mitigating the impact of influenza outbreaks.

Methods

Nasal swabs were collected from individuals with influenza-like symptoms in Jining City during 2021-2022. Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) was used to detect influenza A viruses, followed by isolation using MDCK cells. Additionally, nucleic acid detection was performed to identify influenza A H1N1, seasonal H3N2, B/Victoria, and B/Yamagata strains. Whole-genome sequencing was conducted on 24 influenza virus strains, and subsequent analyses included characterization, phylogenetic construction, mutation analysis, and assessment of nucleotide diversity.

Results

A total of 1,543 throat swab samples were collected. The study revealed the dominance of the B/Victoria influenza virus in Jining during 2021-2022. Whole-genome sequencing showed co-prevalence of B/Victoria influenza viruses in the branches of Victoria clade 1A.3a.1 and Victoria clade 1A.3a.2, with a higher incidence observed in winter and spring. Comparative analysis demonstrated lower similarity in the HA, MP, and PB2 gene segments of the 24 sequenced influenza virus strains compared to the Northern Hemisphere vaccine strain B/Washington/02/2019. Mutations were identified in all antigenic epitopes of the HA protein at R133G, N150K, and N197D, and the 17-sequence antigenic epitopes exhibited more than 4 amino acid variation sites, resulting in antigenic drift. Moreover, one sequence had a D197N mutation in the NA protein, while seven sequences had a K338R mutation in the PA protein.

Conclusion

This study highlights the predominant presence of B/Victoria influenza strain in Jining from 2021 to 2022. The analysis also identified amino acid site variations in the antigenic epitopes, contributing to antigenic drift.

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.

An influenza A(H5N8) virus isolated during an outbreak at a poultry farm in Russia in 2017 has an N294S substitution in the neuraminidase and shows reduced susceptibility to oseltamivir

This study aimed to assess the antiviral susceptibility of influenza A(H5N8) viruses isolated in Russia in 2014-2018. Genetic analysis of 57 Russian isolates with full genome sequences did not find any markers of reduced susceptibility to baloxavir. Only one strain bore an amino acid substitution associated with adamantane resistance (M2-S31N). The neuraminidase of 1 strain had an NA-N293/294S (N8/N2 numbering) substitution associated with reduced inhibition by oseltamivir and normal inhibition by zanamivir, which was confirmed phenotypically. There were no other strains with reduced inhibition by oseltamivir and zanamivir in the phenotypic analysis. In order to estimate the worldwide prevalence of influenza A(H5N8) viruses bearing genetic markers of antiviral resistance, genome sequences deposited in the GISAID database were analyzed (database access: October 2020). The M2 protein of A(H5N8) viruses from the 2.3.4.4c clade had an M2-S31N substitution associated with reduced susceptibility to adamantanes. On the contrary, the majority (94%) of viruses from the 2.3.4.4b clade had the M2-S31 genotype. Fewer than 1% of analyzed viruses had amino acid substitutions associated with reduced susceptibility to baloxavir (PA-E199G, PA-E199E/G) or reduced or highly reduced inhibition by neuraminidase inhibitors (NA-R150/152K, NA-I221/222M, NA-I221/222I/M, NA-I221/222V, NA-I115/117V, NA-G145/147R, NA-R291/292R/K). An NA-N293/294S substitution was not present in sequences from the GISAID database. To the best of our knowledge, influenza A(H5N8) viruses with reduced inhibition by oseltamivir bearing an NA-N293/294S substitution have not been previously reported in epidemiological surveillance studies.

Targeted Profiling of Immunological Genes during Norovirus Replication in Human Intestinal Enteroids

Norovirus is the leading cause of acute gastroenteritis worldwide. The pathogenesis of norovirus and the induced immune response remain poorly understood due to the lack of a robust virus culture system. The monolayers of two secretor-positive Chinese human intestinal enteroid (HIE) lines were challenged with two norovirus pandemic GII.4 Sydney strains. Norovirus RNA replication in supernatants and cell lysates were quantified by RT-qPCR. RNA expression levels of immune-related genes were profiled using PCR arrays. The secreted protein levels of shortlisted upregulated genes were measured in supernatants using analyte-specific enzyme-linked immunosorbent assay (ELISA). Productive norovirus replications were achieved in three (75%) out of four inoculations. The two most upregulated immune-related genes were CXCL10 (93-folds) and IFI44L (580-folds). Gene expressions of CXCL10 and IFI44L were positively correlated with the level of norovirus RNA replication (CXCL10: Spearman’s r = 0.779, p < 0.05; IFI44L: r = 0.881, p < 0.01). The higher level of secreted CXCL10 and IFI44L proteins confirmed their elevated gene expression. The two genes have been reported to be upregulated in norovirus volunteer challenges and natural human infections by other viruses. Our data suggested that HIE could mimic the innate immune response elicited in natural norovirus infection and, therefore, could serve as an experimental model for future virus-host interaction and antiviral studies.

Quantifying mechanistic traits of influenza viral dynamics using in vitro data

When analysing in vitro data, growth kinetics of influenza virus strains are often compared by computing their growth rates, which are sometimes used as proxies for fitness. However, analogous to mathematical models for epidemics, the growth rate can be defined as a function of mechanistic traits: the basic reproduction number (the average number of cells each infected cell infects) and the mean generation time (the average length of a replication cycle). Fitting a model to previously published and newly generated data from experiments in human lung cells, we compared estimates of growth rate, reproduction number and generation time for six influenza A strains. Of four strains in previously published data, A/Canada/RV733/2003 (seasonal H1N1) had the lowest basic reproduction number, followed by A/Mexico/INDRE4487/2009 (pandemic H1N1), then A/Indonesia/05/2005 (spill-over H5N1) and A/Anhui/1/2013 (spill-over H7N9). This ordering of strains was preserved for both generation time and growth rate, suggesting a positive biological correlation between these quantities which have not been previously observed. We further investigated these potential correlations using data from reassortant viruses with different internal proteins (from A/England/195/2009 (pandemic H1N1) and A/Turkey/05/2005 (H5N1)), and the same surface proteins (from A/Puerto Rico/8/34 (lab-adapted H1N1)). Similar correlations between traits were observed for these viruses, confirming our initial findings and suggesting that these patterns were related to the degree of human adaptation of internal genes. Also, the model predicted that strains with a smaller basic reproduction number, shorter generation time and slower growth rate underwent more replication cycles by the time of peak viral load, potentially accumulating mutations more quickly. These results illustrate the utility of mathematical models in inferring traits driving observed differences in in vitro growth of influenza strains.

Modeling Competitive Mixtures With the Lotka-Volterra Framework for More Complex Fitness Assessment Between Strains

With increasing resolution of microbial diversity at the genomic level, experimental and modeling frameworks that translate such diversity into phenotypes are highly needed. This is particularly important when comparing drug-resistant with drug-sensitive pathogen strains, when anticipating epidemiological implications of microbial diversity, and when designing control measures. Classical approaches quantify differences between microbial strains using the exponential growth model, and typically report a selection coefficient for the relative fitness differential between two strains. The apparent simplicity of such approaches comes with the costs of limiting the range of biological scenarios that can be captured, and biases strain fitness estimates to polarized extremes of competitive exclusion. Here, we propose a mathematical and statistical framework based on the Lotka-Volterra model, that can capture frequency-dependent competition between microbial strains within-host and upon transmission. As a proof-of-concept, the model is applied to a previously-published dataset from in-vivo competitive mixture experiments with influenza strains in ferrets (McCaw et al., 2011). We show that for the same data, our model predicts a scenario of coexistence between strains, and supports a higher bottleneck size in the range of 35–145 virions transmitted from donor to recipient host. Thanks to its simplicity and generality, such framework could be applied to other ecological scenarios of microbial competition, enabling a more complex and nuanced view of possible outcomes between two strains, beyond competitive exclusion.

Utilising animal models to evaluate oseltamivir efficacy against influenza A and B viruses with reduced in vitro susceptibility

The neuraminidase (NA) inhibitor (NAI) oseltamivir (OST) is the most widely used influenza antiviral drug. Several NA amino acid substitutions are reported to reduce viral susceptibility to OST in in vitro assays. However, whether there is a correlation between the level of reduction in susceptibility in vitro and the efficacy of OST against these viruses in vivo is not well understood. In this study, a ferret model was utilised to evaluate OST efficacy against circulating influenza A and B viruses with a range of in vitro generated 50% inhibitory concentrations (IC50) values for OST. OST efficacy against an A(H1N1)pdm09 and an A(H1N1)pdm09 virus with the H275Y substitution in neuraminidase was also tested in the macaque model. The results from this study showed that OST had a significant impact on virological parameters compared to placebo treatment of ferrets infected with wild-type influenza A viruses with normal IC50 values (~1 nM). However, this efficacy was lower against wild-type influenza B and other viruses with higher IC50 values. Differing pathogenicity of the viruses made evaluation of clinical parameters difficult, although some effect of OST in reducing clinical signs was observed with influenza A(H1N1) and A(H1N1)pdm09 (H275Y) viruses. Viral titres in macaques were too low to draw conclusive results. Analysis of the ferret data revealed a correlation between IC50 and OST efficacy in reducing viral shedding but highlighted that the current WHO guidelines/criteria for defining normal, reduced or highly reduced inhibition in influenza B viruses based on in vitro data are not well aligned with the low in vivo OST efficacy observed for both wild-type influenza B viruses and those with reduced OST susceptibility.

Triple combination therapy of favipiravir plus two monoclonal antibodies eradicates influenza virus from nude mice

Prolonged treatment of immunocompromised influenza patients with viral neuraminidase (NA) inhibitors is required, because the immune system of such patients fails to eradicate the viruses. Here, we attempted to eradicate influenza virus from the respiratory organs of nude mice, which is a model of immunocompromised hosts, by using combination therapy of the viral polymerase inhibitor favipiravir and monoclonal antibodies (mAbs) against the receptor-binding site (RBS) and stem of viral hemagglutinin (HA). Although monotherapy or combination therapy of two antivirals (two mAbs or favipiravir plus a mAb) suppressed virus replication, they failed to eradicate viruses from nude mice. In contrast, the triple combination therapy of favipiravir plus anti-Stem and anti-RBS mAbs completely stopped virus replication in nude mice, resulting in virus clearance. Triple combination approaches should be considered for the treatment of human immunocompromised patients with severe influenza.

[The model of prediction of the inhibition of neuraminidases of influenza A and B based on a reduced set of energy terms]

The overall model for prediction of IC50 values for inhibitors of neuraminidase influenza virus A and B has been created. It combines data about IC50 values of complexes of 40 variants of neuraminidases of influenza A (7 serotypes) and B and three known inhibitors (oseltamivir, zanamivir, peramivir). The model also uses only data of enthalpy contributions to the potential energy of inhibitor/protein and substrate (MUNANA)/protein complexes. The calculation procedures are ported to use software with support of GPU accelerators, that significant decrease the computation time. The corresponding correlation coefficient (R²) for pIC50 prediction was within 0.45-0.58, the SEM values of around 0.7 (the range of used pIC50 data set is from 4.55 to 10.22).

Influenza B viruses from different genetic backgrounds are variably impaired by neuraminidase inhibitor resistance-associated substitutions

Identifying evolutionary routes to antiviral resistance among influenza viruses informs molecular-based resistance surveillance and clinical decisions. To improve antiviral management and understand whether clinically identified neuraminidase (NA) inhibitor (NAI) resistance-associated markers affect influenza B viruses of the Victoria- or Yamagata-lineages differentially, we generated a panel of NAI-resistant viruses (carrying E105K, G145E, R150K, D197N, I221 L/N/T/V, H273Y, N294S, or G407S substitutions; B numbering) in B/Brisbane/60/2008 (BR/08) and B/Phuket/3073/2013 (PH/13). In both backgrounds, I221 L/N/T/V resulted in reduced or highly reduced inhibition (HRI) by one to three currently available NAIs. D197N reduced inhibition by all NAIs in BR/08 but only by oseltamivir and peramivir in PH/13; R150K caused HRI by all NAIs in PH/13. Although PH/13 generally retained or enhanced NA activity in the presence of the substitutions, enzymatic activity in BR/08 was detrimentally affected. Similarly, substrate affinity and catalysis were relatively stable in PH/13, but not in the BR/08 variants. E105K, R150K, and D197N attenuated replication efficiency of BR/08 in vitro and in mice; only E105K had this effect in PH/13. Notably, the I221 L/N/T/V substitutions did not severely impair replication, particularly in PH/13. Overall, our data show differential effects of NA substitutions in representative Victoria- and Yamagata-lineage viruses, suggesting distinct evolution of these viruses caused variable fitness and NAI susceptibility profiles when similar key NA substitutions arise. Because the viruses harboring the I221 NA substitutions displayed undiminished fitness and are commonly reported, this position is likely to be the most clinically relevant marker for NAI resistance among contemporary influenza B viruses.

Characterization of substitutions in the neuraminidase of A(H7N9) influenza viruses selected following serial passage in the presence of different neuraminidase inhibitors

Avian A(H7N9) infections in humans have been reported in China since 2013 and are of public health concern due to their severity and pandemic potential. Oseltamivir and peramivir are neuraminidase inhibitors (NAIs) routinely used for the treatment of A(H7N9) infections, but variants with reduced sensitivity to these drugs can emerge in patients during treatment. Zanamivir and laninamivir are NAIs that are used less frequently. Herein, we performed in vitro serial passaging experiments with recombinant viruses, containing the neuraminidase (NA) from influenza A/Anhui/1/13 (H7N9) virus, in the presence of each NAI, to determine whether variants with reduced sensitivity would emerge. NA substitutions were characterized for their effect on the NA enzymatic activity and surface expression of the A/Anhui/1/13 (Anhui/1) NA, as well as NAs originating from contemporary A(H7N9) viruses of the Yangtze River Delta and Pearl River Delta lineages. In vitro passage in the presence of oseltamivir, peramivir and laninamivir selected for substitutions associated with reduced sensitivity (E119D, R292K and R152K), whereas passage in the presence of zanamivir did not select for any viruses with reduced sensitivity. All the NA substitutions significantly reduced activity, but not the expression of the Anhui/1 NA. In contemporary N9 NAs, all substitutions tested significantly reduced NA enzyme function in the Yangtze River lineage background, but not in the Pearl River Delta lineage background. Overall, these findings suggest that zanamivir may be less likely than the other NAIs to select for resistance in A(H7N9) viruses and that the impact of substitutions that reduce NAI susceptibility or enzyme function may be less in A(H7N9) viruses from the Pearl River lineage.

In Vitro Properties and Virulence of Contemporary Recombinant Influenza B Viruses Harboring Mutations of Cross-Resistance to Neuraminidase Inhibitors

Three neuraminidase inhibitors (NAIs: Oseltamivir, zanamivir and peramivir) are currently approved in many countries for the treatment of influenza A and B infections. The emergence of influenza B viruses (IBVs) containing mutations of cross-resistance to these NAIs constitutes a serious clinical threat. Herein, we used a reverse genetics system for the current B/Phuket/3073/2013 vaccine strain to investigate the impact on in vitro properties and virulence of H136N, R152K, D198E/N, I222T and N294S NA substitutions (N2 numbering), reported by the World Health Organization (WHO) as clinical markers of reduced or highly-reduced inhibition (RI/HRI) to multiple NAIs. Recombinant viruses were tested by NA inhibition assays. Their replicative capacity and virulence were evaluated in ST6GalI-MDCK cells and BALB/c mice, respectively. All NA mutants (excepted D198E/N) showed RI/HRI phenotypes against ≥ 2 NAIs. These mutants grew to comparable titers of the recombinant wild-type (WT) IBV in vitro, and some of them (H136N, I222T and N294S mutants) induced more weight loss and mortality in BALB/c mice in comparison to the recombinant WT IBV. These results demonstrate that, in contemporary IBVs, some NA mutations may confer RI/HRI phenotypes to existing NAIs without altering the viral fitness. This reinforces the need for development of novel antiviral strategies with different mechanisms of action.