Integrated metabolome and transcriptome revealed the flavonoid biosynthetic pathway in developing Vernonia amygdalina leaves

Integrated transcriptomic, metabolomic, and functional analyses unravel the mechanism of bagging delaying fruit cracking of pomegranate (Punica granatum L.)

The economic impact of fruit cracking in pomegranate products is substantial. In this study, we present the inaugural comprehensive analysis of transcriptome and metabolome in the outermost pericarp of pomegranate fruit in bagging conditions. Our investigation revealed a notable upregulation of differentially expressed genes (DEGs) associated with the calcium signaling pathway (76.92%) and xyloglucan endotransglucosylase/hydrolase (XTH) genes (87.50%) in the fruit peel of non-cracking fruit under bagging. Metabolomic analysis revealed that multiple phenolics, flavonoids, and tannins were identified in pomegranate. Among these, calmodulin-like 23 (PgCML23) exhibited a significant correlation with triterpenoids and demonstrated a marked upregulation under bagging treatment. The transgenic tomatoes overexpressing PgCML23 exhibited significantly higher cellulose content and xyloglucan endotransglucosylase (XET) enzyme activity in the pericarp at the red ripening stage compared to the wild type. Conversely, water-soluble pectin content, polygalacturonase (PG), and β-galactosidase (β-GAL) enzyme activities were significantly lower in the transgenic tomatoes. Importantly, the heterologous expression of PgCML23 led to a substantial reduction in the fruit cracking rate in tomatoes. Our findings highlight the reduction of fruit cracking in bagging conditions through the manipulation of PgCML23 expression.

Metabolomics and transcriptomics analysis revealed the response mechanism of alfalfa to combined cold and saline-alkali stress

Cold and saline-alkali stress are frequently encountered by plants, and they often occur simultaneously in saline-alkali soils at mid to high latitudes, constraining forage crop distribution and production. However, the mechanisms by which forage crops respond to the combination of cold and saline-alkali stress remain unknown. Alfalfa (Medicago sativa L.) is one of the most essential forage grasses in the world. In this study, we analyzed the complex response mechanisms of two alfalfa species (Zhaodong [ZD] and Blue Moon [BM]) to combined cold and saline-alkali stress using multi-omics. The results revealed that ZD had a greater ability to tolerate combined stress than BM. The tricarboxylic acid cycles of the two varieties responded positively to the combined stress, with ZD accumulating more sugars, amino acids, and jasmonic acid. The gene expression and flavonoid content of the flavonoid biosynthesis pathway were significantly different between the two varieties. Weighted gene co-expression network analysis and co-expression network analysis based on RNA-Seq data suggested that the MsMYB12 gene may respond to combined stress by regulating the flavonoid biosynthesis pathway. MsMYB12 can directly bind to the promoter of MsFLS13 and promote its expression. Moreover, MsFLS13 overexpression can enhance flavonol accumulation and antioxidant capacity, which can improve combined stress tolerance. These findings provide new insights into improving alfalfa resistance to combined cold and saline-alkali stress, showing that flavonoids are essential for plant resistance to combined stresses, and provide theoretical guidance for future breeding programs.

A chromosome-scale genome sequence of Aeonium(Aeonium arboreum 'Velour') provides novel insights into the evolution of anthocyanin synthesis

Anthocyanin glycoside is a water-soluble flavonoid compound that colors plants and aids in stress resistance. The mechanism driving the evolution of the anthocyanin synthesis pathway in plants remains unclear. Aeonium plants are highly regarded as model organisms for studying adaptive evolution. These plants can be categorized into various types, each distinguished by the content and distribution of anthocyanins in their leaves. The categories include red leaves, green leaves, black leaves, yellow leaves, and a classification known as the 'spot brocade series. In this study, we successfully assembled and annotated the genome of cultivar 'Aeonium arboreum 'Velour'' at chromosomal level. The genome size is 1,334.85 Mb containing 18 chromosomes in a single set, with a contig N50 of 23.47 Mb and a Scaffold N50 of 25.07 Mb. Through homology prediction, de novo prediction, and transcriptome prediction, we identified 166,228 coding genes, 161,656 of which were successfully annotated in the database. Comparative genomic analysis revealed that Aeonium arboreum 'Velour' underwent an independent genome-wide replication event after differentiating from Sedum album, Kalanchoe laxiflora, and Kalanchoe fedtschenkoi. It also shared a genome-wide replication event with Sedum album and Kalanchoe laxiflora. Aeonium arboreum 'Velour' exhibits a higher number of multi-copy gene families compared to other species. A total of 5,129 gene families unique to Aeonium arboreum 'Velour' were identified, primarily enriched in various metabolic pathways, including monoterpenoid biosynthesis, sesquiterpene and triterpene biosynthesis, cyanamide acid metabolism, flavonoid and flavonol biosynthesis, phosphonate and phosphinate metabolism, fatty acid degradation, biosynthesis of unsaturated fatty acid, ether lipid metabolism, tyrosine metabolism, and isoflavone biosynthesis according to the KEGG pathway analysis. Aeonium arboreum 'Velour' and Sedum album diversion dates back to approximately 43.11 million years ago during the Paleogene period, marked by the expansion of 2,807 gene families. In contrast, the divergence from Kalanchoe laxiflora and Kalanchoe fedtschenkoi began around 57.28 million years ago, with 219 gene families expanding. GO analysis highlighted that most of the expansion or contraction gene families were predominantly enriched in flower organs, leaf organ development, anthocyanin metabolism regulation, and light energy absorption and utilization. Remarkably, anthocyanin metabolism regulation is enriched to 80 expanded genes, including 36 bHLH transcription factors, possibly functioning as photosensitive pigment interaction factors (PIFs). We speculate that flavonoids play a pivotal role in the adaptation of Aeonium arboreum 'Velour' to environmental stress. Moreover, the evolution of the anthocyanin synthesis pathway is potentially driven by the plant's capability to absorb and utilize light energy, especially in high CO2 and high-temperature settings characteristic of the early Paleogene.

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.

In vivo anticoccidial activity and immune response of Moringa oleifera and Vernonia amygdalina leaves against Eimeria magna and Eimeria media in rabbits

Introduction

The aims of this study were to determine the immune response and the anticoccidial activity induced by Moringa oleifera and Vernonia amygdalina leaves in rabbits infected with Eimeria magna and Eimeria media.

Methods

Thirty-five-day-old rabbits, free from coccidia, were infested with 2.10³ oocysts of Eimeria magna and Eimeria media, then received the acetone extract of the leaves of Moringa oleifera and Vernonia amygdalina at different doses by oral gavage.

Results and discussion

The inhibition of the excretion of oocysts was evaluated by the McMaster technique and the levels of cytokines (IL-4 and IL-12) and immunoglobulin IgG were assayed by the ELISA method. The in vivo efficacy on E. magna and E. media oocysts was 95.43% and 96.53% for Moringa oleifera and Vernonia amygdalina at 1000 mg/kg bw against 98% for the positive control. Interestingly the plant extracts increased the production of interleukin (IL) and immunoglobulins (Ig) compared to controls. Plasma IL-4 levels (pg/ml) in rabbits were 128.94 and 131.38; those of IL-12 (pg/ml) were 395.55 and 426.56, and then for those of IgG (μg/ml) were 14.70 and 13.94 respectively with the acetone extracts of Moringa oleifera and Vernonia amygdalina on D14 PT at 1000 mg/kg bw. This study indicates that Moringa oleifera and Vernonia amygdalina can be used as an alternative to synthetic anticoccidials. These plants could be used to increase the resistance of the immune system of rabbits to infestations of Eimeria species in rabbit farms.

Dendrobium Multi-Omics Reveal Lipid Remodeling in Response to Freezing

Freezing damage is a common phenomenon responsible for reduced yields of economic crops. Regulation of lipid metabolism plays an important role in plant growth and adaptation during freezing. We previously carried out transcriptome and untargeted metabolome analyses to determine the regulation of flavonol and anthocyanin biosynthesis during freezing treatment (FT) and post-freezing recovery (FR) in Dendrobium catenatum. However, changes in lipid levels are hard to confirm by untargeted metabolomics analysis alone. Regulation of lipid metabolism in response to freezing is largely unknown in Dendrobium. In this study, a multi-omics strategy was used to offer a better means of studying metabolic flow during FT and FR. To this end, 6976 proteins were identified by the 4D_label-free proteome, including 5343 quantified proteins. For each of the two conditions, we enriched differentially accumulated proteins (DAPs) into 15 gene ontology (GO) terms, including primary metabolism, lipid metabolism, and photosynthesis processes. We also identified 7 lipid categories and 3672 lipid species using lipidome assays. We found significant remodeling occurring in the phospholipid category during FT and FR. We also found that most sphingolipids were significantly upregulated. An integrated multi-omics analysis revealed significant changes in the expression levels of 141 mRNAs and encoding proteins under both FT and FR conditions. During FT, phospholipase A (PLA) and phospholipase D (PLD) were associated with phospholipid editing and galactolipid remodeling. These results provide valuable new insights into how the freezing tolerance of D. catenatum might be improved by genetic engineering.

Transcriptome Analysis and HPLC Profiling of Flavonoid Biosynthesis in Citrus aurantium L. during Its Key Developmental Stages

Citrus aurantium L. (sour orange) is a significant Chinese medicinal and fruit crop rich in flavonoids. However, the pathways and genes involved in flavonoid biosynthesis at the key developmental stages of Citrus aurantium L. are not fully understood. This study found that the total flavonoid concentration gradually decreased as the fruit developed. Additionally, it showed that neohesperidin was the main flavonoid in the early stages of sour orange fruit development. However, as the development stage progressed, naringin content increased rapidly and emerged as the main flavonoid component. From 27 cDNA libraries, RNA sequencing yielded 16.64 billion clean bases, including 8989 differentially expressed genes. We identified 74 flavonoid related unigenes mapped to the phenylalanine, tyrosine, and phenylpropanoid biosynthesis pathways. A total of 152 UDP-glucuronosyltransferase genes (UGTs) were identified from C. aurantium L. transcriptome database, in which 22 key flavonoid-correlated UGTs were divided into five main AtGT groups: E, G, I, L, M. We observed that the ethylene responsive factors (ERF) and myeloblastosis (MYB) family mainly regulated the key genes involved in flavonoid biosynthesis. Overall, our study generated extensive and detailed transcriptome data on the development of C. aurantium L. and characterized the flavonoid biosynthesis pattern during its fruit developmental stages. These results will benefit genetic modification or selection to increase the flavonoid content in sour oranges.

Characterization of the Key Bibenzyl Synthase in Dendrobium sinense

Dendrobium sinense, an endemic medicinal herb in Hainan Island, is rich in bibenzyls. However, the key rate-limited enzyme involved in bibenzyl biosynthesis has yet to be identified in D. sinense. In this study, to explore whether there is a significant difference between the D. sinense tissues, the total contents of bibenzyls were determined in roots, pseudobulbs, and leaves. The results indicated that roots had higher bibenzyl content than pseudobulbs and leaves. Subsequently, transcriptomic sequencings were conducted to excavate the genes encoding type III polyketide synthase (PKS). A total of six D. sinense PKS (DsPKS) genes were identified according to gene function annotation. Phylogenetic analysis classified the type III DsPKS genes into three groups. Importantly, the c93636.graph_c0 was clustered into bibenzyl synthase (BBS) group, named as D. sinense BBS (DsBBS). The expression analysis by FPKM and RT-qPCR indicated that DsBBS showed the highest expression levels in roots, displaying a positive correlation with bibenzyl contents in different tissues. Thus, the recombinant DsBBS-HisTag protein was constructed and expressed to study its catalytic activity. The molecular weight of the recombinant protein was verified to be approximately 45 kDa. Enzyme activity analysis indicated that the recombinant DsBBS-HisTag protein could use 4-coumaryol-CoA and malonyl-CoA as substrates for resveratrol production in vitro. The Vmax of the recombinant protein for the resveratrol production was 0.88 ± 0.07 pmol s⁻¹ mg⁻¹. These results improve our understanding with respect to the process of bibenzyl biosynthesis in D. sinense.

Metabolite analysis in Nymphaea 'Blue Bird' petals reveal the roles of flavonoids in color formation, stress amelioration, and bee orientation

In this study, metabolome of open petals (OP) and closed petals (CP) from Nymphaea 'Blue Bird' was firstly investigated. A total of 455 metabolites was identified in Nymphaea 'Blue Bird' petals, which was mainly composed of 100 flavonoids, 83 phenolic acids, 64 amino acids and derivatives, 60 lipids, 32 alkaloids, 32 organic acids, 24 nucleotides and derivatives, and 12 lignans and coumarins. By differential analysis, 192 metabolites were identified with variable importance in project ≥ 1, among which 83 and 109 metabolites were up- and down-regulated in OP group, respectively. Further analysis (Log₂ fold change ≥ 1) identified 26 and 7 metabolites exhibited significantly lower and higher contents in CP group, relative to OP group. Importantly, KEGG analysis indicated that flavonoid biosynthesis exhibited the most significant enrichment. qRT-PCR analysis indicated that the PAL, CHS, and HCDBR genes showed a significantly higher expression in OP group than in CP group. These data explain the increase of naringenin chalcone and phloretin in OP. However, there was no significant difference of total flavonoids between OP and CP groups. Considering the increase of H₂O₂ content and ultraviolet (UV) absorption peak in OP, our results implied that diurnal stressful conditions induced the degradation of flavonoids, which contributed to environmental stress amelioration. Moreover, a higher absorption peak of 360-380 nm UV was observed in the extract liquor of OP. The sensitivity maximum of the UV-photoreceptor of bees is situated around 340-380 nm UV. This suggested, as noted for the maximum absorption of dihydrokaempferol in 340-370 nm, rhythmic accumulation and loss of these differential flavonoids in Nymphaea 'Blue Bird' petals might enhance UV pattern to some degree, influencing pollinator attraction.