Abstract
Bupleurum falcatum L. is frequently used in traditional herbal medicine in Asia. Saikosaponins (SSs) are the main bioactive ingredients of B. falcatum, but the biosynthetic pathway of SSs is unclear, and the biosynthesis of species-specific phytometabolites is little known. Here we resolved the transcriptome profiles of B. falcatum to identify candidate genes that might be involved in the biosynthesis of SSs. By isoform sequencing (Iso-Seq) analyses of the whole plant, a total of 26.98 Gb of nucleotides were obtained and 124,188 unigenes were identified, and 81,594 unigenes were successfully annotated. A total of 1033 unigenes of 20 families related to the mevalonate (MVA) pathway and methylerythritol phosphate (MEP) pathway of the SS biosynthetic pathway were identified. The WGCNA (weighted gene co-expression network analysis) of these unigenes revealed that only the co-expression module of MEmagenta, which contained 343 unigenes, was highly correlated with the biosynthesis of SSs. Comparing differentially expressed gene analysis and the WGCNA indicated that 130 out of 343 genes of the MEmagenta module exhibited differential expression levels, and genes with the most "hubness" within this module were predicted. Manipulation of these genes might improve the biosynthesis of SSs.