Next Article in Journal
Real-Time Object Tracking Algorithm Based on Siamese Network
Previous Article in Journal
Design and Acceleration of Field Programmable Gate Array-Based Deep Learning for Empty-Dish Recycling Robots
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cluster Thinning Improves Aroma Complexity of White Maraština (Vitis vinifera L.) Wines Compared to Defoliation under Mediterranean Climate

Institute for Adriatic Crops and Karst Reclamation, Put Duilova 11, 21000 Split, Croatia
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(14), 7327; https://doi.org/10.3390/app12147327
Submission received: 21 June 2022 / Revised: 14 July 2022 / Accepted: 19 July 2022 / Published: 21 July 2022
(This article belongs to the Special Issue Polyphenol and Aroma Compounds in Viticulture and Enology)

Abstract

:
Defoliation and cluster thinning are useful canopy management techniques to modulate grapevine carbon distribution and microclimate. Both techniques are directed to achieve the proper balance between fruit and foliage, and to maximize production of well-ripened fruits and quality wines. We performed five canopy treatments on Maraština grapevine grown at a commercial vineyard in the Vrgorac Valley region of Croatia: three different times of basal defoliation, cluster thinning at the veraison, and an untreated control. The effects of the canopy changes on the chemical composition of grapes and wines were studied. The treatments had variable impacts on yield components and basic wine composition. Volatile aroma compounds in produced wines were analyzed using gas chromatography–mass spectrometry coupled with a mass-selective detector. The concentrations of 70 of the 96 individual volatile compounds were significantly influenced by the canopy technique used. The concentrations of 58 of these compounds were different according the timing of defoliation. Cluster thinning at an intensity of 35% produced wines with more terpenes, esters, higher alcohols, other alcohols, volatile phenolic compounds, lactones, and other compounds than other treatments. Among terpenes, cluster thinning increased terpinen-4-ol, linalool, trans-β-farnesen, and geraniol. Odor activity value analysis revealed 16 volatile compounds that contributed to the aroma of cluster-thinned wines.

1. Introduction

Techniques to manipulate the canopy by restricting vegetative and generative growth are useful to interrupt the natural cycle of grape ripening and sugar accumulation temporarily. The two techniques most commonly used to regulate the grape microclimate are defoliation and cluster thinning. Defoliation triggers a general reprogramming of the transcriptome and metabolome [1]. This leads to changes in enzyme activity and/or kinetics of biochemical reactions in berries during development [2]. Two main factors influencing the effects of defoliation and cluster thinning on vine physiology are the intensity and berry development stage, during which the intervention was carried out [3,4,5,6]. Basal defoliation modulates the cluster zone shadow–light ratio, temperature, and airflow [2,5,6]. Basal defoliation significantly disrupted carbon flux between primary and secondary metabolic pathways, delayed sugar accumulation, increased flavonoids, and decreased herbaceous aromas [1,7].
Cluster thinning one week before veraison was shown to increase free and glycosylated terpene concentrations in Sauvignon Blanc berries [4]. Cluster thinning at an intensity of 20% at veraison increased the varietal aroma compounds citronellol, linalool, and vitispirane in sparkling Ribolla Gialla wines, but only in vineyards located in a valley [8]. In Malvasia wines, a limited effect of a low crop level (10 buds per vine) on aroma compounds was reported. A low crop level affected geranyl acetate, but not C13 norisoprenoids, C6 alcohols, and higher alcohols. On the contrary, the concentration of some esters increased with a higher crop level (20 buds per vine) [9].
Maraština is an old white grape variety (Vitis vinifera L.) mostly cultivated in the coastal regions of Croatia (Dalmatia). It produces high-quality wines with fruity flavors and moderate acidity. Maraština is used for dessert-type premium wines called Prošek in Croatia. Bunch rot and powdery mildew can be a problem due to relatively compact clusters and thin berry skin. Maraština is known in Greece and Italy as Pavlos and Malvasia Bianca Lunga (Malvasia del Chianti), respectively [10]. While many studies indicate that canopy-management techniques can result in improved fruit and wine quality, too little fruit or too many leaves can turn the vine toward overcropping and eventually lead to unbalanced wines. There is no information of the impact of canopy techniques on the aroma complexity of Maraština wines. Even more, results on other white varieties grown on dry and shallow karstic soils are scarce, and only a few papers [3,8,9] reported variations in individual wine aroma classes after application of defoliation or cluster-thinning techniques. The effects of three defoliation times (preflowering, after berry set, and at veraison) and cluster thinning at veraison on the aroma profile of white wines have not been compared under the same vineyard conditions previously. The controversial results of the effect of viticulture techniques among other varieties under different climate conditions are not applicable to the improvement of the aroma potential of neutral variety Maraština under Mediterranean conditions.
Here, we fill this information gap and compared the effects of five different canopy techniques: defoliation at preflowering (PF), after berry set (BS), or at veraison (V); cluster thinning at veraison (CT); and an untreated control (C) on the agronomic qualities, cluster architecture, and the physiochemical and aroma composition of cv. Maraština wines. This will provide practical information to optimize cultivation and targeted manipulation of aroma compounds to improve sensory attributes of Maraština wines.

2. Materials and Methods

2.1. Vineyard Site and Climate Conditions

The experiment was conducted in 2019 in a commercial vineyard in the Vrgorac Valley (43°10′9.2′′ N, 17°23′44.2′′ E, 23asl) in the Dalmatian Hinterland wine region (Croatia). Vitis vinifera L. ‘Maraština’ grafted on 1103 Paulsen rootstock was planted in 2010 in fluvisol soil (undeveloped hydromorphic soil) on karst at a spacing of 0.6 m (within row) × 2.2 m (between row) in an east–west orientation. Vines were Guyot-trained to a single-curtain vertical trellis (VSP) and pruned back in winter to a single cane of 10 buds and one spur of two buds.
The climate in the Vrgorac Valley is Mediterranean, based on climate data from the nearest meteorological station with at least one month of drought in summer. The average annual temperature in 2019 was 15.2 °C, and the annual precipitation was 1389.3 mm. Weather conditions during this experiment in Vrgorac (Figure 1) were very hot and dry. Most of the precipitations during the growing season occurred in May (maximum 235.8 mm). The driest month was August (12.9 mm), which was also the month with the hottest average temperature (27.6 °C). In August, the temperature exceeded 30 °C for three consecutive days. Detailed analysis of meteorological data from June to September revealed 78 days with a maximum daily temperature above 30 °C, and the peak of 39.5 °C was on 12 August. Minimum night temperatures were above 20 °C on 54 days, which may be crucial for the balance of photosynthesis and respiration, as well as for flavor development.

2.2. Experimental Setup

A complete randomized block design was used to allow three repetitions of each treatment. Three rows (experimental blocks), each with 15 intermediate spacing (experimental units composed of 12 vines each), were selected in the central part of the vineyard. Canopy treatments were assigned randomly to each experimental units in three rows. In total, 15 vines per each of five treatments were in one row. Four canopy-manipulation techniques were applied at specific phenological development stages according to the Eichhorn and Lorenz (E-L) standard [11]. The three treatments for basal defoliation consisted of manual removal of six basal leaves from shoots at different grapevine phenological stages of the vine: preflowering (E-L 18), at berry set (bunch at right angle to shoot, berries were 2 to 4 mm diameter) (E-L 27), and at veraison (E-L 35). The cluster thinning was carried out at veraison (E-L 35) by cutting off 35% of the clusters, based on cluster counting. The thinned clusters were positioned on the upper nodules. The interval between the first two treatments was nine days, starting on 11 June, and the third treatment was applied on 6 August 2019. A control treatment without leaf or cluster removal was also carried out. Regardless of treatment, the lateral shoots were removed to the proximal node during berry set, as is common in this region. The grapes were harvested manually at the technological ripeness of the untreated control on 3 October. At harvest, yield components, number of clusters, and total mass were determined for each vine per treatment. The total yield per treatment was mixed after weighing to obtain more homogeneous samples for further analysis and winemaking. Immediately after harvest, grapes were transported to the experimental winery with ampelometry of the Institute for Adriatic Crops and Karst Reclamation.

2.3. Analysis of Physicochemical Components of Grape Juice

Ten representative clusters per each of five treatments were randomly selected and weighed, and the number of berries in each cluster was determined. Three randomly selected 100-berry samples per treatment were weighted (g) and processed into juice for rapid physiochemical analysis of basic maturity indicators [12]. Total soluble solids (°Brix) were measured with a portable density meter (RHW-25/Brix (ATC)); the total acid content was determined using titration to the point of equivalence with 0.1 mol L−1 NaOH and the bromthymol blue indicator (expressed as g L−1 tartaric acid); and the pH was measured using a pH meter (Methrom 728, Herisau, Switzerland). Single berry weight (g) was calculated by dividing the total mass of 100 berries by the number of berries.

2.4. Microvinifications

A total of 300 kg of grapes, exactly 60 kg per treatment, was mechanically destemmed, crushed, sulfited to 25 mg L−1 free sulfur dioxide (SO2), and immediately pressed. The free-run must from each treatment was divided in three repetitions and clarified in 10 L stainless steel tanks at 10 °C for 24 h. After settling, the must was decanted and inoculated with 250 mg kg−1 of the commercial Saccharomyces cerevisiae yeast strain (Siha WhiteArome, Eaton, Langenlonsheim, Germany). During rehydration, inactive yeast nutrients (Siha Speedferm, Eaton, Langenlonsheim, Germany) were added to promote formation of active yeast cells and improve fermentation performance. The fermentation temperature was monitored daily and maintained at 18 °C. After six days, 20 g L−1 of a yeast nutrient was added (Siha Profermplus, Eaton, Langenlonsheim, Germany). Once alcoholic fermentation was complete, wines were decanted, sulfited to 30 mg L−1 SO2, and stored at 15 °C. After three months, the wines were decanted a second time, the SO2 was corrected to 20 mg L−1 SO2, and wines were bottled in 750 mL bottles.

2.5. Analysis of Physiochemical Components of Wine

The physiochemical oenological parameters (alcohol, residual sugars, total extract, pH, total acidity, and volatile acidity) of all 15 wines were determined three months after bottling using the official methods of wine and must analysis [12].

2.6. Analysis of Volatile Compounds Using GC-MS

Solid-phase extraction of volatile aroma compounds was performed as previously described [13]. Briefly, cartridges containing 200 mg of LiChrolut EN sorbent (200 mg/3 mL, Merck, Darmstadt, Germany) were preconditioned using successive washing with 3 mL of dichloromethane (UHPLC gradient grade J.T. Baker, Deventer, The Netherlands), methanol (UHPLC gradient grade J.T. Baker, Deventar, The Netherlands), and a 13% aqueous ethanol solution (LiChrosolv, Merck, Darmstadt, Germany). After loading the 50 mL sample onto the column, residual sugars and other polar compounds were washed out with 3 mL of water. The column was dried in a vacuum. The analytes were recovered using elution with 1 mL of dichloromethane. For a quality control, 50 mL of water was injected into the SPE column instead of the sample. Identification and quantification of targeted aromatic compounds in the wines was performed as previously described [14,15]. Gas chromatography–mass spectrometry analysis was performed using an Agilent 6890 system equipped with an Agilent 6890 automatic liquid sampler and a ZB_WAX capillary column (60 m × 0.32 mm i.d., with 0.5 μm film thickness, Phenomenex, Torrance, CA, USA) coupled to a 5973 N mass-selective detector. The injection was performed in splitless mode with volume of 2 μL, while the injector was held at 250 °C. The temperature gradient was as follows: initial oven temperature of 40 °C for 15 min, then increasing to 250 °C in steps of 2 °C/min and held at 250 °C for 15 min. The ion source and MS transfer line were at 250 °C. The carrier gas flow was 1.0 mL/min 5.0 grade helium (Messer, Zaprešić, Croatia) in constant flow mode. The mass detector was operated in electron ionization (EI) mode at an ionization energy of 70 eV with total ion current monitoring. The mass detector ion source and quadrupole were maintained at 250 °C. Compounds were identified by comparing the retention times and mass spectra with those of commercial standards. Six-point calibration was performed for each standard as reported in [14]. Concentration ranges for alcohols, esters, and acids were 5 ug L−1 to 5000 ug L−1; and for terpenes, C13 norisoprenoids, and lactones: 1 ug L−1 to 300 ug L−1. Standard calibration curves were generated using EnhancedChemStation software (Agilent Technologies, Santa Clara, CA, USA) to quantify the compounds. A list of used standards, linear retention indices, and other identification and quantification parameters for the GC-MS are presented in [14]. The results were expressed in µg L−1 as the concentration of the compounds calculated from the peak area of each compound compared to that of the internal standard. In the eluate, 10 µL of internal standard (50 mg L−1, 2-octanol) was added.

2.7. Odor Activity Values of Volatile Aromas (OAV)

To evaluate the potential sensory impact to the wine, the odor activity value of 73 volatile aroma compounds was calculated. The OAV was calculated as the quotients of concentration of a specific compound and the corresponding odor detection threshold from the literature [14,15]. The odor detection threshold values and aroma descriptors were taken from the list of 28 references in the Supplementary Materials and are available in Table S1. When possible, the threshold values determined in wine or winelike medium were used to assure accuracy of the OAV. Volatile aroma compounds with an OAV higher than one had a direct impact on aroma and were considered as active odorants. Those with lower values were considered as contributing to wine complexity through synergistic and additive effects of compounds with similar structure or odor (synergistic compounds) [14].

2.8. Statistical Analysis

A one-way analysis of variance (ANOVA) and mean separation using Tukey’s test (with different letters representing significant differences at p ≤ 0.05) was performed to determine differences between the five different treatments. Analyses of variance were performed in SAS (SAS Institute, Inc., Cary, NC, USA). The results of all physiochemical and aroma parameters were expressed as mean ± standard deviation of three repetitions. A multivariate analysis approach was applied to the volatile metabolites using Statistica version 14 software (TIBCO Software Inc., Palo Alto, CA, USA) [16]. A total of 16 volatile compounds with odor-active compounds >1 were selected as the most likely contributors to the aroma of the wines. Principal component analysis (PCA) was performed on the 16 volatile compounds to differentiate the wines from five canopy treatments and to analyze possible relationships among them.

3. Results

3.1. Yield Components and Physicochemical Composition of Maraština Musts and Wines

Yields components and basic physiochemical composition were determined for Maraština must and wine subjected to defoliation at preflowering (PF), after berry set (BS), or at veraison (V); cluster thinning (CT); and an untreated control (C) (Table 1). Significant differences in the number of clusters per vine were determined between PF and CT. Significant differences in yield were determined in untreated C (1.7 ± 0.8) kg and BS (1.2 ± 0.5) kg, with other treatments in-between. The highest cluster weight and number of berries per cluster was in CT, followed by untreated C and V defoliation. BS caused the greatest reduction in yield per vine and cluster architecture parameters. Defoliation significantly reduced the total soluble solids in must, with a maximum reduction in BS (18.7 ± 0.1) °Brix compared to CT (19.6 ± 0.2) °Brix or the untreated C (19.5 ± 0.2) °Brix. BS significantly reduced total acidity (4.3 ± 0.1) g L−1, followed by treatments performed at veraison (CT and V). The highest pH (3.74 ± 0.01) was in treatments conducted at BS and at veraison, which were significantly higher than in PF and C. As expected, defoliation significantly reduced the alcohol content of the wines, regardless of the timing of treatment. The greatest alcohol reduction, determined as alcoholic strength by volume, of 0.8 vol % was in PF wines, while postponing the treatment reduced differences between those and untreated C or CT wines. In wine, the lowest pH was in early defoliation treatments (PF and BS), followed by C and CT.

3.2. Aroma Composition of Maraština Wines

A total of 96 volatile aroma compounds were identified and quantified in the 15 wines made using the five treatments. These compounds were grouped into 10 chemical classes. The terpenes and C13 norisoprenoids are carriers of varietal aromas (Table 2). Fermentation aromas include esters, higher alcohols, other alcohols, aldehydes, fatty acids, volatile phenols, lactones, and other compounds. Although some aroma compounds were present only in trace amounts, their suppressor or enhancer contributions to the overall wine aroma were very important. The odor activity values with corresponding odor detection thresholds and aroma descriptors of each compound analyzed are listed in the Supplementary Materials (Tables S1 and S2).

3.2.1. Varietal Volatile Composition of Maraština Wine

The CT wines had the most linalool and citronellol: (20.40 ± 0.82) µg L−1 and (19.26 ± 3.05) µg L−1, respectively. Those were the most abundant terpene compounds in Maraština wines. In addition, CT significantly increased the concentrations of terpinen-4-ol, trans-β-farnesene, and geraniol. Consequently, the greatest total sum of terpenes was in CT at 98.44 ± 8.52 µg L−1. There were no significant differences among defoliation treatments (PF, BS, and V) and C in linalool, citronellol, or geraniol concentration or in total sum of terpenes. Defoliation reduced terpene compounds; the lowest total sum of terpene compounds (72.52 ± 3.41 µg L−1) was in PF wines.
Among defoliation treatments, PF wines had a significantly lower concentration of linear sesquiterpene farnesol and 1,8-terpin, although no difference was seen between PF and CT wines in those compounds. Defoliation at berry set led to wines with a significantly higher level of farnesol (6.52 ± 0.46 µg L−1) among all treatments, while a significant increase in terpinen-4-ol (1.60 ± 0.41 µg L−1) in BS wines was seen when defoliation treatments were compared. BS wines had significantly lower concentrations of hotrienol and α-pinene, but had no difference from untreated C. Treatments at veraison (V and CT) had significantly more terpenic hydrocarbon α-pinene. Significantly more cis-linalool oxide and 6,7-dihydro-7-hydroxylinalool were in wines of the untreated C. The treatments did not differ in concentrations of individual and total C13 norisoprenoid compounds.

3.2.2. Fermentation Volatile Composition of Maraština Wine

Wines treated at veraison, CT and V had significantly higher levels of the most abundant ester, isoamyl acetate, and total sum of esters than the other treatments (Table 3). CT wines had significantly higher levels of six individual esters (isobutyl acetate, ethyl lactate, ethyl 2-methylbutanoate, ethyl 3-methylbutanoate, ethyl furoate, and diethyl succinate) than V treatments. Ten individual esters and the total sum of esters reached a minimum in BS wines, with isobutyl acetate, isoamyl acetate, 2-phenylethyl acetate, ethyl butanoate, ethyl-2-methylbutanoate, ethyl-3-methylbutanoate, ethyl-2-hydroxy-3-methylbutanoate, ethyl hexadecanoate, ethyl-3-hydroxyhexanoate, and diethyl succinate being significantly less abundant than in C wines.
The most abundant higher alcohol in Maraština wines was 2-methyl-1-butanol regardless of treatment. Concentrations of this compound ranged from 7,997.82 ± 85.59 μg L−1 in BS wines to 12,452.80 ± 466.25 μg L−1 in CT wines. Wines of early defoliation treatments (PF and BS) had significantly lower concentrations of 3-methyl-1-butanol, isoamyl alcohol, 1-octen-3-ol, and total sum of higher alcohols. CT wines had a significantly higher concentration of 1-hexanol, trans-3-hexen-1-ol, isobutanol, 1-butanol, and 1-pentanol.
The V wines contained significantly more total aldehydes, particularly 2-octenal. In contrast, BS wines had the lowest concentrations of these compounds. Wines from the two veraison treatments did not differ in their benzaldehyde or benzeneacetaldehyde concentration. Both had a significantly higher concentration of benzaldehyde and a lower concentration of benzeneacetaldehyde in comparison to other treatments. The PF wines differed significantly in the total sum of fatty acids compared to the V and CT treatments. In general, the greatest proportion of fatty acids was found in V wines. Delaying the time of defoliation increased the total amount of fatty acids in the wines, reflecting the trend of the most abundant hexanoic acid among them. These wines also had the greatest proportion of 2-methylpropionic acid. Defoliation generally led to significantly lower total volatile phenols and the most abundant among them, 4-vinylguaiacol, which increased with the later defoliation, while the highest concentrations were in CT and C wines. The concentration of γ-butyrolactone reached a maximum of 562.58 ± 27.68 μg L−1 in CT wines. A total sum of lactones followed this trend. Among other compounds, significantly higher concentrations of acetoin were seen in BS wines, contrary to other treatments that did not differ from untreated C. The lowest concentration of C13 norisoprenoid ketone 6-methyl-5-hepten-2-one was in V wines.

3.3. Odor Activity Values of Volatile Compounds in Maraština Wine

A total of 16 volatile compounds out of 96 identified in the Maraština wines were present at concentrations above their thresholds with an OAV > 1, and contributed to the aroma of the wines. Fourteen compounds had OAV > 1 in all wines, regardless of treatment: β-damascenone, isoamyl acetate, ethyl butanoate, ethyl 2-methylbutanoate, ethyl 3-methylbutanoate, ethyl hexanoate, ethyl octanoate, 1-octen-3-ol, 1-heptanol, 2-octenal, benzeneacetaldehyde, hexanoic acid, γ-decalactone, and 6-methyl-5-hepten-2-one. Compounds with the greatest OAV across all treatments were ethyl octanoate (273.61), carrier of pineapple, pear, and floral notes; and 2-octenal (105.08), carrier of herbal and nutty notes. PF wines differed from the others in having higher OAV values for β-damascenone than ethyl hexanoate. Linalool had an OAV > 1 in wines from PF, CT, and untreated C, while 4-vinylguaiacol had an OAV > 1 in CT and C.
The data for 16 volatile variables with OAV > 1 were processed using principal component analysis (PCA) to differentiate among wines from the five different canopy-management treatments according to their volatile compounds. In a projection of 16 volatile variables that defined the principal components F1 and F2, the first two principal components with eigenvalues > 1 explained 73.7% of the variance. The contribution of F1 was 42.93%, and of F2, 30.75% (Figure 2A). According the PCA, wines of different treatments were clearly separated (Figure 2B). Wine sample replicates grouped together for all five canopy-management treatments. V and CT wines were located on the positive side of PC1 and were characterized by 11 OAV values: linalool, isoamyl acetate, four ethyl esters (ethyl hexanoate, ethyl octanoate, ethyl butanoate, and ethyl-3-methylbutanoate), 2-octenal, hexanoic acid, 4-vinylguaiacol, and γ-decalactone. Nine compounds belonged to floral, tropical fruity, sweet, spicy, or nutty aroma odor classes [17], with only two from the green and mineral class, 1-octen-3-ol and hexanoic acid (Flavornet, http://www.flavornet.org/flavornet.html, accessed on 28 February 2022). Veraison samples were located in different quadrants: V1, V2, and V3 in the first; and CT1, CT2, CT3 in the second. Early defoliation treatments (BS and PF) had little or no influence on wine OAV values compared to control wines, all of which were located on the negative side of PC1. Those wines correlated with five OAVs: 6-methyl-5-hepten-2-one, ethyl-2-methylbutanoate, β-damascenone, 1-heptanol, and benzeneacetaldehyde, compounds associated with floral, vegetable, pome, and citrus fruit characteristics. Control and preflowering wines (C1, C2, C3, PF1, PF2, PF3) were located in the third quadrant of the plot, while berry set wines (BS1, BS2, BS3) were in the fourth quadrant. PCA confirmed that the canopy-management treatment affected the volatile profiles of the wine samples.

4. Discussion

The lowest yield per vine and cluster weight was found in BS vines, and not those cluster-thinned at an intensity of 35%, as expected. It is probable that early source removal of six basal leaves immediately after berry set reduced assimilate supply and pressured the vine to select a certain number of berries that would ripen properly [5]. The lack of a significant decrease in yield, despite the low number of clusters left on the vine after cluster thinning, was due to self-compensation increasing the grape cluster and berry weight [18,19]. Those larger berries could decrease the skin-to-pulp ratio and aroma compound concentrations [18].
Regulation of yield with cluster-thinning or defoliation practices depends on vintage [5,19,20], agroecological conditions of location, and variety [19,21,22]. Despite studies describing the impact of defoliation or cluster thinning on aromatic compounds in berries, there are little data concerning their effect on the physiochemical and aroma composition of wine.
Defoliation reduced berry TSS, final alcohol content, and pH, and increased the total acidity of Maraština wines, with the greatest effect with early defoliation treatments (BS and PF). Terpene and C13 norisoprenoid compounds are important carriers of varietal aroma. In this study, defoliation at intensity of six basal leaves, regardless of the timing, reduced the total concentration of terpenes in Maraština wines compared to the untreated C. Despite use of the same training system (Guyot) in this and in an experiment conducted on Malvazija Istarska [23], the reported beneficial impact of early defoliation on linalool, citronellol, nerol, and geraniol was missing. Among the three defoliation treatments, the most individual and total sum of terpene compounds was in V wines, as reported previously for Sauvignon Blanc with 50% or 100% (six leaves) intensity of basal defoliation [6]. Defoliation early in the season had a strong impact on temperature and increase in evaporative and respiration rate of berries [24]. It could be that heat stress during herbaceous stage of berries led to a delay in the onset of accumulation of those compounds and an increase in the degradation and volatilization of already-synthetized compounds [25,26]. Previous research placed the peak of terpene synthesis at two weeks after veraison, and the hypothesis that defoliation at veraison would increase certain terpenes was partly confirmed [27].
The CT wines had the highest total sum of terpenes, with the most abundant being linalool, terpinen-4-ol, trans-β-fernesene, and geraniol. Significantly, the increases in linalool and geraniol in CT wines while others remained constant may have been an indicator of a fully mature grape [28]. One of the most important norisoprenoids is β-damascenone, a carrier of sweet, fruity, floral, and honey aromatic notes. In white wines, the perception threshold of β-damascenone was 0.05 µg L−1, as in hydroalcoholic solution (10–12%, v/v, water/ethanol mixture) [29]. In our study, wines of all treatments had high concentrations of this compound. The esters that carry floral and fruity aromas in wine are yeast-derived compounds. CT wines had the highest concentrations of isoamyl acetate, isobutyl acetate, ethyl lactate, ethyl 2-methylbutanoate, ethyl 3-methylbutanoate, ethyl furoate, diethyl succinate, and total sum of esters. The positive effect of cluster thinning on the overall concentration of esters was reported previously in Ribolla Gialla wines [8]. Similar to terpenes, an increase in esters could not be attributed only to a lower crop load, as we saw the minimum of esters in the BS wines that also had a lowest yield. It was probably a consequence of the enhancement of ripening after the cluster thinning.
The Maraština wines from early defoliation treatments had fewer higher alcohols, contrary to previous reports for Istrian Malvasia [9], which could be explained by a special viticultural condition in certain vintages. V wines did not differ from C in the total sum of higher alcohols, nor in isobutanol or isoamyl alcohol, as found previously for Sauvignon Blanc and Riesling from Continental Croatia [3]. Higher alcohols are the result of sugar or deamination and decarboxylation of certain amino acids during fermentation by the Ehrlich pathway.
Only V wines contained significantly more total aldehydes, particularly 2-octenal, a compound that is considered as one of are markers for oxidative deterioration of wine [30]. Wines from the two veraison treatments had a significantly higher concentration of benzaldehyde (almond) and a lower concentration of benzeneacetaldehyde (2-phenylacetaldehyde; rose) in comparison to other treatments.
CT wines had peak values for 7 among 10 aroma compound classes, with a significant increase in the total sum of terpenes and esters. Cluster thinning improved the aromatic maturity of Maraština wines by enhancing grape ripeness. The odor activity values (OAVs) of volatile compounds were calculated to gain insight into perceptible changes in the concentrations of volatile compounds in wine caused by five different vineyard treatments. According the odor thresholds and description of each aroma compound reported in the literature (Table S1), wines of Maraština have a fruity character. Ethyl octanoate and 2-octenal were two most odor-potent compounds regardless of treatment. A predominance of β-damascenone in PF wines, compared to ethyl hexanoate in other wines, could be important due to its role as an enhancer of other compounds in the wine matrix. CT wines had a more complex aromatic profile, as all 16 compounds in those had an OAV > 1, while in wines from the defoliation treatments, linalool and 4-vinylguaiacol were <1. Previously, high OAV values for linalool and β-damascenone were reported in the aroma of greenhouse-dried Maraština grapes [31].

5. Conclusions

A total of 96 individual volatile compounds were identified and quantified in 15 wines from five different canopy manipulation treatments. Among varietal aromas, terpenes (23) and norisoprenoids (3) were measured. Among fermentation aromas, esters (18), higher alcohols (15), other alcohols (5), aldehydes (8), fatty acids (8), volatile phenols (5), lactones (7), and other compounds (4) were measured. The canopy-manipulation treatment significantly affected the aroma profile of Maraština wines. Early defoliation treatments had a dampening effect on the total concentrations of terpene, esters, higher alcohols, and aldehydes in wines compared to treatments at veraison. Cluster thinning of Maraština at an intensity of 35% at veraison led to wines higher in alcohol, terpenes, esters, higher alcohols, other alcohols, volatile phenolic compounds, lactones, and other compounds. CT treatment stimulated production of varietal aroma compounds: four individual terpenes (terpinen-4-ol, linalool, trans-β-farnesen, and geraniol). An odor activity value analysis showed that cluster thinning improved the quality of Maraština wines by increasing the concentrations of 16 compounds: linalool, β-damascenone, isoamyl acetate, ethyl butanoate, ethyl-2-methylbutanoate, ethyl-3-methylbutanoate, ethyl hexanoate, ethyl octanoate, 1-octen-3-ol, 1-heptanol, 2-octenal, benzeneacetaldehyde, hexanoic acid, 4-vinylguaiacol, γ-decalactone, and 6-methyl-5-hepten-2-one, the carriers of floral and fruity aromas. Cluster thinning led to an improved quality of Maraština wines without a negative impact on yield, while defoliation in our specific region under Mediterranean climate conditions was the least risky canopy technique. In this context, further research on different intensities and timings of cluster thinning is required to understand the effect of this management technique on improvement of the wine quality in a Mediterranean climate.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/app12147327/s1, Table S1. Odor description and odor threshold (ODT) of 73 aroma compounds detected in ‘Maraština’ wines; Table S2. Odor activity values (OAV) for volatile compounds in the wine of ‘Maraština’ from five different canopy-management treatments: preflowering (PF), after berry set (BS), and veraison (V) defoliation; cluster thinning without defoliation; 35% clusters thinned in veraison (CT); and untreated control (C). Refs. [32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58] are cited in Supplementary Materials.

Author Contributions

Conceptualization, methodology, and formal analysis, A.M.; data analysis and preparation of tables, A.M. and G.Z.; preparation of the original draft, A.M.; review and supervision of the manuscript, G.Z., I.B.-L. and K.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We acknowledge Split-Dalmatian County for support of this research through the project: “Implementation of innovative ampelotechnical interventions in the vineyard for the production of high quality wines of native varieties” and the entire staff of the Franić winery for allowing us to host the experiment in their vineyard.

Conflicts of Interest

The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Pastore, C.; Santo, S.D.; Zenoni, S.; Movahed, N.; Allegro, G.; Valentini, G.; Filippetti, I.; Tornielli, G.B. Whole plant temperature manipulation affects flavonoid metabolism and the transcriptome of grapevine berries. Front. Plant Sci. 2017, 8, 929. [Google Scholar] [CrossRef] [PubMed]
  2. Young, P.R.; Eyeghe-Bickong, H.A.; du Plessis, K.; Alexandersson, E.; Jacobson, D.A.; Coetzee, Z.; Deloire, A.; Vivier, M.A. Grapevine plasticity in response to an altered microclimate: Sauvignon Blanc modulates specific metabolites in response to increased berry exposure. Plant Physiol. 2016, 170, 1235–1254. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Kozina, B.; Karoglan, M.; Herjavec, S.; Jeromel, A.; Orlic, S. Influence of basal leaf removal on the chemical composition of Sauvignon Blanc and Riesling wines. J. Food Agric. Environ. 2008, 6, 28. [Google Scholar]
  4. Kok, D. Influences of pre- and post-veraison cluster thinning treatments on grape composition variables and monoterpene levels of Vitis vinifera L. cv. Sauvignon Blanc. J. Food Agric. Environ. 2011, 9, 22–26. [Google Scholar]
  5. Sivilotti, P.; Falchi, R.; Herrera, J.C.; Škvarč, B.; Butinar, L.; Lemut, M.S.; Bubola, M.; Sabbatini, P.; Lisjak, K.; Vanzo, A. Combined effects of early season leaf removal and climatic conditions on aroma precursors in Sauvignon Blanc grapes. J. Agric. Food Chem. 2017, 65, 8426–8434. [Google Scholar] [CrossRef]
  6. Yue, X.; Ma, X.; Tang, Y.; Wang, Y.; Wu, B.; Jiao, X.; Zhang, Z.; Ju, Y. Effect of cluster zone leaf removal on monoterpene profiles of Sauvignon Blanc grapes and wines. Food Res. Int. 2020, 131, 109028. [Google Scholar] [CrossRef]
  7. Šuklje, K.; Antalick, G.; Coetzee, Z.; Schmidtke, L.M.; Česnik, H.B.; Brandt, J.; du Toit, W.J.; Lisjak, K.; Deloire, A. Effect of leaf removal and ultraviolet radiation on the composition and sensory perception of Vitis vinifera L. cv. Sauvignon Blanc wine. Aust. J. Grape Wine Res. 2014, 20, 223–233. [Google Scholar] [CrossRef]
  8. Škrab, D.; Sivilotti, P.; Comuzzo, P.; Voce, S.; Degano, F.; Carlin, S.; Arapitsas, P.; Masuero, D.; Vrhovšek, U. Cluster thinning and vineyard site modulate the metabolomic profile of Ribolla Gialla base and sparkling wines. Metabolites 2021, 11, 331. [Google Scholar] [CrossRef]
  9. Bubola, M.; Rusjan, D.; Lukić, I. Crop level vs. leaf removal: Effects on Istrian Malvasia wine aroma and phenolic acids composition. Food Chem. 2020, 312, 126046. [Google Scholar] [CrossRef]
  10. Maletić, E.; Pejić, I.; Kontić, J.K.; Zdunić, G.; Preiner, D.; Šimon, S.; Andabaka, Ž.; Mihaljević, M.Ž.; Bubola, M.; Marković, Z.; et al. Ampelographic and genetic characterization of Croatian grapevine varieties. Vitis. J. Grapevine Res. 2015, 54, 93–98. [Google Scholar]
  11. Coombe, B.G. Growth Stages of the Grapevine: Adoption of a System for Identifying Grapevine Growth Stages. Aust. J. Grape Wine Res. 1995, 1, 104–110. [Google Scholar] [CrossRef]
  12. Office International de la Vigne et du Vin. Compendium of International Methods of Wine and Must Analysis—Volume I; Organisation Internationale de la Vigne et du Vin (OIV): Paris, France, 2015. [Google Scholar]
  13. López, R.; Aznar, M.; Cacho, J.; Ferreira, V. Determination of minor and trace volatile compounds in wine by solid-phase extraction and gas chromatography with mass spectrometric detection. J. Chromatogr. A 2002, 966, 167–177. [Google Scholar] [CrossRef]
  14. Korenika, A.M.J.; Preiner, D.; Tomaz, I.; Jeromel, A. Volatile profile characterization of Croatian commercial sparkling wines. Molecules 2020, 25, 4349. [Google Scholar] [CrossRef]
  15. Maslov, L.; Tomaz, I.; Žulj, M.M.; Jeromel, A. Aroma characterization of predicate wines from Croatia. Eur. Food Res. Technol. 2017, 243, 263–274. [Google Scholar] [CrossRef]
  16. TIBCO Software Inc. Data Science Workbench, Version 14. 2020. Available online: http://tibco.com (accessed on 2 March 2022).
  17. Flavornet. Available online: http://www.flavornet.org/flavornet.html (accessed on 28 February 2022).
  18. Condurso, C.; Cincotta, F.; Tripodi, G.; Sparacio, A.; Giglio, D.M.L.; Sparla, S.; Verzera, A. Effects of cluster thinning on wine quality of Syrah cultivar (Vitis vinifera L.). Eur. Food Res. Technol. 2016, 242, 1719–1726. [Google Scholar] [CrossRef]
  19. Palliotti, A.; Cartechini, A. Cluster thinning effects on yield and grape composition in different grapevine cultivars. Acta Hortic. 2000, 512, 111–120. [Google Scholar] [CrossRef]
  20. Keller, M.; Mills, L.J.; Wample, R.L.; Spayd, S.E. Cluster thinning effects on three deficit-irrigated Vitis vinifera cultivars. Am. J. Enol. Vitic. 2005, 56, 91–103. [Google Scholar]
  21. Poni, S.; Casalini, L.; Bernizzoni, F.; Civardi, S.; Intrieri, C. Effects of early defoliation on shoot photosynthesis, yield components, and grape composition. Am. J. Enol. Vitic. 2006, 57, 397–407. [Google Scholar]
  22. Hernandez-Orte, P.; Concejero, B.; Astrain, J.; Lacau, B.; Cacho, J.; Ferreira, V. Influence of viticulture practices on grape aroma precursors and their relation with wine aroma. J. Sci. Food Agric. 2015, 95, 688–701. [Google Scholar] [CrossRef]
  23. Bubola, M.; Lukić, I.; Radeka, S.; Sivilotti, P.; Grozić, K.; Vanzo, A.; Bavčar, D.; Lisjak, K. Enhancement of Istrian Malvasia wine aroma and hydroxycinnamate composition by hand and mechanical leaf removal. J. Sci. Food Agric. 2019, 99, 904–914. [Google Scholar] [CrossRef]
  24. Ferrari, V.; Disegna, E.; Dellacassa, E.; Coniberti, A. Influence of timing and intensity of fruit zone leaf removal and kaolin applications on bunch rot control and quality improvement of Sauvignon blanc grapes, and wines, in a temperate humid climate. Sci. Hortic. 2017, 223, 62–71. [Google Scholar] [CrossRef]
  25. Scafidi, P.; Pisciotta, A.; Patti, D.; Tamborra, P.; Di Lorenzo, R.; Barbagallo, M.G. Effect of artificial shading on the tannin accumulation and aromatic composition of the Grillo cultivar (Vitis vinifera L.). BMC Plant Biol. 2013, 13, 175. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Rienth, M.; Torregrosa, L.; Luchaire, N.; Chatbanyong, R.; Lecourieux, D.; Kelly, M.T.; Romieu, C. Day and night heat stress trigger different transcriptomic responses in green and ripening grapevine (Vitis vinifera) fruit. BMC Plant Biol. 2014, 14, 108. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Friedel, M.; Frotscher, J.; Nitsch, M.; Hofmann, M.; Bogs, J.; Stoll, M.; Dietrich, H. Light promotes expression of monoterpene and flavonol metabolic genes and enhances flavour of winegrape berries (Vitis vinifera L. cv. Riesling). Aust. J. Grape Wine Res. 2016, 22, 409–421. [Google Scholar] [CrossRef]
  28. Baron, M.; Prusova, B.; Tomaskova, L.; Kumsta, M.; Sochor, J. Terpene content of wine from the aromatic grape variety “Irsai Oliver” (Vitis vinifera L.) depends on maceration time. Open Life Sci. 2017, 12, 42–50. [Google Scholar] [CrossRef]
  29. Pineau, B.; Barbe, J.C.; van Leeuwen, C.; Dubourdieu, D. Which impact for β-damascenone on red wines aroma? J. Agric. Food Chem. 2007, 55, 4103–4108. [Google Scholar] [CrossRef]
  30. Escudero, A.; Asensio, E.; Cacho, J.; Ferreira, V. Sensory and chemical changes of young white wines stored under oxygen. An assessment of the role played by aldehydes and some other important odorants. Food Chem. 2002, 77, 325–331. [Google Scholar] [CrossRef]
  31. Budić-Leto, I.; Humar, I.; Gajdoš Kljusurić, J.; Zdunić, G.; Zlatić, E. Free and bound volatile aroma compounds of ‘Maraština’ grapes as influenced by dehydration techniques. Appl. Sci. 2020, 10, 8928. [Google Scholar] [CrossRef]
  32. Plotto, A.; Margaría, C.A.; Goodner, K.L.; Goodrich, R.; Baldwin, E.A. Odour and flavour thresholds for key aroma components in an orange juice matrix: Terpenes and aldehydes. Flavour Fragr. J. 2004, 19, 491–498. [Google Scholar] [CrossRef]
  33. Ribereau-Gayon, P.; Boidron, J.N.; Terrier, A. Aroma of Muscat grape varieties. J. Agric. Food Chem. 1975, 23, 1042–1047. [Google Scholar] [CrossRef]
  34. Fariña, L.; Villar, V.; Ares, G.; Carrau, F.; Dellacassa, E.; Boido, E. Volatile composition and aroma profile of Uruguayan Tannat wines. Food Res. Int. 2015, 69, 244–255. [Google Scholar] [CrossRef]
  35. Guth, H. Quantitation and sensory studies of character impact odorants of different white wine varieties. J. Agric. Food Chem. 1997, 45, 3027–3032. [Google Scholar] [CrossRef]
  36. Sánchez-Palomo, E.; Trujillo, M.; García Ruiz, A.; González Viñas, M.A. Aroma profile of Malbec red wines from La Mancha region: Chemical and sensory characterization. Food Res. Int. 2017, 100, 201–208. [Google Scholar] [CrossRef] [PubMed]
  37. Slaghenaufi, D.; Luzzini, G.; Samaniego Solis, J.; Forte, F.; Ugliano, M. Two sides to one story—Aroma chemical and sensory signature of Lugana and Verdicchio wines. Molecules 2021, 26, 2127. [Google Scholar] [CrossRef]
  38. Marais, J. Terpenes in the aroma of grapes: A review. S. Afr. J. Enol. Vitic. 1983, 4, 49–60. [Google Scholar]
  39. Ferreira, V.; López, R.; Cacho, J.F. Quantitative determination of the odorants of young red wines from different grape varieties. J. Sci. Food Agric. 2000, 80, 1659–1667. [Google Scholar] [CrossRef]
  40. Francis, L.; Newton, J. Determining wine aroma from compositional data. Aust. J. Grape Wine Res. 2005, 11, 114–126. [Google Scholar] [CrossRef]
  41. Pardo, E.; Rico, J.; Gil, J.V.; Orejas, M. De novo production of six key grape aroma monoterpenes by a geraniol synthase-engineered S. cerevisiae wine strain. Microb. Cell Fact. 2015, 14, 136. [Google Scholar] [CrossRef] [Green Version]
  42. Sacks, G.L.; Gates, M.J.; Ferry, F.X.; Lavin, E.H.; Kurtz, A.J.; Acree, T.E. Sensory threshold of 1,1,6-trimethyl-1,2-dihydronaphthalene (TDN) and concentrations in young Riesling and non-Riesling wines. J. Agric. Food Chem. 2012, 60, 2998–3004. [Google Scholar] [CrossRef]
  43. Ferreira, V.; Ortín, N.; Escudero, A.; López, R.; Cacho, J. Chemical characterization of the aroma of Grenache rosé wines: Aroma extract dilution analysis, quantitative determination, and sensory reconstitution studies. J. Agric. Food Chem. 2002, 50, 4048–4054. [Google Scholar] [CrossRef]
  44. Gómez-Míguez, M.J.; Cacho, J.F.; Ferreira, V.; Vicario, I.M.; Heredia, F.J. Volatile components of Zalema white wines. Food Chem. 2007, 100, 1464–1473. [Google Scholar] [CrossRef]
  45. Vilanova, M.; Genisheva, Z.; Masa, A.; Oliveira, J.M. Correlation between volatile composition and sensory properties in Spanish Albariño wines. Microchem. J. 2010, 95, 240–246. [Google Scholar] [CrossRef]
  46. Culleré, L.; Escudero, A.; Cacho, J.; Ferreira, V. Gas chromatography− olfactometry and chemical quantitative study of the aroma of six premium quality Spanish aged red wines. J. Agric. Food Chem. 2004, 52, 1653–1660. [Google Scholar] [CrossRef] [PubMed]
  47. Peinado, R.A.; Moreno, J.; Bueno, J.E.; Moreno, J.A.; Mauricio, J.C. Comparative study of aromatic compounds in two young white wines subjected to pre-fermentative cryomaceration. Food Chem. 2004, 84, 585–590. [Google Scholar] [CrossRef]
  48. Gammacurta, M.; Tempere, S.; Marchand, S.; Moine, V.; de Revel, G. Ethyl 2-hydroxy-3-methylbutanoate enantiomers: Quantitation and sensory evaluation in wine. OENO One 2018, 52, 57–65. [Google Scholar] [CrossRef] [Green Version]
  49. Piombino, P.; Moio, L.; Genovese, A. Orthonasal vs. retronasal: Studying how volatiles’ hydrophobicity and matrix composition modulate the release of wine odorants in simulated conditions. Food Res. Int. 2019, 116, 548–558. [Google Scholar] [CrossRef]
  50. Azcarate, S.M.; Cantarelli, M.A.; Marchevsky, E.J.; Camiña, J.M. Evaluation of Geographic Origin of Torrontés Wines by Chemometrics. J. Food Res. 2013, 2, 48. [Google Scholar] [CrossRef]
  51. Moyano, L.; Zea, L.; Moreno, J.; Medina, M. Analytical study of aromatic series in sherry wines subjected to biological aging. J. Agric. Food Chem. 2002, 50, 7356–7361. [Google Scholar] [CrossRef]
  52. Guadagni, D.G.; Buttery, R.G.; Okano, S. Odour thresholds of some organic compounds associated with food flavours. J. Sci. Food Agric. 1963, 14, 761–765. [Google Scholar] [CrossRef]
  53. Sonni, F.; Moore, E.G.; Chinnici, F.; Riponi, C.; Smyth, H.E. Characterisation of Australian Verdelho wines from the Queensland Granite Belt region. Food Chem. 2016, 196, 1163–1171. [Google Scholar] [CrossRef] [Green Version]
  54. Buttery, R.G.; Turnbaugh, J.G.; Ling, L.C. Contribution of volatiles to rice aroma. J. Agric. Food Chem. 1988, 36, 1006–1009. [Google Scholar] [CrossRef]
  55. Cai, J.; Zhu, B.Q.; Wang, Y.H.; Lu, L.; Lan, Y.B.; Reeves, M.J.; Duan, C.Q. Influence of pre-fermentation cold maceration treatment on aroma compounds of Cabernet Sauvignon wines fermented in different industrial scale fermenters. Food Chem. 2014, 154, 217–229. [Google Scholar] [CrossRef] [PubMed]
  56. Escudero, A.; Gogorza, B.; Melús, M.A.; Ortín, N.; Cacho, J.; Ferreira, V. Characterization of the aroma of a wine from Maccabeo. Key role played by compounds with low odor activity values. J. Agric. Food Chem. 2004, 52, 3516–3524. [Google Scholar] [CrossRef] [PubMed]
  57. Ferreira, V.; Jarauta, I.; Ortega, L.; Cacho, J. Simple strategy for the optimization of solid-phase extraction procedures through the use of solid–liquid distribution coefficients: Application to the determination of aliphatic lactones in wine. J. Chromatogr. A 2004, 1025, 147–156. [Google Scholar] [CrossRef] [PubMed]
  58. Liu, H.; Bao, M.L.; Chen, H.L.; Li, Q. Impact of sucrose addition on the physiochemical properties and volatile compounds of “Shuangyou” red wines. J. Food Qual. 2017, 2017, 2926041. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Walter–Lieth climatic diagram for average temperature (red line) and total monthly precipitation (blue line) per month at Vrgorac weather station in 2019. On the left axis, values in black represent the absolute maximum and minimum temperature. At the top right, values in black represent the average values of mean temperature and annual precipitation.
Figure 1. Walter–Lieth climatic diagram for average temperature (red line) and total monthly precipitation (blue line) per month at Vrgorac weather station in 2019. On the left axis, values in black represent the absolute maximum and minimum temperature. At the top right, values in black represent the average values of mean temperature and annual precipitation.
Applsci 12 07327 g001
Figure 2. PCA performed using the concentrations of 16 volatile compounds with OAV > 1 discriminated among wine samples produced from five canopy-management treatments. (A) Projections of 16 volatile characteristics based on latent PC factors 1 and 2, which together explained 73.68% of total variability. (B) Principal component analysis (PCA) of 15 wines (C1, C2, C3: control; V1, V2, V3: defoliation at veraison; CT1, CT2, CT3: cluster thinning; PF1, PF2, PF3: defoliation at preflowering; and BS1, BS2, BS3: defoliation at berry set) on the plane defined by the first two principal components.
Figure 2. PCA performed using the concentrations of 16 volatile compounds with OAV > 1 discriminated among wine samples produced from five canopy-management treatments. (A) Projections of 16 volatile characteristics based on latent PC factors 1 and 2, which together explained 73.68% of total variability. (B) Principal component analysis (PCA) of 15 wines (C1, C2, C3: control; V1, V2, V3: defoliation at veraison; CT1, CT2, CT3: cluster thinning; PF1, PF2, PF3: defoliation at preflowering; and BS1, BS2, BS3: defoliation at berry set) on the plane defined by the first two principal components.
Applsci 12 07327 g002aApplsci 12 07327 g002b
Table 1. Yield components and basic physiochemical composition of Maraština musts and wines subjected to defoliation and cluster-thinning treatments.
Table 1. Yield components and basic physiochemical composition of Maraština musts and wines subjected to defoliation and cluster-thinning treatments.
PFBSVCTC
Yield components
Number clusters/vine8.9 ± 3.2 a7.5 ± 2.9 ab7.5 ± 2.4 ab5.9 ± 1.7 b7.7 ± 2.5 ab
Yield/vine (kg)1.3 ± 0.8 ab1.2 ± 0.5 b1.4 ± 0.6 ab1.4 ± 0.7 ab1.7 ± 0.8 a
Cluster weight (g)165.6 ± 99.2 b161.4 ± 54.2 b183.0 ± 72.0 ab229.1 ± 84.2 a222.5 ± 57.3 a
Berries per cluster225 ± 50 a192 ± 50 a224 ± 53 a253 ± 44 a227 ± 79 a
Must composition
TSS (°Brix) *18.9 ± 0.2 b18.7 ± 0.1 b18.9 ± 0.1 b19.6 ± 0.2 a19.5 ± 0.2 a
Total acidity **4.7 ± 0.0 a4.3 ± 0.1 c4.6 ± 0.1 b4.5 ± 0.0 b4.7 ± 0.1 a
pH3.66 ± 0.01 b3.74 ± 0.01 a3.74 ± 0.01 a3.74 ± 0.01 a3.52 ± 0.01 c
Wine composition
ASv (vol %) ***10.5 ± 0.1 b10.6 ± 0.1 b10.7 ± 0.2 b11.3 ± 0.1 a11.3 ± 0.1 a
Residual sugars (g L−1)1.0 ± 0.0 b1.0 ± 0.1 b1.1 ± 0.2 b1.6 ± 0.3 a1.2 ± 0.2 ab
Total extract (g L−1)16.5 ± 0.6 b16.9 ± 0.2 b17.8 ± 0.4 a18.5 ± 0.6 a19.0 ± 0.9 a
pH3.61 ± 0.02 c3.64 ± 0.04 bc3.76 ± 0.01 a3.69 ± 0.01 ab3.68 ± 0.05 bc
Total acidity (g L−1)4.8 ± 0.1 ab4.9 ± 0.1 a4.5 ± 0.1 b4.6 ± 0.1 b4.7 ± 0.1 ab
Volatile acidity (g L−1) ****0.5 ± 0.0 b0.5 ± 0.0 b0.7 ± 0.0 a0.6 ± 0.0 a0.6 ± 0.0 a
Data were subjected to one-way ANOVA; results are means of biological repetitions ± standard deviation. Means labeled with different letters (a, b, or c) within a row differed significantly according to Tukey’s test at p ≤ 0.05. * TSS: total soluble solids; ** total acidity is expressed as g L−1 tartaric acid equivalents; *** ASv: alcoholic strength by volume; **** volatile acidity is expressed as g L−1 acetic acid equivalents; PF: preflowering defoliation; BS: after berry set defoliation; V: veraison defoliation; CT: cluster thinning at veraison; C: untreated control.
Table 2. The concentrations of varietal aroma compounds of Maraština wines subjected to defoliation and cluster-thinning treatments.
Table 2. The concentrations of varietal aroma compounds of Maraština wines subjected to defoliation and cluster-thinning treatments.
Compound (µg L−1)PFBSVCTC
α-Pinene0.89 ± 0.07 b0.50 ± 0.06 c1.53 ± 0.07 a1.45 ± 0.15 a0.64 ± 0.04 c
Tetrahydrolinalool7.86 ± 0.45 ab5.53 ± 0.46 c5.75 ± 0.32 bc7.97 ± 1.23 a8.52 ± 1.08 a
Ethyl linalyl acetal0.05 ± 0.03 b0.08 ± 0.01 ab0.17 ± 0.20 ab0.09 ± 0.03 ab0.30 ± 0.01 a
cis-Linalool oxide, furan0.18 ± 0.03 b0.26 ± 0.04 b0.14 ± 0.02 b0.17 ± 0.05 b0.46 ± 0.07 a
Linalool14.98 ± 2.35 b12.27 ± 0.95 b12.35 ± 2.14 b20.40 ± 0.82 a15.29 ± 0.81 b
Terpinen-4-ol0.78 ± 0.24 c1.60 ± 0.41 b0.66 ± 0.30 c3.33 ± 0.28 a0.77 ± 0.06 c
Hotrienol2.74 ± 0.29 a0.80 ± 0.06 b2.39 ± 0.09 a2.67 ± 0.28 a0.78 ± 0.20 b
trans-β-Farnesene0.72 ± 0.17 c0.55 ± 0.05 c2.11 ± 0.31 b2.62 ± 0.22 a0.53 ± 0.01 c
cis-β-Farnesene0.95 ± 0.04 a1.00 ± 0.08 a1.00 ± 0.09 a0.92 ± 0.03 a1.07 ± 0.09 a
α-Farnesene0.89 ± 0.02 a0.94 ± 0.06 a0.93 ± 0.09 a0.91 ± 0.12 a1.08 ± 0.06 a
2,6-dimethyl-7-octene-2,6-diol4.32 ± 0.56 a4.04 ± 0.41 a2.50 ± 0.12 b2.36 ± 0.12 b5.23 ± 0.83 a
α-Terpineol1.16 ± 0.24 a1.50 ± 0.36 a1.39 ± 0.29 a1.31 ± 0.16 a1.22 ± 0.10 a
Citronellol12.02 ± 0.62 b12.61 ± 1.74 b12.31 ± 2.12 b19.26 ± 3.05 a15.87 ± 2.17 ab
Nerol0.57 ± 0.07 a0.50 ± 0.13 a0.44 ± 0.06 a0.42 ± 0.03 a0.49 ± 0.08 a
ε-Fenchene0.92 ± 0.11 a0.90 ± 0.04 a1.10 ± 0.11 a1.03 ± 0.09 a1.03 ± 0.11 a
Geraniol3.56 ± 0.41 b3.25 ± 0.11 b3.59 ± 0.40 b4.53 ± 0.38 a3.61 ± 0.17 b
6,7-dihydro-7-hydroxylinalool4.18 ± 0.01 b4.24 ± 0.14 b4.36 ± 0.25 b4.51 ± 0.33 b5.23 ± 0.28 a
Geranyl acetate2.64 ± 0.09 a2.55 ± 0.04 a2.68 ± 0.14 a2.42 ± 0.17 a2.78 ± 0.29 a
8-Hidroxylinalool5.93 ± 0.88 a9.83 ± 0.45 ab13.89 ± 1.04 a12.60 ± 3.88 a7.21 ± 0.43 b
Farnesol3.61 ± 0.45 c6.52 ± 0.46 a5.07 ± 0.09 b4.53 ± 0.19 bc5.21 ± 0.65 b
Neric acid2.57 ± 0.25 a2.40 ± 0.09 a2.50 ± 0.04 a3.78 ± 1.66 a2.62 ± 0.26 a
1,8-Terpin0.35 ± 0.04 d0.70 ± 0.08 bc0.90 ± 0.08 ab0.43 ± 0.11 cd1.13 ± 0.18 a
Neralidol0.66 ± 0.02 a0.70 ± 0.06 a0.82 ± 0.05 a0.73 ± 0.05 a0.70 ± 0.13 a
∑Terpenes72.52 ± 3.41 b73.34 ± 2.73 b78.58 ± 4.63 b98.44 ± 8.52 a81.76 ± 2.59 b
TDN *0.72 ± 0.16 a0.63 ± 0.03 a0.61 ± 0.10 a0.74 ± 0.17 a0.65 ± 0.16 a
β-Damascenone3.61 ± 0.44 a2.38 ± 0.40 a2.43 ± 0.25 a3.10 ± 1.23 a3.46 ± 0.66 a
Dihydroactinidiolide0.83 ± 0.06 a1.16 ± 0.22 a0.82 ± 0.20 a0.93 ± 0.18 a1.02 ± 0.09 a
∑C13-norisoprenoids5.15 ± 0.33 a4.17 ± 0.60 a3.86 ± 0.12 a4.77 ± 1.55 a5.13 ± 0.68 a
Data were subjected to one-way ANOVA; results are means of three biological repetitions ± standard deviation. Means labeled with different letters (a, b, c, or d) within a row differed significantly according to Tukey’s test at p ≤ 0.05. PF: preflowering defoliation; BS: after berry set defoliation; V: veraison defoliation; CT: cluster thinning at veraison; C: untreated control; * TDN: 1,1,6-trimethyl-1,2-dihydronaphthalene.
Table 3. The concentration of fermentation aromas in Maraština wines subjected to defoliation and cluster-thinning treatments.
Table 3. The concentration of fermentation aromas in Maraština wines subjected to defoliation and cluster-thinning treatments.
Compound (µg L−1)PFBSVCTC
Isobutyl acetate22.18 ± 1.60 c16.62 ± 0.47 d24.30 ± 1.06 c27.89 ± 0.95 b31.37 ± 1.89 a
Ethyl butanoate165.94 ± 9.33 c145.79 ± 2.41 d210.23 ± 3.08 a221.89 ± 2.52 a187.91 ± 5.45 b
Ethyl 2-methylbutanoate4.48 ± 0.39 a3.39 ± 0.22 b2.09 ± 0.08 c4.51 ± 0.26 a4.43 ± 0.26 a
Ethyl 3-methylbutanoate9.55 ± 0.44 b7.78 ± 0.81 bc6.12 ± 0.64 c12.06 ± 0.89 a12.16 ± 0.57 a
Isoamyl acetate955.88 ± 37.78 b650.86 ± 20.22 c1373.66 ± 110.61 a1403.59 ± 4.97 a1031.45 ± 22.90 b
Ethyl hexanoate258.46 ± 24.48 c322.49 ± 8.91 bc440.27 ± 36.02 a364.94 ± 28.57 b351.05 ± 11.05 b
Hexyl acetate15.30 ± 1.12 c17.28 ± 0.26 c27.15 ± 1.91 a21.95 ± 1.33 b14.42 ± 1.09 c
Ethyl lactate586.12 ± 15.73 a463.50 ± 29.66 b358.71 ± 13.27 c561.73 ± 36.89 a439.51 ± 33.59 b
Ethyl 2-hydroxy-3-methylbutanoate95.68 ± 3.49 ab61.26 ± 1.12 d101.22 ± 1.31 a90.06 ± 1.97 b81.68 ± 2.55 c
Ethyl octanoate516.94 ± 20.99 ab362.08 ± 32.57 c547.21 ± 47.58 a444.28 ± 25.05 b438.25 ± 11.46 bc
Methyl 2-furoate149.72 ± 0.92 a150.84 ± 0.67 a149.91 ± 0.10 a149.74 ± 1.11 a150.38 ± 0.99 a
Ethyl furoate1.52 ± 0.19 a1.10 ± 0.08 abc0.80 ± 0.15 c1.35 ± 0.25 ab1.05 ± 0.05 bc
Ethyl decanoate110.80 ± 7.23 b82.87 ± 5.16 c150.04 ± 6.29 a145.90 ± 8.78 a97.01 ± 3.79 bc
Diethyl succinate267.11 ± 35.33 bc205.91 ± 4.96 d252.80 ± 21.43 cd412.96 ± 19.04 a319.80 ± 6.12 b
Methyl geranoate2.67 ± 0.25 ab2.37 ± 0.35 b3.40 ± 0.28 a2.52 ± 0.37 b2.56 ± 0.05 b
2-Phenylethyl acetate6.55 ± 0.39 bc5.22 ± 0.30 c8.43 ± 1.14 a7.60 ± 0.37 ab7.04 ± 0.43 ab
Ethyl 3-hydroxyhexanoate1.68 ± 0.34 a0.83 ± 0.14 b1.62 ± 0.14 ab2.34 ± 0.42 a2.02 ± 0.35 a
Ethyl hexadecanoate0.26 ± 0.04 ab0.14 ± 0.04 c0.08 ± 0.02 c0.18 ± 0.01 bc0.29 ± 0.06 a
∑Esters3,170.85 ± 79.85 c2,500.35 ± 80.45 d3,658.03 ± 48.58 b3,875.49 ± 28.09 a3,172.37 ± 42.68 c
Isobutanol2,555.41 ± 126.36 bc2,225.97 ± 18.98 c2,545.03 ± 121.52 bc3,235.33 ± 255.55 a2,752.73 ± 103.98 b
1-Butanol36.84 ± 1.83 b39.20 ± 3.78 b40.86 ± 1.13 b58.69 ± 1.35 a41.10 ± 3.14 b
2-Methyl-1-butanol9,752.26 ± 126.73 c7,997.82 ± 85.59 d10,115.13 ± 112.48 bc12,452.80 ± 466.25 a10,787.12 ± 381.73 b
Isoamyl alcohol2,397.85 ± 232.52 c1,609.71 ± 100.10 c9,411.80 ± 503.45 a7,487.32 ± 408.50 b8,482.01 ± 420.15 a
1-Pentanol5.36 ± 0.46 c5.34 ± 0.03 c8.21 ± 0.26 b11.35 ± 2.12 a6.32 ± 0.50 bc
3-Methyl-1-pentanol132.52 ± 3.26 ab124.94 ± 5.19 b103.70 ± 4.68 c146.56 ± 4.89 a100.66 ± 9.17 c
1-Hexanol0.42 ± 0.03 b0.49 ± 0.09 b0.60 ± 0.08 b1.19 ± 0.23 a0.30 ± 0.03 b
trans-3-Hexen-1-ol22.64 ± 2.66 c24.79 ± 0.94 c32.00 ± 1.06 b44.17 ± 2.94 a28.59 ± 3.10 bc
2-Ethyl-1-hexanol1.16 ± 0.16 a0.98 ± 0.07 a1.23 ± 0.14 a1.35 ± 0.21 a1.00 ± 0.12 a
1-Octen-3-ol1.56 ± 0.27 b1.37 ± 0.24 b2.56 ± 0.18 a2.78 ± 0.25 a2.25 ± 0.19 a
1-Heptanol11.62 ± 0.63 a8.96 ± 0.19 b4.86 ± 0.84 c6.34 ± 0.62 c11.69 ± 0.98 a
1-Octanol7.88 ± 1.29 a7.97 ± 1.11 a10.57 ± 1.34 a10.44 ± 1.07 a11.08 ± 1.45 a
1-Nonanol0.44 ± 0.14 bc0.14 ± 0.02 c0.75 ± 0.15 ab0.90 ± 0.09 a0.91 ± 0.15 a
1-Decanol4.76 ± 0.17 a1.86 ± 0.55 c3.41 ± 0.44 b4.10 ± 0.62 ab3.76 ± 0.18 ab
Phenylethyl alcohol5,126.42 ± 154.61 a4,917.57 ± 69.47 a5,055.14 ± 225.33 a4,861.04 ± 94.37 a5,272.90 ± 184.53 a
∑Higher alcohols20,057.14 ± 201.50 b16,967.12 ± 130.06 c27,335.85 ± 709.84 a28,324.36 ± 215.46 a27,502.42 ± 367.05 a
1,2-Propanediol1.65 ± 0.21 b3.34 ± 0.38 a1.50 ± 0.21 b3.75 ± 0.44 a1.93 ± 0.22 b
2,3-Butanediol31.07 ± 1.51 ab23.52 ± 2.63 c32.46 ± 2.61 ab33.72 ± 1.82 a27.01 ± 1.17 bc
Furfuryl alcohol2.40 ± 0.15 a2.59 ± 0.35 a2.68 ± 0.28 a2.66 ± 0.32 a2.46 ± 0.31 a
Benzyl alcohol8.03 ± 1.35 b7.81 ± 1.73 b10.65 ± 0.72 ab11.79 ± 1.61 a7.75 ± 0.57 b
1,4-Butanediol1.31 ± 0.48 a1.25 ± 0.16 a0.60 ± 0.06 b1.01 ± 0.03 ab1.22 ± 0.13 ab
∑Other alcohols44.45 ± 3.04 bc38.51 ± 4.06 c47.89 ± 3.13 ab52.93 ± 0.48 a40.37 ± 1.81 bc
2-Octenal274.53 ± 19.69 b219.22 ± 8.76 c315.23 ± 6.75 a280.30 ± 8.58 b263.40 ± 12.70 b
2,4-Heptadienal0.38 ± 0.07 ab0.32 ± 0.04 b0.45 ± 0.04 ab0.62 ± 0.20 a0.40 ± 0.03 ab
Furfural1.11 ± 0.11 c0.72 ± 0.04 d1.41 ± 0.03 c2.25 ± 0.21 a1.7 ± 0.16 b
Decanal1.45 ± 0.38 a1.30 ± 0.06 a1.37 ± 0.11 a1.54 ± 0.40 a1.61 ± 0.32 a
Benzaldehyde8.06 ± 1.79 b6.26 ± 0.38 b11.91 ± 0.33 a12.47 ± 0.70 a6.91 ± 0.27 b
Benzeneacetaldehyde7.42 ± 0.37 b8.55 ± 0.29 a6.21 ± 0.60 c5.77 ± 0.22 c8.48 ± 0.13 a
Neral8.84 ± 0.08 a7.13 ± 1.14 ab4.11 ± 0.36 c3.33 ± 0.83 c6.2 ± 0.86 b
2,4-Decadienal0.20 ± 0.03 a0.15 ± 0.02 ab0.08 ± 0.03 bc0.05 ± 0.02 c0.19 ± 0.06 a
∑Aldehydes301.99 ± 19.32 b243.65 ± 8.11 c340.78 ± 5.61 a306.33 ± 7.40 b289.02 ± 13.77 b
Propanoic acid3.01 ± 0.18 c4.88 ± 0.19 a3.84 ± 0.36 bc4.29 ± 0.42 ab4.03 ± 0.35 b
2-Methylpropionic acid3.52 ± 0.40 b2.63 ± 0.23 c4.60 ± 0.45 a2.79 ± 0.17 bc3.19 ± 0.19 bc
Butanoic acid2.71 ± 0.36 c3.80 ± 0.11 a3.02 ± 0.09 bc3.17 ± 0.24 bc3.30 ± 0.18 ab
Isovaleric acid7.67 ± 1.37 a3.02 ± 0.41 b5.94 ± 1.17 ab5.29 ± 0.52 ab8.09 ± 1.58 a
Hexanoic acid1,112.94 ± 103.53 b1,315.46 ± 85.31 ab1,493.72 ± 176.20 a1,455.99 ± 127.64 a1,307.05 ± 95.62 ab
Heptanoic acid3.57 ± 0.33 b8.95 ± 0.85 a3.69 ± 0.20 b8.75 ± 0.33 a9.07 ± 0.83 a
Nonanoic acid8.24 ± 0.16 a7.50 ± 0.77 a8.21 ± 0.26 a8.45 ± 0.46 a8.24 ± 0.43 a
Decanoic acid864.94 ± 51.74 ab737.23 ± 83.28 bc863.64 ± 33.27 ab876.59 ± 22.24 a696.02 ± 43.03 c
∑Fatty acids2,006.59 ± 75.06 b2,083.47 ± 151.65 ab2,386.66 ± 155.65 a2,365.30 ± 126.77 a2,038.99 ± 137.26 ab
Eugenol0.35 ± 0.03 ab0.14 ± 0.05 c0.22 ± 0.07 bc0.39 ± 0.04 a0.49 ± 0.09 a
4-Ethylphenol0.04 ± 0.02 a0.06 ± 0.04 a0.06 ± 0.03 a0.07 ± 0.02 a0.05 ± 0.04 a
4-Vinylguaiacol30.54 ± 6.03 b32.38 ± 3.01 b35.02 ± 1.81 b50.09 ± 0.80 a50.66 ± 4.77 a
Vanillin0.92 ± 0.09 bc1.26 ± 0.14 a0.80 ± 0.07 bc0.75 ± 0.15 c1.12 ± 0.16 ab
Homovanillyl alcohol30.90 ± 6.13 a30.91 ± 5.22 a37.64 ± 4.46 a37.15 ± 7.80 a33.29 ± 6.79 a
∑Volatile phenols62.76 ± 12.16 b64.76 ± 4.53 b73.74 ± 3.41 ab88.45 ± 6.84 a85.61 ± 3.20 a
γ-Butyrolactone311.44 ± 30.80 c380.71 ± 34.38 bc392.35 ± 15.52 b562.58 ± 27.68 a358.31 ± 25.18 bc
γ–Hexalactone1.73 ± 0.15 b1.63 ± 0.43 b2.64 ± 0.32 a2.50 ± 0.29 ab2.34 ± 0.41 ab
γ-Octalactone0.33 ± 0.04 a0.36 ± 0.01 a0.40 ± 0.04 a0.41 ± 0.02 a0.37 ± 0.06 a
γ-Nonalactone5.52 ± 0.39 ab5.07 ± 0.33 b6.34 ± 0.65 ab6.57 ± 0.22 ab6.94 ± 1.02 a
γ–Decalactone1.54 ± 0.16 b1.59 ± 0.13 b1.79 ± 0.19 b2.24 ± 0.15 a1.89 ± 0.15 ab
δ-Decalactone2.11 ± 0.22 a2.03 ± 0.02 a1.96 ± 0.12 a2.17 ± 0.13 a2.17 ± 0.29 a
γ-Undecalactone0.58 ± 0.10 a0.62 ± 0.04 a0.56 ± 0.09 a0.66 ± 0.07 a0.58 ± 0.06 a
∑Lactones323.24 ± 30.85 c392.00 ± 33.76 bc406.04 ± 16.59 b577.13 ± 28.07 a372.59 ± 26.69 bc
2-Pentylfuran228.35 ± 11.34 a224.63 ± 2.79 a201.88 ± 14.39 a213.08 ± 17.82 a207.61 ± 18.96 a
6-Methyl-5-hepten-2-one114.17 ± 4.70 a89.90 ± 5.40 b62.49 ± 4.16 c113.50 ± 14.41 a117.64 ± 4.10 a
Acetoin4.07 ± 0.22 c13.20 ± 1.64 a6.16 ± 0.75 bc7.04 ± 0.50 b5.03 ± 0.99 bc
Acetylfurane1.13 ± 0.20 a1.03 ± 0.14 a1.10 ± 0.14 a1.20 ± 0.21 a1.12 ± 0.06 a
∑Other compounds110.92 ± 15.03 b107.11 ± 7.81 b124.90 ± 6.52 ab145.05 ± 6.58 a129.67 ± 4.89 ab
Data were subjected to one-way ANOVA; results are means of three biological repetitions ± standard deviation. Means labeled with different letters (a, b, c, or d) within a row differed significantly according to Tukey’s test at p ≤ 0.05. PF: preflowering defoliation; BS: after berry set defoliation; V: veraison defoliation; CT: cluster thinning at veraison; C: untreated control.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mucalo, A.; Lukšić, K.; Budić-Leto, I.; Zdunić, G. Cluster Thinning Improves Aroma Complexity of White Maraština (Vitis vinifera L.) Wines Compared to Defoliation under Mediterranean Climate. Appl. Sci. 2022, 12, 7327. https://doi.org/10.3390/app12147327

AMA Style

Mucalo A, Lukšić K, Budić-Leto I, Zdunić G. Cluster Thinning Improves Aroma Complexity of White Maraština (Vitis vinifera L.) Wines Compared to Defoliation under Mediterranean Climate. Applied Sciences. 2022; 12(14):7327. https://doi.org/10.3390/app12147327

Chicago/Turabian Style

Mucalo, Ana, Katarina Lukšić, Irena Budić-Leto, and Goran Zdunić. 2022. "Cluster Thinning Improves Aroma Complexity of White Maraština (Vitis vinifera L.) Wines Compared to Defoliation under Mediterranean Climate" Applied Sciences 12, no. 14: 7327. https://doi.org/10.3390/app12147327

APA Style

Mucalo, A., Lukšić, K., Budić-Leto, I., & Zdunić, G. (2022). Cluster Thinning Improves Aroma Complexity of White Maraština (Vitis vinifera L.) Wines Compared to Defoliation under Mediterranean Climate. Applied Sciences, 12(14), 7327. https://doi.org/10.3390/app12147327

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop