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.

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

A Novel Family of PLAC8 Motif-Containing/pcr Genes Mediates Cd Tolerance and Cd Accumulation in Rice

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

Comparative Transcriptomics Reveals the Key Pathways and Genes of Cadmium Accumulation in the High Cadmium-Accumulating Rice (oryza Sativa L.) Line

Comparison of Transcriptome Differences between Two Rice Cultivars Differing in Cadmium Translocation from Spike-Neck to Grain

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

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

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

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

Genome-wide Identification and Analysis of Japonica and Indica Cultivar-Preferred Transcripts in Rice Using 983 Affymetrix Array Data.

Cadmium Absorption in Various Genotypes of Rice under Cadmium Stress

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

The Relationship between Cadmium-Related Gene Sequence Variations in Rice and Cadmium Accumulation

Genotypic Difference of Germination and Early Seedling Growth in Response to Cd Stress and Its Relation to Cd Accumulation

Time course analysis of gene regulation under cadmium stress in rice

Loci and Natural Alleles for Cadmium-Mediated Growth Responses Revealed by a Genome Wide Association Study and Transcriptome Analysis in Rice

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

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

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

Root-induced Changes to Cadmium Speciation in the Rhizosphere of Two Rice (oryza Sativa L.) Genotypes