Integration of comparative transcriptomics and WGCNA characterizes the regulation of anthocyanin biosynthesis in mung bean (Vigna radiata L.)

Abstract

Mung bean is a dual-use crop widely cultivated in Southeast Asia as a food and medicine resource. The development of new functional mung bean varieties demands identifying new genes regulating anthocyanidin synthesis and investigating their molecular mechanism. In this study, we used high-throughput sequencing technology to generate transcriptome sequence of leaves, petioles, and hypocotyls for investigating the anthocyanins accumulation in common mung bean variety as well as anthocyanidin rich mung bean variety, and to elucidate their molecular mechanisms. 29 kinds of anthocyanin compounds were identified. Most of the anthocyanin components contents were significantly higher in ZL23 compare with AL12. Transcriptome analysis suggested that a total of 93 structural genes encoding the anthocyanin biosynthetic pathway and 273 regulatory genes encoding the ternary complex of MYB-bHLH-WD40 were identified, of which 26 and 78 were differentially expressed in the two varieties. Weighted gene co-expression network analysis revealed that VrMYB3 and VrMYB90 might have enhanced mung bean anthocyanin content by inducing the expression of structural genes such as PAL, 4CL, F3’5’H, LDOX, and F3’H, which was consistent with qRT-PCR results. These findings are envisaged to provide a reference for studying the molecular mechanism of anthocyanin accumulation in mung beans.