key: cord-0006546-30ulsqqu authors: Sha, Jianping; Chen, Xiaowen; Ren, Yajin; Chen, Haijun; Wu, Zuqun; Ying, Dong; Zhang, Zhiruo; Liu, Shelan title: Differences in the epidemiology and virology of mild, severe and fatal human infections with avian influenza A (H7N9) virus date: 2016-02-18 journal: Arch Virol DOI: 10.1007/s00705-016-2781-3 sha: 170eb018c43083eef1333bc6ff6ffd95d19a9a96 doc_id: 6546 cord_uid: 30ulsqqu A novel avian influenza A (H7N9) virus caused 5-10 % mild and 30.5 % fatal human infections as of December 10, 2015. In order to investigate the reason for the higher rate of fatal outcome of this infection, this study compared the molecular epidemiology and virology of avian influenza A (H7N9) viruses from mild (N = 14), severe (N = 50) and fatal (N = 35) cases, as well as from non-human hosts (N = 73). The epidemiological results showed that the average age of the people in the mild, severe and fatal groups was 27.6, 52 and 62 years old, respectively (p < 0.001). Males accounted for 42.9 % (6/14), 58.0 % (29/50), and 74.3 % (26/35) of cases in the mild, severe and fatal group respectively (p = 0.094). Median days from onset to start of antiviral treatment were 2, 5 and 7 days in the mild, severe and fatal group, respectively (p = 0.002). The median time from onset to discharge/death was 12, 40 and 19 days in the mild, severe and fatal group, respectively (p < 0.001). Analysis of whole genome sequences showed that PB2 (E627K), NA (R294K) and PA (V100A) mutations were markedly associated with an increased fatality rate, while HA (N276D) and PB2 (N559T) mutations were clearly related to mild cases. There were no differences in the genotypes, adaptation to mammalian hosts, and genetic identity between the three types of infection. In conclusion, advanced age and delayed confirmation of diagnosis and antiviral intervention were risk factors for death. Furthermore, PB2 (E627K), NA (R294K) and PA (V100A) mutations might contribute to a fatal outcome in human H7N9 infection. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s00705-016-2781-3) contains supplementary material, which is available to authorized users. Since the identification of the first human infection with a novel avian influenza A (H7N9) virus in China in 2013 [6] , 669 confirmed cases have been reported worldwide up to 03 June 2015, according to the China annonmanet for the legislated Infectious diseases (http://www.nhfpc.gov.cn/jkj/ s3578/201509/a67694612ff243a6a0788ad6a949ade0.shtml). The majority of cases were identified in the south and east of China [1, 12] . Almost all infected individuals were hospitalised, and 36 % (158/444) died [12] . The fatality rate was much higher than that observed for seasonal influenza in the United States (*0.04 %) [12] and a previous H7N7 outbreak, but it was lower than that of H5N1 (70.0 %) in China [2] . Numerous studies have indicated that the higher mortality rate is related to delays in diagnosis and antiviral intervention, as well as old age and chronic diseases [17, 19] . Accumulating evidence has indicated that the infection source of H7N9 was presumably poultry and related products; however, human-to-human transmission was implicated in 20.0 % of confirmed cases [3, 10, 21, 26, 36, 38] . Nonetheless, evidence of sustained human-to-human transmission is very limited. Gene sequencing, receptor-bindingsite (RBS) binding assays, and a series of animal studies have revealed that the H7N9 virus has an increased capacity to bind to mammalian respiratory cells. For example, substitutions in the HA at amino acid residues Q226L and G186V, and mutations in PA, PB1 and PB2 (D701N and E627K) are associated with higher virus replication rates, severity of illness and transmissibility [4, 23, 39, 44] . In accordance with this, Wang et al. reported that PB2 E627K mutation slightly increased the case fatality rate [33] . Additional studies by Hai et al. and Sleeman et al. showed that that the emergence of the Arg292Lys mutation in the neuraminidase (NA) gene appears to be associated with failure of response to NA inhibitors and to adverse clinical outcomes [9, 29] . A further important study indicated that H7N9 and H9N2 viruses are co-circulating in some live-bird markets, which might alter host adaptation and virus pathogenicity [14] . This raises concerns regarding the potential for H7N9 to adapt and lead to an infection pandemic. The present study was an epidemiological investigation of 99 avian influenza A (H7N9) cases that occurred from 2013 to 2015. We compared the epidemiological features and signature amino acids, evolution and reassortments of the H7N9 virus with full genome sequences from mild, severe and fatal cases. We attempted to find a potential association between H7N9 mutations and reassortments and clinical outcome. The results might help in improving understanding of the epidemiology of this H7N9 strain, thereby helping to decrease the disease burden and to address the H7N9 challenge. The present study was approved by the Medical Ethical Committee and the National Health and Family Planning Commission. Written informed consent was obtained from all subjects who participated in the study and/or their families. The activities regarding the collection of human samples and the study protocol were approved by ethics committee of the local hospital and local institutional review board located in Zhejiang, Guangdong Province, and Shanghai, China. (1) Definition of a confirmed H7N9 case The case definitions were established on the basis of 'The Diagnosis and Treatment Programs of Human Infections with H7N9 Virus', issued by the National Health and Family Planning Commission of the People's Republic of China [34] . A confirmed H7N9 case was defined as a case in which a patient had influenza-like illness or in which respiratory specimens had tested positive for H7N9 virus via either of the following: isolation of H7N9 virus or positive results by rRT-PCR assay for H7N9, or a fourfold or greater increase in antibody titer for H7N9 virus based on testing of an acute serum specimen (collected 7 days or less after symptom onset) and a convalescent serum specimen collected at least 2 weeks later [27] . (2) Definition of a mild H7N9 case A case in which an individual with confirmed H7N9 virus infection met the respiratory infection criteria was classified as a mild case, presenting with mild respiratory symptoms that did not have any complications throughout the clinical course, such as acute respiratory distress syndrome (ARDS), multi-organ failure, hypoxaemia, etc. [5] . (3) Definition of a severe H7N9 case A case in which an individual with confirmed H7N9 virus infection met any one of the following criteria was classified as a severe case: presenting with severe respiratory symptoms with any complications including ARDS, shock, multi-organ failure, hypoxemia, etc., requiring hospitalization or intensive care unit admission or mechanical ventilation for medical reasons. The objective index was as follows: 1) X-ray showed lesions in multiple lobes or disease progression [50 % within 48 hours; 2) dyspnea with a respiratory rate [ 24 breaths per minute; 3) hypoxemia with oxygen saturation B 92 % on oxygen at a flow rate of 3 to 5 litres/min; 4) shock, ARDS or multiple organ dysfunction syndrome [5, 27, 41] . (4) Definition of a fatal H7N9 case A case in which an individual with confirmed H7N9 virus infection died of this infection was classified as a fatal case. People were predominantly exposed to the virus in two ways. The first of these was via exposure to poultry: individuals who came into direct or indirect contact with, or within proximity of, healthy, sick or dead poultry (including all types of poultry or birds, e.g., chickens, ducks, geese, pet birds and pigeons), having poultry in the neighborhood and eating poultry products that were not properly processed, were all included. The second manner of exposure was human-to-human transmission. This refers to whether infected individuals had had close contact with a confirmed or probable human case (any time from the day before the onset of illness to death or during the period of hospitalization) in the 2 weeks before the onset of his/her symptoms [27] . A total of 552 confirmed cases of avian H7N9 influenza were analyzed. However, only 99 of the 552 confirmed cases were selected for further analysis, as the clinical outcomes of the remaining 453 cases were not known. The selected cases included 14 mild, 50 severe and 35 fatal infections, and whole-genome sequences were generated from all 99 cases. Of these, 84 sequences with numbers are detailed in supplementary Table 1 When a suspected case of H7N9 virus infection had been confirmed, the provincial epidemiologists and local public health doctors conducted the initial field investigations using a standard questionnaire. Field investigators interviewed the patient in confirmed case(s) and/or the relatives of the patient to determine exposure history (poultry exposure and human-to-human transmission), including the dose, duration and intensity of the exposure during the 2 weeks prior to the onset of illness. We calculated the incubation period as the number of days since the date of last exposure to the onset of illness. The questionnaire included data regarding dates, times, frequency and patterns of exposure to poultry and/or other animals, as well as bird environments. All available medical records, including the date of onset of illness, visits to clinical facilities, hospitalization, antiviral treatment, and clinical outcomes, were provided by local clinical health workers. A total of 99 respiratory specimens, including 80 throat swabs and 19 sputum specimens, from patients admitted to different hospitals (located in Zhejiang and Guangdong Province) were collected and shipped at 4°C to the local hospitals and local CDCC for H7N9 laboratory testing by rRT-PCR assay. All of the surveyors and laboratory technicians were strictly trained to ensure that the interviews and laboratory investigations were conducted according to uniform standards. Total viral RNA was extracted from the collected samples using a QIAGEN RNeasy Mini Kit according to the manufacturer's instructions. The specific primer and probe sets were provided by the Chinese CDC. The detection limits of the RT-PCR assays were approximately 100 copies/ml. Whole-genome sequencing of the isolated H7N9 samples was carried out by amplification and sequencing of the eight genomic fragments. Nucleotide sequences of the amplified products were determined directly by dideoxy sequencing, using an ABI PRISM BigDye Terminator Cycle Sequencing Kit [6] . Evolutionary relationships of the generated sequences were identified using the maximum-likelihood method, as implemented in MEGA 6.0 software (http://www.mega software.net). The current addresses of all confirmed cases were geocoded using the Google Map geocoding service (https://googledevelopers.appspot.com/maps/documentation/ javascript/examples/geocoding-simple). After obtaining the X (longitude) and Y (latitude) coordinates, we used ArcGIS (ArcGIS, version10.2). The world basemap used was the publically available map maintained by ESRI (http://www.arcgis.com/home/item.html?id=3864c63872d 84aec91933618e3815dd2). For spatial analysis, a spatial distribution map of the cases was constructed. All statistical analyses were conducted using Statistical Analysis System, version 9.2 (SAS Institute, Cary, NC, USA). Quantitative measurements are presented as the median and range of the observed values, and qualitative measurements are presented as relative and absolute frequencies. An analysis of variance (F test) was applied to the measurement data. Chi-square tests (x 2 ) were used to compare the distribution of the different variables of qualitative measurements between the three groups. Fisher's exact test was used in the analysis of contingency tables when the sample sizes were small (the expected values in any of the cells of a contingency table were below 5; the number of total samples was no more than 40; the data were very unequally distributed among the cells of the table). All p-values given are two-sided and were considered statistically significant at 0.05. Epidemiological characteristics in mild, severe and fatal cases As of 31 January 2015, H7N9 had resulted in 552 confirmed infections, with a fatality rate of 37.7 %. The epidemiological distribution of mild, severe and fatal cases is shown in Fig. 1 and virtual features are shown in Fig. 2 total of 42.9 % (6/14) of the subjects with mild disease were males and 57.1 % (8/14) were females, 58.0 % (29/ 50) of the severe cases were males and 42.0 % (21/50) were females, and 74.3 % (26/35) of the subjects who died were males and 25.7 % (9/35) were females. There was no difference in the gender distribution between the three groups (p = 0.094). However, the average age in the mild cases was 27.6 years (range, 5 months-82 years), versus 52 years (range 1.5-85 years) for the severe cases and 62 years (range 27-91 years) for cases of death. There was a significant difference in the age distribution between the three groups (p \ 0.001). The most predominant age group was 0-, 60-and 60-in the mild, severe and fatal group, respectively (Fig. 3) . The median time from exposure to onset was 4 (1-7), 2 (0-6) and 5 (1-13) days (p = 0.037) in the mild, severe and fatal group, respectively. The median time from onset to consultation was 1 (0-5), 3 (0-13) and 3 (0-7) days (p = 0.059); from onset to laboratory confirmation, it was 3 (0-8), 8 (1--26) and 11 (0-41) days (p = 0.001); from onset to beginning antiviral treatment, it was 2 (0-5), 5 (0-10) and 7 (1-12) days (p = 0.002); and from onset to discharge/ death, it was 12 (0-32), 40 , and 19 (3-85) days (p.001) in the mild, severe and fatal group, respectively (Table 1) . With the exception of the period from onset to consultation, all time relationships were statistically significantly different, as shown in Table 1 and Fig. 4 . Whole genome sequences were generated from 99 H7N9 human isolates and 73 non-human isolates and aligned against 38 reference sequences from databases. Multiple sequence alignments showed 99.1-99.9 % sequence identity between isolates from the three clinical groups and non-human isolates. From 31 March 2013 to 01 October 2015, the H7N9 outbreak occurred in three waves. The first of these lasted from 31 March 2013 to 30 September 2013, the second occurred between 01 October and 30 September 2014, and the third wave lasted from 01 October 2014 to 30 September 2015. Phylogenetic analysis based on the HA segments revealed that the sequences from the mild, severe and fatal groups all clustered in five clades and that clade A varied between the waves and the reported areas. Clade A included the Analysis of mutations in HA sequences indicated that there was no statistical difference in G186V among the three clinical groups and non-human isolates. Q226L was more often found in the severe and fatal groups than in the mild group, but there was no statistical difference (p = 0.080). However, Q226L was identified as the amino acid difference that distinguishes human isolates from non-human isolates (p.001). All of these mutations were associated with increased affinity for a-2,6-linked sialic acid receptors. In contrast, the G228S mutation was not found in any of the isolates from non-human or human cases, irrespective of whether they were mild, severe or fatal. All nonhuman isolates retained the amino acid 'G', which means they had high affinity for human receptors. Conversely, there was an obvious statistical difference between activated C antigen sites (N276D), related to antibody binding sites, among the three clinical groups (p = 0.046) and human vs. non-human isolates (p \ 0.001). However, no difference was found between the other antigen sites, including R57M, R57K, E114K, T133A and G186V, and the three clinical groups and non-human isolates. L168I was identified as statistically different among the mild, severe and fatal groups (p = 0.050) and between the nonhuman isolates and human isolates (p = 0.001), and this mutation had not been reported to date. The HA cleavage site from mild, severe and fatal human cases and nonhuman isolates possessed only a single amino acid R (arginine), indicating that it still had characteristics of a low-pathogenic avian influenza (LPAI) virus ( Table 2 ). The amino acid sequences from the three clinical groups and non-human isolates indicated that aa 69-73 were deleted in the stalk region of the NA protein, which potentially influences virus replication. Although E120V, H276Y and N296S substitutions were absent from sequences from clinical and non-clinical isolates, R294K was found in 2.1 % and 12.1 % of isolates from the severe and fatal cases, respectively. This suggested that resistance to the antiviral NA inhibitor oseltamivir was most often found in the fatal group (p =0.01). Compared to the nonhuman isolates, the R294K mutation rate was much higher in the human isolates (p = 0.01). Of the 35 NA sequences from the patients treated with oseltamivir, three (8.6 %) showed an R294K mutation, while five (14.3 %) exhibited V248I, V251I and N337T substitutions. The E627K substitution in the PB2 gene was more frequently found in the severe (68.2 %) and fatal (86.7 %) groups than in the mild group (33.3 %) (p = 0.006) and non-human isolates (18.6 %) (p \ 0.001). In contrast, N559T was the major mutation in the mild group compared to the other two groups (p = 0.048). No differences were observed among the sequences of the three clinical groups with regard to other site mutations, including L89V, K191E, M570I, Q591K/L and D701N. However, the mutation rates of K191E (p = 0.016), N559T (p = 0.011), M570I (p = 0.026) and E627K (p \ 0.001) differed between the non-human and human isolates, as shown in Table 2 . No significant difference between the three groups was observed in I368V and L598P/M in the PB1 segment. However, the non-human isolates had a higher rate of I368V than the human isolates (p = 0.045), as shown in Table 2 . Full-length PB1-F2 protein (87-90 amino acids) has been associated with increased virulence in mice, and PB1-F2 protein with 90 amino acids was detected in 72.7 %, 80.5 % and 90.3 % of the mild, severe and fatal human cases, respectively, with no statistical difference (p = 0.333). A total of 25.7 % of the chicken isolates and 36.7 % of the environmental virus isolates had PB1-F2 protein with 25 amino acid mutations, which was higher than that observed in the human isolates (p = 0.03). In contrast, 68.6 % of the chicken isolates and 63.3 % of the environmental virus isolates had the 90-aa variant of the PB1-F2 protein. The mutation rates of non-human isolates were much lower than those of human isolates (p = 0.05), as shown in Table 2 . Single letters refer to the amino acid (aa) found in the indicated protein at a specific site. The numbering starts with the first methionine codon for these proteins p1 value: comparison of the mutation frequency among the isolates from mild, severe and fatal cases p2 value: comparison of the mutation frequency between human and non-human isolates Qualitative measurements are presented as relative and absolute frequencies. Chi-square (v 2 ) tests were applied to compare the distribution of the different variables of qualitative measurements. Fisher's exact test was used in the analysis of contingency tables when the sample sizes were small ''-'' denotes no difference between the groups Human infections with avian influenza A (H7N9) virus 1255 the other mutations, I308V, K356R, S409N and T618K, were not statistically significantly different among the three groups. The K356R mutation rate in the non-human isolates was higher than in the human isolates (p = 0.024) ( Table 2) . A total of 11 mild, 41 severe and 32 fatal isolates were analysed at the NP segment. A substitution of N321S and D375E remained conserved in the three clinical groups and the non-human isolates (Table 2) . Similarly, NS segments from 11 mild, 42 severe and 32 fatal isolates were analysed. No differences in P42S and T92A were observed in the three clinical groups and nonhuman isolates. In addition, all H7N9 viruses had a truncated NS gene with a PDZ motif deletion. PDZ deletion might influence the pathogenicity of H7N9 viruses in human or avian hosts ( Table 2 ). It was found that the M1 protein from all of the human and non-human isolates had virulence sites (N30D and T215A) and that the M2 protein from all of the human and non-human isolates had an S31N substitution, indicating resistance to adamantine (amantadine and rimantadine), as shown in Table 2 . Since the spring of 2013, a novel avian-origin influenza A (H7N9) virus has emerged and spread among humans in China, resulting in a high fatality rate [6, 15] . The risk factors that might contribute to a fatal outcome regarding human H7N9 infection include old age, a history of smoking, chronic lung disease, immunosuppressive disorders, chronic drug use, and delayed oseltamivir intervention [19] . The present study described the major differences in epidemiological and viral features between mild, severe and fatal cases during the three epidemic waves of the worldwide H7N9 outbreak. The age and gender characteristics of H7N9 cases are unusual compared to those of H5N1 cases. In H5N1, an average age of 26 years (interquartile range: [19] [20] [21] [22] [23] [24] [25] and a more equal gender distribution were the prominent features [12] . However, the present study indicated that the average age in the mild group (27.6 years) was far younger than that in the severe (52 years) and fatal groups (62 years) [2, 32] , with 50 % of mild cases being found in children \10 years of age. However, the gender distribution has no statistical difference among three groups. Mild cases were frequently identified in family clusters and secondary cases, permitting early detection due to monitoring of the close contacts of index cases, as well as a prompt start to antiviral treatment once symptoms began. However, the predominance of older people among the severe and fatal cases is likely to be due to their having had underlying lung conditions and impaired immune functions, as well as a greater chance and high level of exposure [13, 17] . The age and gender findings were consistent across the first, second and third waves of the epidemic. In accordance with the results reported herein, a strong association has been shown between old age and a fatal outcome in (H1 and H3) seasonal influenza ([65 years old), which is different from what has been observed in H5N1 avian influenza [40] . The median incubation period was 4, 2 and 5 days in the mild, severe and fatal group, respectively, which is comparable to what has been reported previously (3.1 days) [2, 15] . The median number of days from onset to confirmation (11) and onset to start of antiviral treatment (7 days) found in the fatal group was relatively long compared to that observed in the mild and severe groups, and was also in accordance with the findings of Li et al. [15] and Liu et al [20] , which might indicate that delayed time of case confirmation and antiviral intervention are closely associated with a fatal outcome. In particular, the median number of days from the onset to outcome in fatal cases (19 days) was lower than that of onset to outcome or discharge in the severe cases (40 days), but higher than that in the mild cases (12 days) . However, it remains comparable to that of patients infected with H5N1, which is 18.7 days [2] . Detailed characterisation of the genes and protein structure of H7N9 has revealed a number of important features associated with the transmissibility, pathogenicity and severity of this virus in humans [12] . For the HA segment, the H7N9 hemagglutinin outer-surface proteins, from all isolates of the three clinical groups lacked a multibasic cleavage site. Despite the severity of symptoms and fatality rate in human infections, the absence of the multibasic cleavage site in HA indicates that H7N9 would be considered an LPAI virus [30] . Analysis of the protein sequence of the H7N9 viruses infecting humans showed substitutions in the HA at amino acid residues Q226L and G186V. These substitutions have been associated with increased affinity for a-2,6-linked sialic acid receptors, indicating that the virus may adapt more rapidly to infect mammals via an enhanced capacity for mutations [4, 18] . The increased binding affinity of H7N9 compared to that of H5N1, raises concerns regarding a potential adaptation of the virus that would allow it to cause a pandemic [12] . R57K, related to the antigenic site E, and N276D, related to the antigenic site C, were the most frequently observed mutations in mild cases, but R57M, related to the antigenic site E, was the major mutation in severe cases. Interestingly, the present study is the first report of the occurrence of the novel mutations L168I and E114K, near the antibody accessible site. These new mutations were detected in isolates from mild cases; however, their significance remains unclear. Given the severity of the infection in humans, virus susceptibility to antiviral drugs is vitally important for treatment [18, 31] . Genomic sequences of virus isolates showed that 12.1 % of the fatal cases, 2.1 % of the severe cases and none of the mild cases had an R294K mutation, which affects NA inhibitors -oseltamivir, zanamivir, peramivir and laninamivir [9, 11, 29, 31] . In contrast, all H7N9 viruses assessed to date have a mutation in the M gene (S31N), which confers resistance to the commercially available adamantine drugs amantadine and rimantadine [12] . The PB2, PA and PB1 genes code for proteins that form the polymerase enzyme complex, which is necessary for viral replication [7, 12, 16] . Genomic sequencing of H7N9 isolates indicated that the substitution E627K in PB2 was more frequent in the fatal group than in the severe and mild groups, whereas N559T was found in the mild group. The PB2 amino acid change increases virulence and enhances polymerase activity in mice [16] . This, along with other factors, likely contributes to the increased severity of the disease in humans with H7N9 infection. The PB2 from H7N9 isolates from birds retained Glu at position 627 and Asp at 701, strongly suggesting that the mutation is positively selected upon replication in the human host, as reported previously for zoonotic A (H7N7) and A(H5N1) infections [1] . For PA, over 90 % of fatal cases have V100A and K356R mutations, which were present in species-associated signature positions. Interestingly, the PA-356R signature is associated with increased H5N1virulence [37] . Additional markers of adaptation to non-avian hosts or virulence were noted in the PB1-F2 gene, where 90.3 % of fatal isolates showed the 90 amino acid mutations in the PB1 gene -higher than that of non-human isolates and the mild and severe cases. There were no obvious differences in the NS1 and M1 genes among the three clinical groups. All identified virus genotypes contained at least one difference in an internal gene compared with the genotypes detected in the delta of the Yangtze River and Pearl River regions during the past three epidemic waves, indicating that those genotypes may have been generated through inter-provincial poultry trade and further reassortment with local H9N2 viruses [22, 25] . Phylogenetic analysis indicated that a total of six genotypes were circulating in human cases worldwide, but no difference in distribution of genotypes was observed among the three clinical groups in the present study [24, 42] . In the first epidemic wave of 2013, genotype A was detected in all Chinese areas, with the exception of the Pearl River areas (Guangdong Province and Hong Kong) , particularly in the Yangtze River delta, including Zhejiang, Jiangsu, Jiangxi, Shandong, Fujian and Anhui Provinces, as well as Shanghai. These findings support the hypothesis that H7N9 infection originated from this area. Furthermore, some external and internal genes exchanged and generated F genotypes in these areas in the second and third epidemic waves. In contrast, genotypes B, C and E were predominant in the Pearl River (Guangdong Province and Hong Kong) in the first, second and third epidemic waves, which might indicate that other subtypes were circulating in the human and poultry populations. The infectivity, transmissibility and pathogenicity of the different genotypes has not yet been characterised. New genotypes with higher fitness may be generated by genetic tuning in the future [33, 45] . The present study had several limitations, including the small number of sequences from mild cases (14 sequences), the small number of sequences from cases of the third epidemic wave, and the failure to find a statistical association between identified genotypes and clinical outcome due to the limited number of mild cases for each genotype. Therefore, investigations are ongoing, with eyes fully open to any new cases or episodes, in order to cover all cases not included in the present study to analyze, monitor and control this serious challenge. In conclusion, advanced age, as well as a delay in confirmation of diagnosis and start of antiviral treatment, were the greatest contributory factors to a high risk of death. Furthermore, PB2 (E627K), NA(R294K) and PA (V100A) mutations were slightly associated with an increased fatality rate. However, HA (N276D) and PB2 (N559T) were obviously related to a mild outcome. Although H7N9 isolates from all over the world presented with a variety of genotypes, there were no sequence differences in the mild, severe and fatal cases. However, there was an obvious difference in the distribution of genotypes in the regions of the Pearl and Yangtze rivers and between the three different epidemic waves. There was also evidence to suggest that genetic tuning not only mediated species switching but may also have allowed the virus to increase its adaptation to more efficiently infect humans and enable person-to-person transmission [28, 43] . However, the mechanism driving genetic tuning has not yet been fully elucidated, and no vaccine for the prevention of influenza A (H7N9) virus infection is available to date. Thus, an effective avian influenza surveillance network for discovering new infections and monitoring any genetic changes and potential adaptations in the H7N9 virus at the early stages is essential. Canada recently reported two confirmed family cases of influenza A (H7N9), imported from China. This showed that it is possible to transfer the H7N9 virus from China to the rest of the world via travel and transportation [35] . Close cooperation between the different sectors related to human health, poultry and wild birds is necessary to respond to the potential pandemic risk posed by the novel avian influenza A(H7N9) viruses [8] . A/Hong_Kong/734/2014(H7N9)-Death A/Hong_Kong/8113530/2014(H7N9)-Mil d A/Hong_Kong/5942/2013(H7N9)-Severe Death A/Hong_Kong/5581/2014(H7N9)-Mil d A/Hong_Kong/8122430/2014(H7N9)-Severe A/Hong_Kong/4495/2014(H7N9)-Severe A/Hong_Kong/8194273/2014(H7N9)-Severe /11/2013(H7N9)-Death A/Hong_Kong/5731/2014(H7N9)-Mil d A/Beijing/1/2013-Severe A/Hong_Kong/8194273/2014-Severe A/Environment/Wuxi/1/2013(H7N9) A/Zhejiang/1/2014(H7N9)-Death A/Shanghai/05/2013(H7N9)-Mild A/Zhejiang/1/2013(H7N9)-Death A/Zhenjiang/1/2013(H7N9)-Severe A/Hong Kong/5942/2013(H7N9)-Severe A/Zhejiang/7/2014(H7N9)-Mild A/Zhejiang/15/2014(H7N9)-Death A/Beijing/1/2013(H7N9)-Severe A/Zhejiang/2/2014(H7N9)-Severe A/Zhejiang/7/2014(H7N9)-Mild A/Zhejiang/1/2013(H7N9)-Death A/Hong Kong/5731/2014(H7N9)-Mild A/Nanchang/1/2013(H7N9)-Severe A/Shanghai/1/2013(H7N9)-Death A/Zhejiang/15/2014(H7N9)-Death A/Zhejiang/1/2014(H7N9)-Death A/Chicken/Zhejiang/3/2014(H7N9) A/Environment/Wuxi/1/2013(H7N9) A/Nanchang/1/2013(H7N9)-Severe A/Beijing/1/2013(H7N9)-Severe A/Zhejiang/2/2014(H7N9)-Severe A/Hong A/Zhejiang/7/2014(H7N9)-Mild A/Zhenjiang/1/2013(H7N9)-Severe A/Guangdong/02/2013(H7N9)-Mild Human infections with the emerging avian influenza A H7N9 virus from wet market poultry: clinical analysis and characterisation of viral genome Comparative epidemiology of human infections with avian influenza A H7N9 and H5N1 viruses in China: a population-based study of laboratory-confirmed cases A family cluster of three confirmed cases infected with avian influenza A (H7N9) virus in Zhejiang Province of China Adaptation of novel H7N9 influenza A virus to human receptors Clinical severity of human infections with avian influenza A(H7N9) virus, China Human infection with a novel avian-origin influenza A (H7N9) virus Identification of amino acids in HA and PB2 critical for the transmission of H5N1 avian influenza viruses in a mammalian host Importance of Internet surveillance in public health emergency control and prevention: evidence from a digital epidemiologic study during avian influenza A H7N9 outbreaks Influenza A(H7N9) virus gains neuraminidase inhibitor resistance without loss of in vivo virulence or transmissibility Limited human-to-human transmission of avian influenza A(H7N9) virus Association between adverse clinical outcome in human disease caused by novel influenza A H7N9 virus and sustained viral shedding and emergence of antiviral resistance H7N9: preparing for the unexpected in influenza Epidemiological and clinical characteristics and risk factors for death of patients with avian influenza A H7N9 virus infection from Circulation of reassortant influenza A(H7N9) viruses in poultry and humans Epidemiology of human infections with avian influenza A(H7N9) virus in China Molecular basis of replication of duck H5N1 influenza viruses in a mammalian mouse model Risk factors for influenza A(H7N9) disease-China Origin and diversity of novel avian influenza A H7N9 viruses causing human infection: phylogenetic, structural, and coalescent analyses Epidemiological, clinical and viral characteristics of fatal cases of human avian influenza A (H7N9) virus in Zhejiang Province pandemic characteristics and controlling experiences of influenza H1N1 virus 1 year after the inception in Hangzhou One family cluster of avian influenza A(H7N9) virus infection in Shandong Continuing reassortment leads to the genetic diversity of influenza virus H7N9 in Guangdong Amino acid substitutions in polymerase basic protein 2 gene contribute to the pathogenicity of the novel A/H7N9 influenza virus in mammalian hosts Evolution of the H9N2 influenza genotype that facilitated the genesis of the novel H7N9 virus Continuous reassortments with local chicken H9N2 virus underlie the human-infecting influenza A (H7N9) virus in the new influenza season, Guangdong Probable person to person transmission of novel avian influenza A (H7N9) virus in Eastern China, 2013: epidemiological investigation Differences in the epidemiology of human cases of avian influenza A(H7N9) and A(H5N1) viruses infection Infection with possible precursor of avian influenza A(H7N9) virus in a child R292K substitution and drug susceptibility of influenza A(H7N9) viruses The multibasic cleavage site of the hemagglutinin of highly pathogenic A/Vietnam/1203/2004 (H5N1) avian influenza virus acts as a virulence factor in a hostspecific manner in mammals The emergence of influenza A H7N9 in human beings 16 years after influenza A H5N1: a tale of two cities Comparison of patients hospitalized with influenza A subtypes H7N9, H5N1, and 2009 pandemic H1N1 Genetic tuning of the novel avian influenza A(H7N9) virus during interspecies transmission Overview of the emergence and characteristics of the avian influenza A(H7N9) virus Avian influenza (H7N9) virus infection in Chinese tourist in Malaysia Transmission of avian influenza A(H7N9) virus from father to child: a report of limited person-to-person transmission PA-356R is a unique signature of the avian influenza A (H7N9) viruses with bird-to-human transmissibility: potential implication for animal surveillances Family clusters of avian influenza A H7N9 virus infection in Guangdong Province Human H7N9 avian influenza virus infection: a review and pandemic risk assessment Clinical characteristics of 26 human cases of highly pathogenic avian influenza A (H5N1) virus infection in China Human infection with avian influenza A H7N9 virus: an assessment of clinical severity Influenza H7N9 and H9N2 viruses: coexistence in poultry linked to human H7N9 infection and genome characteristics A(H7N9) virus results in early induction of proinflammatory cytokine responses in both human lung epithelial and endothelial cells and shows increased human adaption compared with avian H5N1 virus Dual E627K and D701N mutations in the PB2 protein of A(H7N9) influenza virus increased its virulence in mammalian models Genetic tuning of avian influenza A (H7N9) virus promotes viral fitness within different species Acknowledgements This work was supported by grants from Provincial Medical Research Fund of Zhejiang, China (2014C13G2010106). We acknowledge the authors who submitted the sequences used in this study to GenBank and GISAID. Conflict of interest None declared.