- Research
- Open Access
- Published:

# Modelling hotspots of the two dominant Rift Valley fever vectors (*Aedes vexans* and *Culex poicilipes*) in Barkédji, Sénégal

*Parasites & Vectors*
**volume 9**, Article number: 111 (2016)

## Abstract

### Background

Climatic and environmental variables were used successfully by using models to predict Rift Valley fever (RVF) virus outbreaks in East Africa. However, these models are not replicable in the West African context due to a likely difference of the dynamic of the virus emergence. For these reasons specific models mainly oriented to the risk mapping have been developed. Hence, the areas of high vector pressure or virus activity are commonly predicted. However, the factors impacting their occurrence are poorly investigated and still unknown. In this study, we examine the impact of climate and environmental factors on the likelihood of occurrence of the two main vectors of RVF in West Africa *(Aedes vexans* and *Culex poicilipes)* hotspots.

### Methods

We used generalized linear mixed models taking into account spatial autocorrelation, in order to overcome the default threshold for areas with high mosquito abundance identified by these models. Getis’ *G*i*(*d*) index was used to define local adult mosquito abundance clusters (hotspot).

### Results

For *Culex poicilipes*, a decrease of the minimum temperature promotes the occurrence of hotspots, whereas, for *Aedes vexans*, the likelihood of hotspot occurrence is negatively correlated with relative humidity, maximum and minimum temperatures. However, for the two vectors, proximity to ponds would increase the risk of being in an hotspot area.

### Conclusions

These results may be useful in the improvement of RVF monitoring and vector control management in the Barkedji area.

## Background

Rift Valley fever (RVF) is an emerging arboviral disease considered as a threat to human and animal health in Africa and Arabian Peninsula. RVF outbreaks have led to significant economic impact [1, 2]. This disease affects mainly cattle, sheep, and goats and is responsible for abortion among pregnant females and increasing new-born mortality [3, 4]. In humans, RVF infections are usually asymptomatic, but severe cases can result in complications including haemorrhagic syndromes which are associated with high mortality rates [5]. Humans are exposed to the Rift Valley fever virus (RVFV) by direct contact with infected animals, organs of dead animals, fluids, aerosols, raw milk, and bites of infected mosquitoes [6, 7]. RVFV is transmitted by several species of mosquitoes, mainly of the genera *Aedes* and *Culex*. In Senegal, after the RVF outbreak in Mauritania in 1987, an entomological and animal surveillance program was implemented in several biogeographic zones [8–10]. These studies highlighted several RVFV circulation. Thus, the RVFV was isolated from several mosquito species collected at the end of the rainy season in the Barkedji area [8, 11]. In this area, six mosquito species were found associated to RVFV but the two dominant vectors are *Culex poicilipes* and *Aedes vexans* [12–14]. All these vectors breed in temporary ponds that are flooded after the first rains usually in July. These temporary ponds represent the main source of water for people, livestock and wildlife in this area for at least 6 months per year [14, 15].

The impact of environmental and climatic factors on the heterogeneous spatial distribution of several species of mosquitoes is complex and poorly understood despite the valuable models developed in the recent years [16–19]. This is specially the case for vector-borne diseases such as RVF where knowledge of the relation between host, vectors, and environmental factors is useful for understanding the persistence, emergence, transmission and possible amplification of the virus.

The active circulation of the RVFV underlines the need for efficient surveillance and control tools. To assist prudent decision making in animal and human health resource allocation for surveillance and control, spatio-temporal models for outbreaks prediction have been developed using climate time series data such as rainfall and sea surface temperature and the Normalized Difference of Vegetation Index (NDVI) anomalies map to set up an early warning system for East Africa and Arabian Peninsula [20–22]. However, this system has not been effective in West Africa because the factors involved in the emergence and re-emergence of this virus appear to be different [17]. These included the lack of sensitivity of the system like threshold of NDVI and mosquito vectors involved in the transmission [23]. Further, in contrast in East Africa, above normal rainfall is not always associated to RVFV emergence in West Africa [24]. Thus, for West Africa, several studies have been conducted particularly in Senegal to generate an early warning system [18, 25–28]. One of the main approaches was to focus on the abundance of the main vectors involved in RVFV transmission. Precisely at Barkedji, a Bayesian model was implemented using Markov chain Monte Carlo (MCMC) methods [29]. The main purpose of this previous modelling approach was to define areas of high and low vector abundance, in order to first identify appropriate places for pastoralists to settle and then to minimise contact between vectors and hosts. The other objective was to provide information to the authorities for the implementation and adaptation of targeted vector control strategies based upon small-scale insecticide use after identification of areas to be treated. This statistical model was used to assess the effects of climatic and ecological determinants on the spatiotemporal dynamics of *Ae. vexans* and *Cx. poicilipes* at a local-scale. However, due to the lack of abundance threshold, the previous model developed does not allow to rank among high density areas. This ranking is necessary to reduce costs and surfaces to be treated but also manage the problem of vector resistance to insecticides and the issue of environmental pollution. Indeed, previous studies using satellite images in the Barkédji area have identified during the 2003 rainy season a total of 468 ponds in July [30] and 1354 ponds in August [31]. This large number of ponds demonstrates the need of classification tools in the presence of several sites with high mosquito densities to determine those of primary importance. This is especially important for low income countries where resources are scarce and appropriate guidance is essential for preventing mosquito bites or reducing vector abundance. For that purpose, we developed, in this study, a model to determine the impact of climate and environmental data on the occurrence of hotspots. This model would be an alternative to the lack of threshold. Thus, sites with high densities, predicted by the previous model [29] and with a high probability of being hotspots will be considered in priority as risk areas for any intervention program. We initially determined all hotspots of the two main vectors in our study area using Getis-Ord statistic.

The Getis–Ord Gi* hotspot cluster statistic is one of the possible approaches that can be used for local spatial analysis [32]. This statistics identifies those clusters with higher values (hotspot) in magnitude than expected to be found by random chance. The Gi* statistic measures the degree of spatial clustering of a local sample and indicates how different it is from the expected value which is the mean of the whole dataset. It was used to determine the distribution of multiple mosquito species and to show how they may currently and in the future be affected by climate change, control interventions or other factors [33–35]. It can also be used to help decision maker in identifying areas at risk to target for strategies such as spraying and larval control [36–39].

Statistical species distribution models, including generalized linear models, are commonly used as a tool in decision-making [40–42]. These models also commonly use environmental and climate factors as explanatory variables but their spatial variation influencing predicted distributions are rarely quantified. Hence, the knowledge of climatic and environmental factors affecting the occurrence of these hotspots is required as information-support for decision-making.

The aims of this study were to identify local adult *Ae. vexans* and *Cx. poicilipes* abundance clusters (hotspots) using the Getis’ *G*i*(*d*) index and to model the occurrence of these hotspots by generalized linear mixed models (GLMM). We adopted this modelling approach to quantify the effect of cumulative rainfall, relative humidity, maximum temperature, minimum temperature and the NDVI and the Euclidian distance to larval breeding sites (temporary ponds) on the distribution of vector hotspots. This approach could assist in identifying areas at higher risk and further provide important information supporting decision making including how different intervention actions should be spatially allocated.

## Methods

### Study area

The study was performed around the Barkedji village (14°52’02”W, 15°16’41”N), in the Ferlo area (central north of Senegal) during the 2005 to 2006 rainy seasons (Fig. 1). This area is characterised by a hot dry climate, short rainy season (from June to October), and long dry season (November to May), with annual rainfall ranging from 300 to 500 mm and a number of rainy days around 35.8 [17]. The highest temperatures are usually realised during October with average daily maximum of 45 °C and the lowest temperatures during December with 16 °C. Many temporary ponds flood at the beginning of the rainy season and represent the main source of water for people and livestock. These ponds are the natural habitats of many species of birds, reptiles and rodents and the oviposition and resting sites for mosquito vectors of RVFV. Cattle farming is the main activity of the people in this area with other livestock mainly sheep and goats. The houses are scattered around temporary ponds and are usually composed of small villages with huts.

### Mosquito population sampling and explanatory variables

Mosquitoes were sampled from 79 sites once every 2 weeks from July to December 2005 and 2006 by CO_{2}-CDC light traps (model 512; John W. Hock Company, Gainesville, FL) for a period of 12 h (6 pm to 6 am). Upon collection, mosquitoes were euthanised and identified using available morphological keys [43, 44]. The geographical coordinates of each site were recorded with a Garmin GPSMAP 76S (Olathe, KS, USA) hand-held global positioning unit and projected in UTM Zone 28 N allowing the Euclidian distances from each site to the nearest pond to be estimated using a spatial map of the region.

Rainfall, relative humidity, maximum and minimum temperature, NDVI and distance from the nearest pond were the variables selected to model vector hotspots occurrence. The meteorological variables were collected from automatic weather stations (BWS 200 Campbell Scientific) installed in the Barkedji area. Rainfall represented the cumulative rainfall 15–20 days prior to trapping. We calculated the averages for maximum and minimum temperature, and relative humidity of the sampling period. Remotely-sensed Moderate resolution Imaging Spectroradiometer (MODIS) data set were sourced from the National Aeronautic and Spatial Administration (http://neo.sci.gsfc.nasa.gov/view.php?datasetId=MOD13A2_E_NDVI.) and the NDVI with spatial and temporal resolution of 250 m and 16 days intervals analysed. These explanatory variables chosen are considered likely to be important determinants of the two vectors distribution.

Clustered localities with high or low *Ae. vexans* and *Cx. poicilipes* abundances were identified using the Getis-ord Gi* statistic [32] in R software with “spdep” package [45]. Statistically significant (at a level of 0.05) clusters of sites with high vector abundances were identified as hotspots with Z scores > 1.96, while Z scores of < -1.96 represented coldspots with low vector abundance. The spatial relationship among sites was conceptualised using the inverse distance, which is most appropriate for continuous point datasets because closer sites have larger influences on the computation for each target site than sites that are further away. Data on presence/absence of hotspot was coded as binary factor (presence =1 or absence = 0). At each site, the presence or absence of hotspot (Z scores > 1.96 for presence and Z scores < 1.96 for absence) were then determined.

### Statistical modelling

The generalized linear mixed model (GLMM) with a binomial distribution was fitted taking into account spatial autocorrelation with the Getis-Ord Gi* Z score as the dependent factor. We modelled the likelihood of the presence of hotspot by introducing a random intercept effect between sites [46]. The general form of the logistic regression models is:

Where *p*
_{
st
} is the probability of vector hotspot presence at site s at fortnight t, *β* is a vector of coefficient; *X*
_{st} is a vector of covariates for site s at time t; and *b*
_{s} is the random effect for site s. The parameter *b*
_{s} is distributed normally with mean zero and variance *σ*
^{2}(*b*
_{
s
} ~ *Normal*(0, *σ*
^{2})). GLMM is an alternative for taking into account the spatial autocorrelation within sites [47]. The R package “glmmML” [48] and previously used R source codes were modified and utilised for model fitting [47]. For the GLMM model, we used the methodology described [47].

To verify possible spatial autocorrelation in the parsimonious model residuals, Moran’s I correlograms plot (with a lag interval of 5000 m) from the Pearson residuals of the generalized linear model (GLM) and the GLMM was used [49, 50].

### Model selection

The goodness of fit was assessed for all models tested using the Akaike information criterion (AIC) criteria [51]. AIC criteria is define as:

where L is the maximum log-likelihood of the model and K is the number of parameters in the model. The best approximating model is the model with the lowest AIC.

We tested different combinations of the explanatory variables in the model. Each model were fitted and ranked by their AIC values. We also used Akaike weights *w*
_{
i
} to calculate the relative probability of each model compared to the best model. The Akaike weights for model *i* is defined as [52]:

Where ∆_{
i
} is the difference between the AIC for model *i* and the best model (with the lowest AIC) and the sum in the denominator is over all candidate models (*j* = 1,…,*R*), the denominator is the sum of the relative likelihoods for all candidate models. Finally, for each covariate, the relative importance was evaluated by an indicator obtained by summing the Akaike weights for all models containing the covariate [52]. A 95 % confidence set of models was obtained by ranking model by their Akaike weights and added them successively until the sum exceeds 0.95. We used the area under receiver operating curve (ROC) for the best model to assess its predictive capacity [53, 54]. Values of area under a ROC curve (AUC) greater than 0.7 indicate that the predictive capacity of model may be considered acceptable, excellent for AUC ≥ 0.9, and poor if AUC value is lower than 0.6 [55]. ROC curve were done using ROCR package [56] with the R software [57].

We used graphical methods, quantile-quantile plots with simulation approach and partial residuals plots, to assess the goodness of model fit [58]. Then we applied these methods to the most parsimonious model. The partial residual graph is a plot of the values of a specific covariate in the model versus its partial residuals. Partial residuals *r*
_{
par
} are defined as:

Where y is the observed data (1 or 0), \( \widehat{p} \) is probability estimated by the fitted model, *X* is the model covariate, and \( {\widehat{\beta}}_X \) is the estimated coefficient for the covariate *X* [58]. A linear partial residual plot indicates that linear assumption in the model is adequate. However, a non-linear of partial residual plot suggests that model linear assumption may be not suitable [47]. The partial residual plots were only done for covariates in the most parsimonious model with smoothed curves fitted.

To improve predictive performance of the model, we applied a model averaging method to take into account model uncertainty [52]. Model averaged predictions are more robust than those derived from a single model [47]. Model coefficients were estimated and predictions made taking into account parameter uncertainty [52].

## Results

The localities of high and low clusters of vectors abundances were identified with the Z scores computed by Getis-Ord Gi* (Figs. 2 and 3). For *Aedes vexans*, there were 35 hotspots sites in 2005 and 32 hotspots sites in 2006. For *Culex poicilipes*, there were 19 hotspots sites in 2005 and 21 hotspots sites in 2006 (Table 1). For *Aedes vexans*, in 2005, the first half of August had the highest hotspots numbers (11) whereas in 2006 it was during the second half of September (11). In 2005, *Culex poicilipes* had 15 hotspots in the second half of October and 12 hotspots in 2006 (Table 1).

The most parsimonious model for *Cx. poicilipes* (AIC = 416.6) contained the distance to the nearest pond and the minimum temperature (Tables 2 and 3). The standard deviation of random-effects estimated was 3.07 ± 0.93 indicating relatively large variation amont collection sites.

There was negative correlation between distance to nearest waterpond, minimum temperature and vector hotspot occurrence. However, this model was just 1.26 times better than the model incoporating cumulative rainfall (evidence ratio = 0.155/0.123) and (0.155/0.116) and 1.34 times better than the model with NDVI, distance to the nearest pond and the minimum temperature parameters. The best model for *Ae. vexans* (AIC = 571.3) contained distance to the nearest pond, minimum and maximum temperature and relative humidity (Tables 2 and 3). The standard deviation of random-effects estimated was 0.37 ± 0.27 indicating weak variation among collection sites. These variables were all negatively correlated with this vector hotspot occurrence. Nevertheless, this model was 1.3 (0.368/0.283) times better than the model containing NDVI variable and 2.2 (0.368/0.168) times better than the model with cumulative rainfall.

Significant positive spatial autocorrelation was observed in the model residuals (Fig. 4), but at short lag distances of less than 1.0 km for *Cx. poicilipes*. This indicates that spatial autocorrelation is a problem for sites located close to each other. Observations from these sites can not be considered as independent. However, the spline correlogram plot of the Pearson residuals generated by the GLMM model (Fig. 4) showed absence of spatial correlation. This indicates that the mixed model successfully takes into account the spatial autocorrelation existing within sites. Concerning *Ae. vexans*, the spline correlogram plot demonstrated absence of spatial autocorrelation in Pearson residuals of GLM and GLMM models (Fig. 5).

All graphical methods used to test the goodness of model fit for both species indicated no considerable deviation from the models hypothesis. The first one was the quantile-quantile plots which points were inside the simulated 95 % confidence interval (Additional file 1: Figure S1-A). Partial residuals plots, the second one, with smoothed curve, indicated that the linear hypothesis for each covariate appeared to be suitable (Additional file 2: Figure S1-B and Additional file 3: Figure S1-C).

The cross-validation area under the curve (AUC) was 0.96 and 0.75 for *Cx. poicilipes* and *Ae. vexans,* respectively (Fig. 6), indicating stronger predictive ability and that the framework of the most parsimonious models for each species were appropriate.

For *Ae. vexans*, the most important variable that influenced abundance was maximum temperature (relative importance index = 0.994), followed by minimum temperature (0.992), relative humidity (0.975) and distance to the nearest pond (0.958) (Table 2).

For *Cx. poicilipes*, minimum temperature (0.991) and distance to the nearest pond (0.898) were the most mportant variables that influenced the occurrence of mosquito hotspots (Table 3).

## Discussion

Our study focused on *Cx. poicilipes* and *Ae. vexans* because they are the most abundant and are considered as the main vectors of RVF in the area. They also exhibit a high degree of interaction with the main vertebrate hosts of the virus [11, 14]. Other mosquito species like *Mansonia uniformis, Mansonia africana, Aedes fowleri, Aedes ochraceus*, were found associated with RVF virus. However, they are uncommon and accidently infected. *Mansonia uniformis* and *Ma. africana* are occasionally abundant but their association with aquatic vegetation limits their spatial distribution [59, 60]. Here we developed a statistical model of the RVF vectors hotspot distribution where the probability of hotspot presence was conceptualised as a function of climatic and environmental variables. Our approach demonstrated that climatic and environmental variables were very important for determining the distribution of RVF vectors hotspot in Barkedji area.

There was inverse relationship between distance to pond and predicted probability of hotspot areas for both species. The high abundance of the main RVF vectors near temporary ponds agreed with previous studies that estimated the average flying range less than 600 m and 650 m for *Cx. poicilipes* and *Ae. vexans,* respectively [24]. A study showed also that their abundances decrease linearly up to 843 and 1,394 m from a given pond for *Cx. poicilipes* and *Ae. vexans,* respectively [12]. In addition, *Cx. poicilipes* were generally found near some ponds located in the Ferlo riverbed [11, 14, 61]. The distance of a host to the ponds is an important factor that could affect the risk of being infected by RVFV. It may be an indicator to find the potential vertebrate reservoir host of the virus.

During this study, the lowest and highest temperature ranges recorded were 20.57 °C and 45.39 °C. The negative correlation between minimum temperature and occurrence of vector hotspots indicate that the vectors are very sensitive to lower temperatures that impact negatively on their survival and ability to transmit pathogens such as dengue virus [62, 63]. In addition, vector biology studies showed that higher transmission of RVFV was observed at higher temperatures for *Aedes* and for *Culex* [64–66]. However, the recent emergence of RVFV in Madagascar illustrates that RVFV is also able to circulate under more temperate condition [67, 68].

Previous studies in West Africa have shown that, from 1961 to 2003, periods of RVFV emergence did not coincide with years of high rainfall [16, 18] and that *Ae. vexans* and *Cx. poicilipes* abundance and total rainfall were not correlated [17]. In West Africa, intra seasonal variability of rainfall was suspected to have more impact on mosquitoes dynamics [69]. Therefore, the dynamics of filling ponds was recently identified as a favourable factor in the abundance of the two vectors [18]. RVF outbreaks in Tanzania during 2006–2007 was preceded with periods of above normal rainfall and the occurrence of outbreaks was associated with total amount of rainfall above a threshold (405.4 mm) [70].

The relative performance of several alternative models tested could arise because of correlations between environmental variables, which would consequently limit the amount of additional information derived from varying specific combinations between the variables tested. However, in the model with averaged coefficients, only the two factors (NDVI and cumulative rainfall) had a high likelihood of being included in the parsimonious model. These variables, although not statistically significant in the averaged model, had a great influence in predicting vector hotspot occurrence. Model-averaged could be used to predict the presence of hotspot. The concept of inference is to reduce model selection biais effects on linear regression coefficient estimates [52].

## Conclusions

Vector control is challenging in Barkedji area due to the heterogeneous spatial distribution of RVF vectors across landscape as well as their breeding habitat. The identification of vectors hotspots is a key step in the implementation of most efficient surveillance and control strategies. Thus, this monotoring should target hotspot sites. Pastoralists should be advised to avoid hotspots areas in order to reduce the probability of contact between host and vectors before the government could be able to make more safer and modern water sources for the herds in the future.

## References

Guillaud M, Le Guenno B, Wilson ML, Desoutter D, Gonzalez JP, Digoutte JP. [Prevalence of antibodies against Rift Valley fever virus in sheep and goats in Senegal]. Ann Inst Pasteur Virol. 1988;139(4):455–459.

Wilson ML, Chapman LE, Hall DB, Dykstra EA, Ba K, Zeller HG, et al. Rift Valley fever in rural northern Senegal: human risk factors and potential vectors. Am J Trop Med Hyg. 1994;50(6):663–75.

Botros B, Omar A, Elian K, Mohamed G, Soliman A, Salib A, et al. Adverse response of non-indigenous cattle of European breeds to live attenuated Smithburn Rift Valley fever vaccine. J Med Virol. 2006;78(6):787–91.

Hunter P, Erasmus BJ, Vorster JH. Teratogenicity of a mutagenised Rift Valley fever virus (MVP 12) in sheep. Onderstepoort J Vet Res. 2002;69(1):95–8.

Kamal SA. Pathological studies on postvaccinal reactions of Rift Valley fever in goats. Virol J. 2009;6:94.

Brown JL, Dominik JW, Morrissey RL. Respiratory infectivity of a recently isolated Egyptian strain of Rift Valley fever virus. Infect Immun. 1981;33(3):848–53.

Al-Hazmi M, Ayoola EA, Abdurahman M, Banzal S, Ashraf J, El-Bushra A, et al. Epidemic Rift Valley fever in Saudi Arabia: a clinical study of severe illness in humans. Clin Infect Dis. 2003;36(3):245–52.

Fontenille D, Traore-Lamizana M, Diallo M, Thonnon J, Digoutte JP, Zeller HG. New vectors of Rift Valley fever in West Africa. Emerg Infect Dis. 1998;4(2):289–93.

Zeller HG, Fontenille D, Traore-Lamizana M, Thiongane Y, Digoutte JP. Enzootic activity of Rift Valley fever virus in Senegal. Am J Trop Med Hyg. 1997;56(3):265–72.

Thonnon J, Picquet M, Thiongane Y, Lo M, Sylla R, Vercruysse J. Rift valley fever surveillance in the lower Senegal river basin: update 10 years after the epidemic. Tropical Med Int Health. 1999;4(8):580–5.

Diallo M, Lochouarn L, Ba K, Sall AA, Mondo M, Girault L, et al. First isolation of the Rift Valley fever virus from

*Culex poicilipes*(Diptera: Culicidae) in nature. Am J Trop Med Hyg. 2000;62(6):702–4.Diallo D, Talla C, Ba Y, Dia I, Sall AA, Diallo M. Temporal distribution and spatial pattern of abundance of the Rift Valley fever and West Nile fever vectors in Barkedji, Senegal. J Vector Ecol. 2011;36(2):426–36.

Diallo M, Nabeth P, Ba K, Sall A, Ba Y, Mondo M, et al. Mosquito vectors of the 1998–1999 outbreak of Rift Valley Fever and other arboviruses (Bagaza, Sanar, Wesselsbron and West Nile) in Mauritania and Senegal. Med Vet Entomol. 2005;19(2):119–26.

Ba Y, Sall AA, Diallo D, Mondo M, Girault L, Dia I, et al. Re-emergence of Rift Valley fever virus in Barkedji (Senegal, West Africa) in 2002–2003: identification of new vectors and epidemiological implications. J Am Mosq Control Assoc. 2012;28(3):170–8.

Diallo D, Ba Y, Dia I, Lassana K, Diallo M. Use of insecticide-treated cattle to control Rift Valley fever and West Nile virus vectors in Senegal. Bull Soc Pathol Exot. 2008;101(5):410–7.

Ndione J-A, Diop M, Lacaux JP, Gaye AT. Variabilité intra-saisonnière de la pluviométrie et émergence de la fièvre de la vallée du Rift dans la vallée du fleuve Sénégal : nouvelles considérations. Climatologie. 2008;5:83–97.

Ndione J-A, Besancenot JP, Lacaux JP, Sabatier P. Environnement et épidémiologie de la fièvre de la vallée du Rift (FVR) dans le bassin inférieur du fleuve Sénégal. Environ Risques Santé. 2003;2(3):176–82.

Soti V, Tran A, Degenne P, Chevalier V, Lo Seen D, Thiongane Y, et al. Combining hydrology and mosquito population models to identify the drivers of Rift Valley fever emergence in semi-arid regions of West Africa. PLoS Negl Trop Dis. 2012;6(8):e1795.

Mondet B, Diaïté A, Ndione J-A, Fall AG, Chevalier V, Lancelot R, et al. Rainfall patterns and population dynamics of

*Aedes*(*Aedimorphus*)*vexans arabiensis*Patton, 1905 (Diptera: Culicidae), a potential vector of Rift Valley Fever virus in Senegal. J Vector Ecol J Soc Vector Ecol. 2005;30(1):102–6.Anyamba A, Linthicum KJ, Small J, Britch SC, Pak E, de La Rocque S, et al. Prediction, assessment of the Rift Valley fever activity in East and Southern Africa 2006–2008 and possible vector control strategies. AmJTrop Med Hyg. 2010;83(2 Suppl):43–51.

Anyamba A, Linthicum KJ, Small JL, Collins KM, Tucker CJ, Pak EW, et al. Climate teleconnections and recent patterns of human and animal disease outbreaks. PLoS Negl Trop Dis. 2012;6(1):e1465.

Anyamba A, Chretien J-P, Small J, Tucker CJ, Formenty PB, Richardson JH, et al. Prediction of a Rift Valley fever outbreak. Proc Natl Acad Sci U S A. 2009;106(3):955–9.

Anyamba A, Linthicum KJ, Mahoney R, Tucker CJ, Kelley PW. Mapping potential risk of Rift Valley fever outbreaks in African savannas using vegetation index time series data. Photogramm Eng Remote Sens. 2002;68(2):137–45.

Ba Y, Diallo D, Kebe CMF, Dia I, Diallo M. Aspects of bioecology of two Rift Valley fever virus vectors in Senegal (West Africa):

*Aedes vexans*and*Culex poicilipes*(Diptera: Culicidae). J Med Entomol. 2005;42(5):739–50.Bicout DJ, Sabatier P. Mapping Rift Valley Fever vectors and prevalence using rainfall variations. Vector Borne Zoonotic Dis. 2004;4(1):33–42.

Porphyre T, Bicout D, Sabatier P. Modelling the abundance of mosquito vectors versus flooding dynamics. Ecol Model. 2005;183(2):173–81.

Ndiaye PI, Bicout DJ, Mondet B, Sabatier P. Rainfall triggered dynamics of

*Aedes*mosquito aggressiveness. J Theor Biol. 2006;243(2):222–9.Clements AC, Pfeiffer DU, Martin V, Pittliglio C, Best N, Thiongane Y. Spatial risk assessment of Rift Valley fever in Senegal. Vector Borne Zoonotic Dis. 2007;7(2):203–16.

Talla C, Diallo D, Dia I, Ba Y, Ndione JA, Sall AA, et al. Statistical modeling of the abundance of Vectors of West African Rift Valley Fever in Barkedji. Senegal PLoS One. 2014;9(12):e114047.

Vignolles C, Lacaux JP, Tourre YM, Bigeard G, Ndione JA, Lafaye M. Rift Valley fever in a zone potentially occupied by

*Aedes vexans*in Senegal: dynamics and risk mapping. Geospat Health. 2009;3(2):211–20.Lacaux J, Tourre Y, Vignolles C, Ndione J, Lafaye M. Classification of ponds from high-spatial resolution remote sensing: Application to Rift Valley Fever epidemics in Senegal. Remote Sens Environ. 2007;106(1):66–74.

Ord JK, Getis A. Local spatial autocorrelation statistics: distributional issues and an application. Geogr Anal. 1995;27(4):286–306.

Bisanzio D, Giacobini M, Bertolotti L, Mosca A, Balbo L, Kitron U, et al. Spatio-temporal patterns of distribution of West Nile virus vectors in eastern Piedmont Region, Italy. Parasites Vectors. 2011;4(1):1–11.

Kamdem C, Fouet C, Etouna J, Etoa F-X, Simard F, Besansky NJ, et al. Spatially explicit analyses of anopheline mosquitoes indoor resting density: implications for malaria control. PLoS One. 2012;7(2):e31843.

Barrera R. Spatial stability of adult Aedes aegypti populations. AmJTrop Med Hyg. 2011;85(6):1087–92.

Pinault LL, Hunter FF. New highland distribution records of multiple

*Anopheles*species in the Ecuadorian Andes. Malar J. 2011;10(1):236.Khormi HM, Kumar L, Elzahrany RA. Modeling spatio-temporal risk changes in the incidence of dengue fever in Saudi Arabia: a geographical information system case study. Geospatial Health. 2011;6(1):77–84.

Kelly-Hope LA, Hemingway J, McKenzie FE. Environmental factors associated with the malaria vectors

*Anopheles gambiae*and*Anopheles funestus*in Kenya. Malar J. 2009;8(268):1–8.Kampen H, Medlock JM, Vaux AG, Koenraadt CJ, van Vliet AJ, Bartumeus F, et al. Approaches to passive mosquito surveillance in the EU. Parasites Vectors. 2015;8(1):1–13.

Drake JM, Hassan AN, Beier JC. A statistical model of Rift Valley fever activity in Egypt. J Vector Ecol. 2013;38(2):251–9.

Diuk-Wasser MA, Brown HE, Andreadis TG, Fish D. Modeling the spatial distribution of mosquito vectors for West Nile virus in Connecticut, USA. Vector-Borne Zoonotic Dis. 2006;6(3):283–95.

Wang J, Ogden NH, Zhu H. The impact of weather conditions on

*Culex pipiens*and*Culex restuans*(Diptera: Culicidae) abundance: a case study in Peel region. J Med Entomol. 2011;48(2):468–75.Edwards F. Mosquitoes of the Ethiopian Region. III. Culicine adults and pupae. British. London: Museum (Nat Hist); 1941.

Diagne N, Fontenille D, Konate L, Faye O, Lamizana MT, Legros F, et al. Les anophèles du Sénégal : liste commentée et illustrée. Bull Soc Pathol Exot. 1994;87(4):267–77.

Bivand R, Bernat A, Carvalho M, Chun Y, Dormann C, Dray S, Halbersma R, Lewin-Koh N, Ma J, Millo G: The spdep package. Comprehensive R Arch Network, Version

*03–13*2005.Hosmer Jr DW, Lemeshow S. Applied logistic regression. John Wiley & Sons; 2004.

Rhodes JR, McAlpine CA, Zuur A, Smith G, Ieno E. GLMM applied on the spatial distribution of koalas in a fragmented landscape. In: Mixed effects models and extensions in ecology with R. New York: Springer-Verlag. 2009;469-492.

Broström G, Holmberg H. Generalized linear models with clustered data: Fixed and random effects models. Computational Statistics & Data Analysis. 2011;55(12):3123–3134.

Getis A, Ord JK. The analysis of spatial association by use of distance statistics. Geogr Anal. 1992;24(3):189–206.

BjØrnstad ON, Falck W. Nonparametric spatial covariance functions: estimation and testing. Environ Ecol Stat. 2001;8(1):53–70.

Akaike H: Information theory and an extension of the maximum likelihood principle. In: Second international symposium on information theory1973: Akademinai Kiado: 267–281.

Burnham KP, Anderson DR. Model selection and multimodel inference: a practical information-theoretic approach. Springer Science & Business Media; 2003.

Pearce J, Ferrier S. Evaluating the predictive performance of habitat models developed using logistic regression. Ecol Model. 2000;133(3):225–245.

Fawcett T. An introduction to ROC analysis. Pattern Recogn Lett. 2006;27(8):861–74.

Brooker S, Hay SI, Bundy DA. Tools from ecology: useful for evaluating infection risk models? Trends Parasitol. 2002;18(2):70–4.

Sing T, Sander O, Beerenwinkel N, Lengauer T. ROCR: visualizing classifier performance in R. Bioinformatics. 2005;21(20):3940–1.

R Development Core Team R: R: A Language and Environment for Statistical Computing. 2.11.1 edn: R Foundation for Statistical Computing; 2013: 409–409.

Landwehr JM, Pregibon D, Shoemaker AC. Graphical methods for assessing logistic regression models. J Am Stat Assoc. 1984;79(385):61–71.

Sang R, Kioko E, Lutomiah J, Warigia M, Ochieng C, O’Guinn M, et al. Rift Valley fever virus epidemic in Kenya, 2006/2007: the entomologic investigations. AmJTrop Med Hyg. 2010;83(2 Suppl):28–37.

Chandra G, Ghosh A, Biswas D, Chatterjee S. Host plant preference of

*Mansonia*mosquitoes. J Aquat Plant Manag. 2006;44:142–4.Fontenille D, Traore-Lamizana M, Zeller H, Mondo M, Diallo M, Digoutte JP. Short report: Rift Valley fever in western Africa: isolations from

*Aedes*mosquitoes during an interepizootic period. Am J Trop Med Hyg. 1995;52(5):403–4.Turell MJ, Gargan 2nd TP, Bailey CL.

*Culex pipiens*(Diptera: Culicidae) morbidity and mortality associated with Rift Valley fever virus infection. J Med Entomol. 1985;22(3):332–7.Brady OJ, Johansson MA, Guerra CA, Bhatt S, Golding N, Pigott DM, et al. Modelling adult

*Aedes**aegypti*and*Aedes albopictus*survival at different temperatures in laboratory and field settings. Parasit Vectors. 2013;6:351.Brubaker JF, Turell MJ. Effect of environmental temperature on the susceptibility of

*Culex pipiens*(Diptera: Culicidae) to Rift Valley fever virus. J Med Entomol. 1998;35(6):918–21.Turell MJ. Effect of environmental temperature on the vector competence of

*Aedes taeniorhynchus*for Rift Valley fever and Venezuelan equine encephalitis viruses. Am J Trop Med Hyg. 1993;49(6):672–6.Turell M. Effect of environmental temperature on the vector competence of

*Aedes fowleri*for Rift Valley fever virus. Res Virol. 1989;140:147–54.Carroll SA, Reynes J-M, Khristova ML, Andriamandimby SF, Rollin PE, Nichol ST. Genetic evidence for Rift Valley fever outbreaks in Madagascar resulting from virus introductions from the East African mainland rather than enzootic maintenance. J Virol. 2011;85(13):6162–7.

Jeanmaire EM, Rabenarivahiny R, Biarmann M, Rabibisoa L, Ravaomanana F, Randriamparany T, et al. Prevalence of Rift Valley fever infection in ruminants in Madagascar after the 2008 outbreak. Vector-Borne Zoonotic Dis. 2011;11(4):395–402.

Mondet B, Diaïté A, Fall AG, Chevalier V. Relations entre la pluviométrie et le risque de transmission virale par les moustiques: cas du virus de la Rift Valley fever (RVF) dans le Ferlo (Sénégal). Environ Risques Santé. 2005;4(2):125–9.

Sindato C, Karimuribo ED, Pfeiffer DU, Mboera LE, Kivaria F, Dautu G, et al. Spatial and temporal pattern of Rift Valley fever outbreaks in Tanzania; 1930 to 2007. PLoS One. 2014;9(2):e88897.

## Acknowledgments

This study was partially supported by the EU project QWeCI (Quantifying Weather and Climate Impacts on health in developing countries; funded by the European Commission’s Seventh Framework Research Programme under the grant agreement 243964).

## Author information

### Authors and Affiliations

### Corresponding author

## Additional information

### Competing interests

The authors declare that they have no competing interests.

### Authors’ contributions

Conceived and designed the experiments: CT MD AD. Performed the experiments: CT YB ID DD. Analysed the data: CT MD DD AD. Contributed materials/analysis tools: CT AM AAS JAN DD AD. Wrote the paper: CT MD. All authors read and approved the final version of the manuscript.

## Additional files

### Additional file 1: Figure S1.

Quantile-quantile plot with 95 % pointwise confidence bounds; (A) for *Ae. vexans* and (B) for *Cx. poicilipes*. (TIF 37 kb)

### Additional file 2: Figure S2.

Partial residuals plots for covariates in the parsimonious model for *Cx. poicilipes*; (A) for distance, (B) for maximum temperature, (C) for minimum temperature and (D) for relative humidity. (TIFF 51 kb)

### Additional file 3: Figure S3.

Partial residuals plots for covariates in the parsimonious model for *Ae. vexans*; (A) for distance and (B) for minimum temperature. (TIFF 25 kb)

## Rights and permissions

**Open Access** This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

## About this article

### Cite this article

Talla, C., Diallo, D., Dia, I. *et al.* Modelling hotspots of the two dominant Rift Valley fever vectors (*Aedes vexans* and *Culex poicilipes*) in Barkédji, Sénégal.
*Parasites Vectors* **9**, 111 (2016). https://doi.org/10.1186/s13071-016-1399-3

Received:

Accepted:

Published:

DOI: https://doi.org/10.1186/s13071-016-1399-3

### Keywords

- Rift Valley Fever
- Mixed effects model
- Logistic model
- Spatial distribution model
- Hotspot
- Mosquitoes
- Barkédji