Front. Plant Sci.Frontiers in Plant ScienceFront. Plant Sci.1664-462XFrontiers Media S.A.10.3389/fpls.2021.671286Plant ScienceOriginal ResearchIdentification of a Sulfatase that Detoxifies Glucosinolates in the Phloem-Feeding Insect Bemisia tabaci and Prefers Indolic GlucosinolatesManivannanAbinaya1IsraniBhawana1†LuckKatrin1†GötzMonika2†SeibelElena1EassonMichael L. A. E.1KirschRoy1ReicheltMichael1SteinBeate2WinterStephan2GershenzonJonathan1VassãoDaniel Giddings1*1Max Planck Institute for Chemical Ecology, Jena, Germany2Leibniz Institute DSMZ-German Collection of Microorganisms and Cell Cultures, Braunschweig, Germany
Edited by: Jens Rohloff, Norwegian University of Science and Technology, Norway
Reviewed by: Inga Mewis, Humboldt University of Berlin, Germany; Shaoli Wang, Insititute of Vegetables and Flowers (CAAS), China
*Correspondence: Daniel Giddings Vassão, vassao@ice.mpg.de
†These authors have contributed equally to this work
This article was submitted to Plant Metabolism and Chemodiversity, a section of the journal Frontiers in Plant Science
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Cruciferous plants in the order Brassicales defend themselves from herbivory using glucosinolates: sulfur-containing pro-toxic metabolites that are activated by hydrolysis to form compounds, such as isothiocyanates, which are toxic to insects and other organisms. Some herbivores are known to circumvent glucosinolate activation with glucosinolate sulfatases (GSSs), enzymes that convert glucosinolates into inactive desulfoglucosinolates. This strategy is a major glucosinolate detoxification pathway in a phloem-feeding insect, the silverleaf whitefly Bemisia tabaci, a serious agricultural pest of cruciferous vegetables. In this study, we identified and characterized an enzyme responsible for glucosinolate desulfation in the globally distributed B. tabaci species MEAM1. In in vitro assays, this sulfatase showed a clear preference for indolic glucosinolates compared with aliphatic glucosinolates, consistent with the greater representation of desulfated indolic glucosinolates in honeydew. B. tabaci might use this detoxification strategy specifically against indolic glucosinolates since plants may preferentially deploy indolic glucosinolates against phloem-feeding insects. In vivo silencing of the expression of the B. tabaci GSS gene via RNA interference led to lower levels of desulfoglucosinolates in honeydew. Our findings expand the knowledge on the biochemistry of glucosinolate detoxification in phloem-feeding insects and suggest how detoxification pathways might facilitate plant colonization in a generalist herbivore.
Plants rely on a complex arsenal of toxic chemicals to defend themselves against herbivores and pathogens. One successful defense strategy used to safely accumulate large concentrations of defensive compounds while preventing auto-toxicity is the production of two-component activated defenses, such as cyanogenic glucosides and glucosinolates (Halkier and Gershenzon, 2006; Mithöfer and Boland, 2012). These glucosylated pro-toxins are stable and only become toxic after enzymatic activation coincident with herbivore damage. Multi-component plant defenses like this, however, offer herbivores multiple targets for counter-adaptation, including mechanisms to prevent or redirect activation, or detoxify the activated poisons (Pentzold et al., 2014; Jeschke et al., 2016b). Piercing–sucking, phloem-feeding insects, such as aphids and whiteflies, are thought to lessen plant defensive responses by causing only minimal damage during feeding, stealthily penetrating plant tissues to feed on the sugar-rich phloem sap. Nevertheless, even piercing–sucking feeding leads to some activation of glucosinolates (Kim et al., 2008; Danner et al., 2018), so that detoxification strategies preventing activation of two-component defenses can be greatly advantageous to phloem-feeding insects. Here, we examine the biochemical basis of one such detoxification in the silverleaf whitefly Bemisia tabaci (Gennadius) (Hemiptera): the desulfation of glucosinolates to form derivatives that can no longer be activated.
Glucosinolates are sulfur-containing defensive metabolites restricted to the plant order Brassicales. The family Brassicaceae of this order (mustard family) includes many economically important oilseed, vegetable, condiment, and fodder crops, such as rapeseed, cabbages, and mustards, as well as the model plant Arabidopsis thaliana (L.) Heynh. These plants have distinctive tastes and smells that are conferred in great part by hydrolysis products from glucosinolates, also called mustard oils (Halkier, 2016). The core chemical structure of glucosinolates (Figure 1) consists of a β-thioglucose residue, a sulfonated oxime moiety, and a variable amino acid-derived side chain group R (Halkier and Gershenzon, 2006). A recent review lists up to 137 natural glucosinolates described based on the variability at the R group (Blažević et al., 2020), with these being classifiable according to the precursor amino acid from which they are derived, such as aliphatic (e.g., from Met, Ala, and Val), benzenic (Phe and Tyr), and indolic (Trp).
Glucosinolate activation by myrosinase leads to nitriles and toxic isothiocyanates, whereas desulfation by glucosinolate sulfatase (GSS) forms desulfoglucosinolates that cannot be activated. Below are examples of side chains of aliphatic, benzenic, and indolic glucosinolates derived from different amino acid precursors.
In addition to glucosinolates, cruciferous plants produce thioglucosidase enzymes known as myrosinases that catalyze glucosinolate hydrolysis (Rask et al., 2000). Myrosinases and glucosinolates are physically segregated in the intact plant (Grob and Matile, 1979; Höglund et al., 1991). While the localization of each component is still a topic of study, myrosinases have been reported in protein-accumulating idioblasts called myrosin cells present in the phloem parenchyma. Glucosinolates, on the other hand, are enriched in sulfur-containing S-cells and would thus be separated from myrosinases at a cellular level (Koroleva et al., 2000). However, in other spatial separation models, glucosinolates are confined into the vacuoles, whereas myrosinases are present in the cytoplasm of the same cells, enabling compartmentalization at a subcellular level (Lüthy and Matile, 1984). Upon tissue disruption, such as during herbivory, this compartmentalization is compromised, and myrosinases hydrolyze the glucose moiety of the glucosinolate forming a mixture of toxic products (Figure 1). Isothiocyanates are considered the most toxic components of the mixture. This is attributed to their polarity and reactivity. The lipophilic nature of the isothiocyanate side chains allows them to easily cross cellular membranes, and the electrophilic nature of the –S=C=N group promotes reaction with nucleophilic targets (Kawakishi and Kaneko, 1987). The release of toxic isothiocyanates from glucosinolates is often referred to as the “mustard oil bomb.” By allowing the accumulation of large (millimolar) concentrations of glucosinolates that are only activated into more chemically reactive toxins upon herbivory, this two-component system serves as an important chemical defense against most herbivores (Hopkins et al., 2009; Jeschke et al., 2016a; Schweizer et al., 2017).
Several species of insect herbivores, however, have developed means to detoxify glucosinolates and successfully feed on crucifer plants with relative impunity (Jeschke et al., 2016b). For example, the isothiocyanate hydrolysis products can be converted to hydrophilic glutathione–isothiocyanate conjugates (Schramm et al., 2012; Jeschke et al., 2016a). This reaction, which takes place in many organisms including humans, can occur very efficiently in certain insects adapted to crucifers. On the other hand, the larvae of some crucifer-feeding Pieris species use a “nitrile specifying protein” for directing glucosinolate hydrolysis toward the formation of less toxic nitriles rather than isothiocyanates (Wittstock et al., 2004). Another enzymatic detoxification mechanism is the action of glucosinolate sulfatases (GSSs) in the insect gut. These enzymes are used by crucifer specialists, such as the diamondback moth larvae [Plutella xylostella (L.)] (Ratzka et al., 2002) and the cabbage stem flea beetle [Psylliodes chrysocephala (L.)] (Ahn et al., 2019), as well as the generalist herbivores desert locust [Schistocerca gregaria (Forskål)] (Falk and Gershenzon, 2007) and the silverleaf whitefly (B. tabaci) (Malka et al., 2016). GSSs cleave the sulfate group of glucosinolates to form desulfoglucosinolates (Figure 1), essentially reversing the last step in the biosynthesis of the glucosinolate core structure. As desulfoglucosinolates can no longer be used as substrates by myrosinases, no toxic downstream products are produced.
Among the insects that can desulfate glucosinolates is B. tabaci, a piercing-sucking, phloem-feeder that inhabits mainly tropical and subtropical regions. B. tabaci is a cryptic species complex, existing as a conglomerate of dozens of species (Perring, 2001) that are difficult to distinguish even with molecular markers (Boykin et al., 2012, 2018). The most invasive and geographically widespread species are MED (Mediterranean, formerly named Q) and MEAM1 (Middle East–Asia Minor 1, formerly B). B. tabaci species are altogether extreme generalists with a wide host plant range, colonizing more than 600 plant species including food and ornamental crops (Oliveira et al., 2001), and are, therefore, one of the most serious agricultural pests in the world. B. tabaci has been recently reported to desulfate glucosinolates with the resulting desulfoglucosinolate products excreted in honeydew (Malka et al., 2016). These transformation products were not observed in honeydew from other phloem feeders examined [the green peach aphid Myzus persicae (Sulzer) and the cabbage whitefly Aleyrodes proletella (L.) – unpublished data], suggesting that desulfation is not widespread in hemipteran insects. Indeed, among the insects currently known to detoxify glucosinolates by desulfation, B. tabaci is the only phloem-feeder. This specialized feeding mode minimizes mechanical damage during feeding, but still leads to detectable glucosinolate activation (Danner et al., 2018), especially of indolic glucosinolates (Kim et al., 2008). Pre-emptive deactivation of glucosinolates, such as via GSS activity, may therefore confer benefits to phloem-feeding insects as well.
In this study, we aimed to identify and characterize the enzyme(s) performing glucosinolate desulfation in B. tabaci. Several genes encoding potential GSSs were cloned and heterologously expressed, and the resulting proteins were tested for activity. However, only a single GSS was observed to efficiently desulfate glucosinolates in B. tabaci, and this enzyme had a strong preference for indolic glucosinolates. Silencing of BtGSS gene expression in vivo reduced the proportion of desulfoglucosinolates in the excreted honeydew relative to intact glucosinolates, confirming its role in glucosinolate metabolism in this insect.
Materials and MethodsPlants and Insects
Bemisia tabaci (MEAM1, formerly known as biotype B) was reared in the Leibniz Institute DSMZ (German Collection of Microorganisms and Cell Cultures), Braunschweig, Germany and in the Hebrew University of Jerusalem, Israel on eggplant [Solanum melongena (L.)] and Brussels sprouts [Brassica oleracea (L.)] plants. A. thaliana Col-0 was cultivated in a controlled-environment growth chamber under short day conditions (9.5:13.5 h, L:D, 140 μmol/m2/s photosynthetic photon flux density, 21°C, 50–60% relative humidity).
Chemicals
Solvents used were HPLC grade and were obtained from VWR (Darmstadt, Germany). Glucosinolates were obtained from Sigma-Aldrich (Munich, Germany) and Phytoplan (Heidelberg, Germany). Arylsulfatase from Helix pomatia was obtained from Sigma-Aldrich and further processed as in Graser et al. (2000) and was used for preparation of desulfoglucosinolate standards by incubation of pure glucosinolates.
<italic>Bemisia tabaci</italic> Sulfatase Candidate Selection and Their Tissue-Specific Expression Levels
The gene sequences of the putative B. tabaci MEAM1 GSSs were obtained from the Whitefly Genome Database (whiteflygenomics.org, downloaded January 26, 2017) and NCBI (downloaded July 11, 2018). B. tabaci arylsulfatase gene candidates were picked based on BLAST similarities to the sequence of the previously published P. xylostella GSS (PxGSS, Px0181041). The initial screening of B. tabaci sulfatase candidates included the estimation of ratios of gene expression between dissected B. tabaci gut tissues and whole insect bodies, using previously published data (Ye et al., 2014; Luan et al., 2015) SRR835757 and SRR1523522, previously also used for gene expression analysis in Jing et al. (2016). Expression levels in FPKM were estimated using DNAStar QSeq with the software parameters used in Vogel et al. (2014).
Phylogenetic Analysis of GSSs and Arylsulfatase Candidates
A multiple alignment of the nucleotide sequences of the putative B. tabaci sulfatases and P. xylostella GSS was performed using the MUSCLE algorithm (Edgar, 2004) in the Geneious Prime software (default parameters, eight maximum iterations). The phylogenetic relationships among sequences were inferred using the neighbor-joining method and the Tamura-Nei genetic distance model, with the reliability of the tree branching tested with 1,000 bootstrap replicates, using the Geneious software.
Primer Design, PCR Amplification, and Cloning of Candidate Genes
The gene sequences of the putative B. tabaci MEAM1 GSSs (obtained as described above) were used for the design of primers for the amplification of the respective open reading frames (ORFs). Primer design was carried out using the Geneious software (version 10.0.5; Biomatters Ltd., Auckland, New Zealand). Full-length ORFs were amplified using Phusion HF Polymerase (Thermo Scientific, Darmstadt, Germany) using B. tabaci cDNA as template for amplification. DNA fragments were then extracted from agarose gels using either QIAquick Gel Extraction Kit (Qiagen, Hilden, Germany) or the Zymoclean Gel DNA Recovery Kit (Zymo Research, Freiburg im Breisgau, Germany) as per the manufacturer’s instructions and cloned into the pIB/V5-His TOPO TA expression vector (Thermo Scientific). Final plasmids were transfected into Sf9 cells by lipid-mediated transfection using FuGENE Transfection Reagent (Promega, Mannheim, Germany) according to the manufacturer’s instructions for recombinant protein expression. A pIB/V5-His TOPO TA plasmid containing the corresponding full-length PxGSS from P. xylostella was also transfected and used as a positive expression and activity control. Spodoptera frugiperda Sf9 cells (Life Technologies, Darmstadt, Germany) were cultured in Sf-900 II serum-free medium (Life Technologies). Adherent cell cultures were maintained at 27°C and subcultured every 3–4 days. Collected cells and media were separated by gentle centrifugation at 500 × g for 5 min, with cell pellets then re-suspended and homogenized to generate crude cell-free protein extracts.
Screening <italic>Bemisia tabaci</italic> Arylsulfatase Candidates for GSS Activity
To identify B. tabaci GSSs, in vitro enzyme assays were performed with crude protein extracts and media of cells heterologously expressing B. tabaci arylsulfatases. Reactions were performed in Tris buffer (100 mM, pH 7.5), to which 2–5 mM of glucosinolates being tested was added. The reactions were incubated at 25°C for 30 min and stopped using acetic acid (10% v/v). Assays were performed in duplicates for screening and triplicates for product quantification. Enzymatic assays with recombinant P. xylostella GSS were used as a positive control. The formation of desulfated glucosinolate products was quantified by high-performance liquid chromatography mass spectrometry (HPLC–MS/MS) using external calibration curves.
Polyhistidine-Tagged Protein Purification
The polyhistidine tagged, heterologously produced B. tabaci GSS BtGSS was affinity-purified from the extracellular culture medium. First, the medium was collected and concentrated to 1 ml using centrifugal filter units (Amicon Ultra Centrifugal Filters; Merck, Darmstadt, Germany) with a molecular weight cut-off of 10,000 Da. Purification over Ni-NTA agarose resin (Qiagen, Hilden, Germany) followed the protocol of Eakteiman et al. (2018).
Determination of <italic>B. tabaci</italic> GSS Optimal Temperature and pH
For determination of optimal reaction temperature, BtGSS was incubated with allyl glucosinolate (2 mM) in Tris buffer (100 mM, pH 7.5) for 10 min at 10, 15, 25, 35, 45, 55, 65, and 75°C. The enzyme reactions were stopped using acetic acid (10% v/v). Enzymatic assays were performed in duplicate. The desulfated allyl glucosinolate product formed was quantified by HPLC–MS/MS.
For determination of optimum pH, BtGSS was incubated with allyl glucosinolate (2 mM) in each of the following buffers: broad range: phosphate citrate (0.2 M phosphate, 0.1 M citrate; pH 3.0, 4.0, 5.0, 6.0, 7.0), Tris–HCl (0.1 M; pH 7.0, 8.0), and glycine–NaOH buffer (0.2 M; pH 8.0, 9.0); narrow range: sodium phosphate (0.1 M; pH 6.25, 6.50, 6.75, 7.00, 7.25, 7.50, 7.75, 8.00). Reactions were carried out for 10 min at 25°C and stopped using acetic acid (10% v/v). Assays were performed in duplicate. The desulfated allyl glucosinolate product formed was quantified by HPLC–MS/MS.
Quantitative Real-Time PCR Analysis of Basal <italic>BtGSS</italic> Expression
To check for the expression levels of BtGSS, quantitative real-time PCR (qPCR) analyses were performed using cDNA generated from four biological replicates of B. tabaci reared on eggplant (non-glucosinolate-containing diet) and Brussels sprouts (glucosinolate-containing diet). Ribosomal protein L-13 (Bta04282) was used as a reference gene. Primers were designed with an optimal melting temperature of 60°C. The list of primers is available in Supplementary Table 1. Each reaction was done in triplicate. Expression levels of BtGSS were relatively quantified using the 2–ΔΔCT method (Pfaffl, 2001) and are presented as means ± standard errors.
Silencing of <italic>BtGSS</italic> Expression <italic>in vivo</italic> and Resulting Metabolic Changes
Silencing experiments were carried out using artificial diets containing sucrose, 4msob glucosinolate, and two double-stranded RNases (dsRNAs), one targeting dsRNase2 and the other BtGSS, or a scrambled BtGSS-derived sequence used as a negative control (synthesis and cloning into pUC57 backbone by GenScript, Leiden, Netherlands). Initial amplification was carried out from B. tabaci cDNA (or pUC57 plasmid for the scrambled sequence) as template, using gene-specific primers containing the T7 polymerase site (Supplementary Table 2). The PCR products were used as templates for in vitro transcription reactions using the MEGAscript kit (Thermo Fisher Scientific). The products obtained were quantified using Nanodrop (Thermo Fisher Scientific), and the integrity of the dsRNA was confirmed on a 1.2% denaturing agarose gel. Feeding experiments were carried out in glass vials, with 50–100 whiteflies per vial. Aqueous artificial diets were enclosed between two thin layers of stretched parafilm and contained 0.5 μg/μl BtGSS (or scrambled BtGSS) dsRNA, 0.5 μg/μl dsRNase2 dsRNA, 0.29 M sucrose, and 5 mM 4msob glucosinolate. Silencing was carried out over a period of 48 h. At the end of the feeding period, whiteflies and honeydew were collected for gene expression and metabolite analysis, respectively. Whiteflies were homogenized in 500 μl TRIzol reagent (Thermo Fisher Scientific) for total RNA isolation, and cDNA synthesis was carried out with 0.5 μg total RNA using SuperScript IV Reverse Transcriptase (Thermo Fisher Scientific) according to the manufacturer’s protocol. Honeydew samples were processed by washing the glass vials twice with 500 μl methanol, concentrating under nitrogen flow and re-suspending in 200 μl water. Honeydew samples were used for quantification of 4msob desulfoglucosinolate and non-metabolized 4msob glucosinolate by HPLC–MS/MS.
Glucosinolate Extraction From <italic>Arabidopsis thaliana</italic>
Leaves from A. thaliana Col-0 rosette stage plants were harvested, flash frozen using liquid N2, and freeze-dried (Alpha 1-4 LDplus freeze dryer; Martin Christ, Osterode am Harz, Germany). Intact glucosinolates were extracted with 10 ml 80% (v/v) methanol:water per g dry weight (DW) on ice under continuous shaking for 5 min, followed by the addition of metal beads (3 mm) and vigorous shaking for 10 min in a paint shaker. The supernatant collected after centrifugation (4,000 × g at 4°C for 15 min) was then processed through centrifugal filter units (Amicon Ultra Centrifugal Filters) with a molecular weight cut-off of 10,000 Da to separate glucosinolates from the plant myrosinase. The flow-through was then treated on a rotatory evaporator for removal of the solvent and was resuspended in water.
Tests of <italic>Bemisia tabaci</italic> Sulfatase Substrate Preference With <italic>Arabidopsis</italic> Glucosinolates
To test the substrate preference of the B. tabaci GSS, enzyme assays were performed with the glucosinolate extract obtained from A. thaliana Col-0 leaves. Purified BtGSS was added to Tris buffer (100 mM, pH 7.5) and incubated with glucosinolate extract at 25°C. Aliquots were taken at different time points (0–300 min) and mixed with acetic acid (10% v/v). Enzymatic assays were performed in duplicate, and no-enzyme reactions served as a control for non-enzymatic glucosinolate degradation. The intact glucosinolates remaining in each reaction sample were quantified by HPLC–MS/MS.
Quantitative Tests of BtGSS Substrate Preference With Selected Pure Glucosinolates
To test the substrate preference of the B. tabaci GSS, enzyme assays were performed with selected pure glucosinolates: 4msob-, i3m-, allyl-, and pOHBn-glucosinolates. Purified BtGSS enzyme (1 μg) was added to Tris buffer (0.1 M, pH 7.5, 100 μl total volume) and was incubated with 2 mM glucosinolate extract at 25°C. Reactions proceeded for 10 min and were then stopped with acetic acid (10% v/v). Enzymatic assays were performed in triplicate, and no-enzyme reactions served as a control for non-enzymatic glucosinolate degradation. The desulfoglucosinolates formed from 4msob and i3m glucosinolates were quantified relative to external calibration curves of the respective products by HPLC–MS/MS, as described below.
High-Performance Liquid Chromatography Mass Spectrometry
Analyses of desulfated and intact glucosinolates in honeydew and enzyme assays were performed on an Agilent Technologies (Santa Clara, CA, United States) 1200 Series HPLC using a Nucleodur Sphinx RP column (250 × 4.6 mm × 5 μm; Macherey-Nagel, Düren, Germany) coupled to an API 3200 triple-quadrupole mass spectrometer (Applied Biosystems, Darmstadt, Germany). Formic acid (0.2%) in water and acetonitrile were employed as mobile phases A and B, respectively. The flow rate was 1.1 ml/min. The elution profile was as follows: 0–2.5 min, 1.5% B; 2.5–5 min, 1.5–10% B; 5–12.5 min, 10–40% B; 12.5–17.5 min, 40–70% B; 17.6–20 min, 100% B; and 20.1–24 min, 1.5% B. The ion spray voltage was maintained at 4,500 eV. The turbo gas temperature was set to 700°C. Nebulizing gas was set at 70 psi, curtain gas at 20 psi, heating gas at 60 psi, and collision gas at 10 psi. Multiple reaction monitoring (MRM, Supplementary Tables 3, 4) was used to monitor parent ion to fragment ion conversion (Malka et al., 2016). Analyst 1.5 software (Applied Biosystems) was used for data acquisition and processing.
Analyses of desulfated and intact glucosinolates during RNAi silencing experiments were performed on an identical HPLC setup as those above, but coupled to an API 6500 triple-quadrupole mass spectrometer (Applied Biosystems). Formic acid (0.05%) in water and acetonitrile were employed as mobile phases A and B, respectively. The elution profile was identical to the abovementioned but with a flow rate of 1.0 ml/min. Ion spray voltages were maintained at 4,500 eV (positive mode for desulfoglucosinolates) and -4,500 eV (negative mode for intact glucosinolates). The turbo gas temperature was set to 650°C. Curtain gas was set at 40 psi, collision gas at medium, ion source gas 1 at 60 psi, and ion source gas 2 at 60 psi. MRM was used to monitor parent ion to fragment ion conversion (Supplementary Table 5). Analyst 1.6.3 software (Applied Biosystems) was used for data acquisition and processing.
ResultsIdentification of Putative <italic>Bemisia tabaci</italic> GSS
The protein sequence of the previously identified P. xylostella GSS (Ratzka et al., 2002) was used as a query to identify B. tabaci sulfatases (BtSulf) in the published B. tabaci MEAM1 draft genome (Chen et al., 2016) that might serve as GSSs. Candidates were sought that had high similarity to PxGSS and a high expression level of the encoding genes in B. tabaci guts relative to the whole insect based on publicly available expression data (Table 1 and Figure 2). Most candidates had been previously annotated as arylsulfatases. Selected sequences were cloned from B. tabaci (species MEAM1) whole-insect cDNA and heterologously expressed in Sf9 cells, and the produced protein was subsequently screened for GSS activity.
Predicted Bemisia tabaci sulfatase (“BtSulf”) candidates, ordered based on sequence similarity scores (e-values) to Plutella xylostella GSS (PxGSS).
MEAM
MED
NCBI entry
Bta entry
Name
e-value (vs. PxGSS)
MEAM gut (FPKM)
MEAM adults (FPKM)
Expression ratios, gut vs. whole adults
MED gut (FPKM)
MED adults (FPKM)
XM 019054398
Bta06756
BtSulf1
2,60E-122
0
216
0
0
138
XM 019062599
Bta04898
BtSulf2
4,21E-112
0
81
0
0
20
XM 019061487
Bta03550
BtSulf3
1,18E-107
0
70
0
3
488
XM 019062598
Bta04899
BtSulf4
2,09E-107
0
248
0
1
275
XM 019051255
Bta02222
BtSulf5
1,50E-106
460
675
1,5 lower
370
945
XM 019052592
Bta04774
BtSulf6
1,11E-101
16
277
17,3 lower
5
103
XM 019059032
Bta14665
BtSulf7
5,36E-100
48
272
5,7 lower
81
201
XM 019058912
Bta14669
BtSulf8
1,02E-97
0
20
0
0
93
XM 019059016
Bta14666
BtSulf9
4,87E-96
8526
1203
7,1 higher
7296
135
XM 019049969
Bta01141
BtSulf10
2,88E-83
10
39
3,9 lower
96
62
XM 019059017
Bta14667
BtSulf11
1,76E-65
10
0
higher
2
0
XM 019042567
Bta08750
BtSulf12
5,73E-21
2627
9434
3,6 lower
1587
9732
XM 019054362
Bta06730
BtSulf13
7,94E-17
18
468
26,0 lower
12
746
XM 019049823
Bta01054
BtSulf14
3,13E-13
40
207
5,2 lower
15
55
XM 019053003
Bta05280
BtSulf15
3,59E-08
51
300
5,9 lower
38
360
Genes are denoted according to their accession numbers in NCBI and the whitefly genome databases. Estimated expression levels (FPKM) for these genes in gut tissues and whole insect bodies were calculated from publicly available data for B. tabaci MEAM1 insects and, as a reference, also for B. tabaci MED insects.
Phylogenetic relationships of predicted Bemisia tabaci sulfatases (BtSulf), including BtSulf9 (BtGSS, italicized) and P. xylostella GSS (PxGSS). Gene names are as described in Table 1, and bootstrap values (%) are shown next to each node.
<italic>Bemisia tabaci BtSulf9</italic> Encodes a GSS Enzyme (BtGSS)
Previous reports showed that PxGSS produced in Escherichia coli had no sulfatase activity (Ratzka et al., 2002), suggesting that those cells lacked the post-translational modification machinery necessary to generate an active GSS. Hence, insect cells (S. frugiperda Sf9 cells) were used here for expressing the B. tabaci GSS candidates. The full-length ORFs of selected B. tabaci sulfatases (BtSulf1, 3, 5–9, and 11) were successfully expressed to produce V5-/His-tagged proteins. The resulting B. tabaci enzymes were screened for their activity toward pure glucosinolates: allyl-(sinigrin), p-hydroxybenzyl-(pOHBn, sinalbin), and 4-methylsulfinylbutyl-(4-msob, glucoraphanin) glucosinolates. Heterologously produced PxGSS and commercial H. pomatia sulfatase were used as positive controls for GSS activity. The resulting desulfoglucosinolate products were detected by HPLC–MS/MS. Intact glucosinolates also gave rise to MS signals corresponding to desulfated glucosinolates, due to in-source fragmentation (spontaneous sulfate loss) in the mass spectrometer. However, these artifactual desulfated glucosinolates were readily distinguished from genuine desulfoglucosinolates existing before analysis due to different retention times. These analyses revealed that BtSulf9 had GSS activity toward all tested glucosinolates (Figure 3 and Supplementary Figure 1) and this was henceforth named BtGSS. Both BtGSS and PxGSS were secreted by the Sf9 cells into the extracellular culture medium. BtGSS purified by metal affinity chromatography had an optimum temperature around 50–55°C and an optimum pH between 7.5 and 7.8. None of the other produced sulfatases tested displayed detectable GSS activity in vitro, including BtSulf1 (the closest homolog to PxGSS) and those in the same phylogenetic branch as BtGSS (BtSulf7, 8, 11).
BtGSS has activity with a variety of glucosinolates. (A–D) LC–MS extracted MRM chromatograms of desulfated glucosinolates produced in vitro by BtGSS. The early eluting peaks in each chromatogram indicate genuine desulfoglucosinolates formed during the enzyme assay, whereas the later eluting peaks show desulfoglucosinolates formed from intact glucosinolates via in-source fragmentation during MS analysis and so actually indicate the presence of intact sulfated glucosinolates used as substrates. (E) The time-dependent degradation of a mixture of intact glucosinolates from Arabidopsis thaliana Col-0, with indolic glucosinolates being depleted more quickly than aliphatic ones (amount of each glucosinolate present at 0 min was taken as 100%, N = 2 independent reactions). “ds-” = desulfated glucosinolates; 3msop = 3-methylsulfinylpropyl; 4msob = 4-methylsulfinylbutyl; 4mtb = 4-methylthiobutyl; 5msop = 5-methylsulfinylpentyl; 7msoh = 7-methylsulfinylheptyl; 8msoo = 8-methylsulfinyloctyl; i3m = indol-3-ylmethyl; 4moi3m = 4-methoxyindol-3-ylmethyl.
<italic>BtGSS</italic> Expression Is Independent of Glucosinolate Ingestion
To compare the relative expression levels of BtGSS in B. tabaci reared on non-glucosinolate diet (eggplant) and those reared on a glucosinolate-rich diet (Brussels sprouts), we analyzed the expression of BtGSS by qPCR in B. tabaci MEAM1 adult insects grown on both plants. BtGSS was constitutively and highly expressed in insects grown on both plant species tested (Supplementary Figure 2).
Silencing of BtGSS Supports Its Role as the Major GSS in <italic>Bemisia tabaci</italic>
To examine whether BtGSS acts to desulfate glucosinolates in vivo, we manipulated the expression of its encoding gene BtGSS using an RNAi approach. For this purpose, dsRNA targeting BtGSS was synthesized using in vitro transcription and fed to B. tabaci adults in artificial diets. Since non-specific dsRNases in B. tabaci could cause dsRNA degradation in the whiteflies, artificial diets also contained dsRNA targeting dsRNase2, which had been previously shown to enhance the silencing of genes of interest (Luo et al., 2017). This strategy led to a decrease of 43% in the expression of BtGSS in B. tabaci compared with a scrambled dsRNA control (Figure 4A). After dsRNA treatment, the proportion of 4msob-desulfoglucosinolate relative to non-metabolized intact 4msob-glucosinolate in honeydew samples was 48% lower in BtGSS-silenced adults than in scrambled dsRNA-fed control insects (Figure 4B). Both the extent of BtGSS gene silencing and the differences in desulfoglucosinolate formation indicated a high variability among experimental groups. However, these two factors were linearly correlated, with higher BtGSS expression corresponding to higher proportions of desulfoglucosinolate products in honeydew (Figure 4C).
BtGSS silencing in vivo. (A) Dietary administration of dsRNA (N = 3 per treatment; statistical differences determined via t-tests) targeting BtGSS led to diminished expression of this gene (relative to rpl-13) compared with a treatment using a scrambled dsRNA control. (B) Silencing reduced the proportion of desulfo-4msob glucosinolate relative to intact glucosinolate present in whitefly honeydew. (C)BtGSS expression was positively correlated to the proportions of desulfo-4msob glucosinolate:intact glucosinolate in honeydew.
To compare BtGSS activity toward glucosinolates with different side chains as substrates, the enzyme was incubated for different periods with pure glucosinolates and with a glucosinolate extract obtained from A. thaliana Col-0 leaves, which contained both aliphatic and indolic glucosinolates (Brown et al., 2003). BtGSS desulfated all glucosinolates offered, both as pure compounds and within the A. thaliana extract (Figures 3A–D and Supplementary Figure 1A–H). However, indolic glucosinolates were degraded much more rapidly than aliphatic glucosinolates (Figure 3E). All glucosinolates were stable under non-enzymatic control conditions (Supplementary Figure 3). This confirmed that BtGSS acts on different glucosinolate classes and suggested additionally that it might have distinct substrate affinities/catalytic efficiencies toward them. We confirmed this using two selected pure glucosinolates for in vitro incubations with BtGSS, 4msob- and i3m-glucosinolates. BtGSS desulfated i3m-glucosinolate, the most abundant indolyl glucosinolate in A. thaliana Col-0, with a catalytic efficiency of 466 nmol/mg BtGSS/min. This was 2.93 times more efficient than the conversion of 4msob-glucosinolate (159 nmol/mg BtGSS/min), its most abundant aliphatic counterpart in A. thaliana and several Brassica crops.
Discussion
Although cruciferous plants possess a sophisticated glucosinolate–myrosinase defense system to thwart their herbivores, several insects possess mechanisms to detoxify these plant defensive metabolites. One such strategy is to employ a GSS enzymatic system to convert glucosinolates to their corresponding desulfated forms, and thus prevent them from being broken down to toxic downstream products by myrosinases (Ratzka et al., 2002; Falk and Gershenzon, 2007). Interestingly, while all previously described GSS activities had come from leaf-chewing herbivores, this detoxification mechanism was recently reported from the phloem-feeding herbivore B. tabaci (Malka et al., 2016). It was once believed that phloem-feeding insects, by delicate piercing of plant tissues and stealthy sucking of the plant sap, do not strongly activate the glucosinolate–myrosinase system (Walling, 2008; Pentzold et al., 2014). However, more recent studies showed that phloem feeding indeed leads to glucosinolate activation and production of isothiocyanates (Walling, 2008; Danner et al., 2018). These recent results clarify the advantage of efficiently and pre-emptively detoxifying glucosinolates in piercing–sucking insects, such as whiteflies. As the B. tabaci species complex is a major agricultural pest, it was of interest to identify the GSS(s) performing glucosinolate detoxification in this insect.
We successfully identified one GSS (BtGSS) in B. tabaci MEAM1 that actively desulfated different glucosinolates in vitro. The heterologously produced enzyme displayed optimum GSS activity at pH 7.5–7.8 and was secreted into the extracellular culture medium by Sf9 cells, in agreement with the presence of a native signal peptide guiding BtGSS for secretion. Additionally, available transcriptomic data suggested an elevated expression of this enzyme in the B. tabaci gut tissues relative to the rest of the body in both MEAM1 and MED B. tabaci species (formerly named “B” and “Q” biotypes, respectively). These data together suggest that this enzyme is secreted into the B. tabaci gut lumen, where it can efficiently detoxify potentially harmful glucosinolates after ingestion. This modus operandi mimics that of the P. xylostella GSS, which is secreted from gut epithelial cells and is present in high amounts in the larval gut lumen (Ratzka et al., 2002). Nevertheless, the B. tabaci BtGSS we identified is not the sulfatase with closest homology to the P. xylostella GSS in the whitefly. The sulfatase designated BtSulf1 is phylogenetically closest to PxGSS but did not have GSS activity in vitro. These data suggest that GSS activity evolved independently in whiteflies and lepidopterans from sulfatases performing different innate functions in these distant insect lineages. In fact, BtGSS was the only sulfatase for which GSS activity was detected among all BtSulf candidates tested in this study. These included BtSulf 7, 8, and 11, whose sequences were most closely related to BtGSS, but did not show GSS activity in vitro. Additionally, RNAi-mediated silencing of BtGSS in artificial diet-fed whiteflies led to a decrease in the proportion of desulfoglucosinolates in honeydew relative to unmetabolized intact glucosinolates. This suggests that BtGSS is indeed the major sulfatase responsible for glucosinolate desulfation in B. tabaci MEAM1.
We found that BtGSS exhibits a clear preference in vitro toward indolic glucosinolates rather than their aliphatic counterparts in A. thaliana Col-0. This result is in accordance with the higher proportion of indolic glucosinolates found to be desulfated relative to their aliphatic counterparts in honeydew of B. tabaci feeding on Brussels sprouts and A. thaliana (Malka et al., 2016). Thus, B. tabaci might employ its GSS detoxification machinery to specifically combat the induction and activation of indolic glucosinolates that takes place during phloem sap-feeding (Kim et al., 2008). This substrate preference also complements well another preemptive glucosinolate detoxification mechanism in B. tabaci, the formation of glucosylated glucosinolates. Analysis of honeydew from B. tabaci feeding on A. thaliana showed that glucosylation was found to occur mostly to aliphatic glucosinolates and few glucosylated indolic glucosinolate derivatives were detected (Malka et al., 2020).
The existence of efficient pre-emptive detoxification pathways in B. tabaci can help explain the apparent lack of preference of MEAM1 and MED whiteflies for plants with lower vs. higher glucosinolate contents (Yang et al., 2020); that is, as these B. tabaci species efficiently deactivate glucosinolates, they can feed on plants producing natural concentrations of glucosinolates with impunity. However, when glucosinolate production was increased above natural levels in transgenic lines, B. tabaci development and performance was negatively affected (Elbaz et al., 2012; Markovich et al., 2013). Moreover, myb51, a positive regulator of indolic glucosinolate production, is up-regulated in A. thaliana after infestations by the phloem-feeding insects B. tabaci, M. persicae, and Brevicoryne brassicae (Foyer et al., 2015). Taken together, these findings suggest that plants specifically induce the production of indolic glucosinolates as a response against phloem-feeding insects. These glucosinolates could be activated by one of the recently characterized non-canonical glucosinolate-β-thioglucosidases, such as PEN2 or PYK10 (Bednarek et al., 2009; Clay et al., 2009; Nakano et al., 2017), which may be more specific to indolic glucosinolates. Plants can also convert indol-3-ylmethyl glucosinolate into more bioactive methoxylated forms as a defensive response (Kim and Jander, 2007; Wiesner et al., 2013; Pfalz et al., 2016). The enzymatic breakdown products of indolic glucosinolates are said to have strong antifeedant and toxic effects on crucifer-consuming insects (Kim and Jander, 2007; Agerbirk et al., 2008; Kim et al., 2008).
Quantitative real-time PCR analysis showed that the gene encoding the identified B. tabaci GSS (BtGSS) is highly and constitutively expressed irrespective of whether the insect consumed glucosinolates or not. This suggests that this GSS could have other function(s) for the insect, perhaps serving as a more general arylsulfatase. While the in vivo functions of insect (aryl)sulfatases are still not well-defined, one of their proposed roles is during the molting process, for example, of southern armyworm larvae (Yang et al., 1973). It is known that sulfate conjugates of ecdysone, the steroidal insect molting hormone, are less bioactive than non-sulfated analogs and arylsulfatase titers increase upon molting (King, 1972; Karlson and Koolman, 1973; Yang et al., 1973). Arylsulfatases could, therefore, be involved in converting ecdysone to its active form. However, ecdysteroid localization is typically limited to ecdysteroidogenic glands and hemolymph (Gilbert et al., 2002). According to previous studies of the P. xylostella GSS and S. gregaria GSS, and in accordance with what we observed for B. tabaci, GSS is secreted into the gut lumen (Ratzka et al., 2002; Falk and Gershenzon, 2007). Therefore, it seems unlikely that BtGSS is involved in ecdysone metabolism, and its other potential roles in vivo remain unknown.
In this work, we have characterized the BtGSS gene expression levels and protein activities in B. tabaci of the MEAM1 species. Desulfoglucosinolates had, however, been previously reported in B. tabaci MED whiteflies (Malka et al., 2016), another aggressive, generalist B. tabaci species. Whether MED and other B. tabaci species might have different GSS activity levels, and how this might affect their establishment on crucifers, is not yet known. Nevertheless, it seems clear that B. tabaci uses several different biochemical strategies to deal with glucosinolates. Besides the serial glucosylation of glucosinolates (Malka et al., 2020), a glutathione-S-transferase enzyme has also been reported in B. tabaci (Eakteiman et al., 2018) that can metabolize the hydrolysis products of both major aliphatic and indolic glucosinolates present in A. thaliana. This latter activity might help to mitigate the potential negative effects caused by the small amounts of glucosinolates successfully activated during and after ingestion by neutralizing the toxic hydrolysis products. Whether other strategies might also be employed by this whitefly, as in another generalist phloem-feeding insect, M. persicae (Kim and Jander, 2007; Kim et al., 2008), remains to be determined.
Conclusion
The glucosinolate–myrosinase system is a sophisticated two-component chemical defense strategy present in plants of the order Brassicales. The constituents of this binary defense system are kept spatially segregated, but during herbivory, myrosinases catalyze the hydrolysis of glucosinolates forming toxic breakdown products. In contrast to leaf-chewing insects, which cause extensive damage during herbivory and activate the glucosinolate–myrosinase system, it was long assumed that phloem-feeding insects did not cause sufficient damage for activation. However, glucosinolate hydrolysis products have been detected during attack by aphids; and these hydrolysis products negatively affect phloem-feeding insects. Moreover, plants induce the production of glucosinolates (especially indolic glucosinolates) upon infestation by phloem feeders. Therefore, the strategies used by phloem-feeding insects that feed on glucosinolate-containing plants to circumvent these defensive metabolites deserve further attention.
The whitefly B. tabaci is one such crucifer-colonizing, phloem-feeding insect, which has been reported to efficiently detoxify glucosinolates by conversion to their corresponding desulfated forms. This species complex is a group of serious agricultural pests that causes extensive damage to crops including cruciferous vegetables, and also vectors many plant pathogenic viruses. Here, we have successfully identified and characterized a B. tabaci MEAM1 sulfatase that is responsible for this metabolic reaction. This enzyme exhibits a strong preference in vitro toward indolic glucosinolates in comparison with aliphatic glucosinolates, implying that B. tabaci employs this GSS detoxification mechanism to cope with the negative effects of indolic glucosinolate breakdown. This strategy seems to complement well the other glucosinolate detoxification pathways present in B. tabaci (glucosylation and glutathione conjugation) and helps explain why this insect can feed widely on cruciferous plants. However, the relative importance of the individual pathways, their costs, and how each contributes to the success of B. tabaci still remains to be determined.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
DV designed this study with input from AM, SW, and JG. AM, KL, BI, MG, ES, MLAEE, and BS performed the data collection and analysis. RK and MR provided additional materials, expertise, and technical help. SW, JG, and DV supervised the study. AM and DV drafted the manuscript with input from all coauthors. All authors contributed to the article and approved the submitted version.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We would like to acknowledge the Max Planck Society for funding this research and the publication of this manuscript. We thank Osnat Malka and Shai Morin (The Hebrew University of Jerusalem) for generously providing B. tabaci cDNA.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.671286/full#supplementary-material
ReferencesAgerbirkN.De VosM.KimJ. H.JanderG. (2008). Indole glucosinolate breakdown and its biological effects.Phytochem. Rev.8:101. 10.1007/s11101-008-9098-0AhnS. J.BetzinF.GikonyoM. W.YangZ. L.KollnerT. G.BeranF. (2019). Identification and evolution of glucosinolate sulfatases in a specialist flea beetle.Sci. Rep.9:15725. 10.1038/s41598-019-51749-x31673017BednarekP.Piślewska-BednarekM.SvatošA.SchneiderB.DoubskýJ.MansurovaM. (2009). A glucosinolate metabolism pathway in living plant cells mediates broad-spectrum antifungal defense.Science323101–106. 10.1126/science.116373219095900BlaževićI.MontautS.BurčulF.OlsenC. E.BurowM.RollinP. (2020). Glucosinolate structural diversity, identification, chemical synthesis and metabolism in plants.Phytochemistry169:112100. 10.1016/j.phytochem.2019.11210031771793BoykinL. M.ArmstrongK. F.KubatkoL.De BarroP. (2012). Species delimitation and global biosecurity.Evol. Bioinform.81–37. 10.4137/EBO.S853222267902BoykinL. M.KineneT.WainainaJ. M.SavillA.SealS.MugerwaH. (2018). Review and guide to a future naming system of African Bemisia tabaci species.Syst. Entomol.43427–433. 10.1111/syen.12294BrownP. D.TokuhisaJ. G.ReicheltM.GershenzonJ. (2003). Variation of glucosinolate accumulation among different organs and developmental stages of Arabidopsis thaliana.Phytochemistry62471–481. 10.1016/S0031-9422(02)00549-6ChenW.HasegawaD. K.KaurN.KliotA.PinheiroP. V.LuanJ. (2016). The draft genome of whitefly Bemisia tabaci MEAM1, a global crop pest, provides novel insights into virus transmission, host adaptation, and insecticide resistance.BMC Biol.14:110. 10.1186/s12915-016-0321-y27974049ClayN. K.AdioA. M.DenouxC.JanderG.AusubelF. M. (2009). Glucosinolate metabolites required for an Arabidopsis innate immune response.Science32395–101. 10.1126/science.116462719095898DannerH.DesurmontG. A.CristescuS. M.van DamN. M. (2018). Herbivore-induced plant volatiles accurately predict history of coexistence, diet breadth, and feeding mode of herbivores.New Phytol.220726–738. 10.1111/nph.1442828134434EakteimanG.Moses-KochR.MoshitzkyP.Mestre-RinconN.VassãoD. G.LuckK. (2018). Targeting detoxification genes by phloem-mediated RNAi: a new approach for controlling phloem-feeding insect pests.Insect Biochem. Mol. Biol.10010–21. 10.1016/j.ibmb.2018.05.00829859812EdgarR. C. (2004). MUSCLE: multiple sequence alignment with high accuracy and high throughput.Nucleic Acids Res.321792–1797. 10.1093/nar/gkh34015034147ElbazM.HalonE.MalkaO.MalitskyS.BlumE.AharoniA. (2012). Asymmetric adaptation to indolic and aliphatic glucosinolates in the B and Q sibling species of Bemisia tabaci (Hemiptera: Aleyrodidae.Mol. Ecol.214533–4546. 10.1111/j.1365-294X.2012.05713.x22849567FalkK. L.GershenzonJ. (2007). The desert locust, Schistocerca gregaria, detoxifies the glucosinolates of Schouwia purpurea by desulfation.J. Chem. Ecol.331542–1555. 10.1007/s10886-007-9331-017619221FoyerC. H.VerrallS. R.HancockR. D. (2015). Systematic analysis of phloem-feeding insect-induced transcriptional reprogramming in Arabidopsis highlights common features and reveals distinct responses to specialist and generalist insects.J. Exp. Bot.66495–512. 10.1093/jxb/eru49125540442GilbertL. I.RybczynskiR.WarrenJ. T. (2002). Control and biochemical nature of the ecdysteroidogenic pathway.Annu. Rev. Entomol.47883–916. 10.1146/annurev.ento.47.091201.14530211729094GraserG.SchneiderB.OldhamN. J.GershenzonJ. (2000). The methionine chain elongation pathway in the biosynthesis of glucosinolates in Eruca sativa (Brassicaceae).Arch. Biochem. Biophys.378411–419. 10.1006/abbi.2000.181210860559GrobK.MatileP. (1979). Vacuolar location of glucosinolates in horseradish root-cells.Plant Sci. Lett.14327–335. 10.1016/S0304-4211(79)90281-5HalkierB. A. (2016). General introduction to glucosinolates.Adv. Bot. Res.801–14. 10.1016/bs.abr.2016.07.001HalkierB. A.GershenzonJ. (2006). Biology and biochemistry of glucosinolates.Annu. Rev. Plant Biol.57303–333. 10.1146/annurev.arplant.57.032905.10522816669764HöglundA. S.LenmanM.FalkA.RaskL. (1991). Distribution of myrosinase in rapeseed tissues.Plant Physiol.95213–221. 10.1104/pp.95.1.21316667954HopkinsR. J.van DamN. M.van LoonJ. J. A. (2009). Role of glucosinolates in insect-plant relationships and multitrophic interactions.Annu. Rev. Entomol.5457–83. 10.1146/annurev.ento.54.110807.09062318811249JeschkeV.GershenzonJ.VassãoD. G. (2016a). A mode of action of glucosinolate-derived isothiocyanates: detoxification depletes glutathione and cysteine levels with ramifications on protein metabolism in Spodoptera littoralis.Insect Biochem. Mol. Biol.7137–48. 10.1016/j.ibmb.2016.02.00226855197JeschkeV.GershenzonJ.VassãoD. G. (2016b). Insect detoxification of glucosinolates and their hydrolysis products.Adv. Bot. Res.80199–245. 10.1016/bs.abr.2016.06.003JingX.WhiteT. A.LuanJ.JiaoC.FeiZ.DouglasA. E. (2016). Evolutionary conservation of candidate osmoregulation genes in plant phloem sap-feeding insects.Insect Mol. Biol.25251–258. 10.1111/imb.1221526896054KarlsonP.KoolmanJ. (1973). On the metabolic fate of ecdysone and 3-dehydroecdysone in Calliphora vicina.Insect Biochem.3409–417. 10.1016/0020-1790(73)90074-7KawakishiS.KanekoT. (1987). Interaction of proteins with allyl isothiocyanate.J. Agric. Food Chem.3585–88. 10.1021/jf00073a020KimJ. H.JanderG. (2007). Myzus persicae (green peach aphid) feeding on Arabidopsis induces the formation of a deterrent indole glucosinolate.Plant J.491008–1019. 10.1111/j.1365-313X.2006.03019.x17257166KimJ. H.LeeB. W.SchroederF. C.JanderG. (2008). Identification of indole glucosinolate breakdown products with antifeedant effects on Myzus persicae (green peach aphid.Plant J.541015–1026. 10.1111/j.1365-313X.2008.03476.x18346197KingD. S. (1972). Ecdysone metabolism in insects.Am. Zool.12343–345.KorolevaO. A.DaviesA.DeekenR.ThorpeM. R.TomosA. D.HedrichR. (2000). Identification of a new glucosinolate-rich cell type in Arabidopsis flower stalk.Plant Physiol.124599–608. 10.1104/pp.124.2.59911027710LuanJ. B.ChenW. B.HasegawaD. K.SimmonsA. M.WintermantelW. M.LingK. S. (2015). Metabolic coevolution in the bacterial symbiosis of whiteflies and related plant sap-feeding insects.Genome Biol. Evol.72635–2647. 10.1093/gbe/evv17026377567LuoY. A.ChenQ. G.LuanJ. B.ChungS. H.Van EckJ.TurgeonR. (2017). Towards an understanding of the molecular basis of effective RNAi against a global insect pest, the whitefly Bemisia tabaci.Insect Biochem. Mol. Biol.8821–29. 10.1016/j.ibmb.2017.07.00528736300LüthyB.MatileP. (1984). The mustard oil bomb: rectified analysis of the subcellular organisation of the myrosinase system.Biochem. Physiol. Pflanzen1795–12. 10.1016/s0015-3796(84)80059-1MalkaO.EassonM. L. A. E.PaetzC.GötzM.ReicheltM.SteinB. (2020). Glucosylation prevents plant defense activation in phloem-feeding insects.Nat. Chem. Biol.691795–1797. 10.1038/s41589-020-00658-632989301MalkaO.ShekhovA.ReicheltM.GershenzonJ.VassãoD. G.MorinS. (2016). Glucosinolate desulfation by the phloem-feeding insect Bemisia tabaci.J. Chem. Ecol.42230–235. 10.1007/s10886-016-0675-126961756MarkovichO.KafleD.ElbazM.MalitskyS.AharoniA.SchwarzkopfA. (2013). Arabidopsis thaliana plants with different levels of aliphatic- and indolyl-glucosinolates affect host selection and performance of Bemisia tabaci.J. Chem. Ecol.391361–1372. 10.1007/s10886-013-0358-024190022MithöferA.BolandW. (2012). Plant defense against herbivores: chemical aspects.Annu. Rev. Plant Biol.63431–450. 10.1146/annurev-arplant-042110-10385422404468NakanoR. T.Pislewska-BednarekM.YamadaK.EdgerP. P.MiyaharaM.KondoM. (2017). PYK10 myrosinase reveals a functional coordination between endoplasmic reticulum bodies and glucosinolates in Arabidopsis thaliana.Plant J.89204–220. 10.1111/tpj.1337727612205OliveiraM. R. V.HenneberryT. J.AndersonP. (2001). History, current status, and collaborative research projects for Bemisia tabaci.Crop Protection20709–723. 10.1016/S0261-2194(01)00108-9PentzoldS.ZagrobelnyM.RookF.BakS. (2014). How insects overcome two-component plant chemical defence: plant β-glucosidases as the main target for herbivore adaptation.Biol. Rev.89531–551. 10.1111/brv.1206625165798PerringT. M. (2001). The Bemisia tabaci species complex.Crop Protection20725–737. 10.1016/S0261-2194(01)00109-0PfafflM. W. (2001). A new mathematical model for relative quantification in real-time RT–PCR.Nucleic Acids Res.29:e45.PfalzM.MukhaimarM.PerreauF.KirkJ.HansenC. I. C.OlsenC. E. (2016). Methyl transfer in glucosinolate biosynthesis mediated by indole glucosinolate O-methyltransferase 5.Plant Physiol.1722190–2203. 10.1104/pp.16.0140227810943RaskL.AndreassonE.EkbomB.ErikssonS.PontoppidanB.MeijerJ. (2000). Myrosinase: gene family evolution and herbivore defense in Brassicaceae.Plant Mol. Biol.4293–114. 10.1023/A:1006380021658RatzkaA.VogelH.KliebensteinD. J.Mitchell-OldsT.KroymannJ. (2002). Disarming the mustard oil bomb.Proc. Natl. Acad. Sci. U.S.A.9911223–11228. 10.1073/pnas.17211289912161563SchrammK.VassãoD. G.ReicheltM.GershenzonJ.WittstockU. (2012). Metabolism of glucosinolate-derived isothiocyanates to glutathione conjugates in generalist lepidopteran herbivores.Insect Biochem. Mol. Biol.42174–182. 10.1016/j.ibmb.2011.12.00222193392SchweizerF.Heidel-FischerH.VogelH.ReymondP. (2017). Arabidopsis glucosinolates trigger a contrasting transcriptomic response in a generalist and a specialist herbivore.Insect Biochem. Mol. Biol.8521–31. 10.1016/j.ibmb.2017.04.00428455184VogelH.BadapandaC.KnorrE.VilcinskasA. (2014). RNA-sequencing analysis reveals abundant developmental stage-specific and immunity-related genes in the pollen beetle Meligethes aeneus.Insect Mol. Biol.2398–112. 10.1111/imb.1206724252113WallingL. L. (2008). Avoiding effective defenses: strategies employed by phloem-feeding insects.Plant Physiol.146859–866. 10.1104/pp.107.11314218316641WiesnerM.HanschenF. S.SchreinerM.GlattH.ZrennerR. (2013). Induced production of 1-methoxy-indol-3-ylmethyl glucosinolate by jasmonic acid and methyl jasmonate in sprouts and leaves of pak choi (Brassica rapa) ssp. chinensis.Int. J. Mol. Sci.1414996–15016. 10.3390/ijms14071499623873294WittstockU.AgerbirkN.StauberE. J.OlsenC. E.HipplerM.Mitchell-OldsT. (2004). Successful herbivore attack due to metabolic diversion of a plant chemical defense.Proc. Natl. Acad. Sci. U.S.A.1014859–4864. 10.1073/pnas.030800710115051878YangJ.XieW.LiuB.WangS.WuQ.HeY. (2020). Phenolics, rather than glucosinolates, mediate host choice of Bemisia tabaci MEAM1 and MED on five cabbage genotypes.J. Appl. Entomol.144287–296. 10.1111/jen.12737YangR. S.PellicciaJ. G.WilkinsonC. F. (1973). Age-dependent arylsulphatase and sulphotransferase activities in the southern armyworm: a possible insect endocrine regulatory mechanism?Biochem. J.136817–820. 10.1042/bj13608174780704YeX. D.SuY. L.ZhaoQ. Y.XiaW. Q.LiuS. S.WangX. W. (2014). Transcriptomic analyses reveal the adaptive features and biological differences of guts from two invasive whitefly species.BMC Genomics15:370. 10.1186/1471-2164-15-37024885120
http://iae.fafu.edu.cn/DBM/genePage.php?id=Px018104, downloaded on 01/15/2017, and http://www.insect-genome.com/search/gene_detail1.php?id=790752, downloaded on 05/18/2021