Evaluation of FT-NIR and ATR-FTIR Spectroscopy ... - ACS Publications

their use as diagnostic tools in dairy herd and human health management. Here .... Peter Adewale , Ogan Mba , Marie-Josée Dumont , Michael Ngadi ...
0 downloads 0 Views 985KB Size
Subscriber access provided by DUKE UNIV

Article

EVALUATION OF FT-NIR AND ATR-FTIR SPECTROSCOPY TECHNIQUES FOR DETERMINATION OF MINOR ODD AND BRANCHED CHAIN SATURATED AND TRANS UNSATURATED MILK FATTY ACIDS Ivan Valentinov Stefanov, Vincent Paul Baeten, Ouissam Abbas, Bruno Vlaeminck, Bernard De Baets, and Veerle Fievez J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/jf304515v • Publication Date (Web): 19 Feb 2013 Downloaded from http://pubs.acs.org on March 11, 2013

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Journal of Agricultural and Food Chemistry is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 44

Journal of Agricultural and Food Chemistry

1

EVALUATION OF FT-NIR AND ATR-FTIR SPECTROSCOPY TECHNIQUES FOR

2

DETERMINATION OF MINOR ODD AND BRANCHED CHAIN SATURATED AND

3

TRANS UNSATURATED MILK FATTY ACIDS

4

I. STEFANOV1, V. BAETEN2, O. ABBAS2, B. VLAEMINCK1, B. DE BAETS3 and V.

5

FIEVEZ*1

6

1

7

Bioscience Engineering, Ghent University, Proefhoevestraat 10, B-9090 Melle, Belgium

8

2

9

Agricultural Research Centre. Chaussée de Namur, 24. B-5030 Gembloux, Belgium

Laboratory for Animal Nutrition and Animal Product Quality (LANUPRO), Faculty of

Food and Feed Quality Unit, Valorisation of Agricultural Products Department, Walloon

10

3

11

Engineering, Ghent University, Coupure links 653, B-9000 Ghent, Belgium

Department of Mathematical Modeling, Statistics and Bioinformatics, Faculty of Bioscience

12 13 14 15 16

Corresponding author: Veerle Fievez, Tel: +32 9 264 90 01, Fax: +32 9 264 90 99

17

E-mail address: [email protected]

18 19 20 21 22 23

1 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 2 of 44

24

ABSTRACT

25

Determination of nutritionally important trans MUFA, CLA and OBCFA milk fatty acids (often

26

present in amounts lower than 1.0 g/100g of total fat) using fast and non-destructive analytical

27

methods would enhance their use as diagnostic tools in dairy herd and human health

28

management. Here, PLS regression using ATR/FTIR spectra indicated potential for

29

determination of trans-11 C18:1 and trans-12 C18:1 (Rcv2≥0.80), trans-9 C18:1 in very minor

30

concentration (Rcv2>0.82), as well as for anteiso C15:0 (Rcv2=0.57) and iso C17:0 (Rcv2=0.61).

31

Furthermore, the main cis-9, trans-11 CLA isomer was predicted well despite the high trans

32

MUFA concentration. Differentiation between the CLA and the trans MUFA signals was evident

33

(based on specific cis/trans bands), as well as branched chain saturated fatty acid methyl esters

34

revealed specific iso and anteiso ATR/FTIR absorbance bands. None of the minor FA PLS

35

results with FT-NIR showed interesting potential, except satisfactory predictions for trans-9

36

C18:1 and cis-9, trans-11 CLA. Overall ATR/FTIR resulted in better calibrations and provided

37

more specific information for determination of minor milk fatty acids.

38 39

KEYWORDS: ATR, FTIR, FT-NIR, milk fat, trans fat, CLA, OBCFA, spectroscopy

40 41 42 43 44 45 46

2 ACS Paragon Plus Environment

Page 3 of 44

Journal of Agricultural and Food Chemistry

47

INTRODUCTION

48

In addition to the already established nutritional concerns for trans monounsaturated fatty acids

49

(trans MUFA) and the beneficial biological effects of conjugated linoleic acids (CLA) for human

50

health 1, these two groups of milk fatty acids together with another group of specific odd and

51

branched chain saturated fatty acids (OBCFA) have emerged as important nutritional biomarkers

52

for dairy herds

53

duodenal flow of microbial biomass, as well as the potential of specific OBCFA isomers to

54

monitor the nutrients produced during digestive processes would allow their utilization as

55

excellent diagnostic tools of rumen function

56

MUFA, CLA and OBCFA for early diagnosis of digestive disorders in cows, such as diagnosis

57

of subacute ruminal acidosis (SARA) have recently been identified as isomer specific 3. For

58

example, sudden increase in trans-10 C18:1, which is part of the trans-10, cis-12 C18:2

59

biohydrogenation pathway 7, has been associated with low ruminal pH, thus both fatty acids

60

might be indicative of SARA. In addition, the cis-9, trans-11 C18:2 isomer was recently

61

discovered to be an effective SARA predictor 3. Moreover, drops in the iso C13:0 to iso C16:0

62

concentrations in milk fat were also indicative of SARA induction. The latter is thought to be

63

caused by an OBCFA productivity shrinkage (production of total OBCFA) from ruminal

64

cellulolytic microbiota and/or lower competitiveness of the cellulolytic bacteria relative to the

65

rest of the microbial community (decrease in the cellulolytic/total bacteria ratio).

66

Thus, the ability to follow the variations of these interesting isomers in milk fat would be a

67

valuable non-invasive strategy for detection of altered ruminal environment. Still implementation

68

of the latter requires reliable and effective analytical techniques, which can simultaneously

69

quantify and differentiate minor milk fatty acids.

2, 3

. The microbial origin of OBCFA and their direct relationship with the

4, 5, 6

. Further, the monitoring potential of trans

3 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 4 of 44

70

It would not be trivial to mention that gas chromatography (GC) and high-performance liquid

71

chromatography are the reference methods, as they are still the golden standard and the most

72

widely used analytical techniques for determination of almost any lipid species. However, the

73

effectiveness for disease prevention using this strategy based on quantification of minor fatty

74

acids would be drastically enhanced by the availability of fast, cheap and non-destructive

75

analytical methods, such as spectrophotometry techniques. Developing spectroscopy based

76

methods for the simultaneous determination of several different types of minor fatty acids with

77

different chemical structures would be of major importance for animal as well as for human

78

nutrition.

79

Fourier transform near-infrared (FT-NIR) and attenuated total reflectance infrared (ATR/FTIR)

80

are well established spectroscopy techniques, both of which are widely used in the food industry

81

8

82

official AOCS (Cd 14d-99) method for quantifying total trans fats 9, 10, due to unique property to

83

detect the absorption of the C-H deformation vibration from a trans C=C bond in the vicinity of

84

966 cm-1 11. The sensitivity of the ATR/FTIR technique allows accurate determination of the total

85

monounsaturated trans content in various naturally occurring non-ruminant fats, vegetable and

86

hydrogenated oils and fats (e.g. canola, margarine and lard), but it has not been extensively

87

tested with more complex fat matrices, such as dairy fat.

88

FT-NIR is the more widely used technique in various niches of the food industry, due to its

89

property of a strong absorption signal from many types of molecules and the ability for

90

transmission of the NIR radiation by optical fibres over long distances. However, the former

91

feature has also brought FT-NIR the legacy of low specificity and inability to differentiate

92

between individual components in various homogeneous matrices, due to combinations and

. ATR/FTIR is famous as the only spectroscopy analytical technique to be approved as an

4 ACS Paragon Plus Environment

Page 5 of 44

Journal of Agricultural and Food Chemistry

93

overlaps of broad bands from fundamental and overtone chemical vibrations (generally

94

overlapping vibrational bands that are non specific and non-resolved)

95

NIR was previously reported as a viable technique for determination of the complete fatty acid

96

profile in soybean, olive and flaxseed oils, as well as in different mixtures of shortening and lard

97

14

98

might allow trans FA determinations for regulatory labelling purposes 14. The use of FT-NIR in

99

bovine milk fatty acid analyses has not been investigated.

12, 13

. Nevertheless, FT-

. The latter has permitted FT-NIR to be used for rapid screening/monitoring of fat products and

100

Milk fat is a complex matrix and in addition to the many different types of fatty acids,

101

correlations between biohydrogenation intermediates are inevitable. Thus, having a technique

102

with differentiation ability is of great importance. Here we evaluated both FT-NIR and

103

ATR/FTIR spectroscopy techniques for the quantification of individual and grouped trans

104

MUFA, CLA and OBCFA in a dataset of 75 milk samples sub-selected from a collection of more

105

than 1000 milk samples to cover the full biological concentration span of these milk fatty acids.

106 107

MATERIALS AND METHODS

108

Sample Selection

109

The sample storage and sample selection methodology was previously described

110

total of 100 milk samples were selected from a sample bank (n = 1033) of six different cow

111

feeding experiments aiming at changing the fatty acid (FA) profile of dairy products. The sample

112

subset was selected using a genetic algorithm applied to cover the naturally occurring range of

113

several milk fatty acids of interest (different odd and branched chain saturated fatty acids and

114

different trans-C18:1 and cis/trans-C18:2 unsaturated isomers) in high, mid, and low

15

. Briefly, a

5 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

16

Page 6 of 44

115

concentration ranges

. The milk fat was extracted using a previously described methodology

116

involving dichloromethane-ethanol 17.

117

Fatty Acid Standards

118

Two trans monounsaturated trans-9 and trans-10 C18:1 fatty acid methyl ester (ME) standards,

119

three branched chain saturated iso C14:0, iso C15:0 & anteiso C15:0 ME, one odd chain

120

saturated C15:0 ME and two triacylglyceride (TAG) standards (trielaidin and tristearin) were

121

purchased from Larodan (Larodan Fine Chemicals, Malmo, Sweden) and Nu-Check (Nu-Check

122

Prep, Inc., MN, U.S.A.).

123

GC Reference Data

124

Quantification of trans fatty acids (TFA) and fatty acid groups using spectral data requires

125

precise Gas Liquid Chromatography (GC) reference data for the construction of mathematical

126

models. Identification and quantification of TFA through GC have been greatly improved with

127

new highly polar, long capillary columns, but direct GC without prior fractionation could show

128

overlapping between different trans-n and cis-n C18:1 positional isomers and trans-n C16:1

129

coelution with specific branched chain saturated and cis-n C16:1 mono-unsaturated FAs, which

130

might result in an underestimation of the total TFA content

131

dependency of the polarity of cyanopropyl phases 18 to mathematically deduce concentrations of

132

overlapping fatty acids using two different temperature programs without prior fractionation. A

133

similar approach was described before

134

for 92 of the 100 selected milk fat samples) were methylated

135

(FAME) were analyzed by GC according to Vlaeminck et al.2 (first temperature program) and by

136

an isothermal (T=180oC) (second) temperature program. A different separation with the second

137

temperature program allowed the quantification of most individual FA 15, 20 (Table 1). Due to an

15, 20, 21

18, 15

. Here, we used the temperature

. After extraction, all samples (sufficient quantity 17

and the fatty acid methyl esters

6 ACS Paragon Plus Environment

Page 7 of 44

Journal of Agricultural and Food Chemistry

138

influence from the cows’ diet the total trans mono-unsaturated fatty acids (trans MUFA) were

139

present in very high concentrations and ranged from 1.5 to 21 g/100g in the milk samples

140

Total Conjugated Linoleic fatty acids (CLA) maximally represented 2.0 g/100g of milk fat

141

(Table 1). Further, insufficient GC resolution prevented reliable separation of individual trans-6

142

C18:1, trans-7 C18:1 and trans-8 C18:1 FA isomers, but allowed accurate determination of the

143

trans-6+7+8 C18:1 sum. The cross-correlation results between the individual fatty acids and the

144

major fatty acid groups indicated a high correlation (R2>0.95) between trans-10 & trans-11

145

C18:1 and trans MUFA, between cis-9, trans-11 C18:2 and CLA, as well as between C15:0 and

146

OBCFA and C15:0 and ODD.

147

Vibrational Spectroscopy Analysis

148

ATR/FTIR Spectroscopy

149

All attenuated total reflectance Fourier transform mid-infrared (ATR/FTIR) spectra were

150

acquired on a Vertex 70 – RAM II Bruker spectrometer (Bruker Analytical, Madison, WI)

151

operating with a Golden Gate

152

internal reflection element was a small, non-temperature-controlled Type IIa diamond prism

153

allowing a sampled diameter of approximately 2.0 mm. The optically dense medium was in

154

contact with two ZnSe focusing lenses, one used to focus the incident infrared radiation and the

155

second one to collect the reflected infrared radiation. The optical bench included an

156

interferometer with a RockSolid configuration, KBr substrate beam splitter and a deuterated

157

triglycerin sulfate (RT-DLaTGS) detector. The OPUS 6.5 software for Windows of Bruker

158

Instruments was used for instrument management, spectra acquisition and OPUS JCAMP to

159

JCAMP-DX file transformation. All milk fat samples were stored in vials selected by CRA-W in

160

previous Fourier transform Raman analysis

TM

23

.

diamond ATR accessory (Specac Ltd., Slough, UK). The

22

. Prior to ATR/FTIR analysis, the milk fat vials

7 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 8 of 44

161

were taken out of storage and the milk fat was kept in a melted state at a minimum of 30 minutes

162

prior to spectra acquisition in a 38±1oC water bath, accompanied by gently shaking by hand to

163

ensure homogenization. One drop of each milk fat sample in liquid form was placed on top of the

164

optical medium and spectra were immediately acquired [room temperature (∼24±2oC)], as well

165

as selected FA standards (∼0.1 g/sample) were analyzed after freezing at -80oC. The spectra

166

were collected against air as a background over the wavenumber range of 4498-500 cm-1 with a

167

resolution of 4 cm-1. For each spectrum, 64 scans were co-added and averaged to obtain a good

168

signal-to noise ratio. A total of 2074 data points were recorded from 500 to 4498 cm-1. Because

169

of very low milk fat quantity in 14 and extraction solvent contamination in 3 of all 92 selected

170

samples 15, ATR/FTIR spectra were available for 75 milk fat samples. For each milk fat sample,

171

ATR/FTIR spectra were acquired in duplicate by taking a new milk fat drop and the average of

172

the two spectra was used for chemometrical analysis.

173

FT-NIR Spectroscopy

174

The Fourier transform near-infrared (NIR) spectra of the 75 milk fat samples were acquired

175

using a Bruker Optics Matrix-F instrument equipped with a thermoelectrically cooled InGaAs

176

diode detector and a Bruker Optics N261 diffuse reflection fiber optic probe (Bruker Analytical,

177

Madison, WI). The microprocessor controlled optical bench included an interferometer with a

178

RockSolid configuration, quartz substrate beam splitter with proprietary coating and an air

179

cooled tungsten NIR excitation source (12V, 20W) with a HeNe laser at a frequency of 633 nm.

180

Milk fat samples and pure FA standards (∼0.1-0.5 g/sample) were taken out of storage, left to

181

rest until reaching room temperature (∼21±2oC) and analyzed in vials selected by CRA-W in

182

previous Fourier transform Raman spectroscopy analysis 22 with PE-caps (Klaus Ziemer GmbH,

183

Mannheim, Germany). The spectra of all milk fat samples were collected against air as a

8 ACS Paragon Plus Environment

Page 9 of 44

Journal of Agricultural and Food Chemistry

184

background with a resolution of 8 cm-1 and a total of 780 data points were recorded over the

185

wavenumber range of 4000 to 10,000 cm-1 (2500-1000 nm). For each milk fat sample, FT-NIR

186

spectra were acquired in duplicate and the average of the two spectra was used for

187

chemometrical analysis.

188

Data Treatment

189

Data incorporation, treatment and partial least squares (PLS) regression were performed

190

similarly to previous methodology 15. The OPUS spectra files were converted to the JCAMP-DX

191

file system, which was imported into The UnscramblerTM v. 9.1 software (CAMO, Trondheim,

192

Norway). The ATR/FTIR spectra were reduced from 2074 to 1241 variables, by removing

193

selected regions, e.g. noisy region between 500-560 cm-1, regions identified as carriers of -OH

194

chemical information (1020-1070 cm-1 and 3050-4000 cm-1), as well as minor absorbance from

195

the background (2000-2450 cm-1). The observed -OH absorption signal did not originate from

196

the milk fat samples, but was attributed to moisture/trace of ethanol solvent used to clean the

197

surface of the ATR crystal between analyses. Although a drying fan was used, it appeared that

198

some ethanol moisture/vapours could be trapped between the milk fat and the crystal surface,

199

which was apparent in the ATR/FTIR spectra of six samples. As a precaution, these regions not

200

related to FA information were removed and all mathematical treatments or investigations were

201

performed with the remaining 1241 variables. The pre-processing methods of choice were a

202

combination of linear baseline correction and Multiplicative Scatter Correction (MSC) or a linear

203

baseline correction and Standard Normal Variate (SNV). Both MSC and SNV compensate for

204

undesired variation in the scattering intensity caused by various multiplicative effects, thus the

205

choice of using either pre-treatment would be subjective 24. Standard partial least squares (PLS)

206

regression was carried out in combination with the uncertainty testing “jackknife” procedure

9 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 10 of 44

207

from The UnscramblerTM program, which is used to select significant variables that correlated

208

with the measured parameter 24, 25, 26. Further, PLS regression was performed after ATR/FTIR or

209

FT-NIR spectra transformation according to three different strategies- after MSC or SNV pre-

210

processing only, MSC/SNV followed by 1st Savitzky-Golay derivative or MSC/SNV followed by

211

2nd Savitzky-Golay derivative (using a 1,1,2 derivative treatment, which indicates a smoothing

212

factor with the number of data points in a running average on the left, the secondary smoothing

213

on the right, and the polynomial order of the derivative, respectively). A total of six prediction

214

models for each of the individual trans isolated and conjugated fatty acids (trans-4 C18:1, trans-

215

5 C18:1, trans-6+7+8 C18:1, trans-9 C18:1, trans-10 C18:1, trans-11 C18:1, trans-12 C18:1,

216

trans-15 C18:1, cis-9, trans-11 C18:2 and trans-10, cis-12 C18:2), both trans MUFA and CLA

217

FA groups, as well as for each of the individual OBCFA (iso C13:0, anteiso C13:0, iso C14:0,

218

iso C15:0, anteiso C15:0, C15:0, iso C16:0, iso C17:0,anteiso C17:0 and C17:0), as well as for

219

ODD (sum of C5:0, C7:0, C9:0, C11:0, C13:0, C15:0, C17:0, C19:0, C21:0, C23:0), ISO (sum of

220

iso C13:0, iso C14:0, iso C15:0, iso C16:0, iso C17:0 and iso C18:0), ANTE (sum of anteiso

221

C13:0, anteiso C15:0 and anteiso C17:0), BRANCHED (sum of ISO and ANTE) and OBCFA

222

(sum of ODD and BRANCHED) FA groups were built when using either ATR/FTIR or FT-NIR

223

spectra of all 75 milk samples. In addition, to better evaluate the prediction performance for

224

minor fatty acids in their most common concentration range (below 1.0 g/100g total fat) and to

225

avoid overfitting based on discontinuous dataset distribution, milk fat samples with high amounts

226

of specific individual fatty acids (such as trans-10 C18:1, trans-11 C18:1 and trans-12 C18:1)

227

were removed and new PLS models were constructed using only the remainder. For all PLS

228

models validation was carried out using systematic cross-validation with 3 folds and 25 units per

229

segment. The optimal number of PLS factors used for the regression was determined from the

10 ACS Paragon Plus Environment

Page 11 of 44

Journal of Agricultural and Food Chemistry

230

minimum residual validation variance. The cross-validation coefficient of determination (Rcv2),

231

the root mean square error of cross-validation (RMSECV, g/100g FAME), the number of PLS

232

factors, intercept, slope and bias parameters of the best prediction models for individual and

233

grouped FA are presented in Tables 2 and 3 for ATR/FTIR and FT-NIR, respectively. All other

234

results are available as Supporting Material.

235 236

RESULTS AND DISCUSSION

237

While manufacturers are implementing different dietary strategies for the increase in the amounts

238

of beneficial fats

239

prevention of productivity and quality drops. The latter requires fast and cost effective analytical

240

methods, which could determine fatty acids of interest ex situ or desirably in situ. For this

241

purpose, attenuated total reflectance Fourier transform infrared (ATR/FTIR) and Fourier

242

transform near-infrared (FT-NIR) spectroscopy techniques were considered. Looking at the

243

number of scientific publications it becomes evident that FT-NIR has a predominant impact in

244

the food industry. However, ATR/FTIR would seem a more logical choice, since it possesses

245

enhanced spectral absorbance features, such as from specific unsaturated C=C bonds, and has

246

already been established as an official method for the quantification of total trans fat content in

247

foods and dietary supplements

248

several fatty acids with similar structures present in the complex milk fat matrix. Thus, the

249

evaluation of both most widely used spectrophotometrical methods for their ability to produce

250

fatty acid specific signals and in turn to determine various fatty acids present in different high,

251

mid and low concentration combinations is necessary.

252

Fatty Acid Standards

27

, following the health status of animals is always of importance for quick

11

. Here, the interest is in the simultaneous determination of

11 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

253

Page 12 of 44

trans MUFA standards

254

The general differences between the FT-NIR spectra of pure triolein and trielaidin standards in

255

the vicinity of 4500-4700 cm-1 (2200-2130 nm) and 5600-5900 cm-1 (1785-1695 nm) (attributed

256

to the differences in signals from the cis and trans C=C bonds) were previously reported 14. Here

257

two trans monounsaturated fatty acid methyl ester standards, as well as trielaidin and tristearin

258

pure standards were analyzed using FT-NIR and ATR/FTIR spectroscopy techniques. The 2nd

259

Savitzky-Golay derivative FT-NIR spectra of pure tristearin and trielaidin standards were

260

compared in order to establish specific contributions from the isolated trans C=C bond (Figure

261

1). A striking difference between the spectra of both TAG was in the vicinity of 4705 cm-1,

262

where trielaidin and not tristearin showed increased absorption. The latter is a combination

263

region, which might be specific to absorbance from C-H vibrations in isolated trans C=C bonds.

264

Further, analyses of the pure standards with ATR/FTIR spectra confirmed an enhanced

265

absorbance feature from non-conjugated trans C=C bonds in the vicinity of 966 cm-1 (from the

266

C-H out of a plane deformation band), but no differences in the peak position or absorbance

267

intensity were evident between the two trans-9 and trans-11 methyl ester positional isomers

268

(Figure 2). Other bands around 1245 cm-1 and 1000-1130 cm-1 showed clear absorbance

269

intensity and shift perturbations, which might be isomer specific. In addition, the shape of the

270

966 cm-1 band was very broad when liquid spectra of the standards was considered, thus partial

271

overlap with previously reported absorption bands from conjugated cis and trans C=C bonds in

272

cis-9, trans-11 and trans-10, cis-12 CLA is suspected

273

determinations of CLA isomers in low quantities, especially in the presence of high trans MUFA

274

content (corresponding to an increase in band broadness and intensity) 29.

28

. These features might complicate

12 ACS Paragon Plus Environment

Page 13 of 44

Journal of Agricultural and Food Chemistry

275

Moreover, all standards showed drastic differences between the ATR/FTIR absorption spectra

276

obtained at liquid or at solid physical phases. The crystalline spectra of the trielaidin indicated

277

the appearance of new peaks, peak shifts and absorption increase in specific ATR/FTIR

278

wavebands, which are believed to contain extensive structural and vibrational information

279

(Supporting Materials, Figure SM1). The most striking feature (believed to be of most

280

importance) was the concurrent absorbance increase and peak shift of the isolated trans C=C

281

waveband in the vicinity of 966 cm-1 (shifted to 960 cm-1). Moreover, a closer examination of the

282

latter revealed narrower band shape compared to the liquid phase band. Detailed investigation is

283

required for origin identification of all newly reported vibrational features (Supporting

284

Materials, Table SM1).

285

OBCFA standards

286

Limited information is available for specific FT-NIR and ATR/FTIR absorbance features related

287

to branched (ISO and ANTE) and odd numbered (ODD) saturated fatty acids. Here pure iso

288

C14:0, iso C15:0, anteiso C15:0 and C15:0 fatty acid methyl ester standards were analyzed at

289

room temperature. Overall the ATR/FTIR spectra of the iso C15:0, anteiso C15:0 and C15:0

290

standards were similar, except the appearance of two new peaks in the 700-920 cm-1 region and

291

shifting of bands with absorbance intensity increases in the 1330-1400 and 1475-1435 cm-1

292

regions (Figure 3). Unique signals related to iso and anteiso branching were observed in the

293

vicinity of 920, 1366 & 1466 cm-1 and 770, 1375 & 1462 cm-1, respectively (Figure 3).

294

Literature in regards to the origin of the vibrational modes in the lower wavenumber region

295

(below 920 cm-1) for branched alkanes is scarce, however they are believed to be caused by the

296

FA chain C-C skeletal rocking vibrations [the band position difference between anteiso C15:0

297

and iso C15:0 (770 and 920 cm-1, respectively) might be caused by the difference in the position

13 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 14 of 44

298

of the extra CH3 methyl group at the end of the FA chain, which in turn influences the mobility

299

of the FA chain]

300

1365 cm-1) is known to be associated with the C-H symmetric bending vibrations in -CH3

301

The appearance of this more intense absorbance band in the ATR/FTIR spectra of branched, but

302

not in straight chain standards is due to the presence of an extra methyl group at the end of the

303

fatty acid chain - within branched methyl ester standards, the iso and anteiso configuration

304

position of this extra methyl group causes differences in absorbance intensity and band position

305

at 1366 cm-1 and 1375 cm-1 for iso C15:0 and anteiso C15:0, respectively (Figure 3). The 1475-

306

1435 cm-1 absorption region is related to the asymmetrical deformation vibration of -CH3 in

307

branched alkanes 30, 31. Further, no shifts in these iso bands between the spectra of iso C14:0 and

308

iso C15:0 were evident and only minor differences between the spectra of both fatty acids in C-H

309

(higher in iso C15:0) and C=O (higher in iso C14:0) vibrations due to differences in the fatty

310

acid chain length were present (Supporting Materials, Figure SM2). These results might help

311

in differentiation between anteiso, iso and straight chain fatty acids, but specificity for saturated

312

fatty acids of similar nature within the ISO, ANTE or ODD groups apparently seems low.

313 314 315 316

30, 31

. Absorbance in the vicinity of 1350-1395 cm-1 (often in the vicinity of 30, 31

.

Fatty Acid Determinations Most, trans isolated (trans MUFA) and conjugated (CLA) unsaturated, as well as odd and

317

branched chain saturated (OBCFA) fatty acids (FA) of interest are present in minor

318

concentrations (lower than 1.0 g/100g of total fat) in bovine milk fat. Hence their quantification

319

using a spectroscopy technique is a challenge. This is mainly due to the property of

320

spectrophotometrical

321

physical/molecular characteristics of all FA present in the sample. Thus FA with similar

322

chemical structure characteristics, such as the trans isolated and conjugated FA of interest, could

methods

producing spectral

information

based

strictly on

the

14 ACS Paragon Plus Environment

Page 15 of 44

Journal of Agricultural and Food Chemistry

323

have a significant overlap of spectral bands issued from the wide range of vibrations types from

324

the different fatty acid isomers. These specific absorbances are dependent on the types and

325

concentrations in which the corresponding fatty acids are present. For instance, previous research

326

indicated that both trans-10 C18:1 and cis-9, trans-11 C18:2 CLA are important biomarkers,

327

which might reflect the viability status of the ruminal microbial population. The latter is

328

influenced by drop in ruminal pH values, which in turn could drive the animal into acidotic

329

status 3. For this, it would be important to possess a spectroscopy methodology, which is able to

330

differentiate between both types of fatty acid groups (trans MUFA versus CLA) when they are

331

present in low, as well as in higher than “normal” quantities. Typical average concentrations of

332

CLA in milk fat are 0.3-0.6 g/ 100g with cows fed a total mixed diet, but the levels of CLA in

333

milk can vary widely among herds 32, 33 and could reach an average of 1.1 g/100g of CLA in milk

334

fat with cows consuming pasture only 34. Also, it has been shown that the concentrations of CLA

335

in milk are affected by different dietary fatty acid sources and could drastically vary between

336

individual animals, reaching a maximum of 5.17 g / 100 g in milk fat with cows consuming a

337

sunflower oil based diet

338

trans MUFA (~21 g/100g) and high CLA content (~ 2.5 g/100g), thus we attempted to evaluate

339

the prediction performance within those naturally occurring fatty acid concentration ranges. In

340

such cases, it has already been established that in the presence of at least 1.0 g/100g conjugated

341

unsaturation in total fat, the official trans FA determination method is not applicable anymore 9.

342

However, that does not necessarily mean that quantification is impossible, but rather not

343

comparable to officially applied methods. Thus, determinations with lower performance could

344

still have implications for disease screening or monitoring.

345

PLS accuracy

35

. The current dataset contained both milk fat samples with very high

15 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 16 of 44

346

In general, FT-NIR spectroscopy with partial least squares (PLS) regression resulted in complex

347

prediction models (high number of PLS factors), performed poorly in the determination of

348

individual fatty acids and the ATR/FTIR technique was superior (higher cross-validation

349

coefficient of determination, Rcv2, and lower root mean square error of cross-validation,

350

RMSECV, g/100g FAME).

351

trans MUFA determinations

352

The ATR/FTIR spectroscopy method predicted well (Rcv2≥0.90) all major individual trans

353

monounsaturated fatty acids (trans MUFA) as well as the trans MUFA group in high, mid and

354

low concentration (Table 2). The regression coefficients plots for all models were similar and

355

indicated variables in the 960-970 cm-1 band as main predictors. While derivative transformation

356

did not improve the high concentration range models, the best low trans MUFA model (less than

357

4.6 g/100g in total fat) was with 2nd Savitzky-Golay derivative spectra transformation. No

358

variables in the 930-950 and 980-990 cm-1 wavenumber ranges (associated with absorbance from

359

conjugated C=C) were picked as contributors (Figure 4). It seems that the milk fat matrix does

360

not cause major difficulties for prediction of trans MUFA fatty acids in this concentration range,

361

even in the presence of CLA.

362

Further, satisfactory performance was evident for some of the individual fatty acids in minor

363

concentrations (lower than 1.0 g/100g in total fat), such as trans-11 C18:1 (Rcv2=0.795,

364

RMSECV=0.097 with 5 PLS factors) and trans-12 C18:1(Rcv2=0.814, RMSECV=0.066 with 1

365

PLS factors). While trans-11 C18:1 was mainly determined from a specific signal in the vicinity

366

960-970 cm-1 (Figure 5) without the need for derivative transformation, the trans-12 C18:1

367

prediction was due to a high reference data cross-correlation with the total trans MUFA group.

368

Except for trans-9 C18:1 and the trans-6-7-8 C18:1 sum (Table 2), all other ATR/FTIR models

16 ACS Paragon Plus Environment

Page 17 of 44

Journal of Agricultural and Food Chemistry

369

for fatty acids in very minor concentrations (lower than 0.50 g/100g in total fat) resulted in

370

Rcv2