Complex genome evolution in Anopheles coluzzii associated with increased insecticide usage in Mali

One Health, many approaches: integrated vector management strategies support One Health goals

Integrated pest management (IPM) laid the foundation for integrated vector management (IVM) by introducing strategies that prioritize prevention over intervention and the use of diverse management and control tools for arthropod vectors. Both approaches focus on reducing chemical reliance for management of targeted organisms by employing data-driven decisions and incorporating the use of additional non-chemical based management tactics. While IPM and IVM differ in their specific applications and management contexts, many of their fundamental principles remain the same. By diversifying prevention and control options, these management practices support healthier environments, animals, and populace, the three spheres of the One Health paradigm. One Health principles support and highlight the connections between animal, human, and environmental health and how these spheres influence one another. However, the connections and overlapping goals between IVM and One Health are infrequently discussed in tandem. This review will examine the numerous and varied methods of IVM for reducing global disease burden, detail the benefits of using IVM to achieve One Health goals, as well as benefits and considerations to incorporating a One Health lens with IVM.

Genome wide population genetics and molecular surveillance of insecticide resistance in Anopheles stephensi mosquitoes from Awash Sebat Kilo in Ethiopia

Since the detection of the Asian mosquito Anopheles stephensi in Dijbouti in 2012, it has spread throughout the Horn of Africa. This invasive vector continues to expand across the continent and is a significant threat to malaria control programs. Vector control methods, including insecticide-treated nets and indoor residual spraying, have substantially reduced the malaria burden. However, the increasing prevalence of mosquitoes resistant to insecticides, including An. stephensi populations, undermines ongoing malaria elimination efforts. Understanding population structure, gene flow between populations, and the distribution of insecticide resistance mutations is essential for guiding effective malaria control strategies. Here, we generated whole genome sequencing data for An. stephensi sourced from Awash Sebat Kilo, Ethiopia (n = 27) and compared with South Asian populations (n = 45; India and Pakistan) to assess genomic diversity, population structure, and uncovering insecticide resistance mutations. Population structure analysis using genome-wide single nucleotide polymorphisms (n = 15,533,476) revealed Ethiopian isolates clustering as a distinct ancestral group, separate from South Asian isolates. Three insecticide resistance-associated SNPs (gaba gene: A296S and V327I; vgsc L1014F) were detected. Evidence of ongoing selection was found in several loci, including genes previously associated with neonicotinoids, ivermectin, DDT, and pyrethroid resistance. This study represents the first whole genome population genetics study of invasive An. stephensi, revealing genomic differences from South Asian populations, which can be used for future assessments of vector population dispersal and detection of insecticide resistance mechanisms.

A Comprehensive Review: Biology of Anopheles squamosus, an Understudied Malaria Vector in Africa

With ongoing global efforts to eliminate malaria, several countries have entered a pre-elimination stage, in which populations of the primary mosquito vector are reduced and the rates of malaria transmission are low. Efforts to eliminate malaria have changed the dynamics of disease transmission, causing a shift in anopheline community composition in some regions. One of the most abundant species of Anopheles found in pre-elimination areas, such as southern Africa, is Anopheles squamosus (Theobald, 1901; Diptera: Culicidae). Although An. squamosus has been documented across the African continent, limited information is available about its biology. Knowing what, when, and where the malaria vector bites humans can help develop effective malaria control strategies. The aim of this review is to compile the information available on An. squamosus's biology, ecology, distribution, behavior, and potential role in disease transmission. The review of current literature suggested that An. squamosus is abundant throughout the African continent. Studies indicate that this species occupies the same range and exhibits similar behaviors to primary malaria vectors in southern and eastern African countries, like Zambia. In conclusion, An. squamosus continues to be an understudied species that has circumvented disease control measures and further studies are needed to develop effective control strategies.

A single mutation G454A in the P450 CYP9K1 drives pyrethroid resistance in the major malaria vector Anopheles funestus reducing bed net efficacy

Metabolic mechanisms conferring pyrethroid resistance in malaria vectors are jeopardizing the effectiveness of insecticide-based interventions, and identification of their markers is a key requirement for robust resistance management. Here, using a field-lab-field approach, we demonstrated that a single mutation G454A in the P450 CYP9K1 is driving pyrethroid resistance in the major malaria vector Anopheles funestus in East and Central Africa. Drastic reduction in CYP9K1 diversity was observed in Ugandan samples collected in 2014, with the selection of a predominant haplotype (G454A mutation at 90%), which was completely absent in the other African regions. However, 6 years later (2020) the Ugandan 454A-CYP9K1 haplotype was found predominant in Cameroon (84.6%), but absent in Malawi (Southern Africa) and Ghana (West Africa). Comparative in vitro heterologous expression and metabolism assays revealed that the mutant 454A-CYP9K1 (R) allele significantly metabolizes more type II pyrethroid (deltamethrin) compared with the wild G454-CYP9K1 (S) allele. Transgenic Drosophila melanogaster flies expressing 454A-CYP9K1 (R) allele exhibited significantly higher type I and II pyrethroids resistance compared to flies expressing the wild G454-CYP9K1 (S) allele. Furthermore, laboratory testing and field experimental hut trials in Cameroon demonstrated that mosquitoes harboring the resistant 454A-CYP9K1 allele significantly survived pyrethroids exposure (odds ratio = 567, P < 0.0001). This study highlights the rapid spread of pyrethroid-resistant CYP9K1 allele, under directional selection in East and Central Africa, contributing to reduced bed net efficacy. The newly designed DNA-based assay here will add to the toolbox of resistance monitoring and improving its management strategies.

Insecticide resistance mutations of Anopheles species in the Republic of Korea

The number of reported malaria cases transmitted by Anopheles mosquitoes in the Republic of Korea (ROK) increased from 420 in 2022 to 746 in 2023, a 77.6% increase. Eight Anopheles species are currently reported in the ROK, including six species belonging to the Anopheles Hyrcanus Group and one species each belonging to the Barbirostris Group and Lindesayi Group. However, studies on insecticide resistance in the ROK has predominantly concentrated on Anopheles sinensis or more broadly, members of the Hyrcanus Group. Reported differences in vector competence and ecological characteristics of mosquito species in the ROK highlight the importance for conducting accurate evaluations of insecticide resistance for each of the Anopheles species for informing the potential efficacy of vector control to reduce malaria transmission. All eight species of Anopheles mosquitoes were collected in/near the demilitarized zone (DMZ), a malaria high-risk region in the ROK. Additional specimens were collected in Seoul [Yongsan US Army Garrison (USAG)] and Pyeongtaek (Humphreys USAG) where malaria risks are much lower. Anopheles mosquitoes were identified to species using a multiplex PCR method and then evaluated for the presence of acetylcholinesterase-1 (ace-1) and voltage-gated sodium channel (vgsc) regions to identify mutations linked to insecticide resistance. Analysis of the ace-1 region identified insecticide resistance alleles in four species of the Hyrcanus Group (An. sinensis, An. kleini, An. belenrae, and An. pullus), while ace-1 resistance alleles were not observed in the other four species. The screening of the vgsc gene fragment confirmed the presence of resistant alleles only in An. sinensis (considered a poor malaria vector) and An. kleini (a primary malaria vector) in the ROK. This study represents a preliminary investigation of insecticide resistance mutations across all Anopheles species in the ROK. These findings are crucial in advancing mosquito control strategies to mitigate future malaria infections.

Copy number variants underlie major selective sweeps in insecticide resistance genes in Anopheles arabiensis

To keep ahead of the evolution of resistance to insecticides in mosquitoes, national malaria control programmes must make use of a range of insecticides, both old and new, while monitoring resistance mechanisms. The outdoor-biting malaria vector Anopheles arabiensis is of increasing concern for malaria transmission because it is apparently less susceptible to many indoor control interventions, yet knowledge of its mechanisms of resistance remains limited. Furthermore, comparatively little is known in general about resistance to non-pyrethroid insecticides such as pirimiphos-methyl (PM), which are crucial for effective control in the context of globally high resistance to pyrethroids. We performed a genome-wide association study to determine the molecular mechanisms of resistance to the pyrethroid deltamethrin (commonly used in bednets) and PM (widespread use for indoor spraying), in An. arabiensis from 2 regions in Tanzania. Genomic regions of positive selection in these populations were largely driven by copy number variants (CNVs) in gene families involved in metabolic resistance. We found evidence of a new gene cluster involved in resistance to PM, identifying a strong selective sweep tied to a CNV in the carboxylesterase genes Coeae2g - Coeae6g. Using complementary data from another malaria vector, An. coluzzii, in Ghana, we show that copy number at this locus is significantly associated with PM resistance. Similarly, for deltamethrin, resistance was strongly associated with a novel CNV allele in the Cyp6aa / Cyp6p cluster (Cyp6aap_Dup33). Against this background of metabolic resistance, resistance caused by mutations in the insecticide target sites was very rare or absent. Mutations in the pyrethroid target site Vgsc were at very low frequency in Tanzania, yet combining these samples with 3 An. arabiensis individuals from West Africa revealed a startling evolutionary diversity, with up to 5 independent origins of Vgsc-995 mutations found within just 8 haplotypes. Thus, despite having been first recorded over 10 years ago, Vgsc resistance mutations in Tanzanian An. arabiensis have remained at stable low frequencies. Overall, our results provide a new copy number marker for monitoring resistance to PM in malaria mosquitoes, and reveal the complex picture of resistance patterns in An. arabiensis.

Cryptic population structure and insecticide resistance in Anopheles gambiae from the southern Democratic Republic of Congo

The Democratic Republic of Congo (DRC) suffers from one of the highest malaria burdens worldwide, but information on its Anopheles vector populations is relatively limited. Preventative malaria control in DRC is reliant on pyrethroid-treated nets, raising concerns over the potential impacts of insecticide resistance. We sampled Anopheles gambiae from three geographically distinct populations (Kimpese, Kapolowe and Mikalayi) in southern DRC, collecting from three sub-sites per population and characterising mosquito collections from each for resistance to pyrethroids using WHO tube bioassays. Resistance to each of three different pyrethroids was generally high in An. gambiae with < 92% mortality in all tests, but varied between collections, with mosquitoes from Kimpese being the most resistant. Whole genome sequencing of 165 An. gambiae revealed evidence for genetic differentiation between Kimpese and Kapolowe/Mikalayi, but not between the latter two sample sites despite separation of approximately 800 km. Surprisingly, there was evidence of population structure at a small spatial scale between collection subsites in Kimpese, despite separation of just tens of kilometres. Intra-population (H12) and inter-population (FST) genome scans identified multiple peaks corresponding to genes associated with insecticide resistance such as the voltage gated sodium channel (Vgsc) target site on chromosome 2L, a Cyp6 cytochrome P450 cluster on chromosome arm 2R, and the Cyp9k1 P450 gene on chromosome X. In addition, in the Kimpese subsites, the P450 redox partner gene Cpr showed evidence for contemporary selection (H12) and population differentiation (FST) meriting further exploration as a potential resistance associated marker.

Mitochondrial Variation in Anopheles gambiae and Anopheles coluzzii: Phylogeographic Legacy and Mitonuclear Associations With Metabolic Resistance to Pathogens and Insecticides

Mitochondrial DNA has been a popular marker in phylogeography, phylogeny, and molecular ecology, but its complex evolution is increasingly recognized. Here, we investigated mitochondrial DNA variation in Anopheles gambiae and Anopheles coluzzii, in relation to other species in the Anopheles gambiae complex, by assembling the mitogenomes of 1,219 mosquitoes across Africa. The mitochondrial DNA phylogeny of the Anopheles gambiae complex was consistent with previously reported highly reticulated evolutionary history, revealing important discordances with the species tree. The three most widespread species (An. gambiae, An. coluzzii, and Anopheles arabiensis), known for extensive historical introgression, could not be discriminated based on mitogenomes. Furthermore, a monophyletic clustering of the three saltwater-tolerant species (Anopheles merus, Anopheles melas, and Anopheles bwambae) in the Anopheles gambiae complex also suggested that introgression and possibly selection shaped mitochondrial DNA evolution. Mitochondrial DNA variation in An. gambiae and An. coluzzii across Africa revealed significant partitioning among populations and species. A peculiar mitochondrial DNA lineage found predominantly in An. coluzzii and in the hybrid taxon of the African "far-west" exhibited divergence comparable to the interspecies divergence in the Anopheles gambiae complex, with a geographic distribution matching closely An. coluzzii's geographic range. This phylogeographic relict of the An. coluzzii and An. gambiae split was associated with population and species structure, but not with the rare Wolbachia occurrence. The lineage was significantly associated with single nucleotide polymorphisms in the nuclear genome, particularly in genes associated with pathogen and insecticide resistance. These findings underline potential mitonuclear coevolution history and the role played by mitochondria in shaping metabolic responses to pathogens and insecticides in Anopheles.

Copy number variants underlie the major selective sweeps in insecticide resistance genes in Anopheles arabiensis from Tanzania

To keep ahead of the evolution of resistance to insecticides in mosquitoes, national malaria control programmes must make use of a range of insecticides, both old and new, while monitoring resistance mechanisms. Knowledge of the mechanisms of resistance remains limited in Anopheles arabiensis, which in many parts of Africa is of increasing importance because it is apparently less susceptible to many indoor control interventions. Furthermore, comparatively little is known in general about resistance to non-pyrethroid insecticides such as pirimiphos-methyl (PM), which are crucial for effective control in the context of resistance to pyrethroids. We performed a genome-wide association study to determine the molecular mechanisms of resistance to deltamethrin (commonly used in bednets) and PM, in An. arabiensis from two regions in Tanzania. Genomic regions of positive selection in these populations were largely driven by copy number variants (CNVs) in gene families involved in resistance to these two insecticides. We found evidence of a new gene cluster involved in resistance to PM, identifying a strong selective sweep tied to a CNV in the Coeae2g-Coeae6g cluster of carboxylesterase genes. Using complementary data from An. coluzzii in Ghana, we show that copy number at this locus is significantly associated with PM resistance. Similarly, for deltamethrin, resistance was strongly associated with a novel CNV allele in the Cyp6aa / Cyp6p cluster. Against this background of metabolic resistance, target site resistance was very rare or absent for both insecticides. Mutations in the pyrethroid target site Vgsc were at very low frequency in Tanzania, yet combining these samples with three An. arabiensis individuals from West Africa revealed a startling diversity of evolutionary origins of target site resistance, with up to 5 independent origins of Vgsc-995 mutations found within just 8 haplotypes. Thus, despite having been first recorded over 10 years ago, Vgsc resistance mutations in Tanzanian An. arabiensis have remained at stable low frequencies. Overall, our results provide a new copy number marker for monitoring resistance to PM in malaria mosquitoes, and reveal the complex picture of resistance patterns in An. arabiensis.

Analysis of geometric morphometrics and molecular phylogeny for Anopheles species in the Republic of Korea

Human malaria, transmitted by Anopheles mosquitoes, is the most predominant mosquito-borne disease that is responsible for hundreds of thousands of deaths worldwide each year. In the Republic of Korea (ROK), there are currently several hundred malaria cases annually, mostly near the demilitarized zone (DMZ). Eight species of Anopheles mosquitoes are currently known to be present in the ROK. Similar to other major malaria vectors in Africa and India, it is very challenging to morphologically differentiate Anopheles mosquitoes in the ROK due to their extremely similar morphology. In this study, wing geometric morphometrics (WGM) were used to differentiate the eight Anopheles species collected at six locations near the DMZ, Seoul and Pyeongtaek from April-October 2021. Phylogenetic analysis was also performed using cytochrome c oxidase subunit 1 (COI), internal transcribed spacer 2 (ITS2), and tyrosine hydroxylase (TH) genes for comparison with WGM analysis and to infer evolutionary relationships. The results of cross-validation (overall accuracy = 74.8%) demonstrated that species identification using WGM alone was not possible with a high accuracy for all eight species. While phylogenetic analyses based on the COI region could not clearly distinguish some species, the analysis based on ITS2 and TH was more useful for resolving the phylogenetic correlation of the eight species. Our results may improve Anopheles species identification strategies for effective identification and control of malaria vectors in the ROK.

Genome-wide association studies reveal novel loci associated with pyrethroid and organophosphate resistance in Anopheles gambiae and Anopheles coluzzii

Resistance to insecticides in Anopheles mosquitoes threatens the effectiveness of malaria control, but the genetics of resistance are only partially understood. We performed a large scale multi-country genome-wide association study of resistance to two widely used insecticides: deltamethrin and pirimiphos-methyl, using sequencing data from An. gambiae and An. coluzzii from ten locations in West Africa. Resistance was highly multi-genic, multi-allelic and variable between populations. While the strongest and most consistent association with deltamethrin resistance came from Cyp6aa1, this was based on several independent copy number variants (CNVs) in An. coluzzii, and on a non-CNV haplotype in An. gambiae. For pirimiphos-methyl, signals included Ace1, cytochrome P450s, glutathione S-transferases and the nAChR target site of neonicotinoid insecticides. The regions around Cyp9k1 and the Tep family of immune genes showed evidence of cross-resistance to both insecticides. These locally-varying, multi-allelic patterns highlight the challenges involved in genomic monitoring of resistance, and may form the basis for improved surveillance methods.

The first genome sequence of Anopheles squamous from Macha, Zambia

Despite efforts to minimize the impacts of malaria and reduce the number of primary vectors, malaria has yet to be eliminated in Zambia. Understudied vector species may perpetuate malaria transmission in pre-elimination settings. Anopheles squamosus is one of the most abundantly caught mosquito species in southern Zambia and has previously been found with Plasmodium falciparum sporozoites, a causal agent of human malaria. This species may be a critical vector of malaria transmission, however, there is a lack of genetic information available for An. squamosus. We report the first genome data and the first complete mitogenome (Mt) sequence of An. squamosus. The sequence was extracted from one individual mosquito from the Chidakwa area in Macha, Zambia. The raw reads were obtained using Illumina Novaseq 6000 and assembled through NOVOplasty alignment with related species. The length of the An. squamosus Mt was 15,351 bp, with 77.9 % AT content. The closest match to the whole mitochondrial genome in the phylogenetic tree is the African malaria mosquito, Anopheles gambiae. Its genome data is available through National Center for Biotechnology Information (NCBI) Sequencing Reads Archive (SRA) with accession number SRR22114392. The mitochondrial genome was deposited in NCBI GenBank with the accession number OP776919. The ITS2 containing contig sequence was deposited in GenBank with the accession number OQ241725. Mitogenome annotation and a phylogenetic tree with related Anopheles mosquito species are provided.

Genome-wide association studies reveal novel loci associated with pyrethroid and organophosphate resistance in Anopheles gambiae s.l

Resistance to insecticides in Anopheles mosquitoes threatens the effectiveness of the most widespread tools currently used to control malaria. The genetic underpinnings of resistance are still only partially understood, with much of the variance in resistance phenotype left unexplained. We performed a multi-country large scale genome-wide association study of resistance to two insecticides widely used in malaria control: deltamethrin and pirimiphos-methyl. Using a bioassay methodology designed to maximise the phenotypic difference between resistant and susceptible samples, we sequenced 969 phenotyped female An. gambiae and An. coluzzii from ten locations across four countries in West Africa (Benin, Côte d'Ivoire, Ghana and Togo), identifying single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) segregating in the populations. The patterns of resistance association were highly multiallelic and variable between populations, with different genomic regions contributing to resistance, as well as different mutations within a given region. While the strongest and most consistent association with deltamethrin resistance came from the region around Cyp6aa1 , this resistance was based on a combination of several independent CNVs in An. coluzzii , and on a non-CNV bearing haplotype in An. gambiae . Further signals involved a range of cytochrome P450, mitochondrial, and immunity genes. Similarly, for pirimiphos-methyl, while the strongest signal came from the region of Ace1 , more widespread signals included cytochrome P450s, glutathione S-transferases, and a subunit of the nAChR target site of neonicotinoid insecticides. The regions around Cyp9k1 and the Tep family of immune genes were associated with resistance to both insecticide classes, suggesting possible cross-resistance mechanisms. These locally-varying, multigenic and multiallelic patterns highlight the challenges involved in genomic monitoring and surveillance of resistance, and form the basis for improvement of methods used to detect and predict resistance. Based on simulations of resistance variants, we recommend that yet larger scale studies, exceeding 500 phenotyped samples per population, are required to better identify associated genomic regions.

Outdoor residual spraying for malaria vector-control in Kayin (Karen) state, Myanmar: A cluster randomized controlled trial

Outdoor and early biting by mosquitoes challenge the efficacy of bed nets and indoor residual spraying against malaria in the Greater Mekong Subregion. The objective of this study was to assess the efficacy of outdoor residual spraying (ORS) for malaria vector-control in this region. A cluster randomized controlled trial was conducted between July 2018 and April 2019 in twelve villages in Karen (Kayin) state, Myanmar. Villages were randomly assigned to receive either a single round of ORS with a capsule suspension of lambda-cyhalothrin for two days in October or no intervention (six villages per group). The primary endpoint was the biting rate of malaria mosquitoes assessed with human-landing catch and cow-baited trap collection methods, and was analyzed with a Bayesian multi-level model. In the intervention villages, the proportion of households located within the sprayed area ranged between 42 and 100% and the application rate ranged between 63 and 559 g of active ingredient per hectare. At baseline, the median of Anopheles biting rate estimates in the twelve villages was 2 bites per person per night (inter-quartile range [IQR] 0–5, range 0–48) indoors, 6 bites per person per night (IQR 2–16, range 0–342) outdoors and 206 bites per cow per night (IQR 83–380, range 19–1149) in the cow-baited trap. In intention-to-treat analysis, it was estimated that ORS reduced biting rate by 72% (95% confidence interval [CI] 63–79) from Month 0 to Month 3 and by 79% (95% CI 62–88) from Month 4 to Month 6, considering control villages as the reference. In conclusion, ORS rapidly reduces the biting rates of malaria mosquitoes in a Southeast Asian setting where the vectors bite mostly outdoors and at a time when people are not protected by mosquito bed nets.

Multi-omics analysis identifies a CYP9K1 haplotype conferring pyrethroid resistance in the malaria vector Anopheles funestus in East Africa

Metabolic resistance to pyrethroids is a menace to the continued effectiveness of malaria vector controls. Its molecular basis is complex and varies geographically across Africa. Here, we used a multi‐omics approach, followed‐up with functional validation to show that a directionally selected haplotype of a cytochrome P450, CYP9K1 is a major driver of resistance in Anopheles funestus. A PoolSeq GWAS using mosquitoes alive and dead after permethrin exposure, from Malawi and Cameroon, detected candidate genomic regions, but lacked consistency across replicates. Targeted sequencing of candidate resistance genes detected several SNPs associated with known pyrethroid resistance QTLs. The most significant SNPs were in the cytochrome P450 CYP304B1 (Cameroon), CYP315A1 (Uganda) and the ABC transporter gene ABCG4 (Malawi). However, when comparing field resistant mosquitoes to laboratory susceptible, the pyrethroid resistance locus rp1 and SNPs around the ABC transporter ABCG4 were consistently significant, except for Uganda where SNPs in the P450 CYP9K1 was markedly significant. In vitro heterologous metabolism assays with recombinant CYP9K1 revealed that it metabolises type II pyrethroid (deltamethrin; 64% depletion) but not type I (permethrin; 0%), while moderately metabolising DDT (17%). CYP9K1 exhibited reduced genetic diversity in Uganda underlying an extensive selective sweep. Furthermore, a glycine to alanine (G454A) amino acid change in CYP9K1 was fixed in Ugandan mosquitoes but not in other An. funestus populations. This study sheds further light on the evolution of metabolic resistance in a major malaria vector by implicating more genes and variants that can be used to design field‐applicable markers to better track resistance Africa‐wide.

Multiple Novel Clades of Anopheline Mosquitoes Caught Outdoors in Northern Zambia

Residual vector populations that do not come in contact with the most frequently utilized indoor-directed interventions present major challenges to global malaria eradication. Many of these residual populations are mosquito species about which little is known. As part of a study to assess the threat of outdoor exposure to malaria mosquitoes within the Southern and Central Africa International Centers of Excellence for Malaria Research, foraging female anophelines were collected outside households in Nchelenge District, northern Zambia. These anophelines proved to be more diverse than had previously been reported in the area. In order to further characterize the anopheline species, sequencing and phylogenetic approaches were utilized. Anopheline mosquitoes were collected from outdoor light traps, morphologically identified, and sent to Johns Hopkins Bloomberg School of Public Health for sequencing. Sanger sequencing from 115 field-derived samples yielded mitochondrial COI sequences, which were aligned with a homologous 488 bp gene segment from known anophelines (n = 140) retrieved from NCBI. Nuclear ITS2 sequences (n = 57) for at least one individual from each unique COI clade were generated and compared against NCBI’s nucleotide BLAST database to provide additional evidence for taxonomical identity and structure. Molecular and morphological data were combined for assignment of species or higher taxonomy. Twelve phylogenetic groups were characterized from the COI and ITS2 sequence data, including the primary vector species Anopheles funestus s.s. and An. gambiae s.s. An unexpectedly large proportion of the field collections were identified as An. coustani and An. sp. 6. Six phylogenetic groups remain unidentified to species-level. Outdoor collections of anopheline mosquitoes in areas frequented by people in Nchelenge, northern Zambia, proved to be extremely diverse. Morphological misidentification and underrepresentation of some anopheline species in sequence databases confound efforts to confirm identity of potential malaria vector species. The large number of unidentified anophelines could compromise the malaria vector surveillance and malaria control efforts not only in northern Zambia but other places where surveillance and control are focused on indoor-foraging and resting anophelines. Therefore, it is critical to continue development of methodologies that allow better identification of these populations and revisiting and cleaning current genomic databases.

Assessing cross-resistance within the pyrethroids in terms of their interactions with key cytochrome P450 enzymes and resistance in vector populations

Background

It is important to understand whether the potential impact of pyrethroid resistance on malaria control can be mitigated by switching between different pyrethroids or whether cross-resistance within this insecticide class precludes this approach.

Methods

Here we assess the relationships among pyrethroids in terms of their binding affinity to, and depletion by, key cytochrome P450 enzymes (hereafter P450s) that are known to confer metabolic pyrethroid resistance in Anopheles gambiae (s.l.) and An. funestus, in order to identify which pyrethroids may diverge from the others in their vulnerability to resistance. We then investigate whether these same pyrethroids also diverge from the others in terms of resistance in vector populations.

Results

We found that the type I and II pyrethroids permethrin and deltamethrin, respectively, are closely related in terms of binding affinity to key P450s, depletion by P450s and resistance within vector populations. Bifenthrin, which lacks the common structural moiety of most pyrethroids, diverged from the other pyrethroids tested in terms of both binding affinity to key P450s and depletion by P450s, but resistance to bifenthrin has rarely been tested in vector populations and was not analysed here. Etofenprox, which also lacks the common structural moiety of most pyrethroids, diverged from the more commonly deployed pyrethroids in terms of binding affinity to key P450s and resistance in vector populations, but did not diverge from these pyrethroids in terms of depletion by the P450s. The analysis of depletion by the P450s indicated that etofenprox may be more vulnerable to metabolic resistance mechanisms in vector populations. In addition, greater resistance to etofenprox was found across Aedes aegypti populations, but greater resistance to this compound was not found in any of the malaria vector species analysed. The results for pyrethroid depletion by anopheline P450s in the laboratory were largely not repeated in the findings for resistance in malaria vector populations.

Conclusion

Importantly, the prevalence of resistance to the pyrethroids α-cypermethrin, cyfluthrin, deltamethrin, λ-cyhalothrin and permethrin was correlated across malaria vector populations, and switching between these compounds as a tool to mitigate against pyrethroid resistance is not advised without strong evidence supporting a true difference in resistance.

Laboratory and microcosm experiments reveal contrasted adaptive responses to ammonia and water mineralisation in aquatic stages of the sibling species Anopheles gambiae (sensu stricto) and Anopheles coluzzii

Background

The sibling species of the malaria mosquito, Anopheles gambiae (sensu stricto) and Anopheles coluzzii co-exist in many parts of West Africa and are thought to have recently diverged through a process of ecological speciation with gene flow. Divergent larval ecological adaptations, resulting in Genotype-by-Environment (G × E) interactions, have been proposed as important drivers of speciation in these species. In West Africa, An. coluzzii tends to be associated with permanent man-made larval habitats such as irrigated rice fields, which are typically more eutrophic and mineral and ammonia-rich than the temporary rain pools exploited by An. gambiae (s.s.)

Methods

To highlight G × E interactions at the larval stage and their possible role in ecological speciation of these species, we first investigated the effect of exposure to ammonium hydroxide and water mineralisation on larval developmental success. Mosquito larvae were exposed to two water sources and increasing ammonia concentrations in small containers until adult emergence. In a second experiment, larval developmental success was compared across two contrasted microcosms to highlight G × E interactions under conditions such as those found in the natural environment.

Results

The first experiment revealed significant G × E interactions in developmental success and phenotypic quality for both species in response to increasing ammonia concentrations and water mineralisation. The An. coluzzii strain outperformed the An. gambiae (s.s.) strain under limited conditions that were closer to more eutrophic habitats. The second experiment revealed divergent crisscrossing reaction norms in the developmental success of the sibling species in the two contrasted larval environments. As expected, An. coluzzii had higher emergence rates in the rice paddy environment with emerging adults of superior phenotypic quality compared to An. gambiae (s.s.), and vice versa, in the rain puddle environment.

Conclusions

Evidence for such G × E interactions lends support to the hypothesis that divergent larval adaptations to the environmental conditions found in man-made habitats such as rice fields in An. coluzzii may have been an important driver of its ecological speciation.

A population genomic unveiling of a new cryptic mosquito taxon within the malaria-transmitting Anopheles gambiae complex

The Anopheles gambiae complex consists of multiple morphologically indistinguishable mosquito species including the most important vectors of the malaria parasite Plasmodium falciparum in sub-Saharan Africa. Nine cryptic species have been described so far within the complex. The ecological, immunological and reproductive differences among these species will critically impact population responses to disease control strategies and environmental changes. Here, we examine whole-genome sequencing data from a longitudinal study of putative A. coluzzii in western Burkina Faso. Surprisingly, many specimens are genetically divergent from A. coluzzii and all other Anopheles species and represent a new taxon, here designated Anopheles TENGRELA (AT). Population genetic analysis suggests that the cryptic GOUNDRY subgroup, previously collected as larvae in central Burkina Faso, represents an admixed population descended from both A. coluzzii and AT. AT harbours low nucleotide diversity except for the 2La inversion polymorphism which is maintained by overdominance. It shows numerous fixed differences with A. coluzzii concentrated in several regions reflecting selective sweeps, but the two taxa are identical at standard diagnostic loci used for taxon identification, and thus, AT may often go unnoticed. We present an amplicon-based genotyping assay for identifying AT which could be usefully applied to numerous existing samples. Misidentified cryptic taxa could seriously confound ongoing studies of Anopheles ecology and evolution in western Africa, including phenotypic and genotypic surveys of insecticide resistance. Reproductive barriers between cryptic species may also complicate novel vector control efforts, for example gene drives, and hinder predictions about evolutionary dynamics of Anopheles and Plasmodium.

Spatial Anopheles arabiensis (Diptera: Culicidae) insecticide resistance patterns across malaria-endemic regions of Botswana

Background

Since the advent of the Green Revolution, pesticides have played an important role in the global management of invertebrate pests including vector mosquitoes. Despite optimal efficacy, insects often display insensitivity to synthetic insecticides owing to prolonged exposure that may select for resistance development. Such insecticide insensitivity may regress national and regional coordination in mosquito vector management and indeed malaria control. In Botswana, prolonged use of synthetic insecticides against malaria vectors have been practiced without monitoring of targeted mosquito species susceptibility status.

Methods

Here, susceptibility status of a malaria vector (Anopheles arabiensis), was assessed against World Health Organization-recommended insecticides, across three malaria endemic districts. Adult virgin female mosquitoes (2–5 days old) emerging from wild-collected larvae were exposed to standardized insecticide-impregnated papers with discriminating doses.

Results

The results showed resistance dynamics were variable in space, presumably as a result of spatial differences in insecticide use across malaria endemic districts and the types of insecticides used in the country. Overall, there was a reduced susceptibility of An. arabiensis for the pyrethroid lambda-cyhalothrin and for dichloro diphenyl trichloroethane [DDT], which have similar modes of action and have been used in the country for many years. The Okavango district exhibited the greatest reduction in susceptibility, followed by Ngamiland and then Bobirwa, reflective of national intervention strategy spraying intensities. Vector mosquitoes were, however, highly susceptible to carbamates and organophosphates irrespective of region.

Conclusions

These results provide important findings of vector susceptibility to insecticides recommended for vector control. The results highlight the need to implement insecticide application regimes that more effectively including regionally integrated resistance management strategies for effective malaria control and elimination.

Detecting inversions with PCA in the presence of population structure

Chromosomal inversions can lead to reproductive isolation and adaptation in insects such as Drosophila melanogaster and the non-model malaria vector Anopheles gambiae. Inversions can be detected and characterized using principal component analysis (PCA) of single nucleotide polymorphisms (SNPs). To aid in developing such methods, we formed a new benchmark derived from three publicly-available insect data. We then used this benchmark to perform an extended validation of our software for inversion analysis (Asaph). Through that process, we identified and characterized several problematic test cases liable to misinterpretation that can help guide PCA-based inversion detection. Lastly, we re-analyzed the 2R chromosome arm of 150 An. gambiae and coluzzii samples and observed two inversions (2Rc and 2Rd) that were previously known but not annotated in these particular individuals. The resulting benchmark data set and methods will be useful for future inversion detection based solely on SNP data.

Inversion Genotyping in the Anopheles gambiae Complex Using High-Throughput Array and Sequencing Platforms

Chromosomal inversion polymorphisms have special importance in the Anopheles gambiae complex of malaria vector mosquitoes, due to their role in local adaptation and range expansion. The study of inversions in natural populations is reliant on polytene chromosome analysis by expert cytogeneticists, a process that is limited by the rarity of trained specialists, low throughput, and restrictive sampling requirements. To overcome this barrier, we ascertained tag single nucleotide polymorphisms (SNPs) that are highly correlated with inversion status (inverted or standard orientation). We compared the performance of the tag SNPs using two alternative high throughput molecular genotyping approaches vs. traditional cytogenetic karyotyping of the same 960 individual An. gambiae and An. coluzzii mosquitoes sampled from Burkina Faso, West Africa. We show that both molecular approaches yield comparable results, and that either one performs as well or better than cytogenetics in terms of genotyping accuracy. Given the ability of molecular genotyping approaches to be conducted at scale and at relatively low cost without restriction on mosquito sex or developmental stage, molecular genotyping via tag SNPs has the potential to revitalize research into the role of chromosomal inversions in the behavior and ongoing adaptation of An. gambiae and An. coluzzii to environmental heterogeneities.

The need for new vector control approaches targeting outdoor biting Anopheline malaria vector communities

Since the implementation of Roll Back Malaria, the widespread use of insecticide-treated nets (ITNs) and indoor residual spraying (IRS) is thought to have played a major part in the decrease in mortality and morbidity achieved in malaria-endemic regions. In the past decade, resistance to major classes of insecticides recommended for public health has spread across many malaria vector populations. Increasingly, malaria vectors are also showing changes in vector behaviour in response to current indoor chemical vector control interventions. Changes in the time of biting and proportion of indoor biting of major vectors, as well as changes in the species composition of mosquito communities threaten the progress made to control malaria transmission. Outdoor biting mosquito populations contribute to malaria transmission in many parts of sub-Saharan Africa and pose new challenges as they cannot be reliably monitored or controlled using conventional tools. Here, we review existing and novel approaches that may be used to target outdoor communities of malaria vectors. We conclude that scalable tools designed specifically for the control and monitoring of outdoor biting and resting malaria vectors with increasingly complex and dynamic responses to intensifying malaria control interventions are urgently needed. These are crucial for integrated vector management programmes designed to challenge current and future vector populations.

In Silico Karyotyping of Chromosomally Polymorphic Malaria Mosquitoes in the Anopheles gambiae Complex

Chromosomal inversion polymorphisms play an important role in adaptation to environmental heterogeneities. For mosquito species in the Anopheles gambiae complex that are significant vectors of human malaria, paracentric inversion polymorphisms are abundant and are associated with ecologically and epidemiologically important phenotypes. Improved understanding of these traits relies on determining mosquito karyotype, which currently depends upon laborious cytogenetic methods whose application is limited both by the requirement for specialized expertise and for properly preserved adult females at specific gonotrophic stages. To overcome this limitation, we developed sets of tag single nucleotide polymorphisms (SNPs) inside inversions whose biallelic genotype is strongly correlated with inversion genotype. We leveraged 1,347 fully sequenced An. gambiae and Anopheles coluzzii genomes in the Ag1000G database of natural variation. Beginning with principal components analysis (PCA) of population samples, applied to windows of the genome containing individual chromosomal rearrangements, we classified samples into three inversion genotypes, distinguishing homozygous inverted and homozygous uninverted groups by inclusion of the small subset of specimens in Ag1000G that are associated with cytogenetic metadata. We then assessed the correlation between candidate tag SNP genotypes and PCA-based inversion genotypes in our training sets, selecting those candidates with >80% agreement. Our initial tests both in held-back validation samples from Ag1000G and in data independent of Ag1000G suggest that when used for in silico inversion genotyping of sequenced mosquitoes, these tags perform better than traditional cytogenetics, even for specimens where only a small subset of the tag SNPs can be successfully ascertained.

Whole-genome sequencing reveals high complexity of copy number variation at insecticide resistance loci in malaria mosquitoes

Polymorphisms in genetic copy number can influence gene expression, coding sequence, and zygosity, making them powerful actors in the evolutionary process. Copy number variants (CNVs) are however understudied, being more difficult to detect than single-nucleotide polymorphisms. We take advantage of the intense selective pressures on the major malaria vector Anopheles gambiae, caused by the widespread use of insecticides for malaria control, to investigate the role of CNVs in the evolution of insecticide resistance. Using the whole-genome sequencing data from 1142 samples in the An. gambiae 1000 genomes project, we identified 250 gene-containing CNVs, encompassing a total of 267 genes of which 28 were in gene families linked to metabolic insecticide resistance, representing significant enrichment of these families. The five major gene clusters for metabolic resistance all contained CNVs, with 44 different CNVs being found across these clusters and multiple CNVs frequently covering the same genes. These 44 CNVs are widespread (45% of individuals carry at least one of them) and have been spreading through positive selection, indicated by their high local frequencies and extended haplotype homozygosity. Our results demonstrate the importance of CNVs in the response to selection, highlighting the urgent need to identify the contribution of each CNV to insecticide resistance and to track their spread as the use of insecticides in malaria endemic countries intensifies and as the operational deployment of next-generation bed nets targeting metabolic resistance gathers pace. Our detailed descriptions of CNVs found across the species range provide the tools to do so.

Introgression between Anopheles gambiae and Anopheles coluzzii in Burkina Faso and its associations with kdr resistance and Plasmodium infection

BACKGROUND

Insecticide resistance in Anopheles coluzzii mosquitoes has become widespread throughout West Africa including in Burkina Faso. The insecticide resistance allele (kdr or L1014F) is a prime indicator that is highly correlated with phenotypic resistance in West Africa. Studies from Benin, Ghana and Mali have suggested that the source of the L1014F is introgression of the 2L divergence island via interspecific hybridization with Anopheles gambiae. The goal of this study was to characterize local mosquito populations in the Nouna Department, Burkina Faso with respect to: (i) the extent of introgression between An. coluzzii and An. gambiae, (ii) the frequency of the L1014F mutation and (iii) Plasmodium infection rates.

METHODS

A total of 95 mosquitoes were collected from ten sites surrounding Nouna town in Kossi Province, Burkina Faso in 2012. The species composition, the extent of introgression in An. coluzzii mosquitoes and their Plasmodium infection rates were identified with a modified version of the "Divergence Island SNP" (DIS) genotyping assay.

RESULTS

The mosquito collection contained 70.5% An. coluzzii, 89.3% of which carried a 3 Mb genomic region on the 2L chromosome with L1014F insecticide resistance mutation that was introgressed from An. gambiae. In addition, 22.4% in the introgressed An. coluzzii specimens were infected with Plasmodium falciparum, whereas none of the non-introgressed ("pure") An. coluzzii were infected.

CONCLUSION

This paper is the first report providing divergence island SNP genotypes for natural population of Burkina Faso and corresponding Plasmodium infection rates. These observations warrant further study and could have a major impact on future malaria control strategies in Burkina Faso.

Genome-wide divergence among invasive populations of Aedes aegypti in California

BACKGROUND

In the summer of 2013, Aedes aegypti Linnaeus was first detected in three cities in central California (Clovis, Madera and Menlo Park). It has now been detected in multiple locations in central and southern CA as far south as San Diego and Imperial Counties. A number of published reports suggest that CA populations have been established from multiple independent introductions.

RESULTS

Here we report the first population genomics analyses of Ae. aegypti based on individual, field collected whole genome sequences. We analyzed 46 Ae. aegypti genomes to establish genetic relationships among populations from sites in California, Florida and South Africa. Based on 4.65 million high quality biallelic SNPs, we identified 3 major genetic clusters within California; one that includes all sample sites in the southern part of the state (South of Tehachapi mountain range) plus the town of Exeter in central California and two additional clusters in central California.

CONCLUSIONS

A lack of concordance between mitochondrial and nuclear genealogies suggests that the three founding populations were polymorphic for two main mitochondrial haplotypes prior to being introduced to California. One of these has been lost in the Clovis populations, possibly by a founder effect. Genome-wide comparisons indicate extensive differentiation between genetic clusters. Our observations support recent introductions of Ae. aegypti into California from multiple, genetically diverged source populations. Our data reveal signs of hybridization among diverged populations within CA. Genetic markers identified in this study will be of great value in pursuing classical population genetic studies which require larger sample sizes.

Complete Anopheles funestus mitogenomes reveal an ancient history of mitochondrial lineages and their distribution in southern and central Africa

Anopheles funestus s.s. is a primary vector of malaria in sub-Saharan Africa. Despite its important role in human Plasmodium transmission, evolutionary history, genetic diversity, and population structure of An. funestus in southern and central Africa remains understudied. We deep sequenced, assembled, and annotated the complete mitochondrial genome of An. funestus s.s. for the first time, providing a foundation for further genetic research of this important malaria vector species. We further analyzed the complete mitochondrial genomes of 43 An. funestus s.s. from three sites in Zambia, Democratic Republic of the Congo, and Tanzania. From these 43 mitogenomes we identified 41 unique haplotypes that comprised 567 polymorphic sites. Bayesian phylogenetic reconstruction confirmed the co-existence of two highly divergent An. funestus maternal lineages, herein defined as lineages I and II, in Zambia and Tanzania. The estimated coalescence time of these two mitochondrial lineages is ~500,000 years ago (95% HPD 426,000–594,000 years ago) with subsequent independent diversification. Haplotype network and phylogenetic analysis revealed two major clusters within lineage I, and genetic relatedness of samples with deep branching in lineage II. At this time, data suggest that the lineages are partially sympatric. This study illustrates that accurate retrieval of full mitogenomes of Anopheles vectors enables fine-resolution studies of intraspecies genetic relationships, population differentiation, and demographic history. Further investigations on whether An. funestus mitochondrial lineages represent biologically meaningful populations and their potential implications for malaria vector control are warranted.

Rapid selection of a pyrethroid metabolic enzyme CYP9K1 by operational malaria control activities

Since 2004, indoor residual spraying (IRS) and long-lasting insecticide-impregnated bednets (LLINs) have reduced the malaria parasite prevalence in children on Bioko Island, Equatorial Guinea, from 45% to 12%. After target site-based (knockdown resistance; kdr) pyrethroid resistance was detected in 2004 in Anopheles coluzzii (formerly known as the M form of the Anopheles gambiae complex), the carbamate bendiocarb was introduced. Subsequent analysis showed that kdr alone was not operationally significant, so pyrethroid-based IRS was successfully reintroduced in 2012. In 2007 and 2014-2015, mass distribution of new pyrethroid LLINs was undertaken to increase the net coverage levels. The combined selection pressure of IRS and LLINs resulted in an increase in the frequency of pyrethroid resistance in 2015. In addition to a significant increase in kdr frequency, an additional metabolic pyrethroid resistance mechanism had been selected. Increased metabolism of the pyrethroid deltamethrin was linked with up-regulation of the cytochrome P450 CYP9K1. The increase in resistance prompted a reversion to bendiocarb IRS in 2016 to avoid a resurgence of malaria, in line with the national Malaria Control Program plan.

Genetic variation associated with increased insecticide resistance in the malaria mosquito, Anopheles coluzzii

BACKGROUND

Malaria mortality rates in sub-Saharan Africa have declined significantly in recent years as a result of increased insecticide-treated bed net (ITN) usage. A major challenge to further progress is the emergence and spread of insecticide resistance alleles in the Anopheles mosquito vectors, like An. coluzzii. A non-synonymous mutation in the para voltage-gated sodium channel gene reduces pyrethroid-binding affinity, resulting in knockdown resistance (kdr). Metabolic mechanisms of insecticide resistance involving detoxification genes like cytochrome P450 genes, carboxylesterases, and glutathione S-transferases are also important. As some gene activity is tissue-specific and/or environmentally induced, gene regulatory variation may be overlooked when comparing expression from whole mosquito bodies under standard rearing conditions.

RESULTS

We detected complex insecticide resistance in a 2014 An. coluzzii colony from southern Mali using bottle bioassays. Additional bioassays involving recombinant genotypes from a cross with a relatively susceptible 1995 An. coluzzii colony from Mali confirmed the importance of kdr and associated increased permethrin resistance to the CYP9K1 locus on the X chromosome. Significant differential expression of CYP9K1 was not observed among these colonies in Malpighian tubules. However, the P450 gene CYP6Z1 was overexpressed in resistant individuals following sublethal permethrin exposure and the carboxylesterase gene COEAE5G was constitutively overexpressed.

CONCLUSIONS

The significant P450-related insecticide resistance observed in the 2014 An. coluzzii colony indicates that ITNs treated with the P450 inhibitor piperonyl butoxide (PBO) would be more effective in this region. The known insecticide resistance gene CYP6Z1 was differentially expressed exclusively in the context of sublethal permethrin exposure, highlighting the importance of tissue-specificity and environmental conditions in gene expression studies. The increased activity of the carboxylesterase COEAE5G in the resistant An. coluzzii colony suggests resistance to other insecticides like organophosphates. Additional gene expression studies involving other tissues (e.g. fat body) would provide a more comprehensive view of genes underlying metabolic insecticide resistance in An. coluzzii from Mali. Identifying genetic markers linked to these regulatory alleles is an important next step that would substantially improve insecticide resistance surveillance and population genetic studies in this important vector species.

Mitochondrial genomes of Anopheles arabiensis, An. gambiae and An. coluzzii show no clear species division

Here we report the complete mitochondrial sequences of 70 individual field collected mosquito specimens from throughout Sub-Saharan Africa. We generated this dataset to identify species specific markers for the following Anopheles species and chromosomal forms: An. arabiensis, An. coluzzii (The Forest and Mopti chromosomal forms) and An. gambiae (The Bamako and Savannah chromosomal forms).  The raw Illumina sequencing reads were mapped to the NC_002084 reference mitogenome sequence. A total of 783 single nucleotide polymorphisms (SNPs) were detected on the mitochondrial genome, of which 460 are singletons (58.7%). None of these SNPs are suitable as molecular markers to distinguish among An. arabiensis, An. coluzzii and An. gambiae or any of the chromosomal forms. The lack of species or chromosomal form specific markers is also reflected in the constructed phylogenetic tree, which shows no clear division among the operational taxonomic units considered here.

Mitochondrial genomes of Anopheles arabiensis, An. gambiae and An. coluzzii show no clear species division

Here we report the complete mitochondrial sequences of 70 individual field collected mosquito specimens from throughout Sub-Saharan Africa. We generated this dataset to identify species specific markers for the following Anopheles species and chromosomal forms: An. arabiensis, An. coluzzii (The Forest and Mopti chromosomal forms) and An. gambiae (The Bamako and Savannah chromosomal forms).  The raw Illumina sequencing reads were mapped to the NC_002084 reference mitogenome sequence. A total of 783 single nucleotide polymorphisms (SNPs) were detected on the mitochondrial genome, of which 460 are singletons (58.7%). None of these SNPs are suitable as molecular markers to distinguish among An. arabiensis, An. coluzzii and An. gambiae or any of the chromosomal forms. The lack of species or chromosomal form specific markers is also reflected in the constructed phylogenetic tree, which shows no clear division among the operational taxonomic units considered here.

Human Interventions: Driving Forces of Mosquito Evolution

One of the most common strategies for controlling mosquito-borne diseases relies on the use of chemical pesticides to repel or kill the mosquito vector. Pesticide exposure interferes with several key biological functions in the mosquito and triggers a variety of adaptive responses whose underlying mechanisms are only partially elucidated. The availability of reference genome sequences opens up the possibility of tracking signatures of evolutionary changes, including the most recent, across the genomes of many vector species. In this review, we highlight the recent genomic changes, which contribute to the fascinating adaptation of malaria vectors of the sub-Saharan African region to intensive insecticide-based interventions. We emphasize the operational significance of detailed genomic knowledge for current monitoring and decision-making.

The transcription factor Maf-S regulates metabolic resistance to insecticides in the malaria vector Anopheles gambiae

Background

Malaria control in Africa is dependent upon the use insecticides but intensive use of a limited number of chemicals has led to resistance in mosquito populations. Increased production of enzymes that detoxify insecticides is one of the most potent resistance mechanisms. Several metabolic enzymes have been implicated in insecticide resistance but the processes controlling their expression have remained largely elusive.

Results

Here, we show that the transcription factor Maf-S regulates expression of multiple detoxification genes, including the key insecticide metabolisers CYP6M2 and GSTD1 in the African malaria vector Anopheles gambiae. Attenuation of this transcription factor through RNAi induced knockdown reduced transcript levels of these effectors and significantly increased mortality after exposure to the pyrethroid insecticides and DDT (permethrin: 9.2% to 19.2% (p = 0.015), deltamethrin: 3.9% to 21.6% (p = 0.036) and DDT: 1% to 11.7% (p = <0.01), whilst dramatically decreasing mortality induced by the organophosphate malathion (79.6% to 8.0% (p = <0.01)). Additional genes regulated by Maf-S were also identified providing new insight into the role of this transcription factor in insects.

Conclusion

Maf-S is a key regulator of detoxification genes in Anopheles mosquitoes. Disrupting this transcription factor has opposing effects on the mosquito’s response to different insecticide classes providing a mechanistic explanation to the negative cross resistance that has been reported between pyrethroids and organophosphates.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4086-7) contains supplementary material, which is available to authorized users.

Profiles of soluble proteins in chemosensory organs of three members of the afro-tropical Anopheles gambiae complex

In female mosquitoes, host-seeking and preference as well as several other important behaviors are largely driven by olfaction. Species of the Afrotropical Anopheles gambiae complex display divergent host-preference that are associated with significant differences in their vectorial capacity for human malaria. Olfactory sensitivity begins with signal transduction and activation of peripheral sensory neurons that populate the antennae, maxillary palps and other appendages. We have used shotgun proteomics to characterize the profile of soluble proteins of antennae and maxillary palps of three different species: An. coluzzii, An. arabiensis and An. quadriannulatus that display remarkable differences in anthropophilic behavior. This analysis revealed interspecific differences in the abundance of several proteins that comprise cuticular components, glutathione S-transferase and odorant binding proteins, the latter of which known to be directly involved in odor recognition.

Massive introgression drives species radiation at the range limit of Anopheles gambiae

Impacts of introgressive hybridisation may range from genomic erosion and species collapse to rapid adaptation and speciation but opportunities to study these dynamics are rare. We investigated the extent, causes and consequences of a hybrid zone between Anopheles coluzzii and Anopheles gambiae in Guinea-Bissau, where high hybridisation rates appear to be stable at least since the 1990s. Anopheles gambiae was genetically partitioned into inland and coastal subpopulations, separated by a central region dominated by A. coluzzii. Surprisingly, whole genome sequencing revealed that the coastal region harbours a hybrid form characterised by an A. gambiae-like sex chromosome and massive introgression of A. coluzzii autosomal alleles. Local selection on chromosomal inversions may play a role in this process, suggesting potential for spatiotemporal stability of the coastal hybrid form and providing resilience against introgression of medically-important loci and traits, found to be more prevalent in inland A. gambiae.

A Two-Locus Model of the Evolution of Insecticide Resistance to Inform and Optimise Public Health Insecticide Deployment Strategies

We develop a flexible, two-locus model for the spread of insecticide resistance applicable to mosquito species that transmit human diseases such as malaria. The model allows differential exposure of males and females, allows them to encounter high or low concentrations of insecticide, and allows selection pressures and dominance values to differ depending on the concentration of insecticide encountered. We demonstrate its application by investigating the relative merits of sequential use of insecticides versus their deployment as a mixture to minimise the spread of resistance. We recover previously published results as subsets of this model and conduct a sensitivity analysis over an extensive parameter space to identify what circumstances favour mixtures over sequences. Both strategies lasted more than 500 mosquito generations (or about 40 years) in 24% of runs, while in those runs where resistance had spread to high levels by 500 generations, 56% favoured sequential use and 44% favoured mixtures. Mixtures are favoured when insecticide effectiveness (their ability to kill homozygous susceptible mosquitoes) is high and exposure (the proportion of mosquitoes that encounter the insecticide) is low. If insecticides do not reliably kill homozygous sensitive genotypes, it is likely that sequential deployment will be a more robust strategy. Resistance to an insecticide always spreads slower if that insecticide is used in a mixture although this may be insufficient to outperform sequential use: for example, a mixture may last 5 years while the two insecticides deployed individually may last 3 and 4 years giving an overall 'lifespan' of 7 years for sequential use. We emphasise that this paper is primarily about designing and implementing a flexible modelling strategy to investigate the spread of insecticide resistance in vector populations and demonstrate how our model can identify vector control strategies most likely to minimise the spread of insecticide resistance.

Chromosomal inversions and ecotypic differentiation in Anopheles gambiae: the perspective from whole-genome sequencing

The molecular mechanisms and genetic architecture that facilitate adaptive radiation of lineages remain elusive. Polymorphic chromosomal inversions, due to their recombination-reducing effect, are proposed instruments of ecotypic differentiation. Here, we study an ecologically diversifying lineage of Anopheles gambiae, known as the Bamako chromosomal form based on its unique complement of three chromosomal inversions, to explore the impact of these inversions on ecotypic differentiation. We used pooled and individual genome sequencing of Bamako, typical (non-Bamako) An. gambiae and the sister species Anopheles coluzzii to investigate evolutionary relationships and genomewide patterns of nucleotide diversity and differentiation among lineages. Despite extensive shared polymorphism and limited differentiation from the other taxa, Bamako clusters apart from the other taxa, and forms a maximally supported clade in neighbour-joining trees based on whole-genome data (including inversions) or solely on collinear regions. Nevertheless, F_ST outlier analysis reveals that the majority of differentiated regions between Bamako and typical An. gambiae are located inside chromosomal inversions, consistent with their role in the ecological isolation of Bamako. Exceptionally differentiated genomic regions were enriched for genes implicated in nervous system development and signalling. Candidate genes associated with a selective sweep unique to Bamako contain substitutions not observed in sympatric samples of the other taxa, and several insecticide resistance gene alleles shared between Bamako and other taxa segregate at sharply different frequencies in these samples. Bamako represents a useful window into the initial stages of ecological and genomic differentiation from sympatric populations in this important group of malaria vectors.

A DNA extraction protocol for improved DNA yield from individual mosquitoes

Typical DNA extraction protocols from commercially available kits provide an adequate amount of DNA from a single individual mosquito sufficient for PCR-based assays. However, next-generation sequencing applications and high-throughput SNP genotyping assays exposed the limitation of DNA quantity one usually gets from a single individual mosquito. Whole genome amplification could alleviate the issue but it also creates bias in genome representation. While trying to find alternative DNA extraction protocols for improved DNA yield, we found that a combination of the tissue lysis protocol from Life Technologies and the DNA extraction protocol from Qiagen yielded a higher DNA amount than the protocol using the Qiagen or Life Technologies kit only. We have not rigorously tested all the possible combinations of extraction protocols; we also only tested this on mosquito samples. Therefore, our finding should be noted as a suggestion for improving people’s own DNA extraction protocols and not as an advertisement of a commercially available product.