- Research article
- Open Access
Integrated metabolomics and transcriptome analysis on flavonoid biosynthesis in safflower (Carthamus tinctorius L.) under MeJA treatment
BMC Plant Biology volume 20, Article number: 353 (2020)
Safflower (Carthamus tinctorius L.) is an important cash crop, of which the dried tube flower is not only an important raw material for dyes and cosmetics but also an important herb widely used in traditional Chinese medicine. The pigment and bioactive compounds are composed of flavonoids (mainly quinone chalcones), and studies have reported that MeJA can promote the biosynthesis of quinone chalcones, but the mechanism underlying the effect of MeJA in safflower remains unclear. Here, we attempt to use metabolomics and transcriptome technologies to analyse the molecular mechanism of flavonoid biosynthesis under MeJA treatment in safflower.
Based on a UHPLC-ESI-MS/MS detection platform and a self-built database (including hydroxysafflor yellow A, HSYA), a total of 209 flavonoid metabolites were detected, and 35 metabolites were significantly different after treatment with MeJA. Among them, 24 metabolites were upregulated upon MeJA treatment, especially HSYA. Eleven metabolites were downregulated after MeJA treatment. Integrated metabolomics and transcriptome analysis showed that MeJA might upregulate the expression of upstream genes in the flavonoid biosynthesis pathway (such as CHSs, CHIs and HCTs) and downregulate the expression of downstream genes (such as F3Ms, ANRs and ANSs), thus promoting the biosynthesis of quinone chalcones, such as HSYA. The transcription expressions of these genes were validated by real-time PCR. In addition, the promoters of two genes (CtCHI and CtHCT) that were significantly upregulated under MeJA treatment were cloned and analysed. 7 and 3 MeJA response elements were found in the promoters, respectively.
MeJA might upregulate the expression of the upstream genes in the flavonoid biosynthesis pathway and downregulate the expression of the downstream genes, thus promoting the biosynthesis of quinone chalcones. Our results provide insights and basic data for the molecular mechanism analysis of flavonoid synthesis in safflower under MeJA treatment.
Safflower, Carthamus tinctorius L., is a member of the Asteraceae family and is an important economic plant worldwide. Its dried tubular flowers are an important raw material for dyes and cosmetics and are also an important herb widely used in traditional Chinese medicine. As a traditional Chinese medicine, the dried tubular flowers of safflower have been widely used to improve cerebral blood flow and to treat coronary heart disease, hypertension, and cerebrovascular diseases [1, 2]. Flavonoids, particularly quinone chalcones, are associated with these therapeutic effects. Among them, HSYA is the primary active component. It has the antioxidant activities and the myocardial and cerebral protective effects [3,4,5].
The flavonoid biosynthesis pathway is well understood, especially in some model plants, such as Arabidopsis [6, 7]. Chalcone synthase (CHS) catalyzes the reaction of one molecule of 4-coumaroyl-CoA and three molecules of malonyl-CoA into naringenin chalcone which is required for the flavonoids biosynthesis. It can be hydroxylated or glycosylated to the production of carthamone or HSYA. Chalcone isomerase (CHI) can convert naringenin chalcone into naringenin. In safflower, naringenin can be hydroxylated at the 6 position and glycosylated at the 6 and 7 positions to produce 5,6,7,4′-tetrahydroxyflavanone-6,7-di-O-β-D-glucoside, or it can be hydroxylated at the 6 position and glycosylated at the 5 position to the production of 5,6,7,4′-tetrahydroxyflavanone-5-O-β-D-glucoside (neocarthamin). Additionally, flavanone 3-hydroxylase (F3H) can convert naringenin into dihydrokaempferol. It can be subsequently converted into kaempferol, and kaempferol can then be converted into quercetin . With an emphasis on safflower, the research on flavonoid synthesis has gradually increased. To date, many flavonoid biosynthesis genes have been cloned in safflower, such as chalcone isomerase genes (CHIs) , chalcone synthase genes (CHSs) [10, 11], flavanone 3-hydroxylase genes (F3Hs) , UDP-glucuronosyltransferases genes (UGTs) , and shikimate/quinate hydroxycinnamoyltransferase genes (HCTs) .
As a well-known exogenous inducing factor, methyl jasmonate (MeJA) participates in many plant processes, ranging from plant defence to growth and development . MeJA is of particular interest in plant cell engineering for producing bioactive compounds [16, 17]. It has been reported that flavonoids (mainly quinone chalcones) in safflower can be stimulated under MeJA treatment [13, 18]. However, the molecular mechanism is largely unknown.
Recently, the application of metabolomics to medicinal plants has significantly facilitated the identification of the metabolic pathways of active medicinal compounds in plants. The UHPLC-ESI-MS/MS-based, widely targeted metabolomics method has become very popular in the field of analysis and identification of plant metabolites due to the advantages of high throughput, fast separation, high sensitivity, and wide coverage. Our methodology was based on a multiple reaction monitoring (MRM) approach , with a self-built compounds database according to the results of Chen et al. . This approach has been widely applied in plant metabolite analysis in many plants, such as rice , tomato [21, 22], maize . Furthermore, the integration of transcriptomics and metabolomics has larger advantages in revealing the biosynthetic mechanisms of key metabolic pathways [24,25,26]. Therefore, it is feasible to analyse the flavonoid biosynthesis pathway under MeJA treatment by the two technologies.
Here, metabolic profiling and differential flavonoid metabolites were screened based on a UHPLC-ESI-MS/MS detection platform and a self-built database (including HSYA). In addition, transcriptome sequencing and differential transcripts were analysed. Integrated metabolomics and transcriptome sequencing was analysed based on the KEGG pathway, and the expression of different flavonoid biosynthesis genes with or without MeJA treatment were analysed by real-time PCR. The promoters of genes that were significantly upregulated under MeJA treatment were cloned and analysed. Our results provide insights and basic data for the regulation mechanism analysis of flavonoid synthesis under MeJA treatment.
Metabolic profiling and differential flavonoid metabolite analysis
The flavonoid metabolites in safflower with and without MeJA treatment were investigated based on UHPLC-ESI-MS/MS and a self-built database (including HSYA). A total of 209 flavonoid metabolites were detected, including 62 flavones, 42 flavone C-glycosides, 40 flavonols, 20 flavanones, 18 anthocyanins, 11 isoflavones, 12 flavanols, 2 flavonolignans, 1 quinone chalcone, and 1 alkaloid. (Supplementary Table 1).
Orthogonal partial least squares-discriminant analysis (OPLS-DA) was used for variables with less correlation. As the experiment had biological duplication, the fold change and VIP value of the OPLS-DA model were combined to screen differential metabolites. There were 35 significantly different flavonoid metabolites between MeJA-treated and untreated materials. Among them, 24 metabolites were upregulated upon MeJA treatment, especially hydroxysafflor yellow A (HSYA). Eleven metabolites were downregulated after MeJA treatment. (Fig. 1a,b) The details of the different flavonoid metabolites were provided in Supplementary Table 2. The differential flavonoid metabolites from each comparison group were screened by using the KEGG database. The KEGG classification and enrichment analysis indicated that the differential flavonoid metabolites were mainly involved in the pathway of flavonoid biosynthesis (ko00941). (Fig. 1c).
Transcriptome sequencing and differential transcript analysis
The transcriptomes of the mixed safflower samples were sequenced. In the untreated samples,a 6.50 G clean base was obtained (44,185,480 raw reads and 43,921,238 clean reads), and Q30 was 92.40%. In the MeJA-treated samples, a 6.17G clean base was obtained (42,087,290 raw reads and 41,720,790 clean reads), and Q30 was 91.89%. The different expressed genes were analysed with DESeq2. The total number of differentially expressed genes was 31,822 (20,741 upregulated genes and 11,081 downregulated genes). (Fig. 2).
The differential flavonoid metabolite transcripts from each comparison group were screened by using the KEGG database. The experiments focused on annotations in the phenylpropane metabolic pathway, including flavonoid biosynthesis (ko00941), anthocyanin biosynthesis (ko00942), isoflavonoid biosynthesis (ko00943), and flavone and flavonol biosynthesis (ko00944). The results showed that there were 61 significantly different flavonoid biosynthesis transcripts between MeJA-treated and untreated materials (43 in ko00941, 8 in ko00942, 5 in ko00943, and 5 ko00944). Most of transcripts were involved in flavonoid biosynthesis (ko00941). The details could be viewed in Supplementary Table S3.
Integrated analysis of the transcriptome and metabolome on the KEGG pathway
The metabolic components were mapped onto the pathway of flavonoid metabolism (including ko00941, ko00942, ko00943 and ko00944) by combining metabolic components that had been detected by UHPLC-ESI-MS/MS, thus constructing the metabolic pathway map of integrative flavonoid biosynthesis in safflower. (Fig. 3) At the same time, the differentially annotated metabolites, together with the differentially annotated genes, were indicated on the integrated metabolic map. (Fig. 3) From the integrated metabolic map, it was shown that MeJA might upregulate the expression of the upstream genes in the flavonoid biosynthesis pathway (such as CHSs, CHIs, and HCTs) and might downregulate the expression of downstream genes (such as F3Ms, ANRs, and ANSs), thus promoting the biosynthesis of quinone chalcones, such as HSYA.
Expression analysis of differentially expressed flavonoid biosynthesis genes by real-time PCR
Real-time PCR was used to validate the transcriptome data. Ten differentially expressed genes were selected, including TRINITY_DN28401_c0_g1 (annotated as CHI), TRINITY_DN36537_c1_g1 (annotated as HCT), TRINITY_DN37574_c1_g2 (annotated as HCT), TRINITY_DN39063_c0_g1 (annotated as FLS), TRINITY_DN43344_c1_g1 (annotated as CHI), TRINITY_DN41734_c0_g5 (annotated as CHS), TRINITY_DN42700_c3_g2 (annotated as F3H), TRINITY_DN43120_c5_g1 (annotated as ANR), TRINITY_DN44094_c0_g1 (annotated as ANS), TRINITY_DN44565_c1_g1 (annotated as F3M). The primers used for this experiment could be found in Supplementary Table 4. Among them, HCT contains several transferase enzymes, which include anthranilate N-hydroxycinnamoyl/benzoyltransferase that catalyses the first committed reaction of phytoalexin biosynthesis. F3H can catalyse synthesis N-terminal flavanone 3-hydroxylase. F3M can catalyse flavonoid, NADPH, H+, and O2 into 3′-hydroxyflavonoid and NADP+, and H2O, which functions like F3H. FLS catalyzes the formation of flavonols from dihydroflavonols. CHS catalyses the reaction of one molecule of 4-coumaroyl-CoA and three molecules of malonyl-CoA to form tetrahydroxychalcone. Chalcone isomerase converts tetrahydroxychalcone into naringenin. ANR and ANS are all involved in the synthesis of anthocyanins. The results showed that the expressions of 6 genes, including TRINITY_DN28401_c0_g1 (CHI), TRINITY_DN36537_c1_g1 (HCT), TRINITY_DN37574_c1_g2 (HCT), TRINITY_DN39063_c0_g1 (FLS), TRINITY_DN43344_c1_g1 (CHI) and TRINITY_DN41734_c0_g5 (CHS), were significantly upregulated, and 4 genes, including TRINITY_DN42700_c3_g2 (F3H), TRINITY_DN43120_c5_g1 (ANR), TRINITY_DN44094_c0_g1 (ANS) and TRINITY_DN44565_c1_g1 (F3M), were significantly downregulated. (Fig. 4) From the integrated metabolic map, it could be seen that the expressions of upstream genes in the phenylpropane metabolic pathway were increased with MeJA treatment. And the results of gene expression also suggested that these genes might be involved in metabolic pathways.
Elements analysis of the promoters differentially expressed flavonoid biosynthesis genes
To analyse how MeJA regulates gene expression, the promoters of differentially expressed flavonoid biosynthesis genes were cloned and analysed. Because there were no reference genome sequences from safflower, single oligonucleotide nested PCR (SON-PCR) was used to clone the promoters . In our study, two promoters of MeJA significantly upregulated genes CtCHI and CtHCT (named pCtCHI and pCtHCT, respectively) were successfully cloned, and the fragments were cloned into the T vector and transformed into DH5a bacteria. A positive bacterial solution was selected for sequencing. The sequences were shown in Fig. 5.
The sequences of the promoters were analysed using PlantCARE (http://bioinformatics.psb.ugent.be/webtools/plantCARE) . The sequences upstream of ATG were submitted to a web analysis program. Seven MeJA response elements could be found in the promoters of pCtCHI, while three MeJA response elements could be found in the promoters of pCtHCT. (Fig. 5) The results implied that MeJA could regulate these gene expressions through these elements, thus promoting the synthesis of quinone chalcones.
Many secondary metabolites can be induced by MeJA, such as volatile, stilbene, carotenoids, unsaturated fatty acids, flavonoids, lycopene, among others [29,30,31]. Flavonoids are one of the most reported compounds among them. However, not all the results showed that flavonoids biosynthesis can be upregulated by MeJA. Farag made the analysis of Erythrina lysistemon cell suspension culture in response to MeJA elicitation, and results revealed that triterpene i.e. oleanolic acid and fatty acid i.e. hydroxy-octadecadienoic acid were elicited in response to MeJA, whereas pterocarpans i.e. isoneorautenol showed a decline in response to MeJA elicitation . Our results showed that 24 components of 35 differentially metabolized flavonoids could be induced and 11 components were reduced.
It have been reported that 104 compounds from safflower were isolated and detected . In previous research , 12 flavonoid compounds were detected, including rutin, HSYA, kaempferol-3-O-β-D-rutinoside, kaempferol-3-O-β-D-glucoside, carthamin, luteolin, and quercetin-3-O-β-D-glucoside. We found that rutin, HSYA, kaempferol-3-O-β-D-rutinoside, carthamin, and quercetin-3-O-β-D-glucoside were induced by MeJA treatment. In our results, 209 flavonoid compounds were detected by UHPLC-ESI-MS/MS with a self-built database. We profiled the chemical compositions of safflower in a systematic and comprehensive way, which provided a reference for their sufficient utilization in the future. UHPLC-ESI-MS/MS was an effective method to thoroughly understand plant secondary metabolites.
In previous reports of MeJA regulating flavonoids biosynthesis, many studies focused on the analysis of the genes expression. In Arabidopsis, MeJA similarly induced expression of almost all anthocyanin biosynthetic genes (PAL, CHS, CHI, F3H, F3’H, DFR and UFGT), but PAL, CHS, CHI and F3H were at only low levels. Further analysis show that the late anthocyanin biosynthesis genes (such as DFR and UFGT), were found to be up-regulated strongly by MeJA . While in Gynura bicolor DC., the expression of flavonoid biosynthesis genes, GbCHS, GbCHI, GbDFR and GbANS, was markedly up-regulated. Compared with that in Arabidopsis, the genes, which are classified as up flavonoid biosynthesis genes, were found to be up-regulated strongly by MeJA . From the studies above, it can be indicated that the flavonoid biosynthesis genes expression varied among different plants in response to the MeJA treatment. In our study, MeJA could upregulate the expression of upstream genes in the flavonoid biosynthesis pathway (CHSs, CHIs, and HCTs) and downregulate the expression of downstream genes (F3Ms, ANRs, and ANSs), thus promoting the biosynthesis of quinone chalcones, such as HSYA. It is probably that MeJA could upregulate the expression of TRINITY_DN28401_c0_g1 (CHI), TRINITY_DN36537_c1_g1 (HCT), TRINITY_DN37574_c1_g2 (HCT), TRINITY_DN43344_c1_g1 (CHI), and TRINITY_DN41734_c0_g5 (CHS), and downregulate the expression of TRINITY_DN42700_c3_g2 (F3H), TRINITY_DN43120_c5_g1 (ANR), TRINITY_DN44094_c0_g1 (ANS), and TRINITY_DN44565_c1_g1 (F3M), thus promoting the biosynthesis of quinone chalcones (Fig. 4). Of course, to make more conclusive evidence, some other experiments, such as enzyme activity analysis and transgenic experiments, are required.
In the research of Gynura bicolor DC., some genes were markedly up-regulated under MeJA treatment, such as GbCHS, which is expressed more than 100-fold as compared with that before MJ treatment . However, there was no such high change in our results. The main reason might be the different way of dealing with MeJA. In previous research, MeJA was added into 20 ml fresh MS liquid medium, where Gynura bicolor DC. was planted, while in our research, MeJA was sprayed onto healthy safflower flowers. The amounts for plant exposed to MeJA were different.
Since there is no reference genome in safflower, it is still difficult to clone the fragment. However, we fortunately cloned two promoters of up-regulated gene (HCT and CHI) in our experiment. The sequence analysis results showed that there were MeJA response elements on the promoters, which further proved the RT-PCR results. Athough our results showed that MeJA upregulate the expression of upstream genes in the flavonoid biosynthesis pathway and downregulate the expression of downstream genes, thus promoting the biosynthesis of quinone chalcones, such as HSYA, what genes MeJA regulated to promote the production of HSYA is still unknown, as the genes involved in the biosynthesis of HSYA have not been fully identified. Future researches can be strengthened in this domain.
Here, we used metabolomics and transcriptome technologies to analyse the molecular mechanism of flavonoid biosynthesis under MeJA treatment in safflower. Based on a UHPLC-ESI-MS/MS detection platform and a self-built database (including HSYA), a total of 209 flavonoid metabolites were detected, and 35 metabolites were significantly different. Among them, 24 metabolites were upregulated upon MeJA treatment, especially HSYA. Eleven metabolites were downregulated after MeJA treatment. Integrated metabolomics and transcriptome analysis showed that MeJA might upregulate the expression of upstream genes in the flavonoid biosynthesis pathway (such as CHSs, CHIs and HCTs) and downregulate the expression of downstream genes (such as F3Ms, ANRs and ANSs), thus promoting the biosynthesis of quinone chalcones, such as HSYA. The transcription expressions of these genes were validated by real-time PCR. In addition, the promoters of two genes that were significantly upregulated under MeJA treatment were cloned and analysed. Ten MeJA response elements were found in the promoters. Our results provide insights and basic data for the molecular mechanism analysis of flavonoid synthesis in safflower under MeJA treatment.
Safflower used in this experiment was named as “Chuanhonghua No.1”, which is cultivated by Industrial Crop Research Institute, Sichuan Academy of Agricultural Sciences. It was presented by Renchuan Yao and identified as safflower (Carthamus tinctorius L.) by professor Pei Jin. It was cultivated at the medicinal botanical garden on the Wenjiang Campus of Chengdu University of Traditional Chinese Medicine. The treatment was primarily applied according to the previous report with some modifications . A 100 μM solution of MeJA (Sigma-Aldrich, Switzerland) was sprayed onto healthy safflower flowers 3 days after anthesis (DAA). In the control group, the flowers were sprayed with the same solution but without MeJA. The flowers were then enclosed in clear plastic bags to prevent the emission of volatile phytohormones and to allow for the elicit or solutions to be more highly absorbed. After treatment for 6 h, the plastic bags were removed, and samples of flowers were collected, frozen immediately in liquid nitrogen and stored in a freezer at − 80 °C. For the RNA sequencing, five inflorescences of safflower were mixed as a sample. As the RNA sequencing results were verified by the real-time PCR experiments (the experimental method was listed below) in the study, there were no replicates for the RNA sequencing. For metabolism analysis, ten inflorescences of safflower plants were mixed as one sample. Three biological replicates were performed for the metabolomics analysis. The flowers of safflower with or without MeJA treatment are shown in Fig. 6.
Sample preparation and extraction for metabolomic analysis
The freeze-dried flowers were crushed using a mixer mill (MM 400, Retsch) with a zirconia bead for 1.5 min at 30 Hz. One hundred milligrams of powder was weighed and extracted overnight at 4 °C with 1.0 mL 70% aqueous methanol. Following centrifugation at 10000 g for 10 min, the extracts were absorbed (CNWBOND Carbon-GCB SPE Cartridge, 250 mg, 3 mL; ANPEL, Shanghai, China, www.anpel.com.cn/cnw) and filtered (SCAA-104, 0.22 μm pore size; ANPEL, before LC-MS analysis.
The sample extracts were analysed using an LC-ESI-MS/MS system (HPLC, Shim-pack UFLC SHIMADZU CBM30A system; MS, Applied Biosystems 6500 Q TRAP). The analytical conditions were as follows: UHPLC column, Waters ACQUITY UHPLC HSS T3 C18 (1.8 μm, 2.1 mm × 100 mm); solvent system, A water (0.04% acetic acid): B acetonitrile (0.04% acetic acid); gradient program, 0% B at 0 min, 95% B at 11.0 min, 95% B at 12.0 min,5% B at 12.1 min, and 5% B at 15.0 min; flow rate, 0.40 mL/min; temperature, 40 °C; and injection volume, 2 μL. The UHPLC effluent was connected to an ESI-triple quadrupole-linear ion trap (Q TRAP)-MS.
Mass spectrometry followed the method of Chen et al. . LIT and triple quadrupole (QQQ) scans were acquired on a triple quadrupole-linear ion trap mass spectrometer (Q TRAP), API 6500 QTRAP LC/MS/MS system equipped with an ESI Turbo ion-spray interface, operating in positive ion mode and controlled by Analyst 1.6.3 software (AB Sciex, Waltham, MA, USA). The ESI source operation parameters were as follows: ion source, turbo spray; source temperature, 500 °C; ion spray voltage (IS), 5500 V; the ion source gas I (GSI), gas II (GSII), and curtain gas (CUR) were set at 55, 60, and 25.0 psi, respectively; and the collision gas (CAD) was high. Instrument tuning and mass calibration were performed with 10 and 100 μmol/L polypropylene glycol solutions in QQQ modes. QQQ scans were acquired as MRM experiments with collision gas (nitrogen) set to 5 psi. DP and CE for individual MRM transitions was performed with further DP and CE optimization. A specific set of MRM transitions was monitored for each time period according to the metabolites eluted within that period.
Qualitative and quantitative analysis of metabolites
The flavonoid identification and quantification in our study was made according to a method of scheduled multiple reaction monitoring (MRM), which has been previously described . With this methods, Chen made a genome-wide association analyses between genetic and biochemical, and reported 840 metabolites, from which 277 were identified or annotated . Mainly based on this research, a database was built. All the metabolites in our study were detected or annotated by this database. The database and methods have been wildly used to detected or annotate the metabolites, such as the research of Wang  and Liu . This methods is also called as wild-targetd analysis. As HSYA was not in the database, we use its standard sample (No.: MUST-15072815, Chengdu Manster Biotechnology Co., LTD, China) to build a library and add it to the self-built database. The metabolites of the samples were qualitatively and quantitatively analysed by mass spectrometry. The characteristic ions of each substance were filtered by the triple quadrupole, and the signal strength of the characteristic ions was obtained in the detector. Chromatographic peaks were analyzed with Multi Quant software 3.0.3. The integrated area peak of each compound was used for the PCA and OPLS-DA analysis..
Differential metabolite analysis
Principal component analysis (PCA) can effectively extract main variance information and was used in many other research [29, 38]. In addition, orthogonal partial least squares-discriminant analysis (OPLS-DA) was used for variables with less correlation. The experiment had biological duplication; thus, the fold change and VIP value of the OPLS-DA model were combined to screen differential metabolites. PCA was analyzed by R software built-in functions (www.r-project.org/). The parameter: scale = True. After conversion of the original data by log2, the data was centralized (Mean Centering) and analyzed by OPLSR. Anal of Metabo Analyst in R software. The main steps of our study can be referred to a previouly research . Screening criteria: 1. Metabolites must exhibit fold changes≥2 and fold changes≤0.5. If a metabolite in the experimental groupwas more than 2 times or less than 0.5 times that of the metabolite in the control group, the difference was considered significant. 2. If there is biological duplication in the sample grouping, the metabolites of VIP ≥1 are selected on the basis of the above. The VIP value indicates the influence intensity of the difference between groups of corresponding metabolites in the classification and discrimination of each group of samples in the model.
RNA sequencing and annotation
RNA isolation and purification and cDNA library construction and sequencing were performed as previously described . All tissues were ground on dry ice, and total RNA was prepared by using TRIzol reagent (Invitrogen, CA, USA). To remove DNA, an aliquot of total RNA was treated with DNase (Takara, Dalian, China). RNA quantity and quality were determined by using a NanoDrop 2000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA) and an Agilent 2100 Bioanalyser (Agilent Technologies, CA, USA), respectively. mRNA was isolated from total RNA using magnetic beads with oligo (dT) primer; cDNA was synthesized using a cDNA synthesis kit (TaKaRa, Dalian, China) and linking the sequencing adapter to both ends. The library preparations were sequenced on an Illumina HiSeq 4000 platform, and the unigene sequences obtained from our laboratory transcriptome database by RSEM (RNAseq by Expectation-Maximization) software were integrated for annotation. The whole set of transcript data can be found in the National Center for Biotechnology Information (NCBI) SRA database (SRR10011980 and SRR10011979).
Screening of differential genes
DESeq2 [41, 42] was used for differential expression analysis among sample groups. DESeq2 requires the input of unregulated read counting data of genes rather than RPKM (Reads Per Kilobase Million), FPKM (Fragments Per Kilobase Million) and other standardized data. Read counts of genes are the expected_count outputs calculated using RSEM (RNAseq by Expectation-Maximization). The expected count is generally lower than the read number. After discrepancy analysis, multiple hypothesis tests are needed to correct the hypothesis test probability (P value) with the Benjamini-Hochberg procedure to obtain the false discovery rate (FDR) when |log2fold change| ≥ 1and FDR < 0.05.
Real-time PCR analysis
Real-time PCR analysis was performed as previously described . Total RNA was isolated using an RNA extraction kit (Invitrogen, CA, USA), and reverse transcription was carried out using the prime script reagent kit (Takara, Dalian, China). The primers used to amplify the screened genes (CHS, CHI, HCT, etc.) by real-time PCR were designed by Primer 5.0, and parts of the safflower 28S coding region were used as an internal reference gene. The primer details are listed in Supplementary Table 4. All of the primers were tested for their specificity by agarose gel electrophoresis. RT-PCR was performed using an SYBR prime script RT-PCR kit (Takara) with three replicates, and the cycling conditions were set according to the manual. The Bio-Rad CFX96 real-time PCR detection system (Hercules, CA, USA) was used in our experiment.
Promoter cloning and sequence analysis
Because there is no complete genome sequence for safflower, the 5′ flanking region of the flavonoid gene was isolated from genomic DNA of the safflower plant using SON-PCR. The reactions (50 μL) were performed with 200 μmol of each dNTP, 2 μmol of primer and 2 units of LA Taq DNA polymerase (Takara). The primer details are listed in Supplementary Table 4. The resultant PCR fragment, measuring approximately 2 kb, was cloned into a pMD19-T vector (Takara) and sequenced. The flavonoid gene promoter sequence was analysed with the PLACE Web SignalScan program (http://bioinformatics.psb. ugent.be/webtools/plantCARE ;.
Availability of data and materials
The raw data was uploaded to Sequence Read Archive (SRR10011980 and SRR10011979).
Kyoto Encyclopedia of Genes and Genomes
Hydroxysafflor yellow A
China T. S. P. C. o. Pharmacopoeia of the People’s Republic of China, Part I. Beijing: Chemical Industry Press; 2015. p. 151.
Lou ZQ, Liu ML. A pharmacognostical study of the Chinese drug, Honghua,Carthamitinctor II. Acta Pharm Sin. 1956;4:233–9.
Qu C, Wang LY, Jin WT, Tang YP, Jin Y, Shi QX, Shang LL, Shang EX, Duan JA. Comparative analysis of the effects of Hydroxysafflor yellow a and Anhydrosafflor yellow B in safflower series of herb pairs using prep-HPLC and a selective knock-out approach. Molecules. 2016;21:1480–506.
Wei X, Liu H, Sun X, Fu FH, Zhang XM, Wang J, An J, Ding H. Hydroxysafflor yellow a protects rat brains against ischemia-reperfusion injury byantioxidant action. Neurosci Lett. 2005;386:58–62.
Zhu H, Wang Z, Ma C, Tian J, Fu F, Li C, Guo D, Roeder E, Liu K. Neuroprotective effects of hydroxysafflor yellow a: in vivo and in vitrostudies. Planta Med. 2003;69:429–33.
Winkel-Shirley B. Biosynthesis of flavonoids and effects of stress. Curr Opin Plant Bio. 2002;5:218–39.
Zhang Y, Butelli E, Martin C. Engineering anthocyanin biosynthesis in plants. Curr Opin Plant Bio. 2014;19:81–90.
Huang LL, Yang X, Sun P, Wen T, Hu SQ. The first illuminabased de novo transcriptome sequencing and analysis of safflower flowers. PLoS One. 2012;7:e3865.
Ren CX, Tang XH, Chen CP, Chen J, Pei J, WU YY, Wu QH. Cloning and expression analysis of a new chalcone isomerase gene during flowering in safower. Turkish J Bot. 2019;43. https://doi.org/10.3906/bot-1809-25.
Shinozaki J, Kenmoku H, Nihei K, Masuda K, Noji M, Konno K, Asakawa Y, Kazuma K. Cloning and functional analysis of three chalcone synthases from the flowers of safflowers carthamus tinctorius. Nat Prod Commun. 2016;11(6):787–90.
Guo D, Xue Y, Li D, He B, Jia X, Dong X, Guo M. Overexpression of CtCHS1 increases accumulation of Quinochalcone in safflower. Front Plant Sci. 2017;8:1409–22.
Tu YH, Liu F, Guo DD, Fan LJ, Zhu ZX, Xue YR, Gao Y, Guo ML. Molecular characterization of flavanone 3-hydroxylase gene and flavonoid accumulation in two chemotyped safflower lines in response to methyl jasmonate stimulation. BMC Plant Biol. 2016;16:132–44.
Guo DD, Liu F, Tu YH, He BX, Gao Y, Guo ML. Expression patterns of three UGT genes in different chemotype safflower lines and under MeJA stimulus revealed their potential role in flavonoid biosynthesis. PLoS One. 2016;11:e0158159.
Chen J, Tang XH, Ren CX, Wei B, Wu YY, Wu QH, Pei J. Full-length transcriptome sequences and the identification of putative genes for flavonoid biosynthesis in safflower. BMC Genomics. 2018;19(1):548–60.
Wasternack C. Jasmonates: an update on biosynthesis, signal transduction and action in plant stress response, growth and development. Ann Bot. 2007;100:681–97.
Naoumkina MA, Farag MA, Sumner LW, Tang YH, Liu CJ, Dixon RA. Different mechanisms for phytoalexin induction by pathogen and wound signals in Medicago truncatula. P Natl Acad Sci USA. 2007;104:17909–15.
Suzuki H, Reddy MS, Naoumkina M, Aziz N, May GD, Huhman DV, Sumner LW, Blount JW, Mendes P, Dixon RA. Methyl jasmonate and yeast elicitor induce differential transcriptional and metabolic re-programming in cell suspension cultures of the model legume Medicago truncatula. Planta. 2005;220:696–707.
He BX, Xue YR, Tu YH, Gao Y, Guo ML. CtCHS4 induces the accumulation of safflower quinone chalconesin response to methyl jasmonate induction. Acta Pharm Sin. 2018;53:636–54.
Chen W, Gong L, Guo Z, Wang W, Zhang H, Liu X, Yu S, Xiong L, Luo J. A novel integrated method for large-scale detection, identification, and quantification of widely targeted metabolites: application in the study of Rice metabolomics. Mol Plant. 2013;6:1769–80.
Chen W, Gao Y, Xie W, Gong L, Lu K, Wang W, Li Y, Liu X, Zhang H, Dong H, Zhang W, Zhang L, Yu S, Wang G, Lian X, Luo J. Genome-wide association analyses provide genetic and biochemical insights into natural variation in rice metabolism. Nat Genet. 2014;46:714–21.
Cho K, Cho KS, Sohn HB, Ha IJ, Hong SY, Lee H, Kim YM, Nam MH. Network analysis of the metabolome and transcriptome reveals novel regulation of potato pigmentation. J Exp Bot. 2016;67:1519–33.
Wen W, Li D, Li X, Gao Y, Li W, Li H, Liu J, Liu H, Chen W, Luo J, Yan JB. Metabolome-based genome-wide association study of maize kernel leads to novel biochemical insights. Nat Commun. 2014;5:3438–7.
Zhu G, Wang S, Huang Z, Zhang S, Liao Q, Zhang C, Lin T, Qin M, Peng M, Yang C, Cao X, Han X, Wang X, Knaap E, Zhang Z, Cui X, Klee H, Alisdair RF, Jie L, Huang S. Rewiring of the fruit Metabolome in tomato breeding. Cell. 2018;172:249–61.
Dong T, Han R, Yu J, Zhu M, Zhang Y, Gong Y, Li Z. Anthocyanins accumulation and molecular analysis of correlated genes by metabolome and transcriptome in green and purple asparaguses (Asparagus officinalis, L.). Food Chem. 2019;271:18–28.
Li Y, Fang J, Qi X, Lin M, Zhong Y, Sun L, Cui W. Combined analysis of the fruit Metabolome and Transcriptome reveals candidate genes involved in flavonoid biosynthesis in Actinidia arguta. Int J Mol Sci. 2018;19(5):1471–87.
Meng J, Wang B, He G, Wang Y, Tang X. WangS, Ma Y, Fu C, chai G, Zhou G. metabolomics integrated with Transcriptomics reveals redirection of the Phenylpropanoids metabolic flux in Ginkgo biloba. Agric Food Chem. 2019;67(11):3284–91.
Antal Z, Rascle C, Fèvre M, Bruel C. Single oligonucleotide nested PCR: a rapid method for the isolation of genes and their flankingregions from expressed sequence tags. Curr Genet. 2004;46(4):240–6.
Lescot M, Déhais P, Thijs G, Marchal K, Moreau Y, Peer YVD, Rouzé P, Rombauts S. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002;30(1):325–7.
Saccenti E, Hoefsloot HCJ, Smilde AK, Westerhuis JA, Hendriks MMW. Reflections on univariate and multivariate analysis of metabolomics data. Metabolomics. 2014;10:361–74.
Thakur M, Bhattacharya S, Khosla PK, Puri S. Improving production of plant secondary metabolites through biotic and abiotic elicitation. J Appl Res Med. 2018. https://doi.org/10.1016/j.jarmap.2018.11.004.
Shi J, Ma CY, Qi DD, Lv HP, Yang T, Peng QH, Chen ZM, Lin Z. Erratum to: transcriptional responses and flavor volatiles biosynthesis in methyl jasmonate-treated tea leaves. BMC Plant Biol. 2017;136. https://doi.org/10.1186/s12870-017-1076-5.
Farag MA, Mekky H, El-Masry S. Metabolomics driven analysis of Erythrina lysistemon cell suspension culture in response to methyl jasmonate elicitation. J Adv Res. 2016;7(5):681–9.
Zhou X, Tang L, Xu Y, Zhou G, Wang Z. Towards a better understanding of medicinal uses of Carthamustinctorius L. in traditional Chinese medicine: a phytochemical and pharmacological review. J Ethnopharmacol. 2014;151(1):27–43.
Shan XY, Zhang YS, Peng W, Wang ZL, Xie DX. Molecular mechanism for jasmonate-induction of anthocyanin accumulation in Arabidopsis. J Exp Bot. 2009;60(13):3849–60.
Shimizu Y, Maeda K, Kato M, Shimomura K. Co-expression of GbMYB1 and GbMYC1 induces anthocyanin accumulation in roots of cultured Gynura bicolor DC. Plantlet on methyl jasmonate treatment. Plant Physiol Bioch. 2011;49(2):159–67.
Wang A, Li R, Ren L, Gao X, Zhang Y, Ma Z, Ma D, Luo Y. A comparative metabolomics study of flavonoids in sweet potato with different flesh colors (Ipomoea batatas (L.) lam). Food Chem. 2018;260:124–34.
Liu YH, Lv JH, Liu ZH, Wang J, Yang BZ, Chen WC, Ou LJ, Dai XZ, Zhuqing Z, Xuexiao Z. Integrative analysis of metabolome and transcriptome reveals the mechanism of color formation in pepper fruit (Capsicum annuum L.). Food Chem. 2020:306. https://doi.org/10.1016/j.foodchem.2019.125629.
Li Q, Song J. Analysis of widely targeted metabolites of the euhalophyte Suaeda salsa under saline conditions provides new insights into salt tolerance and nutritional value in halophytic species. BMC Plant Biol. 2019;19:388–98.
Ying H, Shi J, Zhang S, Pingcuo G, Wang S, Zhao F, Cui YN, Zeng XL. Transcriptomic and metabolomic profiling provide novel insights into fruit development and flesh coloration in Prunus mira Koehne, a special wild peach species. BMC Plant Biol. 2019;19:463–78.
Wang Z, Cui Y, Vainstein A, Chen S, Ma H. Regulation of fig (Ficuscarica L.) fruit color: Metabolomic and Transcriptomic analyses of the flavonoid biosynthetic pathway. Front Plant Sci. 2017;8:1990–2004.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-Seq data with DESeq2. Genome Bio. 2014;15:550–70.
Varet H, Loraine BG, Coppée JY, Dillies MA. SARTools: a DESeq2- and EdgeR-based R pipeline for comprehensive differential analysis of RNA-Seq data. PLoS One. 2016;11(6):e0157022.
We acknowledge Dr. Hu Zeng and Yuanyuan Qin (Metware Biotechnology Co., Ltd. Wuhan, China) for the data analysis, American Journal Experts (AJE) for the language editing (verification code: B631-386C-E812-C6A9-351B), and Renchuan Yao for providing us the plant materials.
This project is supported by grants from the National Science Foundation of China (81803669, 81573544), which provide support for the test of metabolic and transcriptional, China Postdoctoral Science Foundation (2018M643790XB), which provide support for the design and the writing of the study.
Ethics approval and consent to participate
This study has not directly involved humans, animals or plants.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
List of the 209 metabolites detected in safflower samples and OPLS-DA results. Two hundred nine flavonoid metabolites details were listed in the table. CK represent treatment without MeJA, M-MeJA represent treatment with MeJA. Three biological repeats were made in the experiment. Table S2. Extract of the 35 metabolites significantly different between MeJA and non-treated materials. CK represents treatment without MeJA. Three biological replicates were used in the experiment. Table S3. Details of the transcripts belonging to the flavonoid biosynthesis which were significantly different after MeJa treatment. M-CK represent treatment without MeJA, M-MeJA represent treatment with MeJA. Table S4. List of primers used for RT-PCR and promoter cloning experiments.
About this article
Cite this article
Chen, J., Wang, J., Wang, R. et al. Integrated metabolomics and transcriptome analysis on flavonoid biosynthesis in safflower (Carthamus tinctorius L.) under MeJA treatment. BMC Plant Biol 20, 353 (2020). https://doi.org/10.1186/s12870-020-02554-6
- Flavonoid biosynthesis
- MeJA treatment
- Molecular mechanism