Next Article in Journal
Integrated Transcriptomics and Metabolomics Analysis Reveals the Effects of Cutting on the Synthesis of Flavonoids and Saponins in Chinese Herbal Medicine Astragalus mongholious
Previous Article in Journal
New Implications of Metabolites and Free Fatty Acids in Quality Control of Crossbred Wagyu Beef during Wet Aging Cold Storage
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Advancing Glucose Conjugated Gibberellins Discovery: A Structure–Oriented Screening and Identification Method for Unraveling Gibberellin Metabolites in Plants

1
School of Bioengineering and Health, Wuhan Textile University, Wuhan 430200, China
2
Department of Chemistry, Wuhan University, Wuhan 430072, China
3
School of Public Health, Wuhan University, Wuhan 430071, China
4
Frontier Science Center for Immunology and Metabolism, Wuhan University, Wuhan 430071, China
*
Authors to whom correspondence should be addressed.
Metabolites 2024, 14(2), 96; https://doi.org/10.3390/metabo14020096
Submission received: 23 December 2023 / Revised: 16 January 2024 / Accepted: 20 January 2024 / Published: 29 January 2024
(This article belongs to the Section Plant Metabolism)

Abstract

:
Gibberellins (GAs) play a pivotal role in modulating plant growth and development. Glucose–conjugated gibberellins (Glc–GAs), a prevalent conjugated form of GAs, regulate intracellular GA levels by the coupling and decoupling of glucose groups. However, the diversity of Glc–GAs identified within individual species remains limited, hinting at a multitude of yet undiscovered gibberellin metabolites. This lacuna poses considerable impediments to research efforts dedicated to comprehensively delineating the GA metabolic pathway. In this study, we developed a structure–oriented screening and identification method for Glc–GAs in plant species by employing LC–MS/MS coupled with chemical derivatization. Through the application of chemical derivatization technique, carboxyl groups on Glc–GAs were labeled which effectively enhanced the sensitivity and selectivity of mass spectrometry detection for these compounds. Concurrently, the integration of mass spectrometry fragmentation and chromatographic retention behavior facilitated the efficient screening and identification of potential Glc–GAs. With this strategy, we screened and identified 12 potential Glc–GAs from six plant species. These findings expand the Glc–GA diversity in plants and contribute to understanding GA metabolic pathways.

1. Introduction

Gibberellins (GAs) are carboxyl–containing plant hormones characterized by a tetracyclic diterpene skeleton. They are closely involved in various plant growth and development processes, including seed germination, hypocotyl elongation, and control of flowering time [1,2]. During the physiological regulation of GAs, free GAs can reversibly bind to small molecular groups through processes like glycosylation, acylation, or esterification, forming conjugated GAs [3,4]. Conjugated GAs also play important roles, mainly by altering the levels of GAs in the plant and regulating the endogenous GA reservoir through the binding and unbinding of conjugated groups [5,6]. During seed germination or plant growth periods, free GAs can be released from conjugated GAs through enzymatic hydrolysis to exert their effects, whereas, in mature plants, free GAs are converted into conjugated GAs and stored, which may be a crucial step in further degradation and metabolism processes [7,8,9]. Moreover, considering the polar properties of some conjugated GAs, they are also believed to be involved in intracellular compartmentalized transport and long–distance transport throughout the plant [10,11]. Therefore, conjugated GAs commonly serve as storage compounds, facilitators of transport, and intermediates in breakdown metabolism, playing an important role in the regulation of GA homeostasis.
Glucose derivatives are the most common conjugated GAs in plants [12]. Glycosylation is a component of both primary and secondary metabolism in plant tissues, and plants can maintain cellular homeostasis by regulating the level, biological activity, and subcellular distribution of glycosylated compounds [3,4,13]. It has been suggested that the polar glucose moiety in some conjugated GAs may cause a twisted orientation of the gibberellin molecule in the membrane, thereby preventing its binding to specific receptors [13,14]. Since the first gibberellin glucoside, GA8–2–O–glucoside, was isolated from mature beans, researchers have identified a series of GA glucose derivatives. Glucose–conjugated gibberellins (Glc–GAs) are mainly divided into two categories: one is the connection of glucose to the hydroxyl group on the GA skeleton, called GA glucosides, in which the glucose component can be connected to the carbon 2–O, 3–O, 11–O, 13–O, or 17–O position of the parent GA, generating a series of isomers. The other is the connection of glucose to the carboxyl group at the C–6 position of the GA skeleton, called GA glucosyl esters [7]. GA glucosyl esters are widely distributed in legumes, and their contents are relatively high in mature seeds [15]. In plants, free GAs are synthesized into gibberellin glucosides and glucosyl esters by glycosyltransferases to reduce the level of free GAs. Glycosyltransferase from mature bean fruit tissue can specifically glycosylate GA3, forming GA3–3–O–glucoside using UDP–glucose as the glucose donor [15]. Piccoli et al. found rapid conversion among GA glucosides, GA glucosyl esters, and free GAs in maize seedlings [16]. Thus, the presence of conjugated gibberellins can effectively maintain GA homeostasis and prevent plants from synthesizing free GAs de novo, which plays a crucial role in the dynamic balance of plant hormones.
As derivatives of free GAs, Glc–GAs typically exist in trace amounts in plant tissues [17]. Analyzing these compounds has been challenging due to their lack of luminescent or chromophoric groups and the complex plant matrix in which they reside [3]. Research on Glc–GAs has mainly focused on the period from the 1960s to the 1990s, during which researchers usually employed activity–directed screening strategies combined with GC–MS analysis to discover and identify Glc–GAs in plant samples [18,19]. Typically, a large amount of plant material is extracted with organic solvents, followed by liquid–liquid extraction for preliminary purification. The resulting fractions obtained from HPLC fractionation were subjected to activity assays, and the active fractions were subsequently analyzed by GC–MS and compared with synthetic standards to obtain qualitative results [18,20,21,22]. This approach is time–consuming, requires large quantities of material, and may overlook the majority of inactive or low–activity Glc–GAs. As a result, the variety of Glc–GAs found in plants is very limited. Compared to the over one hundred species of free–form GAs currently discovered in various organisms, the identified types of conjugated GAs are extremely limited. According to our literature survey, 20 species of conjugated GAs have been reported (Table S1). However, there is still a lack of efficient screening and identification methods for low–content and inactive or low–activity Glc–GAs [23]. Additionally, there is still little knowledge about the mechanism of action of glycosyltransferases in biology, and the relevant genes expressing glycosyltransferases have not been identified. Therefore, achieving accurate and quantitative analysis of trace Glc–GAs will contribute to the discovery of new Glc–GAs in plants and the study of their regulatory mechanism.
In previous studies, we developed a chemically assisted LC–MS technique for GA screening, which had the advantages of strong structure orientation, high sensitivity, and the ability to detect more low–activity GAs [24]. Building upon this, we established a method that combines this technique with structure–oriented Glc–GA screening and identification. By using chemical derivatization to enhance the mass spectrometric detection sensitivity of Glc–GAs, we were able to screen for potential Glc–GAs at low concentrations. Furthermore, the labeled Glc–GA molecules exhibited characteristic fragmentation patterns that could be used for structural inferences, assisting in the identification of potential Glc–GA structures. Ultimately, we screened and identified 12 potential Glc–GAs in six plant samples. This work is of great significance in improving our understanding of GA metabolic pathways and investigating the mechanisms of GA function.

2. Materials and Methods

2.1. Chemicals and Materials

Gibberellin standard GA3 was purchased from Olchemim Co., Ltd. (Olomouc, Czech Republic). Ammonia solution (25%), hydrochloric acid (HCl), and sodium chloride (NaCl) were obtained from Sinopharm Chemical Reagent Co., Ltd. (Shanghai, China). α–D–bromotetraacetylglucopyranose (98%), silver oxide (Ag2O, 99.7%), and sodium methoxide (NaOMe, 97%) were purchased from Aladdin Reagent Co., Ltd. (Shanghai, China). 4A molecular sieves were obtained from Macklin Reagent Co., Ltd. (Shanghai, China). Chromatography–grade methanol (MeOH) was purchased from Merck Reagent Co., Ltd. (Darmstadt, Germany). Dichloromethane (99.9%, anhydrous solvent with molecular sieves) was acquired from J&K Scientific Co., Ltd. (Beijing, China). Analytical–grade formic acid (FA, 88%), 2–chloro–1–methylpyridinium iodide (CMPI), triethylamine (TEA), and N,N–dimethylethylenediamine (DMED) were purchased from Sinopharm Chemical Reagent Co., Ltd. (Shanghai, China). The isotopically labeled reagent d4–DMED was synthesized following the reported methods in our laboratory [25]. Chromatography–grade acetonitrile (ACN) was obtained from Merck Reagent Co., Ltd. (Darmstadt, Germany). SAX–SPE (strong anion exchange–solid–phase extraction column, 3 mL, 200 mg) was purchased from Weltech Co., Ltd. (Wuhan, China). Milli–Q ultrapure water Co., Ltd. (Millipore, Bedford, USA) was used throughout the experiments. Stock solutions of TEA (10 μmol/mL), CMPI (20 μmol/mL), DMED (20 μmol/mL), d4–DMED (20 μmol/mL), and organic acids (1 mg/mL) were prepared in chromatography–grade ACN and stored at −20 °C in a refrigerator.

2.2. Collection of Plant Samples

Mature seeds of Canavalia gladiata, Phaseolus vulgaris, Pisum sativum, and Vicia faba were purchased from local supermarkets in Wuhan, China. Vigna angularis was provided by the Vegetable Breeding Center of the Hubei Academy of Agricultural Sciences (Wuhan, China). Rice (Oryza sativa ssp. japonica cv. Nipponbare) samples were collected during the grain–filling stage, frozen in liquid nitrogen, and stored at −80 °C.

2.3. Extraction and Derivatization of Glc–GAs in Plant Samples

Extraction of Glc–GAs from plant samples: 2 g of plant sample powder was weighed and mixed with 20 mL of MeOH (1 mL/100 mg) as the extraction solvent. The mixture was vortexed for 5 min and left to extract at 4 °C overnight. After centrifugation (5000 rpm, 5 min, 4 °C), the supernatant was collected and purified using SAX–SPE [26]. The purification process involved three steps: activation of the SPE column with 3 mL of MeOH, sample loading, and washing with 3 mL each of MeOH/H2O (1/9, v/v) and H2O/MeOH (1/9, v/v) solutions, followed by elution with 3 mL of the elution solution (MeOH containing 1% FA). The eluted fractions were collected and dried under nitrogen flow with heating at 40 °C.
The purified samples were reconstituted with 200 µL of ACN, divided into two portions, and each portion was separately added with 20 µL of TEA (10 μmol/mL) and CMPI (20 μmol/mL), followed by vortexing. To one portion, 20 µL of DMED (20 μmol/mL) was added, and to the other portion, 20 µL of d4–DMED (20 μmol/mL) was added. The derivatization reactions were conducted at 40 °C with shaking at 1500 rpm for 60 min. After completion of the derivatization, the samples were dried under nitrogen flow with heating at 40 °C. The light– and heavy–labeled samples were reconstituted with 100 µL of ACN/H2O (1/9, v/v) each and combined. The injection volume for subsequent LC–MS analysis was 10 µL.

2.4. Synthesis of Glc–GA3

First, 2 mmol of GA3 (693 mg) standard was weighed, and the GA3 was reacted with 3 mmol of derivatization reagent (TEA, CMPI, DMED) in a molar ratio of reactant to product of 1:1.5 under oscillation at 40 °C for 1 h. The reaction products were purified using a manually packed reverse–phase C18 column to remove excess derivatization reagents. The products were then placed in a 50 mL round–bottom flask and mixed with 1 g of α–D–bromotetraacetylglucopyranose, 1 g of Ag2O, 20 pellets of 4A molecular sieves (previously dried), and 20 mL of anhydrous dichloromethane. The mixture was magnetically stirred at 25 °C for 24 h. The reaction products were extracted with saturated NaCl solution, and the organic phase was dried and dissolved in 10% MeOH prior to mass spectrometry analysis. The software Chemdraw Professional 17.1 was used to predict the hydrophobicity of Glc–GA3 isomers.

2.5. Instrumentation and Analytical Conditions

The instrument used for sample analysis was a UHPLC–MS system. The mass spectrometer employed was an Orbitrap Fusion Tribrid mass spectrometer (Thermo Fisher Scientific, Waltham, MA, USA). The liquid phase system consisted of an UltiMate 3000 UHPLC System (Thermo Fisher Scientific, USA). The chromatographic column used was a Waters ACQUITY UPLC BEH C18 Column (100 mm × 2.1 mm i.d., 1.7 μm, Waters, Milford, MA, USA), maintained at a column temperature of 40 °C. A 0.1% FA aqueous solution was used as mobile phase A and ACN as mobile phase B. The gradient elution for injection chromatography was as follows: 0–2 min, 5% B; 2–10 min, 5–15% B; 10–25 min, 15–50% B; 25–27 min, 50–90% B; 27–32 min, 90% B; 32–35 min, 90–5% B; 35–40 min, 5% B. The flow rate was set at 0.4 mL/min.
Data–dependent acquisition mode (DDA) was used to obtain mass spectrometry information. The DDA–HCD MS/MS data acquisition method employed the following parameters: FTMS (Full–MS): performed in positive ion mode, with a scan range of m/z 300–700; a resolution of 120,000; AGC target set to 2 × 105; precursor priority set to highest charge state and then most intense; FTMS2 (HCD): collision energy (%) at 50–60; a resolution of 30,000; ion source parameters: sheath gas flow at 35 psi, auxiliary gas at 15 psi, spray voltage at 3.5 kV, ion transfer capillary temperature at 320 °C, and heated block temperature at 350 °C.
The synthetic Glc–GA3 was detected in negative mode. The substances derived from DMED were detected in negative mode.

3. Results and Discussion

3.1. Strategies for Screening and Identifying Glc–GAs

Based on the structure–guided GA screening strategy previously reported [24], we established a screening method for Glc–GAs. This strategy utilized the mass spectrometry fragmentation behavior conferred by chemical derivatization to characterize Glc–GAs and screen for potential Glc–GA compounds. It combined mass spectrometry fragmentation and chromatographic retention behaviors for qualitative analysis and aided in structural inference. The specific workflow is illustrated in Figure 1. The plant samples were first pretreated for purification and enrichment. They were then reacted with a pair of derivatizing reagents (DMED/d4–DMED) according to the derivatization reaction equation shown in Figure 2. The labeled product was mixed in a 1:1 ratio and analyzed by UHPLC–MS to obtain the initial data. Table S2 lists the m/z values of labeled Glc–GA products for subsequent peak extraction and peak–pair matching. Additionally, the labeled Glc–GA products exhibited characteristic mass spectrometry fragmentation behavior. We employed five rules to screen the mass spectrometry spectra and identify potential Glc–GA candidates. Furthermore, we utilized mass fragmentation pathway analysis and quantitative structure–retention relationship models for structure elucidation, resulting in the final resolution of the structures.

3.2. Analysis of Mass Fragmentation Behavior of Glc–GAs–Labeled Products

Known GAs’ mass fragmentation behavior provides important clues for structure–oriented screening strategies. Our previous work has studied the fragmentation pathways of DMED–labeled products of existing GA standards, summarizing four distinctive fragmentation behaviors: (1) DMED–related fragmentation behavior; (2) GAs–A ring–related fragmentation behavior; (3) GAs–C–D ring–related fragmentation behavior; and (4) carbon skeleton–related fragmentation behavior [24]. These four characteristic fragmentation behaviors serve as important references for the analysis of the mass fragmentation behavior of Glc–GA labeled products.
Currently, Glc–GAs are not commercially available. In order to obtain Glc–GAs reference standards to explore the fragmentation patterns and retention time patterns using mass spectrometry, we attempted to synthesize Glc–GA3 following the methods described in references [10,27]. The structures of the substrate and potential products are shown in Figure S1. There are three possible Glc–GA3 main products: GA3–3–O–Glucoside, GA3–13–O–Glucoside, and GA3 glucosyl ester. Figure S2 shows the total ion chromatogram and extracted ion chromatograms (XIC) with m/z 507.1872 of the synthesized reaction products. From the XIC, it can be observed that there are three target product peaks with retention times of 8.78 min, 9.40 min, and 10.17 min, respectively. As shown in Figure S3, when observing the mass spectra at the retention times corresponding to the three product peaks, it is found that they exhibit a high degree of similarity in terms of fragment ion types. Among them, Figure S3A has noticeable differences in the relative intensities of the fragment ions compared to Figure S3B,C.
Furthermore, by comparing the hydrophobicity constants predicted by the software for three Glc–GA3 compounds, namely GA3–3–O–Glucoside, GA3–13–O–Glucoside, and GA3 glucosyl ester, it was observed that their hydrophobicity decreases sequentially. As a result, their retention on the chromatographic column decreases, leading to shorter peak elution times. Based on the results of the secondary mass spectrometry (Figure S3), we speculated that the peak at 8.78 min corresponded to GA3 glucosyl ester, the peak at 9.40 min corresponds to GA3–13–O–Glucoside, and the peak at 10.17 min corresponds to GA3–3–O–Glucoside.
DMED can react with the carboxyl group of Glc–GA3. Taking GA3–3–O–glucoside as an example, Figure 2 shows the corresponding reaction equation. Figure 3A displays the total ion chromatogram (TIC) and extracted ion chromatograms (XIC) of the synthetic products after the chemical labeling reaction. In the XIC, the blue peak represents the unreacted GA3 labeled product at m/z 417.2384, while the red peak represents the Glc–GA3 labeled product at m/z 579.2912. Due to the carboxyl site of GA3 glucosyl ester being occupied by the glucose moiety and not reacting with DMED, there are only two red peaks with retention times of 5.57 min and 6.25 min. Considering their hydrophobicity, we speculated that the peak at 5.57 min corresponded to GA3–13–O–Glucoside–DMED and the peak at 6.25 min corresponded to GA3–3–O–Glucoside–DMED. By observing the mass spectra at the retention times corresponding to these peaks, it is found that, apart from the obvious Glucose–related fragmentation behaviors, the Glc–GA3–labeled products exhibit highly similar changes in fragment ion types and relative intensities compared to the substrate GA3 labeled products (Figure 3B–D). Figure 4 shows possible fragmentation pathways for DMED–GA3–3–O–glucoside. In addition to the four characteristic fragmentation behaviors mentioned above, Glc–GAs also exhibit Glucose–related fragmentation behaviors for losing glucose (162.0528), such as the protonated molecule m/z 579.2912 of GA3–3–O–glucoside fragmenting into the product ion m/z 417.2384. This fifth characteristic fragmentation behavior provides a further basis for inferring the structures of Glc–GAs.

3.3. Establishment of Quantitative Structure–Retention Relationship (QSRR) Models

The Quantitative Structure–Retention Relationship (QSRR) model has gained significant attention in recent years. It associates the chromatographic retention properties of substances with their molecular structural features, enabling the prediction of compound retention behaviors and aiding in metabolite identification. In our previous research, we utilized a multivariate linear regression method to construct a QSRR model using a dataset of 13 GA standards. Through various cross–validation analyses, the model demonstrated good predictive ability and successfully predicted the retention index (RI) for 603 GA structures, with a root mean square value of 69.5832, which represents the threshold for matching experimental and predicted RIs [24].
In this research, we utilized the established QSRR model to assist in predicting the RI of potential Glc–GAs in our study. Glc–GAn exhibits a shorter retention time and lower RI compared to the corresponding GAn due to weaker glucose group retention on the C18 reversed–phase column. Figure 3A shows that the retention time of Glc–GA3 is earlier than that of GA3, and the RI of Glc–GA3 (575.25, 592.50) is smaller than that of GA3 (RI = 653.94). Therefore, considering the error threshold of the QSRR, the experimental RI value of potential Glc–GAn should be below the sum of the predicted RI value of GAn plus 69.5832 as shown in Equation (1), and we use this inference to facilitate the identification of potential Glc–GAn.
RIGlcGAn < Predicted RIGAsn + 69.5832

3.4. Screening Potential Glc–GAs in Plant Samples

After overnight extraction of plant samples at a low temperature, they were purified using strong anion exchange solid–phase extraction columns. The samples were then divided into two portions and reacted separately with the derivatizing reagents DMED and d4–DMED. The labeled products were mixed in a 1:1 ratio and subjected to UHPLC–MS analysis to obtain initial data. Peak extraction was performed using the simulated m/z values in Table S2, followed by peak matching based on two criteria: (1) consistent retention time (RTDMED–labeled–RTd4–DMED–labeled: 0–0.013 min) and (2) consistent peak intensity (IntensityDMED–labeled/Intensityd4–DMED–labeled: 0.76–1.50). Peak pairs that satisfied both criteria were considered potential Glc–GA compounds.
Taking sword bean samples as an example, Figure 5A shows the TIC and XIC plots of the labeled products in HRMS full scan mode, while Figure 5B provides information on the precursor ion peaks of the light and heavy labeled products. The peaks with m/z values of 563.2963 and 567.3214 exhibited consistent intensities and retention times, indicating that these compounds are potential Glc–GA compounds. As shown in Table S3, a total of 93 Glc–GA metabolites were detected in the six plant samples. Subsequently, we filtered the secondary mass spectra of 93 potential Glc–GAs using the aforementioned five categories of fragmentation patterns: (1) DMED–related fragmentation behaviors, such as the loss of (CH3)2NH (45.0578), C4H10N+ (72.0808), C3H6ON+ (72.0444), etc.; (2) GAs–A ring–related fragmentation behaviors, such as the loss of H2O (18.0106), HCOOH (46.1055), etc.; (3) GAs–C–D ring–related fragmentation behaviors, involving rearrangement via cleavage of C9–C11 and C8–C14 bonds, and loss of C2H4 (28.0313), etc.; (4) carbon backbone–related fragmentation behaviors, resulting in odd–numbered CxHy+ fragments in the m/z 80–230 range; (5) Glucose–related fragmentation behaviors, including loss of glucose (162.0528). Ultimately, we identified 12 Glc–GA candidates through this screening process. Their molecular weight information is provided in Table 1, and their MS/MS spectra are shown in Figures S4–S5, S9–S11, S13–S16, S18, S21, and S23.

3.5. Identification of Potential Glc–GAs in Plant Samples

After obtaining the 12 Glc–GA candidates, we performed structure identification through a four–step process. The first step involved determining the Glc–GA types corresponding to specific m/z values based on Table S2. In the second step, we analyzed the fragmentation pathways and filtered out Glc–GA candidates with inconsistent structures. The third step involved using a QSRR model to further eliminate Glc–GA candidates with mismatched retention indices (RI). In the fourth step, we compared the MS/MS spectra of the Glc–GA candidates with available standard GA compounds to obtain potential Glc–GA structures. Subsequently, we conducted an identification analysis on the screened potential Glc–GA compounds in the plant samples.
For compound 1 and compound 2, the Glc–GA candidates to which the labeled product with m/z 563.2963 corresponds include Glc–GA5, Glc–GA31, Glc–GA95, Glc–GA96, Glc–GA108, Glc–GA109, Glc–GA117, Glc–GA121, Glc–GA122, Glc–GA7, Glc–GA11, Glc–GA62, Glc–GA88, Glc–GA104, Glc–GA105, Glc–GA106, and Glc–GA107 (Table S2). In the MS/MS spectrum of compound 1 (Figure S4), the labeled product with m/z 563.2963 sequentially loses glucose and (CH3)2NH to yield fragments with m/z 401.2418 and m/z 356.1840, respectively. Additionally, fragment ions with m/z 338.1737 and m/z 310.1791, derived from the loss of H2O and HCOOH from the fragment ion with m/z 356.1840 [M–Glucose–(CH3)2NH+], are observed. Similarly, in the MS/MS spectrum of compound 2 (Figure S5), the labeled product with m/z 563.2963 also exhibits a loss of glucose, (CH3)2NH, H2O, and HCOOH within the high mass range. Based on this information, it could be inferred that both compound 1 and compound 2 possessed a lactone and a hydroxyl group on the A–ring. Considering the structural features of the candidates, the possible structures were Glc–GA7, Glc–GA62, Glc–GA88, Glc–GA104, Glc–GA105, Glc–GA106, and Glc–GA107. The possible fragmentation pathways of Glc–GA7–DMED, Glc–GA62–DMED, Glc–GA88–DMED, Glc–GA104–DMED, Glc–GA105–DMED, Glc–GA106–DMED, and Glc–GA107–DMED are shown in Figure S6, which are consistent with the characteristic fragmentation behaviors of compound 1 and compound 2, making it difficult to make a determination. To assist in the structural identification, a QSRR model was applied as a supplementary tool. From Table S4, the predicted retention indices (RIs) of GA105–DMED, GA62–DMED, GA104–DMED, GA106–DMED, GA107–DMED, GA7–DMED, and GA88–DMED are 535.52, 685.27, 716.94, 796.48, 826.57, 840.94, and 865.36, respectively. The differences between these predicted values and the experimental value of compound 1 (718.97) are −183.45, −33.70, −2.03, 77.51, 107.60, 121.97, and 146.39. Among them, the deviations of GA62–DMED, GA104–DMED, GA106–DMED, GA107–DMED, GA7–DMED, and GA88–DMED from the experimental value are greater than −69.58 but within the threshold of the prediction error of the model. The differences between these predicted values and the experimental value of compound 2 (760.66) are −225.14, −75.39, −43.72, 35.82, 65.91, 80.28, and 104.70. Among them, the deviations of GA104–DMED, GA106–DMED, GA107–DMED, GA7–DMED, and GA88–DMED from the experimental value are within the threshold of prediction error of the model. Furthermore, by comparing the secondary mass spectra of GA7 standards with compound 1 and compound 2 in plant samples, significant differences in the types and relative intensities of MS/MS fragment ions were observed (Figures S7 and S8). Therefore, compound 1 and compound 2 were not Glc–GA7. Based on the above analysis, the possible structures for compound 1 were Glc–GA62, Glc–GA88, Glc–GA104, Glc–GA106, or Glc–GA107, while the possible structures for compound 2 were Glc–GA88, Glc–GA104, Glc–GA106, or Glc–GA107.
For compounds 35, Glc–GA candidates with m/z 565.3119 corresponding to the product ion include Glc–GA20, Glc–GA69, Glc–GA70, Glc–GA84, Glc–GA4, Glc–GA40, Glc–GA51, Glc–GA61, Glc–GA119, Glc–GA45, and Glc–GA12 (Table S2). The MS/MS spectra of compounds 3, 4, and 5 contain fragment ions m/z 340.1907, m/z 312.1963, and m/z 294.1852, generated by the consecutive loss of H2O and HCOOH from the [M–Glucose–(CH3)2NH+] fragment ion (Figures S9–S11), indicating the presence of a lactone and a hydroxyl group on ring A of compounds 35. Based on the structural analysis of the candidates, possible structures were Glc–GA4, Glc–GA40, Glc–GA51, Glc–GA61, and Glc–GA119. Possible fragmentation pathways of Glc–GA40–DMED, Glc–GA51–DMED, Glc–GA61–DMED, Glc–GA119–DMED, and Glc–GA4–DMED are shown in Figure S12, all of which are positional isomers of glucose on ring A with similar fragmentation patterns. Considering the predicted retention indices (RI) (Table S4), it was speculated that the possible structures of compounds 35 were Glc–GA4, Glc–GA40, Glc–GA51, Glc–GA61, or Glc–GA119.
For compounds 69, Glc–GA candidates with m/z 579.2912 corresponding to the product ion include Glc–GA126, Glc–GA3, Glc–GA6, Glc–GA30, Glc–GA80, Glc–GA92, Glc–GA94, Glc–GA22, and Glc–GA68 (Table S2). The MS/MS spectra of compounds 6-9 contain fragment ions m/z 354.1700, m/z 326.1751, and m/z 310.1802, generated by the consecutive loss of H2O and HCOOH from the [M–Glucose–(CH3)2NH+] fragment ion (Figures S13–S16), indicating the presence of a lactone, hydroxyl group, or epoxy group on ring A of compounds 69. Based on the structural analysis of the candidates, possible structures were Glc–GA3, Glc–GA6, Glc–GA30, Glc–GA80, Glc–GA92, and Glc–GA68. Possible fragmentation pathways of Glc–GA92–DMED, Glc–GA80–DMED, Glc–GA68–DMED, Glc–GA3–DMED, Glc–GA6–DMED, and Glc–GA30–DMED are shown in Figure S17, which are consistent with the characteristic fragmentation behaviors of compounds 69. Considering the predicted retention indices (RI) (Table S4) and the fragmentation pathway of Glc–GA3 (Figure 3), compound 6 would be Glc–GA6, Glc–GA30, Glc–GA80, or Glc–GA68; compound 7 would be Glc–GA6, Glc–GA30, Glc–GA80, or Glc–GA68; compound 8 would be Glc–GA6 or Glc–GA30; and compound 9 would be Glc–GA6 or Glc–GA30.
For compound 10, the Glu–GA candidates corresponding to the marked product m/z 581.3068 include Glc–GA77, Glc–GA130, Glc–GA136, Glc–GA1, Glc–GA29, Glc–GA35, Glc–GA58, Glc–GA60, Glc–GA71, Glc–GA81, Glc–GA118, Glc–GA16, Glc–GA34, Glc–GA47, Glc–GA54, Glc–GA90, Glc–GA63, Glc–GA67, and Glc–GA131 (Table S2). The MS/MS spectrum of compound 10 in the high–polymer region contains fragment ions generated by the consecutive loss of H2O, HCOOH, and CO2+2H2O from the [M–Glucose–(CH3)2NH+] fragment ion (374.1962) (Figure S18), indicating the presence of one lactone and two hydroxyl groups on the A–ring of compound 10. Based on the structural characterization of the candidates, possible structures were Glc–GA16, Glc–GA34, Glc–GA47, Glc–GA54, and Glc–GA90. The potential fragmentation pathways of Glc–GA90–DMED, Glc–GA16–DMED, Glc–GA54–DMED, Glc–GA47–DMED, and Glc–GA34–DMED are shown in Figure S19, which are consistent with the characteristic fragmentation behaviors of compound 10. After filtering through the QSRR model, the possible structures were Glc–GA16, Glc–GA34, Glc–GA47, and Glc–GA54 (Table S4). Comparison of the secondary mass spectra of standard GA34 with that of compound 10 (Figure S20) reveals similar types of ion fragments and relative intensity trends. Since the remaining candidates, Glc–GA16, Glc–GA34, Glc–GA47, and Glc–GA54, were isomers with glucose or hydroxyl groups on the A–ring and exhibited similar fragmentation behavior, further determination was not possible. Therefore, the possible structures for compound 10 were Glc–GA16, Glc–GA34, Glc–GA47, or Glc–GA54.
For compound 11, the Glc–GA candidates corresponding to the marked product m/z 583.3225 include Glc–GA2 and Glc–GA82 (Table S2). In the MS/MS spectrum of compound 11, the marked product m/z 583.3225 sequentially loses glucose and (CH3)2NH to obtain m/z 421.2697 and m/z 376.2118, respectively. Meanwhile, the fragment ion m/z 376.2118 exhibits the behavior of losing H2O and HCOOH (Figure S21), indicating the presence of one lactone and one hydroxyl group on the A–ring of compound 11. Based on the structural analysis of the candidates, possible structures were Glc–GA2 and Glc–GA82. The potential fragmentation paths of Glc–GA2–DMED and Glc–GA82–DMED are shown in Figure S22, both of which are consistent with the characteristic fragmentation behaviors of compound 11. Combining the results of the QSRR model filtering (Table S4), we speculated that the possible structures for compound 11 are Glc–GA2 or Glc–GA82.
For compound 12, potential Glc–GA candidates corresponding to the marked product m/z 611.3174 include Glc–GA125, Glc–GA129, Glc–GA23, Glc–GA52, Glc–GA13, Glc–GA17, Glc–GA46, Glc–GA66, Glc–GA99 and Glc–GA102 (Table S2). The MS/MS spectrum of compound 12 includes a series of fragment ions in the higher mass region derived from the loss of 2H2O and HCOOH by the fragment ion [M–Glucose–(CH3)2NH+] (404.2068) (Figure S23), indicating the presence of one lactone and two hydroxyl groups on the A ring of compound 12. Based on the structural analysis of the candidates, the possible structure was Glc–GA52. After filtering through the QSRR model, the possible structure of compound 12 remained Glc–GA52 (Table S4). We analyzed the potential fragmentation pathways of Glc–GA52–DMED (Figure S24), which matched the characteristics of the characteristic fragmentation behaviors of compound 12. Therefore, we speculated that compound 12 was Glc–GA52.

3.6. Possible Locations of Potential Glc–GAs in the Metabolic Pathway of GAs

One hundred and thirty–six species of gibberellin named after GAn have been reported (Figure S25). In most higher plants, the biosynthesis of GAs typically involves three main steps (Figure 6). First, ent–kaurene is synthesized from trans–geranylgeranyl diphosphate (GGPP) as a precursor. The second step leads to the formation of GA12. Finally, GA12 can be converted into intermediate GAs and bioactive GAs through two pathways. One is the non–13–hydroxylation pathway, where GA12 is directly oxidized to produce intermediate GAs and bioactive GA4 and GA7. The other pathway is the 13–hydroxylation pathway. GA12 is initially oxidized to form GA53, which is then further oxidized to generate GA20, GA5, GA6, and other intermediate GAs, as well as bioactive GA1 and GA3. These deactivation and biosynthesis pathways of GAs collectively determine the endogenous content of GAs.
Based on the screening method developed in our study, a total of 12 potential Glc–GAs were screened from six plant samples (Table 1), and their possible structural information is presented in Figure 6. Among these, the positions of the corresponding free–form GAs for the six Glc–GAs (Glc–GA2, Glc–GA4, Glc–GA6, Glc–GA34, Glc–GA51, and Glc–GA61) in the GA metabolic pathway have been clearly elucidated. As shown in Figure 6, within the 13–hydroxylation pathway, intermediates GA2, GA4, GA34, GA51, GA61, and the bioactive GA4 might further undergo glucose conjugation to produce the corresponding glucose conjugated forms (Glc–GA2, Glc–GA6, Glc–GA34, Glc–GA51, and Glc–GA6). It is noteworthy that we have discovered a new Glucose–conjugated form of GA4, GA4–3–O–glucoside, which has only been reported in previous studies for the glucosyl ester form of GA4 (GA4 glucosyl ester) and the Glucose–conjugated product of 17–OH–GA4 (GA4–17–O–glucoside). This indicated that GA4 may also undergo glucose conjugation at the C–3 hydroxyl position. In addition, the intermediate GA6 generated in the non–13–hydroxylation pathway might also undergo glucose conjugation, resulting in Glc–GA6. Except for these six Glc–GA species, the positions in the GA metabolic pathway of other Glc–GAs discovered in this study are not clearly reported. Our results might be helpful for studying the biosynthesis and metabolic pathways of these Glc–GAs. These results suggested the widespread occurrence of Glucose–conjugated GAs in plants, where glucose conjugation and the deconjugation of GAs might be involved in regulating the levels of active GAs and play important roles in mediating the conversion between active GAs and intermediate GAs.

4. Conclusions

In this study, we successfully developed a structure–guided screening and identification strategy for Glc–GAs. This strategy utilizes characteristic MS/MS fragment ion patterns and incorporates QSRR techniques into the structure identification process. The method offers high sensitivity, requires lower plant sample amounts, and has no bias towards specific activities, enabling the efficient screening and structural identification of low–abundance and low–activity GA compounds. Using this approach, we screened 12 potential Glc–GAs in six plant samples and performed structural elucidation. The established analytical method facilitates the screening of more low–abundance and low–activity or inactive Glc–GAs, as the majority of metabolites in GA metabolic networks are inactive or have low activity. The work is highly meaningful for enhancing our understanding of GA metabolic pathways and investigating the mechanisms of GA functions.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/metabo14020096/s1, Table S1. Reported Glc–GAs from plant species [28,29,30,31,32,33,34,35,36,37,38,39,40]; Table S2. The exact mass of reported and potential Glc–GAs and the DMED/d4–DMED derivatization products; Table S3. The screened peak pairs of the proposed m/z of reported and potential Glc–GA–DMED/d4–DMED; Table S4. Comparison of the predicted RI of GAs with the experimental RI of Glc–GA candidates; Figure S1. Chemical structures of GA3 and Glc–GA3; Figure S2. Total ion chromatogram and extracted ion chromatograms (m/z 507.1872) of Glc–GA3 synthesis reaction products; Figure S3. HR MS/MS spectra of Glc–GA3 (m/z 507.18) at (A) 8.78 min, (B) 9.40 min, and (C) 10.17 min under negative mode; Figure S4. HR MS/MS spectra of compound 1; Figure S5. HR MS/MS spectra of compound 2; Figure S6. The proposed fragmentation pathways of Glc–GA7–DMED, Glc–GA62–DMED, Glc–GA88–DMED, Glc–GA104–DMED, Glc–GA105–DMED, Glc–GA106–DMED, and Glc–GA107–DMED; Figure S7. Comparison of HRMS/MS spectra of compound 1 and GA7; Figure S8. Comparison of HRMS/MS spectra of compound 2 and GA7; Figure S9. HR MS/MS spectra of compound 3; Figure S10. HR MS/MS spectra of compound 4; Figure S11. HR MS/MS spectra of compound 5; Figure S12. The proposed fragmentation pathways of Glc–GA40–DMED, Glc–GA51–DMED, Glc–GA61–DMED, Glc–GA119–DMED, and Glc–GA4–DMED; Figure S13. HR MS/MS spectra of compound 6; Figure S14. HR MS/MS spectra of compound 7; Figure S15. HR MS/MS spectra of compound 8; Figure S16. HR MS/MS spectra of compound 9; Figure S17. The proposed fragmentation pathways of Glc–GA92–DMED, Glc–GA80–DMED, Glc–GA68–DMED, Glc–GA3–DMED, Glc–GA6–DMED, and Glc–GA30–DMED; Figure S18. HR MS/MS spectra of compound 10; Figure S19. The proposed fragmentation pathways of Glc–GA90–DMED, Glc–GA16–DMED, Glc–GA54–DMED, Glc–GA47–DMED, and Glc–GA34–DMED; Figure S20. Comparison of HRMS/MS spectra of compound 10 and GA34; Figure S21. HR MS/MS spectra of compound 11; Figure S22. The proposed fragmentation pathways of Glc–GA2–DMED and Glc–GA82–DMED; Figure S23. HR MS/MS spectra of compound 12; Figure S24. The proposed fragmentation pathways of Glc–GA52–DMED; Figure S25. The structures of reported GAs.

Author Contributions

C.Z., conceptualization, methodology, investigation, visualization, writing—original draft.; W.-J.C., conceptualization, methodology, investigation, visualization; L.-C.J., methodology and investigation; T.Y., conceptualization, methodology, investigation, writing—review and editing. Y.-Q.F., supervision, conceptualization, methodology, investigation, writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

The work is supported by the National Natural Science Foundation of China (22204122, 22361132526, and 22274119).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article and Supplementary Materials.

Acknowledgments

We thank Xi-Zhou Hu from the Hubei Academy of Agricultural Sciences (Wuhan, China) for providing Vigna angularis. We thank Guangcun He from the college of life sciences, Wuhan University, for providing the rice seeds.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Hedden, P.; Sponsel, V. A Century of Gibberellin Research. J. Plant Growth Regul. 2015, 34, 740–760. [Google Scholar] [CrossRef] [PubMed]
  2. Plackett, A.R.G.; Wilson, Z.A. Gibberellins and Plant Reproduction. In Annual Plant Reviews; Hedden, P., Thomas, S.G., Eds.; Wiley: Hoboken, NJ, USA, 2016; Volume 49, pp. 323–358. ISBN 978-1-119-21042-9. [Google Scholar]
  3. Sembdner, G.; Atzorn, R. Plant Hormone Conjugation. Plant Mol. Biol. 1994, 26, 1459–1481. [Google Scholar] [CrossRef] [PubMed]
  4. Schneider, G.; Schliemann, W. Gibberellin Conjugates: An Overview. Plant Growth Regul. 1994, 15, 247–260. [Google Scholar] [CrossRef]
  5. Yokota, T.; Murofushi, N.; Takahashi, N.; Katsumi, M. Biological Activities of Gibberellins and Their Glycosides in Pharbitis Nil. Phytochemistry 1971, 10, 2943–2949. [Google Scholar] [CrossRef]
  6. Von Schirach–Szmigiel, L. Alterations in Endogenous Levels of Gibberellin–like Substances during Germination of Phaseolus vulgaris Seeds. Physiol. Plant. 1979, 46, 54–57. [Google Scholar] [CrossRef]
  7. Piotrowska, A.; Bajguz, A. Conjugates of Abscisic Acid, Brassinosteroids, Ethylene, Gibberellins, and Jasmonates. Phytochemistry 2011, 72, 2097–2112. [Google Scholar] [CrossRef]
  8. Knöfel, H.-D.; Schwarzkopf, E.; Müller, P.; Sembdner, G. Enzymic Glycosylation of Gibberellins. J. Plant Growth Regul. 1984, 3, 127–140. [Google Scholar] [CrossRef]
  9. Sakamoto, T.; Miura, K.; Itoh, H.; Tatsumi, T.; Ueguchi-Tanaka, M.; Ishiyama, K.; Kobayashi, M.; Agrawal, G.K.; Takeda, S.; Abe, K.; et al. An Overview of Gibberellin Metabolism Enzyme Genes and Their Related Mutants in Rice. Plant Physiol. 2004, 134, 1642–1653. [Google Scholar] [CrossRef] [PubMed]
  10. Hiraga, K.; Yamane, H.; Takahashi, N. Biological Activity of Some Synthetic Gibberellin Glycosyl Esters. Phytochemistry 1974, 13, 2371–2376. [Google Scholar] [CrossRef]
  11. Yamane, H.; Yamaguchi, I.; Yokota, T.; Murofushi, N.; Takahashi, N.; Katsumi, M. Biological Activities of New Gibberellins A30–A35 and A35 Glycoside. Phytochemistry 1973, 12, 255–261. [Google Scholar] [CrossRef]
  12. Baidsx, G.W.M. Formation of Bound Gibberellins. Planta 1971, 99, 290–301. [Google Scholar] [CrossRef]
  13. Sembdner, G.; Schliemann, W.; Schneider, G. Biochemical and Physiological Aspects of Gibberellin Conjugation. In Gibberellins; Takahashi, N., Phinney, B.O., MacMillan, J., Eds.; Springer: New York, NY, USA, 1991; pp. 249–263. ISBN 978-1-4612-7754-5. [Google Scholar]
  14. Stephan, M.; Bangerth, F.; Schneider, G. Transport and Metabolism of Exogenously Applied Gibberellins to Malus Domestica Borkh. Cv. Jonagold. Plant Growth Regul. 2001, 33, 77–85. [Google Scholar] [CrossRef]
  15. Ostrowski, M.; Jakubowska, A. Udp–Glycosyltransferases of Plant Hormones. Adv. Cell Biol. 2014, 4, 43–60. [Google Scholar] [CrossRef]
  16. Piccoli, P.; Lucangeli, C.D.; Schneider, G.; Bottini, R. Hydrolysis of [17,17–2H2] Gibberellin A20–Glycoside and [17,17–2H2]Gibberellin A20–Glycosyl Ester by Azospirillum Lipoferum Cultured in a Nitrogen–Free Biotin–Based Chemically–Defined Medium. Plant Growth Regul. 1997, 23, 179–182. [Google Scholar] [CrossRef]
  17. Moritz, T.; Monteiro, A.M. Analysis of Endogenous Gibberellins and Gibberellin Metabolites from Dalbergia dolichopetMaby Gas Chromatography–Mass Spectrometry and High–Performance Liquid Chromatography–Mass Spectrometry. Planta 1994, 193, 1–8. [Google Scholar] [CrossRef]
  18. Rivier, L.; Gaskin, P.; Albone, K.S.; MacMillan, J. GC–MS Identification of Endogenous Gibberellins and Gibberellin Conjugates as Their Permethylated Derivatives. Phytochemistry 1981, 20, 687–692. [Google Scholar] [CrossRef]
  19. Binks, R.; MacMillan, J.; Pryce, R.J. Plant Hormones—VIII: Combined Gas Chromatography–Mass Spectrometry of the Methyl Esters of Gibberellins A1 to A24 and Their Trimethylsilyl Ethers. Phytochemistry 1969, 8, 271–284. [Google Scholar] [CrossRef]
  20. Schneider, G.; Sembdner, G.; Jensen, E.; Bernhard, U.; Wagenbreth, D. GC–MS Identification of Native Gibberellin–O–Glycosides in Pea Seeds. J. Plant Growth Regul. 1992, 11, 15–18. [Google Scholar] [CrossRef]
  21. Schneider, G.; Schmidt, J.; Phinney, B.O. GC–MS Identification of GA20–13–O–Glycoside Formed from GA20 in Normal Plants and Dwarf–1 Mutants of Zea mays L. J. Plant Growth Regul. 1987, 5, 217–223. [Google Scholar] [CrossRef]
  22. Senns, B.; Fuchs, P.; Schneider, G. GC–MS Quantification of Gibberellin A20–13–O–Glycoside and Gibberellin A8–2–O–Glycoside in Developing Barley Caryopses. Phytochemistry 1998, 48, 1275–1280. [Google Scholar] [CrossRef]
  23. Moritz, T. The Use of Combined Capillary Liquid Chromatography/Mass Spectrometry for the Identification of a Gibberellin Glycosyl Conjugate. Phytochem. Anal. 1992, 3, 32–37. [Google Scholar] [CrossRef]
  24. Cai, W.-J.; Zeng, C.; Zhang, X.-Y.; Ye, T.; Feng, Y.-Q. A Structure–Guided Screening Strategy for the Discovery and Identification of Potential Gibberellins from Plant Samples Using Liquid Chromatography–Mass Spectrometry Assisted by Chemical Isotope Labeling. Anal. Chim. Acta 2021, 1163, 338505. [Google Scholar] [CrossRef] [PubMed]
  25. Hao, Y.-H.; Zhang, Z.; Wang, L.; Liu, C.; Lei, A.-W.; Yuan, B.-F.; Feng, Y.-Q. Stable Isotope Labeling Assisted Liquid Chromatography–Electrospray Tandem Mass Spectrometry for Quantitative Analysis of Endogenous Gibberellins. Talanta 2015, 144, 341–348. [Google Scholar] [CrossRef]
  26. Zhu, Q.-F.; Yan, J.-W.; Gao, Y.; Zhang, J.-W.; Yuan, B.-F.; Feng, Y.-Q. Highly Sensitive Determination of Fatty Acid Esters of Hydroxyl Fatty Acids by Liquid Chromatography–Mass Spectrometry. J. Chromatogr. B 2017, 1061–1062, 34–40. [Google Scholar] [CrossRef]
  27. Schneider, G.; Sembdner, G.; Schreiber, K. ChemInform Abstract: GIBBERELLINE 26. MITT. SYNTHESE VON GIBBERELLIN–A3–BETA–D–GLYCOPYRANOSIDEN. Chem. Inf. 1975, 6, 226. [Google Scholar] [CrossRef]
  28. Noma, M.; Huber, J.; Ernst, D.; Pharis, R.P. Quantitation of Gibberellins and the Metabolism of [3H]Gibberellin A1 during Somatic Embryogenesis in Carrot and Anise Cell Cultures. Planta 1982, 155, 369–376. [Google Scholar] [CrossRef]
  29. Schliemann, W.; Schaller, B.; Jensen, E.; Schneider, G. Native Gibberellin-O-Glucosides from Mature Seeds of Phaseolus Coccineus. Phytochemistry 1993, 35, 35–38. [Google Scholar] [CrossRef]
  30. Schneider, G.; Fuchs, P.; Schmidt, J. Evidence for the Direct 2β- and 3β-hydroxylation of [2H2]GA20-13-O-[6′-2H2]Glucoside in Seedlings of Phaseolus coccineus. Physiologia Plantarum 2002, 116, 144–147. [Google Scholar] [CrossRef] [PubMed]
  31. Tamura, S.; Takahashi, N.; Yokota, T.; Murofushi, N.; Ogawa, Y. Isolation of Water-Soluble Gibberellins from Immature Seeds of Pharbitis Nil. Planta 1967, 78, 208–212. [Google Scholar] [CrossRef]
  32. Yamaguchi, I.; Yokota, T.; Yoshida, S.; Takahashi, N. High Pressure Liquid Chromatography of Conjugated Gibberellins. Phytochemistry 1979, 18, 1699–1702. [Google Scholar] [CrossRef]
  33. Hasegawa, M.; Nakajima, M.; Takeda, K.; Yamaguchi, I.; Murofushi, N. A Novel Gibberellin Glucoside 16α,` 17-Dihydroxy-16,17-Dihydrogibberellin A4-17-O-β-d-Glucopyranoside, from Rice Anthers. Phytochemistry 1994, 37, 629–634. [Google Scholar] [CrossRef]
  34. Hiraga, K.; Kawabe, S.; Yokota, T.; Murofushi, N.; Takahashi, N. Isolation and Characterization of Plant Growth Substances in Immature Seeds and Etiolated Seedlings of Phaseolus vulgaris. Agric. Biol. Chem. 1974, 38, 2521–2527. [Google Scholar] [CrossRef]
  35. Yokota, T.; Takahashi, N.; Murofushi, N.; Tamura, S. Isolation of Gibberellins A26 and A27 and Their Glucosides from Immature Seeds of Pharbitis Nil. Planta 1969, 87, 180–184. [Google Scholar] [CrossRef] [PubMed]
  36. Yokota, T.; Takahashi, N.; Murofushi, N.; Tamura, S. Structures of new gibberellin glucosides in immature seeds of pharbitis nil. Tetrahedron Lett. 1969, 10, 2081–2084. [Google Scholar] [CrossRef]
  37. Yamane, H.; Yamaguchi, I.; Murofushi, N.; Takahashi, N. Isolation and Structures of Gibberellin A35 and Its Glucoside from Immature Seed of Cytisus scoparius. J. Agric. Chem. Soc. Jpn. 1974, 38, 649–655. [Google Scholar] [CrossRef]
  38. Hiraga, K.; Yokota, T.; Murofushi, N.; Takahashi, N. Isolation and Characterization of a Free Gibberellin and Glucosyl Esters of Gibberellins in Mature Seeds of Phaseolus vulgaris. Agric. Biol. Chem. 1972, 36, 345–347. [Google Scholar] [CrossRef]
  39. Yamaguchi, I.; Kobayashi, M.; Takahashi, N. Isolation and Characterization of Glucosyl Esters of Gibberellin A 5 and A 44 from Immature Seeds of Pharbitis purpurea. Agric. Biol. Chem. 1980, 44, 1975–1977. [Google Scholar] [CrossRef]
  40. Lorenzi, R.; Horgan, R.; Heald, J.K.; Lorenzi, R.; Horgan, R.; Heald, J.K. Gibberellins in Picea Sitchensis Carriere: Seasonal Variation and Partial Characterization. Planta 1975, 126, 75–82. [Google Scholar] [CrossRef]
Figure 1. Overview of the workflow for the screening and identification of potential Glc–GAs from plant species.
Figure 1. Overview of the workflow for the screening and identification of potential Glc–GAs from plant species.
Metabolites 14 00096 g001
Figure 2. Derivatization of GA3–3–O–glucoside with DMED/d4–DMED.
Figure 2. Derivatization of GA3–3–O–glucoside with DMED/d4–DMED.
Metabolites 14 00096 g002
Figure 3. (A). Total ion chromatogram and extracted ion chromatograms of synthetic Glc–GA3 after DMED labeling (red: m/z 579.2912, blue: m/z 417.2384). HR MS/MS spectra of Glc–GA3–DMED (m/z 579) at (B) 5.57 min, (C) 6.25 min, GA3–DMED (m/z 417) at (D) 8.76 min under positive mode.
Figure 3. (A). Total ion chromatogram and extracted ion chromatograms of synthetic Glc–GA3 after DMED labeling (red: m/z 579.2912, blue: m/z 417.2384). HR MS/MS spectra of Glc–GA3–DMED (m/z 579) at (B) 5.57 min, (C) 6.25 min, GA3–DMED (m/z 417) at (D) 8.76 min under positive mode.
Metabolites 14 00096 g003
Figure 4. The proposed fragmentation pathways of DMED–GA3–3–O–Glucoside.
Figure 4. The proposed fragmentation pathways of DMED–GA3–3–O–Glucoside.
Metabolites 14 00096 g004
Figure 5. Chemical isotope–labeled Canavalia gladiata sample analyzed by LC–HRMS under full scan mode. (A) Total ion chromatogram of DMED/d4–DMED labeled Canavalia gladiata sample and extracted ion chromatograms of representative Glc–GA–DMED/d4–DMED channel at m/z 563.2963 and 567.3214. (B) Mass spectra at 10.32 min.
Figure 5. Chemical isotope–labeled Canavalia gladiata sample analyzed by LC–HRMS under full scan mode. (A) Total ion chromatogram of DMED/d4–DMED labeled Canavalia gladiata sample and extracted ion chromatograms of representative Glc–GA–DMED/d4–DMED channel at m/z 563.2963 and 567.3214. (B) Mass spectra at 10.32 min.
Metabolites 14 00096 g005
Figure 6. Stereochemical structures of potential Glc–GA compounds and their locations in the GA biosynthetic and inactive pathways. Compounds in red are possible Glc–GAs in a known metabolic pathway. Compounds in blue are Glc–GAs in an unknown metabolic pathway. Solid lines indicate known GA metabolic pathways. Dotted lines indicate possible GA glycosylation processes. The wavy bonds represent unclear anomers. The superscript number in the compound name (Glc–GAX) is the corresponding number of possible Glucose–conjugated forms of GAX.
Figure 6. Stereochemical structures of potential Glc–GA compounds and their locations in the GA biosynthetic and inactive pathways. Compounds in red are possible Glc–GAs in a known metabolic pathway. Compounds in blue are Glc–GAs in an unknown metabolic pathway. Solid lines indicate known GA metabolic pathways. Dotted lines indicate possible GA glycosylation processes. The wavy bonds represent unclear anomers. The superscript number in the compound name (Glc–GAX) is the corresponding number of possible Glucose–conjugated forms of GAX.
Metabolites 14 00096 g006
Table 1. Summarized information of 12 Glc–GA candidates.
Table 1. Summarized information of 12 Glc–GA candidates.
IDCandidatesTheoretical m/zDerivatized FormulaNative
Formula
RTRIIdentification
DMED
–Labeled
d4–DMED
–Labeled
1Glc–GA5, 31, 95, 96, 108, 109, 117, 121, 122,7, 11, 62, 88, 104, 105, 106, 107563.2963 567.3214 C29H44O9N2C25H32O1010.32 718.97 Glc–GA62, 88, 104,106, 107 a,b
212.10 760.66 Glc–GA62, 88, 104,106, 107 c
3Glc–GA20, 69, 70, 84, 4, 40, 51, 61, 119, 45565.3119 569.3370 C29H46O9N2C25H34O10 8.70 682.35 Glc–GA4, 40, 51, 61, 119 c
410.79 729.98 Glc–GA4, 40, 51, 61, 119 d
513.01 781.97 Glc–GA4, 40, 51, 61, 119 a,b,c,e
6Glc–GA126, 3, 6, 30, 80, 92, 94, 22, 68579.2912 583.3163 C29H44O10N2C25H32O117.47 655.56 Glc–GA6, 30, 80, 68 a
78.20 671.46 Glc–GA6, 30, 80, 68 e
88.94 687.58 Glc–GA6, 30 e
911.12 737.70 Glc–GA6, 30 f
10Glc–GA77, 130, 136, 1, 29, 35, 58, 60, 71, 81, 118, 16, 34, 47, 54, 90, 63, 67, 131 581.3068 585.3320 C29H46O10N2C25H34O1110.77 729.51 Glc–GA16, 34, 47, 54 a
11Glc–GA2, 82583.3225 587.3476 C29H48O10N2C25H36O118.70 682.35 Glc–GA2, 82 e
12Glc–GA125, 129, 23, 52, 13, 17, 46, 66, 99, 102611.3174 615.3425 C30H48O11N2C26H36O129.55 700.94 Glc–GA52 e
Plant species: a Canavalia gladiata; b Phaseolus vulgaris; c Pisum sativum; d Oryza sativa ssp. japonica cv. Nipponbare; e Vicia faba; f Vigna angularis.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zeng, C.; Cai, W.-J.; Jiang, L.-C.; Ye, T.; Feng, Y.-Q. Advancing Glucose Conjugated Gibberellins Discovery: A Structure–Oriented Screening and Identification Method for Unraveling Gibberellin Metabolites in Plants. Metabolites 2024, 14, 96. https://doi.org/10.3390/metabo14020096

AMA Style

Zeng C, Cai W-J, Jiang L-C, Ye T, Feng Y-Q. Advancing Glucose Conjugated Gibberellins Discovery: A Structure–Oriented Screening and Identification Method for Unraveling Gibberellin Metabolites in Plants. Metabolites. 2024; 14(2):96. https://doi.org/10.3390/metabo14020096

Chicago/Turabian Style

Zeng, Chen, Wen-Jing Cai, Liu-Cheng Jiang, Tiantian Ye, and Yu-Qi Feng. 2024. "Advancing Glucose Conjugated Gibberellins Discovery: A Structure–Oriented Screening and Identification Method for Unraveling Gibberellin Metabolites in Plants" Metabolites 14, no. 2: 96. https://doi.org/10.3390/metabo14020096

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop