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

To date, Cd remains a major contaminant in rice production. An in-depth exploration of the mechanism that causes genotypic differences in Cd enrichment in rice is necessary to develop strategies to regulate Cd enrichment in rice. Here, two rice cultivars (low grain Cd, ZZ143; and high grain Cd, YX409) displayed different transcriptomic profile patterns when subjected to 100μmol/L Cd stress. In fact, 18,721(9833 upregulated and 8888 downregulated) and 16,403 (8366 upregulated and 8037 downregulated) differentially expressed genes (DEGs) were found in ZZ143 and YX409, respectively. Gene ontology (GO) classification revealed 28 and 26 terms enriched in ZZ143 and YX409, respectively. ZZ143 had more enriched DEGs than YX409, primarily in cellular processes, metabolic processes, binding terms, catalytic activity, etc. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that 21 and 24 pathways were significantly enriched in ZZ143 and YX409, respectively. Based on the DEGs, ZZ143 had a stronger ability for sulfur assimilation and Cys synthesis, whereas YX409 had a stronger ability to maintain cell wall stability. A series of DEGs involved in metabolic pathways, biosynthesis of secondary metabolites, plant hormone signal transduction pathways, and mitogen-activated protein kinase signaling pathways were identified, which maybe closely related to Cd resistance and the different Cd concentrations between cultivars. The above pathways and the greater number of identified DEGs in more than half of the GO terms and KEGG pathways suggest a higher absorption and stronger tolerance of the roots of ZZ143 than YX409 to Cd.