Frequencies of pyrethroid resistance-associated mutations of Vssc1 and CYP6D1 in field populations of Musca domestica L. in Turkey

House flies were collected from 16 different provinces in the Aegean and Mediterranean regions of Turkey, and the frequencies of pyrethroid resistance-associated mutations in Vssc1 and CYP6D1 in these field-collected populations were studied. Although there is no organized resistance management program for house fly control in Turkey, it is known that different groups of insecticides, including pyrethroids, are used. The frequencies of both Vssc1 and CYP6D1 alleles were weighted toward the susceptibles, with Vssc1-susceptible alleles having higher frequencies in both regions (0.75 in Aegean and 0.69 in Mediterranean populations) than CYP6D1-susceptible alleles (0.65 in Aegean and 0.56 in Mediterranean populations). The frequencies of kdr-his alleles were higher than the frequencies of kdr alleles in these populations. While the frequencies of kdr-his alleles were close to each other in the Aegean (0.23) and Mediterranean (0.17) populations, the frequencies of kdr alleles remarkably differed in these two regions, with values of 0.02 and 0.14, respectively. In contrast to Europe, Asia, and the U.S.A., no super-kdr allele was detected in the samples from both regions. We identified six and eight different Vssc1+CYP6D1 genotype classes in the Aegean and Mediterranean regions, respectively. The three most common genotype classes in the regions were susceptible Vssc1 with heterozygous CYP6D1v1 (29%), sus/kdr-his1 with heterozygous CYP6D1v1 (23%), and susceptible Vssc1 with CYP6D1 (22%). The total frequencies of these three most common genotype classes (approximately 75%) obtained in our study were very close to the value obtained in Florida in a previous study, which was related by the similarity of temperature patterns between Florida and the corresponding regions of Turkey. This may reflect the lack of overwintering fitness cost associated with resistance alleles in both climates.     

Assessing permethrin resistance in the stable fly (Diptera: Muscidae) in Florida by using laboratory selections and field evaluations

Insecticide resistance in the stable fly, Stomoxys calcitrans (L.) (Diptera: Muscidae),has been demonstrated previously, but mostly with insecticides that are no longer used, such as the organochlorines. Resistance to commonly used pyrethroids has been evaluated twice, but only in the midwestern United States. Stable fly susceptibility to a commonly used pyrethroid, permethrin, was determined in Florida to assess the possibility of resistance development. Diagnostic concentration evaluations of three stable fly field strains demonstrated a maximum of 57 and 21% survival to permethrin residues of 3x and 10x the LC99 of a susceptible strain, respectively. Stable flies from an equine facility with no reported insecticide use demonstrated approximately 20% survival with a 3x diagnostic concentration. Despite a distance of 91-km between field collection sites, survival profiles of field-collected stable fly strains were similar. Although an established stable fly colony collected from a local dairy previously expressed low level resistance to permethrin residues, five generations of laboratory permethrin selection increased resistance 15-fold.     

A cis-regulatory sequence driving metabolic insecticide resistance in mosquitoes: Functional characterisation and signatures of selection

Although cytochrome P450 (CYP450) enzymes are frequently up-regulated in mosquitoes resistant to insecticides, no regulatory motifs driving these expression differences with relevance to wild populations have been identified. Transposable elements (TEs) are often enriched upstream of those CYP450s involved in insecticide resistance, leading to the assumption that they contribute regulatory motifs that directly underlie the resistance phenotype. A partial CuRE1 (Culex Repetitive Element 1) transposable element is found directly upstream of CYP9M10, a cytochrome P450 implicated previously in larval resistance to permethrin in the ISOP450 strain of Culex quinquefasciatus, but is absent from the equivalent genomic region of a susceptible strain. Via expression of CYP9M10 in Escherichia coli we have now demonstrated time- and NADPH-dependant permethrin metabolism, prerequisites for confirmation of a role in metabolic resistance, and through qPCR shown that CYP9M10 is >20-fold over-expressed in ISOP450 compared to a susceptible strain. In a fluorescent reporter assay the region upstream of CYP9M10 from ISOP450 drove 10× expression compared to the equivalent region (lacking CuRE1) from the susceptible strain. Close correspondence with the gene expression fold-change implicates the upstream region including CuRE1 as a cis-regulatory element involved in resistance. Only a single CuRE1 bearing allele, identical to the CuRE1 bearing allele in the resistant strain, is found throughout Sub-Saharan Africa, in contrast to the diversity encountered in non-CuRE1 alleles. This suggests a single origin and subsequent spread due to selective advantage. CuRE1 is detectable using a simple diagnostic. When applied to C. quinquefasciatus larvae from Ghana we have demonstrated a significant association with permethrin resistance in multiple field sites (mean Odds Ratio?=?3.86) suggesting this marker has relevance to natural populations of vector mosquitoes. However, when CuRE1 was excised from the allele used in the reporter assay through fusion PCR, expression was unaffected, indicating that the TE has no direct role in resistance and hence that CuRE1 is acting only as a marker of an as yet unidentified regulatory motif in the association analysis. This suggests that a re-evaluation of the assumption that TEs contribute regulatory motifs involved in gene expression may be necessary.     

Identification of a mutation associated with permethrin resistance in the para-type sodium channel of the stable fly (Diptera: Muscidae)

The insect sodium channel is of particular interest for evaluating resistance to pyrethroids because it is the target molecule for this major class of neurotoxic insecticides. The stable fly, Stomoxys calcitrans (L.) (Diptera: Muscidae), sodium channel coding sequence representing domains IS6 through IVS6 was isolated, and the sequence encoding domain II was compared among individuals of a laboratory strain selected for resistance to permethrin and the unselected, parental generation. A point mutation resulting in a leucine-to-histidine amino acid change was identified (Leul014His), and its location corresponded with that observed for knockdown resistance (kdr) mutations in other insects. As a result, the allele was designated kdr-his. A molecular assay was developed to assess the frequency of this mutation in genomic DNA of individual stable flies from the laboratory selections, which provided further evidence that the kdr-his allele accounts for the observed level ofpermethrin resistance in the selected strain. The assay was then used to evaluate the frequency of the mutation from five field-collected populations originating from three horse farms near Ocala, FL; one horse farm near Gainesville, FL; and one dairy farm near Hague, FL. Frequency of the kdr-his allele ranged from 0.46 to 0.78, supporting further investigation of allele prevalence throughout the stable fly season and in response to field insecticide application.     

Co-up-regulation of three P450 genes in response to permethrin exposure in permethrin resistant house flies, Musca domestica

Background

Insects may use various biochemical pathways to enable them to tolerate the lethal action of insecticides. For example, increased cytochrome P450 detoxification is known to play an important role in many insect species. Both constitutively increased expression (overexpression) and induction of P450s are thought to be responsible for increased levels of detoxification of insecticides. However, unlike constitutively overexpressed P450 genes, whose expression association with insecticide resistance has been extensively studied, the induction of P450s is less well characterized in insecticide resistance. The current study focuses on the characterization of individual P450 genes that are induced in response to permethrin treatment in permethrin resistant house flies.

Results

The expression of 3 P450 genes, CYP4D4v2, CYP4G2, and CYP6A38, was co-up-regulated by permethrin treatment in permethrin resistant ALHF house flies in a time and dose-dependent manner. Comparison of the deduced protein sequences of these three P450s from resistant ALHF and susceptible aabys and CS house flies revealed identical protein sequences. Genetic linkage analysis located CYP4D4v2 and CYP6A38 on autosome 5, corresponding to the linkage of P450-mediated resistance in ALHF, whereas CYP4G2 was located on autosome 3, where the major insecticide resistance factor(s) for ALHF had been mapped but no P450 genes reported prior to this study.

Conclusion

Our study provides the first direct evidence that multiple P450 genes are co-up-regulated in permethrin resistant house flies through the induction mechanism, which increases overall expression levels of P450 genes in resistant house flies. Taken together with the significant induction of CYP4D4v2, CYP4G2, and CYP6A38 expression by permethrin only in permethrin resistant house flies and the correlation of the linkage of the genes with resistance and/or P450-mediated resistance in resistant ALHF house flies, this study sheds new light on the functional importance of P450 genes in response to insecticide treatment, detoxification of insecticides, the adaptation of insects to their environment, and the evolution of insecticide resistance.     

Differential expression of CYP6A5 and CYP6A5v2 in pyrethroid-resistant house flies, Musca domestica

Two cytochrome P450 alleles, CYP6A5 and CYP6A5v2, were isolated from a pyrethroid-resistant house fly stain, ALHF. The two alleles shared 98% similarity in amino acid sequence. To understand the importance of these two alleles in resistance and examine the expression profile of the two alleles between resistant and susceptible strains, quantitative real-time PCR (qRT-PCR) was performed and compared with the Northern blot analysis. We found that qRT-PCR was an efficient method to characterize the expression profiles between these two sequence-closely-related P450 genes between resistant and susceptible houses flies. One of them, CYP6A5v2, was constitutively overexpressed in ALHF house flies compared with susceptible house fly strains. Moreover, this gene was predominantly expressed in the abdominal tissues of ALHF, in which the primary detoxification organs of insects are located. However, there was no significant difference in the expression of CYP6A5 between ALHF and susceptible house flies. The genetic linkage analysis was conducted to determine the possible link between the constitutively overexpressed CYP6A5v2 and insecticide resistance. CYP6A5v2 was mapped on autosome 5, which is correlated with the linkage of resistance in ALHF. Taken together, the study suggests the importance of CYP6A5v2 in increasing metabolic detoxification of insecticides in ALHF. The distinct expression of CYP6A5 and CYP6A5v2 in resistant and susceptible house flies implies the functional difference of theses two genes in house flies and suggests that they are two recently diverged P450 genes presented in a single organism.     

Susceptibility of house flies (Diptera: Muscidae) and five pupal parasitoids (Hymenoptera: Pteromalidae) to abamectin and seven commercial insecticides.

Assays of five commercial insecticides applied as residual sprays at label rates to plywood indicated the most toxic insecticide overall for pteromalid parasitoids of house flies, Musca domestica L., was Atroban (permethrin), followed by Ciodrin (crotoxyphos), Rabon (tetrachlorvinphos), Ectrin (fenvalerate), and Cygon (dimethoate). Insecticide-susceptible house flies were susceptible to all five insecticides (mortality, 62-100%). Flies that were recently colonized from populations on dairy farms in New York were susceptible only to Rabon. Urolepis rufipes (Ashmead) was the most susceptible parasitoid species overall to these insecticides, followed by Muscidifurax raptor Girault & Sanders, Nasonia vitripennis Walker, Pachycrepoideus vindemmiae (Rondani), and Spalangia cameroni Perkins. Compared with susceptible flies, newly colonized flies showed moderate resistance to avermectin B1a (abamectin). Abamectin was more toxic to all of the parasitoids except N. vitripennis and S. cameroni than to newly colonized house flies when exposed for 90 min to plywood boards treated with 0.001-0.1% abamectin. Space sprays with Vapona (dichlorvos) killed all of the parasitoids and susceptible flies and 64% of the newly colonized flies when insects were placed directly in the path of the spray; mortality was substantially lower among flies and parasitoids protected under 5 cm of wheat straw. Space sprays with Pyrenone (pyrethrins) killed greater than 86% of all insects exposed to the spray path except for the newly colonized flies (1% mortality); mortality of insects protected under straw was low (less than 12%) except for S. cameroni (76%). Because responses of the five parasitoids to the different insecticides varied considerably, general conclusions about parasitoid susceptibility to active ingredients, insecticide class, or method of application were not possible.     

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.     

The effect of three environmental conditions on the fitness of cytochrome P450 monooxygenase-mediated permethrin resistance in Culex pipiens quinquefasciatus

Background  

The evolution of insecticide resistance and persistence of resistance phenotypes are influenced by the fitness of resistance alleles in the absence of insecticide pressure. Experimental determination of fitness is difficult, but fitness can be inferred by measuring changes in allele frequencies in appropriate environments. We conducted allele competition experiments by crossing two highly related strains of Culex pipiens quinquefasciatus mosquitoes. One strain (ISOP450) was permethrin resistant (due to P450-mediated detoxification) and one was a susceptible strain. Allele and genotype frequencies were examined for 12 generations under three environmental conditions: standard laboratory, temephos exposure (an insecticide to which the P450 detoxification mechanism in ISOP450 confers no resistance and which is commonly used in mosquito control programs) and cold temperature stress (mimics the colder temperatures within the habitat of this mosquito).     

High-resolution melt analysis for the detection of a mutation associated with permethrin resistance in a population of scabies mites

Permethrin as a topical acaricide cream is widely used to treat scabies. The neuronal voltage-sensitive sodium channel (Vssc), necessary for the generation of action potentials in excitable cells, is the target of pyrethroid acaricides such as permethrin. Pyrethroid resistance has been linked to specific mutations in the Vssc gene. Following the partial sequencing of the Vssc gene in the scabies mite Sarcoptes scabiei (L.) (Astigmata: Sarcoptidae), we compared Vssc gene sequences from permethrin-sensitive and -tolerant S. scabiei var. canis Gerlach mites, and identified a G to A single nucleotide polymorphism (SNP) in permethrin-tolerant mites resulting in an amino acid change from glycine to aspartic acid in domain III S6. The mutation is in a region of the gene where mutations have been identified in a range of pyrethroid-resistant arthropods. Results of in vitro permethrin exposure assays showed that survival rates for mites bearing the mutation were similar to those previously reported for mites from human subjects where clinical tolerance to permethrin had been observed. A real-time polymerase chain reaction-high-resolution melt (PCR-HRM) assay was developed to detect this SNP. This assay provides a useful methodology for screening for this and other mutations associated with permethrin resistance in scabies mite populations and thus facilitates surveillance for acaricide resistance.     

Overwintering cost associated with resistance to transgenic cotton in the pink bollworm (Lepidoptera: Gelechiidae)

Fitness costs associated with resistance to transgenic crops producing toxins from Bacillus thuringiensis (Bt) may have important effects on the evolution of resistance. We investigated overwintering costs in pink bollworm, Pectinophora gosypiella (Saunders), strains with different degrees of resistance to Bt cotton. Frequency of resistant individuals in a strain was not associated with induction of diapause or emergence from diapause in early winter. Emergence from diapause in the spring was 71% lower in three highly resistant strains than in two heterogeneous strains from which the resistant strains were derived. This underestimates the overwintering cost because the frequency of the resistance allele was relatively high in the heterogeneous strains. Emergence in the spring in hybrid progeny from crosses between the resistant and heterogeneous strains was greater than in resistant strains but did not differ from susceptible strains, showing that the overwintering cost was recessive to some extent.     

Three mutations identified in the voltage-sensitive sodium channel alpha-subunit gene of permethrin-resistant human head lice reduce the permethrin sensitivity of house fly Vssc1 sodium channels expressed in Xenopus oocytes

Point mutations in the para-orthologous sodium channel alpha-subunit of the head louse (M815I, T917I, and L920F) are associated with permethrin resistance and DDT resistance. These mutations were inserted in all combinations using site-directed mutagenesis at the corresponding amino acid sequence positions (M827I, T929I, and L932F) of the house fly para-orthologous voltage-sensitive sodium channel alpha-subunit (Vssc1(WT)) gene and heterologously co-expressed with the sodium channel auxiliary subunit of house fly (Vsscbeta) in Xenopus oocytes. The double mutant possessing M827I and T929I (Vssc1(MITI)/Vsscbeta) caused a approximately 4.0mV hyperpolarizing shift and the triple mutant, Vssc1(MITILF)/Vsscbeta, caused a approximately 3.2mV depolarizing shift in the voltage dependence of activation curves. Vssc1(MITI)/Vsscbeta, Vssc1(TILF)/Vsscbeta, and Vssc1(MITILF)/Vsscbeta caused depolarizing shifts ( approximately 6.6, approximately 7.6, and approximately 8.8mV, respectively) in the voltage dependence of steady-state inactivation curves. The M827I and L932F mutations reduced permethrin sensitivity when expressed alone but the T929I mutation, either alone or in combination, virtually abolished permethrin sensitivity. Thus, the T929I mutation is the principal cause of permethrin resistance in head lice. Comparison of the expression rates of channels containing single, double and triple mutations with that of Vssc1(WT)/Vsscbeta channels indicates that the M827I mutation may play a role in rescuing the decreased expression of channels containing T929I.     

Allelic Variation of Cytochrome P450s Drives Resistance to Bednet Insecticides in a Major Malaria Vector

Scale up of Long Lasting Insecticide Nets (LLINs) has massively contributed to reduce malaria mortality across Africa. However, resistance to pyrethroid insecticides in malaria vectors threatens its continued effectiveness. Deciphering the detailed molecular basis of such resistance and designing diagnostic tools is critical to implement suitable resistance management strategies. Here, we demonstrated that allelic variation in two cytochrome P450 genes is the most important driver of pyrethroid resistance in the major African malaria vector Anopheles funestus and detected key mutations controlling this resistance. An Africa-wide polymorphism analysis of the duplicated genes CYP6P9a and CYP6P9b revealed that both genes are directionally selected with alleles segregating according to resistance phenotypes. Modelling and docking simulations predicted that resistant alleles were better metabolizers of pyrethroids than susceptible alleles. Metabolism assays performed with recombinant enzymes of various alleles confirmed that alleles from resistant mosquitoes had significantly higher activities toward pyrethroids. Additionally, transgenic expression in Drosophila showed that flies expressing resistant alleles of both genes were significantly more resistant to pyrethroids compared with those expressing the susceptible alleles, indicating that allelic variation is the key resistance mechanism. Furthermore, site-directed mutagenesis and functional analyses demonstrated that three amino acid changes (Val¹⁰⁹Ile, Asp³³⁵Glu and Asn³⁸⁴Ser) from the resistant allele of CYP6P9b were key pyrethroid resistance mutations inducing high metabolic efficiency. The detection of these first DNA markers of metabolic resistance to pyrethroids allows the design of DNA-based diagnostic tools to detect and track resistance associated with bednets scale up, which will improve the design of evidence-based resistance management strategies.     

Fitness consequences of malathion-specific resistance in red flour beetle (Coleoptera: Tenebrionidae) and selection for resistance in the absence of malathion

Malathion resistance in the red flour beetle, Tribolium castaneum (Herbst), is a worldwide problem and is very stable once it becomes widespread in natural populations. In the absence of insecticide the proportion of resistant phenotypes may rapidly decline but the development of resistance does not always involve reduced fitness. Malathion-specific resistance in T. castaneum seems not to involve any loss of fitness in laboratory or field conditions. Susceptible beetles were in competition with resistant beetles at different initial frequencies and modifications of susceptible gene frequency were estimated in these laboratory populations over 10 generations. A significant decrease in susceptible gene frequency was observed in Tribolium populations over time. The selection coefficient of the susceptible allele was estimated and the fitness of susceptible alleles in all tests was observed to range from 0.89 to 0.93 compared with the fitness of resistant genotypes, which was assumed to be 1. Data provided evidence that the resistant strains exhibited fitness advantages in the absence of malathion. We also compared the biotic potential (fecundity and developmental time) of the susceptible strain, the homozygous malathion-specific resistant strain, and their hybrids. Malathion-specific resistant strains showed an 8 -23% increase in biotic potential relative to the susceptible strain. These findings are consistent with those of malathion-specific resistance in T. castaneum; the fitness of the insects seems independent of the genetic background and the fitness of the resistant insects is not affected by this resistance mechanism.     

Analysis of Clines with Variable Selection and Variable Migration

We report a likelihood-based method that estimates both dispersal and natural selection using the rate of change of the shape of a cline when selection and migration are not constant through time. We have investigated the case of local adaptation of the mosquito Culex pipiens to organophosphate insecticides in the Montpellier area in France. We have analyzed the modification of the clinal patterns at two resistance loci during the period from breeding to overwintering. We show that mosquitoes migrate extensively from breeding to overwintering sites at a rate that is markedly different from previous estimates made during the breeding season only. This migration is also strongly asymmetrical, which can be explained by different geographical distributions of breeding and overwintering sites, by variation in mosquito density along the transect, or by behavioral biases. We found that the starting time of overwintering is likely to vary between northern and southern populations and that substantial fitness costs are associated with resistance alleles at the two loci during overwintering. These results illustrate how demography and adaptive microevolution can be studied using selected markers. The method provides a framework to use population genetics and statistical models to reveal ecological and evolutionary processes.     

Influence of permethrin, diazinon and ivermectin treatments on insecticide resistance in the horn fly (Diptera: Muscidae)

The history of insecticide resistance in the horn fly, Haematobia irritans, and the relationship between the characteristics of horn fly biology and insecticide use on resistance development is discussed. Colonies of susceptible horn flies were selected for resistance with six insecticide treatment regimens: continuous single use of permethrin, diazinon and ivermectin: permethrin-diazinon (1:2) mixture; and permethrin-diazinon and permethrin-ivermectin rotation (4-month cycle). Under laboratory conditions, resistance developed during generations 21, 31 and 30 to permethrin, diazinon and ivermectin, respectively. The magnitude of resistance ranged from < 3-fold with ivermectin to 1470-fold with permethrin. Field studies demonstrated that use of a single class of insecticidal ear tag during the horn-fly season resulted in product failure within 3-4 years for pyrethroids and organophosphates, respectively. In laboratory studies, use of alternating insecticides or a mixture of insecticides delayed the onset of resistance for up to 12 generations and reduced the magnitude of pyrethroid resistance. In field studies, yearly alternated use of pyrethroids and organophosphates did not slow or reverse pyrethroid resistance (Barros et al., unpublished data), while a 2-year alternated use with organophosphates resulted in partial reversion of pyrethroid resistance. When pyrethroid and organophosphate ear tags were used in a mosaic strategy at two different locations, efficacy of products did not change during a 3-year period.     

Shifting clinal patterns of stress resistance traits in Drosophila ananassae

Drosophila ananassae, a desiccation and cold sensitive species, is abundant along the latitudinal gradient of the Indian subcontinent. Analysis of seasonally varying wild-caught flies showed two independent patterns of melanisation: (1) narrow to broad melanic stripes on three anterior abdominal segments only; (2) a novel body color pattern (dark vs. light background). We investigated the degree to which these two melanisation systems vary; first with latitude and secondly among seasons. There is a shallow latitudinal cline for percent striped melanisation as well as for frequency of body color alleles during the rainy season. The frequencies of body color alleles vary significantly across seasons in the northern populations i.e. the light allele occur abundantly (>0.94) during the rainy season while the frequency of the dark allele increases (0.22–0.35) during the dry season causing steeper clines during the dry season. By contrast, the low variations in abdominal stripes showed non-significant changes and the cline was similar across seasons. Furthermore, both types of melanisation patterns showed no plasticity with respect to temperature. The present study also investigated clines related to desiccation, heat and cold stress in D. ananassae females across seasons (rainy and dry) from nine latitudinal populations. The clines for stress related traits changes to steeper and non-linear during the dry season. Thus, latitudinal populations of D. ananassae differ in slope values of clines for stress related traits across seasons. This study reports seasonal changes in latitudinal clines of stress resistance traits as seen in a changing frequency of body color alleles of D. ananassae in northern locality, while in southern localities it remains constant. This is presumably the result of only minor seasonal changes in humidity and temperature in the South, whereas in the North seasonal climatic variability is much higher.