MicroRNA390 Is Involved in Cadmium Tolerance and Accumulation in Rice

Cadmium (Cd) is a non-essential heavy metal that is toxic to plants. microRNAs (miRNAs) are 21-nucleotide RNAs that are ubiquitous regulators of gene expression at the post-transcriptional level. Several plant miRNAs, such as miR390, have vital roles in plant growth, development and responses to environmental stresses including heavy metal stress. In this study, the expression of mature miR390 was significantly down-regulated under Cd stress in rice. Consequently, the target gene of miR390, OsSRK was dramatically induced by Cd treatment. Transgenic rice plants overexpressing miR390 displayed reduced Cd tolerance and higher Cd accumulation compared with wild-type plants. Simultaneously, expression of OsSRK was less pronounced in 35S:MIR390 plants than in wild-type. These results indicate that miR390 was a negative regulator involved in Cd stress tolerance in rice.

Introduction

Cadmium is an important environmental pollutant that is highly toxic to plants (Schützendübel et al., 2001). It is widely released into soil, air, and water mainly by effluent from industrial sources, mining, and phosphate fertilization (Pinto et al., 2004). Cd is readily taken up by plants, resulting in toxic symptoms. It damages the photosynthetic apparatus, decreases carbon assimilation and chlorophyll content, and eventually leads to wilting and plant death (Hsu and Kao, 2003; Rodriguez-Serrano et al., 2009). Plants have evolved various molecular mechanisms to respond to heavy metal stress, one of which involves microRNA-guided gene regulation at the post-transcriptional level (Huang et al., 2008, 2010; Zhou et al., 2008, 2012). miRNAs are a class of endogenous non-coding small RNAs that can base pair their target mRNAs to induce their degradation or repress their translation in eukaryotic organisms (Bartel, 2004). In plants, increasing reports have demonstrated that miRNAs have vital regulatory roles in growth, development and plant resistance to abiotic and biotic stresses, including heavy metal stress (Sunkar et al., 2006; Ding and Zhu, 2009; Khraiwesh et al., 2012).

microRNA390 (miR390) is conserved in rice (Oryza sativa), maize (Zea mays), and Arabidopsis thaliana. In rice, it is represented by one member with one locus, whereas in Arabidopsis and poplar (Populus trichocarpa), it is represented by two members with three and six loci, respectively (Sunkar et al., 2005). In Arabidopsis, miR390 exerts its action through the biogenesis of ta-siRNAs that, in turn, leads to the degradation of ARFs (auxin response factors) that play critical roles in lateral root development (Marin et al., 2010; Meng et al., 2010). Recent studies indicate that it also plays an important role in plant stress tolerance because it was up-regulated under drought stress in cowpea (Vigna unguiculata; Barrera-Figueroa et al., 2011) and Brachypodium distachyon (Budak and Akpinar, 2011), while down-regulated under aluminum (Al) treatment in Medicago truncatula (Chen et al., 2012). The rice genome encodes miR390, which regulates the expression of O. sativa stress-responsive leucine-rich repeat receptor-like kinase (OsSRK; LOC_Os02g10100). Cleavage of OsSRK mRNA by miR390 has been confirmed by 5′ rapid amplification of cDNA ends (5′ RACE; Sunkar et al., 2005). RLKs belong to the large RLK/Pelle gene family, and are involved in a range of processes, including regulation of development, disease resistance, and biotic and/or abiotic stress responses (Osakabe et al., 2005; Ouyang et al., 2010). However, the role of miR390 in the rice stress response and the biological consequences of target gene regulation of miR390 remain to be identified.

In this study, a transgenic approach was used to investigate the role of miR390 in the Cd stress response in rice. Whereas, mature miR390 was found to be down-regulated by Cd stress in rice, OsSRK transcripts were dramatically induced under Cd stress. Overexpression of miR390 significantly retarded seedling growth under Cd stress conditions compared with wild-type plants. MiR390-overexpressing plants showed higher accumulation of Cd and a higher degree of Cd-induced oxidative stress than wild-type. These results provided a link between miR390 and Cd stress tolerance in rice.

Materials and Methods

Plant Treatment and RNA Isolation

Seeds of rice Zhonghua 11 (O. sativa L. subsp. japonica) were germinated in an incubator at 37°C in the dark. Uniformly germinated seeds were then transferred to a pot under a 13 h light (29°C)/11 h dark (22°C) photoperiod. For Cd treatment, 2-weeks-old seedlings were transferred to plastic containers containing Yoshida nutrient solution with 60 μM CdCl₂, and the roots were harvested at 1, 3, 6, 12, 36, and 48 h post-treatment. Untreated seedlings were used as controls. The total RNA from roots was extracted using the Trizol reagent (Invitrogen).

Real-Time Quantitative PCR Analysis

Total RNA extracted from rice roots after exposure to 60 μM CdCl₂ for 0, 1, 3, 6, 12, 36, and 48 h was separately isolated using Trizol reagent (Invitrogen) and treated with 5 U of RNase-free DNase I (TaKaRa) to remove DNA contamination from the RNA. The differentially expressed mature miR390 were validated by stem–loop RT-PCR as previously reported (Chen et al., 2005; Varkonyi-Gasic et al., 2007). Real-time PCRs were carried out using a Rotor-Gene Q machine (parameters: 95°C for 1 min, followed by 45 cycles of 95°C for 10 s, 58°C for 15 s, and 72°C for 15 s). Experiments were performed in triplicate and the results were represented by mean ± SE of three replicates. Normalized expression levels were calculated with the 2^-ΔΔC(t) method using U6 as the internal reference gene. Primers used in stem–loop RT-PCR and real-time PCR are listed as follows: miR390 RT primer, 5′-GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACGGCGTA-3′; forward, 5′-GCCGAAGCTCAGGAGG-3′; reverse, 5′-GTGCAGGGTCCGAGGT-3′. U6 RT primer, 5′-ATTTGGACCATTTCTCGATTTGT-3′; forward, 5′-CGATAAAATTGGAACGATACAGA-3′; reverse, 5′-ATTTGGACCATTTCATTTGT-3′. OsACTB forward, 5′-GCCGTCCTCTCTCTGTATGC-3′; reverse, 5′-GGGGACAGTGTGGCTGAC-3′. OsSRK (LOC_Os02g10100) forward, 5′-CCTTCGCAAACTTCTCCG-3′; reverse, 5′-ACTGCCTCCTTCTCAACA-3′.

Sequence Alignment and cis-acting Element Analysis of miR390 and its Target

The pre-miR390 sequence was downloaded from miRBase, miRBase Release 18.0¹. The 1500-bp DNA sequences upstream of the transcription initiation site of the pre-miR390 gene and OsSRK were extracted from the TIGR Rice Genome Annotation site². These sequences were then checked by using PlantCARE³ (Lescot et al., 2001).

Generation and Identification of Transgenic Rice Overexpressing miR390

The fragment of miR390 precursor was amplified from genomic DNA with the following primer pairs: forward, 5′-CGGGGTACCGGAGAGATGTTTTGAGGAAGG-3′; reverse, 5′-AACTGCAGCAGATTTAATTGGTCGTGTGG-3′. The amplified fragments were introduced into p1301-35S-NOS between 35S promoter and NOS terminator with KpnI and PstI enzymes. Cultivars of japonica rice Zhonghua 11 were used for transformation. Transgenic lines were achieved by co-cultivation of rice calli with Agrobacterium tumefaciens strain EHA105 containing p1301-35S:MIR390 (Hiei et al., 1994). All the putative T0 transgenic plants were screened using PCR analysis with genomic DNA from their leaves. The PCR was performed with primers hpt-F and hpt-R in a 20-μL reaction system. Transgenic rice seedlings were selected using 50 mg/L hyg. Expression levels of mature miR390, OsSRK, and metal transporter genes in transgenic rice plants were examined using real-time PCR. OsHMA2 forward, 5′-TGGCACCACAAAAGGCTATT-3′; reverse, 5′-CACCGTCGATTGGAATGACT-3′; OsNRAMP5 forward, 5′-AGTGGTTACAGGGAGGCATC-3′, reverse, 5′-GTCTTCCTCGATAGCACCAAG-3′. Quantification of gene expression was done using the comparative CT method. Experiments were performed in triplicate and the results were represented by mean ± SE of three replicates. GUS staining was also used to confirm the transgenic plants (Jefferson, 1987). Histochemical assay of GUS staining was performed by incubating the transgenic plants in reaction mixture comprising: 50 mM phosphate buffer (pH: 7.0) containing 1 mM 5-bromo-4-chloro-3-indolyl glucuronide (X-gluc), 5% methanol, 10 μg/mL cycloheximide, and 1 mM dithiothreitol. The reaction was done overnight at 37°C, followed by clearing with 70% ethanol.

Effect of Cd Stress on Growth of Rice Seedlings

For the phenotypic analysis of transgenic rice, T1 transgenic rice plants were selected using 50 mg/L hyg. Three weeks after germination, transgenic and wild-type rice seedlings showing consistent growth were transferred to 100 and 200 μM CdCl₂ Yoshida nutrient solution. The phenotypes of the plants were examined and photographed after Cd treatment.

For physiological and biochemical analysis of transgenic rice, T1 transgenic rice seedlings were selected using 50 mg/L hyg and then cultured with Yoshida’s culture solution for 6 weeks. After exposure to 150 μM CdCl₂ for 14 days, the chlorophyll, malonyldialdehyde (MDA), and hydrogen peroxide (H₂O₂) content was determined. All experiments were performed in triplicate. Data points were presented as the mean ± standard deviation (SD) of three replications. Significant differences among wild-type and transgenic lines were analyzed using the Student’s t-test.

The chlorophyll content was measured as described by Knudson et al. (1977). About 100 mg of rice leaves excised to 2–3 cm in length were immersed in 2 ml 95% ethanol for 36 h in the dark. The absorbance of the extracts was read at 645 and 663 nm. The total chlorophyll content was then calculated.

The level of lipid peroxidation in the rice leaves was determined in terms of the peroxidation byproduct MDA following the method of Wang et al. (2009). About 0.1 g of plant leaves was homogenized in 1 ml of 10% trichloroacetic acid (TCA). After the homogenate was centrifuged at 12,000 × g for 10 min at 4°C, 2 ml of the supernatant was added 2 ml of 0.6% thiobarbituric acid (TBA) solution. The mixtures were heated at 95°C for 30 min and then cooled quickly on ice. The resulting mixtures were centrifuged at 10,000 × g for 10 min and the absorbance of the supernatants was measured at 450, 532, and 600 nm.

The H₂O₂ content was determined following the method of Jana and Choudhuri (1982). About 0.1 g of leaf segments were ground into powder and homogenized in 3 mL of pre-cooled 50 mM phosphate buffer (pH 6.5). After centrifugation at 6000 × g for 25 min, 1 mL of 0.1% titanium sulfate in 20% (v/v) H₂SO₄ was added to 3 mL of the resultant suspension. The mixture was immediately centrifuged at 6000 × g for 15 min. Absorbance of the supernatant was recorded at 410 nm. The content of H₂O₂ was confirmed using the extinction coefficient of 0.28 μM^-1 cm^-1, and expressed as μmol g^-1 FW.

Cd Content Determination

For determination of Cd, roots of rice seedlings after 150 μM CdCl₂ treatment for 14 days were immersed in 20 mM disodium ethylenediamine tetra-acetic acid (Na₂-EDTA) for 20 min and then rinsed three times with deionized water (He et al., 2008). The roots, stems (containing the leaf sheath) and leaves were dried at 105°C for 2 h and then at 70°C to achieve a constant weight. 0.1 g DW of plant material was digested with a mixture of HNO₃: HF (6:2, v/v) at 120°C with a Multiwave until they were completely digested. Cd was then quantified using an atomic absorption spectrometer (AA-7000, Shimadzu, Tokyo, Japan). The amount of Cd was expressed on the basis of DW.

Results

Cd Stress Decreases miR390 Level

Rice miRNA expression patterns were examined under Cd stress using a microarray assay containing probes complementary to the miRNA sequences in miRBase Release 11.0⁴ (Ding et al., 2011). 1-week-old rice seedlings were treated with 60 μM CdCl₂ for 6 h, after which small RNA was isolated from roots and used to screen the microarray. Microarray results showed that miR390 expression was substantially down-regulated under Cd stress (Ding et al., 2011). In this study, we performed real time PCR analysis of the expression of mature miR390 in 2-weeks-old rice seedlings. Mature miR390 expression was down-regulated after 6 h exposure to 60 μM CdCl₂, which was consistent with the microarray data (Figure 1B).

FIGURE 1

FIGURE 1. Expression levels of miR390 and its target gene OsSRK as detected by qPCR. (A) Sequence alignment of miR390 and OsSRK. The matched base pairs are indicated by a vertical line. (B) Detection of the transcript of mature miR390 expression from 2-weeks-old rice seedlings treated with 60 μM CdCl₂ for 6 h by qPCR. U6 gene was used as the inner control. (C) QPCR analysis of the mature miR390 and OsSRK transcript levels in roots of 2-weeks-old rice seedlings treated for the indicated times with 60 μM CdCl₂. Quantifications were normalized to the expression of OsACTB. Error bars represent the SE for three independent experiments.

To examine time-dependent changes in the expression of mature miR390, qPCR analysis was performed in 2-weeks-old rice seedlings. In rice roots, miR390 expression decreased after 1 h of treatment with 60 μM Cd, and continued to decrease over 36 h of the same treatment (Figure 1C). By contrast, the target gene of miR390, OsSRK, was induced sharply after exposure to Cd for 36 h (Figure 1C). These results showed that the profiles of miR390 and OsSRK transcripts were complementary, but not exactly opposite, to each other.

Sequence Analysis of miR390 and OsSRK

MicroRNA sequences were derived from miRBase⁵. Sequence alignment of miR390 and OsSRK was as shown in Figure 1A. Upstream cis-acting regulatory element sequences of miR390 and OsSRK were analyzed by plantCARE⁶. Stress- and defense-related elements, such as ARE, TC-rich repeats and MBS, were enriched in these regions, suggesting that miR390 and its target OsSRK were possibly involved in the stress responses in plants (Table 1).

TABLE 1

TABLE 1. Stress-related cis-elements analysis of miR390 and its target OsSRK.

35S:MIR390 Plants are More Sensitive to Cd Stress than are Wild-Type Plants

To determine whether miR390 influences plant growth under heavy metal conditions, transgenic rice plants overexpressing miR390 under control of the CaMV 35S promoter (35S:MIR390 plants) were generated. The hygromycin (hyg) phosphotransferase gene (hpt) was inserted in T-DNA as a selectable marker (Figure 2). The transgenic 35S:MIR390 rice plants were single-locus T-DNA insertion lines. T0 transgenic plants were verified by qPCR and GUS activity assay (Figure 3). Transcript levels of the mature miR390 were substantially increased in transgenic plants, indicating that 35S:MIR390 plants overexpress miR390. OsSRK transcript levels were lower in the transgenic plants than in wild-type controls. These results further validated the negative regulation of miR390 in OsSRK transcripts.

FIGURE 2

FIGURE 2. Expression in plants of the vector of p1301-35S:MIR390. 35S, CaMV 35S promoter; NOS, nopaline synthase terminator; GUS, gene encoding β-glucuronidase; Hpt, gene encoding hygromycin phosphotransferase.

FIGURE 3

FIGURE 3. Confirmation of 35S:MIR390 transgenic rice plants. (A) Genomic PCR analysis of transgenic rice plants. Genomic DNA was used as the template for PCR. The expected size of the product was 621 bp. Key: M, DNA marker; WT, wild-type plant; L1–L5, transgenic plants. (B) GUS activity analysis of transgenic rice plants. (C) qPCR analysis of the mature miR390 and OsSRK levels in wild-type and transgenic rice.

We examined the growth of miR390-overexpressing rice seedlings under Cd stressed conditions (Figures 4A and 5A). When the plants were grown under normal conditions, no noticeable differences in seedling growth were observed between wild-type and transgenic plants. However, the plant growth of transgenic lines was dramatically reduced compared with that of wild-type plants under Cd stress (Figure 4A). Quantitative analyses confirmed that the shoot weight, lateral root weight and root length of 35S:MIR390 plants were significantly lower than those of wild-type plants when grown on Cd (Figures 4B–E). The chlorophyll content was also determined. Under normal conditions, no significant differences in chlorophyll content were observed among the plants. However, 35S:MIR390 lines had less chlorophyll under Cd stress compared with wild-type controls (Figure 5B).

FIGURE 4

FIGURE 4. Decreased Cd resistance of 35S:MIR390 transgenic rice plants. (A) Phenotypic comparison of wild-type and 35S:MIR390 transgenic rice plants under Cd stress. In each growth chamber, WT plants were grown on the left and 35S:MIR390 plants on the right. Upper panel: control, plants grown for 3 weeks under normal conditions before treatment. Middle panel: 3-weeks-old plants treated with 100 μM CdCl₂ (left) and 200 μM CdCl₂ (right) for 20 days. Lower panel: 3-weeks-old plants treated with 100 μM CdCl₂ (left) and 200 μM CdCl₂ (right) for 40 days. (B,C) Shoot fresh weight (FW) and root FW of wild-type and 35S:MIR390 transgenic plants grown under Cd stress for 20 days. (D,E) Shoot FW and root length of wild-type and 35S:MIR390 transgenic plants grown under Cd stress for 40 days. In (B–E), significant differences from the wild-type as determined by the Student’s t-test are indicated (^∗P < 0.05, ^∗∗P < 0.01). Bars indicate SE.

FIGURE 5

FIGURE 5. Performance of wild-type and 35S:MIR390 transgenic plants treated with 150 μM CdCl₂ for 14 days. (A) Phenotypic comparison of WT and 35S:MIR390 transgenic rice. (B) Chlorophyll content. (C) MDA content. (D) H₂O₂ content. Each column represents an average of three replicates, and bars indicate SE. ^∗Indicate significant differences from between wild-type and transgenic plants at P < 0.05.

Malonyldialdehyde is an indicator of oxidative attack on membrane lipids. When subjected to Cd stress, MDA levels increased in wild-type rice seedlings, which attributed to the Cd-induced oxidative stress, in accordance with previous experiments (Shah et al., 2001; Guo et al., 2004). However, in Cd-stressed 35S:MIR390 lines, MDA content increased by 40% compared with wild-type plants (Figure 5C). These results indicated that overexpression of miR390 can increase membrane damage caused by Cd stress.

Stresses usually cause damage in plants via oxidative stress involving the generation of reactive oxygen species (ROS), such as H₂O₂ (Apel and Hirt, 2004). Therefore, the H₂O₂ content in leaves of wild-type and transgenic plants were examined. As shown in Figure 5D, there were no significant differences in H₂O₂ content under normal conditions. After 14 days of 150 μM CdCl₂ treatment, the H₂O₂ content in leaves of wild-type increased by 1.4-fold, whereas those in transgenic plants increased by 2.0-fold, suggesting that overexpression of miR390 cannot efficiently eliminate H₂O₂ produced under Cd stress, resulting in the decreased Cd tolerance of transgenic rice. All these results indicated that the 35S:MIR390 plants were Cd sensitive.

35S:MIR390 Plants Accumulate more Cd than Wild-Type Plants

To assess whether the increase in Cd sensitivity was because of higher metal accumulation, the Cd content in roots and shoots of miR390-overexpressing rice lines was compared with that of the same organs of wild-type plants. The differences in Cd concentration of roots were not significant between wild-type and transgenic lines, whereas the Cd levels were higher in the leaves and stems of transgenic lines than in wild-type plants (Figure 6). Thus, compared with wild-type plants, 35S:MIR390 plants enhanced Cd accumulation, which might result in the sensitivity of 35S:MIR390 plants to Cd stress.

FIGURE 6

FIGURE 6. The contents of Cd in leaves (A), stems (containing leaf sheaths) (B) and roots (C) of wild-type and transgenic plants exposed to 150 mM CdCl₂ for 14 days. Data are presented as the mean ± SD; ^∗∗ indicates significant differences at P < 0.01 between wild-type and transgenic plants.

To estimate the effect of miR390 on Cd accumulation and translocation to shoots, several heavy metal transporter genes were examined in 35S:MIR390 roots after Cd treatment by qPCR (Figure 7). In our analyses, more transcripts of OsHMA2 and OsNRAMP5 were produced in the Cd-stressed roots of 35S:MIR390 than in wild-type plants, indicating that they were positively regulated by miR390 and important for Cd accumulation and translocation in rice.

FIGURE 7

FIGURE 7. Quantitative real-time PCR (QPCR) analysis of metal transporter genes in 35S:MIR390 transgenic plant roots compare with wild-type plants under Cd stress. (A) OsNRAMP5. (B) OsHMA2. Normalized expression levels were calculated with the 2^-ΔΔC(t) method using OsACTB as the internal reference gene.

Discussion

Characterization of the role of miRNAs in the plant stress response is an active field of research. Recently, several methods have been used for the dissection and characterization of diverse miRNAs involved in plant responses and adaptations to heavy metal stress (Sunkar and Zhu, 2004; Zhou et al., 2008; Chen et al., 2012). For example, miRNA microarray assays led to the identification of 19 Cd-responsive miRNAs in rice (Ding et al., 2011). Deep sequencing studies revealed that 52 new miRNA candidates were differentially regulated by mercury Hg(II) in M. truncatula (Zhou et al., 2012). Although, an increasing number of heavy metal-related miRNAs has been identified, their precise roles remain to be verified. Overexpression of miRNAs and/or miRNA targets or knockout mutants of target genes have proved useful for analyzing the function of stress-responsive miRNAs and further uncovering miRNA-regulated transcripts upon application of stress. Transgenic strategies have also been used to improve the heavy metal stress tolerance of plants. Zhang et al. (2013) reported that miR395 overexpressing rapeseed (Brassica napus) showed a lower degree of Cd-induced oxidative stress than wild-type plants, suggesting that miR395 would be involved in detoxification of Cd in B. napus. This study was the first to report that overexpression of a miRNA can affect heavy metal tolerance of rice.

Previous studies suggested that the miR390 family was involved in the plant heavy metal stress response (Mendoza-Soto et al., 2012; Yang and Chen, 2013). miR390 has been reported to be down-regulated by Al, As, Cd, and Hg (Ding et al., 2011; Chen et al., 2012; Zhou et al., 2012). In this study, mature miR390 was decreased in roots of rice under Cd stress, which was consistent to the expression pattern of miR390 precursor in our previous study (Ding et al., 2011). The transgenic plants overexpressing miR390 showed enhanced sensitivity to Cd stress. MiR390-overexpressing plants had higher accumulation of Cd and higher translocation of Cd from roots to shoots relative to the wild-type. Furthermore, expression of miR390 changed the expression pattern of several metal transporter genes. Results showed that OsHMA2 and OsNRAMP5 were increased in miR390-overexpressing transgenic plants compared to wild-type plants under Cd treatment. OsHMA2 was reported to be a major Zn/Cd transporter from roots to shoots (Takahashi et al., 2012). In this study, the expression level of OsHMA2 was increased in 35S:MIR390 plants, in which Cd contents were higher in the shoots compared with wild-type plants. Thus, the increase in Cd translocation of 35S:MIR390 plants could be attributable to the change in the expression level of OsHMA2. In addition, OsNRAMP5, was reported to be a Mn and Cd transporter involved in the root uptake of these metals from the medium (Ishimaru et al., 2012). More OsNRAMP5 transcripts were detected in 35S:MIR390 roots, implying more Cd accumulation in rice.

Plant miRNAs execute biological functions via negative regulation of their specific target genes. In Arabidopsis, miR390 was demonstrated to mediate the regulatory pathway miR390-TAS3-ARFs involved in auxin signaling and the regulation of lateral root development (Marin et al., 2010; Yoon et al., 2010). In this study, we investigated lateral root growth in 35S:MIR390 transgenic lines in response to Cd stress. Quantitative analyses showed that lateral root weight and root length of 35S:MIR390 plants were significantly lower than those of wild-type plants when grown on Cd. However, no reduction in lateral root number was observed in 35S:MIR390 seedlings treated with Cd. These results were partly consistent with Yoon et al.’s (2010) report in which miR390 exerted a negative effect on lateral root formation in Arabidopsis. The role of miR390-TAS3-ARFs pathway on lateral root development in rice remains to be validated and elucidated.

In rice, miR390 has been experimentally validated to target a stress-responsive LRR-like kinase (OsSRK; LOC_Os02g10100) by 5′ RACE and degradome sequencing (Sunkar et al., 2005; Zhou et al., 2010). RLKs are transmembrane proteins containing an extracellular LRR and a Ser/Thr kinase domain, and can activate a complex array of intracellular signaling pathways in response to various plant environmental and developmental signals (Tichtinsky et al., 2003; Osakabe et al., 2005; Chae et al., 2009). RLKs include multigene families. To date, at least 1131 RLK genes have been identified in rice, several of which have been reported to play important roles in plant stress response (Shiu et al., 2004). Transgenic rice plants with overexpression of a putative RLK gene, OsSIK1 showed higher tolerance to salt and drought stresses than wild-type plants (Ouyang et al., 2010). In germinating rice seeds exposed to excess copper (Cu), a RLK protein spot was identified to be up-regulated by MALDI-TOF mass spectrometry (Ahsan et al., 2007). A receptor-like protein gene, OsRMC was reported to be involved in regulation of Fe acquisition in rice (Yang et al., 2013). In this study, the expression of OsSRK (Os02g10100) was up-regulated after exposure to Cd stress in rice seedlings. OsSRK was repressed in 35S:MIR390 plants, in which the expression levels of several metal transporter genes changed. However, the relationship between OsSRK and metal transporter genes needs further research. Gain- and loss-of-function plants targeting OsSRK would be helpful to determine the role of OsSRK in regulation of plant response to heavy metal stress.

Author Contributions

YD and CZ conceived and designed experiment; YY, ZJ, and YW carried out experiments; YD, YY, ZJ, and YW analyzed experimental results; YD and YY wrote the manuscript; All authors read and approved the manuscript.

Conflict of Interest Statement

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.

Acknowledgments

This work was supported by grants from the National Natural Science Foundation of China (31401299, 31170251, 31470368), and Zhejiang Provincial Natural Science Foundation of China (LZ14C020001, LY13C020002).

Abbreviations

CaMV, cauliflower mosaic virus; Cd, cadmium; miRNA, microRNA; qPCR, quantitative real-time PCR; RLK, receptor-like kinase; RT-PCR, reverse transcription PCR.

Footnotes

1. ^http://microrna.sanger.ac.uk/
2. ^http://www.mirbase.org/cgi-bin/mirna_entry.pl?acc=MI0001690
3. ^http://bioinformatics.psb.ugent.be/webtools/plantcare/html/
4. ^http://www.tigr.org/tdb/e2k1/osa1/data_download.shtml
5. ^http://microrna.sanger.ac.uk/sequences/
6. ^http://bioinformatics.psb.ugent.be/webtools/plantcare/html/

References

Ahsan, N., Lee, D. G., Lee, S. H., Kang, K. Y., Lee, J. J., Kim, P. J., et al. (2007). Excess copper induced physiological and proteomic changes in germinating rice seeds. Chemosphere 67, 1182–1193. doi: 10.1016/j.chemosphere.2006.10.075

PubMed Abstract | CrossRef Full Text | Google Scholar

Apel, K., and Hirt, H. (2004). Reactive oxygen species: metabolism, oxidative stress, and signal transduction. Annu. Rev. Plant Biol. 55, 373–399. doi: 10.1146/annurev.arplant.55.031903.141701

PubMed Abstract | CrossRef Full Text | Google Scholar

Barrera-Figueroa, B. E., Gao, L., Diop, N. N., Wu, Z. G., Ehlers, J. D., and Roberts, P. A. (2011). Identification and comparative analysis of drought-associated microRNAs in two cowpea genotypes. BMC Plant Biol. 11:127. doi: 10.1186/1471-2229-11-127

PubMed Abstract | CrossRef Full Text | Google Scholar

Bartel, D. P. (2004). MicroRNAs: genomics, biogenesis, mechanism, and function. Cell 116, 281–297. doi: 10.1016/S0092-8674(04)00045-5

CrossRef Full Text | Google Scholar

Budak, H., and Akpinar, A. (2011). Dehydration stress-responsive miRNA in Brachypodium distachyon: evident by genome-wide screening of microRNAs expression. OMICS 15, 791–799. doi: 10.1089/omi.2011.0073

PubMed Abstract | CrossRef Full Text | Google Scholar

Chae, L., Sudat, S., Dudoit, S., Zhu, T., and Luan, S. (2009). Diverse transcriptional programs associated with environmental stress and hormones in the Arabidopsis receptor-like kinase gene family. Mol. Plant 2, 84–107. doi: 10.1093/mp/ssn083

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, C., Ridzon, D. A., Broomer, A. J., Zhou, Z., Lee, D. H., Nguyen, J. T., et al. (2005). Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res. 33:e179. doi: 10.1093/nar/gni178

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, L., Wang, T. Z., Zhao, M. G., Tian, Q. Y., and Zhang, W. H. (2012). Identification of aluminum-responsive microRNAs in Medicago truncatula by genome-wide high-throughput sequencing. Planta 235, 375–386. doi: 10.1007/s00425-011-1514-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Ding, Y. F., Chen, Z., and Zhu, C. (2011). Microarray-based analysis of cadmium-responsive microRNAs in rice (Oryza sativa). J. Exp. Bot. 62, 3563–3573. doi: 10.1093/jxb/err046

PubMed Abstract | CrossRef Full Text | Google Scholar

Ding, Y. F., and Zhu, C. (2009). The role of microRNAs in copper and cadmium homeostasis. Biochem. Biophy. Res. Commun. 386, 6–10. doi: 10.1016/j.bbrc.2009.05.137

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, T., Zhang, G., Zhou, M., Wu, F., and Chen, J. (2004). Effects of aluminum and cadmium toxicity on growth and antioxidant enzyme activities of two barley genotypes with different Al resistance. Plant Soil 258, 241–248. doi: 10.1023/B:PLSO.0000016554.87519.d6

CrossRef Full Text | Google Scholar

He, J. Y., Zhu, C., Ren, Y. F., Yan, Y. P., Cheng, C., and Jiang, D. A. (2008). Uptake, subcellular distribution, and chemical forms of cadmium in wild-type and mutant rice. Pedosphere 18, 371–377. doi: 10.1016/S1002-0160(08)60027-2

CrossRef Full Text | Google Scholar

Hiei, Y., Ohta, S., Komari, T., and Kumashiro, T. (1994). Efficient transformation of rice (Oryza sativa L.) mediated by Agrobacterium and sequence analysis of the boundaries of the T-DNA. Plant J. 6, 271–282. doi: 10.1046/j.1365-313X.1994.6020271.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Hsu, Y. T., and Kao, C. H. (2003). Role of abscisic acid in cadmium tolerance of rice (Oryza sativa L.) seedlings. Plant Cell Environ. 26, 867–874. doi: 10.1046/j.1365-3040.2003.01018.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, S. Q., Peng, J., Qiu, C. X., and Yang, Z. M. (2008). Heavy metal-regulated new microRNAs from rice. J. Inorg. Biochem. 103, 282–287. doi: 10.1016/j.jinorgbio.2008.10.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, S. Q., Xiang, A. L., Che, L. L., Chen, S., Li, H., and Song, J. B. (2010). A set of miRNAs from Brassica napus in response to sulfate-deficiency and cadmium stress. Plant Biotechnol. J. 8, 887–899. doi: 10.1111/j.1467-7652.2010.00517.x

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Jana, S., and Choudhuri, M. A. (1982). Glycolate metabolism of three submerged aquatic angiosperm during aging. Aquat. Bot. 12, 345–354. doi: 10.1016/0304-3770(82)90026-2

CrossRef Full Text | Google Scholar

Jefferson, R. A. (1987). Assaying chimeric genes in plants: the GUS gene fusion system. Plant Mol. Biol. Rep. 5, 387–405. doi: 10.1007/BF02667740

CrossRef Full Text | Google Scholar

Khraiwesh, B., Zhu, J. K., and Zhu, J. H. (2012). Role of miRNAs and siRNAs in biotic and abiotic stress responses of plants. Biochim. Biophys. Acta 1819, 137–148. doi: 10.1016/j.bbagrm.2011.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Knudson, L. L., Tibbitts, T. W., and Edwards, G. E. (1977). Measurement of ozone injury by determination of leaf chlorophyll concentration. Plant Physiol. 60, 606–608. doi: 10.1104/pp.60.4.606

CrossRef Full Text | Google Scholar

Lescot, M., Déhais, P., Thijs, G., Marchal, K., Moreau, Y., and Van de Peer, Y. (2001). PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 30, 325–327. doi: 10.1093/nar/30.1.325

PubMed Abstract | CrossRef Full Text | Google Scholar

Marin, E., Jouannet, V., Herz, A., Lokerse, A. S., Weijers, D., and Vaucheret, H. (2010). miR390, Arabidopsis TAS3 tasiRNAs, and their AUXIN RESPONSE FACTOR targets define an autoregulatory network quantitatively regulating lateral root growth. Plant Cell 22, 1104–1117. doi: 10.1105/tpc.109.072553

PubMed Abstract | CrossRef Full Text | Google Scholar

Mendoza-Soto, A. B., Sánchez, F., and Hernández, G. (2012). MicroRNAs as regulators in plant metal toxicity response. Front. Plant Sci. 3:105. doi: 10.3389/fpls.2012.00105

PubMed Abstract | CrossRef Full Text | Google Scholar

Meng, Y., Ma, X., Chen, D., Wu, P., and Chen, M. (2010). MicroRNA-mediated signaling involved in plant root development. Biochem. Biophys. Res. Commun. 393, 345–349. doi: 10.1016/j.bbrc.2010.01.129

PubMed Abstract | CrossRef Full Text | Google Scholar

Osakabe, Y., Maruyama, K., Seki, M., Satou, M., Shinozaki, K., and Yamaguchi-Shinozaki, K. (2005). Leucine-rich repeat receptor-like kinase1 is a key membrane-bound regulator of abscisic acid early signaling in Arabidopsis. Plant Cell 17, 1105–1119. doi: 10.1105/tpc.104.027474

PubMed Abstract | CrossRef Full Text | Google Scholar

Ouyang, S. Q., Liu, Y. F., Liu, P., Lei, G., He, S. J., and Ma, B. (2010). Receptor-like kinase OsSIK1 improves drought and salt stress tolerance in rice (Oryza sativa) plants. Plant J. 62, 316–329. doi: 10.1111/j.1365-313X.2010.04146.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Pinto, A. P., Mota, A. M., De Varennes, A., and Pinto, F. C. (2004). Influence of organic matter on the uptake of cadmium, zinc, copper and iron by sorghum plants. Sci. Total Environ. 326, 239–274. doi: 10.1016/j.scitotenv.2004.01.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Rodriguez-Serrano, M., Romero-Puertas, M. C., Pazmino, D. M., Testillano, P. S., Risueno, M. C., and Rio, L. A. (2009). Cellular response of pea plants to cadmium toxicity: cross talk between reactive oxygen species, nitric oxide, and calcium. Plant Physiol. 150, 229–243. doi: 10.1104/pp.108.131524

PubMed Abstract | CrossRef Full Text | Google Scholar

Schützendübel, A., Schwanz, P., Teichmann, T., Gross, K., Langenfeld-Heyser, R., Godbold, D. L., et al. (2001). Cadmium-induced changes in antioxidative systems, hydrogen peroxide content, and differentiation in Scots pine roots. Plant Physiol. 127, 887–898. doi: 10.1104/pp.010318

PubMed Abstract | CrossRef Full Text | Google Scholar

Shah, K., Kumar, R. G., Verma, S., and Dubey, R. S. (2001). Effect of cadmium on lipid peroxidation, superoxide anion generation and activities of antioxidant enzymes in growing rice seedlings. Plant Sci. 161, 1135–1144. doi: 10.1016/S0168-9452(01)00517-9

CrossRef Full Text | Google Scholar

Shiu, S. H., Karlowski, W. M., Pan, R., Tzeng, Y. H., Mayer, K. F., and Li, W. H. (2004). Comparative analysis of the receptor-like kinase family in Arabidopsis and rice. Plant Cell 16, 1220–1234. doi: 10.1105/tpc.020834

PubMed Abstract | CrossRef Full Text | Google Scholar

Sunkar, R., Girke, T., Jain, P. K., and Zhu, J. K. (2005). Cloning and characterization of microRNAs from rice. Plant Cell 17, 1397–1411. doi: 10.1105/tpc.105.031682

PubMed Abstract | CrossRef Full Text | Google Scholar

Sunkar, R., Kapoor, A., and Zhu, J. K. (2006). Posttranscriptional induction of two Cu/Zn superoxide dismutase genes in Arabidopsis is mediated by down-regulation of miR398 and important for oxidative stress tolerance. Plant Cell 18, 2051–2065. doi: 10.1105/tpc.106.041673

PubMed Abstract | CrossRef Full Text | Google Scholar

Sunkar, R., and Zhu, J. K. (2004). Novel and stress-regulated microRNAs and other small RNAs from Arabidopsis. Plant Cell 16, 2001–2019. doi: 10.1105/tpc.104.022830

PubMed Abstract | CrossRef Full Text | Google Scholar

Takahashi, H., Watanabe, A., and Akagi, H. (2012). Mutations in rice (Oryza sativa) Heavy Metal ATPase 2 (OsHMA2) restrict the translocation of zinc and cadmium. Plant Cell Physiol. 53, 213–224. doi: 10.1093/pcp/pcr166

PubMed Abstract | CrossRef Full Text | Google Scholar

Tichtinsky, G., Vanoosthuyse, V., Cock, J. M., and Gaude, T. (2003). Making inroads into plant receptor kinase signalling pathways. Trends Plant Sci. 8, 231–237. doi: 10.1016/S1360-1385(03)00062-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Varkonyi-Gasic, E., Wu, R., Wood, M., Walton, E. F., and Hellens, R. P. (2007). Protocol: a highly sensitive RT-PCR method for detection and quantification of microRNAs. Plant Methods 3, 12. doi: 10.1186/1746-4811-3-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, F. J., Zeng, B., Sun, Z. X., and Zhu, C. (2009). Relationship between proline and Hg²⁺-induced oxidative stress in a tolerant rice mutant. Arch. Environ. Contam. Toxicol. 56, 723–731. doi: 10.1007/s00244-008-9226-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, A., Li, Y., Xu, Y., and Zhang, W. H. (2013). A receptor-like protein RMC is involved in regulation of iron acquisition in rice. J. Exp. Bot. 64, 5009–5020. doi: 10.1093/jxb/ert290

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Z. M., and Chen, J. (2013). A potential role of microRNAs in plant response to metal toxicity. Metallomics 5, 1184–1190. doi: 10.1039/c3mt00022b

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoon, E. K., Yang, J. H., Lim, J., Kim, S. H., Kim, S. K., and Lee, W. S. (2010). Auxin regulation of the microRNA390-dependent transacting small interfering RNA pathway in Arabidopsis lateral root development. Nucleic Acids Res. 38, 1382–1391. doi: 10.1093/nar/gkp1128

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, L. W., Song, J. B., Xx, S., Zhang, Y., and Yang, Z. M. (2013). miR395 is involved in detoxification of cadmium in Brassica napus. J. Hazard. Mater. 250, 204–211. doi: 10.1016/j.jhazmat.2013.01.053

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, M., Gui, L. F., Li, P. C., Song, X. W., Wen, L. Y., and Chen, Z. Y. (2010). Degradome sequencing reveals endogenous small RNA targets in rice (Oryza sativa L. ssp. indica). Front. Biol. 5:90. doi: 10.1007/s11515-010-0007-8

CrossRef Full Text | Google Scholar

Zhou, Z. S., Huang, S. Q., and Yang, Z. M. (2008). Bioinformatic identification and expression analysis of new microRNAs from Medicago truncatula. Biochem. Biophys. Res. Commun. 374, 538–542. doi: 10.1016/j.bbrc.2008.07.083

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Z. S., Zeng, H. Q., Liu, Z. P., and Yang, Z. M. (2012). Genome-wide identification of Medicago truncatula microRNAs and their targets reveals their differential regulation by heavy metal. Plant Cell Environ. 35, 86–99. doi: 10.1111/j.1365-3040.2011.02418.x

PubMed Abstract | CrossRef Full Text | Google Scholar