Quantitative trait loci for grain mineral element accumulation in Vietnamese rice landraces

Hien Linh Tran, Giang Thi Hoang, Nhung Thi Phuong Phung, Ham Huy Le,

Alexandre Grondin, Pascal Gantet

► To cite this version:

Hien Linh Tran, Giang Thi Hoang, Nhung Thi Phuong Phung, Ham Huy Le, Alexandre Grondin, et al.. Quantitative trait loci for grain mineral element accumulation in Vietnamese rice landraces. PLoS ONE, 2024, 19 (12), pp.e0315666. 10.1371/journal.pone.0315666. hal-05080508

HAL Id: hal-05080508

https://hal.science/hal-05080508v1

Submitted on 22 May 2025

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire **HAL**, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Quantitative trait loci for grain mineral element accumulation in Vietnamese rice landraces

Hien Linh Tran¹,²*, Giang Thi Hoang²,³*, Nhung Thi Phuong Phung², Ham Huy Le²,³, Alexandre Grondin¹*, Pascal Gantet¹

1 UMR DIADE, IRD, CIRAD, Université de Montpellier, Montpellier, France, 2 National Key Laboratory for

Plant Cell Biotechnology, LMI RICE, Agricultural Genetics Institute, Hanoi, Vietnam, 3 VNU University of

Engineering and Technology, VNU, Hanoi, Vietnam

© These authors contributed equally to this work.

* alexandre.grondin@ird.fr (AG); nuocngamos@yahoo.com (GTH)

<figure><img src="image_2.png" /></figure>

## Abstract

Rice (*Oryza sativa* L.) is a staple food for half of the world's population, and its biofortification is a key factor in fighting micronutrient malnutrition. However, harmful heavy metals tend to accumulate in rice grains due to soil and water contamination. Therefore, it is important to improve beneficial micronutrient contents and reduce the accumulation of undesirable metals in rice grain. To better characterize the genetic control of mineral accumulation in rice, we conducted association genetics on the ion contents of white and brown grains using a collection of 184 Vietnamese rice landraces. In total, 27 significant associations were identified and delimited into quantitative trait loci associated with macronutrients such as phosphorus, potassium or calcium; micronutrients such as iron or zinc; or toxic heavy metals such as arsenic and cadmium. Several genes related to ion homeostasis or ion transport were identified in the different quantitative trait loci. *LOC_Os10g30610*, present in qRAs10-1 associated with arsenic content in brown rice, encodes an ABC transporter (OsABCG25), which is involved in the silicon-induced formation of the Casparian strip in the rice exodermis and could act as a barrier restricting As diffusion within the root cortex. *LOC_Os05g04330*, present in qRP5-1 and associated with phosphorus content in brown rice, encodes a CHH methylation maintenance protein, and its expression is downregulated in roots in the presence of the *phosphorus uptake 1* (*Pup1*), suggesting a role for epigenetics in the regulation of phosphorus uptake and accumulation in grain. These findings reveal novel quantitative trait loci associated with grain ion content and candidate genes that are potentially valuable for breeding programs aimed at rice grain biofortification and reducing toxic metal accumulation.

## Introduction

Rice (*Oryza sativa* L.) is a staple food for more than half of the world's population, serving as a primary source of nutrients for billions of people worldwide. It is a significant contributor of

<figure><img src="image_3.png" /></figure>

### OPEN ACCESS

**Citation:** Tran HL, Hoang GT, Phung NTP, Le HH, Grondin A, Gantet P (2024) Quantitative trait loci for grain mineral element accumulation in Vietnamese rice landraces. PLoS ONE 19(12): e0315666. https://doi.org/10.1371/journal.pone.0315666

**Editor:** Karthikeyan Adhimoolam, Jeju National University, REPUBLIC OF KOREA

**Received:** May 4, 2024

**Accepted:** November 29, 2024

**Published:** December 23, 2024

**Peer Review History:** PLOS recognizes the benefits of transparency in the peer review process; therefore, we enable the publication of all of the content of peer review and author responses alongside final, published articles. The editorial history of this article is available here: https://doi.org/10.1371/journal.pone.0315666

**Copyright:** © 2024 Tran et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** The data underlying the results presented in the study are available at https://doi.org/10.23708/A76RP9.
**Funding:** Hien Linh Tran is supported by the France Excellence Scholarship Program from the French Embassy in Vietnam. This work was supported by the Institut de Recherche pour le Développement (IRD), the "QUARION" EPPN European Plant Phenotyping Network (EPPN) 2020 project, the "RICE-VN" TRENPLIN-ASEAN prize awarded by the French Academy of Science and the French Ministry of Higher Education and Research (2023–2026). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing interests:** The authors have declared that no competing interests exist.

**Abbreviations:** BR, Brown rice; WR, White rice; qRAs1-1, QTL identifier (qR) for ion arsenic (qRAs) located on chromosome 1 (qRAs1) with specific identification number for this ion on chromosome 1 (qRAs1-1); ICP-MS; Inductively coupled plasma mass spectrometry; BLINK, Bayesian-information and linkage-disequilibrium iteratively nested keyway.

calories while also supplying essential minerals such as calcium (Ca), phosphorus (P), and magnesium (Mg) and, to a lesser extent, copper (Cu), iron (Fe), manganese (Mn) and zinc (Zn) [1]. In some rice-consuming countries, the biofortification of Fe and Zn contents constitutes a crucial objective to reduce micronutrient malnutrition and subsequent effects on human health [2]. Fe deficiency impacts mostly young children and pregnant women in Southeast Asia, West Pacific nations and Africa, whereas Zn deficiency is more widespread in all developing countries [3]. Improvements in the Fe content to over 13 mg/kg and the Zn content to over 28 mg/kg in polished rice grain were set as international targets [4]. Although rice contains many beneficial minerals, it also carries the risk of containing harmful metals for human health, such as arsenic (As) and cadmium (Cd). Rice plants grown in contaminated As and Cd soils are affected in their development, yield and other desired market traits such as the head rice recovery percentage, or consumer traits such as the cooking characteristics [5-7]. Furthermore, accumulations of As and Cd in rice are often associated with reductions in other beneficial ions in grains, as these harmful ions can reduce the accessibility of beneficial ions in the soil or share similar uptake and translocation mechanisms. For instance, As can bind with Fe and Mn oxides and competes with silicon (Si) for uptake [8]. Similarly, NATURAL RESISTANCE-ASSOCIATED MACROPHAGE PROTEIN 5 (OsNRAMP5) acts as a major transporter of Cd and Mn in the roots [9, 10] while OsNRAMP1 is another metal transporter that transports both Cd and Fe [11]. Absorption of Cd can also be mediated by IRON-REGULATED TRANSPORTERS 1 and 2 (OsIRT1 and OsIRT2) [12]. Therefore, when enhancing rice biofortification, caution should be taken to avoid trade-offs caused by the accumulation of harmful elements.

Better agronomical practices combined with varietal improvement successfully improved Zn content by over 24 mg/kg in some released biofortified rice in Bangladesh [13]. The limited variation in the contents of some ions, such as Fe, in polished rice impedes biofortification through conventional breeding. The complex genetic and metabolic networks controlling the ion content in grain, which is highly influenced by the environment, represent another limitation [2, 14]. Transgenic approaches represent an alternative to circumvent these constraints, as constitutive expression of genes involved in ion uptake, transport, storage or metabolism may be less sensitive to the environment and surpass natural variation. For example, the overexpression of nicotianamine synthase 2 (NAS2), which catalyzes the synthesis of the divalent metal chelator nicotianamine acid, resulted in increases of up to 19 mg/kg in Fe and 76 mg/kg in Zn in rice endosperm [15]. However, policy regulations on transgenics limit the use of this approach, and it is important to continue exploring the diversity in ion content in rice grain to identify potential donor parents and quantitative trait loci (QTLs) that could be deployed in breeding. The uncoupling of grain biofortification for beneficial nutrients from the accumulation of harmful metals remains possible, especially through the identification of uptake and transport mechanisms specifically related to these harmful metals [16]. Studies have shown that natural variability in As accumulation in rice varieties exists, with varieties accumulating significantly lower amounts of As than others (up to 20–30 times less) [17–19]. Variation in grain Cd accumulation in a rice mapping population allowed the identification of the *qGCd7* QTL, in which the *HEAVY METAL ATPASE 3* (*OsHMA3*) gene is responsible for limiting root-to-shoot translocation of Cd by its selective sequestration within the vacuole [20]. A better understanding of the genetic and physiological determinants of beneficial and harmful ion accumulation in grains is needed to more efficiently improve rice grain nutritional quality.

Biparental populations (for example, from two contrasted parents or multiparent advanced generation intercross population) can target the identification of genetic determinants of only some ions that present contrasted accumulation in parents [21, 22]. More recently, the capacity to genotype rice diversity collections has been useful for studying the diversity in the grain
ionome and identifying associated QTLs and genes through genome-wide association studies (GWAS) [23, 24]. Studies on "Rice Diversity Panel 1", a collection of accessions representative of *O. sativa* global diversity, identified 330 significant SNPs associated with As, Cu, molybdenum (Mo), and Zn accumulation in grain [25, 26]. By performing a GWAS in a different collection, Yang et al. [27] identified 72 QTLs associated with ion variation and provided evidence for the causal roles of three genes, the sodium transporter gene *HIGH AFFINITY K + TRANSPORTER 1;5 (OsHKT1;5)* for sodium (Na), *MOLYBDATE TRANSPORTER 1;1 (OsMOT1;1)* for Mo and *GRAIN NUMBER, PLANT HEIGHT, AND HEADING DATE7 (Ghd7)* for nitrogen (N). Another study using an *O. sativa* spp. indica panel identified the *NICOTIANAMINE SYNTHASE 3 (OsNAS3)* gene as potentially responsible for the effect of a quantitative trait locus (QTL) associated with Zn content in rice grain [28]. This highlights the usefulness of rice panels to reveal complementary QTLs and genes involved in rice grain ion content. However, most of these association studies were conducted separately for brown rice (BR) and white rice (WR), and only a few studies have compared the genetic basis of both brown and white grains collected from the same set of rice plants [29]. Research on QTLs controlling BR and WR grain ion contents measured in the same plants is becoming increasingly important as brown rice consumption rises.

Vietnamese rice landraces are adapted to a range of environments from North to South Vietnam and present an interesting genetic diversity that has been overlooked. One study using 94 black Vietnamese rice accessions identified original associations for the levels of anthocyanin and flavonoids present in whole grains [30]. Moreover, a genotyped collection of Vietnamese rice landraces originating from various agrosystems across Vietnam was previously established and sequenced [31]. A GWAS of this collection revealed original QTLs associated with root traits [32], panicle development [33], leaf biomass [34], drought resistance [35], jasmonate regulation [36], salt resistance [37] and phosphorous deficiency [38]. However, the diversity in grain ion content is unknown in the Vietnamese rice landraces collection. In this study, we aimed at characterizing the genetic diversity for the contents of 23 ions in brown and white rice grains from plants in this Vietnamese rice landrace collection grown in a field in Vietnam. We further performed a GWAS to identify QTLs governing the ion contents in brown and white rice grains.

## Materials and methods

### Plant materials and genotypic data

In this study, 184 rice genotypes were used, including 181 Vietnamese landraces (113 indica, 62 japonica and 6 admixed accessions) [31] and 3 reference genotypes (Nipponbare, Azucena and IR64; S1 Table). These accessions were collected and are available at the Plant Resources Center (PRC, An Khanh commune, Hoai Duc district, Hanoi city, Vietnam), which acts for the Vietnamese government to preserve the seeds of traditional varieties.

### Field experiment

The field experiment was performed in Hai Phong, Vietnam (20° 51' 59" N and 106° 40' 57" E), between July and November 2019. It was conducted in a field normally used for rice cultivation, therefore no special authorization was required. Seeds were sown in seedling beds for two weeks, and the seedlings were further transplanted into an open field under irrigated low-land conditions (S1 Fig). The rice accessions were grouped according to their flowering time: 67 accessions in the early-flowering group (less than 85 days), 72 accessions in the medium-flowering group (85-119 days), 23 accessions in the late-flowering group (120-154 days), and 22 accessions in the very late-flowering group (more than 155 days). The four flowering groups
were divided into blocks (12 blocks in early, 12 blocks in medium, 4 blocks in late, and 4 blocks in very late) with 50 cm spacing between each pair of blocks and 75 cm spacing between flowing groups. Each accession was planted in three plots within their flowering group, which corresponded to replicates. Each plot was 1 m² in area and contained 25 plants separated by 25 cm between and within rows. The grains from all 25 plants per plot were harvested when the plants reached maturity and were air-dried and stored at 4°C. After drying, the inedible outer hulls of the grains were manually removed to obtain brown rice (BR). White rice (WR), corresponding to the grains without the hull and the bran (including the aleurone layer and the cereal germ), was obtained using a huller machine. Ten BR and WR grains from each plot were finely ground separately using a nonmetal grinder (HARIO Ceramic Coffee Mill Slim Plus Grinder MSS-1DTB) before ion content analysis.

## Grain ion phenotyping

The ion contents were measured at the ionomic platform of the University of Nottingham by inductively coupled plasma-mass spectrometry (ICP-MS) following a similar procedure to [39]. Briefly, ground grains were placed in Pyrex test tubes, and their weights were measured. The tubes were complemented with 1 mL of trace metal grade nitric acid Primar Plus (spiked with an indium internal standard). The samples were subsequently subjected to a preliminary digestion process at room temperature in a fume hood for 20 hours. The samples were transferred to dry block heaters and digested at 115°C for 4 hours. After the samples had cooled to room temperature, 1 mL of hydrogen peroxide was introduced into the tubes, and the samples were subjected to an additional 2-hour digestion in a dry block heater at 115°C. The resulting digested samples were then diluted to a final volume of 10 mL with 18.2 MΩcm Milli-Q Direct water. ICP-MS was performed using a PerkinElmer NexION 2000 in collision mode with Helium (He). Twenty-three elements, including lithium (Li), boron (B), sodium (Na), magnesium (Mg), phosphorus (P), sulfur (S), potassium (K), calcium (Ca), manganese (Mn), iron (Fe), cobalt (Co), nickel (Ni), copper (Cu), zinc (Zn), arsenic (As), selenium (Se), rubidium (Rb), strontium (Sr), molybdate (Mo), cadmium (Cd), lead (Pb), titanium (Ti) and chromium (Cr), were analyzed in both brown and white rice samples. A liquid reference material composed of pooled samples of digested grain powders was prepared before the beginning of the sample run and was used throughout the whole run. It was run after every ninth sample in all the ICP-MS sample sets to correct for variation between and within the ICP-MS analysis runs. Single-element standard solutions were used for calibration. Sample concentrations were calculated using an external calibration method in the instrument's software.

## Statistical analysis

The mean and standard deviation of three repetitions as well as the coefficient of variation for each genotype were measured for all the traits using RStudio version 2022.07.1. A two-way ANOVA (using the aov function) was performed to test the effects of replication and genotype. Broad-sense heritability (H²) was calculated as H² = (F value - 1)/F value, where the F value for the genotype effect was extracted from the ANOVA as described by [34]. Correlations between all traits were analyzed using Pearson's correlation test with a confidence level of 0.95. These correlations were visualized using the corrplot function from the R corrplot package. To perform GWAS, BLUEs per genotype for each ion in BR and WR were calculated in the StatgenSTA package [40] using the spatial analysis of field trials with splines (SpATS) function [41] and considering a resolvable incomplete block design according to the following model:

$$ trait = repID + repID : subBlock + genotype + \varepsilon $$
where repetition (repID) is considered a fixed effect and repID:subBlock is considered a ran-
dom effect.

Genome-wide association studies

GWAS was performed with the GAPIT3 package using a previously published genotypic dataset obtained from 184 accessions from the Vietnamese rice collection, which consisted of 21623 SNP markers [31]. The structure of this population was previously analyzed, and six genotypic groups were identified [31]. Associations between genotypes and phenotypes (BLUEs) were performed using different models: the mixed linear model (MLM), generalized linear model (GLM), multilo- cus mixed model (MLMM), conditional masked language modeling (CMLM), enriched com- pressed mixed linear model (ECMLM) and Bayesian-information and linkage-disequilibrium iteratively nested keyway (BLINK). Among these methods, BLINK is considered to have the high- est statistical power and computing efficiency [42]. Therefore, BLINK was chosen to filter signifi- cant SNPs according to the Bonferroni threshold with a significance level of 0.05. Quantile– quantile plots (QQ plots) and Manhattan plots were generated using GAPIT3 [42].

Linkage disequilibrium heatmap and identification of candidate genes

To delimit QTL regions, pairwise linkage disequilibrium was calculated between significant SNPs and the surrounding SNPs using the LD heatmap package in R [43]. The QTL intervals were restricted to a region where the squared allele frequency correlation (r² values) between markers exceeded 0.4. If the previously delimited region was smaller than 50 kb, the QTL region was extended to include a 50 kb range both upstream and downstream of the significant markers. Genes within the QTL were identified using the MSU Rice Genome Annotation Proj- ect Release 7 (http://rice.uga.edu/) and RAPDB version IRGSP-1.0 2023-09-07 (https://rapdb. dna.affrc.go.jp/). The candidate genes in each QTL region were subsequently selected based on potential links between their predicted function and the ion of interest. The expression levels of the selected candidate genes were retrieved from RiceXPro version 3.0 (https://ricexpro. dna.affrc.go.jp/index.html) [44].

Results

Variation in the grain ion contents of Vietnamese rice landraces

To analyze the variability in ion accumulation in BR and WR in the Vietnamese rice collection, the contents of 23 different ions were analyzed in grains collected from a field experiment performed in Vietnam. ICP-MS analyses revealed the presence of Mg, P, S, K, Ca, Mn, Fe, Co, Ni, Cu, Zn, As, Rb, Sr, Mo and Cd, whereas the Li, B, Na, Ti, Cr, Se and Pb concentrations were below the quantification limits. These latter elements were consequently excluded from further analysis. Most ions presented higher values in BR than in WR, except for Mo and Cd (S2 Fig). In both BR and WR, macronutrients (P, K, Mg, S and Ca) were predominant, as expected, with P having the highest concentration in BR (3369 ppm) and S having the highest concentration in WR (1255 ppm). Among the micronutrients, Co had the lowest concentrations in both BR and WR (0.03 ppm and 0.02 ppm, respectively). In BR, K, Mg, P, and S presented the smallest variations within the collection, with coefficients of variation varying between 11% and 12%. In contrast, Co, As, Mo, and Cd showed the greatest variations, with the coefficients of variation ranging from 33% to 47% in both BR and WR. Significant genotypic effects were observed for all ions with high values of broad-sense heritability (≥ 0.71 for BR and ≥ 0.73 for WR; Table 1), indicating that the grain ion content is genetically controlled in the Vietnamese rice collection.
<table><caption>Table 1. Phenotypic variation and broad sense heritability of grain ion content in brown and white rice in the Vietnamese rice collection.</caption><thead><tr><th rowspan="2">Grain type</th><th rowspan="2">Ion</th><th rowspan="2">Mean (ppm)</th><th rowspan="2">SD</th><th>Min</th><th>Max</th><th>CV</th><th rowspan="2">Replication effect (p-value)</th><th rowspan="2">Genotype effect (p-value)</th><th rowspan="2">H<sup>2</sup></th></tr><tr><th>(ppm)</th><th>(ppm)</th><th>(%)</th></tr></thead><tbody><tr><td rowspan="12">Brown rice</td><td>P</td><td>3369.13</td><td>383.26</td><td>2343.73</td><td>4956.05</td><td>11.38</td><td>0.4430</td><td>&lt;0.001</td><td>0.79</td></tr><tr><td>K</td><td>2653.99</td><td>321.99</td><td>1908.69</td><td>3788.00</td><td>12.13</td><td>0.0298</td><td>&lt;0.001</td><td>0.86</td></tr><tr><td>Mg</td><td>1357.14</td><td>158.92</td><td>916.16</td><td>1927.55</td><td>11.71</td><td>&lt;0.001</td><td>&lt;0.001</td><td>0.82</td></tr><tr><td>S</td><td>1310.55</td><td>166.01</td><td>847.77</td><td>1937.29</td><td>12.67</td><td>&lt;0.001</td><td>&lt;0.001</td><td>0.73</td></tr><tr><td>Ca</td><td>101.47</td><td>19.37</td><td>54.19</td><td>169.55</td><td>19.09</td><td>0.8235</td><td>&lt;0.001</td><td>0.71</td></tr><tr><td>Zn</td><td>35.10</td><td>5.76</td><td>20.76</td><td>64.50</td><td>16.41</td><td>&lt;0.001</td><td>&lt;0.001</td><td>0.85</td></tr><tr><td>Rb</td><td>27.45</td><td>5.82</td><td>13.30</td><td>49.93</td><td>21.20</td><td>0.0323</td><td>&lt;0.001</td><td>0.84</td></tr><tr><td>Mn</td><td>23.98</td><td>4.25</td><td>11.51</td><td>35.68</td><td>17.72</td><td>&lt;0.001</td><td>&lt;0.001</td><td>0.84</td></tr><tr><td>Fe</td><td>10.39</td><td>2.04</td><td>6.63</td><td>21.59</td><td>19.63</td><td>&lt;0.001</td><td>&lt;0.001</td><td>0.74</td></tr><tr><td>Cu</td><td>4.02</td><td>0.78</td><td>2.02</td><td>6.80</td><td>19.40</td><td>&lt;0.001</td><td>&lt;0.001</td><td>0.90</td></tr><tr><td>Mo</td><td>1.02</td><td>0.44</td><td>0.24</td><td>3.17</td><td>43.14</td><td>&lt;0.001</td><td>&lt;0.001</td><td>0.98</td></tr><tr><td>Ni</td><td>0.95</td><td>0.28</td><td>0.24</td><td>1.89</td><td>29.47</td><td>&lt;0.001</td><td>&lt;0.001</td><td>0.90</td></tr><tr><td>Sr</td><td>0.27</td><td>0.08</td><td>0.10</td><td>0.69</td><td>29.63</td><td>0.0151</td><td>&lt;0.001</td><td>0.91</td></tr><tr><td>Cd</td><td>0.21</td><td>0.10</td><td>0.02</td><td>0.53</td><td>47.62</td><td>0.0011</td><td>&lt;0.001</td><td>0.95</td></tr><tr><td>As</td><td>0.16</td><td>0.07</td><td>0.01</td><td>0.35</td><td>43.75</td><td>0.1528</td><td>&lt;0.001</td><td>0.78</td></tr><tr><td>Co</td><td>0.03</td><td>0.01</td><td>0.01</td><td>0.10</td><td>33.33</td><td>0.2419</td><td>&lt;0.001</td><td>0.88</td></tr><tr><td rowspan="12">White rice</td><td>S</td><td>1254.91</td><td>168.63</td><td>766.00</td><td>1958.37</td><td>13.44</td><td>&lt;0.001</td><td>&lt;0.001</td><td>0.75</td></tr><tr><td>P</td><td>1246.63</td><td>260.15</td><td>540.82</td><td>2340.61</td><td>20.87</td><td>&lt;0.001</td><td>&lt;0.001</td><td>0.92</td></tr><tr><td>K</td><td>982.53</td><td>237.79</td><td>499.89</td><td>1944.03</td><td>24.20</td><td>0.0836</td><td>&lt;0.001</td><td>0.94</td></tr><tr><td>Mg</td><td>377.94</td><td>109.77</td><td>148.12</td><td>865.10</td><td>29.04</td><td>0.0779</td><td>&lt;0.001</td><td>0.92</td></tr><tr><td>Ca</td><td>54.11</td><td>13.07</td><td>20.30</td><td>101.37</td><td>24.15</td><td>0.0112</td><td>&lt;0.001</td><td>0.80</td></tr><tr><td>Zn</td><td>28.47</td><td>5.31</td><td>16.94</td><td>47.70</td><td>18.65</td><td>&lt;0.001</td><td>&lt;0.001</td><td>0.91</td></tr><tr><td>Mn</td><td>10.77</td><td>2.09</td><td>6.23</td><td>18.72</td><td>19.41</td><td>0.0059</td><td>&lt;0.001</td><td>0.92</td></tr><tr><td>Rb</td><td>9.70</td><td>2.42</td><td>4.59</td><td>20.69</td><td>24.95</td><td>0.0323</td><td>&lt;0.001</td><td>0.89</td></tr><tr><td>Fe</td><td>4.04</td><td>1.56</td><td>1.18</td><td>11.30</td><td>38.61</td><td>&lt;0.001</td><td>&lt;0.001</td><td>0.73</td></tr><tr><td>Cu</td><td>3.39</td><td>0.64</td><td>1.61</td><td>5.92</td><td>18.88</td><td>&lt;0.001</td><td>&lt;0.001</td><td>0.94</td></tr><tr><td>Mo</td><td>1.01</td><td>0.46</td><td>0.23</td><td>3.25</td><td>45.54</td><td>&lt;0.001</td><td>&lt;0.001</td><td>0.98</td></tr><tr><td>Ni</td><td>0.63</td><td>0.22</td><td>0.18</td><td>1.59</td><td>34.92</td><td>&lt;0.001</td><td>&lt;0.001</td><td>0.90</td></tr><tr><td>Cd</td><td>0.21</td><td>0.10</td><td>0.01</td><td>0.53</td><td>47.62</td><td>&lt;0.001</td><td>&lt;0.001</td><td>0.95</td></tr><tr><td>As</td><td>0.11</td><td>0.05</td><td>0.01</td><td>0.29</td><td>45.45</td><td>0.0495</td><td>&lt;0.001</td><td>0.82</td></tr><tr><td>Sr</td><td>0.09</td><td>0.03</td><td>0.03</td><td>0.38</td><td>33.33</td><td>0.0951</td><td>&lt;0.001</td><td>0.89</td></tr><tr><td>Co</td><td>0.02</td><td>0.01</td><td>0</td><td>0.05</td><td>50.00</td><td>&lt;0.001</td><td>&lt;0.001</td><td>0.90</td></tr></tbody></table>

SD: standard deviation; Min: minimum mean value per genotype observed in the collection; Max: maximum mean value per genotype observed in the collection; CV: coefficient variation; H²: broad sense heritability. p-values from two-way ANOVA tests are represented.

https://doi.org/10.1371/journal.pone.0315666.t001

To check for correlations between ions, Pearson's correlation analysis was performed (Fig 1). Moderate to very strong positive correlations were observed between the P, K and Mg contents for both BR and WR (R = 0.54 to R = 0.87 between combinations in BR and WR). In BR, moderate positive correlations were observed between the S and P contents (R = 0.51) and between the Sr and As contents (R = 0.54). In WR, Mn was moderately positively correlated with Mg, K and Sr (R > 0.50), Rb was also moderately positively correlated with Mg and K (R > 0.55), and Mo was moderately positively correlated with Zn (R = 0.50). For both BR and WR, significant negative correlations were observed between the Ni and As contents (R = -0.34 and R = -0.36, respectively), between the Ni and Mo contents (R = -0.46 and R = -0.39, respectively) and between the Cd and Mo contents (R = -0.45 and R = -0.47, respectively).
<figure><img src="image_4.png" /><figcaption>**Fig 1.** Correlation between ion content in brown rice (A) and white rice (B) grains. Color gradients indicate the Pearson's correlation coefficient. Non-significant correlations at a p-value threshold of 0.05 are indicated with a cross.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0315666.g001

These results indicate strong associations for the ion content in BR and WR (e.g., P, K and Mg).

## Identification of QTLs associated with grain ion content in Vietnamese rice landraces

The high diversity and high heritability observed in the Vietnamese rice collection for ion content in both BR and WR suggest that these traits are amenable for association genetic analysis. To identify genomic regions associated with grain ion content in Vietnamese rice landraces, we carried out a GWAS on the 16 detected ions in BR and WR rice using a previously described genotypic dataset composed of 21623 SNP markers [31]. We identified a total of 27 associations considering the Bonferroni threshold (Table 2 and S3 Fig). Among them, 14 associations passed the significance threshold in BLINK and other GWAS models, and 13 exclusively passed the significance threshold in the BLINK model. Significant associations on chromosome 7 for Cu and on chromosome 10 for As were observed in both BR and WR. Additionally, a similar association was observed on chromosome 11 for both Cu and Fe in BR (Table 2).

Pairwise linkage disequilibrium was calculated between significant SNPs to define the size of each association, which was further referred to as a QTL (Table 2). We investigated whether the QTLs identified in our study colocalized with previously identified QTLs for root and panicle architecture and development, and salt and drought stress tolerance-related traits in the same Vietnamese rice collection [32–35, 37]. We found that qRMo1-1 overlaps with QTL q48, which is related to the number of crown roots per tiller. qRAs1-1 and qRCa4-1 also overlap with QTL_2 and QTL_11, which are related to the ratio of chlorophyll a to chlorophyll b in leaves [37]. qRK10-2 and qRCa4-1 overlap with a QTL associated with leaf mass and yield-related traits [33, 34]. Finally, qRCd11-1 overlaps with QTLs associated with drought tolerance-related traits [35]. In addition, we looked for colocalization between QTLs for grain ion content identified in the Vietnamese rice collection and QTLs for grain ion content identified in other rice populations in the literature. On chromosome 1, qRAs1-1 overlaps with qAs1-3, a QTL for grain As content identified in a rice diversity panel composed of 202 accessions from
<table><caption>Table 2. List and position of QTLs identified in the Vietnamese collection for ion content in brown rice (BR) and white rice (WR) grains.</caption><thead><tr><th>Ion</th><th>QTL</th><th>Grain type</th><th>Chr</th><th>SNP</th><th>SNP position</th><th>QTL position</th></tr></thead><tbody><tr><td>As</td><td>qRAs1-1</td><td>BR</td><td>1</td><td>Sj01_33174934R</td><td>33174934</td><td>33060660-33174934</td></tr><tr><td>As</td><td>qRAs4-1</td><td>WR</td><td>4</td><td><b>Dj04_03138652R</b></td><td>3138652</td><td>2977218-3143723</td></tr><tr><td>As</td><td>qRAs4-2</td><td>WR</td><td>4</td><td>Dj04_32371574F</td><td>32371574</td><td>32347540-32371574</td></tr><tr><td>As</td><td>qRAs10-1</td><td>BR/WR</td><td>10</td><td><b>Sj10_15968777F</b></td><td>15968777</td><td>15873781-16018764</td></tr><tr><td>As</td><td>qRAs11-1</td><td>BR</td><td>11</td><td><b>Dj11_05821648F</b></td><td>5821648</td><td>5747401-5875727</td></tr><tr><td>Ca</td><td>qRCa3-1</td><td>BR</td><td>3</td><td><b>Dj03_29401407R</b></td><td>29401407</td><td>29303078-29439457</td></tr><tr><td>Ca</td><td>qRCa4-1</td><td>BR</td><td>4</td><td><b>Sj04_04361919F</b></td><td>4361919</td><td>4307480-4399327</td></tr><tr><td>Ca</td><td>qRCa6-1</td><td>BR</td><td>6</td><td><b>Sj06_24039978F</b></td><td>24039978</td><td>23989978-24089978</td></tr><tr><td>Ca</td><td>qRCa9-1</td><td>BR</td><td>9</td><td><b>Dj09_15688119R</b></td><td>15688119</td><td>15638119-15738119</td></tr><tr><td>Cd</td><td>qRCd10-1</td><td>BR</td><td>10</td><td>Dj10_02130045F</td><td>2130045</td><td>2081167-2151377</td></tr><tr><td>Cd</td><td>qRCd11-1</td><td>BR</td><td>11</td><td>Dj11_06632630F</td><td>6632630</td><td>6582630-6682630</td></tr><tr><td>Cu</td><td>qRCu7-1</td><td>BR/WR</td><td>7</td><td><b>Sj07_25260997F</b></td><td>25260997</td><td>25232425-25360187</td></tr><tr><td>Cu/Fe</td><td>qRCuFe11-1</td><td>BR</td><td>11</td><td>Dj11_20389966F</td><td>20389966</td><td>20336015-20407398</td></tr><tr><td>Fe</td><td>qRFe3-1</td><td>WR</td><td>3</td><td><b>Dj03_28302202R</b></td><td>28302202</td><td>28109488-28355049</td></tr><tr><td>Fe</td><td>qRFe10-1</td><td>WR</td><td>10</td><td>Sj10_19727114F</td><td>19727114</td><td>19727114-19788660</td></tr><tr><td>K</td><td>qRK5-1</td><td>BR</td><td>5</td><td>Dj05_03417694F</td><td>3417694</td><td>3365689-3504647</td></tr><tr><td>K</td><td>qRK5-2</td><td>BR</td><td>5</td><td>Dj05_05422273F</td><td>5422273</td><td>5390697-5451215</td></tr><tr><td>K</td><td>qRK10-1</td><td>BR</td><td>10</td><td>Dj10_04097517F</td><td>4097517</td><td>4046874-4097937</td></tr><tr><td>K</td><td>qRK10-2</td><td>BR</td><td>10</td><td><b>Sj10_18059798F</b></td><td>18059798</td><td>17948827-18064123</td></tr><tr><td>Mn</td><td>qRMn4-1</td><td>BR</td><td>4</td><td>Sj04_25552228F</td><td>25552228</td><td>25398103-25552228</td></tr><tr><td>Mo</td><td>qRMo1-1</td><td>WR</td><td>1</td><td><b>Sj01_28315371F</b></td><td>28315371</td><td>27868310-28390748</td></tr><tr><td>P</td><td>qRP5-1</td><td>BR</td><td>5</td><td><b>Dj05_01901555F</b></td><td>1901555</td><td>1884166-1970447</td></tr><tr><td>S</td><td>qRS7-1</td><td>BR</td><td>7</td><td>Dj07_21988531R</td><td>21988531</td><td>21896842-21988531</td></tr><tr><td>S</td><td>qRS12-1</td><td>BR</td><td>12</td><td>Dj12_6190592R</td><td>6190592</td><td>6136693-6190592</td></tr><tr><td>Sr</td><td>qRSr3-1</td><td>WR</td><td>3</td><td><b>Sj03_14341377R</b></td><td>14341377</td><td>14266371-14341377</td></tr><tr><td>Sr</td><td>qRSr8-1</td><td>WR</td><td>8</td><td>Sj08_18311304F</td><td>18311304</td><td>18261304-18361304</td></tr><tr><td>Zn</td><td>qRZn1-1</td><td>WR</td><td>1</td><td><b>Sj01_03956268F</b></td><td>3956268</td><td>3956268-4153223</td></tr></tbody></table>

Chr: Chromosome; BR: brown rice; WR: white rice. SNP in bold font passed the significance threshold in BLINK and other GWAS models. SNP in regular font passed the significance threshold in the BLINK model only.

https://doi.org/10.1371/journal.pone.0315666.t002

<table><caption>Table 3. List of candidate genes identified in QTLs associated with P, K and As content in brown rice (BR) and white rice (WR) grains.</caption><thead><tr><th>QTL</th><th>Ion</th><th>Chr</th><th>Grain type</th><th>Gene ID</th><th>Gene annotation -<br>hypothetical function</th><th>Reference</th></tr></thead><tbody><tr><td>qRP5-1</td><td>P</td><td>5</td><td>BR</td><td>LOC_Os05g04330/<br>Os05g0133900</td><td>DNA methyltransferase protein, DRM3/possibly regulated by Pup1 in root under P starvation</td><td>[47]</td></tr><tr><td>qRK5-1</td><td>K</td><td>5</td><td>BR</td><td>LOC_Os05g06660/<br>Os05g0158500</td><td>Putative Serine Carboxypeptidase homologue, OsSCP26/upregulated under K starvation</td><td>[48]</td></tr><tr><td>qRAs10-1</td><td>As</td><td>10</td><td>BR—WR</td><td>LOC_Os10g30610/<br>Os10g0442900</td><td>ABC transporter, OsABCG25/involved in casparian strip formation mediated by Si</td><td>[49]</td></tr><tr><td>qRAs11-1</td><td>As</td><td>11</td><td>BR</td><td>LOC_Os11g10510/<br>Os11g0210500</td><td>alcohol dehydrogenase, OsADH2/indirectly regulates As uptake by roots</td><td>[52]</td></tr></tbody></table>

Chr: Chromosome.

https://doi.org/10.1371/journal.pone.0315666.t003
<figure><img src="image_5.png" /><figcaption>**Fig 2. Genome-wide association mapping of P content in brown rice grain measured in the Vietnamese rice collection.** A. Manhattan plot representing the p-values of associations between SNPs and grain P content using the BLINK model in GAPIT. The solid green line represents the Bonferroni threshold. B. Focus on the association between SNP Dj05_01901555F and grain P content on chromosome 5. C. Definition of QTL qRP5-1 based on blocks of linkage disequilibrium around Dj05_01901555F. The genes present in the QTL region are represented by grey rectangles. The selected candidate gene *LOC_Os05g04330* is highlighted in red. D. Box plots representing the polymorphism distribution at SNP Dj05_01901555F in the Vietnamese collection. The number of genotypes in each haplotype group is represented. E. Expression levels of *LOC_Os05g04330* at different developmental stages. Expression data in different plant tissues were retrieved from the Rice Xpro database [44] as follow; LB_V: Leaf blade at vegetative stage; LB_Rp: Leaf blade at reproductive stage; LB_R: Leaf blade at ripening stage; LS_V: Leaf sheath at vegetative stage; LS_Rp: Leaf sheath at reproductive stage; R_V: Root at vegetative stage; R_Rp: Root at reproductive stage; S_Rp: Stem at reproductive stage; S_R: Stem at ripening stage; I: Inflorescence of 3.0–4.0 mm; A: Anther of 0.7–1.0 mm; Pis: Pistil of 10–14 cm panicle; L: Lemma of 4.0–5.0 mm floret; Pal: Palea of 4.0–5.0 mm floret; O: Ovary 3 days after flowering; Emb: Embryo at 14 days after flowering; End: Endosperm at 14 days after flowering.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0315666.g002

the USDA-ARS Rice MiniCore collection [45]. Furthermore, on chromosome 11, qRCd11-1 (6582630–6682630 bp) is positioned close to qCd11-1 (6233769–6354200 bp), a QTL associated with grain Cd content identified in a subset of 698 varieties present in the 3K rice genome project [46].

## Potential candidate genes related to grain ion content in Vietnamese rice landraces

Genes present in each QTL were investigated (S2 Table) for their putative functions in determining the grain ion content based on their functional annotations and a literature survey. Five potential candidate genes were selected (Table 3). In qRP5-1 associated with grain P content, *LOC_Os05g04330*, which encodes a DNA CHH methylation maintenance protein, is
<figure><img src="image_6.png" /><figcaption>**Fig 3. Genome-wide association mapping of As content in brown and white rice grains measured in the Vietnamese rice collection.** A. Manhattan plot representing the p-values of associations between SNPs and grain As content using the BLINK model in GAPIT. The solid green line represents the Bonferroni threshold and the dotted green line represents the FDR threshold. B. Focus on the association between SNP Sj10_15968777F and grain As content on chromosome 10. C. Definition of QTL qRAs10-1 based on blocks of linkage disequilibrium around Sj10_15968777F. The genes present in the QTL region are represented by grey rectangles. The selected candidate gene *LOC_Os10g30610* is highlighted in red. D. Box plots representing the polymorphism distribution at SNP Sj10_15968777F in the Vietnamese collection. The number of genotypes in each haplotype group is represented. E. Expression levels of *LOC_Os10g30610* at different developmental stages. Expression data in different plant tissues were retrieved from the Rice Xpro database [44] as follow; LB_V: Leaf blade at vegetative stage; LB_Rp: Leaf blade at reproductive stage; LB_R: Leaf blade at ripening stage; LS_V: Leaf sheath at vegetative stage; LS_Rp: Leaf sheath at reproductive stage; R_V: Root at vegetative stage; R_Rp: Root at reproductive stage; S_Rp: Stem at reproductive stage; S_R: Stem at ripening stage; I: Inflorescence of 3.0–4.0 mm; A: Anther of 0.7–1.0 mm; Pis: Pistil of 10–14 cm panicle; L: Lemma of 4.0–5.0 mm floret; Pal: Palea of 4.0–5.0 mm floret; O: Ovary 3 days after flowering; Emb: Embryo at 14 days after flowering; End: Endosperm at 14 days after flowering.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0315666.g003

located 6404 bp downstream of the significant Dj05_01901555F SNP marker (Fig 2A-2C). Haplotype analysis of this SNP revealed that the T allele present in a minority of landraces was associated with increased grain P content (3.6% increase on average; Fig 2D). *LOC_Os05g04330* is highly expressed in roots during the vegetative stage as well as in anthers and pistils under optimal conditions [44] (Fig 2E). Under P-starvation stress, downregulation of this gene has been observed in a near-isogenic line harboring the *phosphorus uptake 1* (*Pup1*) QTL [47]. *LOC_Os05g06660*, located in qRK5-1, also appeared interesting, as it encodes a serine carboxypeptidase homologous to OsSCL26 (Table 3 and S4 Fig). OsSCL26 plays an important role in the homeostasis of P in rice [48], an ion that was highly correlated with K in our study.
Several genes were also linked to the grain As content in qRAs10-1 and qRAs11-1. In qRA10-1, LOC_Os10g30610, located 39981 bp upstream of the significant SNP marker Sj10_15968777F, encodes an ABC transporter (OsABCG25), which is involved in the silicon-induced formation of the Casparian strip in the rice exodermis (Fig 3A–3C and Table 3) [49]. Haplotype analysis revealed that allele T is present in the majority of Vietnamese rice landraces and provides a lower grain As content than the C allele (24% reduction on average; Fig 3D). OsABCG25 is expressed mainly in the roots (Fig 3E). This gene is upregulated under high-As treatment and may be involved in the homeostasis of As in rice [50]. In qRAs11-1, LOC_Os11g10510 encodes ALCOHOL DEHYDROGENASE 2 (OsADH2), which regulates silicate transporters also acting as As transporters in rice (S5 Fig and Table 3) [51, 52]. Overall, some of these genes represent new potential candidates for P uptake and exclusion of heavy metals such as As in rice plants and ultimately grains.

Discussion

Improving the quality of rice grain requires improving their mineral content. However, while some minerals are important for human health, others can be toxic. Many efforts have focused on the biofortification of rice grain with Fe and Zn, but further increases in these nutrients in rice grain are needed to fight malnutrition worldwide. Therefore, it is important to characterize the diversity of ion contents in rice grain, identify the genetic determinants controlling grain ion accumulation and determine the physiological significance of these mechanisms to avoid trade-offs related to the accumulation of toxic heavy metals. To contribute to this objective, we assessed the contents of 16 different ions in BR and WR grains within a collection of 184 sequenced Vietnamese rice landraces and observed a large diversity for most ions. A GWAS identified 27 QTLs associated with the amounts of important macronutrients, such as P or K; micronutrients, such as Fe or Zn; or toxic elements, such as As and Cd. These QTLs contain relevant candidate genes related to ion transport or ion homeostasis adjustment in rice.

In the Vietnamese rice collection, P, K, Mg and S were the most abundant ions in both BR and WR grains, whereas Cd, As, Sr and Co were among the least abundant. The ion contents observed in our study were of the same order of magnitude as those reported in previous rice grain ionome studies [24, 28, 53]. Interestingly, although most ions were more concentrated in BR than in WR, no significant differences in ion content were observed between the two types of rice for Mo and Cd. This was previously observed for Cd [54, 55] and suggests that Cd, and most likely Mo, is not significantly lost during the milling process because of its accumulation in the endosperm rather than in the outer husk and bran layers [55]. The variations in the grain ion contents in the BR and WR samples also aligned with findings from previous studies [24, 28], where K, Mg, and P had the smallest variations, whereas Fe, Zn, As, Mo, and Cd had the largest variations. In our study, Fe and Zn varied from 6.63 to 21.59 ppm and from 20.76 to 64.50 ppm, respectively, in the BR grains. These values are higher than the maximum Fe (16.58 ppm) content observed in BR grains under flooded conditions in the USDA Core Collection by [24]. This indicates that the QTLs identified in our study constitute a good resource for rice grain biofortification toward the international targets (13 to 15 ppm for Fe and 28 ppm for Zn, considering rice daily consumption in countries where rice is a staple food) [4]. Vietnamese landraces with high grain Fe contents may therefore represent interesting donors for breeding programs aimed at rice grain Fe biofortification [56, 57]. Given that grain ion content is influenced by environmental factors whose frequencies are increasing due to climate change [58], it is important to develop climate-resilient rice varieties while ensuring the maintenance of grain quality.
High positive correlations were observed in both BR and WR between the chemical analogs K and Rb, which are known to share similar uptake mechanisms in plants [59]. Ca and Sr, which are also chemical analogs, were not correlated in either BR or WR in our study, possibly because of the low Sr content in the soil. However, high correlations were observed between P, K, and Mg in BR grains, although these ions do not share similar uptake mechanisms and are not tightly chemically coupled in the soil solution [60]. High correlations between P, K, and Mg were previously reported for BR in a set of widely diverse rice germplasms [24] as well as in other species, such as *Arabidopsis thaliana* [61] and maize [62]. The link between these ions in rice grains may be driven by phytate levels, as most of the P in rice grain is found in the form of mixed K-Mg salts of phytic acid in the aleurone layer and germ [24, 63]. Further evidence confirming the robustness of this correlation came from the colocalization of QTLs for P, K and Mg contents in BR grains in a biparental population [64]. Interestingly, we observed that the correlations between these three ions remained present in WR grains, indicating that other mechanisms independent of phytates within the aleurone layer coregulate P, K and Mg contents in the endosperm. Moreover, significant negative correlations between As and Ni contents in BR and WR grains were observed in our study, as previously described by Williams et al. [65]. Significant negative correlations were also observed between Mo and Cd in both BR and WR grains. In line with these results, it has been reported that Mo alleviates Cd uptake, translocation and toxicity in rice [66] and other species, such as *Brassica napus* [67]. The mechanisms involved in this response remain elusive but may be related to the beneficial physiological effects of Mo, which balance Cd stress rather than competition for uptake [66, 67]. Overall, our results support previous observations while providing further indications of the links between ion accumulation within different rice grain tissues, i.e., the aleurone layer versus the endosperm.

The GWAS revealed 27 significant SNPs associated with BR and WR grain ion contents. We extended these associations to QTL regions using linkage disequilibrium as a criterion to group SNPs together. The QTL qRAs1-1 identified in this study for As overlapped with the QTL qAs1-3, which was previously associated with the same traits but was identified in different rice populations [45]. Furthermore, the QTL qRCd11-1 is located close to another QTL associated with Cd in another rice population [46]. This underscores the robustness of these two QTLs and additionally validates our association study to some extent. Interestingly, five QTLs identified in this study (qRMo1-1, qRAs11-1, qRCa4-1, qRK10-2 and qRCd11-1) were found to overlap with QTLs for morphological and functional traits identified in prior studies using the same rice collection. A QTL for Mo (qRMo1-1) was found to overlap with q48 and was associated with the number of crown roots per tiller [32]. Mo uptake occurs predominantly in roots via MO TRANSPORTER 1 (MOT1) [68] and contributes to the structuration and activity of enzymes such as nitrate reductase [69]. The colocalization of qRMo1-1 and q48 suggests that the number of crown roots is important for the accumulation of Mo and ultimately its translocation to grain, likely via better Mo uptake capacity through enhanced soil exploration. This study also identified a QTL for As content in BR (qRAs11-1), in which we identified LOC_Os11g10510, a gene that encodes an alcohol dehydrogenase named OsADH2. The low-arsenic line 3 (las3) mutant, in which OsADH2 is not functional, presents reduced amounts of As in aerial tissues and grains [52]. The decrease in As content in las3 is correlated with the downregulation of genes encoding root silicate transporters (LOW SILICON 1 and 2), which are also known to play major roles in As uptake [70]. It has been proposed that under anaerobic conditions, the loss of function of OsADH2 leads to a decrease in ethanol fermentation, inducing a subsequent decrease in the cytosolic pH of root cells, which represses LSI1 and LSI2 expression and As uptake [52]. It is therefore likely that OsADH2 and its alleles are responsive to the effect of qRA11-1 on As grain accumulation.
Among the new QTLs identified in this study, we identified several interesting candidate genes whose putative functions are related to ion uptake or homeostasis, and which could consequently affect the grain ion content. In qRAs10-1, for example, *LOC_Os10g30610* is upregulated in response to As [50]. This gene encodes an ABC transporter involved in exodermis suberization [49], that act as a barrier restricting ion diffusion such as Fe within the root cortex [71]. In fact, the overexpression of *LOC_Os10g30610* increased the extent of Casparian strips at the exodermis [71] and may contribute to reduce the inflow of As within the root. Therefore, the putative function of qRAs10-1 in As exclusion and qRAs11-1 in reduction of As transport may have a synergistic effect on lowering As uptake in rice and accumulation in grain, if combined. Furthermore, *LOC_Os05g04330* in qRP5-1 is involved in CHH DNA methylation maintenance, and the regulation of its expression in roots is dependent on the major phosphate starvation tolerance *Pup1* QTL [47]. These findings suggest that the response of rice to phosphate deficiency involves the regulation of gene expression via the modulation of DNA methylation. The involvement of these candidate genes in the regulation of grain ion accumulation needs further validation through functional analysis.

In conclusion, our association study on ion contents in BR and WR grains from a Vietnamese rice landrace collection confirmed several previously identified QTLs, provided new insights into the physiological mechanisms of some of these QTLs and identified new QTLs containing interesting potential candidate genes. Vietnamese landraces carrying QTLs for high Fe and Zn grain contents may represent new potential donors for rice grain biofortification, as these contents exceed the values targeted by breeding programs. Moreover, other donors carrying QTLs with complementary effects may be pyramided to reduce the As ion content in grain. Overall, the identification of these QTLs is valuable for improving the mineral ion content of rice grain while reducing the accumulation of toxic heavy metals.

**Supporting information**

**S1 File. Inclusivity in global research.**

(DOCX)

**S1 Fig. Field trial planting in Hai Phong in 2019.**

(TIF)

**S2 Fig. Box plot representing variation in ion content in brown rice (BR) and white rice (WR) in the Vietnamese rice collection.** Ion content is represented as ppm. p-values from ANOVA tests are represented.

(TIF)

**S3 Fig. Localization of QTLs associated with ion content in Vietnamese rice.** The 12 rice chromosomes are represented. The vertical scale indicates the size of chromosomes in mega base pairs (Mb).

(TIF)

**S4 Fig. Genome-wide association mapping of K content in brown rice grain measured in the Vietnamese rice collection.** A. Manhattan plot representing the p-values of associations between SNPs and grain K content using the BLINK model in GAPIT. The solid green line represents the Bonferroni threshold and the dotted green line represents FDR threshold. B. Focus on the association between SNP Dj05_03417694F and grain K content on chromosome 5. C. Definition of QTL qRK5-1 based on blocks of linkage disequilibrium around Dj05_03417694F. The genes present in the QTL region are represented by grey rectangles. The selected candidate gene *LOC_Os05g06660* is highlighted in red. D. Box plots representing the
polymorphism distribution at SNP Dj05_03417694F in the Vietnamese collection. The num-
ber of genotypes in each haplotype group is represented. E. Expression levels of
LOC_Os05g06660 at different developmental stages. Expression data in different plant tissues
were retrieved from the Rice Xpro database [44] as follow; LB_V: Leaf blade at vegetative stage;
LB_Rp: Leaf blade at reproductive stage; LB_R: Leaf blade at ripening stage; LS_V: Leaf sheath
at vegetative stage; LS_Rp: Leaf sheath at reproductive stage; R_V: Root at vegetative stage;
R_Rp: Root at reproductive stage; S_Rp: Stem at reproductive stage; S_R: Stem at ripening
stage; I: Inflorescence of 3.0–4.0 mm; A: Anther of 0.7–1.0 mm; Pis: Pistil of 10–14 cm panicle;
L: Lemma of 4.0–5.0 mm floret; Pal: Palea of 4.0–5.0 mm floret; O: Ovary 3 days after flower-
ing; Emb: Embryo at 14 days after flowering; End: Endosperm at 14 days after flowering.
(TIF)

**S5 Fig. Genome-wide association mapping of As content in brown rice grain measured in the Vietnamese rice collection.**

A. Manhattan plot representing the p-values of associations between SNPs and grain As content using the BLINK model in GAPIT. The solid green line represents the Bonferroni threshold and the dotted green line represents FDR threshold. B. Focus on the association between SNP Dj11_05821648F and grain As content on chromosome 11. C. Definition of QTL qRAs11-1 based on blocks of linkage disequilibrium around Dj11_05821648F. The genes present in the QTL region are represented by grey rectangles. The selected candidate gene *LOC_Os11g10510* is highlighted in red. D. Box plots representing the polymorphism distribution at SNP Dj11_05821648F in the Vietnamese collection. The number of genotypes in each haplotype group is represented. E. Expression levels of *LOC_Os11g10510* at different developmental stages. Expression data in different plant tissues were retrieved from the Rice Xpro database [44] as follow; LB_V: Leaf blade at vegetative stage; LB_Rp: Leaf blade at reproductive stage; LB_R: Leaf blade at ripening stage; LS_V: Leaf sheath at vegetative stage; LS_Rp: Leaf sheath at reproductive stage; R_V: Root at vegetative stage; R_Rp: Root at reproductive stage; S_Rp: Stem at reproductive stage; S_R: Stem at ripening stage; I: Inflorescence of 3.0–4.0 mm; A: Anther of 0.7–1.0 mm; Pis: Pistil of 10–14 cm panicle; L: Lemma of 4.0–5.0 mm floret; Pal: Palea of 4.0–5.0 mm floret; O: Ovary 3 days after flowering; Emb: Embryo at 14 days after flowering; End: Endosperm at 14 days after flowering. (TIF)

**S1 Table.** Passport data of the different landraces present in the Vietnamese rice collection.
I: indica; J: japonica; m: admix; na: no data available.
(XLSX)

**S2 Table.** List of genes present in each QTL identified foreign content in the Vietnamese rice collection.
(XLSX)

**Acknowledgments**

The authors would like to thank Nguyen Van Toan, Vu Thi Huong and Vu Manh An for their
valuable support in field care and rice grain sample harvesting.

**Author Contributions**

**Conceptualization:** Giang Thi Hoang, Alexandre Grondin, Pascal Gantet.

**Data curation:** Hien Linh Tran, Giang Thi Hoang, Nhung Thi Phuong Phung, Ham Huy Le.

**Formal analysis:** Hien Linh Tran, Giang Thi Hoang, Alexandre Grondin, Pascal Gantet.
**Funding acquisition:** Giang Thi Hoang, Alexandre Grondin, Pascal Gantet.

**Investigation:** Hien Linh Tran, Giang Thi Hoang.

**Project administration:** Giang Thi Hoang, Pascal Gantet.

**Supervision:** Giang Thi Hoang, Alexandre Grondin, Pascal Gantet.

**Writing – original draft:** Hien Linh Tran, Giang Thi Hoang, Alexandre Grondin, Pascal Gantet.

**Writing – review & editing:** Hien Linh Tran, Giang Thi Hoang, Nhung Thi Phuong Phung, Ham Huy Le, Alexandre Grondin, Pascal Gantet.

## References

1. Verma D, Srivastav P. Proximate composition, mineral content and fatty acids analyses of aromatic and non-aromatic Indian rice. Rice Science. 2017; 24: 21–31. https://doi.org/10.1016/j.rsci.2016.05.005

2. Senguttuvel P, Padmavathi G, Jasmine C, Sanjeeva Rao D, Neeraja C, Jaldhani V, et al. Rice biofortification: breeding and genomic approaches for genetic enhancement of grain zinc and iron contents. Front Plant Sci. 2023; 14: 1138408. https://doi.org/10.3389/fpls.2023.1138408 PMID: 37332714

3. Ramakrishnan U. Prevalence of micronutrient malnutrition worldwide. Nutrition Reviews. 2002; 60: S46–S52. https://doi.org/10.1301/00296640260130731 PMID: 12035858

4. Trijatmiko KR, Dueñas C, Tsakiralpoglou N, Torrizo L, Arines FM, Adeva C, et al. Biofortified indica rice attains iron and zinc nutrition dietary targets in the field. Sci Rep. 2016; 6: 19792. https://doi.org/10.1038/srep19792 PMID: 26806528

5. Hussain B, Ashraf MN, Shafeeq-ur-Rahman, Abbas A, Li J, Farooq M. Cadmium stress in paddy fields: Effects of soil conditions and remediation strategies. Science of The Total Environment. 2021; 754: 142188. https://doi.org/10.1016/j.scitotenv.2020.142188 PMID: 33254942

6. Moulick D, Santra SC, Ghosh D. Consequences of paddy cultivation in arsenic-contaminated paddy fields of lower indo-gangetic plain on arsenic accumulation pattern and selected grain quality traits: A preliminary assessment. In: Hasanuzzaman M, Nahar K, Fujita M, editors. Mechanisms of arsenic toxicity and tolerance in plants. Singapore: Springer; 2018. pp. 49–78. https://doi.org/10.1007/978-981-13-1292-2_3

7. Moulick D, Ghosh D, Skalicky M, Gharde Y, Mazumder MK, Choudhury S, et al. Interrelationship among rice grain arsenic, micronutrients content and grain quality attributes: An investigation from genotype × environment perspective. Front Environ Sci. 2022; 10: 857629. https://doi.org/10.3389/fenvs.2022.857629

8. Ma JF, Yamaji N, Mitani N, Xu X-Y, Su Y-H, McGrath SP, et al. Transporters of arsenite in rice and their role in arsenic accumulation in rice grain. Proceedings of the National Academy of Sciences. 2008; 105: 9931–9935. https://doi.org/10.1073/pnas.0802361105 PMID: 18626020

9. Ishikawa S, Ishimaru Y, Igura M, Kuramata M, Abe T, Senoura T, et al. Ion-beam irradiation, gene identification, and marker-assisted breeding in the development of low-cadmium rice. Proceedings of the National Academy of Sciences. 2012; 109: 19166–19171. https://doi.org/10.1073/pnas.1211132109 PMID: 23132948

10. Sasaki A, Yamaji N, Yokosho K, Ma JF. Nramp5 is a major transporter responsible for manganese and cadmium uptake in rice. The Plant Cell. 2012; 24: 2155–2167. https://doi.org/10.1105/tpc.112.096925 PMID: 22589467

11. Takahashi R, Ishimaru Y, Senoura T, Shimo H, Ishikawa S, Arao T, et al. The OsNRAMP1 iron transporter is involved in Cd accumulation in rice. J Exp Bot. 2011; 62: 4843–4850. https://doi.org/10.1093/jxb/err136 PMID: 21697258

12. Nakanishi H, Ogawa I, Ishimaru Y, Mori S, Nishizawa NK. Iron deficiency enhances cadmium uptake and translocation mediated by the Fe2+ transporters OsIRT1 and OsIRT2 in rice. Soil Science and Plant Nutrition. 2006; 52: 464–469. https://doi.org/10.1111/j.1747-0765.2006.00055.x

13. Kader MdA Haq ME, Anisuzzaman M, Hore T, Biswas P, Aditya T. Zinc enriched high yielding rice variety BRRI dhan84 for dry season rice growing areas of Bangladesh. Asian Plant Research Journal. 2020; 6: 6–13. https://doi.org/10.9734/APRJ/2020/v6i130117

14. Sharif MK, Butt MS, Anjum FM, Khan SH. Rice bran: A novel functional ingredient. Critical Reviews in Food Science and Nutrition. 2014; 54: 807–816. https://doi.org/10.1080/10408398.2011.608586 PMID: 24345050
15. Johnson AAT, Kyriacou B, Callahan DL, Carruthers L, Stangoulis J, Lombi E, et al. Constitutive overexpression of the OsNAS gene family reveals single-gene strategies for effective iron- and zinc-biofortification of rice endosperm. PLOS ONE. 2011; 6: e24476. https://doi.org/10.1371/journal.pone.0024476 PMID: 21915334

16. Slamet-Loedin IH, Johnson-Beebout SE, Impa S, Tsakirpaloglou N. Enriching rice with Zn and Fe while minimizing Cd risk. Front Plant Sci. 2015; 6: 121. https://doi.org/10.3389/fpls.2015.00121 PMID: 25814994

17. Norton GJ, Pinson SRM, Alexander J, Mckay S, Hansen H, Duan G-L, et al. Variation in grain arsenic assessed in a diverse panel of rice (Oryza sativa) grown in multiple sites. New Phytologist. 2012; 193: 650–664. https://doi.org/10.1111/j.1469-8137.2011.03983.x PMID: 22142234

18. Syu C-H, Huang C-C, Jiang P-Y, Lee C-H, Lee D-Y. Arsenic accumulation and speciation in rice grains influenced by arsenic phytotoxicity and rice genotypes grown in arsenic-elevated paddy soils. Journal of Hazardous Materials. 2015; 286: 179–186. https://doi.org/10.1016/j.jhazmat.2014.12.052 PMID: 25577320

19. Zhang X, Wu S, Ren B, Chen B. Water management, rice varieties and mycorrhizal inoculation influence arsenic concentration and speciation in rice grains. Mycorrhiza. 2016; 26: 299–309. https://doi.org/10.1007/s00572-015-0669-9 PMID: 26585898

20. Sun C, Yang M, Li Y, Tian J, Zhang Y, Liang L, et al. Comprehensive analysis of variation of cadmium accumulation in rice and detection of a new weak allele of OsHMA3. J Exp Bot. 2019; 70: 6389–6400. https://doi.org/10.1093/jxb/erz400 PMID: 31494666

21. Liu S, Zou W, Lu X, Bian J, He H, Chen J, et al. Genome-wide association study using a Multiparent Advanced Generation Intercross (MAGIC) population identified QTLs and candidate genes to predict shoot and grain zinc contents in rice. Agriculture. 2021; 11: 70. https://doi.org/10.3390/agriculture11010070

22. Swamy BPM, Rahman MA, Inabangan-Asilo MA, Amparado A, Manito C, Chadha-Mohanty P, et al. Advances in breeding for high grain zinc in rice. Rice. 2016; 9: 49. https://doi.org/10.1186/s12284-016-0122-5 PMID: 27671163

23. Maganti S, Swaminathan R, Parida A. Variation in iron and zinc content in traditional rice genotypes. Agric Res. 2020; 9: 316–328. https://doi.org/10.1007/s40003-019-00429-3

24. Pinson SRM, Tarpley L, Yan W, Yeater K, Lahner B, Yakubova E, et al. Worldwide genetic diversity for mineral element concentrations in rice grain. Crop Science. 2015; 55: 294–311. https://doi.org/10.2135/cropsci2013.10.0656

25. Eizenga GC, Ali MdL, Bryant RJ, Yeater KM, McClung AM, McCouch SR. Registration of the Rice Diversity Panel 1 for genome-wide association studies. Journal of Plant Registrations. 2014; 8: 109–116. https://doi.org/10.3198/jpr2013.03.0013crmp

26. Norton GJ, Douglas A, Lahner B, Yakubova E, Guerinot ML, Pinson SRM, et al. Genome wide association mapping of grain arsenic, copper, molybdenum and zinc in rice (Oryza sativa L.) grown at four international field sites. PLOS ONE. 2014; 9: e89685. https://doi.org/10.1371/journal.pone.0089685 PMID: 24586963

27. Yang M, Lu K, Zhao F-J, Xie W, Ramakrishna P, Wang G, et al. Genome-wide association studies reveal the genetic basis of ionomic variation in rice. The Plant Cell. 2018; 30: 2720–2740. https://doi.org/10.1105/tpc.18.00375 PMID: 30373760

28. Cu ST, Warnock NI, Pasuquin J, Dingkuhn M, Stangoulis J. A high-resolution genome-wide association study of the grain ionome and agronomic traits in rice Oryza sativa subsp. indica. Sci Rep. 2021; 11: 19230. https://doi.org/10.1038/s41598-021-98573-w PMID: 34584121

29. Bollinedi H, Yadav AK, Vinod KK, Gopala Krishnan S, Bhowmick PK, Nagarajan M, et al. Genome-wide association study reveals novel marker-trait associations (MTAs) governing the localization of Fe and Zn in the rice grain. Front Genet. 2020; 11: 213. https://doi.org/10.3389/fgene.2020.00213 PMID: 32391041

30. Pham CH, Do TD, Nguyen HTL, Hoang NT, Tran TD, Vu MTT, et al. Genome-wide association mapping of genes for anthocyanin and flavonoid contents in Vietnamese landraces of black rice. Euphytica. 2023; 220: 11. https://doi.org/10.1007/s10681-023-03268-0

31. Phung NTP, Mai CD, Mournet P, Frouin J, Droc G, Ta NK, et al. Characterization of a panel of Vietnamese rice varieties using DArT and SNP markers for association mapping purposes. BMC Plant Biol. 2014; 14: 371. https://doi.org/10.1186/s12870-014-0371-7 PMID: 25524444

32. Phung NTP, Mai CD, Hoang GT, Truong HTM, Lavarenne J, Gonin M, et al. Genome-wide association mapping for root traits in a panel of rice accessions from Vietnam. BMC Plant Biology. 2016; 16: 64. https://doi.org/10.1186/s12870-016-0747-y PMID: 26964867

33. Ta KN, Khong NG, Ha LT, Nguyen DT, Mai CD, Hoang GT, et al. A genome-wide association study using a Vietnamese landrace panel of rice (Oryza sativa) reveals new QTLs controlling panicle
morphological traits. BMC Plant Biology. 2018; 18: 282. https://doi.org/10.1186/s12870-018-1504-1
PMID: 30428844

34. Hoang GT, Gantet P, Nguyen KH, Phung NTP, Ha LT, Nguyen TT, et al. Genome-wide association mapping of leaf mass traits in a Vietnamese rice landrace panel. PLoS One. 2019; 14: e0219274. https://doi.org/10.1371/journal.pone.0219274 PMID: 31283792

35. Hoang GT, Lam VD, Nguyen TT, Ta NK, Gathignol F, Mai CD, et al. Genome-wide association study of a panel of Vietnamese rice landraces reveals new QTLs for tolerance to water deficit during the vegetative phase. Rice. 2019; 12: 4. https://doi.org/10.1186/s12284-018-0258-6 PMID: 30701393

36. To HTM, Nguyen HT, Dang NTM, Nguyen NH, Bui TX, Lavarenne J, et al. Unraveling the genetic elements involved in shoot and root growth regulation by jasmonate in rice using a genome-wide association study. Rice. 2019; 12: 69. https://doi.org/10.1186/s12284-019-0327-5 PMID: 31485824

37. Le TD, Gathignol F, Vu HT, Nguyen KL, Tran LH, Vu HTT, et al. Genome-wide association mapping of salinity tolerance at the seedling stage in a panel of Vietnamese landraces reveals new valuable QTLs for salinity stress tolerance breeding in rice. Plants (Basel). 2021; 10: 1088. https://doi.org/10.3390/plants10061088 PMID: 34071570

38. To HTM, Le KQ, Van Nguyen H, Duong LV, Kieu HT, Chu QAT, et al. A genome-wide association study reveals the quantitative trait locus and candidate genes that regulate phosphate efficiency in a Vietnamese rice collection. Physiol Mol Biol Plants. 2020; 26: 2267-2281. https://doi.org/10.1007/s12298-020-00902-2 PMID: 33268928

39. Danku JMC, Lahner B, Yakubova E, Salt DE. Large-scale plant ionomics. In: Maathuis FJM, editor. Plant Mineral Nutrients: Methods and Protocols. Totowa, NJ: Humana Press; 2013. pp. 255-276. https://doi.org/10.1007/978-1-62703-152-3_17 PMID: 23073889

40. van Rossum B. statgenSTA: Single Trial Analysis (STA) of field trials. 2024. Available from: https://github.com/biometris/statgensta.

41. Rodríguez-Álvarez MX, Boer MP, van Eeuwijk FA, Eilers PHC. Correcting for spatial heterogeneity in plant breeding experiments with P-splines. Spatial Statistics. 2018; 23: 52-71. https://doi.org/10.1016/j.spasta.2017.10.003

42. Wang J, Zhang Z. GAPIT Version 3: Boosting power and accuracy for genomic association and prediction. Genomics Proteomics Bioinformatics. 2021; 19: 629-640. https://doi.org/10.1016/j.gpb.2021.08.005 PMID: 34492338

43. Shin J-H, Blay S, McNeney B, Graham J. LDheatmap: An R function for graphical display of pairwise linkage disequilibria between single nucleotide polymorphisms. Journal of Statistical Software. 2006; 16: 1-9. https://doi.org/10.18637/jss.v016.c03

44. Sato Y, Takehisa H, Kamatsuki K, Minami H, Namiki N, Ikawa H, et al. RiceXPro Version 3.0: expanding the informatics resource for rice transcriptome. Nucleic Acids Research. 2013; 41: D1206-D1213. https://doi.org/10.1093/nar/gks1125 PMID: 23180765

45. Pinson SRM, Heuschele DJ, Edwards JD, Jackson AK, Sharma S, Barnaby JY. Relationships among arsenic-related traits, including rice grain arsenic concentration and straighthead resistance, as revealed by genome-wide association. Frontiers in Genetics. 2022; 12: 787767. https://doi.org/10.3389/fgene.2021.787767 PMID: 35371188

46. Zhang G-M, Zheng T-Q, Chen Z, Wang Y-L, Wang Y, Shi Y-M, et al. Joint exploration of favorable haplotypes for mineral concentrations in milled grains of rice (*Oryza sativa* L.). Frontiers in Plant Science. 2018; 9: 447. https://doi.org/10.3389/fpls.2018.00447 PMID: 29706977

47. Kumar S, Seem K, Kumar S, Vinod KK, Chinnusamy V, Mohapatra T. Pup1 QTL regulates gene expression through epigenetic modification of DNA under phosphate starvation stress in rice. Frontiers in Plant Science. 2022; 13: 871890. https://doi.org/10.3389/fpls.2022.871890 PMID: 35712593

48. Lu S, Ye J, Li H, He F, Qi Y, Wang T, et al. The splicing factor OsSCL26 regulates phosphorus homeostasis in rice. Plants. 2023; 12: 2326. https://doi.org/10.3390/plants12122326 PMID: 37375951

49. Hinrichs M, Fleck AT, Biedermann E, Ngo NS, Schreiber L, Schenk MK. An ABC transporter is involved in the silicon-induced formation of casparian bands in the exodermis of rice. Front Plant Sci. 2017; 8: 671. https://doi.org/10.3389/fpls.2017.00671 PMID: 28503184

50. Chauhan R, Awasthi S, Indoliya Y, Chauhan AS, Mishra S, Agrawal L, et al. Transcriptome and proteome analyses reveal selenium mediated amelioration of arsenic toxicity in rice (*Oryza sativa* L.). Journal of Hazardous Materials. 2020; 390: 122122. https://doi.org/10.1016/j.jhazmat.2020.122122 PMID: 32006842

51. Miro B, Longkumer T, Entila FD, Kohli A, Ismail AM. Rice seed germination underwater: Morpho-physiological responses and the bases of differential expression of alcoholic fermentation enzymes. Front Plant Sci. 2017; 8: 1857. https://doi.org/10.3389/fpls.2017.01857 PMID: 29123541
52. Hayashi S, Kuramata M, Abe T, Yamaguchi N, Takagi H, Tanikawa H, et al. Deficiency in alcohol dehydrogenase 2 reduces arsenic in rice grains by suppressing silicate transporters. Plant Physiol. 2021; 186: 611–623. https://doi.org/10.1093/plphys/kiab086 PMID: 33620496

53. Norton GJ, Deacon CM, Xiong L, Huang S, Meharg AA, Price AH. Genetic mapping of the rice ionome in leaves and grain: Identification of QTLs for 17 elements including arsenic, cadmium, iron and selenium. Plant Soil. 2010; 329: 139–153. https://doi.org/10.1007/s11104-009-0141-8

54. Huang D, Fan Y, Hu B, Xiao Y, Chen D, Zhuang J. Assessment and genetic analysis of heavy metal content in rice grain using an *Oryza sativa* × *O. rufipogon* backcross inbred line population. Journal of the Science of Food and Agriculture. 2018; 98: 1339–1345. https://doi.org/10.1002/jsfa.8598 PMID: 28758676

55. Chen X, Chen F, Sun S, Li Y, Li Y, Mo H, et al. Effect of polishing on lead and cadmium bioavailability in rice and its health implications. Foods. 2022; 11: 2718. https://doi.org/10.3390/foods11172718 PMID: 36076903

56. Bouis H, Hotz C, McClafferty B, Meenakshi J, Pfeiffer W. Biofortification: A new tool to reduce micronutrient malnutrition. Food and nutrition bulletin. 2011; 32: S31–40. https://doi.org/10.1177/15648265110321S105 PMID: 21717916

57. Mahender A, Anandan A, Pradhan SK, Pandit E. Rice grain nutritional traits and their enhancement using relevant genes and QTLs through advanced approaches. SpringerPlus. 2016; 5: 2086. https://doi.org/10.1186/s40064-016-3744-6 PMID: 28018794

58. Yang L, Wang Y. Impact of climate change on rice grain quality. In: Bao J, editor. Rice (Fourth Edition). AACC International Press; 2019. pp. 427–441. https://doi.org/10.1016/B978-0-12-811508-4.00013-7

59. Drobner U, Tyler G. Conditions controlling relative uptake of potassium and rubidium by plants from soils. Plant and Soil. 1998; 201: 285–293. https://doi.org/10.1023/A:1004319803952

60. De Datta SK. Principles and practices of rice production. John Wiley and Sons. 1981.

61. Baxter I. Ionomics: Studying the social network of mineral nutrients. Curr Opin Plant Biol. 2009; 12: 381–386. https://doi.org/10.1016/j.pbi.2009.05.002 PMID: 19481970

62. Baxter I. Ionomics: The functional genomics of elements. Briefings in Functional Genomics. 2010; 9: 149–156. https://doi.org/10.1093/bfgp/elp055 PMID: 20081216

63. Bryant RJ, Dorsch JA, Peterson KL, Rutger JN, Raboy V. Phosphorus and mineral concentrations in whole grain and milled low phytic acid (Ipa) 1–1 rice. Cereal Chemistry. 2005; 82: 517–522. https://doi.org/10.1094/CC-82-0517

64. Zhang M, Pinson SRM, Tarpley L, Huang X-Y, Lahner B, Yakubova E, et al. Mapping and validation of quantitative trait loci associated with concentrations of 16 elements in unmilled rice grain. Theor Appl Genet. 2014; 127: 137–165. https://doi.org/10.1007/s00122-013-2207-5 PMID: 24231918

65. Williams PN, Islam S, Islam R, Jahruddin M, Adomako E, Soliaman ARM, et al. Arsenic limits trace mineral nutrition (selenium, zinc, and nickel) in Bangladesh rice grain. Environ Sci Technol. 2009; 43: 8430–8436. https://doi.org/10.1021/es901825t PMID: 19924980

66. Imran M, Hussain S, Rana MS, Saleem MH, Rasul F, Ali KH, et al. Molybdenum improves 2-acetyl-1-pyrroline, grain quality traits and yield attributes in fragrant rice through efficient nitrogen assimilation under cadmium toxicity. Ecotoxicology and Environmental Safety. 2021; 211: 111911. https://doi.org/10.1016/j.ecoenv.2021.111911 PMID: 33453638

67. Han Z, Wei X, Wan D, He W, Wang X, Xiong Y. Effect of molybdenum on plant physiology and cadmium uptake and translocation in rape (Brassica napus L.) under different levels of cadmium stress. Int J Environ Res Public Health. 2020; 17: 2355. https://doi.org/10.3390/ijerph17072355 PMID: 32244320

68. Huang L, Li WC, Tam NFY, Ye Z. Effects of root morphology and anatomy on cadmium uptake and translocation in rice (*Oryza sativa* L.). Journal of Environmental Sciences. 2019; 75: 296–306. https://doi.org/10.1016/j.jes.2018.04.005 PMID: 30473295

69. Mendel RR, Bittner F. Cell biology of molybdenum. Biochimica et Biophysica Acta (BBA)—Molecular Cell Research. 2006; 1763: 621–635. https://doi.org/10.1016/j.bbamcr.2006.03.013 PMID: 16784786

70. Ma JF, Tamai K, Yamaji N, Mitani N, Konishi S, Katsuhara M, et al. A silicon transporter in rice. Nature. 2006; 440: 688–691. https://doi.org/10.1038/nature04590 PMID: 16572174

71. Becker M, Ngo NS, Schenk MKA. Silicon reduces the iron uptake in rice and induces iron homeostasis related genes. Sci Rep. 2020; 10: 5079. https://doi.org/10.1038/s41598-020-61718-4 PMID: 32193423