BackgroundTo combat malaria transmission, the Ugandan government has embarked upon an ambitious programme of indoor residual spraying (IRS) with a carbamate class insecticide, bendiocarb. In preparation for this campaign, we characterized bendiocarb resistance and associated transcriptional variation among Anopheles gambiae s.s. mosquitoes from two sites in Uganda.ResultsGene expression in two mosquito populations displaying some resistance to bendiocarb (95% and 79% An. gambiae s.l. WHO tube bioassay mortality in Nagongera and Kihihi, respectively) was investigated using whole-genome microarrays. Significant overexpression of several genes encoding salivary gland proteins, including D7r2 and D7r4, was detected in mosquitoes from Nagongera. In Kihihi, D7r4, two detoxification-associated genes (Cyp6m2 and Gstd3) and an epithelial serine protease were among the genes most highly overexpressed in resistant mosquitoes. Following the first round of IRS in Nagongera, bendiocarb-resistant mosquitoes were collected, and real-time quantitative PCR analyses detected significant overexpression of D7r2 and D7r4 in resistant mosquitoes. A single nucleotide polymorphism located in a non-coding transcript downstream of the D7 genes was found at a significantly higher frequency in resistant individuals. In silico modelling of the interaction between D7r4 and bendiocarb demonstrated similarity between the insecticide and serotonin, a known ligand of D7 proteins. A meta-analysis of published microarray studies revealed a recurring association between D7 expression and insecticide resistance across Anopheles species and locations.ConclusionsA whole-genome microarray approach identified an association between novel insecticide resistance candidates and bendiocarb resistance in Uganda. In addition, a single nucleotide polymorphism associated with this resistance mechanism was discovered. The use of such impartial screening methods allows for discovery of resistance candidates that have no previously-ascribed function in insecticide binding or detoxification. Characterizing these novel candidates will broaden our understanding of resistance mechanisms and yield new strategies for combatting widespread insecticide resistance among malaria vectors.Electronic supplementary materialThe online version of this article (10.1186/s12864-018-4605-1) contains supplementary material, which is available to authorized users.

The evolution of insecticide resistance mechanisms in natural populations of Anopheles malaria vectors is a major public health concern across Africa. Using genome sequence data, we study the evolution of resistance mutations in the resistance to dieldrin locus (Rdl), a GABA receptor targeted by several insecticides, but most notably by the long-discontinued cyclodiene, dieldrin. The two Rdl resistance mutations (296G and 296S) spread across West and Central African Anopheles via two independent hard selective sweeps that included likely compensatory nearby mutations, and were followed by a rare combination of introgression across species (from A. gambiae and A. arabiensis to A. coluzzii) and across non-concordant karyotypes of the 2La chromosomal inversion. Rdl resistance evolved in the 1950s as the first known adaptation to a large-scale insecticide-based intervention, but the evolutionary lessons from this system highlight contemporary and future dangers for management strategies designed to combat development of resistance in malaria vectors.

Studies of insecticide resistance provide insights into the capacity of populations to show rapid evolutionary responses to contemporary selection. Malaria control remains heavily dependent on pyrethroid insecticides, primarily in long lasting insecticidal nets (LLINs). Resistance in the major malaria vectors has increased in concert with the expansion of LLIN distributions. Identifying genetic mechanisms underlying high-level resistance is crucial for the development and deployment of resistancebreaking tools. Using the Anopheles gambiae 1000 genomes (Ag1000g) data we identified a very recent selective sweep in mosquitoes from Uganda which localized to a cluster of cytochrome P450 genes. Further interrogation revealed a haplotype involving a trio of mutations, a nonsynonymous point mutation in Cyp6p4 (I236M), an upstream insertion of a partial Zanzibar-like transposable element (TE) and a duplication of the Cyp6aa1 gene. The mutations appear to have originated recently in An. gambiae

Insecticide resistance provides both an increasingly pressing threat to the control of vector-borne diseases and insights into the remarkable capacity of natural populations to show rapid evolutionary responses to contemporary selection. Malaria control remains heavily dependent on deployment of pyrethroid insecticides, primarily in long lasting insecticide treated nets (LLINs), but resistance in the major malaria vectors has increased over the last 15 years in concert with dramatic expansion of LLIN distributions. Identifying genetic mechanisms causing high-level resistance in mosquitoes, which may almost entirely overcome pyrethroid efficacy, is crucial for the development and deployment of potentially resistance-breaking tools. Using the Anopheles gambiae 1000 genomes (Ag1000G) data we identified a very recent selective sweep in Ugandan mosquitoes which localized to a cluster of cytochrome P450 genes, including some frequently implicated in resistance. Further interrogation revealed a haplotype involving a trio of mutations, a nonsynonymous point mutation in Cyp6p4 (I236M), an upstream insertion of a partial Zanzibar transposable element (TE) and a duplication of the Cyp6aa1 gene. The mutations appear to have originated recently in An. gambiae from the Kenya-Uganda border region around lake Victoria, with stepwise replacement of the double-mutant (Zanzibar TE and Cyp6p4-236M) with the triple-mutant haplotype (including Cyp6aa1 duplication), which has spread into the Democratic Republic of Congo and Tanzania. The triple-mutant haplotype is strongly associated with increased expression of genes able to metabolise pyrethroids; is strongly predictive of resistance to pyrethroids most notably deltamethrin, a commonly-used LLIN insecticide, but importantly, appears less effective against LLINs co-treated with the synergist piperonyl butoxide (PBO). Frequencies of the triple-mutant haplotype remain spatially variable even within countries, suggesting an effective marker system to guide deployment decisions for limited supplies of PBO-pyrethroid co-treated LLINs across African countries. Duplications of Cyp6aa1 gene are common in An. gambiae across Africa and are likely to be a useful diagnostic for high levels of pyrethroid resistance.

Insecticide resistance provides both a pressing threat to the control of vector-borne diseases and insights into the remarkable capacity of natural populations to show rapid evolutionary responses. Malaria control remains heavily dependent on deployment of insecticides, primarily in long lasting insecticide treated nets (LLINs), but resistance in the major malaria vectors has increased over the last 15 years. Identifying genetic mechanisms causing high-level resistance in mosquitoes, which may almost entirely overcome pyrethroid efficacy, is crucial for the development and deployment of potentially resistance-breaking tools. Using the Anopheles gambiae 1000 genomes data we identified a very recent selective sweep in Ugandan mosquitoes which localized to a cluster of cytochrome P450 genes. Further interrogation revealed a haplotype involving a trio of mutations, a point mutation in Cyp6p4, an insertion of a partial Zanzibar transposable element (TE) and a duplication of the Cyp6aa1 gene. The mutations appear to have originated recently in An. gambiae from the Kenya-Uganda border region, with stepwise replacement of the double-mutant (Zanzibar TE and Cyp6p4-236M) with the triple-mutant haplotype (including Cyp6aa1 duplication), which has spread into the Democratic Republic of Congo and Tanzania. The triple-mutant haplotype is strongly associated with increased expression of genes able to metabolise pyrethroids; is strongly predictive of resistance to pyrethroids but importantly, appears less effective against LLINs co-treated with the synergist piperonyl butoxide (PBO). Frequencies of the triple-mutant haplotype remain spatially variable even within countries, suggesting an effective marker system to guide deployment decisions for limited supplies of PBO-pyrethroid co-treated LLINs across African countries.

Wolbachia is a genus of obligate bacterial endosymbionts that infect a diverse range of arthropod species as well as filarial nematodes, with its single described species, Wolbachia pipientis , divided into several ‘supergroups’ based on multilocus sequence typing. Wolbachia strains in mosquitoes have been shown to inhibit the transmission of human pathogens, including Plasmodium malaria parasites and arboviruses. Despite their large host range, Wolbachia strains within the major malaria vectors of the Anopheles gambiae and Anopheles funestus complexes appear at low density, established solely on PCR-based methods. Questions have been raised as to whether this represents a true endosymbiotic relationship. However, recent definitive evidence for two distinct, high-density strains of supergroup B Wolbachia within Anopheles demeilloni and Anopheles moucheti has opened exciting possibilities to explore naturally occurring Wolbachia endosymbionts in Anopheles for biocontrol strategies to block Plasmodium transmission. Here, we utilize genomic analyses to demonstrate that both Wolbachia strains have retained all key metabolic and transport pathways despite their smaller genome size, with this reduction potentially attributable to degenerated prophage regions. Even with this reduction, we confirmed the presence of cytoplasmic incompatibility (CI) factor genes within both strains, with wAnD maintaining intact copies of these genes while the cifB gene was interrupted in wAnM, so functional analysis is required to determine whether wAnM can induce CI. Additionally, phylogenetic analysis indicates that these Wolbachia strains may have been introduced into these two Anopheles species via horizontal transmission events, rather than by ancestral acquisition and subsequent loss events in the Anopheles gambiae species complex. These are the first Wolbachia genomes, to our knowledge, that enable us to study the relationship between natural strain Plasmodium malaria parasites and their anopheline hosts.

Wolbachia is a genus of obligate bacterial endosymbionts that infect a diverse range of arthropod species as well as filarial nematodes, with its single described species, Wolbachia pipientis, divided into several ‘supergroups’ based on multilocus sequence typing. Wolbachia strains in mosquitoes have been shown to inhibit the transmission of human pathogens including Plasmodium malaria parasites and arboviruses. Despite their large host range, Wolbachia strains within the major malaria vectors of the Anopheles (A.) gambiae and A. funestus complexes appear at low density based solely on PCR-based methods. Questions have been raised as to whether this represents a true endosymbiotic relationship. However, recent definitive evidence for two distinct, high-density strains of supergroup B Wolbachia within A. demeilloni and A. moucheti has opened exciting possibilities to explore naturally occurring Wolbachia endosymbionts in Anopheles for biocontrol strategies to block Plasmodium transmission. Here we utilise genomic analyses to demonstrate that both Wolbachia strains have retained all key metabolic and transport pathways despite their smaller genome size. We further confirm the presence of cytoplasmic incompatibility factor genes, despite noticeably few prophage regions. Additionally, phylogenetic analysis indicates that these Wolbachia strains may have been introduced into these two Anopheles species via horizontal transmission events, and unlikely to be by ancestral acquisition and subsequent loss events in the Anopheles gambiae species complex. These are the first Wolbachia genomes that enable us to study the relationship between natural strains Plasmodium malaria parasites and their Anopheline hosts.Impact statementWolbachia naturally infects a wide range of arthropod species, including insect vectors of human pathogens, where they may play a role in inhibiting their replication. These bacteria have been commonly found within Aedes (Ae.) albopictus and Culex pipiens mosquitoes but have been noticeably absent in the Anopheles mosquito genera, which includes all species responsible for malaria transmission. Recent PCR-based methods have suggested the potential for natural Wolbachia strains within the A. gambiae species complex, which includes major malaria vector species including A. gambiae s.s., A. coluzzii and A. arabiensis. We recently reported the presence of stable Wolbachia strains naturally occurring within two different Anopheles species (A. demeilloni and A. moucheti). In this study, we perform comparative genomic analysis of these two Wolbachia genomes against each other and published Wolbachia strains. The current assemblies are some of the smallest sequenced Wolbachia strains of insects, although their metabolic pathway repertoire is comparable to other strains. Interestingly, prophage fragments were identified within only one of the two strains. The findings of this study will be of significant interest to researchers investigating Wolbachia as a potential malaria biocontrol strategy, giving greater insight into the evolution and diversity of this obligate intracellular endosymbiont.Data summarySequence data generated and used for this analysis are available in the National Centre for Biotechnology Information Sequence Read Archive (NCBI SRA bioproject number PRJNA642000). The two assembled Wolbachia genomes are available with genome accession numbers GCA_018491735.2 and GCA_018491625.2. Additional Wolbachia genomes used for comparative analysis are described in the supplementary material.The authors confirm all supporting data, code and protocols have been provided within the article or through supplementary data files. Additional supplementary data files used to generate several figures can be found at: https://figshare.com/projects/Wolbachia_endosymbionts_in_two_Anopheles_species_indicates_independent_acquisitions_and_lack_of_prophage_elements/126533

Preventable diseases still cause huge mortality in low- and middle-income countries. Research in spatial epidemiology and earth observation is helping academics to understand and prioritise how mortality could be reduced and generates spatial data that are used at a global and national level, to inform disease control policy. These data could also inform operational decision making at a more local level, for example to help officials target efforts at a local/regional level. To be usable for local decision-making, data needs to be presented in a way that is relevant to and understandable by local decision makers. We demonstrate an approach and prototype web application to make spatial outputs from disease modelling more useful for local decision making. Key to our approach is: (1) we focus on a handful of important data layers to maintain simplicity; (2) data are summarised at scales relevant to decision making (administrative units); (3) the application has the ability to rank and compare administrative units; (4) open-source code that can be modified and re-used by others, to target specific user-needs. Our prototype application allows visualisation of a handful of key layers from the Malaria Atlas Project. Data can be summarised by administrative unit for any malaria endemic African country, ranked and compared; e.g. to answer questions such as, ‘does the district with the highest malaria prevalence also have the lowest coverage of insecticide treated nets?’. The application is developed in R and the code is open-source. It would be relatively easy for others to change the source code to incorporate different data layers, administrative boundaries or other data visualisations. We suggest such open-source web application development can facilitate the use of data for public health decision making in low resource settings.