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 screened 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 then gene expression differences were compared 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 the 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.