Anopheles mosquitoes are the principal vectors for malaria and lymphatic filariasis, and evidence for arboviral transmission under laboratory and natural contexts has been demonstrated. Vector management approaches require an understanding of the ecological, epidemiological, and biological contexts of the species in question, and increased interest in gene drive systems for vector control applications has resulted in an increased need for genome assemblies from understudied mosquito vector species. In this study, we present novel genome assemblies for Anopheles crucians, Anopheles freeborni, Anopheles albimanus, and Anopheles quadrimaculatus and examine the evolutionary relationship between these species. We identified 790 shared single-copy orthologs between the newly sequenced genomes and created a phylogeny using 673 of the orthologs, identifying 289 orthologs with evidence for positive selection on at least 1 branch of the phylogeny. Gene ontology terms such as calcium ion signaling, histone binding, and protein acetylation identified as being biased in the set of selected genes. These novel genome sequences will be useful in developing our understanding of the diverse biological traits that drive vectorial capacity in anophelines.

Bacterial content of mosquitoes has given rise to the development of innovative tools that influence and seek to control malaria transmission. This study identified the bacterial microbiota in field-collected female adults of the Anopheles hyrcanus group and three Anopheles species, Anopheles nivipes, Anopheles philippinensis, and Anopheles vagus, from an endemic area in the southeastern part of Ubon Ratchathani Province, northeastern Thailand, near the Lao PDR-Cambodia-Thailand border. A total of 17 DNA libraries were generated from pooled female Anopheles abdomen samples (10 abdomens/ sample). The mosquito microbiota was characterized through the analysis of DNA sequences from the V3-V4 regions of the 16S rRNA gene, and data were analyzed in QIIME2. A total of 3,442 bacterial ASVs were obtained, revealing differences in the microbiota both within the same species/group and between different species/group. Statistical difference in alpha diversity was observed between An. hyrcanus group and An. vagus and between An. nivipes and An. vagus, and beta diversity analyses showed that the bacterial community of An. vagus was the most dissimilar from other species. The most abundant bacteria belonged to the Proteobacteria phylum (48%-75%) in which Pseudomonas, Serratia, and Pantoea were predominant genera among four Anopheles species/group. However, the most significantly abundant genus observed in each Anopheles species/group was as follows: Staphylococcus in the An. hyrcanus group, Pantoea in the An. nivipes, Rosenbergiella in An. philippinensis, and Pseudomonas in An. vagus. Particularly, Pseudomonas sp. was highly abundant in all Anopheles species except An. nivipes. The present study provides the first study on the microbiota of four potential malaria vectors as a starting step towards understanding the role of the microbiota on mosquito biology and ultimately the development of potential tools for malaria control.

Airborne spatial repellency (SR) is characterized and distinguished from other chemical actions including contact locomotor excitation and toxicity. The use of volatile spatial repellents is a potential new intervention class for combatting mosquito-borne pathogen transmission; therefore, continuing investigations on the actions of these chemicals that modify mosquito host-seeking behavior (i.e., bite prevention) is needed. The objective of this study is to characterize the key behavioral avoidance actions of transfluthrin (TFT) to advance spatial repellent development into practical products. Behavioral avoidance responses were observed for adult laboratory strains of Aedes aegypti, Anopheles minimus and An. dirus, and two field populations of An. harrisoni and Ae. aegypti, respectively. Established TFT sublethal (LC50 and LC75), lethal concentrations (LC99) and discriminating concentrations (DCs) were selected corresponding to each mosquito test species. Spatial repellency and contact excitation ('irritancy') responses on adult mosquitoes to TFT were assessed using an excito-repellency assay system. At LC50, TFT exhibited strong avoidance with An. minimus (60.1% escape) and An. dirus (80% escape) laboratory strains, showing between 12 and 16x greater escape response than Ae. aegypti (5% escape). Repellency responses for field collected Ae. aegypti and An. harrisoni were 54.9 and 47.1% escape, respectively. After adjusting the initial contact escape response (a measure of combined irritancy and repellency) to estimate only escape due to contact, the LC50 and LC99 showed moderate escape irritancy with laboratory Ae. aegypti (41.4% escape) and no contact activity against the field population. Adjustment showed only weak contact activity (16.1% escape) in laboratory An. minimus at LC50. Spatial repellency is the predominant mode of action of TFT among colonized and field mosquitoes used in this study. Established baseline (susceptible) dose-response curves assist in optimizing SR products for mosquito control and pathogen transmission prevention.

Anopheles gambiae s.s., Anopheles arabiensis, and Anopheles funestus s.s. are the most important species for malaria transmission. Pyrethroid resistance of these vector mosquitoes is one of the main obstacles against effective vector control. The objective of the present study was to monitor the pyrethroid susceptibility in the 3 major malaria vectors in a highly malaria endemic area in western Kenya and to elucidate the mechanisms of pyrethroid resistance in these species. Gembe East and West, Mbita Division, and 4 main western islands in the Suba district of the Nyanza province in western Kenya were used as the study area. Larval and adult collection and bioassay were conducted, as well as the detection of point mutation in the voltage-gated sodium channel (1014L) by using direct DNA sequencing. A high level of pyrethroid resistance caused by the high frequency of point mutations (L1014S) was detected in An. gambiae s.s. In contrast, P450-related pyrethroid resistance seemed to be widespread in both An. arabiensis and An. funestus s.s. Not a single L1014S mutation was detected in these 2 species. A lack of cross-resistance between DDT and permethrin was also found in An. arabiensis and An. funestus s.s., while An. gambiae s.s. was resistant to both insecticides. It is noteworthy that the above species in the same area are found to be resistant to pyrethroids by their unique resistance mechanisms. Furthermore, it is interesting that 2 different resistance mechanisms have developed in the 2 sibling species in the same area individually. The cross resistance between permethrin and DDT in An. gambiae s.s. may be attributed to the high frequency of kdr mutation, which might be selected by the frequent exposure to ITNs. Similarly, the metabolic pyrethroid resistance in An. arabiensis and An. funestus s.s. is thought to develop without strong selection by DDT.

Despite the medical importance of mosquitoes of the genus Anopheles in the transmission of malaria and other human diseases, its phylogenetic relationships are not settled, and the characteristics of mitochondrial genome (mitogenome) are not thoroughly understood.

Anopheles sundaicus s.l. is an important malaria vector primarily found in coastal landscapes of western and central Indonesia. The species complex has a wide geographical distribution in South and Southeast Asia and exhibits ecological and behavioural variability over its range. Studies on understanding the distribution of different members in the complex and their bionomics related to malaria transmission might be important guiding more effective vector intervention strategies. Female An. sundaicus s.l. were collected from seven provinces, 12 locations in Indonesia representing Sumatra: North Sumatra, Bangka-Belitung, South Lampung, and Bengkulu; in Java: West Java; and the Lesser Sunda Islands: West Nusa Tenggara and East Nusa Tenggara provinces. Sequencing of ribosomal DNA ITS2 gene fragments and two mitochondrial DNA gene markers, COI and cytb, enabled molecular identification of morphologically indistinguishable members of the complex. Findings allowed inference on the distribution of the An. sundaicus s.l. present in Indonesia and further illustrate the phylogenetic relationships of An. epiroticus within the complex. A total of 370 An. sundaicus s.l specimens were analysed for the ITS2 fragment. The ITS2 sequence alignment revealed two consistent species-specific point mutations, a T>C transition at base 479 and a G>T transversion at base 538 that differentiated five haplotypes: TG, CG, TT, CT, and TY. The TG haplotype matched published An. epiroticus-indicative sequences from Thailand, Vietnam and peninsular Malaysia. The previously described insertion event (base 603) was observed in all identified specimens. Analysis of the COI and cytb genes revealed no consistent nucleotide variations that could definitively distinguish An. epiroticus from other members in the Sundaicus Complex. The findings indicate and support the existence of An. epiroticus in North Sumatra and Bangka-Belitung archipelago. Further studies are recommended to determine the full distributional extent of the Sundaicus complex in Indonesia and investigate the role of these species in malaria transmission.

Vector control is a crucial element of anti-malaria campaigns and works best when there is a thorough knowledge of the biology and behaviour of the Anopheles vector species responsible for transmitting malaria within a given locale. With the push to eradicate malaria stronger than ever, there is a growing need to develop and deploy control strategies that exploit the behavioural attributes of local vector species. This is especially true in regions where the vectors are exophagic (i.e., prefer to bite outdoors), exophilic (i.e., prefer to remain outdoors), and zoophagic (i.e., as likely to feed on non-humans as humans). One promising strategy targeting vectors with these behavioural traits is the administration of avermectin-based endectocides, such as ivermectin, to humans and livestock. When ingested in a blood meal, ivermectin has been shown to reduce mosquito survivorship and fecundity in a number of Anopheles species. In this study, the relative toxicity of ivermectin was compared between two zoophagic, exophilic malaria vectors-Anopheles albimanus and Anopheles stephensi.

The Anopheles gambiae sensu lato species complex consists of a number of cryptic species with different habitats and behaviours. These morphologically indistinct species are identified by chromosome banding. Several molecular diagnostic techniques for distinguishing between An. coluzzii and An. gambiae are still under improvement. Although, the current SINE method for identification between An. coluzzii and An. gambiae works reliably, this study describes a refinement of the SINE method to increase sensitivity for identification of An. coluzzii, An. gambiae and An. arabiensis based on amplicon dissociation curve characteristics. Field-collected samples, laboratory-reared colonies and crossed specimens of the two species were used for the design of the protocol. An. gambiae, An. coluzzii, and hybrids of the two species were sampled from Ghana and An. arabiensis from Kenya. Samples were first characterised using conventional SINE PCR method, and further assayed using SYBR green, an intercalating fluorescent dye. The three species and hybrids were clearly differentiated using the melting temperature of the dissociation curves, with derivative peaks at 72°C for An. arabiensis, 75°C for An. gambiae and 86°C for An. coluzzii. The hybrids (An. gambiae / An. coluzzii) showed both peaks. This work is the first to describe a SYBR green real time PCR method for the characterization of An. arabiensis, An. gambiae and An. coluzzii and was purposely designed for basic melt-curve analysis (rather than high-resolution melt-curve) to allow it to be used on a wide range of real-time PCR machines.

Anopheles sinensis is a dominant natural vector of Plasmodium vivax in China, Taiwan, Japan, and Korea. Recent genome sequencing of An. sinensis provides important insights into the genomic basis of vectorial capacity. However, the lack of a physical genome map with chromosome assignment and orientation of sequencing scaffolds hinders comparative analyses with other genomes to infer evolutionary changes relevant to the vector capacity.

Anopheles lesteri and Anopheles sinensis are two major malaria vectors in China and Southeast Asia. They are dramatically different in terms of geographical distribution, host preference, resting habitats, and other traits associated with ecological adaptation and malaria transmission. Both species belong to the Anopheles hyrcanus group, but the extent of genetic differences between them is not well understood. To provide an effective way to differentiate between species and to find useful markers for population genetics studies, we performed a comparative cytogenetic analysis of these two malaria vectors.

Accurate Anopheles species identification is key for effective malaria vector control. Identification primarily depends on morphological analysis of field samples as well as molecular species-specific identifications. During an intra-laboratory assessment (proficiency testing) of the Anopheles funestus group multiplex PCR assay, it was noted that Anopheles arabiensis can be misidentified as Anopheles leesoni, a zoophilic member of the An. funestus group. The aim of this project was, therefore, to ascertain whether other members of the Anopheles gambiae complex can also be misidentified as An. leesoni when using the standard An. funestus multiplex PCR.

Anopheles subpictus and Anopheles sundaicus are closely related species, each comprising several sibling species. Ambiguities exist in the classification of these two nominal species and the specific status of members of these species complexes. Identifying fixed molecular forms and mapping their spatial distribution will help in resolving the taxonomic ambiguities and understanding their relative epidemiological significance.

The so far known species of the Anopheles Claviger Complex, Anopheles claviger s.s. and Anopheles petragnani, can only be distinguished by partial overlapping characteristics of immature stages and by nucleotide sequence variation of the genomic ribosomal DNA (rDNA) internal transcribed spacer 2 (ITS2) region. The known distribution of An. petragnani is so far restricted to the western Mediterranean region, whereas An. claviger s.s. occurs across most of Europe, up to the Middle East and North Africa. In our study, we investigated the larval mosquito fauna in rock pools of the Murg valley (Black Forest, Germany) once a month from April to December 2015.Among other species, larvae belonging to the Anopheles Claviger Complex were found. The fourth instar larvae were morphologically identified by chaetotaxy of the head and abdomen. The results were confirmed by a multiplex PCR and additional sequencing of the amplificates.Of the 1289 collected larvae from the rock pools, seven belonged to the Anopheles Claviger Complex. Five individuals were determined morphologically as An. petragnani and two as An. claviger s.s. The associated mosquito fauna comprised of Aedes japonicus japonicus (548 individuals), Culex pipiens s.l. and Culex torrentium (493 individuals) and Culex hortensis (241 individuals).This is the first record of An. petragnani north of the Alps. Further studies will reveal whether this is an isolated population of An. petragnani and if the investigated rock pool breeding sites represent typical habitats of this species in temperate regions in Central Europe.

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).

Malaria is an important mosquito-borne disease, transmitted to humans by Anopheles mosquitoes. The aim of this study was to gather all records of three main malaria vectors in Iran during the last decades, and to predict the current distribution and the environmental suitability for these species across the country.

Anopheles coluzzii and Anopheles arabiensis belong to the Anopheles gambiae complex and are among the major malaria vectors in sub-Saharan Africa. However, chromosome-level reference genome assemblies are still lacking for these medically important mosquito species.

The age of mosquitoes is a crucial determinant of their ability to transmit pathogens and their resistance to insecticides. We investigated changes to the abundance of proteins found in heads and thoraces of the malaria mosquitoes Anopheles gambiae and Anopheles stephensi as they aged. Protein expression changes were assessed using two-dimensional difference gel electrophoresis and the identity of differentially expressed proteins was determined by using either matrix-assisted laser desorption ionization tandem time-of-flight mass spectrometry or capillary high-pressure liquid chromatography coupled with a linear ion-trap (LTQ)-Orbitrap XL hybrid mass spectrometer. Protein biomarkers were validated by semi quantitative Western blot analysis. Nineteen and nine age dependent protein spots were identified for A. stephensi and A. gambiae, respectively. Among the proteins down-regulated with age were homologs of ADF/Cofilin, cytochome c1, heat shock protein-70 and eukaryotic translation initiation factor 5A (eIF5a). Proteins up-regulated with age included probable methylmalonate-semialdehyde dehydrogenase, voltage-dependent anion-selective channel and fructose bisphosphate aldolase. Semi quantitative Western blot analysis confirmed expression patterns observed by 2-D DIGE for eIF5a and ADF/Cofilin. Further work is recommended to determine whether these biomarkers are robust to infection, blood feeding and insecticide resistance. Robust biomarkers could then be incorporated into rapid diagnostic assays for ecological and epidemiological studies.

To evaluate whether environmental heterogeneity contributes to the genetic heterogeneity in Anopheles triannulatus, larval habitat characteristics across the Brazilian states of Roraima and Pará and genetic sequences were examined. A comparison with Anopheles goeldii was utilised to determine whether high genetic diversity was unique to An. triannulatus. Student t test and analysis of variance found no differences in habitat characteristics between the species. Analysis of population structure of An. triannulatus and An. goeldii revealed distinct demographic histories in a largely overlapping geographic range. Cytochrome oxidase I sequence parsimony networks found geographic clustering for both species; however nuclear marker networks depicted An. triannulatus with a more complex history of fragmentation, secondary contact and recent divergence. Evidence of Pleistocene expansions suggests both species are more likely to be genetically structured by geographic and ecological barriers than demography. We hypothesise that niche partitioning is a driving force for diversity, particularly in An. triannulatus.

Using high-depth whole genome sequencing of F0 mating pairs and multiple individual F1 offspring, we estimated the nuclear mutation rate per generation in the malaria vectors Anopheles coluzzii and Anopheles stephensi by detecting de novo genetic mutations. A purpose-built computer program was employed to filter actual mutations from a deep background of superficially similar artifacts resulting from read misalignment. Performance of filtering parameters was determined using software-simulated mutations, and the resulting estimate of false negative rate was used to correct final mutation rate estimates. Spontaneous mutation rates by base substitution were estimated at 1.00 × 10-9 (95% confidence interval, 2.06 × 10-10-2.91 × 10-9) and 1.36 × 10-9 (95% confidence interval, 4.42 × 10-10-3.18 × 10-9) per site per generation in A. coluzzii and A. stephensi respectively. Although similar studies have been performed on other insect species including dipterans, this is the first study to empirically measure mutation rates in the important genus Anopheles, and thus provides an estimate of µ that will be of utility for comparative evolutionary genomics, as well as for population genetic analysis of malaria vector mosquito species.

Extrinsic environmental factors influence the spatiotemporal dynamics of many organisms, including insects that transmit the pathogens responsible for vector-borne diseases (VBDs). Temperature is an especially important constraint on the fitness of a wide variety of ectothermic insects. A mechanistic understanding of how temperature impacts traits of ectotherms, and thus the distribution of ectotherms and vector-borne infections, is key to predicting the consequences of climate change on transmission of VBDs like malaria. However, the response of transmission to temperature and other drivers is complex, as thermal traits of ectotherms are typically nonlinear, and they interact to determine transmission constraints. In this study, we assess and compare the effect of temperature on the transmission of two malaria parasites, Plasmodium falciparum and Plasmodium vivax, by two malaria vector species, Anopheles gambiae and Anopheles stephensi. We model the nonlinear responses of temperature dependent mosquito and parasite traits (mosquito development rate, bite rate, fecundity, proportion of eggs surviving to adulthood, vector competence, mortality rate, and parasite development rate) and incorporate these traits into a suitability metric based on a model for the basic reproductive number across temperatures. Our model predicts that the optimum temperature for transmission suitability is similar for the four mosquito-parasite combinations assessed in this study, but may differ at the thermal limits. More specifically, we found significant differences in the upper thermal limit between parasites spread by the same mosquito (A. stephensi) and between mosquitoes carrying P. falciparum. In contrast, at the lower thermal limit the significant differences were primarily between the mosquito species that both carried the same pathogen (e.g., A. stephensi and A. gambiae both with P. falciparum). Using prevalence data, we show that the transmission suitability metric ST$$ S(T) $$ calculated from our mechanistic model is consistent with observed P. falciparum prevalence in Africa and Asia but is equivocal for P. vivax prevalence in Asia, and inconsistent with P. vivax prevalence in Africa. We mapped risk to illustrate the number of months various areas in Africa and Asia predicted to be suitable for malaria transmission based on this suitability metric. This mapping provides spatially explicit predictions for suitability and transmission risk.