Next Article in Journal
Use of Bio-Waste of Ilex paraguariensis A. St. Hil. (Yerba mate) to Obtain an Extract Rich in Phenolic Compounds with Preservative Potential
Previous Article in Journal
Strategies to Formulate Value-Added Pastry Products from Composite Flours Based on Spelt Flour and Grape Pomace Powder
Previous Article in Special Issue
Characteristics of Red Deer (Cervus elaphus) Milk: Lactational Changes in Composition and Processing Impacts on Structural and Gelation Properties
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Applications of Solution NMR Spectroscopy in Quality Assessment and Authentication of Bovine Milk

1
Laimburg Research Centre, Laimburg 6—Pfatten (Vadena), 39040 Auer, Italy
2
Faculty of Agricultural, Environmental and Food Sciences, Free University of Bozen-Bolzano, Piazza Unversità 5, 39100 Bolzano, Italy
*
Author to whom correspondence should be addressed.
Foods 2023, 12(17), 3240; https://doi.org/10.3390/foods12173240
Submission received: 6 July 2023 / Revised: 7 August 2023 / Accepted: 18 August 2023 / Published: 28 August 2023

Abstract

:
Nuclear magnetic resonance (NMR) spectroscopy is emerging as a promising technique for the analysis of bovine milk, primarily due to its non-destructive nature, minimal sample preparation requirements, and comprehensive approach to untargeted milk analysis. These inherent strengths of NMR make it a formidable complementary tool to mass spectrometry-based techniques in milk metabolomic studies. This review aims to provide a comprehensive overview of the applications of NMR techniques in the quality assessment and authentication of bovine milk. It will focus on the experimental setup and data processing techniques that contribute to achieving accurate and highly reproducible results. The review will also highlight key studies that have utilized commonly used NMR methodologies in milk analysis, covering a wide range of application fields. These applications include determining milk animal species and feeding regimes, as well as assessing milk nutritional quality and authenticity. By providing an overview of the diverse applications of NMR in milk analysis, this review aims to demonstrate the versatility and significance of NMR spectroscopy as an invaluable tool for milk and dairy metabolomics research and hence, for assessing the quality and authenticity of bovine milk.

1. Introduction

Nuclear magnetic resonance (NMR) spectroscopy has proven to be a powerful tool for the analysis of milk and dairy products, given its unique ability to simultaneously detect and identify small molecules within complex matrices, such as milk [1]. NMR was initially employed by Odeblad and Westing in the 1950s to investigate milk properties and allowed the identification of distinct components of human milk: lactose from milk plasma and milk fat [2]. Since then, NMR spectroscopy has been extensively utilized in various studies, providing valuable insights into the composition, structure, and properties of milk, including the content of key components such as fat, protein, lactose, and minerals. One of the main advantages of NMR spectroscopy in milk analysis is its non-destructive nature, allowing for repeated measurements on the same sample over time and facilitating quality assessment throughout the production process [3]. Additionally, its inherent quantitative nature allows it to be used for quantitative metabolomic analysis of milk.
Metabolomics is an emerging tool in milk research; it gives valuable information about the health status and feeding system of the cow, the milk technological process and shelf-life, and eventually the nutritional quality of the milk. It offers reliable indicators and biomarker compounds that can be used for milk quality assessment and authentication. Over the past two decades, NMR has emerged as one of the most commonly used analytical techniques in milk metabolomic research, alongside gas chromatography coupled with mass spectrometry (GC-MS) and liquid chromatography coupled with single-stage mass spectrometry (LC-MS) [1]. Moreover, NMR-based milk metabolomics is continuously evolving as a promising approach offering high-throughput and comprehensive coverage of metabolites [4,5,6,7,8].
This review aims to highlight the various applications of NMR in the quality assessment and authentication of bovine milk, with particular consideration of metabolomics platforms. While previous literature reviews have covered the applications of NMR in milk and dairy products [4,9] as well as NMR-based metabolomics in their authenticity and quality assessment [5,6,7,10,11], this review places significant emphasis on the most recent advances and applications of NMR, particularly in bovine milk studies. The review will provide a comprehensive summary of commonly used NMR methodologies and NMR data analysis applied in the analysis of bovine milk, along with representative examples from related application fields. Although the primary focus of the review will be on high-resolution NMR, significant examples of low-resolution NMR will also be introduced as a promising tool in milk analysis at the industrial level. Figure 1 outlines the workflow of the review.

2. NMR Principle

Nuclear magnetic resonance (NMR) spectroscopy is considered among the most powerful analytical techniques used to determine the physical, chemical, and biological properties of matter. It allows for the provision of atomic-level information on molecules by measuring the interaction of external magnetic fields with the NMR active nuclei of molecules placed in a strong static magnetic field (B0). In fact, only certain naturally occurring nuclei have intrinsic properties that allow them to be used (and observed) in NMR analysis. Quantum mechanics tells us that those nuclei possess a nuclear spin angular momentum ( I ) characterized by a nuclear spin quantum number (I). Only nuclei with nonzero nuclear spin quantum numbers, I ≠ 0 (meaning isotopes with an odd number of proton and/or neutrons) possess a nuclear magnetic moment, μ , that enables alignment in an external static magnetic field (similar to tiny bar magnets) and the ability to absorb an applied radiofrequency radiation (RF), thus allowing NMR measurements. The most commonly used NMR-active nuclei are those with a half-integer spin quantum number, especially I = ½. Fortunately, most common elements found in living organisms have at least one of the isotopes with spin number I = ½ (i.e., 1H, 13C, 15N, and 31P), thus allowing NMR applications in a wide range of biological samples.
For spin I = ½ nuclei, the application of an external magnetic field B0 aligns the nuclear magnetic moment vector, μ , with or against the B0 field, corresponding to low energy or high energy spin states, respectively. The difference between the two energy levels is:
∆E = γhB0/2π,
Linearly depends on the external magnetic field applied, B0 as well as on the gyromagnetic ratio, which is characteristic of each NMR active nucleus. The Boltzmann distribution defines the difference in population between the two energy levels as:
Nhigh/Nlow = e(∆E/kBT),
where Nhigh and Nlow are the populations of the upper and lower states, respectively, ∆E is the difference between the two energy levels, kB is the Boltzmann constant, and T is the temperature expressed in Kelvin. As for other spectroscopic methods (i.e., UV–vis), the amount of signal observed in NMR Spectroscopy is related to the difference in population between the two energy states. However, differently from UV–vis spectroscopy, at room temperature (T = 298 K), the upper and lower energy spin states are almost equally populated in NMR spectroscopy, with only a very small excess in the lower energy state, explaining NMR’s inherent low sensitivity compared to other techniques.
For practical NMR applications, it is cumbersome to think in terms of individual spins and nuclear magnetic moment vectors ( μ ). Instead, it is more convenient to refer to the sum of all μ in the sample, corresponding to the net magnetization (M). At equilibrium, M is aligned with the main magnetic field (B0). Irradiation of the sample with an RF pulse orthogonal to the main B0 tilts the magnetization M toward the xy plane. In quantum mechanics, irradiation of the sample with an RF pulse (generally applied for a few microseconds) of frequency (in Hz):
v0 = γB0/2π,
which determines the transition of spin between the energy levels (as indicated in Equation (1)). The energy absorbed by the nuclear spins induces a voltage that can be detected by a suitably tuned coil of wire, amplified, and the signal displayed as a free induction decay (FID). By converting the time-domain signal (FID) into a frequency-domain signal using the Fourier transform (FT), the FT provides a high-resolution NMR spectrum that contains valuable information about the molecular structure and dynamics of the sample.
In the absence of any further perturbing RF pulses, two relaxation mechanisms, spin–lattice (T1) and spin–spin (T2) relaxations, eventually return the spin system to thermal equilibrium. The two mechanisms should not be confused and are profoundly different: If the former represents how quickly the net magnetization M returns aligned to B0 (i.e., the system recovers to its ground state), the latter refers to the progressively decaying magnetization M due to dephasing on the xy plane. Both types of relaxation can be measured on several NMR active nuclei by widely used and implemented relaxation-based NMR experiments and provide useful information for molecular structure elucidation and dynamic analyses.
In the following paragraph, the commonly utilized NMR observables in the field of food science, and more specifically in the characterization of milk by NMR spectroscopy, are described.
The chemical shift (δ) determines the position of a signal in the NMR spectrum and conveys information about the nuclei’s chemical environment. The chemical shift is measured in parts per million (ppm) relative to an internal standard, such as tetramethylsilane (TMS) or trimethylsilylpropanoic acid (TSP), which is typically set at 0 ppm. Using ppm instead of Hz makes the chemical shift independent from the spectrometer’s operational frequency, facilitating cross-comparison of data obtained from various spectrometers with different magnetic fields.
Spin–spin scalar J coupling is the result of neighboring spins interacting through bonding electrons, causing signal splitting in NMR spectra. Even though J coupling makes NMR spectra more complicated, it is crucial for structural characterization since its strength depends on the distance and dihedral angle between nuclei. However, splitting due to the J coupling can be removed using “broadband decoupling”, with the benefit of simplifying the spectrum and an increase in the signal-to-noise ratio (SNR).
The area of a given NMR peak is proportional to the number of NMR active nuclei that give rise to the peak, thus making modern NMR spectroscopy a robust method (qNMR) for the quantification of specific compounds. For example, a peak generated by three equivalent protons in a methyl group (CH3) has an area three times greater than the area of a peak caused by a single proton. Thus, one of the greatest advantages of the NMR method applied in food studies is its accuracy in quantitative measurements of the sample constituents because the intensity of the NMR signal is directly proportional to the molar concentration of the constituent [12,13]. Additionally, measurements of two mechanisms involved in the relaxation of magnetization (spin–lattice relaxation, T1, and spin–spin relaxation, T2) and measurements of translation diffusion are based on the determination of signal area.
In the context of this review, we will be mostly focusing on the characterization of milk by 1D and 2D-NMR experiments (see relative section). In the former, the NMR spectrum is typically displayed as a plot of signal intensity (y-axis) versus frequency (x-axis) (expressed as chemical shift and referenced as described previously) measured on one nucleus (generally 1H or 13C), whereas in the latter, cross peaks are defined by two axes (x, y) indicating the chemical shift on two different nuclei. Experiments conducted in two dimensions can either involve the same type of nuclei (homonuclear, i.e., 1H-1H correlation spectroscopy, COSY) or different types of nuclei (heteronuclear, i.e., 1H-13C heteronuclear single quantum coherence spectroscopy, HSQC). Two-dimensional experiments offer improved resolution compared to 1D-NMR since spreading out the spectrum in two dimensions has the advantage for large molecules in that it removes much peak overlap. Additionally, combinations of 2D-NMR experiments (such as 1H-13C heteronuclear multiple bond correlation spectroscopy, HMBC, and 1H-1H correlation spectroscopy, COSY) allow assignment of chemical shifts (i.e., identification of NMR signals attributable to specific atoms).

3. NMR Experiments on Milk

3.1. 1D 1H-NMR Experiments

One-dimensional 1H-NMR spectroscopy is the most preferred NMR method in milk analysis mainly due to the high isotopic natural abundance (∼99%) and short relaxation time of 1H nuclei, leading to higher sensitivity and faster data acquisition [1,4] compared to other NMR experiments. Moreover, 1H-NMR spectra provide quantitative and very detailed structural information about milk constituents; therefore, it is widely performed in metabolomic studies. However, the presence of scalar J coupling between neighboring protons and a narrow chemical shift window (usually from 0 ppm up to ∼15 ppm) leads to an overcrowded spectrum with significantly overlapping peaks . This phenomenon results in low resolution of the 1H-NMR spectra and therefore makes interpretation of the peaks difficult. Resolution is a crucial factor in NMR spectroscopy, and it refers to the ability to distinguish between closely spaced peaks in the NMR spectrum. Furthermore, 1H-NMR experiments carried out on spectrometers with higher magnetic fields resolved this peak overlap by increasing the resolution of the acquired 1H-NMR spectra [1].
1H-NMR has been performed on a very wide range of applications in bovine milk studies. It has demonstrated excellent results on milk chemical composition characterization [14], nutritional quality determination of milk [15,16], discrimination of milk among different feeding regimes [17] and animal species [18], the effects of processing treatments on milk quality and shelf life [19,20,21], and adulteration of milk powder [22,23]. Figure 2 shows an example of the 1D 1H-NMR spectrum of milk that contains about 296 milk metabolites and metabolite species. More detailed information concerning the study is given in Section 7.

3.2. 1D 13C-NMR Experiments

One-dimensional 13C-NMR spectroscopy is commonly used in the analysis of milk lipid composition since it provides better resolution compared to 1H-NMR [9], and thus it is more informative in the identification and quantification of lipids [4]. In fact, 13C-NMR has a broader chemical shift window (∼200 ppm) than 1H-NMR (∼15 ppm), and therefore it has fewer overlapping peaks in the spectrum and, as described before, consequently higher resolution [24]. However, it is well known that quantitative analysis (qNMR) of 13C signals presents some complexities, in particular, related to (a) low sensitivity due to the low isotopic natural abundance of 13C (1%) (compared to other 100% abundance nuclei as 1H, 19F, or 31P) that require long experiment time (~hours), (b) long T1 relaxation time of 13C nuclei, and consequently long recycling delays (D1). Additionally, the degree of signal enhancement caused by the nuclear Overhauser enhancement (NOE) is proportional to the number of hydrogens attached to each 13C (except for quaternary carbons, which do not contribute to the enhancement) [1,25].
Therefore, the application of 13C-NMR, particularly in large-scale metabolomic studies of milk, has been limited. Nevertheless, it has been used as a very promising tool in milk lipid studies to discriminate milk among animal species [26], among different feeding regimes (organic and conventional production) [27,28,29], and among geographical origins [30].

3.3. 1D 31P-NMR Experiments

One-dimensional 31P-NMR experiments allow high resolution and sensitivity due to the 100% natural abundance of the 31P nucleus [25]. 31P-NMR spectroscopy has proved to be an effective technique in milk analysis to study phosphorus-containing compounds such as phospholipids, phosphorylated carbohydrates, inorganic phosphate (Pi), and modified amino acids such as phosphoserine (SerP) and the residue of casein [9]. 31P-NMR has also been applied to study the structure of caseins and casein micelles in milk [31]. It can be used to quantify and distinguish the different phosphorylated compounds present in milk as well as the diphosphates added to milk as stabilizers [32]. However, quantitative 31P-NMR analysis has a significant downside due to the lengthy T1 relaxation of the phosphorous nucleus, thus requiring time-consuming experiments to reach full relaxation of the magnetization, which is crucial for precise signal integration.
31P-NMR spectroscopy in milk analysis provides valuable information about the phospholipid composition of milk, and therefore it has been successfully employed for bovine milk authentication studies. Indeed, 31P-NMR phospholipid fingerprinting [33,34] and the metabolite profile of phosphorylated compounds [31] in bovine milk and milk from other animal species have been widely studied. 31P-NMR spectroscopy was also employed to investigate the effects of heat treatments on milk mineral equilibria and on the interaction between minerals (Pi and Ca) and micellar casein [35]. In general, the structure and structural changes of milk proteins, specifically caseins and casein micelles, can be investigated by 31P-NMR. In fact, this method allows one to observe signals (although broad) from casein micelles since the observable resonances belong to the phosphate of SerP present only in the micelles. Figure 3 represents the 13C and 31P-NMR spectra of bovine milk with the assignments of the most representative milk metabolites.

3.4. Two-Dimensional (2D) NMR Experiments

2D-NMR can potentially cover the main drawbacks of the most commonly used 1D-NMR experiments in milk studies, where the analysis often suffers from spectral congestion due to its low spectral resolution and peak overlapping or extensive signal splitting due to J-coupling. Additionally, 1D-NMR may not always be sufficient for the quantification of complex molecular structures, such as macromolecules in milk. On the other hand, 2D-NMR experiments are very efficient in the identification of both known and unknown metabolites, while 1D-NMR experiments are mostly limited to the identification of known metabolites [1]. Therefore, 2D-NMR seems more suitable for milk metabolomic studies. However, it is important to note that data acquisition and processing of 2D spectra can be more time-consuming and require skilled NMR spectroscopists with expertise in interpreting the obtained multidimensional data as well as in conducting metabolite identification and quantification.
Two-dimensional NMR has been used in milk analysis mainly for structural elucidation [9,38] and more recently also for the quantification of milk constituents [39]. Among 2D-NMR experiments, heteronuclear single quantum coherence spectroscopy (HSQC) and heteronuclear multiple bond correlation spectroscopy (HMBC) are the most frequently conducted in milk studies. In the former, the chemical shift of 1H is correlated to the directly bound 13C nucleus, whereas the latter experiment gives correlations between 1H and 13C nuclei that are separated by multiple bonds, usually up to three. Since HMBC spectra eliminate one-bond C–H correlations, it can be used for the assignments of quaternary and carbonyl carbons that are not detected in HSQC [1]. Indeed, these two 2D-NMR experiments are commonly applied as complementary experiments in complete assignments and chemical structure elucidation. 1H-13C HSQC experiments can adequately resolve signals from more analytes since extra dispersion is gained when the second dimension, 13C, is introduced [40]. Note that single-quantum HSQC experiments are generally preferred over multi-quantum HMQC experiments since they allow higher resolution on the second dimension despite being more susceptible to losses in SNR [41] An example of the 1H-13C HSQC spectrum of whole bovine milk with the most representative metabolite assignments is shown in Figure 4.
Two-dimensional heteronuclear NMR has been applied in milk studies for identification and quantification of milk constituents [40] and milk fatty acid content, aiming to overcome the sensitivity and resolution limitations of 1D-NMR in the detection of minor molecules [34]. Thus, it showed good results on bovine milk studies of authentication [42], geographical origin [43,44] and milk technological properties [45,46] . Although rare, the use of 1H-31P and 1H-15N HMBC was also reported in the literature for the assignment of NMR signals for lecithin [36].

4. Advantages and Disadvantages of NMR in Milk Studies

NMR spectroscopy represents a robust, cost-effective, and versatile technique for providing quantitative information on the metabolic profile at the molecular level. NMR analysis in milk studies does not require extensive sample preparation or chemical derivatization prior to analysis [47,48]. This unique advantage of NMR allows quick analysis, from a few minutes to a few hours for 1D 1H and 2D-NMR experiments, respectively. Additionally, the nondestructive nature of NMR allows repeated experiments with a high level of reproducibility as well as the possibility of performing further experiments on the same samples using additional analytical techniques [49].
Furthermore, NMR became an ideal technique in untargeted metabolomics over other techniques due to its unique strength to identify unknown compounds in milk that can be potentially identified as biomarkers in studies of milk authentication. Thus, it allows for identifying several compounds in milk, such as sugars, organic acids, vitamins, nucleotides, and aromatic compounds [6] that are less tractable in mass spectrometry-based methodologies. Moreover, the minimal sample preparation, without the need for separation steps such as liquid chromatography, allowed direct and reliable determination of analytes in milk samples by NMR and therefore facilitated the application of NMR, especially in large-scale metabolomic studies. Altogether, these unique strengths of NMR propose this technique as a complementary tool to mass spectrometry (MS), which has been widely performed in milk metabolomics coupled with liquid chromatography (LC-MS) or gas chromatography (GC-MS) [1,6].
However, there are several downsides to NMR spectroscopy. The most significant limitation of NMR is its inherently lower sensitivity compared to MS-based techniques. This is a consequence of the low spin polarization of the nuclei being studied (as described previously) as well as the low natural abundance of some nuclei of interest (for example, ~1.1% for 13C). The limit of detection of LC-MS and GC-MS methods is usually considered 10 to 100 times lower compared to NMR [1]. Moreover, even though NMR analysis is generally cheaper than MS analysis (mostly due to the minimal sample preparation steps required), high-field NMR spectrometers are expensive. It is mainly because of the associated running costs (i.e., periodic replenishment of the cryogenic liquids nitrogen and helium, necessary to keep the magnet at ~−271 °C) and the high investment costs of the instruments [6,11].
In addition to these instrumental challenges, the application of different extraction protocols and solvents or different data processing tools precludes any direct comparison of NMR spectra acquired on the same set of foodstuff samples, impeding any further statistical analysis of the NMR data [8,50]. The application of NMR in food analysis is relatively new; therefore, an international collaboration on standardization of data processing, software tools, and database maintenance is still needed. For instance, NMR-based metabolomics is far from having a comprehensive database; in fact, spectral libraries on food are limited with respect to MS-based metabolomics. Table 1 summarizes the advantages and disadvantages of widely used NMR methodologies for the analysis of bovine milk.

5. Experimental Conditions

The experimental workflow, from sample collection to data analysis, for the application of NMR spectroscopy in the study of milk content and composition is described here in detail. These are key steps to follow in the experimental setup to obtain high throughput and reproducible results. A summarized description of the experimental conditions is reported in Table 2.

5.1. Sample Collection and Pre-Treatment

Between collection and sample preparation, milk samples are generally stored at −80 °C. It is recommended to conduct an NMR analysis on the samples immediately after thawing. Even though storing samples at −80 °C after undergoing freeze–thaw cycles did not result in any noticeable impact on the NMR outcomes, it is not the ideal method for ensuring accurate analysis [1].
The quality of NMR experiments is directly affected by sample preparation, which makes it one of the crucial steps for a complete NMR analysis [60]. Several studies have demonstrated that sample preparation has a significant impact on the metabolite profile as well as on the spectral quality [13,29,61,62]. For instance, a study noted that methanol extraction protocols are more efficient concerning both spectral quality and metabolite extraction compared to ultrafiltration, skimming, and unprocessed milk protocols [62]. Another study showed the importance of optimizing sample preparation in milk lipid extraction (i.e., the extraction method, the ratio of the solvents used in the binary mixture, and the volume and quantity of the sample) to achieve an efficient extraction yield [29]. In fact, the application of accurate sample preparation protocols can ensure a more straightforward NMR output, especially for large-scale metabolomic studies where a large number of samples have to be analyzed. Overall, successful sample preparation should produce good-quality spectra with high SNR, low inter-sample chemical shift variability, and the highest number of metabolites detected in the sample [63].
However, the presence of sugars or other macromolecules such as lipids and proteins in milk is challenging due to different concentrations and polarities, which have an impact on the NMR analysis. These molecules give rise to broad signals in NMR spectra, which lead to overlapping or masking the signals of low-concentration metabolites, thus hampering accurate identification and quantification [6]. Therefore, in studies focusing on secondary milk metabolites, macromolecules are typically removed through sample preparation methods. The commonly used method for this purpose is precipitation. Proteins and lipoproteins are precipitated with methanol and ethanol, and lipids with chloroform [13]. Overall, these methods aim to obtain a milk extract with the highest recovery of metabolites and the highest purity from proteins and lipids [6].
In NMR-based milk metabolomic studies, several further sample preparation protocols have been described in the literature. These protocols describe the use of raw milk without any pre-treatment or additive [36,38], freeze-dried milk [29,43,54]or milk powder to be reconstituted prior to NMR analysis [17,18,64].
A quick approach for fatty acid profiling of milk samples requires the use of milk fat extract. Folch and Bligh–Dyer techniques are the most widely used milk lipid extraction methods, based on a chloroform:methanol:water mixture with a chloroform:methanol ratio of 2:1 (v/v) for extraction [54]. After the phase separation, the lower (organic) layer containing milk lipids is collected and evaporated under nitrogen [20].
The aqueous fraction of milk samples is widely studied in NMR analysis. Skimmed milk was obtained by centrifugation [17,21,25,64], by filtration (using a 10 kDa centrifugal filter unit) [14,22] or by ultra-filtration (using a 3 kDa centrifugal filter unit), which removes both lipid and protein compounds [23,24,26]. Alternatively, chemical precipitation of milk fat was achieved by the addition of chloroform [27,28], deuterated chloroform [32], and dichloromethane [29,30,31] followed by centrifugation. Aqueous fractions were also obtained with ultrafiltration [23,24,26], ultracentrifugation at 74,200× g for 60 min at 4 °C, followed by ultrafiltration [29,30], ethanol addition [25,33], or methanol addition [21] for protein denaturation and followed by a centrifugation step. The supernatant is generally dried under a stream of nitrogen and can be stored at −80 °C until further use.

5.2. NMR Sample Preparation and Experimental Setup

NMR analysis typically involves a solution of the sample in an appropriate NMR solvent. In the case of 1H-NMR, a deuterated solvent is used to facilitate locking and shimming and suppress the large solvent signal, which otherwise could dominate the NMR spectrum [65]. The NMR solvent usually contains an internal reference mainly for (a) chemical shift referencing of 1H, 13C, and 31P nuclei and (b) allowing quantitative measurements. The choice of the internal reference strongly depends on the NMR solvent used. The most common internal references for polar solvents such as deuterium oxide (D2O) are 3-(trimethylsilyl) propionic-2,2,3,3-d4 acid, sodium salt (TSP), and 4,4-dimethyl-4-silapentane-1-sulfonic acid (DSS), whereas tetramethylsilane (TMS) is an ideal internal reference for non-polar solvents such as deuterated chloroform (CDCl3) [5].
Therefore, (1) milk fat extract obtained as described previously is dissolved in CDCl3 containing TMS as a chemical shift and intensity reference suitable for NMR measurements. Alternatively, (2) upon protein and/or fat precipitation, skimmed milk samples are prepared in a D2O phosphate buffer solution containing TSP as an internal reference. Finally, adequate volumes (at least 600 µL in a 5 mm tube to ensure good shimming and lineshape and consequently good spectral quality) of samples are transferred into the NMR tube for analysis. Furthermore, NMR samples can be repeatedly analyzed, stored under controlled conditions to avoid sample degradation, and reanalyzed to validate or confirm earlier findings [1].
An internal reference should have the following aspects: No volatility, complete solubility in the solvent, no chemical reaction with the solvent and/or sample, and no overlapping signal (at least one peak should be identified). For example, TSP and DSS have been shown to interact with proteins and fatty acids in the sample, and therefore they must be carefully evaluated, especially in quantitative NMR analysis [34]. If CDCl3 is used as the main solvent, the capillary of the NMR tube should be sealed prior to NMR analysis, e.g., with parafilm, since CDCl3 is volatile (boiling point at 60.9 °C) and evaporation can preclude quantitative measurements. In addition, it is essential to adjust the pH of the milk sample to prevent any possible chemical shift variations in the NMR spectra caused by slight pH differences between the samples. This issue can be solved by carrying out all NMR experiments in phosphate buffer (usually in a range of pH ~7.0 to 7.4) [35,36] and by further optimization during data analysis using an appropriate spectral alignment method [37].

NMR Experimental Setup for Quantitative NMR (qNMR)

Quantitative NMR (qNMR) spectroscopy is a straightforward analytical method with simultaneous identification and quantification of more than one analyte in a complex mixture in a relatively short experimental time (typically a few minutes for 1D 1H experiments) [66]. NMR-based milk metabolomics requires the acquisition of qNMR spectra. In milk analysis, 1H-qNMR is the most widely applied method since it provides high-throughput quantitative information on milk constituents with high accuracy.
The absolute concentration of the component of interest can be directly calculated by the ratio of the peak area of the component of interest and the peak area of the reference material. Besides internal reference, which is most commonly used in NMR analysis, chemical or electronic external reference (acquired under identical conditions with the sample in different NMR tubes) or residual peaks from the solvent can be used as a reference for quantification [67]. Moreover, the reference material in qNMR can also have a different chemical structure from the analyte studied since there is no requirement for identical reference standards, such as in chromatography-based techniques [48,60].
For qNMR spectra, the following experimental conditions during spectra acquisition have been most widely used in the literature: 90° excitation pulse and repetition time (i.e., the sum of the FID acquisition time and any additional relaxation delay inserted prior to any RF pulse) at least five times the longest T1 to ensure sufficient relaxation of the magnetization [12,24]. This strategy produces the highest SNR in a single scan. Additionally, the length of the 90° pulse, the acquisition time (i.e., the length of the FID), the receiver gain, and the quality of shimming and tuning, all significantly influence the overall quality of NMR spectra [68].

5.3. 1H-NMR Spectra Acquisition: Water Signal Suppression

The quality of 1H-NMR spectra in NMR-based metabolomics significantly depends on the water suppression scheme, which is considered one of the critical steps of 1H-NMR signal acquisition. For instance, the presence of a residual water signal in the 1H-NMR spectra might overwhelm the detection of less concentrated milk metabolites that resonate close to H2O as well as determine spectral baseline distortions [69]. Although the water content can be mostly eliminated by lyophilization, this process elongates the sample preparation time [25], and the remaining water in the milk sample is in most cases not negligible for the quality of the NMR spectra [12].
Modern NMR instruments allow water signal suppression by implementing several pulse schemes; the most common is presaturation [69]. WATERGATE and Excitation Sculpting [70,71], and one-dimensional nuclear Overhauser effect spectroscopy (1D-NOESY) [72]. Even though there are several options available, the 1D-NOESY sequence has emerged as the most commonly used in NMR-based metabolomic studies [61]. Using a combination of presaturation and excitation pulses, 1D-NOESY allows minimal baseline distortion in the NMR spectra compared to other methods [25,73]. Various NMR data processing software such as TopSpin (Bruker, San Jose, CA, USA), Delta (JEOL, Tokyo, Japan), Mnova (Mestrelab Research, Santiago, Spain), etc. [74] also allow removing the water signal post-acquisition. Each tool has a special algorithm with different efficiencies for water signal removal. Therefore, the most appropriate technique must be chosen depending on the focus of the research. Spectral processing steps such as the choice of zero filling, apodization, chemical shift referencing protocol, phase, and baseline corrections are of particular importance in NMR spectra treatments [68]. To ensure consistency and standardization of the process, it is essential that all spectra are treated in an identical manner and, in particular, if some steps are carried out manually, by the same operator. Nonetheless, processed NMR spectra may still suffer from chemical shift drift due to instrumental variations such as temperature and magnetic field oscillation, baseline distortion [63], or variations among samples such as differences in pH or ionic strength. While instrumental variations that result in an overall spectrum-to-spectrum chemical shift drift can be solved by using an internal reference peak, it is more difficult to handle sample-borne variations that result in local peak-to-peak chemical shift drifts [75].

6. Spectral Processing and Data Analysis

Following NMR spectra acquisition, proper spectral processing is applied (i.e., phase correction, zero-filling, baseline correction, and chemical shift referencing) to the NMR data. Subsequently, NMR data processing is required prior to multivariate data analysis to correct possible variables in the dataset and extract the most useful information for the study. Several studies have provided detailed guidelines and recommendations on NMR data analysis, covering spectral processing, data processing, and multivariate data analysis [68,75,76,77,78,79]. The present review will briefly highlight the essential steps of NMR data analysis to provide a comprehensive overview of NMR applications.

6.1. Spectral Processing

Phase and baseline corrections are the initial and most important steps in NMR spectral processing. Baseline correction is essential for removing baseline offsets and distortions caused by the imperfect magnetic field homogeneity of the NMR instrument. It is usually performed via semiautomatic tools with manual identification of relevant baseline regions followed by automatic completion of the remaining baseline correction [68]. Frequency domain correction methods have been more widely used compared to time domain methods because they provide more accurate correction by constructing the baseline curves in the spectra directly [80]. The quality of phase and baseline corrections can affect the accuracy of peak integration and quantification. These processing steps are critical, especially in metabolomic studies where low-concentration metabolites significantly contribute to the results. For instance, a proper baseline correction allows accurate assignments of low-intensity peaks that are sensitive to baseline distortions [80]. Moreover, manual phasing is often needed in metabolomics for a more accurate identification of low-intensity peaks [68]. Consequently, it is important to choose appropriate spectral processing parameters to obtain reliable and accurate results.
Nowadays, spectral processing methods are applied automatically to minimize the influence of the operator on the results and simplify the processing of large sets of NMR spectra [81]. A plethora of software packages are available for spectral processing and further data analysis. Previously mentioned TopSpin (Bruker, San Jose, CA, USA) and Delta (JEOL, Tokyo, Japan) software, also used for NMR spectral acquisition, allow optimal spectral processing. Other common NMR processing software packages include NMR Suite (CHENOMX, Edmonton, AB, Canada), Mnova (Mestrelab Research, Santiago, Spain), and Amix toolkits (Bruker, San Jose, CA, USA) [76]. However, milk studies often require more sophisticated multivariate data analysis due to the complexity and variability of milk chemical composition.
The NMR spectra of milk often have overlapping signals that will affect the accuracy of identification and quantification of studied compounds; therefore, conventional peak integration is mostly inadequate [82]. A more advanced spectral deconvolution approach has been alternatively applied for the integration of these peaks, disentangling overlapping signals into their individual spectral patterns [83]. Spectral deconvolution can be performed manually or automatically. Thus, there are various deconvolution tools using different quantification algorithms available for both 1D and 2D-NMR spectra [1].

6.2. Data Processing

NMR spectra are initially aligned to the internal reference in each spectrum to optimize the possible chemical shift drifts, as described previously. Further, a spectral alignment method is needed for an advanced correction of the peak shifting. Binning (also known as bucketing) is the most commonly used spectral alignment method. The binning process involves dividing the NMR spectrum into a set of equally sized intervals (between 0.01 and 0.05 ppm range) or bins, along the chemical shift axis [63,84]. The spectral intensity in each bin is then integrated to generate a set of bin intensities, which are adjusted (i.e., normalized) to achieve a common scale or baseline across the spectra and then used as input for further analysis such as multivariate statistical analysis, machine learning, or database searching. To optimize the binning process, it is important to choose an appropriate bin size that balances spectral resolution and information content and to carefully evaluate the effect of binning on the downstream analysis. In fact, the main drawback of this method is the loss of both spectral resolution and data that may be relevant to the study, due mainly to the uncontrolled data reduction [85]. Therefore, more sophisticated alignment methods such as dynamic time warping (DTW) [80], correlation optimized warping (COW) [86], and interval correlation optimized shifting (icoshift) [84] have emerged to obtain better results from NMR spectra. After NMR spectra are aligned, normalization is applied to correct the variations in signal intensities among the samples. These variations occur due to the possible differences in sample concentrations, dilutions, and instrumental signal instability, which are very common in the analysis of milk samples [29]. Normalization of NMR spectra refers to the process of scaling the spectral intensities to a common reference point to correct differences in sample concentration and instrument settings [12]. The most common method of normalization is based on total spectral area normalization, where the total area under the NMR spectrum is used as the reference point. This method assumes that the total number of protons is constant and that there are no major changes in the sample matrix or experimental conditions [33].
However, in the case of milk, normalization to total intensity would indeed normalize to lactose concentrations, which is the most abundant metabolite in milk [6,62]. This could lead to misleading results if the goal is to compare the relative concentrations of other metabolites in different milk samples. In such cases, alternative normalization strategies can be used to adjust for the differences in the metabolite concentrations. The standard approach is to normalize the concentration of an internal reference that is added to the sample before NMR analysis [87]. The internal standard should have a known concentration and ideally have minimal overlap with the sample metabolites in the NMR spectrum. Normalization to the internal reference can effectively adjust for variations in sample preparation, NMR measurement, and other experimental factors that could affect the spectral intensity.
Prior to data analysis, data scaling is essential to reduce the spread in the variables where the higher intensity peaks have a higher influence in the multivariate data models [88]. The data scaling allows for an increased contribution from low-concentration peaks in the model. This step is crucial in metabolomic studies to ensure accurate results. Among various scaling methods [6,88,89], Pareto scaling is commonly used in NMR-based milk metabolomics due to its ability to maintain the data close to their original values while reducing the misleading effect of large variances [29]. Thus, to achieve optimal results in NMR data processing, a combination of autoscaling, Pareto, and centering scaling methods is recommended for identifying the discriminative metabolites between two groups [90].

6.3. Metabolite Identification

Solution NMR spectroscopy allows for the identification, quantification, and detailed structural information of metabolites present in bovine milk. One-dimensional and 2D-NMR experiments are commonly employed to provide detailed information about the chemical shifts and coupling patterns of the metabolites. The metabolite assignments are confirmed through multiple experiments, including 2D-NMR experiments such as COSY, TOCSY, and HSQC previously described in Section 3.4. The chemical shift values obtained from the NMR spectrum are compared with databases and reference spectra containing known metabolite chemical shifts. The most common databases that provide valuable reference data for metabolite assignments are the Biological Magnetic Resonance Data Bank (BMRB) [91], Human Metabolome Database (HMDB) [92], the Madison-Qingdao Metabolomics Consortium Database (MMCD) [93], and NMRShiftDB2 [94]. Commercial software packages, such as Chenomx NMR Suite (Chenomx, Edmonton, AB, Canada) and KnowItAll Metabolomics (BioRad Corporation, Hercules, CA, USA), rely on standard spectral libraries and proprietary databases that are created and maintained by the respective companies.

6.4. Multivariate Data Analysis (MVDA): Unsupervised and Supervised Methods

NMR-based metabolomic datasets obtained from milk may contain a large number of samples and variables, making it challenging to extract meaningful information without the application of advanced multivariate data analysis (MVDA) tools. Unsupervised MVDA methods, such as principal component analysis (PCA) and hierarchical cluster analysis (HCA), are commonly employed in NMR-based milk metabolomics. These methods are used to identify the inherent patterns and groupings within the samples without prior knowledge of the sample class or treatment. Briefly, PCA can be used to reduce the complexity of large NMR datasets by identifying the most significant variations in the sample set and projecting the data onto a lower-dimensional space. HCA, on the other hand, groups samples into clusters based on similarities in their spectral profiles and can be used to identify outliers or subgroups within the data. Both methods not only reveal patterns and groupings within NMR-based milk metabolomic data, but they also provide valuable insights to predict the values of target variables of interest, such as milk composition, metabolic profiling, and processing quality [7,63]. Contrarily, supervised techniques rely on prior knowledge of class membership to construct classification models that can predict the class membership of unknown samples [88]. Supervised MVDA methods, such as partial least square discriminant analysis (PLS-DA) and orthogonal variants of PLS-DA (O-PLS-DA), are the most commonly applied data analysis methods in milk studies.
In NMR-based studies on milk, PCA is usually performed prior to any other supervised MVDA methods as an initial exploratory tool to gain insights about the dataset [12]. PCA is often combined with PLS-DA, particularly in milk quality assessment and adulteration studies. A typical multivariate data analysis process in milk analysis follows a metabolic approach, where the metabolite content significantly differs between the sample groups according to the NMR spectra and is used as input variables into subsequent multivariate analysis [25,30,95]. Followingly, PCA loading plots are used to determine whether those characteristic variables identified in the NMR spectra can distinguish the sample groups. Variable selection is then performed to identify which variables contribute most to group discrimination. Supervised PLS-DA is subsequently conducted, and variable importance in projection (VIP) scores are calculated to determine the significance of the studied variables. Finally, milk metabolites with a VIP value higher than 1 are usually selected as potential biomarkers for distinguishing the sample groups [54]. Figure 5 shows the PCA score plot, PLS-DA score plot, and OPLS-DA corresponding color-coded plot from a study that investigated bovine milk adulteration in caprine milk [96]. In the study, four NMR signals originating from orotate, citrate, N-acetyl carbohydrates, and one unknown compound were detected as more abundant in bovine milk, and they were identified as potential biomarkers for differentiating the two groups and estimating the adulteration percentage in the bovine–caprine milk mixture. Overall, the multivariate model demonstrated clear clustering between the two groups.

7. NMR Applications on Bovine Milk: Quality Assessment

NMR spectroscopy has been widely applied in studies of milk quality assessment. Various studies have shown the capability of NMR in determining the composition of bovine milk [14,36], distinguishing milk from different animal species [26], different farms [44], and discerning milk from different feeding systems [17,95] Moreover, NMR analysis enables the evaluation of the nutritional value of milk [25,59], monitoring the milk quality during processing and storage [14,25,48,63], and assessing the technological properties of milk [45,46]. Table 2 provides a summarized description of the application fields of NMR studies on bovine milk, including the NMR experiments and types of samples used for the corresponding investigation.

7.1. Milk Composition

The pioneering study by Hu et al. [36] in 2004 holds a significant place in the literature as the first investigation to characterize the chemical composition of whole bovine milk using NMR spectroscopy without any additive or sample pretreatment. In the study, several milk metabolites were assigned to perform both one-dimensional (1H, 13C) and two-dimensional (1H-13C, 1H-15N, 1H-31P) NMR experiments. Notably, they identify two compounds, creatine and N-acetyl carbohydrates, for the first time in milk through NMR analysis. It was observed that the addition of 10% D2O assisted in adjusting the lock and shims, thereby improving the resolution and sensitivity of the NMR spectra. However, the study demonstrated that even without pretreatment, milk can be analyzed for the assignment and quantification of a large number of metabolites [36]. Additionally, the study revealed that non-homogenized milk showed broader signals in the 1H-NMR spectra compared to homogenized milk. This broadening effect was attributed to the presence of larger fat globules in non-homogenized milk, which results in signal broadening and reduced NMR signal sensitivity. This finding later influenced NMR-based milk analysis, allowing for the monitoring of milk fat globule size and stability, as well as the assessment of milk homogenization quality. In their following study published in 2007 [38], fatty acids and various milk compounds were successfully quantified by 2D-NMR analysis of whole milk with the addition of 10% D2O. The authors stated that the study has been reported as the first study using 2D-NMR for the quantification of food components [38]. Overall, these studies have contributed to the development of the application of NMR in milk studies.
Foroutan et al. [16] conducted a comprehensive study on the chemical composition of bovine milk, employing a multi-platform approach that encompassed targeted and quantitative metabolomics through NMR, LC-HRMS, LC-MS/MS, and inductively coupled plasma-mass spectrometry (ICP/MS). This study, conducted in 2019, reported the largest number of metabolites quantified in bovine milk at that time, including the identification of some new compounds in milk such as lysophosphatidylcholine (18:2), phosphatidylcholine (28:1), and triglyceride (48:3) [14]. In total, 296 milk metabolites and metabolite species were quantified and validated using both quantitative metabolomics and literature mining approaches. The authors successfully identified 2355 unique metabolite structures in milk, and the data are readily accessible through a web-accessible database called the Milk Composition Database (MCDB).

7.2. Milk Origin

The pioneering NMR technique has been successfully applied to distinguish milk types according to their animal species origin. Andreotti et al. [26] managed to differentiate bovine milk from sheep and goat milk through fatty acid composition using 13C-NMR spectroscopy combined with fuzzy logic analysis. Moreover, Garcia et al. [33] optimized a 31P-NMR fingerprinting approach (with validation on GC) to differentiate milk samples from different animal species (bovine, human, camel, and mare) through their phospholipid (PL) profiles. This method enabled the identification and quantification of multiple PLs in milk samples, facilitating the discrimination of milk according to their animal origin [33]. Bruschetta et al. [31] also highlighted the effectiveness of 31P-NMR fingerprinting in differentiating milk samples from various animal species (e.g., cow, sheep, goat, mare, and donkey). They observed similar 31P-NMR profiles among milk from the same animal species, irrespective of the breed, age, or lactation of the animal [32]. Moreover, Wei et al. [34] conducted a comparative study on the chemical composition of milk from five different mammals, including bovine, human, goat, yak, and donkey. They investigated the PL composition of milk samples using 31P-NMR and LC-MS, as well as the fatty acid composition using GC. The study revealed variations in the concentrations and composition of PLs and the size of fat globules among milk samples from different animal groups. Interestingly, human milk exhibited a fat globule size similar to that of bovine milk. These findings suggest that NMR analysis could be utilized to optimize infant formulas to more closely resemble human milk in terms of PLs and milk fat globule structure, thereby ensuring the nutritional quality of such formulas [34].

7.3. Milk from Different Feeding Systems

Recent studies have utilized NMR spectroscopy to investigate the relationship between milk cow feeding patterns and milk composition and quality. O’Callaghan et al. [51] conducted a study comparing pasture feeding systems to indoor feeding systems and their effects on the milk metabolome using 1H-NMR metabolomics combined with multivariate data analysis. In the study, milk samples were collected from three commonly practiced feeding systems in Ireland: Pasture perennial ryegrass, perennial ryegrass with white clover, and an indoor total mixed ration consisting of maize and grass silage. The milk samples from cows on the two different pasture feeding systems exhibited similar metabolomes throughout the entire lactation period. In contrast, a significant difference was observed between the milk samples from cows in the indoor silage feeding system and the pasture feeding system. For instance, milk samples from pasture feeding had a higher hippuric acid content, while milk from silage feeding had higher urea content [51]. The increased concentration of urea in milk can have a negative impact on milk protein quality, as it is the main source of non-protein nitrogen content in milk [51]. Thus, the authors stated that hippuric acid could serve as a biomarker for milk from pasture-based feeding regimes since its presence in milk is correlated with forage feeding. Conversely, Lanza et al. [17] reported a higher concentration of hippuric acid in milk from maize silage feeding compared to hay feeding. In the study, both the milk polar fraction and fatty acids of milk samples from three different feeding regimes were analyzed using 1H-NMR and GC-MS, respectively. Milk samples from lucerne hay feeding were found to have considerably higher levels of polyunsaturated fatty acids (PUFA) compared to those from maize silage or maize silage partially replaced with grass–legume silage feeding systems. Furthermore, the NMR milk polar metabolomic profile was found to be less sensitive to alterations in the cow’s diet compared to the milk FA profile. Indeed, the authors noted that the overall metabolomic profile of milk remained largely unaltered when maize silage was partially substituted with a grass–legume silage mixture. However, a significant change was observed only when maize silage was completely replaced with hay [17]. These studies highlight the reliability and precision of NMR in assessing the impact of feeding practices on milk quality.

7.4. Milk Technological Process

Recent NMR techniques have shown significant potential for monitoring milk quality and shelf life during various processing applications and subsequent storage. D. Zhu et al. conducted studies addressing two relatively unexplored aspects in the field: The impact of vat pasteurization [54] and the effect of freeze-drying [19] processes on milk quality under different storage conditions. Both studies employed untargeted metabolomics using 1H-NMR and MS-based techniques, followed by multivariate analysis. In the study on vat pasteurization, the authors observed no significant changes in the milk metabolite profile following the pasteurization process. However, during subsequent refrigerated storage (4 °C), the concentration of free fatty acids and some organic acids (such as succinic acid, ribonic acid, and galactonic/gluconic acid) increased, while the concentration of pantothenic acid (vitamin B5) decreased [54]. In the study on freeze-drying, by using 1H-NMR and 13C-NMR combined with MS, the authors monitored the effect of the freeze-drying process and different storage conditions (varying temperatures and durations) on the milk metabolome over a period of 224 days. Multivariate analysis of the combined data showed that the freeze-drying process had a slight effect on the content of milk metabolites. Additionally, stable metabolome profiles were observed in freeze-dried milk samples stored in the fridge (4 °C) or in the freezer (−20 °C). In contrast, significant changes in the content of some milk metabolites were detected in freeze-dried milk samples stored at room temperature (20 °C), including decreased concentrations of riboflavin, orotic acid, and acetyl carbohydrates, increased concentrations of fatty acids, threonic acid (an oxidized product of ascorbic acid), and uridine [19]. Indeed, the authors highlighted that both vat pasteurization and freeze-drying are effective and mild preservation methods that cause only minor changes in the milk metabolome [19,54].
Innovative processing applications and their effects on milk quality have also been studied using NMR. Lemos et al. [20] focused on the application of NMR-based metabolomics to investigate the effects of hyperbaric preservation on the quality of high-pressure pasteurized (HPP) milk. The study monitored the impact of storage on the metabolic profile of HPP milk stored at an uncontrolled room temperature (~20 °C) with three different pressure levels (50, 75, and 100 MPa) and a storage period of 40 days. For comparison, two control groups were stored at either room temperature or refrigeration (4 °C) under atmospheric pressure. Visual NMR spectra evaluation revealed significant differences between milk samples subjected to hyperbaric storage and the control groups, even after just three days of storage. From the analysis of the NMR spectra, milk samples stored at room temperature showed indications of spoilage, including increased levels of lactate and butyrate; therefore, they were not included in the data analysis. Moreover, milk samples stored under hyperbaric conditions showed broad peaks in both aliphatic and aromatic regions due to an increase in soluble protein levels over time, which was mainly attributed to HPP treatment. This increase in protein content was associated with higher protein digestibility and higher nutritional value. The authors emphasized the potential of NMR-based metabolomics as a powerful tool for monitoring the effects of hyperbaric preservation of HPP milk quality and shelf life [20]. In another study, Yang et al. [21] investigated the effects of single- and double-cycled hydrostatic pressure treatment on whole and skimmed milk samples with high bacterial loads. They utilized 1H-NMR for metabolic profiling, GC-MS for the analysis of fatty acids and volatile compounds, and sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) for protein profiling. The study revealed that skimmed milk is more sensitive to high-pressure treatment compared to whole milk, indicating the possible baroprotective effect of milk fat on milk micelles. Additionally, multi-cycled hydrostatic pressure treatment was more effective in microbial inactivation compared to single-cycle treatment, resulting in better preservation of milk quality [21]. Overall, these studies conclude that NMR-based metabolomics can be used to monitor milk quality during processing and storage, making it a valuable tool for shelf-life studies of milk.

7.5. Milk Nutritional and Technological Quality

The NMR method has been successfully used for the evaluation of the nutritional and technological quality of bovine milk. Sundekilde et al. [46] investigated the relationship between the metabolic profile and technological properties of bovine milk by using NMR-based metabolomics. In the study, skimmed milk samples were classified based on their coagulation properties, which had previously been characterized as either good or poor-coagulating milk. The samples were analyzed using 1H-NMR and 13C-NMR spectroscopy combined with multivariate analysis. The authors found a positive correlation between the coagulation properties of milk and the concentrations of lactose and choline. On the other hand, they observed a negative correlation between carnitine and citrate and the coagulation properties of milk [46]. Therefore, the authors suggested that NMR-based metabolomics could be a rapid classification method to differentiate milk according to its coagulation properties. In a subsequent study, Sundekilde et al. [45] examined the correlation between milk metabolites and their protein content and how they affect the coagulation properties of milk. The study revealed that specific metabolites related to milk protein content, such as choline, creatinine, and carnitine, showed notable differences between non-coagulating and well-coagulating milk samples, indicating their potential use as indicators of milk quality in terms of coagulation properties [45].
The regulation of infant formula requires a detailed declaration of ingredients on the label, as the formulation plays a crucial role in the healthy growth of infants. An untargeted 1H-NMR approach was used by Zhao et al. [23] to determine the nutritional value and compliance of the infant formulas with regulatory requirements. The low-molecular-weight organic compounds in whole bovine milk and commercial infant formulas from different brands were investigated. The study revealed significant differences between milk and infant formulas, as well as the presence of some metabolites that were not declared on the label [23]. Moreover, the authors highlighted the importance of determining choline, creatine, and certain nucleotides/nucleosides by NMR, as these milk compounds have beneficial bioactivities for infants and hence can be used as marker compounds to validate the nutritional value of infant formula [23].

7.6. Milk Cow Health Status

The health and metabolic status of dairy cows strongly influence milk composition and, consequently, the quality of milk. Luangwilai et al. [15] conducted a study in which they characterized and differentiated 46 milk metabolites in raw milk obtained from healthy cows, as well as cows with subclinical and clinical mastitis, using untargeted 1H-NMR. Mastitis inflammation in cows, which leads to an extremely high somatic cell count, was found to be correlated with variations in the milk metabolome. This resulted in significantly increased concentrations of various free amino acids and organic acids in milk. Moreover, the authors reported that increased levels of alanine, valerate, and N-acetylglucosamine in milk can be used as new potential biomarkers for diagnosing mastitis in cows [15]. In another study, Sunds et al. [52] investigated glycosidase enzymatic activities in mastitis milk and their effects on milk composition using untargeted 1H-NMR metabolomics. The presence of certain glycosidases was associated with mastitis inflammation, which had a substantial impact on the milk metabolite profile. The enzymatic activities led to increased levels of free sugars and decreased levels of lactose. Additionally, higher concentrations of fat, protein, and free amino acids were observed, with the latter being attributed to the proteolytic activities of microorganisms due to mastitis [52]. Followingly, the mechanisms of mastitis inflammation were explained by C. Zhu et al. [16] through milk metabolomics using untargeted 1H-NMR. The study provided quantitative information on milk metabolites by comprehensively analyzing both milk and serum samples. They characterized and quantified 54 milk metabolites; a greater number compared to previous reports in the literature. The variations in metabolites between serum and milk samples were analyzed to understand the metabolic pathway of mastitis in milk [16]. In a different study, Xu et al. [53] investigated the biological pathways of negative energy balance in dairy cows by using 1H-NMR combined with LC-MS through integrated analysis. They focused on the metabolite profile of milk serum obtained from cows in the second week of lactation. Milk samples from cows experiencing negative energy balance showed various alterations in the milk metabolic profile. The authors suggested that the combination of NMR with LC-MS as a complementary technique can provide integrated and reliable information about the energy status of dairy cows and, consequently, the nutritional value of milk [53].

8. NMR Applications on Bovine Milk: Authenticity

NMR has emerged as a reliable method for verifying the authenticity of bovine milk. It has been employed as an analytical tool for assessing the authenticity of organic milk [28,29], determining the composition of reconstituted milk [18,55], detecting adulteration in milk powder [22,39], and identifying milk adulteration with melamine [56]. In addition, time-domain NMR has been applied in the assessment of milk quality and authenticity [58,59]. Table 2 offers a concise overview of the application of NMR studies in this field.

8.1. Bovine Milk Authenticity

Erich et al. studied the fat composition of milk to assess the authenticity of organic milk. In the study, 1H-NMR and 13C-NMR methods were used in combination with stable isotope-ratio mass spectrometry (IRMS) and gas chromatography (GC). Visual NMR spectra comparison combined with chemometric models revealed that organic milk exhibited stronger signals on the bis-allyl methylene groups of α-linolenic acid, the carboxyl group of butyric acid, and the long-chain fatty acids than conventional milk samples [28]. Moreover, chemometric analysis (PCA and various classification methods) applied to NMR data combined with stable-isotope data of milk protein and fat (obtained from IRMS) and α-linolenic acid content (quantified using GC) allowed to differentiate organic milk from conventional milk samples [28]. Tsiafoulis et al. [29] analyzed the lipid fraction of lyophilized organic and conventional milk samples using both targeted and untargeted NMR approaches, with a particular interest in the potential health benefits of organic milk for human consumption. From 1D and 2D-NMR spectra, significantly higher concentrations of conjugated linoleic acid (CLA), linoleic and α-linolenic acids, and overall unsaturated fatty acid content were measured in organic milk. Those molecules are associated with the better nutritional value of organic milk, providing human health benefits upon consumption [29].
Renou et al. [30] conducted a study to determine the source of milk from two distinct regions—mountainous and plain—where cows were fed with a combination of pasture and silage. Milk fatty acid composition was studied using 13C-NMR spectroscopy, whereas the milk aqueous fraction was studied by IRMS as a complementary method. Discriminant analyses considering two factors—geographic origin and feeding—were performed based on five variables: Relative proportions of monounsaturated fatty acids (MUFA), polyunsaturated fatty acids (PUFA), saturated fatty acids (SFA), 18O/16O, and deuterium/hydrogen ratios. The study revealed that the level of PUFA in mountain milk was considerably higher than in plain milk, while MUFA and SFA levels did not differ significantly between the two regions. As a result, milk from cows that grazed on pasture had a higher concentration of MUFA compared to milk from cows fed silage. In addition to the impact of feeding practices on milk fatty acid composition, the study conducted by Renou et al. [30] found that oxygen isotope ratios were more influenced by geographical location.
Furthermore, Sacco et al. [43] employed NMR in combination with IRMS to assess the authenticity of milk from Southern Italy. This innovative approach was compared to traditional spectroscopic and chromatographic techniques. The 1H-NMR spectra of milk samples from Southern Italy showed some differences when compared to milk from foreign countries, which were also confirmed from 2D-NMR molecular assignments. The 1H-NMR spectra of southern Italian milk samples displayed two peaks at 1.32 and 4.11 ppm, which represent CH- and CH3-signals of lactic acid, respectively. In contrast, these peaks were absent in foreign milk samples. Thus, foreign milk showed a higher glycerol and sugar content, as well as a different amino acid composition, compared to Southern Italian milk. These differences in milk metabolites were attributed to distinct feeding systems. Additionally, a multivariant statistical analysis using the variables from the most discriminating peaks and isotope ratios confirmed clear differentiation between milk samples from Southern Italy and those from foreign countries [43]. These studies revealed that NMR combined with IRMS provides a more precise approach for geographical origin determination and authenticity assessment of milk compared with traditional techniques [30,43]. Furthermore, Tenori et al. [44] made significant advancements in determining the geographical origin of milk using NMR-based metabolomics. Their approach successfully distinguished milk samples from different farms and brands within the same area and region. However, they emphasized the need to consider seasonal variations that can impact milk composition [44]. Overall, NMR can provide a fast and dependable way to detect the geographical origin of milk, which can be particularly relevant for EU-regulated geographically indicated milk products.

8.2. Bovine Milk Adulteration

NMR spectroscopy has been effectively used in adulteration studies to identify the presence of milk from various other sources in bovine milk. Lamanna et al. [18] conducted a study using 1H-NMR metabolomics to detect the amount of sheep milk in bovine milk as a means to detect adulteration. The study involved preparing a milk mixture of sheep milk and bovine milk at different concentrations. By analyzing the aqueous fraction of milk samples, the study identified spectral differences in metabolites such as citrate, lactate, and protein content, enabling the differentiation of pure bovine milk from pure sheep milk. Moreover, a combination of 1H-NMR metabolic profiles and multilinear regression statistical analysis allowed the determination of the relative concentration of each milk type in the milk mixture [18].
Furthermore, Li et al. [42] applied 1D and 2D-NMR experiments combined with chemometrics to distinguish between bovine, goat, and soy milk. The study identified specific metabolites that differentiate these milk types, such as higher levels of N-acetyl carbohydrates, ethanolamine, citrate, and lecithin in bovine milk compared to goat milk. On the other hand, lower levels of acetate, carnitine, and creatine were found in bovine milk. Additionally, the presence of D-sucrose was identified as a biomarker for the adulteration of bovine milk with soymilk [42]. This method proved effective in detecting bovine milk adulteration, even at low concentrations (as low as 2% (v/v)) [42]. Another study highlights the importance of N-acetyl carbohydrates as a significant biomarker for detecting bovine milk adulteration in caprine milk using 1H-NMR with multivariate data analysis (Figure 5) [55]. Moreover, the study conducted by Cui et al. [39] focused on the adulteration of ultra-high temperature processing (UHT) milk with reconstituted milk using NMR metabolite profiling. The study employed 1D and 2D-NMR experiments with chemometrics to develop a method for identifying UHT and reconstituted milk samples. Through their analysis, the authors identified L-carnitine, succinate, and acetate as biomarkers that could be used to detect the adulteration of UHT with reconstituted milk. Indeed, NMR metabolite profiling can be an effective tool for detecting possible adulteration. This is particularly relevant because UHT milk and reconstituted milk share similarities in their production processes, making it challenging to differentiate them using traditional methods [39].
A study conducted by Lachenmeier et al. [56] delved into the issue of adulteration of bovine milk-based infant formula with melamine, a compound that can have adverse health effects on infants. In order to detect melamine, the study compared the effectiveness of solution NMR and high-resolution magic angle spinning (HR-MAS). The experiments were performed on 400 MHz and 700 MHz NMR spectrometers, respectively. Results obtained from the NMR experiments were then compared with the most commonly used LC-MS/MS method. The study found that both NMR spectroscopy techniques were sensitive in detecting melamine adulteration in infant formula, thanks to the distinct singlet peak of the NH2 groups of melamine (resonating at 5.93 ppm), which do not overlap with other signals from milk metabolites. However, the reference LC-MS/MS procedure was the most sensitive, with LOD = 0.005 mg kg−1 (expressed as the amount of detected melamine per kg of infant formula), followed by HR-MAS and solution-state NMR (LOD = 0.69 and 33.3 mg kg−1, respectively) [58]. Nevertheless, the authors highlighted certain advantages of NMR and HRMAS over the LC-MS/MS procedure. These advantages include the minimal time required for sample preparation compared to the more complex preparation steps involved in LC-MS/MS and the possibility of performing non-targeted analysis.
In another study, the ability of NMR to detect adulterated bovine milk powder samples was exploited by Bergana et al. [22]. The samples were collected worldwide and artificially spiked with the most common adulterants in milk, such as melamine, dicyandiamide, urea, sucrose, maltodextrin, ammonium sulfate, soy protein isolates, and whey protein concentrate. The study found that 1H-NMR combined with conformity index analysis (i.e., a statistical approach used to assess the level of adulteration in a sample) allowed for the detection of low concentrations of all adulterants (0.005–5% w/w), especially for melamine (0.005% w/w), which exhibited a distinct and characteristic signal at 5.95–5.96 ppm in the NMR spectrum [22]. The study concluded that NMR spectroscopy is a powerful and rapid method for ensuring authenticity and detecting adulteration in bovine milk powder.
Beyond high-resolution NMR applications in milk analysis, the use of low-resolution time domain NMR (TD-NMR) has gained attention for online milk processing and quality control in the dairy industry [12,97]. TD-NMR offers practicality, portability, and relatively low instrumental costs, making it suitable for industrial settings. Unlike high-resolution NMR, TD-NMR does not require cryogenic liquids and enables rapid routine NMR analysis in a non-invasive manner. This allows for online NMR analysis directly through the milk packaging, facilitating real-time monitoring and quality control during milk production. Soyler et al. [57] have effectively utilized a low-field benchtop 1H-NMR (43 MHz) instrument to observe the enzymatic hydrolysis of lactose in milk under continuous flow mode. NMR spectra of the reaction were acquired every seven minutes over a total period of 280 min, allowing the tracking of lactose consumption by monitoring the β-lactose peak. However, the signals from the α-lactose site, which include lactose, glucose, and galactose, were concealed due to signal overlap, and the signals from β-galactose, and β-glucose were hidden by the suppressed water signals. Despite those challenges posed by signal overlap at the low magnetic field strength, the researchers were still able to monitor the enzymatic reaction by exploiting the constant anomeric ratio between the α- and β-anomers of lactose throughout the monitored period [57]. This demonstrates the potential of benchtop NMR at the industrial level for online control of residual lactose concentration in lactose-free milk products, which is particularly relevant for meeting regulatory specifications regarding lactose content in milk.
In the field of milk adulteration and quality control, 1H TD-NMR has also been utilized. Santos et al. [58] applied 1H TD-NMR combined with T2 relaxation measurements and chemometric methods to detect adulteration in milk samples. The milk samples were artificially adulterated with common adulterants (water, whey, urea, etc.) at different concentrations in the range from 5% to 50% v/v. The adulteration level in milk samples could be detected through the T2 relaxation time, as higher concentrations of adulterated samples exhibited longer relaxation times compared to control samples. This trend was consistently observed across all milk samples adulterated with various adulterants at the same concentration level [58]. In another study, Coimbra et al. [59] successfully detected formaldehyde adulteration in bovine milk using 1H TD-NMR during refrigerated storage (0 and 48 h). The T2 relaxation curves showed that changes in relaxation times were directly proportional to the concentration of added formaldehyde, with higher concentrations of formaldehyde resulting in increased T2 relaxation times. The authors explained that the addition of formaldehyde caused coagulation of casein in adulterated milk samples during storage, leading to alterations in water and fat mobility and consequently increased T2 relaxation times. Such studies highlight the effectiveness and reliability of TD-NMR in the online assessment of milk quality and authenticity.

9. Conclusions and Future Perspectives

This review focuses on the application of NMR spectroscopy in bovine milk analysis, specifically highlighting studies related to milk animal species origin, nutritional quality, geographical origin, and adulteration of bovine milk. NMR spectroscopy offers unique advantages in milk studies, including minimal to no sample preparation, rapid untargeted analysis, and the ability to identify unknown milk metabolites, which can serve as potential biomarkers for assessing the authenticity and quality of bovine milk.
The high limit of detection and the need for standardization of NMR methods are still considerable limitations compared to MS-based platforms, which are widely used in milk metabolomics. Nevertheless, ongoing advancements in NMR probe and magnet designs, magnet field strength, pulse sequences, and software and database improvements aim to address these limitations and position NMR as a complementary tool to MS platforms.
In milk research, it is recommended to adopt a multiplatform approach that leverages the strengths of each technique. A comprehensive study should employ NMR to achieve extensive metabolic coverage and high-throughput analysis. NMR can be combined with MS as a complementary method, considering its strength as a low-limit detection method. This sequential approach involves untargeted NMR analysis followed by the quantification of especially low-concentration milk metabolites that have been identified as important through NMR analysis.
Overall, NMR spectroscopy serves as a promising technique in milk studies, and with ongoing advancements, it holds great promise for further advancements in milk quality assessment and authentication research.

Author Contributions

Conceptualization, D.E., P.R. and A.C.; writing—original draft preparation, D.E. and A.C.; writing—review and editing, D.E., P.R., M.O. and A.C.; supervision, M.O. and M.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Autonomous Province of Bozen-Bolzano, Department of Innovation, Research and University within the NOI Capacity Building I and II funding frame (Decision n. 1472, 07.10.2013; Decision 864, 04.09.2018). Laimburg Research Centre is funded by the Autonomous Province of Bozen-Bolzano. The PhD-Project of E.D., entitled “Metabolic profiling of alpine food products using NMR”, was funded by the Laimburg Research Centre.

Acknowledgments

The authors thanks the Department of Innovation, Research, University and Museums of the Autonomous Province of Bozen/Bolzano for covering the Open Access publication costs.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Emwas, A.-H.; Roy, R.; McKay, R.T.; Tenori, L.; Saccenti, E.; Gowda, G.A.N.; Raftery, D.; Alahmari, F.; Jaremko, L.; Jaremko, M.; et al. NMR Spectroscopy for Metabolomics Research. Metabolites 2019, 9, 123. [Google Scholar] [CrossRef]
  2. Odeblad, E.; Westin, B. Proton Magnetic Resonance of Human Milk. Acta Radiol. 1958, 49, 389–392. [Google Scholar] [CrossRef]
  3. Chaudhary, V.; Kajla, P.; Dewan, A.; Pandiselvam, R.; Socol, C.T.; Maerescu, C.M. Spectroscopic Techniques for Authentication of Animal Origin Foods. Front. Nutr. 2022, 9, 979205. [Google Scholar] [CrossRef]
  4. Maher, A.; Rochfort, S. Applications of NMR in Dairy Research. Metabolites 2014, 4, 131–141. [Google Scholar] [CrossRef]
  5. Suh, J.H. Critical Review: Metabolomics in Dairy Science—Evaluation of Milk and Milk Product Quality. Food Res. Int. 2022, 154, 110984. [Google Scholar] [CrossRef]
  6. Sundekilde, U.; Larsen, L.; Bertram, H. NMR-Based Milk Metabolomics. Metabolites 2013, 3, 204–222. [Google Scholar] [CrossRef]
  7. Scano, P.; Cusano, E.; Caboni, P.; Consonni, R. NMR Metabolite Profiles of Dairy: A Review. Int. Dairy J. 2019, 90, 56–67. [Google Scholar] [CrossRef]
  8. Qin, C.; Liu, L.; Wang, Y.; Leng, T.; Zhu, M.; Gan, B.; Xie, J.; Yu, Q.; Chen, Y. Advancement of Omics Techniques for Chemical Profile Analysis and Authentication of Milk. Trends Food Sci. Technol. 2022, 127, 114–128. [Google Scholar] [CrossRef]
  9. Belloque, J. Application of NMR Spectroscopy to Milk and Dairy Products. Trends Food Sci. Technol. 1999, 10, 313–320. [Google Scholar] [CrossRef]
  10. Kamal, M.; Karoui, R. Analytical Methods Coupled with Chemometric Tools for Determining the Authenticity and Detecting the Adulteration of Dairy Products: A Review. Trends Food Sci. Technol. 2015, 46, 27–48. [Google Scholar] [CrossRef]
  11. Karoui, R.; Debaerdemaeker, J. A Review of the Analytical Methods Coupled with Chemometric Tools for the Determination of the Quality and Identity of Dairy Products. Food Chem. 2007, 102, 621–640. [Google Scholar] [CrossRef]
  12. Hatzakis, E. Nuclear Magnetic Resonance (NMR) Spectroscopy in Food Science: A Comprehensive Review: NMR Spectroscopy in Food Science. Compr. Rev. Food Sci. Food Saf. 2019, 18, 189–220. [Google Scholar] [CrossRef]
  13. Snytnikova, O.A.; Khlichkina, A.A.; Sagdeev, R.Z.; Tsentalovich, Y.P. Evaluation of Sample Preparation Protocols for Quantitative NMR-Based Metabolomics. Metabolomics 2019, 15. [Google Scholar] [CrossRef]
  14. Foroutan, A.; Guo, A.C.; Vazquez-Fresno, R.; Lipfert, M.; Zhang, L.; Zheng, J.; Badran, H.; Budinski, Z.; Mandal, R.; Ametaj, B.N.; et al. Chemical Composition of Commercial Cow’s Milk. J. Agric. Food Chem. 2019, 67, 4897–4914. [Google Scholar] [CrossRef]
  15. Luangwilai, M.; Duangmal, K.; Chantaprasarn, N.; Settachaimongkon, S. Comparative Metabolite Profiling of Raw Milk from Subclinical and Clinical Mastitis Cows Using 1H-NMR Combined with Chemometric Analysis. Int. J. Food Sci. Technol. 2021, 56, 493–503. [Google Scholar] [CrossRef]
  16. Zhu, C.; Tang, K.; Lu, X.; Tang, J.; Laghi, L. An Untargeted Metabolomics Investigation of Milk from Dairy Cows with Clinical Mastitis by 1H-NMR. Foods 2021, 10. [Google Scholar] [CrossRef]
  17. Lanza, I.; Lolli, V.; Segato, S.; Caligiani, A.; Contiero, B.; Lotto, A.; Galaverna, G.; Magrin, L.; Cozzi, G. Use of GC–MS and 1H NMR Low-Level Data Fusion as an Advanced and Comprehensive Metabolomic Approach to Discriminate Milk from Dairy Chains Based on Different Types of Forage. Int. Dairy J. 2021, 123, 105174. [Google Scholar] [CrossRef]
  18. Lamanna, R.; Braca, A.; Paolo, E.D.; Imparato, G. Identification of Milk Mixtures by 1H NMR Profiling. Magn. Reson. Chem. 2011, 49. [Google Scholar] [CrossRef]
  19. Zhu, D.; Kebede, B.; Chen, G.; McComb, K.; Frew, R. Impact of Freeze-Drying and Subsequent Storage on Milk Metabolites Based on 1H NMR and UHPLC-QToF/MS. Food Control 2020, 116. [Google Scholar] [CrossRef]
  20. Lemos, Á.T.; Goodfellow, B.J.; Delgadillo, I.; Saraiva, J.A. NMR Metabolic Composition Profiling of High Pressure Pasteurized Milk Preserved by Hyperbaric Storage at Room Temperature. Food Control 2022, 134, 108660. [Google Scholar] [CrossRef]
  21. Yang, S.; Liu, G.; Munk, D.M.E.; Qin, Z.; Petersen, M.A.; Cardoso, D.R.; Otte, J.; Ahrné, L. Cycled High Hydrostatic Pressure Processing of Whole and Skimmed Milk: Effects on Physicochemical Properties. Innov. Food Sci. Emerg. Technol. 2020, 63. [Google Scholar] [CrossRef]
  22. Bergana, M.M.; Adams, K.M.; Harnly, J.; Moore, J.C.; Xie, Z. Non-Targeted Detection of Milk Powder Adulteration by 1H NMR Spectroscopy and Conformity Index Analysis. J. Food Compos. Anal. 2019, 78, 49–58. [Google Scholar] [CrossRef]
  23. Zhao, Y.; Chen, H.; Feng, J.; Chen, Z.; Cai, S. 1H NMR-Based Compositional Identification of Different Powdered Infant Formulas. Food Chem. 2017, 230, 164–173. [Google Scholar] [CrossRef]
  24. Koskela, H. Chapter 1 Quantitative 2D NMR Studies. In Annual Reports on NMR Spectroscopy; Elsevier: Amsterdam, The Netherlands, 2009; Volume 66, pp. 1–31. ISBN 978-0-12-374737-2. [Google Scholar]
  25. Tang, F.; Vasas, M.; Hatzakis, E.; Spyros, A. Magnetic Resonance Applications in Food Analysis. In Annual Reports on NMR Spectroscopy; Elsevier, Amsterdam, The Netherlands, 2019; Volume 98, pp. 239–306, ISBN 978-0-08-102910-7.
  26. Andreotti, G.; Lamanna, R.; Trivellone, E.; Motta, A. 13C NMR Spectra of TAG: An Easy Way to Distinguish Milks from Different Animal Species. JAOCS J. Am. Oil Chem. Soc. 2002, 79, 123–127. [Google Scholar] [CrossRef]
  27. Phuenpong, T.; Kongboonkird, M.; Duangmal, K.; Lerdvorasap, W.; Suksawwawimon, M.; Mekboonsonglarp, W.; Nuamchit, J.; Chantaprasarn, N.; Settachaimongkon, S. Molecular Discrimination between Organic and Conventional Liquid Milk Products in Thailand Using 1H-NMR Metabolomics Approach. Trop. Anim. Sci. J. 2021, 44, 478–488. [Google Scholar] [CrossRef]
  28. Erich, S.; Schill, S.; Annweiler, E.; Waiblinger, H.U.; Kuballa, T.; Lachenmeier, D.W.; Monakhova, Y.B. Combined Chemometric Analysis of 1H NMR, 13C NMR and Stable Isotope Data to Differentiate Organic and Conventional Milk. Food Chem. 2015, 188, 1–7. [Google Scholar] [CrossRef] [PubMed]
  29. Tsiafoulis, C.; Papaemmanouil, C.; Alivertis, D.; Tzamaloukas, O.; Miltiadou, D.; Balayssac, S.; Malet-Martino, M.; Gerothanassis, I. NMR-Based Μetabolomics of the Lipid Fraction of Organic and Conventional Bovine Milk. Molecules 2019, 24, 1067. [Google Scholar] [CrossRef]
  30. Renou, J.P.; Deponge, C.; Gachon, P.; Bonnefoy, J.C.; Coulon, J.B.; Garel, J.P.; Vérité, R.; Ritz, P. Characterization of Animal Products According to Geographic Origin and Feeding Diet Using Nuclear Magnetic Resonance and Isotope Ratio Mass Spectrometry: Cow Milk. Food Chem. 2004, 85, 63–66. [Google Scholar] [CrossRef]
  31. Bruschetta, G.; Notti, A.; Lando, G.; Ferlazzo, A. A Promising 31P NMR-Multivariate Analysis Approach for the Identification of Milk Phosphorylated Metabolites and for Rapid Authentication of Milk Samples. Biochem. Biophys. Rep. 2021, 27. [Google Scholar] [CrossRef]
  32. Belloque, J. High-Resolution NMR of Milk and Milk Proteins. In Modern Magnetic Resonance; Webb, G.A., Ed.; Springer: Dordrecht, The Netherlands, 2006; pp. 1631–1635. ISBN 978-1-4020-3894-5. [Google Scholar]
  33. Garcia, C.; Lutz, N.W.; Confort-Gouny, S.; Cozzone, P.J.; Armand, M.; Bernard, M. Phospholipid Fingerprints of Milk from Different Mammalians Determined by 31P NMR: Towards Specific Interest in Human Health. Food Chem. 2012, 135, 1777–1783. [Google Scholar] [CrossRef]
  34. Wei, W.; Li, D.; Jiang, C.; Zhang, X.; Zhang, X.; Jin, Q.; Zhang, X.; Wang, X. Phospholipid Composition and Fat Globule Structure II: Comparison of Mammalian Milk from Five Different Species. Food Chem. 2022, 388. [Google Scholar] [CrossRef] [PubMed]
  35. Boiani, M.; Fenelon, M.; FitzGerald, R.J.; Kelly, P.M. Use of 31P NMR and FTIR to Investigate Key Milk Mineral Equilibria and Their Interactions with Micellar Casein during Heat Treatment. Int. Dairy J. 2018, 81, 12–18. [Google Scholar] [CrossRef]
  36. Hu, F.; Furihata, K.; Ito-Ishida, M.; Kaminogawa, S.; Tanokura, M. Nondestructive Observation of Bovine Milk by NMR Spectroscopy: Analysis of Existing States of Compounds and Detection of New Compounds. J. Agric. Food Chem. 2004, 52, 4969–4974. [Google Scholar] [CrossRef]
  37. Andreotti, G.; Trivellone, E.; Motta, A. Characterization of Buffalo Milk by 31P-Nuclear Magnetic Resonance Spectroscopy. J. Food Compos. Anal. 2006, 19, 843–849. [Google Scholar] [CrossRef]
  38. Hu, F.; Furihata, K.; Kato, Y.; Tanokura, M. Nondestructive Quantification of Organic Compounds in Whole Milk without Pretreatment by Two-Dimensional NMR Spectroscopy. J. Agric. Food Chem. 2007, 55, 4307–4311. [Google Scholar] [CrossRef] [PubMed]
  39. Cui, J.; Zhu, D.; Su, M.; Tan, D.; Zhang, X.; Jia, M.; Chen, G. The Combined Use of 1H and 2D NMR-Based Metabolomics and Chemometrics for Non-Targeted Screening of Biomarkers and Identification of Reconstituted Milk. J. Sci. Food Agric. 2019, 99, 6455–6461. [Google Scholar] [CrossRef]
  40. Wishart, D.S.; Knox, C.; Guo, A.C.; Eisner, R.; Young, N.; Gautam, B.; Hau, D.D.; Psychogios, N.; Dong, E.; Bouatra, S.; et al. HMDB: A Knowledgebase for the Human Metabolome. Nucleic Acids Res. 2009, 37, D603–D610. [Google Scholar] [CrossRef]
  41. Atta-ur-Rahman; Choudhary, M.; Atia-tul-Wahab. Solving Problems with NMR Spectroscopy; Elsevier: Amsterdam, The Netherlands, 1996; ISBN 978-0-12-066320-0. [Google Scholar]
  42. Li, Q.; Yu, Z.; Zhu, D.; Meng, X.; Pang, X.; Liu, Y.; Frew, R.; Chen, H.; Chen, G. The Application of NMR-Based Milk Metabolite Analysis in Milk Authenticity Identification. J. Sci. Food Agric. 2017, 97, 2875–2882. [Google Scholar] [CrossRef] [PubMed]
  43. Sacco, D.; Brescia, M.A.; Sgaramella, A.; Casiello, G.; Buccolieri, A.; Ogrinc, N.; Sacco, A. Discrimination between Southern Italy and Foreign Milk Samples Using Spectroscopic and Analytical Data. Food Chem. 2009, 114, 1559–1563. [Google Scholar] [CrossRef]
  44. Tenori, L.; Santucci, C.; Meoni, G.; Morrocchi, V.; Matteucci, G.; Luchinat, C. NMR Metabolomic Fingerprinting Distinguishes Milk from Different Farms. Food Res. Int. 2018, 113, 131–139. [Google Scholar] [CrossRef]
  45. Sundekilde, U.K.; Gustavsson, F.; Poulsen, N.A.; Glantz, M.; Paulsson, M.; Larsen, L.B.; Bertram, H.C. Association between the Bovine Milk Metabolome and Rennet-Induced Coagulation Properties of Milk. J. Dairy Sci. 2014, 97, 6076–6084. [Google Scholar] [CrossRef]
  46. Sundekilde, U.K.; Frederiksen, P.D.; Clausen, M.R.; Larsen, L.B.; Bertram, H.C. Relationship between the Metabolite Profile and Technological Properties of Bovine Milk from Two Dairy Breeds Elucidated by NMR-Based Metabolomics. J. Agric. Food Chem. 2011, 59, 7360–7367. [Google Scholar] [CrossRef]
  47. Consonni, R.; Cagliani, L.R. The Potentiality of NMR-Based Metabolomics in Food Science and Food Authentication Assessment. Magn Reson. Chem. 2019, 57, 558–578. [Google Scholar] [CrossRef]
  48. Simmler, C.; Napolitano, J.G.; McAlpine, J.B.; Chen, S.-N.; Pauli, G.F. Universal Quantitative NMR Analysis of Complex Natural Samples. Curr. Opin. Biotechnol. 2014, 25, 51–59. [Google Scholar] [CrossRef]
  49. Rocchetti, G.; O’Callaghan, T.F. Application of Metabolomics to Assess Milk Quality and Traceability. Curr. Opin. Food Sci. 2021, 40, 168–178. [Google Scholar] [CrossRef]
  50. 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. [Google Scholar] [CrossRef]
  51. O’Callaghan, T.F.; Vázquez-Fresno, R.; Serra-Cayuela, A.; Dong, E.; Mandal, R.; Hennessy, D.; McAuliffe, S.; Dillon, P.; Wishart, D.S.; Stanton, C.; et al. Pasture Feeding Changes the Bovine Rumen and Milk Metabolome. Metabolites 2018, 8. [Google Scholar] [CrossRef] [PubMed]
  52. Sunds, A.V.; Sundekilde, U.K.; Poulsen, N.A.; Larsen, L.B. Association between Udder Inflammation and Glycosidase Activities and Free Sugar Levels in Bovine Milk. Int. Dairy J. 2021, 120. [Google Scholar] [CrossRef]
  53. Xu, W.; Knegsel, A.V.; Saccenti, E.; Hoeij, R.V.; Kemp, B.; Vervoort, J. Metabolomics of Milk Reflects a Negative Energy Balance in Cows. J. Proteome Res. 2020, 19, 2942–2949. [Google Scholar] [CrossRef]
  54. Zhu, D.; Kebede, B.; Chen, G.; McComb, K.; Frew, R. Effects of the Vat Pasteurization Process and Refrigerated Storage on the Bovine Milk Metabolome. J. Dairy Sci. 2020, 103, 2077–2088. [Google Scholar] [CrossRef]
  55. Rysova, L.; Legarova, V.; Pacakova, Z.; Hanus, O.; Nemeckova, I.; Klimesova, M.; Havlik, J. Detection of Bovine Milk Adulteration in Caprine Milk with N-Acetyl Carbohydrate Biomarkers by Using 1H Nuclear Magnetic Resonance Spectroscopy. J. Dairy Sci. 2021, 104, 9583–9595. [Google Scholar] [CrossRef]
  56. Lachenmeier, D.W.; Eberhard, H.; Fang, F.; Birk, S.; Peter, D.; Constanze, S.; Manfred, S. NMR-Spectroscopy for Nontargeted Screening and Simultaneous Quantification of Health-Relevant Compounds in Foods: The Example of Melamine. J. Agric. Food Chem. 2009, 57, 7194–7199. [Google Scholar] [CrossRef]
  57. Soyler, A.; Cikrikci, S.; Cavdaroglu, C.; Bouillaud, D.; Farjon, J.; Giraudeau, P.; Oztop, M.H. Multi-Scale Benchtop 1H NMR Spectroscopy for Milk Analysis. LWT 2021, 139. [Google Scholar] [CrossRef]
  58. Santos, P.M.; Pereira-Filho, E.R.; Colnago, L.A. Detection and Quantification of Milk Adulteration Using Time Domain Nuclear Magnetic Resonance (TD-NMR). Microchem. J. 2016, 124, 15–19. [Google Scholar] [CrossRef]
  59. Coimbra, P.T.; Bathazar, C.F.; Guimarães, J.T.; Coutinho, N.M.; Pimentel, T.C.; Neto, R.P.C.; Esmerino, E.A.; Freitas, M.Q.; Silva, M.C.; Tavares, M.I.B.; et al. Detection of Formaldehyde in Raw Milk by Time Domain Nuclear Magnetic Resonance and Chemometrics. Food Control 2020, 110, 107006. [Google Scholar] [CrossRef]
  60. Yamazaki, T.; Nakamura, S.; Saito, T. Optimization of Sample Preparation for Accurate Results in Quantitative NMR Spectroscopy. Metrologia 2017, 54, 224–228. [Google Scholar] [CrossRef]
  61. Nagana Gowda, G.A.; Raftery, D. NMR-Based Metabolomics. In Cancer Metabolomics; Hu, S., Ed.; Advances in Experimental Medicine and Biology; Springer International Publishing: Cham, Switzerland, 2021; Volume 1280, pp. 19–37. ISBN 978-3-030-51651-2. [Google Scholar]
  62. Yanibada, B.E.E.; Boudra, H.; Debrauwer, L.; Martin, E.; Morgavi, D.P.; Canlet, E.; Boudra, H. Evaluation of Sample Preparation Methods for NMR-Based Metabolomics of Cow Milk. Heliyon 2018, 4, e00856. [Google Scholar] [CrossRef]
  63. Dais, P.; Hatzakis, E. Quality Assessment and Authentication of Virgin Olive Oil by NMR Spectroscopy: A Critical Review. Anal. Chim. Acta 2013, 765, 1–27. [Google Scholar] [CrossRef]
  64. Zhu, D.; Kebede, B.; McComb, K.; Hayman, A.; Chen, G.; Frew, R. Milk Biomarkers in Relation to Inherent and External Factors Based on Metabolomics. Trends Food Sci. Technol. 2021, 109, 51–64. [Google Scholar] [CrossRef]
  65. Hanson, J.E. NMR Spectroscopy in Nondeuterated Solvents (No-D NMR): Applications in the Undergraduate Organic Laboratory. In ACS Symposium Series; Soulsby, D., Anna, L.J., Wallner, A.S., Eds.; American Chemical Society: Washington, DC, USA, 2013; ISBN 978-0-8412-2794-1. [Google Scholar]
  66. Cao, R.; Liu, X.; Liu, Y.; Zhai, X.; Cao, T.; Wang, A.; Qiu, J. Applications of Nuclear Magnetic Resonance Spectroscopy to the Evaluation of Complex Food Constituents. Food Chem. 2021, 342, 128258. [Google Scholar] [CrossRef]
  67. Bharti, S.K.; Roy, R. Quantitative 1H NMR Spectroscopy. TrAC Trends Anal. Chem. 2012, 35, 5–26. [Google Scholar] [CrossRef]
  68. Emwas, A.-H.; Saccenti, E.; Gao, X.; McKay, R.T.; Dos Santos, V.A.P.M.; Roy, R.; Wishart, D.S. Recommended Strategies for Spectral Processing and Post-Processing of 1D 1H-NMR Data of Biofluids with a Particular Focus on Urine. Metabolomics 2018, 14, 31. [Google Scholar] [CrossRef] [PubMed]
  69. Eads, T.M.; Bryant, R.G. High-Resolution Proton NMR Spectroscopy of Milk, Orange Juice, and Apple Juice with Efficient Suppression of the Water Peak. J. Agric. Food Chem. 1986, 34, 834–837. [Google Scholar] [CrossRef]
  70. Buda, S.; Nawój, M.; Mlynarski, J. Recent Advances in NMR Studies of Carbohydrates. In Annual Reports on NMR Spectroscopy; Elsevier: Amsterdam, The Netherlands, 2016; Volume 89, pp. 185–223. ISBN 978-0-12-804712-5. [Google Scholar]
  71. Piotto, M.; Saudek, V.; Sklen, V. Gradient-Tailored Excitation for Single-Quantum NMR Spectroscopy of Aqueous Solutions. J. Biomol. NMR 1992, 2, 661–665. [Google Scholar] [CrossRef] [PubMed]
  72. Beckonert, O.; Keun, H.C.; Ebbels, T.M.D.; Bundy, J.; Holmes, E.; Lindon, J.C.; Nicholson, J.K. Metabolic Profiling, Metabolomic and Metabonomic Procedures for NMR Spectroscopy of Urine, Plasma, Serum and Tissue Extracts. Nat. Protoc. 2007, 2, 2692–2703. [Google Scholar] [CrossRef] [PubMed]
  73. Mckay, R.T. How the 1D-NOESY Suppresses Solvent Signal in Metabonomics NMR Spectroscopy: An Examination of the Pulse Sequence Components and Evolution. Concepts Magn. Reson. 2011, 38A, 197–220. [Google Scholar] [CrossRef]
  74. Laref, N. Efficiency of different bioinformatics tools in metabolite profiling of whole cow’s milk using synthetically water-removed 1h NMR spectra: A comparative study. J. Microbiol. Biotechnol. Food Sci. 2020, 10, 241–244. [Google Scholar] [CrossRef]
  75. Niklasson, M.; Otten, R.; Ahlner, A.; Andresen, C.; Schlagnitweit, J.; Petzold, K.; Lundström, P. Comprehensive Analysis of NMR Data Using Advanced Line Shape Fitting. J. Biomol. NMR 2017, 69, 93–99. [Google Scholar] [CrossRef]
  76. Izquierdo-García, J.L.; Villa, P.; Kyriazis, A.; Del Puerto-Nevado, L.; Pérez-Rial, S.; Rodriguez, I.; Hernandez, N.; Ruiz-Cabello, J. Descriptive Review of Current NMR-Based Metabolomic Data Analysis Packages. Prog. Nucl. Magn. Reson. Spectrosc. 2011, 59, 263–270. [Google Scholar] [CrossRef]
  77. Kohl, S.M.; Klein, M.S.; Hochrein, J.; Oefner, P.J.; Spang, R.; Gronwald, W. State-of-the Art Data Normalization Methods Improve NMR-Based Metabolomic Analysis. Metabolomics 2012, 8, 146–160. [Google Scholar] [CrossRef]
  78. Guitton, Y.; Tremblay-Franco, M.; Le Corguillé, G.; Martin, J.-F.; Pétéra, M.; Roger-Mele, P.; Delabrière, A.; Goulitquer, S.; Monsoor, M.; Duperier, C.; et al. Create, Run, Share, Publish, and Reference Your LC–MS, FIA–MS, GC–MS, and NMR Data Analysis Workflows with the Workflow4Metabolomics 3.0 Galaxy Online Infrastructure for Metabolomics. Int. J. Biochem. Cell Biol. 2017, 93, 89–101. [Google Scholar] [CrossRef] [PubMed]
  79. Xi, Y.; Rocke, D.M. Baseline Correction for NMR Spectroscopic Metabolomics Data Analysis. BMC Bioinform. 2008, 9, 324. [Google Scholar] [CrossRef] [PubMed]
  80. Tomasi, G.; Van Den Berg, F.; Andersson, C. Correlation Optimized Warping and Dynamic Time Warping as Preprocessing Methods for Chromatographic Data. J. Chemom. 2004, 18, 231–241. [Google Scholar] [CrossRef]
  81. Steimers, E.; Sawall, M.; Behrens, R.; Meinhardt, D.; Simoneau, J.; Münnemann, K.; Neymeyr, K.; Von Harbou, E. Application of a New Method for Simultaneous Phase and Baseline Correction of NMR Signals (SINC). Magn. Reson. Chem. 2020, 58, 260–270. [Google Scholar] [CrossRef] [PubMed]
  82. Klein, M.S.; Buttchereit, N.; Miemczyk, S.P.; Immervoll, A.K.; Louis, C.; Wiedemann, S.; Junge, W.; Thaller, G.; Oefner, P.J.; Gronwald, W. NMR Metabolomic Analysis of Dairy Cows Reveals Milk Glycerophosphocholine to Phosphocholine Ratio as Prognostic Biomarker for Risk of Ketosis. J. Proteome Res. 2012, 11, 1373–1381. [Google Scholar] [CrossRef]
  83. Li, X.; Hu, K. Quantitative NMR Studies of Multiple Compound Mixtures. In Annual Reports on NMR Spectroscopy; Elsevier: Amsterdam, The Netherlands, 2017; Volume 90, pp. 85–143. ISBN 978-0-12-812009-5. [Google Scholar]
  84. Savorani, F.; Tomasi, G.; Engelsen, S.B. Icoshift: A Versatile Tool for the Rapid Alignment of 1D NMR Spectra. J. Magn. Reson. 2010, 202, 190–202. [Google Scholar] [CrossRef]
  85. Mannina, L.; Sobolev, A.P.; Aru, V.; Bellomaria, A.; Bertocchi, F.; Botta, B.; Cagliani, L.R.; Caligiani, A.; Capozzi, F.; Çela, D.; et al. NMR Methodologies in Food Analysis. In Analytical Chemistry—Developments, Applications, and Challenges in Food Analysis; Nova Science Publishers Inc.: Hauppauge, NY, USA, 2017. [Google Scholar]
  86. Nielsen, N.-P.V.; Carstensen, J.M.; Smedsgaard, J. Aligning of Single and Multiple Wavelength Chromatographic Profiles for Chemometric Data Analysis Using Correlation Optimised Warping. J. Chromatogr. A 1998, 805, 17–35. [Google Scholar] [CrossRef]
  87. Giraudeau, P.; Tea, I.; Remaud, G.S.; Akoka, S. Reference and Normalization Methods: Essential Tools for the Intercomparison of NMR Spectra. J. Pharm. Biomed. Anal. 2014, 93, 3–16. [Google Scholar] [CrossRef]
  88. Brennan, L. NMR-Based Metabolomics: From Sample Preparation to Applications in Nutrition Research. Prog. Nucl. Magn. Reson. Spectrosc. 2014, 83, 42–49. [Google Scholar] [CrossRef]
  89. van den Berg, R.A.; Hoefsloot, H.C.J.; Westerhuis, J.A.; Smilde, A.K.; van der Werf, M.J. Centering, Scaling, and Transformations: Improving the Biological Information Content of Metabolomics Data. BMC Genom. 2006, 7, 142. [Google Scholar] [CrossRef]
  90. Liu, X.; Fang, Y.; Ma, H.; Zhang, N.; Li, C. Performance Comparison of Three Scaling Algorithms in NMR-Based Metabolomics Analysis. Open Life Sci. 2023, 18, 20220556. [Google Scholar] [CrossRef] [PubMed]
  91. Ulrich, E.L.; Akutsu, H.; Doreleijers, J.F.; Harano, Y.; Ioannidis, Y.E.; Lin, J.; Livny, M.; Mading, S.; Maziuk, D.; Miller, Z.; et al. BioMagResBank. Nucleic Acids Res. 2008, 36, D402–D408. [Google Scholar] [CrossRef] [PubMed]
  92. Wishart, D.S. NMR Metabolomics: A Look Ahead. J. Magn. Reson. 2019, 306, 155–161. [Google Scholar] [CrossRef] [PubMed]
  93. Cui, Q.; Lewis, I.A.; Hegeman, A.D.; Anderson, M.E.; Li, J.; Schulte, C.F.; Westler, W.M.; Eghbalnia, H.R.; Sussman, M.R.; Markley, J.L. Metabolite Identification via the Madison Metabolomics Consortium Database. Nat. Biotechnol. 2008, 26, 162–164. [Google Scholar] [CrossRef] [PubMed]
  94. Kuhn, S.; Schlörer, N.E. Facilitating Quality Control for Spectra Assignments of Small Organic Molecules: Nmrshiftdb2--a Free in-House NMR Database with Integrated LIMS for Academic Service Laboratories. Magn. Reson. Chem. 2015, 53, 582–589. [Google Scholar] [CrossRef]
  95. O’Callaghan, T.F.; Hennessy, D.; McAuliffe, S.; Kilcawley, K.N.; O’Donovan, M.; Dillon, P.; Ross, R.P.; Stanton, C. Effect of Pasture versus Indoor Feeding Systems on Raw Milk Composition and Quality over an Entire Lactation. J. Dairy Sci. 2016, 99, 9424–9440. [Google Scholar] [CrossRef]
  96. Le Gall, G. Sample Collection and Preparation of Biofluids and Extracts for NMR Spectroscopy. Methods Mol. Biol. 2015, 1277, 15–28. [Google Scholar] [CrossRef]
  97. Van Duynhoven, J.; Voda, A.; Witek, M.; Van As, H. Time-Domain NMR Applied to Food Products. In Annual Reports on NMR Spectroscopy; Elsevier: Amsterdam, The Netherlands, 2010; Volume 69, pp. 145–197. ISBN 978-0-12-381355-8. [Google Scholar]
Figure 1. Workflow of the review.
Figure 1. Workflow of the review.
Foods 12 03240 g001
Figure 2. A representative 700 MHz 1H-NMR spectrum of commercial milk. The milk sample was dissolved in 10% D2O with 1 mM DSS as internal standard. The spectrum represents a comprehensive characterization of milk metabolites; it is dominated by the signals of carbohydrates (primarily lactose, glucose, and galactose; upper inset), organic acids (citrate), amine-containing compounds (creatine and urea; lower inset), orotic acid, and glycerophosphocholine. Additionally, amino acids such as tryptophan, tyrosine, histidine, isoleucine, leucine, and asparagine are also presented. Reprinted from [14] with permission from American Chemical Society.
Figure 2. A representative 700 MHz 1H-NMR spectrum of commercial milk. The milk sample was dissolved in 10% D2O with 1 mM DSS as internal standard. The spectrum represents a comprehensive characterization of milk metabolites; it is dominated by the signals of carbohydrates (primarily lactose, glucose, and galactose; upper inset), organic acids (citrate), amine-containing compounds (creatine and urea; lower inset), orotic acid, and glycerophosphocholine. Additionally, amino acids such as tryptophan, tyrosine, histidine, isoleucine, leucine, and asparagine are also presented. Reprinted from [14] with permission from American Chemical Society.
Foods 12 03240 g002
Figure 3. (a) 500 MHz 13C-NMR spectrum of whole milk in 10% D2O. The spectrum is dominated by the signals of lactose and acyl chains of lipids. Adapted from [36] with permission from American Chemical Society. (b) 400 MHz (operated at 161.97 MHz) 31P-NMR spectrum of milk ultrafiltrate obtained at pH 9.4 in triethylamine/dimethyl-formamide/guanidinium hydrochloride with 10% D2O and phosphocreatine as internal chemical shift reference. The spectrum represents the phosphorous-containing molecules in milk, with inorganic phosphate being the most abundantly present. Adapted from [37] with permission from Elsevier.
Figure 3. (a) 500 MHz 13C-NMR spectrum of whole milk in 10% D2O. The spectrum is dominated by the signals of lactose and acyl chains of lipids. Adapted from [36] with permission from American Chemical Society. (b) 400 MHz (operated at 161.97 MHz) 31P-NMR spectrum of milk ultrafiltrate obtained at pH 9.4 in triethylamine/dimethyl-formamide/guanidinium hydrochloride with 10% D2O and phosphocreatine as internal chemical shift reference. The spectrum represents the phosphorous-containing molecules in milk, with inorganic phosphate being the most abundantly present. Adapted from [37] with permission from Elsevier.
Foods 12 03240 g003
Figure 4. 500 MHz 2D 1H-13C HSQC-NMR spectrum of whole milk in 10% D2O. The spectrum shows resonances from the main constituents of milk (mainly D-lactose, acyl chains of fatty acids, and glycerol backbone of fats) Adapted from [36] with permission from American Chemical Society.
Figure 4. 500 MHz 2D 1H-13C HSQC-NMR spectrum of whole milk in 10% D2O. The spectrum shows resonances from the main constituents of milk (mainly D-lactose, acyl chains of fatty acids, and glycerol backbone of fats) Adapted from [36] with permission from American Chemical Society.
Foods 12 03240 g004
Figure 5. Example of a multivariate analysis procedure applied in NMR milk authentication studies. (a) PCA score plot; (b) OPLS-DA score plot; and (c) OPLS-DA corresponding color/coded coefficient plot of median 1H-NMR spectra of pure bovine (n = 23) and caprine milk (n = 23) represent the contribution of the metabolites toward class discrimination between bovine milk and caprine milk. The color of the bar indicates the correlation loading of the bin p(corr) value. The size of the bar was multiplied by standard deviation. Adapted from [55] with permission from American Chemical Society.
Figure 5. Example of a multivariate analysis procedure applied in NMR milk authentication studies. (a) PCA score plot; (b) OPLS-DA score plot; and (c) OPLS-DA corresponding color/coded coefficient plot of median 1H-NMR spectra of pure bovine (n = 23) and caprine milk (n = 23) represent the contribution of the metabolites toward class discrimination between bovine milk and caprine milk. The color of the bar indicates the correlation loading of the bin p(corr) value. The size of the bar was multiplied by standard deviation. Adapted from [55] with permission from American Chemical Society.
Foods 12 03240 g005
Table 1. Advantages and disadvantages of NMR spectroscopy.
Table 1. Advantages and disadvantages of NMR spectroscopy.
1D 1H-NMR1D 13C-NMR1D 31P-NMR2D-NMR
Advantages DisadvantagesAdvantages DisadvantagesAdvantagesDisadvantagesAdvantagesDisadvantages
High
sensitivity
Narrow
spectral
window
Relatively broader spectral windowLow sensitivityHigh sensitivity
and broad spectral window
Higher sensitivity than 1D 13C experiments
Short relaxation time Low resolutionHigh
resolution
Long relaxation timeHigh resolution Long relaxation timehigh resolution,
more detailed metabolite
assignment,
Short
experimental time
Overcrowded spectra Promising in milk fatty acid analysisLong
experimental time
Promising in milk phospholipid analysisLong
experimental time
Identification of unknown metabolitesLong
experimental time
Suitable for metabolic
fingerprinting
Suitable for metabolic
profiling
Suitable for
metabolic
profiling
Suitable for structural
elucidation
Usually not used for
metabolic
profiling
Table 2. Summary of main characteristics of NMR studies on bovine milk 1.
Table 2. Summary of main characteristics of NMR studies on bovine milk 1.
Application FieldSample
Preparation
pHNMR ExperimentData Pre-ProcessingOther AnalysisMVAReference
Chemical
composition
characterization
Whole and
extract
na1H, 13C-NMR, HSQC, HMBC __[36]
Extract7.01H-NMR_LC-RMS,
LC-MS/MS
ICP-MS
_[14]
Feeding regimeExtractna1H-NMRIntegrationFT-IRHCA, PLS-DA[51]
Extract7.01H-NMR, TOCSYBinning, integrationGC-MSCDA[17]
Milk quality and
animal health
Extract6.0.1H-NMRBinningFT-IRPCA, PLS-DA[15]
Extractna1H-NMRAlignment,
binning
Flow cytometry, MIR, microbiological analysisPCA[52]
Extract7.01H-NMRAlignment, binning, integrationLC-MS_[53]
Extract7.01H-NMRintegration,
normalization
_rPCA[16]
Technological
process and shelf life
Extract of lyophilized milkna1H-NMR_UPLC-QToF/MS, GC-MSPCA, PLS-DA, PLSR[54]
extract of lyophilized milk 7.41H, 13C-NMRIntegrationUPLC-QToF/MSPCA, PLS-DA, PLSR[19]
Extract 7.01H-NMRNormalization,
alignment
_PCA, PLS-DA[20]
Wholena1H, 13C-NMR_FT-IR, GC-MS, SDS-PAGE, microbiological analysis_[21]
Milk technological propertiesExtractna1H-NMR
Alignment, binningFT-IR, flow cytometry, rheometryPCA, PLS, OPLS-DA[45]
Skimmed milkna1H, 13CNMR, 1H-13C HSQC Alignment, binningFT-IR, flow cytometry, rheometryPCA[46]
Infant formula quality assessmentExtract of infant
formula
na1H-NMRIntegration, normalization, binning_PCA, PLS-DA, OPLS-DA[23]
Animal species originExtractna13C-NMRIntegration,
normalization
_Fuzzy logic analysis[26]
Extractna31P-NMR___[33]
Extractna31P-NMRIntegration, normalization_PCA, UNEQ, SIMCA, K-NN[31]
Extractna31P-NMR_GC, LC-MS_[34]
Geographical originExtractna13C-NMR_IRMSDA, SDA[30]
extract of lyophilized milk na1H-NMR, COSY, TOCSY_IRMS, HPIC, ICP-AES, FT-IRPCA, DA[43]
Extract7.41H-NMR,
1H-1H COSY, 1H-13C HSQC
Binning, integration, normalization_PCA, PLS-CA [44]
Organic milk
authentication
Extractna1H, 13C-NMR
Alignment, binning, integrationIRMSPCA, LDA, FDA, PLSDA[28]
Extract of lyophilized milk na1H-NMR, 1D, 2D TOCSY, 1H-13C HSQC, HMBCBinning, integration, normalization_PCA, PLS-DA[29]
Bovine milk
adulteration
Extractna1H-NMR,
1H-13C HSQC, HMBC
Binning_PCA, LDA, ANN[18]
Extractna1H, 13C-NMR, 1H-13C HSQC, HMBCIntegration, normalization_PCA[42]
Extract of reconstituted milk 7.41H-NMR,
1H-13C HSQC, HMBC
Integration, normalization_PCA, PLS-DA[39]
Extract7.41H-NMRBinning, integration, normalizationIR, flow cytometry, microbiological analysisPCA, OPLS-DA[55]
Infant formula fraud with
melamine
Infant
powder
7.01H-NMR,
1H HRMAS NMR
_LC/MS/MS_[56]
Milk powder
adulteration
Extract of milk
powder
na1H-NMRBinning, normalization__[22]
Online quality and adulteration
control
Whole and skimmednaBenchtop
1H-NMR
Integration, normalization_ANN[57]
Whole milkna1H TD-NMR__PCA, SIMCA, k-NN, PLSR[58]
Whole milkna1H TD-NMR_ColorimeterPCA, PLS, SIMCA[59]
1 Extract represents the sample extraction process through sample preparation. Abbreviations are: HSQC, heteronuclear single quantum coherence spectroscopy; HMBC, heteronuclear multiple bond correlation spectroscopy; TOCSY, total correlation spectroscopy; COSY, correlation spectroscopy; HRMAS, high-resolution magic angle spinning; MVA, multivariate statistical data analysis; PCA, principal component analysis; DA, discriminant analysis; UNEQ, unequal class modeling; SIMCA, soft independent modelling of class analogy; K-NN, kernel nearest neighbor; HCA, hierarchical cluster analysis; PLS-DA, partial least squares DA; CDA, canonical DA; PLSR, PLS regression; OPLS-DA, orthogonal PLS-DA; ANN, artificial neural network; rPCA, robust PCA; LDA, linear DA; FDA, factorial DA; SDA, stepwise DA; LC-HRMS, liquid chromatography high-resolution mass spectrometry; LC-MS, liquid chromatography mass spectrometry; ICP-MS, inductively coupled plasma mass spectrometry; GC, gas chromatography; FT-IR, Fourier transform infrared spectroscopy; UPLC-QToF/MS, ultraperformance liquid chromatography-quadrupole time-of-flight MS; SDS-PAGE, sodium dodecyl sulphate-polyacrylamide gel electrophoresis; MIR, mid-infrared spectroscopy; IRMS, isotope ratio MS; HPIC, high-performance ion chromatography; ICP-AES, inductively coupled plasma atomic emission spectroscopy; na, not adjusted.
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

Eltemur, D.; Robatscher, P.; Oberhuber, M.; Scampicchio, M.; Ceccon, A. Applications of Solution NMR Spectroscopy in Quality Assessment and Authentication of Bovine Milk. Foods 2023, 12, 3240. https://doi.org/10.3390/foods12173240

AMA Style

Eltemur D, Robatscher P, Oberhuber M, Scampicchio M, Ceccon A. Applications of Solution NMR Spectroscopy in Quality Assessment and Authentication of Bovine Milk. Foods. 2023; 12(17):3240. https://doi.org/10.3390/foods12173240

Chicago/Turabian Style

Eltemur, Dilek, Peter Robatscher, Michael Oberhuber, Matteo Scampicchio, and Alberto Ceccon. 2023. "Applications of Solution NMR Spectroscopy in Quality Assessment and Authentication of Bovine Milk" Foods 12, no. 17: 3240. https://doi.org/10.3390/foods12173240

APA Style

Eltemur, D., Robatscher, P., Oberhuber, M., Scampicchio, M., & Ceccon, A. (2023). Applications of Solution NMR Spectroscopy in Quality Assessment and Authentication of Bovine Milk. Foods, 12(17), 3240. https://doi.org/10.3390/foods12173240

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