Wing Geometric Morphometrics as a Tool for the Identification of Culex Subgenus Mosquitoes of Culex (Diptera: Culicidae)

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.

Evaluation of the wing cell contour to distinguish between Stomoxys bengalensis and Stomoxys sitiens (Diptera: Muscidae) using outline-based morphometrics

The blood-sucking flies of the genus Stomoxys Geoffroy, 1762 (Diptera: Muscidae) are significant ectoparasites that can cause irritation and transmit pathogens to both animals and humans. Within the genus Stomoxys, two species, Stomoxys bengalensis and Stomoxys sitiens, have similar morphology and coexist in the same habitat. Accurate species identification of these flies is crucial for understanding disease vectors and implementing effective control measures. In this study, we assessed the effectiveness of outline-based geometric morphometrics (GM) by analyzing the wing cell contour of the first posterior cell (R5) to distinguish between species and sexes of S. bengalensis and S. sitiens. Our results demonstrate that the outline-based GM method is highly effective in distinguishing between species and sexes of these flies based on contour shape, with accuracy scores ranging from 90.0% to 97.5%. Therefore, outline-based GM emerges as a promising alternative to landmark-based GM or as a supplementary tool in conjunction with traditional morphology-based methods for species identification.

Geometric Morphometrics and Genetic Diversity Analysis of Chalcidoidea (Diglyphus and Pachyneuron) at Various Elevations

Eulophidae and Pteromalidae are parasitic wasps with a global distribution and import for the biological control of pests. They can be distributed in different altitude regions, but their morphological and genetic adaptations to different altitudes are unclear. Here, we collected specimens that belong to Eulophidae and Pteromalidae from various altitudinal gradients, based on integrated taxonomic approaches to determine the species composition, and we analyzed their body shape and size from different altitudes using geometric morphometrics. Then, we performed an analysis of the D. isaea population's haplotype genes to illustrate their genetic diversity. As a result, eight species that belong to two genera, Diglyphus Walker (Eulophidae) and Pachyneuron Walker (Pteromalidae), were identified, including two newly recorded species from China (D. chabrias and D. sabulosus). Through a geometric morphometrics analysis of body shape, we found that a narrow forewing shape and a widened thorax are the significant characteristics of adaptation to high-altitude environments in D. isaea and P. aphidis. Additionally, the body size studies showed a principal relationship between centroid size and altitude; the size of the forewings and thorax increases at higher altitudes. Next, using haplotype analysis, 32 haplotypes were found in seven geographic populations with high genetic diversity of this species. Our research provides preliminary evidence for the morphological and genetic diversity adaptation of parasitic wasps to extreme environments, and these data can provide important references for investigations on the ecological adaptability of parasitic wasps.

Integrating artificial intelligence and wing geometric morphometry to automate mosquito classification

Mosquitoes (Diptera: Culicidae) comprise over 3500 global species, primarily in tropical regions, where the females act as disease vectors. Thus, identifying medically significant species is vital. In this context, Wing Geometric Morphometry (WGM) emerges as a precise and accessible method, excelling in species differentiation through mathematical approaches. Computational technologies and Artificial Intelligence (AI) promise to overcome WGM challenges, supporting mosquito identification. AI explores computers' thinking capacity, originating in the 1950s. Machine Learning (ML) arose in the 1980s as a subfield of AI, and deep Learning (DL) characterizes ML's subcategory, featuring hierarchical data processing layers. DL relies on data volume and layer adjustments. Over the past decade, AI demonstrated potential in mosquito identification. Various studies employed optical sensors, and Convolutional Neural Networks (CNNs) for mosquito identification, achieving average accuracy rates between 84 % and 93 %. Furthermore, larval Aedes identification reached accuracy rates of 92 % to 94 % using CNNs. DL models such as ResNet50 and VGG16 achieved up to 95 % accuracy in mosquito identification. Applying CNNs to georeference mosquito photos showed promising results. AI algorithms automated landmark detection in various insects' wings with repeatability rates exceeding 90 %. Companies have developed wing landmark detection algorithms, marking significant advancements in the field. In this review, we discuss how AI and WGM are being combined to identify mosquito species, offering benefits in monitoring and controlling mosquito populations.

Unraveling biotypes of the northern house mosquito, Culex pipiens s.l. (Diptera: Culicidae): molecular differentiation and morphometric analysis

Geometric morphometrics was used to determine whether geographic isolation could explain differences in wing size and shape between and within continental (27°S to 41°S) and insular (Rapa Nui) populations of Culex pipiens s.s. Linnaeus and their biotypes (f. pipiens and f. molestus). Molecular protocols based on polymorphisms in the second intron of nuclear locus ace-2 (acetylcholinesterase-2) were used to differentiate Cx. pipiens s.s. from Cx. quinquefasciatus Say, and an assay based on polymorphisms in the flanking region of a microsatellite locus (CQ11) was used to identify biotypes. Culex pipiens f. molestus and hybrids shared larval habitats in all continental sites, while Cx. pipiens f. pipiens was found in 5 of the 10 sites. Only biotype molestus was found in Rapa Nui (Easter Island) Pipiens and molestus biotypes occur sympatrically in aboveground locations, and only molestus was found in the underground site (ME). Biotype molestus was dominant in rural locations and preferably anthropophilic. These results agree with the ecological descriptions previously reported for the biotypes of Cx. pipiens s.s. Procrustes ANOVA only showed differences in centroid size between biotypes in females and males and did not show significant differences in wing shape. However, we found significant differences among the geographic areas in the centroid size and wing shape of both females and males. Particularly, the population of Rapa Nui Island had shorter wings than the continental populations. The results highlight the effects of geographic and environmental processes on morphotypes in vector mosquitoes.

Using geometric wing morphometrics to distinguish Aedes japonicus japonicus and Aedes koreicus

BACKGROUND

Aedes japonicus japonicus (Theobald, 1901) and Aedes koreicus (Edwards, 1917) have rapidly spread in Europe over the last decades. Both species are very closely related and occur in sympatry. Females and males are difficult to distinguish. However, the accurate species discrimination is important as both species may differ in their vectorial capacity and spreading behaviour. In this study, we assessed the potential of geometric wing morphometrics as alternative to distinguish the two species.

METHODS

A total of 147 Ae. j. japonicus specimens (77 females and 70 males) and 124 Ae. koreicus specimens (67 females and 57 males) were collected in southwest Germany. The left wing of each specimen was removed, mounted and photographed. The coordinates of 18 landmarks on the vein crosses were digitalised by a single observer. The resulting two-dimensional dataset was used to analyse the differences in the wing size (i.e. centroid size) and wing shape between Ae. j. japonicus and Ae. koreicus using geometric morphometrics. To analyse the reproducibility of the analysis, the landmark collection was repeated for 20 specimens per sex and species by two additional observers.

RESULTS

The wing size in female Ae. koreicus was significantly greater than in Ae. j. japonicus but did not differ significantly for males. However, the strong overlap in wing size also for the females would not allow to discriminate the two species. In contrast, the wing shape clustering was species specific and a leave-one-out validation resulted in a reclassification accuracy of 96.5% for the females and 91.3% for the males. The data collected by different observers resulted in a similar accuracy, indicating a low observer bias for the landmark collection.

CONCLUSIONS

Geometric wing morphometrics provide a reliable and robust tool to distinguish female and male specimens of Ae. j. japonicus and Ae. koreicus.

Wing geometric morphometrics to distinguish and identify Haematobosca flies (Diptera: Muscidae) from Thailand

The hematophagous flies of the genus Haematobosca Bezzi, 1907 (Diptera: Muscidae) are important ectoparasites in domestic animals and wildlife. Two species of this genus have been recorded in Thailand, viz., Haematobosca sanguinolenta (Austen, 1909) and Haematobosca aberrans (Pont, Duvallet & Changbunjong, 2020). They have a similar morphology and coexist in the same habitat. The correct species identification of these flies is extremely important for understanding disease epidemiology and developing effective control measures. Geometric morphometrics (GM) has been confirmed to be a useful tool for differentiating and identifying morphologically similar insect species. Therefore, GM was used to distinguish and identify H. sanguinolenta and H. aberrans in Thailand. Adult flies of both sexes were collected using Nzi traps, morphologically identified, and analyzed by landmark-based GM of the wing. Results showed that GM was highly effective in distinguishing the two Haematobosca species based on their wing shape, with an overall accuracy score of 99.3%. We also revealed that our study material could be used as reference data to identify new field specimens collected from other geographic locations. We propose that wing GM can be used as a supplement to conventional morphology identification, particularly for Haematobosca specimen that has been damaged or has lost its diagnostic characteristics due to specimen collection and processing in the field.

Detection of Hindwing Landmarks Using Transfer Learning and High-Resolution Networks

Hindwing venation is one of the most important morphological features for the functional and evolutionary analysis of beetles, as it is one of the key features used for the analysis of beetle flight performance and the design of beetle-like flapping wing micro aerial vehicles. However, manual landmark annotation for hindwing morphological analysis is a time-consuming process hindering the development of wing morphology research. In this paper, we present a novel approach for the detection of landmarks on the hindwings of leaf beetles (Coleoptera, Chrysomelidae) using a limited number of samples. The proposed method entails the transfer of a pre-existing model, trained on a large natural image dataset, to the specific domain of leaf beetle hindwings. This is achieved by using a deep high-resolution network as the backbone. The low-stage network parameters are frozen, while the high-stage parameters are re-trained to construct a leaf beetle hindwing landmark detection model. A leaf beetle hindwing landmark dataset was constructed, and the network was trained on varying numbers of randomly selected hindwing samples. The results demonstrate that the average detection normalized mean error for specific landmarks of leaf beetle hindwings (100 samples) remains below 0.02 and only reached 0.045 when using a mere three samples for training. Comparative analyses reveal that the proposed approach out-performs a prevalently used method (i.e., a deep residual network). This study showcases the practicability of employing natural images-specifically, those in ImageNet-for the purpose of pre-training leaf beetle hindwing landmark detection models in particular, providing a promising approach for insect wing venation digitization.

Species Discrimination of Stomoxys Flies S. bengalensis, S. calcitrans, and S. sitiens (Diptera: Muscidae) Using Wing Geometric Morphometrics

The flies of the genus Stomoxys Geoffroy, 1762 (Diptera: Muscidae), are regarded as pests of veterinary and medical importance. In Thailand, Stomoxys calcitrans (Linnaeus, 1758) is the most abundant species and is widely distributed throughout the country. This Stomoxys species can coexist with two other morphologically similar species: Stomoxys bengalensis Picard, 1908, and Stomoxys sitiens Rondani, 1873. Hence, discriminating using morphological characteristics is difficult, especially if the specimen is damaged or loses its diagnostic characteristics. In this study, we evaluated the effectiveness of the landmark-based geometric morphometric (GM) approach to discriminate among the three Stomoxys spp.: S. bengalensis, S. calcitrans, and S. sitiens. Left-wing images of S. bengalensis (n = 120), S. calcitrans (n = 150), and S. sitiens (n = 155) were used for the GM analyses. The results of the wing shape analyses revealed that the GM approach was highly effective for discriminating three Stomoxys, with high accuracy scores ranging from 93.75% to 100%. This study adds to the evidence that landmark-based GM is an excellent alternative approach for discriminating Stomoxys species.

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.

Evaluation of Modern Techniques for Species Identification of Lutzia Mosquitoes (Diptera: Culicidae) in Thailand: Geometric Morphometrics and DNA Barcoding

There are four species of Lutzia mosquitoes in Thailand, including Lutzia chiangmaiensis, Lt. fuscana, Lt. halifaxii, and Lt. vorax. The accurate species identification of adult Lutzia mosquitoes based on morphological features requires many body parts, including the abdominal terga and wing. However, species identification is difficult in the case of damaged specimens when some of their morphological character is missing due to transit or gathering in the field. Thus, we evaluated the efficacy of the landmark-based geometric morphometric (GM) approach for the discrimination of Lutzia species in Thailand. In addition, DNA barcoding was also used in parallel with the GM approach to identify the species. Larvae of Lutzia were collected, raised into adults, and identified based on their morphological characteristics. The validated reclassification test results clearly demonstrated that wing shape resulted in a high level of success in identification (correct identifications ranged from 92.50% to 100%); however, based on the DNA barcoding analyses, our results showed that it was poorly effective in identifying Lt. fuscana and Lt. halifaxii based on an overlap between the intraspecific and interspecific divergence. Moreover, our survey results provide updates on the distribution of Lt. chiangmaiensis and Lt. vorax in Thailand. This research will help medical entomologists more efficiently identify mosquitoes in the genus Lutzia, resulting in more effective mosquito control and surveillance.

Review and redescription of species in the brasiliana group of Smicridea (Rhyacophylax) (Trichoptera, Hydropsychidae, Smicrideinae): exploration of the utility of geometric morphometrics as a method for delimitation and characterization of species in the genus

The Smicrideabrasiliana species group includes five species distributed in northeastern Argentina and Brazil: Smicridea (Rhyacophylax) brasiliana (Ulmer), S. (R.) weidneri Flint, S. (R.) vermiculata Flint, S. (R.) arobasis Flint, and S. (R.) nanda Flint. The original descriptions of these species and their placement in the brasiliana species group were mainly based on the morphology of the male genitalia. However, the fine structure of the internal sclerites of the phallus, which proved to be useful for species delimitation, was not analyzed at the time. In this contribution, we provide a detailed description of the male genitalia and the morphology of the head, and analyze the shape of the wings using geometric morphometrics. The analyzed species can be easily differentiated by the shape of the phallus, especially by the structure of the internal sclerites, the shape of the head in dorsal view, and the shape of the cephalic setose warts. Furthermore, the geometric morphometric approach allowed their separation through the wing shape. The preliminary analysis of these features suggests that the brasiliana species group is not natural but its monophyly should be further tested within the framework of a phylogenetic analysis of all the species of the subgenus Rhyacophylax.

The Influence of the pH and Salinity of Water in Breeding Sites on the Occurrence and Community Composition of Immature Mosquitoes in the Green Belt of the City of São Paulo, Brazil

Simple Summary

Immature mosquitoes are found in natural and artificial aquatic habitats. Variations in physicochemical parameters of water, such as pH, salinity, conductivity, and total dissolved solids, in breeding habitats can influence larval occurrence and drive the proliferation of adult mosquitoes. Herein, we investigated the association between different values of physicochemical parameters in a variety of aquatic habitats and the occurrence and community composition of immature mosquito species in two environmentally protected areas in the city of São Paulo, Brazil. The aquatic habitats surveyed included epiphytic and ground bromeliads, bamboo internodes, ponds, tree hollows, lakes, and artificial containers. Our results revealed a statistically significant relationship between species occurrence and the variables of pH and salinity. The type of aquatic habitat also had a significant influence on mosquito species distribution. Investigating the interactions between immature mosquitoes and the environment in which they develop is important to elucidate the factors driving their occurrence and abundance, and could also be an important tool in planning and implementing immature mosquito control practices.

Abstract

The physicochemical parameters of water, such as pH, salinity, conductivity, and total dissolved solids, can influence mosquito larval development, survival, and abundance. Therefore, it is important to elucidate how these factors influence mosquito occurrence. We hypothesized that the occurrence and community composition of immature mosquito species are driven not only by the availability of suitable aquatic habitats, but also by the physicochemical factors of these habitats. The primary objective of this study was therefore to investigate the influence of the physicochemical parameters of water in different types of aquatic habitats on the occurrence of mosquito species in two remnants of Atlantic Forest in the city of São Paulo, Brazil. Collections of immature mosquitoes and assessment of the physicochemical characteristics of the water in the collection sites were carried out for twelve months. The variation in species composition and occurrence with the measured physicochemical parameters and the type of breeding site was assessed using constrained ordination methods. The results indicate that there was a statistically significant difference in species composition as a function of the different types of aquatic habitats, and that pH had an influence on species occurrence even when the variance explained by the type of aquatic habitat was removed from the analysis. There was a statistically significant association between mosquito species occurrence and pH and salinity, and the former had a significant influence on the mosquito species collected regardless of the type of aquatic habitat, showing that the pH of the breeding site water is an important factor in driving mosquito population dynamics and species distribution.