Annotated Differentially Expressed Salivary Proteins of Susceptible and Insecticide-Resistant Mosquitoes of Anopheles stephensi

Vector control is one of the major global strategies for control of malaria. However, the major obstacle for vector control is the development of multiple resistances to organochlorine, organophosphorus insecticides and pyrethroids that are currently being used in public health for spraying and in bednets. Salivary glands of vectors are the first target organ for human-vector contact during biting and parasite-vector contact prior to parasite development in the mosquito midguts. The salivary glands secrete anti-haemostatic, anti-inflammatory biologically active molecules to facilitate blood feeding from the host and also inadvertently inject malaria parasites into the vertebrate host. The Anopheles stephensi mosquito, an urban vector of malaria to both human and rodent species has been identified as a reference laboratory model to study mosquito—parasite interactions. In this study, we adopted a conventional proteomic approach of 2D-electrophoresis coupled with MALDI-TOF mass spectrometry and bioinformatics to identify putative differentially expressed annotated functional salivary proteins between An. stephensi susceptible and multiresistant strains with same genetic background. Our results show 2D gel profile and MALDI-TOF comparisons that identified 31 differentially expressed putative modulated proteins in deltamethrin/DDT resistant strains of An. stephensi. Among these 15 proteins were found to be upregulated and 16 proteins were downregulated. Our studies interpret that An. stephensi (multiresistant) caused an upregulated expression of proteins and enzymes like cytochrome 450, short chain dehyrdogenase reductase, phosphodiesterase etc that may have an impact in insecticide resistance and xenobiotic detoxification. Our study elucidates a proteomic response of salivary glands differentially regulated proteins in response to insecticide resistance development which include structural, redox and regulatory enzymes of several pathways. These identified proteins may play a role in regulating mosquito biting behavior patterns and may have implications in the development of malaria parasites in resistant mosquitoes during parasite transmission.     

Knockout of Anopheles stephensi immune gene LRIM1 by CRISPR-Cas9 reveals its unexpected role in reproduction and vector competence

PfSPZ Vaccine against malaria is composed of Plasmodium falciparum (Pf) sporozoites (SPZ) manufactured using aseptically reared Anopheles stephensi mosquitoes. Immune response genes of Anopheles mosquitoes such as Leucin-Rich protein (LRIM1), inhibit Plasmodium SPZ development (sporogony) in mosquitoes by supporting melanization and phagocytosis of ookinetes. With the aim of increasing PfSPZ infection intensities, we generated an A. stephensi LRIM1 knockout line, Δaslrim1, by embryonic genome editing using CRISPR-Cas9. Δaslrim1 mosquitoes had a significantly increased midgut bacterial load and an altered microbiome composition, including elimination of commensal acetic acid bacteria. The alterations in the microbiome caused increased mosquito mortality and unexpectedly, significantly reduced sporogony. The survival rate of Δaslrim1 mosquitoes and their ability to support PfSPZ development, were partially restored by antibiotic treatment of the mosquitoes, and fully restored to baseline when Δaslrim1 mosquitoes were produced aseptically. Deletion of LRIM1 also affected reproductive capacity: oviposition, fecundity and male fertility were significantly compromised. Attenuation in fecundity was not associated with the altered microbiome. This work demonstrates that LRIM1’s regulation of the microbiome has a major impact on vector competence and longevity of A. stephensi. Additionally, LRIM1 deletion identified an unexpected role for this gene in fecundity and reduction of sperm transfer by males.     

A mating-induced reproductive gene promotes Anopheles tolerance to Plasmodium falciparum infection

Anopheles mosquitoes have transmitted Plasmodium parasites for millions of years, yet it remains unclear whether they suffer fitness costs to infection. Here we report that the fecundity of virgin and mated females of two important vectors—Anopheles gambiae and Anopheles stephensi—is not affected by infection with Plasmodium falciparum, demonstrating that these human malaria parasites do not inflict this reproductive cost on their natural mosquito hosts. Additionally, parasite development is not impacted by mating status. However, in field studies using different P. falciparum isolates in Anopheles coluzzii, we find that Mating-Induced Stimulator of Oogenesis (MISO), a female reproductive gene strongly induced after mating by the sexual transfer of the steroid hormone 20-hydroxyecdysone (20E), protects females from incurring fecundity costs to infection. MISO-silenced females produce fewer eggs as they become increasingly infected with P. falciparum, while parasite development is not impacted by this gene silencing. Interestingly, previous work had shown that sexual transfer of 20E has specifically evolved in Cellia species of the Anopheles genus, driving the co-adaptation of MISO. Our data therefore suggest that evolution of male-female sexual interactions may have promoted Anopheles tolerance to P. falciparum infection in the Cellia subgenus, which comprises the most important malaria vectors.     

Blood Feeding and Plasmodium Infection Alters the miRNome of Anopheles stephensi

Blood feeding is an integral process required for physiological functions and propagation of the malaria vector Anopheles. During blood feeding, presence of the malaria parasite, Plasmodium in the blood induces several host effector molecules including microRNAs which play important roles in the development and maturation of the parasite within the mosquito. The present study was undertaken to elucidate the dynamic expression of miRNAs during gonotrophic cycle and parasite development in Anopheles stephensi. Using next generation sequencing technology, we identified 126 miRNAs of which 17 were novel miRNAs. The miRNAs were further validated by northern hybridization and cloning. Blood feeding and parasitized blood feeding in the mosquitoes revealed regulation of 13 and 16 miRNAs respectively. Expression profiling of these miRNAs revealed that significant miRNAs were down-regulated upon parasitized blood feeding with a repertoire of miRNAs showing stage specific up-regulation. Expression profiles of significantly modulated miRNAs were further validated by real time PCR. Target prediction of regulated miRNAs revealed overlapping targeting by different miRNAs. These targets included several metabolic pathways including metabolic, redox homeostasis and protein processing machinery components. Our analysis revealed tight regulation of specific miRNAs post blood feeding and parasite infection in An. stephensi. Such regulated expression suggests possible role of these miRNAs during gonotrophic cycle in mosquito. Another set of miRNAs were also significantly regulated at 42 h and 5 days post infection indicating parasite stage-specific role of host miRNAs. This study will result in better understanding of the role of miRNAs during gonotrophic cycle and parasite development in mosquito and can probably facilitate in devising novel malaria control strategies at vector level.     

Stable and heritable gene silencing in the malaria vector Anopheles stephensi

Heritable RNA interference (RNAi), triggered from stably expressed transgenes with an inverted repeat (IR) configuration, is an important tool for reverse genetic studies. Here we report on the development of stable RNAi in Anopheles stephensi mosquitoes, the major vector of human malaria in Asia. Trans genic mosquitoes stably expressing a RNAi transgene, designed to produce intron-spliced double-stranded RNA (dsRNA) targeting the green fluorescent protein EGFP gene, were crossed to an EGFP-expressing target line. EGFP expression was dramatically reduced at both the protein and RNA levels. The levels of gene silencing depended upon the RNAi gene copy number and its site of integration. These results demonstrate that specific RNAi-mediated knockdown of gene function can be achieved with high efficiency in Anopheles. This will be invaluable to systematically unravel the function of Anopheles genes determining the vectorial capacity of the malaria parasite.     

Elizabethkingia anophelis: Molecular Manipulation and Interactions with Mosquito Hosts

Flavobacteria (members of the family Flavobacteriaceae) dominate the bacterial community in the Anopheles mosquito midgut. One such commensal, Elizabethkingia anophelis, is closely associated with Anopheles mosquitoes through transstadial persistence (i.e., from one life stage to the next); these and other properties favor its development for paratransgenic applications in control of malaria parasite transmission. However, the physiological requirements of E. anophelis have not been investigated, nor has its capacity to perpetuate despite digestion pressure in the gut been quantified. To this end, we first developed techniques for genetic manipulation of E. anophelis, including selectable markers, reporter systems (green fluorescent protein [GFP] and NanoLuc), and transposons that function in E. anophelis. A flavobacterial expression system based on the promoter PompA was integrated into the E. anophelis chromosome and showed strong promoter activity to drive GFP and NanoLuc reporter production. Introduced, GFP-tagged E. anophelis associated with mosquitoes at successive developmental stages and propagated in Anopheles gambiae and Anopheles stephensi but not in Aedes triseriatus mosquitoes. Feeding NanoLuc-tagged cells to A. gambiae and A. stephensi in the larval stage led to infection rates of 71% and 82%, respectively. In contrast, a very low infection rate (3%) was detected in Aedes triseriatus mosquitoes under the same conditions. Of the initial E. anophelis cells provided to larvae, 23%, 71%, and 85% were digested in A. stephensi, A. gambiae, and Aedes triseriatus, respectively, demonstrating that E. anophelis adapted to various mosquito midgut environments differently. Bacterial cell growth increased up to 3-fold when arginine was supplemented in the defined medium. Furthermore, the number of NanoLuc-tagged cells in A. stephensi significantly increased when arginine was added to a sugar diet, showing it to be an important amino acid for E. anophelis. Animal erythrocytes promoted E. anophelis growth in vivo and in vitro, indicating that this bacterium could obtain nutrients by participating in erythrocyte lysis in the mosquito midgut.     

Malaria vectors of Papua New Guinea

Understanding malaria transmission in Papua New Guinea (PNG) requires exact knowledge of which Anopheles species are transmitting malaria and is complicated by the cryptic species status of many of these mosquitoes. To identify the malaria vectors in PNG we studied Anopheles specimens from 232 collection localities around human habitation throughout PNG (using CO₂ baited light traps and human bait collections). A total of 22,970 mosquitoes were individually assessed using a Plasmodium sporozoite enzyme-linked immunosorbent assay to identify Plasmodiumfalciparum, Plasmodiumvivax and Plasmodiummalariae circumsporozoite proteins. All mosquitoes were identified to species by morphology and/or PCR. Based on distribution, abundance and their ability to develop sporozoites, we identified five species as major vectors of malaria in PNG. These included: Anophelesfarauti, Anopheleshinesorum (incriminated here, to our knowledge, for the first time), Anophelesfarauti 4, Anopheleskoliensis and Anophelespunctulatus. Anopheleslongirostris and Anophelesbancroftii were also incriminated in this study. Surprisingly, An. longirostris showed a high incidence of infections in some areas. A newly identified taxon within the Punctulatus Group, tentatively called An. farauti 8, was also found positive for circumsporozoite protein. These latter three species, together with Anopheleskarwari and Anophelessubpictus, incriminated in other studies, appear to be only minor vectors, while Anophelesfarauti 6 appears to be the major vector in the highland river valleys (>1500 m above sea level). The nine remaining Anopheles species found in PNG have been little studied and their bionomics are unknown; most appear to be uncommon with limited distribution and their possible role in malaria transmission has yet to be determined.     

Anopheline species and their Plasmodium infection status in Aligarh,India

Malaria is a global issue and India contributes substantially to global malaria incidence. Information related to malaria vectors is very limited in Aligarh. The environmental and climatological situations permit the continual breeding of vectors in permanent breeding sites. This study was designed with the aim to screen all the anophelines species and possible malaria vectors in three different localities of Aligarh. Anopheles mosquitoes were collected from three different localities (Fort, Jalali and Tappal) during peak malaria transmission season (July to November) by using mouth aspirator and CDC light traps. Enzyme-linked immunosorbent assay (ELISA) was done to detect Plasmodium falciparum, Plasmodium vivax-210 and P. vivax-247 circumsporozoite proteins (CSP) from the collected female species. A total of 794 female anopheline mosquitoes belonging to 7 species were collected by different methods. Circumsporozoite protein–enzyme-linked immunosorbent assay was performed with 780 anopheline mosquitoes out of which 13 mosquitoes were positive in CSP–ELISA. Thus, the overall infection rate was 1.66% (13/780). Four (0.51%) mosquitoes belonging to three species were positive for P. falciparum, 7 (0.89%) mosquitoes belonging to three species were positive for VK 210 and 2 (0.25%) mosquitoes belonging to Anopheles culicifacies and Anopheles stephensi species were positive for VK 247. No mixed infection was found in this study. According to species, the highest infection rate was observed in An. culicifacies (7/288, 2.43%) followed by An. stephensi (2.40%) and Anopheles annularis (1.98%). An. culicifacies and An. stephensi were previously incriminated as malaria vectors in Aligarh. There was, however, no previous report in favor of infections in An. annularis in Aligarh. The on-going Malaria Control Program in India needs up to date information on malaria vectors. A major challenge is the lack of knowledge about vectors and their role in malaria transmission. Findings of this study suggested that in the absence of major malaria vectors there is a possibility that other Anopheles species may have been playing a role in malaria transmission in Aligarh.     

Evaluation of Insecticides Susceptibility and Malaria Vector Potential of Anopheles annularis s.l. and Anopheles vagus in Assam,India

During the recent past, development of DDT resistance and reduction to pyrethroid susceptibility among the malaria vectors has posed a serious challenge in many Southeast Asian countries including India. Current study presents the insecticide susceptibility and knock-down data of field collected Anopheles annularis sensu lato and An. vagus mosquito species from endemic areas of Assam in northeast India. Anopheles annularis s.l. and An. vagus adult females were collected from four randomly selected sentinel sites in Orang primary health centre (OPHC) and Balipara primary health centre (BPHC) areas, and used for testing susceptibility to DDT, malathion, deltamethrin and lambda-cyhalothrin. After insecticide susceptibility tests, mosquitoes were subjected to VectorTest^™ assay kits to detect the presence of malaria sporozoite in the mosquitoes. An. annularis s.l. was completely susceptible to deltamethrin, lambda-cyhalothrin and malathion in both the study areas. An. vagus was highly susceptible to deltamethrin in both the areas, but exhibited reduced susceptibility to lambda-cyhalothrin in BPHC. Both the species were resistant to DDT and showed very high KDT₅₀ and KDT₉₉ values for DDT. Probit model used to calculate the KDT₅₀ and KDT₉₉ values did not display normal distribution of percent knock-down with time for malathion in both the mosquito species in OPHC (p<0.05) and An. vagus in BPHC (χ² = 25.3; p = 0.0), and also for deltamethrin to An. vagus in BPHC area (χ² = 15.4; p = 0.004). Minimum infection rate (MIR) of Plasmodium sporozoite for An. vagus was 0.56 in OPHC and 0.13 in BPHC, while for An. annularis MIR was found to be 0.22 in OPHC. Resistance management strategies should be identified to delay the expansion of resistance. Testing of field caught Anopheles vectors from different endemic areas for the presence of malaria sporozoite may be useful to ensure their role in malaria transmission.     

Behavioral Cost & Overdominance in Anopheles gambiae

In response to the widespread use of control strategies such as Insecticide Treated Nets (ITN), Anopheles mosquitoes have evolved various resistance mechanisms. Kdr is a mutation that provides physiological resistance to the pyrethroid insecticides family (PYR). In the present study, we investigated the effect of the Kdr mutation on the ability of female An. gambiae to locate and penetrate a 1cm-diameter hole in a piece of netting, either treated with insecticide or untreated, to reach a bait in a wind tunnel. Kdr homozygous, PYR-resistant mosquitoes were the least efficient at penetrating an untreated damaged net, with about 51% [39-63] success rate compared to 80% [70-90] and 78% [65-91] for homozygous susceptible and heterozygous respectively. This reduced efficiency, likely due to reduced host-seeking activity, as revealed by mosquito video-tracking, is evidence of a recessive behavioral cost of the mutation. Kdr heterozygous mosquitoes were the most efficient at penetrating nets treated with PYR insecticide, thus providing evidence for overdominance, the rarely-described case of heterozygote advantage conveyed by a single locus. The study also highlights the remarkable capacity of female mosquitoes, whether PYR-resistant or not, to locate holes in bed-nets.     

The Genetic Basis of Host Preference and Resting Behavior in the Major African Malaria Vector,Anopheles arabiensis

Malaria transmission is dependent on the propensity of Anopheles mosquitoes to bite humans (anthropophily) instead of other dead end hosts. Recent increases in the usage of Long Lasting Insecticide Treated Nets (LLINs) in Africa have been associated with reductions in highly anthropophilic and endophilic vectors such as Anopheles gambiae s.s., leaving species with a broader host range, such as Anopheles arabiensis, as the most prominent remaining source of transmission in many settings. An. arabiensis appears to be more of a generalist in terms of its host choice and resting behavior, which may be due to phenotypic plasticity and/or segregating allelic variation. To investigate the genetic basis of host choice and resting behavior in An. arabiensis we sequenced the genomes of 23 human-fed and 25 cattle-fed mosquitoes collected both in-doors and out-doors in the Kilombero Valley, Tanzania. We identified a total of 4,820,851 SNPs, which were used to conduct the first genome-wide estimates of “SNP heritability” for host choice and resting behavior in this species. A genetic component was detected for host choice (human vs cow fed; permuted P = 0.002), but there was no evidence of a genetic component for resting behavior (indoors versus outside; permuted P = 0.465). A principal component analysis (PCA) segregated individuals based on genomic variation into three groups which were characterized by differences at the 2Rb and/or 3Ra paracentromeric chromosome inversions. There was a non-random distribution of cattle-fed mosquitoes between the PCA clusters, suggesting that alleles linked to the 2Rb and/or 3Ra inversions may influence host choice. Using a novel inversion genotyping assay, we detected a significant enrichment of the standard arrangement (non-inverted) of 3Ra among cattle-fed mosquitoes (N = 129) versus all non-cattle-fed individuals (N = 234; χ², p = 0.007). Thus, tracking the frequency of the 3Ra in An. arabiensis populations may be of use to infer selection on host choice behavior within these vector populations; possibly in response to vector control. Controlled host-choice assays are needed to discern whether the observed genetic component has a direct relationship with innate host preference. A better understanding of the genetic basis for host feeding behavior in An. arabiensis may also open avenues for novel vector control strategies based on driving genes for zoophily into wild mosquito populations.     

Starvation at the larval stage increases the vector competence of Aedes aegypti females for Zika virus

Aedes aegypti is the primary vector of Zika virus (ZIKV), a flavivirus which typically presents itself as febrile-like symptoms in humans but can also cause neurological and pregnancy complications. The transmission cycle of mosquito-borne arboviruses such as ZIKV requires that various key tissues in the female mosquito get productively infected with the virus before the mosquito can transmit the virus to another vertebrate host. Following ingestion of a viremic blood-meal from a vertebrate, ZIKV initially infects the midgut epithelium before exiting the midgut after blood-meal digestion to disseminate to secondary tissues including the salivary glands. Here we investigated whether smaller Ae. aegypti females resulting from food deprivation as larvae exhibited an altered vector competence for blood-meal acquired ZIKV relative to larger mosquitoes. Midguts from small ‘Starve’ and large ‘Control’ Ae. aegypti were dissected to visualize by transmission electron microscopy (TEM) the midgut basal lamina (BL) as physical evidence for the midgut escape barrier showing Starve mosquitoes with a significantly thinner midgut BL than Control mosquitoes at two timepoints. ZIKV replication was inhibited in Starve mosquitoes following intrathoracic injection of virus, however, Starve mosquitoes exhibited a significantly higher midgut escape and population dissemination rate at 9 days post-infection (dpi) via blood-meal, with more virus present in saliva and head tissue than Control by 10 dpi and 14 dpi, respectively. These results indicate that Ae. aegypti developing under stressful conditions potentially exhibit higher midgut infection and dissemination rates for ZIKV as adults, Thus, variation in food intake as larvae is potentially a source for variable vector competence levels of the emerged adults for the virus.     

Anopheles stephensi in Africa: vector control opportunities for cobreeding An. stephensi and Aedes arbovirus vectors

Anopheles stephensi is an urban malaria vector native in some Asian countries and recently emerged in Africa as an invasive vector competent in transmitting Plasmodium falciparum and P. vivax. The coexistence of An. stephensi and Aedes arboviral vectors offers an optimal opportunity for successful integrated vector management with limited resources.     

Detection of Mayaro virus in Aedes aegypti mosquitoes circulating in Goiânia-Goiás-Brazil

The Arboviral diseases are caused by arthropod-borne viruses, such as Mayaro virus (MAYV), the etiological agent of Mayaro fever. This disease has been drawing the attention of the public health authorities for the increased number of cases likely due to virus adaptation for survival to urban areas as well as infection and multiplication in other vectors insects. Therefore, this work aimed to identify the MAYV infecting Aedes aegypti mosquitoes in Goiânia, the capital of state of Goiás, Brazil. For the development of study, the larvae of A. aegypti were collected in Basic Health Units from different regions of Goiânia then the larvae were grown to adult mosquitoes in controlled laboratory conditions. The female mosquitoes were submitted to the procedure of head and body separation. The RNAs obtained from these samples were analyzed by real-time PCR for identification of arboviruses. We only detect the presence of MAVY in the mosquitoes, in this sense our findings suggest that A. aegypti harbor MAYV in different anatomical sites, and potentially the process of vertical transmission of MAYV can occur in this vector.     

CRISPR/Cas9 modified An. gambiae carrying kdr mutation L1014F functionally validate its contribution in insecticide resistance and combined effect with metabolic enzymes

Insecticide resistance in Anopheles mosquitoes is a major obstacle in maintaining the momentum in reducing the malaria burden; mitigating strategies require improved understanding of the underlying mechanisms. Mutations in the target site of insecticides (the voltage gated sodium channel for the most widely used pyrethroid class) and over-expression of detoxification enzymes are commonly reported, but their relative contribution to phenotypic resistance remain poorly understood. Here we present a genome editing pipeline to introduce single nucleotide polymorphisms in An. gambiae which we have used to study the effect of the classical kdr mutation L1014F (L995F based on An. gambiae numbering), one of the most widely distributed resistance alleles. Introduction of 1014F in an otherwise fully susceptible genetic background increased levels of resistance to all tested pyrethroids and DDT ranging from 9.9-fold for permethrin to >24-fold for DDT. The introduction of the 1014F allele was sufficient to reduce mortality of mosquitoes after exposure to deltamethrin treated bednets, even as the only resistance mechanism present. When 1014F was combined with over-expression of glutathione transferase Gste2, resistance to permethrin increased further demonstrating the critical combined effect between target site resistance and detoxification enzymes in vivo. We also show that mosquitoes carrying the 1014F allele in homozygosity showed fitness disadvantages including increased mortality at the larval stage and a reduction in fecundity and adult longevity, which can have consequences for the strength of selection that will apply to this allele in the field.     

The Role of Oxidative Stress in the Longevity and Insecticide Resistance Phenotype of the Major Malaria Vectors Anopheles arabiensis and Anopheles funestus

Oxidative stress plays numerous biological roles, both functional and pathological. The role of oxidative stress in various epidemiologically relevant biological traits in Anopheles mosquitoes is not well established. In this study, the effects of oxidative stress on the longevity and insecticide resistance phenotype in the major malaria vector species An. arabiensis and An. funestus were examined. Responses to dietary copper sulphate and hydrogen peroxide were used as proxies for the oxidative stress phenotype by determining the effect of copper on longevity and hydrogen peroxide lethal dose. Glutathione peroxidase and catalase activities were determined colorimetrically. Oxidative burden was quantified as protein carbonyl content. Changes in insecticide resistance phenotype were monitored by WHO bioassay. Insecticide resistant individuals showed an increased capacity for coping with oxidative stress, mediated by increased glutathione peroxidase and catalase activity. This effect was observed in both species, as well as in laboratory strains and F₁ individuals derived from wild-caught An. funestus mothers. Phenotypic capacity for coping with oxidative stress was greatest in strains with elevated Cytochrome P450 activity. Synergism of oxidative stress defence enzymes by dietary supplementation with haematin, 3-Amino-1, 2, 4-triazole and Sodium diethyldithiocarbamate significantly increased pyrethroid-induced mortality in An. arabiensis and An. funestus. It is therefore concluded that defence against oxidative stress underlies the augmentation of the insecticide resistance phenotype associated with multiple blood-feeding. This is because multiple blood-feeding ultimately leads to a reduction of oxidative stress in insecticide resistant females, and also reduces the oxidative burden induced by DDT and pyrethroids, by inducing increased glutathione peroxidase activity. This study highlights the importance of oxidative stress in the longevity and insecticide resistance phenotype in malaria vectors.     

Molecular Taxonomy Provides New Insights into Anopheles Species of the Neotropical Arribalzagia Series

Phylogenetic analysis of partial mitochondrial cytochrome oxidase c subunit I (COI) and nuclear internal transcribed spacer 2 (ITS2) sequences were used to evaluate initial identification and to investigate phylogenetic relationships of seven Anopheles morphospecies of the Arribalzagia Series from Colombia. Phylogenetic trees recovered highly supported clades for An. punctimaculas.s., An. calderoni, An. malefactor s.l., An. neomaculipalpus, An. apicimacula s.l., An. mattogrossensis and An. peryassui. This study provides the first molecular confirmation of An. malefactorfrom Colombia and discovered conflicting patterns of divergence for the molecular markers among specimens from northeast and northern Colombia suggesting the presence of two previously unrecognized Molecular Operational Taxonomic Units (MOTUs). Furthermore, two highly differentiated An. apicimacula MOTUs previously found in Panama were detected. Overall, the combined molecular dataset facilitated the detection of known and new Colombian evolutionary lineages, and constitutes the baseline for future research on their bionomics, ecology and potential role as malaria vectors.