1. Introduction
Brassicaceae have a high health potential thanks to their high amounts of several nutraceutical compounds such as glucosinolates and glutathione [
1]. The concentration of these and many other beneficial compounds depends on genetic and environmental factors such as climate, sun exposure and the concentration of nutrients, particularly sulphates, in the culture medium [
2,
3].
Glutathione, the γ-Glu-Cys-Gly tripeptide, is the most abundant thiol compound in plant and animal cells. It is involved in many fundamental cellular functions such as the homeostasis of cellular redox potential [
4], the ascorbate–glutathione cycle that is the main route of detoxification of ROS thanks to its key role in reactions involving enzymes and metabolites with redox properties, avoiding thus oxidative damages in plants [
5,
6] and the regulation of nuclear proteins and transcription factors by S-glutathionylation, affecting chromatin structure and condensation [
7]. Thiols also have a pivotal role in defence and protection from biotic and abiotic stress [
8,
9]. In particular, glutathione can function as a substrate for the synthesis of phytochelatins [
10,
11,
12] involved in the detoxification of heavy metals. It is the main compound involved in the transport and allocation of sulphur during plant growth [
13] and can be used as a source of cysteine in sulphur deprivation [
14]. In plants, glutathione synthesis occurs from cysteine, which is the first organic metabolite formed during sulphur assimilation, through a two-step pathway catalyzed by the γ-glutamylcysteine synthetase to form γ-glutamylcysteine from glutamate and cysteine and glutathione synthetase that add glicine to such compounds to form glutathione [
15].
Glutathione analysis, as well as that of other thiol metabolites, is problematic because the reduced form rapidly oxidizes to the disulphide one during sample preparation [
16], leading to an overestimation of oxidized form and determining a wrong evaluation of physiological cell status. The derivatization reactions, used for their determination, often form unstable products that may affect the analytical results [
17]. In addition, because such compounds may be involved in different reactions, even the HPLC analysis may be difficult [
18].
In the literature, several methods for sulphur compounds analysis in biological samples are described. The first techniques used, which were based on enzymatic [
19], fluorimetric [
20] and colorimetric [
21] assays, were inadequate due to their low sensitivity and low reproducibility.
Recently, several methods based on high-performance liquid chromatography (HPLC) were devised. These methods use: (i) UV absorbance measurement of sulphur compounds following derivatization with reagents that bind the amino group, such as 5,5 γ-dithiobis-(2-nitrobenzoic) acid (DTNB or Ellman’s reagent) [
22] or the sulfhydryl group such as 2,4-dinitrofluorobenzene [
23]; (ii) fluorescence measurement following derivatization with monobromobimane [
24], with
ortho-phthalaldehyde (OPA) derivatives [
25], 1-dimethylaminophthalene-5-sulfonyl chloride (dansyl chloride) [
26], 7-fluorobenzo-2-oxa-1,3-diazol-4-sulfonate (SBD-F) [
27]; and with or with derivatives of
N-(1-pyrenyl) maleimide (NPM) [
28]. There are also electrochemical methods that do not require derivatization [
29,
30].
Furthermore, many reagents used for analysis can support reactions with multiple functional groups or produce unstable derivatives, complicating the determination even via HPLC [
18]. Most of the methods in the literature have complicated protocols, time-wasting procedures, low reproducibility and often insufficient sensitivity for poorly concentrated samples. Finally, most of them still have difficulty measuring the reduced and oxidised thiol forms.
For these reasons, the goal of this study was to develop a method for the determination of thiols that would overcome these problems. We report an improved automated HPLC method for the determination of cysteine (Cys/Cis), γ-glutamylcysteine (γ-GC: oxidised and reduced) and glutathione (GSH/GSSG) based on derivatization with monobromobimane (MBB) in plant samples compared to that previously described by Noctor and Foyer [
31]. We have optimized the sample extraction and handling, the derivatization reaction, the chromatographic separation and the fluorimetric detection. It shows excellent reproducibility and high sensitivity. The method has been tested to determine simultaneously the concentration of cysteine, γ-glutamylcysteine and glutathione in plant extracts of broccoli (
Brassica rapa L. subsp. Sylvestris var. esculenta).
2. Materials and Methods
2.1. Plant Material and Experimental Design
Seeds of ecotype of Brassica rapa L. c.v. sylvestris (Blumen company, Piacenza, Italy) were soaked in 2 mM calcium sulphate (CaSO4) for thirty minutes to stimulate germination. Subsequently, seeds were placed on filter paper imbibed with deionized water in the dark at 22 ± 2 °C for four days. Seedling sprouts were transferred to hydroponics using plastic dishes (15 × 25 × 12 h cm) filled with distilled water. Growth was carried out under the following conditions: photoperiod 16 h, PAR 500 µE, 16/19 °C day/ night, 60–70% RH. After 10 days of cultivation, the young plants were separated into 4 groups with different sulphur nutrition (K2SO4: 0.01, 0.05, 0.5 and 2 mM) in modified Hoagland medium (5 mM KNO3, 0.5 mM KH2PO4, 0.5 mM K2HPO4, 1 mM CaCl2, 1 mM MgCl2, 0.05 mM FeCl3, trace elements solution: 10 mL/L). The nutrient solution was replaced every three days. After 45 days of hydroculture, leaves, petioles, inflorescences and roots were collected separately, frozen in liquid nitrogen and stored at −80 °C. Subsequently, plant material was powdered in a mortar and pestle with liquid nitrogen, and aliquots of 0.1 g of powdered material were stored at −80 °C for further analyses.
Experiments were performed in triplicate.
2.2. Chemicals and Reagents
GSH, GSSG, γ-GC, cysteine (Cys), cistine (cis), monobromobimane (MBB), 3-[4-(2-hydroxyethyl)-1- piperazinyl] propanesulfonic acid (EPPS), Diethylene triamine pentaacetic acid (DTPA), Dithiothreitol (DTT), hydrochloric acid (HCl) and sodium hydroxide were purchased from Sigma (St. Louis, MO, USA), acetic acid from Carlo Erba Reagent SPA (Milan, Italy), methanol and acetonitrile (HPLC grade) from Romil Ltd. (Cambridge, UK). All HPLC buffers were filtered through 0.45 mm Millipore (Billerica, MA, USA) filters prior to use.
2.3. Standard Solutions
The stock solutions of GSH, GSSG, Cys, cis and γ-GC were prepared in 0.1 N HCl to minimize oxidation and stored at −20 °C for 1 month. Working standard solutions were prepared daily by dilution in 0.1 N HCl.
2.4. Derivatization Reagent
MBB stock solution was prepared under dim lighting by dissolving 25 mg of reagent-grade MBB in 614 μL of acetonitrile (final concentration 150 mM) and stored at −20 °C for 1 month. Working 30 mM solutions were prepared daily by dilution in acetonitrile and used within 3 days.
2.5. Chromatographic Conditions
A Hewlett & Packard HPLC system (model 1100; Agilent Technologies Cernusco sul Naviglio MI, Italia S.p.A, Italy), equipped with an autosampler provided the online derivatization of samples and the binary pump, allowing the construction of a discontinuous gradient (
Table 1). The separation was performed using a Beckman Ultrasphere C-18 ODS (250 × 4.6 mm ID 5 μM particle size) protected by a C18 Security Guard pre-column (4 × 3 mm ID; Phenomenex Inc., Torrance, CA, USA). The column was kept at 25 °C. The fluorescence detection of MBB adducts was performed with excitation at 392 nm and emission at 480 nm. The control and data analyses were performed through the HP ChemStation software version A.06.03 Agilent Technologies LC/MSD Cernusco sul Naviglio MI, Italia S.p.A, Italy.
2.6. Extract Preparation and Pre-Column Derivatization and Analysis
2.6.1. Extract Preparation
Frozen powdered samples (0.1 g) were suspended in 1 mL of cold 0.1 M HCl degassed solution, mixed and kept at 4 °C for 2 h. Then, they were centrifuged at 13,000× g for about 2.5 min (Microcentrifuge, 5424R, Eppendorf, Hamburg, Germany), and the resultant supernatants were analyzed by HPLC.
2.6.2. Pre-Column Derivatization Procedure
Briefly, after transferring sample extracts in plugged vials in the HPLC array, the autosampler program allowed us to do the following procedure: (i) keep 30 μL of extract or standard solution by the needle from the sample vial in the HPLC array; (ii) wash the outside of the needle five times in a vial containing a water/methanol 50/50 solution; (iii) keep 70 μL of 200 mM EPPS, 5 mM DTPA buffer pH 11.7 (to have a final pH of 9.0); (iv) wash the outside of the needle; (v) keep 6 μL of 30 mM MBB solutions from a selected vial; (vi) eject the whole solution in a clean reaction vial; vii) mix it five times; (viii) wait for 7.5 min; (xi) keep 70 μL of the derivatized sample and inject it into the column for the analysis.
The MBB reacts only with the reduced thiols (free-SH). In order to also determine the oxidized form, an automated reduction step with DTT was included in the software program before MBB reaction (point v). Therefore, after step iv, the autosampler program allowed us to (v’) add 3.8 μL of 10 mM DTT (freshly daily prepared solution) to the solution vial; (vi’) mix the whole solution five times; (vii’) wait 2 min; (viii’) wash the outside of the needle as above, in order to not contaminate the reagent solution; (xi’) keep 6 μL of MBB 30 mM solution; (x’) add it to the same reaction vial; (xi’) mix five times for the derivatization procedure; (xii’) wait for 7.5 min; (xiii’) keep 70 μL of the derivatized sample and inject it into the column for the analysis.
The software program repeated the procedures of derivatization, separation and analysis for each sample extract. For each sample, a new clean reaction vial for each analysis was programmed in the array, and the selected vials containing reagents and the washing solution were substituted and refilled daily. In addition, to prevent false results due to the eventuality that the needle could be partially clogged by some sample particles that remained in the decanted extract, each analysis was repeated three times.
2.6.3. Chromatographic Conditions and Analysis
The separation and analysis were performed by using the cromatograpic apparatus described in
Section 2.5 (chromatographic conditions). The quantification was performed under the binary elution gradient elution (mobile phase: A = 40 mM sodium acetate and methanol 17% pH 3.9; B = methanol (100%)), with the gradient reported in
Table 1. During the derivatization procedure, the column was furtherly balanced with the eluent A.
The fluorescence detection of MBB adducts was performed with excitation at 392 nm and emission at 480 nm. The control and data analyses were performed through the HP Chemstation software.
2.7. Validation by Recovery Experiments
Authentic GSH, GSSG, Cys and CSSC were added before analysis to the different tissue samples from plants grown with different sulphur nutrition (K2SO4: 0.01, 0.05, 0.5 and 2 mM) or from different plant parts. The added internal standards did not mask endogenous levels because their concentrations were about the same recorded in the samples. Control samples without internal standard addition were also run at the same time. The analytes were analyzed alone and together to verify that there was not any interference in order to validate the method.
2.8. Statistical Methods
Results for oxidized and reduced content of Cys, γ-GC and GSH were expressed as mmol/100 g fresh weight. Data were elaborated and analysed using SigmaPlot v. 14.5, Systat Software Inc., San Jose, CA, USA.
3. Results
3.1. Samples Preparation
An accurate determination of GSH and GSSG levels depends on prevention of GSH oxidation during sample preparation. To reduce this phenomenon, the contact with oxidizing agents was minimized. In order to remove oxygen from the solutions, nitrogen was bubbled into the extraction and assay mixes for 10 min prior to use and nitrogen was flushed onto the samples during their preparation. DTPA was used as a chelant for heavy metal ions. Finally, the oxidation of the thiols during extraction was minimized by using 0.1 N HCl as extracting solution.
3.2. Optimization of Fluorometric Detection
The optimal excitation and detection wavelengths of MBB derivatives were determined by performing a fluorescence spectrum in which the excitation wavelength was set in the range 300–430 nm for excitation and 440–600 nm for emission detection (
Figure 1). The maximum fluorescence value was obtained at 392 nm for the excitation and 480 nm for the emission (
Figure 1).
In this view, the excitation for further analysis was set at 392 nm and emission at 480 nm. These conditions were optimal also for cysteine and γ glutamylcysteine adducts.
3.3. Chromatographic Gradient Optimization
The chromatogram (
Figure 2) was obtained using the chromatographic gradient described in
Table 1.
The gradient was initially very slow (between 2 and 28 min) to resolve the peaks of the tested compounds. Under such conditions, no interference of reagent or adduct peaks were evidenced even at low sample concentrations, which differed from previously reported data [
32]. The elution times for Cys, γ-GC, and GSH were 5.95, 6.96 and 8.64 min, respectively. The identification was made by direct comparison with the pure standards as well as in biological samples in which the internal standards were added.
After the elution of compounds of interest, the gradient was kept at 100% methanol for 5 min to wash the column and then rebalanced to the initial elution buffer gradient with 40 mM sodium acetate in 10% methanol for 10 min. Under these conditions, no significant displacement of the elution times in subsequent analyzes occurred. After analysis, the column was washed with 100% methanol for 30 min.
3.4. Monobromobimane Solution Stability
In order to analyse many samples loaded by the autosampler HPLC system at room temperature, the MBB solution stability over time was tested.
Figure 3 shows the result of the derivatization of a fresh solution of GSH 20 µM by a solution of MBB 30 mM in acetonitrile kept at 25 °C in the dark for up to eight days. MBB solution still retained over 90% reactivity towards GSH after 3 days of storage.
3.5. Derivatization Reaction in Dependence on pH
As shown in
Figure 4, the MBB derivatization reaction occurred at basic pH, and the reaction rate increased up to pH 11.00. However, at higher pH, the MBB reactivity with other functional groups different from thiols increased, leading to the appearance of other peaks that complicated the analyses. In this context, the best condition for a clear chromatogram was to stay at pH 9.00, at which the reaction rate was still 90% of the maximum.
3.6. Reaction Time Course
The long reaction time was responsible for increasing the number of interfering peaks due to the reaction of the MBB with other functional groups or with di-sulphides, tri-sulphides and other compounds generated by multiple functional group thiols.
The reaction time indicated in the literature is 10–15 min [
16,
31].
Figure 5 shows tests carried out by reacting a solution of GSH 20 µM with MBB at different times.
The best compromise to enhance chromatogram quality, thiol selectivity and limit sensitivity was 7.5 min. In these conditions, in fact, more than 98% of GSH reacted, resulting in a clear chromatogram more than that produced in longer reaction times (not shown). In this time, cys and γ-GC also largely reacted.
3.7. Data Analysis
To verify the linearity of the fluorescence response in the overall analysis, serial dilutions of pure standards of GSH, Cys, and γ-GC were injected in sequence.
Figure 6 shows the calibration curves for Cys, γ-GC and GSH. The linear correlation coefficient R
2 = 0.9980, R
2 = 0.9981 and R
2 = 0.9984, for Cys, γ-GC and GSH respectively, ensured a linear response from 0 to 1.2 nmol injected. The high linear correlation coefficient R
2 (greater than 0.9980) provided a linear response at least up to 1.2 nmol of injected adduct (
Figure 6). The fluorescence coefficient for pmoles of injected compounds was 3.2760 for Cys, 3.1685 for γ-GC and 3.7316 for GSH. Method reproducibility was verified by injecting 5 times the same leaf extract of
Brassica rapa L. derivatized with MBB. The per cent coefficient of variation was 1.1 for Cys, 1.6 for γ-GC and 1.4 for GSH.
The recovery of GSH, Cys and γ-GC from Brassica rapa L. tissues was tested by adding an internal standard at concentrations similar to the endogenous ones in the 0.1 M HCl solution used to extract the thiols from the samples. The analyses were performed in triplicate by comparing samples with and without internal standard. Estimated recoveries were 97.4 ± 1.4% for GSH, 96.7 ± 1.6% for γ-GC and 97.6 ± 1.2% for cysteine.
The lower detection limit was determined for GSH by carrying out serial dilutions of GSH. It was evidenced that a signal-to-noise ratio higher than 3 could still be measured when 100 femtomol of GSH in 30 μL were sampled for the derivatization (corresponding to 3 nM).
3.8. Determination of Thiols in Brassica rapa L. Subsp. Sylvestris Var. Esculenta
The method was tested to determine the oxidized and reduced form of Cys, γ-GC and GSH content in florets, stems and leaf blades of
Brassica rapa L. subsp. Sylvestris var. esculenta grown in hydroponics in Hoagland medium supplemented with sufficient sulphate (2 mM) (
Figure 7;
Table 2).
The average contents of GSH/2GSSG were around 840/45, 345/70 and 150/70 neq g
−1 FW for the inflorescences, laminae and petiole, respectively; those of Cys/2Cis were 80/12, 29/12 and 24/6 neq g
−1 FW; and those of γ-GC/2γ-GCCG-γ were 8.0/4.0, 6.0/3.0 and 3.0/2.0 nmol g
−1 FW, respectively. Similar values have been found for Brassicaceae and other plants [
33,
34,
35,
36,
37]. Glutathione showed the highest concentration in the florets, being two- and three-fold higher than that found in the leaf blades and in the stems, respectively. Even though the Cys and γ-GC showed a similar pattern of distribution among the plant parts, this occurred at concentrations much lower than that of glutathione for Cys and γ-GC.
The method was also applied to evaluate total and reduced contents of glutathione in inflorescences, laminae, petioles and roots of
Brassica rapa L. cv Sylvestris plants grown in nutrient solutions containing 0.01, 0.05, 0.5 and 2 mM sulphate. Higher glutathione amounts were recorded in the inflorescences with a content two- and three-fold higher than that found in leaves and petioles, respectively. Total cysteine concentration (cys+cis) in broccoli plants (
Table 3) was about 20-fold lower than the average glutathione content; however, like glutathione, the total cysteine content was also affected by sulphate availability. Cysteine concentration in the inflorescences grown under high sulphur nutrition was tripled compared to that recorded in the inflorescences grown with low sulphur medium and increased ten-fold compared to that recorded in the inflorescences grown under sulphate deficiency. In leaves, petioles and roots, the total cysteine content also increased in response to sulphate, but only up to the 0.5 mM sulphate concentration (
Table 3).
The fraction of reduced glutathione was about 90% in the inflorescences, and it was not influenced by sulphate availability (
Table 3). Conversely, in the other analysed tissues, the glutathione-reduced fraction was dependent on sulphate; in fact, it was reduced to 75% of that recorded in leaves of plants grown with sulphate deficiency. The cysteine reduced fraction (
Table 3) followed the same trend as that observed for glutathione, with lower values recorded in roots and leaves grown under sulphur deficiency (50% and 60% of reduced cysteine, respectively). As for glutathione, reduced cysteine content was not influenced by sulphate availability. The maximum concentration of reduced cysteine (90% of total cystine content) was recorded in the plants grown under high sulphur availability.
4. Discussion
The determination and quantification of glutathione and other thiols in biological samples require a very sensitive, fast and selective method. Their concentrations, in fact, in particular those of glutathione precursors, Cys and γ-GC, are generally low in plant samples. Furthermore, due to their extreme sensitivity to oxidation, the concentration ratio between the reduced and oxidized form (GSH/GSSG), which provides important information on the cellular redox state, can easily be altered during sample handling.
When MBB is used for thiols derivatization, being a highly reactive compound, it is important to avoid the formation of different products that could complicate the chromatographic analysis and lead to an incorrect determination of the tested compounds. The method for the quantification of thiols described in this paper allows obtaining an efficient chromatographic separation in HPLC, accompanied by a high sensitivity guaranteed by the high fluorescent thiol-bimane adducts. Furthermore, having demonstrated the stability of the MBB solution in acetonitrile even at room temperature (
Figure 3), a fully automated sample handling method could be designed, allowing many samples to be run without changing the reagent solution. Another optimization consisted in the grinding of the tissues in liquid nitrogen and the extraction carried out in 0.1 N HCl, which strongly limited the oxidation of the samples. After 2 days of storage at room temperature (data not shown), the samples extracted in 0.1 N HCl still retained 90% of reduced glutathione. Furthermore, the method provided for the immediate injection of the reaction mixture, avoiding the need to add acid to stop the reaction [
31]. On the other hand, the stability of thiol-bimane adducts has been previously demonstrated [
33]. The optimization of the elution gradient in HPLC made it possible to obtain an excellent peaks resolution maintaining a good speed of analysis. The gradient applied (
Table 1) allowed us to elute cysteine, γ-GC and GSH far from the main interfering peaks due to the presence of the reagent, itself fluorescent. The latter were removed from the column with 100% methanol.
The derivatization reaction carried out at pH 9 (
Figure 4), a value higher than that described in the literature [
16,
31], and the reaction time reduced to 7.5 min (
Figure 5), which allowed the reaction of over 97% of the glutathione, rendered clean the chromatogram with a signal-to-noise ratio of 3.3 when 100 fmol of glutathione were injected. In conclusion, the method developed showed an effective separation of the peaks of interest, a high sensitivity (up to 100 fmol), a very good repeatability (maximum 1.6% error), a linear correlation coefficient of the straight calibration (R
2 > 0.99) (
Figure 6) and a high internal standard recovery (>96%), allowing the analysis of about 50 samples per day in total automation without any intervention by the operator.
We tested our optimized method to determine the endogenous content of Cys, γ-GC and GSH in plants under enriched sulphur nutrition. The method was suitable for this evaluation, and our results showed that the S deficiency during growth led to an increase in sulphate in the inflorescences compared to leaf blades. In plants grown under sulphur deficiency and low sulphur availability, in fact, sulphate accumulated mainly in the inflorescences, conversely to that observed in plants grown in medium under high sulphur nutrition (
Table 3). These results suggest that plants are able to redistribute the available sulphate to tissues that have the major need in order to guarantee plant survival and reproduction.
Notwithstanding this representing an optimization of the well-known MBB method [
16] and the fact that other attempts have been done, also recently, to offer improvements of this method by fluorescence detection [
34,
35], this is the first time it has been optimized for a complex tissue like that of
Brassica rapa, in which quenching caused by interference with components present in the plant extract makes it difficult to perform a fluorescent detection.