Spatial characteristics and the epidemiology of human infections with avian influenza A(H7N9) virus in five waves from 2013 to 2017 in Zhejiang Province, China
Spatial characteristics and the epidemiology of human infections with avian influenza A (H7N9) virus in five waves from 2013 to 2017 in Zhejiang Province, China
Haocheng Wu 0 1
XinYi Wang 0 1
Ming Xue 1
Melanie Xue 1
Chen Wu 0 1
Qinbao Lu 0 1
Zheyuan Ding 0 1
Xiaoping Xv 0 1
Junfen Lin 0 1
0 Zhejiang Province Center for Disease Control and Prevention , Hangzhou, Zhejiang Province , China , 2 Key Laboratory for Vaccine, Prevention and Control of Infectious Disease of Zhejiang Province , Hangzhou, Zhejiang Province , China , 3 Hangzhou Centre for Disease Control and Prevention , Hangzhou, Zhejiang, Province , China , 4 Kingston University UK , London , United Kingdom
1 Editor: Florian Krammer, Icahn School of Medicine at Mount Sinai , UNITED STATES
The data for the H7N9 cases from Zhejiang Province between 2013 and 2017 was obtained from the China Information Network System of Disease Prevention and Control. The start date of each wave was 16 March 2013, 1 July 2013, 1 July 2014, 1 July 2015 and 1 July 2016. The F test or Pearson's chi-square test were used to compare the characteristics of the five waves. Global and local autocorrelation analysis was carried out to identify spatial autocorrelations. Ordinary kriging interpolation was analyzed to estimate the number of human infections with H7N9 virus and to infer the extent of infections in the areas with no cases reported before.
Data Availability Statement: All relevant data are
within the paper and its Supporting Information
Funding: This work was supported by grants from
Provincial Medical Research Fund of Zhejiang,
China (grant number were 2015RCB008 and
2017KY131). The document in website were http://
8DD1BAEA38B99297 and http://www.zjmed.org.
The five-wave epidemic of H7N9 in China emerged in the second half of 2016. This study
aimed to compare the epidemiological characteristics among the five waves, estimating the
possible infected cases and inferring the extent of the possible epidemic in the areas that
have not reported cases before.
There were 45, 94, 45, 34 and 80 cases identified from the first wave to the fifth,
respectively. The death rate was significantly different among the five waves of epidemics (χ2 =
10.784, P = 0.029). The age distribution (F = 0.903, P = 0.462), gender (χ2 = 2.674, P =
0.614) and occupation(χ2 = 19.764, P = 0.407) were similar in each period. Most of the
cases were males and farmers. A significant trend (χ2 = 70.328, P<0.001) was identified that
showed a growing proportion of rural cases. There were 31 high-high clusters and 3
highlow clusters at the county level among the five waves and 12, 8, 2, 9 and 3 clusters in each
funders had no role in study design, data collection
and analysis, decision to publish, or preparation of
Competing interests: The authors have declared
that no competing interests exist.
wave, respectively. The total cases infected with the H7N9 virus were far more than those
that have been reported now, and the affected areas continue to expand. The epidemic in
the north of Zhejiang Province persisted in all five waves. Since the second wave, the virus
spread to the south areas and central areas. There was an obvious decline in the infected
cases in the urban areas, and the epidemics mostly occurred in the rural areas after the
fourth wave. The epidemic was relatively dispersed since the third wave had fewer than two
cases in most of the areas and showed a reinforcing trend again in the fifth wave.
The study revealed that there were few differences in the epidemiologic characteristics
among the five waves of the epidemic. However, the areas where the possible epidemic
circulated was larger than reported. The epidemic cross-regional expansion continued and
mostly occurred in rural areas. Continuous closure of the live poultry market (LPM) is
strongly recommended in both rural and urban areas. Illegal and scattered live poultry
trading, especially in rural areas, must be forbidden. It is suggested too that a more rigorous
management be performed on live poultry trade and wholesale across the area. Health
education, surveillance of cases and pathogenicity should also be strengthened.
The H7N9 virus is a virus that reasserts multiple times and whose gene fragments are obtained
from H7N9, H9N2 and H7N3 subtypes of the influenza A virus[1±3]. The first case of human
infections with the H7N9 virus occurred in the spring of 2013 in Eastern China[
epidemic caused great concern due to the increasing number of cases, the expansion of the
affected areas, the high fatality rate, mutation of the virus and the stigma attached to the virus
[6,7]. To date, five waves of H7N9 epidemics have emerged in mainland China. There were
775 laboratory-confirmed infections of A(H7N9) virus across 16 provinces and 3
municipalities as of August 31, 2016[
]. In addition, 23 travelers infected with A(H7N9) exported the
virus to Hong Kong (sixteen cases), Taiwan (four cases), Canada (two cases), and Malaysia
(one case), leading to four deaths [
]. The fifth wave of the epidemic occurred in autumn 2016,
and incidences far exceeded those in the corresponding period of 2015. It should be noted that
Zhejiang is a province located in the Yangtze River Delta region of southeastern China that is well recognized as the original source of the H7N9 outbreaks . This province was the most seriously affected area, accounting for nearly 30% of total cases reported and 25% of the total fatal cases during the last four waves in mainland China.
There are likely to be many more cases than those reported due to there being an unknown
number of mild and subclinical infections [
]. The objectives of this study are to identify the
epidemiological characteristics and the distribution of H7N9 virus in human infection from
2013 to 2017. The Kriging spatial interpolation methods will be applied to estimate the number
of human infections with H7N9 virus and to infer the extent of infections in areas that
currently have no reported cases.
Definition of the five waves
Based on the date of onset, the first wave of H7N9 virus circulation occurred from 13 March to 30 June 2013. In the first wave, the first and last cases occurred on 13 March and 18 April
2 / 17
2013, respectively. The second wave occurred from 1 July 2013 to 30 June 2014, with the first
and last cases occurring on 7 October 2013 and 3 June 2014, respectively. The third wave
occurred from 1 July 2014 to 30 June 2015. In the third wave, the first and last cases occurred
on 17 November 2014 and 28 May 2015, respectively. The fourth wave occurred from 1 July
2015 to 30 June 2016. In the fourth wave, the first and last cases occurred on 18 September
2015 and 24 June 2016, respectively. The fifth wave occurred from 1 July 2016 to 31 March
2017, with the first case occurring on 28 September 2016.
Patient/cluster definition and data collection
The patient definition. The diagnosis of the infections confirmed with the H7N9 virus
was based on the Chinese Guideline of Diagnosis and Treatment for Human Infections with
the Avian Influenza A(H7N9) Virus issued by the National Health and Family Planning
Commission of the People's Republic of China[
]. The definition of the confirmed case is
clinical manifestation with acute influenza (fever, cough, coryza, difficulty breathing) or a
history of contact with a confirmed or suspected case and a laboratory test that includes subtype
confirmation by several parameters. These include PCR, viral isolation or no less than a
fourfold increase in virus-specific serum antibodies isolated from paired positive sera samples
The cluster definition. According to the Chinese Guideline of Epidemic Prevention and
Control for Human Infections with the Avian Influenza A(H7N9) Virus issued by the National
Health and Family Planning Commission of the People's Republic of China, the definition of the cluster is that two confirmed cases emerged within 7 days and in small areas, such as one family or one community.
The data collection. The confirmed case must be reported through the China Informa
tion Network System of Disease Prevention and Control by medical staff. The data of H7N9
cases in Zhejiang Province from 2013 to 2017 were obtained from this network system(S1
Profile of Zhejiang Province. The Zhejiang Province is in the southeast China between
longitudes 118oE-123oE and latitudes 27oN-32oN. There are two sub-provincial cities
Zhou and Ningbo) and nine prefecture-level cities, including Wenzhou, Huzhou, Jiaxing,
Shaoxing, Jinhua, Zhoushan, Quzhou, Taizhou and Lishui, which cover 90 counties(Fig 1) in
Spatial autocorrelation analysis
A global and local autocorrelation analysis using the Moran's I was carried out to identify
spatial autocorrelations. According to the theory, the closer two locations were to each other, the
more likely for the incidence rates have an impact on each other [
]. The global
autocorrelation demonstrates the spatial distribution clusters over areas [
]. The Local Moran's Index
was used to classify the autocorrelations into positive and negative values. When there were
similar high values or low values of the incidence rates, they were defined as having a positive
autocorrelation (represented as High-High or Low-Low autocorrelation). If the attributes held
opposing high and low values, they were considered to have a negative autocorrelation
(represent as High-Low or Low-High autocorrelation) [
Ordinary kriging interpolation
Ordinary kriging (OK) is one of the most commonly used kriging interpolation methods. It has been proven to be the best linear unbiased estimate for variables in regions where the data are spatially autocorrelative. The unsampled value is computed as a linear combination of
3 / 17
Fig 1. Maps of Zhejiang Province, China with area name.
where z (x0,y0) is the predicted value that is spatially located at point (x0,y0), λi is the weight
associated with the measured value of Z at the location (xi,yi), and the weights are derived
from the kriging equation using a semivariance function. The parameters of the semivariance
function and the nugget effect can be estimated by an empirical semivariance function[
xi hx; yi hy
where γ(h) is the semivariance value at the distance interval h, N(h) is the number of sample
pairs within the distance interval h, and Zobs(xi + hx,yi + hy) Zobs(xi,yi) are sample values at two
points separated by the distance interval h.
4 / 17
The F test was employed to compare the mean age of H7N9 cases between the five waves.
Pearson's chi-square or Fisher's exact test were used to compare the occupational and sex
distribution of the cases between the five waves. Pearson's chi-square test and the chi-square test for
trend were used to compare the death rate and the proportion of rural cases between the five
waves. The Bonferroni t-test was employed to further compare the two indexes between each
wave. The methods above were computed by SAS9.2 (SAS Institute Inc., Cary, NC). ArcGIS
software (version 10.1, SERI Inc.; Redlands, CA, USA) was used for the spatial autocorrelation
analysis and the Ordinary kriging interpolation. The county was adopted as the geographic
unit to calculate the spatial autocorrelation. A P value of less than 0.005 was considered to be
statistically significant for the Bonferroni t-test and a P value of less than 0.05 for the other
This study was reviewed and approved by the Ethics Committee of the Zhejiang Provincial
Center for Disease Control and Prevention. All the data for the individuals were kept confidential as requested.
Epidemiologic characteristics in five waves
From March 2013 to March 31, 2017, there were a total of 298 H7N9 human infections and
112 fatal cases reported in the Zhejiang Province, China. There were 45, 94, 45, 34 and 80 cases
identified from the first wave to the fifth epidemic, respectively. In each wave, 10, 39, 24, 13,
and 26 fatalities were found, respectively, with significant differences among them (χ2 =
10.784, P = 0.029). There was no significant difference among the latest four waves (χ2 = 2.662,
P = 0.447), and the significant difference at the 0.005 level was only identified between the first
and the third epidemic (χ2 = 9.265, P = 0.002). The highest death rate was in the third wave
(53.3%) and the lowest was in the first (22.2%). However, there were no significant trend in
the death rate among the five waves (χ2 = 0.034, P = 0.854).
There was a total of 144 rural cases reported, with significant differences in the proportion
among the five waves (χ2 = 73.361, P<0.001). In total, 10, 24, 21, 24 and 65 rural cases were
found in each wave, with a significant trend (χ2 = 70.328, P<0.001) identified showing a
growing proportion of rural cases. The proportion in the fifth wave was significantly different from
that in the first, the second and the third wave (χ2 = 41.811, 53.697, 16.047, P<0.001, P<0.001,
P<0.001, respectively), but similar to the fourth wave (χ2 = 1.584, P = 0.208). All the above
results were shown in Table 1.
An obvious seasonality was present for the occurrence of H7N9 human infection across the
five waves (Fig 2). Cases emerged in March 2013 and grew rapidly, with a peak in April 2013,
and then declined. No case was reported from May to September until a new case remerged in
October 2013. It should be noted that most cases occurred in January in the following three
waves, accounting for an average of 58.4% of the total cases in each wave. The peak incidence
of the fifth wave occurred in December 2016, which was earlier than the other waves, and
continued to January, 2017, followed by a sharp decline in February and March (7 and 3 cases,
respectively). A few cases (9 cases, 4.13% of total) were identified during May and September
in the previous four waves. Of the total cases, male cases predominated, accounting for an
average of 63.1% in the five waves. Gender distribution did not differ significantly among the
five waves (χ2 = 2.674, P = 0.614). The mean age of the cases was 57.2±0.88 in the five waves.
5 / 17
Spatial clustering distribution
According to the global autocorrelation analysis (Table 2), the distribution of the cases infected
with the H7N9 virus in the first, second and fourth waves presented a spatial autocorrelation.
The Moran's I Index was highest in the first wave (Index = 0.8161, P value<0.001), followed
by an index of 0.2658 in the second wave and 0.1468 in the fourth one. Spatial autocorrelation
was not found in the third and fifth waves among the whole province.
Based on the result of the local Moran's I autocorrelation (Table 3), there were 31 high-high
clusters and 3 high-low clusters in total at the county level through the five waves, with 12, 8, 2,
9 and 3 clusters in each wave, respectively. It should be noted that Anji county, Deqing county,
Yuhang county, Xiaoshan county, Shangcheng county, Shaoxing county and Yuecheng county
had long-term high-high clusters throughout the five waves (4,3,3,3,2,2,2 high-high clusters
These counties were all located in the northern areas of Zhejiang Province (Fig 3). There
were 29 counties in this region, which accounted for 85.3% of all the areas with high-high or
high-low clusters located in the northern areas of Zhejiang Province. What is interesting is
that the high clusters emerged in the southern areas (Wencheng and Yueqing county) of
Zhejiang Province for the first time in the fifth wave. It should also be noted that the area of the
counties varied in size, with the largest being 4452 square kilometers and the smallest being
18.17 square kilometers.
The geostatistical parameters for interpolation
It was shown in the semivariogram figure that the value of the semivariation function was rela
tively low over a short range, gradually increased as the distance between the sample points
grew and then fluctuated around an extremum when the distance grew to a particular value.
Similar semivariogram characteristics were shown among the five waves (Fig 4).
6 / 17
Fig 2. The epidemic curve of H7N9 human infections in Zhejiang Province, China, March
2013-Janurary 2017 at monthly interval.
The partial sill (C1) was 0.2582, 7.7863, 0.6320, 0.9723, and 1.9801 in each wave,
respectively. The nugget (C0) was 0.4789, 0.1642, 0.0982, 0.1416 and 0.2161 in each corresponding
wave, respectively. The C1/ C0+ C1 ratios among the last four waves were all above 85%, and
were even above 90% in the second and fifth waves, but was relatively low (35%) in the first
wave, which suggested that the spatial autocorrelation was generally relatively strong. All the
results are shown in Table 4.
7 / 17
Estimation of the extent
The result of the ordinary kriging interpolation showed that there were many unidentified
infections in the areas with no reported cases, indicating that the total number of cases infected
with the H7N9 virus were far more than what has been reported to date (Fig 5). In most of the
areas, there were no more than two cases within one epidemic period, and based on the
estimation of the confidence interval, the number of cases were between 0±2 in the first, the third
and the fourth wave. In contrast, the cases in the second and fifth wave were more than those
in the aforementioned three waves and fluctuated between 0±6 in the second wave and
between 0±4 in the fifth wave in most of the areas (Figs 6 and 7).
Second, the epidemics continuously circulated in the northern areas, including Hangzhou
City, Jiaxing City, Huzhou City and Shaoxing City in Zhejiang Province throughout the five
waves (Figs 5±7), and spread to the southern and central areas since the second wave (Figs 5
and 7). The estimated number of cases in the northern areas were obviously more than those
in the other areas, especially in the second wave, where the cases in Xiaoshan county, Yuhang
county, Deqing county, Anji county and Yinzhou county were all above 7 (Fig 7). Ningbo
City, Taizhou City and Wenzhou City in the southeast areas were affected in the second wave,
and Jinhua City, Quzhou City and Lishui City located in the south-central areas were also
affected in the waves that followed. The affected areas in the fifth wave were larger than those
in previous waves. It was found that most cities in Zhejiang Province were affected by the
epidemics throughout all five waves, but only Zhoushan City, located on the east island, remained
virus free (Figs 5 and 6), and there were perhaps no more than two cases in Zhoushan City,
according to the upper-bound of the estimation (Fig 7). It should also be noted that one
highly-infected region in the east area of Ningbo City and one highly infected region in the
south were identified in the second and the fifth wave, respectively.
Third, according to the estimation maps of H7N9 human infections, there was an interesting result that the cases that emerged in the urban areas declined significantly. Since the fourth wave, most cases occurred in rural areas, while in the previous three waves both urban and rural areas were influenced (Fig 5).
Lastly, the spatial aggregation enhanced stepwise relative to the previous two waves based
on the definition of the cluster[
]. The epidemic presented some relatively dispersed
characteristics since the third wave with cases reported fewer than two in most counties, but the
epidemic showed a reinforcing trend again in the fifth wave (Figs 5 and 7). In general, the affected
areas continued to expand, and the epidemic may have been circulated in some new areas with
no previously reported cases (Figs 5 and 7).
Similar to previous research [3,10], most of the cases infected with H7N9 virus were males and people over 50 years old, and there was a higher proportion of farmers. The epidemic
8 / 17
characteristics, including gender, age and occupational distribution, were similar across the
five waves in the Zhejiang Province in China. Although there were significant differences in
death rate, the fatality rate was similar among the last four waves and there was no significant
trend in death rate across the five waves. The reason for these similarities in the epidemic
characteristics among the five waves may be because the circulating H7N9 virus genotypes were
similar to those in previous waves, although there were a few point mutations discovered in
the amino acids in some samples, which might increase the virus' adaptation to gain
personto-person transmission ability and infect humans more efficiently [21±23]. The conclusion to
date was that the genetic changes have not been sufficient enough to alter the antigenic
characteristics or cause sustained human-to-human transmission[
]. In the fifth wave, there were
more cases occurring in December 2016 than in the corresponding periods in previous years.
Its peak incidence continued in January 2017 and then fell rapidly. Meanwhile, the incidence
9 / 17
Fig 3. Maps of the local autocorrelation analysis of H7N9 human infections in Zhejiang Province among five waves by Local Moran'I.
peak may move up. The reasons for the advance of the peak could be that the temperature in
December 2016 was suitable for the virus and the advance of the Chinese New Year in 2017, which increased live poultry consumption and exposure.
The key assumption of the interpolation analysis was that the distribution of the value was a
]. As we know, most of the infectious diseases were spatially correlative
depending on the geographical environment, meteorological factors, economic levels, local
custom and other factors. The H7N9 epidemic emerged in the Northern areas of the Zhejiang
Province, then the spread of the virus and the expansion of the contaminated region contributed to the dispersion of the epidemic. The similarity of the contaminated environment and the exposure to the virus make the case spatially correlative and the interpolation analysis feasible.
Additionally, we used some parameter to evaluate the assumption of the model. The semivariance parameter is usually used to evaluate the strength of the correlation, except for Moran's I Index. According to the semivariogram figure, the semivariance was relatively low in the short range and gradually grew to the extremum among the five waves. This means that the
10 / 17
Fig 4. The figure of semivariogram of H7N9 human infections in Zhejiang Province among five waves.
more close sample points, the characteristics of the sampleswere more similar. The ratio of the
partial sill to the sill was the key indicator, which represents the extent of variation due to
spatial autocorrelation. In the study, these ratios among the last four waves were all above 85%,
which suggested that the spatial autocorrelation was relatively strong in general and the
variation due to the spatial autocorrelation was the main contribution to the observed difference.
The ratio in the first wave was relatively low(35%), the reason of which maybe the cases in the
first wave almost clustered in the northern areas and the different epidemics in the northern
areas relative to other areas, however the local spatial autocorrelation in the northern areas was
still strong. Therefore, the interpolation analysis was suitable for the incidence of H7N9
infection among the five waves.
Ordinary kriging interpolation was widely used to estimate the extent and geographical dis
tribution for infectious diseases such as HIV, West Nile Disease and Dengue fever [25±27].
Kriging is an effective supplement for regular case surveillance, which depends on the medical personnel report. According to the inference of some previous research, there were mild and subclinical infections occurring at the same time as the epidemic. In addition to the severe cases, the number of cases of symptomatic H7N9 virus infections may be approximately 500±
Ratio(C1/ C0+ C1)
11 / 17
Fig 5. Maps of the estimation of the H7N9 human infections in Zhejiang Province among five waves epidemics by ordinary kriging
1000 times that of the laboratory-confirmed cases [28±30]. Consistent with the previous
research, it was found that areas with possible unidentified cases were larger. The number of
cases could be more than those reported to date and most of them may be the mild or
The epidemic is continuously circulating in the northern areas in Zhejiang Province, China
throughout the five waves. Most of the long-term high-high clusters were identified in the
northern areas, which meant the incidences of h7n9 infection were similarly high among these
areas, and more attention is needed to prevent epidemics in the areas. The possible reason for
the continuous epidemic and clusters in the northern area is that these areas are in the Yangtze
River Delta region along the Taihu Lake, which is recognized as the origin of virus
recombination that lead to the outbreaks [
]. Therefore, the circulation of the virus in the domestic
poultry transmitted by the local wildfowl and contaminated environment sustained, which
contributes to the enzootic H7N9 virus infection in such areas. Although permanent closure of
the live poultry markets was implemented in the main urban areas of Zhejiang Province, live
poultry trade in rural areas and in privacy or wholesale crossing the border still existed. This
contributed to the spread of the virus from the northern contaminated region to other counties
and the expansion of the affected areas, especially the southern and central areas of Zhejiang
Province [9,31]. Based on our study, Wenzhou city, located in the southeast area of Zhejiang,
12 / 17
Fig 6. Maps of the lower-bound of confidence interval(α = 0.05) of the H7N9 human infections in Zhejiang Province
among five waves epidemics by ordinary kriging interpolation.
possibly became another subsource of H7N9 outbreaks. One high-high cluster was found in
the Wencheng county in the fifth wave, which meant similar epidemics in this county and its
neighboring areas. Additionally, 3 high-low clusters were found in the southeast areas.
Although incidences in their surrounding areas were lower, the epidemic may expand to the
neighboring areas where it may become a source of H7N9 outbreak in the future. One
important thing that should be noted is that the continuous trans-regional expansion of the H7N9
viruses has accelerated the reassortment among H7N9 and other subtypes of viruses , which
would make the viruses more mammal-adapted as well as allow them to become
human-tohuman transmissible and lead a pandemic in the future [
It was proven that the closure of the live poultry markets (LPM), which were the main place
of live poultry transaction in urban areas, was highly effective in reducing the risk of H7N9
infection in humans [
]. In our study from the estimation maps, it was found that the
epidemics tended to disappear in urban areas and presented relatively sporadic characteristics,
especially since the fourth wave after the closure of the LPM in the urban areas in 2014.
Meanwhile, according to the total proportion of rural cases, it increased obviously in the fourth and
the fifth waves. However, different from urban areas in Zhejiang Province, live poultry
transactions in rural areas were mainly scattered among the whole area and not centralized in
markets. For example, many farmers sell poultry at the roadside instead of in the market. This may
13 / 17
Fig 7. Maps of the upper-bound of confidence interval(α = 0.05) of the H7N9 human infections in Zhejiang Province among five
waves epidemics by ordinary kriging interpolation.
be one reason for the increasing proportion of rural cases. According to the latest research
], raising backyard poultry at home and having direct contact with backyard poultry have
been identified as two new risk factors associated with increased risk of infection with H7N9
virus, which may be another reason that the epidemic currently circulated primarily in the
rural areas in the Zhejiang Province instead of in urban areas.
Since the third wave, the affected areas expanded and the rural cases increased abruptly,
especially in the fifth wave, although the cases were fewer than two in most of the rural areas.
However, because of the active live poultry trading and the custom of raising backyard poultry,
the virus is very easily transmitted between the village, which leads to the expansion of the
contaminated environment and the circulation of the epidemic in new areas. This may increase
the risk of exposure to the virus and the possibility of a potential cluster. Additionally, this
maybe the reason for the reinforcing trend in the fifth wave.
Several limitations should be noted with our study. First, data, including the clinic informa
tion, environmental specimens and serological surveillance, were not collected, and thus
detailed clinical characteristics, positive rates of H7N9 virus in the environment and
seropositive rates could not be compared. Second, only those cases confirmed to be infected with the
H7N9 virus were analyzed by the ordinary kriging interpolation methods. If the seropositive
14 / 17
rate of the population and the positive rate of H7N9 in the environment were added into the
analysis, the estimation would have been more accurate. Finally, the extent of infection was
estimated only in the spatial as opposed to the temporal dimension. Therefore, the temporal
trend of the epidemic in each area could not be identified. Further research should be carried
out to obtain a more comprehensive conclusion.
In conclusion, this study revealed that there were few differences in the epidemiologic
characteristics among the five waves of the H7N9 infection epidemics. The epidemic was relatively
dispersed, but areas with a possible circulating epidemic were larger than what has been
reported, especially in the fifth wave. The epidemic mostly occurred in rural areas, and
crossregional expansion continued since the fourth wave. Based on the findings, we suggest
continuous closure of the LPM in both rural and urban areas. Illegal and scattered live poultry
trading, especially in rural areas, must be forbidden. In addition, we should strengthen health
education about self-protection from feeding backyard poultry and slaughtering poultry by
hand. Moreover, the management of the live poultry trade and wholesale border crossing
should be more rigorous to prevent the spread of the virus across the area. The surveillance of
the cases, especially the mild ones, should be more heightened. Lastly, ongoing research of the
pathogenicity in the internal genes is needed.
S1 File. The data base of key information of the case infected with H7N9 virus among five
waves in Zhejiang Province.
S2 File. The certification of the manuscript was edited for proper English language.
We thank all the clinical personnel for investigating and reporting the information about the cases infected with the H7N9 virus. We very much appreciate the 90 county-level CDCs for their investigation and correction of the epidemiological information.
Conceptualization: Xiaoping Xv, Junfen Lin.
Data curation: Haocheng Wu, Chen Wu, Qinbao Lu, Zheyuan Ding.
Formal analysis: Haocheng Wu, Chen Wu.
Methodology: Haocheng Wu.
Resources: Chen Wu, Qinbao Lu, Zheyuan Ding.
Supervision: XinYi Wang, Xiaoping Xv, Junfen Lin.
Writing ± original draft: Haocheng Wu.
Writing ± review & editing: XinYi Wang, Ming Xue, Melanie Xue, Junfen Lin.
15 / 17
16 / 17
1. Chen Y , Liang W , Yang S , Wu N , Gao H , Sheng J , et al. Human infections with the emerging avian influenza A H7N9 virus from wet market poultry:clinical analysis and characterisation of viral genome . Lancet . 2013 ; 381 : 1916 ±25. https://doi.org/10.1016/S0140- 6736 ( 13 ) 60903 - 4 PMID: 23623390 .
2. Liu D , Shi W , Shi Y , Wang D , Xiao H , Li W , et al. Origin and diversity of novel avian influenza A H7N9 viruses causing human infection: phylogenetic,structural, and coalescent analyses . Lancet . 2013 ; 381 : 1926 ±32. http://dx.doi.org/10.1016/S0140- 6736 ( 13 ) 60938 - 1 . PMID: 23643111 .
3. He F , Chen EF , Li FD , Wang XY , Wang XX , Lin JF . Human infection and environmental contamination with Avian Influenza A(H7N9) Virus in Zhejiang Province, China:risk trend across the three waves of infection . BMC Public Health , 2015 ; 15 : 931 . https://doi.org/10.1186/s12889-015-2278-0 PMID: 26392274 .
4. Li Q , Zhou L , Zhou M , Chen Z , Li F , Wu H , et al. Epidemiology of human infections with avian influenza A(H7N9) virus in China . N Engl J Med . 2014 ; 370 ( 6 ): 520 ± 32 . https://doi.org/10.1056/NEJMoa1304617 PMID: 23614499 .
5. Cowling BJ , Jin L , Lau EH , Liao Q , Wu P , Jiang H , et al. Comparative epidemiology of human infections with avian influenza A H7N9 and H5N1 viruses in China: a population-based study of laboratory-confirmed cases . Lancet . 2013 ; 382 ( 9887 ): 129 ± 37 .http://dx.doi.org/10.1016/S0140- 6736 ( 13 ) 61171 - X . PMID: 23803488 .
6. Liu J , Wu T , Liu D , Qi X , Shi Y , Gao GF , et al. H7N9: a low pathogenic avian influenza A virus infecting humans . Current Opinion in Virology 2014 ; 5 : 91 ± 97 . https://doi.org/10.1016/j.coviro. 2014 . 03 .001 PMID: 24705093 .
Wang XY , Chai CL , Li FD , He F , Yu Z , Wang XX , et al. Epidemiology of human infections with avian influenza A(H7N9) virus in the two waves before and after October 2013 in Zhejiang province , China. Epidemiol. Infect . 2014 ; 143 ( 9 ): 1839 ±45. https://doi.org/10.1017/S095026881400257X PMID: 25286879 .
8. Xiang NJ , Li XY , Ren RQ , Wang DY , Zhuo SZ , Carolyn M. Greene , et al. Assessing change in avina influenza A(H7N9) virus infections during the fourth epidemic- China , September 2015 - August 2016 . MMWR. 2016 ; 65 ( 49 ): 1390 ± 94 .Available: https://www.altmetric.com/details.php? domain = www.cdc. gov&citation_id=14772368. https://doi.org/10.15585/mmwr.mm6549a2 PMID: 27977644 9.
Wang DY , Yang L , Zhu WF , Zhang Y , Zou SM , Bo H , et al. Two Outbreak Sources of Influenza A (H7N9) Viruses Have Been Established in China . Journal of Virology . 2016 ; 90 ( 12 ): 5561 ± 73 . https:// doi.org/10.1128/JVI.03173-15 PMID: 27030268 .
10. Xiang NJ , Luliano A.D , Zhang YP , Ren RQ , Geng XY , Ye BL , et al. Comparison of the first three waves of avian influenza A(H7N9) virus circulation in the mainland of the People's Republic of China . BMC Infectious Disease . 2016 ; 16 : 734 . https://doi.org/10.1186/s12879-016-2049-2 PMID: 27919225 .
11. Xiang NJ , Bai T , Kang K , Yuan H , Zhou SZ , Ren RQ , et al. Sero-epidemiologic study of influenza A (H7N9) infection among exposed populations , China 2013-2014.Influenza Other Respi Viruses . 2016 . https://doi.org/10.1111/irv.12435 PMID: 27762061 .
12. National Health and Family Planning Commission of the People's Republic of China. Chinese guideline of diagnosis and treatment for human infections with the avian influenza A(H7N9) virus (2nd edn, 2013 ) (http://www.nhfpc.gov.cn/yzygj/s3585u/201304/7e2ad4cdf98b4e2285eab1c15ded8370.shtml). Accessed 10 April 2013 .
13. National Health and Family Planning Commission of the People's Republic of China. Chinese guideline of diagnosis and treatment for human infections with the avian influenza A(H7N9) virus (2014 edition) (http://www .nhfpc.gov.cn/yzygj/s3593g/201401/3f69fe196ecb4cfc8a2d6d96182f8b22.shtml). Accessed 24 January 2014 .
14. National Health and Family Planning Commission of the People's Republic of China. Chinese guideline of epidemic prevention and control for human infections with the avian influenza A(H7N9) virus(3rd edn, 2014 ) (http://www.nhfpc.gov.cn/jkj/s3577/201401/8c1828375a7949cd85454a76bb84f23a.shtml). Accessed 27 January 2014 .
15. Liu XN , Huang F , Wang P . Spatial analysis theories and methods on GIS . Beijing: Science Press. 2008 ; 189 ±194 p.
16. Ashraf M. Dewan1 , Corner Robert1 , Hashizume Masahiro , Ongee Emmanuel T. Typhoid Fever and Its Association with Environmental Factors in the Dhaka Metropolitan Area of Bangladesh: A Spatial and Time-Series Approach . PLOS Neglected Tropical Diseases. 2013 ; 7:e1998 . https://doi.org/10.1371/ journal.pntd. 0001998 PMID: 23359825.
17. Tang FY , Cheng YJ , Bao CJ , Hu JL , Liu WD , Liang Q , et al. Spatio-Temporal Trends and Risk Factors for Shigella from 2001 to 2011 in Jiangsu Province, People's Republic of China . PLOS ONE . 2014 ; 9: e83487 . https://doi.org/10.1371/journal.pone.0083487 PMID: 24416167 .
18. Qu LL , Xiao HY , Zheng NJ , Zhang ZY , Xu Y. Comparison of four methods for spatial interpolation of estimated atmospheric nitrogen deposition in South China . Environ Sci Pollut Res . 2016 ; https://doi.org/ 10.1007/s11356-016-7995-0 PMID: 27826827 .
19. Fu CC , Zhang HB , Tu C , Li LZ , Luo YM . Geostatistical interpolation of available copper in orchard soil as influenced by planting duration . Environ Sci Pollut Res . 2016 ; https://doi.org/10.1007/s11356-016- 7882-8 PMID: 27798802 .
20. Hu Y , Gao J , Meina Chi , Luo C , Henry Lynn , Sun LQ , et al. Spatio-Temporal Patterns of Schistosomiasis Japonica in Lake and Marshland Areas in China:The Effect of Snail Habitats . Am.J.Trop.Med . Hyg. 2014 ; 91 ( 3 ): 547 ± 54 . https://doi.org/10.4269/ajtmh.14-0251 PMID: 24980498 .
21. Chen GW , Kuo SM , Yang SL , Gong YN , Hsiao MR , Liu YC , et al. Genomic Signatures for Avian H7N9 Viruses Adapting to Humans . PLOS ONE . 2016 ; 11 ( 2 ):e0148432. https://doi.org/10.1371/journal.pone. 0148432 PMID: 26845764 .
22. Sun Y , Wang HH , Lou XY , Mao HY , Yan Y , Chen EF , et al. Living poultry markets in rural area: Human infection with H7N9 virus re-emerges in Zhejiang Province, China, in winter 2014 . Journal of Clinical Virology . 2015 ; 70 : 16 ± 22 . http://dx.doi.org/10.1016/j.jcv. 2015 . 06 .092.PMID: 26305813 .
23. Sha JP , Chen XW , Ren YJ , Chen HJ , Wu ZQ , Ying D , et al. Differences in the epidemiology and virology of mild, severe and fatal human infections with avian influenza A (H7N9) virus . Arch Virol . 2016 ; 161 : 1239 ± 1259 . https://doi.org/10.1007/s00705-016-2781-3 PMID: 26887968 .
24. Wang XX , Liu SL , Mao HY , Yu Zhao , Chen EF , Chai CL . Surveillance of Avian H7N9 Virus in Various Environments of Zhejiang Province, China before and after Live Poultry Markets Were Closed in 2013± 2014 . PLOS ONE . 2015 ; 10 ( 8 ):e0135718. https://doi.org/10.1371/journal.pone.0135718 PMID: 26308215 .
25. Carrel Margare , Janko Mark, Melchior Kashamuka Mwandagalirwa Camille Morgan, Fwamba Franck, Muwonga JeÂreÂmie , et al. Changing spatial patterns and increasing rurality of HIV prevalence in the Democratic Republic of the Congo between 2007 and 2013 .Health & Place. 2016 ; 39 : 79 ± 85 . http://dx. doi.org/10.1016/j.healthplace. 2016 . 02 .009.PMID: 26974234 .
26. Carrie F. Nielsen M. Armijos Veronica , Wheeler Sarah, Carpenter Tim E., Boyce Walter M. , Kelley Kara , et al. Risk Factors Associated with Human Infection during the 2006 West Nile Virus Outbreak in Davis, a Residential Community in Northern California . Am J Trop Med Hyg . 2008 ; 78 ( 1 ): 53 ± 62 .PMID: 18187785 .
27. Barbosa Gerson Laurindo, DonalõÂsio Maria Rita, Stephan Celso, LourencËo Roberto Wagner, Andrade Valmir Roberto, de Brito Arduino Marylene, et al. Spatial Distribution of the Risk of Dengue and the Entomological Indicators in SumareÂ, State of SaÄo Paulo, Brazil . PLOS Neglected Tropical Diseases . 2014 ; 8 ( 5 ):e2873. https://doi.org/10.1371/journal.pntd.0002873 PMID: 24831806 .
28. Lin YP , Yang ZF , Liang Y , Li ZT , Helen S. Bond , Chua HY , et al. Population seroprevalence of antibody to influenza A(H7N9) virus, Guangzhou, China . BMC Infectious Disease . 2016 ; 16 : 632 . https://doi.org/ 10.1186/s12879-016-1983-3 PMID: 27814756 .
29. Feng L , Wu JT , Liu X , Yang P , Tsang TK , Jiang H , et al. Clinical severity of human infections with avian influenza A(H7N9) virus , China, 2013 /14. Euro Surveill. 2014 ; 19 ( 49 ):pii = 20984 . Available online: http://www.eurosurveillance.org/ViewArticle.aspx?ArticleId= 20984 .
30. Yu H , Cowling BJ , Feng L , Lau EH , Liao Q , Tsang TK , et al. Human infection with avian influenza A H7N9 virus: an assessment of clinical severity . Lancet . 2013 ; 382 ( 9887 ): 138 ± 45 .http://dx.doi.org/10. 1016/S0140- 6736 ( 13 ) 61207 - 6 .
31. Ling F , Chen EF , Liu Q , Miao Z , Gong Z. Hypothesis On The Source, Transmission and Characteristics of Infection of Avian Influenza A (H7N9) Virus±Based On Analysis of Field Epidemiological Investigation and Gene Sequence Analysis . Zoonoses and Public Health . 2015 ; 62 : 29 ± 37 . https://doi.org/10.1111/ zph.12110 PMID: 24628953 .
32. Wu P , Jiang H , Joseph T .Wu, Chen EF , He JF , Zhou H , et al. Poultry Market Closures and Human Infection with Influenza A(H7N9) Virus , China, 2013 ± 14 . Emerging Infectious Diseases. 2014 ; 20 ( 11 ): 1891 ±93. http://dx.doi.org/10.3201/eid2011.140556.PMID: 25340354 .
33. Zhou L , Ren RQ , Qu JM , Kang M , Wang XX , Fiona Havers , et al. Risk Factors for Influenza A(H7N9) Disease in China, a Matched Case Control Study, October 2014 to April 2015 . Open Forum Infect Dis. 2016 ; 3 ( 3 ):ofw182. https://doi.org/10.1093/ofid/ofw182 PMID: 27704029 .
34. Carrat Fabrice , Valleron Alain Jacques. Epidemiologic Mapping Using the ªKrigingº Method:Application to an InfluenzaÐlike Illness Epidemic in France [J] Am J Epidemiol , 1992 , 135 ( 11 ): 1293 ± 1300 .PMID: 1626546 .