BACKGROUND
fruit are important components in the human diet, serving as sources of vitamins, minerals, fiber and antioxidants  <cit> .
in some parts of the world including equatorial regions of asia, africa and americas, fruits can be among the most significant source of limiting nutrients and food security.
during the final stages of development, fruit undergo the unique process of ripening which is often characterized by dramatic changes in color, texture, flavor and aroma, that are controlled by both external signals  and endogenous hormonal and genetic regulators that render the fruit organ attractive and palatable for seed-dispersing organisms  <cit> .
based on ripening physiology, apple and tomato are classified as climacteric fruit, in which the onset of ripening is accompanied by a rapid increase in respiration rate, normally coincident with elevated ethylene biosynthesis  <cit> .
ethylene is a gaseous hormone able to trigger and coordinate many physiological and response processes in higher plants, including ripening.
in climacteric fruits where the hormone typically plays major regulatory roles  <cit> , it has been shown that ethylene dependent and independent events operate in tandem to regulate overall ripening  <cit> .
the fact that apple and tomato are climacteric suggests that at least some of the regulatory aspects of ripening in both species may be conserved.
tomato has served as a primary model of ripening research due to its short life cycle, ease of transformation, well characterized germplasm http://tgrc.ucdavis.edu/ and availability of extensive molecular resources http://solgenomics.net; http://ted.bti.cornell.edu/;  http://www.pgb.kazusa.or.jp/kaftom/.
in this regard it is an optimal reference system for comparative genomics of climacteric ripening with apple.
in an effort to define comprehensive transcriptome variation with the final aim of identifying candidate ripening genes important for apple and conserved among climacteric species, we employed an expression profiling strategies using both heterologous  and homologous  expression platforms.
the het array tom <dig> was specifically constructed for functional study in tomato http://ted.bti.cornell.edu/cgi-bin/tfgd/array/home.cgi, while the hom array was dedicated for apple analysis as a component of the hidras eu-project http://www.hidras.unimi.it.
microarray technology has received enormous emphasis in recent years by the scientific community, due to its capabilities of analyzing transcription activity in a high throughput fashion  <cit> , especially in those species where large amounts of gene sequence are available  <cit> .
we chose to study the ripening process of apple performing a biological assay with an heterologous tomato array because of its large and well detailed collection of genomic information, and because of the success of a similar approach described in both vertebrate species  <cit>  and plants  <cit> .
this study aimed to improve knowledge about ripening control in apple by identifying new elements involved in this process, keeping in mind that the use of a heterologous cdna array is limited to those genes that retain a minimal degree of sequence homology.
however, a similar phenomena occurs in cdna homologous array hybridization, due to the cross hybridization of members belonging to the same gene family characterized by high sequence similarity  <cit> .
the use of tom <dig> was justified by its greater coverage, while the value of the hom apple array was grounded in the fact that this is a fruit dedicated array, already used to identify genes differentially expressed during fruit development and maturation in apple  <cit> , and in this context represented a valuable tool to confirm the heterologous data.
the use of a common reference genomic tool is an attractive prospect for analysis of non-model plants.
the two platforms were simultaneously assessed to characterize fruit ripening transcription dynamics with an emphasis on ethylene-regulated genes.
earlier high-throughput genomic efforts on fruit ripening and quality have been reported by several groups, with the first being aharoni et al.
who identified genes related to strawberry fruit quality  <cit> .
with regards to prior transcriptomics studies on tree fruit species, the peach microarray μpeach <dig>  was used to study gene expression changes associated with the transition from pre-climacteric to climacteric fruit development  <cit> , while in nectarine it was used to elucidate transcriptome variation in response to the ethylene perception inhibitor 1-mcp  <cit> .
apple has been investigated with a small number of microarray platforms, though most are of limited size, or more focused on fruit development and pre-harvest ripening  <cit> .
large scale statistical analysis of ests in apple have been reported, including an in silico comparison with tomato  <cit> .
gene expression comparisons between apple and tomato, two fleshy-fruited species belonging respectively to the rosaceae and solanaceae, could be very informative in unraveling the unique and common determinants of ripening control.
tomato has been widely used as the primary model species for climacteric fruit ripening  <cit> , and a comprehensive transcriptomic tool kit has been developed to analyze the underlying genetic ripening network.
alba et al.
<cit>  described extensive time-series expression profiling of wild-type tomato fruit using the tom <dig> array .
in this work  <cit> ,  <dig> of approximately  <dig>  genes assayed were differentially expressed during the fruit maturation process, 37% of which were altered in comparison with the never ripe ethylene receptor mutant  <cit> .
here we report our exploration of ethylene dependent and independent trascriptomics of apple fruit maturation and ripening as compared to tomato.
the apple ethylene transcriptome was further characterized in the context of the response caused by the ethylene perception competitor 1-methylcyclopropene   <cit> .
we present a comprehensive biological cross-species genomic comparison between apple and tomato, using and comparing homologous and heterologous cdna microarrays in addition to 2-dimensional protein separation, to highlight conserved and unique gene activities contributing to the complex and important mechanism of climacteric ripening control.
methods
plant material and characterization of fruit ripening physiology
mondial gala fruit were harvested at commercial ripening, and a subset of  <dig> fruits were treated overnight with 1ppm of 1-mcp at 24°c in sealed containers.
ethylene production was monitored by gas-chromatography , on five fruits, three times/week per two fruit batches, analyzed for  <dig> and  <dig> days, respectively, after harvest.
in these two time-course experiments, ethylene was measured for both control and treated sample in order to profile a normal evolution compared to the kinetics affected by the ethylene competitor .
samples were also assessed for fruit firmness with a digital firmness tester  on the two peeled and opposite fruit surfaces of each fruit tested.
transcriptome profiling was carried out with rna isolated from seven tissues  with three biological replicas for each .
three time-points spanned early fruit maturation: green , breaker , and red ripe , and four time-points spanned late ripening stages: t1ctrl  t2ctrl  for the control, and the two corresponding 1-mcp treated samples: t11-mcp and t21-mcp.
the experimental design is thus characterized by two time courses, with the first spanning climacteric ripening  and the second the same but as influenced by 1-mcp treatment .
rna isolation and het  - hom  expression profiling in apple fruit
total rna was isolated from liquid nitrogen frozen flesh  collected from all stages, using a ctab-based extraction buffer  <cit> .
after the first precipitation the procedure followed the protocol reported in the ted database http://ted.bti.cornell.edu/cgi-bin/tfgd/array/total_rna_extraction.cgi and modified  <cit> .
synthesis and cdna labelling were performed according to  <cit> .
first-strand synthesis and purification was obtained with the super-script™ indirect cdna labelling system kit , and microarray hybridization was performed labelling the cdna with cy <dig> for the reference and cy <dig> for the experimental test samples.
labelling and hybridization protocols are detailed in the ted database http://ted.bti.cornell.edu/cgi-bin/tfgd/array/total_rna_label.cgi; http://ted.bti.cornell.edu/cgi-bin/tfgd/array/tom1_hybridization.cgi
<cit> .
our experimental design employed a common reference design, hybridizing three biological replicates for each time point  <cit> .
data processing
microarray slides were processed using a two-channel confocal scanner  and the images were acquired and analyzed with scanarray v <dig>  software , setting the pmt at 65-75%, with a scanning resolution of  <dig> μm.
raw images files were captured and converted to intensity values using imagene software .
data analysis was performed using brb-array tools  <dig>  http://linus.nci.nih.gov/brb-arraytools.html, an integrated software package based on r statistic developed by dr. richard simon and amy peng lam  <cit> .
data were transformed to the log <dig> scale and normalized with the lowess methods to minimize systematic variance.
differentially expressed genes  over the time course were identified using the class comparison tool, performing a paired sample t-test .
the multiple variation tests were used with the maximum false discovery rate set at  <dig>  and 90% of confidence.
expression profile clustering was conduced with the gepas web-based resource for microarray gene expression analysis .
comprehensive heterologous and homologous data were clustered using the sota algorithm , an unsupervised neural network with a binary tree topology  <cit> .
overall graphical representation of the total heterologous transcriptome data was visualized using matlab  <dig>  .
microarray data are available in the array expression database http://www.ebi.ac.uk/microarray-as/ae/ with the number a-mexp- <dig>
quantitative real time pcr for gene specific expression profiling
from the seven samples, five μg of total rna/sample was used to reverse transcribe cdna using superscript ii reverse transcriptase  with oligo dt <dig>  real time pcr was carried out with sybr® green pcr master mix  using the following primers related to aco and pg genes: rt-md-aco1_for: caggcaacgacgcattcat, rt-md-aco1_rev: ggcgtccccagttttcttct and rt-md-pg1_for: accggtgggatagcaacatc, rt-md-pg1_rev: attccctttagctccaaaatcgt.
amplicon detection was performed using an abi prism  <dig> sequencing detection system with the following thermal profile: 95°c for 10' and subsequent  <dig> cycles of 95°c for  <dig> sec, 60°c for  <dig> min and 48°c for  <dig> sec.
amount of target was normalized to an endogenous reference  and expressed as 2-ΔΔct .
2d proteomic analysis
total protein extraction used two grams of frozen cortex collected from two samples of mondial gala apple: t1ctrl and t11-mcp.
the extraction buffer contained  <dig> mm tris-hcl ,  <dig> mm sucrose,  <dig> mm edta,  <dig> mm ascorbic acid,  <dig> mm pmsf,  <dig> % triton x- <dig>   <dig> μm leupeptin and  <dig> mm pefabloc.
sample were solubilized in  <dig> m urea,  <dig> m tiourea, 2% w/v chaps, 2% w/v triton x- <dig>  2% w/v ampholytes ipg buffer , 5% w/v dtt and a trace of bromophenol blue.
the first dimension was carried out on an ettan ipgphor i  at  <dig> kvs.
the strips were subsequently equilibrated for  <dig> min in  <dig> mm tris-hcl ,  <dig> m urea, 30% v/v glycerol, 2% w/v sds and 2% w/v dtt; and for an additional  <dig> min in  <dig> mm tris-hcl ,  <dig> m urea, 30% v/v glycerol, 2% w/v sds and  <dig> % w/v iodacetamide and a trace of bromophenol blue.
after equilibration, the strips were placed on  <dig> % polyacrylamide gels  in the ettan dalt six electrophoresis system , and the second dimension was carried out at  <dig> ma/gel over-night at 10°c.
gel images were analyzed using the image master platinum v. <dig>  software .
spot matching was performed using two synthetic gels by overlapping three gels per sample and placing  <dig> anchors in each gel.
RESULTS
transcriptome dynamics in apple fruit
transcriptome analysis was performed by hybridizing all samples comprising the experimental design  on both the tom <dig> and apple arrays containing  <dig>  and  <dig>  ests, respectively .
the total expression data set, represented by  <dig>  filtered and normalized features, was organized into two structured sota clusters, comparing functional dynamics between normal fruit ripening and that altered by 1-mcp in the following stages: green , breaker , red ripe , t1ctrl/1-mcp  and t2ctrl/1-mcp .
to highlight the functional differences between control and 1-mcp treatment, the entire expression regime was re-plotted considering only the last stages of the experimental design  focusing on the regulatory effects resulting from the ethylene response inhibitor treatment .
four clusters were isolated from the sota tree organization with two down- and two up-regulated profiles during normal ripening .
within these profiles was the transcriptomic variation caused by 1-mcp  impacting genes involved in hormone biosynthesis/response, cell wall metabolism, transcription and secondary metabolism .
similar results in terms of classes of annotated genes were also obtained employing the hom array.
in this case the general profile was represented within a hierarchical clustering , revealing that the 1-mcp impact on general ripening was similar to that detected with the het array .
gene expression dynamics were analyzed over the climacteric time course using a permutation-based paired-sample t-test performed with brb-array tools, which identified  <dig> degs  in the het data set  and  <dig> degs using the hom array .
in silico cross-species nucleotide sequence alignment was performed comparing the total set of  <dig>  annotated tomato unigenes comprising tom <dig> with the public apple est collection http://www.rosaceae.org.
general similarity analysis was performed via blastn with a cut-off value of  <dig> × 10e- <dig>  in this analysis, 52% of the genes were considered homologous, as defined by a nucleotide sequence identity greater than 75% .
out of  <dig> differentially expressed het array unigenes associated with apple ripening,  <dig> matched to a corresponding apple sequence with an identity value higher than the 75% considered necessary for designation as "homologous" .
this analysis yielded results similar to what was reported in a prior solanaceae cross species analysis, where 75% of the available sequences of pepper and eggplant were homologous with those on the tomato array  <cit> .
while apple is more distant from tomato than other members of the solanaceae, it is noteworthy that the tom <dig> array was weighted toward fruit development related genes and these sequences may be more conserved among climacteric fruits.
the comparison between tomato and arabidopsis revealed a lower identity value, further supporting the hypothesis that the tom <dig> array might be enriched for conserved fruit-associated genes and thus for using the tomato array as a reference for gene expression studies in fleshy fruits.
apple is a fleshy, climacteric, indehiscent fruit like tomato, while arabidopsis has non-fleshy, non-climacteric and dehiscent fruit.
it is noteworthy that the difference in array features between the tomato and apple arrays is correspondingly reflected mainly in the number of elements present in each gene category, with  <dig> degs grouped in  <dig> classes, while the  <dig> homologous degs in apple were grouped in  <dig> categories - only 20% less .
to highlight the genes differentially impacted by 1-mcp  we compared the profiles of the control samples with their 1-mcp treated counterparts .
based on pearson correlation  we observed 7% and 25% differential expression using the het and hom arrays, respectively.
it is interesting to note that 1-mcp treatment caused up-regulation of a gene set normally repressed in the presence of ethylene  which contains members functionally associated with nucleic acid metabolism , protein biosynthesis , signal transduction  and gene transcription .
using both arrays we identified elements representing genes up-regulated by 1-mcp treatment including iaa <dig>  ribosomial protein, pyruvate decarboxylase, a ripening-related protein of unknown function and a heat shock protein.
this negative response to ethylene was also observed in a preliminary 2-d proteomic comparison between t1ctrl/1-mcp .
specifically, in silico proteomic comparisons revealed  <dig> common peptide spots  from the total number of  <dig>  and  <dig> spots for the control and 1-mcp treated, respectively.
1-mcp treatment revealed  <dig> unique spots, corresponding to  <dig> % of the total protein pattern.
the difference in the percentage of impacted genes versus proteins likely reflects the fact that the proteomics analysis we performed is likely to over-emphasize the more abundant structural proteins and under-represent the less abundant regulatory peptides suggesting that a greater number of low abundance and putative regulatory proteins may be impacted by 1-mcp.
alternatively  this discrepancy may reflect the fact that multiple gene family members may encode proteins that cannot be distinguished via the 2-d analysis.
more comprehensive proteomics analysis that results in peptide identification would be required to sort out these or other possibilities.
the number of genes and proteins observed to be impacted by 1-mcp treatment is indeed quite significant especially considering that in this work the 1-mcp treatment occurred at harvest  when ripening is well underway.
in the hom array we observed a higher rate of genes influenced by 1-mcp, but it must be noted that while the het array is enriched for fruit-related sequences  <cit>  the hom array is a fully fruit dedicated array, and thus even more genes than in the het array would be anticipated to display differential expression during ripening and as a consequence of 1-mcp treatment.
it is also interesting that we observed an approximately equal ratio of genes either positively or negatively regulated by ethylene, demonstrating that ethylene in apple has an important and complex impact on ripening physiology, regulating both positively or negatively the expression of various genes.
even more interesting is the fact that a considerable number of ripening-related genes were not affected by the hormone, indicating ethylene-independent ripening mechanisms in climacteric apple fruit.
such genes may be especially interesting as candidates for common regulatory control between climacteric and non-climacteric fruits and thus represent a unique set of genes for further investigation.
ethylene dependent and independent genes have been previously reported in melon  <cit>  where it was suggested that in at least some cases members of the same gene family were regulated under these two distinct ripening control processes  <cit> .
transcriptional control of ethylene synthesis, perception and signalling in apple
mondial gala fruit development was characterized at three distinct physiological stages:  green ,  breaker  and  red ripe .
fruit ripening initiation was defined by the induction of ethylene starting from red ripe  and increasing at t1ctrl  and t2ctrl , producing  <dig>  μlkg-1h- <dig> and  <dig>  μlkg-1h- <dig> of ethylene respectively.
following harvest, 1-mcp application resulted in reduced ethylene synthesis with an ethylene synthesis rate of only  <dig>   and  <dig>  μlkg-1h- <dig> , or a reduction in ethylene synthesis of greater than 90% .
to further investigate the efficacy of 1-mcp ethylene repression, we performed a second ethylene assessment extended to  <dig> days after harvest.
at the end of this period the maximum ethylene production in the 1-mcp treated sample was similar to the control , but with a shift of  <dig> days.
specifically, the control sample produced its maximum  after  <dig> days following harvest , while the treated produced  <dig>  μlkg-1h- <dig>  but at day  <dig> .
these two maximum amounts of ethylene were not statistically different .
het microarray analysis identified a set of hormone-related degs whose expression was differentially coordinate by ethylene.
during fruit development, genes involved in the auxin biosynthetic pathway  showed their maximum transcript accumulation at breaker stage .
a similar elevated expression trend was also observed for genes involved in ethylene response including an ethylene inducible protein  and ethylene responsive element .
the tom <dig>  data highlighted the ethylene dependent transcriptional control which was additionally supported by 1-mcp application as reflected in down regulation of genes involved in ethylene synthesis and signal transduction which are normally highly expressed during the ethylene burst and include such genes as: sam 1-2- <dig> , acs , aco , eil , erebp  and erf <dig>  genes.
erfs are especially interesting as they have documented roles in regulating ethylene responsive genes  <cit> , and an ethylene dependent gene expression was confirmed previously through 1-mcp application in apple by wang et al.
<cit> .
identification of differentially expressed genes in the ethylene synthesis pathway and cell wall metabolism was also considered validation of the utility of the heterologous array platform.
thus, to additionally confirm het microarray validity, we assessed the expression of aco and pg, two genes involved in ethylene biosynthesis and cell wall metabolism , respectively, via qpcr in apple.
we also interrogated the ted database for tomato in silico digital expression profiles .
in both cases the expression was consistent with that observed in het and hom profiling.
specifically, positive regulation at the onset of climacteric ripening and down-regulation upon ethylene inhibition .
we note that in this comparison our main limitation was alignment of the apple and tomato developmental time courses where in tomato the ethylene burst occurs at the breaker stage , while in apple this is a post-harvest phenomena.
nevertheless, in the case of the aco and pg controls, maximum gene expression was coincident with the hormone burst and declined thereafter in both species and with both the hom and het arrays suggesting a degree of reliability in both platforms.
in the gene set negatively correlated with climacteric ethylene production  we identified hormone-related genes associated predominantly with plant growth rather than ethylene with the best hits to members involved in auxin  and gibberellin  response pathways.
furthermore, in the hom array we identified seven genes putatively involved in hormone signalling networks.
three were homologous to regulatory genes of the auxin pathway, while the other four were involved in ethylene biosynthesis and signalling: aco, ethylene receptor and md-etr genes .
the expression of both aco and md-etr was ethylene dependent, as both were up-regulated during normal ripening and repressed by 1-mcp.
transcriptional control of fruit softening in apple
to establish ties between gene expression pattern and fruit physiology, we analyzed a number of ripening parameters in the same fruit used for expression profiling including measurement of fruit firmness .
firmness in t <dig> fruit was  <dig>  kg cm- <dig> for the control and  <dig>  kg cm- <dig> for the 1-mcp treated samples and this difference increased in fruit at the t <dig> stage, with  <dig> kg cm- <dig> for t2ctrl and  <dig>  kg cm- <dig> for t21-mcp, respectively.
1-mcp treated samples lost only  <dig>  kg cm- <dig> after  <dig> days of post-harvest ripening, compared to the  <dig>  kg cm- <dig> lost by the control .
using the het array we identified  <dig> degs putatively involved in cell wall metabolism .
we observed two general trends of gene expression related to cell wall enzymes .
the first trend comprise a set of genes whose maximum expression occurs in the breaker/red ripe stages and then decreases during the post-harvest period, such as pectin acetylesterase , cellulose synthase , chitinase  and extensin .
the second category is characterized by genes whose maximum transcript abundance is observed at the end of the time course, coincident with the ethylene burst in apple.
transcripts of xyloglucan endotransglycosylases , xyloglucan endo 1- <dig> glucanase  and polygalacturonase  fall into this category.
our data regarding cell wall gene expression profiling was consistent with the results of others  <cit>  reporting both early and late fruit development enzymatic actions associated with fruit softening.
among the genes involved in later stage cell wall metabolism, polygalacturonase, xyloglucan endotransglycosylases and xyloglucan endo 1- <dig> glucanase, in particular, showed down regulation after 1-mcp treatment, confirming the importance of ethylene and cell wall metabolizing enzymes in fruit softening control.
in the hom array seven cell wall unigenes demonstrated differential expression  including pectin acetylesterase precursor, endoxyloglucan transferase, xylose isomerase and a polygalacturan gene.
according to the het profile, polygalacturonase found on the hom array showed an ethylene dependent profile as well, with a maximum release at t2ctrl, and strong down-regulation following 1-mcp treatment.
transcription factor gene expression in fruit maturation and ripening
fruit development and ripening are highly coordinated by an emerging set of transcription factors which have been defined largely in tomato though shown to have counterparts in other climacteric and non-climacteric species  <cit> .
using tom <dig>   <dig> % of the differentially expressed genes during fruit development and ripening were annotated as putative transcription factors, and 40% of these were up regulated in the preclimacteric phase , while the majority were highly expressed at the onset of ripening , as typical in tomato  <cit> .
within this category, the most common gene families were myb, ap <dig> domain, bzip, mads-box, bhlh and wizz, a set of transcription factors which are generally the most abundant in all eukaryote genomes sequenced to date  <cit> .
1-mcp application affected the expression of 13% of these transcription factor genes .
unigenes belonging to the bzip  group and wizz  were down regulated by 1-mcp treatment , while other elements including mads-box gdef <dig> , myb tmh <dig>  and ap <dig>  genes were stimulated by this treatment, suggesting negative regulation by ethylene in a subset of putative ripening regulators .
five ap <dig> members were identified in the het array, both developmentally and ripening regulated, in agreement with observations in peach  <cit> .
transcription factors that showed the highest change in mrna abundance between the breaker and red ripe stages included bell-like , ccr <dig> , jasmonic acid  <dig> , pti <dig> , yabby <dig> , myb , set domain , zpt <dig> , bzip , bhlh  and ap <dig> .
in addition to those induced during fruit development, others showed a unique profile of being primarily specific to ripening/ethylene expression such as homeodomain , bzip , dof-zinc , tiny , wrky , wizz  and nam .
transcription factors belonging to ap <dig>  bhlh, bzip, homeobox, mads, myb and the squamosa families showed complex and unique fruit development and ripening regulation, with different members of these families induced in both periods.
the role of bzip in fruit ripening has been documented in other species such as watermelon  <cit> , tomato  <cit>  and strawberry  <cit> .
in climacteric fruits the role of mads-box genes has been previously confirmed in tomato via demonstration of the necessity of lemads-rin for ripening  <cit>  and of mdpi in apple seed and fruit development  <cit> .
the relevance of this family in these physiological processes was also supported by interaction studies where  <dig> mads-box elements including tdr <dig> were shown to interact in vitro with lemads-rin  <cit> .
in the current hom apple data set only two transcription factors have been identified : a leucine rich repeats  protein, and a mads box gene further supporting the value of the het data developed here.
both genes were positively regulated by ethylene.
1-mcp application caused dramatic repression of both genes further supporting the role of ethylene in the regulation of these genes.
to have a better picture of fruit transcriptomics, the entire hom transcription profile was organized in a functional dendrogram which revealed two main clusters .
the 1-mcp treated samples clustered in a group closer to the earlier pre-ripening developmental stages .
in fact, within this deg group were identified elements typical of organs in development, involved in energy biosynthesis, nucleic acid metabolism and transcriptional control.
by this functional correlation it is evident that during fruit development and ripening the change between the breaker and ripening stage are determined in large part by differential gene expression.
together with ethylene receptors, transcription factors represent key developmental timing systems  <cit> .
a delayed ethylene burst due to 1-mcp could induce the plant to activate compensatory regulatory machinery in order to re-establish normal ripening physiology.
extending the fruit post-harvest observation up to  <dig> days after harvest we have in fact observed that the sample treated with 1-mcp reached almost the same maximum amount of ethylene suggesting such feed back control mechanisms operate during ripening.
conserved expression dynamics in apple and tomato fruit
one of our primary objectives was to add to the reservoir of ripening knowledge by identifying genes co-associated with ripening and ethylene response in both tomato and apple.
toward this end we compared the  <dig> genes differentially expressed in this work with the  <dig> tomato homologous genes reported by alba et al.
<cit> .
comparison of these two data sets , identified  <dig> genes common to both data sets.
of these,  <dig> genes had apple orthologous sequence counterparts with higher than 75% identity, representing a putative gene set of homologs involved in the ripening of both species .
a similar number of genes  differentially expressed over the course of fruit development of both apple and tomato was reported by janssen et al.
<cit> , however, only 20% of these genes were identified in this study so that a new and larger collection of highly homologous ripening-associated genes is available from tomato and apple.
such genes represent a foundation from which candidate conserved genes among other fruit species can be mined.
the major differences with the janssen report were likely due to the differences between experimental designs.
in their case the samples collected where more representative of the pre-ripening fruit development and initial maturation, while this work focused on the later development and post-harvest ripening.
as such this work focuses on the climacteric stage and the corresponding link with the ethylene production was confirmed through the finding of twelve microarray elements representing genes involved in hormone production and response .
in janssen et al.
<cit>  six genes were defined as ethylene related and all reflected sam synthase genes, one of the early enzymes in the ethylene biosynthetic pathway.
in janssen et al.
<cit>  five genes were found related to cell wall metabolism in common with tomato, of which two are α-expansin, here showed to be more active during the initial fruit softening stage.
our work identified  <dig> cell wall sequences in common with tomato though more related to later stages of fruit softening .
the combination of these two data sets improves the common genomic comparison between apple and tomato, expanding the number of possible genes commonly active in both  species from early development through post-harvest ripening.
twenty-eight genes that are differentially expressed in the ripening of apple and tomato , specifically involved in ethylene biosynthesis and response, cell wall metabolism and transcription factors.
the heterologous expression profiles of the  <dig> common genes for apple were compared with tomato in silico digital expression data retrieved from the ted database http://ted.bti.cornell.edu/.
pearson correlation confirmed that over the course of normal climacteric ripening almost 70% of the genes identified by array analysis were consistent with the digital expression results, again providing validation for the het array assay in apple and suggesting equally useful results could be recovered from other fruit species to expand the comparative fruit genomics base.
analyzing the digital expression profile of the common gene set using the wt vs. nr tomato comparison, we observed differential expression of  <dig> % of genes, consistent with what was reported by alba et al.
<cit> .
in apple, for the same gene set identified with the same array, 1-mcp resulted in differential expression of 30% of these genes, consistent with the tomato results and suggesting that these elements might have a common regulatory role in ripening control of both tomato and apple.
CONCLUSIONS
in this work we presented a heterologous approach to investigate the transcriptome of apple ripening and common genes with tomato that may serve as a base collection of candidates for conserved ripening regulation among diverse fruit species.
the het tomato array was used with the principal goal of identifying candidate genes related to fruit development and ripening that could then be related to previously characterized tomato ripening genes.
gene validation was partially gained by parallel hybridization of the same set of samples with the hom array .
heterologous microarray platforms could provide a useful alternative to explore transcriptome dynamics in the absence of a whole genome array and est data.
here we have demonstrated the use of a tomato array to identify putative apple sequences which are associated with ripening and thus may be targets for further fruit ripening and comparative genomics studies.
at present to the scientific community is presented a growing number of advanced next generation sequencing technologies  that provide a viable alternative to microarray analysis  <cit> .
despite their great potential, these new technologies still present some bias mainly related to technical features of the outputs  <cit>  and the need for strong bioinformatics support to exploit these data.
with the recent availability of the apple genome  <cit>  and the impending release of the tomato genome, these ngst will certainly open new possibilities to target gene expression with high fidelity such as the repertoire of candidate ripening and evolutionarily conserved fruit genes that can be further expanded between tomato and apple and extended to additional important fruit crop species.
authors' contributions
fc designed the experiment, collected the samples, performed the microarray hybridization, analyzed the data and wrote the manuscript.
ra supported the microarray hybridization, contributed to the data analysis and revised the manuscript.
hs provided the apple array.
vs and lg supported the interpretation and annotation of the homologous array.
sserra and sm performed the 2-d proteomic assay.
ssansavini supported the sample collection.
gc supported the 1-mcp treatment and gas-chromatography analysis for ethylene.
zf performed the bioinformatics analysis.
jg contributed to the experiment design, provided important advise during the work and improved the manuscript.
all authors read and approved the final manuscript.
supplementary material
additional file 1
comparison between the hybridization of the tomato array  and the apple array  with apple cdna.
images  and  represent the sub-grid magnification for both arrays.
click here for file
additional file 2
comparative dynamics between the control  and 1-mcp treated  sample.
the 3d plot refers to the up regulated profile.
in the boxes are highlighted specific genes of the functional profile.
click here for file
additional file 3
de genes in apple identified using the heterologous  array tom <dig>
click here for file
additional file 4
de genes in apple identified using the homologous  apple array.
click here for file
additional file 5
tom <dig> unigenes homologous to the apple est dataset collection.
homology is defined by 75% identity.
click here for file
additional file 6
go annotations for de genes identified with the het  and hom  arrays.
tables include go annotations, the number of unigenes identified, and the relative percentage of unigenes identified.
click here for file
additional file 7
comparative dynamics between the control sample  and the 1-mcp treated sample .
the profiles refer to down regulation.
in this particular case a negative regulation is reflected into an up regulation in the positive part of the plot .
click here for file
additional file 8
2d proteomic comparison profile carried out using t1ctrl and t11-mcp samples.
each synthetic gel has been obtained from  <dig> gels per sample.
colored squares represent the anchors used to facilitate the comparison.
data at the bottom of the figure summarize spot numbers and the relative matching values.
click here for file
additional file 9
hierarchical clustering of gene expression patterns identified with the het array.
three functional categories are shown: hormone pathways , transcription factors  and cell wall enzymes .
the three clusters show functional dynamics of late ripening and comparison with 1-mcp.
samples are coded as rr for red ripe: c <dig> and c <dig> for t <dig> and t <dig> control respectively; m <dig> and m <dig> for t <dig> and t <dig> 1-mcp treated respectively.
click here for file
additional file 10
expression profiles for aco and pg in apple and tomato.
shows expression profiles in developing apple fruit, as determined by qpcr.
shows digital expression profiles in tomato fruit, as retrieved from the ted database.
data for aco and pg are shown in red and black, respectively.
the solid line indicates the control samples and the dashed line indicates samples treated with 1-mcp.
abbreviations: dafb, days after full bloom; dap, days after pollination; 1-mcp, 1-methylcyclopropene.
click here for file
additional file 11
expression patterns for genes involved in hormone responses, as determined with the hom array.
the black line indicates the control samples and the red line indicates samples treated with 1-mcp.
abbreviations: 1-mcp, 1-methylcyclopropene.
click here for file
additional file 12
expression patterns for genes involved in cell wall metabolism, as determined with the hom array.
the black line indicates the control samples and the red line indicates samples treated with 1-mcp.
abbreviations: 1-mcp, 1-methylcyclopropene.
click here for file
additional file 13
expression patterns for genes encoding transcription factors, as determined with the hom array.
the black line indicates the control samples and the red line indicates samples treated with 1-mcp.
abbreviations: 1-mcp, 1-methylcyclopropene.
click here for file
additional file 14
expression clustering dendrogram with centered correlation and average linkage.
the cluster was produced using hom array data and shows the expression profile similarity among samples.
click here for file
additional file 15
eighty genes  expressed during the ripening of both apple and tomato.
click here for file
acknowledgements
this paper has been carried out with financial support from the commission of the european communities, specific research program "quality of life and management of living resources", qlrt-2001- <dig> "high-quality disease resistant apples for a sustainable agriculture".
it does not necessarily reflect its views and in no way anticipates the commission's future policy in this area.
the authors thank sara stella for qpcr assay.
the author thanks the reviewers for providing stimulating and constructive comments to the manuscript.

 fruit are important components in the human diet, serving as sources of vitamins, minerals, fiber and antioxidants. in some parts of the world including equatorial regions of asia, africa and americas, fruits can be among the most significant source of limiting nutrients and food security. during the final stages of development, fruit undergo the unique process of ripening which is often characterized by dramatic changes in color, texture, flavor and aroma, that are controlled by both external signals and endogenous hormonal and genetic regulators that render the fruit organ attractive and palatable for seed-dispersing organisms. based on ripening physiology, apple and tomato are classified as climacteric fruit, in which the onset of ripening is accompanied by a rapid increase in respiration rate, normally coincident with elevated ethylene biosynthesis. ethylene is a gaseous hormone able to trigger and coordinate many physiological and response processes in higher plants, including ripening. in climacteric fruits where the hormone typically plays major regulatory roles, it has been shown that ethylene dependent and independent events operate in tandem to regulate overall ripening. the fact that apple and tomato are climacteric suggests that at least some of the regulatory aspects of ripening in both species may be conserved. tomato has served as a primary model of ripening research due to its short life cycle, ease of transformation, well characterized germplasm http://tgrc.ucdavis.edu/ and availability of extensive molecular resources http://solgenomics.net; http://ted.bti.cornell.edu/; http://www.pgb.kazusa.or.jp/kaftom/. in this regard it is an optimal reference system for comparative genomics of climacteric ripening with apple. in an effort to define comprehensive transcriptome variation with the final aim of identifying candidate ripening genes important for apple and conserved among climacteric species, we employed an expression profiling strategies using both heterologous and homologous expression platforms. the het array tom was specifically constructed for functional study in tomato http://ted.bti.cornell.edu/cgi-bin/tfgd/array/home.cgi, while the hom array was dedicated for apple analysis as a component of the hidras eu-project http://www.hidras.unimi.it. microarray technology has received enormous emphasis in recent years by the scientific community, due to its capabilities of analyzing transcription activity in a high throughput fashion, especially in those species where large amounts of gene sequence are available. we chose to study the ripening process of apple performing a biological assay with an heterologous tomato array because of its large and well detailed collection of genomic information, and because of the success of a similar approach described in both vertebrate species and plants. this study aimed to improve knowledge about ripening control in apple by identifying new elements involved in this process, keeping in mind that the use of a heterologous cdna array is limited to those genes that retain a minimal degree of sequence homology. however, a similar phenomena occurs in cdna homologous array hybridization, due to the cross hybridization of members belonging to the same gene family characterized by high sequence similarity. the use of tom was justified by its greater coverage, while the value of the hom apple array was grounded in the fact that this is a fruit dedicated array, already used to identify genes differentially expressed during fruit development and maturation in apple, and in this context represented a valuable tool to confirm the heterologous data. the use of a common reference genomic tool is an attractive prospect for analysis of non-model plants. the two platforms were simultaneously assessed to characterize fruit ripening transcription dynamics with an emphasis on ethylene-regulated genes. earlier high-throughput genomic efforts on fruit ripening and quality have been reported by several groups, with the first being aharoni et al. who identified genes related to strawberry fruit quality. with regards to prior transcriptomics studies on tree fruit species, the peach microarray μpeach was used to study gene expression changes associated with the transition from pre-climacteric to climacteric fruit development, while in nectarine it was used to elucidate transcriptome variation in response to the ethylene perception inhibitor 1-mcp. apple has been investigated with a small number of microarray platforms, though most are of limited size, or more focused on fruit development and pre-harvest ripening. large scale statistical analysis of ests in apple have been reported, including an in silico comparison with tomato. gene expression comparisons between apple and tomato, two fleshy-fruited species belonging respectively to the rosaceae and solanaceae, could be very informative in unraveling the unique and common determinants of ripening control. tomato has been widely used as the primary model species for climacteric fruit ripening, and a comprehensive transcriptomic tool kit has been developed to analyze the underlying genetic ripening network. alba et al. described extensive time-series expression profiling of wild-type tomato fruit using the tom array. in this work, of approximately genes assayed were differentially expressed during the fruit maturation process, 37% of which were altered in comparison with the never ripe ethylene receptor mutant. here we report our exploration of ethylene dependent and independent trascriptomics of apple fruit maturation and ripening as compared to tomato. the apple ethylene transcriptome was further characterized in the context of the response caused by the ethylene perception competitor 1-methylcyclopropene. we present a comprehensive biological cross-species genomic comparison between apple and tomato, using and comparing homologous and heterologous cdna microarrays in addition to 2-dimensional protein separation, to highlight conserved and unique gene activities contributing to the complex and important mechanism of climacteric ripening control. methods plant material and characterization of fruit ripening physiology mondial gala fruit were harvested at commercial ripening, and a subset of fruits were treated overnight with 1ppm of 1-mcp at 24°c in sealed containers. ethylene production was monitored by gas-chromatography, on five fruits, three times/week per two fruit batches, analyzed for and days, respectively, after harvest. in these two time-course experiments, ethylene was measured for both control and treated sample in order to profile a normal evolution compared to the kinetics affected by the ethylene competitor. samples were also assessed for fruit firmness with a digital firmness tester on the two peeled and opposite fruit surfaces of each fruit tested. transcriptome profiling was carried out with rna isolated from seven tissues with three biological replicas for each. three time-points spanned early fruit maturation: green, breaker, and red ripe, and four time-points spanned late ripening stages: t1ctrl t2ctrl for the control, and the two corresponding 1-mcp treated samples: t11-mcp and t21-mcp. the experimental design is thus characterized by two time courses, with the first spanning climacteric ripening and the second the same but as influenced by 1-mcp treatment. rna isolation and het - hom expression profiling in apple fruit total rna was isolated from liquid nitrogen frozen flesh collected from all stages, using a ctab-based extraction buffer. after the first precipitation the procedure followed the protocol reported in the ted database http://ted.bti.cornell.edu/cgi-bin/tfgd/array/total_rna_extraction.cgi and modified. synthesis and cdna labelling were performed according to. first-strand synthesis and purification was obtained with the super-script™ indirect cdna labelling system kit, and microarray hybridization was performed labelling the cdna with cy for the reference and cy for the experimental test samples. labelling and hybridization protocols are detailed in the ted database http://ted.bti.cornell.edu/cgi-bin/tfgd/array/total_rna_label.cgi; http://ted.bti.cornell.edu/cgi-bin/tfgd/array/tom1_hybridization.cgi  . our experimental design employed a common reference design, hybridizing three biological replicates for each time point. data processing microarray slides were processed using a two-channel confocal scanner and the images were acquired and analyzed with scanarray v software, setting the pmt at 65-75%, with a scanning resolution of μm. raw images files were captured and converted to intensity values using imagene software. data analysis was performed using brb-array tools http://linus.nci.nih.gov/brb-arraytools.html, an integrated software package based on r statistic developed by dr. richard simon and amy peng lam. data were transformed to the log scale and normalized with the lowess methods to minimize systematic variance. differentially expressed genes over the time course were identified using the class comparison tool, performing a paired sample t-test. the multiple variation tests were used with the maximum false discovery rate set at and 90% of confidence. expression profile clustering was conduced with the gepas web-based resource for microarray gene expression analysis. comprehensive heterologous and homologous data were clustered using the sota algorithm, an unsupervised neural network with a binary tree topology. overall graphical representation of the total heterologous transcriptome data was visualized using matlab. microarray data are available in the array expression database http://www.ebi.ac.uk/microarray-as/ae/ with the number a-mexp- quantitative real time pcr for gene specific expression profiling from the seven samples, five μg of total rna/sample was used to reverse transcribe cdna using superscript ii reverse transcriptase with oligo dt real time pcr was carried out with sybr® green pcr master mix using the following primers related to aco and pg genes: rt-md-aco1_for: caggcaacgacgcattcat, rt-md-aco1_rev: ggcgtccccagttttcttct and rt-md-pg1_for: accggtgggatagcaacatc, rt-md-pg1_rev: attccctttagctccaaaatcgt. amplicon detection was performed using an abi prism sequencing detection system with the following thermal profile: 95°c for 10' and subsequent cycles of 95°c for sec, 60°c for min and 48°c for sec. amount of target was normalized to an endogenous reference and expressed as 2-ΔΔct. 2d proteomic analysis total protein extraction used two grams of frozen cortex collected from two samples of mondial gala apple: t1ctrl and t11-mcp. the extraction buffer contained mm tris-hcl, mm sucrose, mm edta, mm ascorbic acid, mm pmsf, % triton x- μm leupeptin and mm pefabloc. sample were solubilized in m urea, m tiourea, 2% w/v chaps, 2% w/v triton x- 2% w/v ampholytes ipg buffer, 5% w/v dtt and a trace of bromophenol blue. the first dimension was carried out on an ettan ipgphor i at kvs. the strips were subsequently equilibrated for min in mm tris-hcl, m urea, 30% v/v glycerol, 2% w/v sds and 2% w/v dtt; and for an additional min in mm tris-hcl, m urea, 30% v/v glycerol, 2% w/v sds and % w/v iodacetamide and a trace of bromophenol blue. after equilibration, the strips were placed on % polyacrylamide gels in the ettan dalt six electrophoresis system, and the second dimension was carried out at ma/gel over-night at 10°c. gel images were analyzed using the image master platinum v. software. spot matching was performed using two synthetic gels by overlapping three gels per sample and placing anchors in each gel.  transcriptome dynamics in apple fruit transcriptome analysis was performed by hybridizing all samples comprising the experimental design on both the tom and apple arrays containing and ests, respectively. the total expression data set, represented by filtered and normalized features, was organized into two structured sota clusters, comparing functional dynamics between normal fruit ripening and that altered by 1-mcp in the following stages: green, breaker, red ripe, t1ctrl/1-mcp and t2ctrl/1-mcp. to highlight the functional differences between control and 1-mcp treatment, the entire expression regime was re-plotted considering only the last stages of the experimental design focusing on the regulatory effects resulting from the ethylene response inhibitor treatment. four clusters were isolated from the sota tree organization with two down- and two up-regulated profiles during normal ripening. within these profiles was the transcriptomic variation caused by 1-mcp impacting genes involved in hormone biosynthesis/response, cell wall metabolism, transcription and secondary metabolism. similar results in terms of classes of annotated genes were also obtained employing the hom array. in this case the general profile was represented within a hierarchical clustering, revealing that the 1-mcp impact on general ripening was similar to that detected with the het array. gene expression dynamics were analyzed over the climacteric time course using a permutation-based paired-sample t-test performed with brb-array tools, which identified degs in the het data set and degs using the hom array. in silico cross-species nucleotide sequence alignment was performed comparing the total set of annotated tomato unigenes comprising tom with the public apple est collection http://www.rosaceae.org. general similarity analysis was performed via blastn with a cut-off value of × 10e- in this analysis, 52% of the genes were considered homologous, as defined by a nucleotide sequence identity greater than 75%. out of differentially expressed het array unigenes associated with apple ripening, matched to a corresponding apple sequence with an identity value higher than the 75% considered necessary for designation as "homologous". this analysis yielded results similar to what was reported in a prior solanaceae cross species analysis, where 75% of the available sequences of pepper and eggplant were homologous with those on the tomato array. while apple is more distant from tomato than other members of the solanaceae, it is noteworthy that the tom array was weighted toward fruit development related genes and these sequences may be more conserved among climacteric fruits. the comparison between tomato and arabidopsis revealed a lower identity value, further supporting the hypothesis that the tom array might be enriched for conserved fruit-associated genes and thus for using the tomato array as a reference for gene expression studies in fleshy fruits. apple is a fleshy, climacteric, indehiscent fruit like tomato, while arabidopsis has non-fleshy, non-climacteric and dehiscent fruit. it is noteworthy that the difference in array features between the tomato and apple arrays is correspondingly reflected mainly in the number of elements present in each gene category, with degs grouped in classes, while the homologous degs in apple were grouped in categories - only 20% less. to highlight the genes differentially impacted by 1-mcp we compared the profiles of the control samples with their 1-mcp treated counterparts. based on pearson correlation we observed 7% and 25% differential expression using the het and hom arrays, respectively. it is interesting to note that 1-mcp treatment caused up-regulation of a gene set normally repressed in the presence of ethylene which contains members functionally associated with nucleic acid metabolism, protein biosynthesis, signal transduction and gene transcription. using both arrays we identified elements representing genes up-regulated by 1-mcp treatment including iaa ribosomial protein, pyruvate decarboxylase, a ripening-related protein of unknown function and a heat shock protein. this negative response to ethylene was also observed in a preliminary 2-d proteomic comparison between t1ctrl/1-mcp. specifically, in silico proteomic comparisons revealed common peptide spots from the total number of and spots for the control and 1-mcp treated, respectively. 1-mcp treatment revealed unique spots, corresponding to % of the total protein pattern. the difference in the percentage of impacted genes versus proteins likely reflects the fact that the proteomics analysis we performed is likely to over-emphasize the more abundant structural proteins and under-represent the less abundant regulatory peptides suggesting that a greater number of low abundance and putative regulatory proteins may be impacted by 1-mcp. alternatively this discrepancy may reflect the fact that multiple gene family members may encode proteins that cannot be distinguished via the 2-d analysis. more comprehensive proteomics analysis that results in peptide identification would be required to sort out these or other possibilities. the number of genes and proteins observed to be impacted by 1-mcp treatment is indeed quite significant especially considering that in this work the 1-mcp treatment occurred at harvest when ripening is well underway. in the hom array we observed a higher rate of genes influenced by 1-mcp, but it must be noted that while the het array is enriched for fruit-related sequences the hom array is a fully fruit dedicated array, and thus even more genes than in the het array would be anticipated to display differential expression during ripening and as a consequence of 1-mcp treatment. it is also interesting that we observed an approximately equal ratio of genes either positively or negatively regulated by ethylene, demonstrating that ethylene in apple has an important and complex impact on ripening physiology, regulating both positively or negatively the expression of various genes. even more interesting is the fact that a considerable number of ripening-related genes were not affected by the hormone, indicating ethylene-independent ripening mechanisms in climacteric apple fruit. such genes may be especially interesting as candidates for common regulatory control between climacteric and non-climacteric fruits and thus represent a unique set of genes for further investigation. ethylene dependent and independent genes have been previously reported in melon where it was suggested that in at least some cases members of the same gene family were regulated under these two distinct ripening control processes. transcriptional control of ethylene synthesis, perception and signalling in apple mondial gala fruit development was characterized at three distinct physiological stages: green, breaker and red ripe. fruit ripening initiation was defined by the induction of ethylene starting from red ripe and increasing at t1ctrl and t2ctrl, producing μlkg-1h- and μlkg-1h- of ethylene respectively. following harvest, 1-mcp application resulted in reduced ethylene synthesis with an ethylene synthesis rate of only and μlkg-1h-, or a reduction in ethylene synthesis of greater than 90%. to further investigate the efficacy of 1-mcp ethylene repression, we performed a second ethylene assessment extended to days after harvest. at the end of this period the maximum ethylene production in the 1-mcp treated sample was similar to the control, but with a shift of days. specifically, the control sample produced its maximum after days following harvest, while the treated produced μlkg-1h- but at day. these two maximum amounts of ethylene were not statistically different. het microarray analysis identified a set of hormone-related degs whose expression was differentially coordinate by ethylene. during fruit development, genes involved in the auxin biosynthetic pathway showed their maximum transcript accumulation at breaker stage. a similar elevated expression trend was also observed for genes involved in ethylene response including an ethylene inducible protein and ethylene responsive element. the tom data highlighted the ethylene dependent transcriptional control which was additionally supported by 1-mcp application as reflected in down regulation of genes involved in ethylene synthesis and signal transduction which are normally highly expressed during the ethylene burst and include such genes as: sam 1-2-, acs, aco, eil, erebp and erf genes. erfs are especially interesting as they have documented roles in regulating ethylene responsive genes, and an ethylene dependent gene expression was confirmed previously through 1-mcp application in apple by wang et al.. identification of differentially expressed genes in the ethylene synthesis pathway and cell wall metabolism was also considered validation of the utility of the heterologous array platform. thus, to additionally confirm het microarray validity, we assessed the expression of aco and pg, two genes involved in ethylene biosynthesis and cell wall metabolism, respectively, via qpcr in apple. we also interrogated the ted database for tomato in silico digital expression profiles. in both cases the expression was consistent with that observed in het and hom profiling. specifically, positive regulation at the onset of climacteric ripening and down-regulation upon ethylene inhibition. we note that in this comparison our main limitation was alignment of the apple and tomato developmental time courses where in tomato the ethylene burst occurs at the breaker stage, while in apple this is a post-harvest phenomena. nevertheless, in the case of the aco and pg controls, maximum gene expression was coincident with the hormone burst and declined thereafter in both species and with both the hom and het arrays suggesting a degree of reliability in both platforms. in the gene set negatively correlated with climacteric ethylene production we identified hormone-related genes associated predominantly with plant growth rather than ethylene with the best hits to members involved in auxin and gibberellin response pathways. furthermore, in the hom array we identified seven genes putatively involved in hormone signalling networks. three were homologous to regulatory genes of the auxin pathway, while the other four were involved in ethylene biosynthesis and signalling: aco, ethylene receptor and md-etr genes. the expression of both aco and md-etr was ethylene dependent, as both were up-regulated during normal ripening and repressed by 1-mcp. transcriptional control of fruit softening in apple to establish ties between gene expression pattern and fruit physiology, we analyzed a number of ripening parameters in the same fruit used for expression profiling including measurement of fruit firmness. firmness in t fruit was kg cm- for the control and kg cm- for the 1-mcp treated samples and this difference increased in fruit at the t stage, with kg cm- for t2ctrl and kg cm- for t21-mcp, respectively. 1-mcp treated samples lost only kg cm- after days of post-harvest ripening, compared to the kg cm- lost by the control. using the het array we identified degs putatively involved in cell wall metabolism. we observed two general trends of gene expression related to cell wall enzymes. the first trend comprise a set of genes whose maximum expression occurs in the breaker/red ripe stages and then decreases during the post-harvest period, such as pectin acetylesterase, cellulose synthase, chitinase and extensin. the second category is characterized by genes whose maximum transcript abundance is observed at the end of the time course, coincident with the ethylene burst in apple. transcripts of xyloglucan endotransglycosylases, xyloglucan endo 1- glucanase and polygalacturonase fall into this category. our data regarding cell wall gene expression profiling was consistent with the results of others reporting both early and late fruit development enzymatic actions associated with fruit softening. among the genes involved in later stage cell wall metabolism, polygalacturonase, xyloglucan endotransglycosylases and xyloglucan endo 1- glucanase, in particular, showed down regulation after 1-mcp treatment, confirming the importance of ethylene and cell wall metabolizing enzymes in fruit softening control. in the hom array seven cell wall unigenes demonstrated differential expression including pectin acetylesterase precursor, endoxyloglucan transferase, xylose isomerase and a polygalacturan gene. according to the het profile, polygalacturonase found on the hom array showed an ethylene dependent profile as well, with a maximum release at t2ctrl, and strong down-regulation following 1-mcp treatment. transcription factor gene expression in fruit maturation and ripening fruit development and ripening are highly coordinated by an emerging set of transcription factors which have been defined largely in tomato though shown to have counterparts in other climacteric and non-climacteric species. using tom % of the differentially expressed genes during fruit development and ripening were annotated as putative transcription factors, and 40% of these were up regulated in the preclimacteric phase, while the majority were highly expressed at the onset of ripening, as typical in tomato. within this category, the most common gene families were myb, ap domain, bzip, mads-box, bhlh and wizz, a set of transcription factors which are generally the most abundant in all eukaryote genomes sequenced to date. 1-mcp application affected the expression of 13% of these transcription factor genes. unigenes belonging to the bzip group and wizz were down regulated by 1-mcp treatment, while other elements including mads-box gdef, myb tmh and ap genes were stimulated by this treatment, suggesting negative regulation by ethylene in a subset of putative ripening regulators. five ap members were identified in the het array, both developmentally and ripening regulated, in agreement with observations in peach. transcription factors that showed the highest change in mrna abundance between the breaker and red ripe stages included bell-like, ccr, jasmonic acid, pti, yabby, myb, set domain, zpt, bzip, bhlh and ap. in addition to those induced during fruit development, others showed a unique profile of being primarily specific to ripening/ethylene expression such as homeodomain, bzip, dof-zinc, tiny, wrky, wizz and nam. transcription factors belonging to ap bhlh, bzip, homeobox, mads, myb and the squamosa families showed complex and unique fruit development and ripening regulation, with different members of these families induced in both periods. the role of bzip in fruit ripening has been documented in other species such as watermelon, tomato and strawberry. in climacteric fruits the role of mads-box genes has been previously confirmed in tomato via demonstration of the necessity of lemads-rin for ripening and of mdpi in apple seed and fruit development. the relevance of this family in these physiological processes was also supported by interaction studies where mads-box elements including tdr were shown to interact in vitro with lemads-rin. in the current hom apple data set only two transcription factors have been identified: a leucine rich repeats protein, and a mads box gene further supporting the value of the het data developed here. both genes were positively regulated by ethylene. 1-mcp application caused dramatic repression of both genes further supporting the role of ethylene in the regulation of these genes. to have a better picture of fruit transcriptomics, the entire hom transcription profile was organized in a functional dendrogram which revealed two main clusters. the 1-mcp treated samples clustered in a group closer to the earlier pre-ripening developmental stages. in fact, within this deg group were identified elements typical of organs in development, involved in energy biosynthesis, nucleic acid metabolism and transcriptional control. by this functional correlation it is evident that during fruit development and ripening the change between the breaker and ripening stage are determined in large part by differential gene expression. together with ethylene receptors, transcription factors represent key developmental timing systems. a delayed ethylene burst due to 1-mcp could induce the plant to activate compensatory regulatory machinery in order to re-establish normal ripening physiology. extending the fruit post-harvest observation up to days after harvest we have in fact observed that the sample treated with 1-mcp reached almost the same maximum amount of ethylene suggesting such feed back control mechanisms operate during ripening. conserved expression dynamics in apple and tomato fruit one of our primary objectives was to add to the reservoir of ripening knowledge by identifying genes co-associated with ripening and ethylene response in both tomato and apple. toward this end we compared the genes differentially expressed in this work with the tomato homologous genes reported by alba et al.. comparison of these two data sets, identified genes common to both data sets. of these, genes had apple orthologous sequence counterparts with higher than 75% identity, representing a putative gene set of homologs involved in the ripening of both species. a similar number of genes differentially expressed over the course of fruit development of both apple and tomato was reported by janssen et al., however, only 20% of these genes were identified in this study so that a new and larger collection of highly homologous ripening-associated genes is available from tomato and apple. such genes represent a foundation from which candidate conserved genes among other fruit species can be mined. the major differences with the janssen report were likely due to the differences between experimental designs. in their case the samples collected where more representative of the pre-ripening fruit development and initial maturation, while this work focused on the later development and post-harvest ripening. as such this work focuses on the climacteric stage and the corresponding link with the ethylene production was confirmed through the finding of twelve microarray elements representing genes involved in hormone production and response. in janssen et al. six genes were defined as ethylene related and all reflected sam synthase genes, one of the early enzymes in the ethylene biosynthetic pathway. in janssen et al. five genes were found related to cell wall metabolism in common with tomato, of which two are α-expansin, here showed to be more active during the initial fruit softening stage. our work identified cell wall sequences in common with tomato though more related to later stages of fruit softening. the combination of these two data sets improves the common genomic comparison between apple and tomato, expanding the number of possible genes commonly active in both species from early development through post-harvest ripening. twenty-eight genes that are differentially expressed in the ripening of apple and tomato, specifically involved in ethylene biosynthesis and response, cell wall metabolism and transcription factors. the heterologous expression profiles of the common genes for apple were compared with tomato in silico digital expression data retrieved from the ted database http://ted.bti.cornell.edu/. pearson correlation confirmed that over the course of normal climacteric ripening almost 70% of the genes identified by array analysis were consistent with the digital expression results, again providing validation for the het array assay in apple and suggesting equally useful results could be recovered from other fruit species to expand the comparative fruit genomics base. analyzing the digital expression profile of the common gene set using the wt vs. nr tomato comparison, we observed differential expression of % of genes, consistent with what was reported by alba et al.. in apple, for the same gene set identified with the same array, 1-mcp resulted in differential expression of 30% of these genes, consistent with the tomato results and suggesting that these elements might have a common regulatory role in ripening control of both tomato and apple.  in this work we presented a heterologous approach to investigate the transcriptome of apple ripening and common genes with tomato that may serve as a base collection of candidates for conserved ripening regulation among diverse fruit species. the het tomato array was used with the principal goal of identifying candidate genes related to fruit development and ripening that could then be related to previously characterized tomato ripening genes. gene validation was partially gained by parallel hybridization of the same set of samples with the hom array. heterologous microarray platforms could provide a useful alternative to explore transcriptome dynamics in the absence of a whole genome array and est data. here we have demonstrated the use of a tomato array to identify putative apple sequences which are associated with ripening and thus may be targets for further fruit ripening and comparative genomics studies. at present to the scientific community is presented a growing number of advanced next generation sequencing technologies that provide a viable alternative to microarray analysis. despite their great potential, these new technologies still present some bias mainly related to technical features of the outputs and the need for strong bioinformatics support to exploit these data. with the recent availability of the apple genome and the impending release of the tomato genome, these ngst will certainly open new possibilities to target gene expression with high fidelity such as the repertoire of candidate ripening and evolutionarily conserved fruit genes that can be further expanded between tomato and apple and extended to additional important fruit crop species. authors' contributions fc designed the experiment, collected the samples, performed the microarray hybridization, analyzed the data and wrote the manuscript. ra supported the microarray hybridization, contributed to the data analysis and revised the manuscript. hs provided the apple array. vs and lg supported the interpretation and annotation of the homologous array. sserra and sm performed the 2-d proteomic assay. ssansavini supported the sample collection. gc supported the 1-mcp treatment and gas-chromatography analysis for ethylene. zf performed the bioinformatics analysis. jg contributed to the experiment design, provided important advise during the work and improved the manuscript. all authors read and approved the final manuscript. supplementary material additional file 1 comparison between the hybridization of the tomato array and the apple array with apple cdna. images and represent the sub-grid magnification for both arrays. click here for file additional file 2 comparative dynamics between the control and 1-mcp treated sample. the 3d plot refers to the up regulated profile. in the boxes are highlighted specific genes of the functional profile. click here for file additional file 3 de genes in apple identified using the heterologous array tom click here for file additional file 4 de genes in apple identified using the homologous apple array. click here for file additional file 5 tom unigenes homologous to the apple est dataset collection. homology is defined by 75% identity. click here for file additional file 6 go annotations for de genes identified with the het and hom arrays. tables include go annotations, the number of unigenes identified, and the relative percentage of unigenes identified. click here for file additional file 7 comparative dynamics between the control sample and the 1-mcp treated sample. the profiles refer to down regulation. in this particular case a negative regulation is reflected into an up regulation in the positive part of the plot. click here for file additional file 8 2d proteomic comparison profile carried out using t1ctrl and t11-mcp samples. each synthetic gel has been obtained from gels per sample. colored squares represent the anchors used to facilitate the comparison. data at the bottom of the figure summarize spot numbers and the relative matching values. click here for file additional file 9 hierarchical clustering of gene expression patterns identified with the het array. three functional categories are shown: hormone pathways, transcription factors and cell wall enzymes. the three clusters show functional dynamics of late ripening and comparison with 1-mcp. samples are coded as rr for red ripe: c and c for t and t control respectively; m and m for t and t 1-mcp treated respectively. click here for file additional file 10 expression profiles for aco and pg in apple and tomato. shows expression profiles in developing apple fruit, as determined by qpcr. shows digital expression profiles in tomato fruit, as retrieved from the ted database. data for aco and pg are shown in red and black, respectively. the solid line indicates the control samples and the dashed line indicates samples treated with 1-mcp. abbreviations: dafb, days after full bloom; dap, days after pollination; 1-mcp, 1-methylcyclopropene. click here for file additional file 11 expression patterns for genes involved in hormone responses, as determined with the hom array. the black line indicates the control samples and the red line indicates samples treated with 1-mcp. abbreviations: 1-mcp, 1-methylcyclopropene. click here for file additional file 12 expression patterns for genes involved in cell wall metabolism, as determined with the hom array. the black line indicates the control samples and the red line indicates samples treated with 1-mcp. abbreviations: 1-mcp, 1-methylcyclopropene. click here for file additional file 13 expression patterns for genes encoding transcription factors, as determined with the hom array. the black line indicates the control samples and the red line indicates samples treated with 1-mcp. abbreviations: 1-mcp, 1-methylcyclopropene. click here for file additional file 14 expression clustering dendrogram with centered correlation and average linkage. the cluster was produced using hom array data and shows the expression profile similarity among samples. click here for file additional file 15 eighty genes expressed during the ripening of both apple and tomato. click here for file acknowledgements this paper has been carried out with financial support from the commission of the european communities, specific research program "quality of life and management of living resources", qlrt-2001- "high-quality disease resistant apples for a sustainable agriculture". it does not necessarily reflect its views and in no way anticipates the commission's future policy in this area. the authors thank sara stella for qpcr assay. the author thanks the reviewers for providing stimulating and constructive comments to the manuscript.
fruit development, maturation and ripening consists of a complex series of biochemical and physiological changes that in climacteric fruits, including apple and tomato, are coordinated by the gaseous hormone ethylene. to gain insight into the molecular process regulating ripening in apple, and to compare to tomato, we utilized both homologous and heterologous microarray to profile transcriptome dynamics of genes involved in fruit development and ripening, emphasizing those which are ethylene regulated. to date many reports have been made on the analysis of gene expression in apple. the resulting transcriptomics data revealed coordinated gene expression during fruit ripening of a subset of ripening-related and ethylene responsive genes, further facilitating the analysis of ethylene response during fruit maturation and ripening. comparison with tomato fruit maturation and ethylene responsive transcriptome activity facilitated identification of putative conserved orthologous ripening-related genes, which serve as an initial set of candidates for assessing conservation of gene activity across genomes of fruit bearing plant species. the cross-species analysis had the secondary aim of examining the efficiency of heterologous microarray hybridization for candidate gene identification as related to the ripening process. our combined strategy based on microarray hybridization enabled transcriptome characterization during normal climacteric apple ripening, as well as definition of ethylene-dependent transcriptome changes.
 fruit development, maturation and ripening consists of a complex series of biochemical and physiological changes that in climacteric fruits, including apple and tomato, are coordinated by the gaseous hormone ethylene. these changes lead to final fruit quality and understanding of the functional machinery underlying these processes is of both biological and practical importance. to date many reports have been made on the analysis of gene expression in apple. in this study we focused our investigation on the role of ethylene during apple maturation, specifically comparing transcriptomics of normal ripening with changes resulting from application of the hormone receptor competitor 1-methylcyclopropene.  to gain insight into the molecular process regulating ripening in apple, and to compare to tomato, we utilized both homologous and heterologous microarray to profile transcriptome dynamics of genes involved in fruit development and ripening, emphasizing those which are ethylene regulated. the use of both types of microarrays facilitated transcriptome comparison between apple and tomato and highlighted genes conserved during ripening of both species, which in turn represent a foundation for further comparative genomic studies. the cross-species analysis had the secondary aim of examining the efficiency of heterologous microarray hybridization for candidate gene identification as related to the ripening process. the resulting transcriptomics data revealed coordinated gene expression during fruit ripening of a subset of ripening-related and ethylene responsive genes, further facilitating the analysis of ethylene response during fruit maturation and ripening.  our combined strategy based on microarray hybridization enabled transcriptome characterization during normal climacteric apple ripening, as well as definition of ethylene-dependent transcriptome changes. comparison with tomato fruit maturation and ethylene responsive transcriptome activity facilitated identification of putative conserved orthologous ripening-related genes, which serve as an initial set of candidates for assessing conservation of gene activity across genomes of fruit bearing plant species. 