DNA barcoding and wing morphometrics to distinguish three Aedes vectors in 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.

Ctenocephalides orientis and Ctenocephalides felis in Thailand: Head geometry by species, sex and geography

Fleas in the genus Ctenocephalides serve as biological vectors or intermediate hosts of microorganisms such as bacteria, rickettsia, protozoa and helminths. Ctenocephalides felis has a worldwide distribution, while C. orientis has long been considered as a subspecies of C. felis in Asia. To help the morphological recognition of these two species and further explore their differences, we used the geometric morphometric approach applied to the head. Both sexes were examined. Five anatomical landmarks of the head were used, and to capture the curvature of the front head, 10 semilandmarks were added. There was a consistent difference in species classification accuracy when considering landmarks only versus their combination with semilandmarks, suggesting the importance of the curve of the head as a taxonomic signal. Using or not the labels in the reclassification analyses, the head shape allowed by itself almost perfect recognition of the two species, in both sexes, even after adjustment for prior probabilities. The same approach disclosed a high level of sexual size and shape dimorphism in both species. The contribution of size variation to the discrimination by shape was much more important between sexes (from 27% to 45%) than between species (from 0.7% to 7.1%). Nevertheless, in our data, size never could represent a way to reliably recognise the sex of an individual, even less its species. Geographical variation in head shape could only be explored for the C. orientis sample. No significant correlation of morphometric variation with geography could be detected, which would be consistent with gene flow between Thai provinces. The geometric morphometric approach of the flea head, when it incorporates head curves, is a promising tool for rapid, economical, and accurate species and sex identification. It is, therefore, a useful tool for future epidemiological and demographic studies.

Speciation patterns of Aedes mosquitoes in the Scutellaris Group: a mitochondrial perspective

The Scutellaris Group of Aedes comprises 47 mosquito species, including Aedes albopictus. While Ae. albopictus is widely distributed, the other species are mostly found in the Asia-Pacific region. Evolutionary history researches of Aedes species within the Scutellaris Group have mainly focused on Ae. albopictus, a species that raises significant public health concerns, neglecting the other species. In this study, we aimed to assess genetic diversity and estimate speciation times of several species within the Scutellaris Group. Mosquitoes were therefore collected from various Asia-Pacific countries. Their mitochondrial cytochrome c oxidase subunit 1 (cox1) and subunit 3 (cox3) sequences were analyzed alongside those of other Scutellaris Group species available in the GenBank database. To estimate the divergence time, we analyzed 1849 cox1 gene sequences from 21 species, using three species (Aedes aegypti, Aedes notoscriptus and Aedes vigilax) as outgroups. We found that most of the speciation dates occurred during the Paleogene and the Neogene periods. A separation between the Scutellaris Subgroup and the Albopictus Subgroup occurred approximately 64-61 million years ago (MYA). We also identified a split between species found in Asia/Micronesia and those collected in Melanesia/Polynesia approximately 36-35 MYA. Our findings suggest that the speciation of Aedes species within the Scutellaris Group may be driven by diversity in mammalian hosts, climate and environmental changes, and geological dynamics rather than human migration.

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

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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.

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.

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

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

•

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

•

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

•

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

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

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

Molecular techniques for the taxonomy of Aedes Meigen, 1818 (Culicidae: Aedini): A review of studies from 2010 to 2021

The original description of Aedes Meigen in 1818, written in Latin, was very brief and included a single species, Aedes cinereus. In the last two decades the genus Aedes (Meigen, 1818) has undergone several revisions and reclassifications, with the current proposal being described by Wilkerson in 2015. However, the available keys for morphological identification are still not sufficient to differentiate cryptic species, damaged species, or those with confusing taxonomy. The current study aims to identify and describe the main taxonomic proposals and molecular methodologies available for the identification of the genus Aedes published between the years 2010 and 2021. The main molecular techniques used to identify the genus in the last 10 years, are: Multiplex PCR, DNA barcoding, nuclear and mitochondrial markers, environmental DNA, and bacterial microbiome analysis. This review highlights that there are catalogued data for only a few species of the genus Aedes, being restricted to medically important taxa such as Aedes albopictus and Aedes aegypti. The integrative taxonomy approach is a possibility to reconcile morphological and molecular data to improve species delimitation, contributing to future revisions of the genus.

Species Discrimination of Three Odontomachus (Formicidae: Ponerinae) Species in Thailand Using Outline Morphometrics

Simple Summary

Determination of species in the ant genus Odontomachus, which is a venomous group of ants, may require the use of highly trained entomologists. In Thailand, three species are very similar and difficult to distinguish: O. monticola, O. rixosus, and O. simillimus. In such a situation, a complementary technique not requiring highly specialized entomological knowledge is welcome. The geometric morphometric approach has proven to be this sort of tool, especially powerful for morphologically close or even cryptic species. In its most recent development, the geometric method uses the relative position of some anatomical landmarks. However, in worker ants these landmarks are few in number and can be difficult to assess without dissection. Here, therefore, we use the outline-based approach, an alternative geometric technique that has not yet been tested in ants. We show that the simple outline of the head contains a strong taxonomic signal, much stronger than the one obtained from the pronotum shape. The outline technique therefore represents a promising approach to aid in the determination of ant species.

Abstract

All members of the ant genus Odontomachus Latreille, 1804 are venomous ants. Four species in this genus have been identified from Thailand: Odontomachus latidens Mayr, 1867; O. monticola Emery, 1892; O. rixosus Smith, 1757; and O. simillimus Smith, 1758. The three latter species are available and have been used for an outline morphometric study. They display similar morphology, which makes their distinction very difficult except for highly qualified individuals. A total of 80 worker specimens were studied, exploring the contour shapes of their head and pronotum as possible taxonomic characters. The size of each body part was estimated determining the contour perimeter, the values for which were largely overlapping between O. rixosus and O. simillimus; most O. monticola specimens exhibited a significantly larger size. In contrast to the size, each contour shape of the head or pronotum established O. rixosus as the most distinct species. An exploratory data analysis disclosed the higher taxonomic signal of the head contour relative to the pronotum one. The scores obtained for validated reclassification were much better for the head (99%) than for the pronotum (82%). This study supports outline morphometrics of the head as a promising approach to contribute to the morphological identification of ant species, at least for monomorphic workers.

Wing morphometrics as a tool for the identification of forensic important Lucilia spp. (Diptera: Calliphoridae)

Blow flies of the genera Lucilia Robineau-Desvoidy (Diptera: Calliphoridae) are considered forensically important species across several regions of the world. Due to the similarity of adults, especially females, the usual methods based on morphology or even molecular techniques can experience some limitations; therefore, alternative or supportive tools are required. Recently, the landmark-based geometric morphometric analysis has been applied to discriminate many insects on genus and species level. Herein, we focus on wing morphometric analysis as a tool in classifying five species of Lucilia; three species from Thailand - L. cuprina (Wiedemann, 1830), L. porphyrina (Walker, 1856) and L. sinensis Aubertin, 1933; and two species from Switzerland - L. caesar (Linnaeus, 1758) and L. illustris (Meigen, 1826). Canonical variate analysis of 233 right wings showed four overlapping clusters of L. cuprina, L. sinensis, L. caesar, and L. illustris with one distinct cluster of L. porphyrina. Eighty-eight to 100 percent of correct classification was achieved, with an UPGMA dendrogram analysis revealing clear-cut branch and sub-branch of five species determined. Results from this study suggested that wing morphometric analysis is a useful tool for the identification of adult Lucilia spp.

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&nbsp;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.

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.

Aedes albopictus diversity and relationships in south-western Europe and Brazil by rDNA/mtDNA and phenotypic analyses: ITS-2, a useful marker for spread studies

BACKGROUND

Aedes albopictus is a very invasive mosquito, which has recently colonized tropical and temperate regions worldwide. Of concern is its role in the spread of emerging or re-emerging mosquito-borne diseases. Ae. albopictus from south-western Europe and Brazil were studied to infer genetic and phenetic diversity at intra-individual, intra-population and inter-population levels, and to analyse its spread.

METHODS

Genotyping was made by rDNA 5.8S-ITS-2 and mtDNA cox1 sequencing to assess haplotype and nucleotide diversity, genetic distances and phylogenetic networks. Male and female phenotyping included combined landmark-and outlined-based geometric morphometrics of wing size and shape.

RESULTS

Specimens from seven populations from Spain, France and Brazil provided 12 cox1 and 162 5.8S-ITS-2 haplotypes, with great genetic variability difference between both markers (0.9% vs 31.2%). Five cox1 haplotypes were shared with other countries, mainly Italy, USA and China, but none was shared between Europe and Brazil. The 5.8S-ITS-2 showed 2-7 intra-individual (mean 4.7) and 16-34 intra-/inter-population haplotypes (24.7), including haplotypes shared between Spain, France and Brazil. A 4.3% of ITS-2 haplotypes were shared, mainly with Italy, USA and Thailand, evidencing worldwide spread and introductions from areas where recent outbreaks of Ae. albopictus-transmitted pathogens occurred. Wing size showed sex differences. Wing shape distinguished between Brazilian and European specimens. Both genetic and morphometric markers showed differences between insular Spain and continental Spain, France and Brazil.

CONCLUSIONS

ITS-2 proves to be a useful marker to assess Ae. albopictus spread, providing pronouncedly more information than cox1, including intra-individual, intra-population and inter-population levels, furnishing a complete overview of the evolutionary exchanges followed by this mosquito. Wing morphometry proves to be a useful phenotyping marker, allowing to distinguish different populations at the level of both male and female specimens. Results indicate the need for periodic surveillance monitorings to verify that no Ae. albopictus with high virus transmission capacity is introduced into Europe.

The Unequal Taxonomic Signal of Mosquito Wing Cells

Simple Summary

Mosquitoes of the genus Aedes include important vectors of human disease viruses, including dengue, chikungunya and Zika. Surveillance programs used to detect and control these pests need accurate, fast and low-cost techniques to track the primary target and monitor possible re-infestations. Geometric morphometrics of mosquito wings is a convenient tool in mosquito species identification, but this method requires a complete wing in good condition for maximum accuracy. In this study, we investigate the amount of taxonomic signal provided by shape analysis of the internal cells of the wing. We show that (i) the internal cells of the wing provide differing amounts of taxonomic information, and (ii) the taxonomic signal of a given cell depends on the species under comparison. Since some of these cells are very informative, our study suggests that even damaged wings may provide key taxonomic information to differentiate among species found in mixed species surveillance collections.

Abstract

Accurate identification of mosquito species is critically important for monitoring and controlling the impact of human diseases they transmit. Here, we investigate four mosquito species: Aedes aegypti, Ae. albopictus, Ae. scutellaris and Verrallina dux that co-occur in tropical and subtropical regions, and whose morphological similarity challenges their accurate identification, a crucial requirement in entomological surveillance programs. Previous publications reveal a clear taxonomic signal embedded in wing cell landmark configuration, as well as in the external contour of the wings. We explored this signal for internal cells of the wings as well, to determine whether internal cells could uniformly provide the same taxonomic information. For each cell to be tentatively assigned to its respective species, i.e., to measure the amount of its taxonomic information, we used the shape of its contour, rather than its size. We show that (i) the taxonomic signal of wing shape is not uniformly spread among internal cells of the wing, and (ii) the amount of taxonomic information of a given cell depends on the species under comparison. This unequal taxonomic signal of internal cells is not related to size, nor to apparent shape complexity. The strong taxonomic signal of some cells ensures that even partly damaged wings can be used to improve species recognition.

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.

Current Arboviral Threats and Their Potential Vectors in Thailand

Arthropod-borne viral diseases (arboviruses) are a public-health concern in many regions of the world, including Thailand. This review describes the potential vectors and important human and/or veterinary arboviruses in Thailand. The medically important arboviruses affect humans, while veterinary arboviruses affect livestock and the economy. The main vectors described are mosquitoes, but other arthropods have been reported. Important mosquito-borne arboviruses are transmitted mainly by members of the genus Aedes (e.g., dengue, chikungunya, and Zika virus) and Culex (e.g., Japanese encephalitis, Tembusu and West Nile virus). While mosquitoes are important vectors, arboviruses are transmitted via other vectors, such as sand flies, ticks, cimicids (Family Cimicidae) and Culicoides. Veterinary arboviruses are reported in this review, e.g., duck Tembusu virus (DTMUV), Kaeng Khoi virus (KKV), and African horse sickness virus (AHSV). During arbovirus outbreaks, to target control interventions appropriately, it is critical to identify the vector(s) involved and their ecology. Knowledge of the prevalence of these viruses, and the potential for viral infections to co-circulate in mosquitoes, is also important for outbreak prediction.

Identification of Aedes (Diptera: Culicidae) Species and Arboviruses Circulating in Arauca, Eastern Colombia

The identification of vector species and their natural infection with arboviruses results in important data for the control of their transmission. However, for the eastern region of Colombia, this information is limited. Therefore, this study morphologically and molecularly identified species of the genus Aedes and the detection of arboviruses (Dengue, Chikungunya, Zika, and Mayaro) in female mosquitoes (individually) present in three municipalities (Saravena, Arauquita, and Tame) by amplifying the genetic material using RT-PCR (reverse transcriptase polymerase chain reaction) in the department of Arauca, eastern Colombia. Inconsistencies between morphological and molecular identification were detected in 13 individuals with Aedes albopictus initially determined as Aedes aegypti based on morphology (n = 13). Molecular identification showed the simultaneous presence of A. aegypti (n = 111) and A. albopictus (n = 58) in the urban municipalities of Saravena and Arauquita. These individuals were naturally infected with Dengue virus type 1 (DENV-1) and Chikungunya virus (CHIKV). The most frequent arbovirus was DENV-1 with an infection rate of 40.7% (11/27) for A. aegypti and 39.7% (23/58) for A. albopictus, which was followed by CHIKV with an infection rate of 1.8% for A. aegypti (2/111) and 6.9% for A. albopictus (4/58). Additionally, a mixed infection of DENV-1 and CHIKV was obtained in 4.5% of A. aegypti (5/111). Zika virus (ZIKV) and Mayaro virus (MAYV) infections were not detected. This study found that barcoding (fragment gene COI) is a successful method for identifying Aedes species. Additionally, we recommend the individual processing of insects as a more accurate strategy for arboviruses detection since the infection rate is obtained and co-infection between DENV-1 and CHIKV is also possible.

Seasonal dynamics and molecular differentiation of three natural Anopheles species (Diptera: Culicidae) of the Maculatus group (Neocellia series) in malaria hotspot villages of Thailand

BACKGROUND

Anopheles sawadwongporni Rattanarithikul & Green, Anopheles maculatus Theobald and Anopheles pseudowillmori (Theobald) of the Anopheles maculatus group (Diptera: Culicidae) are recognized as potential malaria vectors in many countries from the Indian subcontinent through Southeast Asia to Taiwan. A number of malaria vectors in malaria hotspot areas along the Thai-Myanmar border belong to this complex. However, the species distribution and dynamic trends remain understudied in this malaria endemic region.

METHODS

Mosquitoes of the Maculatus group were collected using CDC light traps every other week from four villages in Tha Song Yang District, Tak Province, Thailand from January to December 2015. Adult female mosquitoes were morphologically identified on site using taxonomic keys. Molecular species identification was performed by multiplex PCR based on the internal transcribed spacer 2 (ITS2) region of ribosomal DNA (rDNA) and sequencing of the cox1 gene at a DNA barcoding region in a subset of 29 specimens.

RESULTS

A total of 1328 An. maculatus (sensu lato) female mosquitoes were captured with An. maculatus, An. sawadwongporni and An. pseudowilmori accounting for 75.2, 22.1 and 2.7% respectively. The field captured mosquitoes of the Maculatus group were most abundant in the wet season and had a preferred distribution in villages at higher elevations. The phylogenetic relationships of 29 cox1 sequences showed a clear-cut separation of the three member species of the Maculatus group, with the An. pseudowillmori cluster being separated from An. sawadwongporni and An. maculatus.

CONCLUSIONS

This study provides updated information for the species composition, seasonal dynamics and microgeographical distribution of the Maculatus group in malaria-endemic areas of western Thailand. This information can be used to guide the planning and implementation of mosquito control measures in the pursuance of malaria transmission.

New records and DNA barcoding of deer flies, Chrysops (Diptera: Tabanidae) in Thailand

Chrysops spp. or deer flies (Diptera: Tabanidae) are hematophagous flies of medical and veterinary importance and some species are important vectors of Trypanosoma evansi, the causative agent of surra in Thailand. However, data regarding deer fly species and their molecular identification are limited. Accurate species identification will indicate the appropriate control measures. In this study, an entomological survey of deer flies from different sites in Thailand between May 2018 and June 2019 were conducted. In addition, mitochondrial cytochrome oxidase subunit I (COI) barcoding region was used for species identification. A total of 82 females were collected and 6 species were identified. Of these, three species are new records for Thailand: C. designatus, C. fuscomarginalis and C. vanderwulpi bringing the species total found in Thailand to nine. The COI sequences revealed an intraspecific divergence of 0.0%-2.65% and an interspecific divergence of 7.03%-13.47%. Phylogenetic analysis showed that all deer fly species were clearly separated into distinct clusters according to morphologically identified species. These results indicated that COI barcodes were capable in discriminating between deer fly species on the basis of the barcoding gap and phylogenetic analysis. Therefore, DNA barcoding is a valuable tool for species identification of deer flies in Thailand.

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

Culex is the largest subgenus within the genus Culex that includes important vectors of diseases. The correct identification of mosquitoes is critical for effective control strategies. Wing geometric morphometrics (WGM) has been used to identify mosquito species alongside traditional identification methods. Here, WGM was used for eleven Culex species from São Paulo, Brazil, and one from Esquel, Argentina. Adult mosquitoes were collected using CDC (Centers for Disease Control) traps, morphologically identified and analyzed by WGM. The canonical variate analysis (CVA) was performed and a Neighbor-joining (NJ) tree was constructed to illustrate the patterns of species segregation. A cross-validated reclassification test was also carried out. From 110 comparisons in the cross-validated reclassification test, 87 yielded values higher than 70%, with 13 comparisons yielding 100% reclassification scores. Culexquinquefasciatus yielded the highest reclassification scores among the analyzed species, corroborating with the results obtained by the CVA, in which Cx. quinquefasciatus was the most distinct species. The high values obtained at the cross-validated reclassification test and in the NJ analysis as well as the segregation observed at the CVA made it possible to distinguish among Culex species with high degrees of confidence, suggesting that WGM is a reliable tool to identify Culex species of the subgenus Culex.

Molecular Identification and Geometric Morphometric Analysis of Haematobosca aberrans (Diptera: Muscidae)

The genus Haematobosca Bezzi, 1907 (Diptera: Muscidae) contains haematophagous flies of veterinary importance. A new fly species of this genus was recognised from northern Thailand based on morphological characters and described as Haematobosca aberrans Pont, Duvallet & Changbunjong, 2020. In the present study, the mitochondrial cytochrome c oxidase I (COI) gene was used to confirm the morphological identification of H. aberrans. In addition, landmark-based geometric morphometrics was used to determine sexual dimorphism. The molecular analysis was conducted with 10 COI sequences. The results showed that all sequences were 100% identical. The sequence was not highly similar to reference sequences from GenBank and did not match any identified species from Barcode of Life Data Systems (BOLD). Phylogenetic analysis clearly differentiated this species from other species within the subfamily Stomoxyinae. For geometric morphometric analysis, a total of 16 wing pictures were analysed using the landmark-based approach. The results showed significant differences in wing shape between males and females, with a cross-validated classification score of 100%. The allometric analysis showed that wing shape has no correlation with size. Therefore, the COI gene is effective in species identification of H. aberrans, and geometric morphometrics is also effective in determining sexual dimorphism.

Geometric morphometric analysis of the effect of temperature on wing size and shape in Aedes albopictus

Wing geometry helps to identify mosquito species, even cryptic ones. On the other hand, temperature has a well-known effect on insect metric properties. Can such effects blur the taxonomic signal embedded in the wing? Two strains of Aedes albopictus (laboratory and field strain) were examined under three different rearing temperatures (26, 30 and 33 °C) using landmark- and outline-based morphometric approaches. The wings of each experimental line were compared with Aedes aegypti. Both approaches indicated similar associations between wing size and temperature. For the laboratory strain, the wing size significantly decreased as the temperature increased. For the field strain, the largest wings were observed at the intermediate temperature. The two morphometric approaches describing shape showed different sensibilities to temperature. For both strains and sexes, the landmark-based approach disclosed significant wing shape changes with temperature changes. The outline-based approach showed lesser effects, detecting significant changes only in laboratory females and in field males. Despite the size and shape changes induced by temperature, the two strains of Ae. albopictus were always distinguished from Ae. aegypti. The present study confirms the lability of size. However, it also suggests that, despite environmentally-induced variation, the architecture of the wing still provides a strong taxonomic signal.

Effectiveness of landmark- and semi-landmark-based geometric morphometric to identify four species of Culex mosquitoes in Thailand

Objective:

This research aims to study the effectiveness of landmark- and semi-landmark-based geometric morphometric (LMSL) in the identification of four Culex vectors, namely, C. quinquefasciatus, C. visnui, C. sitiens, and C. whitmorei in Thailand and also compared the potential between LMSL and the landmark-based geometric morphometric method (LM).

Materials and Methods:

The overall size of the wing sample was estimated by the centroid size. Shape variables were computed as principal components of the “partial warp” calculated after generalized procrustes analysis of raw coordinates. Discriminant analysis of the canonical variables performed to explore the shape dissimilarity between Culex species has been shown as a factor map and to calculate the Mahalanobis distance. Size and shape differences based on pairwise Mahalanobis distances were tested using non-parametric methods (1,000 cycles) with Bonferroni correction at a p-value of <0.05.

Results:

A total of 120 individuals were used that were divided into 30 individuals per Culex species. The mean CS of C. sitiens had the largest wings followed by C. visnui, C. quinquefasciatus, and C. whitmorei in LM and LMSM. The patterns of statistical difference in CS of both methods were similar and wing shapes among Culex species were different based on a comparison of pairwise Mahalanobis distances (p < 0.05) in both methods. For the cross-validated reclassification test, LM provided Culex species separation ranging from 54% to 84% and 51% to 93% for LMSM.

Conclusion:

Thus, LMSM is another option to use for the identification in mosquito vectors that have a curved line on the wing specific to the species.

Geometric morphometrics approach towards discrimination of three member species of Maculatus group in Thailand

Members of the Maculatus group are important malaria vectors in the border regions of Thailand. However, the role of each species in malaria transmission remains unclear because of their highly similar morphologies, making them difficult to be differentiated. Whereas An. pseudowillmori may be identified by the color pattern of some scales on abdomen and wings, the distinction between An. maculatus and An. sawadwongporni relies on the wings only. Scales are labile structures, as they may be accidentally removed during capture and transportation to the laboratory. To discriminate among the species of this group, we tested the suitability of geometric techniques. Shape variables were used as input for discriminant analyses and validated reclassification. Both landmark- and outline-based geometric techniques disclosed significant differences between the three species. For the delicate An. maculatus - An. sawadwongporni distinction, the outline-based approach appeared as the most promising, with validated reclassification scores reaching 93%, as compared to 77% obtained by landmark data. For An. pseudowillmori, reclassification scores were 100% and 94%, respectively. Geometric morphometrics may provide an alternative and useful complement for discriminating members of the Maculatus group.

Geometric morphometrics of the scutum for differentiation of trombiculid mites within the genus Walchia (Acariformes: Prostigmata: Trombiculidae), a probable vector of scrub typhus

The vectors of scrub typhus are the larval stage of trombiculid mites, termed "chiggers". These vectors are very small - the larvae are approximately 0.2 mm in size - and therefore their morphological identification is difficult. Trombiculid mites are widely distributed across Asia and they can be identified at the genus level by the shape, size and setae/sensilla distribution of their dorsal chitin plate (scutum = shield), while morphological identification at the species level requires more mite characteristics. We recently developed a methodology to ascertain paired matched genotype and morphotype of individual chiggers, based on autofluorescence and brightfield microscopy with subsequent molecular identification using the COI gene (approximately 640bp length). However, based on 20 chigger specimens characterised by paired genotypic and morphological data consisting of the four species [Walchia ewingi with 2 subspecies]: Walchia ewingi lupella (n = 9), W. ewingi ewingi (n = 2), W. alpestris (n = 2), W. kritochaeta (n = 5) and W. minuscuta (n = 2) we found evidence of genetic polymorphism and morphological plasticity within the genus Walchia. The phylogenetic inference of the intra-genus relationships within the Walchia spp., based on COI gene (Blankaartia spp. served as outgroup), delineated the five included species by an average interspecific divergence of mean distance 0.218 (0.126 - 0.323). We therefore applied landmark-based and outline-based geometric morphometric (GM) approaches to differentiate Walchia species using scutum measurements. A total of 261 scutum images of Walchia spp. were examined by landmark-based GM (140 chigger specimens) and outline-based GM (121 specimens) techniques. All Walchia spp. showed significant differences in scutum size and shape. W. minuscuta showed the smallest mean scutum size in both techniques. The largest scutum was found in W. ewingi lupella and W. ewingi ewingi by landmark-based and outline-based GM analysis, respectively. The scutum shapes of W. alpestris and W. minuscuta were clearly distinguished from the other species. Cross-validated classification scores were different depending on species and digitizing techniques and landmark-based GM showed better scores than outline-based GM. We conclude that the morphologically closely-related trombiculid mite species can be further differentiated by their scutum features alone, using GM approaches. This technique is a promising tool for the much-needed characterization studies of chiggers and needs evaluation using matched morphometric and genotyping data for other genera of trombiculids.

Comparison of Landmark- and Outline-Based Geometric Morphometrics for Discriminating Mosquito Vectors in Ratchaburi Province, Thailand

It is often challenging to identify mosquito vectors in the field based on morphological features due to their similar morphologies and difficulties in obtaining undamaged samples but is required for their successful control. Geometric morphometrics (GM) overcomes this issue by analyzing a suite of traits simultaneously and has the added advantages of being easy to use, low cost, and quick. Therefore, this research compared the efficiency and precision of landmark- and outline-based GM techniques for separating species of mosquitoes in Huay Nam Nak village, Ratchaburi Province, Thailand. This research collected 273 individuals belonging to seven species: Anopheles barbirostris, An. subpictus, Culex quinquefasciatus, Cx. vishnui, Cx. whitmorei, Aedes aegypti, and Ae. albopictus. Both landmark-based and outline-based GM techniques could identify malaria vectors in this area to the genus level successfully and were also very effective for identifying the malaria vectors Anopheles spp. and the dengue vectors Aedes spp. to the species level. However, they were less effective for distinguishing between species of Culex. Therefore, GM represents a valuable tool for the identification of mosquito vectors in the field, which will facilitate their successful control.

A potential pitfall in studies of biological shape: Does size matter?

The number of published studies using geometric morphometrics (GM) for analysing biological shape has increased steadily since the beginning of the 1990s, covering multiple research areas such as ecology, evolution, development, taxonomy and palaeontology. Unfortunately, we have observed that many published studies using GM do not evaluate the potential allometric effects of size on shape, which normally require consideration or assessment. This might lead to misinterpretations and flawed conclusions in certain cases, especially when size effects explain a large part of the shape variation. We assessed, for the first time and in a systematic manner, how often published studies that have applied GM consider the potential effects of allometry on shape. We reviewed the 300 most recent published papers that used GM for studying biological shape. We also estimated how much of the shape variation was explained by allometric effects in the reviewed papers. More than one-third (38%) of the reviewed studies did not consider the allometric component of shape variation. In studies where the allometric component was taken into account, it was significant in 88% of the cases, explaining up to 87.3% of total shape variation. We believe that one reason that may cause the observed results is a misunderstanding of the process that superimposes landmark configurations, i.e. the Generalized Procrustes Analysis, which removes isometric effects of size on shape, but not allometric effects. Allometry can be a crucial component of shape variation. We urge authors to address, and report, size effects in studies of biological shape. However, we do not propose to always remove size effects, but rather to evaluate the research question with and without the allometric component of shape variation. This approach can certainly provide a thorough understanding of how much size contributes to the observed shaped variation.

Identification of Muscidae (Diptera) of medico-legal importance by means of wing measurements

Cadavers attract numerous species and genera of Muscidae, both regular elements of carrion insect assemblages, and accidental visitors. Identification of adult Muscidae may be considered difficult, particularly by non-experts. Since species identification is a vital first step in the analysis of entomological material in any forensic entomology orientated experiment and real cases, various alternative methods of species identification have been proposed. We investigated possibility of semiautomated identification by means of wing measurements as an alternative for classic morphology and DNA-based approaches. We examined genus-level identification success for 790 specimens representing 13 genera of the most common European cadavers visiting Muscidae. We found 99.8% of examined specimens correctly identified to the genus-level. Without error, the following were identified: Azelia, Eudasyphora, Graphomya, Hydrotaea, Musca, Muscina, Mydaea, Neomyia, Polietes, Stomoxys and Thricops. Genus-level misidentifications were found only in Helina and Phaonia. Discrimination of examined material on the species level within Hydrotaea (318 specimens representing eight species) and Muscina (163 specimens representing four species) showed lower, yet still high average identification success, 97.2 and 98.8%, respectively. Our results revealed relatively high success in both genus and species identification of Muscidae of medico-legal importance. Semiautomated identification by means of wing measurements can be used by non-experts and does not require sophisticated equipment. This method will facilitate the identification of forensically relevant muscids in comparison to more difficult and more time-consuming identification approaches based on taxonomic keys or DNA-based methods. However, for unambiguous identification of some taxa, we recommend complementary use of identification keys.

The Worldwide Spread of the Tiger Mosquito as Revealed by Mitogenome Haplogroup Diversity

In the last 40 years, the Asian tiger mosquito Aedes albopictus, indigenous to East Asia, has colonized every continent except Antarctica. Its spread is a major public health concern, given that this species is a competent vector for numerous arboviruses, including those causing dengue, chikungunya, West Nile, and the recently emerged Zika fever. To acquire more information on the ancestral source(s) of adventive populations and the overall diffusion process from its native range, we analyzed the mitogenome variation of 27 individuals from representative populations of Asia, the Americas, and Europe. Phylogenetic analyses revealed five haplogroups in Asia, but population surveys appear to indicate that only three of these (A1a1, A1a2, and A1b) were involved in the recent worldwide spread. We also found out that a derived lineage (A1a1a1) within A1a1, which is now common in Italy, most likely arose in North America from an ancestral Japanese source. These different genetic sources now coexist in many of the recently colonized areas, thus probably creating novel genomic combinations which might be one of the causes of the apparently growing ability of A. albopictus to expand its geographical range.