Complex population history of two Anopheles dirus mosquito species in Southeast Asia suggests the influence of Pleistocene climate change rather than human-mediated effects

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.

Island mosquitoes of Thailand: an update on species diversity and DNA barcoding

Mosquitoes (Diptera: Culicidae) are among the most medically significant insects, with several species acting as vectors for human pathogens. Although there are frequent reports of mosquito-borne diseases in the border island areas of Thailand, comprehensive data on the diversity and DNA barcoding of these mosquito species remain limited. This study investigated mosquito diversity in two main archipelagos in Thailand-the Trat archipelago (comprising Chang Island and Kood Island) and the Ranong archipelago (comprising Chang Island and Phayam Island)-and generated DNA barcode data from the mosquitoes found there. The survey across these islands discovered a total of 41 species, highlighting the presence of several species known to be vectors for human diseases. Thirty-seven mosquito species from the island areas were documented to provide reference DNA barcode sequences for mosquitoes in Thailand's island regions. Two species, Aedes fumidus and Finlaya flavipennis, have been added as new COI sequence records in the database. DNA barcoding was highly effective in classifying almost all species by identifying barcoding gaps, except for Anopheles baimaii and Anopheles dirus, which could not be distinguished. Additionally, the study noted that geographical variations might influence certain mosquito species, such as Anopheles barbirostris A3 and Mansonia dives, causing them to be split into two distinct subgroups. The findings of this study are crucial, as they aid in classifying mosquito species using molecular techniques and expand our knowledge of disease vectors in these biodiverse regions.

Insights into the mitochondrial cytochrome oxidase I (mt-COI) gene and wing morphometrics of Anopheles baimaii (Diptera: Culicidae) in malaria-endemic islands of Thailand

Anopheles baimaii (Diptera: Culicidae) significantly contributes to the transmission of parasites causing malaria in Southeast Asia and South Asia. This study examined the morphological (wing shape) and molecular (mitochondrial gene) variations of An. baimaii in four of Thailand's border islands, and also investigated the presence of Plasmodium parasites in these mosquitoes. No Plasmodium infections were detected in the samples. Significant differences in wing shape were observed in most island populations (p < 0.05). A single-linkage tree, constructed using Mahalanobis distances, clustered the populations into two groups based on geographical locations. Genetic variation in An. baimaii was also analyzed through cytochrome c oxidase subunit I (COI) gene sequences. This analysis identified 22 segregating sites and a low nucleotide diversity of 0.004. Furthermore, 18 distinct haplotypes were identified, indicating a high haplotype diversity of 0.825. Neutrality tests for the overall population revealed a significantly negative Fu's Fs value (-5.029), indicating a population expansion. In contrast, Tajima's D yielded a negative value (-1.028) that did not reach statistical significance. The mismatch distribution analysis exhibited a bimodal pattern, and the raggedness index was 0.068, showing no significant discrepancy (p = 0.485) between observed and expected distributions. Pairwise genetic differentiation assessments demonstrated significant differences between all populations (p < 0.05). These findings provide valuable insights into the COI gene and wing morphometric variations in An. baimaii across Thailand's islands, offering critical information for understanding the adaptations of this malaria vector and guiding future comprehensive research.

Co-occurrence of dual lineages within Simulium (Gomphostilbia) atratum De Meijere in the Indonesian Archipelago along Wallace's Line

Mitochondrial cytochrome c oxidase subunit I (COI) sequences were utilized to infer the population genetic structure of Simulium (Gomphostilbia) atratum De Meijere, an endemic simulid species to Indonesia. Both median-joining haplotype network and maximum-likelihood tree revealed two genetic lineages (A and B) within the species, with an overlap distribution in Lombok, which is situated along Wallace's line. Genetic differentiation and gene flow with varying frequencies (FST = 0.02-0.967; Nm = 0.01-10.58) were observed between populations of S. (G.) atratum, of which population pairs of different lineages showed high genetic differentiation. Notably, the high genetic distance of up to 5.92 % observed within S. (G.) atratum in Lombok was attributed to the existence of two genetically distinct lineages. The co-occurrence of distinct lineages in Lombok indicated that Wallace's line did not act as faunistic border for S. (G.) atratum in the present study. Moreover, both lineages also exhibited unimodal distributions and negative values of neutrality tests, suggesting a pattern of population expansion. The expansion and divergence time estimation suggested that the two lineages of S. (G.) atratum diverged and expanded during the Pleistocene era in Indonesia.

Population Genetic Structure and Population History of the Biting Midge Culicoides mahasarakhamense (Diptera: Ceratopogonidae)

Biting midges of the genus Culicoides Latreille are significant pests and vectors of disease agents transmitted to humans and other animals. Understanding the genetic structure and diversity of these insects is crucial for effective control programs. This study examined the genetic diversity, genetic structure, and demographic history of Culicoides mahasarakhamense, a possible vector of avian haemosporidian parasites and Leishmania martiniquensis, in Thailand. The star-like shape of the median joining haplotype network, a unimodal mismatch distribution, and significant negative values for Tajima's D and Fu's F_S tests indicated that populations had undergone recent expansion. Population expansion time was estimated to be 2000-22,000 years ago. Population expansion may have been triggered by climatic amelioration from cold/dry to warm/humid conditions at the end of the last glaciations, resulting in the increased availability of host blood sources. Population pairwise F_ST revealed that most (87%) comparisons were not genetically different, most likely due to a shared recent history. The exception to the generally low level of genetic structuring is a population from the northern region that is genetically highly different from others. Population isolation in the past and the limitation of ongoing gene flows due to large geographic distance separation are possible explanations for genetic differentiation.

Morphological and molecular identification reveals a high diversity of Anopheles species in the forest region of the Cambodia-Laos border

Background

To develop an effective malaria vector intervention method in forested international border regions within the Greater Mekong Subregion (GMS), more in-depth studies should be conducted on local Anopheles species composition and bionomic features. There is a paucity of comprehensive surveys of biodiversity integrating morphological and molecular species identification conducted within the border of Laos and Cambodia.

Methods

A total of 2394 adult mosquitoes were trapped in the Cambodia–Laos border region. We first performed morphological identification of Anopheles mosquitoes and subsequently performed molecular identification using 412 recombinant DNA–internal transcribed spacer 2 (rDNA-ITS2) and 391 mitochondrial DNA–cytochrome c oxidase subunit 2 (mtDNA-COII) sequences. The molecular and morphological identification results were compared, and phylogenetic analysis of rDNA-ITS2 and mtDNA-COII was conducted for the sequence divergence among species.

Results

Thirteen distinct species of Anopheles were molecularly identified in a 26,415 km² border region in Siem Pang (Cambodia) and Pathoomphone (Laos). According to the comparisons of morphological and molecular identity, the interpretation of local species composition for dominant species in the Cambodia–Laos border (An. dirus, An. maculatus, An. philippinensis, An. kochi and An. sinensis) achieved the highest accuracy of morphological identification, from 98.37 to 100%. In contrast, the other species which were molecularly identified were less frequently identified correctly (0–58.3%) by morphological methods. The average rDNA-ITS2 and mtDNA-COII interspecific divergence was respectively 318 times and 15 times higher than their average intraspecific divergence. The barcoding gap ranged from 0.042 to 0.193 for rDNA-ITS2, and from 0.033 to 0.047 for mtDNA-COII.

Conclusions

The Cambodia–Laos border hosts a high diversity of Anopheles species. The morphological identification of Anopheles species provides higher accuracy for dominant species than for other species. Molecular methods combined with morphological analysis to determine species composition, population dynamics and bionomic characteristics can facilitate a better understanding of the factors driving malaria transmission and the effects of interventions, and can aid in achieving the goal of eliminating malaria.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05167-0.

Population structure and population history of the black fly Simulium chumpornense (Diptera: Simuliidae) from Thailand

Understanding genetic structure and diversity of insect vectors is crucial for disease epidemiology. In this study, mitochondrial cytochrome c oxidase I sequences were used to infer genetic diversity, genetic structure and population history of the black fly, Simulium chumpornense Takaoka and Kuvangkadilok, a suspected vector of blood protozoa of the genus Leucocytozoon and Trypanosoma. High intraspecific genetic divergence (max. 3.76%) was found among 142 specimens obtained from 19 locations across Thailand. A median joining network revealed two genetic lineages (A and B) that were geographically associated. Lineage A is representative of central and northeastern regions. Lineage B represents specimens from diverse locations in northern, western, and southern Thailand, including the type locality. Mismatch distribution and the neutrality tests provided signals of past population expansions in both lineages. The expansion time dating back to the end of last glaciations at 12,000 - 15,000 years ago is possibly related to increasing of precipitation at the end of last glacial period. Despite recent population history, population pairwise F_ST analysis revealed that almost all population comparisons were genetically significantly different. The high level of genetic structuring is possibly a result of historical isolation of the population that survived in different refugia sites during the dry conditions during glaciations.

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

BACKGROUND

The malaria vector Anopheles minimus has been influenced by external stresses affecting the survival rate and vectorial capacity of the population. Since An. minimus habitats have continuously undergone ecological changes, this study aimed to determine the population genetic structure and the potential gene flow among the An. minimus populations in Thailand.

METHODS

Anopheles minimus was collected from five malaria transmission areas in Thailand using Centers for Disease Control and Prevention (CDC) light traps. Seventy-nine females from those populations were used as representative samples. The partial mitochondrial cytochrome c oxidase subunit I (COI), cytochrome c oxidase subunit II (COII) and cytochrome b (Cytb) gene sequences were amplified and analyzed to identify species and determine the current population genetic structure. For the past population, we determined the population genetic structure from the 60 deposited COII sequences in GenBank of An. minimus collected from Thailand 20 years ago.

RESULTS

The current populations of An. minimus were genetically divided into two lineages, A and B. Lineage A has high haplotype diversity under gene flow similar to the population in the past. Neutrality tests suggested population expansion of An. minimus, with the detection of abundant rare mutations in all populations, which tend to arise from negative selection.

CONCLUSIONS

This study revealed that the population genetic structure of An. minimus lineage A was similar between the past and present populations, indicating high adaptability of the species. There was substantial gene flow between the eastern and western An. minimus populations without detection of significant gene flow barriers.

Diversity and distribution of cytochrome oxidase I (COI) haplotypes of the brown marmorated stink bug, Halyomorpha halys Stål (Hemiptera, Pentatomidae), along the eastern front of its invasive range in Eurasia

The arrival, establishment and pest status of Halyomorpha halys in Europe and non-native countries in Asia have been well-documented, with thorough characterisation of the genetic diversity and occurrence of cytochrome oxidase I (COI) haplotypes in Switzerland, France, Hungary, Italy and Greece. However, a number of gaps exist in terms of the characterisation of the haplotype diversity and occurrence of H. halys along the invasion front that covers eastern Europe, western and central Asia. To contribute towards filling this gap, the COI haplotype diversity and distribution were investigated for H. halys collected in Serbia, Ukraine, Russia, Georgia and Kazakhstan. A total of 646 specimens were analysed and five haplotypes were found (H1, H3, H8, H33 and H80). Haplotype H1 was present in all five countries investigated and was the only haplotype detected amongst &gt; 500 specimens collected from Ukraine, Russia and Georgia. H1 (82%) was the dominant haplotype found in Kazakhstan, alongside H3 (18%). In contrast to the low or no diversity observed in these four countries, Serbia had higher haplotype diversity and was represented by five haplotypes. Although H3 was dominant (47%) in Serbia, H1 was also prevalent (40%); the remaining haplotypes (H8, H33 and H80) were minor contributors (1–11%) to the haplotype composition. The results are discussed in context with other known populations in neighbouring countries and patterns of haplotype diversity indicate the movement of successful invasive populations in Europe to generate secondary invasions along the eastern front of the invasion in Eurasia. Possible scenarios regarding the spread of particular haplotypes in these regions are discussed, along with suggestions for future research to fill existing gaps.

Highly variable COI haplotype diversity between three species of invasive pest fruit fly reflects remarkably incongruent demographic histories

Distance decay principles predict that species with larger geographic ranges would have greater intraspecific genetic diversity than more restricted species. However, invasive pest species may not follow this prediction, with confounding implications for tracking phenomena including original ranges, invasion pathways and source populations. We sequenced an 815 base-pair section of the COI gene for 441 specimens of Bactrocera correcta, 214 B. zonata and 372 Zeugodacus cucurbitae; three invasive pest fruit fly species with overlapping hostplants. For each species, we explored how many individuals would need to be included in a study to sample the majority of their haplotype diversity. We also tested for phylogeographic signal and used demographic estimators as a proxy for invasion potency. We find contrasting patterns of haplotype diversity amongst the species, where B. zonata has the highest diversity but most haplotypes were represented by singletons; B. correcta has ~7 dominant haplotypes more evenly distributed; Z. cucurbitae has a single dominant haplotype with closely related singletons in a 'star-shape' surrounding it. We discuss how these differing patterns relate to their invasion histories. None of the species showed meaningful phylogeographic patterns, possibly due to gene-flow between areas across their distributions, obscuring or eliminating substructure.

Genetic diversity and population structure of malaria vector mosquitoes Anopheles subpictus, Anopheles peditaeniatus, and Anopheles vagus in five districts of Sri Lanka

Background

Although Sri Lanka is considered as a malaria-free nation, the threat of re-emergence of outbreaks still remains due to the high prevalence and abundance of malaria vectors. Analysis of population genetic structure of malaria vectors is considered to be one of the vital components in implementing successful vector control programmes. The present study was conducted to determine the population genetic structure of three abundant malaria vectors; Anopheles subpictus sensu lato (s.l.), Anopheles peditaneatus and Anopheles vagus from five administrative districts in two climatic zones; intermediate zone (Badulla and Kurunegala districts) and dry zone (Ampara, Batticoloa and Jaffna districts) of Sri Lanka using the mitochondrial gene, cytochrome c oxidase subunit I (COI).

Methods

Adult mosquitoes of An. subpictus s.l., An. peditaeniatus, and An. vagus were collected from five study sites located in five districts using cattle baited traps and backpack aspirators. Representative samples of each species that were morphologically confirmed were selected from each locality in generating COI sequences (> 6 good quality sequences per species per locality).

Results

Anopheles subpictus s.l. specimens collected during the study belonged to two sibling species; An. subpictus ‘A’ (from all study sites except from Jaffna) and An. subpictus ‘B’ (only from Jaffna). The results of haplotype and nucleotide diversity indices showed that all the three species are having high genetic diversity. Although a high significant pairwise difference was observed between An. subpictus ‘A’ and ‘B’ (F_st> 0.950, p < 0.05), there were no significant genetic population structures within An. peditaeniatus, An. vagus and An. subpictus species A (p > 0.05), indicating possible gene flow between these populations.

Conclusions

Gene flow among the populations of An. peditaeniatus, An. vagus and An. subpictus species A was evident. Application of vector control measures against all mosquito species must be done with close monitoring since gene flow can assist the spread of insecticide resistance genes over a vast geographical area.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2419-x) contains supplementary material, which is available to authorized users.

Molecular Identification of Four Members of the Anopheles dirus Complex Using the Mitochondrial Cytochrome C Oxidase Subunit I Gene

Precise mosquito species identification is an essential step for proper management and control of malaria vectors. Misidentification of members in the Anopheles dirus complex, some which are primary malaria vectors in Thailand and mainland Southeast Asia, remains problematic because of indistinguishable or overlapping morphological characters between sibling species. Moreover, there is a need for alternative methods, since the existing molecular techniques in the literature are not entirely satisfactory in differentiating all members in the An. dirus complex. The nucleotide polymorphisms in the mitochondrial cytochrome c oxidase subunit I (COI) sequences were developed to identify the 4 species within the An. dirus complex using an allele-specific (AS) multiplex polymerase chain reaction (PCR). The identified primers amplified and clearly differentiated the 4 members of the complex found in Thailand, Anopheles dirus, An. cracens, An. scanloni, and An. baimaii with PCR products 428/104, 236, 625, and 428 bp, respectively. These results demonstrate that an AS-PCR based on the COI region can accurately identify 4 members of An. dirus complex and would be useful as an alternative PCR-based method for accurate species identification.

Genetic analysis of Bactrocera zonata (Diptera: Tephritidae) populations from India based on cox1 and nad1 gene sequences

The peach fruit fly, Bactrocera zonata, is among the most serious and polyphagous insect pest of fruit crops in many parts of the world under genus Bactrocera. In the present study, the genetic structure, diversity and demographic history of B. zonata in India were inferred from mitochondrial cytochrome oxidase 1 (cox1) and NADH dehydrogenase 1 (nad1) sequences. The efficiency of DNA barcodes for identification of B. zonata was also tested. Genetic diversity indices [number of haplotypes (H), haplotype diversity (Hd), nucleotide diversity (π) and average number of nucleotide differences (k)] of B. zonata populations across India maintain high level of genetic diversity without isolation by distance among the geographic regions. Non-significant negative correlation between pairwise Fst and geographic distance suggests a high level of gene flow among studied populations of B. zonata. The possibility of sudden expansion of B. zonata revealed through mismatch distribution analysis as well as negative Tajima's D and Fu's Fs values further supported by star-like network of haplotypes. DNA barcoding analysis suggests that B. zonata specimens can be clearly differentiated from other species with 100% accuracy of identification. Therefore, cytochrome oxidase 1 (cox1) barcode sequences generated in the present study could be a valuable source for the rapid identification and global population genetic study of B. zonata.

Phylogenetic Relationships Among Malaria Vectors and Closely Related Species in Thailand Using Multilocus DNA Sequences

The evolutionary and taxonomic status is important for understanding speciation events and phylogenetic relationships between closely related vector and nonvector species. This information is useful for targeting important disease vector species groups for the development of novel genetic-based vector and pathogen control methods. In this study, different phylogenetic analyses were performed to reconstruct phylogenetic trees for the primary malaria vectors in Thailand based on sequence information of 4 DNA fragments from the nuclear and mitochondrial regions. The primary Anopheles species in the subgenus Cellia involved in malaria transmission in Thailand separate clearly into 3 distinct clades: the Leucosphyrus group, Minimus subgroup, and Maculatus group. The phylogenetic trees based on different reconstructed algorithms and different gene regions provided congruent phylogenetic status of the mosquito species studied. The phylogenetic relationships of malaria vector species examined followed similar patterns based on morphological characters. An estimate of the divergence time among the Anopheles species infers that they were present during the Eocene and Miocene periods (>41 million years ago). Congruent phylogenetic analysis of malaria vectors is presented with different algorithms and gene regions. The nuclear TOLL6 fragment appears useful for molecular phylogenetic, species DNA barcode, and Anopheles population genetic analyses.

Population genetic structure and demographic history of the black fly vector, Simulium nodosum in Thailand

An understanding of the genetic structure and diversity of vector species is crucial for effective control and management. In this study, mitochondrial DNA sequences were used to examine the genetic structure, diversity and demographic history of a black fly vector, Simulium nodosum Puri (Diptera: Simuliidae), in Thailand. A total of 145 sequences were obtained from 10 sampling locations collected across geographical ranges in the country. Low genetic diversity was found in populations of S. nodosum that could be explained by the recent population history of this species. Demographic history analysis revealed a signature of demographic expansion dating back to only 2600-5200 years ago. Recent population expansion in S. nodosum possibly followed an increase in agriculture that enabled its hosts', humans and domestic animals, densities to increase. Alternatively, the Thai populations could be a derivative of an older expansion event in the more northern populations. Mitochondrial DNA genealogy revealed no genetically divergent lineages, which agrees with the previous cytogenetic study. Genetic structure analyses found that only 27% of the pairwise comparisons were significantly different. The most likely explanation for the pattern of genetic structuring is the effect of genetic drift because of recent colonization.

Cow-baited tents are highly effective in sampling diverse Anopheles malaria vectors in Cambodia

BACKGROUND

The accurate monitoring and evaluation of malaria vectors requires efficient sampling. The objective of this study was to compare methods for sampling outdoor-biting Anopheles mosquitoes in Cambodia.

METHODS

In the Cambodian provinces of Pursat, Preah Vihear, and Ratanakiri, six different mosquito trapping methods were evaluated: human landing collection (HLC), human-baited tent (HBT), cow-baited tent (CBT), CDC miniature light trap (LT), CDC miniature light trap baited with molasses and yeast (LT-M), and barrier fence (F) in a Latin square design during four or six consecutive nights at the height of the malaria transmission season.

RESULTS

Using all traps, a total of 507, 1175, and 615 anophelines were collected in Pursat, Preah Vihear, and Ratanakiri, respectively. CBTs captured 10- to 20-fold more anophelines per night than the other five sampling methods. All 2297 Anopheles mosquitoes were morphologically identified and molecularly typed using standard morphological keys and sequencing the rDNA ITS2 region to distinguish cryptic species, respectively. Overall, an extremely diverse set of 27 known Anopheles species was sampled. CBTs captured the same molecular species that HLCs and the other four traps did, as well as additional species. Nine specimens representing five Anopheles species (Anopheles hyrcanus, Anopheles barbirostris sensu stricto, Anopheles barbirostris clade III, Anopheles nivipes, and Anopheles peditaeniatus) were infected with Plasmodium falciparum and were exclusively captured in CBTs.

CONCLUSIONS

These data indicate that cow-baited tents are highly effective in sampling diverse Anopheles malaria vectors in Cambodia. This sampling method captured high numbers of anophelines with limited sampling effort and greatly reduced human exposure to mosquito bites compared to the gold-standard human landing collection.

Red River barrier and Pleistocene climatic fluctuations shaped the genetic structure of Microhyla fissipes complex (Anura: Microhylidae) in southern China and Indochina

South China and Indochina host striking species diversity and endemism. Complex tectonic and climatic evolutions appear to be the main drivers of the biogeographic patterns. In this study, based on the geologic history of this region, we test 2 hypotheses using the evolutionary history of Microhyla fissipes species complex. Using DNA sequence data from both mitochondrial and nuclear genes, we first test the hypothesis that the Red River is a barrier to gene flow and dispersal. Second, we test the hypothesis that Pleistocene climatic cycling affected the genetic structure and population history of these frogs. We detect 2 major genetic splits that associate with the Red River. Time estimation suggests that late Miocene tectonic movement associated with the Red River drove their diversification. Species distribution modeling (SDM) resolves significant ecological differences between sides of the Red River. Thus, ecological divergence also probably promoted and maintained the diversification. Genogeography, historical demography, and SDM associate patterns in southern China with climate changes of the last glacial maximum (LGM), but not Indochina. Differences in geography and climate between the 2 areas best explain the discovery. Responses to the Pleistocene glacial–interglacial cycling vary among species and regions.

Phylogeography of the black fly Simulium aureohirtum (Diptera: Simuliidae) in Thailand

The black fly Simulium aureohirtum Brunetti (Diptera: Simuliidae) is geographically widespread and utilizes diverse habitats. Previous studies revealed high genetic diversity, suggesting that this species could be a species complex. In this study, we used mitochondrial cytochrome c oxidase II sequences to examine genetic structure and demographic history of S. aureohirtum in Thailand. High level of genetic variation was observed due to the existence of genetically divergent lineages (A, B, and C). Lineages A and C were geographically widespread, while lineage B was geographically restricted to the eastern part of Thailand. The geographic localization of lineage B suggested limitations in gene flow, which is most likely a result of ecologically different habitats. Large genetic differences between individuals of lineages A and C, but geographically sympatric distributions of the members of these lineages, are considered to represent an admixture of the populations that have evolved allopatrically. Historical environmental change during the Pleistocene glaciations is possibly the factor that drove the lineage divergence. Population demographic history analyses revealed recent population expansions in the three lineages dating back to the last glaciations. These observations further highlight the significance of the Pleistocene climatic change on current genetic structure and diversity of living organisms of the Southeast Asian mainland.

DNA barcoding for identification of sand fly species (Diptera: Psychodidae) from leishmaniasis-endemic areas of Peru

Phlebotomine sand flies are the only proven vectors of leishmaniases, a group of human and animal diseases. Accurate knowledge of sand fly species identification is essential in understanding the epidemiology of leishmaniasis and vector control in endemic areas. Classical identification of sand fly species based on morphological characteristics often remains difficult and requires taxonomic expertise. Here, we generated DNA barcodes of the cytochrome c oxidase subunit 1 (COI) gene using 159 adult specimens morphologically identified to be 19 species of sand flies, belonging to 6 subgenera/species groups circulating in Peru, including the vector species. Neighbor-joining (NJ) analysis based on Kimura 2-Parameter genetic distances formed non-overlapping clusters for all species. The levels of intraspecific genetic divergence ranged from 0 to 5.96%, whereas interspecific genetic divergence among different species ranged from 8.39 to 19.08%. The generated COI barcodes could discriminate between all the sand fly taxa. Besides its success in separating known species, we found that DNA barcoding is useful in revealing population differentiation and cryptic diversity, and thus promises to be a valuable tool for epidemiological studies of leishmaniasis.

A review of malaria transmission dynamics in forest ecosystems

Malaria continues to be a major health problem in more than 100 endemic countries located primarily in tropical and sub-tropical regions around the world. Malaria transmission is a dynamic process and involves many interlinked factors, from uncontrollable natural environmental conditions to man-made disturbances to nature. Almost half of the population at risk of malaria lives in forest areas. Forests are hot beds of malaria transmission as they provide conditions such as vegetation cover, temperature, rainfall and humidity conditions that are conducive to distribution and survival of malaria vectors. Forests often lack infrastructure and harbor tribes with distinct genetic traits, socio-cultural beliefs and practices that greatly influence malaria transmission dynamics. Here we summarize the various topographical, entomological, parasitological, human ecological and socio-economic factors, which are crucial and shape malaria transmission in forested areas. An in-depth understanding and synthesis of the intricate relationship of these parameters in achieving better malaria control in various types of forest ecosystems is emphasized.

Evidence for the occurrence of two sympatric sibling species within the Anopheles (Kerteszia) cruzii complex in southeast Brazil and the detection of asymmetric introgression between them using a multilocus analysis

BACKGROUND

Anopheles (Kerteszia) cruzii (Diptera: Culicidae) is a primary vector of human and simian malaria parasites in southern and southeastern Brazil. Earlier studies using chromosome inversions, isoenzymes and a number of molecular markers have suggested that An. cruzii is a species complex.

RESULTS

In this study, a multilocus approach using six loci, three circadian clock genes and three encoding ribosomal proteins, was carried out to investigate in more detail the genetic differentiation between the An. cruzii populations from Florianópolis-Santa Catarina (southern Brazil) and Itatiaia-Rio de Janeiro States (southeastern Brazil). The analyses were performed first comparing Florianópolis and Itatiaia, and then comparing the two putative sympatric incipient species from Itatiaia (Itatiaia A and Itatiaia B). The analysis revealed high FST values between Florianópolis and Itatiaia (considering Itatiaia A and B together) and also between the sympatric Itatiaia A and Itatiaia B, irrespective of their function. Also, using the IM program, no strong indication of migration was found between Florianópolis and Itatiaia (considering Itatiaia A and B together) using all loci together, but between Itatiaia A and Itatiaia B, the results show evidence of migration only in the direction of Itatiaia B.

CONCLUSIONS

The results of the multilocus analysis indicate that Florianópolis and Itatiaia represent different species of the An. cruzii complex that diverged around 0.6 Mya, and also that the Itatiaia sample is composed of two sympatric incipient species A and B, which diverged around 0.2 Mya. Asymmetric introgression was found between the latter two species despite strong divergence in some loci.

An epidemiological overview of malaria in Bangladesh

Bangladesh is one of the four major malaria-endemic countries in South-East Asia having approximately 34% of its population at risk of malaria. This paper aims at providing an overview of the malaria situation in this country. Relevant information was retrieved from published articles and reports in PubMed and Google Scholar. Malaria in Bangladesh is concentrated in 13 districts with a prevalence ranging between 3.1% and 36%, and is mostly caused by Plasmodium falciparum. Geographical conditions pose a potential risk for Plasmodium knowlesi malaria. Resistance to a number of drugs previously recommended for treatment has been reported. Low socio-economic status, poor schooling and close proximity to water bodies and forest areas comprise important risk factors. Despite the significant steps in Long Lasting Insecticide Net (LLIN)/Insecticide Treated Net (ITN) coverage in Bangladesh, there are still many challenges including the extension of malaria support to the remote areas of Bangladesh, where malaria prevalence is higher, and further improvements in the field of referral system and treatment.

Predicted distribution of major malaria vectors belonging to the Anopheles dirus complex in Asia: ecological niche and environmental influences

Methods derived from ecological niche modeling allow to define species distribution based on presence-only data. This is particularly useful to develop models from literature records such as available for the Anopheles dirus complex, a major group of malaria mosquito vectors in Asia. This research defines an innovative modeling design based on presence-only model and hierarchical framework to define the distribution of the complex and attempt to delineate sibling species distribution and environmental preferences. At coarse resolution, the potential distribution was defined using slow changing abiotic factors such as topography and climate representative for the timescale covered by literature records of the species. The distribution area was then refined in a second step using a mask of current suitable land cover. Distribution area and ecological niche were compared between species and environmental factors tested for relevance. Alternatively, extreme values at occurrence points were used to delimit environmental envelopes. The spatial distribution for the complex was broadly consistent with its known distribution and influencing factors included temperature and rainfall. If maps developed from environmental envelopes gave similar results to modeling when the number of sites was high, the results were less similar for species with low number of recorded presences. Using presence-only models and hierarchical framework this study not only predicts the distribution of a major malaria vector, but also improved ecological modeling analysis design and proposed final products better adapted to malaria control decision makers. The resulting maps can help prioritizing areas which need further investigation and help simulate distribution under changing conditions such as climate change or reforestation. The hierarchical framework results in two products one abiotic based model describes the potential maximal distribution and remains valid for decades and the other including a biotic mask easy to update with frequently available information gives current species distribution.

Genetic population structure of the malaria vector Anopheles baimaii in north-east India using mitochondrial DNA

BACKGROUND

Anopheles baimaii is a primary vector of human malaria in the forest settings of Southeast Asia including the north-eastern region of India. Here, the genetic population structure and the basic population genetic parameters of An. baimaii in north-east India were estimated using DNA sequences of the mitochondrial cytochrome oxidase sub unit II (COII) gene.

METHODS

Anopheles baimaii were collected from 26 geo-referenced locations across the seven north-east Indian states and the COII gene was sequenced from 176 individuals across these sites. Fifty-seven COII sequences of An. baimaii from six locations in Bangladesh, Myanmar and Thailand from a previous study were added to this dataset. Altogether, 233 sequences were grouped into eight population groups, to facilitate analyses of genetic diversity, population structure and population history.

RESULTS

A star-shaped median joining haplotype network, unimodal mismatch distribution and significantly negative neutrality tests indicated population expansion in An. baimaii with the start of expansion estimated to be ~0.243 million years before present (MYBP) in north-east India. The populations of An. baimaii from north-east India had the highest haplotype and nucleotide diversity with all other populations having a subset of this diversity, likely as the result of range expansion from north-east India. The north-east Indian populations were genetically distinct from those in Bangladesh, Myanmar and Thailand, indicating that mountains, such as the Arakan mountain range between north-east India and Myanmar, are a significant barrier to gene flow. Within north-east India, there was no genetic differentiation among populations with the exception of the Central 2 population in the Barail hills area that was significantly differentiated from other populations.

CONCLUSIONS

The high genetic distinctiveness of the Central 2 population in the Barail hills area of the north-east India should be confirmed and its epidemiological significance further investigated. The lack of genetic population structure in the other north-east Indian populations likely reflects large population sizes of An. baimaii that, historically, were able to disperse through continuous forest habitats in the north-east India. Additional markers and analytical approaches are required to determine if recent deforestation is now preventing ongoing gene flow. Until such information is acquired, An. baimaii in north-east India should be treated as a single unit for the implementation of vector control measures.

Towards a new role for vector systematics in parasite control

Vector systematics research is being transformed by the recent development of theoretical, experimental and analytical methods, as well as conceptual insights into speciation and reconstruction of evolutionary history. We review this progress using examples from the mosquito genus Anopheles. The conclusion is that recent progress, particularly in the development of better tools for understanding evolutionary history, makes systematics much more informative for vector control purposes, and has increasing potential to inform and improve targeted vector control programmes.

DNA barcodes confirm the presence of a single member of the Anopheles maculipennis group in Morocco and Algeria: An. sicaulti is conspecific with An. labranchiae

Anopheles labranchiae Falleroni is the only member of the Maculipennis Group known to occur in northern Africa; however, confusion exists as to the taxonomic status of its junior synonym, An. sicaulti Roubaud (type locality: near Rabat, Morocco). Based on morphological and behavioural distinctions, it has been suggested that Moroccan populations have been isolated from other North African populations by the Atlas Mountains, and that Moroccan populations may represent An. sicaulti, originally described as a variety of An. maculipennis Meigen. DNA barcodes (658bp of the mitochondrial COI gene) obtained from 89 An. maculipennis s.l. collected in Morocco (n=79) and Algeria (n=10) in 2007 and 2008 were used to determine if Moroccan populations are genetically isolated from those east of the Atlas Mountains (Algeria), and whether there is molecular evidence to support the presence of more than one member of the Maculipennis Group in the region. No evidence for speciation was found between Moroccan and Algerian populations, or within populations in northern Morocco. Moreover shared COI haplotypes between Algeria and Morocco indicate ongoing gene flow between populations in these countries, suggesting that the Atlas Mountains are not a boundary to gene flow in An. labranchiae. The synonymy of An. sicaulti with An. labranchiae is confirmed. That An. labranchiae comprises the same species in these North African countries is important for malaria control.

Genetic and phenotypic variation of the malaria vector Anopheles atroparvus in southern Europe

BACKGROUND

There is a growing concern that global climate change will affect the potential for pathogen transmission by insect species that are vectors of human diseases. One of these species is the former European malaria vector, Anopheles atroparvus. Levels of population differentiation of An. atroparvus from southern Europe were characterized as a first attempt to elucidate patterns of population structure of this former malaria vector. Results are discussed in light of a hypothetical situation of re-establishment of malaria transmission.

METHODS

Genetic and phenotypic variation was analysed in nine mosquito samples collected from five European countries, using eight microsatellite loci and geometric morphometrics on 21 wing landmarks.

RESULTS

Levels of genetic diversity were comparable to those reported for tropical malaria vectors. Low levels of genetic (0.004 <FST <0.086) and phenotypic differentiation were detected among An. atroparvus populations spanning over 3,000 km distance. Genetic differentiation (0.202 <FST <0.299) was higher between the sibling species An. atroparvus and Anopheles maculipennis s.s. Differentiation between sibling species was not so evident at the phenotype level.

CONCLUSIONS

Levels of population differentiation within An. atroparvus were low and not correlated with geographic distance or with putative physical barriers to gene flow (Alps and Pyrenées). While these results may suggest considerable levels of gene flow, other explanations such as the effect of historical population perturbations can also be hypothesized.

Late Pleistocene environmental changes lead to unstable demography and population divergence of Anopheles albimanus in the northern Neotropics

We investigated the historical demography of Anopheles albimanus using mosquitoes from five countries and three different DNA regions, the mitochondrial cytochrome oxidase subunit I gene (COI), the single copy nuclear white gene and the ribosomal internal transcribed spacer two (ITS2). All the molecular markers supported the taxonomic status of a single species of An. albimanus. Furthermore, agreement between the COI and the white genes suggested a scenario of Pleistocene geographic fragmentation (i.e., population contraction) and subsequent range expansion across southern Central America.

Persistent transmission of malaria in Garo hills of Meghalaya bordering Bangladesh, north-east India

Background

Malaria is endemic in Garo hills of Meghalaya, and death cases are reported annually. Plasmodium falciparum is the major parasite, and is solely responsible for each malaria-attributable death case. Garo hills are categorized high-risk for drug-resistant malaria; however, there exists no data on malaria transmitting mosquitoes prevalent in the region. Included in this report are entomological observations with particular reference to vector biology characteristics for devising situation specific intervention strategies for disease transmission reduction.

Methods

The epidemiological data of the West Garo hills have been reviewed retrospectively for 2001-2009 to ascertain the disease transmission profile given the existing interventions. Point prevalence study was conducted in Dalu Community Health Centre that lies in close proximity to international border with Bangladesh to ascertain the true prevalence of malaria, and parasite species. Mosquito collections were made in human dwellings of malaria endemic villages aiming at vector incrimination, and to study relative abundance, resting and feeding preferences, and their present susceptibility status to DDT.

Results

Investigations revealed that the West Garo hill district is co-endemic for Plasmodium falciparum and Plasmodium vivax, but P. falciparum was the predominant infection (> 82%). Malaria transmission was perennial and persistent with seasonal peak during May-July corresponding to months of high rainfall. Entomological collections revealed that Anopheles minimus was the predominant species that was incriminated by detection of sporozoites in salivary glands (infection rate 2.27%), and was ascertained to be fully susceptible to DDT.

Conclusion

For the control of malaria, improved diagnosis and sustained supply of drugs for artemisinin-based combination therapy are strongly advocated, which should be enforced for treatment of every single case of P. falciparum. Greater political commitment is called for organized vector control operations along border/high-risk areas to contain the spread of drug-resistant malaria, and averting impending disease outbreaks.

Mitochondrial DNA variation in the malaria vector Anopheles minimus across China, Thailand and Vietnam: evolutionary hypothesis, population structure and population history

The effects of Pleistocene environmental fluctuations on the distribution and diversity of organisms in Southeast Asia are much less well known than in Europe and North America. In these regions, the combination of palaeoenvironmental reconstruction and inferences about population history from genetic data has been very powerful. In Southeast Asia, mosquitoes are good candidates for the genetic approach, with the added benefit that understanding the relative contributions of historical and current processes to population structure can inform management of vector species. Genetic variation among populations of Anopheles minimus was examined using 144 mtDNA COII sequences from 23 sites in China, Thailand and Vietnam. Haplotype diversity was high, with two distinct lineages that have a sequence divergence of over 2% and exhibit different geographical distributions. We compare alternative hypotheses concerning the origin of this pattern. The observed data deviate from the expectations based on a single-panmictic population with or without growth, or a stable but spatially structured population. However, they can be readily accommodated by a model of past fragmentation into eastern and western refugia, followed by growth and range expansion. This is consistent with the palaeoenvironmental reconstructions currently available for the region.

Partial mitochondrial DNA sequences suggest the existence of a cryptic species within the Leucosphyrus group of the genus Anopheles (Diptera: Culicidae), forest malaria vectors, in northern Vietnam

Background

During the last decade, Southeast Asian countries have been very successful in reducing the burden of malaria. However, malaria remains endemic in these countries, especially in remote and forested areas. The Leucosphyrus group of the genus Anopheles harbors the most important malaria vectors in forested areas of Southeast Asia. In Vietnam, previous molecular studies have resulted in the identification of only Anopheles dirus sensu stricto (previously known as An. dirus species A) among the Leucosphyrus group members. However, Vietnamese entomologists have recognized that mosquitoes belonging to the Leucosphyrus group in northern Vietnam exhibit morphological characteristics similar to those of Anopheles takasagoensis, which has been reported only from Taiwan. Here, we aimed to confirm the genetic and morphological identities of the members of the Leucosphyrus group in Vietnam.

Results

In the molecular phylogenetic trees reconstructed using partial COI and ND6 mitochondrial gene sequences, samples collected from southern and central Vietnam clustered together with GenBank sequences of An. dirus that were obtained from Thailand. However, samples from northern Vietnam formed a distinct clade separated from both An. dirus and An. takasagoensis by other valid species.

Conclusions

The results suggest the existence of a cryptic species in northern Vietnam that is morphologically similar to, but phylogenetically distant from both An. dirus and An. takasagoensis. We have tentatively designated this possible cryptic species as Anopheles aff. takasagoensis for convenience, until a valid name is assigned. However, it is difficult to distinguish the species solely on the basis of morphological characteristics. Further studies on such as karyotypes and polytene chromosome banding patterns are necessary to confirm whether An. aff. takasagoensis is a valid species. Moreover, studies on (1) the geographic distribution, which is potentially spreading along the Vietnam, China, Laos, and Myanmar borders; (2) morphological and ecological characteristics; and (3) vectorial capacity of this newly identified cryptic species of An. dirus, which is one of the most important malaria vectors in the mainland of Southeast Asia, are necessary for planning efficient malaria vector control programs in this region.

Inter-specific gene flow dynamics during the Pleistocene-dated speciation of forest-dependent mosquitoes in Southeast Asia

Tropical forests have undergone repeated fragmentation and expansion during Pleistocene glacial and interglacial periods, respectively. The effects of this repeated forest fragmentation in driving vicariance in tropical taxa have been well studied. However, relatively little is known about how often this process results in allopatric speciation, since it may be inhibited by recurrent gene flow during repeated secondary contact, or to what extent Pleistocene-dated speciation results from ecological specialization in the face of gene flow. Here, divergence times and gene flow between three closely-related mosquito species of the Anopheles dirus species complex endemic to the forests of Southeast Asia, are inferred using coalescent based Bayesian analysis. An Isolation with Migration model is applied to sequences of two mitochondrial and three nuclear genes, and 11 microsatellites. The divergence of An. scanloni has occurred despite unidirectional nuclear gene flow from this species into An. dirus. The inferred asymmetric gene flow may result from the unique evolutionary adaptation of An. scanloni to limestone karst habitat, and therefore the fitness advantage of this species over An. dirus in regions of sympatry. Mitochondrial introgression has led to the complete replacement of An. dirus haplotypes with those of An. baimaii through a recent (approximately 62 kya) selective sweep. Speciation of An. baimaii and An. dirus is inferred to have involved allopatric divergence throughout much of the Pleistocene. Secondary contact and bidirectional gene flow has occurred only within the last 100 000 years, by which time the process of allopatric speciation seems to have been largely completed.

Evidence for pleistocene population divergence and expansion of Anopheles albimanus in Southern Central America

The micro-geographic structure of Anopheles albimanus was studied in southern Central America using partial sequences of the mtDNA cytochrome oxidase subunit I gene (COI). Analysis of molecular variance supported significant genetic structure between populations from Costa Rica and western Panama versus those from central-eastern Panama (Phi(CT) = 0.33), whereas the within group divergence was shallow and statistically insignificant (Phi(ST) = 0.08). Furthermore, a statistical parsimony network depicted three divergent groups of haplotypes that were not evenly distributed across the study area. Our findings are in partial agreement with previous studies, yet they do not support physical barriers to gene flow or contemporary isolation by distance in this region. Instead, three co-occurring groups of An. albimanus may be the result of multiple introductions, most likely caused by historical fragmentation and subsequent secondary contact. In addition, the molecular signature of population expansion of An. albimanus was detected in central-eastern Panama approximately 22,000 years ago (95% confidence interval [CI] 10,183-38,169). We hypothesize that the population structure of An. albimanus, as determined by our COI locus analysis, is the result of late Pleistocene climatic changes in northern South America.

Population structure analyses and demographic history of the malaria vector Anopheles albimanus from the Caribbean and the Pacific regions of Colombia

BACKGROUND

Anopheles albimanus is an important malaria vector in some areas throughout its distribution in the Caribbean and the Pacific regions of Colombia, covering three biogeographic zones of the neotropical region, Maracaibo, Magdalena and Chocó.

METHODS

This study was conducted to estimate intra-population genetic diversity, genetic differentiation and demographic history of An. albimanus populations because knowledge of vector population structure is a useful tool to guide malaria control programmes. Analyses were based on mtDNA COI gene sequences and four microsatellite loci of individuals collected in eight populations from the Caribbean and the Pacific regions of Colombia.

RESULTS

Two distinctive groups were consistently detected corresponding to COI haplotypes from each region. A star-shaped statistical parsimony network, significant and unimodal mismatch distribution, and significant negative neutrality tests together suggest a past demographic expansion or a selective sweep in An. albimanus from the Caribbean coast approximately 21,994 years ago during the late Pleistocene. Overall moderate to low genetic differentiation was observed between populations within each region. However, a significant level of differentiation among the populations closer to Buenaventura in the Pacific region was observed. The isolation by distance model best explained genetic differentiation among the Caribbean region localities: Los Achiotes, Santa Rosa de Lima and Moñitos, but it could not explain the genetic differentiation observed between Turbo (Magdalena providence), and the Pacific region localities (Nuquí, Buenaventura, Tumaco). The patterns of differentiation in the populations from the different biogeographic provinces could not be entirely attributed to isolation by distance.

CONCLUSION

The data provide evidence for limited past gene flow between the Caribbean and the Pacific regions, as estimated by mtDNA sequences and current gene flow patterns among An. albimanus populations as measured by MS loci which may be mainly influenced by semi-permeable natural barriers in each biogeographical region that lead to the genetic differences and effective population sizes detected. The relatively high genetic differentiation in the port city of Buenaventura may be the result of specific ecological conditions, human migration and activities and/or differences in effective population sizes. This knowledge could serve to evaluate and coordinate vector control strategies in these regions of Colombia.

Molecular phylogenetics and biogeography of the Neocellia Series of Anopheles mosquitoes in the Oriental Region

Molecular studies of population divergence and speciation across the Oriental Region are sparse, despite the region's high biodiversity and extensive Pliocene and Pleistocene environmental change. A molecular phylogenetic study of the Neocellia Series of Anopheles mosquitoes was undertaken to identify patterns of diversification across the Oriental Region and to infer the role of Pleistocene and Pliocene climatic change. A robust phylogeny was constructed using CO2 and ND5 mitochondrial genes and ITS2 and D3 nuclear ribosomal markers. Bayesian analysis of mitochondrial genes was used to date divergence events. The repeated contraction and expansion of forest habitat resulting from Pleistocene climatic fluctuations appears to have had a substantial impact on intraspecific diversification, but has not driven speciation within this group. Primarily early to mid Pliocene speciation was detected within the Annularis Group, whereas speciation within the Maculatus and Jamesii Groups occurred during the mid and late Pliocene. Both allopatric divergence driven by late Pliocene environmental changes and ecological adaptation, involving altitudinal replacement and seasonality, are likely to have influenced speciation in the Maculatus Group.