Sucrose as a key nutritional marker distinguishing vegetable and grain soybeans, regulated by GmZF-HD1 via GmSPS17 in seeds

Abstract

Vegetable and grain soybeans are typically distinguished by harvest time and pod size, yet their nutritional differences are often overlooked in breeding programs. This study compared 10 varieties each of vegetable and grain soybeans to find key nutritional markers distinguishing them. Results showed that vegetable soybeans have higher concentrations of sucrose, total soluble sugar, and crude protein, along with lower concentrations of crude oil and total fatty acid. Specifically, vegetable soybeans contain a relatively higher amount of unsaturated fatty acids, particularly oleic acid, at green edible stages. Principal component analysis of 12 nutritional components revealed clear distinctions between vegetable and grain soybeans. Additionally, machine learning algorithms identified sucrose as the most critical nutritional marker for distinguishing these two types. Dynamic RNA-seq analysis combined with weighted gene co-expression network analysis identified a sucrose-related module, highlighting GmSPS17 as a predominant sucrose phosphate synthase encoding gene involved in sucrose accumulation in soybean seeds. Furthermore, we identified GmZF-HD1 as an upstream transcription factor regulating GmSPS17. Yeast one-hybrid, luciferase, and electrophoretic mobility shift assays confirmed that GmZF-HD1 directly activates GmSPS17 transcription. Overexpression experiments in hairy roots validated that GmZF-HD1 enhances GmSPS17 expression, thereby increasing sucrose accumulation. In summary, this study establishes sucrose as a key nutritional marker for distinguishing vegetable soybeans from grain soybeans and elucidates the GmZF-HD1-GmSPS17 regulatory pathway, providing valuable insights into sugar accumulation mechanisms and offering guidance for breeding high-sugar vegetable soybean varieties.

Acknowledgments

This study was funded by the National Key R&D Program of China (2021YFD1201103-03), the Strategic Priority Research Program of Chinese Academy of Sciences (XDA28070402), and Heilongjiang Provincial Natural Science Foundation Outstanding Young Scholar Fund (YQ2024D009).

Data availability

The RNA-Seq data presented in the study are deposited in the Genome Sequence Archive, accession number: PRJCA035555.

Conflict of interest statement

The authors declare that they have no conflicts of interest.

Supplementary data

Supplementary data is available at Horticulture Research online.

References

[1] Silva e Souza R, PAM B, Yassue RM. et al. Combining ability for the improvement of vegetable soybean. Agron J. 2020; 112:3535-48

[2] Mohamed Ali I, Rangappa M. Screening soybean (grain and vegetable) genotypes for nutrients and anti-nutritional factors. Plant Foods Hum Nutr. 1992; 42:87-96

[3] Liu CK, Li YS, Tu BJ. et al. Seed nutritional quality comparison of vegetable soybean genotypes at fresh pod and mature stage. Emir J Food Agric. 2019; 31:405-14

[4] Sonkamble P, Nandanwar RS, Sakhare SB. et al. Assessment of grain and vegetable type soybean genotypes for nutritional composition. Int J Chem Stud. 2020; 8:619-21

[5] Jiang GL, Townsend W, Ren S. Analysis of seed amino acids in vegetable soybeans dried by freeze and thermal drying. Agron-omy. 2023; 13:574

[6] Liu N, Niu Y, Zhang G. et al. Genome sequencing and popu-lation resequencing provide insights into the genetic basis of domestication and diversity of vegetable soybean. Hortic Res. 2022;9:uhab052

[7] Saldivar X, Wang YJ, Chen P. et al. Changes in chemical compo-sition during soybean seed development. Food Chem. 2011; 124: 1369-75

[8] Zhang QY, Li YS, Liu CK. et al. Key components of eating quality and their dynamic accumulation in vegetable soybean varieties [Glycine max (L.) Merr.]. Acta Agron Sin. 2015;41:1692-700 (in Chinese with English abstract)

[9] Dijk A, Kootstra G, Kruijer W. et al. Machine learning in plant science and plant breeding. iScience. 2020; 24:101890

[10] Wang M, Tai C, Weinan E. et al. DeFine: deep convolutional neural networks accurately quantify intensities of transcription factor-DNA binding and facilitate evaluation of functional non-coding variants. Nucleic Acids Res. 2018; 46:e69

[11] Zhu WC, Li WF, Zhang HW. et al. Big data and artificial intelligence-aided crop breeding: progress and prospects. J Integr Plant Biol. 2024; 00:1-18

[12] Çetin N. Machine learning for varietal binary classification of soybean (Glycine max (L.) Merrill) seeds based on shape and size attributes. Food Anal Methods. 2022; 15:2260-73

[13] Teodoro LPR, Silva MO, dos Santos RG. et al. Machine learning for classification of soybean populations for industrial tech-nological variables based on agronomic traits. Euphytica. 2024; 220:40

[14] Wang H, Cimen E, Singh N. et al. Deep learning for plant genomics and crop improvement. Curr Opin Plant Biol. 2020; 54: 34-41

[15] Fong JHC, Wong ASL. Advancing CRISPR/Cas gene editing with machine learning. Curr Opin. Biomed Eng. 2023; 28:100477

[16] Yu X, Yuan F, Fu X. et al. Profiling and relationship of water-soluble sugar and protein compositions in soybean seeds. Food Chem. 2016; 196:776-82

[17] Farrant JM, Hilhorst H. Crops for dry environments. Curr Opin Biotechnol. 2022; 74:84-91

[18] Manan S, Ahmad MZ, Zhang G. et al. Soybean LEC2 regulates subsets of genes involved in controlling the biosynthesis and catabolism of seed storage substances and seed development. Front Plant Sci. 2017; 8:1604

[19] LiuC, ChenH, YuQ. et al. Identification of quantitative trait loci and candidate genes for seed sucrose and soluble sugar concentrations in soybean. Crop Sci. 2023; 63:2976-92

[20] Wang S, Yokosho K, Guo R. et al. The soybean sugar transporter GmSWEET15 mediates sucrose export from endosperm to early embryo. Plant Physiol. 2019; 180:2133-41

[21] Li W, Huang L, Liu N. et al. Key regulators of sucrose metabolism identified through comprehensive comparative transcriptome analysis in peanuts. Int J Mol Sci. 2021; 22:7266

[22] Liu CK, Wang X, Chen H. et al. Nutritional quality of different potassium efficiency types of vegetable soybean as affected by potassium nutrition. Food Qual Saf. 2022; 6:1-11

[23] Morley SA, Ma F, Alazem M. et al. Expression of malic enzyme reveals subcellular carbon partitioning for storage reserve pro-duction in soybeans. New Phytol. 2023; 239:1834-51

[24] Wang R, Peng S, Chi Z. et al. Integrative analyses of metabolome and genome-wide transcriptome reveal the regulatory network governing flavor formation in kiwifruit (Actinidia chinensis). New Phytol. 2022; 233:373-89

[25] Machado FB, Moharana KC, Almeida-Silva F. et al. System-atic analysis of 1298 RNA-Seq samples and construction of a comprehensive soybean (Glycine max) expression atlas. Plant J. 2020; 103:1894-909

[26] Hu Y, Liu Y, Lu L. et al. Global analysis of seed transcriptomes reveals a novel platz regulator for seed size and weight control in soybean. New Phytol. 2023; 240:2436-54

[27] Bilska-Kos A, Mytych J, Suski S. et al. Sucrose phosphate synthase (SPS), sucrose synthase (SUS) and their products in the leaves of Miscanthus × giganteus and Zea Mays at low temperature. Planta. 2020; 252:23

[28] Agyenim-Boateng KG, Zhang S, Zhang S. et al. The nutritional composition of the vegetable soybean (maodou) and its poten-tial in combatting malnutrition. Front Nutr. 2023; 9:1034115

[29] Song J, Liu C, Dajing L. et al. Evaluation of sugar, free amino acid, and organic acid compositions of different varieties of vegetable soybean (Glycine max [L.] Merr). Ind Crop Prod. 2013; 50: 743-9

[30] Moseley D, Mozzoni L, Kaler A. et al. Evaluation of genetic diversity and association mapping for seed weight and size in vegetable soybean germplasm. Crop Sci. 2021; 61:3516-28

[31] Zhang Q, Li Y, Chin KL. et al. Vegetable soybean: seed composi-tion and production research. Ital J Agron. 2017; 12:872

[32] Wszelaki AL, Delwiche JF, Walker SD. et al. Consumer liking and descriptive analysis of six varieties of organically grown edamame-type soybean. Food Qual Prefer. 2005; 16:651-8

[33] TuB, LiuC, TianB. et al. Reduced abscisic acid content is responsible for enhanced sucrose accumulation by potassium nutrition in vegetable soybean seeds. JPlant Res. 2017; 130: 551-8

[34] Kumar V, Rani A, Goyal L. et al. Sucrose and raffinose fam-ily oligosaccharides (RFOs) in soybean seeds as influenced by genotype and growing location. J Agric Food Chem. 2010; 58: 5081-5

[35] Liu CK, Tu BJ, Li YS. et al. Potassium application affects key enzyme activities of sucrose metabolism during seed filling in vegetable soybean. Crop Sci. 2017; 57:2707-17

[36] Singh SK, Barnaby JY, Reddy VR. et al. Varying response of the concentration and yield of soybean seed mineral elements, carbohydrates, organic acids, amino acids, protein, and oil to phosphorus starvation and CO2 enrichment. Front Plant Sci. 2016; 7:1967

[37] Zhang J, Wang X, Lu Y. et al. Genome-wide scan for seed compo-sition provides insights into soybean quality improvement and the impacts of domestication and breeding. Mol Plant. 2018; 11: 460-72

[38] Wang S, Pan C, Zhang C. et al. Key enzymes of sucrose metabolism play an important role in source sink regulation of walnut fruit growth and development. NZJCropHort. 2023; 51: 400-19

[39] Kambhampati S, Aznar-Moreno JA, Bailey SR. et al. Temporal changes in metabolism late in seed development affect biomass composition. Plant Physiol. 2021; 186:874-90

[40] Duan Z, Zhang M, Zhang Z. et al. Natural allelic variation of GmST05 controlling seed size and quality in soybean. Plant Biotechnol J. 2022; 20:1807-18

[41] Hostmark AT, Haug A. Percentage oleic acid is inversely related to percentage arachidonic acid in total lipids of rat serum. Lipids Health Dis. 2013; 12:4

[42] Dhakal KH, Jung KH, Chae JH. et al. Variation of unsaturated fatty acids in soybean sprout of high oleic acid accessions. Food Chem. 2014; 164:70-3

[43] Nan H, Lu S, Fang C. et al. Molecular breeding of a high oleic acid soybean line by integrating natural variations. Mol Breed. 2020; 40:87

[44] Fu M, Chen L, Cai Y. et al. CRISPR/Cas9-mediated mutagenesis of GmFAD2-1A and/or GmFAD2-1 B to create high-oleic-acid soy-bean. Agronomy. 2022; 12:3218

[45] Li Y, Du M, Zhang Q. et al. Greater differences exist in seed protein, oil, total soluble sugar and sucrose content of vegetable soybean genotypes (Glycine max (L.) Merrill) in Northeast China. Aust J Crop Sci. 2012; 6:1681-6

[46] Liu N, Feng W, Zhang G. et al. Genomic and population evidence uncovers divergent improvement of vegetable soybean from grain soybean. Mol Plant. 2025; 18:1094-7

[47] Mariotti M, Lucisano M. Sugar and sweeteners. In: Zhou W, Hui YH, De Leyn I, Pagani MA, Rosell CM, Selman JD, Therdthai N, eds, Bakery Products Science and Technology. 2nd ed. Hoboken, NJ: John Wiley & Sons, Ltd; 2014,199-221

[48] He CY. Identification of Soybean SPS Gene Family and Preliminary Functional Verification of GmSPSA-2 in Roots. Harbin, China: North-east Agricultural University. 2023: (in Chinese with English abstract)

[49] Shen J, Xu Y, Yuan S. et al. Genome-wide identification of GmSPS gene family in soybean and expression analysis in response to cold stress. Int J Mol Sci. 2023; 24:12878

[50] Volkert K, Debast S, Voll LM. et al. Loss of the two major leaf isoforms of sucrose-phosphate synthase in Arabidopsis thaliana limits sucrose synthesis and nocturnal starch degradation but does not alter carbon partitioning during photosynthesis. JExp Bot. 2014; 65:5217-29

[51] Okamura M, Aoki N, Hirose T. et al. Tissue specificity and diurnal change in gene expression of the sucrose phos-phate synthase gene family in rice. Plant Sci. 2011; 181: 159-66

[52] Esparza-Reynoso S, Ruíz-Herrera LF, Pelagio-Flores R. et al. Tri-choderma atroviride-emitted volatiles improve growth of Ara-bidopsis seedlings through modulation of sucrose transport and metabolism. Plant Cell Environ. 2021; 44:1961-76

[53] Mitsuda N, Ikeda M, Takada S. et al. Efficient yeast one-/two-hybrid screening using a library composed only of transcrip-tion factors in Arabidopsis thaliana. Plant Cell Physiol. 2010; 51: 2145-51

[54] Bhattacharjee A, Ghangal R, Garg R. et al. Genome-wide anal-ysis of homeobox gene family in legumes: identification, gene duplication and expression profiling. PLoS One. 2015; 10: e0119198

[55] Fehr WR, Caviness CE. Stages of Soybean Development. Special Report. Ames, IA, USA: Iowa State University of Science and Technology; 1977:87

[56] Qin P, Song W, Yang X. et al. Regional distribution of protein and oil compositions of soybean cultivars in China. Crop Sci. 2014; 54: 1139-46

[57] RenY, YuG, Shi C. et al. Majorbio cloud: a one-stop, compre-hensive bioinformatic platform for multiomics analyses. iMeta. 2022; 1:e12

[58] Lescot M, Déhais P, Thijs G. et al. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002; 30:325-7

[59] Matthews BF, Youssef RM. Agrobacterium rhizogenes-based trans-formation of soybean roots to form composite plants. Bio-protocol. 2016; 6:e1708

[60] Archer E.EricArcher/rfPermute: version 2.5. Geneva, Switzerland: Zenodo; 2021