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).
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).
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.
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).
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.
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.