High-resolution QTL mapping in Tetranychus urticae reveals acaricide-specific responses and common target-site resistance after selection by different METI-I acaricides

Overexpression of a nuclear receptor HR96 contributes to spirodiclofen susceptibility in Panonychus citri (McGregor)

The citrus red mite, Panonychus citri, is one of the most notorious and devastating citrus pests around the world that has developed resistance to multiple chemical acaricides. In previous research, we found that spirodiclofen-resistant is related to overexpression of P450, CCE, and ABC transporter genes in P. citri. However, the regulatory mechanisms of these detoxification genes are still elusive. This study identified all hormone receptor 96 genes of P. citri. 8 PcHR96 genes contained highly conserved domains. The expression profiles showed that PcHR96h was significantly upregulated in spirodiclofen resistant strain and after exposure to spirodiclofen. RNA interference of PcHR96h decreased expression of detoxification genes and increased spirodiclofen susceptibility in P. citri. Furthermore, molecular docking, heterologous expression, and drug affinity responsive target stability demonstrated that PcHR96h can interact with spirodiclofen in vitro. Our research results indicate that PcHR96h plays an important role in regulating spirodiclofen susceptibility and provides theoretical support for the resistance management of P. citri.

Monitoring detoxification enzyme levels and resistance of Tetranychus urticae against some METI-group chemicals in Türkiye cotton fields

Acaricides used against Tetranychus urticae Koch, 1836 (Acari: Tetranychidae) in cotton fields cause control failure over time. To determine the resistance status of T. urticae populations to tebufenpyrad and bifenazate, different populations collected from Aydın (AYD), Adana (ADA), Şanlıurfa (SAN), and Diyarbakır (DIY) provinces of Türkiye, between 2019 and 2020, were subjected to diagnostic dose bioassays. Firstly, the spider mites were eliminated with a discriminating dose. Afterwards, LC50 and LC90 of the remaining populations were determined and the ten highest resistant populations were selected. The highest phenotypic resistance to bifenazate was observed in AYD4 and DIY2 (LC50 57.14 mg L- 1 with 85.01-fold and LC50 30.15 mg L- 1with 44.86-fold, respectively), while the lowest phenotypic resistance was found in SAN6 (LC50 1.5 mg L- 1; 2.28-fold). Considering the phenotypic resistance to tebufenpyrad, the highest resistance was found in AYD4 population (LC50 96.81 mg L- 1; 12.92-fold), while the lowest - in DIY28 population (LC50 21.23 mg L- 1; 2.83-fold). In pharmacokinetic studies, the ADA16 population was compared with the sensitive German Susceptible Strain population and it was determined that carboxylesterase activity was statistically higher (1.46 ± 0.04 nmol/min/mg protein enzyme activation 2.70-fold). The highest activation of glutathione S-transferase was detected in ADA16 (1.49 ± 0.01 nmol/min/mg protein; 2.32-fold). No mutations were found in PSST (METI 1), the point mutation site for tebufenpyrad, and Cytb (METI 3), the point mutation site for bifenazate. In terms of phenotypic resistance, bifenazate was found to be moderately resistant in two populations (85.01 and 44.86-fold), while tebufenpyrad was moderately resistant in one population (12.92-fold). This study showed that both acaricides are still effective against T. urticae populations.

A high-throughput multilocus-amplicon sequencing panel to monitor insecticide resistance in fall armyworm (FAW) Spodoptera frugiperda (Lepidoptera: Noctuidae)

BACKGROUND

Fall armyworm (FAW), Spodoptera frugiperda, is a highly polyphagous crop pest that has spread over the world rapidly and invaded Australia in 2020. Globally, FAW has been reported to be resistant to several insecticides permitted in Australia. Timely resistance diagnosis is critical for integrated pest management-based control of FAW in Australia.

RESULTS

We developed a multi-amplicon panel by next-generation sequencing (multiamplicon-seq) to identify known insecticide resistance mutations in Australian FAW with high throughput and low cost. The panel included nine known mutations causing insecticide resistance in FAW and four gene mutations causing insecticide resistance in several insect species, not yet reported in FAW. We sequenced 36 plates (96-well) in one MiSeq flow cell with easy sequencing library preparation. We found that Australian FAW carried a very high proportion of the F290V mutation in the acetylcholinesterase (AChE) gene that causes resistance to organophosphate and carbamate insecticides. Furthermore, FAW has a GABA-activated chloride channel mutation, A301Q in the RDL gene. The sequencing-based platform provided evidence of a duplication in the AChE gene. Here several single nucleotide polymorphisms (SNPs) within the 476-bp amplicon of the AChE gene demonstrated 100% heterozygosity across samples and some individuals carried two haplotypes with the F290V mutation.

CONCLUSION

Molecular surveillance by multiamplicon-seq will increase capacity for early detection and future resistance monitoring in highly dispersed Australian FAW. It can provide timely resistance information and has the potential to play an important role in the resistance management of FAW. © 2023 Society of Chemical Industry.

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.

CYP4CL2 Confers Metabolic Resistance to Pyridaben in the Citrus Pest Mite Panonychus citri

The citrus red mite Panonychus citri has developed strong resistance to acaricides. Cytochrome P450 monooxygenases (P450s) can detoxify pesticides and are involved in pesticide resistance in many insects. Here, a pyridaben-resistant P. citri strain showed cross-resistance to cyenopyrafen, bifenazate, fenpyroximate, and tolfenpyrad. Piperonyl butoxide, a P450 inhibitor, significantly increased the toxicity of pyridaben to resistant (Pyr_Rs) and susceptible (Pyr_Control) P. citri strains. P450 activity was significantly higher in Pyr_Rs than in Pyr_Control. Analyses of RNA-Seq data identified a P450 gene (CYP4CL2) that is potentially involved in pyridaben resistance. Consistently, it was up-regulated in two field-derived resistant populations (CQ_WZ and CQ_TN). RNA interference-mediated knockdown of CYP4CL2 significantly decreased the pyridaben resistance in P. citri. Transgenic Drosophila melanogaster expressing CYP4CL2 showed increased pyridaben resistance. Molecular docking analysis showed that pyridaben could bind to several amino acids at substrate recognition sites in CYP4CL2. These findings shed light on P450-mediated pyridaben resistance in pest mites.

A glutamate-gated chloride channel as the mite-specific target-site of dicofol and other diphenylcarbinol acaricides

Dicofol has been widely used to control phytophagous mites. Although dicofol is chemically related to DDT, its mode of action has remained elusive. Here, we mapped dicofol resistance in the spider mite Tetranychus urticae to two genomic regions. Each region harbored a glutamate-gated chloride channel (GluCl) gene that contained a mutation-G314D or G326E-known to confer resistance against the unrelated acaricide abamectin. Using electrophysiology assays we showed that dicofol and other diphenylcarbinol acaricides-bromopropylate and chlorobenzilate-induce persistent currents in Xenopus oocytes expressing wild-type T. urticae GluCl3 receptors and potentiate glutamate responses. In contrast, the G326E substitution abolished the agonistic activity of all three compounds. Assays with the wild-type Drosophila GluClα revealed that this receptor was unresponsive to dicofol. Homology modeling combined with ligand-docking confirmed the specificity of electrophysiology assays. Altogether, this work elucidates the mode of action of diphenylcarbinols as mite-specific agonists of GluCl.

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 detoxification pathway initiated by a nuclear receptor TcHR96h in Tetranychus cinnabarinus (Boisduval)

Understanding the mechanism of detoxification initiation in arthropods after pesticide exposure is crucial. Although the identity of transcription factors that induce and regulate the expression of detoxification genes in response to pesticides is beginning to emerge, whether transcription factors directly interact with xenobiotics is unclear. The findings of this study revealed that a nuclear hormone receptor, Tetranychus cinnabarinus hormone receptor (HR) TcHR96h, regulates the overexpression of the detoxification gene TcGSTm02, which is involved in cyflumetofen resistance. The nuclear translocation of TcHR96h increased after cyflumetofen exposure, suggesting direct binding with cyflumetofen. The direct binding of TcHR96h and cyflumetofen was supported by several independent proteomic assays that quantify interactions with small molecules. Together, this study proposes a model for the initiation of xenobiotic detoxification in a polyphagous agricultural pest. These insights not only provide a better understanding of the mechanisms of xenobiotic detoxification and metabolism in arthropods, but also are crucial in understanding adaptation in polyphagous herbivores.

A nuclear receptor HR96-related gene underlies large trans-driven differences in detoxification gene expression in a generalist herbivore

The role, magnitude, and molecular nature of trans-driven expression variation underlying the upregulation of detoxification genes in pesticide resistant arthropod populations has remained enigmatic. In this study, we performed expression quantitative trait locus (eQTL) mapping (n = 458) between a pesticide resistant and a susceptible strain of the generalist herbivore and crop pest Tetranychus urticae. We found that a single trans eQTL hotspot controlled large differences in the expression of a subset of genes in different detoxification gene families, as well as other genes associated with host plant use. As established by additional genetic approaches including RNAi gene knockdown, a duplicated gene with a nuclear hormone receptor HR96-related ligand-binding domain was identified as causal for the expression differences between strains. The presence of a large family of HR96-related genes in T. urticae may enable modular control of detoxification and host plant use genes, facilitating this species' known and rapid evolution to diverse pesticides and host plants.

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.

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.

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.

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.

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.

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

The organotin acaricide fenbutatin oxide (FBO) - an inhibitor of mitochondrial ATP-synthase - has been one of the most extensively used acaricides for the control of spider mites, and is still in use today. Resistance against FBO has evolved in many regions around the world but only few studies have investigated the molecular and genetic mechanisms of resistance to organotin acaricides. Here, we found that FBO resistance is polygenic in two genetically distant, highly resistant strains of the spider mite Tetranychus urticae, MAR-AB and MR-VL. To identify the loci underlying FBO resistance, two independent bulked segregant analysis (BSA) based QTL mapping experiments, BSA MAR-AB and BSA MR-VL, were performed. Two QTLs on chromosome 1 were associated with FBO resistance in each mapping experiment. At the second QTL of BSA MAR-AB, several cytochrome P450 monooxygenase (CYP) genes were located, including CYP392E4, CYP392E6 and CYP392E11, the latter being overexpressed in MAR-AB. Synergism tests further implied a role for CYPs in FBO resistance. Subunit c of mitochondrial ATP-synthase was located near the first QTL of both mapping experiments and harbored a unique V89A mutation enriched in the resistant parents and selected BSA populations. Marker-assisted introgression into a susceptible strain demonstrated a moderate but significant effect of the V89A mutation on toxicity of organotin acaricides. The impact of the mutation on organotin inhibition of ATP synthase was also functionally confirmed by ATPase assays on mitochondrial preparations. To conclude, our findings suggest that FBO resistance in the spider mite T. urticae is a complex interplay between CYP-mediated detoxification and target-site resistance.

freqpcr: Estimation of population allele frequency using qPCR ΔΔCq measures from bulk samples

PCR techniques, both quantitative (qPCR) and nonquantitative, have been used to estimate the frequency of a specific allele in a population. However, the labour required to sample numerous individuals and subsequently handle each sample renders the quantification of rare mutations (e.g., pesticide resistance gene mutations at the early stages of resistance development) challenging. Meanwhile, pooling DNA from multiple individuals as a "bulk sample" combined with qPCR may reduce handling costs. The qPCR output for a bulk sample, however, contains uncertainty owing to variations in DNA yields from each individual, in addition to measurement errors. In this study, we have developed a statistical model to estimate the frequency of the specific allele and its confidence interval when the sample allele frequencies are obtained in the form of ΔΔCq in the qPCR analyses on multiple bulk samples collected from a population. We assumed a gamma distribution as the individual DNA yield and developed an R package for parameter estimation, which was verified using real DNA samples from acaricide-resistant spider mites, as well as a numerical simulation. Our model resulted in unbiased point estimates of the allele frequency compared with simple averaging of the ΔΔCq values. The confidence intervals suggest that dividing the bulk samples into more parts will improve precision if the total number of individuals is equal; however, if the cost of PCR analysis is higher than that of sampling, increasing the total number and pooling them into a few bulk samples may also yield comparable precision.

Cross resistance, inheritance and fitness advantage of cyetpyrafen resistance in two-spotted spider mite, Tetranychus urticae

Cyetpyrafen belonging to mitochondrial electron transport inhibitors of complex II (METI II) has been widely applied to manage pest mites in China. To investigate the adaption of Tetranychus urticae in the evolution of cyetpyrafen resistance, a study of cross resistance, mode of inheritance and fitness comparison of resistance using indoor cyetpyrafen resistant strain (resistance ratio, RR > 2, 000-fold) was executed. Cyet-R showed serious cross resistance to cyenopyrafen (>2500-fold) and cyflumetofen (~190-fold). The number of resistant genes was evaluated via chi-square (χ²) test and the concentration-response curve regarding goodness-of-fit between observed and the expected mortality. The LC₅₀s of F_1RS (Cyet-R^♀ × Tu-YN^♂) and F_1SR (Tu-YN^♀ × Cyet-R^♂) were 3126.30 mg/L and 2743.97 mg/L, respectively, without significance, suggesting autosomal inheritance. The degree of dominance (D) values of F_1RS and F_1SR ranged between 0 and 1, revealing an incompletely dominant inheritance in the tested population of Tetranychus urticae. Plots of concentration-response data for the orthogonal backcross and reverse backcross progenies showed a significant deviation from the expected lines, pointing out a polygenic inheritance. Besides, lifetable analysis showed a fitness advantage of Cyet-R with a significantly decreased adult preadult period and significantly increased total fecundity. This study suggested that cyetpyrafen resistance against T. urticae was inherited as autosomal, incompletely dominant and multigenetic and characterized with serious cross resistance and fitness advantage. Therefore, rational application and preventive strategy should be considered to sustain the efficacy of cyetpyrafen against T. urticae.

Predicting the genomic resolution of bulk segregant analysis

Bulk segregant analysis is a technique for identifying the genetic loci that underlie phenotypic trait differences. The basic approach is to compare two pools of individuals from the opposing tails of the phenotypic distribution, sampled from an interbred population. Each pool is sequenced and scanned for alleles that show divergent frequencies between the pools, indicating potential association with the observed trait differences. Bulk segregant analysis has already been successfully applied to the mapping of various quantitative trait loci in organisms ranging from yeast to maize. However, these studies have typically suffered from rather low mapping resolution, and we still lack a detailed understanding of how this resolution is affected by experimental parameters. Here, we use coalescence theory to calculate the expected genomic resolution of bulk segregant analysis for a simple monogenic trait. We first show that in an idealized interbreeding population of infinite size, the expected length of the mapped region is inversely proportional to the recombination rate, the number of generations of interbreeding, and the number of genomes sampled, as intuitively expected. In a finite population, coalescence events in the genealogy of the sample reduce the number of potentially informative recombination events during interbreeding, thereby increasing the length of the mapped region. This is incorporated into our model by an effective population size parameter that specifies the pairwise coalescence rate of the interbreeding population. The mapping resolution predicted by our calculations closely matches numerical simulations and is surprisingly robust to moderate levels of contamination of the segregant pools with alternative alleles. Furthermore, we show that the approach can easily be extended to modifications of the crossing scheme. Our framework will allow researchers to predict the expected power of their mapping experiments, and to evaluate how their experimental design could be tuned to optimize mapping resolution.

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.

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.

Multiple TaqMan qPCR and droplet digital PCR (ddPCR) diagnostics for pesticide resistance monitoring and management, in the major agricultural pest Tetranychus urticae

BACKGROUND

Decisions on which pesticide to use in agriculture are expected to become more difficult, as the number of available chemicals is decreasing. For Tetranychus urticae (T. urticae), a major pest for which a number of candidate markers for pesticide resistance are in place, molecular diagnostics could support decision-making for the rational use of acaricides.

RESULTS

A suite of 12 TaqMan qPCR assays [G314D (GluCl1), G326E, I321T (GluCl3), G119S, F331W (Ace-1), H92R (PSST), L1024V, F1538I (VGSC), I1017F (CHS1), G126S, S141F, P262T (cytb)], were validated against Sanger-sequencing, and subsequently adapted for use with the ddPCR technology. The concordance correlation coefficient between the actual and ddPCR measured mutant allelic frequencies was 0.995 (95% CI = 0.991-0.998), and no systematic, proportional, or random differences were detected. The achieved Limit of Detection (LoD) was 0.1% (detection of one mutant in a background of 999 wild type mites). The ddPCR assay panel was then assessed in terms of agreement with phenotypic resistance, through a pilot application in field populations from Crete, with strong correlation and thus predictive and diagnostic value of the molecular assays in some cases (e.g., etoxazole and abamectin resistance). Molecular diagnostics were able to capture incipient resistance that was otherwise missed by phenotypic bioassays. The molecular and phenotypic resistance screening of T. urticae field populations from Crete, revealed both multi-resistant and susceptible populations.

CONCLUSION

The highly sensitive T. urticae molecular diagnostic platforms developed in this study could prove a valuable tool for pesticide resistance management. © 2021 Society of Chemical Industry.

Natural tetramic acids elicit multiple inhibitory actions against mitochondrial machineries presiding over oxidative phosphorylation

The mitochondrial machineries presiding over ATP synthesis via oxidative phosphorylation are promising druggable targets. Fusaramin, a 3-acyl tetramic acid isolated from Fusarium concentricum FKI-7550, is an inhibitor of oxidative phosphorylation in Saccharomyces cerevisiae mitochondria, although its target has yet to be identified. Fusaramin significantly interfered with [3H]ADP uptake by yeast mitochondria at the concentration range inhibiting oxidative phosphorylation. A photoreactive fusaramin derivative (pFS-5) specifically labeled voltage-dependent anion channel 1 (VDAC1), which facilitates trafficking of ADP/ATP across the outer mitochondrial membrane. These results strongly suggest that the inhibition of oxidative phosphorylation by fusaramin is predominantly attributable to the impairment of VDAC1 functions. Fusaramin also inhibited FoF1-ATP synthase and ubiquinol-cytochrome c oxidoreductase (complex III) at concentrations higher than those required for the VDAC inhibition. Considering that other tetramic acid derivatives are reported to inhibit FoF1-ATP synthase and complex III, natural tetramic acids were found to elicit multiple inhibitory actions against mitochondrial machineries.

Comparing the efficiency of RNAi after feeding and injection of dsRNA in spider mites

Pesticide resistance in spider mites drives the development of acaricides with novel mode of action, which could benefit from RNAi as a screening tool in search of new molecular targets. RNAi via oral delivery of dsRNA has been frequently reported in spider mites, but injection of dsRNA is rarely reported. We compare here the efficiency of oral delivery versus injection of dsRNA in female adult mites. When comparing silencing efficiency, oral delivery of dsRNAs silenced 40.6 ± 8.9% of CPR, 63.8 ± 6.9% of CHMP2A, and 37.7 ± 5.7% of CHMP3 genes. Similar silencing efficiencies were found for injection (48.6 ± 3.7% of CPR, 70.2 ± 4.1% of CHMP2A, 59.8 ± 2.2% of CHMP3), but with much lower quantities of dsRNAs. Oral delivery of dsRNA failed to silence the expression of the CHMP4B gene, but this could be accomplished by injection of dsRNA (23.1 ± 1.0%). When scoring the phenotypic effects of silencing, both oral delivery and injection of CHMP2A- and CHMP3-dsRNA influenced the locomotion speed of mites significantly. For CPR, silencing could only be accomplished by dsRNA injection, not by feeding. CPR silencing significantly impacted the toxicity of a typical acaricide, pyridaben, as the susceptibility of mites raised 2.75-fold. Last, injection of Eya-dsRNA in adults produced transgenerational phenotypic effects on 3.59% of offspring, as quantified by an observed deviation in eye development, while oral delivery of Eya-dsRNA did not. In conclusion, injection of dsRNA is superior to oral delivery in silencing the expression of the selected genes in this study and could be considered the method of choice to study gene function in reverse genetic approaches.

Co-occurrence of subunit B and C mutations in respiratory complex II confers high resistance levels to pyflubumide and cyenopyrafen in the two-spotted spider mite Tetranychus urticae (Acari: Tetranychidae)

BACKGROUND

Pyflubumide and cyenopyrafen are respiratory complex II (complex II) inhibitors. Previous quantitative trait locus analyses suggested associations of I260V and S56L in complex II subunit B (B-I260V) and subunit C (C-S56L) with pyflubumide and cyenopyrafen resistance, respectively, in Tetranychus urticae. However, although resistant strains had been selected separately by these acaricides, all strains were homozygous for both B-I260V and C-S56L. Hence, the effects of each mutation on resistance development remain unclear.

RESULTS

We established strains homozygous for B-I260V with C-S56 (B-I260V_I260V/C-S56_S56) and for C-S56L with B-I260 (B-I260_I260/C-S56L_S56L). High resistance levels (LC₅₀  > 1000 mg L^-1 ) to pyflubumide and cyenopyrafen was not conferred by B-I260V or C-S56L alone. Next, we prepared intermixed strains by crossing B-I260V_I260V/C-S56_S56 and B-I260_I260/C-S56L_S56L. Selection of the intermixed strains by either acaricide caused very high resistance levels (LC₅₀  ≥ 10 000 mg L^-1 ) to both acaricides and fixed both mutations. Allele-selected recoupling of the mutations without acaricide selection also conferred very high resistance levels to both acaricides in the intermixed strains. Unlike these, B-I260V or C-S56L alone conferred very high and high resistance levels to cyflumetofen, respectively.

CONCLUSION

We conclude that the effect of individual mutations characteristically varies among complex II inhibitors. Moreover, very high resistance levels to pyflubumide and cyenopyrafen is conferred by the co-occurrence of B-I260V and C-S56L mutations, which alone have limited effects on resistance level.

Comparative analysis of the detoxification gene inventory of four major Spodoptera pest species in response to xenobiotics

The genus Spodoptera (Lepidoptera: Noctuidae) comprises some of the most polyphagous and destructive agricultural pests worldwide. The success of many species of this genus is due to their striking abilities to adapt to a broad range of host plants. Superfamilies of detoxification genes play a crucial role in the adaption to overcome plant defense mechanisms mediated by numerous secondary metabolites and toxins. Over the past decade, a substantial amount of expression data in Spodoptera larvae was produced for those genes in response to xenobiotics such as plant secondary metabolites, but also insecticide exposure. However, this information is scattered throughout the literature and in most cases does not allow to clearly identify candidate genes involved in host-plant adaptation and insecticide resistance. In the present review, we analyzed and compiled information on close to 600 pairs of inducers (xenobiotics) and induced genes from four main Spodoptera species: S. exigua, S. frugiperda, S. littoralis and S. litura. The cytochrome P450 monooxygenases (P450s; encoded by CYP genes) were the most upregulated detoxification genes across the literature for all four species. Most of the data was provided from studies on S. litura, followed by S. exigua, S. frugiperda and S. littoralis. We examined whether these detoxification genes were reported for larval survival under xenobiotic challenge in forward and reverse genetic studies. We further analyzed whether biochemical assays were carried out showing the ability of corresponding enzymes and transporters to breakdown and excrete xenobiotics, respectively. This revealed a clear disparity between species and the lack of genetic and biochemical information in S. frugiperda. Finally, we discussed the biological importance of detoxification genes for this genus and propose a workflow to study the involvement of these enzymes in an ecological and agricultural context.

Whole genome sequencing and bulked segregant analysis suggest a new mechanism of amitraz resistance in the citrus red mite, Panonychus citri (Acari: Tetranychidae)

BACKGROUND

Amitraz is a broad-spectrum insecticide/acaricide for the control of aphids, psyllids, ticks and mites. Current evidence suggests that ticks and phytophagous mites have developed strong resistance to amitraz. Previous studies have shown that multiple mechanisms are associated with amitraz resistance in ticks, but very few reports have involved Panonychus citri. We therefore used whole genome sequencing and bulked segregant analysis (BSA) to identify the mechanism underlying P. citri's resistance to amitraz.

RESULTS

High-quality assembly of the whole P. citri genome was completed, resulting in a genome of approximately 83.97 Mb and a contig N50 of approximately 1.81 Mb. Gene structure predictions revealed 11 577 genes, of which 10 940 genes were annotated. Trait-associated regions in the genome were mapped with bulked segregant analysis and 38 candidate SNPs were obtained, of which T752C had the strongest correlation with the resistant trait, located at the 5' untranslated region (UTR) of the β-2R adrenergic-like octopamine receptor gene. The mutation resulted in the formation of a short hairpin loop structure in mRNA and gene expression was down-regulated by more than 50% in the amitraz-resistant strain. Validation of the T752C mutation in field populations of P. citri found that the correlation between the resistance ratio and the base mutation was 94.40%.

CONCLUSION

Our results suggest that this 5' UTR mutation of the β-2R octopamine receptor gene, confers amitraz resistance in P. citri. This discovery provides a new explanation for the mechanism of pest resistance: base mutations in the 5' untranslated region of target gene may regulate the susceptibility of pests to pesticides.

Toxicity of Varronia curassavica Jacq. Essential Oil to Two Arthropod Pests and Their Natural Enemy

This study aimed to evaluate the toxicity of the Varronia curassavica Jacq. essential oil to two cosmopolitan and polyphagous pest species, the two-spotted spider mite (Tetranychus urticae Koch) and the green aphid (Myzus persicae Sulzer). Additionally, we tested the essential oil toxicity to a generalist predatory insect, the green lacewing Ceraeochrysa cubana Hagen. The treatments consisted of four essential oil concentrations (0.25, 0.5, 0.75 and 1.0%) and one control (Tween® 80 + water). At 0.75% concentration, the V. curassavica essential oil caused the highest mortalities for both pests. The concentration of 1.0% of the essential oil interfered mainly in the rates of oviposition and egg hatching and caused a decrease in the instantaneous population growth rate (r_i) of T. urticae. Survival of the C. cubana was not affected by the highest tested concentration of the essential oil (1%). Tested concentrations V. curassavica essential oil were toxic to T. urticae and M. persicae, but not to the predator C. cubana.

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.

Who Is Dermanyssus gallinae? Genetic Structure of Populations and Critical Synthesis of the Current Knowledge

Despite the economic and animal welfare importance of the Poultry Red Mite Dermanyssus gallinae, its genetic structure has been studied in a scattered way so far. The prophylaxis and control of such a globally distributed ectoparasite can be significantly improved by understanding its genetic population structure (composition in species and intraspecific variants). The present study aims to establish a rigorous framework for characterizing the neutral genetic structure of D. gallinae based on a literature review combined with an integrative analysis of the data available in GenBank on population-level nucleotide sequence diversity supplemented by a new dataset. The integrative analysis was conducted on sequence data extracted from GenBank coupled with new sequences of two fragments of the mitochondrial gene encoding Cytochrome Oxidase I (CO1) as well as of an intron of the nuclear gene encoding Tropomyosin (Tpm) from several PRM populations sampled from European poultry farms. Emphasis was placed on using the mitochondrial gene encoding CO1 on which the main universal region of DNA barcoding in animals is located. The species D. gallinae sensu lato is a species complex, encompassing at least two cryptic species, i.e., not distinguishable by morphological characters: D. gallinae sensu stricto and D. gallinae L1. Only D. gallinae s.s. has been recorded among the populations sampled in poultry farms worldwide. Current knowledge suggests they are structured in three mitochondrial groups (haplogroups A, B, and C). Haplogroup A is cosmopolitan, and the other two present slightly contrasted distributions (B rather in the northern part of Europe, C most frequently found in the southern part). Recent data indicate that a dynamic geographic expansion of haplogroup C is underway in Europe. Our results also show that NUMT (nuclear mitochondrial DNA) pseudogenes have generated artifactual groups (haplogroups E and F). It is important to exclude these artifact groups from future analyses to avoid confusion. We provide an operational framework that will promote consistency in the analysis of subsequent results using the CO1 fragment and recommendations for future analyses.

Identification and characterization of striking multiple-insecticide resistance in a Tetranychus urticae field population from Greece

BACKGROUND

Tetranychus urticae is a notorious crop pest with a worldwide distribution that has developed resistance to a wide range of acaricides. Here, we investigated the resistance levels of a T. urticae population collected from an ornamental greenhouse in Peloponnese, Greece, and analyzed its resistance mechanisms at the molecular level.

RESULTS

Toxicological assays showed resistance against compounds with different modes of action, with resistance ratios of: 89-fold for abamectin; > 1000-fold for clofentezine; > 5000-fold for etoxazole; 27-fold for fenpyroximate and pyridaben; 20- and 36-fold for spirodiclofen and spirotetramat, respectively; and 116- and > 500-fold for cyenopyrafen and cyflumetofen, respectively. Bioassays with synergists indicated the involvement of detoxification enzymes in resistance to abamectin, but not to cyflumetofen and spirodiclofen. RNA sequencing (RNA-seq) analysis showed significant over-expression of several genes encoding detoxification enzymes such as cytochrome P450 monooxygenases and UDP-glycosyltransferases, which have been previously associated with acaricide resistance. Known target-site resistance mutations were identified in acetyl-choline esterase, chitin synthase 1 and NDUFS7/psst, but putative novel resistance mutations were also discovered in targets such as glutamate-gated chloride channel subunit 3. Interestingly, target-site resistance mutations against pyrethroids or bifenazate were not identified, possibly indicating a recent reduced selection pressure in Greece, as well as a possible opportunity to rotate these chemistries.

CONCLUSION

We identified and characterized a striking case of multiple acaricide resistance in a field population of T. urticae. Exceptionally strong resistance phenotypes, with accumulation of multiple resistance mutations and over-expression of P450s and other detoxification genes in the same field population are reported.

Fenpyroximate resistance in Iranian populations of the European red mite Panonychus ulmi (Acari: Tetranychidae)

The European red mite, Panonychus ulmi (Koch), is one of the most important apple orchard pests worldwide. Fenpyroximate, a mitochondrial electron transport inhibitor of complex I (METI-I), is a commonly used acaricide to control this pest. In this study, we determined fenpyroximate resistance levels for 11 P. ulmi populations from Iran and a spirodiclofen-resistant strain from Germany (PSR-TK). The LC₅₀ values ranged between 121.8 and 5713.9 mg a.i. L^-1 and the highest resistance ratio (RR) was 47-fold for the Padena population. PBO, TPP and DEM synergist ratios (SRs) were the highest for the PSR-TK (SR = 6.7), Shahin Dej (SR = 6.1) and Semirom3 (SR = 3.6) populations, respectively. In vitro enzyme activity measurements also showed that there was a higher glutathione S-transferases (GSTs) activity in the PSR-TK and Shahin Dej population compared to the most susceptible populations, whereas the esterase and P450 monooxygenase activity were not significantly higher in the resistant populations. Last, we screened all populations for the presence of two mutations previously associated with METI-I resistance in spider mites but none of these mutations could be detected. To conclude, moderate to high levels of fenpyroximate resistance were observed in P. ulmi populations from Iran, with increased detoxification most likely underlying fenpyroximate resistance.

Genetic analysis and screening of detoxification-related genes in an amitraz-resistant strain of Panonychus citri

Panonychus citri (McGregor) is the most common pest in citrus-producing regions. Special low-toxicity acaricides, such as spirocyclic tetronic acids and mite growth inhibitors, have been used for a long time in China. However, pesticide resistance in mites is a growing problem due to the lack of new acaricide development. Wide-spectrum insecticides, such as amitraz have gained acceptance among fruit growers. An amitraz-resistant strain of P. citri was obtained by indoor screening to examine field resistance monitoring of mites to acaricides and to explore the resistant mechanism of mites against amitraz. The amitraz-resistant strain of P. citri had an LC50 value of 2361.45 mg l-1. The resistance ratio was 81.35 times higher in the resistant strain of P. citri compared with the sensitive strain. Crossing experiments between the sensitive and resistant strains of P. citri were conducted, resulting in a D value of 0.11 for F1 SS♀×RS♂ and 0.06 for F1 RS♀×SS♂. Reciprocal cross experiments showed that the dose-mortality curves for the F1 generations coincided, indicating that the resistance trait was not affected by cytoplasmic inheritance. The dose-expected response relationship was evaluated in the backcross generation and a significant difference was observed compared with the actual value. The above results indicate that the inheritance of resistance trait was incompletely dominant, governed by polygenes on the chromosome. Synergism studies demonstrated that cytochrome P450s and esterase may play important roles in the detoxification of amitraz. Based on differential gene analysis, 23 metabolism-related genes of P. citri were identified, consistent with the results of synergism studies. Real-time PCR verification implied that P450s, ABC transporters, and acetylcholinesterase might influence the detoxification of amitraz by P. citri. These results provide the genetic and molecular foundation for the management of pest mite resistance.

Co-Expression of a Homologous Cytochrome P450 Reductase Is Required for In Vivo Validation of the Tetranychus urticae CYP392A16-Based Abamectin Resistance in Drosophila

Simple Summary

The two-spotted spider mite, Tetranychus urticae, is one of the most damaging agricultural pests worldwide, feeding on over 1100 plant species and causing extensive damage to several crops. Chemical acaricides remain the most widely used strategy to control this pest. However, T. urticae has developed significant resistance to numerous acaricide compounds, due to certain features of mite biology and extensive acaricide applications that lead to the selection of resistant pests and subsequently the emergence of resistant populations. Several molecular/genetic mechanisms may contribute to these highly resistant phenotypes. Such mechanisms frequently involve expression of P450 detoxification enzymes, which act together with a partner protein named cytochrome P450 reductase (CPR). In this study, we investigated the potential of a mite P450 enzyme, CYP392A16, to confer resistance to the acaricide abamectin in vivo, when expressed in tissues of the model fruit fly Drosophila melanogaster. We confirmed that expression of this enzyme contributes to abamectin resistance in the fruit fly model, but only when a homologous mite CPR is co-expressed. Our findings indicate that the Drosophila model system can be engineered to facilitate validation of the candidate mite P450s, in order to elucidate resistance mechanisms and their underlying interactions.

Abstract

Overexpression of the cytochrome P450 monooxygenase CYP392A16 has been previously associated with abamectin resistance using transcriptional analysis in the two-spotted spider mite Tetranychus urticae, an important pest species worldwide; however, this association has not been functionally validated in vivo despite the demonstrated ability of CYP392A16 to metabolize abamectin in vitro. We expressed CYP392A16 in vivo via a Gal4 transcription activator protein/Upstream Activating Sequence (GAL4/UAS) system in Drosophila melanogaster flies, driving expression with detoxification tissue-specific drivers. We demonstrated that CYP392A16 expression confers statistically significant abamectin resistance in toxicity bioassays in Drosophila only when its homologous redox partner, cytochrome P450 reductase (TuCPR), is co-expressed in transgenic flies. Our study shows that the Drosophila model can be further improved, to facilitate the functional analysis of insecticide resistance mechanisms acting alone or in combination.

QTL mapping using microsatellite linkage reveals target-site mutations associated with high levels of resistance against three mitochondrial complex II inhibitors in Tetranychus urticae

The acaricides cyflumetofen, cyenopyrafen, and pyflubumide act as inhibitors of the mitochondrial electron transport system at complex II (succinate dehydrogenase; SDH), a new mode of action in arthropods. The development and mechanisms of low-level resistance against cyenopyrafen and cyflumetofen have been previously reported in Tetranychus urticae. In the present study, we investigated high levels of resistance against three SDH inhibitors in T. urticae field populations and clarify the genetic basis of resistance using quantitative trait locus (QTL) analysis. First, we constructed a microsatellite linkage map comprising 64 markers assembled into three linkage groups (LGs) with total length of 683.8 cM and average marker spacing of 11.03 cM. We then used the linkage map to perform QTL mapping, and identified significant QTLs contributing to resistance to cyflumetofen (one QTL on LG1), cyenopyrafen (one QTL on LG3), and pyflubumide (two QTLs on LG1 and LG3). The QTL peaks on LG1 for cyflumetofen and pyflubumide overlapped and included the SdhB locus. For cyenopyrafen resistance, the QTLs on LG3 included the SdhC locus. For cyflumetofen resistance, we found an I260T mutation in SdhB. For pyflubumide and cyenopyrafen resistance, we detected I260V and S56L substitutions in SdhB and SdhC, respectively, by direct sequencing. Both I260 in SdhB and S56 in SdhC were present in highly conserved regions of the ubiquinone binding site formed at the interface among SdhB, SdhC, and SdhD. Mutations at these positions have been implicated in resistance against fungicides that act as Sdh inhibitors in various pathogens. Therefore, we consider these mutations to be target-site resistance mutations for these acaricidal SDH inhibitors.

lincRNA_Tc13743.2-miR-133-5p-TcGSTm02 regulation pathway mediates cyflumetofen resistance in Tetranychus cinnabarinus

Differential expression of metabolic detoxification enzymes is an important mechanism involved in pesticide/acaricide resistance of mite pests. The competing endogenous RNA hypothesis offers a new opportunity to investigate post-transcriptional regulation of those genes. In this study, 4454 long non-coding RNAs were identified in the carmine spider mite Tetranychus cinnabarinus by transcriptome sequencing. Software-based predictions indicated that a long intergenic non-coding RNA, (lincRNA)_Tc13743.2 and a detoxification enzyme gene, TcGSTm02, both contained a microRNA (miR-133-5p) response element. Over-expression of lincRNA_Tc13743.2 and TcGSTm02 were detected in a cyflumetofen-resistant T. cinnabarinus strain (CyR), whereas down-regulation of miR-133-5p was observed in this strain. Conversely, up-regulation of miR-133-5p could inhibit TcGSTm02 expression levels, and both lincRNA_Tc13743.2 and TcGSTm02 were significantly enriched in miR-133-5p biotin-avidin pull-down assays. RNA-binding protein immunoprecipitation assay showed that lincRNA_Tc13743.2 and TcGSTm02 bound to a silencing complex containing miR-133-5p. Moreover, a luciferase reporter assay based on a human cell line revealed that over-expression of lincRNA_Tc13743.2 could significantly reduce the inhibition exerted by miR-133-5p through the TcGSTm02 3'UTR. In addition, co-localization of lincRNA_Tc13743.2 and miR-133-5p was detected using fluorescence in situ hybridization, suggesting that lincRNA_Tc13743.2 interacts directly with miR-133-5p in spider mites. More importantly, silencing the expression of lincRNA_Tc13743.2 significantly reduced the expression levels of TcGSTm02 and increased the sensitivity of spider mites to cyflumetofen. Our data show that lincRNA_Tc13743.2 up-regulates TcGSTm02 expression by competing for miR-133-5p binding, demonstrating that a lincRNA_Tc13743.2-miR-133-5p-TcGSTm02 pathway mediates cyflumetofen resistance in mites.

Using CRISPR/Cas9 genome modification to understand the genetic basis of insecticide resistance: Drosophila and beyond

Chemical insecticides are a major tool for the control of many of the world's most damaging arthropod pests. However, their intensive application is often associated with the emergence of resistance, sometimes with serious implications for sustainable pest control. To mitigate failure of insecticide-based control tools, the mechanisms by which insects have evolved resistance must be elucidated. This includes both identification and functional characterization of putative resistance genes and/or mutations. Research on this topic has been greatly facilitated by using powerful genetic model insects like Drosophila melanogaster, and more recently by advances in genome modification technology, notably CRISPR/Cas9. Here, we present the advances that have been made through the application of genome modification technology in insecticide resistance research. The majority of the work conducted in the field to date has made use of genetic tools and resources available in D. melanogaster. This has greatly enhanced our understanding of resistance mechanisms, especially those mediated by insensitivity of the pesticide target-site. We discuss this progress for a series of different insecticide targets, but also report a number of unsuccessful or inconclusive attempts that highlight some inherent limitations of using Drosophila to characterize resistance mechanisms identified in arthropod pests. We also discuss an experimental framework that may circumvent current limitations while retaining the genetic versatility and robustness that Drosophila has to offer. Finally, we describe examples of direct CRISPR/Cas9 use in non-model pest species, an approach that will likely find much wider application in the near future.

Significance and interpretation of molecular diagnostics for insecticide resistance management of agricultural pests

Insecticide resistant pests become increasingly difficult to control in current day agriculture. Because of environmental and health concerns, the insecticide portfolio to combat agricultural pests is gradually decreasing. It is therefore crucial to make rational decisions on insecticide use to assure effective resistance management. However, resistance monitoring programs that inform on pest susceptibility and resistance are not yet common practice in agriculture. Molecular markers of resistance that are turned into convenient diagnostic tools are urgently needed and will only increase in importance. This review investigates which factors determine the strength, diagnostic value, and success of a diagnostic marker, and in which cases recent technical advances might provide new opportunities for decision making in an operational meaningful way.

RNA interference of NADPH-cytochrome P450 reductase increases susceptibilities to multiple acaricides in Tetranychus urticae

The two-spotted spider mite, Tetranychus urticae, is a polyphagous pest feeding on over 1100 plant species, including numerous highly valued economic crops. The control of T. urticae largely depends on the use of acaricides, which leads to pervasive development of acaricide resistance. Cytochrome P450-mediated metabolic detoxification is one of the major mechanisms of acaricide resistance in T. urticae. NADPH-cytochrome P450 reductase (CPR) plays as a crucial co-factor protein that donates electron(s) to microsomal cytochrome P450s to complete their catalytic cycle. This study seeks to understand the involvement of CPR/P450 in acaricide resistance in T. urticae. The full-length cDNA sequence of T. urticae's CPR (TuCPR) was cloned and characterized. TuCPR was ubiquitously transcribed in different life stages of T. urticae and the highest transcription was observed in the nymph and adult stages. TuCPR was constitutively over-expressed in six acaricide resistant populations compared to a susceptible one. TuCPR transcriptional expression was also induced by multiple acaricides in a time-dependent manner. Down-regulation of TuCPR via RNA interference (RNAi) in T. urticae led to reduced enzymatic activities of TuCPR and cytochrome P450s, as well as a reduction of resistance to multiple acaricides, abamectin, bifenthrin, and fenpyroximate. The outcome of this study highlights CPR as a potential novel target for eco-friendly control of T. urticae and other related plant-feeding pests.

Targeted mutagenesis using CRISPR-Cas9 in the chelicerate herbivore Tetranychus urticae

The use of CRISPR-Cas9 has revolutionized functional genetic work in many organisms, including more and more insect species. However, successful gene editing or genetic transformation has not yet been reported for chelicerates, the second largest group of terrestrial animals. Within this group, some mite and tick species are economically very important for agriculture and human health, and the availability of a gene-editing tool would be a significant advancement for the field. Here, we report on the use of CRISPR-Cas9 in the spider mite Tetranychus urticae. The ovary of virgin adult females was injected with a mix of Cas9 and sgRNAs targeting the phytoene desaturase gene. Natural mutants of this laterally transferred gene have previously shown an easy-to-score albino phenotype. Albino sons of injected virgin females were mated with wild-type females, and two independent transformed lines where created and further characterized. Albinism inherited as a recessive monogenic trait. Sequencing of the complete target-gene of both lines revealed two different lesions at expected locations near the PAM site in the target-gene. Both lines did not genetically complement each other in dedicated crosses, nor when crossed to a reference albino strain with a known genetic defect in the same gene. In conclusion, two independent mutagenesis events were induced in the spider mite T. urticae using CRISPR-Cas9, hereby providing proof-of-concept that CRISPR-Cas9 can be used to create gene knockouts in mites.

Molecular and genetic analysis of resistance to METI-I acaricides in Iranian populations of the citrus red mite Panonychus citri

The citrus red mite, Panonychus citri, is a major pest on citrus all around the world. Mitochondrial Electron Transport Inhibitors of complex I (METI-I) acaricides such as fenpyroximate have been used extensively to control P. citri populations, which resulted in multiple reports of METI-I resistant populations in the field. In this study, biochemical and molecular mechanisms of fenpyroximate resistance were investigated in P. citri. Seven populations were collected from Northern provinces of Iran. Resistance ratios were determined and reached up to 75-fold in comparison to a fenpyroximate susceptible population. Cross-resistance to two additional METI-I acaricides, pyridaben and tebufenpyrad, was detected. PBO synergism experiments, in vivo enzyme assays and gene expression analysis suggest a minor involvement of cytochrome P450 monooxygenases in fenpyroximate resistance, which is in contrast with many reported cases for the closely related Tetranychus urticae. Next, we determined the frequency of a well-known mutation in the target-site of METI-Is, the PSST subunit, associated with METI-I resistance. Indeed, the H92R substitution was detected in a highly fenpyroximate resistant P. citri population. Additionally, a new amino acid substitution at a conserved site in the PSST subunit was detected, A94V, with higher allele frequencies in a moderately resistant population. Marker-assisted back-crossing in a susceptible background confirmed the potential involvement of the newly discovered A94V mutation in fenpyroximate resistance. However, introduction of the A94V mutation in the PSST homologue of D. melanogaster using CRISPR-Cas9 did not result in fenpyroximate resistant flies. In addition, differences in binding curves between METI-Is and complex I measured directly, in isolated transgenic and wildtype mitochondria preparations, could not be found.

Analysis of the Physicochemical Properties of Acaricides Based on Lipinski's Rule of Five

This study aimed to explore the similarity between the physical and chemical properties of different acaricides, determine whether Lipinski's Rule of Five (RO5) used in the design of drug molecules is suitable for screening acaricides, and provide methods for selection of new acaricides. We evaluated and predicted the molecular properties of >180 acaricides using Molinspiration. We calculated physicochemical property parameters, such as log p, molecular weight (MW), and number of hydrogen bond donors (HBDs), hydrogen bond acceptors (HBAs), and rotatable bonds (Rot B). We then conducted qualitative and quantitative analyses of the physicochemical properties of acaricides. The MW of all acaricides ranged from 141 to 663, with an average value of 337.8. The number of HBDs ranged from 0 to 5, with an average value of 0.46. The number of HBAs ranged from 0 to 9, with an average value of 4.07. The log p ranged from -0.79 to 8.74, with an average value of 4.61. The number of Rot B ranged from 0 to 14, with an average value of 5.62. Except for the microbial and plant-derived acaricides, the molecular properties of the remaining acaricides are in accordance with the Lipinski's RO5. Therefore, the Lipinski's RO5 can provide a basis for screening new acaricide drugs.

Trait mapping in diverse arthropods by bulked segregant analysis

Bulked segregant analysis (BSA) is a cross-based method for genetic mapping in sexually reproducing organisms. The method's use of bulked (pooled) samples markedly reduces the genotyping effort associated with traditional linkage mapping studies. Further, it can be applied to species with life histories or physical attributes (as for micro-insects) that render genetic mapping with other methods impractical. Recent studies in both insects and mites have revealed that advanced BSA experimental designs can resolve causal loci to narrow genomic intervals, facilitating follow-up investigations. As high-quality genomes become more widely available, BSA methods are poised to become an increasingly important tool for the rapid mapping of both monogenic and polygenic traits in diverse arthropod species.