RNAseq-based gene expression profiling of the Anopheles funestus pyrethroid resistant strain FUMOZ highlights the predominant role of the duplicated CYP6P9a/b cytochrome P450s

Insecticide-based interventions, notably long-lasting insecticidal nets (LLINs), against mosquito vectors of malaria are currently threatened by pyrethroid resistance. Here, we contrasted RNAseq-based gene expression profiling of laboratory resistant (FUMOZ) and susceptible (FANG) strains of the major malaria vector Anopheles funestus. Cytochrome P450 genes were the predominant over-expressed detoxification genes in FUMOZ, with high expression of the duplicated CYP6P9a (fold-change of 82.23 versus FANG) and CYP6P9b (FC 11.15). Other over-expressed P450s belonged to the same cluster of P450s corresponding to the resistance to pyrethroid 1 (rp1) QTL on chromosome 2R. Several Epsilon class glutathione s-transferases were also over-expressed in FUMOZ, as was the ATP-binding cassette transporter AFUN019220 (ABCA) which also exhibited between-strain alternative splicing events at exon 7. Significant differences in SNP frequencies between strains occurred in resistance QTLs rp1 (CYP6P9a/b, CYP6AA1), rp2 on chromosome 2 L (CYP6Z1, CYP6M7, CYP6Z3) and rp3 on chromosome 3R (CYP9J5, CYP9J4 and CYP9J3). Differences were also detected in CYP4G17 and CYP4G16 genes on the X chromosome, both of which are associated with cuticular resistance in An. gambiae. A close analysis of nonsynonymous diversity at the CYP6P9a/b loci revealed a drastic loss of diversity in FUMOZ with only a single polymorphism and 2 haplotypes vs 18 substitutions and 8 haplotypes in FANG. By contrast, a lowly expressed cytochrome P450 (CYP4C36) did not exhibit diversity differences between strains. We also detected the known pyrethroid resistance conferring amino acid change N384S in CYP6P9b. This study further elucidates the molecular bases of resistance in An. funestus, informing strategies to better manage widespread resistance across Africa.