Toavintsoa Rajonandraina1†
Yoshiaki Ueda2*†
Matthias Wissuwa2,3
Guy J. D. Kirk4
Tovohery Rakotoson1
Hanna Manwaring4
Andry Andriamananjara1Tantely Razafimbelo1- 1Laboratoire des RadioIsotopes (LRI), Université d’Antananarivo, Antananarivo, Madagascar
- 2Crop, Livestock and Environment Division, Japan International Research Center for Agricultural Sciences (JIRCAS), Tsukuba, Japan
- 3PhenoRob Cluster & Institute of Crop Science and Resource Conservation (INRES), University of Bonn, Bonn, Germany
- 4School of Water, Energy and Environment, Cranfield University, Cranfield, United Kingdom
Introduction: Iron (Fe) toxicity is a widespread nutritional disorder in lowland rice causing growth retardation and leaf symptoms referred to as leaf bronzing. It is partly caused by an imbalance of nutrients other than Fe and supply of these is known to mitigate the toxicity. But the physiological and molecular mechanisms involved are unknown.
Methods: We investigated the effect of magnesium (Mg) on Fe toxicity tolerance in a field study in the Central Highlands of Madagascar and in hydroponic experiments with excess Fe (300 mg Fe L-1). An RNA-seq analysis was conducted in a hydroponic experiment to elucidate possible mechanisms underlying Mg effects.
Results and discussion: Addition of Mg consistently decreased leaf bronzing under both field and hydroponic conditions, whereas potassium (K) addition caused minor effects. Plants treated with Mg tended to have smaller shoot Fe concentrations in the field, suggesting enhanced exclusion at the whole-plant level. However, analysis of multiple genotypes showed that Fe toxicity symptoms were also mitigated without a concomitant decrease of Fe concentration, suggesting that increased Mg supply confers tolerance at the tissue level. The hydroponic experiments also suggested that Mg mitigated leaf bronzing without significantly decreasing Fe concentration or oxidative stress as assessed by the content of malondialdehyde, a biomarker for oxidative stress. An RNA-seq analysis revealed that Mg induced more changes in leaves than roots. Subsequent cis-element analysis suggested that NAC transcription factor binding sites were enriched in genes induced by Fe toxicity in leaves. Addition of Mg caused non-significant enrichment of the same binding sites, suggesting that NAC family proteins may mediate the effect of Mg. This study provides clues for mitigating Fe toxicity-induced leaf bronzing in rice.
1 Introduction
Rice is a staple crop for more than half the world’s population and its consumption is steadily increasing with population growth (Arouna et al., 2021). Demand is increasing particularly rapidly in Sub-Saharan Africa (SSA) where economic growth, urbanization and changes in lifestyle have increased people’s preference for rice (Wopereis et al., 2013). Despite increases in production in SSA since the 1960s (Arouna et al., 2021), yields and production are relatively low and more than 40% is imported accounting for a third of the global rice trade (FAO, 2017). Past production increase have been due to expansion of the cultivated area (Tsujimoto et al., 2019) though yields per unit area have also increased recently (Wopereis et al., 2013). Current mean yields (2.1 t ha-1) are far below potential yields under prevailing agro-climatic conditions (7.5–10.8 t ha-1) due to various abiotic and biotic constraints as well as a lack of appropriate technologies (Tsujimoto et al., 2019). Increasing productivity under adverse conditions and closing the gap between actual and potential yields are key to increasing rice productivity in SSA.
Major abiotic constraints to rice in SSA are the prevalence of highly weathered, nutrient depleted soils and the nutritional disorder iron (Fe) toxicity (Becker and Asch, 2005; van Oort, 2018). Iron toxicity occurs exclusively under submerged soil conditions where exclusion of oxygen favors generation of soluble ferrous iron (Fe(II)) from insoluble ferric iron (Fe(III)). The toxicity is exacerbated by low pH and low nutrient status, and hence it is widespread in the highly weathered soils that typify inland valleys in SSA which otherwise have great potential for rice (Becker and Asch, 2005; Kirk et al., 2022). Iron toxicity typically results in >15% rice yield loss and sometimes complete loss (Sahrawat, 2004; Audebert and Fofana, 2009). Although field management such as mid-season drainage and liming may reduce Fe(II) availability and mitigate Fe toxicity (Becker and Asch, 2005; Fageria et al., 2008), such measures are labor intensive and incur high costs. Hence, they are unlikely to be implemented in the resource-poor countries in SSA where Fe toxicity is particularly problematic. Development of tolerant germplasm is a more promising approach.
The presence of excess soluble Fe(II) in soil causes overloading of Fe in plants. Excess Fe disrupts the redox status, causes oxidative stress and foliar damage referred to as leaf bronzing, and retards growth (Aung and Masuda, 2020). Tolerance mechanisms include exclusion of excess Fe from shoots by restricting uptake into roots or by retention in roots, thus reducing root-to-shoot translocation. Genotypes lacking such Fe exclusion mechanisms may have different mechanisms to limit tissue damage despite excessive Fe concentrations (Engel et al., 2012; Matthus et al., 2015; Aung and Masuda, 2020). Such ‘includer’ mechanisms are further classified into ‘avoidance’, which involves compartmentalization of Fe in older tissues, and ‘tissue tolerance’, which involves storage of Fe in less bioactive states such as by sequestration in vacuoles and detoxification of reactive oxygen species (ROS) formed under excess Fe (Wu et al., 2017; Aung and Masuda, 2020). This suggests further genetic improvement of tolerance may be achieved by pyramiding genes for distinct tolerance mechanisms. Quantitative trait loci (QTLs) and candidate genes for Fe toxicity tolerance in rice have been identified (Dufey et al., 2009; Dufey et al., 2015a, Dufey et al., 2015b; Wu et al., 2014; Matthus et al., 2015; Diop et al., 2020; Melandri et al., 2021; Wairich et al., 2021). However, studies so far suggest that many small effect loci underlie each of the above mechanisms, as well as overall Fe toxicity tolerance.
A major problem in understanding tolerance mechanism under field conditions is interactions between Fe toxicity and the nutrient deficiencies that are prevalent in much of SSA, particularly deficiencies of phosphorus (P), potassium (K), calcium (Ca) and magnesium (Mg) (Howeler, 1973; Li et al., 2001; Gao et al., 2014; Suriyagoda et al., 2017). Of these, interaction with K uptake have been most studied. A mutant-based study (Wu et al., 2019) showed that K uptake through the K channel AKT1 restricts excess Fe uptake and mitigates excess Fe-induced leaf bronzing, and a genome-wide association study (Matthus et al., 2015) has also linked AKT1 to genotypic differences in shoot Fe concentration. Evidently K uptake reinforces Fe exclusion mechanisms. However, the interactions of Fe toxicity with other nutrients, such as Mg and Ca, have been less intensively investigated, and genes that potentially mediate the effect of these elements remain to be elucidated.
A recent meta-analysis showed that sub-optimal tissue Mg concentrations are particularly common in rice at Fe toxic sites (Kirk et al., 2022). Highly weathered soils in humid tropical regions are generally depleted in Mg due to high rates of leaching (Sun et al., 2013) and little Mg is added in fertilizers. Root-induced changes in the rhizosphere under Fe toxicity, particularly Fe oxidation and resulting acidification, will tend to decrease Mg solubility and hence its availability for plant uptake (Kirk et al., 2022). Deficiency of Mg tends to increase the levels of toxic reactive oxygen species (ROS) (Hauer-Jákli & Tränkner, 2019) typically generated under Fe toxicity (Wu et al., 2017). Links between Mg status and Fe toxicity response have also been inferred by transcriptome analysis (Kobayashi et al., 2018). This showed Mg deficiency upregulated excess Fe-inducible OsFER2, which encodes a Fe storage protein, and downregulated Fe sufficiency-suppressive OsMIR and OsIRO2, which encode a mitochondrial protein regulating Fe homeostasis and a key transcriptional activator regulating Fe deficiency responses, respectively (Ogo et al., 2007; Ishimaru et al., 2009; Stein et al., 2009; Aung et al., 2018). These pieces of evidence suggest that some key components of the Fe toxicity tolerance response may be inactivated at sub-optimal leaf Mg content. However, our knowledge regarding the effect of Mg supply on the Fe toxicity response and the molecular mechanisms involved remains very incomplete.
Our objectives were first, to characterize the effects of supplying Mg and K on Fe toxicity responses in a diverse set of genotypes at a strongly Fe-toxic field site in the Central Highlands of Madagascar; and second, to dissect possible mechanisms underlying the effects of Mg supply in hydroponic experiments with transcriptomic analyses.
2 Materials and methods
2.1 Field study
A field experiment was conducted with seven rice genotypes differing in their tolerance mechanism to Fe toxicity (Table 1; Rajonandraina et al., 2023). Plants were grown under irrigated conditions at Sambaina, Manjakandriana district in the Central Highlands of Madagascar (18°53’ S, 47°47’ E) during the wet season (December to May) in 2020-2021. The site’s soil is a Gleysol with clay loam texture, pH (H2O) 4.5, pH (KCl) 4.1, organic carbon 62 g kg-1, cation exchange capacity 2.4 cmolc kg-1 and total Fe 57 g kg-1 of which 10 g kg-1 is easily-soluble on soil reduction (Rakotoson et al., 2019).
Table 1 List of rice genotypes used in the field study.
A split-plot design was implemented with eight fertilizer treatments containing all genotypes and four replications. The following treatments were applied: control (i.e. no fertilizer), Mg, NK, NKMg, NP, NPK, NPMg and NPKMg at transplanting at the following dose as in a previous study carried out in the same field (Rajonandraina et al., 2023); 50 kg N ha-1 as urea, 20 kg P ha-1 as triple superphosphate, 20 kg K ha-1 as potassium sulphate, and 26 kg Mg ha-1 as kieserite. The size of a subplot, which contained one genotype, was 1.6 m2 (0.8 m × 2 m). Subplots were randomized in each plot, which was further randomized in blocks of replicates. Four subplots were prepared for each genotype and treatment. Seedlings were transplanted at 21-d old with single plant per hill and a spacing of 20 cm × 20 cm between hills.
2.2 Sampling and measurements-field study
Sampling and measurement were done at different stages of plant growth, namely, tillering, booting, flowering and maturity. Leaf symptoms induced by Fe toxicity were visually scored as leaf bronzing score (LBS) across the whole plant canopy of each subplot as a percentage of leaf area affected on a scale from 0 (no symptoms) to 10 (100% of the leaf area affected) according to Wu et al. (2014). LBS was scored for each time point except at the maturity because of difficulties in differentiation with plant senescence. Plants were collected from two randomly selected hills per subplot, and each plant was dissected into youngest leaf (YL), middle leaves (ML), old leaves (OL) and stem/leaf sheath (ST). Flag leaves (FL) and panicles were additionally sampled at the flowering and maturing stage, respectively. The tissues were oven-dried at 60°C for 48 h until the samples were completely dry and weighed. The samples from two plants from a subplot were pooled prior to element analyses and considered as one replicate. For tissue Fe analysis, oven-dried samples were ground to a powder using a mixer mill (Retsch ZM 200, 0.2-mm sieve).
2.3 Hydroponic study
2.3.1 Experiment 1
Seeds of Nipponbare were sterilized with 5% NaClO for 5 min and rinsed with tap water. Sterilized seeds were placed on a moistened tissue paper on a petri dish and germinated at 28°C in darkness. After 3 d, germinated seeds were placed on top of a stainless mesh tray and grown with deionized water in darkness for 2 d. Subsequently, the seedlings were exposed to natural light and kept with 0.1 mM CaCl2 and 12 µM FeNa-EDTA solution. After further 12 d, seedlings were transplanted to 13-L hydroponic tanks containing 0.5x concentration of Yoshida nutrient solution (Yoshida et al., 1976). After 7 d, the solution was exchanged to 1x concentration of Yoshida solution. After further 7 d, plants were transferred to hydroponic tanks containing control solution (1x Yoshida solution; Fe is provided as 36 µM FeNa-EDTA [=2 mg Fe L-1]) or Fe excess solution (control solution + 300 mg Fe L-1 as FeSO4.7H2O). A low concentration (0.1%) of melted agar (016-15817, Wako Fujifilm) was added to each container (Wang et al., 2008). The full strength (1x) Yoshida solution contained 1.42 mM NH4NO3, 0.1 mM NaH2PO4, 0.5 mM K2SO4, 1 mM CaCl2, 1 mM MgSO4, 36 µM FeNa-EDTA, 9 µM MnCl2, 18.5 µM H3BO3, 0.16 µM CuSO4, 1.5 µM ZnSO4 and 0.07 µM (NH4)6Mo7O24. For combinatorial treatment with different concentrations of Mg and Ca, MgSO4 and CaSO4 was added at different ratio to achieve Ca2+/Mg2+ (mM) = 1/0.2, 1/1, 1/5, and 5/1, together with normal or excess amount of Fe. After 10 d of treatment, aboveground tissues were separated into young leaves (2 most recently fully expanded leaves and emerging leaf) and the rest of tissues (old leaves + stem/leaf sheath) and used for biochemical and element analysis. LBS was measured from the 3 most recently fully expanded leaves from 5 plants in each treatment. Supplemental LED light was installed in the greenhouse due to low solar irradiation during the experiment.
2.3.2 Experiment 2
Plants were grown and treated in the same manner as in Experiment 1. Plants were treated with control, excess Fe (control + 300 mg L-1 Fe) and excess Fe and Mg (control + 4 mM MgSO4 + 300 mg L-1 Fe) conditions. Samples for RNA extraction (n=3) were harvested after 10 d from the onset of the treatment either from the middle part of the top fully expanded leaf or whole roots. The samples were flash-frozen with liquid nitrogen and kept at -80°C until the analysis.
2.4 Mineral analysis
Fe concentration in dry tissue was determined as described previously (Hartmann and Asch, 2018). Briefly, Fe was extracted from pulverized dry tissues using 500 mM sodium dithionite solution, followed by colorimetric assay using 2,2’-dipyridyl. The concentration of Fe in the extract was determined by using a standard curve made with serially diluted Fe standard solution. For the measurement of other mineral elements (Ca, K, Mg), dry samples were acid-digested as reported previously (Ueda and Wissuwa, 2022; Rajonandraina et al., 2023) and measured by an inductively coupled plasma mass spectrometry (NexION 350, PerkinElmer) and inductively coupled plasma atomic emission spectrometer (ICPE-9000, Shimadzu) for field and hydroponic samples, respectively, using a standard curve made with serially diluted standard solution.
2.5 MDA analysis
The content of malondialdehyde (MDA) was quantified according to a previous report (Hodges et al., 1999). Briefly, 0.1% trichloroacetic acid was added to pulverized frozen samples and mixed vigorously. After centrifuging the mixture at 14,000 x g for 10 min at 4°C, 90 µL of the resultant supernatant was mixed either with solution I (0.01% 2,6-di-tert-butyl-4-methylphenol in 20% trichloroacetic acid) or solution II (solution I added with 0.65% 2-thiobarbituric acid) and heated at 95°C for 30 min. The resultant mixture was briefly centrifuged and the absorbance of the supernatant was recorded at 440, 532 and 600 nm using a microplate reader (Multiscan GO, Thermo Scientific).
2.6 RNA extraction and RT-PCR
Total RNA was extracted from frozen pulverized samples using RNeasy Plant Mini Kit (Qiagen) according to manufacturer’s instructions. cDNA was synthesized using the PrimeScript RT Master Mix (Perfect Real Time) kit (Takara). Quantitative RT-PCR was conducted using the CFX96 Real-Time PCR Detection System (Bio-Rad) and TB Green Premix Ex Taq II (Tli RNaseH Plus) (Takara). The expression of target genes was normalized against an internal reference gene OsC3H38 (Höller et al., 2015).
2.7 RNA sequencing
Sequence library was prepared with HiSeq Standard mRNA Library Prep kit using total RNA. The resultant library was sequenced with the HiSeqX ten instrument (Illumina) and 150 bp paired-end reads were generated. Quality of the raw reads were analyzed by the FastQC software (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Reads were trimmed by trimmomatic software (Bolger et al., 2014) using the following options; LEADING:3 TRAILING:3 SLIDINGWINDOW:4:20 MINLEN:25. Gene expression was quantified as reported previously (Pertea et al., 2016). Briefly, obtained reads were mapped to Nipponbare reference transcript (IRGSP-1.0) (Sakai et al., 2013) using HISAT2 software (Kim et al., 2015) with –min-intronlen 20 –max-intronlen 10000 options. Expression of genes was quantified using the Stringtie software (Pertea et al., 2015), yielding transcript per million (TPM) values as well as raw count data.
Differentially expressed genes were analyzed by the DESeq2 software (Love et al., 2014) using count data. PCA plot was created by plotPCA function in DESeq2 after transforming the data by the vst function.
2.8 Cis-element enrichment analysis
Cis-element enrichment analysis was performed as reported previously (Ashrafuzzaman et al., 2020). The promoter element (defined as upstream 1,000-bp from the start codon) for genes significantly affected by excess Fe or combination of excess Fe and Mg was obtained at Ensembl Plants Biomart website (Kinsella et al., 2011). The enrichment analysis was performed by the AME software (McLeay and Bailey, 2010), using the database of Arabidopsis cis-elements (O’Malley et al., 2016) as the reference. The promoter sequences of all the 18,960 genes expressed in leaves (i.e. average TPM >2) were used as control.
2.9 Data analysis
Statistical analyses were performed with the R program (Version 4.2.0; https://www.R-project.org/) with significance level set at P = 0.05. In the field experiment, treatment effects were analyzed with a mixed model ANOVA using the lmer function in the lme4 package (Bates et al., 2015). When evaluating the effect of treatment on all genotypes, treatment was considered as a fixed effect, while genotype, block (replicate), and the interaction between block and treatment factor were considered as random effects, whereas growth stage replaced genotype as random effect when evaluating the treatment effect on each genotype across growth stages. To analyze the effect of Mg, samples from NPK and NP treatments were pooled and considered as -Mg sample group, and samples from NPKMg and NPMg treatments were pooled and considered as +Mg sample group. Similarly, samples from NP and NPMg were pooled and considered as -K sample group, and samples from NPK and NPKMg were pooled and considered as +K sample group. Each of NP, NPMg, NPK, and NPKMg treatment contained 4 replicates of 7 genotypes, thus the sample size was 56 for each of -Mg, +Mg, -K and +K treatment.
In the hydroponic experiment, one-way ANOVA was used to analyze differences among different treatments. Subsequently, Tukey-Kramer post-hoc test was implemented to compare values of each treatment.
3 Results
3.1 Field study
3.1.1 Effects of K and Mg on rice traits
We evaluated the effects of combinatorial application of N, P, K, and Mg on shoot dry weight (SDW), Fe concentration, LBS and yield. Data for all seven genotypes and four growth stages showed that P treatment had a dominant positive effect on SDW and yield, suggesting that the growth of plants was greatly limited by P availability (Table 2; Figure S1). Thus, we only considered data from the plots supplied with P in subsequent analyses.
Table 2 ANOVA results on the effect of genotype, P, Mg and K on Fe toxicity-related traits.
The effects of K and Mg treatments on leaf symptoms were analyzed at different growth stages. The leaf symptoms as assessed by LBS significantly decreased with the application of Mg fertilizer by an average of 18% (Figure 1A). Although the decrease in LBS caused by Mg treatment lessened as the plants grew, the treatment had a consistent suppressive effect on LBS throughout the growth (34 and 6% in tillering and flowering, respectively). However, supply of Mg had no significant effect on SDW in the reproductive stages (i.e. flowering and maturity), resulting in non-significant effect on grain yield (Figures 1B, C; Table 2). The supply of K had a less striking effect on LBS at each growth stage, and the effect throughout the season was not significant (Figure 1D; Table 2). However, K addition significantly increased SDW during the reproductive stage by an average of 17% (Figure 1E) and this effect carried through to grain yield which was also significantly increased (+15%) (Figure 1F; Table 2). These results indicate that increased Mg supply has a beneficial effect on LBS but not on SDW or grain yield, while increased K supply had the opposite effect.
Figure 1 Effects of Mg and K on leaf bronzing, shoot dry weight and grain yield. (A, D) leaf bronzing score; (B, E) shoot DW; (C, F) grain yield for Mg treatment (A-C) and K treatment (D-F). Data are means ± standard errors (n = 56). The result of ANOVA is indicated as follows; ns, P > 0.05; *, P < 0.05; **, P < 0.01; ***, P < 0.001.
3.1.2 Mg effect on shoot Fe uptake and Fe concentrations in different tissues
To get insights into the mechanism of Mg-mediated alleviation of leaf bronzing, we analyzed shoot Fe concentration and content at different growth stages. The supply of Mg significantly increased shoot Fe concentration during the vegetative stages, but the opposite effect was observed during reproductive stages (Figure 2A). Consistent with the effect of Mg supply on SDW and Fe concentration, aboveground Fe content (i.e. the product of SDW and Fe concentration) increased in the +Mg treatment during the vegetative and booting stages (on average by 22%), but at the reproductive growth stage the Fe content was 13% lower in the +Mg treatment (Figure 2B). We further assessed Fe concentrations in different plant parts. Leaf Fe concentrations increased strongly between vegetative and reproductive stages, but the +Mg treatment consistently reduced foliar Fe accumulation, except for old leaves during the vegetative stages (Figure 2C).
Figure 2 Effects of Mg on shoot Fe concentration and Fe content. (A) Whole shoot Fe concentration, (B) whole shoot Fe content, and (C) Fe concentrations in above-ground tissues; from tillering to maturity. Data are means ± standard errors (n = 56). The result of ANOVA is indicated as follows; ns, P > 0.05; *, P < 0.05; **, P < 0.01; ***, P < 0.001.
We investigated if the response to Mg differed between rice genotypes depending on whether they were classified as Fe includers or excluders (Rajonandraina et al., 2023). The analysis revealed that the effect of Mg supply on LBS was genotype-specific and that the +Mg treatment alleviated leaf bronzing in four of the seven genotypes studied, but three includer-type genotypes did not benefit from Mg addition (Figure 3A). Curiously, these 3 includer-type genotypes had substantially reduced Fe concentrations in young leaves (26% decrease on average), whereas the reduction in Fe concentrations of young leaves was less pronounced in excluder genotypes (average decrease below 18%) (Figure 3B). Local rice variety X265, classified as an Fe includer, appears to be an exception, as Mg addition did not lower Fe concentrations in young leaves while significantly decreasing LBS. We analyzed if these changes were due to an increase in tissue Mg concentrations. Mg concentrations in young leaves at the booting stage were not significantly affected by Mg addition in any of the examined genotypes (Figure S2), raising the possibility that the effect of Mg was likely through altered rhizosphere processes or interaction with other nutrients (Figure S2), at least at this growth stage.
Figure 3 Effects of Mg on LBS and young leaves Fe concentrations in different genotypes at different growth stages. (A) leaf bronzing score (LBS) from tillering to flowering stage. (B) Fe concentration in young leaves from tillering to maturity stage. Data are means ± standard errors (n = 8). P values indicate the result of ANOVA analyzing the treatment effect throughout the growth stages.
3.2 Hydroponic study
3.2.1 Mg, but not Ca, alleviates Fe toxicity in hydroponic conditions
We conducted hydroponic studies to reveal potential effects of Mg and reveal associated molecular mechanisms. Plants received normal or excess Fe (2 and 300 mg Fe L-1, respectively), together with different concentrations of Mg. Root morphology and iron plaque formation were not substantially affected by the different Mg concentrations (Figure S3). We also note that additional Mg did not substantially affect the pH or the redox potential of the solution (Figure S4), suggesting that Mg treatment did not affect Fe availability in the hydroponic solution. Supplying additional Mg reduced LBS in the three newest fully expanded leaves, while Mg deficiency had an opposite effect (Figures 4A, B). On the other hand, increasing the concentration of Ca, which like Mg is predominantly present as a divalent cation in the hydroponic solution, did not alleviate leaf bronzing (Figure 4B). The content of malondialdehyde (MDA), which is a biomarker for oxidative stress, was measured in new leaves, which we defined as the emerging and two newest fully expanded leaves. The MDA content increased drastically under Fe toxicity but ANOVA indicated differences between Mg or Ca treatments were not significant (Figure 4C).
Figure 4 Effects of Mg and Ca in hydroponics. (A) Representative image of top fully expanded leaf after 10 d of excess Fe treatment in the presence of different concentrations of Mg. The scale bar indicates 2 cm. The below panels indicate magnified image. (B, C) LBS (B) and MDA content (C) are shown for control (2 mg L-1 Fe) and excess Fe (300 mg L-1 Fe) treatments. In B and C, one-way ANOVA was conducted for each treatment and the resultant P value is shown on the top of each graph. Different alphabets indicate that the values are significantly different. Data are means ± standard deviations (n=4-5).
Element concentrations in shoot tissues were substantially affected by different treatments. Fe concentrations in new leaves were not significantly affected by Ca or Mg treatments (Figure 5A) but in old leaves and stem/leaf sheath the addition of Mg caused a reduction in tissue Fe concentrations whereas the Ca treatment did not have a substantial effect (Figure 5B). Excess Fe drastically reduced tissue Mg concentrations but different levels of Mg significantly affected Mg concentrations both in new leaves and the rest of shoot (i.e. old leaves and stem/leaf sheaths), irrespective of Fe conditions (Figures 5C, D). Ca concentrations were also reduced by excess Fe but drastically increased by additional Ca supply. Increment in Mg supply negatively affected Ca concentrations under normal Fe condition but not under excess Fe (Figures 5E, F). Concentrations of K was rather stable and affected only to a small extent by different levels of Fe, Mg and Ca treatments (Figures 5G, H).
Figure 5 Effects of Fe, Mg and Ca on shoot element concentrations. Concentrations of Fe (A, B), Mg (C, D), Ca (E, F) and K (G, H) treated with different levels of Fe, Mg and Ca. Data for new leaves (A, C, E, G) and old leaves and stem/leaf sheath (B, D, F, H) are shown. One-way ANOVA was conducted for each treatment and the resultant P value is shown on the top of each graph. Different alphabets indicate that the values are significantly different. Data are means ± standard deviations (n=4-5).
3.2.2 Genes and cis-elements involved in Mg-mediated stress alleviation
We performed an RNA-seq analysis using root and leaf tissue of plants grown under the Cont, ++Fe and ++Fe++Mg conditions (Table S1). The PCA analysis clearly separated clusters for Cont, ++Fe and ++Fe++Mg for leaf samples (Figure 6A) but the separation between ++Fe and ++Fe++Mg sample groups was less clear for roots (Figure 6B). Considering that Mg alleviates Fe toxicity, we reasoned that the key genes involved in Mg-mediated effects would exhibit one of the following expression patterns; Cont < ++Fe++Mg < ++Fe or Cont > ++Fe++Mg > ++Fe. Thus, genes that fulfill these criteria in leaf or root samples were identified. Among 18,960 genes expressed in leaves (i.e. average TPM > 2), 219 and 171 genes were down- and up-regulated in ++Fe++Mg compared with ++Fe, respectively (Figure 6C). Among them, 70 and 39 genes exhibited the expression pattern of Cont < ++Fe++Mg < ++Fe and Cont > ++Fe++Mg > ++Fe, respectively (Figure 6C; Table S2). In roots, among 20,375 expressed genes, only a very small number of genes exhibited the pattern of Cont < ++Fe++Mg < ++Fe and Cont > ++Fe++Mg > ++Fe (9 and 0 genes, respectively) (Figure 6D; Table S2). We conclude that the function of Mg to mitigate Fe toxicity is via mechanisms specific to shoot rather than a root tissue.
Figure 6 Effects of Mg on gene expression. (A, B) The result of PCA analysis using gene expression data in leaf (A) and root (B). (C, D) Venn diagrams showing the number of differentially expressed genes in each pair of comparison in leaf (C) and root (D). (E, F) Venn diagrams showing the number of differentially expressed genes and previously suggested ROS-responsive genes. Fisher’s exact test was performed, and P value is indicated for significantly overlapped regions.
To examine if Mg treatment mitigates Fe-induced oxidative stress, we further compared our RNA-seq results with two previously reported transcriptome studies conducted under ozone stress, which causes oxidative stress directly in leaves (Frei et al., 2010; Ashrafuzzaman et al., 2018). 15,870 genes were detected in both previous studies and the current study, and among these, 647 genes were defined as ‘ROS-responsive’ because they were similarly affected by ozone stress in both previous studies (Frei et al., 2010; Ashrafuzzaman et al., 2018). These ROS-responsive genes overlapped with 207 ++Fe-induced and 10 ++Fe-suppressed genes, and the degree of overlap between ROS-responsive and ++Fe-inducible genes was highly significant (P = 1.5 × 10-59 by Fisher’s exact test), confirming that excess Fe causes oxidative stress in leaves (Figure 6E). On the other hand, only 6 ROS-responsive genes (0.9%) overlapped with genes whose expression follows Cont < ++Fe++Mg < ++Fe in leaves, and the significance level of the overlap was much lower (P = 0.028) (Figure 6F).
We next searched for potential regulatory factors that mediate the effect of Mg in leaves. We identified promoter elements that are enriched among differentially expressed genes in leaves under ++Fe or ++Fe++Mg conditions. 52 cis-elements were commonly enriched by both treatments (i.e. among 1,942 and 1,498 genes upregulated by ++Fe and ++Fe++Mg treatment, respectively), such as those for WRKY and homeobox transcription factors (Figure 7A; Table S3). Two homeobox binding motifs and two NAC binding sites were enriched exclusively in ++Fe++Mg and ++Fe treatments, respectively. Among 1,108 and 374 genes downregulated by ++Fe and ++Fe++Mg treatments, only one homeobox binding site was enriched in both treatments, and 5 cis-elements, including the binding site for NAC, bHLH and C2H2 transcription factors, were enriched only by ++Fe treatment, while no element was specifically detected by ++Fe++Mg treatment (Figure 7B; Table S3).
Figure 7 Cis-element enrichment analysis. The number of enriched cis-elements from the genes upregulated (A) and downregulated (B) genes by ++Fe and ++Fe++Mg treatments.
4 Discussion
4.1 Fe exclusion mechanisms conferred by Mg supply
Plant responses to an imbalance in the availability of a nutrient are dependent on the status of other mineral elements. For instance, the strength of P deficiency responses is affected by the availability of N and Fe (Balzergue et al., 2017; Hu et al., 2019; Ueda et al., 2020). Fe toxicity is observed at tissue Fe concentrations ranging from 300 to 2000 ppm (Dobermann and Fairhurst, 2000) and this wide range of critical tissue concentrations is likely due to interactions between excess Fe and deficiencies of other nutrients such as P, K, Ca, Zn and Mg. Such deficiencies may be exacerbated by excess Fe-induced changes in the rhizosphere coupled with unfavorable soil properties such as low cation exchange capacity and pH buffer power (Becker and Asch, 2005; Kirk et al., 2022). Increased availability of these elements consistently mitigates symptoms caused by Fe toxicity (Benckiser et al., 1984; Sahrawat, 2004). For instance, an application of dolomite, which is a mixture of Mg and Ca, effectively mitigated Fe toxicity-triggered reductions in growth and yield of rice in an acid Fe toxic soil in Sri Lanka (Suriyagoda et al., 2017). Similarly, the ratio of available Fe to other divalent cations (i.e. Ca, Mg, Mn) were more closely associated with the degree of yield reduction than Fe concentration alone in Fe toxic acid sulphate soils in Thailand (Moore et al., 1990). These studies would suggest the benefit of increased supply of divalent cations lies in lowering excess Fe uptake, possibly by direct competition, or increasing their supply offsets any negative effects of their deficiency on Fe toxicity responses. Our data suggest that in some genotypes Mg supply helps Fe exclusion mechanisms (Figure 3), which could be explained by multiple factors. Tadano (1975) reported that reduced Mg supply lowered the Fe exclusion capacity of roots and increased the ratio of Fe translocated to shoots in culture solution, suggesting that Mg decreases Fe uptake or excess root-to-shoot translocation or both. Additionally, increased Fe concentration in stems by Mg treatment (Figure 2C) may form another layer of Fe-exclusion mechanism, as suggested previously (Engel et al., 2012). These effects of Mg on Fe uptake, root-to-shoot translocation and retention in stems could be regarded as a general function of Mg under Fe toxicity, as seen by a negative correlation between leaf Fe and Mg concentrations in a pot study using Fe toxic soils (Genon et al., 1994).
4.2 Mg-induced tissue-tolerance mechanisms
Even though Fe exclusion mechanisms are strengthened, some pieces of evidence suggested tissue-tolerance mechanisms conferred by Mg supply. It was suggested that the concentration ratio of Fe to Mg + Ca in flag leaves, rather than the concentration of Fe per se, is a good indicator for leaf bronzing under Fe toxicity (Sylla, 1994). This implies that Ca and Mg may confer at least some tissue tolerance. Our observations in field and hydroponic conditions both indicated attenuated leaf bronzing formation by Mg supply without a concomitant reduction in Fe concentration (X265 in Figures 3–5). A previous study also showed that Mg deficiency did not have a substantial effect on leaf Fe concentration but did induce genetic responses similar to those observed under Fe toxicity (Kobayashi et al., 2018). This suggests that Mg supply may also support tissue-based tolerance mechanisms. The hypothesis that increased Mg supply confers tissue tolerance is also supported by genotypic differences in the benefit provided by Mg treatment. Our data show that the includer genotypes benefitted to a lesser extent from Mg fertilization than the excluder genotypes (Figure 3A). This implies that Mg supply may provide tissue tolerance to plants that already possess exclusion mechanisms, i.e., an additive effect. However, the associated mechanisms underlying increased tissue tolerance conferred by Mg supply might be different in field and hydroponic studies, since we did not observe significant increase in Mg concentrations in young leaves at the booting stage in the field study. It implies that altered rhizospheric chemical processes or secondary effects on the uptake of other nutrients (Figure S2), rather than an increase in foliar Mg concentrations, could explain this phenomenon, at least at this growth stage.
4.3 Possible physiological and genetic factors underlying Mg-induced tissue-tolerance mechanisms
Contrary to the field study, additional Mg supply increased tissue Mg concentrations in hydroponic study (Figure 5), which clearly suggests that tissue tolerance is conferred by increased tissue Mg concentrations. Deficiency of Mg disrupts Fe storage in leaf cell vacuoles (Kobayashi et al., 2018) and additional Mg supply might increase storage of Fe in vacuoles that is otherwise kept in cytosol or other compartments in a toxic state. Interestingly, in our hydroponic experiment, the reduction in leaf bronzing was not accompanied by a decrease of MDA, a biomarker for oxidative stress. These data indicate that the leaves of plants supplied with excess Mg perceive a similar extent of oxidative stress but form less leaf bronzing, suggesting that Mg specifically suppresses the process of cell death caused by oxidative stress. Thus, the underlying mechanism of Mg effects observed in our study is not via reduced oxidative stress as previously suggested (Tränkner et al., 2016; Hauer-Jákli and Tränkner, 2019). This is further supported by the gene expression data showing that most of the genes affected by Mg were not related to oxidative stress (Figures 6E, F).
The promoters of Fe-inducible genes were enriched with cis-elements such as WRKY, NAC and Homeobox-binding motifs as previously suggested (Kakei et al., 2021) (Figure 7A; Table S3). The binding sites for WRKY transcription factors were also enriched among the genes that were induced by the ++Fe++Mg treatment, suggesting that WRKY family proteins may not play pivotal roles in Mg-dependent alleviation of symptoms (Figure 7A; Table S3). In Arabidopsis, half of Mg deficiency-responsive genes are also responsive to ABA (Tanoi and Kobayashi, 2015), suggesting a close link between ABA and Mg status. Interestingly, some members of NAC transcription factor family that likely mediate the effect of Mg are induced by ABA and play central roles in the ABA signal transduction pathway (Todaka et al., 2012). Our results also show that many Mg-responsive genes (i.e. genes that exhibit the expression patterns of Cont <++Fe++Mg < ++Fe or Cont > ++Fe++Mg > ++Fe) are ABA, salinity or drought-related genes; Among 70 genes whose expression pattern follows Cont < ++Fe++Mg < ++Fe in leaf, OsEBP89, bHLH035, ONAC60, MYBR1, PIP1;3, and RAV2 are known to be involved in drought or salinity stress response and tolerance (Fang et al., 2014; Duan et al., 2016; Yin et al., 2017; Chen et al., 2018; Liu et al., 2020; Zhang et al., 2020). These genes also contained some NAC transcription factors (i.e. NAC15, ONAC60, NAC58, NAC121 and NAC103). Further investigations are necessary to identify the key gene(s) mediating the effect of Mg and elucidate the relationship between Mg, ABA and the causal gene under Fe toxicity conditions.
4.4 Beneficial effects of Mg on mineral stresses
Optimal Mg supply is known to counteract several soil-related stresses. For instance, Mg confers tolerance to aluminum (Al) stress in various plant species (Rengel et al., 2015). In rice, Mg alleviates Al stress partly by reducing Al uptake and oxidative stress (Pandey et al., 2013), although the physiological pathways affected by Mg may vary in different plant species (Silva et al., 2001a). Mg also confers tolerance of salinity stress in rice (Chen et al., 2017). In both cases, Mg uptake, rather than the contact of roots with high Mg concentration, is the key factor, as shown by the attenuated tolerance in mutant lines lacking functional MGT1 which is a plasmamembrane-localized major Mg transporter (Chen et al., 2012; Chen et al., 2017).
Contrary to the effect of Mg, additional Ca application in our hydroponic study did not significantly affect leaf bronzing or Fe accumulation. This could be partly due to the similarity of the ionic radii of Fe2+ (77 pm) and Mg2+ (72 pm) compared with Ca2+ (100 pm). Similarity in radii may mean that Mg2+ competes with Fe2+ better than Ca2+ for binding with other molecules and mitigate the deleterious effects of Fe to confer tissue tolerance. Similar Mg-specific effects were also observed under Al stress. A previous investigation in soybean suggested that Al toxicity was specifically mitigated by Mg2+ but not by similar concentration of Ca2+ (Silva et al., 2001a; Silva et al., 2001b), consistent with our observation that the positive effect of Mg2+ is not simply due to the ionic strengths of competing divalent cations. Interestingly, the expression of MGT1 is induced both in roots and leaves under Fe toxicity (Table S1). Thus, it is possible that upregulation of MGT1 and accompanying increase of Mg uptake might be an intrinsic mechanism that plants possess to cope with Fe toxicity, as well as other edaphic stresses such as salinity and Al toxicity. MGT1 is found within a QTL explaining Al stress tolerance (Yamaji et al., 2009). Intriguingly, a QTL for Fe toxicity-induced leaf bronzing (qFETOX1-2) was previously found close to the position of MGT1 (37.7 Mb on chromosome 1) (Wu et al., 2014), which is also close to a QTL for foliar Mg concentration (Norton et al., 2010). Whether allelic variation exists at MGT1 locus, and if that is the cause of genotypic variation of leaf bronzing, needs further investigation.
4.5 Conclusions and future perspectives
We have shown in field and hydroponic experiments that increasing the supply of Mg to rice roots improves Fe toxicity tolerance via exclusion- and tissue tolerance-based mechanisms. Evidence from transcriptome analysis indicated the involvement of several genes expressed in shoot tissue that potentially mediate the Mg effect, and of particular relevance seem to be the altered expression levels of genes containing the binding sites for NAC transcription factors. Previous studies showed more than three-fold variation in leaf Mg concentration among diverse rice accessions grown under control conditions (Chen and Ma, 2013; Yang et al., 2018), suggesting scope for improvement of the Mg status in rice plants through a breeding approach. Whether utilizing potential allelic variation of MGT1, as well as increased foliar Mg concentration effectively reduces Fe toxicity in the field need to be investigated further.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: NCBI SRA accession number: PRJNA936902 under the following link: https://www.ncbi.nlm.nih.gov/sra/PRJNA936902.
Author contributions
TRaj, MW and GK conceived and designed the field experiment. TRaj, MW, GK and TRak conducted field experiment. TRaj and TRak obtained phenotypic data from the field experiment. HM conducted element analyses of field samples. MW, GJDK, TRak, AA and TRaz supervized the field experiment. TRaj, MW and TRak analyzed data from the field experiment. YU conceived and conducted the hydroponic experiments and RNA-seq and analyzed these data. TRaj and YU wrote the manuscript draft and prepared figures. All authors edited and approved the final manuscript.
Funding
This study was partly supported by Biotechnology and Biological Sciences Research Council (Grant number BB/R020388/1) and JIRCAS research project “Africa Rice Farming System”. MW has been partly funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy – EXC 2070 – 390732324.
Acknowledgments
The authors thank local farmers at Madagascar for field management. The authors also thank M. Yonemoto and M. Matsuyama for the assistance in element analysis.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2023.1213456/full#supplementary-material
References
Arouna, A., Fatognon, I. A., Saito, K., Futakuchi, K. (2021). Moving toward rice self-sufficiency in sub-Saharan Africa by 2030: lessons learned from 10 years of the coalition for African rice development. World Dev. Perspect. 21, 100291. doi: 10.1016/j.wdp.2021.100291
Ashrafuzzaman, M., Haque, Z., Ali, B., Mathew, B., Yu, P., Hochholdinger, F., et al. (2018). Ethylenediurea (EDU) mitigates the negative effects of ozone in rice: insights into its mode of action. Plant Cell Environ. 41, 2882–2898. doi: 10.1111/pce.13423
Ashrafuzzaman, M., Ueda, Y., Frei, M. (2020). Natural sequence variation at the OsORAP1 locus is a marker for ozone tolerance in Asian rice. Environ. Exp. Bot. 178, 104153. doi: 10.1016/j.envexpbot.2020.104153
Audebert, A., Fofana, M. (2009). Rice yield gap due to iron toxicity in West Africa. J. Agron. Crop Sci. 195, 66–76. doi: 10.1111/j.1439-037X.2008.00339.x
Aung, M. S., Masuda, H. (2020). How does rice defend against excess iron?: physiological and molecular mechanisms. Front. Plant Sci. 11. doi: 10.3389/fpls.2020.01102
Aung, M. S., Masuda, H., Kobayashi, T., Nishizawa, N. K. (2018). Physiological and transcriptomic analysis of responses to different levels of iron excess stress in various rice tissues. Soil Sci. Plant Nutr. 64, 370–385. doi: 10.1080/00380768.2018.1443754
Balzergue, C., Dartevelle, T., Godon, C., Laugier, E., Meisrimler, C., Teulon, J. M., et al. (2017). Low phosphate activates STOP1-ALMT1 to rapidly inhibit root cell elongation. Nat. Commun. 8, 15300. doi: 10.1038/ncomms15300
Bates, D., Mächler, M., Bolker, B., Walker, S. (2015). Fitting linear mixed-effects models using lme4. J. Stat. Software 67, 1–48. doi: 10.18637/jss.v067.i01
Becker, M., Asch, F. (2005). Iron toxicity in rice–conditions and management concepts. J. Plant Nutr. Soil Sci. 168, 558–573. doi: 10.1002/jpln.200520504
Benckiser, G., Santiago, S., Neue, H. U., Watanabe, I., Ottow, J. C. G. (1984). Effect of fertilization on exudation, dehydrogenase activity, iron-reducing populations and fe++ formation in the rhizosphere of rice (Oryza sativa l.) in relation to iron toxicity. Plant Soil 79, 305–316. doi: 10.1007/BF02184324
Bolger, A. M., Lohse, M., Usadel, B. (2014). Trimmomatic: a flexible trimmer for illumina sequencing data. Bioinformatics 30, 2114–2120. doi: 10.1093/bioinformatics/btu170
Chen, H. C., Cheng, W. H., Hong, C. Y., Chang, Y.S., Chang, M. C. (2018). The transcription factor OsbHLH035 mediates seed germination and enables seedling recovery from salt stress through ABA-dependent and ABA-independent pathways, respectively. Rice 11, 50. doi: 10.1186/s12284-018-0244-z
Chen, Z. C., Ma, J. F. (2013). Magnesium transporters and their role in Al tolerance in plants. Plant Soil 368, 51–56. doi: 10.1007/s11104-012-1433-y
Chen, Z. C., Yamaji, N., Horie, T., Che, J., Li, J., An, G., et al. (2017). A magnesium transporter OsMGT1 plays a critical role in salt tolerance in rice. Plant Physiol. 174, 1837–1849. doi: 10.1104/pp.17.00532
Chen, Z. C., Yamaji, N., Motoyama, R., Nagamura, Y., Ma, J. F. (2012). Up-regulation of a magnesium transporter gene osmgt1 is required for conferring aluminum tolerance in rice. Plant Physiol. 159, 1624–1633. doi: 10.1104/pp.112.199778
Diop, B., Wang, D. R., Drame, K. N., Gracen, V., Tongoona, P., Dzidzienyo, D., et al. (2020). Bridging old and new: diversity and evaluation of high iron-associated stress response of rice cultivated in West Africa. J. Exp. Bot. 71, 4188–4200. doi: 10.1093/jxb/eraa182
Dobermann, A., Fairhurst, T. H. (2000). Rice: nutrient disorders and nutrient management. Int. Rice Res. Instit. Manila Philippines. 1, 1–203.
Duan, Y. B., Li, J., Qin, R. Y., Xu, R. F., Li, H., Yang, Y. C., et al. (2016). Identification of a regulatory element responsible for salt induction of rice OsRAV2 through ex situ and in situ promoter analysis. Plant Mol. Biol. 90, 49–62. doi: 10.1007/s11103-015-0393-z
Dufey, I., Draye, X., Lutts, S., Lorieux, M., Martinez, C., Bertin, P. (2015a). Novel QTLs in an interspecific backcross Oryza sativa × Oryza glaberrima for resistance to iron toxicity in rice. Euphytica 204, 609–625. doi: 10.1007/s10681-014-1342-7
Dufey, I., Hakizimana, P., Draye, X., Lutts, S., Bertin, P. (2009). QTL mapping for biomass and physiological parameters linked to resistance mechanisms to ferrous iron toxicity in rice. Euphytica 167, 143–160. doi: 10.1007/s10681-008-9870-7
Dufey, I., Mathieu, A. S., Draye, X., Lutts, S., Bertin, P. (2015b). Construction of an integrated map through comparative studies allows the identification of candidate regions for resistance to ferrous iron toxicity in rice. Euphytica 203, 59–69. doi: 10.1007/s10681-014-1255-5
Engel, K., Asch, F., Becker, M. (2012). Classification of rice genotypes based on their mechanisms of adaptation to iron toxicity. J. Plant Nutr. Soil Sci. 175, 871–881. doi: 10.1002/jpln.201100421
Fageria, N. K., Santos, A. B., Barbosa Filho, M. P., Guimarães, C. M. (2008). Iron toxicity in lowland rice. J. Plant Nutr. 31, 1676–1697. doi: 10.1080/01904160802244902
Fang, Y., Xie, K., Xiong, L. (2014). Conserved miR164-targeted NAC genes negatively regulate drought resistance in rice. J. Exp. Bot. 65, 2119–2135. doi: 10.1093/jxb/eru072
FAO (2017) Rice market monitor 21. Available at: www.fao.org/economic/est/publications/rice-publications/rice-market-monitor-rmm/.
Frei, M., Tanaka, J. P., Chen, C. P., Wissuwa, M. (2010). Mechanisms of ozone tolerance in rice: characterization of two QTLs affecting leaf bronzing by gene expression profiling and biochemical analyses. J. Exp. Bot. 61, 1405–1417. doi: 10.1093/jxb/erq007
Gao, P. P., Zheng, G. H., Wu, Y. H., Liu, P. (2014). Effect of exogenous potassium on photosynthesis and antioxidant enzymes of rice under iron toxicity. Russ. J. Plant Physiol. 61, 47–52. doi: 10.1134/S1021443714010051
Genon, J. G., de Hepcee, N., Duffy, J. E., Delvaux, B., Hennebert, P. A. (1994). Iron toxicity and other chemical soil constraints to rice in highland swamps of Burundi. Plant Soil 166, 109–115. doi: 10.1007/BF02185487
Hartmann, J., Asch, F. (2018). Micro-method to determine iron concentrations in plant tissues using 2,2′ bipyridine. J. Plant Nutr. Soil Sci. 181, 357–363. doi: 10.1002/jpln.201700433
Hauer-Jákli, M., Tränkner, M. (2019). Critical leaf magnesium thresholds and the impact of magnesium on plant growth and photo-oxidative defense: a systematic review and meta-analysis from 70 years of research. Front. Plant Sci. 10. doi: 10.3389/fpls.2019.00766
Hodges, D. M., DeLong, J. M., Forney, C. F., Prange, R. K. (1999). Improving the thiobarbituric acid-reactive-substances assay for estimating lipid peroxidation in plant tissues containing anthocyanin and other interfering compounds. Planta 207, 604–611. doi: 10.1007/s004250050524
Höller, S., Ueda, Y., Wu, L., Wang, Y., Hajirezaei, M. R., Ghaffari, M. R., et al. (2015). Ascorbate biosynthesis and its involvement in stress tolerance and plant development in rice (Oryza sativa l.). Plant Mol. Biol. 88, 545–560. doi: 10.1007/s11103-015-0341-y
Howeler, R. H. (1973). Iron-induced oranging disease of rice in relation to physico-chemical changes in a flooded oxisol. Soil Sci. Soc Am. J. 37, 898–903. doi: 10.2136/sssaj1973.03615995003700060030x
Hu, B., Jiang, Z., Wang, W., Qiu, Y., Zhang, Z., Liu, Y., et al. (2019). Nitrate–NRT1.1B–SPX4 cascade integrates nitrogen and phosphorus signalling networks in plants. Nat. Plants 5, 401–413. doi: 10.1038/s41477-019-0384-1
Ishimaru, Y., Bashir, K., Fujimoto, M., An, G., Itai, R. N., Tsutsumi, N., et al. (2009). Rice-specific mitochondrial iron-regulated gene (MIR) plays an important role in iron homeostasis. Mol. Plant 2, 1059–1066. doi: 10.1093/mp/ssp051
Kakei, Y., Masuda, H., Nishizawa, N. K., Hattori, H., Aung, M. S. (2021). Elucidation of novel cis-regulatory elements and promoter structures involved in iron excess response mechanisms in rice using a bioinformatics approach. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.660303
Kim, D., Langmead, B., Salzberg, S. L. (2015). HISAT: a fast spliced aligner with low memory requirements. Nat. Methods 12, 357–360. doi: 10.1038/nmeth.3317
Kinsella, R. J., Kähäri, A., Haider, S., Zamora, J., Proctor, G., Spudich, G., et al. (2011). Ensembl BioMarts: a hub for data retrieval across taxonomic space. Database 2011, bar030. doi: 10.1093/database/bar030
Kirk, G. J. D., Manwaring, H. R., Ueda, Y., Semwal, V. K., Wissuwa, M. (2022). Below-ground plant–soil interactions affecting adaptations of rice to iron toxicity. Plant Cell Environ. 45, 705–718. doi: 10.1111/pce.14199
Kobayashi, N. I., Ogura, T., Takagi, K., Sugita, R., Suzuki, H., Iwata, R., et al. (2018). Magnesium deficiency damages the youngest mature leaf in rice through tissue-specific iron toxicity. Plant Soil 428, 137–152. doi: 10.1007/s11104-018-3658-x
Li, H., Yang, X., Luo, A. (2001). Ameliorating effect of potassium on iron toxicity in hybrid rice. J. Plant Nutr. 24, 1849–1860. doi: 10.1081/PLN-100107598
Liu, S., Fukumoto, T., Gena, P., Feng, P., Sun, Q., Li, Q., et al. (2020). Ectopic expression of a rice plasma membrane intrinsic protein (OsPIP1;3) promotes plant growth and water uptake. Plant J. 102, 779–796. doi: 10.1111/tpj.14662
Love, M. I., Huber, W., Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550. doi: 10.1186/s13059-014-0550-8
Matthus, E., Wu, L. B., Ueda, Y., Höller, S., Becker, M., Frei, M. (2015). Loci, genes, and mechanisms associated with tolerance to ferrous iron toxicity in rice (Oryza sativa l.). Theor. Appl. Genet. 128, 2085–2098. doi: 10.1007/s00122-015-2569-y
McLeay, R. C., Bailey, T. L. (2010). Motif enrichment analysis: a unified framework and an evaluation on ChIP data. BMC Bioinf. 11, 165. doi: 10.1186/1471-2105-11-165
Melandri, G., Sikirou, M., Arbelaez, J. D., Shittu, A., Semwal, V. K., Konaté, K. A., et al. (2021). Multiple small-effect alleles of Indica origin enhance high iron-associated stress tolerance in rice under field conditions in West Africa. Front. Plant Sci. 11. doi: 10.3389/fpls.2020.604938
Moore, P. A., Attanandana, T., Patrick, W. H. (1990). Factors affecting rice growth on acid sulfate soils. Soil Sci. Soc Am. J. 54, 1651–1656. doi: 10.2136/sssaj1990.03615995005400060024x
Norton, G. J., Deacon, C. M., Xiong, L., Huang, S., Meharg, A. A., Price, A. H. (2010). Genetic mapping of the rice ionome in leaves and grain: identification of QTLs for 17 elements including arsenic, cadmium, iron and selenium. Plant Soil 329, 139–153. doi: 10.1007/s11104-009-0141-8
Ogo, Y., Nakanishi Itai, R., Nakanishi, H., Kobayashi, T., Takahashi, M., Mori, S., et al. (2007). The rice bHLH protein OsIRO2 is an essential regulator of the genes involved in fe uptake under fe-deficient conditions. Plant J. 51, 366–377. doi: 10.1111/j.1365-313X.2007.03149.x
O’Malley, R. C., Huang, S. S. C., Song, L., Lewsey, M. G., Bartlett, A., Nery, J. R., et al. (2016). Cistrome and epicistrome features shape the regulatory DNA landscape. Cell 165, 1280–1292. doi: 10.1016/j.cell.2016.04.038
Pandey, P., Srivastava, R. K., Dubey, R. S. (2013). Salicylic acid alleviates aluminum toxicity in rice seedlings better than magnesium and calcium by reducing aluminum uptake, suppressing oxidative damage and increasing antioxidative defense. Ecotoxicology 22, 656–670. doi: 10.1007/s10646-013-1058-9
Pertea, M., Pertea, G. M., Antonescu, C. M., Chang, T. C., Mendell, J. T., Salzberg, S. L. (2015). StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 33, 290–295. doi: 10.1038/nbt.3122
Pertea, M., Kim, D., Pertea, G. M., Leek, J. T., Salzberg, S. L. (2016) Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat. Protoc. 11, 1650–1667. doi: 10.1038/nprot.2016-095
Rajonandraina, T., Rakotoson, T., Wissuwa, M., Ueda, Y., Razafimbelo, T., Andriamananjara, A., et al. (2023). Mechanisms of genotypic differences in tolerance of iron toxicity in field-grown rice. Field Crops Res. 298, 108953. doi: 10.1016/j.fcr.2023.108953
Rakotoson, T., Ergezinger, L., Rajonandraina, T., Razafimbelo, T., Wu, L. B., Frei, M. (2019). Physiological investigations of management and genotype options for adapting rice production to iron toxicity in Madagascar. J. Plant Nutr. Soil Sci. 182, 485–495. doi: 10.1002/jpln.201800621
Rengel, Z., Bose, J., Chen, Q., Tripathi, B. N. (2015). Magnesium alleviates plant toxicity of aluminium and heavy metals. Crop Pasture Sci. 66, 1298–1307. doi: 10.1071/CP15284
Sahrawat, K. L. (2004). Iron toxicity in wetland rice and the role of other nutrients. J. Plant Nutr. 27, 1471–1504. doi: 10.1081/PLN-200025869
Sakai, H., Lee, S. S., Tanaka, T., Numa, H., Kim, J., Kawahara, Y., et al. (2013). Rice annotation project database (RAP-DB): an integrative and interactive database for rice genomics. Plant Cell Physiol. 54, e6. doi: 10.1093/pcp/pcs183
Silva, I. R., Smyth, T. J., Israel, D. W., Raper, C. D., Rufty, T. W. (2001a). Magnesium ameliorates aluminum rhizotoxicity in soybean by increasing citric acid production and exudation by roots. Plant Cell Physiol. 42, 546–554. doi: 10.1093/pcp/pce067
Silva, I. R., Smyth, T. J., Israel, D. W., Rufty, T. W. (2001b). Altered aluminum inhibition of soybean root elongation in the presence of magnesium. Plant Soil 230, 223–230. doi: 10.1023/A:1010384516517
Stein, R. J., Ricachenevsky, F. K., Fett, J. P. (2009). Differential regulation of the two rice ferritin genes (OsFER1 and OsFER2). Plant Sci. 177, 563–569. doi: 10.1016/j.plantsci.2009.08.001
Sun, X., Kay, A. D., Kang, H., Small, G. E., Liu, G., Zhou, X., et al. (2013). Correlated biogeographic variation of magnesium across trophic levels in a terrestrial food chain. PloS One 8, e78444. doi: 10.1371/journal.pone.0078444
Suriyagoda, L. D. B., Sirisena, D. N., Somaweera, K. A. T. N., Dissanayake, A., De Costa, W. A. J. M., Lambers, H. (2017). Incorporation of dolomite reduces iron toxicity, enhances growth and yield, and improves phosphorus and potassium nutrition in lowland rice (Oryza sativa l). Plant Soil 410, 299–312. doi: 10.1007/s11104-016-3012-0
Sylla, M. (1994). Soil acidity and salinity: spatial variability and effects on rice production in West africa’s mangrove zone (Wageningen: Wageningen University).
Tadano, T. (1975). Devices of rice roots to tolerate high iron concentration in growth media. Japan Agric. Res. Q. 9, 34–39.
Tanoi, K., Kobayashi, N. I. (2015). Leaf senescence by magnesium deficiency. Plants 4, 756–772. doi: 10.3390/plants4040756
Todaka, D., Nakashima, K., Shinozaki, K., Yamaguchi-Shinozaki, K. (2012). Toward understanding transcriptional regulatory networks in abiotic stress responses and tolerance in rice. Rice 5, 6. doi: 10.1186/1939-8433-5-6
Tränkner, M., Jákli, B., Tavakol, E., Geilfus, C. M., Cakmak, I., Dittert, K., et al. (2016). Magnesium deficiency decreases biomass water-use efficiency and increases leaf water-use efficiency and oxidative stress in barley plants. Plant Soil 406, 409–423. doi: 10.1007/s11104-016-2886-1
Tsujimoto, Y., Rakotoson, T., Tanaka, A., Saito, K. (2019). Challenges and opportunities for improving n use efficiency for rice production in sub-Saharan Africa. Plant Prod. Sci. 22, 413–427. doi: 10.1080/1343943X.2019.1617638
Ueda, Y., Kiba, T., Yanagisawa, S. (2020). Nitrate-inducible NIGT1 proteins modulate phosphate uptake and starvation signalling via transcriptional regulation of SPX genes. Plant J. 102, 448–466. doi: 10.1111/tpj.14637
Ueda, Y., Wissuwa, M. (2022). Physiological evidence that nitrate use positively correlates with internal phosphorus utilization efficiency and phosphorus uptake efficiency in rice (Oryza sativa l.). Plant Soil 481, 547–561. doi: 10.1007/s11104-022-05655-3
van Oort, P. A. J. (2018). Mapping abiotic stresses for rice in Africa: drought, cold, iron toxicity, salinity and sodicity. Field Crops Res. 219, 55–75. doi: 10.1016/j.fcr.2018.01.016
Wairich, A., De Oliveira, B. H. N., Wu, L. B., Murugaiyan, V., Margis-Pinheiro, M., Fett, J. P., et al. (2021). Chromosomal introgressions from Oryza meridionalis into domesticated rice Oryza sativa result in iron tolerance. J. Exp. Bot. 72, 2242–2259. doi: 10.1093/jxb/eraa461
Wang, Y., Frei, M., Wissuwa, M. (2008). An agar nutrient solution technique as a screening tool for tolerance to zinc deficiency and iron toxicity in rice. Soil Sci. Plant Nutr. 54, 744–750. doi: 10.1111/j.1747-0765.2008.00302.x
Wopereis, M. C. S., Diagne, A., Johnson, D. E., Seck, P. A. (2013). Realizing africa’s rice promise: priorities for action. Realizing Africa’s Rice promise 424–436. doi: 10.1079/9781845938123.0424
Wu, L. B., Holtkamp, F., Wairich, A., Frei, M. (2019). Potassium ion channel gene OsAKT1 affects iron translocation in rice plants exposed to iron toxicity. Front. Plant Sci. 10. doi: 10.3389/fpls.2019.00579
Wu, L., Shhadi, M., Gregorio, G., Matthus, E., Becker, M., Frei, M. (2014). Genetic and physiological analysis of tolerance to acute iron toxicity in rice. Rice 7, 8. doi: 10.1186/s12284-014-0008-3
Wu, L.-B., Ueda, Y., Lai, S.-K., Frei, M. (2017). Shoot tolerance mechanisms to iron toxicity in rice (Oryza sativa l.). Plant Cell Environ. 40, 570–584. doi: 10.1111/pce.12733
Yamaji, N., Huang, C. F., Nagao, S., Yano, M., Sato, Y., Nagamura, Y., et al. (2009). A zinc finger transcription factor ART1 regulates multiple genes implicated in aluminum tolerance in rice. Plant Cell 21, 3339–3349. doi: 10.1105/tpc.109.070771
Yang, M., Lu, K., Zhao, F. J., Xie, W., Ramakrishna, P., Wang, G., et al. (2018). Genome-wide association studies reveal the genetic basis of ionomic variation in rice. Plant Cell 30, 2720–2740. doi: 10.1105/tpc.18.00375
Yin, X., Cui, Y., Wang, M., Xia, X. (2017). Overexpression of a novel MYB-related transcription factor, OsMYBR1, confers improved drought tolerance and decreased ABA sensitivity in rice. Biochem. Biophys. Res. Commun. 490, 1355–1361. doi: 10.1016/j.bbrc.2017.07.029
Yoshida, S., Forno, D. A., Cock, J. H., Gomez, K. A. (1976). Laboratory manual for physiological studies of rice (Manila: The International Rice Research Institute).
Keywords: rice, iron toxicity, magnesium, leaf bronzing, fertilizer, RNA-Seq
Citation: Rajonandraina T, Ueda Y, Wissuwa M, Kirk GJD, Rakotoson T, Manwaring H, Andriamananjara A and Razafimbelo T (2023) Magnesium supply alleviates iron toxicity-induced leaf bronzing in rice through exclusion and tissue-tolerance mechanisms. Front. Plant Sci. 14:1213456. doi: 10.3389/fpls.2023.1213456
Received: 28 April 2023; Accepted: 04 July 2023;
Published: 21 July 2023.
Edited by:
Felipe Klein Ricachenevsky, Federal University of Rio Grande do Sul, BrazilCopyright © 2023 Rajonandraina, Ueda, Wissuwa, Kirk, Rakotoson, Manwaring, Andriamananjara and Razafimbelo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yoshiaki Ueda, eW9zaGlha2l1ZWRhQGFmZnJjLmdvLmpw
†These authors have contributed equally to this work