Research Article HRMAS NMR and MRI ... - ACS Publications

were processed by using both Bruker Topspin Software (v 2.1, BrukerBiospin, Rheinstetten,. 161. Germany) and MNOVA .... Concomitantly, the closeness o...
1 downloads 12 Views 1MB Size
Subscriber access provided by READING UNIV

Article

HRMAS NMR and MRI spectroscopies distinguish metabolome and structural properties of maize seeds from plants treated with different fertilizers and arbuscular mycorrhizal fungi Pierluigi Mazzei, Vincenza Cozzolino, and Alessandro Piccolo J. Agric. Food Chem., Just Accepted Manuscript • Publication Date (Web): 11 Jan 2018 Downloaded from http://pubs.acs.org on January 11, 2018

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 32

Journal of Agricultural and Food Chemistry

Research Article

1

HRMAS NMR and MRI spectroscopies distinguish metabolome and structural properties of maize seeds from plants treated with different fertilizers and arbuscular mycorrhizal fungi

2 3 4 5 6

Pierluigi Mazzei†,*, Vincenza Cozzolino‡, Alessandro Piccolo†, ‡, *

7 8 9 10 11 12 13 14 15 16 17 18

†Centro Interdipartimentale per la Risonanza Magnetica Nucleare per l’Ambiente, l’Agro-

Alimentare ed i Nuovi Materiali (CERMANU), Università di Napoli Federico II, Via Università 100, 80055 Portici, Italy ‡

Dipartimento di Agraria, Università di Napoli Federico II, Via Università 100, 80055 Portici, Italy

19 20 21 22 23 24 25 26

Keywords: HRMAS, MRI, Metabolomics, Maize Caryopses, Fertilization, Compost, Arbuscular Mycorrhizal Fungi.

27

1

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

28

ABSTRACT: Both High-Resolution Magic-Angle-Spinning (HRMAS) and Magnetic Resonance

29

Imaging (MRI) NMR spectroscopies were applied here to identify the changes of metabolome,

30

morphology and structural properties induced in seeds (caryopses) of maize plants grown at field

31

level under either mineral or compost fertilization in combination with the inoculation by Arbuscular

32

Mycorrhizal Fungi (AMF). The metabolome of intact caryopses was examined by HRMAS-NMR,

33

while the morphological aspects, endosperm properties and seed water distribution were investigated

34

by MRI. Principal Component Analysis was applied to evaluate 1H CPMG (Carr-Purcel-Meiboom-

35

Gill) HRMAS spectra as well as several MRI-derived parameters (T1, T2 and self-diffusion

36

coefficients) of intact maize caryopses. PCA score-plots from spectral results indicated that both

37

seeds metabolome and structural properties depended on the specific field treatment undergone by

38

maize plants. Our findings show that a combination of multivariate statistical analyses with advanced

39

and non-destructive NMR techniques, such as HRMAS and MRI, enable the evaluation of the effects

40

induced on maize caryopses by different fertilization and management practices at field level. The

41

spectroscopic approach adopted here may become useful for the objective appraisal of the quality of

42

seeds produced under a sustainable agriculture.

43

2

ACS Paragon Plus Environment

Page 2 of 32

Page 3 of 32

44

Journal of Agricultural and Food Chemistry

■ INTRODUCTION

45 46

The identification of efficient and sustainable strategies to promote plants growth is a major

47

challenge in the agrofood sector, especially for the most important food crops such as maize, wheat

48

and rice.1 A traditional practice is the amendment to soils of mineral fertilizers containing the

49

essential macronutrients, whose scarcity may otherwise represent a growth limiting-factor for

50

plants.2 Large if not excessive amount of mineral fertilizers are commonly applied in order to

51

maintain high crop yields, and compensate both the progressive P immobilization3 in soils and the

52

consumption due to microbial activity,4 but leading to costly and environmentally negative

53

conditions.5 An alternative to traditional fertilizers is the amendment of nutrients in their organic

54

form,6 such as in compost that is a biochemically stable organic material rich in organic P and N

55

forms.7-10

56

An increase in plant nutrient uptake is promoted by soil inoculation with beneficial

57

microorganisms, such as Arbuscular Mycorrhizal Fungi (AMF), which enhance the degree of

58

sustainability in industrial agriculture.11 AMF are widespread in most terrestrial ecosystems where

59

they form symbiotic associations with the majority of plants. AMF protect plants against pathogens

60

and facilitate nutrients uptake through several mechanisms, including the formation of an extensive

61

extraradical mycelium, the stimulation of roots elongation, and the induced solubilization of poorly

62

bioavailable nutrients.12 In fact, AMF are regarded as a means to reduce the excessive supply of

63

fertilizers, thus enhancing, directly and indirectly, the natural soil fertility in agricultural

64

ecosystems.13

65

Despite the attempts to introduce sustainable practices in agriculture, an exhaustive knowledge is

66

still lacking on the molecular mechanisms that make the alternative practices effective and how they

67

influence plants metabolic responses and their products. We believe that the use of advanced

68

molecular analytical techniques such as Nuclear Magnetic Resonance (NMR) spectroscopy may 3

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

69

contribute to identify and understand the changes occurring in plants due to field treatments with

70

innovative sustainable practices. In particular, HRMAS (High Resolution Magic Angle Spinning)

71

and MRI (Magnetic Resonance Imaging) represent non-destructive NMR techniques that may

72

provide direct molecular and structural information on plant fresh samples.14-16

73

HRMAS retains the advantages of both traditional solid- and liquid-state NMR techniques

74

and allows the direct examination of whole fresh samples without the need of preliminary

75

extractions.17-18 Moreover, HRMAS 1D spectra were proved useful for the Metabolomics analysis of

76

biological samples and the dynamic evaluation of treatment-dependent metabolic processes

77

occurring in plants.16, 19-22 With regard to maize plants, NMR-based metabolomics was employed to

78

study genetically modified samples,23-25 and to monitor the influence exerted by stress conditions,

79

such as excess salts26 or nitrogen deficiency.27

80

As a complementary NMR technique, MRI enables acquisition of morphological and

81

anatomical images of intact samples at a µm level and permits to extrapolate structural information

82

based on mobility and compartmentalization of constitutive water.28 Currently, MRI is mainly used

83

on mammalians within medical, diagnostic and pharmaceutical applications. In fact, MRI is a

84

consolidated analytical method for the rapid non-invasive evaluation of changes in tissues

85

morphology, and occurrence of specific diseases, such as cancers development.28-29

86

Despite its specific large potential, MRI is still underutilized in the research on plants and

87

food.30-31 Only few works have so far reported on MRI applications for maize plants ad seeds. For

88

example, the effects of steeping processes on maize caryopses were followed by 2D images and MRI

89

moisture profile,32 while 1H MRI was applied to investigate the cell water balance of maize plants

90

during osmotic stress.33 MRI images and relaxation times were used to examine changes in the

91

physical state of water of maize seeds when dried at sub-zero temperatures.34 MRI was applied to

92

design a physical model of heat and mass transfer in maize caryopses in order to assess the

93

dehydration degree and quality degradation of seeds during fluidized-bed drying.35 4

ACS Paragon Plus Environment

Page 4 of 32

Page 5 of 32

Journal of Agricultural and Food Chemistry

94

The aim of this work was thus to apply HRMAS and MRI spectroscopies to evaluate the effects

95

on maize caryopses of field amendments with mineral or compost fertilization and of inoculation

96

with a selected AMF consortium.

97

98

■ MATERIALS AND METHODS

99

Field experiment. Caryopses were collected from maize plants (Zea mays L. cv 3321

100

Limagrain) of field experiments conducted at the Experimental Station of University of Naples

101

Federico II, near Castel Volturno (CE), Campania, Italy. The site is governed by a Mediterranean

102

climate, with moist cool winters and warm dry summers. The soil was classified as Typic Vertic

103

Xerofluvent, with a clayey texture (22% sand, 28% silt, and 50% clay). The field experiment

104

consisted in 5 x 8 meter plots, which received the following treatments: (i) control (M-), (ii) 15 t ha-1

105

of compost made with organic urban wastes (COM), (iii) a mineral fertilizer (MIN) (N250P150K160)

106

composed by 250 kg ha-1 of N (as urea), 150 kg ha-1 of P2O5 (as Ca(H2PO4)2) and 160 kg ha-1 of K

107

(as K2SO4). Moreover, soils were also treated with AMF inoculation (+) (Arbuscular

108

Mycorrhizal Fungi, as a commercial fungal inoculum), either (iv) without (M+) or (v) with mineral

109

fertilizer (MIN+) or (vi) compost (COM+). The commercial inoculum contained 50 AMF infective

110

propagules g-1 of product (information provided by the supplier) and consisted of a mixture of 50%

111

Rhizophagus irregularis and 50% Funneliformis mosseae (previously named Glomus intraradices

112

and Glomus mosseae, respectively) (AEGIS®, Italpollina, Rivoli Veronese, Italy). This microbial

113

mix was defined on the basis of preliminary field trial-tests (data not shown) which selected the

114

combination with the greatest adaptability and performance. Each treatment was replicated three

115

times, resulting in 18 total experimental plots, which were randomly distributed (complete

116

randomized design). Mineral fertilizer and compost were manually broadcasted and incorporated in 5

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

117

soil. The compost was broadcasted three weeks before sowing, while the inocula (20 kg ha−1) were

118

applied during sowing by a seed-drill microgranulator. Maize seeds were linearly sowed at 10 cm

119

distance, while the distance among rows was 75 cm. The plant density corresponded to 7 plants per

120

m2.

121

122

Caryopses collection and sample preparation. A representative number of cobs per plot

123

were harvested during the milky ripeness stage of maize seeds (80 days after sowing) by sampling

124

only plants located in the centre of each plot. A sufficient number (>100) of caryopses were gently

125

removed from cobs of each treatment, and stored at -80°C until NMR analysis. Maize caryopses for

126

MRI analysis were directly loaded into a 10mm MRI tube after a 30 min of defrosting time. Maize

127

samples for HRMAS analysis were prepared by cutting longitudinally the pericarp and aleurone

128

outer layers of each caryopsis and collecting about 15 mg of fresh material with a spongy

129

consistency. In particular, the latter was prevalently composed by the endosperm that typically

130

represents the 69-83% of maize caryopsis dry weight.36 Each HRMAS analytical replicate

131

corresponded to material isolated from a single seed. The material was then packed into a 4 mm

132

zirconia HRMAS rotor, fitted with a perforated Teflon insert, soaked with approximately 15 µL of

133

D2O solution (99.8% D2O/H2O, Armar Chemicals) and sealed with a Kel-F cap (Rototech-Spintech

134

GmbH). A spin rate of 5000 Hz±1 was adopted to rotate samples. A number of 4 replicates for each

135

field parcel were analysed, thus resulting in 12 total replicates for each plant treatment.

136

137

HRMAS NMR Experiments. All HRMAS NMR experiments were conducted at 298±1K on

138

a 400 MHz Bruker Avance magnet (Bruker Biospin, Rheinstetten, Germany), equipped with a

139

Bruker 1H-13C HRMAS probe working at a 1H frequency of 400.13 MHz. A Carr-Purcel-Meiboom6

ACS Paragon Plus Environment

Page 6 of 32

Page 7 of 32

Journal of Agricultural and Food Chemistry

140

Gill (CPMG) NMR pulse sequence was used to acquire 1H spectra of maize caryopses. This

141

sequence was preferred to a conventional 1H acquisition since it allowed application of a T2-filter to

142

enable a selective suppression of overlapping signals due to compounds with short spin-spin (T2)

143

relaxation time. The experiments were acquired by setting 2 seconds of recycle delay, a 90° pulse

144

length ranging within 6.5 and 8.6 µs, 16384 points, a spectral width of 16 ppm (6410.3 Hz) and 256

145

scans. In particular, the CPMG pulse sequence, that is based on a spin-echo method, was performed

146

with a total spin-spin relaxation delay (2nτ) of 320 ms, including single optimal echo times (τ) of 2

147

ms. The signal of residual water was suppressed by an on-resonance pre-saturation during thermal

148

equilibrium delay. Structural identification of maize seeds metabolites was assessed by 2D NMR

149

experiments, consisting in homo-nuclear 1H–1H COSY (Correlation SpectroscopY), TOCSY (Total

150

COrrelationSpectroscopY) and J-RES as well as hetero-nuclear 1H–13C HSQC (Hetero-nuclear

151

Single-Quantum Correlation) and HMBC (Hetero-nuclear Multiple Bond Correlation). All 2D

152

experiments had a spectral width of 16 (6410.3 Hz) and 300 (30186.8 Hz) ppm for 1H and 13C nuclei,

153

respectively, and a time domain of 2048 points (F2) and 256 experiments (F1). In particular, all

154

homo-nuclear 2D spectra consisted in 16 dummy scans and 64 total transients and, in case of

155

TOCSY experiment, a mixing time of 80 ms and a trim pulse length of 2500 ms. HSQC and HMBC

156

hetero-nuclear experiments were acquired with 16 dummy scans, 80 total transients and 0.5 µs of

157

trim pulse length. The optimization of these heteronuclear experiments was conducted by taking into

158

account 145 and 6.5 Hz as optimal 1H –13C short and long range J-couplings, respectively. A phase

159

and baseline correction was applied to all mono- and bi-dimensional spectra. 1H 1D spectra were

160

Fourier Transformed without any apodization and by applying a two-fold zero-filling. NMR data

161

were processed by using both Bruker Topspin Software (v 2.1, BrukerBiospin, Rheinstetten,

162

Germany) and MNOVA Software (v.9.0, Mestrelab Research, Santiago de Compostela, Spain).

7

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

163

Magnetic Resonance Imaging (MRI) Experiments. All MRI experiments were performed at

164

298±1K on a 300 MHz Bruker Avance wide-bore magnet (BrukerBiospin, Rheinstetten, Germany)

165

and equipped with a 10 mm µ-imaging MICRO 5 probe working at a 1H frequency of 300.13 MHz.

166

T2 experiments were performed by applying an echo‐train MSME (Multi‐Slice‐Multi‐Echo) pulse

167

sequence, setting 3 s of recycle delay, 10 scans and acquiring 16 experiments with an increasing

168

number of spin-echoes (total spin-echo time ranging within 10.57 and 169.12 ms). Spin density

169

measurements were performed by considering the images obtained at the shortest echo-delay (10.57

170

ms) in MSME experiments, whereas T2-weighted images were achieved by selecting a total echo-

171

delay of 42.3 ms. T1 measurements were acquired by RARE (Rapid Acquisition with Relaxation

172

Enhancement) experiments which consisted of 10 scans, a rare-factor of 2 and increasing repetition

173

times ranging within 0.188 and 7.5 s in 8 experiments. Diffusion measurements were performed by

174

using a Pulse Field Gradient (PFG) sequence integrated with stimulated echoes. This experiment

175

was set up with 3 s of recycle delay, 5 scans, δ and ∆ diffusion delays of 4 and 13 ms, respectively,

176

and applying gradients at increasing strengths (0, 4.5, 112.0, 362.9, 757.3, 1295.0 and 2168.8 s mm-

177

2

178

gradient pulses in longitudinal (z) or transverse (y) direction, it was assumed an isotropic apparent

179

diffusion in seeds. Consequently, the gradients for PFG experiments have been applied only in the z

180

direction in all cases. Diffusion-weighted images were acquired by selecting the gradient strength of

181

362.9 s mm-2. In all cases the images consisted of a 256x256 matrix produced by acquiring 6

182

interlaced coronal (and sagittal, in some cases) slices with a width of 1 mm, a field of view of

183

1.33x1.0 mm, and 1.07 mm of slice-interdistance.

) in 7 experiments. Since no significant differences in diffusivity were observed by applying

184

3D experiments (rare 3D pulse sequence, 128x128x128 matrix, field of view 13x13x13 mm,

185

rare factor of 8, 24 scans, 3 s of recycle delay and 3.075 ms of echo time) were performed only for

186

several representative samples as well as relaxation and diffusion experiments were repeated by 8

ACS Paragon Plus Environment

Page 8 of 32

Page 9 of 32

Journal of Agricultural and Food Chemistry

187

using also axial slices. ParaVision 5.1 Bruker software was used to process MRI data, attain all of the

188

images, calculate T1 and T2 relaxation times as well as the apparent diffusion coefficients.

189

190

Multivariate data analyses. The spectral region ranging from 0.1 to 9.22 ppm in 1H CPMG

191

HRMAS NMR spectra was equally divided into 228 segments (each single bucket width

192

corresponded to 0.04 ppm). All these segments were integrated, except the region of the baseline

193

distortion resulting from water suppression (4.9-4.66 ppm). The spectral integration produced a data

194

matrix composed by 228 variables x 72 observations (12 replicates for 6 treatments). Data were

195

normalized, by dividing each segment area by the sum of all signals areas, and Pareto-scaled prior to

196

perform Principal Component Analysis (PCA).37-38

197

Also MRI parametric data (apparent diffusion coefficients, T1 and T2 relaxation times) were

198

calculated and evaluated by PCA. In detail, they were selected one large (5 cm2) and two small (1

199

cm2) sized Regions of Interest (ROIs) for both vitreous and floury endosperm. For each MRI

200

parameter and endospermic compartment, it was calculated the average value among the three above

201

mentioned ROIs. The resulting averaged values were then included in the data matrix for MRI-based

202

PCAs which was composed by 6 variables x 36 observations (6 replicates per treatment).

203

In all cases, the variables mostly involved in the differentiation among treatments were

204

identified on the basis of PCA loading-plots, while their significance (Tukey Test, confidence level

205

>95%) was assessed separately by ANOVA (ANalysis Of VAriance) test. Statistical elaborations

206

were carried out by the XLStat software (v. 2012, Addinsoft, Paris, France).

207

208

■ RESULTS AND DISCUSSION

9

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

209

HRMAS. A 1H HRMAS NMR spectrum of a representative maize caryopsis (M-) (Figure 1)

210

was obtained by a CPMG NMR pulse sequence to selectively suppress broad signals, which are

211

typically due to relatively large lipid and protein molecules, such as zein, albumin and glutelin.36

212

Elimination of broad signals improved the proton spectral resolution, thus resulting into easier

213

compound assignment and reproducible quantitative elaboration. Identification of principal signals

214

(Figure 1) was conducted on the basis of previous literature assignments23,26,39 and supported by the

215

homo- and heteronuclear 2D NMR spectra of this study (data not shown). Primary metabolites, such

216

as carbohydrates, amino acids, alcohols and several organic acids, were identified in spectra of fresh

217

caryopses (Figure 1). The most intense proton signals were ascribed to carbohydrates and consisted

218

in overlapped multiplets ranging from 5.5 to 3.5 ppm. This was expected since the analyzed seed part

219

was mostly represented by the endosperm that is prevalently composed by carbohydrates and

220

starch.40 The multiplets detected in the spectral region between 2.35 and 2.01 ppm were assigned to

221

proline\glutamate, as superimposition of proton signals belonging to both glutamate (βCH2 and

222

γCH2) and proline (β’CH2 and βCH) amino acids.23

223

While no specific molecular markers were identified to enable a distinction among

224

treatments, the quantitative evaluation of the relative concentration of several metabolites was found

225

to vary according with plant treatments, thus indicating an induced modulation of the maize

226

metabolism. This became evident by processing the 1H-HRMAS NMR spectra by Principal

227

Component Analysis (PCA). PCA is an unsupervised pattern-recognition technique to evaluate the

228

intrinsic variation within different classes of samples. In fact, this multivariate elaboration calculates

229

linear combinations of a starting set of variables on the basis of their maximum variance and reduces

230

the dimensionality of the original data matrix, while retaining the maximum amount of variability, as

231

well as most of the original information contained in the data set.37-38 Moreover, PCA offers the

232

practical advantage of exploring, in a single output, referred to as score-plot, the molecular response

233

induced in numerous samples on the basis of a large number of variables.38,41 10

ACS Paragon Plus Environment

Page 10 of 32

Page 11 of 32

Journal of Agricultural and Food Chemistry

234

The PCA facilitated the direct comparison among the HRMAS spectral data set of this study,

235

which consisted in a relatively dense matrix composed by 228 variables (buckets) and 72

236

observations (12 replicates x 6 treatments). Three PCA score-plots related to samples deriving from

237

MIN (Figure 2A), COM (Figure 2B) and M+ (Figure 2C) field treatments, indicated significant

238

differences (Tukey Test in ANOVA, confidence level > 95%) as compared to M- samples. Each

239

score-plot shows the direction assumed by the most relevant loading vectors (variables) which are

240

responsible for the differentiation among treatments (Figures S1A, B and C of the Supporting

241

Information, respectively). The almost neat separation among samples of different treatments

242

suggests that field treatments systematically influenced the metabolome of maize caryopses.

243

Concomitantly, the closeness of replicates deriving from the same treatment indicates a good

244

reproducibility of HRMAS technique (Figure 2). Remarkably, the statistical reliability of our results

245

is supported by the fact that the 12 replicates for each treatment were collected from three different

246

subplots (4 field-replicates for each specific parcel).

247

The score-plot of M- versus MIN samples (Figure 2A) results from the combination of PC1

248

and PC3, accounting for 50.6% of the total explained variance. According to the direction of the

249

most relevant loading-vectors (Figure S1A of the Supporting Information), MIN samples were

250

characterized by a larger amount of alanine, proline\glutamate and tryptophan amino acids than for

251

M- (Figure 2A). This abundance of amino acids may be related to a more efficient nitrogen plant

252

uptake, resulting from a larger availability of readily soluble N in soils amended with the mineral

253

fertilizer.42

254

The relationship between M- and COM samples was shown by the PCA score-plot that was

255

obtained by combining PC3 with PC11 (Figure 2B). Despite the relatively lower explained variance,

256

this score-plot showed the differences occurring between M- and COM samples, and revealed that

257

COM samples contained a larger amount of isoleucine, but a lower concentration of alanine,

258

proline\glutamate, leucine and citric acid than for M- (Figure S1B of the Supporting Information). 11

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

259

This result may be again related to the immediate availability of nitrogen42 that, in the case of

260

compost, should be first mineralized from its organic forms, thus leading to a nitrogen plant uptake

261

capacity even lower than for M- (Figure 2B).43

262

The PCA score-plot of Figure 2C highlighted that AMF inoculation significantly influenced

263

the metabolome of maize caryopses. In fact, M+ samples exhibited a larger concentration of sucrose

264

and glucose, accompanied by a lower amount of citric acid, succinate, alanine and proline\glutamate

265

(Figure S1C of the Supporting Information). Assuming that also the plant roots of M- treatment were

266

likely associated with autochthonous mycorrhizal fungi, the addition of a commercial inoculum to

267

M+ plants must have produced a more pronounced variation in the caryopsis metabolome. In

268

particular, the larger abundance of glucose and sucrose observed in M+ suggests a more efficient

269

production and translocation of carbohydrates to seeds. In fact, mycorrizal association elicits the

270

expansion of surface, volume, and roots number, with a consequent enhancement of plant uptake of

271

water and nutrients,

272

relatively lower levels of glucose and sucrose in seeds from AMF inoculated plants, which were

273

accompanied by a larger amount of metabolites related to aerobic glycolysis, such as succinate and

274

citrate. In addition, the larger abundance of several amino acids in M- caryopses also suggests an

275

increased biosynthesis of glycolytic enzymes associated with an increase in respiration and due to a

276

larger availability of nitrogen in M- treatment.42-44

5, 12

and improvement of plant productivity. This explanation agrees with the

277

The effects of the synergy of compost with AFM inoculation (COM+) are worth noting due

278

to its environmental relevance. We report in Figure 3 three PCA score-plots, that differentiate the

279

metabolome of the COM+ treatment from those of COM (Figure 3A), M- (Figure 3B) and M+

280

(Figure 3C). Each score-plot shows the direction assumed by the most relevant loading vectors

281

(variables) which are responsible for the differentiation among treatments (Figures S2A, B and C of

282

the Supporting Information, respectively).The PCA score-plot in Figure 3A was built on the

283

combination of PC3 and PC4 (16.9% of total explained variance) and revealed that the addition of 12

ACS Paragon Plus Environment

Page 12 of 32

Page 13 of 32

Journal of Agricultural and Food Chemistry

284

compost with or without AMF inoculation brought to a different metabolomic response in maize

285

caryopses. In fact, except for the two COM samples associated to negative values of PC3 and PC4

286

(placed in the 3rd quadrant), the differentiation occurred prevalently along the PC3 and was due to a

287

lower amount of glucose in COM, accompanied by a larger amount of phenylalanine (Figure S2A of

288

the Supporting Information). The abundance of glucose in seeds may be related to a more efficient

289

photosynthetic activity in plants under the COM+ treatment. This agrees with previous works that

290

showed that compost may support and promote AMF action, that, in turn, favors plant

291

photosynthesis.45-46

292

The comparison between COM+ and M- treatments is shown in the PCA score-plot of Figure

293

3B. In this case, only the abundance of citric acid (associated to positive PC3) enabled to

294

significantly differentiate the metabolome of COM+ caryopses from that of control (Figure S2B of

295

the Supporting Information). The excess of citric acid in M- samples, as already shown to be

296

responsible for the differentiation between M- and COM (Figure 2B), may be accounted to an

297

enhanced glycolytic process. Figure 3C reports the PCA score-plot of caryopses collected from

298

maize plants that were subjected to AMF inoculation either with (COM+) or without (M+) compost

299

fertilization. The combination of PC2 vs PC4 (20.8% of total explained variance) distinguished the

300

two treatments along the direction intersecting the 2nd and 4th quadrant of the score-plot. In

301

particular, COM+ caryopses were differentiated due to a larger content of citric acid, proline\

302

glutamate, alanine and leucine, and a smaller concentration of sucrose, than for M+ (Figure S2C of

303

the Supporting Information). Contrary to the observation that sucrose and glucose were more

304

abundant in M+ than in M- samples (Figure 2C), no significant difference in glucose content was

305

found between M+ and COM+ caryopses (Figure 3C), although sucrose was very abundant in M+

306

seeds. On the other hand, the excess of several amino acids in COM+ samples (Figure 3C) can be

307

explained by a larger nitrogen availability due to compost amendment, that served not only as N

308

source, but also promoted a beneficial AMF action. 13

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 14 of 32

309

A PCA was also performed on HRMAS spectral data of seeds from maize plants which had

310

been amended with the mineral fertilizer either with (MIN+) or without (MIN) AMF inoculation

311

(Figure S3 of the Supporting Information). In this case, neither the score-plot resulting from the

312

combination of PC1 and PC2 (Figure S3 of the Supporting Information, 49% total variance

313

explained), nor other possible score-plots exploiting the remaining PCs (data not shown), led to a

314

neat differentiation between MIN and MIN+ caryopses. This suggests that the effect exerted on seeds

315

metabolome by the mineral fertilizer was predominant, thus limiting significantly the AMF

316

inoculation.

317

318

Magnetic Resonance Imaging. MRI is a non-invasive NMR spectroscopy technique capable

319

to provide morphological and anatomical details of intact maize caryopses, and indirectly reveal

320

possible effects of field treatments on plant and seed metabolism.

321

images of a intact M- maize caryopsis acquired by applying the MSME pulse sequence at a minimal

322

spin-echo time (10.57 ms, no T2-weighted). The sagittal slice positioned in the center of the seed

323

allowed a neat differentiation among the main regions composing the internal structure of the

324

caryopses, such as the protective pericarp, the embryo region (including the embryonic radicle, the

325

plumule and scutellum) and the endosperm (Figure 4A). The endosperm surrounded by the aleurone

326

layer covers most of the seed volume and is divided into two compartments, referred to as vitreous

327

endosperm (EV) and floury endosperm (EF), which have different composition, density and

328

consistency. Typically, EV is rich in lipids and proteins (including zein) whereas starch, mainly

329

composed by amylose and amylopectine, is predominant in EF. 48 However, the properties and extent

330

of separation of these two endosperm compartments depend strictly on the specific maize cultivar

331

and its quality.49 For example, in the case of the cultivar used here, EV and EF appeared neatly

332

separated from each other (Figure 4A), since EV is concentrated in the seed’s lower part near the 14

ACS Paragon Plus Environment

47

Figure 4 reports three MRI

Page 15 of 32

Journal of Agricultural and Food Chemistry

333

embryo region, whereas EF occupies most of upper seed volume and appears structurally more

334

heterogeneous.49 Two axial slices acquired at different positions of the maize seed (as indicated in

335

Figure 4A) show the core of EF (Figure 4B) and the homogeneous region of EV (Figure 4C).

336

Moreover, the internal structure of maize caryopses are reported in two video clips which display the

337

3D projection of seed (video clip 1 of the Supporting Information) and the bottom-up sequence of

338

axial slices (video clip 2 of the Supporting Information).

339

A coronal and central slice of a control caryopses is shown in Figure 5A, from which it can

340

be excluded the superimposition of the embryo region, thereby making both EF and EV regions more

341

clearly visible. Interestingly, the intrinsic properties of this maize caryopses enabled the systematic

342

isolation of EF and EV regions by acquiring T2-weighted (total spin-echo delay of 42.3 ms) (Figure

343

5B) and diffusion-weighted (gradient strength of 362.95 s mm-2) (Figure 5C) images, respectively. In

344

fact, EF and EV showed a neat difference in both T2 relaxation times and diffusivity, thus suggesting

345

that the detection of specific anatomical characteristics for each endosperm region may lead to a

346

more precise spin-density evaluation. In fact, we compared both conventional (no weighted) and T2-

347

weighted images in order to identify possible peculiar properties induced in maize seeds by the

348

applied field treatments. Although several differences in morphology were observed in several maize

349

seeds, we estimated that the variations could not be related to treatments, but rather to the natural

350

variability in maize seeds.

351

Nevertheless, several other NMR parameters, such as spin-lattice (T1), and spin-spin (T2)

352

relaxation times, and the apparent diffusion coefficients (Diff), were shown to reveal the spatial

353

distribution of the physical status of water in different seed compartments, thus providing detailed

354

and quantifiable information on seeds molecular arrangement.30,33,50 While relatively long T1 and T2

355

relaxation times are attributed to an intrinsic large molecular and structural flexibility in tissues, short

356

relaxation times are due to smaller tumbling rates of water molecules interacting with either

357

macromolecules or solid surfaces. On the other hand, Diff coefficients estimate the apparent 15

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

358

translational self-diffusion (conventionally referred to as ‘diffusion’) of water molecules into tissue

359

compartments, thus indirectly indicating the extent of free diffusion within the tissue matrix.

360

Consequently, Diff coefficients may be decreased if water molecules experience barriers to diffusion,

361

such as rigid tissue structures, interact with cell membranes, or diffuse within small-sized cells.

362

Alternatively, Diff values may be increased if water molecules move more freely into a less compact

363

endosperm (longer free diffusion distance). On these bases, differences in relaxation times and

364

diffusion coefficients measured by MRI are likely to reflect different composition, structure and

365

homogeneity of endosperm samples. In turn, these properties are a function of multiple treatment-

366

dependent factors, such as water/nutrients availability during seed development, extent of seed

367

maturation, compounds partition in the endosperm, distribution of starch and protein reserves for

368

germination.

369

The measurements of relaxation times and diffusion coefficients were conducted in 1 and 5

370

cm2 sized regions of interest (ROIs) selected in endosperms of maize caryopses (Figure S4 of the

371

Supporting Information). In all cases, three ROIs were used for both vitreous (V1, V2 and V10) and

372

floury (F1, F2 and F10) endosperms. The inclusion of one large ROI (V10 and F10, respectively) for

373

each endosperm compartment had the rationale to minimize possible site-specific artifacts resulting

374

from local non-homogeneities in smaller ROIs.

375

The PCA score-plot built up on MRI spectral measurements for the caryopses from the M-,

376

M+, MIN and COM treatments (Figure 6A) was obtained by combining the 1st and the 3rd principal

377

components (49.3 and 14.8 % of explained variance, respectively). A neat differentiation among M-,

378

COM, and M+ samples was observed from positive to negative values along the PC1, that

379

corresponded to the progressive decrease of T1-EF, T2-EF and Diff-EF values (Figure S5A of the

380

Supporting Information). While the caryopses from the MIN treatment could not be differentiated

381

from those of M- along the PC1, they resulted sufficiently separated along the PC3 due to larger and

382

smaller Diff-Ev and T1-EV values than for the M- treatment. The fact that the floury endosperm of M+ 16

ACS Paragon Plus Environment

Page 16 of 32

Page 17 of 32

Journal of Agricultural and Food Chemistry

383

samples exhibited the shortest T1, T2 and smallest Diff values suggests a slower mobility of water

384

molecules in EF, that may be ascribed to a more rigid starchy structure than in caryopses of other

385

treatments. Conversely, although the MRI spectral measurements for the floury endosperm in MIN

386

samples were similar to those for the M- treatment, the vitreous endosperm of maize seeds from the

387

treatment with mineral fertilization provided larger Diff-EV and smaller T1-EV values than for the M-

388

caryopses. This finding related to the MIN treatment may be then explained with a larger molecular

389

rigidity of EV (short T1) compartment and a larger presence of water molecules (high Diff) than for

390

the M- treatment.

391

The PCA score-plot related to maize seeds from compost-treated plants is reported in Figure

392

6B, as based on PC1 and PC2 (39.7 and 22.9 % of variance, respectively). This score-plot revealed a

393

neat horizontal differentiation between COM, COM+ and M- seeds along the PC1 especially due to

394

the characteristics of vitreous endosperm. Moreover, while both COM and COM+ samples provided

395

larger values of Diff-Ev and T1-EV and T2-EV than for M-, the relaxation times and diffusion

396

coefficients were even larger for the caryopses from maize plants which received the combined

397

addition of compost and AMF (COM+) (Figure S5B of the Supporting Information). Conversely, M-

398

seeds revealed the largest diffusion coefficients in the floury endosperm, in line with what is shown

399

in Figure 6A. Also a general separation between COM and COM+ along the PC2 was due to

400

differences in T1-EF values. The MRI spectral measurements thus indicate that field treatments with

401

compost had a significant effect not only on the composition (Figure 3) but also on the

402

morphological structure of caryopses, and this was magnified by the mycorrizal inoculation.

403

The PCA score-plot for the comparison of MRI data between MIN and MIN+ field treatments

404

and those of M- control is shown in Figure 6C. The difference between the two MIN treatments and

405

control was evident along PC1 (33.38% of explained variance), although no differentiation was

406

possible between MIN+ and MIN samples, as already noted from HRMAS spectral results (Figure

407

S3 of the Supporting Information). In line with a previous score-plot (Figure 6A), the difference in 17

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

408

caryopses from plants treated with mineral fertilizers from those of control was represented by both

409

larger Diff-Ev and smaller T1-EV values for MIN treatments (Figure S5C of the Supporting

410

Information). Interestingly, on the basis of HRMAS results, which differentiated MIN from M-

411

samples because of a larger amount of several free amino acids (Figure 2), such a different response

412

in MIN and MIN+ samples may be ascribed to a different proteins/amino acids ratio in vitreous

413

endosperm, that is likely dependent on the different N-availability .

414

In this work we showed that both HRMAS-NMR and MRI spectroscopies enabled the

415

differentiation, respectively, of the molecular composition and morphological structure of caryopses

416

from maize plants, which were treated with either mineral or compost fertilization, as well as with

417

the inoculation by mycorrhizal fungi.

418

1

H HRMAS spectra provided the metabolomic profile of caryopses, whose spectral features

419

enabled the elaboration of PCA and the identification of the metabolite related to the specific field

420

treatments. Concomitantly, the adopted pulse techniques of MRI spectroscopy revealed the

421

morphological and anatomical images of intact maize caryopses. Although no significant differences

422

were found in seeds morphology, the relaxation times (T1, T2) and the diffusion coefficients (Diff)

423

measured by MRI suggested that the physical status of water in the endosperm compartments of

424

seeds were systematically related to the field treatments undergone by the corresponding maize

425

plants. The spectral information on maize obtained here by both HRMAS and MRI spectroscopies

426

were particularly relevant since they relied on samples from a statistically significant experimental

427

design at field level.

428

Our results indicate that the applied NMR techniques can be used to evaluate the molecular

429

composition and structural properties of seeds from maize plants that are subjected to different field

430

management practices. This may be of importance in the increasing quest for food security that

431

requires an objective assessment of the quality of products of sustainable agricultural practices.

18

ACS Paragon Plus Environment

Page 18 of 32

Page 19 of 32

Journal of Agricultural and Food Chemistry

432

433

■ ASSOCIATED CONTENT

434

*S Supporting Information

435

The Supporting Information is available free of charge on the ACS Publications website at DOI:…

436

Supporting information includes a supplementary PCA score-plot evaluating MIN vs MIN+

437

treatments, the most significant loading vectors associated to PCA score-plots shown in Figures 2, 3

438

and 6, the representative location of MRI ROIs, two video clips showing the MRI-detected internal

439

structure of a representative maize seed.

440

441

442

443

444

■ AUTHOR INFORMATION

445

* Corresponding Authors

446

Pierluigi Mazzei: Centro Interdipartimentale per la Risonanza Magnetica Nucleare per l’Ambiente,

447

l’Agro-Alimentare ed i Nuovi Materiali (CERMANU), Università di Napoli Federico II, Via

448

Università 100, 80055 Portici, Italy.E-mail: [email protected]. Phone:39-081-2539448.

449

ORCID (ID)

450

Pierluigi Mazzei: 0000-0002-5312-4969

451

19

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

452

Alessandro Piccolo:

Centro Interdipartimentale per la Risonanza Magnetica Nucleare per

453

l’Ambiente, l’Agro-Alimentare ed i Nuovi Materiali (CERMANU), Università di Napoli Federico II,

454

Via Università 100,80055 Portici, Italy. E-mail: [email protected]. Phone:39-081-

455

2539160

456

Notes

457

The authors declare no competing financial interest.

458

459

■ ACKNOWLEDGEMENTS

460

This work was supported by a project within the 2010 “FARO” program of the University of Naples

461

Federico II.

462

20

ACS Paragon Plus Environment

Page 20 of 32

Page 21 of 32

463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511

Journal of Agricultural and Food Chemistry

■ REFERENCES (01) Khan, N.; Peiqiang, Y.; Mubarak, A.; Cone, J. W.; Hendriksa, W. H. Nutritive value of maize silage in relation to dairy cow performance and milk quality. J. Sci. Food. Agric. 2015, 95, 238– 252. (02) Messias, R. S.; Galli, V.; Silva, S. D. D. A. E.; Schirmer, M. A.; Rombladi, C. V. Micronutrient and functional compounds biofortification of maize grains. CRC Cr. Rev. Food Sci. 2015, 55, 123-139. (03) Walker, V.; Couillerot, O.; Von Felten, A.; Bellvert, F.; Jansa, J.; Maurhofer, M.; Bally, R.; Moënne-Loccoz, Y.; Comte, G. Variation of secondary metabolite levels in maize seedling roots induced by inoculation with azospirillum, pseudomonas and glomus consortium under field conditions. Plant Soil 2012, 356, 151–163. (04) Narsian, V.; Patel, H. H. Aspergillus aculeatus as a rock phosphate solubilizer. Soil Biol. Biochem. 2000, 32, 559-565. (05) Vitousek, P. M.; Aber, J.; Howarth, R. W.; Likens, G. E.; Matson, P. A.; Schindler, D. W.; Schlesinger, W. H.; Tilman, D. G. Human alteration of the global nitrogen cycle: sources and consequences. Ecol. Appl. 1997, 7, 737–750. (06) Kowaljow, E.; Mazzarino, M. J. Soil restoration in semiarid Patagonia: Chemical and biological response to different compost quality. Soil Biol. Biochem. 2007, 39, 1580–1588. (07) Spaccini, R.; Piccolo, A. Molecular characterization of compost at increasing stages of maturity. 1. Chemical fractionation and infrared spectroscopy. J. Agric. Food Chem. 2007, 55, 2293-2302. (08) Spaccini, R., Piccolo, A. Molecular characteristics of humic acids extracted from compost at increasing maturity stages. Soil Biol. Biochem. 2009, 41, 1164–1172. (09) Martinez-Balmori, D.; Spaccini, R.; Aguiar, N. O.; Novotny, E. H.; Olivares, F. L.; Canellas, L. P. Molecular Characteristics of Humic Acids Isolated from Vermicomposts and Their Relationship to Bioactivity. J. Agric. Food Chem. 2014, 62, 11412−11419. (10) Rohlig, R. M.; Engel, K.H. Influence of the Input System (Conventional versus Organic Farming) on Metabolite Profiles of Maize (Zea mays) Kernels. J. Agric. Food Chem. 2010, 58, 3022–3030. (11) Covacevich, F.; Echeverrıa, H. E.; Aguirrezabal, L. A. N. Soil available phosphorus status determines indigenous mycorrhizal colonization of field and glasshouse-grown spring wheat from Argentina. Appl. Soil Ecol. 2007, 35, 1–9. (12) Smith, S.E.; Read, D. J. Mycorrhizal Symbiosis (Third ed.); Elsevier: London, 2008. (13) Cozzolino, V.; Di Meo, V.; Piccolo, A. Impact of arbuscular mycorrhizal fungi applications on maize production and soil phosphorus availability, J. Geochem. Explor. 2013, 129, 40–44. 21

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560

(14) Mazzei, P.; Piccolo, A.; Valentini, M. Intact Food Analysis by Means of HRMAS-NMR Spectroscopy. In Modern Magnetic Resonance, Webb, G.A. (Ed.) Second edition; Springer International Publishing: Switzerland, 2017, pp 1-16, In press, DOI 10.1007/978-3-319-282756_16-1. (15) Doty, F. D.; Entzminger, G.; Yang, A. Y. Magnetism in high-resolution NMR probe design. II: HR MAS. Concept Magnetic Res. 1998, 10, 239–260. (16) Mazzei, P.; Piccolo, A. HRMAS NMR spectroscopy applications in agriculture. Chem. Biol. Technol. Agric. 2017, 4, 1-11. (17) Mazzei, P.; Piccolo, A.; Nugnes, L.; Mascolo, M.; De Rosa, G.; Staibano, S. Metabolic profile of uterine leiomyomas intact tissue using high resolution magic-angle Spinning 1H NMR spectroscopy. NMR Biomed. 2010, 23, 1–10. (18) Mazzei, P.; Piccolo, A. 1H HRMAS-NMR metabolomic to assess quality and traceability of mozzarella cheese from Campania buffalo milk. Food Chem. 2012, 132, 1620–1627. (19) Mazzei, P.; Vinale, F.; Woo, S. L.; Pascale, A.; Lorito, M.; Piccolo, A. Metabolomics by Proton High-resolution magic-angle-spinning nuclear magnetic resonance of tomato plants treated with two secondary metabolites isolated from trichoderma. J. Agric. Food Chem. 2016, 64, 3538−3545. (20) Kruger, N. J.; Troncoso-Ponce, M. A.; Ratcliffe, R. G. 1H NMR metabolite fingerprinting and metabolomic analysis of perchloric acid extracts from plant tissues. Nat. Prot. 2008, 3, 10011012. (21) Schripsema, J. Application of NMR in plant metabolomics: techniques, problems and prospects. Phytochem. Anal. 2010, 21, 14–21. (22) Wu, X.; Li, N.; Li, H., Tang, H. An optimized method for NMR-based plant seed metabolomic analysis with maximized polar metabolite extraction efficiency, signal-to-noise ratio, and chemical shift consistency. Analyst 2014, 139, 1769-1778. (23) Manetti, C.; Bianchetti, C.; Casciani, L.; Castro, C.; Di Cocco, M. E.; Miccheli, A.; Motto, M.; Conti, F. A. A metabonomic study of transgenic maize (Zea mays) seeds revealed variations in osmolytes and branched amino acids. J. Exp. Bot. 2006, 57, 2613–2625. (24) Broyart, C.; Fontaine, J. X.; Molinie, R.; Cailleu, D.; Terce-Laforgue, T.; Dubois, F.; Hirel, B.; Mesnard, F. Metabolic profiling of maize mutants deficient for two glutamine synthetase isoenzymes using 1H-NMR-based metabolomics. Phytochem Anal. 2010, 21, 102-109. (25) Barros, E.; Lezar, S.; Anttonen, M. J.; van Dijk, J. P.; Rohlig, R. M.; Kok, E. J.; Engel, K. H. Comparison of two GM maize varieties with a near isogenic non-GM variety using transcriptomics, proteomics and metabolomics. Plant Biotechnol. 2010, 8, 436–451. (26) Gavaghan, C. L.; Li, J. V.; Hadfield, S. T.; Hole, S.; Nicholson, J. K.; Wilson, I. D.; Howe, P. W. A.; Stanley, P. D.; Holmes, E. Application of NMR-based metabolomics to the investigation of salt stress in maize (Zea mays). Phytochem. Anal. 2011, 22, 214–224. 22

ACS Paragon Plus Environment

Page 22 of 32

Page 23 of 32

561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607

Journal of Agricultural and Food Chemistry

(27) Kusano, M.; Fukushima, A.; Redestig, H.; Saito, K. Metabolomic approaches toward understanding nitrogen metabolism in plants. J. Exp. Bot. 2011, 62, 1439-1453. (28) Shenton, M. E.; Dickey, C. C.; Frumin, M.; McCarly, R. W. A review of MRI findings in schizophrenia. Schizophr. res. 2001, 49, 1-52. (29) Bollineni, V. R.; Kramer, G.; Liu, Y.; Melidis, C.; de Souza, N. M. A literature review of the association between diffusion-weighted MRI derived apparent diffusion coefficient and tumour aggressiveness in pelvic cancer. Cancer Treat. Rev. 2015, 41, 496–502. (30) Schmidt, S.J.; Sun, X.; Litchfield, J. B. Applications of magnetic resonance imaging in food science. Crit. Rev. Food Sci. Nutr. 1996, 36, 357-85. (31) Van As, H.; Scheenen, T.; Vergeldt, F. J. MRI of intact plants. Photosynth Res. 2009, 102, 213–222. (32) Ruan, R.; Litchfield, J. B.; Eckhoff, S. V. Simultaneous and nondestructive measurements of transient moisture profiles and structural changes in corn kernels during steeping using microscopic nuclear magnetic resonance imaging. Cereal Chem. 1996, 69, 600-606. (33) Van Der Weerd, L.; Claessens, M .M. A. E.; Efdè, C.; Van As, H. Nuclear magnetic resonance imaging of membrane permeability changes in plants during osmotic stress. Plant Cell Environ. 2002, 25, 1539–1549. (34) Borompichaichartkul, C.; Moran, G.; Srzednicki, G. S.; Driscoll, R. H. Studies of physical state of water in maize from Northeast China (cv. Huangmo 417) during drying at subzero temperatures. Drying techn. 2004, 22, 295-305. (35) Janas, S.; Boutry, S.; Malumba, P.; Elst, L. V.; Bera, F. Modelling dehydration and quality degradation of maize during fluidized-bed drying. J. Food Eng. 2010, 100, 527–534. (36) Zilic, S.; Milasinovic, M.; Terzic, D.; Barac, M.; Ignjatovic-Micic, D. Grain characteristics and composition of maize specialty hybrids. Span. J. Agric. Res. 2011, 9, 230-241. (37) Todeschini, R. Introduzione alla Chemiometria; EdiSES s.r.l.: Naples, Italy, 1998. (38) Brereton, R. G. Pattern recognition. In Chemometrics: data Analysis for the Laboratory and Chemical Plant J. Wiley & sons (Eds.), Wiley: England (UK), 2003, pp 183−249. (39) Piccioni, F.; Capitani, D.; Zolla, L.; Mannina, L. NMR Metabolic profiling of transgenic maize with the Cry1A(b) Gene. J. Agric. Food Chem. 2009, 57, 6041–6049. (40) Alonso, P.; Val, D. L.; Yair, S. H. Central metabolic fluxes in the endosperm of developing maize seeds and their implications for metabolic engineering. Metab. Eng. 2011, 13, 96–107.

23

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631

(41) Mazzei, P.; Francesca, N.; Moschetti, G.; Piccolo, A. Metabolomic by 1H NMR spectroscopy differentiates “Fiano di Avellino” white wines obtained with different yeast strains. J. Agric. Food Chem. 2013, 61, 10816−10822.

632 633 634 635 636 637 638 639 640 641 642 643 644 645 646

(47) Borisjuk, L.; Rolletschek, H.; Fuchs, J.; Melkus, G.; Neuberger, T. Low and high field magnetic resonance for in vivo analysis of seeds. Materials 2011, 4, 1426-1439.

(42) Amiour, N.; Imbaud, S.; Clement, G.; Agier, N.; Zivy, M.; Valot, B.; Balliau, T.; Armengaud, P.; Quillere, I.; Canas, R.; Tercet-Laforgue, T.; Hirel, B. The use of metabolomics integrated with transcriptomic and proteomic studies for identifying key steps involved in the control of nitrogen metabolism in crops such as maize. J. Exp. Bot. 2012, 63, 5017-5033. (43) Vaccaro, S.; Muscolo, A.; Pizzeghello, D.; Spaccini, R.; Piccolo, A.; Nardi, S. Effect of a Compost and Its Water-Soluble Fractions on Key Enzymes of Nitrogen Metabolism in Maize Seedlings. J. Agric. Food Chem. 2009, 57, 11267–11276. (44) Stitt, M.; von Schaeven, A.; Willmitzer, L. ‘Sink’ regulation of photosynthetic metabolism in transgenic tobacco plants expressing yeast invertase in their cell wall involves a decrease of the Calvin-cycle enzymes and an increase of glycolytic enzymes. Planta 1990, 183, 40–50. (45) Joner, E.J. The effect of long-term fertilization with organic or inorganic fertilizers on mycorrhiza-mediated phosphorus uptake in subterranean clover. Biol. Fertil. Soil 2000, 32, 435– 440. (46) Gosling, P.; Hodge, A.; Goodlass, G.; Bending, G. D. Arbuscular mycorrhizal fungi and organic farming. Agric. Ecosys. Environ. 2006, 113, 17–35.

(48) Gayral, M.; Bakan, B.; Dalgalarrondo, M.; Elmorjani, K.; Delluc, K.; Brunet, S.; Linossier, L.; Morel, M. H.; Marion, D. Lipid partitioning in maize (Zea mays L.) endosperm highlights relationships among starch lipids, amylose, and vitreousness. J. Agric. Food Chem. 2015, 63, 3551-3558. (49) Nuss, E. T.; Tanumihardjo, S. A. Maize: A paramount staple crop in the context of global nutrition. Compr. Rev. Food Sci. Food 2010, 9, 417-436. (50) Sibgatullin, t. A.; Vergeldt, F. J.; Gerkema, E.; Van As, H. Quantitative permeability imaging of plant tissues. Eur. Biophys. J. 2010, 39, 699-710.

24

ACS Paragon Plus Environment

Page 24 of 32

Page 25 of 32

Journal of Agricultural and Food Chemistry

647 648

Figures Captions

649 650 651 652

Figure 1. 1H HRMAS CPMG NMR spectrum of maize seed (M- sample). The proton spectral regions confined within the two grey rectangles (1 and 2) have been magnified on top sides of the figure. The labels indicate assignment of the most intense signals detected in all spectra. (F, Fructose; G, Glucose; LA, Lactic Acid; S, Sucrose; UMP, Uridinmonophosphate; UDP, Uridindiphosphate).

653 654 655 656 657 658 659

Figure 2. Principal Component Analysis score-plots (A, B, and C) of data obtained from 1H CPMG HRMAS spectra of maize seeds deriving from different treatments (M-, no fertilization; MIN, Mineral fertilization, 250 kg ha-1 of N as urea, 150 kg ha-1 of P2O5 as Ca(H2PO4)2 and 160 kg ha-1 of K as K2SO4; COM, 15 t ha-1 compost fertilization; +, AFM inoculation, mix of R. irregularis and F. mosseae). Each score-plot includes the name and the indicative direction assumed (Figure S1 in Supporting Information) by the most significant (ANOVA Tukey test, confidence level > 95%) loading vectors involved in the differentiation among treatments.

660 661 662 663 664 665

Figure 3. Principal Component Analysis score-plots (A, B, and C) of data obtained from 1H CPMG HRMAS spectra of maize seeds deriving from different treatments (M-, no fertilization; COM, 15 t ha-1 compost fertilization; +, AFM inoculation, mix of R. irregularis and F. mosseae). Each scoreplot includes the name and the indicative direction assumed (Figure S2 in Supporting Information) by the most significant (ANOVA Tukey test, confidence level > 95%) loading vectors involved in the differentiation among treatments.

666 667 668

Figure 4. Sagittal (A) and axial (B and C) MRI images of an integer M- maize caryopsis obtained by using the pulse sequence MSME at a total spin-echo time of 10.57 ms. The axial Figures refer to the slices indicated by rectangles (A) and highlight the floury (B) and vitreous (C) endosperm regions.

669 670 671 672

Figure 5. Coronal MRI images of an integer M- maize caryopsis obtained by using the pulse sequence MSME (A and B) and the Diffusion sequence integrated with stimulated echoes (C, gradient strength of 362.9 s mm-2). MSME experiments were set with two total spin-echo times: 10.57 ms (A, No T2-weighted) and 42.3 ms (B, T2- weighted).

673 674 675 676 677 678 679 680

Figure 6. Principal Component Analysis score-plots (A, B, and C) based on MRI parameters (T1, T2, Diff) calculated for vitreous (EV) and floury (EF) endosperm of intact maize caryopses deriving from different treatments (M-, no fertilization; MIN, Mineral fertilization, 250 kg ha-1 of N as urea, 150 kg ha-1 of P2O5 as Ca(H2PO4)2 and 160 kg ha-1 of K as K2SO4; COM, 15 t ha-1 compost fertilization; +, AFM inoculation, mix of R. irregularis and F. mosseae). Each score-plot includes the name and the indicative direction assumed (Figure S5 in Supporting Information) by the of most significant (ANOVA Tukey test, confidence level > 95%) loading vectors involved in the differentiation among treatments.

681 25

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Figure 1. 1H HRMAS CPMG NMR spectrum of maize seed (M- sample). The proton spectral regions confined within the two grey rectangles (1 and 2) have been magnified on top sides of the figure. The labels indicate assignment of the most intense signals detected in all spectra. (F, Fructose; G, Glucose; LA, Lactic Acid; S, Sucrose; UMP, Uridinmonophosphate; UDP, Uridindiphosphate). 500x300mm (96 x 96 DPI)

ACS Paragon Plus Environment

Page 26 of 32

Page 27 of 32

Journal of Agricultural and Food Chemistry

Figure 2. Principal Component Analysis score-plots (A, B, and C) of data obtained from 1H CPMG HRMAS spectra of maize seeds deriving from different treatments (M-, no fertilization; MIN, Mineral fertilization, 250 kg ha-1 of N as urea, 150 kg ha-1 of P2O5 as Ca(H2PO4)2 and 160 kg ha-1 of K as K2SO4; COM, 15 t ha-1 compost fertilization; +, AFM inoculation, mix of R. irregularis and F. mosseae). Each score-plot includes the name and the indicative direction assumed (Figure S1 in Supporting Information) by the most significant (ANOVA Tukey test, confidence level > 95%) loading vectors involved in the differentiation among treatments. 230x664mm (96 x 96 DPI)

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Figure 3. Principal Component Analysis score-plots (A, B, and C) of data obtained from 1H CPMG HRMAS spectra of maize seeds deriving from different treatments (M-, no fertilization; COM, 15 t ha-1 compost fertilization; +, AFM inoculation, mix of R. irregularis and F. mosseae). Each score-plot includes the name and the indicative direction assumed (Figure S2 in Supporting Information) by the most significant (ANOVA Tukey test, confidence level > 95%) loading vectors involved in the differentiation among treatments. 254x665mm (96 x 96 DPI)

ACS Paragon Plus Environment

Page 28 of 32

Page 29 of 32

Journal of Agricultural and Food Chemistry

Figure 4. Sagittal (A) and axial (B and C) MRI images of an integer M- maize caryopsis obtained by using the pulse sequence MSME at a total spin-echo time of 10.57 ms. The axial Figures refer to the slices indicated by rectangles (A) and highlight the floury (B) and vitreous (C) endosperm regions. 289x209mm (96 x 96 DPI)

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Figure 5. Coronal MRI images of an integer M- maize caryopsis obtained by using the pulse sequence MSME (A and B) and the Diffusion sequence integrated with stimulated echoes (C, gradient strength of 362.9 s mm-2). MSME experiments were set with two total spin-echo times: 10.57 ms (A, No T2-weighted) and 42.3 ms (B, T2- weighted). 289x209mm (96 x 96 DPI)

ACS Paragon Plus Environment

Page 30 of 32

Page 31 of 32

Journal of Agricultural and Food Chemistry

Figure 6. Principal Component Analysis score-plots (A, B, and C) based on MRI parameters (T1, T2, Diff) calculated for vitreous (EV) and floury (EF) endosperm of intact maize caryopses deriving from different treatments (M-, no fertilization; MIN, Mineral fertilization, 250 kg ha-1 of N as urea, 150 kg ha-1 of P2O5 as Ca(H2PO4)2 and 160 kg ha-1 of K as K2SO4; COM, 15 t ha-1 compost fertilization; +, AFM inoculation, mix of R. irregularis and F. mosseae). Each score-plot includes the name and the indicative direction assumed (Figure S5 in Supporting Information) by the of most significant (ANOVA Tukey test, confidence level > 95%) loading vectors involved in the differentiation among treatments. 319x695mm (96 x 96 DPI)

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Graphical abstract 84x47mm (96 x 96 DPI)

ACS Paragon Plus Environment

Page 32 of 32