Wing morphometric variability of the malaria vector Anopheles (Cellia) epiroticus Linton et Harbach (Diptera: Culicidae) for the duration of the rainy season in coastal areas of Samut Songkhram, Thailand

Population genetic structure and wing geometric morphometrics of the filarial vector Armigeres subalbatus (Diptera: Culicidae) in Thailand

Armigeres subalbatus (Diptera: Culicidae) is a mosquito species of significant medical and veterinary importance. It is widely distributed across Southeast and East Asia and is commonly found throughout Thailand. This study assessed the genetic diversity and population structure of Ar. subalbatus in Thailand using the cytochrome c oxidase subunit I (COI) gene sequences. Additionally, wing shape variations among these populations were examined using geometric morphometrics (GM). Our results demonstrated that the overall haplotype diversity (Hd) was 0.634, and the nucleotide diversity (π) was 0.0019. Significant negative values in neutrality tests (p < 0.05) indicate that the Ar. subalbatus populations in Thailand are undergoing a phase of expansion following a bottleneck event. The mismatch distribution test suggests that the populations may have started expanding approximately 16,678 years ago. Pairwise genetic differentiation among the 12 populations based on Fst revealed significant differences in 32 pairs (p < 0.05), with the degree of differentiation ranging from 0.000 to 0.419. The GM analysis of wing shape also indicated significant differences in nearly all pairs (p < 0.05), except for between populations from Nakhon Pathom and Samut Songkhram, and between those from Chiang Mai and Mae Hong Son, suggesting no significant difference due to their similar environmental settings. These findings enhance our understanding of the population structure and phenotypic adaptations of mosquito vectors, providing vital insights for the formulation of more efficacious vector control strategies.

Influence of insular conditions on wing phenotypic variation in two dominant mosquito vectors, Aedes albopictus and Armigeres subalbatus (Diptera: Culicidae), in the border archipelagos of Thailand

Insects geographically separated into island and mainland populations often exhibit phenotypic variations, a phenomenon known as insular conditions. These conditions can lead to rapid evolutionary changes that affect the morphological characteristics of mosquito vectors. Nevertheless, studies that specifically examine phenotype differences between island and mainland mosquito populations have been limited. In this study, wing variation in size and shape was investigated using the geometric morphometric (GM) technique in two dominant mosquito vectors, Aedes albopictus and Armigeres subalbatus, in the Ranong and Trat archipelagos of Thailand. Significant differences in average wing centroid size (CS) were found in 6 out of 15 population pairs for Ae. albopictus (p < 0.05) and in 5 population pairs for Ar. subalbatus (p < 0.05). After removing the allometric effect, canonical variate analyses (CVA) based on wing shape analysis revealed overlap across all populations for both Ae. albopictus and Ar. subalbatus. However, the statistical analysis indicated that Ar. subalbatus exhibited wing shape differences across all populations (p < 0.05), and most Ae. albopictus populations also displayed distinct wing shapes (p < 0.05), except for the populations from Chang Island and the mainland of Ranong, which showed no significant differences (p > 0.05). These findings enhance our understanding of mosquito adaptability in island regions and provide valuable data for the surveillance and monitoring of vector evolution.

Geometric morphometrics to differentiate species and explore seasonal variation in three Mansonia species (Diptera: Culicidae) in central Thailand and their association with meteorological factors

Mansonia mosquito species are recognised as a significant vector of human pathogens, primarily transmitting the filarial nematode, Brugia malayi. In central Thailand, the three most prevalent Mansonia species are Mansonia annulifera, Mansonia indiana and Mansonia uniformis. This study explored the influence of seasonal changes on the phenotypic variation of these Mansonia species in central Thailand using the geometric morphometrics (GM). To ensure accurate species identification, we integrated GM techniques with DNA barcoding, examining distinctions in both phenotype and genotype among the species. The intraspecific genetic divergence ranged from 0.00% to 1.69%, whereas the interspecific genetic divergence ranged from 10.52% to 16.36%. The clear distinction between intra- and interspecific distances demonstrated the presence of a barcoding gap, confirming the successful differentiation of the three Mansonia mosquito species through DNA barcoding. Similarly, the interspecies GM assessment for classifying Mansonia species demonstrated a high degree of accuracy, with an overall performance of 98.12%. Exploring seasonal variation in the three Mansonia species revealed wing variations across different seasons, and pronounced variations appearing in the cool season. Regarding their association with meteorological factors, Ma. annulifera and Ma. uniformis showed significant positive correlations with temperature (p < 0.05), and Ma. uniformis also displayed a significant negative correlation with atmospheric pressure (p < 0.05). The insights from this study will deepen our understanding of the adaptive patterns of Mansonia mosquitoes in Thailand's central region, paving the way for enhanced disease surveillance related to these vectors.

Insights into the mitochondrial cytochrome oxidase I (mt-COI) gene and wing morphometrics of Anopheles baimaii (Diptera: Culicidae) in malaria-endemic islands of Thailand

Anopheles baimaii (Diptera: Culicidae) significantly contributes to the transmission of parasites causing malaria in Southeast Asia and South Asia. This study examined the morphological (wing shape) and molecular (mitochondrial gene) variations of An. baimaii in four of Thailand's border islands, and also investigated the presence of Plasmodium parasites in these mosquitoes. No Plasmodium infections were detected in the samples. Significant differences in wing shape were observed in most island populations (p < 0.05). A single-linkage tree, constructed using Mahalanobis distances, clustered the populations into two groups based on geographical locations. Genetic variation in An. baimaii was also analyzed through cytochrome c oxidase subunit I (COI) gene sequences. This analysis identified 22 segregating sites and a low nucleotide diversity of 0.004. Furthermore, 18 distinct haplotypes were identified, indicating a high haplotype diversity of 0.825. Neutrality tests for the overall population revealed a significantly negative Fu's Fs value (-5.029), indicating a population expansion. In contrast, Tajima's D yielded a negative value (-1.028) that did not reach statistical significance. The mismatch distribution analysis exhibited a bimodal pattern, and the raggedness index was 0.068, showing no significant discrepancy (p = 0.485) between observed and expected distributions. Pairwise genetic differentiation assessments demonstrated significant differences between all populations (p < 0.05). These findings provide valuable insights into the COI gene and wing morphometric variations in An. baimaii across Thailand's islands, offering critical information for understanding the adaptations of this malaria vector and guiding future comprehensive research.

Genetic Diversity, Haplotype Relationships, and kdr Mutation of Malaria Anopheles Vectors in the Most Plasmodium knowlesi-Endemic Area of Thailand

Plasmodium knowlesi, a malaria parasite that occurs naturally in long-tailed macaques, pig-tailed macaques, and banded leaf monkeys, is currently regarded as the fifth of the human malaria parasites. We aimed to investigate genetic diversity based on the cytochrome c oxidase subunit I (COI) gene, detect Plasmodium parasites, and screen for the voltage-gated sodium channel (VGSC)-mutation-mediated knockdown resistance (kdr) of Anopheles mosquitoes in Ranong province, which is the most P. knowlesi-endemic area in Thailand. One hundred and fourteen Anopheles females belonging to eight species, including An. baimaii (21.05%), An. minimus s.s. (20.17%), An. epiroticus (19.30%), An. jamesii (19.30%), An. maculatus s.s. (13.16%), An. barbirostris A3 (5.26%), An. sawadwongporni (0.88%), and An. aconitus (0.88%), were caught in three geographical regions of Ranong province. None of the Anopheles mosquitoes sampled in this study were infected with Plasmodium parasites. Based on the sequence analysis of COI sequences, An. epiroticus had the highest level of nucleotide diversity (0.012), followed by An. minimus (0.011). In contrast, An. maculatus (0.002) had the lowest level of nucleotide diversity. The Fu's Fs and Tajima's D values of the Anopheles species in Ranong were all negative, except the Tajima's D values of An. minimus (0.077). Screening of VGSC sequences showed no presence of the kdr mutation of Anopheles mosquitoes. Our results could be used to further select effective techniques for controlling Anopheles populations in Thailand's most P. knowlesi-endemic area.

Geometric morphometric and molecular techniques for discriminating among three cryptic species of the Anopheles barbirostris complex (Diptera: Culicidae) in Thailand

Anopheles members of the Barbirostris complex are important vectors of malaria in Thailand. However, they are morphologically indistinguishable because they are closely related species. In this study, wing geometric morphometrics (GM) and DNA barcoding based on the cytochrome c oxidase subunit 1 (COI) gene were applied to differentiate cryptic species of the Barbirostris complex in Thailand. Three cryptic species of the Barbirostris complex, Anopheles dissidens (19.44%), Anopheles saeungae (24.54%), and Anopheles wejchoochotei (56.02%) were initially identified using the multiplex polymerase chain reaction assay. DNA barcoding analyses showed low intraspecific distances (range, 0.27%–0.63%) and high interspecific distances (range, 1.92%–3.68%), consistent with the phylogenetic analyses that showed clear species groups. While wing size and shape analyses based on landmark-based GM indicated differences between three species (p < 0.05). The cross-validated reclassification revealed that the overall efficacy of wing size analysis for species identification of the Barbirostris complex was less than the wing shape analysis (56.43% vs. 74.29% total performance). Therefore, this study's results are guidelines for applying modern techniques to identify members within the Barbirostris complex, which are still difficult to distinguish by morphology-based identification and contribute to further appropriate malaria control.

•

DNA barcoding based on the cytochrome c oxidase subunit 1 (COI) gene is the most reliable identification tool for the Anopheles barbirostris complex.

•

Analysis of wing size and shape of Anopheles dissidens, An. saeungae and An. wejchoochotei based on geometric morphometrics revealed differences between species (p < 0.05).

•

The efficacy of wing shape analysis for species identification of the Barbirostris complex was moderate levels of performance (74.29% accuracy score).

Mitochondrial DNA barcoding of mosquito species (Diptera: Culicidae) in Thailand

The correct identification of mosquito species is important for effective mosquito vector control. However, the standard morphological identification of mosquito species based on the available keys is not easy with specimens in the field due to missing or damaged morphological features during mosquito collections, often leading to the misidentification of morphologically indistinguishable. To resolve this problem, we collected mosquito species across Thailand to gather genetic information, and evaluated the DNA barcoding efficacy for mosquito species identification in Thailand. A total of 310 mosquito samples, representing 73 mosquito species, were amplified using mitochondrial cytochrome c oxidase subunit I (COI) primers. The average maximum intraspecific genetic variation of the 73 mosquito species was 1% ranged from 0–5.7%. While, average minimum interspecific genetic variation (the distance to the nearest neighbour) of the 73 mosquito species was 7% ranged from 0.3–12.9%. The identification of success rates based on the “Best Match,” “Best Close Match,” and “All Species Barcodes” methods were 97.7%, 91.6%, and 81%, respectively. Phylogenetic analyses of Anopheles COI sequences demonstrated a clear separation between almost all species (except for those between An. baimaii and An. dirus), with high bootstrap support values (97%–99%). Furthermore, phylogenetic analyses revealed potential sibling species of An. annularis, An. tessellatus, and An. subpictus in Thailand. Our results indicated that DNA barcoding is an effective molecular approach for the accurate identification of mosquitoes in Thailand.

Intraspecific variation in wing geometry among Tabanus rubidus (Diptera: Tabanidae) populations in Thailand

Tabanus rubidus (Wiedemann, 1821) (Diptera: Tabanidae) is a hematophagous insect of veterinary and medical importance and is the predominant Tabanus spp. in Thailand. It is a potential mechanical vector of Trypanosoma evansi, which causes surra in domestic and wild animals. Wing geometric morphometrics is widely used as morphological markers for species identification and to assess the insect population structure. Herein, we investigated the intraspecific variation in wing geometry among T. rubidus populations in Thailand using landmark-based geometric morphometric analysis. Tabanus rubidus females were collected from five populations in four geographical regions in Thailand. The left wings of 240 specimens were removed and digitized using 22 landmarks for analysis. While wing size variations were found between some populations, wing shape variations were detected in all. These intraspecific variations in T. rubidus populations indicate an adaptive response to the local environmental conditions.

A novel use of a geometric morphometric technique to distinguish human parasite eggs of twelve different species

Copro-microscopic diagnostic methods are the most common approach for screening patients with parasitic infections. However, expertise is required to identify helminthic eggs from fecal specimens. Consequently, new methods are required. Novel technologies have recently been developed for the classification of organisms, including geometric morphometric (GM) approaches. In this study, the outline-based GM approach was used to distinguish the eggs of 12 common human parasite species, including Ascaris lumbricoides, Trichuris trichiura, Enterobius vermicularis, hookworm, Capillaria philippinensis, Opisthorchis spp., Fasciola spp., Paragonimus spp., Schistosoma mekongi, Taenia spp., Hymenolepis diminuta and Hymenolepis nana. The GM analysis revealed that the size cannot be used as the main variable in the identification of parasite species at the egg stage, producing only 30.18% overall accuracy. However, comparisons of shape based on the Mahalanobis distances between pairs of parasite species showed significant differences in all pairs (p < 0.05). The shape analysis produced 84.29% overall accuracy. This is the first time that outline-based GM has been preliminarily confirmed as a valuable approach to support copro-microscopic analysis, in order to effectively screen helminth eggs. However, further studies with a larger set of helminth eggs and artefacts should be carried out to increase confidence in the identification of parasite species in the absence of local experts.

Wing Phenotypic Variation among Stomoxys calcitrans (Diptera: Muscidae) Populations in Thailand

Simple Summary

The stable fly, Stomoxys calcitrans (Diptera: Muscidae), is the predominant Stomoxys spp. in Thailand and is considered a pest for livestock, pets, wildlife, and occasionally humans. This study investigated the phenotypic variation in the wing size and shape of S. calcitrans populations from different geographical regions in Thailand using landmark-based geometric morphometric analysis. Results showed strong spatial variation in wing shape among S. calcitrans populations and thus suggested the existence of phenotypic plasticity in this fly.

Abstract

Stomoxys calcitrans (Linnaeus, 1758) (Diptera: Muscidae) is a cosmopolitan hematophagous ectoparasite of veterinary and medical importance. It is an important mechanical vector of several animal pathogens and can cause significant economic losses. However, the morphological variation of this species remains unknown. This study aimed to investigate the phenotypic variation in the wing size and shape of S. calcitrans populations in Thailand based on a landmark-based geometric morphometric approach. Specimens were collected from five populations in five geographical regions in Thailand. A total of 490 left wings of S. calcitrans (245 female and 245 male individuals) were used for geometric morphometric analysis. Wing size differences were detected between some populations of S. calcitrans, whereas wing shape differences were found among populations. Therefore, the phenotypic variation in S. calcitrans populations indicated that these populations are adaptive responses to local environmental pressures, suggesting the presence of phenotypic plasticity in this species.

Landmark-based geometric morphometric analysis of wing size and wing shape among Aedes albopictus (Skuse, 1894) populations in Turkey

Aedes albopictus (Skuse, 1894) has a widespread distribution that includes both temperate and tropical climates. It is a potential vector of several important worldwide arboviruses, including chikungunya, zika, and dengue, and its geographical distribution expands rapidly. Despite its vectorial importance, the bioecological features of this species in newly established environments are still not very well known. Thus, this study investigated phenotypic variations among Ae. albopictus populations from different regions in Turkey. This is the first comparative morphological study of Ae. albopictus populations in Turkey. The procrustes distances phenogram derived from wing data indicates that the shape differences among some populations and the population from the Aegean coast differ from the Black Sea populations. Size differences were also detected between some populations.

Geometric morphometric wing analysis represents a robust tool to identify female mosquitoes (Diptera: Culicidae) in Germany

Accurate species identification is the prerequisite to assess the relevance of mosquito specimens, but is often hindered by missing or damaged morphological features. The present study analyses the applicability of wing geometric morphometrics as a low-cost and practical alternative to identify native mosquitoes in Germany. Wing pictures were collected for 502 female mosquitoes of five genera and 19 species from 80 sampling sites. The reliable species identification based on interspecific wing geometry of 18 landmarks per specimen was tested. Leave-one-out cross validation revealed an overall accuracy of 99% for the genus and 90% for the species identification. Misidentifications were mainly due to three pairings of Aedes species: Aedes annulipes vs. Aedes cantans, Aedes cinereus vs. Aedes rossicus and Aedes communis vs. Aedes punctor. Cytochrome oxidase subunit I (COI) gene region was sequenced to validate the morphological and morphometric identification. Similar to the results of the morphometric analysis, the same problematic three Aedes-pairs clustered, but most other species could be well separated. Overall, our study underpins that morphometric wing analysis is a robust tool for reliable mosquito identification, which reach the accuracy of COI barcoding.