Genotype to phenotype, the molecular and physiological dimensions of resistance in arthropods

The dynamic changes of genes revealed that persistently overexpressed genes drive the evolution of cyflumetofen resistance in Tetranychus cinnabarinus

Changes in gene expression are associated with the evolution of pesticide resistance in arthropods. In this study, transcriptome sequencing was performed in 3 different resistance levels (low, L; medium, M; and high, H) of cyflumetofen-resistant strain (YN-CyR). A total of 1 685 genes, including 97 detoxification enzyme genes, were upregulated in all 3 stages, of which 192 genes, including 11 detoxification enzyme genes, showed a continuous increase in expression level with resistance development (L to H). RNA interference experiments showed that overexpression of 7 genes (CYP392A1, TcGSTd05, CCE06, CYP389A1, TcGSTz01, CCE59, and CYP389C2) is involved in the development of cyflumetofen resistance in Tetranychus cinnabarinus. The recombinant CYP392A1 can effectively metabolize cyflumetofen, while CCE06 can bind and sequester cyflumetofen in vitro. We compared 2 methods for rapid screening of resistance molecular markers, including short-term induction and 1-time high-dose selection. Two detoxification enzyme genes were upregulated in the field susceptible strain (YN-S) by induction with 20% lethal concentration (LC₂₀ ) of cyflumetofen. However, 16 detoxification enzyme genes were upregulated by 1-time selection with LC₈₀ of cyflumetofen. Interestingly, the 16 genes were overexpressed in all 3 resistance stages. These results indicated that 1 685 genes that were upregulated at the L stage constituted the basis of cyflumetofen resistance, of which 192 genes in which upregulation continued to increase were the main driving force for the development of resistance. Moreover, the 1-time high-dose selection is an efficient way to rapidly obtain the resistance-related genes that can aid in the development of resistance markers and resistance management in mites.

Trans-driven variation in expression is common among detoxification genes in the extreme generalist herbivore Tetranychus urticae

The extreme adaptation potential of the generalist herbivore Tetranychus urticae (the two-spotted spider mite) to pesticides as well as diverse host plants has been associated with clade-specific gene expansions in known detoxifying enzyme families, and with extensive and rapid transcriptional responses. However, how this broad transcriptional potential is regulated remains largely unknown. Using a parental/F1 design in which four inbred strains were crossed to a common inbred strain, we assessed the genetic basis and inheritance of gene expression variation in T. urticae. Mirroring known phenotypic variation in the progenitor strains of the inbreds, we confirmed that the inbred strains we created were genetically distinct, varied markedly in pesticide resistance, and also captured variation in host plant fitness as is commonly observed in this species. By examining differences in gene expression between parents and allele-specific expression in F1s, we found that variation in RNA abundance was more often explained in trans as compared to cis, with the former associated with dominance in inheritance. Strikingly, in a gene ontology analysis, detoxification genes of the cytochrome P450 monooxygenase (CYP) family, as well as dioxygenases (DOGs) acquired from horizontal gene transfer from fungi, were specifically enriched at the extremes of trans-driven up- and downregulation. In particular, multiple CYPs and DOGs with broad substrate-specificities for pesticides or plant specialized compounds were exceptionally highly upregulated as a result of trans-regulatory variation, or in some cases synergism of cis and trans, in the most multi-pesticide resistant strains. Collectively, our findings highlight the potential importance of trans-driven expression variation in genes associated with xenobiotic metabolism and host plant use for rapid adaptation in T. urticae, and also suggests modular control of these genes, a regulatory architecture that might ameliorate negative pleiotropic effects.

Long-term survey and characterization of cyflumetofen resistance in Tetranychus urticae populations from Turkey

The two-spotted spider mite, Tetranychus urticae Koch (Acari: Tetranychidae) is the most economically important mite pest in agricultural areas and chemical acaricides are widely used to control T. urticae populations. Cyflumetofen is a recently introduced acaricide that inhibits the mitochondrial electron transport chain at complex II (succinate dehydrogenase, SDH), which represents the most recently developed mode of action for mite control worldwide. In the present study, started upon the launch of cyflumetofen in Turkey, a five-year survey was performed to monitor cyflumetofen susceptibility in 28 T. urticae populations collected from agricultural fields across the country. The first resistance case that might cause control failure in practical field conditions was uncovered in 2019, three years after the registration of cyflumetofen. In addition, an extremely resistant population (1722-fold resistance) was also detected towards the end of 2019. Cyflumetofen resistance did not decrease in the laboratory after relaxation of selection pressure for over one year in field-collected populations, suggesting the absence of a fitness cost associated with resistance in these populations. Next to phenotypic resistance, metabolic and physiological mechanisms underlying the decreased susceptibility were also investigated. Synergism assays showed the involvement of P450 monooxygenases in cyflumetofen resistance. Downregulation of carboxylesterases as resistance mechanism, is underpinned by the fact that pre-treatment with esterase inhibitor DEF decreased cyflumetofen toxicity in field-collected strains. Furthermore, a novel H258L substitution in the subunit B of complex II was uncovered in a field population. In silico modeling of the new mutation suggested that the mutation might indeed influence toxicity to complex II inhibitors cyenopyrafen and pyflubumide, but most likely not cyflumetofen. However, further studies are needed to uncover the exact role of this mutation in resistance to this new class of complex II inhibitors.

Polygenic adaptation contributes to the invasive success of the Colorado potato beetle

How invasive species cope with novel selective pressures with limited genetic variation is a fundamental question in molecular ecology. Several mechanisms have been proposed, but they can lack generality. Here, we addressed an alternative solution, polygenic adaptation, wherein traits that arise from multiple combinations of loci may be less sensitive to loss of variation during invasion. We tested the polygenic signal of environmental adaptation of Colorado potato beetle (CPB) introduced in Eurasia. Population genomic analyses showed declining genetic diversity in the eastward expansion of Eurasian populations, and weak population genetic structure (except for the invasion fronts in Asia). Demographic history showed that all populations shared a strong bottleneck about 100 years ago when CPB was introduced to Europe. Genome scans revealed a suite of genes involved in activity regulation functions that are plausibly related to cold stress, including some well-founded functions (e.g., the activity of phosphodiesterase, the G-protein regulator) and discrete functions. Such polygenic architecture supports the hypothesis that polygenic adaptation and potentially genetic redundancy can fuel the adaptation of CPB despite strong genetic depletion, thus representing a promising general mechanism for resolving the genetic paradox of invasion. More broadly, most complex traits based on polygenes may be less sensitive to invasive bottlenecks, thus ensuring the evolutionary success of invasive species in novel environments.

Utility and challenges of using whole-genome resequencing to detect emerging insect and mite resistance in agroecosystems

Arthropods that invade agricultural ecosystems systematically evolve resistance to the control measures used against them, and this remains a significant and ongoing challenge for sustainable food production systems. Early detection of resistance evolution could prompt remedial action to slow the spread of resistance alleles in the landscape. Historical approaches used to detect emerging resistance included phenotypic monitoring of agricultural pest populations, as well as monitoring of allele frequency changes at one or a few candidate pesticide resistance genes. In this article, I discuss the successes and limitations of these traditional monitoring approaches and then consider whether whole-genome scanning could be applied to samples collected from agroecosystems over time for resistance monitoring. I examine the qualities of agroecosystems that could impact application of this approach to pesticide resistance monitoring and describe a recent retrospective analysis where genome scanning successfully detected an oligogenic response to selection by pesticides years prior to pest management failure. I conclude by considering areas of further study that will shed light on the feasibility of applying whole-genome scanning for resistance risk monitoring in agricultural pest species.

Acaricidal Activity and Field Efficacy Analysis of the Potential Biocontrol Agent Bacillus vallismortis NBIF-001 against Spider Mites

In recent years, spider mites have caused considerable economic losses to global agriculture. However, currently available management strategies are limited because of the rapid development of resistance. In this study, Bacillus vallismortis NBIF-001 was isolated and evaluated for its acaricidal activity. NBIF-001 exhibited a significant lethal effect on spider mites within 48 h. The median lethal concentration (LC₅₀) of the culture powders (3.2 × 10¹⁰ CFU/g) was 50.2 µg/mL for Tetranychus urticae (red form), 18.0 µg/mL for T. urticae (green form), and 15.7 µg/mL for Panonychus citri (McGregor). Cultivation optimisation experiments showed that when the number of spores increased, fermentation toxicity also increased. Moreover, field experiments demonstrated that NBIF-001 performed well in the biocontrol of P. citri, which showed a similar corrected field efficacy with the chemical control (67.1 ± 7.9% and 71.1 ± 6.4% after 14 days). Genomics analysis showed that NBIF-001 contains 231 factors and seven gene clusters of metabolites that may be involved in its acaricidal activity. Further bioassays of the fermentation supernatants showed that 50× dilution treatments killed 72.5 ± 5.4% of the mites in 48 h, which was similar with those of the broth. Bioassays of the supernatant proteins confirmed that various proteins exhibited acaricidal activity. Five candidate proteins were expressed and purified successfully. The bioassays showed that the small protein BVP8 exhibited significant acaricidal activity with an LC₅₀ of 12.4 μg/mL (T. urticae). Overall, these findings suggest that B. vallismortis NBIF-001 is a potential biocontrol agent for spider mite management.

The H92R substitution in PSST is a reliable diagnostic biomarker for predicting resistance to mitochondrial electron transport inhibitors of complex I in European populations of Tetranychus urticae

BACKGROUND

Mitochondrial Electron Transport Inhibitors of complex I (METI-I), such as tebufenpyrad and fenpyroximate, are acaricides that have been used extensively to control Tetranychus urticae Koch (Acari: Tetranychidae) for more than 20 years. Because of the ability of this spider mite to rapidly develop acaricide resistance, field (cross-) resistance monitoring and elucidation of resistance mechanisms are extremely important for resistance management (RM). In the present study, 42 European T. urticae field populations were screened for tebufenpyrad and fenpyroximate resistance, and the correlation between resistance and the H92R substitution in PSST was investigated.

RESULTS

According to the calculated lethal concentration values that kill 90% of the population (LC₉₀ ), tebufenpyrad and fenpyroximate would fail to control many of the collected populations at recommended field rates. Six populations exhibited high to very high resistance levels (200- to over 1950-fold) to both METI-Is. Analysis based on the LC₅₀ values displayed a clear correlation between tebufenpyrad and fenpyroximate resistance, further supporting cross-resistance, which is of great operational importance in acaricide RM. The previously uncovered METI-I target-site mutation H92R in the PSST homologue of complex I (NADH:ubiquinone oxidoreductase) was found with high allele frequencies in populations resistant to tebufenpyrad and fenpyroximate. Synergist assays showed this mutation is not the only factor involved in METI-I resistance and additive or synergistic effects of multiple mechanisms most likely determine the phenotypic strength.

CONCLUSIONS

The predictive value of resistance by H92R is very high in European populations and offers great potential to be used as a molecular diagnostic marker for METI-I resistance. © 2022 Society of Chemical Industry.

High frequency of ryanodine receptor and cytochrome P450 CYP9A186 mutations in insecticide-resistant field populations of Spodoptera exigua from China

The beet armyworm, Spodoptera exigua is a global agricultural pest that is polyphagous, highly dispersive, and often difficult to control due to resistance to many insecticides. Previous studies showed that a target site mutation in the S. exigua ryanodine receptor (SeRyR) corresponding to I4743M contributes approximately 20-fold resistance to chlorantraniliprole, whereas a mutation in the cytochrome P450 enzyme CYP9A186 corresponding to F116V confers 200-fold to emamectin benzoate through enhanced metabolic detoxification. Here, high frequencies of mutations were found among six China S. exigua field populations collected from 2016 to 2019 resulting in SeRyR I4743M and CYP9A186 F116V substitutions, with some populations having high levels of resistance to chlorantraniliprole and emamectin benzoate, respectively. Whereas we found a significant correlation between emamectin benzoate resistance level and the allele frequency of CYP9A186 F116V, no significant correlation was found between chlorantraniliprole resistance level and SeRyR I4743M allele frequency in the six field populations. These results suggest that CYP9A186 F116V is a major resistance mechanism for emamectin benzoate in the tested field populations, whereas it is likely that resistance mechanisms other than SeRyR I4743M are responsible for resistance to chlorantraniliprole in the six China field populations. Because of the growing resistance to these two insecticides by S. exigua in China, the use of insecticidal compounds with different modes of action and/or other integrated pest management strategies are needed to further delay the evolution of insecticide resistance and effectively manage S. exigua in China.

Overexpression of Multiple Cytochrome P450 Genes Conferring Clothianidin Resistance in Bradysia odoriphaga

Cytochrome P450 monooxygenases (P450s) play important roles in the detoxification metabolism of xenobiotics and are involved in the resistance of insects to many insecticides. In this study, piperonyl butoxide (PBO), an inhibitor of P450 enzyme activity, significantly increased the toxicity of clothianidin in the clothianidin-resistant (CL-R) population of Bradysia odoriphaga. The enzyme activity of P450 in the CL-R population was significantly higher than that in the SS population. Furthermore, four P450 genes were found to be significantly overexpressed in the CL-R population. Tissue-specific expression analysis indicates that CYP9J57, CYP3828A1, CYP6SX1, and CYP6QE1 were most highly expressed in the midgut and/or Malpighian tubules. After exposure to LC₃₀ of clothianidin, the expression levels of the four P450 genes were significantly upregulated. The RNAi-mediated knockdown of CYP9J57, CYP3828A1, and CYP6QE1 significantly increased the susceptibility of B. odoriphaga to clothianidin. These results suggest that P450 genes are involved in clothianidin resistance in B. odoriphaga. This provides a better understanding of P450-mediated clothianidin resistance in B. odoriphaga and will contribute to the management of insect resistance to insecticides.

RNAi Mediated Gene Silencing of Detoxification Related Genes in the Ectropis oblique

Ectropis oblique is one of the main pests that feed on tea leaves. At present, the main control method is chemical control, but the long-term use of insecticides has been related to the development of insect resistance. One of the resistance mechanisms is the upregulation of relevant detoxification enzymes for defense. In this study, four genes with increased expression were screened from the gene sequences annotated from the transcriptome data of deltamethrin-treated larvae of E. oblique, which are acid phosphatase EoACP138, and cytochrome P450 EoCYP316, carboxylesterase EoCarE592 and acetylcholine esterase EoAchE989, respectively. The fourth instar larvae of E. oblique were stimulated by deltamethrin, chlorpyrifos and fenpropathrin respectively, and the expression levels of the genes were detected by qRT-PCR. The result showed that all four genes’ expression had significantly increased under the stimulation of three insecticides. RNAi technology was used to silence the expression of genes of EoACP138, EoCYP316, EoCarE592 and EoAchE989 in the fourth instar larvae of E. oblique. The change in the expression levels of the above genes in the larvae treated with dsRNA and stimulated with pesticides was determined by qRT-PCR. The target genes have been effectively silenced after feeding on dsRNA and higher sensitivity with higher mortality to pesticides was observed in the larvae interfered with dsRNA. The above genes are related to the detoxification and metabolism of resistance of E. oblique, which lays a foundation for further study on the mechanism of insecticide resistance in E. oblique.

Pesticide resistance in arthropods: Ecology matters too

Pesticide resistance development is an example of rapid contemporary evolution that poses immense challenges for agriculture. It typically evolves due to the strong directional selection that pesticide treatments exert on herbivorous arthropods. However, recent research suggests that some species are more prone to evolve pesticide resistance than others due to their evolutionary history and standing genetic variation. Generalist species might develop pesticide resistance especially rapidly due to pre‐adaptation to handle a wide array of plant allelochemicals. Moreover, research has shown that adaptation to novel host plants could lead to increased pesticide resistance. Exploring such cross‐resistance between host plant range evolution and pesticide resistance development from an ecological perspective is needed to understand its causes and consequences better. Much research has, however, been devoted to the molecular mechanisms underlying pesticide resistance while both the ecological contexts that could facilitate resistance evolution and the ecological consequences of cross‐resistance have been under‐studied. Here, we take an eco‐evolutionary approach and discuss circumstances that may facilitate cross‐resistance in arthropods and the consequences cross‐resistance may have for plant–arthropod interactions in both target and non‐target species and species interactions. Furthermore, we suggest future research avenues and practical implications of an increased ecological understanding of pesticide resistance evolution.

Intradiol ring cleavage dioxygenases from herbivorous spider mites as a new detoxification enzyme family in animals

Background

Generalist herbivores such as the two-spotted spider mite Tetranychus urticae thrive on a wide variety of plants and can rapidly adapt to novel hosts. What traits enable polyphagous herbivores to cope with the diversity of secondary metabolites in their variable plant diet is unclear. Genome sequencing of T. urticae revealed the presence of 17 genes that code for secreted proteins with strong homology to “intradiol ring cleavage dioxygenases (DOGs)” from bacteria and fungi, and phylogenetic analyses show that they have been acquired by horizontal gene transfer from fungi. In bacteria and fungi, DOGs have been well characterized and cleave aromatic rings in catecholic compounds between adjacent hydroxyl groups. Such compounds are found in high amounts in solanaceous plants like tomato, where they protect against herbivory. To better understand the role of this gene family in spider mites, we used a multi-disciplinary approach to functionally characterize the various T. urticae DOG genes.

Results

We confirmed that DOG genes were present in the T. urticae genome and performed a phylogenetic reconstruction using transcriptomic and genomic data to advance our understanding of the evolutionary history of spider mite DOG genes. We found that DOG expression differed between mites from different plant hosts and was induced in response to jasmonic acid defense signaling. In consonance with a presumed role in detoxification, expression was localized in the mite’s gut region. Silencing selected DOGs expression by dsRNA injection reduced the mites’ survival rate on tomato, further supporting a role in mitigating the plant defense response. Recombinant purified DOGs displayed a broad substrate promiscuity, cleaving a surprisingly wide array of aromatic plant metabolites, greatly exceeding the metabolic capacity of previously characterized microbial DOGs.

Conclusion

Our findings suggest that the laterally acquired spider mite DOGs function as detoxification enzymes in the gut, disarming plant metabolites before they reach toxic levels. We provide experimental evidence to support the hypothesis that this proliferated gene family in T. urticae is causally linked to its ability to feed on an extremely wide range of host plants.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01323-1.

Resistance profile of Rhipicephalus (Boophilus) microplus to diazinon and cypermethrin and first report of sodium channel mutation - domain III S6 - T2134A, in field samples from the state of São Paulo, Brazil

ABSTRACT Rhipicephalus (Boophilus) microplus is one of the parasites that has an impact on livestock farming in Brazil. Bioassays using the larval packet test (LPT) impregnated with cypermethrin and diazinon at different concentrations were performed to characterize phenotypic resistance (resistance level, RL), molecular resistance (mutation in domain III - S6 T2134A) and enzymatic metabolism of diazinon and cypermethrin in some populations of R. microplus that were collected in different regions of the state of São Paulo. Among the 40 farms analyzed, 18 of them presented resistance factors for cypermethrin, with RL I (53%) and RL II (47%). Mutation of domain III T2134A was found for the first time in Brazil, on seven farms. Heterozygous larvae were found on six of them and resistant homozygous larvae on four. No differences (p>0.05) in enzymatic activity of α-esterase and β-esterase were found in tests with live and dead larvae at a cypermethrin concentration of 409.6μg/cm2. Ninety percent of the farms showed resistance to the active agent diazinon: RL I in 6%, RL II in 30% and RL III in 64%. There were significant differences (p<0.05) in enzymatic activity at the highest concentration (3.2μg/cm2).

Genomic features of the polyphagous cotton leafworm Spodoptera littoralis

BACKGROUND

The cotton leafworm, Spodoptera littoralis, is a highly polyphagous pest of many cultivated plants and crops in Africa and Europe. The genome of this pest will help us to further understand the molecular mechanisms of polyphagy.

RESULTS

Herein, the high-quality genome of S. littoralis was obtained by Pacific Bioscience (PacBio) sequencing. The assembled genome size of S. littoralis is 436.55 Mb with a scaffold N50 of 6.09 Mb, consisting of 17,207 annotated protein-coding genes. Phylogenetic analysis shows that S. littoralis and its sibling species S. litura diverged about 5.44 million years ago. Expanded gene families were mainly involved in metabolic detoxification and tolerance to toxic xenobiotics based on GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis. Comparative genomics analysis showed that gene families involved in detoxification and chemosensation were significantly expanded in S. littoralis, representing genetic characteristics related to polyphagy and an extensive host range.

CONCLUSIONS

We assembled and annotated the reference genome of S. littoralis, and revealed that this pest has the genetic features of strong detoxification capacity, consistent with it being a significant risk to a wide range of host crops. These data resources will provide support for risk assessment and early warning monitoring of major polyphagous agricultural pests.

A H258Y mutation in subunit B of the succinate dehydrogenase complex of the spider mite Tetranychus urticae confers resistance to cyenopyrafen and pyflubumide, but likely reinforces cyflumetofen binding and toxicity

Succinate dehydrogenase (SDH) inhibitors such as cyflumetofen, cyenopyrafen and pyflubumide, are selective acaricides that control plant-feeding spider mite pests. Resistance development to SDH inhibitors has been investigated in a limited number of populations of the spider mite Tetranychus urticae and is associated with cytochrome P450 based detoxification and target-site mutations such as I260 T/V in subunit B and S56L in subunit C of SDH. Here, we report the discovery of a H258Y substitution in subunit B of SDH in a highly pyflubumide resistant population of T. urticae. As this highly conserved residue corresponds to one of the ubiquinone binding residues in fungi and bacteria, we hypothesized that H258Y could have a strong impact on SDH inhibitors toxicity. Marker assisted introgression and toxicity bioassays revealed that H258Y caused high cross resistance between cyenopyrafen and pyflubumide, but increased cyflumetofen toxicity. Resistance associated with H258Y was determined as dominant for cyenopyrafen, but recessive for pyflubumide. In vitro SDH assays with extracted H258 mitochondria showed that cyenopyrafen and the active metabolites of pyflubumide and cyflumetofen, interacted strongly with complex II. However, a clear shift in IC50s was observed for cyenopyrafen and the metabolite of pyflubumide when Y258 mitochondria were investigated. In contrast, the mutation slightly increased affinity of the cyflumetofen metabolite, likely explaining its increased toxicity for the mite lines carrying the substitution. Homology modeling and ligand docking further revealed that, although the three acaricides share a common binding motif in the Q-site of SDH, H258Y eliminated an important hydrogen bond required for cyenopyrafen and pyflubumide binding. In addition, the hydrogen bond between cyenopyrafen and Y117 in subunit D was also lost upon mutation. In contrast, cyflumetofen affinity was enhanced due to an additional hydrogen bond to W215 and hydrophobic interactions with the introduced Y258 in subunit B. Altogether, our findings not only highlight the importance of the highly conserved histidine residue in the binding of SDH inhibitors, but also reveal that a resistance mutation can provide both positive and negative cross-resistance within the same acaricide mode of action group.

Comparative transcriptome and RNA interference reveal CYP6DC1 and CYP380C47 related to lambda-cyhalothrin resistance in Rhopalosiphum padi

The bird-cherry-oat aphid, Rhopalosiphum padi, is a serious agricultural pest of Triticeae crops, and pyrethroids are the most widely used chemical pesticides for the control of the aphid. Our previous studies found that some R. padi field populations have developed resistance against pyrethroids; an M918L target-site mutation of the voltage gated sodium channel was present in the pyrethroid resistant individuals, while the high-level resistance to lambda-cyhalothrin revealed the presence of other mechanisms in the pest. Here, we conducted genome-wide transcriptional analysis for the lambda-cyhalothrin susceptible (SS) and resistant (LC-RR) strains of R. padi. Results indicated that 2457 genes were differently expressed between the SS and LC-RR strains. In the LC-RR, a total of 1265 and 1192 genes were up- and down-regulated, respectively. KEGG analysis implicated enrichment of P450 involved in insecticide metabolic pathways in the resistant transcriptome. qRT-PCR results confirmed that two P450 genes (CYP6DC1 and CYP380C47) were significantly overexpressed in the LC-RR individuals. Furthermore, RNA interference (RNAi) of CYP6DC1 or CYP380C47 significantly increased mortality of R. padi exposure to lambda-cyhalothrin. These results suggest that the overexpression of CYP6DC1 and CYP380C47 contributed to the lambda-cyhalothrin resistance in the pest. This study provides knowledge for further analyzing the molecular mechanism of resistance to pyrethroids in R. padi.

Composition of mosquito fauna and insecticide resistance status of Anopheles gambiae sensu lato in Itang special district, Gambella, Southwestern Ethiopia

Background

Anopheles arabiensis, member species of the Anopheles gambiae complex, is the primary vector of malaria and is widely distributed in Ethiopia. Anopheles funestus, Anopheles pharoensis and Anopheles nili are secondary vectors occurring with limited distribution in the country. Indoor residual spraying (IRS) and long-lasting insecticidal nets (LLINs) are pillars for the interventions against malaria control and elimination efforts in Ethiopia. However, the emergence and widespread of insecticide resistance in An. gambiae sensu lato (s.l.), might compromise the control efforts of the country. The aim of this study was to investigate composition of mosquito fauna and insecticide resistance status of An. gambiae s.l. in Itang special district ( woreda), Gambella, southwestern Ethiopia.

Methods

Adult mosquitoes were sampled from September 2020 to February 2021 using the CDC light trap and pyrethrum spray catch (PSC). CDC light traps were placed in three selected houses for two consecutive days per month to collect mosquitoes indoor and outdoor from 6:00 P.M. to 06:00 A.M. and PSC was used to collect indoor resting mosquitoes from ten selected houses once in a month from October 2020 to February 2021. Moreover, mosquito larvae were also collected from different breeding sites and reared to adults to assess susceptibility status of populations of An. gambiae s.l. in the study area. Susceptibility tests were conducted on two to three days old non blood fed female An. gambiae s.l. using insecticide impregnated papers with deltamethrin (0.05%), alpha-cypermethrin (0.05%), propoxur (0.1%), pirimiphos-methyl (0.25%) and bendiocarb (0.1%) following World Health Organization (WHO) standard susceptibility test procedure. Molecular diagnostics were done for the identification of member species of An. gambiae s.l. and detection of knockdown resistance (kdr) allele using species specific polymerase chain reaction (PCR) and allele specific PCR.

Results

In total, 468 adult mosquitoes were collected from different houses. Culex mosquitoes were the most dominant (80.4%) followed by Anopheles mosquitoes. Three species of Anopheles (Anopheles coustani, An. pharoensis, and An. gambiae s.l.) were identified, of which An. coustani was the dominant (8.1%) species. Higher number of mosquitoes (231) were collected outdoor by CDC light traps. Out of 468 adult mosquitoes, 294 were blood fed, 46 were half-gravid and gravid whereas the remaining 128 were unfed. WHO bioassay tests revealed that the populations of An. gambiae s.l. in the study area are resistant against alpha-cypermethrin and deltamethrin, but susceptible to bendiocarb, pirimiphos-methyl and propoxur. Of the total 86 An. gambiae s.l. specimens assayed, 79 (92%) successfully amplified and identified as An. arabiensis. West African kdr (L1014F) mutation was detected with high kdr allele frequency ranging from 67 to 88%.

Conclusion

The detection of target site mutation, kdr L1014F allele, coupled with the phenotypic resistance against alpha-cypermethrin and deltamethrin call for continuous resistance monitoring.

Transcriptomic profile of the predatory mite Amblyseius swirskii (Acari: Phytoseiidae) on different host plants

Amblyseius swirskii Athias-Henriot (Acari: Phytoseiidae) is a predatory mite, effective at controlling whiteflies and thrips in protected crops. However, on tomato its efficacy as a biocontrol agent is hindered, most probably by the plant trichomes and their exudates. Our aim was to characterize the response of A. swirskii to the tomato trichome exudates and identify three major detoxification gene sets in this species: cytochromes P450 (CYPs), glutathione S-transferases (GSTs) and carboxyl/cholinesterases (CCEs). Mites were exposed separately to tomato and pepper, a favourable host plant for A. swirskii, after which their transcriptional responses were analysed and compared. The de novo transcriptome assembly resulted in 71,336 unigenes with 66.1% of them annotated. Thirty-nine A. swirskii genes were differentially expressed after transfer on tomato leaves when compared to pepper leaves; some of the expressed genes were associated with the metabolism of tomato exudates. Our results illustrate that the detoxification gene sets CYPs, GSTs and CCEs are abundant in A. swirskii, but do not play a significant role when in contact with the tomato exudates.

Historical samples reveal a combined role of agriculture and public-health applications in vector resistance to insecticides

BACKGROUND

Insecticide resistance is the major threat to vector control and for the prevention of vector-borne diseases. Because almost all insecticides used against insect vectors are or have been used in agriculture, a connection between agricultural insecticide use and resistance in insect vectors has been hypothesized. However, it is challenging to find a causal link between past agricultural use of insecticides and current resistance in vector populations without historical data series. Here we investigated the relative contribution across time of agricultural and public-health insecticide applications in selecting for diflubenzuron (DFB) resistance in Culex pipiens populations. Using DNA sequencing, we looked for DFB resistant mutations in current and historical mosquito samples, dating back to the 1980s-1990s, when DFB was used in agriculture but not yet in mosquito control.

RESULTS

In the samples collected before the introduction of DFB in vector control, we found the resistant mutation I1043M in rural regions but not any of the neighboring urban and natural areas, indicating that the selection pressure was derived by agriculture. However, after the introduction of DFB for vector control, the resistant mutations were found across all study areas showing that the initial selection from agriculture was further boosted by the selection pressure imposed by the mosquito control applications in the 2000s.

CONCLUSIONS

Our findings support a combined role of agricultural and public-health use of insecticides in vector resistance across time and call for specific actions in integrated resistance management, including increased communication between agriculture and health practitioners. © 2022 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

ddRAD Sequencing Identifies Pesticide Resistance-Related Loci and Reveals New Insights into Genetic Structure of Bactericera cockerelli as a Plant Pathogen Vector

(1) Background: Many hemipteran insects transmit plant pathogens that cause devastating crop diseases, while pest management frequently relies primarily on insecticide applications. These intense insecticide applications lead to the development of insecticide resistance, as was the case for potato psyllid, Bactericera cockerelli (Hemiptera: Triozidae), a vector of Candidatus Liberibacter solanacearum, which causes zebra chip disease in potato. (2) Methods: Here, we use double-digest restriction site-associated DNA (ddRAD) to genotype eight psyllid populations (one susceptible and seven resistant to neonicotinoid insecticides). (3) Results: Association tests identified over 400 loci that were strongly segregated between susceptible and resistant populations. Several loci were located within genes involved in insecticide resistance, gene regulation, fertility, and development. Moreover, we explored the genetic structure of these eight populations and discovered that routinely utilized haplotyping was not an accurate predictor of population structure. Pairwise comparisons of the fixation index (FST) of populations of the same haplotype were not different from pairwise FST of populations that belonged to different haplotypes. (4) Conclusions: Our findings suggest that neonicotinoid insecticide resistance has a genetic basis, most likely as a result of similar selection pressure. Furthermore, our results imply that using a single maternally inherited gene marker to designate genetic lineages for potato psyllids should be re-evaluated.

Selectivity and molecular stress responses to classical and botanical acaricides in the predatory mite Phytoseiulus persimilis Athias-Henriot (Acari: Phytoseiidae)

BACKGROUND

Acaricide application remains an integral component of integrated pest management (IPM) for the two-spotted spider mite Tetranychus urticae. Species and strains of phytoseiid predatory mites vary significantly in their response to acaricides. For the success of IPM, it is imperative to identify the determinants of selectivity and molecular stress responses of acaricides in predatory mites.

RESULTS

The three classical acaricides bifenazate, cyflumetofen, and fenbutatin oxide did not affect the survival and fecundity of Phytoseiulus persimilis regardless of the route of exposure. Selectivity of the orange oil and terpenoid blend-based botanical acaricides was low via a combination of direct exposure, acaricide-laced diet, and residual exposure but improved when limiting exposure only to diet. To gain insights into the molecular stress responses, the transcriptome of P. persimilis was assembled. Subsequent gene expression analysis of predatory mites orally exposed to fenbutatin oxide and orange oil yielded only a limited xenobiotic stress response. In contrast, P. persimilis exhibited target-site resistance mutations, including I260M in SdhB, I1017M in CHS1, and kdr and super-kdr in VGSC. Extending the screen using available Phytoseiidae sequences uncovered I136T, S141F in cytb, G119S in AChE, and A2083V in ACC, well-known target-sites of acaricides.

CONCLUSION

Selectivity of the tested botanical acaricides to P. persimilis was low but could be enhanced by restricting exposure to a single route. Differential gene expression analysis did not show a robust induced stress response after sublethal exposure. In contrast, this study uncovered target-site mutations that may help to explain the physiological selectivity of several classical acaricides to phytoseiid predators.

Over-expression in cis of the midgut P450 CYP392A16 contributes to abamectin resistance in Tetranychus urticae

Cytochrome P450 mediated metabolism is a well-known mechanism of insecticide resistance. However, to what extent qualitative or quantitative changes are responsible for increased metabolism, is not well understood. Increased expression of P450 genes is most often reported, but the underlying regulatory mechanisms remain widely unclear. In this study, we investigate CYP392A16, a P450 from the polyphagous and major agricultural pest Tetranychus urticae. High expression levels of CYP392A16 and in vitro metabolism assays have previously associated this P450 with abamectin resistance. Here, we show that CYP392A16 is primarily localized in the midgut epithelial cells, as indicated by immunofluorescence analysis, a finding also supported by a comparison between feeding and contact toxicity bioassays. Silencing via RNAi of CYP392A16 in a highly resistant T. urticae population reduced insecticide resistance levels from 3400- to 1900- fold, compared to the susceptible reference strain. Marker-assisted backcrossing, using a single nucleotide polymorphism (SNP) found in the CYP392A16 allele from the resistant population, was subsequently performed to create congenic lines bearing this gene in a susceptible genetic background. Toxicity assays indicated that the allele derived from the resistant strain confers 3.6-fold abamectin resistance compared to the lines with susceptible genetic background. CYP392A16 is over-expressed at the same levels in these lines, pointing to cis-regulation of gene expression. In support of that, functional analysis of the putative promoter region from the resistant and susceptible parental strains revealed a higher reporter gene expression, confirming the presence of cis-acting regulatory mechanisms.

Chemosensory Proteins Are Associated with Thiamethoxam and Spirotetramat Tolerance in Aphis gossypii Glover

Chemosensory proteins (CSPs) are a class of transporters in arthropods. Deeper research on CSPs showed that CSPs may be involved in some physiological processes beyond chemoreception, such as insect resistance to pesticides. We identified two upregulated CSPs in two resistant strains of Aphis gossypii Glover. To understand their role in the resistance of aphids to pesticides, we performed the functional verification of CSP1 and CSP4 in vivo and in vitro. Results showed that the sensitivity of the thiamethoxam-resistant strain to thiamethoxam increased significantly with the silencing of CSP1 and CSP4 by RNAi (RNA interference), and the sensitivity of the spirotetramat-resistant strain to spirotetramat increased significantly with the silencing of CSP4. Transgenic Drosophila melanogaster expressing CSPs exhibited stronger resistance to thiamethoxam, spirotetramat, and alpha-cypermethrin than the control did. In the bioassay of transgenic Drosophila, CSPs showed different tolerance mechanisms for different pesticides, and the overexpressed CSPs may play a role in processes other than resistance to pesticides. In brief, the present results prove that CSPs are related to the resistance of cotton aphids to insecticides.

The P450 genes of the cat flea, Ctenocephalides felis: a CYPome in flux

The genome of the cat flea, an ectoparasite of major veterinary importance and the first representative of the Siphonaptera, is highly unusual among arthropod genomes in showing a variable size and a very large number of gene duplications (Driscoll et al., 2020). The cat flea is the target of several classes of insecticides, justifying the description of its CYPome, the complement of P450s that are an important family of detoxification enzymes. 103 P450 genes were annotated on the nine chromosomes, with an additional 12 genes on small, extrachromosomal scaffolds. Only 34 genes were found as single sequences, with 47 duplicated two to four-fold. This included duplication of genes that are mostly single copy P450 genes in other arthropods. Large clusters of mitochondrial clan P450s were observed, resulting in a CYP12 bloom within this clan to 34 genes, a number of mitochondrial P450s not seen in other animals so far. The variable geometry of the cat flea CYPome poses a challenge to the study of P450 function in this species, and raises the question of the underlying causes of single copy control versus multicopy licence of P450 genes.

Combination of target site mutation and associated CYPs confers high-level resistance to pyridaben in Tetranychus urticae

Pyridaben is a mitochondrial electron transport complex I inhibitor. The H110R mutation in the PSST subunit has been reported as a major factor in pyridaben resistance in the two-spotted spider mite, Tetranychus urticae. However, backcross experiments revealed that the mutant PSST alone conferred only moderate resistance. In contrast, inhibition of cytochrome P450 (CYP) markedly reduces resistance levels in a number of highly resistant strains. It was reported previously that maternal factors contributed to the inheritance of pyridaben resistance in the egg stage, but the underlying mechanisms have yet to be elucidated. Here, we studied the combined effects of the PSST H110R mutation and candidate CYPs, as metabolic resistance factors, on pyridaben resistance in T. urticae. We found that the maternal effects of inheritance of resistance in the egg stage were associated with CYP activity. Analysis of differential gene expression by RNA-seq identified CYP392A3 as a candidate causal factor for the high resistance level. Congenic strains, where the alleles of both PSST and CYP392A3 were derived from a resistant strain (RR_i; i = 1 or 2) and a susceptible strain (SS_i) in a common susceptible genetic background, were constructed by marker-assisted backcrossing. RR_i showed upregulation of CYP392A3 and high resistance levels (LC₅₀ > 10,000 mg L^-1), while SS_i had LC₅₀ < 10 mg L^-1. To disentangle the individual effects of PSST and CYP392A3 alleles, we also attempted to uncouple these genes in RR_i. We conclude that given the variation in LC₅₀ values and expression levels of CYP392A3 in the congenic and uncoupled strains, it is likely that the high pyridaben resistance levels are due to a synergistic or cumulative effect of the combination of mutant PSST and associated CYPs, including CYP392A3, but other yet to be discovered factors cannot be excluded.

Genome resequencing reveals rapid, repeated evolution in the Colorado potato beetle

Insecticide resistance and rapid pest evolution threatens food security and the development of sustainable agricultural practices, yet the evolutionary mechanisms that allow pests to rapidly adapt to control tactics remains unclear. Here we examine how a global super-pest, the Colorado potato beetle (CPB), Leptinotarsa decemlineata, rapidly evolves resistance to insecticides. Using whole genome resequencing and transcriptomic data focused on its ancestral and pest range in North America, we assess evidence for three, non-mutually exclusive models of rapid evolution: pervasive selection on novel mutations, rapid regulatory evolution, and repeated selection on standing genetic variation. Population genomic analysis demonstrates that CPB is geographically structured, even among recently established pest populations. Pest populations exhibit similar levels of nucleotide diversity, relative to non-pest populations, and show evidence of recent expansion. Genome scans provide clear signatures of repeated adaptation across CPB populations, with especially strong evidence of selection on insecticide resistance genes in different populations. Analyses of gene expression show that constitutive upregulation of candidate insecticide resistance genes drives distinctive population patterns. CPB evolves insecticide resistance repeatedly across agricultural regions, leveraging similar genetic pathways but different genes, demonstrating a polygenic trait architecture for insecticide resistance that can evolve from standing genetic variation. Despite expectations, we do not find support for strong selection on novel mutations, or rapid evolution from selection on regulatory genes. These results suggest that integrated pest management practices must mitigate the evolution of polygenic resistance phenotypes among local pest populations, in order to maintain the efficacy and sustainability of novel control techniques.

Biochemical and molecular mechanisms of acaricide resistance in Dermanyssus gallinae populations from Turkey

The poultry red mite, Dermanyssus gallinae, is the most important blood sucking ectoparasite of egg laying hens and causes economic losses in poultry farms worldwide. Although various management methods exist, the control of poultry red mites (PRMs) mainly relies on acaricides such as pyrethroids and organophosphates (OPs) in many regions of the world. However, repeated use of these synthetic chemicals has resulted in resistance development causing control failure of PRM. In this study, we investigated acaricide resistance mechanisms of Turkish PRM populations. First, we obtained the COI sequence from 30 PRM populations from different regions in Turkey and identified four different COI haplotypes. Toxicity assays showed that four field-collected PRM populations were highly resistant to the pyrethroid alpha-cypermethrin, with resistance ratios (RRs) varying between 100- and 400-fold, while two of these populations had a RR of more than 24-fold against the OP acaricide phoxim. Biochemical assays showed a relatively higher activity of glutathione-S-transferases and carboxyl-cholinesterases, two well-known classes of detoxification enzymes, in one of these resistant populations. In addition, we also screened for mutations in the gene encoding the voltage-gated sodium channel (vgsc) and acetylcholinesterase 1 (ace-1), the target-site of pyrethroids and OPs, respectively. In all but two PRM populations, at least one vgsc mutation was detected. A total of four target-site mutations, previously associated with pyrethroid resistance, M918T, T929I, F1534L, F1538L were found in domain II and III of the VGSC. The T929I mutation was present in the vgsc of almost all PRM populations, while the other mutations were only found at low frequency. The G119S/A mutation in ace-1, previously associated with OP resistance, was found in PRM for the first time and present in fourteen populations. Last, both alive and dead PRMs were genotyped after pesticide exposure and supported the possible role of target-site mutations, T929I and G119S, in alpha-cypermethrin and phoxim resistance, respectively. To conclude, our study provides a current overview of resistance levels and resistance mutations in Turkish PRM populations and might aid in the design of an effective resistance management program of PRM in Turkey.

Pirimicarb resistance and associated mechanisms in field-collected and selected populations of Neoseiulus californicus

The predatory mite Neoseiulus californicus McGregor (Acari: Phytoseiidae) is an important natural enemy of phytophagous mites, and naturally established populations are often found in apple orchards. However, insecticide applications to control pests cause side effects to non-target organisms such as N. californicus. Pirimicarb, a widely used carbamate insecticide in apple orchards, is generally considered a selective aphidicide, however, toxicity to beneficial insects and predatory mites has been reported. Furthermore, the molecular basis for this selectivity, if present in N. californicus, is still largely unknown. In this study, 8 field-collected N. californicus populations were investigated and showed up to 27-fold resistance compared to a susceptible laboratory population. Selection in the laboratory for 5 consecutive generations resulted in a 69-fold pirimicarb resistance. Although there were no significant difference in terms of the acetlycholinesterase (AChE) activity between susceptible and field-collected populations, the selected population exhibited a significantly higher AChE activity. In addition, gene copy number variation of acetylcholinesterase (ace) gene among populations was detected and ranged from 1.6 to 2.1-fold relative to the susceptible population. All field-collected populations, but not the selected population, had a significantly higher ace copy number compared to the susceptible population (t-test, p < 0.05). Molecular analysis of the target-site (AChE) revealed, for the first time, a phenylalanine to tryptophan substition at position 331 in AChE (Torpedo californica numbering), both in field-collected and the selected population, but not in the susceptible population. Last, the selected F5 population consumed significantly more Tetranychusurticae adults than the parental population. Together, the results of this study shed light on the molecular determinants of acaricide selectivity in predatory mites, and will contribute to a better design of an integrated mite management program, including the use of pesticide resistant N. californicus in apple orchards.

Genome evolution in an agricultural pest following adoption of transgenic crops

Evolution of resistance to management approaches in agricultural landscapes is common and results in economic losses. Early detection of pest resistance prior to significant crop damage would benefit the agricultural community. It has been hypothesized that new genomic approaches could track molecular signals of emerging resistance and trigger efforts to preempt widespread damage. We tested this hypothesis by quantifying genomic changes in the pest Helicoverpa zea over a 15-y period concurrent with commercialization of transgenic Bacillus thuringiensis–expressing crops and their subsequent loss of efficacy. Our results demonstrate the complex nature of evolution in agricultural ecosystems and provide insight into the potential and pitfalls of using genomic approaches for resistance monitoring.

Replacing synthetic insecticides with transgenic crops for pest management has been economically and environmentally beneficial, but these benefits erode as pests evolve resistance. It has been proposed that novel genomic approaches could track molecular signals of emerging resistance to aid in resistance management. To test this, we quantified patterns of genomic change in Helicoverpa zea, a major lepidopteran pest and target of transgenic Bacillus thuringiensis (Bt) crops, between 2002 and 2017 as both Bt crop adoption and resistance increased in North America. Genomic scans of wild H. zea were paired with quantitative trait locus (QTL) analyses and showed the genomic architecture of field-evolved Cry1Ab resistance was polygenic, likely arising from standing genetic variation. Resistance to pyramided Cry1A.105 and Cry2Ab2 toxins was controlled by fewer loci. Of the 11 previously described Bt resistance genes, 9 showed no significant change over time or major effects on resistance. We were unable to rule out a contribution of aminopeptidases (apns), as a cluster of apn genes were found within a Cry-associated QTL. Molecular signals of emerging Bt resistance were detectable as early as 2012 in our samples, and we discuss the potential and pitfalls of whole-genome analysis for resistance monitoring based on our findings. This first study of Bt resistance evolution using whole-genome analysis of field-collected specimens demonstrates the need for a more holistic approach to examining rapid adaptation to novel selection pressures in agricultural ecosystems.

Variation of diazinon and amitraz susceptibility of Hyalomma marginatum (Acari: Ixodidae) in the Rabat-Sale-Kenitra region of Morocco

In the present study, the acaricide susceptibility status of Hyalomma marginatum in Morocco was investigated in the Rabat-Sale-Kenitra region using the Larval Packet Test. The overall LC₅₀ value for diazinon and amitraz was 115 ppm (95% CI: [104; 125]) and 22 ppm (95% CI: [21; 23]), respectively. The LC₅₀ values varied significantly between the nine sampled locations (P<0.001) ranging from 75 ppm (95% CI: [65; 84]) in Ouelmes to 179 ppm (95% CI: [139; 201]) in Jorf El Melha for diazinon and from 18 ppm (95% CI: [15; 21]) in Skhirat to 28 ppm (95% CI: [24; 31]) in Ouelmes for amitraz. Sequencing of the target-site of diazinon, acetylcholinesterase 1 (AChE1), indicated that previously reported resistance mutations in AChE1 were absent in ticks from Jorf El Melha surviving 500 ppm diazinon. This study is the first report on the H. marginatum susceptibilty status to the most frequently used acaricides in Morocco and indicates that acaricide tick resistance is emerging.

The G126S substitution in mitochondrially encoded cytochrome b does not confer bifenazate resistance in the spider mite Tetranychus urticae

Several genetic variants of the cd1- and ef-helices of the Q_o site of mitochondrial cytochrome b have been associated with bifenazate resistance in the spider mite Tetranychus urticae, an important crop pest around the world. Maternal inheritance of bifenazate resistance has provided strong evidence for the involvement of many of these mutations alone or in combination. A number of populations highly resistant to bifenazate were uncovered that carried the G126S substitution in combination with other target-site mutations. This G126S mutation has therefore been investigated in several studies in the context of resistance evolution and the development of diagnostic markers. However, experimental data that link bifenazate resistance with the presence of the G126S mutation without additional cd1- and ef-helices mutations, remain very limited. Here, we genotyped 38 T. urticae field populations for cytochrome b and uncovered nine field populations with a fixed or segregating G126S substitution without other target-site mutations in the conserved cd1- and ef-helices of the cytochrome b Q_o pocket. Toxicity bioassays showed that all nine field populations were very susceptible to bifenazate, providing strong evidence that G126S alone does not confer bifenazate resistance. These findings also implicate that previous T. urticae populations with G126S found to be low to moderately resistant to bifenazate, evolved alternative mechanisms of resistance, and more importantly, that this mutation cannot be used as a molecular diagnostic for bifenazate resistance.

Pyrethroid target-site resistance mutations in populations of the honey bee parasite Varroa destructor (Acari: Varroidae) from Flanders, Belgium

The honey bee ectoparasite Varroa destructor is considered the major threat to apiculture, as untreated colonies of Apis mellifera usually collapse within a few years. In order to control this mite, many beekeepers rely on a limited number of approved synthetic acaricides, including the pyrethroids tau-fluvalinate and flumethrin. Due to the intensive use of these products, resistance is now commonplace in many beekeeping regions across the world. In the present study, the occurrence of amino acid substitutions at residue L925 of the voltage-gate sodium channel-the pyrethroid target site-was studied in Varroa populations collected throughout Flanders, Belgium. Dose-response bioassays supported the involvement of the frequently observed L925V substitution in flumethrin resistance, resulting in a 12.64-fold increase of the LC₅₀ in a Varroa population mostly consisting of homozygous 925 V/V mites. With the presence of L925 substitutions in about four out of 10 screened apiaries, the use of pyrethroid-based varroacides in Flanders, including the recently released PolyVar® Yellow, should be carefully considered.

Sequence variation and regulatory variation in acetylcholinesterase genes contribute to insecticide resistance in different populations of Leptinotarsa decemlineata

Although insect herbivores are known to evolve resistance to insecticides through multiple genetic mechanisms, resistance in individual species has been assumed to follow the same mechanism. While both mutations in the target site insensitivity and increased amplification are known to contribute to insecticide resistance, little is known about the degree to which geographic populations of the same species differ at the target site in a response to insecticides. We tested structural (e.g., mutation profiles) and regulatory (e.g., the gene expression of Ldace1 and Ldace2, AChE activity) differences between two populations (Vermont, USA and Belchow, Poland) of the Colorado potato beetle, Leptinotarsa decemlineata in their resistance to two commonly used groups of insecticides, organophosphates, and carbamates. We established that Vermont beetles were more resistant to azinphos-methyl and carbaryl insecticides than Belchow beetles, despite a similar frequency of resistance-associated alleles (i.e., S291G) in the Ldace2 gene. However, the Vermont population had two additional amino acid replacements (G192S and F402Y) in the Ldace1 gene, which were absent in the Belchow population. Moreover, the Vermont population showed higher expression of Ldace1 and was less sensitive to AChE inhibition by azinphos-methyl oxon than the Belchow population. Therefore, the two populations have evolved different genetic mechanisms to adapt to organophosphate and carbamate insecticides.

Genome-Wide Identification and Validation of Target Genes Associated with Insecticide Treatment of the Green Peach Aphid, Myzus persicae

Next-generation sequencing and analyses of whole-genome transcripts can be used to identify genes and potential mechanisms that may be responsible for the development of resistance to insecticides. Such genes can be identified by isolating and sequencing high-quality messenger RNA and identifying differentially expressed genes (DEGs), and gene variants from insecticide-treated and untreated colonies of the Green peach aphid (GPA) or resistant and susceptible GPA populations. Datasets generated would reveal a set of genes whose expression may be associated with the insecticide treatment. The DEGs can then be validated using quantitative PCR assays.

Functional Characterization of Target Genes Associated with Insecticide Resistance of the Green Peach Aphid, Myzus persicae

Identifying genes responsive to insecticide treatment is the first step towards understanding the mechanism(s) of insecticide resistance and the development of effective insecticides against economic insect pests such as the Green peach aphid (GPA). Functional and Reverse Genetics approaches such as the RNA interference (RNAi) technology can be used to assess the possible involvement of genes whose expression is associated with an insecticide treatment. For GPA, this can be done by comparing the behavior and development of the insect following RNAi of a putative gene associated with insecticide treatment and exposure of the RNAi-treated insects to lethal doses of insecticides. In a case where knockdown of a gene or genes increases the susceptibility of RNAi-treated populations compared to controls, the target gene may have a direct role in the development of resistance to the insecticide or the gene may be involved in other metabolic processes that may be required for resilience against the insecticide.

Intraspecific Transcriptome Variation and Sex-Biased Expression in Anopheles arabiensis

The magnitude and functional patterns of intraspecific transcriptional variation in the anophelines, including those of sex-biased genes underlying sex-specific traits relevant for malaria transmission, remain understudied. As a result, how changes in expression levels drive adaptation in these species is poorly understood. We sequenced the female, male, and larval transcriptomes of three populations of Anopheles arabiensis from Burkina Faso. One-third of the genes were differentially expressed between populations, often involving insecticide resistance-related genes in a sample type-specific manner, and with the females showing the largest number of differentially expressed genes. At the genomic level, the X chromosome appears depleted of differentially expressed genes compared with the autosomes, chromosomes harboring inversions do not exhibit evidence for enrichment of such genes, and genes that are top contributors to functional enrichment patterns of population differentiation tend to be clustered in the genome. Further, the magnitude of variation for the sex expression ratio across populations did not substantially differ between male- and female-biased genes, except for some populations in which male-limited expressed genes showed more variation than their female counterparts. In fact, female-biased genes exhibited a larger level of interpopulation variation than male-biased genes, both when assayed in males and females. Beyond uncovering the extensive adaptive potential of transcriptional variation in An. Arabiensis, our findings suggest that the evolutionary rate of changes in expression levels on the X chromosome exceeds that on the autosomes, while pointing to female-biased genes as the most variable component of the An. Arabiensis transcriptome.

Pyrethroid metabolism by eleven Helicoverpa armigera P450s from the CYP6B and CYP9A subfamilies

Lepidopteran P450s of the CYP6B and CYP9A subfamilies are thought to play important roles in host plant adaptation and insecticide resistance. An increasing number of paralogs and orthologs with high levels of sequence identity have been found in these subfamilies by mining recent genome projects. However, the biochemical function of most of them remains unknown. A better understanding of the evolution of P450 genes and of the catalytic competence of the enzymes they encode is needed to facilitate studies of host plant use and insecticide resistance. Here, we focused on the full complement of CYP6B (4 genes) and CYP9A (7 genes) in the generalist herbivore, Helicoverpa armigera. These P450s were heterologously expressed in Sf9 cells and compared functionally. In vitro assays showed that all CYP6B and CYP9A P450s can metabolize esfenvalerate efficiently, except for the evolutionarily divergent CYP6B43. A new 2'-hydroxy-metabolite of esfenvalerate was identified and found to be the main metabolite produced by CYP9A12. All tested P450s showed only low induction responses to esfenvalerate. To put these results from H. armigera P450s in perspective, 158 complete CYP6B and 100 complete CYP9A genes from 34 ditrysian species were manually curated. The CYP9A subfamily was more widespread than the CYP6B subfamily and the latter showed dramatic gains and losses, with ten species lacking CYP6B genes. Two adjacent CYP6B loci were found on chromosome 21, with different fates during the evolution of Lepidoptera. The diversity and functional redundancy of CYP6B and CYP9A genes challenge resistance management and pest control strategies as many P450s are available to insects to cope with chemical stresses they encounter.

Evolution of Adaptive Variation in the Mosquito Culex pipiens: Multiple Independent Origins of Insecticide Resistance Mutations

Simple Summary

The application of insecticides to control insect pests can result in the evolution of resistance. Within a population, the individuals carrying the resistant mutation survive after insecticide application while the others die, leading to the spread of resistance. The finding that the same mutations confer insecticide resistance in different species or populations raises the question how often these mutations arise in natural populations. Resistant mutations can originate once within a population and then spread. Alternatively, multiple origins can occur within the same population or in different geographic areas. Here, we used a phylogenetic approach to investigate the origin of three mutations conferring resistance to diflubenzuron insecticide in Italian and Turkish populations of the mosquito Culex pipiens. Our results support a scenario of multiple origins of the mutations. Resistance is a major threat to mosquito control, and these findings help inform resistance management. At the same time, insecticide resistance is an informative model for studying the origin of adaptive variation. In the words of Rachel Carson, “If Darwin were alive today, he would be astounded and delighted by the impressive verification that his theories of the survival of the fittest are receiving from the insect world” (Silent Spring, 1962)

Abstract

Insecticide resistance is an informative model for studying the appearance of adaptive traits. Simultaneously, understanding how many times resistance mutations originate is essential to design effective resistance management. In the mosquito Culex pipiens, target–site resistance to the insecticide diflubenzuron (DFB) has been recently found in Italian and Turkish populations. Three point mutations confer it at the codon 1043 of the chitin synthase 1 gene (chs-1): I1043L, I1043M, and I1043F. Whether the resistant mutations originated independently from different susceptible alleles or sequentially from resistant alleles and whether resistant alleles from Italy and Turkey have originated once or multiple times remain unresolved. Here, we sequenced a fragment of the chs-1 gene carrying the resistant mutations and inferred the phylogenetic relationships among susceptible and resistant alleles. Confirming previous findings, we found the three mutations in Italy and the I1043M in Turkey. Notably, the I1043F was also found for the first time in Turkish samples, highlighting the need for extensive monitoring activities. Phylogenetic analyses are consistent with an independent origin of the I1043F, I1043M, and I1043L mutations from different susceptible alleles and with multiple independent origins of the Italian and Turkish I1043M and I1043F alleles.

Genome mapping coupled with CRISPR gene editing reveals a P450 gene confers avermectin resistance in the beet armyworm

The evolution of insecticide resistance represents a global constraint to agricultural production. Because of the extreme genetic diversity found in insects and the large numbers of genes involved in insecticide detoxification, better tools are needed to quickly identify and validate the involvement of putative resistance genes for improved monitoring, management, and countering of field-evolved insecticide resistance. The avermectins, emamectin benzoate (EB) and abamectin are relatively new pesticides with reduced environmental risk that target a wide number of insect pests, including the beet armyworm, Spodoptera exigua, an important global pest of many crops. Unfortunately, field resistance to avermectins recently evolved in the beet armyworm, threatening the sustainable use of this class of insecticides. Here, we report a high-quality chromosome-level assembly of the beet armyworm genome and use bulked segregant analysis (BSA) to identify the locus of avermectin resistance, which mapped on 15-16 Mbp of chromosome 17. Knockout of the CYP9A186 gene that maps within this region by CRISPR/Cas9 gene editing fully restored EB susceptibility, implicating this gene in avermectin resistance. Heterologous expression and in vitro functional assays further confirm that a natural substitution (F116V) found in the substrate recognition site 1 (SRS1) of the CYP9A186 protein results in enhanced metabolism of EB and abamectin. Hence, the combined approach of coupling gene editing with BSA allows for the rapid identification of metabolic resistance genes responsible for insecticide resistance, which is critical for effective monitoring and adaptive management of insecticide resistance.

Molecular bases of P450-mediated resistance to the neonicotinoid insecticide imidacloprid in the mosquito Ae. aegypti

Resistance to chemical insecticides including pyrethroids, the main insecticide class used against mosquitoes, has re-kindled interest in the use of neonicotinoids. In this context, the present study aimed to characterize the molecular basis of neonicotinoid resistance in the mosquito Aedes aegypti. Resistance mechanisms were studied by combining transcriptomic and genomic data obtained from a laboratory strain selected at the larval stage after 30 generations of exposure to imidacloprid (Imida-R line). After thirty generations of selection, larvae of the Imida-R line showed an 8-fold increased resistance to imidacloprid and a significant cross-tolerance to the pyrethroids permethrin and deltamethrin. Cross-resistance to pyrethroids was only observed in adults when larvae were previously exposed to imidacloprid suggesting a low but inducible expression of resistance alleles at the adult stage. Resistance of the Imida-R line was associated with a slower larval development time in females. Multiple detoxification enzymes were over-transcribed in larvae in association with resistance including the P450s CYP6BB2, CYP9M9 and CYP6M11 previously associated with pyrethroid resistance. Some of them together with their redox partner NADPH P450 reductase were also affected by non-synonymous mutations associated with resistance. Combining genomic and transcriptomic data allowed identifying promoter variations associated with the up-regulation of CYP6BB2 in the resistant line. Overall, these data confirm the key role of P450s in neonicotinoid resistance in Ae. aegypti and their potential to confer cross-resistance to pyrethroids, raising concerns about the use of neonicotinoids for resistance management in this mosquito species.

High-resolution genetic mapping reveals cis-regulatory and copy number variation in loci associated with cytochrome P450-mediated detoxification in a generalist arthropod pest

Chemical control strategies are driving the evolution of pesticide resistance in pest populations. Understanding the genetic mechanisms of these evolutionary processes is of crucial importance to develop sustainable resistance management strategies. The acaricide pyflubumide is one of the most recently developed mitochondrial complex II inhibitors with a new mode of action that specifically targets spider mite pests. In this study, we characterize the molecular basis of pyflubumide resistance in a highly resistant population of the spider mite Tetranychus urticae. Classical genetic crosses indicated that pyflubumide resistance was incompletely recessive and controlled by more than one gene. To identify resistance loci, we crossed the resistant population to a highly susceptible T. urticae inbred strain and propagated resulting populations with and without pyflubumide exposure for multiple generations in an experimental evolution set-up. High-resolution genetic mapping by a bulked segregant analysis approach led to the identification of three quantitative trait loci (QTL) linked to pyflubumide resistance. Two QTLs were found on the first chromosome and centered on the cytochrome P450 CYP392A16 and a cluster of CYP392E6-8 genes. Comparative transcriptomics revealed a consistent overexpression of CYP392A16 and CYP392E8 in the experimental populations that were selected for pyflubumide resistance. We further corroborated the involvement of CYP392A16 in resistance by in vitro functional expression and metabolism studies. Collectively, these experiments uncovered that CYP392A16 N-demethylates the toxic carboxamide form of pyflubumide to a non-toxic compound. A third QTL coincided with cytochrome P450 reductase (CPR), a vital component of cytochrome P450 metabolism. We show here that the resistant population harbors three gene copies of CPR and that this copy number variation is associated with higher mRNA abundance. Together, we provide evidence for detoxification of pyflubumide by cytochrome P450s that is likely synergized by gene amplification of CPR.

Our understanding of the causal genetic variants that drive the evolution of quantitative traits, such as polygenic pesticide resistance, remains very limited. Here, we followed a high-resolution genetic mapping approach to localize the genetic variants that cause pyflubumide resistance in the two-spotted spider mite Tetranychus urticae. Three well-supported QTL were uncovered and pointed towards a major role for cytochrome P450-mediated detoxification. Cis-regulatory variation for cytochrome P450s was observed, and in vitro cytochrome P450 experiments showed that pyflubumide was metabolized into a non-toxic derivate. A third QTL centered on cytochrome P450 reductase (CPR), which is required for cytochrome P450 activity, and is amplified in pyflubumide resistant populations. Our results indicate that pyflubumide resistance is mediated by cytochrome P450 detoxification that is enhanced by gene amplification at the CPR locus.

Resistance in the Genus Spodoptera: Key Insect Detoxification Genes

The genus Spodoptera (Lepidoptera: Noctuidae) includes species that are among the most important crop pests in the world. These polyphagous species are able to feed on many plants, including corn, rice and cotton. In addition to their ability to adapt to toxic compounds produced by plants, they have developed resistance to the chemical insecticides used for their control. One of the main mechanisms developed by insects to become resistant involves detoxification enzymes. In this review, we illustrate some examples of the role of major families of detoxification enzymes such as cytochromes P450, carboxyl/cholinesterases, glutathione S-transferases (GST) and transporters such as ATP-binding cassette (ABC) transporters in insecticide resistance. We compare available data for four species, Spodoptera exigua, S. frugiperda, S. littoralis and S. litura. Molecular mechanisms underlying the involvement of these genes in resistance will be described, including the duplication of the CYP9A cluster, over-expression of GST epsilon or point mutations in acetylcholinesterase and ABCC2. This review is not intended to be exhaustive but to highlight the key roles of certain genes.

Comparative analysis of mite genomes reveals positive selection for diet adaptation

Diet is a powerful evolutionary force for species adaptation and diversification. Acari is one of the most abundant clades of Arachnida, exhibiting diverse dietary types, while the underlying genetic adaptive mechanisms are not fully understood. Based on comparative analyses of 15 Acari genomes, we found genetic bases for three specialized diets. Herbivores experienced stronger selection pressure than other groups; the olfactory genes and gene families involving metabolizing toxins showed strong adaptive signals. Genes and gene families related to anticoagulation, detoxification, and haemoglobin digestion were found to be under strong selection pressure or significantly expanded in the blood-feeding species. Lipid metabolism genes have a faster evolutionary rate and been subjected to greater selection pressures in fat-feeding species; one positively selected site in the fatty-acid amide hydrolases 2 gene was identified. Our research provides a new perspective for the evolution of Acari and offers potential target loci for novel pesticide development.

Selenium mediated host plant-mite conflict: defense and adaptation

BACKGROUND

Selenium has shown effectiveness in protecting plants from herbivores. However, some insects have evolved adaptability to selenium.

RESULTS

Selenium accumulation in host plants protected them against spider mite feeding. Selenium showed toxic effects on spider mites by reducing growth and interfering with reproduction. After 40 generations on selenium-rich plants, a Tetranychus cinnabarinus strain (Tc-Se) developed adaptability to selenium, with an increased rate of population growth and enhanced ability for selenium metabolism. The high expression of two genes (GSTd07 and SPS1) in the selenium metabolism pathway might be involved in selenium metabolism in spider mites. After GSTd07 and SPS1 were silenced, the selenium adaptability decreased. Recombinant GSTd07 protein promoted the reaction between sodium selenite and glutathione (GSH) and increased the production of sodium selenite metabolites. The results indicated that GSTd07 was involved in the first step of selenium metabolism.

CONCLUSION

Plants can resist spider mite feeding by accumulating selenium. Spider mites subjected to long-term selenium exposure can adapt to selenium by increasing the expression of key genes involved in selenium metabolism. These results elucidate the mechanism of the interaction between mites and host plants mediated by selenium. This study of the interaction between selenium-mediated host plants and spider mites may lead to the development of new and less toxic methods for the prevention and control of spider mites. © 2021 Society of Chemical Industry.

Learnings from over a decade of increasing pesticide resistance in the redlegged earth mite, Halotydeus destructor (Tucker)

BACKGROUND

The redlegged earth mite, Halotydeus destructor (Tucker), is a destructive and economically important pest of winter grain crops and pastures in Australia. It is largely controlled by pesticides, but this mite has evolved resistance to pyrethroid and organophosphate chemicals. A national Resistance Management Strategy has been developed for pro-active management to delay further resistance evolution, though its success is reliant on a detailed understanding of the incidence, patterns of spread, current distribution and the nature of resistance in the field. Here, we report on a long-term resistance surveillance programme undertaken between 2006 and 2019 informed by resistance risk forecasting.

RESULTS

By mapping the Australian distribution of resistance through time, we show that resistance is present across three Australian states and covers more than 3000 km. This current range includes a recently identified population exhibiting organophosphate resistance representing the most easterly location of resistance in H. destructor. Using field history information, we identify associations for the first time between crop management practices employed by farmers and the presence of pyrethroid resistance. Management strategies that could minimize the risk of further resistance include limiting local spread of resistance through farm hygiene practices, crop rotations and reducing pesticide usage.

CONCLUSION

This study highlights the challenges of resistance in H. destructor but also indicates how quantitative resistance risk analysis can be developed to target field surveillance and delay further resistance. The management strategies highlighted in this study can help maintain the effectiveness of control options but will depend on farmer engagement and adoption. © 2021 Society of Chemical Industry.

Genetic analysis and screening of pyrethroid resistance mutations in Varroa destructor populations from Turkey

Varroa destructor is the most common ectoparasite of the Western honey bee (Apis mellifera L.) worldwide and poses a serious threat to bee health. Synthetic acaricides, particularly pyrethroids, are frequently used to control Varroa mites. However, long-term and repeated use of synthetic pyrethroids has led to the development of resistance. In this study, we report on the presence of resistance mutations in the voltage-gated sodium channel in V. destructor populations from Turkish beekeeping areas. Two resistance mutations, L925V and L925I, that were previously associated with pyrethroid resistance, were found in more than 75% of the populations. A general correlation between the presence of mutations and the history of acaricide usage was observed for the sampled hives. In addition, we show there is only a low genetic distance among the sampled V. destructor populations, based on the analysis of three mitochondrial genes: cytochrome b (cytb), ATP synthase subunit 6 (atp6), and cytochrome c oxidase subunit III (cox3). Revealing the presence and geographical distribution of pyrethroid resistance mutations in V. destructor populations from Turkish apiaries will contribute to create more effective mite management programmes.

The role of detoxification enzymes in the susceptibility of Brevipalpus californicus exposed to acaricide and insecticide mixtures

The intense spraying of pesticides to control different arthropod pests has resulted in negative side effects for the management of pests. It was previously discovered that exposure to non-acaricidal insecticides alone or in a mixture, results in lower efficiency of the acaricide spirodiclofen used for Brevipalpus spp. control. We investigate here whether the induced expression of detoxification enzymes by non-lethal insecticides may antagonize spirodiclofen toxicity. Brevipalpus californicus mites exposed to the insecticide phosmet alone or in combination with spirodiclofen showed increased activity of P450 monooxygenases (P450s). No antagonistic effects in mite mortality were observed by the combination of phosmet and spirodiclofen. On the other hand, mites exposed to the insecticide imidacloprid alone or in combination with spirodiclofen showed an increase in the activity of P450s, carboxylcholinesterases (CCE), and glutathione-S-transferases (GST). An antagonistic effect on mite mortality was observed when mites were exposed to the LC₂₅ of spirodiclofen combined with the field rate treatment of imidacloprid. The addition of PBO (a P450 monooxygenase inhibitor) to the mixture of spirodiclofen and imidacloprid resulted in a synergistic effect over mite mortality but the addition of DEM (a GST inhibitor) resulted in an antagonist effect. Taken together, this study showed that the combination of imidacloprid with spirodiclofen is antagonistic for the control of B. californicus, and this results from the induction of detoxification enzymes, such as P450s, CCE, and GST. The use of inhibitors highlights the role of these enzymes in the antagonism of the mixture.

Overexpression of Multiple UDP-Glycosyltransferase Genes Involved in Sulfoxaflor Resistance in Aphis gossypii Glover

UDP-glycosyltransferases (UGTs) are major phase II enzymes involved in the metabolic detoxification of xenobiotics. In this study, two UGT-inhibitors, 5-nitrouracil and sulfinpyrazone, significantly increased sulfoxaflor toxicity against sulfoxaflor-resistant (Sul-R) Aphis gossypii, whereas there were no synergistic effects in susceptible (Sus) A. gossypii. The activity of UGTs in the Sul-R strain was significantly higher (1.35-fold) than that in the Sus strain. Further, gene expression determination demonstrated that 11 of 23 UGT genes were significantly upregulated (1.40- to 5.46-fold) in the Sul-R strain, among which the expression levels of UGT350A2, UGT351A4, UGT350B2, UGT342C2, and UGT343C2 could be induced by sulfoxaflor. Additionally, knockdown of UGT350A2, UGT351A4, UGT350B2, and UGT343C2 using RNA interference (RNAi) significantly increased sensitivity (1.57- to 1.76-fold) to sulfoxaflor in the Sul-R strain. These results suggested that UGTs might be involved in sulfoxaflor resistance in A. gossypii. These findings will facilitate further work to validate the functional roles of these UGT genes in sulfoxaflor resistance.