Next Article in Journal
Estimation of the Enzymatic Activity of Haplic Chernozem under Contamination with Oxides and Nitrates of Ag, Bi, Te and Tl
Next Article in Special Issue
The Efficacy of Micronutrient Fertilizers on the Yield Formulation and Quality of Wheat Grains
Previous Article in Journal
Hyperspectral Vegetation Indices to Assess Water and Nitrogen Status of Sweet Maize Crop
Previous Article in Special Issue
Efficacy of Nitrogen and Zinc Application at Different Growth Stages on Yield, Grain Zinc, and Nitrogen Concentration in Rice
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Zinc Fertilization on Grain Cadmium Accumulation, Gene Expression, and Essential Mineral Partitioning in Rice

by
Michael Tavarez
1,2,
Michael A. Grusak
3 and
Renuka P. Sankaran
1,2,*
1
The Graduate Center, City University of New York, New York, NY 10016, USA
2
Department of Biological Sciences, Lehman College, City University of New York, Bronx, NY 10468, USA
3
USDA-ARS Edward T. Schafer Agricultural Research Center, Fargo, ND 58102, USA
*
Author to whom correspondence should be addressed.
Agronomy 2022, 12(9), 2182; https://doi.org/10.3390/agronomy12092182
Submission received: 29 July 2022 / Revised: 9 September 2022 / Accepted: 10 September 2022 / Published: 14 September 2022
(This article belongs to the Special Issue Micronutrient Homeostasis and Biofortification in Plants)

Abstract

:
Cadmium (Cd) is a toxic heavy metal that can cause severe health issues if ingested. Certain varieties of rice can accumulate high levels of the metal in edible tissues thereby transferring the toxin into the food chain. As chemical analogs, interactions between the essential mineral zinc and the toxic heavy metal cadmium play an important role in regulating the transport of both minerals to rice grains. Understanding these interactions is crucial for limiting cadmium and increasing zinc transfer to the food chain. Previous studies have reported conflicting results suggesting synergistic and antagonistic relationships between the minerals. The goal of this work was to identify the effect of external cadmium and zinc on the uptake and translocation of both minerals from roots to grains of rice that differ in grain cadmium concentrations. The results showed that a higher input of external zinc increased cadmium translocation and accumulation to the grain in two of three varieties, while external cadmium does not influence zinc accumulation. Cadmium synergy and antagonism with other essential minerals were also examined and the effects differed between rice lines. Our results showed that the differential expression of the transport proteins OsNramp5, OsHMA2, and OsHMA3 as well as genes involved in the synthesis of glutathione and phytochelatin could have contributed to differences in grain Cd accumulation. These results add to the knowledge of cadmium and zinc partitioning in one of the most consumed plant foods in the world and can assist fortification efforts to establish rice lines that are both safe and nutritious.

1. Introduction

Cadmium (Cd) is a toxic heavy metal whose industrial uses include batteries, plastic stabilizers, and pigments [1]. It is a contaminant that occurs naturally in some agricultural soils and threatens agricultural productivity, food safety and human health [2,3]. When available, plants can readily accumulate the contaminant in edible tissues and allow it to be transferred through the food chain. Oryza sativa has one of the highest propensities of translocating cadmium to the grains when compared to other cereal crops [4,5,6,7]. Dietary Cd is currently one of the leading causes of Cd poisoning and is linked to bone mineral density loss, renal disease, and cancers [8,9,10]. Since grains contribute 27% of human dietary consumption of cadmium, decreasing cadmium in grains will lower dietary exposure [11].
Essential micronutrients such as zinc (Zn) share a similar atomic structure with Cd, and therefore they can interact similarly with plant proteins. Metal specificity is a challenge in breeding low Cd cereals, because the alteration of pathways to decrease Cd transport could inadvertently decrease Zn transport. The uptake of Cd and Zn from soil and their distribution within the plant involve complex and dynamic processes that may be facilitated by common transporters [12]. As with micronutrients, Cd transport and accumulation in the grains is affected by several processes such as Cd uptake into the root, translocation to aboveground tissues, subcellular sequestration, and phloem transport [13,14,15,16,17].
Cadmium uses proteins from key transport families such as the natural resistance-associated macrophage proteins (NRAMP), the Zrt-Irt-like protein family (ZIP), Yellow stripe-like (YSL), and Heavy Metal p1b-type ATPase family (HMA), which are responsible for transporting essential minerals such as zinc, iron (Fe), and manganese (Mn) [4,15,16,18,19,20,21]. The transport pathway of Zn and Cd from soil begins at the root cells where the mineral enters the plant through either the symplastic or apoplastic pathways [22,23]. Although the exact mechanisms of Zn transport from soil to rice grain are not completely clear, several tightly regulated proteins from the ZIP, YSL, and HMA families have been identified [24,25,26,27,28]. The root uptake of Zn is mediated by OsZIP1 and OsZIP9, which are expressed in the epidermis of rice roots [29,30]. The transport of Zn to aboveground tissues can be attributed to proteins such as OsZIP4 and OsHMA2, which are localized on vascular bundles in nodes [18,24,31,32]. Wong et al. (2009) have also shown that HMA2 and HMA4 are principal transporters contributing to xylem loading of Cd [33,34]. These large complex protein families are often non-specific and can function in the transport of Cd and other divalent ions when available [26,30,35,36,37,38,39,40]. Shared cadmium and zinc transport proteins may influence their transport pathways due to interactions between the minerals.
There are a number of transporters specific for Cd or Cd-conjugated thiol ligands [41,42]. The subcellular sequestration of Cd can be achieved by ligand binding to glutathione or phytochelatin and subsequent transport into the vacuole by ABC-type vacuolar transporters homologous to those found in the yeast [43,44,45,46,47,48]. OsHMA3 is known to play a key role in the root sequestration of Cd, where a single amino acid substitution is associated with high grain Cd [15,32].
Studies on Cd–Zn interactions in rice are limited and remain inconclusive. Some studies have shown that Zn fertilization can ameliorate physiological stresses induced by Cd [49,50]. Zinc supplied to roots or leaves as a foliar spray reduces Cd accumulation in roots and grains [50,51,52,53]. However, other studies have shown that Cd translocation to shoots and accumulation in grains are actually increased as a result of Zn fertilization [54,55,56]. Cadmium and zinc interactions can influence the response of genes in different tissues, resulting in high and low grain Cd accumulators [28]. Tian et al. (2022) have shown that significant differences in the translocation rates of essential minerals under Cd stress were due to competitive interactions between minerals for HMA proteins [57]. Stress response genes, such as γ-glutamyl cysteine ligase (GSH ligase), glutathione synthase (GSHS), and phytochelatin synthase (PCS1), are involved in the biosynthesis of non-protein peptide chelators that bind Cd and Zn, and can also influence cellular transport [19,58,59,60,61]. Because the specifics of these interactions are unknown, there is a need to analyze whole plant mineral transport and gene expression patterns to elucidate the mechanisms responsible for low grain-Cd rice phenotypes. In this study, we sought to understand the dynamics of Cd transport and response to the increased supply of Zn, especially during the grain filling stage. The gene expression of select transport proteins, OsNRAMP5, OsHMA2, and OsHMA3, and stress response genes, GSH ligase, glutathione synthase, and phytochelatin synthase, were compared between the three lines of rice that differed in grain cadmium concentration so that we could investigate the contributions of these genes to the observed phenotypes. Furthermore, we analyzed other essential minerals to elucidate the overall effects of Cd and Zn availability on the transport of essential macro- and micronutrients. A comprehensive understanding of Cd–Zn interactions and the effect of Zn fertilization on Cd accumulation is important to understand the pathway of both Cd and Zn from root to shoot to seed, especially if Zn fertilization is used as a method to reduce grain Cd or increase grain Zn concentration in rice.

2. Materials and Methods

2.1. Sample Selection

Rice accessions (Oryza sativa) were obtained from the USDA-ARS, Dale Bumpers National Rice Research Center (Stuttgart, Arkansas, Genetic Stocks Oryza Collection-GSOR). Based on the grain Cd concentrations published in Pinson et al. (2015) [62], we selected three lines, PI 310546 (herein referred to as line 546), PI 311667 (herein referred to as line 667), and PI 301428 (herein referred to as line 428). Rice lines were grown in hydroponics (details below) containing a full nutrient solution and 1 mM Cd to ensure they could successfully reach reproductive maturity and accumulate Cd in grains. When grown hydroponically, rice lines 546, 667, and 428 accumulated low, medium, and high concentrations of Cd, respectively. As such, these three rice lines were chosen for further experimentation.

2.2. Plant Growth and Treatment

Seeds were surface sterilized with 1% hypochlorite for 10 min, imbibed overnight MilliQ H2O (Millipore, Billerica, MA, USA), and germinated in moistened filter paper lined Petri dishes for 5 days before being transplanted to hydroponics. A total of 96 germinated seedlings per line were transferred to 4 L pots (2 plants per pot/3 pots total for every line × treatment × reproductive stage permutation) (Supplementary Table S1). Each pot represented one of three experimental replicates. Hydroponic solution consisted of a modified Johnson nutrient solution with the following nutrients: 2 mM KNO3, 1 mM Ca(NO3)2, 1 mM KH2PO4, 1 mM MgSO4, 25 µM CaCl2, 25 µM H3BO3, 0.5 µM H2MoO4, 0.1 µM NiSO4, 2 µM MnSO4, 0.5 µM CuSO4, iron supplied as 20 µM Fe(III) HEDTA (N-(2-hydroxymethy) ethylenediaminetriacetic acid), and 2 mM 2-(N-morpholino) ethane-sulfonic acid (MES; Sigma Chemical, St. Louis, MO, USA) to buffer the solution to pH 6.0. Solutions were aerated continuously and replaced every week. The plants were grown in the Lehman College greenhouse with natural lighting plus supplemental metal halide lamps (100W) allowing for a 16 h day and 8 h night photoperiod, 25 ± 3 °C day/23 ± 3 °C night.
A completely randomized factorial design was used to examine Cd and Zn partitioning in three lines (line 546, line 667, and line 428) at different stages of development. When plants grew to the 4 to 5 leaf stage, a fresh solution was given where the plants were assigned to a treatment consisting of different Cd and Zn combinations: 0 µM Cd + 2 µM Zn (c0z2), 1 µM Cd + 2 µM Zn (c1z2), 0 µM Cd + 10 µM Zn (c0z10), or 1 µM Cd + 10 µM Zn (c1z10), respectively. Full solution changes were performed weekly to ensure roots were always submerged in growing solution and mineral availability was consistent over time. Cadmium was supplied as cadmium sulfate for the cadmium treatments and plants were harvested at 1 of 4 predetermined reproductive stages where “vegetative (Veg)” represents 30 days of growth, “Anthesis (Ant)” represents the flowering stage, which was determined to be the first day of heading, and 2 (2WAA) and 4 (4WAA) weeks after anthesis. At harvest, plant material was separated into roots, stems, lower leaves, flag leaves, peduncles, rachis, and grains when available (Supplementary Table S1). Roots were rinsed with deionized water before further processing. Grains were further separated from husks, but low husk weights prevented the tissue from being used in follow-up experiments. Plants from one pot were treated as one biological replicate. A total of three biological replicates were collected and analyzed for each line x treatment x harvest permutation (Supplementary Table S1). All plant material was roughly separated into equal parts, where half of the material was collected for mineral partitioning and the other half was flash frozen in liquid nitrogen to be used for gene expression analyses.

2.3. Tissue Digestion and Mineral Concentration Determination

Tissue for mineral partitioning analysis was collected and oven dried at 60 °C to constant mass and homogenized using a Wiley mill with a 0.2 mm screen (Thomas Scientific, Swedesboro, NJ, USA). For each tissue, ~0.25 g was digested using 3 mL of concentrated nitric acid at 100 °C in a BD50 Digestion Block (Seal Analytical, Kitchener, OT, Canada). Once cooled, 4 mL of 30% hydrogen peroxide was added and temperatures were periodically ramped up to 250 °C over the course of 6 h to dissolve all organic material. Digests were resuspended in 10 mL of 1% nitric acid and filtered using a 0.22 mm pre-syringe filter. The acids used were trace metal-grade (Fisher Scientific, Pittsburg, PA, USA) and the water used was deionized using the MilliQ system. Cadmium sulfate was used as the positive control and blanks were used for every digestion run. Samples were analyzed for nine essential minerals (Ca, Mg, P, S, K, Cu, Fe, Zn, Mn) and Cd by inductively coupled plasma optical emission spectroscopy (ICP-OES; iCAP 7000; Thermo Electron North America LLC, Madison, WI, USA), calibrated with certified standards to detect mineral concentrations as described previously [63]. Negative controls (blanks) were processed and analyzed to monitor for background levels of minerals in the digestions and resuspension solutions. Background mineral levels were subtracted from sample values before concentrations were calculated. The reported averages and standard error are derived from concentrations from three biological replicates (three pots of six plants).

2.4. Gene Expression

Tissue specific gene expression was analyzed to assess the contributions of transport proteins in Zn and Cd transport in different varieties of rice. All harvested tissue for gene expression analysis was instantly flash frozen using liquid nitrogen and crushed using a pre-chilled mortar and pestle. About 0.5 g of tissue was collected and combined with 600 mL of Trizol reagent (Invitrogen. Waltham, MA, USA). After 10 min, 125 mL of chloroform was added, tubes were mixed thoroughly, and chilled on ice for 10 min. Samples were then centrifuged for 20 min at 11,000 RPM. The supernatant was carefully separated and put through the RNeasy Plant RNA Extraction kit (Qiagen. Germantown, MD, USA), where RNA was diluted in 30 mL RNA Free H2O. RNA concentration and quality were analyzed using a Nano drop ND-2000 (Thermo Fisher Scientific, Waltham, MA, USA). DNA impurities were degraded using the RQ1 RNAse-Free DNAse kit (Promega. Madison, WI, USA). RNA (1 mg) was then used for first-strand synthesis RT-PCR with random primers using the GoScript Reverse Transcription kit (Promega. Madison, WI, USA). Primers designed using Integrated DNA Technologies primer design tools were then used (Supplementary Table S2) with SYBR green master mix (Promega. Madison, WI, USA) to perform quantitative real-time polymerase chain reaction. The transport genes OsNramp5, OsHMA2, and OsHMA3 and heavy metal chelator biosynthesis genes phytochelatin synthase (PCS1), γ-glutamyl cysteine ligase (GSH ligase), and glutathione synthase (GSHS) were assessed due to their importance in both zinc and cadmium transport in rice [4,58,59,64,65]. Reactions of two technical replicates and three biological replicates were carried out in a CFX Connect Real Time PCR System (BioRad, Hercules, CA, USA) in 15 mL reactions. Results were calculated as relative gene expression using the delta-delta Ct method with values normalized to the Ubiquitin housekeeping gene [66]. The reference gene was run on each plate to ensure stability. Reported values are average expression and standard error from n = 3 biological replicates.

2.5. Statistics

Analysis of variance and Pearson’s correlation coefficients were performed and determined using GraphPad Prism version 9.1 (GraphPad. San Diego, CA, USA). Two-way analysis of variance (ANOVA) multiple comparisons test with a Tukey correction were used to analyze statistical differences in mean tissue cadmium and mineral concentrations or relative gene expression between lines at a particular treatment and harvest or between treatments in a particular line and harvest (biological replicates n = 3). Significance of p ≤ 0.05 are annotated on graphs and tables. Pearson’s correlation coefficients were also calculated between Cd concentration and essential mineral concentrations. Significant correlations p < 0.05 are annotated on tables.

3. Results

3.1. Mineral Concentrations

3.1.1. Effect of Zinc on Cadmium Concentrations

All three lines, lines 546, 667, and 428, were able to set seed without phytotoxicity in all of the treatments analyzed. Cadmium concentrations in the different tissues of control plants (c0z2) are presented in Supplementary Table S3. Consistently under c1z2 (1 µM Cd + 2 µM Zn) treatments (Figure 1 and Figure 2), line 428 maintained the highest Cd concentrations in all tissues except for flag leaves (Figure 2), where 667 had significantly higher Cd concentrations. Line 546 demonstrated the lowest concentrations of Cd in most tissues compared to 667 and 428. Grain Cd concentrations at full maturity were statistically different (p < 0.05) between all three lines, averaging 0.29, 1.09, and 2.0 mg/g DW in lines 546, 667, and 428, respectively.
When grown with increased Zn levels (c1z10; 1 µM Cd + 10 µM Zn) (Figure 1 and Figure 2), there were no differences between the lines in root Cd concentration except the vegetative stage. The high line (428) had significantly higher Cd concentrations in stems, flag leaves, peduncles, rachis, and grains. There were no differences between the medium and the low Cd lines in most of the tissues. Cadmium concentrations in grains were 4.0, 1.3, and 0.8 μg/g DW in lines 428, 546, and 667, respectively. There were some differences between c1z2 and c1z10 in the vegetative tissues, but differences were not consistent (Table S4). It is interesting to note that the flag leaf Cd concentration was significantly reduced with increased Zn in the medium line (667) when compared to lower Zn treatments (Figure 2; Supplementary Table S5). Zinc fertilization significantly increased grain Cd concentration by 348% in the low line (546) and by 102% in the high line (428) but resulted in a slightly lower Cd concentration in the medium line (667), although the difference was not statistically significant.

3.1.2. Effect of Cadmium and Zinc Fertilization on Zinc Concentrations

Zinc concentrations were comparable between all lines under control (c0z2; 0 mM Cd + 2 mM Zn) treatments (Figure 3 and Figure 4). In the presence of cadmium (c1z2; 1 mM Cd + 2 mM Zn), the high line (428) had significantly higher grain Zn concentrations when compared to the other two lines. However, there were no differences in Zn concentrations between the control and the Cd-treated lines (Supplementary Tables S6 and S7).
Zinc fertilization increased Zn concentration in specific tissues but differed between lines. When extra Zn was supplied in the absence of Cd (c0z10; 0 mM Cd + 10 mM Zn), Zn concentrations were significantly higher when compared to the controls (c0z2; 0 mM Cd + 2 mM Zn) in most tissues (Figure 3 and Figure 4; Supplementary Tables S6 and S7). Although the Zn concentrations were higher in the Zn fertilization treatments, there were no differences between the lines for most tissues.
When Zn fertilization treatments with and without Cd (c0z10 vs. c1z10) were compared, although the root Zn concentrations were significantly lower in the high line (428) in the presence of Cd, there were no differences in most aboveground tissue Zn concentrations between c0z10 and c1z10 treatments (Supplementary Tables S6 and S7).
When grown in the presence of Cd (c1z2 and c1z10), Zn concentration patterns differed between the two treatments (Figure 3 and Figure 4). Zn fertilization (c1z10) led to significant increases in Zn in some tissues in all three lines, and significantly increased in the grains of both the low and high lines. Zinc concentrations were significantly higher in most of the vegetative tissues in the Zn fertilization treatment (c1z10) when compared to c1z2, but that difference was not consistent in the reproductive tissues (Supplementary Tables S6 and S7).

3.1.3. Other Grain Mineral Concentrations

The grain mineral concentrations of nine of the essential minerals were also measured in the three lines in all treatment conditions (Supplementary Table S8). The effects of Zn fertilization were unique to the line and mineral tested. As a general trend, line 428 accumulated the highest concentrations of calcium (Ca), copper (Cu), iron (Fe), magnesium (Mg), and manganese (Mn) under the varying treatments. For all lines, Zn fertilization alone caused no changes in the concentrations of Ca, Cu, Fe, Mg, Mn, Ni, or P, but Cd or Cd and Zn treatments did. Cadmium treatments (c1z2) increased the grain concentrations of Ca, Fe, and Mg in high-line (428) grains compared to controls (c0z2), but there was no difference in these minerals between the control (c0z2) and zinc fertilization (c0z10 and c1z10) treatments. Interestingly, Mg concentrations increased in the presence of Cd in the high line, while Mg concentration was lower in low and medium lines. Pearson’s correlation test revealed significant positive correlations between Cd and Ca and Cd and Ni, while there was significant negative correlation between Cd and P in the medium (667) line (Supplementary Table S9). When the plants were supplied with extra Zn (c1z10), there was strong positive correlation between Cd and Mn in the high Cd line.

3.2. Gene Expression

Gene expression relative to ubiquitin was analyzed in the roots of plants harvested in the vegetative stage (Figure 5) or roots and flag leaves of plants harvested at the anthesis stage (Figure 6). In the vegetative tissues, the gene expression of all transport genes examined were highest in line 667 when cadmium and/or zinc were supplied. At anthesis, however (Figure 6), OsNRAMP5 and OsHMA3 expression were higher in the low line (546) when extra Zn was supplied. In flag leaves, OsHMA3 expression was higher in the medium line (667) under c1z2 treatments (Figure 6).
The genes involved in the synthesis of non-protein thiols glutathione and phytochelatin were also analyzed in tissues harvested at the vegetative and anthesis stages. At anthesis, root GSH ligase expression was significantly higher in the high (428) line when compared to the other two lines (Figure 7A). PCS1 expression was increased in the roots of line 667 as a response to all treatments but only increased in line 546 roots as a response to Zn (Figure 7B). In flag leaves, GSH ligase, GSHS and PCS1 were consistently highest in the medium line (667) (Figure 7C–E) under c1z2 treatments.

4. Discussion

Interactions between minerals can significantly alter their distribution and accumulation in plant tissues [67,68,69,70]. Given the importance of rice as a major food source, it is crucial to produce low Cd cultivars to meet the demand of customers. There is an urgent need to elucidate the mechanisms responsible for differentiating between toxic cadmium and beneficial zinc in uptake, sequestration, and bioavailability. The goal of this work was to assess the influence of Zn levels on Cd and mineral transport and partitioning in three varieties of rice that differed in grain Cd accumulation. To accomplish this goal, we used Cd and Zn treatments that were similar to the environmental concentrations available to plants [71]. All three rice lines were able to successfully reach full reproductive maturity without showing any signs of phytotoxicity.
Our study showed that there was no difference in root Cd concentrations in response to Zn fertilization in all the lines. Some previous studies have shown that Zn fertilization can reduce Cd translocation and accumulation in different tissues and grains [72,73,74,75]. However, our study showed that when increasing the external zinc concentration (c1z10), grain Cd concentrations increased in the high (428) and the low (546) lines. This suggests that either there is a different route of translocation for Cd and Zn or that these concentrations of Zn are not effective in restricting Cd uptake and translocation. This is consistent with some previous studies that showed that Zn does not always competitively inhibit Cd accumulation and it is concentration-dependent [31,76,77,78]. Interestingly, in the medium line, Cd concentration in the flag leaves was significantly lower in response to Zn fertilization and in turn did not cause an increase in grain Cd, which seems to indicate that Zn competitively inhibits Cd transport in this line. Most of the cadmium was accumulated in the stems and lower leaves, consistent with previous studies in cereal crops due to increased transpiration-driven flow in these tissues [48,67,79]. Our study showed only some remobilization of cadmium from the vegetative tissues, which indicates that Cd could be stored in non-labile subcellular compartments such as the vacuoles and thus is not loaded into the phloem [48,67,80].
In terms of the Zn concentration, our results showed that Cd did not influence Zn concentration in any of the lines. It was also interesting to note that in Zn fertilization treatments without Cd (c0z10), there was an increase in Zn concentration in many of the tissues but there was no increase in the grain Zn concentration, especially in the low and medium lines. A previous study has also shown that when Zn was supplied to the leaves of rice during flowering and grain fill, 90% of this Zn was translocated to the vegetative organs and not to the grains [81].
The route that a mineral takes to reach grains is influenced by several factors, including mineral availability, interactions with other minerals, and gene expression, which may be specific to genotype [82]. Transport proteins are critical for the uptake, translocation, and sequestration of minerals in plant tissues [83,84,85]. The regulation of transport proteins at the molecular level can be used to identify the critical mechanisms involved in mineral transport to edible tissues of plants [86,87,88]. The expression of transport proteins can be regulated by mineral availability [89,90,91]. In order to assess Cd and Zn interactions in mineral transport, we analyzed how the gene expression of key transport proteins was affected by Cd, Zn, or combination treatments. While several transport proteins have been identified, we focused on key transport proteins that have been implicated in Cd or Zn transport in rice. In the high-, medium-, and low-grain Cd-accumulating lines, the expressions of OsNRAMP5, OsHMA2, and OsHMA3 were all influenced by treatment (Figure 5). Additionally, expression patterns were unique at reproductive stages. OsNRAMP5 is a root uptake transport protein involved in the transport of cadmium, manganese, and iron [65,92]. In this study, although OsNRAMP5 expression was the highest in the medium line or low line depending on the treatment and developmental stage, there was no difference in root Cd concentrations between the lines. At any given time, minerals are being translocated out of the roots to shoots. Thus, similar Cd and Zn concentrations in roots could be due to the differential expression of other genes as well. One such gene, OsHMA2, is a xylem loading transport protein that has been shown to be essential for Zn homeostasis and root to shoot translocation of both Zn and Cd in rice [32]. Mutations in OsHMA2 were responsible for restricting the translocation of both minerals in rice [93]. The low expression of OsHMA2 in the low Cd-accumulating line, line 546, could explain the reduced Cd in aerial tissues. In rice roots, OsHMA3, is localized to the tonoplast and responsible for vacuolar sequestration of Cd and excess Zn, restricting their transport to aboveground tissues [15,16,18,94,95,96]. In our study, the expression of OsHMA3 was also found to be significantly higher in the flag leaves of the medium accumulating line, 667 (Figure 5). Interestingly, this can probably explain the difference between the high (428) and medium (667) Cd lines, which differed in grain Cd concentration in both c1z2 and c1z10 treatments. Transport proteins including OsZIP3, OsZIP4, Cation Exchange Proteins (CAX), Low Cation Transporters (LCT), Metal Tolerance Proteins (MTP), and ATP Binding Cassette transporters (ABCC) have also been implicated in the transport of Cd and/or Zn, but have been omitted from this study. We did, however, analyze the genes involved in the production of non-protein peptides, glutathione and phytochelatin, which assist in the detoxification of Cd [97]. The three genes tested, GSH ligase, GSHS, and PCS1 were explored in flag leaves and roots at the vegetative and anthesis stages (Figure 6). The expression patterns of all genes were affected by treatment in one or more lines. The differential expression between the three lines suggests that Cd–Zn interactions can influence transport by modifying production of thiolate compounds [98,99,100]. It is interesting to note that the expression of the three genes in the medium-line (667) c1z2 treatment was significantly higher in the flag leaves when compared to the other two lines, which correlates with the high Cd concentration in the flag leaves. This could be the reason for restricted translocation of cadmium from the flag leaf to the reproductive tissues since the concentrations of Cd in the peduncles and rachis were similar to the other two lines. Experiments analyzing the role of glutathione and phytochelatin in cadmium retention suggest that the non-protein peptides play a key role in limiting the translocation of the toxin [61,98].
Concentrations of essential minerals Ca, Cu, Fe, Mg, Mn, Ni, and P were also significantly affected by the Cd and Zn treatments, although effects were specific to the line tested. The high line (428) had the highest concentrations of Ca, Cu, Fe, Mg, and Mn in c1z2 treatments. When Zn fertilization was supplied, Ca, Fe, and P concentrations were lower, but Mn and Ni concentrations in grains increased. We also found strong positive and negative correlations between Cd and essential minerals, but, once again, the specific mineral interactions differed between the lines. Positive or negative correlations between the minerals may be the result of shared transport systems [20,24,36,101]. Our results confirm the importance of assessing the interactions between other essential minerals in attempts to biofortify for specific traits.
Limitations of this work include the inability to measure mineral content or measure precise tissue accumulations or translocation patterns since we used part of the harvested tissue for gene expression studies. We also analyzed the gene expression of only selected transport proteins and stress genes that have been implicated in both cadmium and zinc transport in rice. More robust techniques such as RNA Seq may be able to identify even more genetic differences controlling differential Cd and Zn accumulation.

5. Conclusions

In conclusion, there is a dire need to produce cereals with improved nutritional quality and, to do so, a detailed understanding of mineral interactions on transport is needed. We have shown here that increased Zn supply does not decrease Cd concentration but in fact increases grain Cd concentration and that an increased Zn supply does not increase grain Zn concentrations in these three rice lines. Cadmium and zinc tissue concentrations differed between the lines at different stages of reproductive maturity. These accumulation patterns are likely due to differential gene expression controlling mineral translocation within the plant. The difference in the expression patterns of OsNRAMP5, OsHMA2, OsHMA3, GSH ligase, GSHS, and PCS in our study could be used to explain the difference in Cd translocation and accumulation between the lines. In particular, our study indicates the need to identify genes specifically involved in the export of minerals from flag leaves to grains, which may help us understand the route of mineral translocation into the grains, further assisting in improving strategies to minimize Cd while increasing nutrient density in rice.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/agronomy12092182/s1. Table S1: A and B: Experimental design: List of lines, treatments, reproductive stages, and tissues analyzed in experiment; Table S2: List of primers used; Table S3: Cadmium concentrations in (A) vegetative and (B) reproductive tissues of 3 rice lines grown in the absence of cadmium; Table S4: Effect of zinc fertilization on cadmium concentration in vegetative tissues; Table S5: Effect of zinc fertilization on cadmium concentration in reproductive tissues; Table S6: Effect of cadmium and zinc fertilization on zinc concentration in vegetative tissues; Table S7: Effect of cadmium and zinc fertilization on zinc concentration in reproductive tissues; Table S8: Mineral concentrations in grains of 3 rice cultivars grown c0z2, c1z2, c0z10, and c1z10; Table S9: Correlations between grain Cd and essential minerals.

Author Contributions

R.P.S. conceived and designed the experiments; M.T. performed the experiments and analyzed the data; M.A.G. performed the ICP-OES analysis and contributed to the manuscript; M.T. and R.P.S. wrote the paper. This work is part of the doctoral dissertation work by M.T.: Understanding the influence of zinc on grain cadmium accumulation and bioaccessibility in rice. 2021. Doctoral dissertation, City University of New York. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by a PSC-CUNY Research Award, project TRADA-45-605 to R.P.S. and by the US Department of Agriculture, Agricultural Research Service, project 3060-21650-001-00D to M.A.G.

Acknowledgments

The authors would like to thank Jahnelle Howe, Yardley Martinez, and Thomas Weitzman for their assistance in the greenhouse. The authors also thank Veselina Valkov for assistance with the ICP-OES mineral analyses. The contents of this publication do not necessarily reflect the views or policies of the US Department of Agriculture nor does the mention of trade names, commercial products, or organizations imply endorsement by the US Government.

Conflicts of Interest

The authors declare that this research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Hayat, M.T.; Nauman, M.; Nazir, N.; Ali, S.; Bangash, N. Environmental Hazards of Cadmium: Past, Present, and Future. In Cadmium Toxicity and Tolerance in Plants: From Physiology to Remediation; Elsevier: Amsterdam, The Netherlands, 2018; pp. 163–183. ISBN 9780128148655. [Google Scholar]
  2. Proshad, R.; Kormoker, T.; Mursheed, N.; Monirul Islam, M.; Bhuyan, M.I.; Sazedul Islam, M.; Mithu, T.N. Heavy metal toxicity in agricultural soil due to rapid industrialization in Bangladesh: A review. Int. J. Adv. Geosci. 2018, 6, 83. [Google Scholar] [CrossRef]
  3. Kubier, A.; Wilkin, R.T.; Pichler, T. Cadmium in soils and groundwater: A review. Appl. Geochem. 2019, 108, 104388. [Google Scholar] [CrossRef] [PubMed]
  4. Sui, F.Q.; Chang, J.D.; Tang, Z.; Liu, W.J.; Huang, X.Y.; Zhao, F.J. Nramp5 expression and functionality likely explain higher cadmium uptake in rice than in wheat and maize. Plant Soil 2018, 433, 377–389. [Google Scholar] [CrossRef]
  5. Wu, D.; Yamaji, N.; Yamane, M.; Kashino-Fujii, M.; Sato, K.; Ma, J.F. The HvNramp5 transporter mediates uptake of cadmium and manganese, but not iron. Plant Physiol. 2016, 172, 1899–1910. [Google Scholar] [CrossRef]
  6. Chang, J.D.; Huang, S.; Konishi, N.; Wang, P.; Chen, J.; Huang, X.Y.; Ma, J.F.; Zhao, F.J. Overexpression of the manganese/cadmium transporter OsNRAMP5 reduces cadmium accumulation in rice grain. J. Exp. Bot. 2020, 71, 5705–5715. [Google Scholar] [CrossRef] [PubMed]
  7. Clemens, S.; Ma, J.F. Toxic Heavy Metal and Metalloid Accumulation in Crop Plants and Foods. Annu. Rev. Plant Biol. 2016, 67, 489–512. [Google Scholar] [CrossRef] [PubMed]
  8. Nishijo, M.; Nakagawa, H.; Suwazono, Y.; Nogawa, K.; Kido, T. Causes of death in patients with Itai-itai disease suffering from severe chronic cadmium poisoning: A nested case-control analysis of a follow-up study in Japan. BMJ Open 2017, 7, e015694. [Google Scholar] [CrossRef]
  9. Song, Y.; Wang, Y.; Mao, W.; Sui, H.; Yong, L.; Yang, D.; Jiang, D.; Zhang, L.; Gong, Y. Dietary cadmium exposure assessment among the Chinese population. PLoS ONE 2017, 12, e0177978. [Google Scholar] [CrossRef]
  10. Fujiwara, Y.; Lee, J.Y.; Tokumoto, M.; Satoh, M. Cadmium renal toxicity via apoptotic pathways. Biol. Pharm. Bull. 2012, 35, 1892–1897. [Google Scholar] [CrossRef]
  11. European Food Safety Authority. Cadmium dietary exposure in the European population. EFSA J. 2012, 10, 2551. [Google Scholar] [CrossRef]
  12. Waters, B.M.; Sankaran, R.P. Moving micronutrients from the soil to the seeds: Genes and physiological processes from a biofortification perspective. Plant Sci. 2011, 180, 562–574. [Google Scholar] [CrossRef] [PubMed]
  13. Yoneyama, T.; Gosho, T.; Kato, M.; Goto, S.; Hayashi, H. Xylem and phloem transport of Cd, Zn and Fe into the grains of rice plants (Oryza sativa L.) grown in continuously flooded Cd-contaminated soil. Soil Sci. Plant Nutr. 2010, 56, 445–453. [Google Scholar] [CrossRef]
  14. Ishikawa, S.; Suzui, N.; Ito-Tanabata, S.; Ishii, S.; Igura, M.; Abe, T.; Kuramata, M.; Kawachi, N.; Fujimaki, S. Real-time imaging and analysis of differences in cadmium dynamics in rice cultivars (Oryza sativa) using positron-emitting 107Cd tracer. BMC Plant Biol. 2011, 11, 172. [Google Scholar] [CrossRef]
  15. Ueno, D.; Yamaji, N.; Kono, I.; Huang, C.F.; Ando, T.; Yano, M.; Ma, J.F. Gene limiting cadmium accumulation in rice. Proc. Natl. Acad. Sci. USA 2010, 107, 16500–16505. [Google Scholar] [CrossRef] [PubMed]
  16. Miyadate, H.; Adachi, S.; Hiraizumi, A.; Tezuka, K.; Nakazawa, N.; Kawamoto, T.; Katou, K.; Kodama, I.; Sakurai, K.; Takahashi, H.; et al. OsHMA3, a P1B-type of ATPase affects root-to-shoot cadmium translocation in rice by mediating efflux into vacuoles. New Phytol. 2011, 189, 190–199. [Google Scholar] [CrossRef]
  17. Ma, J.F.; Shen, R.F.; Shao, J.F. Transport of cadmium from soil to grain in cereal crops: A review. Pedosphere 2021, 31, 3–10. [Google Scholar] [CrossRef]
  18. Rasheed, A.; Fahad, S.; Aamer, M.; Hassan, M.U.; Tahir, M.M.; Wu, Z.M. Role of genetic factors in regulating cadmium uptake, transport and accumulation mechanisms and quantitative trait loci mapping in rice. A review. Appl. Ecol. Environ. Res. 2020, 18, 4005–4023. [Google Scholar] [CrossRef]
  19. Guo, J.; Zhang, X.; Ye, D.; Huang, H.; Wang, Y.; Zheng, Z.; Li, T.; Yu, H. Crucial roles of cadmium retention in nodeⅡ for restraining cadmium transport from straw to ear at reproductive period in a grain low-cadmium rice line (Oryza sativa L.). Ecotoxicol. Environ. Saf. 2020, 205, 111323. [Google Scholar] [CrossRef]
  20. Sasaki, A.; Yamaji, N.; Xia, J.; Ma, J.F. OsYSl6 is involved in the detoxification of excess manganese in rice. Plant Physiol. 2011, 157, 1832–1840. [Google Scholar] [CrossRef]
  21. Songmei, L.; Jie, J.; Yang, L.; Jun, M.; Shouling, X.; Yuanyuan, T.; Youfa, L.; Qingyao, S.; Jianzhong, H. Characterization and Evaluation of OsLCT1 and OsNramp5 Mutants Generated Through CRISPR/Cas9-Mediated Mutagenesis for Breeding Low Cd Rice. Rice Sci. 2019, 26, 88–97. [Google Scholar] [CrossRef]
  22. White, P.J.; Whiting, S.N.; Baker, A.J.M.; Broadley, M.R. Does zinc move apoplastically to the xylem in roots of Thlaspi caerulescens? New Phytol. 2002, 153, 201–207. [Google Scholar] [CrossRef]
  23. Ricachenevsky, F.K.; De Araújo Junior, A.T.; Fett, J.P.; Sperotto, R.A. You shall not pass: Root vacuoles as a symplastic checkpoint for metal translocation to shoots and possible application to grain nutritional quality. Front. Plant Sci. 2018, 9, 412. [Google Scholar] [CrossRef] [PubMed]
  24. Chen, W.R.; Feng, Y.; Chao, Y.E. Genomic analysis and expression pattern of OsZIP1, OsZIP3, and OsZIP4 in two rice (Oryza sativa L.) genotypes with different zinc efficiency. Russ. J. Plant Physiol. 2008, 55, 400–409. [Google Scholar] [CrossRef]
  25. Chowdhury, R.; Nallusamy, S.; Shanmugam, V.; Loganathan, A.; Muthurajan, R.; Sivathapandian, S.K.; Paramasivam, J.; Duraialagaraja, S. Genome-wide understanding of evolutionary and functional relationships of rice Yellow Stripe-Like (YSL) transporter family in comparison with other plant species. Biologia 2022, 77, 39–53. [Google Scholar] [CrossRef]
  26. Huang, S.; Yamaji, N.; Feng Ma, J. Zinc transport in rice: How to balance optimal plant requirements and human nutrition. J. Exp. Bot. 2022, 73, 1800–1808. [Google Scholar] [CrossRef] [PubMed]
  27. Bandyopadhyay, T.; Mehra, P.; Hairat, S.; Giri, J. Morpho-physiological and transcriptome profiling reveal novel zinc deficiency-responsive genes in rice. Funct. Integr. Genom. 2017, 17, 565–581. [Google Scholar] [CrossRef]
  28. Adil, M.F.; Sehar, S.; Chen, G.; Chen, Z.H.; Jilani, G.; Chaudhry, A.N.; Shamsi, I.H. Cadmium-zinc cross-talk delineates toxicity tolerance in rice via differential genes expression and physiological / ultrastructural adjustments. Ecotoxicol. Environ. Saf. 2020, 190, 110076. [Google Scholar] [CrossRef]
  29. Ramesh, S.A.; Shin, R.; Eide, D.J.; Schachtman, D.P. Differential metal selectivity and gene expression of two zinc transporters from rice. Plant Physiol. 2003, 133, 126–134. [Google Scholar] [CrossRef]
  30. Huang, S.; Sasaki, A.; Yamaji, N.; Okada, H.; Mitani-Ueno, N.; Ma, J.F. The ZIP transporter family member OsZIP9 contributes to root zinc uptake in rice under zinc-limited conditions. Plant Physiol. 2020, 183, 1224–1234. [Google Scholar] [CrossRef]
  31. Fontanili, L.; Lancilli, C.; Suzui, N.; Dendena, B.; Yin, Y.G.; Ferri, A.; Ishii, S.; Kawachi, N.; Lucchini, G.; Fujimaki, S.; et al. Kinetic Analysis of Zinc/Cadmium Reciprocal Competitions Suggests a Possible Zn-Insensitive Pathway for Root-to-Shoot Cadmium Translocation in Rice. Rice 2016, 9, 16. [Google Scholar] [CrossRef] [Green Version]
  32. Takahashi, R.; Ishimaru, Y.; Shimo, H.; Ogo, Y.; Senoura, T.; Nishizawa, N.K.; Nakanishi, H. The OsHMA2 transporter is involved in root-to-shoot translocation of Zn and Cd in rice. Plant Cell Environ. 2012, 35, 1948–1957. [Google Scholar] [CrossRef]
  33. Wong, C.K.E.; Cobbett, C.S. HMA P-type ATPases are the major mechanism for root-to-shoot Cd translocation in Arabidopsis thaliana. New Phytol. 2009, 181, 71–78. [Google Scholar] [CrossRef]
  34. Wong, C.K.E.; Jarvis, R.S.; Sherson, S.M.; Cobbett, C.S. Functional analysis of the heavy metal binding domains of the Zn/Cd-transporting ATPase, HMA2, in Arabidopsis thaliana. New Phytol. 2009, 181, 79–88. [Google Scholar] [CrossRef]
  35. Liu, X.S.; Feng, S.J.; Zhang, B.Q.; Wang, M.Q.; Cao, H.W.; Rono, J.K.; Chen, X.; Yang, Z.M. OsZIP1 functions as a metal efflux transporter limiting excess zinc, copper and cadmium accumulation in rice. BMC Plant Biol. 2019, 19, 283. [Google Scholar] [CrossRef]
  36. Yang, M.; Li, Y.; Liu, Z.; Tian, J.; Liang, L.; Qiu, Y.; Wang, G.; Du, Q.; Cheng, D.; Cai, H.; et al. A high activity zinc transporter OsZIP9 mediates zinc uptake in rice. Plant J. 2020, 103, 1695–1709. [Google Scholar] [CrossRef]
  37. Lan, H.X.; Wang, Z.F.; Wang, Q.H.; Wang, M.M.; Bao, Y.M.; Huang, J.; Zhang, H.S. Characterization of a vacuolar zinc transporter OZT1 in rice (Oryza sativa L.). Mol. Biol. Rep. 2013, 40, 1201–1210. [Google Scholar] [CrossRef]
  38. Khokhar, J.S.; King, J.; King, I.P.; Young, S.D.; Foulkes, M.J.; De Silva, J.; Weerasinghe, M.; Mossa, A.; Griffiths, S.; Riche, A.B.; et al. Novel sources of variation in grain Zinc (Zn) concentration in bread wheat germplasm derived from Watkins landraces. PLoS ONE 2020, 15, e0229107. [Google Scholar] [CrossRef]
  39. Yang, Y.; Li, Y.; Chen, W.; Wang, M.; Wang, T.; Dai, Y. Dynamic interactions between soil cadmium and zinc affect cadmium phytoavailability to rice and wheat: Regional investigation and risk modeling. Environ. Pollut. 2020, 267, 115613. [Google Scholar] [CrossRef]
  40. Zhang, M.; Liu, B. Identification of a rice metal tolerance protein OsMTP11 as a manganese transporter. PLoS ONE 2017, 12, e0174987. [Google Scholar] [CrossRef]
  41. Mendoza-Cózatl, D.G.; Butko, E.; Springer, F.; Torpey, J.W.; Komives, E.A.; Kehr, J.; Schroeder, J.I. Identification of high levels of phytochelatins, glutathione and cadmium in the phloem sap of Brassica napus. A role for thiol-peptides in the long-distance transport of cadmium and the effect of cadmium on iron translocation. Plant J. 2008, 54, 249–259. [Google Scholar] [CrossRef] [Green Version]
  42. Pike, S.; Patel, A.; Stacey, G.; Gassmann, W. Arabidopsis OPT6 is an oligopeptide transporter with exceptionally broad substrate specificity. Plant Cell Physiol. 2009, 50, 1923–1932. [Google Scholar] [CrossRef]
  43. Salt, D.E.; Rauser, W.E. MgATP-dependent transport of phytochelatins across the tonoplast of oat roots. Plant Physiol. 1995, 107, 1293–1301. [Google Scholar] [CrossRef]
  44. Ortiz, D.F.; Ruscitti, T.; McCue, K.F.; Ow, D.W. Transport of metal-binding peptides by HMT1, a fission yeast ABC-type vacuolar membrane protein. J. Biol. Chem. 1995, 270, 4721–4728. [Google Scholar] [CrossRef]
  45. Li, Z.S.; Szczypka, M.; Lu, Y.P.; Thiele, D.J.; Rea, P.A. The yeast cadmium factor protein (YCF1) is a vacuolar glutathione S-conjugate pump. J. Biol. Chem. 1996, 271, 6509–6517. [Google Scholar] [CrossRef]
  46. Morel, M.; Crouzet, J.; Gravot, A.; Auroy, P.; Leonhardt, N.; Vavasseur, A.; Richaud, P. AtHMA3, a P1B-ATPase allowing Cd/Zn/Co/Pb vacuolar storage in Arabidopsis. Plant Physiol. 2009, 149, 894–904. [Google Scholar] [CrossRef]
  47. Li, Q.; Guo, J.; Zhang, X.; Yu, H.; Huang, F.; Zhang, L.; Zhang, M.; Li, T. Changes of non-protein thiols in root and organic acids in xylem sap involved in cadmium translocation of cadmium-safe rice line (Oryza sative L.). Plant Soil 2019, 439, 475–486. [Google Scholar] [CrossRef]
  48. Hart, J.J.; Welch, R.M.; Norvell, W.A.; Kochian, L.V. Characterization of cadmium uptake, translocation and storage in near-isogenic lines of durum wheat that differ in grain cadmium concentration. New Phytol. 2006, 172, 261–271. [Google Scholar] [CrossRef]
  49. Rizwan, M.; Ali, S.; Adrees, M.; Rizvi, H.; Zia-ur-Rehman, M.; Hannan, F.; Qayyum, M.F.; Hafeez, F.; Ok, Y.S. Cadmium stress in rice: Toxic effects, tolerance mechanisms, and management: A critical review. Environ. Sci. Pollut. Res. 2016, 23, 17859–17879. [Google Scholar] [CrossRef]
  50. Adil, M.F.; Sehar, S.; Han, Z.; Wa Lwalaba, J.L.; Jilani, G.; Zeng, F.; Chen, Z.H.; Shamsi, I.H. Zinc alleviates cadmium toxicity by modulating photosynthesis, ROS homeostasis, and cation flux kinetics in rice. Environ. Pollut. 2020, 265, 114979. [Google Scholar] [CrossRef]
  51. Lv, G.; Wang, H.; Xu, C.; Shuai, H.; Luo, Z.; Zhang, Q.; Zhu, H.; Wang, S.; Zhu, Q.; Zhang, Y.; et al. Effectiveness of simultaneous foliar application of Zn and Mn or P to reduce Cd concentration in rice grains: A field study. Environ. Sci. Pollut. Res. 2019, 26, 9305–9313. [Google Scholar] [CrossRef]
  52. Cai, Y.Y.; Xu, W.; Wang, M.; Chen, W.; Li, X.; Li, Y.; Cai, Y.Y. Mechanisms and uncertainties of Zn supply on regulating rice Cd uptake. Environ. Pollut. 2019, 253, 959–965. [Google Scholar] [CrossRef]
  53. Mutemachani, M.J.; Tichaona, S.; Shafaque, S.; Xin, Z.; Yuqing, H.; Gerald, Z.; Antony, M.; Lwalaba Wa Lwalaba, J.; Haider, S.I. Effects of zinc and silicon on cadmium toxicity and mineral element translocation in two rice (Oryza sativa) genotypes. J. Zhejiang Univ. 2018, 44, 294–310. [Google Scholar] [CrossRef]
  54. Herath, D.; Weerasinghe, A.; Bandara, D.; Wijayawardhana, D. Synergistic Effect of Zinc and Cadmium for Uptake, Accumulation and Growth Responses. Int. J. Chem. Environ. Biol. Sci. 2016, 4, 69–73. [Google Scholar]
  55. Liu, H.J.; Zhang, J.L.; Christie, P.; Zhang, F.S. Influence of external zinc and phosphorus supply on Cd uptake by rice (Oryza sativa L.) seedlings with root surface iron plaque. Plant Soil 2007, 300, 105–115. [Google Scholar] [CrossRef]
  56. Wang, M.; Yang, Y.; Chen, W. Manganese, Zinc, and pH Affect Cadmium Accumulation in Rice Grain under Field Conditions in Southern China. J. Environ. Qual. 2018, 47, 306–311. [Google Scholar] [CrossRef]
  57. Tian, P.; Feng, Y.-X.; Li, C.-Z.; Zhang, P.; Yu, X.-Z. Transcriptional analysis of heavy metal P1B-ATPases (HMAs) elucidates competitive interaction in metal transport between cadmium and mineral elements in rice plants. Environ. Sci. Pollut. Res. 2022. [CrossRef]
  58. Cai, Y.; Cao, F.; Cheng, W.; Zhang, G.; Wu, F. Modulation of exogenous glutathione in phytochelatins and photosynthetic performance against Cd stress in the two rice genotypes differing in Cd tolerance. Biol. Trace Elem. Res. 2011, 143, 1159–1173. [Google Scholar] [CrossRef]
  59. Uraguchi, S.; Tanaka, N.; Hofmann, C.; Abiko, K.; Ohkama-Ohtsu, N.; Weber, M.; Kamiya, T.; Sone, Y.; Nakamura, R.; Takanezawa, Y.; et al. Phytochelatin synthase has contrasting effects on cadmium and arsenic accumulation in rice grains. Plant Cell Physiol. 2017, 58, 1730–1742. [Google Scholar] [CrossRef]
  60. Huang, H.; Li, M.; Rizwan, M.; Dai, Z.; Yuan, Y.; Hossain, M.M.; Cao, M.; Xiong, S.; Tu, S. Synergistic effect of silicon and selenium on the alleviation of cadmium toxicity in rice plants. J. Hazard. Mater. 2021, 401, 123393. [Google Scholar] [CrossRef]
  61. Wang, K.; Yu, H.; Zhang, X.; Ye, D.; Huang, H.; Wang, Y.; Zheng, Z.; Li, T. A transcriptomic view of cadmium retention in roots of cadmium-safe rice line (Oryza sativa L.). J. Hazard. Mater. 2021, 418, 126379. [Google Scholar] [CrossRef]
  62. Pinson, S.R.M.; Tarpley, L.; Yan, W.; Yeater, K.; Lahner, B.; Yakubova, E.; Huang, X.Y.; Zhang, M.; Guerinot, M.L.; Salt, D.E. Worldwide Genetic Diversity for Mineral Element Concentrations in Rice Grain. Crop Sci. 2015, 55, 294–311. [Google Scholar] [CrossRef] [Green Version]
  63. Farnham, M.W.; Keinath, A.P.; Grusak, M.A. Mineral concentration of broccoli florets in relation to year of cultivar release. Crop Sci. 2011, 51, 2721–2727. [Google Scholar] [CrossRef]
  64. Satoh-Nagasawa, N.; Mori, M.; Sakurai, K.; Takahashi, H.; Watanabe, A.; Akagi, H. Functional relationship heavy metal P-type ATPases (OsHMA 2 and OsHMA3) of rice (Oryza sativa) using RNAi. Plant Biotechnol. 2013, 30, 511–515. [Google Scholar] [CrossRef]
  65. Ishimaru, Y.; Takahashi, R.; Bashir, K.; Shimo, H.; Senoura, T.; Sugimoto, K.; Ono, K.; Yano, M.; Ishikawa, S.; Arao, T.; et al. Characterizing the role of rice NRAMP5 in Manganese, Iron and Cadmium Transport. Sci. Rep. 2012, 2, 286. [Google Scholar] [CrossRef]
  66. Rao, X.; Huang, X.; Zhou, Z.; Lin, X. An improvement of the 2^(-delta delta CT) method for quantitative real-time polymerase chain reaction data analysis. Biostat. Bioinforma. Biomath. 2013, 3, 71–85. [Google Scholar]
  67. Tavarez, M.; Macri, A.; Sankaran, R. Cadmium and zinc partitioning and accumulation during grain filling in two near isogenic lines of durum wheat. Plant Physiol. Biochem. 2015, 97, 461–469. [Google Scholar] [CrossRef]
  68. Guttieri, M.J.; Baenziger, S.P.; Frels, K.; Carver, B.; Arnall, B.; Wang, S.; Akhunov, E.; Waters, B.M. Prospects for selecting wheat with increased zinc and decreased cadmium concentration in grain. Crop Sci. 2015, 55, 1712–1728. [Google Scholar] [CrossRef]
  69. Liu, J.; Li, K.; Xu, J.; Liang, J.; Lu, X.; Yang, J.; Zhu, Q. Interaction of Cd and five mineral nutrients for uptake and accumulation in different rice cultivars and genotypes. Field Crops Res. 2003, 83, 271–281. [Google Scholar] [CrossRef]
  70. Sarwar, N.; Ishaq, W.; Farid, G.; Shaheen, M.R.; Imran, M.; Geng, M.; Hussain, S. Zinc-cadmium interactions: Impact on wheat physiology and mineral acquisition. Ecotoxicol. Environ. Saf. 2015, 122, 528–536. [Google Scholar] [CrossRef]
  71. Hart, J.J.; Welch, R.M.; Norvell, W.A.; Clarke, J.M.; Kochian, L.V. Zinc effects on cadmium accumulation and partitioning in near-isogenic lines of durum wheat that differ in grain cadmium concentration. New Phytol. 2005, 167, 391–401. [Google Scholar] [CrossRef]
  72. Harris, N.S.; Taylor, G.J. Remobilization of cadmium in maturing shoots of near isogenic lines of durum wheat that differ in grain cadmium accumulation. J. Exp. Bot. 2001, 52, 1473–1481. [Google Scholar] [CrossRef]
  73. Hart, J.J.; Welch, R.M.; Norvell, W.A.; Kochian, L.V. Transport interactions between cadmium and zinc in roots of bread and durum wheat seedlings. Physiol. Plant. 2002, 116, 73–78. [Google Scholar] [CrossRef]
  74. Hassan, M.J.; Zhang, G.; Wu, F.; Wei, K.; Chen, Z. Zinc alleviates growth inhibition and oxidative stress caused by cadmium in rice. J. Plant Nutr. Soil Sci. 2005, 168, 255–261. [Google Scholar] [CrossRef]
  75. Grant, C.A.; Bailey, L.D. Nitrogen, phosphorus and zinc management effects on grain yield and cadmium concentration in two cultivars of durum wheat. Can. J. Plant Sci. 1998, 78, 63–70. [Google Scholar] [CrossRef]
  76. Green, C.E.; Chaney, R.L.; Bouwkamp, J. Increased zinc supply does not inhibit cadmium accumulation by rice (Oryza sativa L.). J. Plant Nutr. 2017, 40, 869–877. [Google Scholar] [CrossRef]
  77. Zhao, A.Q.; Tian, X.H.; Lu, W.H.; Gale, W.J.; Lu, X.C.; Cao, Y.X. Effect of zinc on cadmium toxicity in winter wheat. J. Plant Nutr. 2011, 34, 1372–1385. [Google Scholar] [CrossRef]
  78. Ishfaq, M.; Kiran, A.; Khaliq, A.; Cheema, S.A.; Alaraidh, I.A.; Hirotsu, N.; Wakeel, A. Zinc biofortified wheat cultivar lessens grain cadmium accumulation under cadmium contaminated conditions. Int. J. Agric. Biol. 2018, 20, 2842–2846. [Google Scholar] [CrossRef]
  79. Harris, N.S.; Taylor, G.J. Cadmium uptake and partitioning in durum wheat during grain filling. BMC Plant Biol. 2013, 13, 103. [Google Scholar] [CrossRef]
  80. Sankaran, R.P.; Ebbs, S.D. Transport of Cd and Zn to seeds of Indian mustard (Brassica juncea) during specific stages of plant growth and development. Physiol. Plant. 2008, 132, 69–78. [Google Scholar] [CrossRef]
  81. Wu, C.Y.; Lu, L.L.; Yang, X.E.; Feng, Y.; Wei, Y.Y.; Hao, H.L.; Stoffella, P.J.; He, Z.L. Uptake, translocation, and remobilization of zinc absorbed at different growth stages by rice genotypes of different Zn densities. J. Agric. Food Chem. 2010, 58, 6767–6773. [Google Scholar] [CrossRef]
  82. Yadav, B.; Jogawat, A.; Lal, S.K.; Lakra, N.; Mehta, S.; Shabek, N.; Narayan, O.P. Plant mineral transport systems and the potential for crop improvement. Planta 2021, 253, 45. [Google Scholar] [CrossRef] [PubMed]
  83. Karley, A.J.; White, P.J. Moving cationic minerals to edible tissues: Potassium, magnesium, calcium. Curr. Opin. Plant Biol. 2009, 12, 291–298. [Google Scholar] [CrossRef] [PubMed]
  84. Ma, J.F.; Tsay, Y.-F. Transport Systems of Mineral Elements in Plants: Transporters, Regulation and Utilization. Plant Cell Physiol. 2021, 62, 539–540. [Google Scholar] [CrossRef] [PubMed]
  85. Tan, L.; Zhu, Y.; Fan, T.; Peng, C.; Wang, J.; Sun, L.; Chen, C. OsZIP7 functions in xylem loading in roots and inter-vascular transfer in nodes to deliver Zn/Cd to grain in rice. Biochem. Biophys. Res. Commun. 2019, 512, 112–118. [Google Scholar] [CrossRef] [PubMed]
  86. Grusak, M.A. Enhancing mineral content in plant food products. J. Am. Coll. Nutr. 2002, 21, 178S–183S. [Google Scholar] [CrossRef]
  87. Liang, C.; Tian, J.; Liao, H. Proteomics dissection of plant responses to mineral nutrient deficiency. Proteomics 2013, 13, 624–636. [Google Scholar] [CrossRef]
  88. Zhang, L.; Gao, C.; Chen, C.; Zhang, W.; Huang, X.Y.; Zhao, F.J. Overexpression of Rice OsHMA3 in Wheat Greatly Decreases Cadmium Accumulation in Wheat Grains. Environ. Sci. Technol. 2020, 54, 10100–10108. [Google Scholar] [CrossRef]
  89. Zlobin, I.E. Current understanding of plant zinc homeostasis regulation mechanisms. Plant Physiol. Biochem. 2021, 162, 327–335. [Google Scholar] [CrossRef]
  90. Zeng, H.; Wu, H.; Yan, F.; Yi, K.; Zhu, Y. Molecular regulation of zinc deficiency responses in plants. J. Plant Physiol. 2021, 261, 153419. [Google Scholar] [CrossRef]
  91. O’rourke, J.A.; Graham, M.A. Gene expression responses to sequential nutrient deficiency stresses in soybean. Int. J. Mol. Sci. 2021, 22, 1252. [Google Scholar] [CrossRef]
  92. Sasaki, A.; Yamaji, N.; Yokosho, K.; Ma, J.F. Nramp5 is a major transporter responsible for manganese and cadmium uptake in rice. Plant Cell 2012, 24, 2155–2167. [Google Scholar] [CrossRef] [PubMed]
  93. Satoh-Nagasawa, N.; Mori, M.; Nakazawa, N.; Kawamoto, T.; Nagato, Y.; Sakurai, K.; Takahashi, H.; Watanabe, A.; Akagi, H. Mutations in rice (Oryza sativa) heavy metal ATPase 2 (OsHMA2) restrict the translocation of zinc and cadmium. Plant Cell Physiol. 2012, 53, 213–224. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  94. Sasaki, A.; Yamaji, N.; Ma, J.F. Overexpression of OsHMA3 enhances Cd tolerance and expression of Zn transporter genes in rice. J. Exp. Bot. 2014, 65, 6013–6021. [Google Scholar] [CrossRef] [PubMed]
  95. Park, W.; Han, K.H.; Ahn, S.J. Differences in root-to-shoot cd and Zn translocation and by HMA3 and 4 could influence chlorophyll and anthocyanin content in arabidopsis ws and col-0 ecotypes under excess metals. Soil Sci. Plant Nutr. 2012, 58, 334–348. [Google Scholar] [CrossRef]
  96. Yan, J.; Wang, P.; Wang, P.; Yang, M.; Lian, X.; Tang, Z.; Huang, C.F.; Salt, D.E.; Zhao, F.J. A loss-of-function allele of OsHMA3 associated with high cadmium accumulation in shoots and grain of Japonica rice cultivars. Plant Cell Environ. 2016, 39, 1941–1954. [Google Scholar] [CrossRef]
  97. Zhang, C.; Ge, Y. Response of Glutathione and Glutathione S-transferase in Rice Seedlings Exposed to Cadmium Stress. Rice Sci. 2008, 15, 73–76. [Google Scholar] [CrossRef]
  98. Wang, K.; Yu, H.; Ye, D.; Wang, Y.; Zhang, X.; Huang, H.; Zheng, Z.; Li, T. The critical role of the shoot base in inhibiting cadmium transport from root to shoot in a cadmium-safe rice line (Oryza sativa L.). Sci. Total Environ. 2021, 765, 142710. [Google Scholar] [CrossRef]
  99. Huang, Y.; Chen, J.; Zhang, D.; Fang, B.; YangJin, T.; Zou, J.; Chen, Y.; Su, N.; Cui, J. Enhanced vacuole compartmentalization of cadmium in root cells contributes to glutathione-induced reduction of cadmium translocation from roots to shoots in pakchoi (Brassica chinensis L.). Ecotoxicol. Environ. Saf. 2021, 208, 111616. [Google Scholar] [CrossRef]
  100. Cao, Z.; Mou, R.; Cao, Z.; Lin, X.; Ma, Y.; Zhu, Z.; Chen, M. Quantitation of glutathione S-transferases in rice (Oryza sativa L.) roots exposed to cadmium by liquid chromatography-tandem mass spectrometry using isotope-labeled wing peptides as an internal standard. Plant Methods 2017, 13, 64. [Google Scholar] [CrossRef]
  101. Zhao, J.; Yang, W.; Zhang, S.; Yang, T.; Liu, Q.; Dong, J.; Fu, H.; Mao, X.; Liu, B. Genome-wide association study and candidate gene analysis of rice cadmium accumulation in grain in a diverse rice collection. Rice 2018, 11, 61. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Effect of zinc fertilization on cadmium concentration in vegetative tissues: Cadmium concentrations (mg/g DW) in vegetative tissues (root, stem, and lower leaves) of rice lines 546, 667, and 428 that were harvested at 4 developmental stages (vegetative (30 days old), anthesis (first day of flowering), 2 weeks after anthesis (2WAA), and 4 weeks after anthesis (4WAA)). Plants were treated with 1 µM Cd + 2 µM Zn (c1z2) or 1 µM Cd + 10 µM Zn (c1z10). Means ± se (n = 3). Variance was compared using 2-way ANOVA with a multiple comparison test and Tukey correction. Different colored letters (abc) indicate significant difference between the three lines tested at a particular treatment and reproductive stage (p < 0.05). No letter indicates no statistical significance.
Figure 1. Effect of zinc fertilization on cadmium concentration in vegetative tissues: Cadmium concentrations (mg/g DW) in vegetative tissues (root, stem, and lower leaves) of rice lines 546, 667, and 428 that were harvested at 4 developmental stages (vegetative (30 days old), anthesis (first day of flowering), 2 weeks after anthesis (2WAA), and 4 weeks after anthesis (4WAA)). Plants were treated with 1 µM Cd + 2 µM Zn (c1z2) or 1 µM Cd + 10 µM Zn (c1z10). Means ± se (n = 3). Variance was compared using 2-way ANOVA with a multiple comparison test and Tukey correction. Different colored letters (abc) indicate significant difference between the three lines tested at a particular treatment and reproductive stage (p < 0.05). No letter indicates no statistical significance.
Agronomy 12 02182 g001
Figure 2. Effect of zinc fertilization on cadmium concentration in reproductive tissues: Cadmium concentrations (mg/g DW) in reproductive tissues (flag leaf, peduncle, rachis, and grains) of rice lines 546, 667, and 428 that were harvested at 4 developmental stages (vegetative (30 days old), anthesis (first day of flowering), 2 weeks after anthesis (2WAA), and 4 weeks after anthesis (4WAA)). Plants were treated with 1 µM Cd + 2 µM Zn (c1z2) or 1 µM Cd + 10 µM Zn (c1z10). Means ± se (n = 3). Variance was compared using 2-way ANOVA with a multiple comparison test and Tukey correction. Different colored letters (abc) indicate significant difference between the three lines tested at a particular treatment and reproductive stage (p < 0.05). No letter indicates no statistical significance.
Figure 2. Effect of zinc fertilization on cadmium concentration in reproductive tissues: Cadmium concentrations (mg/g DW) in reproductive tissues (flag leaf, peduncle, rachis, and grains) of rice lines 546, 667, and 428 that were harvested at 4 developmental stages (vegetative (30 days old), anthesis (first day of flowering), 2 weeks after anthesis (2WAA), and 4 weeks after anthesis (4WAA)). Plants were treated with 1 µM Cd + 2 µM Zn (c1z2) or 1 µM Cd + 10 µM Zn (c1z10). Means ± se (n = 3). Variance was compared using 2-way ANOVA with a multiple comparison test and Tukey correction. Different colored letters (abc) indicate significant difference between the three lines tested at a particular treatment and reproductive stage (p < 0.05). No letter indicates no statistical significance.
Agronomy 12 02182 g002
Figure 3. Effect of cadmium and zinc fertilization on zinc concentration in vegetative tissues: Zinc concentrations (mg/g DW) in vegetative tissues (root, stem, and lower leaves) of rice lines 546, 667, and 428 that were harvested at 4 developmental stages (vegetative (30 days old), anthesis (first day of flowering), 2 weeks after anthesis (2WAA), and 4 weeks after anthesis (4WAA)). Plants were treated with 0 µM Cd + 2 µM Zn (c0z2), 1 µM Cd + 2 µM Zn (c1z2), 0 µM Cd + 10 µM Zn (c0z10), or 1 µM Cd + 10 µM Zn (c1z10). Means ± se (n = 3). Variance was compared using 2-way ANOVA with a multiple comparison test and Tukey correction. Different colored letters (abc) indicate significant difference between the three lines tested at a particular treatment and reproductive stage (p < 0.05). No letter indicates no statistical significance.
Figure 3. Effect of cadmium and zinc fertilization on zinc concentration in vegetative tissues: Zinc concentrations (mg/g DW) in vegetative tissues (root, stem, and lower leaves) of rice lines 546, 667, and 428 that were harvested at 4 developmental stages (vegetative (30 days old), anthesis (first day of flowering), 2 weeks after anthesis (2WAA), and 4 weeks after anthesis (4WAA)). Plants were treated with 0 µM Cd + 2 µM Zn (c0z2), 1 µM Cd + 2 µM Zn (c1z2), 0 µM Cd + 10 µM Zn (c0z10), or 1 µM Cd + 10 µM Zn (c1z10). Means ± se (n = 3). Variance was compared using 2-way ANOVA with a multiple comparison test and Tukey correction. Different colored letters (abc) indicate significant difference between the three lines tested at a particular treatment and reproductive stage (p < 0.05). No letter indicates no statistical significance.
Agronomy 12 02182 g003
Figure 4. Effect of cadmium and zinc fertilization on zinc concentration in reproductive tissues: Zinc concentrations (mg/g DW) in reproductive tissues (flag leaf, peduncle, rachis, and grains) of rice lines 546, 667, and 428 that were harvested at 4 developmental stages (vegetative (30 days old), anthesis (first day of flowering), 2 weeks after anthesis (2WAA), and 4 weeks after anthesis (4WAA). Plants were treated with 0 µM Cd + 2 µM Zn (c0z2), 1 µM Cd + 2 µM Zn (c1z2), 0 µM Cd + 10 µM Zn (c0z10), or 1 µM Cd + 10 µM Zn (c1z10). Means ± se (n = 3). Variance was compared using 2-way ANOVA with a multiple comparison test and Tukey correction. Different colored letters (abc) indicate significant difference between the three lines tested at a particular treatment and reproductive stage (p < 0.05). No letters indicate no statistical significance.
Figure 4. Effect of cadmium and zinc fertilization on zinc concentration in reproductive tissues: Zinc concentrations (mg/g DW) in reproductive tissues (flag leaf, peduncle, rachis, and grains) of rice lines 546, 667, and 428 that were harvested at 4 developmental stages (vegetative (30 days old), anthesis (first day of flowering), 2 weeks after anthesis (2WAA), and 4 weeks after anthesis (4WAA). Plants were treated with 0 µM Cd + 2 µM Zn (c0z2), 1 µM Cd + 2 µM Zn (c1z2), 0 µM Cd + 10 µM Zn (c0z10), or 1 µM Cd + 10 µM Zn (c1z10). Means ± se (n = 3). Variance was compared using 2-way ANOVA with a multiple comparison test and Tukey correction. Different colored letters (abc) indicate significant difference between the three lines tested at a particular treatment and reproductive stage (p < 0.05). No letters indicate no statistical significance.
Agronomy 12 02182 g004
Figure 5. Gene expression of transport proteins OsNramp5 (A), OsHMA2 (B), and OsHMA3 (C), in roots lines 546, 667, and 428 harvested at the vegetative reproductive stages: Quantitative real-time PCR was used to quantify mean relative gene expression (n = 3) using the delta-delta Ct method with standard error. Rice ubiquitin was used as the housekeeping gene. Two-way ANOVA multiple comparisons test with Tukey correction was used to analyze significant differences in gene expression. Statistical significance (p < 0.05) between the three lines at a particular treatment are annotated by (a, b, c) and differences between the four treatments within a rice line are annotated by (w, x, y, z) above bars.
Figure 5. Gene expression of transport proteins OsNramp5 (A), OsHMA2 (B), and OsHMA3 (C), in roots lines 546, 667, and 428 harvested at the vegetative reproductive stages: Quantitative real-time PCR was used to quantify mean relative gene expression (n = 3) using the delta-delta Ct method with standard error. Rice ubiquitin was used as the housekeeping gene. Two-way ANOVA multiple comparisons test with Tukey correction was used to analyze significant differences in gene expression. Statistical significance (p < 0.05) between the three lines at a particular treatment are annotated by (a, b, c) and differences between the four treatments within a rice line are annotated by (w, x, y, z) above bars.
Agronomy 12 02182 g005
Figure 6. Gene expression of transport proteins OsNramp5 (A), OsHMA2 (B), OsHMA3 (C), in roots and OsHMA3 (D) in flag leaves of lines 546, 667, and 428 harvested at the anthesis reproductive stages: Quantitative real-time PCR was used to quantify mean relative gene expression (n = 3) using the delta-delta Ct method with standard error. Rice ubiquitin was used as the housekeeping gene. Two-way ANOVA multiple comparisons test with Tukey correction was used to analyze significant differences in gene expression. Statistical significance (p < 0.05) between the three lines at a particular treatment are annotated by (a, b, c) and differences between the four treatments within a rice line are annotated by (w, x, y) above bars.
Figure 6. Gene expression of transport proteins OsNramp5 (A), OsHMA2 (B), OsHMA3 (C), in roots and OsHMA3 (D) in flag leaves of lines 546, 667, and 428 harvested at the anthesis reproductive stages: Quantitative real-time PCR was used to quantify mean relative gene expression (n = 3) using the delta-delta Ct method with standard error. Rice ubiquitin was used as the housekeeping gene. Two-way ANOVA multiple comparisons test with Tukey correction was used to analyze significant differences in gene expression. Statistical significance (p < 0.05) between the three lines at a particular treatment are annotated by (a, b, c) and differences between the four treatments within a rice line are annotated by (w, x, y) above bars.
Agronomy 12 02182 g006
Figure 7. Relative expression of non-protein thiols γ-Glutamyl cysteine ligase (GSH ligase) (A,C), GSHS (D), and Phytochelatin synthase (PCS1) (B,E) from roots (A,B) and flag leaves (CE) of lines 546, 667, and 428 harvested at anthesis stage: Quantitative real-time PCR was used to quantify mean relative gene expression (n = 3) using the delta-delta Ct method with standard error. Rice ubiquitin was used as the housekeeping gene. Two-way ANOVA multiple comparisons test with Tukey correction was used to analyze significant differences in gene expression. Statistical significance (p < 0.05) between the three lines at a particular treatment are annotated by (a, b, c) and differences between the four treatments within a rice line are annotated by (w, x, y) above bars.
Figure 7. Relative expression of non-protein thiols γ-Glutamyl cysteine ligase (GSH ligase) (A,C), GSHS (D), and Phytochelatin synthase (PCS1) (B,E) from roots (A,B) and flag leaves (CE) of lines 546, 667, and 428 harvested at anthesis stage: Quantitative real-time PCR was used to quantify mean relative gene expression (n = 3) using the delta-delta Ct method with standard error. Rice ubiquitin was used as the housekeeping gene. Two-way ANOVA multiple comparisons test with Tukey correction was used to analyze significant differences in gene expression. Statistical significance (p < 0.05) between the three lines at a particular treatment are annotated by (a, b, c) and differences between the four treatments within a rice line are annotated by (w, x, y) above bars.
Agronomy 12 02182 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tavarez, M.; Grusak, M.A.; Sankaran, R.P. Effects of Zinc Fertilization on Grain Cadmium Accumulation, Gene Expression, and Essential Mineral Partitioning in Rice. Agronomy 2022, 12, 2182. https://doi.org/10.3390/agronomy12092182

AMA Style

Tavarez M, Grusak MA, Sankaran RP. Effects of Zinc Fertilization on Grain Cadmium Accumulation, Gene Expression, and Essential Mineral Partitioning in Rice. Agronomy. 2022; 12(9):2182. https://doi.org/10.3390/agronomy12092182

Chicago/Turabian Style

Tavarez, Michael, Michael A. Grusak, and Renuka P. Sankaran. 2022. "Effects of Zinc Fertilization on Grain Cadmium Accumulation, Gene Expression, and Essential Mineral Partitioning in Rice" Agronomy 12, no. 9: 2182. https://doi.org/10.3390/agronomy12092182

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop