Metabolome and transcriptome profiling reveal regulatory network and mechanism of flavonoid biosynthesis during color formation of Dioscorea cirrhosa L.

Dioscorea cirrhosa is a plant that is used as a dye as well as in medicine. Many metabolites with pharmacological activity exist in the tubers of D. cirrhosa. However, little is known about the mechanism regulating biosynthesis in these metabolites. In this study, transcriptome and metabolome profiling were performed in four color tubers. A total of 531 metabolites, including 62 flavonoids, were identified. Epicatechin and proanthocyanin B2 were the key metabolites that exhibited high content levels in the four tubers. These metabolites were divided into nine classes with distinct change patterns. A total of 22,865 differentially expressed genes (DEGs) were identified by transcriptome analysis. Among these DEGs, we identified 67 candidate genes related to the flavonoid biosynthesis pathway and three genes that played pivotal roles in proanthocyanin (PA) synthesis. A weighted gene co-expression network analysis (WGCNA) revealed that the two modules, “MEblue” and “MEblack,” were two key gene sets strongly associated with phenylpropanoid and flavonoid biosynthesis. We also found that the plant hormone signal transduction biological process exhibited activity in the late stage of tuber color formation. Additionally, we identified 37 hub transcript factors related to flavonoid biosynthesis, of which 24 were found to be highly associated with flavonoid pathway genes. In addition to the MYB-bHLH-WD40 (MBW) genes, we found that the plant hormone gene families exhibited high expression levels. This study provides a reference for understanding the synthesis of D. cirrhosa tuber metabolites at the molecular level and provides a foundation for the further development of D. cirrhosa related plant pigments as well as its further use in the pharmaceutical industry.