Geometric morphometric wing analysis as a tool to discriminate female mosquitoes from different suburban areas of Chiang Mai province, Thailand

Resolution of cryptic mosquito species through wing morphometrics

Mosquitoes are medically important insects, and accurate species identification is crucial to understanding vector biology, forming the cornerstone of successful vector control programs. Identification is difficult owing to morphologically similar species. Wing morphometrics can provide a simple, fast, and accurate way to classify species, and using it as a method to differentiate vector species among its cryptic congeners has been underexplored. Using a total of 227 mosquitoes and 20 landmarks per specimen, we demonstrated the utility of wing morphometrics in differentiating species two groups occurring in sympatry - Culex (Culex) vishnui group and Culex (Lophoceraomyia) subgenus, as well as explored population-level variation in the wing shape of Aedes albopictus across habitats. Cytochrome oxidase subunit I (COI) gene region was sequenced to validate the morphological and morphometric identification. Procrustes ANOVA regression and CVA based on wing shape reflected that the wing landmarks across all species differed significantly, and leave-one-out cross validation revealed an overall high accuracy of >97% for the two Culex groups. Wing morphometrics uncovered population-level variation within Aedes albopictus, but cross validation accuracy was low. Overall, we show that wing geomorphometric analysis is able to resolve cryptic Culex species (including vectors) occurring sympatrically, and is a robust tool for identifying mosquitoes reliably.

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.

Outline-based geometric morphometrics: Wing cell differences for mosquito vector classification in the Tanaosri mountain range, Thailand

Recent studies have revealed taxonomic signals within the wing cells of certain mosquito species. In our study, wing cell differentiation among mosquito vectors from the Tanaosri mountain range in Thailand was evaluated using the outline-based geometric morphometric (GM) approach. Our focus was on four specific wing cells for GM analysis: the wing contour (external cell), the second submarginal cell (internal cell 1), the first posterior cell (internal cell 2), and the third posterior cell (internal cell 3). Before proceeding with the GM approach, the identity of seven mosquito genera and 21 species was confirmed using molecular techniques. Our validated classification tests demonstrated that the performance of mosquito species classification varies according to genus. Notably, three Aedes species exhibited the highest accuracy for both internal cell 2 and internal cell 3, each registering a score of 93.20 %. In the case of two Mansonia species, the wing contour displayed a remarkable accuracy of 98.57 %. Consequently, we suggest the use of the outline-based GM approach, particularly focusing on the wing contour, for differentiating Mansonia annulifera and Mansonia uniformis. In contrast, the highest accuracy for classifying Culex species was found in internal cell 1, at 75.51 %, highlighting the challenges due to similarities in wing cells within this genus. These findings provide a guideline for future applications of the outline-based GM approach, focusing on wing cells, as an alternative method to classify mosquito vector species.

Effect of α-cypermethrin and pirimiphos-methyl on wing morphology of Tribolium castaneum (Herbst) and T. confusum Jacquelin du Val: a comparative study

Tribolium castaneum (Herbst) and Tribolium confusum Jacquelin du Val (Coleoptera: Tenebrionidae) are widespread and serious pests of stored products. Various insecticides are applied aiming to effectively manage both species. Here, two insecticides are tested, the pyrethroid α-cypermethrin and the organophosphate pirimiphos-methyl, hypothesizing that they can lead to morphological changes in the certain body parts of the adult offspring of treated T. castaneum and T. confusum parental female adults. For this purpose, the geometric morphometric method to the elytra and hindwings was applied. Both males and females were included in the analysis. The results showed that adult individuals of T. confusum showed higher tolerance to both insecticides compared to T. castaneum adults. This finding is reflected in analyses of both pairs of wings in T. confusum where changes in shape were negligible. The hindwings of T. castaneum experienced deformations to both insecticides. More significant changes in wing shape were observed in the α-cypermethrin treatment compared to pirimiphos-methyl. In the case of T. castaneum, even the shortest exposure to insecticides (5 min) is enough to provoke shape changes in the hindwings. Deformities in offspring, caused after treatment of their parents with insecticides, could moderate the frequency of insecticidal applications in storages.

A Phylogenetic Morphometric Investigation of Interspecific Relationships of Lyponia s. str. (Coleoptera, Lycidae) Based on Male Genitalia Shapes

The nominate subgenus Lyponia Waterhouse, 1878 from China is reviewed, with two new species described and named L. (s. str.) ruficeps sp. n. (China, Yunnan) and L. (s. str.) zayuana sp. n. (China, Xizang). A distribution map and a key to all species of Lyponia s. str. are provided. Moreover, the phenotypic relationships among the species of Lyponia s. str. are investigated based on phallus shapes using geometric morphometric and phylogenetic morphometric analyses. The topologies demonstrate that the species are divided into two clades. One clade is composed of six species (L. ruficeps sp. n., L. zayuana sp. n., L. kuatunensis, L. shaanxiensis, L. hainanensis, and L. tamdaoensis) and is supported by a stout phallus (less than 3.6 times longer than wide). The other clade includes the remaining species (L. nepalensis, L. debilis, L. cangshanica, L. delicatula, and L. oswai) and is supported by a slender phallus (at least 4.1 times longer than wide). These results provide better understanding of the species diversity and evolution of Lyponia s. str. Nonetheless, more samples and loci are required in the future to verify the present results.

Species Identification of the Major Japanese Encephalitis Vectors within the Culex vishnui Subgroup (Diptera: Culicidae) in Thailand Using Geometric Morphometrics and DNA Barcoding

Japanese encephalitis (JE) is a viral infection of the brain caused by the Japanese encephalitis virus, which spreads globally, particularly in 24 countries of Southeast Asia and the Western Pacific region. In Thailand, the primary vectors of JE are Cx. pseudovishnui, Cx. tritaeniorhynchus, and Cx. vishnui of the Cx. vishnui subgroup. The morphologies of three mosquito species are extremely similar, making identification challenging. Thus, geometric morphometrics (GM) and DNA barcoding were applied for species identification. The results of cross-validation reclassification revealed that the GM technique based on wing shape analysis had relatively high potential for distinguishing Cx. pseudovishnui, Cx. tritaeniorhynchus, and Cx. vishnui (total performance = 88.34% of correctly assigned individuals). While the DNA barcoding yielded excellent results in identifying these Culex species based on the DNA barcode gap (average intraspecific genetic distance = 0.78% ± 0.39% and average interspecific genetic distance = 6.14% ± 0.79%). However, in the absence of the required facilities for DNA barcoding, GM techniques can be employed in conjunction with morphological methods to enhance the reliability of species identification. Based on the results of this study, our approach can help guide efforts to identify members of the Cx. vishnui subgroup, which will be useful for the effective vector control of JE in Thailand.

Wing morphology variations in Culicoides circumscriptus from France

The biting midge Culicoides circumscriptus Kieffer, 1918 is a European widespread vector of avian malaria throughout the continent and is a possible vector of Akabane virus and Bluetongue virus. This species populates a wide range of environments in contrasting ecological settings often exposed to strong seasonal fluctuations. The main goals of this study were to investigate C. circumscriptus phenotypic variation at three departments in France (Corsica Island, Moselle and Var) and to determine if its phenotypes vary with the environment. Culicoides circumscriptus wing phenotypes were analyzed using a geometric morphometric approach based on anatomical landmarks and outlines of the wing. Dendogram trees based on landmarks and the outlines-2 set (cell m4) showed similar topologies and separated populations of C. circumscriptus. In contrast, another set of outlines-1 (covering the r-m cross vein, M, radiale and arculus) presented a different hierarchical clustering tree. The phenotypic variation observed in C. circumscriptus indicated that these populations are exposed to environmental and ecological pressures. Our results suggest the presence of phenotypic plasticity in this species.

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

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

Anopheles (Cellia) dirus Peyton & Harrison and Anopheles baimaii Sallum & Peyton are sibling species within the Dirus complex belonging to the Leucosphyrus group, and have been incriminated as primary vectors of malaria in Thailand. In the present study, DNA barcoding and geometric morphometrics were used to distinguish between An. dirus and An. baimaii in the international border areas, Trat Province, eastern Thailand. Our results revealed that DNA barcoding based on the cytochrome c oxidase subunit I gene could not be used to distinguish An. dirus from An. baimaii. The overlapping values between intra- and interspecific genetic divergence indicated no barcoding gap present for An. dirus and An. baimaii (ranging from 0 to 0.99%). However, the results of the geometric morphometric analysis based on the wing shape clearly distinguished An. dirus and An. baimaii, with 92.42% of specimens assigned to the correct species. We concluded that geometric morphometrics is an effective tool for the correct species identification of these two malaria vectors. Our findings could be used to make entomological surveillance information more accurate, leading to further effective mosquito control planning in Thailand and other countries in Southeast Asia.