Next Article in Journal
Variation Characteristics of Glucosinolate Contents in Leaf Mustard (Brassica juncea)
Previous Article in Journal
UAV Multispectral Data: A Reliable Approach for Managing Phosphate-Solubilizing Bacteria in Common Bean
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Carbon-Isotope Ratio (δ13C) and Phenolic-Compounds Analysis in Authenticity Studies of Wines from Dealu Mare and Cotnari Regions (Romania)

by
Andreea Popîrdă
1,
Camelia Elena Luchian
1,*,
Lucia Cintia Colibaba
1,
Elena Cornelia Focea
1,
Sebastien Nicolas
2,
Laurence Noret
2,
Ionel Bogdan Cioroiu
3,
Régis Gougeon
2 and
Valeriu V. Cotea
1
1
Department of Horticultural Technologies, Faculty of Horticulture, Iasi University of Life Sciences, 3rd M. Sadoveanu Alley, 700490 Iasi, Romania
2
PAM UMR A 02.102, Institut Universitaire de la Vigne et du Vin–Jules Guyot, AgroSup Dijon/Université Bourgogne Franche-Comté, F-21000 Dijon, France
3
Research Centre for Oenology, Romanian Academy-Iasi Branch, 8 Carol I, 700505 Iasi, Romania
*
Author to whom correspondence should be addressed.
Agronomy 2022, 12(10), 2286; https://doi.org/10.3390/agronomy12102286
Submission received: 31 August 2022 / Revised: 17 September 2022 / Accepted: 19 September 2022 / Published: 23 September 2022

Abstract

:
In food quality, authenticity is one of the most important issues in the context of ensuring the safety and security of consumers, but it becomes even more important when wine is involved as this produce has become one of the most counterfeited foods in the world. A batch of 28 wines from Romanian grape varieties obtained in Dealu Mare and Cotnari regions was analysed from a physical–chemical point of view in order to assess the possibility of discriminating them according to geographical origin and variety. The samples were subjected to isotopic analysis, thus determining δ13C isotopic ratio using an elemental analyser EA, VarioMicroCube; while the targeted analysis of phenolic compounds was performed by UHPLC technique, using an Acquity UPLC H-Class. The basic physico-chemical analysis of the samples was carried out using FTIR spectroscopy (OenoFoss). Statistical analyses were performed using the TIBCO STATISTICAL SOFTWARE and the statistical test applied was the Tukey HSD test. Therefore, identified phenolic compounds such as hydroxytyrosol and coumaric acid are considered varietal markers. Tyrosol, dimers B1 and B2 and also catechin and epicatechin are indicators of geographical origin.

1. Introduction

In order to meet the needs of producers, consumers and authorities, and also because of the important value of the wine market, in recent years there has been a global interest in developing analytical methods that provide reliable information on issues related to the authenticity of wines. These methods are meant to diminish and prevent the falsification of wines and to guarantee their quality in the import–export market [1].
Factors such as grape variety, production area, year of harvest, type of soil, technology applied, and climatic conditions are all important in determining the chemical composition of wines [2].
Determining the isotopic ratios of elements such as H, C, N, O has proven, over time, to be a useful tool in determining the geographical origin of food [3,4], a tool that helps to confirm or to deny the information on the label or to verify the conformity of the products with legislation. Regarding wine, the determination of the ratios of stable isotopes (δ2H, δ18O or δ13C, constituents of water or ethanol extracted from the analysed wine) can provide information on the origin of raw materials [5]. The carbon isotope ratio in ethanol after a previous extraction from wine provides relevant information to validate the geographical origin of wines. Therefore, a correlation between the carbon isotope ratio and the phenolic compounds profile can provide an overview of wines of a certain variety or region. In addition, in recent years, untargeted metabolomic analysis (mass spectrometry (MS) hyphenated with separation techniques such as liquid chromatography (LC)) coupled with statistical analysis have allowed good discrimination between wines according to the area of origin, variety, harvest year, etc. [6,7].
The type and concentration of the phenolic compounds are very different depending on grape variety, climatic conditions in each area and the applied wine-making technology [8,9]. Generally, probably the major determinant factor for the variation in the polyphenolic content of different wines throughout the world is the amount of sunlight to which the grapes are exposed during cultivation, but also the average temperature and the water status [10,11].
In Romania, viticulture has an old tradition, over 2000 years old [12], even if at present many of the old, local grape varieties are no longer used in production or they occupy very small areas. Currently, Romania ranks 6th in Europe and 12th in the world in terms of wine production, with more than 180 thousand hectares of vines. In 2021, there was an increase in wine production by 37% compared to the previous year; the Romanian wine market being one of the most dynamic markets in Europe [13]. Two of the local grape varieties that occupy large surfaces and from which wines are highly appreciated by consumers are Fetească neagră and Fetească albă. Busuioacă de Bohotin has recently become a widely cultivated grape variety in Romania and appreciated for its aromatic rosé wines. It is also considered a Romanian variety due to the long time (centuries) it has been cultivated in Romania, although its origins are Greek [14].
The study’s main aim was to correlate, by applying statistical analysis, the chromatographic analysis of phenolic compounds with δ13C isotopic ratio in order to identify markers for the viticultural areas and functions towards discriminating between varietals.

2. Materials and Methods

The studied wines are in fact a set of 28 white, rosé and red commercial wines made in 2019, 2018 and 2017 from the 3 Romanian varieties mentioned above: Fetească albă, Busuioacă de Bohotin and Fetească neagră. These samples have their origins in two of the most important vineyards in Romania, namely: Cotnari (located in the northeast of the country, with 2037 ha of vines) and Dealu Mare (“homeland of red wine”, with 12,327 ha). Cotnari area is characterized by a hilly landscape with slopes of between 3 and 35%, while altitudes vary from 105 to 315 m above sea level. The climate is temperate, the summers are characterised by droughts while winters are cold due to a specific icy wind ”crivăț”, the absolute minimal temperatures can reach −29 °C [15]. In Dealu Mare, conditions are a little bit different. Slopes under vines vary between 8% and 30%. The average annual temperature varies between 10.8 °C and 11.2 °C. Winters are relatively short with cold bursts in January and the first half of February. The average temperature of the warmest month (July) is + 22.4 °C, with variations between + 20.7 °C and + 25.6 °C [16].
The grape varieties, the regions, the wine colours, the years of vintage and the physico-chemical characteristics are listed in Table 1.
The physico-chemical analyses were determined using FTIR spectroscopy (OenoFossTM type 4101, FOSS, Hilleroed, Denmark).
Targeted analysis of phenolic compounds was performed using UHPLC technique, using an Acquity UPLC H-Class (Waters, Milford, MA, USA) equipped with a quaternary pump and an autosampler that was coupled to a fluorimetric detector (FLD) and one with diode (DAD).
The column (Restek brand, model Report ARC-18, particle size 1.8 µm, length 150 mm and internal diameter of 2.1 mm) was thermostated at 35 °C, and the samples were kept at 12 °C. The wine was filtered before injection with PTFE filters with a porosity of 0.45 µm (Restek, Lisses, France). Composition of the mobile phase: an aqueous eluent (A) was prepared, containing 0.28% trifluoroacetic acid and 5% methanol (pH = 1.6). Solvent B was methanol 100%. The flow rate was constant at 0.30 mL min−1. The optimized elution system consisted of a gradient obtained as presented in Table 2. The diode array detector was set at 320 nm for caftaric acid, gentisic acid, and caffeic acid; at 305 nm for coumaric acid; at 280 nm for gallic acid; and at 260 nm for protocatechuic acid. The fluorescence detector was set at λex = 270 nm and λem = 322 nm for hydroxytyrosol, tyrosol, catechin, epicatechin, proanthocyanidin B1, and proanthocyanidin B2, while for trans-resveratrol it was set at λex = 330 nm and λem = 374 nm. Chemically pure standards of phenolic acids (gallic, protocatechuic, caftaric, caffeic, coumaric), stilbenes (trans-resveratrol), flavonols (catechin, epicatechin), hydroxytyrosol and tyrosol, procyanidin dimer B1 and procyanidin dimer B2 were purchased from Sigma-Aldrich (St. Louis, MO, USA). The purity of all phenolic standards was greater than 95%. The individual stock solutions (1000 ppm) were prepared in pure methanol and kept at −20 °C in the dark. The calibration standards were freshly prepared on the day of analysis by diluting the appropriate working solution with the initial medium phase solution. The concentration range was selected based on the sensitivity of UHPLC-PDA for each polyphenol [17].
δ13C measurements were carried out using an elemental analyser (EA, VarioMicroCube, Elementar, F-69623 Villeurbanne, France) which was coupled to isotopic ratio monitoring by mass spectrometry (IRMS, Isoprime/Elementar, F-69623 Villeurbanne, France). Internationally certified reference materials USGS40 (with δ13C org = −26.39‰) and caffeine IAEA-600 (δ13C org = −27.77‰) were used for calibration. As carrier gas, Helium (Air Liquide, 5.6) was used, while carbon dioxide (Air Liquide, 4.8) was used as the reference gas. A total of 10 μL of wine was added to tin cups, and these were then introduced into the oxidation tube (950 °C) under helium flow (200 mL min−1) and oxygen flow (30 mL min−1), the temperature of the reduction furnace being fixed at 550 °C. Tin cups (5 × 9 mm) and copper wire were supplied by Elementar (France). The gases obtained by combustion were dried and eluted in a specific column that physically retains CO2 (60 °C) and then releases it with a higher temperature (110 °C). An open split system allowed adjustment of gas withdrawal to the IRMS; current trap was fixed at 200 μA. The total duration of one measurement cycle was 600 s. The masses measured by IRMS were m/z 44 and 45 corresponding to CO2 without and with 13C, respectively. The isotope ratio is expressed as a relative deviation, δ13C in per mil (‰), from the international standard, Vienna-Pee Dee Belemnite (VPDB), according to δ13C (‰) = 1000 × [(Rs/Rst) − 1], where R corresponds to the carbon 13 isotope ratio of the sample(s) and the standard (st) [18].
Statistical analyses were performed using the TIBCO STATISTICAL SOFTWARE and the statistical test applied was the Tukey HSD test. Tukey’s honest significance test is a single-step multiple comparison procedure for evaluation of statistical significance between multiple parameters of the evaluation of global variance for the same characterization (region or grape variety). In addition, for the discrimination by the variety according to the phenolic content and δ13C ratio, the discriminant function analysis test was applied using the same software [19].

3. Results and Discussions

The study carried out the analyses of phenolic compounds and stable isotopes for the 28 wines from three consecutive years of harvest, from three local grape varieties (Fetească neagra, Busuioacă de Bohotin and Fetească albă) and with two different vineyards: Cotnari and Dealu mare. The results of δ13C analyses are summarized in Table 3, while the phenolic compounds concentrations are presented in Table 4.
It can be observed from Table 4 that the values of the parameter δ13C vary between 28.18 VPDB (‰)—the minimum value reported in the case of a wine from Dealu Mare obtained from the Fetească neagră variety—and −25.66 VPDB (‰)—the maximum value reported in the case of a wine from Busuioacă de Bohotin obtained in the Cotnari vineyard. All other values of this parameter fall within the reported range. A first observation is that the values in this range are generally specific to wines produced in Romania; there are data in the literature to confirm this observation [1,2,20].
Twelve phenolic compounds were targeted and determined for this study. The analyses were performed in duplicate. Table 4 summarizes the results obtained from the analysis of the 12 phenolic compounds grouped as follows: phenolic acids (gallic, protocatechuic, caftaric, caffeic, coumaric), stilbenes (trans-resveratrol), flavonols (catechin, epicatechin), hydroxytyrosol and tyrosol, procyanidin dimer B1 and procyanidin dimer B2. The values thus obtained and expressed in mg/L are comparable to those reported in the literature for Romanian wines [21,22].
From a simple interpretation of the averages in Table 4, variations in the content of phenolic compounds between the tested wines can be observed, depending on: geographical origin, grape variety used or type of wine. Thus, the levels of gallic acid accumulated by the wines obtained from the Fetească neagră variety are considerably higher in the case of the Cotnari vineyard, while the wines obtained from the same variety in Dealu Mare have lower concentrations of the same compound.

3.1. Statistical Tests

3.1.1. ANOVA Evaluation

Considering every categorization group, a normality projection of data was performed in order to verify the distribution of variable results (Figure 1).
Normality of distribution permitted evaluation after one-way ANOVA of the initial data for the variables included in the study considering the grouping of variables region and variety. This revealed that every parameter of characterization was significant at p < 0.05 in the wines that were proposed for analysis and characterization. Evaluation was performed taking Wilks Lamda, which indicated percent variance in the dependent variables not explained by the differences in the levels of categorization. The values of significance for the variables were of 0.00054 for Wilks Lambda (F = 38.75) for the variety as first grouping variable and 0.1072 (F = 7.69) for region as second grouping variable. Since Wilks values are close to zero and F values are positive, the variables contribute to the evaluation model.

3.1.2. Tukey HSD Test

From the statistical analysis performed (Tukey HSD test) whose results are presented in Table 5, it is confirmed that gallic acid can be considered, in the case of the analysed wines, as discriminating by area, p factor being p < 0.05.
Regarding protocatechuic acid, caftaric acid, caffeic acid, coumaric acid, trans-resveratrol and hydroxytyrosol, the differences between the averages of these compounds reported in Table 4 (for Cotnari and Dealu Mare wines) are smaller. Thus, it is difficult to suggest that these compounds may constitute area markers. This hypothesis is also confirmed by the statistical test applied; in their case the values of the p coefficient always exceeded the value of 0.05 (as can be seen in Table 5).
In addition, according to the results of the statistical analyses highlighted in Table 5, the variables tyrosol, dimer B1, dimer B2, catechins and epicatechins have the p factor values p < 0.005, which again constitute statistically significant differences, translating into the possibility of discriminating the two areas analysed using these variables.
The same statistical test was applied in order to discriminate the analysed varieties according to the same variables: the concentrations of phenolic compounds (Table 6). Thus, it was observed that the recorded values of all the analysed phenolic compounds (except tyrosol) showed statistically significant differences between the varieties Fetească neagră and Fetească albă, while only some of them showed statistically significant differences between the other two combined pairs (Fetească albă/Busuioacă de Bohotin and Fetească neagră/Busuioacă de Bohotin). The values of protocatechuic acid also discriminated between the varieties Fetească neagră and Busuioacă de Bohotin, as well as those registered for caffeic acid. The concentrations of coumaric acid as well as hydroxytyrosol can act as discriminants for all three grape varieties. The variable contents of tyrozol recorded statistically significant differences between the varieties Fetească albă and Busuioacă de Bohotin, while the contents of dimers B1 and B2, as well as those of catechin and epicatechin, discriminated between the varieties Fetească neagră and Busuioacă de Bohotin.
A conclusion can be drawn from the obtained results: the isotopic analysis (δ13C) coupled with the statistical analysis makes possible to discriminate between grape varieties and areas. Thus, a clear differentiation could be observed between the wines obtained from the Busuioacă de Bohotin and Fetească neagră varieties and produced in the Cotnari vineyard, as shown in Table 7.
The same observation can be made in the case of wines obtained from the same varieties, but this time in the Dealu Mare vineyard. Within the Dealu Mare vineyard, the wines from the Fetească albă and Fetească neagră varieties could also be discriminated between with the help of this variable (δ13C). The values 0.007517 and 0.001653 of the p parameter in Table 7 show a clear discrimination between the two areas achieved using the δ13C isotopic ratio.
Table 4. Phenolic compound concentrations.
Table 4. Phenolic compound concentrations.
SampleGallic AcidProtocatechuic AcidCaftaric AcidCaffeic AcidCoumaric AcidTrans-ResveratrolHydroxy-TyrosolTyrosolProcyanidin Dimer B1CatechinProcyanidin Dimer B2Epicatechin
mg/Lmg/Lmg/Lmg/Lmg/Lmg/Lmg/Lmg/Lmg/Lmg/Lmg/Lmg/L
Fetească neagră Cotnari (5)Min43.592.2734.5519.026.643.122.5613.4817.7119.709.8112.76
Max98.274.8262.7924.5411.7212.613.2818.4228.2028.6713.0821.09
Average71.383.5046.5222.559.017.862.9115.4221.4923.5111.1716.55
SD21.190.9010.842.191.884.010.302.074.123.561.703.55
Busuioacă de Bohotin Cotnari (7)Min30.441.0813.923.202.650.431.478.924.366.740.362.68
Max89.311.7358.706.825.690961.9311.056.509.650.975.96
Average50.401.3742.454.774.330.751.6810.115.438.120.714.05
SD18.480.2614.801.220.910.190.180.810.671.080.191.21
Fetească albă Cotnari (4)Min1.790.3927.561.250.41Nd0.0014.122.064.320.331.36
Max59.180.7734.122.601.65Nd0.3321.803.326.830.632.77
Average25.190.5429.871.700.82Nd0.2017.122.605.470.442.05
SD24.350.172.910.610.58Nd0.163.660.5711.70.140.62
Fetească neagră Dealu Mare (4)Min33.682.4443.2612.155.513.772.2924.6336.7733.7520.2627.80
Max70.993.6688.6217.308.8512.073.3132.1540.6940.8124.4637.32
Average49.563.1663.9914.447.366.782.7228.4338.9237.0622.7132.39
SD17.270.5719632.201.433.630.463.531.763.291.783.90
Busuioacă de Bohotin Dealu Mare (4)Min0.500.7921.050.711.480.290.797.154.477.660.255.19
Max33.261.9736.196.104.091.793.6910.118.7213.001.528.32
Average10.061.2227.783.963.020.972.338.755.969.660.976.95
SD15.540.546.342.331.150.641.431.401.982.320.541.36
Fetească albă Dealu Mare (4)Min0.790.7419.573.322.11Nd0.2914.651.404.030.192.00
Max12.571.6827.194.813.78Nd1.1622.072.376.830.852.42
Average5.921.2922.553.712.65Nd0.6017.661.825.260.412.26
SD4.950.423.270.730.76Nd0.413.400.471.330.300.19
Table 5. Tukey HSD test results.
Table 5. Tukey HSD test results.
Tukey HSD Test
No.RegionVariable: Gallic acid, MS = 405.39, df = 24.000Variable: Protocatechuic acid, MS = 0.51870, df = 24.000Variable: Caftaric acid, MS = 222.41, df = 24.000Variable: Caffeic acid, MS = 11.421, df = 24.000Variable: coumaric acid, MS = 2.3497, df = 24.000Variable: trans-resveratrol, MS = 718.87, df = 24.000
{1}
50.652
{2}
21.849
{1}
1.8283
{2}
1.8919
{1}
40.575
{2}
38.104
{1}
9.5612
{2}
7.3711
{1}
4.9151
{2}
4.3453
{1}
42.503
{2}
22.892
1Cotnari 0.001133 0.819332 0.668389 0.102783 0.340172 0.067559
2Dealu Mare0.001133 0.819332 0.668389 0.102783 0.340172 0.067559
Tukey HSD Test
No.RegionVariable: Hydroxytyrosol, MS = 0.38400, df = 24.000Variable: Tyrosol, MS = 25.605, df = 24.000Variable: Procyanidin dimer B1, MS = 51.513, df = 24.000Variable: Catechin, MS = 35.256, df = 24.000Variable: Procyanidin dimer B2, MS = 20.830, df = 24.000Variable: Epicatechin, MS = 33.531, df = 24.000
{1}
1.6937
{2}
1.8820
{1}
13.520
{2}
18.277
{1}
9.7428
{2}
15.566
{1}
12.268
{2}
17.327
{1}
3.9133
{2}
8.0264
{1}
7.4569
{2}
13.864
1Cotnari 0.433942 0.021524 0.044233 0.035432 0.026882 0.008065
2Dealu Mare0.433942 0.021524 0.044233 0.035432 0.026882 0.008065
MS—mean sum of squares, df—degrees of freedom, red values—p-value is statistically representative for discrimination; the numbers in parentheses correspond to the regions (in the matrix).
Table 6. Tukey HSD test results.
Table 6. Tukey HSD test results.
Tukey HSD Test
No.VarietyVariable: Gallic acid, MS = 405.39, df = 24.000Variable: Protocatechuic acid, MS = 0.51870, df = 24.000Variable: Caftaric acid, MS = 222.41, df = 24.000Variable: Caffeic acid, MS = 11.421, df = 24.000Variable: coumaric acid, MS = 2.3497, df = 24.000Variable: trans-resveratrol, MS = 718.87, df = 24.000
{1}
55.562
{2}
35.758
{3}
17.666
{1}
3.1021
{2}
1.3759
{3}
0.8286
{1}
50.957
{2}
36.978
{3}
27.159
{1}
17.660
{2}
4.2259
{3}
2.6210
{1}
7.8205
{2}
3.8622
{3}
1.4424
{1}
86.976
{2}
7.7264
{3}
0.0000
1Fetească neagră 0.0828620.002408 0.0001620.0001310.1019110.1019110.009557 0.0001290.000129 0.0001370.000129 0.0001300.000130
2Busuioacă de Bohotin0.082862 0.1726950.000162 0.277098 0.3763520.000129 0.000137 0.0089730.000130 0.823642
3Fetească albă0.0024080.172695 0.0001310.277098 0.0095570.376352 0.0001290.592940 0.0001290.008973 0.0001300.823642
Tukey HSD Test
No.VarietyVariable: Hydroxy-tyrosol, MS = 0.38499, df = 24.000Variable: Tyrosol, MS = 25.605, df = 24.000Variable: Procyanidin dimer B1, MS = 51.513, df = 24.000Variable: Catechin, MS = 35.256, df = 24.000Variable: Procyanidin dimer B2, MS = 20.830, df = 24.000Variable: Epicatechin, MS = 33.531, df = 24.000
{1}
2.8841
{2}
1.6321
{3}
0.4128
{1}
19.797
{2}
10.356
{3}
17.680
{1}
26.774
{2}
5.4234
{3}
2.1837
{1}
27.471
{2}
8.4902
{3}
5.1588
{1}
14.695
{2}
0.8605
{3}
0.3586
{1}
21.897
{2}
4.7014
{3}
2.1414
1Fetească neagră 0.0004190.000129 0.0008560.677037 0.0001300.000129 0.0001290.000129 0.0001300.000131 0.0001300.000130
2Busuioacă de Bohotin0.000419 0.0013490.000856 0.0167870.000130 0.6248980.000129 0.4876310.000130 0.9719990.000130 0.637743
3Fetească albă0.0001290.001349 0.6770370.016787 0.0001290.624898 0.0001290.487631 0.0001310.971999 0.0001300.637743
MS—mean sum of squares, df—degrees of freedom, red values—p-value is statistically representative for discrimination; the numbers in parentheses correspond to the varieties (in the matrix).
Table 7. Tukey HSD test results.
Table 7. Tukey HSD test results.
No.RegionVarietyTukey HSD Test; Variable δ13C; Approximate Probabilities for Post Hoc Tests; Error: Between MS = 0.9002, df = 22.00
{1}
−27.21
{2}
−26.23
{3}
−26.6
{4}
−28.01
{5}
−27.10
{6}
−26.77
1CotnariFetească neagră 0.0002920.0572730.0075170.9935740.285852
2CotnariBusuioacă de Bohotin0.000292 0.3949610.0001440.0016530.079276
3CotnariFetească alba0.0572730.394961 0.0001530.2071910.959551
4Dealu MareFetească neagră0.0075170.0001440.000153 0.0036710.000220
5Dealu MareBusuioacă de Bohotin0.9935740.0016530.2071910.003671 0.636627
6Dealu MareFetească albă0.2858520.0792760.9595510.0002200.636627
MS—mean sum of squares, df—degrees of freedom, red values—p-value is statistically representative for discrimination; the numbers in parentheses correspond to the mic between region and variety (in the matrix).

3.1.3. Statistical Tests—Discriminant Function Analysis (DFA)

The discriminant analysis made it possible to discriminate between varieties, the variables in this case being both the concentrations of the identified phenolic compounds and the results of the isotopic analysis, but not taking into account the origin of the wines because factors are considered as the average of the determined values in each region (Cotnari and Dealul Mare).
The independent contributions of each variable to the variety discrimination model were verified. It was found that all variables contribute differently to the discrimination between grape varieties. The hole model, based on the variables included, gave Wilks Lambda (WL) of 0.00054 for p < 0.05; thus, it was established that there are sufficient variables that influence the model and the discrimination between the functions can be performed. In this situationm only some variables have a greater influence for the determination of a group, namely, caffeic acid (WL-0.0016, p = 0.00052), coumaric acid (WL-0.000934, p = 0.027), trans-resveratrol (WL-0.00096, p = 0.022), tyrosol (WL-0.0014, p = 0.0013) and dimer B1 (WL-0.00097, p = 0.019).
To understand how these variables contribute to the identification of each group (respectively, the grape variety), in this sense, it was necessary to differentiate the functions on the basis of which the groups are individualized, with each function contributing to a greater or lesser extent to the discriminant analysis (Table 8).
Thus, two functions of discrimination were identified, called roots and which are statistically relevant because p < 0.05; therefore, there are two ways of interpreting how the variables (concentrations of phenolic compounds and δ13C ratios) contribute to the separation of groups (grape varieties).
As can be seen in Table 9, the first canonical differentiation function (root 1) is more strongly influenced by variables such as caftaric acid, caffeic acid, and epicatechin; while root 2 is more strongly influenced by coumaric acid, trans-resveratrol, and procyanidin dimer B1, procyanidin dimer B2, hydroxytyrosol, tyrosol, catechin and epicatechin. A study of Montenegrin wines also allowed their discrimination with the help of the profile of polyphenolic compounds [23].
Based on the values of the root functions (as shown in Table 10) it could be shown that the first discrimination function (root 1—which is strongly influenced by the variables caffeic acid, caftaric acid and epicatechin) is the one that discriminates the Fetească neagră variety (and, implicitly, red wines) from the varieties Busuioacă de Bohotin and Fetească albă (meaning rosé and white wines). As can be seen in Figure 2, root 1 cannot make a strong discrimination between the varieties Busuioacă de Bohotin and Fetească albă.
On the other hand, the values of root 2 highlight a good differentiation between all three varieties. Thus, it can discriminate both red from white and rosé wines (respectively, the Fetească neagră variety of Busuioacă de Bohotin and Fetească albă varieties) but also between wines obtained from Busuioacă de Bohotin and Fetească albă varieties.
Table 10. Means of canonical variables.
Table 10. Means of canonical variables.
GroupRoot 1Root 2
Fetească neagră−16.87360.23597
Busuioacă de Bohotin7.3757−3.38780
Fetească albă8.84114.39277
Figure 2. DFA plot of discriminating wine samples by grape variety according to the composition of the phenolic profile.
Figure 2. DFA plot of discriminating wine samples by grape variety according to the composition of the phenolic profile.
Agronomy 12 02286 g002
All the data used for the statistical analyses can be found in the Supplementary Materials.
In the literature were also reported situations when phenolic compounds were markers of the discrimination of varieties. An example is the study conducted on 52 wines of different origins (Romania, Bulgaria and Moldova) which demonstrated the possibility of using phenols, determined by liquid chromatography, to distinguish commercial table wines by variety [24].
The phenolic composition of wines is significantly influenced by genetic and environmental factors, but also an important contribution comes from technology. For this reason, the authentication of wines using phenolic imprinting can sometimes be quite a challenging task [25]. This is the reason why the identification and quantification of phenolic compounds cannot, by itself, answering authenticity issues. It can only be used as an additional verification method.

4. Conclusions

The Tukey Test allowed the discrimination of Romanian wines obtained from three different grape varieties in two separate regions. The discrimination by area and by grape variety was performed according to phenolic composition and δ13C.
The phenolic compounds that could be considered grape variety markers were: hydroxytyrosol and coumaric acid.
The phenolic compounds that could be considered area markers were: gallic acid, tyrosol, B1 and B2 dimers, catechin and epicatechin.
The results of δ13C coupled with statistical analysis made it possible to discriminate between different grape varieties and areas.
The characteristics of these wines were determined by their raw matter, geographical origin as well as technological processes applied specific to each region.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/agronomy12102286/s1, Table S1: Quantitative analysis of phenolic compounds and δ13C values in analysed samples.

Author Contributions

Conceptualization, A.P., R.G. and C.E.L.; methodology, R.G. and C.E.L.; software, A.P. and I.B.C.; validation, L.C.C., V.V.C., R.G. and C.E.L.; formal analysis, A.P., E.C.F. and S.N.; investigation, A.P. and S.N.; resources, V.V.C., L.N. and R.G.; data curation, A.P. and I.B.C.; writing—original draft preparation, A.P. and L.C.C.; writing—review and editing, A.P. and C.E.L.; visualization, L.C.C.; supervision, V.V.C., R.G. and C.E.L. These authors have contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

Funding from the Unesco Chair “Culture et Traditions du Vin” from University of Burgundy, research convention no. 2021PRE00461, received on 21 July 2021, signed on 2 September 2021.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This research was supported by the UNESCO Chair “Culture et Traditions du Vin” of the University of Burgundy, and the IUVV Dijon–PAM Laboratory at University of Burgundy. The samples were acquired by the Iași University of Life Sciences, Romania. We thank Anne-Lise Santoni, technical manager of the GISMO platform from University of Burgundy, for the δ13C measurements.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Geana, E.I.; Popescu, R.; Costinel, D.; Dinca, O.R.; Ionete, R.E.; Stefanescu, I.; Artem, V.; Bala, C. Classification of red wines using suitable markers coupled with multivariatestatistic analysis. Food Chem. 2016, 192, 1015–1024. [Google Scholar] [CrossRef] [PubMed]
  2. Dinca, O.-R.; Ursu, S.O.; Costinel, D.; Miricioiu, M.G.; Radu, G.-L.; Popa, V.; Campeanu, C.B.; Elena, R. Samburesti wines characterization in terms of their stable isotope content. U.P.B. Sci. Bull. Ser. B 2015, 77, 177–188. [Google Scholar]
  3. Versini, G.; Monetti, A.; Reniero, F. Monitoring authenticity and regional origin of wines by natural stable isotope ratios analysis. ACS Symp. Ser. 1991, 661, 113–130. [Google Scholar]
  4. Dinca, O.R.; Ionete, R.E.; Popescu, R.; Costinel, D.; Radu, G.-L. Geographical and Botanical Origin Discrimination of Romanian Honey Using Complex Stable Isotope Data and Chemometrics. Food Anal. Methods 2015, 8, 401–412. [Google Scholar] [CrossRef]
  5. Costinel, D.; Tudorache, A.; Ionete, R.E.; Vremera, R. The Impact of Grape Varieties to Wine Isotopic Characterization. Anal. Lett. 2011, 44, 2856–2864. [Google Scholar] [CrossRef]
  6. Ehlers, M.; Uttl, L.; Riedl, J.; Raeke, J.; Westkamp, I.; Hajslova, J.; Brockmeyer, J.; Fauhl-Hassek, C. Instrument comparability of non-targeted UHPLC-HRMS for wine authentication. Food Control 2022, 144, 109360. [Google Scholar] [CrossRef]
  7. Pan, Y.; Gu, H.-W.; Lv, Y.; Yin, X.-L.; Chen, Y.; Long, W.; Fu, H.; She, Y. Untargeted metabolomic analysis of Chinese red wines for geographical origin traceability by UPLC-QTOF-MS coupled with chemometrics. Food Chem. 2022, 394, 133473. [Google Scholar] [CrossRef]
  8. Yunoki, K.; Yasui, Y.; Naito, A.; Kawaguchi, M.; Ohnishi, M. Comparative analysis of concentrations of trans- and cis-resveratrol in must and wine from eleven different grapes grown in Hokkaido, Northen Japan. J. Oleo Sci. 2001, 50, 895–903. [Google Scholar] [CrossRef]
  9. Rodríguez-Delgado, M.Á.; González-Hernández, G.; Conde-González, J.E.; Pérez-Trujillo, J.P. Principal component analysis of the polyphenol content in young red wines. Food Chem. 2002, 78, 523–532. [Google Scholar] [CrossRef]
  10. Rastija, V.; Srecnik, G.; Medic-Saric, M. Polyphenolic composition of Croatian wines with different geographical origin. Food Chem. 2009, 115, 54–60. [Google Scholar] [CrossRef]
  11. Leeuwen, C.; Destrac-Irvine, A.; Dubernet, M.; Duchêne, E.; Gowdy, M.; Marguerit, E.; Pieri, P.; Parker, A.; de Rességuier, L.; Ollat, N. An Update on the Impact of Climate Change in Viticulture and Potential Adaptations. Agronomy 2019, 9, 514. [Google Scholar] [CrossRef]
  12. Zaldea, G.; Nechita, A.; Damian, D.; Ghiur, A.D.; Cotea, V.V. Climate changes in recent decades, the evolution of the drought phenomenon and their influence on vineyards in north-eastern Romania. Not. Bot. Horti Agrobot. 2021, 49, 12448. [Google Scholar] [CrossRef]
  13. Nedelcu, A.; Privitera, D.; Ivona, A.; Ganusceac, A. Wine tourism as a vector of local and regional development. Case study Prahova county. Udekom Balcan 3rd Int. Themat. Monograph. 2021, 341–357. [Google Scholar] [CrossRef]
  14. Liviu, M.; Nicolae, G.; Carmen, P.; Ramona, C. Tradition and quality of romanian viticulture. Ann. Univ. Craiova 2018, 23, 139–143. [Google Scholar]
  15. Irimia, L.M.; Patriche, C.V.; Quenol, H.; Sfîcă, L.; Foss, C. Shifts in climate suitability for wine production as a result of climate change in a temperate climate wine region of Romania. Theor. Appl. Climatol. 2018, 131, 1069–1081. [Google Scholar] [CrossRef]
  16. Nedelcu, A. Wine Tourism in Romania: Case Study: Dealu Mare Vineyard. Int. J. Sustain. Econ. Manag. 2014, 3, 16–26. [Google Scholar] [CrossRef]
  17. Ballester, J.; Magne, M.; Julien, P.; Noret, L.; Nikolantonaki, M.; Coelho, C.; Gougeon, R.D. Sensory Impact of Polyphenolic Composition on the Oxidative Notes of Chardonnay Wines. Beverages 2018, 4, 19. [Google Scholar] [CrossRef]
  18. Guyon, F.; van Leeuwen, C.; Gaillard, L.; Grand, M.; Akoka, S.; Remaud, G.S.; Sabathié, N.; Salagoïty, M.H. Comparative study of 13C composition in ethanol and bulk dry wine using isotope ratio monitoring by mass spectrometry and by nuclear magnetic resonance as an indicator of vine water status. Anal. Bioanal. Chem. 2015, 407, 9053–9060. [Google Scholar] [CrossRef]
  19. Šorgić, S.; Sredović Ignjatović, I.; Antić, M.; Šaćirović, S.; Pezo, L.; Čejić, V.; Đurović, S. Monitoring of the Wines’ Quality by Gas Chromatography: HSS-GC/FID Method Development, Validation, Verification, for Analysis of Volatile Compounds. Fermentation 2022, 8, 38. [Google Scholar] [CrossRef]
  20. Dinca, O.R.; Ionete, R.E.; Costinel, D.; Geana, I.E.; Popescu, R.; Stefanescu, I.; Radu, G.L. Regional and Vintage Discrimination of Romanian Wines Based on Elemental and Isotopic Fingerprinting. Food Anal. Methods 2016, 9, 2406–2417. [Google Scholar] [CrossRef]
  21. Banc, R.; Loghin, F.; Miere, D.; Ranga, F.; Socaciu, C. Phenolic composition and antioxidant activity of red, rosé and white wines originating from Romanian grape cultivars. Not. Bot. Horti Agrobot. 2020, 48, 716–734. [Google Scholar] [CrossRef]
  22. Geana, E.-I.; Ionete, R.E.; Tudorache, A.; Pasa, R.; Postolache, E.; Ranca, A. Phenolic Contents of Romanian Wines with Different Geographical Origins. Asian J. Chem. 2011, 23, 5. [Google Scholar]
  23. Pajović-Šćepanović, R.; Wendelin, S.; Eder, R. Phenolic composition and varietal discrimination of Montenegrin red wines (Vitis vinifera var. Vranac, Kratošija, and Cabernet Sauvignon). Eur. Food Res. Technol. 2018, 244, 2243–2254. [Google Scholar] [CrossRef]
  24. Geana, E.I.; Popescu, R.; Costinel, D.; Dinca, O.R.; Stefanescu, I.; Ionete, R.E. Verifying the Red Wines Adulteration through Isotopic and Chromatographic Investigations Coupled with Multivariate Statistic Interpretation of the Data. Food Control 2016, 62, 1–9. [Google Scholar] [CrossRef]
  25. Tzachristas, A.; Pasvanka, K.; Calokerinos, A.; Proestos, C. Polyphenols: Natural Antioxidants to be used as a Quality Tool in Wine Authenticity. Appl. Sci. 2020, 10, 5908. [Google Scholar] [CrossRef]
Figure 1. Normal distribution of variables within each group (Fetească neagră, Busuioacă de Bohotin and Fetească alba).
Figure 1. Normal distribution of variables within each group (Fetească neagră, Busuioacă de Bohotin and Fetească alba).
Agronomy 12 02286 g001
Table 1. Physico-chemical characteristics of the studied wines.
Table 1. Physico-chemical characteristics of the studied wines.
Current NumberRegionGrape VarietyWine ColourYearDensity (g/cm3)Ethanol Concentration (% v/v)pHTotal Acidity
(g/L Tartaric Acid)
Sugars (g/L)
1CotnariFetească neagrăRed20170.993413.263.635.575.90
220170.992013.663.724.974.80
320170.993713.083.645.544.90
420170.993813.373.575.796.10
520170.993112.163.735.025.10
6Busuioacă de BohotinRosé20190.989713.633.156.613.70
720191.001212.793.356.3829.50
820190.993712.773.366.1411.40
920190.993512.913.376.1411.20
1020190.993412.983.306.2812.10
1120190.990313.053.266.643.90
1220190.993513.183.345.9112.10
13Fetească albăWhite20180.988512.963.444.852.60
1420181.001812.533.256.1131.50
1520180.990113.193.276.154.20
1620180.992013.523.166.6310.30
17Dealu mareFetească neagrăRed20170.990013.613.445.494.80
1820170.992414.413.645.275.60
1920170.993814.253.835.045.80
2020170.991714.103.814.873.40
21Busuioacă de BohotinRosé20190.990112.123.456.721.40
2220191.005312.303.665.4835.20
2320191.004612.513.107.3638.40
2420190.990813.413.396.285.80
25Fetească albăWhite20180.990511.493.415.493.40
2620180.990013.613.445.494.80
2720180.994713.003.475.7912.30
2820180.988512.963.444.852.60
Table 2. The HPLC gradient elution.
Table 2. The HPLC gradient elution.
White and Rosé WinesRed Wines
Time (min)Eluent B (%)C (%)Time (min)Eluent B (%)C (%)
53970–83–597–95
5–83–597–958–17595
8–1759517–195–995–91
17–195–995–9119–25991
19–2599125–359–14.391–85.7
25–259–14.391–85.735–3814.385.7
35–3614.3–2385.7–6738–4123–2777–73
36–4523–10067–041–4627–5073–50
35–50100046–5250–10050–0
50–55100–30–9752–571000
55–6039757–62100–30–97
62–68397
Table 3. Carbon isotope ratios (δ13C).
Table 3. Carbon isotope ratios (δ13C).
Current NumberRegionGrape Varietyδ13C VPDB (‰)SD
1CotnariFetească neagră−27.300.16
2−26.890.02
3−27.480.16
4−27.360.16
5−27.040.02
6Busuioacă de Bohotin−25.660.06
7−26.130.09
8−26.210.27
9−26.350.00
10−26.440.20
11−26.340.17
12−26.480.06
13Fetească albă−26.120.39
14−26.370.11
15−27.270.11
16−26.630.25
17Dealu MareFetească neagră−27.830.07
18−28.120.00
19−28.180.22
20−27.920.05
21Busuioacă de Bohotin−27.130.25
22−27.140.09
23−26.980.04
24−27.160.08
25Fetească albă−26.740.25
26−26.280.05
27−26.870.14
28−27.200.09
Table 8. Chi-square tests with successive roots removed.
Table 8. Chi-square tests with successive roots removed.
Roots
(Function)
Eigen ValueCanonic RWilks’ LambdaChi-Sqrdfp-Value
1151.44770.9967150.000536143.1067260.000001
211.24490.9582970.08166747.5970120.000004
Table 9. Standardized coefficients for canonical variables.
Table 9. Standardized coefficients for canonical variables.
VariableRoot 1Root 2
Gallic acid−0.71470.28852
Protocatechuic acid0.32250.36014
Caftaric acid1.0130−0.25608
Caffeic acid−2.6101−0.15378
Coumaric acid0.2592−1.23542
Trans-resveratrol−0.73781.36782
Hydroxytyrosol0.89150.18484
Tyrosol0.39682.06080
Procyanidin dimer B1−0.3499−6.83465
Catechin0.23031.35097
Procyanidin dimer B2−0.62675.42857
Epicatechin−3.1831−1.36891
δ13C mean−0.3352−0.32343
Eigenval151.44711.24487
Cum. Prop.0.93091.00000
Red—variables that discriminates with statistical significance (p < 0.05) the first discriminant function (root 1). Green—variables that discriminates with statistical significance (p < 0.05) second discriminant function (root 2).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Popîrdă, A.; Luchian, C.E.; Colibaba, L.C.; Focea, E.C.; Nicolas, S.; Noret, L.; Cioroiu, I.B.; Gougeon, R.; Cotea, V.V. Carbon-Isotope Ratio (δ13C) and Phenolic-Compounds Analysis in Authenticity Studies of Wines from Dealu Mare and Cotnari Regions (Romania). Agronomy 2022, 12, 2286. https://doi.org/10.3390/agronomy12102286

AMA Style

Popîrdă A, Luchian CE, Colibaba LC, Focea EC, Nicolas S, Noret L, Cioroiu IB, Gougeon R, Cotea VV. Carbon-Isotope Ratio (δ13C) and Phenolic-Compounds Analysis in Authenticity Studies of Wines from Dealu Mare and Cotnari Regions (Romania). Agronomy. 2022; 12(10):2286. https://doi.org/10.3390/agronomy12102286

Chicago/Turabian Style

Popîrdă, Andreea, Camelia Elena Luchian, Lucia Cintia Colibaba, Elena Cornelia Focea, Sebastien Nicolas, Laurence Noret, Ionel Bogdan Cioroiu, Régis Gougeon, and Valeriu V. Cotea. 2022. "Carbon-Isotope Ratio (δ13C) and Phenolic-Compounds Analysis in Authenticity Studies of Wines from Dealu Mare and Cotnari Regions (Romania)" Agronomy 12, no. 10: 2286. https://doi.org/10.3390/agronomy12102286

APA Style

Popîrdă, A., Luchian, C. E., Colibaba, L. C., Focea, E. C., Nicolas, S., Noret, L., Cioroiu, I. B., Gougeon, R., & Cotea, V. V. (2022). Carbon-Isotope Ratio (δ13C) and Phenolic-Compounds Analysis in Authenticity Studies of Wines from Dealu Mare and Cotnari Regions (Romania). Agronomy, 12(10), 2286. https://doi.org/10.3390/agronomy12102286

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

Article Metrics

Back to TopTop