Next Article in Journal
Modifying Health Behavior to Prevent Cardiovascular Diseases: A Nationwide Survey among German Primary Care Physicians
Previous Article in Journal
Managing Terrorism or Accidental Nuclear Errors, Preparing for Iodine-131 Emergencies: A Comprehensive Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Japanese Encephalitis Risk and Contextual Risk Factors in Southwest China: A Bayesian Hierarchical Spatial and Spatiotemporal Analysis

1
West China School of Public Health, Sichuan University, Chengdu 610041, China
2
Department of Biostatistics, School of Public Health, University of Washington,Seattle, WA 98195, USA
3
School of Public Health, Xinjiang Medical University, Urumqi 830011, China
4
Office for Disease Control and Emergency Response, Chinese Center for Disease Control and Prevention (China CDC), Beijing 102206, China
5
HSR&D Center of Excellence, VA Puget Sound Health Care System, Seattle, WA 98101, USA
*
Authors to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2014, 11(4), 4201-4217; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph110404201
Submission received: 9 January 2014 / Revised: 12 March 2014 / Accepted: 4 April 2014 / Published: 15 April 2014

Abstract

: It is valuable to study the spatiotemporal pattern of Japanese encephalitis (JE) and its association with the contextual risk factors in southwest China, which is the most endemic area in China. Using data from 2004 to 2009, we applied GISmapping and spatial autocorrelation analysis to analyze reported incidence data of JE in 438 counties in southwest China, finding that JE cases were not randomly distributed, and a Bayesian hierarchical spatiotemporal model identified the east part of southwest China as a high risk area. Meanwhile, the Bayesian hierarchical spatial model in 2006 demonstrated a statistically significant association between JE and the agricultural and climatic variables, including the proportion of rural population, the pig-to-human ratio, the monthly precipitation and the monthly mean minimum and maximum temperatures. Particular emphasis was placed on the time-lagged effect for climatic factors. The regression method and the Spearman correlation analysis both identified a two-month lag for the precipitation, while the regression method found a one-month lag for temperature. The results show that the high risk area in the east part of southwest China may be connected to the agricultural and climatic factors. The routine surveillance and the allocation of health resources should be given more attention in this area. Moreover, the meteorological variables might be considered as possible predictors of JE in southwest China.

1. Introduction

Japanese encephalitis (JE) is a mosquito-borne viral disease, which is widely distributed in South Asia, Southeast Asia, East Asia and the Pacific [1]. It is estimated that about three billion people live in countries where the JE virus is endemic [2]. JE not only causes death, but also induces permanent and psychiatric sequelae [3]. Historically, JE was serious in China. The number of JE cases has declined substantially after the long-term nationwide immunization program starting in the 1970s. However, JE still remains a significant public health issue in China, with approximately a half of the global cases annually [4], and JE virus is one of the four principal arboviruses of public health importance in China [5]. In the past few decades, the Chinese Center for Disease Control and Prevention reported that the annual numbers of death from JE ranked 2–6 among the 26 most serious national notifiable communicable diseases [6,7].

As a mosquito-borne arboviral infection, JE is caused by a flavivirus transmitted by mosquitoes of the Culex species. The spread of JE involves many ecological, environmental, climatic and human behavioral factors, of which two groups are significant; these are climatic factors and agricultural factors. On the one hand, climatic factors can influence the transmission of JE through their effects on mosquitoes. Temperature and precipitation have been reported to be associated with the density of mosquitoes [8]. Although suitable climatic condition could increase JE vector proliferation and longevity, the potential impact of climate change on JE remains to be investigated [2]. Particularly, like other vector-borne diseases, two common questions for JE could be whether the association with meteorological factors exists and if it does exist, what is the time-lagged pattern for the association [9]. On the other hand, two major agricultural factors affecting the JE transmission are pig rearing and the rice-irrigated areas [1,10]. Domestic pigs are the most important reservoir hosts, amplifying the transmission process [3]. The rural population has a greater chance for exposure to mosquitoes compared with the general population [11]. In addition, other socioeconomic variables may affect the transmission process, such as population immunity and housing conditions [12].

There is no specific therapy for JE other than supportive care; therefore, a better understanding of the impact from contextual factors is beneficial to prevent JE, and several studies have been conducted. In [11], it was demonstrated that at the provincial level, the risk of JE seemed to be governed by high proportions of rice-planting area, rural population and the extent of pig rearing. The authors of [13,14,15] studied the impact of the climate variability on the transmission of JE. These studies give evidence on the association between JE cases and the contextual risk factors, including agricultural and climatic factors. However, several limitations and gaps are noticeable in those studies. First, no research has been conducted to assesses the association between JE cases and the contextual risk factors in southwest China, which is the most endemic area in China [16]. The aforementioned single-site studies are in eastern China, and their results may not be valid in southwest China. Meanwhile, a better understanding of the spatiotemporal pattern of JE in southwest China would help identify areas and populations at high risk in southwest China. Second, existing studies focused on only one site to assess the association between JE and climatic variables [17]. There are two drawbacks for this approach. One is that the conclusion is hard to generalize to other sites, resulting in some inconsistent findings between the single-site studies. For example, while [14] reported that a one-month lag is the best time lag for monthly rainfall, [13,18] reported that a two-month lag is the best time lag for monthly rainfall. The other drawback is that it is not sufficient to analyze the association between meteorological factors without adjusting for agricultural variables [13,14,15]. In summary, there is no large-scale analysis (such as over 300 counties) at a small unit (such as the county level) in China to study the association between JE and the contextual factors. To the best of our knowledge, there are few large-scale studies in other countries, with the exception being in Nepal [19]. The study in Nepal analyzed the data in 2005 at the district-level to fit a spatial lag regression model with climatic, agricultural and land-cover variables.

Our aim was to assess the association between JE cases and contextual risk factors in southwest China and to identify the high risk area. This may help in understanding JE epidemiology and provide guidance for JE prevention. Specifically, we applied the Bayesian hierarchical model to explore the spatiotemporal pattern of JE in 438 counties in southwest China from 2004 to 2009. Moreover, the spatial ecological regression was conducted between JE and contextual risk factors (including agricultural and climatic variables) in 2006. This is due to the large number of cases that occurred in that year. Particular attention was given to the lag effect for climatic factors.

2. Materials and Methods

2.1. Study Region

Southwest China (21°14′ to 34°31′N, 97°35′ to 110°19′E) consists of four provinces, Sichuan, Chongqing, Yunnan and Guizhou. The area has a population of 189,977,077 (the sixth national census in 2010) and encompasses 1,137,570 square kilometers. There are 47 cities (autonomous prefectures), which can be further divided into 483 counties (county-level cities and districts).

2.2. Data Collection

JE monthly cases were obtained from 2004 to 2009 for the 483 counties, facilitated by the Chinese Information System for Infectious Diseases Control and Prevention (CISIDCP), which was put in place in 2004 and was more sensitive and efficient than the previous case-reporting system.

Population data for every county from 2004 to 2009 were retrieved from the National Bureau of Statistics of China. The proportion of rural population was from the statistical yearbook of these four provinces. The pig-rearing number came from the census statistics of the four provincial Bureau of Statistics at the prefecture-city level, an administrative division below a province and above a county in China's administrative structure.

County-level Geographical Information System (GIS) shapefiles of southwest China were acquired from the Chinese Institute of Geographic Sciences and Natural Resources Research.

The meteorological data contained monthly precipitation and mean temperatures, including mean minimum and maximum temperatures. They were collected from the 0.5° × 0.5° grid monthly precipitation and monthly mean temperature datasets, publicly available at the Chinese Meteorological Data Sharing Service System “ http://cdc.cma.gov.cn/home.do”, which were constructed by the Chinese National Meteorological Information Center. They are meant to serve as a high resolution source of climate data interpolated from various climate sources [20,21]. We retrieved three monthly climatic variables in 2006 in the format of the 0.5° × 0.5° grid covering the entirety of China, which can be matched to each county. In general, each county only contained one grid point, and all counties took the climatic information of the grid point closest to their administrative centroids.

2.3. Exploratory Analysis

2.3.1. GIS Mapping

All JE cases were geo-coded and matched to the county-level maps by administrative codes, and we calculated annualized incidences of JE at each county over the 6 years. For the descriptive ecological correlation analysis in 2006, the preceding geo-coding was also executed for the agricultural and climatic data. Finally, this information was visualized in the form of maps. All of the implementation above, including the spatial visualization, was accomplished by R. R is a free software programming language and a software environment for statistical computing and graphics [22].

2.3.2. Spatial Autocorrelation Analysis

Global spatial autocorrelation analyses were performed for each year to discern the spatial autocorrelation of JE in southwest China after adjusting for population heterogeneity. The classic and simplest boundary-based neighborhood matrix was used. Thus, the weight is 1 if two counties share a boundary; otherwise, the weight is 0. The analysis was accomplished using R via the Assunção and Reis's global Moran's I statistics [23] in R package “spdep”.

2.3.3. Spearman Correlation Analysis

We examined the relationship between monthly incidences of JE and monthly climatic variables over 438 counties. Spearman's correlation was used to quantify the relationship with lags from 1 to 3 months.

2.4. Bayesian Hierarchical Spatial and Spatiotemporal Modeling

Traditional regression methods assume that the observations are mutually independent, which is not valid, due to the spatial structure. Bayesian hierarchical spatial or spatiotemporal models [24,25] do not make such an assumption, inducing spatial correlation through random effects representing unmeasured (or perhaps unmeasurable) effects not included in the model. It can also take into account the variation of population and overdispersion [26]. These models are extensively applied to epidemiology [27].

2.4.1. Spatiotemporal Model for County-Specific Relative Risks from 2004–2009

The first stage describes the observed data, yijk, the number of cases in county j in month k of the i-th year, as a function of the spatiotemporal-specific relative risk of disease, θijk, and the expected cases, eijk. The expected cases, eijk, were calculated by the internal standardization method, assuming that the disease risk was constant over southwest China. A zero inflated Poisson model was employed to model the observed data, as the county × month combination would lead to excessive zeros compared to the simple Poisson model [28].

y ijk { π + ( 1 π ) × Poisson ( y ijk = 0 | e ijk θ ijk ) if y ijk = 0 , ( 1 π ) × Poisson ( y ijk | e ijk θ ijk ) if y ijk = 0 .
π is the probability of a zero count, and the log-relative risks log θijk are assumed to have the decomposition:
log θ ijk = μ + ϕ j + ψ j + α i + γ i
where:
μ:the intercept quantifying the average Poisson relative risk in the whole of southwest China
ϕj:the non-spatial random effect for the overdispersion for county j
ψj:the spatially structured random effect for county j, accounting for the assumption that geographically close areas are more related than distant areas
αi:the unstructured temporal effect for year i
γi:the temporally structured effect for year i

More generally, we may expand this model to allow for an interaction between space and time, δi,j, which would explain differences in the time trend of JE risk over the six years for different counties.

All random effects were modeled, and default minimally informative hyperpriors were set [29]. Firstly, the spatial effect was modeled as the first order intrinsic Gaussian Markov random field (IGMRF) [30]. The IGMRF accounts for spatial autocorrelation by assuming that the conditional distribution in region i depends on the neighboring regions, j. Secondly, the temporally structured effect was modeled dynamically through a time neighboring structure. Finally, the unstructured spatial and temporal effects were both specified by Gaussian models with a mean of zero. The Gamma (1, 0.0005) was chosen as the prior for the precision of the above Gaussian random effects.

2.4.2. Spatial Model with Contextual Risk Factors in 2006

To examine the association between the JE relative risks and contextual risk factors, data in 2006 was used, due to the large number of cases in that year. As with [13,14], analysis was restricted to endemic months. We restricted to July, August and September, because most of the JE cases in southwest China occurred in these 3 months. The observed number of JE cases in county j in month k, yjk, was modeled as a function of the county-month-specific expected cases, ejk, and the relative risk of disease, θjk, with θjk allowing for the contextual factors:

log θ j k = μ + ϕ j + ψ j + β 1 × x j ( k l 1 ) , 1 + β 2 × x j ( k l 1 ) , 2 + β 3 × x j , 3 + β 4 × x j , 4
where μ, ϕj and ψj are the same as the spatiotemporal model; the others denote the contextual variables as follows:
xj(kl1),1:the temperature of the (kl1) month for county j
xj(kl2),2:the precipitation of the (kl2) month for county j
xj,3:the proportion of rural population for county j
xj,4:the pig-to-human ratio for county j

Here, l1 and l2 represent the lagged effects for the temperature and the precipitation, respectively, and they varied from 1 to 3 months. Moreover, similar to [13], the monthly mean minimum and maximum temperatures were modeled separately, as the two temperatures are highly correlated.

The random effects were exactly the same as the spatiotemporal model. The default highly dispersed Gaussian distribution with the mean equal to 0 and the variance equal to 1,000 was specified for all the fixed effects, β1 to β4.

Two groups of models were fitted using a forward stepwise procedure with an increasing complexity. The first group only included climatic variables to choose the best time lags, while the second model further included the proportion of rural population and the pig-to-human ratio, in order to investigate the association. A Berkson-type measurement error model was used for the pig-to-human ratio. Berkson-type error occurs in biological or epidemiological studies, where averages of exposures in areas are assigned to individuals living or working close-by [31]. Roughly speaking, the method was used to deal with the multiple resolution of the exposure variables. The pig-to-human ratio was observed at a coarser level [32,33], the prefectural level, as opposed to the county level. The prefecture-city level is an administrative division below a province and above a county in China's administrative structure. All covariates were standardized to eliminate the scale effect and to improve the accuracy of parameter estimates.

2.4.3. Computation and Model Choice

All models were computed using integrated nested Laplace approximations (INLA), a method for approximate Bayesian inference within latent Gaussian models. INLA outperforms traditional Markov chain Monte Carlo (MCMC) in terms of computational time, while keeping very precise results [34]. Analyses can be performed using R-INLA, available at “ http://www.r-inla.org/”. Competing models were compared and selected using deviance information criterion (DIC), and a smaller DIC indicates a better trade-off between the model fit and complexity [35].

3. Results

3.1. Descriptive Results

A total of 17,007 JE cases were reported in southwest China from 2004 to 2009. There were 10,489 male and 6,518 female cases. Children under five years old accounted for 48%, 5- to 15-year-old school children for 45.8% and those > 15 years old for 6.2%.

Figure 1 shows a time series of JE cases, a significant seasonal peak with 83.3% of cases occurring in July–August, and 93.8% in July–September. Furthermore, the annual incidence peaked in 2006.

Figure 1. Temporal distribution patterns of Japanese encephalitis (JE) cases in southwest China from 2004 to 2009. (A) The figure shows the epidemic curve of monthly JE cases; (B) the seasonal epidemic patterns of the JE distribution. The bottom and top of the box indicates the lower quartile (P25) and the upper quartile (P75), respectively; the line in the middle of the box represents the median value; whiskers represents 1.5 times the height of the box.
Figure 1. Temporal distribution patterns of Japanese encephalitis (JE) cases in southwest China from 2004 to 2009. (A) The figure shows the epidemic curve of monthly JE cases; (B) the seasonal epidemic patterns of the JE distribution. The bottom and top of the box indicates the lower quartile (P25) and the upper quartile (P75), respectively; the line in the middle of the box represents the median value; whiskers represents 1.5 times the height of the box.
Ijerph 11 04201 g001 1024

Figure 2 demonstrates annualized incidences of JE maps from 2004 to 2009, while Table 1 presents spatial autocorrelation analyses of JE for each year. They show that Moran's I was significant for each year, implying that the distribution of JE was spatially autocorrelated in southwest China. Moreover, Moran's I of the first three years is significantly higher than those of the subsequent three years, implying the autocorrelation became smaller later. This is also evident from Figure 2.

Figure 2. Annualized incidence of JE maps in southwest China from 2004 to 2009. The figure shows the crude annual incidence (number of cases per 100,000) in the following years: (A) 2004; (B) 2005; (C) 2006; (D) 2007; (E) 2008; (F) 2009.
Figure 2. Annualized incidence of JE maps in southwest China from 2004 to 2009. The figure shows the crude annual incidence (number of cases per 100,000) in the following years: (A) 2004; (B) 2005; (C) 2006; (D) 2007; (E) 2008; (F) 2009.
Ijerph 11 04201 g002 1024
Table 1. Spatial autocorrelation analyses for the annualized incidence of JE in southwest China from 2004 to 2009.
Table 1. Spatial autocorrelation analyses for the annualized incidence of JE in southwest China from 2004 to 2009.
YearMoran's Ip-valuePattern
20040.51960.001Clustered
20050.50140.001Clustered
20060.51460.001Clustered
20070.35550.001Clustered
20080.29450.001Clustered
20090.38740.001Clustered

Figure 3 presents the following county-level information in 2006: the crude incidence rates of JE in August, the proportion of the rural population, the pig-to-human ratio, the monthly mean minimum temperature and precipitation, with month lags varying from one to three to August. These graphs show several features. Firstly, the JE cases generally concentrated in the east area. Secondly, there is some degree of similarity between the spatial pattern of JE cases and those of the mean minimum temperature, the proportion of the rural population and the pig-to-human ratio. Thirdly, the geographical pattern for the precipitation did not show a clear correlation with JE cases. Finally, the spatial patterns of the precipitation varied from June to August, while the geographical pattern of the mean minimum temperature was preserved with continuing change in the magnitude. Here, to preserve space, we just reported the monthly mean minimum temperature. Similar results could be found for the monthly mean maximum temperature.

Table 2 shows the Spearman correlation coefficients between monthly incidences (from July to September) and the three climatic variables with different time lags in 2006. The highest positive correlation constantly occurred with a two-month lag for the precipitation, whereas the correlation for the two temperatures does not show a clear pattern.

3.2. Bayesian Hierarchical Spatial and Spatiotemporal Modeling

3.2.1. Results of Spatiotemporal Model for County-Specific Relative Risks from 2004–2009

Two different spatiotemporal models were fitted, with and without a space-time interaction. The model with the interaction had a greater fit (DIC = 26,281.63) than the model without interaction (DIC = 28,529.48). Figure 4A presents the posterior mean of county-specific excessive relative risks, the sum of spatially structured random and non-spatial random effects. This reflects the excessive spatial trend of disease risks for the 438 counties after allowing for sparse counts and correlation effects, and it is evident that high risk areas are located in the eastern and that some are far southern areas. Figure 4B shows the excessive temporal trend posterior mean of year-specific excessive relative risks and the sum of temporally structured and unstructured effects. Overall, JE presented a declining temporal trend in southwest China, with a peak in 2006.

Figure 3. County-level information in southwest China in 2006. The figure shows the following information in 2006 at the county level: (A) the crude monthly incidence in August (number of cases per 100,000); (B) the proportion of rural population; (C) the pig density (the number of pigs per 1,00 people); (D) the monthly mean minimum temperature in July; (E) the monthly mean minimum temperature in June; (F) the monthly mean minimum temperature in May; (G) the precipitation in July; (H) the precipitation in June; and (I) the precipitation in May.
Figure 3. County-level information in southwest China in 2006. The figure shows the following information in 2006 at the county level: (A) the crude monthly incidence in August (number of cases per 100,000); (B) the proportion of rural population; (C) the pig density (the number of pigs per 1,00 people); (D) the monthly mean minimum temperature in July; (E) the monthly mean minimum temperature in June; (F) the monthly mean minimum temperature in May; (G) the precipitation in July; (H) the precipitation in June; and (I) the precipitation in May.
Ijerph 11 04201 g003 1024
Table 2. Spearman correlation coefficients between monthly incidences and climatic variables with different time lags, from July to September in 2006.
Table 2. Spearman correlation coefficients between monthly incidences and climatic variables with different time lags, from July to September in 2006.
Month of incidencePrecipitationMonthly mean minimum temperatureMonthly mean maximum temperature



1-month lag2-month lag3-month lag1-month lag2-month lag3-month lag1-month lag2-month lag3-month lag
July0.233 **0.432 **0.423 **0.331 **0.335 **0.332 **0.258 **0.268 **0.131 **
August0.139 **0.237 **0.187**0.258 **0.259 **0.245 **0.221**0.229 **0.101 **
September0.0630.087 *−0.038−0.0670.074−0.077−0.051−0.055−0.063

**The difference between the correlation coefficient and zero is statistically significant, p < 0.05;*the difference between the correlation coefficient and zero is statistically significant, p < 0.1.

Figure 4. Spatiotemporal trend in southwest China from 2004 to 2009. (A) The posterior mean of county-specific excessive relative risks and the sum of spatially and non-spatial random effects; (B) the excessive temporal trend, posterior mean of year-specific excessive relative risks and the sum of temporally structured and unstructured effects.
Figure 4. Spatiotemporal trend in southwest China from 2004 to 2009. (A) The posterior mean of county-specific excessive relative risks and the sum of spatially and non-spatial random effects; (B) the excessive temporal trend, posterior mean of year-specific excessive relative risks and the sum of temporally structured and unstructured effects.
Ijerph 11 04201 g004 1024

3.2.2. Spatial Model with Contextual Risk Factors in 2006

First, DICs were used to select the best time lags for meteorological variables. Table 3 reports the DICs of models with different month lags for rainfall and temperature. The first column “precipitation” shows the DICs of models with rainfall as the only covariate, and the best fit model is the two-month lag precipitation, with DIC being 8614.96. The result is consistent with the Spearman's correlation analysis. Together with the two-month lag precipitation, different month lags for the monthly mean minimum temperature were modeled. The result is reported in the second column, and the one-month lag presents the best fit, with DIC equal to 8503.39. The monthly mean maximum temperature shows a similar result with the best fit with DIC equal to 8508.47. The minimum temperature shows a slightly better fit than the maximum temperature by the comparison of DICs.

Table 4 presents estimated parameters of models with spatial random effects and all contextual factors, including agricultural variables and the selected time-lagged climatic variables. Monthly mean minimum and maximum temperatures were modeled separately. The two models show a pronounced increase in model fit over the models just with climatic covariates, with the DIC decreasing by about 34% from a comparison of Table 3 and Table 4. Besides, as before, the model with monthly mean minimum temperature has a greater fit (DIC = 5470.98) compared to the model with monthly mean maximum temperature (DIC = 5702.9). These contextual variables all show a statistically significant positive association with JE, as their 95% credible intervals are over zero. Although the monthly mean minimum and maximum temperatures are not directly comparable as they are in different models, the former presents a greater regression coefficient over that of the latter, 1.097 (C.I.= 0.979, 1.222) versus 0.483 (C.I. = 0.408, 0.559). Since all contextual variables are standardized, this may suggest that the variation of the minimum temperature has a greater association with JE cases compared to the maximum temperature. This is consistent with the result of DICs. Furthermore, in both models, the regression coefficients of temperature are greater than the precipitation. For example, in the model with the minimum temperature, the coefficient for the monthly mean minimum temperature is 1.097, while the coefficient for the precipitation is 0.261. Since all independent variables were standardized, this may imply that the variation of temperature is more important for JE variation than precipitation. Moreover, examining Figure 3 shall lead to the same conclusion, as temperature shows a remarkably greater geographic similarity with JE than the precipitation. As expected, coefficients for the proportion of rural population and the pig-to-human ratio are both statistically positively significant.

Table 3. Deviance information criteria (DICs) of models with different month lags for the precipitation and temperature in 2006.
Table 3. Deviance information criteria (DICs) of models with different month lags for the precipitation and temperature in 2006.
Time lagPrecipitationMonthly mean minimum temperature +2-month lag precipitationMonthly mean maximum temperature +2-month lag precipitation
1-month8,707.218,503.398,508.47
2-month8,614.968,614.068,554.03
3-month8,618.338,582.918,616.4
Table 4. Parameter estimates of models with agriculture related variables and the selected time-lagged climatic variables, in addition to the spatial random effects.
Table 4. Parameter estimates of models with agriculture related variables and the selected time-lagged climatic variables, in addition to the spatial random effects.
VariableModel with monthly mean minimum temperatureModel with monthly mean maximum temperature


meanSD95%Credible intervalmeanSD95%Credible interval
temperature1.0970.060.9791.2230.4830.0390.4080.559
precipitation0.2610.0240.2150.3070.2350.0230.190.28
rural0.350.0910.1720.5290.3110.0840.150.48
pig0.2420.0870.0670.410.3330.0810.1770.485
intercept−1.0170.101−1.22−0.819−0.7970.103−1−0.566


DIC5470.985702.9

4. Discussion

Japanese encephalitis still remains one of the major national arboviral diseases in China [5] and is mainly located in southwest China [36], which is our focus. Unlike existing studies in a single site, our study used data from 438 counties in southwest China. Besides, our analysis was relatively comprehensive by including meteorological and agricultural variables.

There are two main findings in the current work. First, spatiotemporal analysis has identified the east part of southwest China as a high risk area and that area may be connected to agricultural and climatic factors. Based on the ecological regression analysis, the high risk area is associated with the high pig-to-human ratio, a high proportion of rural population and weather change. An unpublished thesis [37] had analyzed the association between agricultural factors and JE in the southwest of China at the prefecture-city level, an administrative division below a province and above a county in China's administrative structure, but found no statistically significant factors. Including a comprehensive set of variables and the county-level analysis are the unique characteristics of this study. The model-based approach confirms the association between JE and contextual risk factors. Since southwest China is a less-developed area compared with the other parts in China, it has limited funds for JE control. These identified high risk areas should get the priority in allocating health resources to promote the existing surveillance, infrastructure redevelopment and in-house workforce training. Two strategies may be implemented for pig-rearing, pig immunization and the separation of pig rearing from human settlements. Besides, monitoring the infection rate in pigs may supply information for possible outbreaks.

Second, the time-lagged association between JE incidences and temperature and precipitation are quantified with an adjustment for other agricultural variables. It is important to study the impact of climate on JE transmission, because global warming might change the patterns of temperature and precipitation, which may influence the development of mosquitoes and the virus.

Minimum and maximum temperatures are found to be positively associated with JE incidences, which is biologically plausible. Temperature can affect three aspects of JE transmission: the survival and reproduction rates of mosquitoes, the biting rate of mosquitoes and the development of JE virus within mosquitoes. Thus, higher temperatures, within limits, can lead to the quicker development of larvae, shorter times between blood meals and faster incubation times for viral infections within mosquitoes. As a result, higher temperatures cause mosquito populations to reach a higher density faster and to be maintained for a longer period, thereby increasing the chance for viral transmission [38].

The minimum temperature results in a better model fit and a greater regression coefficient compared to the maximum temperature, implying that the minimum temperature presents a greater effect on JE transmission compared with the maximum temperature. Minimum temperatures might play a more important role for larvae survival in cold conditions. Similar results were reported for other vector-borne diseases [39,40,41].

Precipitation is found to be positively associated with JE incidences. Precipitation would bring a high relative humidity, which may increase the propagation and development of mosquitoes. High rainfall can also provide water support for the development of mosquitoes at the larval and pupal stages.

The one-month lag for temperature is consistent with [14,18], and the two-month lag for precipitation is consistent with [13,18]. There are several periods to be considered for the lagged effect, such as the time for mosquitoes to develop, the development period of JE virus within the mosquito and the incubation period of the virus within the human body. The lag times are distinct between temperature and rainfall, with a shorter lag for temperature. This might result from the difference regarding the biological mechanisms. The temperature can contribute to the transmission of JE through its effect on the development of mosquitoes and the viruses within them, whereas rainfall may just affect the mosquitoes. Hence, temperature may present a faster and stronger effect. The stronger effect can be observed from Figure 3 and the regression coefficient in Table 4, in which the variable temperature factors demonstrate a larger regression coefficient than the precipitation. The interaction between temperature and rainfall was not included in the model for two reasons. First, it is not unreasonable to assume there is no such an interaction, which is the usual assumption from most existing studies for vector-borne diseases [17,42]. Second, the preliminary exploratory analysis did not show a significant interaction.

The relationship between climate and JE may be highly dependent on local environmental factors, and it is not always possible to extrapolate to a broader scale [15]. Therefore, existing results from other sites may not be appropriate for southwest China. Due to the severe consequences of JE, surveillance for the early warning of epidemics based on meteorological factors in combination with agricultural variables may be of paramount importance in southwest China. Our findings give a better understanding for the JE epidemiology in this area. When designing a JE early warning system in southwest China, health departments should consider the lag patterns for meteorological factors, as well as the strongest impact from minimum temperatures.

The limitations of this study should be acknowledged. Firstly, some factors were not included in the ecological regression, such as the population immunity (including vaccination), the area of the irrigated land, the socioeconomic status and the control measures for mosquitoes. Unfortunately, those data were unavailable over the 438 counties in 2006. The hierarchical Bayesian spatial method was used as an effort to better estimate the parameters for the observed variables by allowing for unmeasured covariates [43]. However, the potential bias cannot be eliminated. Secondly, only the data in 2006 was used to investigate the association, as it is very difficult in practice to collect all variables over the 438 counties for the six years.

5. Conclusion

Climate change, especially global warming, has already brought and will continue to bring about challenges for infectious disease control [44]. In our work, a county-level analysis identified the east part of southwest China as the high risk area, which may be connected to agricultural and climatic factors, and therefore, the routine surveillance and the allocation of health resources should be given more attention in this area. Pig immunization and the separation of pig rearing from human settlements may also assist in limiting the transmission. Moreover, the monthly mean minimum and maximum temperatures and monthly rainfall showed lagged correlations with JE. Temperature demonstrated a greater correlation than rainfall, and the minimum temperature presented a greater model fit than the maximum temperature. The climatic factors might be considered as possible predictors of JE in southwest China.

Acknowledgments

This study was supported by the Natural Science Foundation of China (No. 30571618) and China Scholarship Council ( http://en.csc.edu.cn/).

Author Contributions

Xing Zhao performed the statistical analysis and drafted the manuscript. Mingqin Cao, Fei Chen, Heng Fan and Hai-Huan Feng cleared the data. Mingqin Cao conducted some preliminary exploration. Zijian Feng provided the original data. Xiaosong Li conceived of the project concept. Xiao-Hua Zhou gave technical assists. All of the authors have read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Solomon, T. Control of Japanese encephalitis—Within our grasp? N. Engl. J. Med. 2006, 355, 869–871. [Google Scholar]
  2. Erlanger, T.E.; Weiss, S.; Keiser, J.; Utzinger, J.; Wiedenmayer, K. Past, present, and future of Japanese encephalitis. Emerg. Infect. Dis. 2009, 15, 1–7. [Google Scholar]
  3. Diagana, M.; Preux, P.M.; Dumas, M. Japanese encephalitis revisited. J. Neurol. Sci. 2007, 262, 165–170. [Google Scholar]
  4. Campbell, G.L.; Hills, S.L.; Fischer, M.; Jacobson, J.A.; Hoke, C.H.; Hombach, J.M.; Marfin, A.A.; Solomon, T.; Tsai, T.F.; Tsu, V.D.; et al. Estimated global incidence of Japanese encephalitis: A systematic review. Bull. World Health Organ. 2011, 89. [Google Scholar]
  5. Gao, X.; Nasci, R.; Liang, G. The neglected arboviral infections in mainland China. PLoS Negl. Trop. Dis. 2010, 4. [Google Scholar] [CrossRef]
  6. Tao, S. Epidemiological surveillance and prophylaxis of Japanese encephalitis. Chin. J. Vaccines Immun. 2002, 4, 25. [Google Scholar]
  7. Wang, L.; Zhang, W.; Ding, F.; Sun, H.; Yu, S.; Li, S.; Chen, Z.; Huang, L.; Li, C. Deaths associated with Japanese encephalitis, China, 2005–2010. Clin. Infect. Dis. 2013, 56, 752. [Google Scholar] [CrossRef]
  8. Gingrich, J.B.; Nisalak, A.; Latendresse, J.R.; Sattabongkot, J.; Hoke, C.H.; Pomsdhit, J.; Chantalakana, C.; Satayaphanta, C.; Uechiewcharnkit, K.; Innis, B.L. Japanese encephalitis virus in Bangkok: Factors influencing vector infections in three suburban communities. J. Med. Entomol. 1992, 29, 436–444. [Google Scholar]
  9. Rohr, J.R.; Dobson, A.P.; Johnson, P.T.J.; Kilpatrick, A.M.; Paull, S.H.; Raffel, T.R.; Ruiz-Moreno, D.; Thomas, M.B. Frontiers in climate change-disease research. Trends Ecol. Evol. 2011, 26, 270–277. [Google Scholar]
  10. Keiser, J.; Maltese, M.F.; Erlanger, T.E.; Bos, R.; Tanner, M.; Singer, B.H.; Utzinger, J. Effect of irrigated rice agriculture on Japanese encephalitis, including challenges and opportunities for integrated vector management. Acta Trop. 2005, 95, 40–57. [Google Scholar]
  11. Cao, M.; Feng, Z.; Zhang, J.; Ma, J.; Li, X. Contextual risk factors for regional distribution of Japanese encephalitis in the People's Republic of China. Trop. Med. Int. Health 2010, 15, 918–923. [Google Scholar]
  12. Le Flohic, G.; Porphyre, V.; Barbazan, P.; Gonzalez, J.P. Review of climate, landscape, and viral genetics as drivers of the Japanese encephalitis virus ecology. PLoS Negl. Trop. Dis. 2013, 7. [Google Scholar] [CrossRef]
  13. Bi, P.; Zhang, Y.; Parton, K.A. Weather variables and Japanese encephalitis in the metropolitan area of Jinan city, China. J. Infect. 2007, 55, 551–556. [Google Scholar]
  14. Bi, P.; Tong, S.; Donald, K.; Parton, K.A.; Ni, J. Climate variability and transmission of Japanese encephalitis in eastern China. Vector Borne Zoonotic Dis. 2003, 3, 111–115. [Google Scholar]
  15. Lin, H.; Yang, L.; Liu, Q.; Wang, T.; Hossain, S.R.; Ho, S.C.; Tian, L. Time series analysis of Japanese encephalitis and weather in Linyi City, China. Int. J. Public Health 2012, 57, 289–296. [Google Scholar]
  16. Wang, H.; Li, Y.; Liang, X.; Liang, G. Japanese encephalitis in mainland China. Jpn. J. Infect. Dis. 2009, 62, 331–336. [Google Scholar]
  17. Bai, L.; Morton, L.C.; Liu, Q. Climate change and mosquito-borne diseases in China: A review. Glob. Health 2013, 9. [Google Scholar] [CrossRef]
  18. Hsu, S.M.; Yen, A.M.F.; Chen, T.H.H. The impact of climate on Japanese encephalitis. Epidemiol. Infect. 2008, 136, 980–987. [Google Scholar]
  19. Impoinvil, D.E.; Solomon, T.; Schluter, W.W.; Rayamajhi, A.; Bichha, R.P.; Shakya, G.; Caminade, C.; Baylis, M. The spatial heterogeneity between Japanese encephalitis incidence distribution and environmental variables in Nepal. PLoS One 2011, 6. [Google Scholar] [CrossRef]
  20. Shen, Y.; Feng, M.; Hongzheng, Z.; Feng, G. Interpolation methods of China daily orecipitation data. J. Appl. Meterol. Sci. 2010, 21, 279–286. [Google Scholar]
  21. Zhang, Q.; Ruan, X.; Xiong, A. Establishment and assessment of the grid air temperature data sets in China for the past 57 years. J. Appl. Meterol. Sci. 2009, 20, 385–393. [Google Scholar]
  22. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2013. [Google Scholar]
  23. Assunção, R.M.; Reis, E.A. A new proposal to adjust Moran's I for population density. Stat. Med. 1999, 18, 2147–2162. [Google Scholar]
  24. Besag, J.; York, J.; Mollié, A. Bayesian image restoration, with two applications in spatial statistics. Ann. Inst. Stat. Math. 1991, 43, 1–20. [Google Scholar]
  25. Knorr-Held, L. Bayesian modelling of inseparable space-time variation in disease risk. Stat. Med. 2000, 19, 2555–2567. [Google Scholar]
  26. Waller, L.A.; Gotway, C.A. Applied Spatial Statistics for Public Health Data; Wiley-Interscience: Hoboken, NJ, USA, 2004. [Google Scholar]
  27. Elliott, P.; Wartenberg, D. Spatial epidemiology: Current approaches and future challenges. Environ. Health Perspect. 2004, 112, 998–1006. [Google Scholar]
  28. Lambert, D. Zero-inflated Poisson regression, with an application to defects in manufacturing. Technometrics 1992, 34, 1–14. [Google Scholar]
  29. Blangiardo, M.; Cameletti, M.; Baio, G.; Rue, H. Spatial and spatio-temporal models with R-INLA. Spat. Spatiotemporal Epidemiol. 2013, 4, 33–49. [Google Scholar]
  30. Besag, J.; Kooperberg, C. On conditional and intrinsic autoregressions. Biometrika 1995, 82, 733–746. [Google Scholar]
  31. Carroll, R.J.; Ruppert, D.; Stefanski, L.A.; Crainiceanu, CM. Measurement Error in Nonlinear Models: A Modern Perspective, 2nd ed.; Chapman and Hall/CRC: London, UK, 2006. [Google Scholar]
  32. Berkson, J. Are there two regressions? J. Am. Stat. Assoc. 1950, 45, 164–180. [Google Scholar]
  33. Muff, S.; Riebler, A.; Rue, H.; Saner, P.; Held, L. Bayesian Analysis of Measurement Error Models Using INLA. Available online: http://arxiv.org/pdf/1302.3065v2.pdf (accessed on 31 November 2013).
  34. Rue, H.; Martino, S.; Chopin, N. Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations. J. R. Stat. Soc. Ser. B (Stat. Methodol.) 2009, 71, 319–392. [Google Scholar]
  35. Spiegelhalter, D.J.; Best, N.G.; Carlin, B.P.; van Der Linde, A. Bayesian measures of model complexity and fit. J. R. Stat. Soc. Ser. B (Stat. Methodol.) 2002, 64, 583–639. [Google Scholar]
  36. Wang, L.Y.; Zhang, W.Y.; Ding, F.; Hu, W.B.; Soares Magalhaes, R.J.; Sun, H.L.; Li, Y.X.; Zou, W.; Wang, Y.; Liu, Q.Y.; et al. Spatiotemporal patterns of Japanese encephalitis in China, 2002–2010. PLoS Negl. Trop. Dis. 2013, 7. [Google Scholar] [CrossRef]
  37. Cao, M. A Study on the Multiple Membership Multiple Classification Models and Disease Mapping for Analyzing Spatial Dependence and Heterogeneity of Regional Distribution of Japanese Encephalitis in Southwestern China. PhD Thesis, Sichuan University, Sichuan, China, 2009. [Google Scholar]
  38. Russell, R.C. Ross River virus: Ecology and distribution. Annu. Rev. Entomol. 2002, 47, 1–31. [Google Scholar]
  39. Hui, F.M.; Xu, B.; Chen, Z.W.; Cheng, X.; Liang, L.; Huang, H.B.; Fang, L.Q.; Yang, H.; Zhou, H.N.; Yang, H.L.; et al. Spatio-temporal distribution of malaria in Yunnan Province, China. Am J Trop Med Hyg 2009, 81, 503–509. [Google Scholar]
  40. Bi, P.; Tong, S.; Donald, K.; Parton, K.A.; Ni, J. Climatic variables and transmission of malaria: A 12-year data analysis in Shuchen County, China. Public Health Rep. 2003, 118, 65–71. [Google Scholar]
  41. Lu, L.; Lin, H.; Tian, L.; Yang, W.; Sun, J.; Liu, Q. Time series analysis of dengue fever and weather in Guangzhou, China. BMC Public Health 2009, 9. [Google Scholar] [CrossRef]
  42. Mabaso, M.L.H.; Ndlovu, N.C. Critical review of research literature on climate-driven malaria epidemics in sub-Saharan Africa. Public Health 2012, 126, 909–919. [Google Scholar]
  43. Clayton, D.G.; Bernardinelli, L.; Montomoli, C. Spatial correlation in ecological analysis. Int. J. Epidemiol. 1993, 22, 1193–1202. [Google Scholar]
  44. Lafferty, K.D. The ecology of climate change and infectious diseases. Ecology 2009, 90, 888–900. [Google Scholar]

Share and Cite

MDPI and ACS Style

Zhao, X.; Cao, M.; Feng, H.-H.; Fan, H.; Chen, F.; Feng, Z.; Li, X.; Zhou, X.-H. Japanese Encephalitis Risk and Contextual Risk Factors in Southwest China: A Bayesian Hierarchical Spatial and Spatiotemporal Analysis. Int. J. Environ. Res. Public Health 2014, 11, 4201-4217. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph110404201

AMA Style

Zhao X, Cao M, Feng H-H, Fan H, Chen F, Feng Z, Li X, Zhou X-H. Japanese Encephalitis Risk and Contextual Risk Factors in Southwest China: A Bayesian Hierarchical Spatial and Spatiotemporal Analysis. International Journal of Environmental Research and Public Health. 2014; 11(4):4201-4217. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph110404201

Chicago/Turabian Style

Zhao, Xing, Mingqin Cao, Hai-Huan Feng, Heng Fan, Fei Chen, Zijian Feng, Xiaosong Li, and Xiao-Hua Zhou. 2014. "Japanese Encephalitis Risk and Contextual Risk Factors in Southwest China: A Bayesian Hierarchical Spatial and Spatiotemporal Analysis" International Journal of Environmental Research and Public Health 11, no. 4: 4201-4217. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph110404201

Article Metrics

Back to TopTop