A meta-analysis of transcriptomic profiles reveals molecular pathways response to cadmium stress of Gramineae

As a non-essential heavy metal, cadmium (Cd) is toxic to plants. In the last 15 years, over 70 transcriptome studies have been published to decipher the molecular response mechanism against Cd stress in different plants. To extract generalization results from transcriptomic data across different plants and obtain some hub genes that respond to Cd stress, we carried out a meta-analysis of 32 published datasets. Cluster analysis revealed that plant species played a more decisive role than the media used and exposure time in the transcriptome patterns of plant roots response to Cd. The datasets from a Gramineae-like (GL) group were closer in clustering. 838 DEGs were commonly Cd-regulated in at least nine of 18 GL datasets. Gene ontology and KEGG pathway analyses revealed that oxidative stress-related terms and lignin synthesis-related terms were significantly enriched. Mapman analysis revealed that these common DEGs were mainly involved in regulation, cellular response, secondary metabolism, transport, cell wall and lipid metabolism. In Oryza sativa, 15 DEGs were up-regulated in at least four of five HM (As, Cr, Cd, Hg and Pb) groups, such as Os10g0517500 (methionine gamma-lyase) and Os01g0159800 (bHLH107). Moreover, our datasets can be used to retrieve log2FC value of specific genes across 29 studies (48 datasets), which provides data reference for the subsequent selection of HM-related genes. Our results provide the basis for further understanding of Cd tolerance mechanisms in plants.