Early Leaf Removal Increases Berry and Wine Phenolics in Cabernet Sauvignon Grown in Eastern Serbia

Cluster zone leaf removal is a well-established viticulture practice for improving cluster microclimate and wine quality in cooler climates, while its efficacy in warmer conditions is less is known. Here we compared the effect of early (ELR, after fruit set; diameter of berries 3–5 mm) and late (LLR, beginning of veraison) leaf removal on berry composition and wine phenolic profile of grapevine (Vitis vinifera L.) variety Cabernet Sauvignon grown in a temperate, warmer region of Eastern Serbia. Compared to the control (no leaf removal), both leaf removal treatments increased the sugar content in fresh juice and alcohol concentration in wine. Over three consecutive years (2011–2013) markedly different in temperature and rainfall, ELR was clearly most effective in decreasing weights of cluster and of one berry, and in increasing of skin share in a berry. The content of total phenols, tartaric acid esters, anthocyanins, and flavanols in berry skin and wine was the highest in ELR treatment. ELR prominently modified the phenolic profile: Increasing flavanols, myricetin and quercetine in skin and wine, and anthocyanins, peonidin-3-glucoside in skin and delphinidin-3-glucoside in wine. This work demonstrated that early leaf removal positively influenced the chemical composition of berries and wine of Cabernet Sauvignon and might be recommended for practice in the temperate


Introduction
Cluster zone leaf removal (also termed as defoliation) is one of the most effective canopy management practices for altering microclimatic conditions surrounding the grape cluster (e.g., better aeration and sun light exposure of clusters, decreased humidity) and has shown to improve berry and wine quality [1,2]. The canopy leaf area greatly influences the changes in the properties of clusters and berries, such as the contents of sugars, total acids, and colored and aromatic compounds. Accordingly, partial removal of leaves in the cluster zone has been recognized as one of the most important activities in the vineyard. The modification of canopy microclimate and light exposure of clusters by leaf removal in the cluster zone showed the influence on the chemical and organoleptic properties of grapes and wines [3][4][5][6][7][8]. This effect mainly depends on the degree of leaf biomass removal, and on the timing [8,9]. In most vine regions worldwide, leaf removal is usually performed in the period after flowering to the beginning of veraison. Leaf removal from the cluster zone is usually practiced in cooler to moderate climates to increase temperature and light quality in the cluster zone [8,10]. However, in the warm climates a removal of leaves in cluster zone, in particular before or during veraison, may cause sunburns and overheating of berries, which lowers the wine quality [11]. In a cooler climate, early leaf removal usually decreases grape yield due to decreased berry size, but in parallel increases components of wine quality [12]. However, leaf removal in warmer climates has less been studied. In such climates, early leaf removal and subsequently increased cluster exposure to sun light, although less affecting grape yield, increased berry sugar and phenolics, and maintained must acidity [8], but occasionally led to berry discoloration [13].
Phenolic compounds are very important for the quality of red wines, affecting both organoleptic characteristics and benefits for human health [14]. The content and composition of the phenolic compounds depends on the variety, ecological conditions, and vineyard management. Hunter et al. [5] demonstrated that the concentration of anthocyanins increased with partial defoliation at veraison, while the wine quality was considerably improved by this measure regardless of the timing. Late defoliation (post-veraison and pre-harvest) on Sangiovese increased the total flavanols in berry skin; however, it did affect the total anthocyanins [15]. Moreover, in the same study [15], no significant differences between late defoliation and control (no defoliation) treatments were found for total polyphenols, alcohol, titratable acidity, and pH in wine.
Most of phenolic compounds, especially anthocyanins, are concentrated in the berry skin. Early leaf removal (after flowering to the berry set stages), for instance, reduced the size and number of berries and increased the skin/flesh ratio, resulting in an increased content of phenolic compounds in wine [16]. It has further been shown that the smaller size of berries was a result of the reduction of the assimilation surface during the flowering and fruit set stages [17]. Leaf removal also resulted in the higher sugar content in must, increased phenolic content in berry skin and more of stable anthocyanins in wine [18,19].
Despite the extensive research in the past decade, the effect of the timing of cluster zone leaf removal for the specific region and rather continental climatic conditions of Eastern Serbia characterized by warm summertime (July-September), is not well understood. It has been hypothesized that the cluster zone leaf removal will significantly influence the quality of both grapes and wine of Cabernet Sauvignon grown in the specific climate conditions of Eastern Serbia. Therefore, the aim of this study is to evaluate the effects of cluster zone leaf removal timing (early and late, compared to no leaf removal) on the quality of Cabernet Sauvignon grapes and wine with special emphasis on the phenolic profile.

Vineyard Description
This research was conducted over three consecutive seasons (2011-2013) on the Vitis vinifera L. cv. Cabernet Sauvignon. The vineyard was planted in 2001 on rootstock 5BB, with the 3333 plants per hectare (3 × 1 m). The location of vineyard is in Negotinska Krajina wine region, Serbia, South-eastern Europe (44 • 04 05 N and 22 • 33 26 W; altitude of 140 m). Bud break for Cabernet Sauvignon in this area usually occurs by mid April, flowering by mid June; veraison by early August, with harvest by late September, and leaf fall in November. The climate of the experimental area is temperate mesothermal with hot and humid summers (Köppen class Cfa), characterized by both sub-continental and sub-Mediterranean influences [20]. With average annual precipitation about 570 mm (of which 54% is delivered during vegetation period), Lang's rain factor about 45 and de Martonne's aridity index about 24, the research locality is under semi-arid conditions. Mean annual temperature is 11.8 • C and average daily temperature from April to September of 19.2 • C [21].
The vineyard was on a slight slope with the row orientation approximately north-west to south-east. The soil was a Vertisol; pH (in 1 M KCl) 6.07, 1.87% of soil organic matter, 9.7 mg kg −1 soil of ammonium lactate-extractable P and 80.5 mg kg −1 of exchangeable K; no fertilization/irrigation during the entire experiment were applied. Vine training system was a modified two-arm cordon with a trunk height of 1.4 m. The number of buds retained by pruning was 24, distributed on two canes with ten buds and two replacement spurs with two buds each. The locally common soil management and fungicides were applied to all the treatments.

Experimental Design
Leaf removal from the cluster zone was done manually in two period of times: Early leaf removal (ELR) after berry set (diameter of berries 3-5 mm; BBCH 73 growth stage, Lorenz et al., [22] and late leaf removal (LLR) at the beginning of veraison (BBCH 81 growth stage), with the same intensity (removing 5 fully developed basal leaves per shoot), and no leaf removal in the control. The lateral shoots were retained. The experiment was carried out in a random block design with three blocks and three treatments with 10 plants per treatment in each block.

Grape Harvest, Analyses of Grape Juice and Estimation of Vegetative Growth
At commercial harvest time all clusters per vine were harvested (in the stage of "technological maturity" reached at different dates: 25 September 2011, 15 September 2012 and 28 September 2013). Ten vines were harvested in each replicate (block) per treatment. The clusters were weighted to determine the yield and mean cluster weight. The berries were counted from 30 randomly taken clusters per each replicate. In addition, 50 berries were randomly taken from each sample for weighting and then crushed by fingers in order to remove the skins. The skins were rinsed with tap water to completely separate the pulp and, after which they were blotted with filter paper and weighed. The relative skin share (%) was calculated following equation: skin weight × 100/berry weight. The composite grape samples from 10 vines (each containing 12 kg) were divided into two parts; one part (about 2 kg) was frozen (−20 • C) for further chemical analyses, and the other (about 10 kg) was crushed by hand and then subjected to microvinification.
At harvest, total sugar content in the fresh grape juice was determined by portable refractometer (MR200ATC Milwaukee Meters, Neutron Pl Rowville, VIC, Australia) and the titratable acidity (total acid concentration) was determined by titration with 0.1 N NaOH to a pH 8.2 end point.
The vine vegetative growth was estimated by measuring pruning weight per vine during the winter period following the growing season.

Microvinification
The must was treated with potassium metabisulfite (0.1 g kg −1 ) and inoculated with the commercial wine yeast (BDX, Lallemand Oenology, Castel dAzzano, Italy) (0.25 g kg −1 ). After 10 days of maceration and primary fermentation in 15 L plastic pots, the fermenting juice was removed and the skins are pressed by hand to extract the remaining juice. The secondary fermentation was continued in 5 L glass demijohn for the next 10 days and then treated with 0.05 g of potassium metabisulfite. After one month wine was transferred to another 5 L glass demijohn for and after 3 months wine was bottled (0.75 L) and aged for 6 months prior to analyses.

Chemical Analyses of Berry Skin and Wine
Thirty berries were randomly taken from each frozen sample (2 kg collected from 10 vines), the skin was removed by forceps, homogenized in a mortar using liquid nitrogen, and 1 g of the skin homogenate was taken and mixed with 10 mL 50% methanol in the plastic tubes. Tubes were closed, shaken and left for 60 min at room temperature, with shaking every 10 min. The supernatants were subjected to further analysis of phenolic compounds.
Total phenols, total tartaric acid esters and total flavanols were determined following the method of Mazza et al. [23]. The samples (skin extract as described above or 6-month bottle aging wine) were mixed with acidified ethanol (2% HCl in 95% ethanol) in the ratio of 1:10 (v/v), incubated for 15 min at room temperature, filtered through a membrane filter (0.45 µm pore size) and measured spectrophotometrically (Agilent 8453 UV-VIS Spectrophotometer, Santa Clara, CA, USA) at 280 nm (total phenols; gallic acid standard), 320 nm (tartaric acid esters; caffeic acid standard) and 360 nm (total flavanols; quercetin standard). The concentrations were expressed in mg standard equivalent per g of berry skin. Total anthocyanins were determined by the modified method of Di Stefano [24]. After 15-min incubation of 0.5 mL samples (skin extract as described above or 6-month bottle aging wine) in 10 mL solution containing mixture ethanol:d-water:hydrochloric acid (70:29:1 v/v/v), the absorbance was measured spectrophotometrically (Agilent 8453 UV-VIS Spectrophotometer, Santa Clara, CA, USA) at 520 nm, and the total anthocyanins were expressed as malvidin-3-glucoside per g of berry skin.
In addition, the individual phenolic compounds (flavanols and anthocyanins) were determined in the berry skin and wine samples from the year 2012, which happened to be rather typical for the Eastern Serbia terroir (see Table 1 and Supplementary Table S7). The profile of flavanols (e.g., malvidin-3-glucoside, peonidin-3-glucoside, petunid-3-glucoside, cyanidin-3-glucoside and delphidin-3-glucoside) and anthocyanins (e.g., quercetin-glucoside, quercetin, rutin, morin, myricetin, and kaempferol) was determined by the high-performance liquid chromatography (HPLC, Agilent Technologies 1200, Santa Clara, CA, USA) with diode array detector (DAD), whereas the catechin was determined by HPLC with fluorescence detector (275-322 nm). Prior to the HPLC analysis, the extracts were further filtered through membrane filters with a pore size of 0.45 µm. Separation was carried out on the column Agilent-Eclipse XDB C-18 4.6 × 150 mm; the column was calibrated at 30 • C; the solvents were: Formic acid:water (5:95 v/v) − solvent A, and acetonitrile:formic acid:water (80:5:15 v/v) − solvent B. The linear elution gradient was used as follows: From 0 to 10 min, 0.0% B, 10 to 28 min, 10.0% B, 28 to 35 min, 25% B, 35 to 40 min, 50% B, 40 to 45 min, 80% B and for the last 10 min again 0% B. The injection volume was 5 µL and the flow of the mobile phase was 0.9 mL min −1 . For detection wavelengths, 280, 320, 360, and 520 nm for DAD were selected, and 275/322 nm (λ Ex /λ Em ) for fluorescent detection. Different phenolic compounds were identified by comparing their retention times and spectral characteristics with the data of original standard components. Quantification was performed using external calibration with the standards of phenolic compounds (HPLC-grade, ≥98%, Sigma-Aldrich, St. Louis, MO, USA). The calibration curves were linear, with the coefficient of determination R 2 = 0.99. Table 1. The selected weather parameters of the experimental years (2011-2013). Growing degree days (GDD; base 10 • C), rainfall registered during the bud burst to harvest period, annual rainfall, insolation and average daily maximum air temperature (AMDT) during the veraison to harvest period for Cabernet Sauvignon (August-September). Alcohol in the wine samples was determined using pycnometer method. The titratable acidity (total acid concentration) in wine was determined by titration with 0.1 N NaOH to a phenolphthalein end point. Volatile acidity in wines was determined by distillation using a Cazenave-Ferré followed by titration with phenolphthalein.

Statistical Analysis
The data were statistically analyzed by the 2-way ANOVA model with interaction term, after normality of residuals and homoscedasticity tests. The two treatments were leaf removal (control-no leaf removal, late leaf removal, and early leaf removal), and growth season (i.e., experimental year: 2011, 2012, and 2013). The significance of the differences between treatment means were analyzed by Tukey's Honestly significant test, using a conventional level of significance of 5% (α = 0.05). Analyses of variance were performed by the STATISTICA 6 software (StatSoft Inc., Tulsa, OK, USA).
The multivariate structure of the collected data was visualized by Principal Component Analysis (PCA) by PC-ORD 6 software (MjM Software Design, Gleneden Beach, OR, USA), using correlation cross-product matrix, Euclidean distance and no rotation. Correlation cross-product matrix was centered and standardized (relativized by standard deviate), what brought all the analyzed variables to an equal footing. The differences among the chemical profiles of skin and wine, induced by the leaf removal treatment, were analyzed by Multi-Response Permutation Procedure (MRPP; Mielke and Berry [25]). This is a nonparametric procedure for testing the hypothesis of no difference between pre-defined groups; data were relativized by standard deviate to bring all the variables to the same footing prior to analyses and be in accordance with the PCA.

Vine Growth and Grape Yield
Vine vegetative growth, grape yield and yield components were differently affected by leaf removal treatment and experimental year ( Table 2, Supplementary Table S1). Figure 1 shows the main effect of a factor with the strongest influence on each measured parameter (after ANOVA analysis, Supplementary Table S1). Overall, pruning weight ( Figure 1a (Table 1; see also Supplementary Table S7). Vegetative growth (expressed as pruning weight per plant) and number of clusters per vine were even not affected by leaf removal treatment ( Table 2, Supplementary Table S1). Compared to other years, pruning weight per vine was on average by 37% higher in the year with the highest amount of precipitation (2013-see Table 1), Figure 1a). Grape yield ranged from 1.99 kg (ELR treatment in 2012) to 5.08 kg (LLR treatment in 2013) ( Table 2). Across all the canopy management strategies evaluated (no leaf removal, early, and late leaf removal), the grape yield was in 2012 by 43% lower compared to the year 2013, and by 21% lower compared to the year 2011. However, across the three growth seasons, the average grape yield was the lowest in ELR treatment (3.34 ± 1.14 kg per vine); it decreased by 9% in comparison with the LLR treatment, and by 17% compared to the control. The number of clusters per vine ( Figure 1c) followed the same pattern as grape yield; it was the highest in 2013, and decreased by 73% in 2012. Table 2. Effects of leaf removal time and experimental year on the vegetative growth and yield components of Cabernet Sauvignon. LLR, late leaf removal; ELR, early leaf removal. Data are means of three replicates ± sd. Different letters in a column denote significant differences among treatments according to Tukey's test following a two-way ANOVA at α = 0.05.

Experimental
Year   Table S1). For the weight of clusters and weight of berries per cluster, main effect of growth season could not even be detected. These three parameters were the lowest in ELR treatments (see also Table 2). In particular, the weight of a single berry in vines subjected to early leaf removal was consistently the lowest (compared to control and LLR treatments) in each experimental year (Table 2). Across three years, these three parameters did not differ between control and LLR. The ELR treatment however reduced cluster weight by 18%, weight of berries by 23%, and weight of one berry by 20% in comparison to the averages in control and late removal treatments ( Figure 2). With an exception for the year 2011 (which was, incidentally, a very dry one; see Table 1), the results for the other two years suggest that the ELR treatment might be an efficient practice in controlling the yield. Since the yield depends on several factors, such as number of berries and weight and number of clusters, the lower ELR yield in the experiment is probably due to the lower average weight of clusters and berries compared to control and LLR treatments. Present results are in agreement with the results of various studies which reported that early leaf removal (flowering and berry set) is an effective measure in controlling the yield in different climate conditions [9,[26][27][28]. A lower grape yield recorded in ELR treatment was mainly caused by a decreased berry weigh, and thereby weight of clusters (Table 2; Figure 1), which is in accordance with the results of Moreno et al. [26] in Tempranillo grown in a temperate warm vineyard similar to ours (see Table 1 in this study and Table 1 in the study of Moreno et al.). However, Bledsoe et al. [29] reported that defoliation did not affect grape yield components and composition (with exception of an increased TSS in the early leaf removal treatment) of Sauvignon Blanc grown in the climate conditions of Napa (California), with more robust training system and lower vine density than in the present study.

General Quality Parameters of Gape Juice and Wine
Experimental year exhibited the strongest effect on both total sugar content ( Figure 2a) and titratable acidity (Figure 2b) in fresh grape juice (see Supplementary Table S2). The tested canopy management practices, furthermore, did not show any statistically significant effect on the total acidity of grape juice (Supplementary Table S2). On average, across three leaf removal treatments, sugar content (expressed as TSS) in 2012 was by 12% and 15% higher than in 2011 and 2013, respectively. This increase in sugar content in the year 2012 was most prominent in ELR treatment (Figure 2a), which is in agreement with findings of Moreno et al. [26]. Likewise, total juice acidity across all the leaf removal treatments in 2012 (on average 8.2 g L −1 ) was by 9% and by 18% higher than in 2011 and 2013, respectively. Furthermore, across three different growth seasons, TSS in juice was the lowest in the control, and did not significantly differ among the two leaf removal timings. Incidentally, the year 2012 happened to be very favorable for grapevine, with the highest temperatures and insolation, and no drought stress (Table 1). Also, in Tempranillo grown in La Rioja (Spain), basal leaf removal at fruit setting stage led to more ripened grape berries in terms of higher TSS and decreased acidity [30]. However, very recent study on Cabernet Sauvignon grown in the North Serbia terroir showed that early leaf removal decreased total acidity of grape juice, but did not affect TSS [31].
Consistently with the sugars and acids in fresh berry juice (Figure 2), the experimental year exhibited the strongest influence also on quality parameters in wine (Figure 3; Supplementary Table S3). Nevertheless, it has been shown that wine quality might not always be simply a function of TSS and berry acidity [8]. Alcohol concentration in wine was the highest in the year 2012, and in particular in the in ELR treatment (15.7%, Figure 3a). In an apparently most favorable year for wine (2012), alcohol concentration in wine was by 24% higher than in other two experimental years. On average, across three years of the experiment, the concentration of alcohol in wine was the lowest in the control treatments (no leaf removal, Figure 3a). The years 2011 and 2013 conspicuously differed in the rainfall (2011 was rather dry, and 2013 rather humid; Table 1, Supplementary Table S7). and what coincided with a prominent difference in pruning weight, grape yield and number of clusters (Figure 1a-c). Yet no differences between these non-optimal years was observed for the TSS in fresh juice (Figure 2a) nor for the alcohol in wine (Figure 3a). Furthermore, LLR treatment caused, on average across the experimental years, a small, but statistically significant increase of titratable acidity in wine by 4% compared to ELR, and by 7% compared to the control treatments (Figure 3b). In contrast to our findings, wine parameters such as alcohol content and titratable acidity were unaffected by defoliation of Grenache vines grown in La Rioja terroir, although early leaf removal positively affected wine color and sensory parameters [27]. On the other hand, the experimental year did not markedly change the pattern of how the volatile acidity in wine were affected by leaf removal (Figure 3c). Volatile acids in wine were, across the three experimental years, by up to 30% lower in control, then in the treatments involving leaf removal (Figure 3c). These results are in agreement with the findings that wine from shaded clusters of Shiraz showed significantly lower alcohol content, total and titratable acidity [32].

Polyphenolic Composition of Grape Berry Skin and Wine
Interestingly, across three consecutive experimental years markedly differing in precipitation and temperature (Table 1), ELR treatment had a clear, consistently favorable effect on polyphenolic composition of berry and wine. For instance, leaf removal treatment was the strongest modifier of a proportion of skin in a fresh berry (Supplementary Table S4). Share of skin was in fact the only parameter measured in this study not affected by leaf removal x experimental year interaction. In this study, ELR treatment resulted in the highest skin proportion in grape berry in all the experimental years, while there was no difference between LLR and control treatments (Figure 4). Across the three experimental years, average skin share of 6.71% in ELR treatment was 1.14-fold and 1.20-fold higher than in LLR and control treatments, respectively. Early leaf removal (before or after flowering) may also induce a so-called "photosynthetic shock" i.e., reduced overall photosynthetic activity of grapevines, which in turn may cause a reduced influx of assimilates from leaves to blossoms or newly formed clusters and consequently formation of the smaller berries with an increased skin-to-flesh ratio [30,31]. Excessive sunlight exposure of clusters due to an intensive leaf removal, together with high cluster temperature also resulted in smaller berries comparing with those in the shade [13,32]. Increased percentage of the skin (Figure 4) in smaller berries (Figure 1f) was associated with the increased content of phenolic compounds in berry skin (Table 3), especially anthocyanins whose synthesis occurs in the grape berry skin. Reduced mesocarp (flesh) and increased proportion of the skin is crucial for the quality of the wine, since the phenolics from the skin are extracted during the maceration process [17,31]. The highest proportion of the skin in the berry mass (ELR treatment, Figure 4) consequently reflected the highest content of the phenolic compounds in the ELR wine, extracted during the maceration process in each of the three experimental years (Table 4).   Table 3. Effects of leaf removal time and experimental year on the concentrations of total phenols, tartaric acid esters, and anthocyanins and flavanols in berry skin determined by HPLC. LLR, late leaf removal; ELR, early leaf removal. Data are means of three replicates ± sd. Different letters in a column denote significant differences among treatments according to Tukey's test following a two-way ANOVA at α = 0.05.

Experimental
Year In the present study on Cabernet Sauvignon, both leaf removal times changed the phenolic composition of the berry skin (Table 3 and Figure 5a) and wine (Table 4 and Figure 5b) in comparison to the control (no leaf removal). In berry skin, the influence of statistical interaction of leaf removal x experimental year (Supplementary Table S5) strongly modified the four analyzed parameters (Table 3). This interaction had the strongest influence on tartaric acid esters, and the weakest on anthocyanins. The simple main effect for year (across three leaf removal treatments) resulted in the highest content of skin total phenols, tartaric acid esters and flavanols in 2011 and the lowest in 2013, while anthocyanins were the highest in 2012. The simple main effect of defoliation clearly showed a superior effect of ELR on the increase of these parameters, while the differences between LLR and control were inconsistent. On the contrary, leaf removal treatment was the strongest modifier of the concentrations of total phenols and flavanols, and tartaric acid esters in wine (Supplementary Table S6). Only the concentrations of total anthocyanins both in berry skin and in wine were decisively under the strongest influence of the characteristics of the growth season (Supplementary Tables S5 and S6).
Agronomy 2021, 11, x FOR PEER REVIEW 12 of 18 effect of leaf removal treatment). Early leaf removal further caused an increase of the content of total tartaric acid esters in wine by 18% and by 38% relative to the LLR and control treatment, respectively. Finally, although the content of total anthocyanins in wine largely depended on the experimental year (Supplementary Table S6), in the representative 2012 year, it also clearly had the same pattern (ELR > LLR > Control, Table 4). Light stimulates expression and activity of a number enzymes involved in the synthesis of phenolic compounds, and shading reduces these processes [6]. Although light is essential for the production of anthocyanins, direct exposure of clusters to sun light is much more important for the production of flavanols then for the accumulation of anthocyanins; therefore, various authors proposed flavanols as some kind of markers for light exposure and microclimate of clusters [12,28,[33][34][35][36][37][38]. For instance, sun-exposed clusters had greater berry skin concentrations of total and individual (e.g., myricetin, quercetin Despite a strong modifying influence of the experimental year on the effect of leaf removal on skin concentrations of total phenols, tartaric acid esters, flavanols and in particular on anthocyanins and (Table 3, Supplementary Table S5), a general pattern of their increase (across all experimental years) is indicated along the gradient from control, over LLR, to ELR (Figure 5a). Moreover, if these four chemical compounds are used as "fingerprints" to characterize the effect of treatment, treatments involving defoliation are significantly different (by higher concentrations) from the control treatment, with exception of anthocyanins which were more influenced by year and uncorrelated with the main gradient (along the PCA Axis 1, Figure 5). This is in agreement with Moreno et al. [26] in Tempranillo grown in semiarid conditions of Western Spain. These authors reported that early leaf removal did not significantly change the concentration of total anthocyanidins, but increased the concentrations of flavanols. In the same variety (Tempranillo, two-year trial), Diago et al. [30] reported that leaf removal also might significantly increase the accumulation of anthocyanins in berries, but not in every year what is again consistent with our findings (see Supplementary Tables S5 and S6).
In wine, the concentration of total phenols, total tartaric acid esters and total flavanols very regularly increased from no leaf removal (control), over late leaf removal, and were the highest in ELR treatments, in each experimental year (Table 4, Figure 5b). In all treatments the profiles of these four compounds were statistically different (based on MRPP analysis). The total phenols content in ELR wine (average across three years 1924 ± 144 mg L −1 ) was by 15% higher than in LLR wine and about 42% higher than in the control wine. Likewise, total flavanols in ELR wine (average for three years 335 ± 62 mg L −1 ) were 1.22-fold higher than in LLR, and 1.79-fold higher compared to the control wine (main effect of leaf removal treatment). Early leaf removal further caused an increase of the content of total tartaric acid esters in wine by 18% and by 38% relative to the LLR and control treatment, respectively. Finally, although the content of total anthocyanins in wine largely depended on the experimental year (Supplementary Table S6), in the representative 2012 year, it also clearly had the same pattern (ELR > LLR > Control, Table 4).
Light stimulates expression and activity of a number enzymes involved in the synthesis of phenolic compounds, and shading reduces these processes [6]. Although light is essential for the production of anthocyanins, direct exposure of clusters to sun light is much more important for the production of flavanols then for the accumulation of anthocyanins; therefore, various authors proposed flavanols as some kind of markers for light exposure and microclimate of clusters [12,28,[33][34][35][36][37][38]. For instance, sun-exposed clusters had greater berry skin concentrations of total and individual (e.g., myricetin, quercetin and kaempferol) flavanols than in shaded clusters of Merlot [36]. These findings are in agreement with our results (Figures 5 and 6), where anthocyanins were not correlated to the statistically significant PCA axis, while flavanols strongly corresponded to the gradient induced by leaf removal along the PCA axis 1. According to various authors, the shaded clusters showed significantly lower anthocyanins concentration than those exposed to sun [7,28,32,39].
Furthermore, the detailed profile of flavanols and antocyanins in response to leaf removal treatments in berry skin (Table 5, Figure 6a) and in wine (Table 6, Figure 6b) are shown for the year 2012. This year incidentally happened to be rather typical for the Eastern Serbia terroir (see Table 1 and Supplemental Table S7). In fact, 2011 was rather dry and 2013 was cooler and more humid than the average weather conditions of this region. All the treatments (characterized by 12 phenolic compounds) were statistically significantly different (MRPP test for 2012 data, skin and wine, not shown). The skin concentrations of total flavanols measured over three-year experiment were the highest in ELR treatments (Table 3, Figure 5a), and, consistently, all the six detected flavanol compounds (except rutin) had the highest concentration in the ELR treatment of 2012 (Table 5, Figure 6a). Like with total flavanols (Table 3, Figure 5a), the differences in 7 phenolic compounds detected in skin were not consistent between LLR and control. The accumulation of kaempferol, quercetin and morin in particular was 2-fold higher in ELR than in the control in 2012 ( Table 5). The prominent accumulation of flavanols in the berry skin (catechin, quercetin, quercetin-glucoside, morine, myricetine, and kaempferol) in the ELR treatment (Table 5, Figure 6a) is in agreement with findings of Downey et al. [6] in Shiraz,Ćirković et al. [28] in Prokupac, and Spayd et al. [36] in Merlot. The skin concentrations of total anthocyanins were the highest in ELR treatments in each experimental year (Table 3, Figure 5a), but their values highly varied between years (as discussed previously). In the samples from 2012 ( Table 5, Figure 6a), five out of six detected anthocyanins had the highest concentration in the ELR treatment. Malvidin-3-glucoside was the dominant anthocyanin in the berry skin of Cabernet Sauvignon in 2012 (Table 5), which is in accordance with the findings of Ristic et al. [32] in the berry skin of Shiraz. Malvidin-3-glucoside showed the highest concentration in ELR treatment which differed significantly from the other two treatments (Table 5, Figure 6a). On the other hand, the concentrations of the three anthocyanins (peonidin-3-glucoside, petunid-3-glucoside, delphidin-3-glucoside) in berry skin were significantly lower in the LLR treatment than in ELR and control treatments. Relative to the control, the ELR treatment increased the concentration of anthocyanin compounds by up to 40% in 2012 (Table 5) and Spayd et al. [36]. and kaempferol) flavanols than in shaded clusters of Merlot [36]. These findings are in agreement with our results (Figures 5 and 6), where anthocyanins were not correlated to the statistically significant PCA axis, while flavanols strongly corresponded to the gradient induced by leaf removal along the PCA axis 1. According to various authors, the shaded clusters showed significantly lower anthocyanins concentration than those exposed to sun [7,28,32,39].  Tables 5 and 6. All the variables were correlated with the ordination scores by >35%.
Furthermore, the detailed profile of flavanols and antocyanins in response to leaf removal treatments in berry skin (Table 5, Figure 6a) and in wine (Table 6, Figure 6b) are shown for the year 2012. This year incidentally happened to be rather typical for the Eastern Serbia terroir (see Table 1 and Supplemental Table S7). In fact, 2011 was rather dry and 2013 was cooler and more humid than the average weather conditions of this region. All  Tables 5 and 6. All the variables were correlated with the ordination scores by >35%. Table 5. Effects of leaf removal timing on the concentrations (mg g −1 FW) of different phenolic compounds (flavanols and anthocyanins) in the grape berry skin determined by HPLC. LLR, late leaf removal; ELR, early leaf removal. The results from representative year (2012) are shown. Data are means of three replicates ± SD. Different letters in a row denote significant difference between parameter means after Tukey's test following one-way ANOVA (p < 0.05).
of leaf removal timing on total anthocyanins in wine. In samples from 2012, four out of five detected anthocyanins in wine were the highest in the ELR treatments ( Table 6). The strongest increase in ELR treatment was found for delphinidin-3-glucoside (27% relative to the LLR, and 83% relative to the control. The changes in phenolic profiles shown here are only illustrative and have no predicative power. For instance, while in one year we found that ELR increases delfphidine-, malvidine-, cyanindin-, and petunidin-3 glucosides (compared to ELR and control, Table 6), Ivanisevic et al. reported, on the same variety Cabernet, in Serbia, also for one year, a respective effect for delfphidin-and cyanindin-, no effect on petunidin-, and a decrease of malvidine-3 glucosides.
Overall, the general ELR phenolic profile (in the year 2012) was distinguished from the other treatments primarily by the highest concentrations of total quercetin and myricetin in both skin and wine, followed by the highest morin in the skin and the highest petunidinglucoside in the wine ( Figure 6). Interestingly, late defoliation prominently differed from the other two treatments by the lowest concentrations of rutin, peonidin-, cyanidin-, and delphinidine-glucosides in skin tissue (Figure 6a), while in the wine this treatment was set apart by the highest concentrations of rutin and peonidin-glucoside (Figure 6b).

Conclusions
Both early and late leaf removal treatments increased the content of sugars in fresh juice and alcohol in wine, compared to the control. In a three-year field experiment with pronounced interannual variability of temperature and precipitation, the early timing of leaf removal (ELR) showed a very clear effect on a decrease of cluster weight, weight of berries in a cluster and weight of a single berry, and on an increase of skin share in a berry (compared to LLR and control). Also, ELR resulted in the lowest grape yield. Early leaf removal furthermore increased the content of total phenols, tartaric acid esters, anthocyanins and flavanols in both berry skin and wine, with the exception of the total anthocyanins content in 2013, which was atypically humid. Indeed, ELR treatment showed irrespectively of very different growth seasons. In the representative year (2012), the ELR most prominently increased the content of the flavanols, myricetin, and quercetine, in both berry skin and wine, while with respect to anthocyanins the ELR caused the highest increase of peonidin-3-glucoside in skin and the highest increase of delphinidin-3-glucoside in wine.
This work clearly demonstrated the positive impact of early leaf removal on the chemical composition of berries and wine in a rather warm and arid temperature during veraison to harvest period (July-September) of continental climate in Eastern Serbia. Therefore, this measure can be recommended to increase wine quality of Cabernet Sauvignon, but has to be carefully considered together with the other factors that affect wine quality (e.g., climatic and soil conditions, training system, and winemaking technology).

Supplementary Materials:
The following are available online at https://www.mdpi.com/2073-439 5/11/2/238/s1, Table S1: Variance partitioning based on a two-way ANOVA with the interaction term for the vegetative growth and yield components, attributable to treatment factors; Table S2: Variance partitioning based on a two-way ANOVA with the interaction term for the parameters (total soluble solids (TSS)) and total acids in fresh grape juice, attributable to treatment factors; Table S3: Variance partitioning based on a two-way ANOVA with the interaction term for the content of alcohol, total acids and volatile acids in wine, attributable to treatment factors; Table S4: Variance partitioning based on a two-way ANOVA with the interaction term for the share of skin (%) in a fresh berry mass attributable to treatment factors; Table S5: Variance partitioning based on a two-way ANOVA with the interaction term for the concentration of four major groups of chemical compounds in berry skin, attributable to treatment factors; Table S6: Variance partitioning based on a two-way ANOVA with the interaction term for the concentration of four major groups of chemical compounds in wine, attributable to treatment factors; Table S7: Average temperature ( • C), precipitation (mm) and insolation (h) during the bud burst to harvest period (2011-2013).