Junji Inoue
Takamasa Teramoto
Tomohiko Kazama
Takahiro Nakamura- Faculty of Agriculture, Kyushu University, Fukuoka, Japan
Cd is a seriously hazardous heavy metal for both plants and humans and international regulations regarding Cd intake have become stricter in recent years. Three-quarters of the Cd intake comes from plant-based foods, half of which comes from cereals. Therefore, it is anticipated that the Cd uptake efficiency of cereals, including rice, a staple crop in Asia, will be reduced. Natural resistance-associated macrophage protein (Nramp) is the principal transporter involved in the uptake and translocation of metal ions in various plants. In rice, OsNramp5 is a transporter of Mn, which is an essential micronutrient for plant growth, and is responsible for Cd uptake. Although several attempts have been made to engineer the metal uptake characteristics of OsNramp5, in many cases, both Cd and Mn uptake efficiencies are impaired. Therefore, in this study, we engineered OsNramp5 to reduce Cd uptake while retaining Mn uptake efficiency for low-Cd rice production. OsNramp5 was engineered using amino acid substitution(s) at the 232nd Ala and 235th Met of OsNramp5, which have been suggested to be key residues for metal uptake efficiency and/or selectivity by structural analyses of bacterial Nramps. The metal uptake efficiency was first analyzed using a yeast model assay system. Several mutants showed less than 8.6% Cd and more than 64.1% Mn uptake efficiency compared to the original OsNramp5. The improved metal uptake characteristics were confirmed by direct measurement of the metal content in the yeast using inductively coupled plasma optical emission spectroscopy. Notably, several mutants reduced Cd uptake efficiency to the background level while retaining more than 64.7% Mn uptake efficiency under conditions mimicking heavily polluted soils in the world. In addition, computational structural modeling suggested requirements for the spatial and chemical properties of the metal transport tunnel and metal-binding site, respectively, for Cd/Mn uptake efficiency.
Introduction
All organisms require metal ions as micronutrients to maintain their biological activities. Plants absorb biologically useful metals, including Mn2+, Fe2+, Co2+, Ni2+, Cu2+, and Zn2+, from the soil but accidentally incorporate the toxic heavy metals Cd2+, Pb2+, and Hg2+, which pose significant risks to agricultural and human health. Cd is a hazardous heavy metal that can be easily absorbed and accumulated in plant tissues (Gill and Tuteja, 2011). Cd inhibits plant growth by reducing photosynthesis and oxidative stress (Nazir et al., 2022). Moreover, Cd is incorporated into the human body through the food chain and permanently accumulates primarily in the kidneys, with a long biological half-life of 10–35 years (World Health Organization, 2022). High Cd intake can lead to serious health issues, such as kidney dysfunction, cancer, bone fractures, and itai-itai disease, which causes severe pain, softening of bones, and kidney failure (Nogawa and Kido, 1993; Staessen et al., 1999; Tim et al., 2006; Vervaet et al., 2017). Low to moderate Cd exposure is also associated with kidney dysfunction, cancer, and decreased bone density (Järup and Alfvén, 2004; Julin et al., 2012; García-Esquinas et al., 2014; Wallin et al., 2021).
Cd intake has been strictly regulated owing to its hazardous effects. International standards established a permissible maximum intake of Cd at 7 µg/kg/week (1.00 µg/kg/day) in 1988, which was modified to 25 µg/kg/month (0.83 µg/kg/day) in 2010 (World Health Organization, 2011). Europe implemented a stricter standard for intake of 2.5 µg/kg/week (0.36 µg/kg/day) in 2009 (EFSA, 2009). Cd intake is largely caused by seven commodity groups: rice, wheat, roots, tubers, leafy greens, other vegetables, and mollusks (World Health Organization, 2006), six of which are plant-derived. Indeed, it has been reported that Cd intake mostly comes from plant-based foods (74% and 82% in China and Sweden, respectively) and almost half of Cd intake is from cereals: 46% and 48% in China and Sweden, respectively (Julin et al., 2012; Qing et al., 2023). Thus, there is a major challenge in reducing Cd uptake by rice, which is a staple crop for Asian people.
Numerous studies have been conducted to elucidate the molecular mechanisms of Cd uptake in various plants and have found a significant role in natural resistance-associated macrophage proteins (Nramps). Nramp is a principal transporter family for the uptake and/or translocation of metal ions (Socha and Guerinot, 2014; Clemens and Ma, 2016). Genomic analysis has shown that plants possess a significantly larger number of Nramp genes than mammals and bacteria (Chen et al., 2021). In rice, seven Nramp genes have been found and Oryza Sativa Nramp5 (OsNramp5) was identified to be largely responsible for the uptake of Cd and Mn via roots from the surrounding environment (Ishimaru et al., 2012), as evidenced by OsNramp5 knockout lines exhibiting impaired Mn/Cd uptake efficiency (Ishikawa et al., 2012; Sasaki et al., 2012).
Various attempts have been made to reduce Cd uptake efficiency by focusing on Nramps. Simple knockdown of Nramp genes resulted in a reduction in Mn and Cd uptake and had negative effects on growth in Mn-deficient environments in rice (Yang et al., 2014; Honma et al., 2017) and Arabidopsis (Cailliatte et al., 2010; Gao et al., 2018). This is because Mn is an essential micronutrient for plant growth and development, involving several metabolic pathways, including photosynthesis, and as a cofactor of enzymes. A single amino acid substitution mutant of OsNramp5 Q337K mutant identified through targeting induced local lesions in genomes screening showed a 50% reduction in both Mn and Cd uptake. Notably, a field test of rice containing the Q337K mutation exhibited lower Cd uptake than the wild type and higher growth tolerance to Mn deficiency than knockout lines (Kuramata et al., 2022). Another amino acid substitution mutant, A512T OsNramp5 was generated by random mutagenesis. The A512T mutant showed a lower Cd uptake efficiency while maintaining Mn uptake efficiency in a yeast model assay system, although the evaluation in rice plants remains (Qu and Nakanishi, 2023). Studies on Nramps involved in Mn/Cd uptake in other plants have suggested that metal uptake selectivity varies depending on the protein species, for example, differences in a few amino acids (Curie et al., 2000; Thomine et al., 2000; Cailliatte et al., 2010).
Structural analyses of bacterial Nramps provide valuable information regarding the molecular basis of metal binding and selectivity (Ehrnstorfer et al., 2014; Bozzi et al., 2016b; Ray et al., 2023). The Nramp family consists of 11 or 12 transmembrane regions (TMs) and the broken helical regions in TM1 and TM6 comprise the substrate-binding sites. Structural studies also revealed critical residues for metal binding (#2, #5, #8, and #11 in Figure 1) in TM1 and TM6. Several attempts have been made based on structural information and several mutant Nramps with altered metal uptake selectivity (Bozzi et al., 2016a; Li et al., 2018; Lu et al., 2018). However, structure-based engineering of OsNramp5 has not yet been performed.
Figure 1. Structure of OsNramp5 and putative amino acids of the metal-binding site involving metal uptake efficiency and selectivity. (A) Topological diagram of OsNramp5 structure. The diagram was depicted based on the domain architecture of UniProt Q8H4H5, with modifications of Figure 3A in Bozzi and Gaudet, 2021. The OsNramp5 consists of 12 transmembrane regions. The asterisk denotes the location of the metal-binding site composed of TM1 and TM6. The putative key amino acid residues for the metal binding are designated as #2, #5, #8, and #11. (B) Phylogenetic tree and amino acid composition of metal-binding sites of 21 Nramps in various organisms (Os, Oryza sativa; Zm, Zea mays; At, Arabidopsis thaliana; Sca, Staphylococcus capitis; Dra, Deinococcus radiodurans). The amino acid species of #1 to #12 are shown in capital single letter, with the metal uptake characteristics: +, transport; -, no transport;?, promiscuous;/, not confirmed (a: Curie et al., 2000, b: Takahashi et al., 2011, c: Tiwari et al., 2014, d: Chang et al., 2020, e: Ishimaru et al., 2012, f: Ishikawa et al., 2012, g: Sasaki et al., 2012, h: Sui et al., 2018, i: Peris-Peris et al., 2017, j: Xia et al., 2010, k: Li et al., 2022, l: Yamaji et al., 2013, m: Yang et al., 2013, n: Thomine et al., 2000, o: Cailliatte et al., 2009, p: Li et al., 2019, q: Pottier et al., 2015, r: Alejandro et al., 2017, s: Gao et al., 2018, t: Zhao et al., 2018, u: Guo et al., 2022, v: Li et al., 2023, w: Ehrnstorfer et al., 2014, x: Bozzi et al., 2016a, y: Bozzi et al., 2016b).
In this study, we engineered the Mn/Cd uptake efficiencies of OsNramp5 by incorporating amino acid substitution(s) at the putative metal-binding residue(s) to produce rice plants that reduce Cd uptake efficiencies but retain that of Mn. The Mn/Cd uptake efficiency of the mutant OsNramp5 was first evaluated using a yeast model assay system (Podar et al., 2012; Chen et al., 2016). Together with further evaluation by direct measurement of the Mn/Cd uptake efficiencies, we identified several OsNramp5 mutants that exhibited desirable Mn/Cd uptake characteristics and retained more than 64.7% Mn uptake but reduced Cd uptake efficiency to background levels under the conditions mimicking a polluted soil. The molecular mechanisms of Mn/Cd uptake by OsNramp5 and mutant proteins are also discussed based on computational structural modeling.
Materials and methods
Sequence alignment and phylogenetic analysis
The amino acid sequences of the 21 Nramps used in this study were obtained from the NCBI database (https://www.ncbi.nlm.nih.gov) and are listed in Supplementary Table 1, along with their NCBI accession numbers. Protein sequence alignment and phylogenetic analyses were conducted using MUSCLE (Edgar, 2004) and the Maximum likelihood method (Jones et al., 1992), respectively, in the MEGA 11 package (Tamura et al., 2021).
Plasmid construction of OsNramp5 and the mutants
The DNA sequence of the OsNramp5 coding region was chemically synthesized by Geneart (Thermo Fisher Scientific, Waltham, MA, USA). The synthesized DNA was used as a template for polymerase chain reaction amplification using PrimeSTAR Max (Takara Bio, Shiga, Japan) and primers OsNramp5_GA_F and OsNramp5_GA_R. The polymerase chain reaction product was integrated into the expression plasmid of pDR195 (Rentsch et al., 1995; Addgene #36028), using the NEBuilder HiFi DNA Assembly Cloning Kit (New England Biolabs, Ipswich, MA, USA). The resulting plasmid was used for site-directed mutagenesis of OsNramp5 using PrimeSTAR Max (Takara Bio, Shiga, Japan) and the primers listed in Supplementary Table 2 by inverse PCR. The DNA sequence was confirmed using Sanger sequencing.
Metal-sensitivity assay in yeast
The yeast strains used in this study, BY4741(wild type; WT), Cd-sensitive Δycf1 (Y04069), and Mn-sensitive Δpmr1 (Y04534), were obtained from EUROSCARF (Oberursel, Germany). Plasmids containing the OsNramp5 gene or mutants were introduced into yeast using YNB-ura (MP Biomedicals, Irvine, CA, USA) agar plates and the LiAc/SS carrier DNA/PEG method (Gietz and Schiestl, 2007). The obtained transformants were pre-cultured in liquid YNB-ura medium at 30 °C, overnight with gentle shaking. The pre-cultured yeast was replaced with a new YNB-ura liquid medium by adjusting pH 4.5 and an optical density at 600 nm (OD600) = 0.001 in a final volume of 200 µl. CdSO4 was added to be 0 µM to 81.8 µM for the Cd sensitivity test. YNB-ura medium contains 0.04 mg/l (0.25 µM) CuSO4, 0.2 mg/l (1.23 µM) FeCl3, 0.4 mg/l (2.65 µM) MnSO4, and 0.4 mg/l (2.48 µM) ZnSO4 as basic components.
The Cd sensitivity test was conducted using Cd-sensitive Δycf1 by incubating the yeast in 96 plates with rubber lids at 30°C for 3 days with gentle shaking. After incubation, OD600 was measured using EnSight Multimode Plate Reader (Perkin Elmer, Shelton, CT, USA) and IC50 values were calculated using R package “drc ver. 3.0-1” (Ritz et al., 2015). The relative Cd uptake efficiency was estimated using three biological replicates (independent yeast transformants), with the reciprocal IC50 (1/IC50) of the empty plasmid and WT as 0% and 100%, respectively. The Mn sensitivity test was conducted using Mn-sensitive Δpmr1 and a medium containing 0 mM to 1.5 mM MnSO4.
Quantification of the metal uptake characteristics using inductively coupled plasma optical emission spectroscopy
The WT yeast strain (BY4741) containing empty, OsNramp5, or the mutant plasmid was incubated in liquid YNB-ura medium at 30 °C overnight with gentle shaking. The pre-culture medium was transferred to 5 ml YNB-ura liquid medium containing 20 µM MnSO4, CdSO4, FeSO4, CuSO4, ZnSO4, and CoSO4, and adjusted to pH 4.5 and OD600 = 0.1. After 16 hours of incubation at 30 °C with gentle shaking, the cells were harvested via centrifugation, washed once with 20 mM EDTA, and washed three times with deionized water. After drying up the cell at 60 °C for three days, the cell was lysed by 0.5 ml 60% HNO3 at 60 °C for 2 hours. The concentrations of Cd and Mn in the cell lysates were measured using an Agilent 5800 inductively coupled plasma optical emission spectroscopy (ICP-OES; Agilent, Santa Clara, CA, USA) according to the manufacturer’s instructions. The actual Cd and Mn content (ng/mg) was calculated by dividing the metal content of the cell lysate by the dry yeast mass. The apparent metal uptake was estimated from the actual metal uptake by subtracting the value for the empty plasmid as the background. The relative uptake value was calculated using the apparent metal uptake of the empty plasmid and WT as 0% and 100%, respectively. The Cd and Mn uptake efficiencies were analyzed at various Cd concentrations using a medium containing 1.8 mM MnSO4 and 0, 2.5, or 7.5 µM CdSO4. The relative uptake values were calculated using the apparent metal uptake of the empty plasmid and WT as 0% and 100%, respectively, for each CdSO4 concentration. Metal uptake was estimated using three biological replicates (independent yeast transformants).
Computational structural analysis
The 3D structures of OsNramp5 and the mutants were predicted by AlphaFold2 running on the ColabFold v1.5.5 using MMseqs2 server (Mirdita et al., 2019, 2022; Jumper et al., 2021). The highest-ranking model was relaxed by using amber force fields (Eastman et al., 2017). The X-ray crystal structure of Deinococcus radiodurans Nramp (DraNramp)-G223W (Protein Data Bank ID: 8E6N) was chosen as the template for structural modeling (Ray et al., 2023). AlphaFill was used to predict the optimized conformation by incorporating Mn ions into the metal-binding sites of the predicted 3D structures (Hekkelman et al., 2023). Tunnel detection was performed by setting the tunnel bottleneck (minimum radius) to 0.9 Å, with the tunnel search starting point at approximately 3 Å from the Mn ion in CAVER Analyst 2.0 BETA (Jurcik et al., 2018), and the length, radius of the tunnels were calculated. Detection of cavities around the metal-binding sites was performed using pyKVFinder (Eisenberg et al., 1984; da Silva Guerra et al., 2021, Guerra et al., 2023), as well as hydropathy calculations. The cavity detection was performed with the following settings: 0.9 Å probe size, within 3.5 Å to the direction of the metal transport tunnel entrance from Mn, and 3 Å from the center line of the transport tunnel. The predicted 3D structures were visualized using PyMOL v2.5.0 (https://pymol.org/).
Results
Sequence analysis of Nramps
To identify the putative amino acids of plant Nramps involved in metal binding, the amino acid sequence of the OsNramp5 protein was aligned and analyzed phylogenetically, with other 18 Nramps in monocot and dicot plants, as well as bacterial Nramps of Staphylococcus capitis Nramp (ScaNramp) and DraNramp (Figure 1B; Supplementary Figure 1). The key amino acid residues for the metal transporters were designated #1 to #12. Residues #2, #5, #8, and #11 were identified as metal-contacting residues based on structural analyses of bacterial Nramps (Ehrnstorfer et al., 2014). By focusing on these residues, aspartic acid (D) and asparagine (N) at #2 and #5, respectively, in TM1 were found to be highly conserved among all Nramps from bacteria to plants. In contrast, the metal-binding amino acids at positions #8 and #11 in TM6 are diverse. All Nramps transport Cd, which has been characterized so far, containing methionine (M) at position #11 (OsNramp1, 2, and 5; Arabidopsis thaliana Nramp (AtNramp) 1, 3, 4, and 6; Zea mays Nramp (ZmNramp) 2 and 5; ScaNramp; and DraNramp). Nramp, which has been characterized as not transporting Cd, contains different amino acids at #11, such as valine (V) or threonine (T) (OsNramp6, OsNramp4 (also known as OsNrat1), and ZmNrat1). Threonine (T) at position #11 is associated with unique metal uptake characteristics of the trivalent metal ion Al in OsNramp4 and ZmNrat1. No clear relationship was found between the amino acid species at position #8 and the metal uptake characteristics.
These observations suggest that the amino acids in TM1 (#2 and #5) are indispensable for their crucial function as transporters and that amino acids #8 and #11 in TM6 may be involved in metal uptake selectivity. This hypothesis is supported by previous studies on the metal selectivity modification of amino acid substitutions at position #11 in several Nramps (Bozzi et al., 2016a; Li et al., 2018; Lu et al., 2018). Therefore, we modified residues #8 and #11 to engineer OsNramp5 to reduce Cd uptake efficiency.
Analysis of metal uptake characteristics of OsNramp5 mutants using yeast assay system
The yeast assay system is a well-established system for estimating the metal uptake characteristics of transporter proteins of interest. Metal-sensitive mutant strains of yeast show retarded growth in the accumulation of a particular metal ion (Podar et al., 2012; Chen et al., 2016; Li et al., 2018). The extent of growth retardation due to metal accumulation was evaluated as the metal uptake efficiency of OsNramp5 and the mutants in this study. To evaluate the function of residue #11 of OsNramp5, corresponding to 235th Met, the residue was substituted with other 19 amino acid species. The obtained OsNramp5 mutants, designated as M235X, were transformed to the Cd-sensitive yeast strain (Δycf1; Y04069), together with the previously reported mutants of Q337K and A512T (Kuramata et al., 2022; Qu and Nakanishi, 2023). The growth of Cd-sensitive yeast at various Cd concentrations was measured using OD600 to calculate IC50. The IC50 value was used to estimate the metal uptake efficiency (Supplementary Figure 2). All 19 mutant proteins constructed in this study exhibited significantly lower Cd uptake (less than 20%) than that of WT (Figure 2A). Mn uptake efficiency was also analyzed because OsNramp5 is responsible for Mn uptake, which is an essential micronutrient that facilitates the growth and development of plants. The Mn uptake efficiencies were highly variable among the mutants when the assay was conducted using the Mn-sensitive yeast strain (Δpmr1; Y04534). However, the substitution to alanine (A) and cysteine (C) exhibited more than 60% Mn transport compared to WT (M235A and M235C; Figure 2B). The Mn/Cd uptake ratio was calculated by dividing the relative Mn uptake by the relative Cd uptake, to evaluate the metal selectivity of the OsNramp5 mutants (Table 1). The Mn/Cd uptake ratios of M235A and M235C were 11.3 and 7.5, respectively, indicating that these mutants had the expected characteristics of low Cd uptake while retaining Mn uptake. The Mn/Cd uptake ratio of the previously characterized Q337K allele was estimated as 1.0. The Q337K allele transported both Cd and Mn at a medium level compared with WT (36.0% and 37.2%, respectively). This is consistent with a previous report showing that the Q337K allele does not affect the selectivity of Cd or Mn transport (Kuramata et al., 2022). Another previously characterized mutant, A512T, also showed a Mn/Cd uptake ratio of 0.9, due to an approximately 50% reduction in both Mn and Cd uptake. This result was in contrast to that of a previous study in which the mutant showed only a reduction in Cd uptake in yeast (Qu and Nakanishi, 2023). Taken together, M235A and M235C showed the expected metal uptake characteristics, that is, little Cd uptake and efficient Mn uptake, in the yeast assay system; therefore, these mutants were used for subsequent studies.
Table 1. Relative Cd and Mn uptake efficiency and Mn/Cd ratio of OsNramp5 M235X mutants.
Figure 2. Cd and Mn uptake efficiencies of M235X mutants of OsNramp5. (A) Cd uptake efficiency was estimated for 19 M235X mutants in the yeast assay system, together with the previously reported mutants of Q337K and A512T. This test was conducted using Δycf1 in the presence of 0 µM to 81.8 µM CdSO4. The relative Cd uptake was estimated as shown in Supplementary Figure 2. The mean and standard deviation of relative Cd uptake efficiency are shown (N=3, WT=100%). The mutants exhibited <20% Cd and >60% Mn uptake efficiencies are highlighted in red. (B) Mn uptake efficiencies are shown as in (A). This test was conducted using Δpmr1 in the presence of 0 mM to 1.5 mM MnSO4.
Engineering of OsNramp5 by a double mutation
Amino acid substitutions focusing on residue #11 yielded the OsNramp5 mutants, M235A and M235C, with improved Mn/Cd uptake characteristics. The mutants were further subjected to amino acid substitution at residue #8, another putative residue involved in metal uptake selectivity (Figure 1B). Residue #8 for OsNramp5, corresponding to the 232nd Ala, was substituted with cysteine, serine, or valine; these amino acid species are found at residue #8 in other Nramp families. The 232nd Ala residue was also substituted with methionine, which may play a significant role in Cd and Mn uptake at residue #11 (235th Met in OsNramp5). The metal uptake characteristics of the obtained double mutants were analyzed in a yeast assay system, together with mutants containing a single amino acid substitution at position 232. All double mutants showed less than 20% Cd uptake efficiency compared to WT (Figure 3A). When the Mn uptake efficiencies were examined, three double mutants, A232C+M235A, A232S+M235A, and A232S+M235C, showed Mn uptake efficiencies more than 60% (Figure 3B). Two single mutants (A232M and A232V) showed substantially lower Cd uptake efficiencies, although Mn uptake efficiencies were also impaired. Interestingly, the A232M+M235A mutant was simply a swap of the amino acid species between positions 232 and 235 from the WT. However, A232M+M235A displayed extremely low uptake efficiencies for both Cd and Mn, suggesting that both the position and amino acid species determine the efficiency of Cd/Mn uptake.
Figure 3. Cd and Mn uptake efficiencies of OsNramp5 double mutants. (A) Cd uptake efficiency was estimated for eight double mutants in the yeast assay system, together with four single mutants at residue 232. M235A and M235C were also used as references. This test was conducted using Δycf1 in the presence of 0 µM to 81.8 µM CdSO4. The relative uptake value was calculated using the apparent metal uptake of the empty plasmid and WT as 0% and 100%, respectively. The mean and standard deviation of relative Cd uptake efficiency are shown (N=3, WT=100%). The mutants exhibited <20% Cd and >60% Mn uptake efficiencies are highlighted in red. B) Mn uptake efficiencies are shown as in (A). This test was conducted using Δpmr1 in the presence of 0 mM to 1.5 mM MnSO4.
The Mn/Cd uptake ratios of the double mutants were estimated for the single M235X mutants (Table 2). Three double mutants, A232S+M235A, A232S+M235C, and A232C+M235A, showed better Mn/Cd uptake ratios of 17.0, 13.6, and 12.8, respectively, than the single mutants (11.3, M235A; 7.5, M235C). Notably, A232S+M235A showed the most desirable metal uptake characteristics among the mutants in this study, as shown by the Mn/Cd uptake ratio of 17.0, with an equivalent Mn uptake efficiency (96.4%) compared to that of WT, although the Cd uptake efficiency was significantly reduced (5.7%) (Figure 3A; Table 2).
Table 2. Relative Cd and Mn uptake efficiency and Mn/Cd ratio of OsNramp5 single and double mutants.
Direct measurement of the metal uptake of OsNramp5 mutants in yeast
The yeast assay system is a powerful tool that enables high-throughput screening of various mutants. However, the evaluation relies on the growth of yeast, which may be affected by various inter-/intra-cellular conditions. Therefore, the actual metal uptake efficiency of the five selected mutants (M235A, M235C, A232C+M235A, A232S+M235A, and A232S+M235C) were directly determined using inductively coupled plasma optical emission spectrometry (ICP-OES). The mutants were transformed into the WT yeast strain and cultured overnight in the presence of 20 µM each of Cd, Mn, Fe, Cu, Zn, and Co. The incorporated metal content was determined by ICP-OES, and evaluated as metal uptake efficiency in this study.
ICP-OES analysis showed that the Cd uptake of the five mutants was lower than that of the WT, coinciding with the results of the yeast assay system (Table 3). The relative Cd uptake efficiencies were estimated to range from 10.9% to 22.3% against that of WT. When the Mn uptake was examined, the mutants displayed similar or decreased Mn uptake efficiencies compared to the WT, as shown by the relative values from 54.8% to 103.3% (Table 3). The best Mn/Cd ratio of 7.4 was raised by A232C+M235A. Cd uptake was reduced to 13.9% of WT, with no influence on Mn uptake efficiency (103.3%) (Table 3). The Mn/Cd ratio of the single mutant M235A was 3.5, suggesting that the double mutation at residues #8 and #11 successfully improved metal uptake characteristics.
Table 3. Direct measurement of Cd and Mn content for selected OsNramp5 mutants.
When we measured the concentrations of other heavy metals (Co, Cu, Fe, and Zn) under the same experimental conditions, the five mutants exhibited lower Co uptake efficiencies (Supplementary Table 3). This may be due to the role of OsNramp5 in Co uptake and accumulation (Huang et al., 2024). We also found unexpectedly low Cu accumulation in OsNramp5 and all five mutants, although further evaluation is required to elucidate the comprehensive role of OsNramp5 in the uptake and accumulation of various metals.
Taken together, a series of experiments using a yeast assay system and direct measurements by ICP-OES identified several desirable OsNramp5 mutants. Three mutants were selected for further characterization: M235A, the best single mutant in the yeast assay system (Table 1); A232S+M235A, the best double mutant in the yeast assay system (Table 2); and A232C+M235A, the best mutant in the ICP-OES analysis (Table 3).
Cd uptake characterization of the mutants under various Cd concentrations
In natural conditions, the average soil Cd concentration has been reported as 0.36 mg/kg globally, 0.27 mg/kg in the USA, 0.15 mg/kg in Europe, and 0.27 mg/kg in China (Hutchinson et al., 1987; Zhang et al., 2015; Kubier et al., 2019; Ballabio et al., 2023). In Japan, the soil Cd concentration was reported as an average of 0.30 mg/kg and a maximum of 0.80 mg/kg (Asami et al., 1988). Therefore, the Cd uptake of the representative OsNramp5 mutants (M235A, A232S+M235A, and, A232C+M235A) was examined under three conditions: 0 µM for the blank; 2.5 µM (≈ 0.28 mg/kg), mimicking the normal condition; and 7.5 µM (≈ 0.84 mg/kg), mimicking the highly polluted condition. The Mn concentration was determined at 1.8 mM (≈ 98.9 mg/kg) because the previous study of OsNramp5 knockout mutant showed that brown spot disease started to be observed at a soil concentration of 51.9 mg/kg Mn, but not of 220.8 mg/kg Mn (Honma et al., 2017).
This analysis showed that the Cd uptake of all the selected OsNramp5 mutants (A232S+M235A, A232C+M235A, and M235A) was estimated at background levels, even at high Cd concentrations (Figure 4A). A232S+M235A, A232C+M235A, and M235A showed constant Mn uptake under all Cd concentrations: 91.9% to 107.9% of WT in the absence of Cd, 79.8% to 84.0% of WT in the 2.5 µM Cd, and 64.7% to 76.1% of WT at the 7.5 µM Cd (Figure 4B; Table 4). Altogether, this dose-dependent evaluation further strengthened the fact that the three mutants (M235A, A232C+M235A, and A232S+M235A) showed desirable Cd and Mn uptake characteristics under conditions mimicking polluted soil.
Figure 4. Cd and Mn uptake in various Cd concentrations. (A) Apparent Cd content measured using ICP-OES for the WT yeast strain (BY4741) containing OsNramp5 (WT) and the mutants (M235A, A232C+M235A, and A232S+M235A) in the presence of different Cd concentrations (0, 1, 2.5, and 7.5 μM) and fixed 1.8 mM MnSO4. The mean and standard deviation were shown (N = 3). (B) Apparent Mn content measured in the same condition as (A).
Table 4. Apparent Cd and Mn uptake in various Cd concentrations for selected OsNramp5 mutants.
Discussion
Phylogenetic analysis and metal uptake characterization of OsNramp5 mutants using a yeast assay system and direct measurement of metal uptake by ICP-OES successfully identified several OsNramp5 mutants that effectively transported Mn but reduced Cd uptake efficiency compared to the original OsNramp5. The compatibility between the results of the yeast assay system and ICP-OES was evaluated by estimating the correlation coefficient. That was estimated as R=0.95 with p < 0.01 and R=0.85 with p < 0.01 for the Cd and Mn uptake, respectively (Supplementary Figure 3), indicating that both experimental methods were sufficiently worked to evaluate metal uptake characteristics of OsNramp5 mutants. Notably, the Cd dose-dependent test showed that the three mutants (A232S+M235A, A232C+M235A, and M235A) retained more than 64.7% Mn uptake efficiency but background levels of Cd uptake in the conditions mimicking the polluted soil.
Phylogenetic analysis of OsNramp5 and previous studies suggested key residues of TM1 and TM6 for metal uptake efficiency and selectivity. Aspartic acid and asparagine are highly conserved at residues #2 and #5 of TM1, respectively, in various metal ion transporters and are presumably essential for their function as transporters (Figure 1; Chaloupka et al., 2005; Bozzi et al., 2020). Key residues #8 and #11 are degenerated in the Nramp family (Figure 1). The involvement of residues #8 and #11 in metal uptake selectivity was partially demonstrated in this study. All amino acid substitutions at residue #11, M235X of OsNramp5, showed a significant reduction in Cd uptake efficiency compared to the WT (<20%), suggesting that methionine at the metal-binding site is crucial for Cd transport. Mutations #8 and #11 in TM6 frequently reduced the Mn uptake efficiency, whereas efficient Mn transport was achieved when methionine (WT), alanine, and cysteine were positioned at either #8 or #11 (Figures 2B, 3B).
To gain mechanistic insights into the Cd/Mn uptake efficiencies of the OsNramp5 mutants analyzed in this study, their 3D structures were computationally predicted. Several amino acids surrounding the predicted metal transport tunnel of the original OsNramp5 have been reported to affect the metal selectivity of various Nramps (73rd Asn, 225th Asp, 228th Ala, 393th Ile, and 396th Ser; Supplementary Figure 4) (Pottier et al., 2015; Bozzi et al., 2019; Sun et al., 2019), suggesting the reliability of the predicted structure. The predicted structure suggested several requirements that may determine the Cd/Mn uptake efficiencies: the tunnel length, minimum radius as a tunnel width, and hydropathy around the metal-binding sites of the predicted structures (Supplementary Table 4).
For the high Mn uptake efficiency (>60% in Tables 1 and 2), the first requirement was supposed to be tunnel length, as shown that the mutants showed more than 60% Mn uptake efficiency have tunnel lengths between 10.88 and 13.61 Å (Figure 5A and Supplementary Figure 5). In this first node, 17 proteins were eliminated from 34 samples (WT and 33 mutants). A typical example of a change in the tunnel length and shape was observed in the predicted structure of M235P, whose tunnel length was predicted as 10.18 Å (Figure 5B). Proline substitutions are not easily accommodated in transmembrane helices, especially near the middle of the helix (Yohannan et al., 2004). The second requirement was suggested to be a tunnel minimum radius between 1.24 and 1.46 Å. Eight out of 17 proteins were eliminated by this step, as shown in the typical example of the predicted structure of M235V with the tunnel minimum radius of 1.14 Å (Figure 5B). The third requirement is hypothesized to be the hydrophilicity of the metal-binding sites. The mutants that showed more than 60% Mn uptake efficiency have hydrophilicity between -0.30 and -0.03, and six proteins (WT and five mutants) were passed in this step. It is presumed that an appropriate hydrophilicity is required to incorporate metal ions after they pass through the transport tunnel.
Figure 5. Structure and chemical properties that may determine the Cd and Mn uptake efficiency of OsNramp5. (A) Appropriate tunnel length, tunnel minimum radius, and hydropathy around the metal-binding site hypothesized to be required for efficient Mn uptake (>60% of WT). The number of mutants that passed in the node is shown. (B) Structure of the predicted transport tunnel and metal-binding site of WT, M235A, M235P, M235V, and M235S. The predicted tunnels are visualized by blue mesh, together with the key residues. The predicted cavity around the metal-binding site is shown by a small sphere, which is colored according to the Eisenberg & Weiss hydropathy scale. The properties of the length (L) and minimum radius (M) of the transport tunnel, and hydropathy around the metal-binding site (H), are shown.
For the Cd uptake efficiency, methionine at residue #11 was critical for Cd transport. All amino acid substitutions of residue #11 and M235X of OsNramp5 showed a significant reduction of Cd transport efficiency compared to the WT (<20%). The necessity of methionine at residue #11 for Cd uptake has been discussed in a previous report, which provides the significant stabilization necessary for the binding and transport of Cd, which can forge strong covalent-like interactions with sulfur (Bozzi et al., 2016a). The low to moderate Cd uptake efficiency due to the mutation of other positions (A232C, Q337K, and A512T) can be explained by the requirement of tunnel length and minimum radius for efficient Mn uptake. We speculate that these three proteins did not allow Cd ions to reach the metal-binding site because their metal transport tunnels did not match the requirements of either or both spatial properties for metal transport (Supplementary Figure 5). For example, the tunnel length and tunnel minimum radius of the A232C mutant were estimated as 7.44 Å and 1.09 Å, respectively. The sizes did not match the requirements (tunnel lengths between 10.88 and 13.61 Å, and tunnel minimum radius between 1.24 and 1.46 Å; Supplementary Figure 5 and Supplementary Table 4). In brief, structural prediction and mutant analysis in the yeast assay system suggested that metal uptake efficiency may be determined by the tunnel length, tunnel minimum radius, and hydropathy around the metal-binding site. It is interesting to see these observations can be applied to the semi-rational engineering of other Nramps.
This study identified several OsNramp5 mutants that showed desirable metal uptake efficiencies in yeast, including high Mn uptake and low Cd uptake, even under various Mn and Cd concentrations that mimic natural soil conditions (Table 3 and Figure 4). However, it is indispensable to evaluate the Cd/Mn uptake characteristics in plants, because the actual Cd and Mn concentrations of the water fraction in the paddy soil are difficult to estimate, and the amino acid substitution might cause unexpected changes in the physiological or agricultural traits. The metal uptake characteristics of OsNramp5-Q337K have been evaluated in rice, and the Cd and Mn contents in the straw of the mutant rice were estimated as 40% of those in WT rice (Kuramata et al., 2022). This is a comparable value with the estimated Cd and Mn uptake efficiency in this study (40%; Figure 2). In another study, it was reported that in the analysis of amino acid substitution mutants of AtNramp1, Mn concentrations in yeast showed the same tendency as in plants results (Fu et al., 2022). These observations suggest that the characteristics of the OsNramp5 mutant can also be observed in rice plants. The obtained three mutants (M235A, A232C+M235A, and A232S+M235A) which showed Cd uptake of background level even in high concentrations of Cd (7.5 µM; Figure 4) should be the priorities for the evaluation in rice plants.
Future studies should apply a prime editor technique to incorporate M235A, A232C+M235A, or A232S+M235A mutations in planta (Anzalone et al., 2019). This technique enables precise base substitution and eventually produces a final product with no foreign genes in the genome. Several prime editing applications have been reported in various plants (Butt et al., 2020; Lin et al., 2020; Zong et al., 2022). The double mutation in this study can be managed by a single prime editing step, as shown by mutations at two positions of DNA 40 nucleotides apart (Tao et al., 2022).
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Author contributions
JI: Conceptualization, Investigation, Methodology, Writing – original draft, Writing – review & editing. TT: Data curation, Validation, Writing – review & editing. TK: Data curation, Project administration, Validation, Writing – review & editing. TN: Conceptualization, Funding acquisition, Project administration, Supervision, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was partly supported by JSPS KAKENHI (grant number 22H02611).
Acknowledgments
We would like to thank all our colleagues at Nakamura’s laboratory at Kyushu University for their kind advice regarding the experimental protocol and discussions. The ICP-OES (Agilent 5800) measurements were performed at the Center for Advanced Instrumental and Educational Support, Faculty of Agriculture, Kyushu University, Japan.
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.2024.1482099/full#supplementary-material
References
Alejandro, S., Cailliatte, R., Alcon, C., Dirick, L., Domergue, F., Correia, D., et al. (2017). Intracellular distribution of manganese by the trans-golgi network transporter NRAMP2 is critical for photosynthesis and cellular Redox homeostasis. Plant Cell 29, 3068–3084. doi: 10.1105/tpc.17.00578
Anzalone, A. V., Randolph, P. B., Davis, J. R., Sousa, A. A., Koblan, L. W., Levy, J. M., et al. (2019). Search-and-replace genome editing without double-strand breaks or donor DNA. Nature 576, 149–157. doi: 10.1038/s41586-019-1711-4
Asami, T., Kubota, M., Minamisawa, K. (1988). Natural abundance of cadmium, antimony, bismuth and some other heavy metals in Japanese soils. Soil Sci. Plant Nutr. 59, 197–199. doi: 10.20710/dojo.59.2_197
Ballabio, C., Jones, A., Montanarella, L., Toth, G. (2023). Cadmium in the Soils of the EU Analysis of LUCAS Soils data for the review of Fertilizer Directive (Luxembourg: Publications Office of the European Union). doi: 10.2760/3942
Bozzi, A. T., Bane, L. B., Weihofen, W. A., McCabe, A. L., Singharoy, A., Chipot, C. J., et al. (2016a). Conserved methionine dictates substrate preference in Nramp-family divalent metal transporters. Proc. Natl. Acad. Sci. U.S.A. 113, 10310–10315. doi: 10.1073/pnas.1607734113
Bozzi, A. T., Bane, L. B., Weihofen, W. A., Singharoy, A., Guillen, E. R., Ploegh, H. L., et al. (2016b). Crystal structure and conformational change mechanism of a bacterial nramp-family divalent metal transporter. Structure 24, 2102–2114. doi: 10.1016/j.str.2016.09.017
Bozzi, A. T., Gaudet, R. (2021). Molecular mechanism of nramp-family transition metal transport. J. Mol. Biol. 433. doi: 10.1016/j.jmb.2021.166991
Bozzi, A. T., McCabe, A. L., Barnett, B. C., Gaudet, R. (2020). Transmembrane helix 6b links proton and metal release pathways and drives conformational change in an Nrampfamily transition metal transporter. J. Biol. Chem. 295, 1212–1224. doi: 10.1074/jbc.RA119.011336
Bozzi, A. T., Zimanyi, C. M., Nicoludis, J. M., Lee, B. K., Zhang, C. H., Gaudet, R. (2019). Structures in multiple conformations reveal distinct transition metal and proton pathways in an Nramp transporter. Elife 8 (e41124). doi: 10.7554/eLife.41124.001
Butt, H., Rao, G. S., Sedeek, K., Aman, R., Kamel, R., Mahfouz, M. (2020). Engineering herbicide resistance via prime editing in rice. Plant Biotechnol. J. 18, 2370–2372. doi: 10.1111/pbi.13399
Cailliatte, R., Lapeyre, B., Briat, J. F., Mari, S., Curie, C. (2009). The NRAMP6 metal transporter contributes to cadmium toxicity. Biochem. J. 422, 217–228. doi: 10.1042/BJ20090655
Cailliatte, R., Schikora, A., Briat, J. F., Mari, S., Curie, C. (2010). High-affinity manganese uptake by the metal transporter nramp1 is essential for Arabidopsis growth in low manganese conditions. Plant Cell 22, 904–917. doi: 10.1105/tpc.109.073023
Chaloupka, R., Courville, P., Veyrier, F., Knudsen, B., Tompkins, T. A., Cellier, M. F. M. (2005). Identification of functional amino acids in the Nramp family by a combination of evolutionary analysis and biophysical studies of metal and proton cotransport in vivo. Biochemistry 44, 726–733. doi: 10.1021/bi048014v
Chang, J. D., Huang, S., Yamaji, N., Zhang, W., Ma, J. F., Zhao, F. J. (2020). OsNRAMP1 transporter contributes to cadmium and manganese uptake in rice. Plant Cell. Environ. 43, 2476–2491. doi: 10.1111/pce.13843
Chen, X., Li, J., Wang, L., Ma, G., Zhang, W. (2016). A mutagenic study identifying critical residues for the structure and function of rice manganese transporter OsMTP8.1. Sci. Rep. 6 (32073). doi: 10.1038/srep32073
Chen, Y., Zhao, X., Li, G., Kumar, S., Sun, Z., Li, Y., et al. (2021). Genome-wide identification of the nramp gene family in spirodela polyrhiza and expression analysis under cadmium stress. Int. J. Mol. Sci. 22 (6414). doi: 10.3390/ijms22126414
Clemens, S., Ma, J. F. (2016). Toxic heavy metal and metalloid accumulation in crop plants and foods. Annu. Rev. Plant Biol. 67, 489–512. doi: 10.1146/annurev-arplant-043015-112301
Curie, C., Alonso, J. M., Jean, M. L. E., Ecker, J. R., Briat, J.-F. (2000). Involvement of NRAMP1 from Arabidopsis thaliana in iron transport. Biochem. J. 347, 749–755. doi: 10.1042/bj3470749
da Silva Guerra, J. V., Ribeiro-Filho, H. V., Jara, G. E., Bortot, L. O., de Carvalho Pereira, J. G., Lopes-de-Oliveira, P. S. (2021). pyKVFinder: an efficient and integrable Python package for biomolecular cavity detection and characterization in data science. BMC Bioinf. 22 (607). doi: 10.1186/s12859-021-04519-4
Eastman, P., Swails, J., Chodera, J. D., McGibbon, R. T., Zhao, Y., Beauchamp, K. A., et al. (2017). OpenMM 7: Rapid development of high performance algorithms for molecular dynamics. PloS Comput. Biol. 13 (7), e1005659. doi: 10.1371/journal.pcbi.1005659
Edgar, R. C. (2004). MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32, 1792–1797. doi: 10.1093/nar/gkh340
EFSA (2009). Cadmium in food - Scientific opinion of the Panel on Contaminants in the Food Chain. EFSA J. 7 980, 1–139. doi: 10.2903/j.efsa.2009.980
Ehrnstorfer, I. A., Geertsma, E. R., Pardon, E., Steyaert, J., Dutzler, R. (2014). Crystal structure of a SLC11 (NRAMP) transporter reveals the basis for transition-metal ion transport. Nat. Struct. Mol. Biol. 21, 990–996. doi: 10.1038/nsmb.2904
Eisenberg, D., Weiss, R. M., Terwilliger, T. C. (1984). The hydrophobic moment detects periodicity in protein hydrophobicity. Proc. Nadl. Acad. Sci. U.S.A. 81, 140–144. doi: 10.1073/pnas.81.1.140
Fu, D., Zhang, Z., Wallrad, L., Wang, Z., Höller D, S., Ju, C., et al. (2022). Ca 2+-dependent phosphorylation of NRAMP1 by CPK21 and CPK23 facilitates manganese uptake and homeostasis in Arabidopsis. Proc. Nadl. Acad. Sci. U.S.A. 119 (40), e2204574119. doi: 10.1073/pnas.2204574119
Gao, H., Xie, W., Yang, C., Xu, J., Li, J., Wang, H., et al. (2018). NRAMP2, a trans-Golgi network-localized manganese transporter, is required for Arabidopsis root growth under manganese deficiency. New Phytol. 217, 179–193. doi: 10.1111/nph.14783
García-Esquinas, E., Pollan, M., Tellez-Plaza, M., Francesconi, K. A., Goessler, W., Guallar, E., et al. (2014). Cadmium exposure and cancer mortality in a prospective cohort: The strong heart study. Environ. Health Perspect. 122, 363–370. doi: 10.1289/ehp.1306587
Gietz, R. D., Schiestl, R. H. (2007). High-efficiency yeast transformation using the LiAc/SS carrier DNA/PEG method. Nat. Protoc. 2, 31–34. doi: 10.1038/nprot.2007.13
Gill, S. S., Tuteja, N. (2011). Cadmium stress tolerance in crop plants: Probing the role of sulfur. Plant Signal Behav. 6, 215–222. doi: 10.4161/psb.6.2.14880
Guerra, J. V. S., Ribeiro-Filho, H. V., Pereira, J. G. C., Lopes-De-Oliveira, P. S. (2023). KVFinder-web: A web-based application for detecting and characterizing biomolecular cavities. Nucleic Acids Res. 51, W289–W297. doi: 10.1093/nar/gkad324
Guo, J., Long, L., Chen, A., Dong, X., Liu, Z., Chen, L., et al. (2022). Tonoplast-localized transporter ZmNRAMP2 confers root-to-shoot translocation of manganese in maize. Plant Physiol. 190, 2601–2616. doi: 10.1093/plphys/kiac434
Hekkelman, M. L., de Vries, I., Joosten, R. P., Perrakis, A. (2023). AlphaFill: enriching AlphaFold models with ligands and cofactors. Nat. Methods 20, 205–213. doi: 10.1038/s41592-022-01685-y
Honma, T., Shiratori, Y., Ohba, H., Tsuchida, T., Makino, T., Abe, T., et al. (2017). Concentrations of nutrient content in rice variety Koshihikari Kan No.1 and risk estimation of incidence of brown spot disease in different paddy fields. JPNSS 88, 213–220. doi: 10.20710/dojo.88.3_213
Huang, H., Yamaji, N., Huang, S., Ma, J. F. (2024). Uptake and accumulation of cobalt is mediated by osNramp5 in rice. Plant Cell Environ. Advance online publication. doi: 10.1111/pce.15130
Hutchinson, T. C., Meema, K. M., Page, A. L., Chang, A. C., El-Amamy, M. (1987). Lead, Mercury, Cadmium and Arsenic in the Environment. Eds. Hutchinson, T. C., Meema, K. M. (Chichester: Published on behalf of the Scientific Committee on Problems of the Environment of the International Council of Scientific Unions by Wiley).
Ishikawa, S., Ishimaru, Y., Igura, M., Kuramata, M., Abe, T., Senoura, T., et al. (2012). Ion-beam irradiation, gene identification, and marker-assisted breeding in the development of low-cadmium rice. Proc. Natl. Acad. Sci. U.S.A. 109, 19166–19171. doi: 10.1073/pnas.1211132109
Ishimaru, Y., Takahashi, R., Bashir, K., Shimo, H., Senoura, T., Sugimoto, K., et al. (2012). Characterizing the role of rice NRAMP5 in Manganese, Iron and Cadmium Transport. Sci. Rep. 2 (286). doi: 10.1038/srep00286
Järup, L., Alfvén, T. (2004). Low level cadmium exposure, renal and bone effects - the OSCAR study. BioMetals 17, 505–509. doi: 10.1023/B:BIOM.0000045729.68774.a1
Jones, D. T., Taylor, W. R., Thornton, J. M. (1992). The rapid generation of mutation data matrices from protein sequences. Comput. Appl. Biosci. 8, 275–282. doi: 10.1093/bioinformatics/8.3.275
Julin, B., Wolk, A., Johansson, J. E., Andersson, S. O., Andrén, O., Kesson, A. (2012). Dietary cadmium exposure and prostate cancer incidence: A population-based prospective cohort study. Br. J. Cancer 107, 895–900. doi: 10.1038/bjc.2012.311
Jumper, J., Evans, R., Pritzel, A., Green, T., Figurnov, M., Ronneberger, O., et al. (2021). Highly accurate protein structure prediction with AlphaFold. Nature 596, 583–589. doi: 10.1038/s41586-021-03819-2
Jurcik, A., Bednar, D., Byska, J., Marques, S. M., Furmanova, K., Daniel, L., et al. (2018). CAVER Analyst 2.0: Analysis and visualization of channels and tunnels in protein structures and molecular dynamics trajectories. Bioinformatics 34, 3586–3588. doi: 10.1093/bioinformatics/bty386
Kubier, A., Wilkin, R. T., Pichler, T. (2019). Cadmium in soils and groundwater: A review. Appl. Geochem. 108 (104388). doi: 10.1016/j.apgeochem.2019.104388
Kuramata, M., Abe, T., Tanikawa, H., Sugimoto, K., Ishikawa, S. (2022). A weak allele of OsNRAMP5 confers moderate cadmium uptake while avoiding manganese deficiency in rice. J. Exp. Bot. 73, 6475–6489. doi: 10.1093/jxb/erac302
Li, J., Liu, Y., Kong, L., Xu, E., Zou, Y., Zhang, P., et al. (2023). An intracellular transporter OsNRAMP7 is required for distribution and accumulation of iron into rice grains. Plant Sci. 336 (111831). doi: 10.1016/j.plantsci.2023.111831
Li, H., Wang, N., Hu, W., Yan, W., Jin, X., Yu, Y., et al. (2022). ZmNRAMP4 enhances the tolerance to aluminum stress in arabidopsis thaliana. Int. J. Mol. Sci. 23 (15). doi: 10.3390/ijms23158162
Li, J., Wang, Y., Zheng, L., Li, Y., Zhou, X., Li, J., et al. (2019). The intracellular transporter atNRAMP6 is involved in fe homeostasis in arabidopsis. Front. Plant Sci. 10 (5). doi: 10.3389/fpls.2019.01124
Li, J., Wang, L., Zheng, L., Wang, Y., Chen, X., Zhang, W. (2018). A functional study identifying critical residues involving metal transport activity and selectivity in natural resistance-associated macrophage protein 3 in arabidopsis thaliana. Int. J. Mol. Sci. 19. doi: 10.3390/ijms19051430
Lin, Q., Zong, Y., Xue, C., Wang, S., Jin, S., Zhu, Z., et al. (2020). Prime genome editing in rice and wheat. Nat. Biotechnol. 38, 582–585. doi: 10.1038/s41587-020-0455-x
Lu, M., Yang, G., Li, P., Wang, Z., Fu, S., Zhang, X., et al. (2018). Bioinformatic and functional analysis of a key determinant underlying the substrate selectivity of the al transporter, Nrat1. Front. Plant Sci. 9. doi: 10.3389/fpls.2018.00606
Mirdita, M., Schütze, K., Moriwaki, Y., Heo, L., Ovchinnikov, S., Steinegger, M. (2022). ColabFold: making protein folding accessible to all. Nat. Methods 19, 679–682. doi: 10.1038/s41592-022-01488-1
Mirdita, M., Steinegger, M., Söding, J. (2019). MMseqs2 desktop and local web server app for fast, interactive sequence searches. Bioinformatics 35, 2856–2858. doi: 10.1093/bioinformatics/bty1057
Nazir, M. M., Noman, M., Ahmed, T., Ali, S., Ulhassan, Z., Zeng, F., et al. (2022). Exogenous calcium oxide nanoparticles alleviate cadmium toxicity by reducing Cd uptake and enhancing antioxidative capacity in barley seedlings. J. Hazard Mater. 438 (129498). doi: 10.1016/j.jhazmat.2022.129498
Nogawa, K., Kido, T. (1993). Envrolmnental Health Biological monitoring of cadmium exposure in itai-itai disease epidemiology. Int. Arch. Occup. Environ. Health 65, 543–546. doi: 10.1007/BF00381306
Peris-Peris, C., Serra-Cardona, A., Śanchez-Sanuy, F., Campo, S., Ariño, J., Segundo, B. S. (2017). Two NRAMP6 isoforms function as iron and manganese transporters and contribute to disease resistance in rice. Mol. Plant-Microbe Interact. 30, 385–398. doi: 10.1094/MPMI-01-17-0005-R
Podar, D., Scherer, J., Noordally, Z., Herzyk, P., Nies, D., Sanders, D. (2012). Metal selectivity determinants in a family of transition metal transporters. J. Biol. Chem. 287, 3185–3196. doi: 10.1074/jbc.M111.305649
Pottier, M., Oomen, R., Picco, C., Giraudat, J., Scholz-Starke, J., Richaud, P., et al. (2015). Identification of mutations allowing Natural Resistance Associated Macrophage Proteins (NRAMP) to discriminate against cadmium. Plant J. 83, 625–637. doi: 10.1111/tpj.12914
Qing, Y., Li, Y., Cai, X., He, W., Liu, S., Ji, Y., et al. (2023). Assessment of cadmium concentrations in foodstuffs and dietary exposure risk across China: A metadata analysis. Expo. Health 15, 951–961. doi: 10.1007/s12403-022-00530-z
Qu, Z., Nakanishi, H. (2023). Amino acid residues of the metal transporter OsNRAMP5 responsible for cadmium absorption in rice. Plants 12 (24). doi: 10.3390/plants12244182
Ray, S., Berry, S. P., Wilson, E. A., Zhang, C. H., Shekhar, M., Singharoy, A., et al. (2023). High-resolution structures with bound Mn2+ and Cd2+ map the metal import pathway in an Nramp transporter. Elife 12 (e84006). doi: 10.7554/eLife.84006
Rentsch, D., Laloi, M., Rouhara, I., Schmelzer, E., Delrot, S., Frommer, W. B. (1995). NTR1 encodes a high affinity oligopeptide transporter in Arabidopsis. FEBS Lett. 370, 264–268. doi: 10.1016/0014-5793(95)00853-2
Ritz, C., Baty, F., Streibig, J. C., Gerhard, D. (2015). Dose-response analysis using R. PloS One 10 (12). doi: 10.1371/journal.pone.0146021
Sasaki, A., Yamaji, N., Yokosho, K., Ma, J, F. (2012). Nramp5 is a major transporter responsible for manganese and cadmium uptake in rice. Plant Cell 24, 2155–2167. doi: 10.1105/tpc.112.096925
Socha, A. L., Guerinot, M.L. (2014). Mn-euvering manganese: The role of transporter gene family members in manganese uptake and mobilization in plants. Front. Plant Sci. 5. doi: 10.3389/fpls.2014.00106
Staessen, J. A., Roels, H. A., Emelianov, D., Kuznetsova, T., Thijs, L., Vangronsveld, J., et al. (1999). Environmental exposure to cadmium, forearm bone density, and risk of fractures: prospective population study. THE Lancet 353, 1140–1144. doi: 10.1016/s0140-6736(98)09356-8
Sui, F. Q., Chang, J. D., Tang, Z., Liu, W. J., Huang, X. Y., Zhao, F. J. (2018). Nramp5 expression and functionality likely explain higher cadmium uptake in rice than in wheat and maize. Plant Soil 433, 377–389. doi: 10.1007/s11104-018-3849-5
Sun, G. L., Reynolds, E. E., Belcher, A. M. (2019). Designing yeast as plant-like hyperaccumulators for heavy metals. Nat. Commun. 10 (5080). doi: 10.1038/s41467-019-13093-6
Takahashi, R., Ishimaru, Y., Senoura, T., Shimo, H., Ishikawa, S., Arao, T., et al. (2011). The OsNRAMP1 iron transporter is involved in Cd accumulation in rice. J. Exp. Bot. 62, 4843–4850. doi: 10.1093/jxb/err136
Tamura, K., Stecher, G., Kumar, S. (2021). MEGA11: molecular evolutionary genetics analysis version 11. Mol. Biol. Evol. 38, 3022–3027. doi: 10.1093/molbev/msab120
Tao, R., Wang, Y., Jiao, Y., Hu, Y., Li, L., Jiang, L., et al. (2022). Bi-PE: bi-directional priming improves CRISPR/Cas9 prime editing in mammalian cells. Nucleic Acids Res. 50, 6423–6434. doi: 10.1093/nar/gkac506
Thomine, S., Wang, R., Ward, J. M., Crawford, N. M., Schroeder, J. I. (2000). Cadmium and iron transport by members of a plant metal transporter family in Arabidopsis with homology to Nramp genes. Proc. Natl. Acad. Sci. U.S.A. 97, 4991–4996. doi: 10.1073/pnas.97.9.4991
Tim, N., Michelle, P., Janneke, H., Harry, A. R., Hilde, C., Lutgarde, T., et al. (2006). Environmental exposure to cadmium and risk of cancer: a prospective population-based study. Lancet Oncol. 7, 119–126. doi: 10.1016/S1470-2045(06)70545-9
Tiwari, M., Sharma, D., Dwivedi, S., Singh, M., Tripathi, R. D., Trivedi, P. K. (2014). Expression in Arabidopsis and cellular localization reveal involvement of rice NRAMP, OsNRAMP1, in arsenic transport and tolerance. Plant Cell Environ. 37, 140–152. doi: 10.1111/pce.12138
Vervaet, B. A., D’Haese, P. C., Verhulst, A. (2017). Environmental toxin-induced acute kidney injury. Clin. Kidney J. 10, 747–758. doi: 10.1093/ckj/sfx062
Wallin, M., Barregard, L., Sallsten, G., Lundh, T., Sundh, D., Lorentzon, M., et al. (2021). Low-level cadmium exposure is associated with decreased cortical thickness, cortical area and trabecular bone volume fraction in elderly men: The MrOS Sweden study. Bone 143, 115768. doi: 10.1016/j.bone.2020.115768
World Health Organization (2006). Evaluation of Certain Food Contaminants : Sixty-fourth Report of the Joint FAO/WHO Expert Committee on Food Additives (Geneva, Switzerland: Geneva: World Health Organization).
World Health Organization (2011). Evaluation of certain food additives and contaminants : seventy-third report of the Joint FAO/WHO Expert Committee on Food Additives (Geneva, Switzerland: World Health Organization).
World Health Organization (2022). Guidelines for drinking-water quality: Fourth edition incorporating the first and second addenda (World Health Organization).
Xia, J., Yamaji, N., Kasai, T., Ma, J. F. (2010). Plasma membrane-localized transporter for aluminum in rice. Proc. Natl. Acad. Sci. U.S.A. 107, 18381–18385. doi: 10.1073/pnas.1004949107
Yamaji, N., Sasaki, A., Xia, J. X., Yokosho, K., Ma, J. F. (2013). A node-based switch for preferential distribution of manganese in rice. Nat. Commun. 4 (2442). doi: 10.1038/ncomms3442
Yang, M., Zhang, W., Dong, H., Zhang, Y., Lv, K., Wang, D., et al. (2013). OsNRAMP3 is a vascular bundles-specific manganese transporter that is responsible for manganese distribution in rice. PloS One 8 (e83990). doi: 10.1371/journal.pone.0083990
Yang, M., Zhang, Y., Zhang, L., Hu, J., Zhang, X., Lu, K., et al. (2014). OsNRAMP5 contributes to manganese translocation and distribution in rice shoots. J. Exp. Bot. 65, 4849–4861. doi: 10.1093/jxb/eru259
Yohannan, S., Yang, D., Faham, S., Boulting, G., Whitelegge, J., Bowie, J. U. (2004). Proline substitutions are not easily accommodated in a membrane protein. J. Mol. Biol. 341, 1–6. doi: 10.1016/j.jmb.2004.06.025
Zhang, X., Chen, D., Zhong, T., Zhang, X., Cheng, M., Li, X. (2015). Assessment of cadmium (Cd) concentration in arable soil in China. Environ. Sci. pollut. Res. 22, 4932–4941. doi: 10.1007/s11356-014-3892-6
Zhao, J., Yang, W., Zhang, S., Yang, T., Liu, Q., Dong, J., et al. (2018). Genome-wide association study and candidate gene analysis of rice cadmium accumulation in grain in a diverse rice collection. Rice 11 (61). doi: 10.1186/s12284-018-0254-x
Keywords: OsNramp5, rice, cadmium, manganese, transporter, selectivity, protein engineering
Citation: Inoue J, Teramoto T, Kazama T and Nakamura T (2024) Engineering rice Nramp5 modifies cadmium and manganese uptake selectivity using yeast assay system. Front. Plant Sci. 15:1482099. doi: 10.3389/fpls.2024.1482099
Received: 17 August 2024; Accepted: 30 October 2024;
Published: 20 November 2024.
Edited by:
M. J. I. Shohag, University of Florida, United StatesReviewed by:
Ute – Krämer, Ruhr University Bochum, GermanyToshiro Shigaki, The University of Tokyo, Japan
Copyright © 2024 Inoue, Teramoto, Kazama and Nakamura. 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: Takahiro Nakamura, tnaka@agr.kyushu-u.ac.jp