Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

Open Access Research article

Increasing the source/sink ratio in Vitis vinifera (cv Sangiovese) induces extensive transcriptome reprogramming and modifies berry ripening

Chiara Pastore1, Sara Zenoni2, Giovanni Battista Tornielli2, Gianluca Allegro1, Silvia Dal Santo2, Gabriele Valentini1, Cesare Intrieri1, Mario Pezzotti2* and Ilaria Filippetti1

Author affiliations

1 Department of Fruit Tree and Woody Plant Science, University of Bologna, Viale Fanin, 46, 40126, Bologna, Italy

2 Department of Biotechnology, University of Verona, Strada Le Grazie 15, 37134, Verona, Italy

For all author emails, please log on.

Citation and License

BMC Genomics 2011, 12:631  doi:10.1186/1471-2164-12-631

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2164/12/631


Received:24 August 2011
Accepted:23 December 2011
Published:23 December 2011

© 2011 Pastore et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Cluster thinning is an agronomic practice in which a proportion of berry clusters are removed from the vine to increase the source/sink ratio and improve the quality of the remaining berries. Until now no transcriptomic data have been reported describing the mechanisms that underlie the agronomic and biochemical effects of thinning.

Results

We profiled the transcriptome of Vitis vinifera cv. Sangiovese berries before and after thinning at veraison using a genome-wide microarray representing all grapevine genes listed in the latest V1 gene prediction. Thinning increased the source/sink ratio from 0.6 to 1.2 m2 leaf area per kg of berries and boosted the sugar and anthocyanin content at harvest. Extensive transcriptome remodeling was observed in thinned vines 2 weeks after thinning and at ripening. This included the enhanced modulation of genes that are normally regulated during berry development and the induction of a large set of genes that are not usually expressed.

Conclusion

Cluster thinning has a profound effect on several important cellular processes and metabolic pathways including carbohydrate metabolism and the synthesis and transport of secondary products. The integrated agronomic, biochemical and transcriptomic data revealed that the positive impact of cluster thinning on final berry composition reflects a much more complex outcome than simply enhancing the normal ripening process.

Background

Many agronomic practices are employed to maximize grape berry quality in the highly competitive wine industry, including the control of bud load during winter pruning and cluster thinning during berry development. Cluster thinning acts directly to increase the source/sink balance of grapevine plants and the technique is used to prevent overcropping in varieties characterized by excessive bud fertility or in areas where reduced yield is a prerequisite for high-quality wine production. Under such conditions, cluster thinning is performed to obtain a leaf area/yield ratio of 0.8-1.2 m2/kg. Below this threshold value several authors [1,2] reported a positive correlation between berry juice soluble sugars and leaf area/crop weight ratio. Similar results have been achieved following cluster thinning in different varieties and sites [3-7]. Thinning also increases the anthocyanin content of berries, which is an important quality determinant of red wines [8-10]. The anthocyanin composition is also affected, e.g. cluster thinning induced the accumulation of 3',4'-substituted anthocyanins in Sangiovese and Nebbiolo varieties [3,5].

Although it is well known that many different factors influence flavonoid and anthocyanin biosynthesis [11], sugar is likely to play a prominent role because of the concomitant accumulation of soluble solids in the berry flesh and anthocyanins in the skin of red grape varieties. This relationship was first proposed by Pirie and Mullins [12], who suggested that the sugar content of red berry flesh could regulate anthocyanin production, and this was supported by in vitro experiments showing an increase in phenylalanine ammonia-lyase (PAL) activity [13]. Increased anthocyanin accumulation after treatment with sucrose and other sugars has been already demonstrated in grapevine and in a variety of other plant species [14-17]. Several sugars have been shown to induce genes encoding enzymes in the anthocyanin biosynthesis pathway, such as chalcone synthase (CHS), dihydroflavonol reductase (DFR), leucoanthocyanidin dioxygenase (LDOX) [18,19] and flavanone 3-hydroxylase (F3H) [20]. Sucrose boxes have been identified in the promoters of some of these genes [18,19].

Little is currently known about the regulation of gene expression when the source/sink ratio is deliberately altered in the field. We report the first transcriptomic analysis (integrated with agronomic and biochemical data) aiming to determine the mechanisms that control Sangiovese berry composition by comparing gene expression profiles of thinned and control vines. Berry transcriptional profiles were analyzed during ripening using the most comprehensive grapevine microarray available to date, representing 29,549 genes from the most recent 12X grapevine V1 gene prediction http://srs.ebi.ac.uk/ webcite. We observed substantial transcriptomic remodeling in berries from thinned vines which became evident by 2 weeks post-treatment and persisted until ripening, with particular impact on genes involved in carbohydrate metabolism, flavonoid biosynthesis and transport. The results from these studies provide insight into the molecular basis of berry ripening induced by vineyard management techniques.

Results

The effect of cluster thinning on yield and berry ripening

Cluster thinning (CT) was carried out to remove approximately 50% of the bunches on each vine, leaving approximately eight clusters per vine in comparison with 16 on control (C) plants, thus reducing the yield by ~54% (Table 1). The average bunch and berry weight remained the same in CT and C plants (Table 1).

Table 1. Influence of cluster thinning (CT) on yield component and berry composition at harvest.

The leaf area per vine was similar in CT and C plants at harvest, indicating that cluster thinning increased the leaf area/yield ratio from 0.6 m2/kg in C plants to 1.2 m2/kg in CT plants (Table 1). CT berries also accumulated more total soluble solids than C berries from full veraison (8 d after cluster thinning) until harvest, and were less acidic and had higher °Brix values than controls at harvest (Figure 1A and Table 1).

thumbnailFigure 1. Agronomic effects of cluster thinning. (A) Accumulation of soluble sugars (°Brix) and titratable acidity (TA) in berries from C and CT vines (n = 3). (B) Anthocyanin concentration during ripening and (C) anthocyanin composition at harvest in berry skins from C and CT vines (n = 3). Bars represent ± SE. Asterisks indicate significant differences between the treatments at the same date using the Student-Newman-Keuls test (*P < 0.05). Gray background represents the veraison phase.

The anthocyanin content of berry skin was analyzed by HPLC over the same period, showing that anthocyanins accumulated more rapidly in CT berries compared to controls (Figure 1B), and the difference in anthocyanin content between the two samples was already significantly higher by the second sampling date (JD 219) corresponding to full veraison, and gradually declined but was still evident at harvest (Figure 1B). HPLC analysis also revealed that the increase in total anthocyanins was not evenly distributed among the five main glucosylated species that characterize the Sangiovese cultivar (Figure 1C). We observed a significant increase in the levels of 3'4'-OH anthocyanins associated to cluster thinning, which modified the anthocyanin profile of CT berries, particularly increasing levels of the glucosylated form of peonidin compared to C berries.

Transcriptional modulation induced by cluster thinning

To investigate the molecular changes that take place in response to cluster thinning, we carried out a comparative microarray analysis of CT and C berries at time points JD 211, 227 and 266, corresponding to the beginning of veraison (BV), the end of veraison (EV) and harvest (H).

Principal Component Analysis (PCA) of the global transcriptomic data revealed enough uniformity among the three biological replicates to defined associations between treatments (Figure 2). The two principal components, explaining about the 50% of the overall variance, allowed us to clearly separate C and CT at the EV stage, whereas the separation was less clear-cut at the BV and H stages suggesting that the main transcriptomic changes induced by cluster thinning occurred at the EV stage.

thumbnailFigure 2. Principal component analysis (PCA) shows that the most severe changes are at the end of veraison. Biological replicates relative to time point beginning of veraison (BV), end of veraison (EV) and harvest (H) are represented as circles for CT and as triangles for C.

To identify the gene expression profiles with the greatest contribution to the differences between the C and CT transcriptomes, a multiclass comparison analysis was carried out using Significance Analysis of Microarray (SAM) with a false discovery rate (FDR) of 2% (TMev 4.3). We identified 1626 genes modulated during C berry development and 6033 modulated during CT berry development, with a fold change ≥ 2 in at least one comparison (Additional File 1). To evaluate the principal modifications triggered by cluster thinning, we focused on genes with a fold change ≥ 5, narrowing the analysis to 567 genes modulated during C berry development and 2466 genes modulated during CT berry development. A comparison of these datasets indicated three different sets of modulated transcripts. The first grouped 447 genes modulated in both treatments, the second grouped 2019 genes modulated only in CT berries and the third grouped 120 genes modulated only in C berries (Figure 3). For convenience, genes that showed less than a five-fold change in expression were described as 'not highly modulated'.

Additional file 1. Differentially expressed genes along berry development displaying a two-fold or greater change in transcript abundance between EV and BV or EV and H, in CT and C treatments. For each gene the annotation, the EV/BV and the H/BV Fold-change (FC) are indicated. Data obtained for CT and C are listed in two separate worksheets.

Format: XLS Size: 907KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

thumbnailFigure 3. Differentially expressed genes during ripening in CT and C berries. Transcripts are divided into three different gene datasets according to their high modulation in CT, in C or in both. The number of CT highly modulated genes either not expressed in control berries or expressed but not highly modulated is shown. Similarly the number of C highly modulated genes either not highly expressed in CT berries or expressed but not highly modulated is shown. The number of common modulated transcripts with equal or different modulation is specified.

Clustering analysis using Pearson's correlation distance divided the common, CT highly modulated and C highly modulated transcripts (Tables 2, 3 and 4) into eight groups representing the minimum number of profiles required to describe the three sampling time points. Clusters 1-4 represent genes that are downregulated during at least one analyzed time point compared to the BV stage, whereas clusters 5-8 represent genes that are upregulated during at least one analyzed time point compared to the BV stage.

Table 2. Cluster distribution of genes highly modulated both in C and CT berries.

Table 3. Distribution of 192 annotated transcripts (common to C and CT berries but differentially modulated during ripening) into eight clusters according to the gene expression trends.

Table 4. Cluster distribution of highly modulated genes in CT and C berries.

All the transcripts were annotated against the V1 version of the 12X draft annotation of the grapevine genome http://genomes.cribi.unipd.it/DATA/ webcite allowing 70% of the modulated genes to be identified (Additional File 2). To investigate the functional distribution of commonly and specifically modulated transcripts, we distributed them into 18 Gene Ontology (GO) functional categories and determined the percentage of genes in each category for each of the three data sets (Figure 4). The most represented functional categories, shared among the three datasets, were "Transport", "Transcription", "Secondary Metabolic Process", "Response to Stress" and "Cellular Process", which included the main genes involved in the physiology of berry ripening. The functional distribution of modulated genes was similar in the CT and C specific gene sets, suggesting that the CT treatment has a widespread effect on transcription rather than impacting on a specific functional category. The results for the common, CT highly modulated and C highly modulated genes are discussed in more detail below.

Additional file 2. Differentially expressed genes along berry development displaying a five-fold or greater change in transcript abundance between EV and BV or EV and H, in CT and C treatments. For each gene the annotation, the gene ontology, the EV/BV and the H/BV Fold-change (FC), the expression profile and the expression behavior in the other treatment are reported. Data obtained for CT and C are listed in two separate worksheets.

Format: XLS Size: 605KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

thumbnailFigure 4. Distribution of common, CT and C highly modulated transcripts into 18 GO functional categories. Transport: GO:0006810; Transcription: GO:0006350; Signal Transduction: GO:0007165; Secondary Metabolic Process: GO:0019748; Response to Stress: GO:0006950; Response to Hormone Stimulus: GO:0009725; Protein Metabolic Process: GO:0019538; Nucleic Acid Metabolic Process: GO:0090304; Nitrogen Compound Metabolic Process: GO:0006807; Lipid Metabolic Process: GO:0006629; Hormone Metabolic Process: GO:0042445; Generation of precursor metabolites and energy: GO:0006091; Developmental Process: GO:0032502; Cellular Process: GO:0009987; Cellular Homeostasis: GO:0019725; Cellular amino acid and derivative metabolic process: GO:0006519; Cell wall organization or biogenesis: GO:0071554; Carbohydrate metabolic process: GO: 0005975. Percentages were calculated on the total number of modulated genes for each dataset.

Common modulated genes

The first gene set contains 447 genes that are modulated in both C and CT berries, representing 79% of the genes modulated in C berries and 22% of those modulated in CT berries. Approximately equal numbers of genes for each treatment were upregulated in comparison to the BV stage and were therefore distributed into clusters 5-8 (132 for C berries and 130 for CT berries, Table 2). Similarly, there were approximately equal numbers of downregulated genes distributed into clusters 1-4 (315 for C berries and 317 for CT berries) although the nature of the distribution was distinct (Table 2). In particular, there were more genes downregulated throughout ripening (cluster 2) and specifically between BV to EV (cluster 4) in CT berries, whereas there were more genes downregulated specifically between EV and H (cluster 3) in C berries. Interestingly, very few common genes were detected in clusters 1 and 5, which represent EV-specific downregulation (1) and upregulation (5) in both treatments.

We identified 205 genes with similar expression profiles in both treatments and 242 with expression profiles that differed between C and CT berries during ripening (Figure 3, Additional File 2). The cluster distribution for the 192 genes with successful functional annotations belonging to the latter group is shown in Table 3.

Many genes belonging to "Generation of Precursor Metabolites and Energy", the most representative functional category, and to several others such as "Carbohydrate Metabolic Process", "Cell Wall Organization or Biogenesis", "Hormone Metabolic Process", "Response to Hormone Stimulus", "Response to Stress" and "Transport", were allocated to cluster 3 in C berries (downregulated specifically between EV and H) but shifted to cluster 2 in CT berries (downregulated throughout ripening). This pattern was particularly evident among the 24 genes with a role in photosynthesis, the six involved in the cell wall changes associated with berry softening (e.g. pectate lyase, pectinesterase and xyloglucan endotransglucosylase), the three involved in sugar metabolism (fructose-bisphosphate aldolase, the vacuolar invertase GIN1 and galactinol synthase) and the two involved in ethylene biosynthesis, i.e. 1-aminocyclopropane-1-carboxylic acid synthase (ACS) and 1-aminocyclopropane-1-carboxylic acid oxidase (ACO). These data are presented in Additional File 2. For the vacuolar invertase GIN1 (VIT_160022g00670), the expression profiles in C and CT berries were confirmed by real-time RT-PCR (Additional File 3).

Additional file 3. Real time RT-PCR validation of six selected genes. Expression profiles measured by real time RT-PCR are determined by calculating the relative expression ratio value for each stage relative to the BV stage. Real time RT-PCR data are reported as means ± SE of three biological replicates, obtained by using two reference genes. An actin beta/gamma 1 (VIT_12s0178g00200) and an elongation factor 1 (VIT_06s0004g03220) were used as control genes. For each gene the expression profile obtained by microarray analysis is shown on the right side. A-B: Vacuolar invertase 1, GIN 1 (VIT_16s0022g00670); C-D: PAL [Vitis vinifera] (VIT_16s0039g01120); E-F: Flavonol synthase (VIT_18s0001g03430); G-H: β-Fructosidase -invertase (VIT_00s2527g00010); I-J: Flavonoid 3’5’-hydroxylase (VIT_06s0009g02910); K-L: UDP-glucose:flavonoid glucosyltransferase (VIT_04s0023g01290).

Format: PDF Size: 103KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

We also noted several genes that were allocated to cluster 7 in C berries (upregulated specifically between EV to H) but to cluster 6 in CT berries (upregulated throughout ripening) (Table 3). These genes were predominantly involved in secondary metabolic processes, e.g. five PAL isogenes (VIT_00s2508g00010, VIT_16s0039g01120, VIT_16s0039g01320, VIT_16s0039g01300, VIT_00s2849g00010) and two stilbene synthases (STSs, VIT_16s0100g01070, VIT_16s0100g00850). We confirmed the expression profile of one PAL (VIT_16s0039g01120) in C and CT berries by real-time RT-PCR (Additional File 3).

In addition, several of the 205 genes whose expression profile did not change qualitatively after cluster thinning showed marked quantitative differences in expression (Additional File 4). In particular, we observed stronger downregulation of several genes involved in photosynthesis, cell wall metabolism, stress responses and hormone metabolism in CT berries. Notable examples included genes involved in the dynamic remodeling of cell wall polysaccharides, such as a xyloglucan endotransglucosylase/hydrolase (VIT_00s0386g00050), which was downregulated 10-fold more strongly in CT berries compared to controls, and genes related to hormone metabolism, such as a cytokinin dehydrogenase (VIT_18s0001g13200).

Additional file 4. Differentially expressed genes along berry development displaying a similar expression profile between C and CT. For each gene, the annotation, the gene ontology, the profile and the CT/C Fold-change (FC) ratio at EV and at H are reported.

Format: XLS Size: 49KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Whereas many common genes were downregulated more strongly in CT berries compared to controls, very few were induced more strongly after cluster thinning. However, examples included galactinol synthase (VIT_07s0005g0970), which is an important regulator of carbon partitioning, as well as a STS (VIT_09s0018g01490) and a PAL involved in the phenylpropanoid/flavonoid pathway (VIT_16s0039g01170). Remarkably, the glutathione-S-transferase VvGST4 (VIT_04s0079g00690), whose involvement in berry ripening has been reported previously [21], was upregulated throughout ripening in CT berries compared to controls. We validated the expression profile of one gene in this group by real-time RT-PCR, the flavonol synthase VIT_18s0001g03430 (Additional File 3).

Genes highly modulated in response to cluster thinning

The second group of modulated genes comprised 2019 CT highly modulated transcripts that could be divided into eight clusters according to their expression profiles (Table 4). Approximately equal numbers of genes were downregulated (885, clusters 1-4) and upregulated (1134, clusters 5-8). Unlike the cluster distribution of the common genes discussed above, many CT highly modulated genes were assigned to clusters 1 and 5, and most were assigned to clusters 1, 2, 4, 5, 6 and 8, characterized by prompt modulation in response to cluster thinning (i.e. changes already visible at EV). This suggests that physiological and biochemical changes in the berry almost occur immediately after the treatment. We investigated the expression of the CT highly modulated genes in C berries (Table 4), which allowed us to distinguish between those not expressed in C berries and those expressed but not highly modulated in C berries. Among the 885 genes grouped in clusters 1-4, 82 were not expressed in the controls and 803 were expressed but not highly modulated during ripening. The expression of 82 genes in CT berries but not in C berries at the earliest sampling stage is unlikely to reflect a response to treatment and is probably the consequence of inter-vine variability at the BV stage. Certain agronomic parameters such as titratable acidity and soluble sugar content showed higher standard deviations at BV than the subsequent stages (Figure 1A), confirming that samples collected at BV could be more heterogeneous than those collected later. This could likewise explain the presence of C highly downregulated transcripts that are not expressed in CT berries (Table 4). Surprisingly, a significantly majority of the 1134 genes in clusters 5-8 were not expressed in C berries, while only a minority was expressed but not highly modulated in C berries.

The 2019 CT highly modulated genes were assigned to 18 functional categories as shown in Figure 4. Given the major influence of carbohydrate metabolism, secondary metabolism and the transport of carbohydrate and secondary metabolites on berry quality traits, we will focus on genes in the categories "Carbohydrate Metabolic Process", "Secondary Metabolism" and "Transport", which represent 18% of all the genes modulated in CT berries following the treatment.

Genes belonging to the "Carbohydrate Metabolic Process" category

We identified 68 CT highly modulated genes involved in carbohydrate metabolism. These are listed in Table 5, which provides the Gene_ID, annotation, metabolic process, cluster assignment, expression in C berries, and the EV/BV and H/BV fold change. We found that 28 of these genes grouped in clusters 1-4 and 40 (including 33 not expressed in control berries) grouped in clusters 5-8. The modulated genes are involved in various primary metabolic pathways, including sucrose and starch metabolism, glycolysis, the pentose phosphate pathway and the Krebs cycle, indicating a large-scale reprograming of carbohydrate metabolism in response to cluster thinning. This includes the downregulation of two invertases (VIT_18s0072g01040 and VIT_00s2527g00010) and one sucrose-phosphate synthase (VIT_05s0029g01140) probably representing the decline in photosynthesis during berry ripening. The downregulation of VIT_00s2527g00010 was confirmed by real time RT-PCR (Additional File 3).

Table 5. CT highly modulated genes in the "Carbohydrate metabolic process" functional category.

We identified several genes involved in the synthesis and/or degradation of sugars and starch that were specifically induced by cluster thinning, including two genes involved in sucrose re-synthesis (the sucrose synthase VIT_00s1562g00010, and the sucrose-phosphate synthase VIT_18s0075g00330), four genes involved in glycolysis (the glyceraldehyde-3-phosphate dehydrogenases VIT_14s0171g00440 and VIT_18s0089g00590, and the pyruvate kinases VIT_07s0005g00440 and VIT_07s0005g00430), four involved in starch metabolism (the α-amylase VIT_03s0063g00410, the β-amylase, VIT_05s0051g00010, the isoamylase VIT_00s0131g00420 and the granule-bound starch synthase VIT_16s0022g00740), and four involved in the production of osmolytes such as trehalose (VIT_00s0233g00030, VIT_14s0036g01210, VIT_00s0173g00110 and VIT_11s0037g00710).

Microarray analysis also revealed that two CT highly modulated malate dehydrogenase (MDH) isogenes in two different cellular compartments showed opposite expression profiles during ripening, i.e. downregulation of the glyoxysomal isoform compared to control berries (VIT_03s0088g01190) but upregulation of the cytosolic isoform (VIT_07s0005g03360). Interestingly, another cytosolic MDH (VIT_07s0005g03350) was upregulated at EV and thereafter in CT berries but not expressed in controls, along with five NADP-dependent malic enzyme (ME) isogenes (VIT_03s0038g00040, VIT_16s0039g01050, VIT_16s0039g00580, VIT_15s0045g00190, VIT_04s0008g00180) and one phosphoenolpyruvate carboxykinase (VIT_07s0205g00070). Altogether these data suggest that the malate degradation is specifically influenced by cluster thinning, a hypothesis supported by the lower titratable acidity of CT berries compared to controls during ripening (Figure 1A and Table 1).

The oxidative burst that occurs in berries at the onset of ripening [22] seems to be enhanced by the thinning treatment, as supported by the upregulation of two alcohol dehydrogenases (ADHs, VIT_04s0044g01130 and VIT_04s0044g01120), three aldehyde dehydrogenases (ALDHs, VIT_01s0137g00090, VIT_01s0137g00080 and VIT_06s0004g02060) and two aldo/keto reductases (VIT_05s0062g00990 and VIT_08s0070g01010).

Genes belonging to the "Secondary Metabolism Process" category

We identified 108 CT highly modulated genes putatively involved in secondary metabolism (Table 6), 25 of which were downregulated (clusters 1-4), and 83 (including 70 not expressed in control berries) of which were upregulated (clusters 5-8). The downregulated genes included those related to the general phenylpropanoid pathway and phenolic acid metabolism. In particular, we detected two isoforms of 4-coumarate-CoA ligase (4CLs, VIT_17s0000g01790 and VIT_14s0171g00300), one caffeic acid 3-O-methyltransferase (COMT, VIT_08s0007g04520), one cinnamoyl-CoA reductase (VIT_16s0039g01670) and one cinnamyl alchohol dehydrogenase (VIT_19s0014g04980). A small number of genes involved in the flavonoid pathway were also downregulated in CT berries, including a flavonoid 3'-hydroxylase (F3'H, VIT_09s0002g01090), a F3H (VIT_16s0098g00860), a LDOX (VIT_08s0105g00380) and the leucoanthocyanidin reductase LAR1 (VIT_01s0011g02960). The downregulation of such genes may reflect the more robust suppression of proanthocyanidin biosynthesis in CT berries compared to controls, which may also account for the downregulation of a putative UDP-glucose:flavonoid glucosyltrasferase (VIT_16s0115g00340). This gene could be related to the reduced glucosylation of flavonoid compounds such as proanthocyanidin monomers produced during the herbaceous berry growth phase. Glucosylated proanthocyanidin does not accumulate in grapevine tissues but transient glucosylation might be necessary for the vacuolar import of monomers as reported in Medicago truncatula [23]. The downregulation of a geranylgeranyl diphosphatase synthase (VIT_15s0024g00850) and a geranylgeranyl reductase (VIT_17s0000g06280) in CT berries suggests that terpenoid metabolism is also affected by the treatment.

Table 6. CT highly modulated genes in the "Secondary metabolic process" functional category.

The upregulated genes included several transcripts related to the general phenylpropanoid pathway, e.g. four PAL isogenes (VIT_16s0039g01100, VIT_16s0039g01240, VIT_16s0039g01280 and VIT_16s0039g01130) and one 4CL (VIT_08s0007g05050). This suggests more phenolic precursors enter the multibranched phenylpropanoid pathway in CT berries. Indeed, many downstream pathways appeared to be specifically activated after cluster thinning. The induction of stilbene synthesis is suggested by the high upregulation of 19 STSs that are not expressed in control berries (many induced as early as the EV stage), the induction of isoflavone synthesis is indicated by the upregulation of four isoflavone reductase genes (VIT_03s0038g04620, VIT_03s0038g04680, VIT_03s0038g04690 and VIT_03s0038g04630), and the activation of phenolic acid metabolism is indicated by the upregulation of a COMT (VIT_02s0025g02940), a ferulate 5-hydroxylase (VIT_11s0016g01050) and five cinnamyl alcohol dehydrogenase genes (VIT_03s0110g0034, VIT_15s0107g00210, VIT_13s0067g00560, VIT_02s0025g03110 and VIT_13s0064g00290).

A small number of genes from the flavonoid pathway were also induced by cluster thinning, including the flavonoid 3',5'-hydroxylase genes F3'5'Hi (VIT_06s0009g02910), F3'5'Hk (VIT_06s0009g03010), the F3'Hb (VIT_17s0000g07210) [24] and the DFR (VIT_18s0001g12800) described by [25]. The expression profile of F3'5'Hi was further confirmed by real time RT-PCR (Additional File 3). Three putative UDP-glucose:flavonoid glucosyltrasferase transcripts (VIT_12s0134g00620, VIT_06s0004g06400 and VIT_16s0050g01590) were also upregulated, and these could also be involved in anthocyanin synthesis although this has yet to be confirmed.

In addition to the cinnamyl alcohol dehydrogenase genes, a few other genes responsible for the synthesis of flavor compounds were also induced by cluster thinning, including germacrene D-synthase (VIT_19s0014g02550), geraniol 10-hydroxylase (VIT_15s0048g01480), limonoid UDP-glucosyltrasferase (VIT_03s0091g00040) and a carotenoid cleavage dioxygenase (VIT_13s0064g00810), suggesting that the aromatic profile of berries is also affected by the treatment.

A very important consequence of cluster thinning was the high modulation of several members of the GST gene family (Table 6). Twenty different GST genes were activated in CT berries, including 16 tau-type (U), one phi-type (F), one lambda-type (L) and two zeta-type (Z) enzymes, according to the classification devised by Edwards et al. [26]. Five U-type and one F-type CT highly modulated GSTs were expressed but not highly modulated in control berries, whereas 13 U-type GSTs were not expressed in control berries at all. Only one of these transcripts (VIT_19s0015g02610) has already been functionally characterized in grapevine and it corresponds to VvGST5 (VIT_19s0015g02610 [21]). This was strongly upregulated between the EV and H stages in CT berries, which parallels the accumulation of anthocyanins. However, VvGST5 could not induce anthocyanin accumulation in transient assays carried out by Conn et al. [21], indicating that the enzyme is unlikely to have a direct role in anthocyanin biosynthesis. Only two GST transcripts (VIT_16s0039g00990 and VIT_05s0049g01120) were highly downregulated in CT berries.

Genes belonging to the "Transport" category

We identified 175 CT highly modulated genes in the "Transport" category, 36 of which are putatively involved in the transport of carbohydrates (10) or secondary metabolites (26) (Table 7). Of the ten putative carbohydrate transporters, seven were expressed but not highly modulated during ripening in control berries but were highly downregulated after cluster thinning. These included the sucrose transporter VvSUC11 (VIT_18s0001g08220 [27]), the polyol transporter VvPMT5 (VIT_03s0063g02250 [28]), an ERD6-like sugar transporter (VIT_07s0104g00830 [29]), a glucose-6-phosphate/phosphate-translocator required for import into plastids (VIT_11s0052g00430 [30]) and a succinate/fumarate mitochondrial transporter (VIT_08s0217g00010), all of which were downregulated at the EV stage. The remaining two transcripts, representing sucrose transporter VvSUC27 (VIT_18s0076g00250 [27]) and a 2-oxoglutarate/malate carrier protein (VIT_18s0001g07320) were downregulated at both the EV and H stages. In contrast, three dicarboxylate/tricarboxylate carriers (VIT_00s0188g00090, VIT_07s0031g02470 and VIT_08s0007g07270) were strongly upregulated by cluster thinning. Together with the strong induction of malate-degrading enzymes, the induction of dicarboxylate/tricarboxylate carriers reinforces the idea that thinning has a direct and specific impact on malate metabolism.

Table 7. CT highly modulated genes in "Transport" functional category specifically involved in carbohydrate and secondary metabolite translocation.

Among the 26 putative secondary metabolite transporters highly modulated in CT berries, 19 belonged to the ATP-Binding Cassette (ABC) transporter family and seven to the Multidrug and Toxic Compound Extrusion (MATE) transporter family. Five of the transcripts were downregulated in CT berries, while 21 were upregulated (19 of which were not expressed in control berries at all). The ABC and MATE transporters were analyzed by querying the protein sequence of the probe consensus against the NCBI databases using BLASTX (Additional File 5). This was necessary because ABC and MATE transporters play multiple roles at the cellular level. None of the sequences we identified matched the grapevine anthocyanin-acylglucoside MATE transporters that have been functionally characterized [31]. Using Arabidopsis thaliana as a reference organism, we found that most of the transporter sequences matched significantly (e-value < 10-50) with an Arabidopsis protein homolog (Additional File 5). Furthermore, two CT highly upregulated MATE transporters matched TRANSPARENT TESTA 12 (AtTT12), a MATE transporter associated with flavonoid sequestration in vacuoles [32]. These sequences were VIT_18s0001g11760 (e-value = 3 × 10-96), which was not highly modulated in control berries but was upregulated in CT berries in the EV and H stages, and VIT_11s0052g01540 (e-value = 5 × 10-74), which was not expressed in control berries but was upregulated at the EV stage in CT berries. Exploration of the latest release of the Vitis vinifera genome revealed that the two MATE transporters mentioned above are not the closest homologs of AtTT12 (data not shown). However, the potential role of these two TT12-like genes in flavonoid transport cannot be ruled out.

Additional file 5. Differentially expressed genes along CT berry development encoding ABC and MATE transporters. Analysis was performed by querying the protein sequence of the probe consensus against the NCBI databases using BLASTX. For each gene the correspondent Arabidopsis thaliana homologue, the NCBI reference sequence and E-value are specified.

Format: XLS Size: 16KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Genes highly modulated in control berries

The third group of 120 genes was highly modulated in control berries (Table 4). Although, like the common and CT highly modulated genes, these grouped into eight clusters, more than half of them were assigned to clusters 6 and 7, showing that most C highly modulated genes were upregulated during at least one ripening phase. There were 41 downregulated genes (clusters 1-4), six of which were not expressed at all in CT berries and 35 of which were expressed but not highly modulated during ripening after the cluster thinning treatment. There were 79 upregulated genes (clusters 5-8), 20 of which were not expressed in CT berries and 59 of which were expressed, but not highly modulated.

Using the annotation criteria described above, 15 C highly modulated genes were assigned to the "Transport" category, only one of which was involved in sugar transport (Table 8). This was the putative hexose transporter VvHT5 (VIT_05s0020g03140), which was induced strongly between the EV and H stages. Six genes involved in carbohydrate metabolism were highly modulated in control berries. We identified an ADH (VIT_06s0004g04320), an L-idonate dehydrogenase (VIT_16s0100g00290) which is involved in tartaric acid biosynthesis [33], a galactinol synthase (VIT_14s0060g00760) and a trehalose 6-phosphate phosphatase (VIT_18s0001g05300), all of which were downregulated, and an aldo-keto reductase (VIT_05s0077g01300) and a 2-phosphoglycerate kinase (VIT_06s0061g00280) that were upregulated. It is interesting to note that different isoforms of all these enzymes except L-idonate dehydrogenase were also found in the list of CT highly modulated genes. The ADH was highly downregulated in control berries but not highly modulated after cluster thinning, indicating that ADH activity is needed in CT berries as already suggested by the upregulation of two ADH genes in thinned vines.

Table 8. C highly modulated genes in the "Carbohydrate metabolic process", "Secondary metabolic process" and "Transport" (carbohydrate and secondary metabolite) functional categories.

Very few genes related to phenol metabolism were highly modulated in control berries, the exceptions being an isoflavone reductase (VIT_18s0001g12690), a cinnamate-4-hydroxylase (VIT_11s0065g00350) and one transcript putatively involved in anthocyanin modification (VIT_04s0023g01290) which were upregulated. The third of these transcripts was not expressed at all in CT berries as confirmed by real-time RT-PCR (Additional File 3). We did not identify any GST genes that were highly modulated in control berries, or any genes with a putative role in secondary metabolite transport or storage.

Cluster thinning mainly affects the berry transcriptome at the end of veraison

To complete the analysis of our microarray results, we compared the C and CT berry transcriptomes at each time point using a SAM unpaired comparison with a FDR of 2%. This revealed 4123 genes that were differentially expressed in C and CT berries at EV, and 178 at H, in each case with a fold change ≥ 2 (Additional File 6). No genes were found to be differentially expressed in C and CT berries at the BV stage. As anticipated by the previous PCA, these data indicate that C and CT berries at the BV stage are indistinguishable at the transcriptomic level, the main changes occurring at the EV stage, followed by minor further changes at the H stage. As described above, we then focused on genes with a fold change ≥ 5, narrowing the analysis to 1167 differentially-expressed genes at EV and 53 at H, 40 of which were common to both stages. We found that 833 genes were upregulated and 347 were downregulated in CT berries at EV and/or at H.

Additional file 6. Differentially expressed genes displaying a two-fold or greater change in transcript abundance between C and CT at EV and H time points. For each gene the description and the CT/C Fold-change (FC) are indicated. Data obtained for EV and H are listed in two separate worksheets.

Format: XLS Size: 467KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Approximately 80% (940) of these genes had already been identified as differentially expressed either in C or in CT berries by the cluster analysis based on SAM multiclass described above (Additional File 7). Because these genes were mined from the microarray data using two different approaches, we can have greater confidence that they represent a genuine molecular response to thinning and most (893) are indeed included in the list of genes highly modulated in response to CT. The remaining 20% of genes identified only by the direct comparison at each time point represent genes with a less than five-fold difference in expression between C and CT during ripening, and which were therefore excluded from the first cluster analysis.

Additional file 7. Differentially expressed genes displaying a five-fold or greater change in transcript abundance between C and CT at EV and /or H. For each gene the CT/C Fold-change (FC) at EV and /or H, the description, the gene ontology, the presence of the gene among genes identified by cluster analysis as CT specific, C specific, common to C and CT with a different expression profile and common to C and CT with the same expression profile, are reported..

Format: XLS Size: 200KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

As described for the first dataset, we determined the functional categories of the identified genes and focused on those related to carbohydrate metabolism, secondary metabolism and the transport of carbohydrates and secondary metabolites (Table 9). All genes belonging to these categories were differentially expressed in C and CT berries, but only at the EV stage. Only two genes were related to carbohydrate metabolism: a trehalose-phosphate synthase potentially involved in stress tolerance which was downregulated in CT berries (VIT_19s0014g00300), and a pyruvate decarboxylase involved in cellular respiration which was upregulated in CT berries (VIT_15s0024g00630). All the genes related to secondary metabolism were upregulated in CT berries, including two cinnamyl alcohol dehydrogenases (VIT_13s0047g00760 and VIT_00s0371g00010), two GSTs (VIT_07s0104g01810 and VIT_06s0004g05680) and two UDP-glucose:flavonoid 5,3-O-glucosyltransferases (VIT_18s0041g00840 and VIT_18s0041g00950). The induction of these genes reinforces our cluster analysis results, which already showed that the expression of several cinnamyl alcohol dehydrogenase, GST and putative UDP-glucose:flavonoid glucosyl transferase genes was triggered by thinning. The same holds true for the four secondary metabolite transporters, three of which were upregulated (VIT_15s0045g01030, VIT_14s0030g00900, VIT_10s0003g03470) and one of which was downregulated in CT berries (VIT_16s0100g00350).

Table 9. Differentially expressed genes in C and CT berries at the EV stage belonging to "Carbohydrate metabolic process", "Secondary metabolic process" and "Transport" functional categories.

Discussion

Although many previous studies have considered the impact of increasing the source/sink ratio on grape berry composition, ours is the first investigation to look at the consequences of cluster thinning on global gene expression profiles, which is the basis of most of the observed physiological and biochemical changes.

Although there is only limited evidence for a strict relationship between yield and quality [34], in this study the effectiveness of thinning reflected the suboptimal leaf area to yield ratio of control vines (0.6 m2/kg). Increasing this ratio to 1.2 m2/kg by cluster thinning boosted the levels of sugars and anthocyanins and reduced acidity at harvest. Several authors have reported similar effects on sugars and anthocyanins [35,36], but there is disagreement on the impact of thinning on titratable acidity, with reports suggesting that acidity is unaffected [7,37], slightly increased [1] or decreased by yield reduction [2,38]. In all cases, the effect on acidity seems to be related to the impact of cluster thinning on ripening, particularly the soluble solids content, because thinning reduces acidity only when the soluble solids concentration is strongly and positively affected ([38,39] and this study).

The grouping of gene expression profiles by principal component analysis showed that C and CT berries could be distinguished as early as two weeks after treatment. Microarray data were then analyzed by two different approaches. Genes that were differentially expressed in C and CT berries were initially clustered on the basis of their expression profiles, and then we carried out a direct comparison of C and CT transcriptomes at each time point.

The first approach revealed such a large number of genes modulated during ripening that we applied a > 5 fold change cut-off threshold before assigning a gene to the modulated group. This resulted in 2466 transcripts that were considered to be modulated in CT berries and 567 that were modulated in C berries, including 447 that were common to both treatments. Approximately half of the common genes were more strongly modulated in CT berries, including several downregulated genes involved in photosynthesis, carbon utilization, carbohydrate metabolism, cell wall modifications and hormone metabolic processes that are already known to have a role in berry ripening [22,40], and several upregulated genes involved in the normal ripening process (e.g. genes related to secondary metabolite biosynthesis). These data strongly suggest that the entire course of berry ripening is enhanced by the cluster thinning.

An interesting and unexpected result of the microarray analysis was the relatively large number of genes (2019) highly modulated only in CT berries. More than a half of these genes were never expressed in control berries at any ripening stage, and were activated uniquely by the cluster thinning treatment. This shows that thinning is able to trigger the transcription of genes that otherwise would not be activated in untreated berries and therefore that the effect of thinning goes beyond the simple enhancement or acceleration of the normal ripening process. This appears to affect many different metabolic and cellular processes because the CT highly modulated genes are distributed throughout all 18 functional categories we considered.

The CT high modulation of genes involved in carbohydrate metabolism supports the impact of cluster thinning on sugar accumulation commencing at the EV stage. This reflects the achievement of an optimal balance between leaf area and yield in thinned vines (1.2 m2/kg) compared to controls (0.6 m2/kg) as previously reported [2-7]. Interestingly, all CT highly modulated genes involved in sugar transport (including VvSUC11 and VvSUC27) were downregulated in CT berries but expressed at a constant level throughout ripening in controls. The downregulation of sucrose transporters in CT berries contrasts with the increase in sucrose transporter mRNA during berry development reported by Davies et al. [27] and the general enhancing effect of thinning on the entire ripening process. One possible explanation is that the higher sugar concentration triggers negative feedback that affects the sucrose transporters. Indeed the presence of sugar-response elements in the promoters of various sucrose transporter genes, potentially acting as cis-regulatory elements involved in sugar signaling, has recently been reported [41].

Despite the significant decline in starch concentration following veraison, several genes involved in starch biosynthesis and modification are modulated during ripening [40]. We detected both highly up- and down-regulated genes in CT involved in starch degradation, and their role is unclear given that developing CT berries accumulate large amounts of sugar. Simultaneous starch synthesis and degradation may facilitate the unloading and storage of sugars in the ripening fruit [39]. Starch-degrading enzymes might also provide carbon backbones for the biosynthesis of secondary metabolites, which could also act as signaling molecules in the regulation of genes controlling phenylpropanoid synthesis [17,42].

The strong induction of genes encoding malate-degrading enzymes in CT berries, together with the induction of dicarboxylate-tricarboxylate carriers (which transport malate across the mitochondrial membrane thus supplying substrates for the Krebs cycle) supports the specific modulation of malate metabolism by cluster thinning. Malate, whose catabolism is considered responsible of total acidity reduction in grape berry after véraison [43] is liberated from the vacuole in post-véraison becoming available for catabolism through various avenues [44,45] and, with the advance of ripening, malic acid is likely a vital source of carbon.

Although advanced malate degradation is heavily dependent upon the extent to which berry temperature is elevated [46,47] in response to increased sunlight exposure, this condition was not tested in our experiment as the cluster microclimate was not modified by the removal of one cluster of each shoot, maintaining an unchanged canopy structure.

In this respect, our microarray data suggest that malic acid catabolism is accelerated by cluster thinning following the general enhanced ripening process as seen in other researches [48].

Genes involved in the metabolism of phenylpropanoids and aromatic compounds were strongly modulated in CT berries. Several transcripts involved in phenolic acid, stilbene, flavonoid and isoflavonoid metabolism were affected by the treatment, suggesting that these distinct branches of the phenylpropanoid pathway are directly affected by the increased source/sink balance. However it is possible that the synthesis of phenolic compounds such as stilbenes or isoflavonoids could be part of a systemic response to wounding resulting from the removal of berry clusters, since these compounds are normally produced by the plant in response to stress conditions such as wounding or interactions with pathogens [49,50].

The anthocyanin content of CT berries was higher than that of controls (Figure 1B and Table 1) but the only CT highly upregulated transcript related to the flavonoid/anthocyanin pathway that could explain this result was DFR. However, looking specifically for known anthocyanin-related transcripts, we found that VvGST4 (VIT_04s0079g00690) and VvMYBA1 (VIT_02s0033g00410) were also more strongly upregulated in CT berries compared to controls (Additional File 2). The high induction of three flavonoid glucosyltransferases and the large number of CT highly upregulated GSTs and transporters of the ABC and MATE families may also play an important role in triggering anthocyanin biosynthesis, although their precise functions in berry ripening remain to be determined. Several putative flavonoid-related transcripts such as F3H, LDOX, LAR1 and UDP-glucose:flavonoid glucosyltransferase were highly downregulated in CT berries, which may reflect the slowing down of the synthesis of non-anthocyanin flavonoid compounds such as proanthocyanidin. Several reports indicate that agronomic treatments such as water stress and cluster thinning can increase total anthocyanins and also produce a shift in their profile [3,4,51,52]. The anthocyanin composition of CT berries differed from controls in our study, with higher relative levels of peonidin 3-glucoside and total 3'4'-OH anthocyanins, in agreement with earlier reports [3,51]. This might reflect the specific upregulation of F3'Hb in CT berries during the EV and H phases. Although F3'5'Hi and F3'5'Hk were upregulated in CT berries and not expressed at all in controls, there was a significant decrease in the levels of 3'4'5'-OH anthocyanins at harvest, probably because this expression pattern was restricted just to the EV stage. This suggests that the biosynthesis of 3'4'-OH and 3'4'5'-OH anthocyanins is controlled independently [5].

Cluster thinning is known to change the aromatic profile of grape berries and increase the resveratrol content in wine [53,54]. In agreement, our microarray data reveal the induction of many putative aroma-modifying transcripts (e.g. geraniol 10-hydroxylase, germacrene D-synthase, as well as two isoforms each of ADH and ALDH) and the specific upregulation of 19 STSs.

The number of genes highly modulated in control berries (567) was small in comparison with previous studies [22,40] but this reflects our stringent application of the requirement for a fold change > 5. Only 120 of the control highly modulated genes were not highly modulated in CT berries, none of which were involved in carbohydrate or secondary product metabolism. Therefore we assume they play a relatively minor role in the control of berry quality traits at harvest.

The direct comparison of C and CT transcriptomes at each time point indicated that the greatest number of genes differentially expressed in C and CT berries were found at the EV stage, whereas only minor differences were detectable at harvest. This generally supported the cluster analysis results because almost all the differentially expressed genes had already been detected. The novel data obtained by direct comparison related mainly to genes whose expression was only slightly modulated during ripening.

Conclusions

We investigated the effect of cluster thinning on the Sangiovese berry transcriptome in association with agronomic parameters and the biochemical properties of berries during ripening. The increased source-sink ratio achieved by cluster thinning reflected the optimal leaf area/crop weight ratio in thinned vines (1.2 m2/kg) compared to 0.6 m2/kg in controls. This caused the sugar and anthocyanin content to increase from veraison to harvest at an accelerated rate, along with extensive transcriptomic reprogramming involving both the increased modulation of genes that are normally regulated during ripening and the induction of a large group of genes that are not expressed in ripening control berries. Cluster thinning therefore accelerates the normal ripening process but also superimposes a metabolic environment involving the induction of novel processes not found in ripening control berries and the repression of some pathways that are part of the normal ripening process. Although the possibility that the transcriptomic reprogramming was partially an effect and not only the cause of the observed metabolic changes cannot be excluded, our data provide a significant contribution to current understanding of the molecular consequences of cluster thinning, specifically the impact of increasing the source/sink ratio at veraison on carbohydrate and anthocyanin accumulation.

Methods

Plant material

The trial was carried out in 2008 on adult Sangiovese (Vitis vinifera L.) vines, clone 12T grafted to SO4, in a vineyard in Bologna, Italy (44°30'N, 11°24'E), with North-South oriented rows (vine spacing 1 m × 2.8 m, vertical shoot positioned spur pruned cordon training system (12 buds per vine), with cordon height at 1.0 m above ground and a canopy height of 1.3-1.4 m). Pest management was carried out according to Regione Emilia Romagna local practice. Shoots were manually trimmed when most started to outgrow the top wire, which occurred on Julian Day (JD) 192. Nine vines per treatment, with the same cluster number at flowering (16 per vine), were selected in a single uniform row and each vine was randomly assigned in three blocks to the two following groups: control (C, no treatment) and cluster thinning (CT, removal of 50% of total clusters per vine at veraison, JD 211).

Berry sampling

Samples of 40 berries, taken from three vines in each block, were collected at the following stages: a) beginning of veraison (BV, berries softening and °Brix ~8; JD 211); b) full veraison (one week later; JD 219); c) end of veraison (EV, soft and fully-colored berries; JD 227); d) ripening (two weeks after stage c, JD 250); e) full ripening (H, harvest, JD 266). The samples were divided in two and 20 berries were processed immediately in order to monitor berry ripening, while the rest were frozen at -20°C in preparation for HPLC analysis of anthocyanins. We also collect 30 additional samples at stages a, c and e (same method as above), which were immediately frozen in liquid nitrogen and stored at -80°C for subsequent microarray analysis.

Agronomic parameters and biochemical analysis

Final vine leaf area was estimated after the achievement of a linear relationship between shoot length (cm) and leaf area (cm2) for at least 15 shoots for both C and CT samples collected from vines at the end of shoot growth. Leaf area was measured with a leaf area meter (LI-3000A, Li-Cor Biosciences, Lincoln, Nebraska, USA), ensuring the contribution of the main and lateral leaf areas were kept separate. Two separate correlations for the main and lateral areas (y) and shoot lengths (x) were so established (C main: y = 16.212 x, R2 = 0.70; C laterals: y = 14.834x, R2 = 0.70; CT main: y = 13.441x, R2 = 0.71; CT laterals: y = 15.844x, R2 = 0.89). After leaf fall, each shoot per vine was then measured and the shoot lengths were used to estimate the corresponding pooled total leaf area per vine using the linear relationship calculated above.

The must obtained after crushing the first 20 berries per sample was immediately filtered through a strainer and a drop was used for °Brix analysis with a temperature-compensating CR50 refractometer (Maselli Misure Spa, PR, Italy). We then diluted 5 ml of the same must in seven volumes of double distilled water and titrated this against 1 N, 0.5 N or 0.25 N NaOH (Sigma-Aldrich, St. Louis, MO, USA) with a Crison Compact Tritator (Crison, Barcelona, Spain) according to the stage of berry ripening, to obtain pH and titratable acidity data (expressed as g/L of tartaric acid equivalents).

The second 20 berries were used to extract anthocyanins for HPLC analysis as described by Mattivi et al. [55], using a Waters 1525 instrument (Waters, Milford, MA) equipped with a diode array detector (DAD) and using a Phenomenex (Castel Maggiore, BO, Italy) reversed-phase column (RP18, 250 mm × 4 mm, 5 μM). Anthocyanins were quantified at 520 nm using an external calibration curve with malvidin-3-glucoside chloride as the standard (Sigma-Aldrich, St. Louis, MO).

We recorded cluster number and weight, berry weight and total yield per vine for each tested vine at harvest. Statistical analysis of agronomic parameters and biochemical data were carried out using the mixed General Linear Model (GLM) procedure of the SAS statistical package (SAS Institute, Cary, NC), and treatment comparisons were carried out using the Tukey test.

Microarray analysis

Total RNA for microarray analysis was isolated from ~200 mg of the ground berry tissue without seeds using the Spectrum™ Plant Total RNA kit (Sigma-Aldrich, St. Louis, MO). RNA quality and quantity were determined using a Nanodrop 2000 instrument (Thermo Scientific, Wilmington, DE) and a Bioanalyzer Chip RNA 7500 series II (Agilent, Santa Clara, CA). The cDNA synthesis, labeling, hybridization and washing reactions were performed according to the NimbleGen Arrays User's Guide (V 3.2). Each hybridization was carried out on a NimbleGen microarray 090818 Vitis exp HX12 (Roche, NimbleGen Inc., Madison, WI), representing 29,549 predicted genes on the basis of the 12X grapevine V1 gene prediction version http://srs.ebi.ac.uk/ webcite. The chip probe design is available at the following URL: http://ddlab.sci.univr.it/FunctionalGenomics/. webcite

The microarray was scanned using a ScanArray 4000XL (Perkin-Elmer) at 532 nm (Cy-3 absorption peak) and GenePix Pro7 software (Molecular Devices, Sunnyvale, CA, USA) according to the manufacturers' instructions. Images were analyzed using NimbleScan v2.5 software (Roche), which produces Pair Files containing the raw signal intensity data for each probe and Calls Files with normalized expression data derived from the average of the intensities of the four probes for each gene. The normalized gene expression data were finally converted to log2 values. All microarray expression data are available at GEO under the series entry GPL13936.

Pearson Correlation analysis and Principal Component Analysis (PCA; SIMCA P+ Umetrics, Umea, Sweden) were carried out to evaluate the robustness of the three biological replicates in each sample. A gene was considered to be expressed if the normalized expression value was higher than the value obtained by averaging the fluorescence of the negative control present on the chip, for at least two of the three biological replicates. A Significance Analysis of Microarrays (SAM) was implemented using TMeV software http://www.tm4.org/mev webcite, with a false discovery rate of 2%. Cluster analysis was performed by the k-means method with Pearson's correlation distance (TMeV) comparing EV and H gene expression to BV.

Real-time PCR

First strand cDNA was synthesized using 1 μg of total RNA as the template and the Improm-II™ Reverse Transcription system (Promega), according to the manufacturer's instructions. The total RNA was derived from the three biological replicates used for microarray analysis and all RNA samples were first treated with DNase I (Promega). Gene-specific primers were designed for six genes using the sequence information in the 3'-UTR, using actin and elongation factor 1 (EF1) genes as references [56] (Additional File 8). Primers and cDNA were mixed with the Power SYBR® Green PCR Master Mix (Applied Biosystems, Foster City, CA, USA) and the reaction was carried out on an ABI PRISM StepOne Sequence Detection System (Applied Biosystems, Foster City, CA, USA) using the following cycling conditions: 50°C hold for 2 min and a 95°C hold for 10 min followed by 45 cycles at 95°C for 30 s, 55°C for 30 s and 72°C for 20 s. Nonspecific PCR products were identified by the dissociation curves. Amplification efficiency was calculated from raw data using LingRegPCR software [57]. The relative expression ratio was calculated for development time points relative to the first sampling time point (beginning of veraison, BV) according to the Pfaffl equation [58]. Standard error (SE) values were calculated according to Pfaffl et al. [59]. The final data was calculated as previously reported [60].

Additional file 8. List of the primers used for real time RT-PCR validation experiment.

Format: XLS Size: 18KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Authors' contributions

CP set up the RNA extraction, performed microarray and HPLC experiments, participated in data analysis, prepared figures and tables and wrote the initial manuscript draft. SZ set up the real-time RT-PCR protocol, performed data analysis and participated in the drafting of the manuscript. GBT participated in the interpretation of data and supervised and critically revised the manuscript. GA and GV sampled the material and performed the HPLC analysis. SDS set up the microarray experiment and participated in the real time RT-PCR experiments and gene annotation. CI conceived the study and participated in the project design. MP and IF conceived the study, participated in its design and coordination and contributed to and critically revised the manuscript. It is the authors' opinion that CP and SZ contributed equally as first authors to this manuscript. All authors have read and approved the final version of the manuscript.

Acknowledgements

This work was supported by: Fondazione Cariverona (Completamento e attività del Centro di genomica Funzionale Vegetale), Verona, Italy. Special thanks to Alberto Ferrarini and Marianna Fasoli for their support in the hybridization analysis.

References

  1. Kliewer WM, Dokoozlian NK: Leaf Area/Crop Weight Ratios of Grapevines: Influence on Fruit Composition and Wine Quality.

    Am J Enol Vitic 2005, 56(2):170-181. OpenURL

  2. Kliewer WM, Weaver RJ: Effect of Crop Level and Leaf Area on Growth, Composition, and Coloration of 'Tokay' Grapes.

    Am J Enol Vitic 1971, 22(3):172-177. OpenURL

  3. Filippetti I, Ramazzotti S, Centinari M, Brucchetti B, Intrieri C: Esperienze triennali sugli effetti del diradamento dei grappoli sulla qualità delle uve della cultivar Sangiovese.

    Italus Hortus 2007, 14:412-416. OpenURL

  4. Guidoni S, Allara P, Schubert A: Effect of Cluster Thinning on Berry Skin Anthocyanin Composition of Vitis vinifera cv. Nebbiolo.

    Am J Enol Vitic 2002, 53(3):224-226. OpenURL

  5. Guidoni S, Ferrandino A, Novello V: Effects of Seasonal and Agronomical Practices on Skin Anthocyanin Profile of Nebbiolo Grapes.

    Am J Enol Vitic 2008, 59(1):22-29. OpenURL

  6. Intrieri C, Filippetti I, Valentini G, Seghetti L: Il diradamento migliora la qualità del Montepulciano a pergola.

    L'informatore Agrario 2005, 44:59-64. PubMed Abstract | Publisher Full Text OpenURL

  7. Keller M, Mills LJ, Wample RL, Spayd SE: Cluster Thinning Effects on Three Deficit-Irrigated Vitis vinifera Cultivars.

    Am J Enol Vitic 2005, 56(2):91-103. OpenURL

  8. Dokoozlian NK, Hirschfelt DJ: The Influence of Cluster Thinning at Various Stages of Fruit Development on Flame Seedless Table Grapes.

    Am J Enol Vitic 1995, 46(4):429-436. OpenURL

  9. Guidoni S, Argamante N: Influenza del diradamento dei grappoli sull'accumolo di antociani nelle uve.

    Quaderno di Viticoltura ed Enologia Università di Torino 2003, 26:27-42. OpenURL

  10. Palliotti A, Cartechini A: Cluster thinning effetcts on yield and grape composition in different grapevine cultivars.

    Acta Hort 2000, 512:111-120. OpenURL

  11. Downey MO, Dokoozlian NK, Krstic MP: Cultural Practice and Environmental Impacts on the Flavonoid Composition of Grapes and Wine: A Review of Recent Research.

    Am J Enol Vitic 2006, 57(3):257-268. OpenURL

  12. Pirie A, Mullins MG: Changes in anthocyanin and phenolics content of grapevine leaf and fruit tissues treated with sucrose, nitrate, and abscisic Acid.

    Plant Physiol 1976, 58(4):468-472. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Roubelakis-Angelakis KA, Kliewer WM: Effects of Exogenous Factors on Phenylalanine Ammonia-Lyase Activity and Accumulation of Anthocyanins and Total Phenolics in Grape Berries.

    Am J Enol Vitic 1986, 37(4):275-280. OpenURL

  14. Nagira Y, Ozeki Y: A system in which anthocyanin synthesis is induced in regenerated torenia shoots.

    J Plant Res 2004, 117(5):377-383. PubMed Abstract | Publisher Full Text OpenURL

  15. Nagira Y, Shimamura K, Hirai S, Shimanuki M, Kodama H, Ozeki Y: Identification and characterization of genes induced for anthocyanin synthesis and chlorophyll degradation in regenerated torenia shoots using suppression subtractive hybridization, cDNA microarrays, and RNAi techniques.

    J Plant Res 2006, 119(3):217-230. PubMed Abstract | Publisher Full Text OpenURL

  16. Solfanelli C, Poggi A, Loreti E, Alpi A, Perata P: Sucrose-specific induction of the anthocyanin biosynthetic pathway in Arabidopsis.

    Plant Physiol 2006, 140(2):637-646. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Vitrac X, Larronde F, Krisa S, Decendit A, Deffieux G, Merillon JM: Sugar sensing and Ca2+-calmodulin requirement in Vitis vinifera cells producing anthocyanins.

    Phytochemistry 2000, 53(6):659-665. PubMed Abstract | Publisher Full Text OpenURL

  18. Gollop R, Even S, Colova-Tsolova V, Perl A: Expression of the grape dihydroflavonol reductase gene and analysis of its promoter region.

    J Exp Bot 2002, 53(373):1397-1409. PubMed Abstract | Publisher Full Text OpenURL

  19. Gollop R, Farhi S, Perl A: Regulation of the leucoanthocyanidin dioxygenase gene expression in Vitis vinifera.

    Plant Science 2001, 161:579-588. Publisher Full Text OpenURL

  20. Zheng Y, Tian L, Liu H, Pan Q, Zhan J, Huang W: Sugars induce anthocyanin accumulation and flavanone 3-hydroxylase expression in grape berries.

    Plant Growth Regulation 2009, 58(3):251-260. Publisher Full Text OpenURL

  21. Conn S, Curtin C, Bezier A, Franco C, Zhang W: Purification, molecular cloning, and characterization of glutathione S-transferases (GSTs) from pigmented Vitis vinifera L. cell suspension cultures as putative anthocyanin transport proteins.

    J Exp Bot 2008, 59(13):3621-3634. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Pilati S, Perazzolli M, Malossini A, Cestaro A, Dematte L, Fontana P, Dal Ri A, Viola R, Velasco R, Moser C: 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 veraison.

    BMC Genomics 2007, 8:428. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  23. Zhao J, Dixon RA: MATE transporters facilitate vacuolar uptake of epicatechin 3'-O-glucoside for proanthocyanidin biosynthesis in Medicago truncatula and Arabidopsis.

    Plant Cell 2009, 21(8):2323-2340. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Falginella L, Castellarin SD, Testolin R, Gambetta GA, Morgante M, Di Gaspero G: Expansion and subfunctionalisation of flavonoid 3',5'-hydroxylases in the grapevine lineage.

    BMC Genomics 2010, 11:562. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  25. Sparvoli F, Martin C, Scienza A, Gavazzi G, Tonelli C: Cloning and molecular analysis of structural genes involved in flavonoid and stilbene biosynthesis in grape (Vitis vinifera L.).

    Plant Mol Biol 1994, 24(5):743-755. PubMed Abstract | Publisher Full Text OpenURL

  26. Edwards R, Dixon DP, Walbot V: Plant glutathione S-transferases: enzymes with multiple functions in sickness and in health.

    Trends Plant Sci 2000, 5(5):193-198. PubMed Abstract | Publisher Full Text OpenURL

  27. Davies C, Wolf T, Robinson SP: Three putative sucrose transporters are differentially expressed in grapevine tissues.

    Plant Science 1999, 147(2):93-100. Publisher Full Text OpenURL

  28. Klepek YS, Geiger D, Stadler R, Klebl F, Landouar-Arsivaud L, Lemoine R, Hedrich R, Sauer N: Arabidopsis POLYOL TRANSPORTER5, a new member of the monosaccharide transporter-like superfamily, mediates H+-Symport of numerous substrates, including myo-inositol, glycerol, and ribose.

    Plant Cell 2005, 17(1):204-218. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Buttner M: The monosaccharide transporter(-like) gene family in Arabidopsis.

    FEBS Lett 2007, 581(12):2318-2324. PubMed Abstract | Publisher Full Text OpenURL

  30. Niewiadomski P, Knappe S, Geimer S, Fischer K, Schulz B, Unte US, Rosso MG, Ache P, Flugge UI, Schneider A: The Arabidopsis plastidic glucose 6-phosphate/phosphate translocator GPT1 is essential for pollen maturation and embryo sac development.

    Plant Cell 2005, 17(3):760-775. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  31. Gomez C, Terrier N, Torregrosa L, Vialet S, Fournier-Level A, Verries C, Souquet JM, Mazauric JP, Klein M, Cheynier V, et al.: Grapevine MATE-type proteins act as vacuolar H+-dependent acylated anthocyanin transporters.

    Plant Physiol 2009, 150(1):402-415. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Debeaujon I, Peeters AJ, Leon-Kloosterziel KM, Koornneef M: The TRANSPARENT TESTA12 gene of Arabidopsis encodes a multidrug secondary transporter-like protein required for flavonoid sequestration in vacuoles of the seed coat endothelium.

    Plant Cell 2001, 13(4):853-871. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  33. DeBolt S, Cook DR, Ford CM: L-tartaric acid synthesis from vitamin C in higher plants.

    Proc Natl Acad Sci USA 2006, 103(14):5608-5613. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Jackson DI, Lombard PB: Environmental and Management Practices Affecting Grape Composition and Wine Quality - A Review.

    American Journal of Enology and Viticulture 1993, 44(4):409-430. OpenURL

  35. Mazza G, Fukumoto L, Delaquis P, Girard B, Ewert B: Anthocyanins, phenolics, and color of Cabernet Franc, Merlot, and Pinot Noir wines from British Columbia.

    J Agric Food Chem 1999, 47(10):4009-4017. PubMed Abstract | Publisher Full Text OpenURL

  36. Reynolds AG, Price SF, Wardle DA, Watson BT: Fruit Environment and Crop Level Effects on Pinot noir. I. Vine Performance and Fruit Composition in British Columbia.

    American Journal of Enology and Viticulture 1994, 45(4):452-459. OpenURL

  37. Ough CS, Nagaoka R: Effect of Cluster Thinning and Vineyard Yields on Grape and Wine Composition and Wine Quality of Cabernet Sauvignon.

    American Journal of Enology and Viticulture 1984, 35(1):30-34. OpenURL

  38. Di Profio F, Reynolds AG, Kasimos A: Canopy Management and Enzyme Impacts on Merlot, Cabernet franc, and Cabernet Sauvignon. I. Yield and Berry Composition.

    American Journal of Enology and Viticulture 2011, 62(2):139-151. Publisher Full Text OpenURL

  39. Agasse A, Vignault G, Kappel C, Conde C, Gerò H, Delrot S: Sugar transport & sugar sensing in grape. In Grapevine Molecular Physiology Biotechnology. 2nd edition. Edited by Springer. New York; 2009:105-139. OpenURL

  40. 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. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  41. Afoufa-Bastien D, Medici A, Jeauffre J, Coutos-Thevenot P, Lemoine R, Atanassova R, Laloi M: The Vitis vinifera sugar transporter gene family: phylogenetic overview and macroarray expression profiling.

    BMC Plant Biology 2009, 10(1):245. OpenURL

  42. Larronde F, Krisa S, Decendit A, Chèze C, Deffieux G, Mèrrilon JM: Regulation of polyphenols production in Vitis vinifera cell suspension cultures by sugars.

    Plant Cell Reports 1998, 17:946-950. Publisher Full Text OpenURL

  43. Ruffner HP, Possner D, Brem S, Rast DM: The physiological role of malic enzyme in grape ripening.

    Planta 1984, 160(5):444-448. Publisher Full Text OpenURL

  44. Famiani F, Walker RP, Tecsi L, Chen ZH, Proietti P, Leegood RC: An immunohistochemical study of the compartmentation of metabolism during the development of grape (Vitis vinifera L.) berries.

    J Exp Bot 2000, 51(345):675-683. PubMed Abstract | Publisher Full Text OpenURL

  45. Sweetman C, Deluc LG, Cramer GR, Ford CM, Soole KL: Regulation of malate metabolism in grape berry and other developing fruits.

    Phytochemistry 2009, 70(11-12):1329-1344. PubMed Abstract | Publisher Full Text OpenURL

  46. Lakso AN, Kliewer WM: The Influence of Temperature on Malic Acid Metabolism in Grape Berries. II. Temperature Responses of Net Dark CO2 Fixation and Malic Acid Pools.

    Am J Enol Vitic 1978, 29(3):145-149. OpenURL

  47. Ruffner HP, Hawker JS, Hale CR: Temperature and enzymic control of malate metabolism in berries of Vitis vinifera.

    Phytochemistry 1976, 15(12):1877-1880. Publisher Full Text OpenURL

  48. Dokoozlian NK, Kliewer WM: Influence of Light on Grape Berry Growth and Composition Varies during Fruit Development.

    Journal of the American Society for Horticultural Science 1996, 121(5):869-874. OpenURL

  49. Gutierrez-Gonzalez JJ, Guttikonda SK, Tran LS, Aldrich DL, Zhong R, Yu O, Nguyen HT, Sleper DA: Differential expression of isoflavone biosynthetic genes in soybean during water deficits.

    Plant Cell Physiol 2010, 51(6):936-948. PubMed Abstract | Publisher Full Text OpenURL

  50. Xu W, Yu Y, Zhou Q, Ding J, Dai L, Xie X, Xu Y, Zhang C, Wang Y: Expression pattern, genomic structure, and promoter analysis of the gene encoding stilbene synthase from Chinese wild Vitis pseudoreticulata.

    J Exp Bot 2011, 62(8):2745-2761. PubMed Abstract | Publisher Full Text OpenURL

  51. Bucelli P, Giannetti F: Incidenza del diradamento dei grappoli sulla composizione dell'uva e sulla qualità del vino.

    Rivista di Viticoltura ed Enologia 1996, 2:59-67. OpenURL

  52. Castellarin SD, Pfeiffer A, Sivilotti P, Degan M, Peterlunger E, G DIG: Transcriptional regulation of anthocyanin biosynthesis in ripening fruits of grapevine under seasonal water deficit.

    Plant Cell Environ 2007, 30(11):1381-1399. PubMed Abstract | Publisher Full Text OpenURL

  53. Prajitna A, Dami IE, Steiner TE, Ferree DC, Scheerens JC, Schwartz SJ: Influence of Cluster Thinning on Phenolic Composition, Resveratrol, and Antioxidant Capacity in Chambourcin Wine.

    Am J Enol Vitic 2007, 58(3):346-350. OpenURL

  54. Reynolds AG, Schlosser J, Sorokowsky D, Roberts R, Willwerth J, de Savigny C: Magnitude of Viticultural and Enological Effects. II. Relative Impacts of Cluster Thinning and Yeast Strain on Composition and Sensory Attributes of Chardonnay Musque.

    Am J Enol Vitic 2007, 58(1):25-41. OpenURL

  55. Mattivi F, Guzzon R, Vrhovsek U, Stefanini M, Velasco R: Metabolite profiling of grape: Flavonols and anthocyanins.

    J Agric Food Chem 2006, 54(20):7692-7702. PubMed Abstract | Publisher Full Text OpenURL

  56. 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. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  57. Ramakers C, Ruijter JM, Deprez RHL, Moorman AFM: Assumption-free analysis of quantitative real-time polymerase chain reaction (PCR) data.

    Neuroscience Letters 2003, 339(1):62-66. PubMed Abstract | Publisher Full Text OpenURL

  58. Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR.

    Nucleic Acids Res 2001, 29(9):e45. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  59. Pfaffl MW, Horgan GW, Dempfle L: Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real-time PCR.

    Nucleic Acids Res 2002, 30(9):e36. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  60. Mestdagh P, Van Vlierberghe P, De Weer A, Muth D, Westermann F, Speleman F, Vandesompele J: A novel and universal method for microRNA RT-qPCR data normalization.

    Genome Biol 2009, 10(6):R64. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL