Metabolite Profiling of a Zinc-Accumulating Rice Mutant - Journal of

Apr 25, 2017 - Breeding crops with high zinc (Zn) density is an effective way to alleviate human dietary Zn deficiencies. We characterized a mutant Li...
0 downloads 0 Views 724KB Size
Subscriber access provided by UB + Fachbibliothek Chemie | (FU-Bibliothekssystem)

Article

Metabolite profiling of a Zn-accumulating rice mutant Yin Wang, Sha Mei, Zhixue Wang, Zhoulei Jiang, zhangshicang Zhu, Jingwen Ding, Dianxing Wu, and Xiao-Li Shu J. Agric. Food Chem., Just Accepted Manuscript • Publication Date (Web): 25 Apr 2017 Downloaded from http://pubs.acs.org on April 26, 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 26

Journal of Agricultural and Food Chemistry

1

Metabolite profiling of a Zn-accumulating rice mutant

2

Yin Wang, Sha Mei#, Zhixue Wang, Zhoulei Jiang, Zhangshicang Zhu, Jingwen Ding, Dianxing

3

Wu, Xiaoli Shu*

4 5

State Key Laboratory of Rice Biology and Key Laboratory of the Ministry of Agriculture for the

6

Nuclear-Agricultural Sciences, Department of Applied Biosciences, Institute of Nuclear

7

Agricultural Sciences, Zhejiang University, Hangzhou 310029, China

8

# equal contribution as the first author

9

*Corresponding author

10

Email: [email protected]

11 12 13 14 15 16 17 18 19 20 21 22 23 24 1

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 2 of 26

25

Abstract:Breeding crops with high Zn-density is an effective way to alleviate human dietary Zn

26

deficiencies. We characterized a mutant Lilizhi (LLZ) accumulating at least 35% higher Zn

27

concentration in grain than wild-type in hydroponics experiments. The mutant stored less Zn

28

contents in the roots and transported more Zn to the grain. Metabolite profiling demonstrated that

29

with high Zn treatment, the content of proline, asparagine, citric acid and malic acid were

30

enhanced in both LLZ and the WT, which were thought to be involved in Zn transport in rice.

31

Furthermore, the contents of cysteine, allothreonine, alanine, tyrosine, homoserine, β-alanine and

32

nicotianamine that required for the production of many metal-binding proteins were specifically

33

increased in LLZ. LLZ had higher capability of amino acid biosynthesis and metal cation

34

transportation. The current research extends our understanding on the physiological mechanisms

35

of Zn uploading into grain, and provides references for further Zn biofortification breeding in rice.

36

Keywords: :Zinc (Zn); Grain, Amino acid; Nicotianamine; Rice; Biofortification

37 38 39 40 41 42 43 44 45 46 47 48

2

ACS Paragon Plus Environment

Page 3 of 26

Journal of Agricultural and Food Chemistry

49 50

Introduction

51

Zinc is an essential micronutrient for human health, which plays critical roles in diverse

52

biochemical processes, including catalyzing enzymes, facilitating protein folding and regulating

53

gene expression1. The Zn deficiency-induced malnutrition constitutes a major public health

54

problem in the world, affecting more than 25% of the world’s population2. The serious Zn

55

deficiency is a leading risk factor for disability and death, especially for infants, children and

56

pregnant women in developing countries where the cereal-based foods with low concentrations of

57

Zn represent the main diet component3. Thus, breeding new crop cultivars with high grain-Zn

58

concentrations will be an important strategy to alleviate human Zn deficiency.

59

Rice is one of the most important cereals worldwide, providing energy and micronutrients for

60

at least 3 billion people. Unfortunately, rice grain is especially low in Zn content. Among 939

61

rice cultivars analyzed4, Zn concentration in brown rice ranged from 15.9 to 58.4 mg/kg with a

62

mean of 25.4 mg/kg. Due to the narrow range of Zn level in existing germplasm collections, rice

63

seems to be more difficult to make substantial increases in grain Zn content compared to maize

64

and wheat5. Owing to this disadvantageous effect, various efforts have been made to exploit novel

65

Zn-dense germplasm in rice. Jeng et al.6 developed three high Zn-density mutant lines from the

66

rice cultivar IR64, which accumulated considerably higher levels of Zn in their grains (26.58,

67

28.95 and 26.16 mg kg-1, respectively) than the wild type (16.00 mg kg-1). Some transgenic rice

68

lines, which were generated by over-expression of OsNAS genes, could also increase grain Zn

69

concentration7-12. However, none of the above studies has effectively increased the level of Zn

70

content in rice grain to meet the nutritional demand.

71

Due to alkaline pH and high orthophosphate concentration in the phloem, Zn is thought to be

72

transported as complexes with ligands via the phloem while not as free ions13-15. In recent years,

73

advances in high-throughput metabolomics analysis have assisted in studying the multiple

74

physiological processes of Zn homeostasis. A variety of Zn-binding ligands have been identified in

75

plants,

76

metallothioneins16-24. Although a lot of metabolites have been identified and they play diverse

including

organic

acids,

amino

acids,

phytochelatins,

nicotianamine

and

3

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 4 of 26

77

roles in Zn transport, a broad understanding of the metabolic changes in Zn uploading into

78

developing rice grain remains unclear.

79

In this study, a rice mutant, designated as LLZ, which accumulated a large amount of Zn in

80

grains were characterized. We compared the accumulation and distribution of Zn between plant

81

tissues under three Zn conditions in the mutant and wild-type, also we analyzed the metabolic

82

changes in the developing grains of both varieties responding to different Zn treatment. The

83

primary objective of this work is to provide novel metabolic evidence regarding Zn uploading into

84

grain and to assist in elucidating the types of mutation resulting in the enhancement of Zn. The

85

results revealed may lay a foundation for elucidation of the molecular mechanisms associated with

86

Zn transport in rice and digging out genes associated to Zn accumulation in the future and further

87

to develop rice varieties with high Zn concentrations in grain.

88

Materials and Methods

89

Plant materials and growth conditions

90

The rice mutant, Lilizhi (LLZ), was derived from Oryza sativa ssp. Japonica variety

91

Dongbeixiang (DBX) with irradiated by 250 Gy gamma rays (M11 generation). The seeds of both

92

LLZ and WT were surface-sterilized by 3% H2O2 for 20 min and rinsed thoroughly in deionized

93

water, then soaked in deionized water for 48 hours at 30°C. The seeds were transferred to the

94

germination boxes in a growth chamber under photo flux density of 400 µmol m-2 s-1, light/dark

95

period of 16/8 h, day/night temperature of 30/25 °C, and relative humidity of 75-85%. After

96

supplied with only deionized water for the first two weeks, the seedlings were transplanted into

97

48-well plastic buckets (35 L) and concentration of nutrient solution was increased gradually, from

98

quarter of the final concentration for the first three days to half concentration for the next four

99

days, then to the final concentration from the second week. The composition of nutrient solution

100

was the same as that used by Wu et al. 25 except the Zn concentration varied: 1.5 mmol/L NH4NO3,

101

1.0 mmol/L CaCl2, 1.6 mmol/L MgSO4, 1.0 mmol/L, K2SO4, 0.3 mmol/L KH2PO4, 0.05 µmol/L

102

H3BO3, 5.0 µmol/L MnSO4, 0.2 µmol/L CuSO4, 0.5 µmol/LZnSO4, 0.05 µmol/L (NH4)6Mo7O24

103

and 20 µmol/L Na2FeEDTA. The pH of the hydroponic solution was adjusted to 6.0 using 1 M

104

NaOH or HCl. Hydroponics experiments were carried out in a greenhouse facility of Zhejiang 4

ACS Paragon Plus Environment

Page 5 of 26

Journal of Agricultural and Food Chemistry

105

University with natural light and temperature.

106

Treatments and sampling

107

The plants were exposed to three levels of Zn concentration (0.5, 5 and 50 µmol/L). Each

108

treatment was triplet. Whole plant was harvested after seed maturation and separated into grain,

109

glume, rachis, flag leaf, lower stems and leaves, and roots. Lower stems and leaves were collected

110

from the basal leaves and stem internodes in rice. All samples were washed with deionized water

111

to remove superficial nutrient solution. Roots and flag leaves were submerged in a 1.0 L bath

112

containing 1.0 mM LaCl3 and 0.05 mM CaCl2 for 10 min in order to remove apoplastic Zn. Then

113

all tissues were first dried at 105 °C for 30 min and finally dried to constant mass at 70 °C. The

114

dry weight of all tissues was recorded. Dry materials of plant tissues were grounded in the agate

115

mortar, passed through a 60 mesh and stored for elemental analysis. For metabolite profiling

116

analysis, the grain-tissues of LLZ and WT were sampled for 0.5 and 50µmol/L Zn treatment at 15

117

days after flowering (DAF) according to Wu et al. 25 that Zn concentration of rice grains reached

118

the highest level at the stage. Six biological replicates were performed for each sample. All

119

samples were fresh-frozen in liquid nitrogen and stored at -80 °C till metabolite extraction.

120

Elemental analysis

121

The Zn contents of plant tissues were measured as described by Sasaki et al.26. Samples of

122

100±0.01 mg dried flour were hydrolyzed with concentrated HNO3 (60%) at 140 °C. The Zn

123

concentration in the digest solution was determined by inductively coupled plasma-mass

124

spectrometry (ICP-MS 7700X; Agilent Technologies). The key parameters for ICP-MS were

125

plasma RF-power 1500 W, carrier gas flow rate 0.9 L/min, makeup gas flow rate 0.17 L/min,

126

sample depth 8 mm, nebulizer pump seed 0.1 rps and spray chamber temperature 2°C. Each

127

sample was triplet and Student’s t-test was conducted with SPSS 20.0 (IMB Corp. Chicago).

128

Metabolite extraction and metabolite profiling analysis

129

Metabolites were extracted and derivatized according to Lisec et al.

27

. Briefly, frozen

130

grain-tissues samples (50 mg) were homogenized in 700 µL of 100% methanol containing 30 µL

131

of adonitol (2 mg/mL in water) as internal standard. After centrifugation at 11000 g for 15 min, the

132

supernatants were transferred to a fresh 2 mL glass vial. Afterwards, 375 µL of chloroform and 5

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 6 of 26

133

750 µL H2O were added to the samples. After centrifugation at 2200 g for 15 min, 150 µL of upper

134

phases were taken and dried in a vacuum concentrator without heating. Then, 40 µL of

135

methoxyaminhydrochloride (20 mg/mL in pyridine) was added to the dried samples and shaken

136

for 2 h at 37°C. Afterwards, 70 µL of the MSTFA mix (1 ml+20 µL FAMEs) was added to the

137

sample and incubated for 30 min at 37°C. The derivatized samples were then analyzed by a gas

138

chromatograph system coupled with a time-of-flight mass spectrometer (Pegasus HT GC-TOF-MS;

139

LECO). The program of temperature-rise was followed by initial temperature 80 °C for 1 min,

140

10 °C min−1 rate up to 295 °C, then 295 °C for 6 min. The mass spectrometry data were acquired

141

in full-scan mode with the m/z range of 33-600 at a rate of 20 spectra per second after a solvent

142

delay of 300 s. The mixed aliquots of all prepared sample extracts were included into the

143

measurement for quality control. The raw peaks exacting, the data baselines filtering, peak

144

alignment and deconvolution analysis were carried out by Chroma TOF software of LECO

145

Corporation. Peaks were identified and annotated by comparison with retention time indices (RIs)

146

and mass spectra from the LECO/Fiehn metabolomics library databases. The total mass of signal

147

integration area was normalized on the basis of the fresh weight of the sample and the added

148

amount of an internal standard (adonitol).

149

Statistical analysis and metabolic pathway construction

150

The normalized data involving the peak number, sample name, and normalized peak areas

151

was fed to SIMCA14 software package for multivariate statistical analyses comprising of principal

152

component analysis (PCA) and orthogonal projections to latent structures-discriminate analysis

153

(OPLS-DA). Data was log transformed prior to PCA and OPLS-DA. Differences in the metabolite

154

levels were detected using the Limma package in R software platform

155

metabolite showed a significant change (p≤0.05) at a threshold of VIP values greater than 1.0, it

156

was defined as the specific metabolite. Benjamini and Hochberg correction method was also used

157

to account for multiple testing29. The false discovery rate (FDR) for each p-value threshold was

158

estimated by calculating q-values, based on the P-values derived from Limma analysis. Metabolic

28

. Once the level of a

6

ACS Paragon Plus Environment

Page 7 of 26

Journal of Agricultural and Food Chemistry

159

pathway was constructed according to KEGG metabolic database (http://www.genome.jp/kegg/).

160 161

Results

162

Plant growth and Zn accumulation

163

The mutant LLZ exhibited a high Zn density in grain, and the Zn concentration in brown rice

164

of LLZ reached to 52.7 mg/kg in the fields, but only 24.4 mg/kg in WT (Table S1). For comparison

165

of Zn uptake and accumulation in the two genotypes, the mutant LLZ and WT were grown in

166

nutrient solution. When the plants exposed to three Zn concentrations (0.5, 5 and 50 µM), both of

167

them grew healthy and show no symptom of Zn toxicity. Biomass productions of both LLZ and

168

WT were not significantly affected by enhanced Zn in solution (Table 1). While the biomass of

169

shoot and whole-plant of WT was larger than that of the mutant in each hydroponic experiment,

170

but there was no significant difference in root biomass between two varieties (Table 1).

171

With elevated Zn supply, the Zn concentrations increased whether in whole plant or different

172

tissues in both LLZ and WT (Table 1 and Fig 1). Although there was no difference in Zn

173

concentration of whole-plant between LLZ and WT under all three treatments (Table 1), the Zn

174

concentrations in different plant tissues between LLZ and WT varied greatly (Fig 1). Consistent

175

with the previous field data, the average Zn levels of brown rice in LLZ were 1.63-fold, 1.65-fold

176

and 1.35-fold, higher than those in WT under 0.5, 5 and 50 µM Zn conditions respectively (Fig 1).

177

In companied with the Zn accumulation in grains, the average Zn concentration in glume of LLZ

178

was 31.1% higher than that in WT under three Zn conditions (Fig 1). However, the LLZ and WT

179

did not show significant difference in the Zn accumulation of flag leaf and rachis; Furthermore,

180

the mutant LLZ accumulated 34.0% and 34.8% lower Zn concentrations in root and lower stem

181

respectively than those of WT on average of the three Zn conditions (Fig 1). The Zn

182

concentrations of lower leaves between LLZ and WT differed significantly when exposed to 0.5

183

and 5 µM Zn treatments, whereas it showed a similar accumulation pattern in both varieties under

184

50 µM Zn condition (Fig 1). Thus, the mutant displayed lower concentrations of Zn in roots and

185

higher concentrations of Zn in shoots, the shoot to root ratio of Zn concentration (S/R) in LLZ was

186

at least 1.2-fold greater than that in WT under each Zn treatment (Table 1). 7

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

187

Page 8 of 26

Comparison of metabolite levels between LLZ and WT

188

To reveal physiological mechanisms of Zn uploading into rice seed, the metabolite profiles in

189

the grains of LLZ and WT subjected to 0.5 and 50 µM Zn treatments were investigated. With the

190

grains at 15 DAF, 136 metabolites were identified according to the reporting guidelines by Fernie

191

et al.30 among a total of 695 intensive signals (Table S2). These metabolites included amino acids,

192

organic acids, carbohydrates, lipids and other secondary metabolites. The raw GC-TOF-MS

193

metabolomics data were available at the MetaboLights Web site (accession MTBLS448).

194

Principal component analysis (PCA) revealed obvious metabolic differentiation between

195

samples within treatments and genotypes (Fig. 2). The first two principle components can explain

196

46.8% of the total variance. The two genotypes can be separated by the PC1 which represented

197

28.1% of the total detected metabolite variance and the samples with different Zn treatment was

198

clustered along PC2. Among the corresponding loadings, nicotianamine was found to be one of

199

the major contributors for the separation along PC1, organic acids and amino acids, such as citric

200

acid, glutamine and cystine were the dominate metabolites contributors to PC2 (Table S3).

201

To determine the specific different metabolomes between varieties, limma package was

202

further used to dissect the levels of metabolites between the grains of two varieties under 0.5 and

203

50 µM Zn concentrations respectively. The results indicated that there were 29 and 62 metabolites

204

with a significant change between LLZ and WT with the 0.5 µM and 50 µM Zn treatment,

205

respectively (Table S4 and Table S5). And 24 differential metabolites were consistently present in

206

two treatments (Table 2), among which 17 metabolites with higher contents in LLZ and 7

207

metabolites with lower contents in LLZ when compared to WT (Table 2). The consistently

208

accumulated metabolites in LLZ were majorly amino acids and organic acids, including glycine,

209

glutamic acid, asparagine, tryptophan, glutamine, glutaconic acid and pimelic acid, the contents of

210

nicotianamine in LLZ with 0.5 and 50 µM Zn treatment were respectively 1.3-fold and 1.9-fold

211

higher than those in the WT (Table 2), whereas the amounts of 3 carbohydrates (gentiobiose,

212

allose and glucose-6-phosphate) in the mutant were significantly decreased than that in the WT

213

(Table 2).

214

Metabolic changes of LLZ and WT in response to high Zn treatment 8

ACS Paragon Plus Environment

Page 9 of 26

Journal of Agricultural and Food Chemistry

215

To get a comprehensive view for metabolites changes in response to high Zn treatment, we

216

compared the metabolomes of samples with the treatment of 50 µM Zn concentrations to those

217

treated with 0.5 µM Zn. The analysis revealed that the influence of high Zn treatment on

218

metabolism was highly dependent on genotypes. There were 52 and 42 metabolites with a

219

significant change in mutant LLZ and WT, respectively. Among which, 27 metabolites showed

220

identical responses to high Zn treatment in both two genotypes (Fig 3).Among these identical

221

metabolites, 5 amino acids (proline, glycine, asparagine, isoleucine and glutamic acid) 8 organic

222

acids such as threonic acid, pimelic acid, glyceric acid, α-ketoglutaric acid, citric acid, aconitic

223

acid, malic acid and fumaric acid, and 1 carbohydrates (2-deoxy-D-glucose) showed significant

224

increment, and 3 amino acids (glutamine, tryptophan and cystine), 3 organic acids (glyceric

225

acid-3-phosphate, pyruvic acid and succinic acid), and 4 carbohydrates (fructose-6-phosphate,

226

glucose-6-phosphate, fructose-2,6-bisphosphate and glucose-1-phosphate) were significantly

227

down-accumulated in two genotypes(Fig 3).

228

Compared with WT, there were 25 metabolites showing variety-specific responses to high Zn

229

treatment in LLZ (Fig 3). Organic acid and amino acid metabolism appeared to be affected by high

230

Zn treatment in LLZ, resulting in significantly accumulation of 7 amino acids and 7 organic acids,

231

including valine, alanine, tyrosine, tartronic acid and glycolic acid, and down-regulating of

232

β-glutamic acid, 3-hydroxypropionic acid and 2,2-dimethylsuccinic acid. In addition, the level of

233

allose was significantly increased by 1.1-fold in response to high Zn treatment in LLZ (Fig 3). The

234

significant increase of ethanolamine content was also observed in high Zn treated LLZ, indicating

235

that amine biosynthesis was dramatically enhanced (Fig 3). It is noteworthy that nicotianamine

236

and 5-methylthioadenosine were also dramatically enhanced in LLZ, which were biosynthesized

237

from S-adenosyl methionine (Fig 3).

238

Meantime, the contents of 2-ketobutyric acid and cysteine exhibited different responses to

239

high Zn treatment between LLZ and WT, which were significantly down-regulated in WT but

240

up-accumulated in the mutant LLZ (Fig 3). Xylose, which was one of the important carbohydrates

241

in plant, increased by 0.5-fold in the WT but decreased by 0.6-fold in the mutant (Fig 3).

242

Discussion 9

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

243

Page 10 of 26

Tissue-specific accumulation of Zn in LLZ

244

In this study, we characterized a novel Zn-dense mutant LLZ, which showed high efficiency

245

of Zn enrichment in grain. Impa et al. 31 proposed that continued root uptake can fully account for

246

seed Zn allocation under Zn-sufficient conditions. Zn concentration in the whole-plant reflected

247

the absorptive capacity of root uptake32. However, the whole-plant Zn concentrations between the

248

mutant LLZ and WT did not differ significantly at any Zn levels (Fig 1), indicating they have

249

similar absorptive capacity of root uptake.

250

Wu et al. 25 found that large amounts of the Zn in rice seeds was accumulated through

251

remobilization from sources tissues rather than continued root uptake during grain filling. Zinc in

252

the grains may be actively supplied via the phloem after mobilization from the blades of the flag

253

and upper leaves33. The remobilization of Zn from vegetable tissues to grain was major source of

254

grain Zn loading under Zn-deficient conditions31. The mutant LLZ accumulated less Zn in root and

255

the lower stem while more Zn in upper shoot tissues, especially in glumes and grains, presented

256

different root-to-shoot partitioning of Zn from WT (Fig 1), characterized as a higher ratio of

257

shoot/root Zn concentration (Table 1). If the great accumulation of Zn in LLZ grain was majorly

258

transported from the upper leaves as previous reported34, 35, the Zn concentration in the upper leaf

259

tissues of the LLZ would be lower than that in WT. However, there was no significant different in

260

the Zn accumulation of flag leaf and rachis between LLZ and WT (Fig 1). That revealed the

261

mutant LLZ transported more Zn to upper shoot tissues from root and had high efficiency in

262

root-to-shoot Zn transporting which may mainly contribute to the higher Zn accumulation of grain

263

in LLZ.

264

Metabolic changes induced by high Zn treatment

265

Metabolites can reflect the responses of plant to certain environment. Excessive Zn stress

266

caused similar changes of metabolic pathways in LLZ and WT. The glycolysis appeared to be

267

inhibited

268

glucose-6-phosphate, 3-phosphoglycerate and pyruvic acid in both varieties (Fig 3). What’s more,

269

the TCA cycle, associated with glycolysis pathway, was also be affected under high Zn condition

270

in both varieties, shown as higher contents of citric acid, aconitic acid, fumaric acid and malic acid

271

(Fig 3). That suggested Zn might affect the energy production system of rice. Some organic acids

by

Zn

stress,

exhibited

as

decreased

contents

of

fructose-6-phosphate,

10

ACS Paragon Plus Environment

Page 11 of 26

Journal of Agricultural and Food Chemistry

272

involved in TCA cycle, like citric acid and malic acid, were considered as chelator and associated

273

with the transportation of zinc/cadmium in Thlaspi caerulescens, Arabidopsis halleri, and Triticum

274

aestivuml 21, 22, 36. Citric acid and malic acid increased significantly in two varieties under high Zn

275

treatment (Fig 3), being favorable for binding partners of zinc ion into xylem and facilitating the

276

long-distance transport of Zn.

277

Synthesis of some specific amino acids was significantly enhanced when plant were exposed

278

to toxic metals. Amino acids, including proline, isoleucine, glutamine, glycine, and asparagine,

279

were significantly increased in two varieties with high Zn exposure (Fig 3). Proline has been

280

shown to increase rapidly in response to heavy metal exposure16,

281

previously considered as one of candidate ligands for heavy metals during xylem transport40, 41.

282

That the large amounts of proline and asparagine accumulated in response to high Zn treatment

283

was agreement with previous observations (Fig 3). However, the contents of some amino acids,

284

such as serine, threonine, methionine and aspartic acid, were not significantly affected in both

285

varieties by high Zn stress (Fig 3). That might because the buffering effects through the

286

modulation of down-stream amino acids.

287

Varieties dependent metabolic changes induced by high Zn treatment

37-39

; while asparagine was

288

The differential responses of some metabolic biomarkers to high Zn treatment between LLZ

289

and WT indicated that the influences of Zn on metabolites varied between genotypes. LLZ

290

exhibited more active amino acid biosynthesis in response to high Zn stress. Amino acids were the

291

precursors and constituents of proteins. The specific accumulated amino acids, like cysteine,

292

allothreonine, alanine, tyrosine, valine, homoserine and β-alanine in LLZ (Fig 3), were necessary

293

for the production of many metal-binding proteins which participated in the transportation of

294

metal cations, including phytochelatins and metallothioneins18, 42-47. Nicotianamine was thought to

295

be another key ligand involved in long- and short-distance transport of metal cations19. The

296

dramatic enhancement of nicotianamine in grains of LLZ irrespective of Zn treatments might

297

increase the mobility of Zn as the chemical forms of Zn in phloem (Table 2 and Fig 3). Genetic

298

regulation of nicotianamine synthases (NAS), which is responsible for biosynthesis of

299

nicotianamine, could enhance Zn concentration greatly. The Zn concentrations in seed of

300

transgenic tobacco lines constitutively expressing a barley NAS gene increased 1.8-fold48. For rice, 11

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 12 of 26

301

constitutive over-expression of genes encoding nicotianamine synthase also led to elevated

302

nicotianamine levels, greater total content of Zn in grain and better plant tolerance to Zn

303

deficiency7-10. Furthermore, the endosperm of grains from transgenic rice over-expressing NAS

304

genes contained more bio-available Zn9. In reverse, suppression of NAS genes in Arabidopsis

305

halleri rendered a concomitant decrease in root-to-shoot translocation of Zn49, 50. Therefore, the

306

significant higher content of nicotianamine and amino acids in LLZ may contribute to the higher

307

Zn accumulation in LLZ.

308

Enriching Zn in grain is an important goal of micronutrient biofortification in rice. The

309

transportations of zinc from root to shoot and further to grain are the results of interplays among

310

not only different membrane transporters but also the mobility of zinc binding to ligands within

311

intercellular and organs. LLZ accumulated significant higher concentration of Zn in grains due to

312

improving root-to-shoot transport of Zn. Based the different metabolic profiles between the LLZ

313

and WT under normal and high Zn conditions, we proposed the enhancement of metal chelators

314

such as nicotianamine in LLZ was favor in accumulating Zn in grains. The characterization of a

315

novel high Zn accumulation rice mutant LLZ would help for further exploring the physiological

316

mechanisms of Zn uploading into grain, dig out the related genes combined other omics data, and

317

then provide molecular basis for Zn biofortification breeding in rice.

318 319

Acknowledgements

320

The authors gratefully acknowledge the financial support for this work from Zhejiang

321

Provincial Department of Science and Technology (2016YDF101801, 0406) and HarvestPlus

322

Program (2014H8327.ZHU).

323 324 325

Supporting Information Available: Additional Tables. This material is available free of charge via the Internet at http://pubs.acs.org.

326 327 328 329 12

ACS Paragon Plus Environment

Page 13 of 26

Journal of Agricultural and Food Chemistry

330

References

331

1.

768-772.

332 333

2.

3.

Welch, R. M.; Graham, R. D. A new paradigm for world agriculture: meeting human needs Productive, sustainable, nutritious. Field Crop Res. 1999, 60, 1-10.

336 337

Maret, W.; Sandstead, H. H. Zinc requirements and the risks and benefits of zinc supplementation. J. Trace Elem. Med. Biol. 2006, 20, 3-18.

334 335

Saper, R. B.; Rash, R. Zinc: An Essential Micronutrient. Am. Fam. Physician. 2009, 79,

4.

Graham, R.; Senadhira, D.; Beebe, S.; Iglesias, C.; Monasterio, I. Breeding for micronutrient

338

density in edible portions of staple food crops: conventional approaches. Field Crop Res.

339

1999, 60, 57-80.

340

5.

resources perspective. Food Chem. 2003, 80, 589-596.

341 342

Kennedy, G.; Burlingame, B. Analysis of food composition data on rice from a plant genetic

6.

Jeng, T. L.; Lin, Y. W.; Wang, C. S.; Sung, J. M. Comparisons and selection of rice mutants

343

with high iron and zinc contents in their polished grains that were mutated from the indica

344

type cultivar IR64. J. Food Compos. Anal. 2012, 28, 149-154.

345

7.

Johnson, A. A.; Kyriacou, B.; Callahan, D. L.; Carruthers, L.; Stangoulis, J.; Lombi, E.;

346

Tester, M. Constitutive overexpression of the OsNAS gene family reveals single-gene

347

strategies for effective iron- and zinc-biofortification of rice endosperm. PloS one. 2011, 6,

348

e24476.

349

8. Lee, S.; Jeon, U. S.; Lee, S. J.; Kim, Y. K.; Persson, D. P.; Husted, S.; Schjorring, J. K.; Kakei,

350

Y.; Masuda, H.; Nishizawa, N. K.; An, G. Iron fortification of rice seeds through activation of

351

the nicotianamine synthase gene. P Natl Acad Sci USA. 2009, 106, 22014-22019. 13

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

352

9.

Page 14 of 26

Lee, S.; Persson, D. P.; Hansen, T. H.; Husted, S.; Schjoerring, J. K.; Kim, Y. S.; Jeon, U. S.;

353

Kim, Y. K.; Kakei, Y.; Masuda, H.; Nishizawa, N. K.; An, G. Bio-available zinc in rice seeds

354

is increased by activation tagging of nicotianamine synthase. Plant Biotechnol J. 2011, 9,

355

865-873.

356

10. Masuda, H.; Usuda, K.; Kobayashi, T.; Ishimaru, Y.; Kakei, Y.; Takahashi, M.; Higuchi, K.;

357

Nakanishi, H.; Mori, S.; Nishizawa, N. K. Overexpression of the Barley Nicotianamine

358

Synthase Gene HvNAS1 Increases Iron and Zinc Concentrations in Rice Grains. Rice. 2009, 2,

359

155-166.

360

11. Wirth, J.; Poletti, S.; Aeschlimann, B.; Yakandawala, N.; Drosse, B.; Osorio, S.; Tohge, T.;

361

Fernie, A. R.; Gunther, D.; Gruissem, W.; Sautter, C. Rice endosperm iron biofortification by

362

targeted and synergistic action of nicotianamine synthase and ferritin. Plant Biotechnol J.

363

2009, 7, 631-644.

364

12. Zheng, L. Q.; Cheng, Z. Q.; Ai, C. X.; Jiang, X. H.; Bei, X. S.; Zheng, Y.; Glahn, R. P.; Welch,

365

R. M.; Miller, D. D.; Lei, X. G.; Shou, H. X. Nicotianamine, a Novel Enhancer of Rice Iron

366

Bioavailability to Humans. PloS one. 2010, 5, 7.

367 368 369 370 371

13. Fukumorita, T.; Chino, M. Sugar, amino acid and inorganic contents in rice phloem sap. Plant Cell Physiol. 1982, 23, 273-283. 14. Fukumorita, T.; Noziri, Y.; Haraguchi, H.; Chino, M. Inorganic content in rice phloem sap. Soil Sci. Plant Nutr. 1983, 29, 185-192. 15. Nishiyama, R.; Tanoi, K.; Yanagisawa, S.; Yoneyama, T. Quantification of zinc transport via

372

the phloem to the grain in rice plants (Oryza sativa L.) at early grain-filling by a combination

373

of mathematical modeling and Zn-65 tracing. Soil Sci. Plant Nutr. 2013, 59, 750-755. 14

ACS Paragon Plus Environment

Page 15 of 26

Journal of Agricultural and Food Chemistry

374 375 376 377 378 379 380 381 382

16. Bassi, R.; Sharma, S. S. Changes in proline content accompanying the uptake of Zinc and Copper by Lemna minor. Ann. Bot. 1993, 72, 151-154. 17. Bottari, E.; Festa, M. R. Asparagine as a ligand for cadmium(II), lead(II) and zinc(II). Chem. Speciation Bioavail. 1996, 8, 75-83. 18. Clemens, S. Evolution and function of phytochelatin synthases. J Plant Physiol. 2006, 163, 319-32. 19. Clemens, S.; Deinlein, U.; Ahmadi, H.; Horeth, S.; Uraguchi, S. Nicotianamine is a major player in plant Zn homeostasis. Biometals. 2013, 26, 623-632. 20. Curie, C.; Cassin, G.; Couch, D.; Divol, F.; Higuchi, K.; Jean, M.; Misson, J.; Schikora, A.;

383

Czernic, P.; Mari, S. Metal movement within the plant: contribution of nicotianamine and

384

yellow stripe 1-like transporters. Ann. Bot. 2009, 103, 1-11.

385

21. Salt, D. E.; Prince, R. C.; Baker, A. J. M.; Raskin, I.; Pickering, I. J. Zinc ligands in the metal

386

hyperaccumulator Thlaspi caerulescens as determined using X-ray absorption spectroscopy.

387

Environ Sci Technol. 1999, 33, 713-717.

388

22. Sarret, G.; Saumitou-Laprade, P.; Bert, V.; Proux, O.; Hazemann, J. L.; Traverse, A. S.;

389

Marcus, M. A.; Manceau, A. Forms of zinc accumulated in the hyperaccumulator Arabidopsis

390

halleri. Plant Physiol. 2002, 130, 1815-1826.

391

23. Sharma, S. S.; Dietz, K. J. The significance of amino acids and amino acid-derived molecules

392

in plant responses and adaptation to heavy metal stress. J Exp Bot. 2006, 57, 711-726.

393

24. Tennstedt, P.; Peisker, D.; Bottcher, C.; Trampczynska, A.; Clemens, S. Phytochelatin

394

Synthesis Is Essential for the Detoxification of Excess Zinc and Contributes Significantly to

395

the Accumulation of Zinc. Plant Physiol. 2009, 149, 938-948. 15

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

396

Page 16 of 26

25. Wu, C. Y.; Lu, L. L.; Yang, X. E.; Feng, Y.; Wei, Y. Y.; Hao, H. L.; Stoffella, P. J.; He, Z. L.

397

Uptake, Translocation, and Remobilization of Zinc Absorbed at Different Growth Stages by

398

Rice Genotypes of Different Zn Densities. J Agr Food Chem. 2010, 58, 6767-6773.

399

26. Sasaki, A.; Yamaji, N.; Yokosho, K.; Ma, J. F. Nramp5 Is a Major Transporter Responsible

400

for Manganese and Cadmium Uptake in Rice. The Plant cell. 2012, 24, 2155-2167.

401

27. Lisec, J.; Schauer, N.; Kopka, J.; Willmitzer, L.; Fernie, A. R. Gas chromatography mass

402

spectrometry-based metabolite profiling in plants. Nat Protoc. 2015, 10, 1457-1457.

403

28. Smyth, G. K. Linear models and empirical bayes methods for assessing differential

404

expression in microarray experiments. Stat Appl Genet Mol. Biol. 2004, 3, 1-25

405 406

29. Benjamini Y, Hochberg Y. Controlling the false discovery rate: apractical and powerful approach to multiple testing. J R Stat Soc B. 1995, 57, 289-300

407

30. Fernie, A.R.; Aharoni, A.; Willmitzer, L.; Stitt, M.; Tohge, T.; Kopka, J.; Carroll, A.J.; Saito,

408

K.; Fraser, P.D.; DeLuca, V. Recommendations for reporting metabolite data. Plant Cell.

409

2011, 23, 2477-2482

410

31. Impa, S. M.; Morete, M. J.; Ismail, A. M.; Schulin, R.; Johnson-Beebout, S. E. Zn uptake,

411

translocation and grain Zn loading in rice (Oryza sativa L.) genotypes selected for Zn

412

deficiency tolerance and high grain Zn. J Exp Bot. 2013, 64, 2739-2751.

413

32. Jiang, W.; Struik, P. C.; Lingna, J.; van Keulen, H.; Ming, Z.; Stomph, T. J. Uptake and

414

distribution of root-applied or foliar-applied Zn-65 after flowering in aerobic rice. Ann Appl

415

Biol. 2007, 150, 383-391.

416

33. Yoneyama, T.; Gosho, T.; Kato, M.; Goto, S.; Hayashi, H. Xylem and phloem transport of Cd,

417

Zn and Fe into the grains of rice plants (Oryza sativa L.) grown in continuously flooded 16

ACS Paragon Plus Environment

Page 17 of 26

Journal of Agricultural and Food Chemistry

418 419 420 421

Cd-contaminated soil. Soil Sci. Plant Nutr. 2010, 56, 445-453. 34. Stomph, T. J.; Jiang, W.; Struik, P. C. Zinc biofortification of cereals: rice differs from wheat and barley. Trends Plant Sci. 2009, 14, 123-124. 35. Sperotto, R. A.; Ricachenevsky, F. K.; Waldow, V. D.; Muller, A. L. H.; Dressler, V. L.; Fett, J.

422

P. Rice grain Fe, Mn and Zn accumulation: How important are flag leaves and seed number?

423

Plant Soil Environ. 2013, 59, 262-266.

424

36. Gramlich, A.; Tandy, S.; Frossard, E.; Eikenberg, J.; Schulin, R. Availability of Zinc and the

425

Ligands Citrate and Histidine to Wheat: Does Uptake of Entire Complexes Play a Role ? J

426

Agr Food Chem. 2013, 61, 10409-10417.

427 428 429 430 431 432 433 434 435

37. Callahan, D. L.; Baker, A. J.; Kolev, S. D.; Wedd, A. G. Metal ion ligands in hyperaccumulating plants. J Biol Inorg Chem. 2006, 11, 2-12. 38. Rai, V. K. Role of amino acids in plant responses to stresses. Biol Plantarum. 2002, 45, 481-487. 39. Srinivas, V.; Balasubramanian, D. Proline is a protein-compatible hydrotrope. Langmuir. 1995, 11, 2830-2833. 40. Smirnoff, N.; Stewart, G. R. Nitrogen assimilation and Zinc toxicity to Zinc-tolerant and non-tolerant clones of Deschampsia Cespitosa (L) Beauv. New Phytol. 1987, 107, 671-680. 41. White-Monsant, A. C.; Tang, C. Organic acids are not specifically involved in the

436

nitrate-enhanced Zn hyperaccumulation mechanism in Noccaea caerulescens. Environ Exp

437

Bot. 2013, 91, 12-21.

438 439

42. Benatti, M. R.; Yookongkaew, N.; Meetam, M.; Guo, W. J.; Punyasuk, N.; AbuQamar, S.; Goldsbrough, P. Metallothionein deficiency impacts copper accumulation and redistribution 17

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

440 441 442 443

in leaves and seeds of Arabidopsis. New Phytol. 2014, 202, 940-951. 43. Gasic, K.; Korban, S. S. Expression of Arabidopsis phytochelatin synthase in Indian mustard (Brassica juncea) plants enhances tolerance for Cd and Zn. Planta. 2007, 225, 1277-1285. 44. Sheng, J.P.; Liu, K.L.; Fan, B.; Yuan, Y.; Shen, L.; Ru, B.G. Improving zinc content and

444

antioxidant activity in transgenic tomato plants with expression of mouse metallothionein-I

445

by mt-I gene. J Agr Food Chem. 2007, 55, 9846-9849.

446

45. Ren, Y. J.; Liu, Y.; Chen, H. Y.; Li, G.; Zhang, X. L.; Zhao, J. Type 4 metallothionein genes

447

are involved in regulating Zn ion accumulation in late embryo and in controlling early

448

seedling growth in Arabidopsis. Plant Cell Environ. 2012, 35, 770-789.

449

Page 18 of 26

46. Rodriguez-Llorente, I. D.; Perez-Palacios, P.; Doukkali, B.; Caviedes, M. A.; Pajuelo, E.

450

Expression of the seed-specific metallothionein mt4a in plant vegetative tissues increases Cu

451

and Zn tolerance. Plant Sci. 2010, 178, 327-332.

452 453 454

47. Vatamaniuk, O. K.; Bucher, E. A.; Ward, J. T.; Rea, P. A. Worms take the 'phyto' out of "phytochelatins". Trends Biotechnol. 2002, 20, 61-64. 48. Takahashi, M.; Terada, Y.; Nakai, I.; Nakanishi, H.; Yoshimura, E.; Mori, S.; Nishizawa, N. K.

455

Role of nicotianamine in the intracellular delivery of metals and plant reproductive

456

development. The Plant cell. 2003, 15, 1263-1280.

457

49. Cornu, J. Y.; Deinlein, U.; Horeth, S.; Braun, M.; Schmidt, H.; Weber, M.; Persson, D. P.;

458

Husted, S.; Schjoerring, J. K.; Clemens, S. Contrasting effects of nicotianamine synthase

459

knockdown on zinc and nickel tolerance and accumulation in the zinc/cadmium

460

hyperaccumulator Arabidopsis halleri. New Phytol. 2015, 206, 738-750.

461

50. Deinlein, U.; Weber, M.; Schmidt, H.; Rensch, S.; Trampczynska, A.; Hansen, T. H.; Husted, 18

ACS Paragon Plus Environment

Page 19 of 26

Journal of Agricultural and Food Chemistry

462

S.; Schjoerring, J. K.; Talke, I. N.; Kramer, U.; Clemens, S. Elevated nicotianamine levels in

463

Arabidopsis halleri roots play a key role in Zinc hyperaccumulation. The Plant cell. 2012, 24,

464

708-723.

19

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

465

Page 20 of 26

Table 1 Growth and Zn accumulation in LLZ and WT under 0.5, 5 and 50 µM Zn conditions. 5µM

0.5µM LLZ Dry weight (g/plant)

Zn concentration (µg/g)

Whole plant Shoot Root Whole plant Shoot Root S/R

WT

LLZ

5µM WT

LLZ

WT

3.08±0.63**

4.23±0.26

2.89±0.87*

4.23±0.26

2.94±0.36**

4.09±0.39

2.50±0.57** 0.58±0.10

3.72±0.24 0.51±0.04

2.40±0.61** 0.49±0.12

3.72±0.24 0.51±0.04

2.44±0.32** 0.50±0.06

3.54±0.29 0.55±0.08

58.38±2.59

61.15±3.78

84.81±2.59

61.15±3.78

240.86±18.33

248.64±14.05

61.28±1.74 54.37±2.07** 1.13

59.22±2.94 69.81±1.06 0.85

85.77±2.21** 77.3±1.08** 1.11

59.22±2.94 69.81±1.06 0.85

232.86±9.20 320.95±14.40** 0.73

229.06±6.37 442.96±7.63 0.52

466

* and ** indicated there were significant differences between LLZ and WT under the same Zn condition at 0.05 and 0.01 level respectively. The values given are means±SE

467

(n=6).

468 469

20

ACS Paragon Plus Environment

Page 21 of 26

470

Journal of Agricultural and Food Chemistry

Table 2 Significant and consistent different metabolites between LLZ and WT under 0.5 and 50 µM Zn condition.

471

0.5 µM

50 µM

Name WT

LLZ

VIP

P-value

Q-value

Log2(LLZ/WT)

WT

LLZ

VIP

P-value

Q-value

Log2(LLZ/WT)

Amides

malonamide

0.004373

0.006878

1.452720

0.007570

0.037835

0.653303

0.003789

0.006815

1.515390

0.000005

0.000042

0.846681

Amines

putrescine

0.318642

0.195162

1.426440

0.009328

0.044069

-0.707266

0.390817

0.230128

1.413020

0.000291

0.000766

-0.764060

Amino acids

asparagine

0.495972

0.893163

1.638140

0.000448

0.004588

0.848665

0.664333

1.612813

1.515320

0.000005

0.000042

1.279601

citrulline

0.016312

0.040486

1.715440

0.000043

0.000728

1.311509

0.019727

0.053422

1.405370

0.000353

0.000884

1.437237

glutamic acid

0.059774

0.084402

1.652220

0.000322

0.003676

0.497775

0.075415

0.139960

1.431950

0.000189

0.000622

0.892101

glutamine

0.335262

0.566383

1.570330

0.001693

0.013640

0.756490

0.167295

0.262897

1.118920

0.020567

0.021236

0.652109

glycine

0.208921

0.424430

1.717470

0.000043

0.000728

1.022565

0.310111

0.706367

1.522220

0.000004

0.000042

1.187634

lysine

0.068137

0.141847

1.633770

0.000469

0.004588

1.057831

0.079481

0.135788

1.124250

0.020116

0.021236

0.772665

N-acetyl-L-leucine

0.149754

0.057122

1.744570

0.000016

0.000545

-1.390487

0.124004

0.068730

1.293260

0.003427

0.006127

-0.851382

ornithine

0.017881

0.058502

1.777600

0.000001

0.000102

1.710030

0.028683

0.063331

1.535400

0.000001

0.000030

1.142738

tryptophan

0.018536

0.048955

1.606420

0.000876

0.007998

1.401125

0.009145

0.028123

1.448430

0.000113

0.000425

1.620631

0.000411

0.002321

1.784880

0.000001

0.000102

2.496635

0.000596

0.002957

1.516660

0.000005

0.000042

2.309910

allose

0.006561

0.002036

1.726670

0.000033

0.000728

-1.688114

0.006672

0.004416

1.214080

0.008488

0.010901

-0.595473

gentiobiose

0.001811

0.000471

1.752060

0.000016

0.000545

-1.943553

0.001524

0.000478

1.490380

0.000036

0.000189

-1.672321

glucose-6-phosphate

0.438056

0.284239

1.516020

0.003669

0.022964

-0.624011

0.288442

0.210643

1.264330

0.004901

0.007821

-0.453482

sorbitol

0.014946

0.028098

1.515660

0.003669

0.022964

0.910694

0.015113

0.029493

1.335300

0.001693

0.003428

0.964588

Imines

nicotianamine

0.001263

0.003169

1.677260

0.000176

0.002195

1.327817

0.001448

0.005399

1.495980

0.000016

0.000100

1.898905

Lipids

1-monopalmitin

0.000743

0.001280

1.545200

0.003688

0.022964

0.784933

0.001040

0.001513

1.207910

0.012399

0.015184

0.541373

Nucleosides

flavin adenine dinucleotide

0.153576

0.332123

1.558430

0.002017

0.015351

1.112766

0.131737

0.220004

1.026060

0.040507

0.034968

0.739863

Organic acids

2,4-diaminobutyric acid

0.014540

0.039276

1.727370

0.000033

0.000728

1.433589

0.023460

0.046971

1.378910

0.000715

0.001711

1.001562

3-hydroxybutyric acid

0.001747

0.000065

1.497250

0.004493

0.025647

-4.747391

0.001905

0.001493

1.129870

0.011533

0.014460

-0.351935

glutaconic acid

0.000242

0.000465

1.689920

0.005303

0.029062

0.942371

0.000285

0.000898

1.537720

0.000012

0.000087

1.657844

pimelic acid

0.006566

0.012152

1.677470

0.000176

0.002195

0.888071

0.008016

0.021560

1.559120

0.000000

0.000000

1.427481

succinic acid

0.336809

0.217613

1.673370

0.000176

0.002195

-0.630168

0.250803

0.142045

1.418110

0.000260

0.000719

-0.820208

Carbohydrates 2-deoxy-D-glucose

21

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 22 of 26

472

Figure captions

473

Fig 1. The Zn accumulation of different plant tissues under 0.5, 5 and 50 µM Zn conditions. The results are presented as means±SE. The values

474

followed by asterisks are statistically different from the wild type according to a Student’s t-test (n = 6; *P < 0.05; **P < 0.01).

475

Fig 2. Principal component analysis (PCA) of metabolic profiles in grains of LLZ and WT under 0.5 and 50 µM Zn conditions (six biological

476

replicates). Circle and triangle respectively denote metabolomes of LLZ and WT under 0.5µM Zn condition, while inverted triangle and box correspond

477

to LLZ and WT under 50 µM Zn conditions.

478

Fig 3. Metabolic changes of 15 DFA grains in LLZ and WT after 50 µM Zn treatment and the pathway they participated. Box and circle respectively

479

denote significantly changed metabolites of WT and LLZ (VIP>1, P