Abstract
Forests are crucial sustainable sources of natural ecosystems and contribute to human welfare. Cunninghamia lanceolata (Chinese fir) is an economically important conifer and occupies the largest area in China that produces global wood resources. Although Chinese fir has high economic value in China, little information is known regarding its mechanisms of wood formation. Therefore, transcriptome analysis was conducted to study the gene expression patterns and associated timber formation mechanisms in Chinese fir at different stand ages. In the present study, a total of 837,156 unigenes were identified in 84 samples from Chinese fir (pith and root) at different stand ages via RNA-Seq. Among them, most of the differentially expressed genes (DEGs) were significantly enrichment in plant hormone signal transduction, flavonoid metabolism pathway, starch and sucrose metabolism, and MAPK signal transduction pathway, which might be associated with the diameter formation in Chinese fir. The DEGs in these pathways were analyzed in Chinese fir and were related to lignin synthesis, cell wall formation and cell wall reinforcement/thickening. These genes might play an important role in regulating timber formation/growth in Chinese fir. In addition, certain transcriptome factors (TFs) related to Chinese fir timber formation were identified, including WRKY33, WRKY22, PYR/PYL, and MYC2. Weighted co-expression network analysis (WGCNA) showed that glucan endo-1,3-beta-d-glucosidase was a hub gene significantly correlated with the growth-related genes in Chinese fir. Sixteen key genes that related to diameter regulation in Chinese fir were verified by qRT-PCR analysis. These key genes might have a fine regulatory role in timber formation in Chinese fir. Our results pave the way for research on the regulatory mechanisms of wood formation, and provide an insight for improving the quality production of Chinese fir.