Circadian oscillatory transcriptional programs in grapevine ripening fruits
© Carbonell-Bejerano et al.; licensee BioMed Central Ltd. 2014
Received: 28 December 2013
Accepted: 20 March 2014
Published: 25 March 2014
Temperature and solar radiation influence Vitis vinifera L. berry ripening. Both environmental conditions fluctuate cyclically on a daily period basis and the strength of this fluctuation affects grape ripening too. Additionally, a molecular circadian clock regulates daily cyclic expression in a large proportion of the plant transcriptome modulating multiple developmental processes in diverse plant organs and developmental phases. Circadian cycling of fruit transcriptomes has not been characterized in detail despite their putative relevance in the final composition of the fruit. Thus, in this study, gene expression throughout 24 h periods in pre-ripe berries of Tempranillo and Verdejo grapevine cultivars was followed to determine whether different ripening transcriptional programs are activated during certain times of day in different grape tissues and genotypes.
Microarray analyses identified oscillatory transcriptional profiles following circadian variations in the photocycle and the thermocycle. A higher number of expression oscillating transcripts were detected in samples carrying exocarp tissue including biotic stress-responsive transcripts activated around dawn. Thermotolerance-like responses and regulation of circadian clock-related genes were observed in all studied samples. Indeed, homologs of core clock genes were identified in the grapevine genome and, among them, VvREVEILLE1 (VvRVE1), showed a consistent circadian expression rhythm in every grape berry tissue analysed. Light signalling components and terpenoid biosynthetic transcripts were specifically induced during the daytime in Verdejo, a cultivar bearing white-skinned and aromatic berries, whereas transcripts involved in phenylpropanoid biosynthesis were more prominently regulated in Tempranillo, a cultivar bearing black-skinned berries.
The transcriptome of ripening fruits varies in response to daily environmental changes, which might partially be under the control of circadian clock components. Certain cultivar and berry tissue features could rely on specific circadian oscillatory expression profiles. These findings may help to a better understanding of the progress of berry ripening in short term time scales.
The grapevine fruit is a characteristic berry consisting of an external skin surrounding a fleshy pulp that encloses seeds. Grape biochemical composition is crucial for the different uses of grapes such as winemaking, production of juice and liqueur, fresh consumption or elaboration of raisins, and the final composition is mostly achieved during the ripening phase . Grape ripening is triggered once the seeds have developed and radically changes fruit features from frugivore-repulsive to -attractive. This shift comprises sugar accumulation in the vacuoles of mesocarp cells accompanied by organic acid metabolisation and titratable acidity reduction . Mechanisms to protect seeds from biotic and abiotic stress sources are also activated, mainly in the berry skin . Phenolic compounds are accumulated, including phytoalexins, damaging light-absorbing compounds and animal-attractive anthocyanin pigments . Changes in terpenoid composition result in attractive aromatic profiles, while other compounds of the same family accumulate to accomplish protective functions [5, 6].
Although grapevine ripening is regulated by an intrinsic program that is partially triggered by hormonal signals , the process is also strongly modulated by external factors that influence the final berry composition and its commercial quality. Temperature is a major factor altering grape ripening with high temperatures hastening organic acid metabolisation and inhibiting anthocyanin accumulation [8, 9]. Irradiation intensity and quality perceived in the berry skin produces changes in secondary metabolism. Light promotes flavonols and terpenoids accumulation with some effects being more specifically related with ultraviolet radiation, which is also able to enhance stilbenoids accumulation [6, 10–13]. Moreover, temperature and light conditions cyclically fluctuate under field environments due to the Earth's rotation and, in fact, the strength of their daily oscillation has been shown to affect grape ripening. For instance, reduction in temperature fluctuation intensity hastens berry ripening and alters flavonoid partitioning  and high temperatures applied during night-times are also able to reduce anthocyanin accumulation . In contrast, light pulses applied during the night-time enhance anthocyanin accumulation . These effects of environmental variation on grape ripening rely, at least in part, on changes at the level of gene expression [17–21]. Thus, it can be hypothesised that daily oscillations in factors such as temperature, light or humidity could influence grape ripening progression and consequently on its final composition by means of conditioning circadian fluctuations on the grape transcriptome.
An internal molecular clock is another daily cycling element modulating plant physiology. Core clock components are transcription factors that reciprocally regulate their expression resulting in characteristic circadian expression profiles. Altogether, they determine daily rhythms of expression in a great proportion of plant transcriptomes [22–24]. In addition, the circadian clock integrates inputs from fluctuating conditions like light and temperature, which allows for plants to discriminate daily and seasonal changes [25, 26]. In this way, the circadian clock together with environmental factors control several plant physiological and developmental outputs to be triggered at the proper time of day or season such as photosynthesis, starch and nitrogen metabolism, growth, cold acclimation, bud dormancy, flowering time, tuber formation or stomata and organ movements [27–29]. Since grape ripening is modulated by light and temperature inputs, it is possible that part of this modulation could be integrated by circadian clock components as shown in other plant developmental processes. However, the role of circadian clock on fruit ripening control has scarcely been considered before .
In this study, NimbleGen grapevine whole genome microarrays were used to follow transcriptome variations in ripening grapes throughout a 24 h daily cycle. Tempranillo and Verdejo cultivars were analysed as genotypes appreciated in Spanish viticulture for the production of red and white quality wines, respectively. Verdejo is a cultivar with different berry features than Tempranillo bearing white skin colour and a characteristic aromatic profile, giving rise to wines distinguishable from other white wines [31, 32]. Given that berry flesh and skin ripening transcriptional programs significantly differ , and assuming that this divergence could be more pronounced in cultivars bearing black-skinned berries accordingly with the presumably higher gene expression activity related to the phenylpropanoid pathway in its berry skin, both pericarp tissues were separately analysed in Tempranillo to be more precise. This experiment was carried out using potted plants grown under controlled greenhouse conditions that allowed for identifying tissue-shared as well as tissue-specific cyclic expression profiles under monitored variation of environmental parameters. Part of this response was confirmed in a different experiment in Verdejo berries grown under field environmental conditions despite divergences in genotype, berry colour, berry total soluble solids (TSS) or sampling day photocycles and thermocycles.
Circadian oscillatory gene expression profiles in Tempranillo berries grown under controlled conditions
In order to search for those transcripts showing similar expression fluctuation in both pericarp tissues, a similar analysis (5% FDR in a 6-class Limma and 2-fold change) was run considering skin and flesh samples in the same time point as replicates. In that way, 224 DEG were identified, however, this analysis could underestimate transcripts similarly oscillating in expression in both tissues but involving differences in absolute expression levels. This result indicates that a great proportion of daily oscillating transcripts were shared between both pericarp tissues (Additional file 2). In fact, four expression profiles were identified for pericarp DEG (SF1 to SF4) that resembled those identified in individual tissues and shared with them enriched functional categories (Additional file 4). Cluster SF1 grouping transcripts up-regulated at 12:45 pm in both pericarp tissues was enriched in ‘HSP-mediated protein folding’ and ‘HSF family TF’ and moreover in ‘Vacuolar malate transport’ due to the presence of two transcripts putatively coding for tonoplast dicarboxylate transporters (VIT_00s0187g00130 and VIT_00s2188g00010 in Additional file 2). On the other hand, cluster SF2 grouping transcripts down-regulated at 12:45 pm in both pericarp tissues was enriched in ‘JAZ family TF’; while transcripts in cluster SF4, up-regulated in the pericarp during the night-time, included an over-representation of ‘Circadian clock signalling pathway’ as well as in KIP1 SANTA, GRAS and Constans-like families of TFs (Additional files 3 and 4). Thus, a considerable number of transcripts showed similar circadian expression oscillation in grape skin and flesh, including responses to temperature and light cycles as well as circadian clock-related genes.
Transcripts differentially oscillating in expression between Tempranillo berry flesh and skin
Interestingly, the group of transcripts up-regulated during the daytime only in the flesh (cluster SvsF5) was enriched in ‘Secondary metabolism’ classified transcripts, mainly related to the presence of flavonoid biosynthetic transcripts including one phenylalanine ammonia-lyase (VIT_13s0019g04460), three chalcone synthase (VIT_14s0068g00920, VIT_16s0022g01140 and VIT_16s0022g01190), 13 flavonoid hydroxylases (VIT_06s0009g02810, VIT_06s0009g02830, VIT_06s0009g02840, VIT_06s0009g02860, VIT_06s0009g02880, VIT_06s0009g02920, VIT_06s0009g02970, VIT_06s0009g03010, VIT_06s0009g03040, VIT_06s0009g03050, VIT_06s0009g03110, VIT_09s0002g01090 and VIT_16s0098g00860) and two flavonol synthases (VIT_02s0012g00410 and VIT_13s0047g00210). Other phenylpropanoid biosynthetic genes and two monoterpenoid biosynthetic transcripts coding for 10-geraniol hydroxylase (VIT_02s0012g02820 and VIT_03s0097g00460) were also found in cluster SvsF5, together with VvMybA1 (VIT_02s0033g00380) and VvMybA3 (VIT_02s0033g00450) TFs; the former being a key regulator promoting anthocyanin biosynthesis in grapevine [33–35]. Some of these transcripts underwent a net increase in expression in the flesh during the studied 24 h period (Figure 3).
‘Stilbenoid biosynthesis’ was prominently over-represented in both clusters that involved profiles with expression peaking in the flesh at the early morning. This was due to the presence of a set of 21 stilbenoid biosynthetic transcripts among transcripts more specifically up-regulated in the flesh (cluster SvsF6) and another set of 21 transcripts within the group of transcripts that peaked earlier in the skin at the end of the night-time (cluster SvsF7). In this manner, 42 out of 45 stilbenoid biosynthetic genes annotated in the microarray were differentially involved in these two expression profiles and all of them were located in two gene clusters at chromosomes 10 and 16 (Additional file 5), which suggest a functional diversification of grapevine stilbene synthase (STS) promoters before the extensive duplication of these genes. Cluster SvsF6 was also enriched in ‘Nitrogen metabolism’ and ‘Auxin signalling’ resulting both from the presence of seven nitrilase encoding transcripts. An IAA31-like (VIT_05s0020g01070) and an auxin-responsive D22 (VIT_05s0020g04680) encoding transcript were also included within this profile. The enrichment of ‘Oxidative stress response’ in cluster SvsF6 was determined by the presence of six laccases, one lactoylglutathione lyase and one glutathione S-transferase. Concerning cluster SvsF7, it was enriched in ‘WRKY family TF’ and ‘Ammonium transport’ and included other phenylpropanoid biosynthetic transcripts in addition to STSs and the VvMYB14 TF (VIT_07s0005g03340) (Additional file 5). Finally, cluster SvsF8 was enriched in transcripts belonging to the ‘Oxidative stress response’ category. Therefore, circadian gene expression patterns in the flesh and the skin of pre-ripe berries include both tissue-shared and tissue-specific oscillatory genes that in several cases might be related to specific tissue functions.
Circadian oscillatory gene expression profiles in Verdejo berries grown under field conditions
A total of 778 genes significantly oscillating in expression were identified in Verdejo berry pericarps (5% FDR in Limma and 2-fold change; Additional file 7). Almost one third of Verdejo DEG (240 transcripts) was also significant in any of the comparisons carried out in Tempranillo (Additional file 8) indicating that an important part of this response was shared between experiments regardless of genotypic, environmental and developmental differences.
When Verdejo DEG were clustered, only three main expression profile groups were detected because profiles of transcripts up-regulated around noontime could not be distinguished from those up-regulated during the daytime contrasting with the Tempranillo experiment. Such discrepancy could result from differential oscillation in environmental conditions as berries experienced the maximum temperature near noontime in the Tempranillo experiment and close to the end of the daytime in the Verdejo one (Figures 1 and 4B). Moreover, putative responses to night conditions were maintained until the early morning in Verdejo pericarp as observed for Tempranillo berry flesh (Figure 4B). The number of transcripts up-regulated during the daytime in Verdejo was almost half to that of transcripts up-regulated during the night-time. One expression profile included all daytime up-regulated transcripts (cluster V1). Similarly to Tempranillo transcripts up-regulated at 12:45 pm, cluster V1 was enriched in ‘HSP-mediated protein folding’; however it also showed enrichment in ‘ERF subfamily TF’ and ‘Plant photosystem I supercomplex’ among other categories (Additional file 9). Noticeably, cluster V1 included a HY5-like bZIP TF (VIT_04s0008g05210, Additional file 7). Concerning night up-regulated transcripts in Verdejo pericarp, the absence of light in the first time point of this experiment made it possible to discriminate between transcripts putatively up-regulated by low temperatures (cluster V2) or by darkness (cluster V3). Functional categories over-represented among transcripts up-regulated in Tempranillo berries during the night-time were also enriched in these two Verdejo clusters such as ‘GRAS family TF’ in cluster V2 and ‘Constans-like family TF’ in cluster V3. Additionally, cluster V2 was enriched in ‘Biotic stress response’ and ‘ABC multidrug transport’. Two putative 1-deoxy-D-xylulose 5-phosphate synthase (DXS) encoding transcripts were present within cluster V2 (VIT_11s0052g01730 and VIT_11s0052g01780), while the first one was also induced during the night-time in Tempranillo berry skin (cluster SvsF2) together with another putative DXS (VIT_09s0002g02050). Verdejo cluster V3 was enriched in ‘Xyloglucan modification’ because it contained eight xyloglucan endotransglucosylase (XET) encoding transcripts (Additional file 7). ‘Cytokinin-mediated signalling pathway’ was also over-represented in cluster V3, which included three transcripts coding for authentic-response regulator (ARR) type-A (VIT_08s0007g05390, VIT_13s0067g03430 and VIT_13s0067g03480) and one for histidine kinase AHK2-like cytokinin receptor (VIT_12s0057g00690), while ‘Response regulator family TF’ was enriched in both V2 and V3 clusters of night-time up-regulated transcripts also due to the presence of several ARR genes. Noticeably, VvGAI1 (VIT_01s0011g05260), encoding a repressor of gibberellin responses , was also grouped within the V3 cluster.
Although differential expression profiles between Tempranillo and Verdejo could not be clearly attributable to differences in genotype or growth conditions, time series in both cultivars were explicitly compared. Considering that Tempranillo berry skin and flesh shared circadian oscillatory gene expression (Additional file 2) and that a considerable portion of the ripening transcriptional program is common to both pericarp tissues , Tempranillo skin and flesh samples on each time point were considered as replicates to be compared to Verdejo pericarp samples. In this manner, 286 DEG were identified (5% FDR in a maSigPro and 2-fold change, Additional file 10). More transcripts specifically oscillating in expression were identified in Verdejo than in Tempranillo. However, they were enriched in similar functional categories as those regulated in Tempranillo. They involved ‘HSP-mediated protein folding’ and ‘Response regulator family TF’ among transcripts up-regulated only in Verdejo during the daytime or during the night-time, respectively (clusters TvsV3 and 4 in Additional files 11 and 12). Apart from these similarities, differences in secondary metabolism were observed. Three caffeoyl-CoA O-methyltransferase (VIT_01s0010g03470, VIT_01s0010g03490 and VIT_01s0010g03510) and one secoisolariciresinol dehydrogenase (VIT_06s0004g07540) phenylpropanoid biosynthetic transcripts were up-regulated during the daytime only in Tempranillo (cluster TvsV1). Another secoisolariciresinol dehydrogenase (VIT_08s0058g00790) and a flavonoid hydroxylase (VIT_18s0001g11520) were induced during the night-time specifically in Tempranillo (cluster TvsV2). In contrast, terpenoid biosynthetic transcripts specifically oscillated in expression in Verdejo including a linalool synthase (VIT_00s0271g00060), a myrcene synthase (VIT_00s0271g00030), two β-amyrin synthase (VIT_10s0003g03520 and VIT_10s0003g03650) and three lupeol synthase (VIT_10s0003g03600, VIT_10s0003g03610 and VIT_10s0003g03620) that were up-regulated during the daytime; while an isoprene synthase (VIT_12s0134g00020), a β-ocimene synthase (VIT_12s0134g00030) and a germacrene synthase (VIT_18s0001g04560) were up-regulated during the night-time. Interestingly, VvMYBF1 (=VvMYB12, VIT_07s0005g01210), a light inducible TF [19, 37], was up-regulated during the daytime only in the Verdejo experiment (Additional file 10). These results indicate that circadian expression changes in grapes are mostly consistent among different genotypes. Divergences putatively related with genotypic or environmental differences were also found.
Circadian clock homologous genes oscillating in expression in grapevine berries
Grapevine closest homologues to Arabidopsis core clock proteins
Query cover (%)
Significant in clusterb
CCA1 (AT2G46830, 608 aa)
S4, SF4, V2
LHY (AT1G01060, 645 aa)
S4, SF4, V2
RVE1 (AT5G17300, 387 aa)
S4, SF4, V2
PRR7 (AT5G02810, 727 aa)
PRR9 (AT2G46790, 468 aa)
PRR5 ( AT5G24470, 667 aa)
TOC1 ( AT5G61380, 618 aa)
LUX (AT3G46640, 323 aa)
GI (AT1G22770, 1173 aa)
ELF3 (AT2G25930, 695 aa)
ELF4 (AT2G40080, 111 aa)
Berry transcriptome plasticity in the adaptation to daily environmental changes
Plants as sessile organisms have evolved physiological and developmental adaptations to changes in the surrounding environment. Grape transcriptome changes promoted by the soil composition and weather conditions of different years and locations have been described as examples of plasticity displayed by the berry ripening transcriptome . Herein changes are shown in the berry ripening transcriptome taking place within the same plants as a response to sub-daily fluctuations in the surrounding environment. After applying consistent statistical analyses and cut-offs, a total of 1433 and 868 transcripts showing circadian oscillation in expression were respectively identified in Tempranillo and Verdejo experiments (Additional files 2, 5, 7 and 10). Although these daily oscillations could involve a smaller proportion of circadian cycling genes than in other species or tissues [23, 41], they are larger than those transcriptome changes observed between consecutive berry ripening stages near maturity, consistently described to be low under diverse experimental conditions and genotypes [3, 42, 43]. Thus, sub-daily variation of gene expression in the berry could be more decisive for the final berry composition than other developmental cues at late ripening stages.
Most detected cycling profiles were apparently directed to adapt berry cells to changing environmental conditions such as temperature and light. They included the over-expression of HSPs, ERF TFs and BAG5-like Bcl-2 associated athanogene (VIT_01s0146g00150) positively correlating with high temperatures as well as the down-regulation of ABC transporters expression negatively correlating with them as described for thermotolerance responses in Muscat Hamburg berries . The identification of a higher number of HSPs specifically induced during the daytime in Verdejo berries, which experienced more extreme temperatures than the Tempranillo ones, also suggests the existence of plasticity in berry thermotolerance responses. Nonetheless, some of these as well as other circadian transcriptional responses observed herein may not result in physiological responses as they might not always involve changes at the protein level .
On the other hand, transcriptional changes related to biotic stress responses were activated during the night-time and more greatly in samples including berry skin. Induction of pathogen R-gene-mediated resistance before dawn, the time when the pathogen normally disperses the spores, has been described as a mechanism controlled by core clock genes in Arabidopsis to anticipate infection [45, 46]. An analogous mechanism in the grape berry, an organ particularly susceptible to pathogen attacks when approaching ripeness, could have been evolutionarily selected to ensure seeds dispersal and species survival. TFs of the WRKY family were also induced around dawn in all analysed datasets suggesting that they could participate in the regulation of this defence response [47, 48]. STSs were simultaneously induced only in Tempranillo skin and flesh. Higher expression of STS transcripts has also been correlated with low temperatures in Corvina berries grown under different environmental conditions . Circadian regulation of STS genes could be partially controlled by VvMYB14, which was also induced at dawn only in Tempranillo berry skin, since this TF is able to activate STS promoters in grapevine . In contrast, activation of defence pathways against virus at warmest times of day could be characteristic of grape berries in view that RNA silencing ARGONAUTE2 homologues (VIT_10s0042g01150 and VIT_10s0042g01180) were induced around noontime in all experiments. In fact, plant RNA silencing-mediated defence pathways are generally less active under low temperatures and more active at high temperatures concomitantly with plant virus activity [50, 51].
CONSTANS-like genes participate in the integration of light and circadian clock signalling [25, 52], whereas transcripts within this family consistently peaked around dawn both in Tempranillo and Verdejo berries. Relative to circadian clock signalling, CCA1 and LHY are considered the core morning loop genes in the Arabidopsis molecular clock [53–55]. VvLHY was the only ortholog to these genes found in the grapevine reference genome and did not oscillated in expression in our experiments (Table 1). In contrast, in all analysed grape samples, VvRVE1 showed a cyclic expression profile paralleling these of LHY/CCA1 genes in Arabidopsis leaves. This gene is the ortholog of RVE1, another MYB CCA1 subfamily circadian clock-related gene that, like LHY and CCA1, is able to bind to the ‘evening element’ in the promoter of evening-phased circadian rhythmic genes [56, 57]. Taking into account these coincidences, it might be interesting to test whether VvRVE1 could act as a core clock morning gene in grapevine pre-ripe fruits. Absence of significant regulation in other central clock homologs in the grape analyses may suggest that a simplified version of the clock would operate in grapevine late ripening fruits as shown in other non-leaf Arabidopsis organs like roots . The biology of grapevine fruits is extremely different to that of plant tissues and species where the circadian clock signalling has been explained. Therefore, further research is required to confirm how distinct can this pathway be in particular grapevine tissues. On the other hand, hormone signalling was enriched within several expression profiles peaking in the morning in Tempranillo and Verdejo berries. Hormonal responses coordinated by light- and circadian clock-mediated signalling were shown important for daily promotion of growth in Arabidopsis . Thus, it might be worth testing whether similar mechanisms could modulate grape growth and ripening processes that are also under the control of phytohormones [7, 60]. For instance, relative to gibberellins signalling, rhythmic expression of DELLA genes peaking around dawn under circadian clock control has been described in Arabidopsis seedlings  and similar expression profile of VvGAI1 was observed in Verdejo pericarp.
Berry composition features affected by daily changing conditions
Final composition of mature grape berries is highly related with activity of secondary metabolism. The differential daily cycling responses between berry tissues and cultivars included greater variations in the expression of phenylpropanoid biosynthetic transcripts in the Tempranillo experiment despite these results could be underestimated in our approach. VvMybA1 and VvMybA3 TFs were co-induced during the daytime only in Tempranillo berry flesh together with a large group of flavonoid biosynthetic transcripts. The hybridization signal for VvMybA1 and VvMybA3 microarray probe sets was constitutively higher in Tempranillo berry skin and in Verdejo pericarp when compared to Tempranillo flesh samples (Additional file 1). This could indicate that cross hybridization between both TF probe sets and transcripts takes place in the microarray because, as known, VvMybA3 should be expressed in Verdejo berries but not VvMybA1, whose expression is not expected in white cultivars [34, 62–65]. Indeed, VvGT1 (VIT_16s0039g02230), a direct target of VvMybA1 TF and considered to encode the anthocyanidin-glycosyltransferase catalysing the limiting step for anthocyanin accumulation in berry skin [33, 35, 66], was constitutively more highly expressed in Tempranillo berry skin than in flesh and very much less in Verdejo pericarp (Additional file 1). Further research would be required to show whether VvMybA1 or VvMybA3 positively regulate flavonoid biosynthesis in grape mesocarp cells in a circadian cycle-dependent manner.
The activation of terpenoid biosynthetic genes during the daytime, including monoterpenoid synthases, concurrently to that of light signalling components took place only in Verdejo white and aromatic berries. These light signalling genes included a homolog of HY5, a key gene involved in light perception transduction in Arabidopsis [67, 68], and VvMYBF1, which is induced by sunlight and promotes flavonol biosynthesis in grapes [19, 20, 37]. Absolute expression of VvMYBF1 and this HY5 homolog (VIT_04s0008g05210) during the daytime was higher in Verdejo than in Tempranillo samples (Additional file 1). Photosynthetic genes activated during the daytime only in Verdejo could also be indicative of greater light responses in these berries as some of them are induced by HY5 and by the circadian clock during the daytime in Arabidopsis . Considering these results together with the fact that scent terpenoids accumulation in grapes is enhanced by light [6, 12], it would be interesting to test whether HY5 and VvMYBF1 TFs regulate light induced terpenoid biosynthesis in grapes. Additionally, the extension of circadian gene expression analysis to other black and white cultivars under different cycles of controlled temperature and irradiation conditions would help to explain if some of these specific responses are related to the experimental conditions, the berry skin colour or to other cultivar-specific features.
Concerning primary metabolism-related ripening processes, breakdown of malic acid in grapes is hastened by high temperatures [8, 17, 70]. Only a few genes related to this process cycled in expression despite great temperature oscillations took place during the sampling time courses. Interestingly, two dicarboxylate vacuolar transporters, with a likely function in organic acids import into the vacuole , were induced at 12:45 pm in Tempranillo berries opposing high temperature effects on malic acid metabolisation. In general, results suggest that circadian changes of expression in grapes at late ripening stages affect more deeply secondary than primary metabolism pathways. Part of these expression changes may have important consequences on the final berry composition including berry features depending either on the plant genotype or on the circadian fluctuations in the surrounding environment.
Our study demonstrates the existence of daily oscillatory changes in the grapevine berry transcriptome at late ripening stages. These expression profiles involved responses to temperature and light fluctuation as well as oscillation of circadian clock components. Identified expression profiles suggest that regulatory genes such as a VvRVE1 circadian clock homologous gene as well as Constans-like, WRKY and GRAS family TFs might promote night-activated processes in all analysed experiments. Evidence of circadian variations in the expression of these elements in pre-ripe grapes was shown for the first time in this study. Certain tissue-specific oscillatory expression profiles were related to berry skin and flesh features. Circadian expression profiles in secondary metabolism-related transcripts may impact on the progression of ripening and the final composition of berries, which could result in part from the cyclic expression of chief regulatory genes like VvMybA1, VvMybA3, VvMYBF1, VvMYB14 or a HY5 homolog. In summary, our results could contribute to understanding the progress of berry ripening over sub-daily time scales in response to circadian environmental changes. These findings could lead to new considerations in models directed to predict the effect of ripening season environmental conditions on the final berry composition.
Experiment using V. vinifera cv. Tempranillo was conducted on three years old potted plants, grafted on rootstock 110 Richter. Soil inside pots (50 · 45 · 40 cm) was composed by a peat:sand:vegetal sand mix (63:25:12) over a 5 cm stones layer at the bottom. Tempranillo plants were housed in an experimental greenhouse from Bodegas Miguel Torres S.A. located in Pacs del Penedès (Catalonia, Northeast Spain). Greenhouse was roof-covered with glass allowing for sunlight transmission and no other light was supplied. Total irradiation across the glass was reduced in 44 ± 6% (mean ± SD; Additional file 13). Pots were covered with white plastic films and dropper irrigated according to plant requirements that were calculated as: week water supply = WU * PET * watering coefficient (WU, water use by plants in the previous week; PET, potential evapotranspiration; watering coefficient, 90%). WU was calculated by using weighing lysimeters. Temperature inside the chamber was regulated by Exafan coolers (Grup Sabater, Argentona, Spain) to be the same as external ambient temperature. Temperature and humidity inside the chamber were registered every 10 min by MCU Clima sensors (Grup Sabater). Vines were pruned during the vegetative phase and shoot apexes were trimmed throughout the ripening season. Three separated blocks of ten plants inside the greenhouse were used as biological replicates for the Tempranillo experiment.
Experiment using V. vinifera cv. Verdejo was conducted with five years old vines from a grapevine cultivars garden collection belonging to Bodegas Matarromera S.L. located in Valbuena de Duero (Valladolid, Castilla y León, central plateau Spain). Verdejo vines grafted on 110 Richter rootstocks were conducted in a driving trellis system, in a row orientated in E-W direction. The Verdejo vineyard was dropper irrigated twice throughout the ripening season (2 L · h−1) and defoliation practises were not carried out. Three blocks of ten plants in the same row separated by other ten plants in between were delimited as biological replicates for the Verdejo experiment. For this experiment, meteorological data from Valbuena de Duero station obtained from SiAR website (http://eportal.magrama.gob.es/websiar/ were considered.
Experimental design and berry sampling
To identify berry gene expression changes taking place during a daily cycle, grapes were harvested at six time points throughout a 24 hour period. Time point selection was aimed at maximizing differential environmental conditions: 1, dusk; 2, midnight; 3, before dawn; 4, at the early morning; 5, at the time of maximum temperature (near noontime) and 6, before dusk of the second day. A pre-maturity stage was selected for the sampling day to analyse expression changes at the late ripening phase. Berries from three different bunches were harvested for each sample. On each cultivar, berries of the same density were selected to homogenize the ripening state among all sampled time points and blocks. To do this, every berry was separated from its cluster by cutting the pedicel and its density was determined by floatability in a NaCl solution series as a non-invasive indication of the internal sugar concentration [72–74]. The density interval showing higher berry abundance just before the 24-hours sampling cycle onset was selected for every sample of the same cultivar. TSS in berries from the selected density interval was measured by a digital refractometer WM-7 (ATAGO, Tokyo, Japan). For the Tempranillo experiment, berries between 110–130 g NaCl · l−1 (corresponding to a mean ± SD TSS between replicates of 19.3 ± 0.4 °Bx) were collected at 2011/08/09, 18:00 and 23:30 as well as at 2011/08/10 4:00, 7:45, 12:45 and 18:00 time points (Figure 1). For the Verdejo experiment berries between 130–150 g NaCl · l−1 (corresponding to a mean ± SD TSS of 21.1 ± 0.1 °Bx) were collected at 2011/09/12, 19:30 and 23:45 as well as at 2011/09/12 4:30, 8:45, 13:45 and 18:00 time points (Figure 4). Verdejo berries were always obtained from clusters collected from the north side of the vines. The first time point differed between experiments since it was just before dusk for Tempranillo and after the beginning of the dark period for Verdejo. In all cases solar time is indicated, which corresponds to two hours less than local official time. Sampling took place during completely sunny days in the Verdejo experiment and with only some high clouds in the second afternoon of the Tempranillo experiment. Density selected berries were immediately rinsed in distilled H2O, frozen in liquid nitrogen and kept at -80°C until further analyses. Time lapses from cluster collection to berry freezing were between 10 to 20 min.
Gene expression analyses
Tempranillo and Verdejo frozen berries were peeled and de-seeded, respectively. About 10–15 berries were used for RNA extraction on each sample. Total RNA was extracted from frozen Tempranillo berry skin and flesh and Verdejo pericarp as described by Reid et al. . DNase digestion of contaminating DNA in the RNA samples was carried out with the RNase-Free DNase Set (QIAGEN). Final RNA purification was carried out using the Spectrum™ Plant Total RNA kit (Sigma-Aldrich) according to standard protocols.
Microarray hybridization and data processing
RNA integrity of each RNA preparation was tested using an Agilent 2100 Bioanalyzer (Agilent technologies). cDNA was synthesized from 10 μg of total RNA using the cDNA Synthesis System Kit (NimbleGen-Roche). The cDNA preparation (1 μg) was amplified and labelled with Cy3-random nonamers using the One-Color Labeling Kit (NimbleGen-Roche). If the quality control was correct, then 4 μg of labelled cDNA were hybridized on a NimbleGen microarray 090818 Vitis exp HX12 (NimbleGen-Roche), representing 29,549 predicted genes on the basis of the grapevine 12× V1 gene prediction version (platform GPL17894 in Gene Expression Omnibus [GEO] database). Hybridization solution (NimbleGen Hybridization kit) was added to each labelled cDNA and hybridization was performed for 16 h at 42°C in a HS 4 Hybridization station (NimbleGen-Roche). Hybridized microarrays were washed with Wash buffer kit (NimbleGen-Roche) and scanned at 532 nm and 2 μm resolution in a DNA Microarray Scanner with Surescan High-Resolution Technology (Agilent technologies).
After evaluation of hybridization quality by experimental metrics report implemented in the NimbleScan software version 2.6 (NimbleGen-Roche), probe set signal values from all the microarray hybridizations were background corrected and normalized together using the robust microarray average (RMA)  in the NimbleScan software, which produces calls file for each sample with normalized expression data for each gene. The full microarray expression data are available on GEO under the accession number GSE52167. A dataset was generated from normalized data including the expression of the 29,549 annotated genes in the 54 analysed samples (Additional file 1).
Identification of differentially expressed transcripts and functional analysis
Expression changes during a daily period were searched for each analysed series (Tempranillo berry skin, Tempranillo berry flesh and Verdejo pericarp). On each case, in order to identify transcripts showing changes of expression in any of six analysed time points, a six-class Limma comparison was carried out in Babelomics suite [77, 78]. A P value of 0.05 after the Benjamini-Hochberg adjustment for multiple testing and a ≥ 2-fold changes between any two stages were applied as significance cut-offs on each time series. The same cut-offs were applied in Limma after considering Tempranillo skin and flesh samples for the same time point as replicates to identify transcripts with consistent oscillation in expression along the pericarp. To directly contrast gene expression oscillations in a daily period between Tempranillo berry tissues, skin and flesh series were compared using a maSigPro two-class time series comparison , which was also conducted in Babelomics suite. Polynomial degree 3 was selected to be able to identify transcripts cyclically oscillating in expression. P values of 0.05 were chosen as significance levels for gene selection after Benjamini-Hochberg correction and for model variable in maSigPro. In addition, a ≥ 2-fold difference between tissues in at least one time point after expression data normalization to the corresponding 18:00 h time point of the second day on each tissue and besides, a ≥ 2-fold change between any two stages in at least one tissue were applied. These fold change cut-offs were aimed at identifying transcripts oscillating in expression in at least one tissue that showed considerable differences in their profiles independently of differences in absolute expression between tissues. The same cut-offs were used in maSigPro to identify transcripts with differential oscillation in expression between the 24 h periods analysed in Tempranillo and Verdejo experiments. In this case, Tempranillo skin and flesh samples for the same time point were considered replicates to pick out transcripts with consistent expression along Tempranillo pericarp that was the tissue analysed in Verdejo.
DEG identified in all three Limma analysis and both maSigPro comparisons were clustered by shared expression profiles after Log2 expression data normalization to the expression in 18:00 h day-two time point of the corresponding tissue and cultivar using self-organizing maps (SOM) analysis . Euclidean squared metrics and scaled rows were selected for the SOM that were run in Acuity 4.0 (Axon Molecular Devices, http://www.moleculardevices.com). For each SOM analysis, the most informative number of clusters was assessed by gap statistical analyses , which were also carried out in Acuity 4.0. Each identified cluster was analysed on Babelomics suite to search for significant functional enrichment following a grapevine specific functional classification of 12× V1 predicted transcripts . Fisher’s exact test was carried out in a FatiGO analysis  to compare each study list to the list of total transcripts housed in the grapevine 12× V1 gene predictions . Significant enrichment was considered in case of P value ≤0.05 after Benjamini and Hochberg correction for multiple testing. Finally lists of significant transcripts identified from different analysis were compared by Venn diagrams performed in Venny (http://bioinfogp.cnb.csic.es/tools/venny/index.html) .
Search of circadian clock gene homologues
The grapevine genomic sequence was searched for loci encoding homologous proteins to Arabidopsis CCA1, LHY, PRR9, PRR7, PRR5, GI, TOC1, LUX, ELF3 and ELF3, PRR5 core clock proteins . A BLAT alignment from each Arabidopsis protein sequence to the PN40024 12X genomic sequence was carried out in the Genoscope website (http://www.genoscope.cns.fr/blat-server/cgi-bin/vitis/webBlat). For each locus, the corresponding protein in the 12× V1 gene annotations version was identified from Grimplet et al. (2012). Grapevine 12× V1 protein sequences were obtained from the Uniprot website (http://www.uniprot.org/) and were aligned to Arabidopsis protein sequences by blastp (http://blast.ncbi.nlm.nih.gov/).
This work was funded, within the CENIT-DEMETER project, by the Spanish Centre for Technological and Industrial Development (CDTI) from the Spanish Ministry of Science and Innovation (MICINN). The authors express their gratitude to all entities participating in the CENIT-Deméter consortium, and especially to Bodegas Miguel Torres S. A. and Bodegas Matarromera S. L. for their support and for providing the plants and the growing conditions. We also thank Anna Velazquez, Meritxell Salamero and Eudald Crivillés from Bodegas Miguel Torres S. A. as well as Sonia Villanueva, Sara Martín de la Helguera and Rubén Gallo from Bodegas Matarromera S.L. for taking care of plants. We are grateful to Dr. Jérôme Grimplet for the supply of grapevine NimbleGen microarray annotations, Drs Mario Pezzotti and Massimo Delledonne for the design of NimbleGen microarrays, Rafael Torres-Pérez for help formatting gene annotations for FatiGO analysis and Cheo Machín for editing of the English style. Microarray hybridizations were carried out at the Genomics Unit of the National Centre for Biotechnology, CNB-CSIC, Madrid, Spain.
- Dai ZW, Leon C, Feil R, Lunn JE, Delrot S, Gomes E: Metabolic profiling reveals coordinated switches in primary carbohydrate metabolism in grape berry (Vitis vinifera L.), a non-climacteric fleshy fruit. J Exp Bot. 2013, 64: 1345-1355. 10.1093/jxb/ers396.PubMed CentralView ArticlePubMedGoogle Scholar
- Conde C, Silva P, Fontes N, Dias ACP, Tavares RM, Sousa MJ, Agasse A, Delrot S, Gerós H: Biochemical changes throughout Grape Berry development and fruit and wine quality. Food. 2007, 1: 1-22.Google Scholar
- Lijavetzky D, Carbonell-Bejerano P, Grimplet J, Bravo G, Flores P, Fenoll J, Hellin P, Oliveros JC, Martinez-Zapater JM: Berry flesh and skin ripening features in Vitis vinifera as assessed by transcriptional profiling. PLoS One. 2012, 7: e39547-10.1371/journal.pone.0039547.PubMed CentralView ArticlePubMedGoogle Scholar
- Teixeira A, Eiras-Dias J, Castellarin SD, Geros H: Berry phenolics of grapevine under challenging environments. Int J Mol Sci. 2013, 14: 18711-18739. 10.3390/ijms140918711.PubMed CentralView ArticlePubMedGoogle Scholar
- Rodriguez A, Alquezar B, Pena L: Fruit aromas in mature fleshy fruits as signals of readiness for predation and seed dispersal. New Phytol. 2013, 197: 36-48. 10.1111/j.1469-8137.2012.04382.x.View ArticlePubMedGoogle Scholar
- Gil M, Bottini R, Berli F, Pontin M, Silva MF, Piccoli P: Volatile organic compounds characterized from grapevine (Vitis vinifera L. cv. Malbec) berries increase at pre-harvest and in response to UV-B radiation. Phytochemistry. 2013, 96: 148-157.View ArticlePubMedGoogle Scholar
- Davies C, Böttcher C: Hormonal Control of Grape Berry Ripening. Grapevine Molecular Physiology & Biotechnology. Edited by: Roubelakis-Angelakis KA. 2009, Netherlands: Springer, 229-261.View ArticleGoogle Scholar
- Kliewer WM: Effect of temperature on the composition of grapes grown under field and controlled conditions. Proc Am Soc Horti Cult Sci. 1968, 93: 797-806.Google Scholar
- Tarara JM, Lee J, Spayd SE, Scagel CF: Berry temperature and solar radiation alter acylation, proportion, and concentration of anthocyanin in Merlot grapes. Am J Enol Vitic. 2008, 59: 235-247.Google Scholar
- Berli F, D’Angelo J, Cavagnaro B, Bottini R, Wuilloud R, Silva MF: Phenolic composition in grape (Vitis vinifera L. cv. Malbec) ripened with different solar UV-B radiation levels by capillary zone electrophoresis. J Agric Food Chem. 2008, 56: 2892-2898. 10.1021/jf073421+.View ArticlePubMedGoogle Scholar
- Spayd SE, Tarara JM, Mee DL, Ferguson JC: Separation of sunlight and temperature effects on the composition of Vitis vinifera cv. Merlot Berries. Am J Enol Vitic. 2002, 53: 171-182.Google Scholar
- Skinkis PA, Bordelon BP, Butz EM: Effects of sunlight exposure on berry and wine monoterpenes and sensory characteristics of Traminette. Am J Enol Vitic. 2010, 61: 147-156.Google Scholar
- Cortell JM, Kennedy JA: Effect of shading on accumulation of flavonoid compounds in (Vitis vinifera L.) pinot noir fruit and extraction in a model system. J Agric Food Chem. 2006, 54: 8510-8520. 10.1021/jf0616560.View ArticlePubMedGoogle Scholar
- Cohen SD, Tarara JM, Kennedy JA: Diurnal temperature range compression Hastens Berry development and modifies flavonoid partitioning in Grapes. Am J Enol Vitic. 2011, 63: 112-120.View ArticleGoogle Scholar
- Kliewer WM, Torres RE: Effect of Controlled Day and Night Temperatures on Grape Coloration. Am J Enol Vitic. 1972, 23: 71-77.Google Scholar
- Azuma A, Ito A, Moriguchi T, Yakushiji H, Kobayashi S: Light emitting diode irradiation at night accelerates anthocyanin accumulation in grape skin. Acta Horticult. 2012, 956: 341-348.View ArticleGoogle Scholar
- Carbonell-Bejerano P, Santa Maria E, Torres-Perez R, Royo C, Lijavetzky D, Bravo G, Aguirreolea J, Sanchez-Diaz M, Antolin MC, Martinez-Zapater JM: Thermotolerance responses in ripening berries of Vitis vinifera L. cv Muscat Hamburg. Plant Cell Physiol. 2013, 54: 1200-1216. 10.1093/pcp/pct071.View ArticlePubMedGoogle Scholar
- Mori K, Goto-Yamamoto N, Kitayama M, Hashizume K: Loss of anthocyanins in red-wine grape under high temperature. J Exp Bot. 2007, 58: 1935-1945. 10.1093/jxb/erm055.View ArticlePubMedGoogle Scholar
- Matus JT, Loyola R, Vega A, Pena-Neira A, Bordeu E, Arce-Johnson P, Alcalde JA: Post-veraison sunlight exposure induces MYB-mediated transcriptional regulation of anthocyanin and flavonol synthesis in berry skins of Vitis vinifera. J Exp Bot. 2009, 60: 853-867. 10.1093/jxb/ern336.PubMed CentralView ArticlePubMedGoogle Scholar
- Azuma A, Yakushiji H, Koshita Y, Kobayashi S: Flavonoid biosynthesis-related genes in grape skin are differentially regulated by temperature and light conditions. Planta. 2012, 236: 1067-1080. 10.1007/s00425-012-1650-x.View ArticlePubMedGoogle Scholar
- Koyama K, Ikeda H, Poudel PR, Goto-Yamamoto N: Light quality affects flavonoid biosynthesis in young berries of Cabernet Sauvignon grape. Phytochemistry. 2012, 78: 54-64.View ArticlePubMedGoogle Scholar
- Harmer SL: The circadian system in higher plants. Annu Rev Plant Biol. 2009, 60: 357-377. 10.1146/annurev.arplant.043008.092054.View ArticlePubMedGoogle Scholar
- Filichkin SA, Breton G, Priest HD, Dharmawardhana P, Jaiswal P, Fox SE, Michael TP, Chory J, Kay SA, Mockler TC: Global profiling of rice and poplar transcriptomes highlights key conserved circadian-controlled pathways and cis-regulatory modules. PLoS One. 2011, 6: e16907-10.1371/journal.pone.0016907.PubMed CentralView ArticlePubMedGoogle Scholar
- Huang W, Perez-Garcia P, Pokhilko A, Millar AJ, Antoshechkin I, Riechmann JL, Mas P: Mapping the core of the Arabidopsis circadian clock defines the network structure of the oscillator. Science. 2012, 336: 75-79. 10.1126/science.1219075.View ArticlePubMedGoogle Scholar
- Yanovsky MJ, Kay SA: Molecular basis of seasonal time measurement in Arabidopsis. Nature. 2002, 419: 308-312. 10.1038/nature00996.View ArticlePubMedGoogle Scholar
- Salome PA, Weigel D, McClung CR: The role of the Arabidopsis morning loop components CCA1, LHY, PRR7, and PRR9 in temperature compensation. Plant Cell. 2010, 22: 3650-3661. 10.1105/tpc.110.079087.PubMed CentralView ArticlePubMedGoogle Scholar
- McWatters HG, Devlin PF: Timing in plants–a rhythmic arrangement. FEBS Lett. 2011, 585: 1474-1484. 10.1016/j.febslet.2011.03.051.View ArticlePubMedGoogle Scholar
- Hotta CT, Gardner MJ, Hubbard KE, Baek SJ, Dalchau N, Suhita D, Dodd AN, Webb AA: Modulation of environmental responses of plants by circadian clocks. Plant Cell Environ. 2007, 30: 333-349. 10.1111/j.1365-3040.2006.01627.x.View ArticlePubMedGoogle Scholar
- Eriksson ME, Webb AA: Plant cell responses to cold are all about timing. Curr Opin Plant Biol. 2011, 14: 731-737. 10.1016/j.pbi.2011.08.005.View ArticlePubMedGoogle Scholar
- Piechulla B, Gruissem W: Diurnal mRNA fluctuations of nuclear and plastid genes in developing tomato fruits. EMBO J. 1987, 6: 3593-3599.PubMed CentralPubMedGoogle Scholar
- Rodríguez-Nogales JM, Fernández-Fernández E, Vila-Crespo J: Characterisation and classification of Spanish Verdejo young white wines by volatile and sensory analysis with chemometric tools. J Sci Food Agric. 2009, 89: 1927-1935. 10.1002/jsfa.3674.View ArticleGoogle Scholar
- Campo E, Do BV, Ferreira V, Valentin D: Aroma properties of young Spanish monovarietal white wines: a study using sorting task, list of terms and frequency of citation. Aust J Grape Wine Res. 2008, 14: 104-115. 10.1111/j.1755-0238.2008.00010.x.View ArticleGoogle Scholar
- Cutanda-Perez MC, Ageorges A, Gomez C, Vialet S, Terrier N, Romieu C, Torregrosa L: Ectopic expression of VlmybA1 in grapevine activates a narrow set of genes involved in anthocyanin synthesis and transport. Plant Mol Biol. 2009, 69: 633-648. 10.1007/s11103-008-9446-x.View ArticlePubMedGoogle Scholar
- Walker AR, Lee E, Bogs J, McDavid DA, Thomas MR, Robinson SP: White grapes arose through the mutation of two similar and adjacent regulatory genes. Plant J. 2007, 49: 772-785. 10.1111/j.1365-313X.2006.02997.x.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Boss PK, Thomas MR: Association of dwarfism and floral induction with a grape ‘green revolution’ mutation. Nature. 2002, 416: 847-850. 10.1038/416847a.View ArticlePubMedGoogle Scholar
- Czemmel S, Stracke R, Weisshaar B, Cordon N, Harris NN, Walker AR, Robinson SP, Bogs J: The grapevine R2R3-MYB transcription factor VvMYBF1 regulates flavonol synthesis in developing grape berries. Plant Physiol. 2009, 151: 1513-1530. 10.1104/pp.109.142059.PubMed CentralView ArticlePubMedGoogle Scholar
- Takata N, Saito S, Tanaka Saito C, Nanjo T, Shinohara K, Uemura M: Molecular phylogeny and expression of poplar circadian clock genes, LHY1 and LHY2. New Phytol. 2009, 181: 808-819. 10.1111/j.1469-8137.2008.02714.x.View ArticlePubMedGoogle Scholar
- Rawat R, Schwartz J, Jones MA, Sairanen I, Cheng Y, Andersson CR, Zhao Y, Ljung K, Harmer SL: REVEILLE1, a Myb-like transcription factor, integrates the circadian clock and auxin pathways. Proc Natl Acad Sci USA. 2009, 106: 16883-16888. 10.1073/pnas.0813035106.PubMed CentralView ArticlePubMedGoogle Scholar
- Dal Santo S, Tornielli GB, Zenoni S, Fasoli M, Farina L, Anesi A, Guzzo F, Delledonne M, Pezzotti M: The plasticity of the grapevine berry transcriptome. Genome Biol. 2013, 14: r54-10.1186/gb-2013-14-6-r54.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Deluc LG, Grimplet J, Wheatley MD, Tillett RL, Quilici DR, Osborne C, Schooley DA, Schlauch KA, Cushman JC, Cramer GR: Transcriptomic and metabolite analyses of Cabernet Sauvignon grape berry development. BMC Genomics. 2007, 8: 429-10.1186/1471-2164-8-429.PubMed CentralView ArticlePubMedGoogle Scholar
- Guillaumie S, Fouquet R, Kappel C, Camps C, Terrier N, Moncomble D, Dunlevy JD, Davies C, Boss PK, Delrot S: Transcriptional analysis of late ripening stages of grapevine berry. BMC Plant Biol. 2011, 11: 165-10.1186/1471-2229-11-165.PubMed CentralView ArticlePubMedGoogle Scholar
- Li QB, Haskell D, Zhang C, Sung DY, Guy C: Diurnal regulation of Hsp70s in leaf tissue. Plant J. 2000, 21: 373-378. 10.1046/j.1365-313x.2000.00673.x.View ArticlePubMedGoogle Scholar
- Wang W, Barnaby JY, Tada Y, Li H, Tor M, Caldelari D, Lee DU, Fu XD, Dong X: Timing of plant immune responses by a central circadian regulator. Nature. 2011, 470: 110-114. 10.1038/nature09766.View ArticlePubMedGoogle Scholar
- Hua J: Modulation of plant immunity by light, circadian rhythm, and temperature. Curr Opin Plant Biol. 2013, 16: 406-413. 10.1016/j.pbi.2013.06.017.View ArticlePubMedGoogle Scholar
- Eulgem T, Rushton PJ, Schmelzer E, Hahlbrock K, Somssich IE: Early nuclear events in plant defence signalling: rapid gene activation by WRKY transcription factors. EMBO J. 1999, 18: 4689-4699. 10.1093/emboj/18.17.4689.PubMed CentralView ArticlePubMedGoogle Scholar
- Marchive C, Mzid R, Deluc L, Barrieu F, Pirrello J, Gauthier A, Corio-Costet MF, Regad F, Cailleteau B, Hamdi S, Lauvergeat V: Isolation and characterization of a Vitis vinifera transcription factor, VvWRKY1, and its effect on responses to fungal pathogens in transgenic tobacco plants. J Exp Bot. 2007, 58: 1999-2010. 10.1093/jxb/erm062.View ArticlePubMedGoogle Scholar
- Holl J, Vannozzi A, Czemmel S, D'Onofrio C, Walker AR, Rausch T, Lucchin M, Boss PK, Dry IB, Bogs J: The R2R3-MYB transcription factors MYB14 and MYB15 regulate stilbene biosynthesis in Vitis vinifera. Plant Cell. 2013, 25: 4135-4149. 10.1105/tpc.113.117127.PubMed CentralView ArticlePubMedGoogle Scholar
- Szittya G, Silhavy D, Molnar A, Havelda Z, Lovas A, Lakatos L, Banfalvi Z, Burgyan J: Low temperature inhibits RNA silencing-mediated defence by the control of siRNA generation. EMBO J. 2003, 22: 633-640. 10.1093/emboj/cdg74.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang X, Singh J, Li D, Qu F: Temperature-dependent survival of Turnip crinkle virus-infected arabidopsis plants relies on an RNA silencing-based defense that requires dcl2, AGO2, and HEN1. J Virol. 2012, 86: 6847-6854. 10.1128/JVI.00497-12.PubMed CentralView ArticlePubMedGoogle Scholar
- Ledger S, Strayer C, Ashton F, Kay SA, Putterill J: Analysis of the function of two circadian-regulated CONSTANS-LIKE genes. Plant J. 2001, 26: 15-22. 10.1046/j.1365-313x.2001.01003.x.View ArticlePubMedGoogle Scholar
- Alabadi D, Oyama T, Yanovsky MJ, Harmon FG, Mas P, Kay SA: Reciprocal regulation between TOC1 and LHY/CCA1 within the Arabidopsis circadian clock. Science. 2001, 293: 880-883. 10.1126/science.1061320.View ArticlePubMedGoogle Scholar
- Mizoguchi T, Wheatley K, Hanzawa Y, Wright L, Mizoguchi M, Song HR, Carre IA, Coupland G: LHY and CCA1 are partially redundant genes required to maintain circadian rhythms in Arabidopsis. Dev Cell. 2002, 2: 629-641. 10.1016/S1534-5807(02)00170-3.View ArticlePubMedGoogle Scholar
- Alabadi D, Yanovsky MJ, Mas P, Harmer SL, Kay SA: Critical role for CCA1 and LHY in maintaining circadian rhythmicity in Arabidopsis. Curr Biol. 2002, 12: 757-761. 10.1016/S0960-9822(02)00815-1.View ArticlePubMedGoogle Scholar
- Harmer SL, Kay SA: Positive and negative factors confer phase-specific circadian regulation of transcription in Arabidopsis. Plant Cell. 2005, 17: 1926-1940. 10.1105/tpc.105.033035.PubMed CentralView ArticlePubMedGoogle Scholar
- Gong W, He K, Covington M, Dinesh-Kumar SP, Snyder M, Harmer SL, Zhu YX, Deng XW: The development of protein microarrays and their applications in DNA-protein and protein-protein interaction analyses of Arabidopsis transcription factors. Mol Plant. 2008, 1: 27-41. 10.1093/mp/ssm009.PubMed CentralView ArticlePubMedGoogle Scholar
- James AB, Monreal JA, Nimmo GA, Kelly CL, Herzyk P, Jenkins GI, Nimmo HG: The circadian clock in Arabidopsis roots is a simplified slave version of the clock in shoots. Science. 2008, 322: 1832-1835. 10.1126/science.1161403.View ArticlePubMedGoogle Scholar
- Michael TP, Breton G, Hazen SP, Priest H, Mockler TC, Kay SA, Chory J: A morning-specific phytohormone gene expression program underlying rhythmic plant growth. PLoS Biol. 2008, 6: e225-10.1371/journal.pbio.0060225.PubMed CentralView ArticlePubMedGoogle Scholar
- Bottcher C, Burbidge CA, Boss PK, Davies C: Interactions between ethylene and auxin are crucial to the control of grape (Vitis vinifera L.) berry ripening. BMC Plant Biol. 2013, 13: 222-10.1186/1471-2229-13-222.PubMed CentralView ArticlePubMedGoogle Scholar
- Arana MV, Marin-de la Rosa N, Maloof JN, Blazquez MA, Alabadi D: Circadian oscillation of gibberellin signaling in Arabidopsis. Proc Natl Acad Sci USA. 2011, 108: 9292-9297. 10.1073/pnas.1101050108.PubMed CentralView ArticlePubMedGoogle Scholar
- Kobayashi S, Goto-Yamamoto N, Hirochika H: Retrotransposon-induced mutations in grape skin color. Science. 2004, 304: 982-10.1126/science.1095011.View ArticlePubMedGoogle Scholar
- Lijavetzky D, Ruiz-Garcia L, Cabezas JA, De Andres MT, Bravo G, Ibanez A, Carreno J, Cabello F, Ibanez J, Martinez-Zapater JM: Molecular genetics of berry colour variation in table grape. Mol Genet Genomics. 2006, 276: 427-435. 10.1007/s00438-006-0149-1.View ArticlePubMedGoogle Scholar
- This P, Lacombe T, Cadle-Davidson M, Owens CL: Wine grape (Vitis vinifera L.) color associates with allelic variation in the domestication gene VvmybA1. Theor Appl Genet. 2007, 114: 723-730. 10.1007/s00122-006-0472-2.View ArticlePubMedGoogle Scholar
- Castellarin SD, Di Gaspero G: Transcriptional control of anthocyanin biosynthetic genes in extreme phenotypes for berry pigmentation of naturally occurring grapevines. BMC Plant Biol. 2007, 7: 46-10.1186/1471-2229-7-46.PubMed CentralView ArticlePubMedGoogle Scholar
- Boss PK, Davies C, Robinson SP: Expression of anthocyanin biosynthesis pathway genes in red and white grapes. Plant Mol Biol. 1996, 32: 565-569. 10.1007/BF00019111.View ArticlePubMedGoogle Scholar
- Osterlund MT, Hardtke CS, Wei N, Deng XW: Targeted destabilization of HY5 during light-regulated development of Arabidopsis. Nature. 2000, 405: 462-466. 10.1038/35013076.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Lee J, He K, Stolc V, Lee H, Figueroa P, Gao Y, Tongprasit W, Zhao H, Lee I, Deng XW: Analysis of transcription factor HY5 genomic binding sites revealed its hierarchical role in light regulation of development. Plant Cell. 2007, 19: 731-749. 10.1105/tpc.106.047688.PubMed CentralView ArticlePubMedGoogle Scholar
- Ruffner HP, Hawker JS, Hale CR: Temperature and enzymic control of malate metabolism in berries of Vitis vinifera. Phytochemistry. 1976, 15: 1877-1880. 10.1016/S0031-9422(00)88835-4.View ArticleGoogle Scholar
- Hurth MA, Suh SJ, Kretzschmar T, Geis T, Bregante M, Gambale F, Martinoia E, Neuhaus HE: Impaired pH homeostasis in Arabidopsis lacking the vacuolar dicarboxylate transporter and analysis of carboxylic acid transport across the tonoplast. Plant Physiol. 2005, 137: 901-910. 10.1104/pp.104.058453.PubMed CentralView ArticlePubMedGoogle Scholar
- Stein U, Blaich R, Wind R: A novel method for non-destructive determination of the sugar content and for classification of grape berries. Vitis. 1983, 22: 15-22.Google Scholar
- Fournand D, Vicens A, Sidhoum L, Souquet JM, Moutounet M, Cheynier V: Accumulation and extractability of grape skin tannins and anthocyanins at different advanced physiological stages. J Agric Food Chem. 2006, 54: 7331-7338. 10.1021/jf061467h.View ArticlePubMedGoogle Scholar
- Rolle L, Segade SR, Torchio F, Giacosa S, Cagnasso E, Marengo F, Gerbi V: Influence of grape density and harvest date on changes in phenolic composition, phenol extractability indices, and instrumental texture properties during ripening. J Agric Food Chem. 2011, 59: 8796-8805. 10.1021/jf201318x.View ArticlePubMedGoogle Scholar
- Reid KE, Olsson N, Schlosser J, Peng F, Lund ST: An optimized grapevine RNA isolation procedure and statistical determination of reference genes for real-time RT-PCR during berry development. BMC Plant Biol. 2006, 6: 27-10.1186/1471-2229-6-27.PubMed CentralView ArticlePubMedGoogle Scholar
- Irizarry RA, Bolstad BM, Collin F, Cope LM, Hobbs B, Speed TP: Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res. 2003, 31: e15-10.1093/nar/gng015.PubMed CentralView ArticlePubMedGoogle Scholar
- Medina I, Carbonell J, Pulido L, Madeira SC, Goetz S, Conesa A, Tarraga J, Pascual-Montano A, Nogales-Cadenas R, Santoyo J, García F, Marbà M, Montaner D, Dopazo J: Babelomics: an integrative platform for the analysis of transcriptomics, proteomics and genomic data with advanced functional profiling. Nucleic Acids Res. 2010, 38: W210-W213. 10.1093/nar/gkq388.PubMed CentralView ArticlePubMedGoogle Scholar
- Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3 (1): 1-25. Article3Google Scholar
- Conesa A, Nueda MJ, Ferrer A, Talon M: maSigPro: a method to identify significantly differential expression profiles in time-course microarray experiments. Bioinformatics. 2006, 22: 1096-1102. 10.1093/bioinformatics/btl056.View ArticlePubMedGoogle Scholar
- Toronen P, Kolehmainen M, Wong G, Castren E: Analysis of gene expression data using self-organizing maps. FEBS Lett. 1999, 451: 142-146. 10.1016/S0014-5793(99)00524-4.View ArticlePubMedGoogle Scholar
- Tibshirani R, Walther G, Hastie T: Estimating the number of clusters in a data set via the gap statistic. J R Stat Soc B. 2001, 63: 411-423. 10.1111/1467-9868.00293.View ArticleGoogle Scholar
- Grimplet J, Van Hemert J, Carbonell-Bejerano P, Diaz-Riquelme J, Dickerson J, Fennell A, Pezzotti M, Martinez-Zapater JM: Comparative analysis of grapevine whole-genome gene predictions, functional annotation, categorization and integration of the predicted gene sequences. BMC Res Notes. 2012, 5: 213-10.1186/1756-0500-5-213.PubMed CentralView ArticlePubMedGoogle Scholar
- Al-Shahrour F, Diaz-Uriarte R, Dopazo J: FatiGO: a web tool for finding significant associations of Gene Ontology terms with groups of genes. Bioinformatics. 2004, 20: 578-580. 10.1093/bioinformatics/btg455.View ArticlePubMedGoogle Scholar
- VENNY. An interactive tool for comparing lists with Venn Diagrams. [http://bioinfogp.cnb.csic.es/tools/venny/index.html]
- Nagel DH, Kay SA: Complexity in the wiring and regulation of plant circadian networks. Curr Biol. 2012, 22: R648-R657. 10.1016/j.cub.2012.07.025.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.