Variation in mineral elements in grains of 20 brown rice accessions in two environments

Variation in mineral elements in grains of 20 brown rice accessions in two environments

Accepted Manuscript Variation in Mineral Elements in Grains of 20 Brown Rice Accessions in Two Environments Yan Huang, Chuan Tong, Feifei Xu, Yaling C...

445KB Sizes 2 Downloads 29 Views

Accepted Manuscript Variation in Mineral Elements in Grains of 20 Brown Rice Accessions in Two Environments Yan Huang, Chuan Tong, Feifei Xu, Yaling Chen, Caiya Zhang, Jinsong Bao PII: DOI: Reference:

S0308-8146(15)01120-6 http://dx.doi.org/10.1016/j.foodchem.2015.07.087 FOCH 17883

To appear in:

Food Chemistry

Received Date: Revised Date: Accepted Date:

29 April 2015 18 July 2015 20 July 2015

Please cite this article as: Huang, Y., Tong, C., Xu, F., Chen, Y., Zhang, C., Bao, J., Variation in Mineral Elements in Grains of 20 Brown Rice Accessions in Two Environments, Food Chemistry (2015), doi: http://dx.doi.org/ 10.1016/j.foodchem.2015.07.087

This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

1 2

Variation in Mineral Elements in Grains of 20 Brown Rice Accessions in Two Environments

3 Yan Huang1, Chuan Tong1, Feifei Xu1, Yaling Chen1, Caiya Zhang2*, Jinsong

4

Bao1*

5 6 7 8

1

Institute of Nuclear Agricultural Sciences, College of Agriculture and

Biotechnology, Zhejiang University, Hangzhou, 310029, China.

9

2

10

China

Department of Statistics, Zhejiang University City College, Hangzhou 310015,

11 12

*Corresponding author: phone +86-571-88284309, and email [email protected]

13

(C. Zhang); and phone +86-571-86971932; fax +86-571-86971421; and email:

14

[email protected] (J. Bao).

15 16

Running Title: Variation in mineral elements in grains of 20 brown rice

17 18

1

19

Abstract

20

Twenty rice accessions were planted in Hainan province, China, for 2 years to

21

investigate the effects of genotype, environment, and their interactions on the Ca, Mg,

22

Na, K, Fe, Zn, and Cu contents in brown rice. Analysis of variance showed that the Ca,

23

Na and K were mainly affected by the genotypic variance, whereas the Fe, Zn and Cu

24

were mainly affected by the environment variance. The genotype × environment

25

interaction effects for Mg, Na, Zn, and Cu were highly significant (P<0.001), though

26

it only accounted for a small proportion of the total variation (0.5-16.3%). The

27

correlation analyses showed that Mg was significantly positively correlated with K,

28

Fe, and Zn. A total of 9 and 8 single nucleotide polymorphism (SNP) loci were

29

identified in 2011 and 2012, respectively, which were strongly associated with for Ca,

30

Cu, K, Na, and Zn.

31 32

Keywords: Rice; mineral content; genotype and environment interaction, single

33

nucleotide polymorphism

34

2

35

1. Introduction

36

Rice (Oryza sativa) is the most important staple food crop in the world and

37

provides starch, protein, and other essential nutrients for over half of the global

38

population. At least 49 nutrients are required by humans for their normal growth and

39

development (Welch & Graham, 2004), and the demand for most nutrients is generally

40

supplied by cereals, particularly rice as a staple food. Along with an ever growing

41

number of global populations, especially in regions like Asia where rice is regarded as

42

the basic food source, nutritional quality of the grain plays an indispensable role upon

43

human health.

44

Among these nutrients, mineral elements play numerous beneficial roles due to

45

their direct or indirect effects in both plant and human metabolism, otherwise their

46

deficiencies or insufficient intakes may lead to several dysfunctions and diseases in

47

humans (Garcia-Oliveira, Tan, Fu, & Sun, 2009). It was suggested that relatively low

48

intake of Ca in most developing countries was accompanied by the increased risk of

49

osteoporosis (Nordin, 2000; Welch & Graham, 1999). The deficiency of Fe and Zn, two

50

of the major micronutrients, are affecting more than 2 billion people worldwide

51

(Sandstead, 1991); Fe deficiency results in suffering from anemia, while zinc

52

deficiency causes stunted growth and underdeveloped intelligence (Umeta, West,

53

Haidar, Deurenberg, & Hautvast, 2000). Although rice is not considered to be

54

mineral-rich, it still can be an important source for those who eat rice as staple food

55

since it provides human with caloric energy and minerals at the same time (Zhang,

56

Pinson, Tarpley, Huang, Lahner, & Yakubova, 2014). Therefore, even a small increase

57

in the mineral elements concentrations in rice grain, in regions especially in developing

58

countries would be highly significant. Biofortification of mineral nutrients in rice is

59

extremely necessary for rice breeders, which is considered an efficient way for

60

ameliorating mineral deficiency in humans.

61 62

Mineral elements accumulation in rice grains is a complex trait. Quantitative trait locus (QTL) mapping and association analysis are considered two powerful tools for 3

63

discovering the genes/QTLs underlying such complex traits.

In

recent

years,

64

tremendous efforts have been made on the genetic analysis of different mineral

65

elements in rice. Stangoulis, Huynh, Welch, Choi, and Graham (2007) mapped the

66

QTLs for inorganic phosphorus (P), total P, Fe, Zn, Cu and Mn concentrations. Norton,

67

Deacon, Xiong, Huang, Meharg, & Price (2010) identified 41 QTLs for the

68

concentration of 17 elements in rice grain. Garcia-Oliveira et al. (2009) identified 31

69

putative QTLs for Fe, Zn, Mn, Cu, Ca, Mg, P, and K contents of 85 introgression lines

70

derived from a cross between an elite indica cultivar Teqing and the wild rice (Oryza

71

rufipogon), and found that the wild rice contributed favorable alleles for most of the

72

QTLs.

73

In addition, many researchers have explored the impact of environment on the

74

accumulation of mineral elements in rice. Du, Zeng, Wang, Qian, Zheng, and Ling

75

(2013) determined the Ca, Fe, K, Mg, Mn, P, and Zn contents in brown rice and

76

identified 23 and 9 QTLs in two different ecological environments, respectively. Only

77

2 QTLs for Mg content have been detected in both environments, indicating that

78

mineral accumulation QTLs are greatly affected by the environment. Pinson, Tarpley,

79

Yan, Yeater, Lahner, & Yakubova (2015) analyzed 16 mineral elements concentrations

80

in brown rice which were produced over 2 years in Beaumont, TX, under both flooded

81

and unfolded watering regimes, and found that environmental variance was large for

82

some elements.

83

However, it is still unclear how genotypes, environment, and genotypes ×

84

environment interaction influence the mineral accumulation in rice grains. The

85

objectives of this study were (1) to investigate the genetic diversity on the

86

accumulation of mineral elements among 20 diverse rice accessions representing a

87

wide geographical origins in two environments (years), and (2) to explore additional

88

genetic loci (QTLs) for the mineral elements which may be useful for molecular

89

breeding for improving the mineral nutrients in rice grains.

90

2. Materials and method 4

91

2.1. Rice materials

92

Twenty rice accessions (Oryza sativa L.) used for this study, including cultivars,

93

germplasm lines, and landraces from 10 geographical areas, was introduced from the

94

OryzaSNP project (http://www.oryzasnp.org/) (McNally, Childs, Bohnert, Davidson,

95

Zhao, &Ulat, 2009) . All the 20 rice accessions were planted in a randomized block

96

design with two replications in two years at the same site in Lingshui, Hainan

97

province, China. Each accession in the fields was planted into ten rows with six plants

98

per row. The seeds were sowed in early December 2010 and 2011, and harvested in

99

April2011 and 2012, respectively. Hereafter two years (environments) are denoted

100

2011 and 2012 for when the grain was harvested. Field management followed

101

conventional practices including nutrient and pest control procedures to ensure the

102

production in the field. The major environmental conditions, i.e. the mean temperature,

103

total radiation hours and total rainfall for each month during rice growth periods were

104

reported previously (Tong, Chen, Tang, Xu, Huang, Chen, et al. 2014).

105 106

The grains were air-dried to a moisture content of about 12%, stored at room

107

temperature for two months and later at 4 ℃ in the dark until all the grain materials

108

were obtained. The samples were dehusked on a Satake Rice Machine (Satake Co.,

109

Tokyo, Japan) and the dehusked rice samples were ground to flours, which were used

110

for further analyses.

111

2.2. Sample Treatment and Determination of Mineral content

112

Sample preparation and determination of mineral contents were according to

113

Jiang, Wu, Feng, Yang, & Shi (2007). Rice powder (0.5 g) was weighed and placed

114

into a crucible, and then carbonized at 250 °C on an electro thermal plate until the

115

sample turned into black but not smoking completely. The crucibles with samples

116

were dry-ashed by heating in a muffle furnace at 550 °C (about 10–12 h). After

117

sample incineration, a white residue was obtained, which was carefully transferred

118

into a 50 mL volumetric flask, dissolved with 5 mL of 6 M HCl and then diluted to 50 5

119

mL with water. The diluted solutions were subjected to analysis for Ca, Mg, Na, K, Fe,

120

Zn, and Cu by inductively coupled plasma mass spectroscopy (ICP-MS) (Agilent

121

7500A; Agilent Technologies, Inc.).

122

2.3. Statistical Analysis

123

All concentrations of the mineral elements were measured at least in duplicate.

124

Data were presented as mean value ± standard deviation. Data analyses were

125

performed with the SAS program version 9.1.3 (SAS Institute Inc., Cary, NC,U.S.A.).

126

Analysis of variance (ANOVA) was carried out to determine genotypic and

127

environmental variation among the parameters using the general linear model

128

procedure (PROC GLM). Means of different genotypes were determined using the

129

PROC MEANS, followed by least significant difference (LSD) multiple comparison

130

tests at P < 0.05. Duncan’s new multiple-range test was performed to examine

131

significant differences among subpopulations at P < 0.05. The t-test was performed to

132

compare the mean values between two years. Correlation analysis among the different

133

mineral elements in (or between) the two environments was calculated by the PROC

134

CORR procedure.

135

2.4. Genotype data and association mapping

136

The publically available genotype data were downloaded from the Gramene Web

137

site (http://www.gramene.org/) on March 10, 2013. During downloading, we removed

138

the reference genome data firstly, and we then filtered sites with the following

139

parameters: minimum count was 16 out of 20 sequences, and minimum frequency

140

was set at 0.25. Finally, a total of 32 655 SNPs distributed in 12 chromosomes were

141

downloaded (Xu, Tang, Shao, Chen, Tong, & Bao, 2014). Association mapping was

142

performed with GAPIT (Lipka, Tian, Wang, Peiffer, Li, & Bradbury, 2012).

143

Comparison of models and selection of the best model was referred to Xu et al (2014).

144

The results with the best model were presented at the significance level of P < 0.005.

145 6

146

3. Results

147

3.1. Variation in mineral elements concentrations among rice accessions

148

The concentrations of mineral elements in brown rice among 20 rice accessions

149

produced over two years (2011 and 2012) are summarized in Tables 1 and 2. Both

150

macro elements (Ca, Mg, Na, K) and micro-elements (Fe, Zn, Cu) varied in

151

concentrations among different rice accessions. In 2011, the variation for Cu among

152

different rice was largest with ratio of the maximum/minimum content of 5.6-fold

153

ranging from the minimum values of 0.217 mg/kg in rice G03 to the maximum value

154

of 1.216 mg/kg in rice G15. The content of Mg showed the smallest variation with

155

1.35-fold (1.52-2.05 g/kg) (Table 1). In 2012, the mineral contents also varied among

156

20 genotypes with ratios of the maximum/minimum content of 2.52-fold (73.5-185.1

157

mg/kg) for Ca, 1.63-fold (1.42-2.31 g/kg) for Mg, 5.07-fold (7.67-38.86 mg/kg) for Na,

158

1.47-fold (2.55-3.74 g/kg) for K, 3.95-fold (4.01-15.81 mg/kg) for Fe, 2.14-fold

159

(17.98- 38.52 mg/kg) for Zn, 2.36-fold (1.20-2.84 mg/kg) for Cu(Table 2). Moreover,

160

the mean value of Mg, Fe, Zn and Cu were of significant difference between two years

161

(Table 2), suggesting influence of environment on these mineral elements.

162

Of the 20 rice genotypes studied, eight represented the japonica subspecies, eight

163

indica subspecies and four belonged to the aus group (Table 1; Table 2) (Xu et al.

164

2014). The concentrations of the 7 elements also varied among different subgroups in

165

both years. The mean values of Na, K, Fe and Zn in japonica group were higher than

166

those in indica and Aus groups. The aus group had higher mean Ca content than the

167

japonica and indica groups, while the highest mean value of Cu was found in the

168

indica group (Table 3).

169

3.2. Correlation analysis

170

The result of pair-wise correlation analysis among all mineral elements of the 20 rice

171

accessions is listed in Table 4.The correlation coefficients for Ca, Mg, Na, K and Fe

172

between 2 years ranged from 0.56-0.81 (P< 0.001), suggesting that these elements 7

173

were relatively stable in different planting years. The content of Mg was positively

174

correlated with those of K, Fe, and Zn. The content of Fe was positively correlated

175

with the content of Na and K. All these correlations were observed in two years. The

176

content of Cu was negatively correlated with Fe and positively correlated with Zn,

177

which was only observed in 2011.

178

3.3. Analysis of variance for mineral elements

179

Analysis of variance (ANOVA) for mineral elements indicated that the genotype

180

variance for 7 mineral elements were highly significant (P<0.001), but the environment

181

variance for some elements such as Mg, Fe, Zn and Cu were also significant (P<0.001)

182

(Table 5). Additionally, the genotype × environment interaction for most elements

183

except for Ca, K and Fe were significant (P< 0.001), though it only explained a small

184

proportion of the total variation (0.5-16.3%). The genotype effects accounted for more

185

than 70% of the total variation for the elements of Ca, Na and K, indicating that the

186

genotype was the primary factor determining the difference of the mineral elements. In

187

contrast, the environment effects accounted for more than 80% of total variation for the

188

elements of Fe, Zn and Cu (Table 5), indicating that these elements were highly

189

effected by the environment.

190 191

3.4. Association test for mineral elements

192

Association mapping was used to find the genetic factor affecting the genetic

193

diversity in the 7 mineral elements. Seventeen QTLs locating on chromosome 1, 2, 6, 8,

194

10 and 11 were identified for 5 of 7 mineral elements except for Mg and Fe using

195

preliminary association test (Table 6; Supplementary Fig. 1). The Manhattan plots of

196

association mapping results for K and Na were shown in Supplementary Fig. 2 in the

197

two years, respectively.

198

Nine of 17 QTLs were detected in 2011, while the others were identified in the

199

2012. Only one QTL for Cu and Zn content was detected in 2011, indicating that the

200

QTLs for Cu and Zn in rice grains were largely environment-dependent. Five QTLs 8

201

were identified for Ca in the two years, of which two QTLs on chromosome 1 and 6

202

with different positions could be detected in each year. A large number of QTLs for K

203

and Na were detected (Table 6; Supplementary Fig.1 and Fig.2), but none of them

204

could be detected simultaneously in both years.

205

4. Discussion

206

The 20 rice accessions introduced from the OryzaSNP project represent wide

207

genetic diversity in various phenotypes (McNally et al, 2009), such as agronomic traits

208

(Xu et al, 2014), biomass traits (Jahn, McKay, Mauleon, Stephens, McNally, Bush, et

209

al., 2011), bioactive compounds (Heuberger, Lewis, Chen, Brick, Leach, & Ryan, 2010;

210

Shao, Tang, Huang, Xu, Chen, Tong, et al., 2014), and grain quality traits (Tong et al,

211

2014). However, genetic and environmental effects on the mineral elements in 20 rice

212

accessions have not been reported before.

213

In this study, the concentrations of Ca, Mg, Na, K, Fe, Zn, and Cu in brown rice

214

of 20 genotypes were evaluated. Each of these nutrient are essential for the growth

215

and development of plants, animals, and humans. The results showed that 7 mineral of

216

the 20 rice accessions varied greatly among genotypes. For example, the largest

217

variation of mineral elements concentrations were 5.6-fold (0.217-1.216 mg/kg) for Cu

218

in 2011 and 5.07-fold (7.67-38.86 mg/kg) for Na. Genetic variation for mineral

219

contents has been reported by many research groups in different rice accessions

220

(Gregorio, 2002; Jiang et al, 2007; Garcia-Oliveira et al, 2009; Anuradha, Agarwal,

221

Rao, Rao, Viraktamath, & Sarla, 2012). Brar, Jain, Singh, and Jain (2011) showed

222

large variation in Fe (5.1–441.5 µg/g) and Zn (2.12–39.4 µg/g) contents in brown rice

223

in a collection of 220 rice genotypes and identified three rice genotypes with

224

exceptionally higher Fe content (~400µg/g). Jiang et al (2007) indicated that the

225

differences in mineral contents among genotypes might be due to different genetic

226

resources involved. In this study, the genotype effects accounted for more than 70% of

227

variance for the elements of Ca, Na and K, which implied that genotypic variations 9

228

might provide opportunities to select for higher mineral element contents (Gregorio,

229

2002).

230

The mineral contents among different rice accessions were not only influenced

231

by genotypes but also affected by the environment factors. Some previous studies

232

have demonstrated that environment plays a tremendous impact on the accumulation

233

of minerals in rice grain. Du et al (2013) found the mineral concentrations in grains of

234

the two parents grown in Hangzhou and Lingshui, China, were significantly different,

235

especially for Ca, Fe, and Zn. Variations in mineral concentrations in plants may

236

depend on variations in many factors, such as mineral mobilization, uptake,

237

trafficking, and sequestration, which are all relevant processes in the mineral transport

238

pathway from roots to shoots (Ghandilyan, Ilk, Hanhart, Mbengue, Barboza, Schat, et

239

al. 2009; Clemens, 2001). Moreover, the climate may also affect the elemental uptake

240

by modifying plant growth and metabolism (Du et al, 2013). It is straightforward to

241

ascribe the difference in the mean value of most mineral elements except Ca and Na

242

among 20 rice accessions between two years to the effects of environments (Table 2).

243

Under the same environment, different rice accessions also showed different

244

concentrations in 7 mineral elements, which may be ascribed to the genotypic

245

differences (Table 1 and 2). The mean values of most mineral elements among

246

different subgroups showed significant difference in two years, further supporting the

247

influence of genotypes on the mineral accumulation (Table 3). Moreover, ANOVA

248

indicated that the genotype variance for 7 mineral elements were highly significant

249

(P<0.001), while the environment variance for some elements such as Mg, Fe, Zn and

250

Cu were also significant (P<0.001) (Table 4), also suggesting that the mineral elements

251

were affected by hereditary as well as environment. Batten (2002) indicated that the

252

correlation between mineral elements concentrations might be affected by the

253

genotype ×environment, which was supported by our study. The Mg, Na, Zn, and

254

Cu were greatly affected by genotype × environment interaction, though it only

255

explained a small proportion of the total variation. So to improve the mineral nutrients, 10

256

breeders should not ignore the effect of genotype × environment interaction in their

257

breeding lines in different environments.

258

The correlation analyses between 7 mineral elements showed that Mg had a

259

significantly positive correlation with K, Fe, and Zn. Jiang et al, (2007) showed that

260

Mg was negatively correlated with Cu and positively correlated with K, Ca, Na, Fe,

261

Zn and Mn. Ca and was positively correlated Na, Mg, Fe, Zn and Mn. This might be

262

due to the interaction between ions in which chemical properties were equivalently

263

similar, which may compete for the site of absorption, transport, and function in

264

plant tissues (Gussarsson, Adalsteinsson, Jensen, & Asp, 1995; Robson et al, 1983).

265

These correlation results from present and previous studies could help improving the

266

efficiency of early generation selection for rice materials abundant in Mg, K, Fe, and

267

Zn in rice breeding program.

268

Many putative genomic regions associated with the phenotypes using this set of

269

materials have been identified. For example, Xu et al. (2014) identified 23 QTLs for 10

270

agronomic traits in the two years. Tong et al. (2014) identified 22 QTLs for starch

271

physicochemical properties in which some starch biosynthesis related genes such as

272

Wx, SSIIIb, ISA1, and ISA2 were revealed. Shao et al. (2014) found many potential loci

273

for polyphenol and antioxidant capacity traits. Therefore, it is possible to detect the

274

putative genomic regions or QTLs associated with different mineral elements using this

275

set of rice materials. In this study, a total of 17 QTLs for the mineral elements except

276

Mg and Fe were discovered which were located on chromosome 1, 2, 6, 8, 10 and 11.

277

We identified 5 loci for Ca content on chromosomes 2, 6, 8 which were not reported

278

before (Du et al, 2013;Garcia-Oliveira et al, 2009). Only one QTL (qCu6) was

279

identified on the chromosome 6, and it was not the same region as reported by

280

Garcia-Oliveira et al, (2009). Consistent with previous research, qK6.1 was close to

281

the QTLs identified by Du et al, (2013) on chromosome 6, suggesting that the locus

282

may be a hot region for K accumulation in rice grain and more further research on this

283

region is needed. Kumar, Jain, & Jain, (2014) detected one QTL for zinc content on 11

284

chromosome 2 which is close to the position of the qZn2 in the study. Few QTLs were

285

identified for Na content in rice grain, our study identified 5 new QTLs located at

286

different regions from the previous studies, and so further studies are needed to

287

confirm the results.

288

5. Conclusion

289

Wide genetic diversity in the 7 mineral elements in 20 brown rice accessions

290

was discovered in two years. These mineral elements were affected by genotype as

291

well as environment. The effects of genotype × environment interaction should also

292

be considered in different environments. Moreover, we discovered 17 QTLs involved

293

in the accumulation of 7 mineral elements. This study provided new information on

294

the understanding of the contributions of the genotype and environment interaction

295

effects to the complex traits such as mineral contents in the brown rice. The results

296

will also be useful for developing rice varieties to improve the mineral nutrients in

297

rice grain by molecular breeding.

298

Acknowledgement

299

This work was financially supported by the Special Fund for Agro-scientific

300

Research in the Public Interest (201103007) from the Ministry of Agriculture. We

301

thank Dr. Ruaraidh Sackville Hamilton in the International Rice Research Institute for

302

providing the rice materials used in this study.

303

References

304

Anuradha, K., Agarwal, S., Rao, Y. V., Rao, K. V., Viraktamath, B. C., & Sarla, N.

305

(2012). Mapping QTLs and candidate genes for iron and zinc concentrations in

306

unpolished rice of Madhukar  Swarna RILs. Gene, 508, 233-240.

307

Batten, G. D. (2002). Relating minerals in rice shoots and grain to soil tests, yield and

308

grain quality. Rural Industries Research and Development Corporation. New

309

South Wales of Australia, Pp1-32. 12

310

Brar, B., Jain, S., Singh, R., & Jain, R. K. (2011). Genetic diversity for iron and zinc

311

contents in a collection of 220 rice (Oryza sativa L.) genotypes. Indian Journal of

312

Genetics and Plant Breeding, 71, 67-73.

313 314

Clemens, S. (2001). Molecular mechanisms of plant metal tolerance and homeostasis. Planta, 212, 475-486.

315

Du, J., Zeng, D., Wang, B., Qian, Q., Zheng, S., & Ling, H. Q. (201 3). Environmental

316

effects on mineral accumulation in rice grains and identification of ecological

317

specific QTLs. Environmental Geochemistry and Health, 35, 161-170.

318

Garcia-Oliveira, A. L., Tan, L., Fu, Y., & Sun, C. (2009). Genetic identification of

319

auantitative trait loci for contents of mineral nutrients in rice grain. Journal of

320

Integrative Plant Biology, 51, 84-92.

321

Ghandilyan, A., Ilk, N., Hanhart, C., Mbengue, M., Barboza, L., Schat, H., Koornneef,

322

M., El-Lithy, M., Vreugdenhil, D., Reymond, M., & Aarts, M. G. M. (2009). A

323

strong effect of growth medium and organ type on the identification of QTLs for

324

phytate and mineral concentrations in three Arabidopsis thaliana RIL populations.

325

Journal of Experimental Botany, 60, 1409-1425.

326 327

Gregorio, G. B. (2002). Progress in breeding for trace minerals in staple crops. Journal of Nutrition, 132, 500S-502S.

328

Gussarsson, M., Adalsteinsson, S., Jensen, P., & Asp, H. (1995). Cadmium and coppers

329

interactions on the accumulation and distribution of Cd and Cu in birch

330

(Betulapendula Roth) seedlings. Plant and Soil, 171, 185-187.

331

Heuberger, A. L., Lewis, M. R., Chen, M.H., Brick, M. A., Leach, J. E., & Ryan, E. P.

332

(2010). Metabolomic and functional genomic analyses reveal varietal differences

333

in bioactive compounds of cooked Rice. Plos One, 5, e12915.

334

Jahn, C. E., McKay, J. K., Mauleon, R., Stephens, J., McNally, K. L., Bush, D. R.,

335

Leung, H., & Leach, J. E. (2011). Genetic variation in biomass traits among 20

336

diverse rice varieties. Plant Physiology, 155, 157-168.

13

337

Jiang, S. L., Wu, J. G., Feng, Y., Yang, X. E., & Shi, C. H. (2007). Correlation analysis

338

of mineral element contents and quality traits in milled rice (Oryza stavia L.).

339

Journal of Agricultural and Food Chemistry, 55, 9608-9613.

340

Kumar, J., Jain, S., & Jain, R. K. (2014). Linkage Mapping for Grain Iron and Zinc

341

Content in F-2 Population Derived from the Cross between PAU201 and Palman

342

579 in Rice (Oryza sativa L.). Cereal Research Communications, 42, 389-400.

343

Lipka, A. E., Tian, F., Wang, Q., Peiffer, J., Li, M., Bradbury, P. J., Gore, M. A., Buckler,

344

E. S., & Zhang, Z. (2012). GAPIT: genome association and prediction integrated

345

tool. Bioinformatics, 28, 2397-2399.

346

Lu, K. Y., Li, L. Z., Zheng, X. F., Zhang, Z. H., Mou, T. M., & Hu, Z. L. (2008).

347

Quantitative trait loci controlling Cu, Ca, Zn, Mn and Fe content in rice grains.

348

Journal of Genetics, 87, 305–310.

349

McNally, K. L., Childs, K. L., Bohnert, R., Davidson, R. M., Zhao, K., Ulat, V. J., Zeller,

350

G., Clark, R. M., Hoen, D. R., Bureau, T. E., Stokowski, R., Ballinger, D. G.,

351

Frazer, K. A., Cox, D. R., Padhukasahasram, B., Bustamante, C. D., Weigel, D.,

352

Mackill, D. J., Bruskiewich, R. M., Raetsch, G.,

353

Leach, J. E. (2009). Genomewide SNP variation reveals relationships among

354

landraces and modern varieties of rice. Proceedings of the National Academy of

355

Sciences of the United States of America, 106, 12273-12278.

356 357

Buell, C. R., Leung, H., &

Nordin, B. E. C. (2000). Calcium requirement is a sliding scale. American Journal of Clinical Nutrition, 71, 1381-1383.

358

Norton, G. J., Deacon, C. M., Xiong, L., Huang, S., Meharg, A. A., & Price, A. H.

359

(2010). Genetic mapping of the rice ionome in leaves and grain: identification of

360

QTLs for 17 elements including arsenic, cadmium, iron and selenium. Plant and

361

Soil, 329, 139-153.

362

Pinson, S. R. M., Tarpley, L., Yan, W., Yeater, K., Lahner, B., Yakubova, E., Huang,

363

X.Y., Zhang, M., Guerinot, M. L., & Salt, D. E. (2015). Worldwide genetic

14

364

diversity for mineral element concentrations in rice grain. Crop Science, 55,

365

294-311.

366

Robson, A. D., & Pitman, M. G. (1983). Interactions between nutrients in higher

367

plants. In: Lauchli, A. & Bieleski, R. L. (eds), Inorganic Plant Nutrition;

368

Encyclopedia of Plant Physiology, New Series, vol 15A , Pp147-180,

369

Springer-Verlag, Berlin.

370 371

Sandstead, H. H. (1991). Zinc deficiency- A public health problem. American Journal of Diseases of Children, 145, 853-859.

372

Shao, Y., Tang, F., Huang, Y., Xu, F., Chen, Y., Tong, C., Chen, H., & Bao, J. S. (2014).

373

Analysis of genotype  environment interactions for polyphenols and antioxidant

374

capacity of rice by association mapping. Journal of Agricultural and Food

375

Chemistry, 62, 5361-5368.

376

Stangoulis, J. C. R., Huynh, B.L., Welch, R. M., Choi, E.Y., & Graham, R. D. (2007).

377

Quantitative trait loci for phytate in rice grain and their relationship with grain

378

micronutrient content. Euphytica, 154, 289-294.

379

Tong, C., Chen, Y., Tang, F., Xu, F., Huang, Y., Chen, H., & Bao, J. S. (2014). Genetic

380

diversity of amylose content and RVA pasting parameters in 20 rice accessions

381

grown in Hainan, China. Food Chemistry, 161, 239-245.

382

Umeta, M., West, C. E., Haidar, J., Deurenberg, P., & Hautvast, J. (2000). Zinc

383

supplementation and stunted infants in Ethiopia: a randomised controlled trial.

384

Lancet, 355, 2021-2026.

385

Welch, R. M., & Graham, R. D. (1999). A new paradigm for world agriculture: meeting

386

human needs - Productive, sustainable, nutritious. Field Crops Research, 60, 1-10.

387

Welch, R. M., & Graham, R. D. (2004). Breeding for micronutrients in staple food

388

crops from a human nutrition perspective. Journal of Experimental Botany, 55,

389

353-364.

15

390

Xu, F. F., Tang, F. F., Shao, Y. F., Chen, Y. L., Tong, C., & Bao, J. S. (2014). Genotype ×

391

environment interactions for agronomic traits of rice revealed by association

392

mapping. Rice Science, 21, 133-141.

393

Zhang, M., Pinson, S. R. M., Tarpley, L., Huang, X.Y., Lahner, B., Yakubova, E., Baxter,

394

I., Guerinot, M. L., & Salt, D. E. (2014). Mapping and validation of quantitative

395

trait loci associated with concentrations of 16 elements in unmilled rice grain.

396

Theoretical and Applied Genetics, 127, 137-165.

397

16

398 399

Table 1 The mineral contents of 20 rice accessions grown in Hainan in 2011. Rice

Code

Subgroups Ca (mg/kg)

Mg (g/kg)

Na (mg/kg)

K (g/kg)

Fe (mg/kg)

Zn (mg/kg) Cu (mg/kg)

Accessions G01

Azucena

Japonica

111.5±8.6 1.567±0.083 27.68±1.03 3.025±0.151 11.38±1.55

15.68±0.63 0.436±0.046

G02

Dom Sufid

Japonica

113.5±8.0 2.052±0.138 23.37±2.16 3.565±0.196 16.09±2.47

25.77±0.99 0.646±0.016

G03

Dular

Aus

117.1±7.6 1.696±0.092

18.46±0.01 0.217±0.026

G04

FR13A

Aus

115.6±0.9 1.599±0.104 16.29±2.82 2.968±0.232

9.60±2.06

18.97±0.92 0.442±0.034

G05

IR64-21

Indica

71.4±9.0

1.539±0.089

7.84±1.89

13.76±0.35 0.489±0.114

G06

LTH

Japonica

92.2±1.5

1.643±0.085 26.98±1.23 3.308±0.158 12.47±1.34

15.05±0.39 0.513±0.027

G07

M202

Japonica

108.7±5.2 1.636±0.012 29.73±4.94 3.695±0.181 11.74±1.45

11.5±0.14 0.417±0.005

G08

Minghui63

Indica

80.0±7.8

16.2±4.72

7.96±1.74

3.355±0.169 11.06±4.60

2.638±0.173

1.633±0.067

10.51±1.1

3.702±0.201

8.24±1.38

19.22±0.5 0.944±0.082

G09 Moroberekan Japonica

137.9±11.7 1.521±0.112

20.49±0.2

2.871±0.23

11.48±2.38

14.82±0.78 0.302±0.021

G10

N22

Aus

135.0±5.7

1.834±0.09

16.33±1.37 3.195±0.222 13.35±0.16

17.3±0.76 0.221±0.008

G11

Nipponbare

Japonica

77.5±7.2

1.854±0.097 30.98±3.27 3.596±0.245 13.88±2.62

14.81±0.27 0.519±0.033

G12

Pokkali

Indica

111.1±5.3 1.835±0.131 14.67±1.04 3.569±0.319

6.53±2.17

G13

Sadu-Cho

Indica

173.5±16.9 1.889±0.157

21.13±3.5

3.502±0.338

9.58±1.89

G14

SHZ-2

Indica

136.2±1.5

1.573±0.05

8.71±0.7

2.398±0.083

6.96±0.98

14.79±0.34 0.687±0.031

G15

Swarna

Indica

92.1±0.6

1.734±0.082 13.35±0.37 2.707±0.115

7.86±0.32

21.32±0.29 1.216±0.034

G16

Tainung67

Japonica

106.5±12.9 1.578±0.037 31.92±3.58 3.289±0.109

9.55±0.09

19.77±1.08 0.791±0.073

G17 Zhenshan97B

Indica

124.3±5.1 1.788±0.028 13.19±0.22 3.423±0.044

8.58±0.16

14.31±0.94 0.546±0.051

G18

Aswina

Indica

163.7±59.4 1.836±0.126

G19

Cypress

Japonica

G20

Rayada

Aus Average 1

LSD

400 401 402 403 404 405 406 407 408 409 410 411

92.6±1.4

18.71±1.53 0.579±0.059 21.81±2

0.739±0.055

30.3±0.21

3.687±0.353

9.50±1.20

21.29±1.24 0.787±0.077

1.792±0.037 15.76±2.12

3.36±0.097

12.50±0.81

18.43±0.1 0.484±0.023

26.5±1.03

3.06±0.042

10.18±0.88

19.49±2.35 0.451±0.013

160.3±0.9 1.754±0.026 116.0

1.718

20.10

3.246

10.42

17.76

0.571

31.7

0.190

4.92

0.422

3.85

2.07

0.104

Data were presented as mean value ± standard deviation. 1 LSD: least significant difference (P < 0.05).

412 17

413

Table 2

414

The mineral contents of 20 rice accessions grown in Hainan in 2012. Rice

Code

Subgroups Ca (mg/kg)

Mg (g/kg)

Na (mg/kg)

K (g/kg)

Fe (mg/kg) Zn (mg/kg)

Cu (mg/kg)

Accessions G01

Azucena

Japonica

100.6±2.0

1.566±0.023 44.53±2.72 3.121±0.026 9.17±0.61

G02

Dom Sufid

Japonica

136.2±5.0

2.306±0.087

G03

Dular

Aus

90.9±8.0

1.441±0.083 11.39±3.65 2.953±0.153 9.13±0.06

G04

FR13A

Aus

G05

IR64-21

G06

31.49±1.45 2.837±0.153

27.37±7.4 3.736±0.052 15.81±0.84 38.52±0.12 2.398±0.284 25.09±0.98 1.285±0.072

168.7±87.1 1.573±0.004

17.86±1.9 2.917±0.025 7.65±0.97

Indica

73.5±10.3

1.511±0.145

7.87±0.73 2.654±0.233 6.17±1.25

LTH

Japonica

97.8±1.8

1.69±0.028

22.62±3.6 3.385±0.081 8.53±0.69

24.21±0.47 1.772±0.152

G07

M202

Japonica

111.7±13.5

1.53±0.016

19.29±0.28 3.75±0.109

21.67±0.25

G08

Minghui63

Indica

77.1±1.9

1.416±0.014

G09 Moroberekan Japonica

106.6±0.9

1.483±0.002 16.49±1.48 3.046±0.025 9.16±0.54

25.83±1.32

G10

N22

Aus

126.8±4.5

1.86±0.002

18.19±0.64 3.22±0.022

28.42±0.77 1.623±0.029

G11

Nipponbare

Japonica

84.0±3.3

1.889±0.015

18.58±1.5 3.785±0.085 11.37±0.51

26.01±1.4

G12

Pokkali

Indica

120.0±27.6

1.676±0.02

16.96±1.07 3.463±0.062 4.01±0.10

20.93±0.48 1.779±0.064

G13

Sadu-Cho

Indica

122.9±4.2

1.843±0.041 14.65±2.19 3.532±0.039 7.42±0.34

23.6±1.21

2.097±0.117

G14

SHZ-2

Indica

140.7±7.5

1.563±0.062 11.41±0.65 2.549±0.133 5.09±0.43

17.98±1.49

2.757±0.15

G15

Swarna

Indica

81.8±2.6

1.628±0.065

9.93±2.26

4.88±0.35

19.04±0.07 2.416±0.088

G16

Tainung67

Japonica

93.3±6.5

1.257±0.01

17.31±0.5 2.988±0.038 4.47±0.05

21.54±0.59 2.128±0.037

G17 Zhenshan97B

Indica

113.7±5.3

1.931±0.024

7.83±0.65 3.486±0.027 7.73±0.30

18.34±0.09

G18

Aswina

Indica

121.8±40.9 1.691±0.058 35.33±3.08 3.474±0.081 6.82±0.89

20.15±1.37 2.241±0.031

G19

Cypress

Japonica

93.9±10.5

20.17±0.03 1.949±0.021

G20

Rayada

Aus

415 416 417 418 419

1.59±0.027

7.67±1.2

7.56±0.28

3.533±0.007 5.15±0.12

2.811±0.09

8.44±0.17

17.39±2.55 3.286±0.038 9.61±0.59

185.1±76.6 1.616±0.013 38.86±0.98 2.952±0.019 7.23±0.36

22.25±0.89 1.351±0.001 18.33±2

1.844±0.249

1.2±0.03

18.06±0.25 1.375±0.019 1.76±0.02

1.929±0.11

2.36±0.014

22.23±1.5

1.615±0.06

Average

112.3

1.653*

19.08

3.232

7.77*

23.19*

1.936*

LSD1

60.0

0.107

5.27

0.181

1.19

2.13

0.240

Data were presented as mean value ± standard deviation. 1 LSD: least significant difference (P < 0.05). * in the average data indicates significant difference between mean of mineral contents in different years (P < 0.05).

18

420

Table 3

421

Mean of mineral contents among different rice subspecies in different years. Ca Na Subgroup Mg (g/kg) K (g/kg) Fe (mg/kg) Zn (mg/kg) Cu (mg/kg) (mg/kg) (mg/kg) 2011 Aus 132.0a 1.721a 18.83b 3.144a 11.048a 18.55a 0.333c Indica 119.0ab 1.729a 14.98b 3.203a 8.136b 18.15a 0.748a Japonica 105.0b 1.705a 25.86a 3.339a 12.388a 16.98a 0.514b 2012 Aus 142.9a 1.622a 21.58ab 3.011b 8.111a 24.49a 1.468b Indica 106.4b 1.657a 13.96b 3.188ab 5.908b 19.56b 2.109a Japonica 103.0b 1.664a 22.95a 3.387a 9.461a 26.18a 1.997a 422 Different letters in the same column indicate significant difference at P < 0.05. 423

19

Table 4 Correlation analysis of mineral contents in two years. Correlation

424 425 426 427 428

Ca

Mg

Na

K

Fe

Zn

Cu

Ca 0.222 0.200 0.013 -0.014 0.265 -0.119 0.562*** Mg 0.237 0.194 0.618*** 0.495** 0.593*** 0.131 0.713*** Na 0.337* 0.161 0.526*** 0.062 -0.094 0.658*** 0.467** K -0.062 0.543*** 0.147 0.425** 0.255 0.024 0.809*** Fe 0.074 0.661*** 0.318* 0.447** 0.150 -0.356* 0.786*** Zn 0.170 0.558*** 0.485** 0.335** 0.804*** 0.455** 0.295 Cu -0.021 0.300 0.223 -0.132 0.077 0.169 0.383* *,**,*** Indicate significance at P < 0.05, 0.01, and 0.001 levels, respectively. Correlation between different mineral elements above diagonal were for data in 2011, those below diagonal were for data in 2012, and those in diagonal in bold were for data between 2011 and 2012.

20

429

Table 5

430 Mean square values from analysis of variance for mineral contents. Source df Ca Mg Na K Fe Zn Cu Genotype (G) 19 2803.234*** 0.126*** 276.232*** 0.531*** 25.525*** 49.085*** 0.369*** Year (E) 1 268.975 0.084*** 21.061 0.004 140.383*** 589.763*** 37.23*** 19 549.653 0.019*** 58.049*** 0.027 1.447 26.941*** 0.188*** G × E 431 *** indicate significance at P <0.001 levels; df: degree of freedom. 432 433

21

434

Table 6

435 436

Significantly associated QTLs for Ca, Cu, K, Na and Zn detected in 2 years by association test. Trait

Ca

Cu K

Na

Zn

Year

QTL

2011 qCa2.1 qCa6.1 2012 qCa2.2 qCa6.2 qCa8 2011 qCu6 2011 qK1 qK6.1 qK6.2 2012 qK2 qK8 qK11 2011 qNa1 qNa11 2012 qNa6 qNa10 2011 qZn2

Chromosome Position (bp)

2 6 2 6 8 6 1 6 6 2 8 11 1 11 6 10 2

30592700 3450780 19703118 2500449 22222152 15357788 41972397 813231 10689396 24689470 27690577 3086965 43155830 4854990 21112917 21698303 23978212

437 438

22

Major

Miner

Miner allele

allele

allele

frequency

C G A C G C C C T C C C T T G T C

T T G T T T T T C A T A C G A A A

P value

0.0046 0.0048 0.0048 0.0043 0.0025 0.0049 0.0046 0.0038 0.0036 0.0042 0.0044 0.0048 0.0049 0.0048 0.0031 0.0031 0.0036

0.325 0.35 0.375 0.425 0.5 0.3 0.3 0.375 0.325 0.325 0.35 0.3 0.45 0.375 0.45 0.45 0.4

Highlights



Variation in mineral elements in whole grain of 20 rices was investigated.



Ca, Na and K were mainly affected by the genotypic variance.



Fe, Zn and Cu were mainly affected by the environment variance.



17 genetic loci were identified to control the mineral contents.