Skip to main content

ORIGINAL RESEARCH article

Front. Plant Sci., 17 October 2022
Sec. Plant Metabolism and Chemodiversity
This article is part of the Research Topic Biosynthesis and Regulation of Plant Specialized Metabolisms View all 12 articles

Metabolomic and transcriptomic changes in mungbean (Vigna radiata (L.) R. Wilczek) sprouts under salinity stress

  • Department of Plant Science, Gangneung-Wonju National University, Gangneung, South Korea

Mungbean (Vigna radiata) sprouts are consumed globally as a healthy food with high nutritional values, having antioxidant and anticancer capacity. Under mild salinity stress, plants accumulate more secondary metabolites to alleviate oxidative stress. In this study, metabolomic and transcriptomic changes in mungbean sprouts were identified using a reference cultivar, sunhwa, to understand the regulatory mechanisms of secondary metabolites in response to salinity stress. Under salinity conditions, the contents of phenylpropanoid-derived metabolites, including catechin, chlorogenic acid, isovitexin, p-coumaric acid, syringic acid, ferulic acid, and vitexin, significantly increased. Through RNA sequencing, 728 differentially expressed genes (DEGs) were identified and 20 DEGs were detected in phenylpropanoid and flavonoid biosynthetic pathways. Among them, 11 DEGs encoding key enzymes involved in the biosynthesis of the secondary metabolites that increased after NaCl treatment were significantly upregulated, including dihydroflavonol 4-reductase (log2FC 1.46), caffeoyl-CoA O-methyltransferase (1.38), chalcone synthase (1.15), and chalcone isomerase (1.19). Transcription factor families, such as MYB, WRKY, and bHLH, were also identified as upregulated DEGs, which play a crucial role in stress responses in plants. Furthermore, this study showed that mild salinity stress can increase the contents of phenylpropanoids and flavonoids in mungbean sprouts through transcriptional regulation of the key enzymes involved in the biosynthetic pathways. Overall, these findings will provide valuable information for molecular breeders and scientists interested in improving the nutritional quality of sprout vegetables.

Introduction

Mungbean (Vigna radiata (L.) R. Wilczek) is an important legume crop grown over 6 million hectares worldwide (about 8.5% of the global pulse cultivation area) (Gayacharan et al., 2020). Mature mungbean seeds are valuable food sources of protein and starch in many developing countries in Asia (Nair et al., 2013). Mungbean sprouts have long been consumed globally because of their excellent texture and abundant nutrients, including ascorbic acid, dietary fiber, essential amino acids, and vitamins (Gan et al., 2016; Ebert et al., 2017; Singh et al., 2017). In particular, six polyphenolic compounds, including caffeic acid, ferulic acid, gallic acid, p-coumaric acid, catechin, and rutin, have been identified in mungbean sprouts (Gan et al., 2016), which help maintain human health and prevent chronic diseases (Kim et al., 2012; Nderitu et al., 2013; Pudenz et al., 2014; Stagos, 2020). Owing to these nutritional values, mungbean sprout consumption has increased, and the quality of sprout vegetables has attracted consumers’ attention globally.

Plants accumulate polyphenolic compounds as defense metabolites against harmful environments, and this accumulation varies depending on different stress exposures, such as high/low temperature, salinity, ultraviolet, and drought (Ksouri et al., 2008; Yang et al., 2018; Sharma et al., 2019; Chrysargyris et al., 2020). Under salinity conditions, osmotic, ionic, and oxidative stresses affect most plants, inducing morphological and biochemical changes (Yang and Guo, 2018). While severe salinity disrupts general plant development through ionic toxicity and water imbalance (Yang and Guo, 2018), mild NaCl treatment can considerably enhance the polyphenol contents and antioxidant activity in many plant species, such as rice seedlings as well as mungbean, broccoli, buckwheat, and radish sprouts (Yuan et al., 2010; Lim et al., 2012; Guo et al., 2014; Koodkaew, 2019). Similarly, many metabolic studies have shown that anthocyanin and flavonoid metabolisms are significantly influenced by mild NaCl treatment, followed by enhanced production of catechin, chlorogenic acid, ferulic acid, p-coumaric acid, syringic acid, and vitexin (Lim et al., 2012; Hassini et al., 2017; Sarker and Oba, 2018; Banik and Bhattacharjee, 2020; Benincasa et al., 2021). Therefore, salinity stress can play a crucial role in regulating the quality of the nutritional values of plant ingredients.

With RNA-seq-based transcriptome analysis, the gene expression patterns and regulation network in response to salinity have been widely investigated in several plants, including cotton (Wang et al., 2012), rice (Chandran et al., 2019), buckwheat (Ma et al., 2019), soybean (Fan et al., 2013), and medicine plants (Ben Abdallah et al., 2016). The expression of the genes involved in the biosynthetic pathways of flavonoids and phenylpropanoids can be increased by NaCl treatment in solanum nigrum, including phenylalanine ammonialyase, chalcone synthase, and flavonol synthase (Ben Abdallah et al., 2016). Furthermore, various transcription factor (TF) families play a crucial role in salinity stress responses in plants, such as NACs, MYBs, WRKYs, and bHLHs (Eulgem and Somssich, 2007). Overexpression of TaNAC29, AtMYB20, FtMYB10, and OsMYB3R-2 enhanced salt tolerance in transgenic Arabidopsis (Dai et al., 2007, 3; Cui et al., 2013, 3; Xu et al., 2015; Gao et al., 2016). TaMYB56-B, MdSIMYB1, and OsNAC6 were associated with salinity stress responses in wheat, apple, and rice, respectively (Nakashima et al., 2007; Shen et al., 2017).

During the mungbean germination stage, polyphenol accumulation is affected by growth conditions and stresses, such as temperature, salinity, and drought (Kim et al., 2009b). When mungbean sprouts are exposed to NaCl, the total phenolic and flavonoid contents increase, resulting in enhanced antioxidant activities (Koodkaew, 2019). However, only a few transcriptome studies have explored resistance to fungi and drought stress in mungbean seeds/leaves (Kumar et al., 2020), and little is known about the genetic regulation of secondary metabolites for abiotic stresses in mungbean sprouts.

This study’s aim is to investigate the effects of NaCl treatment on morphological, metabolomic, and transcriptomic changes in mungbean sprouts using two cultivars, dahyeon and sunhwa. Dahyeon has been reported to accumulate relatively higher polyphenols in its sprouts compared to other cultivars and sunhwa is a mungbean reference cultivar (Kim et al., 2009a; Habibzadeh and Yagoob, 2014; Kang et al., 2014; Yoseph Ganta et al., 2021). We qualified and quantified the polyphenolic compounds in mungbean sprouts after NaCl treatment and identified the key enzymes involved in the biosynthetic pathways of these compounds. Our findings provide further insights into the molecular mechanisms associated with salinity stress and can help improve the nutritional quality and morphological characteristics of sprout vegetables through genetic engineering in legume crops.

Materials and methods

Plant materials and NaCl treatment

Two mungbean cultivars, dahyeon and sunhwa (VC1973A), were used in our experiments. Seeds were harvested at the Gangneung-Wonju National University Experimental Farm in Gangneung, South Korea (37.77°N, 128.86°E). The seeds were soaked in distilled water for 17 h at 37°C using an incubator (Jeiotech, ISS-4075R, Korea) for germination. The germinated seeds were transplanted into a plant growth chamber (Sundotcom, ST001A, Korea) and cultivated for three days at 28 ± 2°C in the dark. Water spraying automatically occurred for 4 min every 2 h (Kim et al., 2021). During cultivation, mungbean sprouts were soaked in different concentrations (0, 50, 100, and 200 mM) of NaCl solution for 30 min every 12 h (KisanBio, MB-S4636, Korea). Mungbean sprouts were harvested on the third day after germination and stored at −70°C for RNA extraction (Lim et al., 2022). Fifty seeds were cultured for each experimental group, 30 sprouts were randomly picked for further analysis, and five sprouts were used for each biological replication of RNA extraction. The lengths of the root and hypocotyl and the thickness of the hypocotyl were measured using ImageJ (Schneider et al., 2012).

Extraction procedure

Mungbean sprouts were fully dried at 70°C for 24 h using an incubator (Gan et al., 2017). The sample was ground into a fine powder using mortar and pestle and 0.05 g powder of each sample was extracted using 70% ethanol (w/v, 1:10) for 24 h in the dark. Mixtures were vortexed and sonicated for 15 min and centrifuged at 13,000 rpm for 10 min. The supernatant was filtered through syringe filters (SMART I LAB, SPF0213-1, USA, 0.22 um) and diluted with 70% ethanol for further analysis.

Determination of total phenol content

The contents of total phenol were measured using Folin-Ciocalteu colorimetric method with slight modification (Singleton and Rossi, 1965). Gallic acid (0, 50, 100, and 200 mg/L) was used as a standard for quantification. 50 µL Folin–Ciocalteu (Sigma-Aldrich®, 47641, USA) was added to each mungbean sprout extraction (100 µL, 10,000 ppm) in an EP tube. After 5 min, 300 µL of 20% Na2CO3 (FUJIFILM, 199-01605, Japan) solution was added. After 15 min in dark condition, distilled water was added up to 1 mL. After 2 min of centrifuge, 200 ul of the supernatants were transferred to a 96-well plate, and absorbance was measured at 740 nm wavelength using a spectrophotometer (Thermo Scientific MIB, Multiskan FC, Korea, 738 nm), with three replications. The results were expressed as milligrams gallic acid equivalent (GAE) per 1 g of dry weight (GAE mg/g).

Determination of total flavonoid content

The contents of total flavonoids were measured using the aluminum nitrate colorimetric method with slight modifications (Zhishen et al., 1999). Quercetin (Sigma-Aldrich®, Q4951, USA) (0, 50, 100, and 200 mg/L) was used as standard. The mungbean sprouts extract (500 µL, 50,000 ppm) was mixed with a 10% aluminum nitrate (100 µL) solution and 1 M potassium acetate solution (100 µL). After 40 min in dark conditions and 2 min of centrifuge, 150 µL of the supernatants were transferred to a 96-well plate, and absorbance was measured at 429 nm wavelength using a spectrophotometer, with three replications. The results were expressed as milligrams quercetin equivalent (QAE) per 1 g of dry weight (QAE mg/g).

Determination of antioxidant activity through DPPH and ABTS assays

DPPH radical scavenging activity was measured using OxiTec™ DPPH Antioxidant Assay Kit (BIOMAX, BO-DPH-500, Korea) according to the manufacturer’s protocol. Trolox (0, 40, 60, 80, and 100 mg/L) was used as the standard. Each extraction (20 µL, 50,000 ppm) was mixed with assay buffer (80 µL) and DPPH radical (100 µL) in 96-well plates, with three replications. The samples were preserved in the dark for 30 min, and absorbance was measured at 520 nm using a spectrophotometer. The result was expressed as a percentage of scavenging achieved by the DPPH.

ABTS radical scavenging activity was measured as described by Yen (Yen and Chen, 1995). Ascorbic acid (FUJIFILM, 012–04802, Japan) (0, 1, 5, 10, 25, 50, and 100 mg/L) was used as the standard. ABTS (Roche, Cat. No. 10102946001, Swiss) 7.4 mM solution and 2.6 mM potassium persulfate (YAKURI, 28718, Korea) were diluted to a 1:1 ratio (v/v). ABTS solution (180 μL) was added to each extraction (20 µL, 1,000 ppm) in 96-well plates, with three replications. The samples were kept in dark for 10 min, their absorbance was measured at 720 nm using a spectrophotometer. The result was expressed as a percentage of scavenging achieved by the ABTS.

Ultra-high-performance liquid chromatography analysis

Ultra-high-performance liquid chromatography (UHPLC) analysis was conducted using the Shimadzu UHPLC system (Nexera series equipped with MPM-40, SCL-40, SPD-M40, LC-40, SIL-40, and CTO-40 unit from Shimadzu, Kyoto, Japan) with a UV-vis detector. Separation was conducted using a C18 column (Shimadzu, Kyoto, Japan; 250 × 4.6 mm, 5 µm). Ultrapure water (Thermo Fisher, W5-4, Korea) with 0.1% ortho-phthalaldehyde buffer and acetonitrile (Thermo Fisher, 022927 M6, Korea) were used as solvents A and B, respectively. The flow rate was 0.2 mL/min and the column temperature was set at 40°C. The ratio of solvent B proceeded with 5%, increased to 30% for 0–8 min, increased to 36% for 8–20 min, and decreased to 5% for 20–25 min. 14 compounds were used as standard compounds: caffeic acid (CFN99190), catechin (CFN99646), chlorogenic acid (CFN99116), gallic acid (CFN99624), isovitexin (CFN98620), kaempferol (CFN98838), myricetin (CFN98877), neochlorogenic acid (CFN97472), p-coumaric acid (CFN97218), quercetin (CFN99272), resveratrol (CFN98791), syringic acid (CFN98884), t-ferulic acid (CFN99158), and vitexin (CFN98601). The standard compounds were purchased from ChemFaces (Wuhan, China). The results were obtained based on standard calibration curves (10-100 ppm) with a tri-repeat and expressed as milligrams per 100 g of dry weight (mg/100g).

cDNA library construction and transcriptome sequencing

Total RNA was isolated from whole mungbean sprouts using a Ribospin™ Plant RNA extraction Kit (GeneAll, Songpa-gu, South Korea) according to the manufacturer’s protocol. The cDNA libraries for RNA-seq were constructed using a TruSeq Stranded mRNA LT Sample Prep Kit (Illumina, Inc., San Diego, CA, USA). Six libraries were constructed from the treatment and control groups, with three biological replications in each group. The size and quality of the libraries used for sequencing were checked using a 2100 Bioanalyzer (Agilent Technologies Inc., Santa Clara, CA, USA). Sequencing runs were conducted in paired-end mode using a Truseq SBS Kit on the Illumina NovaSeq 6000 platform. RNA sequencing data were deposited in the NCBI SRA database (PRJNA853140).

Identification and analysis of differentially expressed genes

Raw reads were filtered out (q-value > 20) using trimmomatic version 0.36 (http://www.usadellab.org/) with default parameters (Bolger et al., 2014). The reference genome and gene annotation data were downloaded from the Seoul National University Crop Genomics Lab (http://plantgenomics.snu.ac.kr) (Ha et al., 2021). The filtered reads were mapped to the mungbean reference genome using Tophat2 (Kim et al., 2013, 2). Counts per million mapped read values were calculated for 30,999 genes using FeatureCounts (Liao et al., 2014). DEGs were counted using Edge-R (Robinson et al., 2010), which were defined as genes with an absolute log2 fold change (FC) value greater or equal to 1 between the two groups (control over treatment).

Functional annotation of DEGs

Gene ontology (GO) enrichment analysis was performed using BINGO software (http://www.psb.ugent.be/cbd/papers/BiNGO/Home.html) (Maere et al., 2005). The Kyoto Encyclopedia of Genes and Genomes (KEGG) ontology (KO) pathway enrichment analysis was conducted using DAVID (https://david.ncifcrf.gov/) (Sherman et al., 2022). The TFs were identified using PlantTFDB (http://planttfdb.gao-lab.org/) (Tian et al., 2020).

Validation by real-time quantitative reverse transcription PCR

Real-time quantitative reverse transcription PCR (qRT-PCR) was performed using a PrimeScript RT reagent Kit with gDNA Eraser (TaKaRa Bio Inc., San Jose, CA, USA) according to the manufacturer’s protocol. Thirty cycles of PCR amplification were conducted using gene-specific primers. Cytochrome P450 family 20 (CYP20) was used as an internal control gene. The primer sequences are listed in the supplementary information (Table S3).

Statistical analysis

Statistical analysis was performed in R (R Core Team, 2017). The significant differences were calculated by a one-way analysis of variance followed by Duncan’s multiple test; p-value <0.05 was considered statistically significant.

Results

Development of mungbean sprouts under NaCl treatments

We first assessed the effect of NaCl treatment on the development of mungbean sprouts using a gradient of different NaCl concentrations (0, 50, 100, and 200 mM) (Figure 1A). The length and thickness of hypocotyl and the length of the roots were measured in dahyeon and sunhwa (Figure 1B). The length of the roots and hypocotyl significantly decreased in both cultivars when the NaCl concentration increased (Figures 1B, C). While hypocotyl thickness showed no significance in dahyeon, it increased significantly in sunhwa under NaCl treatments (Figure 1D).

FIGURE 1
www.frontiersin.org

Figure 1 Morphological changes in mungbean sprouts of dahyeon and sunhwa cultivars under different concentrations of NaCl treatment (0, 50, 100, and 200 mM). (A) Picture of sprouts, (B) length of root, (C) length of hypocotyl, and (D) thickness of hypocotyl. The color bar indicates the mungbean cultivars, blue represents dahyeon and orange represents sunhwa. The results were expressed with a standard error bar, and different lowercase letters indicate statistical significance (P < 0.05).

Dynamic changes in total phenol and flavonoid contents, and antioxidant activities under NaCl treatment

We also evaluated the total phenol and flavonoid contents and antioxidant activity under different NaCl concentrations (Figure 2). In general, for both cultivars, the phenol contents decreased as the NaCl concentration increased (Figure 2A). However, for total flavonoid content, the two cultivars had opposite responses under low concentrations of NaCl treatment (50 mM, Na50). The flavonoid contents of dahyeon showed a significant decrease when NaCl concentrations increased, while those of sunhwa peaked under Na50 and then decreased under higher concentrations of NaCl (100 and 200 mM) (Figure 2B).

FIGURE 2
www.frontiersin.org

Figure 2 Comparison of (A) total phenolic contents, (B) total flavonoid contents, (C) antioxidant scavenging activity of DPPH, and (D) antioxidant scavenging activity of ABTS in the two mungbean cultivars cultivated on different concentrations of NaCl treatment (0, 50, 100, and 200 mM). The color bar indicates the mungbean cultivars, blue represents dahyeon and orange represents sunhwa. The results were expressed with a standard error bar, and different lowercase letters indicate statistical significance (P < 0.05).

DPPH and ABTS assays were conducted to evaluate the antioxidant activity of the mungbean sprouts (Figures 2C, D). In the DPPH assay, no difference in the antioxidant activities was observed in both cultivars treated with 0, 50, and 100 mM NaCl, and the antioxidant activity significantly declined under 200 mM NaCl treatment (Figure 2C). The antioxidant activity detected through the ABTS assay was gradually reduced in both cultivars as NaCl concentration increased (Figure 2D).

Profile of polyphenolic compounds using UHPLC analysis

Metabolic profiling was conducted on sunhwa sprouts treated with Na50, where flavonoid contents significantly increased (Table 1). To better understand the effects of salinity on polyphenolic compounds, we quantified and qualified 14 compounds in mungbean sprouts, including six flavonoids (catechin, isovitexin, kaempferol, myricetin, quercetin, and vitexin) and eight phenolic acids (caffeic acid, chlorogenic acid gallic acid, p-coumaric acid, neochlorogenic acid, resveratrol, syringic acid, and t-ferulic acid) using UHPLC.

TABLE 1
www.frontiersin.org

Table 1 The effects of 50mM of NaCl treatment on the polyphenol composition of mungbean sprouts of sunhwa cultivar.

Among the 14 compounds, 13 compounds excluding kaempferol were detected in both the control and NaCl treatment samples (Table 1). Chlorogenic acid was the most abundant in both the control and NaCl treatments, with values of 1144.1 mg/100g and 1444.2 mg/100g, respectively (Table 1). Myricetin was the second most abundant compound, with values of 328.6 mg/100g (control) and 233.9 mg/100g (Na50), followed by gallic acid (143.3 and 146.8 mg/g), catechin (114.2 and 160.4 mg/100g), quercetin (132.8 and 113.6 mg/100g), and vitexin (75.5 and 123.8 mg/100g), respectively (Table 1).

The contents of seven compounds (catechin, chlorogenic acid, glycerin, isovitexin, p-coumaric acid, syringic acid, t-ferulic acid, and vitexin) were significantly increased after Na50 treatment (Table 1). In contrast, the contents of four compounds (myricetin, neochlorogenic acid, quercetin, and resveratrol) were significantly reduced under NaCl treatment (Table 1). In the other two compounds (caffeic acid and gallic acid), no significant differences were observed in their contents between the control and Na50 (Table 1).

Transcriptome analysis

To investigate the effects of salinity stress on transcriptional changes, we conducted RNA sequencing on sunhwa sprouts treated with Na50. In total, 96.8 Gb of filtered reads (16.1 Gb per library, on average) were obtained from the control and Na50 samples, with three replications (Table S1). On average, 98.08% of the reads were mapped properly against the mungbean reference genome (Ha et al., 2021) (Table S1). DEGs were identified between the NaCl treatment and control groups. A total of 728 DEGs were identified, and 151 and 577 genes were downregulated and upregulated after NaCl treatment, respectively (Figures 3A, B).

FIGURE 3
www.frontiersin.org

Figure 3 Visualization of identified DEGs using (A) bar plot and (B) volcano plot. The y-axis of the volcano plot indicates log2-fold change values of each DEG, and the x-axis indicates –log10 p-value. The color bar indicates gene regulation, blue represents downregulated genes and red represents upregulated genes.

Functional annotation of DEGs through enrichment analysis

The KEGG pathway and GO enrichment analyses were conducted to understand DEG’s molecular functions. In KEGG pathway analysis, 133 DEGs were assigned to 11 pathways and the top three pathways were “phenylpropanoid biosynthesis” (15), “starch and sucrose metabolisms” (11), and “biosynthesis of various plant secondary metabolites” (8) (Figure 4A). In GO analysis, 464 DEGs were clustered in 104 categories, and the major clusters were “response to stimulus” (137), followed by “response to stress” (100) and “response to chemical stimulus” (83) (Figure 4B). Within the categories involved in secondary metabolites, “phenylpropanoid metabolic process” (18) and “flavonoid metabolic process” (6) were identified as enriched categories. Based on the enrichment analyses, 46 DEGs were identified as genes related to the biosynthesis of phenylpropanoid metabolites. The 46 DEGs were mapped into the phenylpropanoid pathway in KEGG, and 11 DEGs were selected as key genes encoding the enzymes directly involved in the biosynthesis of the secondary metabolites that increased significantly after NaCl treatment (Figure 5; Table S2).

FIGURE 4
www.frontiersin.org

Figure 4 (A) KEGG and (B) GO pathway enrichment analyses of DEGs. The y-axis indicates the pathway or term name, and the x-axis indicates the enriched factor in each pathway and term. The circle size indicates the number of genes. The color bar indicates the adjusted p-value, blue and red represent higher and lower values, respectively. GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes.

FIGURE 5
www.frontiersin.org

Figure 5 Schematic illustration of phenylpropanoid pathway with the expression levels of 11 DEGs (COMT, CCOMT, CHS, CHI, HCT, REF1, and DFR) and the contents of 7 metabolites that increased after NaCl treatment (catechin, chlorogenic acid, glycerin, isovitexin, p-coumaric acid, syringic acid, t-ferulic acid, and vitexin). The color scale indicates the log2-fold change (FC) value of each DEG, green and red represent lower and higher expression levels, respectively. The color blocks indicate FC of expression levels of paralogs of each DEG. The bar plots in the square show the comparison of metabolite contents between the control (gray) and NaCl treatment (black). The different uppercase letters indicate statistical significance (P < 0.05). The dotted line indicates an enzymic reaction that has not been fully identified in plants. COMT, catechol O-methyltransferase; CCOMT, caffeoyl-CoA O-methyltransferase; CHS, chalcone synthase; CHI, chalcone isomerase; HCT, shikimate O-hydroxycinnamoyl transferase; REF1, coniferyl-aldehyde dehydrogenase; DFR, Dihydroflavonol 4-reductase;.

The 11 DEGs encode seven enzymes: catechol O-methyltransferase (COMT: Vradi03g00001673.1, Vradi02g00004009.1), caffeoyl-CoA O-methyltransferase (CCOMT: Vradi09g00001754.1), chalcone synthase (CHS: Vradi09g00003301.1, Vradi09g00003302.1, Vradi09g00003303.1), chalcone isomerase (CHI: Vradi03g00001100.1), shikimate O-hydroxycinnamoyl transferase (HCT: Vradi09g00002363.1), coniferyl-aldehyde dehydrogenase (REF1: Vradi09g00002363.1, Vradi02g00003639.1), and dihydroflavonol 4-reductase (DFR: Vradi07g00001339.1) (Table S3 and Figure 5). In addition, the 11 DEGs were significantly upregulated under salinity conditions, and their Log2 FC values were COMT (1.75, 1.31), CCOMT (1.38), CHI (1.18), CHS (1.14, 1.13, 1.1), HCT (2.35), REF1 (1.85, 1), and DFR (1.46) (Figure 5). For validation of RNA-seq data, we conducted qRT-PCR and found that the result of qRT-PCR well agree with that of RNA-seq, indicating the reliability of the RNA-seq data (Figure 6).

FIGURE 6
www.frontiersin.org

Figure 6 Validation of the RNA expression of DEGs by qRT-PCR analysis.

Discussion

Mungbean sprouts are a good source of natural antioxidants, including catechin, chlorogenic acid, vitexin, isovitexin, and p-coumaric acid (Kim et al., 2009b; Guo et al., 2012; Gan et al., 2016; Singh et al., 2017). These secondary metabolites provide various health benefits for the human body, such as anticancer, antioxidant, antimicrobial, anti-inflammatory, and antitumor effects (Gorzynik-Debicka et al., 2018; Stagos, 2020; Teodor et al., 2020). In plants, abiotic stresses, such as salinity conditions, significantly influence secondary metabolite accumulation (Ksouri et al., 2008; Yang et al., 2018; Sharma et al., 2019; Chrysargyris et al., 2020). Many studies have demonstrated that mild NaCl treatment can improve secondary metabolite contents in various vegetables (Lim et al., 2012; Hassini et al., 2017; Sarker and Oba, 2018; Banik and Bhattacharjee, 2020; Benincasa et al., 2021). In the present study, we investigated the interactive effects of salinity stress on transcriptomic and metabolomic changes in mungbean sprouts.

NaCl treatment significantly suppressed mungbean sprout development (Figure 1). The root and hypocotyl length was reduced in both cultivars as NaCl concentration increased from 50 to 200 mM (Figures 2A, B). These results agree with previous findings that osmotic imbalance caused by salinity results in a significant reduction in root and hypocotyl length in maize (Khodarahmpour et al., 2012), wheat (Akbari et al., 2007), cowpea (Dantas et al., 2005), weed species (Hakim et al., 2011), tomato (Abdel-Farid et al., 2020), and mungbean (Saha et al., 2010). This inhibition of root growth by salinity stress is obviously a drawback in plants (Yin et al., 2022). However, in bean sprout production, long roots are considered a totally undesirable feature to reduce sprout texture, whereas a certain length and thickness of hypocotyl are required (Price, 1988). A previous study suggests that the most suitable features of mungbean sprouts are short roots and hypocotyl longer and thicker than 5 cm in length and 2 mm in diameter (Buescher and Chang, 1982). In our observation, sunhwa treated with Na50 showed an increase in hypocotyl thickness from 2.7 to 3.3 mm, maintaining a hypocotyl length over 8 cm (Figures 1B, C), indicating that mild NaCl treatment can promote the morphological quality of mungbean sprouts in consumers’ favor.

The total flavonoid contents increased in sunhwa at Na50 (Figure 2B), which corroborates previous findings that NaCl treatment can enhance the total flavonoid contents in wheat (Kiani et al., 2021), tomato (Abdel-Farid et al., 2020), and rapeseed sprout (Falcinelli et al., 2017). However, the total flavonoid contents gradually decreased in dahyeon as NaCl concentration increased (Figure 3B). This opposite response might be caused by sensitivity/insensitivity of the two cultivars to salt stress. Sunhwa (VC1973A) is a reference mungbean cultivar developed by the Asian Vegetable Research and Development Center ensuring high yield and insensitivity against diverse stresses (Habibzadeh and Yagoob, 2014; Kang et al., 2014; Yoseph Ganta et al., 2021). Secondary metabolites can be differentially regulated in salt-tolerant or sensitive cultivars under salinity condition (Borghesi et al., 2011; Linić et al., 2019). In the Brassicaceae family, the contents of phenolic acid were higher in more salt-tolerant varieties, such as white cabbage and kale, than in salt-sensitive varieties (Linić et al., 2019). In tomatoes, the accumulation of total anthocyanin contents differed depending on genotypes under NaCl treatment (Borghesi et al., 2011). Our results, along with those from other studies, reveal that the genotype with high salinity tolerance might be a valuable breeding material for high contents of antioxidant compounds in mungbean sprouts.

In this study, the antioxidant activity of mungbean sprouts was measured using two colorimetric methods (Figures 2C, D). DPPH and ABTS are both widely used for measuring the antioxidant activity of plant extraction (Floegel et al., 2011). While ABTS assay is suitable to measure both hydrophilic and lipophilic antioxidants, DPPH assay is more suited to measure hydrophobic systems (Floegel et al., 2011). We found that the antioxidant activities obtained by ABTS and DPPH assays presented slightly different trends (Figures 2C, D). In the ABTS assay, the antioxidant activity gradually declined as NaCl concentration increased, whereas no significant differences were observed in the DPPH assay under mild concentrations of NaCl (Figures 2C, D). This difference between the two assays indicates that the antioxidants mainly affected by mild NaCl treatment (0-100 mM) might be hydrophilic phytochemicals in mungbean sprouts.

Although a few studies have reported the effects of NaCl treatment on the total phenol and flavonoid contents in mungbean sprouts, the changes in the contents of individual phytochemicals along with transcriptomic changes remain unexplored (Koodkaew, 2019; Mankar et al., 2021). We conducted UHPLC analysis using 14 phenylpropanoid compounds in different groups (flavonoids and phenolic acid), which have been reported to be responsible for antioxidant activity in mungbean seeds and sprouts (Kim et al., 2009b; Guo et al., 2012; Gan et al., 2016; Singh et al., 2017). We found that most flavonoid and phenolic acid compounds significantly increased under salinity conditions, including catechin, chlorogenic acid, isovitexin, vitexin, p-coumaric acid, syringic acid, and t-ferulic acid (Table 1). In buckwheat sprouts, the contents of rutin and vitexin were enhanced under 50 mM NaCl treatment (Lim et al., 2012). The contents of catechin were increased in rice treated with 250 mM NaCl, and chlorogenic, p-coumaric, syringic, and t-ferulic acids were also increased by NaCl treatments in the Amaranthus, brassica family, rapeseed, and rice (Hassini et al., 2017; Sarker and Oba, 2018; Banik and Bhattacharjee, 2020; Benincasa et al., 2021). In contrast, salinity treatment can decrease some polyphenols in broccoli and romaine lettuce (Kim et al., 2008; López-Berenguer et al., 2009). In rice, while gallic acid, syringic acid, and catechin increased with 250 mM NaCl treatment, myricetin and quercetin decreased (Banik and Bhattacharjee, 2020). We also observed decreased contents of myricetin and quercetin under salinity conditions in the mungbean sprout (Table 1). Our findings suggest that the biosynthetic pathways of flavonoids and phenolic acids respond differently to salinity stress, which might be regulated by the key enzymes involved in each pathway.

According to the transcriptomic results, some genes (151 downregulated and 528 upregulated) were differentially regulated by NaCl treatment (Figure 3). We found that the expression of the genes related to the biosynthesis of the secondary metabolites was significantly increased by NaCl treatment (Table S2). DFR, REF1, CCOMT, COMT, HCT, CHS, and CHI genes were significantly highly expressed in mungbean sprouts under 50 mM NaCl treatment compared to the control, which corresponded to the changes in the secondary metabolite contents (Figure 5). In a previous study, the expression of TaDFR and TaCHI genes increased following 50 mM NaCl treatment in wheat sprouts (Cuong et al., 2020). In pak choi, DFR, CHS, and CHI genes were upregulated under 250 mM NaCl treatment (Yun et al., 2019). The expression of CHS was enhanced in S. nigrum as NaCl concentrations increased from 50 to 150 mM (Ben Abdallah et al., 2016). The CCOMT and COMT genes showed increased expression under salinity in various plants, including basil leaves, barley root, and tomato root (Sugimoto and Takeda, 2009; Manaa et al., 2011; Rastogi et al., 2019). In addition, many plant TF families were identified as upregulated DEGs, including ERFs, WRKYs, and MYBs (Figure 7). These TF families often stimulate the metabolic pathway under environmental stress, such as drought and salinity, resulting in more accumulated secondary metabolites in plants (Espley et al., 2007; Ilk et al., 2015; Qiu et al., 2016; Wei et al., 2017; Upadhyay et al., 2018; Li et al., 2019). Our TF DEGs were enriched to “response to stress,” “response to osmotic stress,” “response to salt stress,” “response to abiotic stimulus,” among other responses, and the expressions of these TF DEGs were increased up to log2FC 4.5 (log2FC 2.0 on average) (Table S4). In Arabidopsis, the overexpression of MYB111, MYB1D, and Mdmyb10 TFs significantly increased the expression of COMT, CHS, and CHI genes under salinity stress (Gao et al., 2011; Li et al., 2019). The silencing of OscWRKY1 resulted in reduced expression of PAL, COMT, and 4CL transcripts in Arabidopsis (Joshi et al., 2022). In potatoes, CHS, DFR, and ANS genes were upregulated by IbMYB1 overexpression under osmotic stress (Cheng et al., 2013). In this study, four paralogous genes encoding MYB14 were detected as up-regulated DEGs, that had been reported to be positively correlated with key flavonoid synthetic genes, including ANS, ANR, DFR, CHI, and CHS (Liu et al., 2018). WRKY72 that had been reported to be associated with the expression of ANS and DFR was also identified as upregulated TF DEG (Hu et al., 2020). Our transcriptomic results show that the biosynthetic pathways of phenylpropanoids, especially flavonoids, are regulated at the transcription level under salinity stress.

FIGURE 7
www.frontiersin.org

Figure 7 The number of transcription factors identified in differentially expressed genes. Red represents upregulated TF DEGs, and blue represents downregulated TF DEG.

Conclusion

We found that salinity stress can enhance the accumulation of phenylpropanoids and the expression of their responsive genes in mungbean sprouts. Our results indicate that mild NaCl treatment can be a practical approach to increasing the nutritional value of mungbean sprouts. Our comprehensive metabolomic and transcriptomic data provide insights into the complex regulatory mechanisms for salinity tolerance, as well as valuable information for improving the nutritional quality in legume crops.

Data availability statement

The data presented in the study are deposited in the NCBI SRA repository, accession number: https://www.ncbi.nlm.nih.gov/bioproject/PRJNA853140.

Author contributions

IL and MK designed the experiments and wrote the manuscript. IL and MK collected the material data and conducted experiments. BK performed the qualifying and quantifying of the metabolites. IL analyzed the transcriptomic results and contributed to the illustration. JH supervised and revised the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by a grant from the National Research Foundation of Korea (NRF), funded by the Korean government (MSIT; no. 2021R1C1C1004233).

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.

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.

Supplementary material

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

References

Abdel-Farid, I. B., Marghany, M. R., Rowezek, M. M., Sheded, M. G. (2020). Effect of salinity stress on growth and metabolomic profiling of cucumis sativus and solanum lycopersicum. Plants 9, 1626. doi: 10.3390/plants9111626

CrossRef Full Text | Google Scholar

Akbari, G., Sanavy, S. A. M. M., Yousefzadeh, S. (2007). Effect of auxin and salt stress (NaCl) on seed germination of wheat cultivars (Triticum aestivum l.). Pak J. Biol. Sci. 10, 2557–2561. doi: 10.3923/pjbs.2007.2557.2561

PubMed Abstract | CrossRef Full Text | Google Scholar

Banik, N., Bhattacharjee, S. (2020). Complementation of ROS scavenging secondary metabolites with enzymatic antioxidant defense system augments redox-regulation property under salinity stress in rice. Physiol. Mol. Biol. Plants 26, 1623–1633. doi: 10.1007/s12298-020-00844-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Ben Abdallah, S., Aung, B., Amyot, L., Lalin, I., Lachâal, M., Karray-Bouraoui, N., et al. (2016). Salt stress (NaCl) affects plant growth and branch pathways of carotenoid and flavonoid biosyntheses in solanum nigrum. Acta Physiol. Plant 38, 72. doi: 10.1007/s11738-016-2096-8

CrossRef Full Text | Google Scholar

Benincasa, P., Bravi, E., Marconi, O., Lutts, S., Tosti, G., Falcinelli, B. (2021). Transgenerational effects of salt stress imposed to rapeseed (Brassica napus var. oleifera del.) plants involve greater phenolic content and antioxidant activity in the edible sprouts obtained from offspring seeds. Plants 10, 932. doi: 10.3390/plants10050932

PubMed Abstract | CrossRef Full Text | Google Scholar

Bolger, A. M., Lohse, M., Usadel, B. (2014). Trimmomatic: a flexible trimmer for illumina sequence data. Bioinformatics 30, 2114–2120. doi: 10.1093/bioinformatics/btu170

PubMed Abstract | CrossRef Full Text | Google Scholar

Borghesi, E., González-Miret, M. L., Escudero-Gilete, M. L., Malorgio, F., Heredia, F. J., Meléndez-Martínez, A. J. (2011). Effects of salinity stress on carotenoids, anthocyanins, and color of diverse tomato genotypes. J. Agric. Food Chem. 59, 11676–11682. doi: 10.1021/jf2021623

PubMed Abstract | CrossRef Full Text | Google Scholar

Buescher, R., Chang, J. (1982). Production of mung bean sprouts. Arkansas farm Res. Arkansas Agric. Experiment Station 31, 13.

Google Scholar

Chandran, A. K. N., Kim, J.-W., Yoo, Y.-H., Park, H. L., Kim, Y.-J., Cho, M.-H., et al. (2019). Transcriptome analysis of rice-seedling roots under soil–salt stress using RNA-seq method. Plant Biotechnol. Rep. 13, 567–578. doi: 10.1007/s11816-019-00550-3

CrossRef Full Text | Google Scholar

Cheng, Y.-J., Kim, M.-D., Deng, X.-P., Kwak, S.-S., Chen, W. (2013). Enhanced salt stress tolerance in transgenic potato plants expressing IbMYB1, a sweet potato transcription factor. J. Microbiol. Biotechnol. 23, 1737–1746. doi: 10.4014/jmb.1307.07024

PubMed Abstract | CrossRef Full Text | Google Scholar

Chrysargyris, A., Mikallou, M., Petropoulos, S., Tzortzakis, N. (2020). Profiling of essential oils components and polyphenols for their antioxidant activity of medicinal and aromatic plants grown in different environmental conditions. Agronomy 10, 727. doi: 10.3390/agronomy10050727

CrossRef Full Text | Google Scholar

Cui, M. H., Yoo, K. S., Hyoung, S., Nguyen, H. T. K., Kim, Y. Y., Kim, H. J., et al. (2013). An arabidopsis R2R3-MYB transcription factor, AtMYB20, negatively regulates type 2C serine/threonine protein phosphatases to enhance salt tolerance. FEBS Lett. 587, 1773–1778. doi: 10.1016/j.febslet.2013.04.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Cuong, D. M., Kwon, S.-J., Nguyen, B. V., Chun, S. W., Kim, J. K., Park, S. U. (2020). Effect of salinity stress on phenylpropanoid genes expression and related gene expression in wheat sprout. Agronomy 10, 390. doi: 10.3390/agronomy10030390

CrossRef Full Text | Google Scholar

Dai, X., Xu, Y., Ma, Q., Xu, W., Wang, T., Xue, Y., et al. (2007). Overexpression of an R1R2R3 MYB gene, OsMYB3R-2, increases tolerance to freezing, drought, and salt stress in transgenic arabidopsis. Plant Physiol. 143, 1739–1751. doi: 10.1104/pp.106.094532

PubMed Abstract | CrossRef Full Text | Google Scholar

Dantas, B. F., Ribeiro, L., de, S., Aragão, C. A. (2005). Physiological response of cowpea seeds to salinity stress. Rev. Bras. sementes 27, 144–148. doi: 10.1590/S0101-31222005000100018

CrossRef Full Text | Google Scholar

Ebert, A. W., Chang, C.-H., Yan, M.-R., Yang, R.-Y. (2017). Nutritional composition of mungbean and soybean sprouts compared to their adult growth stage. Food Chem. 237, 15–22. doi: 10.1016/j.foodchem.2017.05.073

PubMed Abstract | CrossRef Full Text | Google Scholar

Espley, R. V., Hellens, R. P., Putterill, J., Stevenson, D. E., Kutty-Amma, S., Allan, A. C. (2007). Red colouration in apple fruit is due to the activity of the MYB transcription factor, MdMYB10. Plant J. 49, 414–427. doi: 10.1111/j.1365-313X.2006.02964.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Eulgem, T., Somssich, I. E. (2007). Networks of WRKY transcription factors in defense signaling. Curr. Opin. Plant Biol. 10, 366–371. doi: 10.1016/j.pbi.2007.04.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Falcinelli, B., Sileoni, V., Marconi, O., Perretti, G., Quinet, M., Lutts, S., et al. (2017). Germination under moderate salinity increases phenolic content and antioxidant activity in rapeseed (Brassica napus var oleifera del.) sprouts. Molecules 22, E1377. doi: 10.3390/molecules22081377

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, X.-D., Wang, J.-Q., Yang, N., Dong, Y.-Y., Liu, L., Wang, F.-W., et al. (2013). Gene expression profiling of soybean leaves and roots under salt, saline–alkali and drought stress by high-throughput illumina sequencing. Gene 512, 392–402. doi: 10.1016/j.gene.2012.09.100

PubMed Abstract | CrossRef Full Text | Google Scholar

Floegel, A., Kim, D.-O., Chung, S.-J., Koo, S. I., Chun, O. K. (2011). Comparison of ABTS/DPPH assays to measure antioxidant capacity in popular antioxidant-rich US foods. J. Food Composition Anal. 24, 1043–1048. doi: 10.1016/j.jfca.2011.01.008

CrossRef Full Text | Google Scholar

Gan, R.-Y., Lui, W.-Y., Chan, C.-L., Corke, H. (2017). Hot air drying induces browning and enhances phenolic content and antioxidant capacity in mung bean (Vigna radiata l.) sprouts. J. Food Process. Preservation 41, e12846. doi: 10.1111/jfpp.12846

CrossRef Full Text | Google Scholar

Gan, R.-Y., Wang, M.-F., Lui, W.-Y., Wu, K., Corke, H. (2016). Dynamic changes in phytochemical composition and antioxidant capacity in green and black mung bean (Vigna radiata) sprouts. Int. J. Food Sci. Technol. 51, 2090–2098. doi: 10.1111/ijfs.13185

CrossRef Full Text | Google Scholar

Gao, F., Yao, H., Zhao, H., Zhou, J., Luo, X., Huang, Y., et al. (2016). Tartary buckwheat FtMYB10 encodes an R2R3-MYB transcription factor that acts as a novel negative regulator of salt and drought response in transgenic arabidopsis. Plant Physiol. Biochem. 109, 387–396. doi: 10.1016/j.plaphy.2016.10.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, J.-J., Zhang, Z., Peng, R.-H., Xiong, A.-S., Xu, J., Zhu, B., et al. (2011). Forced expression of Mdmyb10, a myb transcription factor gene from apple, enhances tolerance to osmotic stress in transgenic arabidopsis. Mol. Biol. Rep. 38, 205–211. doi: 10.1007/s11033-010-0096-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Gayacharan., Archak, S., Gupta, K., Gupta, V., Tyagi, V., Singh, K. (2020). “Mungbean genetic resources and utilization,” in The mungbean genome, 9–25. (New York:Springer, Cham). doi: 10.1007/978-3-030-20008-4_2

CrossRef Full Text | Google Scholar

Gorzynik-Debicka, M., Przychodzen, P., Cappello, F., Kuban-Jankowska, A., Marino Gammazza, A., Knap, N., et al. (2018). Potential health benefits of olive oil and plant polyphenols. Int. J. Mol. Sci. 19, 686. doi: 10.3390/ijms19030686

CrossRef Full Text | Google Scholar

Guo, X., Li, T., Tang, K., Liu, R. H. (2012). Effect of germination on phytochemical profiles and antioxidant activity of mung bean sprouts (Vigna radiata). J. Agric. Food Chem. 60, 11050–11055. doi: 10.1021/jf304443u

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, L., Yang, R., Wang, Z., Guo, Q., Gu, Z. (2014). Effect of NaCl stress on health-promoting compounds and antioxidant activity in the sprouts of three broccoli cultivars. Int. J. Food Sci. Nutr. 65, 476–481. doi: 10.3109/09637486.2013.860583

PubMed Abstract | CrossRef Full Text | Google Scholar

Habibzadeh, Y., Yagoob, M. (2014). The effects of water deficit stress on protein yield of mung bean genotypes. Peak J. Agri. Sci. 2, 30–35.

Google Scholar

Hakim, M. A., Juraimi, A. S., Hanafi, M. M., Selamat, A., Ismail, M. R., Karim, S. M. R. (2011). Studies on seed germination and growth in weed species of rice field under salinity stress. J. Environ. Biol. 32, 529–536.

PubMed Abstract | Google Scholar

Ha, J., Satyawan, D., Jeong, H., Lee, E., Cho, K.-H., Kim, M. Y., et al. (2021). A near-complete genome sequence of mungbean (Vigna radiata l.) provides key insights into the modern breeding program. Plant Genome 14, e20121. doi: 10.1002/tpg2.20121

PubMed Abstract | CrossRef Full Text | Google Scholar

Hassini, I., Baenas, N., Moreno, D. A., Carvajal, M., Boughanmi, N., Martinez Ballesta, M. D. C. (2017). Effects of seed priming, salinity and methyl jasmonate treatment on bioactive composition of brassica oleracea var. capitata (white and red varieties) sprouts. J. Sci. Food Agric. 97, 2291–2299. doi: 10.1002/jsfa.8037

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, J., Fang, H., Wang, J., Yue, X., Su, M., Mao, Z., et al. (2020). Ultraviolet b-induced MdWRKY72 expression promotes anthocyanin synthesis in apple. Plant Sci. 292, 110377. doi: 10.1016/j.plantsci.2019.110377

PubMed Abstract | CrossRef Full Text | Google Scholar

Ilk, N., Ding, J., Ihnatowicz, A., Koornneef, M., Reymond, M. (2015). Natural variation for anthocyanin accumulation under high-light and low-temperature stress is attributable to the ENHANCER OF AG-4 2 (HUA2) locus in combination with PRODUCTION OF ANTHOCYANIN PIGMENT1 (PAP1) and PAP2. New Phytol. 206, 422–435. doi: 10.1111/nph.13177

PubMed Abstract | CrossRef Full Text | Google Scholar

Joshi, A., Jeena, G. S., Shikha, Kumar, R. S., Pandey, A., Shukla, R. K. (2022). Ocimum sanctum, OscWRKY1, regulates phenylpropanoid pathway genes and promotes resistance to pathogen infection in arabidopsis. Plant Mol. Biol 2 (3), 30–35. doi: 10.1007/s11103-022-01297-2

CrossRef Full Text | Google Scholar

Kang, Y. J., Kim, S. K., Kim, M. Y., Lestari, P., Kim, K. H., Ha, B.-K., et al. (2014). Genome sequence of mungbean and insights into evolution within vigna species. Nat. Commun. 5, 1–9. doi: 10.1038/ncomms6443

CrossRef Full Text | Google Scholar

Khodarahmpour, Z., Ifar, M., Motamedi, M. (2012). Effects of NaCl salinity on maize (Zea mays l.) at germination and early seedling stage. Afr. J. Biotechnol. 11, 298–304. doi: 10.4314/ajb.v11i2

CrossRef Full Text | Google Scholar

Kiani, R., Arzani, A., Mirmohammady Maibody, S. (2021). Polyphenols, flavonoids, and antioxidant activity involved in salt tolerance in wheat, aegilops cylindrica and their amphidiploids. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.646221

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, D., Choi, J., Lee, Y., Son, D., Moon, J., Oh, Y., et al. (2009a). A new mungbean cultivar, “Dahyeon” with many pod and high yielding. Korean J. Breed. Sci. 41, 36–39.

Google Scholar

Kim, H.-J., Fonseca, J. M., Choi, J.-H., Kubota, C., Kwon, D. Y. (2008). Salt in irrigation water affects the nutritional and visual properties of romaine lettuce (Lactuca sativa l.). J. Agric. Food Chem. 56, 3772–3776. doi: 10.1021/jf0733719

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, D.-K., Jeong, S. C., Gorinstein, S., Chon, S.-U. (2012). Total polyphenols, antioxidant and antiproliferative activities of different extracts in mungbean seeds and sprouts. Plant Foods Hum. Nutr. 67, 71–75. doi: 10.1007/s11130-011-0273-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, B. C., Lim, I., Jeon, S. Y., Kang, M., Ha, J. (2021). Effects of irrigation conditions on development of mungbean (Vigna radiata l.) sprouts. Plant Breed. Biotechnol. 9, 310–317. doi: 10.9787/PBB.2021.9.4.310

CrossRef Full Text | Google Scholar

Kim, D., Pertea, G., Trapnell, C., Pimentel, H., Kelley, R., Salzberg, S. L. (2013). TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 14, R36. doi: 10.1186/gb-2013-14-4-r36

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, D.-K., Son, D.-M., Chon, S.-U., Lee, K.-D., Kim, K.-H., Rim, Y.-S. (2009b). Phenolic compounds content and DPPH, ADH, ALDH activities of mungbean sprout based on growth temperature. Korean J. Of Crop Sci. 54, 1–6.

Google Scholar

Koodkaew, I. (2019). NaCl And glucose improve health-promoting properties in mung bean sprouts. Scientia Hortic. 247, 235–241. doi: 10.1016/j.scienta.2018.12.022

CrossRef Full Text | Google Scholar

Ksouri, R., Megdiche, W., Falleh, H., Trabelsi, N., Boulaaba, M., Smaoui, A., et al. (2008). Influence of biological, environmental and technical factors on phenolic content and antioxidant activities of Tunisian halophytes. Comptes Rendus Biologies 331, 865–873. doi: 10.1016/j.crvi.2008.07.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, S., Ayachit, G., Sahoo, L. (2020). Screening of mungbean for drought tolerance and transcriptome profiling between drought-tolerant and susceptible genotype in response to drought stress. Plant Physiol. Biochem. 157, 229–238. doi: 10.1016/j.plaphy.2020.10.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Liao, Y., Smyth, G. K., Shi, W. (2014). featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30, 923–930. doi: 10.1093/bioinformatics/btt656

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, B., Fan, R., Guo, S., Wang, P., Zhu, X., Fan, Y., et al. (2019). The arabidopsis MYB transcription factor, MYB111 modulates salt responses by regulating flavonoid biosynthesis. Environ. Exp. Bot. 166, 103807. doi: 10.1016/j.envexpbot.2019.103807

CrossRef Full Text | Google Scholar

Lim, I., Kim, B. C., Park, Y., Park, N. I., Ha, J. (2022). Metabolic and developmental changes in germination process of mung bean (Vigna radiata (L.) r. wilczek) sprouts under different water spraying interval and duration. J. Food Qual. 2022, e6256310. doi: 10.1155/2022/6256310

CrossRef Full Text | Google Scholar

Lim, J.-H., Park, K.-J., Kim, B.-K., Jeong, J.-W., Kim, H.-J. (2012). Effect of salinity stress on phenolic compounds and carotenoids in buckwheat (Fagopyrum esculentum m.) sprout. Food Chem. 135, 1065–1070. doi: 10.1016/j.foodchem.2012.05.068

PubMed Abstract | CrossRef Full Text | Google Scholar

Linić, I., Šamec, D., Grúz, J., Vujčić Bok, V., Strnad, M., Salopek-Sondi, B. (2019). Involvement of phenolic acids in short-term adaptation to salinity stress is species-specific among brassicaceae. Plants (Basel) 8, 155. doi: 10.3390/plants8060155

CrossRef Full Text | Google Scholar

Liu, L., Li, Y., She, G., Zhang, X., Jordan, B., Chen, Q., et al. (2018). Metabolite profiling and transcriptomic analyses reveal an essential role of UVR8-mediated signal transduction pathway in regulating flavonoid biosynthesis in tea plants (Camellia sinensis) in response to shading. BMC Plant Biol. 18, 233. doi: 10.1186/s12870-018-1440-0

PubMed Abstract | CrossRef Full Text | Google Scholar

López-Berenguer, C., Martínez-Ballesta, M., del, C., Moreno, D. A., Carvajal, M., García-Viguera, C. (2009). Growing hardier crops for better health: Salinity tolerance and the nutritional value of broccoli. J. Agric. Food Chem. 57, 572–578. doi: 10.1021/jf802994p

PubMed Abstract | CrossRef Full Text | Google Scholar

Maere, S., Heymans, K., Kuiper, M. (2005). BiNGO: a cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics 21, 3448–3449. doi: 10.1093/bioinformatics/bti551

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, W., Kim, J. K., Jia, C., Yin, F., Kim, H. J., Akram, W., et al. (2019). Comparative transcriptome and metabolic profiling analysis of buckwheat (Fagopyrum tataricum (L.) gaertn.) under salinity stress. Metabolites 9, 225. doi: 10.3390/metabo9100225

CrossRef Full Text | Google Scholar

Manaa, A., Ben Ahmed, H., Valot, B., Bouchet, J.-P., Aschi-Smiti, S., Causse, M., et al. (2011). Salt and genotype impact on plant physiology and root proteome variations in tomato. J. Exp. Bot. 62, 2797–2813. doi: 10.1093/jxb/erq460

PubMed Abstract | CrossRef Full Text | Google Scholar

Mankar, G. D., Wayase, U. R., Shelke, D. B., Nikam, T. D., Barmukh, R. B. (2021). Morphological, physiological, and biochemical responses to NaCl-induced salt stress in mungbean (Vigna radiata l.) varieties. Notulae Scientia Biologicae 13, 10936–10936. doi: 10.15835/nsb13210936

CrossRef Full Text | Google Scholar

Nair, R. M., Yang, R.-Y., Easdown, W. J., Thavarajah, D., Thavarajah, P., Hughes, J.D’A, et al. (2013). Biofortification of mungbean (Vigna radiata) as a whole food to enhance human health. J. Sci. Food Agric. 93, 1805–1813. doi: 10.1002/jsfa.6110

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakashima, K., Tran, L.-S. P., Van Nguyen, D., Fujita, M., Maruyama, K., Todaka, D., et al. (2007). Functional analysis of a NAC-type transcription factor OsNAC6 involved in abiotic and biotic stress-responsive gene expression in rice. Plant J. 51, 617–630. doi: 10.1111/j.1365-313X.2007.03168.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Nderitu, A. M., Dykes, L., Awika, J. M., Minnaar, A., Duodu, K. G. (2013). Phenolic composition and inhibitory effect against oxidative DNA damage of cooked cowpeas as affected by simulated in vitro gastrointestinal digestion. Food Chem. 141, 1763–1771. doi: 10.1016/j.foodchem.2013.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Price, T. V. (1988). Seed sprout production for human consumption — a review. Can. Inst Food Sci. Technol. J. 21, 57–65. doi: 10.1016/S0315-5463(88)70718-X

CrossRef Full Text | Google Scholar

Pudenz, M., Roth, K., Gerhauser, C. (2014). Impact of soy isoflavones on the epigenome in cancer prevention. Nutrients 6, 4218–4272. doi: 10.3390/nu6104218

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiu, Z., Wang, X., Gao, J., Guo, Y., Huang, Z., Du, Y. (2016). The tomato hoffman’s anthocyaninless gene encodes a bHLH transcription factor involved in anthocyanin biosynthesis that is developmentally regulated and induced by low temperatures. PloS One 11, e0151067. doi: 10.1371/journal.pone.0151067

PubMed Abstract | CrossRef Full Text | Google Scholar

Rastogi, S., Shah, S., Kumar, R., Vashisth, D., Akhtar, M. Q., Kumar, A., et al. (2019). Ocimum metabolomics in response to abiotic stresses: Cold, flood, drought and salinity. PloS One 14, e0210903. doi: 10.1371/journal.pone.0210903

PubMed Abstract | CrossRef Full Text | Google Scholar

R Core Team (2017). R: A language and environment for statistical computing (Vienna, Austria: R Foundation for Statistical Computing). Available at: https://www.R-project.org/.

Google Scholar

Robinson, M. D., McCarthy, D. J., Smyth, G. K. (2010). edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140. doi: 10.1093/bioinformatics/btp616

PubMed Abstract | CrossRef Full Text | Google Scholar

Saha, P., Chatterjee, P., Biswas, A. K. (2010). NaCl Pretreatment alleviates salt stress by enhancement of antioxidant defense system and osmolyte accumulation in mungbean (Vigna radiata l. wilczek). Indian J. Exp. Biol. 48, 593–600.

PubMed Abstract | Google Scholar

Sarker, U., Oba, S. (2018). Augmentation of leaf color parameters, pigments, vitamins, phenolic acids, flavonoids and antioxidant activity in selected amaranthus tricolor under salinity stress. Sci. Rep. 8, 12349. doi: 10.1038/s41598-018-30897-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Schneider, C. A., Rasband, W. S., Eliceiri, K. W. (2012). NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 9, 671–675. doi: 10.1038/nmeth.2089

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma, A., Shahzad, B., Rehman, A., Bhardwaj, R., Landi, M., Zheng, B. (2019). Response of phenylpropanoid pathway and the role of polyphenols in plants under abiotic stress. Molecules 24, 2452. doi: 10.3390/molecules24132452

CrossRef Full Text | Google Scholar

Shen, X., Guo, X., Guo, X., Zhao, D., Zhao, W., Chen, J., et al. (2017). PacMYBA, a sweet cherry R2R3-MYB transcription factor, is a positive regulator of salt stress tolerance and pathogen resistance. Plant Physiol. Biochem. 112, 302–311. doi: 10.1016/j.plaphy.2017.01.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Sherman, B. T., Hao, M., Qiu, J., Jiao, X., Baseler, M. W., Lane, H. C., et al. (2022). DAVID: a web server for functional enrichment analysis and functional annotation of gene lists, (2021 update). Nucleic Acids Res. 10, gkac194. doi: 10.1093/nar/gkac194

CrossRef Full Text | Google Scholar

Singh, B., Singh, J. P., Kaur, A., Singh, N. (2017). Phenolic composition and antioxidant potential of grain legume seeds: A review. Food Res. Int. 101, 1–16. doi: 10.1016/j.foodres.2017.09.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Singleton, V. L., Rossi, J. A. (1965). Colorimetry of total phenolics with phosphomolybdic-phosphotungstic acid reagents. Am. J. Enol Vitic. 16, 144–158.

Google Scholar

Stagos, D. (2020). Antioxidant activity of polyphenolic plant extracts. Antioxidants 9, 19. doi: 10.3390/antiox9010019

CrossRef Full Text | Google Scholar

Sugimoto, M., Takeda, K. (2009). Proteomic analysis of specific proteins in the root of salt-tolerant barley. Biosci Biotechnol Biochem. 73, 2762–2765. doi: 10.1271/bbb.90456

PubMed Abstract | CrossRef Full Text | Google Scholar

Teodor, E. D., Ungureanu, O., Gatea, F., Radu, G. L. (2020). The potential of flavonoids and tannins from medicinal plants as anticancer agents. Anti-Cancer Agents Med Chemistry- Anti-Cancer Agents) 20, 2216–2227. doi: 10.2174/1871520620666200516150829

CrossRef Full Text | Google Scholar

Tian, F., Yang, D.-C., Meng, Y.-Q., Jin, J., Gao, G. (2020). PlantRegMap: charting functional regulatory maps in plants. Nucleic Acids Res. 48, D1104–D1113. doi: 10.1093/nar/gkz1020

PubMed Abstract | CrossRef Full Text | Google Scholar

Upadhyay, A., Gaonkar, T., Upadhyay, A. K., Jogaiah, S., Shinde, M. P., Kadoo, N. Y., et al. (2018). Global transcriptome analysis of grapevine (Vitis vinifera l.) leaves under salt stress reveals differential response at early and late stages of stress in table grape cv. Thompson seedless. Plant Physiol. Biochem. 129, 168–179. doi: 10.1016/j.plaphy.2018.05.032

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, G., Zhu, Q., Meng, Q., Wu, C. (2012). Transcript profiling during salt stress of young cotton (Gossypium hirsutum) seedlings via solexa sequencing. Acta Physiol. Plant 34, 107–115. doi: 10.1007/s11738-011-0809-6

CrossRef Full Text | Google Scholar

Wei, Q., Zhang, F., Sun, F., Luo, Q., Wang, R., Hu, R., et al. (2017). A wheat MYB transcriptional repressor TaMyb1D regulates phenylpropanoid metabolism and enhances tolerance to drought and oxidative stresses in transgenic tobacco plants. Plant Sci. 265, 112–123. doi: 10.1016/j.plantsci.2017.09.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, Z., Gongbuzhaxi, Wang, C., Xue, F., Zhang, H., Ji, W. (2015). Wheat NAC transcription factor TaNAC29 is involved in response to salt stress. Plant Physiol. Biochem. 96, 356–363. doi: 10.1016/j.plaphy.2015.08.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Y., Guo, Y. (2018). Elucidating the molecular mechanisms mediating plant salt-stress responses. New Phytol. 217, 523–539. doi: 10.1111/nph.14920

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, L., Wen, K.-S., Ruan, X., Zhao, Y.-X., Wei, F., Wang, Q. (2018). Response of plant secondary metabolites to environmental factors. Molecules 23, 762. doi: 10.3390/molecules23040762

CrossRef Full Text | Google Scholar

Yen, G.-C., Chen, H.-Y. (1995). Antioxidant activity of various tea extracts in relation to their antimutagenicity. J. Agric. Food Chem. 43, 27–32. doi: 10.1021/jf00049a007

CrossRef Full Text | Google Scholar

Yin, Y., Xu, J., He, X., Yang, Z., Fang, W., Tao, J. (2022). Role of exogenous melatonin involved in phenolic acid metabolism of germinated hulless barley under NaCl stress. Plant Physiol. Biochem. 170, 14–22. doi: 10.1016/j.plaphy.2021.11.036

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoseph Ganta, T., Mekbib, D. F., Tadele, P. Z. (2021). Genetic diversity, drought tolerance, and genotype by environment interaction of mung bean [vigna radiata (l.) wilczek] genotypes in ethiopia (Ethiopia: University of Haramaya).

Google Scholar

Yuan, G., Wang, X., Guo, R., Wang, Q. (2010). Effect of salt stress on phenolic compounds, glucosinolates, myrosinase and antioxidant activity in radish sprouts. Food Chem. 121, 1014–1019. doi: 10.1016/j.foodchem.2010.01.040

CrossRef Full Text | Google Scholar

Yun, Y. B., Jung, H.-J., Rahim, M. A., Park, J.-I., Kuk, Y. I. (2019). Molecular analysis of genes related to phenylpropanoid and ascorbate biosynthesis in salt and UV-b treated pak choi grown under LEDs. Botany 97, 513–519. doi: 10.1139/cjb-2018-0183

CrossRef Full Text | Google Scholar

Zhishen, J., Mengcheng, T., Jianming, W. (1999). The determination of flavonoid contents in mulberry and their scavenging effects on superoxide radicals. Food Chem. 64, 555–559. doi: 10.1016/S0308-8146(98)00102-2

CrossRef Full Text | Google Scholar

Keywords: mungbean sprout, salinity stress, ultra-high-performance liquid chromatography, phenylpropanoid compound, RNA-seq, gene expression

Citation: Lim I, Kang M, Kim BC and Ha J (2022) Metabolomic and transcriptomic changes in mungbean (Vigna radiata (L.) R. Wilczek) sprouts under salinity stress. Front. Plant Sci. 13:1030677. doi: 10.3389/fpls.2022.1030677

Received: 29 August 2022; Accepted: 27 September 2022;
Published: 17 October 2022.

Edited by:

Moonhyuk Kwon, Gyeongsang National University, South Korea

Reviewed by:

Sang Un Park, Chungnam National University, South Korea
Yue Wang, Zhejiang University, China
Sang Min Kim, Korea Institute of Science and Technology (KIST), South Korea

Copyright © 2022 Lim, Kang, Kim and Ha. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Jungmin Ha, j.ha@gwnu.ac.kr

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.