1. Introduction
Hainan Island, located in the South China Sea and separated from Leizhou Peninsula by the Qiongzhou Strait, makes up the majority of the Hainan Province, China. It was once one of the worst malaria endemic areas in China. In 2003, Hainan Province started implementing the Global Fund Malaria Project, which greatly contributed to the implementation of Hainan’s malaria control program. Since 2010, Hainan had carried out malaria elimination according to the Malaria Elimination Action Plan (2010-2020), the number of malaria cases continued to decline, and realized “No local malaria cases reported” in 2012 [1]. On December 31, 2019, it officially confirmed that Hainan Province had reached the Malaria elimination target.
Anopheles mosquito is the sole vector of malaria. There are 37 species of Anopheles mosquito on Hainan Island. Exophagic, anthropophilic An. dirus and endophilic, anthropophilic An. minimus were the principal vectors in the period of high malaria prevalence in Hainan Island. An. sinensis, An. jeyporiensis candidiensis and several other Anopheles were secondary vectors [2]. With the development of society and economy, large areas of forests have been cut down to grow tropical cash crops such as rubber, areca nut and fruit, destroying breeding habitats for wild inhabiting mosquito species, and with long-term targeted mosquito control measures, the distribution and abundance of An. dirus and An. minimus have been drastically decreased [3]. An. sinensis is a type of exophagic and endophilic mosquito that feeds on both human and livestock blood. It is not an efficient malaria vector compared with An. dirus and An. minimus, but it is widely distributed on Hainan Island and has a large species population. With the development of agricultural mechanization, the drastic reduction of farm cattle has relatively increased the chances of the An. Sinensis mosquito, which used to feed on cattle blood, feeding on human blood, An. sinensis may become the mosquito species responsible for a low prevalence of malaria in Hainan Island, so it is necessary to monitor the spatial and temporal distribution of An. sinensis.
The dry season and the wet season are clearly different in Hainan, as the dry season (from November to April) is warm and dry while the wet season (from May to October) is hot and humid. The peak season of malaria overlapped with wet season during the malaria epidemic [4]. This study took Wanning County as the sample area, and analyzed the spatial distribution characteristics of adult An. sinensis mosquitoes and their relationship with environmental factors in wet season, in order to formulate tailored mosquito control measures to maintain elimination of malaria in Hainan Island.
2. Materials and Methods
2.1. Study Area
Wanning County is located on the southeast coast of Hainan Island with a tropical monsoon climate featuring ample sunshine and precipitation. Wanning is mountainous in the west, hilly in the central and plain in the southeast. The three terrains account for 46%, 32% and 22% respectively, and the average elevation is 123 meters. The topography and geomorphology of Wanning is a good representation of that of Hainan Island.
Wanning used to be a malaria-endemic area. After joining the Malaria Programme of China Global Fund Round One in 2003,the incidence of malaria in Wanning began to decrease significantly. The elimination of malaria has been carried out comprehensively in Wanning since 2010 and reached the standard of elimination in 2017 [5]. The historical malaria hyper-endemic area of Wanning lies in the western and northwestern mountainous regions where An. minimus and An. dirus were highly effective vectors. Although mosquito surveillance has shown that the population density of An. sinensis is significantly higher than that of An. dirus and An. minimus at villages with high malaria incidence, An. sinensis is not considered a major malaria vector and has not drawn the same attention as An. minimus and An. dirus during the high malaria transmission era [6].
2.2. Mosquito Surveillance
In this study, 36 adult mosquito monitoring sites were set up in 36 villages, 26 of which were located in the three historical malaria endemic towns (Nanqiao, Sangengluo and Beida) and 10 in non-malaria-endemic areas (Figure 1).
Figure 1. Study area and monitoring sites.
At each monitoring site, adult mosquitoes were collected outdoors twice a month using human landing catch technique from May to October 2009, once in the middle of the month and once at the end of the month. Due to the time of peak activity of An. sinensis is in the first half of the night [7] [8], mosquito sampling was carried out from 7 p.m. to 12 p.m. Informed consent was obtained from all the participants. Wanning Center for Disease Control and Prevention (CDC) provided medical supervision of participants. All collected mosquitoes were sent to the laboratory for classification. Densities of An. sinensis adult mosquito were measured as the number of An. sinensis adult mosquitoes captured per person-hour.
2.3. Environmental Data Acquisition
GPS receiver was used to record the geographical information (latitude, longitude and altitude) of each sampling site and the distance to the nearest household.
The distances from each sampling site to the nearest river and the nearest highway were calculated using spatial analysis in ArcGIS9.3 software (Environmental Systems Research Institute, Inc.).
The paper map of soil type distribution of Wanning County was obtained from Hainan Administration of Surveying, Mapping and Geoinformation. Soil type of each sampling site was extracted using spatial analysis in ArcGIS9.3 after the paper map being digitized.
The land use types of the sampling sites were classified into three categories: residential area, woodland and grassland respectively.
The temperature, air pressure and relative humidity of the sampling sites were recorded at the beginning of sampling using unified barothermohygrographs.
2.4. Data Analysis
Global Spatial autocorrelation analysis: Global spatial autocorrelation describes the overall distribution of observed objects in space to determine if there was spatial aggregation. Moran’s I index was used to measures the spatial aggregation patterns and aggregation level of attribute variables across the study area. The Moran’s I index was tested using standardized statistic Z values. Values of Moran’s I index ranges from −1 to 1, a positive value indicates tendency toward clustering while a negative value indicates tendency toward dispersion, a value of zero typically indicates no spatial autocorrelation, and the greater the absolute value of Moran’s I is, the higher the degree of autocorrelation in spatial distribution is [9] [10]. In this study, the spatial autocorrelation tool in Arcgis10.2 was used to measure the value of Moran’s I.
Space-time scanning cluster analysis: Scanning statistics are used to discover and evaluate the aggregation of events in pure time, pure space, or space-time dimensions by continuously moving scan window in time and/or space dimensions and recording the number of observations and expectations in the window for each position [11]. SaTScan9.4 software was used to carry out the purely spatial and spatio-temporal scanning cluster analysis of the density distribution of adult An. sinensis mosquitoes.
Because density was continuous data, the normal model was used to conduct the pure space clustering in order to find the high density aggregation area of An. sinensis in each month. For spatio-temporal analysis, the mosquito densities of each sampling site in each month were divided into 1, 2, 3, 4, and 5 levels from low to high, and ordinary model was used to conduct spatio-temporal clustering to clarify the spatial and temporal aggregation of An. sinensis in the whole surveillance period. The unit of space was defined by the coordinates of the sampling sites and the unit of time was one month. The maximum spatial cluster size was set to include up to 50% of the total sampling sites and the maximum size of the time frame was 50% of the surveillance period. For Monte-Carlo inference, 999 replications were performed.
Trend surface analysis: Trend surface analysis is a method based on low-order polynomials of spatial coordinates for estimating a regular grid of points from scattered observations. This approach is aimed to model the overall distribution of properties throughout space. It will sketch the global trend of distribution as a simplified surface. A trend surface model is a particular case of a bivariate regression model with two independent variables, the coordinates X and Y and a dependent variable, the thematic variable Z to be modeled [12] [13]. In this study, ArcGIS 10.2 software was used to analyze trend surface and draw trend analysis graph.
Logistic regression analysis: Bivariate logistic regression models were fitted to identify environmental variables associated with the presence of adult An. sinensis mosquitoes. Then, multivariate logistic regression model, incorporating variables significant at P < 0.05 in the univariate analysis, was developed using a forward-stepwise procedure and assessed using the likelihood ratio test. All regression analyses were carried out with SPSS 20.0.
3. Results
3.1. General
Of the 36 sampling sites, four sites (Jiarun, Toucun, Shanniutian and Minfeng) were stopped sampling from July and one site (Qiaozhong) was stopped sampling in October due to lack of sampling personnel, all the five sites are located in the historical malaria high endemic region. Accumulated 1498 person-hours of sampling were carried out and a total of 619 An. sinensis mosquitoes were caught. The average density of An. sinensis captured in each month from May to October was 0.61, 0.37, 0.39, 0.41, and 0.29 per man-hour, respectively. The density in May was significantly higher than that in other months, and there was no significant difference in the density of other months (Friedman M test, M = 25.273, P = 0.000). Of the 31 sampling sites that have completed all six months of surveillance, the highest average density from May to October was 1.08 per man-hour observed in Shanjia (see Figure 2). The highest density of single
Figure 2. Average mosquito densities of 31 sampling sites during May and October, 2009 (unit: number of Anopheles sinensis captured per person-hour).
surveillance was 2.75 per man-hour observed in Shanjia in mid-May. No An. sinensis was captured during the whole surveillance period at two sampling sites, Yuhou and Liulian, which both lie in the historically low endemic area of malaria.
3.2. Spatial Autocorrelation Analysis
The global autocorrelation analysis showed that the distribution of An. sinensis in each month from May to October 2009 was spatial aggregation rather than random one. The highest level of aggregation of mosquito density was in October, the lowest level of aggregation was in June (Table 1).
3.3. Spatial and Temporal Scanning Clustering
Purely spatial scanning cluster analysis: It was shown that mosquito population was not distributed randomly in Wanning County according to analysis of purely spatial clustering of mosquito density from May to October 2009. With the maximum spatial cluster size of 50% of the total surveillance sites, the spatial clustering identified the most likely cluster with high An. sinensis mosquito density for each of the six months. The secondary cluster was only detected in September and was not statistically significant (Table 2). The cluster covered three sampling sites in May, two sites in June, sixteen sites in July, sixteen sites same as those of July in August, eight sites in September, and eight sites same as those of September in October (Figure 3). These results showed that the relatively high density clusters of each month were not fixed.
Spatio-temporal scanning cluster analysis: Spatio-temporal scanning cluster of An. Sinensis density in Wanning County from May to October, 2009 was tested by ordinary model. The most likely cluster with high mosquito density located in central area of Wanning County (Coordinates: 18.791N, 110.262E; radius: 5.56 km; LLR: 17.636; P = 0.001) where covered three sampling sites adjacent to Wanning Reservoir in non-malaria-endemic areas (Figure 4), and the time frame was from May to July. No secondary cluster was detected.
Table 1. Global spatial autocorrelation of An. sinensis mosquito density in Wanning County from May to October, 2009.
Table 2. Purely spatial An. sinensis clusters in Wanning County from May to October, 2009.
*LLR: Log likelihood ratio, **secondary cluster.
Figure 3. The monthly spatial clusters of An. sinensis in Wanning County from May to October, 2009.
Figure 4. The location of the spatio-temporal cluster of An. sinensis mosquitoes in Wanning from May to October, 2009 (the map was taken from the Google earth. Imagery date: 12/31/2009).
3.4. Trend Surface Analysis
Trend surface analysis was employed to identify the geographic trends of the mosquito density. A coordinate system was created (one axis for each direction with X for West-East and Y for South-North) and projections of densities (Z axis) reflected the variation trend of West-East and South-North orientation.
The result showed that density distribution of An. sinensis each month during May and October presented an obviously inversed “U” shape from west to east, and presented a downward trend from May to August from south to north, a less obviously “U” shape from September to October (see Figure 5).
3.5. Logistic Regression Modeling for Adult An. sinensis Mosquitoes Presence/Absence
Univariate analysis: Univariate two-class logistic regression analysis was used to analyze the correlation between various environmental factors and whether or not An. sinensis was captured.
Univariate analysis was used to select environmental factors that were significantly associated with the presence of adult An. sinensis mosquitoes. Bivariate logistic regression models combining presence of adult An. sinensis mosquitoes and various environmental characteristics of sampling site are presented in Table 3. The results showed that the presence of adult An. sinensis mosquito was more likely in residential areas and sand shale yellow red soil areas, and was negatively correlated with distance to highway, distance to river, temperature, air pressure, and relative humidity, and was positively correlated with the distance to household. The presence of adult An. sinensis mosquito had no correlation with the sampling month or the elevation of sampling site.
Figure 5. Trend surface map of density of An. sinensis in Wanning County from May to October, 2009.
Table 3. Bivariate logistic regression results for presence of adult An. sinensis mosquitoes and the environmental factors in Wanning County.
Multivariate analysis: In order to comprehensively evaluate the effect of environmental factors on the presence of adult An. sinensis population, multivariate logistic regression was used to develop minimum adequate models with all significant variables (P < 0.05) in Table 3. Stepwise logistic regression analysis showed that the five environmental factors (land use type, soil type, and distance to highway, air pressure and relative humidity) were the influencing factors of presence of adult An. sinensis (Table 4). Land use type and soil type had strong correlation with the presence of An. sinensis, and distance to highway, air pressure and relative humidity were negative correlation with the presence of An. sinensis. Three covariates, distance to river, distance to household and temperature, were excluded from the model. The R2 of the regression equation was 0.489, and using the predicted probability of 0.5 as a cut-off point for classification, the correct rate of classification was 80.3%.
4. Discussion
In China, An. sinensis is the main or sole malaria vector in the region north of 34˚ north latitude. It is an important malaria vector in the plains, especially the main vector for Plasmodium vivax in rice growing areas [14]. Hainan island lies south of 21˚ north latitude, the main vectors of malaria during the high-malaria transmission era were An. dirus. and An. minimus. They have a higher susceptibility to Plasmodium vivax and Plasmodium falciparum than An. Sinensis [15] [16]. The spatial and temporal distribution characteristics of An. sinensis on Hainan Island have been rarely studied since the beginning of this century. However, due to long-term targeted mosquito control and environmental changes caused by socioeconomic development, the mosquito populations of An. dirus and An. minimus have been greatly reduced in recent years [17].
Table 4. Stepwise logistic regression results for presence of adult An. sinensis mosquitoes and the environmental variables in Wanning County.
During the surveillance period of this study, only three An. dirus mosquitoes and no An. minimus mosquito were captured. An. sinensis has a significant population advantage over An. dirus and An. minimus in the population structure of malaria vector Anopheles mosquitoes on Hainan Island, the population distribution characteristics and malaria transmission efficiency of the An. sinensis need further attention in order to eliminate malaria on Hainan island.
Anopheles sinensis has a broad adaptability to the breeding environment. The main breeding place is the paddy field and its irrigation systems. The seasonal variation of the population size is closely related to climate change and farming system in different areas. In China, adult An. sinensis are active throughout the year in regions south of 27˚ north latitude, survive the winter in semi-hibernation in areas between 28˚ and 33˚ north latitude, and enter hibernation in winter in regions north of 33˚ north latitude [18]. This survey shows that An. sinensis are active throughout the wet season from May to October in Wanning County, and the density in May is higher than that in the other months. This is consistent with a survey conducted in several districts in Hainan island in 2012, in which the peak of adult mosquito density appeared from April to June [19], but is inconsistent with the observations obtained in China mainland which was from July to August [20]. Perhaps, the reason is that the temperature and humidity conditions in Hainan Island were suitable for breeding and developing of mosquitoes all year round, however, with the increase of precipitation and tropical low-pressure cyclone from July to October each year, mosquito breeding sites and habitats would be destroyed, resulting in a certain drop in mosquito density.
The analysis result showed that the density of adult An. sinensis had significant positive spatial correlation. In other words, mosquito densities of neighboring monitoring sites were closer, indicating the distribution of adult An. sinensis was dependent on the spatial location. Changes of the intensities of correlation in six months suggest that seasonal changes of some environmental factors affecting the density of adult An. sinensis may have an impact on the intensity of this dependence. Positive spatial correlation will inevitably result in spatial aggregation. The monthly space scanning statistics showed that there was a relative high density spatial aggregation area in each month, the changes of location and extent of monthly aggregation areas suggested that the intensity of mosquito density changes of different sampling sites were not consistent in different seasons. The space-time scanning from May to October showed that during the whole surveillance period, the high density aggregation area of adult An. sinensis was located in the northwest of Wanning reservoir and the time frame was from May to July. This area is not the traditional malaria high-endemic area in Wanning and is obviously inconsistent with the spatial distribution of malaria cases in Wanning County in 2009 [21]. Whether this can suggest that An. sinensis is not the main vector of malaria in Wanning County needs more in-depth research, because in addition to mosquito vectors, there are many other factors that affect the malaria epidemic such as the presence of malaria case and whether and how prevention and control measures are taken.
Besides the characteristics of spatial distribution of adult An. sinensis, this study explored the key natural environmental factors related with this distribution feature. Logistic regression analyses showed that land use type, soil type, distance to highway, air pressure, and relative humidity were major factors related with the presence of adult An. sinensis.
Adult An. sinensis was more likely to be found in residential areas than in woodland and grassland, as they are easier to find vertebrate host to suck blood in residential areas. Soil is closely related to insects, it not only indirectly affects insects as the plants grow, but also is places where insects live. Studies shows that the duration of immature life stages of the mosquitoes is inversely related to soil-substrate organic matter content and soil moisture has an impact on the density of mosquitoes in the area [22] [23]. Shililu reported that the absolute moisture content of the soil is not the only factor that determines the successful rate of hatching anopheline eggs, rather, soil moisture, together with the rate of drying, physiological factors associated with the age of the egg, and the type of soil in which the egg rests, likely influences survival [24]. The mechanism how soil type affects distribution of Anopheles mosquitoes in Wanning needs further investigation. This study shows that the closer the monitoring point is to the highway, the greater the possibility of presence of adult An. sinensis, this may be due to some unknown indirect correlations that require further analysis.
Meteorological factors directly affect the mosquito density, because they affect not only reproduction and survival of mosquitoes but also their activities. According to Ma Aimin’s research, total annual precipitation, mean annual air pressure, precipitation of the wettest quarter, and minimum temperature of the coldest month strongly influenced the distribution of An. sinensis [25]. The influence of temperature on mosquito populations is highly complex, as the effects of temperature on various processes are non-linear and difficult to predict. Some surveys have shown a positive correlation between mosquito density and temperature [26] [27] [28]. However, some other surveys showed a negative correlation as we observed in Wanning [29]. The relationship between air pressure and mosquito density also has different conclusions, these different conclusions may be due to differences in mosquito species and basic climatic conditions. For example, a survey conducted in the central Arctic showed that as air pressure rose, the probability of mosquito presence increased and the relationship between air pressure and activity level was also positive [30], while this study shows that the air pressure is inversely related to the presence of An. sinensis and consistent with the findings of several surveys conducted in China which have shown that there is a significant negative correlation between air pressure and mosquito density, and the variation trend of mosquito density can be predicted by using the air pressure fitting model [28] [31] [32]. Hainan Island has high relative humidity, the relative humidity recorded in this study ranged from 68% to 97% with an average of 80%. This survey shows that the relative humidity is inversely related to the presence of An. sinensis, which is consistent with Ferreira’s survey result that humidity above 75% had a negative effect on mosquito abundance [33].
Based on the close relationship between malaria and Anopheles mosquitoes and the high prevalence seasons of malaria in Hainan Island, this study carried out adult Anopheline mosquitos’ surveillance in Wanning County from May to October, 2009, and analyzed the spatial distribution characteristics and environmental factors of An. sinensis population. The results showed that the spatial distribution of the high density areas of adult An. sinensis was inconsistent with the spatial distribution of the high incidence areas of malaria in Wanning. The high prevalence area of malaria has been located in the mountainous forest area of southwest Wanning for many years [21], however, the high density area of An. sinensis in this survey was located in the central hilly area. The density of An. sinensis in May was significantly higher than that in the other months, but the peak of malaria epidemic in Wanning was usually in summer and autumn. Whether this indicates that An. sinensis is not the main vector for maintaining the low level of malaria prevalence in Wanning County needs further study because many socio-economic factors also affect the malaria epidemic in addition to mosquito vectors.
5. Conclusion
The spatial distribution of An. sinensis adult mosquitoes in Wanning shows strong spatial heterogeneity, which is related to geographical and meteorological factors. Using the method of spatial epidemiological analysis, the spatial distribution characteristics of An. sinensis mosquitoes and their influencing factors could be intensively explored to provide a basis for the local anti-mosquito strategies. Hainan Island is adjacent to Southeast Asia where malaria is widespread. As China builds Hainan Island into an international tourism island, population movement between Hainan Island and neighboring countries will be more frequent. This increases the risk of resurgence of malaria in Hainan Island, and monitoring and controlling mosquito vectors remains an important task of maintaining elimination of malaria in Hainan Island. To clarify the ecological habits of An. sinensis in Hainan Island for accurately forecasting the spatial and temporal distribution of An. sinensis, further comprehensive surveys are needed.
Acknowledgements
Special thanks to Chen X, Lin F, Fu JX, Hu MG, Tao Zh, Zhong ShCh, Zhong XH, Xu ChH, Lin DX, Zong ZhH, Liu WX, Zhou ChAn, Lin L, Huang DY, Ji XiY, Li R, Huang G, Tang YP, Li PSh, Shi CB, Huang DY, Lu JH, Huang GP, Chen HSh, Shi LB, Lin XL, Huang HW, Huang HR, Xu ZhW, Wan CQ, Huang ShJ, Huang ZhL, Huang ZhX, Wang DS, Qiu XH, Chen QH, Wang Y, Jiang HY, Wen XZh, Zhang YP, Lan TT, Cui DY, Lin SQ, Chen B, Lu LL, and Wu Sh for their toilsome field mosquitoes sampling.
This study was supported by the National Natural Science Foundation of China (No. 30872151).
NOTES
*These authors contribute equal to this work.