Geometric morphometrics versus DNA barcoding for the identification of malaria vectors Anopheles dirus and An. baimaii in the Thai-Cambodia border (original) (raw)

Introduction

Malaria is a mosquito-borne infectious disease caused by unicellular protozoan parasites belonging to the genus Plasmodium. It is a significant public health problem in many tropical and subtropical countries1. In 2020, the World Health Organization reported an estimated 241 million cases of malaria, leading to 627,000 deaths worldwide[2](/articles/s41598-022-17646-6#ref-CR2 "World Health Organization. WHO malaria report. 2021; https://www.who.int/news-room/fact-sheets/detail/malaria

            ."). Malaria is endemic in Thailand, especially in the international border areas, because the geographical structure of these areas leads to forest malaria[3](/articles/s41598-022-17646-6#ref-CR3 "Parker, D. M., Carrara, V. I., Pukrittayakamee, S., McGready, R. & Nosten, F. H. Malaria ecology along the Thailand-Myanmar border. Malar. J. 14, 388 (2015)."),[4](/articles/s41598-022-17646-6#ref-CR4 "Marasri, N., Overgaard, H. J., Sumarnrote, A., Thanispong, K. & Corbel, V. Abundance and distribution of Anopheles mosquitoes in a malaria endemic area along the Thai-Lao border Abundance and distribution of Anopheles mosquitoes in a malaria endemic area along the Thai-Lao border. J. Vector Ecol. 42, 325–334 (2017)."). Forest malaria is a term that describes the forest areas associated with malaria transmission or the potential for transmission[5](/articles/s41598-022-17646-6#ref-CR5 "Kar, N. P., Kumar, A., Singh, O. P., Carlton, J. M. & Nanda, N. A review of malaria transmission dynamics in forest ecosystems. Parasit. Vectors 7, 265 (2014)."). Some _Anopheles_ mosquitoes (Diptera: Culicidae) are the vectors that transmit _Plasmodium_ pathogens to humans via the bite of an infected female _Anopheles_[6](/articles/s41598-022-17646-6#ref-CR6 "Vantaux, A. et al. Anopheles ecology, genetics and malaria transmission in northern Cambodia. Sci. Rep. 11, 6458 (2021).").

Two Anopheles species complexes, Dirus and Minimus, are the most important malaria vectors in Thailand7. A species complex is a morphologically similar group of organisms, making it difficult to identify their members by gold standard methods based on morphological characteristics8. Accurate identification of Anopheles species is important in order to understand their true bionomic characters and behavioral traits, facilitating the establishment of crucial guidelines for natural mosquito population control and the protection of communities against malaria9.

Malaria is endemic in the Thai–Cambodia border areas10. Anopheles (Cellia) dirus (previously known as species A) and An. baimaii (previously known as species D) are sibling species within the Dirus complex belonging to the Leucosphyrus group in the Neomyzomyia series[11](#ref-CR11 "Harbach, R. E. Anopheles Classification. http://www.mosquito-taxonomic-inventory.info/node/11358.11358

             (2022)."),[12](#ref-CR12 "Green, C. A., Munstermann, L. E., Tan, S. G., Panyim, S. & Baimai, V. Population genetic evidence for species A, B, C and D of the Anopheles dirus complex in Thailand and enzyme electromorphs for their identification. Med. Vet. Entomol. 6, 29–36 (1992)."),[13](/articles/s41598-022-17646-6#ref-CR13 "Sallum, M. A. M., Peyton, E. L. & Wilkerson, R. C. Six new species of the Anopheles leucosphyrus group, reinterpretation of An. elegans and vector implications. Med. Vet. Entomol. 19, 158–199 (2005)."), and have been incriminated as primary malaria vectors in Thailand and some countries in Southeast Asia[14](/articles/s41598-022-17646-6#ref-CR14 "Tainchum, K., Kongmee, M., Manguin, S., Bangs, M. J. & Chareonviriyaphap, T. Anopheles species diversity and distribution of the malaria vectors of Thailand. Trends Parasitol. 31, 109–119 (2015)."). Previously, _An. dirus_ was found in sympatry with _An. baimaii_ in the Thai–Myanmar and Thai–Cambodia border areas[14](/articles/s41598-022-17646-6#ref-CR14 "Tainchum, K., Kongmee, M., Manguin, S., Bangs, M. J. & Chareonviriyaphap, T. Anopheles species diversity and distribution of the malaria vectors of Thailand. Trends Parasitol. 31, 109–119 (2015)."),[15](/articles/s41598-022-17646-6#ref-CR15 "Obsomer, V., Defourny, P. & Coosemans, M. The Anopheles dirus complex: spatial distribution and environmental drivers. Malar. J. 6, 26 (2007)."). However, these species have very similar morphological and anatomical features, making it difficult to separate the species, and leading to errors in local malaria vector information[16](/articles/s41598-022-17646-6#ref-CR16 "World Health Organization. Anopheline Species Complexes in South and South-East Asia. 
              https://apps.who.int/iris/bitstream/handle/10665/204779/B2406.pdf?sequence=1&isAllowed=y
              
             (2007)."). The wings of females of the two malaria vector species differ slightly in morphological features. _Anopheles dirus_ has a presector dark (PSD) spot on the radius wing vein that extends basally beyond the PSD spot on the costa, reaching the humeral dark of the costa vein or at least beyond the middle of the presector pale spot (PSP) of the costa vein, whereas _An. baimaii_ has the PSD spot on the radius vein at the level of the PSD spot on the costa or extending only trivially basally, usually to no more than the middle of the PSP of the costa vein[17](/articles/s41598-022-17646-6#ref-CR17 "Rattanarithikul, R., Harrison, B. A., Harbach, R. E., Panthusiri, P. & Coleman, R. E. (2006) Illustrated keys to the mosquitoes of Thailand IV. Anopheles. Southeast Asian J. Trop. Med. Public Health. 
              https://doi.org/10.1088/0004-637X/709/2/937
              
            "). However, a previous _Anopheles_ survey in Kanchanaburi, western Thailand, reported that morphology-based species identification of _An. dirus_ and _An. baimaii_ is often highly error-prone[18](/articles/s41598-022-17646-6#ref-CR18 "Tananchai, C., Tisgratog, R., Manguin, S. & Chareonviriyaphap, T. Species diversity and biting activity of Anopheles dirus and Anopheles baimaii ( Diptera : Culicidae ) in a malaria prone area of western Thailand. Parasit. Vectors 5, 211 (2012)."). Therefore, molecular methods are widely used for identifying _An. dirus_ and _An. baimaii_[19](/articles/s41598-022-17646-6#ref-CR19 "Huong, N. T., Sonthayanon, P., Ketterman, A. J. & Panyim, S. A. rapid polymerase chain reaction based method for identification of the Anopheles dirus sibling species. Southeast Asian J. Trop. Med. Public Health 32, 615–620 (2001)."). Multiplex species-specific polymerase chain reaction (PCR) assay based on the internal transcribed spacer 2 region (ITS2) of ribosomal DNA is currently used as the standard molecular method for species identification of sibling members within the Dirus complex[19](/articles/s41598-022-17646-6#ref-CR19 "Huong, N. T., Sonthayanon, P., Ketterman, A. J. & Panyim, S. A. rapid polymerase chain reaction based method for identification of the Anopheles dirus sibling species. Southeast Asian J. Trop. Med. Public Health 32, 615–620 (2001)."),[20](/articles/s41598-022-17646-6#ref-CR20 "Manguin, S. et al. SCAR markers and multiplex PCR-based identification of isomorphic species in the Anopheles dirus complex in Southeast Asia. Med. Vet. Entomol. 16, 46–54 (2002)."), but there are several situations in which alternative methods, such as DNA barcoding and geometric morphometrics, are needed.

DNA barcoding is a molecular biology technique commonly used to identify unknown mosquito species, such as specimens which are ambiguous due to damaged morphology caused by field specimen collections21. This molecular technique is based on the nucleotide sequences of a standard short DNA fragment (~ 400–800 bp) with low intragroup genetic divergence but high intergroup genetic divergence22,23. The cytochrome c oxidase subunit I (COI) gene is widely used as the universal barcode locus for the identification of mosquito species24,25. Several studies have reported the success of DNA barcoding in solving the taxonomic problem of mosquito species identification in many countries26,27,28.

Geometric morphometrics (GM)is a modern morphometric measurement technique developed from measurement standards, which relies on geometric principles to assess size and shape variables for species identification or the study of morphological variations29,30. The highlights of the GM technique are that it is inexpensive, rapid, and does not require sophisticated equipment31. This technique has been used to distinguish members of groups of insects that are difficult to identify morphologically, such as horse flies (Tabanidae)32, stable flies (Muscidae)33, muscid flies (Muscidae)34, soldier castes (Rhinotermitidae)35, and mosquitoes (Culicidae)30,36,37. Recently, GM techniques were successful in identifying species members within malaria vector groups in Thailand, including the Maculatus group38 and the Minimus complex39.

In the present study, DNA barcoding and geometric morphometrics were used to assess the effectiveness of identifying An. dirus and An. baimaii, sibling species of the Dirus complex, in the international border areas of the Trat Province in eastern Thailand. Before proceeding with alternative methods, the actual species of all samples was verified using multiplex PCR based on the ITS2 region as a gold standard molecular marker for the identification of members of the Dirus complex. The results of this study provide guidelines for the correct identification of two principal malaria vector species that are difficult to distinguish, making entomological surveillance information more accurate, and leading to further effective mosquito control planning.

Methods

Ethics statement

Our study was conducted in strict accordance with the guidelines and protocols for animal care and use of the Suan Sunandha Rajabhat University, Thailand. All experimental procedures involving mosquitoes were monitored and approved by the Ethics Committee of the Institutional Animal Care and Use Committee of Suan Sunandha Rajabhat University, Bangkok, Thailand (Approval No. IACOC 64-006/2021).

Study sites and mosquito collection

Anopheles mosquitoes were collected from Trat Province, located in the eastern Thai– Cambodia border. Collection was carried out at two sites: the mainland Thai–Cambodia border (12°28′47.2ʺ N 102°40′51.8ʺ E, Dan Chumphon Sub-district, Bo Rai District), and the island Thai–Cambodia border (11°39′21.6ʺ N 102°34′52.6ʺ E, Koh Kood Sub-district, Koh Kood District [Kood island]) (Fig. 1A). Both sites are malaria hot spot areas, according to the malaria report of the Ministry of Public Health of Thailand (accessed at http://malaria.ddc.moph.go.th/). Ten BG-Pro CDC-style mosquito traps (BioGents, Regensbourg, Germany) supplemented with BG-lure cartridges (BioGents) and dry ice (solid carbon dioxide) were used at each border site to catch Anopheles mosquitoes (Fig. 1B). The mosquito collection was carried out simultaneously at the two sites over five consecutive nights monthly from 18:00 to 06:00 during June and August 2021. The field-collected mosquitoes were transported back to the laboratory at the College of Allied Health Sciences, Suan Sunandha Rajabhat University, Thailand for species identification. Before the identification of Anopheles species, all mosquito specimens were kept individually in 1.5 ml Eppendorf tubes and stored in a freezer at − 20 °C.

Figure 1

Map of the mosquito collection sites (A) including the mainland Thai–Cambodia border (red circle) and the island Thai–Cambodia border (orange circle). This map was generated in Google Earth Pro v 7.1.8 (https://earth.google.com). The BG-Pro CDC-style mosquito trap for the mosquito collecting in this study (B).

Mosquito identification

Anopheles mosquitoes were identified using a morphologic taxonomic key[17](/articles/s41598-022-17646-6#ref-CR17 "Rattanarithikul, R., Harrison, B. A., Harbach, R. E., Panthusiri, P. & Coleman, R. E. (2006) Illustrated keys to the mosquitoes of Thailand IV. Anopheles. Southeast Asian J. Trop. Med. Public Health. https://doi.org/10.1088/0004-637X/709/2/937

            "), except for the important malaria vector taxa that are groups and complexes, including the Minimus complex, the Maculatus group, the Dirus complex, and the Barbirostris complex. Member species of the Minimus complex and closely related species (including _An. minimus_ sensu stricto \[s.s.\], _An. harrisoni, An. aconitus, An. varuna_ and _An. pampanai_)[40](/articles/s41598-022-17646-6#ref-CR40 "Garros, C., Koekemoer, L. L., Coetzee, M., Coosemans, M. & Manguin, S. A single multiplex assay to identify major malaria vectors within the African Anopheles funestus and the oriental An. minimus groups. Am. J. Trop. Med. Hyg. 70, 583–590 (2004)."), the Maculatus group (including _An. maculatus_ s.s., _An. pseudowillmori, An. sawadwongporni, An. rampae_ and _An. dravidicus_)[41](/articles/s41598-022-17646-6#ref-CR41 "Walton, C. et al. Genetic diversity and molecular identification of mosquito species in the Anopheles maculatus group using the ITS2 region of rDNA. Infect. Genet. Evol. 7, 93–102 (2007)."), and the Dirus complex (including _An. dirus_ s.s., _An. scanloni, An. cracens, An. baimaii_ and _An. nemophilous_)[42](/articles/s41598-022-17646-6#ref-CR42 "Walton, C. et al. Identification of five species of the Anopheles dirus complex from Thailand, using allele-specific polymerase chain reaction. Med. Vet. Entomol. 13, 24–32 (1999).") were identified and confirmed using the multiplex PCR based on ribosomal DNA ITS2 species-specific primers, following previously published protocols[40](#ref-CR40 "Garros, C., Koekemoer, L. L., Coetzee, M., Coosemans, M. & Manguin, S. A single multiplex assay to identify major malaria vectors within the African Anopheles funestus and the oriental An. minimus groups. Am. J. Trop. Med. Hyg. 70, 583–590 (2004)."),[41](#ref-CR41 "Walton, C. et al. Genetic diversity and molecular identification of mosquito species in the Anopheles maculatus group using the ITS2 region of rDNA. Infect. Genet. Evol. 7, 93–102 (2007)."),[42](/articles/s41598-022-17646-6#ref-CR42 "Walton, C. et al. Identification of five species of the Anopheles dirus complex from Thailand, using allele-specific polymerase chain reaction. Med. Vet. Entomol. 13, 24–32 (1999)."). Members of the Barbirostris complex, including _An. barbirostris_ s.s., _An. dissidens, An. saeungae, An. wejchoochotei,_ and _An. barbirostris_ A3 were identified using species-specific primers based on the _COI_ gene, according to the protocol proposed by Wilai et al. 2020[43](/articles/s41598-022-17646-6#ref-CR43 "Wilai, P. et al. A multiplex PCR based on mitochondrial COI sequences for identification of members of the Anopheles barbirostris complex (Diptera: Culicidae) in Thailand and other countries in the region. Insects 11, 409 (2020)."). Total genomic DNA was extracted from at least two legs of individual adult _Anopheles_ specimens using FavorPrep™ Mini Kits (Favorgen Biotech, Ping-Tung, Taiwan) following the manufacturer’s protocols. Platinum Taq DNA polymerase (Invitrogen, Carlsbad, CA, USA) was used in all PCR multiplex assays. PCR amplicons were separated by electrophoresis on 2% agarose gels stained with Midori Green DNA stain (Nippon Gene, Tokyo, Japan) and species-specific bands of DNA were identified using an ImageQuant LAS 500 imager (GE Healthcare Japan Corp., Tokyo, Japan) for true species identification.

DNA barcoding

To distinguish An. dirus and An. baimaii using DNA barcoding**,** all PCR reactions of 20 μL contained 1 μL genomic DNA, 1 × reaction buffer, 3 mM MgCl2, 0.4 Unit of Platinum Taq DNA polymerase (Invitrogen), 0.2 µM each of forward and reverse primers, 0.2 mM dNTPs, and the remaining volume of distilled water. Universal barcode primer pairs including forward (5'-GGA TTT GGA AAT TGA TTA GTT CCT T-3') and reverse (5'-AAA AAT TTT AAT TCC AGT TGG AAC AGC-3') primers were used to amplify around 707 bp fragments of the COI gene44, under the following PCR conditions: an initial denaturation at 95 °C for 5 min, followed by five cycles of 94 °C for 40 s, 45 °C for 60 s, 72 °C for 1 min, then 35 cycles of 94 °C for 40 s, 54 °C for 60 s, 72 °C for 1 min, and a final extension of 72 °C for 10 min. Negative controls were included in every PCR round. Amplicons were visualized on 1.5% agarose gels with Midori Green DNA stain using an ImageQuant LAS 500 imager. PCR products with clear DNA bands on agarose gels were sent to Solgent Company, Daejeon, South Korea for DNA sequencing. The forward and reverse COI sequence chromatographs obtained were inspected and edited manually using BioEdit software45 to generate a consensus sequence for each individual Anopheles sample. All consensus sequences of An. dirus and An. baimaii were compared with those available in the global sequence database to confirm the species, including the NCBI GenBank using the Basic Local Alignment Search Tool (BLAST) and the Barcode of Life Data Systems (BOLD) using the BOLD Identification System. After that the COI sequences were automatically aligned using the Clustal X46, the genetic distances were computed, including intra- and interspecific genetic divergences based on the Kimura 2-parameter distance model (K2P) and phylogenetic trees constructed using MEGA X software47. The maximum likelihood (ML) and neighbor joining (NJ) methods were used for phylogenetic analyses. The node support of the phylogenetic tree was assessed using bootstrap analysis with 1000 replicates. The Tamura 3-parameter with Gamma distribution (G) was the most appropriate nucleotide substitution model for producing the ML phylogenetic tree.

Wing geometric morphometrics

All complete right wings of female An. dirus and An. baimaii were detached from the thorax, and the wing scales on the veins were removed and mounted on microscope slides under glass coverslips with Hoyer’s mounting medium. The microscope slides of wings were then photographed using a digital camera attached to a Nikon SMZ 800 N stereo-microscope (Nikon Corp., Tokyo, Japan). Subsequently, the coordinates of 18 landmarks (Fig. 2A) were digitized from each wing image. The landmark locations used in this study were based on previous research studies48,49.

Figure 2

Location of 18 landmarks digitized on the right wing of a female Anopheles mosquito (A), superposition of the mean landmark configurations (B), mean superposition line formed by landmark coordinates at positions 17 and 18 (C), and factor map based on discriminant analysis between Anopheles dirus (red) and An. baimaii (blue) (D).

XY Online Morphometrics (XYOM)50, an online tool, was used for landmark digitization, GM analyses, and graphical outputs in this study. Ten images per species were selected at random and digitized twice by the same person to test the repeatability of the landmarks and evaluate the accuracy of landmark collection in a data set51. To assess allometry—the relationship between wing size and wing shape—linear regression of the first (shape derived) discriminant factor on the wing size was performed, and the coefficient of determination was subsequently estimated.

For wing size analyses, the global wing size was estimated from the centroid size (CS), which is defined as the square root of the sum of the squared distances between each landmark and the centroid of the configuration52. The CS was calculated after the Procrustes superimposition of the landmark configurations. The wing CS difference between both species was analyzed using one-way ANOVA. The statistical significance of the one-way ANOVA was estimated using a non-parametric procedure (1000 runs) with a threshold of significance of _p_-value < 0.05. Discriminant analysis (or canonical variate analysis) was used to assess wing shape variation among species using a factor map and then the Mahalanobis distance was computed in order to estimate the divergence in shape between the species. The final shape variables were computed using Generalized Procrustes Analysis and used as input to a discriminant analysis. Wing shape differences based on the Mahalanobis distances between An. dirus and An. baimaii were estimated using a non-parametric procedure (1000 runs) at a _p_-value < 0.05. To assess the correct assignment of An. dirus and An. baimaii and confirm the reliability of the GM tool when distinguishing between the species based on wing shape, a cross-validated reclassification procedure was carried out based on allocating each individual into the closest group53.

Results

Two hundred and sixty-four adult Anopheles mosquito samples, representing 11 species, of which three species belonged to the subgenus Anopheles and eight to the subgenus Cellia, were collected from the Thai–Cambodia border areas of Trat Province, Thailand during June and August 2021 (Table 1). Adult An. jamesii, An. philippinensis, and An. pseudojamesi specimens were morphologically identified. Other Anopheles mosquitoes were identified to the species using multiplex PCR. Ten Anopheles spp. were collected from the mainland border with the exception of An. baimaii, which was found in the island border (16.29% of the total Anopheles captured in this study). The PCR product of the multiplex PCR assay for the identification of member species in the Dirus complex is shown in Fig. 3A.

Table 1 Anopheles species and numbers collected in the Thai–Cambodia border areas, Trat Province, Thailand, and identification methods used.

Full size table

Figure 3

PCR product of the multiplex PCR assay based on the ITS2 region for species identification (A). Lane 1: 3000 bp molecular ladder; lanes 2–6: An. dirus (562 bp); lanes 7–11: An. baimaii (306 bp); lane 12: negative control. Maximum likelihood tree based on the COI sequences of An. dirus (red) and An. baimaii (blue) in this study, and reference Anopheles species sequences obtained from GenBank, with Culex gelidus as the outgroup taxon (B). Bootstrap values (1000 replicates) of maximum likelihood (green) and neighbor joining (purple) higher than 50% are shown at the nodes.

DNA barcoding for species identification

After confirmation of the true species, partial sequences of the COI gene were amplified using universal barcode primers from ten randomly selected specimens of An. dirus and An. baimaii. No insertions, deletions, stop codons, or pseudogenes were found in any of the COI sequences, when the chromatograms were checked. The average nucleotide composition percentages of the COI sequences of An. dirus and An. baimaii were A = 30.3%, T = 39%, C = 15%, and G = 15.7%, and the sequences were AT rich, with an average AT content of 69.3%, higher than the GC content (30.7%). The average intraspecific genetic divergence percentages of An. dirus was 0.51%, ranging from 0% to 1.14%, and that of An. baimaii was 0.35%, ranging from 0% to 0.71%. The average interspecific divergence between the species was 0.47%, ranging from 0% to 0.99%. Thus, the intra- and interspecific divergence showed overlaps between the two Anopheles spp., indicating no barcoding gap. For identification based on the online genetic database, the COI sequences of An. dirus and An. baimaii from this study were compared with those available in the NCBI using BLAST, and in the BOLD database and were shown to have identity overlaps between the species.

Our COI sequences in this study, and several sequences available from GenBank, were used for construction of a phylogenetic tree (Table 2). The ML tree showed the phylogenetic relationships of the members of the Dirus complex (Fig. 3B). All sequences of An. dirus and An. baimaii from this study and GenBank were unambiguously assigned to the same clade. Other members of the Dirus complex, including An. cracens, An. nemophilous, An. elegans, and An. scanloni showed a clear-cut separation between the species. Culex gelidus, as an outgroup species, was clearly separated from clusters of species members within the Dirus complex.

Table 2 GenBank accession numbers of the COI of Anopheles dirus and An. baimaii sequences in this study and reference mosquito species used for construction of the phylogenetic tree.

Full size table

Geometric morphometrics for species identification

Thirty wings from An. dirus and thirty-six wings from An. baimaii were used for GM analyses. The precision and accuracy of the digitization of landmarks in the wing image sets showed good scores, based on their repeatability. The repeatability scores were more than 95% for size and shape. Assessment of the allometric effect showed that our sample exhibited no apparent relationship between size and shape (r2 = 0.2%), which was apparent from the relatively straight line of the linear regression (Fig. 4A).

Figure 4

Scatter plot showing the linear regression prediction (orange dotted line) for examining the allometric effect (A), and boxplots showing the variation of wing centroid sizes (CS) of Anopheles dirus (red) and An. baimaii (blue) (B). Each quartile box displays the median as a vertical line across the middle and the quartiles indicating the 25th and 75th percentiles at its ends.

The wing CS variations of An. dirus and An. baimaii are shown in Fig. 4B. The mean wing CS of An. dirus was 2.89 mm, ranging from 2.68–3.11 mm, while An. baimaii was 2.90 mm, ranging from 2.74–3.07 mm (Table 3). No significant difference in wing CS between species was found (p > 0.05).

Table 3 Mean wing centroid size of Anopheles dirus and An. baimaii, and statistical characteristics_._

Full size table

For wing shape analyses, visualizations of the wing shapes from the superposition of the mean landmark configurations were created before the extraction of the final shape variables (Fig. 2B). The superimposed landmark coordinates at positions 17 and 18 varied the most between An. dirus and An. baimaii (Fig. 2C). The factor map based on discriminant analysis of the wing shape showed that no samples of the species overlapped (Fig. 2D). The Mahalanobis distance between the species was 5.86, and was statistically significant (p < 0.05). These results showed that the wing shapes of An. dirus and An. baimaii were clearly different. A pairwise cross-validated reclassification test showed excellent accuracy scores for correct identification, ranging from 91.67 to 93.33% (Table 4).

Table 4 Cross-validated scores based on the wing shapes of Anopheles dirus and An. baimaii.

Full size table

Discussion

A recent study by Pimnon and Bhumiratana (2018) in the Thai–Cambodia border area of Trat Province, Thailand10 investigated the distribution of six Anopheles spp.: An. dirus, An. minimus, An. pseudowillmori, An. campestris, An. barbirostris, and An. jamesii. However, in this study, 11 Anopheles spp., including An. dirus, An. aconitus, An. baimaii, An. jamesii, An. minimus, An. dissidens, An. maculatus, An. philippinensis, An. pseudojamesi, An. saeungae, and An. wejchoochotei, were collected from the Thai–Cambodia border. An increasing number of Anopheles spp. were reported in this survey, due to the use of molecular biology techniques in the investigation. In particular, sibling members within the Anopheles species complex are difficult to identify species by standard morphological methods43,54. In the present study, three sibling species of the Barbirostris complex, An. dissidens, An. saeungae, and An. wejchoochotei, and two sibling species of the Dirus complex, An. dirus and An. baimaii, were identified using species-specific multiplex PCR. This study is the first time that An. baimaii was found in an island area, Thailand. This Anopheles spp. has been reported to have a distribution restricted to some forested areas of central, southern and western Thailand, especially the Thai–Myanmar border14.

In Thailand, five species of the Dirus complex are reported: An. dirus (species A), An. cracens (species B), An. scanloni (species C), An. baimaii (species D), and An. nemophilous (species F)9,13,[16](/articles/s41598-022-17646-6#ref-CR16 "World Health Organization. Anopheline Species Complexes in South and South-East Asia. https://apps.who.int/iris/bitstream/handle/10665/204779/B2406.pdf?sequence=1&isAllowed=y

             (2007)."). However, only _An. baimaii_ and _An. dirus_ are considered as the primary vectors of malaria in Thailand, and can be found along the Thai–Cambodia border. Tananchai et al. 2012[16](/articles/s41598-022-17646-6#ref-CR16 "World Health Organization. Anopheline Species Complexes in South and South-East Asia. 
              https://apps.who.int/iris/bitstream/handle/10665/204779/B2406.pdf?sequence=1&isAllowed=y
              
             (2007).") studied the biology of two _Anopheles_ spp. in Pu Teuy Village, Kanchanaburi Province, western Thailand, and found that they had different patterns of feeding activity. Their study was the starting point for understanding the behavioral differences between the two vector species, and formed the basis of strategies for malaria control and prevention in areas in which malaria is endemic. However, a variety of effective alternative techniques are needed to aid in the identification of both species, according to the appropriate situation.

DNA barcoding is widely recognized as a modern molecular technique which can help identify unknown mosquito specimens to the species level55. However, the results of this study revealed that DNA barcoding based on the COI gene, as the “gold standard” barcode for animals24,25, could not be used to discriminate An. dirus from An. baimaii. The overlapping values for intra- and interspecific genetic divergence indicated that there was no barcoding gap for An. dirus and An. baimaii. The lack of a barcoding gap means that the genetic distances of our COI sequences cannot be used to identify An. dirus and An. baimaii, due to insufficient nucleotide differences making it impossible to delineate species boundaries56. This finding was consistent with phylogenetic analysis based on the COI sequences, which found that the species were not grouped separately into clades, while the other members of the Dirus complex, including An. cracens, An. nemophilous, An. elegans, and An. scanloni were clearly separated. Therefore, although this technique cannot separate An. dirus and An. baimaii, it can effectively classify the other species members of the Dirus complex. Our results are consistent with previous studies reporting that partial COI gene sequences failed in distinguishing between An. dirus and An. baimaii but can clearly distinguish other species of the Dirus complex and the Leucosphyrus group57. Often, this technique has failed to identify members in species complexes or closely related species, due to overlapping intra- and interspecific variation caused by incomplete lineage sorting between recently diverged species, such as some member species in the Anopheles strodei subgroup58.

In contrast, the results of GM based on wing shape analysis revealed an excellent differentiation between An. dirus and An. baimaii, with 92.42% of specimens correctly assigned based on a cross-validated reclassification procedure. The discriminant analysis identified the differences between the wing shapes of An. dirus and An. baimaii. Visual wing shape based on superposition of the mean landmark configurations showed an incompatibility between An. dirus and An. baimaii at landmark positions 17 and 18. These results suggest that the wing venation of each species is unique, leading to the successful identification based on wing GM analyses. The shape of the wing was recognized as a valuable variable for GM analyses in identifying mosquito species due to its stable character and distinctive structure31,36. These findings are in agreement with previous studies which found that GM techniques based on wing shape analyses can assist in the identification of members among the Anopheles species group and complex in Thailand, such as two species of the Minimus complex, An. minimus and An. harrisoni39, and three species of the Maculatus group, An. maculatus, An. sawadwongporni, and An. pseudowillmori38. Wing size is a variable that is not suitable for species identification by GM because it varies according to environmental influences59,60. Our analyses of the wing size based on CS in this study indicated that An. dirus and An. baimaii were not statistically different. Consistent with the results of assessing allometry, the size and shape of the wings had no relationship. Therefore, the wing size of An. dirus and An. baimaii was not used for species identification.

Conclusions

The findings of this study revealed that the GM technique can be used to effectively distinguish An. dirus from An. baimaii. This modern technique is valuable in laboratories with budget or molecular biology equipment constraints. However, the use of the GM technique should be revalidated each time in every other target area to prevent errors from morphological variations of mosquito specimens in response to environmental factors. In addition, further studies with a large number of specimens are needed to increase the confidence in this technique. For molecular techniques, our results show that DNA barcoding cannot distinguish An. dirus from An. baimaii due to the low difference in the nucleotide sequences of the COI gene of the species. However, the multiplex PCR method based on ITS2 remains the most effective in identifying both species. This study is an important guideline for identifying An. dirus and An. baimaii sibling species within the Dirus complex. These species are important malaria vectors in the border areas, and the ability to identify them will lead to the development of further effective strategies for the control of the Anopheles population in Thailand and other countries in Southeast Asia.

Data availability

All DNA barcode sequences of An. dirus and An. baimaii were submitted and available in the GeneBank Database (https://www.ncbi.nlm.nih.gov/nuccore), accession numbers: OM423762- OM423781 (Table 2).

References

  1. Sato, S. _Plasmodium_—a brief introduction to the parasites causing human malaria and their basic biology. J Physiol Anthropol. 40, 1–13 (2021).
    Article PubMed PubMed Central Google Scholar
  2. World Health Organization. WHO malaria report. 2021; https://www.who.int/news-room/fact-sheets/detail/malaria.
  3. Parker, D. M., Carrara, V. I., Pukrittayakamee, S., McGready, R. & Nosten, F. H. Malaria ecology along the Thailand-Myanmar border. Malar. J. 14, 388 (2015).
    Article PubMed PubMed Central Google Scholar
  4. Marasri, N., Overgaard, H. J., Sumarnrote, A., Thanispong, K. & Corbel, V. Abundance and distribution of Anopheles mosquitoes in a malaria endemic area along the Thai-Lao border Abundance and distribution of Anopheles mosquitoes in a malaria endemic area along the Thai-Lao border. J. Vector Ecol. 42, 325–334 (2017).
    Article PubMed Google Scholar
  5. Kar, N. P., Kumar, A., Singh, O. P., Carlton, J. M. & Nanda, N. A review of malaria transmission dynamics in forest ecosystems. Parasit. Vectors 7, 265 (2014).
    Article PubMed PubMed Central Google Scholar
  6. Vantaux, A. et al. Anopheles ecology, genetics and malaria transmission in northern Cambodia. Sci. Rep. 11, 6458 (2021).
    Article ADS CAS PubMed PubMed Central Google Scholar
  7. Sriwichai, P. et al. Natural human Plasmodium infections in major Anopheles mosquitoes in western Thailand. Parasit. Vectors 9, 17 (2016).
    Article PubMed PubMed Central CAS Google Scholar
  8. Sumruayphol, S. et al. Seasonal dynamics and molecular differentiation of three natural Anopheles species (Diptera: Culicidae) of the Maculatus group (Neocellia series) in malaria hotspot villages of Thailand. Parasit. Vectors 13, 574 (2020).
    Article PubMed PubMed Central Google Scholar
  9. Tananchai, C., Manguin, S., Bangs, M. J. & Chareonviriyaphap, T. Malaria vectors and species complexes in Thailand: implications for vector control. Trends Parasitol. Trends Parasitol. 35, 544–558 (2019).
    Article PubMed Google Scholar
  10. Pimnon, S. & Bhumiratana, A. Adaptation of Anopheles vectors to anthropogenic malaria-associated rubber plantations and indoor residual spraying: Establishing population dynamics and insecticide susceptibility. Can. J. Infect. Dis. Med. Microbiol. 2018, 9853409 (2018).
    Article PubMed PubMed Central Google Scholar
  11. Harbach, R. E. Anopheles Classification. http://www.mosquito-taxonomic-inventory.info/node/11358.11358 (2022).
  12. Green, C. A., Munstermann, L. E., Tan, S. G., Panyim, S. & Baimai, V. Population genetic evidence for species A, B, C and D of the Anopheles dirus complex in Thailand and enzyme electromorphs for their identification. Med. Vet. Entomol. 6, 29–36 (1992).
    Article CAS PubMed Google Scholar
  13. Sallum, M. A. M., Peyton, E. L. & Wilkerson, R. C. Six new species of the Anopheles leucosphyrus group, reinterpretation of An. elegans and vector implications. Med. Vet. Entomol. 19, 158–199 (2005).
    Article CAS PubMed Google Scholar
  14. Tainchum, K., Kongmee, M., Manguin, S., Bangs, M. J. & Chareonviriyaphap, T. Anopheles species diversity and distribution of the malaria vectors of Thailand. Trends Parasitol. 31, 109–119 (2015).
    Article PubMed Google Scholar
  15. Obsomer, V., Defourny, P. & Coosemans, M. The Anopheles dirus complex: spatial distribution and environmental drivers. Malar. J. 6, 26 (2007).
    Article PubMed PubMed Central Google Scholar
  16. World Health Organization. Anopheline Species Complexes in South and South-East Asia. https://apps.who.int/iris/bitstream/handle/10665/204779/B2406.pdf?sequence=1&isAllowed=y (2007).
  17. Rattanarithikul, R., Harrison, B. A., Harbach, R. E., Panthusiri, P. & Coleman, R. E. (2006) Illustrated keys to the mosquitoes of Thailand IV. Anopheles. Southeast Asian J. Trop. Med. Public Health. https://doi.org/10.1088/0004-637X/709/2/937
  18. Tananchai, C., Tisgratog, R., Manguin, S. & Chareonviriyaphap, T. Species diversity and biting activity of Anopheles dirus and Anopheles baimaii ( Diptera : Culicidae ) in a malaria prone area of western Thailand. Parasit. Vectors 5, 211 (2012).
    Article PubMed PubMed Central Google Scholar
  19. Huong, N. T., Sonthayanon, P., Ketterman, A. J. & Panyim, S. A. rapid polymerase chain reaction based method for identification of the Anopheles dirus sibling species. Southeast Asian J. Trop. Med. Public Health 32, 615–620 (2001).
    CAS PubMed Google Scholar
  20. Manguin, S. et al. SCAR markers and multiplex PCR-based identification of isomorphic species in the Anopheles dirus complex in Southeast Asia. Med. Vet. Entomol. 16, 46–54 (2002).
    Article CAS PubMed Google Scholar
  21. Adeniran, A. A. et al. Identification of mosquitoes (Diptera: Culicidae) from Mexico State, Mexico using morphology and COI DNA barcoding. Acta Trop. 213, 105730 (2021).
    Article CAS PubMed Google Scholar
  22. Leite, L. A. R. Mitochondrial pseudogenes in insect DNA barcoding: differing points of view on the same issue. Biota Neotrop. 12, 301–308 (2012).
    Article Google Scholar
  23. Batovska, J., Blacket, M. J., Brown, K. & Lynch, S. E. Molecular identification of mosquitoes (Diptera: Culicidae) in southeastern Australia. Ecol. Evol. 28, 3001–3011 (2016).
    Article Google Scholar
  24. Beebe, N. W. DNA barcoding mosquitoes: advice for potential prospectors. Parasitology 145, 622–633 (2018).
    Article PubMed CAS Google Scholar
  25. Maekawa, Y. et al. DNA barcoding of mosquitoes collected through a nationwide survey in 2011 and 2012 in Malawi Southeast Africa. Acta Trop. 213, 105742 (2021).
    Article CAS PubMed Google Scholar
  26. Hernández-Triana, L. M. et al. DNA barcoding of British mosquitoes (Diptera, Culicidae) to support species identification, discovery of cryptic genetic diversity and monitoring invasive species. Zookeys 19, 57–76 (2019).
    Article Google Scholar
  27. Talaga, S. et al. DNA reference libraries of French Guianese mosquitoes for barcoding and metabarcoding. PLoS ONE 12, e0176993 (2017).
    Article PubMed PubMed Central CAS Google Scholar
  28. Ashfaq, M. et al. Analyzing mosquito (Diptera: Culicidae) diversity in Pakistan by DNA barcoding. PLoS ONE 9, e97268 (2014).
    Article ADS PubMed PubMed Central CAS Google Scholar
  29. Henry, A., Thongsripong, P., Fonseca-Gonzalez, I., Jaramillo-Ocampo, N. & Dujardin, J. P. Wing shape of dengue vectors from around the world. Infect. Genet. Evol. 10, 207–214 (2010).
    Article CAS PubMed Google Scholar
  30. Wilke, A. B. B. et al. Morphometric wing characters as a tool for mosquito identification. PLoS ONE 11, e0161643 (2016).
    Article PubMed PubMed Central CAS Google Scholar
  31. Lorenz, C. et al. Geometric morphometrics in mosquitoes: What has been measured?. Infect. Genet. Evol. 54, 205–215 (2017).
    Article PubMed Google Scholar
  32. Changbunjong, T. et al. Landmark data to distinguish and identify morphologically close Tabanus spp. (Diptera: Tabanidae). Insects 12, 974 (2021).
    Article PubMed PubMed Central Google Scholar
  33. Changbunjong, T. et al. Landmark and outline-based geometric morphometrics analysis of three Stomoxys flies (Diptera: Muscidae). Folia Parasitol. 63, 2016.037 (2016).
    Article Google Scholar
  34. Limsopatham, K. et al. Wing morphometrics of medically and forensically important muscid flies (Diptera: Muscidae). Acta Trop. 222, 106062 (2021).
    Article PubMed Google Scholar
  35. Wikantyoso, B., Tseng, S. P., Himmi, S. K., Yusuf, S. & Yoshimura, T. Morphometric analysis of Coptotermes spp. Soldier caste (blattodea: Rhinotermitidae) in Indonesia and evidence of Coptotermes gestroi extreme head-capsule shapes. Insects 12, 477 (2021).
    Article PubMed PubMed Central Google Scholar
  36. Chaiphongpachara, T. Comparison of landmark- and outline-based geometric morphometrics for discriminating mosquito vectors in Ratchaburi Province Thailand. Biomed Res. Int. 2018, 6170502 (2018).
    Article PubMed PubMed Central Google Scholar
  37. Champakaew, D. et al. Geometric morphometric wing analysis as a tool to discriminate female mosquitoes from different suburban areas of Chiang Mai Province Thailand. PLoS One 16, e0260333 (2021).
    Article CAS PubMed PubMed Central Google Scholar
  38. Chaiphongpachara, T., Sriwichai, P., Samung, Y. & Ruangsittichai, J. Geometric morphometrics approach towards discrimination of three member species of Maculatus group in Thailand. Acta Trop. 192, 66–74 (2019).
    Article PubMed PubMed Central Google Scholar
  39. Chatpiyaphat, K. et al. Geometric morphometrics to distinguish the cryptic species Anopheles minimus and An. harrisoni in malaria hot spot villages, western Thailand. Med. Vet. Entomol. 35, 293–301 (2021).
    Article CAS PubMed Google Scholar
  40. Garros, C., Koekemoer, L. L., Coetzee, M., Coosemans, M. & Manguin, S. A single multiplex assay to identify major malaria vectors within the African Anopheles funestus and the oriental An. minimus groups. Am. J. Trop. Med. Hyg. 70, 583–590 (2004).
    Article CAS PubMed Google Scholar
  41. Walton, C. et al. Genetic diversity and molecular identification of mosquito species in the Anopheles maculatus group using the ITS2 region of rDNA. Infect. Genet. Evol. 7, 93–102 (2007).
    Article CAS PubMed Google Scholar
  42. Walton, C. et al. Identification of five species of the Anopheles dirus complex from Thailand, using allele-specific polymerase chain reaction. Med. Vet. Entomol. 13, 24–32 (1999).
    Article CAS PubMed Google Scholar
  43. Wilai, P. et al. A multiplex PCR based on mitochondrial COI sequences for identification of members of the Anopheles barbirostris complex (Diptera: Culicidae) in Thailand and other countries in the region. Insects 11, 409 (2020).
    Article PubMed Central Google Scholar
  44. Kumar, N. P., Rajavel, A. R., Natarajan, R. & Jambulingam, P. DNA barcodes can distinguish species of Indian mosquitoes (Diptera: Culicidae). J. Med. Entomol. 44, 1–7 (2007).
    Article CAS PubMed Google Scholar
  45. Hall, T. A. BIOEDIT: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/ NT. Nucleic Acids Symp. Ser. 41, 95–98 (1999).
    CAS Google Scholar
  46. Larkin, M. A. et al. Clustal W and Clustal X version 2.0. Bioinformatics 23, 2947–2948 (2007).
    Article CAS PubMed Google Scholar
  47. Kumar, S., Stecher, G., Li, M., Knyaz, C. & Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 35, 1547–1549 (2018).
    Article CAS PubMed PubMed Central Google Scholar
  48. Demari-Silva, B., Suesdek, L., Sallum, M. A. M. & Marrelli, M. T. Wing geometry of Culex coronator (Diptera: Culicidae) from South and Southeast Brazil. Parasit. Vectors 7, 174 (2014).
    Article PubMed PubMed Central Google Scholar
  49. Vidal, P. O. & Suesdek, L. Comparison of wing geometry data and genetic data for assessing the population structure of Aedes aegypti. Infect. Genet. Evol. 12, 591–596 (2012).
    Article PubMed Google Scholar
  50. Dujardin, S. & Dujardin, J. P. Geometric morphometrics in the cloud. Infect. Genet. Evol. 70, 189–196 (2019).
    Article PubMed MATH Google Scholar
  51. Arnqvist, G. & Mårtensson, T. Measurement error in geometric morphometrics: Empirical strategies to assess and reduce its impact on measures of shape. Acta Zool. Acad. Sci. Hungaricae 44, 73–96 (1998).
    Google Scholar
  52. Bookstein, F. L. Morphometric tools for landmark data: geometry and biology. Cambridge University Press (1991).
  53. Dujardin, J. .-P. . Al., Kaba, D. & Henry, A. B. The exchangeability of shape. BMC Res. Notes 3, 266 (2010).
    Article PubMed PubMed Central Google Scholar
  54. Brosseau, L. et al. A multiplex PCR assay for the identification of five species of the Anopheles barbirostris complex in Thailand. Parasit. Vectors 12, 223 (2019).
    Article PubMed PubMed Central Google Scholar
  55. Chan, A. et al. DNA barcoding: Complementing morphological identification of mosquito species in Singapore. Parasit. Vectors 7, 569 (2014).
    Article PubMed PubMed Central CAS Google Scholar
  56. Wiemers, M. & Fiedler, K. Does the DNA barcoding gap exist? - A case study in blue butterflies (Lepidoptera: Lycaenidae). Front. Zool. 4, 8 (2007).
    Article PubMed PubMed Central CAS Google Scholar
  57. Takenaka Takano, K. et al. Partial mitochondrial DNA sequences suggest the existence of a cryptic species within the Leucosphyrus group of the genus Anopheles (Diptera: Culicidae), forest malaria vectors, in northern Vietnam. Parasit. Vectors 3, 41 (2010).
    Article CAS Google Scholar
  58. Bourke, B. P., Oliveira, T. P., Suesdek, L., Bergo, E. S. & Sallum, M. A. M. A multi-locus approach to barcoding in the Anopheles strodei subgroup (Diptera: Culicidae). Parasit. Vectors 6, 111 (2013).
    Article CAS PubMed PubMed Central Google Scholar
  59. Chaiphongpachara, T. & Laojun, S. Annual variability of wing morphology in Culex sitiens Wiedemann (Diptera, Culicidae) mosquito vectors from the coastal area of Samut Songkhram Province Thailand. J. Parasitol. Res. 2019, 3978965 (2019).
    Article PubMed PubMed Central Google Scholar
  60. Chaiphongpachara, T. & Laojun, S. Variation over time in wing size and shape of the coastal malaria vector Anopheles (Cellia) epiroticus Linton and Harbach (Diptera: Culicidae) in Samut Songkhram Thailand. J. Adv. Vet. Anim. Res. 6, 208–214 (2019).
    Article PubMed PubMed Central Google Scholar

Download references

Acknowledgements

This research was kindly supported by the Research Fund of the Suan Sunandha Rajabhat University, Thailand. Furthermore, we would like to thank Mr. Teerayoot Nutepsu for Anopheles collections and sample preparation.

Funding

This research was financially supported by the Suan Sunandha Rajabhat University, Thailand.

Author information

Authors and Affiliations

  1. Department of Public Health and Health Promotion, College of Allied Health Sciences, Suan Sunandha Rajabhat University, Bangkok, 10300, Thailand
    Tanawat Chaiphongpachara, Sedthapong Laojun, Nantana Suwandittakul & Kewarin Kuntawong
  2. Department of Pre-Clinic and Applied Animal Science, Faculty of Veterinary Science, Mahidol University, Nakhon Pathom, 73170, Thailand
    Tanasak Changbunjong
  3. The Monitoring and Surveillance Center for Zoonotic Diseases in Wildlife and Exotic Animals (MoZWE), Faculty of Veterinary Science, Mahidol University, Nakhon Pathom, 73170, Thailand
    Tanasak Changbunjong
  4. Department of Medical Entomology, Faculty of Tropical Medicine, Mahidol University, Bangkok, 10400, Thailand
    Suchada Sumruayphol

Authors

  1. Tanawat Chaiphongpachara
  2. Tanasak Changbunjong
  3. Suchada Sumruayphol
  4. Sedthapong Laojun
  5. Nantana Suwandittakul
  6. Kewarin Kuntawong

Contributions

T.C. (Tanawat Chaiphongpachara): Conceptualization, Investigation, Methodology, Formal Analysis, Visualization, Writing—original draft. T.C. (Tanasak Changbunjong): Methodology, Formal Analysis, Visualization, Writing—Review & Editing. S.S.: Investigation, Visualization, Writing—Review & Editing. S.L.: Investigation, Formal Analysis. N.S.: Visualization, Writing—Review & Editing. K.K.: Investigation, Formal Analysis.

Corresponding author

Correspondence toTanawat Chaiphongpachara.

Ethics declarations

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.

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/.

Reprints and permissions

About this article

Cite this article

Chaiphongpachara, T., Changbunjong, T., Sumruayphol, S. et al. Geometric morphometrics versus DNA barcoding for the identification of malaria vectors Anopheles dirus and An. baimaii in the Thai-Cambodia border.Sci Rep 12, 13236 (2022). https://doi.org/10.1038/s41598-022-17646-6

Download citation