Metabolic and transcriptomic analyses reveal different metabolite biosynthesis profiles between leaf buds and mature leaves in Ziziphus jujuba mill

Gene editing, metabolomics, network pharmacology strategies to explore terpenoid content and anti-TMV activity in NtSPS1 knockout Nicotiana tabacum

The content of terpenoids in tobacco can alter its resistance to TMV. NtSPS1, a pivotal structural gene in tobacco, is capable to regulate the terpenoid content. In this study, we investigated the effect of NtSPS1 knockout in HD on the content of terpenoids and the anti-TMV activity of this mutant using gene editing, widely targeted metabolomics, network pharmacology, and molecular docking. 48 terpenoids (six up-regulated and five down-regulated) in NtSPS1 knockout tobacco compared with WT leaves. Notably, solanesol was remarkable downregulation which was lowered by fourfold and compounds 1 (log2FC = 18.2), 8 (log2FC = 16.7) were significant upregulation between the mutants and wild-type line leaves. The 46 terpenoid's target network encompassed 150 nodes, 509 edges and their underlying mechanisms in the therapeutic management of TMV are discussed. Furthermore, the network pharmacology and molecular docking revealed that compounds 16, 18, 23, 27, and 36 exhibited significant affinity in their respective interactions. Ultimately, five compounds were assayed for their anti-TMV effects, noteworthily, compounds 36 showed potential anti-TMV activity. Above all, we adopted a multifaceted approach to gain a comprehensive understanding of the terpenoid content and anti-TMV properties in NtSPS1 knockout HD. It enlightens the therapeutic potential of NtSPS1 knockout tobacco and it is helpful to find undescribed anti-TMV activity inhibitors, as well as searching for new anti-TMV candidates from the mutants.

Integrated Transcriptomic and Metabolomic Analysis Reveals Tissue-Specific Flavonoid Biosynthesis and MYB-Mediated Regulation of UGT71A1 in Panax quinquefolius

Panax quinquefolius is a globally valued medicinal plant rich in bioactive flavonoids, yet the molecular mechanisms underlying their biosynthesis remain poorly understood. In this study, we integrated transcriptomic and metabolomic analyses to investigate tissue-specific flavonoid accumulation and regulatory networks in roots, leaves, and flowers. Metabolomic profiling identified 141 flavonoid metabolites, with flavones, flavonols, and C-glycosylflavones predominantly enriched in aerial tissues (leaves and flowers), while specific glycosides like tricin 7-O-acetylglucoside showed root-specific accumulation. Transcriptome sequencing revealed 15,551-18,946 DEGs across tissues, and the reliability of the transcriptomic data was validated by qRT-PCR. KEGG and GO annotation analyses suggested that these DEGs may play a crucial role in the biosynthesis and metabolism of secondary metabolites. From the DEGs, UGTs and MYB TFs were identified and subjected to correlation analysis. Functional validation through in vitro enzymatic assays confirmed that PqUGT71A1 catalyzes apigenin and naringenin glycosylation at the 7-OH position. Additionally, subcellular localization and yeast one-hybrid assays demonstrated that PqMYB7 and PqMYB13 interact with the PqUGT71A1 promoter and activate its expression.. This study unveils the spatial dynamics of flavonoid metabolism in P. quinquefolius and establishes a MYB-UGT regulatory axis, providing critical insights for metabolic engineering and bioactive compound optimization in medicinal plants.

Genetic Diversity Analysis and GWAS of Plant Height and Ear Height in Maize Inbred Lines from South-East China

Maize is a critical crop for food, feed, and bioenergy worldwide. This study characterized the genetic diversity and population structure of 212 important inbred lines collected from the Southeast China breeding program using the Maize6H-60K single nucleotide polymorphism (SNP) array. To investigate the genetic architecture of plant height (PH) and ear height (EH), genome-wide association analysis (GWAS) was performed on this population in 2021 and 2022. Cluster analysis and population genetic structure analysis grouped the 212 maize inbred lines into 10 distinct categories. GWAS identified significant associations for PH, EH, and the EH/PH ratio. A total of 40 significant SNP (p < 8.55359 × 10-7) were detected, including nine associated with PH, with phenotypic variation explained (PVE) ranging from 3.42% to 25.92%. Additionally, 16 SNP were linked to EH, with PVE ranging from 2.49% to 38.49%, and 15 SNP were associated with the EH/PH ratio, showing PVE between 3.43% and 16.83%. Five stable SNP, identified across two or more environments, were further analyzed. Three of these SNP loci are reported for the first time in this study: two loci associated with the PH, AX-108020973, and AX-108022922, as well as one new locus, AX-108096437, which was significantly associated with the EH/PH ratio. Additionally, two other significant SNP (AX-247241325 and AX-108097244) were located within a 2 Mb range of previously identified QTL and/or related SNP. Within the 200 kb confidence intervals of these five stable SNP loci, 76 functionally annotated genes were identified. Further functional analysis indicated that 14 of these genes may play a role in regulating plant morphology, which is primarily involved in hormone synthesis, microtubule development, root growth, and cell division regulation. For instance, the homologous genes GRMZM2G375249 and GRMZM2G076029 in maize correspond to OsPEX1 in rice, a protein similar to extension proteins that are implicated in lignin biosynthesis, plant growth promotion, and the negative regulation of root growth through gibberellin-mediated pathways. The candidate gene corresponding to AX-108097244 is GRMZM2G464754; previous studies have reported its involvement in regulating EH in maize. These findings enhance the understanding of QTL associated with maize plant-type traits and provide a foundation for cloning PH, EH-related genes. Therefore, the results also support the development of functional markers for target genes and the breeding of improved maize varieties.

Unraveling site-specific seed formation abnormalities in Picea neoveitchii Mast. trees via widely metabolomic and transcriptomic analysis

Picea neoveitchii Mast. is a rare and threatened species of evergreen coniferous tree in China, commonly facing issues such as damaged seeds, abnormal seed growth, and empty seed shells. These abnormalities vary by location; unfortunately, the reasons behind these inconsistencies are completely unknown. This study compared seeds from two 150-year-old trees located in Taibai (Shaanxi province, TB150) and Zhouqu (Gansu province, ZQ150). The results showed significant differences in 43 metabolites and hormone levels, with higher levels of indole-3-acetic acid (IAA), methyl jasmonate (MeJA), and brassinosteroid (BR) in ZQ150, which were associated with more viable seeds. In contrast, TB150 exhibited more damaged seeds and empty seed shells due to higher abscisic acid (ABA) levels. Moreover, to further investigate these inconsistencies, we performed de-novo transcriptomic assembly and functional annotation of unigenes using high-throughput sequencing. A total of 2,355 differentially expressed unigenes were identified between TB150 and ZQ150, with 1,280 upregulated and 1,075 downregulated. Hormone signaling and sugar metabolism-related unigenes were further examined for their role in seed development. ZQ150 increased the number of normal seeds by enhancing endogenous IAA levels and upregulating auxin signaling and sugar metabolism-related genes. Conversely, TB150 showed more empty seed shells, correlated with elevated ABA levels and the activation of ABA signaling genes. We hypothesize that enhanced IAA levels and the upregulation of sugar metabolism and auxin signaling genes promote normal seed development.

Research Progress and Biosynthetic Mechanisms of Nutritional Compounds Obtained from Various Organs During the Developmental Stages of a Medicinal Plant (Chinese Jujube)

The fruit of the jujube tree is high in nutrients and has various health benefits. China is a major producer of jujube, and it is now cultivated all around the world. Numerous studies have demonstrated the nutritional value and potential health advantages of bioactive compounds found in the jujube tree. Furthermore, the jujube tree has a remarkable 7000-year agricultural history. The jujube plant has developed a rich gene pool, making it a valuable resource for germplasm. Different studies have focused on the developmental stages of jujube fruits to identify the optimal time for harvest and to assess the changes in their bioactive natural compounds or products during the process of development but the molecular mechanism underlying the production of bioactive natural products in Z. jujuba is still poorly understood. Moreover, the potential differential expressed genes (DEGs) identified as responsible for the synthesis of these compounds should be further functionally verified. It has been noticed that the contents of total flavonoids, total phenolic, and vitamin C increase significantly during the ripening process, while the contents of soluble sugars and organic acids decrease gradually. In this review, we have also scrutinized the challenges that hinder the utilization of jujube fruit resources and suggested potential areas for further research. As such, our review serves as a valuable resource for the future development of jujube-based nutritional compounds and the incorporation of their nutritional elements into the functional foods industry.

Integrated metabolomics and transcriptomics analysis provides insights into biosynthesis and accumulation of flavonoids and glucosinolates in different radish varieties

Radish is an important vegetable worldwide, with wide medicinal functions and health benefits. The quality of radish, strongly affected by phytochemicals like flavonoids and glucosinolates, are quite different depending on the radish varieties. However, the comprehensive accumulation profiles of secondary metabolites and their molecular regulatory mechanisms in different radish cultivars remain unclear thus far. Herein, we comprehensively analyzed the secondary metabolite and gene expression profiles of the flesh and skin of four popular radish varieties with different flesh and/or skin colors, using UPLC-MS/MS-based metabolomics and transcriptomics approach combined with RT-qPCR. The results showed that altogether 352 secondary metabolites were identified in radish, of which flavonoids and phenolic acids accounted for 60.51% of the total. The flesh and skin of each variety exhibited distinct metabolic profiles, making them unique in coloration, flavor, taste, and nutritional quality. The differential metabolites were mostly enriched in flavonoid biosynthesis, flavone and flavonol biosynthesis, phenylpropanoid biosynthesis, and glucosinolate biosynthesis pathway. Further, 19 key genes regulating the differential accumulation of flavonoids among different radish varieties were identified, such as RsCHS, RsCCOAMT, RsF3H, RsFLS, RsCYP75B1, RsDFR, and RsANS that were significantly upregulated in red-colored radish tissue. Also, 10 key genes affecting the differential accumulation of glucosinolates among different varieties were identified, such as RsCYP83B1, RsSUR1, and RsST5a that were significantly increased in the skin of green radish. Moreover, systematical biosynthetic pathways of flavonoids and glucosinolates and co-expression networks between genes and metabolites were constructed based on integrative analysis between metabolomics and transcriptomics. Our findings provide a novel insight into the mechanisms of radish quality formation, thereby providing a molecular basis for breeding and cultivation of radish with excellent nutritional quality.

Identification of Key Enzymes and Genes Modulating L-Ascorbic Acid Metabolism During Fruit Development of Lycium chinense by Integrating Metabolome, Transcriptome, and Physiological Analysis

Lycium chinense is acknowledged for its substantial nutritional benefits, particularly attributed to the high levels of ascorbic acid (AsA) found in its fruits. The "Mengqi No.1" variety of L. chinense, which is cultivated in Qinghai, is known for its high yield and exceptional quality. We utilized the "Mengqi No.1" variety as experimental materials and combined metabolomic, transcriptomic, and physiological analyses to investigate the metabolites, genes, and enzymes related to AsA metabolism in L. chinense fruits. The results revealed nine differential metabolites associated with AsA metabolism in L. chinense fruits across three stages, including 1D-Myo-Inositol-1,4-Bisphosphate, D-Fructose, L-(+)-Arabinose, I-Inositol, L-Arabinitol, D-Galactose-1-P, lactose, α-D-Glucose, and D-Glucose-6-P. Notably, the contents of D-Glucose-6-P, D-Galactose-1-P, and D-Fructose were increased as the fruit developed. Additionally, fresh weight, longitudinal length, and radial width were increased, while the contents of AsA and DHA were decreased. GalDH and DHAR are critical enzymes for the accumulation of AsA and DHA, exhibiting positive correlation coefficient. Furthermore, PMM1, PMM5, GME2, and GME3 were identified as key regulatory genes in the L-Galactose pathway of AsA synthesis, influencing D-Galactose-1-P, D-Glucose-6-P, α-D-Glucose, and D-Fructose. DHAR1 and DHAR2 are considered key positive regulator genes of AsA and DHA in the AsA-GSH cycle. However, the majority of genes (nine) act as negative regulators of AsA and DHA. These findings provide a foundation for the understanding of the regulatory mechanism of AsA metabolism in L. chinense fruits and offer insights into the utilization of AsA from L. chinense.

Integrated Metabolomics and Transcriptomics Analyses Reveal the Regulatory Mechanisms of Anthocyanin and Carotenoid Accumulation in the Peel of Coffea arabica

The color of coffee fruits is influenced by several factors, including cultivar, ripening stage, and metabolite composition. However, the metabolic accumulation of pigments and the molecular mechanisms underlying peel coloration during the ripening process of Coffea arabica L. remain relatively understudied. In this study, UPLC-MS/MS-based metabolomics and RNA sequencing (RNA-seq)-based transcriptomics were integrated to investigate the accumulation of anthocyanins and carotenoids in the peel of Coffea arabica at different ripening stages: green peel (GP), green-yellow peel (GYRP), red peel (RP), and red-purple peel (RPP). This integration aimed at elucidating the molecular mechanisms associated with these changes. A total of ten anthocyanins, six carotenoids, and thirty-five xanthophylls were identified throughout the ripening process. The results demonstrated a gradual decrease in the total carotenoid content in the peel with fruit maturation, while anthocyanin content increased significantly. Notably, the accumulation of specific anthocyanins was closely associated with the transition of peel colors from green to red. Integrated metabolomics and transcriptomics analyses identified the GYRP stage as critical for this color transition. A weighted gene co-expression network analysis (WGCNA) revealed that enzyme-coding genes such as 3AT, BZ1, and lcyE, along with transcription factors including MYB, NAC, and bHLH, which interact with PHD and SET TR, may regulate the biosynthesis of anthocyanins and carotenoids, thereby influencing peel pigmentation. These findings provide valuable insights into the molecular mechanisms underlying the accumulation of anthocyanins and carotenoids in Coffea arabica peel during fruit maturation.

Transcriptomics and metabolomics revealed the molecular basis of the color formation in the roots of Panax notoginseng

Panax notoginseng is a traditional Chinese medicine rich in many pharmacological components. The root of the 'Miaoxiang Sanqi No. 2' is yellow or greenish yellow, while a novel cultivar-'Wenyuan Ziqi No. 1' shows purple root and is thought to have high medicinal value. Little information is available about the anthocyanin biosynthesis in P. notoginseng root. In this study, we compared the 'Miaoxiang Sanqi No. 2' and 'Wenyuan Ziqi No. 1' in morphological, transcriptional and metabolic levels. The results showed that purple rich in the periderm, rhizome and phloem around cambium of the 'Wenyuan Ziqi No. 1' root and cyanidin 3-O-galactoside was the main anthocyanin causing purple. Moreover, 'Wenyuan Ziqi No. 1' highly accumulated in 155 metabolites, including flavones, phenylpropanoids and lipids. Transcriptome data showed that phenylpropanoid biosynthesis pathway genes are highly expressed in 'Wenyuan Ziqi No. 1'. Conjoint analysis showed that anthocyanin biosynthesis pathway substances were highly accumulated in 'Wenyuan Ziqi No. 1', and the expression level of structural genes involved in anthocyanin biosynthesis pathway was higher in 'Wenyuan Ziqi No. 1'. Meanwhile, eight R2R3-MYB genes that might be involved in anthocyanin biosynthesis were identified. The comprehensive analysis of two cultivars provides new insights into the understanding of root coloration in P. notoginseng.

A Novel Plant-Derived Biopesticide Mitigates Fusarium Root Rot of Angelica sinensis by Modulating the Rhizosphere Microbiome and Root Metabolome

Fusarium root rot caused by the Fusarium species complex significantly affects the yield and quality of Angelica sinensis, a valuable medicinal herb. Traditional management primarily relies on chemical fungicides, which have led to pathogen resistance, environmental hazards, and concerns regarding public health and the active components in A. sinensis. This study explores the efficacy of a novel plant-derived biopesticide Shi Chuang Zhi Feng Ning (T1; SCZFN), alongside Bacillus subtilis wettable powder (T2) and a chemical fungicide (T3), in controlling root rot and understanding their impacts on the rhizosphere microbial community and root metabolome. Results of the field experiment demonstrated that treatments T1 and T3 achieved control efficiencies of 73.17% and 75.45%, respectively, significantly outperforming T2 (39.99%) and the control. High-throughput sequencing revealed that all treatments altered the diversity and structure of microbial communities, with T1 and T2 reducing the abundance of taxa linked to root rot, such as Muribaculaceae spp., Humicola spp., Fusarium spp., and Mycochlamys spp. Treatment T1 notably enhanced beneficial bacterial taxa, including Acidobacteria spp., Nitrospira spp., and Pedosphaeraceae spp., involved in carbon cycling and plant growth promotion. Metabolomic analysis identified 39, 105, and 45 differentially expressed metabolites (DEMs) across the treatments, demonstrating T1's potential to modulate the root metabolome effectively. Further, a correlation analysis demonstrated a stronger correlation between distinct microorganisms with significant influence and DEMs of T1 treatment compared to other treatments. These findings underscore biopesticide SCZFN's role in enhancing plant health and disease suppression in A. sinensis, providing insights into its biocontrol mechanisms and supporting the development of sustainable disease management strategies in its cultivation.

Integrative analysis of transcriptome and metabolome reveal the differential tolerance mechanisms to low and high salinity in the roots of facultative halophyte Avicennia marina

Mangroves perform a crucial ecological role along the tropical and subtropical coastal intertidal zone where salinity fluctuation occurs frequently. However, the differential responses of mangrove plant at the combined transcriptome and metabolome level to variable salinity are not well documented. In this study, we used Avicennia marina (Forssk.) Vierh., a pioneer species of mangrove wetlands and one of the most salt-tolerant mangroves, to investigate the differential salt tolerance mechanisms under low and high salinity using inductively coupled plasma-mass spectrometry, transcriptomic and metabolomic analysis. The results showed that HAK8 was up-regulated and transported K+ into the roots under low salinity. However, under high salinity, AKT1 and NHX2 were strongly induced, which indicated the transport of K+ and Na+ compartmentalization to maintain ion homeostasis. In addition, A. marina tolerates low salinity by up-regulating ABA signaling pathway and accumulating more mannitol, unsaturated fatty acids, amino acids' and L-ascorbic acid in the roots. Under high salinity, A. marina undergoes a more drastic metabolic network rearrangement in the roots, such as more L-ascorbic acid and oxiglutatione were up-regulated, while carbohydrates, lipids and amino acids were down-regulated in the roots, and, finally, glycolysis and TCA cycle were promoted to provide more energy to improve salt tolerance. Our findings suggest that the major salt tolerance traits in A. marina can be attributed to complex regulatory and signaling mechanisms, and show significant differences between low and high salinity.

The Effect of Heat Stress on Sensory Properties of Fresh Oysters: A Comprehensive Study Using E-Nose, E-Tongue, Sensory Evaluation, HS-SPME-GC-MS, LC-MS, and Transcriptomics

Heat stress has received growing concerns regarding the impact on seafood quality. However, the effects of heat stress on the sensory properties of seafood remain unknown. In this study, the sensory properties of fresh oyster (Crassostrea ariakensis) treated with chronic heat stress (30 °C) for 8 weeks were characterized using electronic nose, electronic tongue, sensory evaluation, HS-SPME-GC-MS, LC-MS and transcriptomics. Overall, chronic heat stress reduced the overall sensory properties of oysters. The metabolic network constructed. based on enrichment results of 423 differential metabolites and 166 differentially expressed genes, showed that the negative effects of chronic heat stress on the sensory properties of oysters were related to oxidative stress, protein degradation, lipid oxidation, and nucleotide metabolism. The results of the study provide valuable insights into the effects of heat stress on the sensory properties of oysters, which are important for ensuring a sustainable supply of high-quality seafood and maintaining food safety.

Integrated physiological, transcriptomics and metabolomics analysis revealed the molecular mechanism of Bupleurum chinense seedlings to drought stress

Drought stress is a prominent abiotic factor that adversely influences the growth and development of Bupleurum chinense during its seedling stage, negatively impacting biomass and secondary metabolite production, thus affecting yield and quality. To investigate the molecular mechanism underlying the response of B. chinense seedlings under drought stress, this study employed comprehensive physiological, transcriptomic, and metabolomic analyses. The results revealed that under drought stress, the root soluble sugar and free proline content in B. chinense seedlings significantly increased, while the activities of SOD, POD, and CAT increased in the leaves. These findings indicate the presence of distinct response mechanisms in B. chinense to cope with drought stress. Integrated analysis further identified significant correlations between genes and metabolites related to amino acid biosynthesis in the leaves, as well as genes and metabolites associated with acetaldehyde and dicarboxylic acid metabolism. In the roots, genes and metabolites related to plant hormone signaling and the tricarboxylic acid (TCA) cycle showed significant correlations. These findings provide vital views into the molecular-level response mechanisms of B. chinense under drought stress. Moreover, this study establishes the groundwork for identifying drought-tolerant genes and breeding drought-resistant varieties, which could improve the drought tolerance of medicinal plants and have broader implications for agriculture and crop production in water-scarce areas.

The Flavor Characteristics and Metabolites of Three Commercial Dried Jujube Cultivars

To understand the flavor and metabolite differences between the three commercial dried jujube cultivars Huizao (HZ), Hamazao 1 (HMZ), and Qiyuexian (QYX), their soluble sugars, organic acids, volatiles, and metabolites were systematically investigated. The results show that sucrose and malic acid were the main soluble sugar and organic acids contained in these dried jujubes, respectively. Sucrose (573.89 mg/g DW) had the highest presence in HZ, and the total sugar content (898.33 mg/g DW) was the highest in QYX. Both of these had a low total acid content, resulting in relatively high sugar-acid ratios (105.49 and 127.86, respectively) compared to that of HMZ (51.50). Additionally, 66 volatile components were detected in the 3 jujubes. These mainly included acids, aldehydes, esters, and ketones (90.5-96.49%). Among them, (E)-2-nonenal, (E)-2-decenal, heptanal, decanal, nonanal, and octanal were identified as the key aromatic substances of the dried jujubes, and their contents were the highest in HMZ. Moreover, 454 metabolites were identified, including alkaloids, amino acids, flavonoids, lipids, nucleotides, and terpenoids. The highest contents of flavonoids (5.6%) and lipids (24.9%) were detected in HMZ, the highest contents of nucleotides (10.2%) and alkaloids (27%) were found in QYX, and the contents of saccharides (5.7%) and amino acids (23.6%) were high in HZ. Overall, HZ, HMZ, and QYX significantly differ in their flavor and nutrition. HZ tastes better, HMZ is more fragrant, and QYX and HMZ possess higher nutritional values.

Carotenoid Accumulation in the Rhododendron chrysanthum Is Mediated by Abscisic Acid Production Driven by UV-B Stress

Rhododendron chrysanthum (R. chrysanthum) development is hampered by UV-B sunlight because it damages the photosynthetic system and encourages the buildup of carotenoids. Nevertheless, it is still unclear how R. chrysanthum repairs the photosynthetic system to encourage the formation of carotenoid pigments. The carotenoid and abscisic acid (ABA) concentrations of the R. chrysanthum were ascertained in this investigation. Following UV-B stress, the level of carotenoids was markedly increased, and there was a strong correlation between carotenoids and ABA. The modifications of R. chrysanthum's OJIP transient curves were examined in order to verify the regulatory effect of ABA on carotenoid accumulation. It was discovered that external application of ABA lessened the degree of damage on the donor side and lessened the damage caused by UV-B stress on R. chrysanthum. Additionally, integrated metabolomics and transcriptomics were used to examine the changes in differentially expressed genes (DEGs) and differential metabolites (DMs) in R. chrysanthum in order to have a better understanding of the role that ABA plays in carotenoid accumulation. The findings indicated that the majority of DEGs were connected to carotenoid accumulation and ABA signaling sensing. To sum up, we proposed a method for R. chrysanthum carotenoid accumulation. UV-B stress activates ABA production, which then interacts with transcription factors to limit photosynthesis and accumulate carotenoids, such as MYB-enhanced carotenoid biosynthesis. This study showed that R. chrysanthum's damage from UV-B exposure was lessened by carotenoid accumulation, and it also offered helpful suggestions for raising the carotenoid content of plants.

Integrative physiological, metabolomic, and transcriptomic analysis reveals the drought responses of two apple rootstock cultivars

BACKGROUND

Drought is considered the main environmental factor restricting apple production and thus the development of the apple industry. Rootstocks play an important role in enhancing the drought tolerance of apple plants. Studies of the physiology have demonstrated that 'ZC9-3' is a strong drought-resistant rootstock, whereas 'Jizhen-2' is a weak drought-resistant rootstock. However, the metabolites in these two apple rootstock varieties that respond to drought stress have not yet been characterized, and the molecular mechanisms underlying their responses to drought stress remain unclear.

RESULTS

In this study, the physiological and molecular mechanisms underlying differences in the drought resistance of 'Jizhen-2' (drought-sensitive) and 'ZC9-3' (drought-resistant) apple rootstocks were explored. Under drought stress, the relative water content of the leaves was maintained at higher levels in 'ZC9-3' than in 'Jizhen-2', and the photosynthetic, antioxidant, and osmoregulatory capacities of 'ZC9-3' were stronger than those of 'Jizhen-2'. Metabolome analysis revealed a total of 95 and 156 differentially accumulated metabolites in 'Jizhen-2' and 'ZC9-3' under drought stress, respectively. The up-regulated metabolites in the two cultivars were mainly amino acids and derivatives. Transcriptome analysis revealed that there were more differentially expressed genes and transcription factors in 'ZC9-3' than in 'Jizhen-2' throughout the drought treatment. Metabolomic and transcriptomic analysis revealed that amino acid biosynthesis pathways play key roles in mediating drought resistance in apple rootstocks. A total of 13 metabolites, including L-α-aminoadipate, L-homoserine, L-threonine, L-isoleucine, L-valine, L-leucine, (2S)-2-isopropylmalate, anthranilate, L-tryptophan, L-phenylalanine, L-tyrosine, L-glutamate, and L-proline, play an important role in the difference in drought resistance between 'ZC9-3' and 'Jizhen-2'. In addition, 13 genes encoding O-acetylserine-(thiol)-lyase, S-adenosylmethionine synthetase, ketol-acid isomeroreductase, dihydroxyacid dehydratase, isopropylmalate isomerase, branched-chain aminotransferase, pyruvate kinase, 3-dehydroquinate dehydratase/shikimate 5-dehydrogenase, N-acetylglutamate-5-P-reductase, and pyrroline-5-carboxylate synthetase positively regulate the response of 'ZC9-3' to drought stress.

CONCLUSIONS

This study enhances our understanding of the response of apple rootstocks to drought stress at the physiological, metabolic, and transcriptional levels and provides key insights that will aid the cultivation of drought-resistant apple rootstock cultivars. Especially, it identifies key metabolites and genes underlying the drought resistance of apple rootstocks.

Metabolomics reveals the importance of metabolites in Mussaenda pubescens for antioxidant properties and quality traits

Mussaenda pubescens (Mp) is a valuable medicinal plant that has traditionally been used for medicinal purposes or as a tea substitute. However, there are few studies on the comprehensive and dynamic evaluation of Mp metabolites. This study used an ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS) approach and biochemical analysis to investigate substance changes in leaves at three different stages and elucidate the relationship between metabolites and antioxidant capacity. The findings showed that Mp leaves contained 957 metabolites, the majority of which were phenolic acids, lipids, and terpenoids. The metabolite profiling of Mp leaves was significantly influenced by their growth and development at different stages. A total of 317 differentially accumulated metabolites (DAMs) were screened, including 150 primary metabolites and 167 secondary metabolites, with 202 DAMs found in bud leaf vs. tender leaf, 54 DAMs in tender leaf vs. mature leaf, and 254 DAMs in bud leaf vs. mature leaf. Total phenolics, flavonoids, and anthocyanin concentrations decreased as Mp leaves grew and developed, whereas terpenoids increased significantly. The secondary metabolites also demonstrated a positive correlation with antioxidant activity. Phenolics, flavonoids, terpenoids, and anthocyanins were the primary factors influencing the antioxidant activity of leaves. These findings provide new insights into the metabolite formation mechanism, as well as the development and utilization of Mp tea.

Rice husk ash based growing media impact on cucumber and melon growth and quality

Rice husk, an agricultural waste from the rice industry, can cause serious environmental pollution if not properly managed. However, rice husk ash (RHA) has been found to have many positive properties, making it a potential replacement for non-renewable peat in soilless planting. Thus, this study investigated the impact of a RHA composite substrate on the growth, photosynthetic parameters, and fruit quality of cucumber (Yuyi longxiang variety) and melon (Yutian yangjiaomi variety). The RHA, peat, vermiculite, and perlite were blended in varying proportions, with the conventional seedling substrate (peat:vermiculite:perlite = 1:1:1 volume ratio) serving as the control (CK). All plants were cultivated in barrels filled with 10L of the mixed substrates. The results from this study found that RHA 40 (RHA:peat:vermiculite:perlite = 4:4:1:1 volume ratio) significantly enhanced substrate ventilation and positively influenced the stem diameter, root activity, seedling index, chlorophyll content, net photosynthetic rate (Pn), stomatal conductance (Gs), and transpiration rate (Tr) of cucumber and melon plants. Additionally, plant planted using RHA 40, the individual fruit weight of cucumber and melon found to increase by 34.62% and 21.67%, respectively, as compared to the control. Aside from that, both cucumber and melon fruits had significantly higher sucrose, total soluble sugar, vitamin C, and soluble protein levels. This subsequently improved the activity of sucrose synthase and sucrose phosphate synthase in both cucumber and melon. In conclusion, the RHA 40 found to best promote cucumber and melon plant growth, increase plant leaf photosynthesis, and improve cucumber and melon fruit quality, making it a suitable substrate formula for cucumber and melon cultivation in place of peat.

Integrated metabolomic and transcriptomic dynamic profiles of endopleura coloration during fruit maturation in three walnut cultivars

BACKGROUND

The color of endopleura is a vital factor in determining the economic value and aesthetics appeal of nut. Walnuts (Juglans) are a key source of edible nuts, high in proteins, amino acids, lipids, carbohydrates. Walnut had a variety endopleura color as yellow, red, and purple. However, the regulation of walnut endopleura color remains little known.

RESULTS

To understand the process of coloration in endopleura, we performed the integrative analysis of transcriptomes and metabolomes at two developmental stages of walnut endopleura. We obtained total of 4,950 differentially expressed genes (DEGs) and 794 metabolites from walnut endopleura, which are involved in flavonoid and phenolic biosynthesis pathways. The enrichment analysis revealed that the cinnamic acid, coniferyl alcohol, naringenin, and naringenin-7-O-glucoside were important metabolites in the development process of walnut endopleura. Transcriptome and metabolome analyses revealed that the DEGs and differentially regulated metabolites (DRMs) were significantly enriched in flavonoid biosynthesis and phenolic metabolic pathways. Through co-expression analysis, CHS (chalcone synthase), CHI (chalcone isomerase), CCR (cinnamoyl CoA reductase), CAD (cinnamyl alcohol dehydrogenase), COMT (catechol-Omethyl transferase), and 4CL (4-coumaroyl: CoA-ligase) may be the key genes that potentially regulate walnut endopleura color in flavonoid biosynthesis and phenolic metabolic pathways.

CONCLUSIONS

This study illuminates the metabolic pathways and candidate genes that underlie the endopleura coloration in walnuts, lay the foundation for further study and provides insights into controlling nut's colour.

Transcriptomic and metabolic analysis unveils the mechanism behind leaf color development in Disanthus cercidifolius var. longipes

Introduction: Leaf coloration in Disanthus cercidifolius var. longipes results from the interplay of various pigments undergoing complex catalytic reactions. Methods: We aimed to elucidate the mechanisms of pigment biosynthesis affecting leaf color transition in D. cercidifolius var. longipes by analyzing variations in pigment accumulation and levels of gene expression. Results: We identified 468, 577, and 215 differential metabolites in green leaves (GL), gradual-color-changing leaves (GCCL), and red leaves (RL), respectively, with 94 metabolites shared across all comparisons. Metabolite accumulation patterns were similar among GL, GCCL, and RL, with flavonoids being the main differential metabolites. Delphinidin, malvidin, and petunidin derivatives were mostly accumulated in GCCL, whereas cyanidin, pelargonidin, and peonidin derivatives accumulated in RL. Transcriptome sequencing was used to identify differentially expressed genes. The expression of anthocyanin biosynthetic pathway genes was associated with anthocyanin accumulation patterns. Discussion: Our findings reveal that the content of delphinidin, malvidin, petunidin, and carotenoids collectively determines the gradual transition of leaf color from green in spring and summer to green, purple, and orange-yellow in early autumn, whereas the content of cyanidin, peonidin, pelargonidin, and carotenoids together causes the autumnal transition to red or orange-red colors as leaves of D. cercidifolius var. longipes age.

Red light alleviates Cd toxicity in Egeria densa by modifying carbon-nitrogen metabolism and boosting energy metabolism

Among the various pollutants detected in aquatic ecosystems, cadmium (Cd) is considered as one of the most hazardous. Freshwater macrophytes have been recognized as possible candidates for eliminating Cd from environment. Nevertheless, the impact of light quality on their ability to tolerate Cd toxicity remains unclear, and the underlying mechanisms have yet to be fully elucidated. In this study, we utilized physiological testing and metabolomics to explore the potential mechanisms by which light quality influences the ability of Egeria densa, a significant Cd hyperaccumulator, to withstand Cd toxicity. The study demonstrated that following Cd treatment, E. densa grown under red light exhibited superior photosynthetic efficiency compared to those grown under blue light, as evidenced by significantly increased photosynthetic rate, higher starch content, and greater activity of photosynthetic enzymes. Moreover, metabolomic analyses revealed that under Cd stress, E. densa grown under red light exhibited an enhanced glycolysis for increased energy production. Sucrose metabolism was also improved to generate sufficient sugar including glucose, fructose and mannose for osmotic adjustment. Moreover, under red light, the heightened production of α-ketoglutarate via tricarboxylic acid (TCA) cycle redirected nitrogen flow towards the synthesis of resilient substances such as γ-Aminobutyric Acid (GABA) and methionine. The production of these substances was ∼2.0 and 1.3 times greater than that of treatment with Cd under blue light, thereby improving E. densa's capacity to withstand Cd stress. This study represents the initial investigation into the possible mechanisms by which light quality influences the ability of E. densa to withstand Cd toxicity through regulating CN metabolism. Furthermore, these findings have the potential to improve phytoremediation strategies aimed at reducing Cd pollution.

Transcriptome Characterization of Pigment-Related Genes in Jujube (Ziziphus Jujuba Mill.) Peel at Different Growth Stages

One of the most important qualities of jujube fruit is its color. Chlorophyll, carotenoid, and anthocyanin all play important roles in the coloring of jujube fruit. However, few studies have focused on the pigment molecular mechanism. In the present study, jujube peels of 'Sanbianhong' in three growth stages were evaluated for their gene expression characteristics and gene regulation related to pigment formation using the transcriptome sequencing analysis. A total of 84.86 Gb of clean data were obtained in the analysis. In the FS1 vs. FS3, FS1 vs. FS5, and FS3 vs. FS5, 4,530, 11,012, and 9,072 differentially expressed genes (DEGs) were identified, respectively. The inter-group screening among the three comparisons yielded 1430 common DEGs. Among these DEGs, 27, 16, and 28 genes were enriched in chlorophyll, carotenoid, and anthocyanin metabolic pathways, respectively. Twelve genes were chosen at random, and the accuracy of the transcriptome data were confirmed using qRT-PCR. The molecular mechanism underlying the pigmentation of jujube fruit was elucidated at the transcriptome level, which would provide a scientific basis for the subsequent functional studies on the color-regulating genes of jujube fruits.

Phenolic Acids and Flavonoids Play Important Roles in Flower Bud Differentiation in Mikania micrantha: Transcriptomics and Metabolomics

Mikania micrantha is a highly invasive vine, and its ability to sexually reproduce is a major obstacle to its eradication. The long-distance dissemination of M. micrantha depends on the distribution of seeds; therefore, inhibiting M. micrantha flowering and seed production is an effective control strategy. The number of blooms of M. micrantha differs at different altitudes (200, 900, and 1300 m). In this study, we used a combination of metabolomics and transcriptomics methods to study the patterns of metabolite accumulation in the flower buds of M. micrantha. Using LC-MS/MS, 658 metabolites were found in the flower buds of M. micrantha at three different altitudes (200, 900, and 1300 m). Flavonoids and phenolic acids were found to be the main differential metabolites, and their concentrations were lower at 900 m than at 200 m and 1300 m, with the concentrations of benzoic acid, ferulic acid, and caffeic acid being the lowest. The biosynthesis pathways for flavonoids and phenolic compounds were significantly enriched for differentially expressed genes (DEGs), according to the results of transcriptome analysis. The production of flavonoid and phenolic acids was strongly linked with the expressions of phenylalanine ammonia-lyase (PAL), caffeoyl-CoA O-methyltransferase (COMT), and 4-coumarate-CoA ligase (4CL), according to the results of the combined transcriptome and metabolome analysis. These genes' roles in the regulation of distinct phenolic acids and flavonoids during M. micrantha bud differentiation are still unknown. This study adds to our understanding of how phenolic acids and flavonoids are regulated in M. micrantha flower buds at various altitudes and identifies regulatory networks that may be involved in this phenomenon, offering a new approach for the prevention and management of M. micrantha.

Characterization of color variation in bamboo sheath of Chimonobambusa hejiangensis by UPLC-ESI-MS/MS and RNA sequencing

BACKGROUND

Chimonobambusa hejiangensis (C.hejiangensis) is a high-quality bamboo species native to China, known for its shoots that are a popular nutritional food. Three C.hejiangensis cultivars exhibit unique color variation in their shoot sheaths, however, the molecular mechanism behind this color change remains unclear.

METHODS

We investigated flavonoid accumulation in the three bamboo cultivar sheaths using metabolomics and transcriptomics.

RESULTS

UPLC-MS/MS identified 969 metabolites, with 187, 103, and 132 having differential accumulation in the yellow-sheath (YShe) vs. spot-sheath (SShe)/black-sheath (BShe) and SShe vs. BShe comparison groups. Flavonoids were the major metabolites that determined bamboo sheath color through differential accumulation of metabolites (DAMs) analysis. Additionally, there were 33 significantly differentially expressed flavonoid structural genes involved in the anthocyanin synthesis pathway based on transcriptome data. We conducted a KEGG analysis on DEGs and DAMs, revealing significant enrichment of phenylpropanoid and flavonoid biosynthetic pathways. Using gene co-expression network analysis, we identified nine structural genes and 29 transcription factors strongly linked to anthocyanin biosynthesis.

CONCLUSION

We identified a comprehensive regulatory network for flavonoid biosynthesis which should improve our comprehension of the molecular mechanisms responsible for color variation and flavonoid biosynthesis in bamboo sheaths.

Integrative analysis of metabolome and transcriptome reveals the different metabolite biosynthesis profiles related to palatability in winter and spring shoot in moso bamboo

Moso bamboo winter shoot has good taste and rich nutritional value. To reveal the causes and regulatory mechanism of palatability deterioration from winter to spring shoot, a conjoint analysis of metabolome and transcriptome was conducted on winter and spring shoots of moso bamboo. Totally 909 metabolites were characterized for the first time. The significant increase of hydrolyzed tannin content intensified the bitterness of spring shoot, which was positively regulated by key metabolite (gallic acid) and genes (DAHPS, DHQS, DHQ, SDH) in the biosynthesis pathway of hydrolyzed tannin. The accumulation of lignified components enhanced the roughness of spring shoot, which was closely connected with the significant changes of important metabolites (cinnamic acid, ferulic acid, UDP-glucose and UDP-xylose) and up-regulation of most enzyme genes involved in the biosynthesis pathways of lignin, cellulose and hemicellulose. The present study provides theoretical support for understanding palatability transition and directional improvement of edible quality of moso bamboo shoots.

Metabolomic and transcriptomice analyses of flavonoid biosynthesis in apricot fruits

Introduction

Flavonoids, as secondary metabolites in plants, play important roles in many biological processes and responses to environmental factors.

Methods

Apricot fruits are rich in flavonoid compounds, and in this study, we performed a combined metabolomic and transcriptomic analysis of orange flesh (JN) and white flesh (ZS) apricot fruits.

Results and discussion

A total of 222 differentially accumulated flavonoids (DAFs) and 15855 differentially expressed genes (DEGs) involved in flavonoid biosynthesis were identified. The biosynthesis of flavonoids in apricot fruit may be regulated by 17 enzyme-encoding genes, namely PAL (2), 4CL (9), C4H (1), HCT (15), C3'H (4), CHS (2), CHI (3), F3H (1), F3'H (CYP75B1) (2), F3'5'H (4), DFR (4), LAR (1), FLS (3), ANS (9), ANR (2), UGT79B1 (6) and CYP81E (2). A structural gene-transcription factor (TF) correlation analysis yielded 3 TFs (2 bHLH, 1 MYB) highly correlated with 2 structural genes. In addition, we obtained 26 candidate genes involved in the biosynthesis of 8 differentially accumulated flavonoids metabolites in ZS by weighted gene coexpression network analysis. The candidate genes and transcription factors identified in this study will provide a highly valuable molecular basis for the in-depth study of flavonoid biosynthesis in apricot fruits.

Integration of Phosphoproteomics and Transcriptome Studies Reveals ABA Signaling Pathways Regulate UV-B Tolerance in Rhododendron chrysanthum Leaves

The influence of UV-B stress on the growth, development, and metabolism of alpine plants, such as the damage to DNA macromolecules, the decline in photosynthetic rate, and changes in growth, development, and morphology cannot be ignored. As an endogenous signal molecule, ABA demonstrates a wide range of responses to UV-B radiation, low temperature, drought, and other stresses. The typical effect of ABA on leaves is to reduce the loss of transpiration by closing the stomata, which helps plants resist abiotic and biological stress. The Changbai Mountains have a harsh environment, with low temperatures and thin air, so Rhododendron chrysanthum (R. chrysanthum) seedlings growing in the Changbai Mountains can be an important research object. In this study, a combination of physiological, phosphorylated proteomic, and transcriptomic approaches was used to investigate the molecular mechanisms by which abiotic stress leads to the phosphorylation of proteins in the ABA signaling pathway, and thereby mitigates UV-B radiation to R. chrysanthum. The experimental results show that a total of 12,289 differentially expressed genes and 109 differentially phosphorylated proteins were detected after UV-B stress in R. chrysanthum, mainly concentrated in plant hormone signaling pathways. Plants were treated with ABA prior to exposure to UV-B stress, and the results showed that ABA mitigated stomatal changes in plants, thus confirming the key role of endogenous ABA in plant adaptation to UV-B. We present a model that suggests a multifaceted R. chrysanthum response to UV-B stress, providing a theoretical basis for further elaboration of the mechanism of ABA signal transduction regulating stomata to resist UV-B radiation.

MutantPpbbx24-delgene positively regulates light-induced anthocyanin accumulation in the red pear.. [DOI: 10.1101/2023.05.19.541476] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Anthocyanins are pigments and nutrients in red pears regulated by BBX family genes. Herein, we characterized a 14-nucleotide deletion mutation in the coding region of thePpBBX24gene from ‘Red Zaosu’ pear (Pyrus pyrifoliaWhite Pear Group), namedPpbbx24-del. Genetic and biochemical approaches were used to compare the roles of PpBBX24 and Ppbbx24-del in anthocyanin accumulation.Ppbbx24-delplayed a positive role in anthocyanin biosynthesis of the ‘Red Zaosu’ pear peel by light treatment. Functional analyses based on overexpression in tobacco and transient overexpression in pear fruit peels showed thatPpbbx24-delpromoted anthocyanin accumulation. Cyanidin and peonidin were major differentially expressed anthocyanins, and transcript levels of some structural genes in the anthocyanin biosynthesis pathway were significantly increased. Protein interaction assays showed that PpBBX24 was located in the nucleus and interacted with PpHY5, whereas Ppbbx24-del was colocalized in the nucleoplasm and did not interact with PpHY5. PpHY5 and Ppbbx24-del had positive regulatory effects on the expression ofPpCHS,PpCHI, andPpMYB10when acting alone, but had cumulative effects on gene activation when acting simultaneously. Alone, PpBBX24 had no significant effect on the expression ofPpCHS,PpCHI, orPpMYB10, whereas it inhibited the activation effects of PpHY5 on downstream genes when it existed with PpHY5. Our study demonstrated that mutant Ppbbx24-del positively regulates the anthocyanin accumulation in pear. The results of this study clarify the mechanism and enrich the regulatory network of anthocyanin biosynthesis, which lays a theoretical foundation forPpbbx24-deluse to create red pear cultivars.

Combined Transcriptome and Metabolome Analysis Reveals Adaptive Defense Responses to DON Induction in Potato

Phytophthora infestans poses a serious threat to potato production, storage, and processing. Understanding plant immunity triggered by fungal elicitors is important for the effective control of plant diseases. However, the role of the potato stress response to Fusarium toxin deoxynivalenol (DON)-induced stress is still not fully understood. In this study, the metabolites of DON-treated potato tubers were studied for four time intervals using UPLC-MS/MS. We identified 676 metabolites, and differential accumulation metabolite analysis showed that alkaloids, phenolic acids, and flavonoids were the major differential metabolites that directly determined defense response. Transcriptome data showed that differentially expressed genes (DEGs) were significantly enriched in phenylpropane and flavonoid metabolic pathways. Weighted gene co-expression network analysis (WGCNA) identified many hub genes, some of which modulate plant immune responses. This study is important for understanding the metabolic changes, transcriptional regulation, and physiological responses of active and signaling substances during DON induction, and it will help to design defense strategies against Phytophthora infestans in potato.

De novo full length transcriptome analysis of a naturally caffeine-free tea plant reveals specificity in secondary metabolic regulation

Tea plants are crops with economic, health and cultural value. Catechin, caffeine and theanine are the main secondary metabolites of taste. In the process of germplasm collection, we found a resource in the Sandu Aquatic Autonomous County of Guizhou (SDT) that possessed significantly different characteristic metabolites compared with the cultivar 'Qiancha 1'. SDT is rich in theobromine and theophylline, possesses low levels of (-)-epicatechin-3-gallate, (-)-epigallocatechin-3-gallate, and theanine content, and is almost free of caffeine. However, research on this tea resource is limited. Full-length transcriptome analysis was performed to investigate the transcriptome and gene expression of these metabolites. In total, 78,809 unique transcripts were obtained, of which 65,263 were complete coding sequences. RNA-seq revealed 3415 differentially expressed transcripts in the tender leaves of 'Qiancha 1' and 'SDT'. Furthermore, 2665, 6231, and 2687 differentially expressed transcripts were found in different SDT tissues. These differentially expressed transcripts were enriched in flavonoid and amino acid metabolism processes. Co-expression network analysis identified five modules associated with metabolites and found that genes of caffeine synthase (TCS) may be responsible for the low caffeine content in SDT. Phenylalanine ammonia lyase (PAL), glutamine synthetase (GS), glutamate synthase (GOGAT), and arginine decarboxylase (ADC) play important roles in the synthesis of catechin and theanine. In addition, we identified that ethylene resposive factor (ERF) and WRKY transcription factors may be involved in theanine biosynthesis. Overall, our study provides candidate genes to improve understanding of the synthesis mechanisms of these metabolites and provides a basis for molecular breeding of tea plant.

ZjFAS2 is involved in the fruit coloration in Ziziphus jujuba Mill. by regulating anthocyanin accumulation

Fruit color is one of the most important traits of jujube (Ziziphus jujuba Mill.). However, the differences in the pigments of different varieties of Jujube are not well studied. In addition, the genes responsible for fruit color and their underlying molecular mechanisms remain unclear. In this study, two jujube varieties, namely "Fengmiguan" (FMG) and "Tailihong" (TLH), were considered. The metabolites from jujube fruits were investigated using ultra-high-performance liquid chromatography/tandem mass spectrometry. Transcriptome was used to screen anthocyanin regulatory genes. The gene function was confirmed by overexpression and transient expression experiments. The gene expression was analyzed by quantitative reverse transcription polymerase chain reaction analyses and subcellular localization. Yeast-two-hybrid and bimolecular fluorescence complementation were used to screen and identify the interacting protein. These cultivars differed in color owing to their respective anthocyanin accumulation patterns. Three and seven types of anthocyanins were found in FMG and TLH, respectively, which played a key role in the process of fruit coloration. ZjFAS2 positively regulates anthocyanin accumulation. The expression profile of ZjFAS2 exhibited its different expression trends in different tissues and varieties. Subcellular localization experiments showed that ZjFAS2 was localized to the nucleus and membrane. A total of 36 interacting proteins were identified, and the possibility of ZjFAS2 interacting with ZjSHV3 to regulate jujube fruit coloration was studied. Herein, we investigated the role of anthocyanins in the different coloring patterns of the jujube fruits and provided a foundation for elucidating the molecular mechanism underlying jujube fruit coloration.

Identification of circular RNAs of Cannabis sativa L. potentially involved in the biosynthesis of cannabinoids

We identified circRNAs in the Cannabis sativa L. genome and examined their association with 28 cannabinoids in three tissues of C. sativa. Nine circRNAs are potentially involved in the biosynthesis of six cannabinoids. Cannabis sativa L. has been widely used in the production of medicine, textiles, and food for over 2500 years. The main bioactive compounds in C. sativa are cannabinoids, which have multiple important pharmacological actions. Circular RNAs (circRNAs) play essential roles in growth and development, stress resistance, and the biosynthesis of secondary metabolites. However, the circRNAs in C. sativa remain unknown. In this study, to explore the role of circRNAs in cannabinoid biosynthesis, we performed RNA-Seq and metabolomics analysis on the leaves, roots, and stems of C. sativa. We identified 741 overlapping circRNAs by three tools, of which 717, 16, and 8 circRNAs were derived from exonic, intronic, and intergenic, respectively. Functional enrichment analysis indicated that the parental genes (PGs) of circRNAs were enriched in many processes related to biological stress responses. We found that most of the circRNAs showed tissue-specific expression and 65 circRNAs were significantly correlated with their PGs (P < 0.05, |r|≥ 0.5). We also determined 28 cannabinoids by High-performance liquid chromatography-ESI-triple quadrupole-linear ion trap mass spectrometry. Ten circRNAs, including ciR0159, ciR0212, ciR0153, ciR0149, ciR0016, ciR0044, ciR0022, ciR0381, ciR0006, and ciR0025 were found to be associated with six cannabinoids by weighted gene co-expression network analysis. Twenty-nine of 53 candidate circRNAs, including 9 cannabinoids related were validated successfully using PCR amplification and Sanger sequencing. Taken together, all these results would help to enhance our acknowledge of the regulation of circRNAs, and lay the foundation for breeding new C. sativa cultivars with high cannabinoids through manipulating circRNAs.

Spatial and Temporal Distribution of Phenolic and Flavonoid Compounds in Sour Jujube (Ziziphus. Acidojujuba Cheng et Liu) and Their Antioxidant Activities

In order to comprehensively analyze the antioxidant substances in sour jujube, total phenolic content (TPC) and total flavonoids contents (TFC) in different organs, including stem, leaf, flower, fruit pulp, and seed were analyzed for their contents and antioxidant activities. The results showed that leaves possessed significantly higher TPC and TFC (20.4 and 20.5 mg/g, respectively) than the other organs and have the highest antioxidant activity, which were also higher than the wild blueberry (A well-known for its high TPC). Subsequently, the variations in the antioxidant content and antioxidant activity of leaves were analyzed during leaf development. TPC in leaves sampled in may and august were significantly higher than that in other months, while the highest one was found in may. The n-hexane, ethyl acetate, n-butanol, and water fractions obtained from the main methanol extract of sour jujube leaves were evaluated for TPC and TFC and their antioxidant activity and it was found that ethyl acetate fraction displayed the highest TPC and TFC (184.5 and 193.3 mg/g, respectively), as well as the best antioxidant activity. In addition, using LC-MS and HPLC, ethyl acetate fraction was analyzed from qualitative and quantitative aspects; 31-one phenolic compounds, including catechin (33.0 mg/g), epigallocatechin (15.3 mg/g), quercetin 3-O-glucoside (11.4 mg/g), naringenin (6.7 mg/g), esculetin (4.8 mg/g), and chlorogenic acid (4.6 mg/g) were identified. Catechin, esculetin, epigallocatechin, chlorogenic acid, quercetin 3-O-glucoside, and naringenin exhibited high antioxidant activity. These results provide a theoretical basis for further study and utilization of flavonoid and polyphenols in sour jujube.

Transcriptome and metabolome response of eggplant against Ralstonia solanacearum infection

Bacterial wilt is a soil-borne disease that represents ubiquitous threat to Solanaceae crops. The whole-root transcriptomes and metabolomes of bacterial wilt-resistant eggplant were studied to understand the response of eggplant to bacterial wilt. A total of 2,896 differentially expressed genes and 63 differences in metabolites were identified after inoculation with Ralstonia solanacearum. Further analysis showed that the biosynthesis pathways for phytohormones, phenylpropanoids, and flavonoids were altered in eggplant after inoculation with R. solanacearum. The results of metabolomes also showed that phytohormones played a key role in eggplant response to bacterial wilt. Integrated analyses of the transcriptomic and metabolic datasets indicated that jasmonic acid (JA) content and gene involved in the JA signaling pathway increased in response to bacterial wilt. These findings remarkably improve our understanding of the mechanisms of induced defense response in eggplant and will provide insights intothe development of disease-resistant varieties of eggplant.

Metabolomics-centered mining of plant metabolic diversity and function: Past decade and future perspectives

Plants are natural experts in organic synthesis, being able to generate large numbers of specific metabolites with widely varying structures that help them adapt to variable survival challenges. Metabolomics is a research discipline that integrates the capabilities of several types of research including analytical chemistry, statistics, and biochemistry. Its ongoing development provides strategies for gaining a systematic understanding of quantitative changes in the levels of metabolites. Metabolomics is usually performed by targeting either a specific cell, a specific tissue, or the entire organism. Considerable advances in science and technology over the last three decades have propelled us into the era of multi-omics, in which metabolomics, despite at an earlier developmental stage than genomics, transcriptomics, and proteomics, offers the distinct advantage of studying the cellular entities that have the greatest influence on end phenotype. Here, we summarize the state of the art of metabolite detection and identification, and illustrate these techniques with four case study applications: (i) comparing metabolite composition within and between species, (ii) assessing spatio-temporal metabolic changes during plant development, (iii) mining characteristic metabolites of plants in different ecological environments and upon exposure to various stresses, and (iv) assessing the performance of metabolomics as a means of functional gene identification , metabolic pathway elucidation, and metabolomics-assisted breeding through analyzing plant populations with diverse genetic variations. In addition, we highlight the prominent contributions of joint analyses of plant metabolomics and other omics datasets, including those from genomics, transcriptomics, proteomics, epigenomics, phenomics, microbiomes, and ion-omics studies. Finally, we discuss future directions and challenges exploiting metabolomics-centered approaches in understanding plant metabolic diversity.

The Jasmine (Jasminum sambac) Genome Provides Insight into the Biosynthesis of Flower Fragrances and Jasmonates

Jasminum sambac (jasmine flower), a world-renowned plant appreciated for its exceptional flower fragrance, is of cultural and economic importance. However, the genetic basis of its fragrance is largely unknown. Here, we present the first de novo genome of J. sambac with 550.12 Mb (scaffold N50 = 40.10 Mb) assembled into 13 pseudochromosomes. Terpene synthase genes associated with flower fragrance are significantly amplified in the form of gene clusters through tandem duplications in the genome. Gene clusters within the salicylic acid/benzoic acid/theobromine (SABATH) and BAHD superfamilies were identified as related to the biosynthesis of phenylpropanoid/benzenoid compounds. Several key genes involved in jasmonate biosynthesis were duplicated, causing increased copy numbers. In addition, multi-omics analyses identified various aromatic compounds and many genes involved in fragrance biosynthesis pathways. Furthermore, the roles of JsTPS3 in β-ocimene biosynthesis, as well as JsAOC1 and JsAOS in jasmonic acid biosynthesis, were functionally validated. The genome assembled in this study for J. sambac offers a basic genetic resource for studying floral scent and jasmonate biosynthesis and provides a foundation for functional genomic research and variety improvements in Jasminum.

Combined Metabolomic and Quantitative RT-PCR Analyses Revealed the Synthetic Differences of 2-Acetyl-1-pyrroline in Aromatic and Non-Aromatic Vegetable Soybeans

Aroma is an important economic trait of vegetable soybeans, which greatly influences their market value. The 2-acetyl-1-pyrroline (2AP) is considered as an important substance affecting the aroma of plants. Although the 2AP synthesis pathway has been resolved, the differences of the 2AP synthesis in the aromatic and non-aromatic vegetable soybeans are unknown. In this study, a broad targeted metabolome analysis including measurement of metabolites levels and gene expression levels was performed to reveal pathways of aroma formation in the two developmental stages of vegetable soybean grains [35 (S5) and 40 (S6) days after anthesis] of the 'Zhexian No. 8' (ZX8, non-aromatic) and ZK1754 (aromatic). The results showed that the differentially accumulated metabolites (DAMs) of the two varieties can be classified into nine main categories including flavonoids, lipids, amino acids and derivatives, saccharides and alcohols, organic acids, nucleotides and derivatives, phenolic acids, alkaloids and vitamin, which mainly contributed to their phenotypic differences. Furthermore, in combination with the 2AP synthesis pathway, the differences of amino acids and derivatives were mainly involved in the 2AP synthesis. Furthermore, 2AP precursors' analysis revealed that the accumulation of 2AP mainly occurred from 1-pyrroline-5-carboxylate (P5C), not 4-aminobutyraldehyde (GABald). The quantitative RT-PCR showed that the associated synthetic genes were 1-pyrroline-5-carboxylate dehydrogenase (P5CDH), ∆¹-pyrroline-5-carboxylate synthetase (P5CS), proline dehydrogenase (PRODH) and pyrroline-5-carboxylate reductase (P5CR), which further verified the synthetic pathway of 2AP. Furthermore, the betaine aldehyde dehydrogenase 2 (GmBADH2) mutant was not only vital for the occurrence of 2AP, but also for the synthesis of 4-aminobutyric acid (GABA) in vegetable soybean. Therefore, the differences of 2AP accumulation in aromatic and non-aromatic vegetable soybeans have been revealed, and it also provides an important theoretical basis for aromatic vegetable soybean breeding.

Photosynthetic carbon and nitrogen metabolism of Camellia oleifera Abel during acclimation to low light conditions

Tea-oil tree (Camellia oleifera Abel) is an important woody oil crop with high economic value. However, it has low photosynthetic production considering the low light intensity of its growth environment. To understand the acclimation mechanism of tea-oil trees to low light conditions, three light intensity treatments were conducted: high light (450-500 μmol. m^-2. s^-1), medium light (180-200 μmol. m^-2. s^-1), and low light (45-50 μmol. m^-2. s^-1). The carbon (C) and nitrogen (N) metabolism network were constructed by investigating the leaf anatomy, photosynthetic characteristics, N partitioning, transcriptome and metabolome. Results demonstrated that a larger proportion light energy was used for photochemical reactions in an environment with lower light intensity, which resulted in an increase in photosystem II photochemical efficiency and instantaneous light use efficiency (LUE) at the leaf level. As the light intensity increased, decreased electron transfer and carboxylation efficiencies, photorespiration and dark respiration rates, LUE at plant level, and N use efficiency (P_NUE) were observed. Leaves trended to harvest more light using higher expression levels of light-harvesting protein genes, higher chlorophyll content, more granum and more tightly stacked granum lamella under lower light intensity. At transcriptional and metabolic levels, the TCA cycle, and the synthesis of starch and saccharides were weakened as light intensity decreased, while the Calvin cycle did not show the regularity between different treatments. Less N was distributed in Rubisco, respiration, and cell wall proteins as light decreased. Storage N was prominently accumulated in forms of amino acids (especially L-arginine) and amino acid derivatives as under medium and low light environments, to make up for C deficiency. Therefore, tea-oil trees actively improve light-harvesting capacity and enlarges the storage N pool to adapt to a low light environment, at the cost of a decrease of photosynthetic C assimilation and P_NUE.

Transcriptomic and metabolomic analyses reveal the main metabolites in Dendrobium officinale leaves during the harvesting period

Dendrobium officinale, which is a medicine food homology plant, contains many metabolites, especially polysaccharides and flavonoids. Unlike flowers and stems, which are the most frequently harvested organs for a variety of uses, leaves tend to be discarded. This study assessed main metabolites in leaves to identify the most appropriate timing of collection during harvest, which was divided into three stages (S1-S3: 8, 10, and 11 months after sprouting, respectively). Metabolomic and transcriptomic analyses of S1-S3 were performed. Water-soluble polysaccharides (WSPs), flavonoids and free amino acids (FAAs) were detected in leaves. WSPs decreased from S1 to S3 but flavonoids and some FAAs (e.g., phophoserine) increased from S1 to S2, then decreased from S2 to S3. In all three stages, mannose was the dominant monosaccharide among WSPs, followed by glucose. In S2, 35 flavonoids were identified, the most abundant being rutin, schaftoside and vitexin, while 34 FAAs were identified in all three stages, the most abundant being tyrosine, phosphoserine and alanine. A total of 2584, 3414 and 2032 differentially expressed genes (DEGs) were discovered in S1 vs S2, S1 vs S3 and S1 vs S3, respectively. Correlation analysis revealed that five DEGs (DoSUS, DoXYLA, DoFRK, DoGMP, and DoCSLA), two DEGs (DoDFR, and DoANS) and a single DEG (DoPGAM) were involved in the metabolism of WSPs, flavonoids and phosphoserine, respectively. The findings of this study lay a foundation for the commercial exploitation of metabolites in the harvested leaves of D. officinale, and the use of detected DEGs in applied genetic studies.

Biosynthetic regulatory network of flavonoid metabolites in stems and leaves of Salvia miltiorrhiza

Flavonoid secondary metabolites can treat and prevent many diseases, but systematic studies on regulation of the biosynthesis of such metabolites in aboveground parts of Salvia miltiorrhiza are lacking. In this study, metabonomic and transcriptomic analyses of different S. miltiorrhiza phenotypes were conducted to explore pathways of synthesis, catalysis, accumulation, and transport of the main flavonoid secondary metabolites regulating pigment accumulation. Tissue localization and quantitative analysis of flavonoid secondary metabolites were conducted by laser scanning confocal microscopy (LSCM). A total 3090 differentially expressed genes were obtained from 114,431 full-length unigenes in purple and green phenotypes, and 108 functional genes were involved in flavonoid biosynthesis. Five key phenylpropane structural genes (PAL, 4CL, ANS, 3AT, HCT) were highly differentially expressed, and four transcription factor genes (MYB, WRKY, bHLH, bZiP) were identified. In addition, six GST genes, nine ABC transporters, 22 MATE genes, and three SNARE genes were detected with key roles in flavonoid transport. According to LSCM, flavonoids were mainly distributed in epidermis, cortex, and collenchyma. Thus, comprehensive and systematic analyses were used to determine biosynthesis, accumulation, and transport of flavonoids in stems and leaves of different S. miltiorrhiza phenotypes. The findings will provide a reference for flavonoid production and cultivar selection.

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.

An integrated metabolome and transcriptome approach reveals the fruit flavor and regulatory network during jujube fruit development

The fruit flavor is a key economic value attribute of jujube. Here we compared metabolomes and transcriptomes of "Mazao" (ST) and "Ping'anhuluzao" (HK) with unique flavors during fruit development. We identified 437 differential metabolites, mainly sugars, acids, and lipids. Fructose, glucose, mannose and citric acid, and malic acid are the determinants of sugar and acid taste of jujube fruit. Based on the transcriptome, 16,245 differentially expressed genes (DEGs) were identified, which were involved in "glucosyltransferase activity," "lipid binding," and "anion transmembrane transporter activity" processes. Both transcriptome and metabolome showed that developmental stages 2 and 3 were important transition periods for jujube maturation. Based on WGCNA and gene-metabolite correlation analysis, modules, and transcription factors (ZjHAP3, ZjTCP14, and ZjMYB78) highly related to sugar and acid were identified. Our results provide new insights into the mechanism of sugar and acid accumulation in jujube fruit and provide clues for the development of jujube with a unique flavor.

Widely targeted metabolomics reveals stamen petaloid tissue of Paeonia lactiflora Pall. being a potential pharmacological resource

Paeonia lactiflflora Pall. has a long edible and medicinal history because of the very high content of biologically active compounds. However, little information is available about the metabolic basis of pharmacological values of P. lactiflora flowers. In this study, we investigated metabolites in the different parts of P. lactiflora flowers, including petal, stamen petaloid tissue and stamen, by widely targeted metabolomics approach. A total of 1102 metabolites were identified, among which 313 and 410 metabolites showed differential accumulation in comparison groups of petal vs. stamen petaloid tissue and stamen vs. stamen petaloid tissue. Differential accumulated metabolites analysis and KEGG pathway analysis showed that the flavonoids were the most critical differential metabolites. Furthermore, difference accumulation of flavonoids, phenolic acids, tannins and alkaloids might lead to the differences in antioxidant activities and tyrosinase inhibition effects. Indeed, stamen petaloid tissue displayed better antioxidant and anti-melanin production activities than petal and stamen through experimental verification. These results not only expand our understanding of metabolites in P. lactiflora flowers, but also reveal that the stamen petaloid tissues of P. lactiflora hold the great potential as promising ingredients for pharmaceuticals, functional foods and skincare products.

Identification of Medicinal Compounds of Fagopyri Dibotryis Rhizome from Different Origins and Its Varieties Using UPLC-MS/MS-Based Metabolomics

Fagopyrum dibotrys, being native to southwest China, is widely distributed in Yunnan, Guizhou Provinces and Chongqing City. However, the quality of medicinal materials growing in different origins varies greatly, and cannot meet the market demand for high-quality F. dibotrys. In this study, 648 metabolites were identified, and phenolic compounds of F. dibotrys from different origins were clearly separated by principal component analysis (PCA). Our results suggested that the medicinal differences of F. dibotrys from different origins can be elucidated via the variations in the abundance of the phenolic and flavonoid compounds. We found that the epicatechin, total flavonoids and total tannin content in Yunnan Qujing (YQ) and Yunnan Kunming (YK) were higher than those in Chongqing Shizhu (CS), Chongqing Fuling (CF) and Guizhou Bijie (GB), suggesting that Yunnan Province can be considered as one of the areas that produce high-quality medicinal materials. Additionally, 1,6-di-O-galloyl-β-D-glucose, 2,3-di-O-galloyl-D-glucose and gallic acid could be used as ideal marker compounds for the quality control of F. dibotrys from different origins caused by metabolites, and the F. dibotrys planted in Yunnan Province is well worth exploiting.

Identification of key genes controlling L-ascorbic acid during Jujube (Ziziphus jujuba Mill.) fruit development by integrating transcriptome and metabolome analysis

Chinese jujube (Ziziphus jujuba) is a vital economic tree native to China. Jujube fruit with abundant L-Ascorbic Acid (AsA) is an ideal material for studying the mechanism of AsA biosynthesis and metabolism. However, the key transcription factors regulating AsA anabolism in jujube have not been reported. Here, we used jujube variety "Mazao" as the experimental material, conducted an integrative analysis of transcriptome and metabolome to investigate changes in differential genes and metabolites, and find the key genes regulating AsA during jujube fruit growth. The results showed that AsA was mostly synthesized in the young stage and enlargement stage, ZjMDHAR gene takes an important part in the AsA recycling. Three gene networks/modules were highly correlated with AsA, among them, three genes were identified as candidates controlling AsA, including ZjERF17 (LOC107404975), ZjbZIP9 (LOC107406320), and ZjGBF4 (LOC107421670). These results provide new directions and insights for further study on the regulation mechanism of AsA in jujube.

UPLC-MS/MS Profile Combined With RNA-Seq Reveals the Amino Acid Metabolism in Zanthoxylum bungeanum Leaves Under Drought Stress

Zanthoxylum bungeanum leaves have a unique taste and incomparable nutritional value and hence are popular as a food item and traditional medicine in China. However, the studies on the metabolites in Z. bungeanum leaves are quite limited, especially for amino acids. Therefore, this study explored the amino acid component in Z. bungeanum leaves and also the accumulation mechanism under drought stress in two Z. bungeanum cultivars using the widely targeted metabolome combined with transcriptome analysis. A total of 56 amino acids and their derivatives were identified in Z. bungeanum leaves, including eight essential amino acids. The total amino acid content with most individual amino acids increased under progressive drought stress. More differentially accumulated amino acids (DAAs) and differentially expressed genes (DEGs) were found in FJ (Z. bungeanum cv. ‘Fengjiao’) than in HJ (Z. bungeanum cv. ‘Hanjiao’). The orthogonal projections to latent structures discriminant analysis identified nine and seven indicator DAAs in FJ and HJ leaves, respectively. The weighted gene co-expression network analysis (WGCNA) showed that the green module was significantly correlated with most indicator DAAs and revealed the important role of FBA3, DELTA-OAT, PROC, and 15 transcription factor genes in regulating the amino acid synthesis. Furthermore, the correlation analysis and redundancy analysis (RDA) identified four candidate synthesis genes (ASNS, AK, ASPS, and PK) in amino acid biosynthesis pathway. This study provided useful information for the development of Z. bungeanum leaves in food and nutrition industry and also laid the foundations for future molecular breeding.

Iminodisuccinic acid enhances antioxidant and mineral element accumulation in young leaves of Ziziphus jujuba

Jujube leaf tea, which is made from the young leaves of Ziziphus jujuba, is a novel functional herb tea or infusion that inhibits the central nervous system. In the current study, the effects of iminodisuccinic acid (IDS), as a metal complexing agent, on mineral element content, oxidative damage, antioxidant enzyme activities, and antioxidant accumulation in the young and mature leaves of Z. jujuba were investigated. Results demonstrated that foliar fertilization with ionic (FeCl2 and ZnCl2) and chelated (Fe-IDS and Zn-IDS) fertilizers could drastically enhance iron and zinc contents, coupled with increased vitamin C level, glutathione accumulation, total phenolic content, and total antioxidant capacity (evaluated based on the Fe3+ reducing power of leaf extracts), compared with the control, particularly in young leaves. However, chelated fertilizers considerably reduced the chlorophyll level, H2O2 content, and lipid peroxidation rate than ionic fertilizers, particularly in young leaves. Compared with the control, chelated fertilizers induced greater superoxide dismutase and catalase activities, particularly in young leaves. Moreover, decreased enzyme activities were observed in the ionic fertilizer-treated leaves compared with the control-treated leaves. Thus, using a chelating agent could improve the accumulation of mineral elements and antioxidants in young leaves by reducing metal-mediated reactive oxygen species toxicity.

Combined Transcriptome and Metabolome Analyses Reveal Candidate Genes Involved in Tangor (Citrus reticulata × Citrus sinensis) Fruit Development and Quality Formation

Tangor, an important citrus type, is a hybrid of orange and mandarin and possesses their advantageous characteristics. Fruit quality is an important factor limiting the development of the citrus industry and highly depends on fruit development and ripening programs. However, fruit development and quality formation have not been completely explored in mandarin-orange hybrids. We sequenced the metabolome and transcriptome of three mandarin-orange hybrid cultivars at the early fruiting [90 days after full bloom (DAFB)], color change (180 DAFB), and ripening (270 DAFB) stages. Metabolome sequencing was performed to preliminarily identify the accumulation patterns of primary and secondary metabolites related to fruit quality and hormones regulating fruit development. Transcriptome analysis showed that many genes related to primary metabolism, secondary metabolism, cell wall metabolism, phytohormones, and transcriptional regulation were up-regulated in all three cultivars during fruit development and ripening. Additionally, multiple key genes were identified that may play a role in sucrose, citric acid and flavonoid accumulation, cell wall modification, and abscisic acid signaling, which may provide a valuable resource for future research on enhancement of fruit quality of hybrid citrus. Overall, this study provides new insights into the molecular basis of pulp growth and development regulation and fruit quality formation in mandarin-orange hybrids.

Comparative metabolomics provides novel insights into the basis of petiole color differences in celery (Apiumgraveolens L.)

Plant metabolites are important for plant development and human health. Plants of celery (Apiumgraveolens L.) with different-colored petioles have been formed in the course of long-term evolution. However, the composition, content distribution, and mechanisms of accumulation of metabolites in different-colored petioles remain elusive. Using ultra-high performance liquid chromatography-tandem mass spectrometry (UHPLC-MS/MS), 1159 metabolites, including 100 lipids, 72 organic acids and derivatives, 83 phenylpropanoids and polyketides, and several alkaloids and terpenoids, were quantified in four celery cultivars, each with a different petiole color. There were significant differences in the types and contents of metabolites in celery with different-colored petioles, with the most striking difference between green celery and purple celery, followed by white celery and green celery. Annotated analysis of metabolic pathways showed that the metabolites of the different-colored petioles were significantly enriched in biosynthetic pathways such as anthocyanin, flavonoid, and chlorophyll pathways, suggesting that these metabolic pathways may play a key role in determining petiole color in celery. The content of chlorophyll in green celery was significantly higher than that in other celery cultivars, yellow celery was rich in carotenoids, and the content of anthocyanin in purple celery was significantly higher than that in the other celery cultivars. The color of the celery petioles was significantly correlated with the content of related metabolites. Among the four celery cultivars, the metabolites of the anthocyanin biosynthesis pathway were enriched in purple celery. The results of quantitative real-time polymerase chain reaction (qRT-PCR) suggested that the differential expression of the chalcone synthase (CHS) gene in the anthocyanin biosynthesis pathway might affect the biosynthesis of anthocyanin in celery. In addition, HPLC analysis revealed that cyanidin is the main pigment in purple celery. This study explored the differences in the types and contents of metabolites in celery cultivars with different-colored petioles and identified key substances for color formation. The results provide a theoretical basis and technical support for genetic improvement of celery petiole color.