Elucidating the genetic basis of metabolic resistance to insecticides in malaria vectors is crucial to prolonging the effectiveness of insecticide-based control tools including long lasting insecticidal nets (LLINs). Here, we show that cis-regulatory variants of the cytochrome P450 gene, CYP6P9b, are associated with pyrethroid resistance in the African malaria vector Anopheles funestus. A DNA-based assay is designed to track this resistance that occurs near fixation in southern Africa but not in West/Central Africa. Applying this assay we demonstrate, using semi-field experimental huts, that CYP6P9b-mediated resistance associates with reduced effectiveness of LLINs. Furthermore, we establish that CYP6P9b combines with another P450, CYP6P9a, to additively exacerbate the reduced efficacy of insecticide-treated nets. Double homozygote resistant mosquitoes (RR/RR) significantly survive exposure to insecticide-treated nets and successfully blood feed more than other genotypes. This study provides tools to track and assess the impact of multi-gene driven metabolic resistance to pyrethroids, helping improve resistance management.

Insecticide resistance in malaria vectors threatens to reverse recent gains in malaria control. Deciphering patterns of gene flow and resistance evolution in malaria vectors is crucial to improve control strategies and preventing malaria resurgence. A genome-wide survey of Anopheles funestus genetic diversity Africa-wide revealed evidences of a major division between southern Africa and elsewhere, associated with different population histories. Three genomic regions exhibited strong signatures of selective sweep, each spanning major resistance loci (CYP6P9a/b, GSTe2 and CYP9K1). However, a sharp regional contrast was observed between populations correlating with gene flow barriers. Signatures of complex molecular evolution of resistance were detected with evidence of copy number variation, transposon insertion and a gene conversion between CYP6P9a/b paralog genes. Temporal analyses of samples before and after bed net scale up suggest that these genomic changes are driven by this control intervention. Multiple independent selective sweeps at the same locus in different parts of Africa suggests that local evolution of resistance in malaria vectors may be a greater threat than trans-regional spread of resistance haplotypes. Malaria control currently relies heavily on insecticide-based vector control interventions. Unfortunately, resistance to insecticides is threatening their continued effectiveness. Metabolic resistance has the greatest operational significance, yet it remains unclear how resistant mosquito populations, evolutionarily respond to the massive selection pressure from control interventions including insecticide-treated nets. Deciphering patterns of gene flow between populations of major malaria vectors such as Anopheles funestus and elucidating genomic signature of resistance evolution are crucial to design resistance management strategies and preventing malaria resurgence. Here, we performed a genome-wide survey of An. funestus genetic diversity from across its continental range using reduced-genome representation (ddRADseq) and whole genome (PoolSeq) approaches revealing evidences of significant barriers to gene flow impacting the spread of insecticide resistance alleles. This study detected signatures of strong selective sweeps occurring in genomic regions controlling cytochrome P450-based and glutathione s-transferase metabolic resistance to insecticides in this species. Fine-scale analysis of the major pyrethroid resistance-associated genomic regions revealed complex molecular evolution with evidence of copy number variation, transposon insertion and a gene conversion highlighting the risk that if this level of selection and spread of resistance continues unabated, our ability to control malaria with current interventions will be compromised.

Galls are plant tissues whose development is induced by another organism for the inducer's benefit. 30,000 arthropod species induce galls, and in most cases the inducing effectors and target plant systems are unknown. Cynipid gall wasps are a speciose monophyletic radiation that induce structurally complex galls on oaks and other plants. We used a model system comprising the gall wasp Biorhiza pallida and the oak Quercus robur to characterise inducer and host plant gene expression at defined stages through the development of galled and ungalled plant tissues, and tested alternative hypotheses for the origin and type of galling effectors and plant metabolic pathways involved. Oak gene expression patterns diverged markedly during development of galled and normal buds. Young galls showed elevated expression of oak genes similar to legume root nodule Nod factor-induced early nodulin (ENOD) genes and developmental parallels with oak buds. In contrast, mature galls showed substantially different patterns of gene expression to mature leaves. While most oak transcripts could be functionally annotated, many gall wasp transcripts of interest were novel. We found no evidence in the gall wasp for involvement of third-party symbionts in gall induction, for effector delivery using virus-like-particles, or for gallwasp expression of genes coding for plant hormones. Many differentially and highly expressed genes in young larvae encoded secretory peptides, which we hypothesise are effector proteins exported to plant tissues. Specifically, we propose that host arabinogalactan proteins and gall wasp chitinases interact in young galls to generate a somatic embryogenesis-like process in oak tissues surrounding the gall wasp larvae. Gall wasp larvae also expressed genes encoding multiple plant cell wall degrading enzymes (PCWDEs). These have functional orthologues in other gall inducing cynipids but not in figitid parasitoid sister groups, suggesting that they may be evolutionary innovations associated with cynipid gall induction. Plant galls are induced by organisms that manipulate host plant development to produce novel structures. The organisms involved range from mutualistic (such as nitrogen fixing bacteria) to parasitic. In the case of parasites, the gall benefits only the gall-inducing partner. A wide range of organisms can induce galls, but the processes involved are understood only for some bacterial and fungal galls. Cynipid gall wasps induce diverse and structurally complex galls, particularly on oaks (Quercus). We used transcriptome and genome sequencing for one gall wasp and its host oak to identify genes active in gall development. On the plant side, when compared to normally developing bud tissues, young gall tissues showed elevated expression of loci similar to those found in nitrogen-fixing root nodules of leguminous plants. On the wasp side, we found no evidence for involvement of viruses or microorganisms carried by the insects in gall induction or delivery of inducing stimuli. We found that gall wasps express many genes whose products may be secreted to the host, including enzymes that degrade plant cell walls. Genome comparisons between galling and non-galling relatives showed cell wall-degrading enzymes are restricted to gall inducers, and hence potentially key components of a gall inducing lifestyle.

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.

Elucidating the complex evolutionary armory that mosquitoes deploy against insecticides is crucial to maintain the effectiveness of insecticide‐based interventions. Here, we deciphered the role of a 6.5‐kb structural variation (SV) in driving cytochrome P450‐mediated pyrethroid resistance in the malaria vector, Anopheles funestus. Whole‐genome pooled sequencing detected an intergenic 6.5‐kb SV between duplicated CYP6P9a/b P450s in pyrethroid‐resistant mosquitoes through a translocation event. Promoter analysis revealed a 17.5‐fold higher activity (p < .0001) for the SV− carrying fragment than the SV− free one. Quantitative real‐time PCR expression profiling of CYP6P9a/b for each SV genotype supported its role as an enhancer because SV+/SV+ homozygote mosquitoes had a significantly greater expression for both genes than heterozygotes SV+/SV− (1.7‐ to 2‐fold) and homozygotes SV−/SV− (4‐to 5‐fold). Designing a PCR assay revealed a strong association between this SV and pyrethroid resistance (SV+/SV+ vs. SV−/SV−; odds ratio [OR] = 2,079.4, p < .001). The 6.5‐kb SV is present at high frequency in southern Africa (80%–100%) but absent in East/Central/West Africa. Experimental hut trials revealed that homozygote SV mosquitoes had a significantly greater chance to survive exposure to pyrethroid‐treated nets (OR 27.7; p < .0001) and to blood feed than susceptible mosquitoes. Furthermore, mosquitoes homozygote‐resistant at the three loci (SV+/CYP6P9a_R/CYP6P9b_R) exhibited a higher resistance level, leading to a far superior ability to survive exposure to nets than those homozygotes susceptible at the three loci, revealing a strong additive effect. This study highlights the important role of structural variations in the development of insecticide resistance in malaria vectors and their detrimental impact on the effectiveness of pyrethroid‐based nets.

Cross-resistance to insecticides in multiple resistant malaria vectors is hampering resistance management. Understanding its underlying molecular basis is critical to implementation of suitable insecticide-based interventions. Here, we established that the tandemly duplicated cytochrome P450s, CYP6P9a/b are driving carbamate and pyrethroid cross-resistance in Southern African populations of the major malaria vector Anopheles funestus. Transcriptome sequencing revealed that cytochrome P450s are the most over-expressed genes in bendiocarb and permethrin-resistant An. funestus. The CYP6P9a and CYP6P9b genes are overexpressed in resistant An. funestus from Southern Africa (Malawi) versus susceptible An. funestus (Fold change (FC) is 53.4 and 17 respectively), while the CYP6P4a and CYP6P4b genes are overexpressed in resistant An. funestus in Ghana, West Africa, (FC is 41.1 and 17.2 respectively). Other up-regulated genes in resistant An. funestus include several additional cytochrome P450s (e.g. CYP9J5, CYP6P2, CYP6P5), glutathione-S transferases, ATP-binding cassette transporters, digestive enzymes, microRNA and transcription factors (FC<7). Targeted enrichment sequencing strongly linked a known major pyrethroid resistance locus (rp1) to carbamate resistance centering around CYP6P9a/b. In bendiocarb resistant An. funestus, this locus exhibits a reduced nucleotide diversity, significant p-values when comparing allele frequencies, and the most non-synonymous substitutions. Recombinant enzyme metabolism assays showed that both CYP6P9a/b metabolize carbamates. Transgenic expression of CYP6P9a/b in Drosophila melanogaster revealed that flies expressing both genes were significantly more resistant to carbamates than controls. Furthermore, a strong correlation was observed between carbamate resistance and CYP6P9a genotypes with homozygote resistant An. funestus (CYP6P9a and the 6.5kb enhancer structural variant) exhibiting a greater ability to withstand bendiocarb/propoxur exposure than homozygote CYP6P9a_susceptible (e.g Odds ratio = 20.8, P<0.0001 for bendiocarb) and heterozygotes (OR = 9.7, P<0.0001). Double homozygote resistant genotype (RR/RR) were even more able to survive than any other genotype combination showing an additive effect. This study highlights the risk that pyrethroid resistance escalation poses to the efficacy of other classes of insecticides. Available metabolic resistance DNA-based diagnostic assays should be used by control programs to monitor cross-resistance between insecticides before implementing new interventions.