Abstract
Objective
Influenza season occurs every year in China, but its presentation was unusual in the period from December 2017 to early 2018. During this period, influenza activity was increasing across the country and was much greater than during the same period in previous years, with great harm to people’s health.
Methods
In this study, we isolated two human influenza virus strains—A/Hebei/F076/2018(H1N1) and B/Hebei/16275B/2018—from patients with severe influenza in Hebei, China, during the flu season in January 2018, and explored their genetic characteristics, pathogenicity, and transmissibility.
Results
A/Hebei/F076/2018(H1N1) belongs to the human-like H1N1 influenza virus lineage, whereas B/Hebei/16275B/2018 belongs to the Victoria lineage and is closely related to the World Health Organization reference strain B/Brisbane/60/2008. Pathogenicity tests revealed that A/Hebei/F076/2018(H1N1) replicated much more strongly in mice, with mice exhibiting 40% mortality, whereas B/Hebei/16275B/2018 was not lethal. Both viruses could be transmitted through direct contact and by the aerosol route between guinea pigs, but the H1N1 strain exhibited higher airborne transmissibility.
Conclusions
These results may contribute to the monitoring of influenza mutation and the prevention of an influenza outbreak.
Keywords
Introduction
Influenza viruses can cause seasonal or pandemic influenza in humans and animals. According to the different antigenicity of their nuclear protein (NP) and matrix protein (M1), they can be separated into four types: A, B, C, and D.1,2 Of these, influenza A virus (IAV) is widespread and can infect humans and a variety of animals. It is the most common type of influenza in human and livestock. 3 Influenza B virus (IBV) has strong host specificity and is currently reported to infect humans and seals.4,5 Influenza C viruses mainly infect humans, pigs, and cattle, 6 and influenza D viruses mainly infect pigs and cattle. 7
Influenza A viruses are currently categorized into 18 hemagglutinin (HA) and 11 neuraminidase (NA) subtypes, but only viruses with HA subtypes H1, H2, and H3, and NA subtypes N1 and N2 are known to cause influenza pandemics in humans.8–10 Since the beginning of the 20th century, there have been five major outbreaks of pandemic influenza worldwide.
11
H1N1pdm09 was a recent human influenza pandemic that caused outbreaks in 214 countries (or regions) and the deaths of at least 18,000 people worldwide.
12
At present, the H1N1pdm09 influenza virus has replaced the H1N1 subtype seasonal influenza virus that dominated before 2009, with H1N1pdm09 spreading widely among the population in the form of new viruses causing seasonal influenza outbreaks every year.13–15 Sequencing and analysis of the genome of the H1N1pdm09 influenza virus showed that it is a triple reassortment virus that carried genes from human, avian, and swine influenza viruses.16,17 Genetic evolution analysis revealed that the
IBV was first obtained from sick children in 1940 and successfully isolated from seals in 2000.4,20 According to differences in the nucleotide sequence of the HA fragment of the virus, IBV is divided into two major lineages; namely, the B-Victoria branch, which is represented by the B/Victoria/2/87 strain, and the B-Yamagata branch, which is represented by the B/Yamagata/16/88 strain. 21 IBV does not pose a pandemic threat, but localized outbreaks cannot be ignored.
Influenza outbreaks occur almost every year, causing respiratory diseases of varying severity in the population and posing a serious threat to global public health.22–24 Virus monitoring is an effective early warning method against influenza virus pandemics. From December 2017 to early 2018, the morbidity and mortality of influenza were much higher than those in previous years (Figure 1), sounding the alarm for a potential epidemic. To understand the characteristics of the strains involved in this outbreak, our laboratory collected 1570 nasopharyngeal swabs from patients with suspected influenza infection in Hebei, China, during the influenza season in January 2018, from which we isolated 17 influenza virus epidemic strains, including 15 IAV (14 H1N1 and 1 H3N2) strains and two IBV strains. In particular, two strains, A/Hebei/F076/2018(H1N1) and B/Hebei/16275B/2018, were isolated from severely ill patients, whereas the others were isolated from mildly symptomatic patients. In this study, we explored the genetic characteristics, pathogenicity, and transmissibility of these two strains and compared them with a previously isolated seasonal influenza A/Hebei/HB17/2017(H1N1) and a H1N1pdm09 strain A/CA/04/2009(H1N1). The results will aid in influenza pandemic preparedness efforts.

Morbidity (blue bars) and mortality (orange curve) of influenza in China 2008 to 2018. Source: National Health Commission of the People’s Republic of China.
Materials and methods
Ethics statement
All animals were adequately cared for, and the animal studies were conducted in strict accordance with the guidelines of animal welfare of the World Organization for Animal Health. 25 Experimental protocols involving animals were approved by the Animal Care and Use Committee of Military Veterinary Institute (Changchun, Jilin, China; approval number: SCXK 20160008; approval date: 7 March 2016). All experiments with the influenza A (H1N1) and influenza B viruses were performed in biosafety level three laboratories approved by the Academy of Military Medical Sciences. This study did not involve human participants.
Viruses
Two human influenza viruses, isolated from severely ill influenza patients in 2018, were used in this study. The viruses were A/Hebei/F076/2018(H1N1) (hereafter F076; GenBank accession numbers MH748636 and MH748645) and B/Hebei/16275B/2018 ( hereafter 16275B; GenBank accession number MH748675). Additionally, A/Hebei/HB17/2017(H1N1) (hereafter HB17), isolated from severely ill influenza patients in 2017, and A/CA/04/2009(H1N1) (hereafter CA04), a pandemic flu virus, were used as control group viruses in this study. Viruses were grown in 9-day-old specific-pathogen-free eggs (Merial Vital Laboratory Animal Technology Company, Beijing, China) and stored at −80°C.
Phylogenetic and sequence analyses
Viral RNA was extracted from allantoic fluid using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) and reverse transcribed into cDNAs using the primer Uni12 (5′-AGCRAAAGCAGG-3′) or Uni9 (5′-AGCAGAAGC-3′). PCR was performed using a previously described method.26,27 The PCR products were subjected to electrophoresis on 1% agarose gels stained with ethidium bromide and visualized with an ultraviolet light transilluminator. The PCR products were purified and sequenced by Comate Bioscience Company Limited (Changchun, Jilin, China). All sequence data were analyzed using the SEQMAN program (DNAStar Inc., Madison, WI, USA). All reference sequences used in this study were obtained from the National Center for Biotechnology Information GenBank database. Phylogenetic analysis was performed by the distance-based neighbor-joining method using MEGA 7 software (DNAStar Inc.).
Receptor-binding specificity assay
The receptor-binding specificities of human influenza viruses were determined by HA assays with 1% chicken and sheep red blood cell (cRBC and sRBC) suspensions. For sialidase treatment, 90 μL of a 10% cRBC suspension was treated with 10 μL of α-2,3-sialidase (50 mU/μL; TaKaRa, Dalian, China) for 10 minutes at 37°C. The sample was then washed twice with PBS, centrifuged at 210 ×
Mouse experiments
Groups of five 6-week-old female BALB/c mice (Merial Vital Laboratory Animal Technology Company) were anesthetized with ether and intranasally inoculated with 50 µL of a 106.0 EID50 (50% embryo infectious dose) solution of strain F076, HB17, CA04, or 16275B.28,29 The mice were monitored for weight loss and mortality daily for 14 days. Mice that lost >30% of their original body weight were humanely euthanized. All animals were euthanized by inhaling excess CO2.
28
To assess the growth characteristics of the two viruses and the pathological changes in the lungs of infected mice, four groups of 20 mice were anesthetized with ether and intranasally instilled with 106.0 EID50 of either F076, HB17, CA04, or 16275B, and another 3 mice were intranasally instilled with PBS as controls. Three mice were euthanized at 3, 5, and 7 days post-infection (dpi). The lungs, hearts, livers, spleens, kidneys, and brains of infected mice were removed to determine the virus titers. Briefly, the tissues were weighed, and 0.1 g of each tissue was placed into 1 mL of PBS containing 100 U/mL penicillin, generating 10% (wt/vol) tissue homogenates. The tissue samples were homogenized by Tissue Lyser (Qiagen, Hilden, Germany) and centrifuged at 1680 ×
Guinea pig experiments
Hartley strain albino female guinea pigs weighing 300 to 350 g (Merial Vital Laboratory Animal Technology Company) were used in this study. In the transmission studies, three guinea pigs per group were intranasally inoculated with 200 μL of the test viruses at 106 EID50 and housed in a cage placed inside an isolator. The next day, three naïve guinea pigs were cohoused (in the same cage) with the three infected guinea pigs to study direct-contact transmission, and another three naïve guinea pigs per group were housed in a wire-frame cage adjacent to the infected guinea pigs to study aerosol transmission. The distance between the infected and aerosol-contact guinea pig cages was 5 cm. To monitor virus shedding, nasal washes were collected and titrated from all animals at 2, 4, 6, and 8 dpi.
Statistical analysis
Significant differences were identified using one-way analysis of variance with GraphPad Prism software (GraphPad Inc., San Diego, CA, USA). All assays were run in triplicate and are representative of at least 3 separate experiments. The error bars represent the standard deviation.
Results
Phylogenetic analysis of surface genes
Full-length sequences of the influenza virus A/Hebei/F076/2018(H1N1) and B/Hebei/16275B/2018 isolates were compared with those of known influenza viruses in the GenBank database. The eight gene segments of F076 and 16275B were found to have the highest nucleotide sequence homologies predominantly with the epidemic strain (H1N1 and influenza B) from recent years, which are available from the influenza sequence database (https://www.ncbi.nlm.nih.gov/genomes/FLU/Database/nph-select.cgi?go=database; Table 1). Phylogenetic trees based on the HA and NA genes of F076 are shown in Figure 2a and 2b. The H1N1 strain in this study (indicated by solid circles) is divided into the human-like H1N1 influenza virus lineage and is closely related to the epidemic strain of recent years. The HA phylogenetic tree for 16275B is shown in Figure 2c; it consists of epidemic strains and WHO reference epidemic strains (marked with solid triangles), as well as reference strains of the Yamagata and Victoria lineages (marked with solid squares). 19 The results showed that B/Hebei/16275B/2018 (marked with solid circles) is a member of the Victoria lineage and is closely related to the reference strain B/Brisbane/60/2008.
Sequence homology of human influenza strains F076 (A/Hebei/F076/2018) and 16275B (B/Hebei/16275B/2018), both associated with severe illness in 2018, with known influenza viruses (https://www.ncbi.nlm.nih.gov/genomes/FLU/Database/nph-select.cgi?go=database).

Phylogenetic analysis of the hemagglutinin (HA) and neuraminidase (NA) genes of strain F076 (A/Hebei/F076/2018, human, associated with severe illness in 2018; panels a and b, respectively) and the HA genes of strain 16275B (B/Hebei/16275B/2018, human, associated with severe illness in 2018; panel c). Phylogenetic trees of the HA and NA genes were constructed by the distance-based neighbor-joining method with 1000 bootstrap replicates using MEGA 7 software (DNAStar Inc., Madison, WI, USA). Horizontal distances are proportional to genetic distances. Influenza virus isolates in this study are marked with solid circles. The WHO reference strains are marked with solid triangles. Reference strains of the Yamagata and Victoria lineage are marked with solid squares.
F076 and 16275B exhibited comparable binding affinities for avian and human sialic acid receptors
We evaluated the receptor-binding specificities of F076 and 16275B using an HA assay, and compared them with HB17 and CA04. The surface of cRBCs contains α-2,3-linked and α-2,6-linked sialic acid receptors. cRBCs treated with α-2,3-sialidase contain only α-2,6-linked sialic acid receptors, cRBCs treated with VCNA contain no receptors, and the surface of sRBCs contain only α-2,3-linked sialic acid receptors. As expected, all viruses could agglutinate untreated cRBCs but not VCNA-treated cRBCs or sRBCs, as shown in Figure 3. The HA titers are shown in Figure 3. These results indicated that the F076, HB17, and 16275B strains maintained affinity for human-like (α-2,6) receptors but exhibited no affinity for avian-like (α-2,3) receptors, suggesting that they are less likely to infect avian species.

Agglutination activity (hemagglutinin [HA] titers) of four influenza strains in red blood cells (RBCs) of different types: (a) chicken RBCs (with α-2,3-linked sialic acid receptors and α-2,6-linked sialic acid receptors); (b) chicken RBCs treated with α-2,3-sialidase (with only α-2,6-linked sialic acid receptors); (c) sheep RBCs (with only α-2,3-linked sialic acid receptors); and (d) chicken RBCs treated with
Pathogenicity of F076 and 16275B in mice
The pathogenicity of F076 and 16275B viruses was evaluated in a mouse model and compared with that of HB17 and CA04. Groups of five 6-week-old female BALB/c mice were used in this experiment. After 24 hours of infection with F076, HB17, CA04, or 16275B, the mice showed clinical manifestations such as coat inversion and loss of appetite. Mice inoculated with viruses rapidly lost approximately 15% of their weight (Figure 4a). Mice infected with F076 exhibited 40% mortality, mice infected with HB17 exhibited 20% mortality, whereas the CA04 and 16275B infections were not lethal (Figure 4b). The change in body weight of mice infected with F076 was more obvious (

Pathogenicity of the isolated viruses in mice. Five mice per group were intranasally inoculated with the viruses at 106 EID50 (50% embryo infectious dose). (a) Body weights during 14 days of monitoring. Values represent the average score of overall body weight loss compared with the initial body weight ± standard deviation (SD). The change in body weight of mice infected with F076 was more obvious (

Histopathological analysis. At 3 days post inoculation (dpi), lungs were collected from mice inoculated with 106.0 EID50 (50% embryo infectious dose) of F076, 16275B, HB17, or CA04, and were fixed with formalin, embedded in paraffin, and stained with hematoxylin and eosin. The images were obtained at a magnification of 200×. Panels a–e show pathological changes in the lungs of mice inoculated with PBS, F076, 16275B, HB17, and CA04, respectively. Arrow a indicates alveolar wall thickening, accompanied by lymphocyte infiltration, neutrophil infiltration; arrow b indicates pulmonary edema, a small amount of eosinophilic serous substance exuding in the alveolar cavity; arrow c indicates local bronchus with necrosis and shedding of epithelial cells and nuclear fragmentation; arrow d indicates alveolar wall thickening, accompanied by a large amount of lymphocytes and neutrophil infiltration; arrow e indicates bleeding. The scale bar represents 100 µm. Influenza strains: F076 (A/Hebei/F076/2018, human, associated with severe illness in 2018), HB17 (A/Hebei/HB17/2017, human, associated with severe illness in 2017; control), CA04 (A/CA/04/2009, pandemic flu virus; control), and 16275B (B/Hebei/16275B/2018, human, associated with severe illness in 2018).
Transmissibility of F076 and 16275B in guinea pigs
Hartley strain albino female guinea pigs weighing 300 to 350 g were used to investigate the contact or airborne transmissibility of the two viruses compared with control strains. Three guinea pigs were intranasally inoculated with 200 μL of the test virus at 106 EID50/mL and housed in a cage placed inside an isolator. Infection was detected by a nasal wash of the influenza A virus (F076) infection group, which had the highest viral titer of 104.2 EID50/mL. F076 was transmitted to two guinea pigs via the direct contact route and to one guinea pig via the aerosol route (three animals were examined in each group; Figure 6a). The highest viral titers in the F076 contact and aerosol infection groups were 102.8 EID50/mL. HB17 was transmitted to two guinea pigs via the direct contact route and to one guinea pig via the airborne route (three animals were examined in each group; Figure 6c). The highest viral titers in the HB17 contact and aerosol infection group were 102.0 EID50/mL and 101.3 EID50/mL, respectively. CA04 was transmitted to three guinea pigs via the direct contact route and to three guinea pigs via the airborne route (three animals were examined in each group; Figure 6d). The highest viral titers in the CA04 contact and aerosol infection group were 103.0 EID50/mL. 16275B was transmitted to three guinea pigs via the direct contact route and to one guinea pig via the airborne route (three animals were examined in each group; Figure 6b). The highest viral titers in the 16275B contact and aerosol infection group were 101.8 EID50/mL. These results suggest that both influenza A and B viruses included in this study could be transmitted via direct contact and by the airborne route.

Horizontal transmission of influenza viruses between guinea pigs. Groups of three guinea pigs were inoculated with the indicated viruses at 106 EID50 (50% embryo infectious dose): (a) F076, (b) 16275B, (c) HB17, and (d) CA04. The next day, the inoculated animals were paired by cohousing with direct-contact guinea pigs; aerosol-transmission animals were also housed in a wire-frame cage adjacent to the infected guinea pigs. Nasal washes were collected from all animals for detection of virus shedding every other day beginning on day 2 after initial infection. Each color bar represents the virus titer in an individual animal. Dashed lines indicate the lower limit of virus detection. Influenza strains: F076 (A/Hebei/F076/2018, human, associated with severe illness in 2018), HB17 (A/Hebei/HB17/2017, human, associated with severe illness in 2017; control), CA04 (A/CA/04/2009, pandemic flu virus; control), and 16275B (B/Hebei/16275B/2018, human, associated with severe illness in 2018).
Discussion
Influenza epidemics and pandemics pose a severe threat to public health and the global economy. Long-term viral surveillance can reflect an outbreak situation in a timely manner and provide early warning of an epidemic or pandemic. In 2018, the proportion of outpatient and emergency influenza-like cases and the positive rate of influenza virus tests reported by the influenza surveillance sentinel hospital in China were higher than those reported during the same period in the previous 3 years. The number of severe cases was also higher.30–32 Therefore, there was an urgent need to study the pathogenicity and transmissibility of the major influenza viruses circulating in 2018.
In this study, we isolated two influenza strains from severely ill patients in Hebei, China, during the 2018 flu season. Phylogenetic analyses revealed that F076 belongs to the human-like H1N1 influenza virus lineage and 16275B belongs to the Victoria lineage. Both are closely related to the epidemic strain present in recent years. In addition, 16275B is closely related to the WHO reference epidemic strain B/Brisbane/60/2008, which in 2009 was recommended by the WHO for inclusion as a vaccine strain.
33
We compared the gene sequences of influenza epidemic strains in recent years and found that F076 has the
The infection caused by influenza viruses begins with the HA glycoprotein binding to sialic acid receptors on the surface of host cells.38,39 Therefore, the receptor-binding specificity of influenza viruses is a key determinant of cross-species transmission.39,40 Avian influenza viruses preferentially bind to α-2,3 sialic acid receptors, whereas mammalian influenza viruses preferentially bind to α-2,6 sialic acid receptors. In this study, both the F076 strain and 16275B strain showed specificity only for α-2,6 sialic acid receptor binding. Therefore, these 2 viruses preferentially infect mammals and are less likely to infect poultry.
In our mouse study, F076 exhibited much higher growth properties and pathogenicity than 16275B, which might be caused by the weak binding affinity for the sialic acid receptor of IBV compared with IAV.41–46 In addition, F076 caused lethal pneumonia in mice without prior adaptation (fatality rate of 40%), demonstrating higher pathogenicity than previous seasonal influenza viruses (Figure 4),
47
which might be associated with the
Airborne transmission capacity is an important indicator of whether an influenza strain can cause an epidemic or a pandemic.51,52 The H5N1, H5N6, and H7N9 strains can infect humans but show limited airborne transmissibility.53,54 Therefore, they are less likely to cause an epidemic. In this study, both F076 and 16275B were transmissible by airborne or direct contact in guinea pigs, but their transmission was not as extensive as that of H1N1pdm09 (Figure 6), indicating that they represent a potential localized outbreak risk.
This study has two limitations. First, the two strains isolated in this study do not and cannot reflect the overall situation of this epidemic. More strains should be characterized in future studies. Second, variations in the virus as they spread from person to person were not considered in this study.
Overall, we assessed the genetic characteristics, pathogenicity, and transmissibility of two influenza strains isolated from severely ill patients during the 2018 flu season. Our results increase our knowledge of current circulating strains and contribute to the prevention and control of influenza epidemic. Continuous monitoring and assessment are needed in the future.
Footnotes
Authors’ contributions
ZDG, CZ, JXL, and LGC conceived and designed the experiments; HC, CZ, JXL, SSD, CMZ, and JXL performed the experiments; and ZDG, ZYW, and CZ wrote the manuscript. All authors approved the content of the final article.
Declaration of conflicting interest
The authors declare that there is no conflict of interest.
Funding
This work was supported by National Key Research and Development Program of China (grant no. 2016YFD0501000). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
