De novo Transcriptome Assembly of Chinese Kale and Global Expression Analysis of Genes Involved in Glucosinolate Metabolism in Multiple Tissues

Integrated metabolomic and transcriptomic analysis reveals the regulatory mechanisms of flavonoid and alkaloid biosynthesis in the new and old leaves of Murraya tetramera Huang

BACKGROUND

Murraya tetramera Huang is a traditional Chinese woody medicine. Its leaves contain flavonoids, alkaloids, and other active compounds, which have anti-inflammatory and analgesic effects, as well as hypoglycemic and lipid-lowering effects, and anti-tumor effects. There are significant differences in the content of flavonoids and alkaloids in leaves during different growth cycles, but the synthesis mechanism is still unclear.

RESULTS

In April 2021, new leaves (one month old) and old leaves (one and a half years old) of M. tetramera were used as experimental materials to systematically analyze the changes in differentially expressed genes (DEGs) and differentially accumulated metabolites (DAMs) with transcriptomics and metabolomics technology. This was done to identify the signaling pathways of flavonoid and alkaloid synthesis. The results showed that the contents of total alkaloids and flavonoids in old leaves were significantly higher than those in new leaves. Thirteen flavonoid compounds, three isoflavone compounds, and nineteen alkaloid compounds were identified, and 125 and 48 DEGs related to flavonoid and alkaloid synthesis were found, respectively. By constructing the KEGG (Kyoto Encyclopedia of Genes and Genomes) network of DEGs and DAMs, it was shown that the molecular mechanism of flavonoid biosynthesis in M. tetramera mainly focuses on the "flavonoid biosynthetic pathway" and the "flavonoid and flavonol biosynthetic pathway". Among them, p-Coumaryl alcohol, Sinapyl alcohol, Phloretin, and Isoquercitrin were significantly accumulated in old leaves, the up-regulated expression of CCR (cinnamoyl-CoA reductase) might promote the accumulation of p-Coumaryl alcohol, upregulation of F5H (ferulate-5-hydroxylase) might promote Sinapyl alcohol accumulation. Alkaloids, including indole alkaloids, pyridine alkaloids, imidazole alkaloids, and quinoline alkaloids, were significantly accumulated in old leaves, and a total of 29 genes were associated with these substances.

CONCLUSIONS

These data are helpful to better understand the biosynthesis of flavonoids and alkaloids in M. tetramera and provide a scientific basis for the development of medicinal components in M. tetramera.

Integrated metabolome and transcriptome analyses reveal the role of BoGSTF12 in anthocyanin accumulation in Chinese kale (Brassica oleracea var. alboglabra)

BACKGROUND

The vivid red, purple, and blue hues that are observed in a variety of plant fruits, flowers, and leaves are produced by anthocyanins, which are naturally occurring pigments produced by a series of biochemical processes occurring inside the plant cells. The purple-stalked Chinese kale, a popular vegetable that contains anthocyanins, has many health benefits but needs to be investigated further to identify the genes involved in the anthocyanin biosynthesis and translocation in this vegetable.

RESULTS

In this study, the purple- and green-stalked Chinese kale were examined using integrative transcriptome and metabolome analyses. The content of anthocyanins such as cyanidin-3-O-(6″-O-feruloyl) sophoroside-5-O-glucoside, cyanidin-3,5-O-diglucoside (cyanin), and cyanidin-3-O-(6″-O-p-hydroxybenzoyl) sophoroside-5-O-glucoside were considerably higher in purple-stalked Chinese kale than in its green-stalked relative. RNA-seq analysis indicated that 23 important anthocyanin biosynthesis genes, including 3 PAL, 2 C4H, 3 4CL, 3 CHS, 1 CHI, 1 F3H, 2 FLS, 2 F3'H, 1 DFR, 3 ANS, and 2 UFGT, along with the transcription factor BoMYB114, were significantly differentially expressed between the purple- and green-stalked varieties. Results of analyzing the expression levels of 11 genes involved in anthocyanin production using qRT-PCR further supported our findings. Association analysis between genes and metabolites revealed a strong correlation between BoGSTF12 and anthocyanin. We overexpressed BoGSTF12 in Arabidopsis thaliana tt19, an anthocyanin transport mutant, and this rescued the anthocyanin-loss phenotype in the stem and rosette leaves, indicating BoGSTF12 encodes an anthocyanin transporter that affects the accumulation of anthocyanins.

CONCLUSION

This work represents a key step forward in our understanding of the molecular processes underlying anthocyanin production in Chinese kale. Our comprehensive metabolomic and transcriptome analyses provide important insights into the regulatory system that controls anthocyanin production and transport, while providing a foundation for further research to elucidate the physiological importance of the metabolites found in this nutritionally significant vegetable.

Identification and in vitro enzymatic activity analysis of the AOP2 gene family associated with glucosinolate biosynthesis in Tumorous stem mustard (Brassica juncea var. tumida)

The major enzyme encoded by the glucosinolate biosynthetic gene AOP2 is involved in catalyzing the conversion of glucoiberin (GIB) into sinigrin (SIN) in Brassicaceae crops. The AOP2 proteins have previously been identified in several Brassicaceae species, but not in Tumorous stem mustard. As per this research, the five identified members of the AOP2 family from the whole genome of Brassica juncea named BjuAOP2.1-BjuAOP2.5 were found to be evenly distributed on five chromosomes. The subcellular localization results implied that BjuAOP2 proteins were mainly concentrated in the cytoplasm. Phylogenetic analysis of the AOP2 proteins from the sequenced Brassicaceae species in BRAD showed that BjuAOP2 genes were more closely linked to Brassica carinata and Brassica rapa than Arabidopsis. In comparison with other Brassicaceae plants, the BjuAOP2 members were conserved in terms of gene structures, protein sequences, and motifs. The light response and hormone response elements were included in the BjuAOP2 genes' cis-regulatory elements. The expression pattern of BjuAOP2 genes was influenced by the different stages of development and the type of tissue being examined. The BjuAOP2 proteins were used to perform the heterologous expression experiment. The results showed that all the five BjuAOP2 proteins can catalyze the conversion of GIB to SIN with different catalytic activity. These results provide the basis for further investigation of the functional study of BjuAOP2 in Tumorous stem mustard glucosinolate biosynthesis.

De novo transcriptome assembly and functional analysis reveal a dihydrochalcone 3-hydroxylase(DHC3H) of wild Malus species that produces sieboldin in vivo

Sieboldin is a specialised secondary metabolite of the group of dihydrochalcones (DHC), found in high concentrations only in some wild Malus species, closely related to the domesticated apple (Malus × domestica L.). To date, the first committed step towards the biosynthesis of sieboldin remains unknown. In this study, we combined transcriptomic analysis and a de novo transcriptome assembly to identify two putative 3-hydroxylases in two wild Malus species (Malus toringo (K. Koch) Carriere syn. sieboldii Rehder, Malus micromalus Makino) whose DHC profile is dominated by sieboldin. We assessed the in vivo activity of putative candidates to produce 3-hydroxyphloretin and sieboldin by de novo production in Saccharomyces cerevisiae. We found that CYP98A proteins of wild Malus accessions (CYP98A195, M. toringo and CYP98A196, M. micromalus) were able to produce 3-hydroxyphloretin, ultimately leading to sieboldin accumulation by co-expression with PGT2. CYP98A197-198 genes of M. × domestica, however, were unable to hydroxylate phloretin in vivo. CYP98A195-196 proteins exerting 3-hydroxylase activity co-localised with an endoplasmic reticulum marker. CYP98A protein model from wild accessions showed mutations in key residues close to the ligand pocket predicted using phloretin for protein docking modelling. These mutations are located within known substrate recognition sites of cytochrome P450s, which could explain the acceptance of phloretin in CYP98A protein of wild accessions. Screening a Malus germplasm collection by HRM marker analysis for CYP98A genes identified three clusters that correspond to the alleles of domesticated and wild species. Moreover, CYP98A isoforms identified in M. toringo and M. micromalus correlate with the accumulation of sieboldin in other wild and hybrid Malus genotypes. Taken together, we provide the first evidence of an enzyme producing sieboldin in vivo that could be involved in the key hydroxylation step towards the synthesis of sieboldin in Malus species.

BocODD1 and BocODD2 Regulate the Biosynthesis of Progoitrin Glucosinolate in Chinese Kale

Progoitrin (2-hydroxy-3-butenyl glucosinolate, PRO) is the main source of bitterness of Brassica plants. Research on the biosynthesis of PRO glucosinolate can aid the understanding of the nutritional value in Brassica plants. In this study, four ODD genes likely involved in PRO biosynthesis were cloned from Chinese kale. These four genes, designated as BocODD1-4, shared 75-82% similarities with the ODD sequence of Arabidopsis. The sequences of these four BocODDs were analyzed, and BocODD1 and BocODD2 were chosen for further study. The gene BocODD1,2 showed the highest expression levels in the roots, followed by the leaves, flowers, and stems, which is in accordance with the trend of the PRO content in the same tissues. Both the expression levels of BocODD1,2 and the content of PRO were significantly induced by high- and low-temperature treatments. The function of BocODDs involved in PRO biosynthesis was identified. Compared with the wild type, the content of PRO was increased twofold in the over-expressing BocODD1 or BocODD2 plants. Meanwhile, the content of PRO was decreased in the BocODD1 or BocODD2 RNAi lines more than twofold compared to the wildtype plants. These results suggested that BocODD1 and BocODD2 may play important roles in the biosynthesis of PRO glucosinolate in Chinese kale.

Transcriptome analysis of genes related to glucoraphanin and sulforaphane synthesis in methyl jasmonate treated broccoli (Brassica oleracea var. italica) hairy roots

Hairy roots obtained by infecting broccoli (Brassica oleracea var. italica) leaves with Agrobacterium rhizogenes (ATCC15834) have the characteristics of phytohormone autonomy, genetic stability and can produce a large amount of the anti-cancer substance Sulforaphane (SF) and the biosynthetic precursor Glucoraphanin (GRA). Under the induction of the exogenous signaling molecule methyl jasmonate (MeJA), the production of SF in broccoli hairy roots was significantly increased. However, the molecular mechanism of GRA and SF synthesis in hairy roots of broccoli treated with MeJA has not been reported. In this study, according to the yield of GRA and SF, the best concentration of MeJA treatment for hairy roots of broccoli was selected. After 18 days of growth, broccoli hairy roots were treated with 10 mmol L^-1 MeJA for 0, 3, 6, 9 and 12 h. Compared with 0 h, the yield of GRA and SF increased under other treatments. The highest yield of GRA and SF occurred at 9 h, which were 2.22-fold and 1.74-fold higher than those at 0 h. Brassica oleracea var. botrytis was used as reference genome, and 5,757 differentially expressed genes (DEG) were observed at 0, 3, 6, 9 and 12 h under 10 mmol L^-1 MeJA treatment, of which 4,673 were down-regulated and 1084 were up-regulated. The key genes regulating GRA synthesis, CYP79F1, CYP83A1, UGT74B1, FMOGS-OX5 and GSL-OH, were up-regulated at 0 and 3 h, and down-regulated the rest of the time; BCAT2 was up-regulated at 6, 9, 12 h, and at 0, 3 h expression was down-regulated, transcription factors MYB28 and MYB29 were down-regulated by exogenous MeJA treatment. A pathway of GRA biosynthesis and transformation pathways in MeJA-treated broccoli hairy roots was simulated and the molecular mechanism of GRA biosynthesis and SF accumulation in broccoli hairy roots under MeJA treatment was revealed.

Non-vernalization requirement in Chinese kale caused by loss of BoFLC and low expressions of its paralogs

We identified the loss of BoFLC gene as the cause of non-vernalization requirement in B. oleracea. Our developed codominant marker of BoFLC gene can be used for breeding program of B. oleracea crops. Many species of the Brassicaceae family, including some Brassica crops, require vernalization to avoid pre-winter flowering. Vernalization is an unfavorable trait for Chinese kale (Brassica oleracea var. chinensis Lei), a stem vegetable, and therefore it has been lost during its domestication/breeding process. To reveal the genetics of vernalization variation, we constructed an F₂ population through crossing a Chinese kale (a non-vernalization crop) with a kale (a vernalization crop). Using bulked segregant analysis (BSA) and RNA-seq, we identified one major quantitative trait locus (QTL) controlling vernalization and fine-mapped it to a region spanning 80 kb. Synteny analysis and PCR-based sequencing results revealed that compared to that of the kale parent, the candidate region of the Chinese kale parent lost a 9,325-bp fragment containing FLC homolog (BoFLC). In addition to the BoFLC gene, there are four other FLC homologs in the genome of B. oleracea, including Bo3g005470, Bo3g024250, Bo9g173370, and Bo9g173400. The qPCR analysis showed that the BoFLC had the highest expression among the five members of the FLC family. Considering the low expression levels of the four paralogs of BoFLC, we speculate that its paralogs cannot compensate the function of the lost BoFLC, therefore the presence/absence (PA) polymorphism of BoFLC determines the vernalization variation. Based on the PA polymorphism of BoFLC, we designed a codominant marker for the vernalization trait, which can be used for breeding programs of B. oleracea crops.

Regulation of Growth and Main Health-Promoting Compounds of Chinese Kale Baby-Leaf by UV-A and FR Light

Chinese kale baby leaves were hydroponically cultured under the basal light (Red: white LEDs = 2:3 at PPFD of 250 μmol·m^-2·s^-1) with different supplemental lighting, including individual ultraviolet-A (UV-A, 380 ± 10 nm, 20 μmol·m^-2·s^-1), far-red (FR, 735 ± 10 nm, 30 μmol·m^-2·s^-1) light, and their combination (UF) radiation in an artificial light plant factory. Effects of supplemental light qualities on morphology and physiology as well as health-promoting compounds of Chinese kale baby leaves were investigated. Application of UV-A and FR presented a positive effect on biomass, with a pronounced increase in petiole length, stem diameter, main stem length, and leaf area. Notably, plants under UF grew more vigorously than under other treatments. Higher levels of FRAP, vitamin C, total phenolic, and flavonoid were observed in plants under UV-A, while no striking changes or a decreasing trend recorded under FR and UF. Moreover, UV-A enhanced the glucosinolates (GLs) accumulation in Chinese kale baby leaves by increasing the predominant GLs (glucoraphanin and glucobrassicin) contents. RT-qPCR results indicated that UV-A upregulated the gene expressions of transcription factors and core structure genes related to GLs biosynthesis. However, downregulated or unchanged gene expressions of GLs biosynthesis-related genes in Chinese kale baby leaves were observed in FR and UF. Therefore, UV-A was benefited for the production of functional substances, while FR was conducive to a significant increase in crop yield. The combination of UV-A and FR, as a balance between yield and production of secondary metabolite, provided a new perspective for the application of artificial light in horticultural crop production.

Genome-Wide Identification and Expression Profiling of 2OGD Superfamily Genes from Three Brassica Plants

The 2-oxoglutarate and Fe(II)-dependent dioxygenase (2OGD) superfamily is the second largest enzyme family in the plant genome, and its members are involved in various oxygenation and hydroxylation reactions. Due to their important biochemical significance in metabolism, a systematic analysis of the plant 2OGD genes family is necessary. Here, we identified 160, 179, and 337 putative 2OGDs from Brassica rapa, Brassica oleracea, and Brassica napus. According to their gene structure, domain, phylogenetic features, function, and previous studies, we also divided 676 2OGDs into three subfamilies: DOXA, DOXB, and DOXC. Additionally, homologous and phylogenetic comparisons of three subfamily genes provided valuable insight into the evolutionary characteristics of the 2OGD genes from Brassica plants. Expression profiles derived from the transcriptome and Genevestigator database exhibited distinct expression patterns of the At2OGD, Br2OGD, and Bo2OGD genes in different developmental stages, tissues, or anatomical parts. Some 2OGD genes showed high expression levels in various tissues, such as callus, seed, silique, and root tissues, while other 2OGD genes were expressed at very low levels in other tissues. Analysis of six Bo2OGD genes in different tissues by qRT-PCR indicated that these genes are involved in the metabolism of gibberellin, which in turn regulates plant growth and development. Our working system analysed 2OGD gene families of three Brassica plants and laid the foundation for further study of their functional characterization.

Improvement of glucosinolates by metabolic engineering in Brassica crops

Glucosinolates (GSLs) are a class of sulfur- and nitrogen-containing, and amino acid-derived important secondary metabolites, which mainly present in plants of Brassicaceae family, including Brassica crops, such as broccoli, cabbage, and oilseed rape. The bioactive GSL metabolites confer benefits to plant defense, human health, and the unique flavor of some Brassica crops. However, certain GSL profiles have adverse effects and are known as anti-nutritional factors. This has attracted mounting attempts to increase beneficial GSLs and reduce detrimental ones in the most commonly consumed Brassica crops. We provide a comprehensive overview of metabolic engineering applied in Brassica crops to achieve this purpose, including modulation of GSL biosynthesis, ablation of GSL hydrolysis, inhibition of GSL transport processes, and redirection of metabolic flux to GSL. Moreover, advances in omics approaches, i.e., genomics, transcriptome, and metabolome, applied in the elucidation of GSL metabolism in Brassica crops, as well as promising and potential genome-editing technologies are also discussed.

Weighted gene co-expression network analysis identifies specific modules and hub genes related to coronary artery disease

This investigation seeks to dissect coronary artery disease molecular target candidates along with its underlying molecular mechanisms. Data on patients with CAD across three separate array data sets, GSE66360, GSE19339 and GSE97320 were extracted. The gene expression profiles were obtained by normalizing and removing the differences between the three data sets, and important modules linked to coronary heart disease were identified using weighted gene co-expression network analysis (WGCNA). Gene Ontology (GO) functional and Kyoto Encyclopedia of Genes and genomes (KEGG) pathway enrichment analyses were applied in order to identify statistically significant genetic modules with the Database for Annotation, Visualization and Integrated Discovery (DAVID) online tool (version 6.8; http://david.abcc.ncifcrf.gov ). The online STRING tool was used to construct a protein-protein interaction (PPI) network, followed by the use of Molecular Complex Detection (MCODE) plug-ins in Cytoscape software to identify hub genes. Two significant modules (green-yellow and magenta) were identified in the CAD samples. Genes in the magenta module were noted to be involved in inflammatory and immune-related pathways, based on GO and KEGG enrichment analyses. After the MCODE analysis, two different MCODE complexes were identified in the magenta module, and four hub genes (ITGAM, degree = 39; CAMP, degree = 37; TYROBP, degree = 28; ICAM1, degree = 18) were uncovered to be critical players in mediating CAD. Independent verification data as well as our RT-qPCR results were highly consistent with the above finding. ITGAM, CAMP, TYROBP and ICAM1 are potential targets in CAD. The underlying mechanism may be related to the transendothelial migration of leukocytes and the immune response.

Seasonal Variation in Transcriptomic Profiling of Tetrastigma hemsleyanum Fully Developed Tuberous Roots Enriches Candidate Genes in Essential Metabolic Pathways and Phytohormone Signaling

Tetrastigma hemsleyanum Diels et Gilg (Sanyeqing, SYQ) is a perennial climbing liana and an endemic plant to southern China. Its tuberous roots (TRs) are used in traditional Chinese medicine for treating some diseases such as high fever, pneumonia, asthma, hepatitis, and cancers. However, the mechanisms underlying the development of TR and the content of flavonoids and phenylpropanoids (FPs) are not well-understood. In this study, we performed a transcriptomic analysis of 12 fully developed TR (FD-TR) samples harvested in four seasons [spring (Sp), summer (Su), autumn (Au), and winter (Wi)] using the RNA-Sequencing (RNA-Seq). We obtained a total of 78.54 Gb raw data and 65,578 unigenes. Then, the unigenes were annotated by using six databases such as non-redundant protein database (NR), Pfam, eggNOG, SWISSProt, Kyoto Encyclopedia of Genes and Genomes (KEGG), and gene ontology (GO). The transcriptomic profiling showed closer relationships between the samples obtained in Su and Au than those obtained in Sp and Wi based on the results of both total unigenes and differentially expressed genes (DEGs). Three pathways, including the biosynthesis of FPs, metabolism of starch and sucrose, and signaling of phytohormones, were highly enriched, suggesting a gene-level seasonal variation. Based on the numbers of DEGs, brassinosteroid (BR) signal transduction factors appeared to play a key role in modulating the development of TRs while most of the auxin signaling genes were mainly activated in Wi and Sp FD-TRs. Most genes in the biosynthesis and biodegradation of starch and biodegradation of cellulose were activated in Wi FD-TRs. As determined by the high performance liquid chromatography (HPLC) and aluminum nitrate colorimetric method, the contents of total flavonoids and most detected FP components increased from Sp to Au but decreased in Wi. Enhanced expression levels of some genes in the biosynthetic pathways of FPs were detected in Su and Au samples, which corroborated well with metabolite content. Our findings provide the first transcriptomic and biochemical data on a seasonal variation in the composition of medically important metabolites in SYQ FD-TRs.

Effects of Supplementary Blue and UV-A LED Lights on Morphology and Phytochemicals of Brassicaceae Baby-Leaves

Brassicaceae baby-leaves are good source of functional phytochemicals. To investigate how Chinese kale and pak-choi baby-leaves in response to different wavebands of blue (430 nm and 465 nm) and UV-A (380 nm and 400 nm) LED, the plant growth, glucosinolates, antioxidants, and minerals were determined. Both agronomy traits and phytochemical contents were significantly affected. Blue and UV-A light played a predominant role in increasing the plant biomass and morphology, as well as the contents of antioxidant compounds (vitamin C, vitamin E, phenolics, and individual flavonols), the antioxidant activity (DPPH and FRAP), and the total glucosinolates accumulation. In particular, four light wavebands significantly decreased the content of progoitrin, while 400 nm UV-A light and 430 nm blue light were efficient in elevating the contents of sinigrin and glucobrassicin in Chinese kale. Meanwhile, 400 nm UV-A light was able to increase the contents of glucoraphanin, sinigrin, and glucobrassicin in pak-choi. From the global view of heatmap, blue lights were more efficient in increasing the yield and phytochemical levels of two baby-leaves.

Cloning and Expression Analysis of the BocMBF1c Gene Involved in Heat Tolerance in Chinese Kale

Chinese kale (Brassica oleracea var. chinensis Lei) is an important vegetable crop in South China, valued for its nutritional content and taste. Nonetheless, the thermal tolerance of Chinese kale still needs improvement. Molecular characterization of Chinese kale's heat stress response could provide a timely solution for developing a thermally tolerant Chinese kale variety. Here, we report the cloning of multi-protein bridging factor (MBF) 1c from Chinese kale (BocMBF1c), an ortholog to the key heat stress responsive gene MBF1c. Phylogenetic analysis showed that BocMBF1c is highly similar to the stress-response transcriptional coactivator MBF1c from Arabidopsis thaliana (AtMBF1c), and the BocMBF1c coding region conserves MBF1 and helix-turn-helix (HTH) domains. Moreover, the promoter region of BocMBF1c contains three heat shock elements (HSEs) and, thus, is highly responsive to heat treatment. This was verified in Nicotiana benthamiana leaf tissue using a green fluorescent protein (GFP) reporter. In addition, the expression of BocMBF1c can be induced by various abiotic stresses in Chinese kale which indicates the involvement of stress responses. The BocMBF1c-eGFP (enhanced green fluorescent protein) chimeric protein quickly translocated into the nucleus under high temperature treatment in Nicotiana benthamiana leaf tissue. Overexpression of BocMBF1c in Arabidopsis thaliana results in a larger size and enhanced thermal tolerance compared with the wild type. Our results provide valuable insight for the role of BocMBF1c during heat stress in Chinese kale.

Comprehensive dissection of transcript and metabolite shifts during seed germination and post-germination stages in poplar

BACKGROUND

Seed germination, a complex, physiological-morphogenetic process, is a critical stage in the life cycle of plants. Biological changes in germinating seeds have not been investigated in poplar, a model woody plant.

RESULTS

In this study, we exploited next-generation sequencing and metabolomics analysis and uncovered a series of significantly different genes and metabolites at various stages of seed germination and post germination. The K-means method was used to identify multiple transcription factors, including AP2/EREBP, DOF, and YABBY, involved in specific seed germination and post-germination stages. A weighted gene coexpression network analysis revealed that cell wall, amino acid metabolism, and transport-related pathways were significantly enriched during stages 3 and 5, with no significant enrichment observed in primary metabolic processes such as glycolysis and the tricarboxylic acid cycle. A metabolomics analysis detected significant changes in intermediate metabolites in these primary metabolic processes, while a targeted correlation network analysis identified the gene family members most relevant to these changing metabolites.

CONCLUSIONS

Taken together, our results provide important insights into the molecular networks underlying poplar seed germination and post-germination processes. The targeted correlation network analysis approach developed in this study can be applied to search for key candidate genes in specific biochemical reactions and represents a new strategy for joint multiomics analyses.

Comparative transcriptome analyses of genes involved in sulforaphane metabolism at different treatment in Chinese kale using full-length transcriptome sequencing

Background

Sulforaphane is a natural isothiocyanate available from cruciferous vegetables with multiple characteristics including antioxidant, antitumor and anti-inflammatory effect. Single-molecule real-time (SMRT) sequencing has been used for long-read de novo assembly of plant genome. Here, we investigated the molecular mechanism related to glucosinolates biosynthesis in Chinese kale using combined NGS and SMRT sequencing.

Results

SMRT sequencing produced 185,134 unigenes, higher than 129,325 in next-generation sequencing (NGS). NaCl (75 mM), methyl jasmonate (MeJA, 40 μM), selenate (Se, sodium selenite 100 μM), and brassinolide (BR, 1.5 μM) treatment induced 6893, 13,287, 13,659 and 11,041 differentially expressed genes (DEGs) in Chinese kale seedlings comparing with control. These genes were associated with pathways of glucosinolates biosynthesis, including phenylalanine, tyrosine and tryptophan biosynthesis, cysteine and methionine metabolism, and glucosinolate biosynthesis. We found NaCl decreased sulforaphane and glucosinolates (indolic and aliphatic) contents and downregulated expression of cytochrome P45083b1 (CYP83b1), S-alkyl-thiohydroximatelyase or carbon–sulfur lyase (SUR1) and UDP-glycosyltransferase 74B1 (UGT74b1). MeJA increased sulforaphane and glucosinolates contents and upregulated the expression of CYP83b1, SUR1 and UGT74b1; Se increased sulforaphane; BR increased expression of CYP83b1, SUR1 and UGT74b1, and increased glucosinolates contents. The desulfoglucosinolate sulfotransferases ST5a_b_c were decreased by all treatments.

Conclusions

We confirmed that NaCl inhibited the biosynthesis of both indolic and aliphatic glucosinolates, while MeJA and BR increased them. MeJA and BR treatments, conferred the biosynthesis of glucosinolates, and Se and MeJA contributed to sulforaphane in Chinese kale via regulating the expression of CYP83b1, SUR1 and UGT74b1.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5758-2) contains supplementary material, which is available to authorized users.

Molecular Cloning and Characterization of Three Glucosinolate Transporter (GTR) Genes from Chinese Kale

Chinese kale is a native vegetable in Southern China and the flowering stalk is the most commonly used edible part due to its high glucosinolate content and other nutritional qualities. The GTR protein played important roles in the glucosinolate transport process. In this study, three BocGTR1 genes were cloned from Chinese kale for the first time. Their gene structure, physicochemical properties, signal peptides, transmembrane structures, functional domains, second and third-order protein structures, and phylogenetic relationships were predicted. The expression levels of BocGTR1a and BocGTR1c were much higher than those of BocGTR1b in various tissues, especially in leaves and buds. In addition, the expression patterns of three genes were examined under various abiotic stresses or hormone treatment, including those induced by wounding, heat stress, methyl jasmonate, salicylic acid, salt, and MgCl₂ treatment. BocGTR1a and BocGTR1c were strongly induced by wounding and heat stress. The expression of BocGTR1a and BocGTR1c was significantly silenced in plants transformed by RNAi technology. Total glucosinolate content was significantly decreased in mature leaves and increased in roots of RNAi-transformed plants compared to wild-type plants. In addition, we found that BocGTR1a and BocGTR1c may participate in glucosinolate accumulation in different tissues with a selection for specific glucosinolates. These results indicated that BocGTR1a and BocGTR1c may be the key genes involved in the glucosinolate accumulation in different tissues of Chinese kale.

Gene co-expression network analysis reveals coordinated regulation of three characteristic secondary biosynthetic pathways in tea plant (Camellia sinensis)

BACKGROUND

The leaves of tea plants (Camellia sinensis) are used to produce tea, which is one of the most popular beverages consumed worldwide. The nutritional value and health benefits of tea are mainly related to three abundant characteristic metabolites; catechins, theanine and caffeine. Weighted gene co-expression network analysis (WGCNA) is a powerful system for investigating correlations between genes, identifying modules among highly correlated genes, and relating modules to phenotypic traits based on gene expression profiling. Currently, relatively little is known about the regulatory mechanisms and correlations between these three secondary metabolic pathways at the omics level in tea.

RESULTS

In this study, levels of the three secondary metabolites in ten different tissues of tea plants were determined, 87,319 high-quality unigenes were assembled, and 55,607 differentially expressed genes (DEGs) were identified by pairwise comparison. The resultant co-expression network included 35 co-expression modules, of which 20 modules were significantly associated with the biosynthesis of catechins, theanine and caffeine. Furthermore, we identified several hub genes related to these three metabolic pathways, and analysed their regulatory relationships using RNA-Seq data. The results showed that these hub genes are regulated by genes involved in all three metabolic pathways, and they regulate the biosynthesis of all three metabolites. It is notable that light was identified as an important regulator for the biosynthesis of catechins.

CONCLUSION

Our integrated omics-level WGCNA analysis provides novel insights into the potential regulatory mechanisms of catechins, theanine and caffeine metabolism, and the identified hub genes provide an important reference for further research on the molecular biology of tea plants.

Evaluation of the Nutritional Quality of Chinese Kale (Brassica alboglabra Bailey) Using UHPLC-Quadrupole-Orbitrap MS/MS-Based Metabolomics

Chinese kale (Brassica alboglabra Bailey) is a widely consumed vegetable which is rich in antioxidants and anticarcinogenic compounds. Herein, we used an untargeted ultra-high-performance liquid chromatography (UHPLC)-Quadrupole-Orbitrap MS/MS-based metabolomics strategy to study the nutrient profiles of Chinese kale. Seven Chinese kale cultivars and three different edible parts were evaluated, and amino acids, sugars, organic acids, glucosinolates and phenolic compounds were analysed simultaneously. We found that two cultivars, a purple-stem cultivar W1 and a yellow-flower cultivar Y1, had more health-promoting compounds than others. The multivariate statistical analysis results showed that gluconapin was the most important contributor for discriminating both cultivars and edible parts. The purple-stem cultivar W1 had higher levels of some phenolic acids and flavonoids than the green stem cultivars. Compared to stems and leaves, the inflorescences contained more amino acids, glucosinolates and most of the phenolic acids. Meanwhile, the stems had the least amounts of phenolic compounds among the organs tested. Metabolomics is a powerful approach for the comprehensive understanding of vegetable nutritional quality. The results provide the basis for future metabolomics-guided breeding and nutritional quality improvement.