The tree that hides the forest: cryptic diversity and phylogenetic relationships in the Palaearctic vector Obsoletus/Scoticus Complex (Diptera: Ceratopogonidae) at the European level
Parasites & Vectors volume 13, Article number: 265 (2020)
Culicoides obsoletus is an abundant and widely distributed Holarctic biting midge species, involved in the transmission of bluetongue virus (BTV) and Schmallenberg virus (SBV) to wild and domestic ruminants. Females of this vector species are often reported jointly with two morphologically very close species, C. scoticus and C. montanus, forming the Obsoletus/Scoticus Complex. Recently, cryptic diversity within C. obsoletus was reported in geographically distant sites. Clear delineation of species and characterization of genetic variability is mandatory to revise their taxonomic status and assess the vector role of each taxonomic entity. Our objectives were to characterize and map the cryptic diversity within the Obsoletus/Scoticus Complex.
Portion of the cox1 mitochondrial gene of 3763 individuals belonging to the Obsoletus/Scoticus Complex was sequenced. Populations from 20 countries along a Palaearctic Mediterranean transect covering Scandinavia to Canary islands (North to South) and Canary islands to Turkey (West to East) were included. Genetic diversity based on cox1 barcoding was supported by 16S rDNA mitochondrial gene sequences and a gene coding for ribosomal 28S rDNA. Species delimitation using a multi-marker methodology was used to revise the current taxonomic scheme of the Obsoletus/Scoticus Complex.
Our analysis showed the existence of three phylogenetic clades (C. obsoletus clade O2, C. obsoletus clade dark and one not yet named and identified) within C. obsoletus. These analyses also revealed two intra-specific clades within C. scoticus and raised questions about the taxonomic status of C. montanus.
To our knowledge, our study provides the first genetic characterization of the Obsoletus/Scoticus Complex on a large geographical scale and allows a revision of the current taxonomic classification for an important group of vector species of livestock viruses in the Palaearctic region.
In 2006, northern Europe faced massive outbreaks of bluetongue disease (BTV), a Culicoides-borne viral infection which affects wild and domestic ruminants. This was followed by the emergence of Schmallenberg virus (SBV) in 2011, another Culicoides-borne virus, which also caused important economic losses for farmers of domestic ruminants . Both, the emergence and massive spread of these diseases over the Palaearctic region raised questions about the vector competence of native Palaearctic biting midges members of the genus Culicoides [2, 3]. Quickly following these epizootics, studies confirmed Culicoides species of the subgenus Avaritia as the most likely vector species and particularly those of the Obsoletus group (see below) [3,4,5]. Culicoides obsoletus, C. scoticus, C. dewulfi and C. chiopterus are reported to be implicated in BTV and SBV transmission in Europe based on virus detection or isolation in field-collected populations . Laboratory experimental infections have also confirmed the vector competence of C. obsoletus and C. scoticus, highlighting species variations in their competence level .
The literature defines the Obsoletus Group as a group of species with a similar morphology, especially for the characters commonly used for the identification of these insect vectors, namely wing spot pattern (poorly defined spotted wings and a second radial cell with a light spot) and distinctive male genitalia [8,9,10]. Adults of species in this Group are known to be abundant, widespread across central and northern Europe, and are characterized by long seasonal occurrence [11, 12]. At present, the group is an artificial taxonomic entity still poorly defined, with no real consensus on the included species, with variable internal groupings and naming. Indeed, the world catalogue of Culicoides does not account for levels below subgenus and does not identify species complexes as groups . At present, the Obsoletus Group is composed of nine valid species: C. obsoletus (Meigen), 1818; C. sinanoensis Tokunaga, 1937; C. scoticus Downes & Kettle, 1952; C. montanus Shakirzjanova, 1962; C. gornostaevae Mirzaeva, 1984; C. abchazicus Dzhafarov, 1964; C. filicinus Gornostaeva & Gachegova, 1972; C. alachua Jamnback & Wirth, 1963; and C. sanguisuga (Coquillet, 1901). The latter two are the only species belonging to this group exclusively present in the Nearctic region, while the others, are sympatric in the Palaearctic region . Culicoides obsoletus is considered Holarctic because it is present in both Nearctic and Palaearctic regions . Combined studies of the geometric morphometry of wings, coupled with molecular analysis have excluded C. chiopterus and C. dewulfi from the Obsoletus Group, despite previously being considered as part of the Group, based only on morphological features [9, 10, 14,15,16]. Life-cycles and trophic behaviors for most of the species in the Obsoletus Group are not well described or vary greatly. For instance, ecological niches of C. chiopterus and C. obsoletus are suspected to be different, although these two species are phylogenetically very close within the subgenus Avaritia. Culicoides obsoletus is a widespread generalist species and occupies a wide range of larval habitats such as forest litter, silage residue, tree holes or manure . Culicoides obsoletus shows opportunistic trophic preferences and is able to take blood meals on various hosts species (man, sheep, goat, cow, horse or rodent) and, occasionally, on birds . Culicoides chiopterus is a more specialist species, found engorged almost exclusively on cattle blood with larvae associated with cattle dung [19, 20].
In addition to the group, Meiswinkel et al. , defined the Obsoletus Complex to close taxa with very similar female adult morphology, including C. obsoletus, C. montanus and C. scoticus. Several authors have recently reported the existence of cryptic diversity within C. obsoletus, namely the C. obsoletus clade ‘O2’ and C. obsoletus clade ‘O3’ in Sweden and Switzerland , and the C. obsoletus clade ‘dark’ in the Netherlands [21, 22]. We will use the term “Obsoletus/Scoticus Complex” here to refer to the cryptic species C. obsoletus, C. scoticus and C. montanus as well as all previously described operational taxonomic units in the literature (see above). The terminology “Obsoletus/Scoticus Complex” is written according to the rules defined by Harbach  for infrasubgeneric categories within the genus Anopheles. Considering that sympatric cryptic species may exhibit different vector competence and may confound epidemiological investigations, it is mandatory to assess the intra- and interspecific diversity within the Obsoletus/Scoticus Complex.
Given the difficulty of separating these species based on morphological identification, various molecular markers have been used to overcome specific identification problems, including 16S ribosomal DNA [24, 25], 28S ribosomal DNA , cytochrome oxidase b (cytb), the internal transcribed spacer region 1 (ITS1)  and ITS2 rDNA . However, the DNA region primarily used to infer phylogenetic relationships in this complex has been the mitochondrial cytochrome c oxidase subunit 1 (cox1) . Most of the diversity within the Obsoletus/Scoticus Complex has been identified using cox1 [29,30,31,32,33]. Despite the diversity of markers used to characterize the phylogeny of the complex, few studies have used a multi-marker approach . This study, therefore, will integrate a multi-marker approach in order to strengthen genetic reconstruction of the Obsoletus/Scoticus Complex.
The Obsoletus/Scoticus Complex, as currently presented in the literature, is confused and needs taxonomic revision. We characterized and mapped the genetic diversity of the Obsoletus/Scoticus Complex along a Palaearctic-Mediterranean transect covering Scandinavia to Canary islands (North to South) and Canary islands to Turkey (West to East). Our main objectives were to identify and describe the cryptic diversity observed within the Obsoletus/Scoticus Complex over a wide geographical area using molecular analyses and to question the taxonomic status of some newly described clades. In order to achieve these objectives, we conducted a molecular analysis which combined multi-marker sequencing, phylogenetic analyses and species delimitation to explore the genetic diversity of the Obsoletus/Scoticus Complex in the western European portion of the Palaearctic region.
Culicoides capture and morphological identification
Biting midges were collected at 68 sites located in 20 countries in the western European portion of the Palaearctic region, between 2009 and 2017, using national surveillance networks for Culicoides populations or local collections (Additional file 1: Table S1). Collections were made overnight with Onderstepoort Veterinary Institute (OVI) light traps set at farms near horses, cattle or sheep and all insects were stored in 70% ethanol. Morphological identification to the species level of adult Culicoides spp. was performed under a binocular microscope using the available identification keys [35, 36].
DNA extraction, amplification and sequencing
DNA was extracted from a total of 3883 adult females belonging to the Obsoletus/Scoticus Complex using the NucleoSpin® DNA kit RapidLyse (Macherey-Nagel, Duren, Germany), following the manufacturer’s instructions. An additional step was added, before extraction, for all individuals (specimens were ground in 50 μl of 1× PBS buffer). DNA samples are available upon request. Fragments of cox1 were amplified for the 3883 individuals. After sequence cleaning, 3763 sequences of cox1 were obtained [dataset cox1 (1)].Fragments of 16S and 28S rDNA were amplified on individuals chosen to be representative of the entire species diversity resulting from cox1 [dataset cox1 (2)] to reinforce mitochondrial gene sequences. All primer sequences as well as the information relating to them are present in Additional file 2: Table S2. PCR’s were performed in a 25 μl reaction volume. The PCR mix contained 1× Qiagen buffer, 1 mM MgCl2, 0.25 mM of each dNTP, 0.2 μM of each primer, 1.25 U Qiagen Taq Polymerase and 0.7 ng/μl genomic DNA for all genes. PCR programs included one-step of 5 cycles before a second step with 35 cycles for 16S rDNA and 28S rDNA. PCR amplification conditions were: an initial denaturation step at 94 °C for 5 min followed by 5 cycles of 94 °C for 30 s; 45 °C for cox1, 42 °C for 16S rDNA or 55 °C for 28S rDNA for 40 s; 72 °C for 1min; 35 cycles of 94 °C for 30 s; 51 °C for cox1, 55 °C for 16S rDNA or 50 °C for 28S rDNA for 30 s; 72 °C for 1 min; and a final extension step at 72 °C for 10 min. For each amplification reaction, negative controls were carried out. The PCR products were visualized on 1.5% agarose gels with a GelRed® Nucleic Acid Gel Stain, staining after migration of 90 min at 130 V by electrophoresis for quality control. After purifications, carried out by the sequencing service provider, the remaining 20 μl were sequenced with the same forward primers used for PCR (https://www.genewiz.com).
The reference sequences of cox1 used to identify individuals to species are available in Additional file 3: Table S3. A total of 3763 cox1 sequences (Additional file 4: Table S4) from female adults morphologically identified as belonging to the Obsoletus/Scoticus Complex were obtained after deletion of short and poor-quality sequences. The cox1 alignment was used to identify all Culicoides to species- or clade-level within the complex using the reference sequences. Among the reference sequences used to specifically assign our Culicoides samples, some sequences previously identified as C. obsoletus O1 and O3  were included in our analysis in order to cover the diversity of the clades described in the literature within the Obsoletus/Scoticus Complex. After comparison with other sequences, it appears that the sequences named C. obsoletus O1 were actually C. obsoletus and that C. obsoletus O3 belongs to C. obsoletus clade dark. For this purpose, cox1 sequences were aligned with reference sequences [dataset cox1 (1)]. A phylogenetic tree based on maximum likelihood method allowed designation of a species name to each sequence if the latter belonged to a monophyletic clade with strong support (bootstrap < 900) that included a reference sequence. Thus, Culicoides were sequenced for markers 16S rDNA and 28S rDNA to support the phylogenetic reconstruction of the complex. All cox1 [dataset cox1 (2)], 16S rDNA and 28S rDNA sequences were independently aligned with the MUSCLE  algorithm available in the software GENEIOUS v.6.0.5 (Biomatters, http://www.geneious.com). Genetic diversity indices, haplotype and nucleotide diversity were evaluated using DNASP v.5.10 . Alignments with gaps were cleaned using the software Gblocks 0.91b . To assess genetic distance between clades and species within, barcoding gap bar chart using R software was performed with ggplot2  and ggthemes packages. Intra- and interspecific genetic differences based on the Kimura 2-Parameter (K2P) distance model [41, 42] were calculated with MEGAX . In order to map the specific diversity of the Obsoletus Group, the R software version 3.6.0 was used with the Leaflet version 2.0.2 and shiny packages version 1.4.0.
Phylogenetic trees were constructed for the three markers using maximum-likelihood (ML) and Bayesian inference (BI). Bayesian inference analyses were conducted on MrBayes version 3.2.6, with tree sampling every 1000 generations in order to calculate posterior probabilities (PP) and 10 million generations. Optimal sequence evolutionary models for each analysis were obtained with Bayesian information criterion (BIC) using jModelTest. Maximum-likelihood analyses were conducted on PhyML 3.0. The ML analyses were conducted with the best model selected using 1000 bootstrap replicates for each dataset to investigate the level of support at each node, with starting tree determined by BioNJ analysis. After independent analysis of each gene, alignment of cox1 and 16S rDNA were concatenated, and analyzed following ML and BI methods.
Species delimitation methods
Two species delimitation methods were applied. The first method was a Bayesian implementation of classical GMYC method, Bayesian General Mixed Yule Coalescent (bGMYC) . The single-locus ultrametric gene trees used for bGMYC methods were created with BEAST 1.8.0  under a strict clock model, a Yule Process Tree Model of speciation, and a random starting tree. This analysis was carried out with default prior distribution, without outgroups and with 10 million generations sampled every 1000 cycles with HKY + G substitution model for 16S rDNA, and with HKY + I for 28S rDNA . The software TREEANNOTATOR v1.8.2 was used to find Ultrametric maximum clade credibility (MCC). Single-threshold GMYC analyses were conducted with splits package in R.
The second species delimitation method used was the Bayesian Poisson Tree Processes (bPTP) through the web server PTP (http://species.h-its.org/ptp/)  with 100,000 MCMC generations and a thinning parameter of 100 on a maximum likelihood phylogenetic tree constructed with cox1 and 16S rDNA genes concatenated with C. dewulfi as outgroup.
In total, 3763 sequences were obtained for cox1, 95 for 16S rDNA and 95 for 28S rDNA (Table 1). All sequences were deposited in GenBank (Additional file 4: Table S4). No stop codons, insertions or deletions were found in any of the cox1 sequences, indicating functional mitochondrial products.
Within the selected samples present in our data set, the most abundant species was C. obsoletus, with 2416 individuals sampled, representing 68% of all Culicoides caught (Fig. 1). Species diversity within the complex varied according to the latitudes of the sampling sites (Fig. 1). The most sampled species in northern Europe (Norway and Finland) was the C. obsoletus clade O2 with 62% (162 individuals) and 58% (147 individuals) of this clade, respectively, in each country. For eastern Europe (Latvia, Poland, Serbia, Bulgaria, Macedonia, Greece and Turkey) the most sampled species was C. obsoletus, representing, for example, 99% (162 individuals) of all Culicoides sampled in Greece. However, a population in Poland (Wronka) appears to be an exception with 100% (35 individuals) of C. obsoletus clade O2. Three individuals, belonging to a phylogenetic clade unidentifiable by our reference sequences close to C. obsoletus dark, were reported from Latvia. The latter sequences are identical to sequences present in the BOLD database (accession numbers: GMGRC1056-13, GMGRC1000-13, GMGRD2587-13) of Culicoides collected in Bavaria, Germany. Finally, western and central Europe (Portugal, Spain, UK, France, Italy, Netherlands, Germany, Switzerland and Denmark) had the higher species diversity of all the species found in the Palaearctic transect. The most sampled species within this area was C. obsoletus with 913 individuals, representing 64% of the Culicoides sampled. However, unlike in eastern Europe, C. scoticus clade 1 was also found in significant numbers with 386 individuals, or 27% of the samples. Culicoides obsoletus clade dark was rarely reported in Europe with only 26 individuals found in France, Denmark, Finland, Italy, Latvia, Norway and Switzerland. Culicoides montanus was found in relatively high proportion in Morocco with 54% (80 individuals) of samples, whereas it was much rarer and more sporadic in European countries.
A total of 228 different cox1 haplotypes were identified. Haplotype diversity varied from 0.829 for C. obsoletus clade O2 to 0.517 for C. scoticus clade 1 and C. obsoletus clade dark (Table 2). In accordance with mitochondrial genes in insects , the GC composition ranged from 16% for 16S rDNA alignment to 39.8% for 28S rDNA, with a strong AT bias (Table 2). Culicoides obsoletus showed very high intraspecific diversity with 94 different haplotypes, representing nearly 43% of the total haplotypic diversity encountered in this study. The maximum interspecific genetic distance (Fig. 2) was reached between C. dewulfi and other taxonomic units, with a minimum of 17% of genetic distance between this outgroup and all other members inside the complex. Culicoides dewulfi is used here as an outgroup. Genetic distances of the same level as the other intraspecific distances were observed between C. scoticus clade 1 and C. scoticus clade 2, and between C. obsoletus and C. montanus, with a maximum of 2% and 4%, respectively. Similar interspecific genetic distances, were observed between all the other clades within the complex, with a minimum of 8% distance. All species had mean intraspecific distances of less than 1%, other than C. scoticus clade 2.
Among the 228 unique cox1 haplotypes, 95 were selected to represent the specific diversity of the complex. Using the same individuals, two alignments of 95 sequences of 16S rDNA (Additional file 5: Figure S1.) and 28S rDNA were constructed. In order to compare tree topology and to concatenate markers, a second cox1 alignment of the 95 sequences was performed (Additional file 6: Figure S2.). The cox1 dataset had a greater haplotype diversity than 16S rDNA and 28S rDNA datasets, with 0.994 vs 0.851 and 0.911 respectively. 28S rDNA was more monomorphic than cox1 and 16S rDNA, with 25, 146 and 51 polymorphic sites, respectively.
Information on the alignments used for the construction of phylogenetic trees is provided in Table 1. Species with confirmed taxonomic validity (C. montanus, C. obsoletus and C. scoticus) and cryptic taxa (C. obsoletus clade O2, C. obsoletus clade dark) highlighted by previous studies were strongly supported (bootstrap > 90%) (Fig. 3). Culicoides obsoletus clade O2 constituted a monophyletic clade with strong support (bootstrap > 90%). Culicoides obsoletus and C. montanus formed a monophyletic clade. Culicoides scoticus, another species considered valid, showed two phylogenetic clades, i.e. C. scoticus clade 1 and C. scoticus clade 2. A monophyletic clade close to C. obsoletus clade dark, was strongly supported by bootstrap values. Topologies of phylogenetics trees constructed via maximum likelihood (Fig. 3a) and Bayesian inference (BI) analyses (Fig. 3b) were congruent for alignment of full haplotype diversity.
Using the cox1 dataset and bGMYC method for species delimitation (Fig. 4), 6 molecular operational taxonomic units (MOTUs) were observed: C. montanus, C. obsoletus, C. scoticus, C. obsoletus clade O2, C. obsoletus clade dark and Culicoides sp. within the Obsoletus/Scoticus Complex. Based on the 16S rDNA dataset analysed with the same delimitation method, three MOTUs were characterised including (i) two clades within C. scoticus; (ii) Culicoides sp., C. obsoletus clade dark and C. obsoletus clade O2; and (iii) C. obsoletus and C. montanus. Species delimitation with the 28S rDNA and 16S rDNA dataset had a lower resolution compared to cox1. The 28S rDNA and 16S rDNA datasets were invaluable in identifying cryptic diversity within the Obsoletus/Scoticus Complex. Indeed, the tree generated with the 28S rDNA dataset showed very low polymorphisms and resolution signal (Additional file 7: Figure S3.). The bPTP method was conducted on concatenated dataset of cox1 and 16S rDNA genes. The best statistical support for species delimitation was for concatenated dataset with bPTP. MOTUs found with the concatenated alignment and cox1 were the same, except for the delimitation of C. montanus and C. obsoletus. Indeed, using molecular delineation based on bPTP analysis, MOTUs were observed for C. scoticus, C. obsoletus clade O2, C. obsoletus clade dark, and Culicoides sp. Using the bPTP method we were not able to distinguish C. obsoletus and C. montanus as well as C. scoticus clades 1 and 2, which were not classified as separate species by either method.
Our study investigated the genetic diversity of the Obsoletus/Scoticus Complex at the European level (68 sampling sites across 20 Palaearctic countries with 3763 cox1 sequences). For the first time, the sample selection covers the whole known western Palaearctic distribution area of these species. The complementary use of the mitochondrial 16S rDNA and nuclear 28S rDNA genes confirms the important level of cryptic diversity found within the Obsoletus/Scoticus Complex. Indeed species delimitation methods allowed us to delineate five MOTUs and: (i) to provide evidence of the taxonomic validity of C. obsoletus clade O2 and C. obsoletus clade dark (ii) to identify individuals belonging to a species not yet described or not present into the databases; and (iii) to question the taxonomic status of C. montanus.
Species assignment at the European scale showed variation in the distribution of the cryptic diversity of Obsoletus/Scoticus Complex. This result confirmed a previous study by Möhlmann et al. , who found a strong latitudinal effect on the relative abundance of species of the Obsoletus/Scoticus Complex. However, the previous study was carried out with few individuals from a relatively small number of countries and sampling sites . This contrasts with our study, the first to be conducted at a European scale with a large data set sufficient to provide a more precise idea of the cryptic diversity within the Obsoletus/Scoticus Complex. The latitudinal variation in the relative abundance of the different cryptic species in the Obsoletus/Scoticus Complex, could be due to a wide range of factors like different ecological niches, or differences in the availability of hosts and breeding sites . For example, a study conducted in Italy showed that C. scoticus collection sites were dominated by areas of natural vegetation or forest, at medium altitudes, preferably in wilder and more pristine environments . However, the heterogeneity of Culicoides collection dates may also explain these variations in specific diversity within the complex. For instance, a study conducted in Sweden, found a seasonal variation in Culicoides community structure . The omnipresence of C. obsoletus makes it the dominant species in Europe, confirming its status as a generalist species, which tolerates a wide range of eco-climatic conditions. The dominant species also varies according to geographical location. For example, C. obsoletus clade O2 is the most sampled species in Nordic countries while C. montanus prevails in Morocco. France, Italy and Spain appear to have the highest specific diversity of Culicoides belonging to the Obsoletus/Scoticus Complex. These three countries bring together all the cryptic diversity known so far, except the new, unidentified, Culicoides taxon that has only been found in Latvia regarding our dataset. This could be due to the significant diversity of ecological niches as well as the high density of hosts in these countries and this species diversity variation is in line with general patterns of latitudinal increase in species richness .
Our phylogenetic analysis allowed us to define seven well supported phylogenetic clades. Some of them correspond to species with taxonomic validity (C. montanus, C. obsoletus and C. scoticus), some more recently described phylogenetic clades (C. obsoletus clade dark and C. obsoletus clade O2) [52, 54] and some clades never described before (Culicoides sp. and C. scoticus clade 2).
Without taking into account C. scoticus, the phylogenetic reconstruction produced herein, confirmed the presence of two divergent groups; one consisting of C. obsoletus and C. montanus and the other of C. obsoletus clade dark, C. obsoletus clade O2 and a clade not yet described in the literature. We were unable to identify these sequences due to the absence of reference sequences identified at the species level in the sequence databases. However, three sequences, from specimens collected in Bavaria in Germany, present in the BOLD database are identical. Further sampling and sequencing of Culicoides from eastern Europe are necessary in order to associate morphological features with this new cryptic species. The large number of Culicoides processed during this study made it difficult to use non-destructive DNA extraction techniques. Non-destructive techniques would be necessary to couple morphological criteria with genetic analysis, in order to identify this species.
We also described a second clade phylogenetically very close to C. scoticus, C. scoticus clade 2 . However, given the small genetic distance observed between these two clades, they can be considered as intraspecific variation within C. scoticus.
According to this study, C. obsoletus clade dark appears as a true cryptic species, with high phylogenetic support. Meiswinkel et al.  hypothesized that C. obsoletus clade dark could be C. gornostaevae Mirzaeva, 1984, but C. gornostaevae was significantly larger and had a distribution restricted to the boreal zone of Siberia. However, this species has been recently reported from Norway, Poland and Sweden  but the lack of C. gornostaevae reference sequences in publicly available databases prevented comparisons with C. obsoletus clade dark.
Culicoides obsoletus clade O2 was also strongly supported by phylogenetic and species delimitation analysis, and thus could be considered a cryptic species within the Obsoletus/Scoticus Complex. However, our sample showed a high abundance of C. obsoletus clade O2 mainly at high latitudes, whereas it had initially been identified for the first time further south in the Swiss Alps [54, 56, 57] and France.
Little is known about the ecology of C. montanus. We found this species over a large geographical area, from Morocco to a few individuals in Norway. In spite of this, we have identified only a very small number of C. montanus except in Morocco where this species constitutes the majority. This latter result is in accordance with the fact that sites in Italy, where C. montanus is predominant are characterised by a high land surface temperature, higher than the values registered in the C. obsoletus and C. scoticus sites . Culicoides obsoletus and C. scoticus are considered to be sibling species because of diagnostic female morphological characters, which are difficult to observe or overlap , as opposed to the morphological identification of males which is simpler . If the genetic analyses by Pages & Sarto  confirm that both C. obsoletus and C. scoticus are distinct species, the question remains open for C. montanus. Previous phylogenetic studies based on cox1 indeed showed limited genetic distance between C. obsoletus and C. montanus and in studies based on ITS sequences C. montanus always appears in one of the subclades of C. obsoletus . In our multi-marker phylogenetic tree, C. montanus and C. obsoletus formed a monophyletic clade. In addition, the genetic distance between C. obsoletus and C. montanus was of the same order of magnitude as some intraspecific distances.
The number of putative species defined within the Obsoletus/Scoticus Complex varied depending on the molecular markers and species delimitation methods used (Fig. 4), particularly pertaining to the status of C. obsoletus and C. montanus. Unlike the other methods, the bGMYc method, based on cox1, distinguished these two species. This can be explained by the fact that some parameters of the analysis (i.e. priors), like differences in population size or speciation rates, can bias the GMYC method by overestimating the number of species [46, 60,61,62]. Insufficient sampling, high gene flow or a recent speciation event are also likely explanations for the differences in results between phylogenetic trees and species delimitation results [63, 64]. Although subject to the same constraints, it has been shown that the bPTP method is significantly more robust . Moreover, methods based on multiple loci improve discovery, resolution and stability of species delimitation [66, 67]. Furthermore, studies have shown that species delineation in insects is more appropriate with multilocus species delimitation methods [68,69,70,71]. These arguments allowed us to validate the species delimitation scheme produced with the bPTP method. This analysis coupled with the low level of genetic distances observed between C. montanus and C. obsoletus led us to question the taxonomic status of C. montanus. This could be the beginning of a speciation process of C. montanus within the clade of C. obsoletus. Hovewer, C. montanus was originally described from Kazakhstan and present in central Asia . Therefore, examination of individuals sampled in this geographical area could strengthen our conclusions.
In the light of our conclusions, in-depth morphological analyses with the deposition of reference individuals will have to be carried out in order to decide whether or not to definitively rule out the taxonomic status of the cryptic species making up this complex. Indeed, although adult females are not morphologically distinguishable on a routine basis, males are easier to identify using their genitalia and pupal differences and can provide evidence of morphological differences between the species.
This study provides clarification of the distribution pattern of species belonging to the Obsoletus/Scoticus Complex, using a dataset based on samples from western Palaeartic and Mediterranean transect. Strong variations in latitudinal cryptic species diversity was observed. This study clarifies the phylogenetic relationships between species belonging to the Obsoletus/Scoticus Complex. We identified and validated five MOTUs, C. obsoletus, C. scoticus, C. obsoletus clade O2, C. obsoletus clade dark and a MOTUs corresponding to an unidentified species. The latter three species have not been formally described but our results confirm that they should be considered as species in their own right. More detailed studies of their morphology and ecology are needed to provide more detailed descriptions of these species. Furthermore, our results raise questions concerning the taxonomic status of C. montanus, which was previously considered as a taxonomically valid species.
Availability of data and materials
All data generated or analysed during this study are included in this published article and its additional files. The newly generated sequences were submitted in the GenBank database under the accession numbers MT170026-MT173788.
- cytb :
cytochrome oxidase b
internal transcribed spacer region 1
internal transcribed spacer region 2
cytochrome c oxidase subunit 1
number of individuals
number of haplotypes
haplotype (gene) diversity
number of polymorphic sites
Bayesian General Mixed Yule Coalescent
molecular operational taxonomic units
Bayesian Poisson Tree Processes
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.
Wilson AJ, Mellor PS. Bluetongue in Europe: past, present and future. Phil Trans R Soc Lond B. 2009;364:2669–81.
Carpenter S, Wilson A, Mellor PS. Culicoides and the emergence of bluetongue virus in northern Europe. Trends Microbiol. 2009;17:172–8.
Hoffmann B, Eschbaumer M, Beer M. Real-time quantitative reverse transcription-PCR assays specifically detecting bluetongue virus serotypes 1, 6, and 8. J Clin Microbiol. 2009;47:2992–4.
Savini G, MacLachlan NJ, Sanchez-Vizcaino JM, Zientara S. Vaccines against bluetongue in Europe. Comp Immunol Microbiol Infect Dis. 2008;31:101–20.
Carpenter S, Lunt HL, Arav D, Venter GJ, Mellor PS. Oral susceptibility to bluetongue virus of Culicoides (Diptera: Ceratopogonidae) from the United Kingdom. J Med Entomol. 2006;43:73–8.
Carpenter S, Mcarthur C, Selby R, Ward R, Nolan DV, Luntz AJ, et al. Experimental infection studies of UK Culicoides species midges with bluetongue virus serotypes 8 and 9. Vet Rec. 2008;163:589–92.
Mathieu B. Les espèces de Culicoides du sous-genre Avaritia (Diptera: Ceratopogonidae) dans le monde: révision systématique et taxonomique des espèces d’intérêt dans la transmission dʼOrbivirus. École doctorale Sciences de la vie et de la santé Université de Strasbourg; 2011.
Meiswinkel R, Gomulski LM, Delecolle JC, Goffredo M, Gasperi G. The taxonomy of Culicoides vector complexes—unfinished business. Vet Ital. 2004;40:151–9.
Hajd Henni L, Sauvage F, Ninio C, Depaquit J, Augot D. Wing geometry as a tool for discrimination of Obsoletus group (Diptera: Ceratopogonidae: Culicoides) in France. Infect Genet Evol. 2014;21:110–7.
Mehlhorn H, Walldorf V, Klimpel S, Jahn B, Jaeger F, Eschweiler J, et al. First occurrence of Culicoides obsoletus-transmitted bluetongue virus epidemic in central Europe. Parasitol Res. 2007;101:219–28.
Savini G, Goffredo M, Monaco F, Di Gennaro A, Cafiero MA, Baldi L, et al. Bluetongue virus isolations from midges belonging to the Obsoletus complex (Culicoides, Diptera: Ceratopogonidae) in Italy. Vet Rec. 2005;157:133–9.
Borkent A. World species of biting midges (Diptera: Ceratopogonidae). Salmon Arm: American Museum of Natural History, and Instituto Nacional de Biodiversidad; 2016.
Gomulski LM, Meiswinkel R, Delecolle JC, Goffredo M, Gasperi G. Phylogenetic relationships of the subgenus Avaritia Fox, 1955 including Culicoides obsoletus (Diptera, Ceratopogonidae) in Italy based on internal transcribed spacer 2 ribosomal DNA sequences. Syst Entomol. 2005;30:619–31.
Ander M, Troell K, Chirico J. Barcoding of biting midges in the genus Culicoides: a tool for species determination. Med Vet Entomol. 2013;27:323–31.
Kluiters G, Pages N, Carpenter S, Gardes L, Guis H, Baylis M, et al. Morphometric discrimination of two sympatric sibling species in the Palaearctic region, Culicoides obsoletus Meigen and C. scoticus Downes & Kettle (Diptera: Ceratopogonidae), vectors of bluetongue and Schmallenberg viruses. Parasit Vectors. 2016;9:262.
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.
Martinez-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.
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.
Kettle DS, Lawson JWH. The early stages of British biting midges Culicoides Latreille (Diptera: Ceratopogonidae) and allied Genera. Bull Entomol Res. 2009;43:421–67.
Meiswinkel R, De Bree F, Bossers-De Vries R, Elbers AR. An unrecognized species of the Culicoides obsoletus complex feeding on livestock in The Netherlands. Vet Parasitol. 2015;207:324–8.
Elbers AR, Meiswinkel R. Culicoides (Diptera: Ceratopogonidae) host preferences and biting rates in the Netherlands: comparing cattle, sheep and the black-light suction trap. Vet Parasitol. 2014;205:330–7.
Harbach RE. Review of the internal classification of the genus Anopheles (Diptera, Culicidae) - the foundation for comparative systematics and phylogenetic research. Bull Entomol Res. 1994;84:331–42.
Debeila TJ. Characterisatoin of selected Culicoides population in South Africa using genetic markers. PhD Thesis, Faculty of Veterinary Science, University of Pretoria; 2010.
Tabachnick W. Culicoides and the global epidemiology of bluetongue virus infection. Vet Ital. 2004;40:144–50.
Bakhoum MT, Labuschagne K, Huber K, Fall M, Mathieu B, Venter G, et al. Phylogenetic relationships and molecular delimitation of Culicoides Latreille (Diptera: Ceratopogonidae) species in the Afrotropical region: interest for the subgenus Avaritia. Syst Entomol. 2018;43:355–71.
Cetre-Sossah C, Baldet T, Delecolle JC, Mathieu B, Perrin A, Grillet C, et al. Molecular detection of Culicoides spp. and Culicoides imicola, the principal vector of bluetongue (BT) and African horse sickness (AHS) in Africa and Europe. Vet Res. 2004;35:325–37.
Gomulski LM, Meiswinkel R, Delecolle JC, Goffredo M, Gasperi G. Phylogeny of the subgenus Culicoides and related species in Italy, inferred from internal transcribed spacer 2 ribosomal DNA sequences. Med Vet Entomol. 2006;20:229–38.
Nolan DV, Carpenter S, Barber J, Mellor PS, Dallas JF, Mordue Luntz AJ, et al. Rapid diagnostic PCR assays for members of the Culicoides obsoletus and Culicoides pulicaris species complexes, implicated vectors of bluetongue virus in Europe. Vet Microbiol. 2007;124:82–94.
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.
Pages N, Munoz-Munoz F, Talavera S, Sarto V, Lorca C, Nunez 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.
Schwenkenbecher JM, Mordue AJ, Piertney SB. Phylogenetic analysis indicates that Culicoides dewulfi should not be considered part of the Culicoides obsoletus complex. Bull Entomol Res. 2009;99:371–5.
Sarto i Monteys V, Ventura D, Pages N, Aranda C, Escosa R. Expansion of Culicoides imicola, the main bluetongue virus vector in Europe, into Catalonia, Spain. Vet Rec. 2005;156:415–7.
Harrup LE, Bellis GA, Balenghien T, Garros C. Culicoides Latreille (Diptera: Ceratopogonidae) taxonomy: current challenges and future directions. Infect Genet Evol. 2015;30:249–66.
Delecolle JC. Nouvelle contribution à lʼétude systématique et iconographique des espèces du genre Culicoides (Diptéra: Cératopogonidae) du Nord-Est de la France. PhD Thesis, University of Strasbourg, Strasbourg; 1985.
Campbell JA, Pelham-Clinton EC. Taxonomic review of the British species of Culicoides Latreille (Diptera, Ceratopogonidae). Proc R Entomol Soc. 1960;67:181–302.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7.
Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25:1451–2.
Castresana J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000;17:540–52.
Wickham H. ggplot2: Elegant graphics for data analysis. J R Stat Soc A Stat. 2009;174:245–50.
Kimura M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980;16:111–20.
Nei M, Rooney AP. Concerted and birth-and-death evolution of multigene families. Annu Rev Genet. 2005;39:121–52.
Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35:1547–9.
Reid NM, Carstens BC. Phylogenetic estimation error can decrease the accuracy of species delimitation: a Bayesian implementation of the general mixed Yule-coalescent model. BMC Evol Biol. 2012;12:196.
Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29:1969–73.
Fujisawa T, Barraclough TG. Delimiting species using single-locus data and the Generalized Mixed Yule Coalescent approach: a revised method and evaluation on simulated data sets. Syst Biol. 2013;62:707–24.
Zhang J, Kapli P, Pavlidis P, Stamatakis A. A general species delimitation method with applications to phylogenetic placements. Bioinformatics. 2013;29:2869–76.
Crozier RH, Crozier YC. The mitochondrial genome of the honeybee Apis mellifera: complete sequence and genome organization. Genetics. 1993;133:97–117.
Möhlmann TWR, Bekendam AM, van Kemenade I, Wennergren U, Favia G, Takken W, et al. Latitudinal diversity of biting midge species within the Obsoletus group across three habitats in Europe. Med Vet Entomol. 2019;33:420–6.
Goffredo M, Meiswinkel R, Federici V, Di Nicola F, Mancini G, Ippoliti C, et al. The ‛Culicoides obsoletus groupʼ in Italy: relative abundance, geographic range, and role as vector for bluetongue virus. Vet Ital. 2016;52:235–41.
Ippoliti C, Goffredo M, Salini R, Bocci M, Quaglia M, Federici V, et al. Eco-climatic indicators for three Culicoides species of the Obsoletus complex in Italy. Vet Ital. 2016;52:213–22.
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.
Allen AP, Gillooly JF. Assessing latitudinal gradients in speciation rates and biodiversity at the global scale. Ecol Lett. 2006;9:947–54.
Wenk CE, Kaufmann C, Schaffner F, Mathis A. Molecular characterization of Swiss Ceratopogonidae (Diptera) and evaluation of real-time PCR assays for the identification of Culicoides biting midges. Vet Parasitol. 2012;184:258–66.
Augot D, Mathieu B, Hadj-Henni L, Barriel V, Zapata Mena S, Smolis S, et al. Molecular phylogeny of 42 species of Culicoides (Diptera, Ceratopogonidae) from three continents. Parasite. 2017;24:23.
Kirkeby C, Dominiak P. Culicoides (Avaritia) gornostaevae Mirzaeva, 1984 (Diptera: Ceratopogonidae)—a possible vector species of the Obsoletus group new to the European fauna. Parasit Vectors. 2014;7:445.
Mathieu B, Delecolle JC, Garros C, Balenghien T, Setier-Rio ML, Candolfi E, et al. Simultaneous quantification of the relative abundance of species complex members: application to Culicoides obsoletus and Culicoides scoticus (Diptera: Ceratopogonidae), potential vectors of bluetongue virus. Vet Parasitol. 2011;182:297–306.
Pages N, Sarto IVM. Differentiation of Culicoides obsoletus and Culicoides scoticus (Diptera: Ceratopogonidae) based on mitochondrial cytochrome oxidase subunit I. J Med Entomol. 2005;42:1026–34.
Kiehl E, Walldorf V, Klimpel S, Al-Quraishy S, Mehlhorn H. The European vectors of Bluetongue virus: are there species complexes, single species or races in Culicoides obsoletus and C. pulicaris detectable by sequencing ITS-1, ITS-2 and 18S-rDNA? Parasitol Res. 2009;105:331–6.
Esselstyn JA, Evans BJ, Sedlock JL, Anwarali Khan FA, Heaney LR. Single-locus species delimitation: a test of the mixed Yule-coalescent model, with an empirical application to Philippine round-leaf bats. Proc R Entomol Soc. 2012;279:3678–86.
Talavera G, Dinca V, Vila R. Factors affecting species delimitations with the GMYC model: insights from a butterfly survey. Methods Ecol Evol. 2013;4:1101–10.
Dellicour S, Flot JF. Delimiting species-poor data sets using single molecular markers: a study of barcode gaps, haplowebs and GMYC. Syst Biol. 2015;64:900–8.
Rosenberg NA. The probability of topological concordance of gene trees and species trees. Theor Popul Biol. 2002;61:225–47.
Nichols R. Gene trees and species trees are not the same. Trends Ecol Evol. 2001;16:358–64.
Pentinsaari M, Vos R, Mutanen M. Algorithmic single-locus species delimitation: effects of sampling effort, variation and nonmonophyly in four methods and 1870 species of beetles. Mol Ecol Resour. 2017;17:393–404.
Fujita MK, Leache AD, Burbrink FT, McGuire JA, Moritz C. Coalescent-based species delimitation in an integrative taxonomy. Trends Ecol Evol. 2012;27:480–8.
Leache AD, Fujita MK. Bayesian species delimitation in West African forest geckos (Hemidactylus fasciatus). Proc R Entomol Soc. 2010;277:3071–7.
Boykin LM, Schutze MK, Krosch MN, Chomič A, Chapman TA, Englezou A, et al. Multi-gene phylogenetic analysis of south-east Asian pest members of the Bactrocera dorsalis species complex (Diptera: Tephritidae) does not support current taxonomy. J Appl Entomol. 2014;138:235–53.
Dinca V, Lukhtanov VA, Talavera G, Vila R. Unexpected layers of cryptic diversity in wood white Leptidea butterflies. Nat Commun. 2011;2:324.
Hsieh CH, Ko CC, Chung CH, Wang HY. Multilocus approach to clarify species status and the divergence history of the Bemisia tabaci (Hemiptera: Aleyrodidae) species complex. Mol Phylogenet Evol. 2014;76:172–80.
Song J-H, Ahn K-J. Species delimitation in the Aleochara fucicola species complex (Coleoptera: Staphylinidae: Aleocharinae) and its phylogenetic relationships. Zool Scr. 2014;43:629–40.
The authors are grateful to all European partners who assisted in Culicoides sampling and shipping. We acknowledge Bruno Mathieu for his friendly cooperation and for allowing us access to use some of his sequences as references for our study.
This study was partially funded by the EU grant H2020-727393 PALE-Blu and by the VectorNet project (OC/EFSA/AHAW/2013/02-FWC1) funded by the European Centre for Disease Prevention and Control (ECDC) and the European Food Safety Authority (EFSA). The contents of this publication are the sole responsibility of the authors and do not necessarily reflect the views of the European Commission, of the ECDC or of the EFSA.
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.
Information on adult female Culicoides sampling sites and results of specific assignation using cox1.
Primers used for PCRs and sequencing in this study.
Reference sequences used for specific assignation.
Information on all Culicoides sequenced and GenBank accession numbers.
Maximum likelihood phylogenetic tree using 16S rDNA. Values at the nodes represent bootstrap values greater than 800 for the most important nodes (1000 replicates).
Maximum likelihood phylogenetic tree using cox1. Values at the nodes represent bootstrap values greater than 800 for the most important nodes (1000 replicates).
Maximum likelihood phylogenetic tree using 28S rDNA. Values at the nodes represent bootstrap values greater than 800 for the most important nodes (1000 replicates).
About this article
Cite this article
Mignotte, A., Garros, C., Gardès, L. 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. Parasites Vectors 13, 265 (2020). https://doi.org/10.1186/s13071-020-04114-1