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