Introduction
The neotropical region is inventoried as a hotspot of mosquito (Diptera: Culicidae) species richness. French Guiana (hereafter FG) is a department of France situated in South America which harbors one of the highest relative species density of mosquitoes in the world [
1,
2]. Currently, 242 mosquito species, classified into 22 genera are known in this French territory [
3]. Some species are medically important because they are proven vectors of human pathogens. For example, several
Anopheles species are vectors of malaria parasitic infection in FG [
4,
5,
6,
7].
Yellow Fever (YFV), dengue (DENV), chikungunya (CHIKV) and Zika (ZIKV) viruses are transmitted by the day-biting
Aedes (
Stegomyia)
aegypti (Linnaeus) in urban area [
8,
9,
10].
Containing 788 species classified into 26 subgenera worldwide, the genus
Culex (abbreviated Cx.) is one of the largest groups of the Culicidae family [
11].
In FG, 104 species belong to the genus Culex [
3]
with a majority of species linked to two subgenus named Culex (abbreviated Cux.) and Melanoconion (abbreviated Mel.). In the Culex subgenus many species are vectors of arboviruses and parasites. Culex (
Cux.)
quinquefasciatus Say was identified as vector of the nematode
Wuchereria bancrofti, agent of lymphatic filariasis in FG [
12] and as primary vector in northeast Brazil [
13]. In a larger scale, this anthropophilic mosquito is also known to be vector of Western equine encephalitis virus (WEEV,
Alphavirus) and many viruses from the
Flavivirus genus, as Saint-Louis encephalitis
virus (SLEV) and West Nile virus (WNV); besides, is also vector of nematode
Dirofilaria immitis, that could infect both humans and dogs [
14,
15,
16,
17]. In South America,
Cx. (
Cux.)
declarator Dyar & Knab is a vector of SLEV and Bussuquara virus (BSQV,
Flavivirus genus) [
18].
Culex (
Cux.)
nigripalpus Theobald allowed the isolation of Cabassou virus (CABV,
Alphavirus genus) in FG [
19] and is a recognized vector for members of
Flavivirus genus with SLEV, WNV and ZIKV and for members of
Alphavirus genus with Venezuelan equine encephalitis virus (VEEV) and Eastern equine encephalitis virus (EEEV) in America [
20,
21].
Culex (
Cux.)
usquatus Dyar belonging to the Coronator Complex is not a recognized as vector, but the females are indistinguishable from
Cx. (
Cux.)
coronator Dyar & Knab, which is known as vector of WNV [
22] and may participate in the transmission of SLEV, VEEV and Mucambo Virus (MUCV,
Alphavirus genus) [
23,
24].
In the
Melanoconion subgenus,
Cx. (
Mel.)
portesi Senevet & Abonnenc and
Cx. (
Mel.)
spissipes Theobald are highly suspected to be natural vectors of Tonate virus (TONV) that belongs to the
Alphavirus genus in FG [
25,
26,
27], but are also found infected by
Bunyavirus (Caru, group Guama, Caraparu and group C) [
19].
Culex (
Mel.)
portesi is a recognized natural vector for MUCV in Trinidad [
28] and is a natural vector for the disease caused by El Huayo virus (a group C
Orthobunyavirus) in Peru [
29]. Ecological studies and the isolation of many arboviral species in
Cx. (
Mel.)
portesi rank it as high potential invertebrate host of numerous viruses [
19].
Culex (
Mel.)
adamesi Sirivanakarn & Galindo,
Cx. (
Mel.)
dunni Dyar, and
Cx. (
Mel.)
pedroi Sirivanakarn & Belkin are natural enzootic vectors of VEEV in Columbia, Peru and Panama [
30,
31,
32].
Culex (
Mel.)
dunni was also found infected with Pacora virus (PCAV,
Bunyavirus-like) in Panama [
33]. A last example is
Cx. (
Mel.)
idottus Dyar which is a caiman-biting mosquito involved in the transmission of
Hepatozoon caimani and suggested as potential vector of WNV [
34].
With global warming, deforestation and urbanization, the distribution areas of mosquitoes change, colonizing new territories. Such changes could promote the emergence of new mosquito vectors due to their exposure to disease agents and hosts, inducing potential outbreaks. To manage and to prevent epidemic emergence, the monitoring and accurate identification of Culex mosquitoes at the species level remains essential. Presently, morphological and molecular methods are the main strategies available for identification and classification of Culex species.
Morphological identification at species level of
Culex mosquitoes requires the careful dissection and mounting of male genitalia by an entomological expert because of the presence of slight morphologic characters that allows the identification of sibling species of the genera [
35]. Females identification is even harder because morphological characters may be either polymorphic or isomorphic among distinct species [
36]. This approach is labor-intensive and time-consuming and, therefore, might not be adapted for routine identification of
Culex females. Furthermore, field-caught specimens are rarely in perfect condition, which preclude any reliable morphological identification. For damaged mosquitoes or for immature stages, molecular tools can be an interesting alternative, especially to distinguish morphologically close species [
37]. Barcoding and metabarcoding markers which combine DNA barcoding with high-throughput sequencing are sometimes used for delimitation and identification of species. However, species-level identification depends heavily on the choice of marker and the selected primer pair, often with a trade-off between successful species amplification and taxonomic resolution. Variation of a partial sequence of the cytochrome
c oxidase subunit I (
COI) gene is often used for the identification of mosquito species and contribute to discover cryptic diversity [
38]. If on the one side, the COI barcode can be successfully used for delimiting and identifying numerous mosquito species, on the other hand recent studies on the
Culex subgenus identification with this fragment indicate a poor resolution in separating species among complexes [
36]. To improve species identification, multiple marker/primer pairs are often recommended and remains an expensive method for screening of large samples [
39].
In the last decade, an innovative and cheaper proteomic tool named MALDI-TOF MS profiling, emerged for rapid mosquito species identification [
40]. The principle is based on the matching of species-specific protein signatures of a specimen with reference spectra database. This promising tool has been reported to distinguish cryptic species of
Anopheles genus with high efficiency [
41]. Already used successfully for different taxonomic groups of insects, including culicids [
42] and phlebotomids [
43], standardization of protocols and optimized procedures could enable sharing all MS spectra reference and lead to the creation of an international MS database [
44].
Thus, the aim of the present study was to assess whether this proteomic tool could distinct species from the neotropical Culex genus. In this way, field-caught mosquitoes from FG of the Culex genus were selected, classified by morphological criteria and confirmed by molecular barcoding. The legs and thoraxes of the selected Culex specimens were then submitted independently to MALDI-TOF MS to assess intra-species reproducibility and inter-species specificity of MS spectra. The establishment of such innovative tool for Culex mosquito identification should improve studies on this genus which are often hampered by the complexity of current methods for accurate classification.
Discussion
MALDI-TOF MS is a method routinely used in microbiological diagnostic laboratory for identification of bacteria and archaea [
59], but the relevance of the technique remains largely underestimated for identification of multicellular organisms. The success of the MALDI-TOF MS for arthropod identification in this last decade, highlights this emerging tool as a relevant alternative for mosquitoes species identification [
43,
60,
61,
62]. As many mosquitoes species from the
Culex genus have been proven to be vectors of pathogenic agents of human and veterinary importance [
63,
64,
65], an accurate classification of the specimens is primordial. The present work demonstrated that MALDI-TOF MS can be used as an alternative to current methods for
Culex female identification at the adult stage, using legs and thorax.
Prior MS submission, a relevant classification of the
Culex mosquitoes per species was compulsory. Firstly, the morphological identification was performed allowing to classify mosquitoes from only five distinct species, two from the
Culex subgenus (
Cx. quinquefasciatus and
Cx. usquatus) and three from the
Melanoconion subgenus (
Cx. portesi, Cx. pedroi and
Cx. spissipes). The accuracy of this morphological classification was confirmed by molecular DNA barcoding of the COI gene. Morphological method failed to identify at the species level 29% (n=60) of
Culex mosquitoes. The limitation of classification of mosquitoes by morphological approach was also reported for others genera [
36,
66,
67]. For instance, in the
Anopheles genus
, the members of the Gambiae Complex and Funestus Group could not be distinguished uniquely by examining external features of their anatomy [
68,
69]. Moreover, the skill and regular training of entomological experts is also another factor playing a role in the capacity of correct specimen identification [
70]. The genus
Culex is recognized as a highly diverse group of mosquitoes for which delimitation and identification of species is particularly difficult. Complementary methods are therefore required for improving species identification.
DNA barcoding using the mitochondrial COI gene is widely used for species identification in molecular taxonomy [
71]. A first COI barcoding database, available in the BOLD system, established based on morphological classification of male
Culex specimens of FG [
3] was used to validate morphological identification and to classify the unidentified females of this study. The COI gene has been recorded as an effective and accessible DNA barcode that provides a high performance in delimiting species within the subgenus
Melanoconion of
Culex species, a subgenus displaying the highest species diversity in tropical regions [
66]. However, others studies also applied COI barcode fragment for identification of species within the subgenus
Culex and concluded that this barcode do not contain enough information to distinguish species within this subgenus [
36]. Based on COI sequences from Brazil and Argentina, Laurito et al. [
36] highlighted that
Cx. declarator cannot be differentiated from
Cx. bidens Dyar and
Cx. tatoi Casal & García. However, these last two species have never been detected in FG [
3], and only
Cx. declarator was considered here. Similarly, COI sequences of
Cx. usquatus can be confused with
Cx. camposi Dyar,
Cx. coronator,
Cx. maxi Dyar and
Cx. surinamensis [
36]. In FG, only
Cx. coronator,
Cx. surinamensis and
Cx. usquatus have already been reported. Nevertheless, a recent review of
Culex species from FG revealed that specimens identified as
Cx. coronator in historical literature were closer to the lectotype of
Cx. usquatus described from Suriname [
3].
Culex usquatus is morphologically distinguishable from
Cx. surinamensis, thus confirming that COI sequences classified as
Cx. usquatus could be attributed to this species in this study.
To improve the identification, molecular tools often develop new algorithms of analysis or add others molecular markers. For
Anopheles species, COI is not enough to discriminate species of the Gambiae Complex. The addition of ITS2 analysis increases greatly the potential of identification of these species [
72]. Indeed, the understanding lack on COI sequences analysis highlight the need to explore complementary or innovative tools.
The proteomic MALDI-TOF MS approach is increasingly used for rapid arthropod identification [
42,
43,
60,
61,
73]. Indeed, to improve the intra-species reproducibility of MS spectra and to share reference protein profiles, protocols were standardized [
56,
74]. In order to develop an accessible mosquito spectra database of
Culex species, the present study employed MALDI-TOF MS double spectra biotyping strategy. Here, for each specimen tested and for DB creation, legs and thorax from the same specimen were independently submitted to MS analysis, whenever possible. Specific MS spectra from the 13
Culex species per body part were obtained. The cluster analysis on MSP dendrogram confirmed this species specificity. Besides, similar ordination of the specimens for both body parts with the COI phylogenetic tree, revealed a relationship between molecular taxonomic classification and MS spectra profiles. However, as MS spectra from closely-related species were grouped in the same branch of MS dendrogram, a mismatching could occur. To prevent misidentification, the legs and thorax were independently query against the DB. This double query allowed to corroborate identification of each specimen using both body parts, enhancing the identification confidence [
55]. In this study, 100% of the samples were correctly identified at the species level on both body parts and the rate of relevant identification (LSVs>1.8) was about 95% for each ones. Interestingly, a higher reproducibility of MS spectra per species was obtained for thorax compared to legs. The lower reproducibility of legs could be attributed to the number of legs available for each specimen. Effectively, legs are breakable and during the catching or storing period some could be lost, which could induced heterogeneity of MS profile intensities among specimens from the same species [
75]. In the cases of all the legs were lost, the specimen could be always identified using thorax compartment. Then, the thorax appears as the more appropriate compartment for mosquito identification, followed by legs accordingly to previous work [
44]. Furthermore, these two body parts does not prevent screening for viruses, parasites, or source of blood feeding which can be researched in the head or abdomen by molecular [
76,
77] or MALDI-TOF MS methods [
78,
79]. These complementary data are primordial in the frame of surveillance programs.
Quality of MS spectra can be altered by various factors such as sample homogenization, quantity of mix buffer, engorgement status, storage conditions or even geographical origin of the collection [
42,
55,
75]. The storage is too critical for field caught mosquitoes. The best method for long storage of arthropods for MALDI-TOF MS analysis is frozen or to maintain them at RT with silica gel [
80], when immediate freezing is not possible. In this study, all the specimens were field-derived, dried and stored frozen at -20°C in tubes or plates from one to three years. The elevate intra-species reproducibility and inter-species specificity of MS spectra for both body parts, suggested that specimens were appropriately preserved. Among the 169 mosquitoes tested, only two
Cx. usquatus failed, on both body parts, to reach the relevant threshold value (LSV>1.8). These low LSVs could be attributed to improper storing of the mosquitoes. It is possible that protein degradation occurred in these specimens conducting to lower quality of MS spectra, as reported in previous studies [
58,
81,
82].
It is noteworthy that all the Culex species selected in this study were correctly classified using MALDI-TOF MS profiling, with concordant results between thorax and legs, and in agreement with molecular identification. Three species from the Melanoconion subgenus, Cx. idottus, Cx. phlogistus and Cx. rabanicolus, presenting isomorphic traits among females and which could not be accurately distinguished morphologically, were unambiguously identified by this proteomic tool.
The present spectral database, created using thoraxes and legs from females of 13
Culex species from FG, represents the first attempt to create a MALDI-TOF MS database for identification of neotropical
Culex species. However, the database encompass only 12.5% (n=13/104) of the total number of
Culex species known in FG and solely two out of the eight subgenera from the
Culex genus of the mosquitoes actually inventoried in the territory [
3]. The widening of this MALDI-TOF MS spectra reference DB with missing
Culex species, become compulsory for application of this tool in the monitoring of
Culex genus vectors in FG. An upgrading of this DB including others mosquito species occurring in FG (n=242) will improve its usage and may be helpful to discriminate cryptic or close-related species. MS spectra for the identification of eight distinct
Anopheles species from FG, among which four are malaria vectors, are already available [
54].
Several works reported the efficiency of the MALDI-TOF MS for identification of mosquitoes at immature stages [
83,
84] and demonstrated it interest for management of Culicidae larval habitats [
85]. As larva can be confusing to distinguish morphologically, it could be interesting to assess MALDI-TOF MS for identification of immature stages. MALDI-TOF MS appears then highly promising for discriminating mosquito fauna [
41,
56,
82].
The present work demonstrated that MALDI-TOF MS could be an alternative to current methods for identification of
Culex female, using dissected legs and thorax. The double body part protein signature of each specimen improved the identification quality. Raw MS spectra from legs and thorax of these 13
Culex species included in the DB are freely available (
Supplementary file S1). The sharing of reference MS spectra is essential in the frame of the creation of an international DB. Except for the expensive cost of the MALDI-TOF MS instrument, this approach is highly competitive economically compared to current molecular methods. It does not require particular skills and could be used for “live” monitoring of vectors due to its short time of handling and obtaining results. This approach appears suitable for the identification of mosquitoes from the
Culex genus, for which morphological and molecular methods remains time consuming and with a potential risk of misidentification. In the near future, the widening of this reference MS spectra DB with specimens from others
Culex species, at both adult and immature stages, could improve knowledge of these species for adapted surveillance and control measures to reduce the risks of pathogen transmission not exclusively to French Guiana.