Discovery of genes involved in anthocyanin biosynthesis from the rind and pith of three sugarcane varieties using integrated metabolic profiling and RNA-seq analysis

Transcriptome and metabolome analyses reveal molecular mechanisms of anthocyanin-related leaf color variation in poplar (Populus deltoides) cultivars

INTRODUCTION

Colored-leaf plants are increasingly popular for their aesthetic, ecological, and social value, which are important materials for research on the regulation of plant pigments. However, anthocyanin components and the molecular mechanisms of anthocyanin biosynthesis in colored-leaf poplar remain unclear. Consequently, an integrative analysis of transcriptome and metabolome is performed to identify the key metabolic pathways and key genes, which could contribute to the molecular mechanism of anthocyanin biosynthesis in the colored-leaf cultivars poplar.

METHODS

In this study, integrated metabolite and transcriptome analysis was performed to explore the anthocyanin composition and the specific regulatory network of anthocyanin biosynthesis in the purple leaves of the cultivars 'Quanhong' (QHP) and 'Zhongshanyuan' (ZSY). Correlation analysis between RNA-seq data and metabolite profiles were also performed to explore the candidate genes associated with anthocyanin biosynthesis. R2R3-MYB and bHLH TFs with differential expression levels were used to perform a correlation analysis with differentially accumulated anthocyanins.

RESULTS AND DISCUSSION

A total of 39 anthocyanin compounds were detected by LC-MS/MS analysis. Twelve cyanidins, seven pelargonidins, five delphinidins, and five procyanidins were identified as the major anthocyanin compounds, which were differentially accumulated in purple leaves of QHP and ZSY. The major genes associated with anthocyanin biosynthesis, including structural genes and transcription factors, were differentially expressed in purple leaves of QHP and ZSY through RNA-sequencing (RNA-seq) data analysis, which was consistent with quantitative real-time PCR analysis results. Correlation analysis between RNA-seq data and metabolite profiles showed that the expression patterns of certain differentially expressed genes in the anthocyanin biosynthesis pathway were strongly correlated with the differential accumulation of anthocyanins. One R2R3-MYB subfamily member in the SG5 subgroup, Podel.04G021100, showed a similar expression pattern to some structural genes. This gene was strongly correlated with 16 anthocyanin compounds, indicating that Podel.04G021100 might be involved in the regulation of anthocyanin biosynthesis. These results contribute to a systematic and comprehensive understanding of anthocyanin accumulation and to the molecular mechanisms of anthocyanin biosynthesis in QHP and ZSY.

Global Responses of Autopolyploid Sugarcane Badila (Saccharum officinarum L.) to Drought Stress Based on Comparative Transcriptome and Metabolome Profiling

Sugarcane (Saccharum spp. hybrid) is frequently affected by seasonal drought, which causes substantial declines in quality and yield. To understand the drought resistance mechanisms of S. officinarum, the main species of modern sugarcane, at a molecular level, we carried out a comparative analysis of transcriptome and metabolome profiling of the sugarcane variety Badila under drought stress (DS). Compared with control group (CG) plants, plants exposed to DS had 13,744 (6663 up-regulated and 7081 down-regulated) differentially expressed genes (DEGs). GO and KEGG analysis showed that the DEGs were enriched in photosynthesis-related pathways and most DEGs had down-regulated expression. Moreover, the chlorophyll content, photosynthesis (Photo), stomatal conductance (Cond), intercellular carbon dioxide concentration (Ci) and transpiration rate (Trmmol) were sharply decreased under DS. These results indicate that DS has a significant negative influence on photosynthesis in sugarcane. Metabolome analysis identified 166 (37 down-regulated and 129 up-regulated) significantly regulated metabolites (SRMs). Over 50% of SRMs were alkaloids, amino acids and their derivatives, and lipids. The five most significantly enriched KEGG pathways among SRMs were Aminoacyl-tRNA biosynthesis, 2-Oxocarboxylic acid metabolism, Biosynthesis of amino acids, Phenylalanine metabolism, and Arginine and proline metabolism (p < 0.05). Comparing CG with DS for transcriptome and metabolome profiling (T_CG/DS and M_CG/DS, respectively), we found three of the same KEGG-enriched pathways, namely Biosynthesis of amino acids, Phenylalanine metabolism and Arginine and proline metabolism. The potential importance of Phenylalanine metabolism and Arginine and proline metabolism was further analyzed for response to DS in sugarcane. Seven SRMs (five up-regulated and two down-regulated) and 60 DEGs (17 up-regulated and 43 down-regulated) were enriched in Phenylalanine metabolism under DS, of which novel.31261, Sspon.04G0008060-1A, Sspon.04G0008060-2B and Sspon.04G0008060-3C were significantly correlated with 7 SRMs. In Arginine and proline metabolism, eight SRMs (seven up-regulated and one down-regulated) and 63 DEGs (32 up-regulated and 31 down-regulated) were enriched, of which Sspon.01G0026110-1A (OAT) and Sspon.03G0002750-3D (P5CS) were strongly associated with proline (r > 0.99). These findings present the dynamic changes and possible molecular mechanisms of Phenylalanine metabolism as well as Arginine and proline metabolism under DS and provide a foundation for future research and sugarcane improvement.

Novel Insights into Anthocyanin Synthesis in the Calyx of Roselle Using Integrated Transcriptomic and Metabolomic Analyses

Roselle (Hibiscus sabdariffa L.) is an annual herbaceous plant of the genus Hibiscus in family Malvaceae. Roselle calyxes are rich in anthocyanins, which play important roles in human health. However, limited information is available on anthocyanin biosynthesis in the roselle calyx. In this study, transcriptomic and metabolomic analyses were performed to identify the key genes involved in anthocyanin biosynthesis in the roselle calyx. Three roselle cultivars with different calyx colors, including FZ-72 (red calyx, R), Baitao K (green calyx, G), and MG5 (stripped calyx, S), were used for metabolomic analyses with UPLC-Q-TOF/MS and RNA-seq. Forty-one compounds were quantified, including six flavonoids and 35 anthocyanins. The calyx of FZ-72 (red calyx) had the highest contents of anthocyanin derivatives such as delphinidin-3-O-sambubioside (955.11 μg/g) and cyanidin-3-O-sambubioside (531.37 μg/g), which were responsible for calyx color, followed by those in MG5 (stripped calyx) (851.97 and 330.06 μg/g, respectively). Baitao K (green calyx) had the lowest levels of these compounds. Furthermore, RNA-seq analysis revealed 114,415 differentially expressed genes (DEGs) in the calyxes at 30 days after flowering (DAF) for the corresponding cultivars FZ-72 (R), Baitao K (G), and MG5(S). The gene expression levels in the calyxes of the three cultivars were compared at different flowering stages, revealing 11,555, 11,949, and 7177 DEGs in R vs. G, R vs. S, and G vs. S, respectively. Phenylpropanoid and flavonoid biosynthesis pathways were found to be enriched. In the flavonoid pathway, 29, 28, and 27 genes were identified in G vs. R, G vs. S, and S vs. R, respectively. In the anthocyanin synthesis pathway, two, two, and one differential genes were identified in the three combinations; these differential genes belonged to the UFGT gene family. After joint analysis of the anthocyanin content in roselle calyxes, nine key genes belonging to the CHS, CHI, UFGT, FLS, ANR, DFR, CCoAOMT, SAT, and HST gene families were identified as strongly related to anthocyanin synthesis. These nine genes were verified using qRT-PCR, and the results were consistent with the transcriptome data. Overall, this study presents the first report on anthocyanin biosynthesis in roselle, laying a foundation for breeding roselle cultivars with high anthocyanin content.

Integrated Analysis of Transcriptome and Metabolome Reveals the Regulation of Chitooligosaccharide on Drought Tolerance in Sugarcane (Saccharum spp. Hybrid) under Drought Stress

Sugarcane (Saccharum spp. hybrid) is an important crop for sugar and biofuels, and often suffers from water shortages during growth. Currently, there is limited knowledge concerning the molecular mechanism involved in sugarcane response to drought stress (DS) and whether chitooligosaccharide could alleviate DS. Here, we carried out a combined transcriptome and metabolome of sugarcane in three different treatment groups: control group (CG), DS group, and DS + chitooligosaccharide group (COS). A total of 12,275 (6404 up-regulated and 5871 down-regulated) differentially expressed genes (DEGs) were identified when comparing the CG and DS transcriptomes (T_CG/DS), and 2525 (1261 up-regulated and 1264 down-regulated) DEGs were identified in comparing the DS and COS transcriptomes (T_DS/COS). GO and KEGG analysis showed that DEGs associated with photosynthesis were significantly enriched and had down-regulated expression. For T_DS/COS, photosynthesis DEGs were also significantly enriched but had up-regulated expression. Together, these results indicate that DS of sugarcane has a significantly negative influence on photosynthesis, and that COS can alleviate these negative effects. In metabolome analysis, lipids, others, amino acids and derivatives and alkaloids were the main significantly different metabolites (SDMs) observed in sugarcane response to DS, and COS treatment reduced the content of these metabolites. KEGG analysis of the metabolome showed that 2-oxocarboxylic acid metabolism, ABC transporters, biosynthesis of amino acids, glucosinolate biosynthesis and valine, leucine and isoleucine biosynthesis were the top-5 KEGG enriched pathways when comparing the CG and DS metabolome (M_CG/DS). Comparing DS with COS (M_DS/COS) showed that purine metabolism and phenylalanine metabolism were enriched. Combined transcriptome and metabolome analysis revealed that pyruvate and phenylalanine metabolism were KEGG-enriched pathways for CG/DS and DS/COS, respectively. For pyruvate metabolism, 87 DEGs (47 up-regulated and 40 down-regulated) and five SDMs (1 up-regulated and 4 down-regulated) were enriched. Pyruvate was closely related with 14 DEGs (|r| > 0.99) after Pearson’s correlation analysis, and only 1 DEG (Sspon.02G0043670-1B) was positively correlated. For phenylalanine metabolism, 13 DEGs (7 up-regulated and 6 down-regulated) and 6 SDMs (1 up-regulated and 5 down-regulated) were identified. Five PAL genes were closely related with 6 SDMs through Pearson’s correlation analysis, and the novel.31257 gene had significantly up-regulated expression. Collectively, our results showed that DS has significant adverse effects on the physiology, transcriptome, and metabolome of sugarcane, particularly genes involved in photosynthesis. We further show that COS treatment can alleviate these negative effects.

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.

Novel Insights into Anthocyanin Metabolism and Molecular Characterization of Associated Genes in Sugarcane Rinds Using the Metabolome and Transcriptome

Saccharum officinarum (sugarcane) is the fifth major cultivated crop around the world. Sugarcane rind is a promising source for anthocyanin pigments; however, limited information is available on the anthocyanin and its biosynthesis in sugarcane rinds. In this study, we have quantified 49 compounds including 6 flavonoids and 43 anthocyanins in the rind of 6 sugarcane cultivars by using LCMS/MS approach. Thirty of them were quantified for the first time in sugarcane. The 43 anthocyanins included 10 cyanidin (Cya), 11 pelargonidin (Pel), 9 peonidin (Peo), 5 malvidin (Mal), 4 delphinidin (Del), and 4 petunidin (Pet) metabolites. High contents of Cya derivatives were observed in the rind of YT71/210 (dark purple rind), such as cya-3-O-(6-O-malonyl)-glu 1283.3 µg/g and cya-3-O-glu 482.67 µg/g followed by ROC22 (red rind) 821.3 µg/g and 409 µg/g, respectively, whereas the YT93/159 (green rind) showed a minimum level of these compounds. Among six cultivars, ROC22 rind has high levels of Peo derivatives such as peo-3-O-glu (197 µg/g), peo-3-O-(6-O-malonyl)-glu (69 µg/g) and peo-3-O-(6-O-p-coumaryl)-glu (55.17 µg/g). The gene expression analysis revealed that some genes, including a MYB(t) gene, were highly associated with the color phenotype. Thus, we cloned and overexpressed the gene in Arabidopsis and found the pinkish brown color in the hypocotyl of all transgenic lines compared with the wild type. Hence, we have quantified a wide range of anthocyanins in major sugarcane cultivars, reported many new anthocyanins for the first time, and concluded that Cya and Peo derivatives are the major contributing factor of dissimilar colors in sugarcane. The finding and the verification of a novel MYB gene involved in anthocyanin biosynthesis have demonstrated that our study was very valuable for gene discovery and genetic improvement of sugarcane cultivars to harvest high anthocyanin contents.

The Transcriptome and Metabolome Reveal Stress Responses in Sulfur-Fumigated Cucumber (Cucumis sativus L.)

Sulfur (S) fumigation is a commonly used sterilization method in horticultural facilities against fungal diseases. S fumigation damaged cucumber leaves, although the response mechanism is unclear. This study analyzes the growth, transcriptome, and metabolomic profiles of young and mature leaves, ovaries, and commercial cucumber fruits to decipher the mechanism of cucumber stress response under S fumigation. S fumigation significantly changed the photosynthetic efficiency and reactive oxygen species (ROS) in leaves, but not fruit development, fruit mass, and peel color. Transcriptome analysis indicated that S fumigation strongly regulated stress defense genes. The weighted gene co-expression network analysis revealed that S fumigation regulated ASPG1, AMC1 defense genes, LECRK3, and PERK1 protein kinase. The abscisic acid (ABA)-mediated model of regulation under S fumigation was constructed. Metabolome analysis showed that S fumigation significantly upregulated or downregulated the contents of amino acids, organic acids, sugars, glycosides, and lipids (VIP > 1 and P-value < 0.05). The opposite Pearson's correlations of these differential metabolites implied that cucumber had different metabolic patterns in short-term and long-term S fumigation. Besides, the elevated levels of proline and triglyceride indicated that stress-responsive mechanisms existed in S-fumigated cucumber. Moreover, the comprehensive analysis indicated that S fumigation elevated secondary S-containing metabolites but decreased sulfate absorption and transportation in cucumber. Overall, our results provided a comprehensive assessment of S fumigation on cucumber, which laid the theoretical foundation for S fumigation in protected cultivation.