Abstract

Bambara groundnut (Vigna subterranea L. Verdc.) is considered an emerging crop for the future and known as a crop for the new millennium. The core intention of this research work was to estimate the variation of landraces of Bambara groundnut considering their 14 qualitative and 27 numerical traits, to discover the best genotype fitted in Malaysia. The findings of the ANOVA observed a highly significant variation () for all the traits evaluated. There was a substantial variation (7.27 to 41.21%) coefficient value, and 14 out of the 27 numerical traits noted . Yield (kg/ha) disclosed positively strong to perfect high significant correlation ( to 1.00; ) with traits like fresh pod weight, dry pod weight, and dry seed weight. The topmost PCV and GCV values were estimated for biomass dry (41.09%) and fresh (40.53%) weight with high heritability (Hb) and genetic advance (GA) %, % and %, %, respectively. The topmost heritability was recorded for fresh pod weight (99.89%) followed by yield (99.75%) with genetic advance 67.95% and 62.03%, respectively. The traits with and suggested the least influenced by the environment as well as governed by the additive genes and direct selection for improvement of such traits can be beneficial. To estimate the genetic variability among accessions, the valuation of variance components, coefficients of variation, heritability, and genetic advance were calculated. To authenticate the genetic inequality, an unweighted pair group produced with arithmetic mean (UPGMA) and principal component analysis was executed based on their measurable traits that could be a steadfast method for judging the degree of diversity. Based on the UPGMA cluster analysis, constructed five distinct clusters and 44 accessions from clusters II and IV consider an elite type of genotypes that produce more than one ton yield per hectare land with desirable traits. This study exposed an extensive disparity among the landraces and the evidence on genetic relatives will be imperative in using the existing germplasm for Bambara groundnut varietal improvement. Moreover, this finding will be beneficial for breeders to choose the desirable numerical traits of V. subterranea in their future breeding program.

1. Introduction

The Bambara groundnut (Vigna subterranea L. Verdc.; Syn: Voandzeia subterranea L. Thouars) is an underutilized grain legume that belongs to the family of Fabaceae and subfamily of Faboidea grown mostly in Africa [1]. Verdcourt [2] suggested the present binomial name Vigna subterranea (L.) Verdc and its chromosome number is [3]. Bambara groundnut is a future emerging legume grown in Africa and Asia, is commonly referred to as a poor man’s crop, mostly known as “Women’s Crop” for family food security [1], and recently noticed as the crop for new millennium [4]. The Bambara groundnut is treated as the 3rd most important legume crop after groundnuts (Arachis hypogaea L) and cowpeas (Vigna unguiculata L. Walp.) in Africa [5] but due to its low rank, it is considered a snack or food supplement but not a lucrative cash crop [6]. The center of origin of Bambara groundnut is believed to be ‘Bambara,’ a place name near Timbuktu in central Mali, West Africa [7], and suffix ‘groundnut’ is because the process of its pods grows under soil, which is equal to peanut/groundnut; hence, its common name is ‘Bambara groundnut.’ The crop has been vastly cultivated in tropical areas, now available in several parts of South America, Asia, and Oceania [8]. Effectively, it is cultivated for human utilization and has been noted as a fully balanced diet due to a high content of carbohydrate (63-65%), protein (18-20%), and oil (17-18%) in its seed [9, 10]. The biochemical investigation reported that it possesses essential amino acids (33.31%) and nonessential amino acids (66.69%) out of the total content [11, 12]. Bambara groundnut has great potential for incorporation into different foods of human where it possesses [11] crude protein (17.5 to 21.1%), crude fat (7.3 to 8.5%), total ash (4-5%), crude fiber (1.8-2.0%), CHO (53.0 to 60.8%), and moisture (7.5 to 12.3%). The mineral content of Bambara groundnut seeds that was accounted for (mg/100 g dry matter) the macro minerals are Ca (37-128), K (1545-2200), Mg (159-335), Na (16-25), P (313-563), and for the micro minerals (ppm) Cu (3.0-13.2), Fe (23.0-150) and Zn (13.9-77.0) reported by Amarteifio et al. [12]. Due to the high content of iron and protein with a significant level of lysine (10.3%), methionine seeds of Bambara were marked as uniformed food sources when compared to other food legume crops [13]. For its rich protein content, Bambara groundnut fulfills the daily demand of proteins for the low-income users and also has the potentiality to develop nutrition and food security and accelerate rural improvement and sustainable use of land where animal proteins are not freely available for service as a staple food for them [14]. In some areas of Nigeria [15], the young fresh seeds were boiled and eaten as a manner of the way like boiled peanut was also made a pudding known as Moi-Moi or Okpa or porridge bean. Bambara groundnut is used for making bread [16], and roasted seeds are eaten as confectionery [17], pounded and mixed with soup [18]. Habitually, the use of raw bean in a manner of chewing and swallowing prevent nausea [19], green leaves are eaten for antivomiting [9, 20], and its fodders are used to feed for animals [21]. Bambara groundnut was treated as a source of fully balanced food [22] and permitted to grow in drought with various agroecological aspects resulting to it becoming an important economic crop for developing countries [23]. Same as other legume crops, it can fix atmospheric nitrogen [24], mostly grown by the female [25] and able to produce a high yield with low input. In Swaziland, it is reported that about 98% of farmers regard Bambara groundnuts as a profitable crop [26] and most of the portion of total yield is used by themselves; the rest the amount (10-40%) was sold to the local market. Globally, the estimated annual production of Bambara groundnut was 160,378 tons in which 111,562 tons are produced in West Africa [27]; besides, Burkina Faso occupied the major portion in Africa with Nigeria leading its production at 100,000 metric tons per annum [28]. However, the existing landraces produce low yield (650 kg/ha) due to lack of following appropriate farming methods, diseases-insect infestation, and lack of improving genotype best adjusted to climate change [1, 29]. On the other hand, some researchers [30] have noted that improved Bambara groundnut genotypes can produce yield 3.0 t/ha to 4.5 t/ha when all factors related to yield are favourable in condition. The adequate knowledge of several genotypes and their evaluation is obligatory for germplasm selection as well as their enhancement approaches [31]. The global population growth rate increases by an estimation of 80 million per year and assume to reach 9.2 billion by 2050 [32]. So, based on the current situation use of potential genetic resources for plant breeding to boost up the production of this crop [33] which can provide as a supplement to meet up a certain defect in the consumption of major crops like rice, wheat and maize also enhance food security in developing countries. Baudoin and Mergeai [8] reported that the Bambara groundnut is an extreme autogamous crop with cleistogamous flowers [34]. Many researchers noted that the effective hybridization between two different lines of Bambara groundnut through traditional breeding has not yet been achieved [13, 35]. The traditional breeding technique for the enhancement of the Bambara groundnut is slow and problematic due to the long generation time and preponderant homozygous nature of the crop [36]. In addition to geotropic pod formation of the Bambara groundnut, it makes its trouble to artificial hybridization [37]. The breeding technique of Bambara groundnut is undefined, and there are no high-yielding cultivars available in Malaysia and local landraces of this crop are still grown. Henceforth, the readily accessible development approach is to apply the selection technique to the already existing variable accessions. The current demand is to discover modern high-yielding cultivars for certain growing regions [38]. Morphological characterization is the first footstep of the germplasm investigation to identify desirable traits of interest [39]. Despite the versatile advantages of Bambara groundnut cultivation, a few findings have experimented on this edible crop in Malaysia compared to the other legumes like sorghum, groundnut, and cowpea. But like other African countries, the Bambara groundnut is vastly grown in Nigeria as a legume crop [40]; besides the African continent, the secondary center of cultivation of Bambara groundnut is in the Asian region like Sri Lanka, Malaysia, Philippines, India, and Brazil [21]. Bambara groundnut can well adapt to the tropical area like Malaysia, where the cultivation of major crops (rice, wheat, maize, etc.) are increasingly challenging due to drought and unpredictable rainfall patterns [41]. This research work emphasize the morphological performance of 150 Bambara groundnut accessions with a view of exploring the variation that exists among the traits with an aim of selection for high yield and further help to the selection of elite genotypes in breeding and agricultural improvement programs. In the current background of global climate change, one of the best approaches to reduce the hereditary erosion of the Bambara groundnut is the germplasm collection and diversity analysis in different growing regions. It could assist to identify the ongoing cultivating landraces as well as to inaugurate the way of management and upgrading of varietal characters. Consequently, this research discovered that the morphological divergence exists in Bambara groundnut accessions in Malaysia. Currently, the available modern strategies are applying for the selection of ongoing cultivated species of Bambara groundnut. After all, the trait improvements that can be made through direct selection are estimated by the availability of heritable variation. Furthermore, heritability influences the magnitude of the selection procedure which would be a powerful tool to improve a certain trait and predict the genetic gain from selection also estimating the comparative effect of genes [42]. Therefore, the pinpoint of this research was to discover the genetic divergence in different qualitative and numerical traits of Bambara groundnut with a certain goal (like high yield) of determination of variance component, heritability, genetic advance, and clustering, based on the selection intensity of the accessions with potentially high-yielding criteria.

2. Materials and Methods

2.1. Experiment Location

This research has been conducted from September 2018 to February 2019 at field-15 at University Putra Malaysia. The research work has experimented under the Institute of Tropical Agriculture and Food Security (ITAFoS), University Agricultural Research Park, University Putra Malaysia (UPM), Malaysia. Based on the Global Positioning System (GPS), the research location was 2°5854.0N latitude and 101°4253.8E longitude. The seeds of accessions were sown in open field conditions during the 2018-2019 cropping season. During the planting season, the mean climatic situation is stated in Table 1. The soil pH is 6.6 to 7.5 with sandy loam to clay loam type (Dept. of land management, UPM).

2.2. Genetic Materials

One hundred fifty accessions of Bambara groundnut were selected for this current research work, all representing the African accessions collected from the local market of Nigeria. The list of Bambara groundnut accession used in this research was displayed in Table 2. To estimate the genetic divergence using the morphophysiological traits, five randomly plants were taken for data investigation among the accessions evaluated [43].

2.3. Experimental Design

The experiment was conducted in a randomized complete block design (RCBD) with three replications. The experimental plot comprised of single rows measuring . The distance between the plant to the plant was 30 cm and row to row distance was 1 m and the distance between replication was 2.0 m according to [43]. During the growing season, the recommended intercultural practices like land preparation, land clearing, weeding, irrigation, and fertilizer were approved. Hand sowing was done into a raised seedbed, with two seeds in a single hole in 3 cm depth, and the seedlings were uprooted to keep in one plant per hole after 2 weeks of sowing when seedlings are completely stable with soil. The recommended fertilizer rates (100% , 100% , and 100% ). The total portion of phosphorus (100% P) and potassium (100% K) was applied during land preparation; hence, 70% N was applied at 5 weeks after sowing [44].

2.4. Parameters Measured for Data Analysis

Both qualitative and quantitative data were taken following the Bambara groundnut descriptors [45]. Twenty-seven measurable characters (Table 3) and fourteen qualitative characters (Table 4) were considered during the morphological characterization. The measurable traits were divided into the following three categories: (1) phenological traits, (2) growth and vegetative traits, and (3) yield traits to easy interpretation. All parameters were visually recorded at different growth stages of five plants in the field and after harvest in the lab as per the following description and descriptors states [45].

2.5. Statistical Analysis

The SAS (statistical analysis software) version 9.3 was followed to test the significant differences using the analysis of variance (ANOVA) procedure at the level of LSD () and to compare among the means of significant traits. The result of this research was expressed as mean, the square of mean, genetic parameter, and lastly, the Pearson correlation was measured to find out the intercorrelation ships among the traits. The correlations between the quantitative variables were determined using the Pearson [46] correlation coefficient formula. A report from [37] noted that correlation is an appropriate guide, particularly for plant breeders who may want to associate a set of traits in their selection programs. Typically, the correlation studies among numerical traits are of great value to plant breeders in selecting elite traits.

3. Genetic Parameter Analysis

3.1. Estimation of Covariance, Broad-Sense Heritability, and Genetic Advance

(1)Analysis of variance (ANOVA) estimation: this analysis was done to discover the uniqueness among the accessions and estimate the effects of environment and genes on several traits(2)The genotypic and phenotypic variation was calculated as per following the formula given by [47]: where is the genotypic variance, is the phenotypic variance, MSG is the genotypic mean square, MSE is the error mean square, and is the replication number(3)The coefficient variation of phenotypic (PCV) and genotypic (GCV) were estimated as per formula given by [48]: (a)(b)(c)where PCV is the phenotypic coefficient of variation, GCV is the genotypic coefficient of variation, is the grand average of the traits, is the phenotypic variance, is the genotypic variance, and RD is the relative difference between PCV and GCV.The estimated values of PCV and GCV were categorized by [48, 49] like as between 0%-10% for low, 10%-20% for intermediate, and greater than (≥20%) for high.(4)Broad sense heritability (): it refers to the proportion of genotypic variance () with the phenotypic variance () multiplied by a hundred. For estimation of (), the formula given by [50] was followed: where = Genotypic variance and = Phenotypic variance. In accordance with [51, 52], the heritability grade was ordered between 0% and 30% for low, 30% and 60% for intermediate, and greater than 60% as high.(5)Genetic advance (GA) (as a percentage of mean) was calculated with a 5% selection intensity () following the method of [51]. Genetic advance is categorized as between 0% and 10% for low, 10% and 20% for intermediate, and more (>20%) than for high, following the formula given by [53]: where is the constant that indicates the intensity of selection. According to Adewale et al. [54], the rate is 2.06 at the point when the is at 5%. is the standard deviation of phenotype, is the broad sense heritability, and is the grand mean values of traits.(6)Genetic gain (%): estimated as 0; it is also categorized [51] as between 0 and 10% for low, 10 and 20% for intermediate, and ≥20% for high genetic advance

Further, to determine the expected gain from selection, estimation of heritability simultaneously with the values of the genetic advance can be an effective tool of crop enhancement program. Valuation of this component is an extensively fundamental step that must be taken into consideration before commencing any breeding program. In our current research, the enhancement of Bambara groundnut yield preference of the selection method was taken based on the magnitude of variation that exists in the gene pool of this crop, measurement of variance component, heritability, and genetic advance.

3.2. Multivariate Analysis

To determine the genetic or hereditary divergence, cluster analysis was used for 27 numerical traits that were considered in this study. Based on the Euclidian distance method, data was analyzed for investigation of genetic diversity. In addition to this, based on the unweighted pair group method using arithmetic average (UPGMA) and following the algorithm and sequential, agglomerative, hierarchic, and nonoverlapping (SAHN) method, the genetic interrelationship among the Bambara groundnut was estimated using the SAS version 9.3 software. Using similar software, the principal component analysis (PCA) was done with subsequent correlation coefficients applied to construct a dendrogram to observe the groupings and relatedness among the accessions of the Bambara groundnut. The NTSYS version 2.1 (Numerical Taxonomy Multivariate Analysis System) and Exeter Software (Setauket, NY, USA software) [55] were used to produce two-dimensional (2D) plots for principal component analysis (PCA). Modena’s [56] stopping rule was followed to select the number of clusters with Milligan and Cooper’s [57] correction.

3.3. Shannon Diversity Index () and Evenness ()

Diversity indexes are statistics regarded to summarize the diversity of a population in which each member belongs to a unique group. Shannon’s diversity index () is another index that is generally used to categorize the species diversity in a certain community. Shannon’s diversity index is an account for both richness and evenness present in the species also used for a wide diversity of fields. It is also known as phylogenetic indices or phylogenetic metrics, which is a numerical estimation that indicates how many types (such as species) of variation are present in a community and simultaneously can consider the phylogenetic relations among the individuals. However, the Shannon evenness index is a synonym for the Shannon equitability index and was calculated [58, 59] using the formula as follows: where is Shannon’s diversity index, is the maximum diversity possible, ln is the natural logarithm of a number, pi is the proportion of the population made up of species ‘’, is the number of species or species richness in a sample, is the number of total samples, and .

4. Result

4.1. Qualitative Diversity

The one hundred fifty landraces of Bambara groundnuts were taken from the local market of Nigeria in different growing regions. The accessions were grouped in 11 units (Table 5) by morphotypes (Figure 1), indicative therefore of the existence of many duplicates in the landraces based on the seed morphological description. The morphological characterization (Figures 2 and 3) summarized the frequency of distribution of some qualitative variables studied in this research. After two weeks, 56.66% of the accessions had greenish stems, 26% had stripped stems, and 17.33% were reddish stems. The terminal leaflets had three different colors: 65.33% accession had a greenish leaflet while the purple and red accessions were both 17.33%. 34.66% of the total accessions had terminal leaflets shaped like lanceolate whereas 48% had oval and 17.33% had elliptic in shape. Among the 150 characterized accessions, three growth habits were found (Figures 2 and 3): bunch-type accessions (26%), semibunch-type accessions (34.66%), and the spreading type (39.33%). Among the landraces, 36.66% had sparse hair on their stems and 17.33% had dense hair while 46% did not have any hair on their stems. Most of the landraces had reddish-brown (45.33%) and brown (35.33%) color pods; some had yellowish-brown (8.66%) and purple (10.66%) color pods. In this current research the seed’s color and texture also taken into consideration and showed significant variation. Maximum accessions were found round (61.33%) shaped and few were oval (38.66%). Seed color had cream and red each of 26%, black-cream and cream-purple each of 19.33%, only 9.33% had black color. 28% of landraces had black eye color and 72% had no eye color. There was no testa pattern (44%) but entire line striped marbled (30%), enough rhomboid spots on both sides (8.66%) of the hilum, entirely dotted spot (8.66%), while little rhomboid spots on both side (8.66%) (Figure 2). Concerning the pod texture, most of the accessions (64%) had smooth little grooves while 26% had enough grooved and only 8.66% had enough folded. The existence of a significant morphological variation was detected for all the qualitative traits such as stem hairiness, terminal leaflet shape, growth habits, pods, and seed color.

4.2. Quantitative Traits
4.2.1. Morphological Diversity

Most of the plant breeders treated yield and other yield contributing traits as high influential parameters for crop improvement. Generally, the traits which are associated directly and indirectly with yield are pod and seed size; shape; quality; plant height; branch number; the total number of pods, plant, and pod biomass weight; resistance to diseases and insect-pest infestation; and biotic and abiotic stresses that also have a significantly remarkable importance in Bambara groundnut breeding programs. However, in this current research, a total number of 27 quantifiable traits of 150 Bambara groundnut accessions were analyzed for the selection of best genotypes with high-yielding capacity. Significant variation, mean, standard error of the mean (SEm), standard deviation (St. Dev), and coefficient of variation (CV%) revealed by analysis of variance among the investigated 27 quantitative traits were displayed in Table 6. The analysis of variance (ANOVA) of 150 landrace’s quantitative traits revealed highly significant () for all the traits studied among the accessions. For replication, all the variables showed a highly significant difference () except the variables no. of branch per plant, fresh biomass weight per plant, immature pod per plant, and shelling percentages observed no significant differences. The minimum and maximum values across overall plants were shown in Table 6, and the observed coefficient of variation (CV%) values ranged from 7.27% (days to maturity) to 41.21% (biomass dry weight per plant). The average days to maturity were found to be days which is statistically significant (). The average no. of branches per plant with the standard error was observed , for total no. of pod per plant , for hundred seed weight and yield  kg/ha. The highest value of standard deviation (SD) was found for the trait yield kg/ha () with standard error (SEm: ±14.93) while the lowest was for internode length (; SEm: ±0.03) (Table 6). We found significant differences in coefficient ( for 14 traits out of the 27 numerical variables studied in this research. Days to 50% flowering varied from 25 to 52 days after sowing (DAS). We noticed that 70.66% of the accessions gave flowers before 40 DAS while 29.33% of the accession produced before 50 DAS. Most of the landraces (90%) had a life cycle of more than 120 days, and it was only one accession S1G141(0.66%) that took an average 160 days to maturity. Hundred seed weight varied between 196.21 g and 364.09 g whereas average highest hundred seed weight was calculated for the accessions S1G143 (329.58 g) and the lowest was S1G13 (203.72 g). As for the yield, ranged from 370.38 to 1679 kg/ha across the plant whereas topmost mean yield was recorded for the accession S1G92 (1635.29 kg/ha) while lowest was S1G28 (380.48 kg/ha). Lowest least significant difference () was noticed for the trait internode length 0.36 while the highest was 35.28 for the trait no. of leaves per plant (Table 6).

4.2.2. Performance of the Several Morphotype’s Elite Accessions

Twenty-three accessions were notified as elite stranded among the whole morphotypes based on a higher yield performance and other yield contributed quantitative traits (Table 7), and the at-a-glance relationship of dry pod weight (g) and a hundred seed weight (g) with yield (kg/ha) is displayed in Figure 4. All the 23 elite accessions gave the best field yield of more than one ton per hectare. The accessions S1G92 (1635.29 kg/ha) produced maximum yield followed by S1G93 (1632.87 kg/ha) and S1G32 (1626.41 kg/ha) among the accessions were studied. In addition to the high yield, the accession S1G108 (1369.01 kg/ha) was identified with a short life cycle of 106 days after sowing (Table 7).

4.3. Analysis of Correlation (Association) Matrix

The phenotypic association among the 27 numerical traits of one hundred fifty Bambara groundnut accessions is given in Table 8. No significant association was found for days to emergence with yield. Fifty percent flowering days had negative and weak () significant correlation with yield (; ). Days to maturity had positive and intermediate () highly significant association with yield (; ). Plant height had positive and weak () highly significant relation with yield (; ) but intermediate correlation with seed length (; ). Significantly positive and intermediate association () was found between biomass fresh weight (; ), biomass dry weight (; ), seed length (; ), and seed width (; ) with yield kg/ha. Correlation values for total no. of pod had positive and moderate () highly significant association with yield (; ) while strong () and significant (; ) relation was found with number of mature pod, dry seed weight (; ), and hundred seed weight (; ). Fresh pod weight detected positive and strong () highly significant correlation with yield (; ) along with dry pod weight (; ) whereas dry pod weight showed positive perfect correlation () with yield. The positive and intermediate () highly significant association with yield was noted for hundred seed weight (; ) per plant (Table 8).

4.4. Genetic Parameter Analysis
4.4.1. Variance and Covariance, Heritability in a Broad Sense, Relative Differences, and Genetic Advances

The investigation of covariance (genotypic and phenotypic), the genotypic coefficient of variation (GCV) and phenotypic coefficient of variation (PCV), relative differences (RD), broad-sense heritability, and genetic (GA) advance (as a percentage of mean) were displayed in Table 9. Briefly, the result exhibited that genotypic variance () varied from internode length (0.25) to seed yield (100025). Unvaryingly, for the trait internode length, the least phenotypic variance (), the value was 0.30 whereas for seed yield kg/ha was reported the topmost (100273.24) value. Seemingly, in the case of all traits, the phenotypic variance () is greater than genotypic variance (). Topmost PCV and GCV values were computed for biomass dry weight per plant and biomass fresh weight with a value of 41.09% and 40.53%, respectively. The lowest estimated values for GCV had 6.74% for shelling percent, whereas the lowest PCV values were found 7.26% for the trait days to maturity. The estimated both GCV and PCV values were found more than 20% for the traits such as number of stem per plant (PCV 23.51% and GCV 22.08%), biomass fresh (PCV 40.83% and GCV 40.53%) and dry (PCV 41.09% and GCV 40.09%) weight, fresh (PCV 33.02% and GCV 33%) and dry (PCV 30.18% and GCV 30.14%) pod weight, dry seed weight (PCV 26.31% and GCV 25.68%), and yield (PCV 30.19% and GCV 30.15%) kg/ha which indicate high variability among these traits and for the improvement of accession further selection could be done since the variation of these traits is due to the effect of additive genes. Rest of the traits showed although the traits days to maturity (PCV 7.26% and GCV 6.82%) and Shelling percent (PCV 8.47% and GCV 6.74%) showed below 10% coefficient of variation which indicates the limited scope of selection based on respected traits due to the effect of environment on their phenotypic expression.

(1) Relative Difference (RD). The relative difference (RD) is referred to as an estimation of the ratio of GCV in association with the respective PCV and the estimated RD values varied from 0.05% (fresh pods weight) to 38.67% for no. of branch per plant (Table 9). The traits like no. of branch per plant (38.67%), no. of node per stem (22.55%), no. of immature pod per plant (27.39%), and pod width (32.93%) discovered a more difference in between their PCV and GCV values compared to other existing traits that indicated these traits had wider genetic variability due to environmental effect and not better feedback to direct selection for the improvement of traits. Oppositely, the traits like day to 50% flowering (3.03%), maturity date (6.08%), biomass fresh (0.74%) and dry (2.43%) weight, fresh (0.05%) and dry (0.12%) pod weight, dry seed weight (2.39%), harvest index (2.55%), and yield kg/ha (0.12%) had lower values of relative difference. Based on this result, it was noticed that the variation present among the traits because of a gene which has a better response to direct selection.

(2) Heritability in a Broad Sense (hb2). Heritability refers to the ratio of the total variation of phenotypic traits in each population between the individuals due to genetic variation. Typically, the estimated values of heritability in the broad sense were high () for almost of the all traits evaluated (Table 9). The range of heritability broad sense for the estimated traits was varied from 37.67% (no. of branches per plant) to 99.89% (fresh pod weight). Generally, moderate () heritability values were marked for the traits like branch number per plant (37.61%), no. of immature pods per plant (52.73%), and pod width (44.98%) whereas the rest of the traits expressed high () heritability values, which indicate that the degree of heritability is less affected by the environment.

(3) Genetic Advance (GA). Genetic advance (as percentage mean) was ranged between the lowest (11.06%) for shelling percent and biomass fresh weight (82.86%) (Table 9). Genetic advance for five traits such as days to maturity (13.19%), no. of branches per plant (16%), pod width per plant (14.34%), hundred seed weight (17.99%), and shelling percent (11.06%) showed intermediate genetic advance () whereas rest of the traits detected high () genetic advance values concurrently with high values of heritability. Superior GCV, alongside high heritability as well as high genetic advance, provides a superior indication of selection than the consideration of individual genetic matrix or measuring unit. The advancement of agronomic variables over selection was carried out based on the degree of the genetic variation that exists among the population. Besides this, it also depends on the strength of the transportability of certain traits which is a scale of estimation of heritability.

4.5. Cluster Analysis

For parental selection, genetic distinction analysis is one of the norms, in which the degree of divergence among obtainable landraces is demonstrated. Clustering provides a very strong and strict clue on the degree, and the nature of genetic divergence is notable for the selection of expected genotype. In this study, the homogenized data was used to calculate the Euclidean distances among the 150 Bambara groundnut accessions and the cluster analysis presented as a dendrogram (Figure 5) using UPGMA (average linkage cluster analysis) revealed numerous clusters depicting associations among these collected accessions. To categorize against the relations in the accessions, accessions were clustered into five major clusters (Figure 5) based on their twenty-seven measurable traits at the dissimilarity of 496.7. In the dendrogram, there was a cut off at the point of 496.7 for picking cluster number and ease of interpretation using Mojena’s stopping rules. Cluster I recorded the highest number (39.33%) accessions with an average yield of 20.67% while maximum average yield was recorded for cluster II (29.45%) which consists of 23 accessions with best agronomic traits (Table 10) followed by cluster V (21% accessions) but lower yield of 13.56%. Cluster IV consists of 14% of the accessions with 2nd average maximum yield (23.5%) while cluster III was constructed by only 10% of the total accessions with lowermost average yield (12.48%). Additionally, we recorded 46.51%, 16.90%, and 2.67% higher (+) mean yield compared to average grand mean yield (1049.19 kg/ha) for cluster II, cluster IV, and cluster I, respectively, while cluster III (36.11%) and cluster V (32.53%) gave lower (-) yield. So, cluster II (23) and cluster IV (21) had accessions adjacent to each other on the aspect of yield traits (Figures 6(a) and 6(b)). For this crop improvement, 44 accessions from cluster II and IV associated with large seed size and high yielding potential were marked as potential accessions.

4.6. Estimation of Principal Component Analysis

Principal component analysis (PCA) has been vastly used in crop research for sorting the traits and grouping of accessions. In the current research first, eight principal components (PC) had accounted for 78.99% of the cumulative variation (Table 11 and Figure 7). The 1st PC gained and recorded for the topmost proportion of the variance in the set of all PCs and rest for gradually smaller and smaller amounts of variation. However, the percent of variation for PC1 and PC2 was 34.29% and 11.63%, respectively, while the 9th PC accounted for 3.28% of the variation. The graphical illustration of 150 accessions (a) and 27 morphological traits (b) exposed by PCA was shown in Figure 8. It is intended from Figure 8(b) 8-9 traits (total no. of pods, no. of mature pods, fresh pods weight, dry pod weight, no. of seed/plant, seed length, dry seed weight, hundred seed weight, and yield) and traits (petiole and leaves number, biomass fresh and dry weight, and seed length and width) had a positive correlation with PC1 and PC2, respectively. The traits contributing to PC1 and PC2 are showing the topmost variability with a high coefficient of variation as also exposed in an analysis of variance. The factor loading of several traits is displayed in Table 11 which were revealed by using PCA. The PC1 allowed loading of traits like total no. of pods, no. of mature pods, fresh and dry pods weight, no. of seed/plant, dry seed weight, hundred seed weight, and field yield indicating the significant for the respective principal components simultaneously, for PC2 (petiole and leaves number, biomass fresh and dry weight), for PC3 (seed length and width, harvest index, andplant height), for PC4 (days to emergence, stem number, pod width, and shelling %), PC5 (Ddays to 50% flowering, maturity date, seed length, and width), for PC6 (days to 50% flowering and maturity date), for PC7 (stems and branch number), for PC8 (plant height and internode length), and for PC9 (no. of nodes per stem and branch and stem number), showed significance to the respective PCs. The relationship between eigenvalues and principal component and their proportion of variation are shown in Figure 7. By considering these 8-9 PCs, it was exposed that these PCs regulate the total variation for all yield contributing traits. Furthermore, the two-dimensional graphical elucidation (Figure 9) demonstrated that most of the accessions were dispersed at low distances whereas the few were dispersed at high distances as reflected by an eigenvector (Table 11). The outermost accession from the centroid was S1G24, S1G13, S1G89, S1G54, S1G121, S1G125, S1G39, S1G147, and S1G4 whereas other accessions were near to the centroid.

4.7. Estimation of Shannon–Weaver Diversity ( Index)

The Shannon–Weaver diversity index was used to assess the phenotypic diversity for each trait. The estimation of the Shannon–Weaver diversity index () and evenness () for the twenty-seven traits is shown in Table 11 using the formula [60]. The estimated Shannon–Weaver diversity index ranged from 4.93 to 5.01 among the traits evaluated. The equitability or evenness was found varied from 0.98 to 1.00. Among the traits, topmost () diversity was estimated for maturity date and shelling percent followed by the traits like fifty percent flowering date, nodes number per stem, internode length, pod width, seed length, and seed width which () indicated that maximum diversity was present among these traits while the lowermost diversity () was noted for the trait biomass fresh and dry weight per plant. Similarly, maximum () values of evenness were marked for almost all the traits whereas minimum () was noted for biomass fresh and dry weight per plant.

5. Discussions

5.1. Qualitative Diversity

Effective selection is done when significant hereditary variation present with high magnitude among the accessions has been reported by Hahn [61] and Adebisi et al. [62]. Several researchers such as Mohammed [63], Sinise and Massawe [64], and Abu and Buah [65] noted a significant level of dissimilarity in numerical traits in Bambara groundnut. An observation among Bambara groundnut genotypes with these 3 types of growth habits was noticed by Ntundu et al. [1] in Tanzania and [30] in Cameroon which supported our research findings on qualitative traits. In the study of 52 landraces of Bambara groundnut by Gbaguidi et al. [66], it was concluded that the significant variation is present among all the qualitative traits, taken into consideration in his study. There was a similar observation by Ntundu et al. [1] that the farmers prefer bunch-type landraces compared to other reasons of it provide more advantages to the farmer especially during harvesting period where the roots and stems were unearthed. Vegetative growth of the Bambara groundnut is varied; our research result was equal with those of [67] who grouped Bambara groundnut into three categories namely, (bunche type, semibunch type, and spreading type) based on its vegetative growth. Similar types of growth habits were observed among the germplasms that were evaluated in this research work. For the Bambara groundnut crop, the qualitative trait growth habit is highly significant to the different cropping patterns [68]. The other qualitative traits which showed low variation did not display a rational distinct identity among the Bambara accessions studied.

5.2. Quantitative Traits

The findings of this work displayed that there was a vast genetic variation that exists among the Bambara groundnut accession estimated for twenty-seven numerical traits listed. The similar descriptive statistics analysis (averages, range, coefficient of variation, and standard deviation) confirmed the genetic diversity of the Vigna subterranea (L.) Verdc [1] and the cowpea (Vigna unguiculata L) [69]. Our research revealed that a high variation coefficient for traits like emergence time, biomass fresh and dry weight, branches number, fresh and dry pod weight, dry seed weight,and yield kg/ha. In the African continent, many researchers conducted intensive researches on Bambara groundnut using several morphological traits and their findings were supported by research findings such as genetic diversity and population structure of the Bambara groundnut improvement program [70]. Relevant observation of high variation coefficients was confirmed by Goli et al. [60] and indicates the existence of a massive heterogeneity among the landraces recorded in Cote d’Ivoire. In the variability test between local and exotic Bambara groundnut [71] in Botswana, Bambara groundnut showed a significant morphological variation reported by Bonny and Dje [72] and Touré et al. [73]. In Cameroon, Bambara groundnut accessions were collected from different locations to contrast their traits for improvement by Ndiang et al. [74] and Sobda et al. [75]. Among the 150 accessions, average days to 50% flowering was found close to 38 days in Malaysia which is lower than the 67.65 days reported by Mohammed [63] out of 101 Ghanaian Bambara landraces. Days to 50% flowering was quite diverse among the Bambara groundnut accessions and the estimated values varied from 25 to 51 days, but the variation of flowering time was reported by Goli et al. [60] ranged from 38 to 68 days for 1384 genotypes. Flowering time for the Bambara groundnut is undefined [76] and plays a role as a vital part of the adaption mechanism of a variety to an environment [73]. The flowering times depends on the various complex proceeding of interaction influenced by genetic and/or environmental element [77] and very related trait for annual farming landraces like Vigna subterranea (L.) Verdc [78]. Identical findings were discovered by Massawe et al. [13] ranging from 64 to 76 days in South Africa, while Masindeni [18] reported 43-80 days in Bloemfontein, South Africa. Ouedraogo et al. [10] observed flowering time differs from 32 to 53 days in Burkina Faso. In addition, flowering happened between 36 and 53 days among twenty Bambara groundnut accessions in Pretoria, state of South Africa, observed by Goli et al. [37], also he stated that flowering can be influenced by several environmental factors such as day length, temperature, altitude, and soil conditions as well as genotypic factors. When Bambara groundnut is planted in long-day flowering, it is either delayed or stopped since it is a short-day plant. Early flowering implies early maturity [79], and in our findings, the early flowering accessions S1G108, S1G56, S1G88, S1G75, S1G28, S1G4, S1G13, S1G26, S1G11, S1G29, S1G148, S1G133, S1G49, S1G14, and S1G118 could be selected for early (below 120 days) maturity. Days to maturity was differed significantly () among the accessions and varies from 103 to 163 days, which was more or less similar range reported by Masindeni [18] and Goli et al. [60]. The maturity duration of the Bambara groundnut depends on the cultivar and climatic situation, from 3 to 6 months [79]. Long photoperiods cause delayed maturity of the Bambara groundnut stated by Linneman et al. [80]. In our study, plant height also showed a significant difference (14.23 to 36.9 cm), while all other morphological traits also exhibited variation significantly. Mohammed [63] stated that average values achieved for this trait in his work are supported by this current study. The outcome founded by Ntundu et al. [1] in Tanzania and Goli et al. [37] in South Africa also supported our findings. A strong significant difference () was found between yield and yield contributing traits such as biomass fresh and dry weight, pod and seed number per plant, fresh and dry pod weight, dry seed weight, pod length and width, mature and immature pods per plant, and hundred seed weight showed high genetic variation among these traits. A similar variation in yield contributing traits was also stated by Goli et al. [37] and suggested that these variations due to genotype by environment (G × E) effect on Bambara groundnut yield. In our assessment hundred seed weight, dry pod weight, and total no. of pods varied from 196.21 g to 364.09 g, 88.89 g to 402.9 g, and 38 to 81, respectively. It has been detected that hundred seed weight was considered a vital tool for the judgment of morphophysiological traits related to yield ([1, 10, 13]; [18, 63, 81]). Our calculated yield ranged from 370.38 to 1679 kg/ha, and this finding was supported by Adebisi et al. [62]; his observation was between 146.6 and 2678.6 kg/ha among 52 landraces in Benin; also the study was carried on by [10] in Burkina Faso. Based on our recorded yield, selected 23 elite accessions were identified. These 23 accessions are fitted to cluster II that produced maximum yield and together with 21 accessions of cluster IV created 44 accessions that gave higher yield compared to other accessions evaluated in this study and considered elite accessions. An average 1537.18 kg/ha (29.45%) yield was recorded for 23 accessions of cluster II while the 21 accessions of cluster IV produced 1226.56 kg/ha, which is 23.5% of the total yield. We observed the minimum (5.95%) yield gap between these two groups related to other cluster (I, III, and V). Moreover, cluster II produced 46.51% higher mean yield than the average grand mean yield of 1049.19 kg/ha followed by cluster IV (16.90%). Considering all the parameter studied, a significant relationship was found between the accessions of cluster II and IV in relation to yield and its contributed traits compared to other groups of accessions and this finding were supported by Onwubiko et al. [35]. On the other hand, FAO in 2014 measured the yield of the Bambara groundnut which is lower than the calculated mean yield (1049.19 kg/ha) in our research. A similar observation was noted inW Africa, which was 703.3 kg/ha from the cultivable land of 158,635 hectares [27]. But for the yield of groundnut (Arachis hypogaea L.) which was introduced from South America, the calculated yield was 1058.8 kg/ha from a cultivated area of 6,207,414 hectare in West Africa; for this reason, Bambara groundnut was treated as a neglected and underutilized crop due to its low production in West Africa [82].

5.3. Correlation (Association) Matrix

The correlations assessment indicated a relationship of some morphological traits with the characters of yield. The association or correlation is an influential tool for the researchers to prefer the traits to be integrated into the genotype selection program [63]. The correlation coefficient is an essential measure of an index in plant breeding; after all, it is the measurement of the magnitude of the correlation between genetic and nongenetic two or more variables. Our findings were supported to those achieved by Mohammed [63] in Cote d’Ivoire, [66] observed among 52 landraces in Benin, and [31] in Cameroon. The accessions with large size seeds fulfill the demand of the farmers as well as consumers and often treated as good commercial qualities [83]. Adebisi et al. [62] declared that for the selection of superior accessions the consideration of correlation values among the variables is a great index of the selection process. Total pod number detected positive and moderate highly significant correlation with yield and positively strong correlation was noted with the number of mature pods, dry seed weight, and hundred seed weight; these results are consistent with the findings of [71], variation correlation [84] studies among yield with its related components. Furthermore, previous researches have reported in the notification by Ntundu et al. [1], Ouedraogo et al. [10], Onwubiko et al. [35], and Goli et al. [37] that the significantly positive correlated traits found in this current research were significantly correlated with seed yield in Bambara groundnut. Consistently, these traits can be directly selected for yield improvement. We observed that there was a positive significant association of plant height and other yield-related traits like the number of total pods, no. of mature pods, fresh pod weight, pod length, seed length, dry seed weight, and hundred seed weight with field yield kg/ha, may be suggested that the selection considering these traits may be useful for yield improvement of Bambara groundnut as well as fodder production.

5.4. Genetic Parameter Analysis
5.4.1. Variance and Covariance, Heritability in a Broad Sense, Relative Differences, and Genetic Advances

In past research findings on heritability, it was reported that the selection constructed for certain trait improvement does not only depend on available genetic variation but also the degrees of heritability for such variations [42, 85]. Besides, the valuation of heritability alongside genetic advance contributes a depth advantage over the sole use of heritability [86, 87]. The analysis of the variance elements viz. phenotypic variance and genotypic variance exhibited that phenotypic values were marginally higher than the respected genotypic values for all the traits, are the indication of the trait’s expression are influenced by the environment. These findings of our research were supported by previous reporters [88]. Further, the coefficient of genotypic and phenotypic variation results was evaluated based on the statement of the research of [46, 47, 89]. They suggested that the values of the genotypic coefficient of variation (GCV) and phenotypic coefficient of variation (PCV) categorized for low (0% to 10%), for intermediate (10%-20%), and high (≥20%) variation. Based on these criteria, our research results noted that both genotypic coefficients of variation (GCV) and phenotypic coefficient of variation (PCV) were medium to high for most of the traits. Generally, the selection will be fruitful for the development of traits associated with the degree of desirable variation [90]. In this current research, almost all the traits related to yield exposed medium to strong heritability and genetic advance values except the trait’s branch number per plant, no. of nodes per stem, and no. of immature pod per plant. So, these traits were significantly remarkable for the selection procedure; however, the traits were controlled by the additive genes with limited response to the environment. This result was supported by Meena et al. [91] and Oladosu et al. [92]. Oppositely, the lower level of genetic advance along with low heritability points out the role of nonadditive genes on these traits, which could be possible to enhance over heterosis breeding [93]. Hence, it is meaningful to prefer those traits with an improved genotypic coefficient of variation (GCV), phenotypic coefficient of variation (PCV), heritability, and genetic advance [94]. Only a powerful selection can be achieved when the effects of additive genes are adequately stronger than the effects of the environment [95]. The greater divergence between the genotypic coefficient of variation (GCV) and phenotypic coefficient of variation (PCV) values is the indication of high effects of the environment to a certain trait whereas the smaller divergence is the indication of strong and significant result of accessions on detectable expression including limited effects of the environment [95].

The estimated result of relative difference (RD) was very high for traits like branch number per plant, no. of nodes per stem, pod width, and no. of immature pod per plant. Bello et al. [96] and Umar et al. [42] confirmed similar findings in their research that the variations were present almost due to the effect of the environment since the improvement of traits cannot be attained by direct selection. Reversely, the characters that had minimum relative differences were due to a genetic effect, indicating that the divergence that exists in these traits can be acquired through direct selection [77].

Estimation of heritability and genetic advance play a vital role in assuming the divergence of phenotypic values which were broadly considered as breeding values. Johnson et al. [49] and Assefa et al. [50] graded the heritability measure as between 0 and 30% for low, 30 and 60% for intermediate, and ≥60% for high. Supposedly, most of the traits considered in this work had high heritability values parallel to high genetic advances. Apparently, in our findings, the traits like biomass fresh (, ) and dry (, ) weight, fresh pod weight (, ), dry pod weight (, ), dry seed weight (, ), and yield kg/ha (, ) had high heritability alongside with high genetic advance and suggested greater additive effect of genes which provide effective selection for traits improvement directly. This result has an uninterrupted background by the previous research of [85, 90, 96]. The traits with intermediate heritability values also considered as the influence of environmental effects [97]. Onwubiko et al. [98] and Jonah et al. [99] estimate the genetic parameters in Bambara Groundnut with the similar findings of my research output. Besides, the traits with low heritability and genetic advance indicated that the estimated result because of non-additive provably (dominance and or epistasis) genes and/or effects of environment or combined effects of these dual factors. It has been declared by Cornelius [100] that the trait selection with low and moderate heritability values together with low genetic advance may be delayed in traits improvement till their genetic effects get high on over the effect of the environment [101]. Finally, it is evident from the current research that the improvement of yield and other yield contributing traits of Vigna subterranea (L.) Verdc. can be obtained through selection by the measurement of heritability and genetic advance.

5.5. Cluster Analysis

In the present, investigation of clustering was supported in the past research, observed by Unigwe et al. [43], Gbaguidi et al. [66], Sobda et al. [75], and Bonny et al. [102] in their studied significant variation regarding morphological characteristics distributed in Bambara groundnut, [96] in chili pepper for high yield and CMT values, [83] in cowpea genotypes, and [103, 104] in Capsicum annuum L. genotypes. The cluster analysis based on the UPGMA model using numerical traits constructed four distinct groups of Bambara groundnut genotypes in south Africa reported by Unigwe et al. [43], and Atoyebi et al. [105] also constructed dendrogram using statistical analysis software (SAS version 9.3) among 300 accessions of Bambara groundnut.

5.6. Estimation of Principal Component Analysis

Typically, the principal component analysis (PCA) is the rejustification tool of cluster analysis. Genetically, identical accessions were clustered into the same group stated by Falconer [50] also genetically dissimilar parents can cover a high degree of heterosis. Johnson [106] noted that principal component analysis intends to determine the total variation that exists in a set of traits which sequentially accounts for the maximum variability in the data. Generally, traits are inter-correlated to varying level and hence all the principal components are not required to summarize the data effectively. The first axes (PC1) elucidate the utmost portion of the total variation in any PCA [107]. Our observation was supported by several types of the research reported by Bello et al. [96], Farhad et al. [108], and Maqbool et al. [109]. Shegro et al. [37] grouped the 20 Bambara groundnut accessions by PCA analysis based on quantitative traits. In our finding’s variation percentages of PC1 and PC2 are 34.29% and 11.63% while [110] identified that PC1 and PC2 highly donated to the total variation at 19% and 14%, respectively, in Bambara groundnut. To cluster the genotypes into groups and subgroups, principal component scores were used because first a few principal components controlled all the information of the original variables [111]. Daudo and Olakojo [112] found the similar output during working on maize genotypes, Mustafa et al. [113] observed comparable findings and decided that selection of characters with greater eigenvalues controlled the diversity among the accessions. Mustafa et al. [114] and Jolliffe [115] also gave attention to using accessions based on component traits.

5.7. Estimation of Shannon–Weaver Diversity ( Index)

In our study, the observed diversity index value was more than 4.93 for most of the traits evaluated. This finding supported by Aliyu et al. [60] reported index varied from 1.60 to 2.07 for twenty quantitative traits of Bambara lines, Olukolu et al. [116] reported index of nineteen qualitative traits of Bambara groundnut varied from 0.1 to 0.15 and twenty-eight numerical traits of 124 accession of Bambara groundnut which showed index values between 0.09 and 0.16 across the four African regions. Bonny et al. [102] evaluated the diversity in qualitative traits of Bambara groundnut landraces (Vigna subterranea L verdc.) in Cốte d’Ivoire of similar findings with me. Nonetheless, the values of index for traits appeared statically more or less similar, suggesting a similar genetic diversity. The report from Alvarez et al. [117], Robert et al. [118], and Thomas et al. [119] showed that Bambara groundnut is a self-pollinated crop; therefore, the diversity level of this crop is influenced by farmers’ agricultural practices as well as seed management techniques such as recycling, storing, exchanging and newly introducing of species.

6. Conclusion

It is noticeable from this current research that the enhancement of yield and other yield-related traits of Bambara groundnut (Vigna subterranea L. Verdc) can be obtained through selection by the determination of different genetic parameters analysis. Additionally, the degree of divergence recorded for almost all the agromorphic variables was studied. However, it can be beneficial to the advancement of agromorphic traits of Bambara groundnut by the plant breeders. The current research also resolved strongly to the perfect association between the morphological traits and the field yield. The yield-related traits like no. of the stem, no. of the petiole, no. of mature pods, biomass fresh and dry weight, fresh and dry pod weight, pod weight, dry seed weight, and yield kg/ha recorded high GCV and PCV values were ≥20% with high genetic advance also a low relative difference. It is evident from this present study that the enhancement of yield and other yield-related traits can be attained through effective selection based on estimates of heritability and genetic advance. The result from the principal component and cluster analyses depict that 23 accessions from cluster II and 21 accessions from cluster IV considered high-yielding accessions and can be suggested as large-sized seed associated with high yield potentials. These 44 elite accessions among 150 of Bambara groundnut lines were suggested to grow for further evaluation via conventional breeding alongside with molecular study for confirming and identifying the best 20 accessions.

Data Availability

All data are provided in full in the results section of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

Authors’ Contributions

Md Mahmudul Hasan Khan and Mohd Y. Rafii did the conceptualization. Md Mahmudul Hasan Khan and Mashitah Jusoh did the data curation. Md Mahmudul Hasan Khan, Mohd Y. Rafii, Shairul Izan Ramlee, and Mashitah Jusoh did the investigation. Md Mahmudul Hasan Khan and Mashitah Jusoh acquired software. Mohd Y. Rafii supervised the study. Md Mahmudul Hasan Khan wrote the original draft. Md Mahmudul Hasan Khan, Mohd Y. Rafii, and Md Al Mamun wrote and reviewed and edited the manuscript.

Acknowledgments

The authors are grateful to the Ministry of Agriculture (MoA), Bangladesh Agricultural Research Council (BARC- Project of NATP Phase-II), Bangladesh Agricultural Research Institute (BARI) of the People's Republic of Bangladesh and Universiti Putra Malaysia (research grant: vote number 6282518), for adequate funding. Thankfulness also goes to Dr. Yussuf Oladosu Adeniyi, postdoctoral fellow, ITAFoS, UPM, Malaysia, for his constructive criticism.