Abstract
Influenza B virus is a major causative agent of respiratory disease in humans. Our study of an outbreak of influenza B virus in Wenzhou from 2011 to 2014 revealed that 163 (5.58 %) of 2921 samples were influenza B positive. Sequencing of the hemagglutinin and neuraminidase genes showed substitutions at the amino acid level. Phylogenetic analysis revealed co-circulation of the B/Victoria and B/Yamagata lineages in the Wenzhou area from 2011 to 2014. Multiple viral introductions from both Chinese and international sources played important roles in endemic co-circulation and transmission in coastal southeastern China.
Introduction
Influenza B virus has been a major pathogen in seasonal influenza outbreaks and has caused respiratory infections in humans globally. It was first isolated in 1940 during an epidemic in the USA (B/Lee/40). Influenza B virus is a member of the family Orthomyxoviridae and is closely related to influenza A viruses, which are similar in viral structure, genome organization and epidemiology [1–4]. Influenza B virus differs from influenza A virus, which has a diversity of subtypes according to surface glycoproteins, in having no subtypes, but it has been separated into two main antigenically distinct lineages, Victoria (B/Victoria/2/87-like) and Yamagata (B/Yamagata/16/88-like), since 1983, based on an analysis of the hemagglutinin gene [5]. Many studies have reported both types to have been predominant during different periods and in different geographic regions worldwide [2, 6, 7].
Wenzhou, a city in southeastern Zhejiang Province, China, includes four districts and 10 counties and is one of the important economic and business centers in Zhejiang. Infectious diseases such as pandemic H1N1 and foot-and-mouth disease have been monitored in Wenzhou, and several outbreaks of these pathogen-caused illnesses were dealt with during the last decade according to surveillance systems established by public health departments in China. Influenza B has now become one of the major public-health problems, as there have been many sporadic cases in recent years. Mutations in both the hemagglutinin (HA) and neuraminidase (NA) genes have allowed influenza B virus to circumvent the immune response in humans, to persist in human populations, to circulate in an endemic environment, and to cause recurrent seasonal epidemics [8–11]. Therefore, by combining the results of molecular and phylogenetic data, we attempted to determine (1) the molecular characteristics of both the hemagglutinin and neuraminidase genes and (2) the phylogenetic pattern of the influenza B virus in the Wenzhou area.
Material and methods
This study was approved by the ethics committee of the Zhejiang Provincial Center for Disease Control and Prevention (ZJCDC), China. Following the ‘Surveillance Program of Influenza in China’, published by the National Health and Family Planning Commission (NHFPC), throat swabs and/or nasopharyngeal samples were collected in local hospitals and delivered to the ZJCDC from 2011 to 2014. In total, 2921 samples were obtained from patients exhibiting flu-like symptoms.
Viral RNA was extracted using an RNeasy Mini Kit (Roche) according to the manufacturer’s instructions. Influenza B virus infection was identified and genotyped by multiplex real-time PCR reactions using an AgPath-IDTM One-Step RT-PCR Kit (Life Technologies) following the protocol for the surveillance program.
Positive specimens were cultured in Madin-Darby canine kidney (MDCK) cells, a gift from the National CDC, for 5 to 7 days. Specific-pathogen-free embryonated chicken eggs were also used for virus isolation. Six 9- to 11-day-old chicken embryos were each inoculated with 300 µl of sample by the chorioallantoic sac route. The eggs were incubated for 48 hours at 35 °C. Cultured supernatants and allantoic fluids were tested by hemagglutination inhibition (HAI). Samples testing negative for hemagglutination were processed a second time.
Positive samples were subjected to RT-PCR amplification and sequencing of the hemagglutinin and neuraminidase genes. RT-PCR reactions for both the hemagglutinin (HA) and neuraminidase (NA) genes were done according to the surveillance program of Takara’s kit (Table S1). Sequencing was performed using an ABI 3730xl DNA Analyzer. All virus sequences have been deposited in the Global Initiative on Sharing All Influenza Data (GISAID) database (EPI630146-EPI630185).
Both the HA and NA gene were assembled and aligned along with additional sequences downloaded from GenBank. Variant positions in the nucleotide and amino acid sequences were checked using Geneious 4.8.5 (http://www.geneious.com). Identical indexes for both HA and NA were calculated using DNAStar Lasergene v7.1 (http://www.dnastar.com). Dataset-specific models that were selected using the Akaike information criterion in Modeltest 3.7 were analyzed [12]. Maximum-likelihood (ML) analysis was done using RAxML v7.2.8 (http://sco.h-its.org/exelixis/software.html). The optimal ML tree and bootstrap percentages (BP) were estimated in the same run. The ML BP values were obtained from 1000 bootstrap replicates using the rapid bootstrap algorithm. BEAST 1.6 (http://beast.bio.ed.ac.uk/Main_Page) was employed to date the divergence of the tree branches in each phylogenetic tree (HA and NA) based on an uncorrelated exponential distributed relaxed-clock model for our sample [13]. A burn-in of 10 % was used, and the convergence of all parameters was assessed using the software TRACER within the BEAST package.
Results and discussion
We collected and tested 2921 samples, including throat swabs and nasopharyngeal swabs from across the Wenzhou area of Zhejiang Province from 2011 to 2014. One hundred sixty-three were positive for influenza B virus, a rate of 5.58 % (163/2921) for all samples (Table S2). This rate was higher than those observed in Thailand and the United States during a similar period of time [7, 14]. Most influenza B virus infections occurred in 2012, at a rate of 8.85 %, while the fewest cases were in 2014 (4.09 %). The majority of samples confirmed influenza B virus infections among children under 10 years of age at a rate of 47.23 % (77/163), which is in line with previous reports [15].
We detected many amino acid substitutions in both the HA and NA segment matrices, since variations in these sequences play a major role in pathogenicity of viral strains (Table 1). Sequence variations in the HA protein were divided into Yamagata and Victoria groups. Twenty-seven substitutions in the HA segment in the Yamagata lineage were detected, while in the Victoria lineage, there were 15 substitutions. A previous study reported four major antigenic epitopes – 120-loop (116-137), 150-loop (141-150), 160-loop (162-167) and 190-helix (194-202) – on the membrane distal domain of the HA1 region of the influenza B virus [11]. We observed many mutations in these important areas: for the Yamagata lineage, we found P123A, N131K and H137Q in the 120-loop area, N141D and R144K in the 150-loop area, R164K and N165S/I in the 160-loop area, and both T196A and K197E in the 190-helix area. We detected only T144N in the 150-loop in the Victoria lineage. A previous study indicated that the 120-loop (116-137) epitope dictates the antigenicity of HA [16]. We found that most mutations were in the 120-loop, especially in the Yamagata lineage, which was in line with several previous studies [11, 17]. N131Y detected in this study was considered a potential alteration in antigenicity of the Wenzhou strains. Variants in the 150-loop, 160-loop and 190-helix also obtained here, may contribute to differences in avoiding antigen recognition by neutralizing antibodies and to maintaining the structural integrity of some key proteins, such as receptor-binding proteins [7, 18]. Forty-two mutations were detected in the NA segment matrix constituted by the Wenzhou sequences when compared to both B/Yamagata/16/1988 and B/Victoria/2/1987. Substitutions found at positions 116, 117 and 402 were considered to be related to sensitivity to NA inhibitors and the function of drug resistance [4, 7]. No biological function has been suggested for many other amino acid changes identified in the NA segment. Pairwise identities were calculated according to different matrices in different segments. In the Victoria lineage, 97.5 % pairwise identity was determined, while 97.7 % was determined in the Yamagata lineage. Pairwise identity was estimated as 93.4 % in the total HA segment and 95.8 % in the NA segment.
Phylogenetic analysis of HA showed that influenza B viruses in the Wenzhou area from 2011 to 2014 separated into two main lineages, Yamagata and Victoria, according to the hemagglutinin gene. In the Yamagata lineage, most sequences, including the Wenzhou-area samples, were clustered into one monophyletic clade with the most common recent ancient time in the year 2004 (Fig. 1 and Table 2). Within this clade, sequences from the Wenzhou area formed two monophyletic clades: Y1 and Y2. Strains from Lishui (Lishui/297/2011) and Hangzhou (Hangzhou/19/2012) showed high homology to Y1 and Y2, respectively. The phylogenetic pattern for the Wenzhou sequences from 2011 to 2014 in the Victoria lineage showed different topological structures, as strains from Wenzhou scattered into different clades in the Victoria part. Strains mainly from 2011 to 2012 formed a monophyletic clade with sister strains from other locations such as Taizhou, Nanjing and Brisbane. The NA phylogenetic tree showed that sequences from the Wenzhou area from 2011 to 2014 clustered into one large monophyletic clade with other NA sequences from outside Wenzhou, such as from Shanghai, Beijing, Chongqing, Hubei, Taiwan, Fujian, Guangdong and Brisbane, with an estimated emergence time of 1994.
Maximum-likelihood tree for the HA and NA segments of influenza B virus. Bootstrap values based on 1000 replicates are indicated above the branches. Blue numbers below the branches show the time of the most recent common ancestor for each node. Blue bars on nodes are 95 % highest posterior density (HPD) intervals. A, HA segment; B, NA segment (color figure online)
Multiple introductions of influenza viruses, such as A/H1N1 and A/H3N2, were reported previously [2, 19–21]. In this study, phylogenetic trees based on both the HA and NA segments and the topological structure of the Yamagata and Victoria lineages were in accordance with the conclusion that the dominant circulating influenza B viruses in the Wenzhou area from 2011 to 2014 were the result of multiple introductions from multiple locations outside of the Wenzhou area in different years or even within a single year. The HA segment in the Yamagata lineage indicated that one dominant circulating influenza B virus in the Wenzhou area during 2013 was introduced from Hangzhou, the capital of Zhejiang Province, while the other monophyletic clade in the same Yamagata lineage showed that another independent introduction in 2012 was from the Lishui area located in southwestern Zhejiang (Fig. 1). These two clades also have different times for the most recent common ancestor as estimated by the MCMC algorithm (Table 2). A similar pattern was observed for the Victoria lineage. Two introduction events related to the Wenzhou strains were from outside of Zhejiang: one was from Nanjing, Jiangsu Province, and the other was from Brisbane, Australia. Additionally, our phylogenetic analysis of the NA segment showed evidence of local persistence of several phylogenetically distinct monophyletic clades circulating in the Wenzhou area from 2011 to 2014. Viruses belonging to those clades may have been introduced from different places. For example, the strain from Brisbane shared high phylogenetic similarity to strains in Wenzhou during 2011 and 2012. There are several possible reasons for this phenomenon. Wenzhou hosts a large number of migrants, as shown in the latest government census report in 2011. They account for one third of the population in the Wenzhou area. Migrants who follow the Chinese custom of returning to their home villages during important festivals throughout the year may facilitate the transmission of viruses from one location to another. This was shown to be the case in previous studies of other viruses [2, 22–24]. Additionally, Wenzhou has connections with other regions in China and, as the third largest city in Zhejiang Province, even internationally through foreign trade, and has a higher risk of other diseases such as chikungunya. Wenzhou, and, in fact, all of Zhejiang, is a popular travel destination for both business travelers and tourists, thereby providing new channels for intercontinental introductions and transmissions. Long-term influenza surveillance is therefore essential for early detection and for providing early warning to the public of the potential for viral transmission and circulation.
This study had several limitations. First, we were unable to investigate the full epidemiological information for the samples we collected. Although we observed that children and young teenagers were the main susceptible population compared to other age groups, which was similar to findings in other studies of influenza B virus, we were unable to identify the proportion of Victoria and Yamagata in the population for each year in the Wenzhou area, which may have reflected background immunity within the population. Second, we failed to obtain the whole genome sequence for all positive samples of influenza B virus. Reassortment analysis was therefore limited. As a result, a blind area exists in the evaluation of the evolutionary processes and endemic speciation of influenza B virus in the Wenzhou area. More samples and detailed epidemiology information are needed to understand the circulation, dynamics and evolution of influenza B virus in this region.
To conclude, this study is the first to highlight the basic epidemiological and molecular characteristics of influenza B virus strains, especially the hemagglutinin and neuraminidase genes, for the 2011 to 2014 outbreak in the Wenzhou area. Amino acid substitutions were identified, and phylogenetic relationships between other strains from different locations were evaluated. Multiple introductions from outside the Wenzhou area were revealed, and potential circulation was detected. Our study revealed the variation, co-circulation and prevalence of both influenza B virus lineages during the period 2011 to 2014 in the coastal zone of southeastern China. Further studies are needed to determine the detailed dynamics and evolution of influenza B virus in this region.
References
Chen R, Holmes EC (2008) The evolutionary dynamics of human influenza B virus. J Mol Evol 66:655–663
Tan Y, Guan WD, Lam TTY, Pan SH, Wu SG, Zhan YQ, Viboud C, Holmes EC, Yang ZF (2013) Differing epidemiological dynamics of influenza B virus lineages in Guangzhou, Southern China, 2009–2010. J Virol 87:12447–12456
Zhu N, Li P, Yu JF, Li YD, Zhao JR, Xia H, Tang S, Zhang Z, Kou Z, Yin SY, Fan ZJ, Li TX (2013) Molecular characterization of influenza B viruses isolated in east-central China in 2009–2010. Virus Genes 46:28–38
Ali G, Amer HM, Ahnajhdi FN (2014) Hemagglutini and neuraminidase genes of influenza B viruses circulating in Riyadh, Saudi Arabia during 2010–2011: evolution and sequence analysis. J Med Virol 86:1003–1016
Kaneage Y, Sugita S, Endo A, Ishida M, Senya S, Osako K, Nerome K, Oya A (1990) Evolutionary pattern of the hemagglutinin gene of influenza B viruses isolated in Japan: cocirculating lineages in the same epidemic season. J Virol 64:2860–2865
Rota PA, Wallis TR, Harmon MW, Rota JS, Kendal AP, Nerome K (1990) Cocirculation of two distinct evolutionary lineages of influenza type B virus since 1983. Virology 175:59–68
Tewawong N, Suwannakarn K, Prachayangprecha S, Korkong S, Vichiwattana P, Vongpunsawad S, Poovorawan Y (2015) Molecular epidemiology and phylogenetic analyses of influenza B virus in Thailand during 2010 to 2014. PLoS One 10(1):e0116302
Burnet SFM (1955) Principles of animal virology. Academic Press, New York
Hay AJ, Gregory V, Douglas AR, Lin YP (2000) The evolution of human influenza viruses. Philos Trans R Soc Lond B Biol Sci 356:1861–1870
Bedford T, Suchard MA, Lerney P, Dudas G, Gregory V, Hay AJ, McCauley JW, Russell CA, Smith DJ, Rambaut A (2014) Integrating influenza antigenic dynamics with molecular evolution. eLife 3:e01914
Fang QX, Gao Y, Chen MF, Guo XL, Yang X, Wei L (2015) Molecular epidemiology and evolution of influenza A and B viruses during winter 2013–2014 in Beijing, China. Arch Virol 160:1083–1095
Posada D, Crandall KA (1998) Modeltest: testing the model of DNA substitution. Bioinformatics 14:817–818
Drummond AJ, Suchard MA, Xie D, Rambaut A (2012) Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol 29:1969–1973
WHO (2014) Review of the 2013–2014 winter influenza season, northern hemisphere. WER 89:245–256
Zaman RU, Alamgir ASM, Rahman M, Azziz-Baumgartner E, Gurley ES, Sharker MAY, Brooks WA, Azim T, Fry AM, Lindstrom S, Gubareva LV, Xu X, Garten RJ, Hossain MJ, Khan SU, Faruque LI, Ameer SS, Klimov AI, Rahman M, Luby SP (2009) Influenza in outpatient ILI case-patients in national hospital-based surveillance, Bangladesh, 2007–2008. PLoS One 29(4):e8452
Wang Q, Cheng F, Lu M, Tian X, Ma J (2008) Crystal structure of unliganded influenza B virus hemagglutinin. J Virol 82:3011–3020
Verhoeyen M, Van Rompuy L, Jou WM, Huylebroeck D, Fiers W (1983) Complete nucleotide sequence of the influenza B/Singapore/222/79 virus hemagglutinin gene and comparison with the B/Lee/40 hemagglutinin. Nucleic Acids Res 11:4703–4712
Ni F, Kondrashkina E, Wang Q (2013) Structural basis of the divergent evolution of influenza B virus hemagglutinin. Virology 446:113–122
Bedford T, Nelson MJ, Chan KH, Chen R, Vijaykrishna D, Halpin RA, Stockwell TB, Lin X, Wentworth DE, Ghedin E, Guan Y, Peiris JS, Riley S, Rambaut A, Holmes EC, Smith GJ (2011) Temporally structured metapopulation dynamics and persistence of influenza A H3N2 virus in human. Proc. Natl Acad Sci USA 108:19359–19364
Cheng X, Tan Y, He M, Lam TT, Lu X, Viboud C, He J, Zhang S, Lu J, Wu C, Fang S, Wang X, Xie X, Ma H, Nelson MJ, Kung HF, Holmes EC (2013) Epidemiological dynamics and phylogeography of influenza in southern China. J Infect Dis 207:106–114
Nelson MI, Edelman L, Spiro DJ, Boyne AR, Bera J, Halpin R, Sengamalay N, Ghedin E, Miller MA, Simonsen L, Viboud C, Holmes EC (2008) Molecular epidemiology of A/H3N2 and A/H1N1 influenza virus during a single epidemic season in the United States. PLoS Pathog 4:e1000133
Sun Y, Yan JY, Mao HY, Zhang L, Lyu QF, Wu ZH, Zheng W, Feng C, Zhang YJ (2013) Characterization of the complete genome of Chikungunya in Zhejiang, China, using a modified virus discovery method based on cDNA-AFLP. PLoS One 8:e83914
Zhang YJ, Mao HY, Yan JY, Zhang L, Sun Y, Wang XY, Chen Y, Lu YY, Chen EF, Lv HK, Gong LM, Li Z, Gao J, Xu CP, Feng Y, Ge Q, Xu BX, Xu F, Yang ZN, Zhao GQ, Han JK, Guus K, Li H, Shu YL, Chen ZP, Xia SC (2014) Isolation and characterization of H7N9 avian influenza A virus from humans with respiratory disease in Zhejiang, China. Virus Res 189:158–164
Mao HY, Guo B, Wang FY, Sun Y, Lou XY, Chen Y, Zhang L, Wang XY, Li Z, Liu SL, Qin SW, Wei JC, Pang ZF, Chen ZP, Zhang YJ (2015) A study of family clustering in two young girls with novel avian influenza A (H7N9) in Dongyang, Zhejiang Province, in 2014. J Clin Virol 63:18–24
Acknowledgments
Dr. David E. Boufford is gratefully acknowledged for editing the manuscript.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Funding
This work was supported by the Zhejiang Provincial Program for the Cultivation of High-level Innovative Health Talents and Monitor Technology Platform of Infectious Diseases of the State Major Science and Technology Special Projects During the 12th Five-Year plan of China (2012ZX10004-210).
Conflict of interest
None.
Ethics standard
This study was approved by the ethics committee of the Zhejiang Provincial Center for Disease Control and Prevention (ZJCDC), China. This study was done in accordance with the Declaration of Helsinki, Good Clinical Practice, and Chinese regulatory requirements.
Additional information
D. Chen and X. Wen contributed equally.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Chen, D., Wen, X., Sun, Y. et al. Multiple viral introductions: molecular characterization of influenza B virus in Wenzhou, Zhejiang, China, from 2011 to 2014 based on hemagglutinin and neuraminidase genes. Arch Virol 161, 1005–1013 (2016). https://doi.org/10.1007/s00705-015-2721-7
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00705-015-2721-7