Research | Open | Published:
Spatial distribution modelling of Culicoides (Diptera: Ceratopogonidae) biting midges, potential vectors of African horse sickness and bluetongue viruses in Senegal
Parasites & Vectorsvolume 11, Article number: 341 (2018)
In Senegal, the last epidemic of African horse sickness (AHS) occurred in 2007. The western part of the country (the Niayes area) concentrates modern farms with exotic horses of high value and was highly affected during the 2007 outbreak that has started in the area. Several studies were initiated in the Niayes area in order to better characterize Culicoides diversity, ecology and the impact of environmental and climatic data on dynamics of proven and suspected vectors. The aims of this study are to better understand the spatial distribution and diversity of Culicoides in Senegal and to map their abundance throughout the country.
Culicoides data were obtained through a nationwide trapping campaign organized in 2012. Two successive collection nights were carried out in 96 sites in 12 (of 14) regions of Senegal at the end of the rainy season (between September and October) using OVI (Onderstepoort Veterinary Institute) light traps. Three different modeling approaches were compared: the first consists in a spatial interpolation by ordinary kriging of Culicoides abundance data. The two others consist in analyzing the relation between Culicoides abundance and environmental and climatic data to model abundance and investigate the environmental suitability; and were carried out by implementing generalized linear models and random forest models.
A total of 1,373,929 specimens of the genus Culicoides belonging to at least 32 different species were collected in 96 sites during the survey. According to the RF (random forest) models which provided better estimates of abundances than Generalized Linear Models (GLM) models, environmental and climatic variables that influence species abundance were identified. Culicoides imicola, C. enderleini and C. miombo were mostly driven by average rainfall and minimum and maximum normalized difference vegetation index. Abundance of C. oxystoma was mostly determined by average rainfall and day temperature. Culicoides bolitinos had a particular trend; the environmental and climatic variables above had a lesser impact on its abundance. RF model prediction maps for the first four species showed high abundance in southern Senegal and in the groundnut basin area, whereas C. bolitinos was present in southern Senegal, but in much lower abundance.
Environmental and climatic variables of importance that influence the spatial distribution of species abundance were identified. It is now crucial to evaluate the vector competence of major species and then combine the vector densities with densities of horses to quantify the risk of transmission of AHS virus across the country.
African horse sickness (AHS) and Bluetongue (BT) are vector-borne viral diseases affecting equids and ruminants (mainly cattle and small ruminants), respectively. These two viruses are biologically transmitted by females of several species of Culicoides (Diptera: Ceratopogonidae) biting midges. Bluetongue is distributed worldwide, affecting almost all continents (except Antarctica) while AHS occurs in sub-Saharan Africa, with rare incursions into Europe and Asia [1,2,3]. The importance of these two arboviral diseases derives from their potential for rapid spread and their major economic impact due to direct mortalities, restriction of animal movements, surveillance and vaccination costs, as shown by the recent outbreaks of BT in Europe [4,5,6] and AHS in Africa [7, 8]. Their importance requires immediate notification to the World Animal Health Organization .
Following the important AHS outbreak in Senegal in 2007 [7, 8], several studies were initiated in the Niayes area (in the western part of the country) in order to better characterize Culicoides diversity and ecology [10,11,12,13,14], and the impact of environmental and climatic variables on the dynamics of major vector species (proven and suspected vectors) . These studies recorded for the first time the presence of Culicoides oxystoma and showed that this species was extremely abundant in the Niayes, both in black-light and in horse-baited traps. Because of its abundance, close contact with horses and its known or suspected role in the transmission of other arboviruses (Akabane, epizootic haemorrhagic disease and bluetongue viruses), these studies concluded that it should be considered as a potential vector of AHS virus. These studies also showed that Culicoides population peak of abundance occurred at end of the rainy season in September and October.
In order to anticipate future events, prevent and better control Culicoides-borne disease outbreaks, it is essential to identify areas with higher risk of Culicoides-borne pathogens transmission at a national or regional level. This allows characterizing vector species distribution in Senegal. This is essential to determine if vectors are present throughout the country, if their abundance is highly variable or not and thus identify high-risk transmission areas. Combined with horse population distribution, this enables to target vaccination in these high-risk areas in the event of an outbreak. The first step towards this aim implies improving our knowledge on the distribution of vectors in Senegal. Thus, a nation-wide entomological trapping campaign was carried out to model the abundance of proven and potential vectors of AHS virus (AHSV) and/or BT virus (BTV) in Senegal and map their distribution.
Several studies demonstrated the ability to predict the presence and/or abundance of Culicoides midges using meteorological and environmental variables mainly derived from satellite imagery. These models were established using statistical techniques such as discriminant analysis [16,17,18,19,20], generalized linear models such as logistic regression [21,22,23,24,25] and more recently, data-mining techniques such as random forests [17, 26, 27].
The aim of this study is to model the spatial distribution of Culicoides species which are known or suspected to be vectors of AHSV and/or which are the most abundant in Senegal. The models could help stakeholders to make decisions for the control of these diseases. For this, entomological data were obtained through a nationwide trapping campaign and models for each species were developed to map their abundance. Three different modeling approaches were compared in order to find a simple yet robust method for each species. The first method consists in a spatial interpolation by ordinary kriging of Culicoides abundance data and does not necessitate environmental data. The other two consist of analyzing the relationship between Culicoides abundance, environmental and climatic data to model abundance in unsampled areas and were carried out by implementing generalized linear models and random forest models.
In 2012, a nationwide Culicoides trapping campaign was organized to collect information on the spatial distribution of Culicoides species in Senegal. Overall, with the help of the veterinary services, 108 sites holding equids were initially selected as follows: 3 sites per department and 3 departments per region in 12 (of the 14) regions of Senegal. Both the Institut Sénégalais de Recherches Agricoles (ISRA) team and veterinary services’ officers carried out trapping over two consecutive collection nights in each site at the end of the rainy season (between September and October). This timeframe was chosen because it corresponds to the peak of abundance of Culicoides in the Niayes region [13,14,15, 28]. Geographical coordinates of sampling sites (Fig. 1) were recorded using Garmin© global positioning system receivers (accurate to within 10 meters). Culicoides species were collected using Onderstepoort black-light suction traps (Onderstepoort Veterinary Institute South Africa) positioned close to equids (for further information on trap set up see [13, 14]). Identification of Culicoides species was carried out as explained by .
Maximum abundance of the two nights was considered as the best estimate of the population present as abundance can decrease rapidly in sub-optimal trapping conditions . If the mean abundance between two nights is used instead of the maximum abundance, then the abundance recorded will be decreased in the sites. Mean is typically meaningful when data are normally distributed. Maximum is much more meaningful here as we are sure that the site holds at least this number of Culicoides corrected by an unknown trap efficiency factor. Entomological data are provided in Additional file 1: Table S1.
Climatic and environmental data
Climatic and environmental variables characterizing favorable habitats for Culicoides were selected based on literature review of presence and abundance models [15,16,17,18, 20, 26, 30,31,32,33,34,35,36,37]. A total of 21 variables were selected belonging to 5 categories (temperature, vegetation index, precipitation, land cover and livestock density). Climatic variables included day-time and night-time land surface temperature (DLST and NLST) and rainfall. Environmental variables included the normalized difference vegetation index (NDVI), percentages of 3 land-cover classes (water, savannah and forest), and ruminant host densities.
Temperature, NDVI and rainfall were extracted from satellite images from September 1, 2011 to October 31, 2012. The DLST and NLST images were obtained from MODIS sensor (http://modis.gsfc.nasa.gov) and NDVI images were derived from Spot-Vegetation (https://rs.vito.be/africa/en/data/Pages/vegetation.aspx). For DLST, NLST and NDVI, the maxima, minima and means were calculated for this 14 months period. Precipitation data were derived from National Oceanic and Atmospheric Administration (http://iridl.ldeo.columbia.edu/expert/SOURCES/.NOAA/.NCEP/.CPC/.FEWS/.Africa/.DAILY/.ARC2/.daily/.est_prcp/). The maximum, mean and cumulative rainfalls were calculated. Land cover data was derived from the Senegal Land Cover Database produced by the Food and Agriculture Organization of the United Nations (FAO) in the framework of Global Land Cover Network activities (http://www.fao.org/geonetwork/srv/en/main.home?uuid=545be438-bc87-480b-83ec-ba3f4e486daf). Percentage of surface covered by water, particularly favorable for Culicoides [15, 26] savannah and forest in the 1 km2 pixels comprising the trap were extracted. Ruminant (cattle, sheep and goats) livestock density data was obtained from the global distribution of livestock maps produced by the FAO (http://www.fao.org/ag/againfo/resources/en/glw/GLW_dens.html).
All data layers were clipped on the Senegalese territory and projected in the same projection system with a spatial resolution of 1 km2 using R2.10.1 statistical language environment  using of the R-package: maps , mapdata , maptools , raster , rgdal  and sp . Climatic data are provided in Additional file 2: Figure S1.
Three methods were used to map the abundances of targeted species of interest: ordinary kriging and two regression techniques: generalized linear models and random forests.
Spatial interpolation by ordinary kriging
Ordinary kriging is an interpolation technique that aims to predict the value of a variable in an unsampled site  taking into account the spatial dependence of the data. It enables converting point data (data measured in specific sample sites) into raster data (images). In spatial statistics, the empirical semivariance is described as follows:
where, z (x i ) are data at a particular location, h is a distance between ordered data, N(h) the number of paired data at a distance of h.
First experimental semi-variograms (a graph of semivariances versus distances between ordered data) are plotted to visualize statistical dependence values. Then a model is adjusted according to semivariogram's points to assess the statistical dependencies between sites and thus determine the maximal interpolation distance (over this distance interpolations can no longer be carried out).
Generalized linear model
One of the distributions commonly used to model count data is the Poisson distribution. However, this distribution assumes equidispersion of count data. Culicoides data generally show strong overdispersion . Overdispersion was tested according to . The statistical approach was planned as follows: in case of significant overdispersion of the residuals of the Poisson model, a negative binomial model was tested, and in case again of significant overdispersion of residuals of the negative binomial model, a hurdle model was tested [47, 48].
Random Forest model
Random Forest (RF) method is a robust ensemble learning technique for data analysis which consists of a set of classifications and regressions trees constructed from sub-samples of the complete data set . This approach can be applied to model either the presence probability by performing classifications for qualitative variables (Random Forest classification) [17, 26, 27] or abundances by performing regressions for quantitative variables (Random Regression Forest) . For more information on RF models and their application to model Culicoides presence and abundance, see . By bootstrapping the data and by randomly changing the predictive variable sets over the different tree, RF increase diversity among regression trees. Each of the k regression trees is grown using another bootstrap subset of the original data set and the nodes are split using the best split variable among a subset of m randomly selected predictive variables . RF parameters, which are the number of trees (k) and the number of predictive variables used to split the nodes (m), are user-defined. In this study, to allow for error convergence, k was set to 500 and m to 4.
The Root means squared error (RMSE) was calculated for internal validation of general linear models (GLM) and RF models. Statistical analysis and modeling were performed with the R2.10.1 statistical language environment  using of the R-package: rgdal , splancs , gstat , maptools  and randomForest .
Between September, 17th and October, 14th 2012, among the 108 sites initially selected, 98 sites were sampled and 10 sites were not prospected because of logistical issues. In two of the 98 sampled sites the trapping results were not exploitable and thus were not considered (one conservation issue and one battery failure).
A total of 1,373,929 specimens of the genus Culicoides belonging to 32 different species were collected in 96 sites during the survey (Table 1). At least 32 species were collected; there may be other species among the 1.8% of captured Culicoides that were not identified. Among those specimens, some belong to the Milnei and Similis groups but could not be identified more precisely as one sex was missing. Some individuals belonging to the Imicola group were shown to world experts of Culicoides but could not be identified; at this stage it is unclear whether they are variants of known species or new species. Two species were recorded for the first time in Senegal: Culicoides murtalai  and Culicoides ochrothorax , increasing the number of described Culicoides species found in Senegal from 53  to 55.
Each of the following five species represented more than 10% of the catch: Culicoides oxystoma, Culicoides enderleini, Culicoides imicola, Culicoides miombo and Culicoides trifasciellus. As found previously in the Niayes region, C. oxystoma (26.9% of catches) was the most abundant species and was present in 93 out of the 96 sites sampled (97% of sites). Percentages of individuals caught compared to the total (and frequencies of sites where the species was present) catches were 23.89% (95.83%) for C. enderleini; 14.38% (95.83%) for C. imicola; 12.87% (82.29%) for C. miombo and 10.76% (56.25%) for C. trifasciellus, respectively (Table 1).
We chose to model the distribution of five species: the first four most abundant (C. oxystoma, C. enderleini, C. imicola and C. miombo) and also Culicoides bolitinos as it is a proven vector of AHSV in South Africa . Culicoides bolitinos was collected in low numbers (0.32% of catches) and in 57.29% of the prospected sites.
For the five Culicoides species of veterinary interest, observed abundance data (maximum abundance of the two consecutive nights) are shown in Fig. 1. Overall, three species had similar (yet not identical) patterns: C. oxystoma, C. imicola and C. enderleini. They were abundant in the west (groundnut basin), in the south-east and in some sites in the north-east (river delta). The distribution of C. miombo differed mainly in that it was very rare in the north-eastern sites while the distribution of C. bolitinos was much sparser.
Regarding C. oxystoma, the sites where abundances were the highest are located in west central (Fatick, Kaolack and Diourbel regions), in the far west (Dakar region) and south of Senegal (Kolda region) (Fig. 2). For C. enderleini, the highest abundances were observed in the south (Kolda and Kedougou), in the groundnut basin (specifically in Fatick and Kaolack) and northwest of the Saint-Louis. For C. imicola, the highest abundances were recorded in the south and southeast (Kolda, Tambacounda and Kedougou) and in the center-west (Thies, Diourbel, Fatick and Kaolack). Culicoides miombo and C. bolitinos were very rare in northern Senegal. The highest abundances of C. miombo were obtained in Kolda, Tambacounda and in the groundnut basin (Fig. 2). Culicoides bolitinos abundances were very low. Culicoides bolitinos highest densities were noticed in Kedougou region and in the groundnut basin (particularly in Fatick, Kaffrine and Kaolack).
The variogram results determined the maximal interpolation distances for each of the five species. These distances were: 20 km for C. oxystoma, 48 km for C. imicola, 13 km for C. enderleini, 25 km for C. miombo and 10 km for C. bolitinos (Fig. 3). Different models were fitted to the variograms: spherical models for C. oxystoma, C. imicola and C. bolitinos, a circular model for C. miombo and an exponential model for C. enderleini. Abundance maps by ordinary kriging were developed respecting the maximal interpolation distance for each species (Fig. 4). For C. oxystoma, C. enderleini, C. imicola and C. miombo, abundances were predicted to be very high in west-central and southern Senegal and decrease gradually towards the north. However, high abundances of C. enderleini were predicted in the northwest of Senegal and, to a lesser extent, for C. oxystoma. For C. bolitinos, the highest abundances were modelled in the extreme south of Senegal and south of the groundnut basin (specifically in Kaffrine region).
To produce abundance maps for the entire Senegalese territory, GLM and RF models were developed. Table 2 presents the 10 most important predictors for each species according to GLM models. Maximum NDVI was among the three most important variables for all species. Average rainfall was also among the 3 most important variables for 3 of the 5 species (C. oxystoma, C. enderleini and C. imicola). Other environmental variables that influence Culicoides abundances were day and night temperatures and percentage cover of water bodies for C. oxystoma and C. imicola; night temperatures for C. enderleini; day and night temperatures for C. miombo; livestock density and percentage cover of savannah for C. bolitinos. Abundance maps based on GLM modeling showed that for all five species of veterinary interest, predicted abundances are very low along the Senegal River Valley and very high in the south (Fig. 5). High abundances of C. imicola and C. oxystoma are predicted on the entire country although for C. imicola areas of high and low abundance are strongly interlinked (areas of high abundance are often close to areas of low abundance). Culicoides imicola and C. enderleini predicted abundances are particularly strong in the southern third of the country and medium in the middle third. The distribution of C. miombo in the south is even more pronounced: abundant in the southern third, it is less abundant in the central and almost absent in the northern third of the country. Not surprisingly, C. bolitinos predicted abundances are lower than those of the other species. The area where it is most abundant is the delta in Casamance in southeast Senegal and in southwest Senegal (Fig. 5).
Table 3 lists the 10 most important predictor variables from RF model of the different Culicoides species. Rainfall and/or NDVI were the most important variables influencing the abundance of the five Culicoides species. This is particularly true for C. imicola, C. enderleini and C. miombo, since the 3 most important variables are average rainfall and minimum and maximum NDVI. Abundance of C. oxystoma was mostly determined by average rainfall and day temperature and less so by NDVI (only ranked 7th by order of importance). Abundance of C. enderleini was mostly determined by average precipitation, NDVI (minimum and maximum) and average day temperature; that of C. imicola by average precipitation and maximum NDVI. Abundance of C. miombo was mostly driven by maximum NDVI followed by average precipitation, minimum NDVI and average night temperature. Culicoides bolitinos had a particular trend; the above environmental and climatic variables had a lesser impact on its abundance.
Predicted Culicoides abundance maps provided by RF model are shown in Fig. 6. For the five species, an increasing abundance gradient from north to south is predicted. For C. oxystoma, C. enderleini and C. imicola, predicted abundances are very high in the southern third of the country, and for C. miombo abundances are high to very high in the southern and middle thirds. For all the species, predicted abundances are very low in the north of the country, particularly in the north of Ferlo and in the Senegal River valley. Culicoides oxystoma is predicted to be particularly abundant along the western coast from Dakar to Ziguinchor and in The Gambia along the river. This species also has very high predicted abundances in the south, southeast and center-west of Senegal (Fig. 6). Culicoides enderleini and C. imicola exhibit a similar pattern with very high abundances predicted in Casamance (Ziguinchor, Sedhiou and Kolda), eastern Senegal (Tambacounda and Kedougou) and the center-west of the groundnut basin (Fatick and surrounding). Regarding C. miombo and C. bolitinos, very low abundances were predicted in the Niayes’ area. Culicoides bolitinos predicted abundances are lower than those of other species; the area where its abundance is predicted to be the highest is the south and southeast, particularly in Kedougou region.
Prediction errors measured by RMSE (on log10(n+1) scale) between observed and predicted values are lower for the RF model than for the GLM model for all five species (Table 4). Among RF models, predictions errors are lowest for C. imicola (RMSE = 0.31).
The identification of suitable areas for five targeted Culicoides species in Senegal was performed with three different modelling techniques: a spatial interpolation method, kriging and two statistical modelling methods using meteorological and environmental variables, GLM and RF models. Interpolation by kriging yielded incomplete prediction maps. To obtain complete maps, more intensive entomological survey would be necessary, especially for C. bolitinos and C. enderleini. Given the logistical efforts required to conduct such survey and identify Culicoides species, this option seems unrealistic.
The comparison of predicted abundance maps from each model showed consistent spatial structures. Although the predictions of maximum abundances provided by kriging were limited by extrapolation distances, abundance maps by kriging and RF models identified similar structures whereas GLM models provided abundance maps slightly different from those resulting from kriging and RF models. GLM outcomes differed for C. oxystoma particularly in Tambacounda and Kedougou, for C. enderleini in Fatick, Kaffrine, Kaolack and western Saint-Louis, for C. imicola in Thies, Fatick, Kaolack, for C. miombo in Dakar, Kaffrine and Tambacounda and for C. bolitinos in Kaolack, Kaffrine and Fatick.
Prediction errors between observed and predicted abundances were lower for the RF model than for the GLM model. Recent studies have also shown that RF approaches provide good modelling results both for regression  and classification [17, 26, 27] trees. The effect of climatic and environmental variables will be discussed for each species based on the results from RF models since they provided better estimates of abundances.
Regarding C. oxystoma, this is, to our knowledge, the first spatial distribution model published. Studies on Culicoides dynamics in Senegal show that this species is very frequent and abundant in the Niayes area [13, 14]. Our study shows that this species is widespread and abundant over the entire country, with an increasing gradient from north to south, and high predicted abundance in the west, near the coast and around river deltas.
Average rainfall was the most important variable influencing C. oxystoma, C. enderleini and C. imicola abundance, and the second most important variable for C. miombo. This variable seems to limit the abundance of these species in the north, more arid (see isohyets representation of Senegal in Fig. 2). In addition to rainfall, day and night temperatures strongly influenced C. oxystoma abundance. Its highest abundance was observed in the groundnut basin and southern Senegal. Although ecological drivers of C. oxystoma remain poorly described, this species is known to be abundant in South Korea  and in Kagoshima, southern Japan . These two areas are characterized by a humid subtropical climate. Culicoides oxystoma was collected from May to October in the south of Korea  and from May to November in southern Japan , i.e. in the warmer and rainy period of the year, showing that rainfall and temperature appear to influence the distribution of this species in both Japan and in Senegal.
Livestock density seems to have a weak influence on C. oxystoma abundances. Generally, livestock density was not identified as a variable with a high impact on Culicoides abundance except for C. bolitinos. This is surprising because for many other Culicoides species, especially those vectors of pathogens to animals, the presence of hosts and their density are key factors known to affect their catch [58,59,60]. Several reasons might explain this outcome. First, host densities data was available at a low spatial resolution (100 ×100 km). Secondly, these consisted in 2005 predictions, and densities could have changed (increased or decreased) since then. Thirdly, host densities included only three ruminant species (cattle, sheep and goats) but did not include horse densities. To our knowledge there is no map of horse densities available for Senegal. The lack of data on horse densities could be an important hindrance since in the field, in Senegal, Fall et al.  found high abundances of C. oxystoma on horses. Finally, some Culicoides species are able to feed on other hosts, including birds and wildlife, and it is possible that the densities of these other host species impact the distribution of that species.
As for C. oxystoma, little information in the literature is available on climatic and environmental variables that could impact C. enderleini. In this study, C. enderleini was very abundant in the southern third of the country. RF models showed that its abundance is significantly associated with rainfall, vegetation index, day and night temperatures and, to a lesser extent, livestock density. This is coherent since in the south of Senegal, rainfall is important, average NDVI, temperatures and livestock density are high.
The distribution of C. imicola and C. enderleini followed a similar pattern with an increasing north-south gradient. Again, average rainfall could be a limiting factor in the north, explaining higher abundances in the south (Ziguinchor, Kolda and Kedougou) and center-west (in the groundnut basin) regions that record the highest rainfall in Senegal. The importance of rainfall on the distribution and abundance of C. imicola was highlighted in Morocco  and Europe [22, 26, 37, 62]. The study of Culicoides dynamics in the Niayes area also showed that peaks of abundance of C. imicola occurred during the end of the rainy season (September-October) .
Maximum and minimum NDVI were the second and third most important variables associated with C. imicola abundance. The link between the distribution of C. imicola and NDVI is consistent with previous studies showing that high NDVI values were associated with increased C. imicola abundance in the Niayes area in Senegal , in north Africa and in Europe [16, 24, 30,31,32,33]. In Senegal, areas where the vegetation index is the highest are also those where rainfall is highest, as both are linked . Studies have also highlighted the link between NDVI and soil moisture . Soil moisture could be an important factor influencing larval development of this species. Temperatures (day and night) were also found associated with C. imicola abundance. This finding has been confirmed by several other studies [16, 18, 19, 25, 29, 32, 35].
Concerning C. miombo, maximum NDVI was by far the most important variable influencing its abundance. Its highest abundances were recorded in areas where the vegetation index is the highest, i.e. in the groundnut basin and in southern Senegal. In the north of the country where the vegetation index is the lowest, C. miombo is almost absent. In the Niayes area (particularly in Dakar and Thies), the abundance of C. miombo was low, confirming the results of a previous study . This is probably due to the low vegetation index in this area. In this study, rainfall and temperatures were found associated with the abundance of C. miombo, in coherence with Meiswinkel’s study, which revealed that this species is sensitive to rainfall and high temperatures .
Regarding C. bolitinos, no single variable was a dominant driver of its abundance. Indeed, the contribution of the best environmental predictor (maximum NDVI) on its abundance remained very low (variable importance = 6.31). The most important variables impacting its abundance were the NDVI, livestock density and temperature. These three types of variables have also been linked to C. bolitinos abundance in South Africa . Livestock density was an important variable which greatly impacted model accuracy. Studies conducted by Fall et al.  showed that C. bolitinos is particularly aggressive on horses with a blood-feeding rate of over 75%. This species is also known to feed on large mammals . Other studies describe the larval habitats of C. bolitinos as being the dung of wild and domestic Bovidae in Africa . The ecology of this species is therefore closely linked to its hosts (through feeding and through its breeding media) and it is thus coherent to find a strong link to livestock densities and weaker environmental drivers of its abundance.
Overall, the highest abundances of these five species, proven or suspected vectors of AHS and BT viruses, were recorded in southern Senegal. If C. oxystoma seems present rather in the southwest of the country, C. enderleini, C. miombo and C. imicola have very high abundances in the center-south (Kolda), in the south-east (Kedougou) and in the groundnut basin. Culicoides abundance maps from RF models were compared with those of AHS outbreaks in 2007 in Senegal , which show significant mortality of horses in the groundnut basin. According to Diouf et al. , the greatest risk of introduction of AHS virus in Senegal is through the northeastern border because of important commercial trade movements and numerous markets to and from which carters transport goods. These authors suggested focusing surveillance of AHS virus introduction in that area in order to prevent the virus from reaching the groundnut basin where the dense network of markets could substantially amplify disease transmission and diffusion. Diouf et al.  considered the south and southeastern part of the country to be less at risk of AHS because of lower horse densities due to the presence of tsetse flies. Because low densities of Culicoides were found in the north-eastern part of Senegal, this study strengthens the hypothesis that the main driver of AHS introduction and spread in 2007 was horse movements. It confirms that the groundnut basin is an area at high risk of AHS transmission, because it combines important vector, host and market densities, which could lead to important epidemics. The very important densities of Culicoides in the south suggest that if the virus was present in the neighboring countries (Mali, the Republic of Guinea and Guinea Bissau), it could be introduced through the spread of infected Culicoides and then be amplified in donkeys.
The two main knowledge gaps to better assess AHS transmission risk in Senegal are (i) vector competence of suspected vectors and of very abundant species; and (ii) a map of equid (horse and donkey) densities. Assessing the vector competence for ASHV of species such as C. oxystoma, C. enderleini and C. miombo is essential to evaluate their roles in disease transmission. Mapping not only equid densities but also equid movement networks in Senegal would be of great help for decision makers, since combining vectors abundance maps with maps of equid densities and movements of equids would enable evaluation of AHS transmission risk. In a second step, combining transmission risk with AHS virus introduction risk would enable making recommendations in terms of early warning systems and vaccination policies.
To our knowledge, this study is the first to map the distribution of five major species of Culicoides of veterinary interest in Senegal. The highest abundances of Culicoides were observed in the south and in the groundnut basin, whereas abundances were the lowest in northern Senegal. Abundance maps were produced using three different modelling approaches. RF models provided better estimates of abundances than GLM models and were not limited by interpolation distances contrary to kriging. Environmental and climatic variables of importance that influence the spatial distribution of species abundance were identified. It is now crucial to evaluate the vector competence of major species and then combine the vector densities with densities of horses and ruminants to quantify the risk of transmission of AHS and BT virus across the country.
African horse sickness virus
Day-time land surface temperature
Generalized Linear Models
Normalized difference vegetation index
Night-time land surface temperature
Bessell PR, Searle KR, Auty HK, Handel IG, Purse BV, Bronsvoort BM. Assessing the potential for Bluetongue virus 8 to spread and vaccination strategies in Scotland. Sci Rep. 2016;6:38940.
Carpenter S, Mellor PS, Fall AG, Garros C, Venter GJ. African horse sickness virus: history, transmission, and current status. Annu Rev Entomol. 2017;62:343–58.
Mellor PS, Hamblin C. African horse sickness. Vet Res. 2004;35:445–66.
Maclachlan NJ, Guthrie AJ. Re-emergence of bluetongue, African horse sickness, and other Orbivirus diseases. Vet Res. 2010;41:35.
Tago D, Hammitt JK, Thomas A, Raboisson D. Cost assessment of the movement restriction policy in France during the 2006 bluetongue virus episode (BTV-8). Prev Vet Med. 2014;117:577–89.
Velthuis AGJ, Saatkamp HW, Mourits MCM, de Koeijer AA, Elbers ARW. Financial consequences of the Dutch bluetongue serotype 8 epidemics of 2006 and 2007. Prev Vet Med. 2010;93:294–304.
Akakpo AJ, Wombou Toukam CM, Mankor A, Ly C. Impact économique de l’épizootie de peste équine de 2007 au Sénégal. Bull Anim Hlth Prod Afr. 2011;59:1–16.
Diouf ND, Etter E, Lo MM, Lo M, Akakpo AJ. Outbreaks of African horse sickness in Senegal and methods of control of the 2007 epidemic. Vet Rec. 2013;172:152.
OIE. OIE-Listed diseases, infections and infestations in force in 2017. 2017. http://www.oie.int/en/animal-health-in-the-world/oie-listed-diseases-2017/. Accessed 21 June 2017.
Bakhoum MT, Fall AG, Fall M, Bassene CK, Baldet T, Seck MT, et al. Insight on the larval habitat of Afrotropical Culicoides Latreille (Diptera: Ceratopogonidae) in the Niayes area of Senegal, West Africa. Parasit Vectors. 2016;9:462.
Bakhoum MT, Fall M, Fall AG, Bellis GA, Gottlieb Y, Labuschagne K, et al. First record of Culicoides oxystoma Kieffer and diversity of species within the Schultzei group of Culicoides Latreille (Diptera: Ceratopogonidae) biting midges in Senegal. PLoS One. 2013;8:e84316.
Bakhoum MT, Fall M, Seck MT, Gardes L, Fall AG, Diop M, et al. Foraging range of arthropods with veterinary interest: new insights for Afrotropical Culicoides biting midges (Diptera: Ceratopogonidae) using the ring method. Acta Trop. 2016;157:59–67.
Diarra M, Fall M, Fall AG, Diop A, Seck MT, Garros C, et al. Seasonal dynamics of Culicoides (Diptera: Ceratopogonidae) biting midges, potential vectors of African horse sickness and bluetongue viruses in the Niayes area of Senegal. Parasit Vectors. 2014;7:147.
Fall M, Diarra M, Fall AG, Balenghien T, Seck MT, Bouyer J, et al. Culicoides (Diptera: Ceratopogonidae) midges, the vectors of African horse sickness virus - a host/vector contact study in the Niayes area of Senegal. Parasit Vectors. 2015;8:39.
Diarra M, Fall M, Lancelot R, Diop A, Fall AG, Dicko A, et al. Modelling the abundances of two major Culicoides (Diptera: Ceratopogonidae) species in the Niayes Area of Senegal. PLoS One. 2015;10:e0131021.
Baylis M, Mellor PS, Wittmann EJ, Rogers DJ. Prediction of areas around the Mediterranean at risk of bluetongue by modelling the distribution of its vector using satellite imaging. Vet Rec. 2001;149:639–43.
Cianci D, Hartemink N, Ibáñez-Justicia A. Modelling the potential spatial distribution of mosquito species using three different techniques. Int J Health Geogr. 2015;14:10.
Conte A, Ippoliti C, Savini L, Goffredo M, Meiswinkel R. Novel environmental factors influencing the distribution and abundance of Culicoides imicola and the Obsoletus Complex in Italy. Vet Ital. 2007;43:571–80.
Purse BV, Baylis M, Tatem AJ, Rogers DJ, Mellor PS, Van Ham M, et al. Predicting the risk of bluetongue through time: climate models of temporal patterns of outbreaks in Israel. Rev Sci Tech. 2004;23:761–75.
Tatem AJ, Baylis M, Mellor PS, Purse BV, Capela R, Pena I, Rogers DJ. Prediction of bluetongue vector distribution in Europe and north Africa using satellite imagery. Vet Microbiol. 2003;97:13–29.
Calistri P, Goffredo M, Caporale V, Meiswinkel R. The distribution of Culicoides imicola in Italy: application and evaluation of current Mediterranean models based on climate. J Vet Med B Infect Dis Vet Public Health. 2003;50:132–8.
Calvete C, Estrada R, Miranda MA, Borras D, Calvo JH, Lucientes J. Modelling the distributions and spatial coincidence of bluetongue vectors Culicoides imicola and the Culicoides obsoletus group throughout the Iberian Peninsula. Med Vet Entomol. 2008;22:124–34.
Conte A, Giovannini A, Savini L, Goffredo M, Calistri P, Meiswinkel R. The effect of climate on the presence of Culicoides imicola in Italy. J Vet Med B Infect Dis Vet Public Health. 2003;50:139–47.
Guis H. Geomatique et epidemiologie: caracterisation des paysages favorables à Culicoides imicola, vecteur de la fievre catarrhale ovine en Corse. These de doctorat es Science de la vie et de la santé. Besançon: Universite de Franche-Comté; 2007.
Wittmann E, Mellor P, Baylis M. Using climate data to map the potential distribution of C ulicoides imicola (Diptera: Ceratopogonidae) in Europe. Rev Sci Tech. 2001;20:731–40.
Ducheyne E, Chueca MAM, Lucientes J, Calvete C, Estrada R, Boender GJ, et al. Abundance modelling of invasive and indigenous Culicoides species in Spain. Geospat Health. 2013;8:241–54.
Peters J, De Baets B, Calvete C, Lucientes J, De Clercq E, Ducheyne E, Verhoest N. Absence reduction in entomological surveillance data to improve niche-based distribution models for Culicoides imicola. Prev Vet Med. 2011;100:15–28.
Fall M, Fall A, Seck M, Bouyer J, Diarra M, Balenghien T, et al. Circadian activity of Culicoides oxystoma (Diptera: Ceratopogonidae), potential vector of bluetongue and African horse sickness viruses in the Niayes area, Senegal. Parasitol Res. 2015;114:3151–8.
Baylis M, el Hasnaoui H, Bouayoune H, Touti J, Mellor PS. The spatial and seasonal distribution of African horse sickness and its potential Culicoides vectors in Morocco. Med Vet Entomol. 1997;11:203–12.
Acevedo P, Ruiz-Fons F, Estrada R, Marquez AL, Miranda MA, Gortazar C, Lucientes J. A broad assessment of factors determining Culicoides imicola abundance: modelling the present and forecasting its future in climate change scenarios. PLoS One. 2010;5:e14236.
Baylis M, Bouayoune H, Touti J, El Hasnaoui H. Use of climatic data and satellite imagery to model the abundance of Culicoides imicola, the vector of African horse sickness virus, in Morocco. Med Vet Entomol. 1998;12:255–66.
Baylis M, Meiswinkel R, Venter GJ. A preliminary attempt to use climate data and satellite imagery to model the abundance and distribution of Culicoides imicola (Diptera: Ceratopogonidae) in southern Africa. J S Afr Vet Assoc. 1999;70:80–9.
Baylis M, Rawlings P. Modelling the distribution and abundance of Culicoides imicola in Morocco and Iberia using climatic data and satellite imagery. Arch Virol Suppl. 1998;14:137–53.
Conte A, Goffredo M, Ippoliti C, Meiswinkel R. Influence of biotic and abiotic factors on the distribution and abundance of Culicoides imicola and the Obsoletus Complex in Italy. Vet Parasitol. 2007;150:333–44.
Ippoliti C, Gilbert M, Vanhuysse S, Goffredo M, Satta G, Wolff E, Conte A. Can landscape metrics help determine the Culicoides imicola distribution in Italy? Geospat Health. 2013;8:267–77.
Purse BV, Falconer D, Sullivan MJ, Carpenter S, Mellor PS, Piertney SB, et al. Impacts of climate, host and landscape factors on Culicoides species in Scotland. Med Vet Entomol. 2011;26:168–77.
Rigot T, Conte A, Goffredo M, Ducheyne E, Hendrickx G, Gilbert M. Predicting the spatio-temporal distribution of Culicoides imicola in Sardinia using a discrete-time population model. Parasit Vectors. 2012;5:270.
R Core Team. R: A language and environment of statistical computing. Vienna: R Foundation for Statistical Computing; 2012. https://www.R-project.org.
Becker RA, Allan RW, Brownrigg R, Minka PT, Deckmyn A. maps: Draw Geographical Maps. R package version 3.2.0. 2017. https://CRAN.R-project.org/package=maps.
Becker RA, Allan RW, Brownrigg R. mapdata: Extra Map Databases. R package version 2.2-6. 2016. https://CRAN.R-project.org/package=mapdata.
Bivand R, Lewin-Koh N. maptools: Tools for Reading and Handling Spatial Objects. R package version 0.9-2. 2017. https://CRAN.R-project.org/package=maptools.
Hijmans RJ. raster: Geographic Data Analysis and Modeling. R package version 2.6-7. 2017. https://CRAN.R-project.org/package=raster.
Bivand R, Keitt T, Rowlingson B. rgdal: Bindings for the Geospatial Data Abstraction Library. R package version 1.1-10. 2016. https://CRAN.R-project.org/package=rgdal.
Pebesma EJ, Bivand RS. Classes and methods for spatial data in R. R news. 2005. https://CRAN.R-project.org/doc/Rnews/.
Cressie N. Statistics for spatial data. Revised edition. New York: Wiley; 1993.
McCullagh P, Nelder J. Generalized linear models. London: Chapman & Hall; 1989.
Mullay J. Specifications and testing of some modified count data model. J Econ. 1986;33:341–65.
Zuur AF, Ieno EN, Walker NJ, Saveliev AA, Smith GM. Mixed effects models and extensions in ecology with R. Statistics for Biology and Health. New York: Springer; 2009.
Breiman L. Random forests. Mach Learn. 2001;45:5–32.
Liaw A, Wiener M. Classification and regression by randomForest. R News. 2002;2:18–22.
Rowlingson B, Diggle P. splancs: Spatial and Space-Time Point Pattern Analysis. R package version 2.01-40. 2017. https://CRAN.R-project.org/package=splancs.
Pebesma EJ. Multivariable geostatistics in S: the gstat package. Comput Geosci. 2004;30:683–91.
Boorman J, Dipeolu OO. A taxonomic study of adult Nigerian Culicoides Latreille (Diptera: Ceratopogoaidae) species. Occ Publ Ent Soc Nigeria. 1979;22:1–121.
Carter HF. New West African Ceratopogoninae. Ann Trop Med Parasitol. 1919;12:289–300.
Meiswinkel R, Paweska JT. Evidence for a new field Culicoides vector of African horse sickness in South Africa. Prev Vet Med. 2003;60:243–53.
Oem JK, Chung JY, Kwon MS, Kim TK, Lee TU, Bae YC. Abundance of biting midge species (Diptera: Ceratopogonidae, Culicoides spp.) on cattle farms in Korea. J Vet Sci. 2013;14:91–4.
Yanase T, Matsumoto Y, Matsumori Y, Aizawa M, Hirata M, Kato T, et al. Molecular identification of field-collected Culicoides larvae in the southern part of Japan. J Med Entomol. 2013;50:1105–10.
Baylis M, Parkin H, Kreppel K, Carpenter S, Mellor PS, McIntyre KM. Evaluation of housing as a means to protect cattle from Culicoides biting midges, the vectors of bluetongue virus. Med Vet Entomol. 2010;24:38–45.
Garcia-Saenz A, McCarter P, Baylis M. The influence of host number on the attraction of biting midges, Culicoides spp. to light traps. Med Vet Entomol. 2010;25:113–5.
Kluiters G, Sugden D, Guis H, Marie McIntyre K, Labuschagne K, Vilar MJ, Baylis M. Modelling the spatial distribution of Culicoides biting midges at the local scale. J Appl Ecol. 2013;50:232–42.
Fall M, Fall AG, Seck MT, Bouyer J, Diarra M, Lancelot R, et al. Host preferences and circadian rhythm of Culicoides (Diptera: Ceratopogonidae), vectors of African horse sickness and bluetongue viruses in Senegal. Acta Trop. 2015;149:239–45.
Calvete C, Estrada R, Miranda MA, Borras D, Calvo JH, Lucientes J. Ecological correlates of bluetongue virus in Spain: predicted spatial occurrence and its relationship with the observed abundance of the potential Culicoides spp. vector. Vet J. 2009;182:235–43.
Rogers DJ, Randolph SE. Mortality rates and population density of tsetse flies correlated with satellite imagery. Nature. 1991;351:739–41.
Narasimha Rao PV, Venkataratnam L, Krishna Rao PV, Ramana KV, Singarao MN. Relation between root zone soil moisture and normalized difference vegetation index of vegetated fields. Int J Remote Sens. 1993;14:441–9.
Meiswinkel R. Afrotropical Culicoides: C (Avaritia) miombo sp. nov., a widespread species closely allied to C. (A.) imicola Kieffer, 1913 (Diptera: Ceratopogonidae). Onderstepoort J Vet Res. 1991;58:155–70.
Eksteen S, Breetzke GD. Predicting the abundance of African horse sickness vectors in South Africa using GIS and artificial neural networks. S Afr J Sci. 2011;107:20–8.
Meiswinkel R, Braack LEO. African horse sickness epidemiology: five species of Culicoides (Diptera: Ceratopogonidae) collected live behind the ears and at the dung of the African elephant in the Kruger National Park, South Africa. Onderstepoort J Vet Res. 1994;61:155–70.
Meiswinkel R. Afrotropical Culicoides: a redescription of C. (Avaritia) imicola Kieffer, 1913 (Diptera: Ceratopogonidae) with description of the closely allied C. (A.) bolitinos sp. nov. reared from the dung of the African buffalo, blue wildebeest and cattle in South Africa. Onderstepoort J Vet Res. 1989;56:23–39.
We thank the veterinary services of Senegal for helping us organize and carry out the trapping campaign. We thank the horse owners for allowing us access to their holdings. We thank Dr Thomas Balenghien and Dr Claire Garros for the Culicoides identification training session they organized at ISRA. We also thank Dr Claire Garros for her constructive comments on the manuscript.
This study was funded by EU grant FP7-261504 EDENext and is catalogued by the EDENext Steering Committee as EDENext389 (http://www.edenext.eu). The contents of this publication are the sole responsibility of the authors and do not necessarily reflect the views of the European Commission.
Availability of data and materials
Data supporting the conclusions of this article are included within the article and its additional files. Raw data are available from the corresponding author upon request.
Ethics approval and consent to participate
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Culicoides data from Senegal throughout a nation-wide trapping campaign in 2012. (XLSX 15 kb)
Figure S1. Climatic and environmental data on Senegalese territory with a spatial resolution on 1 km2. Abbreviations: Av, Average; Min, Minimum; Max, Maximum; Dlst, Day land surface temperature; Nlst, Night land surface temperature; lc, landcover; NDVI, Normalized difference vegetation index; Host, livestock density. (PDF 308 kb)
About this article
- African horse sickness
- Culicoides vectors
- Environmental and climatic data
- Random forest models
- Generalized Linear Models
- Spatial distribution