# 1H-NMR Metabolomics as a Tool for Winemaking Monitoring

Inès Le Mao, Jean Martin-Pernier, Charlyne Bautista, Soizic Lacampagne, Tristan Richard, Gregory da Costa

► To cite this version:

Inès Le Mao, Jean Martin-Pernier, Charlyne Bautista, Soizic Lacampagne, Tristan Richard, et al. 1H-NMR Metabolomics as a Tool for Winemaking Monitoring. Molecules, 2021, 26 (22), pp.6771. 10.3390/molecules26226771. hal-03507798

HAL Id: hal-03507798

https://hal.inrae.fr/hal-03507798

Submitted on 3 Jan 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# ¹H-NMR Metabolomics as a Tool for Winemaking Monitoring

Inès Le Mao, Jean Martin-Pernier, Charlyne Bautista, Soizic Lacampagne, Tristan Richard * and Gregory Da Costa

University of Bordeaux, INRAE, Bordeaux INP, UR OENO, EA 4577, USC 1366, F-33140 Villenave d'Ornon, France; ines.le-mao@u-bordeaux.fr (I.L.M.); jean.martin@u-bordeaux.fr (J.M.-P.); charlynebautista@hotmail.com (C.B.); soizic.lacampagne@u-bordeaux.fr (S.L.); gregory.da-costa@u-bordeaux.fr (G.D.C.)

* Correspondence: tristan.richard@u-bordeaux.fr; Tel.: +33-557575957

**Abstract:** The chemical composition of wine is known to be influenced by multiple factors including some viticulture practices and winemaking processes. ¹H-NMR metabolomics has been successfully applied to the study of wine authenticity. In the present study, ¹H-NMR metabolomics in combination with multivariate analysis was applied to investigate the effects of grape maturity and enzyme and fining treatments on Cabernet Sauvignon wines. A total of forty wine metabolites were quantified. Three different stages of maturity were studied (under-maturity, maturity and over-maturity). Enzyme treatments were carried out using two pectolytic enzymes (E1 and E2). Finally, two proteinaceous fining treatments were compared (vegetable protein, fining F1; pea protein and PVPP, fining F2). The results show a clear difference between the three stages of maturity, with an impact on different classes of metabolites including amino acids, organic acids, sugars, phenolic compounds, alcohols and esters. A clear separation between enzymes E1 and E2 was observed. Both fining agents had a significant effect on metabolite concentrations. The results demonstrate that ¹H-NMR metabolomics provides a fast and robust approach to study the effect of winemaking processes on wine metabolites. These results support the interest to pursue the development of ¹H-NMR metabolomics to investigate the effects of winemaking on wine quality.

**Keywords:** ¹H-NMR; metabolomics; wine; winemaking

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Le Mao, I.; Martin-Pernier, J.; Bautista, C.; Lacampagne, S.; Richard, T.; Da Costa, G. ¹H-NMR Metabolomics as a Tool for Winemaking Monitoring. *Molecules* **2021**, *26*, 6771. https://doi.org/10.3390/molecules26226771

Academic Editor: Fulvio Mattivi

Received: 13 October 2021

Accepted: 5 November 2021

Published: 9 November 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

## 1. Introduction

Over the last decade, the number of studies on wine has grown steadily as quality and authenticity have become major concerns for both producers and consumers. Thus, fine wine analysis has blossomed to study and characterize each step of the production chain from the vineyard to the finished wine. Metabolomics has emerged from the studies that have flourished, which studies all the small molecules of a biological system. Metabolomics science, by adopting a targeted or non-targeted approach, can be discriminating, predictive and informative depending on the subject of the study [1,2].

A well-known challenge in the study of wine is the complexity of the matrix. Indeed, hundreds of compounds of different chemical families make analysis complex. The two most commonly used analytical techniques for wine metabolomics are mass spectrometry and NMR [3]. Although MS is a more sensitive tool, NMR combines many advantages that make it a tool of choice for metabolomic analysis of wines [4–7]: no complex sample preparation, low volume requirements, fast and reproducible analysis and detection of a large number of compounds belonging to different chemical families, such as organic acids, alcohols, sugars, phenolic compounds and amino acids [8,9]. Moreover, NMR is a wine screening tool that can also provide quantitative data thanks to [10].

Since the chemical composition of wine is known to be influenced by multiple factors, many wine NMR studies have been focused on the impact of grape variety [11,12], geographical origin [9,13,14], vintage and aging [15,16]. Comparatively, few NMR technique-based studies have been carried out on viticultural practices and winemaking processes,
which can also influence the metabolite profile of wine, and thus its organoleptic proper-
ties [17,18]. Fermentation processing was explored by establishing the difference between
alcoholic and malolactic fermentation processes on metabolites with a time course evolu-
tion [19], while other studies focused on the impact of using different yeast strains during
white wine fermentation [20], and red wine fermentation [11]. Lee and collaborators
established, by a combination of ¹H-NMR and gas chromatography approaches, the fer-
mentative behaviors and the metabolic variations of five bacteria strains during malolactic
fermentation in Meoru wines, in order to reduce wine acidity [15]. Concerning winemak-
ing technologies, through metabolite profiles, De Pascali and collaborators analyzed the
difference between winemaking processes (ultrasound, cryomaceration and traditional
processes) on Negroamaro red wine [21], while Baiano and collaborators examined the
difference between traditional, cryomaceration and reductive processes and a combination
of the last two procedures on Sauvignon Blanc wine [22].

Wines are subjected to different winemaking processes including the addition of different agents for better wine stability and limpidity, such as finings, used to obtain clarified and limpid wines [23], and enzymatic treatments, used to enhance the degrada- tion of grape berry cell walls and the extraction of aroma and phenolic compounds [24]. The control of these different processes is a crucial point in order to produce quality wines. The impact of finings on wine composition has been studied using different techniques including spectrophotometry or liquid and gas chromatography [23,25]. In the same way, different techniques have been used to study the effect of enzymatic addition on the re- lease of phenols, monosaccharides and polysaccharides, by different chromatographic techniques [24,26].

The aim of this study was to determine the ability of ¹H-NMR spectroscopy combined
with multivariate statistical analyses to discriminate different winemaking parameters and
processes (maturity, fining and enzymatic addition) and to explore the impact of these
processes on wine composition.

**2. Results and Discussion**

2.1. ¹H-NMR Analysis of Wine

The typical ¹H-NMR spectra of wine after solvent suppression are presented in Figure 1. The signals at 0.00 and 8.28 ppm correspond to TMSP and calcium formate, respectively. All other signals correspond to wine constituents. Based on previous studies [8,17], 2D NMR analyses and pure chemical standards, forty-eight constituents were identified (Table 1), including amino acids, organic acids, alcohols, carbohydrates, phenolic compounds, esters, ketones and aldehydes.

## 2.2. *Maturity Stages*

The ability of ¹H-NMR-based metabolomics to discriminate wines produced from grapes harvested at three different stages of maturity (M1: under-maturity; M2: maturity; M3: over-maturity) was assessed using multivariate analysis. Principal component analysis (PCA) followed by orthogonal partial least square discriminant analysis (OPLS-DA) was performed on the quantification data extracted from the ¹H-NMR spectra. The analysis of the maturity effect was carried out using six replicates for each maturity stage.

The PCA score plot (R²X(cum) and Q²(cum) of 0.480 and 0.130, respectively) showed separation between the three stages of maturity (Figure 2a). To confirm this observation, OPLS-DA (number of latent variables 2 + 2 + 0) was performed (Figure 2b). The OPLS-DA score plots confirmed discrimination between wines with high validated predictabil- ity (Q² = 0.873) and goodness of fit values (R²X and R²Y of 0.615 and 0.975, respec- tively). The validity of the OPLS-DA model was confirmed by cross-validated ANOVA (CV-ANOVA) providing a significant *p*-value (3.12 × 10⁻⁵). The recognition ability of the model was 100%, as shown in the misclassification table (Table S2). The under-maturity, ma- turity and over-maturity stages were completely distinguished from one another. This first result indicates that ¹H-NMR metabolomics is able to discriminate wines made from berries
collected at different stages of maturity. To analyze the effect of under-mature and over-
mature grapes on wine quality, the OPLS-DA model was further applied to discriminate
the differential constituents between these wines and those made with grapes harvested at
maturity. The OPLS-DA score and loading plots are shown in Figure 2c,d, respectively, for
comparison of under-maturity/maturity (M1/M2), and in Figure 2e,f, respectively, for com-
parison of over-maturity/maturity (M3/M2). The OPLS-DA model parameters (Table S3)
and cross-validation procedure demonstrated an excellent modeling and predictive ability
in each case.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Typical ¹H-NMR spectrum of wine after water and ethanol suppression (NOESYGPPS1D sequence). Identified constituents are listed in Table 1 (compounds in green are quantified on the ZGPR sequence).</figcaption></figure>
<table><caption>**Table 1.** Chemical shifts and coupling constants used for compound identification. The signals chosen for quantitation are in bold.</caption><thead><tr><th></th><th>Compound</th><th>δ<sup>1</sup>H (Multiplicity, J in Hz, Assignment)</th></tr></thead><tbody><tr><td>1</td><td>leucine</td><td>0.96 (d, 6.2, 2CH<sub>3</sub>), 1.71 (m, CHCH<sub>2</sub>), 3.74 (m, CH)</td></tr><tr><td>2</td><td>isoleucine</td><td>0.93 (t, 7.4, CH<sub>3</sub>), 0.99 (d, 7.0, CH<sub>3</sub>), 1.24 (m, CH<sub>2</sub>), 1.45 (m, CH<sub>2</sub>), 1.97 (m, CH), 3.66 (d, 3.9, CH)</td></tr><tr><td>3</td><td>valine</td><td>0.99 (d, 7.3, CH<sub>3</sub>), 1.04 (d, 7.3, CH<sub>3</sub>), 2.28 (m, CH), 3.66 (d, 4.3, CH)</td></tr><tr><td>4</td><td>2,3-butanediol</td><td>1.13 (d, 6.2, 2CH<sub>3</sub>), 3.61 (m, 2CH)</td></tr><tr><td>5</td><td>ethanol</td><td>1.17 (t, 7.2, CH<sub>3</sub>), 3.65 (q, CH<sub>2</sub>)</td></tr><tr><td>6</td><td>threonine</td><td>1.32 (d, 6.7, CH<sub>3</sub>), 2.58 (d, 4.9, CH), 4.24 (m, CH)</td></tr><tr><td>7</td><td>acetoin</td><td>1.37 (d, 7.0, CH<sub>3</sub>), 2.21 (s, CH<sub>3</sub>), 4.42 (q, CH)</td></tr><tr><td>8</td><td>lactic acid</td><td>1.40 (d, 7.0, CH<sub>3</sub>), 4.31 (q, 7.0, CH)</td></tr><tr><td>9</td><td>alanine</td><td>1.50 (d, 7.2, CH<sub>3</sub>), 3.76 (q, CH)</td></tr><tr><td>10</td><td>isopentanol</td><td>0.88 (d, 6.7, 2CH<sub>3</sub>), 1.44 (q, CH<sub>2</sub>); 1.66 (m, CH), 3.61 (t, 6.7, CH<sub>2</sub>)</td></tr><tr><td>11</td><td>arginine</td><td>1.70 (m, CH<sub>2</sub>), 1.89 (m, CH<sub>2</sub>), 3.23 (t, CH<sub>2</sub>), 3.75 (t, 6.5, CH)</td></tr><tr><td>12</td><td>proline</td><td>1.99 (m, CH<sub>2</sub>), 2.06 (m, CH), 2.33 (m, CH), 3.32 (dt, 14.0, 7.1, CH), 3.42 (dt, 11.6 and 7.0, CH), 4.11 (dd, 8.6 and 6.4, CH)</td></tr><tr><td>13</td><td>ethyl acetate</td><td>1.26 (t, 7.2, CH<sub>3</sub>), 4.12 (q, CH<sub>2</sub>), 2.07 (s, CH<sub>3</sub>)</td></tr><tr><td>14</td><td>acetic acid</td><td>2.08 (s, CH<sub>3</sub>)</td></tr><tr><td>15</td><td>ethanal</td><td>2.23 (d, 3.0, CH<sub>3</sub>), 9.67 (q, CH)</td></tr><tr><td>16</td><td>pyruvic acid</td><td>2.35 (s, CH<sub>3</sub>)</td></tr><tr><td>17</td><td>γ-aminobutyric acid</td><td>1.96 (m, CH<sub>2</sub>), 2.50 (t, 7.3, CH<sub>2</sub>), 3.05 (m, CH<sub>2</sub>)</td></tr><tr><td>18</td><td>succinic acid</td><td>2.65 (s, 2CH<sub>2</sub>)</td></tr><tr><td>19</td><td>malic acid</td><td>2.78 (dd, 16.3 and 7.0, CH), 2.89 (dd, 16.3 and 4.5, CH), 4.53 (dd, CH)</td></tr><tr><td>20</td><td>citric acid</td><td>2.79 (d, 15.6, CH<sub>2</sub>), 2.94 (d, 15.6, CH<sub>2</sub>)</td></tr><tr><td>21</td><td>choline</td><td>3.19 (s, 3CH<sub>3</sub>), 3.51 (dd, CH<sub>2</sub>), 4.05 (m, CH<sub>2</sub>)</td></tr><tr><td>22</td><td>myo-inositol</td><td>3.27 (t, 9.7, CH), 3.52 (dd, 10.0 and 2.8, 2CH), 3.61 (t, 2.8, 2CH), 4.05 (t, 2.8, CH)</td></tr><tr><td>23</td><td>methanol</td><td>3.35 (s, CH<sub>3</sub>)</td></tr><tr><td>24</td><td>isobutanol</td><td>0.87 (d, 6.7, 2CH<sub>3</sub>), 1.73 (m, CH), 3.36 (d, 6.7, CH<sub>2</sub>)</td></tr><tr><td>25</td><td>glycerol</td><td>3.55 (dd, 11.8 and 6.5, CH<sub>2</sub>), 3.64 (dd, CH<sub>2</sub>), 3.77 (m, CH)</td></tr><tr><td>26</td><td>mannitol</td><td>3.65 (dd, 11.7, 6.2 CH<sub>2</sub>), 3.73 (m, CH), 3.77 (d, 9.0, CH), 3.84 (dd, 11.9, 2.8, CH<sub>2</sub>)</td></tr><tr><td>27</td><td>fructose</td><td>3.56 (m, CH<sub>2</sub>), 3.70 (m, 2CH<sub>2</sub>), 3.77 (m, CHCH<sub>2</sub>), 3.87 (dd, 9.9, 3.4, CH), 3.97 (m, CH), 4.00 (dd, 12.8, 1.0 CH<sub>2</sub>), 4.09 (m, 2CH)</td></tr><tr><td>28</td><td>ethyl lactate</td><td>1.28 (t, CH<sub>3</sub>), 1.42 (d, 7.0, CH<sub>3</sub>), 4.22 (q, 7.06, CH), 4.39 (q, 7.0, CH)</td></tr><tr><td>29</td><td>arabinose</td><td>3.51 (dd, CH), 3.68 (m, CHCH<sub>2</sub>), 3.83 (dd, CH), 3.90 (m, CHCH<sub>2</sub>), 3.95(m, CH), 4.02 (m, CHCH<sub>2</sub>), 4.50 (d, 7.7, CH), 5.25 (d, CH)</td></tr><tr><td>30</td><td>glucose</td><td>3.23 (dd, 9.2, 8.0, CH), 3.39 (m, CH), 3.45 (dd, 9.8, 3.7, CH), 3.72 (m, CHCH<sub>2</sub>), 3.82 (m, CHCH<sub>2</sub>), 3.88 (dd, 12.2, 2.1, CH<sub>2</sub>), 4.63 (d, 7.9, CH), 5.22 (d, 3.6, CH)</td></tr><tr><td>31</td><td>tartaric acid</td><td>4.60 (s, 2CH)</td></tr><tr><td>32</td><td>xylose</td><td>3.21 (dd, 9.3, 7.9, CH), 3.31 (t, 11.4, CH<sub>2</sub>), 3.42 (t, 9.25, CH), 3.51 (dd, 9.3, 3.7, CH), 3.63 (m, CHCHCH<sub>2</sub>), 3.91 (dd, 11.5, 5.5, CH<sub>2</sub>), 4.57 (d, 7.9, CH), 5.19 (d, 3.7, CH)</td></tr><tr><td>33</td><td>galacturonic acid</td><td>3.49 (dd, 8.0, 10.0, CH), 3.69 (dd, 9, 3.5, CH), 3.80 (dd, 10.3, 3.8, CH), 3.92 (dd, 10.3, 3.4, CH), 4.24 (dd, 3.6, 1.2, CH), 4.26 (d, 1.2, CH), 4.31 (dd, 3.3, 1.4, CH), 5.32 (d, 3.8, CH)</td></tr><tr><td>34</td><td>glucuronic acid</td><td>3.29 (t, 8.6, CH), 3.51 (m, 2CH), 3.58 (dd, 9.7, 3.7, CH), 3.73 (m, 2CH), 4.08 (d, 10.8, CH), 4.64 (d, 7.9, CH), 5.25 (d, 3.7, CH), 5.55 (d, 4, CH)</td></tr><tr><td>35</td><td>sorbic acid</td><td>1.82 (d, 6.2, CH<sub>3</sub>), 5.78 (d, 15.3, CH), 6.25 (m, 2CH), 7.16 (dd, 15.3, 10.3, CH)</td></tr><tr><td>36</td><td>epicatechin</td><td>2.76 (m, CH<sub>2</sub>), 2.90 (m, CH<sub>2</sub>), 4.32 (m, CH), 4.95 (m, CH), 6.09 (d, 2.0, CH), 6.12 (d, 2.0, CH), 6.93 (m, CH<sub>2</sub>), 7.03 (d, 2.0, CH)</td></tr><tr><td>37</td><td>catechin</td><td>2.53 (dd, CH<sub>2</sub>), 2.85 (m, CH<sub>2</sub>), 4.15 (m, CH), 4.41 (d, 7.0, CH), 5.99 (d, 2.0, CH), 6.08 (d, 2.3, CH), 6.84 (d, 8.6, CH), 6.92 (m, 2CH)</td></tr><tr><td>38</td><td>caffeic acid</td><td>6.33 (d, 16.0, CH), 6.92 (d, 8.0, CH), 7.07 (dd, 8.2, 2.0, CH), 7.14 (d, 2.0, CH), 7.29 (d, CH)</td></tr><tr><td>39</td><td>fumaric acid</td><td>6.78 (s, 2CH)</td></tr><tr><td>40</td><td>shikimic acid</td><td>2.21 (dd, 18.2, 7.0, CH<sub>2</sub>), 2.75 (dd, 18.0, 5.3, CH<sub>2</sub>), 3.74 (dd, 8.6, 4.3, CH), 4.01 (m, CH), 4.42 (t, 4.1, CH), 6.82 (dt, CH)</td></tr><tr><td>41</td><td>tyrosol</td><td>2.77 (t, CH<sub>2</sub>), 3.77 (t, CH<sub>2</sub>), 6.84 (m, 8.4, 2CH), 7.17 (m, 8.4, 2CH)</td></tr><tr><td>42</td><td>tyrosine</td><td>3.02 (dd, CH<sub>2</sub>), 3.17 (dd, CH<sub>2</sub>), 3.92 (dd, CH), 6.86 (m, 8.4, 2CH), 7.17 (m, 8.6, 2CH)</td></tr><tr><td>43</td><td>gallic acid</td><td>7.16 (s, 2CH)</td></tr><tr><td>44</td><td>phenethyl alcohol</td><td>2.85 (t, 6.62, CH<sub>2</sub>), 3.74 (t, CH<sub>2</sub>), 7.33 (m, 5CH)</td></tr><tr><td>45</td><td>syringic acid</td><td>3.84 (s, 2CH<sub>3</sub>), 7.36 (s, 2CH)</td></tr><tr><td>46</td><td>histidine</td><td>3.16 (dd, 15.6, 7.7, CH), 3.23 (dd, 16.0, 5.0, CH), 3.98 (dd, 7.7, 5.0, CH), 7.09 (d, 5.0, CH), 7.90 (d, 1.1, CH)</td></tr><tr><td>47</td><td>trigonelline</td><td>4.42 (s, CH<sub>3</sub>), 8.07 (m, CH), 8.82 (m, 2CH), 9.11 (s, CH)</td></tr></tbody></table>
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Multivariate analysis of ¹H-NMR spectra of wine samples from grapes harvested at three different stages of maturity (M1: under-maturity; M2: maturity; M3: over-maturity): (a) PCA score plot; (b) OPLS-DA score plot; (c) OPLS-DA score showing separation of M1 and M2 samples; (d) loadings from OPLS-DA between M1 and M2 samples; (e) OPLS-DA score showing separation of M2 and M3 samples; (f) loadings from OPLS-DA between M2 and M3 samples (t[1] and to[1]: first predictive and orthogonal components; t[2]: second predictive component; pq[1] and poso[1]: predictive and orthogonal component loadings).</figcaption></figure>

To identify the contribution of wine constituents to the discrimination between different groups, a combination of variable importance in projection (VIP) and p corr values was used [27]. Potential differential compounds were selected with absolute VIP and p corr values greater than 1.0 and 0.5, respectively. The highlighted wine constituents were further controlled by ANOVA followed by Tukey's multiple comparison test, and only those with a *p*-value of <0.05 were conserved. A total of 16 wine constituents including amino acids (alanine, histidine, leucine and tyrosine), organic acids (citric, lactic, tartaric and shikimic acids), sugar derivatives (glucose, galacturonic acids and xylose), phenolic compounds (catechin, epicatechin and gallic acid), an alcohol (methanol) and an ester (ethyl lactate) had a strong influence in distinguishing M2 from M1 wines. Wines made with grapes harvested at under-maturity contained less alanine, histidine, leucine, tyrosine, citric acid, tartaric acid, glucose, galacturonic acid, gallic acid, ethyl lactate and methanol (Figure S1). On the contrary, they contained more lactic acid, shikimic acid, xylose, catechin and epicatechin than wines made with grapes harvested at maturity. Similarly, a total of 10 wine constituents including organic acids (citric acid, lactic acid and shikimic acid), sugar derivatives (glucose, galacturonic acids and xylose) and alcohols (acetoin, 2,3-butanediol, glycerol and methanol) had a strong influence in separating M3 and M2. Wines made with grapes harvested at over-maturity contained more citric acid, glucose, galacturonic acid, 2,3-butanediol, glycerol and methanol (Figure S2). On the contrary, they contained
less lactic acid, shikimic acid, xylose and acetoin than wines made with grapes harvested at maturity.

Most studies explored the ripening effect on the berries' metabolome [28], but few focused on the wine metabolome. Nevertheless, grape maturity generated significant sensory effects on wines [29]. In this study, a significant effect of grape maturity was observed on the Cabernet Sauvignon wine metabolome using ¹H-NMR metabolomics.

The overall trends observed are in agreement with the results obtained by Bindon et al. on Cabernet Sauvignon [30]. Even though the ethanol concentration increased from 13.7 at M1 to 14.2% at M3, this variation was not significant. Nevertheless, a significant increase in other alcohols such glycerol was observed, particularly between M2 and M3 (Figure S2). As expected, the higher lactic acid concentrations in M1 (mean 1.90 ± 0.09 g/L) samples can be explained by the higher initial malic acid concentrations in grapes harvested at under-ripeness (Table S1). Indeed, as the finished wines underwent malolactic fermentation, the initial malate was partly transformed into lactate by lactic acid bacteria [31]. On the contrary, a significant increase in tartaric acid was observed between M1 and M2 samples. This result could be associated with the increase in tartaric acid in grapes between under-maturity and maturity (Table S1). Regarding the other acids, no significant variation was observed except for citric acid and shikimic acid. This observation is in agreement with literature data [30]. Practically all reducing sugars (glucose and fructose) were consumed, indicating the efficiency of fermentation. A higher content was observed in M3 samples (mean 0.89 ± 0.15 g/L for glucose and fructose). As previously reported [32], proline was the main amino acid observed in wines (mean range 0.801–0.962 g/L), and the highest concentration was measured in M3 samples but with no significant difference from other modalities. In addition, a significant increase in four amino acids was observed between M1 and M2 samples (Figure S3). Interestingly, some specific trends were observed such as the significant variation in acetoin between M2 and M3 samples (means 23 ± 2 and 11 ± 3 mg/L, respectively). This compound is a well-known precursor of the biosynthesis of 2,3-butanediol, which contributes to the bouquet of the wine. To conclude, ¹H-NMR metabolomics was able to discriminate wines produced from grapes at different stages of ripeness, providing informative trends on the potential quality of the wine.

## 2.3. *Enzyme Treatments*

Enzymes are used to facilitate hydrolysis of the cell walls, providing a more efficient extraction of phenolic compounds [33]. The effect of two different preparations, E1 and E2, was monitored using ¹H-NMR-based metabolomics and multivariate analysis. Both preparations are pectolytic enzymes. Five replicates for each enzyme were analyzed.

Wine samples were treated by enzymes E1 and E2. The data extracted from ¹H-NMR spectra were compared by PCA and OPLS-DA. The PCA score plot exhibited separation between the three modalities (Figure S3a). To better emphasize the discrimination, OPLS-DA was performed. Score and loading plots are shown in Figure 3a,b, respectively. A clear tendency to discriminate the enzyme treatments was observed (Q², R²X and R²Y of 0.617, 0.869 and 0.518, respectively). To validate the model, a cross-validation procedure was performed, which indicated that all wines were correctly classified (Table S4). ¹H NMR-based metabolomics allows observing the effect of the use of enzymes independent of the grape maturity stage. To analyze the effect of each enzymatic preparation, the OPLS-DA model was further applied to discriminate the differential constituents between these wines and untreated samples. Discriminating constituents were identified based on VIP and absolute p corr) values and confirmed by ANOVA followed by Tukey's multiple comparison test. A total of 11 compounds were identified (Figure 3c).
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Multivariate analysis of ¹H-NMR spectra of wine samples treated by different enzymes (E0: untreated; E1: enzyme 1; E2: enzyme 2): (a) OPLS-DA score plot; (b) loading plot; (c) boxplots of 11 most discriminant wine constituents (t[1] and t[2]: first and second predictive components; pq[1] and pq[2]: first and second predictive component loadings). The significance in the difference was calculated by ANOVA followed by Tukey's multiple comparison test (indicated as \* *p* < 0.05, \*\* *p* < 0.01, \*\*\* *p* < 0.001).</figcaption></figure>

Wine samples treated by enzymes contained significantly more arabinose, galacturonic acid, glucose and catechin and less xylose, myo-inositol and succinic acid. As expected, the enzyme treatments had a significant effect on wine polysaccharide composition [24]. Pectolytic enzymes reduce the molecular weight of polysaccharides, resulting in an increase in free monosaccharides. Significant increases in galacturonic acid, arabinose and glucose have previously been reported in wine samples treated by enzyme preparations [34]. These are clear indicators of the pectolytic activity of enzymes resulting from the degradation of pectic polysaccharides in the cell walls of the grape berry. The increase in galacturonic acid, well above the sensory detection threshold (125 mg/L [35]), could contribute to the wine sourness. Surprisingly, the levels of xylose decreased in both enzymatic treatments. A similar effect was observed in Merlot red wines after enzyme treatments [36]. This could be
due to a competitive effect between the used enzymes and endogenous pectinases that are naturally present in the grape solids. The specific significant increase in catechin confirms the ability of enzymes to enhance phenolic extraction during winemaking [36]. Succinic acid is the main nonvolatile acid produced by yeasts during alcoholic fermentation [37]. Even though this compound could be formed by yeasts until the stationary phase of alcoholic fermentation, the decrease in succinic acid observed with both enzymes cannot be explained by the modulation of yeast activity. An interaction with the products released by the enzymes could be at the origin of this decrease. Interestingly, specific trends of enzymatic preparation E2 were observed (Figure 3c). The use of preparation E2 induced a higher level of arabinose than that from E1. The release of significant quantities of arabinose could be due to an extensive and specific degradation of arabinose-rich polysaccharides by E2. On the contrary, E2 significantly reduced the production of acetoin, isobutanol, acetic acid and leucine in comparison to the E1 and control modalities. As previously mentioned, acetoin is the precursor of 2,3-butanediol, and isobutanol is a higher alcohol involved in wine aroma perception [38]. These results demonstrate the ability of ¹H-NMR metabolomics to monitor the effects of enzyme treatment in winemaking processes.

## 2.4. Fining Treatments

Fining treatments are mainly used in winemaking to clarify and stabilize the wine [39]. They bind to specific compounds (proteins, polymeric phenols or tannins) and form insoluble complexes that can be removed from the wine. Depending on the desired objective, different classes of fining agents can be used. Due to their capacity to bind with tannins, fining agents having a proteinaceous origin are widely used. In this study, the effects of two different fining agents on the wine metabolome were studied by using ¹H-NMR metabolomics. Two fining agents of proteinaceous origin were investigated: a vegetable protein (F1), and a mix of pea protein and PVPP (F2).

Wine samples were treated by finings F1 and F2. The data extracted from ¹H-NMR spectra were compared by PCA (Figure S3b) followed by OPLS-DA to confirm the observed trends. OPLS-DA score and loading plots are shown in Figure 4a,b, respectively. A clear discrimination between the three sets of wine was observed in the OPLS-DA score plot (Q², R²X and R²Y of 0.666, 0.897 and 0.455, respectively). To validate the model, a cross-validation procedure was performed, which indicated that all wines were correctly classified (Table S5). As a first result, these data indicate that ¹H-NMR metabolomics is able to monitor the impact of proteinaceous finings on the wine metabolome. These treatments can induce precipitation of phenolic compounds, particularly polymerized and galloyled tannins [40], affecting the wine organoleptic properties. As shown in Figure 4C, the wines treated by fining F1 presented a higher concentration of amino acids than the untreated wines. These wines contained significantly higher concentrations of alanine and threonine (51 ± 4 and 45 ± 5 mg/L, respectively) in addition to showing an increase in choline (mean 33 ± 3 mg/L) in comparison to the untreated wines and fining F2 treatment. It is surprising to observe an increase in some constituents under the action of the finings. These agents are rather known to precipitate wine constituents. Nevertheless, Maury et al. observed a significant increase in amino acids due to proteinaceous fining treatment in a model wine solution [41]. In addition, a specific action of each fining could be observed. While fining F1 had no significant impact on the main constituents observed by ¹H-NMR, fining F2 induced a significant decrease in some wine constituents such as trigonelline, phenethyl alcohol and isoamyl alcohol (Figure 4c). Phenethyl alcohol is the most important phenolic higher alcohol present in wine [42]. A recent study showed that certain processes can reduce the content of phenyl alcohol in addition to other volatile compounds [43]. Isoamyl alcohol, together with phenethyl alcohol, isobutanol and methionol, is considered one of the most aroma-powerful higher alcohols [38]. In addition, isoamyl alcohol seems to reduce fruity and woody notes in wines. Even though higher alcohol effects depend on the wine aromatic context, ¹H-NMR metabolomics approaches could provide valuable data concerning the effects of fining treatments on wine organic properties.
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Multivariate analysis of ¹H-NMR spectra of wine samples treated by different finings (F0: untreated; F1: fining 1; F2: fining 2): (a) OPLS-DA score plot; (b) loading plot; (c) boxplots of 6 most discriminant wine constituents (t[1] and t[2]: first and second predictive components; pq[1] and pq[2]: first and second predictive component loadings). The significance in the difference was calculated by ANOVA followed by Tukey's multiple comparison test (indicated as \* *p* < 0.05, \*\* *p* < 0.01, \*\*\* *p* < 0.001).</figcaption></figure>

# 3. Materials and Methods

## 3.1. *Wine Samples*

For this study, Cabernet Sauvignon grapes were harvested in 2019 from a vineyard of the Bordeaux Superieur wine region. The grapes were harvested at three different stages of maturity: under-maturity (M1), maturity (M2) and over-maturity (M3), one week between each harvest was observed. Maturity was monitored using a set of different parameters including pH, sugar content, total acidity and malic acid level. These parameters were measured following OIV methods (OIV-OENO 598-2018, OIV-OENO 599-2018, OIV-OENO 600-2018, OIV-MA-AS311-02, OIV-MA-F1-06, OIV-MA-F1-09, OIV-MA-AS313-01) using an Analyzer Y15 (BioSystems, Barcelona, Spain). Data are reported in Table S1.

## 3.2. *Winemaking*

Wines were elaborated at the nanovinification platform (Bordeaux Vinif) of the Institute of Vine and Wine Sciences of Bordeaux (ISVV) following a specific nanovinification method [44]. Briefly, grapes from each of the 3 modalities were destemmed and crushed. The crushed grapes were macerated in 10 L stainless-steel tanks. At the same time, bisulfite (5 mg/L of must) was added. Fermentations were initiated with Actiflore F33 (commercial dry Saccharomyces cerevisiae yeast, Laffort, France) inoculated at 0.1 g/L. Each lot was fermented to completion. The process was monitored by density measurement using the electronic densimeter DMA 35 Basic (Anton Paar France, Ulis, France). At the end of
alcoholic fermentation, each lot was pressed at 1 bar and transferred to 5 L glass tanks. Mal-
olactic fermentations were performed after direct inoculation with commercial Oenoccocus
oeni bacteria SB3 Direct (Laffort, Bordeaux, France). Three weeks after the beginning of
winemaking, the wines were racked, and sulfur dioxide was adjusted to 50 mg/L of wine.
The wines were cold stabilized at 4 °C for 1 week, bottled and stored in the experimental
wine cellar at 15 °C until treatment and analysis. The effects of two enzymes and two
fining agents were compared to untreated samples. For enzyme treatments, two different
pectolytic enzymes were used at the concentration recommended by the provider: enzyme
1 (E1) at 4 g/hL, and enzyme 2 (E2) at 3 mL/hL. The preparations were added at the end
of malolactic fermentation. Samples were kept at 15 °C for 15 days before analysis. For
fining processes, fining agents were added, after malolactic fermentation, to 50 mL of wine
respecting the proportions commonly used in wineries: fining 1 (F1), vegetable protein
at 3 g/hL; fining 2 (F2), pea protein and PVPP (polyvinylpolypyrrolidone) at 80 g/hL.
Samples were kept at 15 °C for 3 days before analysis. The supernatants were collected for
the preparation of the NMR analyses. At least five biological replicates were analyzed for
each treatment.

3.3. *Sample Preparation*

Wines samples were prepared according to our previous work [9]. Briefly, 420 µL wine was mixed with 120 µL phosphate buffer (1 M, pH 2.6), and 60 µL deuterated water containing 0.5 mM sodium trimethylsilylpropanoic acid (TMSP) and 7 mM calcium formate. The pH was adjusted to 3.1 using a small-scale semi-automatic system (BTpH, Bruker BioSpin, Rheinstetten, Germany) using a 1 M HCl solution.

3.4. NMR Analysis

The ¹H-NMR spectra were recorded on a 600 Mhz Avance III NMR spectrometer (Bruker, Wissembourg, France) operating at 600.27 MHz, equipped with a TXI 5 mm probe with z gradient coils. The measurements were performed at 293 K, using Topspin 4.0.8 software (Bruker, Wissembourg, France). Three different sequences were used: the ZG30 and ZGPR sequences were set at 8 ns (number of scans) and NOESYGPPS1D was set at 32 ns for water and ethanol suppressions.

The acquisition parameters were set as follows: The free induction decay (FID) was
collected into a time domain (TD) of 64K data points, with a spectral width (SW) of 16 ppm,
an acquisition time (AQ) of 3.40 s and a relaxation delay (RD) of 5 s per scan. Then,
90° pulse calibration was carried out for each sample automatically, and the shimming was
set manually in gs mode for each spectrum in order to obtain the finest possible line width
(lower than 1 Hz). Water and ethanol signal suppressions were achieved during RD using
a shaped pulse with a multiple-band selective solvent suppression (7 Hz centered on each
ethanol signal, and 20 Hz for water), with a power level for presaturation of 50.37 dB and a
shaped pulse for presaturation of 34.83 dB.

The FID was multiplied by an exponential function corresponding to a 0.3 Hz line-
broadening factor prior to the Fourier transform. Manual phase and baseline correction was
applied to the resulting spectrum, which was then manually phased and aligned to zero
using the TMSP signal. Wine constituents were identified based on previous studies [8,17],
2D NMR analyses and pure chemical standards. Compounds were quantified by targeted
analysis, by the global spectral deconvolution method (GSD) [45], using the simple mixture
analysis (SMA) plugin of MestReNova 12.0 software (Mestrelab Research, Santiago de
Compostela, Spain). Calcium formate was used as an internal standard for quantitation.

3.5. *Multivariate Data Analysis*

Statistical analyses were carried out using SIMCA 16.0 software (Sartorius, Goettingen, Germany). No data preprocessing was performed prior to statistical analysis. Principal component analysis (PCA) was performed to visualize and observe the distribution of the metabolic variance of the dataset. Orthogonal partial least square discriminant analysis
(OPLS-DA) was performed to sharpen the separation between observation groups [27,46]. OPLS-DA was developed as an improvement of the PLS-DA procedure to discriminate classes [47,48]. The results were visualized in the form of score plots and loading plots. Hotelling's T² ellipse, in score plots, defines the 95% confidence region. The model quality was described by R-squared (R²) and Q-squared (Q²) values. R²X and R²Y represent the goodness of fit or explained variation, and Q² represents the goodness of prediction or predicted variation [49]. Results are presented in the form of score plots (each point shows an individual wine sample) and loading plots (each coordinate shows one wine constituent). The variable importance of projection (VIP) represents the total importance of each variable in explaining the model. The p Nature of the data (p-value) indicates the significance of the correlation coefficient between the model and original data (ranged between -1.0 and 1.0) for each loading. OPLS-DA models were validated by analysis of variance testing of cross-validated predictive residuals (CV-ANOVA), with the p-value indicating the model significance. The recognition abilities of the OPLS-DA models were assessed using misclassification tables. The potential discriminant metabolites were identified using VIP and p Nature of the data (p-value) larger than 1.0 and 0.5, respectively, were considered potential discriminant metabolites. One-way analysis of variance (ANOVA) followed by Tukey's multiple comparison procedure was used to identify the discriminant metabolites.

**4. Conclusions**

In our previous studies, we showed that NMR is a valuable tool for ensuring wine traceability complementary to the OIV's official methods [9,50]. NMR combines several advantages that make it a tool of choice for wine metabolomic analysis. In this study, our main goal was to elucidate whether NMR metabolomics could provide valuable data to monitor the effects of winemaking processes including grape maturation and enzyme and fining treatments. The results indicate that NMR metabolomics can discriminate these different winemaking processes. This technique is efficient in quantifying different wine constituents and has the advantages of being fast and reproducible. Combined with multivariate statistical analyses, it allows the identification of some impacted metabolites by each specific winemaking practice. As a complement to other analytical techniques and sensory evaluations, NMR metabolomics could provide valuable new data concerning the impact of winemaking practices on wine organoleptic properties. Additionally, more research is needed to better understand the potential contribution of NMR metabolomics to the investigation of the effects of winemaking practices on wine quality.

**Supplementary Materials:** The following are available online, Figure S1: Boxplots of 16 most discriminant wine constituents based on ¹H-NMR analysis of wines produced from grapes harvested at under-maturity (M1) and at maturity (M2), Figure S2: Boxplots of 10 most discriminant wine constituents based on ¹H-NMR analysis of wines produced from grapes harvested at maturity (M2) and at over-maturity (M3), Figure S3: PCA score plots of ¹H-NMR spectra of wine samples, Table S1: Grape pH, sugar content, total acidity and malic acid contents at the three different harvest times (M1: under-maturity; M2: maturity; M3: over-maturity), Table S2: Cross-validation tables of OPLS analyses of wine samples from grapes at different stages of maturity (M1: under-maturity; M2: maturity; M3: over-maturity), Table S3: OPLS-DA model parameters for discrimination between M1 (under-maturity), M2 (maturity) and M3 (over-maturity). OPLS-DA models were controlled using cross-validation for discrimination between the three groups, and 100 permutation tests of the corresponding validation plots, Table S4: Cross-validation tables of OPLS analyses of wine samples treated by enzymes (E0: control; E1: enzyme 1; E2: enzyme 2), Table S5: Cross-validation tables of OPLS analyses of wine samples treated by fining agents (F0: control; F1: fining 1; F2: fining 2).

**Author Contributions:** Conceptualization, methodology, method validation, formal analysis, investigation, data curation, I.L.M., J.M.-P., C.B., S.L., T.R. and G.D.C.; writing—original draft preparation, I.L.M. and T.R.; writing—review and editing, I.L.M., J.M.-P., S.L., T.R. and G.D.C.; supervision, S.L., T.R. and G.D.C.; funding acquisition, S.L. and T.R. All authors have read and agreed to the published version of the manuscript.
**Funding:** This work was supported by Association Nationale de la Recherche et de la Technologie (Inès Le Mao was the recipient of a CIFRE PhD fellowship from ANRT and Baron Philippe de Rothschild S.A.), Bordeaux Metabolome Facility, MetaboHUB (ANR-11-INBS-0010 project), and the nanovinification platform Bordeaux Vinif. We also want to thank Biolaffort for the financial support during the study.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Data are contained within the article, and they are available from the first corresponding author.

**Conflicts of Interest:** The authors declare no conflict of interest.

**Sample Availability:** Samples are available from the first corresponding author.

## References

1. Cevallos-Cevallos, J.M.; Reyes-De-Corcuera, J.I.; Etxeberria, E.; Danyluk, M.D.; Rodrick, G.E. Metabolomic analysis in food science: A review. *Trends Food Sci. Technol.* **2009**, *20*, 557-566. [CrossRef]

2. Kim, S.; Kim, J.; Yun, E.J.; Kim, K.H. Food metabolomics: From farm to human. *Curr. Opin. Biotechnol.* **2016**, *37*, 16-23. [CrossRef]

3. Valls Fonayet, J.; Loupit, G.; Richard, T. Chapter Ten—MS- and NMR-metabolomic tools for the discrimination of wines: Applications for authenticity. *Adv. Bot. Res.* **2021**, *98*, 297-357.

4. Viskić, M.; Bandić, L.M.; Korenika, A.-M.J.; Jeromel, A. NMR in the service of wine differentiation. *Foods* **2021**, *10*, 120. [CrossRef] [PubMed]

5. Solovyev, P.A.; Fauhl-Hassek, C.; Riedl, J.; Esslinger, S.; Bontempo, L.; Camin, F. NMR spectroscopy in wine authentication: An official control perspective. *Compr. Rev. Food Sci. Food Saf.* **2021**, *20*, 2040-2062. [CrossRef] [PubMed]

6. Sobolev, A.P.; Thomas, F.; Donarski, J.; Ingallina, C.; Circi, S.; Cesare Marincola, F.; Capitani, D.; Mannina, L. Use of NMR applications to tackle future food fraud issues. *Trends Food Sci. Technol.* **2019**, *91*, 347-353. [CrossRef]

7. Tabago, M.K.A.G.; Calingacion, M.N.; Garcia, J. Recent advances in NMR-based metabolomics of alcoholic beverages. *Food Chem. Mol. Sci.* **2021**, *2*, 100009. [CrossRef]

8. Gougeon, L.; Da Costa, G.; Le Mao, I.; Ma, W.; Teissedre, P.-L.; Guyon, F.; Richard, T. Wine analysis and authenticity using 1H-NMR metabolomics data: Application to Chinese wines. *Food Anal. Methods* **2018**, *11*, 3425-3434. [CrossRef]

9. Gougeon, L.; da Costa, G.; Guyon, F.; Richard, T. ¹H NMR metabolomics applied to Bordeaux red wines. *Food Chem.* **2019**, *301*, 125257. [CrossRef] [PubMed]

10. Godelmann, R.; Kost, C.; Patz, C.-D.; Ristow, R.; Wachter, H. Quantitation of compounds in wine using ¹H NMR spectroscopy: Description of the method and collaborative study. *J. AOAC Int.* **2016**, *99*, 1295-1304. [CrossRef]

11. Son, H.-S.; Hwang, G.-S.; Kim, K.M.; Ahn, H.-J.; Park, W.-M.; Van Den Berg, F.; Hong, Y.-S.; Lee, C.-H. Metabolomic studies on geographical grapes and their wines using ¹H NMR analysis coupled with multivariate statistics. *J. Agric. Food. Chem.* **2009**, *57*, 1481-1490. [CrossRef] [PubMed]

12. Fan, S.; Zhong, Q.; Fauhl-Hassek, C.; Pfister, M.K.H.; Horn, B.; Huang, Z. Classification of Chinese wine varieties using ¹H NMR spectroscopy combined with multivariate statistical analysis. *Food Control* **2018**, *88*, 113-122. [CrossRef]

13. Caruso, M.; Galgano, F.; Castiglione Morelli, M.A.; Viggiani, L.; Lencioni, L.; Giussani, B.; Favati, F. Chemical profile of white wines produced from 'Greco bianco' grape variety in different Italian areas by nuclear magnetic resonance (NMR) and conventional physicochemical analyses. *J. Agric. Food. Chem.* **2012**, *60*, 7-15. [CrossRef] [PubMed]

14. Pereira, G.E.; Gaudillère, J.P.; Van Leeuwen, C.; Hilbert, G.; Maucourt, M.; Deborde, C.; Moing, A.; Rolin, D. ¹H-NMR metabolic profiling of wines from three cultivars, three soil types and two contrasting vintages. *J. Int. Sci. Vigne Vin* **2007**, *41*, 103-109.

15. Lee, J.-E.; Hwang, G.-S.; Van Den Berg, F.; Lee, C.-H.; Hong, Y.-S. Evidence of vintage effects on grape wines using ¹H NMR-based metabolomic study. *Anal. Chim. Acta* **2009**, *648*, 71-76. [CrossRef]

16. Cassino, C.; Tsolakis, C.; Bonello, F.; Gianotti, V.; Osella, D. Wine evolution during bottle aging, studied by ¹H NMR spectroscopy and multivariate statistical analysis. *Food Res. Int.* **2018**, *116*, 566-577. [CrossRef]

17. Amargianitaki, M.; Spyros, A. NMR-based metabolomics in wine quality control and authentication. *Chem. Biol. Technol. Agric.* **2017**, *4*, 9. [CrossRef]

18. Alves Filho, E.G.; Silva, L.M.A.; Ribeiro, P.R.V.; de Brito, E.S.; Zocolo, G.J.; Souza-Leão, P.C.; Marques, A.T.B.; Quintela, A.L.; Larsen, F.H.; Canuto, K.M. ¹H NMR and LC-MS-based metabolomic approach for evaluation of the seasonality and viticultural practices in wines from São Francisco River Valley, a Brazilian semi-arid region. *Food Chem.* **2019**, *289*, 558-567. [CrossRef]

19. López-Rituerto, E.; Savorani, F.; Avenoza, A.; Busto, J.H.; Peregrina, J.M.; Engelsen, S.B. Investigations of la Rioja terroir for wine production using ¹H NMR metabolomics. *J. Agric. Food. Chem.* **2012**, *60*, 3452-3461. [CrossRef] [PubMed]

20. Mazzei, P.; Spaccini, R.; Francesca, N.; Moschetti, G.; Piccolo, A. Metabolomic by ¹H NMR spectroscopy differentiates "Fiano Di Avellino" white wines obtained with different yeast strains. *J. Agric. Food. Chem.* **2013**, *61*, 10816-10822. [CrossRef]
21. De Pascali, S.A.; Coletta, A.; Del Coco, L.; Basile, T.; Gambacorta, G.; Fanizzi, F.P. Viticultural practice and winemaking effects on metabolic profile of Negroamaro. *Food Chem.* **2014**, *161*, 112–119. [CrossRef] [PubMed]

22. Baiano, A.; Terracone, C.; Longobardi, F.; Ventrella, A.; Agostiano, A.; Del Nobile, M.A. Effects of different vinification technologies on physical and chemical characteristics of Sauvignon blanc wines. *Food Chem.* **2012**, *135*, 2694–2701. [CrossRef]

23. González-Neves, G.; Favre, G.; Gil, G. Effect of fining on the colour and pigment composition of young red wines. *Food Chem.* **2014**, *157*, 385–392. [CrossRef]

24. Doco, T.; Williams, P.; Cheynier, V. Effect of flash Release and pectinolytic enzyme treatments on wine polysaccharide composition. *J. Agric. Food. Chem.* **2007**, *55*, 6643–6649. [CrossRef]

25. Parish, K.J.; Herbst-Johnstone, M.; Bouda, F.; Klaere, S.; Fedrizzi, B. Pre-fermentation fining effects on the aroma chemistry of Marlborough Sauvignon blanc press fractions. *Food Chem.* **2016**, *208*, 326–335. [CrossRef] [PubMed]

26. Arnous, A.; Meyer, A.S. Discriminated release of phenolic substances from red wine grape skins (*Vitis vinifera* L.) by multicomponent enzymes treatment. *Biochem. Eng. J.* **2010**, *49*, 68–77. [CrossRef]

27. Wheelock, Å.M.; Wheelock, C.E. Trials and tribulations of 'omics data analysis: Assessing quality of SIMCA-based multivariate models using examples from pulmonary medicine. *Mol. BioSyst.* **2013**, *9*, 2589–2596. [CrossRef]

28. Benbouguerra, N.; Valls-Fonayet, J.; Krisa, S.; Garcia, F.; Saucier, C.; Richard, T.; Hornedo-Ortega, R. Polyphenolic characterization of Merlot, Tannat and Syrah skin extracts at different degrees of maturity and anti-inflammatory potential in RAW 264.7 cells. *Foods* **2021**, *10*, 541. [CrossRef]

29. Ferrero-del-Teso, S.; Arias, I.; Escudero, A.; Ferreira, V.; Fernández-Zurbano, P.; Sáenz-Navajas, M.-P. Effect of grape maturity on wine sensory and chemical features: The case of Moristel wines. *LWT* **2020**, *118*, 108848. [CrossRef]

30. Bindon, K.; Varela, C.; Kennedy, J.; Holt, H.; Herderich, M. Relationships between harvest time and wine composition in Vitis vinifera L. cv. Cabernet Sauvignon 1. Grape and wine chemistry. *Food Chem.* **2013**, *138*, 1696–1705. [CrossRef]

31. Lee, J.E.; Hong, Y.S.; Lee, C.H. Characterization of fermentative behaviors of lactic acid bacteria in grape wines through ¹H NMR- and GC-based metabolic profiling. *J. Agric. Food. Chem.* **2009**, *57*, 4810–4817. [CrossRef]

32. Long, D.; Wilkinson, K.L.; Poole, K.; Taylor, D.K.; Warren, T.; Astorga, A.M.; Jiranek, V. Rapid method for proline determination in grape juice and wine. *J. Agric. Food. Chem.* **2012**, *60*, 4259–4264. [CrossRef] [PubMed]

33. Muñoz, O.; Sepúlveda, M.; Schwartz, M. Effects of enzymatic treatment on anthocyanic pigments from grapes skin from chilean wine. *Food Chem.* **2004**, *87*, 487–490. [CrossRef]

34. Li, S.; Bindon, K.; Bastian, S.E.P.; Jiranek, V.; Wilkinson, K.L. Use of winemaking supplements to modify the composition and sensory properties of Shiraz wine. *J. Agric. Food. Chem.* **2017**, *65*, 1353–1364. [CrossRef] [PubMed]

35. Hufnagel, J.C.; Hofmann, T. Quantitative reconstruction of the nonvolatile sensometabolome of a red wine. *J. Agric. Food. Chem.* **2008**, *56*, 9190–9199. [CrossRef]

36. Ducasse, M.-A.; Canal-Llauberes, R.-M.; de Lumley, M.; Williams, P.; Souquet, J.-M.; Fulcrand, H.; Doco, T.; Cheynier, V. Effect of macerating enzyme treatment on the polyphenol and polysaccharide composition of red wines. *Food Chem.* **2010**, *118*, 369–376. [CrossRef]

37. Thoukis, G.; Ueda, M.; Wright, D. The formation of succinic acid during alcoholic fermentation. *Am. J. Enol. Vitic.* **1965**, *16*, 1–8.

38. De-la-Fuente-Blanco, A.; Sáenz-Navajas, M.-P.; Ferreira, V. On the effects of higher alcohols on red wine aroma. *Food Chem.* **2016**, *210*, 107–114. [CrossRef] [PubMed]

39. Marangon, M.; Vincenzi, S.; Curioni, A. Wine Fining with Plant Proteins. *Molecules* **2019**, *24*, 2186. [CrossRef] [PubMed]

40. Gil, M.; Louazil, P.; Iturmendi, N.; Moine, V.; Cheynier, V.; Saucier, C. Effect of polyvinylpolypyrrolidone treatment on rosés wines during fermentation: Impact on color, polyphenols and thiol aromas. *Food Chem.* **2019**, *295*, 493–498. [CrossRef]

41. Maury, C.; Sarni-Manchado, P.; Cheynier, V. Highlighting protein fining residues in a model red wine. *Food Chem.* **2019**, *279*, 272–278. [CrossRef]

42. Jackson, R.S. Chapter 6—Chemical constituents of grapes and wine. In *Wine Science*, 5th ed.; Jackson, R.S., Ed.; Academic Press: Cambridge, MA, USA, 2020; pp. 375–459.

43. Ma, T.-Z.; Gong, P.-F.; Lu, R.-R.; Zhang, B.; Morata, A.; Han, S.-Y. Effect of different clarification treatments on the volatile composition and aromatic attributes of 'Italian Riesling' icewine. *Molecules* **2020**, *25*, 2657. [CrossRef]

44. Rousserie, P.; Lacampagne, S.; Vanbrabant, S.; Rabot, A.; Geny-Denis, L. Wine tannins: Where are they coming from? A method to access the importance of berry part on wine tannins content. *MethodsX* **2020**, *7*, 100961. [CrossRef]

45. Cobas, C.; Seoane, F.; Domínguez, S.; Sykora, S.; Davies, A.N. A new approach to improving automated analysis of proton NMR spectra through Global Spectral Deconvolution (GSD). *Spectrosc. Eur.* **2011**, *23*, 26–30.

46. Tyagi, R.; Maan, K.; Khushu, S.; Rana, P. Urine metabolomics based prediction model approach for radiation exposure. *Sci. Rep.* **2020**, *10*, 16063. [CrossRef] [PubMed]

47. Westerhuis, J.A.; van Velzen, E.J.J.; Hoefsloot, H.C.J.; Smilde, A.K. Multivariate paired data analysis: Multilevel PLSDA versus OPLSDA. *Metabolomics* **2010**, *6*, 119–128. [CrossRef]

48. Trygg, J.; Wold, S. Orthogonal projections to latent structures (O-PLS). *J. Chemom.* **2002**, *16*, 119–128. [CrossRef]

49. Ruiying, C.; Zeyun, L.; Yongliang, Y.; Zijia, Z.; Ji, Z.; Xin, T.; Xiaojian, Z. A comprehensive analysis of metabolomics and transcriptomics in non-small cell lung cancer. *PLoS ONE* **2020**, *15*, e0232272. [CrossRef]

50. Gougeon, L.; da Costa, G.; Richard, T.; Guyon, F. Wine authenticity by quantitative 1 H NMR versus multitechnique analysis: A case study. *Food Anal. Methods* **2019**, *12*, 956–965. [CrossRef]