An integrated transcriptome and metabolome approach reveals the accumulation of taste-related metabolites and gene regulatory networks during watermelon fruit development

Cucurbitacins B, E and I Concentrations and Relationship with Drought Tolerance in Bottle Gourd [Lagenaria siceraria (Molina) Standl.]

Bottle gourd [Lagenaria siceraria (Molina) Standl.]) is a relatively drought-tolerant cucurbit due to the high composition of unique biochemical compositions, including cucurbitacin. The objective of this study was to determine the concentrations of cucurbitacins in bottle gourd and their relationship to drought tolerance. The study assessed 12 bottle gourd accessions grown under two moisture levels (i.e., non-stressed (NS) and drought-stressed (DS)) and three drought stress intensities (i.e., mild, moderate, and severe) using a 12 × 2 × 3 factorial experiment designed in a randomized complete block design with three replications. Control studies were undertaken under glasshouse conditions. The content of cucurbitacins B, E, and I were quantified in leaves and roots using high-performance liquid Cchromatography-mass spectrometry (HPLC-MS). The free radical scavenging activities of pure cucurbitacins B, E, and I were quantified using 2,2-diphenyl-1-picrylhydrazyl (DPPH) and a ferrulic acid power assay (FRAP). Results revealed that cucurbitacins B and I were present in accessions BG-48, BG-58, BG-70, BG-78, BG-79, BG-81, BG-52, and GC in leaves and roots under DS condition. The contents of cucurbitacins B and I were enhanced under increased drought intensity for accessions BG-48, BG-81, and GC. In all the leaf and root samples, cucurbitacin E was not detectable. Based on the DPPH test, pure cucurbitacins I, B, and E reduced free radicals at maximum values of 78, 60, and 66%, respectively. Based on the FRAP assay, pure cucurbitacins I, B, and E had maximum ferric-reducing powers of 67, 62, and 48%. Additionally, cucurbitacin I recorded the highest antioxidant activity compared to cucurbitacins B and E. Increased cucurbitacin accumulation and antioxidant properties indicate their role in minimising cell damage caused by oxidative stress under drought-stressed environments. The present study revealed that cucurbitacins B and I serve as novel biochemical markers for screening drought tolerance in bottle gourd or related cucurbits.

Sugar composition and transcriptome analysis in developing 'Fengtang' plum (Prunus salicina Lindl.) reveal candidate genes regulating sugar accumulation

Sweetness is an important attribute of fruit quality, which directly affects consumers' preference for fresh fruit and is mostly determined by carbohydrate composition. 'Fengtang' plum (Prunus salicina Lindl.) is recognized for its high soluble sugar content, but the sugar composition and the molecular mechanisms underlying sugar overproduction are not fully understood. In this work, the sugar components were analyzed using gas chromatography-mass spectrometry combined with transcription profiles from RNA-sequencing and Quantitative Real-time PCR during fruit development. The target metabolic group showed that sucrose was the dominant sugar component in mature fruit, followed by glucose, fructose, and sorbitol. Based on the transcriptome data and qRT-PCR validation, we identified 12 key structural genes that significantly responded to corresponding component accumulation: sucrose synthase (PsSUS4), sucrose phosphate synthase (PsSPS2), neutral invertase (PsNINV1/3/4), phosphoglucomutase (PsPGM1), UTP-glucose-1-phosphate uridylyl transferase (PsUGP1/2), hexose kinase (PsHXK1/3), sugar transport protein (PsSTP1), and Sugars Will Eventually be Exported Transporter (PsSWEET4). In which PsSUS4 and PsSPS2, whose encoding proteins immediately catalyze sucrose synthesis, were selected to be silenced using the virus-induced gene silencing technology. Silencing of PsSUS4 or PsSPS2 resulted in decreased sucrose content by 27.6% and 8%, respectively, compared with the control, verifying their important roles in sucrose accumulation. Subsequently, sugar metabolism networks in this high-sugar plum were constructed with 12 key structural genes, 72 putative transcription factors, and 4 major sugar components. These results might facilitate a better understanding of the molecular mechanisms of sugar accumulation in 'Fengtang' plum and provide a framework for future fruit quality improvement.

Identification and quantification of cucurbitacin in watermelon frost using molecular networking integrated with ultra-high-performance liquid chromatography-tandem mass spectrometry

Watermelon frost, a traditional Chinese medicine produced using watermelon and Glauber's salt, has been widely used for the therapy of oral and throat disorders. Watermelon contains various phytochemical compounds including cucurbitacins and their glycoside derivatives, which have attracted considerable attention because of their medicinal values. However, whether the composition of cucurbitacins existed in watermelon frost was rarely reported. In this study, three cucurbitacins including cucurbitacin B, isocucurbitacin B, and cucurbitacin E were found from watermelon frost extract assisted by ultra-high-performance liquid chromatography-tandem mass spectrometry and molecular networking guided strategy, and the compounds were verified using standard solutions. Furthermore, a quantification method for simultaneously targeted analysis of cucurbitacins was established using ultra-high-performance liquid chromatography-tandem mass spectrometry operating in the multiple reaction monitoring mode. Among them, cucurbitacin B and cucurbitacin E in watermelon frost samples were determined, and the concentrations were 3.78 ± 0.18 and 0.86 ± 0.19 ng/ml, respectively. While isocucurbitacin B was not detected due to the lower content possibly. In conclusion, ultra-high-performance liquid chromatography-tandem mass spectrometry combined with molecular networking is a very useful technique for the rapid identification of unknown cucurbitacin components in watermelon frost.

Multi-omics integration to explore the molecular insight into the volatile organic compounds in watermelon

A range of volatile organic compounds played an important role in the formation of watermelon fruit aroma, while due to the low content and difficulty in detection, it is often neglected in watermelon breeding programs, resulting in a decline in fruit flavor. VOCs in the flesh of 194 watermelon accessions and seven cultivars at four developmental stages were determined by SPME-GC-MS. Ten metabolites with significant differences in the natural population and positive accumulation during fruit development are considered to be the key metabolite related to watermelon fruit aroma. And the link between metabolite and, flesh color and sugar content by correlation analysis was established. The results of the genome-wide association study showed that (5E)-6,10-dimethylundeca-5,9-dien-2-one, and 1-(4-methylphenyl) ethanone were colocalized with watermelon flesh color on chromosome 4, which may be regulated by LCYB and CCD. (E)-4-(2,6,6-trimethylcyclohexen-1-yl)but-3-en-2-one is the VOC produced by the cleavage of carotenoids, which has a positive correlation with the sugar content of the fruit, and the candidate gene Cla97C05G092490 on chromosome 5 may interact with PSY to influence the accumulation of this metabolite. In addition, Cla97C02G049790 (enol reductase), Cla97C03G051490 (omega-3 fatty acid desaturase gene), LOX, and ADH may play important roles in the synthesis of fatty acids and their derived VOCs. Taken together, our findings provide molecular insights into the accumulation and natural variation of VOCs in watermelon, and give data support for breeding watermelon cultivars with better flavor.

Combined Metabolome and Transcriptome Analysis Elucidates Sugar Accumulation in Wucai (Brassica campestris L.)

Wucai (Brassica campestris L.) is a leafy vegetable that originated in China, its soluble sugars accumulate significantly to improve taste quality during maturation, and it is widely accepted by consumers. In this study, we investigated the soluble sugar content at different developmental stages. Two periods including 34 days after planting (DAP) and 46 DAP, which represent the period prior to and after sugar accumulation, respectively, were selected for metabolomic and transcriptomic profiling. Differentially accumulated metabolites (DAMs) were mainly enriched in the pentose phosphate pathway, galactose metabolism, glycolysis/gluconeogenesis, starch and sucrose metabolism, and fructose and mannose metabolism. By orthogonal projection to latent structures-discriminant s-plot (OPLS-DA S-plot) and MetaboAnalyst analyses, D-galactose and β-D-glucose were identified as the major components of sugar accumulation in wucai. Combined with the transcriptome, the pathway of sugar accumulation and the interact network between 26 DEGs and the two sugars were mapped. CWINV4, CEL1, BGLU16, and BraA03g023380.3C had positive correlations with the accumulation of sugar accumulation in wucai. The lower expression of BraA06g003260.3C, BraA08g002960.3C, BraA05g019040.3C, and BraA05g027230.3C promoted sugar accumulation during the ripening of wucai. These findings provide insights into the mechanisms underlying sugar accumulation during commodity maturity, providing a basis for the breeding of sugar-rich wucai cultivars.

Integrated untargeted metabolome, full-length sequencing, and transcriptome analyses reveal insights into the fruit quality at different harvest times of Chaenomeles speciosa

Chaenomeles speciosa fruit is a homologous medicine and food plant with a long history of multiple uses. It could be harvested near maturity and last for a long time. However, the optimal harvest strategy of Chaenomeles speciosa for various uses is currently unavailable. Here, untargeted metabolome at different harvest times during maturation was investigated for the first time, and 896 metabolites, including sugars, organic acids, amino acids, and phenylpropanoids, were identified. Optimal harvesting methods were proposed for different purposes. During the early maturation stages (before 105 days after full bloom), Ch. speciosa fruit could be harvested as Chinesemedicine. Whereas as snacks and food, Ch. speciosa fruit might be harvested at late maturity (after 120 days after full bloom). In addition, the overall network was revealed by integrating full-length Iso-seq and transcriptomics (RNA-seq) to investigate the association between quality-associated metabolites and Chaenomeles speciosa fruit gene expression during maturation. A few putative genes were captured via screening, dissecting and correlation analysis with the quality-associated metabolites (including d-glucose, catechin, gallocatechin, and succinic acid). Overall, in addition to providing a harvesting strategy for food and medicine, we also investigated the metabolism and gene expression pattern of Chaenomeles speciosa fruit during maturation. This comprehensive data and analyses laid the foundation for further investigating potential regulatory mechanisms during harvest and provided a new possibility for its development and utilization.

Telomere-to-telomere genome assembly of bitter melon (Momordica charantia L. var. abbreviata Ser.) reveals fruit development, composition and ripening genetic characteristics

Momordica charantia L. var. abbreviata Ser. (Mca), known as bitter gourd or bitter melon, is a Momordica variety with medicinal value and belongs to the Cucurbitaceae family. In view of the lack of genomic information on bitter gourd and other Momordica species and to promote Mca genomic research, we assembled a 295.6-Mb telomere-to-telomere (T2T) high-quality Mca genome with six gap-free chromosomes after Hi-C correction. This genome is anchored to 11 chromosomes, which is consistent with the karyotype information, and comprises 98 contigs (N50 of 25.4 Mb) and 95 scaffolds (N50 of 25.4 Mb). The Mca genome harbors 19 895 protein-coding genes, of which 45.59% constitute predicted repeat sequences. Synteny analysis revealed variations involved in fruit quality during the divergence of bitter gourd. In addition, assay for transposase-accessible chromatin by high-throughput sequencing and metabolic analysis showed that momordicosides and other substances are characteristic of Mca fruit pulp. A combined transcriptomic and metabolomic analysis revealed the mechanisms of pigment accumulation and cucurbitacin biosynthesis in Mca fruit peels, providing fundamental molecular information for further research on Mca fruit ripening. This report provides a new genetic resource for Momordica genomic studies and contributes additional insights into Cucurbitaceae phylogeny.

Genome-wide association study provides genetic insights into natural variation in watermelon rind thickness and single fruit weight

Rind thickness and fruit weight are agronomic traits closely related to quality and yield, which have attracted much attention from consumers and breeders. However, the genetic mechanism of these two traits is still not well understood in natural populations. In this study, rind thickness and single fruit weight in 151 watermelon accessions were determined in 2019 and 2020, and genome-wide association analysis was performed by integrating phenotypic and genotype data. Abundant phenotypic variation was found in the test population, and the watermelon with thinner rind thickness tended to have smaller fruit weights. Five significant SNPs were closely associated with rind thickness on chromosome 2 by Genome-wide association study (GWAS), i.e., 32344170, 32321308, 32304738, 32328501, and 32311192. And there were 21 genes were annotated in the candidate interval, most notably, Cla97C02G044160 belonged to the MADS family, and part of the genes in this family played an important role in regulating organ size, further analysis of gene structure, gene expression level, and phylogenetic tree showed that Cla97C02G044160 was a candidate gene for regulating target traits. In conclusion, our study provides molecular insights into the natural variation of watermelon rind thickness and single fruit weight, meanwhile, providing data support for molecular marker-assisted breeding.

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.

Metabolome and Transcriptome Profiling Reveal That Four Terpenoid Hormones Dominate the Growth and Development of Sanghuangporus baumii

Sanghuangporus baumii is a traditional medicinal fungus that produces pharmacological terpenoids, but natural resources are insufficient for applications, and its growth and development mechanisms are poorly understood. Combining metabolomic and transcriptomic analyses, we found four terpenoid hormones and a central gene, isopentenyl diphosphate isomerase (IDI), involved in growth and development. Additionally, an exogenous hormone test was used to further confirm the importance of the four terpenoid hormones. Finally, hormone content determination and qRT−PCR were performed to explore the growth and development mechanism; we found thatcis-zeatin (CZ) plays a major role in the mycelia stage, trans-zeatin (TZ) and gibberellin A4 (GA4) are important in the primordia stage, GA4 is crucial for the fruiting bodies stage, and abscisic acid (ABA) may be a marker of maturity. The IDI gene was also found to affectterpenoid hormone content by regulating the relative gene transcript levels, thereby controlling morphological changes in S. baumii. Our results revealthe growth and development mechanisms of S. baumii and may promote the breeding and utilisation of high-quality varieties.

Integrative analysis of transcriptomic and metabolomic profiles reveals new insights into the molecular foundation of fruit quality formation in Citrullus lanatus (Thunb.) Matsum. & Nakai

In this study, an integrated transcriptome and metabolome analysis was used to explore the molecular foundation of fruit quality in two parent lines of Citrullus lanatus with distinct flesh characteristics, including ‘14-1’ (sweet, red, and soft) and ‘W600’ (bitter, light yellow, and firm), as well as the corresponding F1 population (bitter, light yellow, and firm). Numerous differentially expressed genes (DEGs) were identified in the fruit samples: 3,766 DEGs for ‘14-1’ vs. ‘W600’, 2,767 for ‘14-1’ vs. F1, and 1,178 for F1 vs. ‘W600’ at the transition stage; and 4,221 for ‘14-1’ vs. ‘W600’, 2,447 for ‘14-1’ vs. F1, and 446 for F1 vs. ‘W600’ at the maturity stage. Weighted gene co-expression network analysis (WGCNA) revealed that a gene module including 1,111 DEGs was closely associated with flesh taste and color, and another gene module including 1,575 DEGs contributed significantly to flesh texture. The metabolomic results showed that there were 447 differential metabolites (DMs) for ‘14-1’ vs. ‘W600’ fruits, 394 for ‘14-1’ vs. F1, and 298 for F1 vs. ‘W600’ at the maturity stage. Combining WGNCA and metabolomic results, several DEGs and DMs were further identified as hub players in fruit quality formation: six DEGs with four DMs for flesh sweetness; six DEGs with 13 DMs for bitterness; nine DEGs with 10 DMs for flesh color; and nine DEGs with four DMs for flesh texture. Altogether, these observations not only expand our knowledge of the molecular basis of fruit quality in watermelon, but also provide potential targets for future watermelon improvement.

Metabolome and Transcriptome Analyses of Cucurbitacin Biosynthesis in Luffa (Luffa acutangula)

Cucurbitacins are extremely bitter compounds mainly present in Cucurbitaceae, where Luffa belongs. However, there is no comprehensive analysis of cucurbitacin biosynthesis in Luffa fruit. Therefore, this study analyzed bitter (WM709) and non-bitter (S1174) genotypes of Luffa to reveal the underlying mechanism of cucurbitacin biosynthesis by integrating metabolome and transcriptome analyses. A total of 422 metabolites were detected, including vitamins, essential amino acids, antioxidants, and antitumor substances. Of these, 131 metabolites showed significant differences between bitter (WM709) and non-bitter (S1174) Luffa fruits. The levels of isocucurbitacin B, cucurbitacin D, 23,24-dihydro cucurbitacin E, cucurbitacin F were significantly higher in bitter than in non-bitter Luffa. Transcriptome analysis showed that Bi, cytochromes P450s (CYP450s), and acyltransferase (ACT) of the cucurbitacin biosynthesis pathway, were significantly up-regulated. Moreover, drought stress and abscisic acid (ABA) activated genes of the cucurbitacin biosynthesis pathway. Furthermore, dual-luciferase reporter and yeast one-hybrid assays demonstrated that ABA-response element binding factor 1 (AREB1) binds to the Bi promoter to activate Bi expression. Comparative analysis of the Luffa and cucumber genomes showed that Bi, CYP450s, and ACT are located in the conserved syntenic loci, and formed a cucurbitacin biosynthesis cluster. This study provides important insights into major genes and metabolites of the cucurbitacin biosynthetic pathway, deepening the understanding of regulatory mechanisms of cucurbitacin biosynthesis in Luffa.

Integrated Metabolomics and Transcriptome Analyses Unveil Pathways Involved in Sugar Content and Rind Color of Two Sugarcane Varieties

Metabolic composition can have potential impact on several vital agronomic traits, and metabolomics, which represents the bioactive compounds in plant tissues, is widely considered as a powerful approach for linking phenotype-genotype interactions. However, metabolites related to cane traits such as sugar content, rind color, and texture differences in different sugarcane cultivars using metabolome integrated with transcriptome remain largely inconclusive. In this study, metabolome integrated with transcriptome analyses were performed to identify and quantify metabolites composition, and have better insight into the molecular mechanisms underpinning the different cane traits, namely, brix, rind color, and textures in the stems (S) and leaves (L) of sugarcane varieties FN41 and 165402. We also identified metabolites and associated genes in the phenylpropanoid and flavonoid biosynthesis pathways, starch and sucrose metabolism. A total of 512 metabolites from 11 classes, with the vast majority (122) belonging to flavonoids were identified. Moreover, the relatively high amount of D-fructose 6-p, D-glucose6-p and glucose1-p detected in FN41L may have been transported and distributed by source and sink of the cane, and a majority of them reached the stem of sugarcane FN41L, thereby promoting the high accumulation of sugar in FN41S. Observations also revealed that genes such as C4H, CHS, F3H, F3'H, DFR, and FG2 in phenylpropanoid and flavonoid biosynthesis pathways were the major factors impacting the rind color and contrasting texture of FN41 and 165204. Further analysis revealed that weighted gene co-expression network analysis (WGCNA) hub genes and six transcription factors, namely, Tify and NAC, MYB-related, C2C2-Dof, WRKY, and bHLH play a key role in phenylpropanoid biosynthesis, flavone and flavonol biosynthesis, starch and sucrose metabolism. Additionally, metabolites such as L-phenylalanine, tyrosine, sinapaldehyde, pinobanksin, kaempferin, and nictoflorin were the potential drivers of phenotypic differences. Our finding also demonstrated that genes and metabolites in the starch and sucrose metabolism had a significant effect on cane sugar content. Overall, this study provided valuable insight into the molecular mechanisms underpinning high sugar accumulation and rind color in sugarcane, which we believe is important for future sugarcane breeding programs and the selection of high biomass varieties.