Global Geographical Variation in Elemental and Arsenic Species Concentration in Paddy Rice Grain Identifies a Close Association of Essential Elements Copper, Selenium and Molybdenum with Cadmium

Despite the centrality of staple grains for human well-being, both as a source of nutrients and of toxic ions, there is little understanding of where and how elements vary, and if there are particular elements that correlate. Here, for shop bought polished (white) rice, we comprehensively characterized trace (arsenic species, cadmium, copper, iron, manganese, molybdenum, rubidium and zinc) and macro-nutrients (calcium, chlorine, potassium, phosphorus and sulphur) for grain purchased in 18 countries, across four continents, a total of 1045 samples. This was to investigate if there were any major differences between geographic location and elemental content, and to observe if there were any patterns in elemental distribution. Greatest variation in the median was observed for the non-essential rubidium (15-fold) and arsenic species (fivefold). Rubidium was the highest in the Americas, lowest in Europe, while inorganic arsenic (iAs) and dimethylarsonic acid (DMA) were low for Africa and high in the South American and European continents. The highest concentrations of cadmium were found in Asian samples, and lowest in South America, with variation within these regions. At the extremes of individual counties, China had fivefold higher concentrations than the global median, while Tanzania was fourfold lower than this value. Calcium, potassium, molybdenum and phosphorus were the highest in European and lowest in African grain, though the fold-differences were relatively low, ~ 0.2, while iron was the highest in African grain and lowest in European, Asian and South American grain, with a ~ twofold difference. Selenium was also higher in Africa versus other regions, and copper, manganese and zinc were the highest in American grain. Factor analysis showed that copper, cadmium, molybdenum, rubidium and selenium were strongly associated together, and these element’s factor loadings were diametrically opposed to less tightly associated calcium, chlorine, manganese, potassium, phosphorus and sulphur. Stepwise additions linear region analysis was performed on log-transformed concentrations to investigate cadmium associations in more detail. Selenium was the greatest predictor of cadmium concentration, followed by molybdenum, accounting for over 50% of the contribution to the adjusted R2. Arsenic species were only weakly correlated with other elements. The implications for these findings with respect to dietary nutrition are discussed. Vietnamese rice was notable in being deficient in macro- and micro-nutrients while also being elevated in cadmium at a median of 0.02 mg/kg, with China though still having a median that is ~ 2.5-fold this concentration. These Chinese concentrations are of particular concern as the 75th percentile for China is 0.1 mg/kg, a value that triggers regulatory action for rice products.


Introduction
Hidden hunger is one of the greatest issues facing society globally, defined as where food delivers adequate calories but is deficient in essential nutrients (Gashu et al. 2021). This is nowhere so prominent than in subsistence diets that rely on polished grains, where milling grains, unfortunately, removes a substantial portion of nutrients present in wholegrain (Rahman et al. 2019). Grain mineral nutrition varies due to soil conditions and cultivation and fertilization practice (Nyachoti et al. 2021;Gashu et al. 2021;Welsh and Graham 2004;White and Broadley 2005), and to crop genetics (Cobb et al. 2021;Huang et al. 2015;Norton et al. 2014). Maximizing grain nutrient density, to deliver as great a quantity of each essential element as feasible, given that bioavailability in the gut is also an issue for elements such as iron (Hoekenga et al. 2011), is a major goal for those wishing to improve human health (Gashu et al. 2021;Rahman et al. 2019).
Food security, supplying a population sufficient quantities of nutritionally balanced foods, has to be combined with food safety, minimizing exposure to toxic compounds. White rice is particularly problematic in both regards as it is not overly nutrient dense, while at the same time being widely elevated, ~ tenfold compared to other grains, in cadmium and arsenic (Rahman et al. 2019). There is evidence that for certain elements such as arsenic (Williams et al. 2007), cadmium (Nan et al. 2002) and selenium ) that rice grain concentration is linked to the concentration in the soils in which they are grown. These elements tend to act as analogues of essential nutrients, phosphate and silicic acid for arsenic, dependant on whether arsenic is speciated as arsenate or arsenite, respectively (Meharg and Meharg 2015), nutrient divalent metal cations for cadmium (Schaaf et al. 2004), and sulphur for selenium ). The same is the case for nonessential rubidium, an analogue of essential potassium (Oliferuk et al. 2017). Most nutrient elements are thought to be in homeostasis with plant tissues, given that plant cells have quite stringent requirements for optimum functioning, and all elements/ions can be considered as stressors if outside those optimum activities in cells (Clemens 2001;Ernst et al. 2008;Meharg 1993). Some nutrients can be stored in excess through cellular sub-localization, or the formation of compounds/complexes suitable for storage. Phosphorus, an essential but generally poorly available in soils, undergoes both these storage strategies under luxury uptake, though plants rapidly downregulate phosphate assimilation if leaf tissue becomes too high (Chiou et al. 2006). Quantitative trait loci analysis, also indicates for some elements, notably cadmium, lead, selenium, zinc, that grain accumulation is under genetic control in rice (Cobb et al. 2021;Huang et al. 2015;Norton et al. 2014). Given these complexities, for food security and safety mitigation strategies, it is important to know how elemental concentrations vary regionally for foods, Deficiency or toxicity mitigation for a specific region may include soil fertilization, water management, food supplementation or high grain nutrient directed breed strategies, i.e. food fortification of biofortification (Welsh and Graham 2004;White and Broadley 2005). Global and regional assessments of grain mineral nutrition tend to focus on food fraud, using elemental analysis to try and identify geographical origin (Ariyama et al. 2012;Chung et al. 2018); or geo-nutrition, relating crop grain element concentrations to region of origin (Adomako et al. 2011;Gashu et al. 2021), with a number of studies also focusing on dietary consequences of grain elemental content (Gashu et al. 2021;Joy et al. 2017;Nyachoti et al. 2021). These studies generally use inductively coupled plasma-mass spectrometry (ICP-MS), and this has led to a better characterization of micro-elemental quantification, to which ICP-MS is suited, as compared to macro-elements that are generally too elevated in grains for adequate ICP-MS quantification, with the exception of phosphorus. Of particular interest with respect to ICP-MS quantification are arsenic and cadmium which are of concern in rice but too low to be measured by other techniques such as X-ray fluorescence (XRF) (Rahman et al. 2019). Furthermore, arsenic in rice grain is speciated as iAs (the sum of arsenate and arsenite) and DMA, and it is iAs that is regulated not DMA, yet DMA can constitute up to ~ 50% of grain speciation in some regions, and chromatography interfaced to ICP-MS is required to speciate grain arsenic (Rahman et al. 2019). XRF is better suited to macro-element analysis and is becoming more established in characterizing foods, given that issues with external calibration of XRF need to be overcome by using an appropriate range of Certified Reference Materials (CRMs) (Oliveira et al. 2020;Rahman et al. 2019). Also, for arsenic in particular, there is a need to speciate as toxicity and regulation are driven by iAs content, yet rice grain can have highly variable concentrations of dimethylarsonic acid (DMA) (Adomako et al. 2011;Carey et al. 2020). Here, we report macro-and micro-elemental, along with arsenic species, concentrations of market white rice, chosen to reflect direct dietary consumption of 1045 samples, collected across 18 countries from four major rice growing continents; Africa, South America, Asia and Europe. Rice was analysed by ICP-MS, ion chromatography (IC)-ICP-MS (to ascertain arsenic species concentrations) and XRF. Data were interpreted to examine commonalities in elemental and species patterns, and to identify specific deficiency/toxicity issues for given regions/countries.

Materials and Methods
The sampling strategies are outlined in the publications Carey et al. (2020) and Shi et al. (2020). In brief, polished (white) rice was obtained from 18 countries, and 4 continents, from retail outlets so as to best characterize consumer exposure. In total 1045 samples are reported here, a subset from the 1180 reported in Carey et al. (2020) on arsenic speciation, and 2270 from Shi et al. (2020) on cadmium. The current dataset only reports samples where additional XRF characterization has been conducted, leading to reduced sample coverage, in both number and in countries assessed.
As noted in Carey et al. (2020), the type of market basket survey conducted here has its limitations due to the fact that multiple people conducted the survey, with multiple sampling objectives, with each country sampled having differing geographical size, and the surveyors having differing geographical reaches within a given country. Furthermore, some countries only have one rice season per annum, others have multiple, and there is often within country variation on number of seasons based on geographical setting. Also, genetic variation in rice occurs at the subspecies level (Aus, Japonica, Indica and hybrid etc.), and there is considerable cultivar variation within these subspecies. Markets may directly source rice from neighbouring farmers, or at the other extreme, source rice from larger mills that can operate on a regional scale. The nature of variation in rice sourcing will lead to heterogeneity in retail supply. The data presented here does not include retail outlet type, genetic information, seasonal information, etc., and this must be born in mind in interpreting this study. For a global sampling campaign to include such detail, and the subsequent increase in sample size to account for these multiple variables, were well beyond the resources of this current investigation. However, the utility of the sampling approach used here can be validated by comparison with local studies where details of interest are recorded. Examples of more detailed studies cross-referenced to the current investigation in the discussion are country-specific surveys of grain (Adomako et al. 2011;Joy et al. 2017;Williams et al. 2009); and the role of genetic variation (Huang et al. 2015;Norton et al. 2010Norton et al. , 2012Norton et al. , 2014, with respect to rice grain mineral nutrition. Methodologies for the well-established rice inorganic analysis that we use here are have been published already. A concise overview of the analysis is now presented. Rice samples were freeze dried (Christ Alpha1-4LD) and milled (PM100 Ball Mill, Retsch, Germany) and subsampled for ICP-MS (Thermo Scientific iCap Q), both for arsenic speciation (interfaced with in chromatography, IC), and elemental totals, and for XRF (Nex CG, Rigaku, Japan). ICP-MS and IC-ICP-MS protocols are described in detail in Carey et al. (2020) and Shi et al. (2020), and XRF protocols in Rahman et al. (2019). Certified Reference Material (CRM) rice flour, NIST 1568b, was included in all batches of samples for both elemental totals (ICP-MS and XRF) and for arsenic speciation, and it is certified for both a range of elemental totals and for arsenic speciation. Total elemental analysis was conducted by weighing ~ 100 mg of sample, accurately, into 50 ml polypropylene (PP) tubes, and then adding 2 ml of conc. Aristar nitric acid. The PP tubes were then left to stand for 24 h. Before initiating digestion, 2 ml of undiluted hydrogen peroxide (reagent grade) was pipetted in, then microwaved digested (Mars CEM), using programmes detailed in Carey et al. (2020). Samples were cooled to room temperature on completion of microwaving, and diluted to 30 ml with double distilled deionized water. The protocol for arsenic speciation was identical except that except that 10 ml of one percent nitric acid was used for as the extractant instead of concentrated nitric acid. Prior to IC-ICP-MS speciation, 7 µl of hydrogen peroxide added to 700 µl of the extraction, to convert arsenite to arsenate for subsequent IC chromatography, as arsenite elutes at the solvent front.
ICP-MS analysis was calibrated using eight standards, including one blank. Multi-Element 2 (SPEX CLMS-2 Multi-Element Solution 2, matrix: five percent nitric acid) and Multi-Element 4 (SPEX CLMS-4 Multi-Element Solution 4, matrix: water/Tr-HF) were used to make up all standards in a range of 0-100 µg/l in one percent nitric acid. Rhodium internal standard (Fluka Analytical rhodium internal ICP-MS standard) was added to all standards and all samples also at ten µg/l. For IC-ICP-MS, authentic samples of arsenate, monomethylarsonic acid (MMA) and DMA were used to determine retention times, and the calibration was conducted on DMA. XRF (Rigaku NEX CG) was performed on the ball-milled samples to determine elements that are either problematic on non-determinable by ICP-MS (see Table 1).
Only elements with certified quantification for rice flour CRM NIST 1568b were reported. Reported elements were within the range of 67-109% recovery, with CRM recoveries and Limits of Detection (LoD) for ICP-MS techniques, and Detection Limits (DL), for XRF as this method has an internal calibration, are given in Table 1. As a range of elements can be quantified by both ICP-MS and XRF, and for these elements, the technique with the best recovery of NIST rice flour CRM NIST 1568b was selected for reporting, with the technique used for reported analysis also listed in Table 1.
Statistical analysis considered the data both as grouped by continent or as grouped by country. As the datasets are not normally distributed, Kruskal-Wallis test used to compare medians and stepwise linear regression to look at relationships between elements was conducted in SPSS v.27, and principal component analysis (PCA) and Spearman's rank correlation were conducted using GraphPad Prism v.9.

Results
Medians, 25th and 75th percentiles, are presented for elements and species by continent in Fig. 1. Kruskal-Wallis analysis was used to investigate differences in medians between the four regions. For all elements and species presented here, all were significantly different between continents, all at P < 0.001, with the exception of Mn (P = 0.002) and Zn (P = 0.032). Greatest variation in the median was observed for the non-essential rubidium, being the highest in the Americas, lowest in Europe, with a 15-fold difference. Arsenic species, iAs and DMA, were the next most variable in their medians with both species being low in Africa and high in the South American and European continents, with a ~ fivefold difference, while iAs was relatively high in Asia, but DMA low. Cadmium had a different profile, being the highest in Asia and lowest in South America. A range of macro-nutrient elements, calcium, potassium, molybdenum and phosphorus, were the highest in European, and lowest in African grain, though the fold-differences were relatively low, ~ 0.2. Conversely, iron was high in African grain and lowest in European, Asian and South American grain, ~ twofold difference. Chloride was also high in African and Asian grain, as compared to American and European. Selenium was also higher in Africa versus other regions. Copper, manganese and zinc were the highest in American grain, though there was relatively good homeostasis in concentrations across the continents. Sulphur was also found to have a similar homeostasis.
Presentation of medians for individual countries are given in Table 2, and country medians as a percentage of overall mean of country medians in Fig. 2. Kruskal-Wallis analysis found that all counties differed in their median, all at P < 0.0001 with the exception of zinc where P = 0.046. Tanzania is particularly low in iAs DMA, cadmium and rubidium, but high in chloride. Mali is low in selenium, along with Argentina, Paraguay, Uruguay South Korea, Indonesia and Spain. There is a fourfold difference in selenium concentrations for countries across the globe. For some countries there are notable high or low concentrations for particular elements. South Korea is high in calcium compared to all other countries, while Vietnam is particularly low in this element. Cadmium is specifically highly elevated in China, while other Asian countries also tend to be high in this element, ~ 500% higher than the global mean. Chloride is elevated in Chile, phosphorus in Sri Lanka, and sulphur in Bangladesh. Bangladesh is specifically low in zinc, with a median 70% of that of the overall country medians, but high in iron at 200% of the median (Fig. 2). Countries with notably low macro-and micro-nutrient elements are South Korea for copper, Chile and Sri Lanka for sulphur, zinc in Uruguay and Argentina, and particularly low zinc for Bangladesh, with almost half the concentration of this element compared to other countries. Vietnam is notable as it is deficient in calcium, iron, potassium, manganese and molybdenum, and no other country matches this broad pattern of macro-and micronutrient deficiencies (Fig. 2). For Africa, a depletion in arsenic species is apparent, as is elevation in selenium (Tanzania and Ivory Coast) and iron (universally) (Fig. 2).
There is a specific elevation in Ivory Coast grain cadmium relative to other African countries (Fig. 2). What distinguishes American countries, with the exception of Bolivia, from other regions of the world are that they vary from moderate to highly elevated DMA, ~ 300%, and are relatively depleted in selenium especially in Argentina, Paraguay and Uruguay. Asian countries tend to be elevated in cadmium, but interestingly not South Korea. European   rice was high in DMA and low in rubidium, and Italy was relatively high in cadmium. Turkish rice was notably low, ~ 50%, in selenium. PCA and associated factor analysis for the entire data set are presented in Fig. 3. The PCA score means by country are given with the individual scores due to the high data density. The factor loadings analysis showed that copper, cadmium, molybdenum, rubidium and selenium were strongly associated together, and these element's factor loadings were diametrically opposed to less tightly associated calcium, chlorine, manganese, potassium phosphorus and sulphur, with similar strengths of loadings. iAs and zinc were also associated, but less strongly, with this latter cluster. Iron and DMA had factor loadings close to the origin, suggesting that there were not closely associated with any other elements. Element/species pairs that might be thought to be associated with other are diametrically opposed such as copper and zinc (elevated concentrations of both derived from the same pollutant sources); iAs and cadmium (again being derived in elevated concentrations in soil through pollutant  Fig. 2 Country medians expressed as a percentage relative to the mean of all countries for each element, plotted by continent sources, such as atmospheric deposition of contamination caused by base and precious mining and smelting activities, fossil fuel burning etc.); and the physiological analogues sulphur and selenium. However, potassium and phosphorus factor loadings were similar in direction and strength, which may relate to the fact that these two elements are commonly fertilized together. The American countries, with the exception of Chile, were associated with the copper, cadmium, molybdenum, rubidium and selenium factor loadings, but the countries of other continents were dispersed around the plot, though the nearby countries of Mali and Ivory Coast differed considerably in cadmium content. Vietnam, in particular, has a mean PCA score to the right of the graph but near to the X-axis. This shows, again, that Vietnamese grains are generally nutrient deficient, as shown by plots of the country value of each element expressed as a percentage of the median of country medians given in Table 2 and percentages of medians in Fig. 3. To validate/cross-reference the PCA factor loading observations for cadmium, log-transformed elemental concentrations underwent stepwise additions linear regression  Table 3 Stepwise regression analysis of log cadmium grain concentrations verses the log transformed other elements quantified Section A shows the R 2 and the significance of each model, while B shows the slopes (Beta) and significance for the model containing all elements in the final model with log cadmium concentration as the dependent variable (Table 3). Selenium was the greatest predictor of cadmium, followed by molybdenum, accounting for over 50% of the contribution to the adjusted R 2 . Log iAs in stepwise linear regression was highly correlated with log DMA with an adjusted R 2 of 0.556, but not with other elements. When log DMA was removed as a predictor, the next highest correlating parameter was log molybdenum with an adjusted R 2 of 0.029 (data not tabulated). Figure 4 shows the correlation matrix for the log transformed concentration dataset to enable further comparisons of all elements/species against each other. Besides the strong positive correlation between iAs and DMA, there was also a very strong relationship between potassium and phosphorus (Spearman's R of 0.799). This strong association between potassium and phosphorus is also seen for the factor analysis (Fig. 3). Also notable from the Spearman's rank correlation is that chlorine is negatively correlated with a suite of elements/species including arsenic species, cadmium, copper, molybdenum, phosphorus and rubidium.

Discussion
Plant nutrition is dependent on soil availability and concentration, and plant physiological regulation as determined by genetics (Cobb et al. 2021;Norton et al. 2014). It is considered that most essential elements and ions are in homeostasis for plant tissues in general, and in grain specifically, as outside an individual ion/element's concentration range, the plant will either be in deficiency or toxicity (Bohnert et al. 2006). One essential element that does not appear under tight physiological control, with grain concentration being closely related to soil concentrations is selenium ). This is thought to be because selenium species are analogues of the much more abundant sulphur, and with sulphur ameliorating selenium species toxicity ). Non-essential elements/ions also enter grain as analogues of essential nutrients, such as arsenate, a phosphate analogue (Meharg 1993); arsenite, a silicic acid analogue (Meharg and Meharg 2015); rubidium, a potassium analogue (Oliferuk et al. 2017); and cadmium an analogue of a range of divalent cations (Nan et al. 2002). Thus, a key finding here is that the essential elements copper, molybdenum and selenium are tightly associated with non-essential cadmium and rubidium. Hypotheses regarding the relationship between cadmium with copper, molybdenum, rubidium and selenium are now considered.
Copper and selenium can be elevated in soils due to anthropogenic processes, such as mineral smelting and processing and coal burning, where selenium soil pollution actually supplements plants as it is deficient in most soils ). As noted, selenium can be very variable in plant tissues, while copper is under tight plant homeostasis due its essentiality in a range of cellular functions and due to its toxicity (Clemens 2001). Contrastingly to copper, a divalent metal ion that is also an environmental contaminant, divalent cadmium is also a widely dispersed pollutant that can reach concentrations in grain that cause human health concerns in a wide range of agricultural settings, yet is highly variable in grain (Hu et al. 2016;Nan et al. 2002;Shi et al. 2020). Rubidium differs yet again, being a monovalent ion, and is thought to be preferentially mobilized into plants under acidic soil conditions where soil available potassium is low (Drobner and Tyler 1998). Molybdenum chemistry in soils and plants is complex, with molybdate thought to be a sulphate analogue, and is assimilated by plants more readily at higher pHs due to reduced binding adsorption to amorphous oxides, while sulphur addition also decreased molybdenum assimilation (McGrath et al. 2010). In the context of the close assimilation of the group of elements found to cluster here (cadmium with copper, molybdenum, rubidium and selenium), cadmium is most studied (Hu et al. 2015(Hu et al. , 2016Shi et al. 2020). Cadmium is assimilated under more aerobic conditions of rice cultivation (Hu et al. 2015) and this may be a clue as to why other elements are associated with cadmium, i.e. it could be hypothesized that more aerobic conditions are needed to mobilize these elements in paddy soils effectively. This supposition remains to be tested.
Cadmium, in both soils and plants, is readily complexed by sulphur (Ernst et al. 2008). High grain cadmium is diametrically opposed to high grain sulphur for factor loading plots, and thus when sulphur is high in plants it can be postulated that translocation of cadmium to grain is low Fig. 4 Spearman's rank correlation matrix of log-transformed elements for rice grain concentrations. The scale to the right of the graph is the R-value for analysis as cadmium complexation with sulphur groups of phytochelatins (PCs) are sequestered in the vacuole, retarding their translocation to the plant (Ernst et al. 2008). It is sulphide immobilization under reduced redox conditions that causes cadmium to have poor soil availability under flooded paddy conditions (de Livera et al. 2011;Rinklebe et al. 2016), again emphasizing the intimate link between cadmium and sulphur highlighted in this current study. This sulphur association may also explain grain copper concentrations, with copper like cadmium having a high affinity for sulphur (Ernst et al. 2008;Hartley-Whitaker et al. 2001). Notably, selenium is a sulphur analogue ). For selenium, the ratio of selenium to sulphur increases, maybe as the plant tries to mine plant sulphur, inadvertently assimilating more selenium, i.e. that selenium may be increased in grain under sulphur deficiency. The same could be argued for rubidium and potassium, i.e. that rubidium is enhanced in grain when potassium is deficient, as rubidium is orthogonal to potassium with respect to its grain PCA factor loading plot. Molybdenum is assimilated through plant roots as molybdate (Baxter et al. 2008), an anion that is thought to interact strongly with sulphate (McGrath et al. 2010), i.e. it is postulated that molybdate is a sulphate analogue. Molybdenum is assimilated into grain when sulphur is low, according to PCA analysis, suggesting that under plant sulphur deficiency conditions that molybdenum is more readily assimilated and translocated. Interestingly, Williams et al. (2009) found that grain zinc and selenium was negatively correlated with grain arsenic for field sites across Bangladesh. Thus elemental toxicity, arsenic in this case, is impacting grain nutritional composition. For the current study, Bangladesh had the highest grain arsenic of this survey and the lowest zinc, supporting this finding, but the selenium grain concentration was median.
The global distribution of arsenic species in polished grain has been previously studied for this data-set . What is new here is the understanding that both iAs and DMA are not intimately linked to other elements with respect to co-location in grain. Arsenic species mobilization is strongly driven by soil redox, linked to the dissolution of iron plaque and to the microbial methylation of iAs (Meharg and Zhao 2012). It is notable that grain iron is also not strongly linked to other elements in this current investigation. It appears that other factors, such as potential pollution impact and weathering of soils, have a greater influence on grain arsenic species than co-element nutrition of plant or soils, given that arsenic species can also act as phosphate and silicic acid analogues in both soils and plants (Meharg and Zhao 2012). For example, it is the low pollutant impacted and heavy weathered soils of African and Asian southern hemispheres that are low in grain arsenic species in this current study.
With respect to genetic variation in grain elemental content, Norton et al. (2010) investigating 17 elements in grain through ICP-MS quantification found for an indica-japonica cross that there were 41 QTLs. QTLs for iron, molybdenum, selenium and zinc, in particular, were suggested to be suitably strong to be targets for improving grain concentrations through marker assisted selection. Huang et al. (2015) reported genetic variation in cadmium, iron, lead, selenium, and zinc, for wholegrain rice using Chinese association mapping accessions. Genome-wide association mapping found a strong genotype*environment interaction for loci associated with elements in whole grain rice grown at four locations across the globe, but variation was found for arsenic, copper, molybdenum and zinc (Norton et al. 2014), indicating weak genetic links between the elements. Genetic variation was found to be a large determining factor for grain arsenic concentrations across four field sites globally, Bangladesh, China, India and N. America (Norton et al. 2010). Association mapping has identified root and shoot transporters for cadmium and molybdenum, but these are on different loci and are for different putative transporters (Cobb et al. 2021). Thus, rice genetics, particularly between japonica and indica groups (Huang et al. 2015;Norton et al. 2010), could be an important source of variation in rice's global grain supply chain, but this needs to be systematically investigated in order to validate.
The data sets reported here enable rice from different regions to be grouped as deficient in particular nutrients.
The key country with respect to globally mineral nutrient deficiency in the grain was Vietnam as its grain was particularly deficient in a suite of macro-and micro-nutrients, while also being relatively elevated in cadmium. These deficiencies in grain elemental nutrients suggest that Vietnamese soils may be generally deficient in nutrients. This may be due to a range of factors such as over cropping as many Vietnamese paddies have multiple harvests per year (Ba et al. 2016). Soil acidification is widespread in S.E. Asia due to over use of mineral nitrogen fertilizers (Xie et al. 2019), and acidification will lead to issues with nutrient leaching and availability. Straw biomass is widely burnt in Vietnam leading to poor organic waste-stream recycling, including animal manures produced through straw consumption (Hung et al. 2019;Watanabe et al. 2019). The health consequences for the Vietnamese population for having rice low in nutrient density but high in toxicants needs assessing, given that it is a subsistence rice economy with high rates of per capital consumption (Hung et al. 2019;Watanabe et al. 2017). Another notable health concern is Bangladesh that has high iAs, high cadmium and low zinc, given that zinc is vital in enzymatic roles that are key to the general health of the population ). The low selenium content of rice is of concern for a range of South American countries, and for Spain, South Korea and Indonesia, because selenium deficiency is implicated in a wide range of disease outcomes at a population level ).
The research findings of this current study provide the basis to start to redress dietary deficiencies in populations that have rice as a subsistence diet. This could be through changing agronomic practices or ensuring key nutrients reach the diet through other food groups. The underlying rice grain nutrient deficiencies reported here require further epidemiological studies to investigate human consequences, and agronomic investigations with respect to how to mitigate against deficiencies through soil management and rice breeding.