ORIGINAL RESEARCH article

Front. Nutr., 29 August 2022

Sec. Nutrition and Food Science Technology

Volume 9 - 2022 | https://doi.org/10.3389/fnut.2022.928805

Molecular mechanisms of flavonoid accumulation in germinating common bean (Phaseolus vulgaris) under salt stress

  • 1. Legume Crop Laboratory, Agricultural College, Heilongjiang Bayi Agricultural University, Daqing, China

  • 2. Qiqihar Branch of Heilongjiang Academy of Agricultural Sciences, Qiqihar, China

  • 3. Crop Resources Institute of Heilongjiang Academy of Agricultural Sciences, Harbin, China

  • 4. National Cereals Technology Engineering Research Center, Daqing, China

Article metrics

View details

25

Citations

4,1k

Views

1,3k

Downloads

Abstract

Flavonoids are important secondary metabolites, active biomolecules in germinating beans, and have prominent applications in food and medicine due to their antioxidant effects. Rutin is a plant flavonoid with a wide biological activity range. In this study, flavonoid (rutin) accumulation and its related molecular mechanisms in germinating common bean (Phaseolus vulgaris) were observed at different time points (0–120 h) under salt stress (NaCl). The rutin content increased from germination onset until 96 h, after which a reducing trend was observed. Metabolome analysis showed that salt stress alters flavonoid content by regulating phenylpropanoid (ko00940) and flavonoid (ko00941) biosynthesis pathways, as well as their enzyme activities, including cinnamyl-alcohol dehydrogenase (CAD), peroxidase (POD), chalcone isomerase (CHI), and flavonol synthase (FLS). The RNA-seq and quantitative real-time PCR (qRT-PCR) analyses also showed that these two pathways were linked to changes in flavonoid content following salt treatment. These results reveal that salt stress effectively enhanced rutin content accumulation in germinating beans, hence it could be employed to enhance the functional quality of germinating common beans.

Introduction

The common bean (Phaseolus vulgaris L.) is an annual legume crop that has been cultivated in temperate and semitropical regions for ~8,000 years (1). Common beans constitute an essential part of the human diet not only because of their protein-rich seed but also because they contain active phytochemicals, which are beneficial for health (2), including proteases, galacto-oligosaccharides, phytic acid, and flavonoids (3).

Flavonoids are an abundant and biologically active family of natural plant products that offer a range of health benefits (4). For instance, several flavonoids have anticancer functions, both in vitro and in vivo (5, 6). Some of the reported mechanisms by which flavonoids exert their anticancer functions include (1) preventing new cancer cell development, (2) restraining carcinogens that turn into activation sites, and (3) preventing compound metabolism and reducing their toxicity (7). Recently, with the ever-growing food quality concerns, the international market has shown a growing interest in flavonoid production, with high demand and a growing compound rate of 16.5%.

Plants are the most common material used for flavonoid extraction (7). These flavonoids have a C6–C3–C6 carbon skeleton with two phenyl aromatic rings, along with a heterocyclic ring (8). Flavonoids can be split into several subclasses according to basic skeleton substitution and B-ring attachment, such as chalcones, stilbenes, aurones, flavanones, flavones, isoflavones, phlobaphenes, dihydroflavonols, flavonols, leucoanthocyanidins, proanthocyanidins, and anthocyanins (9, 10). Rutin (C27H30O16) is a relatively large and representative flavonoid in leguminous crops (11–13), which could also be used as a representative flavonoid for studying plant nutrition (14). In plants, certain flavonoids, including rutin, have been proven to play an important role in plant growth, development, hormone signaling, facilitating pollen–tube germination, plant–microorganism interactions, and biotic and abiotic stresses (15, 16).

Legumes are potential raw materials used for producing flavonoids for nutritional and dietary applications (4). However, abiotic stress (e.g., drought, salt, and cold) can affect the flavonoid biosynthetic pathway by regulating flavonoid transport and the expression and accumulation of related proteins, finally affecting the flavonoid content in plants (17). Previous studies have shown that salt stress significantly inhibits the germination of legume crops (including Lathyrus sativus, Vicia sativa, Vigna radiata, and Vigna unguiculata) but increases the flavonoid content in their seeds (18). Besides, another study reported significant differences in the flavonoid content and anti-inflammatory effects of soybean when seeds were germinated under different abiotic stress treatments (19). It has also been reported that abiotic stress could increase the content of total phenolic compounds and total flavonoids in soybean (Glycine max) and mung bean (V. radiata) sprouts (19, 20). Therefore, the regulation of flavonoid biosynthesis is considered an important regulatory mechanism of legume crops in response to salt stress (18). In another study, Gu et al. (21) demonstrated that moderate drought stress treatment could increase flavonoid content in plants. From a nutritional viewpoint, it could be concluded that exposing plants to appropriate stress treatment is an effective strategy for increasing the flavonoid content in food. However, few studies have reported the correlation between germinating common beans and flavonoids. This study determined the accumulation of rutin (used as the representative flavonoid in the common bean) in different common bean tissues at different germination stages under salt stress. Furthermore, the effect and mechanism of salt stress on total flavonoid accumulation in common bean sprouts were evaluated through metabolome and transcriptome analysis. The results presented herein provide a theoretical basis and new insights for improving the flavonoid content of common bean sprouts.

Materials and methods

Common bean materials and treatments

A local common bean cultivar (Longjiang Black Yun) was used as testing material, and its seeds were provided by the National Coarse Cereals Technology Engineering Research Center (Heilongjiang, China). The seeds were surface-sterilized with sodium hypochlorite (NaClO) and washed three times with distilled water. Germination was done in a germination box (BSC-250, Boxun, Shanghai, China) at a constant temperature of 25°C without light. Salt stress was simulated using 70 mmol/L NaCl (22) and was added at different times, as illustrated in Supplementary Table 1. The seeds were exposed to salt stress for 0, 24, 48, 72, 96, and 120 h for sprouting, and the germinating common beans were sampled to measure the rutin content. Different tissues with the highest flavonoid content in the germinating common beans were isolated since they were most suitable for further research. Flavonoid extraction was done using the method described by Liu et al. (23), while the rutin content was determined using HPLC (AB SCIEX, shimadzuLc-20A, MA, USA) from the Customs Quality Inspection Center (Dalian, Liaoning, China), using rutin (C27H30O16) (SR8250, Solarbio, Beijing, China) as the reference (24–26). The standard rutin curve is shown in Supplementary Figure 1.

Metabolome and transcriptome analysis

Bean sprouts, without cotyledons, in the 3d+S treatment were taken at three time points (3d+S-0h, 3d+S-12h, and 3d+S-24h) and used as test samples for metabolome and transcriptome analysis, with three biological replicates at each time point. The hours (0, 12, and 24 h) represent the time salt stress was introduced at 3 days. The broadly targeted metabolome was determined using the UPLC-ESI-MS/MS system (Nexera X2, SHIMADZU, Japan) and ESI-Q TRAP-MS/MS (4500 QTRAP, Applied Biosystems, USA) from Matwell Biotechnology Company (Wuhan, China) (Supplementary Figure 2). The HPLC conditions included liquid-phase chromatographic column (Waters ACQUITY UPLC HSS T3 C18 1.8 μm, 2.1 mm * 100 mm), mobile phase (phase A was ultrapure water with 0.04% acetic acid while phase B was acetonitrile with 0.04% acetic acid), flow rate (0.35 ml/min), column temperature (40°C), and injection volume (4 μl). The mass spectrometry conditions included an electrospray ionization temperature of 550°C, a mass spectrometer voltage of 5,500 V, and a curtain gas pressure of 30 psi. Each ion pair in a triple quadrupole was scanned according to the optimized declustering potential and collision energy (27). The reagents (including methanol, acetonitrile, and ethanol) were bought from Merck (Darmstadt, Germany), while the standards were bought from Sigma (Sigma-Aldrich, Shanghai, China). The blank reagent sample was run during the detection process in order to clean the detection residue. The qualitative analysis of metabolites in samples was based on the MetWare database, so the isotopic and repeating signals (including K+, Na+, and NH) were removed. The metabolite identification method was based on a self-built database, including standard products and experimental samples. The secondary mass spectrometry data were collected based on the mass spectrometry MIM-EPI mode, and the unique Q1 (molecular ion), Q3 (characteristic fragment ion), RT (chromatographic retention time), DP (declustering potential), and CE (collision energy) of each metabolite were constructed by combining with manual analysis of the spectrum. Subsequently, the experimental samples were quantitatively detected using the multiple reaction monitoring (MRM) mode, which could eliminate nontarget ion interference and increase the accuracy and reproducibility of the quantification (27). After screening, the Analyst 1.6.3 software was used to process the mass spectrometry data, and the MultiQuant software to open the sample off-machine mass spectrometry file to integrate and calibrate the chromatographic peaks. The peak area of each chromatographic peak represents the relative content of the corresponding substance. The quality control (QC) in this study was used with mixed samples and internal standards (L-2-chlorophenylalanine, CAS: 103616-89-3) using the processing method applied for the samples. During the instrumental analysis, a QC sample was inserted into every 10 detection and analysis samples to monitor the repeatability of the analysis process. Fold change (FC) and VIP value analysis were used to find differentially altered metabolites (DAMs), and the screening conditions for DAMs were VIP ≥ 1 and FC ≥ 2 or FC ≤ 0.5 (28). The Kyoto Encyclopedia of Genes and Genomes (KEGG) database was used to annotate the information from the DAMs (29), with the p-value set to <0.05.

Also, the bean sprouts without cotyledons at three time points (3d+S-0h, 3d+S-12h, and 3d+S-24h) were used as transcriptome samples with three biological replicates at each time point. The transcriptome of samples was determined using HiSeq X Ten (Illumina, CA, USA) from the Biomarker Company (Beijing, China), and the transcriptome data were analyzed on the BMKCloud platform. The TopHat2 software (30) was used for sequence alignment of clean reads with a reference genome, while fragments per kilobase of transcript per million fragments mapped was used to measure the transcript or gene expression levels (31). The FC and false discovery rate (FDR) were calculated to determine the differentially expressed genes (DEGs). The value of FC was the average value in each group, which was calculated by merging three biological replicates in one treatment. The screening conditions for DEGs were FC ≥ 2 and FDR < 0.01 (28). The KEGG database was used to annotate the information from the DEGs (29), with the p-value set to < 0.05.

Physiology and qRT-PCR analysis

Enzyme activities, including cinnamyl-alcohol dehydrogenase (CAD) (32), peroxidase (POD) (33), chalcone isomerase (CHI) (34), and flavonol synthase (FLS) (35), in the enriched pathways were determined with an ELISA Kit (Michy, Suzhou, China), and data were acquired using a microplate reader (SpectraMax® 190, Molecular Devices, CA, USA). The data were then analyzed using ANOVA at a significance level of p < 0.05 on the SPSS19.0 platform.

The RNA of each sample was extracted with the TreliefTM RNA Pure Plant Kit (TSP412, Tsingke, China), and a NanoDrop (OneC, Thermo, USA) was used to detect quality and concentration. Quality RNA was reverse-transcribed into cDNA with the Evo M-MLV RT Premix (Accurate-Biotechnology, AG11706, Hunan, China). Primers of DEGs enriched in CAD (Phvul.002G144800 and Phvul.003G029500), POD (Phvul.001G143300, Phvul.007G008400 and Phvul.008G249900), and FLS (Phvul.003G216600) were designed using Primer Premier 5.0, while Pvactin-11 was selected as the reference gene; the primer sequences are shown in Supplementary Table 2. The qRT-PCR program was conducted on a Light Cycler 480II (Roche Diagnostics, Switzerland), using a Hieff UNICON® Universal qPCR SYBR Green Master Mix (Yeason, 11184ES03, Shanghai, China). Subsequently, the relative gene expression level was computed using the 2−ΔCt method (36).

Results

The rutin content in germinating common bean

The rutin content of common beans without germination was significantly lower (1.99 mg/ml; p < 0.05) at different germination times compared to common beans with germination. The rutin content of common bean sprouts at different germinating times increased with time and peaked at 96 h, after which a reducing trend was observed. Compared with other germination times, the periods from 72 to 120 h were more suitable for germinating bean sprouts due to the rutin content, while 96 h could be considered optimal (Figure 1A). Bean sprouts at 96 h were divided into two parts for the determination of flavonoid content (Figure 1B); cotyledons and sprouts (including epicotyl and hypocotyl). Sprouts had a significantly (p < 0.05) higher rutin content than cotyledons (Figure 1C). The rutin content in different tissues revealed that sprouts are a suitable target tissue for further research.

Figure 1

Effects of salt stress on rutin content

The rutin contents under different salt exposure times were determined in this study. The result revealed differences in rutin content accumulation at different periods of salt addition. Compared to the treatment without salt (4 days), the rutin content of 2d+S and 3d+S treatments increased, while that of treatments 0d+S and 1d+S decreased (Figure 2). Short-term salt stress, thus, accelerated rutin content accumulation more than long-term stress, and the 3d+S treatment could be used as an optimal treatment to study the mechanism of rutin content accumulation.

Figure 2

Metabolome and transcriptome analysis

The study added a time point (12 h) to enhance broad target metabolome accuracy, with three time points (3d+S-0h, 3d+S-12h, and 3d+S-24h) in the 3d+S treatment, and each time point had three biological replicates. The raw data were uploaded to the Zenodo database (doi.org/10.5281/zenodo.6820497). A total of 732 metabolites were detected; the detailed information, including Q1, Q3, Rt, molecular weight, and ionization model, is shown in Supplementary Table 3. The overlay analysis of the QC-TIC diagram and the sample multipeak detection of MRM analysis diagram (Figure 3) showed that the data had good reproducibility and reliability, hence they could be used for further analysis.

Figure 3

In the comparison of 3d+S-0h vs. 3d+S-24h, 86 DAMs were increased while 55 were decreased (Supplementary Table 4). In the comparison of 3d+S-0h vs. 3d+S-12h, 62 DAMs were increased, while 63 were decreased (Supplementary Table 5 and Table 1). In this study, the DAMs in 3d+S-0h vs. 3d+S-24h and 3d+S-0h vs. 3d+S-12h had been combined together for analysis. A total of 189 DAMs were obtained after removing 77 duplicate DAMs, while 68 DAMs were associated with flavonoids after removing 36 duplicate flavonoid DAMs. In these, 68 DAMs were associated with flavonoids, five types of flavonoids (i.e., flavone, flavonol, flavonoid, flavanone, and isoflavone) had different expression levels (Figure 4), and each compound tested by HPLC-MS is listed in Supplementary Table 6. Rutin made up one DAM metabolome, and its content increased significantly. The KEGG enrichment was analyzed for two comparisons. In 3d+S-0h vs. 3d+S-12h, eight pathways had a p-value < 0.05 and could thus be used as enrichment pathways (Table 2). In 3d+S-0h vs. 3d+S-24h, five pathways had a p-value < 0.05 and could thus be used as enrichment pathways (Table 3). Four candidate pathways were found in the metabolomics analysis, namely, phenylpropanoid biosynthesis (ko00940), flavonoid biosynthesis (ko00941), flavone and flavonol biosynthesis (ko00944), and biosynthesis of secondary metabolites (ko01110).

Table 1

GroupsTotal DAMsUp DAMsDown DAMs
3d+S-0h vs. 3d+S-12h1418655
3d+S-0h vs. 3d+S-24h1256263

The number of differentially altered metabolites (DAMs).

Figure 4

Table 2

No.Ko IDKEGG pathwayP-valueDAMs
1Ko00944Flavone and flavonol biosynthesis0.000pmb0605; pme0199; pme0197; pme3211; pme1622; pme0089
2Ko00941Flavonoid biosynthesis0.001pme0199; pmf0108; pme0089; pme1521; pme3461; pme1201
3Ko00945Stilbenoid, diarylheptanoid and gingerol biosynthesis0.015pme1458
4Ko00940Phenylpropanoid biosynthesis0.018pme3125; pme0030; pma0149; pme3305
5Ko01110Biosynthesis of secondary metabolites0.021pme0199; pme3261; pme0030; pme3125; pme2555; pmf0558; pme1622; pme3266; pme0089; pme1521; pme1692; pme1201; pme0197; pmf0248; pme3211; pme1002; pme1458; pme3233; pme0026; pmb0449; pme3305
6Ko00965Betalain biosynthesis0.027pme0030
7Ko00230Purine metabolism0.027pme1119; pme0193; pme1692; pmb0514; pmb2684; pme2555; pmd0023
8Ko00910Nitrogen metabolism0.033pme0193

The enriched pathways of differentially altered metabolites (DAMs) compared 3d+S-0h vs. 3d+S-12h in KEGG analysis.

Table 3

No.Ko IDKEGG pathwayP-valueDAMs
1Ko00944Flavone and flavonol biosynthesis0.000pme0374; pme1622; pme3211; pme0197
2Ko00941Flavonoid biosynthesis0.000pme1201; pmf0345; pme3509; pme1521; pmf0108
3Ko01110Biosynthesis of secondary metabolites0.001pme1695; pme2801; pme2993; pme0006; pmf0558; pme0303; pme3266; pme3125; pme2122; pme3305; pme0197; pme1622; pme0355; pmf0117; pme2155; pme1201; pmf0345; pme3252; pme1841; pme3233; pme1521; pme3211; pme3261; pme0013; pme3553
4Ko00940Phenylpropanoid biosynthesis0.003pme0303; pmb0242; pme2993; pma0149; pme3125; pme1695; pme3305
5Ko00943Isoflavonoid biosynthesis0.028pme3233; pme3509; pme3208; pme3252; pme0355; pmf0117; pme3266; pme3261

The enriched pathways of differentially altered metabolites (DAMs) compared 3d+S-0h vs. 3d+S-24h in KEGG analysis.

The samples of three time points, i.e., 3d+S-0h, 3d+S-12h, and 3d+S-24h in the 3d+S treatment, were used for RNA-seq. The data quality was sufficient for follow-up analysis (Supplementary Table 7) and data were uploaded to the NCBI database (PRJNA746732) (Supplementary Table 8). In the comparisons of 3d+S-0h vs. 3d+S-12h, and 3d+S-0h vs. 3d+S-24h, a total of 1,230 and 1,450 DEGs were found, of which 812 and 904 DEGs were upregulated, respectively, and 418 and 546 were downregulated, respectively (Table 4). All DEGs were analyzed in KEGG enrichment, respectively. A total of 16 pathways were enriched in 3d+S-0h vs. 3d+S-12h (Table 5, while 13 pathways were enriched in 3d+S-0h vs. 3d+S-24h (Table 6). In the transcriptome analysis of 3d+S-0h vs. 3d+S-12h, and 3d+S-0h vs. 3d+S-24h, seven candidate pathways were found. From the combined metabolome and transcriptome results of enriched KEGG pathways, phenylpropanoid biosynthesis (ko00940) and flavonoid biosynthesis (ko00941) responded to salt stress by accumulating flavonoids.

Table 4

GroupsTotal DEGsUp DEGsDown DEGs
3d+S-0h vs 3d+S-12h1230812418
3d+S-0h vs 3d+S-24h1450904546

The number of differentially expressed genes (DEGs).

Table 5

No.Ko IDKEGG pathwayP-value
1Ko00195Photosynthesis0.000
2Ko00196Photosynthesis—antenna proteins0.000
3Ko00904Diterpenoid biosynthesis0.000
4Ko00360Phenylalanine metabolism0.000
5Ko00940Phenylpropanoid biosynthesis0.000
6Ko00710Carbon fixation in photosynthetic organisms0.001
7Ko00730Thiamine metabolism0.008
8Ko00330Arginine and proline metabolism0.009
9Ko00860Porphyrin and chlorophyll metabolism0.012
10Ko00380Tryptophan metabolism0.015
11Ko00290Valine, leucine and isoleucine biosynthesis0.016
12Ko00908Zeatin biosynthesis0.022
13Ko00906Carotenoid biosynthesis0.026
14Ko00280Valine, leucine and isoleucine degradation0.030
15Ko01220Degradation of aromatic compounds0.033
16Ko00591Linoleic acid metabolism0.047

The enriched pathways of differentially expressed genes (DEGs) compared 3d+S-0h vs. 3d+S-12h in KEGG analysis.

Table 6

No.Ko IDKEGG pathwayP-value
1Ko00195Photosynthesis0.000
2Ko04712Circadian rhythm—plant0.000
3Ko00196Photosynthesis—antenna proteins0.000
4Ko04075Plant hormone signal transduction0.000
5Ko00360Phenylalanine metabolism0.000
6Ko00941Flavonoid biosynthesis0.001
7Ko00904Diterpenoid biosynthesis0.001
8Ko00940Phenylpropanoid biosynthesis0.004
9Ko00906Carotenoid biosynthesis0.008
10Ko00710Carbon fixation in photosynthetic organisms0.014
11Ko00905Brassinosteroid biosynthesis0.043
12Ko00350Tyrosine metabolism0.044
13Ko02010ABC transporters0.046

The enriched pathways of differentially expressed genes (DEGs) compared 3d+S-0h vs. 3d+S-24h in KEGG analysis.

Analysis and verification of pathway mechanisms

The mechanisms of pathways enriched under salt stress are presented in Figure 5. Phenylpropanoid biosynthesis was the upstream pathway in which the metabolites, such as caffeate, scopoletin, sinapic acid, and sinapyl alcohol, were differentially expressed and directly affected the downstream pathway (flavonoid biosynthesis). Also, DEGs were enriched in some enzymes, including 4-coumarate-CoA ligase (4CL), caffeoyl-CoA o-methyltransferase, CAD, and POD. The flavonoid biosynthetic pathway was also enriched as the downstream pathway of phenylpropanoid biosynthesis, and the DEGs were enriched in some enzymes, such as FLS, CHI, and dihydroflavonol-4-reductase. A heatmap of DAMs (such as rutin, cosmosiin, luteolin, garbanzol, and quercetin) showed a significant change in the flavonoid biosynthesis pathway. Analysis of these two pathways revealed that the phenylpropanoid and flavonoid biosynthesis pathways could be crucial in salt stress response by promoting flavonoid accumulation.

Figure 5

Herein, enzyme activity analysis in DEG-enriched points and qRT-PCR analysis was also done. Four enzyme activities that were DEG-enriched were tested. The results showed that salt stress significantly changed the activities of CAD, POD, CHI, and FLS (Figure 6). Besides, the DEGs heatmap showed that DEGs' expression associated with the enzymes differed significantly following salt exposure. Also, six different DEGs in these two pathways, which were enriched in CAD (Phvul.002G144800 and Phvul.003G029500), POD (Phvul.001G143300, Phvul.007G008400 and Phvul.008G249900), and FLS (Phvul.003G216600), were selected for qRT-PCR analysis while their expression and functional analysis are shown in Supplementary Table 9. The results showed that the values of |log2(FoldChange)| expression were all >1 (Figure 7), suggesting that these six genes had a significant expression change at the transcriptional level. Metabolome and transcriptome results revealed that phenylpropanoid and flavonoid biosynthesis responded to salt stress by promoting flavonoid accumulation.

Figure 6

Figure 7

Discussion

Flavonoid trends in germinating common beans

Seed germination is the beginning of most plant life cycles. When seeds (including legumes) begin to germinate, they break dormancy and restore physiological function, which induces metabolism and produces secondary metabolites and other nutrients (37). Rutin is often used as a flavonoid standard in botanical research (26, 38–40) because of its high content in legumes (41, 42). Similarly, the metabolome analysis conducted in this study revealed that germinating common beans have high levels of rutin. Here, the rutin content in germinating common beans initially increased, which illustrates that the beans were activated from hibernation and produced flavonoids (43). The rutin content then decreased after 96 h, which might be related to tissue differentiation after plant germination (44). Other nutrients in germinating beans also show similar trends; for instance, the γ-aminobutyric acid (GABA) content in germinating adzuki beans (Vigna angularis) showed increasing and decreasing, of which the highest content was at 48 h (37). The polyphenol content in germinating mung beans (V. radiata) also showed similar trends, and day 3 was the best time because it had the highest polyphenol content (45). Similarly, Lu and Guo (46) demonstrated that the vitamin C content of beans was at the highest level on day 3 and also showed an increasing and decreasing trend. Flavonoids (represented by rutin in this study) are important intermediates in plants' physiological metabolism, and their content also changes significantly during seed germination (47). Although secondary metabolites show these trends, the time elapsed for the levels to peak varies among species. This phenomenon might be because different nutrient types have different accumulation mechanisms, while different genotypes, or even different legume species, could also affect optimal times. In this study, rutin levels peaked at 96 h, which could be used as the optimal germination time.

Short-term salt stress promotes rutin accumulation

Abiotic stress could influence grain quality and composition in crops (48) and significantly affect the flavonoid content in plants because flavonoids are an important regulator of plant responses to abiotic stress (49). As common abiotic stress, salt stress is important and accelerates the accumulation of certain plant metabolites (50, 51). Several studies have demonstrated that 70 mmol/L NaCl is the appropriate salt concentration for inducing salinity stress in the common bean (22, 52). Rutin is an important class of secondary metabolites widely found in plants and contributes to plant growth and development (16). Previous studies demonstrated that rutin content increases under stress, especially abiotic stress (42, 52, 53). Numerous studies found that flavonoids, including rutin, significantly accumulate under abiotic stress in many plant species, including Arabidopsis, maize (54), tomato (55), and green tea (56). However, flavonoids, including rutin, are stressful substances, and prolonged stress does not necessarily enable continued flavonoid accumulation (57). The antioxidant levels in soybean during drought showed an alternating increasing and decreasing trend, suggesting that short-term drought could rapidly enhance antioxidant intensities, while long-term droughts might cause plants to gradually adapt to these stresses and stop enhancement (58). Also, plant antioxidant enzyme activities showed a similar trend under abiotic stress, which revealed that short-term stress could increase plants' secondary metabolite production, while secondary metabolites under long-term stress might return to a lower level (59). In contrast, stress might inhibit plant germination and reduce flavonoid accumulation. The salt-tolerant common bean germinated better than the salt-sensitive variety, while oxidase enzyme activities of salt-sensitive variety produced insignificant differences at the sprouting stage, limiting plant growth and vitality (52). Here, short-term salt stress was more conducive to rutin accumulation in germinating common beans, and 3d+S could be used to study the mechanism of rutin accumulation under salt stress.

The mechanism of flavonoid (rutin) accumulation under salt stress

In this study, phenylpropanoid and flavonoid biosynthesis were the two pathways that responded to salt stress and promoted flavonoid accumulation. Interestingly, phenylpropanoid biosynthesis was studied by a botanist researching the effects of abiotic stress on plant growth (60). Furthermore, phenylpropanoid biosynthesis is the pathway by which most flavonoids are produced, while enzyme activities in the phenylpropanoid biosynthetic pathway could alter the flavonoid accumulation rate (16). A key enzyme in phenylpropanoid metabolism is 4CL, which responds to abiotic stress, and the 4CL genes have a function in regulating resilience (61). CAD is the final enzyme-catalyzed step in the phenylpropanoid biosynthesis pathway, and CAD activity is affected by salt stress, thus affecting the flavonoid accumulation rate (62). POD is a large family of plant-specific oxidoreductases, and POD activity can change significantly under salt stress (63).

As the downstream pathway of phenylpropanoid biosynthesis, flavonoid biosynthesis is an important pathway for flavonoid accumulation (16). Some enzymes in the flavonoid metabolic pathway, such as CHI and FLS, have also been implicated in flavonoid accumulation and abiotic stress response (64). CHI catalyzes the conversion of bicyclic chalcone into tricyclic (2S)-flavanone and is linked to salt stress (65). FLS mediates the oxidation of dihydroflavonol to flavonol, while FLS gene upregulated and increased flavonoid content under salt stress (66). In this study, enzyme activities and expression levels were determined. The results revealed that salt stress affected the expression level of genes in the phenylpropanoid and flavonoid biosynthesis pathways and influenced the activities of CAD, POD, CHI, and FLS, which might be linked to flavonoid accumulation (Figure 8).

Figure 8

Conclusion

In this study, rutin accumulation trends in germinating common beans were researched, and 96 h was the most suitable germination time because the sprouts had the highest rutin content. Sprouts (including epicotyls and hypocotyls) had more rutin content than cotyledons, hence they could be used as a target research tissue. Short-term salt treatment was more favorable for flavonoid accumulation than long-term salt treatment. Furthermore, metabolome and transcriptome analysis during short-term salt treatment revealed that phenylpropanoid and flavonoid biosynthesis are the enriched pathways that respond to salt stress and promote flavonoid accumulation. This study provides a rationale for germinating common beans and gives a new insight into the molecular mechanism of flavonoid accumulation in common beans under salt stress.

Funding

This study was financially supported by the National Key Research and Development Program (2020YFD1001402), the Research Project of Heilongjiang Bayi Agricultural University (XDB2011-02), and Special Funds from the Central Finance to Support the Development of Local Universities. The funding organizations had no role in experimental design, data collection, analysis, and interpretation of data or writing of the manuscript.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

Data availability statement

The transcriptome raw data was deposited to NCBI database, accession number PRJNA746732. The metabolome raw data was deposited to Zenodo database: doi.org/10.5281/zenodo.6820497.

Author contributions

QZ: data curation and writing original draft. HS, GZ, JixZ, and ZZ: data curation. WZ, QW, and JinZ: conceptualization and methodology. JX and YG: software. Y-lD: formal data analysis and preparation of materials. QZ, ZY, and Y-lD: methodology and revised the manuscript. J-dD: data curation and funding acquisition. All authors contributed to the article and approved the submitted version.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2022.928805/full#supplementary-material

Supplementary Figure 1

The standard reference rutin curve.

Supplementary Figure 2

The proof of metabolome test.

Supplementary Table 1

Detailed information on the different treatments.

Supplementary Table 2

The qRT-PCR primers for differentially expressed genes (DEGs).

Supplementary Table 3

The detailed information of the metabolome analysis in this study.

Supplementary Table 4

The detailed information of differentially altered metabolites (DAMs) in the comparison of 3d+S-0h vs. 3d+S-24h.

Supplementary Table 5

The detailed information of differentially altered metabolites (DAMs) in the comparison of 3d+S-0h vs. 3d+S-12h.

Supplementary Table 6

The order of each compound expression in the five flavonoid kinds analyzed by metabolomics of differentially altered metabolites (DAMs) in the comparison of 3d+S-0h vs. 3d+S-12h or 3d+S-0h vs. 3d+S-24h.

Supplementary Table 7

The quality of the transcriptome results.

Supplementary Table 8

Transcriptome information uploaded to the NCBI database.

Supplementary Table 9

The expression and functional analysis of differentially expressed genes (DEGs) in qRT-PCR analysis.

References

  • 1.

    GanesanKXuB. Polyphenol-rich dry common beans (phaseolus vulgaris l.) and their health benefits. Int J Mol Sci. (2017) 18:2331. 10.3390/ijms181123313

  • 2.

    McDougallGJ. Phenolic-enriched foods: sources and processing for enhanced health benefits. Proc Nutr Soc. (2017) 2:163–71. 10.1017/S0029665116000835

  • 3.

    Díaz-BatallaLWidholmJMFaheyGCCastaño-TostadoEParedes-LópezO. Chemical components with health implications in wild and cultivated mexican common bean seeds (Phaseolus Vulgaris L). J Agric Food Chem. (2006) 6:2045–52. 10.1021/jf051706l

  • 4.

    SajidMChannakesavulaCNStoneSRKaurP. Synthetic Biology Towards Improved Flavonoid Pharmacokinetics. Biomolecules. (2021) 5. 10.3390/biom11050754

  • 5.

    SudhakaranMSardesaiSDoseffAI. Flavonoids: new frontier for immuno-regulation and breast cancer control. Antioxidants (Basel). (2019) 4. 10.3390/antiox8040103

  • 6.

    YeQLiuKShenQLiQHaoJHanFet al. Reversal of multidrug resistance in cancer by multi-functional flavonoids. Front Oncol. (2019) 487. 10.3389/fonc.2019.00487

  • 7.

    AmawiHAshbyCRTiwariAK. Cancer Chemoprevention through Dietary Flavonoids: what's Limiting?Chin J Cancer. (2017) 1:50. 10.1186/s40880-017-0217-4

  • 8.

    NabaviSMÅ amecDTomczykMMilellaLRussoDHabtemariamSet al. Flavonoid biosynthetic pathways in plants: versatile targets for metabolic engineering. Biotechnol Adv. (2020) 38:107316. 10.1016/j.biotechadv.2018.11.005

  • 9.

    Winkel-ShirleyB. Flavonoid biosynthesis. A colorful model for genetics, biochemistry, cell biology, and biotechnology. Plant Physiol. (2001) 2:485–93. 10.1104/pp.126.2.485

  • 10.

    SasakiNNakayamaT. Achievements and perspectives in biochemistry concerning anthocyanin modification for blue flower coloration. Plant Cell Physiol. (2015) 1:28–40. 10.1093/pcp/pcu097

  • 11.

    NakamuraKKoyamaMIshidaRKitaharaTNakajimaTAoyamaT. Characterization of bioactive agents in five types of marketed sprouts and comparison of their antihypertensive, antihyperlipidemic, and antidiabetic effects in fructose-loaded SHRs. J Food Sci Technol. (2016) 53:581–90. 10.1007/s13197-015-2048-0

  • 12.

    SritongtaeBSangsukiamTMorganMRDuangmalK. Effect of acid pretreatment and the germination period on the composition and antioxidant activity of rice bean (Vigna umbellata). Food Chem. (2017) 227:280–8. 10.1016/j.foodchem.2017.01.103

  • 13.

    JohnsonJBSkylasDJManiJSXiangJWalshKBNaikerM. Phenolic profiles of ten australian faba bean varieties. Molecules. (2021) 26:4642. 10.3390/molecules26154642

  • 14.

    TsurunagaYTakahashiTKatsubeTKudoAKuramitsuOIshiwataMMatsumotoS. Effects of UV-B irradiation on the levels of anthocyanin, rutin and radical scavenging activity of buckwheat sprouts. Food Chem. (2013) 141:552–6. 10.1016/j.foodchem.2013.03.032

  • 15.

    TaylorLPGrotewoldE. Flavonoids as developmental regulators. Curr Opin Plant Biol. (2005) 3:317–23. 10.1016/j.pbi.2005.03.005

  • 16.

    LiuWFengYYuSFanZLiXLiJet al. The flavonoid biosynthesis network in plants. Int J Mol Sci. (2021) 23. 10.3390/ijms222312824

  • 17.

    LiQYuHMMengXFLinJSLiYJHouBKet al. Ectopic expression of glycosyltransferase Ugt76e11 increases flavonoid accumulation and enhances abiotic stress tolerance in arabidopsis. Plant Biol (Stuttg). (2018) 1:10–9. 10.1111/plb.12627

  • 18.

    TlahigSBellaniLKarmousIBarbieriFLoumeremMMucciforaSet al. Response to salinity in legume species: an insight on the effects of salt stress during seed germination and seedling growth. Chem Biodivers. (2021) 4:e2000917. 10.1002/cbdv.202000917

  • 19.

    EumHLParkYYiTGLeeJWHaKSChoiIYet al. Effect of germination environment on the biochemical compounds and anti-inflammatory properties of soybean cultivars. PLoS ONE. (2020) 4:e0232159. 10.1371/journal.pone.0232159

  • 20.

    GuoXLiTTangKLiuRH. Effect of Germination on Phytochemical Profiles and Antioxidant Activity of Mung Bean Sprouts (Vigna Radiata). J Agric Food Chem. (2012) 44:11050–5. 10.1021/jf304443u

  • 21.

    GuHWangYXieHQiuCZhangSXiaoJet al. Drought stress triggers proteomic changes involving lignin, flavonoids and fatty acids in tea plants. Sci Rep. (2020) 1:15504. 10.1038/s41598-020-72596-1

  • 22.

    ZhangQGengJDuYZhaoQZhangWFangQet al. Heat shock transcription factor (hsf) gene family in common bean (phaseolus vulgaris): genome-wide identification, phylogeny, evolutionary expansion and expression analyses at the sprout stage under abiotic stress. BMC Plant Biol. (2022) 1:33. 10.1186/s12870-021-03417-4

  • 23.

    LiuJLiCDingGQuanW. Artificial intelligence assisted ultrasonic extraction of total flavonoids from rosa sterilis. Molecules. (2021) 26:3835. 10.3390/molecules26133835

  • 24.

    LouHYuanHYamazakiYSasakiTOkaS. Alkaloids and flavonoids from peanut skins. Planta Med. (2001) 67:345–9. 10.1055/s-2001-14319

  • 25.

    RolimAMacielCPKanekoTMConsiglieriVOSalgado-SantosIMVelascoMV. Validation assay for total flavonoids, as rutin equivalents, from Trichilia catigua Adr. Juss (Meliaceae) and Ptychopetalum olacoides Bentham (Olacaceae) commercial extract. J AOAC Int. (2005) 88:1015–9.

  • 26.

    RolimAOishiTMacielCPZagueVPintoCAKanekoTMet al. Total flavonoids quantification from O/W emulsion with extract of Brazilian plants. Int J Pharm. (2006) 308:107–14. 10.1016/j.ijpharm.2005.10.031

  • 27.

    ChenWGongLGuoZWangWZhangHLiuXet al. A novel integrated method for large-scale detection, identification, and quantification of widely targeted metabolites: application in the study of rice metabolomics. Mol Plant. (2013) 6:1769–80. 10.1093/mp/sst080

  • 28.

    ScholzMSelbigJ. Visualization and analysis of molecular data. Methods Mol Biol. (2007) 87–104. 10.1007/978-1-59745-244-1_6

  • 29.

    WixonJKellD. The Kyoto Encyclopedia of Genes and Genomes—Kegg. Yeast. (2000) 1:48–55. 10.1155/2000/981362

  • 30.

    KimDPerteaGTrapnellCPimentelHKelleyRSalzbergSLet al. Tophat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. (2013) 4:R36. 10.1186/gb-2013-14-4-r36

  • 31.

    FloreaLSongLSalzbergSL. Thousands of exon skipping events differentiate among splicing patterns in sixteen human tissues. (2013) F1000Res: 188. 10.12688/f1000research.2-188.v1

  • 32.

    GoffnerDJoffroyIGrima-PettenatiJHalpinCKnightMESchuchWet al. Purification and characterization of isoforms of cinnamyl alcohol dehydrogenase from eucalyptus xylem. Planta. (1992) 1:48–53. 10.1007/BF01160711

  • 33.

    DoergeDRDiviRLChurchwellMI. Identification of the colored guaiacol oxidation product produced by peroxidases. Anal Biochem. (1997) 1:10–7. 10.1006/abio.1997.2191

  • 34.

    ListerCELancasterJEWalkerJR. Developmental changes in enzymes of flavonoid biosynthesis in the skins of red and green apple cultivars. J Sci Food Agric. (1996). 10.1002/(SICI)1097-0010(199607)71:3&lt;313::AID-JSFA586&gt;3.0.CO;2-N

  • 35.

    HrazdinaGWeedenNF. Enzymatic control of anthocyanin expression in the flowers of pea (Pisum Sativum) mutants. Biochem Genet. (1986) 3–4:309–17. 10.1007/BF00502797

  • 36.

    LivakKJSchmittgenTD. Analysis of relative gene expression data using real-time quantitative pcr and the 2(-Delta Delta C(T)) Method. Methods. (2001) 4:402–8. 10.1006/meth.2001.1262

  • 37.

    JiangXXuQZhangALiuYZhaoLGuLet al. Optimization of Γ-aminobutyric acid (gaba) accumulation in germinating adzuki beans (vigna angularis) by vacuum treatment and monosodium glutamate, and the molecular mechanisms. Front Nutr. (2021) 693862. 10.3389/fnut.2021.693862

  • 38.

    Rice-EvansCAMillerNJPagangaG. Structure-antioxidant activity relationships of flavonoids and phenolic acids. Free Radic Biol Med. (1996) 20:933–56. 10.1016/0891-5849(95)02227-9

  • 39.

    JanjuaNKSiddiqaAYaqubASabahatSQureshiRul HaqueS. Spectrophotometric analysis of flavonoid-DNA binding interactions at physiological conditions. Spectrochim Acta A Mol Biomol Spectrosc. (2009) 74:1135–7. 10.1016/j.saa.2009.09.0225

  • 40.

    YangYZhaoYZuoXWangY. Determination of total flavonoids for paris polyphylla var. yunnanensis in different geographical origins using uv and ft-ir spectroscopy. J AOAC Int. (2019) 2:457–64. 10.5740/jaoacint.18-0188

  • 41.

    KanLNieSHuJWangSCuiSWLiYet al. Nutrients, phytochemicals and antioxidant activities of 26 kidney bean cultivars. Food Chem Toxicol. (2017) PtB:47–77. 10.1016/j.fct.2016.09.007

  • 42.

    XuCWeiLHuangSYangCWangYYuanHet al. Drought resistance in qingke involves a reprogramming of the phenylpropanoid pathway and udp-glucosyltransferase regulation of abiotic stress tolerance targeting flavonoid biosynthesis. J Agric Food Chem. (2021) 13:3992–4005. 10.1021/acs.jafc.0c07810

  • 43.

    KayembeNJansenV. Germination as a processing technique for soybeans in small-scale farming. S Afr J Anim Sci. (2013) 2:167–73. 10.4314/sajas.v43i2.7

  • 44.

    YangRYinYGuoLHanYGuZ. Sequence analysis of diamine oxidase gene from fava bean and its expression related to Γ-aminobutyric acid accumulation in seeds germinating under hypoxia-nacl stress. J Sci Food Agric. (2014) 8:1585–91. 10.1002/jsfa.6461

  • 45.

    KapravelouGMartínezRPerazzoliGSánchez GonzálezCLlopisJCantareroSet al. Germination improves the polyphenolic profile and functional value of mung bean (Vigna radiata L.). Antioxidants (Basel). (2020) 9:746. 10.3390/antiox9080746

  • 46.

    LuYGuoX. The effect of light in vitamin c metabolism regulation and accumulation in mung bean (Vigna Radiata) germination. Plant Foods Hum Nutr. (2020) 1:24–9. 10.1007/s11130-019-00787-x

  • 47.

    YimingZHongWLinlinCXiaoliZWenTXinliSet al. Evolution of nutrient ingredients in tartary buckwheat seeds during germination. Food Chem. (2015) 244–8. 10.1016/j.foodchem.2015.03.115

  • 48.

    FarooqMHussainMUsmanMFarooqSAlghamdiSSSiddiqueKHMet al. Impact of abiotic stresses on grain composition and quality in food legumes. J Agric Food Chem. (2018) 34:8887–97. 10.1021/acs.jafc.8b02924

  • 49.

    NakabayashiRSaitoK. Integrated metabolomics for abiotic stress responses in plants. Curr Opin Plant Biol. (2015) 10–6. 10.1016/j.pbi.2015.01.003

  • 50.

    BorsaiOHassanMANegruierCRaigónMDBoscaiuMSestraREet al. Responses to salt stress in portulaca: Insight into its tolerance mechanisms. Plants (Basel). (2020) 9:1660. 10.3390/plants9121660

  • 51.

    GenzelFDickeMDJunker-FrohnLVNeuwohnerAThieleBPutzAet al. Impact of moderate cold and salt stress on the accumulation of antioxidant flavonoids in the leaves of two capsicum cultivars. J Agric Food Chem. (2021) 69:6431–43. 10.1021/acs.jafc.1c00908

  • 52.

    ZhangQLiMXiaCYZhangWJYinZGZhangYLet al. Transcriptome-based analysis of salt-related genes during the sprout stage of common bean (Phaseolus Vulgaris) under salt stress conditions. Biotechnol Biotechnol Equip. (2021) 1:1086–98. 10.1080/13102818.2021.1954091

  • 53.

    Röhlen-SchmittgenSEllenbergerJGroherTHunscheM. Boosting leaf contents of rutin and solanesol in bio-waste of Solanum lycopersicum. Plant Physiol Biochem. (2020) 155:888–97. 10.1016/j.plaphy.2020.08.035

  • 54.

    QinLSunLWeiLYuanJKongFZhangYet al. Maize Sro1e represses anthocyanin synthesis through regulating the mbw complex in response to abiotic stress. Plant J. (2021) 4:1010–25. 10.1111/tpj.15083

  • 55.

    AtkinsonNJDewTPOrfilaCUrwinPE. Influence of combined biotic and abiotic stress on nutritional quality parameters in tomato (Solanum lycopersicum). J Agric Food Chem. (2011) 17:9673–82. 10.1021/jf202081t

  • 56.

    PrasanthMISivamaruthiBSChaiyasutCTencomnaoT. A review of the role of green tea (Camellia sinensis) in antiphotoaging, stress resistance, neuroprotection, and autophagy. Nutrients. (2019) 11:474. 10.3390/nu11020474

  • 57.

    NakabayashiRMoriTSaitoK. Alternation of flavonoid accumulation under drought stress in arabidopsis thaliana. Plant Signal Behav. (2014) 8:e29518. 10.4161/psb.29518

  • 58.

    Easwar RaoDViswanatha ChaitanyaK. Changes in the Antioxidant Intensities of Seven Different Soybean (Glycine Max (L) Merr) Cultivars During Drought. J Food Biochem. (2020) 2:e.(13118). 10.1111/jfbc.13118

  • 59.

    DuYZhaoQChenLYaoXZhangHWuJet al. Effect of drought stress during soybean R2-R6 growth stages on sucrose metabolism in leaf and seed. Int J Mol Sci. (2020) 21:618. 10.3390/ijms21020618

  • 60.

    CaoLJinXZhangYZhangMWangY. Transcriptomic and metabolomic profiling of melatonin treated soybean (Glycine Max L.) under drought stress during grain filling period through regulation of secondary metabolite biosynthesis pathways. PLoS ONE. (2020) 10:e0239701. 10.1371/journal.pone.0239701

  • 61.

    KhakdanFNasiriJRanjbarMAlizadehH. Water deficit stress fluctuates expression profiles of 4cl, C3h, comt, cvomt and eomt genes involved in the biosynthetic pathway of volatile phenylpropanoids alongside accumulation of methylchavicol and methyleugenol in different iranian cultivars of Basil. J Plant Physiol. (2017) 74–83. 10.1016/j.jplph.2017.07.012

  • 62.

    MouraJCBonineCAde Oliveira Fernandes VianaJDornelasMCMazzaferaP. Abiotic and biotic stresses and changes in the lignin content and composition in plants. J Integr Plant Biol. (2010) 4:360–76. 10.1111/j.1744-7909.2010.00892.x

  • 63.

    ZhuJK. Salt and drought stress signal transduction in plants. Annu Rev Plant Biol. (2002) 247–73. 10.1146/annurev.arplant.53.091401.143329

  • 64.

    ZukMSzperlikJHniteckaASzopaJ. Temporal biosynthesis of flavone constituents in flax growth stages. Plant Physiol Biochem. (2019) 234–45. 10.1016/j.plaphy.2019.07.009

  • 65.

    WangHHuTHuangJLuXHuangBZhengYet al. The expression of millettia pinnata chalcone isomerase in saccharomyces cerevisiae salt-sensitive mutants enhances salt-tolerance. Int J Mol Sci. (2013) 5:8775–86. 10.3390/ijms14058775

  • 66.

    XuNLiuSLuZPangSWangLWangLet al. Gene expression profiles and flavonoid accumulation during salt stress in ginkgo biloba seedlings. Plants (Basel). (2020) 9:1162. 10.3390/plants9091162

Summary

Keywords

flavonoids (rutin) content, common bean, germinating, salt stress, mechanism, phenylpropanoid biosynthesis, flavonoid biosynthesis

Citation

Zhang Q, Zheng G, Wang Q, Zhu J, Zhou Z, Zhou W, Xu J, Sun H, Zhong J, Gu Y, Yin Z, Du Y and Du J (2022) Molecular mechanisms of flavonoid accumulation in germinating common bean (Phaseolus vulgaris) under salt stress. Front. Nutr. 9:928805. doi: 10.3389/fnut.2022.928805

Received

26 April 2022

Accepted

22 July 2022

Published

29 August 2022

Volume

9 - 2022

Edited by

Alberto Valdés, Spanish National Research Council (CSIC), Spain

Reviewed by

Carlos L. Cespedes-Acuña, University of Bio Bio Chillan Chile, Chile; Ya-Fang Shang, Hefei University of Technology, China

Updates

Copyright

*Correspondence: Ji-dao Du

This article was submitted to Nutrition and Food Science Technology, a section of the journal Frontiers in Nutrition

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics