Integrated transcriptomic and metabolomic analyses elucidate the mechanism of flavonoid biosynthesis in the regulation of mulberry seed germination under salt stress

Seed propagation is the main method of mulberry expansion in China, an important economic forest species. However, seed germination is the most sensitive stage to various abiotic stresses, especially salinity stress. To reveal the molecular regulatory mechanism of mulberry seed germination under salt stress, flavonoid metabolomics and transcriptomics analyses were performed on mulberry seeds germinated under 50 and 100 mmol/L NaCl stress. Analysis of the flavonoid metabolome revealed that a total of 145 differential flavonoid metabolites (DFMs) were classified into 9 groups, 40 flavonols, 32 flavones, 16 chalcones and 14 flavanones. Among them, 61.4% (89) of the DFMs accumulated continuously with increasing salt concentration, reaching the highest level at a 100 mmol/L salt concentration; these DFMs included quercetin-3-O-glucoside (isoquercitrin), kaempferol (3,5,7,4'-tetrahydroxyflavone), quercetin-7-O-glucoside, taxifolin (dihydroquercetin) and apigenin (4',5,7-trihydroxyflavone), indicating that these flavonoids may be key metabolites involved in the response to salt stress. Transcriptional analysis identified a total of 3055 differentially expressed genes (DEGs), most of which were enriched in flavonoid biosynthesis (ko00941), phenylpropanoid biosynthesis (ko00940) and biosynthesis of secondary metabolites (ko01110). Combined analysis of flavonoid metabolomic and transcriptomic data indicated that phenylalanine ammonia-lyase (PAL), 4-coumarate-CoA ligase (4CL), chalcone synthase (CHS), flavonol synthase (FLS), bifunctional dihydroflavonol 4-reductase/flavanone 4-reductase (DFR) and anthocyanidin reductase (ANR) were the key genes involved in flavonoid accumulation during mulberry seed germination under 50 and 100 mmol/L NaCl stress. In addition, three transcription factors, MYB, bHLH and NAC, were involved in the regulation of flavonoid accumulation under salt stress. The results of quantitative real-time PCR (qRT‒PCR) validation showed that the expression levels of 11 DEGs, including 7 genes involved in flavonoid biosynthesis, under different salt concentrations were consistent with the transcriptomic data, and parallel reaction monitoring (PRM) results showed that the expression levels of 6 key enzymes (proteins) involved in flavonoid synthesis were consistent with the accumulation of flavonoids. This study provides a new perspective for investigating the regulatory role of flavonoid biosynthesis in the regulation of mulberry seed germination under salt stress at different concentrations.