Skip to main content
Download PDF

Integrated metabolome and transcriptome analyses reveal the role of BoGSTF12 in anthocyanin accumulation in Chinese kale (Brassica oleracea var. alboglabra)

BMC Plant Biology volume 24, Article number: 335 (2024) Cite this article

Abstract

Background

The vivid red, purple, and blue hues that are observed in a variety of plant fruits, flowers, and leaves are produced by anthocyanins, which are naturally occurring pigments produced by a series of biochemical processes occurring inside the plant cells. The purple-stalked Chinese kale, a popular vegetable that contains anthocyanins, has many health benefits but needs to be investigated further to identify the genes involved in the anthocyanin biosynthesis and translocation in this vegetable.

Results

In this study, the purple- and green-stalked Chinese kale were examined using integrative transcriptome and metabolome analyses. The content of anthocyanins such as cyanidin-3-O-(6″-O-feruloyl) sophoroside-5-O-glucoside, cyanidin-3,5-O-diglucoside (cyanin), and cyanidin-3-O-(6″-O-p-hydroxybenzoyl) sophoroside-5-O-glucoside were considerably higher in purple-stalked Chinese kale than in its green-stalked relative. RNA-seq analysis indicated that 23 important anthocyanin biosynthesis genes, including 3 PAL, 2 C4H, 3 4CL, 3 CHS, 1 CHI, 1 F3H, 2 FLS, 2 F3’H, 1 DFR, 3 ANS, and 2 UFGT, along with the transcription factor BoMYB114, were significantly differentially expressed between the purple- and green-stalked varieties. Results of analyzing the expression levels of 11 genes involved in anthocyanin production using qRT-PCR further supported our findings. Association analysis between genes and metabolites revealed a strong correlation between BoGSTF12 and anthocyanin. We overexpressed BoGSTF12 in Arabidopsis thaliana tt19, an anthocyanin transport mutant, and this rescued the anthocyanin-loss phenotype in the stem and rosette leaves, indicating BoGSTF12 encodes an anthocyanin transporter that affects the accumulation of anthocyanins.

Conclusion

This work represents a key step forward in our understanding of the molecular processes underlying anthocyanin production in Chinese kale. Our comprehensive metabolomic and transcriptome analyses provide important insights into the regulatory system that controls anthocyanin production and transport, while providing a foundation for further research to elucidate the physiological importance of the metabolites found in this nutritionally significant vegetable.

Peer Review reports

Introduction

One of the most economically important crop families, Brassicaceae, contains 372 genera and 4,060 species, including Brassica oleracea, B. napus, B. rapa, and the model plant Arabidopsis thaliana [1]. Brassicaceae crops are consumed by people all over the world and are rich in many nutritionally beneficial compounds, including phenolics, vitamins, anthocyanins, glucosinolates, and carotenoids. Due to valuable concentrations of phenolic components, ascorbic acid, tocopherols, and carotenoids, these vegetables are a natural source of antioxidants that help protect the body from damage caused by reactive oxygen species [2,3,4,5,6].

Purple-stalked kale (B. oleracea var. alboglabra) is a variety of B. oleracea. Its color results from the presence of anthocyanins in the stalk epidermal cells [7]. Anthocyanins are water-soluble pigments that occur in a variety of plants and serve a wide range of biological purposes including coloration and stress protection. They give plants their characteristic red, pink, purple, or blue colors, and they also play important roles in resistance to biotic and abiotic stresses [8, 9]. Additionally, anthocyanins have potent antioxidant capabilities that are advantageous for human health [10]. A significant amount of research has been conducted to understand the anthocyanin biosynthesis pathway and its regulation. Key genes in the pathway include phenylalanine ammonia lyase (PAL), chalcone acetylase (CHS), chalcone isomerase (CHI), flavanone 3-hydroxylase (F3H), flavonoid 3’-hydroxylase (F3’H), dihydroflavonol 4-reductase (DFR), anthocyanin synthase (ANS), and UDP-glucose: flavonoid 3-O-glucosyltransferase (UFGT) [11, 12]. In addition, WRKY, basic helix-loop-helix protein (bHLH), MYB, WD40, and NAC transcription factors (TFs) have been identified and studied in higher plants as important in anthocyanin biosynthesis [13, 14]. Several different TFs control the production of anthocyanins in plants. Typically, MYB transcription factors that are involved in anthocyanin production belong to the subgroup MYB-bHLH-WD40 (MBW complex) [15]. Three different TF types make up this complex: WD40 repeat, bHLH, and MYB proteins [15]. Expression of the genes involved in anthocyanin biosynthesis is activated by the combined action of these transcription factors.

Additionally, environmental variables are closely linked to anthocyanin biosynthesis, metabolism, and storage [16], with temperature being a key external signal. In general, during low temperatures, genes involved in anthocyanin synthesis in plants are stimulated, resulting in an increase in anthocyanin content, whereas high temperatures can speed up anthocyanins degradation and result in the fading of plant color. Additionally, it has been demonstrated that some plants accumulate anthocyanins as a defense against adverse environmental conditions, such as cold temperatures [17,18,19,20,21].

Within the Brassicaceae family, genes related to anthocyanin synthesis and regulation have been identified mainly through mapping and transcriptome sequencing. In B. rapa, BrMYB2 [22, 23], BrMYBL2.1 [24], BrMYB114 [25], BrbHLH49 [26], and BrEGL3.2 [27] are associated with the purple trait. In B. juncea, BjTT8 controls the color of purple tumorous stem mustard [28], and BjPl1 is related to its purple leaf color [29]. In (A) thaliana, the transcription factors MYB114, MYB113, and MYB118 are involved in controlling the production of anthocyanins [30]. A number of genes that control some of the leaf colors in (B) oleracea have been identified: it was shown that DFR and Re genes control the red-leaf characteristic [31,32,33,34], and BoMYB2 regulates the production of the purple color in leaves of ornamental kale [35]. Although anthocyanins biosynthesis in B. oleracea is clear, anthocyanin accumulation is less understood.

In this study, using transcriptome and metabolome association analyses, we identified BoGSTF12 gene. Biochemical, genetic and molecular studies showed that BoGSTF12 was a transporter of anthocyanin. These findings will improve our understanding of the regulation of anthocyanin accumulation in Chinese kale and also provide useful gene resources for breeding.

Materials and methods

Plant materials and phenotypic analysis

Inbred lines of purple-stalked (HJJL, R) and green-stalked (ZSJL, G) Chinese kale (B. oleracea var. alboglabra) were used as the experimental materials in this study. The original seeds were provided by the Vegetable Research Institute, Guangdong Academy of Agricultural Sciences. Both varieties of Chinese kale were grown at the Baiyun experimental fields of the Vegetable Research Institute, Guangdong Academy of Agricultural Sciences, Guangzhou, China. The skin of the stalks was peeled off with a blade, then quickly stored in liquid nitrogen. Samples were collected at the same time for metabolome analysis, RNA sequencing (RNA-seq), and qRT-PCR validation. Three biological replicates were taken from each group of samples, each biological replicate consisting of a mixture of five plants.

Metabolite identification

A vacuum freeze-dryer (Scientz-100 F) was used to dry skin. The freeze-dried samples were ground using zirconium beads for 1.5 min at 30 Hz in a mixer mill (MM 400, Retsch). Then, 100 mg of lyophilized powder was dissolved in 1.2 mL of 70% methanol in water, vortexed six times for 30 s each (once every 30 min), and then stored at 4 °C overnight. Before ultra-high performance liquid chromatography-mass spectrometry (UPLC-MS)/MS analysis, the extracts were filtered (SCAA-104, 0.22 μm pore size) after centrifugation at 12,000 rpm for 10 min.

An UPLC-ESI-MS/MS system (UPLC, SHIMADZU Nexera X2; MS, Applied Biosystems 4500 Q TRAP) was used to analyze the sample extracts. The following analytical conditions were used: an Agilent SB-C18 column, 1.8 μm, 2.1 mm×100 mm; Solvent A, sterile clean water with 0.1% formic acid, and Solvent B, acetonitrile with 0.1% formic acid, made up the mobile phase. The starting conditions were 95% A and 5% B, followed by a linear gradient to 5% A, 95% B in less than 9 min, with the composition of 5% A, 95% B maintained for 1 min. Then, within 10 min, a composition of 95% A and 5% B was used and maintained for 2 min. The chosen flow velocity was 0.35 mL/min, the injection volume was 4 µL, and the column oven was adjusted to 40 °C. An alternate connection was made between the effluent and the electrospray ionization (ESI)-triple quadrupole-linear ion trap (QTRAP)-MS.

QTRAP UPLC/MS/MS System was equipped with an ESI Turbo Ion-Spray interface, operating in positive and negative ion modes, and was managed by Analyst 1.6.3 software (AB Sciex); linear ion trap (LIT) and triple quadrupole (QQQ) scans were acquired. The following ESI source operation parameters were used: Turbo spray ion source; 550 °C source temperature; 5500 V (positive ion mode)/−4500 V (negative ion mode) ion spray voltage; 50, 60, and 25 psi, respectively, for the ion source gas I (GSI), gas II (GSII), and curtain gas (CUR); and high collision-activated dissociation (CAD). The instrument calibration and mass calibration were carried out using solutions of 10 and 100 mol/L polypropylene glycol in the QQQ and LIT modes, respectively. For multiple reaction monitoring (MRM) investigations, QQQ scans were recorded with the collision gas (nitrogen) set to medium. Further declustering potential (DP) and collision energy (CE) optimization was used to perform DP and CE for individual MRM transitions. According to the metabolites eluted during each interval, a particular set of MRM transitions were observed.

RNA-seq analysis

Total RNA from stalk skin was isolated using an RNA Extraction Kit (Tiangen, Beijing, China). The Agilent Bioanalyzer 2100 system (Agilent Technologies, Palo Alto, CA, USA) was used to confirm the amount of RNA present. Using the purple- and green-stalked Chinese kale, six cDNA libraries (R1, R2, R3, G1, G2, and G3) were constructed. All samples were sequenced on the Illumina NovaSeq 6000 platform, PE150 model. The manufacturer’s instructions were followed in the construction and sequencing of the RNA-seq library, as previously described [36]. The adapter and low-quality sequences were removed from the raw readings. Clean reads were successfully mapped to the Brassica reference genome. Gene expression levels were calculated using the fragments per kilobases per million fragments (FPKM) technique. The following thresholds were used to identify DEGs: |log2(foldchange)| ≥ 1 and FDR < 0.01. Gene Ontology (GO) and KEGG enrichment analyses of DEGs were performed using the cluster Profiler R package (http://www.geneontology.org/) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) http://www.genome.ad.jp/kegg/.

qRT-PCR analysis

Total RNA was isolated from stalk skin using a Megan RNA Extraction Kit (Guangzhou Magen Biotechnology Co., Ltd.). Table S1 contains a list of the PCR primers used in this work. The control was the actin gene. Each qRT-PCR reaction contained 0.3 µL of the appropriate primers, 1 µL of template cDNA, 3.4 µL of ddH2O, and 5 µL of 2× ChamQ Universal SYBR qPCR Master Mix (Vazyme, Nanjing, China) in a total volume of 10 µL.

Integrated metabolome and transcriptome analyses

Correlation coefficients were computed between the metabolome and transcriptome datasets. These coefficients were derived from the log2-fold changes of individual metabolites and transcripts, using the EXCEL program. Cytoscape version 2.8.2 was used to visualize the relationships between the metabolome and transcriptome.

Statistical evaluation

The statistical analysis of variance was performed using GraphPad Prism8.0.lnk. The mean and standard deviation for three biological replicates of the data were displayed. The threshold point for significant differences was determined at p < 0.05.

Overexpression of BoGSTF12 in the Arabidopsis tt19 mutant

Using gene-specific primers (Table S1), the target gene was cloned and then transferred into the pCAMBIA1301 binary vector. Using the freeze–thaw technique, the vectors generated were introduced into Agrobacterium tumefaciens strain GV3101 before being transformed into the Arabidopsis tt19 mutant. The transgenic seeds were cultivated for 7 days on half-strength Murashige and Skoog plates with hygromycin before being planted. T3 homozygous lines were produced for future investigation. At various developmental stages, the phenotypic traits of the transgenic, mutant, and Arabidopsis wild-type (WT) plants were observed. To determine the anthocyanin contents, we used a method described in a previous study [37]. Samples were frozen in liquid nitrogen and were subsequently ground into a fine powder using a mortar and pestle under liquid nitrogen. After mixing the powder with 1% HCl in methanol, samples were incubated for 24 h at 4 °C. Anthocyanins were quantified by measuring absorbance at 530 and 657 nm.

Phylogenetic analysis

The protein sequences of GSTs from several species were used to construct a phylogenetic tree by MEGA 5 using the neighbor-joining method with a bootstrap value of 1000. Sequences used for phylogenetic tree analysis are listed in Table S2.

Results

Phenotypic analysis and anthocyanin metabolite identification in purple- and green-stalked Chinese kale

The most striking difference between the two varieties was the coloration of their stalks. The purple-stalked Chinese kale displayed a rich and vibrant purple color, whereas the green-stalked variety exhibited a lush and vivid green color (Fig. 1). Detailed genetic and biochemical analysis of these phenotypic differences can offer valuable insights into the physiological bases of these traits and potential implications for breeding and culinary applications of Chinese kale varieties.

Fig. 1
figure 1

Comparison of two Chinese kale phenotypes. (a) Green-stalked Chinese kale and (b) purple-stalked Chinese kale. Photographs show seven -week- old plants grown in the field as described in Methods

Full size image

The anthocyanin levels of the purple- vs. green-stalked Chinese kale differed noticeably. We performed a thorough metabolomic analysis to discover differences in the contents of anthocyanin-related metabolites between purple and green-stalked Chinese kale. A total of 88 differentially occurring metabolites were found, with 81 at higher levels and 7 at lower levels (Fig. 2a) in purple-stalked kale, which had more overall anthocyanin production than the green-stalked variety (Fig. 2b). The levels of anthocyanin-related compounds such as cyanidin-3-O-(6″-O-feruloyl) sophoroside-5-O-glucoside, cyanidin-3,5-O-diglucoside (cyanin), and cyanidin-3-O-(6″-O-p-hydroxybenzoyl) sophoroside-5-O-glucoside were significantly and prominently higher in the purple-stalked plants than in the green-stalked ones (Fig. 2c and Table S3).

Fig. 2
figure 2

Differential levels of anthocyanin metabolites in purple- vs. green-stalked Chinese kale. (a) Volcano plot demonstrating the statistical significance of the differences in metabolite levels in the two Chinese kale varieties. (b) Heatmap comparing metabolite contents in purple and green-stalked Chinese kale. Colors represent differential expression levels after normalization. The left side of the heatmap shows the metabolite classes as identified in the color key to the right. G1–G3 and R1–R3 represent the different samples. (c) Anthocyanins of differential metabolites in two types of Chinese kale

Full size image

Transcriptome analysis

Skin tissues of the purple- and green-stalked Chinese kale stalks were used to construct cDNA libraries, and RNA-seq analysis was carried out to determine the molecular process behind the production of anthocyanins in the samples. More specifically, the clean data for each sample amounted to 5.78 Gb, for a total of 37.22 Gb, and the Q30 baseline percentage was 94.74% and above (Table S4). Clean reads for each sample were sequenced along with the designated reference genome, and the mapped reads were between 91.19% and 92.02% (Table S5). Variable splicing prediction analysis, gene structure optimization analysis, and new gene discovery were done based on the comparison results, and 3,878 new genes were found. Of these, 2,404 genes were found to have functional annotations. In our differential gene expression analysis, we found a total of 5,826 DEGs, of which 2,408 were upregulated and 3,418 were downregulated in the purple-stalked variety relative to the green-stalked variety (Fig. 3, Table S6).

Fig. 3
figure 3

Volcano plot showing DEGs between purple- and green-stalked Chinese kale varieties

Full size image

Functional annotation of DEGs and enrichment analysis

We performed GO, COG, KEGG Orthology (KOG), and KEGG analyses to identify genes involved in the production of anthocyanins in Chinese kale. To predict the molecular function of these genes, we used biological process (BP), cellular component (CC), and molecular function (MF) classifications as the basis for the GO enrichment analysis (Fig. 4, Table S7). We found genes representing a total of 55 enriched terms, and the GO-CC annotation analysis detected 18 enrichments, including extracellular region (GO:0005576), membrane (GO:0016020), nucleoid (GO:0009295), and organelle part (GO:0043226). In GO-MF, we detected 16 enriched terms, including catalytic activity (GO:0003824), binding (GO:0003723), and molecular function regulator (GO:0098772). GO-BP enrichment analysis detected 21 enriched terms, including cellular process (GO:0009987), metabolic process (GO:0008152), response to stimulus (GO:0050896), and single-organism process (GO:0044699).

Fig. 4
figure 4

GO classification of DEGs between the two Chinese kale lines

Full size image

To further explore their biological roles, we mapped the DEGs to particular KEGG pathways. Specifically, 5,826 DEGs were mapped to 130 KEGG pathways. Significantly enriched pathways were plant-pathogen interaction (Ko04626), amino sugar and nucleotide sugar metabolism (Ko00520), plant hormone signal transduction (Ko04075), MAPK signaling pathway–plant (Ko04016), flavone and flavanol biosynthesis (Ko00944), flavonoid biosynthesis (Ko00941), glycerolipid metabolism (Ko00561), and brassinosteroid biosynthesis (Ko00905) (Fig. 5a, b, and Table S8). Notably, pigment-related pathways were significantly enriched among the DEGs, providing important clues about the processes underlying anthocyanin biosynthesis in Chinese kale.

Fig. 5
figure 5

KEGG enrichment pathways of DEGs between purple- and green-stalked Chinese kale. KEGG enrichment pathways of upregulated DEGs (a) and downregulated DEGs (b)

Full size image

Identification of candidate DEGs involved in anthocyanin biosynthesis

Anthocyanin biosynthesis is regulated by a variety of regulatory pathways in plants, with anthocyanin metabolism comprising one of the branches of the flavonoid metabolic pathway. In this study, we found 83 DEGs in this pathway, of which 23 were structural genes related to anthocyanin biosynthesis, including 3 PAL, 3 CHS, 3 4CL, 1 CHI, 2 F3’H, 1 F3H, 2 C4H, 2 FLS, 3 ANS, 1 DFR, and 2 UFGT (Fig. 6), which is consistent with the high-anthocyanin phenotype of purple-stalked kale. Therefore, it is likely that these particular DEGs contribute to the development of variously colored stalks in Chinese kale.

Fig. 6
figure 6

Analysis of the expression of genes involved in anthocyanin biosynthesis in Chinese kale. The illustration shows gene expression levels in the six cDNA libraries created from purple and green-stalked Chinese kale (red = higher expression) in the context of the anthocyanin biosynthesis process using structural genes showing variable expression

Full size image

Identification of relevant transcription factors

Numerous plant TFs control the production of anthocyanins. MYB TFs, including MYB11, MYB12, MYB113, and others, are part of MYB-bHLH-WD40 (MBW) complex [38]. Proteins in a MBW complex work together in a coordinated manner to activate the appropriate target genes involved in anthocyanin biosynthesis [30, 39]. The MYB TFs within the MBW complex regulate the expression of genes encoding key enzymes of the anthocyanin biosynthesis pathway, such as CHI, CHS, and DFR [40,41,42]. Activation of these genes leads to the synthesis of anthocyanin pigments and the subsequent development of colorful and nutritious plant tissues.

In this study, we identified a total of 92 anthocyanin-related TF genes, of which 56 were downregulated, and 36 were upregulated (Table S9). Most of these TFs were MYB, bHLH, and WD40 proteins, components of the MBW complex (Table S9). Of the differentially expressed TF genes related to the anthocyanin biosynthesis pathway, BoMYB114 (Bo6g100940) was shown by qRT-PCR verification to be more highly expressed in purple-stalked Chinese kale vs. the green-stalked variety (Fig. 7). Moreover, supporting our prediction of the importance of BoMYB114 in anthocyanin production in purple-stalked Chinese kale, we observed a higher expression level of three CHS genes, one CHI gene, and one DRF gene in purple-stalked Chinese kale than in green-stalked Chinese kale.

Fig. 7
figure 7

qRT-PCR verification of anthocyanin-related gene expression. Data, from left to right, are represented as relative expression and fragments per kilobase million (FPKM), respectively. Bars show means ± SD of biological replicates data

Full size image

These findings provide the basis to further analyze the role of BoMYB114 in controlling anthocyanin biosynthesis and accumulation in B. oleracea. Furthermore, we identified BoTT8 (Bo9g086910), a bHLH TF homologous to AtTT8, which was reported to be involved in anthocyanin synthesis in Arabidopsis [43]. Transcriptome results showed that Bo9g086910 gene expression in purple-stalked Chinese kale was higher than in green-stalked Chinese kale (Table S9), suggesting its involvement in anthocyanin biosynthesis. On the other hand, the WD40 TF gene Bo7g096780 was not found to be differentially expressed when comparing purple- and green-stalked Chinese kale.

Expression profiles of genes in the anthocyanin biosynthetic pathway

In order to verify the accuracy of the transcriptome data, qRT-PCR of 11 anthocyanin synthesis–related genes were performed, and these verification results were consistent with the transcriptome results (Fig. 7). The expression levels of C4H, 4CL, CHS, CHI, F3H, F3’H, DFR, MYB114, and GSTF12 in purple-stalked Chinese kale were upregulated in RNA-seq as well as in our qRT-PCR analysis, indicating that the transcriptome data were accurate.

Correlation analysis between selected DEGs and anthocyanins

To identify key candidate genes involved in anthocyanin accumulation in Chinese kale, we conducted correlation analyses between identified 25 DEGs (23 structural genes and 2 GSTF genes encoding glutathione S-transferases) and 3 metabolites, including cyanidin-3,5-O-diglucoside (cyanin), cyanidin-3-O-(6″-O-p-hydroxybenzoyl) sophoroside-5-O-glucoside, and cyanidin-3-O-(6″-O-feruloyl) sophoroside-5-O-glucoside. In all cases, we detected significant correlations (correlation coefficient, R2 > 0.8) between our selected DEGs and metabolites (Table S10).

Based on correlation coefficient results, our network analysis showed that all of the selected DEGs were strongly related to these three anthocyanin metabolites (Fig. 8a), demonstrating that these are core genes for anthocyanin accumulation in Chinese kale. Of these, two GSTF subfamily genes GSTF12 (Bo9g161480 and Bo2g013490) were highly correlated with anthocyanin metabolites (Fig. 8b and Table S10). Specifically, of these two genes, GSTF12 (Bo9g161480) showed a strong interaction with all three anthocyanins as compared to GSTF12 (Bo2g013490) (Fig. 8b).

Fig. 8
figure 8

Connection network between core genes and anthocyanin metabolites. (a) Network showing relationships between 25 core genes and 3 anthocyanins. (b) Network for GSTF12 (Bo9g161480 and Bo2g013490) genes and three anthocyanins. Solid lines, stronger interaction; dotted lines, weaker interaction

Full size image

Identification of BoGSTF12

There are 14 different subgroups of GST genes in plants, and the GSTF subfamily is a significant class of genes regulating anthocyanin transport. The results of association analysis indicate that GSTF12 may be involved in the transport of anthocyanins. Therefore, we constructed a phylogenetic tree and sequence alignment between GSTF12 and GSTFs that have been reported to be involved in anthocyanin transport in other species. The result of phylogenetic analysis (Fig. 9a) and sequence alignment (Fig. 9b) showed that BoGSTF12 was homologous with AtGSTF12, which was reported to be involved in anthocyanin transport. Combining the results of transcriptome analysis and qRT-PCR results, we found that the expression level of BoGSTF12 (Bo9g161480) in purple-stalked Chinese kale was higher than that in green-stalked Chinese kale. Because its expression level was consistent with the trend of anthocyanin content in each tissue of Chinese kale (Fig. 9c-d), we speculate that BoGSTF12 (Bo9g161480) regulates anthocyanin transport.

Fig. 9
figure 9

Phylogenetic analysis and expression pattern of BoGSTF12 (a) Phylogenetic analysis, (b) multiple sequence alignment of GST genes in different species (c) phenotypic representation of anthocyanin content in Chinese kale (d), Relative expression of BoGSTF12 in different plant parts of purple-stalked Chinese kale. Scale bar in (c) = 5 cm

Full size image

Overexpression of BoGSTF12 (Bo9g161480) gene in Arabidopsis

Because BoGSTF12 is homologous to Arabidopsis ATGSTF12, to further explore the function of BoGSTF12, it was introduced into the Arabidopsis tt19 mutant (SALK_105779, a deletion mutant of Arabidopsis AtGSTF12). This mutant lacks anthocyanins in the area between the stem and rosette leaves and lacks proanthocyanidins in the seed coat. Two transgenic lines (#6 and #9) were selected for subsequent experiments. BoGSTF12 rescued the anthocyanin-loss phenotype in the area between the stem and rosette leaves of Arabidopsis tt19, but it could not rescue this phenotype in the mutant seeds (Fig. 10a–d). To verify this finding, we measured anthocyanin contents in the Arabidopsis WT, ttl9, #6, and #9 plants. The stalks of the #6, #9, and WT plants contained more anthocyanins than did ttl9 plants (Fig. 10e). The above results indicated that BoGSTF12 was involved in the transport of anthocyanins, but not in the transport of proanthocyanidins, demonstrating a function different from that of AtGSTF12.

Fig. 10
figure 10

Overexpression of BoGSTF12 in the Arabidopsis tt19 mutant and measurements of total anthocyanins contents. (a) Phenotypes of wild-type (WT) Arabidopsis, (b) an Arabidopsis knockout mutant of the anthocyanin transporter AtGSTF12 (tt19), and (c, d) two transgenic lines of 35 S::BoGSTF12-FLAG in the tt19 background. (e) Total contents of anthocyanins as measured in the infiltration patches; data are means ± SD obtained from three biological replicates. The different letters denote significant differences according to one-way analysis of variance (ANOVA) (P < 0.05)

Full size image

Discussion

The goal of the current work was to clarify the molecular processes responsible for the eye-catching purple coloration seen in the stalks of purple-stalked Chinese kale, a popular leafy green vegetable with documented health advantages. Brassica vegetables such as broccoli (Brassica oleracea var. italica), heading Chinese cabbage (Brassica rapa ssp. pekinensis), mizuna (Brassica rapa var. japonica), and ornamental cabbage (Brassica oleracea var. acephala) have attracted a lot of attention because of their high anthocyanin levels [23, 44,45,46]. Using a thorough integration of metabolomic and transcriptome analyses, we have discovered important metabolites and genes involved in the manufacture of the pigments and other secondary metabolites that are responsible for the distinct purple-stalk phenotype.

Metabolomic analysis identified anthocyanins causing the purple-stalk phenotype in Chinese kale

Anthocyanins are some of the most important pigments affecting the color of plant tissues. The majority of research on anthocyanins in Brassica crops has focused on separating and distinguishing among metabolites [47]. The purple variety of cauliflower (Brassica oleracea var. botrytis) contains a cyanidin 3-(coumaryl-caffeyl) glucoside-5-(malonyl) pigment [48]. Additionally, using HPLC-ESI-MS/MS, red cabbage (Brassica oleracea var. capitata) was found to contain over 30 different cyanidin compounds [49].

In this study, our metabolomic analysis revealed a significant accumulation of anthocyanins, flavonoids, and related compounds in purple-stalked Chinese kale compared to the green-stalked variety. These findings are consistent with previous reports that anthocyanins are responsible for the purple coloration in various plant tissues, including leaves, flowers, and stems. Our identification of specific anthocyanin derivatives, such as cyanidin-3-O-(6″-O-feruloyl) sophoroside-5-O-glucoside, cyanidin-3,5-O-diglucoside (cyanin), and cyanidin-3-O-(6″-O-p-hydroxybenzoyl) sophoroside-5-O-glucoside (Fig. 2a-c), has provided insights into the origin of the diversity of pigments contributing to the purple color. Additionally, we identified other secondary metabolites, such as phenolic acids and flavonols, identified in higher quantities in the purple stalks, suggesting a potential interplay between different classes of metabolites in generating the observed pigmentation.

Transcriptome analysis revealed candidate genes involved in anthocyanin biosynthesis and transport

One of products of the flavonoid production pathway are anthocyanins, with many structural genes encoding TFs and enzymes that regulate anthocyanin metabolite production. Studies have shown that the majority of the fundamental genes in the anthocyanin production pathway are more highly expressed during vegetative growth in red cabbage than in green cabbage, resulting in a range of leaf colors [50]. Differently colored mizuna, Arabidopsis, and other plants have been shown to have coordinated expression of DFR, F3H, ANS, UFGT, and F3’H [51,52,53]. Dihydrokaempferol, dihydromyricetin, and dihydroquercetin have a particular substrate bias in DFR derived from different plants [53]. Additionally, ANS, an important enzyme near the end of the anthocyanin synthesis pathway, catalyzes the conversion of monochrome to colored anthocyanins [51]. In our research, comparable outcomes were seen, we found 23 structural genes specifically related to anthocyanin biosynthesis, including three PAL, two C4H, three 4CL, three CHS, one CHI, one DFR, three ANS, two FLS, two F3’H, one F3H, and two UFGT genes (Fig. 6; Table 1). The expression levels of C4H, CHI, CHS, F3’H, 4CL, F3H, and DFR were significantly higher in the purple-stalked Chinese kale than in the green variety. Our transcriptomic analysis complemented and verified the metabolomic data by providing a broader understanding of the genetic regulation underlying the accumulation of anthocyanin in the purple stalks. The upregulation of these key genes suggests an enhanced flux through the anthocyanin pathway, leading to the increased production of anthocyanin precursors and ultimately the pigments responsible for the purple color.

Table 1 DEGs involved in the regulation of anthocyanin levels in purple- and green-stalked Chinese kale
Full size table

Anthocyanins are very prone to degradation and instability. Glucosyltransferase (GST), which determines the position of glycosylation, is crucial for the stability and solubility of plant flower color and anthocyanins. As a result, anthocyanins can function as pigments in vacuoles in a similar manner to flavonoid 3-O-glucosyltransferase and anthocyanidin 3-O-glucosyltransferase [54, 55]. Delphinidin conversion to delphinidin-3-O-glucoside in Clitoria ternatea was also reported to be catalyzed by anthocyanidin 3-O-glucosyltransferase (UGT78K6) [56]. Additionally, Freesia hybrida Fh3GT1, which encodes UF3GT, is essential for the production of anthocyanin glycosides [53]. Similar results were found in our study as well (Fig. 6).

TFs are crucial for regulating the enzymes in the anthocyanin production pathway that are involved in the formation of vegetable, fruit, and flower color [49, 57, 58]. Previous research suggested that, in the flavonoid biosynthesis pathway, the late biosynthetic genes (LBGs) were controlled mostly by the transcriptional complexes made of WD-repeat/MYB/bHLH proteins [59, 60]. Increased expression of early biosynthetic genes (EBGs) and LBGs in Arabidopsis leaves overexpressing PAP1 suggest that the phenylpropanoid pathway is responsible for the increased flavonoid pigment content observed [61]. This shows that LBGs may not be the only genes in the flavonoid pathway that are regulated by transcriptional complexes in Arabidopsis. Similarly, the coordinated expression of TRANSPARENT TESTA8 (TT8) and MYB2-enhanced anthocyanin synthesis may occur by positively activating EBGs and LBGs such as F3H and CHS in red cabbage and F3’H ANS, DFR, LDOX, UFGT, and GST in purple-headed Chinese cabbage [50, 62]. Additionally, the expression of the structural genes F3’H, DFR, LDOX, UGT75C1, and GST12 in Arabidopsis was impacted by the interaction of MYB113 or MYB114 with TTG1 and bHLHs (GL3, EGL3, and TT8) [59]. In our study, we identified the TF BoMYB114, of the MYB family, which is known regulator of anthocyanin biosynthesis. The upregulation of these transcription factors suggests their involvement in orchestrating the transcriptional response leading to purple color development in purple-stalked Chinese kale. Additionally, we noted the increased expression level of three CHS genes, one CHI gene, and one DRF gene in our purple-stalked Chinese kale as compared to green-stalked kale, supporting our results that BoMYB114 is involved in anthocyanin synthesis in purple-stalked Chinese kale. These results provide the basis for further exploration of the functional analysis of BoMYB114 and its role in regulating anthocyanin production and accumulation in B. oleracea.

GST is a pivotal player in anthocyanin transport

Anthocyanins are produced in the cytosol and accumulate in the vacuole. The processes underpinning the intracellular transport of anthocyanins have been partially elucidated in recent years. GST involvement, membrane transport, and/or vesicle trafficking are required for the transfer of anthocyanins from cytosolic production to vacuolar accumulation [63]. The role of GSTs in anthocyanin transport and accumulation has been confirmed in maize [64], petunia [65], Arabidopsis [66], cyclamen [67], perilla [68], grape [69], apple [70], and litchi [71]. In our study, BoGSTF12 expression was highly correlated with anthocyanin content, and further functional analysis revealed that BoGSTF12 rescued the anthocyanin-loss phenotype in Arabidopsis tt19, an anthocyanin transport mutant, in the stem and rosette leaves, but not in the mutant seeds (Fig. 10). Similar to our findings, previously identified genes such as PpGST1 from peach (Prunus persica L. (Batsch)) [72], An9 (Petunia hybrida), LcGST4 (Litchi chinensis), strawberry RAP (Fragaria ananassa), CsGSTF1 (Camellia sinensis), and apple MdGSTF6 (Malus domestica) [70, 71, 73,74,75], could functionally complement the anthocyanin-less phenotype of Arabidopsis tt19 mutant, but not the proanthocyanidin-deficient phenotype in the seed coat.

Implications and future directions

The insights gained from this study have implications for both agricultural and nutritional research. The identification of key genes and metabolites involved in pigment accumulation can inform breeding strategies aimed at enhancing the purple coloration in Chinese kale or related crops. Additionally, the health-promoting properties of the secondary metabolites identified warrant further investigation to confirm their bioactivity and potential contributions to human health. Our integrated metabolomic and transcriptomic analysis provides valuable insights into the molecular basis of the purple-colored stalks in Chinese kale. These findings shed light on the regulatory network governing anthocyanin pigment biosynthesis and lay the groundwork for future studies exploring the functional significance of key metabolites in this nutritionally important vegetable.

Conclusion

The discovery of the function of BoGSTF12 in anthocyanin transport and accumulation using a combination of RNA-seq and metabolomic research marks a substantial advancement in our understanding of the molecular processes driving anthocyanin production in plants. This research provides useful knowledge for both horticultural and agricultural applications by dissecting the complex mechanism of anthocyanin production and transport, enabling the development of crops with improved nutritional value and attractive traits. This work also highlights the ability of contemporary -omics technology to decipher the intricacies of plant biology and has the potential to inspire further advancements in farming and plant breeding.

Data availability

All data generated or analyzed during this study are included in supplementary information.

References

  1. Aires A. Brassica composition and food processing. Processing and Impact on active components in Food. Elsevier; 2015. pp. 17–25.

  2. Cartea ME, Francisco M, Soengas P, Velasco P. Phenolic compounds in Brassica vegetables. Molecules. 2010;16(1):251–80.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Ciska E, Martyniak-Przybyszewska B, Kozlowska H. Content of glucosinolates in cruciferous vegetables grown at the same site for two years under different climatic conditions. J Agric Food Chem. 2000;48(7):2862–7.

    Article  CAS  PubMed  Google Scholar 

  4. Jahangir M, Kim HK, Choi YH, Verpoorte R. Health-affecting compounds in Brassicaceae. CRFSFS. 2009;8(2):31–43.

    CAS  Google Scholar 

  5. Paiva SA, Russell RM. β-carotene and other carotenoids as antioxidants. JACN. 1999;18(5):426–33.

    CAS  Google Scholar 

  6. Park WT, Kim JK, Park S, Lee S-W, Li X, Kim YB, Uddin MR, Park NI, Kim S-J, Park SU. Metabolic profiling of glucosinolates, anthocyanins, carotenoids, and other secondary metabolites in kohlrabi (Brassica oleracea var. gongylodes). J Agric Food Chem. 2012;60(33):8111–6.

    Article  CAS  PubMed  Google Scholar 

  7. Wang Y, Wang H, Gao M, Fan Z, Chen Y, Jin Y. Overexpression of kale (Brassica oleracea L. Var. Acephala) BoMYB increases anthocyanin content in Arabidopsis thaliana. Biotechnol Biotechnol Equip. 2019;33(1):902–10.

    Article  CAS  Google Scholar 

  8. Feild TS, Lee DW, Holbrook NM. Why leaves turn red in autumn. The role of anthocyanins in senescing leaves of red-osier dogwood. Plant Physiol. 2001;127(2):566–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Honda T, Saito N. Recent progress in the chemistry of polyacylated anthocyanins as flower color pigments. Heterocycles: Int J Reviews Commun Heterocycl Chem. 2002;56(1–2):633–92.

    Article  CAS  Google Scholar 

  10. Tena N, Martín J, Asuero AG. State of the art of anthocyanins: antioxidant activity, sources, bioavailability, and therapeutic effect in human health. Antioxidants. 2020;9(5):451.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Koes R, Verweij W, Quattrocchio F. Flavonoids: a colorful model for the regulation and evolution of biochemical pathways. Trends Plant Sci. 2005;10(5):236–42.

    Article  CAS  PubMed  Google Scholar 

  12. Iwashina T. Contribution to flower colors of flavonoids including anthocyanins: a review. Nat Prod Commun. 2015;10(3):1934578X1501000335.

    CAS  Google Scholar 

  13. Ramsay NA, Glover BJ. MYB–bHLH–WD40 protein complex and the evolution of cellular diversity. Trends Plant Sci. 2005;10(2):63–70.

    Article  CAS  PubMed  Google Scholar 

  14. Yang Y, Yao G, Yue W, Zhang S, Wu J. Transcriptome profiling reveals differential gene expression in proanthocyanidin biosynthesis associated with red/green skin color mutant of pear (Pyrus communis L). Front Plant Sci. 2015;6:795.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Zhao L, Gao L, Wang H, Chen X, Wang Y, Yang H, Wei C, Wan X, Xia T. The R2R3-MYB, bHLH, WD40, and related transcription factors in flavonoid biosynthesis. Funct Integr Genomic. 2013;13:75–98.

    Article  CAS  Google Scholar 

  16. Chalker-Scott L. Environmental significance of anthocyanins in plant stress responses. Photochem Photobiol. 1999;70(1):1–9.

    Article  CAS  Google Scholar 

  17. Hughes NM, Reinhardt K, Feild TS, Gerardi AR, Smith WK. Association between winter anthocyanin production and drought stress in angiosperm evergreen species. J Exp Bot. 2010;61(6):1699–709.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. An JP, Wang XF, Zhang XW, Xu HF, Bi SQ, You CX, Hao YJ. An apple MYB transcription factor regulates cold tolerance and anthocyanin accumulation and undergoes MIEL1-mediated degradation. Plant Biotechnol J. 2020;18(2):337–53.

    Article  CAS  PubMed  Google Scholar 

  19. Liang J, He J. Protective role of anthocyanins in plants under low nitrogen stress. Biochem Biophys ResComm. 2018;498(4):946–53.

    Article  CAS  Google Scholar 

  20. Mbarki S, Sytar O, Zivcak M, Abdelly C, Cerda A, Brestic M. Anthocyanins of coloured wheat genotypes in specific response to salstress. Molecules. 2018;23(7):1518.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Zhang Q, Zhai J, Shao L, Lin W, Peng C. Accumulation of anthocyanins: an adaptation strategy of Mikania micrantha to low temperature in winter. Front Plant Sci. 2019;10:1049.

    Article  PubMed  PubMed Central  Google Scholar 

  22. He Q, Wu J, Xue Y, Zhao W, Li R, Zhang L. The novel gene BrMYB2, located on chromosome A07, with a short intron 1 controls the purple-head trait of Chinese cabbage (Brassica rapa L). Hortic Res 2020, 7.

  23. He Q, Lu Q, He Y, Wang Y, Zhang N, Zhao W, Zhang L. Dynamic changes of the anthocyanin biosynthesis mechanism during the development of heading Chinese cabbage (Brassica rapa L.) and Arabidopsis under the control of BrMYB2. Front Plant Sci. 2020;11:593766.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Zhang X, Zhang K, Wu J, Guo N, Liang J, Wang X, Cheng F. QTL-Seq and sequence assembly rapidly mapped the gene BrMYBL2. 1 for the purple trait in Brassica rapa. Sci Rep. 2020;10(1):2328.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Fu M, Guo J, Tang K, Jiang S, Luo S, Luo W, Khan I, Li G. Comparative transcriptome analysis of Purple and Green Flowering Chinese Cabbage and functional analyses of BrMYB114 gene. Int J Mol Sci. 2023;24(18):13951.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Li G-H, Chen H-C, Liu J-L, Luo W-L, Xie D-S, Luo S-B, Wu T-Q, Akram W, Zhong Y-J. A high-density genetic map developed by specific-locus amplified fragment (SLAF) sequencing and identification of a locus controlling anthocyanin pigmentation in stalk of Zicaitai (Brassica rapa L. ssp. chinensis var. Purpurea). BMC Genom. 2019;20:1–13.

    Google Scholar 

  27. Zhang J-Y, Huang S-N, Mo Z-H, Xuan J-P, Jia X-D, Wang G, Guo Z-R. De novo transcriptome sequencing and comparative analysis of differentially expressed genes in kiwifruit under waterlogging stress. Mol Breeding. 2015;35:1–12.

    Article  Google Scholar 

  28. Xie Q, Hu Z, Zhang Y, Tian S, Wang Z, Zhao Z, Yang Y, Chen G. Accumulation and molecular regulation of anthocyanin in purple tumorous stem mustard (Brassica juncea var. Tumida Tsen Et Lee). J Agric Food Chem. 2014;62(31):7813–21.

    Article  CAS  PubMed  Google Scholar 

  29. Zhao Z, Xiao L, Xu L, Xing X, Tang G, Du D. Fine mapping the BjPl1 gene for purple leaf color in B2 of Brassica juncea L. through comparative mapping and whole-genome re-sequencing. Euphytica. 2017;213:1–11.

    Article  Google Scholar 

  30. Petroni K, Tonelli C. Recent advances on the regulation of anthocyanin synthesis in reproductive organs. Plant Sci. 2011;181(3):219–29.

    Article  CAS  PubMed  Google Scholar 

  31. Ren J, Fu W, Du J, Hou A, Liu Z, Feng H. Identification of a candidate gene for Re, the factor determining the red leaf phenotype in ornamental kale using fine mapping and transcriptome analysis. Plant Breed. 2017;136(5):738–48.

    Article  CAS  Google Scholar 

  32. Ren J, Liu Z, Niu R, Feng H. Mapping of Re, a gene conferring the red leaf trait in ornamental kale (Brassica oleracea L. var. Acephala). Plant Breed. 2015;134(4):494–500.

    Article  CAS  Google Scholar 

  33. Liu C, Yao X, Li G, Huang L, Wu X, Xie Z. Identification of major loci and candidate genes for anthocyanin biosynthesis in broccoli using QTL-Seq. Horticulturae. 2021;7(8):246.

    Article  CAS  Google Scholar 

  34. Feng X, Zhang Y, Wang H, Tian Z, Xin S, Zhu P. The dihydroflavonol 4-reductase BoDFR1 drives anthocyanin accumulation in pink-leaved ornamental kale. Theor Appl Genet. 2021;134:159–69.

    Article  CAS  PubMed  Google Scholar 

  35. Yan C, An G, Zhu T, Zhang W, Zhang L, Peng L, Chen J, Kuang H. Independent activation of the BoMYB2 gene leading to purple traits in Brassica oleracea. Theor Appl Genet. 2019;132:895–906.

    Article  CAS  PubMed  Google Scholar 

  36. Wu S, Lei J, Chen G, Chen H, Cao B, Chen C. De novo transcriptome assembly of Chinese kale and global expression analysis of genes involved in glucosinolate metabolism in multiple tissues. Front Plant Sci. 2017;8:92.

    PubMed  PubMed Central  Google Scholar 

  37. Rabino I, Mancinelli AL. Light, temperature, and anthocyanin production. Plant Physiol. 1986;81(3):922–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Li S. Transcriptional control of flavonoid biosynthesis: fine-tuning of the MYB-bHLH-WD40 (MBW) complex. Plant Signal Behav. 2014;9(1):e27522.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Xie Y, Tan H, Ma Z, Huang J. DELLA proteins promote anthocyanin biosynthesis via sequestering MYBL2 and JAZ suppressors of the MYB/bHLH/WD40 complex in Arabidopsis thaliana. Mol Plant. 2016;9(5):711–21.

    Article  CAS  PubMed  Google Scholar 

  40. Matsui K, Oshima Y, Mitsuda N, Sakamoto S, Nishiba Y, Walker AR, Ohme-Takagi M, Robinson SP, Yasui Y, Mori M. Buckwheat R2R3 MYB transcription factor FeMYBF1 regulates flavonol biosynthesis. Plant Sci. 2018;274:466–75.

    Article  CAS  PubMed  Google Scholar 

  41. Liu T, Song S, Yuan Y, Wu D, Chen M, Sun Q, Zhang B, Xu C, Chen K. Improved peach peel color development by fruit bagging. Enhanced expression of anthocyanin biosynthetic and regulatory genes using white non-woven polypropylene as replacement for yellow paper. Sci Hort. 2015;184:142–8.

    Article  CAS  Google Scholar 

  42. Liang W, Ni L, Carballar-Lejarazú R, Zou X, Sun W, Wu L, Yuan X, Mao Y, Huang W, Zou S. Comparative transcriptome among Euscaphis Konishii Hayata tissues and analysis of genes involved in flavonoid biosynthesis and accumulation. BMC Genom. 2019;20(1):1–14.

    Article  Google Scholar 

  43. Baudry A, Caboche M, Lepiniec L. TT8 controls its own expression in a feedback regulation involving TTG1 and homologous MYB and bHLH factors, allowing a strong and cell-specific accumulation of flavonoids in Arabidopsis thaliana. Plant J. 2006;46(5):768–79.

    Article  CAS  PubMed  Google Scholar 

  44. Abdel-Aal E-SM, Young JC, Rabalski I. Anthocyanin composition in black, blue, pink, purple, and red cereal grains. J Agric Food Chem. 2006;54(13):4696–704.

    Article  CAS  PubMed  Google Scholar 

  45. Zhu H-F, Fitzsimmons K, Khandelwal A, Kranz RG. CPC, a single-repeat R3 MYB, is a negative regulator of anthocyanin biosynthesis in Arabidopsis. Mol Plant. 2009;2(4):790–802.

    Article  CAS  PubMed  Google Scholar 

  46. Jin S-W, Rahim MA, Jung H-J, Afrin KS, Kim H-T, Park J-I, Kang J-G, Nou I-S. Abscisic acid and ethylene biosynthesis-related genes are associated with anthocyanin accumulation in purple ornamental cabbage (Brassica oleracea var. acephala). Genome. 2019;62(8):513–26.

    Article  CAS  PubMed  Google Scholar 

  47. Park CH, Bong SJ, Lim CJ, Kim JK, Park SU. Transcriptome analysis and metabolic profiling of green and red mizuna (Brassica rapa L. var. Japonica). Foods. 2020;9(8):1079.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Chiu L-W, Zhou X, Burke S, Wu X, Prior RL, Li L. The purple cauliflower arises from activation of a MYB transcription factor. Plant Physiol. 2010;154(3):1470–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Charron CS, Clevidence BA, Britz SJ, Novotny JA. Effect of dose size on bioavailability of acylated and nonacylated anthocyanins from red cabbage (Brassica oleracea L. Var. capitata). J Agric Food Chem. 2007;55(13):5354–62.

    Article  CAS  PubMed  Google Scholar 

  50. Yuan Y, Chiu L-W, Li L. Transcriptional regulation of anthocyanin biosynthesis in red cabbage. Planta. 2009;230:1141–53.

    Article  CAS  PubMed  Google Scholar 

  51. Zhang H, Gong J, Chen K, Yao W, Zhang B, Wang J, Tian S, Liu H, Wang Y, Liu Y. A novel R3 MYB transcriptional repressor, MaMYBx, finely regulates anthocyanin biosynthesis in grape hyacinth. Plant Sci. 2020;298:110588.

    Article  CAS  PubMed  Google Scholar 

  52. Junqing W, Jing Z, Meiling Q, Yanjing R, Zhang H, Zihui D, Lingyu H, ZHANG L. Genetic analysis and mapping of the purple gene in purple heading Chinese cabbage. Hortic Plant J. 2016;2(6):351–6.

    Article  Google Scholar 

  53. Saito K, Yonekura-Sakakibara K, Nakabayashi R, Higashi Y, Yamazaki M, Tohge T, Fernie AR. The flavonoid biosynthetic pathway in Arabidopsis: structural and genetic diversity. Plant Physiol Biochem. 2013;72:21–34.

    Article  CAS  PubMed  Google Scholar 

  54. Zhao ZC, Hu GB, Hu FC, Wang HC, Yang ZY, Lai B. The UDP glucose: flavonoid-3-O-glucosyltransferase (UFGT) gene regulates anthocyanin biosynthesis in litchi (Litchi Chinesis Sonn.) During fruit coloration. Mol Biol Rep. 2012;39:6409–15.

    Article  CAS  PubMed  Google Scholar 

  55. Offen W, Martinez-Fleites C, Yang M, Kiat‐Lim E, Davis BG, Tarling CA, Ford CM, Bowles DJ, Davies GJ. Structure of a flavonoid glucosyltransferase reveals the basis for plant natural product modification. EMBO J. 2006;25(6):1396–405.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Hiromoto T, Honjo E, Noda N, Tamada T, Kazuma K, Suzuki M, Blaber M, Kuroki R. Structural basis for acceptor-substrate recognition of UDP‐glucose: anthocyanidin 3‐O‐glucosyltransferase from Clitoria ternatea. Protein Sci. 2015;24(3):395–407.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Li Y, Shan X, Tong L, Wei C, Lu K, Li S, Kimani S, Wang S, Wang L, Gao X. The conserved and particular roles of the R2R3-MYB regulator FhPAP1 from Freesia hybrida in flower anthocyanin biosynthesis. Plant Cell Physiol. 2020;61(7):1365–80.

    Article  CAS  PubMed  Google Scholar 

  58. Ni J, Premathilake AT, Gao Y, Yu W, Tao R, Teng Y, Bai S. Ethylene-activated PpERF105 induces the expression of the repressor‐type R2R3‐MYB gene PpMYB140 to inhibit anthocyanin biosynthesis in red pear fruit. Plant J. 2021;105(1):167–81.

    Article  CAS  PubMed  Google Scholar 

  59. Gonzalez A, Zhao M, Leavitt JM, Lloyd AM. Regulation of the anthocyanin biosynthetic pathway by the TTG1/bHLH/Myb transcriptional complex in Arabidopsis seedlings. Plant J. 2008;53(5):814–27.

    Article  CAS  PubMed  Google Scholar 

  60. Zhang F, Gonzalez A, Zhao M, Payne CT, Lloyd A. A network of redundant bHLH proteins functions in all TTG1-dependent pathways of Arabidopsis. 2003.

  61. Tohge T, Nishiyama Y, Hirai MY, Yano M, Nakajima Ji, Awazuhara M, Inoue E, Takahashi H, Goodenowe DB, Kitayama M. Functional genomics by integrated analysis of metabolome and transcriptome of Arabidopsis plants over-expressing an MYB transcription factor. Plant J. 2005;42(2):218–35.

    Article  CAS  PubMed  Google Scholar 

  62. Guo N, Wu J, Zheng S, Cheng F, Liu B, Liang J, Cui Y, Wang X. Anthocyanin profile characterization and quantitative trait locus mapping in zicaitai (Brassica rapa L. ssp. chinensis var. Purpurea). Mol Breeding. 2015;35:1–11.

    Article  CAS  Google Scholar 

  63. Zhao J. Flavonoid transport mechanisms: how to go, and with whom. Trends Plant Sci. 2015;20(9):576–85.

    Article  CAS  PubMed  Google Scholar 

  64. Marrs KA, Alfenito MR, Lloyd AM, Walbot V. A glutathione S-transferase involved in vacuolar transfer encoded by the maize gene Bronze-2. Nature. 1995;375(6530):397–400.

    Article  CAS  PubMed  Google Scholar 

  65. Mueller LA, Goodman CD, Silady RA, Walbot V. AN9, a petunia glutathione S-transferase required for anthocyanin sequestration, is a flavonoid-binding protein. Plant Physiol. 2000;123(4):1561–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Sun Y, Li H, Huang J-R. Arabidopsis TT19 functions as a carrier to transport anthocyanin from the cytosol to tonoplasts. Mol Plant. 2012;5(2):387–400.

    Article  CAS  PubMed  Google Scholar 

  67. Kitamura S, Akita Y, Ishizaka H, Narumi I, Tanaka A. Molecular characterization of an anthocyanin-related glutathione S-transferase gene in cyclamen. J Plant Physiol. 2012;169(6):636–42.

    Article  CAS  PubMed  Google Scholar 

  68. Yamazaki M, Shibata M, Nishiyama Y, Springob K, Kitayama M, Shimada N, Aoki T, Ayabe Si, Saito K. Differential gene expression profiles of red and green forms of Perilla frutescens leading to comprehensive identification of anthocyanin biosynthetic genes. FEBS J. 2008;275(13):3494–502.

    Article  CAS  PubMed  Google Scholar 

  69. Conn S, Curtin C, Bézier A, Franco C, Zhang W. Purification, molecular cloning, and characterization of glutathione S-transferases (GSTs) from pigmented Vitis vinifera L. cell suspension cultures as putative anthocyanin transport proteins. J Exp Bot. 2008;59(13):3621–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Jiang S, Chen M, He N, Chen X, Wang N, Sun Q, Zhang T, Xu H, Fang H, Wang Y. MdGSTF6, activated by MdMYB1, plays an essential role in anthocyanin accumulation in apple. Hortic Res 2019, 6.

  71. Hu B, Zhao J, Lai B, Qin Y, Wang H, Hu G. LcGST4 is an anthocyanin-related glutathione S-transferase gene in Litchi chinensis Sonn. Plant Cell Rep. 2016;35:831–43.

    Article  CAS  PubMed  Google Scholar 

  72. Zhao Y, Dong W, Zhu Y, Allan AC, Lin-Wang K, Xu C. PpGST1, an anthocyanin‐related glutathione S‐transferase gene, is essential for fruit coloration in peach. Plant Biotechnol J. 2020;18(5):1284–95.

    Article  CAS  PubMed  Google Scholar 

  73. Kitamura S, Shikazono N, Tanaka A. TRANSPARENT TESTA 19 is involved in the accumulation of both anthocyanins and proanthocyanidins in Arabidopsis. Plant J. 2004;37(1):104–14.

    Article  CAS  PubMed  Google Scholar 

  74. Luo H, Dai C, Li Y, Feng J, Liu Z, Kang C. Reduced anthocyanins in Petioles codes for a GST anthocyanin transporter that is essential for the foliage and fruit coloration in strawberry. J Exp Bot. 2018;69(10):2595–608.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  75. Wei K, Wang L, Zhang Y, Ruan L, Li H, Wu L, Xu L, Zhang C, Zhou X, Cheng H. A coupled role for CsMYB75 and CsGSTF1 in anthocyanin hyperaccumulation in purple tea. Plant J. 2019;97(5):825–40.

    Article  CAS  PubMed  Google Scholar 

Download references

Funding

This work was financially supported by the 2022 Seed Industry Revitalization Project of Rural Revitalization Strategy Special Fund of Guangdong Province {Yue Cai Nong (2022) No. 184}, 2022 Provincial Rural Revitalization Strategy Special Fund Seed Industry Vitalization Project (2022-NJS-03-001), Guangzhou Science and Technology Planning Project (2023B03J1270).

Author information

Author notes

  1. Kang Tang and Umer Karamat contributed equally to this work.

Authors and Affiliations

  1. College of Horticulture, South China Agricultural University, Guangzhou, 510642, China

    Kang Tang, Shizheng Jiang & Xian Yang

  2. Guangdong Key Laboratory for New Technology Research of Vegetables, Vegetable Research Institute, Guangdong Academy of Agricultural Sciences, Guangzhou, 510642, China

    Kang Tang, Umer Karamat, Guihua Li, Juxian Guo, Shizheng Jiang & Mei Fu

Authors
  1. Kang Tang
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Umer Karamat
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Guihua Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Juxian Guo
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Shizheng Jiang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Mei Fu
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Xian Yang
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Conceived and designed the research: MF and XY. Conducted the experiments: SZJ, KT, and UK. Analyzed the data: GHL and JXG. Contributed to writing and edited the manuscript: KT and UK. All the authors have read and approved the paper.

Corresponding authors

Correspondence to Mei Fu or Xian Yang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1: Table S1:

List of primers

Supplementary Material 2: Table S2:

Sequences used for phylogenetic tree analysis

Supplementary Material 3: Table S3:

Metabolites identifies in both Chinese kale (GREEN and RED)

Supplementary Material 4: Table S4:

Sequencing data statistics

Supplementary Material 5: Table S5:

Mapping rate of sequencing data

Supplementary Material 6: Table S6:

List of DEGs in Reg and green Chinese kale

Supplementary Material 7: Table S7:

List of all significant GO enrichment analyses

Supplementary Material 8: Table S8:

List of all KEGG enrichment pathways

Supplementary Material 9: Table S9:

List of anthocyanin related TFs

Supplementary Material 10: Table S10:

Correlation analysis between DEGs and metabolites

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Tang, K., Karamat, U., Li, G. et al. Integrated metabolome and transcriptome analyses reveal the role of BoGSTF12 in anthocyanin accumulation in Chinese kale (Brassica oleracea var. alboglabra). BMC Plant Biol 24, 335 (2024). https://doi.org/10.1186/s12870-024-05016-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-024-05016-5

Keywords

BMC Plant Biology

ISSN: 1471-2229

Contact us