Cabernet Sauvignon, Sauvignon Blanc and Cabernet Franc are characterized by specific volatile thiols and methoxypyrazines. Enzymes involved in the production of these aromas have been recently characterized. Phenolic compounds play a central role in the physical mouthfeel properties of red wine; recent work relates quality with tannin levels. While the grape genotype has a tremendous impact on tannin content, the environment also plays a very large role in grape composition. The pathway for phenolic biosynthesis is well known, but the mechanisms of environmental influence are poorly understood. Ultimately, there is an interaction between molecular genetics and the environment. Flavor is influenced by climate, topography and viticultural practices. For example, water deficit alters gene expression of enzymes involved in aroma biosynthesis in grapes, which is genotype dependent, plastic gutter and may lead to increased levels of compounds, such as terpenes and hexyl acetate, that contribute to fruity volatile aromas. The grapevine berry can be subdivided into the skin, pulp and seeds.
A thick waxy cuticle covers the epidermis. The hypodermal cells contain chloroplasts, which lose their chlorophyll at veraison and become modified plastids; they are the sites of terpenoid biosynthesis and carotenoid catabolism. Anthocyanins and tannins accumulate in the vacuoles of hypodermal cells . Pulp cells are the main contributors to the sugar and organic acid content of the berries. Pulp cells also have a much higher set of transcripts involved in carbohydrate metabolism, but a lower set of transcripts involved in lipid, amino acid, vitamin, nitrogen and sulfur metabolism than in the skins. Hormones can influence berry development and ripening. Concentrations of auxin, cytokinins and gibberellins tend to increase in early fruit development of the first stage. At veraison, these hormone concentrations have declined concomitant with a peak in abscisic acid concentration just before veraison. Auxin prolongs the Stage 2 lag phase and inhibits anthocyanin biosynthesis and color development in Stage 3. Grapevine, a nonclimacteric fruit, is not very sensitive to ethylene; however, ethylene appears to be necessary for normal fruit ripening.
Ethylene concentration is highest at anthesis, but declines to low levels upon fruit set; ethylene concentrations rise slightly thereafter and peak just before veraison then decline to low levels by maturity. Ethylene also plays a role in the ripening of another non-climacteric fruit, strawberry. ABA also appears to be important in grape berry ripening during veraison when ABA concentrations increase resulting in increased expression of anthocyanin biosynthetic genes and anthocyanin accumulation in the skin. ABA induces ABF2, a transcription factor that affects berry ripening by stimulating berry softening and phenylpropanoid accumulation. In addition, ABA affects sugar accumulation in ripening berries by stimulating acid invertase activity and the induction of sugar transporters. It is not clear whether ABA directly affects flavor volatiles , but there could be indirect effects due to competition for common precursors in the carotenoid pathway. Many grape berry ripening studies have focused on targeted sampling over a broad range of berry development stages, but generally with an emphasis around veraison, when berry ripening is considered to begin. In thisstudy, a narrower focus is taken on the late ripening stages where many berry flavors are known to develop in the skin. We show that that the abundance of transcripts involved in ethylene signaling is increased along with those associated with terpenoid and fatty acid metabolism, particularly in the skin.
Cabernet Sauvignon clusters were harvested in 2008 from a commercial vineyard in Paso Robles, California at various times after veraison with a focus on targeting °Brix levels near maturity. Dates and metabolic details that establish the developmental state of the berries at each harvest are presented in Additional file 1. Berries advanced by harvest date with the typical developmental changes for Cabernet Sauvignon: decreases in titratable acidity and 2- isobutyl-3-methoxypyrazine concentrations and increases in sugar and color . Transcriptomic analysis focused on four harvest dates having average cluster °Brix levels of 22.6, 23.2, 25.0 and 36.7. Wines made in an earlier study from grapes harvested at comparable levels of sugars or total soluble solids to those in the present study showed clear sensory differences. Six biological replicates, comprising two clusters each, were separated into skins and pulp in preparation for RNA extraction and transcriptomic analysis using the NimbleGen Grape Whole-Genome Microarray. Thus, a 4 × 2 factorial experimental design was established. After standard microarray processing and data normalization, two-way ANOVA indicated that the transcript abundance of 16,280 transcripts statistically significantly changed across the °Brix levels below the adjusted p-value of 0.05 , the transcript abundance of 10,581 transcripts changed significantly across Tissue types, and the abundance of 2053 transcripts changed significantly with respect to the °Brix x Tissue interaction term p-value column: adjBrix, adjTissue or adjTissue*Brix. A note of caution must be added here. There are high similarities amongst members in certain Vitis gene families , making it very likely that cross-hybridization can occur with probes on the microarray with high similarity to other genes. We estimate approximately 13,000 genes have the potential for cross-hybridization, with at least one probe of a set of four unique probes for that gene on the microarray potentially cross-hybridizing with probes for another gene on the microarray. Genes with the potential for cross hybridization have been identified and are highlighted in light red in Additional file 2. The rationale to include them is that although individual genes can not be uniquely separated, the probe sets can identify a gene and its highly similar gene family members, thus, providing some useful information about the biological responses of the plant. An additional approach was taken, removing cross-hybridizing probes before quantitative data analysis . Many of the significant genes were unaffected by this processing, but 3600 genes were completely removed from the analysis. Thus, it was felt that valuable information was lost using such a stringent approach. The less stringent approach allowing for analysis of genes with potential cross hybridization was used here in the rest of the analyses. To assess the main processes affected by these treatments, the gene ontologies of significantly affected transcripts were analyzed for statistical significance using BinGO. Based on transcripts that had significant changes in abundance with °Brix level, 230 biological processes were significantly over represented in this group . The three top over represented processes were response to abiotic stress, biosynthetic process, and response to chemical stimulus, a rather generic set of categories. Tissue differences were more revealing at the stage when flavors peak; 4865 transcripts that were significantly higher in skins compared to pulp at 23.2 °Brix were tested for over represented GO functional categories . Some of the top GO categories included photosynthesis, isoprenoid biosynthesis, and pigment biosynthesis . Some of the transcripts with the largest differences between skin and pulp at 23.2 °Brix are β-ketoacyl-CoA synthase , taxane 10-β-hydroxylase , wax synthase, a lipase, an ABC transporter, blueberry container and phenylalanine ammonia-lyase . The abundance of 5716 transcripts was significantly higher in pulp than skin at 23.2 °Brix . Some of the top GO categories over represented were a variety of transport processes and small GTPase mediated signal transduction . Some of the transcripts with the largest differences in abundance with pulp greater than skin at 23.2 °Brix were polygalacturonase , flavonol synthase, stachyose synthase, an amino acid transporter, a potassium channel , and HRE2 . The transcript abundance of 2053 genes had significantly differential expression across °Brix levels and tissues .
The top GOcategories over represented in this set involved photosynthesis and phenylpropanoid metabolism, both associated with the berry skin . Other flavorcentric categories of the 57 categories over represented include aromatic compound biosynthesis, fatty acid metabolism and alcohol catabolism. This transcript set was further analyzed by dividing into 10 clusters using k-means clustering . The over represented GO categories were determined for each cluster . Eight of the 10 clusters had distinct over represented GO categories; two clusters did not have any over represented GO categories, meaning that the genes in these two clusters were assigned to GO categories of expected proportions when compared to the entire NimbleGen array. Clusters 1, 8, 9 and 10 had a large number of over represented categories. Many GO categories within a cluster are subsets of others in that cluster and were grouped together. For example, cluster 4 had four over represented GO categories, oxygen transport, gas transport, heat acclimation and response to heat. The four categories could be grouped into two, as two are subsets of the others; this is how they were listed in Table 1.It would be impossible to discuss here all the transcript abundance changes detected in these berries. As we were interested in compounds associated with berry flavors as they develop or change in the late stages of berry ripening, we took a more targeted approach for analysis with this in mind. Berries at 24° Brix are known to be near-optimal for flavor, thus we took a simple approach to look for genes that were peaking around this stage. We found some significant and large increases in transcript abundance between the 22.6 and 23.2 °Brix levels. A group of VviERF6 transcription factor paralogs represented 6of the top 10 transcripts increasing in transcript abundance from 22.6 to 23.2 °Brix in the skin, but not in the pulp . These VviERF6 TFs were also found in Cluster 8 . This is very interesting since many flavor compounds are derived from the skin and ERF TFs are known to be responsive to ethylene, a known fruit-ripening hormone. These VviERF TFs were named ERF105 in the annotation by Grimplet et al. , however they are more orthologous with AtERF6 as determined by a more comprehensive phylogenetic method using many plant species at Gramene . Annotation details of the V1 gene models of the VviAP2/ERF superfamily can be found in Additional file 8 including updated Vvi symbols according to its closest Arabidopsis ortholog as instructed by the Grapevine Gene Nomenclature System developed by the International Grape Genome Program Supernomenclature committee. This renaming of the AP2/ERF superfamily should facilitate comparative analyses and functions with other species, particularly Arabidopsis. To properly annotate the AP2/ERF superfamily of Vitis vinifera according to the IGGP Supernomenclature committee instructions, a phylogenetic tree was generated for the AP2/ERF superfamily of Arabidopsis thaliana and Vitis vinifera using the TAIR 10 and V1 gene models, respectively . The labeled family classifications were derived from the Arabidopsis naming scheme by Nakano et al.. There are 130 members in the VitisAP2/ERF superfamily in the Pinot Noir reference genome. However, the six paralogs of ERF6 discussed above belong to a Vitis vinifera clade in subfamily IX and are distinctly different or separate from any Arabidopsis subfamily IX ERF TFs . Using k-means clustering, VviERF6L1 fell within Cluster 8 with 369 transcripts, including five additional VviERF6 paralogs. The top GO categories associated with Cluster 8 were genes associated with terpenoid metabolism and pigment biosynthesis . Other interesting flavor associated categories included fatty acid and alcohol metabolism . Representative transcripts from Cluster 8 that were correlated with the transcript abundance profile of VviERF6L1 can be seen in Figure 4. These are ACC oxidase, which is involved in ethylene biosynthesis; a lipoxygenase, part of a fatty acid degradation pathway giving rise to flavor alcohols such as hexenol; α-expansin 1, a cell wall loosening enzyme involved in fruit softening, and two terpene synthases, which produce important terpenes that contribute to Cabernet Sauvignon flavor and aroma. The high similarity of these transcript profiles indicates that ethylene biosynthesis and signaling may be involved in the production of grape aroma. Supporting this argument, two recent studies have shown that a tomato ERF TF , falling in the same ERF IX subfamily, has a strong effect on ethylene signaling and fruit ripening. The transcript abundance of AtERF6 in Arabidopsis is strongly increased by ethylene, which is triggered by the MKK9/MPK3/MPK6 pathway. The transcript abundance of VviMKK9 in the Cabernet Sauvignon berries was higher in the skin than the pulp, but there were no significant differences for VviMPK3 or VviMPK6 . This is not too surprising since AtMKK9 activates AtMPK3 and AtMPK6 by phosphorylation. In addition, the transcript abundance of AtERF6 in Arabidopsis increases with ROS, SA, cold, pathogens, and water deficit. There were no visible signs of pathogen infection in these berries. Additional circumstantial evidence for ethylene signaling in the late stages of berry ripening was that the transcript abundance of many VviERF TFs was significantly affected by berry ripening and/or tissue .