- 1Department of Agricultural, Environmental and Food Sciences, University of Molise, Campobasso, Italy
- 2Department of Agriculture, Food and Environment, University of Catania, Catania, Italy
- 3Department of Agricultural Sciences, Division of Grape and Wine Sciences, University of Naples Federico II, Avellino, Italy
- 4Institute for Mediterranean Agricultural and Forestry Systems, National Reaserch Council, Portici, Italy
- 5Council for Agricultural Research and Economics, Research Centre for Olive, Fruit and Citrus Crops, Acireale, Italy
The cultivated grapevine, Vitis vinifera subsp. vinifera, possesses a rich biodiversity with numerous varieties. Each variety adapts differently to varying pedoclimatic conditions, which greatly influence the terroir expression of wine regions. These conditions impact vine growth, physiology, and berry composition, ultimately shaping the unique characteristics and typicity of the wines produced. Nowadays, the potential of the different adaptation capacities of grape varieties has not yet been thoroughly investigated. We addressed this issue by studying two grape varieties, Aglianico and Cabernet Sauvignon, in two different pedoclimatic conditions of Southern Italy. We evaluated and compared the effect of different pedoclimatic conditions on plant physiology, the microbial quality of grapes using Next-Generation Sequencing (NGS) technology, the expression trends of key genes in ripe berries and the concentration of phenolic compounds in grapes and wines by HPLC-MS, HPLC-DAD, NMR and spectrophotometric analyses. Metabolomic and microbiome data were integrated with quantitative gene expression analyses to examine varietal differences and plasticity of genes involved in important oenological pathways. The data collected showed that the phenotypic response of studied grapes in terms of vigor, production, and fruit quality is strongly influenced by the pedoclimatic conditions and, in particular, by soil physical properties. Furthermore, Aglianico grape variety was more influenced than the Cabernet Sauvignon by environmental conditions. In conclusion, the obtained findings not only reinforce the terroir concept and our comprehension of grape’s ability to adapt to climate variations but can also have implications for the future usage of grape genetic resources.
1 Introduction
Among high-income crops, grapevine (Vitis vinifera subsp. vinifera) is considered one of the most challenged by climate change (Jones and Webb, 2010; Goode, 2012). In the Mediterranean area, a decrease in rainfall associated with an increase in temperature is expected (IPCC, 2022), leading to a shifting of suitable areas able to satisfy the specific thermal request of grapevine cultivars and to a marked effect on soil water availability. The combination of stressor conditions, and the expected water scarcity, will directly affect the grape and wine quality (Duchêne, 2016; Bonfante et al., 2017; Bonfante et al., 2018). Nevertheless, cultivated grapes possess a specific biological capability known as phenotypic plasticity, allowing them to adapt to climate changes; this adaptation is attributed to a comprehensive reorganization of the entire genome. Such reorganization often involves changes at the genomic, epigenetic, transcriptional, and network levels and is triggered by environmental factors. It is best illustrated by the concept of “terroir”, which encompasses the combined influences of varietal attributes, climate, soil conditions, winemaking practices, and the multitude of interactions among them (Bonfante and Brillante, 2022). Previous investigations showed a high level of differentiation among vineyards based on the varieties cultivated or by their area of provenance (Anesi et al., 2015; Paim Pinto et al., 2016; Xie et al., 2017; Dal Santo et al., 2018; Di Vaio et al., 2020). Concerning the varieties, a differentiation can be considered between international and indigenous ones. The former is more uniform and, thus, more resilient to environmental changes than the latter, indeed, international varieties tend to take over indigenous ones across the globe (Wolkowich et al., 2018). This knowledge gap needs to be filled to strengthen solutions towards declining grape diversity richness worldwide. Taken as a whole, the theme of phenotypic plasticity is complex, and researchers did not arrive at the same conclusions, with results sometimes controversial. This is likely linked to the fact that most grapevine studies focused on a single or a few components of the environment and plant responses. For example, Perin et al. (2020) and Vannozzi et al. (2023) observed genotype-dependent responses to soil factors in grapes. They highlighted the importance of specific soil properties, in influencing the metabolite profile of berry skin and found out that Corvina had higher biochemical plasticity compared to Glera, while Glera had higher transcriptional plasticity with respect to Corvina.
In recent years, few studies investigated the effects of various pedoclimatic conditions on berry and wine quality (Bonfante et al., 2018; Dinu et al., 2021; Rodrigues et al., 2023), while Nicolosi et al. (2022) have compared the phenotypic plasticity of indigenous grape varieties over international ones. In this work, we aimed to understand how pedoclimatic conditions can influence berries and wine characteristics, transcriptional regulation and the berries epiphytic fungal community of an international (Cabernet Sauvignon) and an indigenous (Aglianico) red grape variety along a latitudinal climatic gradient.
2 Materials and methods
2.1 Study areas
The study areas are in two hilly environments of southern Italy at different latitudes and elevations, Molise (41° 42’; 606 m a.s.l.), and Sicily (37° 40’; 400-700 m a.s.l.), in farms oriented to the production of high-quality wines. The parent material of soil is different between the experimental sites, clayey marls in Molise, and volcanic material in Sicily that led to the development of Vertisols and Andosols, respectively. The studied grapevines were Aglianico (AGL) and Cabernet Sauvignon (CS) cultivars planted between 2008 and 2010 on 140 Ruggieri rootstocks. In Molise, vines were planted in north-south rows, trained according to espalier Guyot, placed in a flat area. In Sicily, both vineyards were bush trained at 0.5 m, with two to six main branches (each branch was spur-pruned to one spur, with two buds per spur), and placed in anthropic terraces.
2.2 Pedoclimatic characterization and weather monitoring
In each experimental site, the representative soil profile has been described and analyzed. Soil profiles were described according to IUSS (IUSS Working Group WRB, 2014). Chemical analyses were performed according to the official methods of the Italian Ministry of Agriculture and Forestry (Colombo and Miano, 2015). Disturbed and undisturbed soil samples (volume 750 cm3) were collected from each soil horizon to determine chemical and physical soil characteristics. The hydraulic properties were determined in the laboratory according to the procedures reported by Basile et al. (2012) and Bonfante et al. (2010). Climatic characterization of each site was realized on the reference climate dataset 1971-2005 of Italy (spatial resolution about 8 km) realized and optimized on regional high-resolution observational data sets by Bucchignani et al. (2015). This choice has been necessary because there was no available local long-weather dataset reliable near the sites. Weather information (temperature, rainfall, solar radiations, wind speed, air humidity, etc.) was collected using an in-situ weather station (netsense.com) placed in each experimental vineyard within the canopy. Daily reference evapotranspiration (ET0) was estimated according to the equation of Penman-Monteith equation (Monteith, 1965) and the bioclimatic index of Amerine and Winkler (A&W; Amerine and Winkler, 1944) applied to compare the different thermal regimes.
2.3 Berry sample collection
Berries were collected at the full ripening stage (BBCH stage 89, Lorenz et al., 1995) during the 2020 growing season for biochemical, microbiological, and transcriptional analyses, all at the same time. Thirty clusters were harvested from different positions of the vineyard and from random positions on the plant to ensure the representation of the entire vineyard. At least ten berries were randomly selected from different parts of the cluster, avoiding those with visible damage and/or signs of pathogen infection, and pooled with berries from the other plants. For subsequent analyses, three independent pools (biological replicates) of 50 (for transcriptional analysis) or 100 (for biochemical and microbiological analyses) whole berries were selected, immediately frozen in liquid nitrogen, and stored at -80°C for subsequent analyses. For each analysis, four experimental samples were considered: Aglianico Molise (AM), Aglianico Sicily (AS), Cabernet Sauvignon Molise (CM), and Cabernet Sauvignon Sicily (CS).
2.4 Morpho-physiological measurements and maturity index calculation
For each cultivar/terroir combination and time of measurement, 21 fully exposed primary leaves, and 21 lateral leaves were collected. All the measurements have been conducted at the end of flowering, at the berry’s pea-sized stage, and at berry softening (vèraison). The total leaf area (TLA) was estimated through the measurement of leaf area (LA) of the main and lateral leaves realised in the lab with the area meter LI-3100 (Licor, Inc., Lincoln, Nebraska). Plant water status was determined in each experimental site using midday leaf water potential (LWP) measurement realised with a Schöelander pressure chamber (Soil Moisture Equipment Corp., Sta. Barbara, CA, USA; Scholander et al., 1965). Bunches were collected from 5 vines per treatment on each block for yield assessment. A sample of five bunches per treatment on each block was used to determine the bunch weight and both the number and weight of the berry. A 100-berry sample per experimental block was divided into three sub-samples, crushed with a manual press. The free-run juice was utilised to determine the total soluble solids (TSS) measured by a digital refractometer with temperature correction (RX-5000 Atago Co., Ltd., Bellevue, WA, USA). Must pH and titratable acidity (TA) were determined using an automatic titrator (Titrino model 798, Metrohm, Riverview, FL, USA) on 5.0 mL juice samples titrated against 0.1 M NaOH up to pH 8.2. TA was expressed as g/L of tartaric acid equivalents. The maturity index was estimated at harvest in 2020 (soluble solids -Brix° and pH), using the following formula (Lana et al., 2021):
2.5 Chemical analyses and spectrophotometric measurements of grape extracts and wines
Separation and extraction of berry components was carried out in duplicate simulating the maceration process necessary to produce red wines (Mattivi et al., 2002). Monomeric anthocyanins were measured according to the OIV Compendium of International Methods of Analysis of Wine and Musts (https://www.oiv.int/OIV-MA-AS315-11). Total flavanols were determined as vanillin-reactive flavans, according to Gambuti et al. (2015). Total phenolic compounds were dosed as iron-reactive phenolics (Harbertson, 2003). A large-scale extraction (300 g of fresh product) was carried out for each collected sample of red grapes. Destemmed grapes were crushed and extracted twice overnight Skin with 500 mL of H2O:EtOH 2:8 (v/v) mixture. The hydroethanolic extracts from each sample were concentrated under a vacuum and partitioned first against Ethyl Acetate (EtOAc) and then n-butanol 2-mL aliquots of each obtained extract were lyophilized and investigated by NMR. Likewise, 1-mL aliquots from each extract were lyophilized, solubilized in methanol, and subjected to MS analysis. NMR experiments were run on a Varian Unity Inova 700 spectrometer equipped with a 13C Enhanced HCN Cold Probe and by using a Shigemi 5 mm MR tubes. CD3OD (δH 3.31; δC 49.0) was selected as deuterated solvent. Standard Varian pulse sequences were employed for the respective classes of spectra. All NMR data reported in the text were derived from 1D 1H spectra and from 2D experiments, including COSY and TOCSY. Relative quantitation of each identified compound was conducted by selecting representative NMR signals and their areas determined by integration. For Mass Spectrometry analysis, a linear ion trap LTQ Orbitrap XL hybrid Fourier Transform MS (FTMS) instrument equipped with an ESI ION MAX source (Thermo-Fisher, Waltham, MA, USA) and coupled to an Agilent 1100 LC binary system was used. All the employed solvents were of laboratory grade. Following the NMR and High-Resolution MS data for the identified compounds: Catechin: HRESIMS m/z 291.0800 [M+H]+ corresponding to C15H15O6. Epicatechin: HRESIMS m/z 291.0857 [M+H]+ corresponding to C15H15O6. Gallic acid: HRESIMS m/z 169.0169 [M-H]- corresponding to C7H5O5. Syringic acid: HRESIMS m/z 197.0501 [M-H]- corresponding to C9H9O5. Quercetin: HRESIMS m/z 303.1858 [M+H]+ corresponding to C15H11O7. NMR data for the above compounds were consistent with those reported in the literature (Forino et al., 2019; Gambuti et al., 2020).
Analyses of monomeric anthocyanins were performed by an HPLC Shimadzu LC10 ADVP apparatus (Shimadzu Italy, Milan, Italy) equipped with an SCL-10AVP system controller, two LC-10ADVP pumps to create the needed solvent gradient, an SPD-M 10 AVP detector and an injection system full Rheodyne model 7725 (Rheodyne, Cotati, CA, USA). The HPLC solvents were the following: solvent A: water milli- Q (Sigma-Aldrich, Milan, Italy)/formic acid (Sigma- Aldrich ≥ 95%)/acetonitrile (Sigma-Aldrich≥ 99.9%) (87:10:3) v/v; solvent B: water/formic acid/acetonitrile (40:10:50) v/v. The gradient was: zero-time conditions 94% A and 6% B; after 15 min, the pumps were adjusted to 70% A and 30% B, at 30 min to 50% A and 50%B, at 35 min to 40% A and 60% B, at 41 min, end of the analysis, to 94% A and 6% B. 5 minutes re-equilibration time were applied before the subsequent analysis. The column used for the analyses was a water spherisorb column (C 18, Silica particle substrate, ODS2 250 x 4.6 mm, 5 μm particles diameter, 80 Å pore size) with a pre-column was used. 50 µL of calibration standards or wine was injected into the column. The absorbance signals at 520 nm was detected. Detector sensitivity was 0.01 Absorbance units full scale (AUFS). All the samples were filtered through 0.45 µm Durapore membrane filters (Millipore-Ireland) into glass vials and immediately injected into the HPLC system. The calibration curve was obtained by injecting 5 solutions (in triplicate) containing increasing concentrations of malvidin-3-monoglucoside (Extrasynthese, Lyon, France).
2.6 Grapes fungal community analysis
Fungal community composition was analysed by a culture-independent approach using next-generation sequencing (NGS). The grape samples from every vineyard were collected in duplicate, immediately transported to the laboratory and processed. Total genomic DNA was extracted using the Stool DNA Isolation Kit (Norgen, Biotek Corp., Thorold, ON, Canada) according to the manufacturer’s instructions. The concentration and purity of the extracted nucleic acids were visualised and quantified, after agarose gel electrophoresis, by Nanodrop (NanoDropTM 2000/2000c; Thermo Fisher Scientific, Italy). DNA quantity was standardized to a concentration of 10 ng/μl. The ITS2 (internal transcribed spacer) region of the rRNA was amplified using primers 2024F and 2409R. Amplicons were sequenced using the Illumina MiSeq PE300 platform (Illumina, San Diego, CA, USA) at Biodiversa s.r.l. (Rovereto, Italy). Adapter sequences were removed using Cutadapt (Martin, 2011). Read quality was assessed using DADA2 (Callahan et al., 2016). The taxonomic references were assigned using trained OTUs at 99% from UNITE database version 8.2 (https://unite.ut.ee) within Qiime2 tools version 2020.2 (https://qiime2.org) (Bolyen et al., 2019). The raw data have been deposited in Mendeley Data with the accession number DOI: 10.17632/tp2rrdz8wp.2.
2.7 Gene expression analysis in grapes
Total RNA was isolated from 40 mg of ground whole berries, excluding seeds. The Spectrum™ Plant Total RNA kit (Sigma-Aldrich, St. Louis, MO, USA) was used following the manufacturer’s protocol with some modifications (400 mg of fresh tissue powder and an additional Column Wash step). The quantity and quality of the isolated RNA was measured using the NanoDrop ND-1000 spectrophotometer (Thermo Scientific, Wilmington, DE, USA) and Qubit 2.0 fluorometer (Life Technologies, Carlsbad, CA). For cDNA synthesis, 100 ng of each RNA sample was reverse transcribed using the SuperScript® III cDNA Synthesis Kit (Life Technologies), following the manufacturer’s protocol. Expression analysis was conducted by Real-Time quantitative PCR (RT-qPCR) using an SYBR Green method on a 7900HT Fast Real-Time PCR System (Applied Biosystems, Foster City, CA, USA). Each 15 μl PCR reaction contained 330 nM of each primer, 2 μl of 5-fold diluted cDNA and 7.5 μl of SYBR Green Mix (Applied Biosystems, Foster City, CA, USA). The SDS 2.3 and RQ Manager 1.2 software (both Applied Biosystems, Foster City, CA, USA) were used for data elaboration. Primer pairs used in the quantitative analysis were divided into five groups based on their roles: glutathione S-transferase as precursors of volatile thiols, terpene synthase enzymes for terpenoid biosynthesis, O-methyltransferases for the biosynthesis of methoxypyrazines, anthocyanin biosynthetic pathway genes and peroxidase for anthocyanin degradation (Sparvoli et al., 1994; Kobayashi et al., 2002; Dunlevy et al., 2010; Kobayashi et al., 2011; Battilana et al., 2017; Sicilia et al., 2020; Ji et al., 2021; Maritim et al., 2021; Wang et al., 2021). The relative expression value was calculated through the ΔΔCt method using the Sicilian berries as calibrators (Livak and Schmittgen, 2001; Villano et al., 2016). Three biological and three technical replicates were considered.
2.8 Statistical analysis
To evaluate the environmental difference between the experimental sites, the well-known t-test and f-test methods were applied to the weather information collected in 2020: mean temperature, day-night temperature excursion, ET0, maximum wind, the Amerine and Winkler index (A&W), and rainfall. Quantitative grapevine data were compared using Tukey’s least significant differences procedure; all the variance resulted homogeneous. When the variances were not homogeneous, data were analyzed using Kruskall–Wallis test. When results of the Kruskal–Wallis test were significant (p< 0.05), the significance of between-group differences was determined by the Bonferroni–Dunn test (5% significance level). These analyses were performed using XLSTAT (version 2013.6.04; Addinsoft, Paris, France). All the data are expressed as means ± standard deviation of four replicates (two experimental replicates x two analytical replicates).
3 Results
3.1 Pedological and climatic characterization
In the two experimental sites, different soils were identified and described as representative of selected vineyards (Figure 1 and Table 1). In particular, a deep clay soil in Molise site (Gleyic Calcic Vertisols) for both AGL and CS, and two shallow volcanic sandy soils (Vitric Andosols) in the Sicilian site. The available water capacity (AWC) of the first 100 cm differed among the soils, with the highest value expressed in Molise (205 mm) and the lowest in Sicily (127-132 mm). The topsoil layer cation exchange capacity (CEC) was high in both sites except for the Sicilian AGL vineyard. This data agrees with the young nature of the latter soil, less evolved compared to the others (96.1% of volcanic sand), and less enriched in organic matter compared to the Sicilian CS soil. Moreover, a rising of water from a suspended groundwater table was found in the Molise site.
3.2 Weather data and climatic characteristics of experimental sites
The results of climate analysis (RC 1971-2005 dataset) explained the effects of latitude in a Mediterranean climate (cold winters and mild summers), with the high latitude site (Molise) being the coldest compared to the Sicilian site (low latitude). Looking at the grapevine growing season (April-October 2020), the mean temperature was about 7.6% higher in Sicilian than in the Molise site. The behavior of the A&W index reflected the temperature differences, with a gap of 15% between Sicily and Molise. The ET0 for the two sites was quite similar, with a small difference of 9%, while the cumulative rainfall was lower for Sicily compared to Molise (average difference of -8%). Daily mean temperatures and rainfall collected in each experimental site during the monitored growing season (April-October) are reported in Figure 2. The mean temperature over the selected time was 17.2°C ( ± 0.5) for the Molise site, and 19.6°C ( ± 0.4) for the Sicilian one, and the ET0 was 748.4 mm ( ± 1.4) and 714.3 mm ( ± 1.8). The accumulated rainfall was 219.5 mm and 333.4 mm for Molise and Sicily, respectively. The thermal regime showed a higher value of the A&W index in Sicily (2056 GDD compared to 1600 GDD of Molise) with a different trend in the accumulation of GDD during the growing season (Figure 3). Statistical analysis of weather conditions among the experimental sites showed a significant difference between the two sites with a p< 0.05 (p mean temperature = 9,6E-27, p day-night temperature excursion = 0,0036, p ET0 = 0,0325, p maximum wind = 1,96E-51 and p A&W index = 4,29E-55), with an average difference of mean temperature = 2,4 C°, day-night temperature excursion = -0,9 C°, ET0 = -0,2 mm, maximum wind = -2,3 ms-1 and A&W index = 191,4 GDD. Regarding rainfall, there is no statistical difference between the two sites, but from the biological point of view, the difference of 100 mm over the growing season between Molise and Sicily can be considered significant in terms of the effect on plant responses.
Figure 2 Daily mean temperatures (T mean) and daily rainfall (mm) across the selected vineyards in the Molise and Sicily sites from 1 April 2020 to 31 October 2020.
Figure 3 Cumulated Amerine and Winkler index (Growing Degree Days - GDD) trend in the experimental sites (1 April - 31 October 2020).
3.3 Morphological measurements
During early summer, midday LWP was about less than -0.8 MPa in both environments, which is considered the negative threshold for well-irrigated vines (Olivo et al., 2009; Intrigliolo et al., 2012). Indeed, the ψ ≤ -1.0 MPa recorded from midsummer to the end of the season, are considered drought-like conditions occurring from the fruit set throughout the ripening stage. During the latter period in Sicily, the registered values of -1,79 MPa for Aglianico and -1.73 MPa for Cabernet Sauvignon were significantly below -1.2 MPa or even that is considered a condition of severe stress (Nicolosi et al., 2022). The LWP was measured at the end of flowering, at berry’s pea-sized stage and at berry softening. The values were significantly different between the first and third stages for Aglianico. LWP values were -1.31MPa (AS) and -1.79 MPa (CS) in Sicily, whereas in Molise they were higher. During the intermediate stage (berry’s pea-sized stage) no differences were detected. As regards Cabernet Sauvignon, the most stressed vines were observed in Molise during both the end of flowering and berry’s pea-sized stage (ψ= -1,12 MPa vs -0,95 MPa, and ψ= 1,32 MPa vs -1,15 MPa, respectively). At berry softening the most stressed vines were in Sicily (-1,73 MPa vs -1,17 MPa) (Nicolosi et al., 2022).
The total leaf area per vine (TLA) is reported in Table 2. For both cultivars and each phenological stage, the Molise site showed significantly higher values than the Sicilian site. However, in the latter, for Cabernet Sauvignon the TLA slightly increased during the growing season, while in Molise because of the shoot topping (midsummer) a decrease was recorded. For Cabernet Sauvignon, the above-mentioned decrease was observed at the berry softening stage due to the manual late basal leaf removal. The main and lateral shoot behaviour influenced the TLA. In fact, for Aglianico in Sicily, the TLA/main shoot ratio increased during the experimental period (from 0.23 m2 to 0.34 m2), whereas in Molise, a decrease in TLA/main shoot ratio was registered between the end of flowering stage and the berry pea-size stage. Concerning the lateral shoots in Molise, a substantial increase occurred at the berry softening stage. For Cabernet Sauvignon, significant differences between the two experimental sites were observed at the end of flowering and at the softening of berries, both for main and lateral shoots. Regarding the leaf area index, a significant difference was recorded at the berry pea-size stage for Aglianico and Cabernet Sauvignon and the highest value was observed in Sicily. On the contrary, the leaf area index of Aglianico in Molise showed a significant increase at the berry softening stage.
3.4 Berry and wine characteristics
The berries and wine characteristics are reported in Table 3. Aglianico and Cabernet Sauvignon had the highest yield, bunch weight, and berry number in Molise. As regards the berry weight, a significantly high value was detected in Molise for Aglianico, while the berry weight of Cabernet Sauvignon was similar in each region. The rachis weight was always higher in Sicily than in the Molise site. The technological maturity of grapes was evaluated by determining the content of soluble solids, titratable acidity, their ratio, and the pH in 2020. The main classes of phenolic compounds extractable from skins and seeds were also determined (Table 3). Total iron reactive phenolics, total flavanols and anthocyanins are higher in samples grown in Molise, while the content of flavanols strongly differs Aglianico from Cabernet. Data on alcohol content of wines were correlated to the content of sugars in grapes (Table 3). The lower values of titratable acidity of all wines with respect to respective grapes are expected considering tartaric precipitation during fermentation and post-fermentative aging of wines. Comparing data on total phenolics in grapes and wines (Table 3) in AS and AM phenolics were transferred from grapes to wine during winemaking but, especially for AM, a slight loss, probably due to adsorption on grape marcs, was detected. In contrast, in CS higher contents of phenolics compounds in wines with respect to berries were observed. This is easily explained since the method used to determine total phenolic compounds involved the reaction with iron and, thus, it is likely that during winemaking, for this grape variety, the formation of new phenolic structures highly reactive towards iron occurred. This is not surprising given the high reactivity of phenolic compounds among themselves and with fermentation metabolites such as acetaldehyde during winemaking. The lower values of total flavanols and anthocyanins detected in wines with respect to grapes is due to the great number of reactions of oxidation, polymerization, and precipitation that native phenolic compounds underwent. To identify the major polyphenols occurring in the analyzed samples, an NMR-based investigation was conducted. To this aim, both Aglianico and Cabernet Sauvignon grapes were extracted and partitioned as reported in the experimental section. Polyphenols were accumulated in the EtOAc extracts. Gallic acid and syringic acid were identified as phenolic acids, (+)-catechin and (-)-epicatechin and detected as the most abundant flavan-3-ols derivatives, while quercetin emerged as the most abundant flavonol compound. In Figure 4, 1H-NMR spectra of CM and CS EtOAc extracts display all the above proton resonances attributed to the identified polyphenols. Typical NMR signals attributable to hydroxycinnamic acids were identified, but due to their low intensities and occurrence in quite crowded spectral regions, unambiguous identification of specific molecules was prevented. It is to be underlined that anthocyanins were not analyzed by either NMR and MS techniques, as they were investigated by means of DAD-HPLC and spectrophotometric analyses, as reported above. All the above NMR-based structural hypotheses were corroborated by High Resolution Mass Spectrometry analysis (See Materials and methods).
Figure 4 Enlargements of the 1H-NMR spectra of CM (bottom) and CS (top) containing resonances assigned to the identified polyphenols. Red asterisks indicate resonances of quercetin, green asterisks indicate resonances of (+)-catechin, purple asterisks indicate overlapped resonances of (+)-catechin and (-)-epicatechin, magenta asterisk indicates syringic acid, and black asterisk gallic acid.
3.5 Characterization of fungal community
The data on the fungal community present in the harvested grapes is shown graphically in Figure 5. Overall, at harvest time the fungal populations, present in the vineyards of the examined sites, at the phylum level are very similar (Figure 5A). Ascomycota was the most abundant phylum followed by Basidiomycota. Other phylums such as Chytridiomycota are present only in low abundance. The analysis of the fungal community by high-throughput sequencing allowed to identify a total of 13 families and 12 genera. At the family level (Figure 5B), mean values of relative abundance revealed that the population was mainly represented by the Cladosporiaceae family with relative percentages between 11.2 (AS) and 48.6% (CS). The Saccotheciaceae family was widely present on the grapes analyzed with relative percentages between 11.4% (AM) and 33.5% (AS). Saccharomycodaceae was among the predominant family on AS grapes with a relative abundance of 43.9% while on the other grapes sampled this family was little detected. Widespread presence of Pleosporaceae with relative abundance between 2.5% (AS) and 18.9% (CM). At the genus level (Figure 5C), the main contaminants of grapes were fungi belonging to the genera Cladosporium, Aureobasidium, Alternaria and Sporobolomyces. A significant presence of the genus Hanseniaspora was recorded only on AS grapes (43.9%). The diffuse presence of Cladosporium was recorded on all the grapes analyzed, with relative abundances between 11.2% (AS) and 48.4% (CS).
Figure 5 Relative abundance of grape fungal community at phylum (A), family (B) and genus (C) level detected by NGS, Next Generation Sequencing.
3.6 Transcriptional analysis of Aglianico and Cabernet Sauvignon berries
The transcript abundance of representative key genes of phenylpropanoid metabolism is shown in Figure 6 as log2 of the fold change of the samples from Molise normalized to the samples from Sicily. Overall, both varieties showed an overexpression of F3H, LDOX, and UFGT, while DFR and PRX21 were downexpressed. MYBA and CHS were overexpressed in Cabernet Sauvignon and downexpressed in Aglianico. As regards volatile thiol, methoxypyrazine, and terpenoid metabolism in Molise we observed, with respect to Sicilian grapes, an overexpression of GST2, GST3 and GGT for volatile thiols, OMT1 and OMT2 for methoxypyrazine, and TPS56 for terpenoids, in both Aglianico and Cabernet Sauvignon. On the counterpart, OMT3 showed downregulation in both varieties (Figure S1).
Figure 6 Transcript abundance of genes involved in phenylpropanoid (MYBA, CHS, F3H, DFR, LDOX, UFGT, GST, and PRX21) metabolism. The values are expressed in log2 of Fold Change and are normalized on Sicilian samples. Statistically significant differences are asterisked (* p ≤ 0.05, ** p ≤ 0.01, *** p ≤ 0.001).
4 Discussion
The phenotypic plasticity of grape varieties is heavily influenced by the environmental conditions they inhabit, affecting their vigor, production, and fruit quality. Although grapevine plasticity and viticulture practices may provide advantages related to the adaptation of a cultivar to specific growing conditions, it may also cause irregular ripening (Selvaraj et al., 1994) and large inter-seasonal fluctuations (Clingeleffer, 2010), which are undesirable characteristics for wine making (Keller, 2010). In this study we tested how the phenotypic plasticity of two red berries grapes, Aglianico and Cabernet Sauvignon, can be influenced by the pedoclimatic conditions in terms of berries and wine characteristics, transcriptional regulation, and epiphytic fungal community. Overall, we found that both varieties were significantly influenced by the specific pedoclimatic conditions in Molise and Sicily regarding vegetative and productive responses. The sampling sites had a different capacity in soil function expression due to substantial differences in terms of chemical and physical characteristics. This function depends mainly on soil hydraulic properties that can drive grapevine responses and influence the quality of the grapes (Bonfante et al., 2017; Brillante et al., 2017; Brillante et al., 2018). In our experimental sites, the hydraulic properties were higher in Molise compared to Sicily as well as the weather variables affecting the soil water balance with about 9% ET0 in average difference. Such environmental differences between the sites led to different physiological and productive behaviors in both varieties. They can be considered responsible for the growth rate reductions of Aglianico in Sicily from the pea-size berries to the softening of berries stage. Indeed, a depressive effect of drought on shoot elongation and leaf emergence was observed, in agreement with negative water potential found. In contrast, in Molise, the Aglianico growth rate was continuous (mainly for lateral shoots), determining a very long vegetative and reproductive cycle with a high canopy efficiency. By contrast, Cabernet Sauvignon registered a less negative water potential value despite reduction of its leaf area. Furthermore, in Sicily, the growth of both the main and lateral shoots stopped, and the basal leaves lost their functional capabilities. It is crucial to emphasize that the observed discrepancies among sites regarding Cabernet Sauvignon are likely attributed more to the variations in training systems rather than the pedoclimatic conditions. Specifically, the reductions in leaf area index observed in Molise during mid-summer resulted from vine hedging carried out during that period. Consequently, the observed behavior could be due to a better adaptability of this cultivar, which adapted its phenotype in response to various restrictive growth conditions. Concerning Aglianico, the crop load was significantly impacted by the prevailing environmental conditions, leading to notably lower yields in Sicily due to an extended period of severe drought, compared to the Molise site, where abundant water availability fostered vigorous vegetative growth. Consequently, the genotype-environment interaction for Aglianico exhibited its strongest effect during the berry ripening phase. The NMR-based analysis of grape extracts highlighted interesting differences between Molise and Sicily regarding the relative abundance of flavonols, especially in Cabernet Sauvignon. The polyphenolic profile of this variety turned out to be qualitatively superimposable by confronting the relative abundance of quercetin with that of (+)-catechin. Quercetin was more concentrated than (+)-catechin in Sicily than Molise. Furthermore, flavonols concentration was higher in Sicily than in Molise. It might not be surprising since flavonols act in grapes as photo-protectors (Price et al., 1995; Gambuti et al., 2019; Gambuti et al., 2020), and in Sicily, plant tissues were more exposed to sunlight than in Molise. To understand whether this effect detected on grapes is determinant also for wine quality, the corresponding wines were obtained from each grape experimental sample with the same winemaking protocol and evaluated base parameters together with the content of total phenolics, flavonols, and anthocyanins. In regard to the Aglianico, unlike the Cabernet ones, our NMR-based investigation did not bring to light any differences among the identified polyphenols from either qualitative or quantitative standpoints. As already observed in Aglianico grapes (Forino et al., 2019), the content of (+)-catechin (83.1 ± 1.2 mg/L) and (-)-epicatechin (103.4 ± 1.8 mg/L) was comparable, whilst quercetin was not detected under our NMR experimental conditions. In contrast, our results on berries in terms of titratable acidity and pH showed that Aglianico is more sensitive to the environment compared to Cabernet. This datum is relevant and should be confirmed in further years given the relevance of acidic profile for microbial contamination, oxidation risk, and sensory quality of wines. Differences observed in grapes for total flavanols of Cabernet between Molise and Sicily were not confirmed in wines, probably for the extractability of these compounds in winemaking and to the numerous reactions of polymerization and precipitation in which flavanols are involved and that occur just after grape crushing until wine bottling (Forino et al., 2020). Regarding the role of variety on phenolic composition, the amount of flavanol in Aglianico wines was higher than in Cabernet Sauvignon, as already detected for these grape cultivars (Gambuti et al., 2006). For both cultivars, passing from the higher latitude (Molise) to the lowest (Sicily), their content decreases. It is known that the biosynthesis of anthocyanins depends on UV radiation and thermal excursion. The higher the radiation and the thermal excursion, the higher the synthesis of anthocyanins. In contrast, high night temperatures determine the degradation of anthocyanins during berry ripening (Teixeira et al., 2013). Therefore, Sicily’s higher night temperature and lower thermal excursion can easily explain the significant difference between Molise and Sicily for both cultivars. In this sense, the most sensitive cultivar to latitude was the Aglianico, which doubled the total monomeric anthocyanins passing from Sicily to Molise. In contrast, previous studies on Pinot Noir grapevines sampled across a latitudinal gradient in Europe, from southern Spain to central Germany showed that total anthocyanins did not show any correlation with radiation variables (Del-Castillo-Alonso et al., 2016). This confirms the paramount importance of cultivars. The data on anthocyanins in wines also confirmed this dependence on latitude. We observed higher levels of anthocyanins in Cabernet Sauvignon compared to Aglianico. Specifically, when comparing the individual environments, grapes of both varieties cultivated in Molise showed a richer anthocyanin pigment content compared to the corresponding varieties grown in Sicily. This finding aligns with the expression analysis of biosynthetic genes, revealing higher expression levels of F3H, LDOX, GST2, GST3, and GGT in Molise. In Sicily, we recorded reduced anthocyanin accumulation along with low or insignificant levels of biosynthetic gene expression (DFR, UFGT) and a higher degradation rate. This could be attributed to the climatic conditions in Sicily, where heat and drought stressors coincide. Previous studies have shown that drought conditions up-regulate anthocyanin biosynthesis (Castellarin et al., 2007), while heat stress negatively affects anthocyanin concentrations in grapes (Mori et al., 2007; Movahed et al., 2016; Pastore et al., 2017). Water-restricted cultivation has been suggested as a potential countermeasure to mitigate the negative impact of excessive heat on anthocyanin levels, particularly during later growth stages like berry ripening (Scholasch and Rienth, 2019). Recent findings by Tan et al. (2023) also support this, showing the downregulation of anthocyanin biosynthesis genes during combined drought and heat stress in Cabernet Sauvignon. Among the other examined classes, OMT1 and OMT3 have been found to exhibit a significant correlation with the methoxypyrazine content (Zhang et al., 2022). Our study revealed the overexpression of OMT1 in Sicily and varying levels between Aglianico and Cabernet Sauvignon, supporting the hypothesis that methoxypyrazine accumulation in grapes is influenced by light exposure and genotype (Koch et al., 2010; Sanders et al., 2023; Villano et al., 2023). The observed variations in grape characteristics, including soluble solids content, acidity, and pH, highlight the intricate relationship between environmental conditions and grape quality. These factors not only impact the composition of the grapes but also have implications for their fermentation capability and the ultimate sensory attributes of the resulting wines. Indeed, the effect of environmental conditions on grapes’ fermentation capability can be expressed by the ratio between soluble solids content and titratable acidity, which helps to understand if grapes can give problems during fermentation and if they can give equilibrated wines from the gustative point of view. Values too high of this ratio indicate that the fermentation of these grapes is not easy, and several problems such as sluggish or arrest of fermentation, could happen, as well as the appearance of microbial contamination. In addition, the obtained wines cannot be equilibrated and show several off-flavours. In our case studies, this ratio next to 4, well confirms the strong link between grape quality for fermentation and environmental conditions, especially for the Aglianico grapevine. In addition to chemical characteristics of grapes, our results, obtained using NGS method, reinforce the concept that also the grape yeast community composition differs according to geographical location and this biodiversity is a factor that could contribute to regional wine characteristics (Bokulich et al., 2016). In recent years, the contribution of the “native” (non-inoculated) yeasts to improving regional typicality has been enhanced (Liu et al., 2019). In fact, the specific grape-associated microbiota can exert both positive and negative effects on the quality of wine before, during and after fermentation (Spano and Torriani, 2016). This connection between microbial biogeography and regional wine characteristics has been termed “microbial terroir”, a term that hypothesizes a relationship between grapevine microbiology and wine terroir (Gilbert et al., 2014; Morgan et al., 2017).
5 Conclusion
In conclusion, the phenotypic response of Aglianico and Cabernet Sauvignon grapevines, encompassing vigor, production, and fruit quality, is significantly influenced by the prevailing pedoclimatic conditions, particularly soil physical properties. These conditions play a crucial role in shaping grape acidity, notably through the accumulation of organic acids, while their impact on pH expression appears to be comparatively less significant. Furthermore, pedoclimatic factors influence the content of total monomeric anthocyanins in grapes and wines, reflecting their effect on the expression of grapevine genes involved in phenolic and aroma compound metabolism, as well as qualitative traits. The findings presented in this study strongly support the concept of terroir and provide valuable insights into the current adaptation of vineyards and terroirs to climate change. Moreover, these results underscore the importance of comprehending the genetic plasticity of vine cultivars in navigating such environmental shifts.
Data availability statement
The original contributions presented in the study are publicly available. This data can be found here: https://data.mendeley.com/datasets/tp2rrdz8wp/2.
Author contributions
RA, AP, AB and RC conceived and designed research. MI, AS, EN, MF, LP, AV, EM, FF, MS, PT, CV and AB conducted experiments. MI, AS, EN, MF, LP, RA, AP, AV, EM, FF, MS, PT, AG, CV and AB analyzed data. MI, AS, MF, RA, AP, FF, AG, CV and AB wrote the manuscript with comments and input from all authors. RA, AG, CV, AB and RC revised the manuscript. All authors contributed to the article and approved the submitted version.
Acknowledgments
Work in our labs is funded by the Italian Ministry of University and Research, Research Projects of National Interest (PRIN) ADAPT - influence of agro-climatic conditions on the microbiome and genetic expression of grapevines for the production of red wines: a multidisciplinary approach (2017M83XFJ - CUP H34I19000590001). We thank Dr. Arturo Erbaggio for his support of grapevine monitoring, Dr. Roberto De Mascellis and Dr. Maurizio Buoananno for their support of the soil survey and environmental monitoring and Dr. Felice Napolitano for his laboratory assistance.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2023.1250208/full#supplementary-material
References
Amerine, M., Winkler, A. (1944). Composition and quality of musts and wines of California grapes. Hilgardia 15 (6), 493–675. doi: 10.3733/hilg.v15n06p493
Anesi, A., Stocchero, M., dal Santo, S., Commisso, M., Zenoni, S., Ceoldo, S., et al. (2015). Towards a scientific interpretation of the terroir concept: Plasticity of the grape berry metabolome. BMC Plant Biol. 15, 191. doi: 10.1186/s12870-015-0584-4
Basile, A., Buttafuoco, G., Mele, G., Tedeschi, A. (2012). Complementary techniques to assess physical properties of a fine soil irrigated with saline water. Environ. Earth Sci. 66 (7), 1797–1807. doi: 10.1186/s12870-015-0584-4
Battilana, J., Dunlevy, J. D., Boss, P. K. (2017). Histone modifications at the grapevine VvOMT3 locus, which encodes an enzyme responsible for methoxypyrazine production in the berry. Funct. Plant Biol. 44 (7), 655–664. doi: 10.1071/FP16434
Bokulich, N. A., Collins, T. S., Masarweh, C., Allen, G., Heymann, H., Ebeler, S. E., et al. (2016). Associations among wine grape microbiome, metabolome, and fermentation behavior suggest microbial contribution to regional wine characteristics. mBio 7, e00631–e00616. doi: 10.1128/mbio.00631-16
Bolyen, E., Rideout, J. R., Dillon, M. R., Bokulich, N. A., Abnet, C. C., Al-Ghalith, G. A., et al. (2019). Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 37 (8), 852–857. doi: 10.1038/s41587-019-0209-9
Bonfante, A., Alfieri, S. M., Albrizio, R., Basile, A., De Mascellis, R., Gambuti, A., et al. (2017). Evaluation of the effects of future climate change on grape quality through a physically based model application: a case study for the Aglianico grapevine in Campania region, Italy. Agric. Syst. 152, 100–109. doi: 10.1016/j.agsy.2016.12.009
Bonfante, A., Basile, A., Acutis, M., De Mascellis, R., Manna, P., Perego, A., et al. (2010). SWAP, CropSyst and MACRO comparison in two contrasting soils cropped with maize in Northern Italy. Agric. Water Manage. 97 (7), 1051–1062. doi: 10.1016/j.agwat.2010.02.010
Bonfante, A., Brillante, L. (2022). Terroir analysis and its complexity. OENO One 56, 375–388. doi: 10.20870/OENO-ONE.2022.56.2.5448
Bonfante, A., Monaco, E., Langella, G., Mercogliano, P., Bucchignani, E., Manna, P., et al. (2018). A dynamic viticultural zoning to explore the resilience of terroir concept under climate change. Sci. Total Environ. 624, 294–308. doi: 10.1016/j.scitotenv.2017.12.035
Brillante, L., Martínez-Luscher, J., Yu, R., Plank, C. M., Sanchez, L., Bates, T. L., et al. (2017). Assessing spatial variability of grape skin flavonoids at the vineyard scale based on plant water status mapping. J. Agric. Food Chem. 65, 5255–5265. doi: 10.1021/acs.jafc.7b01749
Brillante, L., Mathieu, O., Lévêque, J., van Leeuwen, C., Bois, B. (2018). Water status and must composition in grapevine cv. Chardonnay with different soils and topography and a mini meta-analysis of the δ13C/water potentials correlation. J. Sci. Food Agric. 98, 691–697. doi: 10.1002/jsfa.8516
Brillante, L., Martínez-Lüscher, J., Yu, R., Kurtural, S. K. (2020). Carbon isotope discrimination (δ13 C) of grape musts is a reliable tool for zoning and the physiological ground-truthing of sensor maps in precision viticulture. Front. Environ. Sci. 8 , 561477.
Bucchignani, E., Montesarchio, M., Zollo, A. L., Mercogliano, P. (2015). High-resolution climate simulations with COSMO-CLM over Italy: Performance evaluation and climate projections for the 21st century. Int. J. Climatol. 36, 735–756. doi: 10.1002/joc.4379
Callahan, B. J., McMurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J. A., Holmes, S. P. (2016). DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 13 (7), 581–583. doi: 10.1038/nmeth.3869
Castellarin, S. D., Pfeiffer, A., Sivilotti, P., Degan, M., Peterlunger, E., Di Gaspero, G. (2007). Transcriptional regulation of anthocyanin biosynthesis in ripening fruits of grapevine under seasonal water deficit. Plant Cell Environ. 30 (11), 1381–1399. doi: 10.1111/j.1365-3040.2007.01716.x
Clingeleffer, P. R. (2010). Plant management research: status and what it can offer to address challenges and limitations. Aust. J. Grape Wine Res. 16, 25–32. doi: 10.1111/j.1755-0238.2009.00075.x
Colombo, C., Miano, T. (2015). editor. Metodi di analisi chimicadel suolo – 3 edizione. Pubblicità & Stampa, (Modugno (BA): Italy.
Dal Santo, S., Zenoni, S., Sandri, M., De Lorenzis, G., Magris, G., De Paoli, E., et al. (2018). Grapevine field experiments reveal the contribution of genotype, the influence of environment and the effect of their interaction (G×E) on the berry transcriptome. Plant J. 93 (6), 1143–1159. doi: 10.1111/tpj.13834
Del-Castillo-Alonso, M. Á., Castagna, A., Csepregi, K., Hideg, E., Jakab, G., Jansen, M. A. K., et al. (2016). Environmental factors correlated with the metabolite profile of Vitis vinifera cv. Pinot Noir berry skins along a European latitudinal gradient. J. Agric. Food Chem. 64 (46), 8722–8734. doi: 10.1021/acs.jafc.6b03272
Dinu, D. G., Ricciardi, V., Demarco, C., Zingarofalo, G., De Lorenzis, G., Buccolieri, R., et al. (2021). Climate change impacts on plant phenology: Grapevine (Vitis vinifera) bud break in wintertime in southern Italy. Foods. 10, 2769. doi: 10.3390/FOODS10112769/S1
Di Vaio, C., Villano, C., Lisanti, M. T., Marallo, N., Cirillo, A., Di Lorenzo, R., et al. (2020). Application of anti-transpIrant to control sugar accumulation in grape berries and alcohol degree in wines obtained from thinned and unthinned vines of cv. Falanghina (Vitis vinifera L.). Agronomy 10 (3), 345. doi: 10.3390/agronomy10030345
Duchêne, E. (2016). How can grapevine genetics contribute to the adaptation to climate change? OENO One 50, 113–124. doi: 10.20870/oeno-one.2016.50.3.98
Dunlevy, J. D., Dennis, E. G., Soole, K. L., Perkins, M. V., Davies, C. (2010). Two O-methyltransferases involved in the biosynthesis of methoxypyrazines: grape-derived aroma compounds important to wine flavour. Plant Mol. Biol. 74 (1), 77–89. doi: 10.1007/s11103-013-0012-9
Forino, M., Gambuti, A., Moio, L. (2019). NMR-based systematic analysis of bioactive phytochemicals in red wine. First determination of xanthurenic and oleanic acids. Food Chem. 278, 497–501. doi: 10.1016/j.foodchem.2018.11.103
Forino, M., Picariello, L., Lopatriello, A., Moio, L., Gambuti, A. (2020). New insights into the chemical bases of wine color evolution and stability: The key role of acetaldehyde. Eur. Food Res. Technol. 246 (4), 733–743. doi: 10.1007/s00217-020-03442-x
Gambuti, A., Rinaldi, A., Pessina, R., Moio, L. (2006). Evaluation of aglianico grape skin and seed polyphenol astringency by SDS–PAGE electrophoresis of salivary proteins after the binding reaction. Food Chemistry 97 (4), 614–620.
Gambuti, A., Han, G., Peterson, A. L., Waterhouse, A. L. (2015). Sulfur dioxide and glutathione alter the outcome of microoxygenation. Am. J. Enol Vitic. 66 (4), 411–423. doi: 10.5344/ajev.2015.15005
Gambuti, A., Picariello, L., Moio, L., Waterhouse, A. L. (2019). Cabernet sauvignon aging stability altered by microoxygenation. Am. J. Enol Vitic. 70 (3), 323–331. doi: 10.5344/ajev.2019.18061
Gambuti, A., Picariello, L., Rinaldi, A., Forino, M., Blaiotta, G., Moine, V., et al. (2020). New insights into the formation of precipitates of quercetin in Sangiovese wines. J. Food Sci. Technol. 57, 2602–2611. doi: 10.1007/s13197-020-04296-7
Gilbert, J. A., van der Lelie, D., ZarraonaIndia, I. (2014). Microbial terroir for wine grapes. Proc. Natl. Acad. Sci. 111 (1), 5–6. doi: 10.1073/pnas.1320471110
Goode, J. (2012). Viticulture: Fruity with a hint of drought. Nature. 492, 351–353. doi: 10.1038/492351a
Harbertson, J. F. (2003). The measurement of tannins and polymeric pigments during grape ripening and winemaking using a combined bisulfite bleaching protein precipitation assay (Davis: University of California).
Intrigliolo, D. S., Perez, D., Risco, D., Yeves, A., Castel, J. R. (2012). Yield components and grape composition responses to seasonal water deficits in Tempranillo grapevines. Irrig. Sci. 30, 339–349. doi: 10.1007/s00271-012-0354-0
IPCC (2022). Climate change 2022: mitigation of climate change. Contribution of working group III to the sixth assessment report of the intergovernmental panel on climate change (Cambridge, UK and New York, NY, USA: Cambridge University Press). doi: 10.1017/9781009157926
IUSS Working Group WRB (2014). World Reference Base for Soil Resources 2014. International soil classification system for naming soils and creating legends for soil maps. World Soil Resources Reports No. 106 (Rome: FAO).
Ji, X., Wang, B., Wang, X., Wang, X., Liu, F., Wang, H. (2021). Differences of aroma development and metabolic pathway gene expression between Kyoho and 87-1 grapes. J. Integr. Agric. 20 (6), 1525–1539. doi: 10.1016/S2095-3119(20)63481-5
Jones, G. V., Webb, L. B. (2010). Climate change, viticulture, and wine: Challenges and opportunities. J. Wine Res. 21, 103–106. doi: 10.1080/09571264.2010.530091
Keller, M. (2010). Managing grapevines to optimise fruit development in a challenging environment: a climate change primer for viticulturists. Aust. J. Grape Wine Res. 16, 56–69.
Kobayashi, S., Takase, H., Suzuki, Y., Tanzawa, F., Takata, R., Fujita, K., et al. (2011). Environmental stress enhances biosynthesis of flavor precursors, S-3-(hexan-1-ol)-glutathione and S-3-(hexan-1-ol)-L-cysteine, in grapevine through glutathione S-transferase activation. J. Exp. Bot. 62 (3), 1325–1336. doi: 10.1093/jxb/erq376
Kobayashi, S., Ishimaru, M., Hiraoka, K., Honda, C. (2002). Myb-related genes of the Kyoho grape (Vitis labruscana) regulate anthocyanin biosynthesis. Planta 215, 924–933. doi: 10.1007/s00425-002-0830-5
Koch, A., Doyle, C. L., Matthews, M. A., Williams, L. E., Ebeler, S. E. (2010). 2-Methoxy-3-isobutylpyrazine in grape berries and its dependence on genotype. Phytochem. 71 (17-18), 2190–2198. doi: 10.1016/j.phytochem.2010.09.006
Lana, G., Modica, G., Casas, G. L., Distefano, G., Continella, A., Gentile, S., et al. (2021). Molecular insights into the effects of rootstocks on maturation of blood oranges. Horticulturae 7 (11), 468. doi: 10.3390/horticulturae7110468
Liu, D., Zhang, P., Chen, D., Howell, K. (2019). From the Vineyard to the Winery: How microbial ecology drives regional distinctiveness of wine. Front. Microbiol. 10. doi: 10.3389/fmicb.2019.02679
Livak, K. J., Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2^–ΔΔCT method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262
Lorenz, D. H., Eichhorn, K. W., Bleiholder, H., Klose, R., Meier, U., Weber, E. (1995). Growth stages of the grapevine: phenological growth stages of the grapevine (Vitis vinifera L. ssp. vinifera) - Codes and descriptions according to the extended BBCH scale. Aust. J. Grape Wine Res. 1 (2), 100–103. doi: 10.1111/j.1755-0238.1995.tb00085.x
Maritim, T. K., Masand, M., Seth, R., Sharma, R. K. (2021). Transcriptional analysis reveals key insights into seasonal induced anthocyanin degradation and leaf color transition in purple tea (Camellia sinensis (L.) O. Kuntze). Sci. Rep. 11. doi: 10.1038/s41598-020-80437-4
Martin, M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 17 (1), 10–12. doi: 10.14806/ej.17.1.200
Mattivi, F., Zulian, C., Nicolini, G., Valenti, L. (2002). Wine, biodiversity, technology, and antioxidants. Ann. N. Y. Acad. 957 (1), 37–56. doi: 10.1111/j.1749-6632.2002.tb02904.x
Monteith, J. L. (1965). “Evaporation and environment,” in Symposia of the society for experimental biology, vol. 19. (Cambridge: Cambridge University Press (CUP), 205–234.
Morgan, H. H., Du Toit, M., Setati, M. E. (2017). The grapevine and wine microbiome: insights from high-throughput amplicon sequencing. Front. Microbiol. 8. doi: 10.3389/fmicb.2017.00820
Mori, K., Goto-Yamamoto, N., Kitayama, M., Hashizume, K. (2007). Loss of anthocyanins in red-wine grape under high temperature. J. Exp. Bot. 58 (8), 1935–1945. doi: 10.1093/jxb/erm055
Movahed, N., Pastore, C., Cellini, A., Allegro, G., Valentini, G., Zenoni, S., et al. (2016). The grapevine VviPrx31 peroxidase as a candidate gene involved in anthocyanin degradation in ripening berries under high temperature. J. Plant Res. 129, 513–526. doi: 10.1007/s10265-016-0786-3
Nicolosi, E., Sicilia, A., Ferlito, F., Bonfante, A., Monaco, E., Lo Piero, A. R. (2022). Phenotypic plasticity in bud fruitfulness expressed in two distinct wine grape cultivars grown under three different pedoclimatic conditions. Agric. 12 (10), 1660. doi: 10.3390/agriculture12101660
Olivo, N., Girona, J., Marsal, J. (2009). Seasonal sensitivity of stem water potential to vapor pressure deficit in grapevine. Irrig. Sci. 27, 175–182. doi: 10.1007/s00271-008-0134-z
Paim Pinto, D. L., Brancadoro, L., dal Santo, S., de Lorenzis, G., Pezzotti, M., Meyers, B. C., et al. (2016). The influence of genotype and environment on small RNA Profiles in grapevine berry. Front. Plant Sci. 7. doi: 10.3389/fpls.2016.01459
Pastore, C., Allegro, G., Valentini, G., Muzzi, E., Filippetti, I. (2017). Anthocyanin and flavonol composition response to veraison leaf removal on Cabernet Sauvignon, Nero d’Avola, Raboso Piave and Sangiovese Vitis vinifera L. cultivars. Sci. Hortic. 218, 147–155. doi: 10.1016/j.scienta.2017.01.048
Perin, C., Fait, A., Palumbo, F., Lucchin, M., Vannozzi, A. (2020). The effect of soil on the biochemical plasticity of berry skin in two italian grapevine (V. vinifera L.) cultivars. Front. Plant Sci. 11. doi: 10.3389/fpls.2020.00822
Price, S. F., Breen, P. J., Valladao, M., Watson, B. T. (1995). Cluster sun exposure and quercetin in Pinot noir grapes and wine. Am. J. Enol. Vitic. 46 (2), 187–194. doi: 10.5344/ajev.1995.46.2.187
Rodrigues, R., Gentilesco, G., Coletta, A., Tarricone, L., Alba, V. (2023). Bioclimatic characterization relating to temperature and subsequent future scenarios of vine growing across the Apulia region in southern Italy. Agric. 13, 644. doi: 10.3390/agriculture13030644
Sanders, R. D., Boss, P. K., Capone, D. L., Kidman, C. M., Maffei, S., Jeffery, D. W. (2023). Methoxypyrazine concentrations in the grape bunch rachis of Vitis vinifera L. Cv Shiraz: Influence of rootstock, region and light. Food Chem. 408, 135234. doi: 10.1016/j.foodchem.2022.135234
Scholander, P. F., Hammel, H. T., Bradstreet, E. D., Hemmingsen, E. A. (1965). Sap pressure in vascular plants. Science 148, 339–346. doi: 10.1126/science.148.3668.339
Scholasch, T., Rienth, M. (2019). Review of water deficit mediated changes in vine and berry physiology; consequences for the optimization of irrigation strategies. Oeno One 53 (3). doi: 10.20870/oeno-one.2019.53.3.2407
Selvaraj, Y., Pal, D. K., Singh, R., Roy, T. K. (1994). Biochemistry of uneven ripening in gulabi grape 1. J. Food Biochem. 18 (5), 325–340. doi: 10.1111/j.1745-4514.1994.tb00507.x
Sicilia, A., Scialò, E., Puglisi, I., Lo Piero, A. R. (2020). Anthocyanin biosynthesis and dna methylation dynamics in sweet orange fruit [Citrus sinensis L. (Osbeck)] under Cold Stress. J. Agric. Food Chem. 68 (26), 7024–7031. doi: 10.1021/acs.jafc.0c02360
Spano, G., Torriani, S. (2016). Editorial: Microbiota of grapes: positive and negative role on wine quality. Front. Microbiol. 7. doi: 10.3389/fmicb.2016.02036
Sparvoli, F., Martin, C., Scienza, A., Gavazzi, G., Tonelli, C. (1994). Cloning and molecular analysis of structural genes involved in flavonoid and stilbene biosynthesis in grape (Vitis vinifera L.). Plant Mol. Biol. 24, 743–755. doi: 10.1007/BF00029856
Tan, J. W., Shinde, H., Tesfamicael, K., Hu, Y., Fruzangohar, M., Tricker, P., et al. (2023). Global transcriptome and gene co-expression network analyses reveal regulatory and non-additive effects of drought and heat stress in grapevine. Front. Plant Sci. 14. doi: 10.3389/fpls.2023.1096225
Teixeira, A., Eiras-Dias, J., Castellarin, S. D., Gerós, H. (2013). Berry phenolics of grapevine under challenging environments. Int. J. Mol. Sci. 14 (9), 18711–18739. doi: 10.3390/ijms140918711
Vannozzi, A., Perin, C., Palumbo, F., Sandri, M., Zuccolotto, P., Zenoni, S., et al. (2023). Dissecting the effect of soil on plant phenology and berry transcriptional plasticity in two Italian grapevine varieties (Vitis vinifera L.). Hortic. Res. 10 (05), uhad056. doi: 10.1093/hr/uhad056
Villano, C., Demurtas, O. C., Esposito, S., Granell, A., Rambla, J. L., Piombino, P., et al. (2023). Integrative analysis of metabolome and transcriptome profiles to highlight aroma determinants in Aglianico and Falanghina grape berries. BMC Plant Biol. 23 (1), 1–15. doi: 10.1186/s12870-023-04251-6
Villano, C., Rinaldi, A., Lanzillo, C., Moio, L., Tamburrino, A. J., Carputo, D., et al. (2016). Polyphenol content and differential expression of flavonoid biosynthetic pathway genes in berries of Aglianico. X Int. Symposium Grapevine Physiol. Biotechnol. 1188, 141–148. doi: 10.17660/ActaHortic.2017.1188.19
Wang, W., Feng, J., Wei, L., Khalil-Ur-Rehman, M., Nieuwenhuizen, N. J., Yang, L., et al. (2021). Transcriptomics integrated with free and bound terpenoid aroma profiling during “shine muscat” (Vitis labrusca× V. vinifera) grape berry development reveals coordinate regulation of MEP pathway and terpene synthase gene expression. J. Agric. Food Chem. 69 (4), 1413–1429. doi: 10.1021/acs.jafc.0c06591
Wolkowich, E. M., García de Cortázar-Atauri, I., Morales-Castilla, I., Nicholas, K. A., Lacombe, T. (2018). From Pinot to Xinomavro in the world's future wine-growing regions. Nat. Clim. Change. 8, 29–37. doi: 10.1038/s41558-017-0016-6
Xie, H., Konate, M., Sai, N., Tesfamicael, K. G., Cavagnaro, T., Gilliham, M., et al. (2017). Global dna methylation patterns can play a role in defining terroir in grapevine (Vitis vinifera cv. Shiraz). Front. Plant Sci. 8. doi: 10.3389/fpls.2017.01860
Zhang, Y., Li, X., Guo, X., Wang, N., Geng, K., Li, D., et al. (2022). Comparison of methoxypyrazine content and expression pattern of o-methyltransferase genes in grape berries and wines from six cultivars (Vitis vinifera L.) in the Eastern Foothill of the Helan Mountain. Plants. 11 (12), 1613. doi: 10.3390/plants11121613
Keywords: phenotypic plasticity, grape, soil, metabolites, microbiome, gene expression
Citation: Iorizzo M, Sicilia A, Nicolosi E, Forino M, Picariello L, Lo Piero AR, Vitale A, Monaco E, Ferlito F, Succi M, Tremonte P, Gambuti A, Villano C, Bonfante A, Aversano R and Coppola R (2023) Investigating the impact of pedoclimatic conditions on the oenological performance of two red cultivars grown throughout southern Italy. Front. Plant Sci. 14:1250208. doi: 10.3389/fpls.2023.1250208
Received: 29 June 2023; Accepted: 28 August 2023;
Published: 15 September 2023.
Edited by:
Jesus V. Jorrin Novo, University of Cordoba, SpainReviewed by:
Cinzia Montemurro, University of Bari Aldo Moro, ItalyAnnalisa Marchese, University of Palermo, Italy
Copyright © 2023 Iorizzo, Sicilia, Nicolosi, Forino, Picariello, Lo Piero, Vitale, Monaco, Ferlito, Succi, Tremonte, Gambuti, Villano, Bonfante, Aversano and Coppola. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Clizia Villano, Y2xpemlhLnZpbGxhbm9AdW5pbmEuaXQ=; Antonello Bonfante, YW50b25lbGxvLmJvbmZhbnRlQGNuci5pdA==
†These authors have contributed equally to this work and share first authorship