Next Article in Journal
Trimethylsilylethynyl-Substituted Pyrene Doped Materials as Improved Fluorescent Sensors towards Nitroaromatic Explosives and Related Compounds
Next Article in Special Issue
Hierarchical Modeling to Enhance Spectrophotometry Measurements—Overcoming Dynamic Range Limitations for Remote Monitoring of Neptunium
Previous Article in Journal
Near-Infrared Molecular Photosensitizer Decorated with Quaternary Ammonium for High-Efficiency Photothermal Treatment of Bacterial Infections
Previous Article in Special Issue
Direct Determination of Three PAHs in Drill Cuttings Recycling Products by Solid-Surface 3D Fluorescence Coupled with Chemometrics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Strategy for Differentiating Oak Barrel Aged and Non-Oak Barrel Aged Wines by Using UHPLC–HRMS Combined with Chemometrics

1
Key Laboratory of Quality and Safety of Wolfberry and Wine for State Administration for Market Regulation, Ningxia Food Testing and Research Institute, Yinchuan 750004, China
2
College of Pharmacy, Ningxia Medical University, Yinchuan 750004, China
3
Key Laboratory of Ningxia Minority Medicine Modernization, Ministry of Education, Yinchuan 750004, China
4
Institute of Quality Standard and Testing Technology for Agro-Products of Ningxia, Yinchuan 750002, China
*
Author to whom correspondence should be addressed.
Chemosensors 2023, 11(3), 165; https://doi.org/10.3390/chemosensors11030165
Submission received: 3 February 2023 / Revised: 25 February 2023 / Accepted: 27 February 2023 / Published: 1 March 2023
(This article belongs to the Special Issue Chemometrics for Analytical Chemistry)

Abstract

:
The time involved and the high economic cost of using oak barrels to age wines have led to a significant price difference compared to non-oak barrel aged wines and may lead to some fraudulent sales in the market. In this study, an untargeted metabolomic strategy was developed to detect the metabolite composition of oak barrel aged and non-oak barrel aged wines in both positive and negative ion modes by using UHPLC–HRMS combined with the recently developed chemometric method AntDAS. The results of partial least squares discrimination analysis (PLS-DA) showed that the samples were characterized into two categories. Finally, 51 compounds were identified in positive ion mode, while 26 compounds were identified in negative ion mode. The results indicate that combining UHPLC–HRMS with AntDAS can reveal the material basis of wines and has excellent potential to differentiate between oak barrel aged and non-oak barrel aged wines.

1. Introduction

Wine is a popular alcoholic beverage that has been consumed for hundreds of years [1]. Wine constitutes a reliable and rich source of biologically active phytochemicals, specifically phenolic acids and polyphenols, whose individual and summated actions are believed to provide health benefits [2]. In addition, wine is widely consumed by the public for its pleasant organoleptic properties. According to the literature [3], the organoleptic characteristics of wine, such as smell, taste, and color, largely depend on the wine-making and aging techniques used. The aging of wine is a process widely used in oenology to provide added value to the wine [4]. During the production process, wine-making methods, such as maceration, fermentation, clarification, and aging, can significantly alter the composition and concentration of compounds [5] and substantially enhance the complexity of wine by combining these compounds with important aging characteristic volatiles [6]; this, in turn, affects the wine’s composition and sensory characteristics. During ageing processes wine penetrates oak wood, Maximum sorption capacity depends principally on wood type [7]. Because of the time involved (6 to 18 months in barrels) and the high economic cost [8], the price difference between commercially available oak barrel aged wines and non-oak barrel aged wines is significant, leading to several consumer detriments and fraudulent sales. Therefore, there is a great need to establish a scientific and effective method to distinguish between oak barrel aged and non-oak barrel aged wines. This would not only provide protection for consumers but also provide solid theoretical support for the Chinese wine industry and government supervision.
Nowadays, the emergence of analytical instruments with high separation efficiency, high resolution, and excellent sensitivity coupled with advanced chemometric techniques provides new research perspectives and approaches for analyzing wines and assessing their authenticity. For example, Arapitsas et al. used an untargeted LC–MS method to identify several metabolite markers, including ascorbic acid and tartaric acid, as well as the first identified sulfonated compounds [9]. Gil et al. used ultra-performance liquid chromatography coupled with quadrupole/time of flight mass spectrometry (UPLC–QTOF-MS) to obtain high-resolution mass spectrometric metabolomic fingerprints of rosé wines and then applied a statistical analysis to investigate the discrimination of rosé wines by origin [10]. Ontañón et al. used non-targeted liquid chromatography–mass spectrometry (LC–MS) to analyze the metabolomic space introduced by wine storage under different conditions. Credible markers that may reveal major chemical changes related to storage conditions were also found [11]. In addition to the above studies, other researchers endeavored to explore the organoleptic properties of aged wines, such as Laureati et al., who used gas chromatography–mass spectrometry (GC–MS), gas chromatography–flame ionization detection (GC–FID), and sensory analysis by a wine tasting panel to evaluate the organoleptic characteristics of wines aged for extended periods (up to 15 years) [12].
As a general overview, the existing studies all describe the application of wine-based metabolomics in the discovery of novel compounds, characterization and identification of wine groups [13], prediction of sensory attributes [14], monitoring of storage and aging processes [15,16], and traceability to their varietal and geographical origin [17,18]. It is important to note that little attention has been paid to the effects of the wine barrel aging process on the specific distribution of compounds and the exploration of characteristic compounds between oak barrel aged and non-oak barrel aged wines. Consequently, it is necessary to establish a scientific and effective method to explore the characteristic compounds to distinguish oak barrel aged and non-oak barrel aged wines.
An untargeted metabolomic analysis based on UHPLC–HRMS combined with chemometrics was chosen for this study to detect the metabolite composition of oak barrel aged and non-oak barrel aged wines. This approach aims to compare the metabolite patterns of wines before and after aging to distinguish oak barrel aged wines from non-oak barrel aged wines. The metabolite features were effectively extracted with the automatic target-free data analysis strategy (AntDAS) [19], developed by Yu Yongjie’s research group, to reveal the differences between the metabolites of the two species. In addition, to better exploit the potential of the data, clustering and discriminant analysis were performed using methods including principal component analysis (PCA) and partial least squares discriminant analysis (PLS-DA). The differences in the chemical basis of oak barrel aged and non-oak barrel aged wines were explored and evaluated.

2. Materials and Methods

2.1. Samples and Chemicals

In this study, 66 wine samples from the wine produced in the eastern foothills of Helan Mountain, Ningxia, China, were collected in 2022, including both oak barrel aged (n = 12) and non-oak barrel aged (n = 54) wines. The samples were stored in a refrigerator at −18 °C and protected from light until further use. In addition, the chemical reagents used in this study included HPLC-grade methanol obtained from J.T. Baker, Avantor (Radnor, PA, USA), and HPLC-grade formic acid from Aladdin (Shanghai, China). The LC–MS-grade ammonium acetate that was used as a mobile phase additive was purchased from Fisher Chemical (Waltham, MA, USA). The distilled water was obtained from Watson’s Food & Beverage (Guangzhou, China).

2.2. Sample Preparation

Prior to analysis, the wine samples were removed from the −18 °C refrigerator and thawed. A total of 1.5 mL of wine was placed in a 2 mL centrifuge tube, which was centrifuged at 13,000 r/min for 10 min followed by transferring 1 mL of the supernatant to a chromatography vial for subsequent UHPLC–HRMS analysis. Furthermore, quality control (QC) samples were prepared by mixing all wine samples to be measured in equal volumes under the same conditions to investigate the stability of the method.

2.3. UHPLC–QTOF-MS Analysis

For the untargeted analysis of the wine, an UltiMate 3000 ultra-high-performance liquid chromatograph coupled with a Q Exactive-Orbitrap Mass Spectrometer (Thermo Fisher Scientific, Waltham, MA, USA) was used. The chromatographic separation was performed using an Agilent ZORBAX SB-C18 column (2.1 × 100 mm, 1.8 μm) (Agilent, Santa Clara, CA, USA) with the column temperature set at 40 °C. The mobile phase consisted of 0.1% formic acid in water (solvent A) and 0.1% formic acid in acetonitrile (solvent B). The flow rate was set to 0.2 mL/min, and the gradient elution procedure was as follows: 0–1.00 min, 5–8% B; 1.00–7.00 min, 8–12% B; 7.00–13.00 min, 12–20% B; 13.00–20.00 min, 20–35% B; 20.00–24.00 min, 35–55% B; 24.00–25.00 min, 55–97% B; 25.00–30.00 min, 97–100% B; 30.00–35.00 min, 100–100% B; 35.00–35.50 min, 100–5% B. Sequential sampling analysis was performed with an injection volume of 2μL, and QC was added once for every 10 samples.
Mass spectrometry was performed using an electrospray ionization source (ESI) in positive and negative ion modes. Full scan and data-dependent MS/MS acquisition modes (FULLMS/ddMS2) were selected for the data collection. The optimal source parameters were optimized as follows: For the parameter setting of MS1, the capillary temperature was set to 320 °C; the spray voltage was 3.5 kV; the resolution was 35,000; the AGC target was 5 × 106; and the maximum IT was set to 100 ms for scanning in the range of 120~1000 Da. The auxiliary gas and sheath gas flow rates were 10 and 30 arcs, respectively. The optimized parameters of dd-MS2 were as follows: collision energy 20 eV; resolution 17,500; and AGC target 5× 106, where the minimum AGC target was 5 × 10³ with an intensity threshold 2 × 105 and dynamic exclusion 4.5 s. Unlike MS1, the maximum IT of dd-MS2 was set to 25 ms. The isolation window was set to 0.4 Da.

2.4. Data Analysis

Data processing and analysis were performed with AntDAS, a data processing software. It can directly analyze the raw data collected by the instrument, including extract ion chromatograpy (EIC) automatic peak construction, EIC peak extraction, peak labeling, sample time shift correction, etc. [20]. After data processing, a list of compound information with sample × peak area was obtained. Analysis of variance (ANOVA) was used to screen for differential components and their MS/MS spectra. Principal component analysis (PCA), partial least squares discriminant analysis (PLS-DA), and heat map analysis were performed on the screened differential components to evaluate the differences in chemical basis between the oak barrel aged and non-oak barrel aged wines. Additionally, AntDAS completed the metabolite identification by matching the compound mass spectra to third-party databases (http://prime.psc.riken.jp/compms/msdial/main.html (accessed on 25 February 2023).

3. Results and Discussion

3.1. UHPLC–HRMS Untargeted Metabolomic Analysis of Oak Barrel Aged and Non-Oak Barrel Aged Wines

The total ion chromatograms (TICs) collected in the positive and negative ionization modes are shown in Figure 1. The TIC peaks in both modes are basically well-separated throughout the outflow section. It is also noteworthy that there is a significant difference between the two TICs. Overall, the number of peaks collected in the positive ion mode (39,104) is more abundant than in the negative ion mode (31,944). However, extra consideration needs to be given to the fact that some metabolites may lead to ionization in only one ionization mode due to their unique chemical properties. If only one ionization mode is used for data acquisition, a number of key metabolites may be neglected. Given these considerations, both positive and negative ionization modes were used for the data analysis and screening of the metabolic markers in this study.
To obtain the data matrix, including retention time (RT), m/z, and peak intensity for each sample, the collected raw data were imported into the AntDAS [21] data automation parsing software for data analysis. A list of compound information of 39,104 × 74 was obtained in the positive ionization mode, where 39,104 is the number of EIC peaks and 74 is the number of samples (containing sixty-six samples and eight QCs). For the same samples collected in the negative ionization mode, a data matrix of 31,944 × 74 was obtained. In order to demonstrate the potential of AntDAS peak extraction, [M+H]+ and [M−H] peaks were extracted from the compounds with an accurate molecular weight of 304.0583. The results are shown in Figure 2a, b. Figure 2a shows the extraction of the m/z 305.0655 peaks under positive ions from which it is clear that 21 EIC peaks were well extracted. Another peak extraction at m/z 303.0560 under negative ions is shown in Figure 2b, where 15 EIC peaks were extracted. On the whole, most of the EIC peaks basically achieved baseline separation. As shown in the local magnification of the peak extraction results in Figure 2a, b, AntDAS demonstrates the ability to achieve a more reasonable segmentation of the overlapping signals of the compounds in both positive and negative ion modes. On this basis, it can accurately identify the EIC peaks of the compounds from the background noise of the instrument and accurately estimate the outflow range of the compounds.
To examine the reproducibility and stability of the instrument, the data quality of the QC samples was evaluated. The data quality was assessed by calculating the relative standard deviation (RSD) of each EIC peak area of the QC samples with the positive and negative ion pattern results shown in Figure 2c,d. From the figure, it can be determined that the number of peaks with a RSD ≤ 30% in the positive ion mode is 26,569, accounting for 67.9% of the total number of peaks, compared to 21,332 peaks with a RSD ≤ 30% in the negative ion mode, accounting for 66.8% of the total number of peaks. The number of peaks with a RSD ≤ 20% in the positive ion mode was 18,882, accounting for 48.3% of the total number of peaks, while the number of peaks in this range in the negative ion mode was 16,425, accounting for 51.4% of the total number of peaks. At the level of RSD ≤ 10%, the number of peaks in the negative ion mode still exceeded 32% of the total number of peaks. The data given in the above figures indicate that the method is stationary during this batch injection process, and the quality of the acquired data is reliable for subsequent metabolomics research.

3.2. Chemometric Analysis of High-Throughput Metabolomic Data

To explore the large amount of information implied in the already established data matrix, this study used a rational mathematical and statistical approach to transform the chemical model into a mathematical model in order to permit the maximum collection and extraction of useful information [22]. In this study, the EIC peaks of the registered compounds were screened and analyzed using ANOVA at a p-value of 0.05. Approximately 774 EIC peaks that differed between the groups were screened in the positive ion mode, while 766 differential compounds were screened in the negative ion mode. Based on the differential compounds screened above, the principal component analysis (PCA), which is the unsupervised pattern recognition method [23], was first used in this study to investigate the chemical differences between the oak barrel aged and non-oak barrel aged wines. The PCA plots for the positive and negative ionization pattern data are shown in Figure 3a, b. As can be observed from the figure, the first two principal components (PCs), which are used to explain the matrix information, account for only 36.3% in the ESI+ mode (Figure 3a) and 36.7% in the ESI− mode (Figure 3b). The samples of oak barrel aged wine are not well distinguished from the samples of non-oak barrel aged wine. The data indicate that, according to the unsupervised PCA analysis, no distinction can be made between the oak barrel aged and non-oak barrel aged wines in both the positive and negative ion modes. The analytical results might be due to the large sample difference within the same group, and the difference between the different groups was difficult to find using the unsupervised analysis method.
In view of this, to solve the problems encountered in the unsupervised analysis and observe the overall clustering effect and classification trend of the data matrix, a widely used supervised pattern recognition method, PLS-DA [24], was adopted. The cluster analysis plots obtained with PLS-DA are shown in Figure 3c, d. Initially, in the positive ionization mode, Figure 3c presents a well-clustered effect and a trend of better classification between the oak barrel aged and non-oak barrel aged wine samples. The clustering effect in the negative ionization mode shown in Figure 3d is quite similar to that in Figure 3c. As can be observed, some of the non-oak barrel aged wine samples overlap with the oak barrel aged wine samples. A preliminary judgment was made that the overlapped samples were commercially available adulterated oak barrel aged wines. Overall, there was a better classification trend between the oak barrel aged and non-oak barrel aged wine samples by using PLS-DA compared to PCA. The Monte Carlo simulation method was used to further verify the applicability of PLS-DA. The simulation was performed approximately 10,000 times and, for each simulation, 80% of the samples in a group were selected for building the PLS-DA model to predict the remaining 20% of the samples. Finally, the average prediction accuracies for the positive and negative ion modes were 86.3 ± 7.0% and 84.6 ± 7.5%, respectively. These results indicate that PLS-DA can be used for the classification of oak barrel aged and non-oak barrel aged wines.

3.3. Identification and Analysis of Candidate Metabolites

The compound identification was performed with AntDAS as well. An advantage of AntDAS is that it can identify the insource fragment ions and give detailed information on the MS1 and MS/MS spectra for compound identification [25]. Taking L-tryptophan as an example, the identification process of the compound based on the AntDAS platform is shown in Figure 4. The MS1 spectrum of the compound constructed with AntDAS is shown in Figure 4a, whose spectrum consists mainly of fragment ions formed by cleavage within the source of the compound at m/z 205.0973. From the MS/MS spectrum in Figure 4a, it is observed that the ions form mainly fragment ions with m/z 118.0657, m/z 146.0603, and m/z 188.0710 after the ion break, which is consistent with the ions in MS1 in Figure 4a. Both primary and secondary matches of this compound are above 0.95. Overall, these results indicate that AntDAS can construct accurate MS1 and MS/MS spectra of compounds from complex data. For data visualization, Figure 4b shows the EIC plots corresponding to each fragment ion with highly consistent shapes of the outflow curves in the plots. Identifying the compound using the MS1 and MS/MS spectra constructed with AntDAS gives the best match for the compound. The final identification result of the compound is presented in Figure 4c. Following the above compound identification process, 51 compounds were identified in the positive ion mode (as shown in Table 1), and 26 compounds were identified in the negative ion mode (as shown in Table 2).
A heat map was drawn for the identified compounds to further explore the variation patterns of the compounds between the two wines. Figure 5 shows the heat map visualization of the metabolites identified in the positive and negative ion modes for 66 wine samples. The heat maps show the compounds in the vertical direction and the wine samples in the horizontal direction. Additionally, the red color represents the compounds with a relative content above the average, and the green color represents the compounds with a relative content below the average. The heat map in the positive ion mode is shown in Figure 5a. The result that can be drawn from the graph is that the heat map analysis does not classify the samples into different categories. In other words, the metabolite profiles identified do not allow a clear distinction between oak barrel aged and non-oak barrel aged wines. This may be due to the limited potential of the compound database during the identification process, which led to differential metabolites not being identified through the compound database.
For an in-depth exploration of how these metabolites behave in the two wine samples, the relative amounts of the compounds in each group of wine samples are shown as violin plots in Figure 6 for a comparative analysis. In this study, the main chemical components in the wines were flavonoids, phenolic acids, etc. Therefore, the discussion of the compound content distribution in this section will focus on these compounds. In Figure 6a, four compounds, including two flavonoids (catechin, epicatechin) and two phenolic compounds (caffeic acid, p-coumaric acid), were demonstrated in the positive ion mode. From the information presented in the figure, it is possible to obtain a higher content of flavonoid compounds in oak barrel aged wines even if they did not show a significant difference between the two groups. In previous studies [26], compounds such as catechin and epicatechin decreased in content after aging compared to the initial wines. Such differences may be caused by the characteristics of the oak barrels as well as other parameters such as grape varieties [27]. Alternatively, the graph shows that the levels of coumaric and caffeic acids are low in non-oak barrel aged wines and high in oak barrel aged wines. Studies suggest [28] that this may be due to an increase in coumaric and/or caffeic acids as a result of the contact between the wine and the wood (which may contain p-coumaric and/or caffeic acids) during the aging process in oak barrels. Other studies have explained that the loss of coumarin-based anthocyanins during aging can also increase p-coumaric acid content [29].
The analysis of the four compounds in the negative ion mode is presented in Figure 6b. As can be observed, the levels of 3-phenylactic acid and tartaric acid were higher in the oak barrel aged wine samples than in the non-oak barrel aged wine samples. It is also interesting to note that the level of quercetin was lower in the oak barrel aged wine samples than in the non-oak barrel aged wine samples. According to the literature research [30], oxidation and condensation reactions involving sapogenins and glycoside derivatives could explain the reduction in myricetin and quercetin. Furthermore, the higher content of syringic acid in the oak barrel aged wines is consistent with the results reported by Cadahía et al. [26]. This is because compounds such as gallic acid and its dimers, ellagic acid and other benzoic acids (syringic and vanillic), may be present in oak, and they can be transferred to the wine during aging.

4. Conclusions

A combination of AntDAS and UHPLC–HRMS was used to identify metabolites in wine samples. A total of 51 compounds were identified in the positive ion mode, including eleven amino acids, nine flavonoids, four phenolic acids, and individual compounds including alkaloids and coumarins, while 26 compounds were identified in the negative ion mode, including two amino acids, eight organic acids, four flavonoids, and other compounds. In addition, the unsupervised pattern recognition method PCA and supervised pattern recognition method PLS-DA were used for the clustering analysis of oak barrel aged and non-oak barrel aged wines, and from the results, it is observed that PLS-DA can achieve good clustering and differentiation performance. The application of UHPLC–HRMS in food metabolomics for identification and other purposes is a promising approach.

Author Contributions

Conceptualization, Y.L. and Y.-J.Y.; methodology, Y.Z.; software, Y.-J.Y.; validation, Y.-J.Y., Y.L. and F.-L.M.; formal analysis, Q.G.; investigation, J.-N.W.; resources, Y.Z.; data curation, Y.J. and Q.G.; writing—original draft preparation, F.-L.M. and Y.L.; writing—review and editing, Y.-J.Y. and Y.Z.; visualization, F.-L.M. and J.-N.W.; supervision, Y.J. and Q.G.; project administration, Y.L. and Y.-J.Y.; funding acquisition, Y.L. and Y.-J.Y. All authors have read and agreed to the published version of the manuscript.

Funding

The authors gratefully acknowledge the financial support of the Foundation of the National Natural Science Foundation of China (Grant Number: 21868028), the Key Research and Development Project of Ningxia (Grant Number: 2020BEG03039, 2022BEG03170), research projects of Ningxia Administration for Market Regulation (Grant Number: 2022SJKY0001, 2022SJKY0002, 2022SJKY0003).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Castaldo, L.; Narváez, A.; Izzo, L.; Graziani, G.; Gaspari, A.; Minno, G.D.; Ritieni, A. Red wine consumption and cardiovascular health. Molecules 2019, 24, 3626. [Google Scholar] [CrossRef] [Green Version]
  2. German, J.B.; Walzem, R.L. The health benefits of wine. Annu. Rev. Nutr. 2000, 20, 561–593. [Google Scholar] [CrossRef] [PubMed]
  3. Li, S.Y.; Duan, C.Q. Astringency, bitterness and color changes in dry red wines before and during oak barrel aging: An updated phenolic perspective review. Crit. Rev. Food Sci. 2019, 59, 1840–1867. [Google Scholar] [CrossRef] [PubMed]
  4. Teresa Garde-Cerdán; Carmen Ancín-Azpilicueta. Review of quality factors on wine ageing in oak barrels. Trends Food Sci. Tech. 2006, 17, 438–447. [Google Scholar] [CrossRef]
  5. Casassa, L.F.; Harbertson, J.F. Extraction, evolution, and sensory impact of phenolic compounds during red wine maceration. Annu. Rev. Food. Sci. Technol. 2014, 5, 83–109. [Google Scholar] [CrossRef]
  6. Xu, Q.; Fang, Y.J.; Jian, C.; Ying, S.; Chang, Q.D.; Yi, B.L. Characterization and evolution of volatile compounds of cabernet sauvignon wines from two different clones during oak barrel aging. Foods 2022, 11, 74. [Google Scholar]
  7. Coelho, E.; Domingues, L.; Teixeira, J.A.; Oliveira, J.M.; Tavares, T. Understanding wine sorption by oak wood: Modeling of wine uptake and characterization of volatile compounds retention. Food Res. Int. 2019, 116, 249–257. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Carpena, M.; Pereira, A.G.; Prieto, M.A.; Simal-Gandara, J. Wine aging technology: Fundamental role of wood barrels. Foods 2020, 9, 1160. [Google Scholar] [CrossRef]
  9. Arapitsas, P.; Ugliano, M.; Perenzoni, D.; Angeli, A.; Pangrazzi, P.; Mattivi, F. Wine metabolomics reveals new sulfonated products in bottled white wines, promoted by small amounts of oxygen. J. Chromatogr. A 2016, 1429, 155–165. [Google Scholar] [CrossRef]
  10. Gil, M.; Reynes, C.; Cazals, G.; Enjalbal, C.; Sabatier, R.; Saucier, C. Discrimination of rosé wines using shotgun metabolomics with a genetic algorithm and MS ion intensity ratios. Sci. Rep. 2020, 10, 1170. [Google Scholar] [CrossRef] [Green Version]
  11. Ontañón, I.; Sánchez, D.; Sáez, V.; Mattivi, F.; Ferreira, V.; Arapitsas, P. Liquid chromatography-mass spectrometry-based metabolomics for understanding the compositional changes induced by oxidative or anoxic storage of red wines. J. Agr. Food Chem. 2020, 68, 13367–13379. [Google Scholar] [CrossRef] [PubMed]
  12. Laureati, M.; Cattaneo, C.; Tateo, F.; Bononi, M. Identification of the volatile compounds and sensory attributes of long-term aging Vin Santo Wine from Malvasia di Candia aromatic grapes. Foods 2020, 9, 1736. [Google Scholar] [CrossRef] [PubMed]
  13. Wei, Z.; Liu, X.H.; Huang, Y.; Lu, J.; Zhang, Y.L. Volatile aroma compounds in wines from Chinese wild/hybrid species. J. Food Biochem. 2018, 43, 12684. [Google Scholar] [CrossRef] [PubMed]
  14. Chen, L.; Darriet, P. Strategies for the identification and sensory evaluation of volatile constituents in wine. Compr. Rev. Food Sci. Saf. 2021, 20, 4549–4583. [Google Scholar] [CrossRef] [PubMed]
  15. Arena, E.; Rizzo, V.; Licciardello, F.; Fallico, B.; Muratore, G. Effects of light exposure, bottle colour and storage temperature on the quality of Malvasia delle Lipari sweet wine. Foods 2021, 10, 1881. [Google Scholar] [CrossRef]
  16. Yuan, J.F.; Lai, Y.T.; Chen, Z.Y.; Song, H.X.; Zhang, J.; Wang, D.H.; Gong, M.G.; Sun, J.R. Microwave irradiation: Effects on the change of colour characteristics and main phenolic compounds of Cabernet Gernischt dry Red Wine during storage. Foods 2022, 11, 1778. [Google Scholar] [CrossRef]
  17. Phan, Q.; Tomasino, E. Untargeted lipidomic approach in studying pinot noir wine lipids and predicting wine origin. Food Chem. 2021, 355, 129409. [Google Scholar] [CrossRef]
  18. Onzo, A.; Acquavia, M.A.; Pascale, R.; Iannece, P.; Gaeta, C.; Lelario, F.; Ciriello, R.; Tesoro, C.; Bianco, G.; Di Capua, A. Untargeted metabolomic analysis by ultra-high-resolution mass spectrometry for the profiling of new Italian wine varieties. Anal. Bioanal. Chem. 2022, 414, 7805–7812. [Google Scholar] [CrossRef]
  19. Zhang, Q.; Zhang, Y.Y.; Liu, Z.; Zhang, Y.M.; Lu, N.; Hai, G.Q.; Shao, S.Z.; Zheng, Q.X.; Zhang, X.; Fu, H.Y.; et al. DifferentiatingWestlake Longjing tea from the first- and second-grade producing regions using ultra high performance liquid chromatography with quadrupole time-of-flight mass spectrometry-based untargeted metabolomics in combination with chemometrics. J. Sep. Sci. 2020, 43, 2794–2803. [Google Scholar] [CrossRef]
  20. Zhao, J.-J.; Zhang, Y.; Wang, X.-C.; Wang, X.; Zhang, Q.; Lu, P.; Liu, P.-P.; Yu, Y.-J.; Han, L.; Zhou, H.-N.; et al. A new platform for untargeted UHPLC-HRMS data analysis to address the time-shift problem. Anal. Chim. Acta. 2022, 1193, 339393. [Google Scholar] [CrossRef]
  21. Lu, Y.; Yao, G.; Wang, X.; Zhang, Y.; Zhao, J.; Yu, Y.-J.; Wang, H. Chemometric discrimination of the geographical origin of licorice in China by untargeted metabolomics. Food Chem. 2022, 380, 132235. [Google Scholar] [CrossRef]
  22. El-Gindy, A.; Hadad, G.M. Chemometrics in Pharmaceutical Analysis: An Introduction, Review, and Future Perspectives. J. Aoac. Int. 2012, 95, 609–623. [Google Scholar] [CrossRef]
  23. Huang, Y.P.; Wu, Z.W.; Su, R.H.; Ruan, G.H.; Du, F.Y.; Li, G.K. Current application of chemometrics in traditional Chinese herbal medicine research. J. Chromatogr. B 2016, 1026, 27–35. [Google Scholar] [CrossRef]
  24. Vaclavik, L.; Lacina, O.; Hajslova, J.; Zweigenbaum, J. The use of high performance liquid chromatography–quadrupole time-of-flight mass spectrometry coupled to advanced data mining and chemometric tools for discrimination and classification of red wines according to their variety. Anal. Chim. Acta. 2011, 685, 45–51. [Google Scholar] [CrossRef] [PubMed]
  25. Liu, P.P.; Zhou, H.N.; Zheng, Q.X.; Lu, P.; Yu, Y.J.; Cao, P.J. An automaticUPLC-HRMS data analysis platform for plant metabolomics. Plant Biotechnol. J. 2019, 17, 2038–2040. [Google Scholar] [CrossRef] [PubMed]
  26. Cadahia, E.; Simon BF, D.; Sanz, M.; Poveda, P.; Colio, J. Chemical and chromatic characteristics of Tempranillo, Cabernet Sauvignon and Merlot wines from DO Navarra aged in Spanish and French oak barrels. Food Chem. 2009, 115, 639–649. [Google Scholar] [CrossRef]
  27. Fernández De Simón, B.; Cadahía, E.; Conde, E.; García-Vallejo, M.C. Evolution of phenolic compounds of Spanish oak wood during natural seasoning. First results. J. Agric. Food Chem. 1999, 47, 1687–1694. [Google Scholar] [CrossRef]
  28. Figueiredo-González, M.; Cancho-Grande, B.; Simal-Gándara, J.; Teixeira, N.; Mateus, N.; De Freitas, V. The phenolic chemistry and spectrochemistry of red sweet wine-making and oak-aging. Food Chem. 2014, 152, 522–530. [Google Scholar] [CrossRef]
  29. Monagas, M.; Gomez-Cordoves, C.; Bartolome, B. Evolution of the phenolic content of red wines from L. during ageing in bottle. Food Chem. 2006, 95, 405–412. [Google Scholar] [CrossRef]
  30. Gutierrez, I.H.; Lorenzo, S.P.; Espinosa, A.V. Phenolic composition and magnitude of copigmentation in young and shortly aged red wines made from the cultivars, Cabernet Sauvignon, Cencibel, and Syrah. Food Chem. 2005, 92, 269–283. [Google Scholar] [CrossRef]
Figure 1. (a) Comparison of the total ion chromatograms in the positive (Pos) for oak barrel aged (A) and non-oak barrel aged (N) wine samples. (b) Comparison of the total ion chromatograms in the negative (Neg) for oak barrel aged (A) and non-oak barrel aged (N) wine samples.
Figure 1. (a) Comparison of the total ion chromatograms in the positive (Pos) for oak barrel aged (A) and non-oak barrel aged (N) wine samples. (b) Comparison of the total ion chromatograms in the negative (Neg) for oak barrel aged (A) and non-oak barrel aged (N) wine samples.
Chemosensors 11 00165 g001
Figure 2. (a) EIC extraction and detail diagram example in the positive ion mode. (b) EIC extraction and detail diagram example in the negative ion mode. The data quality assessment diagram in the ESI+ mode (c) and the ESI− mode (d).
Figure 2. (a) EIC extraction and detail diagram example in the positive ion mode. (b) EIC extraction and detail diagram example in the negative ion mode. The data quality assessment diagram in the ESI+ mode (c) and the ESI− mode (d).
Chemosensors 11 00165 g002
Figure 3. Multivariate analysis of metabolomics distinguishes the oak barrel aged and non-oak barrel aged wines: (a) The PCA score plot in the positive ion mode; (b) the PCA score plot in the negative ion mode; (c) the PLS-DA score plot in the positive ion mode; (d) the PLS-DA score plot in the negative ion mode.
Figure 3. Multivariate analysis of metabolomics distinguishes the oak barrel aged and non-oak barrel aged wines: (a) The PCA score plot in the positive ion mode; (b) the PCA score plot in the negative ion mode; (c) the PLS-DA score plot in the positive ion mode; (d) the PLS-DA score plot in the negative ion mode.
Chemosensors 11 00165 g003
Figure 4. Example of the compound identification process with AntDAS. (a) The MS1 and MS/MS spectra of the compounds with m/z of 205.0973 were constructed with AntDAS. (b) The EIC plot corresponding to each fragment ion. (c) The final identification result of the compound.
Figure 4. Example of the compound identification process with AntDAS. (a) The MS1 and MS/MS spectra of the compounds with m/z of 205.0973 were constructed with AntDAS. (b) The EIC plot corresponding to each fragment ion. (c) The final identification result of the compound.
Chemosensors 11 00165 g004
Figure 5. Heatmap visualization of the identified metabolites in 66 wine samples for the ESI+ mode (a) and ESI− mode (b).
Figure 5. Heatmap visualization of the identified metabolites in 66 wine samples for the ESI+ mode (a) and ESI− mode (b).
Chemosensors 11 00165 g005
Figure 6. Violin plots of the relative contents of each compound in the wine samples for the ESI+ mode (a) and ESI− mode (b).
Figure 6. Violin plots of the relative contents of each compound in the wine samples for the ESI+ mode (a) and ESI− mode (b).
Chemosensors 11 00165 g006
Table 1. The results of the compounds identified with MS1 and MS/MS constructed with AntDAS in the positive ion mode.
Table 1. The results of the compounds identified with MS1 and MS/MS constructed with AntDAS in the positive ion mode.
No.Experimental MassRT/minMatch Factor (MS1 + MS/MS)FormulaCompound NameMS/MS Compound
Fragments
IonType
1121.06515.8131.24C8H8OPhenylacetaldehyde53.0395, 77.0393
93.0705, 95.0497
103.0547, 105.0452
[M+H]+
2123.04428.1581.86C7H6O24-Hydroxybenzaldehyde53.0395, 67.0550
77.0394, 95.0497
105.0452
[M+H]+
3124.03961.4491.90C6H5NO2Niacin78.0346, 80.0502
96.0450, 124.0396
[M+H]+
4136.06201.2731.86C5H5N5Adenine67.0299, 94.0406
109.0514, 119.0356
[M+H]+
5138.05511.1811.08C7H7NO2Trigonelline65.0394, 78.0344
92.0500, 94.0656
110.0603
[M+H]+
6138.09141.6911.52C8H11NOTyramine77.0393, 91.0548
95.0497, 103.0547
121.0651
[M+H]+
7144.04812.7371.34C6H9NOS4-Methyl-5-thiazoleethanol112.0218, 113.0299
126.0376,
[M+H]+
8144.10201.2301.14C7H13NO2Proline betaine58.0660, 84.0815
98.0970,
[M+H]+
9146.16550.9951.84C7H19N3Spermidine58.0661, 72.0816
84.0815, 112.1126
[M+H]+
10150.05841.4231.98C5H11NO2SMethionine56.0504, 61.0115
74.0245, 87.0270
133.0322
[M+H]+
11156.07751.0371.98C6H9N3O2Histidine56.0505, 82.0533
83.0611, 93.0454
95.0610, 110.0718
[M+H]+
12161.05987.0011.46C10H8O26-Methylcoumarin79.0550, 103.0548
105.0704, 133.0653
[M+H]+
13162.11241.1481.68C7H15NO3L-Carnitine85.0291, 102.0919
103.0395
[M+H]+
14165.05476.1431.80C9H8O3p-Coumaric acid65.0394, 91.0549
119.0496, 123.0444
147.0443
[M+H]+
15166.086410.0982.00C9H11NO2Phenylalanine79.0545, 93.0704
103.0547, 120.0812
[M+H]+
16170.08131.3572.00C8H11NO3Pyridoxine106.0658, 124.0761
134.0603, 150.0558
152.0709
[M+H]+
17175.11901.1312.00C6H14N4O2Arginine60.0565, 70.0659
72.0816, 116.0711
130.0978
[M+H]+
18179.03397.4221.94C9H6O46,7-Dihydroxycoumarin123.0445, 132.0285
133.0285, 151.0397
[M+H]+
19181.04957.4611.60C9H8O4Caffeic acid83.0392, 117.0340
145.0287, 163.0393
[M+H]+
20182.08152.1381.90C9H11NO3Tyrosine91.0549, 95.0497
119.0495, 136.0759
147.0445
[M+H]+
21183.065412.0081.76C9H10O4Syringaldehyde67.0552, 95.0498
105.0453, 123.0444
125.0239, 140.0470
[M+H]+
22185.154122.8531.96C11H20O2delta-Undecalactone149.1329, 150.1363
167.1435, 168.1468
[M+H]+
23189.13481.1311.90C7H16N4O2Targinine57.0456, 70.0660
88.0764, 116.0711
144.1134, 158.0931
172.1083
[M+H]+
24189.16003.8431.98C9H20N2O2Propamocarb58.0660, 74.0245
86.0971, 144.1022
[M+H]+
25192.07685.2321.98C9H9N3O2Carbendazim132.0560, 160.0508[M+H]+
26199.06038.5161.76C9H10O5Syringic acid95.0498, 125.0236
140.0470,181.0502
[M+H]+
27200.118720.3482.00C12H13N3Pyrimethanil82.0659, 107.0609
125.0710,183.0922
[M+H]+
28204.12381.3231.80C9H18NO4Acetylcarnitine60.0817, 85.0291
145.0498,
[M+H]+
29205.09734.4511.94C11H12N2O2Tryptophan118.0657, 144.0810, 146.0603, 170.0605, 188.0710[M+H]+
30206.081311.5531.92C11H11NOIndolelactic acid118.0656, 130.0656
160.0762, 170.0602
[M+H]+
31211.144211.9421.72C11H18N2O2Cyclo(Leu-Pro)70.0659, 98.0606
114.0918, 138.1279
155.1543, 183.1495
[M+H]+
32217.12781.2391.66C8H16N4O3N-Acetylarginine70.0660, 116.0710
130.0978
[M+H]+
33229.086113.1931.86C14H12O3Resveratrol91.0549, 107.0497
119.0496, 135.0444
211.0750
[M+H]+
34231.102118.2111.08C14H14O3α,β-Dihydroresveratrol91.0549, 107.0497
125.0597, 137.0599
[M+H]+
35239.10282.1021.60C11H14N2O4Gly-Tyr91.0548, 123.0442
136.0759, 147.0444
165.0549, 182.0819
[M+H]+
36261.14433.5771.56C11H20N2O5gamma-Glutamylleucine86.0971, 132.1025[M+H]+
37273.076313.4722.00C15H12O5Naringenin chalcone119.0501, 147.0443
153.0185
[M+H]+
38275.092017.1972.00C15H14O5Phloretin107.0497, 169.0495[M+H]+
39281.11365.0591.72C13H16N2O5Asp-Phe120.0812, 130.0622
166.0864, 235.1076
[M+H]+
40289.070716.5031.72C15H12O6Dihydrokaempferol107.0497, 149.0237
153.0185, 215.0709
[M+H]+
41291.08639.7241.98C15H14O6Epicatechin123.0443, 139.391
147.0441, 207.0654
[M+H]+
42291.08646.0042.00C15H14O6Catechin95.0499, 119.0496
123.0444, 139.0392
147.0444, 207.0656
[M+H]+
43291.08656.6091.98C15H14O6(-)-epicatechin119.0496, 123.0444
139.0392, 147.0445
179.0708, 207.0656
[M+H]+
44291.08674.8862.00C15H14O6CIANIDANOL123.0444, 139.0392
147.0443, 165.0552
207.0657
[M+H]+
45294.15481.8621.78C12H23NO7N-Fructosyl isoleucine86.0972, 144.1024
212.1281, 258.1337
276.1449
[M+H]+
46305.065510.0781.98C15H12O7taxifolin123.0443, 149.0239
153.0183, 167.0344
231.0651, 259.0604
[M+H]+
47319.045216.7942.00C15H10O8Myricetin111.0084, 153.0182
245.0452, 273.0398
301.0350
[M+H]+
48328.13922.8791.84C15H21NO7N-Fructosyl phenylalanine97.0289, 264.1228
292.1178, 310.1290
[M+H]+
49345.14494.1282.00C18H20N2O5Tyr-Tyr119.0498, 136.0761
182.0821
[M+H]+
50377.14579.8192.00C17H20N4O6Riboflavin69.0342, 99.0447
243.0879, 359.1360
[M+H]+
51434.202811.5041.64C19H28O10Sayaendoside87.0447, 115.0393
133.0498, 145.0499
[M+NH4]+
Table 2. The results of the compounds identified with MS1 and MS/MS constructed with AntDAS in the negative ion mode.
Table 2. The results of the compounds identified with MS1 and MS/MS constructed with AntDAS in the negative ion mode.
No.Experimental MassRT/minMatch Factor (MS1 + MS/MS)FormulaCompound NameMS/MS Compound FragmentsIonType
1121.02838.1661.94C7H6O2Benzoic acid93.0332, 108.0203[M−H]
2125.02332.7081.82C6H6O3Benzene-1,2,4-triol69.0331, 95.0122, 97.0281[M−H]
3135.02781.1231.84C4H8O5Threonic acid59.0124, 72.9916
75.0072
[M−H]
4137.02335.4941.92C7H6O3Protocatechuic aldehyde108.0202, 109.0281
119.0126, 136.0155
[M−H]
5147.02891.4501.32C5H8O5Citramalate85.0280, 87.0073
129.0189
[M−H]
6149.00795.9831.48C4H6O6Tartaric acid59.0124, 72.9916
87.0074
[M−H]
7151.02511.5931.96C5H4N4O2Xanthine108.0189, 126.0298[M−H]
8161.04451.9781.98C6H10O5Meglutol57.0332, 59.0125
99.0436, 101.0230
[M−H]
9163.03886.1411.98C9H8O3Phenylpyruvic acid93.0332, 117.0334
119.0489
[M−H]
10165.054711.1661.96C9H10O33-phenyllactic acid72.9917, 91.0537
119.0489, 147.0446
[M−H]
11169.01332.2331.88C7H6O5Gallic acid69.0331, 97.0280
124.0153, 125.0230
[M−H]
12177.01867.4121.84C9H6O4esculetin89.0382, 105.0332
033.0285, 149.0233
[M−H]
13181.05004.6991.92C9H10O43-(4-Hydroxyphenyl)lactic acid72.9916, 119.0489
134.0361, 135.0440
163.0390
[M−H]
14191.01901.4411.98C6H8O7Citric acid57.0331, 85.0280
87.0073, 111.0074
[M−H]
15191.05551.1971.84C7H12O61,3,4,5-Tetrahydroxycyclohexanecarboxylic Acid59.0123, 710.0123
85.0280, 93.0332
[M−H]
16193.03491.1231.88C6H10O72-Keto-L-galactonic acid57.0331, 59.0124
71.0124, 73.0280
85.0280, 101.0230
113.0230, 131.0337
[M−H]
17197.04498.5191.80C9H10O5Syringic Acid95.0120, 123.0077
166.9973, 182.0211
[M−H]
18206.081712.4011.88C11H13NO3N-acetylphenylalanine58.0284, 91.0540
147.0454, 164.0707
[M−H]
19209.02971.1471.74C6H10O8Mucic acid71.0124, 85.0280
133.0128, 191.0177
[M−H]
20263.128919.4111.68C15H20O4Abscisic acid139.0760, 151.0753
201.1281, 204.1145
219.1382
[M−H]
21289.07166.1741.90C15H14O6Catechin151.0389, 179.0340
205.0499, 221.0814
227.0709, 245.0830
271.0612
[M−H]
22301.035319.7551.96C15H10O7Quercetin121.0283, 151.0024
178.9976, 273.0348
[M−H]
23303.051313.5011.94C15H12O7taxifolin125.0232, 175.0390
217.0499, 285.0400
[M−H]
24305.06683.5262.00C15H14O72-(3,4,5-trihydroxyphenyl) chromane-3,5,7-triol125.0231, 167.0339
177.0547, 219.0659
261.0769
[M−H]
25317.030616.8001.98C15H10O8Myricetin107.0126, 137.0232
151.0027, 178.9975
[M−H]
26577.13485.7351.86C30H26O12Procyanidin B1125.0231, 161.0234
287.0555, 407.0769
425.0883
[M−H]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lv, Y.; Ma, F.-L.; Wang, J.-N.; Zhang, Y.; Jiang, Y.; Ge, Q.; Yu, Y.-J. A Strategy for Differentiating Oak Barrel Aged and Non-Oak Barrel Aged Wines by Using UHPLC–HRMS Combined with Chemometrics. Chemosensors 2023, 11, 165. https://doi.org/10.3390/chemosensors11030165

AMA Style

Lv Y, Ma F-L, Wang J-N, Zhang Y, Jiang Y, Ge Q, Yu Y-J. A Strategy for Differentiating Oak Barrel Aged and Non-Oak Barrel Aged Wines by Using UHPLC–HRMS Combined with Chemometrics. Chemosensors. 2023; 11(3):165. https://doi.org/10.3390/chemosensors11030165

Chicago/Turabian Style

Lv, Yi, Feng-Lian Ma, Jia-Nan Wang, Yao Zhang, Yuan Jiang, Qian Ge, and Yong-Jie Yu. 2023. "A Strategy for Differentiating Oak Barrel Aged and Non-Oak Barrel Aged Wines by Using UHPLC–HRMS Combined with Chemometrics" Chemosensors 11, no. 3: 165. https://doi.org/10.3390/chemosensors11030165

APA Style

Lv, Y., Ma, F. -L., Wang, J. -N., Zhang, Y., Jiang, Y., Ge, Q., & Yu, Y. -J. (2023). A Strategy for Differentiating Oak Barrel Aged and Non-Oak Barrel Aged Wines by Using UHPLC–HRMS Combined with Chemometrics. Chemosensors, 11(3), 165. https://doi.org/10.3390/chemosensors11030165

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