Geometric morphometrics in mosquitoes: What has been measured?

Robust mosquito species identification from diverse body and wing images using deep learning

Mosquito-borne diseases are a major global health threat. Traditional morphological or molecular methods for identifying mosquito species often require specialized expertise or expensive laboratory equipment. The use of convolutional neural networks (CNNs) to identify mosquito species based on images may offer a promising alternative, but their practical implementation often remains limited. This study explores the applicability of CNNs in classifying mosquito species. It compares the efficacy of body and wing depictions across three image collection methods: a smartphone, macro-lens attached to a smartphone and a professional stereomicroscope. The study included 796 specimens of four morphologically similar Aedes species, Aedes aegypti, Ae. albopictus, Ae. koreicus and Ae. japonicus japonicus. The findings of this study indicate that CNN models demonstrate superior performance in wing-based classification 87.6% (95% CI: 84.2-91.0) compared to body-based classification 78.9% (95% CI: 77.7-80.0). Nevertheless, there are notable limitations of CNNs as they perform reliably across multiple devices only when trained specifically on those devices, resulting in an average decline of mean accuracy by 14%, even with extensive image augmentation. Additionally, we also estimate the required training data volume for effective classification, noting a reduced requirement for wing-based classification compared to body-based methods. Our study underscores the viability of both body and wing classification methods for mosquito species identification while emphasizing the need to address practical constraints in developing accessible classification systems.

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.

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.

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.

Demographic inference from the mt-DNA COI gene and wing geometry of Culex gelidus (Diptera: Culicidae), an important vector of Japanese encephalitis in Thailand

Culex gelidus (Diptera: Culicidae), an important vector of the Japanese encephalitis virus (JEV), contributes to human viral encephalitis in many Asian countries, including Thailand. This study represents the first investigation of the demographic patterns of Cx. gelidus populations in Thailand using cytochrome c oxidase subunit I (COI) gene analysis and wing geometric morphometrics (GM). Mosquitoes were collected from 10 provinces across six regions of Thailand in 2022. Analysis of the COI sequences (n = 182) indicated high haplotype diversity (0.882) and low nucleotide diversity (0.006), with 72 haplotypes identified. The haplotype network demonstrated no profound splits among the geographic populations. Neutral tests, including Tajima's D and Fu's Fs, displayed negative values, with a significant result observed for Fu's Fs (-33.048, p < 0.05). The mismatch distribution analysis indicated that the population does not statistically deviate from a model of sudden population expansion (SSD = 0.010, p > 0.05; Rg = 0.022, p > 0.05). The estimations suggest that the Cx. gelidus population in Thailand began its expansion approximately between 459,243 and 707,011 years ago. The Mantel test showed no significant relationship between genetic and geographic distances (r = 0.048, p > 0.05). Significant phenotypic differences (based on wing shape) were observed among most populations. Additionally, in this study, we found no significant relationships between phenotypic and genetic distances (r = 0.250, p > 0.05). Understanding the genetic and morphological dynamics of Cx. gelidus is vital for developing targeted surveillance and vector control measures. This knowledge will also help to predict how future environmental changes might affect these populations, thereby informing long-term vector management strategies.

Morphometry of the wings of Anopheles aquasalis in simulated scenarios of climate change

BACKGROUND

Climate change has significant implications on ecosystems. We verified the effects of climate change on the malaria vector Anopheles aquasalis using simulated climate change scenarios (SSCCs).

METHODS

An experimental model was designed for SSCCs, which composed of air-conditioned 25 m3 rooms.

RESULTS

The wing size was significantly different between SSCCs. A colony of Anopheles aquasalis could not be established in extreme scenarios.

CONCLUSIONS

Increases in temperature and CO2 in the atmosphere may modify the global epidemiology of malaria, marking its emergence in currently malaria-free areas.

Wolbachia strains wMel and wAlbB differentially affect Aedes aegypti traits related to fecundity

Two Wolbachia strains, wMel and wAlbB, have been transinfected into Aedes aegypti mosquitoes for population replacement with the aim of reducing dengue transmission. Epidemiological data from various endemic sites suggest a pronounced decrease in dengue transmission after implementing this strategy. In this study, we investigated the impact of the Wolbachia strains wMel and wAlbB on Ae. aegypti fitness in a common genetic background. We found that Ae. aegypti females infected with the wMel strain exhibited several significant differences compared with those infected with the wAlbB strain. Specifically, wMel-infected females laid significantly fewer eggs, ingested a lower amount of blood, had a reduced egg production rate, and exhibited a decreased Wolbachia density at a later age compared with mosquitoes infected with the wAlbB strain. Conversely, the wAlbB strain showed only mild negative effects when compared with Wolbachia-uninfected specimens. These differential effects on Ae. aegypti fitness following infection with either wMel or wAlbB may have important implications for the success of population replacement strategies in invading native Ae. aegypti populations in endemic settings. Further research is needed to better understand the underlying mechanisms responsible for these differences in fitness effects and their potential impact on the long-term efficacy of Wolbachia-based dengue control programs.IMPORTANCEThe transmission of arboviruses such as dengue, Zika, and chikungunya is on the rise globally. Among the most promising strategies to reduce arbovirus burden is the release of one out of two strains of Wolbachia-infected Aedes aegypti: wMel and wAlbB. One critical aspect of whether this approach will succeed involves the fitness cost of either Wolbachia strains on mosquito life history traits. For instance, we found that wMel-infected Ae. aegypti females laid significantly fewer eggs, ingested a lower amount of blood, had a reduced egg production rate, and exhibited a decreased Wolbachia density at a later age compared with mosquitoes infected with the wAlbB strain. Conversely, the wAlbB strain showed only mild negative effects when compared with Wolbachia-uninfected specimens. These differential effects on mosquito fitness following infection with either wMel or wAlbB may have important implications for the success of population replacement strategies in invading native Ae. aegypti populations.

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.

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.

Wing geometric morphometrics and COI barcoding of Culex pipiens subgroup in the Republic of Korea

Two members of the Culex pipiens subgroup, Culex pallens and Culex pipiens f. molestus, are known to occur in the Republic of Korea (ROK). These species exhibit morphologically similar features and are challenging to distinguish below the species level. Therefore, this study utilized wing geometric morphometrics (GM) on the right wing of the Culex pipiens subgroup, alongside sequencing of the cytochrome c oxidase subunit I (COI) region. Mosquitoes were collected from 11 locations between June and October (2020-2022) to minimize regional and seasonal variations. Additionally, Culex pipiens f. pipiens, which is not native to the ROK, was included in the analysis. Culex tritaeniorhynchus, Aedes albopictus, and Anopheles sinensis, the primary vectors in the ROK, were used as outgroups for comparison. All three taxa in the Culex pipiens subgroup could be identified with an 82.4%-97.0% accuracy using GM. However, a comparison of the COI regions of the Culex pipiens subgroup revealed no clear differences between the taxa. These data can be used for accurate identification, contributing to effective mosquito control, in addition to providing a foundation for evolutionary and ecological studies on wing shape differences.

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.

Wing geometric morphometrics and DNA barcoding to distinguish three closely related species of Armigeres mosquitoes (Diptera: Culicidae) in Thailand

Armigeres subalbatus, a mosquito species widely found in Thailand and other Asian countries, serves as a vector for filarial parasites, affecting both humans and animals. However, the surveillance of this vector is complicated because of its morphological similarity to two other species, Armigeres dohami and Armigeres kesseli. To differentiate these morphologically similar species, our study employed both wing geometric morphometrics (GM) and DNA barcoding, offering a comprehensive approach to accurately identify these closely related Armigeres species in Thailand. Our GM analyses based on shape demonstrated significant accuracy in differentiating Armigeres species. Specifically, the outline-based GM method focusing on the 3rd posterior cell exhibited an accuracy rate of 82.61%, closely followed by the landmark-based GM method with 81.54%. Both these GM techniques effectively distinguished Ar. subalbatus from Ar. dohami and Ar. kesseli. Regarding DNA barcoding, our investigation of pairwise intra- and interspecific divergences revealed a "barcoding gap". Furthermore, the results of species confirmation using both species delimitation methods including the automatic barcode gap discovery method (ABGD) and the Multi-rate Poisson tree process (mPTP) were consistent with those of morphological identification, sequence comparisons with the GenBank and Barcode of Life Data System (BOLD) databases, and the neighbor-joining tree construction. These consistent results emphasize the efficacy of DNA barcoding in the precise identification of Armigeres species.

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.

Analysis of geometric morphometrics and molecular phylogeny for Anopheles species in the Republic of Korea

Human malaria, transmitted by Anopheles mosquitoes, is the most predominant mosquito-borne disease that is responsible for hundreds of thousands of deaths worldwide each year. In the Republic of Korea (ROK), there are currently several hundred malaria cases annually, mostly near the demilitarized zone (DMZ). Eight species of Anopheles mosquitoes are currently known to be present in the ROK. Similar to other major malaria vectors in Africa and India, it is very challenging to morphologically differentiate Anopheles mosquitoes in the ROK due to their extremely similar morphology. In this study, wing geometric morphometrics (WGM) were used to differentiate the eight Anopheles species collected at six locations near the DMZ, Seoul and Pyeongtaek from April-October 2021. Phylogenetic analysis was also performed using cytochrome c oxidase subunit 1 (COI), internal transcribed spacer 2 (ITS2), and tyrosine hydroxylase (TH) genes for comparison with WGM analysis and to infer evolutionary relationships. The results of cross-validation (overall accuracy = 74.8%) demonstrated that species identification using WGM alone was not possible with a high accuracy for all eight species. While phylogenetic analyses based on the COI region could not clearly distinguish some species, the analysis based on ITS2 and TH was more useful for resolving the phylogenetic correlation of the eight species. Our results may improve Anopheles species identification strategies for effective identification and control of malaria vectors in the ROK.

Wing morphometrics of biting midges (Diptera: Culicoides) of veterinary importance in Madagascar

Biting midges are vectors of arboviruses such as bluetongue virus, bovine ephemeral fever virus, Akabane virus, African horse sickness virus, epizootic haemorrhagic disease virus and Schmallenberg virus. Fast and accurate identification of biting midges is crucial in the study of Culicoides-borne diseases. Morphological identification of biting midges has revealed the presence of cryptic species. A total of 20 species are reported in Madagascar. In this study, we assessed wing morphometric analysis for identification of seven species namely C. dubitatus Kremer, Rebholtz-Hirtzel and Delécolle, C. enderleini Cornet and Brunhes, C. kibatiensis Goetghebuer, C. miombo Meiswinkel, C. moreli Clastrier, C. nevilli Cornet and Brunhes, and C. zuluensis de Meillon. Culicoides enderleini, C. miombo, C. moreli, C. nevilli and C. zuluensis are vectors diseases. A molecular approach, based on the cytochrome oxidase I gene (Cox1), was used for species delimitation. The molecular analysis presented seven different clades grouped two-by-two according to morphological characters. A total of 179 wing images were digitised. We found morphometric variation among seven species based on 11 landmarks and two outlines. Wing shape variation plots showed that species overlapped with species belonging to the same group. The cross-validation revealed a relatively high percentage of correct classification in most species, ranging from 91.3% to 100% for landmarks; 60% to 82.6% for outlines-1 and 77.1% to 91.3% for outlines-2. Our study suggests that wing geometric morphometric analysis is a robust tool for reliable "Moka Fohy" identification in Madagascar. This inexpensive and simple method is a precise supplement to morphological identification, with reaches the accuracy of Cox1 barcoding.

Genetic diversity and wing geometric morphometrics among four populations of Aedes aegypti (Diptera: Culicidae) from Benin

BACKGROUND

The impact of the arbovirus vector Aedes aegypti is of major concern for global public health as the viruses that it transmits affect millions of people each year worldwide. Originating in Africa, Ae. aegypti has now spread throughout much of the world. While the genetic makeup of Ae. aegypti in the New World has been extensively studied, there is limited knowledge on its genetic diversity in Africa, particularly at a microgeographical level.

METHODS

We investigated mitochondrial cytochrome oxidase I of four Ae. aegypti populations from Benin and employed wing morphometric analyses as a cost-effective and reliable tool to explore population structure. Our sampling encompassed various areas of Benin, from the southern to the northern borders of the country, and included urban, semi-urban, and sylvatic sites.

RESULTS

We observed a notable level of genetic diversity (haplotype diversity of 0.8333) and nucleotide diversity (0.00421986), and identified seven distinct haplotypes. Sylvatic and semi-urban sites exhibited a greater number of haplotypes compared to urban sites. Utilizing 18 wing landmarks, we calculated the centroid size, which revealed significant variation among the three landscape types. However, principal component analysis, employed to assess wing shape variation, did not demonstrate significant differences between populations based on landscape type.

CONCLUSIONS

Our findings indicate substantial genetic and morphological diversity among Ae. aegypti populations in Benin, and provide insight into important biological characteristics of these populations with respect to their potential to transmit viruses. To the best of our knowledge, this is the first study undertaken in Africa to integrate genetics with morphology to analyse the population structure of the major arbovirus vector Ae. aegypti.

Disparification and extinction trade-offs shaped the evolution of Permian to Jurassic Odonata

Owing to their prevalence in nowadays terrestrial ecosystems, insects are a relevant group to assess the impact of mass extinctions on emerged land. However, limitations of the insect fossil record make it difficult to assess the impact of such events based on taxonomic diversity alone. Therefore, we documented trends in morphological diversity, i.e., disparity, using wings of Permian to Jurassic Odonata as model. Our results show a decreasing trend in disparity while species richness increased. Both the Permian-Triassic and Triassic-Jurassic transitions are revealed as important events, associated with strong morphospace restructuring due to selective extinction. In each case, a recovery was assured by the diversification of new forms compensating the loss of others. Early representatives of Odonata continuously evolved new shapes, a pattern contrasting with the classical assertion of a morphospace fulfilled early and followed by selective extinctions and specialization within it.

When size and shape matter: morphometric characterization of two sympatric dragonflies of the genus Perithemis (Odonata: Libellulidae)

Perithemis mooma, Kirby, 1889 and Perithemis icteroptera (Selys in Sagra, 1857) live in sympatry from southern Brazil to central Argentina. The taxonomy of the genus Perithemis Hagen, 1861 has been hampered by the use of characters that are highly variable or show slight differences among species. Our objective was to assess the efficiency of traditional morphometrics (TM) and geometric morphometrics (GM) to discriminate between these species using wing size and shape and vulvar lamina contour, and to analyze the presence of sexual dimorphism in wing size and shape in both species. The TM and landmark-GM methods were applied on the fore and hind wings, while the outline-based GM method was applied on the vulvar lamina. GM allowed species delimitation using shape variables of either wing. The wing and vulvar lamina shapes were confirmed to be good diagnostic characters to separate these species and appear to be promising tools for distinguishing among other species of this genus. Centroid size failed to achieve species separation. Both species exhibited sexual size dimorphism (SSD). In contrast to what would be expected for Perithemis whose males are strongly territorial, P. icteroptera and P. mooma showed female-biased SSD suggesting a common pattern in Perithemis.

Assessment of wing geometric morphometrics of urban Culex quinquefasciatus (Diptera: Culicidae) populations

Culex quinquefasciatus is a cosmopolitan species distributed throughout tropical and subtropical areas of the world. The species is of great epidemiological importance as it is responsible for vectoring the causative agent of lymphatic filariasis and several arboviruses, including West Nile virus. Wing geometric morphometrics has been widely used to assess phenotypic variations in mosquito species. Here, we hypothesize that Cx. quinquefasciatus populations in urban parks in the city of São Paulo, Brazil, have been subjected to anthropogenic selective pressures that are responsible for driving their ecology and behavior. Mosquitoes were collected by CDC traps in five municipal parks in the city of São Paulo. Eighteen anatomical landmark coordinates on each female right wing were digitized. Canonical variate analysis, wireframe graphs, cross-validated reclassification tests and the neighbor-joining method were used to assess phenotypical dissimilarity in wing shape between populations. Centroid size was calculated to assess differences in wing size between populations, which can result from different environmental conditions during immature mosquito development. Moderately heterogeneous wing shape and wing size patterns were found in the populations analyzed, indicating that selective pressures in the urban environment are affecting the wing patterns of Cx. quinquefasciatus populations in the city of São Paulo, Brazil.

Molecular and morphometric differentiation of secondary filariasis vector Coquillettidia mosquitoes (Diptera: Culicidae) in Thailand

Coquillettidia mosquitoes are important nuisance-biting pests and a vector of brugian filariasis in Thailand. However, comprehensive information about these mosquitoes remains unavailable such as molecular and morphometric differences among species. The lack of vector knowledge on Coquillettidia species could affect future disease control. This study aims to investigate differences in molecular variations based on mitochondrial cytochrome oxidase subunit I (COI) gene and wing geometric traits of three Coquillettidia species, namely Cq. crassipes, Cq. nigrosignata, and Cq. ochracea in Thailand. The results of molecular analyses revealed the differences among three Coquillettidia species. The genetic difference measure based on the Kimura two-parameter model among three Coquillettidia species showed low intraspecific distances (0%-3.05%) and large interspecific distances (10.10%-12.41%). The values of intra- and inter-genetic differences of three Coquillettidia species did not overlap which showed the existence of a barcoding gap indicating the efficiency of the identification based on the COI gene. As with molecular analysis, the landmark-based geometric morphometrics approach based on wing shape analysis indicated three distinct species groups which were supported by the high total performance score of cross-validated classification (97.16%). These results provide the first evidence of taxonomic signal based on molecular and wing geometric differences to support species identification and biological variations of Coquillettidia mosquitoes in Thailand for understanding these rare vector mosquitoes in depth and leading to effective further mosquito control.

Sub-Lethal Effects of Pirimiphos-Methyl Are Expressed to Different Levels in Wings of Three Stored-Product Coleopterans: A Geometric Morphometrics Investigation

Tenebrio molitor L. (Coleoptera: Tenebrionidae), Prostephanus truncatus (Horn), and Rhyzopertha dominica (F.) (Coleoptera: Bostrychidae) are noxious insect pests of grains in storages. Pirimiphos-methyl is widely used to protect grains at the post-harvest stage. However, the sub-lethal impact of this active ingredient on the offspring of all three coleopterans remains unknown. Thus, mated females of each species were exposed separately to pirimiphos-methyl at short exposures (30 min, 3, 5, 8, 16, 24, and 36 h), where the elytra and hindwings of the adult progeny were analyzed with the geometric morphometrics method. Males and females of all species were incorporated into the analysis. The results revealed variability among species. Tenebrio molitor was the most sensitive among three species, displaying significant deformations in the elytra and hindwings. Males had more conspicuous morphological changes than females. Prostephanus truncatus hindwings exhibited deformities after 36 h of exposure to pirimiphos-methyl. In contrast, R. dominica offspring were not affected by pirimiphos-methyl. In light of our findings, organophosphorus insecticides may cause variable sub-lethal effects to stored-product insects. This issue may lead to different insecticidal treatments according to the targeted stored-product 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.

Quantitative analysis of the labia minora morphology in 400 Chinese women: A new method for assessing the shape of the labia minora

Purpose

To investigate the shape of labia minora in Chinese adult women.

Methods

Women who visited the Genital Plastic Surgery Center from January 2021 to February 2022 were included and the vulvar regions were photographed. The shape of the labia minora was converted into quantifiable values [left side line segment 1-9 (L1-L9 values) and right side line segment 1-9 (R1-R9 values)]. The 400 enrolled patients were grouped by age, parity, and the types of vulvar surgeries they planned to undergo.

Results

After a graphic-to-digital and digital-to-graphic conversion, the simulated average image of the labia minora in Chinese adult women was constructed based on the mean values of L1-L9 and R1-R9. Comparing the values of Ln and Rn in the four age groups revealed that the mean values of Ln and Rn gradually decreased with age, but only two subgroups showed statistically significant differences (P-value <0.05). When the patients were stratified according to the number of births and whether labiaplasty was planned to perform, there were no significant differences among all groups (P-value >0.05).

Conclusion

To measure and evaluate the labia minora, a new assessing method was used. It is an innovative attempt to transform the simple rough description of the shape of the labia minora into more precise data reflecting the shape features. It allows the comparison of labia minora shape between individuals by comparing quantifiable values. Besides, the average shape of labia minora in Chinese adult women was provided visually. Factors including age, parity, and whether the patient planned to undergo labiaplasty were not significantly associated with the shape of the labia minora.

Automatic acoustic recognition of pollinating bee species can be highly improved by Deep Learning models accompanied by pre-training and strong data augmentation

Introduction

Bees capable of performing floral sonication (or buzz-pollination) are among the most effective pollinators of blueberries. However, the quality of pollination provided varies greatly among species visiting the flowers. Consequently, the correct identification of flower visitors becomes indispensable to distinguishing the most efficient pollinators of blueberry. However, taxonomic identification normally depends on microscopic characteristics and the active participation of experts in the decision-making process. Moreover, the many species of bees (20,507 worldwide) and other insects are a challenge for a decreasing number of insect taxonomists. To overcome the limitations of traditional taxonomy, automatic classification systems of insects based on Machine-Learning (ML) have been raised for detecting and distinguishing a wide variety of bioacoustic signals, including bee buzzing sounds. Despite that, classical ML algorithms fed by spectrogram-type data only reached marginal performance for bee ID recognition. On the other hand, emerging systems from Deep Learning (DL), especially Convolutional Neural Networks (CNNs), have provided a substantial boost to classification performance in other audio domains, but have yet to be tested for acoustic bee species recognition tasks. Therefore, we aimed to automatically identify blueberry pollinating bee species based on characteristics of their buzzing sounds using DL algorithms.

Methods

We designed CNN models combined with Log Mel-Spectrogram representations and strong data augmentation and compared their performance at recognizing blueberry pollinating bee species with the current state-of-the-art models for automatic recognition of bee species.

Results and Discussion

We found that CNN models performed better at assigning bee buzzing sounds to their respective taxa than expected by chance. However, CNN models were highly dependent on acoustic data pre-training and data augmentation to outperform classical ML classifiers in recognizing bee buzzing sounds. Under these conditions, the CNN models could lead to automating the taxonomic recognition of flower-visiting bees of blueberry crops. However, there is still room to improve the performance of CNN models by focusing on recording samples for poorly represented bee species. Automatic acoustic recognition associated with the degree of efficiency of a bee species to pollinate a particular crop would result in a comprehensive and powerful tool for recognizing those that best pollinate and increase fruit yields.

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 Wing Analysis of Avian Malaria Vector, Culiseta longiareolata, from Two Locations in Algeria

The application of geometric morphometry on mosquito wings (Culicidae) is considered a powerful tool for evaluating correlations between the phenotype (e.g., shape) and environmental or genetic variables. However, this has not been used to study the wings of the avian malaria vector, Culiseta longiareolata. Therefore, the goal of this study is to investigate the intra-specific wing variations between male and female Cs. longiareolata populations in different types of larval habitats and climatic conditions in Algeria. A total of 256 Cs. longiareolata mosquito samples were collected from January 2020 to July 2021 in three cities (Annaba, El-Tarf, and Guelma) of northeastern Algeria that have two distinct climatic condition levels (sub-humid and sub-arid) and different types of larval habitats (artificial and natural). Nineteen (19) wing landmarks (LMs) were digitized and analyzed based on geometric morphometry. Our results revealed differences in the wing shape of female and male mosquito populations, indicating sexual dimorphism. Moreover, canonical variance analysis (CVA) showed that factors, such as climatic conditions and type of larval habitats, also affect the wing shape of female and male Cs. longiareolata mosquito populations. Furthermore, the wing shape of male populations was more distinct compared with female populations.

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

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.

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.

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.

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.

Strain Characterisation for Measuring Bioefficacy of ITNs Treated with Two Active Ingredients (Dual-AI ITNs): Developing a Robust Protocol by Building Consensus

Durability monitoring of insecticide-treated nets (ITNs) containing a pyrethroid in combination with a second active ingredient (AI) must be adapted so that the insecticidal bioefficacy of each AI can be monitored independently. An effective way to do this is to measure rapid knock down of a pyrethroid-susceptible strain of mosquitoes to assess the bioefficacy of the pyrethroid component and to use a pyrethroid-resistant strain to measure the bioefficacy of the second ingredient. To allow robust comparison of results across tests within and between test facilities, and over time, protocols for bioefficacy testing must include either characterisation of the resistant strain, standardisation of the mosquitoes used for bioassays, or a combination of the two. Through a series of virtual meetings, key stakeholders and practitioners explored different approaches to achieving these goals. Via an iterative process we decided on the preferred approach and produced a protocol consisting of characterising mosquitoes used for bioefficacy testing before and after a round of bioassays, for example at each time point in a durability monitoring study. We present the final protocol and justify our approach to establishing a standard methodology for durability monitoring of ITNs containing pyrethroid and a second AI.

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

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.

Conserved and Divergent Aspects of Plasticity and Sexual Dimorphism in Wing Size and Shape in Three Diptera

The ability of powered flight in insects facilitated their great evolutionary success allowing them to occupy various ecological niches. Beyond this primary task, wings are often involved in various premating behaviors, such as the generation of courtship songs and the initiation of mating in flight. These specific functions imply special adaptations of wing morphology, as well as sex-specific wing morphologies. Although wing morphology has been extensively studied in Drosophila melanogaster (Meigen, 1830), a comprehensive understanding of developmental plasticity and the impact of sex on wing size and shape plasticity is missing for other Diptera. Therefore, we raised flies of the three Diptera species Drosophila melanogaster, Ceratitis capitata (Wiedemann, 1824) and Musca domestica (Linnaeus, 1758) at different environmental conditions and applied geometric morphometrics to analyze wing shape. Our data showed extensive interspecific differences in wing shape, as well as a clear sexual wing shape dimorphism in all three species. We revealed an impact of different rearing temperatures on wing shape in all three species, which was mostly explained by plasticity in wing size in D. melanogaster. Rearing densities had significant effects on allometric wing shape in D. melanogaster, while no obvious effects were observed for the other two species. Additionally, we did not find evidence for sex-specific response to different rearing conditions in D. melanogaster and C. capitata, while a male-specific impact of different rearing conditions was observed on non-allometric wing shape in M. domestica. Overall, our data strongly suggests that many aspects of wing morphology underly species-specific adaptations and we discuss potential developmental and functional implications of our results.

Landmark Data to Distinguish and Identify Morphologically Close Tabanus spp. (Diptera: Tabanidae)

Tabanus spp., also known as horse flies (Diptera: Tabanidae), are important vectors of several animal pathogens. Adult females of Tabanus megalops and Tabanus striatus, which are members of the T. striatus complex, are morphologically similar and hence difficult to distinguish using morphological characteristics. In addition, molecular identification by DNA barcoding is also unable to distinguish these species. These two species can occur sympatrically with Tabanus rubidus, which is morphologically similar to T. megalops and T. striatus. Wing geometric morphometrics has been widely used in various insects to distinguish morphologically similar species. This study explored the effectiveness of landmark-based geometrics at distinguishing and identifying T. megalops, T. rubidus, and T. striatus in Thailand. Specimens were collected from different geographical regions of Thailand, and only unambiguously identified specimens were used for geometric morphometric analyses. Left wings of females of T. megalops (n = 160), T. rubidus (n = 165), and T. striatus (n = 85) were photographed, and 22 wing landmarks were used for the analysis. Wing shape was able to distinguish among species with high accuracy scores, ranging from 94.38% to 99.39%. We showed that morphologically very close species of Tabanus can be reliably distinguished by the geometry of their wing venation, and we showed how our experimental material could be used as a reference to tentatively identify new field collected specimens.

On the classification of simple and complex biological images using Krawtchouk moments and Generalized pseudo-Zernike moments: a case study with fly wing images and breast cancer mammograms

In image analysis, orthogonal moments are useful mathematical transformations for creating new features from digital images. Moreover, orthogonal moment invariants produce image features that are resistant to translation, rotation, and scaling operations. Here, we show the result of a case study in biological image analysis to help researchers judge the potential efficacy of image features derived from orthogonal moments in a machine learning context. In taxonomic classification of forensically important flies from the Sarcophagidae and the Calliphoridae family (n = 74), we found the GUIDE random forests model was able to completely classify samples from 15 different species correctly based on Krawtchouk moment invariant features generated from fly wing images, with zero out-of-bag error probability. For the more challenging problem of classifying breast masses based solely on digital mammograms from the CBIS-DDSM database (n = 1,151), we found that image features generated from the Generalized pseudo-Zernike moments and the Krawtchouk moments only enabled the GUIDE kernel model to achieve modest classification performance. However, using the predicted probability of malignancy from GUIDE as a feature together with five expert features resulted in a reasonably good model that has mean sensitivity of 85%, mean specificity of 61%, and mean accuracy of 70%. We conclude that orthogonal moments have high potential as informative image features in taxonomic classification problems where the patterns of biological variations are not overly complex. For more complicated and heterogeneous patterns of biological variations such as those present in medical images, relying on orthogonal moments alone to reach strong classification performance is unrealistic, but integrating prediction result using them with carefully selected expert features may still produce reasonably good prediction models.

Geometric morphometrics to distinguish the cryptic species Anopheles minimus and An. harrisoni in malaria hot spot villages, western Thailand

Anopheles minimus Theobald 1901 and An. harrisoni Harbach & Manguin 2007 belong to the same species complex. They are morphologically similar and can exist in sympatry but have blood host preferences. The most accurate method for their identification is based on molecular techniques. Here, we measure the level of interspecific discrimination by geometric morphometry. Sixty-seven An. minimus and 22 An. harrisoni specimens were selected based on their morphological integrity and confirmed by identification polymerase chain reaction of internal transcribed spacer 2. These samples were used as reference data allowing for a morphometric identification based on geometric shape. Despite size overlap between the two species, there was a significant shape divergence allowing for differentiation of An. minimus and An. harrisoni with 90% accuracy. An intraspecific study of An. minimus showed a summer period associated to the reducing of wing size, which did not influence the shape-based differentiation of An. harrisoni. Wing venation geometry can be used to distinguish between these cryptic species mainly based on shaped divergence. This study suggests that geometric morphometrics represent a convenient low-cost method to complement morphological identification, especially concerning damaged specimens, i.e., insects having accidentally lost the anatomical features allowing a reliable morphological identification.

Geometric morphometrics and machine learning as tools for the identification of sibling mosquito species of the Maculipennis complex (Anopheles)

Geometric morphometrics allows researchers to use the specific software to quantify and to visualize morphological differences between taxa from insect wings. Our objective was to assess wing geometry to distinguish four Anopheles sibling species of the Maculipennis complex, An. maculipennis s. s., An. daciae sp. inq., An. atroparvus and An. melanoon, found in Northern Italy. We combined the geometric morphometric approach with different machine learning alghorithms: support vector machine (SVM), random forest (RF), artificial neural network (ANN) and an ensemble model (EN). Centroid size was smaller in An. atroparvus than in An. maculipennis s. s. and An. daciae sp. inq. Principal component analysis (PCA) explained only 33% of the total variance and appeared not very useful to discriminate among species, and in particular between An. maculipennis s. s. and An. daciae sp. inq. The performance of four different machine learning alghorithms using procrustes coordinates of wing shape as predictors was evaluated. All models showed ROC-AUC and PRC-AUC values that were higher than the random classifier but the SVM algorithm maximized the most metrics on the test set. The SVM algorithm with radial basis function allowed the correct classification of 83% of An. maculipennis s. s. and 79% of An. daciae sp. inq. ROC-AUC analysis showed that three landmarks, 11, 16 and 15, were the most important procrustes coordinates in mean wing shape comparison between An. maculipennis s. s. and An. daciae sp. inq. The pattern in the three-dimensional space of the most important procrustes coordinates showed a clearer differentiation between the two species than the PCA. Our study demonstrated that machine learning algorithms could be a useful tool combined with the wing geometric morphometric approach.

Evidence of Wing Shape Sexual Dimorphism in Aedes (Stegomyia) albopictus in Mallorca, Spain

The Asian tiger mosquito Aedes albopictus (Skuse, 1894) is a highly invasive species widely distributed on the Spanish Mediterranean coast and the Balearic archipelago. Most studies involving this species in Spain have been focused on surveillance and control methods. However, micro-evolutionary studies for Ae. albopictus in Spain have been traditionally neglected. Morphological diversity could be the result of long-term evolutionary diversification in responses to selective pressures such as temperature, precipitation, food availability, predation, or competition that may influence flight activity, host-seeking, and blood-feeding behavior. Wing geometric morphometric have been used not only to study micro- and macro-evolution in mosquitoes but also in studies of population structuring and sexual dimorphism. Therefore, the main goal of this study was to investigate the wing shape patterns of Ae. albopictus populations to unveil sexual dimorphism that could provide information about their ecology and behavior. Mosquito eggs were collected using oviposition traps at the main campus of the University of the Balearic Islands (Palma de Mallorca, Spain) and reared under laboratory conditions. In order to study wing shape variation patterns in Ae. albopictus males and females, the left wing of each adult mosquito was removed and analyzed based on 18 landmarks. Our results indicated strong levels of sexual dimorphism between Ae. albopictus males and females. Furthermore, according to the cross-validated reclassification test, males were correctly distinguished from females with an accuracy of 84% and females from males 75%. We observed a significant sexual dimorphism in the wing shape patterns of Ae. albopictus when considering different seasonal patterns (spring vs. autumn). Our results suggested that selective pressures may affect males differently to females. Host-seeking, blood-feeding, and oviposition behavior of females may act as a major driver for wing shape sexual dimorphism. These results should be considered for the development of more effective and targeted mosquito control strategies.

Curcumin/d-mannitol as photolarvicide: induced delay in larval development time, changes in sex ratio and reduced longevity of Aedes aegypti

BACKGROUND

Resistant populations of Ae. aegypti have been a major problem in arboviruses epidemic areas, generating a strong demand for novel methods of vector control. In this regard, our group has demonstrated the use of curcumin as an efficient photoactive larvicide to eliminate Ae. aegypti larvae. This work was aimed to evaluate the Ae. aegypti (Rockefeller) development under sublethal conditions, using a curcumin/d-mannitol (DMC) formulation. The photolarvicidal efficacy under semi-field and field conditions (wild populations) was also analyzed, as well as the photobleaching and residual activity of DMC.

RESULTS

A delay in development time when larvae were exposed to sublethal concentrations of DMC was observed, followed by significant changes in sex ratio and reduction in longevity. DMC also presented a low residual activity when compared to usual larvicides, and had a substantial photolarvicidal activity against wild populations in field trials, achieving 71.3% mortality after 48 h.

CONCLUSIONS

Overall, these findings are of great biological importance for the process of enabling the implementation of DMC as a new product in the control of Ae. aegypti larvae, and contributes to the improvement of new plant-based larvicides. © 2021 Society of Chemical Industry.

Geographical disjunction and environmental conditions drive intraspecific differentiation in the chalk-hill blue butterfly

Drivers of evolution are often related to geographical isolation and/or diverging environmental conditions. Spatial variation in neutral genetic markers mostly reflects past geographical isolation, i.e. long-lasting allopatry, whereas morphology is often driven by local environmental conditions, resulting in more rapid evolution. In Europe, most thermophilic species persisted during the past glacial periods in geographically disjunct refugia, representing long-lasting isolates, frequently with diverging environmental conditions. This situation has driven the evolution of intraspecific signatures in species. Here, we analysed wing shape and wing pigmentation of the chalk-hill blue butterfly, Polyommatus coridon, across its entire distribution range restricted to the western Palaearctic. In addition, we compiled abiotic environmental parameters for each sampling site. Wing colour patterns differentiated a western and an eastern lineage. These lineages might represent two main Pleistocene refugia and differentiation centres, one located on the Italian Peninsula and the other in the Balkan region. The two lineages showed evidence of hybridization across Central Europe, from the Alps and across Germany. The intraspecific differentiation was strongest in the width of the brown band on the outer margin of the wings. The morphological structures obtained are in line with genetic signatures found in previous studies, but the latter are more fine-grained. Current environmental conditions, such as mean temperatures, were only marginally correlated with colour patterns. Our study underlines that Pleistocene range shifts, often resulting in allopatric isolation, shape intraspecific phenotypic structures within species; that pigmentation responds in a more sensitive manner to spatial disjunction than wing shape; and that morphometric and genetic structures in P. coridon provide concordant patterns and thus support identical biogeographical conclusions.

WingBank: A Wing Image Database of Mosquitoes

Mosquito-borne diseases affect millions of people and cause thousands of deaths yearly. Vaccines have been hitherto insufficient to mitigate them, which makes mosquito control the most viable approach. But vector control depends on correct species identification and geographical assignment, and the taxonomic characters of mosquitoes are often inconspicuous to non-taxonomists, which are restricted to a life stage and/or even damaged. Thus, geometric morphometry, a low cost and precise technique that has proven to be efficient for identifying subtle morphological dissimilarities, may contribute to the resolution of these types of problems. We have been applying this technique for more than 10 years and have accumulated thousands of wing images with their metadata. Therefore, the aims of this work were to develop a prototype of a platform for the storage of biological data related to wing morphometry, by means of a relational database and a web system named “WingBank.” In order to build the WingBank prototype, a multidisciplinary team performed a gathering of requirements, modeled and designed the relational database, and implemented a web platform. WingBank was designed to enforce data completeness, to ease data query, to leverage meta-studies, and to support applications of automatic identification of mosquitoes. Currently, the database of the WingBank contains data referring to 77 species belonging to 15 genera of Culicidae. From the 13,287 wing records currently cataloged in the database, 2,138 were already made available for use by third parties. As far as we know, this is the largest database of Culicidae wings of the world.

Is There a Minimum Number of Landmarks That Optimizes the Geometric Morphometric Analysis of Mosquito (Diptera, Culicidae) Wings?

Culicids are the most significant arthropods affecting human health. Thus, their correct identification is critical. The use of Geometric Morphometrics (GM) has been recently incorporated into mosquito taxonomy and has begun to complement classic diagnostic techniques. Since sampling size depends on the number of Landmarks (LMs) used, this study aimed to establish the minimum number of wing LMs needed to optimize GM analysis of mosquito species and/or genera from urban and peri-urban areas of Argentina. Female left wings were used for the optimization phase, in which 17 LMs were reduced to four by iterative LM exclusion. To verify its efficiency, Principal Component Analysis (PCA), Discriminant Analysis (DA), and Canonical Variate Analysis (CVA) were performed. Additionally, a phenogram was constructed to visualize the results. We observed that five LMs for the PCA, CVA, and phenogram and nine for the DA enabled discrimination and/or clustering of almost all species and genera. Therefore, we tested the LM selection by using nine LMs and adding new species. The resulting PCA showed little overlap between species and almost all species clustered as expected, which was also reflected in the phenogram. Significant differences were found between wing shape among all species, together with a low total error rate in the DA. In conclusion, the number of LMs can be reduced and still be used to effectively differentiate and cluster culicids. This is helpful for better exploitation of available material and optimization of data processing time when classic taxonomy methods are inadequate or the material is scarce.

The influence of larval density on triacylglycerol content in Aedes aegypti (Linnaeus) (Diptera: Culicidae)

Given the importance of Aedes aegypti as an arbovirus vector, understanding the impact of population density during larval development is fundamental to improve control methods. The aims of this study were to analyze the effects of larval density on the triacylglycerol (TAG) reserves, fecundity, and the size of mosquitos. To evaluate the influence of density on the amount of TAG in the different stages of development of A. aegypti, the larvae were reared in densities of 0.04, 0.32, and 0.8 larvae/ml, represented by D1, D2, and D3, respectively. Larval rearing density affected the amount of TAG in the different developmental stages of this insect. The group reared with the highest population density (D3) presented greater storage of TAG, except when analyzing the amount of TAG μg/protein in the fat body of blood-fed females. In addition, it was possible to observe higher resistance to fasting and a greater number of eggs laid. Insects reared in D2 and D3 were larger than those in the D1 group. This study provides interesting data for future investigations on how the density of larval rearing affects the lipid metabolism for synthesis and mobilization of TAG stored in A. aegypti and how it influences size and reproduction.

Comprehensive characterisation of Culicoides clastrieri and C. festivipennis (Diptera: Ceratopogonidae) according to morphological and morphometric characters using a multivariate approach and DNA barcode

Biting midges are widespread around the world and play an essential role in the epidemiology of over 100 veterinary and medical diseases. For taxonomists, it is difficult to correctly identify species because of affinities among cryptic species and species complexes. In this study, species boundaries were examined for C. clastrieri and C. festivipennis and compared with six other Culicoides species. The classifiers are partial least squares discriminant analysis (PLS-DA) and sparse partial least squares discriminant analysis (sPLS-DA).The performance of the proposed method was evaluated using four models: (i) geometric morphometrics applied to wings; (ii) morphological wing characters, (iii) "Full wing" (landmarks and morphological characters) and (iv)  "Full model" (morphological characters-wing, head, abdomen, legs-and wing landmarks). Double cross-validation procedures were used to validate the predictive ability of PLS-DA and sPLS-DA models. The AUC (area under the ROC curve) and the balanced error rate showed that the sPLS-DA model performs better than the PLS-DA model. Our final sPLS-DA analysis on the full wing and full model, with nine and seven components respectively, managed to perfectly classify our specimens. The C. clastrieri and C. festivipennis sequences, containing both COI and 28S genes, revealed our markers' weak discrimination power, with an intraspecific and interspecific divergence of 0.4% and 0.1% respectively. Moreover, C. clastrieri and C. festivipennis are grouped in the same clade. The morphology and wing patterns of C. clastrieri and C. festivipennis can be used to clearly distinguish them. Our study confirms C. clastrieri and C. festivipennis as two distinct species. Our results show that caution should be applied when relying solely on DNA barcodes for species identification or discovery.

Interpopulation variations in life history traits and reproductive tactics in Aedes aegypti: A test on populations 50 km apart

The interpopulation variation in life history traits of a species reflects evolutionary adaptation in response to a local environment regime. We examined the life history traits of Aedes aegypti populations from 2 cities in southern Taiwan separated by 50 km. Results revealed a high level of trait differentiation in immature developmental time and survival of Ae. aegypti between the 2 cities. The Kaohsiung populations exhibited total pupation of 40%-60% on day 8; this was significantly lower than that of the Tainan populations and laboratory-reared KHsm mosquitos, which exhibited a pupation rate of 70%-90%. The slow immature development of the Kaohsiung populations was reflected in the low percentage of adult emergence (22%-26%) on day 10. The prolonged immature development did not select larger adults with longer life spans because the Kaohsiung populations had a shorter life span (≈37 d) than that of the Tainan populations (≈42 d). By contrast, immature development and longevity did not differ between populations within each region, indicating weak local differentiation. Three field populations exhibited male-bias sex ratio because of differential mortality of female immatures. The effect of female size on adult life history was nonsignificant. Two reproduction tactics were detected, representing the balanced-mortality hypothesis and the bet-hedging hypothesis. Despite their differential life history strategies and reproductive tactics, these mosquito populations apparently counterbalanced any shortcomings in traits to produce similar population growth. Maintaining optimal population density is essential for Aedes mosquitos to increase the probability of encountering mates and reduce the Allee effect.

Wing geometry and genetic analyses reveal contrasting spatial structures between male and female Aedes aegypti (L.) (Diptera: Culicidae) populations in metropolitan Manila, Philippines

The population genetic structure of Aedes aegypti (Linnaeus, 1762) has been studied in order to understand its role as an efficient vector. Several studies utilized an integrative approach; to combine genetic and phenotypic data to determine its population structure but these studies have only focused on female populations. To address this particular gap, our study compared the population variability and structuring between its male and female populations using phenotypic and genetic data from a highly-urbanized and dengue-endemic region of the Philippines, Metropolitan Manila. Five mosquito populations comprised of female (n = 137) and male (n = 49) adult mosquitoes were used in this study. All mosquito individuals underwent geometric morphometric (26 landmarks), and genetic (11 microsatellite loci) analyses. Results revealed that F_ST estimates (genetic) were 0.055 and 0.009 while Q_ST estimates (phenotypic) were 0.318 and 0.309 in in male and female populations, respectively. Wing shape variation plots showed that male populations were distinctly separated from each other while female populations overlapped. Similarly, discriminant analysis of principal components using genetic data revealed that male populations were also distinctly separated from each other while female populations showed near-overlapping populations. Genetic and phenetic dendrograms showed the formation of two groups in male populations but no groups in female populations. Further analysis indicated a significant correlation (r = 0.68, p = 0.02) between the genetic and phenetic distances of male populations. Bayesian analysis using genetic data also detected multiple clusters in male (K = 3) and female (K = 2) populations, while no clusters were detected using the phenotypic data from both sexes. Our results revealed contrasting phenotypic and genetic patterns between male and female Ae. aegypti, indicating that male populations were more spatially structured than female populations.

Microgeographic Wing-Shape Variation in Aedes albopictus and Aedes scapularis (Diptera: Culicidae) Populations

Aedes albopictus and Aedes scapularis are vectors of several arboviruses, including the dengue, chikungunya, and Rocio virus infection. While Ae. albopictus is a highly invasive species native to Asia and has been dispersed by humans to most parts of the world, Ae. scapularis is native to Brazil and is widely distributed in the southeast of the country. Both species are highly anthropophilic and are often abundant in places with high human population densities. Because of the great epidemiological importance of these two mosquitoes and the paucity of knowledge on how they have adapted to different urban built environments, we investigated the microgeographic population structure of these vector species in the city of São Paulo, Brazil, using wing geometric morphometrics. Females of Ae. albopictus and Ae. scapularis were collected in seven urban parks in the city. The right wings of the specimens were removed and digitized, and eighteen landmarks based on vein intersections in the wing venation patterns were used to assess cross-sectional variation in wing shape and size. The analyses revealed distinct results for Ae. albopictus and Ae. scapularis populations. While the former had less wing shape variation, the latter had more heterogeneity, indicating a higher degree of intraspecific variation. Our results indicate that microgeographic selective pressures exerted by different urban built environments have a distinct effect on wing shape patterns in the populations of these two mosquito species studied here.