Culicoides species community composition and feeding preferences in two aquatic ecosystems in northern Spain
Parasites & Vectors volume 15, Article number: 199 (2022)
Aquatic ecosystems provide breeding sites for blood-sucking insects such as Culicoides biting midges (Diptera: Ceratopogonidae), but factors affecting their distribution and host choice are poorly understood. A study was undertaken at two nature reserves in northern Spain to examine the abundance, species composition, population dynamics and feeding patterns of biting midges between 2018 and 2019.
Culicoides were captured by light suction traps baited with CO2 and by sweep netting vegetation. Blood meals and species identification of blood-fed specimens were determined using cytochrome c oxidase I subunit (COI) DNA barcoding. Multivariate generalized linear models were used to evaluate the associations between the abundance of Culicoides, the species richness and other parameters.
The 4973 identified specimens comprised 28 species of Culicoides. These included two species reported for the first time in northern Spain, thus raising to 54 the number of Culicoides species described in the region. Specimens of all 28 species and 99.6% of the total specimens collected were caught in suction traps, while sweep netting vegetation revealed just 11 species and 0.4% of the total specimens. Midge abundance peaked in June/early July, with five species comprising > 80% of the captures: Culicoides alazanicus (24.9%), Culicoides griseidorsum (20.3%), Culicoides poperinghensis (16.2%), Culicoides kibunensis (10.7%) and Culicoides clastrieri (9.6%). DNA barcode analysis of blood meals from eight Culicoides species revealed that they fed on 17 vertebrate species (3 mammals and 14 birds). Species in the subgenus Avaritia were primarily ornithophilic, except for C. griseidorsum and C. poperinghensis. Host DNA from blood meals was successfully amplified from 75% of blood-fed females. A pictorial blood meal digestion scale is provided to accurately assess the blood-fed status of female Culicoides.
The large number of different blood meal sources identified in the midges captured in this study signals the likely importance of wild birds and mammals (e.g. red deer and wild boar) as reservoir/amplifying hosts for pathogens. Available hosts are more exposed to being bitten by biting midge populations in aquatic ecosystems in late spring and early summer.
Tiny blood-sucking midges (< 4 mm long) of the genus Culicoides (Diptera: Ceratopogonidae) are both biting pests and vectors of several viruses, filarial nematodes and protozoa of veterinary health relevance worldwide [1, 2]. At least 110 species of biting midges are included in the latest key for the Western Palearctic region , of which 84 are known from Spain [4,5,6]. In Europe, biting midges are not a threat to human health , but they do play an important role in transmission of both bluetongue virus (BTV) and Schmallenberg virus (SBV) to wild and domesticated animals [2, 7]. Ornithophilic Culicoides species also transmit avian parasites, such as the avian malaria parasite Plasmodium and the closely related genera Haemoproteus and Leucocytozoon [8, 9].
A thorough understanding of host selection (i.e. feeding preferences) by Culicoides females is necessary to determine the complex relationship between hosts, vectors and pathogens . The feeding preferences of midges in the subgenus Avaritia (Culicoides obsoletus group) and subgenus Culicoides (Culicoides pulicaris group) that frequently feed on livestock are much better studied than species with other feeding preferences [11,12,13,14,15]. In general, females of Culicoides exhibit a plastic feeding behavior determined by host preference, with most species feeding primarily on mammals or birds. Regardless of the preferred host, phylogenetically related species tend to feed on the same class(es) of vertebrates .
The role of Culicoides in the transmission of certain haemosporidian parasites remains poorly known as most surveillance and ecological studies have been conducted in agricultural settings with a focus on viral transmission to cattle or other livestock (goats, sheep and horses). As a result, the community of livestock-associated Culicoides species is relatively known in most of Europe [17, 18]. However, biting midges are also prominent in diverse semi-aquatic habitats, including swamps, marshes, forests, ponds, wet pastures, among others [19,20,21,22]. Wetlands and marshes are particularly valuable ecosystems for many animal and plant species, and serve as refuges and destination sites for native and migratory birds as well as for domestic and wild mammals. Consequently, there is high interest in preserving or expanding protected aquatic ecosystems, which also support many species of biting Diptera , including biting midges .
The propagation and transmission of pathogens via biting insects in wetlands and marshes have received little attention. Therefore, this study addresses this gap by examining Culicoides biting midges in these environments. We report data on the relative abundance and population dynamics of Culicoides and use DNA barcoding to identify the vertebrate species that served as the source of blood meals at various post-ingestion stages. The relevance of these findings is then discussed in relation to the monitoring and transmission of pathogens.
This survey was undertaken in the Basque region of northern Spain, at two nature reserves that are important to native and migratory birds (Fig. 1a). One of these areas is the Salburua wetland (Alava province; 42°51′ N, 02°39′ W) and the other is the Urdaibai marsh (Biscay province; 43°22′ N, 02°40′ E). Because these two sites are important wintering and migratory stopovers for multiple species of birds, they are registered in international programs, including the RAMSAR Convention (designating wetlands of international importance, especially as those of waterfowl habitats) and Natura 2000 (a network of nature protection areas within the European Union). These areas are also important public assets and include the Nature Interpretation Center (Salburua) and the Urdaibai Bird Center (Urdaibai), which together host 40,000–50,000 visitors annually. While the bird, mammal and herpetofauna at both sites have been well studied, prior work [25,26,27] on insects has only targeted the Coleoptera, Lepidoptera and Odonata orders; consequently, nothing is known about blood-sucking dipterans from these sites.
Salburua is a fenced park with a freshwater wetland located on the eastern outskirts of the city of Vitoria-Gasteiz. Spanning 2 km2, it is composed of several large pools, grasslands and oak groves bordered partly by an urbanized area with about 7000 residents. Approximately 290 vertebrate species (including 233 native and 4 exotic bird species) have been recorded at Salburua and a herd of about 120 non-native red deer was introduced for vegetation control (Luis Lobo, Centro de Estudios Ambientales, Vitoria-Gasteiz, personal communication). The climate is transitional between Oceanic and Mediterranean, rainfall is moderate (800 mm) and the average temperature is 11 °C, with more than 30 days of frost per year. The pools are fully recharged during the winter, then gradually become almost dry over the summer.
Urdaibai, a Biosphere Reserve designated by UNESCO in 1984, is located on the Bay of Biscay coast. It covers 220 km2 of small streams that merge into a large salt marsh surrounded by meadows, oak groves, woodlands and conifer plantations. Within the territory, there are 20 municipalities inhabited by about 45,000 people and some small cattle farms. Reflecting its large size and multi-ecosystem landscape, Urdaibai hosts at least 318 vertebrate species (approx. 250 birds), making it the most valuable ecosystem on the northern Spanish coast . The reserve is constantly flooded, but the extent of flooding oscillates during the summer. The marsh ecosystem varies between freshwater and saline depending on location. The climate is Oceanic with abundant annual rainfall (1200 mm) and mild temperatures (average temperature of 14 °C) with less than 1 day of frost per year.
Five distinct environments were sampled from 1 July to 31 October 2018 (two habitats) and from 1 May to 31 October (three habitats) in each reserve. The five sampling sites in Salburua included: aquatic vegetation alongside a large temporary wetland pool (2018); a humid oak grove with streams (2018); grassland near a slow-flowing stream (2019); the margin of small shallow temporary ponds (2019); and a grassy poplar grove (2019). Sites in Urdaibai included: a humid woodland with intermittent ephemeral puddles (2018); vegetation alongside a permanent freshwater marsh (2018); the contact zone between a permanent freshwater habitat vegetated with bulrush; a mixed forest covered by ferns and brambles (2019); a temporary saline marsh with aquatic vegetation, primarily Tamarix spp. (2019); and mixed patches of trees and bushes admixed with streams, puddles and ditches near a livestock farm (2019).
Carbon dioxide-baited CDC-traps
Centers for Disease Control and Prevention (CDC) 6 V battery-powered miniature traps (model 1212; John W Hock Co., Gainesville, FL, USA) equipped with incandescent light and baited with about 1.5 kg of dry ice that is released for 24 h through drilled polyethylene boxes, were employed to collect adult-biting midges as well as other blood-sucking Diptera (Fig. 1b). The traps were suspended at a height of 1–1.5 m and operated for 24 h periods (set up early in the morning and retrieved the next morning) every 2 weeks. Traps were hung on tree branches, where they were protected from sunlight and wind exposure. Collection pots were immediately transported to the laboratory and insects were stored at − 30 °C. A total of 55 carbon dioxide (CO2)-baited CDC traps were examined (16 in 2018 and 39 in 2019) in each wetland.
Biting midges resting on vegetation were collected in the morning (between 9:30 and 11:00 a.m.) using a polyester mesh long-handled net (diameter: 38 cm, mesh size: 0.8 mm; BioQuip Products Inc., Rancho Dominguez, CA, USA) (Fig. 1b) by sweeping through grassy and shrubby vegetation (and on the foliage of bushes and shrubs) in a 25-m radius around each CDC trap. Sweeping was conducted very 2 weeks for approximately 4 min at each location. Sweep netting locations were limited by the type of habitat and time of year because flooding in the spring and densely vegetated areas in the summer restricted the effectiveness of this collection method. Collections were made by the same person each time, and collections were kept at – 30 °C for 20 min to kill the insects. A total of 55 sweep netting collections were examined (16 in 2018 and 39 in 2019) in each wetland.
Once in the laboratory, insects were sorted into major taxonomic groups. Culicoides midges were preserved in ethanol (70%), while other blood-feeding arthropods (data not shown) were preserved frozen for other studies. All biting midges were separated by sex and identified to species or species-group under a stereo microscope based on the wing pattern pigmentation, the shape and size of the third palpal segment (females) and other diagnostic morphological traits . Males and females that presented atypical phenotypes, were damaged or were sibling species, or those whose identity was otherwise uncertain were mounted (after dissection of head, thorax + legs, wings and abdomen) on slides with permanent media (Hoyer´s medium), and species identification was achieved under a compound microscope. All species identifications were validated with the interactive identification key for Western Palaearctic Culicoides . For species within the Obsoletus complex (Culicoides obsoletus and Culicoides scoticus), females were pooled due to the difficulties in discriminating these species morphologically, while males were identified based on diagnostic characters of the male genitalia [30, 31]. Blood-fed (n = 68) and gravid females (n = 12) were also mounted on slides, excluding the abdomen and thorax, which was used for molecular analysis.
Molecular identification of blood-fed Culicoides species and their hosts
To maximize the sample size and assess the limit of detection of the molecular technique, all Culicoides females (n = 80) showing any trace of blood in their abdomen were analyzed. Specimens with successful amplification of host DNA (n = 53) were also DNA barcoded to confirm their species assignment. The specimens were also photographed to illustrate the five stages in blood meal digestion. To achieve these aims, the abdomen and thorax of each Culicoides midge were individually transferred to sterile vials (2 ml) and shipped on dry ice to the Centre for Biodiversity Genomics, University of Guelph (Guelph, Canada) for molecular analysis. Samples were processed following previously established methods [32, 33]. Briefly, DNA was extracted using a modified glass fiber technique . The resulting DNA was used to ascertain the identities of the Culicoides species as well as of the vertebrate hosts upon which they had fed. Culicoides species were identified using standard DNA barcoding techniques, employing universal insect primers (C_LepFolF + C_LepFolR)  followed by Sanger sequencing. Traces were edited in CodonCode Aligner v9.0.1 and uploaded to the Barcode of Life Data System (BOLD). For vertebrate host identification, primers were designed to anneal to vertebrate but not insect DNA (C_BloodmealF1_t1 + Mod.Mamm.R_t1)  followed by next-generation sequencing on an Ion Torrent S5 Sequencer (Thermo Fisher Scientific, Waltham, MA). The resulting sequence reads were processed by first removing reads with a quality score (QV) < 20. Following primer/adapter trimming, reads in the expected size range of 125–250 bp were clustered into operational taxonomic units (OTUs) with a minimum identity of 98%. OTUs represented by at least 10 reads were compared to a reference library consisting of all vertebrate cytochrome c oxidase I gene (COI) barcode records on BOLD. Matches between an OTU sequence and a reference sequence were considered reliable only if at least 100 bp of the query sequence matched a reference with at least 95% homology. For each biting midge, all taxonomic matches were consolidated into a “unique taxonomic hit” table, with each hit supported by a total read count. Any taxonomic hits that occurred in negative control samples were proportionally subtracted from all other read counts, after which hits were only accepted as genuine if supported by at least 100 reads.
Full details for each Culicoides specimen, as well as their sequence information, can be found at the Barcode of Life Database (BOLD) within the “Human Pathogens and Zoonoses Initiative” Working Group 1.4. The Digital Object Identifier (DOI) for the publically available dataset on BOLD is https://doi.org/10.5883/DS-CULSPAIN. Accession numbers for all sequences were obtained from NCBI (OL702716-OL702759). Sequences for the Culicoides was analyzed in MEGA v.6  and a neighbor-joining (NJ) analysis was performed using the Kimura 2-parameter distance. A Barcode Index Number (BIN) was assigned to all sequences longer than 500 bp, and each BIN was mapped according to species.
Data on population dynamics were only plotted for 2019 because data were available for 6 months versus just 4 months for 2018. Generalized linear model (GLM) analysis was performed to evaluate the associations between Culicoides abundance (catch per trap per night) and species richness (S, number of species captured per trap per night) over the sampling period (July–October, shared trapping period for both years) with regards to aquatic ecosystems (Salburua and Urdaibai) and the year (2018 and 2019). Due to overdispersion of the data, a negative binomial GLM (NBGLM) was applied  using the MASS package . For species richness, a GLM with Poisson error distribution and log-link function was used. The best model was selected with the “MuMIn” package of R software using the “dredge” function , which is based on the Akaike information criterion and corrected for sample size (AICc). The overall fit of the model was evaluated with a likelihood ratio test, comparing the best model with the null model. Statistical analyses were performed using R statistical software version 3.6.1 . Relative abundance (RA) was calculated as the number of biting midges captured by one sampling method compared to the total number biting midges captured by both methods.
A total of 4973 Culicoides specimens were collected at the 10 sampling sites in 2018 and 2019 using suction traps baited with CO2 (Table 1) and sweep netting (Table 2). Other blood-sucking arthropods (excluding mosquitoes) were also identified, including eight specimens of Chrysops viduatus (Diptera: Tabanidae) and six specimens of Simulium spp. (Diptera: Simuliidae).
Culicoides species community composition
Among the Culicoides specimens identified in this study, 4748 were female and 225 were male, and collectively they comprised 28 species (Tables 1, 2). Five species of the subgenus Oecacta, all known to feed on birds, dominated the assemblage. They included Culicoides alazanicus (n = 1238, 24.9%), Culicoides griseidorsum (n = 1,012, 20.3%), Culicoides poperinghensis (n = 808, 16.2%), Culicoides kibunensis (n = 531, 10.7%) and Culicoides clastrieri (n = 479, 9.6%). Species in the subgenus Avaritia (C. obsoletus, C. scoticus and C. dewulfi) were scarce (n = 66, 1.3%), while the remaining specimens represented 20 other species (Table 1). The three most common species at Salburua were C. griseidorsum (n = 985, 26.0%), C. alazanicus (n = 976, 25.7%) and C. poperinghensis (n = 793, 20.9%), whereas in Urdaibai they were C. clastrieri (n = 316, 26.8%), C. alazanicus (n = 262, 22.2%) and C. punctatus (n = 165, 14.0%) (Tables 1, 2). The distribution of Culicoides species varied between the two sites as they only shared 17 species (Fig. 2). Morphological traits were found to allow unequivocal species-level identifications for most slide-mounted specimens, excluding sibling species in the Obsoletus complex, which have an indistinct wing pattern. Two species (Culicoides albicans and Culicoides puncticollis) from Urdaibai and Salburua, respectively, represented first records from the Basque region.
Biting midges showed a single abundance peak in early summer (15 June–15 July). Midges from Salburua increased in number throughout the spring (May) until the maximum emergence in early July when 41% of all specimens were collected within an interval of 2 weeks. A similar trend was noted at Urdaibai, but peak abundance was 2 weeks later (i.e. capture of 76% of the total collection was concentrated in mid-late June) (Fig. 3). The most abundant species showed differing seasonal patterns. The earliest emerging species was C. poperinghensis (peaked in May) followed by C. griseidorsum (peaked in June) and C. clastrieri, C. alazanicus and C. kibunensis (peaked in early July) (Fig. 3).
Analysis of variables affecting catch, abundance and species richness of Culicoides
Biting midges were mainly trapped by CO2-baited CDC traps (4945 specimens; relative abundance: 99.4%; species richness: 28) (Table 1) and rarely by sweeping (28 specimens; relative abundance: 0.6%; species richness: 11) (Table 2). Total Culicoides abundance varied between study sites, seasons and years. The abundance of biting midges was significantly higher in Salburua than in the marsh of Urdaibai (Table 3). The number of Culicoides specimens was also associated with the month (Table 3), with the highest catches in July, followed by a sharp decrease in August and a slight rebound in September. An annual variation was also observed, with a higher catch per trap in 2018 than in 2019. Species richness was higher at Salburua than Urdaibai, with more species during July followed by a steady decrease from August to October (Table 3).
Molecular identification of Culicoides and their blood meals
Among the 53 Culicoides females that were barcoded to confirm their identity, 44 yielded a COI DNA barcode sequence. NJ analysis (Additional file 1: Phylogenetic analysis) showed that most barcode sequences clustered as expected based on morphological identifications. However, the COI barcode did not provide good resolution for C. clastrieri versus Culicoides festivipennis and C. griseidorsum versus Culicoides pictipennis, as these species grouped with the same BIN (BOLD: AEB9007 and ACV0334, respectively). Conversely, specimens of C. kibunensis were placed in two BINs that showed a mean divergence of 6.5%.
Among the trapped Culicoides, only 1.6% (n = 68) were scored as blood-fed (containing at least visible remains of blood or red tegument, stages 2–5), while 12 more were scored as gravid. Figure 4 provides a pictorial blood meal digestion scale for Culicoides females. We identified the host DNA source for 75% (51/68) of the females in stages 2 to 5 and for 17% (2/12) of the gravid females. The success of blood meal identification was 100% for fully undigested blood meals (stage 2, n = 9) and early digested blood meals (stage 3, n = 10) but declined to 83% in females in stage 4 (20/24) and to 48% in those in stage 5 (12/25).
In total, 17 vertebrate species (3 mammals and 14 birds) were identified as the hosts for the eight Culicoides species with blood meals (Table 4). Of these specimens, 60% were found to have fed on avian hosts (32/53) and 40% on mammalian hosts (21/53). The most common mammal hosts were Sus scrofa and Cervus elaphus, while Turdus spp. and Sylvia atricapilla were the dominant bird hosts (Table 4). Culicoides midges showed a predominant affinity to feed on birds (subgenus Oecacta) with the exception of C. griseidorsum (Table 4).
Blood-fed Culicoides midges were trapped primarily in June (n = 44) and July (n = 34) and only rarely in other months (May = 1; September = 1). The proportion of mammal versus avian hosts varied between the two sites, with a much higher proportion of mammalian hosts at Salburua.
An in-depth understanding of host-feeding preferences and the species composition of the Culicoides community at a site is a key component of disease surveillance programs and ecosystem health assessments. The investigation reported here revealed a high species richness and abundance of biting midges at two aquatic ecosystems in northern Spain. Interestingly, members of the Obsoletus complex, which are very abundant in most European agricultural settings , were very uncommon in these aquatic ecosystems, where ornithophilic Culicoides species were predominant, as has been reported in other studies . As well, the most abundant species in this study (i.e. C. alazanicus and C. griseidorsum) showed a similar dominance in other aquatic ecosystems [42,43,44]. CDC-traps baited with dry ice captured more Culicoides than sweep netting, but the sampling effort (i.e. time) to collect midges by CDC-traps was much higher (24 h/trap vs. 4 min by sweep netting/sampling site). Although sweep netting is not often used to collect adult Culicoides , it can be effective for diurnal species, swarming aggregations and species that are rarely attracted to light, or to identify their potential resting places [46, 47].
Vector identification is key for the surveillance of arthropod-borne diseases as great differences in transmission capacity have been reported between closely related species . Although a recently updated and illustrated interactive key has helped to identify Culicoides in Europe , the morphological identification of Culicoides species is difficult without training. DNA barcoding has been widely adopted as a tool for the rapid identification of insect species . NJ analysis of the Culicoides barcode sequences generated in this study revealed that most conspecific specimens showed close sequence congruence, while members of different species showed deep divergence. However, specimens of C. clastieri grouped with specimens of C. festivipennis and were placed in the same BIN. Another study also reported low interspecific divergence between these two species , but they can be discriminated by morphometrics (number of spines in the cibarium, distribution of sensillae coeloconia, size of the R5 spot and color of thorax) . Culicoides griseidorsum and C. pictipennis also shared barcode sequences but also could be distinguished by morphological characters. A converse situation was apparent for C. kibunensis, as specimens identified to this species belonged to two sequence clusters showing nearly 7% sequence divergence, suggesting the presence of sibling species.
In Europe, Culicoides midges in farmland settings are well-documented, and a marked variation in species abundances linked to trapping method, latitude and season has been shown . However, much less data are available on Culicoides species in natural habitats. In contrast to farm-associated Culicoides, which show multiple abundance peaks from March to October , our study showed a single pronounced peak in early summer. This difference suggests that the most commonly trapped species (C. alazanicus, C. griseidorsum, C. poperinghensis and C. kibunensis) are univoltine, while species common in agricultural settings (C. festivipennis, Culicoides newsteadi and C. obsoletus) are multivoltine [17, 52]. This difference might be related to breeding sites, which are clearly different between agricultural and natural settings.
The factors responsible for the differences observed between 2018 and 2019 are uncertain, but they could reflect the shift in the sampling sites and/or more favorable climatic conditions in 2018. Threefold more midges were collected in the Salburua wetland, likely due to the abundance of freshwater habitats near the sampling sites [53,54,55], while the Urdaibai marsh was dominated by permanent bodies of saline water. The latter habitats may be less appropriate as developmental sites for the Culicoides species found in Europe, although salt marsh species are common in North America .
The use of blood meal DNA for host identification is constrained by its degradation and by the scarcity of blood-fed specimens. Prior work on mosquitoes has shown that the likelihood of recovering host DNA diminishes as digestion of the blood meal progresses until the formation of eggs [41,42,43]. However, this question is fairly unknown for Culicoides, and Culicoides females have typically been assigned to just two stages, i.e. fully/partly engorged or advanced digestion [11, 57]. To provide greater precision, we classified the digestion status of Culicoides into five categories reflecting the extent of digestion. Host identification was successful in all Culicoides females at stages 2–3, but it also worked, although at a lower efficiency, in specimens showing an advanced degree of blood digestion (stages 4–6) and in gravid Culicoides (stage 7). A South African study recovered host DNA from 19% of parous and 26% of gravid Culicoides specimens , supporting the value of analyzing all abdomens to assess host preferences. This is important because fully engorged females are rarely trapped by light-suction traps. By analyzing blood-fed females at all stages of digestion, we were able to identify the source of the blood meal in 75% of specimens, similar to success rates reported in other studies (44–91%) [15, 42, 59,60,61,62].
To understand and control the spread of pathogens through a community, the identity of both susceptible hosts and insects that vector the pathogens to them must be known. Our blood meal analysis of field-collected Culicoides females revealed differing host preferences, information essential for inferring their vector status . The current study provides new perspectives because it is one of the few undertaken in natural settings as most earlier studies examined the Culicoides feeding preferences on a few domestic livestock host species. As might be expected, Culicoides in natural settings feed on a broad range of avian and mammal hosts. Our results support conclusions reached by Martínez de la Puente et al.  in that we found that members of the subgenus Oecacta (C. alazanicus, C. festivipennis, C. kibunensis, C. clastrieri and Culicoides duddingstoni) only fed on birds, with the exception of C. griseidorsum and C. poperinghensis which were primarily found to attack red deer (C. elaphus), in contrast to results reported in a previous study . This is interesting because Culicoides species that feed on both wild and domestic ruminants act as bridge vectors . Two specimens from the subgenus Avaritia (C. obsoletus) fed on a bird and a mammal. Those species that showed biting preferences for a wide range of birds deserve attention as some of them can transmit avian malaria .
One Culicoides specimen trapped at Urdaibai was found to have fed on cattle from a farm located 150 m from the trapping site and could only have reached the trap by crossing a dense forest patch. This relatively local dispersion after blood ingestion supports previous reports [11, 65]. Although some studies have reported longer dispersal distances , such cases likely reflect wind-mediated passive dispersal . The distinctive differences in blood meals observed between the two study sites might reflect the differing availability of vertebrate hosts at each site. For example, the large population of red deer in Salburua likely provided easy access to hosts for Culicoides. No human-derived blood meals were recorded despite the presence of human settlements near both aquatic settings. This result was not a surprise as local Culicoides species are not attracted to humans, in contrast to the situation in other parts of Europe [60, 64, 67,68,69,70]. Most of the vertebrate hosts identified in this work are common in the study areas but some, such as the golden oriole (Oriolus oriolus), are summer trans-Saharan migrant species. Such instances exemplify how the use of blood meal analysis can be used to monitor rare species within a community [71, 72]. Finally, this survey contributes to knowledge on the species biodiversity of biting midges in the region as it raised the number of Culicoides species from 52  to 54 species. Based on these results, the Basque Country region supports at least 64% of the total Spanish Culicoides fauna.
This study expands our understanding of the Culicoides community at two sites in northern Spain, particularly with regards to host-use. Because this study is one of very few that has examined Culicoides in natural ecosystems as opposed to artificial agricultural settings, the information on host blood meals is more representative of the true host-range preference. The fact that most midges fed primarily on avian hosts has significant implications as it clearly suggests that Culicoides likely play an important role in vectoring avian parasites, particularly over early summer. In addition, some species fed on wild mammals (e.g. red deer and wild boar), which could act to amplify the host–vector cycle of several viruses, which in turn may affect livestock. Further studies should identify pathogens as well as the respective host(s) and vector(s) to help decipher their transmission dynamics.
Availability of data and materials
All data generated and analysed during this study are included in this published article.
Akaike information criterion
Barcode Index Number
Barcode Of Life Database
Centers for Disease Control and Prevention
Cytochrome c oxidase I gene
Digital object identifier
Generalized linear models
Negative binomial generalized linear model
Operational Taxonomic Unit
Carpenter S, Groschup MH, Garros C, Felippe-Bauer ML, Purse BV. Culicoides biting midges, arboviruses and public health in Europe. Antiviral Res. 2013;100:102–13.
Mellor PS, Boorman J, Baylis M. Culicoides biting midges: their role as arbovirus vectors. Annu Rev Entomol. 2000;45:307–40.
Mathieu B, Cetre-Sossah C, Garros C, Chavernac D, Balenghien T, Carpenter S, et al. Development and validation of IIKC: an interactive identification key for Culicoides (Diptera: Ceratopogonidae) females from the Western Palaearctic region. Parasit Vectors. 2012;5:137.
Alarcón-Elbal PM, Lucientes J. Actualización del catálogo de Culicoides Latreille, 1809 (Diptera, Ceratopogonidae) de España. Graellsia. 2012;68:353–62.
Talavera S, Muñoz-Muñoz F, Pagès N. New insights on diversity, morphology and distribution of Culicoides Latreille 1809 (Diptera: Ceratopogonidae) from Northeast Spain. Ann Soc Entomol Fr. 2011;47:214–31.
Sánchez Murillo JM, González M, Martínez Díaz MM, Reyes Galán A, Alarcón-Elbal PM. Primera cita de Culicoides paradoxalis Ramilo & Delécolle, 2013 (Diptera, Ceratopogonidae) en España. Graellsia. 2015;71:e033.
Rasmussen LD, Kristensen B, Kirkeby C, Rasmussen TB, Belsham GJ, Bodker R, et al. Culicoids as vectors of Schmallenberg virus. Emerg Infect Dis. 2012;18:1204–6.
Atkinson CT, van Riper III C. Pathogenicity and epizootiology of avian haematozoa: Plasmodium, Leucocytozoon, and Haemoproteus. In: Loye JE and Zuk M, editors. Bird-Parasite Interactions. Ecology, Evolution, and Behavior. New York: Oxford University Press; 1991. p.19–48.
Valkiunas G. Avian malaria parasites and other Haemosporina. Boca Raton: CRC Press; 2005.
Veiga J, Martínez-de la Puente J, Vaclav R, Figuerola J, Valera F. Culicoides paolae and C. circumscriptus as potential vectors of avian haemosporidians in an arid ecosystem. Parasit Vectors. 2018;11:524.
Bartsch S, Bauer B, Wiemann A, Clausen PH, Steuber S. Feeding patterns of biting midges of the Culicoides obsoletus and Culicoides pulicaris groups on selected farms in Brandenburg. Germany Parasitol Res. 2009;105:373–80.
Garros C, Gardes L, Allene X, Rakotoarivony I, Viennet E, Rossi S, et al. Adaptation of a species-specific multiplex PCR assay for the identification of blood meal source in Culicoides (Ceratopogonidae: Diptera): applications on Palaearctic biting midge species, vectors of Orbiviruses. Infect Genet Evol. 2011;11:1103–10.
Lassen SB, Nielsen SA, Kristensen M. Identity and diversity of blood meal hosts of biting midges (Diptera: Ceratopogonidae: Culicoides Latreille) in Denmark. Parasit Vectors. 2012;5:143.
Martínez-de la Puente J, Martinez J, Ferraguti M, Morales-de la Nuez A, Castro N, Figuerola J. Genetic characterization and molecular identification of the bloodmeal sources of the potential bluetongue vector Culicoides obsoletus in the Canary Islands. Spain Parasit Vectors. 2012;5:147.
Ninio C, Augot D, Delecolle JC, Dufour B, Depaquit J. Contribution to the knowledge of Culicoides (Diptera: Ceratopogonidae) host preferences in France. Parasitol Res. 2011;108:657–63.
Martínez-de la Puente J, Figuerola J, Soriguer R. Fur or feather? Feeding preferences of species of Culicoides biting midges in Europe. Trends Parasitol. 2015;31:16–22.
Cuellar AC, Kjaer LJ, Kirkeby C, Skovgard H, Nielsen SA, Stockmarr A, et al. Spatial and temporal variation in the abundance of Culicoides biting midges (Diptera: Ceratopogonidae) in nine European countries. Parasit Vectors. 2018;11:112.
Mignotte A, Garros C, Gardes L, Balenghien T, Duhayon M, Rakotoarivony I, et al. The tree that hides the forest: cryptic diversity and phylogenetic relationships in the Palaearctic vector Obsoletus/Scoticus Complex (Diptera: Ceratopogonidae) at the European level. Parasit Vectors. 2020;13:265.
Mullen GR, Murphree CS. Biting midges (Ceratopogonidae). In: Mullen GR, Durden LA, editors. Medical and veterinary entomology. San Diego: Academic Press; 2019. p. 213–236.
Purse BV, Carpenter S, Venter GJ, Bellis G, Mullens BA. Bionomics of temperate and tropical Culicoides midges: knowledge gaps and consequences for transmission of Culicoides-borne viruses. Annu Rev Entomol. 2015;60:373–92.
Santiago-Alarcon D, Havelka P, Pineda E, Segelbacher G, Schaefer HM. Urban forests as hubs for novel zoonosis: blood meal analysis, seasonal variation in Culicoides (Diptera: Ceratopogonidae) vectors, and avian haemosporidians. Parasitology. 2013;140:1799–810.
Zimmer JY, Smeets F, Simonon G, Fagot J, Haubruge E, Francis F, et al. Are bogs reservoirs for emerging disease vectors? Evaluation of Culicoides populations in the Hautes Fagnes Nature Reserve (Belgium). PLoS ONE. 2013;8:e66893.
Dale PER, Knight JM. Wetlands and mosquitoes: a review. Wetl Ecol Manag. 2008;16:255–76.
Hendry GAF, Godwin G. Biting midges in Scottish forestry: a costly irritant or a trivial nuisance? Sott For. 1988;42:113–9.
Bahillo de la Puebla P, López-Colón JI, Alonso-Román I. Presencia de Triplax lepida (Faldermann, 1837) (Coleoptera, Erotylidae) en la Reserva de la Biosfera de Urdaibai (Vizcaya, Norte de España). Arquivos Entomolóxicos. 2011;5:31–2.
Vega FJ, García-Criado F, Miguélez D, Valladares LF. Diversidad de odonatos en los humedales rehabilitados del Parque Natural de Salburua (Álava). Est Mus Cienc Nat Álava. 2005;20:107–14.
de Juana F, Monasterio Y, Jiménez R, Albalá J, Belamendia G, Olano I, et al. Los macroheteróceros (Lepidoptera) de los humedales de Salburua (Vitoria-Gasteiz, Araba/Alava, España): un proyecto de ciencia ciudadana. Bol SEA. 2019;64:165–85.
García E, Rebane M. Where to watch birds in Northern & Eastern Spain. London: Helm; 2017.
González M, Goldarazena A. El género Culicoides en el País Vasco: guía práctica para su identificación y control. Vitoria-Gasteiz: Gobierno Vasco-Eusko Jaurlaritza; 2011.
Augot D, Sauvage F, Jouet D, Simphal E, Veuille M, Couloux A, et al. Discrimination of Culicoides obsoletus and Culicoides scoticus, potential bluetongue vectors, by morphometrical and mitochondrial cytochrome oxidase subunit I analysis. Infect Genet Evol. 2010;10:629–37.
Nielsen SA, Kristensen M. Morphological and molecular identification of species of the Obsoletus group (Diptera: Ceratopogonidae) in Scandinavia. Parasitol Res. 2011;109:1133–41.
González MA, Prosser SW, Hernández-Triana LM, Alarcón-Elbal PM, Goiri F, López S, et al. Avian feeding preferences of Culex pipiens and Culiseta spp. along an urban-to-wild gradient in Northern Spain. Front Ecol Evol. 2020;8:352.
Estrada-Franco JG, Fernández-Santos NA, Adebiyi AA, López-López MJ, Aguilar-Duran JA, Hernández-Triana LM, et al. Vertebrate-Aedes aegypti and Culex quinquefasciatus (Diptera)-arbovirus transmission networks: Non-human feeding revealed by meta-barcoding and next-generation sequencing. PLoS Negl Trop Dis. 2020;14:e0008867.
Ivanova NV, Dewaard JR, Herbert PDN. An inexpensive, automation-friendly protocol for recovering high-quality DNA. Mol Ecol Notes. 2006;6:998–1002.
Hernández-Triana LM, Prosser SW, Rodríguez-Pérez MA, Chaverri LG, Hebert PD, Gregory TR. Recovery of DNA barcodes from blackfly museum specimens (Diptera: Simuliidae) using primer sets that target a variety of sequence lengths. Mol Ecol Resour. 2014;14:508–18.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
O’Hara RB, Kotze DJ. Do not log-transform count data. Methods Ecol Evol. 2010;1:118–22.
Venables WN, Ripley BD. Modern applied statistics with S. New York: Springer; 2002.
Barton K. MuMIn: Multi-model inference. R Package version 1.43.17. 2020. https://CRAN.R-project.org/package=MuMIn. Accessed 1 Dec 2021.
R Core Team. A language and environment for statistical computing. 2019. https://www.r-project.org/. Accessed 1 Dec 2021.
Mohlmann TWR, Wennergren U, Talle M, Favia G, Damiani C, Bracchetti L, et al. Community analysis of the abundance and diversity of biting midge species (Diptera: Ceratopogonidae) in three European countries at different latitudes. Parasit Vectors. 2018;11:217.
Tomazatos A, Jost H, Schulze J, Spinu M, Schmidt-Chanasit J, Cadar D, et al. Blood-meal analysis of Culicoides (Diptera: Ceratopogonidae) reveals a broad host range and new species records for Romania. Parasit Vectors. 2020;13:79.
Bobeva A, Zehtindjiev P, Ilieva M, Dimitrov D, Mathis A, Bensch S. Host preferences of ornithophilic biting midges of the genus Culicoides in the Eastern Balkans. Med Vet Entomol. 2015;29:290–6.
Bobeva A, Mathieu B. Preliminary results on species composition of the biting midges-fauna (Culicoides) in a wetland on Lower Danube Flow Bulgaria. ARPHA Conf Abstr. 2019;2:e39697.
McDermott EG, Lysyk TJ. Sampling considerations for adult and immature Culicoides (Diptera: Ceratopogonidae). J Insect Sci. 2020;20:2.
Elbers ARW, Meiswinkel R. Limited attractant range of the black-light suction trap for the capture of Culicoides biting midges (Diptera: Ceratopogonidae). J Appl Entomol. 2016;140:386–94.
González MA, Alarcón-Elbal PM, Venter GJ, López S. Flight and swarming behaviour of Culicoides species (Diptera: Ceratopogonidae) on a livestock farm in Northern Spain. Vet Ital. 2017;53:157–66.
Pagès N, Muñoz-Muñoz F, Talavera S, Sarto V, Lorca C, Núñez JI. Identification of cryptic species of Culicoides (Diptera: Ceratopogonidae) in the subgenus Culicoides and development of species-specific PCR assays based on barcode regions. Vet Parasitol. 2009;165:298–310.
Stur E, Borkent A. When DNA barcoding and morphology mesh: Ceratopogonidae diversity in Finnmark. Norway Zookeys. 2014;463:95–131.
Hadj-Henni L, Djerada Z, Millot C, Augot D. Comprehensive characterisation of Culicoides clastrieri and C. festivipennis (Diptera: Ceratopogonidae) according to morphological and morphometric characters using a multivariate approach and DNA barcode. Sci Rep. 2021;11:521.
González MA, Thierry B, Delecolle JC, López S, Romon P, Goldarazena A. Monitoring of Culicoides Latreille (Diptera: Ceratopogonidae) after BTV outbreaks, in sheep farms and natural habitats from the Basque Country (Northern Spain). Proc Entomol Soc Wash. 2013;115:48–69.
Cuellar AC, Kjaer LJ, Baum A, Stockmarr A, Skovgard H, Nielsen SA, et al. Modelling the monthly abundance of Culicoides biting midges in nine European countries using random forests machine learning. Parasit Vectors. 2020;13:194.
Foxi C, Delrio G. Larval habitats and seasonal abundance of Culicoides biting midges found in association with sheep in northern Sardinia. Italy Med Vet Entomol. 2010;24:199–209.
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.
Uslu U, Dik B. Description of breeding sites of Culicoides species (Diptera: Ceratopogonidae) in Turkey. Parasite. 2007;14:173–7.
Carrasco D, Felippe-Bauer ML, Dumont L, D’Incao F. Abundance of Culicoides (Diptera, Ceratopogonidae) species in salt marshes of the Patos Lagoon estuary, Rio Grande do Sul, Brazil: Influence of climatic variables. Pan-Am J Aquat Sci. 2014;9:8–20.
Hadj-Henni L, De MT, Depaquit J, Noel P, Germain A, Helder R, et al. Comparison of vertebrate cytochrome b and prepronociceptin for blood meal analyses in Culicoides. Front Vet Sci. 2015;2:15.
Riddin MA, Venter GJ, Labuschagne K, Villet MH. Bloodmeal analysis in Culicoides midges collected near horses, donkeys and zebras in the Eastern Cape, South Africa. Med Vet Entomol. 2019;33:467–75.
Calvo JH, Berzal B, Calvete C, Miranda MA, Estrada R, Lucientes J. Host feeding patterns of Culicoides species (Diptera: Ceratopogonidae) within the Picos de Europa National Park in northern Spain. Bull Entomol Res. 2012;102:692–7.
England ME, Pearce-Kelly P, Brugman VA, King S, Gubbins S, Sach F, et al. Culicoides species composition and molecular identification of host blood meals at two zoos in the UK. Parasit Vectors. 2020;13:139.
Talavera S, Muñoz-Muñoz F, Verdún M, Pujol N, Pagès N. Revealing potential bridge vectors for BTV and SBV: a study on Culicoides blood feeding preferences in natural ecosystems in Spain. Med Vet Entomol. 2018;32:35–40.
Videvall E, Bensch S, Ander M, Chirico J, Sigvald R, Ignell R. Molecular identification of bloodmeals and species composition in Culicoides biting midges. Med Vet Entomol. 2013;27:104–12.
Cator LJ, Johnson LR, Mordecai EA, Moustaid FE, Smallwood TRC, LaDeau SL, et al. The role of vector trait variation in vector-borne disease dynamics. Front Ecol Evol. 2020;8:189.
Santiago-Alarcon D, Havelka P, Schaefer HM, Segelbacher G. Bloodmeal analysis reveals avian Plasmodium infections and broad host preferences of Culicoides (Diptera: Ceratopogonidae) vectors. PLoS ONE. 2012;7:e31098.
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.
Kluiters G, Swales H, Baylis M. Local dispersal of palaearctic Culicoides biting midges estimated by mark-release-recapture. Parasit Vectors. 2015;8:86.
Hristescu D, Barbuceanu F, Dascalu L, Nitescu C, Goffredo M, Santilli A, et al. Species composition and relative abundance of the genus Culicoides (Diptera: Ceratopogonidae) in Romania. Parasit Vectors. 2020;13:393.
Kampen H, Werner D. Blood meal analysis in Central European Culicoides biting midge species (Diptera: Ceratopogonidae) by reverse line blot hybridization. Mitt Dtsch Ges Allg Angew Ent. 2020;22:225–32.
Romiti F, Fochetti R, Magliano A, Vinciguerra V, Ermenegildi A, De LC. First report of Culicoides biting midges (Diptera: Ceratopogonidae) attacking people in Italy, with the description of extreme larval breeding sites and diurnal activity of Culicoides riethi. J Med Entomol. 2022;59:772–6.
Vasic A, Zdravkovic N, Anita D, Bojkovski J, Marinov M, Mathis A, et al. Species diversity, host preference and arbovirus detection of Culicoides (Diptera: Ceratopogonidae) in south-eastern Serbia. Parasit Vectors. 2019;12:61.
Brugman VA, Hernandez-Triana LM, England ME, Medlock JM, Mertens PP, Logan JG, et al. Blood-feeding patterns of native mosquitoes and insights into their potential role as pathogen vectors in the Thames estuary region of the United Kingdom. Parasit Vectors. 2017;10:163.
Martínezdela Puente J, Mendez M, Ruiz S, Godoy JA, Soriguer RC, Figuerola J. Individual identification of endangered species using mosquito blood meals: a proof-of-concept study in Iberian lynx. Parasitol Res. 2015;114:1607–10.
The authors thank the “Unidad de Anillo Verde y Biodiversidad y Centro de Estudios Ambientales (CEA) del Ayuntamiento de Vitoria-Gasteiz” and Urdaibai Bird Center for allowing this work to proceed in the Salburua and Urdaibai aquatic ecosystems, respectively.
This work was funded by the Department of Economic Development and Infrastructures of the Basque Government. Mikel A. González was supported by a 3-year postdoctoral fellowship from the Department of Education of the Basque Government. Fatima Goiri was supported by a PhD award from the Department of Economic Development and Infrastructures of the Basque Government (FPI-2019). Aitor Cevidanes is supported by a ‘Juan de la Cierva-Formación’ post-doctoral grant FJC2019-041737-I funded by MCIN/AEI/10.13039/501100011033. Luis Hernández-Triana’s contribution was supported by the Department for Environment, Food and Rural Affairs (Defra), The Scottish Government and Welsh Government through grant SV3045. DNA sequencing and molecular data analysis was enabled by funding to PDNH from the Food from Thought Program with support from the Canada First Research Excellence Fund.
Ethics approval and consect to participate
Consent for publication
The authors declare that they have no competing interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
González, M.A., Goiri, F., Prosser, S.W.J. et al. Culicoides species community composition and feeding preferences in two aquatic ecosystems in northern Spain. Parasites Vectors 15, 199 (2022). https://doi.org/10.1186/s13071-022-05297-5