Variation over time in wing size and shape of the coastal malaria vector Anopheles (Cellia) epiroticus Linton and Harbach (Diptera: Culicidae) in Samut Songkhram, Thailand

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.

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.

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

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

Geometric Morphometric 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.

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.

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

Mosquitoes are hematophagous insects that transmit parasites and pathogens with devastating effects on humans, particularly in subtropical regions. Different mosquito species display various behaviors, breeding sites, and geographic distribution; however, they can be difficult to distinguish in the field due to morphological similarities between species and damage caused during trapping and transportation. Vector control methods for controlling mosquito-borne disease epidemics require an understanding of which vector species are present in the area as well as the epidemiological patterns of disease transmission. Although molecular techniques can accurately distinguish between mosquito species, they are costly and laborious, making them unsuitable for extensive use in the field. Thus, alternative techniques are required. Geometric morphometrics (GM) is a rapid and inexpensive technique that can be used to analyze the size, shape, and shape variation of individuals based on a range of traits. Here, we used GM to analyze the wings of 1,040 female mosquitoes from 12 different species in Thailand. The right wing of each specimen was removed, imaged microscopically, and digitized using 17 landmarks. Wing shape variation among genera and species was analyzed using canonical variate analysis (CVA), while discriminant function analysis was used to cross-validate classification reliability based on Mahalanobis distances. Phenetic relationships were constructed to illustrate the discrimination patterns for genera and species. CVA of the morphological variation among Aedes, Anopheles, Armigeres, Culex, and Mansonia mosquito genera revealed five clusters. In particular, we demonstrated a high percentage of correctly-distinguished samples among Aedes (97.48%), Armigeres (96.15%), Culex (90.07%), and Mansonia (91.67%), but not Anopheles (64.54%). Together, these findings suggest that wing landmark-based GM analysis is an efficient method for identifying mosquito species, particularly among the Aedes, Armigeres, Culex, and Mansonia genera.

Forced Egg Laying Method to Establish F1 Progeny from Field Populations and Laboratory Strains of Anopheles Mosquitoes (Diptera: Culicidae) in Thailand

Successful monitoring of physiological resistance of malaria vectors requires about 150 female mosquitoes for a single set of tests. In some situations, the sampling effort is insufficient due to the low number of field-caught mosquitoes. To address this challenge, we demonstrate the feasibility of using the forced oviposition method for producing F1 from field-caught Anopheles mosquitoes. A total of 430 and 598 gravid Anopheles females from four laboratory strains and five field populations, respectively, were tested. After blood feeding, gravid mosquitoes were individually introduced into transparent plastic vials, containing moistened cotton balls topped with a 4 cm2 piece of filter paper. The number of eggs, hatching larvae, pupation, and adult emergence were recorded daily. The mean number of eggs per female mosquito ranged from 39.3 for Anopheles cracens to 93.6 for Anopheles dirus in the laboratory strains, and from 36.3 for Anopheles harrisoni to 147.6 for Anopheles barbirostris s.l. in the field populations. A relatively high egg hatching rate was found in An. dirus (95.85%), Anopheles minimus (78.22%), and An. cracens (75.59%). Similarly, a relatively high pupation rate was found for almost all test species ranging from 66% for An. minimus to 98.7% for Anopheles maculatus, and lowest for An. harrisoni (43.9%). Highly successful adult emergence rate was observed among 85-100% of pupae that emerged in all tested mosquito populations. The in-tube forced oviposition method is a promising method for the production of sufficient F1 progeny for molecular identification, vector competence, insecticide resistance, and bioassay studies.

Spatial distribution of Plasmodium knowlesi cases and their vectors in Johor, Malaysia: in light of human malaria elimination

Background

Plasmodium knowlesi, a simian malaria parasite infection, increases as Plasmodium falciparum and Plasmodium vivax infections decrease in Johor, Malaysia. Therefore, this study aimed to identify the distribution of vectors involved in knowlesi malaria transmission in Johor. This finding is vital in estimating hotspot areas for targeted control strategies.

Methods

Anopheles mosquitoes were collected from the location where P. knowlesi cases were reported. Cases of knowlesi malaria from 2011 to 2019 in Johor were analyzed. Internal transcribed spacers 2 (ITS2) and cytochrome c oxidase subunit I (COI) genes were used to identify the Leucosphyrus Group of Anopheles mosquitoes. In addition, spatial analysis was carried out on the knowlesi cases and vectors in Johor.

Results

One hundred and eighty-nine cases of P. knowlesi were reported in Johor over 10 years. Young adults between the ages of 20–39 years comprised 65% of the cases. Most infected individuals were involved in agriculture and army-related occupations (22% and 32%, respectively). Four hundred and eighteen Leucosphyrus Group Anopheles mosquitoes were captured during the study. Anopheles introlatus was the predominant species, followed by Anopheles latens. Spatial analysis by Kriging interpolation found that hotspot regions of P. knowlesi overlapped or were close to the areas where An. introlatus and An. latens were found. A significantly high number of vectors and P. knowlesi cases were found near the road within 0–5 km.

Conclusions

This study describes the distribution of P. knowlesi cases and Anopheles species in malaria-endemic transmission areas in Johor. Geospatial analysis is a valuable tool for studying the relationship between vectors and P. knowlesi cases. This study further supports that the Leucosphyrus Group of mosquitoes might be involved in transmitting knowlesi malaria cases in Johor. These findings may provide initial evidence to prioritize diseases and vector surveillance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-021-03963-0.

Infrared light sensors permit rapid recording of wingbeat frequency and bioacoustic species identification of mosquitoes

Recognition and classification of mosquitoes is a critical component of vector-borne disease management. Vector surveillance, based on wingbeat frequency and other parameters, is becoming increasingly important in the development of automated identification systems, but inconsistent data quality and results frequently emerge from different techniques and data processing methods which have not been standardized on wingbeat collection of numerous species. We developed a simple method to detect and record mosquito wingbeat by multi-dimensional optical sensors and collected 21,825 wingbeat files from 29 North American mosquito species. In pairwise comparisons, wingbeat frequency of twenty six species overlapped with at least one other species. No significant differences were observed in wingbeat frequencies between and within individuals of Culex quinquefasciatus over time. This work demonstrates the potential utility of quantifying mosquito wingbeat frequency by infrared light sensors as a component of an automated mosquito identification system. Due to species overlap, wingbeat frequency will need to integrate with other parameters to accurately delineate species in support of efficient mosquito surveillance, an important component of disease intervention.

Preliminary Selection and Analysis of Deltamethrin-Resistant Strains of Aedes albopictus in the Laboratory

Background: Aedes albopictus is a major vector for transmission of many viral pathogens. Deltamethrin resistance was analyzed by catching A. albopictus in the field; the analysis effect is affected by many insecticides that often interact with A. albopictus in the field environment. Materials and Methods: This study examined the development of deltamethrin resistance in A. albopictus mosquitoes under controlled laboratory conditions, focusing on morphological changes, reproductive fitness, and mutation of the knockdown resistance (kdr) gene. Deltamethrin-resistant strains were selected up to the 20th generation. To determine the level of resistance, the lethal concentration 50 (LC₅₀) of deltamethrin in the larvae was obtained, followed by the resistance ratio (RR), and in adult mosquitoes, mortality rates were calculated using the contact tube method. Results: An increase in the LC₅₀ from 0.0070 to 0.0563 mg/L was observed in resistant versus sensitive strains, with an increase of 11.26 in the RR. Overall, the results of the larval resistance bioassay showed that resistant larvae had medium resistance; however, by the 20th generation, adult mosquitoes showed strong resistance. PCR amplification, cloning, and sequencing of sodium channel domain III gene fragments were subsequently carried out using selected resistant and sensitive female mosquitoes. As a result, a number of base mutations were observed in the kdr gene in the resistant strain; however, no amino acid sequence mutations were observed, suggesting that base sequence changes did not affect protein expression. Results of morphological changes between resistant and sensitive strains showed that significant differences in the body, foreleg, mid leg, and hind leg length, as well as wing length and width, antenna length, and proboscis length were observed between 18^th-generation resistant and sensitive strains of A. albopictus. On analysis of reproductive fitness associated with deltamethrin resistance in selection of mosquitoes, observation results showed differences between resistant and sensitive strains; the female/male ratio of mosquitoes decreased after pupa hatching, with more females and fewer males. Conclusions: The model of deltamethrin-resistant selection of A. albopictus was successfully established in the laboratory. The morphological phenotypes of the deltamethrin-resistant population of A. albopictus mosquitoes had changed. The kdr gene of the 19th and 20th generations of deltamethrin-resistant A. albopictus mosquitoes had silent mutations at several sites. After deltamethrin resistance selection, the female/male ratio of mosquitoes increased after pupa hatching, with more females and fewer males, hinting at increased chances of more female mosquitoes transmitting diseases.

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

In Thailand, Anopheles (Cellia) epiroticus Linton et Harbach (Diptera: Culicidae) is the secondary vector of human malaria along coastal regions. While there are some studies of phenotypic variability and population structure of A. epiroticus, more information on morphological variation would enhance epidemiological understanding of medically important mosquito vectors. This research examined morphological variation at three different distances from coastlines of Samut Songkhram Province, Thailand, using landmark-based geometric morphometrics. Wing shape of A. epiroticus was significantly different in the area 0.2 km away from the sea compared to areas 2 and 4 km away from the sea (p < 0.05). Phenotypic variability in wing shape is associated with distance from the sea. Morphological variations in the area closest to the sea were most pronounced, showing a relationship between A. epiroticus and the ecosystem that affects wing geometry. These results provide important information to understand morphological variation of A. epiroticus in coastal areas.