Abstract
Ethylene promotes ripening in fruits as well as the biosynthesis of anthocyanins in plants. However, the question of which ethylene response factors (ERFs) interact with the genes along the anthocyanin biosynthesis pathway is yet to be answered. Herein, we conduct an integrated analysis of transcriptomes and metabolome on fruits of two mulberry genotypes ('Zijin', ZJ, and 'Dashi', DS, with high and low anthocyanin abundance, respectively) at different post-flowering stages. In total, 1035 upregulated genes were identified in ZJ and DS, including MYBA in the MBW complex and anthocyanin related genes such as F3H. A KEGG analysis suggested that flavonoid biosynthesis and plant hormone signaling transduction pathways were significantly enriched in the upregulated gene list. In particular, among 103 ERF genes, the expression of ERF5 showed the most positive correlation with the anthocyanin change pattern across both genotypes and in the post-flowering stages, with a Pearson correlation coefficient (PCC) of 0.93. Electrophoresis mobility shift assay (EMSA) and luciferase assay suggested that ERF5 binds to the promoter regions of MYBA and F3H and transcriptionally activates their gene expression. We elucidated a potential mechanism by which ethylene enhances anthocyanin accumulation in mulberry fruits and highlighted the importance of the ERF5 gene in controlling the anthocyanin content in mulberry species. This knowledge could be used for engineering purposes in future mulberry breeding programs.