BACKGROUND
fruit development is the result of genetically programmed processes influenced by environmental factors.
to identify and characterize genes involved in these processes, different genomic approaches  have been used in several fruit species  <cit>  and the body of information concerning transcriptional networks and regulatory circuits involved in important physiological and developmental processes increased tremendously during the last two decades.
in tomato, large-scale est sequencing projects resulted in a better insight into molecular mechanisms of fruit ripening processes and in the identification of common transcription factors not previously associated with ripening  <cit> .
generation of ests and consequent discovery of genes with potential roles in fruit development have also been reported in grape berry  <cit> .
in apple, an extensive analysis has been made using all est sequences available in public databases to identify genes temporally or spatially regulated during fruit growth and development  <cit> .
other extensive est sequencing projects focusing on fruit development have been set up in peach  <cit> , melon  <cit>  and kiwifruit  <cit> .
sequence information derived from advanced est sequencing is an essential resource for functional genomics studies based on the use of microarray technology and real-time pcr.
following the pioneering work of aharoni and co-workers  <cit>  on strawberry, several papers have now been published on the use of microarrays in different fruit species.
olea europaea l. is an evergreen species, traditionally cultivated in the mediterranean basin.
the oil that results from mechanical extraction of the fruits is a predominant component of the worldwide known 'mediterranean diet', to which increasing attention is being paid for its health benefits and cancer-protective properties  <cit> .
these attributes are closely related to the oil composition and to the concentration of active bio-molecules resulting from the catabolic and anabolic processes taking place throughout olive fruit development which is a long process lasting several months.
the oil content of olives can reach up to 30%  at full ripening  <cit> : it accumulates in the mesocarp and, at a lower extent, in the seed  <cit> .
oil accumulation in the pulp increases slowly, reaching the plateau after veraison.
a marked tryacylglycerol  accumulation in seed and pulp occurs after endocarp lignification, when about  <dig> mg of oil per fruit per week can be synthesized.
the fatty acid profile of the oil accumulating in the fruit is important in relation to its nutritional properties  <cit> .
the main fatty acid is oleic acid , which represents about 75% of total fatty acids, followed by linoleic , palmitic , stearic  and linolenic  acid.
the pattern of fatty acid synthesis and desaturation varies during maturation and ripening, according to cultivars and to environmental conditions  <cit> .
other important metabolites accumulate throughout olive fruit development.
they include polyphenols  <cit> , carotenoids  <cit> , chlorophylls  <cit> , sterols and terpenoids  <cit>  all directly or indirectly affecting olive oil quality and its technological and nutritional properties.
information concerning genetic regulation of these metabolic processes in olive is still very limited.
only few genes involved in fatty acid metabolism have been characterized  <cit> .
a monosaccharide transporter , whose expression increases during fruit maturation, when a massive accumulation of sugars occurs, has been recently cloned  <cit> .
moreover, the gene encoding a geranylgeranyl reductase  has been isolated and its role in organ development and stress response in relation to tocopherol action hypothesized  <cit> .
information at molecular level about polyphenol and triterpenoid metabolism is lacking, as well as the mechanisms involved in olive fruit development and ripening.
among different strategies available for identifying differentially expressed genes, suppression subtractive hybridization  libraries have been successfully used in fruit science to elucidate mechanisms regulating anthocyanin metabolism in grape berries  <cit> , proanthocyanidin biosynthesis in persimmon  <cit> , processes involved in early growth and ethylene-induced ripening in banana  <cit> , and in orange pigmentation  <cit> .
this paper deals with the identification via ssh of large repertories of differentially expressed genes in developing olive fruits, and their computational annotation by means of different bioinformatic software.
the identification and characterization of gene regulatory networks and key metabolic pathways during fruit growth and development represent a prerequisite for improving olive oil quality and its health-related properties.
RESULTS
the study was based on the preparation of cdna libraries using ssh that, likewise to the differential display , represents an efficient strategy to isolate genes with an antagonist expression pattern.
this technique enabled to identify transcripts of genes differentially expressed among the three different developmental stages of olive fruit corresponding to initial fruit set , completed pit hardening  and veraison  .
as far as the composition of the four subtractive libraries is concerned, the number of differentially expressed sequences randomly chosen varied from a minimum of  <dig> to a maximum of  <dig> per library, with a total number of clones equal to  <dig>  .
the average length of the cdna clones was  <dig> bp with a wide range of variation, from  <dig> up to  <dig>  bp.
the redundancy within each single library was relatively low, ranging between  <dig> % and  <dig> %.
taking into account the whole set of sequences, the overall redundancy calculated among the four libraries was equal to  <dig> %.
the sequences of each single library were preliminarily analyzed using the cap <dig> program in order to isolate the singlets and assemble contiguous and overlapping clones into contigs.
this affected the comparative redundancy that increased up to  <dig> %.
querying with cdna sequences the non-redundant ncbi databases allowed the attribution of a blast hit of 79%, 91%, 88% and 78% of the clones belonging to the a, b, c and d libraries, respectively .
the average sequence similarity was around 76%, ranging from 72% to 80%, and the median e-value for each single library ranged from 1e- <dig> to 1e- <dig>
the max e-value was equal to 1e-- <dig> corresponding to the cut-off adopted for go annotation.
around 75% of the blast hits of the olive fruit cdna sequences were homologous to coding sequences present in the rice, arabidopsis and grapevine genomes, with more than  <dig>  hits per species.
it is worthy to note that until now  only  <dig> blast hits for olive could be recorded .
the computational analysis of the whole est collection using the software blast2go allowed the annotation of the expressed sequences according to the terms of the three main gene ontology vocabularies, i.e.
cellular compartment, molecular function and biological process .
as far as cellular compartments are concerned, the most represented are plastids and mitochondria, with more than 50% of the total annotations, followed by cytosol, plasma membrane, endoplasmic reticulum and nucleoplasm, whereas other cellular compartments were represented at a much lower scale .
concerning the molecular function, the most represented categories were those of nucleotide binding proteins, followed by proteins with transport, kinase and enzymatic activities.
the other molecular functions were represented at a lower extent .
more than  <dig> categories were found for the biological process vocabulary, being carbohydrate metabolism, response to biotic and environmental stresses, generation of precursors, metabolites and energy, and catabolic processes the most represented .
although numerically less represented, it is worth to mention the presence of terms related to the secondary metabolites, metabolism of lipids, synthesis of amino acids and derivatives, metabolites and their precursors, and protein modification process.
focusing on the go annotation of each single subtractive library, a number of olive fruit stage-specific go terms were identified .
among the  <dig> and  <dig> go terms found in the a and b libraries,  <dig> and  <dig> were associated to down- and up-regulated genes, respectively.
the most significant go terms were encoding elements of hormone biosynthesis and signal transduction mediated by ethylene, jasmonic acid, salicylic acid, and abscisic acid, as well as biosynthesis of secondary metabolites, such as terpenoids.
in the c and d libraries, a total of  <dig> and  <dig> go terms were recovered,  <dig> and  <dig> of which were related to down- and up-regulated genes, respectively.
among these, there are go terms associated to environmental stress responses, catabolism of secondary metabolites , response to hormones , and auxin signal transduction.
the analysis of go terms shared by pair-wise library combinations allows to determine which genes are continuously or transiently down- or up-regulated during the studied process.
this analysis retrieved only  <dig> terms for the down-regulated genes, whereas as many as  <dig> were those collected among the up-regulated ones.
on the contrary, comparable numbers  of go terms associated to transiently up- and down-regulated genes among the three fruit developmental stages were found .
quantitative real-time pcr experiments were carried out to corroborate the expression patterns of a subset of sequences , corresponding to  <dig> different genes.
expression patterns related to the selected gene sequences and estimated in pair-wise comparisons between the three different fruit stages are reported in figure  <dig>  the real-time pcr analyses validated the results from ssh experiments for  <dig> out of  <dig> genes .
the validated genes were grouped according to different expression patterns .
the kyoto encyclopaedia of genes and genomes  was queried for sequences encoding enzymes and the deduced gene products were associated to specific metabolic and/or biosynthetic pathways related to carbohydrates, fatty acid and secondary metabolism.
the most represented kegg maps associated with carbohydrate and fatty acid biosynthesis and metabolism were organized in a simple network analysis representative of major communication ways among retrieved metabolic pathways .
as expected, several genes encoding enzymes related to carbohydrate and fatty acid compounds were transcriptionally up- or down-regulated during olive fruit development.
several components of the carbon fixation in photosynthetic organisms , and starch and sucrose metabolism  were modulated in their expression in the early stages of fruit development .
kegg maps analysis pointed out an intense up-regulation of the majority of enzymes related to the pentose phosphate pathway , glycolysis and gluconeogenesis  along with starch and sucrose metabolism  .
table  <dig> reports the enzymes involved in starch and sucrose metabolism, glycolysis and gluconeogenesis.
transcripts of enzymes involved in the synthesis of pyruvate from β-d-fructose-6-p, such as 6-phosphofructokinase , glyceraldehyde-3-phosphate dehydrogenase , phosphoglycerate kinase , phosphopyruvate hydratase  and pyruvate kinase  were up-regulated from  <dig> to  <dig> daf and down-regulated from  <dig> daf to  <dig> daf .
* g: glycolysis and gluconeogenesis; c: citrate cycle ; p1: pentose phosphate pathway; p2: pentose and glucuronate interconversions; f: fructose and mannose metabolism; s: starch and sucrose metabolism; py: pyruvate metabolism; cf: carbon fixation; a: ascorbate and aldarate metabolism
consistent with the accumulation of oil in the drupe that starts around pit hardening and reaches the highest rate well before ripening begins, several transcripts encoding enzymes leading to fatty acid biosynthesis  from glycolysis , through the pyruvate metabolism  and citrate cycle  were up-regulated .
it is worthy to mention that several enzymes were found to be up-regulated from  <dig> to  <dig> daf, when the lipid accumulation shows an increasing rate.
a similar dynamics in terms of specific transcript accumulation was recorded for several components of the fa biosynthesis as acetyl-coa carboxylase  and enoyl-  reductase , actively involved in the synthesis of malonylcoa and fa chain elongation, respectively .
moreover, enzymes controlling fa chain elongation were up-regulated throughout fruit development, whereas the synthesis of short chain fas was impaired during early development.
the malonylcoa-acp transacylase , the enzyme involved in one of the early steps of fa biosynthesis, and the acyl-  desaturase , related to the synthesis of hexadecenoyl-  and octadecenoyl- , the precursors of palmitic and oleic acid, respectively, were up-regulated at veraison as well as a number of genes involved in the fa metabolism .
* b: fatty acid biosynthesis; m: fatty acid metabolism; g: glycerolipid metabolism; gp: glycerophospholipid metabolism; e: ether lipid metabolism; a: arachidonic acid metabolism; l: linoleic acid metabolism; al: alpha-linolenic acid metabolism; u: biosynthesis of unsaturated fatty acids.
kegg maps were also produced for secondary metabolic pathways.
considering the flavonoids as an example, the isoforms of four different enzymes  controlling flavone and flavonol as well as anthocyanin biosynthesis were mapped .
transcripts encoding these enzymes were more abundant during the transition from  <dig> to  <dig> daf, thus suggesting an increased accumulation of the related metabolites during late fruit development.
a transient down-regulation of the dihydrokaempferol 4-reductase gene  was recorded at  <dig> daf .
two isoforms of a peroxidase specifically involved in four different steps of the phenylpropanoid biosynthesis as well as few enzymes involved in alkaloid biosynthesis , limonene and pinene degradation  and caffeine metabolism  were also identified .
the analysis of the olive est dataset with the mapman software enabled to reconstruct overview metabolism maps  and to group sequences in the main regulatory networks .
as far as the in silico expression analysis of all olive est clones linked to the arabidopsis genechip sequences , most genes proved to be associated to constantly and transiently regulated genes during fruit development involved in cell wall synthesis and breakdown, fatty acid biosynthesis and lipid breakdown, starch and sucrose metabolism, glycolysis, and secondary metabolism .
gene products involved in amino acid biosynthesis and metabolism were also identified .
furthermore, several genes encoding transcription factors were mainly down-regulated throughout fruit development, while some others were related to protein modification and degradation .
genes related to hormone biosynthesis and action appeared to be differently regulated according to the type of hormone.
a down-regulation of genes involved in auxin biosynthesis and metabolism, as an oxido/reductase and an iaa-amino acid synthase, occurred between  <dig> and  <dig> daf, while during late development the expression of auxin responsive factors such as arf <dig> and arf <dig>  were up- and down-regulated, respectively.
by contrast, the synthesis of abscisic acid  was clearly stimulated throughout development, since two key enzymes encoded by aba <dig>  and aao <dig> , were up-regulated at completed pit-hardening and veraison, respectively.
in addition, a down-regulation of a gene encoding a ga-regulated protein occurred between completed pit hardening and veraison.
during the same developmental phase, an up-regulation of the arr <dig>  a protein involved in the citokinin response, as well as a down-regulation of an ugt <dig>  encoding a udp-glucosyl transferase, actively involved in the metabolism of the hormone, were observed.
genes involved in jasmonic acid  metabolism were up-regulated throughout fruit development.
a transcriptional up-regulation of an ethylene receptor of the ers type as well as a down-regulation of genes involved in brassinosteroid  biosynthesis and action were observed during early development.
taking into account genes related to biotic and abiotic stress responses, dual multiple contingency tests were performed to identify go terms associated to sequences significantly and antagonistically distributed between libraries.
the results of this analysis evidenced the abundance of sequences related to light stimulus  and radiation , as well as to biotic and abiotic stresses .
sequences related to light stimulus and radiation resulted to be up-regulated from  <dig> to  <dig> daf, while those associated with response to oxidative stresses were down-regulated during the same developmental period.
discussion
olea europaea l. is a common tree species of the mediterranean basin that plays a peculiar role in the landscape characterization and represents a major agricultural commodity as source of olive oil.
olives are not only a significant food source, but also contribute to human health and are becoming popular in health-conscious diets far beyond the mediterranean area of olive oil traditional use.
taking into account the increasing worldwide commercial interest of olive oil and the lack of information on its genomic features, a functional genomic approach, able to gain insights into the genetic and molecular aspects controlling fruit development and ripening, may be considered of primary interest.
a number of physiological and biochemical data is available on growth, development and ripening of the olive drupe but, unlike other fruit species such as peach, apple and grape, information on olive gene sequences and gene products is very limited in the main public gene databases.
this study was carried out on cv.
leccino, one of the most widespread italian varieties, characterized by a relatively short fruit developmental cycle, and a high degree of synchronization of processes defining ripening.
the ssh approach allowed the identification of  <dig>  differentially expressed gene sequences in three selected developmental stages of the olive fruit, named initial fruit set , completed pit hardening  and veraison .
all the sequences were deposited in the ncbi databases on august  <dig>   <dig> .
this contribution represented a major implementation of the public olive est repertories that currently include a total number of  <dig>  sequences .
the ssh technique confirmed its suitability for transcriptome studies and for identifying genes differentially expressed along fruit development in olive.
in fact, of the  <dig>  differentially expressed sequences as many as  <dig> demonstrated to be unique sequences.
among these,  <dig>  corresponding to  <dig> different key genes, further investigated by quantitative real-time pcr, pointed out a validation of ssh findings as high as 69%
olive fruits showed significant differences for the distribution of go terms among ontological vocabularies and categories in relation to the developmental stages.
go terms involved in carbohydrate, fatty acid and flavonoid metabolism were analyzed by setting up kegg maps.
as far as carbohydrates are concerned, genes involved in carbohydrate metabolism were modulated in their expression and, in particular, the expression pattern of genes related to starch metabolism was coherent with a temporary role of starch as storage compound during fruit early development.
sequences of enzymes involved within the pentose phosphate pathway, glycolysis and gluconeogenesis, along with starch and sucrose metabolism, were up-regulated.
the up-regulation of genes encoding enzymes involved in the synthesis of pyruvate from different substrates occurring during early development, may indicate that the fruit at this stage is highly energy demanding.
there are two sources of assimilates for fruit growth in olive.
the major source is certainly the sugars translocated in the phloem from leaves or sites of storage, comprising mannitol, raffinose, stachyose, and sucrose  <cit> .
the secondary source is sugars formed by photosynthesis in developing fruits themselves that remain green for a considerable period and retain active chlorophyll even when they change colour as approaching maturity.
while chlorophyll is mostly in the exocarp, the mesocarp has been shown to contain significant amounts of phosphoenol pyruvate  <cit> , the co <dig> fixation enzyme of the cam and c <dig> photosynthetic pathway.
genes related to fatty acid  biosynthesis appeared to be up-regulated throughout development although transcripts of specific enzymes accumulate at different extent depending on the developmental stage.
enzymes controlling fa chain elongation were up-regulated throughout fruit development: the synthesis of short chain fas was impaired during early development whereas the synthesis of the precursors of palmitoleic and oleic acids was up-regulated during late development.
this finding is consistent with the oil accumulation pattern in olive mesocarp, that starts around 40- <dig> daf and reaches at ripening the highest amount, confirming what previously observed  <cit> .
a number of genes involved in the fa metabolism appeared to be differentially expressed throughout development with a significant up-regulation at veraison.
this might be interpreted as a homeostatic reaction to the large fa accumulation occurring at this developmental stage.
taking into account the secondary metabolites, genes involved in phenylpropanoid and alkaloid biosynthesis and caffeine, limonene and pinene metabolism appeared to be differentially expressed throughout fruit development.
dihydrokaempferol 4-reductase, flavanone 3-dioxygenase, naringenin-chalcone synthase, and leucocyanidin oxygenase, four enzymes controlling flavone and flavonol, as well as anthocyanin biosynthesis, were up-regulated from  <dig> to  <dig> daf, thus suggesting an increased accumulation of the related metabolites during late development.
flavonoids are important secondary metabolites precursor of flavonols and anthocyanins, the latter being responsible of the color development occurring at ripening.
key genes related to the anthocyanin biosynthetic pathway  proved to be up-regulated at veraison stage.
a similar coordinated up-regulation of these genes has been observed in grape berries at the onset of ripening, at the time of pigmentation changes  <cit> .
the relatively high abundance of transcripts related to hormones supports their key regulatory role in olive fruit development, as demonstrated in several other fruits  <cit> .
the synthesis of aba is clearly stimulated throughout fruit development, as demonstrated by the up-regulation of two key enzymes of its biosynthetic pathway.
the pattern of expression changed according to the type of hormone.
the observed down-regulation of genes involved in auxin biosynthesis and metabolism is consistent with a lowering of the auxin content reported in other fruits throughout development.
the different regulation of auxin responsive factors, such as arf <dig> and arf <dig>  which act as negative and positive regulators of iaa responsive genes, respectively  <cit> , might imply a decreased sensitivity to the hormone during late fruit development.
this is consistent with a negative regulation of auxin on the onset of the ripening syndrome observed in non climacteric fruits.
at late developmental stage, a down regulation of the zeatin o-glucosyltransferase  <dig>  actively involved in ck metabolism, as well as an up-regulation of arr <dig>  a protein involved in the ck response, have been observed.
arr <dig> is a type-b arr transcription factor involved in ck-responsive phenomena.
it has been proposed that arr <dig>  together with arr <dig> and arr <dig>  redundantly play pivotal roles in the ahk-dependent phosphorelay signaling in response to ck  <cit> .
taking into account that ck metabolic enzymes are up-regulated by the hormone, these transcriptional changes may reflect a lowering of ck concentration, along with an increase of fruit sensitivity to ck occurring during late development.
the veraison is a developmental stage characterized by a strengthening of the fruit sink action, that, also in olive, might be regulated by a complex hormone cross-talk and interaction, as demonstrated in other non climacteric fruits as grape berries  <cit> .
in spite of any signal variation in terms of expression of genes related to ethylene biosynthesis, an up-regulation of an ers type ethylene receptor has been observed during early fruit development.
this might imply an increase in sensitivity of the fruitlet to ethylene, the hormone involved in the regulation of the immature fruit physiological drop that in olive occurs between fruit set and pit hardening  <cit> .
genes involved in the jasmonates  metabolism and brassinosteroids  biosynthesis were up- and down regulated, respectively.
if the gene expression pattern is mirrored by a similar evolution of ja and br concentration in fruit tissues, this would indicate that the role of jas and brs in olive is different from other fruit types such as peach and grape.
in fact, in peach it has been demonstrated that jas delay fruit ripening  <cit> , while brs stimulate the onset of veraison in grape berry  <cit> .
the dual multiple contingency tests allowed the identification of go terms related to biotic and abiotic stress responses significantly and antagonistically distributed between libraries.
the transcriptional profile of genes related to light, temperature and biotic stresses parallels the dynamics of day length and light intensity, both peaking at the pit hardening stage.
moreover, these data may indicate that genes responding to environmental stimuli are transcriptionally regulated up to pit hardening, and post-transcriptionally up to the veraison.
CONCLUSIONS
as a concluding remark, the ssh technique allowed to identify a set of  <dig> differentially expressed unique sequences.
among these,  <dig>  corresponding to  <dig> different key genes were further investigated by real-time pcr, pointing out a validation of the ssh results as high as 69%.
the bioinformatic annotation of all gene sequences was useful to shed light on metabolic pathways and to understand specific regulatory networks.
in fact, data here reported represent a significant contribution to the elucidation of transcriptional aspects related to carbohydrates, fas, secondary metabolites, transcription factors and hormones as well as response to biotic and abiotic stresses.
particularly interesting are data related to hormones, pointing out the complexity of the role played by these compounds in olive fruit development and ripening.
these molecular and bioinformatic data represent a first step toward both functional genomics and systems biology research for understanding the gene functions and regulatory networks in olive fruit growth, development and ripening.

 fruit development is the result of genetically programmed processes influenced by environmental factors. to identify and characterize genes involved in these processes, different genomic approaches have been used in several fruit species and the body of information concerning transcriptional networks and regulatory circuits involved in important physiological and developmental processes increased tremendously during the last two decades. in tomato, large-scale est sequencing projects resulted in a better insight into molecular mechanisms of fruit ripening processes and in the identification of common transcription factors not previously associated with ripening. generation of ests and consequent discovery of genes with potential roles in fruit development have also been reported in grape berry. in apple, an extensive analysis has been made using all est sequences available in public databases to identify genes temporally or spatially regulated during fruit growth and development. other extensive est sequencing projects focusing on fruit development have been set up in peach, melon and kiwifruit. sequence information derived from advanced est sequencing is an essential resource for functional genomics studies based on the use of microarray technology and real-time pcr. following the pioneering work of aharoni and co-workers on strawberry, several papers have now been published on the use of microarrays in different fruit species. olea europaea l. is an evergreen species, traditionally cultivated in the mediterranean basin. the oil that results from mechanical extraction of the fruits is a predominant component of the worldwide known'mediterranean diet', to which increasing attention is being paid for its health benefits and cancer-protective properties. these attributes are closely related to the oil composition and to the concentration of active bio-molecules resulting from the catabolic and anabolic processes taking place throughout olive fruit development which is a long process lasting several months. the oil content of olives can reach up to 30% at full ripening: it accumulates in the mesocarp and, at a lower extent, in the seed. oil accumulation in the pulp increases slowly, reaching the plateau after veraison. a marked tryacylglycerol accumulation in seed and pulp occurs after endocarp lignification, when about mg of oil per fruit per week can be synthesized. the fatty acid profile of the oil accumulating in the fruit is important in relation to its nutritional properties. the main fatty acid is oleic acid, which represents about 75% of total fatty acids, followed by linoleic, palmitic, stearic and linolenic acid. the pattern of fatty acid synthesis and desaturation varies during maturation and ripening, according to cultivars and to environmental conditions. other important metabolites accumulate throughout olive fruit development. they include polyphenols, carotenoids, chlorophylls, sterols and terpenoids all directly or indirectly affecting olive oil quality and its technological and nutritional properties. information concerning genetic regulation of these metabolic processes in olive is still very limited. only few genes involved in fatty acid metabolism have been characterized. a monosaccharide transporter, whose expression increases during fruit maturation, when a massive accumulation of sugars occurs, has been recently cloned. moreover, the gene encoding a geranylgeranyl reductase has been isolated and its role in organ development and stress response in relation to tocopherol action hypothesized. information at molecular level about polyphenol and triterpenoid metabolism is lacking, as well as the mechanisms involved in olive fruit development and ripening. among different strategies available for identifying differentially expressed genes, suppression subtractive hybridization libraries have been successfully used in fruit science to elucidate mechanisms regulating anthocyanin metabolism in grape berries, proanthocyanidin biosynthesis in persimmon, processes involved in early growth and ethylene-induced ripening in banana, and in orange pigmentation. this paper deals with the identification via ssh of large repertories of differentially expressed genes in developing olive fruits, and their computational annotation by means of different bioinformatic software. the identification and characterization of gene regulatory networks and key metabolic pathways during fruit growth and development represent a prerequisite for improving olive oil quality and its health-related properties.  the study was based on the preparation of cdna libraries using ssh that, likewise to the differential display, represents an efficient strategy to isolate genes with an antagonist expression pattern. this technique enabled to identify transcripts of genes differentially expressed among the three different developmental stages of olive fruit corresponding to initial fruit set, completed pit hardening and veraison. as far as the composition of the four subtractive libraries is concerned, the number of differentially expressed sequences randomly chosen varied from a minimum of to a maximum of per library, with a total number of clones equal to. the average length of the cdna clones was bp with a wide range of variation, from up to bp. the redundancy within each single library was relatively low, ranging between % and %. taking into account the whole set of sequences, the overall redundancy calculated among the four libraries was equal to %. the sequences of each single library were preliminarily analyzed using the cap program in order to isolate the singlets and assemble contiguous and overlapping clones into contigs. this affected the comparative redundancy that increased up to %. querying with cdna sequences the non-redundant ncbi databases allowed the attribution of a blast hit of 79%, 91%, 88% and 78% of the clones belonging to the a, b, c and d libraries, respectively. the average sequence similarity was around 76%, ranging from 72% to 80%, and the median e-value for each single library ranged from 1e- to 1e- the max e-value was equal to 1e-- corresponding to the cut-off adopted for go annotation. around 75% of the blast hits of the olive fruit cdna sequences were homologous to coding sequences present in the rice, arabidopsis and grapevine genomes, with more than hits per species. it is worthy to note that until now only blast hits for olive could be recorded. the computational analysis of the whole est collection using the software blast2go allowed the annotation of the expressed sequences according to the terms of the three main gene ontology vocabularies, i.e. cellular compartment, molecular function and biological process. as far as cellular compartments are concerned, the most represented are plastids and mitochondria, with more than 50% of the total annotations, followed by cytosol, plasma membrane, endoplasmic reticulum and nucleoplasm, whereas other cellular compartments were represented at a much lower scale. concerning the molecular function, the most represented categories were those of nucleotide binding proteins, followed by proteins with transport, kinase and enzymatic activities. the other molecular functions were represented at a lower extent. more than categories were found for the biological process vocabulary, being carbohydrate metabolism, response to biotic and environmental stresses, generation of precursors, metabolites and energy, and catabolic processes the most represented. although numerically less represented, it is worth to mention the presence of terms related to the secondary metabolites, metabolism of lipids, synthesis of amino acids and derivatives, metabolites and their precursors, and protein modification process. focusing on the go annotation of each single subtractive library, a number of olive fruit stage-specific go terms were identified. among the and go terms found in the a and b libraries, and were associated to down- and up-regulated genes, respectively. the most significant go terms were encoding elements of hormone biosynthesis and signal transduction mediated by ethylene, jasmonic acid, salicylic acid, and abscisic acid, as well as biosynthesis of secondary metabolites, such as terpenoids. in the c and d libraries, a total of and go terms were recovered, and of which were related to down- and up-regulated genes, respectively. among these, there are go terms associated to environmental stress responses, catabolism of secondary metabolites, response to hormones, and auxin signal transduction. the analysis of go terms shared by pair-wise library combinations allows to determine which genes are continuously or transiently down- or up-regulated during the studied process. this analysis retrieved only terms for the down-regulated genes, whereas as many as were those collected among the up-regulated ones. on the contrary, comparable numbers of go terms associated to transiently up- and down-regulated genes among the three fruit developmental stages were found. quantitative real-time pcr experiments were carried out to corroborate the expression patterns of a subset of sequences, corresponding to different genes. expression patterns related to the selected gene sequences and estimated in pair-wise comparisons between the three different fruit stages are reported in figure the real-time pcr analyses validated the results from ssh experiments for out of genes. the validated genes were grouped according to different expression patterns. the kyoto encyclopaedia of genes and genomes was queried for sequences encoding enzymes and the deduced gene products were associated to specific metabolic and/or biosynthetic pathways related to carbohydrates, fatty acid and secondary metabolism. the most represented kegg maps associated with carbohydrate and fatty acid biosynthesis and metabolism were organized in a simple network analysis representative of major communication ways among retrieved metabolic pathways. as expected, several genes encoding enzymes related to carbohydrate and fatty acid compounds were transcriptionally up- or down-regulated during olive fruit development. several components of the carbon fixation in photosynthetic organisms, and starch and sucrose metabolism were modulated in their expression in the early stages of fruit development. kegg maps analysis pointed out an intense up-regulation of the majority of enzymes related to the pentose phosphate pathway, glycolysis and gluconeogenesis along with starch and sucrose metabolism. table reports the enzymes involved in starch and sucrose metabolism, glycolysis and gluconeogenesis. transcripts of enzymes involved in the synthesis of pyruvate from β-d-fructose-6-p, such as 6-phosphofructokinase, glyceraldehyde-3-phosphate dehydrogenase, phosphoglycerate kinase, phosphopyruvate hydratase and pyruvate kinase were up-regulated from to daf and down-regulated from daf to daf. * g: glycolysis and gluconeogenesis; c: citrate cycle; p1: pentose phosphate pathway; p2: pentose and glucuronate interconversions; f: fructose and mannose metabolism; s: starch and sucrose metabolism; py: pyruvate metabolism; cf: carbon fixation; a: ascorbate and aldarate metabolism consistent with the accumulation of oil in the drupe that starts around pit hardening and reaches the highest rate well before ripening begins, several transcripts encoding enzymes leading to fatty acid biosynthesis from glycolysis, through the pyruvate metabolism and citrate cycle were up-regulated. it is worthy to mention that several enzymes were found to be up-regulated from to daf, when the lipid accumulation shows an increasing rate. a similar dynamics in terms of specific transcript accumulation was recorded for several components of the fa biosynthesis as acetyl-coa carboxylase and enoyl- reductase, actively involved in the synthesis of malonylcoa and fa chain elongation, respectively. moreover, enzymes controlling fa chain elongation were up-regulated throughout fruit development, whereas the synthesis of short chain fas was impaired during early development. the malonylcoa-acp transacylase, the enzyme involved in one of the early steps of fa biosynthesis, and the acyl- desaturase, related to the synthesis of hexadecenoyl- and octadecenoyl-, the precursors of palmitic and oleic acid, respectively, were up-regulated at veraison as well as a number of genes involved in the fa metabolism. * b: fatty acid biosynthesis; m: fatty acid metabolism; g: glycerolipid metabolism; gp: glycerophospholipid metabolism; e: ether lipid metabolism; a: arachidonic acid metabolism; l: linoleic acid metabolism; al: alpha-linolenic acid metabolism; u: biosynthesis of unsaturated fatty acids. kegg maps were also produced for secondary metabolic pathways. considering the flavonoids as an example, the isoforms of four different enzymes controlling flavone and flavonol as well as anthocyanin biosynthesis were mapped. transcripts encoding these enzymes were more abundant during the transition from to daf, thus suggesting an increased accumulation of the related metabolites during late fruit development. a transient down-regulation of the dihydrokaempferol 4-reductase gene was recorded at daf. two isoforms of a peroxidase specifically involved in four different steps of the phenylpropanoid biosynthesis as well as few enzymes involved in alkaloid biosynthesis, limonene and pinene degradation and caffeine metabolism were also identified. the analysis of the olive est dataset with the mapman software enabled to reconstruct overview metabolism maps and to group sequences in the main regulatory networks. as far as the in silico expression analysis of all olive est clones linked to the arabidopsis genechip sequences, most genes proved to be associated to constantly and transiently regulated genes during fruit development involved in cell wall synthesis and breakdown, fatty acid biosynthesis and lipid breakdown, starch and sucrose metabolism, glycolysis, and secondary metabolism. gene products involved in amino acid biosynthesis and metabolism were also identified. furthermore, several genes encoding transcription factors were mainly down-regulated throughout fruit development, while some others were related to protein modification and degradation. genes related to hormone biosynthesis and action appeared to be differently regulated according to the type of hormone. a down-regulation of genes involved in auxin biosynthesis and metabolism, as an oxido/reductase and an iaa-amino acid synthase, occurred between and daf, while during late development the expression of auxin responsive factors such as arf and arf were up- and down-regulated, respectively. by contrast, the synthesis of abscisic acid was clearly stimulated throughout development, since two key enzymes encoded by aba and aao, were up-regulated at completed pit-hardening and veraison, respectively. in addition, a down-regulation of a gene encoding a ga-regulated protein occurred between completed pit hardening and veraison. during the same developmental phase, an up-regulation of the arr a protein involved in the citokinin response, as well as a down-regulation of an ugt encoding a udp-glucosyl transferase, actively involved in the metabolism of the hormone, were observed. genes involved in jasmonic acid metabolism were up-regulated throughout fruit development. a transcriptional up-regulation of an ethylene receptor of the ers type as well as a down-regulation of genes involved in brassinosteroid biosynthesis and action were observed during early development. taking into account genes related to biotic and abiotic stress responses, dual multiple contingency tests were performed to identify go terms associated to sequences significantly and antagonistically distributed between libraries. the results of this analysis evidenced the abundance of sequences related to light stimulus and radiation, as well as to biotic and abiotic stresses. sequences related to light stimulus and radiation resulted to be up-regulated from to daf, while those associated with response to oxidative stresses were down-regulated during the same developmental period. discussion olea europaea l. is a common tree species of the mediterranean basin that plays a peculiar role in the landscape characterization and represents a major agricultural commodity as source of olive oil. olives are not only a significant food source, but also contribute to human health and are becoming popular in health-conscious diets far beyond the mediterranean area of olive oil traditional use. taking into account the increasing worldwide commercial interest of olive oil and the lack of information on its genomic features, a functional genomic approach, able to gain insights into the genetic and molecular aspects controlling fruit development and ripening, may be considered of primary interest. a number of physiological and biochemical data is available on growth, development and ripening of the olive drupe but, unlike other fruit species such as peach, apple and grape, information on olive gene sequences and gene products is very limited in the main public gene databases. this study was carried out on cv. leccino, one of the most widespread italian varieties, characterized by a relatively short fruit developmental cycle, and a high degree of synchronization of processes defining ripening. the ssh approach allowed the identification of differentially expressed gene sequences in three selected developmental stages of the olive fruit, named initial fruit set, completed pit hardening and veraison. all the sequences were deposited in the ncbi databases on august. this contribution represented a major implementation of the public olive est repertories that currently include a total number of sequences. the ssh technique confirmed its suitability for transcriptome studies and for identifying genes differentially expressed along fruit development in olive. in fact, of the differentially expressed sequences as many as demonstrated to be unique sequences. among these, corresponding to different key genes, further investigated by quantitative real-time pcr, pointed out a validation of ssh findings as high as 69% olive fruits showed significant differences for the distribution of go terms among ontological vocabularies and categories in relation to the developmental stages. go terms involved in carbohydrate, fatty acid and flavonoid metabolism were analyzed by setting up kegg maps. as far as carbohydrates are concerned, genes involved in carbohydrate metabolism were modulated in their expression and, in particular, the expression pattern of genes related to starch metabolism was coherent with a temporary role of starch as storage compound during fruit early development. sequences of enzymes involved within the pentose phosphate pathway, glycolysis and gluconeogenesis, along with starch and sucrose metabolism, were up-regulated. the up-regulation of genes encoding enzymes involved in the synthesis of pyruvate from different substrates occurring during early development, may indicate that the fruit at this stage is highly energy demanding. there are two sources of assimilates for fruit growth in olive. the major source is certainly the sugars translocated in the phloem from leaves or sites of storage, comprising mannitol, raffinose, stachyose, and sucrose. the secondary source is sugars formed by photosynthesis in developing fruits themselves that remain green for a considerable period and retain active chlorophyll even when they change colour as approaching maturity. while chlorophyll is mostly in the exocarp, the mesocarp has been shown to contain significant amounts of phosphoenol pyruvate, the co fixation enzyme of the cam and c photosynthetic pathway. genes related to fatty acid biosynthesis appeared to be up-regulated throughout development although transcripts of specific enzymes accumulate at different extent depending on the developmental stage. enzymes controlling fa chain elongation were up-regulated throughout fruit development: the synthesis of short chain fas was impaired during early development whereas the synthesis of the precursors of palmitoleic and oleic acids was up-regulated during late development. this finding is consistent with the oil accumulation pattern in olive mesocarp, that starts around 40- daf and reaches at ripening the highest amount, confirming what previously observed. a number of genes involved in the fa metabolism appeared to be differentially expressed throughout development with a significant up-regulation at veraison. this might be interpreted as a homeostatic reaction to the large fa accumulation occurring at this developmental stage. taking into account the secondary metabolites, genes involved in phenylpropanoid and alkaloid biosynthesis and caffeine, limonene and pinene metabolism appeared to be differentially expressed throughout fruit development. dihydrokaempferol 4-reductase, flavanone 3-dioxygenase, naringenin-chalcone synthase, and leucocyanidin oxygenase, four enzymes controlling flavone and flavonol, as well as anthocyanin biosynthesis, were up-regulated from to daf, thus suggesting an increased accumulation of the related metabolites during late development. flavonoids are important secondary metabolites precursor of flavonols and anthocyanins, the latter being responsible of the color development occurring at ripening. key genes related to the anthocyanin biosynthetic pathway proved to be up-regulated at veraison stage. a similar coordinated up-regulation of these genes has been observed in grape berries at the onset of ripening, at the time of pigmentation changes. the relatively high abundance of transcripts related to hormones supports their key regulatory role in olive fruit development, as demonstrated in several other fruits. the synthesis of aba is clearly stimulated throughout fruit development, as demonstrated by the up-regulation of two key enzymes of its biosynthetic pathway. the pattern of expression changed according to the type of hormone. the observed down-regulation of genes involved in auxin biosynthesis and metabolism is consistent with a lowering of the auxin content reported in other fruits throughout development. the different regulation of auxin responsive factors, such as arf and arf which act as negative and positive regulators of iaa responsive genes, respectively, might imply a decreased sensitivity to the hormone during late fruit development. this is consistent with a negative regulation of auxin on the onset of the ripening syndrome observed in non climacteric fruits. at late developmental stage, a down regulation of the zeatin o-glucosyltransferase actively involved in ck metabolism, as well as an up-regulation of arr a protein involved in the ck response, have been observed. arr is a type-b arr transcription factor involved in ck-responsive phenomena. it has been proposed that arr together with arr and arr redundantly play pivotal roles in the ahk-dependent phosphorelay signaling in response to ck. taking into account that ck metabolic enzymes are up-regulated by the hormone, these transcriptional changes may reflect a lowering of ck concentration, along with an increase of fruit sensitivity to ck occurring during late development. the veraison is a developmental stage characterized by a strengthening of the fruit sink action, that, also in olive, might be regulated by a complex hormone cross-talk and interaction, as demonstrated in other non climacteric fruits as grape berries. in spite of any signal variation in terms of expression of genes related to ethylene biosynthesis, an up-regulation of an ers type ethylene receptor has been observed during early fruit development. this might imply an increase in sensitivity of the fruitlet to ethylene, the hormone involved in the regulation of the immature fruit physiological drop that in olive occurs between fruit set and pit hardening. genes involved in the jasmonates metabolism and brassinosteroids biosynthesis were up- and down regulated, respectively. if the gene expression pattern is mirrored by a similar evolution of ja and br concentration in fruit tissues, this would indicate that the role of jas and brs in olive is different from other fruit types such as peach and grape. in fact, in peach it has been demonstrated that jas delay fruit ripening, while brs stimulate the onset of veraison in grape berry. the dual multiple contingency tests allowed the identification of go terms related to biotic and abiotic stress responses significantly and antagonistically distributed between libraries. the transcriptional profile of genes related to light, temperature and biotic stresses parallels the dynamics of day length and light intensity, both peaking at the pit hardening stage. moreover, these data may indicate that genes responding to environmental stimuli are transcriptionally regulated up to pit hardening, and post-transcriptionally up to the veraison.  as a concluding remark, the ssh technique allowed to identify a set of differentially expressed unique sequences. among these, corresponding to different key genes were further investigated by real-time pcr, pointing out a validation of the ssh results as high as 69%. the bioinformatic annotation of all gene sequences was useful to shed light on metabolic pathways and to understand specific regulatory networks. in fact, data here reported represent a significant contribution to the elucidation of transcriptional aspects related to carbohydrates, fas, secondary metabolites, transcription factors and hormones as well as response to biotic and abiotic stresses. particularly interesting are data related to hormones, pointing out the complexity of the role played by these compounds in olive fruit development and ripening. these molecular and bioinformatic data represent a first step toward both functional genomics and systems biology research for understanding the gene functions and regulatory networks in olive fruit growth, development and ripening. 
these results represent a first step toward both functional genomics and systems biology research for understanding the gene functions and regulatory networks in olive fruit growth and ripening. the bioinformatic annotation of all gene sequences was useful to shed light on metabolic pathways and transcriptional aspects related to carbohydrates, fatty acids, secondary metabolites, transcription factors and hormones as well as response to biotic and abiotic stresses throughout olive drupe development. a total of out of differentially expressed unique sequences was further investigated by real-time pcr, showing a validation of the ssh results as high as 69%. library-specific cdna repertories were annotated according to the three main vocabularies of the gene ontology: cellular component, biological process and molecular function. bioinformatic analysis pointed out a significantly different distribution of the annotated sequences for each go category, when comparing the three fruit developmental stages. differently from other fruit tree species, little is known about the physiological and molecular basis of the olive fruit development and a few sequences of genes and gene products are available for olive in public databases. leccino sampled at three different stages was used for the identification of differentially expressed genes putatively involved in main processes along fruit development. olea europaea l. is a traditional tree crop of the mediterranean basin with a worldwide economical high impact. this study deals with the identification of large sets of differentially expressed genes in developing olive fruits and the subsequent computational annotation by means of different software. the olive fruit-specific transcriptome dataset was used to query all known kegg metabolic pathways for characterizing and positioning retrieved est records. blastx analysis, go terms mapping and annotation analysis were performed using the blast2go software, a research tool designed with the main purpose of enabling go based data mining on sequence sets for which no go annotation is yet available. all sequenced clones were analyzed through blastx against non-redundant ncbi databases and about 60% of them showed similarity to known proteins.
 olea europaea l. is a traditional tree crop of the mediterranean basin with a worldwide economical high impact. differently from other fruit tree species, little is known about the physiological and molecular basis of the olive fruit development and a few sequences of genes and gene products are available for olive in public databases. this study deals with the identification of large sets of differentially expressed genes in developing olive fruits and the subsequent computational annotation by means of different software.  mrna from fruits of the cv. leccino sampled at three different stages was used for the identification of differentially expressed genes putatively involved in main processes along fruit development. four subtractive hybridization libraries were constructed: forward and reverse between stage and, and and. all sequenced clones were analyzed through blastx against non-redundant ncbi databases and about 60% of them showed similarity to known proteins. a total of out of differentially expressed unique sequences was further investigated by real-time pcr, showing a validation of the ssh results as high as 69%. library-specific cdna repertories were annotated according to the three main vocabularies of the gene ontology: cellular component, biological process and molecular function. blastx analysis, go terms mapping and annotation analysis were performed using the blast2go software, a research tool designed with the main purpose of enabling go based data mining on sequence sets for which no go annotation is yet available. bioinformatic analysis pointed out a significantly different distribution of the annotated sequences for each go category, when comparing the three fruit developmental stages. the olive fruit-specific transcriptome dataset was used to query all known kegg metabolic pathways for characterizing and positioning retrieved est records. the integration of the olive sequence datasets within the mapman platform for microarray analysis allowed the identification of specific biosynthetic pathways useful for the definition of key functional categories in time course analyses for gene groups.  the bioinformatic annotation of all gene sequences was useful to shed light on metabolic pathways and transcriptional aspects related to carbohydrates, fatty acids, secondary metabolites, transcription factors and hormones as well as response to biotic and abiotic stresses throughout olive drupe development. these results represent a first step toward both functional genomics and systems biology research for understanding the gene functions and regulatory networks in olive fruit growth and ripening. 