Literature DB >> 31013598

Effects of Two Different Irrigation Systems on the Amino Acid Concentrations, Volatile Composition and Sensory Profiles of Godello Musts and Wines.

José Manuel Mirás-Avalos1,2, Yolanda Bouzas-Cid3, Emiliano Trigo-Córdoba4, Ignacio Orriols5,6, Elena Falqué7.   

Abstract

The concentrations of amino acids and volatile compounds of a given grapevine cultivar may be modified by climate variability between years and by management practices, such as irrigation, that may alter the typicality of its wines. The current study aimed at assessing the amino acid profile of musts and wines, volatile composition and sensory profile of wines from Vitis vinifera (L.) cultivar Godello under rain-fed and two drip irrigation systems (above, drip irrigation (DI), and under the soil surface, subsurface drip irrigation (SDI)) over three consecutive years. Irrigation tended to increase must and wine total acidity; however, it did not alter must amino acid concentrations significantly. Irrigation reduced the concentrations of acetaldehyde and methanol in Godello wines. Moreover, irrigation tended to decrease the concentrations of compounds giving fruity aromas, such as acetaldehyde (by 31% in SDI) and isoamyl acetate (by 21% in SDI), when compared to rain-fed conditions. Sensory analysis revealed slight differences between treatments. Rain-fed and SDI were the treatments showing the greatest differences. Weather conditions affected more must and wine composition than in-season effects caused by irrigation.

Entities:  

Keywords:  aroma; global change; nitrogen compounds; vine water status; white wine

Year:  2019        PMID: 31013598      PMCID: PMC6517885          DOI: 10.3390/foods8040135

Source DB:  PubMed          Journal:  Foods        ISSN: 2304-8158


1. Introduction

Wine aroma is defined by its volatile composition, which depends on many factors including grape variety, climate, soil and vineyard management, amongst others [1,2]. Amino acids in grapes constitute a nitrogen source for yeasts and are responsible for the formation of volatile compounds that define wine aroma, including higher alcohols, volatile fatty acids and ethyl esters [3,4]. These compounds accumulate in grapes during ripening, although their concentrations depend on temperature and water availability [1,5]. Climate change alters the temporal distribution of rainfall and increases the intensity of drought events and heatwaves, raising a great concern in viticulture regions worldwide [6]. In order to counteract these negative effects, irrigation use in vineyards is increasing, even in cool-humid regions [7,8]. Irrigation management can be a tool for modulating berry growth and quality [9]. A great research effort has been devoted to assess the effect of water stress and irrigation scheduling on grapevine (Vitis vinifera L.) yield and grape and wine composition, mainly on red cultivars [9,10]. Recently, the effects of water stress on the amino acid and volatile composition of grapes and wines have been determined on white cultivars [11,12,13,14], in which aroma plays an essential role. However, these studies offered contrasting results depending on grapevine cultivar, intensity of water stress, irrigation strategy and climate conditions. In the Northwest of the Iberian Peninsula (Galicia and North of Portugal), white grapevine cultivars are predominantly grown and the volatile composition of their grapes and wines might be altered by increasing temperatures and drought events over the growing cycle. Amongst them, Godello provides monovarietal wines recognized worldwide. Volatile composition of Godello wines has already been described [15,16]. Recently, our research group assessed the effect that irrigation might exert on Godello amino acid profile in the Ribeiro Designation of Origin (DO) [17], although only one irrigation strategy was tested and its effects on the wine sensory profile have not been determined. Therefore, the aim of the current study was to assess the effects of two irrigation systems on the amino acid composition of musts and wines, and the aromatic and sensory profiles of wines from Godello during three years in Northwest Spain. This is a comprehensive approach for studying wine quality since it involves both chemical and sensory evaluations. The obtained results should be of interest to grape growers, winemakers and other wine professionals.

2. Materials and Methods

2.1. Studied Vineyard and Experimental Design

The study was carried out from 2012 to 2014 in a Godello vineyard planted in 1997 onto 110-Richter rootstock and located in A Rúa, Ourense, Spain (42° 23’ 59’’ N, 7° 7’ 15’’ W, elevation 320 m). Vines were spaced 1.35 × 1.95 m (3800 vines ha−1) and vertically trellised on a double cordon system. Rows were oriented in the north-south direction. Soil is loamy-textured, very acidic and available water capacity is 170 mm m−1. According to the multicriteria climatic classification system [18], this area is temperate, moderately dry with cool nights. The experiment consisted of three treatments on a randomized-block design with three replicates (three rows each, the middle row was used for measurements and sampling, whereas the external rows acted as buffers). The treatments were: Rain-fed (R), surface (DI) and subsurface (SDI) drip irrigation, which reduces direct evaporation from soil and, therefore, more water is available for the vines. The DI pipes were installed on the vineyard row 40 cm above the soil, while SDI pipes were buried 40 cm deep into the soil. Both systems were equipped with 2 L h−1 emitters [7], one emitter per vine in the case of DI and one emitter per meter in the case of SDI. Irrigation treatments started on June 1st and finished in mid-August in 2012; during 2013, irrigation began in July and ended by late August; however, due to the pumping system malfunctioning, irrigation started in mid-July and finished by the end of August. In summary, water was applied for 59, 46 and 34 days in 2012, 2013 and 2014, respectively, at a rate of 1.5 h per day. Average irrigation depths were 1.14 mm and 1.54 mm for DI and SDI, respectively [19]. Agricultural practices were the same for all treatments, except for irrigation. Several variables were measured over the study period for characterizing the water status, vegetative growth and yield of vines in each treatment. Both methods and results about these determinations are out of the scope of this report and have been described elsewhere [7,19]. However, we summarized the most important findings in the Results and Discussion section.

2.2. Sampling and Winemaking

Harvest were performed on the same date, each year, for the three treatments. Samples of 20 kg per field replication were collected for winemaking, one tank was used for two replications; therefore, there were two wines per treatment and year. Two 250 mL must samples from each replicate were collected for determining the general attributes of musts (soluble solids, total acidity, pH) according to official methods after grape pressing [20]. Must samples for determination of amino acids were stored at −4 °C until analysis. The winemaking process was the same as that reported previously [13]. In summary, grapes from each replicate were separately destemmed, crushed and pressed (yielding 50% must). Pectolytic enzyme was added (4 g per 100 kg of grape) to favor settling and SO2 (50 mg L−1) was added to avoid oxidation. After 24 h, musts were racked and, then, fermented in 35 L stainless steel tanks. Commercial yeast (Excellence FW, Lamothe-Abiet, Bordeaux, France) was added (20 g h L−1). Once fermentation finished, wines were racked and SO2 was added to 35 mg L−1 free sulfur dioxide. A natural clarification was performed at 4 °C for one month. Finally, wines were filtered, bottled and stored for five months at 10 °C until analysis. The main attributes of wines were determined according to the official methods [20].

2.3. Analytical Determinations

2.3.1. Chemical Reagents

Milli-Q equipment (Millipore, Bedford, MA, USA) allowed for obtaining ultra-pure water. Super-gradient High-Performance Liquid-Chromatography (HPLC) grade acetonitrile and methanol were from Scharlau (Sentmenat, Spain). Ammonium chloride was from Merck (Darmstadt, Germany). Individual L-α amino acids and diethylethoxymethylenemalonate were from Acros Organics (New Jersey, NJ, USA). Amino acid solutions were made with HCl 0.1 N. Dichloromethane, n-pentane and anhydrous sodium sulfate were from Scharlau (Sentmenat, Spain). Standards for volatile compounds were from Merck (Madrid, Spain), Aldrich (Madrid, Spain), Fluka (Seelze, Germany), Alfa Aesar (Barcelona, Spain) and Sigma (Madrid, Spain). All the standards were prepared in 50% hydroalcoholic solutions.

2.3.2. Determination of Free Amino Acids

The determination of amino acids in Godello musts and wines from each treatment was performed by HPLC according to previous works [13,21], using Agilent 1100 series equipment (Agilent Technologies, Palo Alto, CA, USA). Extraction method, reagents, elution conditions and chromatographic separation were performed as described in a previous report [13]. Determinations were carried out in triplicate.

2.3.3. Determination of Volatile Compounds

Concentrations of major volatile compounds were quantified by direct injection in a 7890A gas chromatograph (Agilent Technologies, Palo Alto, CA, USA) [12,16]. Terpenes, C6 alcohols, volatile fatty acids, ethyl esters of fatty acids and acetates of higher alcohols were extracted and injected into a 6890 gas chromatograph (Agilent Technologies, Palo Alto, CA, USA) using the conditions reported in a previous work [13]. Identification was performed using the National Institute of Standards and Technology (NIST) Mass Spectral library by comparing mass spectra and retention times with those of pure standard compounds. All determinations were carried out in triplicate. The ratio between concentration of an individual compound and its perception threshold was computed to obtain the odor activity value (OAV) of volatile compounds [22,23].

2.4. Sensory Evaluation

A panel of nine judges (five male, four female, age span 30–64) evaluated, each year, the wines from this study. These professionals were oenologists and technicians from Galician wineries, with expertise on assessing white wines, because experienced tasters describe better what they like than consumers or trained panelists [24,25]. For this evaluation, we used a scorecard consisting of 21 descriptors (four for color, 10 for aroma, and seven for palate), which were specifically chosen for Galician white wines and were scored from 0 (not present) to 9 (very intense) [26,27]. Apart from the fruity descriptors such as dry fruit or tropical fruit, quantitative descriptors such as persistence (duration of the perception of aroma) and intensity were included in this scorecard. Furthermore, tasters scored the global quality of the wine. Wine samples (30 mL) coded with three random numbers were presented in clear tulip-shaped glasses. The tasting sessions took place in different days (April and May from the year after the corresponding vintage).

2.5. Statistical Analysis

The effects of year, irrigation treatment and their interaction on amino acids and volatile compound concentrations were assessed by analysis of variance (ANOVA). When required, mean separation was performed using the Tukey’s test. A principal component analysis (PCA) was carried out to separate must samples according to their amino acid concentrations and wines according to their volatile and sensory profiles. Statistical analyses were performed using R v3.2.2 [28]. Data from the sensory evaluation were processed using Big Sensory Soft 1.02 (Centro Studi Assagiatori, Brescia, Italy), using the non-parametric Friedman test to discern which descriptors differed between treatments since these data do not comply with the assumption of normality [29].

3. Results and Discussion

3.1. Weather Conditions and Vine Water Status

A summary of the weather conditions over the study period (2012–2014) is presented in Table 1. Annual rainfall was lower in 2012 than in the other studied years. However, rainfall over the growing season (April to harvest) was 260 mm, 331 mm and 239 mm in 2012, 2013 and 2014, respectively (Table 1). During maturation (veraison to harvest), temperatures were similar among the studied years; however, the rainfall pattern differed. Annual mean temperature was higher in 2014 than in 2012 and 2013. Mean temperatures over the growing season increased from year to year, and reference evapotranspiration (ETo) was greater in 2013 than in the other studied years (Table 1).
Table 1

Mean temperature, potential evapotranspiration (ETo) and rainfall over the growing season (April to harvest), and annual rainfall and mean temperature of each year (2012–2014) in the studied vineyard.

YearGrowing Season Rainfall(mm)Annual Rainfall(mm)Growing Season Mean Temperature(°C)Annual Mean Temperature(°C)Growing Season ETo(mm)
201226054316.812.7706
201333192617.112.7741
201423982517.413.4698
Previous studies reported data on the effects of irrigation on physiology and yield response of vines in this site [7,19]. In summary, rain-fed vines showed more negative midday stem water potentials than those irrigated; whereas yield and cluster weight were similar among treatments (Table 2). Mild water stress was observed in both rain-fed and irrigated vines in 2012; while rain-fed vines suffered from mild water stress also in 2013 (Table 2). No stress conditions were observed in 2014. In summary, rain-fed vines showed more negative stem water potentials than those irrigated, especially from veraison onwards; these differences were more marked in 2013, the warmest growing season [19]. This explains the absence of significant differences among treatments in yield. In fact, only the number of clusters per vine differed significantly among treatments in 2013, when the vines from SDI had a greater number of clusters than those from the other treatments (Table 2).
Table 2

Water status, number of clusters, yield and cluster weight of Godello vines subjected to rain-fed and irrigation conditions during three consecutive seasons (2012, 2013, and 2014).

YearTreatmentMidday Stem Water Potential from Veraison to Harvest (MPa)Clusters per VineYield (kg vine−1)Cluster Weight (g cluster−1)
2012R−0.93 b22.22.8130.3
DI−0.85 ab19.42.9145.0
SDI−0.80 a21.63.5152.5
2013R−0.91 b21.3 a3.2142.7
DI−0.78 a20.8 a3.7171.0
SDI−0.71 a26.7 b4.2152.8
2014R−0.72 b22.92.9122.1
DI−0.62 a21.42.7125.0
SDI−0.58 a19.32.4117.7

For each year, different letters indicate significant differences between treatments at p < 0.05. R = rain-fed; DI = drip irrigation; SDI = subsurface drip irrigation.

3.2. General Parameters of Musts and Wines

Year exerted a significant effect on tartaric acid, whereas treatment did not affect any of the considered attributes, although a trend (p < 0.1) to higher acidities in musts from the irrigated treatments was observed (Table 3). No significant interactions between factors were detected. Acid-balanced wines usually have refreshing or crisp sensory undertones because organic acids, such as tartaric, play important roles in the development of specific flavor compounds [30].
Table 3

General attributes of the musts (mean ± standard error) of Godello vines under three irrigation regimes over three seasons (2012, 2013, and 2014). The significances of the factors year and treatment, as well as their interaction are shown.

Parameter201220132014TYT × Y
RDISDIRDISDIRDISDI
Total soluble solids (°Brix) 23.3 ± 0.122.8 ± 0.122.2 ± 0.424.3 ± 0.224.1 ± 0.624.0 ± 0.322.0 ± 0.021.7 ± 0.121.7 ± 0.2nsnsns
Total acidity (g L−1) 6.5 ± 0.77.0 ± 0.27.9 ± 0.66.2 ± 0.26.1 ± 0.16.9 ± 0.57.1 ± 0.67.4 ± 0.37.9 ± 0.2nsnsns
pH 3.20 ± 0.023.18 ± 0.013.14 ± 0.013.33 ± 0.013.33 ± 0.023.26 ± 0.063.17 ± 0.053.13 ± 0.043.06 ± 0.02nsnsns
Malic acid (g L−1) 4.0 ± 1.43.4 ± 0.24.1 ± 0.32.8 ± 0.32.9 ± 0.13.1 ± 0.13.0 ± 0.33.1 ± 0.23.2 ± 0.1nsnsns
Tartaric acid (g L−1) 5.7 ± 0.14.7 ± 0.86.3 ± 0.46.6 ± 0.06.3 ± 0.26.8 ± 0.48.2 ± 0.48.2 ± 0.38.6 ± 0.4ns***ns

For each year, different letters indicate significant differences between treatments at p < 0.05. The significance of the year, treatment and their interaction is expressed as ns = not significant; *** p < 0.001. R = rain-fed; DI = drip irrigation; SDI = subsurface drip irrigation; T = Treatment; Y = Year.

Vines in our study experienced a weak level of water stress (Table 2) and this prevents to observe changes in must attributes [31]. However, total acidity showed a great sensitivity to irrigation, as previously observed for other cultivars [8,9,14]. These slight changes might be not sufficient to detect differences among treatments in the analytical determinations carried out; however, a year-to-year variability existed in the general attributes of musts. For instance, up to 2 °Brix difference between grapes harvested in 2013 and 2014, accompanied by a 1 g L−1 difference in total acidity and about 1.5 g L−1 in tartaric acid concentration (Table 3). This different maturation can be explained by rainfall and temperature patterns over the ripening period in each year, exerting a great effect on musts and wine quality [1]. In wines, both year and treatment significantly altered total acidity and tartaric acid concentrations, following the trends detected in musts. Wines from SDI showed the highest TA and the lowest pH values (Table 4). This response is common in white cultivars under weak or mild stress conditions [7,8,32], and might have been caused by higher water availability over the maturation period in the irrigated treatments, especially in SDI.
Table 4

General attributes of the wines (mean ± standard error) from Godello vines under three irrigation regimes over three seasons (2012, 2013, and 2014). The significances of the factors year and treatment, as well as their interaction are shown.

Parameter201220132014TYT × Y
RDISDIRDISDIRDISDI
Alcohol (%vol.) 14.1 ± 0.114.0 ± 0.213.6 ± 0.514.0 ± 0.214.2 ± 0.114.1 ± 0.014.3 ± 0.214.3 ± 0.114.0 ± 0.0nsnsns
Total acidity (g L−1) 7.0 ± 0.1 ab 6.7 ± 0.1 a 7.6 ± 0.2 b 7.1 ± 0.2 a 7.5 ± 0.2 a 8.0 ± 0.3 b 7.9 ± 0.38.1 ± 0.28.7 ± 0.3******ns
pH 3.26 ± 0.023.26 ± 0.023.12 ± 0.063.37 ± 0.083.39 ± 0.013.26 ± 0.06 3.18 ± 0.01 b 3.09 ± 0.09 ab 2.89 ± 0.01 a nsnsns
Malic acid (g L−1) 2.4 ± 0.32.4 ± 0.22.7 ± 0.12.5 ± 0.32.9 ± 0.12.7 ± 0.12.3 ± 0.22.4 ± 0.12.2 ± 0.2nsnsns
Tartaric acid (g L−1) 2.2 ± 0.32.4 ± 0.12.8 ± 0.42.9 ± 0.32.7 ± 0.13.3 ± 0.3 4.2 ± 0.2 a 4.3 ± 0.4 a 5.6 ± 0.3 b ****ns

For each year, different letters indicate significant differences between treatments at p < 0.05. The significance of the year, treatment and their interaction is expressed as ns = not significant; * p < 0.05; *** p < 0.001. Bold letters indicate significant differences among treatments for a given parameter and year. R = rain-fed; DI = drip irrigation; SDI = subsurface drip irrigation; T = Treatment; Y = Year.

3.3. Amino Acids Profiles of Musts and Wines

The amino acid concentrations in Godello musts (Table 5) were within previously reported levels, except for proline that appeared at concentrations lower than those previously reported for other cultivars [3]. Proline is accumulated in grapes when vines suffer from water stress [12,33], which was not the case in our study (Table 2). The concentrations of amino acids observed in the current study were very similar to those reported for Godello in a previous work despite the different rootstock, soil and weather conditions between the experiments in both studies [17]. This agrees with the hypothesis that amino acid profile is relatively constant for a given variety [21].
Table 5

Irrigation effects on the amino acid concentrations (mean ± standard error, mg L−1) of Godello musts. The significances of the factors year and treatment, as well as their interaction are shown. R = rain-fed, DI = drip irrigation, SDI = subsurface drip irrigation, T = Treatment, Y = Year.

Compound201220132014TYT × Y
RDISDIRDISDIRDISDI
Aspartic acid 42.5 ± 8.544.5 ± 5.954.1 ± 0.151.7 ± 0.859.8 ± 3.158.7 ± 1.922.2 ± 3.824.7 ± 1.725.3 ± 2.1ns*ns
Glutamic acid 75.7 ± 10.060.3 ± 4.842.7 ± 1.995.7 ± 2.3101.9 ± 4.4100.4 ± 3.0102.8 ± 10.796.5 ± 2.989.6 ± 2.3ns**ns
Asparagine 3.2 ± 0.73.0 ± 1.32.8 ± 1.03.1 ± 0.63.4 ± 0.24.2 ± 0.51.3 ± 0.31.3 ± 0.01.7 ± 0.4ns*ns
Serine 44.4 ± 3.847.9 ± 7.540.8 ± 7.336.7 ± 3.739.5 ± 4.638.3 ± 0.526.6 ± 5.328.4 ± 0.926.7 ± 1.2ns***ns
Glutamine 61.3 ± 8.969.5 ± 24.555.1 ± 11.043.6 ± 6.643.7 ± 7.041.3 ± 2.354.9 ± 17.856.0 ± 1.549.6 ± 4.1nsnsns
Histidine 18.7 ± 3.317.5 ± 4.112.9 ± 2.710.0 ± 1.310.9 ± 1.610.0 ± 0.48.9 ± 1.99.5 ± 0.18.8 ± 1.0ns**ns
Glycine 2.5 ± 0.32.6 ± 0.32.0 ± 0.62.7 ± 0.32.8 ± 0.52.7 ± 0.02.6 ± 0.33.0 ± 0.02.8 ± 0.0nsnsns
Threonine 79.8 ± 13.185.6 ± 19.774.4 ± 13.755.3 ± 7.254.6 ± 7.554.4 ± 2.344.0 ± 11.546.7 ± 1.444.9 ± 3.3ns***ns
Arginine 153.9 ± 23.8245.6 ± 88.3148.7 ± 41.9120.4 ± 21.3127.9 ± 16.8126.5 ± 2.4157.7 ± 47.1170.0 ± 2.9159.6 ± 20.6nsnsns
Alanine 62.4 ± 15.978.6 ± 27.857.1 ± 16.551.0 ± 11.157.8 ± 12.856.5 ± 4.146.5 ± 13.451.9 ± 3.743.5 ± 3.4nsnsns
γ-Aminobutyric acid (GABA) 102.1 ± 4.3111.3 ± 4.782.8 ± 18.930.0 ± 3.732.8 ± 6.928.6 ± 1.726.6 ± 3.026.5 ± 1.220.6 ± 1.8ns***ns
Proline 4.9 ± 1.13.5 ± 1.92.2 ± 1.36.4 ± 0.06.1 ± 0.46.7 ± 1.03.2 ± 0.42.7 ± 0.22.7 ± 0.2nsnsns
Tyrosine 3.0 ± 0.33.2 ± 0.52.8 ± 0.12.5 ± 0.62.3 ± 0.31.9 ± 0.52.0 ± 0.22.4 ± 0.12.2 ± 0.1ns*ns
Ammonium ion 124.3 ± 19.4124.1 ± 21.9120.9 ± 1.1142.4 ± 3.4148.6 ± 10.9149.1 ± 13.2117.0 ± 23.0132.1 ± 6.4127.3 ± 9.4nsnsns
Valine 18.8 ± 1.517.9 ± 1.412.0 ± 2.213.6 ± 0.714.7 ± 2.713.4 ± 0.413.9 ± 1.414.8 ± 0.113.2 ± 0.6nsnsns
Methionine 1.6 ± 0.21.9 ± 0.61.2 ± 0.21.7 ± 0.21.5 ± 0.01.7 ± 0.01.2 ± 0.21.5 ± 0.21.3 ± 0.1nsnsns
Cysteine 2.5 ± 0.22.3 ± 0.22.6 ± 0.02.2 ± 0.12.3 ± 0.02.1 ± 0.11.9 ± 0.01.9 ± 0.11.9 ± 0.1ns***ns
Isoleucine 8.8 ± 1.07.9 ± 0.96.0 ± 1.05.7 ± 0.36.1 ± 1.15.8 ± 0.17.0 ± 0.87.4 ± 0.16.7 ± 0.4nsnsns
Tryptophan 5.7 ± 0.25.3 ± 2.35.5 ± 0.3 4.3 ± 0.2 b 3.7 ± 0.0 ab 2.7 ± 0.3 a 2.1 ± 0.32.2 ± 0.21.5 ± 0.1ns***ns
Leucine 10.3 ± 0.911.1 ± 1.56.9 ± 1.46.0 ± 0.67.4 ± 1.47.3 ± 0.49.2 ± 0.79.2 ± 0.08.8 ± 0.8nsnsns
Phenylalanine 8.0 ± 0.38.3 ± 1.86.1 ± 1.25.5 ± 0.75.8 ± 0.95.5 ± 0.17.1 ± 0.97.6 ± 0.07.1 ± 0.1nsnsns
Ornithine 1.1 ± 0.31.6 ± 0.71.0 ± 0.20.7 ± 0.20.7 ± 0.10.7 ± 0.00.5 ± 0.10.5 ± 0.00.5 ± 0.0ns**ns
Lysine 3.7 ± 0.63.9 ± 1.03.1 ± 0.72.6 ± 0.52.8 ± 0.42.8 ± 0.12.1 ± 0.32.2 ± 0.12.1 ± 0.2ns**ns
Sum of amino acids 714.9 ± 92.4833.3 ± 201.1622.6 ± 123.8551.1 ± 62.7588.4 ± 66.7572.0 ± 11.6544.2 ± 114.2566.5 ± 6.9520.7 ± 39.0ns*ns

Different letters indicate significant differences between treatments at p < 0.05 for a given amino acid in a given year; for a rapid identification, these values appear in bold. The significance of the year, treatment and their interaction is expressed as ns = non-significant; * p < 0.05; ** p < 0.01; *** p < 0.001.

In this sense, Godello is an arginine accumulator as this amino acid accounted for 25% of the total free amino acids in musts, similarly to other cultivars, such as Syrah, Petit Verdot and Merlot [21]. Other amino acids present at relevant concentrations in Godello musts were glutamic acid (14%), alanine (10%) and glutamine (9%); abundant in other Spanish white varieties, such as Verdejo and Albariño [11,34]. In our study, irrigation did not modify the amino acid profile of this variety over the three years studied. In the current study, irrigation did not affect the amino acid concentrations of Godello musts except for reducing that of tryptophan in 2013 (Table 5). Although not significantly, SDI tended to reduce glutamic acid (p = 0.07) and valine (p = 0.06) concentrations when compared to the rain-fed control in 2012. In contrast, the effect of year was significant for aspartic and glutamic acids, asparagine, serine, histidine, threonine, γ-aminobutyric acid (GABA), tyrosine, cysteine, tryptophan, ornithine and lysine, and the sum of amino acid concentrations in Godello musts (Table 5). No significant interactions between treatment and year were detected. As previously stated, rainfall and temperature regimes between veraison and harvest differed among the studied years, being 2013 a drier year than 2012 and 2014, likely causing significant variations in the concentrations of amino acids in Godello musts, as these compounds are strongly affected by weather conditions during berry maturation [1]. In contrast, the weak water deficit suffered by rain-fed vines between veraison and harvest compared to irrigated vines (Table 2) was not enough to alter must amino acid profiles. This was confirmed by a PCA ran on the covariance matrix of the concentrations of the free amino acids in must. The two first principal components (PC) explained 80.8% of the total variance in the dataset (Figure 1). Samples separated according to year, whereas irrigation treatment had a lower influence. Musts from 2012 were on the positive side of PC1, due to high concentrations of most amino acids, including asparagine, GABA, histidine, valine, alanine, etc. Musts from 2013 were on the positive side of PC2 but on the negative side of PC1 due to their high concentrations of proline. Samples from 2014 were on the negative side of both PC due to their high concentrations of glycine. In the case of 2012 samples, PCA showed a separation according to treatment, samples from R treatment had higher concentrations of GABA, histidine, tyrosine, alanine and ornithine, whereas samples from DI had high concentrations of valine and those from SDI had higher concentrations of asparagine and cysteine (Figure 1).
Figure 1

Principal component analysis (PCA) of Godello musts: Biplot for the first two principal components (PC) for the concentrations of free amino acids in musts. R = Rain-fed, DI = Drip irrigation, SDI = Subsurface drip irrigation.

In agreement with a previous report on the cultivar Verdejo [11], our results showed that irrigation did not affect amino acid concentrations in Godello musts, except for some individual compounds in a given year, being weather conditions and maturation stage the main drivers for determining the concentrations of amino acids in the musts. In our study, the highest amino acid concentrations in Godello musts were observed in 2012 [7], the year with the most negative midday stem water potentials (Table 2). Weather conditions can cause a great variation in the concentrations of individual amino acids from year to year, although their relative proportions over the total free amino acids remain constant [11,34]. Although previous research showed that the percentage of amino acids allowed for a differentiation of grapes according to the variety and cultivation system [21], regardless of the year, other studies revealed a great inter-annual variation on the amino acid composition of grapes and musts [11]. This is directly related to the weather conditions of each year. In the current study, mean temperature increased from year to year; however, evapotranspiration was higher in 2013. In contrast, rainfall distribution was different among years. In general, 2012 could be considered a cooler and drier year (Table 1). Likely, this has caused that samples were more different between years than between irrigation treatments (Figure 2).
Figure 2

Aroma descriptors for Godello wines (from 2012 to 2014) as a function of rain-fed and irrigation conditions. R = Rain-fed, DI = Drip irrigation, SDI = Subsurface drip irrigation. Asterisks indicate significant differences among treatments according to the Friedman’s test.

In the case of wines, those from SDI showed the lowest concentrations of serine, glutamine, glycine and isoleucine, as well as the lowest sum of free amino acids (Supplementary Table S1). Year exerted a significant influence on the concentration of 15 individual compounds and on the total concentration of free amino acids in Godello wines. No significant interactions between treatment and year were detected. These results were similar to those from a previous work [17].

3.4. Volatile Composition of Godello Wines

Volatile compounds in the wines from the current study (Table 6) were present at similar concentrations to those previously reported for Godello wines [15,16]. Irrigation tended to reduce the concentrations of methanol, acetaldehyde, acetoine, acetol and citronellol (Table 6). In addition, year affected most of the volatile compounds determined in Godello wines (Table 6). Significant interactions between treatment and year were observed for benzyl alcohol, acetol and hodiol I.
Table 6

Irrigation effects on the concentrations of volatile compounds (mean ± standard error) of Godello wines. The significances of the year, treatment as well as their interaction are also shown. R = rain-fed, DI = drip irrigation, SDI = subsurface drip irrigation, T = Treatment, Y = Year.

201220132014TYT × Y
RDISDIRDISDIRDISDI
Methanol (mg L−1)21 ± 118 ± 117 ± 019 ± 120 ± 118 ± 022 ± 319 ± 118 ± 1*nsns
Ethyl acetate (mg L−1)59 ± 1156 ± 559 ± 126 ± 226 ± 127 ± 337 ± 634 ± 138 ± 0ns*ns
Acetaldehyde (mg L−1)58 ± 1239 ± 128 ± 231 ± 335 ± 330 ± 231 ± 628 ± 025 ± 0***ns
Higher alcohols (mg L−1)
1-propanol12 ± 013 ± 112 ± 2 22 ± 1 ab 24 ± 0 b 20 ± 0 a 23 ± 424 ± 417 ± 2ns**ns
2-methyl-1-propanol61 ± 560 ± 369 ± 735 ± 237 ± 134 ± 028 ± 624 ± 345 ± 8ns***ns
2-methyl-1-butanol67 ± 464 ± 167 ± 164 ± 172 ± 164 ± 342 ± 444 ± 450 ± 1ns***ns
3-methyl-1-butanol252 ± 23261 ± 19275 ± 4 289 ± 6 a 313 ± 2 b 286 ± 0 a 221 ± 8219 ± 8253 ± 15nsnsns
∑ Higher alcohols392 ± 21398 ± 23423 ± 12 409 ± 8 a 446 ± 2 b 403 ± 3 a 314 ± 14311 ± 11365 ± 21ns**ns
Other alcohols (mg L−1)
1-hexanol1.8 ± 0.11.7 ± 0.82.3 ± 0.41.8 ± 0.01.8 ± 0.21.8 ± 0.02.0 ± 0.22.0 ± 0.12.1 ± 0.1nsnsns
trans-3-hexen-1-ol0.26 ± 0.060.34 ± 0.080.16 ± 0.020.14 ± 0.000.13 ± 0.000.12 ± 0.000.15 ± 0.010.15 ± 0.000.12 ± 0.03ns*ns
cis-3-hexen-1-ol0.54 ± 0.100.52 ± 0.010.27 ± 0.030.21 ± 0.020.22 ± 0.010.23 ± 0.030.22 ± 0.000.23 ± 0.010.18 ± 0.05ns**ns
Benzyl alcohol 3.82 ± 0.06 b 2.73 ± 0.03 a 2.77 ± 0.04 a 2.28 ± 0.142.47 ± 0.132.15 ± 0.11 2.12 ± 0.08 a 2.48 ± 0.03 ab 2.79 ± 0.11 b ns****
2-phenylethanol56 ± 250 ± 443 ± 548 ± 457 ± 150 ± 4 34 ± 1 ab 38 ± 0 b 33 ± 1 a ns**ns
Other compounds (mg L−1)
Ethyl lactate14 ± 110 ± 012 ± 25 ± 06 ± 06 ± 0 4 ± 0 a 5 ± 0 ab 6 ± 0 b ns***ns
Acetoine 8 ± 1 b 4 ± 1 a 3 ± 0 a 3 ± 03 ± 03 ± 04 ± 23 ± 2<LOD**ns
Acetol 92 ± 0 b 34 ± 6 a 30 ± 2 a 33 ± 335 ± 128 ± 110 ± 48 ± 011 ± 1********
2,3-butanediol levo729 ± 116776 ± 51618 ± 54712 ± 11691 ± 29692 ± 131565 ± 3671552 ± 422733 ± 110ns**ns
2,3-butanediol meso332 ± 17332 ± 15317 ± 13352 ± 3350 ± 6347 ± 8318 ± 81314 ± 90130 ± 28nsnsns
Methionol0.78 ± 0.050.82 ± 0.140.67 ± 0.100.41 ± 0.040.49 ± 0.060.44 ± 0.060.25 ± 0.050.22 ± 0.040.43 ± 0.07ns***ns
Acetates of higher alcohols (mg L−1)
Isoamyl acetate0.71 ± 0.030.67 ± 0.040.51 ± 0.101.25 ± 0.381.34 ± 0.350.76 ± 0.010.74 ± 0.090.83 ± 0.030.86 ± 0.06nsnsns
Hexyl acetate0.54 ± 0.210.36 ± 0.080.22 ± 0.160.26 ± 0.020.37 ± 0.040.42 ± 0.050.29 ± 0.040.22 ± 0.020.25 ± 0.01nsnsns
2-phenylethyl acetate 0.04 ± 0.00 ab 0.05 ± 0.00 b 0.03 ± 0.00 a 0.25 ± 0.010.22 ± 0.030.20 ± 0.010.08 ± 0.030.09 ± 0.010.06 ± 0.01nsnsns
∑ Acetates1.29 ± 0.201.08 ± 0.180.76 ± 0.141.76 ± 0.331.93 ± 0.351.38 ± 0.161.11 ± 0.201.14 ± 0.231.17 ± 0.24nsnsns
Esters (mg L−1)
Ethyl butyrate0.04 ± 0.000.04 ± 0.000.03 ± 0.02< LOD< LOD< LOD0.05 ± 0.000.05 ± 0.000.06 ± 0.00ns*ns
Ethyl hexanoate0.26 ± 0.010.29 ± 0.010.30 ± 0.090.44 ± 0.030.36 ± 0.010.39 ± 0.020.24 ± 0.020.25 ± 0.000.22 ± 0.00nsnsns
Ethyl octanoate0.34 ± 0.000.41 ± 0.010.38 ± 0.080.71 ± 0.070.54 ± 0.070.60 ± 0.110.29 ± 0.000.34 ± 0.010.32 ± 0.04nsnsns
Ethyl decanoate0.08 ± 0.020.09 ± 0.010.08 ± 0.010.30 ± 0.030.28 ± 0.050.29 ± 0.030.11 ± 0.010.12 ± 0.000.12 ± 0.02nsnsns
∑ Ethyl esters C6-C100.72 ± 0.030.83 ± 0.030.79 ± 0.211.45 ± 0.131.18 ± 0.131.29 ± 0.160.70 ± 0.020.76 ± 0.010.72 ± 0.06nsnsns
Volatile fatty acids (mg L−1)
Isobutyric acid3.40 ± 0.472.78 ± 0.172.76 ± 0.113.87 ± 0.304.15 ± 0.403.65 ± 0.15 1.76 ± 0.17 a 1.53 ± 0.10 a 2.48 ± 0.14 bns*ns
Butyric acid0.98 ± 0.110.92 ± 0.000.87 ± 0.081.48 ± 0.031.46 ± 0.031.37 ± 0.040.85 ± 0.080.86 ± 0.060.86 ± 0.00nsnsns
Isovaleric acid1.41 ± 0.011.35 ± 0.181.16 ± 0.152.59 ± 0.032.84 ± 0.202.81 ± 0.011.25 ± 0.101.29 ± 0.051.30 ± 0.03nsnsns
∑ Volatile fatty acids C4-C55.79 ± 0.595.05 ± 0.354.79 ± 0.197.94 ± 0.358.45 ± 0.647.84 ± 0.19 3.86 ± 0.01 a 3.68 ± 0.10 a 4.64 ± 0.11 bnsnsns
Hexanoic acid1.96 ± 0.272.27 ± 0.122.12 ± 0.362.94 ± 0.042.78 ± 0.152.75 ± 0.232.05 ± 0.222.22 ± 0.041.95 ± 0.09nsnsns
Octanoic acid1.84 ± 0.242.23 ± 0.221.88 ± 0.242.75 ± 0.312.36 ± 0.012.71 ± 0.321.97 ± 0.062.33 ± 0.062.05 ± 0.28nsnsns
Decanoic acid0.39 ± 0.150.43 ± 0.060.34 ± 0.661.22 ± 0.080.99 ± 0.031.12 ± 0.100.62 ± 0.060.69 ± 0.050.59 ± 0.02nsnsns
∑ Volatile fatty acids C6-C104.19 ± 0.664.93 ± 0.414.33 ± 0.656.91 ± 0.436.14 ± 0.196.58 ± 0.654.64 ± 0.225.24 ± 0.084.58 ± 0.39nsnsns
Free terpenes (µg L−1)
trans-linalool oxide (furan) †7.0 ± 0.85.4 ± 0.96.2 ± 0.15.2 ± 0.93.9 ± 0.34.1 ± 0.57.9 ± 0.27.8 ± 0.18.1 ± 0.7nsnsns
cis-linalool oxide (furan) †0.7 ± 0.00.7 ± 0.20.8 ± 0.10.9 ± 0.00.6 ± 0.10.7 ± 0.11.1 ± 0.00.7 ± 0.10.8 ± 0.2nsnsns
trans-linalool oxide (pyran) †4.5 ± 0.14.1 ± 0.14.0 ± 0.0 4.7 ± 0.2 b 4.2 ± 0.1 ab 3.4 ± 0.3 a 4.9 ± 2.54.9 ± 0.14.6 ± 0.7nsnsns
cis-linalool oxide (pyran) †1.2 ± 0.20.9 ± 0.00.9 ± 0.11.1 ± 0.21.1 ± 0.20.9 ± 0.21.1 ± 0.10.9 ± 0.11.1 ± 0.5nsnsns
Linalool (L)1.5 ± 0.01.4 ± 0.01.3 ± 0.23.6 ± 0.63.7 ± 0.73.4 ± 0.26.3 ± 0.26.2 ± 0.65.3 ± 0.5ns***ns
Hotrienol †0.5 ± 0.00.6 ± 0.00.6 ± 0.00.9 ± 0.21.1 ± 0.31.8 ± 0.61.3 ± 0.21.0 ± 0.11.3 ± 0.2ns*ns
A-terpineol (αT)2.2 ± 0.32.4 ± 0.12.3 ± 0.03.2 ± 0.02.8 ± 0.13.3 ± 0.34.5 ± 0.13.9 ± 0.73.1 ± 1.3ns**ns
Citronellol (C)<LOD<LOD<LOD 8.1 ± 1.4 b 5.2 ± 0.4 a 4.2 ± 0.6 a 7.5 ± 0.6 b 6.3 ± 0.1 ab 5.7 ± 0.3 a *nsns
Nerol (N)<LOD<LOD<LOD0.2 ± 0.20.2 ± 0.10.7 ± 0.32.2 ± 0.12.6 ± 0.42.7 ± 0.6ns***ns
Geraniol (G)<LOD<LOD<LOD5.6 ± 1.15.5 ± 0.26.1 ± 0.56.8 ± 0.47.3 ± 0.96.3 ± 0.1nsnsns
Σ Free terpenes (L + αT + C + N + G)3.6 ± 0.33.8 ± 0.43.6 ± 0.420.7 ± 1.717.2 ± 0.817.8 ± 1.127.3 ± 1.226.3 ± 1.123.1 ± 0.9ns***ns
Hodiol I †<LOD<LOD<LOD5.4 ± 0.55.5 ± 0.55.9 ± 0.6 5.7 ± 0.0 c 4.3 ± 0.0 b 3.7 ± 0.2 a ns**
2,7-dimethyloctane-4,5-diol †81.2 ± 6.4111.0 ± 17.891.5 ± 12.5160.9 ± 2.7169.3 ± 15.1152.7 ± 27.686.3 ± 17.376.0 ± 4.364.5 ± 5.3ns***ns

Different letters indicate significant differences between treatments at p < 0.05 for a given compound or sum of compounds on a given year; for a rapid identification, these values appear in bold. The significance of the year, treatment and their interaction is expressed as ns = not significant; * p < 0.05; ** p < 0.01; *** p < 0.001. LOD means limit of detection. † indicates that the compound is expressed in µg L−1 of internal standard.

The effect of the treatments differed among years. In 2012, higher concentrations of benzyl alcohol, acetoine and acetol appeared in wines from R, whereas 2-phenylethyl acetate was at lower concentrations in wines from SDI. In 2013, 1-propanol, 3-methyl-1-butanol and the sum of higher alcohols were greater in wines from DI, whereas the concentrations of trans-linalool-oxide (pyran) and citronellol were higher in wines from R. In 2014, wines from SDI had higher concentrations of benzyl alcohol, ethyl lactate, isobutyric acid and the sum of short-chained fatty acids, whereas wines from R had higher concentrations of citronellol and hodiol I (Table 6). Previous reports on the effect of irrigation on wine volatile composition offered contrasting results depending on grapevine variety, soil type, climate conditions and severity and time of occurrence of water stress [10,14,17]. For instance, Talaverano et al. [10] reported higher concentrations of alcohols, C6 compounds and volatile phenols in Tempranillo wines from rain-fed vines. Similarly, Vilanova et al. [14] observed a reduction in C6 compounds in Verdejo wines from irrigated vines. These findings contrasted with the results from the current study likely because of the different cultivar/rootstock combination and climate conditions. In fact, soil water availability was adequate for bud-break and the development of leaf surface in our study [7,19], so no differences in grapevine metabolism were expected. In this study, 19 volatile compounds appeared in wines at concentrations higher than their odor thresholds (Table 7), likely contributing to wine aroma [22,23]. Acetaldehyde, 1-propanol, isoamyl acetate, ethyl hexanoate, ethyl octanoate and isovaleric acid were present at concentrations more than 15 times higher than their respective odor thresholds, being consistent with previous studies [16,17]. Irrigation did not affect OAV of any compound, except on specific years (Table 7). For instance, OAV for 1-propanol was higher for the DI treatment in 2013, whereas that of benzyl alcohol was higher for SDI in 2014 (Table 7).
Table 7

Irrigation effects on the odor activity values (mean ± standard error) of Godello wines. Odor thresholds and descriptors for each compound are displayed. The significances of the year, treatment as well as their interaction are also shown. R = rain-fed, DI = drip irrigation, SDI = subsurface drip irrigation, T = Treatment, Y = Year.

Odor Threshold (µg L−1)Odor Descriptor201220132014TYT × Y
RDISDIRDISDIRDISDI
Ethyl acetate7500Pineapple8 ± 17 ± 18 ± 04 ± 03 ± 04 ± 05 ± 15 ± 05 ± 0ns*ns
Acetaldehyde500Fruity116 ± 2577 ± 256 ± 462 ± 671 ± 761 ± 462 ± 1255 ± 149 ± 1***ns
Higher alcohols
1-propanol750Alcohol16 ± 017 ± 116 ± 3 29 ± 1 ab 32 ± 0 b 26 ± 0 a 31 ± 632 ± 522 ± 3ns**ns
2-methyl-1-propanol40000Alcohol2 ± 01 ± 02 ± 01 ± 01 ± 01 ± 01 ± 01 ± 01 ± 0ns***ns
3-methyl-1-butanol30000Alcohol8 ± 19 ± 19 ± 010 ± 010 ± 010 ± 07 ± 07 ± 08 ± 0nsnsns
Other alcohols
cis-3-hexen-1-ol400Grass1 ± 01 ± 01 ± 01 ± 01 ± 01 ± 01 ± 01 ± 0<1ns**ns
Benzyl alcohol620Blackberry 6 ± 0 b 4 ± 0 a 4 ± 0 a 4 ± 04 ± 03 ± 0 3 ± 0 a 4 ± 0 ab 5 ± 0 b ns****
2-phenylethanol14000Rose4 ± 04 ± 03 ± 03 ± 04 ± 04 ± 0 2 ± 0 a 3 ± 0 b 2 ± 0 a ns**ns
Acetates of higher alcohols
Isoamyl acetate30Banana24 ± 122 ± 117 ± 342 ± 1345 ± 1226 ± 125 ± 328 ± 129 ± 2nsnsns
Esters
Ethyl butyrate20Fruity2 ± 02 ± 02 ± 1<1<1<13 ± 03 ± 03 ± 0nsnsns
Ethyl hexanoate14Fruity19 ± 121 ± 121 ± 631 ± 226 ± 128 ± 218 ± 215 ± 015 ± 0nsnsns
Ethyl octanoate5Fruity67 ± 182 ± 276 ± 16141 ± 13108 ± 14120 ± 2259 ± 167 ± 165 ± 9nsnsns
Ethyl decanoate200Grape<1<1<12 ± 01 ± 01 ± 01 ± 01 ± 01 ± 0nsnsns
Volatile fatty acids
Isobutyric acid2300Cheese1 ± 01 ± 01 ± 02 ± 02 ± 02 ± 0 <1 a <1 a 1 ± 0 b nsnsns
Butyric acid173Cheese6 ± 15 ± 05 ± 09 ± 08 ± 08 ± 05 ± 05 ± 05 ± 0nsnsns
Isovaleric acid33Cheese43 ± 041 ± 535 ± 578 ± 186 ± 685 ± 038 ± 339 ± 239 ± 1nsnsns
Hexanoic acid3000Cheese1 ± 01 ± 01 ± 01 ± 01 ± 01 ± 01 ± 01 ± 01 ± 0nsnsns
Octanoic acid500Rancid4 ± 04 ± 04 ± 06 ± 15 ± 05 ± 14 ± 05 ± 04 ± 1nsnsns
Decanoic acid1000Rancid<1<1<11 ± 01 ± 01 ± 01 ± 01 ± 01 ± 0nsnsns

Different letters indicate significant differences between treatments at p < 0.05 for a given compound on a given year; for a rapid identification, these values appear in bold. The significance of the year, treatment and their interaction is expressed as ns = not significant; * p < 0.05; ** p < 0.01; *** p < 0.001. Odor thresholds have been taken from Guth [22] and Ferreira et al. [23].

3.5. Sensory Profiles of Godello Wines

No clear differences among treatments were detected for color (Supplementary Figure S1). Aroma descriptors (Figure 2) were similar for the three treatments in 2012 and 2013; however, wines from R received greater marks to fresh, dry and tropical fruit descriptors than wines from the irrigated treatments, especially those from SDI (Figure 2). Palate descriptors did not reveal a clear preference of the tasters for wines from a given treatment (Supplementary Figure S1). As reported for other white varieties [8], Godello wines from irrigated vines were more acidic but showed an aroma complexity similar to that of wines from the rain-fed treatment. Overall, the slight differences observed in the sensory characteristics were likely caused by the low level of water stress experienced by grapevines [7], which did not affect the Godello vine performance. When applying a PCA to the aroma descriptors and the families of volatile compounds quantified in Godello wines (Figure 3), the first two PC explained 66% of the variability in the samples. PC1 explained 41% of the variance and was positively related to the citric and herbaceous descriptors, and negatively associated to aroma intensity and the concentration of volatile fatty acids, higher alcohols and esters. The PC2 explained 25% of the variance and was positively related to tropical fruit and floral descriptors, and negatively related to the concentration of terpenes. The projection of the samples on this plane allowed for a clear differentiation of the studied years, but not by treatments. Wines from 2012 were on the positive side of both PC due to their marks for tropical fruit and floral descriptors. In contrast, wines from 2013 had greater concentrations of higher alcohols and volatile fatty acids, appearing on the negative side of both PC. Finally, wines from 2014 were on the positive side of PC1 due to their greater concentrations of free terpenes and higher marks for citric and herbaceous descriptors. Aggrupation of wines from the same year is consistent with other reports [35], proving that in-season effects caused by irrigation were less important than the influence of weather conditions.
Figure 3

Principal component analysis (PCA) of Godello wines: Biplot for the first two principal components (PC) for aroma descriptors (in upper-case letters) and the different families of volatile compounds quantified (in lower-case letters). R = Rain-fed, DI = Drip irrigation, SDI = Subsurface drip irrigation.

4. Conclusions

Under the climate conditions of Northwest Iberian Peninsula, irrigation had a low incidence on amino acid and volatile composition of Godello musts and wines, because of the weak level of water stress experienced by vines from veraison to harvest. Despite this, rain-fed vines produced wines with lower acidity and higher pH than those from the treatment under subsurface drip irrigation. Moreover, irrigation decreased the concentration of volatile compounds providing fruity aromas, such as acetaldehyde and isoamyl acetate, leading to slight differences in wine perception by tasters.
  2 in total

1.  Revealing the Usefulness of Aroma Networks to Explain Wine Aroma Properties: A Case Study of Portuguese Wines.

Authors:  Sílvia Petronilho; Ricardo Lopez; Vicente Ferreira; Manuel A Coimbra; Sílvia M Rocha
Journal:  Molecules       Date:  2020-01-09       Impact factor: 4.411

2.  Amino Acid Profiles to Differentiate White Wines from Three Autochtonous Galician Varieties.

Authors:  José Manuel Mirás-Avalos; Yolanda Bouzas-Cid; Emiliano Trigo-Córdoba; Ignacio Orriols; Elena Falqué
Journal:  Foods       Date:  2020-01-21
  2 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.