We’re sorry, something doesn't seem to be working properly.
Please try refreshing the page. If that doesn't work, please contact support so we can address the problem.
Spatio-temporal pattern analysis for evaluation of the spread of human infections with avian influenza A(H7N9) virus in China, 2013–2014
BMC Infectious Diseasesvolume 17, Article number: 704 (2017)
A large number (n = 460) of A(H7N9) human infections have been reported in China from March 2013 through December 2014, and H7N9 outbreaks in humans became an emerging issue for China health, which have caused numerous disease outbreaks in domestic poultry and wild bird populations, and threatened human health severely. The aims of this study were to investigate the directional trend of the epidemic and to identify the significant presence of spatial-temporal clustering of influenza A(H7N9) human cases between March 2013 and December 2014.
Three distinct epidemic phases of A(H7N9) human infections were identified in this study. In each phase, standard deviational ellipse analysis was conducted to examine the directional trend of disease spreading, and retrospective space-time permutation scan statistic was then used to identify the spatio-temporal cluster patterns of H7N9 outbreaks in humans.
The ever-changing location and the increasing size of the three identified standard deviational ellipses showed that the epidemic moved from east to southeast coast, and hence to some central regions, with a future epidemiological trend of continue dispersing to more central regions of China, and a few new human cases might also appear in parts of the western China. Furthermore, A(H7N9) human infections were clustering in space and time in the first two phases with five significant spatio-temporal clusters (p < 0.05), but there was no significant cluster identified in phase III.
There was a new epidemiologic pattern that the decrease in significant spatio-temporal cluster of A(H7N9) human infections was accompanied with an obvious spatial expansion of the outbreaks during the study period, and identification of the spatio-temporal patterns of the epidemic can provide valuable insights for better understanding the spreading dynamics of the disease in China.
Since the first avian influenza A (H7N9) outbreak in human occurred in March 2013 in Zhejiang, H7N9 outbreaks became an emerging issue for China health: they have caused numerous disease outbreaks in domestic poultry and wild bird populations, and threatened human health severely [1,2,3]. Between March 2013 and December 2014, a total of 460 A(H7N9) human infections with severe respiratory illness and mortality have been reported to the World Health Organization .
Phylogenetic analysis of H7N9 demonstrated that the novel avian-origin influenza A(H7N9) virus was a novel triple reassortant of avian influenza A H7N3, A H7N9, and A H9N2 viruses [5,6,7]. The virus was previously found only in wild bird populations, but in this outbreak it has been proved to be a novel strain, the first of its kind ever found to infect humans . Gene sequence studies have showed that the virus combines with human receptors more easily than with birds . The virus can infect birds without symptoms, but it can cause severe respiratory illness and result in high mortality in human populations. In spite of the long presence of A(H7N9) human infections in China, the long-term spatiotemporal dynamics and geographic spread mechanisms of H7N9 outbreaks in China have not been well characterised. So it is necessary to investigate the transmission pattern of A(H7N9) human infections and identify the significant clustering of influenza A(H7N9) human cases. Understanding how A(H7N9) human infections spread in China can help controlling the epidemic in the future.
Standard deviational ellipse analysis has been used in some studies to investigate the directional trend and model the spread of outbreak over time [10,11,12,13,14], but the research on the directional trend of H7N9 is very limited, and so far there was only one former study used this method to explore the directional distribution of influenza A(H7N9) human cases in the Yangtze River Delta between 19 February 2013 and 18 April 2013 . In addition, as an helpful adjunct for assessing the spatiotemporal distribution of an epidemic, the SatScan method has been undertaken in some previous studies to explore the spatio-temporal clustering of H7N9 outbreaks in China [3, 16]. Although some spatio-temporal clusters of H7N9 have been identified in these researches, the study periods are not long enough neither (the longest study in these papers was conducted between 19 February 2013 and 31 March 2014 ). However, the spreading of influenza A(H7N9) human infections is a dynamic development process and the spread mechanisms of the epidemic are still unclear till now, consequently these researches may not be adequate for seasonal H7N9 outbreaks in different epidemic phases, and further studies with more case data covering all seasons of a full year are still necessary to investigate the periodical spatiotemporal spreading pattern of H7N9.
In order to examine the ever-changing spreading pattern of H7N9 outbreaks, by using information on all influenza A(H7N9) human cases reported in China between March 2013 and December 2014, GIS-based standard deviational ellipse analyses and SaTScan-based retrospective space-time permutation scan statistics were conducted in different epidemic phases of China. The main purposes of this study were: (1) to investigate the epidemic characteristics and the directional trend of the H7N9 outbreaks, which can provide clues for understanding the spreading dynamics of the epidemic that have occurred in China in different phases. And (2) to identify the significant presence of spatial-temporal clustering of influenza A(H7N9) human cases during each phase, and this is helpful in confirming geographical locations and times of the year in which H7N9 outbreaks are most likely to occur, and then the government can apply control strategies specifically. As a complement to epidemiological studies, results of this study can provide valuable insights for better understanding the spreading dynamics of the disease, and hence improve the efficiency of future strategic planning in the prevention and control of H7N9 outbreaks.
In order to fully examine the spread of the H7N9 outbreaks in humans, the disease data consisted of 460 influenza A(H7N9) human cases reported in China from March 2013 through December 2014, which were aggregated at the municipal level, and the spatio-temporal attributes of each H7N9 human case were used, including the location (latitude, longitude) and the date of influenza A(H7N9) human cases.
We obtained information about the confirmed influenza A(H7N9) human cases from the World Health Organization website . A Chinese topographic map (1:4,000,000 scales) available from the National Fundamental Geographic Information System of China  was used to visualize the reported influenza A(H7N9) human cases, the identified standard deviational ellipses and the spatio-temporal clusters of H7N9 outbreaks by using ArcGIS 10.1 software (ESRI Inc., Redlands, CA, USA).
Based on the daily and monthly data of influenza A(H7N9) human cases reported in China from March 2013 to December 2014, the corresponding epidemic curves of H7N9 were created by R 3.1.3 software  to display the magnitude and temporal trend of H7N9 outbreak. According to the epidemic curves, we regarded a full time period as an epidemic phase, in which the number of influenza A(H7N9) human cases peaked and then decreased until the cases finally disappeared. Accordingly, a quantitative description of H7N9 transmission could be presented by displaying the staged dynamics of the standard deviational ellipse and the spatio-temporal clusters of the epidemic.
Directional trend analysis
The standard deviational ellipse analysis can well measure the spatial distribution of geographic features around their geometric center, and provide useful information about feature dispersion and orientation . How standard deviational ellipse works has been described in detail in ArcGIS Help 10.1 . In this study, standard deviational ellipse analysis was carried out with one standard deviation in ArcGIS 10.1 (ESRI Inc., Redlands, CA, USA) to investigate the size and the movement of the H7N9 epidemic, and then to measure whether the distribution of influenza A(H7N9) human cases displays a directional trend during the three epidemic phases.
This analysis can create ellipse to summarize the spatial characteristics of the influenza A(H7N9) human cases: central tendency, dispersion, and directional trends. The ellipse shows the spatial spread of a set of point locations, and helps us to determine whether the distribution of these points is elongated and hence they spread with a particular orientation.
Spatio-temporal cluster analysis
For identifying the spatio-temporal clusters of A(H7N9) human infections in each epidemic phase, retrospective space-time permutation scan statistics were applied by using SaTScan 9.4.2 . The space-time permutation model only requires case data with exact spatial location and time, but no need about the background population at risk [22,23,24].
The scanning window in the space-time permutation model was a cylinder with the base representing space and the height representing time, respectively . In this research the window moves both in time and in space to cover each possible time interval for each possible geographic location, and the model parameters for maximum spatial and temporal window sizes were set such that a cluster could include a maximum of 50% of all influenza A(H7N9) human cases. During the analysis, spatio-temporal clusters occur when an excess number of H7N9 outbreaks are observed, and the log-likelihood ratio statistic [26,27,28] was used to evaluate whether the cylinder contained a cluster or not. The scanning window with the maximum log-likelihood ratio was defined be the most likely cluster, and the other windows with statistically significant log-likelihood ratio were considered as secondary clusters, which were ranked by the value of their log-likelihood ratio . The P-value for the detected clusters was evaluated by using Monte Carlo simulation . Ensuring adequate precision, the number of Monte Carlo simulations was set at 999, and the time units of 5 days were used to reflect potential disease spread. For all analyses, the most likely cluster and non-overlapping secondary clusters with statistical significance of P < 0.05 were reported.
Spacial distribution analysis
During the H7N9 epidemic period (March 2013 to December 2014, 658 days) studied in this paper, a total of 460 influenza A(H7N9) human cases occurred in 18 provinces, municipalities and autonomous regions of China. The spatial distribution of influenza A(H7N9) human cases is shown in Fig. 1. It clearly indicates that the distribution of A(H7N9) human infections was heterogeneous at the provincial level. The results show that the largest number of influenza A(H7N9) human cases was reported from Zhejiang (142), followed by Guangdong (104), Jiangsu (62) and Shanghai (42), all of them were located in the southeast coast of China, which accounted for about 76.1% of influenza A(H7N9) human cases (350/460).
The epidemic curve shows the temporal pattern of H7N9 outbreaks in humans from March 2013 to December 2014 in China, which consisted of three distinct phases (Fig. 2): phase I was from March 13, 2013 to May 31, 2013 (n = 130), with peak observed in April 16, 2013; phase II was from June 1, 2013 to May 31, 2014 (n = 294), with peak observed in January 20, 2014; phase III was from June 1, 2014 to December 31, 2014 (n = 36), with peak observed in December 21, 2014. From Fig. 2 it also can be seen that there was a similar seasonal pattern showed by different time scales of day and month: influenza A(H7N9) season generally began in early winter (November to December), and then the epidemic reached its peak between late winter and late spring (January to April), but declined in incidence between May and June and then entered a sustained low incidence period (July to October).
Directional trend analysis
Table 1 summarizes the characteristics of standard deviational ellipses: X coordinate (CenterX) and Y coordinate (CenterY) for the mean center; the standard distance of X-axis (XStdDist) and of Y-axis (YStdDist), in which the larger one is long axis and the smaller one is short axis; the orientation of the ellipse (Rotation). Fig. 3 illustrates clearly the mean centers and the directional trends of influenza A(H7N9) human cases in the three phases by standard deviational ellipses (1 standard deviation).
In phase I, the mean center of the ellipse (x = 1,383,592.750160 m, y = 3,949,191.040550 m) was located in Xuancheng, Anhui province, and the ellipse was confined in eastern China, mainly including Jiangsu, Anhui, Zhejiang and Shanghai. Moreover, the ellipse showed a northwest-southeast elongation (long axis = 325.181481980 km, short axis = 277.821569251 km), with an angular rotation of 152.621015°.
In phase II, the spatial distribution of A(H7N9) human infections spread to southeast coastal areas, and the mean center of the ellipse (x = 1,184,775.846750 m, y = 3,516,687.088910 m) moved southwestward and finally located in Fuzhou, Jiangxi province, which was about 476.01 km southwest of the mean center of phase I. Consequently, the coverage of the ellipse in phase II also spread from the east to southeast coastal areas of China, mainly including Jiangsu, Anhui, Zhejiang, Shanghai, Jiangxi, Fujian, Guangdong, Hunan, and so on. In addition, the ellipse presented a southwest-northeast elongation (long axis = 785.662350587 km, short axis = 294.358161767 km), and the angular rotation was 26.947956°.
The H7N9 virus then spread rapidly from southeast coastal areas to central China, even dispersed to a slice of areas of western China during the third phase. In phase III, the mean center of the ellipse (x = 656,953.359159 m, y = 4,089,979.224960 m) moved northwestward and finally located in Nanyang, Henan province, which was about 779.27 km northwest of the mean center of phase II. Accordingly, the extent of the ellipse in phase III reached the maximum of the three phases, mainly including Taiwan, Fujian, Jiangxi, Zhejiang, Jiangsu, Anhui, Hubei, Henan, Shaanxi, Shanxi, Ningxia, Gansu, Hunan, and so on. Furthermore, the ellipse showed a northwest-southeast elongation (long axis = 1976.727605830 km, short axis = 539.929271709 km), with an angular rotation of 123.918912°.
Spatio-temporal cluster analysis
Using the surveillance data of influenza A(H7N9) human cases collected in the three epidemic phases, spatio-temporal cluster analysis identified two most-likely clusters and three secondary clusters during the study period.
Spatio-temporal cluster in phase I
There were totally 31 localities and 130 influenza A(H7N9) human cases contained in the first epidemic phase. The identified statistically significant spatio-temporal clusters of A(H7N9) human infections in this period was displayed in Fig. 4, which overlaid with the correlated influenza A(H7N9) human cases in phase I. The scan test identified two significant spatio-temporal clusters of H7N9 outbreaks (P < 0.05) between March 13, 2013 and May 31, 2013 (79 days), a most likely cluster and a secondary cluster (Fig. 4 and Table 2), subsequently the null hypothesis could be rejected. The most likely cluster (P = 0.0000055) including eight localities was observed in the period from April 27, 2013 to May 11, 2013 with a radius of 596.94 km, which was detected in Fujian, Jiangxi, Taiwan, Zhejiang and Hunan provinces, southeast of China. Based on the time period and the area of the cluster, this cluster showed 10 influenza A(H7N9) human cases, which accounted for 7.70% of all the human cases of phase I (10/130), when only 1.38 human cases could be expected theoretically. Other, the secondary cluster 2 (P = 0.044) including 12 localities was located in Jiangsu, Anhui provinces and Shanghai municipality, Yangtze river delta region of eastern China, with a radius of 271.46 km. This cluster showed 31 influenza A(H7N9) human cases between March 13, 2013 and April 11, 2013, accounting for 23.85% of the total human cases occurred in this phase (31/130), whereas only 17.36 could be expected.
Spatio-temporal cluster in phase II
There were totally 72 localities and 294 influenza A(H7N9) human cases contained in the second epidemic phase. The identified statistically significant spatio-temporal clusters of A(H7N9) human infections from June 1, 2013 to May 31, 2014 (364 days) were displayed in Fig. 5, in which the correlated influenza A(H7N9) human cases in this period were overlaid too. In the first 6 months of the second H7N9 epidemic phase (between June 1, 2013 and December 31, 2013), there was no statistically significant cluster detected in China, and then the space-time test identified three significant clusters (P < 0.05), one most likely cluster and two secondary clusters (Fig. 5 and Table 3). The most likely cluster (P = 0.0016) including 10 localities was identified in Zhejiang province, Jiangsu province and Shanghai municipality between January 12, 2014 and January 31, 2014 with a radius of 194.90 km, which showed 61 influenza A(H7N9) human cases accounting for 20.75% of all the human cases of phase II (61/294), when only 33.80 were expected theoretically. The secondary cluster 2 (P = 0.011) including 18 localities was identified in Anhui, Jiangsu, Shandong and Jiangxi provinces between April 22, 2014 and May 31, 2014, with a radius of 453.14 km, including 4.76% of the total human cases of this phase (14/294) whereas only 3.58 could be expected. The secondary cluster 3 (P = 0.028) including 4 localities was identified in Guangdong province between February 16, 2014 and March 2, 2014, with a radius of 80.49 km. This cluster included 6.46% of the total human cases of this phase (19/294) when only 6.53 could be expected theoretically.
Spatio-temporal cluster in phase III
Although there were still 36 influenza A(H7N9) human cases occurred in 21 localities during the third epidemic phase (213 days), there was no statistically significant spatio-temporal cluster of A(H7N9) human infections (P < 0.05) identified in phase III (June 1, 2014 to December 31, 2014) using the space-time permutation model.
By using the space-time permutation model, a total of five statistically significant spatio-temporal clusters (P < 0.05) were identified during the study period, in which two clusters were observed in phase I, another three clusters were observed in phase II, respectively, and there was no significant cluster occurred in phase III. In the three phases, around 29.35% of the total influenza A(H7N9) human cases were located in the five significant spatio-temporal clusters identified in this study (135/460), the radius of these clusters varied between 80.49 km and 596.94 km.
The emergence of A(H7N9) human infections in China has posed a great threat to both the society and the public. How H7N9 outbreaks spread is not clear till now. Although some H7N9 family clusters have been reported in China [30,31,32], so far there is no evidence of sustained or efficient human-to-human H7N9 virus transmission [33,34,35]. As a complementary instrument of public health surveillance, geographical analysis can investigate spatial distribution of a disease within a territory . In order to thoroughly examine the spatial-temporal spreading process of A(H7N9) human infections between March 2013 and December 2014, we investigated the spatial and temporal distribution, the directional trend and the statistically significant spatio-temporal cluster of A(H7N9) human infections by using a range of GIS and space-time statistical methods in this study.
A former research showed that some areas of Eastern China have a higher risk of infection with non-uniform distribution of H7N9 outbreaks . However, our results suggested that the high incidence areas of H7N9 epidemic were mainly located in the eastern and the southeast coast of China, such as Shanghai, Zhejiang, Guangdong and Jiangsu. This difference may be due to a longer study period in our paper and the epidemic has changed as time went on. It is noteworthy that the similar living and geomorphological environments in Shanghai, Zhejiang, Guangdong and Jiangsu may be suitable for the virus spreading. Since a host of the affected cities were the most densely developed middle- and large-sized cities, the population density may associate with the spread of H7N9, consisting with a previous research . In addition, we also noted that there were 9 new influenza A(H7N9) human cases occurred in Xinjiang between July 2014 and December 2014, where away from the eastern and the southeast coast of China and no case had reported before, this phenomenon is of concern and merits further investigation.
Influenza A(H7N9) human infections may have strong epidemic characteristics due to the seasonal pattern showed in Fig. 2. The epidemic was mainly concentrated between January and April, which suggested there was a possible seasonally higher risk of A(H7N9) human infections in China. This also provided some useful insights into H7N9 outbreaks in humans that the high incidence stage might coincide with the dry season months of January to April in China, and the occurrence of the epidemic peak was likely related to increased virus survival at rising temperature in this period. Some previous researches have highlighted that relative humidity and temperature may contribute to the occurrence of A(H7N9) human infections [15, 38], and an optimal temperature between 7 °C and 15 °C may be one of the driving forces for H7N9 . Consequently, three epidemic phases were identified in this study (Fig. 2): (1) phase I was the initial outbreak stage of A(H7N9) human infections which lasted 79 days, and the sudden outbreak of H7N9 epidemic then resulted in a large number of influenza A(H7N9) human cases (n = 130). However, it is notable that few human cases were reported in the summer of 2013 after some suitable public health measures such as closure of live poultry markets which announced and implemented by the government [40, 41], and this confirmed that suspension of live poultry trading was an effective measure to contain H7N9 outbreaks in humans . Furthermore, the existing researches also showed a possibility that live poultry markets contaminated by H7N9 virus were the most likely sources of the epidemic [6, 38, 43, 44]. And (2) in the beginning of phase II, there was an obvious reduction in the number of influenza A(H7N9) human cases and the low occurrence of the desease lasted nearly 6 months. But H7N9 outbreaks in humans resurged in the winter of 2013 , from January 2014 the number of the human cases began to sharply increase. And last (3) in phase III, only a slight rising but no obvious decline tendency of influenza A(H7N9) human cases’ number was observed due to the limitation of the research data in this study. This might because that the third phase was in an incubation period of the next epidemic curve, and the recurring H7N9 outbreaks in humans called for more improved preventive and control measures in China in the future.
In the first epidemic phase, the ellipse of influenza A(H7N9) human cases was the smallest one during the study period, which suggested that there were only some local H7N9 outbreaks in humans appeared in limited affected areas in these early days of the epidemic, especially some areas in eastern China. Moreover, the first epidemic phase displayed a more regular circular distribution than the others with a slight northwest-southeast elongation, and there was no obvious difference between the long axis and the short axis, which revealed that the directional trend of the epidemic was not apparent at the beginning. Although the spread of the outbreaks was limited then, there still appeared two significant spatio-temporal clusters of the epidemic. The secondary cluster 2 (P = 0.044) first appeared in phase I, with more influenza A(H7N9) human cases (n = 31), closer geographical distance between the human cases (radius = 271.46 km) and longer duration time (29 days) than the most likely cluster of phase I. An army of influenza A(H7N9) human cases gathered within small geographic scope in this cluster between March 13, 2013 and April 11, 2013, this was probably a consequence of the sudden outbreak of H7N9 and the lack of prevention measures in these early days of the epidemic. Subsequently, the most likely cluster (P = 0.0000055) appeared, which was the largest identified cluster (radius = 596.94 km) during the study period. Furthermore, over time the cluster shifted from the east to the southeast, which inferred that the epidemic might have a possible trend of spreading to the southeast coast of China.
In the second epidemic phase, both directional trend and local spatio-temporal clusters of H7N9 outbreaks in humans changed greatly. From the beginning of phase II (June 2013), influenza A(H7N9) human cases appeared sporadically for almost 6 months. But then a large number of influenza A(H7N9) human cases started appearing in the east and southeast of China between January 2014 and May 2014, which forming a larger ellipse and three significant spatio-temporal clusters of A(H7N9) human infections in this phase. The mean center (phase II) moved southwestward indicating that the H7N9 virus had spread along the southeast coast of China and hence the endemic area had shifted southwesterly. Furthermore, the standard deviational ellipse around the areas of A(H7N9) human infections substantially increased, suggesting that the infected areas had become larger than before. We also noted that the orientation of the ellipse had decreased from 152.621015 (phase I) to 26.947956 (phase II), indicating that the directional trend of A(H7N9) human infections had transformed from northwest-southeast to southwest-northeast direction. The long axis grew more significantly with an extending in the southwest-northeast directions and its length was nearly 2.7 times the short axis, which implied that the epidemic mainly spread along the southeast coast of China with a strong coastal growth. In addition, three statistically significant spatio-temporal clusters were observed in this phase, which shifted from the east to the southeast, and then back to the east-central China, and this suggested a possible epidemic trend of spreading to central China. The most likely cluster of phase II (P = 0.0016) contained the most influenza A(H7N9) human cases of all clusters detected in the three phases (n = 61). Then the secondary cluster 3 (P = 0.028) was observed in Zhaoqing, Foshan, Jiangmen and Guangzhou, and these four cities all belong to Guangdong province, suggesting that Guangdong was a high incidence area of A(H7N9) human infections in space - time dimension at that time. As the season progressed, the largest cluster of phase II developed in Anhui, Jiangsu, Shandong and Jiangxi provinces (Secondary cluster 2). Compared with the previous two clusters identified in phase II, this cluster showed the least influenza A(H7N9) human cases (n = 14), the largest geographical range (radius = 453.14 km) and the longest duration time (39 days) in this phase. This indicated that the number of the human cases decreased, but there was still an obvious trend of spatial diffusion of the epidemic.
However, there was no significant spatio-temporal cluster of A(H7N9) human infections identified during the third epidemic phase. This may be caused by a slice of sporadic influenza A(H7N9) human cases reported in the second half of 2014 and the relatively long period (213 days) with low H7N9 virus activity in this period. In phase III, the mean center of standard deviational ellipse moved from southeast to central regions, suggesting a further disease spreading throughout central China. Moreover, the size of the ellipse increased continuously, accordingly the geographical size of H7N9 outbreaks in humans reached the maximum value of the three phases. We observed that the orientation of the ellipse increased to 123.918912°, meanwhile, the long axis increased obviously with a northwest-southeast elongation and its length was nearly 3.66 times the short axis, all of this indicated that the directional trend of A(H7N9) human infections was more focused on the northwest-southeast with a strengthened spreading force, signaling a strong inland growth. Consequently, the significant reduction of the number of A(H7N9) human infections and the further expansion of the epidemic all contributed to the disappearance of the spatio-temporal cluster in phase III.
Although A(H7N9) human infections had not spread across the whole country, the epidemic moved from east to southeast coast, and thence to central China as the time changed, suggesting a future epidemiological trend that H7N9 outbreaks in humans would likely to continue and disperse to more central regions, and a few new human cases might appear in parts of the western China too. Meanwhile, the continuous increase in the size of the standard deviational ellipse signifying the relentless expansion of the infected area warned that the disease transmission cycle had not been interrupted yet. Surveillance and regional disease control preparedness should be strengthened accordingly. Furthermore, the disappearance of the cluster in phase III confirmed that the spatial-temporal clustering of A(H7N9) human infections in China between 2013 and 2014 presented different forms in different areas and times .
There are three potential limitations in this paper. Firstly, our analyses did not include the data of influenza A(H7N9) human cases after 2014 because which were either not available or imprecise. Secondly, this study has mainly examined the spatiotemporal spread patterns of A(H7N9) human infections in China between March 2013 and December 2014, further researches are still required to investigate the ecological risk and delineate the epidemiological nature of A(H7N9) human infections in the future. Thirdly, because the retrospective space-time permutation scan statistic model for detecting of disease outbreaks only uses case numbers and does not require the population-at-risk data, so the association of the spatio-temporal clusters of A(H7N9) human infections with population density has not been assessed in this study. However, the spatio-temporal clusters of the epidemic might still be associated with the population density of the study areas, which will be further studied in our future research. In addition, although we cannot get the complete and accurate data of influenza A(H7N9) human cases after 2014, the case reports in World Health Organization website  showed that H7N9 outbreaks in humans between January 2015 and August 2017 presented obvious seasonality and regionality: Winter–spring season was the most prevalent period for H7N9, and most cases were concentrated in the eastern and the southeast coast of China; with the spreading of the epidemic, human cases in central China have increased significantly, and a small number of new human cases have also been reported in parts of the western regions of China. These phenomena are consistent with the epidemiological trend predicted in this study, therefore our results are credible and reasonable, which can lay a foundation for not only risk-based surveillance but also investigations into risk factors associated with A(H7N9) human infections, and also can provide some meaningful information to future researches on assessing the transmission mechanism of H7N9 in China.
A(H7N9) human infections has essentially become a severe public health threat in China for several years. In order to explore the spreading dynamics of the epidemic, a phased assessment of the potential directional trend and the significant spatio-temporal clustering of H7N9 outbreaks in humans were investigated in this study. Consequently, three standard deviational ellipses and five significant spatio-temporal clusters (P < 0.05) of the epidemic were identified during the study period. The results showed that A(H7N9) human infections started in eastern China, persisted for a long period in southeast coastal areas, and then dispersed throughout the central regions, with a trend of continuous spatial expansion. Furthermore, the epidemic was space-time clustered during the first two phases in China, but it was randomly distributed in phase III. In order to limit the spread of the outbreaks, surveillance activities, closure of live poultry markets, and health education should be selectively taken within highly endemics areas in the epidemic season (especially between January and April). In short, this study can provide valuable references for establishing effective monitoring and preventive measures and controlling the virus spread during future H7N9 epidemic in China.
X coordinate for the mean center
Y coordinate for the mean center
The standard distance of X-axis
The standard distance of Y-axis
Cohen J. New flu virus in China worries and confuses. Science. 2013;340(6129):129–30.
Huang Y, Xu K, Ren DF, Ai J, Ji H, Ge AH, Bao CJ, Shi GQ, Shen T, Tang FY, et al. Probable longer incubation period for human infection with avian influenza a(H7N9) virus in Jiangsu Province, China, 2013. Epidemiol Infect. 2014;142(12):2647–53.
Zhang Y, Shen Z, Ma C, Jiang C, Feng C, Shankar N, Yang P, Sun W, Wang Q. Cluster of human infections with avian influenza a (H7N9) cases: a temporal and spatial analysis. Int J Environ Res Public Health. 2015;12(1):816–28.
World Health Organization, 2015. (http://www.who.int/csr/don/archive/country/chn/en/ ). Accessed 15 Jan 2015.
Gao R, Cao B, Hu Y, Feng Z, Wang D, Hu W, Chen J, Jie Z, Qiu H, Xu K. Human infection with a novel avian-origin influenza a (H7N9) virus. N Engl J Med. 2013;368(20):1888–97.
Chen Y, Liang W, Yang S, Wu N, Gao H, Sheng J, Yao H, Wo J, Fang Q, Cui D. Human infections with the emerging avian influenza a H7N9 virus from wet market poultry: clinical analysis and characterisation of viral genome. Lancet. 2013;381(9881):1916–25.
Li X-L, Liu K, Yao H-W, Sun Y, Chen W-J, Sun R-X, de Vlas SJ, Fang L-Q, Cao W-C. Highly pathogenic avian influenza H5N1 in mainland China. Int J Environ Res Public Health. 2015;12(5):5026–45.
Ding H, Xie L, Sun Z, Kao Q-J, Huang R-J, Yang X-H, Huang C-P, Wen Y-Y, Pan J-C, Pu X-Y, et al. Epidemiologic characterization of 30 confirmed cases of human infection with avian influenza A(H7N9) virus in Hangzhou, China. BMC Infect Dis. 2014;14(1):1-8.
Liu D, Shi W, Shi Y, Wang D, Xiao H, Li W, Bi Y, Wu Y, Li X, Yan J. Origin and diversity of novel avian influenza a H7N9 viruses causing human infection: phylogenetic, structural, and coalescent analyses. Lancet. 2013;381(9881):1926–32.
Ward MP, Maftei D, Apostu C, Suru A. Geostatistical visualisation and spatial statistics for evaluation of the dispersion of epidemic highly pathogenic avian influenza subtype H5N1. Vet Res. 2008;39(3).
Svensson K, Back E, Eliasson H, Berglund L, Granberg M, Karlsson L, Larsson P, Forsman M, Johansson A. Landscape epidemiology of tularemia outbreaks in Sweden. Emerg Infect Dis. 2009;15(12):1937–47.
Carnes A, Ogneva-Himmelberger Y. Temporal variations in the distribution of West Nile virus within the United States; 2000-2008. Appl Spat Ana Policy. 2012;5(3):211–29.
Yang K, Li W, Sun LP, Huang YX, Zhang JF, Wu F, Hang DR, Steinmann P, Liang YS. Spatio-temporal analysis to identify determinants of Oncomelania hupensis infection with Schistosoma Japonicum in Jiangsu province, China. Parasit Vectors. 2013;6:138.
Osmani MG, Ward MP, Giasuddin M, Islam MR, Kalam A. The spread of highly pathogenic avian influenza (subtype H5N1) clades in Bangladesh, 2010 and 2011. Prev Vet Med. 2014;114(1):21–7.
Qiu J, Li RD, XJ X, Hong XC, Xia X, CH Y. Spatiotemporal pattern and risk factors of the reported novel avian-origin influenza a(H7N9) cases in China. Prev VetMed. 2014;115(3–4):229–37.
Liu W, Yang K, Qi X, Xu K, Ji H, Ai J, Ge A, Wu Y, Li Y, Dai Q. Spatial and temporal analysis of human infection with avian influenza a(H7N9) virus in China, 2013. Euro surveillance: bulletin europeen sur les maladies transmissibles = European communicable disease bulletin. 2013;18(47):31–8.
National Fundamental Geographic Information System of China. 2015. (http://ngcc.sbsm.gov.cn/). Accessed 16 Feb 2015.
R Development Core Team, 2015. R Foundation for Statistical Computing, Vienna, . (http://www.R-project.org). Accessed 19 June 2015.
Scott LM, Janikas MV. Spatial statistics in ArcGIS. In: Fischer MM, Berlin GA, editors. Handbook of Applied Spatial Analysis: Software Tools, Methods and Applications. Edn. Heidelberg: Springer Berlin Heidelberg; 2010. p. 27–41.
ArcGIS Help 10.1. How Directional Distribution (Standard Deviational Ellipse) works. (http://resources.arcgis.com/en/help/main/10.1/index.html#/How_Directional_Distribution_Standard_Deviational_Ellipse_works/005p0000001q000000/). Accessed 26 Aug 2015.
Kulldorff, M., Information management services inc., 2015. (http://www.satscan.org/). Accessed 10 Aug2015.
Kulldorff M, Heffernan R, Hartman J, Assunçao R, Mostashari F. A space–time permutation scan statistic for disease outbreak detection. PLoS Med. 2005;2(3):216–24.
Liu Y, Wang X, Liu Y, Sun D, Ding S, Zhang B, Du Z, Xue F. Detecting spatial-temporal clusters of HFMD from 2007 to 2011 in Shandong Province, China. PLoS One. 2013;8(5):e63447.
Tuia D, Ratle F, Lasaponara R, Telesca L, Kanevski M. Scan statistics analysis of forest fire clusters. Commun Nonlinear Sci Numer Simul. 2008;13(8):1689–94.
Kulldorff M, Athas WF, Feurer EJ, Miller BA, Key CR. Evaluating cluster alarms: a space-time scan statistic and brain cancer in los Alamos, New Mexico. Am J Public Health. 1998;88(9):1377–80.
Yang K, Pan A, Yang Y, Zhang S, Ong S, Tang H. Remote sensing image registration using multiple image features. Remote Sensing. 2017;9(6):581.
Wei Z, Han Y, Li M, Yang K, Yang Y, Luo Y, Ong S-H. A small UAV based multi-temporal image registration for dynamic agricultural terrace monitoring. Remote Sensing. 2017;9(9):904.
Ma J, Chen C, Li C, Huang J. Infrared and visible image fusion via gradient transfer and total variation minimization. Information Fusion. 2016;31:100–9.
Wang J, Cao Z, Zeng DD, Wang Q, Wang X, Qian H. Epidemiological analysis, detection, and comparison of space-time patterns of Beijing hand-foot-mouth disease (2008–2012). PLoS One. 2014;9(3):e92745.
Xinhua News, 2013. Expert review: the clinical symptoms between the father and the son with H7N9 virus infection in Shandong province were obviously different, difficult to support the probability of human-to-human transmission of H7N9 virus. (http://news.china.com/focus/h7n9/11138618/20130430/17808678.html). Accessed 11 Aug 2013.
Ding H, Chen Y, Yu Z, Horby PW, Wang F, Hu J, Yang X, Mao H, Qin S, Chai C, et al. A family cluster of three confirmed cases infected with avian influenza a (H7N9) virus in Zhejiang Province of China. BMC Infect Dis. 2014;14(1):1–14.
Liu T, Bi Z, Wang X, Li Z, Ding S, Bi Z, Wang L, Pei Y, Song S, Zhang S, et al. One family cluster of avian influenza a(H7N9) virus infection in Shandong, China. BMC Infect Dis. 2014;14(1):1–6.
Nishiura H, Mizumoto K, Ejima K. How to interpret the transmissibility of novel influenza a (H7N9): an analysis of initial epidemiological data of human cases from China. Theor Biol Med Model. 2013;10(1):30.
Cowling BJ, Jin L, Lau EH, Liao Q, Wu P, Jiang H, Tsang TK, Zheng J, Fang VJ, Chang Z. 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.
Parry J. H7N9 avian flu kills seven and infects 23 in China. BMJ. 2013;346
Gomez-Barroso D, Martinez-Beneito M, Flores V, Amorós R, Delgado C, Botella P, Zurriaga O, Larrauri A. Geographical spread of influenza incidence in Spain during the 2009 a (H1N1) pandemic wave and the two succeeding influenza seasons. Epidemiol Infect. 2014;142(12):2629–41.
Pan Y-N, Lou J-J, Han X-P. Outbreak patterns of the novel avian influenza (H7N9). Physica A: Stat Mech Appl. 2014;401:265–70.
Fang L-Q, Li X-L, Liu K, Li Y-J, Yao H-W, Liang S, Yang Y, Feng Z-J, Gray GC, Cao W-C. Mapping spread and risk of avian influenza a (H7N9) in China. Sci Rep. 2013;3:2722.
Hu W, Zhang W, Huang X, Clements A, Mengersen K, Tong S. Weather variability and influenza a (H7N9) transmission in shanghai, China: a Bayesian spatial analysis. Environ Res. 2015;136:405–12.
Gong ZY, Lv HK, Ding H, Han JK, Sun JM, Chai CL, Cai J, Yu Z, Chen EF. Epidemiology of the avian influenza a (H7N9) outbreak in Zhejiang Province, China. BMC Infect Dis. 2014;14
Goodwin R, Sun S. Early responses to H7N9 in southern mainland China. BMC Infect Dis. 2014;14(1):1–7.
Abdelwhab EM, Veits J, Mettenleiter TC. Prevalence and control of H7 avian influenza viruses in birds and humans. Epidemiol Infect. 2014;142(5):896–920.
Chen E, Chen Y, Fu L, Chen Z, Gong Z, Mao H, Wang D, Ni M, Wu P, Yu Z. Human infection with avian influenza a (H7N9) virus re-emerges in China in winter 2013. Euro Surveill. 2013;18(43)
Ma X, Liao Q, Yuan J, Liu Y, Liu Y, Chen J, Liu J, Cai W, Cowling BJ, Di B, et al. Knowledge, attitudes and practices relating to influenza a(H7N9) risk among live poultry traders in Guangzhou City, China. BMC Infect Dis. 2014;14(1):1–12.
Availability of data and materials
The datasets analysed during the current study are available from the first author Wen Dong (email address: email@example.com) on reasonable request.
The study was supported by two projects under grant Number 2012AA1214002 and 2012AA1214003, the Hi-Tech Research and Development Program of China (grant No. 2012AA1214000), the Specialized Research Fund for the Doctoral Program of Higher Education (grant No. 20115303110002), the Yunnan Provincial Science and Technology Project (grant No. 2011XX2005), the Yunnan Provincial Department of Education Research Fund (grant No. 2011Y307), the National Natural Science Foundation of China (grant No. 41461038) and the National Natural Science Foundation of China (grant No. 41661087).
Ethics approval and consent to participate
This study was undertaken with the retrospective analysis of data collected from the World Health Organization website, and the data used in the analysis were anonymized and no additional contact was made with patients to gain further information as a part of this study, therefore, no ethical approval was required.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.