Microgeographic genetic variation of the malaria vector Anopheles darlingi root (Diptera: Culicidae) from Cordoba and Antioquia, Colombia

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.

Molecular diversity of genes related to biological rhythms (period and timeless) and insecticide resistance (Na V and ace-1) in Anopheles darlingi

BACKGROUND

Malaria is a public health concern in the Amazonian Region, where Anopheles darlingi is the main vector of Plasmodium spp. Several studies hypothesised the existence of cryptic species in An. darlingi, considering variations in behaviour, morphological and genetic aspects. Determining their overall genetic background for vector competence, insecticide resistance, and other elements is essential to better guide strategies for malaria control.

OBJECTIVES

This study aimed to evaluate the molecular diversity in genes related to behaviour and insecticide resistance, estimating genetic differentiation in An. darlingi populations from Amazonian localities in Brazil and Pacific Colombian region.

METHODS

We amplified, cloned and sequenced fragments of genes related to behaviour: timeless (tim) and period (per), and to insecticide resistance: voltage-gated sodium channel (Na V ) and acetylcholinesterase (ace-1) from 516 An. darlingi DNA samples from Manaus, Unini River, Jaú River and Porto Velho - Brazil, and Chocó - Colombia. We discriminated single nucleotide polymorphisms (SNPs), determined haplotypes and evaluate the phylogenetic relationship among the populations.

FINDINGS

The genes per, tim and ace-1 were more polymorphic than Na V . The classical kdr and ace-1 R mutations were not observed. Phylogenetic analyses suggested a significant differentiation between An. darlingi populations from Brazil and Colombia, except for the Na V gene. There was a geographic differentiation within Brazilian populations considering per and ace-1.

CONCLUSIONS

Our results add genetic data to the discussion about polymorphisms at population levels in An. darlingi. The search for insecticide resistance-related mechanisms should be extended to more populations, especially from localities with a vector control failure scenario.

Humidity - The overlooked variable in the thermal biology of mosquito-borne disease

Vector-borne diseases cause significant financial and human loss, with billions of dollars spent on control. Arthropod vectors experience a complex suite of environmental factors that affect fitness, population growth and species interactions across multiple spatial and temporal scales. Temperature and water availability are two of the most important abiotic variables influencing their distributions and abundances. While extensive research on temperature exists, the influence of humidity on vector and pathogen parameters affecting disease dynamics are less understood. Humidity is often underemphasized, and when considered, is often treated as independent of temperature even though desiccation likely contributes to declines in trait performance at warmer temperatures. This Perspectives explores how humidity shapes the thermal performance of mosquito-borne pathogen transmission. We summarize what is known about its effects and propose a conceptual model for how temperature and humidity interact to shape the range of temperatures across which mosquitoes persist and achieve high transmission potential. We discuss how failing to account for these interactions hinders efforts to forecast transmission dynamics and respond to epidemics of mosquito-borne infections. We outline future research areas that will ground the effects of humidity on the thermal biology of pathogen transmission in a theoretical and empirical framework to improve spatial and temporal prediction of vector-borne pathogen transmission.

Dispersion routes of the main vectors of human malaria in the Americas: Genetic evidence from the mitochondrial COI gene

BACKGROUND AND OBJECTIVES

In America, of the 44 species of Anopheles, nine are main vectors of malaria and, of these, genetic information exists for seven. Hence, this study sought to know the gene flow and diversity of the seven principal vectors of malaria at the Americas level.

METHODS

For the seven species and the sequences of the mitochondrial cytochrome c oxidase I (COI) gene obtained from the GenBank and Bold System, genetic analyzes of populations and genetic structure were performed and haplotype networks and phylogenetic trees were obtained.

RESULTS

For the seven species, 1440 sequences were analyzed and 519 haplotypes were detected. The Hd and π values were higher within a continental context than by countries. Neutrality tests indicated positive and negative values with most of these being significant (p < 0.05). Phylogenetic analyses for all the species recovered three clades with no geographic pattern among them.

INTERPRETATION & CONCLUSION

Studies suggest that native species of Anopheles from the Americas have greater haplotype diversity and low genetic differentiation due to the lack of physical barriers to impede gene flow among these populations. Moreover, all the species are interconnected by roadways. This scenario complicates the epidemiological picture of malaria in the Americas.

Population genetic structure of the malaria vector Anopheles fluviatilis species T (Diptera: Culicidae) in India

Anopheles fluviatilis James (Diptera: Culicidae) represents a complex that comprises four sibling species (S, T, U, and V). Among these, species T is widely distributed in India. Chromosomal inversion polymorphism exists among different geographic populations of An. fluviatilis species T; however, population genetic structure is not understood. This study inferred a genetic structure among six geographically diverse populations of species T using a panel of microsatellite markers. Analyses indicated a significant but low genetic differentiation among the majority of the studied populations. A significant correlation was observed between genetic and geographic distances, exhibiting stepwise migration patterns among populations.

Genetic and morphological polymorphisms of Aedes scapularis (Diptera: Culicidae), vector of filariae and arboviruses

BACKGROUND

Aedes scapularis is a neotropical mosquito that is competent to vector viruses and filariae. It is reputed to be highly morphologically and genetically polymorphic, facts that have raised questions about whether it is a single taxonomic entity. In the last five decades, authors have posed the hypothesis that it could actually be a species complex under incipient speciation. Due to its epidemiological importance, its taxonomic status should be determined with confidence.

AIM AND METHOD

Our objective was to investigate more deeply the polymorphism of Ae. scapularis to detect any evidence of incipient speciation of cryptic species. We then compared populational samples from the Southeastern, Northern and Northeastern regions of Brazil. The biological markers used in the comparison were: the complete mitochondrial DNA, the isolated mitochondrial gene cytochrome oxidase subunit I (COI) and wing geometry.

RESULTS AND DISCUSSION

As expected, high morphological/genetic polymorphism was observed in all Ae. scapularis populations, however it was not indicative of segregation or incipient speciation. There was no correlation between wing shape and the geographical origin of the populations analysed. A congruent observation resulted from the analysis of the COI gene, which revealed a high number of haplotypes (51) and no clusterization of populational samples according to the original biomes. In the phylogenetic analysis of the 13 mitochondrial protein-coding genes, the Ae. scapularis clade clustered with maximum support (100% bootstrap support and posterior probability of 1). No significant internal structure was observed in the Ae. scapularis clade, which was nearly a polytomy. Taken together, our results indicate that this species is not a species complex.

CONCLUSION

We conclude that there was no indication, in the analysed regions, of the occurrence of more than one taxon in the species Ae. scapularis, despite it being highly polymorphic. By ruling out the former species complex hypothesis, our phylogenetic results reinforce that Ae. scapularis is a single taxonomic unit and should be monitored with standardized surveillance and control methods.

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.

Natural malaria infection in anophelines vectors and their incrimination in local malaria transmission in Darién, Panama

Background

More than 85% of the malaria cases in Panama occur in poor, rural and indigenous regions like Darien Province. Vector diversity, infection rate and spatial distribution are important entomological parameters of malaria transmission dynamics. Their understanding is crucial for the development of effective disease control strategies. The objective of this study was to determine the composition of Anopheles species, their natural infection rate and their geographic distribution to better understand the malaria transmission dynamics in Darién, Panama.

Methods

Anophelines mosquitoes were captured during the rainy and dry season of 2016. We selected five communities where adult anophelines were collected using CDC light-traps, and through protective human-baited traps. Detection of natural infection and Plasmodium genotype were detected via nested PCR through the amplification of ssrRNA and the circumsporozoite protein gene (csp), respectively.

Results

A total of 1,063 mosquitoes were collected mosquitoes were collected for the detection of natural infection with Plasmodium spp. Nine Anophelines species were identified, with the predominant species being: An. (Nys.) darlingi (45.0%) and An. (Nys.) albimanus (42.6%). Natural infection in An. (Nys.) albimanus with P. vivax was detected in one mosquito pool from the community Pueblo Tortuga (0.6%), three from Marraganti (1.7%), two from Bajo Chiquito (1.1%) and three pools from Alto Playona 3 (1.7%). For An. (Nys.) darlingi mosquitoes, we detected seven positive pools from the community Bajo Chiquito (4.0%), two pools from Marraganti (1.1%) and two pools from Alto Playona (1.1%). The P. vivax allelic variant VK210 was detected in infected mosquitoes.

Conclusion

The results from this study provide new information on the transmission dynamics associated with anophelines vectors in the Darién region. This is the first report of natural P. vivax infection in An. (Nys.) darlingi and its incrimination as a potential malaria vector in this region of Panama. Additional studies are necessary to expand our knowledge and determine crucial parameters in malaria transmission in Darién, which in turn will aid the National Malaria Program in attaining an adequate malaria control strategy towards malaria elimination.

Antibody Responses Against Anopheles darlingi Immunogenic Peptides in Plasmodium Infected Humans

Introduction: Malaria is still an important vector-borne disease in the New World tropics. Despite the recent decline in malaria due to Plasmodium falciparum infection in Africa, a rise in Plasmodium infections has been detected in several low malaria transmission areas in Latin America. One of the main obstacles in the battle against malaria is the lack of innovative tools to assess malaria transmission risk, and the behavioral plasticity of one of the main malaria vectors in Latin America, Anopheles darlingi.

Methods: We used human IgG antibodies against mosquito salivary gland proteins as a measure of disease risk. Whole salivary gland antigen (SGA) from Anopheles darlingi mosquitoes was used as antigen in Western blot experiments, in which a ~65 kDa protein was visualized as the main immunogenic band and sent for sequencing by mass spectrometry. Apyrase and peroxidase peptides were designed and used as antigens in an ELISA-based test to measure human IgG antibody responses in people with different clinical presentations of malaria.

Results: Liquid chromatography–mass spectrometry revealed 17 proteins contained in the ~65 kDa band, with an apyrase and a peroxidase as the two most abundant proteins. Detection of IgG antibodies against salivary antigens by ELISA revealed a significant higher antibody levels in people with malaria infection when compared to uninfected volunteers using the AnDar_Apy1 and AnDar_Apy2 peptides. We also detected a significant positive correlation between the anti-peptides IgG levels and antibodies against the Plasmodium vivax and P. falciparum antigens PvMSP1 and PfMSP1. Odd ratios suggest that people with higher IgG antibodies against the apyrase peptides were up to five times more likely to have a malaria infection.

Conclusion: Antibodies against salivary peptides from An. darlingi salivary gland proteins may be used as biomarkers for malaria risk.

Identification and Pilot Evaluation of Salivary Peptides from Anopheles albimanus as Biomarkers for Bite Exposure and Malaria Infection in Colombia

Insect saliva induces significant antibody responses associated with the intensity of exposure to bites and the risk of disease in humans. Several salivary biomarkers have been characterized to determine exposure intensity to Old World Anopheles mosquito species. However, new tools are needed to quantify the intensity of human exposure to Anopheles bites and understand the risk of malaria in low-transmission areas in the Americas. To address this need, we conducted proteomic and bioinformatic analyses of immunogenic candidate proteins present in the saliva of uninfected Anopheles albimanus from two separate colonies—one originating from Central America (STECLA strain) and one originating from South America (Cartagena strain). A ~65 kDa band was identified by IgG antibodies in serum samples from healthy volunteers living in a malaria endemic area in Colombia, and a total of five peptides were designed from the sequences of two immunogenic candidate proteins that were shared by both strains. ELISA-based testing of human IgG antibody levels against the peptides revealed that the transferrin-derived peptides, TRANS-P1, TRANS-P2 and a salivary peroxidase peptide (PEROX-P3) were able to distinguish between malaria-infected and uninfected groups. Interestingly, IgG antibody levels against PEROX-P3 were significantly lower in people that have never experienced malaria, suggesting that it may be a good marker for mosquito bite exposure in naïve populations such as travelers and deployed military personnel. In addition, the strength of the differences in the IgG levels against the peptides varied according to location, suggesting that the peptides may able to detect differences in intensities of bite exposure according to the mosquito population density. Thus, the An. albimanus salivary peptides TRANS-P1, TRANS-P2, and PEROX-P3 are promising biomarkers that could be exploited in a quantitative immunoassay for determination of human-vector contact and calculation of disease risk.

Behavior and abundance of Anopheles darlingi in communities living in the Colombian Amazon riverside

In the past few years, relative frequencies of malaria parasite species in communities living in the Colombian Amazon riverside have changed, being Plasmodium vivax (61.4%) and Plasmodium malariae (43.8%) the most frequent. Given this epidemiological scenario, it is important to determine the species of anophelines involved in these parasites' transmission. This study was carried out in June 2016 in two indigenous communities living close to the tributaries of the Amazon River using protected human bait. The results of this study showed a total abundance of 1,085 mosquitos, of which 99.2% corresponded to Anopheles darlingi. Additionally, only two anopheline species were found, showing low diversity in the study areas. Molecular confirmation of some individuals was then followed by evolutionary analysis by using the COI gene. Nested PCR was used for identifying the three Plasmodium species circulating in the study areas. Of the two species collected in this study, 21.0% of the An. darlingi mosquitoes were infected with P. malariae, 21.9% with P. vivax and 10.3% with Plasmodium falciparum. It exhibited exophilic and exophagic behavior in both study areas, having marked differences regarding its abundance in each community (Tipisca first sampling 49.4%, Tipisca second sampling 39.6% and Doce de Octubre 10.9%). Interestingly, An. mattogrossensis infected by P. vivax was found for the first time in Colombia (in 50% of the four females collected). Analysis of An. darlingi COI gene diversity indicated a single population maintaining a high gene flow between the study areas. The An. darlingi behavior pattern found in both communities represents a risk factor for the region's inhabitants living/working near these sites. This highlights the need for vector control efforts such as the use of personal repellents and insecticides for use on cattle, which must be made available in order to reduce this Anopheline's abundance.

Genetic diversity and population structure of Anopheles triannulatus s. l. in the department of Córdoba, Colombia, using DNA barcoding

Introduction: Anopheles triannulatus is not incriminated as a vector of malaria transmission in Colombia despite recent reports of infection with Plasmodium spp. in populations related to the northwestern and southeastern lineages. Genetic diversity can delimit information about gene flow and population differentiation in localities with malaria. Objective: To estimate the genetic diversity of An. triannulatus in five municipalities with high and low incidence of malaria in the department of Córdoba. Materials and methods: The entomological collections were done between August and November, 2016, in Tierralta, Puerto Libertador, Montelíbano, Sahagún, and Planeta Rica. We used the COI barcoding fragment as molecular marker. The genetic analysis included the estimation of genetic parameters such as the diversity haplotype, the genetic structure, the gene flow, the Tajima’s D test, the haplotype network, and the phylogenetic relationship. Results: We obtained 148 sequences with a length of 655 nucleotides of the COI gene, from which we derived 44 haplotypes. The H2 and H21 haplotypes were the most frequent in the populations. The values of the Tajima’s D test were negative and not significant (p>0.10). The genetic structure index (FST=0.01427) and the gene flow (Nm=17.27) evidenced no differentiation between sampled populations due to the high exchange of migrants. Using phylogenetic inferences and the haplotype network, we identified one single species without geographic differentiation or lineages in the geographic range studied. Conclusions: The genetic diversity calculated for An. triannulatus in this context indicated stable populations in constant exchange.

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.

Genetic Differentiation of Colombian Populations of Anopheles darlingi Root (Diptera: Culicidae)

Anopheles darlingi Root is a primary vector of malaria in the neotropic region, a species not just highly anthropophilic but very efficient in transmitting Plasmodium species and considered the most important vector in the Amazon region. The main goal of this study was to determine the genetic structure of the A. darlingi populations using microsatellites (STR) in western and eastern regions of Colombia. DNA extraction was done with the cited protocol of band using the Genomic Prep™ cell and tissue isolation commercial kits. We used the STR reported by Conn et al (Mol Ecol Notes 1: 223-225, 2001). The analysis with STR proved there was a high genetic diversity and significant alterations of the Hardy-Weinberg equilibrium. The greatest genetic diversity was recorded in Mitu (Vaupes) (Na = 14, Ho = 0.520). The lowest was in Pueblo Nuevo (Cordoba) (Na = 12, Ho = 0.457). The eastern region and the Mitu (Vaupes) populations presented the highest number of primer alleles (Ap = 30; Ap = 13; Ap = 9), with variations between 0.010 and 0.097. The AMOVA revealed that the whole population underwent moderate genetic differentiation (F _ST = 0.063, p < 0.05). The same differentiation was noticed (0.06 < F _ST > 0.06, p < 0.05) with five of the six populations included in this job, and there was a low differentiation in the Las Margaritas (Santander) area (F _ST = 0.02s3, p < 0.05). Our results suggest a slight positive correlation, which does not show a statistical significance between the geographic and genetic distances, probably suggesting that the moderate genetic differentiation found between pairs of populations does not need to be explained for the hypothesis of separation by distance.

[Anopheles darlingi (Diptera: Culicidae) Rood 1926: Morphometric variations in wings and legs of populations from Colombia]

INTRODUCTION

Natural populations of Anopheles darlingi, the main malaria vector in Colombia, have shown phenotypic variations in some of their diagnostic characters.

OBJECTIVE

To characterize morphometric variations in wing spot patterns and diagnostic characters of the hind leg of adult females of An. darlingi collected in areas where malaria is endemic.

MATERIALS AND METHODS

Patterns of the costal vein spots of wild females of An. darling populations collected in the departments of Chocó, Guaviare, Meta and Vichada were analyzed using linear and morphometric approaches. The second tarsomere of the hind leg of females was characterized by linear morphometric analysis.

RESULTS

We found 19 patterns of spots in the costal vein; patterns I (n=118/240, 49%) and VI (n=66, 28%) were the most frequent. The proportion of the basal dark area of hind tarsomere II and the length of hind tarsomere II (DSIII2/Ta-III2) constituted a robust diagnostic character as it represented 89% (n=213/240) of the total specimens studied. Significant differences were found in the wing shape (F=1.65, df =50, p<0.001) and the wing size (F=3.37, df=5, p=0.005) among populations from different locations. The smallest centroid size (2.64 mm) was found in populations from Chocó.

CONCLUSIONS

We registered 11 new wing spot patterns in the costal vein and the dominance of the patterns I and VI for populations of An. darlingi from Colombia. We confirmed DSIII2/TaIII2 ratio as a robust diagnostic character for the taxonomy of this species. We found differences between the size and shape of the wings of An. darlingi populations in accordance to their geographical distribution, which constitute important bionomic aspects for this malaria vector.

Genetic diversity and population structure of the primary malaria vector Anopheles sinensis (Diptera: Culicidae) in China inferred by cox1 gene

Background

Anopheles sinensis is a primary vector for Plasmodium vivax malaria in most regions of China. A comprehensive understanding of genetic variation and structure of the mosquito would be of benefit to the vector control and in a further attempt to contribute to malaria elimination in China. However, there is only inadequate population genetic data pertaining to An. sinensis currently.

Methods

Genetic variations and structure among populations of An. sinensis was examined and analyzed based on the nucleotide sequences of a 662 nt variable region of the mitochondrial cox1 gene among 15 populations from 20 collection sites in China.

Results

A total of 453 individuals in 15 populations were analyzed. The cox1 gene sequences were aligned, and 247 haplotypes were detected, 41 of these shared between populations. The range of haplotype diversity was from 0.709 (Yunnan) to 0.998 (Anhui). The genealogic network showed that the haplotypes were divided into two clusters, cluster I was at a high level of homoplasy, while cluster II included almost all individuals from the Yunnan population. The Yunnan population displayed a significantly high level of genetic differentiation (0.452−0.622) and a restricted gene flow with other populations. The pairwise F_ST values among other populations were lower. The AMOVA result showed that the percentage of variation within populations (83.83%) was higher than that among populations (16.17%). Mantel test suggested that geographical distance did not significantly contribute to the genetic differentiation (R² = 0.0125, P = 0.59). Neutral test and mismatch analysis results showed that the An. sinensis population has undergone demographic expansions.

Conclusions

Anopheles sinensis populations showed high genetic polymorphism by cox1 gene. The weak genetic structure may be a consequence of low genetic differentiation and high gene flow among populations, except the Yunnan samples. The Yunnan population was isolated from the other populations, gene flow limited by geographical distance and barriers. These findings will provide a theoretical basis for vector surveillance and vector control in China.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-017-2013-z) contains supplementary material, which is available to authorized users.

Use of DNA barcoding to distinguish the malaria vector Anopheles neivai in Colombia

A reference 535 bp barcode sequence from a fragment of the mitochondrial gene cytochrome oxidase I (COI), acquired from specimens of An. neivai Howard, Dyar & Knab, 1913 from its type locality in Panama, was used as a tool for distinguishing this species from others in the subgenus Kerteszia. Comparisons with corresponding regions of COI between An. neivai and other species in the subgenus (An. bellator Dyar & Knab 1906, An. homunculus Komp 1937, An cruzii Dyar & Knab, 1908 and An. laneanus Corrêa & Cerqueira, 1944) produced K2P genetic distances of 8.3-12.6%, values well above those associated with intraspecific variation. In contrast, genetic distances among 55 specimens from five municipalities in the Colombian Pacific coastal state of Chocó were all within the range of 0-2.5%, with an optimized barcode threshold of 1.3%, the limit for unambiguous differentiation of An. neivai. Among specimens from the Chocó region, 18 haplotypes were detected, two of which were widely distributed over the municipalities sampled. The barcode sequence permits discrimination of An. neivai from sympatric species and indicates genetic variability within the species; aspects key to malaria surveillance and control as well as defining geographic distribution and dispersion patterns.

Population dynamics of Anopheles nuneztovari in Colombia

Anopheles nuneztovari is an important Colombian malaria vector widespread on both sides of the Andean Mountains, presenting morphological, behavioral and genetic heterogeneity throughout the country. The aim of this study was to evaluate whether the population structure and distribution of An. nuneztovari in Colombia are associated with ecological and physical barriers present in a heterogeneous landscape. Further, differences in behavior were addressed. A total of 5392 specimens of An. nuneztovari were collected. Mitochondrial and nuclear marker analyses detected subdivision among the northwest-west, northeast and east populations. For both markers, isolation by distance (~53%) and isolation by resistance (>30%) were determinants of population genetic differentiation. This suggests that physical barriers, geographical distance and ecological differences on both sides of the Andean Mountains promoted the genetic differentiation and population subdivision of An. nuneztovari in Colombia. This species showed the highest biting activity after 20:00h; indoor and outdoor preferences were found in all localities. These results indicated that the most effective interventions for controlling vector populations on both sides of the Andes need to be region-specific.

Behavior and population structure of Anopheles darlingi in Colombia

Anopheles darlingi is a widely distributed and important malaria vector in Colombia. Biogeographical and ecological heterogeneity across the Colombian distribution led to the hypothesis of behavioral and genetic differentiation among A. darlingi populations. A total of 2017 A. darlingi specimens were collected during 222 h of sampling. This vector was the most abundant anopheline species in most of the localities sampled. Subdivision between samples collected west and east of the Andes was indicated by 1) mitochondrial COI and nuclear CAD sequences from NW-W and CE-S populations (COI ΦST=0.48761-0.81974, CAD FST=0.11319-0.21321), 2) a COI haplotype network, and 3) SAMOVA. Endo- and exophagy were detected in populations west of the Andes, whereas exophagy was evident in PTG, a locality east of the Andes. Isolation by resistance was significant for COI and explained 26% of the genetic differentiation. We suggest that at a macrogeographic scale, the Andes influence the differentiation of A. darlingi in Colombia and may drive divergence, and, at a microgeographic scale, ecological differences have a significant impact on structure. These data could constitute a baseline for the design of effective vector interventions, locality-specific for the east and similar for panmictic populations west of the Andes.

Evidence for temporal population replacement and the signature of ecological adaptation in a major Neotropical malaria vector in Amazonian Peru

BACKGROUND

The major Neotropical malaria vector, Anopheles darlingi, was reintroduced into the Iquitos, Loreto, Peru area during the early 1990s, where it displaced other anophelines and caused a major malaria epidemic. Since then, case numbers in Loreto have fluctuated, but annual increases have been reported since 2012.

METHODS

The population genetic structure of An. darlingi sampled before and after the introduction of long-lasting insecticidal nets (LLINs) was investigated to test the hypothesis of temporal population change (2006 vs. 2012). Current samples of An. darlingi were used to test the hypothesis of ecological adaptation to human modified (highway) compared with wild (riverine) habitat, linked to forest cover. In total, 693 An. darlingi from nine localities in Loreto, Peru area were genotyped using 13 microsatellite loci. To test the hypothesis of habitat differentiation in An. darlingi biting time patterns, HBR and EIR, four collections of An. darlingi from five localities (two riverine and three highway) were analysed.

RESULTS

Analyses of microsatellite loci from seven (2006) and nine settlements (2012-2014) in the Iquitos area detected two distinctive populations with little overlap, although it is unclear whether this population replacement event is associated with LLIN distribution or climate. Within the 2012-2014 population two admixed subpopulations, A and B, were differentiated by habitat, with B significantly overrepresented in highway, and both in near-equal proportions in riverine. Both subpopulations had a signature of expansion and there was moderate genetic differentiation between them. Habitat and forest cover level had significant effects on HBR, such that Plasmodium transmission risk, as measured by EIR, in peridomestic riverine settlements was threefold higher than in peridomestic highway settlements. HBR was directly associated with available host biomass rather than forest cover.

CONCLUSIONS

A population replacement event occurred between 2006 and 2012-2014, concurrently with LLIN distribution and a moderate El Niño event, and prior to an increase in malaria incidence. The likely drivers of this replacement cannot be determined with current data. The present-day An. darlingi population is composed of two highly admixed subpopulations, which appear to be in an early stage of differentiation, triggered by anthropogenic alterations to local habitat.

Brazilian Anopheles darlingi Root (Diptera: Culicidae) Clusters by Major Biogeographical Region

The major drivers of the extensive biodiversity of the Neotropics are proposed to be geological and tectonic events together with Pliocene and Pleistocene environmental and climatic change. Geographical barriers represented by the rivers Amazonas/Solimões, the Andes and the coastal mountain ranges in eastern Brazil have been hypothesized to lead to diversification within the primary malaria vector, Anopheles (Nyssorhynchus) darlingi Root, which primarily inhabits rainforest. To test this biogeographical hypothesis, we analyzed 786 single nucleotide polymorphisms (SNPs) in 12 populations of An. darlingi from across the complex Brazilian landscape. Both model-based (STRUCTURE) and non-model-based (Principal Components and Discriminant Analysis) analysis of population structure detected three major genetic clusters that correspond with newly described Neotropical biogeographical regions: 1) Atlantic Forest province (= southeast population); 2) Parana Forest province (= West Atlantic forest population, with one Chacoan population - SP); and 3) Brazilian dominion population (= Amazonian population with one Chacoan population - TO). Significant levels of pairwise genetic divergences were found among the three clusters, allele sharing among clusters was negligible, and geographical distance did not contribute to differentiation. We infer that the Atlantic forest coastal mountain range limited dispersal between the Atlantic Forest province and the Parana Forest province populations, and that the large, diagonal open vegetation region of the Chacoan dominion dramatically reduced dispersal between the Parana and Brazilian dominion populations. We hypothesize that the three genetic clusters may represent three putative species.

Anopheles species composition explains differences in Plasmodium transmission in La Guajira, northern Colombia

Malaria in La Guajira, the most northern state of Colombia, shows two different epidemiological patterns. Malaria is endemic in the municipality of Dibulla whereas in Riohacha it is characterised by sporadic outbreaks. This study aimed to establish whether differences in transmission patterns could be attributed to different vector species. The most abundant adult female species were Anopheles aquasalis, exclusive to Riohacha, and Anopheles darlingi, restricted to Dibulla. Anopheles mosquitoes were identified using morphology and the molecular markers internal transcribed spacer 2 and cytochrome c oxidase I. All specimens (n = 1,393) were tested by ELISA to determine natural infection rates with Plasmodium falciparum and Plasmodium vivax. An. darlingi was positive for P. vivax 210, with an infection rate of 0.355% and an entomological inoculation rate of 15.87 infective bites/person/year. Anopheles albimanus larvae were the most common species in Riohacha, found in temporary swamps; in contrast, in Dibulla An. darlingi were detected mainly in permanent streams. Distinctive species composition and larval habitats in each municipality may explain the differences in Plasmodium transmission and suggest different local strategies should be used for vector control.

High morphological and genetic variabilities of Ochlerotatus scapularis, a potential vector of filarias and arboviruses

Background

Ochlerotatus scapularis is a potential vector of filarias and arboviruses in the Neotropics. This species was once typically associated with sylvatic environments; however, cases of synanthropy and urbanization of this species have been increasingly reported in southeast Brazil. Despite the medical relevance of Oc. scapularis, its populational variability is not yet known. To our knowledge, this is the first report describing the morphological and genetic variabilities of this species.

Methods

Population samples were characterized using the cytochrome oxidase subunit I (COI) mitochondrial gene and wing geometrics. Adult mosquitoes were collected from five sampling sites from remnants of the Atlantic forest embedded in the urban or rural areas of southeast Brazil.

Results

In the 130 individuals analyzed, 46 COI haplotypes were detected. Haplotype diversity was high and ranged from 0.66 to 0.97. Six haplotypes were present in 61% of the individuals, whereas the remaining haplotypes were less frequent (39%). Wing shape was also highly polymorphic. Differentiation of populations across sampling sites according to genetic distances (F_st = −0.009 to 0.060) and morphological distances (Q_st = 0.47) indicated that populations were not identical. No correlations were noted for phenetic and genetic diversities (p = 0.19) or for genetic or phenetic distances with geographical distances (p = 0.2 and p = 0.18, respectively).

Conclusions

Our study results suggest that Oc. scapularis has a rich genetic patrimony, even though its habitat is fragmented. Implications of such genetic richness with respect to vectorial competence, plasticity, and ability to exploit urbanized areas need to be further investigated.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-015-0740-6) contains supplementary material, which is available to authorized users.

Geometric morphometric analysis of Colombian Anopheles albimanus (Diptera: Culicidae) reveals significant effect of environmental factors on wing traits and presence of a metapopulation

Anopheles albimanus is a major malaria mosquito vector in Colombia. In the present study, wing variability (size and shape) in An. albimanus populations from Colombian Maracaibo and Chocó bio-geographical eco-regions and the relationship of these phenotypic traits with environmental factors were evaluated. Microsatellite and morphometric data facilitated a comparison of the genetic and phenetic structure of this species. Wing size was influenced by elevation and relative humidity, whereas wing shape was affected by these two variables and also by rainfall, latitude, temperature and eco-region. Significant differences in mean shape between populations and eco-regions were detected, but they were smaller than those at the intra-population level. Correct assignment based on wing shape was low at the population level (<58%) and only slightly higher (>70%) at the eco-regional level, supporting the low population structure inferred from microsatellite data. Wing size was similar among populations with no significant differences between eco-regions. Population relationships in the genetic tree did not agree with those from the morphometric data; however, both datasets consistently reinforced a panmictic population of An. albimanus. Overall, site-specific population differentiation is not strongly supported by wing traits or genotypic data. We hypothesize that the metapopulation structure of An. albimanus throughout these Colombian eco-regions is favoring plasticity in wing traits, a relevant characteristic of species living under variable environmental conditions and colonizing new habitats.

Temporal genetic structure of major dengue vector Aedes aegypti from Manaus, Amazonas, Brazil

In recent years, high levels of Aedes aegypti infestation and several dengue outbreaks with fatal outcome cases have been reported in Manaus, State of Amazonas, Brazil. This situation made it important to understand the genetic structure and gene flow patterns among the populations of this vector in Manaus, vital pieces of information for their management and development of new control strategies. In this study, we used nine microsatellite loci to examine the effect of seasonality on the genetic structure and gene flow patterns in Ae. aegypti populations from four urban neighborhoods of Manaus, collected during the two main rainy and dry seasons. All loci were polymorphic in the eight samples from the two seasons, with a total of 41 alleles. The genetic structure analyses of the samples from the rainy season revealed genetic homogeneity and extensive gene flow, a result consistent with the abundance of breeding sites for this vector. However, the samples from the dry season were significantly structured, due to a reduction of Ne in two (Praça 14 de Janeiro and Cidade Nova) of the four samples analyzed, and this was the primary factor influencing structure during the dry season. Genetic bottleneck analyses suggested that the Ae. aegypti populations from Manaus are being maintained continuously throughout the year, with seasonal reduction rather than severe bottleneck or extinction, corroborating previous reports. These findings are of extremely great importance for designing new dengue control strategies in Manaus.

Phylogeography of the neotropical Anopheles triannulatus complex (Diptera: Culicidae) supports deep structure and complex patterns

BACKGROUND

The molecular phylogenetic relationships and population structure of the species of the Anopheles triannulatus complex: Anopheles triannulatus s.s., Anopheles halophylus and the putative species Anopheles triannulatus C were investigated.

METHODS

The mitochondrial COI gene, the nuclear white gene and rDNA ITS2 of samples that include the known geographic distribution of these taxa were analyzed. Phylogenetic analyses were performed using Bayesian inference, Maximum parsimony and Maximum likelihood approaches.

RESULTS

Each data set analyzed septely yielded a different topology but none provided evidence for the seption of An. halophylus and An. triannulatus C, consistent with the hypothesis that the two are undergoing incipient speciation. The phylogenetic analyses of the white gene found three main clades, whereas the statistical parsimony network detected only a single metapopulation of Anopheles triannulatus s.l. Seven COI lineages were detected by phylogenetic and network analysis. In contrast, the network, but not the phylogenetic analyses, strongly supported three ITS2 groups. Combined data analyses provided the best resolution of the trees, with two major clades, Amazonian (clade I) and trans-Andean + Amazon Delta (clade II). Clade I consists of multiple subclades: An. halophylus + An. triannulatus C; trans-Andean Venezuela; central Amazonia + central Bolivia; Atlantic coastal lowland; and Amazon delta. Clade II includes three subclades: Panama; cis-Andean Colombia; and cis-Venezuela. The Amazon delta specimens are in both clades, likely indicating local sympatry. Spatial and molecular variance analyses detected nine groups, corroborating some of subclades obtained in the combined data analysis.

CONCLUSION

Combination of the three molecular markers provided the best resolution for differentiation within An. triannulatus s.s. and An. halophylus and C. The latest two species seem to be very closely related and the analyses performed were not conclusive regarding species differentiation. Further studies including new molecular markers would be desirable to solve this species status question. Besides, results of the study indicate a trans-Andean origin for An. triannulatus s.l. The potential implications for malaria epidemiology remain to be investigated.

Molecular evidence for a single taxon, Anopheles nuneztovari s.l., from two endemic malaria regions in Colombia

To elucidate the Anopheles nuneztovari s.l. taxonomic status at a microgeographic level in four malaria endemic localities from Antioquia and Córdoba, Colombia, fragments of the cytochrome oxidase subunit I (COI) and the white gene were used. The COI analysis showed low genetic differentiation with fixation index (F(ST)) levels between -0.02-0.137 and Nm values between 3-∞, indicating the presence of high gene flow among An. nuneztovari s.l. populations from the four localities. The COI network showed a single most common haplotype, type 1 (n = 55), present in all localities, as the likely ancestral haplotype. Analysis of the white gene showed that An. nuneztovari s.l. populations from both departments grouped with haplotypes 19 and 20, which are part of lineage 3 reported previously. The results of the present study suggest that An. nuneztovari s.l. is a single taxon in the area of the present study.

Genetic diversity of Anopheles triannulatus s.l. (Diptera: Culicidae) from northwestern and southeastern Colombia

Anopheles triannulatus s.l. is a species complex, however in Colombia its taxonomic status is unclear. This study was conducted to understand the level of genetic differentiation or population structure of specimens of An. triannulatus s.l. from northwestern and southeastern Colombia. Cytochrome oxidase subunit I (COI) and internal transcribed spacer (ITS2) sequence analyses suggested high genetic differentiation between the NW and SE populations. A TCS network and Bayesian inference analysis based on 814 bp of COI showed two main groups: group I included samples from the NW and group II samples from the SE. Two main ITS2-polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) patterns were found. Pattern I is present in both the NW and SE, and pattern II is found in the SE specimens. To further elucidate the taxonomic status of An. triannulatus s.l. in Colombia and how these COI lineages are related to the Triannulatus Complex species, the evaluation of immature stages, male genitalia, and additional mitochondrial and nuclear markers will be needed.

Ecology of Anopheles darlingi Root with respect to vector importance: a review

Anopheles darlingi is one of the most important malaria vectors in the Americas. In this era of new tools and strategies for malaria and vector control it is essential to have knowledge on the ecology and behavior of vectors in order to evaluate appropriateness and impact of control measures. This paper aims to provide information on the importance, ecology and behavior of An. darlingi. It reviews publications that addressed ecological and behavioral aspects that are important to understand the role and importance of An. darlingi in the transmission of malaria throughout its area of distribution. The results show that Anopheles darlingi is especially important for malaria transmission in the Amazon region. Although numerous studies exist, many aspects determining the vectorial capacity of An. darlingi, i.e. its relation to seasons and environmental conditions, its gonotrophic cycle and longevity, and its feeding behavior and biting preferences, are still unknown. The vector shows a high degree of variability in behavioral traits. This makes it difficult to predict the impact of ongoing changes in the environment on the mosquito populations. Recent studies indicate a good ability of An. darlingi to adapt to environments modified by human development. This allows the vector to establish populations in areas where it previously did not exist or had been controlled to date. The behavioral variability of the vector, its adaptability, and our limited knowledge of these impede the establishment of effective control strategies. Increasing our knowledge of An. darlingi is necessary.

Mitochondrial genetic differentiation across populations of the malaria vector Anopheles lesteri from China (Diptera: Culicidae)

Background

Anopheles lesteri is a primary vector of Plasmodium spp. in central China. A complete understanding of vector population structure and the processes responsible for the differentiation is important to the vector-based malaria control programmes and for identifying heterogeneity in disease transmission as a result of discrete vector populations. There is no adequate An. lesteri population genetic data available.

Methods

Polymorphism of sequence variations in mitochondrial COII and Cytb genes were assessed to explore the level of genetic variability and differentiation among six populations of An. lesteri from China.

Results

There were 30 (4.37%) and 21 (5.33%) polymorphic sites for mtDNA-COII and Cytb gene, respectively. Totally 31 COII and 30 Cytb haplotypes were obtained. The range of F_ST values was from 0.101 to 0.655 by mtDNA-COII, and 0.029 to 0.231 by Cytb gene. The analysis of molecular variance (AMOVA) showed that the percentage of variation within populations (65.83%, 88.48%) was greater than that among populations (34.17%, 11.52%) using both genes. The Tajima's D and Fu's Fs values were all negative, except Tajima's D values of YN and HNB populations, which suggest a large number of low-frequency mutations in populations and the populations were in expansion proceeding.

Conclusions

Levels of genetic variation within An. lesteri populations were higher than among them. While these results may suggest considerable levels of gene flow, other explanations, such as the effect of historical population perturbations can also be hypothesized.

Malaria vector species in Colombia: a review

Here we present a comprehensive review of the literature on the vectorial importance of the major Anopheles malaria vectors in Colombia. We provide basic information on the geographical distribution, altitudinal range, immature habitats, adult behaviour, feeding preferences and anthropophily, endophily and infectivity rates. We additionally review information on the life cycle, longevity and population fluctuation of Colombian Anopheles species. Emphasis was placed on the primary vectors that have been epidemiologically incriminated in malaria transmission: Anopheles darlingi, Anopheles albimanus and Anopheles nuneztovari. The role of a selection of local, regional or secondary vectors (e.g., Anopheles pseudopunctipennis and Anopheles neivai) is also discussed. We highlight the importance of combining biological, morphological and molecular data for the correct taxonomical determination of a given species, particularly for members of the species complexes. We likewise emphasise the importance of studying the bionomics of primary and secondary vectors along with an examination of the local conditions affecting the transmission of malaria. The presence and spread of the major vectors and the emergence of secondary species capable of transmitting human Plasmodia are of great interest. When selecting control measures, the anopheline diversity in the region must be considered. Variation in macroclimate conditions over a species' geographical range must be well understood and targeted to plan effective control measures based on the population dynamics of the local Anopheles species.

Population structure of the malaria vector Anopheles sinensis (Diptera: Culicidae) in China: two gene pools inferred by microsatellites

BACKGROUND

Anopheles sinensis is a competent malaria vector in China. An understanding of vector population structure is important to the vector-based malaria control programs. However, there is no adequate data of A. sinensis population genetics available yet.

METHODOLOGY/PRINCIPAL FINDINGS

This study used 5 microsatellite loci to estimate population genetic diversity, genetic differentiation and demographic history of A. sinensis from 14 representative localities in China. All 5 microsatellite loci were highly polymorphic across populations, with high allelic richness and heterozygosity. Hardy-Weinberg disequilibrium was found in 12 populations associated with heterozygote deficits, which was likely caused by the presence of null allele and the Wahlund effect. Bayesian clustering analysis revealed two gene pools, grouping samples into two population clusters; one includes six and the other includes eight populations. Out of 14 samples, six samples were mixed with individuals from both gene pools, indicating the coexistence of two genetic units in the areas sampled. The overall differentiation between two genetic pools was moderate (F(ST) = 0.156). Pairwise differentiation between populations were lower within clusters (F(ST) = 0.008-0.028 in cluster I and F(ST) = 0.004-0.048 in cluster II) than between clusters (F(ST) = 0.120-0.201). A reduced gene flow (Nm = 1-1.7) was detected between clusters. No evidence of isolation by distance was detected among populations neither within nor between the two clusters. There are differences in effective population size (Ne = 14.3-infinite) across sampled populations.

CONCLUSIONS/SIGNIFICANCE

Two genetic pools with moderate genetic differentiation were identified in the A. sinensis populations in China. The population divergence was not correlated with geographic distance or barrier in the range. Variable effective population size and other demographic effects of historical population perturbations could be the factors affecting the population differentiation. The structured populations may limit the migration of genes under pressures/selections, such as insecticides and immune genes against malaria.