Next Article in Journal
New Halogen-Containing Drugs Approved by FDA in 2021: An Overview on Their Syntheses and Pharmaceutical Use
Next Article in Special Issue
Metabonomics Analysis of Stem Extracts from Dalbergia sissoo
Previous Article in Journal
Understanding the Biosynthesis of Paxisterol in Lichen-Derived Penicillium aurantiacobrunneum for Production of Fluorinated Derivatives
Previous Article in Special Issue
Use of Gas Chromatography-Mass Spectrometry Techniques (GC-MS, GC-MS/MS and GC-QTOF) for the Characterization of Photooxidation and Autoxidation Products of Lipids of Autotrophic Organisms in Environmental Samples
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of a Nuclear Magnetic Resonance Method and a Near Infrared Calibration Model for the Rapid Determination of Lipid Content in the Field Pea (Pisum sativum)

1
Bioresource Engineering Department, McGill University, Macdonald Campus, Ste-Anne-de-Bellevue, QC H9X 3V9, Canada
2
Drug Discovery Platform, Research Institute McGill University Health Centre, Montreal, QC H4A 3J1, Canada
3
Plant Science Department, McGill University, Macdonald Campus, Ste-Anne-de-Bellevue, QC H9X 3V9, Canada
*
Author to whom correspondence should be addressed.
Submission received: 30 January 2022 / Revised: 23 February 2022 / Accepted: 28 February 2022 / Published: 2 March 2022

Abstract

:
Pisum sativum is a leguminous crop suitable for cultivation worldwide. It is used as a forage or dried seed supplement in animal feed and, more recently, as a potential non-traditional oilseed. This study aimed to develop a low-cost, rapid, and non-destructive method for analyzing pea lipids with no chemical modifications that would prove superior to existing destructive solvent extraction methods. Different pea accession seed samples, prepared as either small portions (0.5 mm2) of endosperm or ground pea seed powder for comparison, were subjected to HR-MAS NMR analyses and whole seed samples underwent NIR analyses. The total lipid content ranged between 0.57–3.45% and 1.3–2.6% with NMR and NIR, respectively. Compared to traditional extraction with butanol, hexane-isopropanol, and petroleum ether, correlation coefficients were 0.77 (R2 = 0.60), 0.56 (R2 = 0.47), and 0.78 (R2 = 0.62), respectively. Correlation coefficients for NMR compared to traditional extraction increased to 0.97 (R2 = 0.99) with appropriate correction factors. PLS regression analyses confirmed the application of this technology for rapid lipid content determination, with trends fitting models often close to an R2 of 0.95. A better robust NIR quantification model can be developed by increasing the number of samples with more diversity.

1. Introduction

For more than 50 years, society has depended on the organic chemistry industry to produce basic necessities, such as food, pharmaceuticals, and chemicals. Today, 92% of organic products are synthesized in the petrochemical industry using either crude oil or natural gas. Recent issues related to climate change, due to the increase in greenhouse gases, as a result of the petrochemical industry, has led to more research into bio-products.
Although legume or pulse crops are mostly cultivated to enhance protein content for human consumption and animal feed, it may be possible to improve lipid concentration through cross-breeding and genetic engineering to increase their value as enhanced bio-products. Legume protein content ranges from 25–32% [1,2] and lipid content ranges between 0.9–46% [3,4].
Field pea (Pisum sativum) is an important pulse crop that is mostly grown in temperate regions. It is used as either forage or a dried seed supplement in animal feed, and more recently, has been considered, along with other pulse crops, as a potential non-traditional oilseed crop that can help meet global demands for vegetable oil and biodiesel [5]. In addition, field peas are a good source of protein, carbohydrates, phenolic compounds, phytic acid, fiber, lipids, tannins, and amino acids [6]. Peas contain approximately 87% total digestible nutrients, increasing the nutritional value of most diets [7].
Lipids are comprised of fatty acids that are present in the endosperm of the field pea seed [8]. Fatty acid composition is influenced by agro-climatic conditions as the enzymes in fatty acid biosynthesis depend on these factors [9]. Fatty acid content ranges between 0.9–5.0% of dry matter [5]. The low quantity of saturated fatty acids found in peas in comparison to unsaturated fatty acids can be explained by the conversion of the polyunsaturated fatty acids to carbohydrates by lipid peroxidation [10]. This is because the fatty acids contain multiple double bonds and the methyl group that lies within is prone to hydrogen atom abstraction. Research is ongoing to increase the lipid content in field pea using genetic engineering and traditional breeding [11].
Various studies have been conducted using traditional methods for extraction and quantification of lipids, as well as for fatty acid profiling of plant material [12,13,14]. Traditional methods for fatty acid determination involve triacylglycerol hydrolysis, followed by the esterification of the liberated fatty acyl groups with methanol or by direct transesterification and, finally, gas chromatography to separate, identify and quantify fatty acids that are present. This process is required because the majority of the fatty acids exist in the form of esters and amides, and only a small portion constitutes free fatty acids. These traditional methods are labor-intensive and time- and chemical-consuming. They expose the oil to lipid oxidation and the process requires fatty acid standards [15]. Previous solvent extractions have demonstrated percentage lipid content in pea samples ranging from 1.3% ± 0.07% to 2.6% ± 0.11% [16]. More recently, fatty acid determination has evolved with methods that allow for qualitative and quantitative analyses. Modern methods include nuclear magnetic resonance (NMR) and near infrared (NIR) spectroscopies, where spectroscopy is defined as the interaction between matter and electromagnetic radiation.
NMR is based on the absorption of electromagnetic radiation in the radio frequency region by atomic nuclei [17]. This technology has been studied for decades and it has become an important and powerful analytical tool used by researchers [18]. NMR has applications in a wide range of disciplines such as the food, pharmaceutical, chemical and petroleum industries, as well as the agricultural and disease control sectors. In the pharmaceutical industry, it is a useful technique for drug discovery and determining the conformations of the compounds bound to enzymes, receptors, and other proteins [19]. In the agricultural and food industries, NMR is mostly used for quantitative and qualitative analyses, such as structure determination, chemical composition, product purity and adulterations, kinetic reaction monitoring, in situ tests for ripening, etc. NMR can be used to compare the different profiles of genetically modified foods to their unmodified counterparts [20].
NIR spectroscopy is based on the region of the electromagnetic spectrum between visible light and infrared radiation, occurring from approximately 780 nm to 2600 nm [21]. It involves changes in vibration energies of the molecular bonds O-H, C-H, C-O, and N-H when subjected to NIR frequencies, wherein bonds undergo vibrational stretching and bending [22]. Calibration methods have been developed using NIR spectroscopy for rapid determination of seed contents, including lipids, starch, carbohydrates, proteins, moisture, fatty acid composition, phytonutrients, and many other components [23,24]. In the food industry, NIR is used to determine variations in raw materials and finished food products due to different processing procedures [25]. This study aimed to investigate how modern spectroscopy methods, NMR and NIR, could be used to study fatty acids in pea, while developing calibration models to predict pea lipid content in this non-traditional oilseed. Data presented here could advance the development of a low-cost, rapid, and non-destructive method for lipid profiling and determining lipid content in leguminous crops that is superior to traditional destructive solvent extraction methods.

2. Results and Discussion

2.1. HR-MAS NMR Spectroscopy

A review of the literature shows that the NMR spectrum of the pea is similar to those of various oilseeds analyzed using 1H high-resolution magic angle spinning (HR-MAS) NMR under similar conditions [26,27,28]. The 1H NMR spectrum of the pea oil contains 11 distinct signals of variable intensity (Figure 1). Different signals observed in the spectrum are a result of protons in the same chemical environment, where neighbouring protons cause peak splitting. However, the intensity of the signal depends on the proportion of different acyl groups in any given sample [26]. The use of the magic angle spinning property of the HR-MAS probe was efficient at averaging residual dipolar interactions and variations in the bulk magnetic susceptibility of the samples. NMR spectra displayed a resolution approaching that of liquid samples. Some samples achieved lower spectral resolution than others where the splitting of multiplets could not be distinguished.
Signal A (=HC-CH2-CH=) represents bis-allylic hydrogen atoms and is caused by overlapping signals from α-methylene protons in two double bonds. Bis-allylic hydrogens are used to determine the nature of the unsaturated fatty acids present in each sample. Signals B (-OCO-CH2-) and D (-OCO-CH2-CH2-) represent methylene hydrogen atoms in the β and α positions in relation to the carbonyl group, respectively. Signal C (-CH2-CH=CH-) represents allylic hydrogen atoms, caused by α methylene hydrogen atoms in a single double bond (allylic carbon). All methyl protons in saturated or unsaturated acyl groups containing a double bond further from the terminal methyl group give rise to a triplet signal (Signal F). Signal F (-CH3) is produced by the overlap of triplet signals from the methyl group protons of saturated, oleic acyl and linoleic acyl groups [26,27,28].
Signal E (-CH3) represents the methyl protons of ω-3 acyl groups. The difference in chemical shifts between methyl proton signals is due to their proximity to the double bond of the fatty acid chain. Signal G (-(CH2)n-) represents the hydrogen atoms of those methylene groups, either in position β or further from the olefinic groups, or in position γ or further from the carboxylic groups inside the triglyceride molecule. The α hydrogens of the glycerol backbone on carbons 1 and 3 cause signals α1 and α2 (-CH2OCOR). The β (-CHOCOR) signal is due to the β hydrogen atom on carbon 2 of the glycerol backbone. Lastly, signal H (-CH=CH-) overlaps slightly with signal β due to the olefinic hydrogen atoms of the different acyl groups [26,27,28].

2.1.1. Fatty Acid Composition of Pea Seeds Using 1H HR-MAS NMR

The observed chemical shifts of the signals were not significantly (Figure 1) (p > 0.05) different from those reported in the literature [26,27,28]. Importantly, data indicate that the different levels of fatty acid types are not affected by sample preparation into small pea seed portions or as ground pea seed powder (Table 1). The linoleic acid content of the tested pea seed varieties ranged between 13.89–43.98%. Oleic acid, a monounsaturated omega-9 fatty acid, was detected in all field pea varieties tested; the highest content measured was 55.42% (Sample 29579) and the lowest content measured was 21.71% (Sample 42819). Linolenic acid, an essential fatty acid, belongs to the group of omega-3 fatty acids and is highly concentrated in most plant oils. It reportedly inhibits prostaglandin synthesis and reduces the risk of chronic diseases [29]. In this study, the highest linolenic fatty acid content measured was 26.62% (Sample F1) and the lowest content measured was 20.81% (Sample 45760). Saturated fatty acids most commonly found in field peas are arachidic, behenic, palmitic, stearic, and lignoceric fatty acids, with palmitic acid being the most dominant [30]. In this study, saturated fatty acids ranged between 4.61–28.24%. The results are similar to published data [4,8].

2.1.2. Quantification of Total Lipid Content Using 1H HR-MAS NMR

Non-extraction methods for determining total lipid content are indirect methods based on the measurement of a physical or chemical property of a sample. This includes density, dielectric measurements, NIR spectroscopy, NMR spectroscopy, ultrasound, colorimetry and X-ray absorption [31]. Extraction methods are based on the water insolubility property of lipids. This is achieved by using a solvent that is highly lipid-soluble, and less or insoluble for other cell components such as proteins and carbohydrates [32]. Diethyl ether, pentane and methylene chloride are some of the solvents commonly used for extraction of lipids from pea seeds.
Analysis of variances for different extraction methods showed no significant differences (p > 0.05) between Soxhlet extraction with hexane-isopropanol and butanol solvents. In contrast to petroleum ether, significant differences (p > 0.05) existed when a larger amount of lipid was extracted with the hexane-isopropanol and butanol solvents. This is because non-polar lipids exhibit better solubility in hexane-isopropanol and butanol [33]. In addition, these data indicate that a larger fraction of the lipid concentration in pea seeds is non-polar. In Table 2, quantification by both NMR and Soxhlet extraction with petroleum ether underestimated the total lipid content present in field pea seeds compared to Soxhlet extraction with hexane-isopropanol and butanol solvents.
The NMR method employed in this study is rapid and non-destructive compared to traditional solvent extraction methods, and a correlation between the NMR method and other traditional methods was observed. Correlation coefficients of 0.77 (R2 = 0.60), 0.56 (R2 = 0.47), and 0.78 (R2 = 0.62) were observed in the Khodapanahi et al. (2012), Kalia (2016), and Ossowski (2017) studies, respectively (Table 2). Correlation coefficients were determined by comparing the results from the NMR and the three studies. The correlation coefficients for each method increased to 0.97 (R2 of 0.99) when NMR data were corrected with respect to the other methods using a correction factor, namely 3.08, 1.97 and 0.85 for Khodapanahi et al., (2012), Kalia, (2016), and Ossowski, (2017) respectively (Figure 2). The correction factor used was the average ratio of the data acquired for each traditional method to the NMR data. Statistical analyses showed no significant differences (p > 0.05) between the corrected NMR data for NMR and the corresponding quantitative data acquired with traditional methods. This demonstrates that NMR can effectively quantify total lipids in field pea seeds.
The NMR method was used to quantify lipids in soybean to validate the model. Soybean pieces were tested and a total lipid content of 4.60% was recorded by NMR (Table 2). The total lipid content was 14.17% after the correction factor was applied, which was not significantly different when compared to total lipid content determined using the Soxhlet extraction with butanol at 13.90%.

2.1.3. Seed Viability

Seed viability was verified by subjecting all remaining pea seed portions from each variety to a germination test. The germination rate for all pea seeds from each variety (27 in total) was 100% (Table S1). Successful germination was as defined the emergence of the radicle (day four), plumule (day five), foliage leaves (day seven), and secondary roots (day eight). These data indicate that this NMR sampling method may be used for the non-destructive analytical testing of oilseed lipids.

2.2. NIR Spectroscopy

The entire spectrum of absorbance was used to generate the partial least squares (PLS) model (Figure S3) for the near infrared analyses. The PLS was calculated over the selected wavelength range of 1100 nm–2500 nm [21,34]. The near infrared spectra were pre-treated mathematically and derivatized. PLS regression generated a set of b-vector coefficients for each of these points; they were used as coefficients for calculating new predictions from similarly treated spectra. By the nature of the PLS model, quantification of lipids in each pea seed sample was achieved by using spectral features in the entire range from 1100 nm–2500 nm, rather than using the individual peaks or regions in the spectrum.
From the PLS regression analyses (Table 3), different approaches for pre-processing data were assessed for the ability to produce models that best fit the acquired data. Results describing the validation models produced are summarized in Table 3. When using NIR spectroscopy for oil determination in pea seeds [21], the use of PLS regression is recommended over other multivariate tests, such as principal component regression. Here, the effectiveness of PLS is demonstrated through the coefficient of determination (R2) and the standard error of cross-validation (SECV), which is the error of cross validation when corrected for bias using the CAMO Unscrambler® X Software (CAMO Analytics, Oslo, Norway).
A basic PLS regression analysis of the data without manipulation beyond conversion to absorbance demonstrates a four-factor model with a reasonable SECV of 0.13 (R2 = 0.96). Improvement to the model’s fitness is demonstrated when the data is pre-processed. An interesting metric of the effectiveness of a model is the calculation of RPD, or the ratio of standard error of performance to standard deviation. RPD is calculated as the ratio of the standard derivative to the SECV. RPD values can be used to predict the approximate suitability of data for use at different levels. The higher the RPD, the more reliable is the calibration. The model is suitable for rough scans and general approximation of oil content as the values approach or exceed 4.5 RPD. There is no doubt that the models produced are insufficiently robust, as only 18 samples were used for the study. Performance of this caliber with such a small sample set would suffer from over-fitting. Modest results using NIR for component quantification with larger sample counts have reported lower RPD scores [35]. A greater number of samples would lend to a more robust model. The overall deviation between the reference values versus the calibration and validation models is seen in Figure 3.

3. Materials and Methods

3.1. HR-MAS NMR Spectroscopy

3.1.1. Pea Accessions and Sample Preparation

Pea (Pisum sativum) oil, soybean seeds (Glycine max, cv. Champion), and nine wild and domestic pea accessions catalogued as 29579, 112351, 42819, F2 (G611 x Wando), F1 (112351 × 43016), 29600, 43016, 45760, and 29526 were used in this study. Field pea seeds were harvested in 2015 from the Harbin Seed Farm (Rivercourse, AL, Canada) and McGill University’s Macdonald campus (Sainte-Anne-de-Bellevue, QC, Canada). Soybean seeds were obtained from plants grown in 2015 at the Belcan Agro Centre (Sainte-Marthe, QC, Canada).
Samples were prepared as two different forms. For pea seed portions, a small fragment (0.5 mm2) of the endosperm of each seed was manually removed to serve as the pea seed sample for NMR analysis. Germination tests were performed on the remaining pea seed fraction to determine viability for each variety and confirm that the method used was non-destructive. Ground pea samples were achieved by grinding dry pea seeds to a fine powder with a mortar and pestle prior to NMR analysis. Serving as the external reference for total lipid quantification, one pea oil sample was extracted from pea accession 29,526 using the Soxhlet extraction method with 12 mL hexane-isopropanol (3:2, v/v) at room temperature.

3.1.2. Instrumental Setup

Fatty acid profiling and total lipid content of pea seeds were determined using the high-resolution magic angle spinning (HR-MAS) NMR equipment at the Research Institute of the McGill University Health Centre (Montreal, QC, Canada). These experiments were performed using a Bruker AVANCE III HD 600 spectrometer (Bruker Cooperation, Billerica, MA, United States) operating at 600.17 MHz for 1H. A 4-mm triple-resonance 1H/31P/13C high-resolution MAS probe with a Z-gradient directed along the magic angle axis was employed. All seed samples were spun at 5000 Hz in zirconia rotors and experiments were run at room temperature (~20 °C). Approximately 15 μL deuterium oxide (D2O) containing 5 mM 3-(trimethylsilyl) propionic-2, 2, 3, 3-d4 acid (D2O + TSP; Sigma-Aldrich, Oakville, ON, Canada) was pipetted into the bottom of the rotor to serve as an internal standard for chemical shift references. Approximately 7 mg of each sample (pea seed portion or ground sample) were placed into a Bruker KelF insert and combined with 15μL D2O. The insert was then closed with a plug and a sealing screw. The insert was placed in the rotor and 1H spectra were acquired with a 1-D pulse sequence and water pre-saturation (Bruker pulse sequence zgpr). Each spectrum was acquired with 128 scans, a 6 μs 90° pulse length, 8417.5 Hz spectral width, 16k data points, a 0.97 s acquisition time, and a 1 s relaxation delay. Solution-type NMR experiments were achieved for the study by spinning at the magic angle using the high-resolution magic angle spinning (HR-MAS) probe. Only 1H experiments were acquired, as is common when using HR-MAS probes. The experiment was conducted in triplicate for each accession.

3.1.3. HR-MAS NMR Spectra Processing

The NMR method used for lipid quantification was determined by knowing that all fatty acid acyl chains were esterified to a common moiety, a glycerol backbone that forms a triacylglyceride [36]. The area under each signal in the proton NMR spectrum is proportional to the number of hydrogen atoms for the different types of aliphatic groups of the fatty acid chain [37]. Integration of the signals was done using Topspin 4.0.3 software (Bruker Cooperation, Billerica, MA, USA) to find the area values under the signals. Preliminary tests showed that TSP was not a reliable reference for the quantification of fatty acids. The accuracy of the TSP quantification method can be affected by the binding of TSP to fatty acids and the displacement of TSP into small spaces inside the spinning rotor that are outside of the radio frequency coil’s detectable volume [38,39]. The pea oil sample was used as the quantification reference and for normalizing the sum of the integrals to 100. The spectrum for each field pea seed sample (Figure S1) was compared to the pea oil sample using the multiple integration NMR spectra function, generating integral values (Table S2). This was used to calculate proportions of the different acyl groups in each sample using the following equations:
Linolenic fatty acid (C18:3) (%) = 100% [E/(F + E)]
Linoleic fatty acid (C18:2) (%) = 100% [(A/B) − 2[E/(F + E)]
Oleic fatty acid(C18:1) (%) = 100% [(C/2B) − (A/B) + [E/(F + E)]
Saturated fatty acids (%) = 100% [1 − (C/2B)]
where A represents HC-CH2-CH acyl groups; B represents -OCO-CH2- acyl groups; C represents -CH2-CH=CH- acyl groups; D represents -OCO-CH2-CH2- acyl groups; E represents a -CH3 linolenic acyl group; and F represents -CH3 saturated, oleic and linoleic acyl groups [37].

3.1.4. Statistical Analyses of HR-MAS NMR Data

SAS University Edition vApp (SAS Institute Inc., Cary, NC, US) was used for data analyses. A completely random design model of two factors (variety and processed sample form) and Student–Newman–Keuls pairwise comparison was used for this study.

3.2. NIR Spectroscopy

Six wild pea accessions, 29526, 29579, 42819, 43016, 45760, and 112351, were tested for NIR analyses. Pea plants were field grown and harvested in the summers of 2014 and 2015 at McGill University’s Macdonald campus and dried in an oven at 40 °C for 72 h before long-term storage.

3.2.1. Instrumental Setup and Scan Parameters

NIR scans were taken on a FieldSpec 4 Standard-Res portable spectroradiometer (ASD Incorporated, Longmont, CO, United States). The FieldSpec 4 was equipped with an ASD Contact Probe (ASD Incorporated, Longmont, CO, USA), a 10-mm scan spot size and a halogen bulb light source (2900 K colour temperature). Scans were taken over the full instrument spectrum of 350 nm to 2500 nm. The room where scans took place was environmentally lit with standard fluorescent lighting, rated at a colour temperature of 3500 K. The instrument was calibrated with a Spectralon fluoropolymer calibration standard with a diffuse reflectance of >95% between 250 nm and 2500 nm. Clear Petri dishes (100 mm × 15 mm) were filled and levelled with whole pea samples. The experiment was conducted in triplicate for each accession.
A turntable was designed from the rotating head of a tube rotator adjusted to sit horizontally (Figure S2). A Styrofoam tray was modified and attached to the head of the rotator to center and secure Petri dishes containing each whole pea sample. Samples were rotated at 22 rpm, with the instrument’s probe fixated approximately 30 mm from the center of each dish, 5 mm above the sample. For each sample dish, the instrument was programmed to take average reflectance values of 100 scans over a 17-s period as the sample dish rotated. Each Petri dish containing one sample was scanned and analyzed chemometrically using UCal 4TM Custom Calibration Software (Unity Scientific, Milford, MA, USA).

3.2.2. Data Pre-Processing

NIR spectra were analyzed and processed using UCal 4TM Custom Calibration Software. A data matrix was configured to associate the reference lipid values determined previously [16] with their corresponding spectra. Reflectance spectra were converted to absorbance by taking the logarithm of the inverse of the reflectance values. Applicable wavelengths were consolidated to the NIR range of 1100–2500 nm. Standard normal variate and math treatments were performed on the full sample set to assess their ability to remove additive baseline and scattering effects. Multiplicative scattering correction was applied to reduce the effects of light scattering produced by the sampling method. This eliminated deleterious spectral effects that might have been introduced by the instrument setup and not the samples, such as slight changes in ambient lighting or slight movements of the scanning instrument in three-dimensional space during and between scans.

3.2.3. Calibration and Validation Model Design

Partial least squares (PLS) regression analyses were run for datasets that were pre-processed in the preceding section using UCal 4TM Custom Calibration Software. The calibration model design relied on data previously obtained as reference values for lipid content in the scanned pea samples obtained by solvent extraction using hexane-isopropanol (3:2, v/v) in a modified version of a method described previously [4]. In the absence of an independent validation sample set due to the small sample size available for this study, a cross-validation method was used for validating the calibration model. A leave-one-out cross-validation (LOOCV) was performed to generate the standard error of cross-validation (SECV) statistics, and each sample was used as a singular test set against the model built from the remaining training set.

4. Conclusions

Traditional methods for lipid quantification require several grams of sample, costly extraction solvents, and are time-consuming. The results showed that NMR and NIR are low-cost, rapid, and non-destructive methods for analyzing lipids in pea. With this study, we confirmed that NMR technology can be used for quantitative and qualitative analyses of the fatty acids present in oilseeds. These data show a high correlation between the NMR method and the traditional method in predicting the lipid content present in seeds. In comparing traditional extraction methods, however, an appropriate correction factor must be used with the measured NMR values. PLS regression analyses confirm the possible application of NIR analysis for the rapid determination of pea lipid content. With trends fitting models often close to R2 = 0.96, more samples and more plant cultivars will improve the robustness of the calibration model and an effective system may be developed to profile fatty acids and quantify lipids in various samples. Due to the versatility of NMR and NIR, future studies should explore the use of these technologies for the nutritional assessment of food products and the identification of adulterations or impurities.

Supplementary Materials

The following are available online, Table S1. Germination rate of pea seeds, Table S2. Average integral values of the peaks used for fatty acid quantification, Figure S1. 1H NMR spectrum of the pea seed, and Figure S2. Detailed view of the contact probe, light source, sample, and turntable sample trya for NIR analyses. Figure S3. NIR spectra from six pea cultivars.

Author Contributions

Conceptualization, M.L., M.-J.D., V.O. and J.S.; methodology, P.W.A., P.O., A.E.G. and R.K.; software, A.E.G., P.W.A. and P.O.; formal analysis, P.W.A. and P.O.; visualization, P.W.A.; investigation, P.W.A. and P.O.; writing—original draft preparation, P.W.A. and P.O.; writing—review and editing, S.M. and V.H.-V.; supervision, M.L.; funding acquisition, M.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Natural Sciences and Engineering Research Council, grant number NSERC CRDPJ 477365-14.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We extend our sincere gratitude to Harbin Seed Farm, Rivercourse, Alberta, Canada and Belcan Agro Centre, Sainte-Marthe, Quebec, Canada for their support. We are grateful to Andrée Gravel for her laboratory support.

Conflicts of Interest

The authors declare that they have no conflict of interest.

Sample Availability

Samples of the compounds are available from the authors.

References

  1. Iqbal, A.; Khalil, I.A.; Ateeq, N.; Khan, M.S. Nutritional quality of important food legumes. Food Chem. 2006, 97, 331–335. [Google Scholar] [CrossRef]
  2. Jezierny, D.; Mosenthin, R.; Bauer, E. The use of grain legumes as a protein source in pig nutrition: A review. Anim. Feed Sci. Technol. 2010, 157, 111–128. [Google Scholar] [CrossRef]
  3. Karn, A.; Heim, C.; Flint-Garcia, S.; Bilyeu, K.; Gillman, J. Development of rigorous fatty acid near-infrared spectroscopy quantitation methods in support of soybean oil improvement. J. Am. Oil Chem. Soc. 2017, 94, 69–76. [Google Scholar] [CrossRef]
  4. Ryan, E.; Galvin, K.; O’Connor, T.P.; Maguire, A.R.; O’Brien, N.M. Phytosterol, squalene, tocopherol content and fatty acid profile of selected seeds, grains, and legumes. Plant Foods Hum. Nutr. 2007, 62, 85–91. [Google Scholar] [CrossRef] [PubMed]
  5. Khodapanahi, E.; Lefsrud, M.; Orsat, V.; Singh, J.; Warkentin, T.D. Study of pea accessions for development of an oilseed pea. Energies 2012, 5, 3788–3802. [Google Scholar] [CrossRef]
  6. Sarwar, M.F.; Sarwar, M.H.; Sarwar, M.; Qadri, N.A.; Moghal, S. The role of oilseeds nutrition in human health: A critical review. J. Cereal Oil 2013, 4, 97–100. [Google Scholar] [CrossRef] [Green Version]
  7. Kocer, A.; Albayrak, S. Determination of forage yield and quality of pea (Pisum sativum L.) mixtures with oat and barley. Turkish J. Field Crop. 2012, 17, 96–99. [Google Scholar]
  8. Solis, M.I.V.; Patel, A.; Orsat, V.; Singh, J.; Lefsrud, M. Fatty acid profiling of the seed oils of some varieties of field peas (Pisum sativum) by RP-LC/ESI-MS/MS: Towards the development of an oilseed pea. Food Chem. 2013, 139, 986–993. [Google Scholar] [CrossRef]
  9. Yoshida, H.; Tomiyama, Y.; Tanaka, M.; Mizushina, Y. Characteristic profiles of lipid classes, fatty acids and triacylglycerol molecular species of peas (Pisum sativum L.). Eur. J. Lipid Sci. Technol. 2007, 109, 600–607. [Google Scholar] [CrossRef]
  10. Ayala, A.; Muñoz, M.F.; Argüelles, S. Lipid peroxidation: Production, metabolism, and signaling mechanisms of malondialdehyde and 4-hydroxy-2-nonenal. Oxid. Med. Cell. Longev. 2014, 2014, 1–32. [Google Scholar] [CrossRef]
  11. Ahmad, S.; Kaur, S.; Lamb-Palmer, N.D.; Lefsrud, M.; Singh, J. Genetic diversity and population structure of Pisum sativum accessions for marker-trait association of lipid content. Crop J. 2015, 3, 238–245. [Google Scholar] [CrossRef] [Green Version]
  12. Caprioli, G.; Giusti, F.; Ballini, R.; Sagratini, G.; Vila-Donat, P.; Vittori, S.; Fiorini, D. Lipid nutritional value of legumes: Evaluation of different extraction methods and determination of fatty acid composition. Food Chem. 2016, 192, 965–971. [Google Scholar] [CrossRef] [PubMed]
  13. Costa, C.S.; Vicenti, J.R.; Moron-Villarreyes, J.A.; Caldas, S.; Cardoso, L.V.; Freitas, R.F.; D’oca, M.G. Extraction and characterization of lipids from Sarcocornia ambigua meal: A halophyte biomass produced with shrimp farm effluent irrigation. An. Acad. Bras. Ciênc. 2014, 86, 935–943. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Padhi, E.M.; Liu, R.; Hernandez, M.; Tsao, R.; Ramdath, D.D. Total polyphenol content, carotenoid, tocopherol and fatty acid composition of commonly consumed Canadian pulses and their contribution to antioxidant activity. J. Funct. Food. 2017, 38, 602–611. [Google Scholar] [CrossRef]
  15. Mazurek, B.; Chmiel, M.; Górecka, B. Fatty acids analysis using gas chromatography-mass spectrometer detector (GC/MSD)-method validation based on berry seed extract samples. Food Anal. Method. 2017, 10, 2868–2880. [Google Scholar] [CrossRef] [Green Version]
  16. Kalia, R. Correlation of Lipid Content and Phenotypic Markers of Canadian Field Peas (Pisum sativum); McGill University Libraries: Montreal, QC, Canada, 2016. [Google Scholar]
  17. Kleinberg, R.L.; Jackson, J.A. An introduction to the history of NMR well logging. Concepts Magn. Res. 2001, 13, 340–342. [Google Scholar] [CrossRef]
  18. Moser, E.; Laistler, E.; Schmitt, F.; Kontaxis, G. Ultra-High Field NMR and MRI-The Role of Magnet Technology to Increase Sensitivity and Specificity. Front. Phys. 2017, 5, 33–48. [Google Scholar] [CrossRef] [Green Version]
  19. Baias, M.; Dumez, J.-N.; Svensson, P.H.; Schantz, S.; Day, G.M.; Emsley, L. De novo determination of the crystal structure of a large drug molecule by crystal structure prediction-based powder NMR crystallography. J. Am. Chem. Soc. 2013, 135, 17501–17507. [Google Scholar] [CrossRef] [Green Version]
  20. Sobolev, A.P.; Testone, G.; Santoro, F.; Nicolodi, C.; Iannelli, M.A.; Amato, M.E.; Ianniello, A.; Brosio, E.; Giannino, D.; Mannina, L. Quality traits of conventional and transgenic lettuce (Lactuca sativa L.) at harvesting by NMR metabolic profiling. J. Agric. Food Chem. 2010, 58, 6928–6936. [Google Scholar] [CrossRef]
  21. Baianu, I.; You, T.; Costescu, D.; Lozano, P.; Prisecaru, V.; Nelson, R.L. Determination of soybean oil, protein and amino acid residues in soybean seeds by high resolution nuclear magnetic resonance (NMRS) and near infrared (NIRS). Nature 2012, 1, 1–62. [Google Scholar] [CrossRef]
  22. Cen, H.; He, Y. Theory and application of near infrared reflectance spectroscopy in determination of food quality. Trends Food Sci. Technol. 2007, 18, 72–83. [Google Scholar] [CrossRef]
  23. Shi, H.; Lei, Y.; Prates, L.L.; Yu, P. Evaluation of near-infrared (NIR) and Fourier transform mid-infrared (ATR-FT/MIR) spectroscopy techniques combined with chemometrics for the determination of crude protein and intestinal protein digestibility of wheat. Food Chem. 2019, 272, 507–513. [Google Scholar] [CrossRef] [PubMed]
  24. Zhu, Z.; Chen, S.; Wu, X.; Xing, C.; Yuan, J. Determination of soybean routine quality parameters using near-infrared spectroscopy. Food Sci. Nutr. 2018, 6, 1109–1118. [Google Scholar] [CrossRef] [PubMed]
  25. Porep, J.U.; Kammerer, D.R.; Carle, R. On-line application of near infrared (NIR) spectroscopy in food production. Trends Food Sci. Technol. 2015, 46, 211–230. [Google Scholar] [CrossRef]
  26. Barison, A.; Pereira da Silva, C.W.; Campos, F.R.; Simonelli, F.; Lenz, C.A.; Ferreira, A.G. A simple methodology for the determination of fatty acid composition in edible oils through 1H NMR spectroscopy. Magn. Reson. Chem. 2010, 48, 642–650. [Google Scholar] [CrossRef] [PubMed]
  27. Salinero, C.; Feás, X.; Mansilla, J.P.; Seijas, J.A.; Vázquez-Tato, M.P.; Vela, P.; Sainz, M.J. 1H-nuclear magnetic resonance analysis of the triacylglyceride composition of cold-pressed oil from Camellia japonica. Molecules 2012, 17, 6716–6727. [Google Scholar] [CrossRef] [Green Version]
  28. Guillén, M.D.; Ruiz, A. Monitoring of heat-induced degradation of edible oils by proton NMR. Eur. J. Lipid Sci. Technol. 2008, 110, 52–60. [Google Scholar] [CrossRef]
  29. Lanzmann-Petithory, D. Alpha-linolenic acid and cardiovascular diseases. J. Nutr. Health Aging 2001, 5, 179–183. [Google Scholar]
  30. Carta, G.; Murru, E.; Banni, S.; Manca, C. Palmitic acid: Physiological role, metabolism and nutritional implications. Front. Physiol. 2017, 8, 902. [Google Scholar] [CrossRef] [Green Version]
  31. Xu, J.L.; Riccioli, C.; Sun, D.W. An overview on nondestructive spectroscopic techniques for lipid and lipid oxidation analysis in fish and fish products. Compr. Rev. Food Sci. Food Saf. 2015, 14, 466–477. [Google Scholar] [CrossRef] [Green Version]
  32. Moreau, R.A.; Powell, M.J.; Singh, V. Pressurized liquid extraction of polar and nonpolar lipids in corn and oats with hexane, methylene chloride, isopropanol, and ethanol. J. Oil Am. Chem. Soc. 2003, 80, 1063–1067. [Google Scholar] [CrossRef]
  33. Ossowski, P. An Overview of Results Associated with the Experiments Performed for the Development of an Oilseed Pea; McGill University Libraries: Montreal, QC, Canada, 2017. [Google Scholar]
  34. Baianu, I.; Guo, J. NIR Calibrations for Soybean Seeds and Soy Food Composition Analysis: Total Carbohydrates, Oil, Proteins and Water Contents. Nat. Preced. 2011. [Google Scholar] [CrossRef]
  35. Armstrong, P. Rapid single-kernel NIR measurement of grain and oil-seed attributes. Appl. Eng. Agric. 2006, 22, 767–772. [Google Scholar] [CrossRef] [Green Version]
  36. Prajapati, H.N.; Dalrymple, D.M.; Serajuddin, A.T. A comparative evaluation of mono-, di-and triglyceride of medium chain fatty acids by lipid/surfactant/water phase diagram, solubility determination and dispersion testing for application in pharmaceutical dosage form development. Pharm. Res. 2012, 29, 285–305. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Guillén, M.D.; Ruiz, A. 1H nuclear magnetic resonance as a fast tool for determining the composition of acyl chains in acylglycerol mixtures. Eur. J. Lipid Sci. Technol. 2003, 105, 502–507. [Google Scholar] [CrossRef]
  38. Albers, M.J.; Butler, T.N.; Rahwa, I.; Bao, N.; Keshari, K.R.; Swanson, M.G.; Kurhanewicz, J. Evaluation of the ERETIC method as an improved quantitative reference for 1H HR-MAS spectroscopy of prostate tissue. Magn. Reson. Med. 2009, 61, 525–532. [Google Scholar] [CrossRef] [Green Version]
  39. Bharti, S.K.; Roy, R. Quantitative 1H NMR spectroscopy. TrAC-Trend. Anal. Chem. 2012, 35, 5–26. [Google Scholar] [CrossRef]
Figure 1. 1H NMR spectrum of the pea oil sample.
Figure 1. 1H NMR spectrum of the pea oil sample.
Molecules 27 01642 g001
Figure 2. Comparing of the actual results of the traditional methods to the actual and corrected results of the NMR of field pea.
Figure 2. Comparing of the actual results of the traditional methods to the actual and corrected results of the NMR of field pea.
Molecules 27 01642 g002
Figure 3. Graph of predicted (NIR) versus reference values demonstrating fitness of the model.
Figure 3. Graph of predicted (NIR) versus reference values demonstrating fitness of the model.
Molecules 27 01642 g003
Table 1. Sample preparation comparison and percentage (% dry mass) of different fatty acids in field pea seeds.
Table 1. Sample preparation comparison and percentage (% dry mass) of different fatty acids in field pea seeds.
SampleFormLinolenic Fatty Acid (%)Linoleic Fatty Acid (%)Oleic Fatty Acid (%)Saturated Fatty Acid (%)
F1Seed portion26.62 ± 1.65 a23.57 ± 4.60 c44.49 ± 3.86 a5.32 ± 4.26 b
Ground24.56 ± 0.56 a20.28 ± 5.07 c45.64 ± 5.52 a9.51 ± 0.90 b
112351Seed portion23.84 ± 0.56 ab31.42 ± 1.42 b32.33 ± 0.65 bc12.42 ± 1.42 b
Ground22.35 ± 2.69 ab30.81 ± 0.77 b34.89 ± 0.54 bc11.96 ± 2.63 b
F2Seed portion22.29 ± 1.15 b34.16 ± 2.40 ab31.98 ± 1.10 bc11.56 ± 3.05 b
Ground22.95 ± 2.15 b36.50 ± 3.60 ab29.28 ± 2.31 bc11.28 ± 0.20 b
42819Seed portion23.42 ± 0.52 ab43.98 ± 7.73 a21.71 ± 7.72 c10.89 ± 0.21 b
Ground23.86 ± 0.26 ab35.79 ± 3.45 a30.59 ± 3.77 c9.76 ± 0.79 b
43016Seed portion20.88 ± 0.86 b35.30 ± 5.24 ab34.01 ± 4.90 bc9.81 ± 0.69 b
Ground21.47 ± 0.92 b35.78 ± 5.32 ab32.29 ± 6.18 bc10.46 ± 0.45 b
29600Seed portion23.27 ± 1.35 ab31.13 ± 3.05 b40.99 ± 3.42 ab4.61 ± 2.14 b
Ground22.85 ± 1.61 ab28.97 ± 7.28 b41.03 ± 7.25 ab7.15 ± 1.25 b
45760Seed portion21.24 ± 1.64 b36.63 ± 1.33 ab31.37 ± 1.30 bc10.76 ± 1.63 b
Ground20.81 ± 2.72 b34.81 ± 3.13 ab31.71 ± 3.81 bc12.66 ± 0.94 b
29579Seed portion24.08 ± 0.58 ab20.19 ± 0.52 c44.63 ± 3.75 a11.09 ± 3.24 b
Ground22.61 ± 2.85 ab13.89 ± 3.45 c55.42 ± 1.17 a8.08 ± 4.74 b
29526Seed portion25.52 ± 1.04 a21.00 ± 7.40 c25.24 ± 3.57 bc28.24 ± 4.50 a
Ground25.45 ± 7.46 a18.42 ± 7.85 c34.37 ± 3.52 bc21.76 ± 2.83 a
Values in the same column with the same letter superscript are not significantly different (p < 0.05). Values are the means of three replicates ± standard deviation.
Table 2. Comparative table of the actual total fatty acid (% dry mass) vs the corrected total fatty acid (% dry mass) recorded for soybean and field pea seed accessions.
Table 2. Comparative table of the actual total fatty acid (% dry mass) vs the corrected total fatty acid (% dry mass) recorded for soybean and field pea seed accessions.
Samples Actual NMR
Data
Khodapanahi et al. (2012)Kalia (2016)Ossowski (2017)
Actual Soxhlet (Butanol) Corrected NMR Actual Soxhlet (Hexane-Isopropanol)Corrected NMR Actual Soxhlet (Petroleum Ether)Corrected NMR
Soybean4.60 ± 0.23 b13.9014.17NRNRNRNR
F10.80 ± 0.06 aNRNRNRNRNRNR
1123510.67 ± 0.05 a1.602.061.401.320.610.57
F21.09 ± 0.05 aNRNRNRNRNRNR
428190.77 ± 0.10 a2.502.371.701.520.670.65
430161.09± 0.09 a3.603.362.202.15NRNR
296001.12 ± 0.08 a2.603.451.702.21NRNR
457601.07 ± 0.12 a3.503.302.102.111.050.91
295791.02 ± 0.06 a3.703.142.602.010.740.87
295260.93 ± 0.10 a3.102.861.301.830.690.79
Columns followed by the same letter are not significantly different (p < 0.05). Values represent the means of three replicates ± standard deviation. NR: not reported.
Table 3. PLS regression model cross validation data.
Table 3. PLS regression model cross validation data.
Pre-Processing Step(s)FactorsR2SECVReference Std. Dev.RPD
[a]40.960.130.503.97
[b]30.980.110.504.62
[b] [c]40.990.110.504.62
[a] No math treatment; [b] Standard normal variate and detrend treatment, forward derivatization (1st order); [c] Multiplicative Scattering Correction; SECV-Standard of error for cross-validation; RPD-Residual predictive deviation.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Addo, P.W.; Ossowski, P.; MacPherson, S.; Gravel, A.E.; Kaur, R.; Hoyos-Villegas, V.; Singh, J.; Orsat, V.; Dumont, M.-J.; Lefsrud, M. Development of a Nuclear Magnetic Resonance Method and a Near Infrared Calibration Model for the Rapid Determination of Lipid Content in the Field Pea (Pisum sativum). Molecules 2022, 27, 1642. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules27051642

AMA Style

Addo PW, Ossowski P, MacPherson S, Gravel AE, Kaur R, Hoyos-Villegas V, Singh J, Orsat V, Dumont M-J, Lefsrud M. Development of a Nuclear Magnetic Resonance Method and a Near Infrared Calibration Model for the Rapid Determination of Lipid Content in the Field Pea (Pisum sativum). Molecules. 2022; 27(5):1642. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules27051642

Chicago/Turabian Style

Addo, Philip Wiredu, Philip Ossowski, Sarah MacPherson, Andrée E. Gravel, Rajvinder Kaur, Valerio Hoyos-Villegas, Jaswinder Singh, Valérie Orsat, Marie-Josée Dumont, and Mark Lefsrud. 2022. "Development of a Nuclear Magnetic Resonance Method and a Near Infrared Calibration Model for the Rapid Determination of Lipid Content in the Field Pea (Pisum sativum)" Molecules 27, no. 5: 1642. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules27051642

Article Metrics

Back to TopTop