Recent advances in the transcriptional regulation of the flavonoid biosynthetic pathway

Modulation of anthocyanin accumulation in storage roots of sweetpotato by transcription factor IbMYB1-2 through direct binding to anthocyanin biosynthetic gene promoters

The storage roots of purple-fleshed sweetpotato rich in anthocyanins are considered nutrient-rich foods with health effects. However, the molecular mechanism underlying anthocyanin biosynthesis and regulation remains to be revealed. In this study, IbMYB1-2 was isolated from purple-fleshed sweetpotato "Xuzishu8". The phylogenetic and sequence analysis indicated that IbMYB1-2 belongs to the SG6 subfamily with a conserved bHLH motif. Subcellular localization analysis and transcriptional activity assay revealed that IbMYB1-2 is a key transcriptional activator and is specific to the nucleus. Agrobacterium rhizogenes-mediated overexpression of IbMYB1-2 in sweetpotato through in vivo root transgenic system led to an increase in anthocyanins in the root of sweetpotato. qRT-PCR and transcriptome analysis depicted that the transcript levels of IbMYB1-2, IbbHLH42, and eight structural genes that are associated with the synthesis of anthocyanin were upregulated in overexpressed IbMYB1-2 transgenic roots. Dual-luciferase reporter (DLR) assay and yeast one-hybrid (Y1H) assay demonstrated IbMYB1-2 binding to the promoter regions of IbbHLH42 and other anthocyanin biosynthetic genes, including IbCHS, IbCHI, IbF3H, IbDFR, IbANS, IbGSTF12, IbUGT78D2, and IbUF3GT. Moreover, IbbHLH42 was shown to be an active enhancer for the formation of MYB-bHLH-WD40 (MBW) complex, which strongly supports the promoter activities of the IbCHS, IbANS, IbUGT78D2, and IbGSTF12 genes to induce anthocyanin accumulation. Taken together, our findings not only revealed the underlying regulatory molecular mechanism of IbMYB1-2 for anthocyanin accumulation in the storage roots of sweetpotato but also uncovered a potential mechanism by which IbbHLH42 modulated anthocyanin biosynthesis through a positive feedback regulatory loop.

Transcriptomic-based analysis to identify candidate genes for blue color rose breeding

Analysis of the flower color formation mechanism of 'Rhapsody in Blue' by BF and WF transcriptomes reveals that RhF3'H and RhGT74F2 play a key role in flower color formation. Rosa hybrida has colorful flowers and a high ornamental value. Although rose flowers have a wide range of colors, no blue roses exist in nature, and the reason for this is unclear. In this study, the blue-purple petals (BF) of the rose variety 'Rhapsody in Blue' and the white petals (WF) of its natural mutant were subjected to transcriptome analysis to find genes related to the formation of the blue-purple color. The results showed that the anthocyanin content was significantly higher in BF than in WF. A total of 1077 differentially expressed genes (DEGs) were detected by RNA-Seq analysis, of which 555 were up-regulated and 522 were down-regulated in the WF vs. BF petals. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses of the DEGs revealed that a single gene up-regulated in BF was related to multiple metabolic pathways including metabolic process, cellular process, protein-containing complex, etc. Additionally, the transcript levels of most of the structural genes related to anthocyanin synthesis were significantly higher in BF than in WF. Selected genes were analyzed by qRT-PCR and the results were highly consistent with the RNA-Seq results. The functions of RhF3'H and RhGT74F2 were verified by transient overexpression analyses, and the results confirmed that both affect the accumulation of anthocyanins in 'Rhapsody in Blue'. We have obtained comprehensive transcriptome data for the rose variety 'Rhapsody in Blue'. Our results provide new insights into the mechanisms underlying rose color formation and even blue rose formation.

Comparative transcriptome and genome analysis unravels the response of Tatary buckwheat root to nitrogen deficiency

Tartary buckwheat (Fagopyrum tataricum Garetn.), a dicotyledonous herbaceous crop, has good adaptation to low nitrogen (LN) condition. The plasticity of roots drives the adaption of Tartary buckwheat under LN, but the detailed mechanism behind the response of TB roots to LN remains unclear. In this study, the molecular mechanism of two Tartary buckwheat genotypes' roots with contrasting sensitivity in response to LN was investigated by integrating physiological, transcriptome and whole-genome re-sequencing analysis. LN improved primary and lateral root growth of LN-sensitive genotype, whereas the roots of LN-insensitive genotype showed no response to LN. 2, 661 LN-responsive differentially expressed genes (DEGs) were identified by transcriptome analysis. Of these genes, 17 N transport and assimilation-related and 29 hormone biosynthesis and signaling genes showed response to LN, and they may play important role in Tartary buckwheat root development under LN. The flavonoid biosynthetic genes' expression was improved by LN, and their transcriptional regulations mediated by MYB and bHLH were analyzed. 78 transcription factors, 124 small secreted peptides and 38 receptor-like protein kinases encoding genes involved in LN response. 438 genes were differentially expressed between LN-sensitive and LN-insensitive genotypes by comparing their transcriptome, including 176 LN-responsive DEGs. Furthermore, nine key LN-responsive genes with sequence variation were identified, including FtNRT2.4, FtNPF2.6 and FtMYB1R1. This paper provided useful information on the response and adaptation of Tartary buckwheat root to LN, and the candidate genes for breeding Tartary buckwheat with high N use efficiency were identified.

OsCCRL1 is Essential for Phenylpropanoid Metabolism in Rice Anthers

Phenylpropanoid metabolism and timely tapetal degradation are essential for anther and pollen development, but the underlying mechanisms are unclear. In the current study, to investigate this, we identified and analyzed the male-sterile mutant, osccrl1 (cinnamoyl coA reductase-like 1), which exhibited delayed tapetal programmed cell death (PCD) and defective mature pollen. Map-based cloning, genetic complementation, and gene knockout revealed that OsCCRL1 corresponds to the gene LOC_Os09g32020.2, a member of SDR (short-chain dehydrogenase/reductase) family enzyme. OsCCRL1 was preferentially expressed in the tapetal cells and microspores, and localized to the nucleus and cytoplasm in both rice protoplasts and Nicotiana benthamiana leaves. The osccrl1 mutant exhibited reduced CCRs enzyme activity, less lignin accumulation, delayed tapetum degradation, and disrupted phenylpropanoid metabolism. Furthermore, an R2R3 MYB transcription factor OsMYB103/OsMYB80/OsMS188/BM1, involved in tapetum and pollen development, regulates the expression of OsCCRL1. Finally, the osmyb103 osccrl1 double mutants, exhibited the same phenotype as the osmyb103 single mutant, further indicating that OsMYB103/OsMYB80/OsMS188/BM1 functions upstream of OsCCRL1. These findings help to clarify the role of phenylpropanoid metabolism in male sterility and the regulatory network underlying the tapetum degradation.

Comparative transcriptomic and metabolomic analyses reveal differences in flavonoid biosynthesis between PCNA and PCA persimmon fruit

The fruit of the persimmon (Diospyros kaki.) has high economic and nutritional value and is rich in flavonoids. Flavonoids are essential secondary metabolisms in plants. The association between persimmon astringency and changes in the proanthocyanidins (a flavonoid subclass) content is well-known. However, information on the relationships between different astringency types and other flavonoid subclasses and biosynthetic genes is more limited. In this study, an initial correlation analysis between total flavonoids and fruit astringency type, and KEGG analysis of metabolites showed that flavonoid-related pathways were linked to differences between mature pollination-constant non-astringent (PCNA) varieties ('Jiro' and 'Yohou') and pollination-constant astringent (PCA) fruit varieties ('Zhongshi5' and 'Huojing'). Based on these findings, variations in the expression of genes and metabolites associated with flavonoid biosynthesis were investigated between typical PCNA ('Jiro') and PCA ('Huojing') persimmons during fruit development. The flavonoid concentration in 'Huojing' fruit was significantly higher than that of 'Jiro' fruit, especially, in levels of proanthocyanin precursor epicatechin and anthocyanin cyanidin derivatives. Combined WGCNA and KEGG analyses showed that genes such as PAL, C4H, CHI, CHS, F3H, F3'5'H, FLS, DFR, ANR, ANS, and UF3GT in the phenylpropanoid and flavonoid biosynthesis pathways may be significant factors impacting the proanthocyanin precursor and anthocyanin contents. Moreover, interactions between the R2R3MYB (evm.TU.contig7272.598) and WD40 (evm.TU.contig3208.5) transcription factors were found to be associated with the above structural genes. These findings provide essential information on flavonoid biosynthesis and its regulation in the persimmon and lay a foundation for further investigation into how astringency types affect flavor components in PCNA and PCA persimmons.

Genome-Wide Characterization and Analysis of bHLH Transcription Factors Related to Anthocyanin Biosynthesis in Cinnamomum camphora ('Gantong 1')

Cinnamomum camphora is one of the most commonly used tree species in landscaping. Improving its ornamental traits, particularly bark and leaf colors, is one of the key breeding goals. The basic helix-loop-helix (bHLH) transcription factors (TFs) are crucial in controlling anthocyanin biosynthesis in many plants. However, their role in C. camphora remains largely unknown. In this study, we identified 150 bHLH TFs (CcbHLHs) using natural mutant C. camphora 'Gantong 1', which has unusual bark and leaf colors. Phylogenetic analysis revealed that 150 CcbHLHs were divided into 26 subfamilies which shared similar gene structures and conserved motifs. According to the protein homology analysis, we identified four candidate CcbHLHs that were highly conserved compared to the TT8 protein in A. thaliana. These TFs are potentially involved in anthocyanin biosynthesis in C. camphora. RNA-seq analysis revealed specific expression patterns of CcbHLHs in different tissue types. Furthermore, we verified expression patterns of seven CcbHLHs (CcbHLH001, CcbHLH015, CcbHLH017, CcbHLH022, CcbHLH101, CcbHLH118, and CcbHLH134) in various tissue types at different growth stages using qRT-PCR. This study opens a new avenue for subsequent research on anthocyanin biosynthesis regulated by CcbHLH TFs in C. camphora.

NtMYB12 requires for competition between flavonol and (pro)anthocyanin biosynthesis in Narcissus tazetta tepals

The color of flowers is one of the main characteristics adopted for plants to attract pollinators to ensure the reproductive success of the plant, they are also important in their ornamental appeal in Narcissus plant. In this study, we identified a NtMYB12 locus encoding an R2R3-MYB transcription factor. Comparative transcriptome analysis of loss- and gain- of NtMYB12 tissue relative to wild-type narcissus showed NtMYB12 was mainly involved in flavonol and phenylpropanoid metabolic pathways. Biochemical evidences of dual-luciferase activity and chromatin immunoprecipitation assay supported that MYB12 directly bound to promoters of NtFLS, NtLAR, and NtDFR that were cloned by genome walking assay, and activated NtFLS and NtLAR expression but repressed NtDFR expression. More interestingly, NtMYB12 can interact with NtbHLH1 and NtWD40-1 proteins via R3 domain that were selected by transcriptome-based WGCNA and confirmed by yeast two hybrid, bimolecular fluorescence complementation and coimmunoprecipitation assay. Interaction of NtMYB12 with NtbHLH1 and NtWD40-1 forming MYB-bHLH-WD40 triplex specially activated NtDFR and NtANS expression and promoted (pro)anthocyanin accumulation, while NtMYB12 alone activated NtFLS and NtLAR expression and accumulated flavonols, but repressed NtDFR expression. These results indicated that NtMYB12 alone or NtMYB12-bHLH1-WD40-1 triplex requires for competition of metabolism fluxes between flavonol and (pro)anthocyanin biosynthesis. NtMYB12 dually functions on flavonol and proanthocyanin biogenesis via physically binding to NtFLS and NtLAR promoter activating their expression and on (pro)anthocyanin biosynthesis via NtMYB12-NtWD40-NtbHLH (MBW) triplex activating NtDFR and NtANS expression. Requirement of NtMYB12 alone or MBW complex for the competition between flavonol and anthocyanin biosynthesis results in narcissus colorized petal traits.

Genome-Wide Identification and Characterization of R2R3-MYB Provide Insight into Anthocyanin Biosynthesis Regulation Mechanism of Ananas comosus var. bracteatus

The R2R3-MYB proteins comprise the largest class of MYB transcription factors, which play an essential role in regulating anthocyanin synthesis in various plant species. Ananas comosus var. bracteatus is an important colorful anthocyanins-rich garden plant. The spatio-temporal accumulation of anthocyanins in chimeric leaves, bracts, flowers, and peels makes it an important plant with a long ornamental period and highly improves its commercial value. We conducted a comprehensive bioinformatic analysis of the R2R3-MYB gene family based on genome data from A. comosus var. bracteatus. Phylogenetic analysis, gene structure and motif analysis, gene duplication, collinearity, and promoter analysis were used to analyze the characteristics of this gene family. In this work, a total of 99 R2R3-MYB genes were identified and classified into 33 subfamilies according to phylogenetic analysis, and most of them were localized in the nucleus. We found these genes were mapped to 25 chromosomes. Gene structure and protein motifs were conserved among AbR2R3-MYB genes, especially within the same subfamily. Collinearity analysis revealed four pairs of tandem duplicated genes and 32 segmental duplicates in AbR2R3-MYB genes, indicating that segmental duplication contributed to the amplification of the AbR2R3-MYB gene family. A total of 273 ABRE responsiveness, 66 TCA elements, 97 CGTCA motifs, and TGACG motifs were the main cis elements in the promoter region under response to ABA, SA, and MEJA. These results revealed the potential function of AbR2R3-MYB genes in response to hormone stress. Ten R2R3-MYBs were found to have high homology to MYB proteins reported to be involved in anthocyanin biosynthesis from other plants. RT-qPCR results revealed the 10 AbR2R3-MYB genes showed tissue-specific expression patterns, six of them expressed the highest in the flower, two genes in the bract, and two genes in the leaf. These results suggested that these genes may be the candidates that regulate anthocyanin biosynthesis of A. comosus var. bracteatus in the flower, leaf, and bract, respectively. In addition, the expressions of these 10 AbR2R3-MYB genes were differentially induced by ABA, MEJA, and SA, implying that these genes may play crucial roles in hormone-induced anthocyanin biosynthesis. Our study provided a comprehensive and systematic analysis of AbR2R3-MYB genes and identified the AbR2R3-MYB genes regulating the spatial-temporal anthocyanin biosynthesis in A. comosus var. bracteatus, which would be valuable for further study on the anthocyanin regulation mechanism of A. comosus var. bracteatus.

Foxtail millet MYB-like transcription factor SiMYB16 confers salt tolerance in transgenic rice by regulating phenylpropane pathway

R2R3-MYB transcription factors play an important role in the synthesis of phenylpropanoid-derived compounds, which in turn provide salt tolerance in plant. In this study, we found that the expression of foxtail millet R2R3-MYB factor SiMYB16 can be induced by salt and drought. SiMYB16 is localized in the nucleus and acts as a transcriptional activator. Phylogenetic analysis indicates that SiMYB16 belongs to the R2R3-MYB transcription factor family subgroup 24. Transgenic rice expressing SiMYB16 (OX16) had a higher survival rate, lower malondialdehyde content, and heavier fresh weight compared with type (WT) under salt stress conditions. The transgenic plants also had a higher germination rate in salt treatment conditions and higher yield in the field compared with wild-type plants. Transcriptome analysis revealed that the up-regulated differential expression genes in the transgenic rice were mainly involved in phenylpropanoid biosynthesis, fatty acid elongation, phenylalanine metabolism, and flavonoid biosynthesis pathways. Quantitative real-time PCR analysis also showed that the genes encoding the major enzymes in the lignin and suberin biosynthesis pathways had higher expression level in SiMYB16 transgenic plants. Correspondingly, the content of flavonoid and lignin, and the activity of fatty acid synthase increased in SiMYB16 transgenic rice compared with wild-type plants under salt stress treatment. These results indicate that SiMYB16 gene can enhance plant salt tolerance by regulating the biosynthesis of lignin and suberin.

Light-Induced Flavonoid Biosynthesis in Sinopodophyllum hexandrum with High-Altitude Adaptation

Sinopodophyllum hexandrum is a perennial alpine herb producing the anti-cancer metabolite podophyllotoxin (PPT). Although the adaptation of S. hexandrum to high altitudes has been demonstrated and the effects of temperature, precipitation, and UV-B light on plant growth and metabolite accumulation have been studied, knowledge on the role of flavonoid biosynthesis in adapting to high altitudes is limited. In this study, light intensity, amount and type of flavonoids, and differentially expressed proteins (DEPs) and genes (DEGs) at 2300 and 3300 m were analyzed by HPLC, proteomic, transcriptomic, and qRT-PCR analysis. We found that higher light intensity correlated with greater flavonoid, flavonol, and anthocyanin content as well as higher anthocyanin to total flavonoid and flavonol ratios observed at the higher altitude. Based on proteomic and transcriptomic analyses, nine DEPs and 41 DEGs were identified to be involved in flavonoid biosynthesis and light response at 3300 m. The relative expression of nine genes (PAL, CHS1, IFRL, ANS, MYB4, BHLH137, CYP6, PPO1, and ABCB19) involved in flavonoid biosynthesis and seven genes (HSP18.1, HSP70, UBC4, ERF5, ERF9, APX3, and EX2) involved in light stress were observed to be up-regulated at 3300 m compared with 2300 m. These findings indicate that light intensity may play a regulatory role in enhancing flavonoid accumulation that allows S. hexandrum to adapt to elevated-altitude coupled with high light intensity.

Physiological and transcriptomic study reveal SeNPs-mediated AsIII stress detoxification mechanisms involved modulation of antioxidants, metal transporters, and transcription factors in Glycine max L. (Merr.) roots

Physiological changes and genome-wide alteration in gene expression were performed in soybean (Glycine max [L.] Merr.) roots exposed to AsⅢ (25 μmol/L) alone and supplemented with selenium nanoparticles (SeNPs) at the concentration of 10 and 25 μmol/L at the V2 growth stage. Excessive arsenic in the root zone poses a potential threat to soybean yield, particularly to roots, due to the limited translocation of AsIII from root to shoot in the case of soybean. We hypothesized that SeNPs can relieve AsⅢ toxicity to soybean root by reducing the AsⅢ uptake and regulating the internal tolerance mechanism of the plants. Results accomplished that SeNPs had positive impact on soybean dry weight and roots parameters under AsⅢ stress. Then, we further evaluated physiological indexes, whole genome transcriptomic analysis and quantitative real-time PCR to elucidate the underlying mechanism of AsⅢ tolerance under SeNPs supplementation. Under the condition of AsⅢ-stress, SeNPs exposure significantly reduced the electrolyte leakage, O₂^-•, H₂O₂ and MDA accumulation while increasing the antioxidants level. The RNA-seq dataset revealed total of 5819 up and 7231 down expressed DEGs across all libraries. The number of exclusively regulated genes were higher under As + SeNP10 (4909) treatment than in the AsⅢ-alone (4830) and As + SeNP25 (3311) treatments. The KEGG and GO analyses revealed that stress responsive DEGs such as glutathione S-transferase, glutathione peroxidase, ascorbate, glutaredoxin, thioredoxin, and phytochelatins synthase are responsible for AsⅢ tolerance under the SeNPs supplementation. Similarly, sulfate transporter, and ABC transporters (ATP-binding cassettes) expression were induced, and aquaporin channels related DEGs expression were reduced under SeNPs application in AsⅢ exposure condition. Furthermore, the expression of molecular chaperones (HSP) and transcription factors (MYB, bZIP, bHLH, and HSFs) were increased in SeNPs treatment groups. These results provide vital information of AsⅢ tolerance mechanism in response to SeNPs in soybean. We suggest that functional characterization of these genes will help us learn more about the SeNPs responsive arsenic tolerance mechanism in soybean.

Integrative Analysis of Metabolome and Transcriptome Reveals the Mechanism of Color Formation in Yellow-Fleshed Kiwifruit

During the development of yellow-fleshed kiwifruit (Actinidia chinensis), the flesh appeared light pink at the initial stage, the pink faded at the fastest growth stage, and gradually changed into green. At the maturity stage, it showed bright yellow. In order to analyze the mechanism of flesh color change at the metabolic and gene transcription level, the relationship between color and changes of metabolites and key enzyme genes was studied. In this study, five time points (20 d, 58 d, 97 d, 136 d, and 175 d) of yellow-fleshed kiwifruit were used for flavonoid metabolites detection and transcriptome, and four time points (20 d, 97 d, 136 d, and 175 d) were used for targeted detection of carotenoids. Through the analysis of the content changes of flavonoid metabolites, it was found that the accumulation of pelargonidin and cyanidin and their respective anthocyanin derivatives was related to the pink flesh of young fruit, but not to delphinidin and its derivative anthocyanins. A total of 140 flavonoid compounds were detected in the flesh, among which anthocyanin and 76% of the flavonoid compounds had the highest content at 20 d, and began to decrease significantly at 58 d until 175 d, resulting in the pale-pink fading of the flesh. At the mature stage of fruit development (175 d), the degradation of chlorophyll and the increase of carotenoids jointly led to the change of flesh color from green to yellow, in addition to chlorophyll degradation. In kiwifruit flesh, 10 carotenoids were detected, with none of them being linear carotenoids. During the whole development process of kiwifruit, the content of β-carotene was always higher than that of α-carotene. In addition, β-cryptoxanthin was the most-accumulated pigment in the kiwifruit at 175 d. Through transcriptome analysis of kiwifruit flesh, seven key transcription factors for flavonoid biosynthesis and ten key transcription factors for carotenoid synthesis were screened. This study was the first to analyze the effect of flavonoid accumulation on the pink color of yellow-fleshed kiwifruit. The high proportion of β-cryptoxanthin in yellow-fleshed kiwifruit was preliminarily found. This provides information on metabolite accumulation for further revealing the pink color of yellow-fleshed kiwifruit, and also provides a new direction for the study of carotenoid biosynthesis and regulation in yellow-fleshed kiwifruit.

Basic Helix-Loop-Helix Transcription Factors: Regulators for Plant Growth Development and Abiotic Stress Responses

Plant basic helix-loop-helix (bHLH) transcription factors are involved in many physiological processes, and they play important roles in the abiotic stress responses. The literature related to genome sequences has increased, with genome-wide studies on the bHLH transcription factors in plants. Researchers have detailed the functionally characterized bHLH transcription factors from different aspects in the model plant Arabidopsis thaliana, such as iron homeostasis and abiotic stresses; however, other important economic crops, such as rice, have not been summarized and highlighted. The bHLH members in the same subfamily have similar functions; therefore, unraveling their regulatory mechanisms will help us to identify and understand the roles of some of the unknown bHLH transcription factors in the same subfamily. In this review, we summarize the available knowledge on functionally characterized bHLH transcription factors according to four categories: plant growth and development; metabolism synthesis; plant signaling, and abiotic stress responses. We also highlight the roles of the bHLH transcription factors in some economic crops, especially in rice, and discuss future research directions for possible genetic applications in crop breeding.

HDAC19 recruits ERF4 to the MYB5a promoter and diminishes anthocyanin accumulation during grape ripening

DNA acetylation alters the expression of responsive genes during plant development. In grapes (Vitis vinifera), however, little is known about this regulatory mechanism. In the present study, 'Kyoho' grapes treated with trichostatin A (TSA, a deacetylase inhibitor) were used for transcriptome sequencing and quantitative proteomics analysis. We observed that acetylation was associated with anthocyanin accumulation and gene expression. Acetylation positively regulated phenylalanine metabolism and flavonoid biosynthesis pathways. Using omics analysis, we detected an increase in the levels of the AP2/EREBP transcription factor family after TSA treatment, indicating its association with acetylation-deacetylation dynamics in grapes. Furthermore, ethylene response factor 4 (ERF4) physically interacted with VvHDAC19, a histone deacetylase, which synergistically reduced the expression of target genes involved in anthocyanin biosynthesis owing to the binding of VvERF4 to the GCC-box cis-regulatory element in the VvMYB5a promoter. VvHDAC19 and VvERF4 also controlled anthocyanin biosynthesis and accumulation by regulating acetylation levels of histones H3 and H4. Therefore, alterations in histone modification can significantly regulate the expression of genes involved in anthocyanin biosynthesis and affect grape ripening.

A transcriptional complex of FtMYB102 and FtbHLH4 coordinately regulates the accumulation of rutin in Fagopyrum tataricum

Tartary buckwheat is rich in flavonoids, which not only play an important role in the plant-environment interaction, but are also beneficial to human health. Rutin is a therapeutic flavonol which is massively accumulated in Tartary buckwheat. It has been demonstrated that transcription factors control rutin biosynthesis. However, the transcriptional regulatory network of rutin is not fully clear. In this study, through transcriptome and target metabolomics, we validated the role of FtMYB102 and FtbHLH4 TFs at the different developmental stages of Tartary buckwheat. The elevated accumulation of rutin in the sprout appears to be closely associated with the expression of FtMYB102 and FtbHLH4. Yeast two-hybrid, transient luciferase activity and co-immunoprecipitation demonstrated that FtMYB102 and FtbHLH4 can interact and form a transcriptional complex. Moreover, yeast one-hybrid showed that both FtMYB102 and FtbHLH4 directly bind to the promoter of chalcone isomerase (CHI), and they can coordinately induce CHI expression as shown by transient luciferase activity assay. Finally, we transferred FtMYB102 and FtbHLH4 into the hairy roots of Tartary buckwheat and found that they both can promote the accumulation of rutin. Our results indicate that FtMYB102 and FtbHLH4 can form a transcriptional complex by inducing CHI expression to coordinately promote the accumulation of rutin.

Metabolomics and transcriptomics provide insights into the flavonoid biosynthesis pathway in the roots of developing Aster tataricus

Aster tataricus (L.) is an important medicinal plant in China. Its roots are rich in flavonoids, the main medicinal components. However, the molecular basis of flavonoid biosynthesis in the roots of A. tataricus remains unclear. In this study, the content of total flavonoid of A. tataricus roots at different developmental stages was measured first, and the results showed that the content of total flavonoid gradually decreased from September to November, which may be caused by the stagnation of A. tataricus growth due to the decrease in temperature after September. Then, an integrated analysis of transcriptome and metabolome was conducted on five developing stages of A. tataricus roots to identify flavonoid compositions and potential genes involved in flavonoid biosynthesis. A total of 80 flavonoid metabolites, of which 75% were flavonols and flavonoids, were identified in metabolomic analyses, among which isorhamnetin, kaempferol, quercetin, and myricetin were the main skeletons of these flavonoids. Cluster analysis divided these 80 flavonoids into 3 clusters. The compounds in cluster I mainly accumulated in S1, S3, and S5. In cluster II, the relative content of the flavonoid metabolites showed an upward trend from S2 to S4. In cluster III, the flavonoids decreased from S1 to S5. A total of 129 structural genes, including 43 PAL, 23 4CL, 9 C4H, 4 CHS, 18 CHI, 3 F3H, 5 F3'H, 1 F3'5'H, 21 FLS, and 2 FSII, and 65 transcription factors, including 22 AP2/ERF, 7 bHLH, 5 bZIP, 8 MYB, 11 NAC, and 12 WRKY, showed significant correlation with total flavonoid content. Eighteen genes (7 4CL, 5 C4H, 2 CHI, 1 F3H, and 3 FLS) and 30 genes (5 PAL, 9 4CL, 1 C4H, 2 CHI, 1 F3H, 1 DFR, 7 3AT, 1 BZ1, and 3 UGT79B1) were identified as key structural genes for kaempferol and anthocyanins biosynthesis, respectively. Our study provides valuable information for understanding the mechanism of flavonoid biosynthesis in A. tataricus root.

Integrative analysis of metabolome and transcriptome reveals a dynamic regulatory network of potato tuber pigmentation

Potatoes consist of flavonoids that provide health benefits for human consumers. To learn more about how potato tuber flavonoid accumulation and flesh pigmentation are controlled, we analyzed the transcriptomic and metabolomic profile of potato tubers from three colored potato clones at three developmental phases using an integrated approach. From the 72 flavonoids identified in pigmented flesh, differential abundance was noted for anthocyanins, flavonols, and flavones. Weighted gene co-expression network analysis further allowed modules and candidate genes that positively or negatively regulate flavonoid biosynthesis to be identified. Furthermore, an R2R3-MYB repressor StMYB3 and an R3-MYB repressor StMYBATV involved in the modulation of anthocyanin biosynthesis during tuber development were identified. Both StMYB3 and StMYBATV could interact with the cofactor StbHLH1 and repress anthocyanin biosynthesis. Our results indicate a feedback regulatory mechanism of a coordinated MYB activator-repressor network on fine-tuning of potato tuber pigmentation during tuber development.

Genome-Wide Characterization and Analysis of R2R3-MYB Genes Related to Fruit Ripening and Stress Response in Banana (Musa acuminata L. AAA Group, cv. 'Cavendish')

MYB is an important type of transcription factor in eukaryotes. It is widely involved in a variety of biological processes and plays a role in plant morphogenesis, growth and development, primary and secondary metabolite synthesis, and other life processes. In this study, bioinformatics methods were used to identify the R2R3-MYB transcription factor family members in the whole Musa acuminata (DH-Pahang) genome, one of the wild ancestors of banana. A total of 280 MaMYBs were obtained, and phylogenetic analysis indicated that these MaMYBs could be classified into 33 clades with MYBs from Arabidopsis thaliana. The amino acid sequences of the R2 and R3 Myb-DNA binding in all MaMYB protein sequences were quite conserved, especially Arg-12, Arg-13, Leu-23, and Leu-79. Distribution mapping results showed that 277 MaMYBs were localized on the 11 chromosomes in the Musa acuminata genome. The MaMYBs were distributed unevenly across the 11 chromosomes. More than 40.0% of the MaMYBs were located in collinear fragments, and segmental duplications likely played a key role in the expansion of the MaMYBs. Moreover, the expression profiles of MaMYBs in different fruit development and ripening stages and under various abiotic and biotic stresses were investigated using available RNA-sequencing data to obtain fruit development, ripening-specific, and stress-responsive candidate genes. Weighted gene co-expression network analysis (WGCNA) was used to analyze transcriptome data of banana from the above 11 samples. We found MaMYBs participating in important metabolic biosynthesis pathways in banana. Collectively, our results represent a comprehensive genome-wide study of the MaMYB gene family, which should be helpful in further detailed studies on MaMYBs functions related to fruit development, postharvest ripening, and the seedling response to stress in an important banana cultivar.

Transcriptome and metabolome analyses reveal anthocyanins pathways associated with fruit color changes in plum (Prunus salicina Lindl.)

Plum (Prunus salicina Lindl.) is one of the most widely cultivated and important fruit trees in temperate and cold regions. Fruit color is a significant trait relating to fruit quality in plum. However, its development mechanism has not been studied from the aspects of transcriptional regulation and metabolomic progress. To reveal the mechanism of fruit color developments in plums, we selected the fruits of two plum cultivars, 'Changli84' (Ch84, red fruit) and 'Dahuangganhe' (D, yellow fruit) as plant materials for transcriptome sequencing and metabolomic analysis were performed. Based on the data of transcriptome and metabolome at three fruit developmental stages, young fruit stage, color-change stage, and maturation stage, we identified 2,492 differentially expressed genes (DEGs) and 54 differential metabolites (DMs). The KEGG analysis indicated that "Flavonoid biosynthesis" was significantly enriched during three fruit development stages. Some DEGs in the "Flavonoid biosynthesis" pathway, had opposite trends between Ch84 and D, including chalcone synthase (CHS), dihydroflavonol 4-reductase (DFR) and flavonol synthase (FLS). Also, the genes encoding MYB-bHLH-WD (MBW) protein complexes, especially MYBs and bHLHs, showed a close relationship with plum fruit color. In the current study, DMs like procyanidin B1, cyanidin 3-glucoside, and cyanidin-3-O-alpha-arabinopyranoside were key pigments (or precursors), while the carotene and carotenoids did not show key relationships with fruit color. In conclusion, the anthocyanins dominate the color change of plum fruit. Carotenes and carotenoids might be related to the color of plum fruit, but do not play a dominate role.

Classification and Prediction by Pigment Content in Lettuce (Lactuca sativa L.) Varieties Using Machine Learning and ATR-FTIR Spectroscopy

Green or purple lettuce varieties produce many secondary metabolites, such as chlorophylls, carotenoids, anthocyanins, flavonoids, and phenolic compounds, which is an emergent search in the field of biomolecule research. The main objective of this study was to use multivariate and machine learning algorithms on Attenuated Total Reflectance Fourier Transform Infrared Spectroscopy (ATR-FTIR)-based spectra to classify, predict, and categorize chemometric attributes. The cluster heatmap showed the highest efficiency in grouping similar lettuce varieties based on pigment profiles. The relationship among pigments was more significant than the absolute contents. Other results allow classification based on ATR-FTIR fingerprints of inflections associated with structural and chemical components present in lettuce, obtaining high accuracy and precision (>97%) by using principal component analysis and discriminant analysis (PCA-LDA)-associated linear LDA and SVM machine learning algorithms. In addition, PLSR models were capable of predicting Chla, Chlb, Chla+b, Car, AnC, Flv, and Phe contents, with R2P and RPDP values considered very good (0.81−0.88) for Car, Anc, and Flv and excellent (0.91−0.93) for Phe. According to the RPDP metric, the models were considered excellent (>2.10) for all variables estimated. Thus, this research shows the potential of machine learning solutions for ATR-FTIR spectroscopy analysis to classify, estimate, and characterize the biomolecules associated with secondary metabolites in lettuce.

CmHY5 functions in apigenin biosynthesis by regulating flavone synthase II expression in chrysanthemum flowers

The predominant flavones in the ray florets of chrysanthemum flowers are apigenin and its derivatives. CmHY5 participates in apigenin biosynthesis by directly regulating the expression of FNSII-1 in chrysanthemum. Chrysanthemum (Chrysanthemum morifolium) flowers have been used for centuries as functional food and in herbal tea and traditional medicine. The chrysanthemum flower contains significant amounts of the biologically active compound flavones, which has medicinal properties. However, the mechanism regulating flavones biosynthesis in chrysanthemum flowers organs is still unclear. Here, we compared the transcriptomes and metabolomes of different floral organs between two cultivars with contrasting flavone levels in their flowers. We identified 186 flavonoids by metabolome analysis. The predominant flavones in the ray florets of chrysanthemum flowers are apigenin and its derivatives, of which the contents are highly correlated with the expression of flavones synthase II gene CmFNSII-1. We also determined that CmHY5 is a direct upstream regulator of CmFNSII-1 transcription. We showed that CmHY5 RNAi interference lines in chrysanthemum have lower contents of apigenin compared to wild-type chrysanthemum. Our results demonstrated that CmHY5 participates in flavone biosynthesis by directly regulating the expression of FNSII-1 in chrysanthemum.

Anthocyanins Profiling Analysis and RNA-Seq Revealed the Dominating Pigments and Coloring Mechanism in Cyclamen Flowers

Pigments in cyclamen (Cyclamen purpurascens) endows flowers with great ornamental and medicinal values. However, little is known about the biosynthetic pathways of pigments, especially anthocyanins, in cyclamen flowers. Herein, anthocyanins profiling and RNA-Seq were used to decipher the molecular events using cyclamen genotypes of red (HXK) or white (BXK) flowers. We found that red cyclamen petals are rich in cyanidin-3-O-rutinoside, cyanidin-3-O-glucoside, delphinidin-3-O-glucoside, malvidin-3-O-glucoside, peonidin-3-O-rutinoside, quercetin-3-O-glucoside, and ruti. In addition, our transcriptomics data revealed 3589 up-regulated genes and 2788 down-regulated genes comparing the BXK to HXK. Our rich dataset also identified eight putative key genes for anthocyanin synthesis, including four chalcone synthase (CHS, g13809_i0, g12097_i0, g18851_i0, g36714_i0), one chalcone isomerase (CHI, g26337_i0), two flavonoid 3-hydroxylase (F3'H, g14710_i0 and g15005_i0), and one anthocyanidin synthase (ANS, g18981_i0). Importantly, we found a 2.5 order of magnitude higher expression of anthocyanin 3-O-glucosyltransferase (g8206_i0), which encodes a key gene in glycosylation of anthocyanins, in HXK compared to BXK. Taken together, our multiomics approach demonstrated massive changes in gene regulatory networks and anthocyanin metabolism in controlling cyclamen flower color.

Transcriptome and Metabolome Analysis Revealed That Exogenous Spermidine-Modulated Flavone Enhances the Heat Tolerance of Lettuce

Lettuce is sensitive to high temperature, and exogenous spermidine can improve heat tolerance in lettuce, but its intrinsic mechanism is still unclear. We analyzed the effects of exogenous spermidine on the leaf physiological metabolism, transcriptome and metabolome of lettuce seedlings under high-temperature stress using the heat-sensitive lettuce variety 'Beisansheng No. 3' as the material. The results showed that exogenous spermidine increased the total fresh weight, total dry weight, root length, chlorophyll content and total flavonoid content, increased the activities of antioxidant enzymes such as superoxide dismutase (SOD), peroxidase (POD) and catalase (CAT), and decreased malondialdehyde (MDA) content in lettuce under high temperature stress. Transcriptome and metabolome analyses revealed 818 differentially expressed genes (DEGs) and 393 metabolites between water spray and spermidine spray treatments under high temperature stress, and 75 genes from 13 transcription factors (TF) families were included in the DEGs. The Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analysis of DEG contains pathways for plant-pathogen interactions, photosynthesis-antennal proteins, mitogen-activated protein kinase (MAPK) signaling pathway and flavonoid biosynthesis. A total of 19 genes related to flavonoid synthesis were detected. Most of these 19 DEGs were down-regulated under high temperature stress and up-regulated after spermidine application, which may be responsible for the increase in total flavonoid content. We provide a possible source and conjecture for exploring the mechanism of exogenous spermidine-mediated heat tolerance in lettuce.

A naturally-occurring phenomenon of flower color change during flower development in Xanthoceras sorbifolium

Xanthoceras sorbifolium (yellowhorn) is originated in China and is a unique tree in northern China. Yellowhorn is very popular because of the gradual color change of its flower at different flower developmental stages. During flowering development, the color at the base of yellowhorn flower petals gradually changes from yellow to purple. The mechanism of this miraculous phenomenon is still unclear. Here we show that anthocyanin accumulation during flowering development is the main reason for this color change. RT-PCR results show that the expression level of a variety of anthocyanin biosynthesis genes changes in different flower developmental stages. Realtime results show that the expression changes of these anthocyanin biosynthesis genes are positively regulated by a cluster of R2R3-MYB transcription factor genes, XsMYB113s. Furthermore, the DNA methylation analysis showed that CHH methylation status on the transposon element near the XsMYB113-1 influence its transcript level during flowering development. Our results suggest that dynamic epigenetic regulation of the XsMYB113-1 leads to the accumulation of anthocyanins during yellowhorn flower color change. These findings reemphasize the role of epigenetic regulation in flower development and provide a foundation for further studies of epigenetic regulation in long-lived woody perennials.

DhMYB2 and DhbHLH1 regulates anthocyanin accumulation via activation of late biosynthesis genes in Phalaenopsis-type Dendrobium

Phalaenopsis-type Dendrobium is a popular orchid with good ornamental and market value. Despite their popularity, molecular regulation of anthocyanin biosynthesis during flower development remains poorly understood. In this study, we systematically investigated the regulatory roles of the transcription factors DhMYB2 and DhbHLH1 in anthocyanins biosynthesis. Gene expression analyses indicated that both DhMYB2 and DhbHLH1 are specifically expressed in flowers and have similar expression patterns, showing high expression in purple floral tissues with anthocyanin accumulation. Transcriptomic analyses showed 29 differentially expressed genes corresponding to eight enzymes in anthocyanin biosynthesis pathway have similar expression patterns to DhMYB2 and DhbHLH1, with higher expression in the purple lips than the yellow petals and sepals of Dendrobium 'Suriya Gold'. Further gene expression analyses and Pearson correlation matrix analyses of Dendrobium hybrid progenies revealed expression profiles of DhMYB2 and DhbHLH1 were positively correlated with the structural genes DhF3'H1, DhF3'5'H2, DhDFR, DhANS, and DhGT4. Yeast one-hybrid and dual-luciferase reporter assays revealed DhMYB2 and DhbHLH1 can bind to promoter regions of DhF3'H1, DhF3'5'H2, DhDFR, DhANS and DhGT4, suggesting a role as transcriptional activators. These results provide new evidence of the molecular mechanisms of DhMYB2 and DhbHLH1 in anthocyanin biosynthesis in Phalaenopsis-type Dendrobium.

A chromosome-level genome assembly for Dracaena cochinchinensis reveals the molecular basis of its longevity and formation of dragon's blood

Dracaena, a remarkably long-lived and slowly maturing species of plant, is world famous for its ability to produce dragon's blood, a precious traditional medicine used by different cultures since ancient times. However, there is no detailed and high-quality genome available for this species at present; thus, the molecular mechanisms that underlie its important traits are largely unknown. These factors seriously limit the protection and regeneration of this rare and endangered plant resource. Here, we sequenced and assembled the genome of Dracaena cochinchinensis at the chromosome level. The D. cochinchinensis genome covers 1.21 Gb with a scaffold N50 of 50.06 Mb and encodes 31 619 predicted protein-coding genes. Analysis showed that D. cochinchinensis has undergone two whole-genome duplications and two bursts of long terminal repeat insertions. The expansion of two gene classes, cis-zeatin O-glucosyltransferase and small auxin upregulated RNA, were found to account for its longevity and slow growth. Two transcription factors (bHLH and MYB) were found to be core regulators of the flavonoid biosynthesis pathway, and reactive oxygen species were identified as the specific signaling molecules responsible for the injury-induced formation of dragon's blood. Our study provides high-quality genomic information relating to D. cochinchinensis and significant insight into the molecular mechanisms responsible for its longevity and formation of dragon's blood. These findings will facilitate resource protection and sustainable utilization of Dracaena.

Integrated metabolomic and transcriptomic analysis of the anthocyanin and proanthocyanidin regulatory networks in red walnut natural hybrid progeny leaves

Background

Walnuts are among the most important dry fruit crops worldwide, typically exhibiting green leaves and yellow-brown or gray-yellow seed coats. A specific walnut accession with red leaves and seed coats, 'RW-1', was selected for study because of its high anthocyanin and proanthocyanidin (PA) contents. Anthocyanins and PAs are important secondary metabolites and play key roles in plant responses to biotic and abiotic stresses. However, few studies have focused on the molecular mechanism of anthocyanin biosynthesis in walnuts.

Methods

In this study, we determined the anthocyanin and PA components and their contents in different color leaves of 'RW-1' natural hybrid progenies at various developmental stages. Integrated transcriptome and metabolome analyses were used to identify the differentially expressed genes (DEGs) and differentially accumulated metabolites (DAMs). We also performed conjoint analyses on DEGs and DAMs to ascertain the degree pathways, and explore the regulation of anthocyanin and PA biosynthesis.

Results

The results of widely targeted metabolome profiling and anthocyanin detection revealed 395 substances, including four PAs and 26 anthocyanins, in red (SR) and green leaves (SG) of 'RW-1' natural hybrid progenies. From the research, the contents of all anthocyanin components in SR were higher than that in SG. Among them, the contents of delphinidin 3-O-galactoside, cyanidin 3-O-galactoside, delphinidin 3-O-arabinoside and cyanidin 3-O-glucoside were significantly higher than others, and they were considered as the main types of anthocyanins. However, nine anthocyanins were detected only in SR. For PAs, the content of procyanidin C1 was higher in SR compared with SG, while procyanidin B1 and procyanidin B3 were higher in SR-1 and SR-3 but downregulated in SR-2 compared with the controls. Furthermore, transcriptome analysis revealed that the expressions of structural genes (C4H, F3H, F3'5'H, UFGT, LAR and ANR), three MYBs predicted as the activators of anthocyanin and PA biosynthesis, two MYBs predicted as the repressors of anthocyanin biosynthesis, and five WD40s in the anthocyanin and PA biosynthetic pathways were significantly higher in the SR walnuts. Gene-metabolite correlation analyses revealed a core set of 31 genes that were strongly correlated with four anthocyanins and one PA metabolites. The alteration of gene coding sequence altered the binding or regulation of regulatory factors to structural genes in different color leaves, resulting in the effective increase of anthocyanins and PAs accumulation in red walnut.

Conclusions

This study provides valuable information on anthocyanin and PA metabolites and candidate genes for anthocyanin and PA biosynthesis, yielding new insights into anthocyanin and PA biosynthesis in walnuts.

Selenium Application Enhances the Accumulation of Flavones and Anthocyanins in Bread Wheat (Triticum aestivum L.) Grains

Selenium (Se) biofortification in wheat reduces the risk of Se deficiency in humans. Se biofortification increases the concentration of Se and anthocyanins in wheat grains. However, it is unknown whether Se biofortification can enhance flavonoids other than anthocyanins and the mechanism underlying flavonoid accumulation in wheat grains. Here, foliar application of selenite solution in wheat was conducted 10 days after flowering. Metabolite profiling and transcriptome sequencing were performed in Se-treated grains. A significant increase in the total contents of Se, anthocyanins, and flavonoids was observed in Se-treated mature grains. Twenty-seven significantly increased flavonoids were identified in Se-treated immature grains. The significant accumulation of flavones (tricin, tricin derivatives, and chrysoeriol derivatives) was detected, and six anthocyanins, dihydroquercetin (the precursor for anthocyanin biosynthesis) and catechins were also increased. Integrated analysis of metabolites and transcriptome revealed that Se application enhanced the biosynthesis of flavones, dihydroquercetin, anthocyanins, and catechins by increasing the expression levels of seven key structural genes in flavonoid biosynthesis (two TaF3Hs, two TaDFRs, one TaF3'5'H, one TaOMT, and one TaANR). Our findings shed new light on the molecular mechanism underlying the enhancement in flavonoid accumulation by Se supplementation and pave the way for further enhancing the nutritional value of wheat grains.

Tropaeolum majus R2R3 MYB Transcription Factor TmPAP2 Functions as a Positive Regulator of Anthocyanin Biosynthesis

Anthocyanins are an important group of water-soluble and non-toxic natural pigments with antioxidant and anti-inflammatory properties that can be found in flowers, vegetables, and fruits. Anthocyanin biosynthesis is regulated by several different types of transcription factors, including the WD40-repeat protein Transparent Testa Glabra 1 (TTG1), the bHLH transcription factor Transparent Testa 8 (TT8), Glabra3 (GL3), Enhancer of GL3 (EGL3), and the R2R3 MYB transcription factor Production of Anthocyanin Pigment 1 (PAP1), PAP2, MYB113, and MYB114, which are able to form MYB-bHLH-WD40 (MBW) complexes to regulate the expression of late biosynthesis genes (LBGs) in the anthocyanin biosynthesis pathway. Nasturtium (Tropaeolum majus) is an edible flower plant that offers many health benefits, as it contains numerous medicinally important ingredients, including anthocyanins. By a comparative examination of the possible anthocyanin biosynthesis regulator genes in nasturtium varieties with different anthocyanin contents, we found that TmPAP2, an R2R3 MYB transcription factor gene, is highly expressed in “Empress of India”, a nasturtium variety with high anthocyanin content, while the expression of TmPAP2 in Arabidopsis led to the overproduction of anthocyanins. Protoplast transfection shows that TmPAP2 functions as a transcription activator; consistent with this finding, some of the biosynthesis genes in the general phenylpropanoid pathway and anthocyanin biosynthesis pathway were highly expressed in “Empress of India” and the 35S:TmPAP2 transgenic Arabidopsis plants. However, protoplast transfection indicates that TmPAP2 may not be able to form an MBW complex with TmGL3 and TmTTG1. These results suggest that TmPAP2 may function alone as a key regulator of anthocyanin biosynthesis in nasturtiums.

RsTTG1, a WD40 Protein, Interacts with the bHLH Transcription Factor RsTT8 to Regulate Anthocyanin and Proanthocyanidin Biosynthesis in Raphanus sativus

MBW complexes, consisting of MYB, basic helix–loop–helix (bHLH), and WD40 proteins, regulate multiple traits in plants, including anthocyanin and proanthocyanidin (PA) biosynthesis and the determination of epidermal cell fate. Here, a WD40 gene from Raphanus sativus, designated TRANSPARENT TESTA GLABRA 1 (RsTTG1), was cloned and functionally characterized. Heterologous expression of RsTTG1 in the Arabidopsis thaliana mutant ttg1-22 background restored accumulation of anthocyanin and PA in the mutant and rescued trichome development. In radish, RsTTG1 was abundantly expressed in all root and leaf tissues, independently of anthocyanin accumulation, while its MBW partners RsMYB1 and TRANSPARENT TESTA 8 (RsTT8) were expressed at higher levels in pigment-accumulating tissues. In yeast two-hybrid analysis, the full-length RsTTG1 protein interacted with RsTT8. Moreover, transient protoplast co-expression assays demonstrated that RsTTG1, which localized to both the cytoplasm and nucleus, moves from the cytoplasm to the nucleus in the presence of RsTT8. When co-expressed with RsMYB1 and RsTT8, RsTTG1 stably activated the promoters of the anthocyanin biosynthesis genes CHALCONE SYNTHASE (RsCHS) and DIHYDROFLAVONOL 4-REDUCTASE (RsDFR). Transient expression of RsTTG1 in tobacco leaves exhibited an increase in anthocyanin accumulation due to activation of the expression of anthocyanin biosynthesis genes when simultaneously expressed with RsMYB1 and RsTT8. These results indicate that RsTTG1 is a vital regulator of pigmentation and trichome development as a functional homolog of AtTTG1.

Methylation level of potato gene OMT30376 regulates tuber anthocyanin transformations

After anthocyanin synthesis, a variety of anthocyanin compounds are produced through further methylation, glycosylation, and acylation. However, the effect of the potato methylase gene on anthocyanin biosynthesis has not been reported. Red and purple mutation types appear in tubers of the potato cultivar 'Purple Viking' with chimeric skin phenotypes. In this study, transcriptome and anthocyanin metabolome analyses were performed on skin of Purple Viking tubers and associated mutants. According to the metabolome analysis, the transformation of delphinidin into malvidin-3-O-glucoside and petunidin 3-O-glucoside and that of cyanidin into rosinidin O-hexoside and peonidin-3-O-glucoside were hindered in red tubers. Expression of methyltransferase gene OMT30376 was significantly lower in red tubers than in purple ones, whereas the methylation level of OMT30376 was significantly higher in red tubers. In addition, red skin appeared in tubers from purple tuber plants treated with S-adenosylmethionine (SAM), indicating the difference between purple and red was caused by the methylation degree of the gene OMT30376. Thus, the results of the study suggest that the OMT30376 gene is involved in the transformation of anthocyanins in potato tubers. The results also provide an important reference to reveal the regulatory mechanisms of anthocyanin biosynthesis and transformation.

Genome-wide transcriptome analysis reveals the regulatory network governing terpene trilactones biosynthesis in Ginkgo biloba

Ginkgo biloba L. is currently the only remaining gymnosperm of the Ginkgoaceae Ginkgo genus, and its history can be traced back to the Carboniferous 200 million years ago. Terpene trilactones (TTLs) are one of the main active ingredients in G. biloba, including ginkgolides and bilobalide. They have a good curative effect on cardiovascular and cerebrovascular diseases because of their special antagonistic effect on platelet-activating factors. Therefore, it is necessary to deeply mine genes related to TTLs and to analyze their transcriptional regulation mechanism, which will hold vitally important scientific and practical significance for quality improvement and regulation of G. biloba. In this study, we performed RNA-Seq on the root, stem, immature leaf, mature leaf, microstrobilus, ovulate strobilus, immature fruit and mature fruit of G. biloba. The TTL regulatory network of G. biloba in different organs was revealed by different transcriptomic analysis strategies. Weighted gene co-expression network analysis (WGCNA) revealed that the five modules were closely correlated with organs. The 12 transcription factors, 5 structural genes and 24 Cytochrome P450 (CYP450) were identified as candidate regulators for TTL accumulation by WGCNA and cytoscape visualization. Finally, 6 APETALA2/ethylene response factors, 2 CYP450s and bHLH were inferred to regulate the metabolism of TTLs by correlation analysis. This study is the comprehensive in authenticating transcription factors, structural genes and CYP450 involved in TTL biosynthesis, thereby providing molecular evidence for revealing the comprehensive regulatory network involved in TTL metabolism in G. biloba.

Bioactivity and Therapeutic Potential of Kaempferol and Quercetin: New Insights for Plant and Human Health

Plant secondary metabolites, especially flavonoids, are major metabolites widely found in plants that play several key roles in plant defence and signalling in response to stress conditions. The most studied among these flavonoids are kaempferol and quercetin due to their anti-oxidative potential and their key roles in the defence system, making them more critical for plant adaptation in stress environments. Kaempferol and quercetin in plants have great therapeutic potential for human health. Despite being well-studied, some of their functional aspects regarding plants and human health need further evaluation. This review summarizes the emerging potential of kaempferol and quercetin in terms of antimicrobial activity, bioavailability and bioactivity in the human body as well as in the regulation of plant defence in response to stresses and as a signalling molecule in terms of hormonal modulation under stress conditions. We also evaluated the safe use of both metabolites in the pharmaceutical industry.

Multi-omics analysis reveals the molecular changes accompanying heavy-grazing-induced dwarfing of Stipa grandis

Heavy grazing significantly reduces Stipa grandis growth. To enhance our understanding of plant responses to heavy grazing, we conducted transcriptomic, proteomic, and metabolic analyses of the leaves of non-grazed plants (NG) and heavy-grazing-induced dwarf plants (HG) of S. grandis. A total of 101 metabolites, 167 proteins, and 1,268 genes differed in abundance between the HG and NG groups. Analysis of Kyoto Encyclopedia of Genes and Genomes pathways among differentially accumulated metabolites (DAMs) revealed that the most enriched pathways were flavone and flavonol biosynthesis, tryptophan metabolism, and phenylpropanoid biosynthesis. An integrative analysis of differentially expressed genes (DEGs) and proteins, and DAMs in these three pathways was performed. Heavy-grazing-induced dwarfism decreased the accumulation of DAMs enriched in phenylpropanoid biosynthesis, among which four DAMs were associated with lignin biosynthesis. In contrast, all DAMs enriched in flavone and flavonol biosynthesis and tryptophan metabolism showed increased accumulation in HG compared with NG plants. Among the DAMs enriched in tryptophan metabolism, three were involved in tryptophan-dependent IAA biosynthesis. Some of the DEGs and proteins enriched in these pathways showed different expression trends. The results indicated that these pathways play important roles in the regulation of growth and grazing-associated stress adaptions of S. grandis. This study enriches the knowledge of the mechanism of heavy-grazing-induced growth inhibition of S. grandis and provides valuable information for restoration of the productivity in degraded grassland.

Integration of metabolome and transcriptome analyses reveals the mechanism of anthocyanin accumulation in purple radish leaves

Anthocyanins are natural pigments and play significant roles in multiple growth, development, and stress response processes in plants. The vegetables with high anthocyanin content have better colours, higher antioxidant activity than green vegetables and are potent antioxidants with health benefits. However, the mechanism of anthocyanin accumulation in purple and green leaves of Raphanus sativus (radish) is poorly understood and needs further investigation. In the present study, the pigment content in a green leaf cultivar "RA9" and a purple-leaf cultivar "MU17" was characterized and revealed that the MU17 had significantly increased accumulation of anthocyanins and reduced content of chlorophyll and carotenoid compared with that in RA9. Meanwhile, these two cultivars were subjected to a combination of metabolomic and transcriptome studies. A total of 52 massively content-changed metabolites and 3463 differentially expressed genes were discovered in MU17 compared with RA9. In addition, the content of significantly increased flavonoids (such as pelargonidin and cyanidin) was identified in MU17 compared to RA9 using an integrated analysis of metabolic and transcriptome data. Moreover, the quantitative real-time polymerase chain reaction results also confirmed the differences in the expression of genes related to pathways of flavonoids and anthocyanin metabolism in MU17 leaves. The present findings provide valuable information for anthocyanin metabolism and further genetic manipulation of anthocyanin biosynthesis in radish leaves.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01245-w.

Wall-associated kinase-like gene RL1 contributes to red leaves in sorghum

Red leaves are common in trees but rare in cereal crops. Red leaves can be used as raw materials for anthocyanin extraction and may have some adaptive significance for plants. In this study, we discovered a red leaf phenotype in the F₁ hybrids derived from a cross between two sorghum accessions with green leaf. Histological analysis of red leaves and green leaves showed that red compounds accumulate in mesophyll cells and gradually spreads to the entire leaf blade. In addition, we found chloroplasts degraded more quickly in red leaves than in green leaves based on transmission electron microscopy. Metabolic analysis revealed that flavonoids including six anthocyanins are more abundant in red leaves. Moreover, transcriptome analysis revealed that expression of flavonoid biosynthesis genes was upregulated in red leaves. These observations indicate that flavonoids and anthocyanins in particular, are ideal candidates for the red compounds accumulating in red leaves. Segregation analysis of the red leaf phenotype suggested a genetic architecture consisting of three dominant genes, one (RL1 for RED LEAF1) of which we mapped to a 55-kb region on chromosome 7 containing seven genes. Sequencing, reverse transcription-polymerase chain reaction, and transcriptome analysis suggested Sobic.007G214300, encoding a wall-associated kinase, as the most likely candidate for RL1. Fine mapping the red leaf gene and identifying the metabolites that cause red leaf in sorghum provide us with a better understanding of the red leaf phenotype in the natural population of sorghum.

Integrated analysis of multi-omics and fine-mapping reveals a candidate gene regulating pericarp color and flavonoids accumulation in wax gourd (Benincasa hispida)

Wax gourd (Benincasa hispida), a popular fruit of the Cucurbitaceae (cucurbits) family, contains many nutrients with health benefits and is widely grown in China and other tropical areas. In this study, a wax gourd mutant hfc12 with light-color pericarp was obtained through ethane methylsulfonate (EMS) mutagenesis. Integrative analysis of the metabolome and transcriptome identified 31 differentially accumulated flavonoids (DAFs; flavonoids or flavonoid glycosides) and 828 differentially expressed genes (DEGs) between the hfc12 mutant and wild-type 'BWT'. Furthermore, BSA-seq and kompetitive allele specific PCR (KASP) analysis suggested that the light-color pericarp and higher flavonoid content was controlled by a single gene BhiPRR6 (Bhi12M000742), a typical two-component system (TCS) pseudo-response regulator (PRR). Genetic analysis detected only one nonsynonymous mutation (C-T) in the second exon region of the BhiPRR6. Weighted correlation network analysis (WGCNA) identified the downstream target genes of BhiPRR6, probably regulated by light and were intermediated in the regulatory enzyme reaction of flavonoid biosynthetic pathway. Thus, these results speculated that the transcription factor BhiPRR6, interacting with multiple genes, regulates the absorption of light signals and thereby changes the pericarp color and synthesis of flavonoids in wax gourd.

Integrative Analysis of the Metabolome and Transcriptome Provides Insights into the Mechanisms of Flavonoid Biosynthesis in Quinoa Seeds at Different Developmental Stages

Quinoa (Chenopodium quinoa Willd.) is a crop with high nutritional and health benefits. Quinoa seeds are rich in flavonoid compounds; however, the mechanisms behind quinoa flavonoid biosynthesis remain unclear. We independently selected the high-generation quinoa strain ‘Dianli-3260′, and used its seeds at the filling, milk ripening, wax ripening, and mature stages for extensive targeted metabolome analysis combined with joint transcriptome analysis. The results showed that the molecular mechanism of flavonoid biosynthesis in quinoa seeds was mainly concentrated in two pathways: “flavonoid biosynthesis pathway” and “flavone and flavonol biosynthesis pathway”. Totally, 154 flavonoid-related metabolites, mainly flavones and flavonols, were detected in the four development stages. Moreover, 39,738 genes were annotated with KEGG functions, and most structural genes of flavonoid biosynthesis were differentially expressed during grain development. We analyzed the differential flavonoid metabolites and transcriptome changes between the four development stages of quinoa seeds and found that 11 differential flavonoid metabolites and 22 differential genes were the key factors for the difference in flavonoid biosynthesis. This study provides important information on the mechanisms underlying quinoa flavonoid biosynthesis, the screening of potential quinoa flavonoid biosynthesis regulation target genes, and the development of quinoa products.

Transcriptome analysis reveals chrysanthemum flower discoloration under high-temperature stress

Temperature is an important environmental factor affecting plant anthocyanin synthesis. High temperatures are associated with decreased anthocyanin pigmentation in chrysanthemum. To reveal the effects of high temperature on anthocyanin biosynthesis in chrysanthemum, ray florets of the heat-sensitive cultivar "Nannong Ziyunying" (ZYY) were subjected to RNA sequencing. A total of 18,286 unigenes were differentially expressed between the control and treatment groups. Functional annotation and enrichment analyses of these unigenes revealed that the heat shock response and flavonoid pathways were significantly enriched, suggesting that the expression of these genes in response to high temperature is associated with the fading of chrysanthemum flower color. In addition, genes related to anthocyanin synthesis and heat shock response were differentially expressed under high-temperature stress. Finally, to further investigate the molecular mechanism of discoloration under high-temperature stress and facilitate the use of marker-assisted breeding for developing novel heat-tolerant cultivars, these results were used to mine candidate genes by analyzing changes in their transcription levels in chrysanthemum.

Transcriptional regulation of proanthocyanidin biosynthesis pathway genes and transcription factors in Indigofera stachyodes Lindl. roots

BACKGROUND

Proanthocyanidins (PAs) have always been considered as important medicinal value component. In order to gain insights into the PA biosynthesis regulatory network in I. stachyodes roots, we analyzed the transcriptome of the I. stachyodes in Leaf, Stem, RootI (one-year-old root), and RootII (two-year-old root).

RESULTS

In this study, a total of 110,779 non-redundant unigenes were obtained, of which 63,863 could be functionally annotated. Simultaneously, 75 structural genes that regulate PA biosynthesis were identified, of these 6 structural genes (IsF3'H1, IsANR2, IsLAR2, IsUGT72L1-3, IsMATE2, IsMATE3) may play an important role in the synthesis of PAs in I. stachyodes roots. Furthermore, co-expression network analysis revealed that 34 IsMYBs, 18 IsbHLHs, 15 IsWRKYs, 9 IsMADSs, and 3 IsWIPs hub TFs are potential regulators for PA accumulation. Among them, IsMYB24 and IsMYB79 may be closely involved in the PA biosynthesis in I. stachyodes roots.

CONCLUSIONS

The biosynthesis of PAs in I. stachyodes roots is mainly produced by the subsequent pathway of cyanidin. Our work provides new insights into the molecular pathways underlying PA accumulation and enhances our global understanding of transcriptome dynamics throughout different tissues.

Rice β-Glucosidase 4 (Os1βGlu4) Regulates the Hull Pigmentation via Accumulation of Salicylic Acid

Salicylic acid (SA) is a stress hormone synthesized in phenylalanine ammonia-lyase (PAL) and the branching acid pathway. SA has two interconvertible forms in plants: SAG (SA O-β-glucoside) and SA (free form). The molecular mechanism of conversion of SA to SAG had been reported previously. However, which genes regulate SAG to SA remained unknown. Here, we report a cytoplasmic β-glucosidase (β-Glu) which participates in the SA pathway and is involved in the brown hull pigmentation in rice grain. In the current study, an EMS-generated mutant brown hull 1 (bh1) displayed decreased contents of SA in hulls, a lower photosynthesis rate, and high-temperature sensitivity compared to the wild type (WT). A plaque-like phenotype (brown pigmentation) was present on the hulls of bh1, which causes a significant decrease in the seed setting rate. Genetic analysis revealed a mutation in LOC_Os01g67220, which encodes a cytoplasmic Os1βGlu4. The knock-out lines displayed the phenotype of brown pigmentation on hulls and decreased seed setting rate comparable with bh1. Overexpression and complementation lines of Os1βGlu4 restored the phenotype of hulls and normal seed setting rate comparable with WT. Subcellular localization revealed that the protein of Os1βGlu4 was localized in the cytoplasm. In contrast to WT, bh1 could not hydrolyze SAG into SA in vivo. Together, our results revealed the novel role of Os1βGlu4 in the accumulation of flavonoids in hulls by regulating the level of free SA in the cellular pool.

Decoding the formation of diverse petal colors of Lagerstroemia indica by integrating the data from transcriptome and metabolome

Lagerstroemia indica has great economic value due to its ecological, medicinal, and ornamental properties. Because its bloom color is one of the most essential characteristics, research into its color development is a hot topic. In this study, five representative colored cultivars were chosen, each representing a different color, such as white, red, pink, violet, and purple. Fully bloomed flowers were used to detect flavonoids in the petals. Anthocyanin is the main factor for the color formation of L. indica. 14 anthocyanins were discovered among the 299 flavonoids. Among 14 anthocyanins, malvidin-3,5-di-O-glucoside varied greatly among four colored samples and is the main contributor to color diversity. Transcriptome sequencing revealed that compared to white flowers, Anthocyanin pathway genes appear to be more active in colored samples. Analyzing the correlation network between metabolites and differential expressed genes, 53 key structural genes, and 24 TFs were detected that may play an essential role in the formation of color in L. indica flowers. Among these, the differential expression of F3'5'H and F3'H between all samples are contributors to color diversity. These findings lay the foundation for discovering the molecular mechanism of L. indica flower color diversity.

The MdBBX22-miR858-MdMYB9/11/12 module regulates proanthocyanidin biosynthesis in apple peel

Proanthocyanidins (PAs) have antioxidant properties and are beneficial to human health. The fruit of apple (Malus × domestica Borkh.), especially the peel, is rich in various flavonoids, such as PAs, and thus is an important source of dietary antioxidants. Previous research on the regulation of PAs in apple has mainly focussed on the transcription level, whereas studies conducted at the post-transcriptional level are relatively rare. In this study, we investigated the function of mdm-miR858, a miRNA with multiple functions in plant development, in the peel of apple fruit. We showed that mdm-miR858 negatively regulated PA accumulation by targeting MdMYB9/11/12 in the peel. During fruit development, mdm-miR858 expression was negatively correlated with MdMYB9/11/12 expression and PA accumulation. A 5'-RACE experiment, GUS staining assays and transient luminescent assays indicated that mdm-miR858 cleaved and inhibited the expression of MdMYB9/11/12. Overexpression of mdm-miR858 in apple calli, tobacco and Arabidopsis reduced the accumulation of PAs induced by overexpression of MdMYB9/11/12. Furthermore, we found that MdBBX22 bound to the mdm-miR858 promoter and induced its expression. Overexpression of MdBBX22 induced the expression of mdm-miR858 to inhibit the accumulation of PAs in apple calli overexpressing MdMYB9/11/12. Under light stress, MdBBX22 induced mdm-miR858 expression to inhibit PA accumulation and thereby indirectly enhanced anthocyanin synthesis in the peel. The present results revealed that the MdBBX22-miR858-MdMYB9/11/12 module regulates PA accumulation in apple. The findings provide a reference for further studies of the regulatory mechanism of PA accumulation and the relationship between PAs and anthocyanins.

Comparative transcriptomics of Pinus massoniana organs provides insights on terpene biosynthesis regulation

Terpenoids are the most important natural products collected from conifer species. However, the molecular mechanisms and core factors underlying terpenoid biosynthesis in Pinus massoniana remain unclear. To clarify these mechanisms, this study aimed to identify potential genes that might participate in the terpenoid biosynthesis of P. massoniana. In this study, single molecule real-time (SMRT) sequencing and expression analysis were used to confirm the expression patterns of genes involved in the cones, immature needles, mature needles, immature branches, and mature branches of P. massoniana. A total of 31,331 lncRNAs and 71,240 mRNAs were identified from these organs, and the greatest number of differentially expressed genes (DEGs) was discovered between needles and branches. Weighted gene coexpression network analysis (WGCNA) classified all expressed genes into nine typical modules with 11 kinds of transcription factors (TFs), namely, AP2-ERF, ARF, AUX-IAA, C2H2, Dof, F-box, SBP, WRKY, bHLH, bZIP, and GRAS, and seven kinds of functional genes, namely, ABC transporter, cellulose synthase (CesA), leucine-rich repeats (LRR), cytochrome P450 (CYT P450), pathogenesis-related protein (PR), terpene synthase (TPS), and chlorophyllase enzyme. A molecular network was constructed for hub genes, TFs, and functional genes in three modules. The potential function of eight candidate genes, including PmbHLH2, PmERF1, PmRGA, PmGAI, PmbZIP1, PmLOB1, PmMADS1, and PmMYB1, was validated through correlation analysis between terpenoid contents and expression levels, subcellular localization, and transcriptional activation activity, which provides us with probable regulators of terpenoid biosynthesis in conifers.

Anthocyanins identification and transcriptional regulation of anthocyanin biosynthesis in purple Brassica napus

The main anthocyanin components were identified, and the transcriptional regulation pattern of anthocyanin related genes in leaves and stem bark was elucidated in a purple B. napus. Brassica napus is one of the most important oil crops planted worldwide, and developing varieties of dual-purpose for oil and vegetable is beneficial to improve economic benefits. Anthocyanins are a class of secondary metabolites that not only make plants present beautiful colors, but have a variety of important physiological functions and biological activities. Therefore, increasing the accumulation of anthocyanin in vegetative organs can improve vegetable value of rapeseed. However, anthocyanin enriched varieties in vegetative organs are rare, and there are few studies on category identification and accumulation mechanism of anthocyanin, which limits the utilization of anthocyanins in B. napus. In this study, 157 anthocyanin biosynthesis related genes (ABGs) were identified in B. napus genome by homology comparison and collinearity analysis of genes related to anthocyanin synthesis and regulation in Arabidopsis. Moreover, five anthocyanins were identified in the stem bark and leaves of the purple B. napus PR01 by high performance liquid chromatography-mass spectrometry (HPLC-MS), and the expression characteristics of ABGs in the leaves and stem bark of PR01 were analyzed and compared with the green cultivar ZS11 by RNA-Seq. Combining further weighted gene co-expression network analysis (WGCNA), the up-regulation of transcript factors BnaA07. PAP2 and BnaC06. PAP2 were identified as the key to the up-regulation of most of anthocyanin synthesis genes that promoted anthocyanin accumulation in PR01. This study is helpful to understand the transcriptional regulation of anthocyanin biosynthesis in B. napus and provides the theoretical basis for breeding novel varieties of dual-purpose for oil and vegetable.

Transcriptome profiles of yellowish-white and fuchsia colored flowers in the Rheum palmatum complex reveal genes related to color polymorphism

Flower color variation is ubiquitous in many plant species, and several studies have been conducted to elucidate the underlying molecular mechanism. There are two flower color variants (yellowish-white and fuchsia) in the Rheum palmatum complex, however, few studies have investigated this phenomenon. Here, we used transcriptome sequencing of the two color variants to shed light on the molecular and biochemical basis for these color morphs. Comparison of the two transcriptomes identified 9641 differentially expressed unigenes (DEGs), including 6477 up-regulated and 3163 down-regulated genes. Functional analyses indicated that several DEGs were related to the anthocyanin biosynthesis pathway, and the expression profiles of these DEGs were coincident with the qRT-PCR validation results, indicating that expression levels of structural genes have a profound effect on the color variation in the R. palmatum complex. Our results suggested that the interaction of transcription factors (MYB, bHLH and WRKY) also regulated the anthocyanin biosynthesis in the R. palmatum complex. Estimation of selection pressures using the dN/dS ratio showed that 1106 pairs of orthologous genes have undergone positive selection. Of these positively selected genes, 21 were involved in the anthocyanin biosynthetic pathway, indicating that they may encode the proteins for structural alteration and affect flower color in the R. palmatum complex.

Transcriptomic and metabolomic profiling of flavonoid biosynthesis provides novel insights into petals coloration in Asian cotton (Gossypium arboreum L.)

BACKGROUND

Asian cotton (Gossypium arboreum L.), as a precious germplasm resource of cotton with insect resistance and stress tolerance, possesses a broad spectrum of phenotypic variation related to pigmentation. Flower color affects insect pollination and the ornamental value of plants. Studying flower color of Asian cotton varieties improves the rate of hybridization and thus enriches the diversity of germplasm resources. Meanwhile, it also impacts the development of the horticultural industry. Unfortunately, there is a clear lack of studies concerning intricate mechanisms of cotton flower-color differentiation. Hereby, we report an integrative approach utilizing transcriptome and metabolome concerning flower color variation in three Gossypium arboreum cultivars.

RESULTS

A total of 215 differentially accumulated metabolites (DAMs) were identified, including 83 differentially accumulated flavonoids (DAFs). Colorless kaempferol was more abundant in white flowers, while gossypetin-fer showed specificity in white flowers. Quercetin and gossypetin were the main contributors to yellow petal formation. Pelargonidin 3-O-beta-D-glucoside and cyanidin-3-O-(6''-Malonylglucoside) showed high accumulation levels in purple petals. Quercetin and gossypetin pigments also promoted purple flower coloration. Moreover, 8178 differentially expressed genes (DEGs) were identified by RNA sequencing. The correlation results between total anthocyanins and DEGs were explored, indicating that 10 key structural genes and 29 transcription factors promoted anthocyanin biosynthesis and could be candidates for anthocyanin accumulation. Ultimately, we constructed co-expression networks of key DAFs and DEGs and demonstrated the interactions between specific metabolites and transcripts in different color flowers.

CONCLUSION

This study provides new insights into elucidating the regulatory mechanisms of cotton flower color and lays a potential foundation for generate cotton varieties with highly attractive flowers for pollinators.

Integrated metabolomics and transcriptomic analysis of the flavonoid regulatory networks in Sorghum bicolor seeds

Background

The objective of this study was to reveal the flavonoid biosynthesis pathway in white (Z6), red (Z27) and black (HC4) seeds of the sweet sorghum (Sorghum bicolor) using metabolomics and transcriptomics, to identify different flavonoid metabolites, and to analyze the differentially expressed genes involved in flavonoid biosynthesis.

Results

We analyzed the metabolomics and transcriptomics data of sweet sorghum seeds. Six hundred and fifty-one metabolites including 171 flavonoids were identified in three samples. Integrated analysis of transcriptomics and metabolomics showed that 8 chalcone synthase genes (gene19114, gene19115, gene19116, gene19117, gene19118, gene19120, gene19122 and gene19123) involved in flavonoid biosynthesis, were identified and play central role in change of color. Six flavanone including homoeriodictyol, naringin, prunin, naringenin, hesperetin and pinocembrin were main reason for the color difference.

Conclusions

Our results provide valuable information on the flavonoid metabolites and the candidate genes involved in the flavonoid biosynthesis pathway in sweet sorghum seeds.