Substrate specificity and promiscuity of horizontally transferred UDP-glycosyltransferases in the generalist herbivore Tetranychus urticae

UDP-Glycosyltransferases UGT350C3 and UGT344L7 Confer Tolerance to Neonicotinoids in Field Populations of Aphis gossypii

The cotton aphid, Aphis gossypii, is a polyphagous pest that stunts host plant growth via direct feeding or transmitting plant virus. Due to the long-term application of insecticides, A. gossypii has developed different levels of resistance to numerous insecticides. We found that five field populations had evolved multiple resistances to neonicotinoids. To explore the resistance mechanism mediated by uridine diphosphate glycosyltransferases (UGTs), two upregulated UGT genes in these five strains, UGT350C3 and UGT344L7, were selected for functional analysis of their roles in neonicotinoid detoxification. Transgenic Drosophila bioassay results indicated that compared with the control lines, the UGT350C3 and UGT344L7 overexpression lines were more tolerant to thiamethoxam, imidacloprid, and dinotefuran. Knockdown of UGT350C3 and UGT344L7 significantly increased A. gossypii sensitivity to thiamethoxam, imidacloprid, and dinotefuran. Molecular docking analysis demonstrated that these neonicotinoids could bind to the active pockets of UGT350C3 and UGT344L7. This study provides functional evidence of neonicotinoid detoxification mediated by UGTs and will facilitate further work to identify strategies for preventing the development of neonicotinoid resistance in insects.

Duplications and Losses of the Detoxification Enzyme Glycosyltransferase 1 Are Related to Insect Adaptations to Plant Feeding

Insects have developed sophisticated detoxification systems to protect them from plant secondary metabolites while feeding on plants to obtain necessary nutrients. As an important enzyme in the system, glycosyltransferase 1 (GT1) conjugates toxic compounds to mitigate their harm to insects. However, the evolutionary link between GT1s and insect plant feeding remains elusive. In this study, we explored the evolution of GT1s across different insect orders and feeding niches using publicly available insect genomes. GT1 is widely present in insect species; however, its gene number differs among insect orders. Notably, plant-sap-feeding species have the highest GT1 gene numbers, whereas blood-feeding species display the lowest. GT1s appear to be associated with insect adaptations to different plant substrates in different orders, while the shift to non-plant feeding is related to several losses of GT1s. Most large gene numbers are likely the consequence of tandem duplications showing variations in collinearity among insect orders. These results reveal the potential relationships between the evolution of GT1s and insect adaptation to plant feeding, facilitating our understanding of the molecular mechanisms underlying insect-plant interactions.

UDP-glycosyltransferases act as key determinants of host plant range in generalist and specialist Spodoptera species

Phytophagous insects have evolved sophisticated detoxification systems to overcome the antiherbivore chemical defenses produced by many plants. However, how these biotransformation systems differ in generalist and specialist insect species and their role in determining insect host plant range remains an open question. Here, we show that UDP-glucosyltransferases (UGTs) play a key role in determining the host range of insect species within the Spodoptera genus. Comparative genomic analyses of Spodoptera species that differ in host plant breadth identified a relatively conserved number of UGT genes in generalist species but high levels of UGT gene pseudogenization in the specialist Spodoptera picta. CRISPR-Cas9 knockouts of the three main UGT gene clusters of Spodoptera frugiperda revealed that UGT33 genes play an important role in allowing this species to utilize the poaceous plants maize, wheat, and rice, while UGT40 genes facilitate utilization of cotton. Further functional analyses in vivo and in vitro identified the UGT SfUGT33F32 as the key mechanism that allows generalist S. frugiperda to detoxify the benzoxazinoid DIMBOA (2,4-dihydroxy-7-methoxy-2H-1,4-benzoxazin-3(4H)-one), a potent insecticidal phytotoxin produced by poaceous plants. However, while this detoxification capacity is conserved in several generalist Spodoptera species, Spodoptera picta, which specializes on Crinum plants, is unable to detoxify DIMBOA due to a nonfunctionalizing mutation in SpUGT33F34. Collectively, these findings provide insight into the role of insect UGTs in host plant adaptation, the mechanistic basis of evolutionary transitions between generalism and specialism and offer molecular targets for controlling a group of notorious insect pests.

Uridine Diphosphate-Glycosyltransferase RpUGT344D38 Contributes to λ-Cyhalothrin Resistance in Rhopalosiphum padi

Uridine diphosphate-glycosyltransferase (UGT) is a key phase II enzyme in the insect detoxification system. Pyrethroids are commonly used to control the destructive wheat aphid Rhopalosiphum padi. In this study, we found a highly expressed UGT gene, RpUGT344D38, in both λ-cyhalothrin (LCR)- and bifenthrin (BTR)-resistant strains of R. padi. After exposure to λ-cyhalothrin and bifenthrin, the expression levels of RpUGT344D38 were significantly increased in the resistant strains. Knockdown of RpUGT344D38 did not affect the resistance of BTR, but it did significantly increase the susceptibility of LCR aphids to λ-cyhalothrin. Molecular docking analysis demonstrated that RpUGT344D38 had a stable binding interaction with both bifenthrin and λ-cyhalothrin. The recombinant RpUGT344D38 was able to metabolize 50% of λ-cyhalothrin. This study provides a comprehensive analysis of the role of RpUGT344D38 in the resistance of R. padi to bifenthrin and λ-cyhalothrin, contributing to a better understanding of aphid resistance to pyrethroids.

Structural and functional studies reveal the molecular basis of substrate promiscuity of a glycosyltransferase originating from a major agricultural pest

The two-spotted spider mite, Tetranychus urticae, is a major cosmopolitan pest that feeds on more than 1100 plant species. Its genome contains an unprecedentedly large number of genes involved in detoxifying and transporting xenobiotics, including 80 genes that code for UDP glycosyltransferases (UGTs). These enzymes were acquired via horizontal gene transfer from bacteria after loss in the Chelicerata lineage. UGTs are well-known for their role in phase II metabolism; however, their contribution to host adaptation and acaricide resistance in arthropods, such as T. urticae, is not yet resolved. TuUGT202A2 (Tetur22g00270) has been linked to the ability of this pest to adapt to tomato plants. Moreover, it was shown that this enzyme can glycosylate a wide range of flavonoids. To understand this relationship at the molecular level, structural, functional, and computational studies were performed. Structural studies provided specific snapshots of the enzyme in different catalytically relevant stages. The crystal structure of TuUGT202A2 in complex with UDP-glucose was obtained and site-directed mutagenesis paired with molecular dynamic simulations revealed a novel lid-like mechanism involved in the binding of the activated sugar donor. Two additional TuUGT202A2 crystal complexes, UDP-(S)-naringenin and UDP-naringin, demonstrated that this enzyme has a highly plastic and open-ended acceptor-binding site. Overall, this work reveals the molecular basis of substrate promiscuity of TuUGT202A2 and provides novel insights into the structural mechanism of UGTs catalysis.

The genome sequence of a spider mite, Tetranychus truncatus, provides insights into interspecific host range variation and the genetic basis of adaptation to a low-quality host plant

The phytophagous mite Tetranychus truncatus is a serious pest in East Asia but has a relatively narrower host range than the pest mite Tetranychus urticae, which can feed on over 1200 plant species. Here, we generated a high-quality chromosomal level genome of T. truncatus and compared it with that of T. urticae, with an emphasis on the genes related to detoxification and chemoreception, to explore the genomic basis underlying the evolution of host range. We also conducted population genetics analyses (in 86 females from 10 populations) and host transfer experiments (in 4 populations) to investigate transcription changes following transfer to a low-quality host (Solanum melongena, eggplant), and we established possible connections between fitness on eggplant and genes related to detoxification and chemoreception. We found that T. truncatus has fewer genes related to detoxification, transport, and chemoreception than T. urticae, with a particularly strong reduction in gustatory receptor (GR) genes. We also found widespread transcriptional variation among T. truncatus populations, which varied in fitness on eggplant. We characterized selection on detoxification-related genes through ω values and found a negative correlation between expression levels and ω values. Based on the transcription results, as well as the fitness and genetic differences among populations, we identified genes potentially involved in adaptation to eggplant in T. truncatus. Our work provides a genomic resource for this pest mite and new insights into mechanisms underlying the adaptation of herbivorous mites to host plants.

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.

Assessment of the role of an ABCC transporter TuMRP1 in the toxicity of abamectin to Tetranychus urticae

The rapid evolution of pest resistance threatens the sustainable utilization of bioinsecticides such as abamectin, and so deciphering the molecular mechanisms affecting toxicity and resistance is essential for their long-term application. Historical studies of abamectin resistance in arthropods have mainly focused on mechanisms involving the glutamate-gated chloride channel (GluCl) targets, with the role of metabolic processes less clear. The two-spotted spider mite, Tetranychus urticae, is a generalist herbivore notorious for rapidly developing resistance to pesticides worldwide, and abamectin has been widely used for its control in the field. After reanalyzing previous transcriptome and RNA-seq data, we here identified an ABC transporter subfamily C gene in T. urticae named multidrug resistance-associated protein 1 (TuMRP1), whose expression differed between susceptible and resistant populations. Synergism bioassays with the inhibitor MK-571, the existence of a genetic association between TuMRP1 expression and susceptibility to abamectin, and the effect of RNA interference mediated silencing of TuMRP1 were all consistent with a direct role of this transporter protein in the toxicity of abamectin. Although ABC transporters are often involved in removing insecticidal compounds from cells, our data suggest either an alternative role for these proteins in the mechanism of action of abamectin or highlight an indirect association between their expression and abamectin toxicity.

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.

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.

Contribution of UDP-glycosyltransferases to chlorpyrifos resistance in Nilaparvata lugens

As a multigene superfamily of Phase II detoxification enzymes, uridine diphosphate (UDP)-glycosyltransferases (UGTs) play important roles in the metabolism of xenobiotics including insecticides. In this study, 5-nitrouracil, an inhibitor of UGT enzyme activity, effectively increased the toxicity of chlorpyrifos to the chlorpyrifos-resistant strain of Nilaparvata lugens, one of the most resistant rice pests. The enzyme content of UGT in the resistant strain was significantly higher than that in the susceptible strain. Among 20 identified UGT genes, UGT386H2, UGT386J2, UGT386N2 and UGT386P1 were found significantly overexpressed in the resistant strain and can be effectively induced by chlorpyrifos. These four UGT genes were most highly expressed in the midgut and/or fat body, two main insect detoxification tissues. Amino acid sequence alignments revealed that these four UGTs contained a variable N-terminal substrate-binding domain and a conserved C-terminal sugar donor-binding domain. Furthermore, homology modeling and molecular docking analyses showed that these UGTs could stably bind to chlorpyrifos and chlorpyrifos oxon, with the binding free energies from -19.4 to -110.62 kcal mol^-1. Knockdown of UGT386H2 or UGT386P1 by RNA interference dramatically increased the susceptibility of the resistant strain to chlorpyrifos. These findings suggest that overexpression of these two UGT genes contributes to chlorpyrifos resistance in N. lugens.

Transcriptome analysis revealed detoxification gene expression changes in Tetranychus cinnabarinus challenged with ethyl oleate

Natural acaricides are potential biorational mite control alternatives to conventional chemical acaricides. However, little is known about the molecular mechanism of defense response to natural acaricides in mites. We previously reported significant acaricidal properties of ethyl oleate (EO) against Tetranychus cinnabarinus (here referred to as a sibling species of two-spotted spider mite, Tetranychus urticae), a highly polyphagous pest devastating crops in fields and greenhouses worldwide. In this study, we explored the molecular responses of T. cinnabarinus exposed to EO using RNA-Seq and differentially expressed gene (DEG) analysis. A total of 131, 185, and 154 DEGs were identified in T. cinnabarinus after 1, 6, and 24 h of EO treatment. In addition, 36 putative detoxification-related DEGs, including 10 cytochrome P450s (P450s), three glutathione S-transferases (GSTs), nine UDP-glycosyltransferases (UGTs), eight esterases (ESTs), and six ATP-binding cassette transporters (ABC transporters), were identified. Interestingly, the upregulation of these detoxification-related genes might be the main defense response of T. cinnabarinus exposed to EO. A quantitative real-time PCR analysis indicated that the expression profiles of 19 random DEGs were consistent with the RNA-Seq results. These findings serve as valuable information for a better understanding of the acaricide-mite interaction and molecular mechanisms involved in the defense response of T. cinnabarinus against EO.

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.

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.

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

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

The Evolution of Glycoside Hydrolase Family 1 in Insects Related to Their Adaptation to Plant Utilization

Insects closely interact with plants with multiple genes involved in their interactions. β-glucosidase, constituted mainly by glycoside hydrolase family 1 (GH1), is a crucial enzyme in insects to digest plant cell walls and defend against natural enemies with sequestered plant metabolites. To gain more insights into the role of this enzyme in plant-insect interactions, we analyzed the evolutionary history of the GH1 gene family with publicly available insect genomes. We found that GH1 is widely present in insects, while the gene numbers are significantly higher in insect herbivores directly feeding on plant cell walls than in other insects. After reconciling the insect GH1 gene tree with a species tree, we found that the patterns of duplication and loss of GH1 genes differ among insect orders, which may be associated with the evolution of their ecology. Furthermore, the majority of insects' GH1 genes were tandem-duplicated and subsequently went through neofunctionalization. This study shows the evolutionary history of an important gene family GH1 in insects and facilitates our understanding of the evolution of insect-plant interactions.

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

Background

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

Results

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

Conclusion

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

Supplementary Information

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

Comparative Genomics Reveals Insights into the Divergent Evolution of Astigmatic Mites and Household Pest Adaptations

Highly diversified astigmatic mites comprise many medically important human household pests such as house dust mites causing ∼1–2% of all allergic diseases globally; however, their evolutionary origin and diverse lifestyles including reversible parasitism have not been illustrated at the genomic level, which hampers allergy prevention and our exploration of these household pests. Using six high-quality assembled and annotated genomes, this study not only refuted the monophyly of mites and ticks, but also thoroughly explored the divergence of Acariformes and the diversification of astigmatic mites. In monophyletic Acariformes, Prostigmata known as notorious plant pests first evolved, and then rapidly evolving Astigmata diverged from soil oribatid mites. Within astigmatic mites, a wide range of gene families rapidly expanded via tandem gene duplications, including ionotropic glutamate receptors, triacylglycerol lipases, serine proteases and UDP glucuronosyltransferases. Gene diversification after tandem duplications provides many genetic resources for adaptation to sensing environmental signals, digestion, and detoxification in rapidly changing household environments. Many gene decay events only occurred in the skin-burrowing parasitic mite Sarcoptes scabiei. Throughout the evolution of Acariformes, massive horizontal gene transfer events occurred in gene families such as UDP glucuronosyltransferases and several important fungal cell wall lytic enzymes, which enable detoxification and digestive functions and provide perfect drug targets for pest control. This comparative study sheds light on the divergent evolution and quick adaptation to human household environments of astigmatic mites and provides insights into the genetic adaptations and even control of human household pests.

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

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

Knockdown of Uridine Diphosphate Glucosyltransferase 86Dg Enhances Susceptibility of Tribolium castaneum (Coleoptera: Tenebrionidae) to Artemisia vulgaris (Asterales: Asteraceae) Essential Oil

Uridine diphosphate glucosyltransferases (UGTs), which are phase II detoxification enzymes, are found in various organisms. These enzymes play an important role in the detoxification mechanisms of plant allelopathy and in insects. Artemisia vulgaris L. (Asterales: Asteraceae: Artemisia) essential oil has strong contact toxicity to Tribolium castaneum Herbst (Coleoptera: Tenebrionidae) larvae. However, the effect of A. vulgaris essential oil on UGTs is unclear. In this study, A. vulgaris essential oil was shown to significantly induce the expression of the TcUgt86Dg transcript. Furthermore, treatment of TcUgt86Dg-silenced individuals with A. vulgaris essential oil resulted in higher mortality than for the control individuals, indicating that TcUgt86Dg is involved in detoxification of A. vulgaris essential oil in T. castaneum. The developmental expression profile showed that the expression of TcUgt86Dg in late adults was higher than in other developmental stages. Furthermore, the expression profile in adult tissues revealed higher expression of TcUgt86Dg in the head, antenna, fat body, and accessory gland than in other tissues. These data show that TcUgt86Dg may be involved in the metabolism of exogenous toxins by T. castaneum; thus, our results have elucidated one possible mechanism of resistance to A. vulgaris essential oil and provide a theoretical basis for a control scheme for T. castaneum.

Rapid specialization of counter defenses enables two-spotted spider mite to adapt to novel plant hosts

Genetic adaptation, occurring over a long evolutionary time, enables host-specialized herbivores to develop novel resistance traits and to efficiently counteract the defenses of a narrow range of host plants. In contrast, physiological acclimation, leading to the suppression and/or detoxification of host defenses, is hypothesized to enable broad generalists to shift between plant hosts. However, the host adaptation mechanisms used by generalists composed of host-adapted populations are not known. Two-spotted spider mite (TSSM; Tetranychus urticae) is an extreme generalist herbivore whose individual populations perform well only on a subset of potential hosts. We combined experimental evolution, Arabidopsis thaliana genetics, mite reverse genetics, and pharmacological approaches to examine mite host adaptation upon the shift of a bean (Phaseolus vulgaris)-adapted population to Arabidopsis. We showed that cytochrome P450 monooxygenases are required for mite adaptation to Arabidopsis. We identified activities of two tiers of P450s: general xenobiotic-responsive P450s that have a limited contribution to mite adaptation to Arabidopsis and adaptation-associated P450s that efficiently counteract Arabidopsis defenses. In approximately 25 generations of mite selection on Arabidopsis plants, mites evolved highly efficient detoxification-based adaptation, characteristic of specialist herbivores. This demonstrates that specialization to plant resistance traits can occur within the ecological timescale, enabling the TSSM to shift to novel plant hosts.

Transcriptome analysis reveals important candidate gene families related to oligosaccharides biosynthesis in Morindaofficinalis

Morinda officinalis How (MO) is one of the best-known traditional herbs and is widely cultivated in subtropical and tropical areas for many years, especially in southern China. Oligosaccharides are the major constituents in the roots of MO, which is well known for its therapeutic effects with anti-depression, anti-osteoporosis, memory-enhancing, ect. To date, the main gene families that regulate the biosynthetic pathway of MO oligosaccharides metabolism yet have been published. In our study, six cDNA libraries generated from six plants of MO were sequenced utilizing an Illumina HiSeq 4000 platform. Corresponding totals of more than 132.60 million clean reads were obtained from the six libraries and assembled into 25,812 unigenes with an average length of 1288 bp. Moreover, 6036 unigenes were found to be allocated to 26 pathways maps using several public databases, and 2538 differential expression genes (DEGs) were screened. Among them, 25 genes from three families were selected as the mainly candidate genes related to MO oligosaccharides biosynthesis. Then, the expression patterns of six DEGs closely related to MO oligosaccharides biosynthesis were verified by quantitative real-time PCR (qRT-PCR). Besides, the MO was clustered more closely to Coffea arabica of Rubiaceae. In summary, the transcriptomic analysis was used to investigate the differences in expression genes of oligosaccharides biosynthesis, with the notable outcome that several key gene families were closely linked to oligosaccharides biosynthesis.

Delta class glutathione S-transferase (TuGSTd01) from the two-spotted spider mite Tetranychus urticae is inhibited by abamectin

GSTs (Glutathione S-transferases) are known to catalyze the nucleophilic attack of the sulfhydryl group of reduced glutathione (GSH) on electrophilic centers of xenobiotic compounds, including insecticides and acaricides. Genome analyses of the polyphagous spider mite herbivore Tetranychus urticae (two-spotted spider mite) revealed the presence of a set of 32 genes that code for secreted proteins belonging to the GST family of enzymes. To better understand the role of these proteins in T. urticae, we have functionally characterized TuGSTd01. Moreover, we have modeled the structure of the enzyme in apo form, as well as in the form with bound inhibitor. We demonstrated that this protein is a glutathione S-transferase that can conjugate glutathione to 1-chloro-2,4-dinitrobenzene (CDNB). We have tested TuGSTd01 activity with a range of potential substrates such as cinnamic acid, cumene hydroperoxide, and allyl isothiocyanate; however, the enzyme was unable to process these compounds. Using mutagenesis, we showed that putative active site variants S11A, E66A, S67A, and R68A mutants, which were residues predicted to interact directly with GSH, have no measurable activity, and these residues are required for the enzymatic activity of TuGSTd01. There are several reports that associate some T. urticae acaricide resistance with increased activity of GSTs . However, we found that TuGSTd01 is not able to detoxify abamectin; in fact, the acaricide inhibits the enzyme with K_i = 101 μM. Therefore, we suggest that the increased GST activity observed in abamectin resistant T. urticae field populations is a part of the compensatory feedback loop. In this case, the increased production of GSTs and relatively high concentration of GSH in cells allow GSTs to maintain physiological functions despite the presence of the acaricide.

Comparative analyses of transcriptional responses of Dectes texanus LeConte (Coleoptera: Cerambycidae) larvae fed on three different host plants and artificial diet

Dectes texanus is an important coleopteran pest of soybeans and cultivated sunflowers in the Midwestern United States that causes yield losses by girdling stems of their host plants. Although sunflower and giant ragweed are primary hosts of D. texanus, they began colonizing soybeans approximately 50 years ago and no reliable management method has been established to prevent or reduce losses by this pest. To identify genes putatively involved when feeding soybean, we compared gene expression of D. texanus third-instar larvae fed soybean to those fed sunflower, giant ragweed, or artificial diet. Dectes texanus larvae differentially expressed 514 unigenes when fed on soybean compared to those fed the other diet treatments. Enrichment analyses of gene ontology terms from up-regulated unigenes in soybean-fed larvae compared to those fed both primary hosts highlighted unigenes involved in oxidoreductase and polygalacturonase activities. Cytochrome P450s, carboxylesterases, major facilitator superfamily transporters, lipocalins, apolipoproteins, glycoside hydrolases 1 and 28, and lytic monooxygenases were among the most commonly up-regulated unigenes in soybean-fed larvae compared to those fed their primary hosts. These results suggest that D. texanus larvae differentially expressed unigenes involved in biotransformation of allelochemicals, digestion of plant cell walls and transport of small solutes and lipids when feeding in soybean.

Glycosyltransferases: the multifaceted enzymatic regulator in insects

Glycosyltransferases (GTs) catalyse the reaction of glyco-conjugation of various biomolecules by transferring the saccharide moieties from an activated nucleotide sugar to nucleophilic glycosyl acceptor. In insects, GTs show diverse temporal and site-specific expression patterns and thus play significant roles in forming the complex biomolecular structures that are necessary for insect survival, growth and development. Several insects exhibit GT-mediated detoxification as a key defence strategy against plant allelochemicals and xenobiotic compounds, as well as a mechanism for pesticide cross-resistance. Also, these enzymes act as crucial effectors and modulators in various developmental processes of insects such as eye development, UV shielding, cuticle formation, epithelial development and other specialized functions. Furthermore, many of the known insect GTs have been shown to play a fundamental role in other physiological processes like body pigmentation, cuticular tanning, chemosensation and stress response. This review provides a detailed overview of the multifaceted functionality of insect GTs and summarizes numerous case studies associated with it.

Whitefly hijacks a plant detoxification gene that neutralizes plant toxins

Plants protect themselves with a vast array of toxic secondary metabolites, yet most plants serve as food for insects. The evolutionary processes that allow herbivorous insects to resist plant defenses remain largely unknown. The whitefly Bemisia tabaci is a cosmopolitan, highly polyphagous agricultural pest that vectors several serious plant pathogenic viruses and is an excellent model to probe the molecular mechanisms involved in overcoming plant defenses. Here, we show that, through an exceptional horizontal gene transfer event, the whitefly has acquired the plant-derived phenolic glucoside malonyltransferase gene BtPMaT1. This gene enables whiteflies to neutralize phenolic glucosides. This was confirmed by genetically transforming tomato plants to produce small interfering RNAs that silence BtPMaT1, thus impairing the whiteflies' detoxification ability. These findings reveal an evolutionary scenario whereby herbivores harness the genetic toolkit of their host plants to develop resistance to plant defenses and how this can be exploited for crop protection.

The UDP-Glycosyltransferase Family in Drosophila melanogaster: Nomenclature Update, Gene Expression and Phylogenetic Analysis

UDP-glycosyltransferases (UGTs) are important conjugation enzymes found in all kingdoms of life, catalyzing a sugar conjugation with small lipophilic compounds and playing a crucial role in detoxification and homeostasis. The UGT gene family is defined by a signature motif in the C-terminal domain where the uridine diphosphate (UDP)-sugar donor binds. UGTs have been identified in a number of insect genomes over the last decade and much progress has been achieved in characterizing their expression patterns and molecular functions. Here, we present an update of the complete repertoire of UGT genes in Drosophila melanogaster and provide a brief overview of the latest research in this model insect. A total of 35 UGT genes are found in the D. melanogaster genome, localized to chromosomes 2 and 3 with a high degree of gene duplications on the chromosome arm 3R. All D. melanogaster UGT genes have now been named in FlyBase according to the unified UGT nomenclature guidelines. A phylogenetic analysis of UGT genes shows lineage-specific gene duplications. Analysis of anatomical and induced gene expression patterns demonstrate that some UGT genes are differentially expressed in various tissues or after environmental treatments. Extended searches of UGT orthologs from 18 additional Drosophila species reveal a diversity of UGT gene numbers and composition. The roles of Drosophila UGTs identified to date are briefly reviewed, and include xenobiotic metabolism, nicotine resistance, olfaction, cold tolerance, sclerotization, pigmentation, and immunity. Together, the updated genomic information and research overview provided herein will aid further research in this developing field.

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.

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.

Genome-Wide Identification and Expression Analysis of Udp-Glucuronosyltransferases in the Whitefly Bemisia Tabaci (Gennadius) (HemipterA: Aleyrodidae)

Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae) is an important agricultural pest worldwide. Uridine diphosphate (UDP)-glucuronosyltransferases (UGTs) are one of the largest and most ubiquitous groups of proteins. Because of their role in detoxification, insect UGTs are attracting increasing attention. In this study, we identified and analyzed UGT genes in B. tabaci MEAM1 to investigate their potential roles in host adaptation and reproductive capacity. Based on phylogenetic and structural analyses, we identified 76 UGT genes in the B. tabaci MEAM1 genome. RNA-seq and real-time quantitative PCR (RT-qPCR) revealed differential expression patterns of these genes at different developmental stages and in association with four host plants (cabbage, cucumber, cotton and tomato). RNA interference results of selected UGTs showed that, when UGT352A1, UGT352B1, and UGT354A1 were respectively silenced by feeding on dsRNA, the fecundity of B. tabaci MEAM1 was reduced, suggesting that the expressions of these three UGT genes in this species may be associated with host-related fecundity. Together, our results provide detailed UGTs data in B.tabaci and help guide future studies on the mechanisms of host adaptation by B.tabaci.

Genome streamlining in a minute herbivore that manipulates its host plant

The tomato russet mite, Aculops lycopersici, is among the smallest animals on earth. It is a worldwide pest on tomato and can potently suppress the host's natural resistance. We sequenced its genome, the first of an eriophyoid, and explored whether there are genomic features associated with the mite's minute size and lifestyle. At only 32.5 Mb, the genome is the smallest yet reported for any arthropod and, reminiscent of microbial eukaryotes, exceptionally streamlined. It has few transposable elements, tiny intergenic regions, and is remarkably intron-poor, as more than 80% of coding genes are intronless. Furthermore, in accordance with ecological specialization theory, this defense-suppressing herbivore has extremely reduced environmental response gene families such as those involved in chemoreception and detoxification. Other losses associate with this species' highly derived body plan. Our findings accelerate the understanding of evolutionary forces underpinning metazoan life at the limits of small physical and genome size.

Geographical distribution and molecular insights into abamectin and milbemectin cross-resistance in European field populations of Tetranychus urticae

BACKGROUND

Milbemectin and abamectin are frequently used to control the spider mite Tetranychus urticae. The development of abamectin resistance in this major pest has become an increasing problem worldwide, potentially compromising the use of milbemectin. In this study, a large collection of European field populations was screened for milbemectin and abamectin resistance, and both target-site and metabolic (cross-)resistance mechanisms were investigated.

RESULTS

High to very high levels of abamectin resistance were found in one third of all populations, while milbemectin resistance levels were low for most populations. The occurrence of well-known target-site resistance mutations in glutamate-gated chloride channels (G314D in GluCl1 and G326E in GluCl3) was documented in the most resistant populations. However, a new mutation, I321T in GluCl3, was also uncovered in three resistant populations, while a V327G and L329F mutation was found in GluCl3 of one resistant population. A differential gene-expression analysis revealed the overexpression of detoxification genes, more specifically cytochrome P450 monooxygenase (P450) and UDP-glycosyltransferase (UGT) genes. Multiple UGTs were functionally expressed, and their capability to glycosylate abamectin and milbemectin, was tested and confirmed.

CONCLUSIONS

We found a clear correlation between abamectin and milbemectin resistance in European T. urticae populations, but as milbemectin resistance levels were low, the observed cross-resistance is probably not of operational importance. The presence of target-site resistance mutations in GluCl genes was confirmed in most but not all resistant populations. Gene-expression analysis and functional characterization of P450s and UGTs suggests that also metabolic abamectin resistance mechanisms are common in European T. urticae populations. © 2020 Society of Chemical Industry.

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.