Zhoupeng Ren

Predicting malaria vector distribution under climate change scenarios in China: Challenges for malaria elimination

Projecting the distribution of malaria vectors under climate change is essential for planning integrated vector control activities for sustaining elimination and preventing reintroduction of malaria. In China, however, little knowledge exists on the possible effects of climate change on malaria vectors. Here we assess the potential impact of climate change on four dominant malaria vectors (An. dirus, An. minimus, An. lesteri and An. sinensis) using species distribution models for two future decades: the 2030 s and the 2050 s. Simulation-based estimates suggest that the environmentally suitable area (ESA) for An. dirus and An. minimus would increase by an average of 49% and 16%, respectively, under all three scenarios for the 2030 s, but decrease by 11% and 16%, respectively in the 2050 s. By contrast, an increase of 36% and 11%, respectively, in ESA of An. lesteri and An. sinensis, was estimated under medium stabilizing (RCP4.5) and very heavy (RCP8.5) emission scenarios. in the 2050 s. In total, we predict a substantial net increase in the population exposed to the four dominant malaria vectors in the decades of the 2030 s and 2050 s, considering land use changes and urbanization simultaneously. Strategies to achieve and sustain malaria elimination in China will need to account for these potential changes in vector distributions and receptivity.

Spatio-temporal analysis of malaria vectors in national malaria surveillance sites in China

BackgroundTo reveal the spatio-temporal distribution of malaria vectors in the national malaria surveillance sites from 2005 to 2010 and provide reference for the current National Malaria Elimination Programme (NMEP) in China.MethodsA 6-year longitudinal surveillance on density of malaria vectors was carried out in the 62 national malaria surveillance sites. The spatial and temporal analyses of the four primary vectors distribution were conducted by the methods of kernel k-means and the cluster distribution of the most widely distribution vector of An.sinensis was identified using the empirical mode decomposition (EMD).ResultsTotally 4 species of Anopheles mosquitoes including An.sinensis, An.lesteri, An.dirus and An.minimus were captured with significant difference of distribution as well as density. An. sinensis was the most widely distributed, accounting for 96.25% of all collections, and its distribution was divided into three different clusters with a significant increase of density observed in the second cluster which located mostly in the central parts of China.ConclusionThis study first described the spatio-temporal distribution of malaria vectors based on the nationwide surveillance during 2005–2010, which served as a baseline for the ongoing national malaria elimination program.

Geographically weighted regression-based determinants of malaria incidences in northern China

Geographically weighted regression (GWR) is an important local method to explore spatial non-stationarity in data relationships. It has been repeatedly used to examine spatially varying relationships between epidemic diseases and predictors. Malaria, a serious parasitic disease around the world, shows spatial clustering in areas at risk. In this article, we used GWR to explore the local determinants of malaria incidences over a 7-year period in northern China, a typical mid-latitude, high-risk malaria area. Normalized difference vegetation index (NDVI), land surface temperature (LST), temperature difference, elevation, water density index (WDI) and gross domestic product (GDP) were selected as predictors. Results showed that both positively and negatively local effects on malaria incidences appeared for all predictors except for WDI and GDP. The GWR model calibrations successfully depicted spatial variations in the effect sizes and levels of parameters, and also showed substantially improvements in terms of goodness of fits in contrast to the corresponding non-spatial ordinary least squares (OLS) model fits. For example, the diagnostic information of the OLS fit for the 7-year average case is R2 = 0.243 and AICc = 837.99, while significant improvement has been made by the GWR calibration with R2 = 0.800 and AICc = 618.54.

Spatial-Temporal Variation and Primary Ecological Drivers of Anopheles sinensis Human Biting Rates in Malaria Epidemic- Prone Regions of China

Background Robust malaria vector surveillance is essential for optimally selecting and targeting vector control measures. Sixty-two vector surveillance sites were established between 2005 and 2008 by the national malaria surveillance program in China to measure Anopheles sinensis human biting rates. Using these data to determine the primary ecological drivers of malaria vector human biting rates in malaria epidemic-prone regions of China will allow better targeting of vector control resources in space and time as the country aims to eliminate malaria. Methods We analyzed data from 62 malaria surveillance sentinel sites from 2005 to 2008. Linear mixed effects models were used to identify the primary ecological drivers for Anopheles sinensis human biting rates as well as to explore the spatial-temporal variation of relevant factors at surveillance sites throughout China. Results Minimum semimonthly temperature (β = 2.99; 95% confidence interval (CI) 2.07- 3.92), enhanced vegetation index (β =1.07; 95% CI 0.11–2.03), and paddy index (the percentage of rice paddy field in the total cultivated land area of each site) (β = 0.86; 95% CI 0.17–1.56) were associated with greater An. Sinensis human biting rates, while increasing distance to the nearest river was associated with lower An. Sinensis human biting rates (β = −1.47; 95% CI −2.88, −0.06). The temporal variation (σt02=1.35) in biting rates was much larger than the spatial variation (σs02=0.83), with 19.3% of temporal variation attributable to differences in minimum temperature and enhanced vegetation index and 16.9% of spatial variance due to distance to the nearest river and the paddy index. Discussion Substantial spatial-temporal variation in An. Sinensis human biting rates exists in malaria epidemic-prone regions of China, with minimum temperature and enhanced vegetation index accounting for the greatest proportion of temporal variation and distance to nearest river and paddy index accounting for the greatest proportion of spatial variation amongst observed ecological drivers. Conclusions Targeted vector control measures based on these findings can support the ongoing malaria elimination efforts in China more effectively.

Spatial patterns of the congenital heart disease prevalence among 0- to 14-year-old children in Sichuan Basin, P. R China, from 2004 to 2009

BackgroundCongenital heart disease (CHD) is the most common type of major birth defects in Sichuan, the most populous province in China. The detailed etiology of CHD is unknown but some environmental factors are suspected as the cause of this disease. However, the geographical variations in CHD prevalence would be highly valuable in providing a clue on the role of the environment in CHD etiology. Here, we investigate the spatial patterns and geographic differences in CHD prevalence among 0- to 14-year-old children, discuss the possible environmental risk factors that might be associated with CHD prevalence in Sichuan Basin from 2004 to 2009.MethodsThe hierarchical Bayesian model was used to estimate CHD prevalence at the township level. Spatial autocorrelation statistics were performed, and a hot-spot analysis with different distance thresholds was used to identify the spatial pattern of CHD prevalence. Distribution and clustering maps were drawn using geographic information system tools.ResultsCHD prevalence was significantly clustered in Sichuan Basin in different spatial scale. Typical hot/cold clusters were identified, and possible CHD causes were discussed. The association between selected hypothetical environmental factors of maternal exposure and CHD prevalence was evaluated.ConclusionsThe largest hot-spot clustering phenomena and the CHD prevalence clustering trend among 0- to 14-year-old children in the study area showed a plausibly close similarity with those observed in the Tuojiang River Basin. The high ecological risk of heavy metal(Cd, As, and Pb)sediments in the middle and lower streams of the Tuojiang River watershed and ammonia–nitrogen pollution may have contribution to the high prevalence of CHD in this area.

A spatial and temporal analysis of Japanese encephalitis in mainland China, 1963-1975: a period without Japanese encephalitis vaccination.

More than a million Japanese encephalitis (JE) cases occurred in mainland China from the 1960s to 1970s without vaccine interventions. The aim of this study is to analyze the spatial and temporal pattern of JE cases reported in mainland China from 1965 to 1973 in the absence of JE vaccination, and to discuss the impacts of climatic and geographical factors on JE during that period. Thus, the data of reported JE cases at provincial level and monthly precipitation and monthly mean temperature from 1963 to 1975 in mainland China were collected. Local Indicators of Spatial Association analysis was performed to identify spatial clusters at the province level. During that period, The epidemic peaked in 1966 and 1971 and the JE incidence reached up to 20.58/100000 and 20.92/100000, respectively. The endemic regions can be divided into three classes including high, medium, and low prevalence regions. Through spatial cluster analysis, JE epidemic hot spots were identified; most were located in the Yangtze River Plain which lies in the southeast of China. In addition, JE incidence was shown to vary among eight geomorphic units in China. Also, the JE incidence in the Loess Plateau and the North China Plain was showed to increase with the rise of temperature. Likewise, JE incidence in the Loess Plateau and the Yangtze River Plain was observed a same trend with the increase of rainfall. In conclusion, the JE cases clustered geographically during the epidemic period. Besides, the JE incidence was markedly higher on the plains than plateaus. These results may provide an insight into the epidemiological characteristics of JE in the absence of vaccine interventions and assist health authorities, both in China and potentially in Europe and Americas, in JE prevention and control strategies.

Using spatial multilevel regression analysis to assess soil type contextual effects on neural tube defects

The rate of neural tube defects (NTDs) in Shanxi Province is the highest world widely. Both human and environmental factors can induce NTDs, but various studies ignored contextual effects. This research examines whether there are significant soil type contextual effects on the rate of NTDs. A spatial two-level regression model is used to quantify the magnitude of contextual effects. Spatial autocorrelated errors structure is used to control autocorrelation of residuals. The results suggest that the spatial multilevel model fit the data better than non-spatial multilevel models. Our findings indicate that there are significant soil type contextual effects on the rate of NTDs, even after taking into account of fertilizer and net income. More attentions should be focused on how characteristics of each soil type may affect the rates of NTDs in further studies, which is a relevant issue for understanding etiology of NTDs.

Maternal exposure to ambient PM 10 during pregnancy increases the risk of congenital heart defects: Evidence from machine learning models

Previous research suggested an association between maternal exposure to ambient air pollutants and risk of congenital heart defects (CHDs), though the effects of particulate matter ≤10μm in aerodynamic diameter (PM10) on CHDs are inconsistent. We used two machine learning models (i.e., random forest (RF) and gradient boosting (GB)) to investigate the non-linear effects of PM10 exposure during the critical time window, weeks 3-8 in pregnancy, on risk of CHDs. From 2009 through 2012, we carried out a population-based birth cohort study on 39,053 live-born infants in Beijing. RF and GB models were used to calculate odds ratios for CHDs associated with increase in PM10 exposure, adjusting for maternal and perinatal characteristics. Maternal exposure to PM10 was identified as the primary risk factor for CHDs in all machine learning models. We observed a clear non-linear effect of maternal exposure to PM10 on CHDs risk. Compared to 40μgm-3, the following odds ratios resulted: 1) 92μgm-3 [RF: 1.16 (95% CI: 1.06, 1.28); GB: 1.26 (95% CI: 1.17, 1.35)]; 2) 111μgm-3 [RF: 1.04 (95% CI: 0.96, 1.14); GB: 1.04 (95% CI: 0.99, 1.08)]; 3) 124μgm-3 [RF: 1.01 (95% CI: 0.94, 1.10); GB: 0.98 (95% CI: 0.93, 1.02)]; 4) 190μgm-3 [RF: 1.29 (95% CI: 1.14, 1.44); GB: 1.71 (95% CI: 1.04, 2.17)]. Overall, both machine models showed an association between maternal exposure to ambient PM10 and CHDs in Beijing, highlighting the need for non-linear methods to investigate dose-response relationships.

Co-endemicity of Pulmonary Tuberculosis and Intestinal Helminth Infection in the People's Republic of China.

Both pulmonary tuberculosis (PTB) and intestinal helminth infection (IHI) affect millions of individuals every year in China. However, the national-scale estimation of prevalence predictors and prevalence maps for these diseases, as well as co-endemic relative risk (RR) maps of both diseases’ prevalence are not well developed. There are co-endemic, high prevalence areas of both diseases, whose delimitation is essential for devising effective control strategies. Bayesian geostatistical logistic regression models including socio-economic, climatic, geographical and environmental predictors were fitted separately for active PTB and IHI based on data from the national surveys for PTB and major human parasitic diseases that were completed in 2010 and 2004, respectively. Prevalence maps and co-endemic RR maps were constructed for both diseases by means of Bayesian Kriging model and Bayesian shared component model capable of appraising the fraction of variance of spatial RRs shared by both diseases, and those specific for each one, under an assumption that there are unobserved covariates common to both diseases. Our results indicate that gross domestic product (GDP) per capita had a negative association, while rural regions, the arid and polar zones and elevation had positive association with active PTB prevalence; for the IHI prevalence, GDP per capita and distance to water bodies had a negative association, the equatorial and warm zones and the normalized difference vegetation index had a positive association. Moderate to high prevalence of active PTB and low prevalence of IHI were predicted in western regions, low to moderate prevalence of active PTB and low prevalence of IHI were predicted in north-central regions and the southeast coastal regions, and moderate to high prevalence of active PTB and high prevalence of IHI were predicted in the south-western regions. Thus, co-endemic areas of active PTB and IHI were located in the south-western regions of China, which might be determined by socio-economic factors, such as GDP per capita.

Spatial pattern and variations in the prevalence of congenital heart disease in children aged 4–18 years in the Qinghai-Tibetan Plateau

PURPOSE This study aimed to investigate the spatial distribution pattern of the prevalence of congenital heart disease (CHD) in children in Qinghai-Tibetan Plateau (QTP), a high-altitude region in China. METHODS Epidemiological data from a survey on the prevalence of CHD in Qinghai Province including 288,066 children (4-18 years) were used in this study. The prevalence and distribution pattern of CHD was determined by sex, CHD subtype, and nationality and altitude. Spatial pattern analysis using Getis-Ord Gi⁎ was used to identify the spatial distribution of CHD. Bayesian spatial binomial regression was performed to examine the relationship between the prevalence of CHD and environmental risk factors in the QTP. RESULTS The prevalence of CHD showed a significant spatial clustering pattern. The Tibetan autonomous prefecture of Yushu (average altitude > 4000 m) and the Mongolian autonomous county of Henan (average altitude > 3600 m) in Huangnan had the highest prevalence of CHD. Univariate analysis showed that with ascending altitude, the total prevalence of CHD, that in girls and boys with CHD, and that of the subtypes PDA and ASD increasing accordingly. Thus, environmental factors greatly contributed to the prevalence of CHD. CONCLUSIONS The prevalence of CHD shows significant spatial clustering pattern in the QTP. The CHD subtype prevalence clustering pattern has statistical regularity which would provide convenient clues of environmental risk factors. Our results may provide support to make strategies of CHD prevention, to reduce the incidence of CHD in high altitude regions of China.