Next Article in Journal
Methyl Donor Micronutrients: A Potential Dietary Epigenetic Target in Systemic Lupus Erythematosus Patients
Previous Article in Journal
Study of Effector CD8+ T Cell Interactions with Cortical Neurons in Response to Inflammation in Mouse Brain Slices and Neuronal Cultures
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 24
Issue 4
10.3390/ijms24043170
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multiple-P450 Gene Co-Up-Regulation in the Development of Permethrin Resistance in the House Fly, Musca domestica

by
Ming Li
1,2,
Xuechun Feng
1,3,
William R. Reid
1,4,
Fang Tang
1,5 and
Nannan Liu
1,*
1
Department of Entomology and Plant Pathology, Auburn University, Auburn, AL 36849, USA
2
Department of Entomology, University of California, San Diego, CA 92093, USA
3
Institute of Infectious Diseases, Shenzhen Bay Laboratory, Shenzhen 518000, China
4
Department of Biochemistry, State University of New York at Buffalo, Buffalo, NY 14203, USA
5
College of Forestry, Nanjing Forestry University, Nanjing 210037, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(4), 3170; https://doi.org/10.3390/ijms24043170
Submission received: 5 December 2022 / Revised: 17 January 2023 / Accepted: 24 January 2023 / Published: 6 February 2023
(This article belongs to the Section Molecular Toxicology)
Browse Figures
Review Reports Versions Notes

Abstract

:
This paper reports a study conducted at the whole transcriptome level to characterize the P450 genes involved in the development of pyrethroid resistance, utilizing expression profile analyses of 86 cytochrome P450 genes in house fly strains with different levels of resistance to pyrethroids/permethrin. Interactions among the up-regulated P450 genes and possible regulatory factors in different autosomes were examined in house fly lines with different combinations of autosomes from a resistant house fly strain, ALHF. Eleven P450 genes that were significantly up-regulated, with levels > 2-fold those in the resistant ALHF house flies, were in CYP families 4 and 6 and located on autosomes 1, 3 and 5. The expression of these P450 genes was regulated by trans- and/or cis-acting factors, especially on autosomes 1 and 2. An in vivo functional study indicated that the up-regulated P450 genes also conferred permethrin resistance in Drosophila melanogaster transgenic lines. An in vitro functional study confirmed that the up-regulated P450 genes are able to metabolize not only cis- and trans-permethrin, but also two metabolites of permethrin, PBalc and PBald. In silico homology modeling and the molecular docking methodology further support the metabolic capacity of these P450s for permethrin and substrates. Taken together, the findings of this study highlight the important function of multi-up-regulated P450 genes in the development of insecticide resistance in house flies.

Graphical Abstract

1. Introduction

Efforts to characterize the molecular mechanisms involved in the development of insecticide resistance in insect pests have generally focused on building a better understanding of the mechanisms involved, providing vital information for practical applications such as the design of novel strategies to prevent or minimize the spread and evolution of resistance development and the control of insect pests [1]. The interactions among multiple mechanisms or genes in response to the development of insecticide resistance have been extensively explored; moreover, the transcriptional up-regulation of the detoxification machinery responsible for increasing the metabolism of insecticides into less harmful substances and facilitating insecticide excretion is known to play an important role in allowing insects to defend themselves against insecticides [2,3]. Cytochrome P450s are known to be important components in the detoxification machinery of insects; their transcriptional up-regulation in insects increases both protein production and enzymatic activities, which, in turn, enhance the metabolic detoxification of insecticides and plant toxins, thus leading to the development of insecticide resistance [1,4]. Multiple P450s have also been shown to work in combination to boost insecticide resistance through co-up-regulation in resistant insects [1,5,6,7,8,9,10]. Recent advances in genome/whole transcriptome sequencing technology represent a major milestone, opening up new ways for researchers to explore precisely how many P450 genes are involved in insecticide resistance in individual insect pests [11,12,13,14,15,16,17].
The house fly, Musca domestica, is a major domestic, medical and veterinary pest that can transmit over 200 human and animal pathogens [18]. Control efforts have mainly relied on the application of insecticides, primarily pyrethroids [19], but the major barrier hampering house fly control, as in other insect species, is their remarkable ability to develop not only resistance to the insecticide used against them but also cross-resistance to unrelated classes of insecticides [4,12]. The house fly is widely used as a useful model with which to study and predict insecticide resistance because of its ability to rapidly develop resistance and cross-resistance to insecticides, its well-established linkage map for five autosomes and two sex chromosomes (X and Y), its relatively well-understood biochemistry and genetics of insecticide resistance, and its readily available genomic and transcriptome data analyses [11,12]. Building on our earlier adult house fly transcriptome data analysis [11], where cytochrome P450s were found to be the major functional group up-regulated in resistant house flies, the current study moved on to characterize the total number of P450 genes directly involved in resistance in a single resistant insect. A systematic study conducted at the whole transcriptome level was used to develop an in-depth understanding of these P450 genes, with a focus in how their co-up-regulation, interaction, metabolism and function contribute to the development of insecticide resistance in house flies.

2. Results

2.1. Cytochrome P450 Genes and Their Transcriptome Expression Profiles in M. domestica

There are 146 P450 genes identified in the genome database of house flies [12], of which 86 are expressed in ALHF adult house flies based on our previous whole transcriptome analysis [11]. These cytochrome P450s were selected to pinpoint the most important P450s involved in the development of permethrin resistance for this study. The majority of these 86 P450 genes were assembled into clans 3 and 4, with thirty-nine P450 genes in clan 3, of which twenty-seven were in family 6, six were in family 9, three were in family 28, one was in family 310 and one was in family 317 (Table 1 and Table S2). Of the thirty P450 genes in clan 4, twenty-three were in family 4, two were in family 313, one was in family 311 and one was in family 318. Of the remaining P450 genes, four were in families 304, 305, 306 and 18 of clan 2, and thirteen P450 genes were found to be in five families in the mitochondrial clan, namely 12, 301, 302, 314 and 315. Interestingly, the detectable P450 genes in the ALHF house flies showed a clear expansion in families 4, 6 and 9 compared to the other families. P450 genes in these three families have previously been implicated in environmental response, xenobiotic metabolism and the development of resistance in insects [20,21,22,23].
The relative expression profile of the 86 P450 genes was examined in resistant ALHF and two susceptible aabys and CS house fly strains to eliminate false positive results. The expression levels of 11 P450 genes were found to be significantly (p < 0.05) up-regulated, with levels > 2-fold higher in the ALHF house flies compared to the levels in both the aabys and CS house flies (Table 2 and Table S3). These genes were distributed into two clans (clans 3 and 4), with five genes in family 4 and six in family 6. Fourteen P450 genes were found to be significantly down-regulated in the ALHF strain compared to both the aabys and CS strains (Table S3). The remaining 59 P450 genes showed no significant difference in expression levels between ALHF and either or both the aabys and CS strains.

2.2. Autosome Location and Co-Up-Regulation of P450 Gene Expression in Insecticide-Resistant House Flies

Autosomal location analyses were conducted for all 11 of the up-regulated P450 genes using allele-specific PCR (AS-PCR) determination utilizing five house fly BC1 lines, and allele-specific primer pairs designed based on the specific sequence of the genes from ALHF, by placing a specific nucleotide polymorphism at the 3′ end of one primer of each primer pair to facilitate the preferential amplification of specific alleles from ALHF [24]. Our results show that the ALHF allele-specific primer sets for CYP6D3 and CYP6D10 amplified specific DNA fragments only in flies containing the autosome 1 wild-type marker from ALHF (Figure 1), confirming that these two P450 genes are indeed on autosome 1. Similarly, CYP4G13, CYP4G99 and CYP4S24 were found to be located on autosome 3, and CYP4E10, CYP4E11, CYP6A36, CYP6A40, CYP6A52 and CYP6A58 on autosome 5 (Figure 1).
To better understand the cis-/trans-regulation of the up-regulated P450 genes in ALHF house flies, we next determined the effects and co-regulation of various factors on the expression of the up-regulated P450 genes by analyzing the gene expression changes resulting from autosome replacement in ALHF. This enabled us to evaluate the role that the factors on each autosome plays in P450 gene overexpression in ALHF. The P450 gene expression levels were examined in resistant ALHF, susceptible aabys and five house fly homozygous lines: A2345, A1345, A1245, A1235 and A1234; these lines represent the ALHF strains where autosomes 1, 2, 3, 4 and 5, respectively, have been replaced by the autosome from aabys. We found no significant change in the level of expression in any of the P450 genes tested when autosome 4 of ALHF was replaced with that from aabys (i.e., line A1235; Figure 2, Table 2), strongly suggesting that the factors/genes on autosome 4 do not play a major role in the up-regulation of P450 genes in ALHF. This agrees with the findings reported in previous studies by Tian et al. [25], Liu and Scott [26], and Liu and Yue [27].
Clear changes in the P450 gene expression levels were observed when any one of autosomes 1, 2, 3 and 5 in ALHF was replaced by its aabys counterpart. Apart from one P450 gene, CYP6A52, whose up-regulation was solely controlled by cis factor(s) on autosome 5, where the P450 gene itself is located, the expression of the remaining 10 up-regulated P450 genes were all linked to factors on more than one autosome (Figure 2, Table 2). These results suggest that cis/trans factors are capable of co-regulating the P450 gene expression responsible for insecticide resistance. Of the five P450 genes (CYP4E10, CYP4E11, CYP6A36, CYP6A40 and CYP6A58) located on autosome 5, the expression of three (CYP4E10, CYP6A36, and CYP6A58) was controlled by trans-regulatory factors on autosome 1 and 2, while that of the other two (CYP4E11 and CYP6A40) was controlled by trans-regulatory factors on autosome 2 only. Of the three genes located on autosome 3, two (CYP4G13 and CYP4S24) were found to be regulated by the trans-regulatory factors on autosomes 1 and 2, and the other (CYP4G99) by trans-regulatory factors on autosome 1. One of the two genes located on autosome 1 (CYP6D10,) was up-regulated by trans-regulatory factors on autosomes 2 and 5, while the other (CYP6D3) was regulated by trans-regulatory factors on autosome 2 (Figure 2, Table 2). Clearly, factors on autosomes 2 and 5, especially those on autosome 2, are commonly involved in the up-regulation of P450 gene expression associated with the development of resistance in house flies.

2.3. Up-Regulated P450 Genes in Transgenic D. melanogaster—An Alternative/Independent Approach to Defining P450 Gene Functions

To further validate the function of the up-regulated P450 genes in the development of insecticide resistance in ALHF house flies, three P450 genes, CYP4S24, CYP6A36, and CYP6D10, were selected for a transgenic study using the GAL4-UAS enhancer trap system of D. melanogaster. These three genes represent the up-regulated P450 genes on each of three different autosomes: the expression of CYP4S24, located on autosome 3, is regulated by trans-regulatory factors on autosomes 1 and 2; the expression of CYP6A36, located on autosome 5, is regulated by trans-regulatory factors on autosomes 1 and 2; and the expression of CYP6D10, located on autosome 1, is regulated by trans-regulatory factors on autosomes 2 and 5.
The first step was to confirm the presence of house fly P450 genes in the transgenic lines of D. melanogaster via RT-PCR; the results show that all the transgenic lines did indeed express the target transgenes (Figure 3A). To test whether the GAL4-UAS expression system in D. melanogaster increases the expression of target house fly P450 genes, qRT-PCR was further employed to detect the expression levels among crosses of the transgenic lines including: a control (Bloomington stock #24484 D. melanogaster line containing an empty pUAST vector), GAL4 (the ubiquitous Act5C driver line), the control crossed with GAL4 (the F1 progeny produced by a cross between the control females and the GAL4 males), the two homozygous transgenic P450 lines, and the P450 transgenic lines crossed with the GAL4 driver line (the F1 progeny from a cross between the P450 homozygous transgene line females and the GAL4 males). There was no sign of the expression of any of the three P450 genes in the non-transgenic D. melanogaster lines (control, GAL4 and control × GAL4), but it was detected in both of the transgenic lines (P450 and P450 × GAL4). The expression of these P450 genes was enhanced by the GAL4-UAS expression system of D. melanogaster (P450 × GAL4 > P450) (Figure 3B).
The next step was to characterize the sensitivity of the non-transgenic and transgenic D. melanogaster lines to permethrin. The bioassay results showed no significant difference in permethrin toxicity among the three non-transgenic lines based on overlapping 95% confidence intervals. However, transgenic D. melanogaster lines with house fly P450s showed significantly increased levels of resistance to permethrin compared to non-transgenic lines (Figure 3C). Of the three P450s tested, transgenic flies with the house fly CYP6A36 gene had the highest level of resistance to permethrin (a ~6.5-fold increase compared with the control lines), followed by flies with CYP6D10 (a ~4-fold increase) and CYP4S24 (a ~2-fold increase). These results indicate that these up-regulated P450 genes in the house fly boosted permethrin resistance in D. melanogaster, highlighting the important function of these P450 genes in the permethrin resistance of M. domestica.

2.4. Functional Expression of P450s and In Vitro Metabolism Studies

Functional validation to confirm that metabolic enzymes encoded by the overexpressed metabolic genes do indeed metabolize the insecticides in resistant mosquitoes is crucial in efforts to pin down the gene’s precise function in the metabolism of insecticides and in the development of insecticide resistance. We therefore heterogeneously expressed CYP4S24, CYP6A36 and CYP6D10 proteins with house fly cytochrome P450 reductase (CPR) in a baculovirus expression system. Each of the P450s was co-expressed with CPR via simultaneous infection of Sf9 cells through two recombinant viruses (P450-recombinant baculovirus (P450rbv) with an MOI of 0.5 and CPR-recombinant baculovirus (CPRrbv) with an MOI of 0.05) in a standard suspension culture [28]. The functional catalytic activity of these recombinant P450/CPR proteins was confirmed via an ECOD assay, where 7-ethoxycoumarin was hydrolyzed to 7-hydroxycoumarin by recombinant P450 and CPR proteins with different rates of metabolism (Table S4). CYP6D10/CPR exhibited the highest ECOD activity, followed by CYP6A36, and then, CYP4S24, confirming that all three of these recombinant P450/CPR proteins can indeed convert significant amounts of the substrate 7-ethoxycoumarin to the fluorescent product 7-hydroxycoumarin. No significant P450 activity was observed in proteins isolated from the original Sf9 cells.
Permethrin metabolic activity was assayed using microsomal proteins of P450/CPR in the presence or absence of NADPH. The degradation of the substrate was monitored via reversed-phase HPLC. Our in vitro study demonstrated that all three P450s and CPR co-expressed microsomal proteins and metabolized both trans- and cis- permethrin. Significant levels of the permethrin were metabolized by the three P450s and CPR co-expressed microsomal proteins in a 20 µM permethrin solution containing both cis- and trans-isomers compared to a no-NADPH control after a 120 min incubation period. Of the three P450s, CYP6A36 showed the highest activity toward both trans- and cis-permethrin (64.27% and 59.81% of the trans- and cis-permethrin, respectively, was eliminated in 120 min), closely followed by CYP6D10 (50.55% and 43.01% of the trans- and cis-permethrin, respectively). Meanwhile, CYP4S24 exhibited the lowest activity toward trans- and cis-permethrin (37.53% and 30.61% of the trans- and cis-permethrin, respectively) (Figure 4). These results clearly demonstrate that all three up-regulated P450s are capable of metabolizing permethrin in vitro.
It has been reported that pyrethroids can be hydrolyzed by carboxylesterases to produce PBalc and PBald, and these metabolites can be further processed to 4′-OH PBalc and PBald by P450s in mammals and mosquitos [29,30,31,32,33,34,35]. We therefore conducted a metabolism assay of PBalc and PBald using house fly P450s and CPR co-expressed microsomal proteins. All three P450s metabolized both PBalc and PBald, with CYP4S24 exhibiting the highest ability to metabolize both PBalc and PBald (61.20% and 92.97% of the PBalc and PBald, respectively, was eliminated in 120 min) (Figure 4). Our results confirm that CYP6A36, CYP6D10 and CYP4S24 are not only able to metabolize cis- and trans-permethrin, but can also metabolize permethrin’s metabolites, PBalc and PBald.

2.5. Homology Modeling and Substrate Docking

To understand the underlying metabolic mechanisms and the differences between the responses of these selected P450s toward permethrin, PBalc and PBald, homology modeling and molecular docking simulations were conducted. Not surprisingly, these P450s were found to contain a number of conserved P450 characteristics, including the oxygen binding motif of helix I (A/GGxE/DTT/S), the heme-binding signature (FXXGXRXCXG) near the C-terminal, the EXXR motif which stabilizes the heme core, and six predicted substrate recognition sites (SRSs) (Figures S1A–C and S2).
The active site for P450 is buried deep within the enzyme structure. It is connected to the surrounding environment by a network of channels that serve as access/egress paths and may determine the substrate specificity of P450s [36,37]. The most probable candidate substrate access/egress channels (solvent and family 2) were therefore investigated in this study. An analysis of the geometry involved showed that eight channels (2a, 2ac, 2b, 2c, 2d, 2e, 2f and S) were present in CYP6A36 and CYP6D10 (Figure 5). These channels merge together to form a large opening to the heme prosthetic group, and thus, provide a large active cavity volume. The seven channels (2a, 2ac, 2b, 2c, 2d, 2f and S) detected in CYP4S24 formed a much more restricted narrow opening to the heme prosthetic group and a small active cavity, likely due to two amino acids (Glutamate 313 and Threonine 317) that protrude into the helix I structure (Figure 5).
A comparison of permethrin, PBalc and PBald binding conformations at the active sites of the P450 models revealed several possible binding and metabolic mechanisms of P450s to these chemicals. The predicted active site of these three P450s easily accommodates both cis- and trans-permethrin. Although permethrin can dock with the gem-dimethyl group, with 5-phenoxybenzyl carbon and 4′-phenoxybenzyl carbon closer to the heme center of the active sites of CYP6A36 and CYP6D10, docking simulations suggest that 4′-hydroxylation is a major route for metabolism (Figure 6A–D, Table S5), which agrees with the findings reported in other studies [29,38,39,40,41]. Only one binding conformation has the gem-dimethyl groups of the cyclopropane moiety approaching the heme, which, when modeled, produces cis- or trans-hydroxymethyl-permethrin at the active site of CYP4s24 (Figure 6E,F, Table S5). This has previously been identified as a minor route of pyrethroid metabolism [38].
The high binding affinity and the short distance from the substrate metabolic sites to the P450 heme iron are responsible for the high electrophilic character of this position [42] and the high catalytic reactivity [40]. Interestingly, our docking study not only showed that the catalytic reactivity of these three P450s is potentially greater with trans-permethrin compared to cis-permethrin, as trans-permethrin docked closer to the heme center and presented a higher binding affinity, but also demonstrated that CYP6A36 has the strongest ability to metabolize permethrin compared with CYP4S24, which had the lowest ability to metabolize permethrin (Table S5). These findings are consistent with those of our transgenic D. melanogaster study and in vitro permethrin metabolism study.
In CYP6A36 and CYP6D10, PBalc and PBald docked with productive poses such that the 4′-phenoxybenzyl carbon was oriented above the heme (Figure 7A,B,D,E, Table S5), creating two products (4′-hydroxylated PBalc (4′-OH-PBalc) [34] and 4′-hydroxylated PBald (4′-OH-PBald), respectively), as predicted. However, in CYP4S24, PBalc and PBald docked with the carbon atom in the hydroxyl and the aldehyde group directly over the heme, with high binding affinity and only a short distance from the heme iron (Figure 7C,F, Table S5). This positioning may lead to the production of 3-phenoxybenzoic acid (PBAcid), which has been identified as the final metabolite of pyrethroids [29,31,34,35].

3. Discussion

Previous studies have suggested that the constitutively increased expression and induction of P450 genes in insects is responsible for increased levels of total P450 and P450 activities, which significantly affect the disposition of xenobiotics or endogenous compounds in the tissues of organisms, and thus, alter their toxicological effects. Multiple P450s are commonly thought to share roles associated with the development of insecticide resistance through co-regulation in resistant insects [7,8,24,43,44,45,46]. Using coupled CRISPR-Cas9-based genome editing techniques and in vitro metabolism analysis, Wang et al. [47] approved a cluster of nine P450 genes involved in the metabolism of xenobiotics in cotton bollworms, Helicoverpa armigera. While many studies have revealed that the overexpression of P450 genes is regulated by unknown trans- or cis-regulatory factors [1,11,26,27,43,44,48,49,50,51,52,53,54], mutations in P450 proteins have also been proven to impact the xanthotoxin metabolism of P450s [55]. Taken together, these studies have demonstrated not only that insecticide resistance in a single insect is conferred by the co-contribution of multiple P450 genes, but also that the interaction with regulatory factors plays an important role. In view of the above, the current paper presents a systematic study conducted at the whole transcriptome level to characterize the P450 genes known to be involved in the development of pyrethroid resistance, by analyzing the expression profiles of a total of 86 cytochrome P450 genes [11] in house fly strains with different levels of resistance to pyrethroid/permethrin. A total of 11 P450 genes that were found to be significantly up-regulated with levels > 2-fold those in resistant ALHF house flies were in families 4 and 6, located on autosomes 1, 3 and 5. Except for CYP6A52, whose up-regulation was solely controlled by the cis factor(s) on a single autosome, the up-regulated P450 genes were shown to be up-regulated through interactions with factors on more than one autosome (primarily autosomes 1 and 2) suggesting that the expression of up-regulated P450 genes is controlled by the cis or trans factors in resistant house flies. Hence, further characterizing the regulatory factors that are involved in resistance P450 gene expression will not only support efforts to elucidate the mechanisms of insecticide resistance development, but also facilitate work being carried out to diagram the entire transcriptional regulatory network and select new, more efficient targets for insect pest control.
While our knowledge regarding the regulation of insect P450 gene expression in insecticide resistance and the identification of regulatory factors in insect P450 expression is still far from complete for most insect P450s, several recent studies have clearly demonstrated the ambiguity of the transcriptional regulation mechanisms involved in the regulation of resistance P450 gene expression through cis elements and trans-regulatory factors [2,56,57,58,59,60,61,62,63]. A few of the regulatory factors known to be involved in insect P450 gene expression in insecticide resistance, including orphan nuclear receptors, cap ‘n’ collar C (CncC)/muscle aponeurosis fibromatosis (Maf), G-protein-coupled receptors (GPCRs), and the cAMP-response element binding protein, have been identified in several insect species [54,57,59,60,64]. Hu et al. [65] have reported the constitutive up-regulation of transcription factors CncC/Maf to be at least partially responsible for the upregulated expression of CYP321A8 in a resistant strain of Spodoptera exigua and demonstrated that CncC/Maf enhanced the expression of CYP321A8 by binding to specific sites in the promoter. The authors also found that cis-regulatory elements resulting from a mutation in the CYP321A8 promoter in a resistant strain facilitate the binding of the orphan nuclear receptor, Knirps, and enhance promoter activity. Yang et al. [64] discovered that the overexpression of CYP6CM1 in white flies, known to be responsible for the development of resistance to neonicotinoid insecticides, is regulated via the trans-regulatory factor cAMP-responsive element modulator (CREM), the activation of which is mediated by the mitogen-activated protein kinase (MAPK) signaling pathway. A GPCR regulatory pathway through Gαs, AC, PKA and cAMP in mosquitoes that has been implicated in the regulation of multiple resistance P450 gene expression levels, eventually leading to the enhanced detoxification of insecticides in the mosquito Culex quinquefasciatus, has been systematically studied and modeled [2,54,58,59]. The importance of GPCR genes and GPCR signaling pathways in the regulation of resistance P450 genes has also been reported in other resistant insect species, including Lymantria dispar [66], Culex pipiens [67] and M. domestica [68]. Li et al. [2,54] conducted functional studies of rhodopsin-like GPCR using transgenic lines of D. melanogaster and found that not only did the tolerance to permethrin insecticide increase in the transgenic lines of D. melanogaster, but the expression of Drosophila’s resistance the P450 361 genes CYP12d1 and CYP6a8 also increased. Similar functional studies using Drosophila transgenic lines have further confirmed the involvement of GPCRs in the regulation of P450 gene expression and insecticide resistance in L. dispar [66] and M. domestica [68]. Transposable elements (TEs), DNA sequences with the ability to change their position within a genome, have been reported in the evolutionary and adaptive processes of a variety of organisms [69,70] by introducing new cis-acting factors or renovating transcriptional networks to regulate the expression of genes, including metabolic detoxification-related cytochrome P450 genes in insects [71,72,73,74,75]. By screening the association of ETs with P450 genes in both Helicoverpa zea and D. melanogaster, Chen and Li [71] reported that several TEs were inserted into 5’ or 3′ flanking regions, the exons or the introns of several xenobiotic-metabolizing cytochrome P450 genes in both insect species, and proposed that TEs were selectively retained within or in close proximity to xenobiotic-metabolizing P450 genes. Chung et al. [72] reported the insertion of a long terminal repeat (LTR) of an Accord retrotransposon in the upstream region of the Cyp6g1 gene, the gene involved in DDT resistance in D. melanogaster. The authors demonstrated that the Accord LTR carried regulatory sequences that increased the expression of Cyp6g1 in detoxification tissues, implying their role in resistance P450 gene regulation and insecticide resistance. A similar result was reported by Schlenke and Begun [73], in which the presence of the insertion of a Doc transposable element was correlated with increased abundance of Cyp6g1 transcription. In the mosquito Culex quinquefasciatus, a partial CuRE1 (Culex Repetitive Element 1) transposable element has been found to be inserted directly upstream of CYP9M10, a cytochrome P450 gene involved in mosquito resistance to permethrin [74]. A study by Marsano et al. [75] provided an example of a cis-regulatory mechanism through the insertion of TE(s) at the 3’ untranslated region of genes. The above examples provide evidence of TE-mediated innovation in gene regulation. Further characterization of molecular mechanisms involved in this cis- or trans-regulation of resistance P450 genes will be essential in shedding light on the development of P450 detoxification-mediated insecticide resistance in insects.
Techniques such as in silico homology modeling and the molecular docking of proteins are allowing researchers to construct three-dimensional models of proteins, supporting efforts to develop a deeper understanding of the relationship between the structure of proteins such as P450s and substrates such as insecticides, and thus, provide reasonable explanations for the substrate specificities and metabolic specificities caused by allelic variations or mutations [1,76]. In our study, we employed this technique to further confirm the findings of our transgenic and in vitro metabolism study of CYP6A36, CYP6D10 and CYP4S24. Comparisons of the 3D models for permethrin and the potential PBalc and PBald docking positions among CYP6A36, CYP6A52, CYP6D10 and CYP4S24 not only support the metabolic capacity of these P450s for these substrates, but also indicates that CYP6A36, CYP6A52 and CYP6D10 are likely to produce similar metabolites. In CYP4S24, the low metabolic ability of permethrin observed may be due to the narrow channel opening to the heme iron and the small active site cavity (attributed to the protrusion of two amino acids in the helix I structure), which prevent the access of large molecules such as permethrin into the catalytic hotspot [41,77]. However, this structure does provide a good match for small molecules such as PBalc and PBald. This is probably the main reason why CYP4S24 exhibits the highest metabolic ability toward PBalc and PBald among these three P450s. The differences in metabolic activities among the various P450 enzymes of house flies may well be attributed to structural differences in the P450 proteins that facilitate substrate binding, conduction and activation.

4. Materials and Methods

4.1. House Fly Strains and Lines

Three house fly strains, five backcross (BC1) lines and five homozygous lines were used in this study. ALHF is a wild-type multiple-insecticide-resistant strain collected from a poultry farm in Alabama in 1998 [78]. This strain was further selected with permethrin, a pyrethroid insecticide, for six generations after collection to reach a high level of resistance, and has been maintained under biannual selection with permethrin [25,78]. Aabys is an insecticide-susceptible strain with recessive morphological markers including ali-curve (ac), aristapedia (ar), brown body (bwb), yellow eyes (ye) and snipped wings (snp) on autosomes 1, 2, 3, 4 and 5, respectively. CS is a wild-type insecticide-susceptible strain that has been maintained in the laboratory for more than five decades. Both the aabys and CS strains were originally obtained from Dr. J. G. Scott (Cornell University).
Five back-cross (BC1) lines were generated from a cross of ALHF females and aabys males with the following genotypes: ac/ac, +/ar, +/bwb, +/ye, +/sw (A2345); +/ac, ar/ar, +/bwb, +/ye, +/sw (A1345); +/ac, +/ar, bwb/bwb, +/ye, +/sw (A1245); +/ac, +/ar, +/bwb, ye/ye, +/sw (A1235); and +/ac, +/ar, +/bwb, +/ye, sw/sw (A1234). Five homozygous lines were also generated: ac/ac, +/+, +/+, +/+, +/+ (A2345); +/+, ar/ar, +/+, +/+, +/+ (A1345); (+/+, +/+, bwb/bwb, +/+, +/+ (A1245); +/+, +/+, +/+, ye/ye, +/+ (A1235); and +/+, +/+, +/+, +/+, snp/snp (A1234) [11,25]. The name of each line indicates which of its autosomes bear wild-type markers from ALHF. For instance, the A2345 strain has wild-type markers on autosomes 2, 3, 4 and 5 from ALHF and the mutant marker on autosome 1 from aabys.

4.2. Quantitative Real-Time PCR (qRT-PCR)

A total of 20 3-day-old adult female house flies from each of the three house fly strains and house fly lines had their RNA extracted for each experiment using the acidic guanidine thiocyanate-phenol-chloroform method [43]. Total RNA (0.5 µg/sample) from each house fly sample was reverse-transcribed using SuperScript II reverse transcriptase (Stratagene, San Diego, CA, USA) in a total volume of 20 μL. The quantity of cDNAs was measured using a spectrophotometer prior to qRT-PCR, which was performed using the SYBR Green Master Mix Kit and the ABI 7500 Real-Time PCR system (Applied Biosystems, Waltham, MA, USA). Each qRT-PCR reaction (15 μL final volume) contained 1× SYBR Green Master Mix, 1 μL of cDNA, and a P450 gene-specific primer pair (Table S1) designed according to each of the P450 gene sequences at a final concentration of 0.3–0.5 μM. All sampling, including the ‘no-template’ negative control, was performed in triplicate. The specificity of the PCR reactions was assessed by a melting curve analysis for each PCR reaction using Dissociation Curves software [79]. Relative expression levels for the P450 genes were calculated via the 2−ΔΔCT method using SDS RQ software [80]. The β-actin gene and ribosomal protein S3 (RPS3) were used as endogenous controls, to normalize the expression of the target genes [81]. Each experiment was repeated either three or four times with different preparations of RNA samples. The statistical significance of the gene expressions was calculated using a Student’s t-test for all pairwise sample comparisons, and one-way analysis of variance (ANOVA) was utilized for multiple sample comparisons (SPSS 19.0 software); a value of p ≤ 0.05 was considered statistically significant.

4.3. Autosomal Linkage Analysis of the P450 Genes in M. domestica

Five back-cross (BC1) lines (ac/ac, +/ar, +/bwb, +/ye, +/sw (A2345); +/ac, ar/ar, +/bwb, +/ye, +/sw (A1345); +/ac, +/ar, bwb/bwb, +/ye, +/sw (A1245); +/ac, +/ar, +/bwb, ye/ye, +/sw (A1235); and +/ac, +/ar, +/bwb, +/ye, sw/sw (A1234)) were used in this study. The genotype of each BC1 line was homozygous for the recessive mutant allele from aabys and heterozygous for the dominant wild-type alleles from ALHF. Allele-specific PCR was conducted using a two-round PCR strategy as described by Liu and Scott [26]. For the first PCR reaction, allele-independent primer pairs (Table S1) were used to generate P450 cDNA fragments as follows: the first PCR solution with the cDNA template and a primer pair were heated to 95 °C for 3 min, followed by 35 cycles of 95 °C for 30 s, 60 °C for 30 s and 72 °C for 1 min, then, 72 °C for 10 min. The second PCR reaction utilized 0.5 µL of the first-round PCR reaction solution and the allele-specific primer pair (Table S1), designed based on the specific sequence of the genes from ALHF by placing a specific nucleotide polymorphism at the 3′ end of the primer to permit preferential amplification of the allele from ALHF. This second PCR reaction sample was heated to 95 °C for 3 min, followed by 35 cycles of 95 °C for 30 s, 62 °C for 30 s and 72 °C for 30 s, then, 72 °C for 10 min. Each experiment was repeated three times with different mRNAs, and the PCR products were sequenced at least once for each gene.

4.4. Construction of Transgenic D. melanogaster and Toxicity of Permethrin for the Transgenic Lines

The full lengths of the candidate up-regulated genes from the ALHF strain of M. domestica were amplified from cDNA of M. domestica using platinum Taq DNA polymerase High Fidelity (Invitrogen) with specific primer pairs (Table S1) based on the 5′ and 3′ end sequences of the genes. The PCR products of the full-length genes were purified using a QIAquick Gel Extraction Kit (Qiagen). The purified PCR products were ligated into pCR 2.1 vector using the Original TA Cloning kit (Invitrogen) as described by the manufacturer. The full lengths of the genes were cloned in One Shot TOP10F’ cells using the One Shot TOP10F’ Chemically Competent E. coli kit (Invitrogen). Cloning and sequence analyses of the cDNAs were repeated at least three times and three TA clones from each replication were verified via sequencing. The clones were then sub-cloned into a pUASTattB vector (a gift from Dr. Johannes Bischof, University of Zurich). The plasmid of each pUASTattB-up-regulated gene was transformed into the germ line of the M{vas-int.Dm}ZH-2A, M{3×P3-RFP.attP’}ZH-58A strain of D. melanogaster (Bloomington stock #24484), resulting in site-specific integration on chromosome 2R (Rainbow Transgenic Flies Inc., Camarillo, CA, USA). These flies were then reciprocally crossed against a W1118 strain to obtain transformants with the orange eye phenotype, then, balanced against the D. melanogaster balancer strain w[1118]/Dp(1;Y)y[+]; sna[Sco]/CyO, P{ry[+t7.2] = sevRas1.V12}FK1 (Bloomington stock #6312) to generate a homozygous line containing the cytochrome P450 transgene. The insertion of the up-regulated genes in the transgenic fruit fly lines was further confirmed using RT-PCR. The resulting homozygous line was crossed with the GAL4-expressing D. melanogaster strain P{Act5C-GAL4}17bFO1 (Bloomington stock #3954), which expresses GAL4 under control of the Act5C promoter, resulting in ubiquitous non-tissue-specific expression. The F1 generation of these crosses expressed GAL4 and contained a single copy of the cytochrome P450 transgene that was under the control of the UAS enhancer.
Permethrin toxicity bioassays were then conducted 2–3 days post-eclosion of female Drosophila of the F1 UAS-GAL4 crosses to examine the toxicity of permethrin to transgenic flies. Briefly, serial concentrations of permethrin solution in acetone ranging from 3 ng/µL to 50 ng/µL, which gave >0 and <100% mortality to the tested insects, were prepared, and 200 µL of each permethrin solution was evenly coated on the inside of individual 20 mL glass scintillation vials. Fifteen female flies were transferred to each of the prepared vials, which were then plugged with cotton balls soaked with 5% sucrose. The vials for the control groups were coated with acetone alone and plugged with identical 5% sucrose-soaked cotton balls. Mortality was scored after 24 h exposure to permethrin. Once the concentration range had been determined, a single concentration (25 ng/µL) that showed a significant difference in mortality between the transgenic and non-transgenic flies was selected for further permethrin bioassays. For each bioassay, a total of four vials each containing 15 2–3-day-old adult female D. melanogaster were treated with 200 µL of 25 ng/µL permethrin solution for each strain. Each bioassay was independently replicated three times using only flies that exhibited the correct morphological marker (orange eyes). The D. melanogaster strain M{vas-int.Dm}ZH-2A, M{3×P3-RFP.attP’}ZH-58A containing the empty pUAST vector-donated insert, but no transgene from M. domestica, was used as the control reference strain. All tests were run at 27 °C and mortality was assessed after 24 h. Bioassay data were pooled and probit analysis was conducted. Significant differences in the resistance levels of the D. melanogaster lines were determined based on the non-overlap of 95% confidence intervals. All D. melanogaster were reared on Jazz-Mix Drosophila food (Fisher Scientific, Kansas City, MO, USA) at 25 ± 2 °C under a photoperiod of 12:12 (L:D) h.

4.5. Recombinant Protein Production

The recombinant P450s and NADPH-cytochrome P450 reductase (CPR) proteins used in this study were produced in Spodoptera frugiperda (Sf9) insect cells using a BaculoDirectTM baculovirus-mediated expression system (Invitrogen). Briefly, the full-lengths of CYP6A36, CYP6D10, CYP4S24 and NADPH-cytochrome P450 reductase (CPR) were amplified from the house fly ALHF strain with specific primers (Table S1) and ligated with the BaculoDirectTM Linear DNA via an in vitro Lambda Recombination (LR) reaction. The constructed recombinant baculovirus for each gene was then transfected to Sf9 cells using CellfectinR II reagents (Invitrogen) to produce the baculovirus stock solution for further infection. The titers of baculoviruses were later measured via a plaque-forming assay, and a titer of ~2 × 108 pfu/mL P2 virus was used as a final stock to infect cells for the large-scale expression of targeted proteins. According to protocols employed in Culex mosquitoes for their P450 characterizations [28,75], the co-expression of house fly P450s and CPR (P450/CPR) in Sf9 cells was performed with an optimal multiplicity of infection (MOI) ratio of 10:1, supplemented with 1 μg/mL hemin and 0.1 mM 5-ALA, to achieve the highest P450 and catalytic activities [28,75]. The cell lysate proteins were harvested 72 h post-infection and centrifuged at 1000× g for 10 min at 4 °C. The cell pellets were washed twice using ice-cold PBS buffer (pH 7.4), and then, re-suspended in homogenization buffer (0.1 M phosphate (pH 7.4), 1 mM EDTA, 0.5 mM PMSF and 0.25 M sucrose). Subsequently, the dissolved cell lysate was then homogenized using a sonicator and ultracentrifuged at 37,000× g for 60 min at 4 °C. The microsomal protein pellets were isolated and dissolved in re-suspension buffer (0.1 M phosphate buffer, pH = 7.4, 1 mM EDTA, 0.1 mM DTT, 1 mM PMSF and 20% glycerol stock) and stored at −80 °C for use in in vitro studies.

4.6. P450/CPR Activity Determination

The activities of CYP6A36, CYP6D10 and CYP4S24, with co-expressed CPR as their electron donor, were measured via a 7-ethoxycoumarin O-deethylation (ECOD) assay with modifications [28]. Briefly, a mixture of 5 μL microsomal protein of each P450/CPR, 1 μL 7-ethoxycoumarin (50 mM dissolved in acetone), and 85 μL sodium phosphate buffer (0.1 M, pH = 7.5) was added to each well of a 96-well microplate, and 10 μL of NADPH was added to initiate the reaction; then, it was incubated at 30 °C for 30 min with orbital shaking, after which 10 μL oxidized glutathione (30 mM) and 10 μL glutathione reductase (5U) were added to remove the NADPH. After 10 min incubation at room temperature, the reaction was stopped by adding 120 μL of 50% acetonitrile in TRIZMA-base buffer (0.05 M, pH = 10). The amount of 7-hydroxycoumarin (7-OH) released during incubation was quantified using a Cytation 3 imagine reader (BioTekUSA) with 390 nm excitation and 450 nm emission filters. A well containing Sf9 cell microsomal protein only, and a well without microsomes, served as controls. A standard curve for 7-hydroxycoumarin (with substrate concentrations ranging from 0.1 to 1.2 mM) was used to calculate the kinetic parameters.

4.7. P450/CPR-Mediated In Vitro Metabolism of Permethrin, 3-Phenoxybenzyl Alcohol (PBalc) and 3-Phenoxybenzaldehyde (PBald)

The method and procedure used for the metabolism experiments were as described in our previous studies [28,82,83]. The in vitro metabolism of permethrin and its metabolites, PBalc and PBald, was conducted by incubating chemicals (20 μM) dissolved in Tris-HCl buffer (0.2 M, pH 7.4)/acetonitrile (1:1, v:v) with 0.2 μM of different recombinant P450/CPR microsomes and an NADPH-regenerating system (1 mM NADP+, 0.25 mM MgCl2) in a total volume of 700 μL; the mixture was incubated at 30 °C with 100 rpm orbital shaking and quenched after 120 min by adding 700 μL ice-cold acetonitrile. Samples were then incubated with shaking for an additional 20 min and centrifuged at 16,000× g for 10 min. The control group contained the same components, without NADPH added. The supernatant from each sample was then filtered through 0.45 μm membranes and the amounts of the chemicals were analyzed via HPLC using a Waters 2695 High-Performance Liquid Chromatograph System (Milford, MA, USA) with the conditions employed in previous studies [28]. The permethrin amount was measured as the integrated area under the two peaks of trans-permethrin and cis-permethrin [38], and calculated based on the standard curves of each chemical. Reactions were performed in triplicate and a paired t-test of sample reactions (with NADPH) vs. the control (without NADPH) was performed for the statistical measurements of substrate depletion following the method described by Gong et al. [28,75].

4.8. In Silico Modeling and Substrate Docking Analysis

In silico 3D structural modeling was performed by the I-TASSER server using a combination of threading and ab initio modeling [84,85,86]. Five models were predicted by I-TASSER for each P450, with the top-scoring model being submitted to the FG-MD server for fragment guided molecular dynamics structure refinement [87]. The model quality was controlled by Ramachandran plots, generated using Procheck (http://services.mbi.ucla.edu/SAVES/; accessed on 15 October 2021) [87] and ProSA-web (https://prosa.services.came.sbg.ac.at/prosa.php; accessed on 15 October 2021) [88,89]. The P450 channels were calculated using CAVER 3.0 (http://caver.cz/index.php; accessed on 15 October 2021) [88,89]. The P450 channels were calculated using CAVER 3.0 (http://caver.cz/index.php; accessed on 15 October 2021) [90,91]; those permitting the passage of a sphere with a maximal radius greater than 1.2 Å were identified, tabulated and named according to the nomenclature proposed by Cojocaru et al. [36]. The volume of the substrate binding cavity was characterized by VOIDOO using a 1.4 Å probe [92], and the proteins and ligands were prepared for docking using Autodock Tools v1.5.6 (http://mgltools.scripps.edu/downloads; accessed on 15 October 2021). Molecular docking was performed using Autodock 4.2. [93]. The ligand permethrin structures were retrieved from the ZINC database [94]. For all dockings, a search space with a grid box of 60 × 60 × 60 Å, centered on the heme iron, was set that corresponded to the substrate recognition sites (SRSs), following those for the CYP2 family proposed by Gohoh [95]. The figures were produced using Pymol (http://www.pymol.org/; accessed on 15 October 2021) [96].

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms24043170/s1.

Author Contributions

Conception and design of the study: M.L. and N.L. Performance of the experiments: M.L., X.F., W.R.R. and F.T. Preparation of the materials: N.L. Writing and discussion of the paper: M.L., X.F., W.R.R., F.T. and N.L. All authors reviewed the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This project was supported by the following awards: the AAES Hatch/Multistate Grants ALA015-1-10026 and ALA015-1-16009 to N.L.

Institutional Review Board Statement

Not applicable. The study did not require ethical approval.

Informed Consent Statement

Not applicable. This article does not contain any studies with human participants or animals performed by any of the authors.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author, N.L., upon request.

Acknowledgments

The authors thank Jeff Scott for his valuable comments and discussion on the project.

Conflicts of Interest

The authors declare no competing interest.

References

  1. Liu, N.; Li, M.; Gong, Y.; Liu, F.; Li, T. Cytochrome P450s—Their expression, regulation, and role in insecticide resistance. Pestic. Biochem. Phys. 2015, 120, 77–81. [Google Scholar] [CrossRef] [PubMed]
  2. Li, T.; Liu, L.; Zhang, L.; Liu, N. Role of g-protein-coupled receptor-related genes in insecticide resistance of the mosquito, Culex quinquefasciatus. Sci. Rep. 2014, 4, 6474. [Google Scholar] [CrossRef] [PubMed]
  3. Liu, N. Insecticide resistance in mosquitoes: Impact, mechanisms, and research directions. Annu. Rev. Entomol. 2015, 60, 537–559. [Google Scholar] [CrossRef]
  4. Scott, J.G. Cytochromes P450 and insecticide resistance. Insect Biochem. Mol. Biol. 1999, 29, 757–777. [Google Scholar] [CrossRef] [PubMed]
  5. Brun-Barale, A.; Héma, O.; Martin, T.; Suraporn, S.; Audant, P.; Sezutsu, H.; Feyereisen, R. Multiple P450 genes overexpressed in deltamethrin-resistant strains of Helicoverpa armigera. Pest Manag. Sci. 2010, 66, 900–909. [Google Scholar] [CrossRef]
  6. Yang, T.; Liu, N. Genome analysis of cytochrome p450s and their expression profiles in insecticide resistant mosquitoes, Culex quinquefasciatus. PLoS ONE 2011, 6, e29418. [Google Scholar] [CrossRef]
  7. Gong, Y.; Li, T.; Zhang, L.; Gao, X.; Liu, N. Permethrin induction of multiple cytochrome p450 genes in insecticide resistant mosquitoes, Culex quinquefasciatus. Int. J. Biol. Sci. 2013, 9, 863–871. [Google Scholar] [CrossRef]
  8. Riveron, J.M.; Irving, H.; Ndula, M.; Barnes, K.G.; Ibrahim, S.S.; Paine, M.J.I.; Wondji, C.S. Directionally selected cytochrome P450 alleles are driving the spread of pyrethroid resistance in the major malaria vector Anopheles funestus. Proc. Natl. Acad. Sci. USA 2013, 110, 252–257. [Google Scholar] [CrossRef]
  9. Ma, K.; Tang, Q.; Zhang, B.; Liang, P.; Wang, B.; Gao, X. Overexpression of multiple cytochrome P450 genes associated with sulfoxaflor resistance in Aphis gossypii Glover. Pestic. Biochem. Phys. 2019, 157, 204–210. [Google Scholar] [CrossRef]
  10. Yang, T.; Li, T.; Feng, X.; Li, M.; Liu, S.; Liu, N. Multiple cytochrome P450 genes: Conferring high levels of permethrin resistance in mosquitoes, Culex quinquefasciatus. Sci. Rep. 2021, 11, 9041. [Google Scholar] [CrossRef]
  11. Li, M.; Reid, W.R.; Zhang, L.; Scott, J.G.; Gao, X.; Kristensen, M.; Liu, N. A whole transcriptomal linkage analysis of gene co-regulation in insecticide resistant house flies, Musca domestica. BMC Genom. 2013, 14, 803. [Google Scholar] [CrossRef] [PubMed]
  12. Scott, J.G.; Warren, W.C.; Beukeboom, L.W.; Bopp, D.; Clark, A.G.; Giers, S.D.; Hediger, M.; Jones, A.K.; Kasai, S.; Leichter, C.A.; et al. Genome of the house fly, Musca domestica L., a global vector of diseases with adaptations to a septic environment. Genome Biol. 2014, 15, 466. [Google Scholar] [CrossRef] [PubMed]
  13. Pavlidi, N.; Vontas, J.; Van Leeuwen, T. The role of glutathione S-transferases (GSTs) in insecticide resistance in crop pests and disease vectors. Curr. Opin. Insect Sci. 2018, 27, 97–102. [Google Scholar] [CrossRef] [PubMed]
  14. Rosen, R.; Lebedev, G.; Kontsedalov, S.; Ben-Yakir, D.; Ghanim, M. A de novo transcriptomics approach reveals genes involved in thrips tabaci resistance to spinosad. Insects 2021, 12, 67. [Google Scholar] [CrossRef] [PubMed]
  15. Main, B.J.; Everitt, A.; Cornel, A.J.; Hormozdiari, F.; Lanzaro, G.C. Genetic variation associated with increased insecticide resistance in the malaria mosquito, Anopheles coluzzii. Parasites Vectors 2018, 11, 225. [Google Scholar] [CrossRef]
  16. Simma, E.A.; Dermauw, W.; Balabanidou, V.; Snoeck, S.; Bryon, A.; Clark, R.M.; Yewhalaw, D.; Vontas, J.; Duchateau, L.; Van Leeuwen, T. Genome-wide gene expression profiling reveals that cuticle alterations and P450 detoxification are associated with deltamethrin and DDT resistance in Anopheles arabiensis populations from Ethiopia. Pest Manag. Sci. 2019, 75, 1808–1818. [Google Scholar] [CrossRef]
  17. Mugenzi, L.M.J.; Menze, B.D.; Tchouakui, M.; Wondji, M.J.; Irving, H.; Tchoupo, M.; Hearn, J.; Weedall, G.D.; Riveron, J.M.; Wondji, C.S. Cis-regulatory CYP6P9b P450 variants associated with loss of insecticide-treated bed net efficacy against Anopheles funestus. Nat. Commun. 2019, 10, 4652. [Google Scholar] [CrossRef]
  18. Nayduch, D.; Burrus, R.G. Flourishing in filth: House fly–microbe interactions across life history. Ann. Entomol. Soc. Am. 2017, 110, 6–18. [Google Scholar] [CrossRef]
  19. Freeman, J.C.; Ross, D.H.; Scott, J.G. Insecticide resistance monitoring of house fly populations from the United States. Pestic. Biochem. Phys. 2019, 158, 61–68. [Google Scholar] [CrossRef]
  20. Feyereisen, R. Insect cytochrome P450. In Comprehensive Molecular Insect Science; Elsevier: Amsterdam, The Netherlands, 2005; Volume 4, pp. 1–77. [Google Scholar]
  21. Li, X.; Schuler, M.A.; Berenbaum, M.R. Molecular mechanisms of metabolic resistance to synthetic and natural xenobiotics. Annu. Rev. Entomol. 2007, 52, 231–253. [Google Scholar] [CrossRef]
  22. Pavlidi, N.; Monastirioti, M.; Daborn, P.; Livadaras, I.; Van Leeuwen, T.; Vontas, J. Transgenic expression of the Aedes aegypti CYP9J28 confers pyrethroid resistance in Drosophila melanogaster. Pestic. Biochem. Phys. 2012, 104, 132–135. [Google Scholar] [CrossRef]
  23. Yang, Y.; Yue, L.; Chen, S.; Wu, Y. Functional expression of Helicoverpa armigera CYP9A12 and CYP9A14 in Saccharomyces cerevisiae. Pestic. Biochem. Phys. 2008, 92, 101–105. [Google Scholar] [CrossRef]
  24. Zhu, F.; Liu, N. Differential expression of CYP6A5 and CYP6A5v2 in pyrethroid-resistant house flies, Musca domestica. Arch. Insect Biochem. Physiol. 2008, 67, 107–119. [Google Scholar] [CrossRef] [PubMed]
  25. Tian, L.; Cao, C.; He, L.; Li, M.; Zhang, L.; Zhang, L.; Liu, H.; Liu, N. Autosomal interactions and mechanisms of pyrethroid resistance in house flies, Musca domestica. Int. J. Biol. Sci. 2011, 7, 902–911. [Google Scholar] [CrossRef]
  26. Liu, N.; Scott, J.G. Genetics of resistance to pyrethroid insecticides in the house fly, Musca domestica. Pestic. Biochem. Phys. 1995, 52, 116–124. [Google Scholar] [CrossRef]
  27. Liu, N.; Yue, X. Genetics of pyrethroid resistance in a strain (ALHF) of house flies (Diptera: Muscidae). Pestic. Biochem. Phys. 2001, 70, 151–158. [Google Scholar] [CrossRef]
  28. Gong, Y.; Li, T.; Feng, Y.; Liu, N. The function of two P450s, CYP9M10 and CYP6AA7, in the permethrin resistance of Culex quinquefasciatus. Sci. Rep. 2017, 7, 587. [Google Scholar] [CrossRef] [PubMed]
  29. Chandor-Proust, A.; Bibby, J.; Régent-Kloeckner, M.; Roux, J.; Guittard-Crilat, E.; Poupardin, R.; Riaz, M.A.; Paine, M.; Dauphin-Villemant, C.; Reynaud, S.; et al. The central role of mosquito cytochrome P450 CYP6Zs in insecticide detoxification revealed by functional expression and structural modelling. Biochem. J. 2013, 455, 75–85. [Google Scholar] [CrossRef]
  30. Godin, S.J.; Scollon, E.J.; Hughes, M.F.; Potter, P.M.; DeVito, M.J.; Ross, M.K. Species differences in the in vitro metabolism of deltamethrin and esfenvalerate: Differential oxidative and hydrolytic metabolism by humans and rats. Drug Metab. Dispos. 2006, 34, 1764–1771. [Google Scholar] [CrossRef]
  31. Nakamura, Y.; Sugihara, K.; Sone, T.; Isobe, M.; Ohta, S.; Kitamura, S. The in vitro metabolism of a pyrethroid insecticide, permethrin, and its hydrolysis products in rats. Toxicology 2007, 235, 176–184. [Google Scholar] [CrossRef] [PubMed]
  32. Ross, M.K.; Borazjani, A.; Edwards, C.C.; Potter, P.M. Hydrolytic metabolism of pyrethroids by human and other mammalian carboxylesterases. Biochem. Pharmacol. 2006, 71, 657–669. [Google Scholar] [CrossRef] [PubMed]
  33. Somwang, P.; Yanola, J.; Suwan, W.; Walton, C.; Lumjuan, N.; Prapanthadara, L.; Somboon, P. Enzymes-based resistant mechanism in pyrethroid resistant and susceptible Aedes aegypti strains from northern Thailand. Parasitol. Res. 2011, 109, 531–537. [Google Scholar] [CrossRef] [PubMed]
  34. Takaku, T.; Mikata, K.; Matsui, M.; Nishioka, K.; Isobe, N.; Kaneko, H. In vitro metabolism of trans-permethrin and its major metabolites, PBalc and PBacid, in humans. J. Agric. Food Chem. 2011, 59, 5001–5005. [Google Scholar] [CrossRef] [PubMed]
  35. Tange, S.; Fujimoto, N.; Uramaru, N.; Sugihara, K.; Ohta, S.; Kitamura, S. In vitro metabolism of cis- and trans-permethrin by rat liver microsomes, and its effect on estrogenic and anti-androgenic activities. Environ. Toxicol. Pharmacol. 2014, 37, 996–1005. [Google Scholar] [CrossRef] [PubMed]
  36. Cojocaru, V.; Winn, P.J.; Wade, R.C. The ins and outs of cytochrome P450s. Biochim. Biophys. Acta 2007, 1770, 390–401. [Google Scholar] [CrossRef]
  37. Otyepka, M.; Berka, K.; Anzenbacher, P. I Is there a relationship between the substrate preferences and structural flexibility of cytochromes P450? Curr. Drug Metab. 2012, 13, 130–142. [Google Scholar] [CrossRef]
  38. Stevenson, B.J.; Bibby, J.; Pignatelli, P.; Muangnoicharoen, S.; O’Neill, P.M.; Lian, L.-Y.; Müller, P.; Nikou, D.; Steven, A.; Hemingway, J.; et al. Cytochrome P450 6M2 from the malaria vector Anopheles gambiae metabolizes pyrethroids: Sequential metabolism of deltamethrin revealed. Insect Biochem. Mol. Biol. 2011, 41, 492–502. [Google Scholar] [CrossRef]
  39. Zhu, F.; Moural, T.W.; Shah, K.; Palli, S.R. Integrated analysis of cytochrome P450 gene superfamily in the red flour beetle, Tribolium castaneum. BMC Genom 2013, 14, 174. [Google Scholar] [CrossRef]
  40. Zimmer, C.T.; Bass, C.; Williamson, M.S.; Kaussmann, M.; Wölfel, K.; Gutbrod, O.; Nauen, R. Molecular and functional characterization of CYP6BQ23, a cytochrome P450 conferring resistance to pyrethroids in European populations of pollen beetle, Meligethes aeneus. Insect Biochem. Mol. Biol. 2014, 45, 18–29. [Google Scholar] [CrossRef]
  41. Ibrahim, S.S.; Riveron, J.M.; Stott, R.; Irving, H.; Wondji, C.S. The cytochrome P450 CYP6P4 is responsible for the high pyrethroid resistance in knockdown resistance-free Anopheles arabiensis. Insect Biochem. Mol. Biol. 2016, 68, 23–32. [Google Scholar] [CrossRef] [Green Version]
  42. Karunker, I.; Morou, E.; Nikou, D.; Nauen, R.; Sertchook, R.; Stevenson, B.J.; Paine, M.J.; Morin, S.; Vontas, J. Structural model and functional characterization of the Bemisia tabaci CYP6CM1vQ, a cytochrome P450 associated with high levels of imidacloprid resistance. Insect Biochem. Mol. Biol. 2009, 39, 697–706. [Google Scholar] [CrossRef]
  43. Liu, N.; Scott, J.G. Phenobarbital induction of CYP6D1 is due to a trans acting factor on autosome 2 in house flies, Musca domestica. Insect Mol. Biol. 1997, 6, 77–81. [Google Scholar] [CrossRef] [PubMed]
  44. Carino, F.A.; Koener, J.F.; Plapp Jr, F.W.; Feyereisen, R. Constitutive overexpression of the cytochrome P450 gene CYP6A1 in a house fly strain with metabolic resistance to insecticides. Insect Biochem. Mol. Biol. 1994, 24, 411–418. [Google Scholar] [CrossRef]
  45. Pavek, P.; Dvorak, Z. Xenobiotic-induced transcriptional regulation of xenobiotic metabolizing enzymes of the cytochrome P450 superfamily in human extrahepatic tissues. Curr. Drug Metab. 2008, 9, 129–143. [Google Scholar] [CrossRef] [PubMed]
  46. Zhu, F.; Parthasarathy, R.; Bai, H.; Woithe, K.; Kaussmann, M.; Nauen, R.; Harrison, D.A.; Palli, S.R. A brain-specific cytochrome P450 responsible for the majority of deltamethrin resistance in the QTC279 strain of Tribolium castaneum. Proc. Natl. Acad. Sci. USA 2010, 107, 8557–8562. [Google Scholar] [CrossRef]
  47. Wang, H.; Shi, Y.; Wang, L.; Liu, S.; Wu, S.; Yang, Y.; Feyereisen, R.; Wu, Y. CYP6AE gene cluster knockout in Helicoverpa armigera reveals role in detoxification of phytochemicals and insecticides. Nat. Commun. 2018, 9, 4820. [Google Scholar] [CrossRef] [PubMed]
  48. Liu, N.; Scott, J.G. Genetic analysis of factors controlling high-level expression of cytochrome P450, CYP6D1, cytochrome b5, P450 reductase, and monooxygenase activities in LPR house flies, Musca domestica. Biochem. Genet. 1996, 34, 133–148. [Google Scholar] [CrossRef] [PubMed]
  49. Maitra, S.; Dombrowski, S.M.; Waters, L.C.; Ganguly, R. Three Second Chromosome-Linked Clustered Cyp6 Genes Show Differential Constitutive and Barbital-Induced Expression in DDT-Resistant and Susceptible Strains of Drosophila Melanogaster. Gene 1996, 180, 165–171. [Google Scholar] [CrossRef]
  50. Maitra, S.; Dombrowski, S.M.; Basu, M.; Raustol, O.; Waters, L.C.; Ganguly, R. Factors on the third chromosome affect the level of Cyp6a2 and Cyp6a8 expression in Drosophila melanogaster. Gene 2000, 248, 147–156. [Google Scholar] [CrossRef]
  51. Fotoukkiaii, S.M.; Wybouw, N.; Kurlovs, A.H.; Tsakireli, D.; Pergantis, S.A.; Clark, R.M.; Vontas, J.; Leeuwen, T.V. High-resolution genetic mapping reveals cis-regulatory and copy number variation in loci associated with cytochrome P450-mediated detoxification in a generalist arthropod pest. PLoS Genet. 2021, 17, e1009422. [Google Scholar] [CrossRef]
  52. Papapostolou, K.M.; Riga, M.; Samantsidis, G.-R.; Skoufa, E.; Balabanidou, V.; Van Leeuwen, T.; Vontas, J. Over-expression in cis of the midgut P450 CYP392A16 contributes to abamectin resistance in Tetranychus urticae. Insect Biochem. Mol. Biol. 2022, 142, 103709. [Google Scholar] [CrossRef] [PubMed]
  53. Chung, H.; Boey, A.; Lumb, C.; Willoughby, L.; Batterham, P.; Daborn, P.J. Induction of a detoxification gene in Drosophila melanogaster requires an interaction between tissue specific enhancers and a novel cis-regulatory element. Insect Biochem. Mol. Biol. 2011, 41, 863–871. [Google Scholar] [CrossRef] [PubMed]
  54. Li, T.; Cao, C.; Yang, T.; Zhang, L.; He, L.; Xi, Z.; Bian, G.; Liu, N. A G-protein-coupled receptor regulation pathway in cytochrome P450-mediated permethrin-resistance in mosquitoes, Culex quinquefasciatus. Sci. Rep. 2016, 5, 17772. [Google Scholar] [CrossRef] [PubMed]
  55. Shi, Y.; Sun, S.; Zhang, Y.; He, Y.; Du, M.; O’Reilly, A.O.; Wu, S.; Yang, Y.; Wu, Y. Single amino acid variations drive functional divergence of cytochrome P450s in Helicoverpa species. Insect Biochem. Mol. Biol. 2022, 146, 103796. [Google Scholar] [CrossRef]
  56. McDonnell, C.M.; Brown, R.P.; Berenbaum, M.R.; Schuler, M.A. Conserved regulatory elements in the promoters of two allelochemical-inducible cytochrome P450 genes differentially regulate transcription. Insect Biochem. Mol. Biol. 2004, 34, 1129–1139. [Google Scholar] [CrossRef]
  57. Mohammed, B.R.; Kawe, M.S.; Nwachukwu, O.M.; Jegede, O.C.; Agbede, R.I.; Finn, R.D. Understanding the mechanisms involved in the regulation of cytochrome P450 gene expression in Drosophila melanogaster (Diptera: Drosophilidae). Entomol. Ornithol. Herpetol. 2017, 6, 1. [Google Scholar] [CrossRef]
  58. Li, T.; Liu, N. Regulation of P450-mediated permethrin resistance in Culex quinquefasciatus by the GPCR/Gαs/AC/cAMP/PKA signaling cascade. Biochem. Biophys. Rep. 2017, 12, 12–19. [Google Scholar] [CrossRef]
  59. Li, T.; Liu, N. The function of g-protein-coupled receptor-regulatory cascade in southern house mosquitoes (Diptera: Culicidae). J. Med. Entomol. 2018, 55, 862–870. [Google Scholar] [CrossRef]
  60. Kalsi, M.; Palli, S.R. Transcription Factors, CncC and Maf, Transcription factors, CncC and Maf, regulate expression of CYP6BQ genes responsible for deltamethrin resistance in Tribolium castaneum. Insect Biochem. Mol. Biol. 2015, 65, 47–56. [Google Scholar] [CrossRef]
  61. Pu, J.; Sun, H.; Wang, J.; Wu, M.; Wang, K.; Denholm, I.; Han, Z. Multiple cis-acting elements involved in up-regulation of a cytochrome P450 gene conferring resistance to deltamethrin in smal brown planthopper, Laodelphax striatellus (Fallén). Insect Biochem. Mol. Biol. 2016, 78, 20–28. [Google Scholar] [CrossRef]
  62. Palli, S.R. CncC/Maf-mediated xenobiotic response pathway in insects. Arch. Insect Biochem. Physiol. 2020, 104, e21674. [Google Scholar] [CrossRef]
  63. Smith, L.B.; Tyagi, R.; Kasai, S.; Scott, J.G. CYP-mediated permethrin resistance in Aedes aegypti and evidence for trans-regulation. PLoS Negl. Trop. Dis. 2018, 12, e0006933. [Google Scholar] [CrossRef] [PubMed]
  64. Yang, X.; Deng, S.; Wei, X.; Yang, J.; Zhao, Q.; Yin, C.; Du, T.; Guo, Z.; Xia, J.; Yang, Z.; et al. MAPK-directed activation of the whitefly transcription factor CREB leads to P450-mediated imidacloprid resistance. Proc. Natl. Acad. Sci. USA 2020, 117, 10246–10253. [Google Scholar] [CrossRef] [PubMed]
  65. Hu, B.; Huang, H.; Hu, S.; Ren, M.; Wei, Q.; Tian, X.; Elzaki, M.E.A.; Bass, C.; Su, J.; Palli, S.R. Changes in both trans- and cis-regulatory elements mediate insecticide resistance in a lepidopteron pest, Spodoptera exigua. PLoS Genet. 2021, 17, e1009403. [Google Scholar] [CrossRef] [PubMed]
  66. Cao, C.; Sun, L.; Du, H.; Moural, T.W.; Bai, H.; Liu, P.; Zhu, F. Physiological functions of a methuselah-like G protein coupled receptor in Lymantria dispar Linnaeus. Pestic. Biochem. Phys. 2019, 160, 1–10. [Google Scholar] [CrossRef] [PubMed]
  67. Sun, Y.; Zou, P.; Yu, X.-Y.; Chen, C.; Yu, J.; Shi, L.-N.; Hong, S.-C.; Zhou, D.; Chang, X.-L.; Wang, W.-J.; et al. Functional characterization of an arrestin gene on insecticide resistance of Culex pipiens pallens. Parasites Vectors 2012, 5, 134. [Google Scholar] [CrossRef]
  68. Ma, Z.; Zhang, Y.; You, C.; Zeng, X.; Gao, X. The role of G protein-coupled receptor-related genes in cytochrome P450-mediated resistance of the house fly, Musca domestica (Diptera: Muscidae), to imidacloprid. Insect Mol. Biol. 2020, 29, 92–103. [Google Scholar] [CrossRef]
  69. Sundaram, V.; Wysocka, J. Transposable elements as a potent source of diverse cis-regulatory sequences in mammalian genomes. Phil. Trans. R. Soc. 2020, 375, 20190347. [Google Scholar] [CrossRef]
  70. Moschetti, R.; Palazzo, A.; Lorusso, P.; Viggiano, L.; Marsano, R.M. “What you need, baby, I got it”: Transposable elements as suppliers of cis-operating sequences in Drosophila. Biology 2020, 9, 25. [Google Scholar] [CrossRef]
  71. Chen, S.; Li, X. Transposable elements are enriched within or in close proximity to xenobiotic-metabolizing cytochrome P450 genes. BMC Evol. Biol. 2007, 7, 46. [Google Scholar] [CrossRef] [Green Version]
  72. Chung, H.; Bogwitz, M.R.; McCart, C.; Andrianopoulos, A.; Ffrench-Constant, R.H.; Batterham, P.; Daborn, P.J. Cis-regulatory elements in the Accord retrotransposon result in tissue-specific expression of the Drosophila melanogaster insecticide resistance gene Cyp6g1. Genetics 2007, 175, 1071–1977. [Google Scholar] [CrossRef] [PubMed]
  73. Schlenke, T.A.; Begun, D.J. Strong selective sweep associated with a transposon insertion in Drosophila simulans. Proc. Natl. Acad. Sci. USA 2004, 101, 1626–1631. [Google Scholar] [CrossRef] [PubMed]
  74. Wilding, C.S.; Smith, I.; Lynd, A.; Yawson, A.E.; Weetman, D.; Paine, M.J.; Donnelly, M.J. A cis-regulatory sequence driving metabolic insecticide resistance in mosquitoes: Functional characterization and signatures of selection. Insect Biochem. Mol. Biol. 2012, 42, 699–707. [Google Scholar] [CrossRef] [PubMed]
  75. Marsano, R.M.; Caizzi, R.; Moschetti, R.; Junakovic, N. Evidence for a functional interaction between the Bari1 transposable element and the cytochrome P450 cyp12a4 gene in Drosophila melanogaster. Gene 2005, 357, 122–128. [Google Scholar] [CrossRef]
  76. Schuler, M.A.; Berenbaum, M.R. Structure and function of cytochrome P450s in insect adaptation to natural and synthetic toxins: Insights gained from molecular modeling. J. Chem. Ecol. 2013, 39, 1232–1245. [Google Scholar] [CrossRef]
  77. Lertkiatmongkol, P.; Jenwitheesuk, E.; Rongnoparut, P. Homology modeling of mosquito cytochrome P450 enzymes involved in pyrethroid metabolism: Insights into differences in substrate selectivity. BMC Res. Notes 2011, 4, 321. [Google Scholar] [CrossRef]
  78. Liu, N.; Yue, X. Insecticide resistance and cross-resistance in the house fly (Diptera: Muscidae). J. Econ. Entomol. 2000, 93, 1269–1275. [Google Scholar] [CrossRef]
  79. Wittwer, C.T.; Herrmann, M.G.; Moss, A.A.; Rasmussen, R.P. Continuous fluorescence monitoring of rapid cycle DNA amplification. Biotechniques 1997, 22, 130–138. [Google Scholar] [CrossRef]
  80. Livak, K.J.; Schmittgen, T.D. A Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  81. Zhu, F.; Feng, J.-N.; Zhang, L.; Liu, N. Characterization of two novel cytochrome P450 genes in insecticide-resistant house-flies. Insect Mol. Biol. 2008, 17, 27–37. [Google Scholar] [CrossRef]
  82. Gong, Y.; Li, T.; Li, Q.; Liu, S.; Liu, N. The central role of multiple P450 genes and their co-factor CPR in the development of permethrin resistance in the mosquito Culex quinquefasciatus. Front. Physiol. 2022, 12, 802584. [Google Scholar] [CrossRef]
  83. Feng, X.; Liu, N. Functional analyses of house fly carboxylesterases involved in insecticide resistance. Front. Physiol. 2020, 11, 595009. [Google Scholar] [CrossRef] [PubMed]
  84. Roy, A.; Kucukural, A.; Zhang, Y. I-TASSER: A unified platform for automated protein structure and function prediction. Nat. Protoc. 2010, 5, 725–738. [Google Scholar] [CrossRef]
  85. Zhang, Y. I-TASSER server for protein 3D structure prediction. BMC Bioinform. 2008, 9, 40. [Google Scholar] [CrossRef] [PubMed]
  86. Zhang, J.; Liang, Y.; Zhang, Y. Atomic-level protein structure refinement using fragment-guided molecular dynamics conformation sampling. Structure 2011, 19, 1784–1795. [Google Scholar] [CrossRef] [PubMed]
  87. Laskowski, R.A.; MacArthur, M.W.; Moss, D.S.; Thornton, J.M. PROCHECK: A program to check the stereochemical quality of protein structures. J. Appl. Crystallogr. 1993, 26, 283–291. [Google Scholar] [CrossRef]
  88. Sippl, M.J. Recognition of errors in three-dimensional structures of proteins. Proteins 1993, 17, 355–362. [Google Scholar] [CrossRef] [PubMed]
  89. Wiederstein, M.; Sippl, M.J. ProSA-web: Interactive web service for the recognition of errors in three-dimensional structures of proteins. Nucleic Acids Res. 2007, 35, W407–W410. [Google Scholar] [CrossRef]
  90. Chovancova, E.; Pavelka, A.; Benes, P.; Strnad, O.; Brezovsky, J.; Kozlikova, B.; Gora, A.; Sustr, V.; Klvana, M.; Medek, P. CAVER 3.0: A tool for the analysis of transport pathways in dynamic protein structures. PLoS Comput. Biol. 2012, 8, e1002708. [Google Scholar] [CrossRef]
  91. Medek, P.; Beneš, P.; Sochor, J. Computation of tunnels in protein molecules using Delaunay triangulation. J. WSCG 2007, 15, 107–114. [Google Scholar]
  92. Kleywegt, G.J.; Zou, J.Y.; Kjeldgaard, M.; Jones, T.A.; Around, O. International Tables for Crystallography; Kluwer: Dordrecht, The Netherlands, 2001; Volume F, Chapter 17.1; pp. 353–356, 366–367. [Google Scholar]
  93. Morris, G.M.; Huey, R.; Lindstrom, W.; Sanner, M.F.; Belew, R.K.; Goodsell, D.S.; Olson, A.J. AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J. Comput. Chem. 2009, 30, 2785–2791. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  94. Irwin, J.J.; Sterling, T.; Mysinger, M.M.; Bolstad, E.S.; Coleman, R.G. ZINC: A Free Tool to Discover Chemistry for Biology. J. Chem. Inf. Model. 2012, 52, 1757–1768. [Google Scholar] [CrossRef] [PubMed]
  95. Gotoh, O. S Substrate recognition sites in cytochrome P450 family 2 (CYP2) proteins inferred from comparative analyses of amino acid and coding nucleotide sequences. J. Biol. Chem. 1992, 267, 83–89. [Google Scholar] [CrossRef] [PubMed]
  96. Delano, W.L. The PyMOL Molecular Graphics System. 2002. Available online: https://www.pymol.org (accessed on 15 May 2022).
Ijms 24 03170 g001 550
Figure 1. Allele-specific RT-PCR autosomal mapping of Musca domestica P450 genes. PCR fragments were generated using an allele-specific primer set according to the sequence for each gene from ALHF. The absence of a PCR product in a house fly line indicates that the gene was located on the corresponding autosome from aabys (i.e., the absence of a band in the A1234 line indicates that the gene was present on autosome 5).
Figure 1. Allele-specific RT-PCR autosomal mapping of Musca domestica P450 genes. PCR fragments were generated using an allele-specific primer set according to the sequence for each gene from ALHF. The absence of a PCR product in a house fly line indicates that the gene was located on the corresponding autosome from aabys (i.e., the absence of a band in the A1234 line indicates that the gene was present on autosome 5).
Ijms 24 03170 g001
Ijms 24 03170 g002 550
Figure 2. Autosomal interaction in up-regulated P450 gene expression in houseflies. Genes that are located on Autosome 5 are shown in red, on autosome 1 in blue, and on autosome 3 in Yellow. NO P450 genes are identified on either autosome 2 or 4. The autosomes that include factors involved in the interaction and regulation of P450 gene expression are shown with shaded highlights.
Figure 2. Autosomal interaction in up-regulated P450 gene expression in houseflies. Genes that are located on Autosome 5 are shown in red, on autosome 1 in blue, and on autosome 3 in Yellow. NO P450 genes are identified on either autosome 2 or 4. The autosomes that include factors involved in the interaction and regulation of P450 gene expression are shown with shaded highlights.
Ijms 24 03170 g002
Ijms 24 03170 g003 550
Figure 3. Transgenic study of CYP4S24, CYP6A36 and CYP6D10 in D. melanogaster. (A) RT-PCR amplification of the P450 genes in transgenic D. melanogaster lines. “−” represents the non-transgenic empty-vector control D. melanogaster line. “+” represents the transgenic D. melanogaster lines containing house fly P450 genes. * The GelPilot 1Kb (+) ladder was used as a molecular size reference, with the numbers on the figure indicating the DNA band size in bp. (B) Transgenic expression of CYP4S24, CYP6A36 and CYP6D10 in D. melanogaster. The relative expressions of the three transgenes were quantified via qRT-PCR. “P450” represents the homozygous transgenic D. line with house fly P450 genes, and “P450 × GAL4” represents the F1 generation of the homozygous transgenic Drosophila line crossed with GAL4 driver line. Data shown as mean + S.E.M (n = 3). (C) Toxicity of permethrin to non-transgenic and transgenic D. melanogaster lines. Resistance ratios = LD50 of D. melanogaster lines/LD50 of control line. There were no significant differences in expression levels labeled with the same letter (p < 0.05).
Figure 3. Transgenic study of CYP4S24, CYP6A36 and CYP6D10 in D. melanogaster. (A) RT-PCR amplification of the P450 genes in transgenic D. melanogaster lines. “−” represents the non-transgenic empty-vector control D. melanogaster line. “+” represents the transgenic D. melanogaster lines containing house fly P450 genes. * The GelPilot 1Kb (+) ladder was used as a molecular size reference, with the numbers on the figure indicating the DNA band size in bp. (B) Transgenic expression of CYP4S24, CYP6A36 and CYP6D10 in D. melanogaster. The relative expressions of the three transgenes were quantified via qRT-PCR. “P450” represents the homozygous transgenic D. line with house fly P450 genes, and “P450 × GAL4” represents the F1 generation of the homozygous transgenic Drosophila line crossed with GAL4 driver line. Data shown as mean + S.E.M (n = 3). (C) Toxicity of permethrin to non-transgenic and transgenic D. melanogaster lines. Resistance ratios = LD50 of D. melanogaster lines/LD50 of control line. There were no significant differences in expression levels labeled with the same letter (p < 0.05).
Ijms 24 03170 g003
Ijms 24 03170 g004 550
Figure 4. HPLC profiles of permethrin and substrate metabolism mediated by the expressed P450/CPR isomers. Reactions were performed with 20 µM chemical (cis-/trans-permethrin, PBalc or PBald), 1 mM NADPH, 0.25 mM MgCl2 and 0.2 μM P450s/CPR microsomal proteins incubated for 2 h at 30 °C. Each P450/CPR protein with NADPH had a control reaction consisting of the same microsome without NADPH. A negative control with Sf9 cell microsomes (0.2 μM) and without P450/CPR expression was used to calculate substrate depletion. (A) Percentage of trans-permethrin depletion. (B) Percentage of cis-permethrin depletion. (C) Percentage of PBalc depletion. (D) Percentage of PBald depletion. Data shown as mean ± S.E.M (n = 3).
Figure 4. HPLC profiles of permethrin and substrate metabolism mediated by the expressed P450/CPR isomers. Reactions were performed with 20 µM chemical (cis-/trans-permethrin, PBalc or PBald), 1 mM NADPH, 0.25 mM MgCl2 and 0.2 μM P450s/CPR microsomal proteins incubated for 2 h at 30 °C. Each P450/CPR protein with NADPH had a control reaction consisting of the same microsome without NADPH. A negative control with Sf9 cell microsomes (0.2 μM) and without P450/CPR expression was used to calculate substrate depletion. (A) Percentage of trans-permethrin depletion. (B) Percentage of cis-permethrin depletion. (C) Percentage of PBalc depletion. (D) Percentage of PBald depletion. Data shown as mean ± S.E.M (n = 3).
Ijms 24 03170 g004
Ijms 24 03170 g005 550
Figure 5. Solvent channel, family 2 channels and active site of P450s. Channels are shown as connected spheres with different colors: cyan (2a), orange (2ac), purple (2b), blue (2c), wheat (2d), red (2e) and yellow (2f). Helix I and heme of P450s are labeled. (A) CYP6A36, (B) CYP6D10 and (C) CYP4S24. Glutamate 313 and Threonine 317 are also labeled.
Figure 5. Solvent channel, family 2 channels and active site of P450s. Channels are shown as connected spheres with different colors: cyan (2a), orange (2ac), purple (2b), blue (2c), wheat (2d), red (2e) and yellow (2f). Helix I and heme of P450s are labeled. (A) CYP6A36, (B) CYP6D10 and (C) CYP4S24. Glutamate 313 and Threonine 317 are also labeled.
Ijms 24 03170 g005
Ijms 24 03170 g006 550
Figure 6. Permethrin binding models for the active site of P450s. The heme group is represented by green sticks and permethrin by red sticks. (A) Predicted binding model for cis-permethrin in CYP6A36. (B) Predicted binding model for trans-permethrin in CYP6A36. (C) Predicted binding model for cis-permethrin in CYP6D10. (D) Predicted binding model for trans-permethrin in CYP6D10. (E) Predicted binding model for cis-permethrin in CYP4S24. (F) Predicted binding model for trans-permethrin in CYP4S24.
Figure 6. Permethrin binding models for the active site of P450s. The heme group is represented by green sticks and permethrin by red sticks. (A) Predicted binding model for cis-permethrin in CYP6A36. (B) Predicted binding model for trans-permethrin in CYP6A36. (C) Predicted binding model for cis-permethrin in CYP6D10. (D) Predicted binding model for trans-permethrin in CYP6D10. (E) Predicted binding model for cis-permethrin in CYP4S24. (F) Predicted binding model for trans-permethrin in CYP4S24.
Ijms 24 03170 g006
Ijms 24 03170 g007 550
Figure 7. PBalc and PBald binding models for the active site of P450s. The heme group is represented by green sticks, PBalc by yellow sticks and PBald by pink sticks. (A) Predicted binding model for PBalc in CYP6A36. (B) Predicted binding model for PBalc in CYP6D10. (C) Predicted binding model for PBalc in CYP4S24. (D) Predicted binding model for PBald in CYP6A36. (E) Predicted binding model for PBald in CYP6D10. (F) Predicted binding model for PBald in CYP4S24.
Figure 7. PBalc and PBald binding models for the active site of P450s. The heme group is represented by green sticks, PBalc by yellow sticks and PBald by pink sticks. (A) Predicted binding model for PBalc in CYP6A36. (B) Predicted binding model for PBalc in CYP6D10. (C) Predicted binding model for PBalc in CYP4S24. (D) Predicted binding model for PBald in CYP6A36. (E) Predicted binding model for PBald in CYP6D10. (F) Predicted binding model for PBald in CYP4S24.
Ijms 24 03170 g007
Table
Table 1. Clan, family and number distribution of selected cytochrome P450 genes in ALHF Musca domestica.
Table 1. Clan, family and number distribution of selected cytochrome P450 genes in ALHF Musca domestica.
CYP ClanCYP Family *Number of CYP Members
Clan 2CYP181
CYP3041
CYP3051
CYP3061
Clan 3CYP627
CYP96
CYP283
CYP3101
CYP3171
CYP438A41
Clan 4CYP423
CYP3111
CYP30733
CYP3132
CYP3181
Mitochondrial CYPCYP129
CYP3011
CYP3021
CYP3141
CYP3151
* The P450 gene sequence information was generated from the RNA-seq (NCBI:SRR521288, NCBI:SRR521289) and genome (PRJNA210139, PRJNA176013) analysis of Musca domestica.
Table
Table 2. Relative expression levels and the predicted autosomal interactions of the up-regulated P450 genes among ALHF strain and other lines of M. domestica.
Table 2. Relative expression levels and the predicted autosomal interactions of the up-regulated P450 genes among ALHF strain and other lines of M. domestica.
Relative Gene Expression in Strains/Lines §
GeneALHFA2345A1345A1245A1235A1234
CYP4E10274.42 ± 22.40125.67 ± 10.39 *144.62 ± 12.62 *257.81 ± 23.61268.61 ± 25.866.6 ± 6.74 *
CYP4E115.86 ± 0.556.31 ± 0.482.61 ± 0.17 *5.99 ± 0.715.80 ± 0.362.71 ± 0.26 *
CYP4G132.13 ± 0.131.04 ± 0.10 *1.22 ± 0.09 *1.04 ± 0.04 *2.01 ± 0.121.91 ± 0.08
CYP4G994.52 ± 0.312.12 ± 0.22 *4.28 ± 0.411.38 ± 0.13 *4.43 ± 0.454.43 ± 0.23
CYP4S242.94 ± 0.081.03 ± 0.05 *1.14 ± 0.06 *0.83 ± 0.03 *2.83 ± 0.493.04 ± 0.53
CYP6A367.03 ± 0.642.60 ± 0.21 *2.80 ± 0.41 *6.63 ± 0.627.01 ± 1.041.72 ± 0.19 *
CYP6A402.95 ± 0.222.90 ± 0.321.16 ± 0.10 *2.99 ± 0.412.68 ± 0.361.17 ± 0.05 *
CYP6A522.74 ± 0.152.50 ± 0.222.70 ± 0.172.54 ± 0.083.01 ± 0.240.81 ± 0.03 *
CYP6A584.32 ± 1.021.32 ± 0.23 *1.86 ± 0.36 *5.63 ± 2.174.19 ± 0.220.64 ± 0.15 *
CYP6D32.40 ± 0.131.38 ± 0.07 *1.34 ± 0.06 *2.25 ± 0.132.40 ± 0.372.34 ± 0.18
CYP6D106.24 ± 0.632.55 ± 0.37 *3.06 ± 0.55 *5.99 ± 0.256.19 ± 0.512.62 ± 0.32 *
§ The relative levels of gene expression are shown as a ratio in comparison with those in aabys flies; the data are shown as the mean ± SEM. * Gene expression value within a given M. domestica autosomal line was significantly lower than the expression in the parental ALHF strain at a p < 0.05 level of significance.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, M.; Feng, X.; Reid, W.R.; Tang, F.; Liu, N. Multiple-P450 Gene Co-Up-Regulation in the Development of Permethrin Resistance in the House Fly, Musca domestica. Int. J. Mol. Sci. 2023, 24, 3170. https://doi.org/10.3390/ijms24043170

AMA Style

Li M, Feng X, Reid WR, Tang F, Liu N. Multiple-P450 Gene Co-Up-Regulation in the Development of Permethrin Resistance in the House Fly, Musca domestica. International Journal of Molecular Sciences. 2023; 24(4):3170. https://doi.org/10.3390/ijms24043170

Chicago/Turabian Style

Li, Ming, Xuechun Feng, William R. Reid, Fang Tang, and Nannan Liu. 2023. "Multiple-P450 Gene Co-Up-Regulation in the Development of Permethrin Resistance in the House Fly, Musca domestica" International Journal of Molecular Sciences 24, no. 4: 3170. https://doi.org/10.3390/ijms24043170

APA Style

Li, M., Feng, X., Reid, W. R., Tang, F., & Liu, N. (2023). Multiple-P450 Gene Co-Up-Regulation in the Development of Permethrin Resistance in the House Fly, Musca domestica. International Journal of Molecular Sciences, 24(4), 3170. https://doi.org/10.3390/ijms24043170

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop