Characterization of Influenza A and B Viruses Circulating in Southern China During the 2017–2018 Season

The trivalent seasonal influenza vaccine was the only approved and available vaccine during the 2016–2018 influenza seasons. It did not include the B/Yamagata strain. In this study, we report an acute respiratory disease outbreak associated with influenza B/Yamagata infections in Guangzhou, Southern China (January through March, 2018). Among the 9914 patients, 2241 (22.6%) were positive for the influenza B virus, with only 312 (3.1%) positive for the influenza A virus. The influenza B/Yamagata lineage dominated during this period in Southern China. The highest incidence of influenza A virus infection occurred in the children aged 5–14 years. In contrast, populations across all age groups were susceptible to the influenza B virus. Phylogenetic, mutations, and 3D structure analyses of hemagglutinin (HA) genes were performed to assess the vaccine-virus relatedness. The recommended A/H1N1 vaccine strain (A/Michigan/45/2015) during both 2017–2018 and 2018–2019 was antigen-specific for these circulating isolates (clade 6B.1) in Spring 2018. An outbreak of influenza B/Yamagata (clade 3) infections in 2018 occurred during the absence of the corresponding vaccine during 2016–2018. The recommended influenza B/Yamagata vaccine strain (B/Phuket/3073/2013) for the following season (2018–2019) was antigen-specific. Although there were only a few influenza B/Victoria infections in Spring 2018, five amino acid mutations were identified in the HA antigenic sites of the 19 B/Victoria isolates (clade 1A), when compared with the 2016–2018 B/Victoria vaccine strain. The number was larger than expected and suggested that the influenza B HA gene may be more variable than previously thought. One of the mutations (K180N) was noted to likely alter the epitope and to potentially affect the viral antigenicity. Seven mutations were also identified in the HA antigenic sites of 2018–2020 B/Victoria vaccine strain, of which some or all may reduce immunogenicity and the protective efficacy of the vaccine, perhaps leading to more outbreaks in subsequent seasons. The combined epidemiological, phylogenetic, mutations, and 3D structural analyses of the HA genes of influenza strains reported here contribute to the understanding and evaluation of how HA mutations affect vaccine efficacy, as well as to providing important data for screening and selecting more specific, appropriate, and effective influenza vaccine candidate strains.

Introduction

Influenza viruses are the major cause of acute respiratory diseases in humans, causing several serious global pandemics because of their transmission dynamics and great antigenic variability. Compounding this is the potential of emergent strains originating in swine and avian hosts (Mostafa et al., 2018), as well as from reverse zoonosis from human to animals (Morens et al., 2013; Nelson and Vincent, 2015). For both directions of transmission, knowing the distribution of strains in human populations is important for perhaps predicting strains for the coming season.

Influenza A virus presents as a permutation of 18 hemagglutinin (HA) and 11 neuraminidase (NA) possible subtypes, with a wide range of hosts (Cauldwell et al., 2014). The latest two subtypes H17N10 and H18N11 are of bat origins (Tong et al., 2012, 2013). Influenza B virus, including the Victoria and Yamagata lineages, is hosted by humans and seals (Osterhaus et al., 2000). Although there is no direct evidence that influenza B is transmitted between humans and animals, it has been observed that influenza B can infect seal, swine, ferret, guinea pig, pheasant, dog and horse (Osterhaus et al., 2000; Pica et al., 2012; Huang et al., 2014; Ran et al., 2015; Pascua et al., 2016). Pandemic influenza is usually caused by influenza A virus, due to its rapid antigenic variation, strong replication capacity, and transmission ability associated with genetic reassortment (Chan et al., 2010; White and Lowen, 2018). Influenza B virus, often neglected, has been circulating and, in some seasons, has predominated over influenza A, particularly among children (Paul Glezen et al., 2013; Tewawong et al., 2015; Furuse et al., 2016). Vaccines against influenza A/H1N1 and A/H3N2 have been used effectively for more than 20 years for controlling and preventing potential pandemics. However, influenza vaccines against the B/Victoria and B/Yamagata lineages were used interchangeably. In the 2016–2018 seasons, the trivalent influenza vaccine did not include the B/Yamagata lineage.

A high incidence of influenza-like illnesses was reported in January 2018 by the China CDC¹. To investigate whether this was due to the unavailability of B/Yamagata vaccine, we collected 9914 nasal swabs from patients presenting with influenza-like symptoms in Guangzhou, Southern China from January through March in 2018. Most of the patients were outpatients (88.7%). Influenza-positive samples were further subtyped by sequencing and phylogenetic analysis of HA and NA genes. The data provide a view into the molecular epidemiology of the influenza viruses circulating in the current population, which may, in turn, confirm the efficacy of the influenza vaccines administered as well as provide insights for the future development and deployment of effective subtype or lineage-specific influenza vaccines.

Materials and Methods

Influenza Diagnostic Tests

From January through March 2018, nasal swab specimens from 9914 patients with influenza-like symptoms were collected and tested for influenza viral antigens (Li et al., 2019) using the Clearview^® Exact Influenza A&B kit with the colloidal gold detection method (Jennings et al., 2009)². Briefly, each nasal swab was eluted with 192 μL lysis buffer and the immunochromatographic strip in the kit was inserted into the eluent and incubated for 15 min at room temperature. Then the results were determined according to the manufacturer’s instructions. The interaction between the colloidal gold conjugated antibodies and the influenza antigens will indicate red in the test line. This study was approved by the Institutional Review Board of the Affiliated Cancer Hospital of Guangzhou Medical University in accordance with the Declaration of Helsinki, with the patient consent for using left-over specimens waived.

RNA Extraction, RT-PCR and Sequencing

Each nasal swab was eluted with 1 mL DMEM. The eluent was immediately stored at −80°C. Subsequently, 1 mL Trizol was used to extract viral RNA from the 200 μL of nasal swab eluent, and 200 μL chloroform was added immediately and centrifuged at 12,000 × g. The upper aqueous layer was transferred to a tube, with 1 mL 100% isopropanol for RNA precipitation. Following centrifuge at 12,000 × g, 1 mL of 75% cold ethanol to the pellet as a wash. The RNA was eluted into 20 μL RNase-free water. This RNA was reverse transcribed into cDNA using the PrimerScript^TM RT Reagent Kit with gDNA Eraser (TAKARA), according to the manufacturer’s protocol. All primers used in this experiment are based on the “WHO protocols for molecular diagnosis of influenza virus”³ (Table 1). The cDNA was amplified using the Premix Taq^TM (TaKaRa^TM Taq Version 2.0), under the following conditions: Initial denaturation at 94°C for 2 min; amplification in 40 cycles of 94°C for 30 s/50°C for 30 s/72°C for 2 min; and final extension for 10 min at 72°C. These amplified sequences were confirmed by Sanger DNA sequencing using the primers listed in Table 1 as sequencing primers.

TABLE 1

Table 1. Primers used for RT and PCR amplification of influenza A and B viruses.

Phylogenetic Sequence Analyses

Phylogenetic analysis was performed by Molecular Evolutionary Genetics Analysis (MEGA X)⁴ (Kumar et al., 2018). The phylogenetic tree was generated by the maximum-likelihood method with 1,000 bootstrap replicates using the MEGA X by applying default parameters (Zhang et al., 2017, 2019). Sequences of HA and NA genes were retrieved from Influenza Research Database⁵ are used for comparison. The gene reference sequences of H1N1, H3N2, and Victoria and Yamagata lineages were extracted from published sequences, using BLAST, obtained from strains that circulated worldwide in recent years, as well as China strains and northern hemisphere vaccine strains recommended by WHO and the reference strains of known clades as reported by WHO.

Amino Acid Mutation Analysis

Amino acid mutations in the HA genes of 23 influenza A/H1N1 isolates from residues 7 to 566 were compared to homologous sequences in the 2017–2019 vaccine strain (A/Michigan/45/2015). In addition, amino acid mutations of the HA genes of 80 influenza B/Yamagata isolates from residues 96–215 were compared to the equivalent sequences from the 2018–2019 vaccine strain (B/Phuket/3073/2013). Amino acid mutations in the HA genes from the 19 influenza B/Victoria isolates and the 2018–2020 B/Victoria vaccine strain from residues 93 to 180 were also compared to the counterparts in the 2016–2018 northern hemisphere vaccine strain (B/Brisbane/60/2008). These analyses were compared using the software CLUSTAL.

3D Protein Structure Analysis

The SWISS-MODEL server and web-tool⁶ (Waterhouse et al., 2018) was used to predict the 3D structure of the HA protein of the two northern hemisphere vaccine strains (A/Michigan/45/2015 and B/Phuket/3073/2013). SWISS-MODEL searches for related protein structure templates from the SMTL database by BLAST and HHblits were conducted to ensure accuracy and sensitivity as much as possible. The crystal structure of the B/Brisbane/60/2008 vaccine strain of the HA protein was retrieved from NCBI. The PDB ID: 4FQM. PyMOL software⁷ as used to overlay the positions of the antigenic sites and mutation sites on the HA crystal structure of influenza A and B isolates.

HA and NA Gene Sequences Annotation and GenBank Accession Numbers

Sequence data were assembled with the SeqMan Pro software 7.0.1 (DNASTAR, Inc., Madison, WI, United States). Nucleotide and amino acid sequences of the HA and NA genes were aligned using the CLUSTAL and BLAST software. These HA and NA gene sequences from influenza A and B were archived in GenBank with the following accession numbers: (1) MN653601-MN653624 (influenza A, HA gene); (2) MN653549-MN653572 (influenza A, NA gene); and (3) MN653573-MN653600 and MT123908-MT123978 (influenza B, HA gene).

Results and Discussion

The Epidemiology of Influenza A and B Infections

Type-specific colloidal gold test revealed that 2547 of 9914 patients (25.7%) during this outbreak were positive for influenza viruses, of which 312 (3.1%) were identified as influenza A and 2241 (22.6%) as influenza B. The dominant type circulating was influenza B (87.8%). Outpatients of influenza A and B positive accounted for 297 (95.2%) and 2141 (95.5%), respectively. Positive specimens collected during the outbreak were selected randomly (100 specimens for each month) for further subtyping by Sanger sequencing DNA sequencing of the HA and NA genes. The sequence data resulted in the typing of the influenza A strains as 23 H1N1 with one H3N2, and the typing of influenza B as identified 80 Yamagata and 19 Victoria strains.

Influenza cases peaked initially from January 8 to 21, followed by another peak from February 5 to March 4 (Figure 1). This was different from the Northern China epidemic profile (one peak, January 1–21). During the two peaks, there was a total of 1776 influenza positive cases. Among these, influenza B accounted for 1579 cases, with only 203 as influenza A positive. These include six cases testing as both influenza A and B positive.

FIGURE 1

Figure 1. Timeline and type distribution of influenza A and B cases of onset in Guangzhou, Southern China, January–March 2018. The number of cases of each influenza A and B are shown encompassing the two peaks of the epidemic.

During the outbreak, influenza A and B viruses circulated in all age groups. Median ages of the patients were 12 and 17, respectively, with the male to female ratio of 1.1:1. The highest positive rate of influenza A and B infections appeared in children between the ages of 5–14 years, with incidences of 4.2 and 31.7%, respectively (Table 2). All age groups were susceptible to influenza B virus, accounting for 17.9–31.7% of the flu-like cases.

TABLE 2

Table 2. The age distribution of influenza A and B patients in Guangzhou, Southern China, January–March 2018.

The Subtype Distribution and Phylogenetic Analysis of Influenza A and B Viruses

To reveal the subtypes circulating during this outbreak, the HA and NA genes were further analyzed using their sequences. Phylogenetic analysis showed that the HA and NA genes of the 23 H1N1 isolates had close sequence relationships to the strains circulating in Northern China, Japan, Italy, England, and the USA during 2017–2018, as well as to the 2016–2017 vaccine strain (A/California/7/2009/vaccine) and 2017–2019 vaccine strain (A/Michigan/45/2015/vaccine). The phylogenetic analysis indicated that there were few HA and NA genetic variations between the A/H1N1 isolates circulating in 2018 and the 2016–2019 vaccine strains. All of the A/H1N1 strains circulating during January to March, as well as the 2017–2019 vaccine strain (A/Michigan/45/2015), fell into clade 6B.1, which is defined by the HA1 amino acid substitutions S84N, S162N and I216T. The 2016–2017 vaccine strain (A/California/7/2009) and the other four China strains circulating in 2016 formed a clade noted as 6B.2 (Figure 2).

FIGURE 2

Figure 2. Phylogenetic analysis of the HA (A) and NA (B) genes from influenza A/H1N1 viruses. The phylogenetic tree was generated by the maximum-likelihood method with 1,000 bootstrap replicates using the MEGA X by applying default parameters. Bootstrap numbers shown at the nodes mean the percentages of 1000 replications producing the clade. The scale bar indicates the units of nucleotide substitutions per site. Nucleotide sequences of HA and NA genes retrieved from Influenza Research Database (https://www.fludb.org/brc/influenza_sequence_search_segment_display.spg?method=ShowCleanSearch&decorator=influenza) are used as references. Sequences of the isolates obtained in this study are noted (▲), as well as sequences of the northern hemisphere vaccine strains of known clades as reported by WHO as reference (•).

The majority (80/99) of the influenza B isolates belonged to the Yamagata lineage, which was further divided into two sub-branches: One was phylogenetically close to strains circulating in China, Japan, Italy, South Korea, German, Canada, and the USA during 2016–2018 and the other was close to strains B/Georgia/7276/2018, B/Florida/31/2017, B/Idaho/03/2017, B/New Mexico/6335/2017, B/Turkey/5879/2017, and B/West Virginia/13/2017 (Figure 3). Influenza B Yamagata-lineage was separated previously into two major antigenically distinct clades (clades 2 and 3), based on phylogenetic analysis of its HA and NA genes (Barr et al., 2010, 2014). All of the B/Yamagata strains circulating during January to March and the 2018–2019 vaccine strain (B/Phuket/3073/2013) parsed into clade 3.

FIGURE 3

Figure 3. Phylogenetic analysis of the HA genes of influenza B/Yamagata isolates (A) and B/Victoria isolates (B). The phylogenetic trees were generated by the maximum-likelihood method with 1,000 bootstrap replicates using the MEGA X by applying default parameters. Bootstrap numbers shown at the nodes mean the percentages of 1000 replications producing the clade. The scale bar indicates the units of nucleotide substitutions per site Nucleotide sequences of the HA genes retrieved from the Influenza Research Database (https://www.fludb.org/brc/influenza_sequence_search_segment_display.spg?method = ShowCleanSearch& decorator = influenza) are used as references. Sequences of the isolates obtained in this study are noted (▲), as well as sequences of the northern hemisphere vaccine strains of known clades as reported by WHO as reference (•). Δ2 means two amino acid deletions (177K and 178N).

During the 2016–2018 seasons, the trivalent vaccine didn’t include B/Yamagata lineage, so its effectiveness could not be determined. However, it may be speculated that the absence of B/Yamagata vaccine might have contributed to an influenza B outbreak due to the low herd immunity.

The 19 B/Victoria strains clustered within two close sub-branches, 16 of which were close phylogenetically to the 2018–2020 vaccine strain (B/Colorado/06/2017) and three of which were close to the 2016–2018 vaccine strain (B/Brisbane/60/2008) (Figure 3). All of the 19 B/Victoria strains were members of clade 1A.

Amino Acid Mutations and 3D Structure Analysis of the HA Antigenic Sites Indicate the Potential Effect on Vaccines

Current influenza vaccines provide important protection in humans by inducing strain-specific neutralizing antibodies targeting highly variable antigenic epitopes in the globular domain of the HA protein that is in the virus envelope. They play a critical role in host cell recognition, viral binding, as well as the subsequent fusion and entry processes (Shen et al., 2013; Qi et al., 2018). Amino acid mutations of the HA genes of circulating influenza isolates and the predicted 3D structure were further analyzed, and compared to the northern hemisphere vaccine strains.

When compared with the A/H1N1 2017–2019 vaccine strain (A/Michigan/45/2015), there were three major amino acid mutations within the HA gene of all 23 isolates (S91R, S181T, and I312V), of which one mutation, S91R, located in the Cb antigenic site was identified in all 23 isolates (Figure 4A). Another mutation P154S in the Ca2 antigenic site was also found in two isolates (SMU214 and SMU941). Except for these mutations, no additional mutations were located within these antigenic sites of HA gene: Sa, Sb, Ca1, Ca2, Cb (Huang et al., 2013; Matsuzaki et al., 2014). S181T was identified in all the 23 isolates, although the meaning of the mutation is unknown (Table 3). To be clear, it is not known if these patients were vaccinated for influenza. Whether these mutations actually changed the immunogenicity of these isolates needs to be investigated further. However, the incidence of influenza A was significantly lower than influenza B in Jan-Mar 2018. The use of the influenza A vaccine appears to have been effective in controlling the potential influenza A outbreak.

FIGURE 4

Figure 4. Structural modeling of the HA mutations. The hemagglutinin structures of influenza A/H1N1 (A), B/Yamagata (B), and B/Victoria (C) vaccine strains are presented, in which the antigenic epitopes of the protein are mapped and marked in yellow. Amino acid mutations of the isolates are either marked in pink, if not located in antigenic sites, or in red, if located in antigenic sites. The frequency was also noted.

TABLE 3

Table 3. Amino acid mutations in the HA genes of A/H1N1 isolates.

Because there was no influenza B/Yamagata vaccine used during the 2017–2018 season, the 2018–2019 influenza B/Yamagata vaccine strain (B/Phuket/3073/2013) was used as reference. The HA gene of most of the isolates was conserved. Only five out of 80 isolates contained one non-synonymous substitution, located in the 120-loop, 150-loop, 160-loop, and 190-helix sites (Wang et al., 2008), respectively (marked in red in Table 4). One isolate (SMU926) harbored four mutations in the 160-loop site. Except for these mutations, no other changes were found in the HA antigenic sites (Table 4). The mutation L187Q was noted in all 80 isolates, but it was not found in any antigenic sites and therefore may not affect the antigenicity. At amino acid residue 179, three of eighty isolates had amino acid mutations in the antigenic epitopes (Figure 4B). For the other six mutations, they were only found in one isolate for each mutation. These indicated that the 2018–2019 B/Yamagata vaccine strain was antigen-specific and appropriate for the 2018 circulating B/Yamagata strains.

TABLE 4

Table 4. Amino acid mutations in the HA genes of B/Yamagata isolates.

When compared with the B/Victoria 2016–2018 vaccine strain (B/Brisbane/60/2008), all 19 isolates had one common amino acid mutation (K180N), located in the HA 160-loop (Table 5). Five mutations located in the HA antigenic epitopes were identified in the 19 B/Victoria isolates: K180N (19/19), I132V (16/19), N144D (15/19), I161V (3/19), H137Q or H137N (1/19), respectively (Figure 4C). Note that the mutations K180N, I132V, and N144D occurred in most of the isolates. In the 150-loop, three isolates contained a mutation I161V. Additionally, 17 of 19 isolates had three mutations located in the three antigenic sites (120-loop, 150-loop, and 160-loop). When compared with the B/Victoria 2018–2020 vaccine strain (B/Colorado/06/2017), surprisingly, all 19 isolates had K180N mutation in the 160-loop and 15 of 19 isolates had a G144D mutation in the 120-loop. A previous study found that the amino acid substitution N180K altered the epitope and affected the viral antigenicity. Human antibodies did not substantially inhibit the hemagglutination reaction in the hemagglutination inhibition tests (Nakagawa et al., 2005). Such variants may be important in future epidemics. Except for the B/Colorado/06/2017/vaccine strain, which had two deletions (Δ2) in the HA antigenic sites (177K and 178N), the other isolates had no amino acid deletions at any of the antigenic sites (Table 5). A larger than expected number of mutations found within HA multiple antigenic epitopes of B/Victoria isolates suggest that the HA gene may be more variable than previously thought (Lindstrom et al., 1999). The isolates accumulated mutations at several antigenic sites, which may reduce the protective efficacy of the vaccine. Therefore, the B/Vitoria vaccine strain used in 2018–2020 might not be appropriately selected for Southern China should this be the case. According to the report of the National Health Commission of China, the B/Victoria strain re-dominated in China in 2019⁸. In the context of our study, this indicates that the protective efficacy of 2018–2019 B/Victoria vaccine is likely reduced due to the mutations in the antigenic sites of the HA protein, particularly if the vaccination rate in China remained at the same level in 2019 as in prior years.

TABLE 5

Table 5. Amino acid mutations in the HA genes of B/Victoria isolates.

Conclusion

According to the report of the National Health Commission of China, from January through March of 2018, 487,773 people were infected with influenza, of which 119 died. During the outbreak, both influenza A and B dominated in Northern China. In contrast, our study found that influenza B, rather than influenza A, dominated in Southern China. Throughout the years, influenza A viruses have attracted a great deal of attention and have caused several global pandemics due to their strong transmission capacity and great variability. In contrast, there is less literature focused on the epidemiology and societal burdens of influenza B, particularly outside the United States and Europe (Caini et al., 2015). Nevertheless, influenza B remains an critical respiratory pathogen that imparts an important and large public health impact (Paul Glezen et al., 2013), and studies such as this provide important context for both influenza A and B, as well as insights into the dynamics of their epidemiology. It should be noted that a domination of influenza B/Yamagata in both Northern and Southern China has occurred twice, in the past 2007–2008 and in 2014–2015 (Yu et al., 2013; Yang et al., 2018). In this study, we found that, 3 years later, influenza B/Yamagata dominated once again in Southern China.

The highest incidence of influenza infection occurred among the 5–14-year-old children, which indicates that children and adolescents were more susceptible to influenza viruses. Outpatients accounted for more than 95% of influenza cases, which suggested that the illness was not particularly severe during this outbreak. H1N1 was the major subtype of influenza A (95.8%) and the Yamagata lineage was the major lineage of influenza B (80.8%). The outbreak of influenza viruses in Southern China in 2018 was mainly caused by Yamagata lineage, which was phylogenetically close to other influenza B strains circulating worldwide in the same timeframe.

With respect to the phylogenetic, mutation, and 3D structure analysis of HA genes, the A/H1N1 vaccine strain used in 2017–2018 was antigen-specific and appeared to be effective in this season. However, the low influenza vaccination rate in China might have contributed to the emergence of H1N1 cases⁹. Due to the lack of the influenza B/Yamagata vaccine during 2016–2018, the outbreak of the influenza B/Yamagata virus in 2018 might have been expected. According to our analysis, the influenza B/Yamagata vaccine strain used during the subsequent season (2018–2019) was antigen-specific for the circulating strains. An outbreak of influenza B/Yamagata infections did not occur in 2019 perhaps due to the use of the inactivated B/Yamagata vaccine. However, the HA gene of the 2018 B/Victoria isolates had accumulated mutations in several antigenic sites, which may have potentially changed the immunogenicity. Furthermore, when compared with the B/Vitoria 2018–2019 vaccine strain, more HA mutations were found in the circulating isolates, which may have reduced the protective efficacy of B/Vitoria 2018–2019 vaccine. This is an indication that the vaccine strain might not have been the best candidate vaccine for Southern China in future influenza seasons. According to the reports of the United States CDC, although the influenza vaccination rate for adults was as high as 45.3% in 2018–2019, an increase of 8.2% from the 2017–2018 influenza season, a severe influenza infection still broke out in 2019–2020, causing more than 10,000 deaths, primarily caused by influenza B/Victoria lineage (Owusu et al., 2020)¹⁰. HA gene variations may and should be considered for evaluating potential vaccine efficacy.

Previous studies demonstrated that vaccination with the live and attenuated influenza vaccine elicited lung CD4⁺ and virus-specific CD8⁺ T cell responses, similar in phenotype to those generated by influenza virus infection, and ultimately established lung Tissue-resident memory T cells (TRM) capable of providing long-term, hetero-subtypic protection to multiple, non-vaccine influenza strains. In contrast, vaccination with inactivated influenza vaccine generated durable, strain-specific humoral immunity but failed to elicit T cell responses (Zens et al., 2016). In China, all the influenza vaccines approved and used are inactivated vaccines. Therefore, the effect of mutations on HA T cell epitopes was not analyzed.

The epidemiological, phylogenetic, mutation, and 3D structural analyses presented in this study, contribute to a better understanding of circulating influenza strains by revealing critical mutations, allowing for an evaluation of vaccine efficacy, and providing a basis for the improved selection of more specific and effective influenza vaccine candidate strains. A continuing surveillance of the sequence variations of the HA genes is important for managing, controlling, and limiting future influenza outbreaks and pandemics. These surveys and characterizations of circulating influenza strains are especially important should there be emergent or re-emergent influenza viral pathogens crossing the human-animal interface.

Data Availability Statement

All datasets generated for this study are included in the article.

Author Contributions

QZ, YY, and ZY conceived and designed the experiments. ZY collected and provided clinical samples. YY performed the experiments. YY, JO, SZ, KM, WL, WG, XW, JZ, and QZ analyzed the data. BZ, WZ, CW, WS, JW, DS, ZY, and QZ contributed to the preparation of the manuscript.

Funding

This work was supported by grants from the National Key Research and Development Program of China (2018YFE0204503), National Natural Science Foundation of China (31570155, 81701995 and 31370199), Natural Science Foundation of Guangdong Province (2018B030312010) as well as from the Guangzhou Healthcare Collaborative Innovation Major Project (201803040004 and 201803040007).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank Dr. Huachen Zhu at the University of Hong Kong, Dr. Yuhai Bi at the Center for Influenza Research and Early-Warning, Chinese Academy of Sciences (CASCIRE), Dr. Sooksan Wong at Guangzhou Medial University, and Lirong Zou at Guangdong CDC for their constructive suggestions on the data analysis.

Footnotes

1. ^ http://www.chinacdc.cn/jkzt/crb/bl/lxxgm/zstd/201801/t20180108_158017.html
2. ^ https://www.alere.com/zh/zh/product-details/clearview-exact-influenza-a-and-b.html
3. ^ https://www.who.int/influenza/gisrs_laboratory/molecular_diagnosis/en
4. ^ http://www.megasoftware.net/index.html
5. ^ https://www.fludb.org/brc/influenza_sequence_search_segment_display.spg?method=ShowCleanSearch&decorator=influenza
6. ^ https://swissmodel.expasy.org/
7. ^ https://pymol.org/2/
8. ^ http://www.chinaivdc.cn/cnic/zyzx/lgzb/
9. ^ http://www.chinacdc.cn/mtbd_8067/201810/t20181030_196735.html
10. ^ https://www.cdc.gov/flu/fluvaxview/coverage-1819estimates.htm

References

Barr, I. G., McCauley, J., Cox, N., Daniels, R., Engelhardt, O. G., Fukuda, K., et al. (2010). Epidemiological, antigenic and genetic characteristics of seasonal influenza A(H1N1), A(H3N2) and B influenza viruses: basis for the WHO recommendation on the composition of influenza vaccines for use in the 2009-2010 northern hemisphere season. Vaccine 28, 1156–1167. doi: 10.1016/j.vaccine.2009.11.043

PubMed Abstract | CrossRef Full Text | Google Scholar

Barr, I. G., Russell, C., Besselaar, T. G., Cox, N. J., Daniels, R. S., Donis, R., et al. (2014). WHO recommendations for the viruses used in the 2013-2014 Northern Hemisphere influenza vaccine: epidemiology, antigenic and genetic characteristics of influenza A(H1N1)pdm09, A(H3N2) and B influenza viruses collected from October 2012 to January 2013. Vaccine 32, 4713–4725. doi: 10.1016/j.vaccine.2014.02.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Caini, S., Huang, Q. S., Ciblak, M. A., Kusznierz, G., Owen, R., Wangchuk, S., et al. (2015). Epidemiological and virological characteristics of influenza B: results of the global influenza B study. Influenza Other Respir. Viruses 9(Suppl. 1), 3–12. doi: 10.1111/irv.12319

PubMed Abstract | CrossRef Full Text | Google Scholar

Cauldwell, A. V., Long, J. S., Moncorge, O., and Barclay, W. S. (2014). Viral determinants of influenza A virus host range. J. Gen. Virol. 95(Pt 6), 1193–1210. doi: 10.1099/vir.0.062836-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Chan, K.-H., Zhang, A. J. X., To, K. K. W., Chan, C. C. S., Poon, V. K. M., Guo, K., et al. (2010). Wild type and mutant 2009 pandemic influenza A (H1N1) viruses cause more severe disease and higher mortality in pregnant BALB/c Mice. PLoS One 5:e13757. doi: 10.1371/journal.pone.0013757

PubMed Abstract | CrossRef Full Text | Google Scholar

Furuse, Y., Odagiri, T., Tamaki, R., Kamigaki, T., Otomaru, H., Opinion, J., et al. (2016). Local persistence and global dissemination play a significant role in the circulation of influenza B viruses in Leyte Island, Philippines. Virology 492, 21–24. doi: 10.1016/j.virol.2016.02.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, P., Yu, S. Y., Wu, C. Y., and Liang, L. J. (2013). Highly conserved antigenic epitope regions of hemagglutinin and neuraminidase genes between 2009 H1N1 and seasonal H1N1 influenza: vaccine considerations. J. Transl. Med. 11:47. doi: 10.1186/1479-5876-11-47

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, S. S. H., Banner, D., Paquette, S. G., Leon, A. J., Kelvin, A. A., and Kelvin, D. J. (2014). Pathogenic influenza B virus in the ferret model establishes lower respiratory tract infection. J. Gen. Virol. 95(Pt 10), 2127–2139. doi: 10.1099/vir.0.064352-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Jennings, L. C., Skopnik, H., Burckhardt, I., Hribar, I., Del Piero, L., and Deichmann, K. A. (2009). Effect of rapid influenza testing on the clinical management of paediatric influenza. Influenza Other Respir. Viruses 3, 91–98. doi: 10.1111/j.1750-2659.2009.00079.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, S., Stecher, G., Li, M., Knyaz, C., and Tamura, K. (2018). MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 35, 1547–1549. doi: 10.1093/molbev/msy096

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, W., Liu, L., Chen, L., and Shang, S. (2019). Evaluation of a commercial colloidal gold assay for detection of influenza A and B virus in children’s respiratory specimens. Fetal Pediatr. Pathol. 39, 93–98. doi: 10.1080/15513815.2019.1639088

PubMed Abstract | CrossRef Full Text | Google Scholar

Lindstrom, S. E., Hiromoto, Y., Nishimura, H., Saito, T., Nerome, R., and Nerome, K. (1999). Comparative analysis of evolutionary mechanisms of the hemagglutinin and three internal protein genes of influenza B virus: multiple cocirculating lineages and frequent reassortment of the NP, M, and NS genes. J. Virol. 73, 4413–4426.

PubMed Abstract | Google Scholar

Matsuzaki, Y., Sugawara, K., Nakauchi, M., Takahashi, Y., Onodera, T., Tsunetsugu-Yokota, Y., et al. (2014). Epitope mapping of the hemagglutinin molecule of A/(H1N1)pdm09 influenza virus by using monoclonal antibody escape mutants. J. Virol. 88, 12364–12373. doi: 10.1128/jvi.01381-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Morens, D. M., Taubenberger, J. K., and Fauci, A. S. (2013). H7N9 avian influenza A virus and the perpetual challenge of potential human pandemicity. mBio 4:e00445-13. doi: 10.1128/mBio.00445-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Mostafa, A., Abdelwhab, E. M., Mettenleiter, T. C., and Pleschka, S. (2018). Zoonotic potential of influenza A viruses: a comprehensive overview. Viruses 10:497. doi: 10.3390/v10090497

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakagawa, N., Kubota, R., and Okuno, Y. (2005). Variation of the conserved neutralizing epitope in influenza B virus victoria group isolates in Japan. J. Clin. Microbiol. 43, 4212–4214. doi: 10.1128/JCM.43.8.4212-4214.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Nelson, M. I., and Vincent, A. L. (2015). Reverse zoonosis of influenza to swine: new perspectives on the human-animal interface. Trends Microbiol. 23, 142–153. doi: 10.1016/j.tim.2014.12.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Osterhaus, A. D., Rimmelzwaan, G. F., Martina, B. E., Bestebroer, T. M., and Fouchier, R. A. (2000). Influenza B virus in seals. Science 288, 1051–1053.

PubMed Abstract | Google Scholar

Owusu, D., Hand, J., Tenforde, M. W., Feldstein, L. R., DaSilva, J., Barnes, J., et al. (2020). Early season pediatric influenza B/victoria virus infections associated with a recently emerged virus subclade – Louisiana, 2019. MMWR Morb. Mortal. Wkly. Rep. 69, 40–43. doi: 10.15585/mmwr.mm6902e1

PubMed Abstract | CrossRef Full Text | Google Scholar

Pascua, P. N., Marathe, B. M., Burnham, A. J., Vogel, P., Webby, R. J., Webster, R. G., et al. (2016). Competitive fitness of influenza B viruses possessing E119A and H274Y neuraminidase inhibitor resistance-associated substitutions in Ferrets. PLoS One 11:e0159847. doi: 10.1371/journal.pone.0159847

PubMed Abstract | CrossRef Full Text | Google Scholar

Paul Glezen, W., Schmier, J. K., Kuehn, C. M., Ryan, K. J., and Oxford, J. (2013). The burden of influenza B: a structured literature review. Am. J. Public Health 103, e43–e51. doi: 10.2105/ajph.2012.301137

PubMed Abstract | CrossRef Full Text | Google Scholar

Pica, N., Chou, Y. Y., Bouvier, N. M., and Palese, P. (2012). Transmission of influenza B viruses in the guinea pig. J. Virol. 86, 4279–4287. doi: 10.1128/jvi.06645-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Qi, Z. L., Hu, H. Y., Wang, Z. H., Wang, G. H., Li, Y., Zhao, X. R., et al. (2018). Antibodies against H1N1 influenza virus cross-react with -cells of pancreatic islets. J. Diabet. Invest. 9, 265–269. doi: 10.1111/jdi.12690

PubMed Abstract | CrossRef Full Text | Google Scholar

Ran, Z., Shen, H., Lang, Y., Kolb, E. A., Turan, N., Zhu, L., et al. (2015). Domestic pigs are susceptible to infection with influenza B viruses. J. Virol. 89, 4818–4826. doi: 10.1128/jvi.00059-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, X. T., Zhang, X. X., and Liu, S. W. (2013). Novel hemagglutinin-based influenza virus inhibitors. J. Thor. Dis. 5, S149–S159. doi: 10.3978/j.issn.2072-1439.2013.06.14

PubMed Abstract | CrossRef Full Text | Google Scholar

Tewawong, N., Suwannakarn, K., Prachayangprecha, S., Korkong, S., Vichiwattana, P., Vongpunsawad, S., et al. (2015). Molecular epidemiology and phylogenetic analyses of influenza B virus in Thailand during 2010 to 2014. PLoS One 10:e0116302. doi: 10.1371/journal.pone.0116302

PubMed Abstract | CrossRef Full Text | Google Scholar

Tong, S., Li, Y., Rivailler, P., Conrardy, C., Castillo, D. A., Chen, L. M., et al. (2012). A distinct lineage of influenza A virus from bats. Proc. Natl. Acad. Sci. U.S.A. 109, 4269–4274. doi: 10.1073/pnas.1116200109

PubMed Abstract | CrossRef Full Text | Google Scholar

Tong, S., Zhu, X., Li, Y., Shi, M., Zhang, J., Bourgeois, M., et al. (2013). New world bats harbor diverse influenza A viruses. PLoS Pathog. 9:e1003657. doi: 10.1371/journal.ppat.1003657

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Q., Cheng, F., Lu, M., Tian, X., and Ma, J. (2008). Crystal structure of unliganded influenza B virus hemagglutinin. J. Virol. 82, 3011–3020. doi: 10.1128/jvi.02477-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Waterhouse, A., Bertoni, M., Bienert, S., Studer, G., Tauriello, G., Gumienny, R., et al. (2018). SWISS-MODEL: homology modelling of protein structures and complexes. Nucleic Acids Res. 46, W296–W303. doi: 10.1093/nar/gky427

PubMed Abstract | CrossRef Full Text | Google Scholar

White, M. C., and Lowen, A. C. (2018). Implications of segment mismatch for influenza A virus evolution. J. Gen. Virol. 99, 3–16. doi: 10.1099/jgv.0.000989

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J., Lau, Y. C., Wu, P., Feng, L., Wang, X., Chen, T., et al. (2018). Variation in influenza B virus epidemiology by Lineage, China. Emerg. Infect. Dis. 24, 1536–1540. doi: 10.3201/eid2408.180063

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, H., Feng, L., Viboud, C. G., Shay, D. K., Jiang, Y., Zhou, H., et al. (2013). Regional variation in mortality impact of the 2009 A(H1N1) influenza pandemic in China. Influenza Other Respir. Viruses 7, 1350–1360. doi: 10.1111/irv.12121

PubMed Abstract | CrossRef Full Text | Google Scholar

Zens, K. D., Chen, J. K., and Farber, D. L. (2016). Vaccine-generated lung tissue-resident memory T cells provide heterosubtypic protection to influenza infection. JCI Insight 1:e85832. doi: 10.1172/jci.insight.85832

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Kang, J., Dehghan, S., Sridhar, S., Lau, S., Ou, J., et al. (2019). A survey of recent adenoviral respiratory pathogens in Hong Kong reveals emergent and recombinant human adenovirus type 4 (HAdV-E4) circulating in civilian populations. Viruses 11:129. doi: 10.3390/v11020129

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Q., Jing, S., Cheng, Z., Yu, Z., Dehghan, S., Shamsaddini, A., et al. (2017). Comparative genomic analysis of two emergent human adenovirus type 14 respiratory pathogen isolates in China reveals similar yet divergent genomes. Emerg. Microb. Infect. 6:e92. doi: 10.1038/emi.2017.78

PubMed Abstract | CrossRef Full Text | Google Scholar