Global spread and genetic variants of the two CYP9M10 haplotype forms associated with insecticide resistance in Culex quinquefasciatus Say

Potential key genes involved in metabolic resistance to malathion in the southern house mosquito, Culex quinquefasciatus, and functional validation of CYP325BC1 and CYP9M12 as candidate genes using RNA interference

BACKGROUND

Metabolic detoxification is one of the major mechanisms contributing to the development of resistance in mosquitoes, including the southern house mosquito, Culex quinquefasciatus. The three major detoxification supergene families, cytochrome P450s, glutathione S-transferases and general esterases, have been demonstrated to play an important role in metabolic resistance. In this study, we performed differential gene expression analysis based on high-throughput transcriptome sequencing on samples from four experimental groups to give insight into key genes involved in metabolic resistance to malathion in Cx. quinquefasciatus. We conducted a whole transcriptome analysis of field captured wild Cx. quinquefasciatus from Harris County (WI), Texas and a malathion susceptible laboratory-maintained Sebring colony (CO) to investigate metabolic insecticide resistance. Field captured mosquitoes were also phenotypically classified into the malathion resistant and malathion susceptible groups following a mortality response measure conducted using a Centers for Disease Control and Prevention (CDC) bottle assay. The live (MR) and dead (MS) specimens from the bottle assay, along with an unselected WI sample and a CO sample were processed for total RNA extraction and subjected to whole-transcriptome sequencing.

RESULTS

We demonstrated that the genes coding for detoxification enzymes, particularly cytochrome P450s, were highly up-regulated in the MR group compared to the MS group with similar up-regulation observed in the WI group compared to the CO group. A total of 1,438 genes were differentially expressed in comparison between MR and MS group, including 614 up-regulated genes and 824 down-regulated genes. Additionally, 1,871 genes were differentially expressed in comparison between WI and CO group, including 1,083 up-regulated genes and 788 down-regulated genes. Further analysis on differentially expressed genes from three major detoxification supergene families in both comparisons resulted in 16 detoxification genes as candidates potentially associated with metabolic resistance to malathion. Knockdown of CYP325BC1 and CYP9M12 using RNA interference on the laboratory-maintained Sebring strain significantly increased the mortality of Cx. quinquefasciatus after exposure to malathion.

CONCLUSION

We generated substantial transcriptomic evidence on metabolic detoxification of malathion in Cx. quinquefasciatus. We also validated the functional roles of two candidate P450 genes identified through DGE analysis. Our results are the first to demonstrate that knockdown of CYP325BC1 and CYP9M12 both significantly increased malathion susceptibility in Cx. quinquefasciatus, indicating involvement of these two genes in metabolic resistance to malathion.

Multiple cytochrome P450 genes: conferring high levels of permethrin resistance in mosquitoes, Culex quinquefasciatus

Insecticides, especially pyrethroids, are the most important in the insect pest control and preventing insect vector-borne human diseases. However, insect pests, including mosquitoes, have developed resistance in the insecticides that used against them. Cytochrome P450s are associated with insecticide resistance through overexpression and detoxification mechanisms in insect species. In this study, we utilized a powerful tool, the RNAi technique, to determine the roles of key P450 genes overexpressed in permethrin resistant mosquitoes that confer insecticide resistance to unravel the molecular basis of resistance mechanisms in the mosquito Culex quinquefasciatus. The results showed that knockdown of 8 key P450 genes using RNAi techniques significantly decreased resistance to permethrin in resistant mosquitoes. In silico modeling and docking analysis further revealed the potential metabolic function of overexpressed P450 genes in the development of insecticide resistance in mosquitoes. These findings not only highlighted the functional importance of these P450 genes in insecticide resistance, but also revealed that overexpression of multiple P450 genes was responsible for the high levels of insecticide resistance in a mosquito population of Culex quinquefasciatus.

From Global to Local-New Insights into Features of Pyrethroid Detoxification in Vector Mosquitoes

The threat of mosquito-borne diseases continues to be a problem for public health in subtropical and tropical regions of the world; in response, there has been increased use of adulticidal insecticides, such as pyrethroids, in human habitation areas over the last thirty years. As a result, the prevalence of pyrethroid-resistant genetic markers in natural mosquito populations has increased at an alarming rate. This review details recent advances in the understanding of specific mechanisms associated with pyrethroid resistance, with emphasis on features of insecticide detoxification and the interdependence of multiple cellular pathways. Together, these advances add important context to the understanding of the processes that are selected in resistant mosquitoes. Specifically, before pyrethroids bind to their targets on motoneurons, they must first permeate the outer cuticle and diffuse to inner tissues. Resistant mosquitoes have evolved detoxification mechanisms that rely on cytochrome P450s (CYP), esterases, carboxyesterases, and other oxidation/reduction (redox) components to effectively detoxify pyrethroids to nontoxic breakdown products that are then excreted. Enhanced resistance mechanisms have evolved to include alteration of gene copy number, transcriptional and post-transcriptional regulation of gene expression, as well as changes to cellular signaling mechanisms. Here, we outline the variety of ways in which detoxification has been selected in various mosquito populations, as well as key gene categories involved. Pathways associated with potential new genes of interest are proposed. Consideration of multiple cellular pathways could provide opportunities for development of new insecticides.

Novel CONCOMITANT mutations L932F and I936V in the Voltage-Gated Sodium Channel and Its Association With Pyrethroid Resistance in Culex quinquefasciatus (Diptera: Culicidae)

Highly residual pyrethroids such as permethrin have been used for controlling mosquitoes that transmit infectious diseases. However, the selective pressure from such insecticides may result in cross-resistance against other pyrethroids used for household insecticides. In this study, we investigated the susceptibility of Culex quinquefasciatus Say collected from Brazil and Myanmar to permethrin in addition to four types of household pyrethroids. Both strains exhibited high resistance against all pyrethroids tested, indicating cross-resistance. Furthermore, we detected the knockdown resistance (kdr) mutations L932F+I936V in the voltage-gated sodium channel gene (VGSC) in the Brazilian strain. Notably, the L932F+I936V haplotype has previously been observed in in silico data, but it should be detected not directly from living insects. In comparison, a common kdr mutation, L1014F, was detected from the Myanmar strain. Although L1014F was also detected from the Brazilian strain, the allele frequency was too low to affect resistance. Both strains harbored the resistance-associated haplotypes of the cytochrome P450 gene, CYP9M10. The Brazilian strain demonstrated comparable resistance against pyrethroids as that of the Myanmar strain even when a cytochrome P450 inhibitor, piperonyl butoxide was added to the bioassay. Our results suggested that the L932F+I936V mutations confer the Brazilian strain of Cx. Quiquefasciatus with resistance at a comparable level to that conferred by the well-recognized kdr mutation L1014F in the Myanmar strain. The identification of unexplored mutations may improve the diagnosis and understanding of resistance of this medically important species.

Essential Oils of Zingiber Species from Vietnam: Chemical Compositions and Biological Activities

Mosquito-borne diseases are a large problem in Vietnam as elsewhere. Due to environmental concerns regarding the use of synthetic insecticides as well as developing insecticidal resistance, there is a need for environmentally-benign alternative mosquito control agents. In addition, resistance of pathogenic microorganisms to antibiotics is an increasing problem. As part of a program to identify essential oils as alternative larvicidal and antimicrobial agents, the leaf, stem, and rhizome essential oils of several Zingiber species, obtained from wild-growing specimens in northern Vietnam, were acquired by hydrodistillation and investigated using gas chromatography. The mosquito larvicidal activities of the essential oils were assessed against Culex quinquefasciatus, Aedes albopictus, and Ae. aegypti, and for antibacterial activity against a selection of Gram-positive and Gram-negative bacteria, and for activity against Candida albicans. Zingiber essential oils rich in α-pinene and β-pinene showed the best larvicidal activity. Zingiber nudicarpum rhizome essential oil showed excellent antibacterial activity against Enterococcus faecalis, Staphylococcus aureus, and Bacillus cereus, with minimum inhibitory concentrations (MIC) of 2, 8, and 1 μg/mL, respectively. However, the major components, α-pinene and β-pinene, cannot explain the antibacterial activities obtained.

Using targeted next-generation sequencing to characterize genetic differences associated with insecticide resistance in Culex quinquefasciatus populations from the southern U.S

Resistance to insecticides can hamper the control of mosquitoes such as Culex quinquefasciatus, known to vector arboviruses such as West Nile virus and others. The strong selective pressure exerted on a mosquito population by the use of insecticides can result in heritable genetic changes associated with resistance. We sought to characterize genetic differences between insecticide resistant and susceptible Culex quinquefasciatus mosquitoes using targeted DNA sequencing. To that end, we developed a panel of 122 genes known or hypothesized to be involved in insecticide resistance, and used an Ion Torrent PGM sequencer to sequence 125 unrelated individuals from seven populations in the southern U.S. whose resistance phenotypes to permethrin and malathion were known from previous CDC bottle bioassay testing. Data analysis consisted of discovering SNPs (Single Nucleotide Polymorphism) and genes with evidence of copy number variants (CNVs) statistically associated with resistance. Ten of the seventeen genes found to be present in higher copy numbers were experimentally validated with real-time PCR. Of those, six, including the gene with the knock-down resistance (kdr) mutation, showed evidence of a ≥ 1.5 fold increase compared to control DNA. The SNP analysis revealed 228 unique SNPs that had significant p-values for both a Fisher’s Exact Test and the Cochran-Armitage Test for Trend. We calculated the population frequency for each of the 64 nonsynonymous SNPs in this group. Several genes not previously well characterized represent potential candidates for diagnostic assays when further validation is conducted.

The Effect of Permethrin Resistance on Aedes aegypti Transcriptome Following Ingestion of Zika Virus Infected Blood

Aedes aegypti (L.) is the primary vector of many emerging arboviruses. Insecticide resistance among mosquito populations is a consequence of the application of insecticides for mosquito control. We used RNA-sequencing to compare transcriptomes between permethrin resistant and susceptible strains of Florida Ae. aegypti in response to Zika virus infection. A total of 2459 transcripts were expressed at significantly different levels between resistant and susceptible Ae. aegypti. Gene ontology analysis placed these genes into seven categories of biological processes. The 863 transcripts were expressed at significantly different levels between the two mosquito strains (up/down regulated) more than 2-fold. Quantitative real-time PCR analysis was used to validate the Zika-infection response. Our results suggested a highly overexpressed P450, with AAEL014617 and AAEL006798 as potential candidates for the molecular mechanism of permethrin resistance in Ae. aegypti. Our findings indicated that most detoxification enzymes and immune system enzymes altered their gene expression between the two strains of Ae. aegypti in response to Zika virus infection. Understanding the interactions of arboviruses with resistant mosquito vectors at the molecular level allows for the possible development of new approaches in mitigating arbovirus transmission. This information sheds light on Zika-induced changes in insecticide resistant Ae. aegypti with implications for mosquito control strategies.

Genome-wide patterns of copy number variations in Spodoptera litura

Spodoptera litura is a polyphagous pest and can feed on more than 100 species of plants, causing great damage to agricultural production. The SNP results showed that there were gene exchanges between different regions. To explore the variations of larger segments in S. litura genome, we used genome resequencing samples from 14 regions of China, India, and Japan to study the copy number variations (CNVs). We identified 3976 CNV events and 1581 unique copy number variation regions (CNVRs) occupying the 108.5 Mb genome of S. litura. A total of 5527 genes that overlapped with CNVRs were detected. Selection signal analysis identified 19 shared CNVRs and 105 group-specific CNVRs, whose related genes were involved in various biological processes in S. litura. We constructed the first CNVs map in S. litura genome, and our findings will be valuable for understanding the genomic variations and population differences of S. litura.

Molecular identification of four novel cytochrome P450 genes related to the development of resistance of Spodoptera exigua (Lepidoptera: Noctuidae) to chlorantraniliprole

BACKGROUND

The beet armyworm, Spodoptera exigua, is an omnivorous insect that damages a variety of crops worldwide. Chlorantraniliprole is a new diamide insecticide that acts on the ryanodine receptors in insects. The aim of this study was to explore key genes related to the development of chlorantraniliprole resistance in S. exigua.

RESULTS

Transcriptomes were compared between beet armyworms from a susceptible laboratory strain (Sus-Lab) and Sus-Lab screened with LC₂₅ sublethal doses of chlorantraniliprole for six generations (SE-Sel). Ten of 11 cytochrome P450 genes with upregulated expression verified by quantitative reverse transcription polymerase chain reaction (qRT-PCR) in SE-Sel strains were detected in SE-PZ16 and SE-ZY16, two extremely resistant field populations. In addition, expression of four new cytochrome P450 genes, CYP9A21v1, CYP9A21v2, CYP9A21v3 and CYP9A21v4, was shown in the two field populations and was significantly higher in the SE-Sel strain than in the Sus-Lab strain (P < 0.05). Their full-length and protein tertiary structures were also cloned and predicted. The function of CYP9A21v3 was analysed by RNA interference, and the relative expression of CYP9A21v3 in the SE-ZY16 population after feeding on dsRNA was lower than in the control group. Moreover, mortality rates in insects treated at the LC₅₀ of chlorantraniliprole after dsRNA feeding were significantly higher than in the control group 24 h after treatment (P < 0.05).

CONCLUSION

Overexpression of CYP9A21v3 may be a primary factor in the development of chlorantraniliprole resistance in beet armyworms. © 2018 Society of Chemical Industry.

Genomic exaptation enables Lasius niger adaptation to urban environments

Background

The world is rapidly urbanizing, and only a subset of species are able to succeed in stressful city environments. Efficient genome-enabled stress response appears to be a likely prerequisite for urban adaptation. Despite the important role ants play in the ecosytem, only the genomes of ~13 have been sequenced so far. Here, we present the draft genome assembly of the black garden ant Lasius niger – the most successful urban inhabitant of all ants – and we compare it with the genomes of other ant species, including the closely related Camponotus floridanus.

Results

Sequences from 272 M Illumina reads were assembled into 41,406 contigs with total length of 245 MB, and N50 of 16,382 bp, similar to other ant genome assemblies enabling comparative genomic analysis. Remarkably, the predicted proteome of L. niger is significantly enriched relative to other ant genomes in terms of abundance of domains involved in nucleic acid binding, DNA repair, and nucleotidyl transferase activity, reflecting transposable element proliferation and a likely genomic response. With respect to environmental stress, we note a proliferation of various detoxification genes, including glutatione-S-transferases and those in the cytochrome P450 families. Notably, the CYP9 family is highly expanded with 19 complete and 21 nearly complete members - over twice as many compared to other ants. This family exhibits the signatures of strong directional selection, with eleven positively selected positions in ligand-binding pockets of enzymes. Gene family contraction was detected for several components of the olfactory system, accompanied by instances of both directional selection and relaxation.

Conclusions

Our results suggest that the success of L. niger in urbanized areas may be the result of fortuitous coincidence of several factors, including the expansion of the CYP9 cytochrome family due to coevolution with parasitic fungi, the diversification of DNA repair systems as an answer to proliferation of retroelements, and the reduction of olfactory system and behavioral preadaptations from non-territorial subdominant life strategies found in natural environments. Diversification of cytochromes and DNA repair systems along with reduced odorant communication are the basis of L. niger pollutant resistance and polyphagy, while non-territorial and mobilization strategies allows more efficient exploitation of large but patchy food sources.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0867-x) contains supplementary material, which is available to authorized users.

Testing the causality between CYP9M10 and pyrethroid resistance using the TALEN and CRISPR/Cas9 technologies

Recently-emerging genome editing technologies have enabled targeted gene knockout experiments even in non-model insect species. For studies on insecticide resistance, genome editing technologies offer some advantages over the conventional reverse genetic technique, RNA interference, for testing causal relationships between genes of detoxifying enzymes and resistance phenotypes. There were relatively abundant evidences indicating that the overexpression of a cytochrome P450 gene CYP9M10 confers strong pyrethroid resistance in larvae of the southern house mosquito Culex quinquefasciatus. However, reverse genetic verification has not yet been obtained because of the technical difficulty of microinjection into larvae. Here, we tested two genome editing technologies, transcription activator-like effector nucleases (TALEN)s and clustered regularly interspaced short palindromic repeats (CRISPR/Cas9), to disrupt CYP9M10 in a resistant strain of C. quinquefasciatus. Additionally, we developed a novel, effective approach to construct a TALE using the chemical cleavage of phosphorothioate inter-nucleotide linkages in the level 1 assembly. Both TALEN and CRISPR/Cas9 induced frame-shifting mutations in one or all copies of CYP9M10 in a pyrethroid-resistant strain. A line fixed with a completely disrupted CYP9M10 haplotype showed more than 100-fold reduction in pyrethroid resistance in the larval stage.

Genetic characterization of field-evolved resistance to phosphine in the rusty grain beetle, Cryptolestes ferrugineus (Laemophloeidae: Coleoptera)

Inheritance of resistance to phosphine fumigant was investigated in three field-collected strains of rusty grain beetle, Cryptolestes ferrugineus, Susceptible (S-strain), Weakly Resistant (Weak-R) and Strongly Resistant (Strong-R). The strains were purified for susceptibility, weak resistance and strong resistance to phosphine, respectively, to ensure homozygosity of resistance genotype. Crosses were established between S-strain×Weak-R, S-strain×Strong-R and Weak-R×Strong-R, and the dose mortality responses to phosphine of these strains and their F1, F2 and F1-backcross progeny were obtained. The fumigations were undertaken at 25°C and 55% RH for 72h. Weak-R and Strong-R showed resistance factors of 6.3× and 505× compared with S-strain at the LC50. Both weak and strong resistances were expressed as incompletely recessive with degrees of dominance of -0.48 and -0.43 at the LC50, respectively. Responses of F2 and F1-backcross progeny indicated the existence of one major gene in Weak-R, and at least two major genes in Strong-R, one of which was allelic with the major factor in Weak-R. Phenotypic variance analyses also estimated that the number of independently segregating genes conferring weak resistance was 1 (nE=0.89) whereas there were two genes controlling strong resistance (nE=1.2). The second gene, unique to Strong-R, interacted synergistically with the first gene to confer a very high level of resistance (~80×). Neither of the two major resistance genes was sex linked. Despite the similarity of the genetics of resistance to that previously observed in other pest species, a significant proportion (~15 to 30%) of F1 individuals survived at phosphine concentrations higher than predicted. Thus it is likely that additional dominant heritable factors, present in some individuals in the population, also influenced the resistance phenotype. Our results will help in understanding the process of selection for phosphine resistance in the field which will inform resistance management strategies. In addition, this information will provide a basis for the identification of the resistance genes.

A single nucleotide change in a core promoter is involved in the progressive overexpression of the duplicated CYP9M10 haplotype lineage in Culex quinquefasciatus

Although the importance of cis-acting mutations on detoxification enzyme genes for insecticide resistance is widely accepted, only a few of them have been determined as concrete mutations present in genomic DNA till date. The overexpression of a cytochrome P450 gene, CYP9M10, is associated with pyrethroid resistance in the southern house mosquito Culex quinquefasciatus. The haplotypes of CYP9M10 exhibiting overexpression (resistant haplotypes) belong to one specific phylogenetic lineage that shares high nucleotide sequence homology and the same insertion of a transposable element. Among the resistant haplotypes, allelic progression involving an additional cis-acting mutation and gene duplication evolved a CYP9M10 haplotype associated with extremely high transcription and strong pyrethroid resistance. Here we show that a single nucleotide substitution G-27A, which is located near the transcription start site of CYP9M10, is involved in the progression of the duplicated haplotype lineage. The deletion of a 7-bp AT-rich sequence that includes nucleotide -27 inhibited the initiation of transcription from the original transcriptional initiation site. The mutation was suspected to reside within a core promoter, TATA-box, of CYP9M10.

Designing drugs that overcome antibacterial resistance: where do we stand and what should we do?

INTRODUCTION

In recent years, infections caused by multidrug-resistant bacterial pathogens have become a huge issue to public healthcare systems. Indeed, the misuse of antibiotics has led to, over the past 30 years, the emergence of a number of resistant bacterial strains including Staphylococcus aureus, Neisseria gonorrhoeae, Escherichia coli and Mycobacterium tuberculosis. Unfortunately, efforts to produce new antibiotics have not been sufficient to cope with the emergence of these new antibiotic-resistant (AR) strains.

AREAS COVERED

There is an urgent need to invent and employ unconventional strategies for antimicrobial drug development to tackle the rising global threats imposed by the spread of antimicrobial resistance. Herein, the authors discuss these novel design strategies and provide their expert perspective on the subject.

EXPERT OPINION

To deal with the growing threat of AR, it is important to cut down the use of antibiotics to the very minimum to diminish the risk of unknown drug-resistant bacteria and increase antibacterial vaccination programs. Furthermore, it is important to develop new classes of antibiotics that can deal with multidrug-resistant bacterial pathogens.

Pyrethroid resistance in Culex pipiens mosquitoes

Mosquitoes within the Culex pipiens complex are widely distributed and important in the transmission of many human diseases. Insecticides, pyrethroids in particular, remain a mainstay for control of these important vectors. In this paper we review what is known about the levels, mechanisms and fitness costs of pyrethroid resistance in Cx. pipiens. Pyrethroid resistance in Cx. pipiens is a global problem, and resistance ratios of up to 7000-fold have been found in larvae of field collected mosquitoes. However, there is considerable variation between populations, indicating significant geographic heterogeneity of the resistance. The two major mechanisms of resistance to pyrethroids in Culex are mutations in Vssc (target site insensitivity) and overexpression of cytochrome P450(s) (increased detoxification). The most frequently reported Vssc mutation is L1014F (i.e. kdr), which has been found throughout the world. The L1014S mutation has been found in Cx. p. pallens from Japan and China, and in Cx. p. pipiens from China. The L1014C mutation has only been reported for Cx. p. pipens molestus from China and the V1016G mutation has only been reported from Saudi Arabia. Studies on the P450s of Cx. pipiens have identified several that are overexpressed (measured as transcript levels) in pyrethroid resistant strains. CYP9M10 is consistently overexpressed in pyrethroid resistant Cx. pipiens from at least seven countries, suggesting this P450 might be of global importance in resistance. Both CYP9M10-mediated pyrethroid resistance and kdr have fitness costs in the absence of insecticides under certain environmental conditions. Research needs and future directions are discussed.

Permethrin resistance variation and susceptible reference line isolation in a field population of the mosquito, Culex quinquefasciatus (Diptera: Culicidae)

This study examines the genetic variations and mechanisms involved in the development of permethrin resistance in individual mosquitoes from a field population of Culex quinquefasciatus, HAmCq(G0) , and characterizes susceptible reference lines of mosquitoes with a similar genetic background to the field HAmCq(G0) strain. Six upregulated cytochrome P450 genes, CYP9M10, CYP9J34, CYP6P14, CYP9J40, CYP6AA7, and CYP4C52v1, previously identified as being upregulated in the larvae of resistant HAmCq(G8) mosquitoes were examined in the larvae of 3 strains (susceptible S-Lab, parental HAmCq(G0) and permethrin-selected highly resistant HAmCq(G8) ) and 8 HAmCq(G0) single-egg raft colonies, covering a range of levels of susceptibility/resistance to permethrin and exhibiting different variations in the expression of A and/or T alleles at the L-to-F kdr locus of the sodium channel. The 2 lines with the lowest tolerance to permethrin and bearing solely the susceptible A allele at the L-to-F kdr locus of the sodium channels, from colonies Cx_SERC5 and Cx_SERC8, showed lower or similar levels of all 6 of the P450 genes tested compared with the S-Lab strain, suggesting that these 2 lines could be used as the reference mosquitoes in future studies characterizing insecticide resistance in HAmCq mosquitoes. This study also provides a detailed investigation of the mechanisms involved in insecticide resistance in individuals within a population: individuals with elevated levels of resistance to permethrin all displayed one or more potential resistance mechanisms-either elevated levels of P450 gene expression, or L-to-F mutations in the sodium channel, or both.

Mechanisms of pyrethroid resistance in the dengue mosquito vector, Aedes aegypti: target site insensitivity, penetration, and metabolism

Aedes aegypti is the major vector of yellow and dengue fevers. After 10 generations of adult selection, an A. aegypti strain (SP) developed 1650-fold resistance to permethrin, which is one of the most widely used pyrethroid insecticides for mosquito control. SP larvae also developed 8790-fold resistance following selection of the adults. Prior to the selections, the frequencies of V1016G and F1534C mutations in domains II and III, respectively, of voltage-sensitive sodium channel (Vssc, the target site of pyrethroid insecticide) were 0.44 and 0.56, respectively. In contrast, only G1016 alleles were present after two permethrin selections, indicating that G1016 can more contribute to the insensitivity of Vssc than C1534. In vivo metabolism studies showed that the SP strain excreted permethrin metabolites more rapidly than a susceptible SMK strain. Pretreatment with piperonyl butoxide caused strong inhibition of excretion of permethrin metabolites, suggesting that cytochrome P450 monooxygenases (P450s) play an important role in resistance development. In vitro metabolism studies also indicated an association of P450s with resistance. Microarray analysis showed that multiple P450 genes were over expressed during the larval and adult stages in the SP strain. Following quantitative real time PCR, we focused on two P450 isoforms, CYP9M6 and CYP6BB2. Transcription levels of these P450s were well correlated with the rate of permethrin excretion and they were certainly capable of detoxifying permethrin to 4'-HO-permethrin. Over expression of CYP9M6 was partially due to gene amplification. There was no significant difference in the rate of permethrin reduction from cuticle between SP and SMK strains.