- Research article
- Open Access
Genome-wide transcriptional analysis of grapevine berry ripening reveals a set of genes similarly modulated during three seasons and the occurrence of an oxidative burst at vèraison
BMC Genomicsvolume 8, Article number: 428 (2007)
Grapevine (Vitis species) is among the most important fruit crops in terms of cultivated area and economic impact. Despite this relevance, little is known about the transcriptional changes and the regulatory circuits underlying the biochemical and physical changes occurring during berry development.
Fruit ripening in the non-climacteric crop species Vitis vinifera L. has been investigated at the transcriptional level by the use of the Affymetrix Vitis GeneChip® which contains approximately 14,500 unigenes. Gene expression data obtained from berries sampled before and after véraison in three growing years, were analyzed to identify genes specifically involved in fruit ripening and to investigate seasonal influences on the process. From these analyses a core set of 1477 genes was found which was similarly modulated in all seasons. We were able to separate ripening specific isoforms within gene families and to identify ripening related genes which appeared strongly regulated also by the seasonal weather conditions. Transcripts annotation by Gene Ontology vocabulary revealed five overrepresented functional categories of which cell wall organization and biogenesis, carbohydrate and secondary metabolisms and stress response were specifically induced during the ripening phase, while photosynthesis was strongly repressed. About 19% of the core gene set was characterized by genes involved in regulatory processes, such as transcription factors and transcripts related to hormonal metabolism and signal transduction. Auxin, ethylene and light emerged as the main stimuli influencing berry development. In addition, an oxidative burst, previously not detected in grapevine, characterized by rapid accumulation of H2O2 starting from véraison and by the modulation of many ROS scavenging enzymes, was observed.
The time-course gene expression analysis of grapevine berry development has identified the occurrence of two well distinct phases along the process. The pre-véraison phase represents a reprogramming stage of the cellular metabolism, characterized by the expression of numerous genes involved in hormonal signalling and transcriptional regulation. The post-véraison phase is characterized by the onset of a ripening-specialized metabolism responsible for the phenotypic traits of the ripe berry. Between the two phases, at véraison, an oxidative burst and the concurrent modulation of the anti-oxidative enzymatic network was observed. The large number of regulatory genes we have identified represents a powerful new resource for dissecting the mechanisms of fruit ripening control in non-climacteric plants.
Grape is among the most ancient, widely cultivated fruit crops. The keen interest in the understanding of grape berry ripening is justified by the economic relevance of the quality of grapes and their processed products, such as wine, juice and dried fruit. The onset of the genomic era has brought unprecedented progress into our knowledge of Vitis molecular biology, and powerful genetic and genomic tools are now available or are being developed, such as mapping populations, genetic and physical maps [1, 2], an extensive ESTs collection (342,576 sequences deposited at the NCBI EST database-March 19, 2007, ) and several tools for gene expression analysis. These resources enable us to identify grape as a model plant for studies on non-climacteric anthocyanin accumulating fruits, beside tomato, which is currently the model plant for climacteric carotenoid accumulating fruits. Grape berry development from anthesis to ripening is classically divided into three phases on the basis of chemical and morphological traits : berry formation, characterized by exponential growth of the berry and accumulation of organic acids, mainly malate, in the vacuole; véraison, a transition phase during which growth declines and berries start to change colour and soften; ripening, characterized by an increase in pH, marked berry growth and accumulation of sugars, anthocyanins and flavour-enhancing compounds. Published studies based on EST analysis [5–7], differential display , cDNA-arrays  and oligo arrays  have shown distinct and extensive changes in the transcriptome occurring during ripening. These studies report up to 740 developmentally regulated genes mainly involved in defence, stress response, primary and secondary metabolism, berry growth and hormonal metabolism. The very recent study of Grimplet et al. , provides an extensive description of these functional classes investigating their tissue-specific expression in the ripe berry by means of the Affymetrix Vitis GeneChip®.
With the aim of understanding the transcriptional changes underlying grapevine berry ripening and the regulatory mechanisms acting at this level, we were able to identify 1477 genes modulated during berry development, comprising 282 regulatory factors in addition to the structural genes, by the use of the Affymetrix Vitis GeneChip®. The results, obtained from replicated experiments carried out in three growing seasons, provide also insights into the relationship between genetic and environmental effects on berry development control.
Results and Discussion
Grape berries phenotypic evaluations
Figure 1 shows the changes in average berry weight, organic acids, total soluble solids and anthocyanin content of grape berries sampled weekly during the 2003, 2005 and 2006 growing seasons. The date of véraison was set at approximately 8 weeks post-flowering (wpf) in all years. In order to characterize berry ripening, we chose three distinct moments according to : the last stage of berry formation, time-point A (TP A, stage E-L 33), characterized by small hard green berries still accumulating organic acids; and two stages during berry ripening, TP B (stage E-L 34), just before véraison, characterized by berries in the green state with maximum acidic content starting to soften and TP C (stage E-L 36), when the ripening process is well established and berries are growing fast, colouring and accumulating sugars. The identification of the corresponding developmental stages in the three years was based on the biochemical profiles. As average temperatures during 2003 remained constantly unusually higher from April to October [see Additional file 1], flowering and harvest occurred approximately one week earlier in 2003.
RNA extracted in triplicate from berries harvested at three time points was used to hybridize 27 Affymetrix Vitis GeneChips®. Datasets from each season were analyzed separately. The quality of the biological replicates was evaluated by the R-squared coefficient [see Additional file 2], which ranged from 0.89 to 0.99. Minimum R-squared values were found for TP B samples, suggesting a higher sample heterogeneity approaching véraison, when dramatic transcriptional changes occur in less than 48 hours .
Principal Component Analysis (PCA) confirmed the uniformity of the biological replicates as the nine groups of replicates clusterised tightly (Fig. 2). The greatest variance in gene expression was found between samples from the three growing seasons, as they were separated along the first component (51% of the total variance). In particular, the largest difference was evident between the 2003 expression data compared to the 2005 and 2006, suggesting that the environmental conditions greatly affect the berry transcriptome. [see Additional file 1]. On the other hand, the second principal component (24% of the total variance) separated neatly the pre-véraison stages (TP A and TP B) from the post-véraison one (TP C), confirming previous observations of extensive transcriptional changes occurring from berry formation to berry ripening .
Class comparison analysis by means of Significance Analysis of Microarray , imposing a minimum fold change of 2 in at least one comparison (TP A vs. TP B and TP B vs. TP C), produced three sets of 2632, 4238 and 3856 differentially expressed probesets for the 2003, 2005 and 2006 samples, respectively. The size difference between the three sets likely reflects inter-seasonal biological differences. The comparison of the three datasets identified a common set of 1700 modulated transcripts and among these, 1477 transcripts with a conserved profile in the three years (Pearson's correlation coefficient >0.5 in any pairwise comparison). This 1477 core dataset, which corresponds to 9% of the chip probesets, is listed in Additional file 3. It represents a highly conserved expression network reasonably involved in an internal basic functional program, such as berry ripening.
Annotation of the differentially expressed genes
Automatic annotation of all the Vitis GeneChip® sequences was performed using the Gene Ontology (GO) classification . The core set of 1477 differentially expressed genes was subsequently manually checked and integrated with additional GO 'biological process' terms. Transcripts were then grouped into 17 GO functional categories [see Additional file 3 and 4] whose distribution among the entire chip and the core set was statistically compared (Fig. 3). Nine categories were differentially represented at significant levels during berry development and five of these were overrepresented in the regulated gene core set. These were "carbohydrate metabolism", "cell wall organization and biogenesis," "response to stimulus," "secondary metabolism" and "photosynthesis," all processes that ought to be specifically regulated during grape berry development. The four underrepresented categories were "nucleobase, nucleoside, nucleotide and nucleic acid metabolism," "protein biosynthesis," "protein metabolism" and "transport." These classes may either play a minor role in the developmental program or their involvement could be restricted to a small subset of genes.
The functional class distribution of the Vitis GeneChip® was also compared to that of the Arabidopsis genome available at TAIR (, [Additional file 5]). The overall agreement in the class distributions confirms the Vitis GeneChip® as a reliable toolkit for the investigation of the grape expressed genome.
Cluster analysis and functional categories modulation during berry ripening
Cluster analysis of the gene core set was based on the k-means method using Pearson's correlation distance calculated on the gene expression profiles. Transcripts were divided into eight groups representing the minimum number of profiles that could be obtained with three time-points. We observed good agreement between clustering in the three gene sets (80% of the transcripts fell in the same cluster in all seasons). We decided to use the 2005 expression data for the cluster representation (Fig. 4), due to the smaller variance among replicates [see Additional file 2]. As expected from PCA, the two most populated clusters were number 4 and 8, composed by genes modulated positively and negatively from TP B to TP C, respectively. Conversely, clusters 2 and 6, consisting of genes modulated only during the pre-véraison interval, were the less populated ones. Approximately the same number of genes was positively or negatively modulated along the whole study interval (cluster 3 and 7). Genes specifically induced around véraison should fall in cluster1 while genes with roles in early and late stages of berry development should fall in cluster 5.
Functional class distribution frequency was then calculated for each cluster and represented as histogram [Additional file 6]. To better understand the modulation of the functional classes during berry ripening, their frequencies were determined after splitting the experimental interval into two phases: one from TP A to TP B and one from TP B to TP C. The genes modulated during the first phase derived from clusters 1, 2, 3 (induced) and 5, 6, 7 (repressed), while the genes modulated during the second phase derived from clusters 3, 4, 5 (induced) and 1, 7, 8 (repressed). During the pre-véraison phase, the number of the induced genes (408) exceeded that of the repressed ones (270), while the opposite occurred during the following phase (651 induced vs. 1118 repressed). As depicted in Figure 5, all but two categories showed a general tendency towards an opposite modulation between the first and the second phase. Furthermore in the first phase a general bias towards induction for the GO categories involved in regulatory mechanisms, namely "nucleobase, nucleoside, nucleotide and nucleic acid metabolism" (mainly transcription factors), "response to endogenous stimulus" (mainly hormone metabolism), "signal transduction," and "protein metabolism" was observed. This suggests a strong cell re-programming taking place in berry cells up to véraison. During the second phase, a marked negative regulation is evident for categories involved in cell division, such as "organelle organization and biogenesis" and "protein biosynthesis," and "photosynthesis." This behaviour is coherent with the slowing down of cell replication and the loss of photosynthetic capacity. On the other hand, the classes "amino acid and derivative metabolism," "carbohydrate metabolism," "cell wall organization and biogenesis," "development," "lipid metabolism," "primary" and "secondary metabolism" display a positive regulation. This indicates the prevalence after véraison of metabolic processes involved in cell wall loosening and synthesis, sugar accumulation and synthesis and transport of metabolites responsible of grape colour and flavour. Finally, the class "response to stimulus" appeared induced in both phases, suggesting that throughout the entire berry development process the cell devotes considerable effort to face different kinds of stimuli, presumably including osmotic, oxidative and biotic stress.
Gene regulation of grape berry ripening
Our study of the berry ripening process identified 125 transcription factors, 65 genes involved in hormone metabolism and response and 92 involved in signal transduction, corresponding to 8.5%, 4.4% and 6.2% of the regulated core set, respectively. These figures and similar data recently obtained for tomato and Citrus fruit [16, 17], suggest that a consistent fraction of the modulated fruit transcriptome is devoted to the control of the developmental program. In this section we present our main contribution to the knowledge of gene regulation during berry ripening in grapevine. A restricted list of selected genes involved in the process and discussed in the text is reported in Table 1 with the specific AffyID, whereas the complete list is provided in part A of Additional file 7.
Hormones metabolism and signalling
Amongst the genes related to hormone metabolism in the core set, those related to auxin and ethylene were the most represented, followed by those related to abscisic acid and brassinosteroids. Few genes related to cytokinines, gibberellins or jasmonic acid were found. These observations are consistent with our knowledge of the relative importance of these signalling pathways during grape berry development [18, 19].
In fruits from both climacteric and non-climacteric species, auxins are known to be more abundant immediately after ovary fertilization and to mediate important signals for the onset of fruit development [20, 21]. In grape, it has been generally accepted that indole-3-acetic acid (IAA) content reaches its maximal level just after anthesis and then declines to very low levels in the ripe fruit , even if a recent study on Vitis vinifera cv Cabernet Sauvignon did not confirm it . The expression of genes involved in cell wall metabolism (GRIP4), anthocyanin synthesis (CHS and UFGT) and sugar storage (GIN1) and the accumulation of the correlated metabolites were all shown to be retarded by auxin treatments at véraison .
In the present study, we have identified a gene homologous to the Arabidopsis amidase AtAMI1, which in vitro synthesizes IAA from indole-3-acetamide , and its descending profile during ripening agrees well with the decline in IAA levels observed in Vitis labrusca berries . We have also found a gene homologous to an Arabidopsis IAA-amino acid synthetase which may contribute to IAA intracellular homeostasis via amino acid conjugation of excess IAA . Two auxin carriers (an AUX1-like and a PIN1-like), putatively mediating auxin efflux, were also expressed before véraison.
We detected numerous auxin-regulated genes belonging to the Aux/IAA, Small Auxin-Up RNA and GH3 gene families. These families are characterized by the presence of an Auxin Responsive Element in the promoter region of their members, recognized by specific Auxin Response Factors (ARFs). Two ARFs, ARF5 and ARF18 and an auxin receptor of the ABP family were found expressed during the pre-véraison stage and then steeply repressed during ripening. With the exception of two Aux/IAA proteins, which were induced during ripening, all the other transcripts presented profiles with either a peak of expression around véraison or a decreasing profile after it. The auxin responsive genes seem therefore timely co-ordinated with the cellular auxin level, reflecting its prevalent role in the first stages of berry development. We have also detected a gene (1612001_s_at) homologous to a Capsicum chinense GH3 gene, which has been proposed as point of convergence between auxin and ethylene signals in non-climacteric fruit development . While in pepper this GH3 transcript is up-regulated during ripening due to its ethylene-inducibility, our analyses showed a steep repression of the grape homologue after TP B.
The role of ethylene during development and ripening of some non-climacteric fruits was suggested since the 1970s and has recently been well documented [27–29]. In grape, a small and transient increase of endogenous ethylene production occurs just before véraison, in concomitance with an increase in 1-aminocyclopropane-1-carboxylic acid (ACC) oxidase activity . Our results well suit these reports, as we observed a peak of expression of ACC oxidase around véraison, while ACC synthase and a putative ethylene-forming-enzyme dioxygenase, detected with similar profile also by others , were expressed until véraison and then repressed.
We were not able to identify any receptor or kinase involved in ethylene signal transduction in grape, within the modulated gene set. Instead, we detected four transcription factors of the Ethylene Responsive Factors (ERFs) and AP2/EREBP families. The transcript homologous to tomato LeERF1 and a putative AP2/EREBP transcription factor were induced during berry ripening, while the homologous of ERF3b and another putative AP2/EREBP transcription factor were repressed. Considering that in tomato, tobacco and Arabidopsis, ERF1 genes are induced by ethylene and wounding and positively regulate downstream genes involved in biotic and abiotic stresses , the induced transcription factors that we have identified could be involved in ethylene mediated responses to plant stresses in grape as well [31, 32]. The Vitis homologues to osmotin, PR-1 and PR-4, downstream genes known to be regulated by ERFs in the species above mentioned, were indeed found to be induced in our study, suggesting that they could be under the same regulatory mechanism. The homologous of tobacco nitrilase NtNIT4B , another gene related to the transcriptional control of ethylene responsive defence genes, was found to be induced during ripening in our analysis as well. We could not thus exclude that an additional role of ethylene during grape berry development is related to biotic stress response.
Sequence annotation revealed the modulation of some transcripts homologous to the tomato genes ER6 and ER33, whose functions are still unknown but are reported to be induced by ethylene during tomato ripening . In our experiment, however, these genes did not show strong modulation around véraison and they were all repressed during ripening. On the other hand, a transcript homologous to tomato E8, was found induced in grape, in agreement with tomato and strawberry gene expression data [35, 36]. Although E8 resembles ACC synthase structure, it apparently exerts a negative control on ethylene biosynthesis  and its expression seems affected by the homeotic transcription factor LeMADS-RIN . E8 similar modulation in both climacteric and non-climacteric fruits raises the possibility of a common regulatory mechanism for ethylene homeostasis control in fleshy fruit ripening.
The phytohormone abscisic acid (ABA) regulates processes such as floral transition, embryo maturation, seed development and tolerance to abiotic and biotic stress . In grape berries, ABA begins to accumulate at véraison and is thus considered a good candidate for triggering berry ripening . The present analysis found three ABA-related transcripts induced after véraison: one is homologous to the pea nine-cis-epoxycarotenoid dioxygenase 4 and thus putatively involved in ABA synthesis; the other two are homologous to two ABA responsive proteins. Recently, a Vitis dehydrin, DHN1 , and some Arabidopsis plasma membrane aquaporins (PIPs) , have been reported to be regulated by ABA and involved in stress and senescence, respectively. Interestingly, in our analysis we found VvDHN1 and numerous VvPIP isoforms modulated in a ripening specific way (clusters 4 and 8).
Brassinosteroids accumulate during fruit development and seem to play a key role in determining the onset of ripening in fleshy fruits [23, 41]. Three Vitis genes involved in brassinosteroid biosynthesis and sensing have been recently cloned and transcriptionally characterized during berry development . Our analysis identified VvBR6OX1, which converts 6-deoxocastasterone to castasterone, the only bioactive brassinosteroid detected in grape , and VvBRI1, a brassinosteroid receptor. Their expression profiles were characterized by a peak of induction at TP B, in agreement with previous data .
We identified numerous transcripts putatively involved in the circadian rhythm oscillatory system. In particular, one member of the His-to-Asp two-component signal transduction family, homologous to Arabidopsis APRR2, and one transcription factor of the Constans-like family, showed a decreasing profile after véraison, while a kinase homologous to WNK1 , and two transcription factors, belonging to the Constans-like and Early flowering family respectively , displayed a positive regulation, more marked after véraison. A transcript homologous to Flowering locus T, putatively coordinating the circadian rhythm with the flowering switch , appeared expressed until TP B, thus suggesting its involvement in signalling circuits during early berry development.
Some transcripts putatively involved in light and hormonal signalling cross-talk were also identified. IAA26/phytochrome associated protein 1 and IAA27/phytochrome associated protein 2, involved in light and auxin signal transduction, were down regulated after véraison. Members of the His-to-Asp family, such as Type A response regulators, pseudo response regulators and histidine-containing phosphotransfer proteins, are known to mediate the responses to light, osmotic, cytokinin and ethylene stimuli (reviewed in ) and were found modulated during the process of berry development.
In our study, we identified 125 transcription factors, among which members of the myb, MADS-box, NAC and WRKY families and some homeotic and development-specific genes.
Within the myb family, we found two previously isolated Vitis genes: VvmybA1, which showed a ripening-specific profile highly correlated to those of UDP-glucose:flavonoid 3-O-glucosyltransferase, glutathione-S-transferase 4 and the caffeoyl-CoA O-methyltransferase as described in , and a VlmybB-like gene, which was induced during berry development as reported for Vitis labruscana . The recently characterized grape Vvmyb5A, proposed to be involved in the regulation of the phenylpropanoid pathway , was found to be expressed invariantly at low levels in all years during this study interval (1621471_s_at), thus not supporting a ripening-specific function for this gene.
Concerning the MADS-box family of transcription factors, only one previously characterized grape gene was present in the regulated set, namely VvMADS5, which was highly induced during early berry development, in agreement with . Both VvMADS1 (1619742_at) and VvMADS4 (1614965_at), which were previously reported as abundant transcription factors involved in berry development [49, 50], were highly expressed at all time-points in our experiment and thus not included in the regulated set.
1613748_at could represent the homologue of tomato LeMADS-rin (E-value = 8.2e-50, 68% identity), which is part of the developmental signalling system that initiates ripening in this fruit . However, while LeMADS-rin is induced from breaker stage to ripeness, its putative grape homologue was slightly induced until véraison and then showed a decreasing profile. This behaviour could either suggest that the grape transcript is not the homologue of LeMADS-rin or that its role is different in non-climacteric fruits.
Five transcripts with homology to NAC transcription factors appeared modulated, all in a positive way in the study interval. This is of interest as recent transcriptome studies report the regulation of NAC transcription factors family not only during biotic and abiotic stress responses, but also in many other processes, such as fruit development and ABA signalling (reviewed in ).
Three transcripts homologous to the Arabidopsis WRKY family were also modulated. In plants, the WRKY family is quite large (74 members in Arabidopsis and 90 in rice) and its members participate in numerous cellular processes, such as defence and hormonal signalling, developmental programs and fruit maturation (reviewed in ).
Nine putative homeotic genes previously identified in floral development studies and 11 transcripts involved in other morphological and/or developmental processes were modulated in our study although their connection to fruit development has not yet been reported and will need further investigation. 1606591_at and 1619334_at correspond to two transcriptional regulators, SHORT ROOT and SCARECROW, which in Arabidopsis are both implicated in developmental root layers specification . As they were co-expressed in our analysis, they might interact in berry as well.
Though not present in the regulated gene set, 1616455_s_at, corresponding to the grape VvMSA transcription factor ought to be mentioned. VvMSA is induced by sugar and ABA and positively regulates VvHT1 expression [55, 56]. In our study, 1616455_s_at was very highly expressed at all three time-points, in agreement with previous proteomic and gene expression studies [7, 55, 57].
Genes responsible for berry phenotypic traits and stress response
Berry development involves major transcriptional changes of genes directly responsible for the biochemical and phenotypical characteristics of each growth stage. Several protein families linked to these phenotypic traits have been expansively described in past studies focussed on single or small gene sets  and recently in a genome-wide study . The present work provides an extensive description of gene expression during ripening, reporting for the first time the occurrence of an oxidative burst immediately after véraison. The complete list of the genes discussed in the text is in part B of Additional file 7, ordered by GO functional category; for clarity, they are further divided according to more specific GO metabolisms, when possible.
Cell wall metabolism
The previously described involvement of specific gene families in berry growth and softening  was confirmed in the present study. Among these, there were pectin modifying enzymes (methylesterase, invertase/pectin methylesterase inhibitor, pectate lyase, pectinesterase), extensins and expansins, grape ripening-induced proteins (GRIP 3, 13, 15, 22 and 28) and xyloglucan endotransglycosylases (XETs) [8, 59]. In particular, 10 members of the XET family were modulated and four of these very strongly during ripening (1613415_at, 1615809_at, 1619355_at, 1621251_s_at, cluster 4). These induced isoforms are homologous to the tomato SIXTH5, 6 and 8 XET isoforms, recently defined as Group 3 XETs . Although from the analysis of the tomato fruit cDNA libraries only SIXTH5 was found expressed , our array-based gene expression analysis suggests that all the three isoforms of group 3 could be involved in fruit development sharing a common function. Vitis XET isoforms have been proposed to be implied in brassinosteroids mediated berry cell wall softening .
The present analysis identified 270 transcripts related to plant stress perception and response which fall into the "response to stimulus" GO category.
Sixty-nine transcripts were involved in plant biotic stress response, and they were mostly accumulating during ripening. The induction of defence genes including GRIPs (GRIP 22, 28, 31, 32, 51, 61), thaumatins, lipid transfer proteins, chitinases, γ-thionins and pathogenesis related proteins during berry ripening, confirmed previous results .
A particular kind of stress response which has been reported during fruit development in several species such as tomato , strawberry , avocado , pineapple  and pear , is the response to oxidative stress. In recent years a dual role for reactive oxygen species (ROS) in plants has been identified: toxic by-products of the aerobic metabolism and important regulators of growth, development and defence reactions. A complex network of genes finely tuning ROS concentration as the result of their production and scavenging has been documented .
The occurrence of oxidative stress during grape berry development has been rather controversial: catalase was identified among the abundant proteins in the ripe berry , but at the transcriptional level the typical oxidative stress markers seemed absent or negatively regulated . The present study clearly indicated that oxidative processes take place also during berry ripening and that an enzyme-mediated scavenging system is activated. Berry H2O2 content increased significantly in correspondence of the véraison stage, reaching its maximum one-two weeks after, and then decreasing at a slower pace towards harvest time (Fig. 6). At least 32 transcripts of the modulated core set appeared involved in the enzymatic detoxification from H2O2 and the other ROS species which accumulate in the cell starting from véraison (Fig. 7). Some of them, such as ascorbate peroxidase, glutathione peroxidase and peroxiredoxins, catalyze the direct reduction of hydrogen peroxide to water; the others (thioredoxins, glutaredoxins, glutathione-S-transferases and metallothioneins) regulate the balance of the oxidized and reduced forms of the antioxidants ascorbic acid and glutathione. Although the regulation of the detoxifying enzymes must be very complex, isozyme specific, and occurring at different levels (transcriptional- and post-transcriptional regulation, subcellular compartimentalization, etc) as previously shown for some of them , most of the enzymes shown in Figure 7 appeared to be unchanged from TP A to TP B and be strongly modulated from TP B to TP C. Further studies are needed to address how ROS molecules are taking part to these regulatory circuits and what is the main reason of ROS accumulation at véraison. The superoxide dismutase (SOD) and catalase (CAT) transcripts which also codify for two important detoxifying enzymes were not significantly modulated in our study. The CAT gene (1618705_s_at) was highly expressed in all the three TPs, the SOD encoding transcripts (5 probes) showed very different expression levels but no one fell into the modulated core set.
The Nudix gene family may also be involved in the oxidative stress response, as three transcripts were found modulated (one in cluster 3 and two in cluster 8). In Arabidopsis, Nudix hydrolases appear to be involved in the sanitization of the oxidized nucleotide pool accumulating during oxidative stress .
Ten genes involved in programmed cell death (PCD) appeared to be modulated. In particular, three transcripts with putative function of PCD suppression are grouped in cluster 1, while a putative death associated protein is induced after véraison (cluster 4). In addition, listed in the protein metabolism category there are two cysteine proteases (cluster 4) and a cystatin-like gene (cluster 5). Our data are insufficient to draw any conclusion, however they may suggest the occurrence of a PCD process in grape berry during ripening, as highlighted for strawberry .
Sugar/acid ratio and water balance
Our results support the model of sugar phloem unloading via an apoplastic pathway by means of plasma membrane sucrose and hexoses transporters as proposed by  and recently confirmed by [68, 69]. Four sugar transporters were found modulated: the hexose transporter VvHT6 (1615257_at and 1619691_at) was positively modulated during ripening (cluster 3), as previously described ; interestingly, a transcript homologous to a sugar transporter isolated in Citrus fruit was also induced while the forth transcript, which is homologous to a sugar transporter from tomato root is repressed. The sucrose transporter VvSUC27 appeared dramatically down-regulated after véraison, as already reported in the cv. Syrah , while, VvSUC11 and VvSUC12, transporters with a ripening specific expression profile , were highly expressed at all time points in our experiment.
There is evidence that during berry ripening, sucrose is hydrolysed in the cell wall and then re-synthesised in the cytoplasm following hexose uptake . Our data corroborate this mechanism as three sucrose-phosphate synthases (cluster 3 and 4), a sucrose synthase (cluster 4), a sucrose-6-phosphate phosphatase (cluster 4), a UDP-sugar pyrophosphorylase (cluster 2) and a beta-phosphoglucomutase (cluster 4) were all up-regulated during ripening. Vacuolar invertases, VvGIN1 and VvGIN2, both resulted highly expressed until véraison (cluster 8), as previously described [57, 71]. The existence of 'futile cycles' of sucrose synthesis and breakdown has been demonstrated in other fruits, such as tomato, where the relative rate of synthesis and breakdown cycles controls the rate of sugar and starch accumulation .
The early expression of a phosphoenolpyruvate carboxylase isoforms (cluster 8) well correlated with the rate of malate accumulation during the first phase of berry development. Furthermore, a vacuolar pyrophosphatase (1614834_at) putatively involved in malate vacuolar uptake was expressed early and then induced during ripening (cluster 5), as previously observed . At véraison, when malate decompartimentalization and breakdown take place , induction of malate degrading enzymes, such as malic enzyme and phosphoenolpyruvate carboxykinase (cluster 1) and mitochondrial malate dehydrogenase (cluster 8) was observed. The profiles observed for grape alcohol dehydrogenase 2 and a putative short chain alcohol dehydrogenase (cluster 4 and 5, respectively) and three isoforms of aldehyde dehydrogenase (one in cluster 4 and two in cluster 8) may be indicative of a shift to an aerobic fermentative metabolism during ripening . The alpha subunit of the pyruvate dehydrogenase complex was induced after véraison (cluster 5); however, this enzyme is known to be highly regulated post-translationally, allosterically by NADH and covalently by reversible phosphorylation. Our analysis identified two mitochondrial pyruvate dehydrogenase kinase isoforms, one of which induced in a ripening-specific way. Two enzymes of the glycolysis/gluconeogenesis pathways, such as enolase and glyceraldehyde-3-phosphate dehydrogenase (cluster 5 and 4, respectively) were induced during ripening, suggesting their involvement into the acid/sugar interconversion and/or the fermentative metabolism.
About 90% of ripe berries fresh weight is comprised of water whose transport into the berry vacuole is mediated by aquaporins. In our study we detected 13 transcripts coding for aquaporin isoforms: eight belong to the plasma membrane intrinsic protein family, four to the tonoplast intrinsic protein family and one to the NOD26-like intrinsic protein family. All isoforms were highly expressed in the pre-véraison phase and three of them, coding for the putative aquaporin PIP2-1, were also induced during ripening. A similar modulation has been recently observed in Citrus, although the induced isoforms during fruit development, belong to the TIP family .
The "secondary metabolism" GO functional category included genes involved in the synthesis of aromatic and volatile compounds, such as terpenes and benzenoid compounds; antioxidant compounds, such as polyphenols and vitamin E, and genes of the phenylpropanoid pathway. Due to its importance for berry quality traits, we restricted our analysis to the latter pathway. Most structural genes of this branched pathway, which starts with phenylalanine ammonia lyase, were present among the regulated genes of the current study. During ripening, the genes coding for the enzymes acting in flavonols synthesis (flavonol synthase, cluster 3), stilbenes synthesis (three isoforms of stilbene synthase, cluster 4) and anthocyanins synthesis (UDP-glucose:flavonoid 3-O-glucosyltransferase, cluster 4) were induced, while those involved in tannins synthesis (anthocyanidin reductase and leucoanthocyanidin reductase) were repressed. The profiles of genes coding for phenylalanine ammonia lyase, chalcone isomerase 1 and flavonoid 3'5'-hydroxylase showed a decrease until véraison followed by an induction during ripening, confirming previous data . Two isoforms of 4-coumarate-CoA ligase with opposite profiles were detected (clusters 3 and 8) and their pattern of expression is consistent with an involvement into the anthocyanin and lignin synthesis pathways, respectively.
Eleven members of the ATP-binding cassette (ABC) and multidrug and toxic compound extrusion (MATE) transporters families which are involved in several cellular transport processes (reviewed in ), were modulated. Among them, the three transporters positively modulated during berry development (1610275_at, 1607763_at, 1616929_at) represent possible candidates for secondary metabolites transport.
Gene expression profiles highly influenced by seasonal variation
The influence on gene expression by environmental conditions is well reported in literature  and was evident in our PCA analysis (Fig. 2). The comparison of the three sampling years highlighted a large number of genes modulated only during one or two seasons, and thus excluded by the core set of ripening-specific genes: 938 genes in 2003, 2530 in 2005 and 2143 in 2006. Two examples are the isoform 2 of the phenylpropanoid biosynthetic genes chalcone synthase and chalcone isomerase whose expression profiles are depicted in Fig. 8E and 8F; such discrepancy in their modulation is indeed expected for non-ripening specific isoforms. Although these genes are probably highly affected by changes in the environment their GO functional categories distribution was not significantly different to that of the core set genes (data not shown). This result suggests that, on a long time frame, the plant reacts to the seasonal variations by adjusting the whole metabolism, and not just a part of it, to maintain homeostasis.
The group of genes modulated in all three seasons with different transcriptional profiles, deserves in our opinion a special discussion as they are likely involved into the ripening process but in the meantime highly influenced by the climatic conditions. These 223 genes represent about 13% of the modulated gene set (1700 transcripts). Among them, at least four transcripts are putatively involved in light perception and signalling: cryptochrome 1 (1615199_at) and pseudo-response regulator 9 (APRR9, 1616872_at) showed opposite behaviours in 2003 vs. 2005 and 2006 (Fig. 8A and 8B); pseudo-response-regulator 7 (APRR7, 1608006_at) displayed a "V" profile in 2003 and 2005 and the opposite pattern during 2006; finally LONG HYPOCOTYL 5 (HY5, 1609930_at) showed a different profile from TP A to TP B in 2003 and 2005 vs. 2006 while it was repressed in all three seasons from TP B to TP C. Interestingly, APRR7 and APRR9 seem involved in the temperature responsiveness of the Arabidopsis circadian clock, synchronizing photo- and thermo-cycles  and HY5 is a positive regulator of light signalling which acts downstream of several photoreceptors and mediates the response to different types of light (red/far red, blue and UV-B [79, 80]).
Out of the three previously characterized grape hexose transporters present on the chip, two (Vvht1 and Vvht2) seemed influenced by the season. Vvht1 (1616083_at) showed a "V" profile in all the three years, but with different slopes, suggesting a strong environmental conditioning along with the reported cultivar specificity [10, 81, 82]. Unlike Vvht1, Vvht2 is considered a ripening specific isoform [10, 82, 83]; however, the present work highlighted a certain degree of variability in its induction time in the three years (1614764_at, Fig. 8C). A similar behaviour was observed for another ripening specific gene, namely the isoform 3 of the chalcone synthase family [75, 84] (Fig. 8D).
To ensure that the observed differences among seasons and among gene isoforms were not technical artefacts, the expression profile of six genes was validated by real time RT-PCR (Fig. 8, insets). In all but one case (Vvht2 year 2003) there was complete agreement with the array data.
The time-course gene expression analysis of Pinot Noir berry development, based on a systematic and genome-wide approach, has identified the occurrence of two well distinct phases along the process. The pre-véraison phase represents a reprogramming stage of the cellular metabolism, characterized by the expression of numerous genes involved in hormonal signalling and transcriptional regulation. The post-véraison phase is characterized by the onset of a ripening-specialized metabolism responsible for the phenotypic traits of the ripe berry. A similar switch in fruit metabolism has been recently reported also in Citrus . Five functional categories are overrepresented in this ripening-specialized metabolism: cell wall, carbohydrate and secondary metabolisms, stress response and photosynthesis; the latter being obviously switched off. An oxidative burst, previously not detected in grapevine and characterized by a rapid accumulation of H2O2 starting from véraison, was demonstrated. A concurrent modulation of the enzymatic ROS scavenging network was also highlighted, suggesting a possible role of ROS species in the regulation of berry ripening. We have found a very large number of transcription factors and transcripts related to hormonal metabolism and signal transduction that appear to be regulated during berry development. These represent 19% of the modulated gene set, corroborating the widely-accepted hypothesis that fruit development is under tight transcriptional control. These data provide the basis for subsequent correlation analysis aimed at highlighting gene co-regulation and metabolic networks . Our strategy to perform the experiment during three seasons allowed us to restrict the dataset of modulated genes to those more likely to be ripening-specific. Nonetheless, we were able to separate ripening specific isoforms within gene families and to identify ripening related genes which appeared strongly regulated also by the seasonal weather conditions.
In perspective, the wealth of information we present here will provide a new extended platform to study the complex process of grapevine berry development and more in general of non-climacteric fruits.
Plant Material and Experimental Design
Ten grape clusters of Vitis vinifera cv. Pinot Noir were collected weekly between 8 a.m and 9 a.m, from flowering to over-ripening during 2003, 2005 and 2006 seasons at the IASMA study site 'ai Molini' (San Michele all'Adige-TN Italy). To be more representative of the canopy five clusters were collected from the north side (shady) and the other five from the south side (sunny). Berries were rapidly detached with their petioles and divided randomly into groups of 50 that were either immediately frozen for subsequent transcriptional analysis or subjected to biochemical analyses, including must acid and sugar content, average berry weight determination and anthocyanin concentration . Based on these parameters, the peak of acidity just before véraison was estimated and two additional time-points (two weeks before and three weeks after véraison) were chosen for gene expression analysis. The three TPs correspond to the 33, 34 and 36 of the modified E-L system . For each time-point three biological replicates were obtained by performing an independent RNA extraction from three sub-pools of 8 berries derived from the group of 50 berries. Consequently nine arrays were utilized for each season.
RNA Preparation, Array Hybridization, Data Analysis
Total RNA was extracted from frozen deseeded berries according to  and further purified using RNeasy spin columns (Qiagen). RNA quality and quantity were assessed by gel-electrophoresis and by absorbance measurements. Ten μg of each sample were sent to the IFOM (Milan-I) Affymetrix platform facility for probe synthesis, GeneChip® Vitis arrays hybridization, washing, staining and scanning with GeneChip Scanner 3000, according to the 'Affymetrix GeneChip® Expression Analysis Technical Manual'. Data referred to 2003, 2005 and 2006 were analyzed independently in order to obtain one set of differentially expressed genes for each season. Bioconductor packages  in R v.2.1.0  were used for the statistical analyses. Data were pre-processed using the GCRMA method , which performs background correction, quantile normalization and summarization, and returns log2 converted expression values. Biological replicates quality was assessed by means of R-squared coefficient, which represents the fraction of variance explained by a linear model. For variation assessment of the gene expression values between all 27 data sets, data from the three seasons were scaled and Principal Components Analysis (PCA) was performed. In order to identify genes modulated during berry development, a filtering step based on the inter-quantile range method (IQR = 0.25) was used. This restricted dataset, enriched in differentially expressed genes, was suited to run a multiclass comparison method of Significance Analysis of Microarrays (SAM), with false discovery rate =0.15% . SAM output was further restricted to genes with fold-change greater than 2 in at least one of the two comparisons (TP A vs. TP B and TP B vs. TP C). The final datasets relative to the three seasons were intersected obtaining a dataset of 1700 common modulated genes. Pearson's correlation coefficients between genes profiles relative each pair of seasons were then calculated and a threshold of 0.5 was fixed to obtain a berry ripening-modulated core dataset of 1477 genes. Expression profiles clustering was carried out by the k-means method (k = 12) with Pearson's correlation distance using the T-MeV software . The number of clusters was manually reduced to eight, which represents the minimum number of expression profiles considering three time-points.
Affymetrix Vitis Array
The GeneChip® Vitis vinifera genome array (Affymetrix) consists of 16,436 probesets: 14,496 derived from V. vinifera transcripts and 1,940 derived from other Vitis species or hybrids transcripts. It can interrogate 12,908 GenBank accessions of V. vinifera and 1,547 of other Vitis species or hybrids. Sequences used in the design of the Vitis GeneChip® were selected from GenBank, dbEST, and RefSeq. The sequence clusters were created from the UniGene database (Build 7, Oct 2003). V. vinifera sequences represented on the chip correspond to 10,042 TIGR Tentative Consensus and 1,940 Singletons (Release 4, Sept 2004), while 102 GenBank accessions are not present in the TIGR database. Overall chip redundancy is estimated to be 16.6%.
The consensus sequence of each Vitis GeneChip® probeset (available at  was analyzed to identify the coding sequence (cds) by means of an in-house developed method, which combines the results of two cds predictors: Estscan  and FrameFinder (at ESTate, ). 13,923 cds were obtainded and used to interrogate the Uniprot database (Sept. '05) by means of 'blastp', while the remaining 2,513 sequences were used to interrogate the same database by means of 'blastx' . Blast results (E-value <e-10) with GO associated terms () were analyzed by the in-house developed program 'GORetriever' for annotation and 12,242 sequences were annotated with high confidence (> 80%). The automatic annotation results for the 1700 modulated genes were manually inspected and integrated with GO 'biological process' terms supported by literature evidences.
Genes were grouped into 17 functional categories based on GO 'biological process' terms by means of the in-house developed program 'GOSlimer' [see Additional file 4 for the list of the GO categories and Additional file 3 for genes annotation]. Functional categories distribution in the modulated and chip sequences were compared by means of Chi square and Fisher statistical tests (p-value < 0.001).
Determination of the H2O2 content
H2O2 content was measured by a fluorimetric assay based on the substrate 10-acetyl-3,7-dihydroxyphenoxazine which is oxidized in the presence of H2O2 and peroxidase, using the Amplex Red Hydrogen Peroxide/Peroxidase Assay (Molecular Probes, Eugene, USA) following the manufacturer's instructions. Ten frozen berries for each time point considered were ground in liquid nitrogen and 0.2 g of fine powder were then solubilized, for each replicate, with 0.5 ml of 50 mM phosphate buffer (pH 7.4) and kept 5 min on ice. After centrifugation at 20,000 g for 15 min, the cleared supernatant was extracted with an equal volume of 2:1 (v/v) chloroform:methanol mixture and centrifuged at 12,000 g for 5 min. To perform the H2O2 assay 50 μl of the aqueous phase were added to 50 μl of the reagent working solution. After 30 min of incubation at 25°C, relative fluorescence (excitation at 544 nm and emission at 590 nm) was measured in a black 96 wells plate. Absolute quantification was determined by the use of a H2O2 standard curve.
Real time RT-PCR
First strand cDNA synthesis was performed on one mRNA sample of the three biological replicates used for the microarray experiments using the SuperScript™ III Reverse Transcriptase kit (Invitrogen) according to the manufacturer instructions. Primers [see Additional file 8 for sequences] and cDNA were mixed with the Platinum® SYBR® Green qPCR SuperMix-UDG (Invitrogen) and the reaction was carried out on an ABI PRISM 7000 Sequence Detection System (Applied Biosystems). Cycling conditions were: 50°C for 2 minutes, 95°C for 2 minutes, then 40 cycles of 95°C for 15 seconds and 60°C for 1 minute. Raw data were analyzed with ABI PRISM 7000 SDS software to extract Ct values and with the LinReg software to calculate the reaction efficiency . Relative expression of each gene (target) was then calculated according to the equation by  using actin for normalization (reference) and centered on the mean Ct calculated on the three seasons and on the three time points (control): RelExp = E(target)exp[ΔCt(target)(sample-control)]/E(reference)exp[ΔCt(reference)(sample-control)].
All microarray expression data produced by this work are available at EBI ArrayExpress  under the series entry E-MEXP-1282.
Expressed Sequence Tag
National Center for Biotechnology Information
Significance Analysis of Microarrays
The Institute for Genomic Research
The Arabidopsis Information Resource
Doligez A, Adam-Blondon AF, Cipriani G, Di Gaspero G, Laucou V, Merdinoglu D, Meredith CP, Riaz S, Roux C, This P: An integrated SSR map of grapevine based on five mapping populations. Theor Appl Genet. 2006, 113: 369-382. 10.1007/s00122-006-0295-1.
IASMA Genomics. [http://genomics.research.iasma.it]
NCBI dbest. [http://www.ncbi.nlm.nih.gov/dbEST/]
Coombe BG, Mccarthy MG: Dynamics of Grape Berry Growth and Physiology of Ripening. Aust J Grape Wine Res. 2000, 6: 131-135.
Ablett E, Seaton G, Scott K, Shelton D, Graham MW, Baverstock P, Lee LS, Henry R: Analysis of grape ESTs: global gene expression patterns in leaf and berry. Plant Sci. 2000, 159: 87-95. 10.1016/S0168-9452(00)00335-6.
Moser C, Segala C, Fontana P, Salakhudtinov I, Gatto P, Pindo M, Zyprian E, Toepfer R, Grando MS, Velasco R: Comparative analysis of expressed sequence tags from different organs of Vitis vinifera L. Funct Integr Genomics. 2005, 5: 208-217. 10.1007/s10142-005-0143-4.
da Silva FG, Iandolino A, Al-Kayal F, Bohlmann MC, Cushman MA, Lim H, Ergul A, Figueroa R, Kabuloglu EK, Osborne C, Rowe J, Tattersall E, Leslie A, Xu J, Baek J, Cramer GR, Cushman JC, Cook DR: Characterizing the grape transcriptome. Analysis of expressed sequence tags from multiple Vitis species and development of a compendium of gene expression during berry development. Plant Physiol. 2005, 139: 574-597. 10.1104/pp.105.065748.
Davies C, Robinson SP: Differential screening indicates a dramatic change in mRNA profiles during grape berry ripening. Cloning and characterization of cDNAs encoding putative cell wall and stress response proteins. Plant Physiol. 2000, 122: 803-812. 10.1104/pp.122.3.803.
Waters DL, Holton TA, Ablett EM, Lee LS, Henry RJ: cDNA microarray analysis of developing grape (Vitis vinifera cv. Shiraz) berry skin. Funct Integr Genomics. 2005, 5: 40-58. 10.1007/s10142-004-0124-z.
Terrier N, Glissant D, Grimplet J, Barrieu F, Abbal P, Couture C, Ageorges A, Atanassova R, Leon C, Renaudin JP, Dedaldechamp F, Romieu C, Delrot S, Hamdi S: Isogene specific oligo arrays reveal multifaceted changes in gene expression during grape berry (Vitis vinifera L.) development. Planta. 2005, 222: 832-847. 10.1007/s00425-005-0017-y.
Grimplet J, Deluc LG, Tillett RL, Wheatley MD, Schlauch KA, Cramer GR, Cushman JC: Tissue-specific mRNA expression profiling in grape berry tissues. BMC Genomics. 2007, 8: 187-10.1186/1471-2164-8-187.
Coombe BG: Adoption of a system for identifying grapevine growth stages. Aust J Grape Wine Res. 1995, 1: 104-110.
Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA. 2001, 98: 5116-5121. 10.1073/pnas.091062498.
Gene Ontology home. [http://www.geneontology.org/]
TAIR:The Arabidopsis Information Resource. [http://www.arabidopsis.org]
Cercos M, Soler G, Iglesias DJ, Gadea J, Forment J, Talon M: Global analysis of gene expression during development and ripening of citrus fruit flesh. A proposed mechanism for citric acid utilization. Plant Mol Biol. 2006, 62: 513-527. 10.1007/s11103-006-9037-7.
Alba R, Payton P, Fei Z, McQuinn R, Debbie P, Martin GB, Tanksley SD, Giovannoni JJ: Transcriptome and selected metabolite analyses reveal multiple points of ethylene control during tomato fruit development. Plant Cell. 2005, 17: 2954-2965. 10.1105/tpc.105.036053.
Coombe BG, Hale CR: The hormone content of ripening grape berries and the effects of growth substance treatments. Plant Physiol. 1973, 51: 629-634.
Davies C, Boss PK, Robinson SP: Treatment of grape berries, a nonclimacteric fruit with a synthetic auxin, retards ripening and alters the expression of developmentally regulated genes. Plant Physiol. 1997, 115: 1155-1161.
Aharoni A, Keizer LC, Van Den Broeck HC, Blanco-Portales R, Munoz-Blanco J, Bois G, Smit P, De Vos RC, O'Connell AP: Novel insight into vascular, stress, and auxin-dependent and -independent gene expression programs in strawberry, a non-climacteric fruit. Plant Physiol. 2002, 129: 1019-1031. 10.1104/pp.003558.
Gillaspy G, Ben-David H, Gruissem W: Fruits: a developmental perspective. Plant Cell. 1993, 5: 1439-1451. 10.1105/tpc.5.10.1439.
Cawthon DL, Morris JR: Relationship of seed number and maturity to berry development, fruit maturation, hormonal changes, and uneven berry ripening in 'Concord' (Vitis labrusca) grapes. J Am Soc Hortic Sci. 1982, 1097-1104.
Symons GM, Davies C, Shavrukov Y, Dry IB, Reid JB, Thomas MR: Grapes on steroids. Brassinosteroids are involved in grape berry ripening. Plant Physiol. 2006, 140: 150-158. 10.1104/pp.105.070706.
Pollmann S, Neu D, Weiler EW: Molecular cloning and characterization of an amidase from Arabidopsis thaliana capable of converting indole-3-acetamide into the plant growth hormone, indole-3-acetic acid. Phytochemistry. 2003, 62: 293-300. 10.1016/S0031-9422(02)00563-0.
Staswick PE, Serban B, Rowe M, Tiryaki I, Maldonado MT, Maldonado MC, Suza W: Characterization of an Arabidopsis enzyme family that conjugates amino acids to indole-3-acetic acid. Plant Cell. 2005, 17: 616-627. 10.1105/tpc.104.026690.
Liu K, Kang BC, Jiang H, Moore SL, Li H, Watkins CB, Setter TL, Jahn MM: A GH3-like gene, CcGH3, isolated from Capsicum chinense L. fruit is regulated by auxin and ethylene. Plant Mol Biol. 2005, 58: 447-464. 10.1007/s11103-005-6505-4.
Hale CR, Coombe BG, Hawker JS: Effects of ethylene and 2-chloroethylphosphonic acid on the ripening of grapes. Plant Physiol. 1970, 45: 620-623.
Trainotti L, Pavanello A, Casadoro G: Different ethylene receptors show an increased expression during the ripening of strawberries: does such an increment imply a role for ethylene in the ripening of these non-climacteric fruits?. J Exp Bot. 2005, 56: 2037-2046. 10.1093/jxb/eri202.
Chervin C, El-Kereamy A, Roustan JP, Latche A, Lamon J, Bouzayen M: Ethylene seems required for the berry development and ripening in grape, a non-climacteric fruit. Plant Sci. 2004, 167: 1301-1305. 10.1016/j.plantsci.2004.06.026.
Tournier B, Sanchez-Ballesta MT, Jones B, Pesquet E, Regad F, Latche A, Pech JC, Bouzayen M: New members of the tomato ERF family show specific expression pattern and diverse DNA-binding capacity to the GCC box element. FEBS Lett. 2003, 550: 149-154. 10.1016/S0014-5793(03)00757-9.
Lorenzo O, Piqueras R, Sanchez-Serrano JJ, Solano R: ETHYLENE RESPONSE FACTOR1 integrates signals from ethylene and jasmonate pathways in plant defense. Plant Cell. 2003, 15: 165-178. 10.1105/tpc.007468.
Ohme-Takagi M, Shinshi H: Ethylene-inducible DNA binding proteins that interact with an ethylene-responsive element. Plant Cell. 1995, 7: 173-182. 10.1105/tpc.7.2.173.
Xu P, Narasimhan ML, Samson T, Coca MA, Huh GH, Zhou J, Martin GB, Hasegawa PM, Bressan RA: A nitrilase-like protein interacts with GCC box DNA-binding proteins involved in ethylene and defense responses. Plant Physiol. 1998, 118: 867-874. 10.1104/pp.118.3.867.
Zegzouti H, Jones B, Frasse P, Marty C, Maitre B, Latch A, Pech JC, Bouzayen M: Ethylene-regulated gene expression in tomato fruit: characterization of novel ethylene-responsive and ripening-related genes isolated by differential display. Plant J. 1999, 18: 589-600. 10.1046/j.1365-313x.1999.00483.x.
Kneissl ML, Deikman J: The tomato E8 gene influences ethylene biosynthesis in fruit but not in flowers. Plant Physiol. 1996, 112: 537-547.
Aharoni A, O'Connell AP: Gene expression analysis of strawberry achene and receptacle maturation using DNA microarrays. J Exp Bot. 2002, 53: 2073-2087. 10.1093/jxb/erf026.
Giovannoni JJ: Genetic regulation of fruit development and ripening. Plant Cell. 2004, 16 (Suppl): S170-S180. 10.1105/tpc.019158.
Zeevaart JAD, Creelman RA: Metabolism and physiology of abscisic-acid. Annu Rev Plant Physiol Plant Mol Biol. 1988, 39: 439-473. 10.1146/annurev.pp.39.060188.002255.
Xiao H, Nassuth A: Stress- and development-induced expression of spliced and unspliced transcripts from two highly similar dehydrin 1 genes in V. riparia and V. vinifera. Plant Cell Rep. 2006, 25: 968-977. 10.1007/s00299-006-0151-4.
Jang JY, Kim DG, Kim YO, Kim JS, Kang H: An expression analysis of a gene family encoding plasma membrane aquaporins in response to abiotic stresses in Arabidopsis thaliana. Plant Mol Biol. 2004, 54: 713-725. 10.1023/B:PLAN.0000040900.61345.a6.
Montoya T, Nomura T, Yokota T, Farrar K, Harrison K, Jones JD, Kaneta T, Kamiya Y, Szekeres M, Bishop GJ: Patterns of Dwarf expression and brassinosteroid accumulation in tomato reveal the importance of brassinosteroid synthesis during fruit development. Plant J. 2005, 42: 262-269. 10.1111/j.1365-313X.2005.02376.x.
Nakamichi N, Murakami-Kojima M, Sato E, Kishi Y, Yamashino T, Mizuno T: Compilation and characterization of a novel WNK family of protein kinases in Arabiodpsis thaliana with reference to circadian rhythms. Biosci Biotechnol Biochem. 2002, 66: 2429-36. 10.1271/bbb.66.2429.
Kikis EA, Khanna R, Quail PH: ELF4 is a phytochrome-regulated component of a negative-feedback loop involving the central oscillator components CCA1 and LHY. Plant J. 2005, 44: 300-313. 10.1111/j.1365-313X.2005.02531.x.
Kardailsky I, Shukla VK, Ahn JH, Dagenais N, Christensen SK, Nguyen JT, Chory J, Harrison MJ, Weigel D: Activation tagging of the floral inducer FT. Science. 1999, 286: 1962-1965. 10.1126/science.286.5446.1962.
Hwang I, Chen HC, Sheen J: Two-component signal transduction pathways in Arabidopsis. Plant Physiol. 2002, 129: 500-515. 10.1104/pp.005504.
Ageorges A, Fernandez L, Vialet S, Merdinoglu D, Terrier N, Romieu C: Four specific isogenes of the anthocyanin metabolic pathway are systematically co-expressed with the red colour of grape berries. Plant Sci. 2006, 170: 372-383. 10.1016/j.plantsci.2005.09.007.
Kobayashi S, Ishimaru M, Hiraoka K, Honda C: Myb-related genes of the Kyoho grape (Vitis labruscana) regulate anthocyanin biosynthesis. Planta. 2002, 215: 924-933. 10.1007/s00425-002-0830-5.
Deluc L, Barrieu F, Marchive C, Lauvergeat V, Decendit A, Richard T, Carde JP, Merillon JM, Hamdi S: Characterization of a grapevine R2R3-MYB transcription factor that regulates the phenylpropanoid pathway. Plant Physiol. 2006, 140: 499-511. 10.1104/pp.105.067231.
Boss PK, Buckeridge EJ, Poole A, Thomas MR: New insights into grapevine flowering. Funct Plant Biol. 2003, 30: 593-606. 10.1071/FP02112.
Boss PK, Vivier M, Matsumoto S, Dry IB, Thomas MR: A cDNA from grapevine (Vitis vinifera L.), which shows homology to AGAMOUS and SHATTERPROOF, is not only expressed in flowers but also throughout berry development. Plant Mol Biol. 2001, 45: 541-553. 10.1023/A:1010634132156.
Vrebalov J, Ruezinsky D, Padmanabhan V, White R, Medrano D, Drake R, Schuch W, Giovannoni J: A MADS-box gene necessary for fruit ripening at the tomato ripening-inhibitor (rin) locus. Science. 2002, 296: 343-346. 10.1126/science.1068181.
Olsen AN, Ernst HA, Leggio LL, Skriver K: NAC transcription factors: structurally distinct, functionally diverse. Trends Plant Sci. 2005, 10: 79-87. 10.1016/j.tplants.2004.12.010.
Ulker B, Somssich IE: WRKY transcription factors: from DNA binding towards biological function. Curr Opin Plant Biol. 2004, 7: 491-498. 10.1016/j.pbi.2004.07.012.
Paquette AJ, Benfey PN: Maturation of the ground tissue of the root is regulated by gibberellin and SCARECROW and requires SHORT-ROOT. Plant Physiol. 2005, 138: 636-640. 10.1104/pp.104.058362.
Cakir B, Agasse A, Gaillard C, Saumonneau A, Delrot S, Atanassova R: A grape ASR protein involved in sugar and abscisic acid signaling. Plant Cell. 2003, 15: 2165-2180. 10.1105/tpc.013854.
Atanassova R, Leterrier M, Gaillard C, Agasse A, Sagot E, Coutos-Thevenot P, Delrot S: Sugar-regulated expression of a putative hexose transport gene in grape. Plant Physiol. 2003, 131: 326-334. 10.1104/pp.009522.
Sarry JE, Sommerer N, Sauvage FX, Bergoin A, Rossignol M, Albagnac G, Romieu C: Grape berry biochemistry revisited upon proteomic analysis of the mesocarp. Proteomics. 2004, 4: 201-215. 10.1002/pmic.200300499.
Robinson SP, Davies C: Molecular biology of grape berry ripening. Aust J Grape Wine Res. 2000, 6: 175-188.
Nunan KJ, Davies C, Robinson SP, Fincher GB: Expression patterns of cell wall-modifying enzymes during grape berry development. Planta. 2001, 214: 257-264.
Saladie M, Rose JK, Cosgrove DJ, Catala C: Characterization of a new xyloglucan endotransglucosylase/hydrolase (XTH) from ripening tomato fruit and implications for the diverse modes of enzymic action. Plant J. 2006, 47: 282-295. 10.1111/j.1365-313X.2006.02784.x.
Jimenez A, Creissen G, Kular B, Firmin J, Robinson S, Verhoeyen M, Mullineaux P: Changes in oxidative processes and components of the antioxidant system during tomato fruit ripening. Planta. 2002, 214: 751-758. 10.1007/s004250100667.
Meir S, Philosoph-Hadas S, Zauberman G, Fuchs Y, Akerman M, Aharoni N: Increased formation of fluorescent lipid-peroxidation products in avocado peels precedes other signs of ripening. J Amer Soc Hort Sci. 1991, 116: 823-826.
Moyle R, Fairbairn DJ, Ripi J, Crowe M, Botella JR: Developing pineapple fruit has a small transcriptome dominated by metallothionein. J Exp Bot. 2005, 56: 101-112.
Brennan T, Frenkel C: Involvement of hydrogen peroxide in the regulation of senescence in pear. Plant Physiol. 1977, 59: 411-416.
Mittler R, Vanderauwera S, Gollery M, Van Breusegem F: Reactive oxygen gene network of plants. Trends Plant Sci. 2004, 9: 490-498. 10.1016/j.tplants.2004.08.009.
Ogawa T, Ueda Y, Yoshimura K, Shigeoka S: Comprehensive analysis of cytosolic Nudix hydrolases in Arabidopsis thaliana. J Biol Chem. 2005, 280: 25277-25283. 10.1074/jbc.M503536200.
Terrier N, Romieu C: Grape berry acidity. 2001, Kluwer Academic Publisher, 35-57.
Cheng GW, Breen PJ: Activity of phenylalanine ammonia-lyase (PAL) and concentrations of anthocyanins and phenolics in developing strawberry fruit. J Am Soc Hortic Sci. 1991, 116: 865-869.
Zhang X-Y, Wang X-L, Wang X-F, Xia G-H, Pan Q-H, Fan R-C, Wu F-Q, Yu X-C, Zhang D-P: A shift of phloem unloading from symplasmic to apoplasmic pathway is involved in developmental onset of ripening in grape berry. Plant Physiol. 2006, 142: 220-232. 10.1104/pp.106.081430.
Davies C, Wolf T, Robinson SP: Three putative sucrose transporters are differentially expressed in grapvine tissues. Plant Sci. 1999, 147: 93-100. 10.1016/S0168-9452(99)00059-X.
Davies C, Robinson SP: Sugar accumulation in grape berries. Cloning of two putative vacuolar invertase cDNAs and their expression in grapevine tissues. Plant Physiol. 1996, 111: 275-283. 10.1104/pp.111.1.275.
Nguyen-Quoc B, Foyer CH: A role for 'futile cycles' involving invertase and sucrose synthase in sucrose metabolism of tomato fruit. J Exp Bot. 2001, 52: 881-889. 10.1093/jexbot/52.358.881.
Terrier N, Sauvage FX, Ageorges A, Romieu C: Changes in acidity and in proton transport at the tonoplast of grape berries during development. Planta. 2001, 213: 20-8. 10.1007/s004250000472.
Mellema S, Eichenberger W, Rawyler A, Suter M, Tadege M, Kuhlemeier C: The ethanolic fermentation pathway supports respiration and lipid biosynthesis in tobacco pollen. Plant J. 2002, 30: 329-336. 10.1046/j.1365-313X.2002.01293.x.
Boss PK, Davies C, Robinson SP: Analysis of the expression of anthocyanin pathway genes in developing Vitis vinifera L. cv Shiraz grape berries and the implications for pathway regulation. Plant Physiol. 1996, 111: 1059-1066.
Martinoia E, Klein M, Geisler M, Bovet L, Forestier C, Kolukisaoglu U, Muller-Rober B, Schulz B: Multifunctionality of plant ABC transporters – more than just detoxifiers. Planta. 2002, 214: 345-355. 10.1007/s004250100661.
Wissel K, Pettersson F, Berglund A, Jansson S: What affects mRNA levels in leaves of field-grown aspen? A study of developmental and environmental influences. Plant Physiol. 2003, 133: 1190-1197. 10.1104/pp.103.028191.
Salome PA, McClung CR: PSEUDO-RESPONSE REGULATOR 7 and 9 are partially redundant genes essential for the temperature responsiveness of the Arabidopsis circadian clock. Plant Cell. 2005, 17: 791-803. 10.1105/tpc.104.029504.
Ulm R, Baumann A, Oravecz A, Mate Z, Adam E, Oakeley EJ, Schafer E, Nagy F: Genome-wide analysis of gene expression reveals function of the bZIP transcription factor HY5 in the UV-B response of Arabidopsis. Proc Natl Acad Sci USA. 2004, 101: 1397-1402. 10.1073/pnas.0308044100.
Tepperman JM, Zhu T, Chang HS, Wang X, Quail PH: Multiple transcription-factor genes are early targets of phytochrome A signaling. Proc Natl Acad Sci USA. 2001, 98: 9437-9442. 10.1073/pnas.161300998.
Conde C, Agasse A, Glissant D, Tavares R, Geros H, Delrot S: Pathways of glucose regulation of monosaccharide transport in grape cells. Plant Physiol. 2006, 141: 1563-1577. 10.1104/pp.106.080804.
Fillion L, Ageorges A, Picaud S, Coutos-Thevenot P, Lemoine R, Romieu C, Delrot S: Cloning and expression of a hexose transporter gene expressed during the ripening of grape berry. Plant Physiol. 1999, 120: 1083-1094. 10.1104/pp.120.4.1083.
Vignault C, Vachaud M, Cakir B, Glissant D, Dedaldechamp F, Buttner M, Atanassova R, Fleurat-Lessard P, Lemoine R, Delrot S: VvHT1 encodes a monosaccharide transporter expressed in the conducting complex of the grape berry phloem. J Exp Bot. 2005, 56: 1409-1418. 10.1093/jxb/eri142.
Goto-Yamamoto N, Wan GH, Masaki K, Kobayashi S: Structure and transcription of three chalcone synthase genes of grapevine (Vitis vinifera). Plant Sci. 2002, 162: 867-872. 10.1016/S0168-9452(02)00042-0.
Harmer SL, Hogenesch JB, Straume M, Chang HS, Han B, Zhu T, Wang X, Kreps JA, Kay SA: Orchestrated transcription of key pathways in Arabidopsis by the circadian clock. Science. 2000, 290: 2110-2113. 10.1126/science.290.5499.2110.
Moser C, Gatto P, Moser M, Pindo M, Velasco R: Isolation of functional RNA from small amounts of different grape and apple tissues. Mol Biotechnol. 2004, 26: 95-100. 10.1385/MB:26:2:95.
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5: R80-10.1186/gb-2004-5-10-r80.
Ihaka R, Gentleman R: R: a language for data analysis and graphics. J Comput Graph Stat. 1996, 5: 299-314. 10.2307/1390807.
Wu Z, Irizarry RA: Preprocessing of oligonucleotide array data. Nat Biotechno. 2004, 22: 656-658. 10.1038/nbt0604-656b. author reply 658
TM4: Microarray software suite. [http://www.tm4.org/]
Affymetrix home site. [http://www.affymetrix.com/index.affx]
Nadershahi A, Fahrenkrug SC, Ellis LB: Comparison of computational methods for identifying translation initiation sites in EST data. BMC Bioinformatics. 2004, 5: 14-. 10.1186/1471-2105-5-14.
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.
Ramakers C, Ruijter JM, Deprez RH, Moorman AF: Assumption-free analysis of quantitative real-time polymerase chain reaction (PCR) data. Neurosci Lett. 2003, 339: 62-66. 10.1016/S0304-3940(02)01423-4.
Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29: e45-10.1093/nar/29.9.e45.
ArrayExpress home. [http://www.ebi.ac.uk/microarray-as/aer/?#ae-main]
We wish to thank Dr. Davide Marchignoli for informatic support and stimulating discussion, Dr. Massimo Bertamini, Dr. Roberto Zorer for physiological data on berry samples and seasonal climate data (IASMA Research Center), Dr. Mauro Dallaserra (FBK & CNR-Istituto di Biofisica, Trento-Italy) for H2O2 measurements, Dr. Etti Or (Volcani Center, Israel) and Dr. David Neale (University of California, Davis – USA) for critical reading of the manuscript and Dr. Enrico Blanzieri for help with data analysis.
This work was supported in part by grants from the Fondo Unico of the Provincia Autonoma di Trento, post doc project "Profiles" (S.P.), and from the Advanced Biology Project funded by the Fondazione delle Casse di Risparmio di Trento e Rovereto.
SP has made substantial contribution to conception, data analysis and manuscript drafting. MP participated in the GO annotation of the sequences, performed profiles clustering and revised critically the manuscript. AM has contributed to data analysis and revised critically the manuscript. AC, LD and PF contributed to sequence annotation. ADR carried out samples collection, RNA isolation, the biochemical assays and the real time PCR experiments. RV contributed to critically revise the manuscript. RV participated to project's design and coordination. CM has made substantial contribution to conception, project coordination and helped in drafting the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 3: Core set of modulated genes during berry ripening. Log2 average expression data for 2003, 2005 and 2006 of the 1477 modulated genes during berry ripening. Some more information such as the corresponding TC or singleton of the TIGR gene index v 4.0,, the expression profile cluster, the sequence annotation and the corresponding GO functional category are also included. (XLS 606 KB)
Additional file 7: Complete list of genes belonging to the functional categories discussed in the article. Genes of the core set involved in the regulation of berry development (part A) and responsible for berry phenotypic traits (part B), ordered by major GO functional categories (in red). For clarity, transcripts discussed in the manuscript were further classified according to GO more specific terms. (XLS 212 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.