Deciphering the Biosynthetic Pathways and Regulatory Networks of the Active Components of Cibotium barometz by Transcriptomic Analysis

Cibotium barometz (L.) J. Sm., a medicinally significant fern in traditional Chinese medicine, is little explored at the genomic level regarding its bioactive compounds. Using an integrated approach combining Illumina and PacBio sequencing technologies, we profiled its root, rachis, and pinna transcriptomes, identifying 12,718, 21,341, and 11,441 unigenes, respectively. Our analysis systematically characterized the transcriptional features of transcription factors (TFs), simple sequence repeats (SSRs), long non-coding RNAs (lncRNAs), and differentially expressed genes (DEGs). Enrichment analyses highlighted the roles of highly expressed unigenes in secondary metabolism. Seventeen key enzymes involved in polysaccharide biosynthesis showed tissue-specific expression patterns. Notably, total polysaccharide content correlated positively with UDP-arabinose 4-epimerase (UXE) expression but negatively with phosphoglucomutase (PGM) and 3,5-epimerase/4-reductase (UER1). Flavonoid accumulation inversely correlated with chalcone synthase (CHS) expression. Two lignin pathways (H-lignin and G-lignin) were characterized, with phenylalanine ammonia-lyase (PAL), cinnamate-4-hydroxylase (C4H), and cinnamyl alcohol dehydrogenase (CAD) as key genes. The absence of ferulate-5-hydroxylase (F5H) explains the undetected S-lignin pathway. Regulatory network analysis revealed positive correlations between PAL expression and NAC72/NAC78/WRKY35 and C4H expression and WRKY65/WRKY69/WRKY71, while a negative correlation was revealed between flavonoid 3′,5′-hydroxylase (F3′5′H) and MYB3R4. This study provides comprehensive transcriptomic insights into C. barometz bioactive compound biosynthesis, serving as a foundation for mechanistic research.