Abstract
Broussonetia papyrifera is rich in flavonoids, which have significant antioxidant, antibacterial, and anti-inflammatory activities and certain pharmacological activities. Nevertheless, scarce transcriptome resources of B. papyrifera have impeded further study regarding the process of its production and accumulation. In this study, RNA-seq was utilized to evaluate the gene expression of B. papyrifera leaves at three distinct developmental phases (T1: young leaves, T3: immature leaves, T4: matured leaves). We obtained 2447 upregulated and 2960 downregulated DEGs, 4657 upregulated and 4804 downregulated DEGs, and 805 upregulated and 484 downregulated DEGs from T1 vs. T3, T1 vs. T4, and T3 vs. T4, respectively. Further research found that the following variables contributed to the formation of flavonoids in the leaves of B. papyrifera: Several important enzyme genes involved in flavonoid production pathways have been discovered. The results demonstrated that the dynamic changing trend of flavonoid contents is related to the expression pattern of the vast majority of essential genes in the biosynthetic pathway. Genes associated in energy and glucose metabolism, polysaccharide, cell wall and cytoskeleton metabolism, signal transduction, and protein and amino acid metabolism may affect the growth and development of B. papyrifera leaves, and eventually their flavonoid content. This study's results offer a strong platform for future research into the metabolic pathways of B. papyrifera.