- Research article
- Open Access
Metabolite and transcript profiling of berry skin during fruit development elucidates differential regulation between Cabernet Sauvignon and Shiraz cultivars at branching points in the polyphenol pathway
© Degu et al.; licensee BioMed Central Ltd. 2014
- Received: 21 March 2014
- Accepted: 11 July 2014
- Published: 26 July 2014
Grapevine berries undergo complex biochemical changes during fruit maturation, many of which are dependent upon the variety and its environment. In order to elucidate the varietal dependent developmental regulation of primary and specialized metabolism, berry skins of Cabernet Sauvignon and Shiraz were subjected to gas chromatography–mass spectrometry (GC-MS) and liquid chromatography–mass spectrometry (LC-MS) based metabolite profiling from pre-veraison to harvest. The generated dataset was augmented with transcript profiling using RNAseq.
The analysis of the metabolite data revealed similar developmental patterns of change in primary metabolites between the two cultivars. Nevertheless, towards maturity the extent of change in the major organic acid and sugars (i.e. sucrose, trehalose, malate) and precursors of aromatic and phenolic compounds such as quinate and shikimate was greater in Shiraz compared to Cabernet Sauvignon. In contrast, distinct directional projections on the PCA plot of the two cultivars samples towards maturation when using the specialized metabolite profiles were apparent, suggesting a cultivar-dependent regulation of the specialized metabolism. Generally, Shiraz displayed greater upregulation of the entire polyphenol pathway and specifically higher accumulation of piceid and coumaroyl anthocyanin forms than Cabernet Sauvignon from veraison onwards. Transcript profiling revealed coordinated increased transcript abundance for genes encoding enzymes of committing steps in the phenylpropanoid pathway. The anthocyanin metabolite profile showed F3′5′H-mediated delphinidin-type anthocyanin enrichment in both varieties towards maturation, consistent with the transcript data, indicating that the F3′5′H-governed branching step dominates the anthocyanin profile at late berry development. Correlation analysis confirmed the tightly coordinated metabolic changes during development, and suggested a source-sink relation between the central and specialized metabolism, stronger in Shiraz than Cabernet Sauvignon. RNAseq analysis also revealed that the two cultivars exhibited distinct pattern of changes in genes related to abscisic acid (ABA) biosynthesis enzymes.
Compared with CS, Shiraz showed higher number of significant correlations between metabolites, which together with the relatively higher expression of flavonoid genes supports the evidence of increased accumulation of coumaroyl anthocyanins in that cultivar. Enhanced stress related metabolism, e.g. trehalose, stilbene and ABA in Shiraz berry-skin are consistent with its relatively higher susceptibility to environmental cues.
- Metabolite profiling
- Grape berry metabolism
- Transcript analysis
Grape cultivars evolved from many origins  and were distributed to geographical regions with different climates  leading to the evolution of varietal specific physiology, morphology, growth and development ,. Grapevine vegetative growth and canopy structure affects the amount of light reaching the grape bunch and modifies its microclimate -. Moreover, the partitioning of photoassimilates and other metabolites to the bunch will also depend on the extent of vegetative growth during berry development and eventually it will determine the metabolite profile of the flesh and skin of the berry ,,. Hence, a complex interaction between genotype, macro- and micro-climate and plant vigor will define berry metabolism and eventually wine quality .
Grape berry, a non-climacteric fruit , follows a double sigmoidal growth pattern which is composed of three major phases , each characterized by changes in size, color, texture, metabolic composition and gene expression ,. The hard green stage, phase I, is characterized by rapid cell growth and elongation with a concomitant massive accumulation of organic acids, hydroxyl cinnamic acid, amino acids and tannins towards the end of this phase -. In phase two, berry enlargement ceases, sugars start to accumulate at the expense of acids eventually leading to a rise in pH . During veraison, berries enter a second rapid growth phase with a progressive decline in acidity and increase in softness. Veraison is a stage that marks the beginning of phase III, it is associated with berry softening and fundamental metabolite changes such as the accumulation of sugars, aroma precursors and pigments, , as well as other classes of specialized metabolites that will contribute to berry quality . The rate of biochemical changes during berry development is varietal dependent - as well as influenced by climatic factors -, soil , irrigation - and viticultural practices -.
The role of secondary metabolism in berries is not understood completely. Plants with fleshy and tasty fruits were likely favored during evolution by animal driven dispersal  and the pericarp and skin possibly evolved as protective tissues for developing seed . As such, polyphenols, including numerous flavonoids as well as non-flavonoid compounds, play a central protective role in plant adaptation to the environment . Flavonols protect berries from UV radiation , anthocyanins play a role in seed dispersal, provide protection from high temperature and UV radiation, and proanthocyanidins posses antifungal properties -. Flavonoids are biologically active compounds due to their antioxidant activities , and contribute to organoleptic properties of the fruit and wine ,.
Since the draft sequence of the Pinot Noir grapevine genome , the regulation of berry development and the associated metabolic processes were the focus of an increasing number of whole-genome gene expression studies ,,-. Nevertheless, there are extensive gaps in our understanding of the regulation of secondary metabolism and the differences in secondary metabolism amongst different cultivars. For example, recent transcriptomic and genomic studies on different commercial clones discovered hundreds of genes missing from the Pinot Noir reference genome -. Thus, more comparative studies are needed to fully elucidate the regulation of grape metabolism during berry development. Recently developed RNA-seq transcriptome profiling approaches, provides higher resolution and capability of detecting different isoforms of a transcript compared to the microarray based methods . To date, there are only a few studies in grape berry development that were performed using RNAseq - and none in comparative studies. Moreover, given the significant post-transcriptional regulation of fruit metabolism -, transcriptional studies are increasingly being supported by- and integrated with metabolic and proteomic data. Large-scale metabolite profiling has been applied to grape only recently ,,- and mostly using targeted analysis of primary or secondary metabolism.
In the present study we used GC-MS and LC-MS based metabolic profiling of the berry skin to study the developmental processes characteristic of two red wine grape varieties Shiraz (SH) and Cabernet Sauvignon (CS) grown in the field in a semiarid environment in southern Israel. The gas chromatography–mass spectrometry (GC-MS) and liquid chromatography–mass spectrometry (LC-MS) generated dataset was augmented with RNAseq-transcript profiles. The study aimed to examine varietal differences in the metabolism of developing berries and the underlying regulatory mechanism.
To the best of our knowledge, this is the first attempt to integrate metabolomics and RNAseq-based transcriptome data in a comparative study of berry development on two physiologically contrasting cultivars in grape.
Plant material and experimental conditions
The field experiment was conducted on four-year old commercial vines of two red wine grape varieties of Vitis vinifera L. (Shiraz and Cabernet Sauvignon) grafted on Ruggeri during the 2011 growing season. Both varieties were planted with 3 m spacing between rows and 1.5 m between vines. The vines were trained onto vertically oriented trellis in a East–west orientation. The vineyard is located at 30.68 N 34.80E; with an average rainfall of 86 mm and an average potential evaporation of over 2000 mm. Soil water content was measured on an hourly basis and it is given in Additional file 1: Figure S1. Irrigation was applied every four days based on evaporation from evaporation pan Class A (ETpan) using an irrigation factor of 0.25. Plants were drip irrigated using in-line emitters at a flow rate of 2.2 L h−1. To avoid any edge effects, only plants located in the middle of each row were used for the experiment. The experimental plot was bordered with vine row receiving a similar irrigation amount.
Berry quality related traits
Prior to each sampling, total soluble solids (TSS) in oBRIX were assayed from the berry juice of five different vines using a hand refractometer (Palette PR-100, Atago, USA). Berry weight was averaged from at least 15 berries selected from the middle part of representative bunches from eight plants of each cultivar on both sides (N and S) of the row. Titratable acidity measurement (g/L tartaric acid equivalents) was performed on the same samples according to the standard procedures used in Guymon and Ough .
Sampling and metabolite extraction
Samples were taken for metabolic analysis at four developmental stages in both cultivars (number of days are given in reference to veraison stage): 30 days before veraison (DTV) (pre-veraison); veraison; 20 days after veraison (DAV) (post-veraison) and harvest 37 DAV for Shiraz and 53 DAV for CS. Shiraz and CS were differentially sampled at harvest in relation to the oBRIX measured. At all sampling dates, five berries representative of the middle portion of the bunch and at similar developmental stage were collected on both N and S side of a given vine-plant and pooled together. Six independent biological replicates, each consisting of the pooled berries from six separate vines, were collected in the middle section of the row in both cultivars. Skin was carefully peeled from the berries, gently removing the remaining flesh from the skin with a clean cloth. Immediately after, the tissue was snap-frozen in liquid nitrogen and kept at −80°C until further analysis. Ahead of extraction samples were freeze dried in a lyophilizer (VIRTIS GARDINER, N.Y. R525, Model 10-MR-TR). The freeze-dried samples were extracted for parallel metabolite profiling (LC and GC/MS) using a slightly modified version of the method described previously . All chemicals were purchased from Sigma-Aldrich if not indicated otherwise. Berry skin tissue was ground using a RETCH-mill (Retsch Gmbh, 42787 Haan, Germany) with pre-chilled steel holders and grinding beads. For metabolite extraction, 70 mg of frozen powder were weighed and extracted in a pre-chilled methanol: chloroform: water extraction solution (2.5:1:1 v/v). Internal standards, (i.e. 0.2 mg/ml ribitol in water, 1 mg/ml ampicillin in water and 1 mg/ml corticosterone in methanol), were subsequently added. The mixture was then briefly vortexed, centrifuged for 2 min at 14000 RPM (microcentrifuge 5417R) and the supernatant was decanted into the new tubes. The supernatant was mixed with 300 μl of chloroform (LC/MS grade) and 300 μl of UPLC-grade water and then centrifuged at 14,000 RPM for 2 min. After that, 100 μl of the water/methanol phase was dried in a vacuum concentrator (Eppendorf Concentrator Plus) for derivatization  for GC-MS analysis. The remaining water/methanol phase was transferred to UPLC vials for LC-MS analysis.
GC-MS derivatization and data processing
GC-MS samples from the above extraction were re-dissolved and derivatized. Eight microliters of a retention time standard mixture (0.029% v/v n-dodecane, n-pentadecane, n-nonadecane, n-docosane, n-octacosane, n-dotracontane, and n-hexatriacontane dissolved in pyridine) was added. The sample set also included a reference quality control of authentic metabolite standards (1 mg ml−1 each) (Additional file 2: Table S1A). Volumes of 1 μL were then injected onto 30-m VF-5 ms GC column with 0.25 mm i.d., film thickness of 0.25 μ m, and + 10 m EZ-Guard (Agilent) in splitless and split mode (32:1) allowing a more accurate comparison of highly abundant metabolites (e.g. tartarate, sugars, and inositol). The GC-MS system consisted of an AS 3000 autosampler, a TRACE GC ULTRA gas chromatograph, and a DSQII quadrupole mass spectrometer (Thermo-Fisher ltd). The parameters of the machine were exactly as described in . Spectral searching was done by consulting the National Institute of Standards and Technology (NIST, Gaithersburg, USA) algorithm incorporated in the Xcalibur® data software (version 2.0.7) against RI libraries from the Max-Planck Institute for Plant Physiology in Golm, Germany (http://www.mpimp-golm.mpg.de/mms-library/) and finally normalized by the total metabolites and correcte d for the dilution factor.
For LC-MS analysis, 4 μl of extracted sample was injected onto a UPLC-QTOF-MS system equipped with an ESI interface (Waters Q-TOF XevoTM: Waters MS Technologies, Manchester, UK) operating in negative and positive ion modes. Chromatographic separation was carried out on an Acquity UPLC BEH C18 column (100 mm × 2.1 mm, 1.7 μm). The column and autosampler were maintained at 40°C and 10°C, respectively. During each sample running, the mobile phase comprised 95% water, 5% acetonitrile, 0.1% formic acid (phase A), and 0.1% formic acid in acetonitrile (phase B). The solvent gradient program was conditioned exactly as described previously . All analyses were acquired using leucine enkephalin for lock mass calibration to ensure accuracy and reproducibility, at a concentration of 0.4 ng L−1, in 50/50 of acetonitrile/H2O with 0.1% v/v formic acid. The MS conditions were set essentially as described previously .
UPLC data processing
MassLynxTM software (Waters) version 4.1 was used as the system controlling the UPLC and for data acquisition as described previously . The raw data acquired were processed using MarkerLynx application manager (Waters) essentially as described previously . To verify metabolite identification, representative samples of different developmental stages from each cultivar were run using the same instruments and under the same operating conditions at the metabolomics facility of the Edmund Mach Foundation in San Michele all’Adige – Italy, where an in-house standard library described in details in  was used to validate the annotation of the identified metabolites based on retention time order of commercial standards (Additional file 2: Table S1B), and metabolites were also identified based on a fragmentation pattern searched against the Chemspider metabolite database (http://www.chemspider.com/) and further confirmed with previous metabolite annotations ,-.
The normalized data set (to tissue dry weight and internal standards) was subjected to Kruskal-Wallis and Wilcoxon rank sum tests to identify metabolites changing significantly during berry development. Comparisons between cultivars at each sampling date were performed using the Student’s t-test. Pairwise correlations to all annotated metabolites were calculated using the Pearson correlation algorithm. The corresponding p-values were computed using the cor.test function of R as previously described . To avoid false positives, false discovery rate (FDR) at a Q value of 0.05 was performed. The correlation matrices were used as a data for network construction and analysis using “igraph” package (http://igraph.org/r/#docs). All the statistical analyses and network properties were calculated using the R-software environment R 3.0.1 (http://cran.r-project.org/).
To facilitate parallel comparisons to the metabolite data, frozen, ground skin tissues from the same samples used for metabolite extraction were used. Total RNA was extracted from three biological replicates of 70 mg of berry skin tissue at veraison and 20 DAV in the two cultivars essentially as described by Japelaghi et al. . The quality and concentration of extracted RNA was determined using Bioanalyzer Chip RNA 7500 series II (Agilent, Santa Clara, CA) and a Nanodrop 2000 spectrophotometer (Thermo Scientific, Wilmington, DE). Following quality assessment, poly (A) mRNA preparation and sequencing with an Illumina HiSeq 1000 sequencer (Illumina Inc., San Diego, CA, USA) were performed as described previously . The resultant reads were aligned to the reference Vitis vinifera genome using TopHat software (version 2.0.6), which received as input the Vitis vinifera GTF file (option “-G”) and also using the following parameters: “–b2-very-sensitive –r 150 –mate-std-dev 50” . Subsequently, the Cufflinks software ; version 2.0.2) was used to assemble aligned RNA-seq reads into transcripts with the parameters “--min-intron-length 10”; these assemblies were processed into a full transcriptome set by CuffMerge . Finally, their abundance was estimated with CuffDiff . In order to facilitate easy comparison and visualization of transcriptomic changes between developmental stages and the two cultivars, the data were normalized to the CS veraison stage transcript level.
Functional enrichment analysis of differentially expressed transcripts of CS and Shiraz was carried out based on VitisNet molecular networks ). The complete set of networks and their corresponding genes was downloaded from http://www.sdstate.edu/ps/research/vitis/upload/AdditionalFile2_2.xlsx. Enrichment against these networks was performed using a hyper-geometric test in Expander , Enrichment p-values were adjusted for multiple testing using the Bonferroni correction.
Sampling stages and berry quality traits
GC/MS based metabolic changes across berry development
Most amino acids showed progressive reduction towards maturation in both cultivars. However, pyruvate derived amino acids, alanine and valine, had a constant abundance during berry development for both cultivars (Figure 3). Coupled to a pronounced decrease in glycerate, amino acids, serine and glycine, that are derived from glycerate progressively decreased in relative abundance in both cultivars towards harvest but the reduction was higher in Shiraz compared to CS. The 2-oxoglutarate derived amino acids glutamate, arginine, ornithine and glutamine, and the oxaloacetate derived amino acids aspartate and threonine had a progressive decrease in relative abundance in both cultivars. Among shikimate-derived amino acids, phenylalanine showed maximum accumulation at veraison followed by a progressive reduction in relative abundance. The extent of reduction at harvest compared to its level at veraison was somewhat higher in Shiraz compared to CS. The other shikimate-derived amino acids, tyrosine and tryptophan, remained at relatively steady levels following veraison. In contrast to most amino acids, proline relative abundance increased towards maturity in both cultivars (Figure 3).
TCA cycle intermediates, citrate, succinate, fumarate and malate showed similar patterns of change in relative abundance declining progressively towards maturity. However, the extent of decrease was variety dependent. For example, the major organic acid, malate decreased by 6-fold in Shiraz compared to 4-fold reduction in CS. Reduction in citrate level was also higher in Shiraz (~9-fold) compared to CS (7-fold). The lesser abundant organic acids, succinate and maleate , decreased toward maturity by 18-fold and 24-fold in CS compared to 4-fold and 6-fold respectively, in Shiraz (Figure 3).
Quinate and shikimate links the primary metabolism to the phenylpropanoid biosynthesis pathways. These metabolites exhibited high reduction across berry development. Quinate and shikimate decreased on average by 14- and 20-fold respectively, at harvest in both cultivars. Cinnamate, the immediate precursor of the phenylpropanoid pathway, showed higher accumulation at pre-veraison followed by a progressive reduction towards maturity in both cultivars.
LC-MS based metabolic changes across berry development
Stilbenes are a metabolite class that branches from the phenylpropanoid pathway. Stilbenes accumulate in the berry in response to biotic - and abiotic stresses -. The two cultivars have similar patterns of accumulation of the stilbene, resveratrol, which exhibited a progressive accumulation in the berry skin towards maturity. In contrast, the stilbene, viniferin, did not change in abundance during development in both cultivars, while the pattern of change in the glucoside form of resveratrol, piceid, differed between cultivars; it gradually accumulated in Shiraz while keeping at comparable levels in CS following veraison (Figure 5).
Flavonols, the immediate competitors of the precursors for the anthocyanin pathway, showed mild induction following veraison stage with the exception of quercetin glucoside, which accumulated at pre-veraison and post veraison stages with a transient reduction in relative abundance at veraison. Unlike the other flavonols, kaempferol glucuronide and rutin decreased in relative abundance progressively towards maturity (Figure 5). All the flavanol classes (catechin, epicatechin and procyanidin dimers) decreased progressively in abundance towards maturity in both cultivars but the extent of reduction was more intense in Shiraz compared to CS (Figure 5) likely due to competition with the massive anthocyanin accumulation.
Anthocyanins, as expected, displayed progressive accumulation across development in both cultivars. Dihydrokaempferol represents the branching point between different anthocyanin downstream pathways (Figure 5, ). With flavonoid 3′-hydroxylase (F3′H), dihydrokaempferol is converted to dihydroquercetin and supports the biosynthesis of glucoside, acetyl and coumaroyl forms of cyanidin and peonidin (cyanidin-type anthocyanins), while with flavonoid 3′,5′ hydroxylase (F3′5′H) it can be converted to dihydromyricetin, that gives rise to the glucoside, acetyl, coumaroyl and caffeoyl forms of delphinidin, petunidin and malvidin (delphinidin-type anthocyanins). Cyanidin and peonidin glucoside biosynthesis had a mild but significant induction in ripening stages after veraison. In both cultivars the results indicated (Figure 5) that all anthocyanins derived from the dihydromyricetin pathway were highly induced compared with the parallel pathway derived from dihydroquercetin (Figure 5). The two cultivars exhibited similar pattern of anthocyanin accumulation with the exception of malvidin 3-O-glucoside, which progressively accumulated in CS, and petunidin, which significantly accumulated towards maturity in Shiraz (Additional file 2: Table S3). The two initial glucoside forms, cyanidin glucoside and delphinidin glucoside reached maximum level at mid-ripening and reduced at harvest (Figure 5). Comparing the coumaroyl forms of anthocyanins, Shiraz exhibited remarkably higher levels from veraison onwards (Additional file 1: Figure S3).
Metabolite correlation and network analysis
To gain an understanding of how metabolites coordinately changed during the season and to identify cultivar-specific metabolite relations, a pairwise correlation analysis was performed among all primary and specialized metabolites as described previously in Hochberg et al. . Metabolites in the same biosynthetic pathway mostly displayed a positive correlation (Additional file 2: Table S4A & S4B). The highest correlations were observed between primary and specialized metabolite classes. Anthocyanins and flavonols exhibited mostly similar trends of correlation, reflected in positive r values between the two classes and in negative r values with organic acids and most of other identified primary metabolites. Particularly, anthocyanins displayed strong negative relationship with flavanols, flavanones and hydroxycinamic acids and other organic acids. The Shiraz data show a greater number of significant correlations between metabolites than CS (Additional file 2: Table S4A & S4B, Additional file 1: Figure S4A).
Metabolic networks generated from matrices of significant correlations (Q < 0.05 and 0.7 < r < −0.7) built for Shiraz and CS indicated an overall higher number of edges (1957 vs 1771, respectively), average nodal degree (39.5 vs 34.7, respectively) (Additional file 2: Table S4C) and frequency of the nodes with large degree in Shiraz network compared to CS (Additional file 1: Figure S4B). Taking together, the lower diameter of the Shiraz network (Additional file 2: Table S4C) suggests a tighter coordinated behavior of metabolism during the development of Shiraz compared to CS. This suggestion is also supported by higher transitivity score of the Shiraz network (Additional file 2: Table S4C), i.e. the probability of a network to form clusters with stronger interconnections.
Transcriptomic changes during berry development
Sampling berries with similar developmental stages can be difficult due to variability even between clusters on the same vine and among berries within clusters. Precautions were taken to collect representative and average (in size and appearance) berries from the middle of the cluster to reduce such developmental variations. Veraison is a transitional phase from green stage to ripening when massive metabolite changes towards ripening begin ,,, and this transition brings a fundamental transcriptomic reprograming of the entire vine  that determines the final metabolite composition of the berry. Taking into consideration the results of the whole-berry weight and biochemical parameters and to minimize the developmental element in the variation of gene expression between the two cultivars, samples for RNA-seq data analysis were harvested at veraison and 20 days after veraison (DAV), i.e. maximal anthocyanin content (Figure 1). At both developmentally and biochemically defined stages, the two cultivars displayed similar oBRIX levels, which made the comparative investigation of dynamics of transcriptional and metabolite changes possible . Following identification of significantly changed transcripts, to facilitate the comparison and visualization of transcriptomic changes in both cultivars, the data were normalized to the CS veraison stage transcript level to visualize the change between stages and between cultivars.
Transcripts putatively encoding for TCA cycle enzymes such as citrate synthase (VIT_13s0156g00110) aconitase (VIT_05s0049g01980), isocitrate dehydrogenase (VIT_14s0066g00950), succinyl-CoA synthetase (VIT_10s0042g00950 and VIT_01s0127g00260) and fumarate hydratase (VIT_07s0005g00820) were downregulated following veraison (Figure 6), consistent with the metabolite profiles of these TCA intermediates.
Differentially accumulated transcripts involved in the phenylpropanoid and flavonoid pathway are depicted in Figure 7. Interestingly genes involved in the hydroxylation of dihydrokaempferol via putative flavonoid 3′,5′ hydroxylase (F3′5′H) (VIT_06s0009g02970 and VIT_06s0009g02840) that lead to anthocyanin precursor, delphinidin, significantly increased as the berry developed from veraison to post-veraison stage in both cultivars. In contrast, the parallel biosynthetic pathway that utilizes the same precursor, dihydrokaempferol to another anthocyanin precursor (cyanidin) via the enzyme flavonoid 3′ hydroxylase (F3′H) did not show significant increased expression in Shiraz while being downregulated in CS from veraison to mid-ripening stage (Additional file 2: Table S5). Previous studies showed that high expression of F3′H and F3′5′H correlated with the composition of flavonols, flavanols and some anthocyanins ,. Hence, F3′H as an early acting flavonoid gene can provide precursor for flavanol and cyanidin type of anthocyanin whereas F3′5′H might divert the flux to delphinidin derivatives at later stages of berry development. Thus, the significantly higher expression of F3′5′H transcripts at postharvest is presumably involved in polyphenol maturation. Our analysis of metabolites and transcripts suggest a consistent association between metabolites and related gene transcripts in the polyphenol biosynthesis (Figures 5 and 7).
The cuticular membrane of grape berry, which is composed of a wax layer that protects the fruit from desiccation and injuries. Three putative wax synthase-coding transcripts (VIT_15s0046g00590, VIT_15s0046g00520 and VIT_03s0063g00050) decreased significantly in abundance from veraison, while three other wax synthase-coding transcripts (VIT_15s0046g00470, VIT_19s0090g01340 and VIT_19s0090g01350) were induced significantly as the berry developed from veraison to the mid-ripening stage in CS; this was not so in Shiraz berry skins, where a down regulation of two of the wax-synthesis related genes (VIT_15s0046g00590 and VIT_15s0046g00520) was detected and only one was significantly upregulated (VIT_15s0046g00470) (Additional file 2: Table S5). Our results partly contrast with previous microarray transcriptomic studies where maximum wax synthase-coding transcripts were shown at veraison ,. The two cultivars exhibited differential expression in two different transcripts coding for universal stress protein (USP). Shiraz exhibited greater changes (3.2- and 2.4-fold changes for VIT_06s0004g05670 and VIT_08s0007g01360, respectively) compared to changes in CS (2.4- and 1.7-fold changes for VIT_07s0005g01290 and VIT_08s0007g01360, respectively) (Additional file 2: Table S5). These data are in line with previous results where we showed a more pronounced drought related response in Shiraz to progressive water deficit compared to CS .
Expectedly, many phytohormones related genes involved in the berry ripening process showed changes in transcription level (Additional file 2: Table S6) in both cultivars as the berry developed from veraison to post-veraison stage. For example, a total of 19 genes related to ABA metabolism showed differential expression in development and cultivar specificity, i.e. an upregulation of ABA related genes was shown in Shiraz, whereas a repression of the same process was evident in CS (Additional file 2: Table S6). The expression of genes related with ethylene biosynthesis also highly changed in expression during development while showing similar trends in the two cultivars (Additional file 2: Table S6). Similarly, brassinosteroid biosynthesis related genes were shown to change significantly in the two cultivars but to different extents. Namely, in CS, two genes upregulated and eight genes downregulated while in Shiraz three genes were upregulated while five were downregulated. In relation to jasmonate biosynthesis, the two cultivars upregulated one gene, while showing a greater downregulation of the process in CS (five genes) than Shiraz (two genes, Additional file 2: Table S6).
Gene enrichment analysis
The results of functional classification on classified genes using Expander enrichment analysis tool are shown in Additional file 1: Figure S5. The enrichment analysis yielded 44 and 37 significantly enriched networks in CS and Shiraz, respectively. The detail list of significantly enriched networks (Bonferroni FDR corrected p- value cutoff: 0.05) and member of gene lists with a further partitioning into upregulated and downregulated genes, is given in Additional file 2: Table S7. Among the classified genes, Shiraz exhibited higher enrichment in carbohydrate and specialized metabolic processes, whereas lipid, developmental and DNA/RNA processes were over represented in CS berries (Additional file 1: Figure S5).
Cultivar differences in berry metabolism and its regulation by endogenous and environmental cues during development give rise to the distinct varieties in wine flavor ,,-. Light, temperature and water demand are significant factors affecting plant growth and development - . In this respect, Cabernet Sauvignon and Shiraz were recently shown to harbor different hydraulic behavior and water use , confirmed in the field in the present study measuring the volumetric soil water content (Additional file 1: Figure S1), and responded differently to terminal drought in respect to leaf metabolism and vegetative physiology . Previously, separate investigations on berry skin microarray analysis of Cabernet Sauvignon  and RNAseq of Shiraz  reported transcriptional profile and dynamics of gene expression during the course of berry development. Here a comparative parallel analysis of metabolites and RNA transcripts by RNAseq of the two cultivars enabled a comprehensive description of the commonalities and differences in the regulation of the metabolic processes during berry development.
During the early stages of berry development the berry skin is entirely green and photosynthetically active ,, contributing locally to the biosynthesis of sugars and other organic acids , that will later be used as substrate in the specialized metabolism. Following veraison, however, a progressive decline in the sugar content, as here shown, is known to occur in support of a continuous use of carbohydrate resources produced via photosynthesis as recently reported by Dai et al. . Most of the identified organic acids showed progressive decrease in both cultivars from veraison to berry ripening. The metabolite profiles were supported by the transcript profiles indicating a down regulation of the TCA cycle in both CS and Shiraz as the berry developed from veraison to 20 DAV. Our results are in accordance with previous findings  indicating that the progressive decline of TCA cycle activity likely matches the increased demand of supply for precursors in the synthesis of carbohydrates, amino acids and subsequent flavonoids . Nevertheless a larger carbon demand may be channeled into the specialized metabolism of Shiraz and specifically in the enhanced production of anthocyanins, particularly coumaroyl forms (discussed below). For example greater sucrose accumulation was exhibited in CS than in Shiraz from veraison to maturation; lower accumulation of quinate and shikimate, precursors of the polyphenol metabolism, were measured in Shiraz during later stages of maturation, suggesting an increased integration into downstream processes (anthocyanins production). Correlation analysis further supported the evidence of highly coordinated metabolic changes between central and specialized metabolism. The occurrence of a greater number of significant correlations in the Shiraz metabolite profiling, the stronger properties of the metabolic network of Shiraz compared with CS and the relatively higher expression of flavonoid genes suggest that Shiraz metabolism has a relatively tighter coordination of the metabolic pathways investigated.
The progressive development to maturation was characterized by the increased expression of a high number of transcripts putatively encoding for phenylalanine ammonia lyase (PAL), stilbene synthase, chalcone synthase and chalcone isomerase; these were upregulated in both cultivars, but more significantly in Shiraz. Concomitantly, at the metabolite level, shikimate content was substantially reduced. An intermediate metabolite that links primary metabolism to the polyphenol biosynthesis ,, shikimate is the precursor of aromatic amino acid, phenylalanine, which was decreased in relative abundance towards maturity in coordination with higher content of downstream polyphenols. Changes in metabolite (Phe) and related transcripts in Shiraz were greater than CS, indicating a coordinated shift of central metabolism fueling the biosynthesis of specialized compounds in both cultivars, but to a greater extent in Shiraz. The continuous supply of nutrients during ripening via phloem from matured leaves fuels berry metabolic processes . Nevertheless, it is tempting to hypothesize that the commencement of specialized metabolism in the skin at veraison is supported by the biosynthesis and accumulation of primary substrates, e.g. sugars and organic acids including shikimate, via localized production of photoassimilates during the pre-veraison stages. The role of photosynthesis in the fleshy fruits was addressed recently in developing tomatoes . The results showed that at least under the conditions used fruit photosynthesis was not necessary for fruit energy metabolism or development, but it was necessary for properly timed seed development. The role of photosynthesis in organs shifting from autotrophic to hetrotrophic tissue was investigated in green seeds too. Here the results showed that photosynthesis contributed to maintain energy and oxygen levels over a metabolic-promotive threshold in a tissue that subsequently would become dense and hypoxic . These contradicting results prompt for more work to elucidate the role of localized photosynthesis to tissue development. Berry maturation was also accompanied by enhanced stress related metabolism such as trehalose, stilbene and ABA suggesting for a more susceptible Shiraz to environmental cues, likely related with the conditions in the field, reflected in high radiation and temperature. In support of this hypothesis in a recent study under progressive drought conditions Shiraz responded with greater alterations in the central metabolism of the leaf and by the accumulation of stress related compounds .
The two cultivars exhibited similar patterns in transcripts and metabolites of the flavonoid pathway, but differed in their magnitude of change
Anthocyanins make up a significant portion of berry skin flavonoids  and five classes (cyanidin, peonidin, delphinidin, petunidin and malvidin) are commonly found in Vitis vinifera grape with their glucosides, acetyl-glucosides, coumaroyl-glucosides, peonidin and malvidin caffeoyl-glucosides . Dihydrokaempferol is a common substrate for F3′H and F3′5′H to produce their corresponding cyanidin-type and delphinidin-type anthocyanins . The two enzymes, F3′H and F3′5′H, belong to the cytochrome p450 protein family  and compete a common precursor for the synthesis of red and blue anthocyanins, respectively . Previous phylogenetic analysis indicated that F3′5′H was recruited from F3′H before the divergence of angiosperms and gymnosperms . In the current study transcripts coding for F3′5′H increased significantly in relative abundance from veraison to post veraison stages in both cultivars. The increase in F3′5′H transcript abundance was coupled to a decrease of F3′H transcript abundance and with higher accumulation of delphinidin-type anthocyanins in the downstream steps of this pathway. Presumably F3′5′H plays a key role in driving the phenylpropanoid substrate flux towards delphinidin-type anthocyanin pathways. Mattivi et al.  also suggested that higher F3′5′OH activity would direct anthocyanin modification towards delphinidin-type anthocyanin forms. Liu et al.  demonstrated the greater increase in transcript abundance of F3′5′H genes under UV treatment compared with the changes in transcript abundance of F3′H genes in Antarctic moss. Perhaps, F3′5′H genes are highly associated with strategies to cope with increased stress during berry ripening.
In the present study the coumaroyl form of anthocyanins were accumulated in Shiraz to a significantly higher extent compared to CS (Additional file 2: Table S8) towards the later stage of berry development and this coincided with the increase in the level of phenylpropanoid-related gene transcripts. The acylation step of anthocyanin biosynthesis is mediated by anthocyanin acyltransferase which links aromatic constituents to the C6′ position of the glucosyl group . This step of biosynthesis promotes color stability and intensity of anthocyanins -. Our results are in agreement with substantial cultivar variability in coumaroyl forms of anthocyanins as previously reported by Mazza et al.  and Boss et al. .
Shiraz showed enhanced stilbene metabolism
A significant cultivar specificity was observed in the stilbene metabolism and regulation. Specifically, Shiraz displayed a significant upregulation of many STS genes during ripening, which was not observed in CS. The transcript profile was in accordance with the content of resveratrol and piceid, which accumulated in Shiraz more than in CS. These findings indicate a differential regulation of stilbene metabolism between the two cultivars and confirm the high correlation of STS transcript level and stilbene compounds abundance in grape berry skin ,. Having said that, the branching between stilbene and flavonoid metabolism may imply a competition for the same precursors. However the concomitant accumulation of anthocyanins in Shiraz suggests that the entire polyphenol pathway is upregulated in this cultivar compared with CS, possibly as a result to local climate conditions (expressed by the stress related genes and metabolites) and at the expense of primary C substrates (e.g. sucrose).
Shiraz exhibited enhanced transcription of phytohormone related genes
The role of phytohormones in relation to berry ripening with their coordinated interaction with other signal molecules has been reported previously . Exogenous ABA application induces the ripening process in climacteric fruits such as tomato , and banana . Endogenous ABA was associated with sugar metabolism in water stressed grapevine . Previous work demonstrated that the level of ABA in the berry increase towards ripening ,, and it is likely involved in enhanced phenolic compound biosynthesis -. Our data indicates upregulation of ABA related genes in Shiraz and indeed associated with varietal specific enhanced anthocyanin accumulation and specifically coumaroyl anthocyanin forms. Cramer et al.  demonstrated that among the genes related to hormone metabolism, ABA related transcripts changed the most in response to water and salinity stresses. Given the recently established differences in stress response between Shiraz and CS ,, a possible explanation for the enhanced ABA-related genes and other stress related processes, could be related to magnitude of stress experienced by that cultivar under the condition of growth in the experimental field, e.g. high evapotranspiration and radiation and low humidity.
The present comparative study showed a tight coordination between polyphenol metabolism and respective transcripts in the developing berries of the two cultivars. The inverse regulation of the two committing steps of phenylpropanoid metabolism, F3′H and F3′5′H, sharing a common substrate, was reflected by the metabolite profiles of the two parallel downstream anthocyanin pathways. Correlation analysis showed a “source sink”- like relation between central and specialized metabolism, also supported by the transcript analysis. At the variety level, Shiraz was characterized by a more synchronized developmental metabolism as well as by a greater upregulation of the entire polyphenol pathway, stress related processes. Taken together these data indicate that the varietal susceptibility to the environment may be influenced by the vines hydraulic behavior and physiology, and may determine important metabolic traits at fruit harvest.
Availability of supporting data
Supporting data for list of metabolite standards, principle component analyses, metabolite correlation analysis, list of significantly changed transcripts, list of phytohormones related genes, and enrichment analysis, metabolite pairwise comparison of samples and statistical analysis of coumaroyl anthocyanin forms are available in Additional file 2: Tables S1–S8. The change in volumetric water content of the soil, relative abundance of hydroxycinamate, relative abundance of coumaroyl anthocyanin forms, summary of correlation analysis and functional category distribution are presented in Additional file 1: Figures S1–S5.
AF and GC are indebted with the BARD (The United States – Israel Binational Agriculture and Development Fund) foundation for supporting this study (grant number 814512). AD, UH and RG would like to thank the GRCN ((NSF grant number DBI 0741876)) and the Goldinger foundations for financial support. AF and FM acknowledge the support of the Blaustein Center for Scientific Cooperation (BCSC) to their collaboration.
- Arroyo-Garcia R, Ruiz‐Garcia L, Bolling L, Ocete R, Lopez M, Arnold C, Ergul A, Uzun H, Cabello F, Ibáñez J: Multiple origins of cultivated grapevine (Vitis vinifera L. ssp. sativa) based on chloroplast DNA polymorphisms. Mol Ecol. 2006, 15 (12): 3707-3714.PubMedGoogle Scholar
- Webb L, Whetton P, Barlow E: Climate change and winegrape quality in Australia. Climate Res. 2008, 36 (2): 99-Google Scholar
- This P, Lacombe T, Thomas MR: Historical origins and genetic diversity of wine grapes. Trends Genet. 2006, 22 (9): 511-519.PubMedGoogle Scholar
- Alonso-Blanco C, Aarts MG, Bentsink L, Keurentjes JJ, Reymond M, Vreugdenhil D, Koornneef M: What has natural variation taught us about plant development, physiology, and adaptation?. Plant Cell. 2009, 21 (7): 1877-1896.PubMed CentralPubMedGoogle Scholar
- Mabrouk H, Sinoquet H: Indices of light microclimate and canopy structure of grapevines determined by 3D digitising and image analysis, and their relationship to grape quality. Aust J Grape Wine Res. 1998, 4 (1): 2-13.Google Scholar
- Percival D, Fisher K, Sullivan J: Use of fruit zone leaf removal with Vitis vinifera L. cv. Riesling grapevines. I. Effects on canopy structure, microclimate, bud survival, shoot density, and vine vigor. Am J Enol Vitic. 1994, 45 (2): 123-132.Google Scholar
- Louarn G, Dauzat J, Lecoeur J, Lebon E: Influence of trellis system and shoot positioning on light interception and distribution in two grapevine cultivars with different architectures: an original approach based on 3D canopy modelling. Aust J Grape Wine Res. 2008, 14 (3): 143-152.Google Scholar
- Delgado R, Martín P, del Álamo M, González M: Changes in the phenolic composition of grape berries during ripening in relation to vineyard nitrogen and potassium fertilisation rates. J Sci Food Agric. 2004, 84 (7): 623-630.Google Scholar
- Downey MO, Harvey JS, Robinson SP: The effect of bunch shading on berry development and flavonoid accumulation in Shiraz grapes. Aust J Grape Wine Res. 2004, 10 (1): 55-73.Google Scholar
- de Orduña MR: Climate change associated effects on grape and wine quality and production. Food Res Int. 2010, 43 (7): 1844-1855.Google Scholar
- Coombe BG, Hale C: The hormone content of ripening grape berries and the effects of growth substance treatments. Plant Physiol. 1973, 51 (4): 629-634.PubMed CentralPubMedGoogle Scholar
- Coombe B: Research on development and ripening of the grape berry. Am J Enol Vitic. 1992, 43 (1): 101-Google Scholar
- Waters DLE, Holton TA, Ablett EM, Lee LS, Henry RJ: cDNA microarray analysis of developing grape (Vitis vinifera cv. Shiraz) berry skin. Funct Integr Genomics. 2005, 5 (1): 40-58.PubMedGoogle Scholar
- Deluc L, Bogs J, Walker AR, Ferrier T, Decendit A, Merillon JM, Robinson SP, Barrieu F: The transcription factor VvMYB5b contributes to the regulation of anthocyanin and proanthocyanidin biosynthesis in developing grape berries. Plant Physiol. 2008, 147 (4): 2041-PubMed CentralPubMedGoogle 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 (1): 429-PubMed CentralPubMedGoogle 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 (Global Science books). 2007, 1 (1): 1-22.Google Scholar
- Kennedy JA, Matthews MA, Waterhouse AL: Changes in grape seed polyphenols during fruit ripening. Phytochemistry. 2000, 55 (1): 77-85.PubMedGoogle Scholar
- Kennedy JA, Troup GJ, Pilbrow JR, Hutton DR, Hewitt D, Hunter CR, Ristic R, Iland PG, Jones GP: Development of seed polyphenols in berries from Vitis vinifera L. cv. Shiraz. Aust J Grape Wine Res. 2000, 6 (3): 244-254.Google Scholar
- Kennedy JA, Hayasaka Y, Vidal S, Waters EJ, Jones GP: Composition of grape skin proanthocyanidins at different stages of berry development. J Agric Food Chem. 2001, 49 (11): 5348-5355.PubMedGoogle Scholar
- Robinson SP, DAVIES C: Molecular biology of grape berry ripening. Aust J Grape Wine Res. 2000, 6 (2): 175-188.Google Scholar
- 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.PubMedGoogle Scholar
- Katalinić V, Možina SS, Skroza D, Generalić I, Abramovič H, Miloš M, Ljubenkov I, Piskernik S, Pezo I, Terpinc P: Polyphenolic profile, antioxidant properties and antimicrobial activity of grape skin extracts of 14 < i > Vitis vinifera varieties grown in Dalmatia (Croatia). Food Chem. 2010, 119 (2): 715-723.Google Scholar
- Liang Z, Wu B, Fan P, Yang C, Duan W, Zheng X, Liu C, Li S: Anthocyanin composition and content in grape berry skin in < i > Vitis germplasm. Food Chem. 2008, 111 (4): 837-844.Google Scholar
- Jones GV, White MA, Cooper OR, Storchmann K: Climate change and global wine quality. Clim Change. 2005, 73 (3): 319-343.Google Scholar
- Kliewer WM: Influence of environment on metabolism of organic acids and carbohydrates in Vitis vinifera. I Temperature. Plant Physiol. 1964, 39 (6): 869-PubMed CentralPubMedGoogle Scholar
- White MA, Diffenbaugh N, Jones GV, Pal J, Giorgi F: Extreme heat reduces and shifts United States premium wine production in the 21st century. Proc Natl Acad Sci. 2006, 103 (30): 11217-11222.PubMed CentralPubMedGoogle Scholar
- Pomar F, Novo M, Masa A: Varietal differences among the anthocyanin profiles of 50 red table grape cultivars studied by high performance liquid chromatography. J Chromatogr A. 2005, 1094 (1): 34-41.PubMedGoogle Scholar
- Bravdo B, Naor A, Zahavi T, Gal Y, Snyder R: The effect of water stress applied alternately to part of the wetting zone along the season (PRD-partial rootzone drying) on wine quality, yield and water relations of red wine grapes. Acta Hortic. 2004, 664: 101-109.Google Scholar
- Hardie W, Considine J: Response of grapes to water-deficit stress in particular stages of development. Am J Enol Vitic. 1976, 27 (2): 55-61.Google Scholar
- Matthews MA, Anderson MM: Am J Enol Vitic. 1989, 40 (1): 52-60.Google Scholar
- Ojeda H, Andary C, Kraeva E, Carbonneau A, Deloire A: Influence of pre-and postveraison water deficit on synthesis and concentration of skin phenolic compounds during berry growth of Vitis vinifera cv. Shiraz. Am J Enol Vitic. 2002, 53 (4): 261-267.Google Scholar
- Hunter J, Ruffner H, Volschenk C, Le Roux D: Partial defoliation of Vitis vinifera L. cv. Cabernet Sauvignon/99 Richter: effect on root growth, canopy efficiency; grape composition, and wine quality. Am J Enol Vitic. 1995, 46 (3): 306-314.Google Scholar
- Ollat N, Gaudillere J: The effect of limiting leaf area during stage I of berry growth on development and composition of berries of Vitis vinifera L. cv. Cabernet Sauvignon. Am J Enol Vitic. 1998, 49 (3): 251-258.Google Scholar
- 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.Google Scholar
- Zohary D: Unconscious selection and the evolution of domesticated plants. Econ Bot. 2004, 58 (1): 5-10.Google Scholar
- Lemaire-Chamley M, Petit J, Garcia V, Just D, Baldet P, Germain V, Fagard M, Mouassite M, Cheniclet C, Rothan C: Changes in transcriptional profiles are associated with early fruit tissue specialization in tomato. Plant Physiol. 2005, 139 (2): 750-769.PubMed CentralPubMedGoogle Scholar
- Hichri I, Barrieu F, Bogs J, Kappel C, Delrot S, Lauvergeat V: Recent advances in the transcriptional regulation of the flavonoid biosynthetic pathway. J Exp Bot. 2011, 62 (8): 2465-2483.PubMedGoogle Scholar
- Kolb CA, Kopecký J, Riederer M, Pfündel EE: UV screening by phenolics in berries of grapevine (Vitis vinifera). Funct Plant Biol. 2003, 30 (12): 1177-1186.Google Scholar
- Dixon RA, Achnine L, Kota P, Liu C, Reddy M, Wang L: The phenylpropanoid pathway and plant defence—a genomics perspective. Mol Plant Pathol. 2002, 3 (5): 371-390.PubMedGoogle Scholar
- Winkel-Shirley B: Biosynthesis of flavonoids and effects of stress. Curr Opin Plant Biol. 2002, 5 (3): 218-223.PubMedGoogle Scholar
- Adams DO: Phenolics and ripening in grape berries. Am J Enol Vitic. 2006, 57 (3): 249-256.Google Scholar
- Rice-Evans CA, Miller NJ, Paganga G: Structure-antioxidant activity relationships of flavonoids and phenolic acids. Free Radic Biol Med. 1996, 20 (7): 933-956.PubMedGoogle Scholar
- Rice-Evans C, Miller N, Paganga G: Antioxidant properties of phenolic compounds. Trends Plant Sci. 1997, 2 (4): 152-159.Google Scholar
- Braidot E, Petrussa E, Bertolini A, Peresson C, Ermacora P, Loi N, Terdoslavich M, Passamonti S, Macrì F, Vianello A: Evidence for a putative flavonoid translocator similar to mammalian bilitranslocase in grape berries (Vitis vinifera L.) during ripening. Planta. 2008, 228 (1): 203-213.PubMedGoogle Scholar
- Jaillon O, Aury J, Noel B, Policriti A, Clepet C, Casagrande A, Choisne N, Aubourg S, Vitulo N, Jubin C: The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007, 449 (7161): 463-467.PubMedGoogle Scholar
- Ali M, Howard S, Chen S, Wang Y, Yu O, Kovacs L, Qiu W: Berry skin development in Norton grape: distinct patterns of transcriptional regulation and flavonoid biosynthesis. BMC Plant Biol. 2011, 11 (1): 7-PubMed CentralPubMedGoogle Scholar
- Lijavetzky D, Carbonell-Bejerano P, Grimplet J, Bravo G, Flores P, Fenoll J, Hellín P, Oliveros JC, Martínez-Zapater JM: Berry flesh and skin ripening features in Vitis vinifera as assessed by transcriptional profiling. PLoS One. 2012, 7 (6): e39547-PubMed CentralPubMedGoogle Scholar
- Jeong S, Goto-Yamamoto N, Hashizume K, Esaka M: Expression of multi-copy flavonoid pathway genes coincides with anthocyanin, flavonol and flavan-3-ol accumulation of grapevine. Vitis. 2008, 47 (3): 135-140.Google Scholar
- Pilati S, Perazzolli M, Malossini A, Cestaro A, Demattè 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 (1): 428-PubMed CentralPubMedGoogle Scholar
- Terrier N, Glissant D, Grimplet J, Barrieu F, Abbal P, Couture C, Ageorges A, Atanassova R, Léon C, Renaudin J: Isogene specific oligo arrays reveal multifaceted changes in gene expression during grape berry (Vitis vinifera L.) development. Planta. 2005, 222 (5): 832-847.PubMedGoogle Scholar
- Da Silva C, Zamperin G, Ferrarini A, Minio A, Dal Molin A, Venturini L, Buson G, Tononi P, Avanzato C, Zago E, Boido E, Dellacassa E, Gaggero C, Pezzotti M, Carrau F, Delledonne M: The high polyphenol content of grapevine cultivar tannat berries is conferred primarily by genes that are not shared with the reference genome. Plant Cell. 2013, 25 (12): 4777-4788.PubMed CentralPubMedGoogle Scholar
- Venturini L, Ferrarini A, Zenoni S, Tornielli GB, Fasoli M, Dal Santo S, Minio A, Buson G, Tononi P, Zago ED: De novo transcriptome characterization of Vitis vinifera cv. Corvina unveils varietal diversity. BMC Genomics. 2013, 14 (1): 41-PubMed CentralPubMedGoogle Scholar
- Di Genova A, Almeida AM, Mun C, Vizoso P, Travisany D, Moraga C, Pinto M, Hinrichsen P, Orellana A, Maass A: Whole genome comparison between table and wine grapes reveals a comprehensive catalog of structural variants. BMC Plant Biol. 2014, 14 (1): 7-PubMed CentralPubMedGoogle Scholar
- Wang Z, Gerstein M, Snyder M: RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009, 10 (1): 57-63.PubMed CentralPubMedGoogle Scholar
- Zenoni S, Ferrarini A, Giacomelli E, Xumerle L, Fasoli M, Malerba G, Bellin D, Pezzotti M, Delledonne M: Characterization of transcriptional complexity during berry development in Vitis vinifera using RNA-Seq. Plant Physiol. 2010, 152 (4): 1787-1795.PubMed CentralPubMedGoogle Scholar
- Fasoli M, Dal Santo S, Zenoni S, Tornielli GB, Farina L, Zamboni A, Porceddu A, Venturini L, Bicego M, Murino V: The grapevine expression atlas reveals a deep transcriptome shift driving the entire plant into a maturation program. Plant Cell. 2012, 24 (9): 3489-3505.PubMed CentralPubMedGoogle Scholar
- Sweetman C, Wong DC, Ford CM, Drew DP: Transcriptome analysis at four developmental stages of grape berry (Vitis vinifera cv. Shiraz) provides insights into regulated and coordinated gene expression. BMC Genomics. 2012, 13 (1): 691-PubMed CentralPubMedGoogle Scholar
- Mattoo AK, Chung SH, Goyal RK, Fatima T, Solomos T, Srivastava A, Handa AK: Overaccumulation of higher polyamines in ripening transgenic tomato fruit revives metabolic memory, upregulates anabolism-related genes, and positively impacts nutritional quality. J AOAC Int. 2007, 90 (5): 1456-1464.PubMedGoogle Scholar
- Bramley PM: Regulation of carotenoid formation during tomato fruit ripening and development. J Exp Bot. 2002, 53 (377): 2107-2113.PubMedGoogle Scholar
- Rosati C, Aquilani R, Dharmapuri S, Pallara P, Marusic C, Tavazza R, Bouvier F, Camara B, Giuliano G: Metabolic engineering of beta-carotene and lycopene content in tomato fruit. Plant J. 2000, 24 (3): 413-420.PubMedGoogle Scholar
- Castellarin SD, Matthews MA, Di Gaspero G, Gambetta GA: Water deficits accelerate ripening and induce changes in gene expression regulating flavonoid biosynthesis in grape berries. Planta. 2007, 227 (1): 101-112.PubMedGoogle Scholar
- Castellarin SD, Di Gaspero G, Marconi R, Nonis A, Peterlunger E, Paillard S, Adam-Blondon A, Testolin R: Colour variation in red grapevines (Vitis vinifera L.): genomic organisation, expression of flavonoid 3′-hydroxylase, flavonoid 3′, 5′-hydroxylase genes and related metabolite profiling of red cyanidin-/blue delphinidin-based anthocyanins in berry skin. BMC Genomics. 2006, 7 (1): 12-PubMed CentralPubMedGoogle Scholar
- Pereira G, Gaudillere J, Pieri P, Hilbert G, Maucourt M, Deborde C, Moing A, Rolin D: Microclimate influence on mineral and metabolic profiles of grape berries. J Agric Food Chem. 2006, 54 (18): 6765-6775.PubMedGoogle Scholar
- Pereira GE, Gaudillere J, Leeuwen C, Hilbert G, Maucourt M, Deborde C, Moing A, Rolin D: < sup> 1 H NMR metabolite fingerprints of grape berry: Comparison of vintage and soil effects in Bordeaux grapevine growing areas. Anal Chim Acta. 2006, 563 (1): 346-352.Google Scholar
- Dai ZW, Léon C, Feil R, Lunn JE, Delrot S, Gomès 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 (5): 1345-1355.PubMed CentralPubMedGoogle Scholar
- Guymon J, Ough C: A uniform method for total acid determination in wines. Am J Enol Vitic. 1962, 13 (1): 40-45.Google Scholar
- Weckwerth W, Wenzel K, Fiehn O: Process for the integrated extraction, identification and quantification of metabolites, proteins and RNA to reveal their co-regulation in biochemical networks. Proteomics. 2004, 4 (1): 78-83.PubMedGoogle Scholar
- Lisec J, Schauer N, Kopka J, Willmitzer L, Fernie AR: Gas chromatography mass spectrometry–based metabolite profiling in plants. Nat Protoc. 2006, 1 (1): 387-396.PubMedGoogle Scholar
- Hochberg U, Degu A, Toubiana D, Gendler T, Nikoloski Z, Rachmilevitch S, Fait A: Metabolite profiling and network analysis reveal coordinated changes in grapevine water stress response. BMC Plant Biol. 2013, 13 (1): 184-PubMed CentralPubMedGoogle Scholar
- Arapitsas P, Scholz M, Vrhovsek U, Di Blasi S, Bartolini AB, Masuero D, Perenzoni D, Rigo A, Mattivi F: A metabolomic approach to the study of wine micro-oxygenation. PLoS One. 2012, 7 (5): e37783-PubMed CentralPubMedGoogle Scholar
- Monagas M, Suárez R, Gómez-Cordovés C, Bartolomé B: Simultaneous determination of nonanthocyanin phenolic compounds in red wines by HPLC-DAD/ESI-MS. Am J Enol Vitic. 2005, 56 (2): 139-147.Google Scholar
- Monagas M, Núñez V, Bartolomé B, Gómez-Cordovés C: Anthocyanin-derived pigments in Graciano, Tempranillo, and Cabernet Sauvignon wines produced in Spain. Am J Enol Vitic. 2003, 54 (3): 163-169.Google Scholar
- Sánchez‐Rabaneda F, Jauregui O, Lamuela‐Raventós RM, Viladomat F, Bastida J, Codina C: Qualitative analysis of phenolic compounds in apple pomace using liquid chromatography coupled to mass spectrometry in tandem mode. Rapid Commun Mass Spectrom. 2004, 18 (5): 553-563.PubMedGoogle Scholar
- Iijima Y, Nakamura Y, Ogata Y, Tanaka K, Sakurai N, Suda K, Suzuki T, Suzuki H, Okazaki K, Kitayama M: Metabolite annotations based on the integration of mass spectral information. Plant J. 2008, 54 (5): 949-962.PubMed CentralPubMedGoogle Scholar
- Hanhineva K, Rogachev I, Kokko H, Mintz-Oron S, Venger I, Kärenlampi S, Aharoni A: Non-targeted analysis of spatial metabolite composition in strawberry (<i > Fragaria × < i > ananassa) flowers. Phytochemistry. 2008, 69 (13): 2463-2481.PubMedGoogle Scholar
- Moco S, Bino RJ, Vorst O, Verhoeven HA, de Groot J, van Beek TA, Vervoort J, De Vos CR: A liquid chromatography-mass spectrometry-based metabolome database for tomato. Plant Physiol. 2006, 141 (4): 1205-1218.PubMed CentralPubMedGoogle Scholar
- Toubiana D, Semel Y, Tohge T, Beleggia R, Cattivelli L, Rosental L, Nikoloski Z, Zamir D, Fernie AR, Fait A: Metabolic profiling of a mapping population exposes new insights in the regulation of seed metabolism and seed, fruit, and plant relations. PLoS Genet. 2012, 8 (3): e1002612-PubMed CentralPubMedGoogle Scholar
- Japelaghi RH, Haddad R, Garoosi G: Rapid and efficient isolation of high quality nucleic acids from plant tissues rich in polyphenols and polysaccharides. Mol Biotechnol. 2011, 49 (2): 129-137.PubMedGoogle Scholar
- Surget-Groba Y, Montoya-Burgos JI: Optimization of de novo transcriptome assembly from next-generation sequencing data. Genome Res. 2010, 20 (10): 1432-1440.PubMed CentralPubMedGoogle Scholar
- Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L: Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010, 28 (5): 511-515.PubMed CentralPubMedGoogle Scholar
- Roberts A, Trapnell C, Donaghey J, Rinn JL, Pachter L: Improving RNA-Seq expression estimates by correcting for fragment bias. Genome Biol. 2011, 12 (3): R22-PubMed CentralPubMedGoogle Scholar
- Grimplet J, Van Hemert J, Carbonell-Bejerano P, Díaz-Riquelme J, Dickerson J, Fennell A, Pezzotti M, Martínez-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 (1): 213-PubMed CentralPubMedGoogle Scholar
- Ulitsky I, Maron-Katz A, Shavit S, Sagir D, Linhart C, Elkon R, Tanay A, Sharan R, Shiloh Y, Shamir R: Expander: from expression microarrays to networks and functions. Nat Protoc. 2010, 5 (2): 303-322.PubMedGoogle Scholar
- Paley SM, Karp PD: The pathway tools cellular overview diagram and omics viewer. Nucleic Acids Res. 2006, 34 (13): 3771-3778.PubMed CentralPubMedGoogle Scholar
- Dharmadhikari M: Composition of grapes. Vineyard Vintage View Mo State Univ. 1994, 9 (7/8): 3-8.Google Scholar
- Fung RW, Gonzalo M, Fekete C, Kovacs LG, He Y, Marsh E, McIntyre LM, Schachtman DP, Qiu W: Powdery mildew induces defense-oriented reprogramming of the transcriptome in a susceptible but not in a resistant grapevine. Plant Physiol. 2008, 146 (1): 236-249.PubMed CentralPubMedGoogle Scholar
- Schnee S, Viret O, Gindro K: Role of stilbenes in the resistance of grapevine to powdery mildew. Physiol Mol Plant Pathol. 2008, 72 (4): 128-133.Google Scholar
- Langcake P, Pryce R: The production of resveratrol by < i > Vitis vinifera and other members of the Vitaceae as a response to infection or injury. Physiol Plant Pathol. 1976, 9 (1): 77-86.Google Scholar
- Adrian M, Jeandet P, Veneau J, Weston LA, Bessis R: Biological activity of resveratrol, a stilbenic compound from grapevines, against Botrytis cinerea, the causal agent for gray mold. J Chem Ecol. 1997, 23 (7): 1689-1702.Google Scholar
- Bézier A, Lambert B, Baillieul F: Study of defense-related gene expression in grapevine leaves and berries infected with Botrytis cinerea. Eur J Plant Pathol. 2002, 108 (2): 111-120.Google Scholar
- Langcake P, McCarthy W: The relationship of resveratrol production to infection ofgrapevine leaves by Botrytis cinerea. Vitis. 1979, 18: 244-253.Google Scholar
- Chiron H, Drouet A, Lieutier F, Payer H, Ernst D, Sandermann H: Gene induction of stilbene biosynthesis in Scots pine in response to ozone treatment, wounding, and fungal infection. Plant Physiol. 2000, 124 (2): 865-872.PubMed CentralPubMedGoogle Scholar
- Wang W, Tang K, Yang H, Wen P, Zhang P, Wang H, Huang W: Distribution of resveratrol and stilbene synthase in young grape plants (<i > Vitis vinifera L. cv. Cabernet Sauvignon) and the effect of UV-C on its accumulation. Plant Physiol Biochem. 2010, 48 (2): 142-152.PubMedGoogle Scholar
- Rosemann D, Heller W, Sandermann H: Biochemical plant responses to ozone II. Induction of stilbene biosynthesis in Scots pine (Pinus sylvestris L.) seedlings. Plant Physiol. 1991, 97 (4): 1280-1286.PubMed CentralPubMedGoogle Scholar
- Teixeira A, Eiras-Dias J, Castellarin SD, Gerós H: Berry phenolics of grapevine under challenging environments. Int J Mol Sci. 2013, 14 (9): 18711-18739.PubMed CentralPubMedGoogle Scholar
- Toffali K, Zamboni A, Anesi A, Stocchero M, Pezzotti M, Levi M, Guzzo F: Novel aspects of grape berry ripening and post-harvest withering revealed by untargeted LC-ESI-MS metabolomics analysis. Metabolomics. 2011, 7 (3): 424-436.Google Scholar
- Castellarin SD, Pfeiffer A, Sivilotti P, Degan M, Peterlunger E, Di Gaspero G: Transcriptional regulation of anthocyanin biosynthesis in ripening fruits of grapevine under seasonal water deficit. Plant Cell Environ. 2007, 30 (11): 1381-1399.PubMedGoogle Scholar
- Jeong S, Goto-Yamamoto N, Hashizume K, Esaka M: Expression of the flavonoid 3′-hydroxylase and flavonoid 3′, 5′-hydroxylase genes and flavonoid composition in grape (<i > Vitis vinifera). Plant Sci. 2006, 170 (1): 61-69.Google Scholar
- Hardie W, O’brien T, Jaudzems V: Morphology, anatomy and development of the pericarp after anthesis in grape, Vitis vinifera L. Aust J Grape Wine Res. 1996, 2 (2): 97-142.Google Scholar
- Vilanova M, Masa A, Tardaguila J: Evaluation of the aroma descriptors variability in Spanish grape cultivars by a quantitative descriptive analysis. Euphytica. 2009, 165 (2): 383-389.Google Scholar
- Gatto P, Vrhovsek U, Muth J, Segala C, Romualdi C, Fontana P, Pruefer D, Stefanini M, Moser C, Mattivi F: Ripening and genotype control stilbene accumulation in healthy grapes. J Agric Food Chem. 2008, 56 (24): 11773-11785.PubMedGoogle Scholar
- Yang J, Martinson TE, Liu RH: Phytochemical profiles and antioxidant activities of wine grapes. Food Chem. 2009, 116 (1): 332-339.Google Scholar
- Havaux M: Stress tolerance of photosystem II in vivo antagonistic effects of water, heat, and photoinhibition stresses. Plant Physiol. 1992, 100 (1): 424-432.PubMed CentralPubMedGoogle Scholar
- Chaves M, OsoÂrio M, OsoÂrio J, Pereira J: Non-stomatal limitation of photosynthesis under high temperature and water deficits in Lupinus albus. Photosynthica. 1992, 27: 521-528.Google Scholar
- Cornic G, Massacci A: Leaf photosynthesis under drought stress. Photosynthesis and the Environment. Edited by: Baker NR. 1996, Kluwer Academic Publishers, New York, 347-366.Google Scholar
- Hochberg U, Degu A, Fait A, Rachmilevitch S: Near isohydric grapevine cultivar displays higher photosynthetic efficiency and photorespiration rates under drought stress as compared with near anisohydric grapevine cultivar. Physiol Plant. 2012, 147 (4): 443-453.PubMedGoogle Scholar
- Blanke M, Lenz F: Fruit photosynthesis. Plant Cell Environ. 1989, 12 (1): 31-46.Google Scholar
- 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): 1329-1344.PubMedGoogle Scholar
- Agudelo-Romero P, Erban A, Sousa L, Pais MS, Kopka J, Fortes AM: Search for transcriptional and metabolic markers of grape Pre-ripening and ripening and insights into specific aroma development in three Portuguese cultivars. PLoS One. 2013, 8 (4): e60422-PubMed CentralPubMedGoogle Scholar
- Castellarin S, Bavaresco L, Falginella L, Gonçalves M, Di Gaspero G, Gerós H, Chaves M, Delrot S: Phenolics in Grape Berry and Key Antioxidants. The Biochemistry of Grape Berry. Edited by: Gerós HM, Chaves M, Delrot S. 2012, Bentham Science Publishers, The Netherlands, Bussum, 89-110.Google Scholar
- Hunter JK, Ruffner HP: Assimilate transport in grapevines-effect of phloem disruption. Aust J Grape Wine Res. 2001, 7 (3): 118-126.Google Scholar
- Lytovchenko A, Eickmeier I, Pons C, Osorio S, Szecowka M, Lehmberg K, Arrivault S, Tohge T, Pineda B, Anton MT, Hedtke B, Lu Y, Fisahn J, Bock R, Stitt M, Grimm B, Granell A, Fernie AR: Tomato fruit photosynthesis is seemingly unimportant in primary metabolism and ripening but plays a considerable role in seed development. Plant Physiol. 2011, 157 (4): 1650-1663.PubMed CentralPubMedGoogle Scholar
- Tschiersch H, Borisjuk L, Rutten T, Rolletschek H: Gradients of seed photosynthesis and its role for oxygen balancing. Biosystems. 2011, 103 (2): 302-308.PubMedGoogle Scholar
- Souquet J, Cheynier V, Brossaud F, Moutounet M: Polymeric proanthocyanidins from grape skins. Phytochemistry. 1996, 43 (2): 509-512.Google Scholar
- Bogs J, Ebadi A, McDavid D, Robinson SP: Identification of the flavonoid hydroxylases from grapevine and their regulation during fruit development. Plant Physiol. 2006, 140 (1): 279-291.PubMed CentralPubMedGoogle Scholar
- 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-2164-11-562-Google Scholar
- Seitz C, Eder C, Deiml B, Kellner S, Martens S, Forkmann G: Cloning, functional identification and sequence analysis of flavonoid 3′-hydroxylase and flavonoid 3′, 5′-hydroxylase cDNAs reveals independent evolution of flavonoid 3′, 5′-hydroxylase in the Asteraceae family. Plant Mol Biol. 2006, 61 (3): 365-381.PubMedGoogle Scholar
- Liu S, Ju J, Xia G: Identification of the flavonoid 3′-hydroxylase and flavonoid 3′, 5′-hydroxylase genes from Antarctic moss and their regulation during abiotic stress. Gene. 2014, 543 (1): 145-152.PubMedGoogle Scholar
- Mazza G, Francis F: Anthocyanins in grapes and grape products. Crit Rev Food Sci Nutr. 1995, 35 (4): 341-371.PubMedGoogle Scholar
- Anderson D, Gueffroy D, Webb A, Kepner R: Identification of acetic acid as an acylating agent of anthocyanin pigments in grapes. Phytochemistry. 1970, 9 (7): 1579-1583.Google Scholar
- Gueffroy D, Kepner R, Webb A: Acylated anthocyanin pigments in < i > Vitis vinifera grapes: Identification of malvidin-3-(6-p-coumaroyl) glucoside. Phytochemistry. 1971, 10 (4): 813-819.Google Scholar
- Tamura H, Hayashi Y, Sugisawa H, Kondo T: Structure determination of acylated anthocyanins in Muscat Bailey a grapes by homonuclear Hartmann-Hahn (HOHAHA) spectroscopy and liquid chromatography-mass spectrometry. Phytochem Anal. 1994, 5 (4): 190-196.Google Scholar
- 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.PubMedGoogle Scholar
- Boss PK, Davies C, Robinson SP: Anthocyanin composition and anthocyanin pathway gene expression in grapevine sports differing in berry skin colour. Aust J Grape Wine Res. 1996, 2 (3): 163-170.Google Scholar
- Versari A, Parpinello GP, Tornielli GB, Ferrarini R, Giulivo C: Stilbene compounds and stilbene synthase expression during ripening, wilting, and UV treatment in grape cv. Corvina. J Agric Food Chem. 2001, 49 (11): 5531-5536.PubMedGoogle Scholar
- Deluc LG, Decendit A, Papastamoulis Y, Mérillon J, Cushman JC, Cramer GR: Water deficit increases stilbene metabolism in Cabernet Sauvignon berries. J Agric Food Chem. 2010, 59 (1): 289-297.PubMed CentralPubMedGoogle Scholar
- Davies C, Böttcher C: Hormonal control of grape berry ripening. Grapevine Molecular Physiology & Biotechnology. Edited by: Roubelakis-Angelakis K. 2009, Springer, Heidelberg, London, New York, 229-261. 2Google Scholar
- Zhang M, Yuan B, Leng P: The role of ABA in triggering ethylene biosynthesis and ripening of tomato fruit. J Exp Bot. 2009, 60 (6): 1579-1588.PubMed CentralPubMedGoogle Scholar
- Alexander L, Grierson D: Ethylene biosynthesis and action in tomato: a model for climacteric fruit ripening. J Exp Bot. 2002, 53 (377): 2039-2055.PubMedGoogle Scholar
- Jiang Y, Joyce DC, Macnish AJ: Effect of abscisic acid on banana fruit ripening in relation to the role of ethylene. J Plant Growth Regul. 2000, 19 (1): 106-111.PubMedGoogle Scholar
- Deluc L, Quilici D, Decendit A, Grimplet J, Wheatley M, Schlauch K, Mérillon J, Cushman J, Cramer G: Water deficit alters differentially metabolic pathways affecting important flavor and quality traits in grape berries of Cabernet Sauvignon and Chardonnay. BMC Genomics. 2009, 10 (1): 212-PubMed CentralPubMedGoogle Scholar
- Inaba A, Ishida M, Sobajima Y: Changes in endogenous hormone concentrations during berry development in relation to the ripening of Delaware grapes. J Jap Soc Hort Sci. 1976, 45 (3): 245-252.Google Scholar
- Cawthon DL, Morris J: Relationship of seed number and maturity to berry development, fruit maturation, hormonal changes, and uneven ripening of’Concord’ (Vitis labrusca L.) grapes. J Jap Soc Hort Sci. 1982, 107: 1097-1104.Google Scholar
- Peppi MC, Fidelibus MW, Dokoozlian N: Abscisic acid application timing and concentration affect firmness, pigmentation, and color of flame seedless’ grapes. Hort Science. 2006, 41 (6): 1440-1445.Google Scholar
- Jeong S, Goto-Yamamoto N, Kobayashi S, Esaka M: Effects of plant hormones and shading on the accumulation of anthocyanins and the expression of anthocyanin biosynthetic genes in grape berry skins. Plant Sci. 2004, 167 (2): 247-252.Google Scholar
- Ban T, Ishimaru M, Kobayashi S, Shiozaki S, Goto-Yamamoto N, Horiuchi S: Abscisic acid and 2, 4-dichlorophenoxyacetic acid affect the expression of anthocyanin biosynthetic pathway genes in’Kyoho’grape berries. J Hortic Sci Biotechnol. 2003, 78 (4): 586-589.Google Scholar
- Gambetta GA, Matthews MA, Shaghasi TH, McElrone AJ, Castellarin SD: Sugar and abscisic acid signaling orthologs are activated at the onset of ripening in grape. Planta. 2010, 232 (1): 219-234.PubMed CentralPubMedGoogle Scholar
- Cramer GR, Ergül A, Grimplet J, Tillett RL, Tattersall EAR, Bohlman MC, Vincent D, Sonderegger J, Evans J, Osborne C: Water and salinity stress in grapevines: early and late changes in transcript and metabolite profiles. Funct Integr Genomics. 2007, 7 (2): 111-134.PubMedGoogle 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/4.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.