key: cord-0000822-zss38mct authors: Fang, Li-Qun; Wang, Li-Ping; de Vlas, Sake J.; Liang, Song; Tong, Shi-Lu; Li, Yan-Li; Li, Ya-Pin; Qian, Quan; Yang, Hong; Zhou, Mai-Geng; Wang, Xiao-Feng; Richardus, Jan Hendrik; Ma, Jia-Qi; Cao, Wu-Chun title: Distribution and Risk Factors of 2009 Pandemic Influenza A (H1N1) in Mainland China date: 2012-05-01 journal: Am J Epidemiol DOI: 10.1093/aje/kwr411 sha: 043aa68f2c784899f71225e8eb233150760a6a54 doc_id: 822 cord_uid: zss38mct Data from all reported cases of 2009 pandemic influenza A (H1N1) were obtained from the China Information System for Disease Control and Prevention. The spatiotemporal distribution patterns of cases were characterized through spatial analysis. The impact of travel-related risk factors on invasion of the disease was analyzed using survival analysis, and climatic factors related to local transmission were identified using multilevel Poisson regression, both at the county level. The results showed that the epidemic spanned a large geographic area, with the most affected areas being in western China. Significant differences in incidence were found among age groups, with incidences peaking in school-age children. Overall, the epidemic spread from southeast to northwest. Proximity to airports and being intersected by national highways or freeways but not railways were variables associated with the presence of the disease in a county. Lower temperature and lower relative humidity were the climatic factors facilitating local transmission after correction for the effects of school summer vacation and public holidays, as well as population density and the density of medical facilities. These findings indicate that interventions focused on domestic travel, population density, and climatic factors could play a role in mitigating the public health impact of future influenza pandemics. In early April 2009, human cases of infection with 2009 pandemic influenza A (H1N1) virus were first identified in the United States and Mexico (1) . The virus then spread rapidly to other regions of the world. As of January 24, 2010, laboratory-confirmed cases of pandemic influenza (H1N1-2009) were being reported in more than 209 countries or regions worldwide, with 14,711 deaths among confirmed cases (2) . Pandemic influenza was introduced to mainland China on May 9, 2009 (3, 4) , and then spread across the whole country. By the end of 2009, more than 120,000 confirmed cases were reported to the Chinese Center for Disease Control and Prevention (CCDC), including 648 deaths (5) . Information on reported cases was released daily by the Ministry of Health of the People's Republic of China in the early stages of the pandemic and then twice weekly later on. Analyzing the information gathered and unearthing underlying risk factors provide an opportunity to identify epidemic characteristics and transmission patterns of the pandemic in China, thereby producing useful information for prevention and control measures during future epidemics. In this study, we aimed to characterize the temporal and spatial distribution of pandemic influenza in mainland China, to understand the diffusion pattern of the disease within the country, and to identify risk factors for invasion and local transmission of this disease. We used a database that included all cases of pandemic influenza (H1N1-2009) reported to the China Information System for Disease Control and Prevention (CISDCP) from May 9, 2009 , when the first confirmed case in China was reported, to December 31, 2009 (6) . The CISDCP covers all provincial, prefectural, and county centers for disease control and prevention, 95.3% of the provincial, prefectural, and county hospitals (9,084 in total), and 84.0% of township clinics (38,175 in total) across mainland China. After the World Health Organization issued an alert about the novel influenza virus (H1N1-2009), pandemic influenza was classified as a class B notifiable infectious disease on April 30, 2009 , by the Ministry of Health but was managed according to the criteria for class A notifiable infectious diseases. According to the Law for Prevention and Control of Infectious Diseases in China, information regarding each patient, once identified, must be reported to the CCDC within 2 hours through the Web-based CISDCP system. A suspected case was defined as a person with influenza-like symptoms who had had close contact with a confirmed case within the past 7 days, had a history of travel to affected areas within the past 7 days, or tested positive for influenza A virus, excluding other known subtypes of influenza A. A laboratory-confirmed case was defined according to World Health Organization criteria-that is, a person with influenza-like symptoms and laboratory-confirmed pandemic influenza A virus infection by one or more of the following tests: reverse-transcriptase polymerase chain reaction, real-time reverse-transcriptase polymerase chain reaction, viral culture, or a 4-fold rise in specific antibodies to pandemic influenza A virus (7) . Influenzalike symptoms were defined according to World Health Organization criteria: sudden onset of fever greater than 38°C, cough or sore throat, and absence of other diagnoses (8) . All laboratory-confirmed cases from 2009 were included in our database, including information about age, sex, occupation, residence address, work address, onset date and location, hospital admission date and address, and clinical outcome. Furthermore, census information was obtained from the National Bureau of Statistics of China (9) . Epidemic curves were created by plotting the daily number of newly confirmed cases and the number of deaths. Mainland China is divided into 2,925 counties, which are political subdivisions of provinces, usually containing several townships. The incidences for different sex and age groups were calculated using 2009 census data. As was done previously (10, 11) , each case was linked to a digital map of China (1:100,000) according to its onset location. The incidence for each county was calculated and standardized using direct standardization for age and sex according to the overall composition of the 2009 Chinese population, using 5-year age-group categories. To explore the spatial and temporal diffusion trend of pandemic influenza in mainland China, a map of pandemic influenza spread was developed using trend surface analysis, which is a spatial smoothing method that uses polynomials with geographic coordinates, as defined by the central point of each county (12) (13) (14) . The time delay of the first confirmed case for each county was defined as the duration of time (in days) since May 9, 2009, the date of the first confirmed case in mainland China. A trend surface on these durations was created in ArcGIS 9.2 (ESRI Inc., Redlands, California) using a second-order trend surface model with a local polynomial method to explore the diffusion patterns of pandemic influenza over time. In addition, the epidemic curves were also plotted for 5 selected cities/provinces in different regions of mainland China: Guangdong (southern China), Shanghai (eastern China), Shaanxi (central China), Tibet (western China), and Beijing (northern China). To assess the association between the invasion of pandemic influenza and different means of domestic travel, we performed a survival analysis (Cox analysis) of the time delay of the first confirmed case for each affected county, considering unaffected counties as right-censored. The time delay of the first confirmed case was defined as the duration (in days) since May 9, 2009, the date of the first confirmed case in mainland China. Domestic travel was expressed using 4 variables: distance (from the midpoint of each county) to the nearest civil airport and whether or not a county was intersected by national highways, freeways, or railways. This information was obtained from the National Bureau of Statistics of China (9, 15) . Spatial analyses were used to extract data on these variables in ArcGIS 9.2 (ESRI Inc.). Since population density is also linked to human activities and may facilitate influenza transmission (16, 17) and since density of medical facilities could be linked to patient reporting, we adjusted for the effects of these variables. In this study, the population density for each county was obtained from the National Bureau of Statistics of China (9, 15) , and the densities of medical facilities were based on the CISDCP, including all reporting sectors for notifiable infectious diseases, comprising provincial, prefectural, and county centers for disease control and prevention and provincial, prefectural, and county hospitals and township clinics. Hazard ratios and their 95% confidence intervals and P values were estimated using maximum likelihood methods. Hazard ratios for the continuous variables were calculated for the following units: distance to the nearest airport (in 50-km increments), population density (in 1,000 persons per km 2 ), and density of medical facilities (in number of reporting sectors for pandemic influenza per 10,000 persons). To explore the effect of climatic factors on local transmission within counties, we performed multilevel Poisson regression. Climatic data (temperature, relative humidity, and precipitation) during May-December 2009 were obtained from the National Meteorological Bureau of China (18) . Owing to probable time lags, the climatic variables were processed by calculating the average value for the current day and a lag of 1-3 days, which is the observed incubation period of pandemic influenza (19) . Poisson regression deals with the daily number of laboratoryconfirmed cases per county. The inclusion of the population size for each county as an offset makes it an analysis of incidence. To account for possible confounding, we included school summer vacation and public holidays, the proportion of the school-age population (ages 6-19 years), population density, and the density of medical facilities as correction factors in the analysis. The percentage change in incidence in response to the change of the variable by a given amount (10°C for temperature, 10% for relative humidity, 1 mm for precipitation, 10% for school-age population, 1,000 persons per km 2 for population density, and number of facilities per 10,000 persons) was used to reflect the impact of each variable. The 95% confidence intervals and corresponding P values were estimated after correcting for overdispersion because of the nature of infectious diseases with spatial clustering patterns (20, 21) . For temperature, we also included a quadratic term in the analysis. For all analyses, univariate analysis was performed first to examine the effect of each variable separately. Multivariate analysis was then performed by including all variables with P values less than 0.20 in univariate analysis and exclusion of variables with P values greater than 0.10, using a standard backward likelihood ratio method. For all continuous variables, we also presented categorical results in 3-5 categories to allow inspection of the data and determine whether or not the assumption regarding continuous variables (quadratic for temperature) was justified. Statistical analyses were performed using the Stata package (StataCorp LP, College Station, Texas) (20) . Readers interested in further research can contact the corresponding author to obtain the full data set used in this study. A total of 121,805 cases of pandemic influenza (H1N1-2009), distributed in all 31 provinces in mainland China, were reported from May 9, 2009, to December 31, 2009. There was much variation in the numbers of confirmed cases in different provinces, ranging from 881 to 12,748, with a median of 2,958 cases, and in the incidence of confirmed cases in different provinces, ranging from 3.94 per 100,000 population to 71.72 per 100,000, with a median of 8.41 per 100,000. From the time profile, we found that the number of confirmed cases increased rapidly beginning at the end of August, when a new term began for school students, and peaked by the end of November. The first death caused by pandemic influenza was reported on October 4, 2009. The number of deaths eventually rose to 648 by the end of the year, and peaked in early December (Figure 1 ). The age-and sex-standardized incidence map shows that the epidemic spanned a large geographic area, and the most affected areas were in western China (see Web Figure 1 , which appears on the Journal's website (http://aje.oxfordjournals.org/)). Significant differences in incidence were found among age groups, with incidences peaking in school-age groups (Web Figure 1 ). Boys showed a higher incidence than girls (ages <20 years). Web Figure 2 shows the trend of the spatial spread of pandemic influenza over time in mainland China and indicates that the epidemic areas during the first 120 days after May 9, 2009, were limited to the circumferences of cities with international airports, such as Beijing, Shanghai, Guangzhou, Shenzhen, Chengdu, and Changchun. Thereafter, it spread to the rest of mainland China, roughly from southeast to northwest. The largest-scale spread took place 150-180 days after the first case (Web Figure 2) . Figure 2 shows the large variation in the temporal patterns of pandemic influenza for the 5 selected cities/provinces, but there was a marked drop in incidence during the first week of October for all locations. Survival analysis of the duration of time to the first confirmed case in each county indicated that all 4 factors related to domestic travel or human mobility were significantly associated with the invasion of pandemic influenza in the Cox univariate analysis (Table 1) . Population density and the density of medical facilities also showed a significant association. The significant effect of being intersected by railways disappeared, and the density of medical facilities showed borderline significance after correction for other factors in multivariate analysis, whereas being intersected by national highways and freeways and proximity to airports and higher population den-sity remained as significant factors, all showing a positive association (Table 1) . Table 2 shows that all climatic factors (except precipitation), school summer vacation and public holidays, proportion of the school-age population, population density, and the density of medical facilities were significantly associated with the extent of local transmission in univariate multilevel Poisson regression. School summer vacation and public holidays showed a significant negative association with the incidence of pandemic influenza. The significant effect of the proportion of schoolage children disappeared after correction for other factors; thus, temperature, relative humidity, school summer vacation and public holidays, population density, and the density of medical facilities remained as significant factors in multivariate analysis. Temperature showed a peak pattern, with the highest incidences for the range from 0°C to 10°C, which was also reflected in the statistically significant quadratic term. Our study provides a complete overview of the spatial and temporal characteristics of the pandemic influenza (H1N1-2009) epidemic in mainland China in 2009. The epidemic spanned a large geographic area and presented spatial and temporal heterogeneity in different regions of mainland China. Our analyses of the invasion of pandemic influenza indicated that domestic travel by air and by national highways and freeways and population density contributed to the spread of the epidemic. Lower temperatures and lower relative humidity were climatic factors that facilitated local transmission after correction for the effects of school summer vacation and public holidays, as well as population density and the density of medical facilities. The density of medical facilities could have influenced pandemic influenza patient reporting, and this effect seemed more important for the reporting of local transmission (a highly significant positive association) than for reporting of the invasion (borderline significance). This indicates that the CISDCP can be further improved. In the initial phase of the epidemic, the Chinese government took measures to prevent and control the spread of the novel influenza virus, declaring it a notifiable infectious disease in order to strengthen national surveillance and find newly confirmed cases quickly. In addition, quarantine measures were implemented at the international airports (e.g., Beijing, Shanghai, Guangzhou, and Fuzhou) to identify and isolate prob-able cases and close contacts in order to decrease the risk of local transmission caused by imported cases. It is possible that these measures were effective in achieving a slower pace of spread in the early stages of the epidemic, but the current evidence is inconclusive (22) . However, following the development of the global pandemic and the beginning of the new school term, a rapid spread of the epidemic occurred in mainland China, and local outbreaks increased at the end of August 2009. A reduction in incidence was observed during the first week of October, when there was an 8-day public holiday during the National Days from October 1 to October 8. This drop in incidence was largely due to a lower tendency for patients to visit medical facilities at that time, together with the fact that many hospitals had reduced the open hours of their outpatient clinics during the holidays. This is clearly visible in Figure 1 , where we see the drop beginning on September 28, 3 days (i.e., the average duration between onset and seeing a physician) before the start of the holiday period. Thus, the holiday period led to a reduction in the number of people being diagnosed with pandemic influenza. Apparently, most undiagnosed patients recovered in the following days, because there was no marked compensation visible in the days after the holiday period. Additionally, there Abbreviations: CI, confidence interval; NS, not significant. a Results were adjusted for school summer vacation and public holidays, population density, and the density of medical facilities. b For all continuous variables, categorical results are also reported to allow inspection of the data and assessment of whether or not the assumption regarding continuous variables was justified. may have been some reduced transmission because of school closure, as was observed in Japan, where transmission was substantially reduced during school closure (23) . In addition, the temporal death curve could reflect the massive rise in confirmed cases with approximately 1 week's delay following the peak of confirmed cases at the end of November. The direction of the spread of pandemic influenza was from the southeast to the northwest, indicating how the virus benefited from entering the country through international airports in the coastal areas and spreading further along routes of long-distance domestic travel. With the fast-growing public transportation infrastructure and increasing socioeconomic activities, travel has become an important issue in the prevention of emerging airborne infectious diseases such as influenza, especially during the introduction period. Obviously, the presence of airports and high densities of transportation routes coincide with more developed areas (i.e., those with a higher population density and more medical facilities); however, after we corrected for these factors, proximity to airports and the presence of national highways or freeways remained significantly associated with the spread of the infection. Geoinformatics plays an important role in the study and control of infectious disease outbreaks, and it includes techniques such as geographic mapping and location-based alert services (16, (24) (25) (26) . As was recognized previously, the international spread of pandemic influenza and severe acute respiratory syndrome was largely related to air travel (27) (28) (29) . Our study confirms that air travel and transportation routes accelerated the spread of pandemic influenza between counties in mainland China. Air travel and travel by national highways and freeways especially appeared to play a role, whereas railways were less important. In mainland China, trains are mainly used for occasional long-distance travel, whereas highways are more often used for daily or weekly commuting, especially because bus schedules are more flexible. Our previous study on the geographic spread of the severe acute respiratory syndrome epidemic in China also demonstrated that domestic travel along national highways played a more important role than travel by railway (30) . Data from the National Bureau of Statistics of China show that passenger traffic by highways in 2009 was 18.2 times that of railway travel (9) . Transportation by highway remains an important mode of travel between Chinese cities and provinces and therefore is a potential target for controlling any future emerging airborne infections. We also showed that lower temperature and lower relative humidity create a higher risk of local transmission of pandemic influenza. However, much lower temperature (e.g., <0°C) did not facilitate local transmission, as indicated by the daily incidence over categorical temperature in Table 2 , and required inclusion of a quadratic variable in the model. The observation of an influence of temperature and relative humidity on pandemic influenza is in accordance with animal experiments on seasonal influenza virus (31) . In addition, recent studies have suggested that absolute humidity could also play an important role in the transmission of pandemic influenza and seasonal variations in influenza epidemics in temperate regions (32) (33) (34) . As expected, population density further facilitated both invasion and local transmission, whereas holiday periods reduced spread (16, 17, 35, 36) . We also used population size as a correction factor instead of population density, which led to similar results (not shown). In conclusion, this is the first complete documentation of pandemic influenza in mainland China, to our knowledge. The findings indicate that interventions focused on domestic travel, population density, and climatic factors could play a major role in mitigating the public health impact of future influenza pandemics. Emergence of a novel swineorigin influenza A (H1N1) virus in humans. Novel Swine-Origin Influenza A (H1N1) Virus Investigation Team Pandemic (H1N1) 2009-Update 85 Clinical and epidemiologic characteristics of 3 early cases of influenza A pandemic (H1N1) 2009 virus infection, People's Republic of China Clinical features of the initial cases of 2009 pandemic influenza A (H1N1) virus infection in China. National Influenza A Pandemic (H1N1) 2009 Clinical Investigation Group of China Ministry of Health of the People's Republic of China. Announcement of the Situation of Prevention and Control of Pandemic Influenza A (H1N1) in China Ministry of Health of the People's Republic of China Chinese Center for Disease Control and Prevention. China Information System for Disease Control and Prevention [database Chinese Center for Disease Control and Prevention Ministry of Health of the People's Republic of China. Protocol of Diagnosis and Treatment of Pandemic Influenza A (H1N1) (Preliminary Version) Ministry of Health of the People's Republic of China World Health Organization China Statistical Yearbook Spatiotemporal dynamics of hemorrhagic fever with renal syndrome, Beijing, People's Republic of China Environmental factors contributing to the spread of H5N1 avian influenza in mainland China Spatial analysis of sleeping sickness, southeastern Uganda Spatial diffusion of raccoon rabies in Pennsylvania, USA Estimation and prediction in the spatial linear model China Data Online [database Synchrony, waves, and spatial hierarchies in the spread of influenza Monitoring the impact of influenza by age: emergency department fever and respiratory complaint surveillance in New York City China Meteorological Administration The transmissibility and control of pandemic influenza A (H1N1) virus A Handbook of Statistical Analyses Using STATA Spatiotemporal trends and climatic factors of hemorrhagic fever with renal syndrome epidemic in Shandong Province, China China's tough flu measures appear to be effective Influenza (H1N1) 2009 outbreak and school closure Improving imperfect data from health management information systems in Africa using space-time geostatistics Need for improved methods to collect and present spatial epidemiologic data for vectorborne diseases Large-scale spatial-transmission models of infectious disease Spread of a novel influenza A (H1N1) virus via global airline transportation The severe acute respiratory syndrome: impact on travel and tourism Experience of severe acute respiratory syndrome in Singapore: importation of cases, and defense strategies at the airport Geographical spread of SARS in mainland China Influenza virus transmission is dependent on relative humidity and temperature Absolute humidity and pandemic versus epidemic influenza Absolute humidity and the seasonal onset of influenza in the continental United States Absolute humidity modulates influenza survival, transmission, and seasonality Strategies for mitigating an influenza pandemic Estimating the impact of school closure on influenza transmission from sentinel data The authors thank the medical staff who contributed to the ascertainment of pandemic influenza (H1N1-2009) cases, including variable reporting coverage across time, age groups, and locations. They also thank Dr. Caspar W. N. Looman at the Department of Public Health of Erasmus Medical Center, University Medical Center Rotterdam, for his comments on the statistical analyses and interpretation of the results.Conflict of interest: none declared.