Geo-statistical analysis of Culicoides spp. distribution and abundance in Sicily, Italy
- Valeria Blanda†1,
- Marcellocalogero Blanda†1,
- Francesco La Russa1,
- Rossella Scimeca1,
- Salvatore Scimeca1,
- Rosalia D’Agostino1,
- Michelangelo Auteri1 and
- Alessandra Torina1Email authorView ORCID ID profile
© The Author(s). 2018
Received: 25 May 2017
Accepted: 18 January 2018
Published: 1 February 2018
Biting midges belonging to Culicoides imicola, Culicoides obsoletus complex and Culicoides pulicaris complex (Diptera: Ceratopogonidae) are increasingly implicated as vectors of bluetongue virus in Palaearctic regions. Culicoides obsoletus complex includes C. obsoletus (sensu stricto), C. scoticus, C. dewulfi and C. chiopterus. Culicoides pulicaris and C. lupicaris belong to the Culicoides pulicaris complex. The aim of this study was a geo-statistical analysis of the abundance and spatial distribution of Culicoides spp. involved in bluetongue virus transmission. As part of the national bluetongue surveillance plan 7081 catches were collected in 897 Sicilian farms from 2000 to 2013.
Onderstepoort-type blacklight traps were used for sample collection and each catch was analysed for the presence of Culicoides spp. and for the presence and abundance of Culicoides vector species (C. imicola, C. pulicaris / C. obsoletus complexes). A geo-statistical analysis was carried out monthly via the interpolation of measured values based on the Inverse Distance Weighted method, using a GIS tool. Raster maps were reclassified into seven classes according to the presence and abundance of Culicoides, in order to obtain suitable maps for Map Algebra operations.
Sicilian provinces showing a very high abundance of Culicoides vector species were Messina (80% of the whole area), Palermo (20%) and Catania (12%). A total of 5654 farms fell within the very high risk area for bluetongue (21% of the 26,676 farms active in Sicily); of these, 3483 farms were in Messina, 1567 in Palermo and 604 in Catania. Culicoides imicola was prevalent in Palermo, C. pulicaris in Messina and C. obsoletus complex was very abundant over the whole island with the highest abundance value in Messina.
Our study reports the results of a geo-statistical analysis concerning the abundance and spatial distribution of Culicoides spp. in Sicily throughout the fourteen year study. It provides useful decision support in the field of epidemiology, allowing the identification of areas to be monitored as bases for improved surveillance plans. Moreover, this knowledge can become a tool for the evaluation of virus transmission risks, especially if related to vector competence.
Culicoides biting midges (Diptera: Ceratopogonidae) are small hematophagous insects widely spread through different geographical areas, from America to Europe, Asia and Oceania . Breeding sites for Culicoides include a wide range of habitats often located nearby their hosts, in and around farm holdings, such as decaying vegetation, dung, pond borders and moist soils. Furthermore, their larvae can survive within aquatic sites characterized by different ranges of acidity and salinity . Blood-feeding is required for egg production by females, which bite hosts such as amphibians, birds and mammals including humans and domestic animals .
Different species included in the genus Culicoides are of veterinary importance as vectors of different arboviruses causing severe animal diseases. Bluetongue virus (BTV) and African horse sickness virus (AHSV), listed as causative agents of globally important diseases by the OIE (Office International des Epizooties), are transmitted by Culicoides species . Moreover, the species of the Culicoides obsoletus complex have been recently considered the potential vectors of the Schmallenberg virus (SVB), a virus that was not previously detected in Europe . Culicoides spp. can also transmit filarial diseases such as onchocercosis and mansonellosis, affecting various species including humans [3, 5].
Concerning BTV infection, just about 30 Culicoides species have fulfilled the criteria required to be considered as BTV vectors, due to a series of barrier systems preventing virus replication within the biting midges. Recognized BTV vectors include C. imicola (the main vector in Africa, the Middle East, Southeast Asia and areas of southern Europe), C. sonorensis (the main vector in North America) and C. brevitarsis (the main vector in Australia) . Wind-borne transportation of Culicoides or accidental importation of infected hosts can contribute to BTV entry in novel areas . Spreading of species such as Culicoides imicola was possibly favoured by ongoing climate changes, allowing its diffusion in northern Europe as well as an increase of the seasonal activity period, vector density and virus infections susceptibility .
In Europe, BTV vectors belong to three main species or species complexes . Species within the Culicoides obsoletus complex, i.e. Culicoides obsoletus (sensu stricto), Culicoides scoticus, Culicoides dewulfi and Culicoides chiopterus, are putative BTV vectors in northern and central Europe [10, 11]. Other relevant BTV vector species belong to the Culicoides pulicaris complex . Finally, Culicoides imicola is the most important vector in the Mediterranean basin and it may have a role in northern Europe . Within the Mediterrean, Sicily (Italy) has geographical and climatic features particularly suitable for introduction and spread of Culicoides vectors. Sicily has a typical Mediterranean climate, with mild winters, warm autumns/springs and hot summers; however, temperatures can vary among areas depending on the distance from the sea and the presence of mountains.
Studies on the abundance and distribution of Culicoides species involved in vector-borne diseases are essential to define infection risk areas, as well as to identify possible local factors favouring the diffusion of Culicoides and associated pathogens. For such purposes, geographical information system (GIS) technologies have become useful tools for disease mapping, ecological analyses, prediction of parasite occurrence/seasonality and surveillance of parasitic diseases . Thus, GIS technology is increasingly used to carry out a systematic analysis of spatial distribution of vectors and related diseases, usable for the development of effective countermeasures against arthropod-borne diseases.
Our study was focused on a geo-statistical analysis of the abundance and spatial distribution of Culicoides involved in the bluetongue virus transmission, elaborating data obtained from a long monitoring period (2000–2013) in Sicily. Our results will be of importance for improving current surveillance plans in Sicily and may represent a basis for a systematic use of GIS in Culicoides-trasmitted disease risk analysis.
The insects, attracted by the UV light, were sucked into plastic containers containing 200–300 ml of water in which a few drops of detergent were added to reduce surface tension. Collected insects were then placed in vials filled with 70% ethanol.
Culicoides spp. were first separated from other insects under a stereomicroscope. Afterwards, Culicoides species/complexes were identified by the characteristic wing patterns and spots [16–18]. In particular, regarding the species most involved in pathogen transmission, C. imicola specimens were recorded from 2000 to 2013, while the occurrence of species belonging to the C. obsoletus and C. pulicaris complexes was analysed from 2005 to 2013.
Results were recorded in a database and the monthly average abundance was calculated for each farm. The monthly average was chosen as a reliable parameter for estimating Culicoides abundance, since it may vary in the Sicilian territories throughout a month, in a negative or positive way due to transitory adverse or favourable weather conditions. For every month, a geo-statistical analysis via a GIS tool was performed by interpolating measured values with the IDW (inverse distance weighted) method, one of the most important tools in GIS-based spatial interpolation. The IDW method was used for data analysis considering its easy implementation and its suitability to the characteristics of our data, compared to other interpolation methods. This methodology allowed the estimation of values in unsampled points focusing on the basic assumption that values in unsampled points are the weighted average of values in nearby sampled points. Weights are elaborated as the inverse squares of distances. The final objective of IDW analysis is the creation of a smooth surface where values in unsampled points are more similar to nearby points than distant points .
In detail, 12 monthly maps were processed for each species. Once processed, each raster map was reclassified into 7 classes in relation to Culicoides presence and abundance (0, absent to 6, very high), giving a score for each class to each pixel of the image. The abundance classes were the same for each month and for each Culicoides species, and they were defined using the same threshold values.
Through the sum of the monthly maps, it was possible to draw up a map of abundance for each of the three Culicoides species/complexes (C. imicola, C. pulicaris and C. obsoletus). Such maps were further processed to produce a general map comprehensive of all the three species, from which a final map reporting abundance classes in all the Sicilian provinces was elaborated with the evidence of municipalities at risk.
Abundance of Culicoides spp., C. imicola, C. obsoletus complex and C. pulicaris complex based on 7081 catches during 2000–2013
No. of specimens
C. obsoletus complex
C. pulicaris complex
Percentages of areas with different Culicoides abundance in Sicilian provinces
Our work reports the results of a geo-statistical approach applied to the study of the abundance of some relevant Culicoides vector species, analysing a significant number of catches distributed all over the seasons during 14 years. Our approach allows integration of the field data with the ones resulting from the GIS-based IDW method, permitting the elaboration of comprehensive maps including areas where no measurements were available.
As reported, Culicoides spp. are spread all over Sicilian territories, with C. imicola particularly abundant in Palermo, and C. pulicaris complex and C. obsoletus complex being particularly abundant in Messina. Considering all the three species together (Fig. 7), Culicoides are particularly abundant in Messina province and less abundant in Ragusa territories.
Several factors may influence the spread of Culicoides spp., including both biotic (e.g. vegetation, human presence, other animal presence) and abiotic (light, soil, water, air, climatic factors etc.) factors as well as other factors that can vary in terms of macrohabitat or microhabitat. All of these factors are implicated in producing substantial differences in the presence and abundance of Culicoides spp., and local variations occur also at short distance. As an example, C. imicola, is widely distributed across the world, from South Africa to the Mediterranean basin and the Middle and Far East . There are, however, areas with either low abundance or even completely free of C. imicola , underlining a patchy distribution of this species. Other studies have confirmed that the spatial distribution of Culicoides spp. could be influenced both by different climatic zones  as well as by local factors, such as flock management systems. Influence of seasonal and metereological parameters on Culicoides activity was also investigated [22–25] and other studies reported risk maps obtained via the analysis of habitat characteristics [26–28]. As an example, an association between the bluetongue disease diffusion and some landscape metrics, such as woodland and open prairies, was revealed .
Accordingly, even in our study we observed an irregular distribution of Culicoides abundance in Sicily that could also be explained by a possible correlation to specific biotic or abiotic factors in the different provinces. For instance, the high number of farms falling in high risk territories in Messina province could be explained by the farm type, mostly represented by family-owned farms with a small number of animals. Conversely, in Catania province industrial farms with a great number of animals are common. Other factors could influence the distribution pattern of Culicoides in Sicily, including different vegetation, animal presence and microclimatic factors such as altitude, temperature and precipitation. In this study, we provided a picture of Culicoides presence in the Sicilian territories; we related the risk for Culicoides-transmitted diseases with the Culicoides abundance in the different provinces. However, several facets of the data can be further examined and developed, as indeed the correlation with environmental factors and seasonality, which would be our future perspectives of analysis. At present, we have provided an overview about the annual overall abundance of each species, reporting maps obtained by the sum of monthly maps. However, in the perspective of seasonality studies, the specific abundance values within each season/month could be used for correlation with seasonal/monthly environmental parameters. In this view, the implementation of modelling analysis would provide a useful tool to enhance data analysis and to uncover the multiform facets connecting entomological field data to environmental factors, population dynamics, and pathogen diffusion.
The significance of the maps we obtained should also be evaluated in relation to the vector competence of a species towards a pathogen. In Sicily, for example, the bluetongue outbreaks that occurred from 2000 to 2013 severely affected Trapani province, as well as the western part of Palermo province. These occurrences can be explained analysing the map of C. imicola that is a very competent vector of bluetongue virus (Fig. 3). Similarly, to evaluate the risk of Schmallenberg virus spreading, the map reporting C. obsoletus complex distribution has to be considered.
However, the transmission of the same pathogen may be due to several Culicoides species with different vector competence. Thus, correction factors relative to the competence of a Culicoides species towards the examined pathogen should be considered for map elaboration. However, studies on vector competence are complex, requiring highly specialized laboratories, pathogen availability and laboratory colonies of the vector. Nevertheless, for some vectors and pathogens, this information is available. For example, a study conducted on vector competence of C. sonorensis  showed that this species has a competence of 32% towards epizootic hemorragic disease virus and this information certainly allows a precise assessment of the disease risk correlated to this species. Studies providing novel data on competence for the other Culicoides species involved in pathogen transmission could lead to the elaboration of improved maps that better define the risk for Culicoides-related diseases.
In conclusion, this study uses a GIS-based geo-statistical approach to perform a comprehensive analysis of abundance and spatial distribution of Culicoides spp. most involved in virus transmission. We reported data derived from a long time surveillance period (2000–2013). Our study confirmed the high abundance of Culicoides species in Sicily, underlining the presence of C. imicola in Palermo province and of C. pulicaris complex and C. obsoletus complex in Messina province. The analysis provides a useful tool for decision support in the field of epidemiology, allowing the identification of areas to focus on for surveillance purposes. Indeed, the data showed the presence of medium- to high-risk areas in the whole Sicily and of a particularly high-risk area in Messina province. Moreover, if related to vector competence, these data can become an instrument for the prevention of pathogen transmission risks.
We would like to thank Pippo Bono, Gaspare Lo Bue and Nicola Galati (Istituto Zooprofilattico Sperimentale della Sicilia, Palermo, Sicily, Italy) for technical assistance and their contribution in the realization of the catches.
The project was funded by Italian Ministry of Health (Project RC IZS SI 01/13).
Availability of data and materials
All data generated or analysed during this study are included in this published article. The raw datasets generated and/or analysed during the current study are available from the corresponding author upon reasonable request.
VB created and filled the data base for data analysis and wrote the paper. MB elaborated the maps and performed the spatial analysis by IDW (Inverse Distance Weighted) method. MA carried out data analysis. FLR, RS, SS and RDA carried out morphological Culicoides spp. identification. AT organized monitoring activities and project plan. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Borkent A. World species of biting midges (Diptera: Ceratopogonidae), 2012. http://wwx.inhs.illinois.edu/files/8413/4219/9566/CeratopogonidaeCatalog.pdf. Accessed 17 May 2017.
- González M, López S, Mullens BA, Baldet T, Goldarazena A. A survey of Culicoides developmental sites on a farm in northern Spain, with a brief review of immature habitats of European species. Vet Parasitol. 2013;191:81–93.View ArticlePubMedGoogle Scholar
- Mellor PS, Boorman J, Baylis M. Culicoides biting midges: their role as arbovirus vectors. Annu Rev Entomol. 2000;45:307–40.Google Scholar
- Rasmussen LD, Kristensen B, Kirkeby C, Rasmussen TB, Belsham GJ, Bødker R, Bøtner A. Culicoids as vectors of Schmallenberg virus. Emerg Infect Dis. 2012;18:1204–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Linley JR, Hoch AL, Pinheiro FP, et al. J Med Entomol. 1983;20:347–64.View ArticlePubMedGoogle Scholar
- Wilson AJ, Mellor PS. Bluetongue in Europe: past, present and future. Phil Trans R Soc B. 2009;364:2669–81.View ArticlePubMedPubMed CentralGoogle Scholar
- Barros SC, Ramos F, Luis TM, Vaz A, Duarte M, Henriques M, et al. Molecular epidemiology of bluetongue virus in Portugal during 2004–2006 outbreak. Vet Microbiol. 2007;124:25–34.View ArticlePubMedGoogle Scholar
- Purse BV, Mellor PS, Rogers DJ, Samuel AR, PPC M, Baylis M. Climate change and the recent emergence of bluetongue in Europe. Nat Rev Microbiol. 2005;3:171–81.View ArticlePubMedGoogle Scholar
- Wilson AJ, Mellor PS. Bluetongue in Europe: vectors, epidemiology and climate change. Parasitol Res. 2008;103(Suppl. 1):S69–77.View ArticlePubMedGoogle Scholar
- Meiswinkel R, Baldet T, De Deken R, Takken W, Delecolle JC, Mellor PS. The 2006 outbreak of bluetongue in northern Europe - the entomological perspective. Prev Vet Med. 2008;87:55–63.View ArticlePubMedGoogle Scholar
- De Liberato C, Scavia G, Lorenzetti R, Scaramozzino P, Amaddeo D, Cardeti G, et al. Identification of Culicoides obsoletus (Diptera: Ceratopogonidae) as a vector of bluetongue virus in central Italy. Vet Rec. 2005;156:301–4.Google Scholar
- Caracappa S, Torina A, Guercio A, Vitale F, Calabro A, Purpari G, et al. Identification of a novel bluetongue virus vector species of Culicoides in Sicily. Vet Rec. 2003;153:71–4.Google Scholar
- Cringoli G, Rinaldi L, Veneziano V, Musella V. Disease mapping and risk assessment in veterinary parasitology: some case studies. Parassitologia. 2005;47:9–25.PubMedGoogle Scholar
- Venter GJ, Meiswinkel R. The virtual absence of Culicoides imicola (Diptera: Ceratopogonidae) in a light-trap survey of the colder, high-lying area of the eastern Orange free state, South Africa, and implications for the transmission of arboviruses. Onderstepoort J Vet Res. 1994;61:327–40.PubMedGoogle Scholar
- Venter GJ, Labushagne K, Hermanides KG, SNB B, Majatladi DM. Comparison of the efficiency of five suction light traps underfield condition in South Africa for the collection of Culicodes species. Vet Parasitol. 2009;166:299–307.View ArticlePubMedGoogle Scholar
- Goffredo M, Meiswinkel R. Entormological Surveillance of bluetongue in Italy: method of capture, catch analysis and identification of Culicoides biting midges. Vet Ital. 2004;40:260–5.PubMedGoogle Scholar
- Delécolle JC. Nouvelle contribution a l’etude systematique et iconographique des especes du genre Culicoides (Diptera: Ceratopogonidae) du Nord-Est de la France. Strasbourg: Universite Louis Pasteur de Strasbourg, “Vie et Terre”; 1985.Google Scholar
- Rawlings P. A key, based on wing patterns of biting midges (genus: Culicoides Latreille - Diptera: Ceratopogonidae) in the Iberian peninsula, for use in epidemiological studies. Graellsia. 1996;52:57–71.View ArticleGoogle Scholar
- Longley PA, Goodchild MF, Maguire DJ, Rhind DW. Geographic information system and science. Chichester: John Wiley & Sons, Ltd; 2005. p. 333–6.Google Scholar
- Meiswinkel R, Nevill EM, Venter GJ. Vectors: Culicoides spp. In: Coetzer JAW, Thomson GR, Tustin RC, editors. Infectious diseases of livestock with special reference to southern Africa, vol. 1. Cape Town: Oxford University Press; 1994. p. 68–89.Google Scholar
- Brugger K, Rubel F. Characterizing the species composition of European Culicoides vectors by means of the Köppen-Geiger climate classification. Parasit Vectors. 2013;6:333.View ArticlePubMedPubMed CentralGoogle Scholar
- Brugger K, Rubel F. Bluetongue disease risk assessment based on observed and projected Culicoides obsoletus spp. vector densities. PLoS One. 2013;8:e60330.View ArticlePubMedPubMed CentralGoogle Scholar
- Sanders CJ, Shortall CR, Gubbins S, Burgin L, Gloster J, Harrington R, et al. Influence of season and meteorological parameters on flight activity of Culicoides biting midges. J Appl Ecol. 2011;48:1355–136.Google Scholar
- Racloz V, Venter G, Griot C, KDC S. Estimating the temporal and spatial risk of bluetongue related to the incursion of infected vectors into Switzerland. BMC Vet Res. 2008;4:42.View ArticlePubMedPubMed CentralGoogle Scholar
- Ander M, Meiswinkel R, Chirico J. Seasonal dynamics of biting midges (Diptera: Ceratopogonidae: Culicoides), the potential vectors of bluetongue virus, in Sweden. Vet Parasitol. 2012;184:59–67.View ArticlePubMedGoogle Scholar
- Conte A, Ippoliti C, Calistri P, Pelini S, Savini L, Salini R, et al. Towards the identification of potential infectious sites for bluetongue in Italy: a spatial analysis approach based on the distribution of Culicoides imicola. Vet Ital. 2004;40:311–5.Google Scholar
- Takken W, Verhulst N, Scholte EJ, Jacobs F, Jongema Y, Conte A, et al. The phenology and population dynamics of Culicoides spp. in different ecosystems in the Netherlands. Prev Vet Med. 2008;87:41–54.View ArticlePubMedGoogle Scholar
- Caligiuri V, Giuliano GA, Vitale V, Chiavacci L, Travaglio S, Manelli L, et al. Bluetongue surveillance in the Campania region of Italy using a geographic information system to create risk maps. Vet Ital. 2004;40:385–9.PubMedGoogle Scholar
- Guis H, Tran A, de La Roque S, Baldet T, Gerbier G, Barraguè B, et al. Use of high spatial resolution satellite imagery to characterize landscapes at risk of bluetongue. Vet Res. 2007;38:669–83.View ArticlePubMedGoogle Scholar
- Ruder MG, Howerth E, Stallknecht DE, Allison AB, Carter DL, Drolet BS, et al. Vector competence of Culicoides sonorensis (Diptera: Ceratopogonidae) to epizootic hemorrhagic disease virus serotype 7. Parasit Vectors. 2012;5:236.Google Scholar