Comparative transcriptome analysis of the hyperaccumulator plant Phytolacca americana in response to cadmium stress

To study the molecular mechanism of the hyperaccumulator plant Phytolacca americana against cadmium (Cd) stress, the leaves of P. americana treated with 400 μM Cd for 0, 2, 12, and 24 h were harvested for comparative transcriptome analysis. In total, 110.07 Gb of clean data were obtained, and 63,957 unigenes were acquired after being assembled. Due to the lack of P. americana genome information, only 24,517 unigenes were annotated by public databases. After Cd treatment, 5054 differentially expressed genes (DEGs) were identified. KEGG pathway enrichment analysis of DEGs showed that genes involved in the flavonoid biosynthesis and antenna proteins of photosynthesis were significantly down-regulated, while genes related to the lignin biosynthesis pathway were remarkably up-regulated, indicating that P. americana could synthesize more lignin to cope with Cd stress. Moreover, genes related to heavy metal accumulation, sulfur metabolism and glutathione metabolism were also significantly up-regulated. The gene expression pattern of several key genes related to distinct metabolic pathways was verified by qRT-PCR. The results indicated that the immobilization of lignin in cell wall, chelation, vacuolar compartmentalization, as well as the increase of thiol compounds content may be the important mechanisms of Cd detoxification in hyperaccumulator plant P. americana.