Abstract

Bambara groundnut (BGN) fits the bill when it comes to an acceptable level of nutrient and mineral composition. BGN is a balanced food that can help eradicate food and nutritional insecurity if it is incorporated into the major food system. However, there is a large degree of variation in nutrient composition and antinutritional factors among BGN accessions. Here, we show the degree of variability of nutrient and antinutrient components such as percentage ash, moisture, protein, fat, tryptophan, tannin, and phytate contents in seeds of 95 accessions of BGN. Data were subjected to analysis of variance (ANOVA), followed by correlation and principal component analysis. Clustering was done to show the relatedness between the accessions in response to the various traits. A high level of heterogeneity was observed among the accessions for the various traits studied. PC1 and PC2 show 41.2% of the total observed variations. Cluster analysis grouped accessions into four main clusters. This study was able to confirm the high level of diversity in the components of nutrients and antinutrients previously reported in BGN. The results of this study are expected to aid in identifying parent lines for improved breeding programs.

1. Introduction

The world’s population is predicted to exceed 9 billion by 2050 [1]. This, along with the impact of climate change, the unavailability of agricultural land, and the reliance on limited species of crop, has hampered the global achievement of food and nutrition security. In Sub-Saharan Africa and parts of Asia, the effect of climate change is erratic. Flooding, drought, and an increase in temperature are experienced in these areas [2, 3]. Farmlands are affected by environmental impacts from climate change effects such as temperature, humidity, rainfall pattern, and changes in light intensity and duration. Erratic yield and, in some cases, total loss have been the results [4]. Although these cannot be attributed to climate change, biotic factors such as pests, rodents, birds, and other small mammals have their fair share of recorded losses. The food being produced is merely meeting the demands of the ever-increasing population. Food and nutrition insecurity is taking its toll on less-developed nations. According to FAO, we will need to produce 60% more food if we are to sustain the population by 2050 (https://www.fao.org/sustainability/en/).

The nutritional and mineral components are vital to the survival of humans and livestock. Lack of an adequate amount of this can result in various fatal illnesses. The developing world is the most vulnerable to nutrient and mineral deficiencies. Nutrients are necessary for proper development, and an adequate supply is present in most crops, but most of these crops are significantly underutilized in the world food system. Most of these crops are an important part of their immediate local food system. One of these food crops is Bambara groundnut (BGN). Due to its high protein content (9.60–40.0%) [5, 6] and a good balance of essential amino acids [7], it is regarded as a complete food [8, 9]. Hence, it can be a good alternative to meat in terms of protein source [10]. In terms of nutrient components, it competes well with the major food crops, with very few boasting higher nutritional and mineral contents [11]. The high fiber and protein content also makes it an ideal component for animal feed. With all the promising attributes of BGN, it is reported to contain some antinutrient components such as phytates and tannins [12, 13]. These are not good for consumption in large quantities.

Variations in nutrient compositions occupy a unique role in achieving food and nutritional security in the developing world. Diversities in nutrient and mineral composition enhance the reputation of BGN as a complete food. However, the presence of antinutrient factors has a significant effect on the selection and use of some of the accessions. The high level of these factors in these accessions should be taken seriously, as antinutrient factors have been reported to affect the bioavailability of nutrients and minerals by chelating with the required minerals [14], thereby limiting the concentration of nutrients available for use. However, the situation can be taken care of by fermentation of the seeds and other processing methods as reported in some studies [1517].

In terms of trait improvement, no known varieties of BGN have been developed for improvement. The major hindrance to trait development includes a lack of a complete and well-annotated reference genome for the crop. However, other approaches can be used to improve the crop based on its nutrient and antinutrient composition, such as cooking processes [17]. One such approach is the use of a genome-wide association study (GWAS) approach. GWAS has been used to develop small nucleotide polymorphism (SNP) markers for regions of simple and complex traits in crops [18, 19]. Due to the success of this technology in other crops [19, 20], it would be of great interest if BGN could also benefit from improving its nutrient composition, especially by region mapping the antinutrient components. Locating the region mapping with the antinutrient traits can aid the removal of these traits using advanced technologies like gene editing, thereby producing a safe and nutrient-enriched complete food crop.

How well can BGN improve nutritional security, especially in Sub-Saharan Africa? Most people, especially in the developed world prefer the major crops for their source of nutrients and global acceptability, however, most underutilized crops reportedly contain enough of these nutrients as is found in the major crops. Therefore, in this study, we will the estimate the amount of nutrients and antinutrient components in selected accessions of BGN. Having this information will be key for nutrient-improvement programs in BGN, which will facilitate its incorporation into the food system in ensuring food and nutrition security in Sub-Saharan Africa and the world at large.

2. Materials and Methods

2.1. Source and Preparation of Bambara Groundnut Seed Materials

Ninety-five accessions of BGN seeds obtained from the genetic resource center of the International Institute of Tropical Agriculture (IITA), Ibadan, Nigeria, were selected. The seeds from each accession were milled into powder in the laboratory using a grinder [21]. 20 g of milled seeds were used for the analysis.

2.2. Nutrient Components
2.2.1. Moisture Content

3 g of the milled sample was weighed in duplicates and placed in moisture canisters. The samples were then placed in an oven (Memmert, GmbH, Model-30-750) for 16 hrs at a temperature of 105°C. The weight was taken every hour until it became constant. The percentage moisture content was calculated according to the formulawhere wt = weight.

2.2.2. Protein Content

Protein content was estimated from nitrogen content, which was determined by using the micro-Kjeldahl N distillation method followed by multiplication by 6.25. The value of 6.25 is the universally accepted protein content estimation factor. 0.2 g of each powdered sample was weighed into the digestion tubes of a Foss Tecator™ Digestor, and 4 ml of concentrated 98% H2SO4 and H2O2 were added. The copper tablet was used as the catalyst. The tubes were heated to 420°C in the digestion block until a clear solution was obtained. The resultant was distilled using Kjeltec 2200 distillation apparatus, and automated titer equipment was used for titration. The percentage of protein content displayed on the screen was recorded [22].

2.2.3. Fat Content Determination

The procedure of Official Methods of Analysis (AOAC, 1990) was followed using Foss Soxtec 2055 fat extraction and autoanalyzer equipment. 3.0 g of powdered sample was placed in a tube extracting petroleum ether and heated between 100 and 145°C; the reflux continues until the oil in the samples is leached out of the paper in a cup, leaving the water to evaporate. The oil was now allowed to cool down, and the volume of oil (fat) was obtained by the differential weight after evaporation as compared with the initial weight before evaporation occurred, as shown in the formula:where W3 = weight of the tube with the extracted oil, W2 = weight of the empty tube, and W1 = weight of the sample.

2.2.4. Ash Content Determination

The procedure for determining the ash content involves weighing 2 g of the powdered sample in an already known weighed crucible. They were then oven-dried for 4 hours at 105°C, after which they were put in a muffle furnace (Vulcan, Model-3-1750) at 550°C until white ash was obtained [22]. Then, the samples were cooled in a bench-top desiccator and reweighed to estimate the amount of ash. The percentage ash content was calculated:

2.2.5. Tryptophan Analysis

The determination of the amino acid tryptophan was done as a separate analysis because of its differing hydrolysis conditions. The samples (400 mg) are hydrolyzed under alkaline conditions with saturated barium hydroxide solution heated to 110°C for 20 hours. The hydrolysates were analyzed by reverse-phase liquid chromatography with UV detection at 285 nm using a Waters Breeze HPLC with Empower software (Waters, Millipore Corp., Milford, MA) [23].

2.2.6. Carbohydrate

It was derived from the subtraction of the results obtained from ash, protein, moisture content, and fat from 100, that is, 100 – (ash + protein + moisture content + fat) [24].

2.3. Antinutritional Components
2.3.1. Determination of Phytic Acid

This was analyzed with 1 g of the sample. Extraction and precipitation were done according to the method of Wheeler and Ferrel [25]. The phytic content of the samples was then determined by a 4 : 6 Fe/P atomic ratio.

2.3.2. Determination of Tannins

This was done according to the method of Adegunwa and Alamu [26]. A mixture of 0.5 g samples in 50 ml of distilled water was allowed to stand for 30 mins at 28°C and filtered using Whatman No. 42 filter paper. After dispensing approximately 2 ml of the extract into a 50 ml volumetric flask, standard tannin solution (2 ml) and 2 ml of distilled water were placed in separate volumetric flasks to serve as standards. Folin’s reagent and 2.5 ml of saturated NA2CO3 solution were added to each flask. The contents of each flask were then made up to 50 ml with distilled water and incubated at 28°C for 90 min. Their absorbance was measured at 260 nm on a spectrophotometer (ThermoFisher Scientific, Model G10 UV-VIS). The blank reagent was used to calibrate the instrument at zero.

2.4. Statistical Analysis

All analysis was done using the R statistical package [27]. Descriptive and summary statistics were computed using the stat.dev function in the pastecs package. One-way analysis of variance (ANOVA) was performed using the lmer function to determine significant differences between the means of the nutrients and antinutrient components. Fischer’s least significant difference (F-LSD) at a probability level of 5% was used to separate the means that were significantly different. The principal component analysis was done with subsequent cluster analysis using the fviz function in the FactoMineR package and correlation coefficients to view the groupings and relatedness among these collections of accessions on the different traits.

3. Results

3.1. Descriptive Statistics and ANOVA Analysis

Descriptive statistics showed a significant coefficient of variation (CV > 20) in ash content, moisture content, tryptophan, and tannin contents (Table 1, Figure S1). The ranges of 50.6–69.3, 0.76–6.39, 2.88–13.55, 12.51–26.72, 3.19–9.88, 0.09–0.47, 0.08–0.42, and 3.17–7.42 were recorded for carbohydrate, ash, moisture, protein, fat, tryptophan, tannin, and phytate content, respectively. Among the nutrient components, carbohydrate has the highest range but lowest variation, while phytate has the highest range between the two antinutrient components and the lowest variation. In the ANOVA result (Table 2), all components were highly significant (). The resulting means separation showed a high level of heterogeneity among the accessions for the traits, especially for ash, moisture, phytate, and tryptophan contents (Table 3).

3.2. Principal Component Analysis

When we looked at the variation in nutrient and antinutrient content among the 95 BGN accessions, the first three principal components (PC) with eigenvalues ≥1 described around 62% of the variation (Tables 4, and . The first, second, and third PCs accounted for 29.74%, 17.55%, and 14.45% of the variations, respectively, in the dataset under consideration. The principal components (PC) are represented as dimensions (dim.). PC1 (Dim.1) was mostly influenced by the following characteristics: CHO (40.97%), moisture content (18.58%), protein (16.87%), ash (10.51%), and fat (9.67%). In the second PC (Dim.2), the variables ash, moisture content, tannin, and protein each accounted for 21.33% of the variations, 16.91% of the variations, 42.66% of the variations, and 11.56% of the variations. Tryptophan (48.70%) and phytate (40.95%) were the most important contributors to the variation described by PC3 (Dim.3).

The various links between accessions and attributes are depicted in Figure 1. These accessions whose vectors are plotted close to the vector for a specific nutrient or anti-nutrient component are highly correlated to that nutrient or antinutrient component, whereas the length of a vector for an accession will estimate its mean for that specific nutrient or antinutrient component. An accession that is drawn along with a lengthy vector implies that it has a high mean value for the nutrient or antinutrient that it relates to. In the study of CHO, three accessions were among those identified as being strongly connected with the vector: TVSu-2065, TVSu-2071, and TVSu-2086. A few accessions associated with ash were identified, including TVSu-1859, TVSu-1823, and TVSu-2046, among others. In addition, TVSu-1943, TVSu-1930, and TVSu-1979 were among those identified as being near the vector for moisture content, and the accessions associated with protein were identified as being TVSu-2048, TVSu-2021, and TVSu-1923, among others. There were other accessions in close proximity to the fat vector, including TVSu-1905, TVSu-1952, and TVSu-1899, all of which were identified as such. In terms of tryptophan content, accessions TVSu-2031 and TVSu-1764 were among those found to be associated, whereas accessions TVSu-1918, TVSu-1892, and TVSu-2025 were among those found to relate to tannin content. Phytate is associated with the accession TVSu-1863, which is one of the accessions associated with it.

3.3. Clustering Based on Analyzed Traits

The 95 accessions may be divided into four distinct clusters based on their nutritional and antinutritional features, according to hierarchical clustering (Figure 2). Three clusters with a similar number of accessions were identified. The first cluster had thirteen accessions, whereas the second cluster, with thirty accessions, was the largest. With twenty-five and twenty-seven accessions, the third and fourth clusters followed closely behind the second cluster. The first cluster was defined by accessions with high CHO mean values but low MC mean values. Accessions in the second cluster had high mean values for ash, fat, and phytate content, whereas accessions in the third cluster had low mean values for phytate and fat but a high mean value for MC, and accessions in the fourth cluster had high mean values for tryptophan and tannin (Figure 2).

3.4. Correlation among Nutritional and Antinutritional Components

In the bivariate correlations, it was discovered that the nutrients had varying degrees of association with one another (Figure 3). Asterisks indicate that the correlations are statistically significant at a value of less than 0.05. Positive correlations are represented by the blue hues, whilst negative correlations are represented by the red colors. The deeper the hue, the greater the degree of negativity or positivity represented by the matrix. Phytate was strongly linked with protein (r = 0.58, ), fat (r = 0.28, ), and ash (r = 0.58, ) in the correlation matrix (Figure 3). Tryptophan correlated negatively with protein content (r = -0.17, ), tannin (r = -0.10, ), and moisture content (r = 0.43, ). Fat correlated positively with moisture content (r = 0.43, ) and tannin (r = 0.16, ). The moisture content and CHO (r = −0.83) had the largest negative connection, followed by protein and CHO (r = −0.77), which had the second highest negative correlation. A positive correlation was found between phytate and protein (r = 0.23, ) as well as tryptophan (r = 0.08, ), while a negative correlation was found between phytate and tannin (r = 0.21, ). This means that the amount of phytate is not affected by the amount of tannin present.

4. Discussion

Food insecurity in the developing world often results from the low nutritional quality of some of the traditional foodstuffs and the high cost of protein foods. Fortified foods are not easily available for most families in the developing world [28] due to the high level of poverty. Legumes have therefore come under immense importance in these regions as they have been proven to be alternative sources to fortified foods for nutrients. Hence, it is of necessity to increase their level of consumpton.

4.1. Genotypic Variability in Nutrient and Antinutrient Contents in BGN

The genetic, physiological, and morphological variations of plants have a direct impact on breeding outcomes. Therefore, the importance of variations in different traits cannot be overemphasized. There are significant variations in nutrient and antinutrient contents among the accessions studied (Tables 2 and 3, Fig S1), therefore, they are excellent resources for developing or selecting parental lines for better nutritional performance.

Carbohydrates are needed for energy and are present in the form of starch and sugars, varying from simple sugars to complex sugars, i.e., the monosaccharides, disaccharides, and polysaccharides. According to Liu, Tang [29], pulse seeds are high in carbohydrate content. Therefore, it is expected that BGN should contain a high amount of carbohydrate. In this study, the carbohydrate content of 58.34% is similar to the result reported in the study of Atoyebi and Osilesi [30] (range of 42.77%–62.76%) and higher than 50.2% reported by Yao and Kouassi [7]. The higher the moisture content, the higher the susceptibility of the seeds to attacks from microorganisms [28]. Therefore, those seeds with low moisture content will have a longer shelf life than those with high moisture content. Pulse seeds contain 20%–40% protein [31, 32], a value twice as much as that present in cereals. The value for protein content in this study, which is below 25%, is similar to that obtained by Nti [33], who reported a range of 19.3% to 27.1%, and Sirivongpaisal [34], who reported 15.48% protein content. Compared with soybean, BGN is lower in protein content, as reported in the study of Wijewardana and Reddy [35] on soybean seeds under well-watered conditions. However, these values fall in the range obtained for protein content in other pulses such as peas, chickpeas, beans, lupin, and cowpea [32, 36]. Within the Vigna genus, BGN mean concentrations of the nutrient components are similar to those reported in cowpea and mungbean [13, 37, 38]. Tryptophan is a precursor of auxin, which is an essential growth hormone in plants [3941], and bioactive compounds such as nicotinamide, melatonin, and kynurenine in humans [42, 43]. The tryptophan levels recorded in this study show that BGN can help get the required amount of tryptophan in the quantity recommended by WHO, which is 1.70, 0.85, and 0.66 for infants, children, and adolescents, respectively [44]. The tannin content reported in this study (0.09%–0.41%) is lower than 0.25%–2.27% reported by Ofori and Kumaga [45]. However, the mean value of 0.20 is greater than 0.046 reported by Mazahib and Nuha [46] in Sudan and 0.039 reported by Ijarotimi and Esho [47] in Nigeria. However, it is lower than the tannin content in faba beans, which is between 8 and 9% [14]. Tannins are known to be present in food products and to inhibit the activities of trypsin, chymotrypsin, amylase, and lipase, as well as lower protein quality and interfere with iron absorption [48]. Tanning can be detrimental to microbial enzyme activities, including cellulose digestion and intestinal digestion, if the tannin concentration in the diet becomes too high [49]. The utilization efficiency and ease of digestion of nutrients are affected by the tannin content present [50]. The strong antinutritional compound phytate is present in nearly all legumes and seeds [48]. Phytic acid inhibits dietary tyrosinase, trypsin, pepsin, and lipases [51]. The average phytic acid level reported by Mazahib and Nuha [46] in Egyptian landraces was 14.78 mg/100 g, which is higher than that reported in this study. The mean phytic acid content reported in this study is higher than that reported in the study of Yao and Kouassi [7], who reported a mean value of 1.1 mg/100 g. Phytic acid binds to ion complexes, resulting in a decrease in the concentrations of these ions [52]. Most of the phosphorus in phytic acid is unavailable to monogastric animals because they do not have the enzyme phytase [53]. Furthermore, dehulling in some pulses increases their phytic acid content but decreases their tannin content because tannins are mainly found in the seed coats [54].

4.2. Association among Accession, Nutrient, and Antinutrient Components

In addition, the PCA biplot showed four clear separations in the biplot according to the nutrient components (Figure 1). Accessions high in carbohydrates are also high in tryptophan, while those high in tannin content are also high in moisture and fat content, but they are low in carbohydrate, tryptophan, protein, ash, and phytate content. Likewise, those high in protein content are also high in ash and phytate content but low in the other components. Considering the antinutrients, those high in phytates are low in tannins and vice versa. Tryptophan is an amino acid constituting total proteins. They are therefore part of the total protein; hence, accessions high in tryptophan are expected to be high in protein and correlate positively with total protein. In this study, however, accessions high in tryptophan were low in total protein and vice versa.

The hierarchical clustering (Figure 2) agreed with the biplot result (Figure 1), which showed that accessions with comparable nutritional and antinutrient contents were clustered together. Divergence among BGN accessions has been successfully identified using clustering based on nutritional and antinutrient components [7, 4547, 5558]. To avoid inbreeding depression, the selection of accessions for crossing should focus on divergent accessions, while proposing a few accessions for production might reduce the risk posed by biotic and abiotic pressures, particularly in stress-prone areas such as Sub-Saharan Africa.

4.3. Association between Nutrient and Antinutrient Components

Finally, the correlation between the nutrient and antinutrient components showed that breeding for nutritional components in these accessions will require careful selection processes. Although there are many significant correlations, the numbers of positive correlations are however similar to those of negative correlations. The negative relationship between carbohydrates and protein aligns with the result reported in the study of Ndidi and Ndidi [59]. However, the opposite was reported for the relationship between ash and phytate. The correlation between tannin and phytate reported in this study is supported in the study of Unigwe and Doria [60].

The environment interacts with different plants to give different responses, which are observed in the nutrient, antinutrient components, and other traits, such as reported in the studies of Ribeiro and Kläsener [61], Fauziah and Mas’ udah [62], Vogelsang-O’Dwyer and Petersen [63], Olanrewaju and Oyatomi [64], and Olanrewaju and Oyatomi [65]. Because these traits are regulated by specific genes. The impact of the environment on these genes to affect the nutrient components was not considered in this study. Therefore, nutrient and antinutrient components from plants grown in various environments should be studied to ascertain the extent of the environmental influence on the nutrient and antinutrient components. Secondly, the emergence of NGS has revolutionized molecular studies in various areas; therefore, gene locus regulating these traits can be identified using technologies such as GWAS for improvement breeding programs.

5. Conclusion

In this study, two of the traits studied, viz. carbohydrate and tryptophan, are high in more than half of the accessions, while the others have high amounts of the remaining traits, including the antinutrient factors. Therefore, the selection of accessions having low or minimal amounts of antinutrient components can be made possible for future improvement. The variation in nutrient and antinutrient components in BGN observed here agrees with previous findings. The results of this study offer opportunities for the selection and breeding of the crop for beneficial nutrient compositions and, at the same time, those with low amounts of antinutrient factors. Therefore, BGN can successfully improve nutritional security even more than some major crops.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The first author is gratefully acknowledged for school bursaries and research fellowships, respectively, from North-West University and IITA. The assistance rendered by the staff of the Genetic Resource Center, IITA, is greatly appreciated. The work of GRC was funded by the Global Crop Diversity Trust and CGIAR.

Supplementary Materials

Figure S1: nutrient and antinutrient compositions of each accession. Table S1: contribution of the traits to each PCA component. (Supplementary Materials)