Transcriptomic Analysis Reveals the Flavonoid Biosynthesis Pathway Involved in Rhizome Development in Polygonatum cyrtonema Hua

Polygonatum cyrtonema Hua (P. cyrtonema) rhizomes are rich in flavonoids and other secondary metabolites, exhibiting remarkable antioxidant, anti-tumor, and immunomodulatory effects. Polygonatum flavonoid-biosynthesis-related genes have been characterized already. However, a comprehensive overview of Polygonatum flavonoid biosynthesis pathways is still absent. To articulate the accumulation of the flavonoid biosynthesis pathways, we examined transcriptome changes using Illumina HiSeq from five different tissues and the RNA-seq of 15 samples had over 105 Gb of a clean base, generating a total of 277,955 unigenes. The cDNA libraries of the fruits (F), leaves (L), roots (R), stems (S), and rhizomes (T) of three-year-old P. cyrtonema plants generated 57,591, 53,578, 60,321, 51,530, and 54,935 unigenes. Comparative transcriptome analysis revealed that 379 differentially expressed genes (DEGs) were in the group of F _vs_ T, L _vs_ T, R _vs_ T, and S _vs_ T, and the transcripts of flavonoid-biosynthesis-related DEGs were principally enriched in rhizomes. In addition, combined with WGCNA and the FPKM of five tissues' transcription, nine differentially expressed transcription factor families (MYB, WRKY, AP2/ERF, etc.) were characterized in the red module, the red module positively correlated with rhizome flavonoid accumulation. Quantitative real-time PCR (qRT-PCR) further indicated that BZIP1, C3H31, ERF114, and DREB21 are differentially expressed in rhizomes, accompanied in rhizome development in P. cyrtonema. Therefore, this study provides a foundation for further research into uncovering the accumulation of flavonoid biosynthesis in the rhizomes of P. cyrtonema.