Mouse adaptation of influenza B virus increases replication in the upper respiratory tract and results in droplet transmissibility in ferrets.

To investigate the molecular changes that allow influenza B viruses to adapt to new mammalian hosts, influenza B/Florida/04/2006 was serially passaged in BALB/c mice until highly virulent. The viral factors underlying this transition were then investigated in mice and ferrets. Five viruses, including the wild-type virus (P0), three intermediate viruses (P5, P9, and P12), and a lethal mouse-adapted virus (P17 (MA)), harbored one to five amino acid substitutions in the hemagglutinin, M, NP, and PA segments suggesting that these mutations enhance virulence. The P17 (MA) virus replicated significantly more efficiently than the P0 virus both in vitro and in vivo (P < 0.0001), and was highly virulent (MLD50: 10(5.25)TCID50) while the P0, P5, and P9 viruses did not kill any infected mice (MLD50 > 10(6.0)TCID50). Furthermore, the P17 (MA) virus grew to greater titers in the ferret upper respiratory tract compared with the P0 and intermediate viruses, and only the P17 (MA) virus was transmissible between ferrets via both direct and aerosol contact. To our knowledge, this is the first study to demonstrate ferret-to-ferret transmission of influenza B virus and to delineate factors that may affect its transmission.

Influenza A and B viruses cause global outbreaks of mild-to-severe respiratory disease each year and are among the most significant public health problems worldwide. The influenza A viruses can be divided into subtypes depending on the antigenicity of their two surface glycoproteins, hemagglutinin (HA) and neuraminidase (NA)1. In contrast, the influenza B virus is not divided into subtypes and has only two genetically distinct lineages in humans, B/Vitoria/2/87-like (Victoria lineage) and B/Yamagata/16/88-like (Yamagata lineage). The B/Victoria lineage was predominant during the 1980s, while the B/Yamagata lineage was predominant throughout most of the world during the 1990s2. However, a recent study revealed that one of each influenza B linage viruses has been predominantly circulating approximately once every three years with genetic shift and antigenic drift3. While influenza B viruses are believed to cause milder clinical disease in mammalian hosts than influenza A viruses, these viruses are essential components of influenza vaccines administered to susceptible groups, such as children and young adults4. In most years, the vaccination process is successful; however, antigenic drift and mismatches between circulating viruses and selected vaccine strains have occurred with significant consequences5, and little cross-protection between lineages is considered to exist4. In this regard, during the 2007–2008 flu seasons, the selected vaccine strains protected against infection with the B/Victoria lineage, but over 95% of the circulating influenza B viruses were of the B/Yamagata lineage6. Thus, one proposed alternative to current vaccine strategies is the production of a quadrivalent influenza vaccine (QIV) that includes both influenza B virus lineages7.

Previous reports of influenza B virus isolation from dogs, a pig, and harbor seals have not met established criteria for determining whether an infection has occurred, and the virus has shown little-to-no adaptation to these hosts. In fact, influenza B viruses isolated from non-human species have been found to be closely related to concurrently circulating human strains8910. A natural reservoir of influenza B viruses outside humans has not been identified and thus, these viruses are considered to be exclusively human pathogens. Although both influenza A and B viruses have caused widespread human epidemics in the past, the genetic determinants of the virulence and transmissibility of influenza B viruses have not been studied in depth. In fact, the limited known host range of the influenza B virus and the paucity of available animal models have precluded investigation of influenza B pathogenicity and transmission factors, as well as the evaluation of antiviral compounds and vaccine efficacies in animals. In general, mice have been considered a good mammalian model for investigating pathogenic mechanisms and host range determinants of influenza viruses11. In this regard, mouse adaptation via serial passage in the lungs has often been used to study influenza A virulence factors in mice1213141516. In addition, ferrets are a suitable small animal model with which to study influenza infection because they manifest flu-like symptoms and immune responses similar to those observed in humans1718192021. While ferrets are used routinely to study the transmission of influenza A viruses2223242526, prior to this study, their utility for influenza B virus transmission studies had not been established. Limited previous studies have shown that ferrets infected with B/Victoria-like and B/Yamagata-like viruses display mild clinical signs of infection and have lower elevations in core body temperature than do animals infected with influenza A viruses27; however, ferret-to-ferret transmission of influenza B viruses has not been evaluated to date.

In this study, we provide insight into the mammalian adaptation of influenza B viruses through molecular modifications that increase pathogenicity. Furthermore, we report that mouse-adapted influenza B viruses not only replicate in the upper respiratory tract of ferrets, but are also transmitted efficiently from infected to naïve animals by both direct and aerosol contact. Our results suggest that mutations in multiple genes encoding different viral proteins are associated with the growth and virulence of influenza B viruses.

Results

Genetic characterization of mouse-adapted B/Florida/04/06 influenza virus

To generate the mouse-adapted (MA) influenza B strain, B/Florida/04/06 was serially passaged in the lungs every three dpi until high virulence with 100% mortality was observed in BALB/c mice. After seventeen passages, infection with the B/Florida/ma04/06 virus resulted in severe influenza-like symptoms with 100% mortality in mice at 7 dpi. To investigate the molecular changes that occurred during the mouse adaptation of the virus, the P0 (B/Florida/04/05), three intermediate-passage viruses (P5, P9, and P12), and the seventeenth passage virus P17 (MA) were isolated by plaque purification, and the full-length sequences of each virus were analyzed. Compared with the P0 strain, the P17 (MA) strain possessed a total of seven nucleotide changes in four of the eight gene segments and five of the ten genes; these changes resulted in five amino acid mutations (Table 1). These amino acid mutations developed sequentially during the mouse adaptation process (Table 2). The first mutation that was observed was in the HA gene (D424G) in the P5 (MAHA) viruses (three of eight plaques). The second and third mutations occurred simultaneously in M1N240T and BM2N14S in the P9 (MAM) viruses, while six of eight plaques lost the HA gene (D424G) mutation. Increased virulence with 80% mortality was observed for P12 (MAHA/M), and full-length sequencing revealed that the P12 (MAHA/M) virus had all three of the mutations that were detected in the P5 and P9 viruses (HAD424G, M1N240T, and BM2 N14S). Notably, the fully mouse-adapted P17 (MA) virus expressed the three aforementioned mutations together with an additional two mutations in the NPK294N and PAV609I proteins. Together, these results demonstrate that mutations occurred sequentially during the mouse lung-to-lung passage of the B/Florida/04/06 virus and suggest that multiple genetic mutations might be required for the acquisition of high virulence by influenza B in mice.

Table 1

Nucleotide and amino acid changes in wild-type and mouse-adapted influenza virus B/Florida/04/2006.

Genesegment	Changed nucleotidesa			Changed amino acidsb
	Position	Original	Mouse- adapted	Position	Original	Mouse-adapted
HA	951	T	A	317	P	—
	1271	A	G	424	D	G
M1c	719	A	C	240	N	T
BM2c	27	C	T	9	I	—
NP	40	A	G	14	N	S
	882	G	C	294	K	N
PA	1825	G	A	609	V	I

aNumbering of nucleotides is relative to B/Lee/40 and begins at the start of the non-coding region of the RNA segment.

bNumbering of amino acids is relative to B/Lee/40 and begins at the start of the coding region.

cRNA segment 7 encodes two polypeptides, M1 and BM2, using tandem cistrons.

Table 2

Sequences of amino acid changes during mouse adaptation.

Viral proteins(amino acid)	P0	P5	P9	P12	P17
HA (424)	Asp	Gly	—	Gly	Gly
M1 (240)	Asn	—	Thr	Thr	Thr
M2 (14)	Asn	—	Ser	Ser	Ser
NP (294)	Lys	—	—	—	Asn
PA (609)	Val	—	—	—	Ile

Note: P0, parent virus B/Florida/04/2006; P5, 5th mouse passage; P9, 9th mouse passage; P12, 12th mouse passage; and MA, mouse-adapted virus after 17 passages.

Growth characteristics of mouse-adapted influenza B viruses in vitro

To investigate whether sequential mutations of the B/Florida/04/06 virus affect its growth kinetics, each plaque-purified representative mutant, P0, P5 (MAHA), P9 (MAM), P12 (MAHA/M), and P17(MA), was inoculated onto MDCK cell monolayers, and the growth characteristics at 12, 24, 48, and 72 hours post infection (hpi) were assessed. The P0 virus reached a peak titer (5.0 log10 TCID50/mL) at 48 hpi, while the P17 MA virus showed a significantly higher peak viral titer (6.25 log10 TCID50/mL) than did the P0 virus (P < 0.001) at 48 hpi (Fig. 1A). Although the P5 (MAHA) and P9 (MAM) viruses did not display significantly higher peak titers in infected MDCK cells compared with the P0 virus, the P12 (MAHA/M) virus displayed significantly greater peak viral titers at 48 hpi (P < 0.001; Fig. 1A). These results clearly demonstrate that the serially passaged, mouse-adapted viruses have a considerable growth advantage in vitro.

Figure 1

Replication kinetics of mouse-adapted virus.

(A) Cells were infected with four influenza B viruses and mouse-adapted viruses at a multiplicity of infection55 of 0.001. Cell culture supernatants were collected at 12, 24, 48, and 72 hours post infection, and viral titers were measured in units of log10 TCID50/mL. (B) Determination of plaque phenotype. *indicates P < 0.0001, and †indicates P < 0.001 compared with the P0 (B/Florida/04/06). The error bars show the standard errors of the mean (SEMs), which were derived from three independent experiments.

Next, we compared the phenotypes of plaques formed in MDCK cells for each mouse-adapted virus. The P0 virus formed homogeneous, small, punctiform plaques (mean size, 1.0 ± 0.13 mm). During the course of mouse adaptation, the plaque sizes of each intermediate virus (P5, P9, and P12) became larger and more heterogeneous in size, ranging from ~1.5 mm to ~2.3 mm. Strikingly, the P17 (MA) virus produced large, homogeneous, circular plaques with an average size of 3.5 ± 0.11 mm (Fig. 1B). These results indicate that the P17 (MA) virus was homogeneous and readily replicated in MDCK cells.

Pathogenicity of mouse-adapted influenza B viruses in vivo

To determine the pathogenicities of the mouse-adapted influenza B viruses in mice, we infected ten mice per group with each MA influenza B strain and monitored survival and weight loss for 14 dpi. The three cohorts of mice infected with P0 (parental virus), P5 (MAHA), and P9 (MAM) demonstrated 100% survival for 14 dpi, a very modest body weight loss (<1%), and no clinical symptoms (Fig. 2A). Mice infected with the P12 (MAHA/M) virus exhibited up to a 20% reduction in body weight by 8 dpi, and only two mice survived the infection (20% survival rate; Fig. 2B). In contrast, all P17 (MA)-infected mice became severely ill, lost body weight starting at 2 dpi (>25% from baseline), and eventually died or were humanly euthanized by 7 dpi (0% survival). These results suggest that mutations in the P12 (MAHA/M) and P17 (MA) viruses contribute to viral pathogenicity in mice.

Figure 2

Pathogenicities of mouse-adapted influenza B viruses in mice.

All groups of mice were observed for 14 days for (A) morbidity and (B) mortality. (C) Virus growth in mouse lungs. Five groups were infected intranasally with 6.0 log10 TCID50 of each virus. Lung virus titers on days 1, 3, 5, and 7 dpi are expressed as log10 TCID50 per gram of tissue. *indicates P < 0.0001, and † indicates P < 0.001 compared with the P0 (B/Florida/04/06); mouse death is indicated by ‡. The lower limit of detection (1.5 log10 TCID50) is represented by the dotted line.

We next assessed the viral replication kinetics in mouse lungs for each MA influenza B mutant to determine the association between lung viral titers and virulence. The parental virus (P0) could be recovered from mouse lungs until 5 dpi and reached peak titers (2.5 log10 TCID50/g) at 3 dpi. Although the P5 (MAHA) virus showed higher lung viral titers (3.5 log10 TCID50/g) than did P9 (MAM) at 1 dpi, the P9 (MAM) virus displayed greater viral titers compared with P5 (MAHA) at 3 and 5 dpi (Fig. 2C). However, the P12 (MAHA/M) virus had higher lung viral titers compared with those in the P0, P5 (MAHA), and P9 (MAM) groups until 7 dpi (5.25 log10 TCID50/mL peak titer). Finally, the P17 (MA) mutant-infected mice showed the highest viral lung titers until 5 dpi (5.5 TCID50/mL peak titer), and all mice succumbed to infection by 7 dpi (Fig. 2). These data suggest that the P12 (MAHA/M) and P17 (MA) viruses replicate efficiently in mouse lungs and that the P17 (MA) virus is better adapted to mouse lungs than the parental and intermediate viruses, as evidenced by its improved growth kinetics.

Analysis of proinflammatory cytokine expression levels

The dysregulation of proinflammatory cytokines, which is sometimes referred to as a cytokine storm, is associated with the pathogenesis of severe influenza virus infections3536. Therefore, we investigated whether the infection of mice with mouse-adapted viruses results in cytokine responses related to pathogenicity in influenza B virus infections. As shown in Fig. 3, levels of the proinflammatory cytokines TNF-α, IL-4, IL-5 and IL-6 were higher in BALF from the P12 (MAHA/M) and P17 (MA) virus-infected mice relative to mice infected with the other viruses from 3dpi (p < 0.0001) through 5 dpi (p < 0.001; Fig. 3A–D). The levels of IL-12p70 were also elevated in the P12- (MAHA/M) and P17 (MA)-infected mice from 1dpi to 3 dpi compared with the P0 virus-infected mice; however, by 5dpi the levels of IL-12p70 decreased to levels comparable to those observed for the other groups (Fig. 3E). Only the P12 (MAHA/M)-infected group showed a significant upregulation of IFN-γ compared with the P0 virus-infected mice from 3 dpi (p < 0.05) to 5dpi (p < 0.001; Fig. 3F). These data suggest that the virulence of each mouse-adapted virus was closely correlated with the elevated proinflammatory cytokine expression levels observed in the lungs of infected mice.

Figure 3

Cytokine expression associated with mouse-adapted influenza B virus.

The mice were euthanized (n = 5 per each time point), and lung BALF was harvested. Cytokine production in lung BALF was analyzed, including the levels of TNF-α, IL-4, IL-5, IL-6, IL-12p70, and IFN-γ, by BioPlex analysis. Statistical significance compared to mice infected with the P0 (B/Florida/04/06) virus was determined by the t-test. *indicates P < 0.0001 and † indicates P < 0.001 compared with the P0 (B/Florida/04/06).

Replication and transmission of mouse-adapted influenza B viruses in ferrets

Although the ferret is the most relevant small-animal model for studying human influenza A virus infections, a preferred transmission model for recent influenza B virus strains has not yet been described. To determine whether the P17 (MA) virus is transmissible as well as pathogenic in a ferret model, groups of ferrets were inoculated with each of virus and 24 h later, these ferrets were caged with or adjacent to naïve ferrets. The parental virus (P0) grew moderately well in the upper respiratory tract of the infected ferrets, reaching a mean peak titer of 4.75 log10 TCID50/mL at 3 dpi, and the virus was recovered in nasal washes until 5 dpi (Fig. 4A). There were no clinical signs of influenza infection in the ferrets infected with P0 (B/Florida/04/06) (<2% weight loss; Fig. 4A,B). The P17 (MA)-infected ferrets exhibited high viral titers (>10 times greater than that observed with P0) at all time points in the upper respiratory tract, with maximum nasal wash titers of 5.5 log10 TCID50/mL observed at 3 dpi. In addition, all ferrets inoculated with the P17 (MA) virus showed mild signs of influenza, marked by infrequent sneezing, a 2–8% reduction in mean body weight, and temperature elevations of up to 2.5 °C between 3 and 5 dpi (Fig. 4D–F).

Figure 4

Weight loss, body temperature, and viral shedding in ferrets inoculated with mouse-adapted influenza virus.

Two groups of ferrets (n = 12/group) were inoculated intranasally with 106 TCID50 of virus. (A–C) P0 (B/Florida/04/06) and (D, E, and F) P17 (MA) virus. (A,D) Viral shedding from the nasal wash samples of inoculated ferrets. (B,E) Ferret weight changes during the mouse-adapted virus transmission experiments. Weight loss of inoculated animals is represented as a percentage of body weight. (C,F) Temperature change in inoculated animals is represented as a °C. After 24 h, each two inoculated animals in each group were individually paired and cohoused with direct contact ferrets (n = 2) and indirect aerosol contact ferrets (n = 2). The daily weights and temperature changes (as percentages) of each ferret were compared with the weights measured prior to the initiation of the experiments. Viral replication was monitored by titrating the viruses in the nasal washes that were collected every other day from both the inoculated (Left) and contact-exposed (Right) ferrets. The limit of detection was 1.5 log10TCID50/mL and is indicated by the dotted line for each representation.

Notably, only the P17 (MA) virus was readily transmitted to naïve ferrets via both direct contact and aerosol contact, as shown by the viral titers at 3dpc (4.0 log10 TCID50/mL at 5dpc; Fig. 4D). In contrast, no virus was recovered from the ferrets exposed either directly or through aerosol contact to P0 virus-infected ferrets (Fig. 4A). Furthermore, both direct and aerosol contact with P17 (MA) virus-infected ferrets resulted in high HI titers (160 HI GMT for direct contact and 40 HI GMT for aerosol contact) at 17 dpc (Table 3). Therefore, the P17 (MA) virus both maintained pathogenicity and exhibited enhanced airborne droplet transmission. A comparison of the viral titers in the ferret respiratory organs between the P0 (parent) and P17 (MA) viruses revealed increased replication of the P17 (MA) virus (up to 5.25 log10 TCID50/mL in the nasal turbinate of the inoculated ferrets), which was detected until 7 dpi. In contrast, infection with the P0 virus resulted in significantly lower viral titers at all time points (P < 0.05; Fig. 5). In addition, only the P17 (MA) virus was recovered from the trachea and lungs of infected ferrets.

Table 3

Seroconversion associated with mouse-adapted influenza B viruses in ferrets exposed to direct and aerosol contact.

Virus	Direct contact	Aerosol contact
	Serum (HI GMT*)†
P0 (Parent)	<10	<10
P5 (MAHA)	<10	<10
P9 (MAM)	<10	<10
P12 (MAHA/M)	40	<10
P17 (MA)	160	40

†The parent virus was used in HI assays; sera were collected at 17 dpc.

*Note: GMT, geometric mean titers.

Figure 5

Mouse-adapted influenza B virus replication in the respiratory organs of ferrets.

Ferrets (n = 2/day) were euthanized for virus titration on days 1, 3, 5, and 7 after infection. (A) P0 (B/Florida/04/06). (B) P17 (MA) virus. The virus titers in the nasal turbinates, trachea, and lung were determined in terms of the TCID50 in MDCK cells. Data are means ± SM of virus titration from three independent experiments. Data are representative of three independent experiments. The lower limit of detection (1.5 log10TCID50) is indicated by the dotted line.

Because the P17 (MA) virus exhibited efficient airborne droplet transmission to ferrets housed in adjacent cages, the remaining three intermediate mutants (P5 (MAHA), P9 (MAM), and P12 (MAHA/M)) were also assessed to determine which gene(s) or mutation(s) play critical roles in this process. We found that the P5 (MAHA) and P9 (MAM) viruses had viral replication properties similar to those of the wild-type P0 (parent) virus and failed to be transmitted either directly or through aerosol contact (Fig. 6A,B). Notably, the P12 (MAHA/M) virus showed increased viral peak titers (4.75 log10 TCID50/mL) compared to those of the wild-type P0 (parent), P5 (MAHA), and P9 (MAM) viruses and was transmitted to ferrets through direct contact, but not through aerosol contact (Fig. 6C). Only direct contact with the P12 (MAHA/M) virus-infected ferrets resulted in a modest serum HI titer (40 HI GMT) on day 17 dpc (Table 3). Moreover, viral titers approximately five to ten times higher (3.25 to 4.75 log10 TCID50/g) were observed in P12 (MAHA/M) virus-infected nasal turbinates compared with those infected with the P5 (MAHA) and P9 (MAM) viruses, and these elevated levels persisted until 5 dpi (Fig. 7). In contrast, no virus was detected in the lungs or trachea, with the exception of infection with the P12 (MAHA/M) virus at 1 dpi, suggesting that the ability of each intermediate mutant to replicate in the respiratory tracts of infected ferrets was impaired compared with the P17 (MA) variant. Taken together, our results indicate that although the clinical symptoms caused by the P17 (MA) virus were not significantly different from those caused by the wild-type P0 (parent) virus, mutations in various genes generated during mouse adaptation collectively contributed to the transmissibility of the influenza B virus in the ferret model. Furthermore, the finding that the P17 (MA) virus replicates efficiently in the upper respiratory tract supports the finding that this virus is transmitted efficiently between ferrets, even by indirect means.

Figure 6

Transmission of mouse-adapted intermediate-passage influenza viruses in ferrets.

Three groups of ferrets (n = 12/group) were inoculated intranasally with 106 TCID50 of the following viruses: (A) P5 (MAHA), (B) P9 (MAM), and (C) P12 (MAHA/M). After 24 h, each two inoculated animals in each group were individually paired and cohoused with direct contact ferrets (n = 2) and indirect aerosol contact ferrets (n = 2). Viral replication was monitored by the titration of nasal washes collected every other day from both the inoculated (Left) and contact (Right) animals. The limit of detection was 1.5 log10TCID50/mL and is indicated by the dotted line.

Figure 7

Mouse-adapted, intermediate-passage influenza virus replication in respiratory organs.

Ferrets (n = 2/day) were euthanized on days 1, 3, 5, and 7 after infection for titration of the following viruses: (A) P5 (MAHA), (B) P9 (MAM), and (C) P12 (MAHA/M) viruses. The virus titers in the nasal turbinates, trachea, and lung were determined by measuring in terms of the TCID50 in MDCK cells. Data are representative of three independent experiments. The lower limit of detection (1.5 log10TCID50) is indicated by the dotted line.

Polymerase activity is enhanced in a mouse-adapted influenza B virus

To determine to what extent the PA and NP mutations were responsible for the high titer of the P17 (MA) virus observed in the upper respiratory tract (turbinate), we assessed polymerase activity using a luciferase reporter assay34. The PA (V609I) and/or NP (K294N) mutations were incorporated into parental virus, individually or simultaneously, using site-directed mutagenesis. At 33 °C, a temperature representative of the upper respiratory tract, polymerase complexes with the individual PA and NP mutants showed approximately 105% and 146% polymerase activity, respectively, compared with that of the parent polymerase complex (Fig. 8). Notably, the simultaneous PA and NP mutants showed the highest polymerase activity (163%) out of all of the viral polymerase complexes tested at 33 °C (P < 0.0001)(Fig. 8). However, the polymerase activity of the simultaneous PA and NP mutant (124%) was lower than that of the single NP (294N) mutant (150%) at 37 °C. Nevertheless, these results indicate that the NP (294N) mutation increases polymerase activity at both 33 °C and 37 °C.

Figure 8

Enhanced polymerase activity of mouse-adapted influenza B virus.

The polymerase activities of reconstituted RNP complexes composed of the PB2, PB1, PA, and NP plasmids of the virus strains P0 (B/Florida/04/06) and P17 (MA). The indicated combinations and single-point mutants are shown. Note that PB2 and PB1 are the same in all complexes because no mutation was observed in this gene. The luciferase activity values are the mean values of at least three assays; (A) 33 °C and (B) 37 °C. *indicates P < 0.0001, and ‡indicates P < 0.01 compared with the parent virus (P0).

Discussion

Viral adaptation to a specific host is considered to be one of the primary mechanisms underlying influenza virus evolution. The resulting zoonotic influenza viruses have the potential to become novel pandemic influenza strains that can spread throughout immunologically-compromised human populations373839. Although the influenza B virus can induce severe disease and even death in human populations, most human influenza B virus strains are avirulent in mouse models, only transiently infecting mice producing low titers in the lungs.

In this study, we selected a recent influenza B vaccine strain (B/Florida/04/06) and generated a mouse-adapted, lethal B/Florida/ma04/06 virus. Of note, full-length sequencing of intermediate strains revealed the mutations that were observed in four genes (the HA (P5 and P12), M1 (P9 and P12), M2 (P9 and P12), NP (P17), and PA (P17) genes; Table 2) accumulated serially. However, increased virulence was observed between the P12 (MAHA/M) and P17 (MA) passages of B/Florida/04/06. A relatively high degree of virulence (80% case fatality rate and an MLD50 of 6.5 TCID50 was first observed in the P12 (HAN424G, M1N240V, and BM2N107S) virus, and the P17 (MA) mutant strain showed 100% mortality with a low MLD50 of 5.5 TCID50. These results suggest that multigenic factors are associated with the influenza B virus virulence and demonstrate that the HA, NP, M, and PA genes play important roles in the virulence of B/Florida/04/06.

Previous studies have focused on the roles of the M1 and NS1 proteins in increased virulence following the adaptation of influenza viruses in mouse lungs32404142. Similarly, McCullers et al. showed that the serial passage of an initially avirulent influenza B virus (B/Memphis/12/97) results in the selection of a lethal mouse-adapted strain, and suggested a critical role for a single amino acid mutation (N221S) in the C-terminal domain of the M1 gene product40. However, no mutation was commonly found between the current study and the previous B/Memphis/12/97 study. In fact, only one similar mutation (N240V) in the C-terminal domain of the M1 gene was acquired as early as P9 in the current study, but in vivo experiments revealed no critical role of this mutation in increased virulence (Fig. 2). The influenza B virus NS1 has been shown to bind to the interferon-stimulated gene 15 protein (ISG15) and subsequently to prevent the ISGylation of target proteins, thereby blocking downstream antiviral effects434445. Of note, we did not detect any NS1 mutations in the B/Florida/04/06 strain until passage 17. Interestingly, we identified a consensus change in the HA gene (HAN424G) in the P5 and P12 viruses. A single HA mutation in the P5 virus was associated with increased viral growth kinetics in infected mouse lungs at early time points, although there was no increase in body weight loss or mortality caused by this virus (Fig. 2). In contrast, a virus with a combination of mutations in M1, BM2, and HA (P12, MAHA/M) caused significantly more body weight loss and mortality and enhanced viral growth kinetics compared with the parental strain. Subsequently, additional NP and PA mutations were found only in the fully mouse-adapted, lethal B/Florida/ma04/06 P17 (MA) strain. These results suggest that the HA mutation may increase the susceptibility of mice to the B/Florida/04/06 virus and rather than a single amino acid mutation, it is the combination of this HA mutation with mutations in internal genes (M1, BM2, NP, and PA genes) that is critical for high virulence in mice.

Inflammatory cytokines and chemokines are produced during influenza virus infections4647, which is the first consequence of the activation of innate immune cells. However, the multiple functions of cytokines can be either beneficial or detrimental to virus-infected hosts. Previous studies have suggested that IL-6 is an early inflammatory marker, and its level has been shown to be well-correlated with infection severity48 and a Th2 cytokine (IL-4 and IL-5)-mediated delay in viral clearance4950. In addition, IL-12 has been reported to contribute to early NK cell-dependent IFN-γ production and to inhibit early influenza virus replication in mice51. In this study, the evaluation of various mouse-adapted B/Florida/04/06 virus infections revealed notable patterns in cytokine expression in the P12 (MAHA/M) and P17 (MA) groups relative to the parental (P0) group (Fig. 3). We found that proinflammatory cytokine (IL-5, IL-4, and IL-6) production was significantly increased later in the P17 (MA) viral infection process (at 5 dpi) compared with the parent (P0) virus infection. In addition, the immunoregulatory interleukin IL-12p70, which is involved in the differentiation of naïve T cells into Th1 cells, was robustly induced early in the course of infection (at 1 and 3 dpi) in the P12 (MAHA/M) and P17 (MA) groups. This increased expression of inflammatory cytokines may be predominantly reflective of the enhanced virulence of these viruses in infected mice. Thus, the induction of proinflammatory cytokines was positively correlated with viral loads in the lung at the early stages of infection, and their prolonged expression may be the cause of the high mortality.

Ferrets are a suitable model in which to study influenza infection because they manifest human flu-like symptoms and immune responses1752. In fact, the utility of the ferret model for studying the clinical and pathogenic features of influenza B virus infection have been demonstrated previously. Stephen S. H. et al. found that certain influenza B strains, such as B/Brisbane/60/2008, establish persistent lower respiratory tract infections in ferrets, suggesting that the influenza B virus may cause moderate to severe clinical illness in these animals53. In this study, we demonstrate that mouse-adapted B/Florida/04/06 viruses not only replicate in the upper and lower respiratory tracts of ferrets but they were also efficiently transmitted from infected to naïve animals under both contact and aerosol droplet transmission conditions (Fig. 4). Interestingly, only the P17 (MA) virus transmitted readily to direct contact- and aerosol contact-naïve ferrets at 3 dpc, while no virus was recovered from either contact group exposed to the P0 virus-infected ferrets. To our knowledge, this is the first report of ferret-to-ferret aerosol droplet transmission of influenza B virus.

In addition, the results with three intermediate passage mutants (P5, P9, and P12) suggest that the limited replication of the influenza B virus in the upper respiratory tract could explain why influenza B virus transmission has not been previously documented in ferrets. Nevertheless, we hypothesize that the P12 virus (with the HAN424G, M1N240V, and BM2N14S mutations) possesses an increased ability to replicate in ferrets and that the mutations present in this virus are responsible for its ability to transmit to ferrets via direct contact. Furthermore, additional NP and PA gene mutations, which are closely associated with polymerase activity, are crucial for efficient viral replication in the upper respiratory tract and for aerosol transmission in the ferret model, as observed in P17 (MA) virus-infected ferrets. Interestingly, the NP and PA genes are core components of the viral ribonucleoprotein (vRNP) complex of influenza viruses and are also known to harbor host range determinants5455. In particular, the PB2 and PA polymerase subunits of influenza A viruses are essential for the efficient replication of avian influenza viruses in mammals145657 and also affect the host range and pathogenicity of avian influenza viruses58. A recent study has found that an arginine to lysine substitution at position 185 of an H5N1 virus PA protein significantly affects virulence and pathogenicity in mice59. Interestingly, an NP protein mutation similar to that observed in the B/Florida/ma04/06 virus (K294N) has been identified in the avian A/H7N1 NP protein (V284M) after passage in ferrets, resulting in the acquisition of the ability to transmit via aerosols. This mutation is in the C-terminal PB2 interaction domain and occurs in parallel with a mutation in the PB2 N-terminal domain (PB2T81I)60, indicating that together these mutations may enhance polymerase activity. Further, a polymerase activity assay indicated that a point mutation at the NP 294N residue is sufficient to enhance polymerase activity in mammalian cells at both temperatures tested (33 °C and 37 °C), while a PA mutation (609I) is additionally required for high polymerase activity at a temperature representative of the upper respiratory tract, such as that in the nasal turbinate. Therefore, we hypothesize that the increased polymerase activity of the P17 (MA) mutant, which contains both the PA and NP mutations, led to the high viral titer observed in the nasal turbinate and that this may have contributed to the elevated rate of transmission to ferrets exposed via direct and aerosol contact in our study.

In conclusion, this study demonstrates that mouse adaptation is a useful strategy for efficiently generating influenza B viruses for the infection of various mammalian hosts. In addition, multigenic virulence factors, including mutations in HA, M, and the vRNP complex, are required for efficient replication in the upper respiratory tract, which contributes to transmission of type B influenza viruses in ferrets. Thus, these results provide a foundation for the investigation of influenza B virus transmission in an animal model, and also emphasize the persistent surveillance for containment purposes and reduce opportunities for further genetic evolution. Moreover, this model may also prove useful for investigating influenza B virus biology, including virulence markers, as well as for vaccine development, vaccine efficacy, and antiviral compound testing in in vivo experiments.

Materials and Methods

Cells

Madin-Darby canine kidney (MDCK) cells obtained from the American Type Culture Collection (ATCC) were maintained in Eagle’s minimal essential medium (EMEM; Lonza, Inc., Allendale, NJ) supplemented with 10% fetal bovine serum (Lonza, Inc.), 1% penicillin/streptomycin (Gibco-Invitrogen, Inc., Carlsbad, CA), and 1% nonessential amino acids (Gibco-Invitrogen, Inc.). 293T human embryonic kidney cells were grown in Dulbecco’s modified Eagle’s medium containing 10% FBS and 1% penicillin/streptomycin (Gibco-Invitrogen, Inc., Carlsbad, CA).

Viruses

Influenza virus B/Florida/04/2006 (Yamagata lineage) was obtained from the Green Cross Corporation and passaged in mice. It was serially passaged 17 times in mice until a mouse-adapted, highly virulent strain was generated, which was termed the mouse-adapted (P17 (MA)) variant. The parent virus (P0), 3 intermediate passages (P5, P9, and P12), and the 17th passage (P17 (MA)) were isolated by plaque purification, as described below, and were selected based on sequence differences relative to the parental wild-type virus.

Plaque purification and selection

To isolate single-phenotype viruses that cause mortality at rates similar to those observed with the mouse-adapted virus, we plaque-purified lung isolates of the virulent mouse-adapted strain from MDCK cells, as described previously28. Briefly, supernatants of lung tissue homogenates were serially diluted 10-fold in appropriate media. MDCK cells were infected with the diluted samples in six-well plates. After 1 h of incubation, the cells were washed with PBS and overlaid with a 0.7%-agarose-medium mixture containing 1 μg/mL of L-1-tosylamide-2-phenylmethyl chloromethyl ketone (TPCK)-treated trypsin. Seventy-two hours later, eight single-plaque colonies were selected from each plate, resuspended in medium, and used to infect MDCK cells. After 48 h of incubation, viruses were harvested, and the 50% tissue culture infectious dose (TCID50) of each strain was calculated by the method of Reed and Muench29. After the colonies were selected from plates, the MDCK cells on the plates were fixed with 4% ice-cold formalin and stained with 1% crystal violet in order to count and measure plaques. The average plaque size was calculated as an arithmetic mean of the sizes of fifteen or more individual plaques. Viruses that formed plaques of ≤1.0 mm, 1.0–1.9 mm or ≥2.0 mm in diameter on average were designated as having small-size, medium-size or large-size plaque phenotypes, respectively30.

Genomic sequencing and phylogenetic analyses

Viral RNA was extracted from cell culture isolates using an RNeasy Mini Kit (Qiagen, Valencia, CA). RT-PCR and full-genome sequencing were performed under standard conditions using influenza-specific primers as previously described3132. DNA sequences were compiled and edited using the DNAStar Lasergene sequence analysis software package version 5.0 (Madison, WI).

In vitro replication assays

MDCK cells were prepared in 6-well plates and allowed to grow for 24 h, after which they were infected with B/Florida/04/2006 (P0), three intermediate-passage viruses (P5, P9, and P12), and one 17th passage (P17 (MA))virus at a multiplicity of infection of 0.001. Infected cells were incubated at 33 °C in appropriate medium containing 0.2% bovine serum albumin and TPCK-treated trypsin. Supernatants were collected at 12, 24, 48, and 72 h post infection (p.i.), and the viruses were titrated (log10 TCID50) in MDCK cells by the method of Reed and Muench29.

Studies in mice

Five groups of BALB/c (H-2d) mice (5-week-old females weighing ≥18 g/mouse; Samtaco, Seoul, Republic of Korea) were anesthetized with an intraperitoneal injection of a Zoletil/xylazine mixture (Zoletil 50®, 80 mg/kg, Virbac, France; Rompun®, 20 mg/kg, Bayer HealthCare, Germany) and inoculated intranasally (i.n.) with 6.0 log10 TCID50 of virus. The mice (n = 10/group) were monitored for weight loss and mortality daily for fourteen days. A mouse was euthanized and, for the purpose of data analysis, was considered to have died if its body weight was reduced by >25% compared with its original body weight. The experimental animals were used and handled according to the principles described in the Guide for the Health and Conditions of Mice33. Each group of mice (n = 5) were euthanized on days 1, 3, 5, and 7 p.i. Lungs were collected and homogenized (1 g/mL) in cold phosphate-buffered saline (PBS) containing antibiotics (0.1% penicillin/streptomycin; Gibco). The supernatants were serially diluted 10-fold and inoculated into MDCK cells for virus titration (log10 TCID50 per gram of tissue). After 72 h, hemagglutination assays were performed using 0.5% turkey red blood cells. To determine the 50% mouse lethal dose (MLD50) of the viruses, we inoculated groups of ten mice i.n. with 10-fold serial dilutions containing 101 to 106 TCID50 of virus in a 30 μL volume. The MLD50 was expressed in terms of log10 TCID50. All TCID50 and MLD50 calculations were performed according to Reed and Muench29.

Studies in ferrets

Sixteen- to 18-week-old adult outbred female ferrets (Mustela putorius furo) weighing 0.5 to 0.8 kg (Wuxi Sangosho Pet Park Co., Wuxi, China) were tested for the absence of antibodies to currently circulating influenza viruses (H5, H7, H9, pH1N1, human seasonal H1N1 and H3N2, and all of the viruses used in this study) using hemagglutination inhibition (HI) assays. For pathogenesis and transmission experiments, 106.0 TCID50 of the mouse-adapted influenza B viruses (P0, P5, P9, P12, and P17) in 1 mL of sterile PBS was instilled intranasally (i.n.) (500 μL/each nostril) into five groups of ferrets (n = 12/group) under anesthesia. After 24 h, inoculated animals in each group were individually paired and cohoused with direct contact ferrets (n = 2). For aerosol contact, naïve ferrets (n = 2) were housed in the other half of a cage containing inoculated ferrets, and the ferrets were separated by two stainless steel grids, allowing for airborne virus transmission but preventing transmission via direct contact. The ferrets were weighed and observed for clinical signs of infection daily for up to 11 days post infection (dpi). Nasal washes were collected every other day from the inoculated ferrets and daily from the contact ferrets for virus titration. Ferrets (n = 2) from each inoculated group were euthanized at 1, 3, 5, and 7 dpi, and the viral titers in the lungs, trachea, and nasal turbinate were assessed using the TCID50 method. All remaining inoculated and contact animals were euthanized at 18 dpi (17 days post contact, dpc), and their blood samples were tested for specific antibodies to a homologous virus with the HI assay using 0.5% turkey erythrocytes.

BALF samples and cytokine detection by multiplex microbead immunoassay

Bronchoalveolar lavage fluid (BALF) samples were collected from mouse lungs at 1, 3, and 5 dpi or from uninfected control mouse lungs (five mice/group). The BALF samples were centrifuged at 12 000 rpm for 5 min at 4 °C, aliquoted, and stored at −70 °C until analysis. A multiplex biometric immunoassay using fluorescently dyed microspheres conjugated with monoclonal antibodies specific for target proteins was used for cytokine measurements according to the manufacturer’s instructions (Panomics Mouse Cytokine Assay; Affymetrix Inc., Fremont, CA). Levels of the following cytokines were measured: IL-4, IL-5, IL-6, IL-12p70, interferon-gamma (INF-γ), and tumor necrosis factor-alpha (TNF-α). Briefly, 20 μL of BALF was diluted 1:4 and incubated with antibody-coupled beads. Complexes were washed, incubated with a biotinylated detection antibody and labeled with streptavidin–phycoerythrin before cytokine concentrations were measured. Concentrated mouse recombinant cytokines were obtained from Panomics Laboratories. A range of 9350–25 000 pg/mL of recombinant cytokines was used to establish standard curves and to maximize the sensitivity and dynamic range of the assay. The cytokine levels were determined using a Luminex™ Instrumentation Systems multiplex array reader (Bio-Plex Workstation from Bio-Rad Laboratories). Cytokine concentrations were calculated using software provided by the manufacturer (Bio-Plex Manager Software).

Luciferase mini-genome reporter assays

A luciferase reporter plasmid (pHW72-B/Luc) driven by a pol I transcription unit containing the non-coding region of the influenza B virus M gene was constructed. Luciferase activity was measured according to Salomon et al.34 to compare the activities of the recombinant viral polymerase complexes. 293T cells were prepared in 24-well plates 24 h before use and were transfected with 0.1 μg each of pHW72-B/Luc, pHW-B/2000-PB2, pHW-B/2000-PB1, pHW-B/2000-PA, pHW-B/2000-NP, and pCMV-β-galactosidase plasmids using the TransIT-LT1 transfection reagent as directed. After 4 h, the transfection medium was replaced with fresh DMEM (Gibco) with 5% FBS. After 24 h, the cells were washed with PBS and lysed for 30 min with 100 μL of lysis buffer (Promega). The cell lysates were then harvested, and luciferase activity was assayed in triplicate using a Promega Luciferase Assay System. The results were normalized to the β-galactosidase activity level of the cells.

Statistical analyses

To assess significant differences in viral titers and cytokine expression levels, one-tailed unpaired Student’s t-tests were used. Welch’s correction was applied when variances were found to be inhomogeneous. All statistical analyses were performed using Prism 5 (GraphPad Software).

Ethics statement

All animals were housed at an appropriate temperature under a standard light/dark cycle. The research protocols involving the use of all animals in this study were conducted in strict accordance with the animal handling policies mandated by the Guidelines for Animal Use and Care of the Korea Center for Disease Control (K-CDC) and were approved by the Medical Research Institute of Chungbuk National University in Korea (approval number CBNU-IRB-2012-GM01 and CBNUA-856-15-01). All of the experiments were progressed in an Enhanced Animal BioSafety Level 3 (ABSL3+) facility in Chungbuk National University permitted by the Center for Disease Control in the Republic of Korea.

Additional Information

How to cite this article: Kim, E.-H. et al. Mouse adaptation of influenza B virus increases replication in the upper respiratory tract and results in droplet transmissibility in ferrets. Sci. Rep. 5, 15940; doi: 10.1038/srep15940 (2015).