- Research
- Open access
- Published:
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)
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.
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.
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).
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).
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).
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.
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.
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.
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).
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.
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.
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.
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
Aires A. Brassica composition and food processing. Processing and Impact on active components in Food. Elsevier; 2015. pp. 17–25.
Cartea ME, Francisco M, Soengas P, Velasco P. Phenolic compounds in Brassica vegetables. Molecules. 2010;16(1):251–80.
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.
Jahangir M, Kim HK, Choi YH, Verpoorte R. Health-affecting compounds in Brassicaceae. CRFSFS. 2009;8(2):31–43.
Paiva SA, Russell RM. β-carotene and other carotenoids as antioxidants. JACN. 1999;18(5):426–33.
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.
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.
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.
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.
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.
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.
Iwashina T. Contribution to flower colors of flavonoids including anthocyanins: a review. Nat Prod Commun. 2015;10(3):1934578X1501000335.
Ramsay NA, Glover BJ. MYB–bHLH–WD40 protein complex and the evolution of cellular diversity. Trends Plant Sci. 2005;10(2):63–70.
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.
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.
Chalker-Scott L. Environmental significance of anthocyanins in plant stress responses. Photochem Photobiol. 1999;70(1):1–9.
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.
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.
Liang J, He J. Protective role of anthocyanins in plants under low nitrogen stress. Biochem Biophys ResComm. 2018;498(4):946–53.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Petroni K, Tonelli C. Recent advances on the regulation of anthocyanin synthesis in reproductive organs. Plant Sci. 2011;181(3):219–29.
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.
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.
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.
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.
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.
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.
Rabino I, Mancinelli AL. Light, temperature, and anthocyanin production. Plant Physiol. 1986;81(3):922–4.
Li S. Transcriptional control of flavonoid biosynthesis: fine-tuning of the MYB-bHLH-WD40 (MBW) complex. Plant Signal Behav. 2014;9(1):e27522.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Yuan Y, Chiu L-W, Li L. Transcriptional regulation of anthocyanin biosynthesis in red cabbage. Planta. 2009;230:1141–53.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Zhao J. Flavonoid transport mechanisms: how to go, and with whom. Trends Plant Sci. 2015;20(9):576–85.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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
Authors and Affiliations
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
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.
About this article
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
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12870-024-05016-5