Plant Defensive β-Glucosidases Resist Digestion and Sustain Activity in the Gut of a Lepidopteran Herbivore.

Two-component activated chemical defenses are a major part of many plants' strategies to disrupt herbivory. The activation step is often the β-glucosidase-catalyzed removal of a glucose moiety from a pro-toxin, leading to an unstable and toxic aglycone. While some β-glucosidases have been well studied, several aspects of their roles in vivo, such as their precise sites of enzymatic activity during and after ingestion, and the importance of particular isoforms in plant defense are still not fully understood. Here, plant defensive β-glucosidases from maize, white mustard and almonds were shown to resist digestion by larvae of the generalist lepidopteran Spodoptera littoralis, and the majority of the ingested activities toward both general and plant pro-toxic substrates was recovered in the frass. Among other proteins potentially involved in defense, we identified specific plant β-glucosidases and a maize β-glucosidase aggregating factor in frass from plant-fed insects using proteomic methods. We therefore found that, while S. littoralis larvae efficiently degraded bulk food protein during digestion, β-glucosidases were among a small number of plant defensive proteins that resist insect digestive proteolysis. These enzymes remain intact in the gut lumen and frass and can therefore further catalyze the activation of plant defenses after ingestion, especially in pH-neutral regions of the digestive system. As most of the ingested enzymatic activity persists in the frass, and only particular β-glucosidases were detected via proteomic analyses, our data support the involvement of specific isoforms (maize ZmGlu1 and S. alba MA1 myrosinase) in defense in vivo.

Introduction

Many plants utilize activated plant chemical defenses to fend off herbivore attacks. These small molecular weight compounds are stored as biologically inactive pro-toxins in the intact tissue, and are then enzymatically activated to form bioactive toxic compounds after tissue damage by an attacking herbivore (Halkier and Gershenzon, 2006; Morant et al., 2008). This two-component system avoids the problem of auto-toxicity, as these defensive molecules can be stored without ill effects to the plants. A very common activation strategy is the enzymatic removal of a protecting glucose group by a β-glucosidase (Morant et al., 2008). The resulting aglycone is typically more reactive than the parent glucoside and accordingly displays increased toxicity. Furthermore, its lower polarity facilitates diffusion through cell membranes and penetration into the herbivore cells where toxicity is often exerted. In order to avoid uncontrolled and unnecessary activation, the two components of the defensive system (the glucosylated pro-toxin and the activating β-glucosidase) are kept spatially separated. These meet only when tissue integrity is lost, leading to mixing of the contents of different cells and cellular compartments. The increased water-solubility conferred by glycosylation allows for storage of millimolar concentrations of compounds such as glucosinolates, benzoxazinoids, and cyanogenic glucosides. While much is known about the chemistry and compartmentalization of activated defenses, it is not clear whether they are activated only on plant damage or also later in the herbivore gut.

Among activated plant chemical defenses, one of the best-studied systems involves glucosinolates and myrosinases that are mostly restricted to brassicaceous plants such as broccoli, mustard, and rapeseed (Halkier and Gershenzon, 2006). Hydrolysis of glucosinolates results in the liberation of glucose and an unstable N-(sulfooxy)alkylimidothioic acid aglycone that undergoes rearrangement in aqueous media to form compounds of varying toxicity. These products include isothiocyanates (mustard oils), nitriles, epithionitriles, and thiocyanates (Figure ) depending on the reaction conditions (e.g., pH, concentrations of metal cations and plant-derived specifier proteins) and the nature of the glucosinolate side chain (Bones and Rossiter, 2006; Wittstock and Burow, 2010; Jeschke et al., 2016b). Among those, isothiocyanates are considered particularly toxic, reacting with biological nucleophiles such as glutathione (GSH) and certain amino acid side-chains in proteins (Bruggeman et al., 1986; Mi et al., 2011; Mithöfer and Boland, 2012; Schramm et al., 2012; Jeschke et al., 2016a). The sugar moiety in glucosinolates is connected to the aglycone via an S-β-glucosidic bond. The β-glucosidases that activate glucosinolates are therefore β-thioglucoside glucohydrolases, commonly named myrosinases.

Schematic representation of some reactions catalyzed by plant β-glucosidases. (A) Activation of glucosinolates by myrosinase. (B) Step-wise hydrolysis of the cyanogenic diglucoside amygdalin. (C) Hydrolytic activation of the benzoxazinoid glucoside DIMBOA-Glc and formation of the open-ring aglycone form in solution.

Cyanogenic glucosides are similarly activated by removal of their sugar group(s), leading to α-cyanohydrin aglycones. For example, one sugar molecule from the diglucoside amygdalin, found especially in fruits of the Rosaceae, is removed by amygdalin hydrolases to give prunasin, and this cyanogenic monoglucoside is further hydrolyzed into mandelonitrile and glucose by prunasin hydrolases (Figure ). The mandelonitrile aglycone thus formed may then dissociate into benzaldehyde and the very toxic hydrogen cyanide (HCN) non-enzymatically, especially under high pH (Fomunyam et al., 1985), or with enzymatic catalysis by mandelonitrile lyase (Swain and Poulton, 1994a,b). Emission of HCN, or cyanogenesis, occurs in a surprisingly large number of plant species of importance to humans, being a feature of grain crops such as oats, sorghum and maize, fruit crops such as apples and mangoes, and in the sub-Saharan staple food crop cassava.

Finally, benzoxazinoids (BXDs) are indole-derived defense compounds produced by a small number of agriculturally relevant grass (Poaceae) crops, including maize, wheat and rye (Frey et al., 1997; Niemeyer, 2009), and in scattered dicotyledonous species. They are mostly accumulated in intact plant tissues as the corresponding BXD glucosides. These lead, after glucose hydrolysis, to cyclic hemiacetal aglycones that equilibrate to α-oxoaldehydes by ring opening (Figure ) and are reactive toward a wide range of biological nucleophiles (Maresh et al., 2006; Dixon et al., 2012). The major BXD glucoside in young maize leaves is (2R)-2-β-D-glucopyranosyloxy-4-hydroxy-7-methoxy-2H-1,4-benzoxazin-3(4H)-one (DIMBOA-Glc), whose aglycone DIMBOA is considered to be one of the primary maize defensive chemicals against generalist leaf-chewing herbivores (Cambier et al., 2000; Glauser et al., 2011).

Glycosidases are a generally well-studied class of enzymes, due in part to their industrial utility. Plant β-glucosidases are classified as part of the Family 1 of glycoside hydrolases, catalyzing the hydrolysis of a terminal β-D-glucose residue via a mechanism in which the anomeric configuration of the glucose is retained. These β-glucosidases generally have a subunit molecular mass of 55–65 kDa, slightly acidic pH optima (pH 5–6), and an absolute specificity toward β-glucosides (Fowler, 1993). The tertiary structure of plant β-glucosidases belonging to family 1 glycoside hydrolases is highly conserved, as are the active site amino acids involved in binding of the sugar moiety. The residues that define the aglycone binding pocket are found at conserved positions in the active site, but are variable. The regions involved in oligomerization also vary among β-glucosidases. This is reflected in the diversity of quaternary structures observed in plant β-glucosidases, where the different active enzymes can be monomers, dimers, tetramers, or even higher oligomers. While the glucoside pro-toxins are typically stored in the vacuole, the subcellular localization of activating β-glucosidases differs in monocotyledonous and dicotyledonous plants. Monocotyledonous β-glucosidases contain an N-terminal transit peptide leading to localization to plastids. In dicotyledons, these activating enzymes contain an N-terminal signal that results in co-translational glycosylation and secretion (Morant et al., 2008). Protein glycosylation promotes stability and solubility and mediates protein localization and protein–protein interactions (Nastruzzi et al., 1996; Burmeister et al., 1997), but is not essential for full enzymatic activity (Zhou et al., 2012).

Myrosinases (E.C. 3.2.1.147), the thioglucosidases responsible for hydrolytic activation of glucosinolates, are commonly present in multigene families, with members possibly having different roles in different tissues. The Arabidopsis thaliana genome encodes for six myrosinase genes (tgg1-6), with two apparently being pseudogenes, and two genes being expressed each in aerial parts (tgg1/2) and underground tissues (tgg4/5) (Andersson et al., 2009; Zhou et al., 2012). More recently, additional A. thaliana proteins (PEN2 and PYK10) have been described to have myrosinase-like activities restricted toward indolic glucosinolates (Bednarek et al., 2009; Nakano et al., 2014), but their specific functions in plant defense against herbivores are still unclear.

The widely commercially available β-glucosidase from almonds, which produce cyanogenic glycosides, has been studied biochemically in great detail. While it was shown to artificially induce herbivory-like volatile emission in cabbages (Mattiacci et al., 1995), its potential roles in planta also remain to be better characterized. It can be isolated in large amounts from sweet almonds as a mixture of isozymes, and it has been massively used as a model enzyme in the study of β-glucosidase mechanism, kinetics, and chemical inactivation [for example (Shulman et al., 1976; Grover and Cushley, 1977; He and Withers, 1997)]. Its stability also allows its use in the synthesis of glucosides in vitro in high concentrations of organic solvents, exploiting the reversibility of the glucoside hydrolysis reactions [for example (Vic et al., 1995; Ducret et al., 2006)]. Nevertheless, while this enzyme accepts many compounds in vitro, information about whether it hydrolyzes cyanogenic glucosides in vivo is lacking.

In maize, BXD β-glucosidase reactions can be performed by two enzymes, ZmGlu1 and ZmGlu2 (Czjzek et al., 2000), which share 90% sequence identity. ZmGlu1 and ZmGlu2 hydrolyze a broad spectrum of artificial and natural compounds in addition to the benzoxazinoid DIMBOA-Glc (Czjzek et al., 2000, 2001). The maize β-glucosidase Zm-p60.1, of which ZmGlu1 is the corresponding allozyme, has been implicated in the release of active cytokinins from inactive stored cytokinin-O-glucosides (Campos et al., 1993; Zouhar et al., 1999; Zouhar et al., 2001).

Some maize proteins have been shown to interact directly with β-glucosidases resulting in the formation of large aggregates. Biochemical and immunological studies demonstrated that β-glucosidase activity could not be detected in zymograms of certain maize genotypes (so-called “nulls”) due to the enzyme occurring as large insoluble or poorly soluble quaternary complexes (Esen and Blanchard, 2000). This complexation is mediated by a β-glucosidase-aggregating factor (BGAF), a chimeric maize protein consisting of an N-terminal dirigent (disease response) domain and a C-terminal jacalin-related lectin domain (Kittur et al., 2007). An analogous jacalin-like lectin protein from wheat (Triticum aestivum), HFR1, has direct antifeedant effects and is involved in resistance against the Hessian fly (Mayetiola destructor) (Subramanyam et al., 2008), but its individual protein binding partners were not examined. The enzymatic activities of ZmGlu1 and ZmGlu2 are unaffected by BGAF binding, suggesting that interaction with BGAF does not change the conformation of the enzymes (Blanchard, 2001). Similar myrosinase-binding proteins (MBP) have been found in rapeseed (Falk et al., 1995) and in A. thaliana (Takechi et al., 1999; Takeda et al., 2008), where they also lead to formation of higher-molecular weight active complexes. However, whether these protein aggregates are stable in herbivore guts after ingestion and whether they play any roles in the hydrolysis of activated defenses in vivo is not yet known.

Nitrogen in the form of proteins and essential amino acids is considered the most limiting macronutrient in herbivorous diets, justifying the evolution of efficient and flexible cocktails of digestive proteases in lepidopteran herbivores (Terra and Ferreira, 1994; Kuwar et al., 2015). Accordingly, most bulk plant protein is effectively hydrolyzed after ingestion by these insects (Chen et al., 2005) and can be used for insect growth. However, a growing number of plant proteins have been described in insect frass, several of which remaining catalytically or biologically active after passing through the insect digestive tract and acting in plant defense both during and after ingestion (Chen et al., 2005, 2007; Chuang et al., 2013; Ray et al., 2016).

In studying plant two-component activated defenses, it is important to learn if plant β-glucosidases and other hydrolyzing enzymes critical to defense survive initial proteolysis and are active in the insect gut. Recent experiments using insects fed on maize leaves showed that maize β-glucosidases involved in benzoxazinoid metabolism are active even after going through insect digestion (Glauser et al., 2011; Wouters et al., 2014). Hence we have now examined the stability of plant defensive β-glucosidases after digestion by a generalist insect herbivore (Spodoptera littoralis) in order to better understand their roles in defense. We investigated three different β-glucosidases said to be involved in benzoxazinoid, glucosinolate, and cyanogenic glycoside activation, respectively, by determining their presence in protein extracts made from insect frass and comparing their catalytic activity in non-ingested material and after passage through the gut.

Materials and Methods

Materials

Experiments with maize (Zea mays) used seeds of the variety Badischer Gelber from Kiepenkerl (Bad Marienberg, Germany) while white mustard (Sinapis alba) seeds were obtained from N. L. Chrestensen (Erfurt, Germany). Maize and white mustard were sown in plastic pots (9 cm × 9 cm) with two plants per pot under controlled light and temperature conditions (16:8 h light/dark, day-time temperature 22°C, night-time 20°C). White mustard thioglucosidase (180 U/g toward sinigrin), almond (Prunus dulcis) β-glucosidase (5.2 U/mg toward salicin), sinigrin (allyl glucosinolate), amygdalin and p-nitrophenyl-β-D-glucoside (pNPG) were obtained from Sigma-Aldrich (Munich, Germany). Larvae of S. littoralis (Lepidoptera: Noctuidae) were reared from eggs kindly provided by Syngenta Crop Protection (Stein, Switzerland) on artificial diet modified from Bergomaz and Boppré (1986) at 18°C under 12:12 h light cycle, and used at the fourth instar. The diet was prepared as follows: 2.4 g agar-agar was boiled with 80 mL tap water and added to a mixture of 25 g white bean flour, 450 mg ascorbic acid, 250 mg 4-ethylbenzoic acid and 25 mL tap water. The diet was stirred until hand-warm (∼40°C) and then supplemented with 45 mg α-tocopherol in 625 μL canola oil and 200 μL 3.7% formaldehyde.

Maize β-Glucosidase Extraction

About 100 maize kernels were washed with 80% ethanol for 2 min and 3% sodium hypochlorite with 0.5% Tween 20 for 15 min, then left in a Petri dish with water for 2 days to germinate, and sown as described above. Six- to seven-day-old seedlings were harvested, frozen and ground in liquid nitrogen and stored at -80°C prior to extraction. Enzyme extraction and cryoprecipitation were performed according to (Esen, 1992) with the final supernatant concentrated with Centricon Plus-70 regenerated cellulose membrane 30,000 NMWL ultracentrifugation filters (Merck Millipore, Hessen, Germany) for about 40 min. All operations were carried out between 0 and 4°C. The enzymatic activity of the extract was determined using the non-specific substrate pNPG as described below. The final enzymatic extract was adjusted to 1 U/mL.

DIMBOA-Glc Extraction

About 70 maize kernels were germinated as above and grown in the dark. Seven-day-old seedlings were harvested, weighed, and quickly frozen and ground in liquid nitrogen. Plant material was extracted with 3 mL of methanol per gram of ground plant material, quickly mixing the ground material in the solvent to avoid degradation. The suspension was centrifuged at 3,400 g for 10 min and the supernatant was collected. The extraction and centrifugation steps were repeated two more times, and the combined supernatants were filtered. Ammonium sulfate was added to the filtrate and left overnight. The supernatant was then filtered and evaporated under vacuum, and the solution was lyophilized. After resuspension at 200 mg/mL in 0.5% aqueous formic acid (FA)/MeOH (1:1) the compound was fractionated and purified using semi-preparative HPLC. Purification of DIMBOA-Glc used a Supelcosil LC-18-DB semi-preparative column (25 cm × 10 mm I.D., 5 μm particle size, Sigma-Aldrich, maintained at 25°C) with a flow rate of 4 mL/min, using FA (0.005%) in water and acetonitrile as mobile phases A and B, respectively. The gradient was as follows: 10–20% B (10 min), 100% B (3 min), 10% B (3 min) with detection by UV absorption at 254 nm. The solvents where evaporated partially by using a rotary evaporator, while the remaining water was removed by lyophilization, with 37.2 mg of DIMBOA-Glc obtained as a white powder.

Caterpillar Feeding and Collection of Diet and Frass Samples

In artificial diet feeding experiments, S. littoralis caterpillars previously fed on artificial diet were switched to a cube of artificial diet (0.5 cm × 0.5 cm) to which pure enzyme (10 μL 0.026 U/mL P. dulcis β-glucosidase or 10 μL 10 U/mL S. alba myrosinase) or the semi-purified maize enzyme extract (10 μL 1 U/mL) was applied (Figure ). Negative controls used 10 μL of the corresponding buffer. Each caterpillar was kept in an individual plastic cup and allowed to feed until the cube was consumed. An additional feeding experiment used bovine β-lactoglobulin (20 μL 10 mg/mL) added to the diet. Non-ingested samples used diet cubes with added enzyme kept under identical temperature conditions but in the absence of S. littoralis larvae. The resulting diet and frass samples were collected and extracted by adding two metal beads (diameter 3 mm) and 500 μL of the respective buffer used for each assay, followed by vigorous agitation with a paint shaker for 4 × 4 min, and centrifugation at 4,300 g at 4°C for 40 min. The supernatants were collected and assayed for enzymatic activity.

Schematic representation of enzyme feeding strategies from divided leaves (left) and in artificial diet cubes (right).

Feeding experiments with leaves instead of artificial diets were performed with maize and white mustard. For maize, one piece (about 3 cm × 1.5 cm) from a half of a leaf (divided longitudinally) from a 2-week-old plant was fed, with the other corresponding half frozen and used to determine non-ingested enzyme activity (Figure ). For white mustard, five leaves from a 1-week-old plant were fed per caterpillar, with a set containing another five leaves used to quantify the initial myrosinase activity being fed. Six replicates were made for each feeding experiment. After complete ingestion of the glucosidase-containing diet, a small amount of control diet was offered to the insect, to allow for excretion of the material still present in the insect gut.

Maize β-Glucosidase Activity Assay

In initial assays the non-specific β-glucosidase substrate p-nitrophenyl-β-D-glucopyranoside (pNPG) was used. For that, 500 μL of enzyme extract in 10 mM phosphate buffer at pH 7.1 were incubated with 500 μL 5 mM aqueous pNPG at 25°C for 5 min. The reaction was stopped by adding 1 mL 400 mM aqueous Na2CO3 and freezing in liquid nitrogen until further use, and the p-nitrophenol liberated was measured at 400 nm in a UV-Vis spectrophotometer shortly after thawing. Other enzyme assays used DIMBOA-Glc as substrate. The standard assay mixture for DIMBOA-Glc consisted of 5 μL enzyme extract (1 U/mL), 2 μL DIMBOA-Glc in DMSO (200 μM), and 43 μL 10 mM MOPS buffer pH 7.0. After incubation at 25°C for 20 min, 50 μL MeOH/FA (1:1) were added, and the samples were analyzed by LC-MS/MS as in (Wouters et al., 2014). Assays were performed in triplicate. Statistical testing for all enzyme activity results was performed with unpaired t-tests using the software GraphPad (La Jolla, CA, United States). β-Glucosidase activity was evident from both the decrease of the glucoside and concomitant appearance of the corresponding aglycone, which were quantified using external calibration curves.

Almond β-Glucosidase Enzymatic Assay

The assay mixture consisted of 250 μL enzyme extract, 500 μL 0.1 M acetate buffer (pH 5.0) and 250 μL 20 mM aqueous pNPG. The reaction mixture was incubated at 37°C for 15 min and analyzed as above. Assays were performed in triplicate.

S. alba Myrosinase Enzymatic Assay

The standard assay mixture consisted of 30 μL protein extract, 20 μL 50 mM aqueous sinigrin, and 200 μL of 50 mM citrate buffer (pH 5.5). After incubation for 60 min at 37°C, the reaction was stopped by addition of 25 μL glacial acetic acid and cooling on ice. The mixture was extracted with 250 μL CH2Cl2 containing 0.01% (v/v) hexyl-isothiocyanate as internal standard. A 100 μL aliquot of the CH2Cl2 fraction was removed and left overnight on 20 mg ammonium sulfate. Assays were performed in sextuplicate. The samples were analyzed by GC-MS, with enzymatic activities measured by quantification of released isothiocyanates using an external calibration curve of allyl-isothiocyanate.

Analysis of Assay Products

LC-MS/MS analyses of DIMBOA-Glc β-glucosidase activity were performed as described in (Wouters et al., 2014). GC-MS analyses were performed in an Agilent 6890 series gas chromatograph (Agilent Technologies, Waldbronn, Germany) using an Agilent 19091S-433 capillary column (30 m × 0.25 mm × 0.25 μm), using 2 μL injection on splitless mode at 200°C. The temperature program was: 40°C for 3 min, a 10°C∗min−1 ramp to 130°C, a 60°C∗min−1 ramp to 300°C, held for 3 min. The total running time was 17.83 min.

Spectrophotometry Analysis

A double-beam UV-Vis spectrophotometer (UV-2401PC, Shimadzu Corporation, Kyoto, Japan) equipped with integrating sphere was used for all spectrophotometry assays.

Proteomic Analysis of Gel-Separated Proteins

Frass extracts for proteomic analyses were prepared as described above. Protein was quantified using the Bradford method, and about 200 μg of protein per sample were concentrated by lyophilization and separated using an Any kD Mini-Protean TGX Gel (Bio-Rad).

In-Gel Digestion of Proteins

Protein bands between 20 and 120 kDa were cut out from the Coomassie-stained gels, cut into small pieces, washed several times with 25 mM NH4HCO3 and destained with 50% ACN/25 mM NH4HCO3. The proteins were then reduced in-gel with 10 mM dithiothreitol at 50°C for 1 h and alkylated with 55 mM iodoacetamide at room temperature in the dark for 45 min. Next, destained, washed, dehydrated gel pieces were rehydrated for 60 min in 0.5 μM solution of porcine trypsin (Promega) in 25 mM NH4HCO3 at 4°C and incubated overnight at 37°C. The tryptic peptides were extracted from gel pieces using 75% ACN/5% FA, and dried down in a SpeedVac. For LC-MS, analysis samples were reconstructed in 10 μL aqueous 1% FA (Shevchenko et al., 2007).

Nano-UPLC-MS Analysis of Peptides

One to 8 μL of the peptide mixture was injected onto a nanoAcquity nanoUPLC system (Waters) online coupled to a Q-ToF HDMS mass spectrometer (Waters). Samples were initially transferred with 0.1% aqueous FA for desalting onto a Symmetry C18 trap-column (20 mm × 0.18 mm, 5 μm particle size) at a flow rate of 15 μL/min (0.1% aqueous FA), and peptides were subsequently eluted onto a nanoAcquity C18 analytical column (200 mm × 75 μm I.D., BEH 130 material, 1.7 μm particle size) at a flow rate of 350 nL/min with the following gradient: 1–30% B over 13 min, 30–50% B over 5 min, 50–95% B over 5 min, 95% B for 4 min, and a return to 1% B over 1 min (phases A and B composed of 0.1% aqueous FA and 100% acetonitrile containing 0.1% FA, respectively). The analytical column was re-equilibrated for 9 min prior to the next injection. LC-MS data were acquired in positive ESI mode under data-independent acquisition (MSE) controlled by MassLynx v4.1 software. The collision energy was set at 4 eV in low energy (MS) scans, and ramped from 15 to 40 eV in elevated energy (MSE) scans. The mass range (m/z) for both scans was 300–1900 and 50–1700 Da, respectively. The scan time was set at 1.5 s for both modes of acquisition with an inter-scan delay of 0.2 s. A reference compound, human glu-fibrinopeptide B [650 fmol/mL in 0.1% aqueous FA/acetonitrile (v/v, 1:1)], was infused continuously through a reference sprayer for external calibration.

Data Processing and Protein Identification

ProteinLynx Global Server (PLGS) version 2.5.2 (Waters) was used for processing of raw files and for database searching. The continuum LC-MSE data were lock-mass corrected, smoothed, background subtracted, centered, deisotoped, and charge state reduced. Thresholds for low/high energy scan ions and peptide intensity were set at 150, 30, and 750 counts, respectively. Processed data were searched against the Swissprot database[1] (downloaded on July 27, 2011) combined with protein sequences of Z. mays[2] (downloaded on June 11, 2013) according to the described Ion Accounting Algorithm (Li et al., 2009). Database searches were performed at 4% FDR, using the following parameters: minimum number of product ion matches per peptide (5), minimum number of product ion matches per protein (7), minimum number of peptide matches (2), and maximum number of missed tryptic cleavage sites (1). Searches were restricted to tryptic peptides with a fixed carbamidomethyl modification for Cys residues, along with variable oxidation of Met. To avoid protein interference, identified hits were classified using a script based on the classification described for PAnalyzer software (Prieto et al., 2012). According to this algorithm, proteins were divided into four groups: conclusive, non-conclusive, indistinguishable, and ambiguous. Only proteins classified as conclusive hits (proteins with at least one unique peptide) were considered as confident matches.

Estimation of Gene Expression Levels

Gene expression was estimated using the maize eFP browser of the Bio-Analytic Resource at the University of Toronto[3] using the Sekhon et al. atlas setting (Winter et al., 2007; Sekhon et al., 2011) referring to data from B73 maize. Absolute expression values in the topmost leaf at growth stage V3 [“V3_Topmost leaf” (Sekhon et al., 2011)] are reported.

Results

Plants produce β-glucosidases that activate anti-herbivore defenses by cleaving various β-glucoside pro-toxins. In order to analyze the stability of these β-glucosidases during digestion by larvae from a generalist lepidopteran pest, we fed known quantities of these enzymes to S. littoralis caterpillars (Figure ) and quantified the enzymatic activities in non-ingested material and post-ingestion. Initial in vitro β-glucosidase assays using the general substrate p-nitrophenyl-β-D-glucoside (pNPG) revealed that non-specific β-glucosidases are present in frass from untreated larvae that arise from their diet and from the insect’s own endogenous activities, but these enzymes were not capable of cleaving specific plant defense compounds [such as DIMBOA-Glc (Glauser et al., 2011)]. Therefore, when using pNPG, the non-specific activities were measured (and corresponded to 5–30% of the added enzyme activities, depending on the enzyme being combined) and subtracted from activities measured in samples containing the added β-glucosidases.

Maize β-Glucosidase (Involved in Benzoxazinoid Hydrolysis) Was Active After Insect Digestion

Previous studies had shown that maize BXD-Glc β-glucosidases are active throughout the insect gut (Glauser et al., 2011), including within the hindgut and frass (Wouters et al., 2014). In order to quantify these observations, we first fed caterpillars on an artificial diet to which we applied a maize β-glucosidase enzyme extract in a similar amount/concentration to that in leaves. The resulting frass was then collected and extracted. Assays using pNPG indicated that maize β-glucosidases added to the diet were indeed excreted in active form, with 69.6% of the initial activity fed to the larvae being detected in the frass extracts (Figure ).

Maize β-glucosidases were recovered in active form after digestion by S. littoralis larvae. Enzymatic activity levels were determined in non-ingested material and after digestion of diet cubes spiked with semi-purified maize β-glucosidase (A), or after feeding on maize leaf tissue (B). Activities were measured either using the non-specific β-glucosidase substrate pNPG (A) or the endogenous plant substrate DIMBOA-Glc (B). Shown are the means ± standard errors [N = 3 in both analyses; P = 0.2 (A) and P = 0.3 (B)].

To assess how other plant-derived proteins might influence the digestion of these maize defensive β-glucosidases, caterpillars were subsequently fed a controlled amount of young maize leaf tissue. Leaves were divided along the mid-vein, and the matching leaf pieces (i.e., from the opposing sides of the leaf blade, see Figure ) were (a) fed to larvae and (b) used to quantify the initial β-glucosidase activity being fed. Quantification of activity was done with DIMBOA-Glc, a more specific and biologically relevant substrate than pNPG. We recovered 75.9% of the enzyme activity initially fed to the larvae, further confirming that the maize BXD β-glucosidases maintain most of their activity throughout the caterpillar gut (Figure ).

Myrosinase (Involved in Glucosinolate Hydrolysis) Was Active After Insect Digestion

Similarly, we quantified the activity of S. alba myrosinase remaining after digestion by S. littoralis larvae, using pure enzyme added to an artificial diet and by feeding with S. alba leaves. In both cases, enzymatic activities in non-ingested material and post-digestion were measured using a specific myrosinase substrate, sinigrin (allyl glucosinolate), with quantification of its main hydrolysis end-product, allyl isothiocyanate.

Myrosinase was also recovered in active form after digestion by the larvae. When the pure commercial enzyme was fed in artificial diet, 92.6% of the myrosinase activity administered was successfully recovered in the insect feces (Figure ), while 59.1% was recovered when comparing the activity in frass to that present in the corresponding unfed leaf (Figure ). The potential effects of the presence of ascorbate on the measured myrosinase activities were not tested, but this compound is known to increase S. alba myrosinase activity (Pihakaski and Pihakaski, 1978). Variation in its concentration between non-ingested material and post-digestion samples from plants and artificial diets might therefore account for some of the contrast in the apparent recovered activities.

Sinapis alba myrosinase activity was resistant to digestion by S. littoralis larvae and detected in larval frass. Enzymatic activity levels were determined in non-ingested material and after digestion of diet cubes spiked with pure enzyme (A) or of S. alba leaf material (B). Activities were measured using the myrosinase-specific substrate sinigrin. Shown are the means ± standard errors [N = 6 in both analyses; P = 0.2 (A) and P < 0.0001 (B)].

Almond β-Glucosidase Was Active After Insect Digestion

For quantification of the enzymatic activity of P. dulcis β-glucosidase after digestion by S. littoralis, pure commercial enzyme added to an artificial diet was used for the feeding experiment. Enzymatic assays using amygdalin did not lead to detectable amounts of its expected hydrolysis products. Enzymatic assays performed with the non-specific substrate pNPG, on the other hand, indicated that activity was quantitatively recovered after ingestion by the larvae (Figure ). The measured post-digestion activities were slightly higher than those for the undigested control, but the differences were not statistically significant (p = 0.09) and are likely due to natural variation in levels of endogenous insect β-glucosidase activities.

Prunus dulcis glucosidase activity was resistant to digestion by S. littoralis larvae and detected in larval frass. Enzymatic activity levels were determined in non-ingested material and after digestion of diet cubes spiked with pure enzyme. Activities were measured using the general β-glucosidase substrate pNPG. Shown are the means ± standard errors (N = 3; P = 0.09).

Proteomic Analyses

Protein extracts for proteomic analyses were prepared from maize and white mustard leaves, both in non-ingested material and post-ingestion. After separation by SDS–PAGE, gel regions were excised and digested with trypsin. The resulting peptides were analyzed by LC-MS using data-independent acquisition (DIA). In this approach, referred to as MSE, data were comprehensively acquired with an alternating mode of acquisition: a low energy (MS) mode that provided information on the mass of eluted intact peptide ions (precursors), and an elevated energy mode (MSE) when fragmentation of these ions occurred. The ion accounting algorithm implemented in the software PLGS (Waters) was applied for searching parallel fragmentation of multiple precursor ions, matching fragment ions to their corresponding precursor peptides while taking into account their chromatographic coelution in LC-MS/MS.

These analyses led, in the case of maize, to detection of six BXD β-glucosidase-derived peptides in frass material (Figure ). Two maize β-glucosidases have been identified to act on DIMBOA-Glc (ZmGlu1 and 2), but we were only able to confidently identify peptides derived from ZmGlu1. Seven peptides derived from the maize β-glucosidase aggregating factor (BGAF) protein were observed in frass extracts. Eleven peptides derived from S. alba myrosinase were found in frass extract samples, with unambiguous peptides indicating that the enzyme detected here is the isoform MA1. In contrast to the peptides belonging to the two maize proteins, the peptides corresponding to myrosinase were spread over a broad region of the SDS–PAGE gel (data not shown). This could result from partial deglycosylation/digestion of the proteins during passage through the gut.

LC-MSE-based proteomic identification of peptides (marked with an underline) derived from maize β-glucosidase ZmGlu1, maize β-glucosidase aggregating factor (BGAF), and S. alba myrosinase in frass samples.

Peptides belonging to 13 additional maize proteins were also detected in S. littoralis frass extracts (Supplementary Table ), including some with predicted functions matching those previously observed in the frass of maize-fed Spodoptera frugiperda larvae (Chuang et al., 2013). Expression levels of the corresponding genes in maize leaves (Sekhon et al., 2011) are scattered widely (Supplementary Table ) and do not correlate with numbers of detected peptides or sequence coverage, supporting a differential resistance to insect proteolysis among plant proteins. Many differences were observed in the general profiles of plant-derived polypeptides in total leaf extracts (Supplementary Table ) relative to those isolated from frass (Supplementary Table ). One of the major distinctions was in the large subunit (rbcL) of Rubisco (UniProt P00874) considered a marker for bulk leaf protein (Chen et al., 2005, 2007). In contrast to the good (48%) coverage of rbcL primary sequence in leaf extracts, our proteomic analyses did not detect any peptides derived from rbcL in frass extracts, supporting the degradation of bulk leaf protein during passage through the larval gut. In order to confirm this finding, we fed diet cubes containing bovine β-lactoglobulin in similar concentrations to those above for β-glucosidases, and analyzed non-ingested material and post-digestion samples. We could not detect a band corresponding to β-lactoglobulin with SDS–PAGE and Coomassie staining in insect-digested samples, and no peptides matching this protein could be found in frass proteomic data within the size range analyzed (∼7–35 kDa), again indicating that efficient proteolysis of bulk protein takes place during digestion by this herbivore.

Discussion

Previous work had identified defensive plant proteins in insect frass following passage through the insect digestive system (Chen et al., 2005, 2007; Chuang et al., 2013; Ray et al., 2016), but their biological activities were assessed only in a few specific cases. Here, we fed three β-glucosidases presumed to be involved in activating plant defenses to S. littoralis caterpillars and found that they remained active after passage through the gut. A maize β-glucosidase involved in benzoxazinoid hydrolysis, a S. alba myrosinase involved in glucosinolate activation, and an almond β-glucosidase associated with cyanogenic glycoside accumulation were all detected in frass using both in vitro enzymatic assays and proteomic analyses. Our proteomic analyses identified several additional proteins in extracts from frass of maize-fed larvae that could participate in plant defense, such as oxidases, proteases, and other glycoside hydrolases. Interestingly, both arginine decarboxylase and agmatine deiminase were observed, suggesting that these may function sequentially to deplete the essential amino acid arginine from the food bolus as part of a defensive strategy. Other predicted glycosyl hydrolases were also detected, and their utilization of benzoxazinoids as substrates and roles in defense remain to be explored.

The persistence of these proteins and corresponding activities in the frass means they are at least partially resistant to the proteases, high pH (9.5–11) and detergents present in S. littoralis guts. These results are in agreement with previous reports showing that other β-glucosidases, such as cassava cyanogenic β-glucosidase isozymes (linamarases) and maize and almond β-glucosidases are stable against short treatments with trypsin, high temperature and chemical denaturants (Mkpong et al., 1990; Esen and Gungor, 1993). This unusual resistance may be a general trait of plant β-glucosidases and is an important aspect of their commercial value as industrial biocatalysts. In contrast to β-glucosidases which remained active in the resulting frass, bulk dietary proteins like Rubisco and artificially added bovine β-lactoglobulin were efficiently degraded in the S. littoralis gut.

The biophysical proprieties that contribute to the high stability of β-glucosidases remain to be fully understood, but some structural features common among them may be responsible. Alignment and structural superimposition have revealed that, despite relatively low primary-structure similarities (17–44%), the tertiary structures of family 1 β-glucosidases are remarkably similar. Several β-glucosidases, including those studied here, have a tight and stable folded core, presumably enabling activity over a wide range of conditions. Some amino acids present in the binding and catalytic sites of these glucosidases, e.g., the peptide motifs TFNEP and ITENG, are also highly conserved and could help stabilize the protein core. In contrast to the tertiary structures, a diversity of β-glucosidase quaternary structures has been observed; oligomerization may lead to variable functions or allow specific regulation of bioactivation of defense compounds.

Enzyme stability might also be influenced by glycosylation. We detected two myrosinase-derived peptides (those containing N90 and N482) that contain residues proposed (Burmeister et al., 1997) to be modified by sugars. Their presence in our frass extracts coupled with the relatively broad migration of the parent protein in SDS–PAGE suggest that these proteins which are glycosylated may have been partially cleaved and deglycosylated during digestion. The actual contribution of these modifications to protein activity is still unclear, as previous studies have demonstrated that myrosinase glycosylation is not essential for activity (Zhou et al., 2012), while monocotyledonous glucosidases are not glycoproteins. However, glycosylation may contribute to stability with deglycosylated proteins being more susceptible to further attack by insect gut proteases, oxidants, electrophiles (including their enzymatic products), and to variation in temperature and pH.

Further features responsible for protein stability in insect guts have been revealed by investigations of maize ZmGlu1. The dimerization of this protein is stabilized by a disulfide bridge that shields a cluster of hydrophobic residues in the active site from the solvent (Rotrekl et al., 1999; Zouhar et al., 2001). ZmGlu1 also has a high number of proline residues as well as several intramolecular ion pairs and hydrogen and electrostatic bonds (Czjzek et al., 2001), which give the protein thermal stability and resistance to denaturing agents. The average number of hydrogen bonds per residue is rather high (>1) in ZmGlu1 and S. alba myrosinase MA1 (Czjzek et al., 2001).

Our proteomic analyses also support the assignment of more specific in vivo roles to individual β-glucosidase isoforms in spite of their overlapping substrate preferences in vitro. For example, ZmGlu1 but not ZmGlu2, and S. alba MA1 but not the other annotated S. alba myrosinases, were unambiguously identified in insect frass extracts. As most of the catalytic activity of the ingested foliage was still present in these frass extracts, these particular enzymes may therefore be the isoforms most relevant in activating their respective plant defensive substrates. However, this is in contrast to a previous proteomic analysis of frass from maize-fed S. frugiperda larvae, where ZmGlu2 was abundantly detected (Ray et al., 2016). Future studies are needed to clarify the comparative stabilities and substrate preferences of these close homologs and dissect their functions. The degraded isoforms may fulfill other roles, such as hydrolysis of cytokinins by ZmGlu enzymes (Brzobohatý et al., 1993). The almond enzyme, on the other hand, was found not to hydrolyze the cyanogenic diglucoside amygdalin in vitro, but may still act to hydrolyze prunasin (the corresponding monoglucoside), although this possibility was not examined and its activity was measured using the general glucosidase substrate pNPG.

The maize β-glucosidase aggregating factor (BGAF) whose aggregation of ZmGlu might protect the latter from insect proteases (Kittur et al., 2007) also resists the digestive system of S. littoralis. However, the roles of such aggregation factors have not been fully determined. Our recovery of maize BXD β-glucosidase activity did not differ between frass derived from feeding whole plant tissues (containing BGAF) and that from feeding semi-purified β-glucosidase preparations (presumably without BGAF), suggesting that BGAF did not protect against digestive inactivation. However, strongly interacting protein partners could have remained bound to ZmGlu during our simple purification procedure.

Insect gut pH can influence enzyme stability and also have a direct impact on catalytic activity. While larval Orthoptera, Hemiptera, and the larvae of most coleopteran families have slightly acidic to neutral midguts, many larvae of Lepidoptera, Diptera, and scarab beetles (Coleoptera) have highly alkaline midguts (Terra and Ferreira, 1994). In previous work, we observed that the midgut lumen of S. frugiperda fed on maize leaves is alkaline, but becomes neutral toward the hindgut (Wouters et al., 2014). While maize DIMBOA-Glc β-glucosidases were active under both neutral and basic conditions, their enzymatic activity was much higher at pH 7.0 (closer to their slightly acidic pH optima) than at pH 10.0 (Wouters et al., 2014). This resulted in slow hydrolysis of DIMBOA-Glc in the alkaline fore- and midgut, but extensive activation in the hindgut and neutral rectum where absorption of water and salts takes place. In the case of Zygaena filipendula, the cyanogenic β-glucosidases of Lotus spp. also had lower activity in the highly alkaline gut significantly reducing cyanogenesis (Pentzold et al., 2014). Therefore, by lowering the activity of these β-glucosidases, the high pH observed in the midgut (but not hindgut) lumen of these and other lepidopteran herbivores may help the insect to partially counteract the resistance of these defensive β-glucosidases to proteolytic inactivation in the gut.

Conclusion

Here, we have shown that selected plant β-glucosidases catalyzing the activation of chemical defenses survive digestion in a generalist lepidopteran herbivore, and are detectable in frass using in vitro activity assays and proteomic analyses. Although the structural features responsible for persistence in the digestive tract remain to be fully demonstrated, these enzymes seem particularly resistant because of their compact tertiary structure and the glycosylation of several residues, which have optimized them for resisting the protease-rich and high pH gut of lepidopteran herbivores. In parallel these features have allowed the use of such plant β-glucosidases in industrial applications even in the presence of detergents and organic solvents. The resistance to insect gut proteases even under high pH is presumably a result of natural selection on these enzymes to catalyze the continuous release of toxic aglycones in the herbivore gut, especially in the hindgut and rectum where the lumen reverts to more neutral conditions. Further investigations are needed to establish this extended role of β-glucosidases in the herbivore gut, and to determine which features of the proteins or their binding partners are responsible.

Author Contributions

DV and NW designed the experiments. DV, NW, AG, SG-J, and YH acquired and analyzed the data. DV, NW, AG, AS, and JG interpreted the data and drafted the article.

Conflict of Interest Statement

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.