Original Research ARTICLE
Organ-Specific Quantitative Genetics and Candidate Genes of Phenylpropanoid Metabolism in Brassica oleracea
- 1Group of Genetics, Breeding and Biochemistry of Brassicas, Misión Biológica de Galicia - Consejo Superior de Investigaciones Científicas (MBG-CSIC), Pontevedra, Spain
- 2Department of Horticulture, Faculty of Agriculture, Ain Shams University, Cairo, Egypt
- 3Research Group on Quality, Safety and Bioactivity of Plant Foods, Department of Food Science and Technology, Centro de Edafología y Biología Aplicada del Segura - Consejo Superior de Investigaciones Científicas (CEBAS-CSIC), Murcia, Spain
Phenolic compounds are proving to be increasingly important for human health and in crop development, defense and adaptation. In spite of the economical importance of Brassica crops in agriculture, the mechanisms involved in the biosynthesis of phenolic compounds presents in these species remain unknown. The genetic and metabolic basis of phenolics accumulation was dissected through analysis of total phenolics concentration and its individual components in leaves, flower buds, and seeds of a double haploid (DH) mapping population of Brassica oleracea. The quantitative trait loci (QTL) that had an effect on phenolics concentration in each organ were integrated, resulting in 33 consensus QTLs controlling phenolics traits. Most of the studied compounds had organ-specific genomic regulation. Moreover, this information allowed us to propose candidate genes and to predict the function of genes underlying the QTL. A number of previously unknown potential regulatory regions involved in phenylpropanoid metabolism were identified and this study illustrates how plant ontogeny can affect a biochemical pathway.
A wide range of secondary metabolites that are synthesized by plants are not required in the primary processes of growth and development but are of vital importance for plant interaction with the environment, for their defense mechanism and for their reproductive strategy (Cheynier et al., 2013). Phenolic compounds are the most widely distributed secondary metabolites, ubiquitously present in the plant kingdom. They are synthesized from phenylalaline via the shikimate/phenylpropanoid pathway (Figure 1). Research in Arabidopsis has shown that the pathway starts with the conversion of phenylalanine into cinnamic acid by phenylalanine ammonia lyase (PAL). The cinnamate is hydroxylated by cinnamate 4-hydroxylase (C4H) to form p-coumaric acid. The 4-coumarate:CoA ligase enzyme (4CL) then converts p-coumaric acid into p-coumaroyl CoA by addiction of a CoA thioester, which is the precursor of variuos phenylpropanoid derivatives, including flavonoids, lignins, and isoflavonoids (Fraser and Chapple, 2011). More than 8000 phenolic compounds have been isolated and reported from plant sources which have diverse structural configurations and polarities (Robbins, 2003).
Figure 1. Proposed phenylpropanoid pathway in the Brassicaceae family. ALDH, aldehyde dehydrogenase; CAD, cinnamyl alcohol dehydrogenase; CCR, cinnamyl CoA reductase; CCoAMT, caffeoyl CoA O-methyltransferase; C4H, cinnamate-4-hydroxylase; CHS, chalcone synthase; CHI, chalcone isomerase; COMT, caffeic acid O-methyltransferase; C3H, coumarate 3-hydroxylase; 4CL, 4-coumaroyl CoA-ligase; F3GT, flavonol 3-O-glucosyltransferase; F7GT, flavonol 7-O-glucosyltransferase; FLS, flavonol synthase; F3H, flavanone 3-hydroxylase; F3'H, flavonoid-3′ hydroxylase; F5H, ferulate 5-hydroxylase; HCT, hydroxycinnamoyl-CoA:shikimate hydroxycinnamoyl transferase; OMT, O-methyltransferase; PAL, phenylalanine ammonia-lyase; glc, glucose.
Phenolic compounds in plants are involved in a variety of functions. They can modulate essential physiological processes, such as transcriptional regulation but also they are involved in processes such us growth, development, adaptation, symbiosis, diseases, and responses to pathogen attack (Bhattacharya et al., 2010; Mandal et al., 2010). In humans, phenolic compounds have received considerable attention for being potentially protective factors against cancer and heart diseases (Crozier et al., 2009), in part because of their potent antioxidative properties and their ubiquity in a wide range of commonly consumed foods of plant origin.
The Brassicaceae family includes horticultural species consumed by people all over the world considered important food crops in China, Japan, India, and European countries. The main vegetable species of this family is Brassica oleracea, which includes many well-known crops such as kale, cabbage, broccoli, Brussels sprouts, and cauliflower. These vegetables are known for their characteristic and high content of secondary metabolites such as glucosinolates (GLS) and phenolic compounds widely studied for its beneficial properties (Cartea et al., 2010; Dinkova-Kostova and Kostov, 2012). For a better understanding of metabolic pathways and their regulations it is important to mention the significance of interactions and possible competition for precursors between GLS pathway and the phenylpropanoid pathways, with respect to environmental effects and genetics (Francisco et al., 2012; Kim et al., 2015).
The most common phenolic compounds presents in Brassica are flavonoids derivatives (flavonols and anothocyannins) and hydroxycinnamic acids. The flavonoids consist of 15 carbons with two aromatic rings connected by a three-carbon bridge (C6-C3-C6 carbon skeleton). Among the flavonoids, flavonols of quercetin, kaempferol, and isorhamnetin glycosilated and/or acylated by different hydroxycinnamic acids are the most common in Brassica crops. Hydroxycinnamic acids, on the other hand, are aromatic compounds with a three-carbon side chain (C6–C3). Caffeic, ferulic, p-coumaric, and sinapic acids are appearing in high quantities in Brassica, often found in conjugation with sugar moieties or other hydroxycinnamics (Cartea et al., 2010).
In a green part of Brassica vegetables it can be found up to 40 structurally different flavonoids and hydroxycinnamic acids. Moreover, the phenolic profile and concentration may vary between different tissues of the plant and within different populations of the same plant species (Ferreres et al., 2006; Ayaz et al., 2008; Sousa et al., 2008; Francisco et al., 2009). Several studies have found that the regulation of phenol production occurs principally through changes in the transcription rate of the biosynthetic genes (Takos et al., 2006; Pereira et al., 2009; Carbonell-Bejerano et al., 2014; Koyama et al., 2014). During the last years, there is an increasing interest in phenolics regulation. Extensive studies on model organisms like tobacco, Arabidopsis thaliana, and Medicago truncatula facilitate the understanding of regulation system of the phenolic pathway (Tamagnone et al., 1998; Pang et al., 2008; Zhao and Dixon, 2009; Verdier et al., 2012). However, in other plant species, the genetic control of phenolics production has been the subject of very few genetic studies. In Brassica vegetables the majority of them are focused in the molecular regulation of sinapine metabolism in oilseed rape (B. napus L. var. napus) (Milkowski et al., 2004; Baumert et al., 2005). Seeds of oilseed rape accumulate high amounts of sinapate esters. Given the antinutritive properties that these metabolites confer on the protein fraction, low sinapate ester content is one major aim of breeding programs designed to increase the nutritional value of seeds from B. napus. Based on the absence of soluble sinapate esters in the A. thaliana mutant sin1, the homologous gene from B. napus was identified and the sequence information was used for antisense suppression strategy to obtain lines with reduced levels of sinapate esters (Nair et al., 2000). Comprehensive analysis of transgenic seeds with other silenced genes involved in sinapate ester biosynthesis revealed that both the amount of sinapoylglucose and that of the minor sinapate esters can be drastically reduced to trace amounts (Hüsken et al., 2005). In a more recent study, Rezaeizad et al. (2011) detected SSR markers linked to phenolic compounds also in oilseed rape using an association mapping approach, which could be employed in a marker assisted selection.
In spite of the economical importance of Brassica species in agriculture, the mechanisms involved in the biosynthesis of defense and/or health-related phenolic compounds presents in these plants remain unknown. This study aims to identify genomic regions controlling phenolic composition and content in three different organs (leaves, flower buds, and seeds) in a double haploid (DH) population of B. oleracea. We utilized a quantitative trait locus (QTL) mapping approach to better understand the genetic basis of phenolic accumulation within plant ontogeny. For major QTL regions, candidate genes were proposed using two different approaches, by in silico search in the B. oleracea available sequences and by genome comparison with A. thaliana through analysis of syntenic regions.
Materials and Methods
A double haploid (DH) mapping population (BolTBDH) was employed in this work. The population was created from an F1individual from a cross between a DH rapid cycling of Chinese kale (TO1000DH3, P1) and a DH broccoli line “Early Big” (P2) (Iniguez-Luy et al., 2009). Parents and 155 DH inbred lines were sown in a greenhouse under 16 h of daylight and a temperature of 24 ± 2°C, and 8 h of darkness with 18 ± 2°C at night, and a relative humidity of 55%. Plants were sown in a completely randomized experiment with two replications and four plants per replication and DH line. From each line, leaf samples were taken at the 4 leaves stage and flower buds were taken differentially depending on the flowering time of each plant. One bulk was taken from each replication by mixing the four samples of leaves and flower buds. Samples were immediately conserved at −80°C, and afterwards they were lyophilized for 48 h (Christ® Beta 2-8 LD Plus Freeze Dryer, Germany). The dried material was powdered using an IKA-A10 (IKA-Werke GmbH & Co.KG) mill and the powder was used for analysis. Besides, two bulks of 50 mg of seed for each line were prepared for phenolic analysis.
Fifty milligrams of each sample were extracted in 500 μL 70% methanol and sonicated during 1 h (model: 3510E-MTH, Bransonic®, Mexico) to facilitate the extraction. The suspensions were allowed to stand overnight at 4°C and afterwards they were sonicated again for 1 h. Then, samples were centrifugated at 3700 rpm for 15 min. Supernatants were recovered and were filtered through 0.22 μm PTFE filters (Multi Screen®, Ireland).
Identification of Phenolic Compounds by HPLC-DAD-ESI/MSn Analysis
For phenolics identification, three samples from each parent at each plant organ were analyzed. Besides, we identified these compounds in samples from one inbreed line from the mapping population (TODH101). The UV spectra analysis was carried out on a Kinetex column (5 μm, C18, 100 A, 150 × 4.6 mm; Phenomenex, Macclesfield, UK). The mobile phase consisted of two solvents: water-acetic acid (1%) (A) and methanol (B), starting with 20% B and using a gradient to obtain 40% B at 20 min and 60% B at 30 min. The flow rate was 1 mL/min, and the injection volume was 20 μL. Spectral data from all peaks were accumulated in the range of 240–400 nm, and chromatograms were recorded at 330 nm. The HPLC-DAD-ESI/MSn analyses were carried out in an Agilent HPLC 1100 series equipped with a diode array detector and mass detector in series (Agilent Technologies, Waldbronn, Germany). The HPLC consisted of a binary pump (model G1312A), an auto sampler (model G1313A), a degasser (model G1322A), and a photodiode array detector (model G1315B). The HPLC system was controlled by ChemStation software (Agilent, v. 08.03). The mass detector was an ion trap spectrometer (model G2445A) equipped with an electro spray ionization interface and was controlled by LCMSD software (Agilent, v. 4.1). The ionization conditions were adjusted at 350°C and 4 kV for capillary temperature and voltage, respectively. The nebulizer pressure and flow rate of nitrogen were 65.0 psi and 11 L/min, respectively. The full scan mass covered the range from m/z 100 up to m/z 1600. Collision-induced fragmentation experiments were performed in the ion trap using helium as the collision gas, with voltage ramping cycles from 0.3 up to 2 V. Mass spectrometry data were acquired in the negative ionization mode. Individual phenolic compounds were identified based on the data obtained from the standard substances or published literature including tR, λmax, ([M–H]−), and major fragment ions (Vallejo et al., 2004; Ferreres et al., 2006; Francisco et al., 2009; Velasco et al., 2011).
Quantification of Phenolic Compounds by UHPLC
Quantification of phenolic compounds was carried out in the two parents and in the 155 DH inbred lines of the mapping population. An Ultra-High-Performance Liquid-Chromatograph (UHPLC Nexera LC-30AD; Shimadzu) equipped with a Nexera SIL-30AC injector and one SPD-M20A UV/VIS photodiode array detector was used for the quantification. The UHPLC column was a Kinetex™ 2.6 μm C18 82–102 Å, LC Column 100 × 4.6 mm, protected with a C18 guard cartridge. The flow rate was 0.4 mL/min and the oven temperature was set at 30°C. The mobile phase consisted of two solvents: water-acetic acid (1%) (A) and methanol (B), starting with 10% B and using a gradient to obtain 40% B at 15 min and 60% B at 24 min. The injection volume was 5 μL. Chromatograms were recorded at 330 nm and data was processed on a computer with the LabSolutions software (Shimadzu). Caffeoyl quinic and p-coumparoyl quinic acids derivatives were quantified as chlorogenic acid (Sigma–Aldrich Chemie GmbH, Steinheim, Germany), flavonoids as kaempferol-3-O- glucoside and sinapic acid and derivatives as sinapic acid (Sigma–Aldrich Chemie GmbH, Steinheim, Germany).
A combined analysis of variance across organs and individual analyses of variance for each organ were made for individual and total phenolics. Lines and organs were considered as fixed factors and replications were considered as random factor. Comparisons between parents for each of the studied traits were accomplished with the Student's t-test. All the statistical analyses were performed with the PROC GLM of SAS (SAS Institute Inc., 2011).
The genetic map employed for the QTL analysis was created by Iniguez-Luy et al. (2009). It has 279 markers (SSRs and RFLPs) distributed along nine linkage groups assigned to chromosomes C1–C9 of B. oleracea with a total distance of 891.4 cM and a marker density of 3.2 cM/marker. Quantitative trait locus mapping was carried out through a composite interval mapping method (Zeng, 1994) by using the software PLABQTL (version 1.2) (Utz, 2003). Individual analyses were carried out for each trait and plant developmental stage (leaf, flower bud and seed). Empirical thresholds were build with permutation tests, with N = 1000 (Van Ooijen, 1999). The confidence intervals were set to 95%. The analysis and cofactor election were carried out by following PLABQTL's recommendations. The proportion of phenotypic variance explained for a specific trait was determined by the adjusted coefficient of determination of regression (R2) fitting a model which includes all detected QTLs (Papst et al., 2004). Fivefold cross-validation of QTLs was performed by following the procedures described by Utz et al. (2000). The frequency of QTL detection gives us an estimation of the precision of QTL localization. Significant QTLs for phenolic compounds were integrated by using a QTL meta-analysis with BioMercator 4.2 software in order to give consensus QTLs (Goffinet and Gerber, 2000). An Akaike-type statistical criterion (AIC value) indicated the model which best fitted the data, including the number and the consensus QTLs positions. The aim of performing a meta-analysis was to find if a genomic region could determine the content of different phenolic compounds.
Candidate Gene Identification
In order to identify candidate genes which may directly account for QTLs in B. oleracea we used two different approaches, by in silico search in the B. oleracea available sequences and by genome comparison with A. thaliana through analysis of syntenic regions. DNA sequences codes of the SSRs and RFLPs markers of the BolTBDH mapping population were collected from http://www.brassica.info/resource/markers.php and DNA sequences were obtained from NCBI data base. A BLAST search at the B. oleracea Genomics Project web site, Bolbase was done (http://www.ocri-genomics.org/bolbase/index.html) using the DNA sequences of the markers. Two sequences were declared putative homologous if the E value of their alignment was = 10−4. Then, genes related to phenylpropanoid pathway were searched in those regions corresponding to the confidence interval of each QTL (Table S1). This approach was not possible to use for all the identified QTLs since the B. oleracea genome for download is yet incomplete. Then we also perform a comparative mapping study of B. oleracea and A. thaliana. Iniguez-Luy et al. (2009) identified collinear genomic blocks between the mapping population used in this study (BolTBDH) and A. thaliana by using a synteny analysis. Following this approach, we tried to locate genes involved in phenylpropanoid metabolism in A. thaliana which were obtained from TAIR (The Arabidopsis Information Resource) on the BolTBDH map by in silico mapping (Table S2).
Metabolite Profiling and Phenotype Variation in the B. oleracea BolTBDH Mapping Population
A total of 36, 38, and 7 phenolic compounds, including flavonoids and hydroxycinnamic acids, were identified in the leaves, flower buds and seeds of the BolTBDH mapping population, respectively (Figure 2). The flavonoids were identified as O-glycosides containing a substituent at 3− and/or 7− positions of quercetin, kaempferol, and isorhamnetin and/or in conjugation with different hydroxycinnamic acids. Acylation were most frequently at the 3-position. The identified hydroxycinnamic acids included caffeoyl and p-coumaroyl quinic acids, feruloyl glycoside, and different sinapoyl derivatives glycosilated with at least two sugar moieties (gentiobiosides).
Figure 2. Phenolic profile in leaves, flower buds and seeds of a double haploid population of Brassica oleracea; P1 (TO1000DH3) and P2 (“Early Big” broccoli). See Table S4 for compounds abbreviations.
The major individual phenolics in each plant organ under study were then quantified in both parents and in the 155 inbred lines of the DH population. Fifteen individual phenolic compounds were quantified in the leaves samples, 15 in the flower buds, and five in the seeds. The mean of each individual phenolic compound quantified in the two parents and in 155 lines from the mapping population at the three plant organs under study are shown in Table 1. Besides individual compounds, total flavonoids, total hydroxycinnamic acids and total phenolics contents were also computed. At each plant stage, the two parents exhibited a similar phenolic profile but differed on their content. In leaves and flower buds, total phenolic content was higher in P2 than in P1. In contrast, P1 showed higher phenolics content than P2 in the seeds.
Table 1. Individual and total phenolics content of parents and in the 155 lines of BolTBDH mapping population in three different plant organs.
The frequency distributions of individual compounds as well as of total flavonoids and total hydroxycinnamic acids content in three plant organs were relatively normal (Figure 3). Transgressive distributions were observed for almost all the individual phenolics quantified in leaves and seeds, meaning that values for the DH mapping population were substantially larger and/or smaller than both of the parents. For instance, the content of the flavonoid F25 in leaves was 0.26 μmol g−1 dw in P1 and 0.12 μmol g−1 dw in P2. The average content of this kaempferol derivative in the mapping population was 1.36 μmol g−1 dw and ranged from 0 to 8.16 μmol g−1 dw (Table 1).
Figure 3. Frequency distributions of total hydroxycinnamic acids and total flavonoids in the 155 lines of the BolTBDH mapping population and their parental lines in three plant organs (leaves, flower buds, and seeds). The vertical axis of each figure represents the number of DH lines. The parental values are the mean of two replicates, indicates as P1 (TO1000DH3) and P2 (“Early Big” broccoli).
Phenolic Composition Differs Across Plant Organs
The phenolic profile and content of the mapping population varied depending on the plant organ analyzed (Figure 4). The average content of total phenolics in the population was 35.61 μmol g−1 dw in leaves, 53.49 μmol g−1 dw in flower buds and 75.09 μmol g−1 dw in seeds. In leaves, 26% of total phenolics content was due to flavonoids and 73% to hydroxycinnamic acids. In flower buds, 36% of total phenolics content was due to flavonoids and 64% to hydroxycinnamic acids. On the other hand, seeds showed a completely different phenolic profile with total hydroxycinnamic acids accounting for almost 100% of the total phenolic content.
Figure 4. Phenolic composition differs across plant organs in Brassica oleracea. Average content of flavonoids (kaempferol and isorhamnetin derivatives) and hydroxycinnamic acids (sinapoyl, p-coumaroyl, and caffeoyl derivatives) in leaves, flower buds, and seeds of the 155 lines from BolTBDH mapping population.
Across plant organs, flavonoids reached the higher content at flower bud stage (19.02 μmol g−1 dw) rather than leaves (9.28 μmol g−1 dw). Seeds showed very low quantities of flavonoids, less than 1 μmol g−1 dw. The flavonoid profile was also different between plant organs (Figure 2; Table 1). None of the identified flavonoids was coincident between leaves and flower buds. In leaves, all the quantified flavonoids were kaemperol-3-O-di/triglucoside-7-O-glucoside acylated with one or more hydroxycinnamic acids (F1, F2, F5, F7, F8, F25, F34, and F37). In flower buds, the major flavonoids were kaempferol-3-O-diglucoside, kaempferol-3-O-glucoside-7-O-glucoside and isorhamnetin-3-O-diglucoside-7-O-glucoside, and in contrast with the flavonoids found in leaves, none of them was in conjugation with hydroxycinnamic acids (F11, F15, F16, F18, F23, F28, and F44). Besides, isorhamnetin derivatives were not present in leaves. Flavonoids of quercetin were identified in both plant organs, but as they appeared in small quantities were not quantified.
Total hydroxycinnamic acids reached the highest average levels in seeds (74.72 μmol g−1 dw) and they were mainly represented by only one compound; a putative sinapoyl derivative (SA2) (Table 1). This compound was found exclusively in this plant organ, representing 91% of total phenolic content. The hydroxycinnamic acids profile was similar between leaves and flower buds although differed by content. The sinapoyl gentiobiosides derivatives SA9, SA10, and SA11 were the major phenolic acids in both plant organs. The average content of SA9 was higher in leaves while the content of SA10 and SA11 was higher in flower buds (Table 1).
Organ-Specific QTLs and Phenylpropanoid Pathway Candidate Genes
To investigate the genetic control of quantitative variation of phenolics accumulation in leaves, flower buds and seeds of the DH population of B. oleracea, a QTL analysis was performed. Detailed information of identified QTLs is given in Table S3 and Figure 5. Seventy nine significant QTLs were detected for the 29 traits analyzed which included the individual phenolics quantified as well as the sum of flavonoids (FlavT), the sum of hydroxycinnamic acids (HydroxT), and total phenolics (PhenolT) in each plant organ. The QTLs were spread all over the 9 linkage groups. These QTLs were not equally distributed over the B. oleracea linkage groups, as hotspots for the genetic control of metabolite content could be identified (Figure 5). The number of QTLs by linkage group ranged between two in C1 and 18 in C9. QTLs were detected for all the traits analyzed unless for F2 and F44. In leaves, 30 significant QTLs controlling phenolics content were found. The value of R2 ranged between 6.23% for F37 in C1 and 46.05% for 3CQAc in C4. Seventeen QTLs had a frequency of cross-validation higher than 50% (Table S3). In flower buds, 33 QTLs were found. The value of R2 ranged between 1.32% for SA11 in C5 and 36.79% for 5CQAc in C4. Fifteen QTLs had a frequency of cross-validation higher than 50% (Table S3). Sixteen QTLs were found in seeds. The value of R2 ranged between 3.18% for total hydroxycinnamic acids in C9 and 30.74% for SA1 in C3. Five QTLs had a frequency of cross-validation higher than 50% (Table S3).
Figure 5. Linkage map of the BolTBDH mapping population. Marker locus names are on the right side of each linkage group and distances (cM) are shown for each marker interval. Colors indicate the organ where the QTL was found: green for leaves; blue for flower buds; red for seeds. See Table S4 for compounds abbreviations.
The meta-analysis successfully reduced the total QTL by 58%. In total, 33 independent consensus QTL were identified at the three plant organ under study spread all over the nine linkage groups (Table 2). There were six on C3; five on C2 and C9; four on C8; three on C4, C5, C6, and C7; and two on C1.
Table 2. Description of organ-specific consensus quantitative trait loci (QTL) found in BolTBDH mapping population.
Twenty three of the consensus QTLs were plant organ-specific, 13 were found exclusively in flower buds, six in leaves, and four in seeds (Table 2). Some of them were also phenolic kind specific, i.e., QTL-2.5 and QTL-6.1 regulated only flavonoid content in leaves. Moreover, there were found genomic regions were the same QTL regulate the content of more than one compound in a specific organ. In leaves, the consensus QTL-9.3 was found to be responsible for the 8–15% of the phenotypic variation of six phenolic traits. In this region of B. oleracea genome a cluster of three putative hydroxycinnamoyl-Coenzyme A shikimate/quinate hydroxycinnamoyltransferases (HCT) was located (Table S1). In flower buds, the consensus QTL-6.3, QTL-7.1, and QTL-9.4 explained 8–12% of the phenotypic variation for different traits including total phenolic content. In the confidence interval of QTL-7.1 two candidate genes codifying for two UDP-glucosyltransferases were found (Table S1). In seeds, the consensus QTL-8.4 was exclusively related with the accumulation of the major hydroxycinnamic acids presents in this plant organ and was responsible for 11–18% of these compounds phenotypic variation.
QTL-7.3 and QTL-9.5, overlapped between the three plant organs. These QTLs were related with up to seven traits, including individual flavonoids and hydroxycinnamic acids as well as total phenolics (Table 2). These QTLs explained 6–30% of the phenolics accumulation among the three plant organs. A putative caffeoyl CoA O-methyltransferase was proposed as candidate gene for QTL-7.3 phenotypic variation (Table S1).
Leaves and flower buds showed other three co-located consensus QTLs in common in the linkage groups C4, C7, and C8 (Table 2). The QTL-4.1 was exclusively related with the accumulation of the chlorogenic acids (3CQAc, 5CQAc) and 3-p-coumaroyl quinic acid (3pCoQAc). This QTL was found to be responsible for a large proportion of phenotypic variation, 36–46%. Alleles for increasing the content of the chlorogenic acids are given by P1, while alleles for increasing the content of the 3pCoQAc are given by P2. In the confidence interval of this QTL we found synteny with chromosome 2 of A. thaliana. In this region of the model plant was found a key gene of the phenylopropanoid pathway (AT2g40890) which codifies a coumarate 3-hydroxylase enzyme (C3H) (Table S2).
Other consensus QTLs co-located between plant organs. The QTL-2.2, QTL-3.5, and QTL-4.2 were coincident for leaves and seeds. The QTL-2.3 and QTL-8.3 were related with the accumulation of different phenolics in both flower buds and seeds (Table 2).
The phenylpropanoid pathway serves as a rich source of metabolites in plants, being required for the biosynthesis of lignin, and serving as a starting point for the production of many other important compounds, such as the flavonoids and hydroxycinnamic acids. These compounds perform a variety of functions in the plant, generally centered on responses to pathogen attach and UV protection, as well as contributing toward the color and sensory characteristics of vegetables (Pereira et al., 2009). Besides, diets rich in foods containing phenolic compounds, such as B. oleracea crops, have been reported to possess many useful properties for human health including anti-inflammatory, enzyme inhibition, antimicrobial, antiallergic, vascular and cytotoxic antitumor activity, but the most important action of phenolics is their antioxidant activity (Crozier et al., 2009; Jahangir et al., 2009). In recent years, there is growing interest to increase the content of phytonutrients in crop species. Both traditional plant breeding and biotechnology-based techniques involving modifications of key enzymes catalyzing the synthesis of target compounds are needed to produce plants with the desired quality traits. Therefore, understanding phenolics regulation is one of the great challenges of modern plant biology. In the present work, we identified genomic regions controlling individual and total phenolic compounds accumulation in the B. oleracea TO1000DH3 × “Early Big” broccoli DH population. To gain insight into the biological relationship between the metabolism of phenolic compounds and plant development, the accumulation of phenolic compounds were comprehensively studied in three different organs (leaves, flower buds and seeds). To the best of our knowledge, this the first report dissecting genomic regions regulating individual and total phenolic compounds in different plant organs in B. oleracea. Moreover, for some of the described QTLs it was possible to propose candidate genes involved in phenylpropanoid metabolism that could control the phenotypic variation of the QTL.
Metabolite Profiling and Phenotype Variation in the B. oleracea BolTBDH Mapping Population
The main flavonoids identified were kaempferol, isorhamnetin and quercetin derivatives. The chemical structure more common in our samples were flavonol-3-O-(cinnamoyl)glycoside-7-O-glycoside derivatives. On the other hand, the predominant hydroxycinnamic acids have been identified as 3CQAc, 5CQAc, 3pCoQAc and different sinapic acid derivatives. Similar phenolic profile has already been characterized in other Brassica vegetables such as broccoli, turnips, kale and tronchuda cabbage (Vallejo et al., 2004; Ferreres et al., 2006; Francisco et al., 2009; Velasco et al., 2011).
Quantification of the major B. oleracea phenolic compounds revealed that these secondary metabolites varied greatly in concentration among lines in the DH mapping population in the three different plant organs under study. Distributions of individual and sums of phenolics were in most of cases transgressive. Therefore, new combination of additive alleles or epistatic interactions among loci for phenolics accumulation could be the cause to exceed the parental phenotypes. These types of segregation have been described before in the same DH population for GLS compounds (Sotelo et al., 2014b), and in other B. oleracea mapping population in a flavonoid genomic study (Lee et al., 2014).
Phenolic Composition Differs Across Plant Organs
The accumulation of some kind of phenolics was distributed in an organ-dependent manner. This different metabolic composition among organs may depend on de novo biosynthesis and catabolism or can also result from re-allocation of compounds. The ability to synthesize different phenolic compounds during plant developmental stages seems to be related with specific needs, allowing plants to cope with the constantly changing environmental challenges over their life cycle, thus obtaining protective adaptation to environmental stresses.
In this respect, the high diversity in the profile of flavonoids found in leaves and flower buds is not surprising. In each plant organ more than 20 structurally different flavonoids were found. On the other hand, seeds had less diversity, only one flavonoid in small quantities was found. In concordance with that, it has been described that flavonoid content and antioxidant activity in Brassica crops might differ completely form young to mature leaves as well as between external to internal leaves (Ferreres et al., 2006; Soengas et al., 2012). This particular accumulation of flavonoids suggests specialized functions in different stages of the plant development. It has been reported that flavonoids could have a role as internal physiological regulators or chemical messengers within the plant (Buer et al., 2010). Quercetin and kaempferol flavonoids can act through a specific receptor (NPA) in the plant cell plasma membrane blocking the polar auxin transport which could influence plant architecture (Brown et al., 2001; Buer and Muday, 2004). Besides, flavonoids have different roles in plant defense against pathogens, herbivores, and stress-related oxidative pressure (Buer et al., 2010). Therefore, any plant organ or temporal variation in the quality and quantity of flavonoids may have a role in plant organ differentiation, development, defense and adaptation.
In seeds, total phenolics content was mainly represented by only one hydroxycinnamic acid, a synapoyl derivative that was exclusively found in this plant organ. Hydroxycinnamic acids are precursors of lignin biosynthesis, important in the first plant stages to rigidifying cell walls and rendering them impermeable to water. Seeds of many species of Brassica genus accumulate high amount of sinapic acid derivaties such as O-sinapoylcholine (sinapine) and synapoyl gentiobiosides as the predominating phenolic compounds (Milkowski et al., 2004). During seed germination these compounds are transformed into sinapoylmalate which is accumulated in vacuoles of the leaf epidermal cell layer serving as UV screen. This could explain why seeds showed higher concentration of these hydroxycinnamic acids. So far, the biological role of synapoyl derivatives accumulation during seed development has not been elucidated, although it has been shown that overexpression of genes involved in sinapine metabolism influence the metabolism, morphology and physiology of developing seedlings (Clauß et al., 2011).
Organ-Specific QTLs and Phenylpropanoid Pathway Candidate Genes
In this study, we found 79 QTLs controlling phenolic content at three different organs. As it was expected, most of the detected QTLs were found in leaves and flower buds as they showed the highest variability. The meta-analysis reduced the total QTL number of independent consensus QTLs to 33 and they were spread all over the nine linkage groups. Comparative analysis of the consensus QTLs among plant organs revealed that only two of the consensus QTLs were located in the same genomic region for the three organs under study; while 23 of them were found within a single plant organ. Thus, most of the studied compounds had organ-specific genomic control, probably due to a selective expression of distinct fractions of the genome in response to developmental and environmental cues.
It has been shown that organ type is the most important factor in determining differentially expressed genes in A. thaliana (Schmid et al., 2005). The consensus QTL-9.3 seems to control most of the phenolics traits only in leaves. In silico mapping approaches located two candidate genes in this region, a putative HCT and a CCR (Table S1). These genes are in the core pathway of the phenylpropanoid metabolism influencing the accumulation of flavonoids, chlorogenic and sinapic acids derivatives in A. thaliana (Fraser and Chapple, 2011) (Figure 1). In flower buds, the consensus QTL-6.3, QTL-7.1, and QTL-9.4 regulated the accumulation of different traits including total phenolic content, individual flavonoids and hydroxycinnamic acids content. Two putative UDP-glucosyltransferases seem to be involved on the phenotypic variation related with QTL-7.1. In seeds, the QTL-8.4 regulated the content of the major hydroxycinnamic acids exclusively in this plant organ; unfortunately it was not possible to do in silico mapping in this region. Our findings, suggest a different gene regulation in the accumulation of phenolic compounds depending on the developmental stage of the plant or environmental signals. Moreover, this could explain the different metabolic profiles of phenolics accumulation among organs and will help to elaborate the genetic mechanism of phenolic biosynthesis in B. oleracea.
Only in two genomic regions, at C7 and C9 linkage groups, co-located consensus QTLs were found for leaves, flower buds and seeds (QTL-7.3 and QTL-9.5). They were related with up to seven traits including individual flavonoids, hydroxycinnamics as well as total phenolic compounds and explained between 6 and 20% of the phenotypic variation for those traits. This suggests that the genes that underlie these consensus QTLs could be involved in the core pathway of phenolics biosynthesis. By in silico analysis within the B. oleracea genome it was proposed a candidate gene (Bol043270) that could explain the phenotypic variation of the QTL-7.3 (Table S1). This gene has O-methyltransferase activity and codifies for a putative CCoAMT. In A. thaliana, this enzyme catalyzes the methylation of caffeoyl-CoA to feruloyl-CoA (in vitro and in vivo) and 5-hydroxyferuloyl-CoA to sinapoyl-CoA (at least in vitro) and, together with COMT, are involved in redundant functions for lignin, flavonoids and sinapoyl malate biosynthesis (Do et al., 2007). Disrupting the expression of both CCoAOMT and COMT dramatically reduces G and S lignin and 3-O-methoxylated soluble metabolites in A. thaliana indicating a main role of those enzymes in the phenylpropanoid pathway (Raes et al., 2003). The other consensus QTL shared by the three organs under study, the QTL-9.5, mapped to regions previously identified to contain a major QTL for individual and total GLS content in the same B. oleracea population (Sotelo et al., 2014b). Unfortunately it was not possible to use the in silico approach to localize candidate genes in this region. The fact that a single QTL affects the content of different secondary metabolites suggests that there may be a regulatory loci underlying the QTL or that there is a tight linkage of distinct genes involved in those metabolites accumulation. Co-localization of different QTLs might also be a first indication that some loci have a pleiotropic effect, due to a common mechanistic basis. Although GLSs and phenylpropanoids are synthesized through distinct biosynthetic pathways and have unique functions, Kim et al. (2015) showed that there is crosstalk between the two pathways. Significant contribution of the region underlying QTL-9.5 to the regulation of the level of secondary metabolites in B. oleracea opens the possibility for application of metabolic engineering in Brassica crops. However, it is important to realize that correlation between different pathways can exist, preventing the identification of a hotspot for regulation of specific metabolites. In the present study, since the markers were not uniformly distributed, large gaps appeared with low marker density on the region of the QTL-9.5 implying that more markers should be developed among these gaps and the authenticity of those common QTL for secondary metabolism should be further clarified.
Consensus QTLs were also detected for a specific kind of compounds. Some of the QTLs related with individual hydroxycinnamic acids explained a large proportion of the observed variation suggesting that hydroxycinnamics accumulation in B. oleracea is under the genetic control of a few additive loci. The QTL-4.3 regulated caffeoyl and p-coumaroyl quinic acids content in leaves and flower buds. This QTL showed the highest LOD scores and explained a large proportion of the variation for those traits (Table S3). Presence of a specific allele in this region of C4 resulted in higher levels of p-coumaroyl quinic acids and lower levels of caffeoyl quinic acids within both organs. In the confidence interval of this QTL we found synteny with a region of A. thaliana that mapped for a key gene of the phenylpropanoid pathway which codifies C3H enzyme. The C3H is a cytochrome P450 which transfers another hydroxyl group to the 3-position of the aromatic ring in p-coumaric acid required for the production of caffeoyl quinic acids (Figure 1). It is described as one of the main enzymes for the cinnamic intermediates in lignin biosynthesis (Rosler et al., 1997; Berner et al., 2006). The position of this QTL on the top of C4 linkage group co-localized with a previous QTL described by Sotelo et al. (2014a) related with antioxidant capacity in this B. oleracea population. This finding strongly supports further validation by fine mapping of C3H gene in the B. oleracea genome.
The hypothesis that regulation of flavonoid accumulation is a complex trait affected by many small effect alleles is also supported by the fact that at least 10 consensus QTLs were exclusively related with this kind of phenolics, being most of them also organ-specific. Moreover, most of those QTLs explained less than 10% of the phenotypic variation for a specific flavonoid trait. While over 60 genes have now been characterized to be involved in flavonoid metabolism, we could not find any known genes in the flavonoid metabolic pathway by in silico approaches. Recently, a study on A. thaliana seed flavonoid QTLs carried out by Routaboul et al. (2012) had similar results meaning that most loci identified did not co-localize with any gene previously characterized in the phenylpropanoid pathway. These findings suggest that the regulation of flavonoid accumulation may involve higher order transcriptional control or other yet uncharacterized limiting step such as transporter activity.
Since plant phenolics can modulate essential biological processes in the plant and are involved in a variety of defense response, a better understanding of the genomic regions that regulate phenol accumulation in Brassica vegetables within plant organs would be useful in terms of plant physiological ecology. The availability of the whole B. oleracea genome sequence together with possible comparative alignment with the related model species A. thaliana will facilitate fine mapping and cloning of candidate genes underlying the desired QTL. In addition, this information is also of great importance to both consumers and producers, as it can accelerate breeding programs to enhance plant defense and health protecting capacity of Brassica crops while decreasing the antinutritional properties of others.
Conception and design of the work: MF, PV, and PS; acquisition, analysis, and data interpretation: MF, MA, FF, DM, PV, and PS; wrote the paper: MF and PS.
Research supported by the Spanish Ministry of Economy and Competitiveness through the project AGL2012-35539 supported in part by the European Union (FEDER funds).
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.
The authors thank Dr. Federico Iñiguez-Luy for supplying seeds of the DH population and Rosaura Abilleira, Salomé Ruíz and Juan Carlos Fernández for laboratory help and field work and Dr. Ana María Butrón for all their help with QTLs analysis.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2015.01240
Ayaz, F. A., Hayırlıoglu-Ayaz, S., Alpay-Karaoglu, S., Grúz, J., Valentová, K., Ulrichová, J., et al. (2008). Phenolic acid contents of kale (Brassica oleraceae L. var. acephala DC.) extracts and their antioxidant and antibacterial activities. Food Chem. 107, 19–25. doi: 10.1016/j.foodchem.2007.07.003
Baumert, A., Milkowski, C., Schmidt, J., Nimtz, M., Wray, V., and Strack, D. (2005). Formation of a complex pattern of sinapate esters in Brassica napus seeds, catalyzed by enzymes of a serine carboxypeptidase-like acyltransferase family? Phytochemistry 66, 1334–1345. doi: 10.1016/j.phytochem.2005.02.031
Berner, M., Krug, D., Bihlmaier, C., Vente, A., Müller, R., and Bechthold, A. (2006). Genes and enzymes involved in caffeic acid biosynthesis in the actinomycete Saccharothrix espanaensis. J. Bacteriol. 188, 2666–2673. doi: 10.1128/JB.188.7.2666-2673.2006
Bhattacharya, A., Sood, P., and Citovsky, V. (2010). The roles of plant phenolics in defence and communication during Agrobacterium and Rhizobium infection. Mol. Plant Pathol. 11, 705–719. doi: 10.1111/j.1364-3703.2010.00625.x
Brown, D. E., Rashotte, A. M., Murphy, A. S., Normanly, J., Tague, B. W., Peer, W. A., et al. (2001). Flavonoids act as negative regulators of auxin transport in vivo in Arabidopsis. Plant Physiol. 126, 524–535. doi: 10.1104/pp.126.2.524
Buer, C. S., and Muday, G. K. (2004). The transparent testa4 mutation prevents flavonoid synthesis and alters auxin transport and the response of Arabidopsis roots to gravity and light. Plant Cell 16, 1191–1205. doi: 10.1105/tpc.020313
Carbonell-Bejerano, P., Diago, M.-P., Martinez-Abaigar, J., Martinez-Zapater, J., Tardaguila, J., and Nunez-Olivera, E. (2014). Solar ultraviolet radiation is necessary to enhance grapevine fruit ripening transcriptional and phenolic responses. BMC Plant Biol. 14:183. doi: 10.1186/1471-2229-14-183
Cheynier, V., Comte, G., Davies, K. M., Lattanzio, V., and Martens, S. (2013). Plant phenolics: recent advances on their biosynthesis, genetics, and ecophysiology. Plant Physiol. Biochem. 72, 1–20. doi: 10.1016/j.plaphy.2013.05.009
Clauß, K., von Roepenack-Lahaye, E., Böttcher, C., Roth, M. R., Welti, R., Erban, A., et al. (2011). Overexpression of sinapine esterase BnSCE3 in oilseed rape seeds triggers global changes in seed metabolism. Plant Physiol. 155, 1127–1145. doi: 10.1104/pp.110.169821
Do, C.-T., Pollet, B., Thévenin, J., Sibout, R., Denoue, D., Barrière, Y., et al. (2007). Both caffeoyl Coenzyme A 3-O-methyltransferase 1 and caffeic acid O-methyltransferase 1 are involved in redundant functions for lignin, flavonoids and sinapoyl malate biosynthesis in Arabidopsis. Planta 226, 1117–1129. doi: 10.1007/s00425-007-0558-3
Ferreres, F., Sousa, C., Vrchovská, V., Valentão, P., Pereira, J., Seabra, R., et al. (2006). Chemical composition and antioxidant activity of tronchuda cabbage internal leaves. Eur. Food Res. Technol. 222, 88–98. doi: 10.1007/s00217-005-0104-0
Francisco, M., Cartea, M. E., Butrón, A. M., Sotelo, T., and Velasco, P. (2012). Environmental and genetic effects on yield and secondary metabolite production in Brassica rapa crops. J. Agric. Food Chem. 60, 5507–5514. doi: 10.1021/jf301070q
Francisco, M., Moreno, D. A., Cartea, M. E., Ferreres, F., García-Viguera, C., and Velasco, P. (2009). Simultaneous identification of glucosinolates and phenolic compounds in a representative collection of vegetable Brassica rapa. J. Chromatogr. A 1216, 6611–6619. doi: 10.1016/j.chroma.2009.07.055
Hüsken, A., Baumert, A., Strack, D., Becker, H. C., Möllers, C., and Milkowski, C. (2005). Reduction of sinapate ester content in transgenic oilseed rape (Brassica napus) by dsRNAi-based suppression of BnSGT1 gene expression. Mol. Breed. 16, 127–138. doi: 10.1007/s11032-005-6825-8
Iniguez-Luy, F., Lukens, L., Farnham, M., Amasino, R., and Osborn, T. (2009). Development of public immortal mapping populations, molecular markers and linkage maps for rapid cycling Brassica rapa and B. oleracea. Theor. Appl. Genet. 120, 31–43. doi: 10.1007/s00122-009-1157-4
Kim, J. I., Dolan, W. L., Anderson, N. A., and Chapple, C. (2015). Indole glucosinolate biosynthesis limits phenylpropanoid accumulation in Arabidopsis thaliana. Plant Cell 27, 1529–1546. doi: 10.1105/tpc.15.00127
Koyama, K., Numata, M., Nakajima, I., Goto-Yamamoto, N., Matsumura, H., and Tanaka, N. (2014). Functional characterization of a new grapevine MYB transcription factor and regulation of proanthocyanidin biosynthesis in grapes. J. Exp. Bot. 65, 4433–4449. doi: 10.1093/jxb/eru213
Lee, R. H., Malchev, I., Rajcan, I., and Kott, L. (2014). Identification of putative quantitative trait loci associated with a flavonoid related to resistance to cabbage seedpod weevil (Ceutorhynchus obstrictus) in canola derived from an intergeneric cross, Sinapis alba × Brassica napus. Theor. Appl. Genet. 127, 419–428. doi: 10.1007/s00122-013-2228-0
Milkowski, C., Baumert, A., Schmidt, D., Nehlin, L., and Strack, D. (2004). Molecular regulation of sinapate ester metabolism in Brassica napus: expression of genes, properties of the encoded proteins and correlation of enzyme activities with metabolite accumulation. Plant J. 38, 80–92. doi: 10.1111/j.1365-313X.2004.02036.x
Nair, R. B., Joy, R. W., Kurylo, E., Shi, X., Schnaider, J., Datla, R. S. S., et al. (2000). Identification of a CYP84 family of cytochrome P450-dependent mono-oxygenase genes in Brassica napus and perturbation of their expression for engineering sinapine reduction in the seeds. Plant Physiol. 123, 1623–1634. doi: 10.1104/pp.123.4.1623
Pang, Y., Peel, G. J., Sharma, S. B., Tang, Y., and Dixon, R. A. (2008). A transcript profiling approach reveals an epicatechin-specific glucosyltransferase expressed in the seed coat of Medicago truncatula. Proc. Natl. Acad. Sci. U.S.A. 105, 14210–14215. doi: 10.1073/pnas.0805954105
Papst, C., Bohn, M., Utz, H. F., Melchinger, A. E., Klein, D., and Eder, J. (2004). QTL mapping for European corn borer resistance (Ostrinia nubilalis Hb.), agronomic and forage quality traits of testcross progenies in early-maturing European maize (Zea mays L.) germplasm. Theor. Appl. Genet. 108, 1545–1554. doi: 10.1007/s00122-003-1579-3
Raes, J., Rohde, A., Christensen, J. H., Van de Peer, Y., and Boerjan, W. (2003). Genome-wide characterization of the lignification toolbox in Arabidopsis. Plant Physiol. 133, 1051–1071. doi: 10.1104/pp.103.026484
Rezaeizad, A., Wittkop, B., Snowdon, R., Hasan, M., Mohammadi, V., Zali, A., et al. (2011). Identification of QTLs for phenolic compounds in oilseed rape (Brassica napus L.) by association mapping using SSR markers. Euphytica 177, 335–342. doi: 10.1007/s10681-010-0231-y
Routaboul, J.-M., Dubos, C., Beck, G., Marquis, C., Bidzinski, P., Loudet, O., et al. (2012). Metabolite profiling and quantitative genetics of natural variation for flavonoids in Arabidopsis. J. Exp. Bot. 63, 3749–3764. doi: 10.1093/jxb/ers067
Sotelo, T., Soengas, P., Velasco, P., Rodríguez, V. M., and Cartea, M. E. (2014b). Identification of metabolic QTLs and candidate genes for glucosinolate synthesis in Brassica oleracea leaves, seeds and flower buds. PLoS ONE 9:e91428. doi: 10.1371/journal.pone.0091428
Sousa, C., Pereira, D. M., Pereira, J. A., Bento, A., Rodrigues, M. A., Dopico-García, S., et al. (2008). Multivariate analysis of tronchuda cabbage (Brassica oleracea L. var. costata DC) phenolics: influence of fertilizers. J. Agric. Food Chem. 56, 2231–2239. doi: 10.1021/jf073041o
Takos, A. M., Ubi, B. E., Robinson, S. P., and Walker, A. R. (2006). Condensed tannin biosynthesis genes are regulated separately from other flavonoid biosynthesis genes in apple fruit skin. Plant Sci. 170, 487–499. doi: 10.1016/j.plantsci.2005.10.001
Tamagnone, L., Merida, A., Parr, A., Mackay, S., Culianez-Macia, F. A., Roberts, K., et al. (1998). The AmMYB308 and AmMYB330 transcription factors from antirrhinum regulate phenylpropanoid and lignin biosynthesis in transgenic tobacco. Plant Cell 10, 135–154. doi: 10.1105/tpc.10.2.135
Utz, H. F. M. A. E. (2003). PLABQTL, a Computer Program to Map QTL, Version 1.2. 2006-06-01. Stuttgart, Institute of Plant Breeding, Seed Science, and Population Genetics. Stuttgart: University of Hohenheim.
Utz, H. F., Melchinger, A. E., and Schön, C. C. (2000). Bias and sampling error of the estimated proportion of genotypic variance explained by quantitative trait loci determined from experimental data in maize using cross validation and validation with independent samples. Genetics 154, 1839–1849.
Vallejo, F., Tomás-Barberán, F. A., and Ferreres, F. (2004). Characterisation of flavonols in broccoli (Brassica oleracea L. var. italica) by liquid chromatography-UV diode-array detection-electrospray ionisation mass spectrometry. J. Chromatogr. A 1054, 181–193. doi: 10.1016/j.chroma.2004.05.045
Velasco, P., Francisco, M., Moreno, D. A., Ferreres, F., García-Viguera, C., and Cartea, M. E. (2011). Phytochemical fingerprinting of vegetable Brassica oleracea and Brassica napus by simultaneous identification of glucosinolates and phenolics. Phytochem. Anal. 22, 144–152. doi: 10.1002/pca.1259
Verdier, J., Zhao, J., Torres-Jerez, I., Ge, S., Liu, C., He, X., et al. (2012). MtPAR MYB transcription factor acts as an on switch for proanthocyanidin biosynthesis in Medicago truncatula. Proc. Natl. Acad. Sci. U.S.A. 109, 1766–1771. doi: 10.1073/pnas.1120916109
Keywords: Brassica oleracea, phenylpropanoid metabolism, QTL mapping, candidate gene, flavonoid, hydroxycinnamic acid
Citation: Francisco M, Ali M, Ferreres F, Moreno DA, Velasco P and Soengas P (2016) Organ-Specific Quantitative Genetics and Candidate Genes of Phenylpropanoid Metabolism in Brassica oleracea. Front. Plant Sci. 6:1240. doi: 10.3389/fpls.2015.01240
Received: 23 October 2015; Accepted: 20 December 2015;
Published: 28 January 2016.
Edited by:Antonio Encina García, University of León, Spain
Reviewed by:Akifumi Sugiyama, Kyoto University, Japan
Jens Rohloff, Norwegian University of Science and Technology, Norway
Copyright © 2016 Francisco, Ali, Ferreres, Moreno, Velasco and Soengas. 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: Pilar Soengas, firstname.lastname@example.org