Effects of closing and reopening live poultry markets on the epidemic of human infection with avian influenza A virus

Jian Lu, Wendong Liu, Rui Xia, Qigang Dai, Changjun Bao, Fenyang Tang, Yefei Zhu,3,?, Qiao Wang,?

1School of Information Science and Engineering, Southeast University, Nanjing, Jiangsu 211111, China;

2 Department of Acute Infectious Disease Control and Prevention, Jiangsu Provincial Center for Disease Control and Prevention, Nanjing, Jiangsu 210009, China;

3 Center of Laboratory Medicine, The Second Affiliated Hospital of Nanjing Medical University, Nanjing, Jiangsu 210011, China.

Effects of closing and reopening live poultry markets on the epidemic of human infection with avian influenza A virus

Jian Lu1,△, Wendong Liu2,△, Rui Xia1, Qigang Dai2, Changjun Bao2, Fenyang Tang2, Yefei Zhu2,3,?, Qiao Wang1,?

1School of Information Science and Engineering, Southeast University, Nanjing, Jiangsu 211111, China;

2 Department of Acute Infectious Disease Control and Prevention, Jiangsu Provincial Center for Disease Control and Prevention, Nanjing, Jiangsu 210009, China;

3 Center of Laboratory Medicine, The Second Affiliated Hospital of Nanjing Medical University, Nanjing, Jiangsu 210011, China.

Live poultry markets (LPMs) are crucial places for human infection of influenza A (H7N9 virus). In Yangtze River Delta, LPMs were closed after the outbreak of human infection with avian influenza A (H7N9) virus, and then reopened when no case was found. Our purpose was to quantify the effect of LPMs’ operations in this region on the transmission of influenza A (H7N9) virus. We obtained information about dates of symptom onset and locations for all human influenza A (H7N9) cases reported from Shanghai, Jiangsu and Zhejiang provinces by May 31, 2014, and acquired dates of closures and reopening of LPMs from official media. A two-phase Bayesian model was fitted by Markov Chain Monte Carlo methods to process the spatial and temporal influence of human cases. A total of 235 cases of influenza A (H7N9) were confirmed in Shanghai, Jiangsu and Zhejiang by May 31, 2014. Using these data, our analysis showed that, after LPM closures, the influenza A (H7N9) outbreak disappeared within two weeks in Shanghai, one week in Jiangsu, and one week in Zhejiang, respectively. Local authorities reopened LPMs when there was no outbreak of influenza A (H7N9), which did not lead to reemergence of human influenza A (H7N9). LPM closures were effective in controlling the H7N9 outbreak. Reopening of LPM in summer did not increase the risk of human infection with H7N9. Our findings showed that LPMs should be closed immediately in areas where the H7N9 virus is confirmed in LPM. When there is no outbreak of H7N9 virus, LPMs can be reopened to satisfy the Chinese traditional culture of buying live poultry. In the long term, local authorities should take a cautious attitude in permanent LPM closure.

live poultry market, avian influenza A (H7N9)

Introduction

In March 2013, cases of human infection with a novel avian influenza A (H7N9) virus emerged in China[1-7]. Until May 31, 2014, two waves of human influenza A (H7N9) epidemic affected China, and caused 399 cases. The confirmed infection cases were concentrated in eastern China, including the Yangtze River Delta region and Beijing. Most H7N9 virus-infected patients developed severe pneumonia and acute respiratorydistress syndrome[8]with an approximate mortality of 30%[9-10]. Chinese health authorities reported that H7N9 virus was found in domestic poultries, such as chickens, pigeons, and ducks[11-12]. Since there is no evidence of human-to-human transmission so far[13-14], the main focus of intervention measure is to reduce transmission from poultry to humans[15-17].

Live poultry was known as a potential source of human infection of H7N9 virus as 75% of H7N9 virus-infected patients reported recent exposures to poultry[6]. The most likely locations for poultry-to- human transmission of H7N9 were considered as live poultry markets (LPMs)[4-6,18-20]. In China, LPMs are common places. A survey in Guangzhou in 2006 showed that 80% of households reported visit to LPMs at least once a year[21]. Even in 2013, 47% of respondents in Guangzhou reported visiting LPMs more than one time per year[22]. Studies in Hong Kong indicated that LPMs were high risk regions of human infection with H5N1 virus[23-24]. Recently, a research on 8,943 LPMs in China has shown a strong positive association between market-level H7N9 presence and local density of LPMs[25].

In April 2013, closures of LPMs in four cities were tried as an effective measure in curtailing the spread of H7N9 virus[20]. However, closures of LPMs were only temporary, and LPMs were reopened in June 2013. The influence of reopening LPMs on H7N9 virus transmission is still lack of analysis. Since H7N9 virus was initially detected in the Yangtze River Delta, we aimed to assess the effects of closing and reopening LPMs around this area. In this study, we proposed a two-phase Bayesian network, which may provide scientific evidence for authorities to employ more comprehensive measures to LPMs in those regions affected by H7N9 virus.

Materials and methods

Data source

Shanghai and two adjacent provinces, Jiangsu and Zhejiang, were selected as the region for the current study. As of 31 May 2014, a total of 235 laboratory-confirmed human cases of H7N9 had occurred in the region. For each case, we collected information about the dates of symptom onset and locations from H7N9 epidemic report released by the National Health and Family Planning Commission of the People's Republic of China (NHFPC). By the same way, we obtained the dates of closing and reopening of LPMs in the region.

Spatial and temporal analysis

To assess the effect of LPMs' operations on the epidemic of H7N9, we proposed a two-phase Bayesian model. We segmented the human cases of each province into time series. For convenience, we treated the first week in 2013 as the first week in our study. Let Xi,tdenote the number of cases at the tthweek in province i, and Zi,tdenote the phase of epidemic in city or province i, where Zi,t= 0,1 corresponds to the non-outbreak phase and outbreak phase, respectively. Index i = 1, 2, 3 denotes Shanghai, Jiangsu, and Zhejiang, respectively.

To measure the effect of LPM closure and reopening in each specified city or province, we firstly assumed that the weekly numbers of human cases obey the conditional Poisson distribution under the given phase Z for each specified city or province. Then, we introduced the spatial and temporal dependence on the phases of all provinces.

For Shanghai, we supposed that the observation X1,tobeys a Poisson distribution under the given phase Z1,t. Given Z1,t= 0, which denotes that there is no outbreak of influenza A (H7N9):

where λ1,tdenotes the mean number of H7N9 cases when there is no outbreak of H7N9, where the mean number refers to the mean number of H7N9 cases every week. Since Shanghai closed LPMs in the 14thand 53thweek, and reopened in the 25thand 70thweek, respectively, we segmented λ1,tinto five different time periods as below:

When Z1,t= 1, which means the presence of outbreak of H7N9, we assumed:

where ρ1,tdenotes the infection force of H7N9 transmission. Similarly, we segmented ρ1,tlike λ1,t:

In the same way as (1)-(4), we constructed the conditional Poisson distribution of X2,tand X3,tfor Jiangsu and Zhejiang, respectively.

Fig. 1 Weekly distribution of human influenza A (H7N9) cases. The weekly number of H7N9 in Shanghai, Jiangsu, and Zhejiang are plotted, where we treat the first week in 2013 as the 1stweek of our model.

To investigate the patterns of H7N9 outbreak, we adopted the spatial and temporal dependence into the model, where the phase Z1,tfor Shanghai at week t was not only dependent on its previous phase Z1,t-1, but also on its adjacent provinces. This pattern can be described by a generalized linear model[26]. For province i at week t, the probability of Zi,tis defined as follows:

where Θ and Φ were parameter matrices that denoted the temporal and spatial influences. Since phase Z takes only two possible values, the Θ is a 2 × 2 matrix. Θm,nmeasures the possibility that phase Zi,ttakes value n given the previous phase Zi,t-1equals to m. The matrix Φ is also a 2 × 2 matrix, in which Φl,kdenotes the possibility that phase Zi,ttakes value k given that the phase of the adjacent location Zj,tequals to l. Each element in Θ and Φ obeys a zero-mean Gaussian prior distribution.

The sequence of state Z for a specific city reflects the local effect of LPM closure or reopening. This Bayesian model is similar to models previously used to detect the influenza from Twitter[26]. We compute all the parameters by Markov Chain Monte Carlo methods in R software (Version 3.2.0).

Results

As of 31 May 2014, there were 235 laboratoryconfirmed human cases of H7N9 reported in the study areas, 41 cases in Shanghai, 55 cases in Jiangsu, and 139 cases in Zhejiang, respectively (Fig. 1). The first epidemic wave occurred between19 February 2013 and 5 April 2013 and led to 107 cases, and the second wave caused 128 cases between 7 October 2013 and 22 May 2014.

For prevention of the first epidemic wave in 2013, local authorities closed LPMs in Shanghai on April 6, in Jiangsu on April 16, and in Zhejiang on April 22, respectively. When this wave ended, LPMs reopened in Shanghai on June 20, Jiangsu on June 1, and Zhejiang on May 27, respectively. For the second epidemic wave in 2014, the three provinces have taken different strategies on LPMs based on their epidemic situations. Shanghai closed and reopened LPMs on January 3 and April 30 in 2014, respectively. Jiangsu did not close LPMs during the second epidemic wave. Zhejiang closed all LPMs on February 8 in 2014, which never reopened. The details of dates of closing and opening of LPMs are shown in Table 1.

The analysis of phase and estimated number of H7N9 cases

Based on the two-phase Bayesian model, we estimated the phases for each province. The Kruskal-Wallis test with chi-squared = 88.728, df = 1, P-value＜2.2e-16 showed a significant difference between the outbreak-phase and the non-outbreak phase.

The number of cases and phases of H7N9 epidemic for Shanghai are shown in Fig. 2. Outbreak phases in Shanghai only occurred among the 13thand 15thweek.The date of first closure of LPMs in Shanghai was the 14thweek. For the first closure of LPMs, the number of cases decreased from 8.5 at the 14thweek to 6.2 at the 15thweek. From the 16thweek to 25thweek, the number of H7N9 cases was estimated at 0.07. The estimated mean number of H7N9 cases before and after the firstreopening of LPMs was 0.07 and 0.14, respectively. The estimated mean number of H7N9 cases before and after the second reopening of LPMs was 0.18 and 0.19, respectively. The box-plot analysis and details of the estimated mean number of H7N9 cases are shown in Fig. 3 and Table 2, respectively.

Table 1 The dates of closing and reopening of live poultry markets (LPMs) in Shanghai, Jiangsu and Zhejiang provinces.

Fig. 2 Number of human influenza A (H7N9) cases and phases of epidemic in Shanghai. The blue line denotes the number of human cases, and the red circles and green circles denotes the outbreak phase and non-outbreak phase of Shanghai, respectively. In our model, Shanghai only occurred one wave of H7N9 epidemic from the 13thweek to the 15thweek.

Fig. 3 The box-plot of mean number of human influenza A (H7N9) cases in Shanghai. The mean numbers of H7N9 cases before and after live poultry markets were open or closed in Shanghai were analyzed.

Table 2 Estimated mean number of human influenza A (H7N9) cases in different time period segmented by the closures and reopening LPMs of Shanghai, Jiangsu and Zhejiang.

Fig. 4 Number of human influenza A (H7N9) cases and phases of epidemic in Jiangsu. The blue line denotes the number of human cases, and the red circles and green circles denotes the outbreak phase and no-outbreak phase of Jiangsu, respectively. In our model, Jiangsu only occurred one wave of H7N9 epidemic from the 12thweek to the 15thweek.

The numbers of cases and phases of H7N9 epidemic for Jiangsu are shown in Fig. 4, where the outbreak phases only occurred between the 12thand 15thweek. The da e of he firs closure of LPMs in Jiangsu was the 15thweek. The number of H7N9 cases before and after the first closure was 6.5 and 0.57, respectively. The estimated mean number of H7N9 cases before and after the first reopening was 0.57 and 0.54, respectively. The box-plot analysis and details of the estimated mean number of H7N9 cases are shown in Fig. 5 and Table 2, respectively.

The number of cases and phases of H7N9 cases for Zhejiang are shown in Fig. 6. The outbreak phases in Zhejiang occurred among the 14th-16thweek, and the 53th-58thweek. The dates of the first and second closure of LPMs were the 16thand 58thweek, respectively. For the first closure, the estimated mean number of H7N9 cases decreased from 9.9 to 0.14. The estimated mean number of H7N9 cases before and after the second closure was 10.1 and 0.25, respectively. The estimated mean number of H7N9 cases before and after the first reopening was 0.14 and 0.43, respectively. The boxplot analysis and details of the estimated mean number of H7N9 cases are shown in Fig. 7 and Table 2, respectively.

Relationships between temperature and H7N9 cases

The relationships between the number of H7N9 cases and temperature are shown in Fig. 8, which plots the daily lowest and highest temperatures as well as thenumber of cases in Shanghai, Jiangsu, and Zhejiang between February 2013 and May 2014. The highest temperature ranged from 3°C to 29°C, and the lowest temperature ranged from -4°C to 15°C.

Fig. 5 The box-plot of mean number of human influenza A (H7N9) cases in Jiangsu. The mean numbers of H7N9 cases before and after LPMs were closed and opened in Jiangsu were analyzed.

Spatial and temporal dependence

Our model analyzed the spatial and temporal dependence for the cases of H7N9 in Shanghai, Jiangsu, and Zhejiang provinces. We estimated the temporal dependence of phase Z. When there was no outbreak of H7N9 at the (t-1)thweek, the probability of no outbreak of H7N9 at the tthweek was 0.97 (95% CI: 0.9671-0.9686). When there happened an outbreak of H7N9 at the (t-1)thweek, the probability of the outbreak of H7N9 at the tthweek was 0.72 (95% CI: 0.7143-0.7279). The spatial dependence of phase Z was also considered. The province had a high probability of 0.96 (95% CI: 0.9607-0.9617) of the no-outbreak phase when there was no outbreak of H7N9 at adjacent provinces. The probability of outbreak phase was 0.21 (95% CI: 0.2049-0.2099) when adjacent provinces had outbreak of H7N9.

Discussion

Our model systematically evaluated the effects of intervention measures to the LPMs on the epidemic of H7N9 in the study areas. According to our findings, Shanghai and Jiangsu had one wave of H7N9 epidemic, while Zhejiang had two waves of epidemic. The H7N9 epidemic in Shanghai started at the 13thweek, and LPMs were shut down on the next week. Then, the H7N9 epidemic faded after the 15thweek. Our model showed that the closure of LPMsreduced the mean number of H7N9 cases from 8.5277 at the 13th-14thweek to 0.07 at the 16th-25thweek. Furthermore, the effectiveness of closures of LPMs in Jiangsu and Zhejiang can also be found by our model in a similar way. This result supported the idea that LPMs have a major role in poultry-to-human transmission of H7N9 virus[20,27]. From the public health standpoint, once signs of outbreak of human infection with H7N9 virus are detected by routine monitoring, especially in winter and spring, local authorities should decisively close LPMs as soon as possible to block the epidemic.

Fig. 6 Number of human influenza A (H7N9) cases and phases of epidemic in Zhejiang. The blue line denotes the number of human cases, and the red circles and green circles denotes the outbreak phase and no-outbreak phase of Zhejiang, respectively. In our model, Zhejiang had two waves of H7N9 epidemic. The first wave happened from the 14thweek to 16thweek, and the second wave from the 53thweek to the 58thweek.

Conversely, we found that reopening LPMs in the non-outbreak phase did not increase the risk of human infections with H7N9 virus. As mentioned before, Shanghai reopened the LPMs twice, while Jiangsu and Zhejiang reopened LPMs once. The date of the first reopening of LPMs in Shanghai was June 20, 2013. In our model, Shanghai remained in the non-outbreak phases of H7N9 epidemic after reopening. The mean number of H7N9 cases before and after the reopening was estimated at 0.07 and 0.14, respectively. Similarly, the second reopening of LPMs in Shanghai, and reopening of LPMs in Jiangsu and Zhejiang did not lead to outbreak of H7N9 epidemic. Climate conditions in the hot summer may not be suitable for H7N9 virus to survive. Many experimental studies on H5N1 have revealed that temperature can affect the survival time of the H5N1 virus[28-29]. It is also reported that H7N9 virusis inactivated in high temperature[30]. Analysis of H7N9 cases and their corresponding temperature showed that there were no cases when the lowest daily temperature were less than 15°C or the highest daily temperature exceeded over 29°C. In this case, closing LPMs would not be a wise decision to the local authorities, taking into account the heavy cost.

Fig. 7 The box-plot of mean number of human influenza A (H7N9) cases in Zhejiang. The mean numbers of H7N9 cases before and after the LPMs’ operations in Zhejiang were analyzed.

The spatial and temporal dependence of the outbreak of H7N9 epidemic was discussed in our model. Analysis of temporal dependence showed that, in the same province, the phases of the adjacent weeks showed strong correlations. That is, if the outbreak phase of H7N9 epidemic was found at the present, the probability of occurring of H7N9 at the next week was as high as 0.97. This result suggested that, when the H7N9 epidemic was in the outbreak phase, LPMs should be rapidly closed to protect against poultry-to-human transmission. Spatial analysis revealed that, in the same week, the phases of the adjacent provinces showed mild correlations. This result was from the analysis of all cases between 2013 and 2014. In 2013, the first wave of H7N9 occurred almost on the same week in Shanghai, Jiangsu and Zhejiang. This phenomenon in 2013 reflected a strong spatial correlation on the epidemic outbreak among adjacent provinces. But after effective precautionary measures were taken in 2014, only Zhejiang detected the outbreak of H7N9 cases in the 53rd-58thweek. No outbreak of H7N9 occurred in Shanghai and Jiangsu. This phenomenon in 2014 showed that the spatial correlation in epidemic outbreak among adjacent provinces was weak. Combining all outbreak phases in 2013 and 2014, the spatial analysis of our model showed a mild spatial correlation.

Fig. 8 The relationship between the number of human influenza A (H7N9) cases and the daily lowest temperature in (A), and daily highest temperature in (B).

Currently, the considerations of permanent LPM closure and central slaughtering for the poultry supply chain are fiercely discussed. Although the evidence has shown that LPM closures can dramatically cut down H7N9 cases, the policymakers should also consider the social and economic impacts of LPM closures. At the end of June 2013, the first wave of epidemic caused 60 billion RMB loss in China's poultry industry, and more than 50 million enterprises and farmers suffered losses[10]. Based on our findings, we propose that local authorities should be more reasonable in using the strategy of closing LPMs.

Limitations of our study have to be addressed here. Firstly, we made a simple assumption that LPMs of each province were closed on the same week. In fact, even in the same province, the dates of LPM closures in different cities may have one week's bias. In Jiangsu, the provincial government started Level 4 Emergency Response to avian influenza on April 3, 2013[31]. Based on this response, LPMs were shut down on April 8 (the 14thweek) in Nanjing and Suzhou, on April 9 (the 14thweek) in Zhenjiang, on April 16 (the 15thweek) in Yancheng and Xuzhou. Our model assumed 15thweek as the date of LPM closures for Jiangsu. This assumption may have certain impact on effective analysis of LPM closures. To deal with this situation, a more complex model needs to be considered. Secondly, we targeted our surveys to three provinces in the Yangtze River delta region, where most H7N9 cases occurred. We were unable to examine the effect of LPM closure and reopening in other Chinese provinces which had small numbers of H7N9 cases.

Despite these limitations, it is noteworthy that we have analyzed the effectiveness of closing and reopening of LPMs by using the H7N9 epidemic phases which were built through a two-phase Bayesian model. The result can provide a basis for local authorities for taking proper control strategies on LPMs.

References

[1] WHO/GAR (2013). Human infection with influenza A (H7N9) virus in China[R/OL]. Available: http://www.who. int/csr/don/2013_04_01/en. Accessed 1 April 2013.

[2] Yang F, Wang J, Jiang L, et al. A fatal case caused by novel H7N9 avian influenza A virus in China[J/OL]. Emerg Microbes Infect, 2013,2(4):el9.

[3] Gao R, Cao B, Hu Y, et al. Human infection with a novel avian-origin influenza A (H7N9) virus[J]. N Engl J Med, 2013,368(20):1888–1897.

[4] Li Q, Zhou L, Zhou M, et al. Epidemiology of human infections with avian influenza A (H7N9) virus in China[J]. N Engl J Med, 2014,370:520–532.

[5] Chen Y, Liang W, Yang S, et al. Human infections with the emerging avian influenza A H7N9 virus from wet market poultry: clinical analysis and characterisation of viral genome[J]. Lancet, 2013,381(9881):1916–1925.

[6] Cowling BJ, Jin L, Lau EHY, 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[J]. Lancet, 2013, 382(9887):129–137.

[7] Zhuang Q, Wang S, Wu M, et al. Epidemiological and risk analysis of the H7N9 subtype influenza outbreak in China at its early stage[J]. Chin Sci Bull, 2013,58:1403–1406.

[8] Gao H, Lu H, Cao B, Du B, Shang H, et al. Clinical findings in 111 cases of influenza A (H7N9) virus infection[J]. N Engl J Med, 2013,368(24):2277–2285.

[9] Li Q, Zhou L, Zhou M, Chen Z, et al. Preliminary Report: Epidemiology of the Avian Influenza A (H7N9) Outbreak in China[J]. N Engl J Med, 2013,370:520–532.

[10] China’s Ministry of Agriculture (2014). National plan for Eradication of avian influenza A (H7N9)[R/OL]. Available: http://www.chinairn.com/news/20140221/13542832.html.

[11] Kahn RE, Richt JA. The Novel H7N9 Influenza A Virus: Its Present Impact and Indeterminate Future[J]. Vector-Borne Zoonot, 2013,13(6):347–348.

[12] Shi J, Deng G, Liu P, et al. Isolation and characterization of H7N9 viruses from live poultry markets-Implication of the source of current H7N9 infection in humans[J]. Chin Sci Bull, 2013,58:1857–1863.

[13] Li Q, Zhou L, Zhou M, et al. Preliminary report: Epidemiology of the avian influenza A (H7N9) outbreak in China[J]. N Engl J Med, 2014,370:520–532.

[14] Zhu H, Wang D, Kelvin DJ, et al. Infectivity, transmission, and pathology of human-isolated H7N9 influenza virus in ferrets and pigs[J]. Science, 2013,341:183–186.

[15] Horby P. H7N9 is a virus worth worrying about[J]. Nature, 2013,496(7446):399.

[16] Uyeki TM, Cox NJ. Global concerns regarding novel influenza A (H7N9) virus infections[J]. N Engl J Med, 2013,368(20):1862–1864.

[17] Han J, Jin M, Zhang P, et al. Epidemiological link between exposure to poultry and all influenza A (H7N9) confirmed cases in Huzhou city, China, March to May 2013[J]. Euro Surveill, 2013,18(20):20481.

[18] He Y, Liu P, Tang S, et al. Live poultry market closure and control of avian influenza A (H7N9), Shanghai, China[J]. Emerg Infect Dis, 2014,20(9):1565–1566.

[19] Bao C, Cui L, Zhou M, et al. Live-animal markets and influenza A (H7N9) virus infection[J]. N Engl J Med, 2013,368(24):2337–2339.

[20] Yu H, Wu J, Cowling BJ, et al. Effect of closure of live poultry markets on poultry-to-person transmission of avian influenza A H7N9 virus: an ecological study[J]. Lancet, 2014,383(9916):541–548.

[21] Liao Q, Lam WT, Leung GM, et al. Live poultry exposure, Guangzhou, China, 2006[J]. Epidemics, 2009, 1(4):207–212.

[22] Wang L, Cowling BJ, Wu P, et al. Human exposure to live poultry and psychological and behavioral responses to influenza A(H7N9), China[J]. Emerg Infect Dis, 2014,20(8):1296–1305.

[23] Shortridge KF, Gao P, Guan Y, et al. Interspecies transmission of influenza viruses: H5N1 virus and a Hong Kong SAR perspective[J]. Vet Microbiol, 2000,74(1–2):141–147.

[24] Mounts AW, Kwong H, Izurieta HS, et al. Case-control study of risk factors for avian influenza A(H5N1) disease, Hong Kong, 1997[J]. J Infect Dis, 1999,180(2):505–508.

[25] Gilbert M, Golding N, Zhou H, et al. Predicting the risk of avian influenza A H7N9 infection in live-poultry markets across Asia[J]. Nat Commun, 2014,5:4116.

[26] Li J, Cardie C. Early stage influenza detection from twitter. arXiv:1309.7340[cs.SI]. Available: http://arxiv.org/ abs/1309.7340.

[27] Gao G. Influenza and the Live Poultry Trade[J]. Science, 2014,344(6181):235.

[28] Paek MR, Lee Y J, Yoon H, et al. Survival rate of H5N1 highly pathogenic avian influenza viruses at different temperatures[J]. Poult Sci, 2010,89(8):1647–1650.

[29] Tiensin T, Nielen M, Songserm T, et al. Geographic and temporal distribution of highly pathogenic avian influenza A virus (H5N1) in Thailand, 2004–2005: an overview[J]. Avian Dis, 2007,51(1 Suppl):182–188.

[30] Zhang Z, Xia Y, Lu Y, et al. Prediction of H7N9 epidemic in China[J]. Chin Med J, 2014,127(2):254–260.

[31] Jiangsu Province (2013). Jiangsu starts the 4th level Emergency Response to avian influenza[EB/OL]. Available: http://www.chinanews.com/gn/2013/04-03/4703079.shtml.

△These authors contributed equally to this work.

? Yefei Zhu, Department of Acute Infectious Disease Control and Prevention, Jiangsu Provincial Center for Disease Control and Prevention, Nanjing, China. Email: jszyf@jscdc.cn, Tel: 02558509923, Fax: 02558509923; Qiao Wang, School of Information Science and Engineering, Southeast University, Nanjing, Jiangsu China.

Email: qiaowang@seu.edu.cn, Tel: 02552091680, Fax: 02552091680. Received 03 April 2015, Revised 03 July 2015, Accepted 26 October 2015, Epub 29 December 2015

R513, Document code: A

The authors reported no conflict of interests.