QTL mapping suggests that both cytochrome P450-mediated detoxification and target-site resistance are involved in fenbutatin oxide resistance in Tetranychus urticae

The cytochrome P450 subfamilies CYP392A and CYP392D are key players in acaricide metabolism in Tetranychusurticae

The spider mite Tetranychus urticae is a major agricultural pest with a global distribution, extremely diverse host range and a remarkable ability to develop resistance to a wide variety of acaricides. P450 mono-oxygenases have been frequently associated with resistance development in this species. In particular enzymes of the CYP392A-subfamily were shown to metabolize a number of key acaricides, including abamectin, amitraz, fenpyroximate and the active metabolite of pyflubumide. However, transcriptomic studies comparing highly resistant and susceptible populations have often revealed high expression of members of the CYP392D-subfamily, but these have been only poorly studied. Here, we conducted a meta-analysis of gene expression data of 20 populations and identified two key enzymes of this family, CYP392D2 and CYP392D8, whose expression is associated with resistance. We subsequently functionally expressed these enzymes, together with CYP392A11 and CYP392A16 as known metabolizers, and compared their potential to accept a wide diversity of acaricides as substrate. This study overall confirms previous discovered substrates for CYP392A11 and CYP392A16, but also reveals unreported metabolic activity towards new acaricides. These include carbaryl, chlorpyrifos and etoxazole for CYP392A16 and carbaryl, chlorpyrifos and NNI-0711-NH pyflubumide for CYP392A11. For the newly studied CYP392D-family, we show that CYP392D2 metabolizes pyridaben, fenpyroximate, etoxazole and chlorpyrifos, while CYP392D8 metabolizes carbaryl, fenazaquin and tebufenpyrad. Last, we observed that both CYP392A- and CYP392D-subfamily enzymes activate chlorpyrifos to its corresponding oxon. Our study indicates that there is both overlap and specificity in the activity of A- and D-subfamily enzymes against acaricides and model substrates. With the recent advent of highly efficient CRISPR/Cas9 gene editing protocols in T. urticae, the way is now paved to conduct further genetic experiments revealing and quantifying the role of these enzymes in the resistance phenotype in field populations.

Cell Painting unravels insecticidal modes of action on Spodoptera frugiperda insect cells

The "Cell Painting" technology utilizes multiplexed fluorescent staining of various cell organelles, to produce high-content microscopy images of cells for multidimensional phenotype assessment. The phenotypic profiles extracted from those images can be analyzed upon perturbations with biologically active molecules to annotate the mode of action or biological activity by comparison with reference profiles of already known mechanisms of action, ultimately enabling the determination of on-target and off-target effects. This approach is already described in various human cell cultures, the most commonly used being the U2OS cell line, yet allows broad applications in additional areas of chemical-biological research. Here we describe for the first time the application and adaptation of Cell Painting to an insect cell line, the Sf9 cells from Spodoptera frugiperda. By adjusting image acquisition and analysis models, specific phenotypic profiles were obtained in a dose-dependent manner for 20 reference compounds, including representatives for the most relevant insecticidal modes of action categories (nerve & muscle, respiration and growth & development). Through a dimensionality-reduction method, both calculations of phenotypic half maximal inhibition concentration (IC50) values as well as similarity analysis of the obtained profiles by hierarchical clustering were performed. By Cell Painting effects on the phenotype could be obtained at higher sensitivity than in other assay formats, such as cytotoxicity assessments. More importantly, these analyses provide insight into mechanistic determinants of biological activity. Compounds with similar modes of action showed a high degree of proximity in a hierarchical clustering analysis while being distinct from actives with an unrelated mode of action. In essence, we provide strong evidence on the impact of Cell Painting mechanistic understanding of insecticides with regards to determinants of efficacy and safety utilizing an insect cell model system.

Functional analysis of SDR112C1 associated with fenpropathrin tolerance in Tetranychus cinnabarinus (Boisduval)

Short-chain dehydrogenases/reductases (SDRs) are ubiquitously distributed across diverse organisms and play pivotal roles in the growth, as well as endogenous and exogenous metabolism of various substances, including drugs. The expression levels of SDR genes are reportedly upregulated in the fenpropathrin (FEN)-resistant (FeR) strain of Tetranychus cinnabarinus. However, the functions of these SDR genes in acaricide tolerance remain elusive. In this study, the activity of SDRs was found to be significantly higher (2.26-fold) in the FeR strain compared to the susceptible strain (SS) of T. cinnabarinus. A specific upregulated SDR gene, named SDR112C1, exhibited significant overexpression (3.13-fold) in the FeR population compared with that in the SS population. Furthermore, the expression of SDR112C1 showed a significant increase in the response to FEN induction. Additionally, knockdown of the SDR112C1 gene resulted in decreased SDR activity and reduced mite viability against FEN. Importantly, heterologous expression and in vitro incubation assays confirmed that recombinant SDR112C1 could effectively deplete FEN. Moreover, the overexpression of the SDR112C1 gene in Drosophila melanogaster significantly decreased the toxicity of FEN to transgenic fruit flies. These findings suggest that the overexpression of SDR SDR112C1 is a crucial factor contributing to FEN tolerance in T. cinnabarinus. This discovery not only enhances our understanding of SDR-mediated acaricide tolerance but also introduces a new family of detoxification enzymes to consider in practice, beyond cytochrome P450s, carboxyl/choline esterases and glutathione S-transferases.

SYNCAS: Efficient CRISPR/Cas9 gene-editing in difficult to transform arthropods

The genome editing technique CRISPR/Cas9 has led to major advancements in many research fields and this state-of-the-art tool has proven its use in genetic studies for various arthropods. However, most transformation protocols rely on microinjection of CRISPR/Cas9 components into embryos, a method which is challenging for many species. Alternatively, injections can be performed on adult females, but transformation efficiencies can be very low as was shown for the two-spotted spider mite, Tetranychus urticae, a minute but important chelicerate pest on many crops. In this study, we explored different CRISPR/Cas9 formulations to optimize a maternal injection protocol for T. urticae. We observed a strong synergy between branched amphipathic peptide capsules and saponins, resulting in a significant increase of CRISPR/Cas9 knock-out efficiency, exceeding 20%. This CRISPR/Cas9 formulation, termed SYNCAS, was used to knock-out different T. urticae genes - phytoene desaturase, CYP384A1 and Antennapedia - but also allowed to develop a co-CRISPR strategy and facilitated the generation of T. urticae knock-in mutants. In addition, SYNCAS was successfully applied to knock-out white and white-like genes in the western flower thrips, Frankliniella occidentalis. The SYNCAS method allows routine genome editing in these species and can be a game changer for genetic research in other hard to transform arthropods.

Contrasting roles of cytochrome P450s in amitraz and chlorfenapyr resistance in the crop pest Tetranychus urticae

The molecular mechanisms of amitraz and chlorfenapyr resistance remain only poorly understood for major agricultural pests and vectors of human diseases. This study focusses on a multi-resistant field strain of the crop pest Tetranychus urticae, which could be readily selected in the laboratory to high levels of amitraz and chlorfenapyr resistance. Toxicity experiments using tralopyril, the active toxophore of chlorfenapyr, suggested decreased activation as a likely mechanism underlying resistance. Starting from the same parental strain, transcriptome profiling revealed that a cluster of detoxifying genes was upregulated after amitraz selection, but unexpectedly downregulated after chlorfenapyr selection. Further functional validation associated the upregulation of CYP392A16 with amitraz metabolism and the downregulation of CYP392D8 with reduced activation of chlorfenapyr to tralopyril. Genetic mapping (QTL analysis by BSA) was conducted in an attempt to unravel the genetic mechanisms of expression variation and resistance. This revealed that chlorfenapyr resistance was associated with a single QTL, while 3 QTLs were uncovered for amitraz resistance. Together with the observed contrasting gene expression patterns, we argue that transcriptional regulators most likely underly the distinct expression profiles associated with resistance, but these await further functional validation.

The complex II resistance mutation H258Y in succinate dehydrogenase subunit B causes fitness penalties associated with mitochondrial respiratory deficiency

BACKGROUND

The acaricides cyflumetofen, cyenopyrafen and pyflubumide inhibit the mitochondrial electron transport chain at complex II [succinate dehydrogenase (SDH) complex]. A target site mutation H258Y was recently discovered in a resistant strain of the spider mite pest Tetranychus urticae. H258Y causes strong cross-resistance between cyenopyrafen and pyflubumide, but not cyflumetofen. In fungal pests, fitness costs associated with substitutions at the corresponding H258 position that confer resistance to fungicidal SDH inhibitors have not been uncovered. Here, we used H258 and Y258 near-isogenic lines of T. urticae to quantify potential pleiotropic fitness effects on mite physiology.

RESULTS

The H258Y mutation was not associated with consistent significant changes of single generation life history traits and fertility life table parameters. In contrast, proportional Sanger sequencing and droplet digital polymerase chain reaction showed that the frequency of the resistant Y258 allele decreased when replicated 50:50 Y258:H258 experimentally evolving populations were maintained in an acaricide-free environment for approximately 12 generations. Using in vitro assays with mitochondrial extracts from resistant (Y258) and susceptible (H258) lines, we identified a significantly reduced SDH activity (48% lower activity) and a slightly enhanced combined complex I and III activity (18% higher activity) in the Y258 lines.

CONCLUSION

Our findings suggest that the H258Y mutation is associated with a high fitness cost in the spider mite T. urticae. Importantly, while it is the most common approach, it is clear that only comparing life history traits and life table fecundity does not allow to reliably estimate fitness costs of target site mutations in natural pest populations. © 2023 Society of Chemical Industry.

The host plant strongly modulates acaricide resistance levels to mitochondrial complex II inhibitors in a multi-resistant field population of Tetranychus urticae

The two-spotted spider mite Tetranychus urticae is a polyphagous pest with an extraordinary ability to develop acaricide resistance. Here, we characterize the resistance mechanisms in a T. urticae population (VR-BE) collected from a Belgian tomato greenhouse, where the grower was unsuccessful in chemically controlling the mite population resulting in crop loss. Upon arrival in the laboratory, the VR-BE population was established both on bean and tomato plants as hosts. Toxicity bioassays on both populations confirmed that the population was highly multi-resistant, recording resistance to 12 out of 13 compounds tested from various mode of action groups. DNA sequencing revealed the presence of multiple target-site resistance mutations, but these could not explain resistance to all compounds. In addition, striking differences in toxicity for six acaricides were observed between the populations on bean and tomato. The highest difference was recorded for the complex II inhibitors cyenopyrafen and cyflumetofen, which were 4.4 and 3.3-fold less toxic for VR-BE mites on tomato versus bean. PBO synergism bioassays suggested increased P450 based detoxification contribute to the host-dependent toxicity. Given the involvement of increased detoxification, we subsequently determined genome-wide gene expression levels of VR-BE on both hosts, in comparison to a reference susceptible population, revealing overexpression of a large set of detoxification genes in VR-BE on both hosts compared to the reference. In addition, a number of mainly detoxification genes with higher expression in VR-BE on tomato compared to bean was identified, including several cytochrome P450s. Together, our work suggests that multi-resistant field populations can accumulate a striking number of target-site resistance mutations. We also show that the host plant can have a profound effect on the P450-associated resistance levels to cyenopyrafen and cyflumetofen.

A review of the molecular mechanisms of acaricide resistance in mites and ticks

The Arachnida subclass of Acari comprises many harmful pests that threaten agriculture as well as animal health, including herbivorous spider mites, the bee parasite Varroa, the poultry mite Dermanyssus and several species of ticks. Especially in agriculture, acaricides are often used intensively to minimize the damage they inflict, promoting the development of resistance. Beneficial predatory mites used in biological control are also subjected to acaricide selection in the field. The development and use of new genetic and genomic tools such as genome and transcriptome sequencing, bulked segregant analysis (QTL mapping), and reverse genetics via RNAi or CRISPR/Cas9, have greatly increased our understanding of the molecular genetic mechanisms of resistance in Acari, especially in the spider mite Tetranychus urticae which emerged as a model species. These new techniques allowed to uncover and validate new resistance mutations in a larger range of species. In addition, they provided an impetus to start elucidating more challenging questions on mechanisms of gene regulation of detoxification associated with resistance.

Increased metabolism in combination with the novel cytochrome b target-site mutation L258F confers cross-resistance between the Qo inhibitors acequinocyl and bifenazate in Tetranychus urticae

Acequinocyl and bifenazate are potent acaricides acting at the Q_o site of complex III of the electron transport chain, but frequent applications of these acaricides have led to the development of resistance in spider mites. Target-site resistance caused by mutations in the conserved cd1- and ef-helices of the Q_o pocket of cytochrome b has been elucidated as the main resistance mechanism. We therefore monitored Q_o pocket mutations in European field populations of Tetranychus urticae and uncovered a new mutation, L258F. The role of this mutation was validated by revealing patterns of maternal inheritance and by the independently replicated introgression in an unrelated susceptible genetic background. However, the parental strain exhibited higher resistance levels than conferred by the mutation alone in isogenic lines, especially for acequinocyl, implying the involvement of strong additional resistance mechanisms. This was confirmed by revealing a polygenic inheritance pattern with classical genetic crosses and via synergism experiments. Therefore, a genome-wide expression analysis was conducted that identified a number of highly overexpressed detoxification genes, including many P450s. Functional expression revealed that the P450 CYP392A11 can metabolize bifenazate by hydroxylation of the ring structure. In conclusion, the novel cytochrome b target-site mutation L258F was uncovered in a recently collected field strain and its role in acequinocyl and bifenazate resistance was validated. However, the high level of resistance in this strain is most likely caused by a combination of target-site resistance and P450-based increased detoxification, potentially acting in synergism.

Intraspecific diversity in the mechanisms underlying abamectin resistance in a cosmopolitan pest

Pesticide resistance relies on a myriad of mechanisms, ranging from single mutations to a complex and polygenic architecture, and it involves mechanisms such as target-site insensitivity, metabolic detoxification, or a combination of these, with either additive or synergistic effects. Several resistance mechanisms against abamectin, a macrocyclic lactone widely used in crop protection, have been reported in the cosmopolitan pest Tetranychus urticae. However, it has been shown that a single mechanism cannot account for the high levels of abamectin resistance found across different mite populations. Here, we used experimental evolution combined with bulked segregant analyses to map quantitative trait loci (QTL) associated with abamectin resistance in two genetically unrelated populations of T. urticae. In these two independent QTL mapping experiments, three and four QTLs were identified, of which three were shared between experiments. Shared QTLs contained genes encoding subunits of the glutamate-gated chloride channel (GluCl) and harboured previously reported mutations, including G314D in GluCl1 and G326E in GluCl3, but also novel resistance candidate loci, including DNA helicases and chemosensory receptors. Surprisingly, the fourth QTL, present only in only one of the experiments and thus unique for one resistant parental line, revealed a non-functional variant of GluCl2, suggesting gene knock-out as resistance mechanism. Our study uncovers the complex basis of abamectin resistance, and it highlights the intraspecific diversity of genetic mechanisms underlying resistance in a cosmopolitan pest.

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.

High-Resolution Genetic Mapping Combined with Transcriptome Profiling Reveals That Both Target-Site Resistance and Increased Detoxification Confer Resistance to the Pyrethroid Bifenthrin in the Spider Mite Tetranychus urticae

Pyrethroids are widely applied insecticides in agriculture, but their frequent use has provoked many cases of resistance, in which mutations in the voltage-gated sodium channel (VGSC), the pyrethroid target-site, were shown to play a major role. However, for the spider mite Tetranychus urticae, it has also been shown that increased detoxification contributes to resistance against the pyrethroid bifenthrin. Here, we performed QTL-mapping to identify the genomic loci underlying bifenthrin resistance in T. urticae. Two loci on chromosome 1 were identified, with the VGSC gene being located near the second QTL and harboring the well-known L1024V mutation. In addition, the presence of an L925M mutation in the VGSC of a highly bifenthrin-resistant strain and its loss in its derived, susceptible, inbred line indicated the importance of target-site mutations in bifenthrin resistance. Further, RNAseq experiments revealed that genes encoding detoxification enzymes, including carboxyl/choline esterases (CCEs), cytochrome P450 monooxygenases and UDP-glycosyl transferases (UGTs), were overexpressed in resistant strains. Toxicity bioassays with bifenthrin (ester pyrethroid) and etofenprox (non-ester pyrethroid) also indicated a possible role for CCEs in bifenthrin resistance. A selection of CCEs and UGTs were therefore functionally expressed, and CCEinc18 was shown to metabolize bifenthrin, while teturUGT10 could glycosylate bifenthrin-alcohol. To conclude, our findings suggest that both target-site and metabolic mechanisms underlie bifenthrin resistance in T. urticae, and these might synergize high levels of resistance.

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.