New insight into the molecular basis of cadmium stress responses of wild paper mulberry plant by transcriptome analysis

Background: Heavy metal contamination is becoming a limitation to the utilization of soil and the distribution of vegetation. In particular, cadmium (Cd) pollution has had a serious impact on the food chain. Broussonetia papyrifera is a widely distributed pioneer tree species of heavy metal contaminated areas with important economic value. However, little is known about the genomic background of the Cd-tolerance mechanism in B. papyrifera. Results: The CdCl2 responsive physiology was evaluated and proved to be involved in antioxidase activity and active oxygen species (ROS) accumulation. The leaf and root transcriptomes derived from B. papyrifera grown under normal and CdCl2 stress conditions were systematically investigated using the Illumina HiSeq method. A total of 180,678,660 bp (27.1 GB) clean reads were assembled into 589,487 high-quality unigenes, of which 256,025 (43.43% of the total) and 250,251 (42.45% of the total) were aligned in Gene Ontology (GO) and Protein family (Pfam), respectively. A total of 24,414 differentially expressed genes (DEGs) were GO-annotated into 53, 23, 55, and 60 terms from the transcriptomes of root and leaf tissues under Cd stress and control conditions. A total of 117,547 Kyoto Encyclopedia of Genes and Genomes (KEGG) Orthology (KO)-annotated DEGs were enriched in at least 47 KEGG pathway terms among the four comparisons. Many genes encoding important transcription factors (e.g., auxin/indole-3-acetic acid (AUX/IAA), basic helix-loop-helix (bHLH), DNA-binding one zinc finger (Dof), and MYB) and proteins involved in plant-pathogen interactions, phenylpropanoid biosynthesis, plant hormone signal transduction, oxidative phosphorylation, carbon fixation, peroxisomes, flavonoid biosynthesis, and glutathione metabolism, among others, were substantially upregulated under CdCl2 stress. Conclusions: These genes represent important candidates for studying Cd-response mechanisms and molecular biology of B. papyrifera and related species. Our findings provide a genomic sequence resource for functional genetic assignments in B. papyrifera, which will help elucidate the molecular mechanisms of its Cd-stress responses and facilitate the bioremediation of heavy metal contaminated areas via breeding of new stress-tolerant cultivars.