Genetic analysis and selection of Bambara groundnut (Vigna subterranea [L.] Verdc.) landraces for high yield revealed by qualitative and quantitative traits

As a crop for the new millennium Bambara groundnut (Vigna subterranea [L.] Verdc.) considered as leading legumes in the tropical regions due to its versatile advantages. The main intent of this study was to find out the high yielding potential genotypes and considering these genotypes to develop pure lines for commercial cultivation in Malaysia. Considering the 14 qualitative and 27 quantitative traits of fifteen landraces the variation and genetic parameters namely, variability, heritability, genetic advance, characters association, and cluster matrix were determined. ANOVA revealed significant variation for all the agronomic traits (except plant height). Among the accessions, highly significant differences (P ≤ 0.01) were found for almost all the traits excluding fifty percent flowering date, seed length, seed width. The 16 traits out of the 27 quantitative traits had a coefficient of variation (CV) ≥ 20%. A positive and intermediate to perfect highly significant association (r = 0.23 to 1.00; P < 0.00) was found between yield and its related traits. The trait dry seed weight per plant (g) had the highest GCV = 59.91% and PCV = 59.57% whereas the trait fresh pod weight (99.55%), dry seed weight (98.86%), and yield (98.10%) were highly heritable. The genetic advance recorded the highest for dry seed weight (122.01%) and lowest (3.97%) for plant height. To validate the genetic disparity, an unweighted pair-group produce with arithmetic mean (UPGMA), principal component analysis (PCA), heatmap, and H’-index was performed considering 27 quantitative traits. The constructed dendrogram showed five distinct groups of accessions. Genotypes G2, G3, and G9 from Group IV consider as promising lines which gave 70.05% higher mean yield compared to grand mean yield (1180 kg ha−1) with desirable traits. Group II had a maximum number of accessions while group III and group V had one of each. However, findings declared that the availability of genetic variance will be beneficial for this crop improvement through cross breeding as well as the plant breeders to prefer desirable traits in V. subterranea L. Verdc. for further breeding purposes.

ground soil, hence jointly its common name is 'Bambara groundnut' 7 . Due to rich in carbohydrate (63-65%), protein (18-20%) and oil (17-18%), Bambara groundnut has been defined as a fully well-adjusted food for human feeding 8 . According to 9 it contains 32.72% essential amino acids and 66.10% non-essential amino acids. The seed is regarded as a balanced food because when compared to other food legumes, it is rich in iron and its protein contains high level of lysine and methionine 10 . Lysine is the major essential amino acid and represents 10.3% of the total essential amino acid 10 . Bambara groundnut can fulfill the regular demand of protein for the marginal users where the animal protein is badly available due to its high cost 11 . Seeds of Bambara groundnut contain considerable amount of minerals such as Ca: 260 mg; K: 1723.25 mg; Fe: 3.6 mg and Na: 75.25 mg of each 100 g dry (weight) seed, moreover due to rich in potassium (K) it also have ability to lessen diabetes by prompting the insulin hormone 12 . The undeveloped fresh seeds of Bambara groundnut can be boiled to make pudding (Moi-Moi or Okpa) 13 ; fodders to feed animals 14 and extract of leaves has medicinal values as anti-vomiting agent 15 . In many developing countries where cultivation of other major crops is difficult, but Bambara groundnut can be accommodated nicely due to its drought tolerant and low diseases-insects infestation nature 16 and as legume it can fix atmospheric nitrogen via nodulation 17 . This crop has the capacity to give high yield with low input and mostly grown by female in sole culture without any modern techniques 18 and 10-40% of their total yield they sold in market rest is used by themselves 19 . According to FAOSTAT 20 , the annual production is about 42,023 kg ha −1 , of which Africa produces half, with Burkina Faso occupied the major producing country. Bambara groundnut can well adapt to the tropical area like Malaysia, where cultivation of major crops (Rice, Wheat, Maize, etc.) are increasingly challenging due to drought and unpredictable rainfall patterns 21 . The average yield of Bambara groundnut is low 0.650-0.850 t ha −122 , 28.96 g plant −112 , and 0.38-1.6 t ha −123 . The production of Bambara groundnut is mainly limited due to lack of improved cultural techniques 24 and improvement of Bambara groundnut was neglected for many years by researchers because of the lack of available fund and unprivileged effort on its improvement 25 . Although numerous efforts have unsuccessful for the varietal improvement of Bambara groundnut through hybridization due to cleistogamous and autogamous nature of flower and knowledge gap on its reproductive biology as under developing crop 26 . On the aspect of Malaysia, the breeding approach of Bambara groundnut is undetermined, and no commercial high yielding variety is available, so the requirement is to discover commercially high yielding cultivar for certain rising areas 27 . Germplasm screening considering the agronomic variables is the initial attempt to identify the targeted characters of interest 28 . This current research reveals the genetic divergence of fifteen Bambara groundnut accessions to discover the existing variation and the selection to develop high yielding pure lines for this crop improvement. Accordingly, this study provides an evidence on Bambara groundnut diversity among the landraces that introduced from Africa (Nigeria) to Asia (Malaysia). All the modern applicable techniques may be applied for the betterment of this ongoing cultivated crops, but the dual approaches like conventional breeding linked with molecular breeding is highly successive over the solely use of one approach. However, traits improvement can be possible through direct selection with valuation of different genetic parameter analysis. The extent of selection approach exceedingly inspire by heritability and genetic gain estimation is the commanding tools for the enhancement of a certain traits 29 . Hence, the core intent of this study was to determine the inherent variation of Bamabara groundnut landraces using both qualitative and quantitative traits via valuation of characters association, variance component and different genetic parameters, resulting the identification of high yielding potentials from which pure lines will be developed for commercial cultivation. The soil PH is 6.6 to 7.5 with sandy loam to clay loam type (Dept. of land management, UPM). Fifteen accessions of Bambara groundnut were selected for this current research work, all representing the African accessions collected from the local market of Nigeria. Land races of Bambara groundnut used in this research was listed in Table 1. Randomly five plants were taken into consideration to evaluate genetic variability based on the agronomic traits 30 .

Materials and methods
Experimental design. The experiment was conducted in a randomized complete block design (RCBD) with three replications. The experimental plot comprised of two rows measuring 1.6 m × 0.80 m. The distance between plant to the plant 30 cm, row to row 50 cm, plot to plot 1.5 m and the distance between replication was 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; P ≤ 0.05 and to compare among the mean of significant of traits. The correlations between the quantitative variables were determined using Pearson 33 correlation coefficient formula. The genotypic and phenotypic variation were calculated as per following the formula given by Singh and Choudhary 34 . The coefficient variation of phenotypic (PCV) and genotypic (GCV): were estimated as per formula given by Khan et al. 23 also relative differences was estimated using the formula (RD) = Relative difference between PCV and GCV. The estimated values of PCV and GCV were categorized by Robinson et al. 35 and Khan et al. 23 , like as between 0 and 10% for low, 10-20% for intermediate and greater than (≥ 20%) for high. Broad sense heritability ( h 2 b ) was estimated using the formula given by Falconer 36 and Khan et al. 23 . In accordance with Johnson et al. 37 and Khan et al. 23 , the heritability was categorized as between 0 and 30% for low, 30-60% for intermediate and greater than 60% as high. Genetic Advance (GA) (as a percentage of mean): was calculated with 5% selection intensity (K) following the method of 37 . Genetic advance is categorized as between 0 and 10% for low, 10% to 20% for intermediate and more (> 20%) than for high, following the formula given by Khan et al. 23 . K for constant also indicates the intensity of selection.

Results
Assessment of qualitative variation. The frequency of distribution of some qualitative variables are summarized in Fig. 1 Table 3. We found 16 out of the 27 quantitative traits had coefficient of variation (CV) ≥ 20% which ranged from 5.96% (Shelling%) to 58.55% (dry seed weight per plant (g). The average days to maturity was found (131.56 ± 1.18) days which is statistically significant (p ≤ 0.01). The highest value of standard deviation (SD) was found for yield kg/ha (650.98) with standard error (SEm ± 97.04) while the lowest was for internode length (SD = 0.49; SEm ± 0.07) ( Table 3). Standard error (SE) is the indication of consistency of the average values, lowest SE values indicate the sample mean is more precise reflection of the real population mean. The maximum and minimum values for overall accessions and mean comparison with least significant difference (LSD = 0.05) were shown in the Tables 4 and 5. Fig. 3 showing the graphical relationship of DSW(g) and HSW(g) with yield (kg ha −1 ). The days to 50% flowering varied from 31 to 49 days after sowing (DAS) while 66.67% of the accessions gave 50% flower before 40 days after sowing. Most of the landraces (80%) took more than 120 days to maturity which varied from 122 to 141 DAS. The genotype G1 marked as short duration line with a total day to maturity of 119 DAS ( Table 4). The genotype G2 recorded highest values for the traits like-BFW (614.67 g); BDW(361.57 g) ( Analysis of correlation matrix. The phenotypic correlation among the 27 numerical traits of fifteen Bambara groundnut accessions is given in Table 6. Days to 50% flowering showed negative and intermediate (0.25 ≤ r < 0.75) significant association with NB ( r = −0.29; P = 0.04), NNS(r = 0.33; P = 0.02), TNP (r = −0.38; P = 0.00), NMP (r = −0.35; P = 0.01) and NSP (r = −0.38; P = 0.00). A positive and highly strong (0.75 ≤ r < 1) significant association was found for the traits like NMP (r = 0.93; P ≤ 0.00), DPW (r = 0.76; P ≤ 0.00), NSP (r = 0.99; P ≤ 0.00), DSW (r = 0.77; P ≤ 0.00) and yield (r = 0.76; P ≤ 0.00) with the total number of pods. The trait yield (kg ha −1 ) revealed positive and perfect (r = 1.00) highly significant association with DPW (r = 1.00; P ≤ 0.00) while positive and highly strong ( 0.75 ≤ r < 1) significant association was found with DSW (r = 0.99; P ≤ 0.00), NSP (r = 0.76; P ≤ 0.00) and HSW (r = 0.76; P ≤ 0.00) per plant ( Table 6). Correlation networking and pattern search plot are more visual illustration of correlation matrix. In correlation network (Fig. 4a) each node represents a variable with its colour based on the defined geographical population (Gombe-7 accessions, Kwami-2 accessions, Akko-2 accessions, Alkaleri-2 accessions, and Sokoto-2 accessions) group and its size is based on number of correlations to that variable. The traits yield, DPW, FPW, HSW, NP, and NL showed larger node size then other traits. Two variables are connected by an edge if the correlation between the two variables meet the p-value (0.05) and thresholds. The edge size also reflects the magnitude of the correlation. Helps in identifying biologically meaningful relationship or associations between group or features. On other hand, pattern search ( Fig. 4b)  Abundance (richness) analysis for traits. Abundance sketching (Stacked Bar) provides overall as well as comparative abundance profile across landraces at different variable levels (Fig. 5a). Variation present in data of multiple variable level can be visualized for all individual genotypes wise. It also summarized and compared the abundance of different variables based on the multiple data. The pie chart provides an exact clue for all variables at multiple variation level present in data. Pie chart helps in visualizing the degree of diversity or abundance of variables compositions for different BG samples. Provides exact composition of each group through direct quantitative comparison of abundances in percentages (Fig. 5b). In the "Stacked bar or Area plot" each bar represents an individual genotype under the respective location (at the top of area olot). Among the traits, yield, NL, and BFW possess the higher area with showing different colour for each genotype however, according to pic chart these traits also occupied 28%, 21% and 7%, respectively.
Analysis of genetic components. Variance and covariance, heritability in broad sense, relative differences, and genetic advances. The output of genetic components analysis was compiled in Table 7. Apparently, the   .29%) and shell% (PCV 6.05% & GCV 4.50%) showed below 20% of phenotypic coefficient of variation (PCV) and genotypic coefficient of variation (GCV). For the improvement of this crop further selection could be done considering the traits having GCV ≥ 20% (BFW, BDW, FPW, DPW, DSW and grain yield) which indicated high degree of variability among these traits although the variation is due to the effect of additive genes. Due to the lower GCV values (≤ 10%) the vegetative traits (PH, D50%F, and DTM) and yield traits (SL, SW, and Shel%) indicated the limited chance of selection based on respective traits due to the effect of environment on their phenotypic expression.
The relative difference (RD) is the ratio of GCV in association with the respective PCV and the estimated RD values varied from 0.22% (fresh pods weight) to 57.46% for plant height (Table 7). Relatively low difference value between GCV and PCV was recorded for the traits like DTM(9.74%), NB(2.59%), BFW(1.17%), BDW(5.18%), TNP(4.91%), FPW(0.22%), DPW(0.95%), PL(3.15%), DSW(0.57%), HSW(4.71%), and yield (0.95%) kg ha −1 and noticed that the variation present among the traits due to the effect of gene and have a better response to direct selection. On the other hand, the traits with higher difference in between their PCV and GCV values indicated the wider genetic variability due to environmental effect and not better feedback to direct selection for the improvement of traits.
Clustering patterns. In this study, the homogenized data was used to calculate the Euclidean distances among the 15 Bambara groundnut accessions and a UPGMA dendrogram was designed (Fig. 6). To discriminate against the relations in the population, the dendrograms of the 15 Bambara groundnut accessions were clustered into five major groups based on their twenty-seven measurable traits at 1.16 dissimilarity coefficients. In the dendrogram, there was a cut off at the point of 1.16 coefficient for ease of interpretation. The Table 8 showed the mean performances of the selections according to each class. Group I represented by G1. G6 and G12 are characterized by early germinating (6 days) but need more time to flowering close to 44 days among the group also took medium time to maturity and plant height (24.45 cm). Internode length was maximum while hundred seed weight and yield kg/ha showed medium values within the five groups. Group II is formed by maximum number of (G4, G5, G7, G11, G14 and G15) accessions characterized by the following traits: maximum maturation date (137 days) and minimum pod length ( (Table 8). Moreover, we estimated 70.07% higher ( +) mean yield over the average grand mean yield (1180 kg ha −1 ) for cluster IV whereas the cluster I (9.34%), cluster II (31.44%), cluster III (50.08%) and cluster V ( 12.95%) produced lower (−) yield.

Heatmap analysis for genotypes and agro-morphological traits.
A heatmap is a data imagining practice that displays extent of a phenomenon as color in two dimensions. The variation in color may be by hue or intensity, giving noticeable visual indications to the reader about how the phenomenon is clustered or varies over space. It visualizes the relative patterns of high-abundance features against a background of features that are mostly lowabundance or absent. Heatmap analysis of agro-morphological descriptors were carried out to show a chromatic evaluation of the Bambara groundnut genotypes. The heatmap analysis constructed double dendrograms, the 1st dendrogram on the vertical direction, an arrangement that represent the Bambara groundnut accessions, and the 2nd dendrogram on the horizontal direction representing traits that influenced this diffusion. Dendrogram 1 showed two major group, group (a) linked to two genotypes Maikai (G2) and Giiwa (G9) while the group (b) comprised of rest 13 genotypes (Fig. 7). Dendrogram 2 also displayed two major groups, group  The two-dimensional (Fig. 8a) and three dimensional (Fig. 8b) graphical elucidation demonstrated that most of the accessions were dispersed at low distances whereas the few were dispersed at high distances as reflected by eigenvector ( Table 9). The farthest accession from the centroid was G2, G3, G8, G9, G12 and G13 whereas other accessions were near to centroid. The proportion of variation for principal component (PC1) and (PC2) were 45.88% and 10.64% respectively, in which the first principal component occupied the topmost position of the total variation existed (Table 9). PCA biplot loaded the both variables and cases (accessions) at the same time shows how strongly each trait influences a PC and correlated to each other it also shows the how distances the genotype from each other. The lesser angle between two vectors (Fig. 9) indicate higher and positive correlation (e.g. DSW & Yield), when angle between two vectors form 90° indicate no correlation while it goes more than 90° to near 180°, indicate negative correlation between the traits (e.g. HSW & NL) .

Valuation of Shannon-Weaver diversity (H' Index). The Shannon-Weaver diversity index was used
to assess the phenotypic diversity for each trait. The estimation of the Shannon-Weaver diversity index (H) and Evenness (E H ) for the twenty-seven traits shown in

Selection of Bambara groundnut accessions based on Rank summation index (RSI). Pearson's
correlation coefficients were estimated for all traits to define the traits of a positive or negative association with yield. The traits which were strongly positively associated with yield were used to compute Rank Summation Index (RSI) to lead high yielding Bambara groundnut selection from the population. Furthermore, correlations within a couple of traits are probably significant when the ideal values of the coefficient (%) are higher than 0.20 43 . Using Mulumba and Mock 44 reported formula of RSI method the accessions were first ranked using average values of the respective traits of positively correlated with yield (here, 1 = topmost and 15 = lowermost), consequently the ranked values of traits were summed to estimate the total performance of individual accessions. In this way, the accession with the lowest RSI score indicates the superior yielding potential. The values of RSI ranking of Bambara groundnut with the traits of positive and strong significantly correlated traits with yield namely fresh pod weight, dry pod weight, no. of seed per plant and dry seed weight per plant are represented in Table 10. Based on yield, the topmost genotype G2 (2991 kg ha −1 ), G9 (2226 kg ha −1 ), G3 (1557 kg ha −1 ), G6 (1414 kg ha −1 ), G10 (1250 kg ha −1 ), G1(1101 kg ha −1 ) were identified as high yielding accessions among the

Discussion
Qualitative disparity. The existence of a significant qualitative variation was found for all the qualitative traits, supported by Gbaguidi et al. 45 he found significant variation among all the qualitative traits. We recorded three types of growth habit and similar observation were noticed by Ntundu et al. 1 in Tanzania, Khan et al. 23 in Malaysia and Azam-Ali et al. 46 in Cameroon. We categorized the vegetative growth of Bambara groundnut namely, bunches type, semi bunches type and spreading type which is matched by the result of Doku 47 and highly significant difference among the qualitative trait was noted by Egbadzor et al. 48 .
Quantitative traits. The estimated 27 quantitative traits showed a massive genetic variation and similar variation was confirmed by Ntundu et al. 1 and Aliyu et al. 49 in Vigna subterranea (L.) Verdc and the cowpea (Vigna unguiculata L) 50 . The estimated high coefficients of variation (CV) in our study is the indication of vast scale of heterogeneity confirmed by Goli et al. 51 in Bambara groundnut. We found D50%F close to 39 days but 52 noticed close to 68 days in Ghana. The indeterminate 53 nature of flower bearing make it vital issue for adjustment mechanism to an environment 54  Several climatic issues photoperiod, temperature, altitude and soil structure as well as genotypic nature is responsible to bearing flower in Bambara groundnut 28 and reported flowering happened between 36 to 53 days. In our study, genotype G1, G2, G3, G8, G9 and G10 identified as early flowering lines; early flowering ensures early maturity 56 . A significant difference (P ≤ 0.01) was recorded for maturity (119.67 to 141.33) days is supported by Goli et al. 51 and Masindeni 55 and due to diverse cultivar along with multi-environmental factors maturation time varied from 90 to180 days 57 . Plant hight had no significant variation, supported by Ntundu et al. 1 in Tanzania and Shegro et al. 28 in south Africa. The yield and yield related traits like TNP, NSP, FPW, DPW, DSW, PL, PW, NMP, NIP and HSW showed high genetic discrepancy, similar variation stated by Shegro et al. 28 with a recommendation of variation happened due to effect of genotype by environment interaction Bambara yield. Hundred seed weight varied from 177.52 g to 360.15 g, is a vital factor for the measurement of morphological traits linked to yield 23,52,55,58 it also influences the yield directly. The yield of Bambara groundnut was recorded from 146.6 to 2678.6 kg ha −1 by Gbaguidi et al. 45  Correlation coefficient. In plant breeding correlation matrix is a prominent approach for the judgement of degree of the association between two or more variables, is supported by Mohammed 52 . For superior genotype's selection programme consideration of correlation matrix can be a great scale of measurement 60 . Strong and positive significant correlation for total number of pods (TNP) was identified with the traits NMP, DPW, NSP, DSW and Yield this result is consistent with the study of Pranesh et al. 61 and Jonah et al. 62 . We got moderate and positively high significant association of plant hight (PH) with TNP, NMP, FPW and yield can be proposed the selection based on these traits may be beneficial for yield enhancement of this crop as well as fodder production for animal feeding. Similar recommendation was stated by Mohammed 52 in Cote d'Ivoire and 63 in Cameroon. Nankar et al. 64 described correlation network in tomato phenotypic diversity and pattern search plot that was also supportive to my findings. Relative abundance study using area plot and pie chart revealed richness of traits diversity for each genotype. The traits possessing higher extent in both area plot and pie chart indicates that these traits are highly contributed to governing the diversity among the landraces. www.nature.com/scientificreports/ Genetic components. For the selection program variation presents among the traits was taken into consideration which depends on the degree of heritability. To know the projected gain from selection, valuation of genetic advance with heritability can be a significant approach of crop improvement. Various research findings reported that the selection may be effective for a specific trait improvement using available genetic variation with the degree of heritability 29,65 . Consideration of both heritability and genetic advance is more effective over the uniquely use of heritability 66,67 . Like the previous reporters Adebola et al. 68 findings, we disclose higher phenotypic variance values than genotypic variance for all traits, indicates the trait expression govern by the environment. The obtained GCV and PCV value was categorized based on the suggested index of 0%-10% for low, 10-20% for moderate and ≥ 20% for high variation 23,69,70 . Intermediate to strong genetic advance with heritability was found for all yield related traits except seed width and shelling% is the indication of the traits have significant potential in the selection process due to low environmental influences, supported by Meena et al. 71 .
The improvement of the traits with low heritability and genetic advance can be boost over heterosis breeding this is supported by Usman et al. 70 . The value of relative differences between GCV and PCV had higher for the trait plant hight, seed width, days to 50% flowering, and seed length is the sign of higher environmental effect and the improvement of these traits are tough via direct selection whereas the trait with lower difference is the symbol of lower influence by the environment which may give desirable strong and significant output in crop improvement program, is supported by Umar et al. 29 and Usman et al. 70 . Direct selection can be effective considering the traits Table 8.   65,73 . Low to moderate heritability and genetic advance values may hindrance in the trait's betterment due to high environmental effects over the genetic effects on its stated by Ridzuan et al. 74 . So, only an effective selection can be gained picking the traits with higher GCV, PCV, HB, and GA meaning that effect of additive genes is sufficiently robust than environmental effect 70 .
Clustering patterns. Five clusters were constructed based on the 27 quantitative traits at 1.16 of the distant coefficients that indicates a degree of diversity among the genotypes. The cluster V considered as potential group of genotypes for the crop betterment associated with high yielding capacity. The findings of previous researchers 30,45,75,76 stated that they constructed same type of cluster and found significant variation regarding morphological traits in Bambara groundnut. The study of Unigwe et al. 30 explored the four distinct groups of Bambara groundnut genotypes in south Africa using UPGMA model. The timing of flowering duration is a motivational factor for the final yield also play a positive role to the best yield of the group and selection could be effective from this class noted by Tourél et al. 77 . Flowering in Bambara groundnut is indeterminate up harvesting stage explained by Kumaga et al. 78 . However, early flowering has been considered as a well agronomic trait of crops to quick maturity, uniform yield as well as generally crop production 78 thus, accessions that have early flowering criteria should be treated as best to production of Bambara groundnut 79 . The groups achieved from the cluster analysis of quantitative characteristics illustrate the performances of Bambara groundnut accession cultivated in Benin would be the future guideline for this crop improvement 22,80 . The clustering and characterization of accessions considering their agro-morphological traits and genetic similarity would be the crucial issue to identification and selection of the best parents for hybridisation 81 . Additionally, cluster IV produced 70.05% higher mean yield than the average grand mean yield of 1180 kg ha −1 while the other groups gave lower yield and this finding were supported by Onwubiko et al. 73 . Therefore, current research represents significant information to the plant breeders based on their similarity and grouping of accessions through univariate and multivariate methods. The heatmap analysis depicted the depth of correspondence among morphological traits evaluated of Bambara groundnut genotypes and this result was constantly supported by Virga 28 grouped the 20 Bambara groundnut accessions by PCA analysis using quantitative traits. For yield improvement the selection PC1 was revealed as the most powerful criterion concluded by the work of Adéoti et al. 87 and Mih et al. 88 . In my research total pod numbers, mature pods number, seed number, dry seed weight and yield kg/ha occupied high values in PC1. This finding supported by Stoilova & Pereira 80 described that the most significant components for yield are the pods number and seeds number per plant. The cluster analysis together with principal component analysis explored the common association among landraces in terms of seed yield and related agronomic traits. www.nature.com/scientificreports/

Shannon diversity index (H) and evenness (E). Shannon's diversity index (H) 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. The estimated H' Index varied from 2.57 for dry seed weight per plant to 2.71 for plant height, maturity date and shelling percent among the phenotypic traits. In our study the observed diversity index value was more than 2.50 for most of the traits evaluated and highest the value indicates higher diversity, though H' index ranges typically from 1.5 to 3.5 but rare can be reaches 4.5 89,90 . Olukolu et al. 22 reported H' Index of nineteen qualitative traits (0.1 to 0.15) and twenty-eight numerical traits (0.09 to 0.16) of Bambara groundnut that supported our findings. Bonny et al. 76 evaluated the diversity in qualitative traits of Bambara groundnut landraces of similar findings with our result.
Selection of high yielding accession using RSI method. The result from correlation studies revealed the traits that have a positive significant correlation with yield also give valuable information in a breeding program for the selection of high yielding Bambara groundnut. Similar findings were concluded by Onwubiko et al. 43 and Ajala et al. 91 , he got two high yielding genotypes from 33 accessions using the RSI method. Other

Conclusion
From the present study an evident has been established that the improvement of Bambara groundnut (Vigna subterranea [L.] Verdc.) yield and it related traits can be gained via selection with the valuation of different genetic parameters analysis like GCV, PCV, HB, and GA. Based on the recorded data and considering the supplementary analysis (heatmap study, correlation network, abundance analysis) it can be state that a considerable degree of variation exist in almost all the agronomic traits evaluated in this study. Moderate to perfect significant association was noted between the yield and its related traits. Additionally, this research also depicts selection criteria using the traits which had strongly positive correlation with grain yield. More than 20% PCV and GCV values was estimated for all traits excluding the traits like TNP, PW, NSP, SL, SW, and shell% beside this, the six traits like DSW, DPW, FPW, BDW, BFW, and Grain yield showed high (≥ 20%) genetic advance (as percentage mean) with high heritable values. However, it can be declared that a higher extent of divergence was detected among tested landraces based on H'-index (1.57-2.71) as well as Euclidian distance clustering (into five group). Considering all statistical findings, the genotype G2, G3, G8 and G9 identified as high yielding promising lines and can be use as distance parents for hybridization program. We suggested that further research can be conducted to gain the homogeneity of genotypes based on yield and its contributed traits improvement. Concurrently, we must provide emphasis on intensive research of these potentially high yielding lines, together with conventional breeding and molecular approaches.