Abstract:Background The migration of cadmium (Cd) from contaminated soil to rice is a cause for concern. However, the molecular mechanism underlying the response of rice roots to various Cd stresses remains to be clarified from the viewpoint of the co-expression network at a system-wide scale. Results We employed a comparative RNAseq-based approach to identify early Cd-responsive differentially expressed genes (DEGs) in rice ‘Nipponbare’ seedling roots after 1 h of high-Cd treatment. A multiplicity of the identified 1772 DEGs were implicated in hormone signaling and transcriptional regulation, particularly NACs and WRKYs were all upregulated under Cd stress. All of the 6 Cd-upregulated ABC transporters were pleiotropic drug resistance proteins (PDRs), whereas all of the 6 ZRT/IRT-like proteins (ZIPs) were consistently downregulated by Cd treatment. To further confirm our results of this early transcriptomic response to Cd exposure, we then conducted weighted gene co-expression network analysis (WGCNA) to re-analyze our RNAseq data in combination with other 11 previously published RNAseq datasets for rice roots exposed to diverse concentrations of Cd for extended treatment periods. This integrative approach identified 271 transcripts as universal Cd-regulated DEGs that are key components of the Cd treatment coupled co-expression module. A global view of the 164 transcripts with annotated functions in pathway networks revealed several Cd-upregulated key functional genes, including transporter ABCG36/OsPDR9, 12-oxo-phytodienoic acid reductases (OPRs) for JA synthesis, and ZIM domain proteins JAZs in JA signaling, as well as OsWRKY10 , NAC, and ZFP transcription factors. More importantly, 104 of these, including ABCG36/OsPDR9 , OsNAC3 , as well as several orthologs in group metalloendoproteinase, plastocyanin-like domain containing proteins and pectin methylesterase inhibitor, may respond specifically to various Cd pressures, after subtracting the 60 general stress-responsive genes reported to be commonly upregulated following multiple stresses. Conclusion An integrative approach was implemented to identify DEGs and co-expression network modules in response to various Cd pressures, and 104 of the 164 annotatable universal Cd-responsive DEGs may specifically respond to various Cd pressures. These results provide insight into the universal molecular mechanisms beneath the Cd response in rice roots, and suggest many promising targets for improving the rice acclimation process against Cd toxicity.

Transcriptome Profiling of Abiotic Stress-Responsive Genes During Cadmium Chloride-Mediated Stress in Two Indica Rice Varieties

Transcriptome analysis of rice (Oryza sativa L.) shoots responsive to cadmium stress

Time course analysis of gene regulation under cadmium stress in rice

Multi-omics analysis reveals differential molecular responses to cadmium toxicity in rice root tip and mature zone

Co-expression Network Analysis of the Transcriptomes of Rice Roots Exposed to Various Cadmium Stresses Reveals Universal Cadmium-Responsive Genes

A transcriptomic (RNA-seq) analysis of genes responsive to both cadmium and arsenic stress in rice root.

Transcriptomic and physiological comparsion of the short-term responses of two Oryza sativa L. varieties to cadmium

Unique Transcriptome and Gene Expression Analysis of Rice Seedling Reveals Different Cadmium Response Regulatory Mechanisms Between Indica and Japonica Rice

Cadmium Absorption in Various Genotypes of Rice under Cadmium Stress

Comparative Analysis of Cd-responsive Maize and Rice Transcriptomes Highlights Cd Co-Modulated Orthologs.

Comparative transcriptomic analysis reveals the important process in two rice cultivars with differences in cadmium accumulation

Comparative transcriptome analysis reveals key genes and coordinated mechanisms in two rice cultivars differing in cadmium accumulation

Integrated morphological, physiological and transcriptomic analyses reveal response mechanisms of rice under different cadmium exposure routes

Identification of key genes and modules in response to Cadmium stress in different rice varieties and stem nodes by weighted gene co-expression network analysis

Comprehensive physiology and proteomics analysis revealed the resistance mechanism of rice (Oryza sativa L) to cadmium stress

Cadmium-zinc Cross-Talk Delineates Toxicity Tolerance in Rice Via Differential Genes Expression and Physiological / Ultrastructural Adjustments.

Phenomic and Metabolomic Responses of Roots to Cadmium Reveal Contrasting Resistance Strategies in Two Rice Cultivars (oryza Sativa L.)

Physiological Responses Involved in Cadmium Tolerance in a High-Cadmium-accumulating Rice (oryza Sativa L.) Line

Cadmium tolerance is associated with the root-driven coordination of cadmium sequestration, iron regulation, and ROS scavenging in rice

Profiling of rice Cd-tolerant genes through yeast-based cDNA library survival screening

A transcriptomic view of cadmium retention in roots of cadmium-safe rice line (Oryza sativa L.)