Identification of the diterpenoid biosynthesis genes and their expression status in relation to oleoresin yield of masson pine

Abstract

Conifer oleoresin is a highly precious raw material for industrial sectors and a precursor for biofuels. Diterpenoid, particularly diterpenoid resin acids (DRAs), is one of the most abundant oleoresin components. Masson pine (Pinus massoniana L.) is one of the most important native species for oleoresin yield in China. To better understand the biosynthesis of diterpenoid in high oleoresin-yield masson pine and develop methods for its production in needle and trunk xylem, we used the combination of SMRT full-length transcriptome and Illumina RNA sequencing technologies to establish a transcriptome dataset. Approximately 20,560 and 17,950 full-length non-redundant transcripts were assembled in trunk xylem and needle, respectively. Of those, 29 and 20 DEGs were identified involving in the Terpenoid Backbone Biosynthesis pathway (TBB) (KO 00900) and the Diterpenoid Biosynthesis pathway (DB) (KO 00904) in trunk xylem, as well as 35 and 12 DEGs in the TBB and the DB in the needle, respectively. Compared with the low oleoresin-yield germplasm, 13 up-regulated DEGs with complete open reading frame (ORF) sequence were obtained by Venn analysis in trunk xylem and needle of a high oleoresin-yield clone. Subsequently, the expression patterns of the 16 diterpenoid biosynthesis-related genes (including three types of DXS genes) were validated by quantitative real-time PCR in both 12-year-old and 1-year-old high and low oleoresin-yield masson pine clones. The results showed that the 13 genes in the higher yield oleoresin trees were proved to be higher expression in comparison with the lower ones irrespectively of the tree age; as well as MCS (2-C-methyl-d-erythritol 2, 4-cyclodiphosphate synthase), HDS (1-hydroxy-2-methyl-2-(E)-butenyl4-diphosphate synthase), and IDI (Isopentenyl-diphosphate Delta-isomerase) of the methylerythritol 4-phosphate (MEP) pathway demonstrated a higher expression compared to that of other TBB genes. Both cytosol-localized MVA and plastid-localized MEP pathways gave significant contributions to diterpenoid accumulation, and the latter played a more dominant role. To clarify the potential relationships between oleoresin yield and the expression patterns of 16 genes from the TBB and DB pathways, the Pearson correlation analysis was performed using three tissues (young needle, old needle, and trunk xylem) from the 60-year-old tree. Remarkably, HMGR (Hydroxy methylglutaryl-CoA reductase) and ent-CDS (Ent-copalyl diphosphate synthase) were positive correlation with oleoresin yield in the new needle, whereas FPPS (Farnesyl pyrophosphate synthase) was a negative correlation with that in the old needle; furthermore, HMGS (Hydroxy methylglutaryl-CoA synthase), IDI and CDS (Copalyl diphosphate synthase) were positively correlated with oleoresin synthesis in the trunk xylem. Moreover, the Pearson correlation analysis was also carried out to determine the potential relationships between the 13 genes of the TBB pathways and ent-CDS, CDS, and LDS ((13E)-labda-7 13-dien-15-ol synthase) from DB pathways within three tissues of the 60-year-old masson pine. Furthermore, the diterpenoid and diterpene resin acids (DRAs) demonstrated highly positive correlations with oleoresin yield, and DRAs were detected as the most abundant diterpenoids in oleoresin. Promisingly, HMGR, MCS, HDS, IDI, ent-CDS, and CDS may be used as the key target genes of the molecular-assisted selection for high-oleoresin-yielding germplasm in this species.