Transcriptomic and metabolomic analyses reveal several critical metabolic pathways and candidate genes involved in resin biosynthesis in Pinus massoniana

Pine resin, which typically consists of terpenoids, is a natural product used in various industrial applications. Oleoresin can be obtained from the xylem tissue by wounding the stem bark. Pinus massoniana (masson pine) is an important resin-tapping tree species that originated in southern China. Masson pines with different genetic backgrounds typically have different resin-yielding capacities (RYCs). However, the mechanisms underlying high resin yield in masson pines are unclear. The aim of this study was to identify the possible genetic regulation pathways and functional genes that influence the resin yield. In this study, we conducted transcriptomic and metabolomic studies of masson pine secondary xylem with high, medium, and low RYCs. A total of 230,068 unigenes and 3894 metabolites were identified from the tissue of the secondary xylem. Several differentially expressed regulation factors, including WRKY, bHLH, and ERF, and functional genes such as PKc and LRR-RLKs, were identified among these masson pines. The Kyoto Encyclopedia of Genes and Genomes pathways were mainly focused on diterpenoid biosynthesis, plant hormone signal transduction, and ABC transporters. Furthermore, integration of the transcriptomic and metabolomic data indicated that the PKc- and LRR-RLK-related regulatory and metabolic pathways may play critical roles in the biosynthesis of terpenoids. These above results improve our understanding of the biosynthesis mechanism of oleoresin in P. massoniana and facilitate further research work into the functional analysis of these candidate genes.