Agronomic and metabolomics analysis of rice-Tartary buckwheat (Fagopyrum tataricum Gaertn) bred by hybridization

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Scientific Reports volume 12, Article number: 11986 (2022 ) Cite this article Chinese Buckwheat Rice

Agronomic and metabolomics analysis of rice-Tartary buckwheat (Fagopyrum tataricum Gaertn) bred by hybridization | Scientific Reports

Tartary buckwheat (TB) is an edible pseudocereal with good health benefits, but its adhering thick shell and bitter taste inhibit its consumption. In this study, the first hybrid rice-Tartary buckwheat (RTB) variety Mikuqiao18 (M18), bred by the pedigree selection of crossbreeding ‘Miqiao’ (MQ) with ‘Jingqiaomai2’ (JQ2), was selected for an agronomic and metabolomics analysis. Compared with JQ2, M18 demonstrated a significantly lower yield per plant owing to the decreased grain weight and similar full-filling grain number per plant. However, M18 had a similar kernel weight per plant because of the thinner shell. The sense organ test suggested that M18 had higher taste quality regardless of partial replacement of rice through the improvement of preponderant indicators related to cereal taste quality, including lower values of total protein, albumin, glutelin, globulin, pasting temperature, cool paste viscosity, and setback. Meanwhile, M18 contained high levels of flavonoids, including rutin and quercetin, but presented a positive summary appraisal of cooking with 25% rice. Additionally, 92 metabolites were positively identified by GC–MS, including 59 differentially expressed metabolites (DEMs) between M18 and JQ2. Typically, M18 exhibited lower levels of 20 amino acids and higher levels of 6 sugars and 4 polyols. These DEMs might partly explain the superior eating quality of M18. In addition, M18 was abundant in 4-aminobutyric acid, which is beneficial to human health. The current findings offer a theoretical foundation for breeding rice-Tartary buckwheat with high yield and quality and promoting the cultivation and consumption of rice-Tartary buckwheat as a daily functional cereal.

Tartary buckwheat (TB, Fagopyrum tataricum Gaertn) and common buckwheat (F. esculentum Moench), native to southwest China and the Himalayas, are two dominant cultispecies of Fagopyrum Mill1,2. Common buckwheat is called ‘Tianqiao’ due to its good taste and has been suggested as a source of rutin in the daily diet for a long time3. Conventional Tartary buckwheat (CTB) is called ‘Kuqiao’ due to its bitter taste. Unlike common buckwheat, the shell of CTB is thick and adheres to the testa layer, contributing to increasing the difficulty of dehulling1,4. Traditional mechanical methods for the dehulling of CTB are complicated. Before dehulling, stewing and drying are two necessary steps, which not only raise the cost of processing but also lower the active ingredients of caryopsis4,5. Therefore, CTB is not as popular for daily consumption around the world despite containing several times more rutin than common buckwheat6,7,8,9,10.

‘Miqiao’ (MQ) is a special natural germplasm of TB and is also known as rice-Tartary buckwheat (RTB) given its thin shell for ready dehulling and cooking, similar to rice11. MQ is a local variety in Yunnan Province, China, with a long vegetative period and low yield characteristics. The poor agronomic traits of MQ were confirmed to be inherited when planted in many TB production areas11,12. In RTB breeding, genetic improvement for ecological adaptability and yield is recommended. Initially, Wang and Campbell reported the hybridization of RTB and found that the difficulty for crossbreeding was derived from the tiny flowers and low seed fertility11. Then, hot-water emasculation for TB was utilized in the crossbreeding of RTB13. Nonetheless, these hybrid RTB germplasms are rarely cultivated around the world due to defects such as poor ecological adaptability and low yield.

TB is classified as a cereal considering many similarities to other cereals. Starch and protein are the dominant storage components in CTB seeds. As the starch pasting viscosity and physicochemical properties of CTB are unique and different from most cereals, it has been suggested as an available source of retrograded starch14. Meanwhile, CTB seeds are abundant in proteins that are important active ingredients for reducing plasma total cholesterol15. CTB is famous for its high levels of flavonoids, which are recognized as the major active components for health care uses16,17,18. The cereal quality of RTB directly affects food development and processing. Compared with the CTB cultivar ‘Jinqiaomai2’ (JQ2), MQ exhibits similar levels of storage components, such as starch and proteins, while presenting significantly lower levels of flavonoids12. Considering that JQ2 is a high yield CTB cultivar with high levels of flavonoids, a series of RTB lines were created by crossing JQ2 with MQ19,20. Among them, Mikuqiao18 (M18) was identified and authorized as a high-yield variety. Therefore, the RTB variety M18 was selected to investigate the yield, major storage components, metabolites, starch pasting viscosity, and taste quality. The current findings will provide fresh insights into the ecological adaptability and cereal quality of RTB to popularize its planting, processing, and consumption as a daily food.

An effective method for TB hybridization was reported, and a considerable number of hybridization combinations were performed in RTB breeding19. Afterwards, hundreds of RTB lines were planted for pedigree selection. Among them, an excellent line from crossbreeding JQ2 with MQ, designated Mikuqiao18 (M18), was identified and authorized as a high-yield RTB variety by 6 generations of pedigree selection. M18 typically demonstrated an easy dehulling property inherited from its mother parent MQ, and it displayed a visible shell phenotype differing from that of JQ2. Notably, the shell of M18 was thin and partly cracked without a longitudinal furrow that could not cover the kernel completely (Fig. 1a).

Seed phenotype (a), grain weight (b) and weight ratio of kernel and shell (c) in grains of Tartary buckwheat. The photo was spliced by PowerPoint in Microsoft Office 2003. Blue arrows point to the longitudinal furrows of CTB seed, and red arrows point to the cracked shell of RTB seed. JQ2 indicates Jinqiaomai 2, M18 indicates Mikuqiao 18. Bar = 2 cm. * and ** indicate significant differences at 0.05 and 0.01, respectively, according to independent samples t-test (n = 3).

The RTB variety MQ exhibited a long vegetative period and low yield characteristics when planted in many production areas11,12. As M18 presented a longer period of growth, inherited from its mother parent, M18 was sown 10 days before JQ2 and achieved a similar maturing status in this experiment. At harvest, M18 exhibited extremely significantly higher plant height than JQ2 at harvest in two tested years12. Although the effective branching number was undifferentiated, the grain number per plant in M18 was extremely significantly increased. However, no significant difference in the full-filling grain number per plant between JQ2 and M18 was observed owing to the higher flat and empty grains in M18 (Table 1).

Compared with JQ2, the 1000-grain weight of M18 was extremely significantly lower (Fig. 1b), while the kernel weight ratio in grain (85.0%) was extremely significantly higher, and the shell weight ratio in grain (15.0%) was extremely significantly lower (Fig. 1c). M18 obtained a parallel yield per plant, although the yield per plant was remarkably lower (Table 1).

To identify the cooking and eating quality of RTB, the starch pasting viscosity parameters were detected using a rapid visco analyser (RVA). M18 demonstrated a significantly lower value of pasting temperature (PT) (61.0 °C) and an extremely significant lower cool paste viscosity (CPV) value, leading to significantly lower values of setback (SB) and consistency (CS) compared with JQ2 (Table 2). Nevertheless, no significant difference was detected in the peak viscosity (PV), hot paste viscosity (HPV), and breakdown (BD) between M18 and JQ2.

Currently, people prefer to cook coarse cereals with rice together as a staple food to achieve nutritional and health benefits. Therefore, it might be acceptable to consume RTB kernel with rice daily. To further evaluate the taste of cooked kernel of RTB, the appearance, aroma, viscosity, and taste of different formulated M18 with rice were investigated according to the method of rice taste technology and compared with those of the famous rice variety ‘Koshihikari’. Consistently, both M18 and JQ2 exhibited lower values of appearance, aroma, viscosity, and taste when individually cooked (Table 3). Nevertheless, M18 exhibited higher values of appearance, taste, viscosity, and summary than JQ2. Moreover, the positive appraisal of the sense organ test was identified when TB kernel was mixed with rice. Note that the positive summary appraisal was tested when the additive amount of M18 reached 75%.

As storage components are closely related to crop quality, the major storage components in kernels of M18 were determined. Starch was the dominant storage component in TB seeds (Fig. 2). Compared with JQ2, M18 demonstrated similar levels of starch consisting of amylose and amylopectin but significantly lower levels of total protein due to the significantly decreased albumin, glutelin, and globulin (Fig. 2b). Flavonoids are the dominant biologically active components in TB seeds, such as rutin and quercetin. Although no significant difference in total flavonoids was detected between M18 and JQ2, M18 presented extremely significantly lower contents of rutin and significantly higher contents of quercetin (Fig. 2c).

The levels of starch (a), proteins (b) and flavonoids (c) in grain of Tartary buckwheat. * and ** indicate significant differences at 0.05 and 0.01, respectively, according to independent samples t-test (n = 3).

To identify the metabolites in M18, the GC–MS metabolomics platform was used. A total of 92 metabolites were positively identified, mainly including 27 amino acids, 26 organic acids, 10 sugars, 6 phosphoric acids, 6 polyols, and 3 fatty acids (Fig. 3a). The principal component analysis (PCA) indicated that the metabolomes were different among M18, JQ2, and the quality control, with R2X[1] and R2X[2] values of 0.59 and 0.125, respectively (Fig. 3b). OPLS-DA was performed to maximize the distinction between M18 and JQ2, with R2X[1] and R2X[2] values of 0.59 and 0.101, respectively (Fig. 3c), indicating that the established model was appropriate for differential expressed metabolites (DEMs) screening.

The groups (a), PCA (b) and PLS-DA (c) of metabolites in Tartary buckwheat grains by GC–MS. Multivariate statistical analysis was performed using the software package SIMCA-P (V13.0) and R language ropls package.

The heatmap exhibited a great difference in DEMs between M18 and JQ2 (Fig. 4a). A total of 22 elevated and 37 reduced DEMs were identified in M18 (p value ≤ 0.05 and VIP value ≥ 1) compared with JQ2 (Fig. 4a, Table S1). Among them, M18 demonstrated elevated levels of 4-hydroxyproline (2.69-fold), l -homoserine (2.49-fold), and l -methionine (1.3-fold), but reduced levels of 20 amino acids, such as glycine (0.53-fold), l -tryptophan (0.53-fold), tyrosine (0.53-fold), and ornithine (0.53-fold), indicating that the metabolism of amino acids in M18 was weaker than that in JQ2 (Fig. 4b, Table S1). By contrast, M18 showed elevated levels of 6 sugars (arabinose (1.79-fold), fucose (1.25-fold), galactose (1.68-fold), glucose (1.47-fold), isomaltose (1.30-fold), and ribose (1.39-fold)) and 4 polyols (erythritol (2.73-fold), myo-inositol (1.31-fold), threitol (2.48-fold), and xylitol (2.02-fold)), implying that the metabolism of sugars and polyols was improved in M18 (Fig. 4b, Table S1). Meanwhile, M18 displayed elevated levels of metabolites associated with the tricarboxylic acid cycle, consisting of pyruvic acid (1.76-fold), a-ketoglutaric acid (1.62-fold), and malic acid (1.82-fold). Additionally, M18 showed a 2.3-fold increase in 4-aminobutyric acid compared with JQ2 (Table S1).

The heatmap (a) of DEMs and the boxplot-visualizations of partial DEMs (b) by GC–MS. The scale of data set and the bidirectional clustering of samples and metabolites by hierarchical clustering were performed using the pheatmap package in R (V1.0.8). S1 indicates arabinose, S2 indicates fucose, S3 indicates galactose, S4 indicates glucose, S5 indicates isomaltose, S6 indicates ribose, P1 indicates erythritol, P2 indicates myo-inositol, P3 indicates threitol, P4 indicates xylitol, A1 indicates 2-aminoadiopic acid, A2 indicates 2-aminobutyric acid, A3 indicates beta-Alanine, A4 indicates glycine, A5 indicates l -asparagine, A6 indicates l -aspartic acid, A7 indicates l -cysteine, A8 indicates l -glutamic acid, A9 indicates l -glutamine, A10 indicates l -isoleucine, A11 indicates l -lysine, A12 indicates l -phenylalanine, A13 indicates l -proline, A14 indicates l -serine, A15 indicates l -threonine, A16 indicates l -tryptophan, A17 indicates l -valine, A18 indicates ornithine, A19 indicates pyroglutamic acid, A20 indicates tyrosine. (n = 6).

The KEGG analysis demonstrated that the DEMs covered 50 pathways or metabolisms. The most diverse metabolic pathways included aminoacyl-tRNA biosynthesis; arginine and proline metabolism; glycine, serine and threonine metabolism; and alanine, aspartate, and glutamate metabolism (Table S2).

Person’s correlation coefficient analysis revealed 1434 significant correlation coefficients (p < 0.01, r2 ≥ 0.49, FDR ≤ 0.01) between DEMs, including 822 positive and 612 negative significant correlations (Fig. 5, Table S3). Among them, l -homoserine, ornithine, erythritol, threitol, xylitol, and kaempferol were the most associated metabolites and all contained 57 cytonodes (Table S3). Negative significant correlations were observed between metabolites from different metabolic pathways, such as xylitol and phosphoric acid (− 0.993), xylitol and glycine (− 0.989), and xylitol and kaempferol (− 0.979). Meanwhile, positive significant correlations were discovered between metabolites belonging to the same metabolic pathways, such as l -aspartic acid and glycine (0.997), threitol and erythritol (0.995), and glucose and galactose (0.994) (Table S3).

Pearson’s correlation coefficient analysis of DEMs by GC–MS. The correlation coefficient was calculated by Pearson correlation coefficient using the cor() function in R (V3.1.3). The significance statistical test of metabolite correlation analysis was conducted using the cor.test() function in R language package (|R| ≥ 0.7 and FDR p value ≤ 0.05).

Breeding adaptive RTB cultivars with easy dehulling seeds through hybridization using ‘Miqiao’ as a parent has long been considered. However, the crossbreeding of RTB is challenging due to the tiny flowers and low seed fertility11. Although the success of hot-water emasculation for RTB has been reported13, crossbreeding RTB is not popular with other research teams given some unknown mechanism that interferes with fertilization. Thus, RTB has not been widely planted around the world. Our previous work reported the success of hybridization of RTB by artificial emasculation on newly opened flowers using sharp materials such as toothpicks and tweezers19. Hereafter, hundreds of lines were created by crossing high yield CTB varieties with RTB germplasm MQ. Finally, M18, crossbred by MQ with JQ2, was identified and authorized as a high-yield RTB variety by 6 generations of pedigree selection. M18 exhibited an easy dehulling property inherited from its mother parent MQ12, and its shell was thin and accounted for only 15% of the grain weight. Additionally, M18 exhibited a significantly higher grain number per plant than JQ2, suggesting that the ‘sink’ was improved by crossbreeding. However, the full-filling grain number per plant of M18 was not significantly increased, and the grain weight was significantly decreased, leading to a significantly lower yield per plant compared with JQ2. Moreover, the kernel weight per plant of M18 was similar to that of JQ2 because of the significantly thinner shell. Although M18 showed a longer period of growth than JQ2, the poor yield characteristics of MQ were overcome in M18 as similar yield levels were achieved in two tested years compared with the high-yield variety JQ212. Therefore, M18 was proposed as a potential RTB variety for cultivation in TB-producing areas.

Buckwheat is classified as a pseudocereal for its similar chemical composition as conventional cereals. CTB contains several times more flavonoids than common buckwheat and has long been suggested as a dietary source of rutin21. CTB is rarely directly cooked for daily consumption around the world because of its bitter taste, despite its excellent health benefits22. Therefore, the cooking and eating quality remains to be improved, even though the problem of dehulling is not present in RTB. Nonetheless, the taste quality of crops is a complicated trait controlled by multitudinous factors, such as seed storage components and metabolites.

Starch is the dominant storage component in CTB seeds and is composed of amylose and amylopectin. The pasting viscosity of cereal is dramatically influenced by the starch level, especially the amylose level. RVA has been widely employed to evaluate the starch pasting properties, cooking and taste quality of cereals23,24,25. In this study, the levels of starch, amylose, and amylopectin in M18 were similar to those in JQ2. Hence, no significant difference in the values of RV, HPV and BD between M18 and JQ2 was detected by RVA. In addition, the values of SB and CPV are significantly negative correlated with the cooking and taste quality of cereals26. Consistent results were revealed in TB, as higher values of appearance, viscosity, and taste and lower values of SB and CPV were discovered in M18 compared with JQ2.

Proteins are the major nutritional component in CTB seeds. High levels of proteins are proposed in high-quality breeding of cereals, while excessive levels of proteins negatively affect taste quality27. This relationship probably remained in TB considering that lower levels of most proteins (albumin, glutelin, prolamine) were detected in M18 while higher taste quality was found compared with JQ2. Therefore, the contradiction between nutritional quality and taste quality might exist and should be further improved in the crossbreeding of RTB for high yield and quality.

Flavonoids, especially rutin and quercetin, are the major active components in CTB seeds16,17,18. Compared with JQ2, M18 exhibited extremely significantly lower levels of rutin and significantly higher levels of quercetin. No significant difference in total flavonoid levels between JQ2 and M18 was detected, although the total rutin and quercetin in M18 was slightly lower than that of JQ2. Thus, the RTB variety M18 had a high capacity for health efficacy, similar to JQ2. High levels of rutin in CTB can be hydrolysed into quercetin during the cooking process and have been verified to be responsible for its bitter taste21,28. Furthermore, 87.5% of the rutin RTB flour degraded 10 min after the addition of water, while the rutin in RTB dehulled whole seed (whole groats) was stable even when it was immersed in water21. Currently, an increasing number of people prefer to eat rice mixed with coarse cereals daily to achieve health benefits. In this study, the sense organ test of palatability characteristics was investigated among different ratios of dehulled whole seeds of RTB mixed with rice. The positive appraisal of the sense organ test, including appearance, aroma, viscosity, and summary, was identified at the maximum level of 75% M18 mixed with 25% rice variety ‘cv. Koshihikari’. Therefore, it is feasible to directly utilize the dehulled whole seed of RTB cooked with rice.

Metabolomics is a powerful strategy for determining metabolites and explaining macroscopic biological phenomena, providing new insight into crop improvement29. In this study, a total of 59 DEMs were identified between M18 and JQ2 by GC–MS. Amino acids are not only precursors for protein synthesis but also crucial nutritional components for cereals. Compared with JQ2, M18 had only three elevated amino acids (4-hydroxyproline, l -homoserine and l -methionine) belonging to the aspartic acid family. Most amino acids of the five families were reduced in M18, including the a-ketoglutaric acid family, pyruvic acid family, glycerate-3-phosphate family and aromatic amino acid family amino acids. Nonetheless, a-ketoglutaric acid and pyruvic acid were elevated in M18, suggesting that the synthesis of amino acids was limited by the deficiency of the nitrogen source instead of the carbon source. This was consistent with the results of the significantly lower protein contents in M18, which weakened the nutritional quality but might contribute to the better taste quality. M18 is a 4-aminobutyric acid-rich variety, the level of which is 2.3-fold that of JQ2, bringing many benefits such as improving sleep quality and lowering blood pressure30.

Tartary buckwheat is a well-known functional cereal called ‘kuqiao’ owing to its bitter taste2,22. Currently, the demand for high-quality Tartary buckwheat is gradually increasing in modern life. Sugars are not only energy sources but also important sweeteners in cereals. Compared with JQ2, M18 had 6 elevated sugars. Among them, arabinose was the most elevated sugar in M18, and it benefits the hypoglycemic effect by reducing disaccharide hydrolysis and lipid lowering by reducing the synthesis of fat31. Meanwhile, M18 showed elevated fucose, the conjugation of which plays an important role in cancer biology32. Polyols are well-known sugar-free sweeteners. Through GC–MS, 6 polyols were identified, while M18 presented 4 elevated polyols. Among them, the levels of erythritol, threitol, and xylitol were elevated more than twofold. As reported in rice, polyols vary in cultivars with different tastes33. Hence, these differentially accumulated metabolites might partly explain the superior eating quality of M18 by the sense organ test.

The cereal quality of Tartary buckwheat is a complex system consisting of nutritional quality, health care quality, and cooking and taste quality. It is difficult for traditional testing methods to comprehensively evaluate cereal quality32. In this study, the combined metabolomics and traditional testing methods provided an effective strategy for a comprehensive assessment of cereal quality, which would significantly accelerate the high-quality breeding process and promote the consumption of RTB as a daily functional cereal.

The Tartary buckwheat cultivars ‘Jingqiaomai2’ (JQ2) and ‘Mikuqiao18’ (M18) were planted at the Experimental Station of the Research Center of Buckwheat Engineering Technology in Guizhou Province, China (1146.0 m, 26°50′ N, 106°58′ E) during the autumn growing season in 2018 and 2019. The average temperature was 21.79 °C/14.91 °C and 24.00 °C/16.55 °C, and the amount of precipitation was 329 mm and 244.2 mm during the growing period in 2018 and 2019, respectively (Supplement Table 2). JQ2 is a CTB high yield variety authorized by the National Identification Committee of Minor Grain and Bean in China. M18 is a high-yield RTB variety obtained by crossbreeding Miqiao with JQ2.

M18 was sown in early August, and JQ2 was sown in mid-August. However, they were synchronously harvested in late November. The soil in the test pot was yellow loam containing 12.15 g kg−1 organic matter, 0.81 g kg−1 total nitrogen, 96.37 mg kg−1 alkaline nitrogen, 12.42 mg kg−1 valid phosphorus, and 118.57 mg kg−1 valid potassium. The fertilizer included 60 kg N ha−1, 60 kg P2O5 ha−1, and 30 kg K2O ha−1, and was applied as a base fertilizer in each plot before sowing.

A randomized block design with three replications was used in each region. The area for each test plot was 2 m × 5 m, with row width of 33 cm. The planting density was approximately 1,000,000 plants per hectare. Consistent management measures were taken, and there were no significant pests, diseases, or weeds during the growth period.

A total of 12 separate plants from each plot in the Experiment Station of the Research Center of Buckwheat Engineering Technology were sampled at harvest to measure the plant height, effective branching number, grain number per plant, full-filling grain number per plant, grain weight, kernel weight, shell weight, yield per plant and kernel yield per plant. The rest of the plants in each plot were harvested to investigate the yield and quality.

A total of 3000 mg of sample and 25 mL of distilled water were added to a tube and detected by a rapid visco analyser (RVA, Newport Scientific, Australia). Parameters were determined according to Tone et al.30. Then, the values of PV, HPV, CPV, BD, SB, CS, and PT were analysed by software kit of Thermal Cycle Window according to Tone et al.34.

Regarding the identification of taste quality, M18 and JQ2 were mixed with rice variety ‘Koshihikari’ at 0, 25%, 50%, 75% and 100% ratios, respectively, and cooked in an electric rice cooker (MB-40EASY202, Midea, China). Twenty-two trained panelists participated in evaluating the taste quality from different aspects of appearance, aroma, taste, and viscosity, compared with Koshihikari and provided a summary appraisal following the method of Zhao et al.35.

TB kernels were previously crushed and sieved with a 100-mesh sieve for the further measurement of storage components.

A total of 100 mg of sample was eluted with 2 mL of 95% ethanol to remove soluble sugars three times. The details for amylose and amylopectin measurement were performed using the methods of Liang et al.20. The starch content was calculated with amylopectin and amylose.

A total of 100 mg of sample was used to extract albumin, glutelin, prolamine and globulin by adding 5 mL of 10 mmol/L Tris–HCl (pH 7.5), 5 mL of 0.5 mol/L NaCl, 5 mL of 60% normal propyl alcohol and 5 mL of mixed liquor (containing 0.5% sodium tartrate, 0.24% copper sulfate, 1.68% KOH and 50% normal propyl alcohol) three times. The albumin, glutelin, prolamine and globulin measurement were performed following the methods of Wang et al.32. The total protein content was calculated with the total albumin, glutelin, prolamine and globulin.

A total of 100 mg of sample was used to extract flavonoids by adding 2 mL of 80% ethanol to a hot bath at 70 °C for 5 h and ultrasonic extraction for 10 min. The flavonoids measurement was conducted following the methods of Wang et al.36.

A total of 100 mg sample was used to extract rutin by adding 20 mL methanol for ultrasonic extraction at 30 °C. The rutin measurement was achieved using the methods of Zhang et al.37.

A total of 600 mg of sample was used to extract quercetin by adding a moderate amount of ether 3 times for ultrasonic extraction at 30 °C for 30 min and 2 mL of ethyl acetate 5 times for ultrasonic extraction at 30 °C for 30 min to collect the supernatant. Then, the supernatant was dried at 80 °C and washed with sterile water. The sediment was collected by centrifugation and redissolved to a constant volume of 10 mL by ethanol. The quercetin measurement was performed following the methods of Lei et al.38.

A total of 100 mg of fresh sample was ground by a high flux organization grinding apparatus under liquid nitrogen conditions. The sample preparation was conducted following the methods of Lisec et al.39 and Sangster et al.40. Metabolites were detected using GC (Agilent GC 7890A-5975C, USA) using a a 280 °C injection temperature, 1:20 split ratio, 150 °C interface, 230 °C ion source and temperature control (60 °C for 2 min, 10 °C/min rate up to 300 °C and holding there for 5 min). Then, mass spectrometry was recorded from 35 to 750 (m/z).

The data were processed by G1701 MSD ChemStation and R (v3.3.2) and analysed by AMDIS (v2.71) using NIST, Wiley Registry and GMD following the method of Smith et al.41. Multivariate data analysis was performed by SIMCA-P (v13.0, Umetrics, Sweden) and the R package ropls. The PCA and OPLS-DA were validated by the method of Thévenot et al.42. The potential biomarkers and variable contributions were evaluated through the analysis of VIP values (VIP > 1) and ANOVA (p value < 0.05) using SPSS (v19.0, IBM Corporation, Armonk, NY, USA). The heatmap was performed using the pheatmap package in R (V1.0.8, https://CRAN.R-project.org/package=pheatmap)43 KEGG analysis of DEMs was performed using the Kyoto Encyclopedia of Genes and Genomes database (http://www.kegg.jp/kegg/kegg1.html)44,45.

MS-Excel was used to classify Data. The means and differences of the data from three replications were determined by SPSS 17.0 (SPSS Inc., Chicago, IL, USA).

The permission of this study was obtained from Research Center of Buckwheat Engineering Technology in Guizhou Province. The experimental research and field studies on plants (either cultivated or wild) in this study were complied with relevant institutional, national, and international guidelines and legislation.

The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.

Zhang, Q. & Xu, J. G. Determining the geographical origin of common buckwheat from China by multivariate analysis based on mineral elements, amino acids and vitamins. Sci. Rep. 7(1), 9696. https://doi.org/10.1038/s41598-017-08808-y (2017).

Article ADS PubMed PubMed Central Google Scholar

Weisskopf, A. & Fuller, D. Q. Buckwheat: Origins and development. In Encyclopedia of Global Archaeology 1025–1028 (Springer, 2014). https://doi.org/10.1007/978-1-4419-0465-2_2310.

Couch, J. F., Naghski, J. & Krewson, C. F. Buckwheat as a source of rutin. Science 103(2668), 97–198. https://doi.org/10.1126/science.103.2668.197 (1946).

Zhu, X., Fan, W., Li, Z. & Guo, W. Test and development of a non-thermal tartary buckwheat huller. J. Chin. Agric. Univ. 22(12), 146–155. http://en.cnki.com.cn/Article_en/CJFDTotal-NYDX201712018.htm (2017).

Liu, Y., Du, W. & Wu, Y. Experimental study on the optimization of the hulling method and main parameters of Tartarian buckwheat. J. Agric. Mech. Res. 12, 137–139. http://en.cnki.com.cn/article_en/cjfdtotal-njyj200812044.htm (2008).

Wei, Y. M. Buckwheat remains from the Late Neolithic Site of Donghuishan, Gansu Province, China. Cereal Chem. 96, 332–337. https://doi.org/10.1002/cche.10130 (2019).

Fabjan, N. et al. Tartary buckwheat (Fagopyrum tataricum Gaertn.) as a source of dietary rutin and quercitrin. J. Agric. Food Chem. 51(22), 6452–6455. https://doi.org/10.1021/jf034543e (2003).

Article CAS PubMed Google Scholar

Ladan, M. K., Straus, J., Benković, E. T. & Kreft, S. FT-IR-based method for rutin, quercetin and quercitrin quantification in different buckwheat (Fagopyrum) species. Sci. Rep. 7, 7226. https://doi.org/10.1038/s41598-017-07665-z (2017).

Article ADS CAS Google Scholar

Juan Antonio, G. B. & Henryk, Z. Buckwheat as a functional food and its effects on health. J. Agric. Food Chem. 63(36), 7896–7913. https://doi.org/10.1021/acs.jafc.5b02498 (2015).

Tomotake, H. et al. Hypolipidemic activity of common (Fagopyrum esculentum Moench) and tartary (Fagopyrum tataricum Gaertn.) buckwheat. J. Sci. Food Agric. 95(10), 1963–1967. https://doi.org/10.1002/jsfa.6981 (2015).

Article CAS PubMed Google Scholar

Wang, Y. & Campbell, C. G. Tartary buckwheat breeding (Fagopyrum tataricum L. Gaertn.) through hybridization with its Rice-Tartary type. Euphytica 156(3), 399–405. https://doi.org/10.1007/s10681-007-9389-3 (2007).

Liao, K., Wang, Y., Yu, W., Wei, C. & Liang, C. Characteristics of starch synthesis, yield and quality provide insights into rice-Tartary Buckwheat utilization and improvement. Agric. Biotechnol. 2(9), 31-34+38. https://doi.org/10.19759/j.cnki.2164-4993.2020.02.009.html (2020).

Mukasa, Y., Suzuki, T. & Honda, Y. Suitability of rice-Tartary buckwheat for crossbreeding and for utilization of rutin. Jpn. Agric. Res. Q. 43(3), 199–206. https://doi.org/10.6090/jarq.43.199 (2009).

Gao, J. et al. Tartary buckwheat (Fagopyrum tataricum Gaertn.) starch, a side product in functional food production, as a potential source of retrograded starch. Food Chem. 190, 552–558. https://doi.org/10.1016/j.foodchem.2015.05.122 (2016).

Article CAS PubMed Google Scholar

Zhang, C., Zhang, R., Li, Y. M., Liang, N. & Chen, Z. Y. Cholesterol-lowering activity of Tartary buckwheat protein. J. Agric. Food Chem. 65(9), 1900–1906. https://doi.org/10.1021/acs.jafc.7b00066 (2017).

Article CAS PubMed Google Scholar

Tao, B. et al. Antioxidant and antidiabetic properties of Tartary buckwheat rice flavonoids after in vitro digestion. J. Zhejiang Univ. Sci. B 17(12), 941. https://doi.org/10.1631/jzus.B1600243 (2016).

Lee, L. S. et al. Contribution of flavonoids to the antioxidant properties of common and tartary buckwheat. J. Cereal Sci. 68, 181–186. https://doi.org/10.1016/j.jcs.2015.07.005 (2016).

Song, C. et al. Changes in seed growth, levels and distribution of flavonoids during tartary buckwheat seed development. Plant Prod. Sci. 19(4), 1–10. https://doi.org/10.1080/1343943X.2016.1207485 (2016).

Chen, Q., Chen, Q., Shi, T., Huang, K., Liang, L., Rao, Q., Tang, L. & Zhang, Q. Inheritance of Tartary buckwheat thick shell character and its relationships with yield factors. Crops 2, 27–31. http://en.cnki.com.cn/Article_en/CJFDTOTAL-ZWZZ201502006.htm (2015).

Liang, C. et al. Characteristics of the grain-filling process and starch accumulation of high-yield common buckwheat “cv. Fengtian 1” and Tartary buckwheat “cv. Jingqiao 2”. Cereal Res. Commun. 3(44), 393–403. https://doi.org/10.1556/0806.44.2016.005 (2016).

Mukasa, Y., Suzuki, T. & Morishita, T. Rutin degradation in dehulled seeds of rice-Tartary buckwheat. Fagopyrum 29, 39–43 (2009).

Suzuki, T. et al. Breeding of “Manten-Kirari”, a non-bitter and trace-rutinosidase variety of Tartary buckwheat (Fagopyrum tataricum Gaertn.). Breed. Sci. 64(4), 344–350. https://doi.org/10.1270/jsbbs.64.344 (2014).

Article CAS PubMed PubMed Central Google Scholar

Ragaee, S. & Abdel-Aal, E. S. M. Pasting properties of starch and protein in selected cereals and quality of their food products. Food Chem. 95(1), 9–18. https://doi.org/10.1016/j.foodchem.2004.12.012 (2006).

Balindong, J. L. et al. Rice grain protein composition influences instrumental measures of rice cooking and eating quality. J. Cereal Sci. https://doi.org/10.1016/j.jcs.2017.09.008 (2017).

Liu, S. et al. Behaviors of starches evaluated at high heating temperatures using a new model of Rapid Visco Analyzer—RVA 4800. Food Hydrocoll. 94, 217–228. https://doi.org/10.1016/j.foodhyd.2019.03.015 (2019).

Balet, S., Guelpa, A., Fox, G. & Manley, M. Rapid visco analyser (RVA) as a tool for measuring starch-related physiochemical properties in cereals: A review. Food Analyt. Methods 12(3), 1–17. https://doi.org/10.1007/s12161-019-01581-w (2019).

Zhang, X., Shi, L., Ding, D., Wang, S. & Cui, J. The relationship between rice protein related character and the RVA characteristic profile and palatability character. Food Sci. Technol. 39(10), 188–191. https://doi.org/10.13684/j.cnki.spkj.2014.10.040.html (2014).

Lukšič, L. et al. Hydrothermal treatment of Tartary buckwheat grain hinders the transformation of rutin to quercetin. J. Cereal Sci. 72, 131–134. https://doi.org/10.1016/j.jcs.2016.10.009 (2016).

Parray, J. A. & Shameem, N. Metabolomics for crop improvement: Quality and productivity. Sustain. Agric. https://doi.org/10.1016/B978-0-12-817109-7.00001-8 (2020).

Pol, K., Graaf, K. D., Bruin, D. D., Balvers, M. & Mars, M. The effect of replacing sucrose with l -arabinose in drinks and cereal foods on blood glucose and plasma insulin responses in healthy adults. J. Funct. Foods. 73, 104114. https://doi.org/10.1016/j.jff.2020.104114 (2020).

Keeley, T. S., Yang, S. & Lau, E. The diverse contributions of fucose linkages in cancer. Cancers 11(9), 1241. https://doi.org/10.3390/cancers11091241 (2019).

Article CAS PubMed Central Google Scholar

Gong, R. et al. Comparative metabolomics analysis reveals the variations of eating quality among three high-quality rice cultivars. Mol. Breed. 40, 112. https://doi.org/10.1007/s11032-020-01192-y (2020).

Shimada, M., Hasegawa, T., Nishimura, C., Kan, H. & Matsubayashi, T. Anti-hypertensive effect of gamma-aminobutyric acid (gaba)-rich chlorella on high-normal blood pressure and borderline hypertension in placebo-controlled double blind study. Clin. Exp. Hypertens. 31(4), 342–354. https://doi.org/10.1080/10641960902977908 (2009).

Article CAS PubMed Google Scholar

Tong, C. et al. Genetic diversity of amylose content and RVA pasting parameters in 20 rice accessions grown in Hainan, China. Food Chem. 161, 239–245. https://doi.org/10.1016/j.foodchem.2014.04.011 (2008).

Zhao, J., Kusutani, A., Cui, J., Chen, X., Liu, J. & Niu, J. Sense organ test method of palatability characteristics on Japonica varieties. Tianjin Agric. Sci. 1(9), 12–14. http://www.en.cnki.com.cn/Article_en/CJFDTotal-HNNB201603007.htm (2003).

Wang, Y. et al. The yield and quality of tartary buckwheat varieties and the response to low nitrogen. J. Guizhou Normal Univ. (Nat. Sci.) 6(35), 66–73. https://doi.org/10.16614/j.cnki.issn1004-5570.2017.06.006.html (2017).

Zhang, H., Guan, Z. & Tian, X. Determination of rutin in Fagopyrum by spectrometric method. Chin. J. Spectr. Lab. 3, 24–27. https://kns.cnki.net/KCMS/detail/detail.aspx?dbcode=CJFD&filename=GPSS603.009 (1996).

Lei, T., Meng, F. & Wu, Y. Quantitative determination of quercetin of cortex Juglandis mandshuricae in Mudanjiang. Lab. Sci. 15(1), 88–90. http://en.cnki.com.cn/Article_en/CJFDTotal-YSKT201201030.htm (2012).

Lisec, J., Schauer, N., Kopka, J., Willmitzer, L. & Fernie, A. R. Gas chromatography mass spectrometry-based metabolite profiling in plants. Nat. Protoc. 1(1), 387–396. https://doi.org/10.1038/nprot.2006.59 (2006).

Article CAS PubMed Google Scholar

Sangster, T., Major, H., Plumb, R., Wilson, A. J. & Wilson, I. D. A pragmatic and readily implemented quality control strategy for HPLC–MS and GC–MS-based metabonomic analysis. Analyst 131(10), 1075–1088. https://doi.org/10.1039/b604498k (2006).

Article ADS CAS PubMed Google Scholar

Smith, C. A., Want, E. J., O’Maille, G., Abagyan, R. & Siuzdak, G. XCMS: Processing mass spectrometry data for metabolite profiling using nonlinear peak alignment, matching, and identification. Anal. Chem. 78(3), 779–787. https://doi.org/10.1021/ac051437y (2006).

Article CAS PubMed Google Scholar

Thévenot, E. A., Roux, A., Xu, Y., Ezan, E. & Junot, C. Analysis of the human adult urinary metabolome variations with age, body mass index, and gender by implementing a comprehensive workflow for univariate and OPLS statistical analyses. J. Proteome Res. 14(8), 3322–3335. https://doi.org/10.1021/acs.jproteome.5b00354 (2015).

Article CAS PubMed Google Scholar

Kolde, R. pheatmap: Pretty Heatmaps.https://CRAN.R-project.org/package=pheatmap (2015).

Kanehisa, M. & Goto, S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 28, 27–30. https://doi.org/10.1093/nar/28.1.27 (2000).

Article CAS PubMed PubMed Central Google Scholar

Kanehisa, M., Sato, Y., Kawashima, M., Furumichi, M. & Tanabe, M. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 44, 457–462. https://doi.org/10.1093/nar/gkv1070 (2016).

This research was funded by the Co-project of Karst Research Center of Guizhou Province and National Natural Science Foundation of China (Grant number: U1812401), the National Natural Science Foundation of China (Grant number: 31660366 and 32160428) and Natural Science Foundation of Guizhou Province (QKH ZC[2019]2298). We thank the Fund for construction of the Key Laboratory for Conservation and Innovation of Buckwheat Germplasm in Guizhou (QJH KY Zi [2017]002) and the researchers worked in the comprehensive test stations of China Agriculture Research System for ecological adaptability test and the panelists for sense organ test.

These authors contributed equally: Yan Wang and Zhixiu Guan.

Research Center of Buckwheat Industry Technology, Guizhou Normal University, Guiyang, 550001, China

Yan Wang, Zhixiu Guan, Chenggang Liang, Kai Liao, Juan Huang, Chunyu Wei, Taoxiong Shi & Qingfu Chen

College of Bioengineering, Chengdu University, Chengdu, 610106, China

National Research and Development Center for Coarse Cereal Processing, Chengdu, 610106, China

You can also search for this author in PubMed Google Scholar

Y.W. and Z.G. analysed the results and drafted the manuscript. K.L., C.W., and D.X. conducted the experiment(s). J.H., T.X. and Q.C. partly funded the experiments and revised the manuscript. C.L. and Q.C. conceived the experiment(s). All authors reviewed the manuscript.

Correspondence to Chenggang Liang or Qingfu Chen.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Wang, Y., Guan, Z., Liang, C. et al. Agronomic and metabolomics analysis of rice-Tartary buckwheat (Fagopyrum tataricum Gaertn) bred by hybridization. Sci Rep 12, 11986 (2022). https://doi.org/10.1038/s41598-022-16001-z

DOI: https://doi.org/10.1038/s41598-022-16001-z

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Scientific Reports (Sci Rep) ISSN 2045-2322 (online)

Baked Cassia Tea Sign up for the Nature Briefing: Translational Research newsletter — top stories in biotechnology, drug discovery and pharma.

Agronomic and metabolomics analysis of rice-Tartary buckwheat (Fagopyrum tataricum Gaertn) bred by hybridization | Scientific Reports