Spatial abundance and clustering of Culicoides (Diptera: Ceratopogonidae) on a local scale
© Kirkeby et al.; licensee BioMed Central Ltd. 2013
Received: 24 August 2012
Accepted: 24 January 2013
Published: 22 February 2013
Biting midges, Culicoides, of the Obsoletus group and the Pulicaris group have been involved in recent outbreaks of bluetongue virus and the former was also involved in the Schmallenberg virus outbreak in northern Europe.
For the first time, here we investigate the local abundance pattern of these two species groups in the field by intensive sampling with a grid of light traps on 16 catch nights. Neighboring trap catches can be spatially dependent on each other, hence we developed a conditional autoregressive (CAR) model framework to test a number of spatial and non-spatial covariates expected to affect Culicoides abundance.
The distance to sheep penned in the corner of the study field significantly increased the abundance level up to 200 meters away from the sheep. Spatial clustering was found to be significant but could not be explained by any known factors, and cluster locations shifted between catch nights. No significant temporal autocorrelation was detected. CAR models for both species groups identified a significant positive impact of humidity and significant negative impacts of precipitation and wind turbulence. Temperature was also found to be significant with a peak at just below 16 degrees Celcius. Surprisingly, there was a significant positive impact of wind speed. The CAR model for the Pulicaris group also identified a significant attraction to the smaller groups of sheep placed in the field. Furthermore, a large number of spatial covariates which were incorrectly found to be significant in ordinary regression models were not significant in the CAR models. The 95% C.I. on the prediction estimates ranged from 20.4% to 304.8%, underlining the difficulties of predicting the abundance of Culicoides.
We found that significant spatial clusters of Culicoides moved around in a dynamic pattern varying between catch nights. This conforms with the modeling but was not explained by any of the tested covariates. The mean abundance within these clusters was up to 11 times higher for the Obsoletus group and 4 times higher for the Pulicaris group compared to the rest of the field.
Since the incursion of bluetongue virus into northern Europe and the subsequent discovery of Schmallenberg virus in the same region, Culicoides populations on farms have become important for epidemiological research. Species of the Obsoletus group and the Pulicaris group are suspected to play an important role in north European outbreaks of bluetongue and are found throughout northern Europe [1–5]. Recently, it was confirmed that species in the Obsoletus group can replicate Schmallenberg virus . Many large-scale studies and transmission models have included spatial estimates of the abundance of Culicoides in Europe ([7–16]), but few studies have investigated the spatial pattern of Culicoides abundance on a local scale: In 1951, Kettle  found that the abundance of C. impunctatus decreased proportionally with distance to their breeding sites. This species is not dominant on farms but frequently associated with bogs (e.g. ). Later, Kettle  found indication of higher abundances of C. impunctatus and C. pulicaris L. near hosts (cattle, horses and humans). Garcia-Saenz et al. found a positive correlation between the number of sheep near a light trap, and the number of female C. obsoletus caught in the trap. Rigot et al. found that the abundance of different species of Culicoides were positively correlated with closeness to farms in Belgium. In a large scale study, Purse et al. found that the abundance of adult C. pulicaris sensu stricto was correlated with vegetation indices, land use and elevation above sea level; C. punctatus abundance was correlated with the presence of sheep, temperature, land use and vegetation; and the abundance of C. obsoletus was only correlated with temperature. Also, the abundance of adult C. impunctatus was found to have a negative correlation with the presence of cattle, which might be because of their breeding sites (bogs) that are often located away from cattle. Remote sensing can be used to estimate the abundance of Culicoides (e.g. [8, 18, 22]), but provides only estimates of Culicoides abundance on a rough scale. In this study we take a novel approach, using local-scale abundance data to investigate possible spatial and temporal covariates for prediction of Culicoides abundance within a field.
Neighboring insect traps can be spatially dependent on each other (e.g. [23, 24]), and Rigot et al.  found significant overlapping catching areas between 8 W Onderstepoort traps situated 50 meters apart. Thus it is necessary to take spatial autocorrelation into account. We developed conditional autoregressive (CAR) models for the abundance of two Culicoides vector species groups in order to account for the spatial dependency. For the first time, spatial autocorrelation is incorporated in a prediction model for Culicoides on a local scale, making trap catches spatially independent by including information from neighboring traps. Using this approach, a number of spatial covariates which have a significant impact in ordinary regression modeling, no longer appear significant and some temporal covariates become significant. At the same time we provide a method to deal with a lot of missing data in a spatial dataset by including second order neighbors when first order neighbors are missing. Furthermore, we estimate the spatial autocorrelation between trap catches and demonstrate the need to take it into account by incorporating it into statistical models. Lastly, we examine the abundance pattern not explained by the systematic part of the CAR models through cluster analysis.
The total number of Culicoides caught on each successful catch night and the temporal covariates
Date of sampling night
Obsoletus group total
Obsoletus group min/max
Pulicaris group total
Pulicaris group min/max
Wind speed (m/s)
Mean humidity (%RH)
Date of sampling night
Obsoletus group total
Obsoletus group min/max
Pulicaris group total
Pulicaris group min/max
Wind speed (m/s)
Mean humidity (%RH)
Four potential breeding sites (A-D on Figure 1) for the Pulicaris group were subjectively identified on the field : Site A was a shallow assembly of water without any boundary vegetation and a 1-3 meter broad mud zone; Site B was an old marl pit with shallow water, heavily shaded by dense thicket with trees; Site C was a small pond with reed along the steep edges; site D was a muddy area on the field with small temporary water bodies. Throughout the study period, twenty to fifty (according to area size) mud samples (97 mm in diameter) were taken weekly from each potential breeding site and kept in emergence chambers at an indoor facility (following ) to confirm breeding. Outside the potential breeding sites, an additional 50 soil samples were taken weekly at randomly generated coordinates to screen for Obsoletus group breeding sites and for unexpected breeding sites of the Pulicaris group. We did not target breeding sites of the Obsoletus group as they are poorly investigated. They are associated with factors that are difficult to include in a model such as dung heaps and leaf litter, but this topic is still largely uncovered [28–30].
On three nights, sheep were placed in a transect in the middle of the field to test the attraction effect of a few sheep compared to the flock. In three transect points, two sheep were placed together in a 3 by 3 meter enclosure under a light trap (see Figure 1). The distance between transect points was 150 meters. During the study period, a weather station (Davis Vantage Pro 2) with a data logger was set up to record temperature, precipitation, humidity, wind speed and wind direction at 5 minute intervals. It was placed in the middle of the field to keep away from interfering vegetation. Light traps were emptied at dawn, and the caught Culicoides were preserved in 70% ethanol. The samples were analyzed under a dissection microscope, and sorted to species group and sex following Campbell and Pelham-Clinton . Only females of the Obsoletus group (comprising C. obsoletus, C. scoticus, C. chiopterus and C. dewulfi) and the Pulicaris group (here comprising C. pulicaris and C. punctatus), were included in this analysis. We only considered the two dominant species in the latter group since other members of this group are rare in farm areas (pers. obs.) and not identified as a disease vector in this region. Due to time constraints, on 8 of the catch nights we only counted 50% of the trap catches. On these catch nights (the dates are underlined in Table 1), every second sample, chosen in a checkerboard pattern, was analyzed. All 16 nights were included in the models.
To deal with a high number of low catches we stabilized the observations by transforming the numbers with the natural logarithm prior to analysis, log(x+1). Thus for low numbers the observations will converge towards 1 instead of zero. For simplicity, we here denote the transformation as log(X) in the equations.
Only weather records during the flight periods of Culicoides (assumed to be one hour before to three hours after sunset and two hours before to one hour after sunrise) were used in the analysis because we assume that the trap catches were only directly affected by the weather in this time interval. Mean temperature, humidity and wind speed measurements recorded on the field during the flight periods were included directly as covariates. Precipitation was summed over each flight period and included as a covariate. As an estimate of the wind turbulence, changes in wind direction was defined in steps as a minimum change of wind direction of 22.5 degrees. The highest number of steps that the wind direction changed in either 5 or 10 minute intervals, measured within each flight period, was calculated. As each catch night consisted of two flight periods, the mean of the two highest step change numbers for each flight period was used as the turbulence covariate for each catch night.
Mean, standard error and ranges of spatial covariates in the models
38 - 653
Breeding site effect
4.7 · 10−4
9.6 · 10−7
3.6 · 10−5 - 5.1 ∗ 10−3
0.001 - 1.383
Sheep scent effect
0.000 - 0.041
Ordinary regression modeling
Where the log of the abundance of Culicoides (X) is determined by covariates Z and their coefficients β (where β T signifies the matrix transpose of β), and a residual error term ∈. Model reduction was performed with the likelihood ratio method, and covariates that did not contribute significantly (p ≥ 0.05) were excluded. After model reduction, all excluded covariates were tested again by forward selection, with the test sequence defined through the Akaike Information Criterion (AIC) . These models treated the trap catches as stochastically independent of each other and hence ignored potential spatial autocorrelation (clustering). All regression modeling was carried out in R 2.14.2 (http://www.r-project.org).
To account for the spatial autocorrelation within the dataset, a conditional autoregressive model (CAR) model was constructed for each species group by assuming spatial dependence in the model (1) as described in the following. The model estimation and test procedure is described in Additional file 1.
To evaluate the performance of the CAR models, the models were examined for significant spatial clusters in the residuals using a normal distribution model in SaTScan v. 91.1.1 (http://www.satscan.org). For each catch night and each species group, the model residuals were tested for circular or elliptic hotspots or coldspots, without penalty for elliptic clusters and allowing multiple hotspots. Each scan was run for 9999 iterations, testing for significant clusters at the 5% level. To investigate the spatial autocorrelation pattern not explained by the systematic part of the CAR models, we adjusted the observations for the significant effects found in the CAR models and then tested each catch night for significant clusters in SaTScan. Ordinary regression model and CAR model fit were tested by plotting the distribution of residuals and by quantile-quantile-plots.
No sheep were harmed in this study. Permission to move the sheep was obtained by the owners.
During the study period, successful catches from 16 nights, consisting of 530 trap catches, were included in the analysis. A total of 19,654 female Culicoides were counted: 15,166 from the Pulicaris group and 4,488 from the Obsoletus group (Table 1). The parameter ranges within the active period were: mean temperature: 12.1 – 19.9 degrees Celsius, mean wind speed: 0.08 – 5.47 m/s, precipitation: 0-3.6 mm, Relative humidity: 54-100%. The distance from each trap to the sheep was 38 – 653 meters and the distance to the nearest breeding sites 1-45 meters. Catches were excluded if the sheep broke through the enclosure during the night; the trap was damaged or not operating properly.
A total number of 208 Culicoides spp. hatched from the emergence chambers, of which 16 were from the Pulicaris group and none were from the Obsoletus group. The other species that hatched were mostly C. pictipennis and C. festivipennis. From breeding site A, 24 Culicoides spp. (none were from the Pulicaris group) emerged from 350 soil samples. From the shaded breeding site B (Figure 1) no Culicoides but many Psychodidae spp. emerged from 140 soil samples. From breeding site C, 152 Culicoides spp. (of which 13 were from the Pulicaris group) emerged from 140 samples. From breeding site D, 32 Culicoides spp. (of which 3 were from the Pulicaris group) emerged from 140 samples. No Culicoides emerged from the 350 random samples on the field, indicating that Pulicaris group breeding sites were confined to the identified breeding sites and that the Obsoletus group did not emerge on the field during the study period. Distance to breeding sites A, C and D was included in the modeling procedure as they were found to be breeding sites for Culicoides.
Significant coefficients from the models
Distance to sheep
Distance to sheep 2
Wind speed 2
Windbreaks * Sheep scent
Wind speed * Sheep scent
Catch night 21.07
Catch night 27.07
Catch night 03.08
Catch night 04.08
Catch night 06.08
Catch night 17.08
Catch night 18.08
Catch night 21.08
Catch night 24.08
Catch night 25.08
Catch night 27.08
Catch night 28.08
Catch night 31.08
Catch night 03.09
Catch night 04.09
The ordinary regression models without spatial autocorrelation identified more significant spatial covariates than the CAR models did, and the CAR models identified more temporal covariates than the ordinary regression models (Table 3).
For both species groups we found a significant positive effect of humidity and a significant negative effect of turbulence and precipitation (Table 3). There was also a significant effect of temperature and wind speed including their squared terms. The temperature effect showed peak abundance at just below 16 degrees Celcius and the wind speed surprisingly showed a positive effect with increasing wind speed between 1 and 5 m/s (Figure 3).
Only the CAR model for the Pulicaris group identified a significant effect of the transect of sheep on three catch nights. The effect of the inverse Euclidean distance to the small enclosures with pairs of sheep is positive, meaning that the Pulicaris group abundance is higher close to the pairs of sheep.
We tested for clusters in the residuals of the CAR models to check if the spatial autocorrelation was fully extracted in the models. In the residuals of the Obsoletus group CAR model we found two clusters (p=0.0045 and 0.0006) on the nights of 28.08 and 31.08. We also found two clusters (p=0.0375 and 0.0427) in the Pulicaris group CAR model on the nights of 06.08 and 31.08.
We tested the observations for clustering without the effect of host animals to investigate the spatial clustering pattern not explained by the systematic covariates in the CAR models. It revealed a dynamic pattern with higher Culicoides abundance in different places, varying between catch nights, so clusters were moving around on the field (Figures 4 and 5). This is consistent with the CAR modeling, and implies that one or more yet unidentified factors influenced the Culicoides abundance in a spatial pattern that changes each night. The ratios of the significant hotspots show that the mean abundance of the Obsoletus group in a significant hotspot was 0.62-10.82 times higher than the rest of the field (Figure 4), and 1.75-4.16 times higher for the Pulicaris group (Figure 5). This result is striking and can seriously impact field studies of Culicoides abundance. Since no known factor could explain this dynamic pattern, it will cause noise in abundance studies. The best way to to take account for this is to conduct large-scale studies with many traps and locations reducing the noise from the spatial clustering. It is not possible to obtain a reliable measure of the level of abundance in an area by using a single trap. However, this does not mean that national or regional scale predictive abundance models are invalid if they are based on just one trap per farm. If a large number of farms are sampled, the general relationship between environmental factors and the mean abundance can still be quantified. Such models may therefore be able to predict a mean trap collection on farms associated with a specific combination of environmental covariates (e.g. [8, 11, 12]). But if the same models are used to predict catch sizes in a trap at a specific farm it may result in very large residuals as a result of the large spatial variation in abundance on the same farm.
In this study we found that the spatial autocorrelation between traps was highly significant. This means that if a trap catches more than expected, another trap close by is also likely to catch more than expected. For the Obsoletus group the spatial autocorrelation was 0.41. We explored this further (using equation (5), see Additional file 1), to have a look at the relation between two traps, A and B, with 50 m distance. For an expected level of abundance at 100% in both traps, if trap A catches 20% more than expected, then trap B will be expected to catch 7.8% more. If trap A catches 50% more than expected, trap B will be expected to catch 18.1% more. For the Pulicaris group the spatial autocorrelation was 0.235. Using the same scenarios, trap B would catch 4.4% and 10.0% more than expected, respectively. The spatial autocorrelation means that traps placed close to each other do not provide independent estimates of abundance. The true variance in abundance will therefore be underestimated unless traps are widely separated. This has to be taken into account when using more than a single trap at a site.
The CAR models should extract the spatial clustering from the data and therefore leave no significant clusters in the residuals. However, we found two clusters (p=0.0045 and 0.0006) in the residuals of the Obsoletus group CAR model. The first cluster is on the night of the 28.08.2009 where only one specimen from the Obsoletus group was caught on the entire field, and thus we ascribe this cluster as an artefact. The second cluster in the Obsoletus group CAR model residuals on the night of the 31.08.2009 is also highly significant. We performed the parameter estimation again without this catch night and found similar estimates of the effects (data not shown). Thus we conclude that this model violation does not influence the general validity of the model. Two clusters were found in the Pulicaris group CAR model residuals with p-values only just below the significance level (p=0.0375 and 0.0427). Therefore we do not doubt the general applicability of this model either. Furthermore, the SaTScan analysis used to detect clusters is not able to deal with the varying variance included in the model residuals created by differing neighbor configurations, making this test very rigid.
From the two CAR models, the residual variance was estimated to be 0.69 and 0.65 (Table 3). We can use this variance to estimate the general 95% confidence intervals of the abundance estimates. Thus the 95% interval for Obsoletus group CAR model ranged from 20.4% to 304.8% of the predicted catch size. For the Pulicaris group CAR model the interval ranged from 22.6% to 289.4%. This highlights the huge variation in the catches. Estimates of vector abundance based on single traps are expected to vary dramatically depending on the exact position chosen for the trap. This high uncertainty associated with abundance estimates based on single traps needs to be taken into account when modeling the abundance of Culicoides on a greater scale and in simulation models of vectorborne disease that rely on vector abundance estimates.
The estimates of the significant effects in the models are fairly similar between the two species groups (Figure 3, Table 3). This supports the results of the models and indicates that the effects found may be general for species of Culicoides. Especially the significant temporal covariates, which may be general for Culicoides because they are not influenced by host preferences.
The dynamic pattern is also fairly similar between the two species groups. Surprisingly, three of the significant hotspots for the Obsoletus group and two for the Pulicaris group were found in the southern part of the field, away from the sheep. A possible explanation for this is swarming behavior. Downes observed in 1955  that different species of Culicoides swarm above certain markers such as cow dung, a dark cloth or other conspicuous objects. Both the Obsoletus group and the Pulicaris group have been observed swarming, and it is likely that swarming can blur the general abundance pattern. Very few males were caught in the light traps in this study, and they seemed to be correlated with high female abundance (data not shown), which could also indicate swarming behaviour.
Similar to the results from other studies [20, 21, 36], we found a significant effect of the vicinity of host animals for both the Obsoletus group and the Pulicaris group. In a study of Calvete et al., traps were placed within 30 m from each farm to obtain estimates of the abundance of Culicoides, and Goffredo and Meiswinkel  pointed out that when monitoring Culicoides, light traps should be placed in the near vicinity of vertebrate hosts. This is supported by the present study where we quantified the effect of host animals. We found that traps placed near host animals increased the overall vector abundance with approximately 50% - 100% compared with 300-400 m away from the host animals. However, we also found an increased level of abundance for both species groups in the southern part of the field. This could be an artefact in the simple two-parameter model construction, or it could indicate a depletion of Culicoides abundance around the host animals. In the latter case, the abundance level is normal again at 650 m distance from the host animals. An alternative explanation could be that this effect is caused by the small forest area in the southern part of the field.
This pattern is relevant for other studies of the abundance of Culicoides. Traditionally, Culicoides monitoring programmes are carried out running a single trap on each farm near host animals. Calvete et al. mentions that traps were placed within 30 m from the hosts to ensure a high catch. Goffredo and Meiswinkel  suggest that traps are placed in the vicinity of hosts for monitoring programmes. We suggest, that the trap placement should be standardized or adjusted with regards to the distance to host animals because the distance to the hosts impacts directly on the trap catch. For instance, if placement of the traps just next to the host animals is impossible, all traps in a study should be placed at the same distance to obtain comparable measures at different farms. Alternatively, if one trap is placed sub-optimally at for instance 300 meters distance from the host animals, catches of the Pulicaris group made here should be adjusted up by 150%.
We also found a significant effect of the sheep placed in transects on the field for the Pulicaris group. This emphasizes that this species group is more abundant where the host animals are, and that even two sheep can have an impact on the abundance of this species group as found by Garcia-Saenz et al.. It also underlines the fact that Culicoides can find any small group of host animals regardless of other groups of hosts nearby, which makes them very efficient disease vectors.
The temperature was significant for both species groups with peak abundance at 16 degrees Celcius and no effect below 14 degrees or above 18 degrees (Figure 3). This is in concordance with Conte et al. who found that the minimum temperature for activity of the Obsoletus Complex was 14.2 (13.9–14.6) degrees Celcius. Garcia-Saenz et al. found no significant effect of temperature on the abundance of Culicoides, but Carpenter et al. found a peak biting rate at 21 degrees. The latter study included catches at temperature up to 29 degrees, which was not possible to include in the present study.
The humidity was found to have a positive significant effect on the abundance of both species groups. This is in concordance with Carpenter et al. who found a positive correlation between humidity and Culicoides abundance. Carpenter et al. and Baylis et al. also found a positive effect of humidity on the abundance of the Obsoletus group. Turbulence had a significant negative effect in the CAR models for both species groups. Carpenter et al. also found this significant effect. In the present study we also found that precipitation had a negative effect on the Culicoides abundance. This contrasts with the findings of Blackwell  who found a positive effect of rain on catches of C. impunctatus.
We found a significant effect of wind speed and its quadratic term for both species groups. When plotting with confidence intervals, the abundance increases with the wind speed in the investigated interval (Figure 3). In contrast, Carpenter et al. found decreasing abundance for wind speeds exceeding 3 m/s. A possible explanation of the findings in our study is that if the wind is weak and the Culicoides therefore have difficulties in determining the direction of hosts by scent, they are reluctant to waste energy on flying. Thus, within the investigated range of windspeed, higher wind speeds yield a higher abundance of active Culicoides.
No Culicoides emerged from breeding site B (Figure 1). This could be due to the thicket and trees shading the pond, which prevents the sun from heating up the mud to the necessary temperature for Culicoides to breed. The other three sunlit breeding sites were expected as breeding sites for Culicoides spp. In this study we used light traps to measure Culicoides abundance. Therefore the results may be influenced by bias of the trapping method such as variation in attraction for different species and for different lifestages of Culicoides[38, 41–43]. Future trapping studies should ideally distinguish specimens to species level in order to determine the differences in the behaviour between species with regards to light traps.
The spatial autocorrelation, ρ, was found significant, meaning that it is necessary to take spatial clustering into account on this scale. Even on a larger scale, spatial clustering is important to incorporate in the modeling framework as shown by . The temporal autocorrelation, θ, was found non-significant. This was expected since the intervals between catch nights ranged from 0 to 10 nights. The ordinary regression models identified more significant spatial covariates than the CAR models, effects which the CAR models discarded through the inclusion of local dependence given by the spatial correlation (Table 3). A possible explanation for the extra significant spatial covariates included in the ordinary regression models is that they compensate for the spatial clustering by including more explanatory covariates, and it should be noted that given the validity of the CAR model, these significances are type 1 errors, ie. false significances. This interpretation is further supported by the fact that the significant covariates shared by the CAR models and the ordinary regression models are fairly alike (Table 3). In the present study, the systematic effect of each catch night may have overtaken the effect of some of the covariates when few catch nights are sampled because non-spatial covariates will covary with catch night, which is a drawback of this type of model. However, the advantage is that we obtain more precise estimates of significant covariates corrected for the effect of spatial autocorrelation.
We used Besag’s block design to build the CAR models in this study . Formulating the spatial autocorrelation as an exponentially decreasing correlation between neighboring traps we were able to include data points where all first order neighbors were missing by taking second order neighbors into account. This approach is useful in studies of grid measurements where many missing data are present.
The spatial autocorrelation between trap catches, ρ, accounts for other potential unknown covariates which were not spatially consistent between catch nights. However, if an unknown, spatially fixed factor influenced the abundance of Culicoides, the temporal autocorrelation, θ, would tend to be significant, indicating that some traps consistently caught higher numbers of Culicoides. But since the temporal autocorrelation was found insignificant and the spatial autocorrelation was found significant, there is no evidence for the presence of unknown spatially fixed covariates.
We revealed a spatially varying pattern of abundance that varies between catch nights, where unpredictable hotspots caused the mean trap catch to be up to 11 times higher for the Obsoletus group and 4 times higher for the Pulicaris group. From the residual variance of the models we calculated that the 95% C.I. on the prediction of abundance is approximately 20% to 300%, which is important to consider when conducting large-scale studies. We found no significant spatial covariates determining the abundance of the studied species groups other than the distance to host animals and for the Pulicaris group this also included pairs of sheep placed in small enclosures on the field. Thus no low risk areas for placing host animals susceptible to bluetongue or Schmallenberg virus were identified on this scale because the abundance of Culicoides was indeed determined by the presence of host animals. We have demonstrated the importance of placing traps near the hosts when monitoring Culicoides, as we see a significantly increased abundance of Culicoides (up to 100%) in a radius of approximately 200 meters from the hosts. We also found significant positive effects of humidity and wind speed, significant negative effects of precipitation and turbulence. The optimum temperature for abundance of both species groups was found to be just below 16 degrees Celcius.
Thanks to Frank Hansen and Rune Ploug from Vallø Lam for their invaluable help with experimental setup, Simon Haarder and Peter Iversen for help with species determination and Kaare Græsbøll for commenting on the script. This study was partially funded by the Danish Ministry of Food, Agriculture and Fisheries and by EU grant GOCE-2003-010284 EDENext and is catalogued by the EDENext Steering Committee as EDENext 065 (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.
- Caracappa S, Torina A, Guercio A, Vitale F, Calabro A, Purpari G, Ferrantelli V, Vitale M, Mellor P: Identification of a novel bluetongue virus vector species of Culicoides in Sicily. Vet Rec. 2003, 153 (3): 71-74. 10.1136/vr.153.3.71.View ArticlePubMedGoogle Scholar
- De Liberato C, Scavia G, Lorenzetti R, Scaramozzino P, Amaddeo D, Cardeti G, Scicluna M, Ferrari G, Autorino G: Identification of Culicoides obsoletus (Diptera: Ceratopogonidae) as a vector of bluetongue virus in central Italy. Vet Rec. 2005, 156 (10): 301-304.View ArticlePubMedGoogle Scholar
- Carpenter S, Lunt H, Arav D, Venter G, Mellor P: Oral susceptibility to bluetongue virus of Culicoides (Diptera: Ceratopogonidae) from the United Kingdom. J Med Entomol. 2006, 43: 73-78. 10.1603/0022-2585(2006)043[0073:OSTBVO]2.0.CO;2.View ArticlePubMedGoogle Scholar
- Carpenter S, Wilson A, Mellor P: Culicoides and the emergence of bluetongue virus in northern Europe. Trends Microbiol. 2009, 17 (4): 172-178. 10.1016/j.tim.2009.01.001.View ArticlePubMedGoogle Scholar
- Nielsen S, Nielsen B, Chirico J: Monitoring of biting midges (Diptera: Ceratopogonidae: Culicoides Latreille) on farms in Sweden during the emergence of the 2008 epidemic of bluetongue. Parasitol Res. 2010, 106 (5): 1197-1203. 10.1007/s00436-010-1791-4.View ArticlePubMedGoogle Scholar
- Rasmussen L, Kristensen B, Kirkeby C, Rasmussen T, Belsham G, Bødker R, Bøtner A: Culicoids as vectors of schmallenberg virus. Emerg Infect Dis. 2012, 18 (7): 1204-1206.PubMed CentralPubMedGoogle Scholar
- Rawlings P, Pro M, Pena I, Ortega M, Capela R: Spatial and seasonal distribution of Culicoides imicola in Iberia in relation to the transmission of African horse sickness virus. Med Vet Entomol. 1997, 11: 49-57. 10.1111/j.1365-2915.1997.tb00289.x.View ArticlePubMedGoogle Scholar
- Purse B, Tatem A, Caracappa S, Rogers D, Mellor P, Baylis M, Torina A: Modelling the distributions of Culicoides bluetongue virus vectors in Sicily in relation to satellite-derived climate variables. Med Vet Entomol. 2004, 18 (2): 90-101. 10.1111/j.0269-283X.2004.00492.x.View ArticlePubMedGoogle Scholar
- Scavia G, Autorino G, De Liberato C, Farina F, Ferrari G, Guidoni M, Magliano A, Miceli M, Scholl F, Scicluna M: Association between the 2001-2003 bluetongue epidemic in Lazio and Tuscany (central Italy) and distribution and abundance of Culicoides imicola and C. obsoletus vectors. Vet Ital. 2004, 40 (3): 274-PubMedGoogle Scholar
- Acevedo P, Ruiz-Fons F, Estrada R, Márquez A, Miranda M, Gortázar 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 (12): e14236-10.1371/journal.pone.0014236.PubMed CentralView ArticlePubMedGoogle Scholar
- Calvete C, Estrada R, Miranda M, Borrás D, Calvo J, 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 (2): 235-243. 10.1016/j.tvjl.2008.06.010.View ArticlePubMedGoogle Scholar
- 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 (4): 333-344. 10.1016/j.vetpar.2007.09.021.View ArticlePubMedGoogle Scholar
- Guis H, Labuschagne K, Sugden D, McIntire J, Vilar M, Baylis M: Spatial variability in catches of Culicoides across farms in the Bala area of North Wales. Revue d’élevage et de médecine vétérinaire des pays tropicaux 2009: 2008 Medreonet Proceedings. 2008, 62 (2-4): 103-Google Scholar
- Racloz V, Venter G, Griot C, Stärk K: Estimating the temporal and spatial risk of bluetongue related to the incursion of infected vectors into Switzerland. BMC Vet Res. 2008, 4: 42-10.1186/1746-6148-4-42.PubMed CentralView ArticlePubMedGoogle Scholar
- Sedda L, Brown H, Purse B, Burgin L, Gloster J, Rogers D: A new algorithm quantifies the roles of wind and midge flight activity in the bluetongue epizootic in northwest Europe. Proc R Soc B: Biol Sci. 2012, 279 (1737): 2354-2362. 10.1098/rspb.2011.2555.View ArticleGoogle Scholar
- Searle K, Blackwell A, Falconer D, Sullivan M, Butler A, Purse B: Identifying environmental drivers of insect phenology across space and time: Culicoides in Scotland as a case study. Bull Entomol Res. 2012, 1: 1-16.Google Scholar
- Kettle D: The spatial distribution of Culicoides impunctatus Goet. under woodland and moorland conditions and its flight range through woodland. Bull Entomol Res. 1951, 42 (02): 239-291. 10.1017/S0007485300025323.View ArticleGoogle Scholar
- Purse B, Falconer D, Sullivan M, Carpenter S, Mellor P, Piertney S, Albon S, Gunn G, Blackwell A: Impacts of climate, host and landscape factors on Culicoides species in Scotland. Med Vet Entomol. 2011, 26 (2): 168-177.View ArticlePubMedGoogle Scholar
- Kettle D: The flight of Culicoides impunctatus Goetghebuer (Diptera, ceratopogonidae) over moorland and its bearing on midge control. Bull Entomol Res. 1960, 51 (03): 461-489. 10.1017/S0007485300055115.View ArticleGoogle Scholar
- 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. 2011, 25: 113-115. 10.1111/j.1365-2915.2010.00904.x.View ArticlePubMedGoogle Scholar
- Rigot T, Vercauteren Drubbel M, Delécolle J, Gilbert M: Farms, pastures and woodlands: the fine-scale distribution of Palearctic Culicoides spp. biting midges along an agro-ecological gradient. Med Vet Entomol. 2012, [http://dx.doi.org/10.1111/j.1365--2915.2012.01032.x ]Google Scholar
- Guis H, Tran A, de La Rocque S, Baldet T, Gerbier G, Biteau-Coroller F, Roger F, Viel J, Mauny F, Barragué B: Use of high spatial resolution satellite imagery to characterize landscapes at risk for bluetongue. Vet Res. 2007, 38 (5): 669-683. 10.1051/vetres:2007025.View ArticlePubMedGoogle Scholar
- Papadopoulos N, Katsoyannos B, Nestle D: Spatial autocorrelation analysis of a Ceratitis capitata (Diptera: Tephritidae) adult population in a mixed deciduous fruit orchard in northern Greece. Environ Entomol. 2003, 32 (2): 319-326. 10.1603/0046-225X-32.2.319.View ArticleGoogle Scholar
- Cocu N, Harrington R, Rounsevell M, Hullé M: Spatial autocorrelation as a tool for identifying the geographical patterns of aphid annual abundance. Agric Forest Entomol. 2005, 7: 31-43. 10.1111/j.1461-9555.2005.00245.x.View ArticleGoogle Scholar
- Rigot T, Gilbert M: Quantifying the spatial dependence of Culicoides midge samples collected by Onderstepoort-type blacklight traps: an experimental approach to infer the range of attraction of light traps. Med Vet Entomol. 2011, 26 (2): 152-161.View ArticlePubMedGoogle Scholar
- Nielsen S: Forurening og formering af mitter i Egeløkke Lung: larvetæthed og aktiviteter af voksne mitter. 1998, Denmark: Miljø-og Energiministeriet, MiljøstyrelsenGoogle Scholar
- Kirkeby C, Bødker R, Stockmarr A, Enøe C: Association between land cover and Culicoides (Diptera: Ceratopogonidae) breeding sites on four Danish cattlefarms. Entomologica Fennica. 2009, 20 (4): 228-232.Google Scholar
- Zimmer J, Haubruge E, Francis F, Bortels J, Simonon G, Losson B, Mignon B, Paternostre J, De Deken R, De Deken G: Breeding sites of bluetongue vectors in northern Europe. Vet Rec. 2008, 162 (4): 131-131.View ArticlePubMedGoogle Scholar
- Zimmer J, Saegerman C, Losson B, Haubruge E: Breeding sites of bluetongue virus vectors, Belgium. Emerg Infect Dis. 2010, 16 (3): 575-10.3201/eid1603.091311.PubMed CentralView ArticlePubMedGoogle Scholar
- Ninio C, Augot D, Dufour B, Depaquit J: Emergence of Culicoides obsoletus from indoor and outdoor breeding sites. Vet Parasitol. 2011, 183 (1-2): 125-10.1016/j.vetpar.2011.07.020.View ArticlePubMedGoogle Scholar
- Campbell J, Pelham-Clinton E: X.A taxonomic review of the british species of Culicoides Latreille (Diptera, Ceratopogonidæ). Proc R Soc Edinburgh. Section B. Biol. 1960, 67 (03): 181-302. 10.1017/S0080455X00000758.View ArticleGoogle Scholar
- Akaike H: Information theory and an extension of the maximum likelihood principle. Second International Symposium on Information Theory, Volume 1. Edited by: Petrov BN, Csaki F. 1973, Budapest, Akademiai Kiado: Springer Verlag, 267-281.Google Scholar
- Besag J: Spatial interaction and the statistical analysis of lattice systems. J R Stat Soc Ser B (Methodological). 1974, 36 (2): 192-236.Google Scholar
- Fisher R: On the mathematical foundations of theoretical statistics. Proc R Soc Lond B Biol Sci. 1922, 222 (594-604): 309-368.Google Scholar
- Downes J: Observations on the swarming flight and mating of Culicoides (Diptera: Ceratopogonidae) 1. Trans Royal Entomol Soc Lond. 1955, 106 (5): 213-236.View ArticleGoogle Scholar
- Bellis G, Reid D: Sampling bias in determining the parous rate of collections of Culicoides brevitarsis Kieffer and C. wadai Kitaoka (Diptera: Ceratopogonidae). Aust J Entomol. 1996, 35 (4): 319-322. 10.1111/j.1440-6055.1996.tb01411.x.View ArticleGoogle Scholar
- Goffredo M, Meiswinkel R: Entomological surveillance of bluetongue in Italy: methods of capture, catch analysis and identification of Culicoides biting midges. Vet Ital. 2004, 40 (3): 260-265.PubMedGoogle Scholar
- Carpenter S, Szmaragd C, Barber J, Labuschagne K, Gubbins S, Mellor P: An assessment of Culicoides surveillance techniques in northern Europe: have we underestimated a potential bluetongue virus vector?. J Appl Ecol. 2008, 45 (4): 1237-1245.Google Scholar
- Baylis M, Parkin H, Kreppel K, Carpenter S, Mellor P, McIntyre K: 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. 10.1111/j.1365-2915.2009.00842.x.View ArticlePubMedGoogle Scholar
- Blackwell A: Diel flight periodicity of the biting midge Culicoides impunctatus and the effects of meteorological conditions. Med Vet Entomol. 1997, 11 (4): 361-367. 10.1111/j.1365-2915.1997.tb00423.x.View ArticlePubMedGoogle Scholar
- Gerry A, Monteys V, Vidal J, Francino O, Mullens B: Biting rates of Culicoides midges (Diptera: Ceratopogonidae) on sheep in northeastern Spain in relation to midge capture using UV light and carbon dioxide-baited traps. J Med Entomol. 2009, 46 (3): 615-624. 10.1603/033.046.0329.View ArticlePubMedGoogle Scholar
- Viennet E, Garros C, Lancelot R, Allène X, Gardès L, Rakotoarivony I, Crochet D, Delécolle J, Moulia C, Baldet T: Assessment of vector/host contact: comparison of animal-baited traps and UV-light/suction trap for collecting Culicoides biting midges (Diptera: Ceratopogonidae), vectors of Orbiviruses. Parasites Vectors. 2011, 4: 1-12. 10.1186/1756-3305-4-1.View ArticleGoogle Scholar
- Venter G, Hermanides K, Boikanyo S, Majatladi D, Morey L: The effect of light trap height on the numbers of Culicoides midges collected under field conditions in South Africa. Vet Parasitol. 2009, 166 (3): 343-345. 10.1016/j.vetpar.2009.09.003.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.