Multi-omics analysis reveals the biosynthesis of flavonoids during the browning process of Malus sieversii explants

Malus sieversii is a precious apple germplasm resource. Browning of explants is one of the most important factors limiting the survival rate of plant tissue culture. In order to explore the molecular mechanism of the browning degree of different strains of Malus sieversii, we compared the dynamic changes of Malus sieversii and Malus robusta Rehd. during the whole browning process using a multi-group method. A total of 44 048 differentially expressed genes (DEGs) were identified by transcriptome analysis on the DNBSEQ-T7 sequencing platform. KEGG enrichment analysis showed that the DEGs were significantly enriched in the flavonoid biosynthesis pathway. In addition, metabonomic analysis showed that (-)-epicatechin, astragalin, chrysin, irigenin, isoquercitrin, naringenin, neobavaisoflavone and prunin exhibited different degrees of free radical scavenging ability in the tissue culture browning process, and their accumulation in different varieties led to differences in the browning degree among varieties. Comprehensive transcriptome and metabonomics analysis of the data related to flavonoid biosynthesis showed that PAL, 4CL, F3H, CYP73A, CHS, CHI, ANS, DFR and PGT1 were the key genes for flavonoid accumulation during browning. In addition, WGCNA analysis revealed a strong correlation between the known flavonoid structure genes and the selected transcriptional genes. Protein interaction predictions demonstrated that 19 transcription factors (7 MYBs and 12 bHLHs) and 8 flavonoid structural genes had targeted relationships. The results show that the interspecific differential expression of flavonoid genes is the key influencing factor of the difference in browning degree between Malus sieversii and Malus robusta Rehd., providing a theoretical basis for further study on the regulation of flavonoid biosynthesis.