Next Article in Journal
BMS-986020, a Specific LPA1 Antagonist, Provides Neuroprotection against Ischemic Stroke in Mice
Next Article in Special Issue
Does Herbal and/or Zinc Dietary Supplementation Improve the Antioxidant and Mineral Status of Lambs with Parasite Infection?
Previous Article in Journal
Non-Coding RNAs as Sensors of Oxidative Stress in Neurodegenerative Diseases
Previous Article in Special Issue
Chemico-Biological Characterization of Torpedino Di Fondi® Tomato Fruits: A Comparison with San Marzano Cultivar at Two Ripeness Stages
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Changes in Phenolics and Fatty Acids Composition and Related Gene Expression during the Development from Seed to Leaves of Three Cultivated Cardoon Genotypes

1
Department of Pharmacy, University of Naples Federico II, Via Domenico Montesano 49, 80131 Naples, Italy
2
Institute of Bioscience and Bioresources, Consiglio Nazionale delle Ricerche, via Università 133, 80055 Portici, Italy
3
Institute of Agricultural Biology and Biotechnology, Consiglio Nazionale delle Ricerche, Via E. Bassini 15, 20133 Milan, Italy
4
Unesco Chair for Health Education and Sustainable Development, 80131 Naples, Italy
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Antioxidants 2020, 9(11), 1096; https://doi.org/10.3390/antiox9111096
Submission received: 19 October 2020 / Revised: 2 November 2020 / Accepted: 5 November 2020 / Published: 8 November 2020
(This article belongs to the Special Issue The Potential of Dietary Antioxidants)

Abstract

:
Cultivated cardoon (Cynara cardunculus var. altilis) has long been used as a food and medicine remedy and nowadays is considered a functional food. Its leaf bioactive compounds are mostly represented by chlorogenic acids and coumaroyl derivatives, known for their nutritional value and bioactivity. Having antioxidant and hepatoprotective properties, these molecules are used for medicinal purposes. Apart from the phenolic compounds in green tissues, cultivated cardoon is also used for the seed oil, having a composition suitable for the human diet, but also valuable as feedstock for the production of biofuel and biodegradable bioplastics. Given the wide spectrum of valuable cardoon molecules and their numerous industrial applications, a detailed characterization of different organs and tissues for their metabolic profiles as well as an extensive transcriptional analysis of associated key biosynthetic genes were performed to provide a deeper insight into metabolites biosynthesis and accumulation sites. This study aimed to provide a comprehensive analysis of the phenylpropanoids profile through UHPLC-Q-Orbitrap HRMS analysis, of fatty acids content through GC-MS analysis, along with quantitative transcriptional analyses by qRT-PCR of hydroxycinnamoyl-quinate transferase (HQT), stearic acid desaturase (SAD), and fatty acid desaturase (FAD) genes in seeds, hypocotyls, cotyledons and leaves of the cardoon genotypes “Spagnolo”, “Bianco Avorio”, and “Gigante”. Both oil yield and total phenols accumulation in all the tissues and organs indicated higher production in “Bianco Avorio” and “Spagnolo” than in “Gigante”. Antioxidant activity evaluation by DPPH, ABTS, and FRAP assays mirrored total phenols content. Overall, this study provides a detailed analysis of tissue composition of cardoon, enabling to elucidate value-added product accumulation and distribution during plant development and hence contributing to better address and optimize the sustainable use of this natural resource. Besides, our metabolic and transcriptional screening could be useful to guide the selection of superior genotypes.

1. Introduction

Cynara cardunculus L. var. altilis DC., the cultivated cardoon, belongs to the Asteraceae family and, with its sister species globe artichoke (Cynara cardunculus var. scolymus) and their common ancestor wild cardoon (Cynara cardunculus var. sylvestris), participates in the small Cynara genus. Cardoon and artichoke originate from the Mediterranean area [1], where they are mostly cultivated. Artichoke is famous as a food crop for its edible immature inflorescences (heads) and is also produced in America and Asia, while cardoon, although having a similar composition, has a restricted traditional food use in Mediterranean countries, though its cultivation for industrial purposes is increasing. Indeed, cultivated cardoon is gaining interest as a multipurpose crop, due to the versatility of its main products, which offer a wide spectrum of applications [2]. Due to its high biomass and related content of cellulose, emicellulose, and lignin, cardoon is starting to be exploited for energy production and green chemistry [3,4,5,6]. Concomitantly, cardoon extracts maintain a prominent role in the nutraceutical and pharmaceutical sectors, since numerous bioactive molecules are biosynthesized and accumulated in different parts of the plant. In particular, most of the health-beneficial compounds are polyphenols, which are present in all plant organs, though in higher amounts in leaves and seeds [7].
Cardoon biosynthesizes specialized metabolites, which are distinctive of Asteraceae [8], namely flavonoids, such as apigenin and luteolin, and hydroxycinnamic derivatives such as mono- and di-caffeoylquinic acids (CGA) [9,10,11]. Extraction of these cardoon bioactive compounds has been reported in many works [7,12,13] and several in vitro and in vivo studies have demonstrated the health-promoting effects related to its polyphenols. Among these polyphenols, chlorogenic acid (CGA) represents the major bioactive component because of its abundance and numerous bioactive functions, especially as a free-radical scavenger and antioxidant. These functions seem to be associated with a well-known dual role as a protectant against oxidative injury caused by free radicals [14,15] and as a substrate for both chemical and enzymatic browning reactions [9]. Apart from CGA, most cardoon phenolics confer several medicinal properties to cardoon leaf extracts, such as cholesterol-lowering, hepatoprotective, and diuretic activities, but also antifungal and antibacterial properties [16,17,18,19,20]. The biosynthesis of general phenylpropanoids and CGA proceeds through several sequential enzymatic reactions, starting from the amino-acid phenylalanine that, being the substrate of phenylalanine-ammonia-lyase (PAL), is the entry point for the formation of a plethora of phenylpropanoids. Afterward, cinnamate 4-hydroxylase (C4H) and 4-coumarate-CoA ligase generate p-coumaroyl-CoA, which is the precursor of both CGA and flavonoids. Three possible routes of CGA biosynthesis have been proposed in several plant species [21,22]. The route characterized in artichoke and other Asteraceae [23,24] proceeds through the catalytic activity of hydroxycinnamoyl-CoA quinate transferase (HQT) on caffeoyl-CoA and quinic acid to form CGA. Several HQT genes with catalytic activity towards p-coumaroyl-CoA have been characterized in Asteraceae species, and it cannot be excluded that several other HQT isoforms with a possible tissue-specific activity and specific regulation [25] are still undiscovered. In this regard, the study of the spatio-temporal expression pattern of HQT along with dynamic changes in CGA accumulation and distribution can provide insights into possible biosynthetic correlations and any developmental regulatory factors acting on the accumulation of this specialized metabolite.
Indeed, cardoon use is not limited to phenylpropanoids, but also oils. As oleaginous species, cardoon produces seed oils, which have been investigated both for human nutrition [26] and as a potential renewable source for the production of biofuels and bioplastics [27]. In plants, fatty acids and derived molecules represent an energy source but are also fundamental in many metabolic and signaling processes, other than being structural components of membranes and thus directly involved in defenses responses [28,29]. In human and animal nutrition, fatty acids such as oleic and linoleic are essential major nutrients since mammals cannot synthesize but only uptake them through the diet. In plants, fatty acids are synthesized in the plastids starting from acetyl-CoA, and, after passing into the cytosol, are modified and assembled in lipids in the endoplasmic reticulum (ER) [30,31]. The desaturation of stearic acid (C18:0) to oleic acid (C18:1) is catalyzed by stearoyl-acyl carrier protein desaturase (SAD). Further desaturation of oleic to linoleic acid (18:2) is catalyzed by FAD2 in the ER and FAD6 in the plastid [32]. As oleic (C18:1Δ9) and linoleic acids (C18:2Δ9,12) are the main determinants for oil quality, their content is very important to assess the nutritional properties of edible oils or direct their technological applications [33,34]. Oleic and linoleic acids have both the capability of lowering total serum cholesterols, but oleic acid, having one double bond less than linoleic, has higher oxidative stability. For this reason, high oleic acid content has also a great potential for industrial uses, for example for making products such as biodiesel and lubricants, needing high oxidative stability. In this regard, the knowledge of active genes in oleic and linoleic biosynthesis is important for selection of MUFA- (mono-unsaturated fatty acids) and/or PUFA- (poly unsaturated fatty acids) enriched tissues and genotypes.
In our study, we followed the distribution and accumulation of phenylpropanoids and oils during the growth and development of cultivated cardoon from seeds to leaf formation. We characterized the metabolic composition of seeds, hypocotyls, cotyledons, and leaves of the three cultivated cardoon genotypes “Bianco Avorio”, “Spagnolo”, and “Gigante” by GC and high-resolution mass spectrometric analysis, along with the tissue and developmental expression of key genes in the biosynthesis of chlorogenic acid and oleic and linoleic acid by qRT-PCR. To minimize the effects of environment and growth conditions, that are known to influence the level and composition of these metabolites [35], we grew the three cardoon genotypes in controlled greenhouse conditions, thus allowing to specifically highlight genotype differences and assist the selection of the most suitable tissues and genotypes for different purposes. Moreover, to characterize tissues and genotypes based on the biological properties, we analyzed the antioxidant activities by ABTS, DPPH, and FRAP for all the tissues.

2. Materials and Methods

2.1. Reagents and Materials

Polyphenolic standards (purity > 98%), including luteolin, apigenin, diosmin, apigenin-8-C-glucoside, naringenin, quercetin-3-O-glucoside, quercetin, kaempferol, myricetin, naringin, kaempferol-3-O-glucoside, luteolin-7-O-glucoside, p-coumaric acid, chlorogenic acid, and 4-hydroxybenzoic acid were purchased from Sigma-Aldrich Chemical Co. (St. Louis, MO, USA). The reference standards were accurately weighed and dissolved in methanol to prepare stock solutions (1 mg mL−1). For the construction of calibration curves, a multi-compound working solution was prepared by appropriate dilution of stock solutions with methanol–water (80:20 v/v) to make concentrations in the range 0.02–5 µg·mL−1. For the antioxidant tests, gallic acid, 6-hydroxy-2,5,7,8-tetramethylchromane-2-carboxylic acid (Trolox), 1,1-diphenyl-2-picrylhydrazyl (DPPH), 2,3,5-triphenyltetrazolium chloride (TPTZ), anhydrous ferric chloride, hydrochloric acid, potassium persulfate (K2S2O8), 2,2’-azino-bis(3-ethylbenzothiazoline-6-sulfonic acid) diammonium salt (ABTS) and sodium acetate were purchased from Sigma-Aldrich (MO, United States). Methanol (MeOH), ethanol, and water (LC-MS grade) were acquired from Carlo Erba reagents (Milan, Italy), whereas formic acid (98–100%) was purchased from Fluka (Milan, Italy).

2.2. Plant Growth and Sample Collection

Seeds of C. cardunculus L. var. altilis DC. “Gigante” and “Spagnolo” genotypes were provided by ARCA 2010 s.c.a.r.l., Acerra (Italy), and “Bianco Avorio” seeds were provided by La Semiorto Sementi, Lavorate di Sarno (Italy). Seeds were placed in germination pots in the greenhouse of the Institute of Bioscience and Bioresources, Portici, (Italy), and grown for the collection of the different tissues at different vegetative stages (Supplementary Figure S1). Part of the germinated seeds were used for the collection of hypocotyls and cotyledons. Hypocotyls were cut with a razor blade, approx. 5 mm below the cotyledons, about 7–10 days post-germination. Other germinated seeds were used for the collection of cotyledons, and leaves were collected 30 days after emergence. Seeds, hypocotyls, cotyledons, and leaves were flash frozen at −80 °C and stored for RNA extraction and biochemical characterization.

2.3. RNA Extraction and qRT-PCR

Seeds and tissues of cultivated cardoon were harvested in biological triplicates and ground in liquid nitrogen. Total RNA was extracted from 100 mg of tissue using the RNeasy kit (Qiagen, CA, USA). After verification of RNA integrity and NanoDrop ND-8000 spectrophotometer (Fisher Scientific, DE, USA) quantification, 1 μg of total RNA was reversed transcribed using the QuantiTec Reverse Transcription Kit (Quiagen, CA, USA), according to the manufacturer’s instructions. Real-time qRT-PCR was performed with Platinum® SYBR® Green qPCR SuperMix (Applied Biosystem, CA, USA) in an ABI7900 HT (Life Technologies, CA, USA). Each PCR reaction (20 μL) contained 10 μL real-time qRT-PCR Mix, 4 μL of a 1:25 dilution of cDNA, and 0.25 μM of each specific primer designed on sequences retrieved from NCBI. Thermal cycling conditions were as follows: 50 °C for 2 min, 95 °C for 2 min, followed by 40 cycles of 15 s at 95 °C and 30 s at 60 °C. The dissociation stage was performed for each primer pairs and a single melting curve for each analyzed gene indicated the amplification of a single band. All reactions were performed on biological triplicates and technical duplicates and the 2−ΔΔCT method [36,37] was used for fold change measurements. For the relative expression of each gene in each genotype, the tissue with the lowest expression level was used as an internal calibrator by setting its expression equal to 1. C. cardunculus actin gene, as previously reported [38] was used to normalize gene expression. In Supplementary Table S1, the accession number and sequences of primers used for real-time qRT-PCR are reported.

2.4. Ultrasound-Assisted Extraction of Polyphenolic Compounds

Lyophilized samples were extracted using the method reported in the literature [39] with few modifications. In particular, 3 g of dried sample were extracted with 30 mL of ethanol/water (50:50 v/v) by sonication at room temperature for 30 min. Samples were centrifuged to 4000 rpm at 4 °C, filtered through 0.45 mm nylon syringe membranes, and then used for high-resolution mass spectrometry analysis and antioxidant activity assay. The ultrasonic plant material extraction procedure was repeated three times.

2.5. Antioxidant Activity: ABTS Assay

The antioxidant capacity assay was conducted based on the method described by [40]. Briefly, 44 µL of aqueous 2.45 mM potassium persulfate were added to aqueous 7mM ABTS and incubated in the dark at room temperature (23 °C) for 12–16 h. After dilution (1:88) with ethanol, the absorbance of this ABTS+ working solution should be 0.700 ± 0.050 at 734 nm. The assay was performed by adding 0.1 mL of filtered and suitably diluted sample to 1 mL of ABTS+ working solution and the absorbance was monitored at 734 nm after 2.5 min. Results were expressed as Trolox® equivalent antioxidant capacity (TEAC, mmol Trolox® equivalents Kg-1 dry weight of plant). All determinations were performed in triplicate.

2.6. Antioxidant Activity: DPPH Assay

The DPPH radical-scavenging activity was determined using the method proposed by [41] with minor modifications. Briefly, the DPPH methanolic solution 100 μM (1 mL) was added to polyphenol extract (0.2 mL), and the decrease in absorbance of the resulting solution was monitored at 517 nm after 10 min. The results were corrected for dilution and expressed as Trolox® equivalent antioxidant capacity (TEAC, mmol Trolox® equivalents Kg−1 dry weight of plant). All determinations were performed in triplicate.

2.7. Antioxidant Activity: FRAP Assay

The FRAP assay was conducted according to the method reported by [42] with slight adaptations. Briefly, the FRAP reagent contained a solution of 10 μM TPTZ in 40 μM HCl, 20 μM of aqueous FeCl3 and acetate buffer (300 μM, pH 3.6) at 1:1:10 (v/v/v). The FRAP reagent (300 μL) and sample solutions (10 μL) were mixed and the absorbance was monitored at 593 nm after 10 min. The results were expressed in mmol Trolox® Kg−1 dry weight (dw). The results were corrected for dilution and expressed as Trolox® equivalent antioxidant capacity (TEAC, mmol Trolox® equivalents Kg−1 dry weight of plant). All determinations were performed in triplicate.

2.8. HRMS Orbitrap Analysis of Bioactive Polyphenols

An Ultra-High-Pressure Liquid Chromatograph (UHPLC, Dionex UltiMate 3000, Thermo Fisher Scientific, Ma, USA) coupled with a Q-Exactive Orbitrap mass spectrometer (UHPLC, Thermo Fischer Scientific, Ma, USA) was used to investigate the quali-quantitative profile of polyphenolic compounds applying conditions reported in our previous work [43]. A Kinetex 2.6 µm Biphenyl (100 × 2.1 mm, Phenomenex) column was applied for chromatographic separation of polyphenols with a column temperature set at 25 °C. The mobile phase consisted of water containing 0.1% of formic acid (eluent A) and methanol containing 0.1% of formic acid (eluent B). Polyphenolic compounds were eluted using the following gradient program: 0–1.3 min 5% B, 1.3–9.3 min 5–100% B, 9.3–11.3 min 100% B, 11.3–13.3 min 100–5% B, 13.3–20 min 5% B. The flow rate was 0.2 mL min−1 and the injection volume was 2 μL. The mass spectrometer was operated in negative ion mode (ESI-) setting two scan events (Full ion MS and All ion fragmentation, AIF) for all compounds of interest. Full scan data were acquired setting a resolving power of 35,000 FWHM (full width at half maximum) at m/z 200. The key parameters were as follows: spray voltage −2.8 kV, sheath gas flow rate 35 arbitrary units, auxiliary-gas flow rate, 10 arbitrary units, capillary temperature 310 °C, auxiliary gas heater temperature 350 °C, S-lens RF level 50. For the scan event of AIF, the resolving power was set at 17,500 FWHM, the collision energies were 10, 20, and 45 eV, and the scan range was m/z 80–1200. Data acquisition and processing were performed with Quan/Qual Browser Xcalibur software, v. 3.1.66.10 (Xcalibur, Thermo Fisher Scientific, (Thermo Fisher Scientific, Waltham, MA, USA). Parameters of UHPLC-HRMS Orbitrap validation are reported in Supplementary Table S2.

2.9. Oil Extraction

The oil from the samples was extracted using solvent extraction according to the procedure described in the literature [44]. In particular, 100 mL of hexane were added to 10 g of powdered sample and shaken in a conical flask for 24 h. Then the solution was filtered under reduced pressure through filter paper (Whatman No. 1) and concentrated to remove completely hexane on a rotary evaporator (Rotavapor RE 120; Büchi, Flavil, Sweden). The final weight of solid material left after evaporation was recorded and used for the fatty acid profile study.

2.10. Fatty Acids Analysis

Preparation of fatty acid methyl ester derivatives (FAMEs) was carried out dissolving 0.1 g of oil sample in 1 mL hexane and 500 μL of methanolic KOH (2N). Afterward, the mixture was vortexed for 2 min at room temperature and the upper layer (hexane) was directly injected into GCMS. The fatty acid composition was investigated by GC-MS (Agilent Technologies, Santa Clara, CA, USA) using a capillary column, Rxi-5MS 5% Phenyl 95% Dimethylpolysiloxane (29.7 m × 0.25 mm i.d., 0.25 µm) (Restek, Bellefonte, PA, USA). One µL of the sample was injected in splitless mode into the injection port. Helium (grade 6.0) was used as carrier gas with a flow rate of 1.4 mL min−1. The temperature program of the column was as follows: the initial oven temperature = 50 °C for 1.5 min then raising from 50 °C to 80 °C with a program ramp rate of 2.5 °C/min, the column was kept at 80 °C for 1 min then up to 300 °C at 10 °C min−1, making a total run time of 36 min. The injector temperature was 280 °C.
The mass spectrometer was tuned according to the manufacturer’s recommendation using perfluorotributylamine (PFTBA). The raw GC- TOF-MS data were processed by the ChromaTOF software version 5.03.09.0 (LECO Corporation, Saint Joseph, MI, USA) that deconvolutes mass spectra and performs peak identification using the NIST11 library (NIST, MD, USA). Results were expressed as mass response areas in relative percentages.

2.11. Statistical Analysis

Statistical analyses were performed using SigmaStat version 11.0 software (Sigma Stat, Statcon, Witzenhausen, Germany, SigmaStat for Windows). Mean and standard deviation were calculated on three biological and two technical replicates for all the biochemical and qRT-PCR data. Statistical differences were evaluated through a one-way analysis of variance (ANOVA). Tukey’s post hoc test was used for mean separation and the statistical significance of the comparisons was defined as p < 0.05.

3. Results and Discussion

3.1. Polyphenol Profiling of Cardoon Genotypes

Qualitative and quantitative profiles of polyphenolic compounds were obtained using an ultra-high-performance liquid chromatography method coupled with electrospray ionization hybrid linear trap quadrupole orbitrap mass spectrometry (UHPLC-Q-Orbitrap-HRMS). Data were analyzed using Qual Browser Xcalibur 3.0 (Thermo Fisher Scientific), and identification of individual compounds was supported by retention times, exact mass spectra data, and MS/MS spectra. Individual phenolic compounds were quantified using the calibration curves of the respective reference compounds as described in Section 2.1.
The bibliographic data indicated that the chemical profile of cardoon depends on the variety [45,46] and geographical origin [47]. The UHPLC technique coupled to Orbitrap-HRMS showed that the different analyzed cardoon tissues and organs contained the same phenolic qualitative profile, but their relative abundances were considerably different between samples and genotypes. A typical full-scan MS chromatogram of a cardoon sample (cotyledons) is reported in Supplementary Figure S2. The studied cardoon parts, namely seeds, hypocotyls, cotyledons, and leaves, showed significant differences in phenolic compounds content, with leaves having the highest total phenolic content and reaching an average value of 3.99 mg g−1 on a dry weight (dw) basis. Table 1 summarizes the peak characteristics of fifteen detected phenolic compounds. Two phenolic acids (peaks 1 and 3), nine non-anthocyanin flavonoids (peaks 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, and 15), and one caffeoylquinic acid (peak 2) were identified in C. cardunculus samples.
Flavonoids are represented by the greatest number of compounds, among which the most abundant were apigenin-8-C-glucoside (vitexin), luteolin glucoside, naringin, and luteolin, which is in accordance with the review by [48]. Peak 1 and 3 were positively identified as 4-hydroxy benzoic and coumaric acid, respectively, by comparing their retention time and mass characteristics with commercial standards. Regarding peak 2, caffeoylquinic acid, ([M − H] at m/z 353) was positively identified as 5-O-caffeoylquinic acid, yielding the base peak at m/z 191 and a secondary ion at m/z 179, characteristic of 5-acylchlorogenic acids, as reported by [49]. Peak 4 ([M − H] at m/z 463) released a unique MS2 unit at m/z 301 (−162 μ), corresponding to the loss of a glucosyl unit, being tentatively identified as quercetin 3-O-glucoside. Peak 5 was identified as apigenin 8-C-glucoside, presenting pseudomolecular ion [M − H] at m/z 431 and an MS2 fragment at m/z 269, corresponding to the loss of glucosyl unit. Peaks 8 and 10 were identified respectively as luteolin and kaempferol glucosides both showing pseudomolecular ion [M − H] at m/z 447 (MS2 unit at m/z 285). Peak 6 was identified as diosmin, characterized by the typical mass spectra in negative mode with the pseudomolecular ion [M − H] at m/z 607. Peaks 7, 9, 11, 12, 13, 14, and 15 were identified as naringin, myricetin, naringenin, quercetin, luteolin, kaempferol, and apigenin, presenting pseudomolecular ions [M − H] at m/z 579, 317, 271, 301, 285, 285, and 269, respectively.
The qualitative and quantitative profile of identified compounds (Table 2) varied between seeds, hypocotyls, cotyledons, and leaves of the three C. cardunculus analyzed genotypes. Vitexin (apigenin-8-C-glucoside) concentration in tissues and genotypes showed an increase during plant development in the three genotypes, reaching a maximum in leaves with the highest values in “Bianco Avorio” (2107.906 µg g−1) (dw), followed by “Spagnolo” (2100.344 µg g−1) (dw) (Table 2 and Supplementary Table S3). Similar behavior was also observed for luteolin 7-O-glucoside, quercetin 3-glucoside, luteolin, and kaempferol (Supplementary Table S3), possibly suggesting a developmental regulation of their production. In all tissues and genotypes, chlorogenic acid was the most represented compound, reaching concentrations between 786.031 and 3735.790 µg g−1 (dw). Leaves accumulated high levels of chlorogenic acid in all genotypes, especially in “Spagnolo”, though the highest concentration for “Gigante” and “Bianco Avorio” was detected in hypocotyls (Supplementary Table S3). Phenylpropanoids represent essential constitutive and inducible defenses and are crucial for the safe growth and development of the plant. Therefore, these metabolites are usually allocated to protect tissues more exposed to biotic and abiotic stressors [50]. Consistently, the accumulation of total phenols was found to be high in leaves. The high accumulation of CQA precursors and other phenylpropanoids detected in hypocotyls and, to a lower extent, in cotyledons may be associated with the physiological function of these tissues. In particular, hypocotyls could be actively involved in metabolites transport to new growing tissues, actively translocating phenolics to new leaves, whereas the cotyledons could go through metabolic events to mobilize storage reserves [51]. In our tissues and genotypes, we did not detect di-caffeoyl quinic acids. Therefore, in our conditions di-caffeoyl quinic acid isomers, if any, were below the detection limit. We hypothesize that the accumulation of these compounds might depend on the physiological state of the plant and be developmentally regulated, as already reported [12]. These interesting findings remain to be further elucidated and will be the subject of future investigations.

3.2. Qualitative and Quantitative Profile of (Seeds, Hypocotyls, Cotyledons, and Leaves) Oil

Regarding the oil content, the results showed that the highest oil concentration is reached in the seeds, with an average level of 25.2%, in agreement with previous reports [52] but without substantial differences between the three analyzed genotypes (Figure 1). The lowest average level was found in the leaves for all genotypes, with the exception of “Gigante”, for which the lowest average oil content was detected in cotyledons.
The fatty acid profile of the oil represents the main factor that determines the choice of using the oil either for nutritional or industrial purposes. [53] found that cardoon oil is a rich source of unsaturated fatty acids such as linoleic and oleic acids (44.5% and 42.6%, respectively), whereas saturated fatty acids such as palmitic and stearic acid were detected in lower amounts (9.8% and 3.1%, respectively). Since the oil composition is deeply influenced by genotype, climatic factors, and geographical origin [54], growing the three cardoon genotypes in the same controlled greenhouse conditions should have allowed to highlight only differences due to genetic variability. Our results are in agreement with those reported by [53,55,56], highlighting a considerable variation in fatty acid profile between the analyzed genotypes, as shown in Table 3.
Palmitic acid was the most representative fatty acid in all the analyzed tissues and genotypes (Table 3 and Supplementary Table S4). In contrast to what was reported in the literature [57], in all the studied tissues and genotypes we did not detect stearic acid. Oleic acid content was found to vary from 6.54% to 19.27% and reached the highest amount in the oil extracted from the seeds of the genotype “Gigante”, but also in the other genotypes its content was higher in seeds than in all the other analyzed tissues (Table 3 and Supplementary Table S4). The content of linoleic acid varied from 9.20% to 45.43%, with hypocotyls of “Bianco Avorio” showing the highest concentration. Although the hypocotyls of the three genotypes accumulated the highest content of linoleic acid, it is worth noting that a considerable amount of linoleic acid was detected also in “Spagnolo” leaves (Table 3 and Supplementary Table S4). Finally, the content of palmitic acid varied from 34.33 to 63.32% and the highest amount was observed for the oil extracted from the cotyledons of the genotype “Gigante”. The fatty acids profile differences had a significant impact on the ratio of polyunsaturated/saturated fatty acids, ranging from 0.38 to 1.42, and highlighted the highest ratio for the oil extracted from hypocotyls of all genotypes analyzed. Surprisingly, the same was found for the oil extracted from the leaves of the genotype “Spagnolo”. This value is an important parameter in determining the nutritional value and functional properties of food products and oils as oils with a polyunsaturated/saturated fatty acids ratio greater than one are considered valuable edible oils.

3.3. Transcriptional Analysis of Key Biosynthetic Genes in Chlorogenic Acid and Monounsaturated Fatty Acids

To gain insight into the metabolic activities of cardoon tissues and organs, we analyzed transcripts abundance for HQT (hydroxycinnamoyl quinate transferase; Figure 2), as a key gene in the biosynthesis of chlorogenic acid, and for SAD (stearoyl acid desaturase) and FAD2 (fatty acid desaturase) genes (Figure 3), involved in oleic and linoleic acid formation, respectively. We identified the genes coding for enzymes involved in the biosynthesis of CGAs and mono-unsaturated fatty acids from deposited sequences. HQT was previously isolated and characterized in C. cardunculus. var altilis and globe artichoke [23,24], whereas for putative SAD and FAD2 genes, sequences were retrieved from BLAST searches in the C. cardunculus database (PRJNA453787) and chosen based on the high nucleotidic and proteic identity to orthologous genes in closely related Asteraceae and other oleaginous species. Regarding SAD, we found only two cardoon sequences, highly similar to each other, and coding for proteins closely related to characterized Lactuca sativa (XM 023873428.1) and Olea europaea var. sylvestris (XM 023023317.1) SAD2. As regards FAD genes, two main oleate desaturase genes occur in oleaginous plants, namely the microsomal and plastidial types [32]. In detail, two microsomial genes are characterized by tissue specificity, the seed type FAD2.1 and the house-keeping FAD2.2, along with the plastidial one, known as FAD6. We selected three FAD2 representatives: CcFAD2.1 (70.08% to GmFAD2-l, AAB00859, and 77.89% to HaFAD2-1, AAL68981.1), CcFAD2.2 (90.00% to CpFAD2-2, CAA76157 and 90.82% to HaFad2.2, AAL68982.1), and FAD6 (92.88% to AaFAD6, PWA37615 and 83% to TcFAD6, GEU42198). Although gene expression analyses were carried out for only one gene within the phenylpropanoid pathway and for SAD and three FAD2 isoforms for fatty acids biosynthesis, the selected genes represent key steps in the respective pathways. Gene transcripts were detected for all the tissues analyzed, and transcript abundances varied in tissues and genotypes underlining a tight regulation for the accumulation of these metabolites.

3.3.1. HQT Expression Analysis

In details, transcriptional analysis of HQT gene expression (Figure 2) was consistent with CQA accumulation among the genotypes (Table 2), which was higher in “Bianco Avorio” than “Gigante” and “Spagnolo”. Nevertheless, within each genotype, HQT expression did not reflect CQA accumulation for most of the analyzed tissues, in particular seeds and cotyledons. In all genotypes, HQT expression was barely detectable in seeds and much higher in cotyledons, while CGA levels were comparable to other tissues in seeds and the lowest in cotyledons (Figure 2 and Table 2). Due to their quiescent state, mature seeds might have a reduced transcriptional activity not temporally related with CGA accumulation, whose biosynthesis is highly regulated during seed development in order to guarantee a safe storage of energy in form of chlorogenic acid at seed maturity. Instead, in cotyledons, the high HQT transcriptional activity might be needed to provide high levels of CGA, which is degraded to provide lignin-precursors necessary for cotyledonary cell wall formation [58]. Regarding HQT expression in the other tissues we cannot exclude that additional pathways can contribute to CGA synthesis, as reported in other species [59]. Moreover, expression analysis can be only indicative of the transcriptional dynamics linked to CQA accumulation, which is mostly dependent on the enzymatic activity, as reported in other species [59,60]. We believe that HQT expression is also the result of a fine regulatory scheme, since this key gene is rate limiting and we cannot exclude that a feedback regulation is involved in its transcriptional rate, as recently reported [61].

3.3.2. SAD and FAD2 Expression Analysis

Transcriptional analysis of key biosynthetic genes in oleic and linoleic acid formation revealed that SAD expression was low in seeds of all three genotypes. In “Gigante” and “Spagnolo” it was low also in the other analyzed tissues, whereas in “Bianco Avorio” SAD transcription was remarkable in cotyledons but not consistent with quantitative oleic acid detection (Figure 3). The weak correlation between oleic acid accumulation and SAD expression might be related to the developmental seed stage investigated in this study. In mature seed, oleic acid might have reached an accumulation threshold, driving a feedback inhibition of SAD expression in favor of oleic to linoleic conversion. This hypothesis is consistent with the higher content of linoleic acid in almost all the analyzed tissues and genotypes. Regarding the expression analysis of the FAD isogenes, we found that the expression levels of the three FAD genes were higher in hypocotyls than in the other tissues in the “Gigante” genotype. Consistently, high linoleic acid accumulation was detected in hypocotyls of this genotype (Table 3), possibly resulting from the contribution of all the FAD isoforms. In “Spagnolo”, FAD2.1 and FAD2.2 showed a clear tissue-specificity, being highly expressed in hypocotyls and leaves, respectively, and we might hypothesize that each isoform contributed to the higher linoleic acid content of these two tissues in this genotype. As for FAD6, transcripts accumulation was higher in “Spagnolo” seeds and hypocotyls than in cotyledons and leaves. Regarding “Bianco Avorio”, high FAD2.1 expression was detected in cotyledons, where also FAD2.2 accumulated at high levels, similarly to leaves, whereas FAD6 was higher in seeds and hypocotyls. The distribution of FAD transcripts, in particular for FAD2.2, was in accordance with published literature in other oleaginous plants, indicating it as constitutive gene [59]. Indeed, FAD2.2 was expressed in all tissues, except for dormant seeds, and, according to transcriptional data, could play a main role in linoleic acid formation in “Spagnolo” and “Bianco Avorio”. Nevertheless, functional studies are necessary to ascertain the role of each gene in this species and whether it might be affected by genotype differences. Overall, our results reflect the harmonic interplay between SAD and FAD genes, the lower level of SAD being seemingly compensated by the higher expression of FAD genes, which finally regulates the MUFA/ PUFA ratio in all the cardoon tissues.

3.4. Antioxidant Activity of Polyphenolic Extracts

The results obtained for antioxidant activity of cultivated cardoon samples by the ABTS, FRAP and DPPH assays are presented in Table 4 and expressed as TEAC (mmol kg-1 dw). The measurement of the antioxidant activity of food is carried out using several methods since it is influenced by various factors. The most commonly used for the evaluation of the antioxidant potential are DPPH radical scavenging, ABTS decolorization and FRAP assays due to their simplicity, stability, accuracy and reproducibility [62].
As for seeds, FRAP values ranged from 49.83 ± 1.28 to 107.73 ± 1.30 mmol kg−1 dw, whereas ABTS values ranged from 87.24 ± 53.31 to 148.04 ± 1.13 mmol kg-1 dw and DPPH values ranged from 49.49 ± 0.27 to 143.87 ± 0.94 mmol kg-1 dw. Antioxidant activity showed significant differences between the tested genotypes, with “Bianco Avorio” having the overall highest antioxidant activity. Our results are in accordance with the results obtained by [11], reporting significant differences between wild and cultivated cardoon genotypes, as well as between hydromethanolic extracts and its residues. [63], similarly to our data, estimated that seed extracts of wild cardoons showed a high DPPH scavenging activity and a considerable ABTS•+ radical scavenging capacity. Moreover, the correlation coefficient between phenolic contents and TEAC was highly significant especially in the case of ABTS and FRAP (r2 = 0.82 and 0.94, respectively), indicating that polyphenolics may play an important role in free radical cations scavenging and ferric reducing antioxidant power.
Similar results were obtained for hypocotyls, where the correlation coefficients between phenolic contents and TEAC values were highly significant (r2 = 1, 0.94, and 0.95, respectively for DPPH, ABTS and FRAP). The results obtained showed a remarkable antioxidant activity of the genotype “Bianco Avorio” also for the hypocotyls and the highest antioxidant activity of polyphenols especially as free radical cation scavengers, being the antioxidant activity evaluated with the ABTS assay the highest. For cardoon hypocotyls, in fact, FRAP values ranged from 57.94 ± 1.043 to 79.76.73 ± 2.69 mmol kg−1 dw, ABTS values ranged from 317.21 ± 2.071 to 667.37 ± 1. mmol kg−1 dw and finally DPPH values ranged from 37.54 ± 0.311 to 60.87 ± 0.106 mmol kg-1 dw. These results highlight a significant variability in nutraceutical potential between genotypes.
For cotyledons, the results showed a significant variability in antioxidant potential between genotypes (Table 4). FRAP values ranged from 47.78 ± 0.987 to 79.148 ± 2.876 mmol kg−1 dw, whereas ABTS values ranged from 412.113 ± 0.112 to 543.012 ± 0.765 mmol kg-1 dw and finally DPPH values ranged from 19.192 ± 0.123 to 63.951 ± 0.543 mmol kg−1 dw. The highest antioxidant capacities were recorded for the “Spagnolo” genotype using the DPPH and FRAP methods while, in analogy with the results described so far, “Bianco Avorio” showed the highest scavenger capacity against the ABTS radical cations. In literature, the deposition of phenolic polymers was reported for Coffea arabica cotyledons during development, suggesting that lignification may occur via the utilization of the stored 5-O-caffeoylquinic acid. This result would explain the lower antioxidant activity of the cotyledons also associated with a low level of chlorogenic acid. Finally, for cardoon leaves, FRAP values ranged from 62.45 ± 0.432 to 91.675 ± 0.761 mmol kg−1 dw, whereas ABTS values ranged from 510.315 ± 0.766 to 771.412 ± 0.665 mmol kg−1 dw and finally DPPH values ranged from 27.167 ± 0.551 to 91.856 ± 0.116 mmol kg−1 dw.
The nutraceutical properties of leaves of artichoke and cardoon are linked to their specific chemical composition. Caffeoylquinic acids, followed by flavonoids, are the main phenolic compounds in leaves, among which chlorogenic acid is the most abundant. The same order of activity was found in the DPPH and FRAP assays, but the absolute values were lower than those found in the ABTS assay. According to literature data [64] ABTS assay is the most sensitive assay, as the working solution is soluble in aqueous and organic solvents, at several pH values, and the quenching reaction of free radicals is faster than other methods. The antioxidant capacity of cardoon leaf extracts was reported [12] and different artichoke genotypes were studied in terms of their antioxidant capacity [65]. Antioxidant activity measured on the leaf polyphenolic extracts was remarkably higher than that reported in literature, on the other hand little change in experimental conditions can widely influence the results making the comparison only qualitative. For leaves, the correlation coefficients between phenolic contents and TEAC values were mildly significant (r2 = 0.98, 0.83 and 0.70, respectively for DPPH, ABTS and FRAP), with the genotype “Bianco Avorio” showing the highest antioxidant activity especially as ABTS radical scavenger.

4. Conclusions

In conclusion, the results presented in this work provide a detailed analysis of tissues composition of cardoon, enabling to elucidate value-added products accumulation and distribution during plant development and hence contributing to optimize the sustainable use of this natural resource. Moreover, transcriptional and metabolic screening could be useful for a careful selection of the genotype.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-3921/9/11/1096/s1.

Author Contributions

Conceptualization, G.G., T.D.; methodology, G.G., T.D., M.D.P. and L.I.; software, G.G., T.D.; validation, G.G., T.D. and M.D.P.; formal analysis, T.D., G.G. and M.D.P.; investigation, G.G., T.D., L.I. and M.D.P.; resources, M.T. and A.R.; data curation, G.G. and T.D.; writing—original draft preparation, G.G., T.D., M.D.P., M.T. and R.T.; writing—review and editing, G.G., T.D., M.D.P., M.T., A.R.; supervision, A.R. and M.T.; project administration, A.R., F.S. and M.T.; funding acquisition, A.R., F.S. and M.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the following project: BOBCAT Biotecnologie per la produzione sostenibile di bio-materiali e prodotti specialistici utilizzando colture cellulari di cardo come bioraffineria. Fondazione Cariplo (2018_0955). In addition, part of this work was also carried out with the support of the CNR project NUTR-AGE (FOE-2019, DSB.AD004.271, Subtask 2.1.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bianco, V.V. Present situation and future potential of artichoke in the mediterranean basin. Acta Hortic. 2005, 681, 39–58. [Google Scholar] [CrossRef]
  2. Pappalardo, H.D.; Toscano, V.; Puglia, G.D.; Genovese, C.; Raccuia, S.A. Cynara cardunculus L. as a Multipurpose Crop for Plant Secondary Metabolites Production in Marginal Stressed Lands. Front. Plant Sci. 2020, 11, 240. [Google Scholar] [CrossRef] [PubMed]
  3. Ciancolini, A.; Alignan, M.; Pagnotta, M.A.; Vilarem, G.; Crinò, P. Selection of Italian cardoon genotypes as industrial crop for biomass and polyphenol production. Ind. Crop. Prod. 2013, 51, 145–151. [Google Scholar] [CrossRef]
  4. Foti, S.; Mauromicale, G.; Raccuia, S.A.; Fallico, B.; Fanella, F. Possible alternative utilization of Cynara spp. I. Biomass, grain yield and chemical composition of grain. Ind. Crops Prod. 1999, 10, 219–228. [Google Scholar] [CrossRef]
  5. Petropoulos, S.; Fernandes, Â.; Pereira, C.; Tzortzakis, N.; Vaz, J.; Soković, M.; Barros, L.; Ferreira, I.C. Bioactivities, chemical composition and nutritional value of Cynara cardunculus L. seeds. Food Chem. 2019, 289, 404–412. [Google Scholar] [CrossRef] [PubMed]
  6. Ierna, A.; Mauro, R.P.; Mauromicale, G.; Catania, S.; Lancia, S.V.; Industriale, Z.; Palma, B.; Catania, I. Biomass, grain and energy yield in Cynara cardunculus L. as affected by fertilization, genotype and harvest time. Biomass Bioenergy 2011, 36, 404–410. [Google Scholar] [CrossRef]
  7. Falleh, H.; Ksouri, R.; Chaieb, K.; Karray-Bouraoui, N.; Trabelsi, N.; Boulaaba, M.; Abdelly, C. Phenolic composition of Cynara cardunculus L. organs, and their biological activities. Comptes Rendus Biol. 2008, 331, 372–379. [Google Scholar] [CrossRef]
  8. Calabria, L.M.; Emerenciano, V.P.; Ferreira, M.J.P.; Scotti, M.T.; Mabry, T.J. A Phylogenetic Analysis of Tribes of the Asteraceae Based on Phytochemical Data. Nat. Prod. Commun. 2007, 2, 277–285. [Google Scholar] [CrossRef]
  9. Lattanzio, V.; Cardinali, A.; Di Venere, D.; Linsalata, V.; Palmieri, S. Browning phenomena in stored artichoke (Cynara scolymus L.) heads: Enzymic or chemical reactions? Food Chem. 1994, 50, 1–7. [Google Scholar] [CrossRef]
  10. Williamson, G.; Holst, B. Dietary reference intake (DRI) value for dietary polyphenols: Are we heading in the right direction? Br. J. Nutr. 2008, 99, S55–S58. [Google Scholar] [CrossRef] [Green Version]
  11. Durazzo, A.; Foddai, M.; Temperini, A.; Azzini, E.; Venneria, E.; Lucarini, M.; Finotti, E.; Maiani, G.; Crinò, P.; Saccardo, F.; et al. Antioxidant Properties of Seeds from Lines of Artichoke, Cultivated Cardoon and Wild Cardoon. Antioxidants 2013, 2, 52–61. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Gouveia, S.C.; Castilho, P.C. Phenolic composition and antioxidant capacity of cultivated artichoke, Madeira cardoon and artichoke-based dietary supplements. Food Res. Int. 2012, 48, 712–724. [Google Scholar] [CrossRef]
  13. Mandim, F.; Petropoulos, S.A.; Dias, M.I.; Pinela, J.; Kostic, M.; Soković, M.; Santos-Buelga, C.; Ferreira, I.C.F.R.; Barros, L. Seasonal variation in bioactive properties and phenolic composition of cardoon (Cynara cardunculus var. altilis) bracts. Food Chem. 2020, 336, 127744. [Google Scholar] [CrossRef]
  14. Racchi, M.; Daglia, M.; Lanni, C.; Papetti, A.; Govoni, S.; Gazzani, G. Antiradical Activity of Water Soluble Components in Common Diet Vegetables. J. Agric. Food Chem. 2002, 50, 1272–1277. [Google Scholar] [CrossRef]
  15. Miller, N.J.; Rice-Evans, C.A. Spectrophotometric determination of antioxidant activity. Redox Rep. 1996, 2, 161–171. [Google Scholar] [CrossRef]
  16. Cho, A.S.; Jeon, S.M.; Kim, M.J.; Yeo, J.; Seo, K.I.; Choi, M.S.; Lee, M.K. Chlorogenic acid exhibits anti-obesity property and improves lipid metabolism in high-fat diet-induced-obese mice. Food Chem. Toxicol. 2010, 48, 937–943. [Google Scholar] [CrossRef]
  17. Dos Santos, M.D.; Almeida, M.C.; Lopes, N.P.; de Souza, G.E.P. Evaluation of the anti-inflammatory, analgesic and antipyretic activities of the natural polyphenol chlorogenic acid. Biol. Pharm. Bull. 2006, 29, 2236–2240. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Naveed, M.; Hejazi, V.; Abbas, M.; Kamboh, A.A.; Khan, G.J.; Shumzaid, M.; Ahmad, F.; Babazadeh, D.; FangFang, X.; Modarresi-Ghazani, F.; et al. Chlorogenic acid (CGA): A pharmacological review and call for further research. Biomed. Pharmacother. 2018, 97, 67–74. [Google Scholar] [CrossRef]
  19. Santana-Gálvez, J.; Cisneros-Zevallos, L.; Jacobo-Velázquez, D. Chlorogenic Acid: Recent Advances on Its Dual Role as a Food Additive and a Nutraceutical against Metabolic Syndrome. Molecules 2017, 22, 358. [Google Scholar] [CrossRef] [Green Version]
  20. Lou, Z.; Wang, H.; Zhu, S.; Ma, C.; Wang, Z. Antibacterial Activity and Mechanism of Action of Chlorogenic Acid. J. Food Sci. 2011, 76, M398–M403. [Google Scholar] [CrossRef]
  21. Niggeweg, R.; Michael, A.J.; Martin, C. Engineering plants with increased levels of the antioxidant chlorogenic acid. Nat. Biotechnol. 2004, 22, 746–754. [Google Scholar] [CrossRef]
  22. Mahesh, V.; Million-Rousseau, R.; Ullmann, P.; Chabrillange, N.; Bustamante, J.; Mondolot, L.; Morant, M.; Noirot, M.; Hamon, S.; de Kochko, A.; et al. Functional characterization of two p-coumaroyl ester 3′-hydroxylase genes from coffee tree: Evidence of a candidate for chlorogenic acid biosynthesis. Plant Mol. Biol. 2007, 64, 145–159. [Google Scholar] [CrossRef]
  23. Sonnante, G.; D’Amore, R.; Blanco, E.; Pierri, C.L.; De Palma, M.; Luo, J.; Tucci, M.; Martin, C. Novel hydroxycinnamoyl-coenzyme A quinate transferase genes from artichoke are involved in the synthesis of chlorogenic acid. Plant Physiol. 2010, 153, 1224–1238. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Comino, C.; Hehn, A.; Moglia, A.; Menin, B.; Bourgaud, F.; Lanteri, S.; Portis, E. The isolation and mapping of a novel hydroxycinnamoyltransferase in the globe artichoke chlorogenic acid pathway. BMC Plant Biol. 2009, 9, 30. [Google Scholar] [CrossRef] [Green Version]
  25. Cheevarungnapakul, K.; Khaksar, G.; Panpetch, P.; Boonjing, P.; Sirikantaramas, S. Identification and Functional Characterization of Genes Involved in the Biosynthesis of Caffeoylquinic Acids in Sunflower (Helianthus annuus L.). Front. Plant Sci. 2019, 10, 968. [Google Scholar] [CrossRef] [Green Version]
  26. Raccuia, S.A.; Melilli, M.G. Biomass and grain oil yields in Cynara cardunculus L. genotypes grown in a Mediterranean environment. Field Crop. Res. 2007, 101, 187–197. [Google Scholar] [CrossRef]
  27. Gominho, J.; Lourenc, A.; Palma, P.; Lourenc, M.E.; Curt, M.D.; Fernández, J.; Pereira, H. Large scale cultivation of Cynara cardunculus L. for biomass production—A case study. Ind. Crops Prod. 2011, 33, 1–6. [Google Scholar] [CrossRef]
  28. Meï, C.; Michaud, M.; Cussac, M.; Albrieux, C.; Gros, V.; Maréchal, E.; Block, M.A.; Jouhet, J.; Rébeillé, F. Levels of polyunsaturated fatty acids correlate with growth rate in plant cell cultures. Sci. Rep. 2015, 5, 1–9. [Google Scholar] [CrossRef] [Green Version]
  29. Upchurch, R.G. Fatty acid unsaturation, mobilization, and regulation in the response of plants to stress. Biotechnol. Lett. 2008, 30, 967–977. [Google Scholar] [CrossRef]
  30. Somerville, C.; Browse, J. Plant Lipids: Metabolism, Mutants, and Membranes. Science 1991, 252, 80–87. [Google Scholar] [CrossRef]
  31. Browse, J.; Somerville, C. Glycerolipid Synthesis: Biochemistry and Regulation. Annu. Rev. Plant Biol. 1991, 42, 467–506. [Google Scholar] [CrossRef]
  32. Dar, A.A.; Choudhury, A.R.; Kancharla, P.K.; Arumugam, N. The FAD2 Gene in Plants: Occurrence, Regulation, and Role. Front. Plant Sci. 2017, 8, 1–16. [Google Scholar] [CrossRef] [Green Version]
  33. Kumar, A.; Sharma, A.; Upadhyaya, K.C. Vegetable Oil: Nutritional and Industrial Perspective. Curr. Genom. 2016, 17, 230–240. [Google Scholar] [CrossRef] [Green Version]
  34. Orsavova, J.; Misurcova, L.; Ambrozova, J.; Vicha, R.; Mlcek, J. Fatty Acids Composition of Vegetable Oils and Its Contribution to Dietary Energy Intake and Dependence of Cardiovascular Mortality on Dietary Intake of Fatty Acids. Int. J. Mol. Sci. 2015, 16, 12871–12890. [Google Scholar] [CrossRef] [PubMed]
  35. Clé, C.; Hill, L.M.; Niggeweg, R.; Martin, C.R.; Guisez, Y.; Prinsen, E.; Jansen, M.A.K.K. Modulation of chlorogenic acid biosynthesis in Solanum lycopersicum; consequences for phenolic accumulation and UV-tolerance. Phytochemistry 2008, 69, 2149–2156. [Google Scholar] [CrossRef]
  36. Pfaffl, M.W. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29, e45. [Google Scholar] [CrossRef]
  37. Pfaffl, M. Quantification strategies in real-time PCR. AZ Quant. PCR 2004, 1, 87–112. [Google Scholar]
  38. Docimo, T.; De Stefano, R.; Cappetta, E.; Lisa Piccinelli, A.; Celano, R.; De Palma, M.; Tucci, M. Physiological, biochemical, and metabolic responses to short and prolonged saline stress in two cultivated cardoon genotypes. Plants 2020, 9, 554. [Google Scholar] [CrossRef]
  39. Tungmunnithum, D.; Garros, L.; Drouet, S.; Renouard, S.; Lainé, E.; Hano, C. Green Ultrasound Assisted Extraction of trans Rosmarinic Acid from Plectranthus scutellarioides (L.) R.Br. Leaves. Plants 2019, 8, 50. [Google Scholar] [CrossRef] [Green Version]
  40. Re, R.; Pellegrini, N.; Proteggente, A.; Pannala, A.; Yang, M.; Rice-Evans, C. Antioxidant activity applying an improved ABTS radical cation decolorization assay. Free Radic. Biol. Med. 1999, 26, 1231–1237. [Google Scholar] [CrossRef]
  41. Brand-Williams, W.; Cuvelier, M.E.; Berset, C. Use of a free radical method to evaluate antioxidant activity. LWT-Food Sci. Technol. 1995, 28, 25–30. [Google Scholar] [CrossRef]
  42. Benzie, I.F.F.; Strain, J.J. The ferric reducing ability of plasma as a measure of antioxodant. Anal. Biochem. 1996, 239, 70–76. [Google Scholar] [CrossRef] [Green Version]
  43. Graziani, G.; Ritieni, A.; Cirillo, A.; Cice, D.; Di Vaio, C. Effects of Biostimulants on Annurca Fruit Quality and Potential Nutraceutical Compounds at Harvest and during Storage. Plants 2020, 9, 775. [Google Scholar] [CrossRef]
  44. Kostić, M.D.; Joković, N.M.; Stamenković, O.S.; Rajković, K.M.; Milić, P.S.; Veljković, V.B. Optimization of hempseed oil extraction by n-hexane. Ind. Crop. Prod. 2013, 48, 133–143. [Google Scholar] [CrossRef]
  45. Pandino, G.; Lombardo, S.; Mauromicale, G.; Williamson, G. Phenolic acids and flavonoids in leaf and floral stem of cultivated and wild Cynara cardunculus L. genotypes. Food Chem. 2011, 126, 417–422. [Google Scholar] [CrossRef]
  46. Dias, M.I.; Barros, L.; Barreira, J.C.M.; Alves, M.J.; Barracosa, P.; Ferreira, I.C.F.R. Phenolic profile and bioactivity of cardoon (Cynara cardunculus L.) inflorescence parts: Selecting the best genotype for food applications. Food Chem. 2018, 268, 196–202. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Ramos, P.A.B.; Santos, S.A.O.; Guerra, Â.R.; Guerreiro, O.; Freire, C.S.R.; Rocha, S.M.; Duarte, M.F.; Silvestre, A.J.D. Phenolic composition and antioxidant activity of different morphological parts of Cynara cardunculus L. var. altilis (DC). Ind. Crop. Prod. 2014, 61, 460–471. [Google Scholar] [CrossRef]
  48. Lattanzio, V.; Kroon, P.A.; Linsalata, V.; Cardinali, A. Globe artichoke: A functional food and source of nutraceutical ingredients. J. Funct. Foods 2009, 1, 131–144. [Google Scholar] [CrossRef]
  49. Clifford, M.N.; Johnston, K.L.; Knight, S.; Kuhnert, N. Hierarchical scheme for LC-MSn identification of chlorogenic acids. J. Agric. Food Chem. 2003, 51, 2900–2911. [Google Scholar] [CrossRef] [PubMed]
  50. Tattini, M.; Loreto, F.; Fini, A.; Guidi, L.; Brunetti, C.; Velikova, V.; Gori, A.; Ferrini, F. Isoprenoids and phenylpropanoids are part of the antioxidant defense orchestrated daily by drought-stressed Platanus × acerifolia plants during Mediterranean summers. New Phytol. 2015, 207, 613–626. [Google Scholar] [CrossRef]
  51. Penfield, S.; Graham, S.; Graham, I.A. Storage reserve mobilization in germinating oilseeds: Arabidopsis as a model system. Biochem. Soc. Trans. 2005, 33, 380–383. [Google Scholar] [CrossRef]
  52. Gominho, J.; Curt, M.D.; Lourenço, A.; Fernández, J.; Pereira, H. Cynara cardunculus L. as a biomass and multi-purpose crop: A review of 30 years of research. Biomass Bioenergy 2018, 109, 257–275. [Google Scholar] [CrossRef]
  53. Raccuia, S.A.; Piscioneri, I.; Sharma, N.; Melilli, M.G. Genetic variability in Cynara cardunculus L. domestic and wild types for grain oil production and fatty acids composition. Biomass Bioenergy 2011, 35, 3167–3173. [Google Scholar] [CrossRef]
  54. Francisco, M.; Elena Cartea, M.; Maria Butron, A.; Sotelo, T.; Velasco, P. Environmental and Genetic Effects on Yield and Secondary Metabolite Production in Brassica rapa Crops. J. Agric. Food Chem. 2012, 60, 5507–5514. [Google Scholar] [CrossRef]
  55. Maccarone, E.; Fallico, B.; Fanella, F.; Mauromicale, G.; Raccuia, S.A. Possible alternative utilization of Cynara spp. II. Chemical characterization of their grain oil. Ind. Crops Prod. 1999, 10, 229–237. [Google Scholar] [CrossRef]
  56. Curt, M.D.; Sánchez, G.; Fernández, J. The potential of Cynara cardunculus L. for seed oil production in a perennial cultivation system. Biomass Bioenergy 2002, 23, 33–46. [Google Scholar] [CrossRef]
  57. Angelova, V.; Nemska, M.P.; Krustev, L. Chemical composition of cardoon (Cynara cardunculus L.) Grown in south bulgaria. AGROFOR Int. J. 2019, 4, 100–110. [Google Scholar]
  58. Volpi E Silva, N.; Mazzafera, P.; Cesarino, I. Should I stay or should I go: Are chlorogenic acids mobilized towards lignin biosynthesis? Phytochemistry 2019, 166, 112063. [Google Scholar] [CrossRef]
  59. Payyavula, R.S.; Navarre, D.; Kuhl, J.C.; Pantoja, A.; Pillai, S.S. Differential effects of environment on potato phenylpropanoid and carotenoid expression. BMC Plant Biol. 2012, 12, 39. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Payyavula, R.S.; Shakya, R.; Sengoda, V.G.; Munyaneza, J.E.; Swamy, P.; Navarre, D.A. Synthesis and regulation of chlorogenic acid in potato: Rerouting phenylpropanoid flux in HQT-silenced lines. Plant Biotechnol. J. 2015, 13, 551–564. [Google Scholar] [CrossRef]
  61. Li, Y.; Kong, D.; Bai, M.; He, H.; Wang, H.; Wu, H. Correlation of the temporal and spatial expression patterns of HQT with the biosynthesis and accumulation of chlorogenic acid in Lonicera japonica flowers. Hortic. Res. 2019, 6, 73. [Google Scholar] [CrossRef] [Green Version]
  62. Vijaya Kumar Reddy, C.; Sreeramulu, D.; Raghunath, M. Antioxidant activity of fresh and dry fruits commonly consumed in India. Food Res. Int. 2010, 43, 285–288. [Google Scholar] [CrossRef]
  63. Kollia, E.; Markaki, P.; Zoumpoulakis, P.; Proestos, C. Antioxidant activity of Cynara scolymus L. and Cynara cardunculus L. extracts obtained by different extraction techniques. Nat. Prod. Res. 2017, 31, 1163–1167. [Google Scholar] [CrossRef]
  64. Arnao, M.B. Some methodological problems in the determination of antioxidant activity using chromogen radicals: A practical case. Trends Food Sci. Technol. 2000, 11, 419–421. [Google Scholar] [CrossRef]
  65. Colla, G.; Rouphael, Y.; Cardarelli, M.; Svecova, E.; Rea, E.; Lucini, L. Effects of saline stress on mineral composition, phenolic acids and flavonoids in leaves of artichoke and cardoon genotypes grown in floating system. J. Sci. Food Agric. 2013, 93, 1119–1127. [Google Scholar] [CrossRef]
Figure 1. Percent oil content (w:w) in the n-exan extract in cotyledons, leaves, hypocotyls, and seeds of cardoon genotypes. Each value represents the mean of three biological and two technical replicates. Different letters denote significant differences among tissues by analysis of variance [ANOVA]. Statistical significance was defined as p < 0.05, using Tukey’s post hoc test for mean separation.
Figure 1. Percent oil content (w:w) in the n-exan extract in cotyledons, leaves, hypocotyls, and seeds of cardoon genotypes. Each value represents the mean of three biological and two technical replicates. Different letters denote significant differences among tissues by analysis of variance [ANOVA]. Statistical significance was defined as p < 0.05, using Tukey’s post hoc test for mean separation.
Antioxidants 09 01096 g001
Figure 2. Transcriptional analysis of the CcHQT gene by qRT-PCR in seeds, hypocotyls, cotyledons and leaves of “Gigante”, “Spagnolo” and “Bianco Avorio” cultivated cardoon genotypes. Results are expressed as fold changes relatively to seeds, used as internal calibrator of gene expression for each genotype. Each value represents the mean ± SD of three biological and two technical replicates. Different letters denote a significant difference between tissues by analysis of variance [ANOVA]. Statistical significance was defined as p < 0.05, using the Tukey’s post hoc test for mean separation.
Figure 2. Transcriptional analysis of the CcHQT gene by qRT-PCR in seeds, hypocotyls, cotyledons and leaves of “Gigante”, “Spagnolo” and “Bianco Avorio” cultivated cardoon genotypes. Results are expressed as fold changes relatively to seeds, used as internal calibrator of gene expression for each genotype. Each value represents the mean ± SD of three biological and two technical replicates. Different letters denote a significant difference between tissues by analysis of variance [ANOVA]. Statistical significance was defined as p < 0.05, using the Tukey’s post hoc test for mean separation.
Antioxidants 09 01096 g002
Figure 3. Transcriptional profiling of SAD, FAD2.1, FAD2.2 and FAD6 genes by qRT-PCR in seeds, hypocotyls, cotyledons and leaves of cardoon genotypes. (a) “Gigante”, (b) “Spagnolo”, (c) “Bianco Avorio”. Results are expressed as fold changes relatively to the tissue with the lowest expression level, used as internal calibrator for each gene within each genotype. Each value represents the mean ± SD of three biological and two technical replicates. For each gene, different letters denote a significant difference between tissues within each genotype by analysis of variance [ANOVA]. Statistical significance was defined as p < 0.05, using the Tukey’s post hoc test for mean separation.
Figure 3. Transcriptional profiling of SAD, FAD2.1, FAD2.2 and FAD6 genes by qRT-PCR in seeds, hypocotyls, cotyledons and leaves of cardoon genotypes. (a) “Gigante”, (b) “Spagnolo”, (c) “Bianco Avorio”. Results are expressed as fold changes relatively to the tissue with the lowest expression level, used as internal calibrator for each gene within each genotype. Each value represents the mean ± SD of three biological and two technical replicates. For each gene, different letters denote a significant difference between tissues within each genotype by analysis of variance [ANOVA]. Statistical significance was defined as p < 0.05, using the Tukey’s post hoc test for mean separation.
Antioxidants 09 01096 g003
Table 1. Retention time and exact mass spectra data of cultivated cardoon polyphenols investigated by UHPLC-HRMS Orbitrap.
Table 1. Retention time and exact mass spectra data of cultivated cardoon polyphenols investigated by UHPLC-HRMS Orbitrap.
Peak NumberPolyphenolRetention Time (min) Chemical Formula Theoretical Mass (m/z) Measured Mass (m/z) Accuracy (ppm)
13Luteolin10.99C15H10O6285.04046285.040710.88
15Apigenin11.45C15H10O5269.04555269.046051.86
6Diosmin9.96C28H32O15607.16684607.168212.26
5Apigenin-8-C-glucoside9.89C21H20O10431.09837431.099733.15
11(+/−) Naringenin10.58C15H12O5271.06120271.061651.66
4Quercitin-3-O-glucoside9.71C21H19O12463.08820463.089262.29
12Quercetin10.71C15H10O7301.03538301.035400.07
14Kaempferol11.31C15H10O6285.04046285.040771.09
9Myricetin10.06C15H10O8317.03029317.031192.84
7Naringin10.01C27H32O14579.54123579.541620.67
10Kaempferol-3-O-glucoside10.11C21H20O11447.09328447.093961.52
8Luteolin-7-O-glucoside10.05C21H20O12447.09328447.09311−0.38
3p-Coumaric acid8.66C9H8O3163.04007163.040563.00
2Chlorogenic acid7.85C16H18O9353.08781353.088562.12
14-Hydroxybenzoic acid6.76C7H6O3137.02442137.024893.43
Table 2. Polyphenols content in seeds, hypocotyls, cotyledons, and leaves of “Gigante”, “Spagnolo” and “Bianco Avorio” genotypes detected by HRMS-Orbitrap. Values are expressed in mg g−1 (dw). Each value represents the mean of three biological and two technical replicates. Different letters denote a significant difference between genotypes within each tissue by analysis of variance [ANOVA]. Statistical significance was defined as p < 0.05, using Tukey’s post hoc test for mean separation. nd: not detected.
Table 2. Polyphenols content in seeds, hypocotyls, cotyledons, and leaves of “Gigante”, “Spagnolo” and “Bianco Avorio” genotypes detected by HRMS-Orbitrap. Values are expressed in mg g−1 (dw). Each value represents the mean of three biological and two technical replicates. Different letters denote a significant difference between genotypes within each tissue by analysis of variance [ANOVA]. Statistical significance was defined as p < 0.05, using Tukey’s post hoc test for mean separation. nd: not detected.
Phenolic CompoundsSEEDSHYPOCOTYLSCOTYLEDONSLEAVES
GiganteSpagnoloB. AvorioGiganteSpagnoloB. AvorioGiganteSpagnoloB. AvorioGiganteSpagnoloB. Avorio
4-hydroxy benzoic acid2.750a2.510b2.040c5.620a3.760c3.910b12.592b8.795c16.771a1.221a0.661b0.329c
Vitexin107.300a91.400b83.700c121.700b111.710c250.110a350.432c876.529b1182.892a528.396c2100.344b2107.906a
luteolin-7-O-glucoside1.470b2.760a0.910c1.610c8.610b12.700a18.092c34.299a24.086b81.922c122.687a88.773b
naringin10.900c11.700b14.400a5.710c51.700b67.100a14.599b438.433a436.022a0.251b0.357a0.038c
chlorogenic acid1036.310c1201.990b1430.460a3461.130b1006.780c3735.790a786.031a769.301a463.390b1468.968b2467.679a2637.733a
coumaric acid1.790a1.840a1.820a3.620a3.410a3.720a0.180a0.080b0.080b0.080c0.560a0.200b
quercetin-3-glucoside0.090c0.270b1.160a0.700b0.180c1.180a1.587b2.670a1.970b3.386a3.618a2.050b
diosmin0.050andndndndndndndndndndnd
kaempferol-3-O-glucoside2.180a2.140a2.200a4.910ab4.460b5.000a0.420b0.780a0.860a9.960a9.960a6.720b
myricetin0.440a0.440a0.450a0.880a0.890a0.880a2.470c3.110b8.470a6.884a6.063b5.493c
naringenin1.730a1.610b1.680ab3.130a3.110a3.110a0.120and0.040b0.160c0.320b0.500a
luteolin0.320a0.020b0.020b0.040b0.200and11.700c37.137a31.464b56.040c84.640a72.671b
kaempferol0.340a0.020b0.020b0.040b0.240and3.600c4.300b7.900a4.800c7.360a5.780b
quercetin0.820ndndnd1.640ndnd0.020a0.020a36.580a21.300b14.160c
apigenin0.240a0.040b0.020cnd0.080a0.040bnd4.580b7.040a2.260c3.710b8.080a
Total polyphenols1166.730c1316.740b1538.880a3609.090b1196.770c4083.540a1201.823b2180.034a2181.005a2200.908b4829.259a4950.433a
Table 3. GC/MS analysis of fatty acids composition in cardoon genotypes (expressed as % of the total fatty acid composition). Each value represents the mean of three biological and two technical replicates. Different letters denote a significant difference between genotypes within each tissue by analysis of variance [ANOVA]. Statistical significance was defined as p < 0.05, using Tukey’s post hoc test for mean separation. n.d.: not detected.
Table 3. GC/MS analysis of fatty acids composition in cardoon genotypes (expressed as % of the total fatty acid composition). Each value represents the mean of three biological and two technical replicates. Different letters denote a significant difference between genotypes within each tissue by analysis of variance [ANOVA]. Statistical significance was defined as p < 0.05, using Tukey’s post hoc test for mean separation. n.d.: not detected.
Fatty AcidsSeedsHypocotylsCotyledonsLeaves
GiganteSpagnoloB. AvorioGiganteSpagnoloB. AvorioGiganteSpagnoloB. AvorioGiganteSpagnoloB. Avorio
Pentadeconoic (C15:0)0.123c0.283a0.173b0.172b0.182a0.143c1.852a0.578c0.711b0.985a0.441c0.611b
Palmitic (C16:0) 59.868a44.265b56.708a49.361a47.722a40.791b63.320a44.594c50.348b44.484b34.330c61.643a
Margaric (C17:0) 0.694b1.25a0.797b0.109a0.117a0.098b2.253a0.639b0.642b0.882a0.571b0.608b
Nonadecanoic (C19:0) 0.013b0.024a0.01b4n.d.n.d.n.d.n.d.0.252n.d.0.289n.d.0.187
Arachidic (C20:0) 4.446ab4.427a3.4bn.d.n.d.0.659n.d.8.176a5.157b7.635a3.852c5.309b
Behenic (C22:0) 0.634a0.635a0.489b0.271b0.431a0.315b1.976b3.767a1.459c4.162a2.103c3.463b
Lignoceric (C24:0) 0.457a0.246b0.432a0.326n.d.n.d.2.508b5.408a0.979c5.219bn.d.5.706a
Cerotic (C26:0) 0.052n.d.n.d.n.d.n.d.n.d.0.465a4.414an.d.3.977n.d.n.d.
Melissic (C30:0) 0.0150.03n.d.n.d.n.d.0.282n.d.n.d.0.557n.d.n.d.n.d.
Palmitoleic (C16:1) 0.372b0.816a0.348bn.d.0.135n.d.0.607an.d.0.439b1.452a0.420c0.516b
Hexadicadienoic (C16:2)0.006b1.389a0.009bn.d.n.d.n.d.n.d.n.d.n.d.n.d.0.088n.d.
Hexadicatrienoic (C16:3)n.d.n.d.n.d.n.d.n.d.n.d.9.774n.d.n.d.n.d.n.d.n.d.
Oleic (C18:1)19.271a18.395b15.828c6.54c7.869b12.278an.d.9.737a8.240b10.524b15.519a6.592c
Linoleic (C18:2)13.684c27.34a21.212b37.892b38.954b45.434a8.935c12.997b24.714a12.806b37.460a9.197c
Linolenic (C18:3) 0.355a0.441a0.021b5.259a4.588bn.d.8.312ab9.438a6.644b7.584a5.215c6.169b
Nonadecenoic (C19:1) 0.010.018n.d.n.d.n.d.n.d.n.d.n.d.n.d.n.d.n.d.n.d.
Gadoleic (C20:1) n.d.0.441b0.566a0.07n.d.n.d.n.d.n.d.0.110n.d.n.d.n.d.
Table 4. Antioxidant activity in seeds, hypocotyls, cotyledons and leaves of three cultivated cardoon genotypes. Each value represents the mean ± SD of three biological and two technical replicates. Different letters denote a significant difference between genotypes within each tissue by analysis of variance [ANOVA] Statistical significance was defined as p < 0.05, using the Tukey’s post hoc test for mean separation.
Table 4. Antioxidant activity in seeds, hypocotyls, cotyledons and leaves of three cultivated cardoon genotypes. Each value represents the mean ± SD of three biological and two technical replicates. Different letters denote a significant difference between genotypes within each tissue by analysis of variance [ANOVA] Statistical significance was defined as p < 0.05, using the Tukey’s post hoc test for mean separation.
GENOTYPEDPPHABTSFRAP
mmol trolox kg−1mmol trolox kg−1mmol trolox kg−1
SeedsGIGANTE49.492 ± 0.275c87.237 ± 3.312c49.833 ± 1.283c
SPAGNOLO143.871 ± 0.942a135.413 ± 1.119b59.848 ± 3.340b
BIANCO AVORIO137.946 ± 0.265b148.037 ± 1.127a107.728 ± 1.303a
HypocotylsGIGANTE57.149 ± 0.377b667.365 ± 1.036a72.000 ± 0.954b
SPAGNOLO37.543 ± 0.311c317.206 ± 2.071c57.939 ± 1.043c
BIANCO AVORIO60.866 ± 0.106a631.482 ± 8.285b79.758 ± 2.694a
CotyledonsGIGANTE19.192 ± 0.123c487.277 ± 2.111b47.776 ± 0.987c
SPAGNOLO63.951 ± 0.543a412.113 ± 0.112c79.148 ± 2.876a
BIANCO AVORIO55.126 ± 0.156b543.012 ± 0.765a57.122 ± 1.832b
LeavesGIGANTE27.167 ± 0.551c510.315 ± 0.765c62.450 ± 0.432c
SPAGNOLO81.543 ± 0.111b654.216 ± 0.981b97.659 ± 0.876a
BIANCO AVORIO91.856 ± 0.116a771.412 ± 0.665a91.675 ± 0.761b
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Graziani, G.; Docimo, T.; Palma, M.D.; Sparvoli, F.; Izzo, L.; Tucci, M.; Ritieni, A. Changes in Phenolics and Fatty Acids Composition and Related Gene Expression during the Development from Seed to Leaves of Three Cultivated Cardoon Genotypes. Antioxidants 2020, 9, 1096. https://doi.org/10.3390/antiox9111096

AMA Style

Graziani G, Docimo T, Palma MD, Sparvoli F, Izzo L, Tucci M, Ritieni A. Changes in Phenolics and Fatty Acids Composition and Related Gene Expression during the Development from Seed to Leaves of Three Cultivated Cardoon Genotypes. Antioxidants. 2020; 9(11):1096. https://doi.org/10.3390/antiox9111096

Chicago/Turabian Style

Graziani, Giulia, Teresa Docimo, Monica De Palma, Francesca Sparvoli, Luana Izzo, Marina Tucci, and Alberto Ritieni. 2020. "Changes in Phenolics and Fatty Acids Composition and Related Gene Expression during the Development from Seed to Leaves of Three Cultivated Cardoon Genotypes" Antioxidants 9, no. 11: 1096. https://doi.org/10.3390/antiox9111096

APA Style

Graziani, G., Docimo, T., Palma, M. D., Sparvoli, F., Izzo, L., Tucci, M., & Ritieni, A. (2020). Changes in Phenolics and Fatty Acids Composition and Related Gene Expression during the Development from Seed to Leaves of Three Cultivated Cardoon Genotypes. Antioxidants, 9(11), 1096. https://doi.org/10.3390/antiox9111096

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop