Effects of Cultivar and Nitrogen Nutrition on the ... - ACS Publications

MS/MS and GC-FID. 440. Figure S2 Starch damage analysis of the six cultivars in three nitrogen conditions. 441. Figure S3 Relevant means on the log sc...
0 downloads 15 Views 1MB Size
Subscriber access provided by TUFTS UNIV

Article

Effects of Cultivar and Nitrogen Nutrition on the Lipid Composition of Wheat Flour Byoung Min, Irene Gonzalez-Thuillier, Stephen J. Powers, Peter Wilde, Peter R. Shewry, and Richard P. Haslam J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.7b01437 • Publication Date (Web): 14 Jun 2017 Downloaded from http://pubs.acs.org on June 20, 2017

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 31

Journal of Agricultural and Food Chemistry

1

Effects of Cultivar and Nitrogen Nutrition on the Lipid Composition of Wheat

2

Flour

3

Byoung Min1, Irene González-Thuillier1, Stephen J. Powers2, Peter Wilde3, Peter R

4

Shewry1* and Richard P Haslam1

5

6

1

Plant Science, Rothamsted Research, Harpenden, AL5 2JQ, United Kingdom

2

Computational and Analytical Sciences, Rothamsted Research, Harpenden, AL5

7

2JQ, United Kingdom

8

4

9

NR4 7UA, United Kingdom

10

Quadram Institute Bioscience, Institute of Food Research, Norwich Research Park,

*Author for correspondence: [email protected]

11

12

13

14

15

16

17

1 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

18

ABSTRACT: Despite being minor components of flour, wheat (Triticum aestivum L.)

19

lipids contribute to grain processing. They are particularly important for bread making

20

where they adsorb to the surface of gas bubbles formed during the proving stage of

21

bread making, stabilizing the gas cells and improving gas retention within the dough.

22

This contributes to the volume and texture of the loaf. However, little is understood

23

about how their amount, composition and properties vary in response to genotype

24

(G), environment (E) or G x E interactions. Six wheat lines were therefore grown at

25

three levels of nitrogen supply at Rothamsted Research and 48 lipid species across

26

six lipid classes were identified and quantified in white flour using electrospray

27

ionization-tandem mass spectrometry (ESI-MS/MS). This showed clear differences

28

in the contents and compositions of lipids between cultivar, as well as effects of

29

nitrogen fertilization, which would be expected to have impacts on the processing

30

properties of the samples.

31

32

KEY WORDS: wheat grain, lipids, bread making, lipidomics, genotype, environment

33

34

35

2 ACS Paragon Plus Environment

Page 2 of 31

Page 3 of 31

Journal of Agricultural and Food Chemistry

36

INTRODUCTION

37

Wheat (Triticum aestivum L.) flour is widely used for food production, including

38

breads, other baked products, pasta, noodles and as an ingredient in many

39

processed goods such as sauces and processed meats. White flour is derived from

40

the starchy endosperm storage tissue of the grain and comprises mainly starch

41

(about 75-85%) and protein (about 10%). These two major components largely

42

determine the processing properties, with the gluten proteins determining the visco-

43

elastic properties of dough that underpin many of its uses, including bread making.

44

However, wheat flour also contains a range of other components including about 2.0

45

to 2.5% lipids.1-3

46

Despite their low concentrations, lipids have an impact on the quality of wheat flour

47

for bread making, affecting the volume and crumb structure of the loaf.4 These

48

effects are not completely understood, but they are thought to include indirect effects

49

by binding to and plasticising the gluten network, and direct effects by stabilizing the

50

structure of gas cells which are formed during dough mixing.5-6 These gas cells

51

entrap the carbon dioxide released by yeast during fermentation leading to the

52

expansion of dough mass, increasing its volume and determining the crumb

53

structure. The importance of lipids in bread making has also been shown by recent

54

studies in which specific lipases have been used to optimise lipid composition and

55

improve processing quality.7-8 Therefore, understanding the processes that define

56

grain lipid content and composition is critical to the development of wheat varieties

57

with optimized quality for bread making.

58

Lipids are present throughout the wheat grain and include phospholipids (PL) which

59

are structural components of membranes in all tissues, and triacylglycerols (TAG)

3 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

60

which are stored in oil bodies in the aleurone layer and embryo. PLs and TAGs are

61

acyl lipids in that they contain fatty acids esterified to a glycerol backbone, while fatty

62

acids also occur in their free form (FFA). Wheat flour lipids are often divided into

63

starch and non-starch lipids. Starch lipids are intrinsic components of starch granules

64

and as such they can only be extracted when the starch granules are broken down

65

or damaged. They do not play a significant role in determining the properties of

66

dough, and therefore do not affect the final quality of the bread.9 By contrast, non-

67

starch lipids have significant effects on dough properties and product quality.6,9,10,11

68

Lipidomic analysis has previous been used to show that lipid remodelling takes place

69

in vegetative tissues of resilient wheat cultivars when exposed to temperature

70

extremes.12-13 However, although two studies have used detailed lipidomic

71

approaches to determine the compositions of lipids in wheat flour,14-15 currently little

72

is known about the extent to which lipid composition varies between genotypes, or is

73

affected by environmental factors (notably temperature and water availability), or by

74

crop nutrition, in particular by the application of nitrogen fertilizer which varies greatly

75

between production systems. This information is crucial to enable grain processors

76

to predict and exploit differences in lipid composition to improve the quality and

77

consistency of products. We have therefore applied a lipidomic approach to

78

determine the contents and compositions of the major lipid classes (PLs, TAGs,

79

diacylglycerols (DAGs) and FFAs) in five UK hard bread making wheat cultivars and

80

one soft feed wheat, comparing material grown at three levels of nitrogen

81

fertilization: 100 kg/Ha (low input), 200 kg/Ha (typical of UK intensive production) and

82

350 kg/Ha (luxury supply).

83

4 ACS Paragon Plus Environment

Page 4 of 31

Page 5 of 31

Journal of Agricultural and Food Chemistry

84

MATERIALS AND METHODS

85

Samples. Six UK wheat cultivars were grown on the experimental farm at

86

Rothamsted Research, Harpenden, UK in 2013 as part of the Wheat Genetic

87

Improvement Network (WGIN) nitrogen use efficiency trial.16 These were five modern

88

hard bread making cultivars (Avalon, Cadenza, Crusoe, Hereward and Malacca) and

89

one soft feed cultivar (Istabraq).

90

Cultivars were grown in triplicate 9 x 3 m plots at three levels of nitrogen fertilization,

91

100, 200 and 350 kg/Ha using ammonium nitrate. 200 kg/Ha (200N) is typical for

92

intensive wheat production in the UK, while 100 kg/Ha (100N) is more typical of less

93

intensive systems used in many other countries. The highest application level of 350

94

kg/Ha (350N) was above those typically used in crop production but was included to

95

explore the effects of unusually high nitrogen availability.

96

Plots were randomized within main plots of the nitrogen treatments. Grain was

97

pooled from the three plots of each combination to obtain representative samples. 50

98

g samples of grain per combination were then milled using a Chopin CD1 mill, which

99

complies with the NF EN ISO 27971 standard, after being conditioned to a 14%

100

moisture content. Straight grade flour was obtained for each sample, after

101

undergoing a break and reduction stage and sifting through a centrifugal sifter. Flour

102

samples were stored at -20°C. Starch damage was measured using the Megazyme

103

starch damage kit (Megazyme, Bray, Ireland).

104

105

Lipid Extraction. Non-starch lipids were extracted from 5 replicate samples of each

106

flour as described by Finnie et al14 with some modifications. 150 mg of flour in a

107

glass tube was heated in boiling water for 12 mins to inactivate lipid-hydrolysing 5 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 6 of 31

108

enzymes. Samples were then extracted sequentially with petroleum ether (PEt),

109

water-saturated butan-1-ol (WSB), and propan-2-ol/water (90:10) (IW). The PEt and

110

WSB extracts were washed by shaking with 0.88% KCl in a 1:1 ratio for 2 mins at

111

18000 rpm. The supernatants from the PEt and WSB extractions were combined,

112

evaporated under nitrogen at 40 °C, re-suspended in an equal volume of chloroform

113

and washed with 0.88% KCl. The lower phase was retained, filtered through a

114

0.45µm Millex-FH filter (Merck Millipore, Germany), dried under a stream of nitrogen

115

gas and re-suspended in 2mL of chloroform. These were then flushed with nitrogen

116

and stored at -20 °C.

117

Quantitative Lipid Analysis. Electrospray ionisation tandem triple-quadrupole mass

118

spectrometry (4000 QTRAP; SCIEX; ESI-MS/MS) was used to quantify the major

119

types

120

phosphatidylethaloamine (PE), lysophosphatidylcholine (LPC), and the neutral lipids

121

(NLs) FFA, DAG and TAG. Monoacylglycerols (MAG) were not analysed due to their

122

low levels in flour. Internal standards for polar lipids were obtained from Avanti

123

(Alabaster, AL, USA) and were incorporated as 0.857 nmol of 13:0-LPC, 0.086 nmol

124

of di18:0-PI, and 0.080 nmol of di14:0-PE, 0.800nmol of di18:0 PI, and 0.080 nmol of

125

di14:0-PG, dissolved in chloroform. 25 µL of sample dissolved in chloroform was

126

combined with the standard and chloroform: methanol:400mM ammonium acetate

127

(300:665:3.5 v/v) to make a final volume of 1ml. The lipid extracts were infused at

128

15uL/min with an autosampler (HTS-xt PAL, CTC-PAL Analytics AG, Switzerland).

129

Data acquisition and acyl group identification were as described by Gonzalez-

130

Thuillier et al15 with modifications

131

FFAs were analysed by combining 25 µl of sample, 0.607 nmol of 15:0-FFA (Sigma-

132

Aldrich,

of

St

lipid

Louis,

in

wheat

MO,

USA)

flour:

the

and

PLs

phosphatidylcholine

propan-2-ol/methanol/50mM 6

ACS Paragon Plus Environment

(PC),

ammonium

Page 7 of 31

Journal of Agricultural and Food Chemistry

133

acetate/dichloromethane (4:3:2:1; v/v) to a final volume of 1ml, quantified by the Q1

134

ESI-MS negative mode. DAG and TAG contents were quantified using the ESI-

135

MS/MS method described by Li et al

136

sample and 0.857 nmol of tri15:0-TAG (Nu-Chek Prep, Elysian, MN, USA) was

137

combined with chloroform/methanol/300 mM ammonium acetate (24:24:1.75: v/v),

138

and for DAG, 25uL of sample and 0.857 nmol of 18:0-20:4-DAG (Sigma-Aldrich)

139

were

140

acetate:dichloromethane (4:3:2:1; v/v), to the final volume of 1mL for direct infusion

141

to the mass spectrometer. TAG and DAG were detected as [M+NH4]

142

series of different neutral loss scans, targeting the losses of fatty acids. The scans as

143

well as the parameters used for the three NLs are shown in Table S1 of the

144

Supporting Information. The data were processed using the Lipid View Software

145

(Sciex, Framingham, MA, USA) where isotope corrections were applied. The peak

146

area for each lipid was normalized to the internal standard and further normalized to

147

the weight of the initial sample.

148

The ionization efficiency varies among acyl glycerol species with different fatty acyl

149

groups, and response factors for individual species were not determined in this

150

study. Consequently, the values are not directly proportional to the TAG/DAG

151

contents of each species. However, the approach does provide a valid comparison

152

of TAG/DAG species between samples.

153

Determination of Total Fatty Acid Methyl Esters (FAMEs). Methyl ester

154

derivatives of total fatty acids extracted were analysed by GC-FID (flame ionization

155

detection) using a gas chromatography-flame ionization detector (GC-FID, Agilent

156

6890, Palo Alto, CA, USA) with an AT-225 capillary column of fused silica (30 m

157

length, 0.25 mm id, 0.20 µm film thickness). The oven temperature cycle was set

combined

with

17

with some modifications. For TAG, 15 µl of

propan-2-ol:

methanol:50

7 ACS Paragon Plus Environment

mM

ammonium

+

ions by a

Journal of Agricultural and Food Chemistry

158

with a start temperature of 50 °C which was held for 1 min to allow vaporized

159

samples and solvent (hexane) to condense at the front of the column. The oven

160

temperature was then raised rapidly to 190 °C at a rate of 40 °C/min followed by a

161

slower increase to 220 °C, which was held for 1 min, giving a total run time of 25 min

162

and 50 s per sample. Hydrogen was used as the carrier gas. FAMEs were identified

163

by comparison with known standards (37 FAMEs, Sigma, St Louis MO) and

164

confirmed by GC-mass spectrometry (Agilent 6890N, Palo Alto, CA, USA).

165

Validation of Quantification. Phospholipids are more precisely determined by MS

166

than FFA and TAG, which are overestimated and underestimated, respectively.

167

Correction factors were therefore determined and applied to the data for FFA and

168

TAG. Equal volumes of lipid extractions from 28 flour samples were taken and

169

quantified using two systems: ESI-MS/MS as described above, and by thin layer

170

chromatography – gas chromatography (TLC-GC). TLC-GC was carried out by

171

adding 50 µg of 45:0 TAG/15:0 FFA standard to each sample, and separating out the

172

lipid extracts on a silica gel TLC (plate thickness = 0.25 mm) using the solvent

173

hexane / diethyl ether/acetic acid (150:50:2 by volume). The individual lipid classes

174

were identified under UV light after spraying with primuline (0.05% w/v in

175

acetone/water, 80:20 v/v), and the TAG and FFA were scraped from the plate and

176

directly methylated for FAME analysis.

177

The quantity of each fatty acid was calculated in comparison to the internal standard

178

and normalized for sample weight (g of flour). The values (nmol/g flour) from the two

179

procedures were compared in a scatter plot. The data for FFA were directly

180

compared, whereas for TAG, the data for molecular species determined by ESI-

181

MS/MS systems were compared to the sum of FAMEs determined by GC-FID. The

182

sum of FAMEs provides a valid representation of the quantities of TAG comparable 8 ACS Paragon Plus Environment

Page 8 of 31

Page 9 of 31

Journal of Agricultural and Food Chemistry

183

to that measured by ESI-MS/MS. (see Figure S1A and S1B of the Supporting

184

Information)

185

The GenStat statistics package (2016, 18th edition © VSN International Ltd, Hemel

186

Hempstead, UK) was used to fit an asymptotic exponential model using the method

187

of non-linear least squares to estimate two parameters (the asymptote and the

188

exponential rate) with standard error for both GC-FID (FAMES) FFA and TAG in

189

terms of the ESI-MS/MS (QTRAP) equivalents. The SigmaPlot package (13th edition,

190

Systat Software Inc, San Jose, California, USA) was used to produce the picture of

191

the fitted curve on the scatter plot (Figure S1A, B). The model allowed prediction of

192

values of FFA and TAG for one method given values for the other.

193

Statistical Analysis. The dataset included a small number of missing values.

194

Estimates of these values were therefore initially made based on the lipid means

195

within each cultivar by nitrogen combination using the values for units that had no

196

missing values for any lipid. The estimates were then recalculated as the fitted

197

values from the multiple regression of that lipid on all the other lipids and the lipid

198

means recalculated. If any of the means differed from the previous mean by more

199

than a proposed tolerance (the initial standard error divided by 1000) the process

200

was repeated.18-19

201

This allowed the calculation of the percentage data per lipid for analysis, using total

202

quantity per sample lipid profile without the issue of missing values biasing the

203

calculation of the percentages.

204

Following this, Canonical Variates Analysis (CVA) was performed to give a linear

205

discrimination between the cultivar by N combinations and to provide a low

206

dimensional representation of the differences between them, so that (in two

9 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

207

dimensional plots) biologically meaningful differences could be tentatively assigned

208

using non-overlapping 95% confidence circles around the means of CV scores per

209

treatment combination, assuming of a multivariate Normal distribution for the data,

210

although the differences cannot be assigned statistical significance because the

211

replicates were technical rather than biological. However, with five technical,

212

replicates the total number of units, 90, minus the number of treatment combinations,

213

18, is 72 which is greater than the number of variables (lipids), 48, required for a

214

statistically valid CVA. The distributional assumption was acceptable having applied

215

a natural log (to base e) transformation to the data.

216

The magnitudes of the CV loadings on the variables (quantified lipids) indicated the

217

relative importance of the original lipids in the discrimination observed (see

218

Supplementary Data S6 for loadings). This allowed only the most relevant results of

219

univariate analysis of variance (ANOVA) to be interpreted, testing the main effects

220

and interactions between the factors of cultivar and nitrogen (F-tests) for each lipid

221

on the log scale, but again with the caveat of the replication being technical rather

222

than biological. The GenStat statistics package (18th Edition, © VSN International

223

Ltd, Hemel Hempstead, UK) was used for this analysis.

224

10 ACS Paragon Plus Environment

Page 10 of 31

Page 11 of 31

Journal of Agricultural and Food Chemistry

225

RESULTS AND DISCUSSION

226

In order to determine the effects of genotype and nutrition on grain lipid composition,

227

six UK wheat cultivars, comprising five hard bread making cultivars (Avalon,

228

Cadenza, Crusoe, Hereward and Malacca) and the soft cultivar Istabraq (used for

229

making biscuits and livestock feed), were analysed from field trials grown with three

230

levels of nitrogen fertilization.

231

Mature grain from the replicated field trial was milled using a Chopin CD1 mill to

232

produce white flour. This mill gives low levels of starch damage (determined as

233

ranging from 4 to 6% in the present study, Supporting Information Figure S2) which

234

would minimise the extraction of intrinsic starch granule lipids which do not affect

235

bread making quality.9 ESI-MS/MS

236

classes of NLs (DAG, TAG and FFA) and PL (PC, LPC and PE). The notation for

237

FFA species reflects the number of carbons and double bonds (for example, oleic

238

acid is C18:1) while the notation for species of neutral lipids and PL reflects the total

239

numbers of double carbons and double bonds for the acyl chains.

15

was used to identify and quantify the major

240

241

Relative abundance of lipids. An ESI-MS/MS approach was used to compare the

242

amount of each lipid molecular species in the different grain samples. The lipid

243

molecular species were identified by precursor or neutral loss scanning, and the

244

lipids in each head group class were quantified in comparison with internal standards

245

of that class. Acquisition of MS data requires long periods of sample infusion.

246

Therefore, a quality-controlled approach was employed to remove any instrument or

247

analytical variation from the data acquisition process. For example, the intensity of

248

lipid species in samples was normalized to internal standards run with the biological

11 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 12 of 31

249

extractions and separately. Furthermore, each sample infusion was replicated and

250

the order of survey scans (Supporting Information Table S1) was changed, for

251

example, the position of the PC 184 m/z head group scan was varied between the

252

first and last positions within the analytical run. Lipids account for about 1.4-2.0% of

253

the dry weight of white wheat flour20 with FFA, TAG, DAG, PC, LPC and PE

254

accounting for over 90% of the total (15 and authors’ unpublished results). Minor lipid

255

components, comprising galactolipids and minor phospholipids, were not analysed.

256

The lipid data are presented as a mol% of the total lipids, obtained by normalizing

257

the peak area of each lipid to the internal standard which was further normalized to

258

the weight of the initial flour sample.

259

Of the lipids extracted from white flour in this study, FFA were generally the most

260

abundant class, accounting for an average of 31% of the lipids analysed. These

261

were followed by LPC, TAG and DAG, which accounted for averages of 27%, 23%

262

and 11% of the total lipids, respectively (Figure 1). The observed percentages agree

263

with Chung et al

264

used samples grown at about 200 kgN/Ha), although in the present study TAGs

265

were more abundant than FFA or LPC for Cadenza and Crusoe grown at 350N.

266

Nevertheless, other authors have reported similar results, with TAGs being more

267

abundant than FFAs in four wheat cultivars20.

268

The proportions of individual species within lipid classes also varied. Among the ten

269

FFA species, C16:0 and C18:2 comprised about 22% and 48% of the total FFA,

270

respectively, while LPC species LPC 16:0 and LPC 18:2 represented 86% to 88% of

271

total LPC. Among the ten species of PC, PC 34:2 and PC 36:4 comprised 37% and

272

32%, of the total, respectively. TAGs comprised 13 species with TAGs 52:4, 54:5

273

and 54:6 comprising about 20%, 18% and 24% of the total, respectively.

21

and the more recent study by Gonzalez-Thuillier et al

12 ACS Paragon Plus Environment

15

(which

Page 13 of 31

Journal of Agricultural and Food Chemistry

274

275

Effects of Genotype and Nitrogen on Lipid Classes. In order to test the main

276

effects and interactions between the six cultivars and three nitrogen conditions, firstly

277

the total amounts of six lipid classes (LPC, PE, PC, FFA, TAG and DAG), together

278

with total PL (LPC, PC, PE) and total NL (FFA, TAG and DAG) were analysed by

279

Analysis of Variance (ANOVA) (Table 1).

280

No main effects of cultivar and nitrogen, or of cultivar x nitrogen interactions, were

281

observed for total LPC, total DAG or total PL. However, cultivar strongly affected the

282

levels of PC, PE, TAG and total NL, while nitrogen treatment had the greatest effects

283

on total NL, FFA, TAG and PE. Although interactions between cultivar and nitrogen

284

were less significant than the main effects of cultivar and nitrogen alone, they did

285

significantly (p < 0.05, F-test) affect the levels of PE, FFA and total NL (Table 1).

286

Although cultivar had the greatest effects on PC and PE, there was no effect on LPC

287

and hence on total PL. By contrast, total NL was significantly (p < 0.05, F-test)

288

affected by both factors and the interactions between them. The effects on NL

289

probably resulted mainly from effects on TAG species, which were greatly affected

290

by both cultivar and nitrogen (see Table 1). However, there was no effect of cultivar x

291

nitrogen on total TAG. Total FFAs were not affected by cultivar but were significantly

292

(p < 0.05, F-test) affected by nitrogen, and by cultivar x nitrogen.

293

294

Effects of Genotype and Nitrogen on Lipid Species. Differences in PL

295

composition between cultivars were observed with Cadenza having lower

296

proportions of total PC and PE species than the other cultivars (Figure 2). The

297

proportions of individual PE species also varied between cultivars and nitrogen 13 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

298

inputs (Supporting Information Figure S4B). Statistical analyses (p-values from

299

ANOVA F-tests) of the effects of cultivar, nitrogen and cultivar x nitrogen on the

300

individual PL species are shown in Supporting Information Table S2. All nine PC and

301

all three PE species were significantly (p < 0.05, F-test) affected by cultivar, with four

302

PC species and all three PE species also being affected by cultivar x nitrogen. Only

303

one out of five LPC species (LPC 18:0) was affected by cultivar and cultivar x

304

nitrogen.

305

The compositions of NLs differed to a limited extent between cultivars, but greater

306

differences were observed depending on nitrogen supply, particularly for Crusoe

307

grown at N100 and N200 (Figure 2). Statistical analysis of the NL species showed

308

that only five of the twelve TAG species (TAG 50:1, 52:2, 52:5, 54:3 and 56:5) were

309

significantly (p < 0.05, F-test) affected by cultivar and three of these (TAG 50:1, 52:5,

310

56:5) by nitrogen and cultivar x nitrogen (Supporting Information Table S4). Six of

311

the nine DAG species were also affected by either one or both main effects, with

312

only DAG 32:0 being affected by the interaction. Although ten species of FFA were

313

determined, only four of these showed significant (p < 0.05, F-test) effects of

314

nitrogen, with none being affected by cultivar or cultivar x nitrogen.

315

In broad terms, these analyses show that PLs were generally more affected by

316

cultivar and NLs by nitrogen.

317

318

Multivariate Statistical Analysis. Canonical Variates Analysis (CVA) was used to

319

obtain a linear discrimination between the cultivar by nitrogen combinations and

320

allowed a low-dimensional representation of the differences to be made (in two-

321

dimensional plots). Bearing in mind the technical nature of the replication, this

14 ACS Paragon Plus Environment

Page 14 of 31

Page 15 of 31

Journal of Agricultural and Food Chemistry

322

allowed tentative biological differences to be assigned by non-overlapping 95%

323

confidence circles around the means of the CV scores for each treatment

324

combination, assuming a multivariate Normal distribution for the dataset on the log-

325

scale. The magnitudes of the CV loadings on the variables (quantified lipids) indicate

326

the relative importance of the lipids in the discrimination observed. In this way, the

327

effects of cultivar, nitrogen and cultivar x nitrogen on the lipid composition could be

328

identified.

329

The first three CVs accounted for 66.17% of the total variance and the possible

330

discrimination in the dataset. DAG 36:4, DAG 36:3, FFA 16:1 and PC 34:1 were the

331

most important lipids for separation in the first CV dimension, TAG 52:5, FFA 18:2,

332

TAG 54:7, TAG 54:3, FFA 24:0, TAG 50:2, LPC 16:0, DAG 36:4 and DAG 36:5 for

333

CV2 and LPC 16:0, LPC 18:2, DAG 36:3, PC 34:1, TAG 52:5, DAG 36:5, FFA 18:2,

334

FFA 24:0, DAG 32:0, PC 36:2 and FFA 18:1 for CV3. The full data regarding the

335

loading values for the dataset, which determines the “important” lipids, are given in

336

Supporting Information, Table S3.

337

CV1 clearly separates Istabraq at all three nitrogen levels and Avalon grown at 100

338

and 200kg/N/Ha from the other cultivars and nitrogen treatments (Figure 4).

339

Cadenza is separated from the other cultivars in the CV2 dimension, particularly the

340

samples grown at 200 and 350kg/N/Ha. In broad terms, the separations in CV1 and

341

CV2 are determined more by differences between cultivars than between nitrogen

342

treatments. Istabraq is particularly well separated from the other cultivars by CV3

343

with good separation also of Hereward grown at 350N and Crusoe grown at 200N.

344

345

GENERAL DISCUSSION

15 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

346

Comparison of six cultivars grown at three levels of applied nitrogen showed that

347

cultivar and nitrogen had effects on the lipid composition of the white wheat flour,

348

with stronger effects on PC, PE and TAG species than on other lipid species.

349

PLs were most affected by cultivar, in terms of total classes and individual PC and

350

PE species, but little effect on LPC was observed. The proportions of LPC

351

determined here were also higher than in other studies.22,14 This may have resulted

352

from partial hydrolysis of PC by phospholipase A2 during storage of the flour,

353

although we failed to detect increases in PC in white flour stored over a period of 3

354

months (authors’ unpublished results). We therefore consider that it is more likely to

355

have resulted from endogenous enzyme activity in the grain before milling, as LPC is

356

produced from PC during germination.23

357

PLs have been shown to have a positive effect on loaf volume,24 acting as

358

emulsifiers to improve the baking performance of wheat dough.25-26 They have also

359

been shown to interact with gluten proteins, which may contribute to increased loaf

360

volume and improved dough strength.27 The effects of cultivar and nitrogen on PC

361

and PE demonstrated here may therefore have implications for processing quality.

362

The amount of TAG was significantly (p < 0.05, F-test) affected by both cultivar and

363

nitrogen. TAG is hydrolysed to DAG and FFA during germination and this may also

364

occur during flour storage at room temperature. High levels of FFA are undesirable

365

for food production as oxidation may lead to rancidity while NLs have been reported

366

to destabilize gas cells and negatively affect loaf volume.4 Unsaturated FFA in

367

particular have been reported to reduce bread loaf volume.28-29 FFA do not usually

368

accumulate in healthy living tissues and it is therefore likely that the high levels of

369

FFA reported here result largely from TAG degradation, either in the grain or stored

16 ACS Paragon Plus Environment

Page 16 of 31

Page 17 of 31

Journal of Agricultural and Food Chemistry

370

flour. The degradation of TAG into FFA has been widely reported in various

371

systems30 while our own studies have shown increases in FFA content together with

372

decreases in TAG in wheat flour.31

373

Multivariate analysis showed tentative differences between cultivars, with Cadenza

374

and Istabraq being the most distinct, followed by Avalon. The other three cultivars,

375

Crusoe, Hereward and Malacca showed no differences, except for Crusoe and

376

Malacca at N200 and Hereward at N350. In broad terms, the variation in lipid

377

composition between the cultivars was greatest at a nitrogen input of 200kg/N/Ha,

378

which is the typical level of fertilization for intensive wheat production in the UK.

379

The separation of Istabraq is not surprising as it is the only soft wheat in this study.

380

Several studies have shown a correlation between endosperm texture and the

381

amount and/or type of lipid. For example, hard Australian wheat cultivars have higher

382

levels of hexane-extractable free lipids than soft cultivars

383

relationship was reported between hardness and total PL content in British

384

cultivars.33 Cadenza is a hard bread making wheat but differs from the other five

385

cultivars in being a spring type (though it is winter hardy and generally sown in

386

autumn in the UK). Hence, the differences in composition between Cadenza and the

387

other cultivars could relate to their pedigree, as spring and winter breeding programs

388

generally use different sets of germplasm. Other authors have also reported

389

differences between the lipid profiles of cultivars, with Hargin and Morrison

390

comparing four bread wheat varieties and Beleggia et al (2013) four durum wheat

391

cultivars.35

392

We have therefore demonstrated effects of genotype and nitrogen nutrition on the

393

content and composition of PL, TAG and FFA, which could have impacts on dough

17 ACS Paragon Plus Environment

32

while a strong inverse

34

Journal of Agricultural and Food Chemistry

Page 18 of 31

394

properties and bread making performance. However, although this is the most

395

detailed study so far reported, it should be noted that only a single harvest year was

396

studied.

397

environmental factors can be expected, including those demonstrated here related to

398

genotype and agronomy. Consequently, the exploitation of differences in processing

399

quality resulting from lipid composition by wheat breeders and wheat processors

400

poses a challenge.

Hence,

year-to-year

differences

between

grain

samples

due

to

401

402

AUTHOR INFORMATION

403

Corresponding author

404

*Email: [email protected]

405

Funding

406

This work was supported by the BBSRC Crop Improvement Club (BB/J019526/1

407

“The role of lipids in determining gas bubble retention in wheat dough”) Rothamsted

408

Research and the IFR recieve strategic funding from the Biotechnological and

409

Biological Sciences Research Council (BBSRC).

410

Notes

411

The authors declare no competing financial interest.

412

ABBREVIATIONS USED

413

AV, Avalon; CA, Cadenza; CR, Crusoe; CVA, Canonical Variates Analysis; DAG,

414

diacylglycerols; ESI-MS/MS, electrospray ionization tandem triple quadrupole mass

415

spectrometry; FAME, fatty acid methyl ester; FFA, free fatty acid; GxE, genetic x

18 ACS Paragon Plus Environment

Page 19 of 31

Journal of Agricultural and Food Chemistry

416

environmental; HE, Hereward; IS, Istabraq; KCl, potassium chloride; LPC,

417

lysophosphatidylcholine; MA, Malacca; N, nitrogen; NL, neutral lipid; PC;

418

phosphatidylcholine; PE, phosphatidylethanolamine; PEt, petroleum ether; PL,

419

phospholipids; TAG, triacylglycerol; TLC-GC-FID, thin layer chromatography-gas

420

chromatography; WGIN, Wheat Genetic Improvement Network; WSB, water

421

saturated butal-1-ol

422

423

ACKNOWLEDGEMENTS

424

The work reported here was supported by the Crop Improvement Research Club by

425

a PhD studentship to BM as part of the grant BB/J019526/1 “The role of lipids in

426

determining gas bubble retention in wheat dough”). Rothamsted Research and IFR

427

receive strategic funding from the Biotechnology and Biological Sciences Research

428

Council (BBSRC).

429

430

ASSOCIATED CONTENT

431

Supporting Information

432

Table S1 The ESI-MS/MS methods and parameters to identify the molecular species

433

of the neutral lipids FFA, DAG and TAG.

434

Table S2 The p-values for F-tests of cultivar, nitrogen and cultivar x nitrogen

435

interactions on the lipids analysed in this study.

436

Table S3 The latent (loading) vector values on the variables (quantified lipids).

19 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

437

Figure S1A Fitted asymptotic exponential curve for the FFA data comparing ESI-

438

MS/MS and GC-FID.

439

Figure S1B Fitted asymptotic exponential curve for the TAG data comparing ESI-

440

MS/MS and GC-FID.

441

Figure S2 Starch damage analysis of the six cultivars in three nitrogen conditions.

442

Figure S3 Relevant means on the log scale for the WGIN lipid totals data.

443

Figure S4A Lipid composition of selected minor neutral lipids (mol% of total lipids).

444

Figure S4B Comparison of LPC and PE (mol% of total lipids) between the six

445

cultivars.

446

The Supporting Information is available free of charge on the ACS Publications

447

website at DOI

448

449

REFERENCES

450 451

1)

Kim, W. S.; Seib, P. A. Lipids in some commercial wheat flours 1. Cereal

452

Chem. 1993, 70, 367-372.

453

2)

454

Chem. 1998, 75, 826-829.

455

3)

456

sources, interactions, and impact on bread quality. J. Cereal Sci. 2011, 54, 266-279.

457

4)

458

50, 292-302.

Jun, W. J.; Chung, O. K.; Seib, P. A. Lipids in Japanese noodle flours. Cereal

Pareyt, B.; Finnie, S. M.; Putseys, J. A.; Delcour. J. A. Lipids in bread making:

MacRitchie, F.; Gras, P. W. Role of flour lipids in baking. Cereal Chem. 1973,

20 ACS Paragon Plus Environment

Page 20 of 31

Page 21 of 31

Journal of Agricultural and Food Chemistry

459

5)

Köhler, P. Study of the effect of DATEM. 3: Synthesis and characterization of

460

DATEM components. LWT-Food Science and Technology. 2001, 34, 359-366.

461

6)

462

Cereal Chem. 1978, 55, 598-618.

463

7)

464

the role of wheat lipids in bread making. Food Chem. 2014, 156,190-196.

465

8)

466

analysis and functional effects of lipid reaction products. J Agric. Food Chem. 2014,

467

62, 8229-8237.

468

9)

469

wheat starch granules, Starch‐Stärke 1978, 30, 119-125.

470

10)

Morrison, W. R. Wheat lipid composition. Cereal Chem. 1978, 55, 548-558.

471

11)

Carr, N. O.; Daniels N. W. R.; Frazier, P.J. Lipid interactions in breadmaking.

472

Crit. Rev. Food Sci. Nutr. 1992. 31, 237-258.

473

12)

474

lipids during heat stress: I. High day and night temperatures result in major lipid

475

alterations. Plant, Cell Environ. 2016, 39, 787-803.

476

13)

477

Lipids experiencing coordinated metabolism are detected by analysis of lipid co‐

478

occurrence. Plant, Cell Environ. 2016, 39, 608-617.

479

14)

480

polar lipids from wheat whole meal, flour, and starch. Cereal Chem. 2009, 86, 637-

481

645.

482

15)

483

Powers, S. J.; Ward, J. L.; Wilde, P.; Shewry, P. R.; Haslam, R. P. Distribution of

Chung, O. K.; Pomeranz, Y.; Finney, K. F. Wheat flour lipids in breadmaking.

Gerits, L. R.; Pareyt, B.; Delcour, J. A. A lipase based approach for studying

Schaffarczyk, M.; Østdal, H.; Koehler, P. Lipases in wheat breadmaking:

Meredith, P.; Dengate, H. N.; Morrison, W. R. The lipids of various sizes of

Narayanan, S.; Tamura, P. J.; Roth, M. R.; Prasad, P. V.; Welti, R. Wheat leaf

Narayanan, S.; Prasad, P. V.; Welti, R. Wheat leaf lipids during heat stress: II.

Finnie, S. M.; Jeannotte, R.; Faubion, J. M. Quantitative characterization of

Gonzalez-Thuillier, I.; Salt, L.; Chope, C.; Penson, S.; Skeggs, P.; Tosi, P.;

21 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

484

lipids in the grain of wheat (cv. Hereward) determined by lipidomic analysis of milling

485

and pearling fractions. J. Agric. Food Chem. 2015, 63, 10705-10716.

486

16)

487

Shepherd, C. E.; Hawkesford, M. J. Nitrogen efficiency of wheat: genotypic and

488

environmental variation and prospects for improvement. Eur. J. Agron. 2010, 33, 1-

489

11.

490

17)

491

profiling and pattern analysis of triacylglycerol species in Arabidopsis seeds by

492

electrospray ionization mass spectrometry. Plant J. 2014, 77,160-172.

493

18)

494

Stat. Soc. Ser. B. 1975, 37,129-145.

495

19)

496

applications. In Proceedings of the 6th Berkeley Symposium on Mathematical

497

Statistics and Probability Vol. 1, University of California Press, Berkeley, CA, 1972;

498

pp. 697-715.

499

20)

500

using a 3-step solvent extraction and acetic acid fractionation techniques. Cereals

501

2005 2005, 1, 273-276.

502

21)

503

In Wheat: Chemistry and Technology, 5th ed.; Khan, K., Shewry, P.R., Eds. AACC

504

International: St. Paul, MN, USA, 2009, pp. 363-399.

505

22)

506

analysis of lipids in cereal grains and similar tissues. J. Sci. Food Agric. 1980, 31,

507

329-340.

Barraclough, P. B.; Howarth, J. R.; Jones, J.; Lopez-Bellido, R.; Parmar, S.;

Li, M.; Baughman, E.; Roth, M. R.; Han, X.; Welti, R.; Wang, J. Quantitative

Beale, E. M. L.; Little, R. J. A. Missing values in multivariate analysis. J. R. l.

Orchard, T.; Woodbury, M. A. A missing information principle: theory and

McCann, T., Small, D. M.; Day, L. Study of lipid-protein interactions in gluten

Chung, O. K.; Ohm, J. B.; Ram, M. S.; Park, S. O.; Howitt, C. A. Wheat lipids.

Morrison, W. R.; Tan, S. L.; Hargin, K. D. Methods for the quantitative

22 ACS Paragon Plus Environment

Page 22 of 31

Page 23 of 31

Journal of Agricultural and Food Chemistry

508

23)

De la Roche, I. A.; Andrews, C. J.; Pomeroy, M. K.; Weinberger, P.; Kates, M.

509

Lipid changes in winter wheat seedlings (Triticum aestivum) at temperatures

510

inducing cold hardiness. Can. J. Bot. 1972, 50, 2401-2409.

511

24)

512

Barnes, P.J.; Ed.; London, Academic Press: London, UK, 1983, pp.11-32.

513

25)

514

effects on bread quality. Food Technol. 1968, 22 ,1897.

515

26)

516

Wheat germ in breadmaking. II. Improving breadmaking properties by physical and

517

chemical methods. Cereal Chem. 1970, 47, 429-437.

518

27)

519

interactions in gluten elucidated using acetic acid fractionation. Food Chem. 2009,

520

115, 105-112.

521

28)

522

wheat-flour nonpolar lipids. Cereal Chem. 1976, 53, 636-642.

523

29)

524

breadmaking. II. The secondary liquid lamellae. J. Cereal Sci. 2009, 49, 41-46.

525

30)

526

and degradation. Cell Mol Life Sci CMLS. 2006, 1344-1369.

527

31)

528

P.; Skeggs, P. K.; Shewry, P. R.; Wilde, P. J. Intrinsic wheat lipid composition effects the

529

interfacial and foaming properties of dough liquor. Food Hydrocolloids. 2017, in press

MacRitchie, F. Role of flour lipids in baking. In Lipids in Cereal Technology,

Pomeranz, Y.; Shogren, M.; Finney, K. F. Natural and modified phospholipids

Pomeranz, Y.; Carvajal, M. J.; Shogren, M. D.; Hoseney, R. C.; Finney. K. F.

McCann, T. H.; Small, D. M.; Batey, I. L.; Wrigley, C. L.; Day, L. Protein–lipid

De Stefanis, V. A.; Ponte, J. G. Studies on the breadmaking properties of

Sroan, B. S.; MacRitchie, F. Mechanism of gas cell stabilization in

Athenstaedt, K.; Daum, G. The life cycle of neutral lipids: synthesis, storage

Salt, L.J.; González-Thuillier, I.; Chope, G.; Penson, S.; Tosi, P.; Haslam, R.

530

23 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

531

32)

Panozzo, J. F.; O'Brien, L.; MacRitchie, F.; Bekes., F. Baking quality of

532

Australian wheat cultivars varying in their free lipid composition. J. Cereal Sci. 1990,

533

11, 51-57.

534

33)

535

effect of group 5 chromosomes on the free polar lipids and breadmaking quality of

536

wheat. J. Cereal Sci. 1989, 9, 41-51.

537

34)

538

aleurone, starch and non‐starch endosperm of four wheat varieties. J. Sci. Food

539

Agric. 1980, 31, 877-888.

540

35)

541

genotype, environment and genotype-by-environment interaction on metabolite

542

profiling in durum wheat (Triticum durum Desf.) grain. J. Cereal Sci. 2013, 57,183-

543

192.

Morrison, W. R.; Law, C. N.; Wylie, L. J.; Coventry, A. M.; Seekings, J. The

Hargin, K. D.; Morrison, W. R. The distribution of acyl lipids in the germ,

Beleggia, R.; Platani, C.; Nigro, F.; De Vita, P.; Cattivelli, L.; Papa, R. Effect of

544

545

24 ACS Paragon Plus Environment

Page 24 of 31

Page 25 of 31

Journal of Agricultural and Food Chemistry

Table 1 The p-values for F-tests of cultivar, N and cultivar by N interactions. Lipids

Cultivar

N

Cultivar x Nitrogen (GxE)

LPC Total

0.651

0.950

0.997

PC Total