Spatio-temporal analysis of malaria vectors in national malaria surveillance sites in China
- Ji-Xia Huang†^{1, 2, 3},
- Zhi-Gui Xia†^{4},
- Shui-Sen Zhou^{4},
- Xiao-Jun Pu^{5},
- Mao-Gui Hu^{3, 6},
- Da-Cang Huang^{7},
- Zhou-Peng Ren^{3, 6},
- Shao-Sen Zhang^{4},
- Man-ni Yang^{4},
- Duo-Quan Wang^{4}Email author and
- Jin-Feng Wang^{3, 6}Email author
https://doi.org/10.1186/s13071-015-0741-5
© Huang et al.; licensee BioMed Central. 2015
Received: 11 September 2014
Accepted: 13 February 2015
Published: 7 March 2015
Abstract
Background
To 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.
Methods
A 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).
Results
Totally 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.
Conclusion
This 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.
Keywords
Malaria vectors Spatio-temporal distribution Cluster Surveillance sitesBackground
Malaria, HIV/AIDS and tuberculosis (TB) are considered to be the three major public health problems in the world. In order to timely monitor the malaria epidemic and provide evidence for strategies development as well as evaluation of control activities, China government issued in 2005 the malaria surveillance program which was trialed out until 2010, and totally 62 townships in 18 provinces were selected as the sentinel sites where the malaria vectors were surveyed every 15 days from June to October [1]. However, the spatio-temporal information of continuous vectors surveillance has not been mined particularly for the NMEP from the perspective of vector control.
With continuous reduction of malaria incidence, a strong political commitment was made by China through issuing the National Action Plan for Malaria Elimination (2010–2020) in 2010. The goal is to eliminate locally acquired malaria by the end of 2015 except for the bordering areas in Yunnan Province, and eliminate the disease by 2020 nationwide. Though well-described relationships between the parasite, man and vector will provide an empirical basis for eradication of this disease, until recently there has been a paucity of systematic data on malaria vectors distribution in the malaria epidemic-prone regions of China. Undoubtedly, the efficiency of malaria elimination interventions will depend largely on information of the distribution of primary vectors. It is evident from a series of studies in China as well as in other countries [2-7] that the temporal and spatial distribution of malaria vectors is very dynamic. Thus, knowledge of the spatio-temporal malaria vector distribution will facilitate the malaria elimination efforts and overall mapping is needed to focus on priority areas for additional surveillance and response.
Therefore, the purpose of this study is to map the spatio-temporal distribution of malaria vectors to provide a dataset that would aid informed decision making on malaria elimination program, and to provide a basis for implementation of the elimination strategy in China.
Methods
Study areas
- 1.
Unstable endemic areas: 30 townships in 6 provinces (Hainan, Yunnan, Anhui, Hubei, Henan, Jiangsu) [5 counties/ province, 1 township/county].
- 2.
Low endemic areas: 24 townships in 8 provinces (Sichuan, Chongqing, Guizhou, Guangdong, Guangxi, Hunan, Jiangxi, Fujian) [3 counties/province, 1 township/county].
- 3.
Pre-elimination areas: 8 townships in 4 provinces (Shanghai, Zhejiang, Shandong, Liaoning) [2 counties/province, 1 township/county].
Malaria vector surveillance methods
According to WHO recommendations [8], human landing catches were made by two adult volunteers from the local population working beside the bed-net with one sleeping person. Mosquitoes coming to bite the collectors or sleeping person were detected using a flashlight, collected with glass tubes (CDC backpack aspirator: John W. Hock Co., Florida, USA) and placed in the screened pint-sized containers. Collections were conducted from 18:00 to 06:00 overnight during every 15 days from June to October, 2005–2010. Collectors worked in pairs for 6-h shifts. Mosquitoes were taken to the laboratory and killed by suffocation with chloroform vapor. They were counted and identified morphologically using taxonomic keys [9], and the density was calculated as the number of female adults per man per night.
Spatial and temporal clustering
- (a)
Initializing the centers m _{ k }(k =1,…,K) for those clusters in feature space.
- (b)Assigning each observed value φ(x _{ i })(i = 1, …, N) to the nearest center, and the distance between observed value and center is computed through Euclidean distance:$$ d= \arg \min \left\Vert \varphi \left({x}_i\right)\right.-{\left.{m}_k\right\Vert}^2 $$
- (c)Update the centers and recalculate the sum of the within-cluster variation E;$$ {m}_k=\frac{1}{N_k}{\displaystyle \sum_{x_j\in classk}\varphi \left({x}_j\right)},E={{\displaystyle \sum_{k=1}^K{\displaystyle \sum_{j=1}^{N_k}\left\Vert \varphi \left({x}_j\right)-{m}_k\right\Vert}}}^2 $$
- (d)
Repeat step (b) and (c), until the E is stable or less than a threshold value.
B denotes variances of between-clusters, W are the variances of within-cluster, k is the clusters’ count, n indicates the count of monitoring station. The high CH value indicates the well clustering results. In this research, we chose the count of clusters when the CH index is highest. The kernel k-means function was used the “kern lab” package [13] in R [14] in our study.
Temporal period and trend decomposition
Empirical mode decomposition (EMD), an adaptive time series data analysis model [15,16], was used to identify the temporal trend of malaria vectors in terms of density in the clustered surveillance sites. In EMD analysis, the original data is decomposed into a series of modes without requiring prior knowledge. Comparing with Fourier and wavelet decomposition, EMD has its own advantage. The Fourier analyses could transform data into the combination of sine and cosine functions with different frequencies, while wavelet analysis needs wavelet splines. These types of decomposition involve many spurious components due to serious restriction of the harmonic nature of the basis function [10,17]. The EMD method decomposes the data into several oscillatory components, corresponding to some physical phenomenon underlying the data, and the residual of the decomposition is the trend of the data [10].
The model decomposed from EMD methods is called the intrinsic mode function (IMF) [15,16]. An IMF should satisfy two conditions [10,15]: First, the difference between the number of extrema and the number of zero-crossings in the whole data span should be less than or equal to one. Second, the mean value of the envelope defined by the local maxima and the envelope defined by the local minima is zero at any point. We calculated the IMFs and the residuals of the density of the malaria vectors using the sifting method by taking the following steps:
Initialize r (t) = x(t), i = 0, k = 1, where x(t) is the time series of density of the malaria vectors. We define a threshold δ , the count of extreme points N.
Step1: Calculate the local maximum value and the local minimum value of the time-series r(t);
Step2: Calculate the upper envelope e _{max}(t) through interpolating the local maximum value, and the lower envelope e _{min}(t) through interpolating the local minimum value;
Step3: Calculate the local average value using the formula: m(t) = (e _{max}(t) + e _{min}(t))/2;
Step4: Appoint i = i + 1, then the proto-mode function (PMF) is calculated through the following formula: p _{ i }(t) = r(t) − m(t), r(t) = p _{ i }(t);
If SD < δ, then imf _{ k }(t) = p _{ i }(t), go to the following step; otherwise, go back to step 1–4;
Step 6: Appoint r(t) = r(t) − imf _{ k }(t). If the count of extreme points of r(t) is greater than N, then k = k + 1, i = 0, and go back to step 1; otherwise, stop the sifting process.
Where K is the number of IMFs and r is the trend of the original data. The first IMF (IMF1) is the component with the highest frequency. All of the IMFs were statistically tested against the null hypothesis [18]. In this research, the EMD method was processed with the code written by Rilling et al. [19]
Results
Overview of the dataset
An. lesteri was found at 13 surveillance sites in 7 provinces (Henan, Hubei, Guangdong, Hainan, Sichuan, Guizhou and Yunnan), with a peak of abundance each year from June to July. The density per site averaged 0.45 and ranged from 0.03 to 2.29 bites per man per night, the three highest density was observed at Anlu (2.29 bites per man per night) in Hubei Province, followed by Tongbai (0.80 bites per man per night) in Henan Province and Baisha (0.78 bites per man per night) in Hainan Province (Figure 3B).
An. minimus was captured at 13 national malaria surveillance sites in 5 provinces (Fujian, Guangdong, Hainan, Guizhou and Yunnan),with a peak of abundance each year from July to September. The density of An. minimus per site averaged 1.26 and ranged from 0.20 to 3.57 bites per man per night, the two highest density was observed at Menglian (3.57 bites per man per night) and Jinghong (3.50 bites per man per night) from Yunnan Province (Figure 3C).
An. dirus was captured at 8 national malaria surveillance sites in 3 provinces (Guangdong, Guizhou and Hainan). The density of An. dirus per site averaged 1.21 and ranged between 0.11 and 4.22 bites per man per night, the highest density was observed at Zengcheng (4.22 bites per man per night) in Guangdong Province, followed by Congjiang (1.76 bites per man per night) in Guizhou Province (Figure 3D).
Clustering distribution
The density distribution in three clusters of An. sinensis during 2005-2010
Clusters | Mean(SD) | Minimum | Percentile | Maximum | ||
---|---|---|---|---|---|---|
25% | 50% | 75% | ||||
Low density | 0.21 | 0.05 | 0.29 | 0.48 | 0.58 | 0.87 |
Medium density | 0.95 | 0.16 | 0.83 | 1.67 | 2.26 | 4.24 |
High density | 4.00 | 0.29 | 2.38 | 5.03 | 7.56 | 17.20 |
The temporal trend of identified clusters
Discussion
Like many former studies [20-25] indicating four primary malaria vector species including An.sinensis, An.lesteri, An.dirus and An.minimus existed in China, this study also showed the same malaria vectors existing with significant discrepancy of spatial distribution in the surveillance regions.
Although An.sinensis is an inefficient vector mainly because of its zoophilic habit, it is still considered an important vector of P. vivax malaria due to its wide distribution as well as high density. Malaria transmission had often taken place due to its high density under suitable conditions such as the open plains in the greater part of China between 25 ~ 33 north latitude [25-27].
An. lesteri was mainly distributed in low hill, hillock and shallow hilly plain between 22 ~ 33 north latitude which had still been found in 11 provinces with 205 counties through the survey conducted from 1998 to 2001 [3]. However, An. lesteri had been captured only in 13 counties from 7 provinces (Henan, Hubei, Guangdong, Hainan, Sichuan, Guizhou and Yunnan) in the surveillance sites during 2005–2010, indicating that the spatial distribution of this vector became smaller than before. The extensive agricultural use of pesticides and indoor insecticide spraying may have substantially reduced its distribution [23,28,29].
An.minimus had long been regarded as a principal malaria vector in China. In the past, mosquitoes of the An. minimus group were recorded from 16 provinces in southern China including Sichuan, Chongqing, Hubei, Henan, Anhui, Zhejiang, Yunnan, Guizhou, Guangxi, Hainan, Hunan, Jiangxi, Guangdong, Fujian, Hong Kong and Taiwan throughout southern China from Yunnan Province eastward and from Hainan Island northward to approximately 32.5 north latitude [30-33]. However, An. minimus was captured only in 5 provinces (Fujian, Guangdong, Hainan, Guizhou and Yunnan) from the surveillance sites during 2005–2010. It is still a major malaria vector in the hilly bordering regions of Yunnan Province near Myanmar [33].
In China, An.leucosphyrus was found for the first time by the writer in 1941 in Yunnan [34]. It was then rectified to An.dirus, a new species described by Peyton and Harrison [35], based on the detailed morphological comparison of specimens from Thailand. However, different from previous research that An.dirus had been captured from four provinces of Hainan, Yunnan, Guangxi and Tibet [36], it was only captured at 8 national malaria surveillance sites during 2005–2010 from 3 provinces (Guangdong, Guizhou and Hainan). An.dirus has so far been known as a malaria vector only on Hainan Island.
Further, this is the first study to systematically explore the temporal-spatial dynamics of An.sinensis in China using kernel k-means based on CH value [10]. The significant difference of population density may have resulted from available resources that more rice fields in the southern part of China than the others contributing to the suitable breeding sites of An.sinensis [37]. Different from the first and third cluster that a significant decrease of An. sinensis density was observed, the second cluster displayed a significant increase during 2005–2010 in the central part of China, which the increasing vectorial capacity leading to malaria re-emergence in Huang-Huai River of central China may be related to this change [38,39].
Conclusion
As malaria transmission relies mainly on Anopheles mosquitoes, the present results will serve as a baseline for the ongoing NMEP in China, which will help to address the knowledge gaps with respect to the malaria profile under different settings.
Notes
Declarations
Acknowledgments
This study was supported by the Fundamental Research Funds for the Central
Universities (no. BLX2014-08) and the National S & T Major Program (Grant No. 2012ZX10004-220 and 2012CB955503). We would like to thank all the staff enrolled in this national surveillance program for their excellent cooperation.
Authors’ Affiliations
References
- Ministry of Health. Malaria surveillance project in China. 2005.Google Scholar
- Lu B. The present status of mosquito-borne diseases in China. Program and Abstracts of 4th Intern Meeting of Tropical Medicine and parasitology. 1998;18-20.Google Scholar
- Guo C. Advance in the research of Anopheles anthropophagus in China. China Tropical Med. 2005;5(2):396–9.Google Scholar
- Qu F. Research of mosquitoes during recent fifty years in China. Chinese J Parasitological Parasitic Dis. 1999;17(5):264–6.Google Scholar
- Hay SI, Guerra CA, Tatem AJ, Noor AM, Snow RW. The global distribution and population at risk of malaria:past, present, and future. Lancet Infect Dis. 2004;4(6):327–36.View ArticlePubMed CentralPubMedGoogle Scholar
- Greenwood BM, Fidock DA, Kyle DE, Kappe SH, Alonso PL, Collins FH, et al. Malaria: progress, perils, and prospects for eradication. J Clin Invest. 2008;118(4):1266–76.View ArticlePubMed CentralPubMedGoogle Scholar
- Zhou G, Minakawa N, Githeko A, Yan G. Spatial distribution patterns of malaria vectors and sample size determination in spatially heterogeneous environments: a case study in the west Kenyan highland. J Med Entomol. 2004;41(6):1001–9.View ArticlePubMedGoogle Scholar
- Wang DQ, Tang LH, Liu HH, Gu ZC, Zheng X. Application of Structural Equation Models for Elucidating the Ecological Drivers of Anopheles sinensis in the Three Gorges Reservoir. PLoS One. 2013;8(7):e68766. doi:10.1371/journal.pone. 0068766.View ArticleGoogle Scholar
- Wang DQ, Tang LH, Gu ZC, Zheng X, Jiang WK. Malaria transmission potential in the three gorges reservoir of the Yangtze River, China. Biomed Environ Sci. 2013;26(1):54–62.PubMedGoogle Scholar
- Hu MG, Jia L, Wang JF, Pan YP. Spatial and temporal characteristics of particulate matter in Beijing, China using the empirical model decomposition method. Sci Total Environ. 2013;458–460:70–80.View ArticlePubMedGoogle Scholar
- Dhillon I, Guan Y, Kulis B. A unified view of kernel k-means, spectral clustering and graph partitioning. UTCS Technical Report. 2004.Google Scholar
- Caliński T, Harabasz J. A dendrite method for cluster analysis. Commun Stat. 1974;3:1–27.View ArticleGoogle Scholar
- Karatzoglou A, Smola A, Hornik K, Zeileis A. kernlab - an S4 package for kernel methods in R. J Stat Softw. 2004;11:1–20.Google Scholar
- R Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2012.Google Scholar
- Huang NE, Shen Z, Long SR, Wu MC, Shih HH, Zheng Q, et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc R SocLond A. 1998;454:903–95.View ArticleGoogle Scholar
- Huang NE, Wu MC, Long SR, Shen SSP, Qu W, Gloersen P, et al. A confidence limit for the empirical mode decomposition and Hilbert spectral analysis. Proc R SocLond A. 2003;459:2317–45.View ArticleGoogle Scholar
- Klionski DM, Oreshko NI, Geppener VV, Vasiljev AV. Applications of empirical mode de-composition for processing nonstationary signals. Pattern Recognition Image Anal. 2004;18:390–9.View ArticleGoogle Scholar
- Wu Z, Huang NE. A study of the characteristics of white noise using the empirical modedecomposition method. Proc R SocLond A. 2004;460:1597–611.View ArticleGoogle Scholar
- Rilling G, Flandrin P, Goncalves P. On empirical mode decomposition and its algorithms. Proceedings of IEEE-EURASIP workshop on nonlinear signal and image processing; 2003.Google Scholar
- Chow CY. Malaria vector in China. Chin J Entomol. 1991;Special Publ(6):67–79.Google Scholar
- Ministry of Health Disease Prevention and Control Bureau. Handbook for malaria control and prevention. Beijing: People’s Hygiene Publishing House Press; 2007.Google Scholar
- Sinka EM, Bangs JM, Manguin S, Chareonviriyaphap T, Patil PA, Temperley HW, et al. The dominant Anopheles vectors of human malaria in the Asia-Pacific region: occurrence data, distribution maps and bionomic précis. Parasi Vectors. 2011;4:89.View ArticleGoogle Scholar
- Tang LH, Qian HL, Xu SH. Malaria and its control in the People’s Republic of China. Southeast Asian J Trop Med Public Health. 1991;22:467–76.PubMedGoogle Scholar
- Rueda LM, Ma Y, Song GH, Gao Q. Notes on the distribution of Anopheles (Anopheles) sinensis Wiedemann (Diptera: Culicidae) in China and the status of some Anopheleshyrcanus group type specimens from China. Proc Entomol Soc Wash. 2005;107:235–8.Google Scholar
- Sleigh AC, Liu XL, Jackson S, Li P, Shang LY. Resurgence of vivax malaria in Henan Province. China Bull World Health Organ. 1998;76:265–70.Google Scholar
- Xu JJ, Feng LC. Studies on the Anopheles hyrcanus group of mosquitoes in China [in Chinese]. Acta Entomol Sinica. 1975;18:77–98. 76 pls.Google Scholar
- Ma SF. Studies on the Anopheles (A.) sinensis group of mosquitoes in China, including four new sibling species [in Chinese]. Sinozoologia. 1981;1:59–70. 54 pls.Google Scholar
- Ho C, Chou TC, Ch’En TH, Hsueh AT. The Anopheles hyrcanus group and its relation to malaria in east China. Chin Med J. 1962;81:71–8.PubMedGoogle Scholar
- Liu C. Comparative studies on the role of Anopheles anthropophagus and Anopheles sinensis in malaria transmission in China [in Chinese]. Zhonghua Liu Xing Bing Xue Za Zhi. 1990;11:360–3.PubMedGoogle Scholar
- Ho C, Feng LC. Studies on malaria in new China. Chin Med J. 1958;77:533–51.PubMedGoogle Scholar
- Dong XS. The malaria vectors and their ecology in Yunnan Province. Chinese J Parasitic Dis Control. 2000;13(2):144–7.Google Scholar
- Wang XZ, Du ZW, Lu YY, Zhu GJ, Huang R, Gu Y, et al. Bionomics of An. minimus and its role in malaria transmission in south parts of Yunnan. Chinese J Vector Biology Control. 1999;10(5):343–6.Google Scholar
- Lu BL. Fauna Sinica, Insecta, Vol. 9, Diptera. CulicidaeI I. Beijing: Science Press; 1997.Google Scholar
- Shi W, Zhou X, Zhang Y, Zhou X, Hu L, Wang X, et al. Investigation on malaria vectors in western part of China-Myanmar border. Chin J Parasitol Parasit Dis. 2011;9(2):134–7.Google Scholar
- Sweet WC, Feng LC, Chow CY, Hsu SC. Anophelines of Southwestern Yunnan and their relation to malaria. J Natl Mal Soc. 1942;1:25–32.Google Scholar
- Colless DH. Further notes on the systematic of the Anopheles leucosphyrus group. Proc Roy Ent Soc London. 1957;26:131–9.Google Scholar
- Harrison. PA Medical entomology studies XII. Contrib Am Ent Inst. 1980;17:195.Google Scholar
- Xiao X, Boles S, Liu J, Zhuang D, Frolking S, Lia C, et al. Mapping paddy rice agriculture in southern China using multi-temporal MODIS images. Remote Sens Environ. 2005;95(4):480–92.View ArticleGoogle Scholar
- Pan JY, Zhou SS, Zheng X, Huang F, Wang DQ, Shen YZ, et al. Vector capacity of An. sinensisin malaria outbreak areas of central China. Parasit Vectors. 2012;5:136.View ArticlePubMed CentralPubMedGoogle Scholar
Copyright
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.