Olive phenolic compounds: metabolic and transcriptional profiling during fruit development
© Alagna et al.; licensee BioMed Central Ltd. 2012
Received: 26 April 2012
Accepted: 30 August 2012
Published: 10 September 2012
Olive (Olea europaea L.) fruits contain numerous secondary metabolites, primarily phenolics, terpenes and sterols, some of which are particularly interesting for their nutraceutical properties. This study will attempt to provide further insight into the profile of olive phenolic compounds during fruit development and to identify the major genetic determinants of phenolic metabolism.
The concentration of the major phenolic compounds, such as oleuropein, demethyloleuropein, 3–4 DHPEA-EDA, ligstroside, tyrosol, hydroxytyrosol, verbascoside and lignans, were measured in the developing fruits of 12 olive cultivars. The content of these compounds varied significantly among the cultivars and decreased during fruit development and maturation, with some compounds showing specificity for certain cultivars. Thirty-five olive transcripts homologous to genes involved in the pathways of the main secondary metabolites were identified from the massive sequencing data of the olive fruit transcriptome or from cDNA-AFLP analysis. Their mRNA levels were determined using RT-qPCR analysis on fruits of high- and low-phenolic varieties (Coratina and Dolce d’Andria, respectively) during three different fruit developmental stages. A strong correlation was observed between phenolic compound concentrations and transcripts putatively involved in their biosynthesis, suggesting a transcriptional regulation of the corresponding pathways. OeDXS, OeGES, OeGE10H and OeADH, encoding putative 1-deoxy-D-xylulose-5-P synthase, geraniol synthase, geraniol 10-hydroxylase and arogenate dehydrogenase, respectively, were almost exclusively present at 45 days after flowering (DAF), suggesting that these compounds might play a key role in regulating secoiridoid accumulation during fruit development.
Metabolic and transcriptional profiling led to the identification of some major players putatively involved in biosynthesis of secondary compounds in the olive tree. Our data represent the first step towards the functional characterisation of important genes for the determination of olive fruit quality.
KeywordsOlea europaea Phenolics Secoiridoids RT-qPCR Transcriptome Secondary metabolism
The olive fruit mesocarp accumulates a wide range of secondary metabolites. The main category of secondary metabolites is represented by secoiridoids, a group of monoterpenoids with a cleaved methylcyclopentane skeleton, which are typical of the Oleaceae and other few dicotyledonous families. Secoiridoids are abundant in olives as phenol-conjugated compounds that might contain a glycoside moiety. The most important secoiridoids of the olive fruit and virgin olive oil are oleuropein, demethyloleuropein, oleuroside, ligstroside, nüzhenide and their aglycon forms, such as the dialdehydic form of decarboxymethyl elenolic acid linked to either 3,4-DHPEA or p-HPEA (3,4-DHPEA-EDA and p-HPEA-EDA, respectively), an isomer of oleuropein aglycon (3,4-DHPEA-EA) and the ligstroside aglycon (p-HPEA-EA) . The secoiridoid compound p-HPEA-EDA, also called oleocanthal, was identified in the extra-virgin olive oil but has never been observed in the fruits, probably as a consequence of post-harvest enzymatic activity . Other olive phenolics include phenolic acids, phenolic alcohols (e.g., hydroxytyrosol (3,4-DHPEA) and tyrosol (p-HPEA), flavonoids and lignans . These compounds are observed in all parts of the drupe, with the highest concentration in the pulp . The fruits of several olive cultivars also contain high amounts of verbascoside  and other phenolics present at maturity, such as homovanillic alcohol, 3,4-dihydroxyphenylacetic acid (DHPAC), caffeic acid, p-coumaric acid, phloretic acid, vanillic acid  and low amounts of comselogoside .
Studies of the phenolic profiles in the mesocarp, exocarp, seed, stone and leaf of the olive have demonstrated that different tissues contain distinct compounds. For example, nüzhenide and salidroside are only observed in the olive seed , while the flavonoids luteolin-7-glucoside, rutin and quercetin are exclusively present in the fruit peel .
Squalene, an intermediate of the sterol pathway, is the precursor of α- and β-amyrins and the triterpenic diols erythrodiol and uvaol. Squalene is another important compound with recognised effects on human health that is present in consistent amount only in olive and other vegetable oils . Sterols, such as β-sitosterol, cycloartenol and 24-methylenecycloartanol, progressively accumulate when the olive fruit reaches its final size and veraison begins .
Secoiridoids are not soluble in oil and, after the process of mechanical extraction, only a small portion is recovered in the oil, representing the most important microconstituents of virgin olive oils for their health and sensory proprieties . Indeed, olive secoiridoids play a role in the prevention of atherosclerosis and the inhibition of low-density lipoprotein peroxidation . Numerous studies have clearly demonstrated that these compounds exhibit cancer preventive activities  and can contribute to the nutritional prevention of osteoporosis . In particular, oleuropein, hydroxytyrosol  and oleocanthal  have shown effects on human health. The secoiridoids contribute to the quality of olive oil, influencing the oil taste, being responsible for bitter and pungency sensory notes and as primary antioxidants, secoiridoids are involved in oil oxidative stability .
Phenolics play a crucial role in the plant response to environmental cues, being the most important defence compounds against defoliating insects . They also affect shoot branching  and have been hypothesised to protect cells and prevent fungal penetration into the cambial zone [19, 20]. Some data support the idea that the resistance to specific pathogens might also be related to certain types of phenolics [21, 22]. Oleuropein is responsible for the release of phytoalexins , and it is also a multivalent alkylator that functions as an ideal protein cross-linker, exhibiting the strongest activity reported for a plant metabolite, which adversely affects herbivores by decreasing the nutritive value of dietary proteins .
To date, secoiridoid metabolism has not been well clarified, but a pathway has been proposed for some Oleaceae species [25, 26]. Secoiridoid accumulation is a controlled process with expression and composition varying considerably among varieties, tissues, developmental stages and in response to different environmental conditions .
The key genes that modulate the synthesis and degradation of secondary compounds in olive fruits have not been characterised, with the exception of a few genes involved in triterpene biosynthesis [28, 29], due to the lack of information for the olive genome sequence. The first olive fruit transcriptome data were recently released [30, 31], representing an important resource for the identification of genes involved in fruit metabolism.
Functional genetic studies are difficult to perform in perennial woody species due to the lack of efficient protocols for mutagenesis, transformation and in vitro regeneration; therefore, understanding the natural variations for traits of interest represents a valuable tool. In plant science, the integration of gene expression and metabolic data sets is currently being attempted to study metabolic pathways .
The aim of this work was to provide further insight into the evolution of the olive fruit phenolic compounds and to identify their major genetic determinants.
Results and discussion
Concentration of phenolics according to variety and fruit developmental stage
The cultivars chosen for the analyses represent a high level of variation in the fruit phenolics content, based on a number of studies previously performed on subgroups of these compounds or empirical information, as is the case for Dolce d’Andria and Tendellone, which are traditionally used as table olives and are directly edible without undergoing the debittering process. To our knowledge, our work represents the first attempt to directly compare phenolic profiles of olive cultivars grown under the same environmental conditions. These cultivars, in fact, are cultivated in different regions of Italy, under different climate conditions, soils, water availabilities and agro-techniques that greatly affect the fruit phenolic content. A previous molecular characterisation of these varieties has established that they have originated independently, and no close genetic relationships have been observed among either the high- or low-phenolics cultivars (Baldoni, unpub. data).
In the high phenolics (HP) cultivars Coratina and Rosciola, oleuropein represented the most abundant compound (up to 82% of the total), decreasing in concentration after fruit set, as previously observed in other olive varieties . In the low phenolics (LP) cultivars, Tendellone and Dolce d’Andria, the main phenolic compound was not oleuropein, lower than 15 mg g-1 dw during all stages of the maturation process (Figure 1B), but 3–4 DHPEA-EDA, which accounted for up to 94% of the total phenolic contents at the beginning of fruit development.
Demethyloleuropein showed an opposite trend, accumulating during the last period of fruit development, after 105 DAF, with the highest levels observed in the cvs. Coratina, Rosciola, Frantoio, Canino, Moraiolo, and Leccino. In contrast, only trace amounts of this compound (lower than 3 mg g-1 dw) were detected in Tendellone, Bianchella, and Dritta, and demethyloleuropein expression was completely absent in Dolce d’Andria, Nocellara del Belice and Nocellara Etnea, suggesting a block in the reaction leading to the synthesis or accumulation of this compound (Figure 1C). In HP cultivars the progressive decrease in the oleuropein concentration corresponded with an increase of demethyloleuropein content. The highest concentration of demethyloleuropein was observed at the stages with the lowest concentration of oleuropein. This observation supports the hypothesis of the derivation of demethyloleuropein from the demethylation of oleuropein, as previously envisaged by Obied et al.. In addition, the observation that in LP cultivars and Nocellara del Belice, Nocellara Etnea, Bianchella and Dritta, oleuropein expression and demethyloleuropein content remain constant during fruit development is in accordance with this hypothesis.
The compound 3–4 DHPEA-EDA represented a consistent portion of the total phenolics during the early stages of fruit development in the cvs. Tendellone, Dritta, Nocellara del Belice, and Nocellara Etnea, (95%, 31.5% 29.7% and 42.3%, respectively), but the its concentration decreased to values nearing 0 at 90 DAF for most varieties (Figure 1D). Interestingly, in Dolce d’Andria it was the only secoiridoid detected at the early stages of fruit development, suggesting that in this cultivar secoiridoid synthesis is not totally blocked and also that this compound could be positioned upstream along the pathway before oleuropein, as reported by Ryan et al..
In all cultivars expressing ligstroside, the concentration of this compound increased with fruit development (Figure 1E). This trend was opposite to that observed for oleuropein, suggesting that this compound might be derived from oleuropein rather than acting as its precursor.
Among phenolics, the tyrosol concentration showed a decreasing pattern from less than 2 mg g-1 dw, whereas the hydroxytyrosol concentration remained below 4 mg g-1 dw, with minor variation among cultivars and fruit developmental stages; notably, a slight increase in hydroxytyrosol concentration was observed after 120 DAF in Canino, Leccino, Dritta and Bianchella (Figure 2A-B, respectively).
Verbascoside, a molecule comprising tyrosol/hydroxytyrosol and cinnamyl-derivative moieties , was detected in consistent amounts (up to 90 mg g-1 dw) in Coratina and Rosciola at 90 DAF, while its concentration remained lower than 30 mg g-1 dw in the other cultivars during all developmental phases (Figure 2C), in accordance with the results of Jemai et al.. These authors detected verbascoside only in earlier maturation stages in the cv. Dhokar. Based on our data, it seems that the two major phenolic compounds of olive fruits, oleuropein and verbascoside, with common tyrosol/hydroxytyrosol moieties, do not show any clear relationship.
The lignans acetoxypinoresinol and pinoresinol (Figure 2D) showed a modest increase in expression during fruit development, with the former being greater than the latter. The strongest difference in the relative content of these compounds was observed in Dolce d’Andria, for which the acetoxypinoresinol contents were approximately ten times higher than those of pinoresinol (Additional file 4).
The metabolic analysis allowed us to select HP and LP cultivars for further transcriptional analyses, as the contrasting phenotypes of these cultivars might reflect differences in the expression of the genes involved in phenolic metabolism.
Identification of transcripts putatively involved in the synthesis of the main secondary metabolites identified from fruit EST datasets
Transcripts putatively involved in the secondary metabolism of olive fruits
1-deoxy-d-xylulose 5-phosphate synthase (EC:126.96.36.199)
7.9E-124, 92% (H. brasiliensis)
1-deoxy-d-xylulose-5-phosphate reductoisomerase (EC:188.8.131.527)
2.5E-129, 92% (H. brasiliensis)
2-C-methyl-D-erythritol 4-phosphate cytidyltransferase (EC:184.108.40.206)
4.0E-141, 87% (S. miltiorrhiza)
4-diphosphocytidyl-2-C-methyl-D-erythritol kinase (EC:220.127.116.11)
5.1E-111, 76% (S. miltiorrhiza)
2-C-methyl-D-erythritol 2,4-cyclodiphosphate synthase (EC:18.104.22.168)
5.8E-89, 75% (A. adenophora)
4-hydroxy-3-methylbut-2-en-1-yl diphosphate synthase-like (EC 22.214.171.124)
9.0E-32, 95% (V. vinifera)
4-hydroxy-3-methylbut-2-enyl diphosphate reductase (EC:126.96.36.199)
0.0, 87% (S. lycopersicum)
Isopentenyl diphosphate isomerase (EC:188.8.131.52)
2.3E-156, 95% (N. tabacum)
Mevalonate pathway (MVA)
3-hydroxy-3-methyl glutaryl CoA reductase (EC:184.108.40.206)
7.6E-40, 89% (S. nigrum)
Mevalonate kinase (EC:220.127.116.11)
4.1E-142, 84% (C. roseus)
Phosphomevalonate kinase (EC:18.104.22.168)
6.2E-24, 67% (V. vinifera)
Mevalonate diphosphate decarboxylase (EC:22.214.171.124)
1.0E-14, 85% (C. roseus)
Synthesis of monoterpenic moiety of secoiridoids
geraniol synthase (EC:4.2.3.-)
4.0E-166, 76% (P.dulcis)
Geraniol 10-hydroxylase (EC:1.14.13.B15)
0.0, 80% (C. roseus)
NADH dehydrogenase I (EC:126.96.36.199)
1.0E-92, 77% (P. trichocarpa x P. deltoides)
8.0E-126, 69% (N. tabacum)
Secologanin synthase-like (EC:188.8.131.52)
1.2E-145, 67% (V. vinifera)
Secologanin synthase-like (EC:184.108.40.206)
1.1E-120, 45% (G. max)
Secologanin synthase-like (EC:220.127.116.11)
1.8E-79, 44% (G. max)
Secologanin synthase (EC:18.104.22.168)
1.3E-44, 38% (C. roseus)
S-adenosylmethionine-dependent methyltransferase (EC:22.214.171.124)
1.9E-65, 53% (V. vinifera)
Synthesis of phenolic moieties
Arogenate dehydrogenase (EC:126.96.36.199)
6.0E-70, 84% (S. pennellii)
Copper amine oxidase (EC:188.8.131.52)
0.0, 89% (R. communis)
Polyphenol oxidase (EC:184.108.40.206)
0.0, 60% (S. indicum)
Tyrosine/dopa decarboxylase (EC:220.127.116.11)
0.0, 78% (P. somniferum)
Alcohol dehydrogenase (EC:18.104.22.168)
2.0E-163, 98% (O. europaea)
Alcohol dehydrogenase class-3 (EC:22.214.171.124)
1.0E-95, 93% (G.max)
Phenylalanine ammonia-lyase (EC:126.96.36.199)
0.0, 90% (P. frutescens)
4-coumarate coenzyme A ligase (EC: 188.8.131.52)
9.1E-50, 98% (P. fortunei)
Sterol and terpene biosynthesis
Limonene synthase like (EC: 184.108.40.206)
6.0E-17, 67% (A. thaliana)
Farnesyl diphosphate synthase (EC:220.127.116.11)
0.0, 89% (G. uralensis)
Squalene synthase (EC: 18.104.22.168)
1.2E-111, 95% (B. monnieri)
Geranylgeranyl pyrophosphate synthase (EC:22.214.171.124)
2.1E-159, 90% (C. roseus)
Degradation of phenolics
6.0E-85, 85% (R. communis)
Peroxidase (EC: 126.96.36.199)
4.0E-140, 82% (R. communis)
The unigenes showing best tBLASTn scores to protein sequences functionally characterised in other species were selected for subsequent transcriptional analyses. To restrict the number of candidates, only those with the highest levels of expression, deduced by the number of ESTs per contig, were considered. In fact, taking into account that phenolic compounds represent the most abundant secondary metabolites of olive fruit, we assumed that the metabolism should be strongly oriented to their biosynthesis. The E-value and the percentage identity to known proteins have been reported for the selected genes (Table 1, Additional file 5). Although this method has a certain level of uncertainty and further studies are necessary to confirm the proposed functions, the identification of candidate genes considering the similarity to known proteins is a powerful approach, as demonstrated for many other species.
The transcripts were resequenced for the cvs. Coratina and Dolce d’Andria and named according to their putative function.
Among the selected transcripts, some were involved in plastidial 2-C-methyl-d-erythritol 4-phosphate (MEP) and cytosolic mevalonate (MVA) pathways, whereas other candidate transcripts were putatively involved in the synthesis of secoiridoids (monoterpenic and phenolic moieties), other phenolics, terpenoids and sterols (Table 1).
Moreover, an olive gene encoding for a lupeol synthase (LUPS) [Genbank: AB025343], which resulted implicated in triterpenoid biosynthesis in yeast , was considered for the expression analyses.
Identification of other transcripts from cDNA-AFLP analysis
To identify genes that might play a role in secondary metabolism but not included in the fruit EST dataset, a cDNA-AFLP analysis was performed. Comparing the fruit transcriptional profiles of the HP cvs. Coratina and Rosciola and LP cvs. Tendellone and Dolce d’Andria in three developmental stages (45, 90 and 165 DAF) allowed the identification of 93 fragments showing differential expression patterns. For 59 fragments, high quality sequences were obtained and blasted against protein databases. Forty-five transcripts (approximately 76% of selected sequences) showed significant similarity to proteins with known function. The BLAST results and expression patterns of these transcripts are shown in Additional file 6. The sequences of cDNA-AFLP fragments matched entries in the fruit EST database, with low differences attributable to cultivar allele variations, EST sequence uncertainties, or different genes of the same family.
The results of the Blast2GO analysis allowed the annotation of the expressed sequences according to the terms of three main Gene Ontology vocabularies. The majority of the transcripts encoded for putative proteins with binding or catalytic activities (44% and 37%, respectively) (Additional file 7).
For 24 transcripts, the expression profiles were confirmed using semi-quantitative PCR (sqPCR) analysis, and 23 transcripts showed a similar expression pattern to that observed in the cDNA-AFLP analysis (Additional file 8).
Seven transcripts (17,7%) were implicated in the synthesis or degradation of secoiridoids (Table 1). The sequences of the transcripts for OeDXS, OeNDHD, OeGT, OeADH, OeCuAO, OeGLU and OePOX were extended to 574, 643, 1050, 790, 1843, 857 and 1,182 bp, respectively, using RACE-PCR. The sequences of OeNDHD and OePOX include the complete coding and partial 5’ and 3’ UTR regions.
The OeDXS sequence exhibited high identity (86% identity, 92% similarity) to DXS type II [Genbank:CAD22531] of M. truncatula, which is putatively involved in secondary metabolism, and a lower BLASTP score (79% identity, 86% similarity) to DXS type I [Genbank:CAD22530], which is proposed to play a role in primary metabolism .
Quantitative expression analyses for olive fruit mRNAs
The expression of the 28 candidate genes identified from EST datasets and the seven selected ones from the cDNA-AFLP analysis, putatively encoding for enzymes involved in secoiridoid, phenolic, terpene and sterol metabolism, was characterised using RT-qPCR to detect a possible correlation with the metabolic data.
The analysis was performed at 45, 90 and 165 DAF using the fruits of HP (Coratina) and LP (Dolce d’Andria) cultivars.
Relative expression of transcripts putatively involved in MEP and MVA pathways
In higher plants, the five-carbon building blocks of all terpenoids, isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP) may derive from the plastid-localised MEP pathway and the cytosolic MVA pathway. Enzymes and related genes in both pathways are well known and have been characterised.
The four transcripts putatively encoding enzymes involved in the MVA pathway (OeHMGR, OeMVAK, OeMVAPK, OeMVAPPD) showed completely different profiles compared with the MEP transcripts. In particular, OeHMGR and OeMVAK did not display strong differences among the developmental stages or varieties, except for a weak expression increase at 165 DAF in cv. Coratina (Figure 4).
Two isoforms of the enzyme HMGR were identified in Coffea arabica fruits and transcripts of isoform CaHMGR1 were expressed only at the initial stages of fruit development, while the isoform CaHMGR2 was constitutively expressed . It is also possible that these two isoforms are expressed in the olive fruit, and it is likely that we have identified the constitutive form.
Coratina OeMVAPK was more expressed at fruit pit hardening (90 DAF), whereas in Dolce d’Andria, it was expressed at the same level in all three developmental stages. OeMVAPPD mRNA increased during fruit development and was always more highly expressed (up to 13-fold) in Coratina than in Dolce d’Andria.
Unlike transcripts involved in the MVA pathway, the expression patterns of transcripts involved in the MEP pathway, in accordance with secoiridoid, decrease during fruit development, suggesting that the MEP pathway is also present in olives and might contribute to the terpenoid portion of secoiridoids. These results are consistent with those reported in other plant species, where MEP and MVA pathways produce different terpenoid classes. In fact, it is generally accepted that geranyl diphosphate (GDP) and geranylgeranyl diphosphate (GGDP), deriving from the MEP pathway, are used in plastids as substrates for monoterpene, diterpene and tetraterpene synthesis, whereas farnesyl diphosphate (FDP), obtained from the MVA pathway in the cytosol, is involved in triterpene and sesquiterpene biosynthesis . However, cross talk between these two different IPP biosynthetic pathways has been documented, and the relative contribution of each pathway to the biosynthesis of the various classes of terpenes remains uncertain . In olive, only a contrasting report has been published , suggesting the involvement of the MVA pathway in secoiridoid formation in the Oleaceae family, but neither enzymatic nor molecular data were provided.
Relative expression of transcripts putatively involved in biosynthesis of the terpenic moiety of secoiridoids
Biosynthetic steps leading to formation of the terpenic and phenolic portions of secoiridoids are still not well clarified and, as a consequence, enzymes involved in these pathways remain uncertain.
GES is a monoterpene synthase involved in the synthesis of geraniol , and GE10H is a cytochrome P450 monooxygenase that hydroxylates the monoterpenoid geraniol at the C-10 position to generate 10-hydroxygeraniol. This enzyme was reported to be involved in the biosynthesis of iridoid monoterpenoids and several classes of monoterpenoid alkaloids observed in a diverse range of plant species . Feeding experiments on Syringa and Fraxinus (other genera within the Oleaceae family) showed that the biosynthesis of oleoside-type secoiridoids proceeds via iridodial [25, 26, 52], which presumably originates from geraniol and 10-hydroxygeraniol [1, 49, 53] (Additional file 10). Interestingly, it has been recently shown using in vitro enzymatic assays that the same OeGES1 that we have analysed is involved in the synthesis of geraniol , demonstrating the effectiveness of our approach in identifying candidate transcripts for phenolics biosynthesis.
The G10H of C. roseus (CrG10H) also catalyses the 3'-hydroxylation of naringenin to produce eriodictyol with a catalytic activity efficiency that is 10 times lower compared to geraniol hydroxylation. These studies demonstrated that G10H plays an important role in the biosynthesis of flavonoids, in addition to its previously described role in the metabolism of terpenoids . Based on these studies, we cannot exclude that OeG10H might also be involved in flavonoid biosynthesis in the olive.
Despite the fact that NDHI and GT are enzymes required in different pathways, the strong differential expression of OeNDHI and OeGT in LP and HP varieties and during different stages of fruit development, suggested a putative role for these transcripts in secoiridoid metabolism. These enzymes might play a role in various steps of the pathway. NDHI might work in both, the conversion of iridotrial to deoxyloganic acid aglucone and the conversion of 7-epi-loganic acid to 7-ketologanic acid. GTs are required for terpenic and phenolic fractions, transferring glucosylic groups to deoxyloganic aglucones for the formation of deoxyloganic acid, the conversion of oleoside 11-methyl ester to 7-ß-1-D-glucopyranosyl 11-methyloleoside and catalysing the formation of oleuropein from oleuropein aglycone, which is the last step of the pathway (Additional file 10).
The position of 3,4-DHPEA-EDA along the secoiridoid pathway remains controversial. It has been considered either as a derivative of oleuropein, produced by its enzymatic degradation by endogenous β-glucosidases  (Additional file 2), or as the intermediate compound of an alternative biosynthetic pathway leading to oleuropein formation [1, 8, 55]. The high level of expression of genes encoding enzymes promoting the conversion of 3,4-DHPEA-EDA to oleuropein, such as OeGT, during the early stages of fruit development in the HP cv. Coratina supports the hypothesis of 3,4-DHPEA-EDA as an oleuropein precursor. This hypothesis finds a further confirmation in the observation that 3,4-DHPEA-EDA is the only secoiridoid detected in cv. Dolce d’Andria at the early stages of fruit development, suggesting that a downstream block may prevent the formation of the other secoiridoids.
In O. europaea both epoxides of secologanin and secoxyloganin could be precursors of oleuropein [49, 53]. The oxidation required for the conversion of 7-ketologanin to oleoside-11- methyl ester is similar to the mechanism taking place for the conversion of loganin, an epimer of 7-ketologanin, to secologanin. Therefore, we assumed that the gene encoding this enzymatic function might have high sequence similarity with the secologanin synthase (a cytochrome P450 enzyme). Moreover, the role played by loganin and secologanin in secoiridoid biosynthesis in Olea europaea remains controversial, and it cannot be excluded that these compounds might be intermediates of secoiridoid biosynthesis , further supporting the involvement of a secologanin synthase.
In the fruit EST database, we identified four transcripts (OeSLS1 OeSLS2 OeSLS3 OeSLS4) showing high similarities to secologanin synthase. Using RT-qPCR analyses, OeSLS2 OeSLS3 OeSLS4 showed a pattern similar to the other transcripts that were included in the secoiridoid synthesis, supporting the hypothesis of their involvement in this pathway. However, OeSLS1 showed a low level of expression at the first sampling, which increased later. This result could be explained by the putative involvement of OeSLS1 in the biosynthesis of other secondary compounds, such as terpene indole alkaloids, as reported in other plant species , rather than playing a major role in oleuropein synthesis.
The methylation of 7-ketologanic acid might be catalysed by an enzyme similar to loganic acid O-methyltransferase (LAMT), which converts loganic acid to loganin (epimers of 7-ketologanic acid and 7-ketologanin, respectively), as indicated by the functional characterization in C. roseus (CrLAMT), showing that this enzyme exhibits high specificity for the loganic acid substrate . Therefore, we identified a putative olive homolog (OeLAMT) of the LAMT gene.
Relative expression of transcripts putatively involved in phenolics biosynthesis
The phenolic moiety of secoiridoids is presumably derived from tyrosine and proceeds through tyrosol [1, 8] (Additional file 11). In most plants, tyrosine is synthesised from arogenate decarboxylated by arogenate dehydrogenase (ADH) [57, 58], and hydroxytyrosol is synthesised from tyrosine through DOPA and dopamine. Recently, the biosynthesis of hydroxytyrosol was clarified in Olea europaea using cultured cells . Another pathway for the tyrosol formation has been reported in other plant species, where tyrosol might be produced from a p-coumaric acid precursor, which is derived primarily from phenylalanine . However, in olive, the presence of this alternative pathway has never been demonstrated. Instead, two alternative routes from tyrosol to oleuropein have been proposed: one considering ligstroside as direct oleuropein precursor  and the other proceeding via oleuropein aglycone  (Additional file 11). Alternative biosynthetic pathways are proposed for verbascoside, e.g., from tyramine via dopamine or from tyrosol via hydroxytyrosol .
A schematic representation of the putative olive phenolics biosynthetic pathway is proposed in Figure 5, according to the previous findings discussed above. Based on these findings, we searched for transcripts putatively involved in the biosynthesis of the phenolic portion of secoiridoids and other related phenolic compounds, such as tyrosol, hydroxytyrosol and verbascoside.
We selected and analysed the following genes that might function in the considered pathway: OeADH OeCuAO OePPO OeTYRD OeALDH1 and OeALDH2. ADH decarboxylates the arogenate producing tyrosine . CuAO deaminates various compounds with biologically active amines producing their corresponding aminoaldehydes, H2O2 and NH3. PPO catalyses the o-hydroxylation of monophenols to o-diphenols, and it might be involved in different steps of phenolics metabolism . An enzyme similar to a tyrosine/dopa decarboxylase (TYRD) is required for both the conversion of tyrosine in tyramine and DOPA in dopamine [64, 65]. The conversion of both tyramine to tyrosol and dopamine to hydroxytyrosol requires an amino-oxidase (AO) and an alcohol dehydrogenase (ALDH) [59, 60].
The expression profiles of four of the six genes (OeADH, OePPO, OeTYRD and OeALDH1) correlated with the secoiridoid content, decreasing during fruit development, similarly to the genes involved in the terpenic moiety (Figure 5B). mRNAs of OeADH, OePPO and OeTYRD were exclusively present at 45 DAF, and OeALDH1 and OeALDH2 were more highly expressed at the first sampling and strongly decreased with fruit development. OeTYRD was more expressed in Coratina at 45 DAF compared to Dolce d’Andria. Only OeCuAO was more highly expressed at 90 DAF, decreasing at 165 DAF in Coratina, whereas OeCuAO expression remained constant during development in Dolce d’Andria.
In phenylpropanoid metabolism, two key genes have been investigated: PAL and 4CL. They are involved in the synthesis of large groups of compounds, such as flavonoids, lignans and verbascoside. The expression of OePAL increased at 165 DAF in both cultivars, whereas Oe4CL showed the highest expression only at the first sampling in both cultivars (Figure 5C). The strong increase of OePAL mRNA levels at 165 DAF is consistent with the accumulation of anthocyanins and the change in fruit colour to purple-black that occurs in this phase . The opposite trend observed for Oe4CL did not correlate with the content of lignans and flavonoids. In fact, the concentration of lignans first increased and then decreased during fruit growth in all cultivars, whereas the main flavonoids, rutin and luteolin 7-O-glucoside, are reported to increase during olive fruit ripening . It is possible that Oe4CL might not be associated with lignans or flavonoid formation but with that of unique phenylpropanoid end-products, as reported for some species, such as Arabidopsis, aspen and soybean .
Relative expression of transcripts putatively involved in phenolics degradation
The expression of putative OeGLU and OePOX orthologs was analysed. GLU and POX enzymes are involved in phenolic degradation; moreover, GLU plays a role in the formation of oleuropein and ligstroside derivatives  (Additional file 2).
OeGLU was almost exclusively expressed during the early stages of fruit development (45 DAF) of cultivar Coratina, whereas OePOX was almost exclusively expressed at 45 DAF in both cultivars (Figure 5D). These profiles are similar to those observed for transcripts putatively involved in secoiridoid synthesis. These results confirm the role of these enzymes in processes that lead to the decrease in phenolic concentration observed at 90 and 165 DAF, and their expression might be down-regulated whenever a lower availability of oleuropein, their main substrate, occurs . A similar mechanism might explain the lower expression observed for both types of transcripts in the LP cultivar, where the lack of oleuropein could be due to differences in the regulation of enzymes involved in its biosynthesis rather than in its degradation.
Relative expression of transcripts putatively involved in terpenoid and sterol biosynthesis
OeFPPS and OeSQS mRNA levels, putatively coding for two key enzymes leading to the synthesis of squalene, were also analysed, and their expression reached a peak at 90 DAF in Coratina (Figure 6B). It is noteworthy that important compounds are synthesised from squalene-oxide at the onset of fruit development, such as maslinic and oleanolic acids and α- and β-amyrins . The expression profiles observed for OeFPPS and OeSQS of Coratina were consistent with these findings. In contrast, differences among developmental stages were not observed in Dolce d’Andria.
The expression profile in developing fruits of the gene for lupeol synthase, catalysing the formation of triterpenic lupeol and functionally characterised in yeast , has never been analysed. Our OeLUPS showed a different mRNA profile between the two cultivars at 45 DAF (Figure 6B), which was more highly expressed in Coratina than in Dolce d’Andria, and, in the last two stages, they shared the same pattern. These results could be explained by a different accumulation of squalene and triterpenes during fruit development in Dolce d’Andria; however, the accumulation of these compounds in this variety has not yet been investigated.
The mRNA levels of GGPS, catalysing the synthesis of geranyl geranyl pyrophosphate (GGPP), an important intermediate for diterpenes and carotenoids formation, were higher at 45 DAF compared with the other developmental stages in both cultivars, and, at the same stage, the expression was higher in Coratina (Figure 6C). Our results are consistent with the higher accumulation of carotenoids during the first stages of olive fruit development .
Our study represents an effort to characterise the transcriptional profile of candidate genes putatively involved in secondary metabolism in olive fruits by assessing the content of major phenolic compounds in a comprehensive number of cultivars at different fruit developmental stages. These methods allowed us to identify candidate genes for secondary metabolites.
We observed a strong variation in overall metabolite content, in the level of specific compounds and in the accumulation trends among genotypes and fruit developmental stages. Coratina and Rosciola were characterised by the highest oleuropein concentration at 45 DAF. In contrast, Dolce d’Andria and Tendellone showed the lowest content at all stages. Interestingly, we identified genotypes characterised by the complete absence of other compounds. In particular, Nocellara del Belice, Dolce d’Andria and Nocellara Etnea did not accumulate demethyloleuropein, and Dolce d’Andria was the only variety showing undetectable amounts of all secoiridoids, excepting 3,4-DHPEA-EDA at the early stages of fruit development. These findings might result from a block in the enzymatic steps leading to the biosynthesis of these compounds.
The evidence of an accumulation of ligstroside occurring when oleuropein decreases, leaves space to the hypothesis that ligstroside may derive from oleuropein instead of the contrary.
The observed differences in the phenolic profiles among the cultivars might reflect great variability in the modulation of their biosynthesis and accumulation. This variability can be exploited in breeding programmes to increase the fruit composition of important phenolic compounds.
The recent development of EST datasets for olive fruits allowed the extrapolation of gene information using sequence-similarity-based approaches. We used protein sequences that were previously characterised or had an assigned function in other species. The high level of similarity detected allowed us to predict the enzyme classes and, with some degree of approximation, the substrate specificity. As a complementary approach, the cDNA-AFLP analysis proved to be an efficient technique for the isolation of differentially expressed transcripts without any sequence similarity assumption.
Based on these approaches, 36 transcripts were identified and their expression profiles were characterised and associated with corresponding metabolite profiles.
The strong correlation observed between the content of specific metabolites during fruit development and the expression of transcripts putatively involved in their biosynthesis, suggests that metabolite content is regulated at transcriptional level and strengthens the involvement of the candidate genes in the proposed pathways. The differences observed in the expression of some genes between Coratina and Dolce d’Andria might indicate a different regulation of the transcripts involved in the secondary metabolism among olive genotypes.
In particular, the levels of most of transcripts putatively involved in the biosynthesis of secoiridoids (both terpenic and phenolic moieties) showed a strong decrease during fruit development, according to the decrease in oleuropein concentration at the same stages, and many genes involved in secoiridoid pathway were more highly expressed in HP cv. Coratina than in LP cv. Dolce d’Andria.
We observed a correlation between the MEP pathway and secoiridoid transcriptional profiles, supporting the hypothesis that this pathway, rather than the MVA pathway, primarily contributes to secoiridoid biosynthesis.
Interestingly, some key genes of monoterpenoid and phenolic biosynthesis, such as OeDXS, OeGES, OeG10H and OeADH, were exclusively expressed in the earliest sampling, when the highest secoiridoid concentration was detected. These genes might regulate the accumulation of these compounds during fruit development.
A strong correlation between metabolic and transcriptional data was also identified for the biosynthesis of limonene and GGPP. OeLS was almost exclusively expressed at the developmental stage when the highest level of this compound was detected. OeGGPPS, putatively involved in carotenoid biosynthesis, was expressed only in young fruits when the highest levels of these compounds occur, as reported in other studies.
These data provide useful information for functional genetic studies of this crop species and for the identification of functional markers related to the accumulation of compounds and metabolites affecting the nutraceutical and organoleptic properties of olive fruits and oil.
Based on previous information on the phenolic profile of their oils, the following cultivars were chosen, putatively representing a high level of variation in fruit phenolics content: Bianchella, Canino, Coratina, Dolce d’Andria, Dritta, Frantoio, Leccino, Moraiolo, Nocellara del Belice, Nocellara Etnea, Rosciola and Tendellone. Fruits used for the phenolic composition and cDNA-AFLP analyses were harvested from 45 to 165 days after full bloom (DAF) every 15 days from plants of an olive cultivar collection at the experimental farm of the CRA–OLI (Collececco, Spoleto, Perugia) in central Italy (42° 48’ 48”N, 12° 39’ 15”E, 356 m above sea level). The phenological stage of the fruits at sampling dates was recorded. Plants were grown under the same environmental and agronomical conditions. To avoid possible effects of different levels of water availability on the phenolic content among trees, their water status was periodically monitored and occasional irrigation was applied as needed to maintain all plants at similar values of pre-dawn water potential during the sampling period (data not shown). Immediately after harvesting, the olive fruits were frozen in liquid nitrogen and stored at −80°C until further analysis.
Evaluation of phenolic compounds
The extraction of phenolic compounds was performed according to Servili et al., with few modifications, Briefly, 3 g of olive fruit mesocarp and exocarp were homogenised using a Homogeniser A/S N (Foss Electric, Denmark) in a 100-ml solution of methanol/water 80:20%, followed by two further homogenisation in 50 ml of methanol/water 80:20 using an Ultra-Turrax T 25 (IKA, Staufen, Germany). After methanol evaporation in vacuum under a nitrogen flow at 37°C, solid-phase extraction (SPE) was performed to separate the phenolics from the aqueous extract. During the SPE, a 900-mg Extraclean high load C18 cartridge (Alltech Italia s.r.l., Sedriano, Italy) was loaded with 1 ml of olive extract using 50 ml of methanol as the eluting solvent.
The HPLC analysis was performed according to Selvaggini et al., using a Spherisorb ODS-1 250 mm x 4.6 mm column with a particle size of 5 μm (Phase Separation Ltd., Deeside, UK).
The phenolic separation was performed using semi-preparative high-performance liquid chromatography (HPLC) analysis with a 9.4 mm i.d. 500-mm Whatman Partisil 10 ODS-2 semipreparative column; the mobile phase was 0.2% acetic acid in water (pH 3.1) (A)/methanol (B) at a flow rate of 6.5 mL/min and phenol detection was performed using a diode array detector (DAD). The purity of all compounds obtained from direct extraction was tested using HPLC, and their chemical structure was verified using nuclear magnetic resonance (NMR) with the same operative conditions reported in previous studies . The following phenolic compounds were considered: oleuropein, demethyloleuropein, 3,4-DHPEA-EDA, ligstroside, tyrosol, hydroxytyrosol, verbascoside and lignans.
Identification of mRNAs putatively involved in the metabolism of phenolic compounds
Sequences derived from the OLEA EST database , a collection of over 102,000 Olea europaea L. fruit EST reads generated through 454 massive sequencing technology , were analysed to identify transcripts putatively involved in the pathways of terpenoids, phenolics and other secondary metabolite synthesis. Amino acid sequences of genes involved in the pathways of interest and published in the Kyoto Encyclopaedia of Genes and Genomes (KEGG) databases [74, 75] were used to search olive homologs in the OLEA EST database using basic local alignment (tBLASTn). For those genes represented using multiple unique transcripts, only the largest contig and/or that representing the highest number of singletons was used for the analysis. Sequences obtained from this collection include candidates for the isoprenoid (both MVP and MEP pathways), phenylpropanoid, terpene (monoterpenes, secoiridoids, diterpenes and triterpenes), phenol, sterol, lignan and flavonoid biosynthesis, leading to the main secondary metabolites present in the olive fruit. All transcripts of interest were resequenced in the HP and LP cultivars, Coratina and Dolce D’Andria, respectively, using gene-specific primers (Additional file 12).
Total RNA was isolated from fruit mesocarp and exocarp using the RNeasy Plant Mini Kit (Qiagen), and contaminating genomic DNA were removed with DNase I (Qiagen) treatment. Samples of two HP (Coratina and Rosciola) and two LP (Tendellone and Dolce d’Andria) cultivars, at three samplings (45, 90 and 165 DAF), were considered.
PolyA RNA was isolated from approximately 50 μg of total RNA using oligo dT Dynabeads (Invitrogen). All purified mRNA was used to synthesise first strand cDNA using SuperScript III (Invitrogen) according to the manufacturer’s 18-bp oligo d(T) protocol. Double-stranded (ds) cDNA was synthesised by incubating the first-strand product for 2 h at 16°C with 30 U of DNA polymerase I (Invitrogen) and 3 U of RNase H (Invitrogen) in a reaction mixture containing 20 mM Tris–HCl, 75 mM KCl, 10 mM (NH4)2SO4, 5 mM MgCl2 and 1 mM DTT. Subsequently, a 10-μl aliquot of each sample was assessed on a 1% agarose gel, and a clear DNA smear was visible between 500 and 4,000 bp. The samples were purified using a phenol/chloroform procedure and quantified using a spectrophotometer at a wavelength of 260 nm.
The cDNA-AFLP procedure was conducted according to Bachem and coworkers , with some modifications. A total of 500 ng of cDNA was digested in a 50-μl volume using 10 U of MseI and 20 U of EcoRI for 2 h at 37°C. The digestion mix was ligated to 5 pmol of EcoRI adapter and 50 pmol of MseI adapter using 68 U of T4 Ligase (New England Biolabs). The reaction was performed for 2 h at 37°C. The PCR reaction solution (50 μl) for preamplification contained 10 μl of the digestion mix, 75 ng of each primer constructed using the adaptors, 0.2 mM dNTP mix, 1.5 mM Mg2+ and 1 U of Taq polymerase (Invitrogen). The PCR reaction was conducted using the following conditions: 94°C for 2 min, 25 cycles at 94°C for 30 sec, 56°C for 1 min and 72°C for 1 min, followed by 72°C for 7 min. The selective amplification was performed using 5 μl of preamplification product diluted 1:10, 60 ng of each selective primer, 1.5 mM Mg2+, 0.2 mM dNTP mix and 1 U of Taq polymerase in a 20 μl total reaction volume. A total of 20 combinations of the selective primers containing 3 and 2 selective bases for the primer were designed using the EcoRI and MseI adapters, respectively, and the Eco primers were labelled with fluorescein. Touch-down PCR was performed using the following cycling parameters: 94°C for 2 min, 13 cycles at: 94°C for 30 sec, 65°C (−0.7°C/cycle) for 30 sec and 72°C for 1 min, 23 cycles at 94°C for 30 sec, 56°C for 30 sec and 72°C for 1 min, followed by 72°C for 7 min. The selected PCR product (5.7 μl) was denatured at 96°C for 12 min and separated on a 4% denatured polyacrylamide gel, pre-run at 6 Watts for 10 min and run at 125 Watts for 1.75 h at 50°C in Genomyx GX 100 (Beckman Instruments, Fullerton, CA). The gel was scanned, and the separated fragments were visualised on a gel using fluorescent emission. Any differentially expressed fragment was excised from the gel, treated in a speed vacuum for 30 sec and incubated in a 100 μl volume of water overnight at 4°C. A total of 5 μl of each sample was reamplified with 50 ng of the same primers used in the preamplification procedure, 1.5 mM Mg2+, 0.2 mM dNTP mix and 1 U of Taq polymerase in a 20 μl total reaction volume. The PCR reaction was conducted using the following conditions: 94°C for 2 min; 35 cycles at: 94°C for 30 sec, 56°C for 1 min and 72°C for 1 min, followed by 72°C for 7 min. A 5 μl aliquot of the PCR mixture was run on an agarose gel, and the remaining reactions were purified with PolyEthylene Glycol (PEG) 8000. The cDNA fragments were eluted in 10 μl of H2O and either directly sequenced or sequenced after subcloning into the pGEM-T Easy Vector (Promega) followed by plasmid purification with QIAprep Spin Miniprep kit (Qiagen). The fragments were sequenced using the automated sequencer 3130 Genetic Analyser (Applied Biosystems, Foster City, CA).
The sequences were aligned to the non-redundant protein sequences (nr) database at the NCBI  using BLASTX, which also calculated protein similarity and predicted the hypothetical function of the unknown fragment sequences.
All sequences differentially expressed in cDNA-AFLP analysis were annotated according to the three main gene ontological categories (cellular component, biological process and molecular function) using the Blast2GO software v1.3.3 [78, 79] according to Galla and coworkers , with minor modifications. All sequences were also aligned to the OLEA EST database using the BlastN algorithm to detect possible homologous transcripts.
Enzyme mapping of the annotated sequences was performed using a direct GO to Enzyme annotation, and this programme was also used to query the KEGG maps to define the main metabolic pathways involved. The conserved domains were identified using the Conserved Domain Database (CDD) [80, 81].
Sequences showing expression profiles compatible with the content of corresponding metabolites were further characterised. Fragments were extended at 3’ and 5’ ends using a RACE-PCR method to obtain long partial or full-length cDNA clones. The resulting sequences were used to confirm the putative function of each compound. RACE-PCR was performed using the SMART RACE cDNA Amplification Kit (Clontech), with some modifications. Briefly, total RNA (2 μg) was reverse transcribed using the 3’CDS primer II A and the SMART II A oligonucleotide (Clontech). The reactions were incubated in a solution containing 50 mM Tris–HCl pH 8.3, 75 mM KCl, 6 mM MgCl2 and 2 mM DTT for 1.5 h at 42°C with Superscript III (Invitrogen). For the second strand synthesis, touch-down PCR was performed using a small aliquot (1/20 volume) of the primary template with EX Taq (Takara), the 5’ PCR primer II A (Clontech) and a gene specific primer. A thermal cycling programme was performed using the following conditions: an initial denaturation step of 94°C for 2 min, followed by 16 cycles of 94°C for 30 sec, 68°C for 30 sec with a 0.5°C decrease per cycle and 72°C for 3 min. An additional 20 cycles were performed at 94°C for 30 sec, 60°C for 30 sec and 72°C for 3 min, with a final elongation step at 72°C for 5 min. Nested PCR was required in some cases to improve the specificity of the amplification. The samples were purified using PEG 8000 and either directly sequenced or sequenced after being cloned into the pGEM-T Easy vector with the 3130 Genetic Analyser (Applied Biosystems, Foster City, CA).
To validate the differential expression patterns of cDNA-AFLP selected clones, preliminary semi-quantitative PCR experiments were performed on 24 transcripts using different biological replicates of the samples. The RNeasy Plant Mini Kit (Qiagen) was used to isolate 2 μg of total RNA from biological replicates of the olive samples used in the cDNA-AFLP analysis. The cDNA synthesis was performed using SuperScript III according to the manufacturer’s protocol. Elongation Factor (EF1α) was used as the endogenous gene for sample normalisation . Specific EF1α primers were used for cDNA sample normalisation with the following conditions: a 25 μl total reaction volume containing 10 μl of first strand cDNA diluted 1:50, 1.5 mM Mg2+, 0.2 mM dNTP mix and 1.25 U of Taq polymerase (Invitrogen). The following PCR amplification conditions were used: 94°C for 3 min, 20 to 45 cycles (depending on the amplicon signal intensity obtained under non-saturating PCR conditions) at 94°C for 30 sec, 60°C for 30 sec and 72°C for 50 sec, followed by 72°C for 7 min. The cDNAs were normalised based upon band intensity. After sample normalisation, a PCR was performed using the previously described conditions with 16-bp specific primers that included the restriction site sequences (EcoI and MseI) used for the cDNA-AFLP analysis and a fragment specific region (Additional file 13) to validate both the expression and the allelic polymorphisms of any genes of interest.
The RT-qPCR experiments were performed on the most representative candidates of main branches of secondary metabolites maps. RT-qPCR was performed on the subset of most robust 37 putative transcripts to verify if their expression profiles corresponded to the patterns of secoiridoid, phenolics and other secondary metabolites synthesis or degradation. Coratina and Dolce d’Andria were used as reference HP and LP varieties, respectively. Total RNA was extracted from 0.2 g of fruit mesocarp and exocarp with the RNeasy Plant Mini Kit (Qiagen) and treated with DNase I (Qiagen). Reverse transcription of 2 μg of RNA was performed using oligo(dT)18 and the SuperScript III Reverse Transcriptase kit (Invitrogen) according to the manufacturer’s instructions. Quantitative real-time PCR was performed on a PCR Real Time 7300 (Applied Biosystems, Foster City, CA) according to the manufacturer's protocol and using the Reagent kit for SYBR Green analysis (Applied Biosystems) and gene-specific primers (Additional file 14). Primers were verified by the presence of a single PCR product band after agarose gel electrophoresis. All reactions were performed in triplicate. After each assay, a dissociation kinetics analysis was performed to verify the specificity of the amplification products. Relative amounts of all mRNAs were calculated using the 2-ΔΔCT method , where ΔCt = Cttarget gene - Ctreference gene using Dolce D’Andria at 165 DAF as a control sample. The housekeeping Elongation Factor 1-α gene was used as an endogenous reference gene for cDNA normalisation. The data for three biological replicates were analysed using an analysis of variance (ANOVA) followed by Bonferroni’s post hoc test (P < 0.05) with R software (version 2.14.0) .
Days after flowering
Real-time reverse transcription polymerase chain reaction
2-C-methyl-D-erythritol 4-phosphate cytidyltransferase
2-C-methyl-D-erythritol 2,4-cyclodiphosphate synthase
Isopentenyl diphosphate delta isomerase
Mevalonate diphosphate decarboxylase
NADH dehydrogenase I
Loganic acid methyltransferase
Copper amine oxidase
4-coumarate coenzyme A ligase
Farnesyl diphosphate synthase
Geranyl geranyl pyrophosphate
Geranyl geranyl pyrophosphate synthase
The research was supported by the Project OLEA – Genomics and Breeding of Olive, funded by MIPAF, Italy and by the Project FISR - Improving flavour and nutritional properties of plant food after first and second transformation, Italian Ministry of Research.
- Obied HK, Prenzler PD, Ryan D, Servili M, Taticchi A, Esposto S, Robards K: Biosynthesis and biotransformations of phenol-conjugated oleosidic secoiridoids from Olea europaea L. Nat Prod Rep. 2008, 25 (6): 1167-1179. 10.1039/b719736e.PubMedView ArticleGoogle Scholar
- Beauchamp GK, Keast RS, Morel D, Lin J, Pika J, Han Q, Lee CH, Smith AB, Breslin PA: Phytochemistry: ibuprofen-like activity in extra-virgin olive oil. Nature. 2005, 437 (7055): 45-46. 10.1038/437045a.PubMedView ArticleGoogle Scholar
- Oliveras López MJ, Innocenti M, Ieri F, Giaccherini C, Romani A, Mulinacci N: HPLC/DAD/ESI/MS detection of lignans from Spanish and Italian Olea europaea L. fruits. Journal of Food Composition and Analysis. 2008, 21: 62-70. 10.1016/j.jfca.2007.04.012.View ArticleGoogle Scholar
- Owen RW, Haubner R, Mier W, Giacosa A, Hull WE, Spiegelhalder B, Bartsch H: Isolation, structure elucidation and antioxidant potential of the major phenolic and flavonoid compounds in brined olive drupes. Food Chem Toxicol. 2003, 41 (5): 703-717. 10.1016/S0278-6915(03)00011-5.PubMedView ArticleGoogle Scholar
- Servili M, Selvaggini R, Esposto S, Taticchi A, Montedoro G, Morozzi G: Health and sensory properties of virgin olive oil hydrophilic phenols: agronomic and technological aspects of production that affect their occurrence in the oil. J Chromatogr A. 2004, 1054 (1–2): 113-127.PubMedView ArticleGoogle Scholar
- Boskou G, Salta FN, Chrysostomou S, Mylona A, Chiou A, Andrikopoulos NK: Antioxidant capacity and phenolic profile of table olives from the Greek market. Food Chem. 2006, 94: 558-564. 10.1016/j.foodchem.2004.12.005.View ArticleGoogle Scholar
- Jerman T, Trebše P, Mozeticˇ Vodopivec B: Ultrasound-assisted solid liquid extraction (USLE) of olive fruit (Olea europaea) phenolic compounds. Food Chem. 2010, 123 (123): 175-182.View ArticleGoogle Scholar
- Ryan D, Antolovich M, Herlt T, Prenzler PD, Lavee S, Robards K: Identification of phenolic compounds in tissues of the novel olive cultivar hardy's mammoth. J Agric Food Chem. 2002, 50 (23): 6716-6724. 10.1021/jf025736p.PubMedView ArticleGoogle Scholar
- Servili M, Baldioli M, Selvaggini R, Macchioni A, Montedoro G: Phenolic compounds of olive fruit: one- and two-dimensional nuclear magnetic resonance characterization of Nuzhenide and its distribution in the constitutive parts of fruit. J Agric Food Chem. 1999, 47 (1): 12-18. 10.1021/jf9806210.PubMedView ArticleGoogle Scholar
- Goulas V, Manganaris GA: Towards an Efficient Protocol for the Determination of Triterpenic Acids in Olive Fruit: A Comparative Study of Drying and Extraction Methods. Phytochem Anal. 2012, 23 (5): 444-449. 10.1002/pca.1378.PubMedView ArticleGoogle Scholar
- Dabbou S, Issaoui M, Esposto S, Sifi S, Taticchi A, Servili M, Montedoro GF, Hammami M: Cultivar and growing area effects on minor compounds of olive oil from autochthonous and European introduced cultivars in Tunisia. J Sci Food Agric. 2009, 89: 1314-1325. 10.1002/jsfa.3588.View ArticleGoogle Scholar
- Waterman E, Lockwood B: Active Components and Clinical Applications of Olive Oil. Altern Med Rev. 2007, 12 (4): 331-342.PubMedGoogle Scholar
- Stiti N, Triki S, Hartmann MA: Formation of triterpenoids throughout Olea europaea fruit ontogeny. Lipids. 2007, 42 (1): 55-67. 10.1007/s11745-006-3002-8.PubMedView ArticleGoogle Scholar
- Servili M, Montedoro GF: Contribution of phenolic compounds to virgin olive oil quality. European Journal of Lipid Science and Technology. 2002, 104: 602-613. 10.1002/1438-9312(200210)104:9/10<602::AID-EJLT602>3.0.CO;2-X.View ArticleGoogle Scholar
- Omar SH: Oleuropein in olive and its pharmacological effects. Sci Pharm. 2010, 78 (2): 133-154. 10.3797/scipharm.0912-18.PubMedPubMed CentralView ArticleGoogle Scholar
- Puel C, Mathey J, Agalias A, Kati-Coulibaly S, Mardon J, Obled C, Davicco MJ, Lebecque P, Horcajada MN, Skaltsounis AL, et al: Dose–response study of effect of oleuropein, an olive oil polyphenol, in an ovariectomy/inflammation experimental model of bone loss in the rat. Clin Nutr. 2006, 25 (5): 859-868. 10.1016/j.clnu.2006.03.009.PubMedView ArticleGoogle Scholar
- Mumm R, Hilker M: Direct and indirect chemical defence of pine against folivorous insects. Trends Plant Sci. 2006, 11 (7): 351-358. 10.1016/j.tplants.2006.05.007.PubMedView ArticleGoogle Scholar
- Umehara M, Hanada A, Yoshida S, Akiyama K, Arite T, Takeda-Kamiya N, Magome H, Kamiya Y, Shirasu K, Yoneyama K, et al: Inhibition of shoot branching by new terpenoid plant hormones. Nature. 2008, 455 (7210): 195-200. 10.1038/nature07272.PubMedView ArticleGoogle Scholar
- Beckman CH: Phenolic-storing cells: keys to programmed cell death and periderm formation in wilt disease resistance and in general defence responses in plants?. Physiological and Molecular Plant Pathology. 2000, 57: 101-110. 10.1006/pmpp.2000.0287.View ArticleGoogle Scholar
- Franceschi VR, Krokene P, Christiansen E, Krekling T: Anatomical and chemical defenses of conifer bark against bark beetles and other pests. New Phytol. 2005, 167 (2): 353-375. 10.1111/j.1469-8137.2005.01436.x.PubMedView ArticleGoogle Scholar
- Ockels FS, Eyles A, McPherson BA, Wood DL, Bonello P: Phenolic chemistry of coast live oak response to Phytophthora ramorum infection. J Chem Ecol. 2007, 33 (9): 1721-1732. 10.1007/s10886-007-9332-z.PubMedView ArticleGoogle Scholar
- Wallis C, Eyles A, Chorbadjian R, McSpadden Gardener B, Hansen R, Cipollini D, Herms DA, Bonello P: Systemic induction of phloem secondary metabolism and its relationship to resistance to a canker pathogen in Austrian pine. New Phytol. 2008, 177 (3)): 767-778.PubMedView ArticleGoogle Scholar
- Kubo I, Matsumoto A, Takase I: A multichemical defense mechanism of bitter olive Olea europaea (Oleaceae). Is oleuropein a phytoalexin precursor?. J Chem Ecol. 1985, 11: 251-263. 10.1007/BF00988207.PubMedView ArticleGoogle Scholar
- Konno K, Hirayama C, Yasui H, Nakamura M: Enzymatic activation of oleuropein: a protein crosslinker used as a chemical defense in the privet tree. Proc Natl Acad Sci. 1999, 96 (16): 9159-9164. 10.1073/pnas.96.16.9159.PubMedPubMed CentralView ArticleGoogle Scholar
- Damtoft S, Franzyk H, Jensen SR: Biosynthesis of secoiridoids in Fontanesia. Phytochemistry. 1995, 38 (3): 615-621. 10.1016/0031-9422(94)00683-K.View ArticleGoogle Scholar
- Damtoft S, Franzyk H, Jensen SR: Biosynthesis of iridoids in Syringa and Fraxinus: carbocyclic iridoid precursors. Phytochemistry. 1995, 40 (3): 785-792. 10.1016/0031-9422(95)00210-X.View ArticleGoogle Scholar
- Malik NSA, Bradford JM: Recovery and stability of oleuropein and other phenolic compounds during extraction and processing of olive (Olea europaea L.) leaves. Journal of Food Agriculture and Environment. 2008, 6: 8-13.Google Scholar
- Shibuya M, Zhang H, Endo A, Shishikura K, Kushiro T, Ebizuka Y: Two branches of the lupeol synthase gene in the molecular evolution of plant oxidosqualene cyclases. Eur J Biochem. 1999, 266 (1): 302-307. 10.1046/j.1432-1327.1999.00875.x.PubMedView ArticleGoogle Scholar
- Saimaru H, Orihara Y, Tansakul P, Kang YH, Shibuya M, Ebizuka Y: Production of triterpene acids by cell suspension cultures of Olea europaea. Chem Pharm Bull(Tokyo). 2007, 55 (5): 784-788. 10.1248/cpb.55.784.View ArticleGoogle Scholar
- Alagna F, D'Agostino N, Torchia L, Servili M, Rao R, Pietrella M, Giuliano G, Chiusano ML, Baldoni L, Perrotta G: Comparative 454 pyrosequencing of transcripts from two olive genotypes during fruit development. BMC Genomics. 2009, 10: 399-10.1186/1471-2164-10-399.PubMedPubMed CentralView ArticleGoogle Scholar
- Galla G, Barcaccia G, Ramina A, Collani S, Alagna F, Baldoni L, Cultrera NG, Martinelli F, Sebastiani L, Tonutti P: Computational annotation of genes differentially expressed along olive fruit development. BMC Plant Biol. 2009, 9: 128-10.1186/1471-2229-9-128.PubMedPubMed CentralView ArticleGoogle Scholar
- Saito K, Hirai MY, Yonekura-Sakakibara K: Decoding genes with coexpression networks and metabolomics - 'majority report by precogs'. Trends Plant Sci. 2008, 13 (1): 36-43. 10.1016/j.tplants.2007.10.006.PubMedView ArticleGoogle Scholar
- Morello JR, Romero MP, Motilva MJ: Effect of the maturation process of the olive fruit on the phenolic fraction of drupes and oils from Arbequina, Farga, and Morrut cultivars. J Agric Food Chem. 2004, 52 (19): 6002-6009. 10.1021/jf035300p.PubMedView ArticleGoogle Scholar
- Jemai H, Bouaziz M, Sayadi S: Phenolic composition, sugar contents and antioxidant activity of Tunisian sweet olive cultivar with regard to fruit ripening. J Agric Food Chem. 2009, 57 (7): 2961-2968. 10.1021/jf8034176.PubMedView ArticleGoogle Scholar
- The OLEA EST database: The OLEA EST database. http://188.8.131.52/oleaestdb/,
- Collu G, Unver N, Peltenburg-Looman AM, van der Heijden R, Verpoorte R, Memelink J: Geraniol 10-hydroxylase, a cytochrome P450 enzyme involved in terpenoid indole alkaloid biosynthesis. FEBS Lett. 2001, 508 (2): 215-220. 10.1016/S0014-5793(01)03045-9.PubMedView ArticleGoogle Scholar
- Irmler S, Schroder G, St-Pierre B, Crouch NP, Hotze M, Schmidt J, Strack D, Matern U, Schroder J: Indole alkaloid biosynthesis in Catharanthus roseus: new enzyme activities and identification of cytochrome P450 CYP72A1 as secologanin synthase. Plant J. 2000, 24 (6): 797-804. 10.1046/j.1365-313x.2000.00922.x.PubMedView ArticleGoogle Scholar
- Murata J, Roepke J, Gordon H, De Luca V: The leaf epidermome of Catharanthus roseus reveals its biochemical specialization. Plant Cell. 2008, 20 (3): 524-542. 10.1105/tpc.107.056630.PubMedPubMed CentralView ArticleGoogle Scholar
- Walter MH, Hans J, Strack D: Two distantly related genes encoding 1-deoxy-d-xylulose 5-phosphate synthases: differential regulation in shoots and apocarotenoid-accumulating mycorrhizal roots. Plant J. 2002, 31 (3): 243-254. 10.1046/j.1365-313X.2002.01352.x.PubMedView ArticleGoogle Scholar
- Cordoba E, Salmi M, Leon P: Unravelling the regulatory mechanisms that modulate the MEP pathway in higher plants. J Exp Bot. 2009, 60 (10): 2933-2943. 10.1093/jxb/erp190.PubMedView ArticleGoogle Scholar
- Estevez JM, Cantero A, Reindl A, Reichler S, Leon P: 1-Deoxy-D-xylulose-5-phosphate synthase, a limiting enzyme for plastidic isoprenoid biosynthesis in plants. J Biol Chem. 2001, 276 (25): 22901-22909. 10.1074/jbc.M100854200.PubMedView ArticleGoogle Scholar
- Xing S, Miao J, Li S, Qin G, Tang S, Li H, Gu H, Qu LJ: Disruption of the 1-deoxy-D-xylulose-5-phosphate reductoisomerase (DXR) gene results in albino, dwarf and defects in trichome initiation and stomata closure in Arabidopsis. Cell Res. 2010, 20 (6): 688-700. 10.1038/cr.2010.54.PubMedView ArticleGoogle Scholar
- Zulak KG, Bohlmann J: Terpenoid biosynthesis and specialized vascular cells of conifer defense. J Integr Plant Biol. 2010, 52 (1): 86-97. 10.1111/j.1744-7909.2010.00910.x.PubMedView ArticleGoogle Scholar
- Battilana J, Costantini L, Emanuelli F, Sevini F, Segala C, Moser S, Velasco R, Versini G, Stella Grando M: The 1-deoxy-D: -xylulose 5-phosphate synthase gene co-localizes with a major QTL affecting monoterpene content in grapevine. Theor Appl Genet. 2009, 118 (4)): 653-669.PubMedView ArticleGoogle Scholar
- Emanuelli F, Battilana J, Costantini L, Le Cunff L, Boursiquot JM, This P, Grando MS: A candidate gene association study on muscat flavor in grapevine (Vitis vinifera L.). BMC Plant Biol. 2010, 10: 241-10.1186/1471-2229-10-241.PubMedPubMed CentralView ArticleGoogle Scholar
- Tiski I, Marraccini P, Pot D, Vieira LG, Pereira LF: Characterization and expression of two cDNA encoding 3-Hydroxy-3-methylglutaryl coenzyme A reductase isoforms in coffee (Coffea arabica L.). Omics. 2011, 15 (10)): 719-727.PubMedView ArticleGoogle Scholar
- Aharoni A, Jongsma MA, Kim TY, Ri MB, Giri AP, Verstappen FWA, Schwab W, Bouwmeester HJ: Metabolic engineering of terpenoid biosynthesis in plants. Phytochemistry Reviews. 2006, 5: 49-58. 10.1007/s11101-005-3747-3.View ArticleGoogle Scholar
- Dudareva N, Andersson S, Orlova I, Gatto N, Reichelt M, Rhodes D, Boland W, Gershenzon J: The nonmevalonate pathway supports both monoterpene and sesquiterpene formation in snapdragon flowers. Proc Natl Acad Sci U S A. 2005, 102 (3): 933-938. 10.1073/pnas.0407360102.PubMedPubMed CentralView ArticleGoogle Scholar
- Damtoft S, Franzyk H, Jensen SR: Biosynthesis of secoiridoid glucosides in Oleaceae. Phytochemistry. 1993, 34 (5): 1291-1299. 10.1016/0031-9422(91)80018-V.View ArticleGoogle Scholar
- Iijima Y, Gang DR, Fridman E, Lewinsohn E, Pichersky E: Characterization of geraniol synthase from the peltate glands of sweet basil. Plant Physiol. 2004, 134 (1): 370-379. 10.1104/pp.103.032946.PubMedPubMed CentralView ArticleGoogle Scholar
- Sung PH, Huang FC, Do YY, Huang PL: Functional expression of geraniol 10-hydroxylase reveals its dual function in the biosynthesis of terpenoid and phenylpropanoid. J Agric Food Chem. 2011, 59 (9): 4637-4643. 10.1021/jf200259n.PubMedView ArticleGoogle Scholar
- Jensen SR, Franzyk H, Wallander E: Chemotaxonomy of the Oleaceae: iridoids as taxonomic markers. Phytochemistry. 2002, 60 (3): 213-231. 10.1016/S0031-9422(02)00102-4.PubMedView ArticleGoogle Scholar
- Soler-Rivas C, Espìn JC, Wichers HJ: Oleuropein and related compounds. J Sci Food Agric. 2000, 80: 1013-1023. 10.1002/(SICI)1097-0010(20000515)80:7<1013::AID-JSFA571>3.0.CO;2-C.View ArticleGoogle Scholar
- Vezzaro A, Krause ST, Nonis A, Ramina A, Degenhardt J, Ruperti B: Isolation and characterization of terpene synthases potentially involved in flavor development of ripening olive (Olea europaea) fruits. J Plant Physiol. 2012, 169 (9): 908-914. 10.1016/j.jplph.2012.01.021.PubMedView ArticleGoogle Scholar
- Ryan D, Prenzler PD, Lavee S, Antolovich M, Robards K: Quantitative changes in phenolic content during physiological development of the olive (Olea europaea) cultivar Hardy's Mammoth. J Agric Food Chem. 2003, 51 (9): 2532-2538. 10.1021/jf0261351.PubMedView ArticleGoogle Scholar
- O'Connor SE, Maresh JJ: Chemistry and biology of monoterpene indole alkaloid biosynthesis. Nat Prod Rep. 2006, 23 (4): 532-547. 10.1039/b512615k.PubMedView ArticleGoogle Scholar
- Legrand P, Dumas R, Seux M, Rippert P, Ravelli R, Ferrer JL, Matringe M: Biochemical characterization and crystal structure of Synechocystis arogenate dehydrogenase provide insights into catalytic reaction. Structure. 2006, 14 (4): 767-776. 10.1016/j.str.2006.01.006.PubMedView ArticleGoogle Scholar
- Rippert P, Puyaubert J, Grisollet D, Derrier L, Matringe M: Tyrosine and phenylalanine are synthesized within the plastids in Arabidopsis. Plant Physiol. 2009, 149 (3): 1251-1260. 10.1104/pp.108.130070.PubMedPubMed CentralView ArticleGoogle Scholar
- Saimaru H, Orihara Y: Biosynthesis of acteoside in cultured cells of Olea europaea. J Nat Med. 2010, 64 (2): 139-145. 10.1007/s11418-009-0383-z.PubMedView ArticleGoogle Scholar
- Zhang JX, Ma LQ, Yu HS, Zhang H, Wang HT, Qin YF, Shi GL, Wang YN: A tyrosine decarboxylase catalyzes the initial reaction of the salidroside biosynthesis pathway in Rhodiola sachalinensis. Plant Cell Rep. 2011, 30 (8): 1443-1453. 10.1007/s00299-011-1053-7.PubMedView ArticleGoogle Scholar
- Rippert P, Matringe M: Molecular and biochemical characterization of an Arabidopsis thaliana arogenate dehydrogenase with two highly similar and active protein domains. Plant Mol Biol. 2002, 48 (4): 361-368. 10.1023/A:1014018926676.PubMedView ArticleGoogle Scholar
- Sebela M, Frebort I, Lemr K, Brauner F, Pec P: A study on the reactions of plant copper amine oxidase with C3 and C4 aliphatic diamines. Arch Biochem Biophys. 2000, 384 (1): 88-99. 10.1006/abbi.2000.2081.PubMedView ArticleGoogle Scholar
- Ortega-Garcia F, Blanco S, Peinado MA, Peragon J: Polyphenol oxidase and its relationship with oleuropein concentration in fruits and leaves of olive (Olea europaea) cv. 'Picual' trees during fruit ripening. Tree Physiol. 2008, 28 (1)): 45-54.PubMedView ArticleGoogle Scholar
- Guillet G, Poupart J, Basurco J, De Luca V: Expression of tryptophan decarboxylase and tyrosine decarboxylase genes in tobacco results in altered biochemical and physiological phenotypes. Plant Physiol. 2000, 122 (3): 933-943. 10.1104/pp.122.3.933.PubMedPubMed CentralView ArticleGoogle Scholar
- Facchini PJ, De Luca V: Differential and tissue-specific expression of a gene family for tyrosine/dopa decarboxylase in opium poppy. J Biol Chem. 1994, 269 (43): 26684-26690.PubMedGoogle Scholar
- Agati G, Pinelli P, Cortes Ebner S, Romani A, Cartelat A, Cerovic ZG: Nondestructive evaluation of anthocyanins in olive (Olea europaea) fruits by in situ chlorophyll fluorescence spectroscopy. J Agric Food Chem. 2005, 53 (5)): 1354-1363.PubMedView ArticleGoogle Scholar
- Gomez-Rico A, Salvador MD, Fregapane G: Effect of cultivar and ripening on minor components in Spanish olive fruits and their corresponding virgin olive oils. Food Res Int. 2008, 41: 433-440. 10.1016/j.foodres.2008.02.003.View ArticleGoogle Scholar
- Kumar A, Ellis BE: 4-coumarate:CoA ligase gene family in Rubus idaeus: cDNA structures, evolution, and expression. Plant Mol Biol. 2003, 51 (3): 327-340. 10.1023/A:1022004923982.PubMedView ArticleGoogle Scholar
- Sivakumar G, Uccella N, Briccoli Bati C: Demethyloleuropein and beta-glucosidase activity in olive fruits. Biotechnol J. 2007, 2 (3): 381-385. 10.1002/biot.200600118.PubMedView ArticleGoogle Scholar
- Flamini G, Cioni PL, Morelli I: Volatiles from leaves, fruits, and virgin oil from Olea europaea Cv. Olivastra Seggianese from Italy. J Agric Food Chem. 2003, 51 (5)): 1382-1386.PubMedView ArticleGoogle Scholar
- Roca M, Minguez-Mosquera MI: Carotenoid levels during the period of growth and ripening in fruits of different olive varieties (Hojiblanca, Picual and Arbequina). J Plant Physiol. 2003, 160: 451-459. 10.1078/0176-1617-00759.PubMedView ArticleGoogle Scholar
- Servili M, Badioli M, Selvaggini R, Miniati E, Macchioni A, Montedoro GF: High-Performance Liquid Chromatography evaluation of phenols in olive fruits, virgin olive oil, vegetation waters, and pomace and ID- and 2D-Nuclear Magnetic Resonance characterization. J Am Oil Chem Soc. 1999, 7 (76): 873-882.View ArticleGoogle Scholar
- Selvaggini R, Servili M, Urbani S, Esposto S, Taticchi A, Montedoro G: Evaluation of phenolic compounds in virgin olive oil by direct injection in high-performance liquid chromatography with fluorometric detection. J Agric Food Chem. 2006, 54 (8): 2832-2838. 10.1021/jf0527596.PubMedView ArticleGoogle Scholar
- The Kyoto Encyclopaedia of Genes and Genomes (KEGG). http://www.genome.jp/kegg,
- Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T, et al: KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008, 36 (Database issue)): D480-D484.PubMedPubMed CentralGoogle Scholar
- Bachem CW, van der Hoeven RS, de Bruijn SM, Vreugdenhil D, Zabeau M, Visser RG: Visualization of differential gene expression using a novel method of RNA fingerprinting based on AFLP: analysis of gene expression during potato tuber development. Plant J. 1996, 9 (5): 745-753. 10.1046/j.1365-313X.1996.9050745.x.PubMedView ArticleGoogle Scholar
- The National Center for Biotechnological Information (NCBI). http://blast.ncbi.nlm.nih.gov,
- Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21 (18): 3674-3676. 10.1093/bioinformatics/bti610.PubMedView ArticleGoogle Scholar
- The Blast2GO software,. http://www.blast2go.org]
- The Conserved Domain Database (CDD): The Conserved Domain Database (CDD). http://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi]
- Marchler-Bauer A, Anderson JB, Chitsaz F, Derbyshire MK, DeWeese-Scott C, Fong JH, Geer LY, Geer RC, Gonzales NR, Gwadz M, et al: CDD: specific functional annotation with the Conserved Domain Database. Nucleic Acids Res. 2009, 37 (Database issue)): D205-D210.PubMedPubMed CentralView ArticleGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.PubMedView ArticleGoogle Scholar
- The R project. www.r-project.org,
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.