Original Research ARTICLE
Influence of Light and Temperature on Gene Expression Leading to Accumulation of Specific Flavonol Glycosides and Hydroxycinnamic Acid Derivatives in Kale (Brassica oleracea var. sabellica)
- Leibniz Institute of Vegetable and Ornamental Crops Großbeeren/Erfurt e.V., Großbeeren, Germany
Light intensity and temperature are very important signals for the regulation of plant growth and development. Plants subjected to less favorable light or temperature conditions often respond with accumulation of secondary metabolites. Some of these metabolites have been identified as bioactive compounds, considered to exert positive effects on human health when consumed regularly. In order to test a typical range of growth parameters for the winter crop Brassica oleracea var. sabellica, plants were grown either at 400 μmol m−2 s−1 or 100 μmol m−2 s−1 at 10°C, or at 400 μmol m−2 s−1 with 5 or 15°C. The higher light intensity overall increased flavonol content of leaves, favoring the main quercetin glycosides, a caffeic acid monoacylated kaempferol triglycoside, and disinapoyl-gentiobiose. The higher temperature mainly increased the hydroxycinnamic acid derivative disinapoyl-gentiobiose, while at lower temperature synthesis is in favor of very complex sinapic acid acylated flavonol tetraglycosides such as kaempferol-3-O-sinapoyl-sophoroside-7-O-diglucoside. A global analysis of light and temperature dependent alterations of gene expression in B. oleracea var. sabellica leaves was performed with the most comprehensive Brassica microarray. When compared to the light experiment much less genes were differentially expressed in kale leaves grown at 5 or 15°C. A structured evaluation of differentially expressed genes revealed the expected enrichment in the functional categories of e.g. protein degradation at different light intensities or phytohormone metabolism at different temperature. Genes of the secondary metabolism namely phenylpropanoids are significantly enriched with both treatments. Thus, the genome of B. oleracea was screened for predicted genes putatively involved in the biosynthesis of flavonoids and hydroxycinnamic acid derivatives. All identified B. oleracea genes were analyzed for their most specific 60-mer oligonucleotides present on the 2 × 105 K format Brassica microarray. Expression differences were correlated to the structure-dependent response of flavonoid glycosides and hydroxycinnamic acid derivatives to alterations in either light or temperature. The altered metabolite accumulation was mainly reflected on gene expression level of core biosynthetic pathway genes and gave further hints to an isoform specific functional specialization.
Kale (Brassica oleracea var. sabellica) is an increasingly popular leafy vegetable of the Brassicaceae, the plant family with many important food crops. Eating these vegetables is supposed to be health-supporting since a number of epidemiological studies have identified a positive association between its consumption and human health (Verkerk et al., 2009). These positive effects have been mainly attributed to secondary metabolites like glucosinolates and flavonoids present in plant tissues (Calderón-Montaño et al., 2011; Maalik et al., 2014). Kale as a typical winter crop is mainly cultivated in Northern and Central Europe, e.g., in Germany as well as in North America. During cultivation kale is exposed to low or even freezing temperature conditions and it also has to grow at low levels of photosynthetically active radiation (PAR).
It is undoubtable that light and temperature are most important signals for development, growth and metabolism. The light environment is perceived by several receptors, controlling acclimation and avoidance responses to limiting or excess light (Chen et al., 2004; Li et al., 2009) thus creating a complex network of light response (Casal, 2013). Upon temperature changes signaling pathways trigger lower temperature response (Yamaguchi-Shinozaki and Shinozaki, 2006) finally contributing to cold adaptation and freezing tolerance (Knight and Knight, 2012). Furthermore, light and temperature signaling networks are integrated in an interconnected system that may be evolutionarily favorable (Franklin et al., 2014).
When subjected to stress plants often respond with accumulation of secondary metabolites, especially flavonoids (Falcone Ferreyra et al., 2012). These metabolites may play a role in adaptation to the environment and in overcoming stress conditions (Ramakrishna and Ravishankar, 2011). Mutant phenotypes can be detected through altered pigmentation and corresponding genes involved in the flavonoid biosynthetic pathway have been characterized in several plant species (Lepiniec et al., 2006). A scheme of the phenylpropanoid pathways and Arabidopsis genes coding for respective enzyme activities are listed in Figure 1.
Figure 1. Phenylpropanoid metabolism and related pathways. Identified Arabidopsis genes coding for respective enzyme activities are listed. 4CL, 4-coumarate-CoA ligase (CL1, At1g51680; 4CL2, At3g21240; 4CL3, At1g65060; 4CL5, At3g21230); BAHD-ATs, acyl-CoA dependent acyltransferase (e.g., AT5MAT, At3g29590); C3H, coumarate 3′-hydroxylase (REF8, At2g40890); C4H, cinnamate 4-hydroxylase (C4H, At2g30490); CAD, cinnamylalcohol dehydrogenase (CAD4, At3g19450; CAD5, At4g34230); CCoAOMT, caffeoyl-CoA O-methyltransferase (CCoAOMT1, At4g34050; CCoAOMT7, At4g26220); CCR, cinnamoyl-CoA reductase (CCR1, At1g15950; CCR2, At1g15950); CHI, chalcone isomerase (TT5, At3g55120); CHS, chalcone synthase (TT4, At5g13930); COMT, flavone 3′-O-methyltransferase (COMT1, At5g54160); DFR, dihydroflavonol 4-reductase (DFR, At5g42800); F3GTs, flavonol 3-O-glucosyltransferase (e.g., UGT73B2, At 4g34135); F3H, flavanone 3-hydroxylase (TT6, At3g51240); F3′H, flavonoid 3′-hydroxylase (TT7, At5g07990); F3′5′H, flavonoid 3′5′-hydroxylase (not present in A. thaliana); F7GTs, flavonol 7-O-glucosyltransferase (e.g., UGT73C6, At2g36790); FAH1, ferulate 5-hydroxylase (FAH1, At4g36220); FLS, flavonol synthase (FLS1, At5g08640; FLS3, At5g63590); FNSII, flavone synthase (not present in A. thaliana); HCT, shikimate O-hydroxycinnamoyltransferase (HCT, At5g48930); PAL, phenylalanine ammonia lyase (PAL1, At2g37040; PAL2, At3g53260); REF1, sinapaldehyde dehydrogenase (ALDH2C4, At3g24503); SCPL-ATs, acylglucose dependent acyltransferases (e.g., SCPL10, AtSAT, At2g23000); SGT1, sinapate 1-glucosyltransferase (UGT84A2, At3g21560; UGT84A1, At4g15480; UGT84A3, At4g15490; UGT84A4, At4g15500); SST, sinapoylglucose O-sinapoyltransferase (SCPL9, At2g23010; SCPL13, At2g22980).
It is known that flavonoids and hydroxycinnamic acid derivatives protect plants from damage caused by high levels of sunlight. Besides acting as shielding components due to their absorption maxima phenolic compounds are also scavengers of reactive oxygen species (ROS) (Edreva, 2005). The reduction of sunlight has been associated with lower concentrations of flavonoids for the perennial woody species Ligustrum vulgare and Phillyrea latifoilia (Agati and Tattini, 2010) as well as for the annual leafy species B. rapa and B. juncea (Fallovo et al., 2011). The effects of different light conditions on the control of flavonoid biosynthesis is mainly exerted through R2R3 MYB transcription factors, regulating the biosynthesis of distinct flavonoids for example in fruits (Zoratti et al., 2014). It was demonstrated that low temperature enhances phenolic compounds and total flavonoids due to an enzymatic repair inhibition combined with higher quantities of ROS (Bilger et al., 2007; Klimov et al., 2008; Moheb et al., 2011; Lavola et al., 2013). Moreover, a decrease in quercetin glycosides and caffeic acid derivatives with a temperature increase of 2°C was shown for willow (Salix myrsinifolia) (Nybakken et al., 2012).
Compared to other leafy vegetables such as rocket salads or lettuce, kale leaves have high concentrations of a large number of naturally occurring flavonol glycosides and hydroxycinnamic acid derivatives (Ferreres et al., 2009; Lin and Harnly, 2009; Olsen et al., 2009a; Schmidt et al., 2010b; Fiol et al., 2012). Up to now, a total number of 71 flavonol glycosides has been identified in kale by HPLC-DAD-MSn. So far, 27 non-acylated, 30 monoacylated and 14 diacylated glycosides have been described based on the flavonol aglycones quercetin, kaempferol and isorhamnetin (Schmidt et al., 2010b). The main flavonol glycosides in kale are non-acylated and monoacylated quercetin and kaempferol glucosides, with the majority of flavonol glucosides being acylated with hydroxycinnamic acids (Schmidt et al., 2010b). Unlike Arabidopsis thaliana containing only non acylated flavonol glycosides (Saito et al., 2013), the acylation pattern makes kale a very useful Brassicaceae model plant to analyze the biosynthesis and function of hydroxycinnamic acid acylated flavonol glycosides. The main free hydroxycinnamic acids derivatives in kale are caffeoylquinic acid (chlorogenic acid), disinapoyl-gentiobiose and sinapoyl-feruloyl-gentiobiose (Zietz et al., 2010; Fiol et al., 2012). This is also very different from A. thaliana that mainly contains sinapoylcholine, sinapoylglucose, and sinapoylmalate as free hydroxycinnamic acid derivatives (Vogt, 2010; Fraser and Chapple, 2011).
Based on field experiments in a typical growing season in Germany multiple regressions with the factors temperature and radiation were calculated for kale's flavonol glycosides and underlined a structure-specific response (Schmidt et al., 2010a; Neugart et al., 2012a). In this field experiment the effect of PAR remained remarkably low thus indicating that the influence of PAR on the investigated structures might be suppressed due to the decreasing temperatures. In a consecutive experiment with young kale plants exposed for 1 week to low temperatures (5 and 10°C) with diverse PAR levels (200 to 800 μmol m−2 s−1) in growth chambers the specific effect of low temperature and moderate PAR levels on structurally different flavonoid glycosides was established further (Neugart et al., 2013). Under these constant conditions the specific impact of PAR on the accumulation of flavonoid glycosides was mainly exerted on the levels of quercetin triglycosides, however, expression of respective genes was unaffected assuming an adaptation response of the rather young plants (Neugart et al., 2013).
The aim of this study was to (I) distinguish the effect of ecologically relevant PAR (400 μmol m−2 s−1 or 100 μmol m−2 s−1 at 10°C) and temperature levels (5 or 15°C at 400 μmol m−2 s−1) on structurally different flavonoid glycosides and hydroxycinnamic acid derivatives in adult kale plants under defined conditions in climate chambers. These conditions are comparable to conditions in the field when kale is cultivated during central European winter thus resembling typical light and temperature conditions in a growing season. We want to (II) further differentiate the plant's response to PAR and temperature on expression level using the most comprehensive 2 × 104 k format Brassica Array. Consequently, the focus of the transcriptome analysis will be directed toward genes putatively involved in the phenylpropanoid pathway.
Materials and Methods
Plant Material and Experimental Design
Kale (B. oleracea var. sabellica) plants of the cultivar Winterbor (by Bruno Nebelung, Norken, Germany) were grown in pots under greenhouse conditions for 12 weeks. For both experiments, the plants were cultivated in potting soil (Standard soil Type 1, Fitz Kausek GmbH & Co.KG, Mittenwalde, Germany). Water supply was given as needed by the plants and fertilization was given as nutrient solution (pH 6.4, P 41 mg L−1, Mg 43 mg L−1, NO3154 mg L−1, Ca 205 mg L−1, K 230 mg L−1, SO4 280 mg L−1, Zn 0.82 mg L−1) every 4 weeks. The adult kale plants were then transferred to climate chambers. The plants were treated either with (I) different PAR, with 100 μmol m−2 s−1 and 400 μmol m−2 s−1, measured with light meter Li-250 (Li-COR, USA), at 10°C or with (II) different temperatures, at 5 and 15°C with 400 μmol m−2 s−1. The PAR radiation was applied from 8 am to 6 pm each day. These conditions are comparable to conditions in the field during cultivation in central European winter season. The temperature was measured with a Pt1000 sensor integrated in the climate chamber with a fluctuation range of ±0.3°C independent from the PAR level. The relative humidity in the climate chambers was 70%. The duration of each experiment was 4 weeks. Leaves of six fully developed plants were harvested using three repetitions per treatment.
Sample Preparation and Extraction of Flavonol Glycosides and Hydroxycinnamic Acid Derivatives
A mixed sample of leaves (without petiole and midrib) taken from six fully developed plants per replication was frozen (−40°C), lyophilized for 72 h and subsequently ground to a fine powder. The samples were either stored in darkness at room temperature less than 6 months until required for analysis or they were immediately used for total RNA extraction. To analyze flavonol glycosides and hydroxycinnamic acid derivatives samples were extracted with 60% aqueous methanol as described (Neugart et al., 2012b).
HPLC-DAD-ESI-MSn Measurement of Flavonol Glycosides and Hydroxycinnamic Acid Derivatives
An HPLC series 1100 by Agilent (Waldbronn, Germany), consisting of a degaser, binary pump, autosampler, column oven and photodiode array detector, was used to quantify the flavonol glycosides and hydroxycinnamic acid derivatives. An ion trap (Agilent series 1100 MSD) with an electrospray ionization (ESI) ion source in negative ionization mode was used as the mass spectrometer. Nitrogen was used as the dry gas (12 L min−1, 350°C) and nebulizer gas (40 psi). Helium was the inert, collision gas in the ion trap. The flavonol glycosides and hydroxycinnamic acid derivatives were separated on a Phenomenex Prodigy column (150 × 3.0 mm, ODS 3, 5 ⌈m, 100 A°) with a C18 security guard column (4 × 3.0 mm, ODS 3, 5 ⌈m, 100 A°) at a temperature of 30°C using a water/acetonitrile gradient (Schmidt et al., 2010a). Solvent A consisted of 99.5% water and 0.5% acetic acid; solvent B was 100% acetonitrile. The following gradient was used for eluent B: 5–7% (0–12 min), 7–9% (12–25 min), 9–12% (25–45 min), 12–15% (45–100 min), 15% isocratic (100–150 min), 15–50% (150–155 min), 50% isocratic (155–165 min), 50–5% (165–170 min) and 5% isocratic (170–175 min). Separation was performed using a flow of 0.4 ml min−1; the measured detector wavelength for quantification was set at 370 nm for non-acylated flavonol glycosides, 330 nm for acylated flavonol glycosides and 320 nm for hydroxycinnamic acid derivatives (Olsen et al., 2009a; Schmidt et al., 2010a). The standards quercetin-3-O-glucoside, kaempferol-3-O-glucoside and chlorogenic acid (Carl Roth GmbH, Karlsruhe, Germany) were used to obtain an external calibration curve ranging from 0.1 to 10 mg 100 ml−1. Mass optimization for the ion optics of the mass spectrometer was performed for quercetin m/z 301. In addition, arbitrary m/z 1000 was used as the target mass in the auto mode to include higher mass fragments. The MSn experiments were performed in auto or manual mode up to MS4 in a scan from m/z 200 to 2000.
The microarray analysis was performed as described (Mewis et al., 2012). Lyophilized and ground leaf material was immediately used for total RNA extraction with the RNeasy Plant Mini Kit including the on-column DNase digestion step (Qiagen GmbH, Hilden, Germany). Microarray analysis was done with 1 mg of total RNA isolated from each of three to four replicates of plants grown at 10°C and 100 μmol m−2 s−1 (LL, low light) or at 10°C and 400 μmol m−2 s−1 (HL, high light), and plants grown at 400 μmol m−2 s−1 and 5°C (LT, low temperature) or at 400 μmol m−2 s−1 and 15°C (HT, high temperature). Agilent One-Color Gene Expression Microarray analysis following the recommendation of MIAME (www.mged.org) was performed at Beckman Coulter Genomics (Morrisville, NC, United States; www.beckmangenomics.com/) using their 2 × 105 K format Brassica Array (Custom Gene Expression Microarray G4503A, 2 × 105 K; (Trick et al., 2009); brassica.bbsrc.ac.uk/). Raw data files are available with accession number E-MTAB-2677 at ArrayExpress (www.ebi.ac.uk/arrayexpress/). The Open Source Microarray Processing Software Robin was used (Lohse et al., 2010; mapman.gabipd.org/web/guest/home) including quality control tools using the default settings. Here the data were shifted to be zero centered and scaled to have unit variance. The PCA was performed by singular value decomposition. The analysis settings were as follows: Normalization with background correction “subtraction” and between array normalization method “scale” with the analysis strategy of linear models (package limma). Results were given as log fold change of expression in relation to different light or temperature treatment. The assignment of the different genes was done by comparison of the translated protein sequences of the 95 K Brassica unigene set with the Arabidopsis TAIR9 database using the newly developed Mercator pipeline for automated sequence annotation (Lohse et al., 2014; mapman.gabipd.org/web/guest/app/mercator). For each identifier the gene with the highest homology was provided with identifier and description. The respective bitscores were classified as follows: very weakly similar (bitscore smaller than 100); weakly similar (bitscore 101–200); moderately similar (bitscore 201–500); highly similar (bitscore greater than 500). Assignment of the different genes represented by identifiers to respective bins and visualization of data sets was realized using MapMan 3.5.1 (Thimm et al., 2004; http://mapman.gabipd.org/web/guest/mapman-download). Complete lists of significant alterations of mRNA abundance induced by either LL, or HL, and either LT, or HT are given as Supplementary Tables 1–4.
Statistical analyses were conducted with the software STATISTICA 9 for Windows (version 9, Statsoft Inc.). Differences due to temperature or radiation were statistically analyzed by a one-factorial ANOVA by a comparison of the factor levels using Tukey's HSD test. The individual values of flavonol gylcosides and hydroxycinnamic acid derivative represent the mean of 3 biological replicates (each biological replicate contained 6 fully developed plants). All tests were conducted at a significance level of 5%.
Results and Discussion
As a distinctive winter crop, kale has to handle low levels of PAR and also low temperature. However, under field conditions there was a certain interaction of PAR and temperature as previously shown for the formation of flavonol glycosides (Neugart et al., 2012a). Therefore, the present study used defined conditions in climate chambers to elucidate the specific effect of either low to moderate PAR levels or typical winter temperature. In order to gain a deeper insight into the adaptation processes when exposed to different PAR levels (100 μmol m−2 s−1 and 400 μmol m−2 s−1 at 10°C) and temperatures (5 and 15°C, at 400 μmol m−2 s−1), the plant's response was evaluated on expression level using microarrays. In addition to the global analysis of transcriptomic alterations, a specific focus of the expression analysis was directed toward genes putatively involved in the phenylpropanoid pathway of B. oleracea.
The comprehensive gene expression analysis was performed using the 2 × 105 K format Brassica Array based on the 95 K Brassica unigene set (Trick et al., 2009) and the microarray processing software Robin (Lohse et al., 2010; mapman.gabipd.org/web/guest/home) using default settings for quality assessment and normalization of data. The principal component analysis of all expression data demonstrates that 54% of expression variation can be explained by the two components PC1 and PC2 (Figure 2). The clear separation of the LL data (100 μmol m−2 s−1, 10°C) and HL data (400 μmol m−2 s−1, 10°C) indicated that light is the factor that mainly determines PC1. On the other hand the rather good separation of the LT data (5°C, 400 μmol m−2 s−1) and the HT data (15°C, 400 μmol m−2 s−1) indicated that temperature is the factor that mainly determines PC2. This was further supported by the fact that the HL data of the other experiments at 10°C and 400 μmol m−2 s−1 positioned between LT and HT data. This data reduction shows that in the two experimental setups with different PAR or different temperature typical expression patterns can be received for the respective situation.
Figure 2. Principal components analysis of expression data. Score plot of the principal components analysis of normalized Brassica Array expression values analyzed with open source microarray processing software Robin. PC1 with standard deviation of 134 explaining a proportion of variance of 0.372, and PC2 with standard deviation of 89, explaining a proportion of variance of 0.164 and a cumulative proportion of 0.536 are shown. PC3, standard deviation 75, proportion of variance 0.115, and PC4, standard deviation 60, proportion of variance 0.075 with a cumulative proportion of 0.726 are not shown.
Differences in Gene Expression in Kale Leaves at Low and Moderate Levels of Photosynthetic Active Radiation
Since light is a very important environmental factor it does not astonish that the expression of almost 18% of the genes represented on the Brassica chip were significantly altered in kale leaves when plants were grown either at 400 μmol m−2 s−1 or 100 μmol m−2 s−1. Gene expression was reliably detectable for 89,621 identifier represented on the 105 K format Brassica Array. Among these, 16,000 reliably detectable elements showed significant differences in expression between growth at 400 μmol m−2 s−1 and 100 μmol m−2 s−1, with 12,666 upregulated at 400 PAR and 3334 upregulated at 100 PAR. The pure numbers of differentially expressed genes in this experiment showed that higher light intensities require almost 4 times more genes to be induced than at low light intensity.
At low light conditions, plants obviously need to absorb sufficient light for photosynthesis, survival, and growth, while at excess light conditions plants have to prevent photo-oxidative damages. Photosynthetic organisms have developed direct and indirect mechanisms for sensing and responding to limiting or excess light (Li et al., 2009; Casal, 2013). For a structured evaluation of differentially expressed genes a functional gene enrichment analysis was performed using the Wilcoxon Rank Sum Test including Benjamini Hochberg correction as integrated in MapMan (Thimm et al., 2004). This analysis revealed that most significant differences in gene expression between plant leaves grown at 400 μmol m−2 s−1 and 100 μmol m−2 s−1 were in the functional categories of protein degradation, transport processes, amino acid metabolism, and in secondary metabolism (Table 1). The huge size of the MapMan bin containing members involved in protein degradation indicates the most important role of E3 ubiquitin ligation in controlling plant processes including vegetative growth, light response and stress tolerance (Mazzucotelli et al., 2006). Up to now only a few of these members have been characterized with a known biological function. Amongst them is COP1, the conserved regulator that has established the ubiquitin-proteasome system as a central subject in light signal transduction (Lau and Deng, 2012). Additionally, the significant enrichment of differentially expressed genes belonging to the categories of sugar and amino acid transport and amino acid metabolism was well expected since tissue and phloem concentrations are not static; these metabolite concentrations are modulated by factors such as light and nutritional and developmental stage (Coruzzi et al., 2015).
Table 1. Significantly altered gene expression in functional categories at 100 PAR-400 PAR light intensity difference.
Light as an environmental factor with strong effects on particular metabolites expects expression differences of genes belonging to the category of secondary metabolism in plants grown at different light intensities. Only the bins with genes of phenylpropanoid and flavonoid metabolism were significantly enriched within this large group of genes coding for enzymes of various secondary pathways (Table 1). This clearly shows that light intensity is a key regulator of gene expression of flavonoid biosynthesis. This further supports the hypothesis that flavonoids and hydroxycinnamic acid derivatives may act as shielding components and scavengers of ROS, thus protecting plants from damage caused by high levels of sunlight (Edreva, 2005).
Temperature Effects on Overall Gene Expression in Kale Leaves
Less genes are differentially expressed in kale leaves grown at 5 or 15°C when compared to the light experiment. Among the reliably detectable 89,621 identifier represented on the 105 K format Brassica Array 3332 elements or 3.7% showed significant differences in expression between growth at 5 and 15°C, with 1525 upregulated at 15°C and 1807 upregulated at 5°C. The pure numbers of differentially expressed genes in this experiment showed that lower temperature require a few more genes to be induced. Functional gene enrichment analysis revealed that most significant differences in gene expression between plant leaves at 5 and 15°C were in the functional categories of phytohormone metabolism, amino acid synthesis, stress and again secondary metabolism (Table 2). Differentially expressed genes were also found in the functional category RNA.regulation of transcription with specific enrichment of certain classes of transcriptional regulators (Table 2).
Table 2. Significantly altered gene expression in functional categories at 5–15°C temperature difference.
Plants from temperate regions show a variable degree of chilling tolerance and can increase their freezing tolerance during exposure to chilling and non-freezing temperatures. This process is known as cold acclimation and involves a complex and interactive relationship between different pathways (Miura and Furumoto, 2013). The significant enrichment of differentially expressed genes belonging to the functional category of jasmonate hormone metabolism supports recent findings that plant hormones are also involved in responses to cold stress (Shi et al., 2015). Additionally, the significant alterations in the category of amino acid synthesis especially aromatic amino acid synthesis point to alterations in phenolic compound production, a prerequisite of the increased lignification and suberin deposition at cold temperatures (Ramakrishna and Ravishankar, 2011). This is further supported by significant enrichment of altered gene expression between plant leaves at 5 and 15°C in the functional category of phenylpropanoid secondary metabolism. Differentially expressed genes were also found in the functional category RNA.regulation of transcription with specific enrichment of transcriptional regulators of the C2H2 and the C3H zink finger families, which are mainly induced at 5°C and bHLH basic helix-loop-helix family members, which are mainly induced at 15°C.
Content of Flavonoid Glycosides and Hydroxycinnamic Acid Derivatives in Kale Leaves Specifically Influenced by Light or Temperature
Kale as a leafy Brassica vegetable contains a high number of structurally different flavonol glycosides and hydroxycinnamic acid derivatives. It is therefore a suitable model plant for investigating the structure-specific formation of these phenolic compounds with respect to temperature and PAR, especially for acylated flavonol glycosides. In the present study defined conditions in climate chambers can clarify the specific effect PAR or temperature can have on the formation of structurally different flavonol glycosides (17 compounds) and hydroxycinnamic acid derivatives (3 compounds) (Table 3).
Table 3. Influence of light or temperature on the content of flavonoid glycosides and hydroxycinnamic acid derivatives in kale plants.
Higher sunlight due to an 8-fold higher reflection caused by a reflective mulch resulted in higher concentrations of flavonoids and anthocyanins in apple fruit (Overbeck et al., 2013). Li and co-worker reported higher concentrations of quercetin glycosides independent of the glycosylated sugar moiety in sun-exposed peel of apple fruits (Li et al., 2013). In kale we found structure-dependent responses of flavonoid glycosides between low and moderate light intensities (Table 3). The main quercetin glycosides quercetin-3-O-sinapoyl-sophoroside-7-O-D-glucoside/quercetin-3-O-sophoroside-7-O-sinapoyl-diglucoside and quercetin-3-O-sophoroside-7-O-D-glucoside were increased with increasing PAR while the main kaempferol glycosides (kaempferol-3-O-sophoroside-7-O-D-glucoside, kaempferol-3-O-feruloyl-sophoroside-7-O-D-glucoside, kaempferol-3-O-hydroxyferuloyl-sophoroside-7-O-D-glucoside, kaempferol-3-O-sinapoyl-sophoroside-7-O-D-glucoside, kaempferol-3-O-hydroxyferuloyl-sophoroside-7-O-diglucose, and kaempferol-3-O-sinapoyl-sophoroside-7-O-diglucoside) were decreased or not affected. The concentrations of other investigated minor flavonol glycosides were also decreased or unchanged at higher PAR levels. It is known that quercetin glycosides have a higher antioxidant activity than their corresponding kaempferol glycosides (Zietz et al., 2010) and are better scavengers of ROS (Edreva, 2005). Out of these feruloyl monoacylated quercetin triglycoside and sinapic acid diacylated quercetin tetraglycoside showed high anti-oxidant activity with 2.54 and 4.19, respectively (Zietz et al., 2010). This implies that even low dosages of PAR result in ROS production and consequences a shift to B-ring polyhydroxylated flavonoids (Agati and Tattini, 2010; Majer et al., 2014). In conclusion, this underlines that ROS concentration presumably at low levels in the present experiment is not damaging the plant's cells. In addition, kaempferol-3-O-caffeoyl-sophoroside-7-O-D-glucoside was increased with higher PAR level whereas caffeoylquinic acid was decreased. This reflects previous investigations where increased UV-B (Neugart et al., 2012b, 2014) lead to the same results. Flavonoids are well-known antioxidants in plants but also act as shielding components (Edreva, 2005). We assume that kaempferol-3-O-caffeoyl-sophoroside-7-O-D-glucoside has a higher shielding activity in the plant than caffeoylquinic acid due to the bigger π-electron system. In contrast, disinapoyl-gentiobiose, which has a high antioxidative activity (4.51), increased at 400 μmol m−2 s−1 while sinapoyl-feruloyl-gentiobiose remained unaffected. In summary, the higher light intensity resulted in structure-dependent responses favoring the main quercetin glycoside quercetin-3-O-sinapoyl-sophoroside-7-O-D-glucoside/quercetin-3-O-sophoroside-7-O-sinapoyl-diglucoside, the caffeic acid monoacylated kaempferol triglycoside kaempferol-3-O-caffeoyl-sophoroside-7-O-D-glucoside and the disinapoyl-gentiobiose.
Lower temperatures resulted in higher concentrations of quercetin glycosides in the investigated kale. Especially the minor compounds quercetin-3-O-hydroxyferuloyl-sophoroside-7-O-D-glucoside and quercetin-3-O-disnapoyl-triglucoside-7-O-D-glucoside that have a high antioxidant activity (2.54 and 4.19, respectively) (Zietz et al., 2010) were increased assuming a higher level of ROS in kale plants at 5°C compared to 15°C. This reflects the shift to B-ring polyhydroxylated glycosides in different plant species due to higher ROS levels in plants caused by lower temperatures (Bilger et al., 2007; Klimov et al., 2008; Moheb et al., 2011; Lavola et al., 2013). Lower temperatures were shown to increase the flavonoids and free phenolic acids in cucumber leaves grown at day/night temperatures of 15/10°C under defined conditions (Chen et al., 2013). Concomitantly, the activity of polyphenol oxidase and different peroxidases was increased at lower day/night temperatures (Chen et al., 2013). However, in kale, different specific flavonol glycosides were increased at 5°C suggesting different ROS to be produced or different plant defense mechanisms to interact with the flavonol glycosides and respective enzymes (catalase, polyphenol oxidase, peroxidase and superoxid dismutase) and vitamins (vitamin C and vitamin E). Additionally, Steindal and co-workers reported higher concentrations of quercetin and kaempferol in broccoli grown at lower day/night temperatures of 15/9°C under defied conditions (Steindal et al., 2013). Interestingly, in kale low temperature of 5°C resulted in higher concentrations of kaempferol tetraglycosides (kaempferol-3-O-hydroxyferuloyl-sophoroside-7-O-diglucoside, kaempferol-3-O-sinapoyl-sophoroside-7-O-diglucoside, kaempferol-3-O-disinapoyl-trigluco-side-7-O-D-glucoside) independent of the acylation pattern. In general, glycosylation causes better water solubility and higher concentrations of sugars act as cryo-protectans in plants at lower temperatures (Theocharis et al., 2012). These kaempferol tetraglycosides can act as antioxidants comparably to their corresponding kaempferol triglycosides that are main flavonol glycosides in kale. The main hydroxycinnamic acid derivatives of kale were not affected or decreased (disinapoyl-gentiobiose) with decreasing temperatures. This was also shown for caffeoylquinic acid in box tree that was not affected by higher altitude including lower temperatures and higher UV-B radiation (Bernal et al., 2013). Even though disinapoyl-gentiobioside which has a high antioxidant potential (4.51) was decreased at lower temperature the total antioxidant activity in the kale plants was increased at lower temperatures (Zietz et al., 2010). Nevertheless, Dong and co-worker found for Eucommia ulmoides a higher correlation of flavonoids and chlorogenic acids with the annual average temperature than with the annual average sunshine duration (Dong et al., 2011). To summarize, in response to low temperatures highly complex sinapic acid acylated flavonol glycosides are favored (quercetin-3-O-disinapoyl-triglucoside-7-O-D-glucoside, kaempferol-3-O-sinapoyl-sophoroside-7-O-diglucoside, kaempferol-3-O-disinapoyl-triglucoside-7-O-D-glucoside) whereas disinapoyl-gentiobiose was decreased at the same time assuming them to be better cryo-protectants.
Identification of the B. oleracea Genes Putatively Involved in the Biosynthesis of Flavonoids and Hydroxycinnamic Acids
Thanks to the identification of genes involved in the phenylpropanoid biosynthetic pathway in several plant species (Milkowski and Strack, 2010; Vogt, 2010; Fraser and Chapple, 2011) the B. oleracea genome could be screened for respective homologs. All relevant genes (Figure 1) were selected by sequence comparison of predicted genes from the B. oleracea genome using the respective coding sequences of A. thaliana and the genome annotation system in e!EnsemblGenomes (Kersey et al., 2014). Genes coding for flavone synthase (FNSII) of Glycine max and flavonoid 3′5′-hydroxylase (F3′5′H) of Vitis vinifera were also included, however, putative homologs to FNSII and F3′5′H were not found in the B. oleracea genome. Predicted B. oleracea genes putatively involved in the biosynthesis of flavonol aglycons (gray and blue pathways in Figure 1) are presented in Table 4, and genes putatively involved in the biosynthesis of phenylpropanoids (orange pathway in Figure 1) are shown in Table 5. Because during evolution the Brassica genomes experienced whole genome triplication followed by massive gene loss and frequent reshuffling of triplicated genomic blocks (Liu et al., 2014), we expected to find in most cases more than one gene homolog in B. oleracea when compared to A. thaliana. In case of flavonoid 3′-hydroxylase F3′H and sinapaldehyde dehydrogenase REF1 only single genes could be found in the B. oleracea genome, while for cinnamate 4-hydroxylase C4H and isoform 5 of 4-coumarate-CoA ligase 4CL5 five homologous genes are present in the B. oleracea genome.
Table 4. Expression difference of genes putatively involved in the biosynthesis of flavonol aglycons in Brassica oleracea.
Table 5. Expression difference of genes putatively involved in the phenylpropanoid biosynthesis in Brassica oleracea.
The Brassica community microarray resource was built with the 95K Brassica unigene set assembled with sequences originating from B. napus, B. rapa, and B. oleracea (Trick et al., 2009). Therefore, all identified B. oleracea genes putatively involved in flavonoid biosynthesis were analyzed for their most specific 60-mer oligonucleotides present on the 2 × 105 K format Brassica Array. In some cases, e.g. for the three CHS genes, gene specific microarray elements could be obtained for all individual isogenes (Table 4). While in other cases, only a single microarray element was present that is highly identical to all selected homologs thus hampering differentiation, e.g. between the four different genes of F3H (Table 5). In addition, there are individual genes in B. oleracea where there is no highly homologous, specific 60-mer oligonucleotide present on this most comprehensive 105 K Brassica Array.
Specific Expression Changes of Genes Putatively Involved in the Biosynthesis of Flavonoids Induced by Light or Temperature
As already shown, different light intensities lead to significant enrichment of differentially expressed genes of secondary metabolism (Table 1), especially predicted genes putatively involved in the biosynthesis of flavonoids and phenylpropanoids. As evident from Table 4 most of the predicted genes putatively involved in the core biosynthesis of flavonol aglycons in kale were strongly increased with increasing light. It is known from various studies that light is one of the most important environmental factors affecting flavonoid biosynthesis, involving R2R3 MYB transcription factors and a COP1 mediated signaling pathways (Zoratti et al., 2014). In fact, a closer analysis of array elements with high homology to MYB transcription factors revealed light-increased expression of CX267652 (log2-difference LL-HL -2.30), specific for the predicted two B. oleracea MYB1 genes, Bo1g144680 and Bo7g117880. This underpins a coordinated light-increased expression of at least one gene coding for each step of flavonoid biosynthesis in kale, resulting in accumulation of higher levels of flavonoids (Table 3). Moreover, the higher light intensity caused a structure-dependent response favoring the main quercetin glycosides quercetin-3-O-sinapoyl-sophoroside-7-O-D-glucoside/ quercetin-3-O-sophoroside-7-O-sinapoyl-diglucoside. This preferential formation of quercetin as the flavonoid aglycone at higher light intensity is further supported by the very strong induction of Bo9g174880, coding for the flavonoid 3′-hydroxylase. This gene is among the most significantly, differentially expressed genes in this analysis. It further underlines the light dependent increased synthesis of shielding components and scavengers protecting plants from damage caused by high levels of sunlight (Edreva, 2005).
As shown previously, different temperature lead to significant enrichment of differentially expressed genes of secondary metabolism (Table 2), especially predicted genes putatively involved in the biosynthesis of phenylpropanoids. As evident from Table 4, significant temperature specific responses were only found in predicted genes putatively involved in the first three steps of core biosynthesis, phenylalanine ammonia lyase (PAL), cinnamate 4-hydroxylase (C4H) and 4-coumarate-CoA ligase (4CL) (Table 4). These first three steps are in common in flavonoid and phenylpropanoid biosynthesis. Thus, temperature effects on gene expression are also exerted on the level of genes involved in hydroxycinnamic acid synthesis, what is discussed in more detail in the next chapter.
Specific Expression Changes of Genes Putatively Involved in the Biosynthesis of Hydroxycinnamic Acids Induced by Light or Temperature
At higher light intensities, the main quercetin glycosides and kaempferol-3-O-caffeoyl-sophoroside-7-O-D-glucoside were significantly accumulating (Table 3), together with disinapoyl-gentiobiose, which has a high antioxidant potential (4.51). This elevated formation of hydroxycinnamic acid derivatives is accompanied by an increased expression of genes presumably involved in the synthesis of hydroxycinnamic acid derivatives at higher PAR levels (Table 5). Up to now, only a few genes involved in these biosynthetic pathways are known from research in Arabidopsis (Fraser and Chapple, 2011). For example, genes directly involved in disinapoyl-gentiobiose biosynthesis, still needs to be discovered, as well as isoforms catalyzing the acylation of flavonoids. Gene identification becomes so difficult, because the model organism Arabidopsis is lacking disinapoyl-gentiobiose, and also acylation of flavonols (Saito et al., 2013). In addition, candidate genes involved in the biosynthetic steps all belong to huge gene families that need further detailed exploration (Sasaki et al., 2014).
As shown previously, different temperature lead to significant enrichment of differentially expressed genes putatively involved in the biosynthesis of phenylpropanoids (Table 2). As evident from Tables 4, 5, significant temperature specific responses were found in predicted genes putatively involved in the first three steps of biosynthesis (PAL, C4H, and 4CL) (Table 4), and hydroxycinnamoyl-transferase (HCT) (Table 5).
Within the identified genes there is evidence of an isoform specific response to different stimuli: while both phenylalanine ammonia lyase1 genes (PAL1, Bo4g03091, Bo4g186590) are strongly induced with higher light and unaffected by temperature changes, one or both of the PAL2 genes (Bo6g067250, Bo8g082620) were unaffected by light and stronger expressed at low temperature. A functional specialization has been shown previously for PAL1 and PAL2 in Arabidopsis for an abiotic environmental-triggered flavonoid biosynthesis (Olsen et al., 2008, 2009b). The isoform specific expression changes to different stimuli proved functional specialization of PAL also in kale. Furthermore, a specific metabolic function of individual isoforms has recently also been discovered for 4-coumarate-CoA ligase 4CL (Li et al., 2015). The four isoforms of 4CL in Arabidopsis have overlapping yet distinct roles in phenylpropanoid metabolism, with 4CL1 and 4CL2 predominant in lignifying cells (Ehlting et al., 1999; Saito et al., 2013), and 4CL3 with a distinct role in flavonoid metabolism (Li et al., 2015). Comparable roles might also be expected in kale, as 4CL1 and 4CL2 isoform expression is unaffected by different light and temperature, while expression of the single 4CL3 isoform is significantly upregulated with high light, thus co-expressing with the other genes of core biosynthesis of flavonol aglycons. Two isoforms of 4CL5 were significantly increased at high temperature and light. It might be speculated that this is related to increased levels of disinapoyl-gentiobiose at high temperature and light. It was demonstrated previously in Arabidopsis, that a specific 4CL gene product is able to activate sinapic acid (Hamberger and Hahlbrock, 2004). Together with specific isoforms of cinnamate 4-hydroxylase C4H and hydroxycinnamoyl-transferase HCT co-expression of the hydroxycinnamic acid pathway at high temperature and high light was detectable.
Expression Pattern of Gene Family Members Putatively Involved in Flavonoid Modification
In kale flavonoid glycosylation occurs preferentially at the C-3 and the C-7 position of the aglycons, with only glucose residues being detected as conjugates (Table 3). So far no arabinose or rhamnose was found to be conjugated with flavonols in kale leaves. In addition, diversity is increased by a variety of acylations (Schmidt et al., 2010b), a modification not present in A. thaliana. Another differentiation to Arabidopsis is the accumulation of other hydroxycinnamic acid derivatives: while disinapoyl-gentiobiose and sinapoyl-feruloyl-gentiobiose are highly abundant in kale, the major hydroxycinnamic acid derivative in Arabidopsis, the sinapoylmalate is absent in kale.
Huge gene families of glycosyltransferases (Bowles et al., 2005) and acyltransferases (Sasaki et al., 2014) are involved in the synthesis of chemically diverse flavonoids and hydroxycinnamic acid derivatives. Dependent on their subcellular location within the plant cell glycosylation is either performed by UDP-sugar dependent glycosyltransferases (UGT) in the cytosol, or by acyl-glucose dependent glucosyltransferases in the vacuole (glycoside hydrolase family 1 type glucosyltransferases, GH1-GT). Also dependent on localization acylation is either performed by acyl-CoA dependent acyltransferases in the cytosol (BAHD type family of acyltransferases, BAHD-AT), or by acyl-glucose dependent acyltransferases (serine carboxypeptidase-like family of acyltransferases, SCPL-AT) in the vacuole. In the Arabidopsis genome 116 UGT genes are present while the B. oleracea genome harbors at least 189 genes coding for predicted UDP-sugar glycosyltransferases. The gene family of vacuolar glycoside hydrolase family 1 type glucosyltransferases is a little smaller with 48 members in Arabidopsis and 75 GH1-GT present in the B. oleracea genome. While for the BAHD type family of acyltransferases in Arabidopsis 64 genes are known, at least 99 genes encoding putative BAHD-AT can be extracted from the B. oleracea genome. Finally the 54 genes coding for serine carboxypeptidase-like acyltransferases in Arabidopsis are exceeded by the 64 SCPL-AT present in the B. oleracea genome. Up to now a couple of these glycosyltransferases and acyltransferases have been determined with their functional preferences, and very few of these can be linked to flavonoid and anthocyanin modifications or hydroxycinnamic acid derivative biosynthesis (Fraser and Chapple, 2011; Saito et al., 2013).
Since these gene families are so large and most of their members still lacking functional characterization a differentiation or distinct assignment to certain pathways or enzymatic steps will not be possible at that point. However, extraction of some isoforms with correlated expression patterns might be possible and give hints to some promising candidates. For example, candidates of the UGT73B and UGT73C family with predicted functions as flavonol 3-O-glucosyltransferase and flavonol 7-O-glucosyltransferase (Table 6) are among the most induced UGT family members at high light intensities. Furthermore, the two microarray elements EV108083 and JCVI_27911 are among the top 100 most cold responsive elements. These microarray elements correspond two adjacent genes on chromosome 4 predicted as UGT family members of UGT73C5/6. Whether one or both of these genes are responsible for the enhanced accumulation of kaempferol tetraglycosides catalyzing the transfer of the second glucose residue to the flavonol 7-O position in kale at low temperature will be the focus of further studies.
Table 6. Expression difference of genes putatively involved in glucosylation and acylation of flavonoids in Brassica oleracea.
In general, additional experimental evidence will be needed to unequivocally show the functional involvement of specific genes and gene products in flavonoid related pathways in kale. This can be done only with some degree of success in in vitro experiments using recombinant gene products analyzed for their enzyme activity with only limiting substrate availability. An in vivo functional analysis will be complicated by the fact that knock-out mutants of kale are not available and gene editing techniques in Brassicaceae are still in its infancies (Lawrenson et al., 2015). Furthermore, related knock-out mutants in Arabidopsis are not very useful, because specific flavonoid modifications present in kale are absent in Arabidopsis. However, this latter point might also be adventitious, since it allows functional in vivo characterization by overexpressing the respective B. oleracea gene candidates in transgenic Arabidopsis plants. A subsequent analysis of the flavonoids in these plants may result in more kale like flavonoid profiles and may give valuable hints on the enzyme specificities.
Kale as a typical winter crop is exposed to less favorable growth conditions. Having to deal with the disadvantageous light and/or temperature conditions kale leaves accumulate high concentrations of a large number of naturally occurring flavonol glycosides and hydroxycinnamic acid derivatives. This makes kale an increasingly important vegetable, since its secondary metabolites are considered as health promoting substances. To understand the underlying processes and expression networks, the temperature and light dependent alterations of gene expression were analyzed in B. oleracea var. sabellica using the most comprehensive Brassica microarray. Expression differences were correlated to the structure-dependent response of flavonoid glycosides and hydroxycinnamic acid derivatives to alterations in either temperature or PAR. Altered metabolite accumulation was mainly reflected on gene expression level of core biosynthetic pathway genes and gave hints to an isoform specific functional specialization. However, additional experimental evidence will be needed to unequivocally show the functional involvement of specific genes and gene products in flavonoid modification pathways. On the one hand such an analysis will be complicated by the fact that knock-out mutants are not available and some flavonoid modifications are specific in kale and absent in Arabidopsis. On the other hand this will be adventitious since it allows functional in vivo characterization by overexpressing the respective B. oleracea gene candidates in transgenic Arabidopsis plants with subsequent flavonoid profiling.
SN and AK designed the study, carried out the metabolite analysis, helped in data interpretation, and made a draft of parts of the manuscript. RZ carried out the expression analysis, data interpretation and wrote the manuscript. All authors have given final approval for this version to be published.
We sincerely thank the DFG (Deutsche Forschungsgemeinschaft) for funding (Projects KR-2066/3-2). The Leibniz Association is gratefully acknowledged for supporting this work by an additional grant.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank Andrea Maikath (Leibniz Institute of Vegetable and Ornamental Crops) for extracting the RNA for microarray studies and the excellent technical support, and Ursula Zentner, Elke Büsch and Uwe Kunert and his team at the Leibniz Institute of Vegetable and Ornamental Crops for cultivating the kale plants. We also thank Dr. Marc Lohse (Max-Planck-Institute of Molecular Plant Physiology, Potsdam, Germany) for helpful advices regarding the analysis of microarray data.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2016.00326
Bernal, M., Llorens, L., Julkunen-Tiitto, R., Badosa, J., and Verdaguer, D. (2013). Altitudinal and seasonal changes of phenolic compounds in Buxus sempervirens leaves and cuticles. Plant Physiol. Biochem. 70, 471–482. doi: 10.1016/j.plaphy.2013.06.012
Bilger, W., Rolland, M., and Nybakken, L. (2007). UV screening in higher plants induced by low temperature in the absence of UV-B radiation. Photochem. Photobiol. Sci. 6, 190–195. doi: 10.1039/b609820g
Calderón-Montaño, J. M., Burgos-Morón, E., Pérez-Guerrero, C., and López-Lázaro, M. (2011). A review on the dietary flavonoid kaempferol. Mini Rev. Med. Chem. 11, 298–344. doi: 10.2174/138955711795305335
Chen, S., Jin, W., Liu, A., Zhang, S., Liu, D., Wang, F., et al. (2013). Arbuscular mycorrhizal fungi (AMF) increase growth and secondary metabolism in cucumber subjected to low temperature stress. Sci. Hortic. 160, 222–229. doi: 10.1016/j.scienta.2013.05.039
Coruzzi, G., Last, R., Dudareva, N., and Amrhein, N. (2015). “Amino Acids,” in Biochemistry and Molecular Biology of Plants, ed W. G. Bob, B. Buchanan, Russell, and L. Jones 2nd Edn. (Blackwell: John Wiley & Sons. Ltd.), 289–376.
Dong, J., Ma, X., Wei, Q., Peng, S., and Zhang, S. (2011). Effects of growing location on the contents of secondary metabolites in the leaves of four selected superior clones of Eucommia ulmoides. Ind. Crops Prod. 34, 1607–1614. doi: 10.1016/j.indcrop.2011.06.007
Ehlting, J., Büttner, D., Wang, Q., Douglas, C. J., Somssich, I. E., and Kombrink, E. (1999). Three 4-coumarate: coenzyme a ligases in Arabidopsis thaliana represent two evolutionarily divergent classes in angiosperms. Plant J. 19, 9–20. doi: 10.1046/j.1365-313X.1999.00491.x
Fallovo, C., Schreiner, M., Schwarz, D., Colla, G., and Krumbein, A. (2011). Phytochemical changes induced by different nitrogen supply forms and radiation levels in two leafy Brassica species. J. Agric. Food Chem. 59, 4198–4207. doi: 10.1021/jf1048904
Ferreres, F., Fernandes, F., Oliveira, J. M. A., Valentão, P., Pereira, J. A., and Andrade, P. B. (2009). Metabolic profiling and biological capacity of Pieris brassicae fed with kale (Brassica oleracea L. var. acephala). Food Chem. Toxicol. 47, 1209–1220. doi: 10.1016/j.fct.2009.02.014
Fiol, M., Adermann, S., Neugart, S., Rohn, S., Mügge, C., Schreiner, M., et al. (2012). Highly glycosylated and acylated flavonols isolated from kale (Brassica oleracea var. sabellica) — structure–antioxidant activity relationship. Food Res. Int. 47, 80–89. doi: 10.1016/j.foodres.2012.01.014
Hamberger, B., and Hahlbrock, K. (2004). The 4-coumarate: CoA ligase gene family in Arabidopsis thaliana comprises one rare, sinapate-activating and three commonly occurring isoenzymes. Proc. Natl. Acad. Sci. U.S.A. 101, 2209–2214. doi: 10.1073/pnas.0307307101
Kersey, P. J., Allen, J. E., Christensen, M., Davis, P., Falin, L. J., Grabmueller, C., et al. (2014). Ensembl Genomes 2013: scaling up access to genome-wide data. Nucleic Acids Res. 42, D546–D552. doi: 10.1093/nar/gkt979
Klimov, S. V., Burakhanova, E. A., Dubinina, I. M., Alieva, G. P., Sal'nikova, E. B., Olenichenko, N. A., et al. (2008). Suppression of the source activity affects carbon distribution and frost hardiness of vegetating winter wheat plants. Russ. J. Plant Physiol. 55, 308–314. doi: 10.1134/S1021443708030035
Lavola, A., Nybakken, L., Rousi, M., Pusenius, J., Petrelius, M., Kellomaki, S., et al. (2013). Combination treatment of elevated UVB radiation, CO2 and temperature has little effect on silver birch (Betula pendula) growth and phytochemistry. Physiol. Plant. 149, 499–514. doi: 10.1111/ppl.12051
Lawrenson, T., Shorinola, O., Stacey, N., Li, C., Ostergaard, L., Patron, N., et al. (2015). Induction of targeted, heritable mutations in barley and Brassica oleracea using RNA-guided Cas9 nuclease. Genome Biol. 16:258. doi: 10.1186/s13059-015-0826-7
Lepiniec, L., Debeaujon, I., Routaboul, J. M., Baudry, A., Pourcel, L., Nesi, N., et al. (2006). Genetics and biochemistry of seed flavonoids. Ann. Rev. Plant Biol. 57, 405–430. doi: 10.1146/annurev.arplant.57.032905.105252
Li, Y., Kim, J. I., Pysh, L., and Chapple, C. (2015). Four isoforms of Arabidopsis thaliana 4-coumarate: CoA ligase (4CL) have overlapping yet distinct roles in phenylpropanoid metabolism. Plant Physiol. 169, 2409–2421. doi: 10.1104/pp.15.00838
Liu, S., Liu, Y., Yang, X., Tong, C., Edwards, D., Parkin, I. A., et al. (2014). The Brassica oleracea genome reveals the asymmetrical evolution of polyploid genomes. Nat. Commun. 5:3930. doi: 10.1038/ncomms4930
Lohse, M., Nagel, A., Herter, T., May, P., Schroda, M., Zrenner, R., et al. (2014). Mercator: a fast and simple web server for genome scale functional annotation of plant sequence data. Plant Cell Environ. 37, 1250–1258. doi: 10.1111/pce.12231
Lohse, M., Nunes-Nesi, A., Krüger, P., Nagel, A., Hannemann, J., Giorgi, F. M., et al. (2010). Robin: an intuitive wizard application for R-based expression microarray quality assessment and analysis. Plant Physiol. 153, 642–651. doi: 10.1104/pp.109.152553
Maalik, A., Khan, F. A., Mumtaz, A., Mehmood, A., Azhar, S., Atif, M., et al. (2014). Pharmacological applications of quercetin and its derivatives: a short review. Trop. J. Pharmaceut. Res. 13, 1561–1566. doi: 10.4314/tjpr.v13i9.26
Majer, P., Neugart, S., Krumbein, A., Schreiner, M., and Hideg, É. (2014). Singlet oxygen scavenging by leaf flavonoids contributes to sunlight acclimation in Tilia platyphyllos. Environ. Exp. Bot. 100, 1–9. doi: 10.1016/j.envexpbot.2013.12.001
Mazzucotelli, E., Belloni, S., Marone, D., De Leonardis, A. M., Guerra, D., Di Fonzo, N., et al. (2006). The E3 ubiquitin ligase gene family in plants: regulation by degradation. Curr. Genomics 7, 509–522. doi: 10.2174/138920206779315728
Mewis, I., Schreiner, M., Nguyen, C. N., Krumbein, A., Ulrichs, C., Lohse, M., et al. (2012). UV-B irradiation changes specifically the secondary metabolite profile in broccoli sprouts: induced signaling overlaps with defense response to biotic stressors. Plant Cell Physiol. 53, 1546–1560. doi: 10.1093/pcp/pcs096
Neugart, S., Fiol, M., Schreiner, M., Rohn, S., Zrenner, R., Kroh, L. W., et al. (2013). Low and moderate photosynthetically active radiation affects the flavonol glycosides and hydroxycinnamic acid derivatives in kale (Brassica oleracea var. sabellica) dependent on two low temperatures. Plant Physiol. Biochem. 72, 161–168. doi: 10.1016/j.plaphy.2013.04.002
Neugart, S., Fiol, M., Schreiner, M., Rohn, S., Zrenner, R., Kroh, L. W., et al. (2014). Interaction of moderate UV-B exposure and temperature on the formation of structurally different flavonol glycosides and hydroxycinnamic acid derivatives in Kale (Brassica oleracea var. sabellica). J. Agric. Food Chem. 62, 4054–4062. doi: 10.1021/jf4054066
Neugart, S., Kläring, H.-P., Zietz, M., Schreiner, M., Rohn, S., Kroh, L. W., et al. (2012a). The effect of temperature and radiation on flavonol aglycones and flavonol glycosides of kale (Brassica oleracea var. sabellica). Food Chem. 133, 1456–1465. doi: 10.1016/j.foodchem.2012.02.034
Neugart, S., Zietz, M., Schreiner, M., Rohn, S., Kroh, L. W., and Krumbein, A. (2012b). Structurally different flavonol glycosides and hydroxycinnamic acid derivatives respond differently to moderate UV-B radiation exposure. Physiol. Plant. 145, 582–593. doi: 10.1111/j.1399-3054.2012.01567.x
Nybakken, L., Hörkkä, R., and Julkunen-Tiitto, R. (2012). Combined enhancements of temperature and UVB influence growth and phenolics in clones of the sexually dimorphic Salix myrsinifolia. Physiol. Plant. 145, 551–564. doi: 10.1111/j.1399-3054.2011.01565.x
Olsen, H., Aaby, K., and Borge, G. I. A. (2009a). Characterization and quantification of flavonoids and hydroxycinnamic acids in curly kale (Brassica oleracea L. Convar. acephala Var. sabellica) by HPLC-DAD-ESI-MSn. J. Agric. Food Chem. 57, 2816–2825. doi: 10.1021/jf803693t
Olsen, K. M., Lea, U. S., Slimestad, R., Verheul, M., and Lillo, C. (2008). Differential expression of four Arabidopsis PAL genes; PAL1 and PAL2 have functional specialization in abiotic environmental-triggered flavonoid synthesis. J. Plant Physiol. 165, 1491–1499. doi: 10.1016/j.jplph.2007.11.005
Olsen, K. M., Slimestad, R., Lea, U. S., Brede, C., Løvdal, T., Ruoff, P., et al. (2009b). Temperature and nitrogen effects on regulators and products of the flavonoid pathway: experimental and kinetic model studies. Plant Cell Environ. 32, 286–299. doi: 10.1111/j.1365-3040.2008.01920.x
Saito, K., Yonekura-Sakakibara, K., Nakabayashi, R., Higashi, Y., Yamazaki, M., Tohge, T., et al. (2013). The flavonoid biosynthetic pathway in Arabidopsis: structural and genetic diversity. Plant Physiol. Biochem. 72, 21–34. doi: 10.1016/j.plaphy.2013.02.001
Schmidt, S., Zietz, M., Schreiner, M., Rohn, S., Kroh, L. W., and Krumbein, A. (2010a). Genotypic and climatic influences on the concentration and composition of flavonoids in kale (Brassica oleracea var. sabellica). Food Chem. 119, 1293–1299. doi: 10.1016/j.foodchem.2009.09.004
Schmidt, S., Zietz, M., Schreiner, M., Rohn, S., Kroh, L. W., and Krumbein, A. (2010b). Identification of complex, naturally occurring flavonoid glycosides in kale (Brassica oleracea var. sabellica) by high-performance liquid chromatography diode-array detection/electrospray ionization multi-stage mass spectrometry. Rapid Commun. Mass Spectrom. 24, 2009–2022. doi: 10.1002/rcm.4605
Steindal, A. L. H., Mølmann, J., Bengtsson, G. B., and Johansen, T. J. (2013). Influence of day length and temperature on the content of health-related compounds in broccoli (Brassica oleracea L. var. italica). J. Agric. Food Chem. 61, 10779–10786. doi: 10.1021/jf403466r
Thimm, O., Bläsing, O., Gibon, Y., Nagel, A., Meyer, S., Krüger, P., et al. (2004). MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 37, 914–939. doi: 10.1111/j.1365-313X.2004.02016.x
Trick, M., Cheung, F., Drou, N., Fraser, F., Lobenhofer, E. K., Hurban, P., et al. (2009). A newly-developed community microarray resource for transcriptome profiling in Brassica species enables the confirmation of Brassica-specific expressed sequences. BMC Plant Biol. 9:50. doi: 10.1186/1471-2229-9-50
Verkerk, R., Schreiner, M., Krumbein, A., Ciska, E., Holst, B., Rowland, I., et al. (2009). Glucosinolates in Brassica vegetables: the influence of the food supply chain on intake, bioavailability and human health. Mol. Nutr. Food Res. 53(Suppl. 2), S219. doi: 10.1002/mnfr.200800065
Yamaguchi-Shinozaki, K., and Shinozaki, K. (2006). Transcriptional regulatory networks in cellular responses and tolerance to dehydration and cold stresses. Ann. Rev. Plant Biol. 57, 781–803. doi: 10.1146/annurev.arplant.57.032905.105444
Zietz, M., Weckmüller, A., Schmidt, S., Rohn, S., Schreiner, M., Krumbein, A., et al. (2010). Genotypic and climatic influence on the antioxidant activity of flavonoids in kale (Brassica oleracea var. sabellica). J. Agric. Food Chem. 58, 2123–2130. doi: 10.1021/jf9033909
Keywords: Brassica, flavonol glycosides, temperature, photosynthetically active radiation (PAR), gene expression, microarray
Citation: Neugart S, Krumbein A and Zrenner R (2016) Influence of Light and Temperature on Gene Expression Leading to Accumulation of Specific Flavonol Glycosides and Hydroxycinnamic Acid Derivatives in Kale (Brassica oleracea var. sabellica). Front. Plant Sci. 7:326. doi: 10.3389/fpls.2016.00326
Received: 10 December 2015; Accepted: 03 March 2016;
Published: 30 March 2016.
Edited by:Antonio Encina García, University of León, Spain
Reviewed by:Dayakar Badri, Hill's Pet Nutrition Inc., USA
Biswapriya Biswavas Misra, University of Florida, USA
Copyright © 2016 Neugart, Krumbein and Zrenner. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Rita Zrenner, firstname.lastname@example.org