- Open Access
MALDI-TOF mass spectrometry identification of mosquitoes collected in Vietnam
Parasites & Vectors volume 15, Article number: 39 (2022)
Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) is a tool that has revolutionised clinical microbiology and has recently been described as an innovative and effective approach to arthropod identification.
In this study, mosquitoes were captured in Vietnam using four different methods (human landing catch, CDC light traps, BG-Sentinel traps, animal-baited net traps). A total of 4215 mosquitoes were captured and morphologically identified as belonging to three genera: Aedes, Anopheles and Culex. We randomly selected 1253 mosquitoes, including 662 specimens of 14 Anopheles species, 200 specimens of two Aedes species and 391 morphologically unidentified Culex specimens, for molecular and MALDI-TOF MS analysis. The DNA from 98 mosquitoes (69 Anopheles specimens, 23 Culex specimens and six Aedes sp. specimens) was subjected to molecular analysis, either to confirm our morphological identification or the MALDI-TOF MS results, as well as to identify the Culex species that were morphologically identified at the genus level and to resolve the discrepancies between the morphological identification and the MALDI-TOF MS identification.
High-quality MS spectra were obtained for 1058 of the 1253 specimens (84%), including 192/200 for Aedes, 589/662 for Anopheles and 277/391 for Culex. The blind test showed that 986/997 (99%) of the specimens were correctly identified by MALDI-TOF MS, with log score values ranging from 1.708 to 2.843. Eleven specimens of Culex could not be identified based on morphological features, MALDI-TOF MS or molecular analysis.
This study enabled us to identify several species of mosquitoes from Vietnam using MALDI-TOF MS, and to enrich our database of MALDI-TOF MS reference spectra.
Mosquitoes are one of the most important insect species in the world and approximately 3500 species are currently recognised . They pose a major threat to human public health and economic development due to their ability to transmit and spread infectious pathogens, such as those causing malaria, dengue fever [dengue virus (DENV)], yellow fever [yellow fever virus (YFV)], Japanese encephalitis [Japanese encephalitis virus (JEV)], chikungunya [chikungunya virus (CHIKV)] and filariasis . Malaria, caused by Plasmodium malarial parasites (i.e. P. falciparum, P. vivax, P. malariae, P. ovale and P. knowlesi), is transmitted to humans by female Anopheles spp. mosquitoes [3,4,5]. The WHO estimated that there were 229 million cases of malaria in 2019, resulting in 409,000 deaths worldwide . Dengue infection, caused by any of the four DENV serotypes (DENV-1, -2, -3, and -4), is transmitted to humans via the bites of infected female Aedes aegypti or Aedes albopictus mosquitoes . Dengue is epidemic in 128 countries, with four billion people living in areas at risk of dengue infection . There are 100 million cases of dengue annually around the world, resulting in approximately 20,000 deaths .
In Vietnam, 255 mosquito species have been described, belonging to 42 subgenera and 21 genera . The majority fall into three genera, namely Aedes, Anopheles and Culex, and include vectors of agents of infectious pathogens . The number of species being discovered is steadily increasing, some of which are difficult to distinguish from others or are morphologically indistinguishable using traditional systematic methods . Molecular biology techniques, such as enzyme electrophoresis and PCR analyses, have been used to discriminate between homogeneous species [13, 14]. However, these techniques are limited due to being time-consuming and expensive and by requiring amplified gene sequences, which are not available in the GenBank databases for certain species [15, 16].
Recent studies show that matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) may be a rapid and accurate tool for the identification of arthropods, their blood meal sources and potentially associated microorganisms [17,18,19]. Compared to sequencing, MALDI-TOF MS is an easy-to-use, time- and cost-saving technique in entomology, especially if a machine that is generally purchased for microbiology is available [20, 21]. It is a technique based on the identification of proteins contained in a sample of interest. The mass/charge ratio of each protein is measured as it crosses an electromagnetic field after propulsion of the protein molecules by an ultraviolet laser. The mass/charge of the proteins thus measured generates a specific protein mass spectrum of the sample known as the “protein signature”. This spectrum is then compared to a reference database containing the spectra of species which have been formally identified morphologically and molecularly . In medical entomology, the choice of the body part of the arthropod and the method of preservation are parameters that can influence the quality of the spectrum . However, the choice of the body part to be used for MALDI-TOF MS analysis depends on the type of arthropod, such as the legs for mosquitoes and ticks [20, 22, 23], the head for bedbugs  and the cephalothorax for lice . This proteomic technique has overcome the limitations of traditional morphological and molecular identification methods because it does not require specific entomological knowledge and considerably reduces the time and cost involved [14, 21, 26].
The aim of the present study was to apply MALDI-TOF MS to identify some mosquito species collected in Vietnam in order to first create a reference database for rapid and accurate characterisation.
Mosquito collection and morphological identification
All mosquitoes were collected by an entomological team from the Institute of Malariology, Parasitology and Entomology in Quy Nhon, Vietnam (IMPE-QN), at the beginning of the rainy season (May) and the beginning of the dry season (January) between 2018 and 2020 in the Central Highlands of Vietnam. The collection area included eight provinces: Hoa Vang (16°03′N; 108°01′E) District, Da Nang Province; Nam Giang (15°65′N; 107°50′E) District, Quang Nam Province; Van Canh (13°37′N; 108°59′E) District, Binh Dinh Province; Dong Xuan (13°22′N; 109°02′E) District, Phu Yen Province; Ea Kar (12°49′N; 108°27′E) District, Dak Lak Province; Krong Pa (13°15′N; 108°45′E) District, Gia Lai Province; Khanh Vinh (12°16′N; 108°53′E) District, Khanh Hoa Province; and Ham Thuan Nam (11°09′N; 108°03′E) District, Binh Thuan Province (Fig. 1).
Collections were performed on eight consecutive nights/days using four different methods: human landing catch (HLC), CDC light trap (CDC-LT), animal-baited net trap and Biogents Sentinel mosquito trap (BG-Sentinel trap; BioGents AG, Regensburg, Germany). Collections took place in the forests, homes and plot huts and around local people’s homes in order to capture mosquitoes representing a greater number of genera and species. HLCs were used for night/day collection both indoors and outdoors. The mosquito collectors were local adults who volunteered to capture mosquitoes using HLCs and subsequently trained as mosquito collectors. They were regularly observed for signs of malaria or dengue and were treated if infected. A group of four people worked from 6:00 pm until 12:00 pm sampling Anopheles spp. mosquitoes, and from 5:30–7:30 am to 5:30–7:30 pm sampling Aedes spp. mosquitoes. Four CDC-LTs and two BG-Sentinel traps were also used for indoor and outdoor overnight captures. They were positioned outdoors with a minimum between-trap distance of 100 m. Indoor traps were placed near bed nets around human sleeping areas. All traps were suspended in the areas targeted from 6:00 pm until 6:00 am. The animal-baited net traps were positioned from 6:00 pm until 6:00 am.
Anopheles and Aedes mosquitoes were identified morphologically at the species level and Culex mosquitoes were identified at the genus level using Vietnamese mosquito identification keys [27,28,29,30]. Only female mosquitoes were used in this study, and captured specimens were stored in silica gel (Carl Roth GmbH, Karlsruhe, Germany) for the most part, or frozen at − 20 °C, before being sent to the Institut Hospitalo-Universitaire (IHU) Méditerranée Infection in Marseille, France for analysis. Samples were stored and subsequently analysed over a period of between 1 and 2.5 years.
Molecular identification of mosquitoes and phylogenetic analysis
Total DNA was extracted using a BioRobot EZ1 instrument (QIAGEN, Hilden, Germany) according to the manufacturer’s instructions. Cephalothoraxes from individual mosquitoes were placed in tubes containing 180 µl of G2 buffer and 20 µl of proteinase K and the samples incubated at 56 °C overnight. DNA extraction was performed from 200 µl of the incubation solution, eluted in 100 μl and stored at − 20 °C.
The DNA from 98 mosquitoes, including 43 mosquitoes the spectra of which were used to create our MALDI-TOF MS database and 55 mosquitoes selected for quality control of our MALDI-TOF MS identification, was subjected to molecular identification. The group of 55 mosquitoes selected for quality control of the MALDI-TOF MS identification method consisted of 34 mosquitoes whose identification by MALDI-TOF MS agreed with that by our morphological identification and 21 mosquitoes whose identification by MALDI-TOF MS showed discrepancies with that by our morphological identification. Standard PCR that specifically amplified a partial 720-bp cytochrome C oxidase I (COI) gene [forward primer COI_1 (LCO1490): 5′-GGT CAA ATC ATA AGA TAT TGG-3′; reverse primer (HC02198): 5′-TAA ACT TCA GGG TGA CCA AAA AAT CA-3′] was carried out and followed by sequencing . We also sequenced the 720-bp COI fragment of 34 specimens of different species of mosquitoes to confirm our MALDI-TOF MS identification using the same primers. Samples for which we did not successfully sequence the COI_1 gene were then subjected to another standard PCR and sequencing using the gene known as COI_2 [forward primer (CI-J-1632): 5′-TGATCAAAAAATTTATAAT-3′; reverse primer (CI-N-2191): 5′-GGGGTAAAAAAAATATAAACTTC-3′] amplifying a 560-bp partial sequence . We used the acetylcholinesterase 2 (Ace2) gene [direct primer (ACEpip): 5′-GGAAAACAAACGACGT ATGTACT-3′; reverse primer (B1246s): 5′-TGGAGCCTCCTCTTCACGGG-3′] amplifying a 610-bp fragment of Culex pipiens and a 274-bp fragment of Culex quinquefasciatus to distinguish the two species .
For specimens of Anopheles maculatus sensu lato (s.l.) that were not identified by MALDI-TOF MS, we performed standard PCR and sequencing targeting a partial 459-bp internal transcribed spacer 2 (ITS2) [forward primer (5.8F): 5′-TGTGAACTG CAGGACACATG-3′; reverse primer (28R): 5′-ATGCTTAAATTTAGGGGGTA-3′] .
Standard PCR and sequencing were performed as described in previous studies . Negative controls (single PCR mix and sterile water) and positive controls (DNA extracted from identified Anopheles gambiae and Ae. albopictus from our laboratory colony) were included in each PCR run. The sequences obtained were assembled and analysed using ChromasPro software (version 1.77) (Technelysium Pty Ltd., Tewantin, Australia) and submitted for analysis to the NCBI BLAST website (http://blast.ncbi.nlm.nih.gov). The phylogenetic tree was built using Mega 7.0 software (https://www.megasoftware.net) . Analyses were performed using the maximum likelihood method proposed by MEGA 7 and bootstrap analyses were performed using 500 replicates.
Preparation of mosquitoes for MALDI-TOF MS analysis
Mosquitoes with at least three legs were selected for the MALDI-TOF MS analysis. Each mosquito was treated individually by grinding the legs in a 1.5-ml Eppendorf tube containing 15 μl of 70% (v/v) formic acid (Sigma-Aldrich, St. Louis, MO, USA) and 15 μl of 50% (v/v) acetonitrile (Fluka, Buchs, Switzerland), with 1.0-mm-diameter glass beads (Sigma France, Lyon, France) using a tissue lyser machine (Qiagen); the optimal parameters had been previously established . The crushed legs were centrifuged, and a 1-μl aliquot of the supernatant of each sample was deposited in quadruplicate on a MALDI-TOF MS steel plate (Bruker Daltonics, Wissembourg, France) and covered after drying with a matrix solution composed of 1 μl of saturated alpha-cyano-4-hydroxycinnamic acid (Sigma France), 50% acetonitrile (v/v), 2.5% trifluoroacetic acid (v/v) (Sigma-Aldrich Co. Ltd., Gillingham Dorset, UK) and high-performance liquid chromatography-grade water [19, 37]. The target plate was then air-dried for a few minutes at room temperature before being introduced into the Microflex LT MALDI-TOF MS apparatus (Bruker Daltonics) for analysis. The quality of the matrix, the sample loading and the performance of the MALDI-TOF MS apparatus were controlled using the leg of an Ae. albopictus from our laboratory as a positive control.
MALDI-TOF MS parameters
The spectral profiles obtained from the mosquito legs were visualised using a Microflex LT MALDI-TOF mass spectrometer with FlexControl software (version 3.3; Bruker Daltonics). The spectra were acquired in a positive linear mode at a laser frequency of 50 Hz. The accelerating voltage was 20 kV, and the extraction delay time was 200 ns. Each spectrum corresponded to the ions obtained from the 240 laser shots performed in six regions at the same location and acquired automatically using the AutoXecute Flex Control software (version 2.4; Bruker Daltonics). The MS profiles were visualised using the FlexAnalysis v3.4 software, MALDI Biotyper Compass Explorer v4.1.70 (Bruker Daltonics), and ClinProTools v3.0 software (Bruker Daltonics) was used for data processing.
Spectral analysis and reference database creation
All spectra were visualised using the FlexAnalysis v.3.3 software to check their quality, reproducibility, and specificity. Spectra of poor quality, i.e. those which had low intensity (< 3000 AU), which were non-reproducible and which featured background noise, were excluded from the study. Finally, cluster analysis was performed to visualise the reproducibility and specificity of the spectra using the Main Spectrum Profile (MSP) dendrogram function of the MALDI Biotyper Compass Explorer v4.1.70 software. The dendrogram was performed using two to three spectra per species, excluding Anopheles barbirostris and Culex pallidothorax, for which there was only one spectrum for each species. The MS dendrogram is the result of comparing the MSPs produced by the MALDI Biotyper software and clustered according to the mass profile of the proteins (i.e. their mass signals and intensities), and this MS dendrogram illustrates how the samples are related to one another. The MS dendrograms helped us select representative reference spectra of the species . Between one and four reference spectra per species were then added to our in-house arthropod database of MS reference spectra . However, the reference spectra of Anopheles kochi and Anopheles jeyporiensis were not added to the database as the identification of An. kochi had not been molecularly confirmed and the spectrum of An. jeyporiensis was of poor quality.
Blind test for mosquito identification
The MS spectra from Ae. albopictus, Ae. aegypti and Culex spp. were blind tested against our pre-existing database containing spectra from Ae. albopictus, Ae. aegypti, Cx. quinquefasciatus and Culex sitiens species, the reference spectra of which were already available in our MALDI-TOF MS database. For specimens of Culex spp. that were not identified by this first blind test and the Anopheles species, a blind test was performed with the remaining spectra after creating a database of reference spectra of the different species confirmed by molecular biology. A total of 43 specimens, including between one and four of each species, were used to create the MALDI-TOF MS spectral database of Vietnamese mosquitoes. The level of reliability of species’ identification was determined using log score values (LSVs) of between 0 and 3 using the MALDI-TOF Biotyper software, which corresponded to a degree of similarity between the queried spectra and the reference spectrum in the database. The LSV threshold score set for correct identification of mosquito species was 1.70 in this study as this is the optimal threshold value applied for arthropod identification using MALDI-TOF MS [38, 39].
Mosquito collection and morphological identification
A total of 4215 individual mosquitoes were captured using different methods, including 1301 specimens (30.9%) by HLCs, 1335 (31.7%) by CDC-LTs, 891 (21.1%) by animal-baited net traps and 688 (16.3%) by BG-Sentinel traps. Mosquitoes collected by CDC-LTs, HLCs and animal-baited net traps were more abundant in quantity and number of species than those caught using the BG-Sentinel traps. The number of specimens collected per province was: Da Nang (n = 902), Quang Nam (n = 220), Binh Dinh (n = 1098), Gia Lai (n = 473), Dak Lak (n = 233), Phu Yen (n = 899), Khanh Hoa (n = 209) and Binh Thuan (n = 181) (Table 1; Fig. 1; Additional file 1: Table S1). Two Aedes species were identified morphologically: Ae. aegypti (n = 507; 12%) and Ae. albopictus (n = 640; 15.2%). Fourteen Anopheles species were identified based on morphological characteristics: Anopheles aconitus (n = 15; 0.4%), Anopheles annularis (n = 3; 0.1%), Anopheles barbirostris (n = 7; 0.2%), Anopheles dirus s.l. (n = 281; 6.7%), An. kochi (n = 8; 0.2%), An. maculatus s.l. (n = 1062; 25.2%), Anopheles minimus s.l. (n = 183; 4.2%), An. jeyporiensis (n = 1; 0.02%), Anopheles peditaeniatus (n = 81; 1.9%), Anopheles sinensis (n = 11; 0.2%), Anopheles splendidus (n = 36; 0.8%), Anopheles jamesii (n = 49; 1.2%), Anopheles vagus (n = 35; 0.8%) and Anopheles varuna (n = 9; 0.2%). The Culex mosquitoes, representing 1287 mosquitoes (30.7%), were not identified to the species level when specific morphological characters could not be identified with certitude due to being damaged during collection, transport or storage.
MALDI-TOF MS identification of the mosquitoes
Of the 4215 mosquitoes collected, 1253 (30%) mosquitoes, including 662 specimens of Anopheles spp., 391 specimens of Culex spp. and 200 specimens of Aedes spp., were randomly selected for MALDI-TOF MS analysis (Table 2). Overall, good-quality MS spectra were obtained from the legs of 1058 of the 1253 (84%) mosquitoes subjected to MALDI-TOF MS analysis. The percentage of good-quality spectra was highest for Aedes spp., with 96% (192/200) good-quality spectra, followed by Anopheles spp. (88.97%; 589/662) and Culex spp. (70.84%; 277/391). Visually, the comparison of MS profiles of legs of different species using FlexControl analysis revealed intra-species reproducibility and inter-species specificity (Additional file 2: Fig. S1). The MS protein profiles of between two and three specimens of each mosquito species used to generate an MSP dendrogram showed a clustering of specimens of the same species on the same branch (Fig. 2). The average discriminant peak intensity of 19 mosquito species, excluding those with low numbers, is presented in Additional file 3: Table S2.
Querying 192 good-quality spectra of Aedes spp. against our in-house MALDI-TOF MS database showed that 94 spectra were Ae. albopictus and 98 were Ae. aegypti, with LSVs ranging from 1.729 to 2.39 (Table 2). For Anopheles spp., we added 28 reference spectra from between one and four spectra of each species, confirmed by molecular biology analysis, to our MALDI-TOF MS database, which did not contain spectra of the concerned species. The remaining 561 spectra were then blind tested against the updated database, revealing that 96.8% (543/561) had been correctly identified, i.e. presented a match between our morphological identification and MALDI-TOF MS results. The LSVs of these Anopheles spp. specimens ranged from 1.708 to 2.821 (Table 2). For the remaining 18 specimens (3.2%) there were discrepancies between our morphological identification and the MALDI-TOF MS identification. Molecular biology analysis confirmed the MALDI-TOF MS identification of all species, with the exception of specimens of An. maculatus s.l. The dendrogram visually produced using the spectra of mosquitoes identified as An. maculatus s.l. shows that they form two distinct branches (Fig. 3).
Querying the 277 legs of Culex spp. against our database, which contained the spectra of 17 Culex species, showed that 28.2% (78/277) were identified, including 48 Cx. quinquefasciatus and 30 Cx. sitiens, with LSVs ranging from 1.725 to 2.401, and from 1.737 to 2.144, respectively. From the remaining 199/277 legs of Culex spp. which were unidentified, we randomly selected 15 of the 199 Culex specimens to create our database with reference spectra of four Culex tritaeniorhynchus, four Culex pseudovishnui, four Culex vishnui s.l., two Culex fuscocephala, and one Culex pallidothorax identified by molecular biology. The second blind test of the remaining 184 spectra against the updated database showed that 92.4% (170/184) were identified, with LSVs ranging from 1.739 to 2.843. Of the 173 Culex mosquitoes which were identified, 58 specimens were Cx. tritaeniorhynchus, 46 were Cx. pseudovishnui, 44 were Cx. vishnui s.l. and 25 were Cx. fuscocephala (including 3 specimens which showed discrepancies between the morphological and the MALDI-TOF MS results, which were subsequently confirmed by molecular analysis). Finally, 11 specimens with good spectra were not identified by MALDI-TOF MS and could not be identified molecularly because we could not sequence them using the three types of genes (Additional file 4: Fig. S2).
Molecular identification of the mosquitoes
Sequencing and BLAST of the COI_1 or COI_2 gene sequences from the Anopheles and Aedes species selected for the database revealed reliable and consistent species’ identification according to the morphological identification, with between 98.4 and 100% identity with GenBank sequences (Table 3). The COI_1 gene was used to identify species of 15 Culex spp. randomly selected to enter in our database. BLAST showed that four were between 99.3% and 100% identical to Cx. tritaeniorhynchus (MF179219; MF179221; MF179231; MF179232), four were between 97.9% and 98.2% identical to Cx. pseudovishnui (AB738092; LC054483), four were between 99.2% and 100% identical to Cx. vishnui (MF179240), two were 100% identical to Cx. fuscocephala (HQ398887) and one was 98% identical to Cx. pallidothorax (MF179281). Only one Cx. quinquefasciatus species could not be identified using COI_1 DNA barcoding; this species was then identified using the Ace2 gene, which was shown to share between 99.6% and 100% identity with GenBank (MN299021; MK575480) (Table 3).
Similarly, for specimens that showed discrepancies between our morphological identification and MALDI-TOF MS identification, i.e., An. aconitus identified as An. varuna by MALDI-TOF MS, An. dirus identified as An. sinensis, An. minimus identified as An. varuna, An. jamesi identified as An. maculatus, An. pediteniatus identified as An. sinensis, and Culex sp. identified as An. minimus, An. jamesii, and An. dirus. We found consistent identification between MALDI-TOF MS and molecular tools after querying the sequences against the GenBank database. However, for some specimens that were morphologically identified as An. maculatus s.l., the COI_1 sequences were 93.44% identical to An. maculatus (KT382822). We had no sequence for these specimens with the ITS gene generally used to differentiate between species of the An. maculatus complex (Table 3).
Finally, interrogation of the sequences from the samples selected for confirmation showed a perfect match between MALDI-TOF MS and molecular identification (Tables 2 and 3).
The COI_1, COI_2 or Ace2 sequences of some specimens were deposited in GenBank, and the names and accession numbers are listed in Additional file 5: Table S3. The phylogenetic tree constructed using the COI_1 sequences and their homologues available on GenBank shows that these 22 mosquito species formed distinct groups (Fig. 4) while correlated with morphological and MALDI-TOF MS identifications.
MALDI-TOF MS is an innovative tool for the rapid and accurate identification of arthropods . Here, we report the first MALDI-TOF MS reference spectra of a selection of protein extractions from Vietnamese mosquito legs. These species are widely distributed, medically important mosquitoes, and were most frequently captured during entomological surveillance and management programmes in the Central Highlands of Vietnam. Among them, An. minimus s.l., An. dirus s.l. and An. maculatus s.l. are considered to be the main vectors of human malaria in Vietnam [40,41,42] and other parts of Southeast Asia, such as Laos, Cambodia and Thailand [28, 43, 44]. Aedes albopictus and Ae. aegypti, are vectors of DENV, and recent epidemics of dengue fever in 2017 and 2019 resulted in 183,287 and 200,000 cases, with 32 and 50 deaths, respectively [45,46,47]. These two Aedes species are also important vectors of the CHIKV and Zika virus, both of which are present in the arbovirus burden of disease in Vietnam . Mosquitoes of the genus Culex are the main vectors responsible for the transmission of JEV throughout Vietnam . Despite the introduction of a JEV vaccine in 1997 and mandatory vaccination of children under the age of 5 years in Vietnam  over the past 10 years, the average annual number of JEV cases is between 1000 and 1200, with 20 to 50 deaths .
In terms of entomological surveillance, all mosquito species collected using the four methods in this study did not really have much of an impact on the quality of the MALDI-TOF MS profile, despite the duration of exposure to the mosquitoes according to temperature and humidity varying widely. After 12 h of use of the animal-baited net trap, the CDC-LT and the BG-sentinel trap, compared to the HLC method, the LSV median values remained at the optimal threshold score (i.e. 1.7). The stability of the MALDI-TOF MS results could likely be explained by the mosquitoes being handled carefully before being stored in dry tubes with silica gel.
Two species of Aedes and 14 species of Anopheles were identified based on morphological characteristics. In contrast, Culex specimens were only identified to the genus level based on morphological features due to the integrity of the mosquitoes, which did not allow for discrimination between morphologically identical sibling mosquito species and cryptic species groups. Although morphological identification is the “gold standard” for distinguishing mosquitoes based on external characteristics, it is unable to separate species belonging to the same complex. In the present study, almost all of the Culex species that were identified belong to species complexes that are not only difficult to distinguish morphologically  but also molecularly [23, 39, 52].
High-quality MS spectra were obtained for 84% of mosquito samples analysed by MALDI-TOF MS. The quality of the spectra obtained from Aedes spp. was higher than that obtained for MS spectra from Anopheles spp. and Culex spp. However, our percentage of good-quality spectra was lower than that found in a study of mosquitoes stored at − 20 °C . This is due to the difference in storage methods, as demonstrated by previous studies which have highlighted that storage at − 20 °C is the best method of preserving arthropods for MALDI-TOF MS analysis for periods not exceeding 6 months  while for MALDI-TOF MS analyses that are to be performed after periods exceeding 6 months, the best method of storage is − 80 °C . However, many laboratories do not have a − 80 °C freezer for mosquito storage, and in these cases − 20 °C freezing is the optimal method. In our study, conducted under field conditions, the percentage of good-quality spectra was higher than that of another study on mosquitoes which were also preserved in silica gel .
The MS spectra generated from the legs of the different mosquito species and the dendrogram show intra-species reproducibility and inter-species specificity or grouping of specimens of the same species on the same branch (Fig. 2; Additional file 2: Figure S1). This specificity of the MALDI-TOF MS spectra of mosquitoes has already been reported [20, 23, 53].
MALDI-TOF MS analysis allowed us to identify 22 species of mosquitoes, with LSV scores ranging from 1.708 to 2.843, including 18 species that were not in our MALDI-TOF MS database. Two species of Anopheles (An. kochi and An. jeyporiensis) were identified based on morphological features and 11 specimens of Culex sp. were not identified by MALDI-TOF MS because their spectra were not in our database. The latter is due to An. kochi and Culex not being molecularly identified in our study because we did not obtain DNA sequences from them, despite having good-quality spectra. An arthropod spectrum is added to our database only when the specimen is unambiguously identified by morphology and molecular biology analysis . We had only one specimen of An. jeyporiensis for which the MS spectrum was of poor quality [lower intensity: < 3000 AU; which was non-reproducible and which contained background noise). In addition to being able to identify closely related mosquito species, we found that MALDI-TOF MS can distinguish between subgroups of species of the same complex, which are impossible to distinguish morphologically, as was the case for the An. minimus complex (An. minimus and An. harrisoni formerly referred to as An. minimus A and C) [55, 56] and for the An. maculatus s.l. complex (Figs. 2, 3). Interestingly, the first two top hits for species classification were Cx. quinquefasciatus and Cx. pipiens, both ranging from 99 to 100% of sequence correspondence. These Cx. quinquefasciatus specimens were then classified by Ace_2 (Table 3). Preliminary studies had already reported that the MALDI-TOF MS technique is capable of identifying arthropods at the species, species subgroup or complex level [20, 57].
The morphological identification of all the specimens introduced in our database was confirmed by molecular biology using at least one gene. We found consistency between our morphological and molecular identifications, as our sequences had identities ranging from 96.6 to 100% with their GenBank homologues. For specimens that showed inconsistencies between morphological identification and MALDI-TOF MS identification, molecular biology analysis agreed with the MALDI-TOF MS identification. Morphological identification errors thus occurred, showing the limitations of this method. These errors were partly identified by MALDI-TOF MS, which we consider to be reliable for the identification of arthropods. Sequencing, which is the last technique we used to either confirm morphological and MALDI-TOF MS identification or to separate the two, did not allow us to identify 11 Culex specimens for which we had no morphological and MALDI-TOF MS identification of the species. This shows that all techniques have limitations and that molecular identification requires a reliable genomic sequence and the presence in GenBank of the homologous sequences of the specimen to be identified . Despite the difficulties and limitations of morphological identification, it remains the gold standard when it comes to identifying an unknown specimen, especially when the molecular method cannot identify it due to the absence of its sequence in GenBank or due to the lack of a sequence.
This study is the first to apply MALDI-TOF MS to the identification of field-caught mosquitoes in Vietnam, where the burden of disease transmitted by these vectors is high. It confirms that misidentification of mosquitoes by morphology can be prevented with MALDI-TOF MS, provided that the reference spectrum of that mosquito species is present in the MALDI-TOF MS database. The study allowed us to enrich our in-house MALDI-TOF MS database by adding reference spectra of new species. MALDI-TOF MS is a promising entomological surveillance tool that can significantly contribute to the control of mosquito-borne diseases in Vietnam.
Availability of data and materials
All relevant data are presented in the paper and its supporting ESM files. Our MALDI-TOF MS database is publicly available on our laboratory website and can be downloaded with the following DOI number: https://doi.org/10.35081/hwtr-5224.
Biogents Sentinel mosquito trap
Basic Local Alignment Search Tool
US Centers for Disease Control and Prevention miniature light traps
Log score values
- MALDI-TOF MS:
Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry
Xia H, Wang Y, Atoni E, Zhang B, Yuan Z. Mosquito-associated viruses in China. Virol Sin. 2018;33:5–20.
World Health Organization (WHO). Vector-borne diseases. 2017. https://www.who.int/news-room/fact-sheets/detail/vector-borne-diseases. Accessed 06 Jan 2021.
Nakazawa S, Marchand RP, Quang NT, Culleton R, Manh ND, Maeno Y. Anopheles dirus co-infection with human and monkey malaria parasites in Vietnam. Int J Parasitol. 2009;39:1533–7.
Van Long B, Allen G, Brauny M, Linh LTK, Pallerla SR, Huyen TTT, et al. Molecular surveillance and temporal monitoring of malaria parasites in focal Vietnamese provinces. Malar J. 2020;19:458. https://doi.org/10.1186/s12936-020-03561-6.
World Health Organization (WHO). World Malaria Report 2016. 2016. https://www.who.int/publications/i/item/9789241511711. Accessed 11 Apr 2021.
World Health Organization (WHO). Malaria. 2021. https://www.who.int/news-room/fact-sheets/detail/malaria. Accessed 18 Jan 2021.
Moncayo AC, Fernandez Z, Ortiz D, Diallo M, Sall A, Hartman S, et al. Dengue emergence and adaptation to peridomestic mosquitoes. Emerg Infect Dis. 2004;10:1790–6.
Brady OJ, Gething PW, Bhatt S, Messina JP, Brownstein JS, Hoen AG, et al. Refining the global spatial limits of dengue virus transmission by evidence-based consensus. PLoS Negl Trop Dis. 2012;6: e1760.
World Health Organization (WHO). Global strategy for dengue prevention and control, 2012–2020. https://www.who.int/denguecontrol/9789241504034/en/. Accessed 16 May 2021.
Van Dung N. Checklist of the known species of Culicidae found in Vietnam. Hoi Nghi Khoa Hoc Toan Quoc Ve Sinh Thai va Tai Nguyen Sinh Vat Lan. 2021, http://iebr.ac.vn/database/HNTQ6/504.pdf. Accessed 18 Jan 2021.
Tandina F, Doumbo O, Yaro AS, Traoré SF, Parola P, Robert V. Mosquitoes (Diptera: Culicidae) and mosquito-borne diseases in Mali, West Africa. Parasit Vectors. 2018;11:467.
Nadler SA, DE León GPP. Integrating molecular and morphological approaches for characterizing parasite cryptic species: implications for parasitology. Parasitology. 2011;138:1688–709.
Green CA, Gass RF, Munstermann LE, Baimai V. Population-genetic evidence for two species in Anopheles minimus in Thailand. Med Vet Entomol. 1990;4:25–34.
Lauri A, Mariani PO. Potentials and limitations of molecular diagnostic methods in food safety. Genes Nutr. 2009;4:1–12.
Hebert PDN, Cywinska A, Ball SL, De Waard JR. Biological identifications through DNA barcodes. Proc Biol Sci. 2003;270:313–21.
Vega-Rúa A, Pagès N, Fontaine A, Nuccio C, Hery L, Goindin D, et al. Improvement of mosquito identification by MALDI-TOF MS biotyping using protein signatures from two body parts. Parasit Vectors. 2018;11:574.
Laroche M, Jean-Michel Bérenger JM, Pascal Delaunay P, Didier Raoult D, Parola P, et al. Medical entomology: a reemerging field of research to better understand vector-borne infectious diseases. Clin Infect Dis. 2017;65:S30-8.
Sevestre J, Diarra AZ, Laroche M, Almeras L, Parola P. Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry: an emerging tool for studying the vectors of human infectious diseases. Future Microbiol. 2021;16:323–40.
Yssouf A, Socolovschi C, Flaudrops C, Sokhna CS, Raoult D, Parola P, et al. Matrix-assisted laser desorption ionization—time of flight mass spectrometry: an emerging tool for the rapid identification of mosquito vectors. PLoS ONE. 2013. https://doi.org/10.1371/journal.pone.0072380.
Diarra AZ, Laroche M, Berger F, Parola P. Use of MALDI-TOF MS for the identification of chad mosquitoes and the origin of their blood meal. Am J Trop Med Hyg. 2019;100:47–53.
Yssouf A, Almeras L, Raoult D, Parola P. Emerging tools for identification of arthropod vectors. Future Microbiol. 2016;11:549–66.
Diarra AZ, Almeras L, Laroche M, Doumbo O, Raoult D, Parola P, et al. Molecular and MALDI-TOF identification of ticks and tick-associated bacteria in Mali. PLoS Negl Trop Dis. 2017;11:e0005762.
Lawrence AL, Jana Batovska J, Webb CE, Šlapeta J, Parola P, Laroche M, et al. Accurate identification of Australian mosquitoes using protein profiling. Parasitology. 2019;146:462–71.
Benkacimi L, Gazelle G, El Hamzaoui B, Bérenger JM, Parola P, Laroche M. MALDI-TOF MS identification of Cimex lectularius and Cimex hemipterus bedbugs. Infect Genet Evol. 2020;85:104536. https://doi.org/10.1016/j.meegid.2020.104536.
Ouarti B, Laroche M, Righi S, Meguini MN, Benakhla A, Raoult D, et al. Development of MALDI-TOF mass spectrometry for the identification of lice isolated from farm animals. Parasite. 2020. https://doi.org/10.1051/parasite/2020026.
Favrot C. Polymerase chain reaction: advantages and drawbacks. In: 3. Congresso Latinoamericano de Dermatologia Veterinaria [Third Latin American Congress of Veterinary Dermatology], 26–27 Nov 2015; Buenos Aires. https://www.zora.uzh.ch/id/eprint/116536. Accessed 18 Jan 2021.
National Institute of Malariology, Parasitology and Entomology (NIMPE) Vietnam. Bang Dinh Loai muoi Anophelinae o Vietnam. 2006. http://nimpe.vn/TinChiTiet.aspx?id=33&cat=4. Accessed 21 Dec 2020.
Rattanarithikul R, Harrison BA, Harbach RE, Panthusiri P, Coleman RE. Illustrated keys to the mosquitoes of Thailand. IV. Anopheles. Southeast Asian J Trop Med Public Health. 2006;37:1–128.
Reuben R, Tewari SC, Hiriyan J, Akiyama J. Illustrated keys to species of Culex (Culex) associated with Japanese encephalitis in Southeast Asia (Diptera: Culicidae). Mosquito Systematics. 1994;26:75–96.
Rueda LM. Pictorial keys for the identification of mosquitoes (Diptera: Culicidae) associated with dengue virus transmission. Zootaxa. 2004;586:1–60.
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.
Kambhampati S. A phylogeny of cockroaches and related insects based on DNA sequence of mitochondrial ribosomal RNA genes. Proc Natl Acad Sci USA. 1995;92:017–20.
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.
Collins FH, Paskewitz SM. A review of the use of ribosomal DNA (rDNA) to differentiate among cryptic Anopheles species. Insect Mol Biol. 1996;5:1–9.
Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4.
Nebbak A, Willcox AC, Bitam I, Raoult D, Parola P, Almeras L. Standardization of sample homogenization for mosquito identification using an innovative proteomic tool based on protein profiling. Proteomics. 2016;16:3148–60.
Yssouf A, Parola P, Lindström A, Lilja T, Berenger JM, Raoult D, et al. Identification of European mosquito species by MALDI-TOF MS. Parasitol Res. 2014;113:2375–8.
Nabet C, Kone AK, Dia AK, Sylla M, Gautier M, Yattara M, et al. New assessment of Anopheles vector species identification using MALDI-TOF MS. Malar J. 2021;20:33.
Rakotonirina A, Pol M, Kainiu M, Barsac E, Tutagata J, Kilama S, et al. MALDI-TOF MS: optimization for future uses in entomological surveillance and identification of mosquitoes from New Caledonia. Parasit Vectors. 2020;13:359.
Do MC, Beebe NW, Nguyen TVV, Nguyen XT, Le NA, Cooper RD, et al. Vectors and malaria transmission in deforested, rural communities in north-central Vietnam. Malar J. 2010;9:259.
Garros C, Van Nguyen C, Trung HD, Bortel WV, Coosemans M, Manguin S. Distribution of Anopheles in Vietnam, with particular attention to malaria vectors of the Anopheles minimus complex. Malar J. 2008;7:11.
Ngo CT, Dubois G, Sinou V, Parzy D, Le HQ, Harbach RE, Manguin S. Diversity of Anopheles mosquitoes in Binh Phuoc and Dak Nong provinces of Vietnam and their relation to disease. Parasit Vectors. 2014;7:316.
Marasri N, Overgaard HJ, Sumarnrote A, Thanispong K, Corbel V, Chareonviriyaphap T. Abundance and distribution of Anopheles mosquitoes in a malaria endemic area along the Thai–Lao border. J Vector Ecol. 2017;42:325–34.
Vantaux A, Samreth R, Piv E, Siv S, Taylor WR, Ménard D, et al. Contribution to malaria transmission of symptomatic and asymptomatic parasite carriers in Cambodia. J Infect Dis. 2018;217:1561–8.
World Health Organization (WHO). Dengue situation update 533, 2 January 2018. 2018. https://reliefweb.int/report/viet-nam/dengue-situation-update-533-2-january-2018. Accessed 17 Jan 2021.
Quy N. Dengue fever outbreak kills 50 across Vietnam. VnExpress International. 31 Oct 2019. https://e.vnexpress.net/news/news/dengue-fever-outbreak-kills-50-across-vietnam-4005827.html. Accessed 17 Jan 2021.
Dang TT, Pham HM, Bui VH, Le VD. Whole genome sequencing and genetic variations in several dengue virus type 1 strains from unusual dengue epidemic of 2017 in Vietnam. Virol J. 2020;17:7.
Quyen NTH, Kien DTH, Rabaa M, Van Vinh Chau N, Wills B, Simmons CP, et al. Chikungunya and zika virus cases detected against a backdrop of endemic dengue transmission in Vietnam. Am J Trop Med Hyg. 2017;97:146–50.
Nguyen TY, Duffy MR, Nguyen MH, Nguyen TH, Fischer M, Hills SL. Surveillance for Japanese encephalitis in Vietnam, 1998–2007. Am J Trop Med Hyg. 2010;83:816–9.
Ngan D. Nang nong, canh giac cao voi benh viem nao Nhat Ban. HQ Online. 2020. https://haiquanonline.com.vn/nang-nong-canh-giac-cao-voi-benh-viem-nao-nhat-ban-127028-127028.html. Accessed 17 Jan 2021.
Shaikevich EV, Vinogradova EB, Bouattour A, de Almeida APG. Genetic diversity of Culex pipiens mosquitoes in distinct populations from Europe: contribution of Culex quinquefasciatus in Mediterranean populations. Parasit Vectors. 2016;9:47.
Batovska J, Blacket MJ, Brown K, Lynch SE. Molecular identification of mosquitoes (Diptera: Culicidae) in southeastern Australia. Ecol Evol. 2016;6:3001–11.
Tandina F, Niaré S, Laroche M, Doumbo OK, Raoult D, Parola P, et al. Using MALDI-TOF MS to identify mosquitoes collected in Mali and their blood meals. Parasitology. 2018;145:1170–82.
Nebbak A, El Hamzaoui B, Berenger J-M, Bitam I, Raoult D, Almeras L, et al. Comparative analysis of storage conditions and homogenization methods for tick and flea species for identification by MALDI-TOF MS. Med Vet Entomol. 2017;31:438–48.
Harbach RE, Garros C, Manh CN, Manguin S. Formal taxonomy of species C of the Anopheles minimus sibling species complex (Diptera: Culicidae). Zootaxa. 2007;1654:1. https://doi.org/10.11646/zootaxa.1654.1.3.
Harbach RE, Parkin E, Chen B, Butlin RK. Anopheles (Cellia) minimus theobald (Diptera: Culicidae): neotype designation, characterization, and systematics. Proc Entomol Soc Wash. 2006;108:198–209.
Niare S, Tandina F, Davous B, Raoult D, Parola P, Lionel Almeras L, et al. Accurate identification of Anopheles gambiae Giles trophic preferences by MALDI-TOF MS. Infect Genet Evol. 2018;63:410–9.
We are sincerely grateful to the staff of the Entomological Department of the Institute of Malariology, Parasitology and Entomology Quy Nhon (IMPE-QN), Vietnam, especially to the staff in the experimental entomology group, for their support with specimen collection and transportation. We would also thank the insectarium team of the Institut Hospitalo-Universitaire Méditerranée Infection, Marseille, France, for their assistance and the provision of the MALDI-TOF MS and molecular biology materials.
This study was supported by the Institut Hospitalo-Universitaire (IHU) Méditerranée Infection, the National Research Agency under the “Investissements d’avenir” programme, reference ANR-10-IAHU-03, the Région Provence Alpes Côte d’Azur and European ERDF PRIMI funding.
Ethics approval and consent to participate
This study was approved by the Research Ethics Committee of the Institute of Malariology, Parasitology and Entomology in Quy Nhon, on behalf of the Vietnam Ministry of Health (Approval Number 308/VSR-CT-2018). For ethical approval regarding human landing catches (HLC), we obtained consent from local adult volunteers. With regard to animals that were used for animal-baited net traps, such as domestic cows, we obtained permission from the owners of the cows.
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.
Additional file 1: Table S1.
Distribution of the adult mosquito species collected in different regions of Vietnam.
Additional file 2: Figure S1.
Comparison of MALDI-TOF MS spectra from the legs of 22 mosquito species collected in Vietnam between May 2018 and January 2020. MS spectra revealed intra-species reproducibility and inter-species specificity. Abbreviations: a.u., Arbitrary units; m/z, mass to charge ratio
Additional file 3: Table S2.
List of the top ten average mass peaks of the different mosquito species. The top ten average mass peaks by mosquito species are shown in bold. The average intensity of the mass peaks is indicated using a colour scale from green to red, indicating moderate and high average peak intensity, respectively. The list of mass peaks used to distinguish different mosquito species is based on the analysis of the ClinProTools genetic algorithm model. Abbreviations: Da, Daltons
Additional file 4: Figure S2.
Flow diagram of mosquito specimens which were included, analysed, and added to our in-house database using MALDI-TOF MS and molecular tools
Additional file 5: Table S3.
Accession numbers of sequences of different species of mosquitoes deposited in GenBank
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.
About this article
Cite this article
Huynh, L.N., Diarra, A.Z., Nguyen, H.S. et al. MALDI-TOF mass spectrometry identification of mosquitoes collected in Vietnam. Parasites Vectors 15, 39 (2022). https://doi.org/10.1186/s13071-022-05149-2
- MALDI-TOF MS
- Molecular identification