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

Phenotypic and genetic variation of Aedes albopictus (Diptera: Culicidae) in Thailand and its global relationships: Insights from wing morphometric and mitochondrial COI gene analyses

Aedes albopictus (Diptera: Culicidae), commonly known as the Asian tiger mosquito, is an important vector transmitting dangerous arboviruses to humans. This study investigated the phenotypic and genetic variation of this species in Thailand through wing geometric morphometric (GM) and mitochondrial cytochrome c oxidase subunit I (COI) gene sequence analyses. A total of 236 Ae. albopictus specimens from 12 populations in Thailand and 89 specimens from invasive populations in Florida, Hawaii and Brazil underwent wing GM analysis. The centroid size (CS) of Ae. albopictus populations in Thailand ranged from 2.00 mm in Bangkok to 2.36 mm in Chanthaburi, while in invasive populations, CS varied from 2.25 mm in Brazil to 2.47 mm in Florida. Pairwise comparisons of wing shape revealed significant differences for most population pairs, with distances ranging from 1.63 to 10.02. The clustering tree indicated distant relationships in wing shape between native and invasive populations. Additionally, partial COI gene sequences were amplified from 108 specimens, revealing a mean haplotype diversity of 0.842 ± 0.025 and a mean nucleotide diversity of 0.002 ± 0.001. The results from neutral Tajima's D and Fu's Fs tests indicated negative and statistically significant values (-2.159 and -33.846, respectively), suggesting population expansion. Further examination of haplotype relationships between Thailand and other countries identified two distinct groups: a Southeast Asia group, with Thai haplotypes clustered exclusively within it, and a non-Southeast Asia group. These findings highlight the phenotypic and genetic variation of Ae. albopictus in Thailand, providing essential insights for disease control strategies and tracing the mosquito's origins across regions.

Reliability of wing morphometrics for species identification of human-biting black flies (Diptera: Simuliidae) in Thailand

BACKGROUND

Fast and reliable species identification of black flies is essential for research proposes and effective vector control. Besides traditional identification based on morphology, which is usually supplemented with molecular methods, geometric morphometrics (GM) has emerged as a promising tool for identification. Despite its potential, no specific GM techniques have been established for the identification of black fly species.

METHODS

Adult female black flies collected using human bait, as well as those reared from pupae, were used in this study. Here, landmark-based GM analysis of wings was assessed for the first time to identify human-biting black fly species in Thailand, comparing this approach with the standard morphological identification method and DNA barcoding based on the mitochondrial cytochrome c oxidase subunit I (COI) gene. To explore genetic relationships between species, maximum likelihood (ML) and neighbor-joining (NJ) phylogenetic trees were built. Additionally, three different methods of species delimitation, i.e., assemble species by automatic partitioning (ASAP), generalized mixed yule coalescent (GMYC), and single Poisson tree processes (PTP), were utilized to identify the morphologically defined species. The effectiveness of a COI barcode in identifying black fly species was further examined through the best match (BM) and best close match (BCM) methods.

RESULTS

Seven black fly species, namely Simulium tenebrosum Takaoka, Srisuka & Saeung, 2018 (complex), S. doipuiense Takaoka & Choochote, 2005 (complex), S. nigrogilvum Summers, 1911, S. nodosum Puri, 1933, S. asakoae Takaoka & Davies, 1995, S. chamlongi Takaoka & Suzuki, 1984, and S. umphangense Takaoka, Srisuka & Saeung, 2017 were morphologically identified. Compared with the standard method, the GM analysis based on wing shape showed high success in separating species, achieving an overall accuracy rate of 88.54%. On the other hand, DNA barcoding surpassed wing GM for species identification with a correct identification rate of 98.57%. Species delimitation analyses confirmed the validity of most nominal species, with an exception for S. tenebrosum complex and S. doipuiense complex, being delimited as a single species. Moreover, the analyses unveiled hidden diversity within S. asakoae, indicating the possible existence of up to four putative species.

CONCLUSIONS

This study highlights the potential of wing GM as a promising and reliable complementary tool for species identification of human-biting black flies in Thailand.

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.

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.

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.