Integrated transcriptomic and metabolomic analyses of anthocyanin biosynthesis in the Leaf of quinoa(Chenopodium quinoa Willd.)leaf
In this study,we analyzed the differentially expressed genes(DEGs)and metabolites(DAMs)between different varieties at different growth stages to investigate the mechanism of anthocyanin synthesis in emerald green and pink quinoa(Chenopodium quinoa Willd.)leaves.11 DEGs were closely related to anthocyanin biosynthesis including 4CL,C3'H,HCT,CHS,CHI,ANR,CYP75B1,UGT79B1,FG3,FG2 and CYP73A.4 DEGs were MYC2,BHLH14,HY5 and TGA transcription factors.The DEGs enriched in seven GO Terms(GO:0009812(Flavonoid metabolic process),GO:0010468(Regulation of gene expression),GO:0051553(Flavone biosynthetic process),GO:0009813(Flavonoid biosynthetic process),GO:0009411(Response to UV),GO:0043169(Cation binding)and GO:0016703(Oxidoreductase activity))that are closely related to anthocyanin biosynthesis.KEGG enrichment analysis showed that there were 6 metabolic pathways significantly related to anthocyanin biosynthesis including Phenylalanine metabolism(ko00360),Phenylpropanoid biosynthesis(ko00940),Flavonoid biosynthesis(ko00941),Plant hormone signal transduction(ko04075),Flavone and flavonol biosynthesis(ko00944),and Circadian rhythm-plant(ko04712).Metabolomics analysis confirmed that Cyanidin 3-O-(3″,6″-O-dimalonyl glucoside)and Naringenin were the key DAMs formed by pink leaves.Transcriptomics-metabolomics joint analysis indicated that flavonoid biosynthesis pathway(ko00941)and anthocyanin biosynthesis pathway(ko00942)were the enrichment pathway of anthocyanin biosynthesis.10 DEGs were verified by Qrt-PCR,and the results showed that the verification results were consistent with the results of transcriptome.