RNA-Seq Analysis of Rice Roots Reveals the Involvement of Post-Transcriptional Regulation in Response to Cadmium Stress

Widely-spread cadmium (Cd) pollution in the soil threatens both crop production and human health. How plants deal with the excess Cd are largely unknown. To evaluate the molecular mechanism by which plants respond to Cd stress, rice seedlings were treated with two concentrations of Cd and subjected to deep RNA sequencing. Comprehensive RNA-Seq analysis of rice roots under two gradients of Cd treatment revealed 1169 Cd toxicity-responsive genes. These genes were involved in the reactive oxygen species scavenging system, stress response, cell wall formation, ion transport, and signal transduction. Nine out of 93 predicted long non-coding RNAs (lncRNAs) were detected as Cd-responsive lncRNAs due to their high correlation with the Cd stress response. In addition, we analyzed alternative splicing (AS) events under different Cd concentrations. Four hundred and seventy-six differential alternatively spliced genes with 542 aberrant splicing events were identified. GO analysis indicated that these genes were highly enriched in oxidation reduction and cellular response to chemical stimulus. Real-time qRT-PCR validation analysis strengthened the reliability of our RNA-Seq results. The results suggest that post-transcriptional AS regulation may also be involved in plant responses to high Cd stress.