Comparative Study of the Nutritional, Phytochemical Composition and Antioxidant Potential of Soybean [Glycine max (L.) Merrill] from Burkina Faso

P
Pierre Alexandre Eric Djifaby Sombié3
S
Samson Guenne4
H
Hemayoro Sama4,5
I
Ibié Gilles Thio2
I
Issaka Ouedraogo4
O
Oumar Boro6
N
Nerbéwendé Sawadogo1
1Biosciences Laboratory, Joseph KI-ZERBO University, Ouagadougou, Burkina Faso.
2Institute of Environment and Agricultural Research (INERA), Ouagadougou, Burkina Faso.
3Institute of Environment and Agricultural Research/Farako-Ba (INERA/Farako-Ba), Bobo-Dioulasso, Burkina Faso.
4Laboratory of Biochemistry and Chemistry Applied (LABIOCA), Joseph KI-ZERBO University, Ouagadougou, Burkina Faso.
5Laboratory of Biochemistry, Biotechnology, Food Technology and Nutrition, Joseph KI-ZERBO University, Ouagadougou, Burkina Faso.
6West Africa Centre for Crop Improvement, School of Agriculture, University of Ghana, Legon, Ghana.
  • Submitted24-12-2025|

  • Accepted25-02-2026|

  • First Online 23-03-2026|

  • doi 10.18805/LRF-927

Background: The soybean crop plays a crucial role in Burkina Faso’s agricultural sector and contributes significantly to food security and nutrition. The study aimed to evaluate the nutritional composition and antioxidant potential of soybean cultivated in Burkina Faso.

Methods: A total of 40 soybean varieties cultivated in 2022 at two sites in Burkina Faso were analyzed using standard methods for six parameters, including total phenolic, total flavonoid, total protein, soluble sugars, phytic acid, DPPH radical scavenging activity and ferric reducing antioxidant power (FRAP). Total phenolic ranged from 1.769 to 2.923 µg GAE/100 mg, total flavonoid from 0.008 to 0.031 µg QE/100 mg, total protein from 1.526 to 2.021 mg/100 mg, soluble sugar from 107.6 to 228.1 µg GE/100 mg and phytic acid concentration from 130.9 to 378.0 µg APE/100 mg. Antioxidant activity varied from 0.786 to 2.282 µg AAE/100 mg for DPPH and from 0.019 to 0.079 µg AAE/100 mg dw for FRAP.

Result: Four soybean varieties, GPV12, GPV17, GPV22 and GMV19 are characterized by high nutrient content, particularly total protein and soluble sugars. The GPV3, GPV6 and GPV8 lines are distinguished by their high concentration of bioactive compounds, such as total phenolic compounds and total flavonoids. In contrast, the GPV4 and GPV5 varieties have relatively high levels of antinutrients. The varieties with high nutritional potential could help strengthen the resilience of poor households against food insecurity in Sub-Saharan Africa.

Soybean [Glycine max (L.) Merrill] is a legume widely consumed worldwide due to its high nutritional value and health benefits (Amol et al., 2021). In terms of protein, the nutritional value of soybean amino acids is significantly higher than that of most legumes (Messina, 2022). The seed is an excellent source of essential nutrients, providing proteins, essential amino acids, oil, carbohydrates, dietary fiber, vitamins and minerals for both human and animal nutrition (Alghafis and Raouf, 2020). Moreover, soybean, as a nitrogen-fixing legume, helps reduce the use of chemical sources for nitrogen fertilizer production (Patil et al., 2024). However, the relatively low genetic variation among soybean cultivars limits the improvement of these nutritional characteristics. Determining the phytochemical diversity of soybean cultivars is of utmost importance (Azam et al., 2021). This will help breeders to develop new cultivars with enhanced nutrition and expand the genetic diversity of soybean cultivars (Singer et al., 2023). The nutritional and phytochemical characteristics of soybean are mainly influenced by various factors such as genotype, location, climate, water availability and maturity group (Jiang et al., 2018). Multi-location trials are necessary for the effective selection of different genotypes with higher levels of desirable nutritional traits (Azam et al., 2021). The present study was conducted to identify the variability of total phenolic, total flavonoid, total protein, soluble sugars, phytic acid, 2,2-diphenyl-1-picrylhydrazyl (DPPH) radical scavenging activity and ferric reducing antioxidant power (FRAP) in Burkina Faso’ soybean genotype.
Plant material
 
This research was conducted in 2024 at the Laboratory of Biochemistry and Chemistry Applied (LaBioCa), University Joseph Ki-Zerbo, Burkina Faso. The study material (40) consisted of 37 soybean lines (Table 1) developed by the International Institute of Tropical Agriculture (IITA) and 3 soybean lines from the Institute of Environment and Agricultural Research (INERA). The seeds used were harvested from plants at full maturity in two different agroecological zones of Burkina Faso: the northern Sudanian zone (Gampela: 12°15' North latitude, 1°12' West longitude) and the southern Sudanian zone (Farako-Ba: 11°12'0" North latitude, 4°18'0" West longitude).

Table 1: List of soybean lines used in this study.


 
Determination of nutrients and antinutrients contents
 
Total protein content
 
Total protein content was evaluated at 595 nm by using the methods of Bradford (1976). The means of three values were obtained, expressed as mg/100 mg dw.
 
Soluble sugar content
 
Total soluble sugars was evaluated using the phenol-acid sulfuric method of  DuBois et al., (1956) and the absorbance was read at 490 nm. The total sugar content was expressed as µg GE/100 mg dw.
 
Phytic acid content
 
Phytic acid content (PAC) was determined using the method of Latta and Eskin (1980). The PAC was determined using phytic acid as the standard and expressed in µg APE/100 mg dw.
 
Bioactive compounds quantification
 
Extraction method of bioactive compounds
 
The dried soybean seeds were ground. Seeds powders (0.5 g each) from different soybean varieties were extracted with 10 mL of ethanol 80% (v/v). The mixture was shaken for 24 h and centrifuged at 4500 rpm for 10 minutes. The supernatant was used for the quantification of total phenolics, total flavonoids, DPPH and FRAP.
 
Total phenolic content
 
Total phenolic content of the soybean seed extracts was determined at 760 nm according to the method described by Singleton et al., (1999), using gallic acid as a standard. Total phenolic was expressed in µg of gallic acid equivalent per 100 g of dry weight.
 
Total flavonoid content
 
The total flavonoid content of the soybean seed extracts was determined at 415 nm using the method described by Arvouet-Grand  et al. (1994). The total flavonoid content was determined on a quercetin calibration curve and expressed as µg of quercetin equivalents (QE) per 100 mg of dry seed weight.
 
DPPH radical scavenging activity
 
The ability of soybean seed extracts to scavenge the DPPH radical was evaluated at 517 nm as described by Sombié  et al. (2011). The means of three values were obtained, expressed as µg of ascorbic acid equivalent per 100 g of dry seeds weight (µg AAE/100 mg dw).
 
Ferric reducing antioxidant power (FRAP) assay
 
The ability of the cowpea seed extracts to reduce iron (III) to iron (II) was measured at 700 nm following the procedure described by Hinneburg et al., (2006). Iron (III) reducing activity was determined as µg of ascorbic acid equivalent (AAE) per 100 mg of dry seed weight (µg AAE)/100 mg dw).
 
Statistical analysis
 
Before analyzing variance (ANOVA), the assumptions of normality and homoscedasticity were verified using the Shapiro-Wilk and Levene’s tests, respectively. An analysis of variance was performed to estimate the effects of genotype, environment and genotype x environment (GxE) interaction on the measured traits at the 1% threshold level. Analysis of variance across the two sites was conducted for each trait and the LSD (least significant difference) test was applied to identify groups of genotypes that actually differ from each other. Pearson’s correlation test was used to determine the relation between phenolic content, flavonoid content and antioxidant activities using the R software package (version 4.4.0).
Variation in nutrient and antinutrient contents, bioactive compounds and antioxidant activity
 
Soybean is an important but underutilized legume contributing to food and nutritional security in Burkina Faso. This study evaluated variation in nutrient, antinutrient, phytochemical composition and antioxidant activity among 40 soybean genotypes grown under two contrasting environments. Analysis of variance revealed highly significant effects (p≤0.001) of genotype (G), environment (E) and genotype x environment (G x E) interaction for all traits (Table 2). Genotype explained the largest proportion of total variation for most traits, notably total phenolic content (57.59%), soluble sugars (63.17%), FRAP (58.35%) and phytic acid content (75.04%), consistent with previous reports (Njoroge and Oyoo, 2020; Wan, 2023). In contrast, total flavonoid content, total protein and DPPH activity were more strongly influenced by G×E interaction, indicating differential genotype responses across environments, as also reported by Hailemariam (2022).

Table 2: Combined analysis of variance for 7 traits and the percentage of mean squares for genotype (G), environment (E) and genotype x environment interaction (G x E).


       
Significant genetic variability was observed among genotypes for all studied traits (Table 3). Total flavonoid content ranged from 0.008 to 0.031 µg QE/100 mg, while total phenolic content exceeded the grand mean (2.233 µg GAE/100 mg) in 55% of the genotypes. Genotypes GPV3 and GPV8 exhibited high phenolic content and strong antioxidant activity (DPPH), confirming the close association between phenolics and antioxidant capacity. Similar environment-dependent variability in phytochemical composition has been attributed to biochemical regulation of antioxidant enzymes and phenolic biosynthesis in response to growing conditions (Kokebie et al., 2024).

Table 3: Mean values of soybean phytochemical traits of 40 genotypes evaluated at gampela and farako-Ba, burkina faso, during the 2022 growing season.


 
Performance of the genotypes in each environment
 
Genotype performance varied markedly between environments (Fig 1). Farako-Ba, characterized by higher rainfall, generally promoted greater accumulation of nutritional and antioxidant compounds compared with the drier Gampela site, in agreement with earlier findings on the influence of climatic conditions on phenolic content and antioxidant capacity (Tolic et al., 2017). Several genotypes exhibited high levels of phytic acid and soluble sugars, traits of functional and nutritional interest. Although phytic acid is considered an antinutrient due to mineral chelation (Lott et al., 2000), its antioxidant and anticancer properties have been widely documented (Vucenik and Shamsuddin, 2006). Similarly, soluble sugars play a key role in glycemic regulation and intestinal health (Fu et al., 2022), with marked genetic variability previously reported in tropical soybean environments (Jiang et al., (2018). Overall, the wide genetic diversity observed for nutritional, phytochemical and antioxidant traits highlights the importance of genotype selection and environment-specific breeding strategies for improving soybean nutritional quality and sustainability in sub-Saharan Africa.

Fig 1: Biochemical compound content of soybean varieties grown at gampela (E1) and farako-Ba (E2), burkina faso, during the 2022 growing season.


 
Relationship between antioxidant activity and the studied traits
 
To identify the phytochemical compounds likely contributing to the antioxidant activity of soybeans, Pearson’s correlation coefficient was calculated (Table 4). Total phenolic content is positively correlated with DPPH (r = 0.519) and FRAP (r = 0.304). The antioxidant activities show a strong positive correlation with each other (r = 0.508). However, no significant relationship was found between the antioxidant activity and total flavonoids content. Similar results were reported by Priastomo et al., (2024). The total protein content is negatively correlated with FRAP (r = -0.237). Similarly, the phytate content is negatively and weakly correlated with FRAP (r = -0.164). The positive correlation between DPPH, FRAP and total phenolic content in soybean seeds indicates that higher phenolic levels enhance antioxidant activity. This finding is consistent with a previous study that reported that the phenolic compounds are responsible for the level of antioxidant activity (Arifin et al., 2021). Studies show that as total phenolic content increases, both DPPH scavenging and FRAP values improve, reflecting a synergistic relationship where phenolics contribute significantly to the overall antioxidant capacity of soybean seeds (Alsuwayt, 2025). Flavonoids present in soybeans play a crucial role in improving protein quality (Oliveira et al., 2024). They can enhance protein digestibility and nutritional value by acting as antioxidants, which protect amino acids from oxidative damage (Prodić et al., 2023). By incorporating flavonoids into the diet, it is possible to potentially improve the overall quality of soybean proteins (Jia et al., 2022).

Table 4: Pearson’s correlation coefficient.


 
Association between nutritional and phytochemical contents
 
The first three principal components explained 69.60% of the total variance (Table 5; Fig 2), indicating a substantial representation of the biochemical diversity among soybean genotypes. The first principal component (Dim.1; eigenvalue = 2.44) accounted for 34.85% of the total variance and was strongly associated with antioxidant activities (DPPH, FRAP) and total phenolic content, reflecting the bioactive potential of the genotypes. Genotypes positioned in the positive quadrant of Dim.1 exhibited high antioxidant capacity and phenolic content and are therefore of interest for nutraceutical and functional applications. In this regard, GMV21, GPV20, GPV3, GPV8 and GMV20 displayed the greatest variability and represent promising candidates. The second principal component (Dim.2; eigenvalue = 1.30), explaining 18.62% of the total variance, was mainly influenced by total flavonoid and protein contents, indicating a nutritional and functional dimension. Genotypes located in the positive quadrant of Dim.2 combined high flavonoid and protein levels, making them suitable for breeding programs targeting both nutritional quality and health benefits. The most variable genotypes along this axis were GPV5, GPV12, GMV18, GMV16 and GMV9. The third principal component (Dim.3; eigenvalue = 1.12) explained an additional 16.12% of the total variance and was characterized by strong contributions from soluble sugars and phytic acid, highlighting variability related to carbohydrate content and specific antinutritional compounds. Genotypes GMV13, GMV12, GMV7, GMV11 and GPV1 showed the highest variability along this component. While genotypes located in the negative quadrants of Dim.1 and/or Dim.2 exhibited lower antioxidant or functional compound contents, some displayed high soluble sugar or tannin levels, indicating alternative valorization pathways.

Table 5: Eigenvalues and proportion of the variance explained for the three principal components of the 40 soybean genotypes based on phytochemical components.



Fig 2: PCA showing the spatial distribution of 40 soybean genotypes based on their Nutritional, phytochemical and antioxidant characteristics.


       
Globally, several genotypes (GMV5, GMV15, GMV21 and GPV8) showed positive loadings in at least two principal components, demonstrating their potential for targeted phytochemical improvement. The wide dispersion of genotypes confirms high biochemical diversity within the collection, allowing effective discrimination based on antioxidant, nutritional and phytochemical attributes and providing a robust basis for selection in breeding and utilization programs.
 
Cluster and heatmap analysis
 
A combined cluster and heatmap analysis were conducted to assess the variability among soybean genotypes concerning the total flavonoid content, total phenolic content, total protein, soluble sugars, phytic acid content and antioxidant activity. All genotypes were divided into three main groups (Table 6 and 7). In particular, Fig 3 shows the cluster analysis according to the average linkage method. Cluster 1 included only three varieties (GPV5, GPV4 and GPV7), which have high values in phytic acid content and low content in bioactive compounds (DPPH, FRAP, TPC, TFC). Cluster 2 also contained only two varieties (GPV6, GMV7) that have the highest total flavonoid content and total proteins. Cluster 2 is intermediate in antioxidant compounds and total flavonoids. Cluster 3 consisted of three varieties (GPV3, GPV8 and GPV20), which are prominent in the total phenol content and antioxidant activity. This cluster is rich in flavonoids and polyphenols and has high antioxidant activity (DPPH, FRAP), but low phytic acid content.

Table 6: Average cluster values of nutritional components, phytochemicals and antioxidant activities of 40 soybean genotypes.



Table 7: Composition and characteristics of the different groups.



Fig 3: Cluster and heatmap analysis based on significantly correlated variables.

In this study, for all the traits analyzed, genetic variation was primarily influenced by genotype, followed by the effects of the environment and the genotype × environment interaction (GEI). The Gampela locality was found to be less productive in terms of nutritional compounds, phytochemicals and antioxidant activities, whereas Farako-Ba was more productive. The variability among soybean varieties was also described through cluster analysis, which divided them into three groups. Moreover, correlation analysis indicated that total phenolic content was correlated with antioxidant activities. These results could contribute to the improvement of soybean breeding programs through an approach based on nutritional, phytochemical and antioxidant criteria. The soybean lines GPV12, GPV17, GPV22 and GMV19 are characterized by high nutrient content, particularly total protein and soluble sugars. The GPV3, GPV6 and GPV8 lines are distinguished by their high concentration of bioactive compounds, such as total phenolic compounds and total flavonoids. In contrast, the GPV4 and GPV5 lines have relatively high levels of antinutrients. Furthermore, this study makes a significant contribution to the mapping of Burkina Faso’s plant genetic resources.
 
The present study was supported by the International Center of Genetic Engineering and Biotechnology.
The authors have declared that no conflict of interest exists.

  1. Alghafis, A. and Raouf, E.A. (2020). Optimization of injection timing and injection duration of a diesel engine running on pure biodiesel SME (Soya methyl ester). Open J. Appl. Sci. 10: 486.

  2. Alsuwayt, B. (2025). Evaluation of antioxidant, anti-inflammatory and analgesic potentials of combined polyphenol-rich fractions from Ziziphus mauritiana and Ziziphus spina- christi leaves through modulation of inflammatory and oxidative stress markers in Sprague Dawley rats model. Inflammopharmacology. 33(9): 5509-5523. https:// doi.org/10.1007/s10787-025-01686-1.

  3. Amol, V., Bhati, K.R. and Bhati, K.R. (2021). Nutritive benefits of soybean (Glycine max). Indian J. Nutr. Diet. pp 522- 533. https://doi.org/10.21048/IJND.2021.58.4.27339.

  4. Arifin, H.A., Hashiguchi, T., Nagahama, K., Hashiguchi, M., Muguerza, M., Sakakibara, Y., Tanaka, H. and Akashi, R. (2021). Varietal differences in flavonoid and antioxidant activity in Japanese soybean accessions. Biosci. Biotechnol. Biochem. 85: 916-922.

  5. Arvouet-Grand, A., Vennat, B., Pourrat, A. and Legret, P. (1994). Standardization of propolis extract and identification of principal constituents. J. Pharm. Belg. 49: 462-468.

  6. Azam, M., Zhang, S., Qi, J., Abdelghany, A.M., Shaibu, A.S., Ghosh, S., Feng, Y., Huai, Y., Gebregziabher, B.S., Li, J., Li, B. and Sun, J. (2021). Profiling and associations of seed nutritional characteristics in Chinese and USA soybean cultivars. J. Food Compos. Anal. 98: 103803. https:// doi.org/10.1016/j.jfca.2021.103803.

  7. Bradford, M.M. (1976). A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal. Biochem. 72: 248-254.

  8. DuBois, M., Gilles, K.A., Hamilton, J.K., Rebers, P.T. and Smith, F. (1956). Colorimetric method for determination of sugars and related substances. Anal. Chem. 28: 350-356.

  9. Fu, J., Zheng, Y., Gao, Y. and Xu, W. (2022). Dietary fiber intake and gut microbiota in human health. Microorganisms. 10: 2507. https://doi.org/10.3390/microorganisms10122507.

  10. Hailemariam, H.M. (2022). Adaptability and stability for soybean yield by AMMI and GGE models in Ethiopia. Front. Plant Sci. 13: 950992.

  11. Hinneburg, I., Dorman, H.D. and Hiltunen, R. (2006). Antioxidant activities of extracts from selected culinary herbs and spices. Food Chem. 97: 122-129.

  12. Hungria, M. and Mendes, I.C. (2015). Nitrogen fixation with soybean: The perfect symbiosis? Biol. Nitrogen Fixat. pp 1009- 1024.

  13. Jia, Y., Fu, Y., Man, H., Yan, X., Huang, Y., Sun, S., Qi, B. and Li, Y. (2022). Comparative study of binding interactions between different dietary flavonoids and soybean â- conglycinin and glycinin: Impact on structure and function of the proteins. Food Res. Int. 161: 111784. https:// doi.org/10.1016/j.foodres.2022.111784.

  14. Jiang, G.L., Chen, P., Zhang, J., Florez-Palacios, L., Zeng, A., Wang, X., Bowen, R.A., Miller, A. and Berry, H. (2018). Genetic analysis of sugar composition and its relationship with protein, oil and fiber in soybean. Crop Sci. 58: 2413-2421.

  15. Kokebie, D., Enyew, A., Masresha, G., Fentie, T. and Mulat, E. (2024). Morphological, physiological and biochemical responses of three different soybean (Glycine max L.) varieties under salinity stress conditions. Front. Plant Sci. 15. https://doi.org/10.3389/fpls.2024.1440445.

  16. Latta, M. and Eskin, M. (1980). A simple and rapid colorimetric method for phytate determination. J. Agric. Food Chem. 28: 1313-1315.

  17. Lott, J.N., Ockenden, I., Raboy, V. and Batten, G.D. (2000). Phytic acid and phosphorus in crop seeds and fruits: A global estimate. Seed Sci. Res. 10: 11-33.

  18. Messina, M. (2022). Perspective: Soybeans can help address the caloric and protein needs of a growing global population. Front. Nutr. 9: 909464. https://doi.org/10.3389/fnut. 2022.909464.

  19. Njoroge, J.N. and Oyoo, M.E. (2020). Stability and adaptability of soybean [Glycine max (L.) Merrill] genotypes for yield, protein and oil content using AMMI analysis in Kenya. Agric. Sci. Dig. 40(2): 122-128. doi: 10.18805/ag.D-218.

  20. Oliveira, I.C., Santana, D.C., de Oliveira, J.L.G., Silva, E.V.M., da Silva, C.S.A.C., Blanco, M., Teodoro, L.P.R., da Silva Júnior, C.A., Baio, F.H.R., Alves, C.Z. and Teodoro, P.E. (2024). Flavonoids and their relationship with the physiological quality of seeds from different soybean genotypes. Sci. Rep. 14: 17008. https://doi.org/10.1038/ s41598-024-68117-z.

  21. Patil, S.S., Hiwale, S.D., Patil, P.P., Yadav, K.R., Patil, S.R. and Bisarya, D. (2024). Impact of levels of nutrient management on economics of soybean-onion and soybean-potato cropping system. Indian J. Agric. Res. 59(5): 725-730. doi: 10.18805/IJARe.A-6144.

  22. Priastomo, M., Adlia, A.R., Lumbantobing, V. and Adnyana, I.K. (2024). Determination of total phenols, total flavonoids and antioxidant activity of watermelon peel and rind from several cultivation areas in Indonesia. Indian J. Agric. Res. 58(5): 865-871. doi: 10.18805/IJARe.AF-872.

  23. Prodic, I., Krstic, R.M. and Smiljanic, K. (2023). Antioxidant properties of protein-rich plant foods in gastrointestinal digestion- peanuts as our antioxidant friend or foe in allergies. Antioxidants. 12: 886. https://doi.org/10.3390/antiox12040886.

  24. Singer, W.M., Lee, Y., Shea, Z., Vieira, C.C., Lee, D., Li, X., Cunicelli, M., Kadam, S.S., Khan, M.A.W., Shannon, G., Mian, M.A.R., Nguyen, H.T. and Zhang, B. (2023). Soybean genetics, genomics and breeding for improving nutritional value and reducing antinutritional traits in food and feed. Plant Genome. 16: e20415. https://doi.org/10.1002/tpg2.20415.

  25. Singleton, V.L., Orthofer, R. and Lamuela-Raventós, R.M. (1999). Analysis of total phenols and other oxidation substrates and antioxidants by means of folin-ciocalteu reagent, In: Methods in Enzymology. Elsevier. pp. 152-178.

  26. Sombié, P., Hilou, A., Mounier, C., Coulibaly, A., Kiendrebeogo, M., Millogo, J. and Nacoulma, O. (2011). Antioxidant and anti-inflammatory activities from galls of Guiera senegalensis JF Gmel (Combretaceae). Research Journal of Medicinal Plants. 5(4): 448-4461.

  27. Tolic, M.T., Krbavcic, I.P., Vujevic, P., Milinovic, B., Jurcevic, I.L. and Vahcic, N. (2017). Effects of weather conditions on phenolic content and antioxidant capacity in juice of chokeberries (Aronia melanocarpa L.). Pol. J. Food Nutr. Sci. 67(1): 67-74.

  28. Vucenik, I. and Shamsuddin, A.M. (2006). Protection against cancer by dietary IP6 and inositol. Nutr. Cancer. 55: 109-125.

  29. Wan, Z. (2023). Exploring the effects of genotype, environment and genotype by environment interactions on field pea (Pisum sativum L.) protein and amino acid contents using near-infrared reflectance spectroscopy. Food and Human Nutritional Sciences.

Comparative Study of the Nutritional, Phytochemical Composition and Antioxidant Potential of Soybean [Glycine max (L.) Merrill] from Burkina Faso

P
Pierre Alexandre Eric Djifaby Sombié3
S
Samson Guenne4
H
Hemayoro Sama4,5
I
Ibié Gilles Thio2
I
Issaka Ouedraogo4
O
Oumar Boro6
N
Nerbéwendé Sawadogo1
1Biosciences Laboratory, Joseph KI-ZERBO University, Ouagadougou, Burkina Faso.
2Institute of Environment and Agricultural Research (INERA), Ouagadougou, Burkina Faso.
3Institute of Environment and Agricultural Research/Farako-Ba (INERA/Farako-Ba), Bobo-Dioulasso, Burkina Faso.
4Laboratory of Biochemistry and Chemistry Applied (LABIOCA), Joseph KI-ZERBO University, Ouagadougou, Burkina Faso.
5Laboratory of Biochemistry, Biotechnology, Food Technology and Nutrition, Joseph KI-ZERBO University, Ouagadougou, Burkina Faso.
6West Africa Centre for Crop Improvement, School of Agriculture, University of Ghana, Legon, Ghana.
  • Submitted24-12-2025|

  • Accepted25-02-2026|

  • First Online 23-03-2026|

  • doi 10.18805/LRF-927

Background: The soybean crop plays a crucial role in Burkina Faso’s agricultural sector and contributes significantly to food security and nutrition. The study aimed to evaluate the nutritional composition and antioxidant potential of soybean cultivated in Burkina Faso.

Methods: A total of 40 soybean varieties cultivated in 2022 at two sites in Burkina Faso were analyzed using standard methods for six parameters, including total phenolic, total flavonoid, total protein, soluble sugars, phytic acid, DPPH radical scavenging activity and ferric reducing antioxidant power (FRAP). Total phenolic ranged from 1.769 to 2.923 µg GAE/100 mg, total flavonoid from 0.008 to 0.031 µg QE/100 mg, total protein from 1.526 to 2.021 mg/100 mg, soluble sugar from 107.6 to 228.1 µg GE/100 mg and phytic acid concentration from 130.9 to 378.0 µg APE/100 mg. Antioxidant activity varied from 0.786 to 2.282 µg AAE/100 mg for DPPH and from 0.019 to 0.079 µg AAE/100 mg dw for FRAP.

Result: Four soybean varieties, GPV12, GPV17, GPV22 and GMV19 are characterized by high nutrient content, particularly total protein and soluble sugars. The GPV3, GPV6 and GPV8 lines are distinguished by their high concentration of bioactive compounds, such as total phenolic compounds and total flavonoids. In contrast, the GPV4 and GPV5 varieties have relatively high levels of antinutrients. The varieties with high nutritional potential could help strengthen the resilience of poor households against food insecurity in Sub-Saharan Africa.

Soybean [Glycine max (L.) Merrill] is a legume widely consumed worldwide due to its high nutritional value and health benefits (Amol et al., 2021). In terms of protein, the nutritional value of soybean amino acids is significantly higher than that of most legumes (Messina, 2022). The seed is an excellent source of essential nutrients, providing proteins, essential amino acids, oil, carbohydrates, dietary fiber, vitamins and minerals for both human and animal nutrition (Alghafis and Raouf, 2020). Moreover, soybean, as a nitrogen-fixing legume, helps reduce the use of chemical sources for nitrogen fertilizer production (Patil et al., 2024). However, the relatively low genetic variation among soybean cultivars limits the improvement of these nutritional characteristics. Determining the phytochemical diversity of soybean cultivars is of utmost importance (Azam et al., 2021). This will help breeders to develop new cultivars with enhanced nutrition and expand the genetic diversity of soybean cultivars (Singer et al., 2023). The nutritional and phytochemical characteristics of soybean are mainly influenced by various factors such as genotype, location, climate, water availability and maturity group (Jiang et al., 2018). Multi-location trials are necessary for the effective selection of different genotypes with higher levels of desirable nutritional traits (Azam et al., 2021). The present study was conducted to identify the variability of total phenolic, total flavonoid, total protein, soluble sugars, phytic acid, 2,2-diphenyl-1-picrylhydrazyl (DPPH) radical scavenging activity and ferric reducing antioxidant power (FRAP) in Burkina Faso’ soybean genotype.
Plant material
 
This research was conducted in 2024 at the Laboratory of Biochemistry and Chemistry Applied (LaBioCa), University Joseph Ki-Zerbo, Burkina Faso. The study material (40) consisted of 37 soybean lines (Table 1) developed by the International Institute of Tropical Agriculture (IITA) and 3 soybean lines from the Institute of Environment and Agricultural Research (INERA). The seeds used were harvested from plants at full maturity in two different agroecological zones of Burkina Faso: the northern Sudanian zone (Gampela: 12°15' North latitude, 1°12' West longitude) and the southern Sudanian zone (Farako-Ba: 11°12'0" North latitude, 4°18'0" West longitude).

Table 1: List of soybean lines used in this study.


 
Determination of nutrients and antinutrients contents
 
Total protein content
 
Total protein content was evaluated at 595 nm by using the methods of Bradford (1976). The means of three values were obtained, expressed as mg/100 mg dw.
 
Soluble sugar content
 
Total soluble sugars was evaluated using the phenol-acid sulfuric method of  DuBois et al., (1956) and the absorbance was read at 490 nm. The total sugar content was expressed as µg GE/100 mg dw.
 
Phytic acid content
 
Phytic acid content (PAC) was determined using the method of Latta and Eskin (1980). The PAC was determined using phytic acid as the standard and expressed in µg APE/100 mg dw.
 
Bioactive compounds quantification
 
Extraction method of bioactive compounds
 
The dried soybean seeds were ground. Seeds powders (0.5 g each) from different soybean varieties were extracted with 10 mL of ethanol 80% (v/v). The mixture was shaken for 24 h and centrifuged at 4500 rpm for 10 minutes. The supernatant was used for the quantification of total phenolics, total flavonoids, DPPH and FRAP.
 
Total phenolic content
 
Total phenolic content of the soybean seed extracts was determined at 760 nm according to the method described by Singleton et al., (1999), using gallic acid as a standard. Total phenolic was expressed in µg of gallic acid equivalent per 100 g of dry weight.
 
Total flavonoid content
 
The total flavonoid content of the soybean seed extracts was determined at 415 nm using the method described by Arvouet-Grand  et al. (1994). The total flavonoid content was determined on a quercetin calibration curve and expressed as µg of quercetin equivalents (QE) per 100 mg of dry seed weight.
 
DPPH radical scavenging activity
 
The ability of soybean seed extracts to scavenge the DPPH radical was evaluated at 517 nm as described by Sombié  et al. (2011). The means of three values were obtained, expressed as µg of ascorbic acid equivalent per 100 g of dry seeds weight (µg AAE/100 mg dw).
 
Ferric reducing antioxidant power (FRAP) assay
 
The ability of the cowpea seed extracts to reduce iron (III) to iron (II) was measured at 700 nm following the procedure described by Hinneburg et al., (2006). Iron (III) reducing activity was determined as µg of ascorbic acid equivalent (AAE) per 100 mg of dry seed weight (µg AAE)/100 mg dw).
 
Statistical analysis
 
Before analyzing variance (ANOVA), the assumptions of normality and homoscedasticity were verified using the Shapiro-Wilk and Levene’s tests, respectively. An analysis of variance was performed to estimate the effects of genotype, environment and genotype x environment (GxE) interaction on the measured traits at the 1% threshold level. Analysis of variance across the two sites was conducted for each trait and the LSD (least significant difference) test was applied to identify groups of genotypes that actually differ from each other. Pearson’s correlation test was used to determine the relation between phenolic content, flavonoid content and antioxidant activities using the R software package (version 4.4.0).
Variation in nutrient and antinutrient contents, bioactive compounds and antioxidant activity
 
Soybean is an important but underutilized legume contributing to food and nutritional security in Burkina Faso. This study evaluated variation in nutrient, antinutrient, phytochemical composition and antioxidant activity among 40 soybean genotypes grown under two contrasting environments. Analysis of variance revealed highly significant effects (p≤0.001) of genotype (G), environment (E) and genotype x environment (G x E) interaction for all traits (Table 2). Genotype explained the largest proportion of total variation for most traits, notably total phenolic content (57.59%), soluble sugars (63.17%), FRAP (58.35%) and phytic acid content (75.04%), consistent with previous reports (Njoroge and Oyoo, 2020; Wan, 2023). In contrast, total flavonoid content, total protein and DPPH activity were more strongly influenced by G×E interaction, indicating differential genotype responses across environments, as also reported by Hailemariam (2022).

Table 2: Combined analysis of variance for 7 traits and the percentage of mean squares for genotype (G), environment (E) and genotype x environment interaction (G x E).


       
Significant genetic variability was observed among genotypes for all studied traits (Table 3). Total flavonoid content ranged from 0.008 to 0.031 µg QE/100 mg, while total phenolic content exceeded the grand mean (2.233 µg GAE/100 mg) in 55% of the genotypes. Genotypes GPV3 and GPV8 exhibited high phenolic content and strong antioxidant activity (DPPH), confirming the close association between phenolics and antioxidant capacity. Similar environment-dependent variability in phytochemical composition has been attributed to biochemical regulation of antioxidant enzymes and phenolic biosynthesis in response to growing conditions (Kokebie et al., 2024).

Table 3: Mean values of soybean phytochemical traits of 40 genotypes evaluated at gampela and farako-Ba, burkina faso, during the 2022 growing season.


 
Performance of the genotypes in each environment
 
Genotype performance varied markedly between environments (Fig 1). Farako-Ba, characterized by higher rainfall, generally promoted greater accumulation of nutritional and antioxidant compounds compared with the drier Gampela site, in agreement with earlier findings on the influence of climatic conditions on phenolic content and antioxidant capacity (Tolic et al., 2017). Several genotypes exhibited high levels of phytic acid and soluble sugars, traits of functional and nutritional interest. Although phytic acid is considered an antinutrient due to mineral chelation (Lott et al., 2000), its antioxidant and anticancer properties have been widely documented (Vucenik and Shamsuddin, 2006). Similarly, soluble sugars play a key role in glycemic regulation and intestinal health (Fu et al., 2022), with marked genetic variability previously reported in tropical soybean environments (Jiang et al., (2018). Overall, the wide genetic diversity observed for nutritional, phytochemical and antioxidant traits highlights the importance of genotype selection and environment-specific breeding strategies for improving soybean nutritional quality and sustainability in sub-Saharan Africa.

Fig 1: Biochemical compound content of soybean varieties grown at gampela (E1) and farako-Ba (E2), burkina faso, during the 2022 growing season.


 
Relationship between antioxidant activity and the studied traits
 
To identify the phytochemical compounds likely contributing to the antioxidant activity of soybeans, Pearson’s correlation coefficient was calculated (Table 4). Total phenolic content is positively correlated with DPPH (r = 0.519) and FRAP (r = 0.304). The antioxidant activities show a strong positive correlation with each other (r = 0.508). However, no significant relationship was found between the antioxidant activity and total flavonoids content. Similar results were reported by Priastomo et al., (2024). The total protein content is negatively correlated with FRAP (r = -0.237). Similarly, the phytate content is negatively and weakly correlated with FRAP (r = -0.164). The positive correlation between DPPH, FRAP and total phenolic content in soybean seeds indicates that higher phenolic levels enhance antioxidant activity. This finding is consistent with a previous study that reported that the phenolic compounds are responsible for the level of antioxidant activity (Arifin et al., 2021). Studies show that as total phenolic content increases, both DPPH scavenging and FRAP values improve, reflecting a synergistic relationship where phenolics contribute significantly to the overall antioxidant capacity of soybean seeds (Alsuwayt, 2025). Flavonoids present in soybeans play a crucial role in improving protein quality (Oliveira et al., 2024). They can enhance protein digestibility and nutritional value by acting as antioxidants, which protect amino acids from oxidative damage (Prodić et al., 2023). By incorporating flavonoids into the diet, it is possible to potentially improve the overall quality of soybean proteins (Jia et al., 2022).

Table 4: Pearson’s correlation coefficient.


 
Association between nutritional and phytochemical contents
 
The first three principal components explained 69.60% of the total variance (Table 5; Fig 2), indicating a substantial representation of the biochemical diversity among soybean genotypes. The first principal component (Dim.1; eigenvalue = 2.44) accounted for 34.85% of the total variance and was strongly associated with antioxidant activities (DPPH, FRAP) and total phenolic content, reflecting the bioactive potential of the genotypes. Genotypes positioned in the positive quadrant of Dim.1 exhibited high antioxidant capacity and phenolic content and are therefore of interest for nutraceutical and functional applications. In this regard, GMV21, GPV20, GPV3, GPV8 and GMV20 displayed the greatest variability and represent promising candidates. The second principal component (Dim.2; eigenvalue = 1.30), explaining 18.62% of the total variance, was mainly influenced by total flavonoid and protein contents, indicating a nutritional and functional dimension. Genotypes located in the positive quadrant of Dim.2 combined high flavonoid and protein levels, making them suitable for breeding programs targeting both nutritional quality and health benefits. The most variable genotypes along this axis were GPV5, GPV12, GMV18, GMV16 and GMV9. The third principal component (Dim.3; eigenvalue = 1.12) explained an additional 16.12% of the total variance and was characterized by strong contributions from soluble sugars and phytic acid, highlighting variability related to carbohydrate content and specific antinutritional compounds. Genotypes GMV13, GMV12, GMV7, GMV11 and GPV1 showed the highest variability along this component. While genotypes located in the negative quadrants of Dim.1 and/or Dim.2 exhibited lower antioxidant or functional compound contents, some displayed high soluble sugar or tannin levels, indicating alternative valorization pathways.

Table 5: Eigenvalues and proportion of the variance explained for the three principal components of the 40 soybean genotypes based on phytochemical components.



Fig 2: PCA showing the spatial distribution of 40 soybean genotypes based on their Nutritional, phytochemical and antioxidant characteristics.


       
Globally, several genotypes (GMV5, GMV15, GMV21 and GPV8) showed positive loadings in at least two principal components, demonstrating their potential for targeted phytochemical improvement. The wide dispersion of genotypes confirms high biochemical diversity within the collection, allowing effective discrimination based on antioxidant, nutritional and phytochemical attributes and providing a robust basis for selection in breeding and utilization programs.
 
Cluster and heatmap analysis
 
A combined cluster and heatmap analysis were conducted to assess the variability among soybean genotypes concerning the total flavonoid content, total phenolic content, total protein, soluble sugars, phytic acid content and antioxidant activity. All genotypes were divided into three main groups (Table 6 and 7). In particular, Fig 3 shows the cluster analysis according to the average linkage method. Cluster 1 included only three varieties (GPV5, GPV4 and GPV7), which have high values in phytic acid content and low content in bioactive compounds (DPPH, FRAP, TPC, TFC). Cluster 2 also contained only two varieties (GPV6, GMV7) that have the highest total flavonoid content and total proteins. Cluster 2 is intermediate in antioxidant compounds and total flavonoids. Cluster 3 consisted of three varieties (GPV3, GPV8 and GPV20), which are prominent in the total phenol content and antioxidant activity. This cluster is rich in flavonoids and polyphenols and has high antioxidant activity (DPPH, FRAP), but low phytic acid content.

Table 6: Average cluster values of nutritional components, phytochemicals and antioxidant activities of 40 soybean genotypes.



Table 7: Composition and characteristics of the different groups.



Fig 3: Cluster and heatmap analysis based on significantly correlated variables.

In this study, for all the traits analyzed, genetic variation was primarily influenced by genotype, followed by the effects of the environment and the genotype × environment interaction (GEI). The Gampela locality was found to be less productive in terms of nutritional compounds, phytochemicals and antioxidant activities, whereas Farako-Ba was more productive. The variability among soybean varieties was also described through cluster analysis, which divided them into three groups. Moreover, correlation analysis indicated that total phenolic content was correlated with antioxidant activities. These results could contribute to the improvement of soybean breeding programs through an approach based on nutritional, phytochemical and antioxidant criteria. The soybean lines GPV12, GPV17, GPV22 and GMV19 are characterized by high nutrient content, particularly total protein and soluble sugars. The GPV3, GPV6 and GPV8 lines are distinguished by their high concentration of bioactive compounds, such as total phenolic compounds and total flavonoids. In contrast, the GPV4 and GPV5 lines have relatively high levels of antinutrients. Furthermore, this study makes a significant contribution to the mapping of Burkina Faso’s plant genetic resources.
 
The present study was supported by the International Center of Genetic Engineering and Biotechnology.
The authors have declared that no conflict of interest exists.

  1. Alghafis, A. and Raouf, E.A. (2020). Optimization of injection timing and injection duration of a diesel engine running on pure biodiesel SME (Soya methyl ester). Open J. Appl. Sci. 10: 486.

  2. Alsuwayt, B. (2025). Evaluation of antioxidant, anti-inflammatory and analgesic potentials of combined polyphenol-rich fractions from Ziziphus mauritiana and Ziziphus spina- christi leaves through modulation of inflammatory and oxidative stress markers in Sprague Dawley rats model. Inflammopharmacology. 33(9): 5509-5523. https:// doi.org/10.1007/s10787-025-01686-1.

  3. Amol, V., Bhati, K.R. and Bhati, K.R. (2021). Nutritive benefits of soybean (Glycine max). Indian J. Nutr. Diet. pp 522- 533. https://doi.org/10.21048/IJND.2021.58.4.27339.

  4. Arifin, H.A., Hashiguchi, T., Nagahama, K., Hashiguchi, M., Muguerza, M., Sakakibara, Y., Tanaka, H. and Akashi, R. (2021). Varietal differences in flavonoid and antioxidant activity in Japanese soybean accessions. Biosci. Biotechnol. Biochem. 85: 916-922.

  5. Arvouet-Grand, A., Vennat, B., Pourrat, A. and Legret, P. (1994). Standardization of propolis extract and identification of principal constituents. J. Pharm. Belg. 49: 462-468.

  6. Azam, M., Zhang, S., Qi, J., Abdelghany, A.M., Shaibu, A.S., Ghosh, S., Feng, Y., Huai, Y., Gebregziabher, B.S., Li, J., Li, B. and Sun, J. (2021). Profiling and associations of seed nutritional characteristics in Chinese and USA soybean cultivars. J. Food Compos. Anal. 98: 103803. https:// doi.org/10.1016/j.jfca.2021.103803.

  7. Bradford, M.M. (1976). A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal. Biochem. 72: 248-254.

  8. DuBois, M., Gilles, K.A., Hamilton, J.K., Rebers, P.T. and Smith, F. (1956). Colorimetric method for determination of sugars and related substances. Anal. Chem. 28: 350-356.

  9. Fu, J., Zheng, Y., Gao, Y. and Xu, W. (2022). Dietary fiber intake and gut microbiota in human health. Microorganisms. 10: 2507. https://doi.org/10.3390/microorganisms10122507.

  10. Hailemariam, H.M. (2022). Adaptability and stability for soybean yield by AMMI and GGE models in Ethiopia. Front. Plant Sci. 13: 950992.

  11. Hinneburg, I., Dorman, H.D. and Hiltunen, R. (2006). Antioxidant activities of extracts from selected culinary herbs and spices. Food Chem. 97: 122-129.

  12. Hungria, M. and Mendes, I.C. (2015). Nitrogen fixation with soybean: The perfect symbiosis? Biol. Nitrogen Fixat. pp 1009- 1024.

  13. Jia, Y., Fu, Y., Man, H., Yan, X., Huang, Y., Sun, S., Qi, B. and Li, Y. (2022). Comparative study of binding interactions between different dietary flavonoids and soybean â- conglycinin and glycinin: Impact on structure and function of the proteins. Food Res. Int. 161: 111784. https:// doi.org/10.1016/j.foodres.2022.111784.

  14. Jiang, G.L., Chen, P., Zhang, J., Florez-Palacios, L., Zeng, A., Wang, X., Bowen, R.A., Miller, A. and Berry, H. (2018). Genetic analysis of sugar composition and its relationship with protein, oil and fiber in soybean. Crop Sci. 58: 2413-2421.

  15. Kokebie, D., Enyew, A., Masresha, G., Fentie, T. and Mulat, E. (2024). Morphological, physiological and biochemical responses of three different soybean (Glycine max L.) varieties under salinity stress conditions. Front. Plant Sci. 15. https://doi.org/10.3389/fpls.2024.1440445.

  16. Latta, M. and Eskin, M. (1980). A simple and rapid colorimetric method for phytate determination. J. Agric. Food Chem. 28: 1313-1315.

  17. Lott, J.N., Ockenden, I., Raboy, V. and Batten, G.D. (2000). Phytic acid and phosphorus in crop seeds and fruits: A global estimate. Seed Sci. Res. 10: 11-33.

  18. Messina, M. (2022). Perspective: Soybeans can help address the caloric and protein needs of a growing global population. Front. Nutr. 9: 909464. https://doi.org/10.3389/fnut. 2022.909464.

  19. Njoroge, J.N. and Oyoo, M.E. (2020). Stability and adaptability of soybean [Glycine max (L.) Merrill] genotypes for yield, protein and oil content using AMMI analysis in Kenya. Agric. Sci. Dig. 40(2): 122-128. doi: 10.18805/ag.D-218.

  20. Oliveira, I.C., Santana, D.C., de Oliveira, J.L.G., Silva, E.V.M., da Silva, C.S.A.C., Blanco, M., Teodoro, L.P.R., da Silva Júnior, C.A., Baio, F.H.R., Alves, C.Z. and Teodoro, P.E. (2024). Flavonoids and their relationship with the physiological quality of seeds from different soybean genotypes. Sci. Rep. 14: 17008. https://doi.org/10.1038/ s41598-024-68117-z.

  21. Patil, S.S., Hiwale, S.D., Patil, P.P., Yadav, K.R., Patil, S.R. and Bisarya, D. (2024). Impact of levels of nutrient management on economics of soybean-onion and soybean-potato cropping system. Indian J. Agric. Res. 59(5): 725-730. doi: 10.18805/IJARe.A-6144.

  22. Priastomo, M., Adlia, A.R., Lumbantobing, V. and Adnyana, I.K. (2024). Determination of total phenols, total flavonoids and antioxidant activity of watermelon peel and rind from several cultivation areas in Indonesia. Indian J. Agric. Res. 58(5): 865-871. doi: 10.18805/IJARe.AF-872.

  23. Prodic, I., Krstic, R.M. and Smiljanic, K. (2023). Antioxidant properties of protein-rich plant foods in gastrointestinal digestion- peanuts as our antioxidant friend or foe in allergies. Antioxidants. 12: 886. https://doi.org/10.3390/antiox12040886.

  24. Singer, W.M., Lee, Y., Shea, Z., Vieira, C.C., Lee, D., Li, X., Cunicelli, M., Kadam, S.S., Khan, M.A.W., Shannon, G., Mian, M.A.R., Nguyen, H.T. and Zhang, B. (2023). Soybean genetics, genomics and breeding for improving nutritional value and reducing antinutritional traits in food and feed. Plant Genome. 16: e20415. https://doi.org/10.1002/tpg2.20415.

  25. Singleton, V.L., Orthofer, R. and Lamuela-Raventós, R.M. (1999). Analysis of total phenols and other oxidation substrates and antioxidants by means of folin-ciocalteu reagent, In: Methods in Enzymology. Elsevier. pp. 152-178.

  26. Sombié, P., Hilou, A., Mounier, C., Coulibaly, A., Kiendrebeogo, M., Millogo, J. and Nacoulma, O. (2011). Antioxidant and anti-inflammatory activities from galls of Guiera senegalensis JF Gmel (Combretaceae). Research Journal of Medicinal Plants. 5(4): 448-4461.

  27. Tolic, M.T., Krbavcic, I.P., Vujevic, P., Milinovic, B., Jurcevic, I.L. and Vahcic, N. (2017). Effects of weather conditions on phenolic content and antioxidant capacity in juice of chokeberries (Aronia melanocarpa L.). Pol. J. Food Nutr. Sci. 67(1): 67-74.

  28. Vucenik, I. and Shamsuddin, A.M. (2006). Protection against cancer by dietary IP6 and inositol. Nutr. Cancer. 55: 109-125.

  29. Wan, Z. (2023). Exploring the effects of genotype, environment and genotype by environment interactions on field pea (Pisum sativum L.) protein and amino acid contents using near-infrared reflectance spectroscopy. Food and Human Nutritional Sciences.
In this Article
Published In
Legume Research

Editorial Board

View all (0)