Response of Chinese fir seedlings to low phosphorus stress and analysis of gene expression differences

Abstract

Chinese fir (Cunninghamia lanceolata) is an excellent fast-growing timber species occurring in southern China and has significant value in the forestry industry. In order to enhance the phosphorus utilization efficiency in Chinese fir, four clones named X6, S3, S39 and FK were used, and low phosphorus (LP) stress experiments were performed to analyze the response of different clones to phosphorus deficiency. According to the results on seedling height, maximum root length, leaf blade aspect ratio, root ratio, malondialdehyde content, acid phosphates activity, proline content, soluble protein level, and chlorophyll a and b levels of the tested clones, compared to the control groups (CK), the phosphorus high efficiency clone X6 was screen out for transcriptome sequencing experiments. De novo RNA-seq was then used to sequence the root transcriptomes of X6 under LP stress and CK, and we then compared the gene expression differences under the two conditions. A total of 3416 SDEGs were obtained by comparing the LP and CK groups, among which 1742 were up-regulated and 1682 were down-regulated. All SDEGs obtained from the LP and CK treated samples were subjected to KEGG annotation and classification. Through classification statistical analysis using WEGO software, 607 SDEGs obtained KEGG pathway annotations, which were related to 206 metabolic pathways. In Chinese fir subjected to LP stress, 53 SDEGs related with phosphorus metabolism, and phosphate uptake and transport were obtained from our transcriptome data. Based on the phosphorus metabolism pathway obtained by KEGG classification, combined with previously report on gene annotation related with phosphorus metabolism, the enzymes encoded by SDEG related with phosphorus metabolism and their expression pattern were mapped onto phosphorus metabolism pathway.

Keywords

Chinese fir / Low phosphorus stress / Root transcriptomes / SDEG / Phosphorus metabolism

References

[1] Audic S, Claverie J-M. The significance of digital gene expression profiles. Genome Res, 1997, 7: 986-995.

[2] Benjamini Y, Yekutieli D. The control of the false discovery rate in multiple testing under dependency. Ann Stat, 2001, 29: 1165-1188.

[3] Bian L, Zheng R, Su S, Lin H, Xiao H, Wu HX, Shi J. Spatial analysis increases efficiency of progeny testing of chinese fir. J For Res, 2017, 28(3): 445-452.

[4] Bozzo GG, Raghothama KG, Plaxton WC. Purification and characterization of two secreted purple acid phosphatase isozymes from phosphate-starved tomato (Lycopersicon esculentum) cell cultures. Eur J Biochem, 2003, 269: 6278-6286.

[5] Chen FJ, Liu XS, Guo-Hua MI. Varietal differences in plant growth, phosphorus uptake and yield formation in two maize inbred lines grown under field conditions. J Integr Agric, 2012, 11: 1738-1743.

[6] Coelho GTCP, Carneiro NP, Karthikeyan AS, Raghothama KG, Schaffert RE, Brandão RL, Paiva LV, Souza IRP, Alves VM, Imolesi A. A phosphate transporter promoter from Arabidopsis thaliana AtPHT1;4 gene drives preferential gene expression in transgenic maize roots under phosphorus starvation. Plant Mol Biol Report, 2010, 28: 717-723.

[7] Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics, 2005, 21: 3674-3676.

[8] Duan H, Hu D, Li Y, Zheng H. Characterization of a collection of chinese fir elite genotypes using sequence-related amplified polymorphism markers. J For Res, 2016, 27(5): 1-6.

[9] Dubey NC, Tripathi BP, Müller M, Stamm M, Ionov L. Enhanced activity of acetyl CoA synthetase adsorbed on smart microgel: an implication for precursor biosynthesis. ACS Appl Mater Interface, 2015, 7: 1500-1507.

[10] Fan F, Ding G, Wen X. Proteomic analyses provide new insights into the responses of Pinus massoniana seedlings to phosphorus deficiency. Proteomics, 2016, 16(3): 504-515.

[11] Fort F, Cruz P, Stroia C, Jouany C, Catrice O, Delbrut A, Luzarreta M. Root functional trait syndromes and plasticity drive the ability of grassland Fabaceae to tolerate water and phosphorus shortage. Environ Exp Bot, 2015, 110: 62-72.

[12] Gilbert GA, Knight JD, Vance CP, Allan DL. Acid phosphatase activity in phosphorus-deficient white lupin roots. Plant Cell Environ, 1999, 22(7): 801-810.

[13] Guo B, Irigoyen S, Fowler TB, Versaw WK. Differential expression and phylogenetic analysis suggest specialization of plastid-localized members of the PHT4 phosphate transporter family for photosynthetic and heterotrophic tissues. Plant Signal Behav, 2008, 3: 784-790.

[14] Hammond JP, Broadley MR, Bowen HC, Spracklen WP, Hayden RM, White PJ. Gene expression changes in phosphorus deficient potato (Solanum teubrosum L.) leaves and the potential for diagnostic gene expression markers. PLoS ONE, 2011 6 9 e24606

[15] He P, Qin H, Wu M, Wu B, Wei J, Wang D. Identification of genes differentially expressed in the roots of rubber tree (Hevea brasiliensis Muell. Arg.) in response to phosphorus deficiency. Mol Biol Rep, 2013, 40: 1397-1405.

[16] Lapis-Gaza HR, Jost R, Finnegan PM. Arabidopsis PHOSPHATE TRANSPORTER1 genes PHT1;8 and PHT1;9 are involved in root-to-shoot translocation of orthophosphate. BMC Plant Biol, 2014, 14: 1-19.

[17] Li L, Liu C, Lian X. Gene expression profiles in rice roots under low phosphorus stress. Plant Mol Biol, 2010, 72: 423-432.

[18] Li S, Ha SJ, Kim HJ, Galazka JM, Cate JHD, Jin YS, Zhao H. Investigation of the functional role of aldose 1-epimerase in engineered cellobiose utilization. J Biotechnol, 2013, 168: 1-6.

[19] Lu Yong-quan, Jia Qing, Tong Zai-kang. Cloning and sequence analysis of nine novel MYB genes in Taxodiaceae plants. J For Res, 2014, 25(4): 795-804.

[20] Ma Z, Huang B, Xu S, Chen Y, Cao G, Ding G, Lin S. Ion flux in roots of Chinese fir (Cunninghamia lanceolata (Lamb.) Hook) under aluminum stress. PLoS ONE, 2016, 11: e0156832.

[21] Marschner H. Mineral nutrition of higher plants, 1995 2 London: Academic Press.

[22] Miguel MA, Widrig A, Vieira RF, Brown KM, Lynch JP. Basal root whorl number: a modulator of phosphorus acquisition in common bean (Phaseolus vulgaris). Ann Bot, 2013, 112: 973-982.

[23] Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods, 2008, 5: 621-628.

[24] Orwa C, Mutua A, Kindt R, Jamnadass R, Simons A (2015) Agroforestree database: a tree reference and selection guide version 4.0. 2009. http://www.worldagroforestry.org/af/treedb/. Accessed 15 Feb 2011

[25] Rausch C, Zimmermann P, Amrhein N, Bucher M. Expression analysis suggests novel roles for the plastidic phosphate transporter Pht2; 1 in auto- and heterotrophic tissues in potato and Arabidopsis. Plant J, 2004, 39: 13-28.

[26] Richardson AE, Lynch JP, Ryan PR, Delhaize E, Smith FA, Smith SE, Harvey PR, Ryan MH, Veneklaas EJ, Lambers H. Plant and microbial strategies to improve the phosphorus efficiency of agriculture. Plant Soil, 2011, 349: 121-156.

[27] Rouached H, Arpat AB, Poirier Y. Regulation of phosphate starvation responses in plants: signaling players and cross-talks. Mol Plant, 2010, 3: 288-299.

[28] Saldanha AJ. Java Treeview—extensible visualization of microarray data. Bioinformatics, 2004, 20: 3246-3248.

[29] Sharma V, Kumar A, Archana G, Kumar GN. Ensifermeliloti overexpressing Escherichia coli phytase gene (appA) improves phosphorus (P) acquisition in maize plants. Sci Nat, 2016, 103: 76.

[30] Simpson RJ, Oberson A, Culvenor RA, Ryan MH, Veneklaas EJ, Lambers H, Lynch JP, Ryan PR, Delhaize E, Smith FA. Strategies and agronomic interventions to improve the phosphorus-use efficiency of farming systems. Plant Soil, 2011, 349: 89-120.

[31] Vance CP, Uhde-Stone C, Allan DL. Phosphorus acquisition and use: critical adaptations by plants for securing a nonrenewable resource. New Phytol, 2003, 157: 423-447.

[32] Wu P, Ma L, Hou X, Wang M, Wu Y, Liu F, Deng XW. Phosphate starvation triggers distinct alterations of genome expression in Arabidopsis roots and leaves. Plant Physiol, 2003, 132: 1260-1271.

[33] Yan W, Chen GH, Yang LF, Gai JY, Zhu YL. Overexpression of the rice phosphate transporter gene OsPT6 enhances tolerance to low phosphorus stress in vegetable soybean. Sci Hortic, 2014, 177: 71-76.

[34] Ye J, Fang L, Zheng H, Zhang Y, Chen J, Zhang Z, Wang J, Li S, Li R, Bolund L. WEGO: a web tool for plotting GO annotations. Nucleic Acids Res, 2006, 34: W293-W297.