Transcriptome and Difference Analysis of Fenpropathrin Resistant Predatory Mite, Neoseiulus barkeri (Hughes)

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.

Combined PacBio Iso-Seq and Illumina RNA-Seq Analysis of the Tuta absoluta (Meyrick) Transcriptome and Cytochrome P450 Genes

Tuta absoluta (Meyrick) is a devastating invasive pest worldwide. The abamectin and chlorantraniliprole complex have become an alternative option for chemical control because they can enhance insecticidal activity and delay increased drug resistance. Notably, pests are inevitably resistant to various types of insecticides, and compound insecticides are no exception. To identify potential genes involved in the detoxification of abamectin and chlorantraniliprole complex in T. absoluta, PacBio SMRT-seq transcriptome sequencing and Illumina RNA-seq analysis of abamectin and chlorantraniliprole complex-treated T. absoluta were performed. We obtained 80,492 non-redundant transcripts, 62,762 (77.97%) transcripts that were successfully annotated, and 15,524 differentially expressed transcripts (DETs). GO annotation results showed that most of these DETs were involved in the biological processes of life-sustaining activities, such as cellular, metabolic, and single-organism processes. The KEGG pathway enrichment results showed that the pathways related to glutathione metabolism, fatty acid and amino acid synthesis, and metabolism were related to the response to abamectin and chlorantraniliprole complex in T. absoluta. Among these, 21 P450s were differentially expressed (11 upregulated and 10 downregulated). The qRT-PCR results for the eight upregulated P450 genes after abamectin and chlorantraniliprole complex treatment were consistent with the RNA-Seq data. Our findings provide new full-length transcriptional data and information for further studies on detoxification-related genes in T. absoluta.

Biochemical and Molecular Analysis of Field Resistance to Spirodiclofen in Panonychus citri (McGregor)

Spirodiclofen is one of the most widely used acaricides in China. The citrus red mite, Panonychus citri (McGregor) (Acari: Tetranychidae), is one of the most destructive citrus pests worldwide and has developed a high resistance to spirodiclofen. However, the molecular mechanism of spirodiclofen resistance in P. citri is still unknown. In this study, we identified a field spirodiclofen-resistant strain (DL-SC) that showed 712-fold resistance to spirodiclofen by egg bioassay compared to the susceptible strain. Target-site resistance was not detected as non-synonymous mutations were not found by amplification and sequencing of the ACCase gene of resistant and susceptible strains; in addition, the mRNA expression levels of ACCase were similar in both resistant and susceptible strains. The activity of detoxifying enzymes P450s and CCEs in the resistant strain was significantly higher than in the susceptible strain. The transcriptome expression data showed 19 xenobiotic metabolisms genes that were upregulated. Stage-specific expression profiling revealed that the most prominent upregulated gene, CYP385C10, in transcriptome data was significantly higher in resistant strains in all stages. Furthermore, functional analysis by RNAi indicated that the mortality caused by spirodiclofen was significantly increased by silencing the P450 gene CYP385C10. The current results suggest that overexpression of the P450 gene, CYP385C10, may be involved in spirodiclofen resistance in P. citri.

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

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

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.

Detrimental Impact of λ-Cyhalothrin on the Biocontrol Efficacy of Eocanthecona furcellata by Affecting Global Transcriptome and Predatory Behavior

Whether and how insecticide exposure will affect the biological control efficacy of predatory arthropods is critical in insecticide toxicology research but largely unexplored. In the current study, reduced biocontrol efficacy was observed in a predatory stink bug─Eocanthecona furcellata─after insecticide application in the field. Thus, we constructed a comparative transcriptome analysis and identified a total of 4364 upregulated and 1043 down regulated differentially expressed genes following the sublethal exposure of λ-cyhalothrin. The reduced juvenile hormone (JH) titer and increased trehalose content were observed. The predation capacity and theoretical maximum predation of predators were decreased by 31.08 and 48.90% in response to λ-cyhalothrin, respectively. Furthermore, JH supplementation after λ-cyhalothrin treatment could significantly stimulate trehalase and detoxification enzyme activities, as well as restore the predatory ability of E. furcellata. Our results help to understand the toxicological mechanism of predatory stink bug species in responding to insecticides, benefit predators' ecological services, and optimize the insecticide selection.

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

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

Next-generation biological control: the need for integrating genetics and genomics

Biological control is widely successful at controlling pests, but effective biocontrol agents are now more difficult to import from countries of origin due to more restrictive international trade laws (the Nagoya Protocol). Coupled with increasing demand, the efficacy of existing and new biocontrol agents needs to be improved with genetic and genomic approaches. Although they have been underutilised in the past, application of genetic and genomic techniques is becoming more feasible from both technological and economic perspectives. We review current methods and provide a framework for using them. First, it is necessary to identify which biocontrol trait to select and in what direction. Next, the genes or markers linked to these traits need be determined, including how to implement this information into a selective breeding program. Choosing a trait can be assisted by modelling to account for the proper agro‐ecological context, and by knowing which traits have sufficiently high heritability values. We provide guidelines for designing genomic strategies in biocontrol programs, which depend on the organism, budget, and desired objective. Genomic approaches start with genome sequencing and assembly. We provide a guide for deciding the most successful sequencing strategy for biocontrol agents. Gene discovery involves quantitative trait loci analyses, transcriptomic and proteomic studies, and gene editing. Improving biocontrol practices includes marker‐assisted selection, genomic selection and microbiome manipulation of biocontrol agents, and monitoring for genetic variation during rearing and post‐release. We conclude by identifying the most promising applications of genetic and genomic methods to improve biological control efficacy.

Improving the compatibility of pesticides and predatory mites: recent findings on physiological and ecological selectivity

Integrated pest management relies upon the application of selective pesticides that do not hinder biological control. Phytoseiid mites represent an interesting case-study: they are amongst the most frequently used biological control agents and often are less affected by pesticides than their prey by natural tolerance or by developing resistance. The selectivity of a pesticide is determined by physiological processes that include metabolism, transport, and the affinity to the target-site. Genomic and transcriptomic studies start to elucidate the genetic and molecular mechanisms of differential toxicity in some phytoseiid species, such as a mutation in the sodium channel conferring pyrethroid resistance. Ecological selectivity is achieved by smart applications of pesticides and management practices that influence the persistence of phytoseiid mites on plants. Although modern pesticides often show lower acute toxicity, there is a need for robust assays and procedures that quantify lethal and sublethal effects, through different routes and times of exposure.

Laboratory selection, resistance risk assessment, multi-resistance, and management of Tetranychus urticae Koch to bifenthrin, bifenazate and cyflumetofen on cowpea

BACKGROUND

Tetranychus urticae (T. urticae) Koch is an important pest of vegetable crops worldwide. In this study, bioassays were carried out to analyze the resistance risk, multi-resistance and management of T. urticae Koch to bifenthrin, bifenazate and cyflumetofen on cowpea.

RESULTS

The resistance ratios of the adult T. urticae population to bifenthrin (G16), bifenazate (G12) and cyflumetofen (G12) were 31.29, 9.38 and 5.81, respectively. Realized heritability (h ² ) analysis showed that, under a selection pressure of 50-90% mortality, the generations needed to increase 10-fold LC₅₀ values of bifenthrin, bifenazate and cyflumetofen were 3.64-8.05, 5.75-12.71, and 10.93-24.15, respectively. No obvious multi-resistance among these three acaricides was observed. Synergist bioassay results showed that microsomal multifunctional oxidase (MFO) was involved in bifenthrin resistance of T. urticae, with a synergistic ratio of 22.38. However, MFO and GSTs were not the main factors conferring the resistance to bifenazate. MFO, glutathione S-transferases(GSTs), together with esterase contributed to the development of the resistance to cyflumetofen. Additionally, the toxicity selection index test showed that bifenazate was safe to the natural enemy Neoseiulus barkeri (N. barkeri) with a toxicity selection index (TSI) >484.85, while bifenthrin was the least safe (TSI = 0.92).

CONCLUSIONS

These results demonstrated the T. urticae developed higher resistance risk to bifenthrin compared to bifenazate and cyflumetofen and no obvious multi-resistance among these three acaricides, providing guidance for designing appropriate strategies for the effective application of bifenthrin, bifenazate and cyflumetofen in the field and delaying the development of insecticide resistance. © 2019 Society of Chemical Industry.

De novo RNA-seq and functional annotation of Ornithonyssus bacoti

Ornithonyssus bacoti (Hirst) (Acari: Macronyssidae) is a vector and reservoir of pathogens causing serious infectious diseases, such as epidemic hemorrhagic fever, endemic typhus, tularemia, and leptospirosis. Its genome and transcriptome data are lacking in public databases. In this study, total RNA was extracted from live O. bacoti to conduct RNA-seq, functional annotation, coding domain sequence (CDS) prediction and simple sequence repeats (SSRs) detection. The results showed that 65.8 million clean reads were generated and assembled into 72,185 unigenes, of which 49.4% were annotated by seven functional databases. 23,121 unigenes were annotated and assigned to 457 species by non-redundant protein sequence database. The BLAST top-two hit species were Metaseiulus occidentalis and Ixodes scapularis. The procedure detected 12,426 SSRs, of which tri- and di-nucleotides were the most abundant types and the representative motifs were AAT/ATT and AC/GT. 26,936 CDS were predicted with a mean length of 711 bp. 87 unigenes of 30 functional genes, which are usually involved in stress responses, drug resistance, movement, metabolism and allergy, were further identified by bioinformatics methods. The unigenes putatively encoding cytochrome P450 proteins were further analyzed phylogenetically. In conclusion, this study completed the RNA-seq and functional annotation of O. bacoti successfully, which provides reliable molecular data for its future studies of gene function and molecular markers.

Gene cloning and difference analysis of vitellogenin in Neoseiulus barkeri (Hughes)

Neoseiulus barkeri (HUGHES) is the natural enemy of spider mites, whiteflies and thrips. Screening for chemically-resistant predatory mites is a practical way to balance the contradiction between the pesticide using and biological control. In this study, the number of eggs laid by fenpropathrin-susceptible and resistant strains of N. barkeri was compared. Additionally, we cloned three N. barkeri vitellogenin (Vg) genes and used quantitative real-time polymerase chain reaction to quantify Vg expression in susceptible and resistant strains. The total number of eggs significantly increased in the fenpropathrin-resistant strain. The full-length cDNA cloning of three N. barkeri Vg genes (NbVg1, NbVg2 and NbVg3) revealed that the open reading frames of NbVg1, NbVg2 and NbVg3 were 5571, 5532 and 4728 bp, encoding 1856, 1843 and 1575 amino acids, respectively. The three N. barkeri Vg possessed the Vitellogenin-N domain (or lipoprotein N-terminal domain (LPD_N)), von Willebrand factor type D domain (VWD) and the domain with unknown function 1943 (DUF1943). The NbVg1 and NbVg2 expression levels were significantly higher in the resistant strain than in the susceptible strain, while the NbVg3 expression level was lower in the resistant strain. Thus, we speculate that the increased number of eggs laid by the fenpropathrin-resistant strain of N. barkeri may be a consequence of changes in Vg gene expression.

Antioxidant Role of PcGSTd1 in Fenpropathrin Resistant Population of the Citrus Red Mite, Panonychus citri (McGregor)

The citrus red mite, Panonychus citri, a major citrus pest distributed worldwide, has evolved severe resistance to various classes of chemical acaricides/insecticides including pyrethroids. It is well known that the resistance to pyrethroids is mainly caused by point mutations of voltage-gated sodium channel gene in a wide range of pests. However, increasing number of evidences support that pyrethroids resistance might also be resulted from the integrated mechanisms including metabolic mechanisms. In this study, firstly, comparative analysis of RNA-seq data showed that multiple detoxification genes, including a GSTs gene PcGSTd1, were up-regulated in a fenpropathrin-resistant population compared with the susceptible strain (SS). Quantitative real time-PCR results showed that the exposure of fenpropathrin had an induction effect on the transcription of PcGSTd1 in a time-dependent manner. In vitro inhibition and metabolic assay of recombinant PcGSTd1 found that fenpropathrin might not be metabolized directly by this protein. However, its antioxidant role in alleviating the oxidative stress caused by fenpropathrin was demonstrated via the reversely genetic experiment. Our results provide a list of candidate genes which may contribute to a multiple metabolic mechanisms implicated in the evolution of fenpropathrin resistance in the field population of P. citri. Furthermore, during the detoxification process, PcGSTd1 plays an antioxidant role by detoxifying lipid peroxidation products induced by fenpropathrin.