Zhenggang Xu1, Meng Dong2, Xiaoyun Peng2, Wenzhen Ku2, Yunlin Zhao3, Guiyan Yang4. 1. Hunan Research Center of Engineering Technology for Utilization of Environmental and Resources Plant, Central South University of Forestry and Technology, 498 Shaoshan South Road, 410004 Changsha, Hunan Province, China; School of Material and Chemical Engineering, Hunan City University, 518 Yingbin Road, 413000 Yiyang, Hunan Province, China. 2. School of Material and Chemical Engineering, Hunan City University, 518 Yingbin Road, 413000 Yiyang, Hunan Province, China. 3. Hunan Research Center of Engineering Technology for Utilization of Environmental and Resources Plant, Central South University of Forestry and Technology, 498 Shaoshan South Road, 410004 Changsha, Hunan Province, China. Electronic address: zyl8291290@163.com. 4. Laboratory of Walnut Research Center, College of Forestry, Northwest A & F University, Yangling 712100 Shaanxi, China. Electronic address: yangguiyan@nwsuaf.edu.cn.
Abstract
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.
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.