Comparative transcriptome analysis between two different cadmium-accumulating genotypes of soybean (Glycine max) in response to cadmium stress

Cadmium (Cd) is extremely toxic and non-essential for plants. Different soybean varieties differ greatly in their Cd accumulation ability, but little is known about the underlying molecular mechanisms. Here, we performed transcriptomic analysis using Illumina pair-end sequencing on root tissues from two soybean varieties (su8, high-Cd-accumulating (HAS) and su7, low Cd-accumulating (LAS)) grown with 0 or 50 μM CdSO4. A total of 18.76 million clean reads from the soybean root samples were obtained after quality assessment and data filtering. After Cd treatment, 739 differentially expressed genes (DEGs; 265 up and 474 down) were found in HAS; however, only 259 DEGs (88 up and 171 down) were found in LAS, and 64 genes were same between the two varieties. Pathway enrichment analysis suggested that after cadmium treatment, the DEGs between LAS and HAS were mainly enriched in glutathione metabolism and plant-pathogen interaction pathways. KEGG analysis showed that phenylalanine metabolism responding to cadmium stress in LAS, while ABC transporters responding to cadmium stress in HAS. Besides we found more differential expressed heavy metal transporters such as ABC transporters and zinc transporters in HAS than LAS, and there were more transcription factors differently expressed in HAS than LAS after cadmium treatment in two soybean varieties, eg. bHLH transcription factor, WRKY transcription factor and ZIP transcription factor. Findings from this study will shed new insights on the underlying molecular mechanisms behind the Cd accumulation in soybean.