Integrative analysis of transcriptome and metabolome reveals flavonoid biosynthesis regulation in Rhododendron pulchrum petals

Comprehensive transcriptome and metabolome analysis deciphers the mechanism underlying rapid xylem growth in the dominant hybrid poplar QB3

MAIN CONCLUSION

Compared with its parents, the heterosis in growth of QB3 is primarily attributed to the upregulation of auxin and brassinosteroid-related genes, as well as the induced expression of numerous xylem and phloem synthesis genes, particularly the accumulation of lignin. Interestingly, QB3 significantly increased resistance to gray mold, which may be related to anthocyanin accumulation. Our findings illuminate the complex interplay of biological mechanisms that govern the regulation of wood growth and resistance. Poplar, as a fast-growing energy species widely distributed in the northern hemisphere, has important ecological and economic value. The hybridization of poplars is very common and often can bring to the progeny superior growth and resilience traits, but the molecular mechanism of heterosis remains to be studied. Through decades of crossbreeding work, a high-growth rate hybrid offspring named QinBai3 (QB3) was selected from P. alba × (P. alba × P. glandulosa), which provided an ideal model for investigating the molecular mechanism of heterosis. We found that the plant height, ground diameter, and xylem thickness of QB3 were much higher than those of I101 and 84 K. Through transcriptome and qRT-PCR analyses, we found that the expression levels of poplar regulatory genes associated with vegetative growth, brassinosteroid (BR), and auxin hormone signaling were significantly elevated in July compared to February. Meanwhile, compared to its parents, QB3 exhibited more specifically up-regulated genes in the processes of xylem and phloem synthesis, notably PalOPS and PalPRX52. However, in response to certain abiotic stresses, such as water deprivation and UV-B exposure, more down-regulated genes were identified. Metabolome analyses indicated that QB3 significantly increased the levels of lignin and anthocyanin, a result that aligns with the transcriptome data. Additionally, chemical assays confirmed the substantial accumulation of lignin and anthocyanin in QB3, suggesting that increased lignin accumulation may enhance the stem growth rate of QB3. Surprisingly, QB3 significantly increased resistance to Botrytis cinerea B05.10, which was accompanied by anthocyanin accumulation. In addition, our study offers detailed insights into the molecular mechanisms underlying rapid growth and stress resistance in hybrid poplar, thereby providing a new theoretical foundation and practical guidance for forest genetic breeding.

Joint analysis of transcriptional metabolism for flavonoid synthesis during different developmental periods in oil palm exocarp

To identify candidate genes for breeding oil palm varieties with high flavonoid content through molecular biotechnology, this study analyzed the metabolomes and transcriptomes of oil palm exocarp at different developmental stages using LC-MS/MS and RNA-Seq techniques. The green fruiting type (FS) oil palm exocarp at 95 days (FS1), 125 days (FS2), and 185 days (FS3) after pollination served as the materials. The enzyme genes F3H, CHS, ANS, and DFR were positively correlated with Quercetin-3-O-sambubioside. DFR also showed positive correlations with Afzelechin, Epiafzelechin, and Baimaside. In contrast, F3H, CHS, and ANS were negatively correlated with Hesperetin-7-O-glucoside. Additionally, CYP73A, UGT73C6, FG2-1, and FG2-2 were negatively correlated with Afzelechin, Epiafzelechin, Quercetin-3-O-sambubioside, and Baimaside, while CYP75A was negatively correlated with Epiafzelechin, Quercetin-3-O-sambubioside, and Baimaside. These results suggest that F3H, CHS, ANS, and DFR play a role in promoting Quercetin-3-O-sambubioside* synthesis, with DFR further enhancing the production of Afzelechin, Epiafzelechin, and Baimaside. On the other hand, F3H, CHS, and ANS may inhibit Hesperetin-7-O-glucoside synthesis. Meanwhile, CYP73A, UGT73C6, FG2-1, and FG2-2 appear to suppress the synthesis of multiple flavonoids, including Afzelechin, Epiafzelechin, Quercetin-3-O-sambubioside*, and Baimaside. Lastly, CYP75A is implicated in suppressing Epiafzelechin, Quercetin-3-O-sambubioside*, and Baimaside synthesis. These findings provide a foundation for future molecular breeding efforts targeting flavonoid-rich oil palm varieties.

Integrated Metabolomics and Transcriptome Analysis of Anthocyanin Biosynthetic Pathway in Prunus serrulata

Prunus serrulata is an important landscape tree species whose flower color has high ornamental value. However, the molecular mechanisms regulating flower color in P. serrulata remain unclear. By studying the metabolomics and transcriptomics of three different color varieties under the species lineage of P. serrulata, 'Eigeng' (EG, white), 'Albo-rosea' (AR, pink), and 'Grandiflora' (GF, green), the biosynthetic mechanisms of different flower colors in P. serrulata were revealed. The results showed that the different colors of the petals were related to the content of chlorophyll and anthocyanins. Among these, cyanidin-3-O-glucoside and cyanidin-3-O-(6-O-malonyl-β-D-glucoside) were highly expressed in AR. A combined transcriptomic analysis revealed that five flavonoid structural genes, including two DFR genes and three UFGT genes, were specifically expressed. In addition, three key transcription factors, PsMYB77, PsMYB17, and PsMYB105, were identified as regulators of the structural genes DFR and UFGT and participants in the forward synthesis of anthocyanin. This study provides convincing evidence elucidating the regulatory mechanisms of anthocyanin synthesis of P. serrulata and provides a theoretical basis for the breeding and development of new varieties and germplasm resource innovation for cherry blossom.

Comparative global profiling of Perilla leaf and stem via transcriptomics and metabolomics

Perilla (Perilla frutescens L.) is a time-honored herbal plant with widespread applications in both medicine and culinary practices around the world. Profiling the essential organs and tissues with medicinal significance on a global scale offers valuable insights for enhancing the yield of desirable compounds in Perilla and other medicinal plants. In the present study, genome-wide RNA-sequencing (RNA-seq) and assessing the global spectrum of metabolites were carried out in the two major organs/tissues of stem (PfST) and leaf (PfLE) in Perilla. The results showed a total of 18,490 transcripts as the DEGs (differentially expressed genes) and 144 metabolites as the DAMs (differentially accumulated metabolites) through the comparative profiling of PfST vs PfLE, and all the DEGs and DAMs exhibited tissue-specific trends. An association analysis between the transcriptomics and metabolomics revealed 14 significantly enriched pathways for both DEGs and DAMs, among which the pathways of Glycine, serine and threonine metabolism (ko00260), Glyoxylate and dicarboxylate metabolism (ko00630), and Glucagon signaling pathway (ko04922) involved relatively more DEGs and DAMs. The results of qRT-PCR assays of 18 selected DEGs confirmed the distinct tissue-specific characteristics of all identified DEGs between PfST and PfLE. Notably, all eight genes associated with the flavonoid biosynthesis/metabolism pathways exhibited significantly elevated expression levels in PfLE compared to PfST. This observation suggests a heightened accumulation of metabolites related to flavonoids in Perilla leaves. The findings of this study offer a comprehensive overview of the organs and tissues in Perilla that have medicinal significance.

Integrated Transcriptome and Metabolome Analysis Reveals Mechanism of Flavonoid Synthesis During Low-Temperature Storage of Sweet Corn Kernels

Sweet corn is a globally important food source and vegetable renowned for its rich nutritional content. However, post-harvest quality deterioration remains a significant challenge due to sweet corn's high sensitivity to environmental factors. Currently, low-temperature storage is the primary method for preserving sweet corn; however, the molecular mechanisms involved in this process remain unclear. In this study, kernels stored at different temperatures (28 °C and 4 °C) for 1, 3, and 5 days after harvest were collected for physiological and transcriptomic analysis. Low temperature storage significantly improved the PPO and SOD activity in sweet corn kernels compared to storage at a normal temperature. A total of 1993 common differentially expressed genes (DEGs) were identified in kernels stored at low temperatures across all three time points. Integrated analysis of transcriptomic and previous metabolomic data revealed that low temperature storage significantly affected flavonoid biosynthesis. Furthermore, 11 genes involved in flavonoid biosynthesis exhibited differential expression across the three storage periods, including CHI, HCT, ANS, F3'H, F3'5'H, FLS, and NOMT, with Eriodictyol, Myricetin, and Hesperetin-7-O-glucoside among the key flavonoids. Correlation analysis revealed three AP2/ERF-ERF transcription factors (EREB14, EREB182, and EREB200) as potential regulators of flavonoid biosynthesis during low temperature treatment. These results enhance our understanding of the mechanisms of flavonoid synthesis in sweet corn kernels during low-temperature storage.

Recent advancements in the physiological, genetic, and genomic research on Rhododendrons for trait improvement

High species diversity, hybridization potential, broad geographical dispersal range and ornamental characteristics (i.e., attractive size, shape, structure, flowers, and evergreen) have fetched a good international market for Rhododendron. However, most species are restricted to specific geographic areas due to their habitat specificity in acidic soil and cold climates, resulting many species being classified under threat categories of the IUCN. In this review, advances in research on Rhododendron for improvement to floral display quality and stress resistance have been described. The low genetic barrier among species has created opportunities for extensive hybridization and ploidy alteration for introducing quality and adaptive traits during the development of new varieties. Recent technological advances have supported investigations into the mechanism of flower development, as well as cold tolerance and pathogen resistance mechanisms in the Rhododendron. However, most of the species have limited adaptability to drought, line-tolerance, pathogen resistance, and high-temperature conditions and this resistance ability present in few species largely remains unexplored. Additionally, the available genetic diversity and genomic information on species, and possibilities for their application in molecular breeding have been summarized. Overall, genomic resource data are scarce in the majority of the members of this genus. Finally, various research gaps such as genetic mapping of quality traits, understanding the molecular mechanism of quality-related traits and genomic assortment in Rhododendron members have been discussed in the future perspective section.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-04006-6.

Metabolite analysis reveals flavonoids accumulation during flower development in Rhododendron pulchrum sweet (Ericaceae)

The azalea (Rhododendron simsii Planch.) is an important ornamental woody plant with various medicinal properties due to its phytochemical compositions and components. However little information on the metabolite variation during flower development in Rhododendron has been provided. In our study, a comparative analysis of the flavonoid profile was performed in Rhododendron pulchrum sweet at three stages of flower development, bud (stage 1), partially open flower (stage 2), and full bloom (stage 3). A total of 199 flavonoids, including flavone, flavonol, flavone C-glycosides, flavanone, anthocyanin, and isoflavone were identified. In hierarchical clustering analysis (HCA) and principal component analysis (PCA), the accumulation of flavonoids displayed a clear development stage variation. During flower development, 78 differential accumulated metabolites (DAMs) were identified, and most were enriched to higher levels at the full bloom stage. A total of 11 DAMs including flavone (chrysin, chrysoeriol O-glucuronic acid, and chrysoeriol O-hexosyl-O-pentoside), isoflavone (biochanin A), and flavonol (3,7-di-O-methyl quercetin and isorhamnetin) were significantly altered at three stages. In particular, 3,7-di-O-methyl quercetin was the top increased metabolite during flower development. Furthermore, integrative analyses of metabolomic and transcriptomic were conducted, revealing that the contents of isoflavone, biochanin A, glycitin, and prunetin were correlated with the expression of 2-hydroxyisoflavanone dehydratase (HIDH), which provide insight into the regulatory mechanism that controls isoflavone biosynthesis in R. pulchrum. This study will provide a new reference for increasing desired metabolites effectively by more accurate or appropriate genetic engineering strategies.

Integration of Metabolomic and Transcriptomic Analyses Reveals the Molecular Mechanisms of Flower Color Formation in Prunus mume

Flower color is an important trait that affects the economic value of Prunus mume, a famous ornamental plant in the Rosaceae family. P. mume with purple-red flowers is uniquely charming and highly favored in landscape applications. However, little is known about its flower coloring mechanism, which stands as a critical obstacle on the path to innovative breeding for P. mume flower color. In this study, transcriptomic and targeted metabolomic analyses of purple-red P. mume and white P. mume were performed to elucidate the mechanism of flower color formation. In addition, the expression patterns of key genes were analyzed using an RT-qPCR experiment. The results showed that the differential metabolites were significantly enriched in the flavonoid synthesis pathway. A total of 14 anthocyanins emerged as the pivotal metabolites responsible for the differences in flower color between the two P. mume cultivars, comprising seven cyanidin derivatives, five pelargonium derivatives, and two paeoniflorin derivatives. Moreover, the results clarified that the metabolic pathway determining flower color in purple-red P. mume encompasses two distinct branches: cyanidin and pelargonidin, excluding the delphinidin branch. Additionally, through the integrated analysis of transcriptomic and metabolomic data, we identified 18 key genes responsible for anthocyanin regulation, thereby constructing the gene regulatory network for P. mume anthocyanin synthesis. Among them, ten genes (PmCHI, PmGT2, PmGT5, PmGST3, PmMYB17, PmMYB22, PmMYB23, PmbHLH4, PmbHLH10, and PmbHLH20) related to anthocyanin synthesis were significantly positively correlated with anthocyanin contents, indicating that they may be the key contributors to anthocyanin accumulation. Our investigation contributes a novel perspective to understanding the mechanisms responsible for flower color formation in P. mume. The findings of this study introduce novel strategies for molecular design breeding aimed at manipulating flower color in P. mume.

Integrated transcriptomic and WGCNA analyses reveal candidate genes regulating mainly flavonoid biosynthesis in Litsea coreana var. sinensis

Litsea coreana Levl. var. sinensis (Allen) Yang et P. H. Huang is a popular ethnic herb and beverage plant known for its high flavonoid content, which has been linked to a variety of pharmacological benefits and crucial health-promoting impacts in humans. The progress in understanding the molecular mechanisms of flavonoid accumulation in this plant has been hindered due to the deficiency of genomic and transcriptomic resources. We utilized a combination of Illumina and Oxford Nanopore Technology (ONT) sequencing to generate a de novo hybrid transcriptome assembly. In total, 126,977 unigenes were characterized, out of which 107,977 were successfully annotated in seven public databases. Within the annotated unigenes, 3,781 were categorized into 58 transcription factor families. Furthermore, we investigated the presence of four valuable flavonoids-quercetin-3-O-β-D-galactoside, quercetin-3-O-β-D-glucoside, kaempferol-3-O-β-D-galactoside, and kaempferol-3-O-β-D-glucoside in 98 samples, using high-performance liquid chromatography. A weighted gene co-expression network analysis identified two co-expression modules, MEpink and MEturquoise, that showed strong positive correlation with flavonoid content. Within these modules, four transcription factor genes (R2R3-MYB, NAC, WD40, and ARF) and four key enzyme-encoding genes (CHI, F3H, PAL, and C4H) emerged as potential hub genes. Among them, the R2R3-MYB (LcsMYB123) as a homologous gene to AtMYB123/TT2, was speculated to play a significant role in flavonol biosynthesis based on phylogenetic analysis. Our findings provided a theoretical foundation for further research into the molecular mechanisms of flavonoid biosynthesis. Additionally, The hybrid transcriptome sequences will serve as a valuable molecular resource for the transcriptional annotation of L. coreana var. sinensis, which will contribute to the improvement of high-flavonoid materials.

Integrated metabolomics and transcriptomics reveal molecular mechanisms of corolla coloration in Rhododendron dauricum L

Rhododendron dauricum L. is a semi-evergreen shrub of high ornamental and medicinal values in Northeast China. To study the molecular mechanisms of corolla coloration in R. dauricum, integrated metabolomics and transcriptomics were performed in R. dauricum featuring purple flowers and R. dauricum var. album featuring white flowers. Comparative metabolomics revealed 25 differential metabolites in the corolla of the two distinct colors, enriched in flavonoids that are closely related to pigmentation in the flower. Differential analysis of the transcriptomics data revealed enrichment of structural genes for flavonoid biosynthesis (99 up- and 58 down-regulated, respectively, in purple corollas compared to white ones). Significantly, CHS and CHI, key genes in the early stage of anthocyanin synthesis, as well as F3H, F3'H, F3'5'H, DFR, ANS, and UFGT that promote the accumulation of pigments in the late stage of anthocyanin synthesis, were up-regulated in R. dauricum (purple color). In R. dauricum var. album, FLS were key genes determining the accumulation of flavonols. In addition, transcriptome-metabolome correlation analysis identified 16 R2R3 MYB transcription factors (out of 83 MYBs) that are important for corolla coloration. Five negative and four positive MYBs were further identified by integrated transcriptional and metabolic network analysis, revealing a key role of MYBA and MYB12 in regulating anthocyanins and flavonols, respectively. Moreover, we validated the function of RdMYBA by creating stable transgenic plants and found that RdMYBA promotes anthocyanin biosynthesis. In summary, we systematically characterized the transcriptome and metabolome of two R. dauricum cultivars with different flower colors and identified MYBs as key factors in modulating corolla coloration.

Integrated metabolome, full-length sequencing, and transcriptome analyses unveil the molecular mechanisms of color formation of the canary yellow and red bracts of Bougainvillea × buttiana 'Chitra'

Bougainvillea is a typical tropical flower of great ornamental value due to its colorful bracts. The molecular mechanism behind color formation is not well-understood. Therefore, this research conducted metabolome analysis, transcriptome analysis, and multi-flux full-length sequencing in two color bracts of Bougainvillea × buttiana 'Chitra' to investigate the significantly different metabolites (SDMs) and differentially expressed genes (DEGs). Overall, 261 SDMs, including 62 flavonoids and 26 alkaloids, were detected, and flavonols and betalains were significantly differentially accumulated among the two bracts. Furthermore, the complete-length transcriptome of Bougainvillea × buttiana was also developed, which contained 512 493 non-redundant isoforms. Among them, 341 210 (66.58%) displayed multiple annotations in the KOG, GO, NR, KEGG, Pfam, Swissprot, and NT databases. RNA-seq findings revealed that 3610 DEGs were identified between two bracts. Co-expression analysis demonstrated that the DEGs and SDMs involved in flavonol metabolism (such as CHS, CHI, F3H, FLS, CYP75B1, kaempferol, and quercetin) and betacyanin metabolism (DODA, betanidin, and betacyanins) were the main contributors for the canary yellow and red bract formation, respectively. Further investigation revealed that several putative transcription factors (TFs) might interact with the promoters of the genes mentioned above. The expression profiles of the putative TFs displayed that they may positively and negatively regulate the structural genes' expression profiles. The data revealed a potential regulatory network between important genes, putative TFs, and metabolites in the flavonol and betacyanin biosynthesis of Bougainvillea × buttiana 'Chitra' bracts. These findings will serve as a rich genetic resource for future studies that could create new color bracts.

Transcriptomic Analysis Reveals the Regulatory Mechanism of Color Diversity in Rhododendron pulchrum Sweet (Ericaceae)

Rhododendron pulchrum Sweet is a traditional ornamental plant cultivated in China and presents a great variation in petal coloration. However, few studies have been performed to reveal the genes involved and the regulatory mechanism of flower color formation in this plant. In this study, to explore the underlying genetic basis of flower color formation, transcriptome analysis was performed by high-throughput sequencing techniques on four petal samples of different colors: purple, pink, light pink, and white. Results show that a total of 35.55 to 40.56 million high-quality clean reads were obtained, of which 28.56 to 32.65 million reads were mapped to the reference genome. For their annotation, 28,273, 18,054, 24,301, 19,099, and 11,507 genes were allocated to Nr, Swiss-Prot, Pfam, GO, and KEGG databases, correspondingly. There were differentially expressed genes among the four different petal samples, including signal-transduction-related genes, anthocyanin biosynthesis genes, and transcription factors. We found that the higher expressed levels of genes associated with flavonol synthase (FLS) might be the key to white formation, and the formation of red color may be related to the higher expression of flavanone 4-reductase (DFR) families. Overall, our study provides some valuable information for exploring and understanding the flower color intensity variation in R. pulchrum.

The DREB transcription factor, a biomacromolecule, responds to abiotic stress by regulating the expression of stress-related genes

Abiotic stress is a crucial factor that affects plant survival and growth and even leads to plant death in severe cases. Transcription factors can enhance the ability of plants to fight against various stresses by controlling the expression of downstream genes. The dehydration response element binding protein (DREB) is the most extensive subfamily of AP2/ERF transcription factors involved in abiotic stress. However, insufficient research on the signal network of DREB transcription factors has limited plant growth and reproduction. Furthermore, field planting of DREB transcription factors and their roles under multiple stress also require extensive research. Previous reports on DREB transcription factors have focused on the regulation of DREB expression and its roles in plant abiotic stress. In recent years, there has been new progress in DREB transcription factors. Here, the structure and classification, evolution and regulation, role in abiotic stress, and application in crops of DREB transcription factors were reviewed. And this paper highlighted the evolution of DREB1/CBF, as well as the regulation of DREB transcription factors under the participation of plant hormone signals and the roles of subgroups in abiotic stress. In the future, it will lay a solid foundation for further study of DREB transcription factors and pave the way for the cultivation of resistant plants.

Intracellular metabolomics and microRNAomics unveil new insight into the regulatory network for potential biocontrol mechanism of stress-tolerant Tricho-fusants interacting with phytopathogen Sclerotium rolfsii Sacc

The present study employed microRNA (miRNA) sequencing and metabolome profiling of Trichoderma parental strains and fusants during normal growth and interaction with the phytopathogen Sclerotium rolfsii Sacc. In-vitro antagonism indicated that abiotic stress-tolerant Tricho-fusant FU21 was examined as a potent biocontroller with mycoparasitic action after 10 days. During interaction with the test pathogen, the most abundant uprising intracellular metabolite was recognized as l-proline, which corresponds to held-down l-alanine, associated with arginine and proline metabolism, biosynthesis of secondary metabolites, and nitrogen metabolism linked to predicted genes controlled by miRNAs viz., cel-miR-8210-3p, hsa-miR-3613-5p, and mml-miR-7174-3p. The miRNAs- mml-miR-320c and mmu-miR-6980-5p were found to be associated with phenylpropanoid biosynthesis, transcription factors, and signal transduction pathways, respectively, and were ascertained downregulated in potent FU21_IB compared with FU21_CB. The amino benzoate degradation and T cell receptor signaling pathways were regulated by miRNAs cel-miR-8210 and tca-miR-3824 as stress tolerance mechanisms of FU21. The intracellular metabolites l-proline, maleic acid, d-fructose, Myo-inositol, arabinitol, d-xylose, mannitol, and butane were significantly elevated as potential biocontrol and stress-tolerant constituents associated with miRNA regulatory pathways in potent FU21_IB. A network analysis between regulatory miRNA-predicted genes and intracellular metabolomics acknowledged possible biocontrol pathways/mechanisms in potent FU21_IB to restrain phytopathogen.

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.