Not just shrivelling: time-series profiling of the biochemical changes in Corvina (Vitis vinifera L.) berries subjected to post-harvest withering

The grape berry withering process is often seen as a means to concentrate the constituents of the berry via water removal, however, recent molecular studies have indicated that the reprogramming of biosynthetic pathways that impart the unique aroma and flavour of the final wine also occur. Metabolic analysis was performed using GC-MS and LC-MS on Corvina berry dehydrated for up to 108 days. The temporal pattern of metabolic changes and time-series relationship between various classes of metabolites were investigated on berries sampled at 23-time points. Principal component analysis of both GC-MS and LC-MS datasets revealed three distinct phases of post-harvest withering; early (day 0–28), mid (day 35–66) and late (day 69–108) stages. Stress-associated amino acids such as proline, serine, ethanolamine, and leucine accumulated significantly during the mid and late stages of the drying process while phosphorylated glucose and fructose, sucrose, kestose and resveratrol increased massively across time. Unlike most of the identified metabolites, low molecular weight flavanols exhibited a consistent pattern of decline during the withering period. Our network analysis revealed that increased metabolic network connectivity occurred during the middle stage of withering, thus reflecting more coordinated metabolite changes while reduced network connectivity at early and late stages indicates minimal metabolite perturbation during these stages. The current prolonged berry withering experiment revealed the timing of metabolite interconversions in the central metabolism and provided critical clues that link the concentration effect, protein degradation, and the onset of stress-like conditions in drying berries.


INTRODUCTION
The production of wines from partially dried grapes in the Veneto dates back to the times of the Ancient Romans. The structure of the berry renders low resistance to dehydration and enables berries of some selected grape varieties to shrivel under a slow withering process and leads to the development of the typical organoleptic characteristics of Amarone wine (Accordini, 2013). The appellation of origin controlled and guaranteed (DOCG) "Amarone della Valpolicella" is reserved for wines produced from autochthonous local cultivars including Corvina and/or Corvinone (from 45 % to 95 %), and Rondinella (from 5 % to 30 %), while up to 25 % can come from other red cultivars (MIPAAF, 2019). The rate of berry's water loss and the extent of metabolic changes including aroma compounds in Corvina, Corvinone and Rondinella grapes are quite different during the postharvest withering process (D'Onofrio et al., 2019;Bellincontro et al., 2016).
The postharvest withering of Corvina grape berry is a common practice for the production of the Amarone wine. The monovarietal wines from Corvina produced without withering have been described to be dominated by an unripe astringency and, at the same time, by a high relative contribution for surface smoothness in their astringency profile , and they have lower polyphenol levels and ethanol in comparison with most Italian red wines (Arapitsas et al., 2020;Giacosa et al., 2021). These characteristics of the cultivar are ideal for the production of a wine, Amarone, produced from partially dried grapes. The use of this traditional process dates back to the sixth century AD (Paronetto and Dellaglio, 2011). Berry withering process is usually done after harvest as an offvine drying method (Mencarelli et al., 2010;Panceri et al., 2013). In the current Amarone wine production, a combination of natural and artificial ventilation practices transforms berries into partially dried fruit that can conserve and release desirable quality attributes during the vinification process (Bellincontro et al., 2004). Indeed the biosynthesis of aroma and flavour of the final wine involves a series of complex biosynthetic pathways (Sanz et al., 1996).
The post-harvest withering process of grapes lasts up to three months during which many witheringrelated metabolic changes undergo over time due to water loss (Kays, 1991). The final wine quality is mainly influenced by the accumulation of precursor metabolites during berry ripening and post-harvest withering. However, it has been shown that the speed of dehydration has a significant impact on the profile of volatile release in different grapevine varieties (Bellincontro et al., 2004). Among the factors that may modulate the rate of water loss are the structural arrangement and composition of the grape berry skin wax layer (Possingham et al., 1967;Dimopoulos et al., 2020), but also environmental factors such as temperature, wind speed, and relative humidity (Kays, 1991;Roy et al., 1999;Schirra et al., 2000;Mencarelli et al., 1994). All these elements could lead to metabolite compositional changes with a potentially detrimental impact on the final wine quality. Bellincontro et al. (2004) showed a reduction in volatile compounds while polyphenols show an increasing trend under a short-term withering experiment. Slaghenaufi et al. (2020) showed that the withering process is associated with the modulation of both fermentative and grape-related organic volatiles, including esters, higher alcohols, acids, C6 alcohols, terpenes and norisoprenoids. At the molecular level, water loss is associated with changes in cell wall enzyme activity (Hsiao, 1973), respiration and ethylene production (Kays, 1991). In the berry, these changes are known to take place (Costantini et al., 2006;Zamboni et al., 2008;Zenoni et al., 2016) alongside transcriptome-wide reprogramming of metabolism pathways involving phytohormones, secondary metabolites, carbohydrate, and cell wall, as well as stress-protection, respiration, and transcriptional regulatory genes (Zenoni et al., 2016). In the current study, MS-based profiling of primary and secondary metabolites was performed to characterize the withering-associated temporal pattern of metabolic changes in Corvina grapes.

Berry sampling
Fully matured V. vinifera cv. Corvina grape berries were harvested on 25/09/2013 at Valpolicella (province of Verona, Italy), in a commercial vineyard of Monte dei Galli, San Pietro in Cariano, and withered for 108 days under controlled ventilation in the warehouse of Allegrini winery, located a few kilometres (< 5 km) away from the vineyard. A total of 23 samples were taken for metabolic analysis at different stages of the withering process and for each time-point of sampling, three biological replicates were collected.
For berry weight loss measurement, 7 trays (about 40 kg of grapes) were weighted at 20 sampling points (Supplementary Table S1).

Sample preparation and extraction for metabolite profiling
At each sampling stage during the withering process, 20 berries were randomly sampled, immediately snap-frozen with liquid nitrogen, and kept at −80 °C until further analysis. The grape berries (skin, flesh and seed together) were ground under liquid nitrogen using an IKA analytical mill (Staufen, Germany) to obtain a frozen powder and ahead of metabolite extraction, the frozen powder was freeze-dried in a lyophilizer. The freeze-dried samples were used for both primary and secondary metabolite analysis, according to a previously validated protocol (Theodoridis et al., 2012).

Sample preparation and LC-MS targeted metabolomics.
An aliquot of 2 g of powder from each sample was extracted in sealed glass vials by adding 2400 µL of MeOH/H 2 O (2:1), 1600 µL of chloroform and 40 µL of internal standard (gentisic acid and rosmarinic acid, 500 mg/L). The rosmarinic acid standard was used to adjust for the recovery, while gentisic acid was used to verify the stability within the sequence. The vials were put in an orbital shaker for 15 min at room temperature. Samples were centrifuged at 4000 g and 4 °C for 5 min, and the upper phase consisting of aqueous methanol was collected in a volumetric flask of 10 mL. Extraction was repeated by adding another 2400 µL of MeOH/H 2 O) to the pellets and chloroform fractions. The upper phases from the second extraction were also added to the flask, the final volume of the combined extracts was brought to 10 mL, and filtered through a 0.2 μm PTFE filter before LC-MS analysis.
Quantitative, targeted analysis of phenolic metabolites was performed with a previously validated protocol by UPLC/QqQ-MS/MS , with few modifications. A Waters TM Acquity UPLC system (Milford, MA) consisting of a binary pump, an online vacuum degasser, an autosampler and a column compartment was used. Phenolic compounds were separated at a temperature of 40 °C on a Waters TM Acquity HSS T3 column 1.8 μm, 150 mm × 2.1 mm (Milford, MA, USA). The mobile phase was composed of eluent A (0.1 % formic acid in water) and eluent B (0.1 % formic acid in acetonitrile), with a flow rate set to 0.4 mL/min.
The multi-segment linear gradient elution conditions had the following steps: 0 min, 5 % B; from 0 to 3 min, to 20 % B; from 3 to 4.3 min, isocratic 20 % B; from 4.3 to 9 min, to 45 % B; from 9 to 11 min, to 100 % B; from 11 to 13 min, hold at 100 % B; from 13.01 to 15 min, the column was re-equilibrated to the initial conditions of 5 % B. The injection volume was 2 μL for both sample and standard solutions. Each sample was analyzed in triplicate. After each injection, the needle was rinsed with 600 μL of a weak washing solution (water/methanol, 90:10) and 200 μL of a strong washing solution (methanol/water, 90:10). Samples were kept at 6 °C during the analysis.
Mass spectrometry detection was performed on a Waters TM Xevo TQMS (Milford, MA, USA) instrument equipped with an electrospray (ESI) source. The capillary voltage was 3.5 kV in positive mode and −2.5 kV in negative mode; the source was kept at 150 °C ; the desolvation temperature was 500 °C ; cone gas flow, 50 L/h; and desolvation gas flow, 800 L/h. The MS parameters for MRM quantification are detailed in Vrhovsek et al. (2012). Quantification was done using Waters TM MassLynx 4.1 and TargetLynx software against a calibration curve made with pure standard.

GC-MS derivatization, run conditions and data processing
Tissue powder of 30 mg was transferred to a 2 mL tube, and metabolites were extracted in a pre-chilled methanol/chloroform/water extraction solution (2.5/1/1 v/v/v). From the upper water/ methanol phase, 100 μL of the extract was dried in a vacuum concentrator (Eppendorf Concentrator Plus) for derivatization (Lisec et al., 2006). The dried GC/MS samples were derivatized as described in Hochberg et al. (2013) with a similar retention time standard mix. The sample set also included a reference quality control of authentic metabolite standards (1 mg ml −1 , each) (Table S2). GC-MS run conditions were maintained exactly as described previously (Degu et al., 2014). 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 more accurate comparison of highly abundant metabolites (e.g., sugars). The GC-MS system and the parameters of the machine were set exactly as described in (Degu et al., 2014). 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 (Kopka et al., 2004) and normalized by the internal standard ribitol.

Data analysis and statistical treatment
Principal component analysis (PCA) was applied to the entire data set after mean-centring and unit variance scaling using R (version 4.0.2) (R Core Team, 2020). Pairwise correlations to all identified metabolites for the three withering stages were performed separately using Pearson's correlation algorithm. The corresponding p-values were computed using the "cor.test" function of R as previously described (Toubiana et al., 2012). Analysis of variance (ANOVA) was employed to analyse the effect of withering on metabolite accumulation among the three stages using R (version 4.0.2) (R Core Team, 2020).
Correlation-based network analysis was performed using a previously published method (Batushansky et al., 2016). First, data from both platforms (GC-and LC-MS) was united, inspected on the missing values and divided into three groups based on PCA clustering: early (T1-T7), middle (T8-T14) and late (T15-T21). The first (T0) and the last (T22) time-points were excluded from the network analysis to equilibrate dimensions of the three groups, and avoid artificial differences between them during the comparative analysis. Next, Spearman's correlation matrices for each group separately were computed using "psych" package (Revelle, 2020) for R (Team, 2013) as described in Batushansky et al. (Batushansky et al., 2016). Only correlation coefficients that passed the selected threshold (r > |±0.5|, p < 0.05) were transformed to network format for visualization using Cytoscape (Shannon et al., 2003). Graph theory-based network properties were calculated using the "igraph" package for R (Csardi and Nepusz, 2006).

Metabolite profiling during post-harvest withering
GC-MS and LC-MS analysis identified a number of metabolites being altered in content during the postwithering process of Corvina berries. Metabolite profiling of both GC-and LC-MS based PCA plot clearly discriminated the long-term withering dataset into three distinct groups (Figure 1). Based on the PCA trajectories, the temporal pattern of metabolite changes during the withering process can be categorized into three phases; early (day 0-28), middle (day 35-66) and late (day 69-108) stage. Berry lost weight up to 16 %, 18-29 % and 30-36 % during the early, middle and late stages of withering, respectively (Table S3). The first two principal components account for 44.6 % of the total variance (Table S3). The first principal component separated the early stage from the remaining two stages (middle and late), and the clustering was mainly due to sugars (Table  S4). The latter two withering stages (middle and late) had similar trajectories along PC1 while they resolved on PC2, mainly due to differences in organic acids and hydroxycinamic acids (Table S4).
To clearly show the post-withering metabolite alteration, a fold-change of metabolites at the different drying stages relative to the initial time point (T i ) (harvest T 0 ) is displayed in Figures 2A, 2B and 3. Identified metabolites are grouped into sub-classes representing amino acids, organic acids, sugar acids, sugars, stilbenes, flavanols, flavonols, sugar alcohols and hydroxycinnamic acids. The levels to which the different metabolite classes accumulated during the withering process vary between these three mentioned phases. The largest magnitude of changes was observed within the mid and late stage of the dehydration process when the berry underwent 18-29 % and 30-36 % weight loss, respectively (Figures 2A,B and 3).

Changes in amino acids and their derivatives during the withering process
Amino acids showed a mild change during the early stage of withering, with the exception of proline which showed an early increase up to 2.7-fold-change at 10 % weight loss (Figure 2A). During the mid (18-29 % weight loss) and late withering period (30-36 % weight loss), amino acids showed different patterns of changes. A significant difference in the accumulation of all identified amino acid occurred among the three withering stages (Table S5-S7). Pyroglutamate had mild change during early withering, accumulated moderately up to 1.8-fold change during mid-stage at 25 % weight loss, and showed mild change during the late withering phase. Parallel to the progressive water loss, proline level massively accumulated (up to 24-fold change at 31 % weight loss) a process that could be ascribed to stress-like conditions developing in the berry. Proline is one of the main water stressrelated molecules in grapevine operating in the leaves (Hochberg et al., 2013;Ju et al., 2018), berries (Savoi et al., 2017) and its progressive accumulation during the withering process is likely related to osmotic protection (Hong et al., 2000;Bandurska, 2000), stabilizing membranes and preventing oxidative burst activities (Hayat et al., 2012). Tholkappian et al. (2001) pointed out that water stress-associated conditions inhibit the utilization of amino acid in the downstream biosynthesis pathways and their build-up into protein. Due to such under-utilization, other amino acids such as serine, ethanolamine and leucine accumulated progressively during the mid and late stages of the drying process.
Plants synthesized ethanolamine directly from serine by serine decarboxylase enzyme in one reaction step and it plays an important role in the formation of cell membranes (Liu et al., 2018). These phospholipid classes are essential membrane components in most biological organisms (Vance, 2015). The dehydration process accelerates the rate of respiration (Kays, 1991;Costantini et al., 2006) and induces stress-like condition which facilitates the disintegration of biological membranes and could explain the accumulation of ethanolamine. Putrescine, the naturally occurring polyamine, is a well-known stress-associated amino alcohol, which notably increased (> 8 fold) during the mid and late stage of the withering. Similarly, environmental stress such as elevated temperature (Sweetman et al., 2014) and water deficit (Savoi et al., 2017) are also known to upregulate putrescine biosynthesis in grape berries. In other plants, Flores and Galston (1982) also pointed out a sharp increase in putrescine biosynthesis under osmotic stress. This result supports the previous prediction suggesting the role of polyamine in the re-generation of H 2 O 2 for ROS-mediated drought stress responses (Alcázar et al., 2011).
The other amino acids such as alanine, isoleucine, valine, beta-alanine, pyroglutamate, glutamate and arginine showed a moderate gradual increase. The general increase in amino-acid levels during the dehydration period could be due to the degradation of existing large protein molecules facilitated by proteases and peptidase released from damaged cells (Mandelstam, 1963;Yao et al., 2006). Indeed, expression of enzymes encoding various aspartate-, subtilisin-, cysteine-, and serine proteases/peptidases, among others were not only highly expressed but also extensively modulated in several berry transcriptomes (e.g., Corvina, Merlot, Sangiovese) subjected to postharvest dehydration (Zenoni et al., 2016). This observation supports the hypothesis that as the protein degraded, amino acids will be available for the biosynthesis of new protein necessary for survival under stress conditions (Willetts, 1967). Alternatively, the cell might use the protein degradation process as a means to eliminate abnormal proteins formed as a result of stress. The increase in branched-chain amino acids (isoleucine, valine) presumably contributes to the peculiar aroma of Amarone wine as previously suggested in melon fruit (Gonda et al., 2010).

Changes in organic acids and their derivatives during the withering process
Organic acids play an important role in the organoleptic properties of the final wine (Bauer and Dicks, 2004). Organic acids also influence wine stability, colour, and pH (Vilela, 2019). Similar to the amino acids, all the identified organic acids and their derivatives exhibited a significant difference in accumulation among the three stages (Table S5 -S7). During the early stage of withering (up to 16 % berry weight loss), levels of maleic acid, malonic acid, malic acid, erythronic acid, citric acid, tartaric acid, and fumaric acid showed only mild changes in the content; all exhibited moderate accumulation (up to 3.2 fold change at 29 % berry weight loss) during mid-stage of withering; during the late stage of withering maleic acid, malonic acid and fumaric acid showed progressive decline while malic acid, erythronic acid, citric acid, tartaric acid maintained moderate levels ( Figure 2B).
Similarly, cinnamic acid derivatives such as cinnamate, caffeate and synapate also showed a transit peak during the dehydration process. Hydroxylated benzoate also showed a transit peak followed by a significant reduction towards the end of the dehydration period. Organic acid biosynthesis is central to carbon metabolism producing precursors in the biosynthetic route of diverse compounds such as amino acids, vitamins, polyphenols and terpenic aroma volatiles. Therefore, the observed mild changes during withering probably reflect a reduction in the utilization of metabolites in the aforementioned downstream biosynthetic pathways any further.

Changes in sugars and their derivatives during the withering process
As expected, progressive water loss caused a significant concomitant sugar accumulation, particularly sugar phosphates, fructose 1,6-bisphosphate (F1,6 BP) and glucose 1,6-bisphosphate (G1,6-BP) which increased up to 9-fold changes towards the end of withering (36 % berry weight loss) while sucrose increased moderately during the mid-stage (29 % berry weight loss) and up to 4.5-fold-change during the late stage of withering (30 % berry weight loss) ( Figure 2B). The general increase in sugar levels during withering is likely associated with the dehydration process though the accumulation of sugar in the berry is regulated by complex mechanisms and poorly understood (Delrot et al., 2000). The general increase in sugar contents in our study ( Figure 2B) is in accordance with Reščič et al. (2016), suggesting the changes in soluble sugar accumulation is the foremost change during the post-harvest withering process. Unlike pre-harvest accumulation of soluble carbohydrates in berries, which are imported into the fruit as photoassimilate, the progressive increase in sugars during the post-harvest withering period mostly associated with the loss of water during withering and anaerobic metabolism as previously suggested by Costantini et al. (2006). Thus, the withering process favoured a sugar concentration, which may lead to fine alcoholic Amarone wine with a characteristic flavour.

Changes in stilbene and other healthpromoting compounds during the withering process
Another general effect of the post-harvest withering process is the significant accumulation of stilbene compounds including cis-piceide, trans-piceide and trans-resveratrol (Figures 2A and 3). Corvina grapes are rich in resveratrol, and a strong increase in these compounds in both the grapes and the wines can be elicited by the UV-C radiation during the withering (Paronetto and Mattivi, 1999). The modulation of this pathway has emerged as one of the key effects also in our trial. Particularly, trans-resveratrol had a massive induction (> 100-fold change at 36 % berry weight loss) from the initial time point (Figure 2A), suggesting improved nutritional status of raisin grape during the period of withering. Similarly, several reports (e.g., Versari et al. 2001;Mencarelli et al. 2010;Zenoni et al. 2016) also reported that post-harvest berry withering induced the accumulation of stilbene compounds. It is well established that multiple factors including grape cultivar, genotype, location, environmental conditions attribute for the variation of grape berry stilbene compound accumulation. Nevertheless, extensive transcriptome analysis of post-harvest withering grapes (Zamboni et al., 2008;Zenoni et al., 2016) indicated induction of a large suite (> 40) of stilbene synthase genes (Vannozzi et al., 2012) as well as their transcriptional regulators such as MYB14/15 (Höll et al., 2013) and WRKY03/24 (Vannozzi et al., 2018). In addition, the withering process massively induced the fermentation of alcoholic product, triazole, 4-methyl-5hydroxyethyl, over time. This result likely indicates a progressive decline in oxygen level in the drying berry and leads to the accumulation of anaerobic fermentation products. Triazole and its derivatives possess health-promoting properties including antimicrobial, antihypertensive, antiviral, antioxidant, antitumor and antiinflammatory (Rauf and Farshori, 2011).

Changes in flavonoids during the withering process
Unlike the other metabolite classes, most of the flavanol and flavonol metabolite classes did not show a significant difference among the three withering stages (Table S5-S7). A consistent decrease in the levels of metabolites in flavanols class including catechin, gallocatechin, epicatechin, epigallocatechin (5TMS) and epigallocatechin (6TMS) procyanidin B1 and B2 is one of the prominent change of the post-harvest withering (Figures 2A and 3). This finding is in agreement with a previous report showing a progressive decline in flavanols during Corvina berry withering (Versari et al., 2001). An increase in polyphenol oxidase during the postharvest withering process (Bonghi et al., 2012;Rolle et al., 2009) might be associated with the decline in flavanol compounds. Metabolites in the flavonol group including kaempferol glucoside, quercetin glucoside, dihydrokaempferol, quercetin galactoside, kaempferol glucoside, syringetin galactoside, isorhamnetin galactoside and rutin exhibited a general increase during the withering process whereas myricetin exhibited a consistent declining trend. Previous sensory analysis suggested that flavonols were among the compounds responsible for the velvety astringency of Amarone wines, having recognition threshold concentrations in the range 0.2 to 4.8 μmol/L (Hufnagel and Hofmann, 2008a), and the velvety character is a distinctive sub-quality of the taste of Corvina wines .
In the field time-series experiment, Shojaie et al. (2016) reported that the amount of flavonol was significantly higher under severe drought stress than its amount under control conditions or mild drought stress, suggesting that flavonols are potential mitigators for drought stress (Shojaie et al., 2016). Reščič et al. (2016) reported that the increase in phenolics due to dehydration of grape berries resulted in higher levels of these compounds in the final wines compared to the control. The post-harvest withering process may advance the phenolic ripening, leading to higher levels of phenolic compounds in withered berries. Nevertheless, compared to the pre-harvest massive flavonoid accumulation where the upstream genes regulate the pathway in a coordinated pattern (Chen et al., 2017), the general moderate increase in our study could be ascribed to the concentration effect caused by the withering of berries. Each data point represents the means of three samples obtained from twenty berries. Values in each box represent fold change relative to day 0 or when the first detection level of the metabolite is recorded. «-» represents the metabolites that exist below the detection level.

Metabolite network and correlation analysis
To investigate systematic changes occurring during berry withering a correlation-based network analysis was performed on three sets of data profiles from the three phases of withering. It is well known that dramatic physiological, metabolic and textural changes characterize berry ripening following veraison. The classical analysis of metabolite profile demonstrated distinct patterns of changes in different compound classes during the long-term berry withering. Correlation-based network analysis revealed that extensive topological differences exist among the three networks of withering stages. Nodes correspond to the different metabolite classes; node colours correspond to compound classes as detailed in the figure legend. Edges between nodes are correlation coefficients that passed the selected threshold (r > |±0.5|, p < 0.05), where pink and dark edges correspond to negative and positive correlations, respectively. As the withering process advances from early (up to 16 % berry weight loss) to middle phase (18-29 % berry weight loss), it increased the number of edges from 159 to 267, and network density from 0.09 to 0.16 (Table 1). This reflects that the berry underwent minimal metabolic rearrangements during the early stage may be due to mild water loss (Supplementary Table S1). Further water loss towards the late withering stage (30-36 % berry weight loss) markedly dropped the number of edges from 267 to 143 and network density, from 0.16 to 0.08. During the middle stage of withering, however, the berry underwent a massive coordinated metabolic rearrangement perhaps due to larger water loss leads to protein degradation and activation of stress-related biosynthetic pathways. During the late stage of withering, the decrease in network connectivity and network density might be associated with a reduction in the aforementioned biological processes. Diameter (the longest/ shortest pathways between any two nodes in the graph) of the metabolic networks at early and late phases were remarkably higher compared to the middle (9 vs. 6, respectively) supporting weak interaction during early and late stages of withering. Interestingly, early and late stages networks demonstrated remarkable similarity in main properties; including the proportion of positive and negative integrations between metabolites (~60:40), while the middle stage network was characterized by the overwhelmingly high number of positive interactions (88 %).
The current networks analysis also revealed that strong connectivity between compounds were obtained in the networks of primary biosynthetic pathways. Surprisingly, specialized metabolites like flavonoids and phenylpropanoids are more connected with primary metabolites than with each other. Several studies have shown that environmental stress significantly influences the network connectivity of metabolites in various grapevine tissues. For example, water deficit stress was previously shown to increase the network connectivity of primary and secondary metabolites in grape berries (Savoi et al., 2017) and leaves (Hochberg et al., 2013). Nonetheless, the phenomenon observed here distinguishes the postharvest dehydration metabolism from the general pre-harvest berry ripening process (Degu et al., 2014) and deficit irrigation in leaves (Hochberg et al., 2013) that mainly links primary metabolism from the downstream specialized metabolites (Hochberg et al., 2013).
To gain insight into the relationship between different metabolite classes during the berry withering process both the LC-and GC-MS dataset were used in the correlation analysis (Table S8A -9C). Furthermore, correlation coefficient heat map of metabolites for the three withering stages were generated using the R corrplot package ( Figure S1 -S3). There were highly significant positive (Table S8A -S9C) correlations between ethanolamine with sugars like fructose (r = 0.9), glucose MP (r = 0.85), glucose BP (r = 0.92) and sucrose (r = 0.89) during the late stage of withering. Similar correlation values were observed during the mid-stage of withering. The correlation of tartrate with sugars also exhibited similar correlation trends during the mid and late stages of withering. These metabolites belong to different metabolic classes but showed a progressive increase during the withering process (Figure 2A, 2B) suggesting a dehydration-driven concentration of metabolites in the berries.

CONCLUSIONS
In conclusion, Corvina grape berry undergoes a pronounced shift in metabolism in response to the prolonged withering period with the accumulation of stress-related amino acids, sugars, organic acids, flavonols, stilbene compounds and ethanolamine while the flavanols declined over time. These changes define the body and flavour of the final wine. The general increase in sugars and organic acids during the withering process perhaps is likely due to the concentration effect through the decrease in berry weight. Alternatively, stress markers like resveratrol and proline exhibited a massive increase more than the concentration effect and perhaps associated with the stimulation of different genes involved in stress protection mechanisms due to water loss. The metabolite data presented here reveal the timing of important switches in both central and specialized metabolites in the withering process, which can potentially be related to wine flavour composition and production of wines with distinct styles.