Skip to main content

A DNA barcode reference library of Croatian mosquitoes (Diptera: Culicidae): implications for identification and delimitation of species, with notes on the distribution of potential vector species

Abstract

Background

Mosquitoes pose a risk to human health worldwide, and correct species identification and detection of cryptic species are the most important keys for surveillance and control of mosquito vectors. In addition to traditional identification based on morphology, DNA barcoding has recently been widely used as a complementary tool for reliable identification of mosquito species. The main objective of this study was to create a reference DNA barcode library for the Croatian mosquito fauna, which should contribute to more accurate and faster identification of species, including cryptic species, and recognition of relevant vector species.

Methods

Sampling was carried out in three biogeographical regions of Croatia over six years (2017–2022). The mosquitoes were morphologically identified; molecular identification was based on the standard barcoding region of the mitochondrial COI gene and the nuclear ITS2 region, the latter to identify species within the Anopheles maculipennis complex. The BIN-RESL algorithm assigned the COI sequences to the corresponding BINs (Barcode Index Number clusters) in BOLD, i.e. to putative MOTUs (Molecular Operational Taxonomic Units). The bPTP and ASAP species delimitation methods were applied to the genus datasets in order to verify/confirm the assignment of specimens to specific MOTUs.

Results

A total of 405 mosquito specimens belonging to six genera and 30 morphospecies were collected and processed. Species delimitation methods assigned the samples to 31 (BIN-RESL), 30 (bPTP) and 28 (ASAP) MOTUs, with most delimited MOTUs matching the morphological identification. Some species of the genera Culex, Aedes and Anopheles were assigned to the same MOTUs, especially species that are difficult to distinguish morphologically and/or represent species complexes. In total, COI barcode sequences for 34 mosquito species and ITS2 sequences for three species of the genus Anopheles were added to the mosquito sequence database for Croatia, including one individual from the Intrudens Group, which represents a new record for the Croatian mosquito fauna.

Conclusion

We present the results of the first comprehensive study combining morphological and molecular identification of most mosquito species present in Croatia, including several invasive and vector species. With the exception of some closely related species, this study confirmed that DNA barcoding based on COI provides a reliable basis for the identification of mosquito species in Croatia.

Graphical Abstract

Background

Mosquitoes are one of the most studied groups of insects in the world, mainly because of their medical and veterinary role as vectors of pathogens causing millions of deaths per year.. Although most data on mosquito-borne diseases mainly refer to tropical and subtropical countries, quite a few cases of such infections have also recently been recorded in Europe [1,2,3,4]. Due to increasing trade and travel, invasive mosquito species are introduced and spread very easily [5,6,7], such as Aedes albopictus and Aedes aegypti, which are vectors of many viruses that cause various infectious diseases (e.g. dengue, chikungunya, yellow fever, Japanese encephalitis, West Nile and Zika viruses) [8,9,10]. In Croatia, well established populations of Aedes albopictus are present [11], while Aedes aegypti has not yet been recorded. In 2010, autochthonous cases of dengue fever caused by the dengue virus were registered for the first time in Croatia [12]. Also, several smaller dengue epidemics were continuously registered in Europe, with a number of autochthonous cases in France, Spain and Italy [13,14,15]. Although invasive species pose a serious medical problem, native mosquito species such as those of the genus Culex should not be ignored. This has been demonstrated by recent cases of West Nile virus (WNV) outbreaks in Romania, Spain, the Netherlands, Italy, Hungary, Germany, and Serbia [3, 16,17,18]. In Croatia, the first clinical cases of WNV infection were reported in 2012 [19], associated with Culex pipiens complex, and possibly some other species [20]. After 2012, WNV infections occurred continuously in Croatia [21]. The main vector of Usutu virus (USUV) infection is Culex pipiens, although the virus has also been found in several other species [22, 23]. Recent studies show that USUV has become endemic in northwestern Croatia [24]. Malaria used to be a prevalent disease in Europe, but currently only imported cases are reported; no deaths occurred following autochthonous infections in 2000–2019 [25]. However, numbers of imported malaria cases in Europe have increased, leading to the re-emergence of indigenous cases in Greece, Spain, Italy, and France [26,27,28,29]. Within the Anopheles maculipennis complex, several species are considered vectors of malaria parasites in Europe: Anopheles atroparvus, Anopheles labranchiae, Anopheles messeae and Anopheles sacharovi [30, 31].

The correct identification of mosquitoes is an important part of implementing effective vector management strategies. Conventional identification with dichotomous keys is essential, but has many shortcomings and is not always sufficient to identify mosquito specimens. Important features required for accurate morphological identification of mosquitoes often fall off or are damaged during sampling (scales, legs, wings), or these defining differences are only visible at a certain stage of development or are related to sex. In addition, mosquitoes often occur in complexes of closely related species, and morphological identification proved to be insufficient in most of these cases [32].

Molecular identification by DNA barcoding is an accurate method of species identification, independent of the developmental stage and condition of the specimens examined [33,34,35]. Many studies have confirmed that the standard mitochondrial COI barcoding region is a suitable marker for identification of mosquito species and recognition of cryptic species [36,37,38,39,40,41]. Nevertheless, for several genera and species complexes (Culex, Aedes and Anopheles), the use of additional molecular markers is necessary to increase identification accuracy and to distinguish between closely related species, forms and hybrids. For example, ITS2 (nuclear ribosomal internal transciber spacer 2), Ace2 gene and polymorphisms of various microsatellite loci (e.g. CQ11) are now routinely used [42,43,44,45,46,47]. However, despite major advances in molecular methods, accurate identification requires a multidisciplinary approach to taxonomy that includes morphological, molecular, distributional and ecological data [36, 48].

To date, 52 species of mosquitoes have been recorded in Croatia, two of which are invasive: Aedes albopictus and Aedes japonicus [11]. Traditionally, the majority of research has been based on morphological identification, but several recently published papers have used molecular methods to confirm the presence of certain species. By these, Cx. torrentium was proven for the first time to belong to the mosquito fauna of Croatia in 2018 [49], and the presence of several other species was confirmed for certain regions of Croatia [50, 51]. However, no systematic survey has been conducted to date that would provide an overview of the barcodes of the Croatian mosquito fauna, and there is also very little data for neighbour countries [52, 53]. The importance of such data lies in the potential presence and spread of new invasive species that are potential arbovirus vectors, such as Aedes koreicus, which is currently established in surrounding countries but has not yet been recorded in Croatia [54, 55]. Recent DNA barcoding studies in Europe [56,57,58,59,60,61,62,63,64] established national DNA barcode libraries for several countries, and resulted in new species findings and confirmation of unrecorded species in certain areas, as well as identification of cryptic taxonomic units [65]. Nevertheless, the number of barcode sequences of European mosquitoes with accurate country and species designation is still limited and accounts for less than 10% of all Culicidae records in BOLD (accessed December 15, 2023).

The aim of this study was to create a DNA barcode library for the Croatian mosquito fauna and to gain insight into the genetic diversity and geographical distribution of mosquito species through DNA barcoding. The data obtained will contribute to public databases (BOLD, NCBI GenBank) and help create a platform for easier, faster and more accurate identification of mosquitoes of the Croatian fauna. The results will also serve as a basis for projects on the surveillance of invasive and vector mosquitoes and control of mosquito-borne diseases in the studied region.

Methods

Study area

With its geographical location in Central Europe and the Mediterranean region, Croatia is one of the "hotspots" of European biodiversity. The high level of biodiversity in Croatia is a consequence of the diverse composition of habitats, climatic and hydrological characteristics and the complex geological history of the region.

According to Bertić et al. [66], Croatia was divided into three biogeographic regions for the purpose of this study:the Pannonian-Peripannonian region in the north and east, the central mountainous region in the middle and the Mediterranean region in the south (Fig. 1). From 2017 to 2022, 50 sites were sampled in the Pannonian-Peripannonian region, 38 in the central mountainous region and 52 in the Mediterranean region of Croatia (Fig. 1; Additional file: Table S1).

Fig. 1
figure 1

Sampling localities in ​​three biogeographic regions of Croatia: Pannonian-Peripannonian region—green dots, Mediterranean region—blue dots, and mountainous region—orange dots; high-altitude localities in the Pannonian-Peripannonian region are also marked with orange dots. Some localities may represent multiple sampling sites (details for each sample are available in the public BOLD dataset DS-CROCU2)

Mosquito collection and a priori morphological identification

Adult mosquitoes were caught with CDC traps baited with CO2, BG-Sentinel traps baited with BG-Lure and CO2 and human landing catch (HLC) method. Larvae were caught individually using a dipper. The mosquitoes were morphologically identified using the identification keys of Becker et al. [67]. All specimens were stored in 96% ethanol at -20° C. The vouchers of the barcoded specimens are kept at the Josip Juraj Strossmayer University of Osijek, Department of Biology, and at the Andrija Štampar Teaching Institute of Public Health in Zagreb.

DNA extraction and amplification

DNA was extracted from single legs of adult specimens and from entire larvae using the GenEluteTM Mammalian Genomic DNA Miniprep Kit (Sigma, Taufkirchen, Germany). The protocol for rodent tail preparation as provided with the kit was followed with slight modifications (incubation in proteinase K overnight; final DNA elution in 100 µl elution solution).

For all samples, the standard barcoding region of the mitochondrial COI gene [34] was amplified using the universal primers LCO1490 and HCO2198 [68]. For certain species complexes, the COI gene marker does not provide sufficient resolution, so other genomic regions were used. For the identification of species within the Anopheles maculipennis complex, the ITS2 region was amplified with the primers 5.8S (forward) and 28S (reverse) [69]. The amplification mixtures and PCR reaction conditions for COI and ITS2 are described in Bušić et al. [50]. The PCR products were enzymatically purified using the ExoI-rSAP system (NEB, Ipswich, MA, USA) according to the manufacturer's protocol and bidirectionally sequenced at Macrogen Inc. (Amsterdam, The Netherlands) using the amplification primers.

Data analysis

The sequences of the COI and ITS2 regions were verified and edited in Geneious 8.1.4. (https://www.geneious.com) and subsequently deposited in the NCBI GenBank and BOLD databases (GenBank accession numbers PP694665-PP694812, BOLD ID numbers in Additional file 1: Table S1; additional information available in public BOLD dataset DS-CROCU2). The percentage identity of the newly obtained COI and ITS2 sequences with the records in the GenBank database was checked with the BLAST tool using the Megablast algorithm (https://blast.ncbi.nlm.nih.gov/Blast.cgi; accessed December 15, 2023). The BOLD identification tool (http://www.boldsystems.org/index.php/IDS_OpenIdEngine; accessed December 15, 2023) was used to compare COI barcode sequences amplified from our samples with publicly available data in the BOLD database. Available published COI sequences of conspecific and congeneric mosquito specimens were selected from the BOLD database of public records (accessed December 15, 2023) and used for subsequent analyses.

COI and ITS2 sequences were analysed in datasets corresponding to specific genera. Multiple sequence alignments were performed with MAFFT version 7 using the "auto" strategy [70] (https://mafft.cbrc.jp/alignment/server/index.html; final alignments in Additional file 2: Dataset S1). Intraspecific and interspecific p-distances were calculated in MEGA 7.0.25 [71]. Maximum likelihood (ML) trees were generated on the PhyML 3.0 web server [72] (http://www.atgc-montpellier.fr/phyml/), with automatic model selection by SMS (Smart Model Selection; model determined by the AIC selection criterion) [73] and aLRT SH-like support [74].

Species delimitation methods (SDMs) bPTP [75] (https://species.h-its.org/ptp/) and ASAP [76] (https://bioinfo.mnhn.fr/abi/public/asap/asapweb.html) were applied to genus datasets to verify/confirm the assignment of specimens to specific Molecular Operational Taxonomic Units (MOTUs). As input for bPTP, the inferred ML tree was used, while for the ASAP method, the MAFFT alignment was used. In addition, the BIN-RESL algorithm (Barcode Index Number system in BOLD) [77] assigned the sequences to the corresponding BINs in BOLD. The results of the SDMs are presented in a combined ML tree that includes all newly barcoded specimens. The tree was graphically processed in FigTree v.1.4.3. (http://tree.bio.ed.ac.uk/software/figtree/).

Results

In this study, a total of 405 mosquito specimens were collected and processed; 181 specimens in the Pannonian-Peripannonian region, 122 specimens in the Mediterranean region and 102 specimens in the mountainous region of Croatia. Based on morphological identification, a total of 30 species previously recorded in Croatia were sampled, belonging to six genera (Aedes—16 species, Anopheles—5 species, Culex—5 species, Culiseta—2 species, Coquillettidia—1 species and Orthopodomyia—1 species) (Table 1). For some specimens, such as those of the Aedes annulipes/cantans/excrucians and Aedes cinereus/rossicus groups or some species complexes such as the Culex pipiens and Anopheles maculipennis complexes, it was not possible to determine morphologically to which species they belonged.

Table 1 Mean and maximum intraspecific distance (intra-SP) between individuals of the same BIN (species / MOTUs) and distance to nearest neighbour (NN) BIN in BOLD

The largest number of species (28) was sampled from the Pannonian-Peripannonian region, followed by the Mediterranean region (19) and the mountainous region (18). Ten species were recorded in all three regions (Additional file 1: Table S1). One individual from the Intrudens Group (CROCU102-21) was recorded in the Pannonian-Peripannonian region (Klokočevci – Fig. 1) and represents new data not only for this region, but also for the Croatian fauna in general.

The success rate of COI amplification and sequencing from both directions in relation to the number of samples processed was 72.7% (197 / 271 specimens), while the percentage of successful amplification reactions of ITS2 was 98.5% (132 / 134 specimens of the Anopheles maculipennis complex). The remaining specimens were excluded from further analyses due to failed amplification or low sequence quality. A minimum of one and a maximum of 21 COI sequences were obtained per species, with an average of six.

The BIN-RESL algorithm in BOLD assigned the specimens to 31 BINs (Table 1), with most delineated MOTUs matching the morphological identification. Several BINs were discordant, with two or more species placed together. BIN discordance was present in the genus Aedes and the genus Anopheles, mainly in the species that were difficult to distinguish morphologically. These are the BINs with samples identified as Ae. annulipes, Ae. cantans, Ae. riparius and Ae. excrucians (BOLD:AAB1098); Ae. rossicus and Ae. cinereus (BOLD:AAP8897); Ae. caspius and Ae. zammitii (BOLD:AAB7911); An. messeae and An. daciae (BOLD:ABY8239). Barcode gap analysis in BOLD confirmed the presence of a clear barcoding gap within the public BOLD dataset DS-CROCU2 (Fig. 2). The concordance with the morphological identification was 98.48%.

Fig. 2
figure 2

Frequency histogram of p-distances within and between investigated mosquito species (MOTUs) for the public BOLD dataset DS-CROCU2

In 13 species, the genetic distances of the species BIN to the nearest neighbour (NN) (i.e. non-specific specimen with the lowest interspecific distance in BOLD) were higher than the maximum intraspecific genetic distances, while other species had a lower value (Table 1). The mean intraspecific p-distance ranged from 0% to 1.68%, calculated for MOTUs delimited with BIN-RESL. The largest intraspecific divergence was observed in BIN BOLD:AAA7067 for the species Ae. vexans (7.22%), followed by BIN BOLD:AAB2483 for An. hyrcanus (5.32%) and BIN BOLD:AAA4751 for Cx. pipiens s.l. (5.02%). Only one species, Or. pulcripalpis (BOLD:AAW9535), showed no intraspecific variability in BOLD (Table 1).

The ASAP method delimited 28 MOTUs (Fig. 3). Some species of the genera Culex, Aedes and Anopheles were assigned to the same MOTUs, mostly again the species that are difficult to distinguish morphologically (Fig. 3). A very similar delimitation was achieved using the bPTP method, resulting in 30 MOTUs being assigned to our dataset. Using bPTP and BIN-RESL analysis, An. claviger s.s. was subdivided into two MOTUs. Also, Cx. modestus was subdivided into two MOTUs, but only using the bPTP delimitation method. The BIN-RESL method divided Cx. torrentium and Cx. pipiens s.l. into two different BINs, while ASAP and bPTP analyses classified them to the same MOTU.

Fig. 3
figure 3figure 3

ML tree based on COI DNA sequences of the specimens sequenced in this study. Numbers on branches denote bootstrap / aLRT support values (values lower than 70% are not shown). The results of species delimitation methods (bPTP, ASAP and BIN-RESL) are shown as vertical bars on the right

Most of the morphologically identified species formed well-supported monophyletic groups in the ML tree, which mostly corresponded to the assigned BINs, with support values between 81 and 100% bootstrap / > 0.9 aLRT (Fig. 3). Several species appeared to be paraphyletic, including Cx. pipiens s.l. in relation to Cx. torrentium and An. daciae and An. messeae in relation to An. maculipennis s.s.

Of the 132 specimens morphologically identified as Anopheles maculipennis complex, three species were confirmed using the ITS2 sequence: An. maculipennis s.s., An. messeae and An. daciae. An. maculipennis s.s. was the predominant species with 56.06% of the samples, followed by An. messeae with 36.36% and An. daciae with 7.57% of the samples.

Based on both morphological and molecular identification, a total of 34 species were identified in this study (Table 1).

Discussion

This work represents the first comprehensive study aimed at generating DNA barcodes for the establishment of a DNA barcode library for the Croatian mosquito fauna. Based on morphological characteristics, 30 mosquito species were identified, divided into six genera, and four more species were confirmed based on molecular data, which in total represents about 63.5% of the total number of mosquito species recorded in Croatia [11]. The validity and accuracy of the data obtained was examined using species delimitation algorithms, with various clustering methods additionally supporting the MOTUs obtained. This study shows that, in addition to morphological classification, molecular taxonomy can also be a suitable tool for the identification and delimitation of mosquito species in Croatia, as has recently been shown for many other countries [36,37,38, 57, 62, 78, 79]. The concordance of morphological and molecular identification in this study is high, which is also consistent with recent similar studies [38, 62, 80]. Morphological and molecular data discrepancies might be due to misidentifications, resulting in inaccurate species names for sequences in public databases [81]. This study found minimal differences in species delimitation methods, mostly for morphologically indistinguishable species. The effectiveness of barcoding was confirmed by the presence of a clear barcoding gap, which is necessary for the effective application of DNA barcodes to identify specimens and delimit species [62, 82]. According to Meyer & Paulay [83], a barcoding gap can be defined as the difference between the mean intraspecific sequence variability and the interspecific variability for congeneric COI sequences.

Based on DNA barcode analysis, most of the species morphologically identified in this study were categorised into established groups in ML tree analysis. A few individuals in certain taxa of the genera Aedes and Anopheles, such as Ae. annulipes/cantans/excrucians, Ae. rossicus/cinereus, Ae. caspius/zammitii, and An. daciae/messeae, were grouped within the same BINs, suggesting that the COI gene may not be informative enough for distinguishing among these species. This is also confirmed by the fact that these taxa occur in similar combinations in BOLD BINs. According to recent studies, these taxa are phylogenetically very closely related, although they have different morphological characteristics [38, 57, 62].

The Anopheles maculipennis complex consists of several species, seven of which have been recorded to belong to the Croatian fauna [11]. In this study three of them were identified, namely An. maculipennis s.s., An. messeae and An. daciae. In the COI sequence-based species delimitation, the Anopheles maculipennis complex was divided into two different MOTUs, An. maculipennis s. s. and An. messeae/daciae, by the results of the BIN-RESL algorithm. These two BINs are mutually nearest neighbours in BOLD, with only 1.33% distance between them, while in both BINs the highest intraspecific distance is about 4% (Table 1). The other two delimitation methods do not support this partitioning, as they group all An. maculipennis/messeae/daciae specimens into a single MOTU, with An. messeae/daciae specimens appearing as a paraphyletic group in the ML tree with respect to the An. maculipennis s.s. samples. This confirms the highly entangled situation in this species group which requires further investigation involving other molecular markers, similar as it was observed in Sedaghat et al. [84].

A rare species within the Anopheles maculipennis complex, Anopheles melanoon, previously recorded in the Croatian areas of southern Dalmatia and Istria [85], could not be confirmed in this study, possibly due to environmental modification of its distribution area. The occurrence of this species is closely linked to areas with horse and cow stables [86, 87], and such habitats have become quite rare. Other species within this complex previously detected in Croatia (An. atroparvus, An. sacharovi and An. labranchiae) [88, 89] were also not found in other recent surveys [50, 85]. Future studies should definitely focus on the detection of these species to confirm or exclude their current distribution in Croatia, as they are the most important malaria vectors in Europe [30].

The sibling species Cx. torrentium and Cx. pipiens s.l. can only be distinguished morphologically based on the characteristics of the male genitalia [67, 90]. The results of this study are consistent with the results of studies in Belgium [57] and show that these two species are separated as distinct MOTUs based on the BIN-RESL method, which was not the case in another recent study [62].

In contrast, the ASAP and bPTP methods did not show sufficient discrimination between species within genetically closely related species groups such as the Culex pipiens complex or the Anopheles maculipennis complex. For each of these two groups, the species are grouped into a single MOTU according to ASAP and bPTP.

The two sibling species Anopheles claviger s.s. (Meigen) and Anopheles petragnani Del Vecchio belong to the Anopheles claviger species complex [67], but An. petragnani has not yet been recorded in the Croatian fauna [11]. Within this species complex, COI has proven to be a sufficient tool to distinguish the two species [62]. Our two An. claviger s.s. specimens formed a strongly supported clade in the ML tree according to ASAP method, but were identified as separate MOTUs using bPTP and BIN-RESL methods. The specimen CROCU199-21 from a mountainous region was grouped in a BIN BOLD:AAM4220, where a wide European distribution with larger genetic distance (max. 3.56%) is observed. Another specimen, CROCU077-21 from the Mediterranean region, was grouped in BIN BOLD:AAE3979, with a narrower distribution (mainly Iran and Kosovo) and somewhat lower distances (max. 3.15%) (Table 1). A third BIN in BOLD with specimens from Tajikistan and China suggests that the species, An. claviger s.s., may be divided into several cryptic species. However, An. petragnani is placed in a separate BIN, including specimens from western Europe, some of them misidentified as An. claviger.

Culex modestus is considereda potential vector of WNV [91]. In this study, two individuals of Cx. modestus were detected in the Pannonian-Peripannonian region of Croatia near wetlands. This species is considered as the main vector of WNV in similar areas of southern France [92]. Although its dispersal ability is low, it could serve as an important enzootic and, given its ornithophilic and mammalophilic biting behaviour, bridging vector in natural/rural wetlands across Europe [67]. Therefore, in order to gain insights into the genetic structure of populations in this region, it should be studied in more detail in the future, using a larger number of individuals and localities.

Currently, there is no comprehensive database for all European mosquito species of the subgenus Ochlerotatus. According to Becker et al. [67], species groups associated with European species were classified according to region (Palaearctic, Fennoscandia, Germany and former USSR), which are also mentioned in our study (see below). The differences in the classification of species within each region will remain as long as there is no global analysis of the subgenus [67].

Within the Caspius Group there are six species, Aedes berlandi, Ae. caspius, Aedes dorsalis, Aedes mariae, Ae. pulcritarsis and Ae. zammitii [67], of which only Ae. berlandi has never been recorded in Croatia. According to all species delimitation algorithms used in this study, Ae. caspius and Ae. zammitii are grouped in the same MOTU, although they can be clearly distinguished morphologically. A similar case was observed in a recent study on the territory of Spain [80] for Ae. mariae and Ae. caspius which were grouped in the same MOTU. In our study the minimum interspecific distance between Ae. zammitii and Ae. caspius was 1.15%, similar to that in the aforementioned study [80], which is value lower than values typical between species [33]. Clearly, additional analyses of other loci that may better support species delimitation are required for a more accurate separation of species within this group.

Of the Annulipes Group, Ae. annulipes, Aedes behningi, Ae. cantans, Ae. excrucians, Aedes flavescens and Ae. riparius were recorded in Croatia [11, 67]. In general, distinction between species Ae. annulipes and Ae.cantans cannot always be accomplished with certainty based on morphological characters. Some sequences of the specimens morphologically identified as Ae. annulipes, Ae. cantans or Ae. excrucians matched the sequences in BIN: BOLD:AAB1098, which includes samples identified as Ae. annulipes/cantans or Ae. excrucians/cantans. One specimen, which had been assigned to the species Ae. behningi according to morphological characteristics, turned out to be Ae. cantans based on molecular analysis. In our ML tree, all specimens morphologically assigned to Ae. annulipes and Ae. cantans clustered into a single, well-supported clade, and all three species delimitation algorithms grouped the samples into a single MOTU. Consequently, COI proved to be an insufficient marker for distinguishing species within this species, as previously mentioned in other studies [38, 57].

In this study, the species Ae. cinereus and Ae. rossicus are grouped in the same MOTU according to all methods of species delimitation. Some samples morphologically identified as Ae. cinereus or Ae. rossicus were also confirmed based on COI sequences, while some could no be differentiated and had to be labelled as Ae. rossicus/cinereus (Fig. 3, Table 1). All individuals were collected in the Pannonian-Peripannonian region of Croatia. The taxonomic status of the species of the subgenus Aedes occurring in the Palaearctic is still unclear. Ae. cinereus, Ae. rossicus and Aedes esoensis were considered subspecies of the nominate form Ae. cinereus by Gutsevich et al. [93]. This view is not generally accepted, as Ae. cinereus and Ae. rossicus overlap strongly in Europe, the larvae often occur at the same breeding sites (which was also the case for our samples CROCU107-21, 108–21 and 109–21 from the same location and habitat—Fig. 3, Table 1), and transitional forms are unknown [67]. This is also supported by the fact that these two species are often found together, with samples from this study being assigned to the same BIN (BOLD:AAP8897).

A possible new species for the Croatian fauna, detected in this work, belongs to the Intrudens Group, which includes the species Ae. diantaeus, Ae. intrudens and Ae. pullatus [67]. Morphologically, the respective sample, CROCU102-21, was incorrectly determined as Aedes punctor, but after BIN-RESL analysis it was placed in BIN BOLD:AAF2904, which includes samples mainly identified as Ae. intrudens and Ae. diantaeus. This sequence matches 99.06% with the sequence of a specimen from Russia (KC855601) and 98.96% with a specimen from Sweden (JX040505), both of which were identified as Ae. intrudens. According to our photographs (available in the BOLD dataset DS-CROCU2) and the identification key according to Becker et al. [67], our specimen also corresponded to the description of the species Ae. intrudens. The eggs of this species overwinter, and the larvae are found from early spring onwards until the beginning of summer [67], which is consistent with our finding in April. In addition, the larvae of this species develop in temporary forest pools with dead leaves on the bottom [67], supporting our observations of its habitat. Considering that one single specimen was found, future studies should focus on the rediscovery and morphological and molecular confirmation of this species so that it can be reliably included in the list of Croatian mosquito species.

This six-year study covered 63.46% of the total mosquito fauna in Croatia. The distribution of species in Croatia by region determined in this study essentially corresponds to the results of previous studies [11, 50, 51]. It should be emphasized that some species included in the currently valid Croatian mosquito species list [11] were last recorded a long time ago, such as species of the Anopheles maculipennis complex: An. atroparvus, An. sacharovi and An. labranchiae [88, 89, 94]. Efforts have been made to eliminate malaria mosquitoes over the past centuries, and habitat changes and pollution may have contributed to reducing their population or causing extinction in this area [95,96,97]. However, as these species are the main vectors of malaria in Europe [30], future research should focus on confirming their presence (or absence) in the area. Species that were not recorded in this study but had been recorded in some other recent studies [98,99,100] are monocyclic, rare or uncommon [67], such as Aedes cataphylla, Aedes leucomelas, Ae. behningi, Ae. riparius, Culex martinii, Culiseta morsitans and Culiseta subochrea. In addition, some species are linked exclusively to certain habitats and periods of occurrence, and there was a high probability that they would not be recorded in this type of research. It is necessary to focus on these rare and unrecorded species in the future so that they can be sampled specifically and thus expand the database of Croatian mosquito barcodes in BOLD.

Conclusions

With barcoding sequences for 34 Culicidae species, the results of this study represent the basis for the establishment of a reference DNA barcode library for mosquitoes in Croatia. It has been proven that barcoding is an appropriate tool for the additional identification and delimitation of species of the Croatian mosquito fauna, even for closely related species such as Cx. pipiens s.l. and Cx. torrentium, while ITS2 proved a suitable marker to differentiate closely related species within the Anopheles maculipennis complex. The inability of COI to distinguish certain biotypes or other closely related species (such as Cx. pipiens biotype molestus; Ae. annulipes/cantans; An. messeae/ daciae) needs to be compensated for by additional molecular markers. In the future, studies should focus on confirming the species from the Croatian mosquito fauna list that were not included in this study. For the discovered potentially new species of the Intrudens Group, targeted investigations involving sampling of individuals in all life stages are needed to substantiate their inclusion in the list of Croatian mosquito fauna. Rapid and accurate identification is a crucial step in mosquito surveillance and control, so the data here are important not only for the assessment of biodiversity and of the geographical distribution of potential vector species in the studied area. The reference barcode sequences will contribute to future research on the mosquito fauna in Croatia, neighbouring countries and Europe in general, facilitate the identification or detection of potentially misidentified or cryptic species and provide a basis for invasive and vector species surveillance and monitoring projects.

Availability of data and materials

All relevant data are within the manuscript and its additional files. All information about barcoded specimens can be found in the public BOLD dataset DS-CROCU2. The vouchers of barcoded specimens are kept in Josip Juraj Strossmayer University of Osijek, Department of Biology, and in Andrija Štampar Teaching Institute of Public Health in Zagreb. DNAs of barcoded specimens are kept in Josip Juraj Strossmayer University of Osijek, Department of Biology.

References

  1. Lindh E, Argentini C, Remoli ME, Fortuna C, Faggioni G, Benedetti E, et al. The Italian 2017 outbreak chikungunya virus belongs to an emerging Aedes albopictus - adapted virus cluster introduced from the Indian subcontinent. Open Forum Infect Dis. 2018. https://doi.org/10.1093/ofid/ofy321.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Martinet J-P, Ferté H, Failloux A-B, Schaffner F, Depaquit J. Mosquitoes of north-western Europe as potential vectors of arboviruses: a review. Viruses. 2019;11:1059.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Bakonyi T, Haussig JM. West Nile virus keeps on moving up in Europe. Euro Surveill. 2020;25:46.

    Article  Google Scholar 

  4. Lazzarini L, Barzon L, Foglia F, Manfrin V, Pacenti M, Pavan G, et al. First autochthonous dengue outbreak in Italy, August 2020. Euro Surveill. 2020;25:2001606.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Kilpatrick AM. Globalization, land use and the invasion of West Nile virus. Science. 2011;334:323–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Vilà M, Hulme PE. Impact of biological invasions on ecosystem services. Berlin/Heidelberg, Germany: Springer International Publishing; 2017. https://doi.org/10.1007/978-3-319-45121-3. Accessed 18 Dec 2023.

  7. Dahmana H, Mediannikov O. Mosquito-borne diseases emergence/resurgence and how to effectively control it biologically. Pathogens. 2020;9:310.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Bhatt S, Gething PW, Brady OJ, Messina JP, Farlow AW, Moyes CL, et al. The global distribution and burden of dengue. Nature. 2013;496:504–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. de Wispelaere M, Desprès P, Choumet V. European Aedes albopictus and Culex pipiens are competent vectors for Japanese encephalitis virus. PLoS Negl Trop Dis. 2017;11:e0005294.

    Article  PubMed  PubMed Central  Google Scholar 

  10. European Centre for Disease Prevention and Control. Mosquito factsheets [Internet]. 2017. https://www.ecdc.europa.eu/en/disease-vectors/facts/mosquito-factsheets. Accessed 10 Jan 2024.

  11. Merdić E, Klobučar A, Žitko T, Sudarić Bogojević M, Vrućina I, Turić N, et al. Updated checklist of the mosquitoes (Diptera: Culicidae) of Croatia. J Vector Ecol. 2020;45:135–9.

    Article  PubMed  Google Scholar 

  12. Gjenero-Margan I, Aleraj B, Krajcar D, Lesnikar V, Klobučar A, Pem-Novosel I, et al. Autochthonous dengue fever in Croatia, August-September 2010. Euro Surveill. 2011;16:19805.

    Article  PubMed  Google Scholar 

  13. Riccò M, Peruzzi S, Balzarini F, Zaniboni A, Ranzieri S. Dengue fever in Italy: The “eternal return” of an emerging arboviral disease. Trop Med Infect Dis. 2022;7:10.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Brem J, Elankeswaran B, Erne D, Hedrich N, Lovey T, Marzetta V, et al. Dengue “homegrown” in Europe (2022 to 2023). New Microbes New Infect. 2024;56:101205.

    Article  PubMed  Google Scholar 

  15. Autochthonous vectorial transmission of dengue virus in mainland EU/EEA, 2010-present. European Centre for Disease Prevention and Control. 2018. https://www.ecdc.europa.eu/en/all-topics-z/dengue/surveillance-and-disease-data/autochthonous-transmission-dengue-virus-eueea. Accessed 28 Feb 2024.

  16. Busquets N, Laranjo-González M, Soler M, Nicolás O, Rivas R, Talavera S, et al. Detection of West Nile virus lineage 2 in north-eastern Spain (Catalonia). Transbound Emerg Dis. 2019;66:617–21.

    Article  PubMed  Google Scholar 

  17. Petrić D, Petrović T, Hrnjaković Cvjetković I, Zgomba M, Milošević V, Lazić G, et al. West Nile virus “circulation” in Vojvodina, Serbia: mosquito, bird, horse and human surveillance. Mol Cell Probes. 2017;31:28–36.

    Article  PubMed  Google Scholar 

  18. Rapid risk assessment: early large increase in West Nile virus infections reported in the EU/EEA and EU neighbouring countries 2018. European Centre for Disease Prevention and Control.https://www.ecdc.europa.eu/en/publications-data/rapid-risk-assessment-early-large-increase-west-nile-virus-infections-reported. Accessed 10 Jan 2024.

  19. Merdić E, Perić L, Pandak N, Kurolt IC, Turić N, Vignjević G, et al. West Nile virus outbreak in humans in Croatia, 2012. Coll Antropol. 2013;37:943–7.

    PubMed  Google Scholar 

  20. Merdić E 2013 Mosquitoes—vectors of West Nile virus in Croatia Rad HAZU 517 Med Sci 39 115 21.

    Google Scholar 

  21. Vilibic-Cavlek T, Janev-Holcer N, Bogdanic M, Ferenc T, Vujica Ferenc M, Krcmar S, et al. Current status of vector-borne diseases in Croatia: challenges and future prospects. Life. 2023;13:1856.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Nikolay B. A review of West Nile and Usutu virus co-circulation in Europe: how much do transmission cycles overlap? Trans R Soc Trop Med Hyg. 2015;109:609–18.

    Article  PubMed  Google Scholar 

  23. Vilibic-Cavlek T, Petrovic T, Savic V, Barbic L, Tabain I, Stevanovic V, et al. Epidemiology of Usutu virus: the European scenario. Pathogens. 2020;9:699.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Klobucar A, Savic V, Curman Posavec M, Petrinic S, Kuhar U, Toplak I, et al. Screening of mosquitoes for West Nile virus and Usutu virus in Croatia, 2015–2020. Trop Med Infect Dis. 2021;6:45.

    Article  PubMed  PubMed Central  Google Scholar 

  25. World Health Organization. World malaria report 2020: 20 years of global progress and challenges 2020. https://www.who.int/publications-detail-redirect/9789240015791. Accessed 10 Jan 2024.

  26. Kampen H, Maltezos E, Pagonaki M, Hunfeld K-P, Maier WA, Seitz HM. Individual cases of autochthonous malaria in Evros Province, northern Greece: serological aspects. Parasitol Res. 2002;88:261–6.

    Article  PubMed  Google Scholar 

  27. Doudier B, Bogreau H, DeVries A, Ponçon N, Stauffer WM, Fontenille D, et al. Possible autochthonous malaria from Marseille to Minneapolis. Emerg Infect Dis. 2007;13:1236–8.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Santa-Olalla Peralta P, Vazquez-Torres MC, Latorre-Fandos E, Mairal-Claver P, Cortina-Solano P, Puy-Azón A, et al. First autochthonous malaria case due to Plasmodium vivax since eradication, Spain, October 2010. Euro Surveill. 2010;15:19684.

    Article  CAS  PubMed  Google Scholar 

  29. Malaria - annual epidemiological report for 2021 2023. European Centre for Disease Prevention and Control. 2021 https://www.ecdc.europa.eu/en/publications-data/malaria-annual-epidemiological-report-2021. Accessed 10 Jan 2024.

  30. Sinka ME, Bangs MJ, Manguin S, Coetzee M, Mbogo CM, Hemingway J, et al. The dominant Anopheles vectors of human malaria in Africa, Europe and the Middle East: occurrence data, distribution maps and bionomic précis. Parasit Vectors. 2010;3:117.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Meibalan E, Marti M. Biology of malaria transmission. Cold Spring Harb Perspect Med. 2017;7:a025452.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Beebe NW. DNA barcoding mosquitoes: advice for potential prospectors. Parasitology. 2018;145:622–33.

    Article  PubMed  Google Scholar 

  33. Hebert PDN, Cywinska A, Ball SL, deWaard JR. Biological identifications through DNA barcodes. Proc R Soc B Biol Sci. 2003;270:313–21.

    Article  CAS  Google Scholar 

  34. Hebert PDN, Ratnasingham S, deWaard JR. Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species. Proc R Soc B Biol Sci. 2003;270:S96–9.

    Article  CAS  Google Scholar 

  35. Remigio EA, Hebert PDN. Testing the utility of partial COI sequences for phylogenetic estimates of gastropod relationships. Mol Phylogenet Evol. 2003;29:641–7.

    Article  CAS  PubMed  Google Scholar 

  36. Cywinska A, Hunter FF, Hebert PDN. Identifying Canadian mosquito species through DNA barcodes. Med Vet Entomol. 2006;20:413–24.

    Article  CAS  PubMed  Google Scholar 

  37. Kumar NP, Rajavel AR, Natarajan R, Jambulingam P. DNA barcodes can distinguish species of Indian mosquitoes (Diptera: Culicidae). J Med Entomol. 2007;44:1–7.

    Article  CAS  PubMed  Google Scholar 

  38. Hernández-Triana LM, Brugman VA, Nikolova NI, Ruiz-Arrondo I, Barrero E, Thorne L, et al. DNA barcoding of British mosquitoes (Diptera, Culicidae) to support species identification, discovery of cryptic genetic diversity and monitoring invasive species. ZooKeys. 2019;832:57–76.

    Article  Google Scholar 

  39. Adeniran AA, Hernández-Triana LM, Ortega-Morales AI, Garza-Hernández JA, de la Cruz-Ramos J, Chan-Chable RJ, et al. Identification of mosquitoes (Diptera: Culicidae) from Mexico State, Mexico using morphology and COI DNA barcoding. Acta Trop. 2021;213:105730.

    Article  CAS  PubMed  Google Scholar 

  40. Laurito M, Arias-Builes Dl Almirón WR. Improving the DNA barcode library of mosquito species with new identifications and discoveries in North-Central Argentina J Med Entomol. J Med Entomol. 2021https://doi.org/10.1093/jme/tjab160.

    Article  PubMed  Google Scholar 

  41. Viveros-Santos V, Hernández-Triana LM, Ibáñez-Bernal S, Ortega-Morales AI, Nikolova NI, Pairot P, et al. Integrated approaches for the identification of mosquitoes (Diptera: Culicidae) from the volcanoes of Central America physiographic subprovince of the state of Chiapas. Mexico Vector-Borne Zoonot Dis. 2022;22:120–37.

    Article  Google Scholar 

  42. Smith JL, Fonseca DM. Rapid assays for identification of members of the Culex (Culex) pipiens complex, their hybrids, and other sibling species (Diptera: Culicidae). Am J Trop Med Hyg. 2004;70:339–45.

    Article  CAS  PubMed  Google Scholar 

  43. Vinogradova EB Shaikevich EV 2007 Morphometric, physiological and molecular characteristics of underground populations of the urban mosquito Culex pipiens Linnaeus f. molestus Forskål (Diptera: Culicidae) from several areas of Russia Eur Mosq Bull 22 17 24.

    Google Scholar 

  44. Khoshdel-Nezamiha F, Vatandoost H, Oshaghi MA, Azari-Hamidian S, Mianroodi RA, Dabiri F, et al. Molecular characterization of mosquitoes (Diptera: Culicidae) in northwestern Iran by using rDNA-ITS2. Jpn J Infect Dis. 2016;69:319–22.

    Article  CAS  PubMed  Google Scholar 

  45. Shaikevich EV Vinogradova EB Bouattour A Gouveia de Almeida AP 2016 Genetic diversity of Culex pipiens mosquitoes in distinct populations from Europe: contribution of Cx. quinquefasciatus in Mediterranean populations Parasit Vectors 9 47.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Akiner MM Öztürk M Beriş FŞ Karacaoğlu Ç Şimşek FM Akgeyik AU 2022 Distribution and molecular differentiation of Culex pipiens complex species in the middle and Eastern Black sea regions of Turkey Turk J Zool 46 207 19.

    Article  CAS  Google Scholar 

  47. Vanderheyden A Smitz N Wolf De K Deblauwe I Dekoninck W Meganck K 2022 DNA Identification and diversity of the vector mosquitoes Culex pipiens s.s. and Culex torrentium in Belgium (Diptera: Culicidae) Diversity 14 486.

    Article  CAS  Google Scholar 

  48. Chan A, Chiang L-P, Hapuarachchi HC, Tan C-H, Pang S-C, Lee R, et al. DNA barcoding: complementing morphological identification of mosquito species in Singapore. Parasit Vectors. 2014;7:569.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Merdić E, Kujavec M, Kovačević M, Žulj M, Graovac N, Vignjević G, et al. Culex torrentium (Martini), a new species in Croatian mosquito fauna. Nat Croat. 2018;27:323–9.

    Article  Google Scholar 

  50. Bušić N, Kučinić M, Merdić E, Bruvo-Mađarić B. Diversity of mosquito fauna (Diptera, Culicidae) in higher-altitude regions of Croatia. J Vector Ecol. 2021;46:65–75.

    Article  PubMed  Google Scholar 

  51. Bušić N, Modrušan M, Vilc H, Merdić E. Contribution to the mosquito fauna (Diptera, Culicidae) of Lika, central Croatia, with special reference to invasive species. Nat Croat. 2021;30:231–42.

    Article  Google Scholar 

  52. Kavran M, Zgomba M, Weitzel T, Petric D, Manz C, Becker N. Distribution of Anopheles daciae and other Anopheles maculipennis complex species in Serbia. Parasitol Res. 2018;117:3277–87.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Török E, Kolcsár L-P, Keresztes L. New records and faunistic data of mosquitoes (Diptera, Culicidae) from Albania, Hungary, Macedonia, Montenegro, and Serbia. Turk J Zool. 2019;43:123–30.

    Article  Google Scholar 

  54. Kurucz K, Kiss V, Zana B, Schmieder V, Kepner A, Jakab F, et al. Emergence of Aedes koreicus (Diptera: Culicidae) in an urban area, Hungary, 2016. Parasitol Res. 2016;115:4687–9.

    Article  PubMed  Google Scholar 

  55. Kalan K, Šušnjar J, Ivović V, Buzan E. First record of Aedes koreicus (Diptera, Culicidae) in Slovenia. Parasitol Res. 2017;116:2355–8.

    Article  PubMed  Google Scholar 

  56. Engdahl C, Larsson P, Näslund J, Bravo M, Evander M, Lundström JO, et al. Identification of Swedish mosquitoes based on molecular barcoding of the COI gene and SNP analysis. Mol Ecol Resour. 2014;14:478–88.

    Article  CAS  PubMed  Google Scholar 

  57. Versteirt V, Nagy ZT, Roelants P, Denis L, Breman FC, Damiens D, et al. Identification of Belgian mosquito species (Diptera: Culicidae) by DNA barcoding. Mol Ecol Resour. 2015;15:449–57.

    Article  CAS  PubMed  Google Scholar 

  58. Talaga S, Leroy C, Guidez A, Dusfour I, Girod R, Dejean A, et al. DNA reference libraries of French Guianese mosquitoes for barcoding and metabarcoding. PLoS ONE. 2017;12:e0176993.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Muja-Bajraktari N, Zhushi-Etemi F, Dikolli-Velo E, Kadriaj P, Gunay F. The composition, diversity, and distribution of mosquito fauna (Diptera: Culicidae) in Kosovo. J Vector Ecol. 2019;44:94–104.

    Article  PubMed  Google Scholar 

  60. Deblauwe I, Ibáñez-Justicia A, Wolf De K, Smitz N, Schneider A, Stroo A. First detections of Culiseta longiareolata (Diptera: Culicidae) in Belgium and the Netherlands. J Med Entomol 2021;58:2524–2532.

    Article  CAS  PubMed  Google Scholar 

  61. Hohmeister N Werner D Kampen H 2021 The invasive Korean bush mosquito Aedes koreicus (Diptera: Culicidae) in Germany as of 2020 Parasit Vectors https://doi.org/10.1186/s13071-021-05077-7.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Madeira S, Duarte A, Boinas F, Costa OH. A DNA barcode reference library of Portuguese mosquitoes. Zoonoses Public Hlth. 2021;68:926–36.

    Article  CAS  Google Scholar 

  63. Bakran-Lebl K, Jerrentrup H, Daroglou E, Pfitzner WP, Fuehrer H-P, Allerberger F. First records of Aedes pulcritarsis (Rondani, 1872) (Diptera: Culicidae) in Austria. Parasitol Res. 2022;121:765–8.

    Article  PubMed  PubMed Central  Google Scholar 

  64. Kirik H, Tummeleht L, Kurina O. Rediscovering the mosquito fauna (Diptera: Culicidae) of Estonia: an annotated checklist with distribution maps and DNA evidence. Zootaxa. 2022;5094:261–87.

    Article  PubMed  Google Scholar 

  65. Lilja T, Troell K, Kirik H, Lindström A. A distinct group of north European Aedes vexans as determined by mitochondrial and nuclear markers. Med Vet Entomol. 2018;32:282–9.

    Article  CAS  PubMed  Google Scholar 

  66. Bertić I Lampek Pavčnik I. Republika Hrvatska-Prirodna obilježja, stanovništvo i geografske regije Satelitski Atlas Hrvatske Zagreb Naknada Ljevak & Gis Data. 2001

    Google Scholar 

  67. Becker N, Petrić D, Zgomba M, Boase C, Madon MB, Dahl C, et al. Mosquitoes: Identification, Ecology and Control. Cham: Springer International Publishing; 2020.

    Book  Google Scholar 

  68. Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994;3:294–9.

    CAS  PubMed  Google Scholar 

  69. Luca Di M Boccolini D Marinuccil M Romi R 2004 Intrapopulation polymorphism in Anopheles messeae (An. maculipennis complex) inferred by molecular analysis J Med Entomol 41 582 6.

    Article  PubMed  Google Scholar 

  70. Katoh K, Rozewicki J, Yamada KD. MAFFT online service: multiple sequence alignment, interactive sequence choice and visualization. Brief Bioinform. 2019;20:1160–6.

    Article  CAS  PubMed  Google Scholar 

  71. Kumar S Stecher G Tamura K 2016 MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets Mol Biol Evol 33 1870 4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Guindon S, Dufayard J-F, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010;59:307–21.

    Article  CAS  PubMed  Google Scholar 

  73. Lefort V, Longueville J-E, Gascuel O. SMS: smart model selection in PhyML. Mol Biol Evol. 2017;34:2422–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Anisimova M, Gascuel O. Approximate likelihood-ratio test for branches: a fast, accurate, and powerful alternative. Syst Biol. 2006;55:539–52.

    Article  PubMed  Google Scholar 

  75. Zhang J, Kapli P, Pavlidis P, Stamatakis A. A general species delimitation method with applications to phylogenetic placements. Bioinformatics. 2013;29:2869–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  76. Puillandre N, Brouillet S, Achaz G. ASAP: assemble species by automatic partitioning. Mol Ecol Resour. 2021;21:609–20.

    Article  PubMed  Google Scholar 

  77. Ratnasingham S, Hebert PDN. A DNA-based registry for all animal species: the barcode index number (BIN) system. PLoS ONE. 2013;8:e66213.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  78. Wang G, Li C, Guo X, Xing D, Dong Y, Wang Z, et al. Identifying the main mosquito species in China based on DNA barcoding. PLoS ONE. 2012;7:e47051.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  79. Ashfaq M, Hebert PDN, Mirza JH, Khan AM, Zafar Y, Mirza MS. Analyzing mosquito (Diptera: Culicidae) diversity in Pakistan by DNA barcoding. PLoS ONE. 2014;9:e97268.

    Article  PubMed  PubMed Central  Google Scholar 

  80. Delgado-Serra S, Viader M, Ruiz-Arrondo I, Miranda MÁ, Barceló C, Bueno-Marí R, et al. Molecular characterization of mosquito diversity in the Balearic Islands. J Med Entomol. 2021;58:608–15.

    Article  CAS  PubMed  Google Scholar 

  81. Krzywinski J, Besansky NJ. Molecular systematics of Anopheles: from subgenera to subpopulations. Annu Rev Entomol. 2003;48:111–39.

    Article  CAS  PubMed  Google Scholar 

  82. Meier R, Zhang G, Ali F. The use of mean instead of smallest interspecific distances exaggerates the size of the “barcoding gap” and leads to misidentification. Syst Biol. 2008;57:809–13.

    Article  PubMed  Google Scholar 

  83. Meyer CP, Paulay G. DNA barcoding: error rates based on comprehensive sampling. PLoS Biol. 2005;3:e422.

    Article  PubMed  PubMed Central  Google Scholar 

  84. Sedaghat MM, Linton Y-M, Nicolescu G, Smith L, Koliopoulos G, Zounos AK, et al. Morphological and molecular characterization of Anopheles (Anopheles) sacharovi Favre, a primary vector of malaria in the Middle East. Syst Entomol. 2003;28:241–56.

    Article  Google Scholar 

  85. Vignjević G. Molekularna identifikacija i distribucija vrsta komaraca Anopheles maculipennis kompleks u Hrvatskoj. Osijek, Croatia; 2014.

  86. Linton Y, Smith L, Koliopoulos G, Zounos AK, Samanidou-Voyadjoglou A, Patsoula E, et al. The Anopheles (Anopheles) maculipennis complex (Diptera: Culicidae) in Greece. J Nat Hist. 2007;41:2683–99.

    Article  Google Scholar 

  87. Nicolescu G, Linton Y-M, Vladimirescu A, Howard TM, Harbach RE. Mosquitoes of the Anopheles maculipennis group (Diptera: Culicidae) in Romania, with the discovery and formal recognition of a new species based on molecular and morphological evidence. Bull Entomol Res. 2004;94:525–35.

    Article  CAS  PubMed  Google Scholar 

  88. Adamović Ž. Anophelinae populations (Diptera, Culicidae) in Srem. Serbia Acta Vet. 1982;32:8.

    Google Scholar 

  89. Adamović Ž. Anophelinae mosquitoes in Podravina. Croatia Acta Entomol Jug. 1982;19:9.

    Google Scholar 

  90. Schaffner F, Angel G, Geoffroy B, Hervy J-P, Rhaiem A, Brunhes J. The Mosquitoes of Europe. An identification and training programme. CD-ROM. Paris: IRD editions. 2001.

  91. Ribeiro H Ramos HC Pires CA Capela RA 1988 An annotated checklist of the mosquitoes of continental Portugal (Diptera, Culicidae). Actas do III Congresso Iberico de Entomologia Granada Servicio de Publicaciones de la Universidad.

    Google Scholar 

  92. Balenghien T, Vazeille M, Grandadam M, Schaffner F, Zeller H, Reiter P, et al. Vector competence of some French Culex and Aedes mosquitoes for West Nile virus. Vector Borne Zoonotic Dis. 2008;8:589–95.

    Article  PubMed  Google Scholar 

  93. Gutsevich AV, Monchadskii AS, Shtakel’berg AA. Fauna of the USSR. Diptera. Mosquitoes of the family Culicidae. Leningrad, Leningrad Akad. Nauk SSSR Zool. Inst. N S No. 100, English translation: Israel Program for Scientific Translations: Jerusalem, 1974.

  94. Trausmiller O. Problem malarije u NR Hrvatskoj. Nar Zdrav. 1946;1:28–35.

    Google Scholar 

  95. Takken W Geene R Adam W Jetten TH Velden van der JA 2002 Distribution and dynamics of larval populations of Anopheles messeae and A. atroparvus in the delta of the rivers Rhine and Meuse, The Netherlands Ambio 31 212 8.

    Article  PubMed  Google Scholar 

  96. Bueno Marí R, Bernués Bañeres A, Jiménez PR. Biodiversity and larval habitat heterogeneity of mosquitoes (Diptera: Culicidae) in northern Spain. Acta Zool Bulg. 2014;66:547–54.

    Google Scholar 

  97. Kronefeld M, Werner D, Kampen H. PCR identification and distribution of Anopheles daciae (Diptera, Culicidae) in Germany. Parasitol Res. 2014;113:2079–86.

    Article  PubMed  Google Scholar 

  98. Vrućina I. Utjecaj prirodne i umjetne svjetlosti na ponašanje i aktivnost komaraca (Diptera, Culicidae) na urbanom i šumskom području istočne Hrvatske. Osijek, Croatia. https://www.bib.irb.hr:8443/index.php/891116. Accessed 11 Jan 2024. 2010.

  99. Merdić E, Sudarić Bogojević M, Boca I, Turić N. Determined and estimated mosquito (Diptera, Culicidae) fauna in the city of Osijek, Croatia, using dry-ice baited CDC traps. Period Biol. 2010;112:201–5.

    Google Scholar 

  100. Klobucar A, Benic N, Krajcar D, Kosanovic-Licina ML, Tesic V, Merdic E, et al. An overview of mosquitoes and emerging arboviral infections in the Zagreb area. Croatia J Infect Dev Ctries. 2016;10:1286–93.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We would like to thank Ivana Vrućina, mr.sc. and Željko Zahirović, mr.sc. (Josip Juraj Strossmayer University of Osijek) for assistance in the field while collecting samples and to Mario Dunić (Josip Juraj Strossmayer University of Osijek) for graphical assistance. We also thank two anonymous reviewers for constructive comments and suggestions.

Funding

This work has been supported by the Croatian Science Foundation under the project IP-2016–06-9988, DNA Barcoding of Diversity of Croatian Fauna (P.I. Prof. Dr. Sc. Mladen Kučinić) and by Josip Juraj Strossmayer University of Osijek project ZUP-2018–55 (P.I. Prof. Dr. Sc. Enrih Merdić) and internal project 3105–5 (P.I. Prof. Dr. Sc. Enrih Merdić).

Author information

Authors and Affiliations

Authors

Contributions

Conception and design of the work: NB, BBM, EM, MK; samples acquisition: NB, AK, NL, TŽ, GV, NT, MSB; analyses: NB, BBM; interpretation of data: NB, BBM, EM; manuscript draft: NB, BBM; manuscript revision: AK, EM; funding acquisition: MK, EM, NB.

Corresponding authors

Correspondence to Nataša Bušić or Branka Bruvo Mađarić.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

13071_2024_6291_MOESM1_ESM.xlsx

Supplementary Material 1: Table S1. All samples analysed in this study with their BOLD IDs, divided by regions of Croatia where they were found.

Supplementary Material 2. Dataset S1. Nucleotide alignments of COI and ITS2 datasets for all analysed specimens.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Bušić, N., Klobučar, A., Landeka, N. et al. A DNA barcode reference library of Croatian mosquitoes (Diptera: Culicidae): implications for identification and delimitation of species, with notes on the distribution of potential vector species. Parasites Vectors 17, 216 (2024). https://doi.org/10.1186/s13071-024-06291-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-024-06291-9

Keywords