Metabolite Profiling and Transcriptome Analyses Provide Insights into the Flavonoid Biosynthesis in the Developing Seed of Tartary Buckwheat (Fagopyrum tataricum)

Buckwheat OMICS: present status and future prospects

Buckwheat (Fagopyrum spp.) is an underutilized resilient crop of North Western Himalayas belonging to the family Polygonaceae and is a source of essential nutrients and therapeutics. Common Buckwheat and Tatary Buckwheat are the two main cultivated species used as food. It is the only grain crop possessing rutin, an important metabolite with high nutraceutical potential. Due to its inherent tolerance to various biotic and abiotic stresses and a short life cycle, Buckwheat has been proposed as a model crop plant. Nutritional security is one of the major concerns, breeding for a nutrient-dense crop such as Buckwheat will provide a sustainable solution. Efforts toward improving Buckwheat for nutrition and yield are limited due to the lack of available: genetic resources, genomics, transcriptomics and metabolomics. In order to harness the agricultural importance of Buckwheat, an integrated breeding and OMICS platforms needs to be established that can pave the way for a better understanding of crop biology and developing commercial varieties. This, coupled with the availability of the genome sequences of both Buckwheat species in the public domain, should facilitate the identification of alleles/QTLs and candidate genes. There is a need to further our understanding of the molecular basis of the genetic regulation that controls various economically important traits. The present review focuses on: the food and nutritional importance of Buckwheat, its various omics resources, utilization of omics approaches in understanding Buckwheat biology and, finally, how an integrated platform of breeding and omics will help in developing commercially high yielding nutrient rich cultivars in Buckwheat.

A decade of advances in the study of buckwheat for organic farming and agroecology (2013-2023)

During the last decade, research has shown the environment and human health benefits of growing buckwheat (Fagopyrum spp.). This comprehensive review aims to summarize the major advancements made in the study of buckwheat from 2013 to 2023, focusing on its agronomic characteristics, nutritional value, and potential applications in sustainable agriculture. The review examines the diverse applications of buckwheat in organic and agroecological farming systems, and discusses the ability of buckwheat to control weeds through allelopathy, competition, and other sustainable farming methods, such as crop rotation, intercropping and green manure, while improving soil health and biodiversity. The review also explores the nutritional value of buckwheat. It delves into the composition of buckwheat grains, emphasizing their high protein content, and the presence of essential amino acids and valuable micronutrients, which is linked to health benefits such as lowering cholesterol levels, controlling diabetes and acting against different types of cancer, among others. Finally, the review concludes by highlighting the gaps in current knowledge, and proposing future research directions to further optimize buckwheat production in organic or agroecological farming systems. It emphasizes the need for interdisciplinary collaboration, and the integration of traditional knowledge with modern scientific approaches to unlock the full potential of buckwheat as a sustainable crop.

Tartary buckwheat rutin: Accumulation, metabolic pathways, regulation mechanisms, and biofortification strategies

Rutin is a significant flavonoid with strong antioxidant property and various therapeutic effects. It plays a crucial role in disease prevention and human health maintenance, especially in anti-inflammatory, antidiabetic, hepatoprotective and cardiovascular effects. While many plants can synthesize and accumulate rutin, tartary buckwheat is the only food crop possessing high levels of rutin. At present, the rutin content (RC) is regarded as the key index for evaluating the nutritional quality of tartary buckwheat. Consequently, rutin has become the focus for tartary buckwheat breeders and has made considerable progress. Here, we summarize research on the rutin in tartary buckwheat in the past two decades, including its accumulation, biosynthesis and breakdown pathways, and regulatory mechanisms. Furthermore, we propose several strategies to increase the RC in tartary buckwheat seeds based on current knowledge. This review aims to provide valuable references for elevating the quality of tartary buckwheat in the future.

Metabolic profiling provides insights into the accumulation patterns of flavonoids and phenolic acids in different parts of Lactuca indica L

Lactuca indica L. contains high amounts of flavonoids and phenolic acids. However, there is limited information on the composition of these compounds in different parts of the plant. The present study analyzed the secondary metabolite profiles of the stem, leaf, flower, and seed of Lactuca indicaL.cv. Mengzao (LIM) using a widely targeted metabolomic approach. A total of 576 secondary metabolites were identified, including 218, 267, 232, 286, 302, and 308 differentially accumulated metabolites (DAMs) in the stem_vs_leaf, stem_vs_flower, stem_vs_seed, leaf_vs_flower, leaf_vs_seed, and flower_vs_seed comparisons. In particular, considerable differences were detected in the flavonoids and phenolic acids, five flavonoids, five phenolic acids, one triterpenoid and one alkaloid being differentially accumulated in the four parts. Compared to the stem and flower, the leaf and seed had higher total flavonoid content and total phenolic content. Thses findings provide comprehensive insights into utilizing different parts of LIM in developing functional food products.

Germplasm Resources and Metabolite Marker Screening of High-Flavonoid Tartary Buckwheat (Fagopyrum tataricum)

Tartary buckwheat is an annual minor cereal crop with a variety of secondary metabolites, endowing it with a high nutritional and medicinal value. Flavonoids constitute the primary compounds of Tartary buckwheat. Recently, metabolomics, as an adjunct breeding method, has been increasingly employed in crop research. This study explores the correlation between the total flavonoid content (TFC) and antioxidant capacity in 167 Tartary buckwheat varieties. Ten Tartary buckwheat varieties with significant differences in flavonoid content and antioxidant capacity were selected by cluster analysis. With the use of liquid chromatography-mass spectrometry, 58 flavonoid compounds were identified, namely, 42 flavonols, 10 flavanols, 3 flavanones, 1 isoflavone, 1 anthocyanidin, and 1 proanthocyanidin. Different samples were clearly separated by employing principal component analysis and partial least-squares discriminant analysis. Eight differential flavonoid compounds were further selected through volcano plots and variable importance in projection. Differential metabolites were highly correlated with TFC and antioxidant capacity. Finally, metabolic markers of kaempferol-3-O-hexoside, kaempferol-7-O-glucoside, and naringenin-O-hexoside were determined by the random forest model. The findings provide a basis for the selection and identification of Tartary buckwheat varieties with high flavonoid content and strong antioxidant activity.

Tartary Buckwheat (Fagopyrum tataricum) FtTT8 Inhibits Anthocyanin Biosynthesis and Promotes Proanthocyanidin Biosynthesis

Tartary buckwheat (Fagopyrum tataricum) is an important plant, utilized for both medicine and food. It has become a current research hotspot due to its rich content of flavonoids, which are beneficial for human health. Anthocyanins (ATs) and proanthocyanidins (PAs) are the two main kinds of flavonoid compounds in Tartary buckwheat, which participate in the pigmentation of some tissue as well as rendering resistance to many biotic and abiotic stresses. Additionally, Tartary buckwheat anthocyanins and PAs have many health benefits for humans and the plant itself. However, little is known about the regulation mechanism of the biosynthesis of anthocyanin and PA in Tartary buckwheat. In the present study, a bHLH transcription factor (TF) FtTT8 was characterized to be homologous with AtTT8 and phylogenetically close to bHLH proteins from other plant species. Subcellular location and yeast two-hybrid assays suggested that FtTT8 locates in the nucleus and plays a role as a transcription factor. Complementation analysis in Arabidopsis tt8 mutant showed that FtTT8 could not recover anthocyanin deficiency but could promote PAs accumulation. Overexpression of FtTT8 in red-flowering tobacco showed that FtTT8 inhibits anthocyanin biosynthesis and accelerates proanthocyanidin biosynthesis. QRT-PCR and yeast one-hybrid assay revealed that FtTT8 might bind to the promoter of NtUFGT and suppress its expression, while binding to the promoter of NtLAR and upregulating its expression in K326 tobacco. This displayed the bidirectional regulating function of FtTT8 that negatively regulates anthocyanin biosynthesis and positively regulates proanthocyanidin biosynthesis. The results provide new insights on TT8 in Tartary buckwheat, which is inconsistent with TT8 from other plant species, and FtTT8 might be a high-quality gene resource for Tartary buckwheat breeding.

Targeted metabolic reveals different part of maize in polyphenolic metabolites during germination and hypoglycemic activity analysis

In this study, qualitative and quantitative analyses of phenolic compounds in the maize germinating seed embryo, radicle, and germ were performed at 0, 48, and 96 h post-germination, followed by the evaluation of their hypoglycemic activity. The results revealed the accumulation of 80 phenolics in different parts of germinated maize, of which 47, 48, and 53 were present in the seed embryo, radicle, and germ. After germination 22, 26, and 34 polyphenols were found to differential accumulate in the seed embryo, radicle, and germ. At 96 h post-germination, the content of monomeric phenols in the germ was higher than that in the radicle and seed embryo. Moreover, the inhibitory activity of polyphenols in the germ towards α-glucosidase and α-amylase was higher than that in the radicle and seed embryo. These results indicate that germination can effectively improve the type and content of phenolic compounds in different parts of maize.

Integration of comparative transcriptomics and WGCNA characterizes the regulation of anthocyanin biosynthesis in mung bean (Vigna radiata L.)

Mung bean is a dual-use crop widely cultivated in Southeast Asia as a food and medicine resource. The development of new functional mung bean varieties demands identifying new genes regulating anthocyanidin synthesis and investigating their molecular mechanism. In this study, we used high-throughput sequencing technology to generate transcriptome sequence of leaves, petioles, and hypocotyls for investigating the anthocyanins accumulation in common mung bean variety as well as anthocyanidin rich mung bean variety, and to elucidate their molecular mechanisms. 29 kinds of anthocyanin compounds were identified. Most of the anthocyanin components contents were significantly higher in ZL23 compare with AL12. Transcriptome analysis suggested that a total of 93 structural genes encoding the anthocyanin biosynthetic pathway and 273 regulatory genes encoding the ternary complex of MYB-bHLH-WD40 were identified, of which 26 and 78 were differentially expressed in the two varieties. Weighted gene co-expression network analysis revealed that VrMYB3 and VrMYB90 might have enhanced mung bean anthocyanin content by inducing the expression of structural genes such as PAL, 4CL, F3'5'H, LDOX, and F3'H, which was consistent with qRT-PCR results. These findings are envisaged to provide a reference for studying the molecular mechanism of anthocyanin accumulation in mung beans.

Multi-Omics Elucidates Difference in Accumulation of Bioactive Constituents in Licorice (Glycyrrhiza uralensis) under Drought Stress

Licorice is a frequently applied herb with potential edible and medicinal value based on various flavonoids and triterpenes. However, studies on detailed flavonoid and triterpene metabolism and the molecular basis of their biosynthesis in licorice are very limited, especially under drought conditions. In the present study, we carried out transcriptome, proteome, and metabolome experiments. To ultimately combine three omics for analysis, we performed a bioinformatics comparison, integrating transcriptome data and proteome data through a Cloud platform, along with a simplified biosynthesis of primary flavonoids and triterpenoids in the KEGG pathway based on metabolomic results. The biosynthesis pathways of triterpenes and flavonoids are enriched at both gene and protein levels. Key flavonoid-related genes (PAL, 4CL, CHS, CHI, CYP93C, HIDH, HI4OMT, and CYP81E1_7) and representative proteins (HIDH, CYP81E1_7, CYP93C, and VR) were obtained, which all showed high levels after drought treatment. Notably, one R2R3-MYB transcription factor (Glyur000237s00014382.1), a critical regulator of flavonoid biosynthesis, achieved a significant upregulated expression as well. In the biosynthesis of glycyrrhizin, both gene and protein levels of bAS and CYP88D6 have been found with upregulated expression under drought conditions. Most of the differentially expressed genes (DEGs) and proteins (DEPs) showed similar expression patterns and positively related to metabolic profiles of flavonoid and saponin. We believe that suitable drought stress may contribute to the accumulation of bioactive constituents in licorice, and our research provides an insight into the genetic study and quality breeding in this plant.

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.

Identification of common buckwheat (Fagopyrum esculentum Moench) adulterated in Tartary buckwheat (Fagopyrum tataricum (L.) Gaertn) flour based on near-infrared spectroscopy and chemometrics

Near-infrared spectroscopy (NIRS) presents great potential in the identification of food adulteration due to its advantages of nondestructive, simple, and easy to operate. In this paper, a method based on NIRS and chemometrics was proposed to predict the content of common buckwheat (Fagopyrum esculentum Moench) flour in Tartary buckwheat (Fagopyrum tataricum (L.) Gaertn) flour. Partial least squares regression (PLSR) and support vector regression (SVR) models were used to analyze the spectrum data of adulterated samples and predict the adulteration level. Various preprocessing methods, parameter-optimization methods, and competitive adaptive reweighted sampling (CARS) wavelength-selection methods were used to optimize the model prediction accuracy. The results of PLSR and SVR modeling for predicting of Tartary buckwheat adulteration content were satisfactory, and the correlation coefficients of the optimum identification models were above 0.99. In conclusion, the combinations of NIRS and chemometrics indicated excellent predictive performance and applicability to analyze the adulteration of common buckwheat flour in Tartary buckwheat flour. This work provides a promising method to identify the adulteration of Tartary buckwheat flour and results obtained can give theoretical and data support for adulteration identification of agro-products.

The Tartary buckwheat bHLH gene ALCATRAZ enables Arabidopsis thaliana silique dehiscence

Tartary buckwheat is popular because of its rich nutrients. However, the difficulty in shelling restricts food production. The gene ALCATRAZ (AtALC) plays a key role in silique dehiscence in Arabidopsis thaliana. In this study, an atalc mutant was obtained by CRISPR/Cas9, and a FtALC gene homologous to AtALC was complemented into the atalc mutant to verify its function. Phenotypic observations showed that three atalc mutant lines did not dehiscence, while ComFtALC lines recovered the dehiscence phenotype. The contents of lignin, cellulose, hemicellulose, and pectin in the siliques of all the atalc mutant lines were significantly higher than those in the wild-type and ComFtALC lines. Moreover, FtALC was found to regulate the expression of cell wall pathway genes. Finally, the interaction of FtALC with FtSHP and FtIND was verified by yeast two-hybrid, bimolecular fluorescent complimentary (BIFC) and firefly luciferase completion imaging assays (LCIs). Our findings enrich the silique regulatory network and lay the foundation for the cultivation of easily shelled tartary buckwheat varieties.

Metabolomics and transcriptomics provide insights into the molecular mechanisms of anthocyanin accumulation in the seed coat of differently colored mung bean (Vigna radiata L.)

Black mung bean is rich in anthocyanin, however, the accumulation and the molecular mechanism of anthocyanin synthesis in black mung bean are unclear. In this study, anthocyanin metabolomics and transcriptomics on the seed coats of two different colors of mung bean were performed to clarify the composition of anthocyanins, and identify transcription factors involved in regulating anthocyanin biosynthesis. In the mature stage, 23 kinds of anthocyanin compounds were identified. All anthocyanin components contents were significantly higher in seed coat of black mung bean compare with green mung bean. Transcriptome analysis suggested that most of the structural genes for anthocyanin biosynthesis and some potential regulatory genes were significantly differentially expressed. WGCNA suggested VrMYB90 was an important regulatory gene in anthocyanin biosynthesis. Arabidopsis thaliana overexpressing VrMYB90 showed significant accumulation of anthocyanins. PAL, 4CL, DFR, F3'5'H, LDOX, F3'H and UFGT were up-regulated in 35S:VrMYB90 Arabidopsis thaliana. These findings provide valuable information for understanding the synthesis mechanism of anthocyanins in black mung bean seed coats.

Proteomic and metabolomic analyses showing the differentially accumulation of NnUFGT2 is involved in the petal red-white bicolor pigmentation in lotus (Nelumbo nucifera)

Bicolor flower lotus is rare with high ornamental value. During the long history of breeding and artificial selection, a very famous lotus cultivar 'Da Sajin' with red and white picotee bicolor petals were obtained. In order to reveal the mechanism underlying the formation of its picotee bicolor pattern in the petal, an integrative metabolomics and proteomics analyses were conducted between red and white parts of its petals. The results showed that the defect of anthocyanidin 3-O-glucosyltransferases (UFGTs) accumulation resulted in the failure of the glycosylation of anthocyanidin, the last step of anthocyanin biosynthesis in white part of the petals. And proteomic data and biochemical analysis showed that the defect of UFGTs accumulation is not related to their transcription, but because of their degradation. Function of one differentially accumulated NnUFGT were proven being involved in anthocyanin biosynthesis through both in-vitro enzyme assay and in-vivo transgenic analyses. This regulation on the protein accumulation of structural genes in anthocyanin biosynthesis was not explored in any other plants, and hence supposed to be a novel mechanism for the formation of picotee bicolor pattern flower. The results not only provide some new insights into the understanding of lotus flower coloration, but also might assist the breeding of flower lotus.

Tartary Buckwheat Grain as a Source of Bioactive Compounds in Husked Groats

Tartary buckwheat (Fagopyrum tataricum Gaertn.) originates in mountain regions of Western China, and is cultivated in China, Bhutan, Northern India, Nepal, and Central Europe. The content of flavonoids in Tartary buckwheat grain and groats is much higher than in common buckwheat (Fagopyrum esculentum Moench), and depends on ecological conditions, such as UV-B radiation. Buckwheat intake has preventative effects in chronic diseases, such as cardiovascular diseases, diabetes, and obesity, due to its content of bioactive substances. The main bioactive compounds in Tartary buckwheat groats are flavonoids (rutin and quercetin). There are differences in the bioactivities of buckwheat groats obtained using different husking technologies, based on husking raw or pretreated grain. Husking hydrothermally pretreated grain is among the traditional ways of consuming buckwheat in Europe and some parts of China and Japan. During hydrothermal and other processing of Tartary buckwheat grain, a part of rutin is transformed to quercetin, the degradation product of rutin. By adjusting the humidity of materials and the processing temperature, it is possible to regulate the degree of conversion of rutin to quercetin. Rutin is degraded to quercetin in Tartary buckwheat grain due to the enzyme rutinosidase. The high-temperature treatment of wet Tartary buckwheat grain is able to prevent the transformation of rutin to quercetin.

Widely targeted metabolomics analysis of different parts of kudzu

Pueraria lobata is a traditional medicinal and edible plant. Its root is often used as a nutritional supplement, but its stems and leaves are often discarded. In this study, the types and contents of compounds in roots, stems and leaves of kudzu were studied by ultra-high-performance liquid chromatography electrospray ionization tandem mass spectrometry. A total of 446 metabolites were identified, which were mainly divided into eight categories; 40 unique compounds were detected in roots. Multivariate statistical analysis showed that there were significant differences in the contents of metabolites in the roots, stems and leaves of kudzu, and the most significant differences were found in the roots and stems. Most flavonoids and terpenoids accumulated in roots, while alkaloids, flavonoid glycosides, lignans and coumarins accumulated in stems and leaves. In addition, isoflavones accumulated least in stems and most in roots. These results improved our understanding of the accumulation of metabolites in kudzu and provided a reference for the full study of the medicinal value of different parts of kudzu.

Metabolomics and transcriptomics provide insights into the flavonoid biosynthesis pathway in the roots of developing Aster tataricus

Aster tataricus (L.) is an important medicinal plant in China. Its roots are rich in flavonoids, the main medicinal components. However, the molecular basis of flavonoid biosynthesis in the roots of A. tataricus remains unclear. In this study, the content of total flavonoid of A. tataricus roots at different developmental stages was measured first, and the results showed that the content of total flavonoid gradually decreased from September to November, which may be caused by the stagnation of A. tataricus growth due to the decrease in temperature after September. Then, an integrated analysis of transcriptome and metabolome was conducted on five developing stages of A. tataricus roots to identify flavonoid compositions and potential genes involved in flavonoid biosynthesis. A total of 80 flavonoid metabolites, of which 75% were flavonols and flavonoids, were identified in metabolomic analyses, among which isorhamnetin, kaempferol, quercetin, and myricetin were the main skeletons of these flavonoids. Cluster analysis divided these 80 flavonoids into 3 clusters. The compounds in cluster I mainly accumulated in S1, S3, and S5. In cluster II, the relative content of the flavonoid metabolites showed an upward trend from S2 to S4. In cluster III, the flavonoids decreased from S1 to S5. A total of 129 structural genes, including 43 PAL, 23 4CL, 9 C4H, 4 CHS, 18 CHI, 3 F3H, 5 F3'H, 1 F3'5'H, 21 FLS, and 2 FSII, and 65 transcription factors, including 22 AP2/ERF, 7 bHLH, 5 bZIP, 8 MYB, 11 NAC, and 12 WRKY, showed significant correlation with total flavonoid content. Eighteen genes (7 4CL, 5 C4H, 2 CHI, 1 F3H, and 3 FLS) and 30 genes (5 PAL, 9 4CL, 1 C4H, 2 CHI, 1 F3H, 1 DFR, 7 3AT, 1 BZ1, and 3 UGT79B1) were identified as key structural genes for kaempferol and anthocyanins biosynthesis, respectively. Our study provides valuable information for understanding the mechanism of flavonoid biosynthesis in A. tataricus root.

Discovery of the key active compounds in Citri Reticulatae Pericarpium (Citrus reticulata "Chachi") and their therapeutic potential for the treatment of COVID-19 based on comparative metabolomics and network pharmacology

Edible herbal medicines contain macro- and micronutrients and active metabolites that can take part in biochemical processes to help achieve or maintain a state of well-being. Citri Reticulatae Pericarpium (CRP) is an edible and medicinal herb used as a component of the traditional Chinese medicine (TCM) approach to treating COVID-19 in China. However, the material basis and related mechanistic research regarding this herb for the treatment of COVID-19 are still unclear. First, a wide-targeted UPLC-ESI-MS/MS-based comparative metabolomics analysis was conducted to screen for the active metabolites of CRP. Second, network pharmacology was used to uncover the initial linkages among these metabolites, their possible targets, and COVID-19. Each metabolite was then further studied via molecular docking with the identified potential SARS-CoV-2 targets 3CL hydrolase, host cell target angiotensin-converting enzyme II, spike protein, and RNA-dependent RNA polymerase. Finally, the most potential small molecule compound was verified by in vitro and in vivo experiments, and the mechanism of its treatment of COVID-19 was further explored. In total, 399 metabolites were identified and nine upregulated differential metabolites were screened out as potential key active metabolites, among which isorhamnetin have anti-inflammatory activity in vitro validation assays. In addition, the molecular docking results also showed that isorhamnetin had a good binding ability with the key targets of COVID-19. Furthermore, in vivo results showed that isorhamnetin could significantly reduced the lung pathological injury and inflammatory injury by regulating ATK1, EGFR, MAPK8, and MAPK14 to involve in TNF signaling pathway, PI3K-Akt signalling pathway, and T cell receptor signaling pathway. Our results indicated that isorhamnetin, as screened from CRP, may have great potential for use in the treatment of patients with COVID-19. This study has also demonstrated that comparative metabolomics combined with network pharmacology strategy could be used as an effective approach for discovering potential compounds in herbal medicines that are effective against COVID-19.

Comparative proteomic analyses of Tartary buckwheat (Fagopyrum tataricum) seeds at three stages of development

Tartary buckwheat is among the valuable crops, utilized as both food and Chinese herbal medicine. To uncover the accumulation dynamics of the main nutrients and their regulatory mechanism of Tartary buckwheat seeds, microscopic observations and nutrient analysis were conducted which suggested that starch, proteins as well as flavonoid gradually accumulated among seed development. Comparative proteomic analysis of rice Tartary buckwheat at three different developmental stages was performed. A total of 78 protein spots showed differential expression with 74 of them being successfully identified by MALDI-TOF/TOF MS. Among them, granule bound starch synthase (GBSS1) might be the critical enzyme that determines starch biosynthesis, while 11 S seed storage protein and vicilin seemed to be the main globulin and affect seed storage protein accumulation in Tartary buckwheat seeds. Two enzymes, flavanone 3-hydroxylase (F3H) and anthocyanidin reductase (ANR), involved in the flavonoid biosynthesis pathway were identified. Further analysis on the expression profiles of flavonoid biosynthetic genes revealed that F3H might be the key enzyme that promote flavonoid accumulation. This study provides insights into the mechanism of nutrition accumulation at the protein level in Tartary buckwheat seeds and may facilitate in the breeding and enhancement of Tartary buckwheat germplasm.

Integrated metabolite profiling and transcriptome analysis reveal candidate genes involved in the formation of yellow Nelumbo nucifera

As a worldwide major ornamental flower and a edible plant, lotus (Nelumbo nucifera) is also used as medicine and tea beverage. Here, transcriptome and metabolites of yellow (MLQS) and white (YGB) lotus cultivars during five key flower coloration stages were profiled. 2014 differentially expressed genes were detected with 11 carotenoids in lotus were identified for the first time. Then, regulatory networks between and within functional modules was reconstructed, and the correlation between module-metabolites and gene-metabolites was conducted within 3 core modules. 18 candidate genes related to the formation of yellow flower were screened out and a gene regulatory model for the flower color difference between MLQS and YGB were speculated as follows: The substrate competition between F3'H and F3'5'H and substrate specificity of FLS, together with differential expression of CCD4a and CCD4b were contribute to the differences in flavonoids and carotenoids. Besides, UGT73C2, UGT91C1-2 and SGTase, and regulation of UGTs by transcription factors PLATZ, MADS, NAC031, and MYB308 may also play a role in the upstream regulation. The following verification results indicated that functional differences existed in the coding sequences of NnCCD4b and promoters of NnCCD4a of MLQS and YGB. In all, this study preliminarily reveals the mechanism of yellow flower coloration in lotus and provides new ideas for the study of complex ornamental characters of other plants.

Metabolomic and transcriptomic exploration of the uric acid-reducing flavonoids biosynthetic pathways in the fruit of Actinidia arguta Sieb. Zucc

Flavonoids from Actinidia arguta Sieb. Zucc. can reduce uric acid in mice. However, the molecular basis of its biosynthesis is still unclear. In this paper, we used a combination of extensively targeted metabolomics and transcriptomics analysis to determine the types and differences of flavonoids in the fruit ripening period (August to September) of two main cultivated varieties in northern China. The ethanol extract was prepared, and the potential flavonoids of Chrysin (Flavone1), Rutin (Flavone2), and Daidzein (Flavone3) in Actinidia arguta Sieb. Zucc. were separated and purified by HPD600 macroporous adsorption resin and preparative liquid chromatography. The structure was identified by MS-HPLC, and the serum uric acid index of male Kunming mice was determined by an animal model test.125 flavonoids and 50 differentially regulated genes were identified. The contents of UA (uric acid), BUN (urea nitrogen), Cr (creatinine), and GAPDH in mouse serum and mouse liver glycogen decreased or increased in varying degrees. This paper reveals the biosynthetic pathway of uric acid-reducing flavonoids in the fruit of Actinidia arguta Sieb. Zucc., a major cultivar in northern China, provides valuable information for the development of food and drug homologous functional foods.

Niche differentiation modulates metabolites abundance and composition in silicon fertilizer amended soil during sugarcane growth

BACKGROUND

As one of the vital crops globally, sugarcane (Saccharum officinarum L.) has been one of model crops for conducting metabolome research. Although many studies have focused on understanding bioactive components in specific sugarcane tissues, crucial questions have been left unanswered about the response of metabolites to niche differentiation such as different sugarcane tissues (leaf, stem and root), and soil regions (rhizosphere and bulk) under silicon (Si) amended soils. Here, nontargeted metabolite profiling method was leveraged to assess the similarities and differences in the abundance and community composition of metabolites in the different sugarcane and soil compartments. Identify the compartment-specific expression patterns of metabolites, and their association with cane agronomic traits and edaphic factors. We also investigated the response of sugarcane agronomic traits and edaphic factors to Si amended soil.

RESULTS

We found that Si fertilizer exhibited the advantages of overwhelmingly promoting the height and theoretical production of cane, and profoundly increased soil Si content by 24.8 and 27.0%, while soil available potassium (AK) was enhanced by 3.07 and 2.67 folds in the bulk and rhizosphere soils, respectively. It was also noticed that available phosphorus (AP) in the rhizosphere soil tremendously increased by 105.5%. We detected 339 metabolites in 30 samples using LC-MS/MS analyses, 161 of which were classified and annotated, including organooxygen compounds (19.9%), carboxylic acids and derivatives (15.5%), fatty acyls (15.5%), flavonoids (4.4%), phenols (4.4%), and benzene and substituted derivatives (3.7%). In addition, the total percentages covered by these core metabolites in each compartment ranged from 94.0% (bulk soil) to 93.4% (rhizosphere soil), followed by 87.4% (leaf), 81.0% (root) and 80.5% (stem), suggesting that these bioactive compounds may have migrated from the belowground tissues and gradually filtered in various aboveground niches of the plant. We also observed that the variations and enrichment of metabolites abundance and community were compartment-specific. Furthermore, some key bioactive compounds were markedly associated with plant growth parameters and soil edaphic.

CONCLUSION

Taken together, we hypothesized that Si utilization can exhibit the advantage of enhancing edaphic factors and cane agronomic traits, and variations in metabolites community are tissue-specific.

Selenium Application Enhances the Accumulation of Flavones and Anthocyanins in Bread Wheat (Triticum aestivum L.) Grains

Selenium (Se) biofortification in wheat reduces the risk of Se deficiency in humans. Se biofortification increases the concentration of Se and anthocyanins in wheat grains. However, it is unknown whether Se biofortification can enhance flavonoids other than anthocyanins and the mechanism underlying flavonoid accumulation in wheat grains. Here, foliar application of selenite solution in wheat was conducted 10 days after flowering. Metabolite profiling and transcriptome sequencing were performed in Se-treated grains. A significant increase in the total contents of Se, anthocyanins, and flavonoids was observed in Se-treated mature grains. Twenty-seven significantly increased flavonoids were identified in Se-treated immature grains. The significant accumulation of flavones (tricin, tricin derivatives, and chrysoeriol derivatives) was detected, and six anthocyanins, dihydroquercetin (the precursor for anthocyanin biosynthesis) and catechins were also increased. Integrated analysis of metabolites and transcriptome revealed that Se application enhanced the biosynthesis of flavones, dihydroquercetin, anthocyanins, and catechins by increasing the expression levels of seven key structural genes in flavonoid biosynthesis (two TaF3Hs, two TaDFRs, one TaF3'5'H, one TaOMT, and one TaANR). Our findings shed new light on the molecular mechanism underlying the enhancement in flavonoid accumulation by Se supplementation and pave the way for further enhancing the nutritional value of wheat grains.

Different adaptive patterns of wheat with different drought tolerance under drought stresses and rehydration revealed by integrated metabolomic and transcriptomic analysis

Wheat as a staple food crop is enduring ever-frequent intermittent and changing drought with the climate change. It is of great significance to highlight the adaptive approaches under such variable conditions at multiple levels to provide a comprehensive understanding of drought tolerance and facilitate the genetic breeding of wheat. Therefore, three wheat lines with different drought tolerance (drought-tolerant mutant Mu > common wheat CK > drought susceptible mutant mu) were analyzed under moderate and severe drought stresses as well as rehydration. Samples were subjected to transcriptomic and metabolomic profiling in combination with physiological and biochemical determination. The moderate drought stress rendered 198 and 115 differentially expressed metabolites (DEMs) in CK and Mu, respectively. The severe drought stress rendered 166, 151 and 137 DEMs in CK, Mu and mu, respectively. The rehydration rendered 150 and 127 DEMs in CK and Mu. 12,557 and 10,402 differentially expressed genes (DEGs) were identified for CK and Mu under moderate drought stress, respectively. 9,893, 7,924, and 9,387 DEGs were identified for CK, Mu, and mu under severe drought stress, respectively. 13,874 and 14,839 were identified in CK and Mu under rehydration, respectively. Metabolomics results showed that amino acid was the most differentially expressed metabolites, followed by phenolic acids. Flavonoids played an important role in drought tolerance. Most enriched pathways under drought included biosynthesis of secondary metabolites, metabolic pathways and photosynthesis. Metabolites and genes involved in osmotic regulation, antioxidase activities, and ABA signaling were more enriched in Mu than in CK and mu. Various drought-responsive genes and metabolites in Mu showed different trends with those in CK and mu. Increased amino acids biosynthetic capability and ROS scavenging ability resulted from higher antioxidase activities and increased flavonoids may be the mechanisms underlying the drought tolerance characteristic of Mu. Recovery from reversible ROS damage and rapid amino acid biosynthesis may contribute to the rapid recovery of Mu. The present study provides new insights for mechanisms of wheat under complex drought conditions.

Comparative physiological, transcriptomic, and WGCNA analyses reveal the key genes and regulatory pathways associated with drought tolerance in Tartary buckwheat

Drought stress is one of the major abiotic stress factors that affect plant growth and crop productivity. Tartary buckwheat is a nutritionally balanced and flavonoid-rich pseudocereal crop and also has strong adaptability to different adverse environments including drought. However, little is known about its drought tolerance mechanism. In this study, we performed comparative physiological and transcriptomic analyses of two contrasting drought-resistant Tartary buckwheat genotypes under nature drought treatment in the reproductive stage. Under drought stress, the drought-tolerant genotype XZSN had significantly higher contents of relative water, proline, and soluble sugar, as well as lower relative electrolyte leakage in the leaves than the drought-susceptible LK3. A total of 5,058 (2,165 upregulated and 2,893 downregulated) and 5,182 (2,358 upregulated and 2,824 downregulated) potential drought-responsive genes were identified in XZSN and LK3 by transcriptome sequencing analysis, respectively. Among the potential drought-responsive genes of XZSN, 1,206 and 1,274 genes were identified to be potential positive and negative contributors for XZSN having higher drought resistance ability than LK3. Furthermore, 851 out of 1,206 positive drought-resistant genes were further identified to be the core drought-resistant genes of XZSN based on WGCNA analysis, and most of them were induced earlier and quicker by drought stress than those in LK3. Functional annotation of the 851 core drought-resistant genes found that a large number of stress-responsive genes were involved in TFs, abscisic acid (ABA) biosynthesis, signal transduction and response, non-ABA signal molecule biosynthesis, water holding, oxygen species scavenging, osmotic adjustment, cell damage prevention, and so on. Transcriptional regulatory network analyses identified the potential regulators of these drought-resistant functional genes and found that the HD-ZIP and MYB TFs might be the key downstream TFs of drought resistance in Tartary buckwheat. Taken together, these results indicated that the XZSN genotype was more drought-tolerant than the LK3 genotype as evidenced by triggering the rapid and dramatic transcriptional reprogramming of drought-resistant genes to reduce water loss, prevent cell damage, and so on. This research expands our current understanding of the drought tolerance mechanisms of Tartary buckwheat and provides important information for its further drought resistance research and variety breeding.

Integrative Analysis of the Metabolome and Transcriptome Provides Insights into the Mechanisms of Flavonoid Biosynthesis in Quinoa Seeds at Different Developmental Stages

Quinoa (Chenopodium quinoa Willd.) is a crop with high nutritional and health benefits. Quinoa seeds are rich in flavonoid compounds; however, the mechanisms behind quinoa flavonoid biosynthesis remain unclear. We independently selected the high-generation quinoa strain ‘Dianli-3260′, and used its seeds at the filling, milk ripening, wax ripening, and mature stages for extensive targeted metabolome analysis combined with joint transcriptome analysis. The results showed that the molecular mechanism of flavonoid biosynthesis in quinoa seeds was mainly concentrated in two pathways: “flavonoid biosynthesis pathway” and “flavone and flavonol biosynthesis pathway”. Totally, 154 flavonoid-related metabolites, mainly flavones and flavonols, were detected in the four development stages. Moreover, 39,738 genes were annotated with KEGG functions, and most structural genes of flavonoid biosynthesis were differentially expressed during grain development. We analyzed the differential flavonoid metabolites and transcriptome changes between the four development stages of quinoa seeds and found that 11 differential flavonoid metabolites and 22 differential genes were the key factors for the difference in flavonoid biosynthesis. This study provides important information on the mechanisms underlying quinoa flavonoid biosynthesis, the screening of potential quinoa flavonoid biosynthesis regulation target genes, and the development of quinoa products.

Transcriptomic and metabolomic profiling of flavonoid biosynthesis provides novel insights into petals coloration in Asian cotton (Gossypium arboreum L.)

BACKGROUND

Asian cotton (Gossypium arboreum L.), as a precious germplasm resource of cotton with insect resistance and stress tolerance, possesses a broad spectrum of phenotypic variation related to pigmentation. Flower color affects insect pollination and the ornamental value of plants. Studying flower color of Asian cotton varieties improves the rate of hybridization and thus enriches the diversity of germplasm resources. Meanwhile, it also impacts the development of the horticultural industry. Unfortunately, there is a clear lack of studies concerning intricate mechanisms of cotton flower-color differentiation. Hereby, we report an integrative approach utilizing transcriptome and metabolome concerning flower color variation in three Gossypium arboreum cultivars.

RESULTS

A total of 215 differentially accumulated metabolites (DAMs) were identified, including 83 differentially accumulated flavonoids (DAFs). Colorless kaempferol was more abundant in white flowers, while gossypetin-fer showed specificity in white flowers. Quercetin and gossypetin were the main contributors to yellow petal formation. Pelargonidin 3-O-beta-D-glucoside and cyanidin-3-O-(6''-Malonylglucoside) showed high accumulation levels in purple petals. Quercetin and gossypetin pigments also promoted purple flower coloration. Moreover, 8178 differentially expressed genes (DEGs) were identified by RNA sequencing. The correlation results between total anthocyanins and DEGs were explored, indicating that 10 key structural genes and 29 transcription factors promoted anthocyanin biosynthesis and could be candidates for anthocyanin accumulation. Ultimately, we constructed co-expression networks of key DAFs and DEGs and demonstrated the interactions between specific metabolites and transcripts in different color flowers.

CONCLUSION

This study provides new insights into elucidating the regulatory mechanisms of cotton flower color and lays a potential foundation for generate cotton varieties with highly attractive flowers for pollinators.

Visualizing the spatial distribution of functional metabolites in Forsythia suspensa at different harvest stages by MALDI mass spectrometry imaging

As a traditional Chinese medicine, Forsythia suspensa (F. suspensa) has attracted much attention due to its significant pharmacological activity. Revealing the spatial distribution of metabolites during F. suspensa development is important for understanding its biosynthesis rules and improving the quality of medicinal materials. However, there is currently a lack of information on the spatial distribution of F. suspensa metabolites. In this work, the spatial distribution and growth metabolism patterns of important metabolites of F. suspensa were studied for the first time using matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI). Using 2,5-dimethylnaphthalene (DAN) as the matrix and detecting in negative ion mode, the spatial distribution and growth patterns of 11 metabolites obtained from longitudinal sections of F. suspensa included pinoresinol, phillygenin, forsythoside A, forsythoside E, rutin, caffeic acid, malic acid, citric acid, stearic acid, oleic acid, and linoleic acid. These results showed the mesocarp and endosperm tissues of F. suspensa were important for storing important functional metabolites. Changes in mesocarp and endosperm growth and development tissues caused large changes in the content of important functional metabolites in F. suspensa. These results provide a basis for understanding the spatial distribution of metabolites in F. suspensa tissues and the significant changes that occur during growth and development, exploring the mechanism of important synthesis of metabolites, regulating the harvest of F. suspensa, and improving the quality of medicinal herbs.

Dynamic transcriptome analysis suggests the key genes regulating seed development and filling in Tartary buckwheat (Fagopyrum tataricum Garetn.)

Tartary buckwheat is highly attractive for the richness of nutrients and quality, yet post-embryonic seed abortion greatly halts the yield. Seed development is crucial for determining grain yield, whereas the molecular basis and regulatory network of Tartary buckwheat seed development and filling is not well understood at present. Here, we assessed the transcriptional dynamics of filling stage Tartary buckwheat seeds at three developmental stages by RNA sequencing. Among the 4249 differentially expressed genes (DEGs), genes related to seed development were identified. Specifically, 88 phytohormone biosynthesis signaling genes, 309 TFs, and 16 expansin genes participating in cell enlargement, 37 structural genes involved in starch biosynthesis represented significant variation and were candidate key seed development genes. Cis-element enrichment analysis indicated that the promoters of differentially expressed expansin genes and starch biosynthesis genes are rich of hormone-responsive (ABA-, AUX-, ET-, and JA-), and seed growth-related (MYB, MYC and WRKY) binding sites. The expansin DEGs showed strong correlations with DEGs in phytohormone pathways and transcription factors (TFs). In total, phytohormone ABA, AUX, ET, BR and CTK, and related TFs could substantially regulate seed development in Tartary buckwheat through targeting downstream expansin genes and structural starch biosynthetic genes. This transcriptome data could provide a theoretical basis for improving yield of Tartary buckwheat.

Transcriptome and Metabolome Analyses Reveal Differences in Terpenoid and Flavonoid Biosynthesis in Cryptomeria fortunei Needles Across Different Seasons

Cryptomeria fortunei (Chinese cedar) has outstanding medicinal value due to its abundant flavonoid and terpenoid contents. The metabolite contents of C. fortunei needles differ across different seasons. However, the biosynthetic mechanism of these differentially synthesized metabolites (DSMs) is poorly understood. To improve our understanding of this process, we performed integrated non-targeted metabolomic liquid chromatography and gas chromatography mass spectrometry (LC-MS and GC-MS), and transcriptomic analyses of summer and winter needles. In winter, the C. fortunei needle ultrastructure was damaged, and the chlorophyll content and F _v/F _m were significantly (p < 0.05) reduced. Based on GC-MS and LC-MS, we obtained 106 and 413 DSMs, respectively; based on transcriptome analysis, we obtained a total of 41.17 Gb of clean data and assembled 33,063 unigenes, including 14,057 differentially expressed unigenes (DEGs). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that these DSMs/DEGs were significantly (p < 0.05) enriched in many biosynthesis pathways, such as terpenoids, photosynthates, and flavonoids. Integrated transcriptomic and metabonomic analyses showed that seasonal changes have the greatest impact on photosynthesis pathways, followed by terpenoid and flavonoid biosynthesis pathways. In summer Chinese cedar (SCC) needles, DXS, DXR, and ispH in the 2-methyl-pentaerythritol 4-phosphate (MEP) pathway and GGPS were highly expressed and promoted the accumulation of terpenoids, especially diterpenoids. In winter Chinese cedar (WCC) needles, 9 genes (HCT, CHS, CHI, F3H, F3'H, F3'5'H, FLS, DFR, and LAR) involved in flavonoid biosynthesis were highly expressed and promoted flavonoid accumulation. This study broadens our understanding of the metabolic and transcriptomic changes in C. fortunei needles caused by seasonal changes and provides a reference regarding the adaptive mechanisms of C. fortunei and the extraction of its metabolites.

Transcriptome analysis reveals regulation mechanism of methyl jasmonate-induced terpenes biosynthesis in Curcuma wenyujin

Curcuma wenyujin is the source plant of three traditional Chinese medicines, which have been widely used in clinical treatment over 1000 years. The content of terpenes, the major medicinal active ingredients, is relatively low in this plant. Studies have shown that MeJA can promote terpenes biosynthesis in plants. However, the mechanism underlying the effect of MeJA in C. wenyujin remains unclear. In this work, the transcriptome of C. wenyujin leaves with MeJA treatment was analyzed to elucidate the regulation mechanism of MeJA-mediated terpene biosynthesis. Based on the RNA-seq data, 7,246 unigenes were differentially expressed with MeJA treatment. Expression pattern clustering of DEGs revealed that unigenes, related to JA biosynthesis and signal transduction, responded to exogenous MeJA stimulation on the early stage and maintained throughout the process. Subsequently, unigenes related to terpene biosynthesis pathway showed a significant up-regulation with 6 h treatment. The analysis results suggested that MeJA induced the expression of JA biosynthesis genes (such as LOXs, AOSs, AOCs, OPRs, and MFPs) and JA signal transduction core genes (JAZs and MYCs) to activate JA signaling pathway. Meanwhile, downstream JA-responsive genes presented up-regulated expression levels such as AACT, HMGSs, HMGRs, DXSs, DXRs, MCTs, HDSs, and HDRs, thus promoting terpenes biosynthesis. The transcriptional expressions of these genes were validated by qRT-PCR. In addition, six CwTPS genes in response to MeJA were identified. With MeJA treatment, the expression levels of CwTPSs were increased as well as those of the transcription factors MYB, NAC, bZIP, WRKY, AP2/ERF, and HLH. These TFs might potentially regulate terpenes biosynthesis. These results provide insights for regulation mechanism of terpenes biosynthesis.

FtUGT79A15 is responsible for rutinosylation in flavonoid diglycoside biosynthesis in Fagopyrum tataricum

Tartary buckwheat shows health benefits with its high antioxidant activity and abundant flavonoid content. However, glycosylated flavonoid accumulation patterns and their molecular basis remain unidentified in Tartary buckwheat. Here, our metabolomics analysis revealed that F3'H branching was the major flavonoid metabolic flux in Tartary buckwheat. Interestingly, metabolome results also showed that the most abundant flavonoids were mainly in the glycosylated form, including flavonoid glycosides and flavonoid diglycosides in Tartary buckwheat. However, the flavonoid glycosides glycosyltransferase (GGT) gene catalyzing the second glycosylation step of flavonoid diglycoside has not been discovered yet in Tartary buckwheat. Thus, we explored GGT genes in the transcriptome-metabolome correlation network and confirmed that FtUGT79A15 showed the rhamnosyltransferase activity to catalyze quercetin 3-O-glucoside to rutin invitro and inplanta. Overall, FtUGT79A15 was identified to involve in the flavonoid diglycoside biosynthesis pathway in Tartary buckwheat.

Metabolome and Transcriptome Analysis Reveals the Transcriptional Regulatory Mechanism of Triterpenoid Saponin Biosynthesis in Soapberry (Sapindus mukorossi Gaertn.)

Soapberry (Sapindus mukorossi Gaertn.) pericarps are rich in valuable bioactive triterpenoid saponins. However, the saponin content dynamics and the molecular regulatory network of saponin biosynthesis in soapberry pericarps remain largely unclear. Here, we performed combined metabolite profiling and transcriptome analysis to identify saponin accumulation kinetic patterns, investigate gene networks, and characterize key candidate genes and transcription factors (TFs) involved in saponin biosynthesis in soapberry pericarps. A total of 54 saponins were tentatively identified, including 25 that were differentially accumulated. Furthermore, 49 genes putatively involved in sapogenin backbone biosynthesis and some candidate genes assumed to be responsible for the backbone modification, including 41 cytochrome P450s and 45 glycosyltransferases, were identified. Saponin-specific clusters/modules were identified by Mfuzz clustering and weighted gene coexpression network analysis, and one TF-gene regulatory network underlying saponin biosynthesis was proposed. The results of yeast one-hybrid assay and electrophoretic mobility shift assay suggested that SmbHLH2, SmTCP4, and SmWRKY27 may play important roles in the triterpenoid saponin biosynthesis by directly regulating the transcription of SmCYP71D-3 in the soapberry pericarp. Overall, these findings provide valuable information for understanding the molecular regulatory mechanism of saponin biosynthesis, enriching the gene resources, and guiding further research on triterpenoid saponin accumulation in soapberry pericarps.

Comparative metabolomics study of Tartary (Fagopyrum tataricum (L.) Gaertn) and common (Fagopyrum esculentum Moench) buckwheat seeds

Tartary buckwheat has higher health-promoting value than common buckwheat. However, the related metabolites information except flavonoids is largely deficient. Here, we compared the seed metabolomes of the two species using a UHPLC-QqQ-MS-based metabolomics approach. In total, 722 metabolites were obtained, of which 84 and 78 were identified as the key active ingredients of Traditional Chinese Medicines and the active pharmaceutical ingredients for six major diseases-resistance, respectively. Comparative analysis showed there were obviously difference in metabolic profiles between the two buckwheat species, and further found 61 flavonoids and 94 non-flavonoids metabolites displayed significantly higher contents (≥2 fold) in Tartary buckwheat than in common buckwheat. Our results suggest that Tartary and common buckwheat seeds are rich in metabolites beneficial to human health, and non-flavonoids metabolites also contributed to Tartary buckwheat's higher health-promoting value than common buckwheat. This study provides valuable information for the development of new functional foods of Tartary buckwheat.

Chemical Profile, Antimicrobial and Antioxidant Activity Assessment of the Crude Extract and Its Main Flavonoids from Tartary Buckwheat Sprouts

The purpose of this study was to investigate the major flavonoids content and bioactivities of Tartary buckwheat sprouts. The crude methanol extract (ME) of Tartary buckwheat sprouts was abundant in flavonoids, and six major flavonoids, including isoorientin, vitexin, isovitexin, rutin, quercetin, and kaemferol were successfully determined from the sprouts by the high-performance liquid chromatography (HPLC) method. Generally, the flavonoid content of buckwheat sprouts was in the order of rutin > quercetin > isovitexin > vitexin> isoorientin > kaemferol. The highest rutin content of the ME and sprout cultures was 89.81 mg/g and 31.50 mg/g, respectively. Antibacterial activity results indicated the ME displayed notable inhibitory activity against the five tested bacteria, and its minimum inhibitory concentration (MIC) values ranged from 0.8 mg/mL to 3.2 mg/mL. Among the six flavonoids, quercetin was the most active compound, which exhibited strong activity against all tested bacteria except for E. coli and S. epidermidis, with its MIC values ranging from 0.2 mg/mL to 0.4 mg/mL. For the antifungal activity assay, the ME of Tartary buckwheat sprouts and four flavonoids could significantly inhibit the spore germination of two pathogenic fungi, and their inhibitory efficiency was concentration dependent. Quercetin was the most active one, which significantly inhibited the spore germination of F. oxysporum f. sp. vasinfectum and F. oxysporum f. sp. cucumerinum, and its median effective inhibitory concentration (IC50) value was 42.36 and 32.85 µg/mL, respectively. The antioxidant activity results showed that quercetin, kaemferol, and rutin displayed excellent antioxidant activity in the DPPH radical scavenging test, and their IC50 value was calculated as 5.60, 16.23, and 27.95 µg/mL, respectively. This is the first report on the antimicrobial activity of the crude extract of Tartary buckwheat sprouts. These results indicated that the methanol extract of Tartary buckwheat sprouts could be used as a potential antimicrobial or antioxidant agent in the future.

Integrated Transcriptomics and Widely Targeted Metabolomics Analyses Provide Insights Into Flavonoid Biosynthesis in the Rhizomes of Golden Buckwheat (Fagopyrum cymosum)

Golden buckwheat (Fagopyrum cymosum) is used in Traditional Chinese Medicine. It has received attention because of the high value of its various medicinal and nutritional metabolites, especially flavonoids (catechin and epicatechin). However, the metabolites and their encoding genes in golden buckwheat have not yet been identified in the global landscape. This study performed transcriptomics and widely targeted metabolomics analyses for the first time on rhizomes of golden buckwheat. As a result, 10,191 differentially expressed genes (DEGs) and 297 differentially regulated metabolites (DRMs) were identified, among which the flavonoid biosynthesis pathway was enriched in both transcriptome and metabolome. The integration analyses of the transcriptome and the metabolome revealed a network related to catechin, in which four metabolites and 14 genes interacted with each other. Subsequently, an SG5 R2R3-MYB transcription factor, named FcMYB1, was identified as a transcriptional activator in catechin biosynthesis, as it was positively correlated to eight flavonoid biosynthesis genes in their expression patterns and was directly bound to the promoters of FcLAR2 and FcF3'H1 by yeast one hybrid analysis. Finally, a flavonoid biosynthesis pathway was proposed in the rhizomes of golden buckwheat, including 13 metabolites, 11 genes encoding 9 enzymes, and 1 MYB transcription factor. The expression of 12 DEGs were validated by qRT-PCR, resulting in a good agreement with the Pearson R ranging from 0.83 to 1. The study provided a comprehensive flavonoid biosynthesis and regulatory network of golden buckwheat.

Transcriptome sequencing and metabolomics analyses provide insights into the flavonoid biosynthesis in Torreya grandis kernels

The kernel of Torreya grandis (T. grandis) is a rare nut with a variety of bioactive compounds. Flavonoids are a very important class of bioactive compounds with high antioxidant activity in T. grandis kernels. However, the flavonoid compositions which mainly contribute to antioxidant capacity and the molecular basis of flavonoid biosynthesis in T. grandis remain unclear. Here, transcriptome sequencing and metabolomics analysis for kernels were performed. In total, 124 flavonoids were identified. Among them, 9 flavonoids were highly correlated with antioxidant activity. Furthermore, unigenes encoding CHS, DFR and ANS showed significant correlation with the 9 flavonoids. Transient overexpression of TgDFR1 in tobacco leaves resulted in increased antioxidant activity. Moreover, several transcription factors from MYB, bHLH and bZIP families were identified by co-expression assay, suggesting that they may regulate flavonoid biosynthesis. Our findings provide a molecular basis and new insights into the flavonoid biosynthesis in T. grandis kernels.

Metabolite Profiling and Transcriptome Analysis Explains Difference in Accumulation of Bioactive Constituents in Licorice (Glycyrrhiza uralensis) Under Salt Stress

Salinity stress significantly affects the contents of bioactive constituents in licorice Glycyrrhiza uralensis. To elucidate the molecular mechanism underlying the difference in the accumulation of these constituents under sodium chloride (NaCl, salt) stress, licorice seedlings were treated with NaCl and then subjected to an integrated transcriptomic and metabolite profiling analysis. The transcriptomic analysis results identified 3,664 differentially expressed genes (DEGs) including transcription factor family MYB and basic helix-loop-helix (bHLH). Most DEGs were involved in flavonoid and terpenoid biosynthesis pathways. In addition, 121 compounds including a triterpenoid and five classes of flavonoids (isoflavone, flavone, flavanone, isoflavan, and chalcone) were identified, and their relative levels were compared between the stressed and control groups using data from the ultrafast liquid chromatography (UFLC)-triple quadrupole-time of flight-tandem mass spectrometry (TOF-MS/MS) analysis. Putative biosynthesis networks of the flavonoids and triterpenoids were created and combined with structural DEGs such as phenylalanine ammonia-lyase (PAL), 4-coumarate-CoA ligase [4CL], cinnamate 4-hydroxylase [C4H], chalcone synthase [CHS], chalcone-flavanone isomerase [CHI], and flavonoid-3',5' hydroxylase (F3',5'H) for flavonoids, and CYP88D6 and CYP72A154 for glycyrrhizin biosynthesis. Notably, significant upregulation of UDP-glycosyltransferase genes (UGT) in salt-stressed licorice indicated that postmodification of glycosyltransferase may participate in downstream biosynthesis of flavonoid glycosides and triterpenoid saponins. Accordingly, the expression trend of the DEGs is positively correlated with the accumulation of glycosides. Our study findings indicate that key DEGs and crucial UGT genes co-regulate flavonoid and saponin biosynthesis in licorice under salt stress.

Metabolomics and Transcriptomics Provide Insights into Anthocyanin Biosynthesis in the Developing Grains of Purple Wheat (Triticum aestivum L.)

Purple wheat is thought to have beneficial effects on humans owing to its high anthocyanin content. However, a systematic understanding of the anthocyanin biosynthesis process in developing wheat grain is lacking. Here, the dynamic changes in anthocyanin components and transcripts in the grain of purple wheat ZNM168 at five developmental stages (10, 15, 20, 25, and 30 DAF) were characterized. Compared with other anthocyanins, four components, cyanidin 3-O-rutinoside, cyanidin 3-O-glucoside, cyanidin 3,5-O-diglucoside, and malvidin 3-O-glucoside, were significantly accumulated with grain development. In particular, the considerable accumulation of cyanidin 3-O-rutinoside indicated that it was the pivotal pigment for the purple grain. Transcriptome analysis revealed that the nine differentially expressed genes related to anthocyanin biosynthesis belonged to the BZ1 group, the homologous enzyme encoded by the maize Bronze-1 locus, which may primarily serve to glucosylate anthocyanidins. By constructing a gene coexpression network based on weighted gene coexpression network analysis, the TaBZ1 UniGene (TraesCS1D02G019200) was predicted as a core gene in anthocyanin biosynthesis. In addition, correlation analysis between the metabolites and transcripts suggested that TraesCS2A01G527700 (TaCHS) and TraesCS6B01G006200 (TaANS) were considered critical structural genes in the anthocyanin biosynthesis pathway. This study provides insights to exploit genes pinpointed as genetic engineering targets, thereby breeding anthocyanin-enriched wheat.

Genome-wide analysis of BBX gene family in Tartary buckwheat (Fagopyrum tataricum)

BBX (B-box), a zinc finger transcription factor with one or two B-box domains, plays an important role in plant photomorphogenesis, growth, and development as well as response to environmental changes. In this study, 28 Tartary buckwheat BBX (FtBBX) genes were identified and screened using a comparison program. Their physicochemical properties, gene structures, conserved motifs, distribution in chromosomal, and phylogeny of the coding proteins, as well as their expression patterns, were analyzed. In addition, multiple collinearity analysis in three monocots and three dicot species illustrated that the BBX proteins identified from monocots clustered separately from those of dicots. Moreover, the expression of 11 candidate BBX genes with probable involvement in the regulation of anthocyanin biosynthesis was analyzed in the sprouts of Tartary buckwheat during light treatment. The results of gene structure analysis showed that all the 28 BBX genes contained B-box domain, three genes lacked introns, and these genes were unevenly distributed on the other seven chromosomes except for chromosome 6. The 28 proteins contained 10 conserved motifs and could be divided into five subfamilies. BBX genes of Tartary buckwheat showed varying expression under different conditions demonstrating that FtBBXs might play important roles in Tartary buckwheat growth and development. This study lays a foundation for further understanding of Tartary buckwheat BBX genes and their functions in growth and development as well as regulation of pigmentation in Tartary buckwheat.

Joint Transcriptomic and Metabolic Analysis of Flavonoids in Cyclocarya paliurus Leaves

Flavonoids are a class of commonly occurring natural compounds in the plant kingdom with various biological activities. This study compares the content of flavonoids in Cyclocarya paliurus at different developmental stages to better inform the selection of the optimal picking period. Thus, we analyzed the transcriptome and metabolome of C. paliurus at different developmental stages. The transcriptome analysis revealed 44 genes involved in the biosynthesis of flavonoids in C. paliurus, with 10 differentially expressed genes across the four different developmental stages. The metabolites were separated and identified by a combination of chromatography and mass spectrometry, followed by multi-reaction monitoring mode analysis of triple quadrupole mass spectrometry for complete metabolite quantification. In the flavonoid synthesis pathway, a total of 137 differential flavonoids were detected. The joint transcriptome and metabolome analysis showed that the expression trends in differential metabolites and genes were significantly related. Four MYB transcription factors and two bHLH transcription factors that are closely related to flavonoid biosynthesis were identified. The regulation network of flavonoid biosynthesis in C. paliurus was thus established, providing guidance for follow-up research.

Integrated microRNA and transcriptome profiling reveal key miRNA-mRNA interaction pairs associated with seed development in Tartary buckwheat (Fagopyrum tataricum)

BACKGROUND

Tartary buckwheat seed development is an extremely complex process involving many gene regulatory pathways. MicroRNAs (miRNAs) have been identified as the important negative regulators of gene expression and performed crucial regulatory roles in various plant biological processes. However, whether miRNAs participate in Tartary buckwheat seed development remains unexplored.

RESULTS

In this study, we first identified 26 miRNA biosynthesis genes in the Tartary buckwheat genome and described their phylogeny and expression profiling. Then we performed small RNA (sRNA) sequencing for Tartary buckwheat seeds at three developmental stages to identify the miRNAs associated with seed development. In total, 230 miRNAs, including 101 conserved and 129 novel miRNAs, were first identified in Tartary buckwheat, and 3268 target genes were successfully predicted. Among these miRNAs, 76 exhibited differential expression during seed development, and 1534 target genes which correspond to 74 differentially expressed miRNAs (DEMs) were identified. Based on integrated analysis of DEMs and their targets expression, 65 miRNA-mRNA interaction pairs (25 DEMs corresponding to 65 target genes) were identified that exhibited significantly opposite expression during Tartary buckwheat seed development, and 6 of the miRNA-mRNA pairs were further verified by quantitative real-time polymerase chain reaction (qRT-PCR) and ligase-mediated rapid amplification of 5' cDNA ends (5'-RLM-RACE). Functional annotation of the 65 target mRNAs showed that 56 miRNA-mRNA interaction pairs major involved in cell differentiation and proliferation, cell elongation, hormones response, organogenesis, embryo and endosperm development, seed size, mineral elements transport, and flavonoid biosynthesis, which indicated that they are the key miRNA-mRNA pairs for Tartary buckwheat seed development.

CONCLUSIONS

Our findings provided insights for the first time into miRNA-mediated regulatory pathways in Tartary buckwheat seed development and suggested that miRNAs play important role in Tartary buckwheat seed development. These findings will be help to study the roles and regulatory mechanism of miRNAs in Tartary buckwheat seed development.

Mapping QTLs for 1000-grain weight and genes controlling hull type using SNP marker in Tartary buckwheat (Fagopyrum tataricum)

BACKGROUND

Tartary buckwheat (Fagopyrum tataricum), an important pseudocereal crop, has high economic value due to its nutritional and medicinal properties. However, dehulling of Tartary buckwheat is difficult owing to its thick and tough hull, which has greatly limited the development of the Tartary buckwheat processing industry. The construction of high-resolution genetic maps serves as a basis for identifying quantitative trait loci (QTLs) and qualitative trait genes for agronomic traits. In this study, a recombinant inbred lines (XJ-RILs) population derived from a cross between the easily dehulled Rice-Tartary type and Tartary buckwheat type was genotyped using restriction site-associated DNA (RAD) sequencing to construct a high-density SNP genetic map. Furthermore, QTLs for 1000-grain weight (TGW) and genes controlling hull type were mapped in multiple environments.

RESULTS

In total, 4151 bin markers comprising 122,185 SNPs were used to construct the genetic linkage map. The map consisted of 8 linkage groups and covered 1444.15 cM, with an average distance of 0.35 cM between adjacent bin markers. Nine QTLs for TGW were detected and distributed on four loci on chromosome 1 and 4. A major locus detected in all three trials was mapped in 38.2-39.8 cM region on chromosome 1, with an LOD score of 18.1-37.0, and explained for 23.6-47.5% of the phenotypic variation. The genes controlling hull type were mapped to chromosome 1 between marker Block330 and Block331, which was closely followed by the major locus for TGW. The expression levels of the seven candidate genes controlling hull type present in the region between Block330 and Block336 was low during grain development, and no significant difference was observed between the parental lines. Six non-synonymous coding SNPs were found between the two parents in the region.

CONCLUSIONS

We constructed a high-density SNP genetic map for the first time in Tartary buckwheat. The mapped major loci controlling TGW and hull type will be valuable for gene cloning and revealing the mechanism underlying grain development and easy dehulling, and marker-assisted selection in Tartary buckwheat.

Widely Targeted Metabolomics Analysis of Different Parts of Salsola collina Pall

Salsola collina Pall has a long history of being used as a traditional medicine to treat hypertension, headache, insomnia, constipation and vertigo. However, only a few biologically active substances have been identified from S. collina. Here, the shoots and roots of S. collina, namely L-Sc and R-Sc, were studied. The primary and secondary metabolites were investigated using ultrahigh-performance liquid chromatography-electrospray ionization-tandem mass spectrometry (UPLC-ESI-MS/MS). A total of 637 putative metabolites were identified and these metabolites were mainly classified into ten different categories. Correlation analysis, hierarchical clustering analysis, principal component analysis and orthogonal partial least squares discriminant analysis of metabolites showed that the L-Sc samples could be clearly separated from the R-Sc samples. Differential accumulated metabolite analysis revealed that most of differential primary metabolites were significantly lower in the L-Sc than in the R-Sc. Conversely, the major differential secondary metabolites had higher levels in the L-Sc than in the R-Sc. Further analysis indicated that the flavonoids were the major putative antioxidant components and most of putative antioxidant components exhibited higher relative concentrations in the L-Sc than the R-Sc. These results improve our understanding of metabolite accumulation and provide a reference for the study of medicinal value in S. collina.

Metabolomics Analyses of Cotyledon and Plumule Showing the Potential Domestic Selection in Lotus Breeding

Lotus (Nelumbo nucifera) seeds are widely consumed as functional food or herbal medicine, of which cotyledon (CL) is the main edible part, and lotus plumule (LP) is commonly utilized in traditional Chinese medicine. However, few studies have been conducted to investigate the chemical components of CL and LP in dry lotus seeds, not to mention the comparison between wild and domesticated varieties. In this study, a widely targeted metabolomics approach based on Ultra Performance Liquid Chromatography-electrospray ionization-Tandem mass spectrometry (UPLC-ESI-MS/MS) was utilized to analyze the metabolites in CL and LP of China Antique (“CA”, a wild variety) and Jianxuan-17 (“JX”, a popular cultivar). A total of 402 metabolites were identified, which included flavonoids (23.08% to 27.84%), amino acids and derivatives (14.18–16.57%), phenolic acids (11.49–12.63%), and lipids (9.14–10.95%). These metabolites were classified into ten clusters based on their organ or cultivar-specific characters. Most of these metabolites were more abundant in LP than in CL for both varieties, except for metabolites belonging to organic acids and lipids. The analysis of differentially accumulated metabolites (DAMs) demonstrated that more than 25% of metabolites detected in our study were DAMs in CL and LP comparing “JX” with “CA”, most of which were less abundant in “JX”, including 35 flavonoids in LP, 23 amino acids and derivatives in CL, 7 alkaloids in CL, and 10 nucleotides and derivatives in LP, whereas all of 11 differentially accumulated lipids in LP were more abundant in “JX”. Together with the fact that the seed yield of “JX” is much higher than that of “CA”, these results indicated that abundant metabolites, especially the functional secondary metabolites (mainly flavonoids and alkaloids), were lost during the process of breeding selection.

Transcriptomic and Metabolomic Analyses Reveals That Exogenous Methyl Jasmonate Regulates Galanthamine Biosynthesis in Lycoris longituba Seedlings

The Amaryllidaceae alkaloid galanthamine (Gal) in Lycoris longituba is a secondary metabolite that has been used to treat Alzheimer's disease. Plant secondary metabolism is affected by methyl jasmonate (MeJA) exposure, although the regulatory mechanisms of MeJA on L. longituba seedlings remains largely unknown. In the present study, 75, 150, and 300 μM MeJA were used as treatments on L. longituba seedlings for 7, 14, 21, and 28 days, while 0 μM MeJA was used as the control (MJ-0). The effect of exogenous MeJA on Gal synthesis in L. longituba was then investigated using transcriptomic sequencing and metabolite profiling via GC-MS and LC-MS analysis. Galanthamine (Gal), lycorine (Lyc), and lycoramine (Lycm) abundances were 2. 71-, 2. 01-, and 2.85-fold higher in 75 μM MeJA (MJ-75) treatment plants compared to MJ-0 treatment plants after 7 days of cultivation. Transcriptomic analysis further showed that MJ-75 treatment significantly induced the expression of norbelladine synthase (NBS) and norbelladine 4'-O-methyltransferase (OMT), which are involved in the Gal biosynthesis pathway. In addition, increased expression was observed in MJ-75 treatment plants for genes in the JA synthesis and JA signaling pathways including those of allene oxide cyclase (AOC), 12-oxo-phytodienoic acid reductase (OPR), jasmonic acid amino acid synthase (JAR), and transcription factor MYC. The L. longituba tyrosine decarboxylase (LlTYDC) enzyme was identified and proposed to be involved in the Gal biosynthetic pathway. Metabolomics results demonstrated that the accumulation of Amaryllidaceae alkaloids, and especially alkaloids in the Gal biosynthesis pathway, could be induced by MJ-75 treatment. Interestingly, metabolites in the JA synthesis pathway were also affected by MeJA treatment. Overall, this multi-omics study suggests that both the JA synthesis/JA signaling and Gal biosynthesis pathways were affected by exogenous MeJA treatment. This comprehensive study of gene expression and metabolite contents can help us better understand the molecular mechanisms underlying MeJA-mediated Gal biosynthesis in L. longituba.

Seeing the forest for the trees: Retrieving plant secondary biochemical pathways from metabolome networks

Over the last decade, a giant leap forward has been made in resolving the main bottleneck in metabolomics, i.e., the structural characterization of the many unknowns. This has led to the next challenge in this research field: retrieving biochemical pathway information from the various types of networks that can be constructed from metabolome data. Searching putative biochemical pathways, referred to as biotransformation paths, is complicated because several flaws occur during the construction of metabolome networks. Multiple network analysis tools have been developed to deal with these flaws, while in silico retrosynthesis is appearing as an alternative approach. In this review, the different types of metabolome networks, their flaws, and the various tools to trace these biotransformation paths are discussed.