Metabolomic and transcriptomic profiling reveals distinct nutritional properties of cassavas with different flesh colors

Transcriptomic study of the role of MeFtsZ2-1 in pigment accumulation in cassava leaves

MeFtsZ2-1 is a key gene for plant plastid division, but the mechanism by which MeFtsZ2-1 affects pigment accumulation in cassava (Manihot esculenta Crantz) through plastids remains unclear. We found that MeFtsZ2-1 overexpression in cassava (OE) exhibited darker colors of leaves, with increased levels of anthocyanins and carotenoids. Further observation via Transmission Electron Microscopy (TEM) revealed no apparent defects in chloroplast structure but an increase in the number of plastoglobule in OE leaves. RNA-seq results showed 1582 differentially expressed genes (DEGs) in leaves of OE. KEGG pathway analysis indicated that these DEGs were enriched in pathways related to flavonoid, anthocyanin, and carotenoid biosynthesis. This study reveals the role of MeFtsZ2-1 in cassava pigment accumulation from a physiological and transcriptomic perspective, providing a theoretical basis for improving cassava quality.

Metabolic GWAS-based dissection of genetic basis underlying nutrient quality variation and domestication of cassava storage root

BACKGROUND

Metabolites play critical roles in regulating nutritional qualities of plants, thereby influencing their consumption and human health. However, the genetic basis underlying the metabolite-based nutrient quality and domestication of root and tuber crops remain largely unknown.

RESULTS

We report a comprehensive study combining metabolic and phenotypic genome-wide association studies to dissect the genetic basis of metabolites in the storage root (SR) of cassava. We quantify 2,980 metabolic features in 299 cultivated cassava accessions. We detect 18,218 significant marker-metabolite associations via metabolic genome-wide association mapping and identify 12 candidate genes responsible for the levels of metabolites that are of potential nutritional importance. Me3GT, MeMYB4, and UGT85K4/UGT85K5, which are involved in flavone, anthocyanin, and cyanogenic glucoside metabolism, respectively, are functionally validated through in vitro enzyme assays and in vivo gene silencing analyses. We identify a cluster of cyanogenic glucoside biosynthesis genes, among which CYP79D1, CYP71E7b, and UGT85K5 are highly co-expressed and their allelic combination contributes to low linamarin content. We find MeMYB4 is responsible for variations in cyanidin 3-O-glucoside and delphinidin 3-O-rutinoside contents, thus controlling SR endothelium color. We find human selection affects quercetin 3-O-glucoside content and SR weight per plant. The candidate gene MeFLS1 is subject to selection during cassava domestication, leading to decreased quercetin 3-O-glucoside content and thus increased SR weight per plant.

CONCLUSIONS

These findings reveal the genetic basis of cassava SR metabolome variation, establish a linkage between metabolites and agronomic traits, and offer useful resources for genetically improving the nutrition of cassava and other root crops.

Metabolite profiling and transcriptomic analyses demonstrate the effects of biocontrol agents on alkaloid accumulation in Fritillaria thunbergii

BACKGROUND

During Fritillaria thunbergii planting, pests and diseases usually invade the plant, resulting in reduced yield and quality. Previous studies have demonstrated that using biocontrol agents can effectively control grubs and affect the steroid alkaloids content in F. thunbergii. However, the molecular regulatory mechanisms underlying the differences in the accumulation of steroid alkaloids in response to biocontrol agents remain unclear.

RESULTS

Combined transcriptomic and metabolic analyses were performed by treating the bulbs of F. thunbergii treated with biocontrol agents during planting. Otherwise, 48 alkaloids including 32 steroid alkaloids, 6 indole alkaloids, 2 scopolamine-type alkaloids, 1 isoquinoline alkaloid, 1 furoquinoline alkaloid, and 6 other alkaloids were identified. The content of steroidal alkaloids particularly peimine, peiminine, and veratramine, increased significantly in the group treated with the biocontrol agents. Transcriptome sequencing identified 929 differential genes using biocontrol agents, including 589 upregulated and 340 downregulated genes. Putative biosynthesis networks of steroid alkaloids have been established and combined with differentially expressed structural unigenes, such as acetyl-CoA C-acetyl-transferase, acelyl-CoAC-acetyltransferase3-hydroxy-3-methylglutaryl-coenzyme A synthase, 1-deoxy-D-xylulose-5-phosphate reductor-isomerase, 2-C-methyl-D-erythritol-4-phosphate cytidylyltransferase and 4-hydroxy-3-methylbut-2-enyl diphosphate reductase. In addition, biological processes such as amino acid accumulation and oxidative phosphorylation were predicted to be related to the synthesis of steroid alkaloids. Cytochrome P450 enzymes also play crucial roles in the steroid alkaloid synthesis. The transcription factor families MYB and bHLH were significantly upregulated after using biocontrol agents.

CONCLUSIONS

Biocontrol agents increased the steroid alkaloids accumulation of steroid alkaloids by affecting key enzymes in the steroid alkaloid synthesis pathway, biological processes of oxidative phosphorylation and amino acid synthesis, cytochrome P450 enzymes, and transcription factors. This study revealed the mechanism underlying the difference in steroidal alkaloids in F. thunbergii after using biocontrol agents, laying the groundwork for future industrial production of steroid alkaloids and ecological planting of medicinal materials in the future.

Genomic-wide identification and expression analysis of R2R3-MYB transcription factors related to flavonol biosynthesis in Morinda officinalis

BACKGROUND

The R2R3-MYB transcription factors are a crucial and extensive gene family in plants, which participate in diverse processes, including development, metabolism, defense, differentiation, and stress response. In the Lingnan region of China, Morinda officinalis is extensively grown and is renowned for its use as both a medicinal herb and food source. However, there are relatively few reports on the R2R3-MYB transcription factor family in M.officinalis.

RESULTS

In this study, we identified 97 R2R3-MYB genes in the genome of Morinda officinalis and classified them into 32 subgroups based on phylogenetic comparison with Arabidopsis thaliana. The lack of recent whole-genome duplication events in M.officinalis may be the reason for the relatively few members of the R2R3-MYB family. We also further analyzed the physical and chemical characteristics, conserved motifs, gene structure, and chromosomal location. Gene duplication events found 21 fragment duplication pairs and five tandem duplication event R2R3-MYB genes in M.officinalis may also affect gene family expansion. Based on phylogenetic analysis, cis-element analysis, co-expression analysis and RT-qPCR, we concluded that MoMYB33 might modulate flavonol levels by regulating the expression of 4-coumarate-CoA ligase Mo4CL2, chalcone isomerase MoCHI3, and flavonol synthase MoFLS4/11/12. MoMYB33 and AtMYB111 showed the highest similarity of 79% and may be involved in flavonol synthase networks by the STRING database. Moreover, we also identified MoMYB genes that respond to methyl Jasmonate (MeJA) and abscisic acid (ABA) stress by RT-qPCR.

CONCLUSIONS

This study offers a thorough comprehension of R2R3-MYB in M.officinalis, which lays the foundation for the regulation of flavonol synthesis and the response of MoMYB genes to phytohormones in M.officinalis.

Integrative analysis of metabolome and transcriptome reveals the mechanism of color formation in cassava (Manihot esculenta Crantz) leaves

Cassava (Manihot esculenta Crantz) leaves are often used as vegetables in Africa. Anthocyanins possess antioxidant, anti-inflammatory, anti-cancer, and other biological activities. They are poor in green leaves but rich in the purple leaves of cassava. The mechanism of anthocyanin's accumulation in cassava is poorly understood. In this study, two cassava varieties, SC9 with green leaves and Ziyehuangxin with purple leaves (PL), were selected to perform an integrative analysis using metabolomics and transcriptomics. The metabolomic analysis indicated that the most significantly differential metabolites (SDMs) belong to anthocyanins and are highly accumulated in PL. The transcriptomic analysis revealed that differentially expressed genes (DEGs) are enriched in secondary metabolites biosynthesis. The analysis of the combination of metabolomics and transcriptomics showed that metabolite changes are associated with the gene expressions in the anthocyanin biosynthesis pathway. In addition, some transcription factors (TFs) may be involved in anthocyanin biosynthesis. To further investigate the correlation between anthocyanin accumulation and color formation in cassava leaves, the virus-induced gene silencing (VIGS) system was used. VIGS-MeANR silenced plant showed the altered phenotypes of cassava leaves, partially from green to purple color, resulting in a significant increase of the total anthocyanin content and reduction in the expression of MeANR. These results provide a theoretical basis for breeding cassava varieties with anthocyanin-rich leaves.

Integrative Transcriptomic and Metabolomic Analysis Reveals the Molecular Mechanism of Red Maple (Acer rubrum L.) Leaf Coloring

This study employed a combination of ultraviolet spectrophotometry, LC-ESI-MS/MS system, and RNA-sequencing technology; the extracts and isolation of total RNA from the red and yellow leaf strains of red maple (Acer rubrum L.) at different developmental stages were subjected to an intercomparison of the dynamic content of chlorophyll and total anthocyanin, flavonoid metabolite fingerprinting, and gene expression. The metabonomic results indicated that one hundred and ninety-two flavonoids were identified, which could be classified into eight categories in the red maple leaves. Among them, 39% and 19% were flavones and flavonols, respectively. The metabolomic analysis identified 23, 32, 24, 24, 38, and 41 DAMs in the AR1018r vs. AR1031r comparison, the AR1018r vs. AR1119r comparison, the AR1031r vs. AR1119r comparison, the AR1018y vs. AR1031y comparison, the AR1018y vs. AR1119y comparison, and the AR1031y vs. AR1119y comparison, respectively. In total, 6003 and 8888 DEGs were identified in AR1018r vs. AR1031r comparison and in the AR1018y vs. AR1031y comparison, respectively. The GO and KEGG analyses showed that the DEGs were mainly involved in plant hormone signal transduction, flavonoid biosynthesis, and other metabolite metabolic processes. The comprehensive analysis revealed that caffeoyl-CoA 3-O-methyltransferase (Cluster-28704.45358 and Cluster-28704.50421) was up-regulated in the red strain but down-regulated in the yellow strain, while Peonidin 3-O-glucoside chloride and Pelargonidin 3-O-beta-D-glucoside were up-regulated in both the red and yellow strains. By successfully integrating the analyses on the behavior of pigment accumulation, dynamics of flavonoids, and differentially expressed genes with omics tools, the regulation mechanisms underlying leaf coloring in red maple at the transcriptomic and metabolomic levels were demonstrated, and the results provide valuable information for further research on gene function in red maple.

Metabolomics analysis unveils important changes involved in the salt tolerance of Salicornia europaea

Salicornia europaea is one of the world's salt-tolerant plant species and is recognized as a model plant for studying the metabolism and molecular mechanisms of halophytes under salinity. To investigate the metabolic responses to salinity stress in S. europaea, this study performed a widely targeted metabolomic analysis after analyzing the physiological characteristics of plants exposed to various NaCl treatments. S. europaea exhibited excellent salt tolerance and could withstand extremely high NaCl concentrations, while lower NaCl conditions (50 and 100 mM) significantly promoted growth by increasing tissue succulence and maintaining a relatively stable K⁺ concentration. A total of 552 metabolites were detected, 500 of which were differently accumulated, mainly consisting of lipids, organic acids, saccharides, alcohols, amino acids, flavonoids, phenolic acids, and alkaloids. Sucrose, glucose, p-proline, quercetin and its derivatives, and kaempferol derivatives represented core metabolites that are responsive to salinity stress. Glycolysis, flavone and flavonol biosynthesis, and phenylpropanoid biosynthesis were considered as the most important pathways responsible for salt stress response by increasing the osmotic tolerance and antioxidant activities. The high accumulation of some saccharides, flavonoids, and phenolic acids under 50 mM NaCl compared with 300 mM NaCl might contribute to the improved salt tolerance under the 50 mM NaCl treatment. Furthermore, quercetin, quercetin derivatives, and kaempferol derivatives showed varied change patterns in the roots and shoots, while coumaric, caffeic, and ferulic acids increased significantly in the roots, implying that the coping strategies in the shoots and roots varied under salinity stress. These findings lay the foundation for further analysis of the mechanism underlying the response of S. europaea to salinity.

Transcriptomic and metabolic analyses reveal the potential mechanism of increasing steroidal alkaloids in Fritillaria hupehensis through intercropping with Magnolia officinalis

Fritillaria hupehensis, a well-known medicinal perennial herb, is used as an antitussive and an expectorant. Continuous cropping and monoculture cultivation usually negativly affect the growth of F. hupehensis. Compared with the monoculture system, the F. hupehensis-Magnolia officinalis intercropping system significantly increases the yield of F. hupehensis. However, changes in steroidal alkaloid metabolites (the most important bioactive components) and their molecular regulatory mechanisms in F. hupehensis intercropping system remain unclear. We performed comparative transcriptomic and metabolomic analyses of F. hupehensis bulbs grown in monocropping and intercropping systems. A total of 40 alkaloids were identified, including 26 steroidal alkaloids, 4 plumeranes, 3 phenolamines, 1 pyridine alkaloid, and 6 other alkaloids. The results showed that intercropping significantly increased the levels of peimine, peiminine, hupehenine, korseveridine, verticinone N-oxide, delafrine, tortifoline, pingbeinone, puqienine B, puqienine E, jervine, ussuriedine, hydroxymandelonitrile, N-feruloylputrescine, and N-benzylmethylene isomethylamine in F. hupehensis, but decreased the levels of indole, p-coumaroylputrescine, and N-benzylformamide. Transcriptome sequencing identified 11,466 differentially expressed unigenes in F. hupehensis under the intercropping system, of which 5,656 genes were up-regulated and 5,810 genes were down-regulated. We proposed a possible steroidal alkaloid biosynthesis pathway, in which 12 differentially expressed genes were identified. The higher expressions of these genes in the intercropping system positively correlated with the high accumulation of peimine, peiminine, and hupehenine, further validating our proposal. Moreover, the biological processes of oxidative phosphorylation and plant hormone signal transduction, cytochrome P450 enzymes, ATP-binding cassette transporters, and transcription factors may play pivotal roles in the regulation of steroidal alkaloid biosynthesis. This study revealed the underlying molecular mechanisms of intercropping in improving steroidal alkaloids in F. hupehensis at the transcriptome and metabolome levels. These findings provided a theoretical foundation for sustainable development of this ecological planting method.

Identification of Key Gene Networks Controlling Soluble Sugar and Organic Acid Metabolism During Oriental Melon Fruit Development by Integrated Analysis of Metabolic and Transcriptomic Analyses

Oriental melon (Cucumis melo var. acidulus) is one of the most economically important fruit crops worldwide. To elucidate the molecular basis related to soluble sugar and organic acid metabolism in the fruits of two oriental melon cultivars with different sweetness, we performed integrated metabolomic and transcriptomic analyses of the fruits of 'Tianbao' (A) with high sweetness and 'Xiaocuigua' (B) with low sweetness at different ripening stages. The high accumulation of sucrose, D-glucose, D-(+)-raffinose, and the relatively lower citric acid and malic acid might contribute to the sweet taste of A. By screening the differentially expressed genes (DEGs) and correlation analysis of the DEGs and differentially accumulated metabolites, we deduced that the B cultivar might promote the conversion of glucose and fructose into intermediate compounds for downstream processes such as glycolysis. The tricarboxylic acid (TCA) cycle might also be enhanced compared to A, thus resulting in the differential accumulation of soluble sugars and organic acids, ultimately causing the taste difference between the two oriental melon cultivars. Our finding provides important information for further exploring the metabolic mechanisms of soluble sugars and organic acids in oriental melon.

Integrated Metabolomic and Transcriptomic Analyses Reveal Novel Insights of Anthocyanin Biosynthesis on Color Formation in Cassava Tuberous Roots

Yellow roots are of higher nutritional quality and better appearance than white roots in cassava, a crucial tropical and subtropical root crop. In this work, two varieties with yellow and white cassava roots were selected to explore the mechanisms of color formation by using comparative metabolome and transcriptome analyses during seven developmental stages. Compared with the white-rooted cassava, anthocyanins, catechin derivatives, coumarin derivatives, and phenolic acids accumulated at higher levels in yellow-rooted cassava. Anthocyanins were particularly enriched and displayed different accumulation patterns during tuberous root development. This was confirmed by metabolic comparisons between five yellow-rooted and five white-rooted cassava accessions. The integrative metabolomic and transcriptomic analysis further revealed a coordinate regulation of 16 metabolites and 11 co-expression genes participating in anthocyanin biosynthesis, suggesting a vital role of anthocyanin biosynthesis in yellow pigmentation in cassava tuberous roots. In addition, two transcriptional factors, i.e., MeMYB5 and MeMYB42, were also identified to co-express with these anthocyanin biosynthesis genes. These findings expand our knowledge on the role of anthocyanin biosynthesis in cassava root color formation, and offer useful information for the genetic breeding of yellow-rooted cassava in the future.

Comparison of Metabolome and Transcriptome of Flavonoid Biosynthesis in Two Colors of Coreopsis tinctoria Nutt

Coreopsis tinctoria Nutt. (C. tinctoria) has a long history of application and high economic and medicinal value. Flavonoids, the main active components of C. tinctoria, are widely studied in pharmacology and food development. However, the flavonoid biosynthesis pathway in C. tinctoria is unclear. In this study, we comprehensively compared the transcriptomes and metabolite profiles of two colors of C. tinctoria flowers (LS and JS) at different flowering stages. A total of 165 flavonoids (46 flavonoids, 42 flavonols, 22 anthocyanins, 18 chalcones, 12 dihydroflavonols, nine isoflavones, eight dihydroflavonoids, six flavanols, and two tannins) were identified in LS and JS at different flowering stages. Thirty-three metabolites (11 anthocyanins, 11 flavonols, seven flavonoids, two dihydroflavonols, one dihydroflavone, and one chalcone) were found to be statistically significantly different in the LS vs. JS groups. LS flowers accumulated higher levels of 10 anthocyanins (seven cyanidins and three pelargonidins) than JS flowers. Furthermore, candidate genes related to the regulation of flavonoid and anthocyanin synthesis were identified and included 28 structural genes (especially F3H, Cluster-28756.299649, and 3GT, Cluster-28756.230942) in LS and JS, six key differentially expressed transcription factors (especially MYB90a, Cluster-28756.143139) in LS and JS, and 17 other regulators (mainly including transporter proteins and others) in LS. Our results provide valuable information for further studies on the mechanism underlying flavonoid biosynthesis in C. tinctoria.

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.