Comparative analysis of Cd-responsive maize and rice transcriptomes highlights Cd co-modulated orthologs

Background: Metal tolerance is often an integrative result of metal uptake and distribution, which are fine-tuned by a network of signaling cascades and metal transporters. Thus, with the goal of advancing the molecular understanding of such metal homeostatic mechanisms, comparative RNAseq-based transcriptome analysis was conducted to dissect differentially expressed genes (DEGs) in maize roots exposed to cadmium (Cd) stress. Results: To unveil conserved Cd-responsive genes in cereal plants, the obtained 5166 maize DEGs were compared with 2567 Cd-regulated orthologs in rice roots, and this comparison generated 880 universal Cd-responsive orthologs groups composed of 1074 maize DEGs and 981 rice counterparts. More importantly, most of the orthologous DEGs showed coordinated expression pattern between Cd-treated maize and rice, and these include one large orthologs group of pleiotropic drug resistance (PDR)-type ABC transporters, two clusters of amino acid transporters, and 3 blocks of multidrug and toxic compound extrusion (MATE) efflux family transporters, and 3 clusters of heavy metal-associated domain (HMAD) isoprenylated plant proteins (HIPPs), as well as all 4 groups of zinc/iron regulated transporter protein (ZIPs). Additionally, several blocks of tandem maize paralogs, such as germin-like proteins (GLPs), phenylalanine ammonia-lyases (PALs) and several enzymes involved in JA biosynthesis, displayed consistent co-expression pattern under Cd stress. Out of the 1074 maize DEGs, approximately 30 maize Cd-responsive genes such as ZmHIPP27, stress-responsive NAC transcription factor (ZmSNAC1) and 9-cis-epoxycarotenoid dioxygenase (NCED, vp14) were also common stress-responsive genes reported to be uniformly regulated by multiple abiotic stresses. Moreover, the aforementioned three promising Cd-upregulated genes with rice counterparts were identified to be novel Cd-responsive genes in maize. Meanwhile, one maize glutamate decarboxylase (ZmGAD1) with Cd co-modulated rice ortholog was selected for further analysis of Cd tolerance via heterologous expression, and the results suggest that ZmGAD1 can confer Cd tolerance in yeast and tobacco leaves. Conclusions: These novel findings revealed the conserved function of Cd-responsive orthologs and paralogs, which would be valuable for elucidating the genetic basis of the plant response to Cd stress and unraveling Cd tolerance genes.