Understanding the gene regulatory basis of plant response to heavy metals (HMs) is fundamental for the management of food safety and security. However, a comprehensive and comparative view of the plant responses to different HMs is still lacking. Here, we compared root transcriptomes in common bean under 9 HM treatments at 50 μM for three time points each. Cd, Cr, Co, Ni, and Pb caused most severe morphological and/or biochemical retardations. A total of 448 genes were found to be responsive to all nine HMs, which were mostly involved in photosynthesis, oxidization-reduction, and ion binding. Cd and Cu triggered the greatest number of unique differentially expressed genes (DEG)s, which were predominantly related to cellular transport/localization in the former and RNA binding in the latter. Short-term and prolonged HM treatments shaped very different DEG patterns. Weighted gene co-expression network analysis identified six co-expression modules showing exceptionally high transcripts abundance in specific HM × time scenarios. We experimentally verified the promoter activity of the gene GIP1 and the novel function of XTH23 under Cu/Cd stress. Collectively, the transcriptomic atlas provides valuable resources for better understanding the common and unique mechanisms of plant response to different HMs and offers a mass of candidate target genes/promoters for genetic engineering.
Keywords: common bean; hairy root; heavy metal; promoter; transcriptome.