Population genetic structure of the malaria vector Anopheles minimus in Thailand based on mitochondrial DNA markers

Fine-scale adaptive divergence and population genetic structure of Aedes aegypti in Metropolitan Manila, Philippines

BACKGROUND

The adaptive divergence of Aedes aegypti populations to heterogeneous environments can be a driving force behind the recent expansion of their habitat distribution and outbreaks of dengue disease in urbanized areas. In this study, we investigated the population genomics of Ae. aegypti at a regional scale in Metropolitan Manila, Philippines.

METHODS

We used the Pool-Seq double digestion restriction-site association DNA sequencing (ddRAD-Seq) approach to generate a high number of single nucleotide polymorphisms (SNPs), with the aim to determine local adaptation and compare the population structure with 11 microsatellite markers. A total of 217 Ae. aegypti individuals from seven female and seven male populations collected from Metropolitan Manila were used in the assays.

RESULTS

We detected 65,473 SNPs across the populations, of which 76 were non-neutral SNPs. Of these non-neutral SNPs, the multivariate regression test associated 50 with eight landscape variables (e.g. open space, forest, etc.) and 29 with five climate variables (e.g. air temperature, humidity, etc.) (P-value range 0.005-0.045) in female and male populations separately. Male and female populations exhibited contrasting spatial divergence, with males exhibiting greater divergence than females, most likely reflecting the different dispersal abilities of male and female mosquitoes. In the comparative analysis of the same Ae. aegypti individuals, the pairwise FST values of 11 microsatellite markers were lower than those of the neutral SNPs, indicating that the neutral SNPs generated via pool ddRAD-Seq were more sensitive in terms of detecting genetic differences between populations at fine-spatial scales.

CONCLUSIONS

Overall, our study demonstrates the utility of pool ddRAD-Seq for examining genetic differences in Ae. aegypti populations in areas at fine-spatial scales that could inform vector control programs such as Wolbachia-infected mosquito mass-release programs. This in turn would provide information on mosquito population dispersal patterns and the potential barriers to mosquito movement within and around the release area. In addition, the potential of environmental adaptability observed in Ae. aegypti could help population control efforts.

Population genetic structure and wing geometric morphometrics of the filarial vector Armigeres subalbatus (Diptera: Culicidae) in Thailand

Armigeres subalbatus (Diptera: Culicidae) is a mosquito species of significant medical and veterinary importance. It is widely distributed across Southeast and East Asia and is commonly found throughout Thailand. This study assessed the genetic diversity and population structure of Ar. subalbatus in Thailand using the cytochrome c oxidase subunit I (COI) gene sequences. Additionally, wing shape variations among these populations were examined using geometric morphometrics (GM). Our results demonstrated that the overall haplotype diversity (Hd) was 0.634, and the nucleotide diversity (π) was 0.0019. Significant negative values in neutrality tests (p < 0.05) indicate that the Ar. subalbatus populations in Thailand are undergoing a phase of expansion following a bottleneck event. The mismatch distribution test suggests that the populations may have started expanding approximately 16,678 years ago. Pairwise genetic differentiation among the 12 populations based on Fst revealed significant differences in 32 pairs (p < 0.05), with the degree of differentiation ranging from 0.000 to 0.419. The GM analysis of wing shape also indicated significant differences in nearly all pairs (p < 0.05), except for between populations from Nakhon Pathom and Samut Songkhram, and between those from Chiang Mai and Mae Hong Son, suggesting no significant difference due to their similar environmental settings. These findings enhance our understanding of the population structure and phenotypic adaptations of mosquito vectors, providing vital insights for the formulation of more efficacious vector control strategies.

Molecular characterization of anopheline mosquitoes from the goat malaria-endemic areas of Thailand

Despite the fact that over a 100 anopheline mosquito species have been identified as human malaria vectors, little is known about ungulate malaria vectors. Consequently, we focused on investigating the bionomics and genetic characterizations of anopheline mosquitoes in goat malaria-endemic regions. We also attempted to screen for ungulate malaria potential vectors. A total of 1019 female anopheline mosquitoes were collected from six goat farms in four provinces of Thailand from 2020 to 2021. Mosquitoes were morphologically identified and subsequently confirmed using the mitochondrial DNA barcoding region-cytochrome oxidase c subunit I (MtDNA-COI), mitochondrial DNA-cytochrome c oxidase subunit II (MtDNA-COII), and ribosomal DNA internal transcribed spacer 2 (rDNA-ITS2) sequences. The current study reveals the genetic characteristics and distribution of nine mosquito species within the Anopheles and Cellia subgenera. Four dominant species, including Anopheles peditaeniatus, Anopheles subpictus, Anopheles vagus, and Anopheles aconitus exhibited significant intraspecific gene flow within their corresponding species. Although malaria parasites were not found in 126 mosquito pools, meaning more investigation is necessary, the current study adds to the existing DNA barcoding data collection and improves the current understanding of the genetic structure and distribution of anopheline mosquito species, which could be useful for effective control of mosquito-borne diseases.

Genetic Diversity and Population Structure of Anopheles funestus in Western Kenya Based on Mitochondrial DNA Marker COII

The mitochondrial marker, COII, was employed to assess the genetic structure and diversity of Anopheles funestus, a very important malaria vector in Africa that adapt and colonize different ecological niches in western Kenya. Mosquitoes were collected using mechanical aspirators in four areas (Bungoma, Port Victoria, Kombewa, and Migori) in western Kenya. Following morphological identification, PCR was used to confirm the species. The COII gene was amplified, sequenced, and analyzed to determine genetic diversity and population structure. A total of 126 (Port Victoria-38, Migori-38, Bungoma-22, and Kombewa-28) sequences of COII were used for population genetic analysis. Anopheles funestus had a high haplotype diversity (Hd = 0.97 to 0.98) but low nucleotide diversity (Π = 0.004 to 0.005). The neutrality test revealed negative Tajima's D and Fs values indicating an excess of low-frequency variation. This could be attributed to either population expansion or negative selection pressure across all the populations. No genetic or structural differentiation (Fst = -0.01) and a high level of gene flow (Gamma St, Nm = 17.99 to 35.22) were observed among the populations. Population expansion suggests the high adaptability of this species to various ecological requirements, hence sustaining its vectorial capacity and malaria transmission.

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.

Multidisciplinary Investigations of Sustained Malaria Transmission in the Greater Mekong Subregion

In the course of malaria elimination in the Greater Mekong Subregion (GMS), malaria epidemiology has experienced drastic spatiotemporal changes with residual transmission concentrated along international borders and the rising predominance of Plasmodium vivax. The emergence of Plasmodium falciparum parasites resistant to artemisinin and partner drugs renders artemisinin-based combination therapies less effective while the potential spread of multidrug-resistant parasites elicits concern. Vector behavioral changes and insecticide resistance have reduced the effectiveness of core vector control measures. In recognition of these problems, the Southeast Asian International Center of Excellence for Malaria Research (ICEMR) has been conducting multidisciplinary research to determine how human migration, antimalarial drug resistance, vector behavior, and insecticide resistance sustain malaria transmission at international borders. These efforts allow us to comprehensively understand the ecology of border malaria transmission and develop population genomics tools to identify and track parasite introduction. In addition to employing in vivo, in vitro, and molecular approaches to monitor the emergence and spread of drug-resistant parasites, we also use genomic and genetic methods to reveal novel mechanisms of antimalarial drug resistance of parasites. We also use omics and population genetics approaches to study insecticide resistance in malaria vectors and identify changes in mosquito community structure, vectorial potential, and seasonal dynamics. Collectively, the scientific findings from the ICEMR research activities offer a systematic view of the factors sustaining residual malaria transmission and identify potential solutions to these problems to accelerate malaria elimination in the GMS.

Malaria Research for Tailored Control and Elimination Strategies in the Greater Mekong Subregion

The malaria landscape in the Greater Mekong Subregion has experienced drastic changes with the ramp-up of the control efforts, revealing formidable challenges that slowed down the progress toward malaria elimination. Problems such as border malaria and cross-border malaria introduction, multidrug resistance in Plasmodium falciparum, the persistence of Plasmodium vivax, the asymptomatic parasite reservoirs, and insecticide resistance in primary vectors require integrated strategies tailored for individual nations in the region. In recognition of these challenges and the need for research, the Southeast Asian International Center of Excellence for Malaria Research has established a network of researchers and stakeholders and conducted basic and translational research to identify existing and emerging problems and develop new countermeasures. The installation of a comprehensive disease and vector surveillance system at sentinel sites in border areas with the implementation of passive/active case detection and cross-sectional surveys allowed timely detection and management of malaria cases, provided updated knowledge for effective vector control measures, and facilitated the efficacy studies of antimalarials. Incorporating sensitive molecular diagnosis to expose the significance of asymptomatic parasite reservoirs for sustaining transmission helped establish the necessary evidence to guide targeted control to eliminate residual transmission. In addition, this program has developed point-of-care diagnostics to monitor the quality of artemisinin combination therapies, delivering the needed information to the drug regulatory authorities to take measures against falsified and substandard antimalarials. To accelerate malaria elimination, this program has actively engaged with stakeholders of all levels, fostered vertical and horizontal collaborations, and enabled the effective dissemination of research findings.