Next Article in Journal
Korean Traditional Medicine (Jakyakgamcho-tang) Ameliorates Colitis by Regulating Gut Microbiota
Previous Article in Journal
Metabolomics Approach Reveals the Effects of Breed and Feed on the Composition of Chicken Eggs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Transcriptome and Metabolic Profiling Analysis of Buckwheat (Fagopyrum Tataricum (L.) Gaertn.) under Salinity Stress

1
Laboratory of Natural Medicine and Molecular Engineering, Department of Medicinal Plant, College of Plant Science and Technology, Huazhong Agriculture University, Wuhan 430070, China
2
Division of Life Sciences and Bio-Resource and Environmental Center, Incheon National University, Yeonsu-gu, Incheon 406-772, Korea
3
Key Laboratory of Environment Correlative Dietology (Ministry of Education), College of Food Science and Technology, Huazhong Agricultural University, Wuhan 430070, China
4
Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan 430070, China
*
Author to whom correspondence should be addressed.
Submission received: 12 September 2019 / Revised: 9 October 2019 / Accepted: 11 October 2019 / Published: 14 October 2019

Abstract

:
Tartary buckwheat (Fagopyrum tataricum (L.) Gaertn.) is a nutritional crop, which has high flavonoid content. However, buckwheat is a salt sensitive glycophyte cereal crop and the growth and grain yield of buckwheat are significantly affected by soil salinity. In this study, we performed a comprehensive analysis of the transcriptome and metabolome of salt treated-buckwheat to understand the effects of salinity on buckwheat. A total of 50,681,938 clean reads were acquired from all samples. We acquired 94,950 unigenes with a mean length of 1133 bp and N50 length of 1900 bp assembly. Of these, 63,305 unigenes (66.7%) were matched in public databases. Comparison of the transcriptome expression patterns between control and salt treated groups showed that 4098 unigenes were up-regulated and 3292 unigenes were down-regulated significantly. Further, we found that genes involved with amino acid, lipid and nucleotide metabolism were most responsive to salt stress. Additionally, many genes involved in secondary metabolite biosynthesis changed significantly following treatment. Those affected included phenylpropanoid biosynthesis and flavonoid biosynthesis. Chromatographic analysis was used to examine the differences in concentration of flavonoids, carotenoids, amino acids and organic acids in the samples following treatment. There was a significant increase in rutin (12.115 mg/g dry weight), following salt stress; whereas, six carotenoids (lutein, zeaxanthin, 13Z-β-carotene, α-carotene, E-β-carotene and 9Z-β-carotene) did not significantly respond to salt stress. Ultimately, our data acts as a valuable resource for future research on buckwheat and can be used as the basis for future analysis focused on gene-to-metabolite networks in buckwheat.

Graphical Abstract

1. Introduction

Salinity is one of the environmental factors that has the greatest effect on food production and quality. High salinity can limit the growth, development and yield of crops [1]. There are three salt-induced stress pathways—osmotic stress, ionic stress and secondary stress. High salinity can change salt and ion concentrations, with toxic effects and alterations to important metabolic pathways [2]. Most crops are salt sensitive, high Na+ concentration in the cytosol inhibits many metabolic and physiological processes [3]. High concentration of Na+ ions can also cause secondary stress in plant, including the accumulation of toxic compounds and disruption of nutrient balances [2]. Systematic research of salt-tolerance mechanisms has been reported in Arabidopsis, canola, rice and several others [4,5,6,7]. Previous studies have shown that canola responds to high salinity at the morphological, physiological and biochemical levels along with changes at the molecular level [4]. Secondary metabolites play an important role in the response to biotic and abiotic stresses in plants [8,9,10]. Many reports have shown that the synthesis of secondary metabolites, including carotenoids, anthocyanins and flavonoids, is heavily influenced by salt stress [11,12], these compounds could serve as the nonenzymatic scavengers to reduce reactive oxygen species (ROS) damage and protect cellular structures and macromolecules in plant [5,13,14].
Transcription factors play an important role in plant development and response to environment stress [15,16]. There are many transcription factors has been reported and involved in stress-responsive in plant. The transcription factor TaMYB73 has been reported to as responding salinity stress improving ionic resistance partly via the regulation of stress-responsive genes in wheat [17]. AtMYB20 was shown to negatively regulate type 2C serine/threonine protein phosphatases resulting in enhanced salt tolerance in Arabidopsis [18]. Zhong et al. reported two NAC transcription factors involved in response to high-salinity in Brassica napus [19]. Similar studies on the regulation of transcription factors in response to biotic or abiotic stress, including NAC, WRKY and bZIP amongst others, were also reported in rice, wheat, tomato and other plants [20,21,22,23,24].
Tartary buckwheat (Fagopyrum tataricum (L.) Gaertn., Polygonaceae family) is both a food and medicinal plant. It is predominantly grown in the semi-arid and arid regions of the mountainous areas of southwest China, India, Nepal, Japan and Korea [25,26]. The chemical composition of tartary buckwheat was comprehensively investigated, this included analysis of its flavonoid [27], phenylpropanoid glycosides [28], organic acids [29], production of various other important compounds [30,31]. In buckwheat, the major flavonoid is rutin, whose concentration can reach 81 mg/g in the groats of tartary buckwheat [32]. Rutin-rich buckwheat is known for reducing cholesterol, blood clots and high blood pressure [33,34,35,36], in addition to showing anti-oedema effects, antioxidant activity [37,38,39,40], anti-inflammatory activity [41] and antifatigue properties [42]. In addition, some proteins in buckwheat are believed to have anti-tumour properties [43,44,45].
The synthesis of secondary metabolites in buckwheat is strongly influenced by environmental factors. It has been reported that the content of flavonoids was mainly up-regulated and accumulated in buckwheat sprouts when treated with exogenous methyl jasmonate (MeJA) [46,47]. The accumulation of phenolic compounds and anthocyanins were obviously increased with the treatment of L-phenylalanine and light-emitting diodes (LEDs) lights in tartary buckwheat sprouts [48]. Buckwheat is a salt sensitive glycophyte cereal crop and the growth and grain yields are significantly affected by soil salinity [49,50]. In addition, the metabolic synthesis of buckwheat is also affected by salt stress. Gamma-aminobutyric acid (GABA) and flavonoid production are reported to alter when buckwheat sprouts are treated to high salinity conditions [51,52]. Two flavonol compounds, quercetin and kaempferol, were up-regulated in buckwheat cultivar “Hokkai T10” when treated with NaCl [53]. However, the comprehensive effect of salinity on transcriptional and metabolic activities in buckwheat has not been widely reported. [54].
In order to provide a global view of salinity-tolerant regulatory mechanisms in buckwheat, in this study, the influence of salt-induced salinity stress was investigated in buckwheat using transcriptome and metabolite profiling analysis. To do this at a molecular level, the Illumina HiSeq™ 2000 platform was employed to screen differential gene expression under salt stress. This data were then analyzed and the main metabolites, especially, flavonoids and carotenoids, were detected using chromatographic analysis. Our research provided foundational information which may lead to a better understanding of the molecular mechanisms of salt-resistance in buckwheat.

2. Results

2.1. De Novo Transcriptome Assembly of F. Tataricum and Sequence Analysis

We performed RNA-sequencing analysis using total RNA from plants exposed to normal and high salinity conditions. After removal of adaptor sequences and reads with low quality, a total of 50,681,938 clean reads were acquired from all samples. Subsequently, after assembling of these clean reads, 94,950 unigenes, with the average length of 1133 nucleotides (nt) and the N50 is 1900 nt, have been identified between control (C) and salt-treated (T) samples (Table 1). Among them, 42,127 unigenes have a length between 300 nt and 1000 nt and 15,263 unigenes are longer than 1000 nt (as shown in Supplementary Figure S1).

2.2. Functional Annotation and Classification of the Buckwheat Transcriptome

To obtain gene function information, the assembled transcripts were subjected to a search through six public databases, including the NCBI non-redundant protein sequences database (Nr), the NCBI non-redundant nucleotide database (Nt), the Clusters of Orthologous Groups of protein databases (COG), the Swiss-Prot protein database (Swiss-Prot), the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Finally, 63,305 unigenes (66.67% of the assembled transcripts) were well annotated in at least one of the public databases (Supplementary Table S1, Figure S2). In detail, 60,148 unigenes were found in the Nr database, 52,067 in the Nt database, 43,129 in the Swiss-Prot database, 38,296 in the KEGG database and 28,107 in the COG database (Supplementary Table S1). Subsequently, using E-value frequency distribution analysis, 65.62% of the Nr annotated unigenes had an e-value < 1.0e−30 (Figure 1A). While sequence similarity distribution indicated that 48.50% of the unigenes showed more than 60% similarity with matched genes (Figure 1B). The species homologous distribution analysis of unigenes is shown in Figure 1C. These unigenes matched with the sequences from Vitis vinifera (21,470, 35.70%), Arabidopsis thaliana (8175, 13.59%), Glycine max (3.978, 6.61%), Medicago truncatula (3238, 5.38%) and Sorghum bicolor (2877, 4.78%) Figure 1C). The function of the unigenes from F. tataricum were divided into 3 major categories—‘biological process,’ ‘cellular component’ and ‘molecular function’ according to Gene Ontology (GO) analysis. Among the 57 sub-categories, the highest number of unigenes were putatively ascribed to involvement in the ‘cell’ (45,305), ‘cell part’ (45,017) and ‘organelle’ (39,256) under the category of cellular components. In the biological process category, the following were most strongly represented, ‘cellular process’ (40,743), ‘metabolic process’ (39,980) and ‘single-organism process’ (37,819) In the molecular function category enrichment for ‘binding’ (29,090) and ‘catalytic activity’ (28,351) was observed (Supplementary Figure S3). Classification of orthologous group (COG) analysis showed that 28,107 of the total transcripts could be placed in 25 categories. The top homologous categories were the cluster for ‘General function prediction only’ followed by ‘Transcription’ and ‘Replication, recombination and repair’ (Supplementary Figure S4).

2.3. Differentially Expressed Genes (DEGs) of Buckwheat in Response to Salt Stress

The differentially expressed genes between salt treated and non-treated samples of buckwheat were analyzed. The results showed that significantly more transcripts (5350) were significantly up-regulated in salt treated samples than down-regulated (4207) when compared with control (p-value < 0.001 and |log2FoldChange| > 2). (Supplementary Figure S5). These differentially expressed genes were annotated using GO analysis and the results indicated that top clusters of differentially expressed genes (DEGs) belong to ‘cellular process,’ ‘single-organism process’ and ‘metabolic process’ in the biological process category. They fell into ‘cell,’ ‘cell part’ and ‘organelle’ in the cellular component category and ‘catalytic activity’ and ‘binding’ in the molecular function category (Figure 2). According to pathway annotation, the DEGs mainly fell into ‘Metabolic pathway’ (1414, 33.26%), ‘Biosynthesis of secondary metabolites’ (897, 21.1%) and ‘plant hormone signal transduction’ (419, 9.86%) ((Figure 3, Supplementary Table S2).
According to GO analysis, 39,980 unigenes were linked to ‘metabolic processes,’ of which genes in ‘carbohydrates metabolism,’ ‘amino acid metabolism,’ ‘lipid and nucleotide metabolism’ were significantly differentially regulated under stress (Supplementary Figures S6 and S7, Table S2). In detail, 536 unigenes were up-regulated and 349 unigenes were down-regulated in carbohydrate metabolism under salt-stress (Supplementary Figure S6). Further, 167 unigenes from the starch and sucrose metabolism pathways were found to strongly respond to salt stress. The unigenes involved in glycolysis/gluconeogenesis, amino sugar and nucleotide metabolism, pyruvate metabolism, citrate cycle (TCA cycle), galactose metabolism, fructose and mannose metabolism were also predominantly up and not down regulated following salt stress (Supplementary Figure S6). While the unigenes involved in amino acid metabolism were largely up-regulated, especially in cysteine and methionine metabolism and the valine, leucine and isoleucine degradation process (Supplementary Figure S7). Secondary metabolites play an important role in stress resistance in plants. In our study, differentially expressed genes involved in 14 biosynthetic pathways were identified from the transcriptome data (Figure 4). The top salt-response genes belong to phenylpropanoid biosynthesis and flavonoid biosynthesis. In all, our results indicate that primary and secondary metabolic biosynthetic pathway were directly altered in response to salt-stress in buckwheat.

2.4. The Content of Metabolites in Tartary Buckwheat in Response to Salt Stress

Primary metabolites, such as carbohydrates, amino acids, lipids and nucleotides, are all predominantly involved in physiological processes in living organisms. The comprehensive metabolic profiling of the primary metabolites was detected in tartary buckwheat in response to salt stress treatment. In our study, 39 hydrophilic metabolites and 14 liphophilic metabolites were identified as differentially regulated in buckwheat using GC-TOFMS, these included amino acids, organic acids and carbohydrates (Supplementary Tables S3 and S4). The sugars, including fructose, galactose, glucose, xylose and sucrose, content was increased in most of the salt-stress samples when compared with control. The level of galactose was 4.38- and 2.67-fold accumulated in the first two days treated samples, while, the content of glucose and fructose were 3.01-fold and 3.07-fold increased after the two days of treatment (Figure 5). Organic acids and amino acids including pyruvic acid, pyroglutamic acid, succinic acid and methionine, also showed a similar trend (Figure 5), indicating that they may contribute to the differential accumulation of secondary metabolites between salt-treated and control buckwheat samples. The concentration of lipophilic metabolites, including isoprenoid compounds α-tocopherol and β-sitosterol, had minor fluctuations, while no significant change was detected (Supplementary Tables S3 and S4).
Correlation analysis is a statistical method used to analyze the strength of the relationship between two quantitative samples; it can be applied to demonstrate the associations between metabolites belonging to specific biological systems [55]. To examine the relationship between lipophilic metabolites and hydrophilic metabolites in buckwheat, Pearson’s correlation analysis was performed (Figure 6). The most distinct relationships were observed between sugars and their derivatives (galactose, inositol, fructose, glucose, mannitol, xylose, sucrose), organic acids and the products of the tricarboxylic acid cycle (fumaric acid, malic acid, citric acid and succinic acid) and amino acids and their derivatives (β-alanine, phenylalanine, methionine, 4-aminobutyric acid).

2.5. The Transcriptional and Metabolic Regulation of Phenylpropanoid Biosynthetic Pathway of Buckwheat in Response to Salt Stress

A total of 122 unigenes from phenylpropanoid biosynthesis were significantly differentially expressed between control and salt-treated (p < 0.01), of these, 63 unigenes were up-regulated and 59 were down-regulated (Supplementary Table S5). Interestingly, according to transcriptome data, a total of 92 unigenes from flavonoid biosynthesis were significantly differentially expressed under salt stress (p < 0.01), 61 of these unigenes were significantly up-regulated, while, only 31 were down-regulated (Supplementary Table S5). More unigenes in flavone and flavonol biosynthesis, isoflavonoid biosynthesis and anthocyanin biosynthesis pathway were down-regulated after salt stress when compared with the control (Figure 4, Supplementary Table S5). In this study, the gene expression profiles for the flavonoid biosynthesis pathway was established using quantification real-time PCR (qRT-PCR) (Figure 7). Our results showed that the main genes involved in flavonoid biosynthesis, including the genes from FtPAL to FtRT, were down-regulated after the first day of treatment and then significantly up-regulated in samples taken five days post treatment (Figure 7).
At the same time, we analyzed the main metabolites in buckwheat using high-performance liquid chromatography (HPLC). In total, eight main metabolites from flavonoid biosynthesis pathway, including catechin, chlorogenic acid, caffeic acid, p-coumaric acid, benzoic acid, rutin, quercetin, kaempferol, were detected and are shown in Table 2. The major metabolite in buckwheat, rutin, was significantly up-regulated in salt-treated samples, the majority of the rutin was accumulated during day one of treatment, where an increase of 12.24% was recorded when compare with the control. Rutin concentration initially decreased and then underwent significant upregulation in the following five-days post treatment, which was in agreement with the transcriptional data collected. Benzoic acid also showed a concentration spike at one day post treatment, while, catechin and quercetin were most enriched at four-days post treatment. Quercetin was 10-fold higher at four days post treatment when compared with the control. Other compounds did not show significant changes after salt treatment.

2.6. The Transcriptional and Metabolic Regulation of Carotenoid Biosynthesis of Buckwheat in Response to Salt Stress

Plant secondary metabolites, specifically carotenoids and phenolics, are stress inducible. These compounds play important roles in growth and contribute to the nutritional value of vegetables and crops [8]. Our transcriptome data shows that, 23 genes from the carotenoid biosynthesis pathway were up-regulated and 23 were down-regulated (Supplementary Table S6). As for the metabolites, 6 carotenoids including lutein, zeaxanthin, 13Z-β-carotene, α-carotene, E-β-carotene and 9Z-β-carotene were detected in this study, our data indicate that the predominant composition of carotenoids in tartary buckwheat seedlings were lutein and E-β-carotene, which was consistent with previous reports [56]. The accumulation of carotenoids changed with seedling age and the content of lutein, E-β-carotene and 13Z-β-carotene reached their highest level at four-days control group rather than salt-treated group. In detail, 299.29 μg/g and 249.16 μg/g of lutein, 235.54 μg/g and 164.47 μg/g of E-β- carotene, along with 47.27 μg/g and 37.24 μg/g of 13Z-β-carotene were detected at four days post treatment in the control and salt stress groups respectively. However, only zeaxanthin was found to increase significantly in four-days salt-treated buckwheat. Other carotenoids also did not show any significant accumulation following salt stress (Table 3).

3. Discussion

Plants have evolved sophisticated physiological adaptations to help them cope with various environmental stressors, including increased salinity [57]. This stress resistance is usually accomplished through various primary and secondary metabolites [8,10,12]. In previous reports, Jeon et al., reported that some organic acids, sugars and their derivatives, anthocyanins and proanthocyanidins were significantly increased in tartary buckwheat seedlings by cold treatment [58]. Anthocyanins accumulation was primarily detected in the epidermal and cortex cells of hypocotyls in tartary buckwheat sprout in response to cold stress [59]. According to metabolic analysis of wheat, five and eleven metabolites involved in glycolysis, TCA cycle and shikimate pathway, were significantly increased or decreased in the leaves of wheat seedlings during saline stress [60]. Various metabolites, including ascorbic acid, alkaloids, carotenoids, flavonoids, phenolics, tocopherol amongst others, have been reported as having expression profiles that are closely correlated with various stressors in plants, these metabolites mostly function as the nonenzymatic scavengers to protect plants from salt stress [9,12,61,62,63].
Transcriptomic analysis using next-generation sequencing (RNA-seq) has been successfully and widely applied to investigate the gene expression of plants and other organisms [64]. This technology has been used to elucidate the stress-response mechanisms in plants at the transcriptional level [2,4,65]. Lu et al. [54] reported that there was 385 DEGs with a fold change ≥ 2 and a false discovery rate (FDR) of <  0.05 in 15-days-old-seedlings in common buckwheat in response to salt stress (100 mM NaCl, 72 h). DEGs were highly enriched for ‘Protein processing in endoplasmic reticulum,’ ‘Plant hormone signal transduction,’ ‘Carbon metabolism,’ ‘Pentose phosphate pathway,’ ‘Biosynthesis of amino acids’ and ‘Citrate cycle (TCA cycle)’ [54]. Even the nutritution and biological characteristics are different from common buckwheat, the similar result was also detected in this tartary buckwheat research and the DEGs from control versus 100 mM NaCl treated groups mainly belong to ‘Metabolic pathway,’ ‘Biosynthesis of secondary metabolites,’ ‘Plant hormone signal transduction,’ ‘Starch and sucrose metabolism,’ ‘Amino sugar and nucleotide sugar metabolism’ and so on according to pathway annotation (Figure 3, Table S2); this result indicated that the primary and secondary metabolites biosynthetic pathway were obviously activated at the transcriptional level in buckwheat in response to salt stress.
According to the transcriptome data, there were a number of transcripts (5350 up- and 4207 down-regulated) differentially regulated following salt treatment (p-value < 0.001 and |log2FoldChange| > 2) (Figure S5). A total of 683 unigenes were annotated, 432 were up-regulated and 251 down-regulated in buckwheat during salt stress. These results reveal that more genes are involved in buckwheat in response to salt stress than previous studies had identified [66]. In the previous study 14-day-old Tartary buckwheat seedlings (Chuanqiao No.1) were exposed to 200 mM NaCl for 24 hours. Here they identified 455 DEGs, these DEGs were mainly involved the biochemical pathways including ‘Antigen processing and presentation,’ ‘Protein processing in endoplasmic reticulum,’ ‘Estrogen signaling pathway’ and ‘Plant-pathogen interaction’ when analyzed using KEGG enrichment analysis [66]. In contrast, the DEGs identified in this study predominantly featured in the following pathways ‘Metabolic pathways’ (33.26%), ‘Biosynthesis of secondary metabolites’ (21.1%) and ‘plant hormone signal transduction’ (9.86%) (Figure 3, Table S2). These differences may be the result of differences in the cultivars, seedling age of buckwheat and the severity and timing of the salt stress [66]. According to previous reports, the accumulation of primary or secondary metabolites induced by salinity stress may vary and base on different species, cultivars or treatment way [10,11,53]. It has been reported that the phenolic acids are species-specific among Brassicaceae and some phenolic acids may especially participate in stress tolerance in more salt-tolerant varieties or species, such as white cabbage and kale [67]. Interestingly, the NaCl treatment lead to different increased degree of carotenoids in one nonanthocyanin-accumulating and three anthocyanin-accumulating genotypes tomato, however, the accumulation of total anthocyanins was opposite in different tomato cultivars [11]. The result also suggested that a considerable number of genes and metabolites were differentially regulated in salt-tolerant or not tolerant varieties or species. This further indicated that it is important to study and exploit the viability of saline soil suitable buckwheat species or cultivars with increased levels of functional compounds.
Combining transcriptomic and metabolic data is one approach to understanding the complex biological processes involved in plants [68,69,70,71]. Using transcriptomic and metabolic analysis, Jeon et al., reported that the metabolites including some organic acids, sugars and their derivatives, anthocyanins and proanthocyanidins were significantly increased and some amino acids and their derivatives were decreased in tartary buckwheat seedlings by cold treatment and the transcription levels of genes in phenylpropanoid biosynthetic pathway were correspondingly accumulated in response to cold stress [58]. Even some researchers reported that salinity effected the accumulation of phenolic compounds, carotenoids and flavonoids in buckwheat and the salt-related genes were analyzed by de novo assembly method [54,66], however, the interactive effects of salinity on both transcriptional and metabolic regulation of buckwheat was unclear [53,54,66,72]. In this study, we investigated the impact of salt stress on the transcriptional and metabolism response of buckwheat through combining transcriptomic and metabolic profiling approach. The content of primary and secondary metabolites, including carbohydrates, organic acids or amino acids, flavonoids, carotenoids and tocopherols were detected, while the accumulations of above metabolites were differentially modulated in buckwheat in response to salt stress treatment.
According to our results, most of the sugars, organic acids and amino acids increased during salt-stress, which was a little different from the accumulation pattern of metabolites in response to cold stress [58]. This result further indicated that environment stress could affect the nutritional constituents of buckwheat, on the other side, buckwheat may modulate the metabolites accumulation to increase tolerance to environmental stress [46,58,72]. In addition, this metabolites result was consistent with the transcriptomic data, in which DEGs in ‘metabolic pathways’ (33.26%), ‘carbohydrate metabolism,’ ‘amino acid metabolism,’ ‘lipid metabolism’ and ‘nucleotide metabolism’ underwent significant expressional changes during salt stress (Figures S6 and S7). Flavonoids are the major secondary metabolites with diverse biological activities in F. tataricum. According to the transcriptomic data, 61 unigenes in flavonoid biosynthesis were up-regulated and 31 unigenes were down-regulated following salt treatment (Table S5). Correspondingly, metabolic analysis showed massive accumulation of rutin following salt stress, the result was corresponded with the most of genes expression level of rutin biosynthesis pathway in buckwheat. Carotenoids biosynthetic metabolites plays a major role in photosynthesis, which is one of the most severely affected pathways in plants response to salt stress [6,73]. Here, 23 genes from the carotenoid biosynthesis pathway were up-regulated and 23 down-regulated, however, the content of the six main carotenoids were not strongly affected by salt stress, which may be the result of the limited treatment course as discussed above. Other factors also could affect the carotenoids synthesis in buckwheat, it has been reported that the light-grown buckwheat sprouts exhibited higher and increased content of carotenoids during the mid-cultivation when compared with dark condition, while, the decreased trend was monitored under the longer culture time, the accumulation pattern of carotenoid was also fluctuated with the seed development stages, different organs and cultivars of tartary buckwheat [56].
In addition, it has been shown that many transcription factors (TFs), including NACs, MYBs and WRKYs, are involved in the regulation of the salt stress response in plants [74,75]. TaMYB56-B in wheat and MdSIMYB1 in apple were reported to be involved in the salt stress response [76] and the over-expression of TaMYB73, AtMYB20 and OsMYB3R-2 in Arabidopsis can enhance salt tolerance [17,18,77]. It has been reported that eight R2R3-MYB genes from buckwheat serve a functional role in abiotic stress tolerance [78]. FtMYB10 from buckwheat has been reported to act as a negative regulator of salt and drought stress, this mechanism has been linked to ABA signaling feedback regulation [79]. Our study has identified a number of transcription factors, including NAC, bHLH and MYB genes, which respond to salt stress in buckwheat seedlings. Of the 2664 identified transcription factors, the bHLHs, ERFs, MYBs, NACs and WRKYs, were significantly differentially regulated during salt stress. Totally, 30 and 25 MYBs and MYBs-related genes, 32 and 4 NACs, 8 and 33 bHLHs, 25 and 6 WRKYs were up or down-regulated to response to salt stress, respectively (Figure 8, Table S7), which may play an important role in regulating salt stress in buckwheat. It has been demonstrated that transcription factors are involved in modulating metabolic pathways in response to abiotic or biotic stressors [20,21,77]. Understanding these processes at the molecular level will play an important role in stabilizing crop performance under high salt conditions [65].

4. Materials and Methods

4.1. Plant Materials, Culture Conditions and Salinity Treatment

F. tataricum seeds (Yu Qiao 2) were germinated and cultured in soil under light/dark (16/8 h,13000 flx) conditions at 25 °C with approximately 50% humidity. 10-days-old buckwheat seedlings were treated with exogenously NaCl solution (100 mM). The control and salt-treated samples were harvested after 1, 2, 4 and 5 days without roots. Samples were rapidly frozen in liquid nitrogen after collection and stored at −80 °C until analysis.

4.2. RNA Isolation, Library Preparation and Transcriptome Sequencing

For the transcriptome sequencing analysis, the 4 days samples of three repeats were precooled and ground into powder in a mortar using liquid nitrogen. Total RNA was extracted from samples using a Trizol reagent (Invitrogen). Finally, the purified double-stranded cDNA samples were further enriched by polymerase chain reaction (PCR) to construct the final cDNA libraries that were sequenced using Hiseq 2000 (Illumina Technologies) to generate 150 bp paired ends.

4.3. Illumina Sequencing and De Novo ASSEMBLY

The raw reads were filtered low quality and removal of adapter sequences, then get clean reads do de novo assembly by Trinity software package [80]. The Trinity program combines overlapping reads into longer contig, through the relationship between reads. Then confirm different contig from same sample and each distance, connect those contig to get unigenes. In order to annotate the assembled transcripts, blast searches were performed against databases including the Nr, Nt, Swiss-Prot, KEGG and COG with the threshold of e-value < 0.00001), along with function annotation of each sequence. The GO annotations were using Blast2GO program and using WEGO software for gene function distributions. The genome information of buckwheat was used as reference information to annotate the transcripts data [81].

4.4. Availability of Supporting Data

The sequencing data in this study were stored in the NCBI Sequence Read Archive (SRA, http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/Traces/sra/) with the accession number of PRJNA528524.

4.5. Total RNA Extraction and cDNA Synthesis

Total RNA was isolated from F. tataricum frozen samples using the RNeasy Plant Mini Kit (Qiagen, Valencia, CA, USA). First-strand cDNA was synthesized from 1 μg of high-quality total RNA with a ReverTra Ace-R reverse transcription (RT) kit using oligo(dT)18 primer (Toyobo Co., Ltd., Osaka, Japan). A 20-fold dilution of 20 μL of the resulting cDNA was used as a template for quantitative real-time PCR.

4.6. Quantitative Real-Time PCR (qRT-PCR) Analysis

The quantitative real-time PCR (qRT-PCR) was performed and the primers were designed and used according to previous reference [82]. The templates were amplified under following conditions—95 °C for 5 min, followed by 35 cycles at 95 °C for 15 s, annealing for 15 s at 56 °C and elongation for 20 s at 72 °C. Real-time PCR reactions were run on a MiniOpticon system (Bio-Rad Laboratories, Hercules, CA, USA) using the SYBR Green Real-Time PCR Master Mix (Toyobo, Osaka, Japan). The housekeeping gene FtActin (KC571237) was used as the internal control to normalize gene expression, as described previously [56]. All reactions were performed in triplicate, addition each run contained a series of standards and a negative control (using water instead of cDNA). The relative expression levels were calculated as differences in cycle threshold (Ct) between genes and actin and analyzed using the 2−△△Ct method [83].

4.7. Extraction and Quantitative HPLC Analysis for Flavonoids

The materials were freeze dried and then ground to a powder before extraction using a method described previously with some modifications [53]. 0.10 g of each powdered sample was extracted with 5 mL of 80% methanol (MeOH), after extraction with sonicator at 25 °C for 1 h, the centrifugation was performed at 12,000 × g for 10 min, the supernatant was filtered and then used for analysis. The HPLC analysis of flavonoids was performed on a Shimadzu class-VP HPLC system (Kyoto, Japan) equipped with a Capcell PAK ODS column (250 mm × 4.6 mm, 5 μm; Shiseido, Tokyo, Japan). The mobile phase consisted of (A) MeOH and (B) water—0.1% formic acid (v/v). The gradient condition was as follows—0 min, 95% A; 4 min, 95–85% A; 9 min, 85% A; 14 min, 85–80% A; 24 min, 80% A; 54 min, 80–70% A; 55 min, 70–55% A; 65 min, 55% A; 75 min, 55–44% A; 77.0 min, 44–40% A; 79 min, 40% A; 80 min, 40–20% A; 90 min, 20% A; 91.0 min, 20–95% A; and 98.0 min, 95% A. The column was maintained at 30 °C, the flow rate was 1.0 mL/min, the injection volume was 20 μL and the detected wavelength was 280 nm. Each compound was quantified according to the peak areas. Results were expressed as milligrams per gram dry weight.

4.8. Liphophilic and Hydrophilic Metabolites Analysis Using GC-TOF-MS

The extraction and analysis of hydrophilic and liphophilic metabolites were performed as described previously [55,70,84,85]. For extraction of hydrophilic compounds, a total of 20 mg of ground sample was extracted with 1 mL of a mixed solvent of methanol/water/chloroform (2.5:1:1, v/v/v). Ribitol (120 μL, 0.2 mg/mL) was added as an internal standard (IS). Extraction was performed at 37 °C using a thermomixer compact (Eppendorf AG, Germany). For methoxime derivatization, 80 μL of methoxyamine hydrochloride (20,000 ppm) in pyridine was added and shaken for 90 min at 30 °C. After the addition of 80 μL of N-methyl-N-trimethylsilyltrifluoroacetamide (MSTFA), the mixtures were incubated for 30 min at 37 °C. The metabolites were analyzed by GC-TOF-MS using an Agilent 7890A gas chromatograph (Agilent, Atlanta, GA, USA) coupled to a Pegasus HT TOF mass spectrometer (LECO, St. Joseph, MI). Derivatized sample (1 μL) was separated on a fused-silica capillary column (30 m × 0.25 mm id) coated with 0.25 mm of CPSIL 8 CB low bleed (Varian Inc., Palo Alto, CA, USA). The parameters of instrument and temperature program was set as mentioned in reference [70]. The quantification of each compounds was based on the peak area ratio relative to that of the IS [55,70]. To extract the lipophilic metabolites, fine-ground samples (0.05 g) were added to 3 mL of ethanol containing 0.1% ascorbic acid (w/v) and 5 µL of 5 α-cholestane (100 µg/mL), which acted as IS and placed in a water bath at 85 °C for 5 min [84]. The mixture was added to 120 µL potassium hydroxide (80%, w/v) for saponification. The samples were then directly treated in ice for 5 min and 1.5 mL of deionized water and hexane were added to each sample and this was mixed for 20 s by vortex and centrifugation (1200× g, 5 min, 4 °C). The supernatant was separated in fresh tubes. The pellet was subjected to re-extraction using 1.5 mL hexane. The hexane fraction was evaporated under a stream of N2 to a volume of 200 µL and dried in a centrifugal concentrator (CC-105). For derivatization, 30 µL of MSTFA and 30 µL of pyridine were added to each sample and the mixture was incubated (60 °C, 30 min, 1200 rpm) using a Thermomixer Comfort. The samples (1.0 µL each) were injected in split mode (10:1 ratio) and gas chromatography was performed using an Rtx-5MS column (30 m length, 0.25 mm diameter and 0.25 µm thickness). The GCMS-QP2010 Ultra system, with autosampler AOC-20i (Shimadzu), was used for the separation of the lipophilic compounds. The parameters of the instrument and temperature program were set as mentioned in References [84,85]. Calibration samples were prepared by mixing appropriate aliquots of each stock solution. Quantification was performed by means of three-point calibration curves, for which the concentrations of a mixture of authentic standards ranged from 0.25 to 5.0 µg, whereas the amount of the IS was constant at 0.5 µg [84,85].

4.9. Carotenoid Extraction and Analysis

Carotenoids were extracted and measured using HPLC as described previously [56]. Briefly, the carotenoids were extracted from cauliflower samples (0.1 g) by adding 3 mL of ethanol containing 0.1% ascorbic acid (w/v), vortex mixing for 20 s and placing in a water bath at 85 °C for 5 min. The potassium hydroxide (120 μL, 80% w/v) was used for carotenoid saponification (85 °C for 10 min). The cold deionized water (1.5 mL) was added after saponification. β-Apo-8′-carotenal (0.2 mL, 25 μg/mL) was added as an internal standard (IS). The carotenoids were extracted twice with 1.5 mL hexane and followed by centrifugation at 1200g. Aliquots of the extracts were dried under a stream of nitrogen and redissolved in 50:50 (v/v) dichloromethane/methanol before HPLC analysis. The samples were analyzed by an Agilent 1100 HPLC instrument (Massy, France) equipped with a photodiode array detector. The C30 YMC column (250 × 4.6 mm, 3 μm; YMC Co., Kyoto, Japan) was used. The wavelength of monitor was 450 nm. Solvent A consisted of methanol/water (92:8 v/v) with 10mM ammonium acetate; solvent B consisted of 100% methyl tert-butyl ether. The elution rate was 1 mL/min with a gradient program to analyze samples. Gradient elution was performed as follows—0 min, 90% A–10% B; 20 min, 83% A–17% B; 29 min, 75% A–25% B; 35 min, 30% A–0% B; 40 min, 30% A–70% B; 42 min, 25% A–75% B; 45 min, 90% A–10% B; and 55 min, 90% A–10% B. The quantification of individual carotenoid was calculated with carotenoid standards (Carote Nature, Lupsingen, Switzerland) according to the peak area ratios with internal standard (IS).

4.10. Data Processing and Statistical Analysis

All analyses were performed at least three times. Experimental data were processed by analysis of variance (ANOVA) and significant differences among the means were determined by Duncan’s multiple-range test (SPSS 22.0, Chicago: SPSS Inc. IL, USA). Pearson correlation analysis was calculated by SPSS 22.0, then cluster analysis and heatmap were drawn by R (http://www.r-project.org/).

5. Conclusions

To our knowledge, no previous studies have provided a comprehensive description of salt-resistance in buckwheat using transcriptome and metabolite profiling analysis and the interactive effects of salinity on both transcriptional and metabolic regulation of buckwheat was unclear. In this study, we performed a comprehensive analysis of the transcriptome and metabolome of salt treated-buckwheat to understand the effects of salinity on buckwheat. Many genes and metabolites involved in primary and secondary metabolic biosynthesis were changed significantly following treatment. Those effects especially included amino acids and organic acids biosynthesis, phenylpropanoid biosynthesis and flavonoid biosynthesis in buckwheat. In all, this research illuminated more information to uncover the regulation mechanism of salt tolerance in buckwheat at the transcriptional and metabolic level. Our study could be helpful for the better understanding of the salt tolerance mechanisms in this important crop.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2218-1989/9/10/225/s1, Figure S1: Length distribution of all unigenes from F. tataricum transcriptome. Figure S2: Annotation of all unigenes from F. tataricum transcriptome. Number of genes were annotated by five representatively databases (Nr, KEGG, COG, Swissprot and Nt). Figure S3: GO analysis of all unigenes based on biological process, cellular component and molecular function categories. Figure S4: COG function classification of all unigenes sequences from F. tataricum transcriptome. Figure S5: Distribution of differentially expressed genes (DEGs) between control and salt treated buckwheat. DEGs was selected by p-value < 0.001 and |log2FoldChange| > 2. Figure S6: DEGs involved in carbohydrate metabolism between control and salt treated buckwheat. Figure S7: DEGs involved in amino acid metabolism between control and salt treated buckwheat. Table S1: Statistics of annotations for assembled unigenes of F. tataricum in different public databases. Table S2: DEGs genes with KEGG pathway annotation. Table S3: The relative comparation of hydrophilic metabolites in control (Con.) and salt-treated (Tr) buckwheat samples. Table S4: The content of liphophilic metabolites in control (Con) and salt-treated (Tr) buckwheat samples (μg/g dry weight). Table S5: Different expression of phenylpropanoid biosynthetic genes under salt stress. Table S6: The differentially expressed genes of carotenoids biosynthetic pathway in buckwheat under salt stress. Table S7: Identification of differentially expressed transcription factors (TFs) in buckwheat to response salt stress.

Author Contributions

X.L. designed the experiments, X.L., W.M., J.K.K. and H.J.K. conducted the experiments, X.L. and W.M. and F.Y. analyzed the data and wrote the manuscript. C.J., X.H. and W.A. revised and edited the manuscript.

Funding

This work was funded by the National Natural Science Foundation of China No. 31600234 (X.L.), the Fundamental Research Funds for the Central Universities Program No. 2662015QD005 (X.L.).

Acknowledgments

We thank Xiao Jun Han from Agricultural Bureau of Jingyuan province, Shanxi, China, for providing materials.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

GOGene ontology
KEGGKyoto Encyclopedia of Genes and Genomes
KOGEukaryotic Ortholog Groups
NrNCBI non-redundant protein sequences
NtNCBI non-redundant nucleotide sequences
DEGsDifferential expressed genes
HPLCHigh-performance liquid chromatography
GC-TOF-MSGas Chromatography Time-Of-Flight Mass Spectrometry

References

  1. Tester, M.; Davenport, R. Na+ tolerance and Na+ transport in higher plants. Ann. Bot. 2003, 91, 503–527. [Google Scholar] [CrossRef] [PubMed]
  2. Yang, Y.; Guo, Y. Elucidating the molecular mechanisms mediating plant salt-stress responses. New Phytol. 2018, 217, 523–539. [Google Scholar] [CrossRef] [PubMed]
  3. Xue, G.P. Characterisation of the DNA-binding profile of barley HvCBF1 using an enzymatic method for rapid, quantitative and high-throughput analysis of the DNA-binding activity. Nucleic Acids Res. 2002, 30, e77. [Google Scholar] [CrossRef] [PubMed]
  4. Shokri-Gharelo, R.; Noparvar, P.M. Molecular response of canola to salt stress: Insights on tolerance mechanisms. PeerJ 2018, 6, e4822. [Google Scholar] [CrossRef] [PubMed]
  5. Golldack, D.; Li, C.; Mohan, H.; Probst, N. Tolerance to drought and salt stress in plants: Unraveling the signaling networks. Front. Plant Sci. 2014, 5, 151. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Stoeva, N.; Kaymakanova, M. Effect of salt stress on the growth and photosynthesis rate of bean plants (Phaseolus vulgaris L.). J. Cent. Eur. Agric. 2008, 9, 385–391. [Google Scholar]
  7. Zhang, H.; Han, B.; Wang, T.; Chen, S.; Li, H.; Zhang, Y.; Dai, S. Mechanisms of plant salt response: Insights from proteomics. J. Proteome Res. 2012, 11, 49–67. [Google Scholar] [CrossRef] [PubMed]
  8. Payyavula, R.S.; Navarre, D.A.; Kuhl, J.C.; Pantoja, A.; Pillai, S.S. Differential effects of environment on potato phenylpropanoid and carotenoid expression. BMC Plant Biol. 2012, 12, 39. [Google Scholar] [CrossRef]
  9. Petridis, A.; Therios, I.; Samouris, G.; Tananaki, C. Salinity-induced changes in phenolic compounds in leaves and roots of four olive cultivars (Olea europaea L.) and their relationship to antioxidant activity. Environ. Exp. Bot. 2012, 79, 37–43. [Google Scholar] [CrossRef]
  10. Sreenivasulu, N.; Grimm, B.; Wobus, U.; Weschke, W. Differential response of antioxidant compounds to salinity stress in salt-tolerant and salt-sensitive seedlings of foxtail millet (Setaria italica). Physiol. Plant 2000, 109, 435–442. [Google Scholar] [CrossRef]
  11. Borghesi, E.; Gonzà lez-Miret, M.L.; Escudero-Gilete, M.L.; Malorgio, F.; Heredia, F.J.; Meléndez-Martínez, A.J. Effects of salinity stress on carotenoids, anthocyanins and color of diverse tomato genotypes. J. Agric. Food Chem. 2011, 59, 11676–11682. [Google Scholar] [CrossRef] [PubMed]
  12. Abdallah, S.B.; Aung, B.; Amyot, L.; Lalin, I.; Lachâal, M.; Karray-Bouraoui, N.; Hannoufa, A. Salt stress (NaCl) affects plant growth and branch pathways of carotenoid and flavonoid biosyntheses in Solanum nigrum. Acta Physiol. Plant. 2016, 38, 1–13. [Google Scholar] [CrossRef]
  13. Tuteja, N. Mechanisms of high salinity tolerance in plants. Methods Enzymol. 2007, 428, 419–438. [Google Scholar] [PubMed]
  14. Agati, G.; Azzarello, E.; Pollastri, S.; Tattini, M. Flavonoids as antioxidants in plants: Location and functional significance. Plant Sci. 2012, 196, 67–76. [Google Scholar] [CrossRef] [PubMed]
  15. Zhao, Y.; Tian, X.; Wang, F.; Zhang, L.; Xin, M.; Hu, Z.; Yao, Y.; Ni, Z.; Sun, Q.; Peng, H. Characterization of wheat MYB genes responsive to high temperatures. BMC Plant Biol. 2017, 17, 208. [Google Scholar] [CrossRef] [PubMed]
  16. Ma, J.; Lu, J.; Xu, J.; Duan, B.; He, X.; Liu, J. Genome-wide identification of WRKY genes in the desert Poplar Populus euphratica and adaptive evolution of the genes in response to salt stress. Evol. Bioinform. 2015, 11, 47–55. [Google Scholar]
  17. He, Y.; Li, W.; Lv, J.; Jia, Y.; Wang, M.; Xia, G. Ectopic expression of a wheat MYB transcription factor gene, TaMYB73, improves salinity stress tolerance in Arabidopsis thaliana. J. Exp. Bot. 2012, 63, 1511–1522. [Google Scholar] [CrossRef] [PubMed]
  18. Cui, M.H.; Yoo, K.S.; Hyoung, S.; Nguyen, H.T.; Kim, Y.Y.; Kim, H.J.; Ok, S.H.; Yoo, S.D.; Shin, J.S. An Arabidopsis R2R3-MYB transcription factor, AtMYB20, negatively regulates type 2C serine/threonine protein phosphatases to enhance salt tolerance. FEBS Lett. 2013, 587, 1773–1778. [Google Scholar] [CrossRef]
  19. Zhong, H.; Guo, Q.Q.; Chen, L.; Ren, F.; Wang, Q.Q.; Zheng, Y.; Li, X.B. Two Brassica napus genes encoding NAC transcription factors are involved in response to high-salinity stress. Plant Cell Rep. 2012, 31, 1991–2003. [Google Scholar] [CrossRef]
  20. Nakashima, K.; Tran, L.S.; Van Nguyen, D.; Fujita, M.; Maruyama, K.; Todaka, D.; Ito, Y.; Hayashi, N.; Shinozaki, K.; Yamaguchi-Shinozaki, K. Functional analysis of a NAC-type transcription factor OsNAC6 involved in abiotic and biotic stress-responsive gene expression in rice. Plant J. 2007, 51, 617–630. [Google Scholar] [CrossRef]
  21. Xu, Z.; Zhaxi, G.; Wang, C.; Xue, F.; Zhang, H.; Ji, W. Wheat NAC transcription factor TaNAC29 is involved in response to salt stress. Plant Physiol. Biochem. 2015, 96, 356–363. [Google Scholar] [CrossRef] [PubMed]
  22. Zhu, M.; Chen, G.; Zhang, J.; Zhang, Y.; Xie, Q.; Zhao, Z.; Pan, Y.; Hu, Z. The abiotic stress-responsive NAC-type transcription factor SlNAC4 regulates salt and drought tolerance and stress-related genes in tomato (Solanum lycopersicum). Plant Cell Rep. 2014, 33, 1851–1863. [Google Scholar] [CrossRef] [PubMed]
  23. Cheng, P.; Gao, J.; Feng, Y.; Zhang, Z.; Liu, Y.; Fang, W.; Chen, S.; Chen, F.; Jiang, J. The chrysanthemum leaf and root transcript profiling in response to salinity stress. Gene 2018, 674, 161–169. [Google Scholar] [CrossRef] [PubMed]
  24. Golldack, D.; Luking, I.; Yang, O. Plant tolerance to drought and salinity: Stress regulating transcription factors and their functional significance in the cellular transcriptional network. Plant Cell Rep. 2011, 30, 1383–1391. [Google Scholar] [CrossRef] [PubMed]
  25. Park, N.I.; Li, X.; Suzuki, T.; Kim, S.J.; Woo, S.H.; Park, C.H.; Park, S.U. Differential expression of anthocyanin biosynthetic genes and anthocyanin accumulation in tartary buckwheat cultivars ‘Hokkai t8’ and ‘Hokkai t10’. J. Agric. Food Chem. 2011, 59, 2356–2361. [Google Scholar] [CrossRef]
  26. Mendler-Drienyovszki, N.; Cal, A.J.; Dobranszki, J. Progress and prospects for interspecific hybridization in buckwheat and the genus Fagopyrum. Biotechnol. Adv. 2013, 31, 1768–1775. [Google Scholar] [CrossRef] [PubMed]
  27. Kim, S.J.; Maeda, T.; Sarker, M.Z.; Takigawa, S.; Matsuura-Endo, C.; Yamauchi, H.; Mukasa, Y.; Saito, K.; Hashimoto, N.; Noda, T.; et al. Identification of anthocyanins in the sprouts of buckwheat. J. Agric. Food Chem. 2007, 55, 6314–6318. [Google Scholar] [CrossRef]
  28. Zheng, C.; Hu, C.; Ma, X.; Peng, C.; Zhang, H.; Qin, L. Cytotoxic phenylpropanoid glycosides from Fagopyrum tataricum (L.) Gaertn. Food Chem. 2012, 132, 433–438. [Google Scholar] [CrossRef]
  29. Matsui, T.; Kudo, A.; Tokuda, S.; Matsumoto, K.; Hosoyama, H. Identification of a new natural vasorelaxatant compound, (+)-osbeckic acid, from rutin-free tartary buckwheat extract. J. Agric. Food Chem. 2010, 58, 10876–10879. [Google Scholar] [CrossRef]
  30. Yang, N.; Ren, G. Determination of D-chiro-Inositol in tartary buckwheat using high-performance liquid chromatography with an evaporative light-scattering detector. J. Agric. Food Chem. 2008, 56, 757–760. [Google Scholar] [CrossRef]
  31. Ren, Q.; Wu, C.; Ren, Y.; Zhang, J. Characterization and identification of the chemical constituents from tartary buckwheat (Fagopyrum tataricum Gaertn) by high performance liquid chromatography/photodiode array detector/linear ion trap FTICR hybrid mass spectrometry. Food Chem. 2013, 136, 1377–1389. [Google Scholar] [CrossRef] [PubMed]
  32. Wijngaard, H.H.; Arendt, E.K. Buckwheat. Cereal Chem. 2006, 83, 391–401. [Google Scholar] [CrossRef]
  33. Middleton, E., Jr.; Kandaswami, C.; Theoharides, T.C. The effects of plant flavonoids on mammalian cells: Implications for inflammation, heart disease and cancer. Pharmacol. Rev. 2000, 52, 673–751. [Google Scholar] [PubMed]
  34. Panwar, A.; Gupta, N.; Chauhan, R.S. Biosynthesis and accumulation of flavonoids in Fagopyrum spp. Eur. J. Plant Sci. Biotechnol. 2012, 6, 17–26. [Google Scholar]
  35. Suzuki, T.; Morishita, T.; Mukasa, Y.; Takigawa, S.; Yokota, S.; Ishiguro, K.; Noda, T. Breeding of ‘Manten-Kirari,’ a non-bitter and trace-rutinosidase variety of Tartary buckwheat (Fagopyrum tataricum Gaertn.). Breed. Sci. 2014, 64, 344–350. [Google Scholar] [CrossRef] [PubMed]
  36. Nishimura, M.; Ohkawara, T.; Sato, Y.; Satoh, H.; Suzuki, T.; Ishiguro, K.; Noda, T.; Morishita, T.; Nishihira, J. Effectiveness of rutin-rich Tartary buckwheat (Fagopyrum tataricum Gaertn.) ‘Manten-Kirari’ in body weight reduction related to its antioxidant properties: A randomised, double-blind, placebo-controlled study. J. Funct. Foods 2016, 26, 460–469. [Google Scholar] [CrossRef]
  37. Watanabe, M. Catechins as antioxidants from buckwheat (Fagopyrum esculentum Moench) groats. J. Agric. Food Chem. 1998, 46, 839–845. [Google Scholar] [CrossRef]
  38. Holasova, M.; Fiedlerova, V.H.; Orsak, M.; Lachman, J. Vavreinova Buckwheat—The source of antioxidant activity in functional foods. Food Res. Int. 2002, 35, 207–211. [Google Scholar] [CrossRef]
  39. Krkošková, B.; Mrázová, Z. Prophylactic components of buckwheat. Food Res. Int. 2005, 38, 561–568. [Google Scholar] [CrossRef]
  40. Park, J.W.; Kang, D.B.; Kim, C.W.; Koh, S.H.; Yum, H.Y.; Kim, K.E.; Hong, C.S.; Lee, K.Y. Identification and characterization of the major allergens of buckwheat. Allergy 2015, 55, 1035–1041. [Google Scholar] [CrossRef]
  41. Karki, R.; Park, C.H.; Kim, D.W. Extract of buckwheat sprouts scavenges oxidation and inhibits pro-inflammatory mediators in lipopolysaccharide-stimulated macrophages (RAW264.7). J. Integr. Med. 2013, 11, 246–252. [Google Scholar] [CrossRef] [PubMed]
  42. Jin, H.M.; Wei, P. Anti-fatigue properties of tartary buckwheat extracts in mice. Int. J. Mol. Sci. 2011, 12, 4770–4780. [Google Scholar] [CrossRef] [PubMed]
  43. Ren, W.; Qiao, Z.; Wang, H.; Zhu, L.; Zhang, L.; Lu, Y.; Cui, Y.; Zhang, Z.; Wang, Z. Tartary buckwheat flavonoid activates caspase 3 and induces HL-60 cell apoptosis. Methods Find. Exp. Clin. Pharmacol. 2001, 23, 427–432. [Google Scholar] [CrossRef] [PubMed]
  44. Ren, W.; Qiao, Z.; Wang, H.; Zhu, L.; Zhang, L.; Lu, Y.; Zhang, Z.; Wang, Z. Molecular basis of fas and cytochrome c pathways of apoptosis induced by tartary buckwheat flavonoid in HL-60 cells. Methods Find. Exp. Clin. Pharmacol. 2003, 25, 431–436. [Google Scholar] [CrossRef] [PubMed]
  45. Guo, X.; Zhu, K.; Hui, Z.; Yao, H. Anti-tumor activity of a novel protein obtained from tartary buckwheat. Int. J. Mol. Sci. 2010, 11, 5201–5211. [Google Scholar] [CrossRef] [PubMed]
  46. Kim, H.-J.; Park, K.-J.; Lim, J.-H. Metabolomic analysis of phenolic compounds in buckwheat (Fagopyrum esculentum M.) sprouts treated with methyl Jasmonate. J. Agric. Food Chem. 2011, 59, 5707–5713. [Google Scholar] [CrossRef] [PubMed]
  47. Yang, H.J.; Lim, J.H.; Park, K.J.; Kang, S.; Kim, D.S.; Park, S. Methyl jasmolate treated buckwheat sprout powder enhances glucose metabolism by potentiating hepatic insulin signaling in estrogen-deficient rats. Nutrition 2016, 32, 129–137. [Google Scholar] [CrossRef] [PubMed]
  48. Seo, J.-M.; Arasu, M.V.; Kim, Y.-B.; Park, S.U.; Kim, S.-J. Phenylalanine and LED lights enhance phenolic compound production in Tartary buckwheat sprouts. Food Chem. 2015, 177, 204–213. [Google Scholar] [CrossRef] [PubMed]
  49. Chen, L.H.; Zhang, B.; Xu, Z.Q. Salt tolerance conferred by overexpression of Arabidopsis vacuolar Na(+)/H (+) antiporter gene AtNHX1 in common buckwheat (Fagopyrum esculentum). Transgenic Res. 2008, 17, 121–132. [Google Scholar] [CrossRef]
  50. Qin, P.; Wang, Q.; Shan, F.; Hou, Z.; Ren, G. Nutritional composition and flavonoids content of flour from different buckwheat cultivars. Int. J. Food Sci. Technol. 2010, 45, 951–958. [Google Scholar] [CrossRef]
  51. Wan, Y.; Xiang, D.B.; Zeng, X.L.; Wang, X.; Yuan, Z.H.; Tan, M.L.; Zou, L.; Zhao, G. Salt stress influence on yield of tartary buckwheat sprouts and flavonoids content. Sci. Technol. Food Ind. 2016, 37, 328–332. [Google Scholar]
  52. Zhu, Y.; Guo, Y. Optimization of culture conditions for accumulating gamma-aminobutyric acid (GABA) in germinated tartary buckwheat under salt stress by response surface methodology. Food Sci. 2015, 36, 96–100. [Google Scholar]
  53. Li, X.; Kim, Y.B.; Kim, Y.; Zhao, S.; Kim, H.H.; Chung, E.; Lee, J.H.; Park, S.U. Differential stress-response expression of two flavonol synthase genes and accumulation of flavonols in tartary buckwheat. J. Plant Physiol. 2013, 170, 1630–1636. [Google Scholar] [CrossRef] [PubMed]
  54. Lu, Q.H.; Wang, Y.Q.; Song, J.N.; Yang, H.B. Transcriptomic identification of salt-related genes and de novo assembly in common buckwheat (F. esculentum). Plant Physiol. Biochem. 2018, 127, 299–309. [Google Scholar] [CrossRef] [PubMed]
  55. Park, S.Y.; Lim, S.H.; Ha, S.H.; Yeo, Y.; Park, W.T.; Kwon, D.Y.; Park, S.U.; Kim, J.K. Metabolite profiling approach reveals the interface of primary and secondary metabolism in colored cauliflowers (Brassica oleracea, L. ssp. botrytis). J. Agric. Food Chem. 2013, 61, 6999–7007. [Google Scholar] [CrossRef] [PubMed]
  56. Tuan, P.A.; Thwe, A.A.; Kim, J.K.; Kim, Y.B.; Lee, S.; Park, S.U. Molecular characterisation and the light-dark regulation of carotenoid biosynthesis in sprouts of tartary buckwheat (Fagopyrum tataricum Gaertn.). Food Chem. 2013, 141, 3803–3812. [Google Scholar] [CrossRef] [PubMed]
  57. Suzuki, N.; Rivero, R.M.; Shulaev, V.; Blumwald, E.; Mittler, R. Abiotic and biotic stress combinations. New Phytol. 2014, 203, 32–43. [Google Scholar] [CrossRef]
  58. Jeon, J.; Kim, J.K.; Wu, Q.; Park, S.U. Effects of cold stress on transcripts and metabolites in tartary buckwheat (Fagopyrum tataricum). Environ. Exp. Bot. 2018, 155, 488–496. [Google Scholar] [CrossRef]
  59. Li, S.-J.; Bai, Y.-C.; Li, C.-L.; Yao, H.-P.; Chen, H.; Zhao, H.-X.; Wu, Q.J.A.P.P. Anthocyanins accumulate in tartary buckwheat (Fagopyrum tataricum) sprout in response to cold stress. Acta Physiol. Plant. 2015, 37, 159. [Google Scholar] [CrossRef]
  60. Guo, R.; Yang, Z.; Li, F.; Yan, C.; Zhong, X.; Liu, Q.; Xia, X.; Li, H.; Zhao, L. Comparative metabolic responses and adaptive strategies of wheat (Triticum aestivum) to salt and alkali stress. BMC Plant Biol. 2015, 15, 170. [Google Scholar] [CrossRef] [PubMed]
  61. Hernández, J.A.; Jiménez, A.; Mullineaux, P.; Sevilia, F. Tolerance of pea (Pisum sativum L.) to long-term salt stress is associated with induction of antioxidant defences. Plant Cell Environ. 2000, 23, 853–862. [Google Scholar] [CrossRef]
  62. Pascale, S.D.; Maggio, A.; Angelino, G.; Graziani, G. Effect of salt stress on water relations and antioxidant activity in tomato. Acta Hortic. 2003, 613, 39–46. [Google Scholar] [CrossRef]
  63. Huang, C.; He, W.; Guo, J.; Chang, X.; Su, P.; Zhang, L. Increased sensitivity to salt stress in an ascorbate-deficient Arabidopsis mutant. J. Exp. Bot. 2005, 56, 3041–3049. [Google Scholar] [CrossRef] [PubMed]
  64. Zhang, H.; He, L.; Cai, L. Transcriptome Sequencing: RNA-Seq. Methods Mol. Biol. 2018, 1754, 15–27. [Google Scholar]
  65. Upadhyay, A.; Gaonkar, T.; Upadhyay, A.K.; Jogaiah, S.; Shinde, M.P.; Kadoo, N.Y.; Gupta, V.S. Global transcriptome analysis of grapevine (Vitis vinifera L.) leaves under salt stress reveals differential response at early and late stages of stress in table grape cv. Thompson Seedless. Plant Physiol. Biochem. 2018, 129, 168–179. [Google Scholar] [CrossRef] [PubMed]
  66. Wu, Q.; Bai, X.; Zhao, W.; Xiang, D.; Wan, Y.; Yan, J.; Zou, L.; Zhao, G. De novo assembly and analysis of tartary buckwheat (Fagopyrum tataricum Garetn.) transcriptome discloses key regulators involved in salt-stress response. Genes 2017, 8, 255. [Google Scholar] [CrossRef] [PubMed]
  67. Linic, I.; Samec, D.; Gruz, J.; Vujcic Bok, V.; Strnad, M.; Salopek-Sondi, B. Involvement of phenolic acids in short-term adaptation to salinity stress is species-specific among Brassicaceae. Plants 2019, 8, 155. [Google Scholar] [CrossRef]
  68. Enfissi, E.M.A.; Barneche, F.; Ahmed, I.; Lichtlé, C.; Gerrish, C.; McQuinn, R.P.; Giovannoni, J.J.; Lopez-Juez, E.; Bowler, C.; Bramley, P.M.; et al. Integrative transcript and metabolite analysis of nutritionally enhanced DE-ETIOLATED1 downregulated tomato fruit. Plant Cell 2010, 22, 1190–1215. [Google Scholar] [CrossRef] [PubMed]
  69. Ayenew, B.; Degu, A.; Manela, N.; Perl, A.; Shamir, M.O.; Fait, A. Metabolite profiling and transcript analysis reveal specificities in the response of a berry derived cell culture to abiotic stresses. Front. Plant Sci. 2015, 6, 728. [Google Scholar] [CrossRef] [Green Version]
  70. Park, C.H.; Morgan, A.M.A.; Park, B.B.; Lee, S.Y.; Lee, S.; Kim, J.K.; Park, S.U. Metabolic analysis of four cultivars of Liriope platyphylla. Metabolites 2019, 9, 59. [Google Scholar] [CrossRef]
  71. Park, C.H.; Yeo, H.J.; Park, Y.E.; Baek, S.A.; Kim, J.K.; Park, S.U. Transcriptome analysis and metabolic profiling of Lycoris Radiata. Biology 2019, 8, 63. [Google Scholar] [CrossRef] [PubMed]
  72. Lim, J.H.; Park, K.J.; Kim, B.K.; Jeong, J.W.; Kim, H.J. Effect of salinity stress on phenolic compounds and carotenoids in buckwheat (Fagopyrum esculentum M.) sprout. Food Chem. 2012, 135, 1065–1070. [Google Scholar] [CrossRef] [PubMed]
  73. Sudhir, P.; Murthy, S.D.S. Effects of salt stress on basic processes of photosynthesis. Photosynthetica 2004, 42, 481–486. [Google Scholar] [CrossRef]
  74. Eulgem, T.; Somssich, I.E. Networks of WRKY transcription factors in defense signaling. Curr. Opin. Plant Biol. 2007, 10, 366–371. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Jin, H.; Dong, D.; Yang, Q.; Zhu, D. Salt-responsive transcriptome profiling of suaeda glaucavia RNA sequencing. PLoS ONE 2016, 11, e0150504. [Google Scholar] [CrossRef] [PubMed]
  76. Shen, X.; Guo, X.; Guo, X.; Zhao, D.; Zhao, W.; Chen, J.; Li, T. PacMYBA, a sweet cherry R2R3-MYB transcription factor, is a positive regulator of salt stress tolerance and pathogen resistance. Plant Physiol. Biochem. 2017, 112, 302–311. [Google Scholar] [CrossRef] [PubMed]
  77. Dai, X.; Xu, Y.; Ma, Q.; Xu, W.; Wang, T.; Xue, Y.; Chong, K. Overexpression of an R1R2R3 MYB Gene, OsMYB3R-2, Increases tolerance to freezing, drought and salt stress in transgenic arabidopsis. Plant Physiol. 2007, 143, 1739–1751. [Google Scholar] [CrossRef]
  78. Gao, F.; Zhao, H.X.; Yao, H.P.; Li, C.L.; Chen, H.; Wang, A.H.; Park, S.U.; Wu, Q. Identification, isolation and expression analysis of eight stress-related R2R3-MYB genes in tartary buckwheat (Fagopyrum tataricum). Plant Cell Rep. 2016, 35, 1385–1396. [Google Scholar] [CrossRef]
  79. Gao, F.; Yao, H.; Zhao, H.; Zhou, J.; Luo, X.; Huang, Y.; Li, C.; Chen, H.; Wu, Q. Tartary buckwheat FtMYB10 encodes an R2R3-MYB transcription factor that acts as a novel negative regulator of salt and drought response in transgenic Arabidopsis. Plant Physiol. Biochem. 2016, 109, 387–396. [Google Scholar] [CrossRef]
  80. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q.; et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 2011, 29, 644–652. [Google Scholar] [CrossRef] [Green Version]
  81. Zhang, L.; Li, X.; Ma, B.; Gao, Q.; Du, H.; Han, Y.; Li, Y.; Cao, Y.; Qi, M.; Zhu, Y.; et al. The tartary buckwheat genome provides insights into rutin biosynthesis and abiotic stress tolerance. Mol. Plant 2017, 10, 1224–1237. [Google Scholar] [CrossRef] [PubMed]
  82. Li, X.; Kim, J.K.; Park, S.Y.; Zhao, S.; Kim, Y.B.; Lee, S.; Park, S.U. Comparative analysis of flavonoids and polar metabolite profiling of Tanno-original and Tanno-high rutin buckwheat. J. Agric. Food Chem. 2014, 62, 2701–2708. [Google Scholar] [CrossRef] [PubMed]
  83. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  84. Kim, T.J.; Choi, J.; Kim, K.W.; Ahn, S.K.; Ha, S.-H.; Choi, Y.; Park, N.I.; Kim, J.K. Metabolite profiling of peppers of various colors reveals relationships between tocopherol, carotenoid and phytosterol content. J. Food Sci. 2017, 82, 2885–2893. [Google Scholar] [CrossRef] [PubMed]
  85. Baek, S.A.; Jung, Y.H.; Lim, S.H.; Park, S.U.; Kim, J.K. Metabolic profiling in chinese cabbage (Brassica rapa, L. subsp. pekinensis) cultivars reveals that glucosinolate content is correlated with carotenoid content. J. Agric. Food Chem. 2016, 64, 4426–4434. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Species statistics and distribution of unigenes in F. tataricum transcriptome. (A) The E-value frequency distribution of Nr annotated unigenes using E-value frequency distribution analysis; (B) Sequence similarity distribution of the Blast hits for each unigenes; (C) The species homologous distribution analysis. It is shown as the percentage of the total homologous sequences in the Genebank non redundant (Nr) database.
Figure 1. Species statistics and distribution of unigenes in F. tataricum transcriptome. (A) The E-value frequency distribution of Nr annotated unigenes using E-value frequency distribution analysis; (B) Sequence similarity distribution of the Blast hits for each unigenes; (C) The species homologous distribution analysis. It is shown as the percentage of the total homologous sequences in the Genebank non redundant (Nr) database.
Metabolites 09 00225 g001
Figure 2. GO analysis of differentially expressed genes (DEGs) based on biological process, cellular component and molecular function categories.
Figure 2. GO analysis of differentially expressed genes (DEGs) based on biological process, cellular component and molecular function categories.
Metabolites 09 00225 g002
Figure 3. The differentially expressed genes (DEGs) with Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment (top 20) in tartary buckwheat in response to salt stress.
Figure 3. The differentially expressed genes (DEGs) with Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment (top 20) in tartary buckwheat in response to salt stress.
Metabolites 09 00225 g003
Figure 4. DEGs involved in secondary metabolites biosynthesis in buckwheat in response to salt stress.
Figure 4. DEGs involved in secondary metabolites biosynthesis in buckwheat in response to salt stress.
Metabolites 09 00225 g004
Figure 5. Overview of the primary metabolic changes of buckwheat in response to salt stress. The fold changes of each metabolite are given by salt-treated relative to control group in 1, 2, 4 and 5 day (cells from left to right, respectively). Different colors (from red, orange, yellow, grey yellow to blue) represent the fold change of the metabolites as indicated in the legend.
Figure 5. Overview of the primary metabolic changes of buckwheat in response to salt stress. The fold changes of each metabolite are given by salt-treated relative to control group in 1, 2, 4 and 5 day (cells from left to right, respectively). Different colors (from red, orange, yellow, grey yellow to blue) represent the fold change of the metabolites as indicated in the legend.
Metabolites 09 00225 g005
Figure 6. Correlation matrix and cluster analysis of results obtained from data on 53 primary metabolites of buckwheat. Each square indicates the Pearson’s correlation coefficient of a pair of compounds and the value for the correlation coefficient is represented by the intensity of the blue or red color, as indicated on the color scale.
Figure 6. Correlation matrix and cluster analysis of results obtained from data on 53 primary metabolites of buckwheat. Each square indicates the Pearson’s correlation coefficient of a pair of compounds and the value for the correlation coefficient is represented by the intensity of the blue or red color, as indicated on the color scale.
Metabolites 09 00225 g006
Figure 7. Expression of the key genes in phenylpropanoid biosynthetic pathway of tartary buckwheat in response to salt stress. PAL, phenylalanine ammonium lyase; CHI, chalcone isomerase; CHS, chalcone synthase; 4CL, 4-coumarate-CoA ligase; C4H, cinnamic acid 4-hydroxylase; F3H, flavanone-3-hydroxylase; F3′H, flavonoid-3′-hydroxylase; FLS, flavonol synthase; DFR, dihydroflavonol reductase; ANS, anthocyanin synthase. The asterisk indicates a significant difference between control (Con.) and salt-treated (Tr.) samples (*, p < 0.05; **, p < 0.01).
Figure 7. Expression of the key genes in phenylpropanoid biosynthetic pathway of tartary buckwheat in response to salt stress. PAL, phenylalanine ammonium lyase; CHI, chalcone isomerase; CHS, chalcone synthase; 4CL, 4-coumarate-CoA ligase; C4H, cinnamic acid 4-hydroxylase; F3H, flavanone-3-hydroxylase; F3′H, flavonoid-3′-hydroxylase; FLS, flavonol synthase; DFR, dihydroflavonol reductase; ANS, anthocyanin synthase. The asterisk indicates a significant difference between control (Con.) and salt-treated (Tr.) samples (*, p < 0.05; **, p < 0.01).
Metabolites 09 00225 g007
Figure 8. The differentially expressed transcription factors (TFs) of buckwheat in response to salt stress. The up and down regulated TFs were represented by the green or blue color, respectively.
Figure 8. The differentially expressed transcription factors (TFs) of buckwheat in response to salt stress. The up and down regulated TFs were represented by the green or blue color, respectively.
Metabolites 09 00225 g008
Table 1. Statistic of sequencing and de novo assembling of transcriptome in F. tataricum.
Table 1. Statistic of sequencing and de novo assembling of transcriptome in F. tataricum.
SampleTotal NumberTotal Length (nt)Mean Length (nt)N50Total Consensus SequencesDistinct ClustersDistinct Singletons
ContigsC197,01748,783,4745031133---
C2103,96151,398,9104941104---
C398,58849,113,0214981130---
T1117,18252,926,2624521007---
T2104,42949,990,9584791081---
T3118,86853,985,9594541001---
UnigenesC158,78960,570,9981030168058,78921,10537,684
C262,55964,289,6881028168662,55922,44540,114
C358,64460,081,2021025166858,64420,63738,007
T164,47062,513,110970163964,47021,29143,179
T260,14660,524,7421006167560,14620,49839,648
T365,74963,160,040961163165,74920,94844,801
All94,950107,581,1091133190094,95038,24256,708
C1-3 and T1-3 mean control (C) and salt-treated (T) samples with three repeats, respectively.
Table 2. The contents of flavonoids in tartary buckwheat in response to salt stress (μg/g dry weight).
Table 2. The contents of flavonoids in tartary buckwheat in response to salt stress (μg/g dry weight).
SamplesCatechinChlorogenic AcidCaffeic Acidp-Coumaric AcidBenzoic AcidRutinQuercetinKaempferol
Con.1D12.04 ± 2.36a625.66 ± 28.28cd46.32 ± 0.03ab12.73 ± 2.72abc176.59 ± 8.75c12115.99 ± 132.05f54.84 ± 10.24c34.28 ± 4.42bc
2D11.56 ± 0.72a596.23 ± 3.06c37.29 ± 6.07ab14.2 ± 2.06bc162.52 ± 8.19b10406.33 ± 41.03c36.55 ± 14.71bc32.74 ± 6.92ab
4D19.21 ± 1.99bc776.19 ± 26.62d49.46 ± 0.02b10.32 ± 0.72a144.29 ± 3.23a9506.46 ± 212.82a10.52 ± 4.88a39.45 ± 1.21bc
5D14.61 ± 0.07ab535.18 ± 18.28ab38.83 ± 1.13ab10.24 ± 0.27a148.43 ± 4.77a9906.37 ± 78.64b52.76 ± 6.32c26.81 ± 1.09a
Tr.1D13.62 ± 5.58ab543 ± 31.06b48.32 ± 5.86b15.85 ± 2.44c200.01 ± 3.07d13599.51 ± 112.49g22.76 ± 4.03ab38.02 ± 4.34bc
2D21.66 ± 2.00c516.64 ± 10.10ab44.26 ± 0.43ab10.61 ± 0.21ab184.06 ± 2.38c11743.61 ± 110.46e46.73 ± 19.37c32.67 ± 1.83ab
4D23.04 ± 3.32c639.35 ± 6.33c39.99±14.99ab11.6 ± 0.48ab146.85 ± 1.43a9975.95 ± 93.02b111.04 ± 11.74d41.14 ± 5.04c
5D13.47 ± 5.87ab503.38 ± 1.84a33.46 ± 9.32a13.07 ± 3.55abc167.23 ± 5.67b11370.11 ± 0.01d47.5 ± 15.36c37.12 ± 1.94bc
Values with different letters indicate a significant difference between control (Con.) and salt-treated (Tr.) samples at p < 0.05, applying a Duncan test (SPSS 22.0).
Table 3. The contents of carotenoid in tartary buckwheat in response to salt stress (μg/g dry weight).
Table 3. The contents of carotenoid in tartary buckwheat in response to salt stress (μg/g dry weight).
SampleLuteinZeaxanthin13Z-β-Caroteneα-CaroteneE-β-Carotene9Z-β-Carotene
Con.1D244.38 ± 7.96a1.01 ± 0.25a38.58 ± 4.24a5.46 ± 0.46a193.78 ± 9.61b26.63 ± 1.88ab
2D236.83 ± 7.12a1.16 ± 0.19a36.34 ± 1.91a4.32 ± 0.59a172.21 ± 12.27ab22.57 ± 0.73a
4D299.29 ± 20.97b2.68 ± 0.41d47.27 ± 2.02b4.88 ± 0.09a235.54 ± 25.51c29.11 ± 3.31b
5D225.37 ± 17.38a1.80 ± 0.15bc33.19 ± 1.53a4.61 ± 0.72a156.30 ± 14.33a21.86 ± 2.42a
Tr.1D242.01 ± 14.83a1.42 ± 0.25abc38.46 ± 0.67a4.86 ± 0.29a179.58 ± 19.13ab25.89 ± 2.49ab
2D253.28 ± 36.49a1.27 ± 0.13ab37.51 ± 5.69a5.11 ± 1.57a177.76 ± 25.53ab26.02 ± 6.26ab
4D249.16 ± 31.94a3.24 ± 0.55e37.24 ± 1.20a5.54 ± 0.19a164.47 ± 23.28ab29.88 ± 3.85b
5D241.95 ± 9.48a1.89 ± 0.22c34.61 ± 0.52a5.02 ± 0.06a157.60 ± 11.81a26.12 ± 2.43ab
Values with different letters indicate a significant difference between control (Con.) and salt-treated (Tr.) samples at p < 0.05, applying a Duncan test (SPSS 22.0).

Share and Cite

MDPI and ACS Style

Ma, W.; Kim, J.K.; Jia, C.; Yin, F.; Kim, H.J.; Akram, W.; Hu, X.; Li, X. Comparative Transcriptome and Metabolic Profiling Analysis of Buckwheat (Fagopyrum Tataricum (L.) Gaertn.) under Salinity Stress. Metabolites 2019, 9, 225. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo9100225

AMA Style

Ma W, Kim JK, Jia C, Yin F, Kim HJ, Akram W, Hu X, Li X. Comparative Transcriptome and Metabolic Profiling Analysis of Buckwheat (Fagopyrum Tataricum (L.) Gaertn.) under Salinity Stress. Metabolites. 2019; 9(10):225. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo9100225

Chicago/Turabian Style

Ma, Weibo, Jae Kwang Kim, Caihua Jia, Feifan Yin, Hyo Jin Kim, Waheed Akram, Xuebo Hu, and Xiaohua Li. 2019. "Comparative Transcriptome and Metabolic Profiling Analysis of Buckwheat (Fagopyrum Tataricum (L.) Gaertn.) under Salinity Stress" Metabolites 9, no. 10: 225. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo9100225

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop