Hostname: page-component-586b7cd67f-2plfb Total loading time: 0 Render date: 2024-11-24T17:39:18.560Z Has data issue: false hasContentIssue false

Genetic variation in Bambara ground nuts as revealed by agro-morphological and DArTseq markers and selection for improved yield performance

Published online by Cambridge University Press:  01 August 2023

Michael M. Chipeta*
Affiliation:
Lilongwe University of Agriculture and Natural Resources, Bunda College, P.O BOX 219, Lilongwe, Malawi
Davis Gimode
Affiliation:
International Crops Research Institute for the Semi-Arid Tropics- Kenya, PO BOX 39063-00623, Nairobi, Kenya
*
Corresponding author: Michael M. Chipeta; Email: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

Bambara groundnut (Vigna subterranea (L.) Verdc) has been neglected in terms of variety selection and development which has resulted in farmers growing a mixture of landraces that are not genetically characterized and are low yielding. With the need to set up a breeding programme in Malawi, it was necessary to thoroughly understand the genetic diversity (GD) present in the available germplasm. The objectives of the study were to assess Bambara genotypes GD using agro-morphological traits and SNP markers, and to identify and select high yielding Bambara genotypes. Field trials were conducted for two seasons at Bunda College. Later, genotypes were genotyped using DartSeqLD SNP markers. All data were analysed using R Package. Significant genetic variations (P < 0.001) were observed for most traits including grain yield, which suggests that genetic variability exists in Bambara groundnuts which can be exploited in breeding programmes aimed at developing high performing varieties. Based on grain yield, the study identified 18 top performing genotypes across the evaluation seasons which will be tested under farmers’ fields’ conditions. DArTseqLD grouped the genotypes into three clusters. It was noted that majority of the genotypes from the same origin clustered together. High genetic distances were observed between genotypes from Southern African and West African regions and this has important implications in parental selection for the genetic improvement of Bambara. Our results provide valuable insights about the extent of genetic variability and how parental lines can be selected for improved genetic gain in Bambara groundnuts.

Type
Research Article
Copyright
Copyright © The Author(s), 2023. Published by Cambridge University Press on behalf of National Institute of Agricultural Botany

Introduction

Bambara groundnut (Vigna subterranea (L.) Verdc) is an underutilized and research neglected leguminous crop (Massawe et al., Reference Massawe, Mwale, Azam-Ali, Roberts, Ochatt and Jain2007) and in Malawi, it is grown in most areas especially by rural women. It is regarded as a complete food as it contains sufficient quantities of protein (18–20%), carbohydrate (60–65%) and fat (Goli, Reference Goli, Begemann and Mushonga1995; National Research Council, 2006). Therefore, Bambara groundnut is one of the legumes with the potential of improving nutrition and addressing malnutrition levels reportedly at 47% among the children under the age of five in Malawi.

The world holds large collections of Bambara germplasm characterized by wide genetic diversity (GD). Several studies of GD analysis in Bambara groundnut have been conducted using both agro-morphological traits (Karikari and Tabona, Reference Karikari and Tabona2004; Ntundu et al., Reference Ntundu, Shillah, Marandu and Christiansen2006; Ouedraogo et al., Reference Ouedraogo, Ouedraogo, Tignere, Bilma, Dabire and Konate2008; Abu and Buah, Reference Abu and Buah2011; Olukolu et al., Reference Olukolu, Mayes, Stadler, Ng, Fawole, Dominique, Azam-Ali, Abbott and Kole2012; Aliyu and Massawe, Reference Aliyu and Massawe2013; Molosiwa et al., Reference Molosiwa, Aliyu, Stadler, Mayes, Massawe, Kilian and Mayes2015; Olanrewaju et al., Reference Olanrewaju, Oyatomi, Babalola and Abberton2021) and molecular markers (Olukolu et al., Reference Olukolu, Mayes, Stadler, Ng, Fawole, Dominique, Azam-Ali, Abbott and Kole2012; Aliyu and Massawe, Reference Aliyu and Massawe2013; Mayes et al., Reference Mayes, Basu, Molosiwa, Redjek, Ahmad, Khan, Zehra, Noah, Mayes, Roberts, Stadler, Massawe, Kilian and Azam-Ali2013; Molosiwa et al., Reference Molosiwa, Aliyu, Stadler, Mayes, Massawe, Kilian and Mayes2015). The molecular markers that have been used in GD analysis in Bambara include; RAPD (Amadou et al., Reference Amadou, Bebeli and Kaltsikes2001; Massawe et al., Reference Massawe, Roberts and Davey2003; Rungnoi et al., Reference Rungnoi, Suwanprasert, Somta and Srinives2012), AFLP (Massawe et al., Reference Massawe, Dickinson, Roberts and Azam-Ali2002; Ntundu et al., Reference Ntundu, Bach, Christiansen and Andersen2004), SSR (Somta et al., Reference Somta, Chankaew, Rungnoi and Srinives2011; Olukolu et al., Reference Olukolu, Mayes, Stadler, Ng, Fawole, Dominique, Azam-Ali, Abbott and Kole2012; Aliyu and Massawe, Reference Aliyu and Massawe2013; Molosiwa et al., Reference Molosiwa, Aliyu, Stadler, Mayes, Massawe, Kilian and Mayes2015), DArT Array markers (Stadler, Reference Stadler2009; Olukolu et al., Reference Olukolu, Mayes, Stadler, Ng, Fawole, Dominique, Azam-Ali, Abbott and Kole2012; Molosiwa et al., Reference Molosiwa, Aliyu, Stadler, Mayes, Massawe, Kilian and Mayes2015; Ahmad et al., Reference Ahmad, Chai, Basu, Sri Redjeki, Moreton, Mayes, Ho, Massawe and Mayes2016).

GD analysis helps to align the breeding objectives with identified diversity more especially during selection of parental materials to initiate a breeding programme. No comprehensive agro-morphological or molecular marker GD analysis has been reported in Malawi on Bambara groundnuts except for one genotype that was reported by Molosiwa et al. (Reference Molosiwa, Aliyu, Stadler, Mayes, Massawe, Kilian and Mayes2015). This suggests that Bambara groundnut has been grossly neglected in Malawi in terms of research and variety selection. To date, Malawi lacks an active Bambara groundnut breeding programme. As such there is a dearth of improved varieties for farmers, resulting in the cultivation of low yielding landraces with an average yield of less than 400 kg/ha among smallholder farmers (Pungulani et al., Reference Pungulani, Kadyampakeni, Nsapato and Kachapila2012). However, this is lower than the Africa average of 0.85 t/ha. Yields well over 3 t/ha have been reported in Africa (Azam-Ali et al., Reference Azam-Ali, Sesay, Karikari, Massawe, Aguilar-Manjarrez, Bannayan and Hampson2001; Hillocks et al., Reference Hillocks, Bennett and Mponda2012; Berchie et al., Reference Berchie, Dapaah, Agyeman, Sarkodie-Addo, Addo, Addy and Blankson2016; Mubaiwa et al., Reference Mubaiwa, Fogliano, Chidewe, Jan Bakker and Linnemann2018; Tan et al., Reference Tan, Azam-Ali, Goh, Mustafa, Chai, Ho, Mayes, Mabhaudhi, Azam-Ali and Massawe2020). This suggests there is potential to develop high yielding Bambara groundnuts for improved food and nutrition security.

With renewed interest in setting up a Bambara groundnut improvement programme for Malawi, there was need for a thorough understanding of the GD present in the existing germplasm using both morphological characterization and genetic markers. Therefore, this research was formulated to characterize Bambara genotypes using agro morphological traits and DartSeqLd SNP markers. Further, to identify and select high yielding Bambara groundnuts lines that can be released as new varieties to the farming communities or used as parents in a Bambara groundnut breeding programme.

Materials and methods

Study site description

The study was conducted at Bunda College of Agriculture Research farm in the 2019/2020 and 2021/2022 seasons. The site receives a mean annual rainfall of about 930 mm with a temperature range of 17.2°C to 19.1°C.

Plant materials

A diverse panel of 56 Bambara groundnut germplasm was assembled. This included materials collected locally from farmers and markets in Malawi as well as materials from different countries sourced through the IITA Gene bank Nigeria.

Field trials experimental designs and data collection

The field agronomic trials were carried out under rain-fed conditions in 2019/2020 and 2021/2022 seasons. For the 2019/2020 trials, 56 genotypes were planted using a 7 by 8 alpha lattice design in which there were 7 blocks of size 8 and replicated three times. The plot constituted of 2 ridges each measuring 1.5 m with an intra-row spacing of 15 cm and 0.75 m between ridges. One seed was planted per planting station and the trial was planted on 30 November 2019. For the 2021/2022 trial, 60 genotypes were planted on 3rd January, 2022 using a 5 by 12 alpha lattice design. The trial had 3 replications with 12 blocks of size 5. The plot had 2 ridges spaced at 0.75 m, each 2.1 m long. Data were collected following the procedure described in the Bambara groundnut descriptor (Ipgri/Iita/Bamnet, 2000).

Genotyping sample preparation, processing and analysis

For genotyping purposes, six plants for each of the 60 genotypes were tagged in preparation for leaf sampling. The six plants within the genotypes were sampled to detect the within genotype variation to facilitate single plant selection.

Leaf tissue sampling for DNA analysis

Three weeks after germination, a young leaf from each of the six tagged plants from each genotyped was sampled for DNA analysis. The sampling was done according to the required procedure by Intertek-Agritech laboratory (Intertek-Agritech, 2016). First, unique sample identifications were generated using breeding management system, which were used to track the samples throughout the genotyping process. A single hole-puncher (6.0 mm diameter) was used to punch and collect four leaf discs per sample from young and healthy newly developed leaf. The punched leaf discs from each sample were transferred with forceps directly into a single well of a 96-well plate (1.2 ml AbGene Storage Plate, number AB0564, Thermo Fisher Scientific). Wells G12 and H12 of the sampling plates were left blank as control wells according to Intertek-Agritech lab standard protocol. The forceps were cleaned with 75% alcohol before and after placing each sample in a well to avoid cross contamination. The sample plates were later oven dried for 24 h at 48 degrees Celsius. The 96-well sample plates were sealed with sealing mats (AB0674, Thermo Fisher Scientific), wrapped in plastic bags, secured firmly and shipped for DNA extraction to Intertek-Agritech laboratory in Sweden.

DNA isolation and genotyping

After the extraction of High Molecular Weight DNA from the tissue samples and after quality check, DNA was sent for high-throughput DArTseqLD SNP genotyping at the Diversity Arrays Technology Pty Ltd. Canberra, Australia. The DNA samples were subjected to complexity reduction by digestion/ligation reactions as described by Kilian et al. (Reference Kilian, Wenzl, Huttner, Carling, Xia, Blois, Caig, Heller-Uszynska, Jaccoud, Hopper, Aschenbrenner-Kilian, Evers, Peng, Cayla, Hok and Uszynski2012) but replacing the single PstI-compatible adaptor with two adaptors corresponding to PstI and MseI restriction enzymes’ overhangs. The PstI-compatible adaptors were designed to include the Illumina flowcell attachment sequence, barcode region, while the reverse adaptor contained a flowcell attachment region and MseI-compatible overhang sequence. Only mixed fragments (PstI-MseI) were effectively amplified in PCR. Subsequently, equimolar amounts of PCR product from each sample were bulked and sequenced by the Hiseq2500/ Novaseq6000 (Illumina® Inc., San Diego, CA, USA). Thereafter, the generated sequences were processed using proprietary DArT analytical pipelines. Approximately 250,000 sequences per barcode/sample were used in marker call. Identical sequences were collapsed into FASTQCOL files followed by SNPs calling using the software package DArTsoft14.

Data analysis

Agro-morphological data

An exploratory analysis using density plots was done on quantitative data to check for ANOVA assumptions which included Shapiro–Wilk test for normality and Levine's test for equal variance. Thereafter, all analyses were done in the R software environment (R Core Team, 2021). ANOVA was done by fitting the model with the lmer function of the lme4 package (Bates et al., Reference Bates, Mächler, Bolker and Walker2015) in R and P values were adjusted using Tukey method to compare means at 0.05 significance level. Principal component analysis (PCA) was done using the FactoMineR package (Lê et al., Reference Lê, Josse and Husson2008) and Pearson correlation was performed using the corr. function of the stats package in R. A hierarchical cluster analysis using Euclidean distance metric and the associations between groups were done using the ward D2 method in a cluster R package (Maechler et al., Reference Maechler, Rousseeuw, Struyf, Hubert and Hornik2019).

Analysis of Bambara groundnut genotyping data

Genotyping data were filtered using a minor allele frequency (MAF) threshold of 0.05 which eliminated 660 SNPs. Diversity metrics which included MAF, Neis GD, polymorphic information content and inbreeding coefficient (Fi) were estimated using the snpReady package (Granato et al., Reference Granato, Galli, de Oliveira Couto, e Souza, Mendonça and Fritsche-Neto2018). Principal coordinate analysis was also conducted using dartR package in R (Gruber et al., Reference Gruber, Unmack, Berry and Georges2018) to determine the contributions of each component to the variation that existed in the germplasm. AMOVA was performed to quantify genetic variation among individuals, among populations and within individuals using Genalex ver. 6.502 (Peakall and Smouse, Reference Peakall and Smouse2012), with statistical significance based on 999 permutations. Bayesian information criterion was used to estimate the optimum number of clusters (K), followed by examination of population structure by discriminant analysis of principal components (DAPC) using the adegenet package in R (Jombart et al., Reference Jombart, Kamvar, Collins, Lustrik, Beugin, Knaus and Jombart2020).

Results

Genotypic variation in Bambara groundnuts

Tables 1 and 2 show the mean square values for the different traits scored in 2019/2020 and 2021/2022 seasons, respectively. For the 2019/2020 season, there were significant genotype effects for grain yield (P < 0.001), pod length (P < 0.001), pod width (P < 0.001), number of pods per plant (P < 0.001) and number of seeds per pod (P < 0.001). For the 2021/2022 season, a similar trend was observed where significant differences were observed for grain yield (P < 0.01), 100 seed weight (P < 0.05), seed length (P < 0.001), number of pods per plant (P < 0.001), pod width (P < 0.001), pod length (P < 0.001), plant height (P < 0.001), canopy spread (P < 0.001), number of stems per plant (P < 0.001), number of branches per stem (P < 0.001), number of nodes per stem (P < 0.001), internode length (P < 0.001), petiole length (P < 0.001), terminal leaflet width (P < 0.001), terminal leaflet length (P < 0.001), number of leaves per plant and peduncle length (P < 0.001). However, significant genotypic effects, were not observed for number of seeds per pod (P > 0.05), number of days to 50% flowering (P > 0.05) and number of days from sowing to first flowering (P > 0.05).

Table 1. Mean square values for traits scored in 2019/2020 season

SOV, Source of variation; DF, Degrees of freedom; GYLD, Grain yield (kg/ha); PODL, Pod length; PODW, Pod width. ***highly significant at P < 0.001, **highly significant at P < 0.01, *Significant at P < 0.05. P values adjusted using Tukey method at 0.95 Confidence level.

Table 2. Mean square values for traits scored in 2021/2022 season

GYLD, Grain yield (kg/ha); S100W, Weight of 100 seeds; SDL, Seed length; NSP, Number of seeds per pod; NPP, Number of pods per plant; PDW, Pod width; PDL, Pod length; D50F, Number of days from first flowering to 50% flowering; PHT, Plant height; CNS, Canopy spread; DFF, Number of days from sowing to first flowering; NSTP, Number of stems per plant; NBS, Number of branches per stem; NNS, Number of nodes per stem; INL, Internode length; PTL, Petiole length; TLW, Terminal leaflet width; TLL, Terminal leaflet length; NLP, Number of leaves per plant; PEDL, Peduncle length. ***highly significant at P < 0.001, **highly significant at P < 0.01, *Significant at P < 0.05. P values adjusted using Tukey method at 0.95 Confidence level.

Performance of Bambara groundnuts for yield and related traits

The results on the mean performance of Bambara genotypes for the 2019/2020 growing season are presented in Table S1. Grain yield ranged from 252 kg/ha (TVSu-1831) to 1410 kg/ha (TVSu-1321). For pod length, values ranged from 15.4 mm to 21.3 mm for TVSu-1831 and TVSu-368, respectively. In terms of pod width, genotypes TVSu-1448 and TVSu-1733 had the highest pod width of 11.3 mm and 11.0 mm, respectively and genotypes TVSu-1258 (7.6 mm) and TVSu-1831 (6.7 mm) recorded lowest pod width. On average, the Bambara genotypes produced 18.4 pods per plant. For the 2021/2022 growing season, (Tables S2 and S3), grain yield ranged from 104 kg/ha (TVSu-250) to 1828 kg/ha (Mazira-a-mpheta). In terms of weight of 100 seeds, the highest weights were obtained from genotypes TVSu-1557 (90.1 g) followed by BC13 (85.1 g). The lowest 100 seed weight was recorded on genotype TVSu-1888 (46.5 g).

Principal component and correlation analysis for agro-morphological traits

PCAs grouped all the agronomic traits scored in 2021/2022 season under five components (Eigen values ≥1) which accounted for 71.8% of the variation (Table S4). All the traits exhibited a negative association with PC1 which accounted for 32.2% of the total variation. Number of branches per stem, number of stems per plant and number of leaves per plant were correlated with PC2 and PC3, while grain yield, pod length, pod width and peduncle length were correlated with PC4.

The correlation analysis exhibited several significant correlations among the traits scored (Fig. 1). The distribution of each trait is shown on the diagonal while on the bottom of the diagonal, the bivariate scatter plots with a fitted line are displayed. On the top of the diagonal, the correlation coefficients and the significance level (P-values) are indicated. Grain yield (kg/ha) was significantly correlated with 100 seed weight (r = 0.32, P < 0.05), number of pods per plant (r = 0.53, P < 0.001), plant height (r = 0.41, P < 0.01), canopy spread (r = 0.58, P< 0.001), petiole length (r = 0.41, P< 0.01), number of leaves per plant (r = 0.44, P< 0.001), terminal leaflet width (r = 0.35, P< 0.01) and peduncle length (r = 0.29, P< 0.05). Apart from grain yield, several other traits were also highly correlated showing a positive association between the traits.

Figure 1. A scatter plot matrix showing the correlation coefficients among traits and the significance levels scored in 2021/2022 season.

Hierarchical clustering of genotypes based on agro-morphological traits

Cluster analysis grouped the genotypes into five clusters based on the agro-morphological traits (Fig. 2). Cluster V had the highest number of genotypes (17) and lowest number was observed in cluster I (8) (Table S5). In terms of origin of the genotypes, clusters II and IV were dominated with genotypes from Malawi while cluster I was dominated by genotypes from Cameroon and cluster III was dominated by genotypes from Burkina Faso and Malawi. Cluster V had a unique distribution in which genotypes origin varied remarkably but it was dominated by genotypes from West Africa. Clusters I and II were dominated by high yielding genotypes, high number of pods per plant, and most of these originated from Malawi. Clustering analysis showed genotypes from the same geographical areas were generally grouped together. A similar trend was observed with clustering based on SNP markers.

Figure 2. Hierarchical clustering dendrogram analysis; Euclidean distance was used and the associations between groups were done by the Ward method for the morphological traits.

Molecular diversity and population structure

The average Nei genetic distances (Table 3) among the genotypes within each population revealed that the genotypes in population 1 (mostly West and some East African genotypes) were closely related to population 3 (Southern African genotypes) which exhibited the smallest genetic distance value (0.030). Though both population 1 and 2 were dominated by West African genotypes, they exhibited the greatest genetic distance (0.137). Similarly, population 3 (Southern African genotypes) were distantly (0.115) related to population 2 (West African genotypes). The populations 3 and 1 were closely related probably due to the presence of East African genotypes (Sudan and Kenya) which geographically are in close proximity to the Southern African region, or they might have shared many common alleles and are closely related to each other.

Table 3. Pairwise Nei genetic distances among three subpopulations

Pop 1, Cluster 1 = > West African genotypes (including some East African Genotypes from Kenya & Sudan); Pop 2, Cluster 2 = > West African genotypes (Tanzania is the outlier, and UK); Pop 3, Cluster 3 = > Mostly Southern African genotypes.

The AMOVA results (Table 4) showed that the total genetic variance found among populations was 10% while among the individuals contributed to 90% of the genetic variance and no variation was detected within the individuals. The high percentage value of GD within-population (among individuals) obtained from AMOVA could be due to natural adaptation.

Table 4. Analysis of molecular variance (AMOVA) for the three geographical regions of the Bambara groundnut Genotypes

The cumulative percentage contributions of the first two principal components (PCs) to the total variations in the populations were 33%, with 22.6% and 10.4% respective contributions to PC1 and PC2 (Figure S1). The membership clustering of DAPC revealed three major groups of genotypes (Figure S2). DArTseqLD showed that 358 samples derived from 60 genotypes of Bambara groundnut were grouped into three clusters (Figure S3). Cluster I had 60 samples, cluster II had 110 samples and cluster III had 188 samples.

Discussion

The significant variations observed in the genotypes for a number of traits suggest that genetic variability exists in Bambara groundnuts which could be exploited in breeding programmes aimed at developing high performing varieties. Significant genotypic variability has extensively been reported (Ntundu et al., Reference Ntundu, Shillah, Marandu and Christiansen2006; Vurayai et al., Reference Vurayai, Emongor and Moseki2011; Shegro et al., Reference Shegro, Jansen Van Rensburg and Adebola2013; Mohammed, Reference Mohammed2014; Aliyu et al., Reference Aliyu, Massawe and Mayes2016; Gbaguidi et al., Reference Gbaguidi, Dansi, Dossou-Aminon, Gbemavo, Orobiyi, Sanoussi and Yedomonhan2018; Valombola et al., Reference Valombola, Akundabweni, Awala and Hove2019; Olanrewaju et al., Reference Olanrewaju, Oyatomi, Babalola and Abberton2021; Onwubiko, Reference Onwubiko2021) on a number of traits in Bambara groundnut including grain yield, 100 seed weight, pod length, pod width and many others. The existence of this genetic variability signifies that Bambara groundnut could be responsive to selection and that superior genotypes could be developed despite reported difficulties in making crosses (Basu et al., Reference Basu, Roberts, Azam-Ali, Mayes and Kole2007) due to strict self-pollination in the crop.

The average grain yield reported here is 714 kg/ha which is above the national average that smallholder farmers realize in their farms in Malawi (Pungulani et al., Reference Pungulani, Kadyampakeni, Nsapato and Kachapila2012). The differential performance could be attributed to different management conditions under farmers’ fields and heterogeneity of the farm lands. This can also be attributed to better genetic potential of the genotypes used in this study than landraces used by farmers which tend to be inherently low yielding compounded by poor management practices. Other studies have also reported grain yield values within this margin, for example, Akpalu et al. (Reference Akpalu, Sarkodie-Addo and Akpalu2012) reported Bambara groundnut grain yield in the average range of 650 kg/ha–850 kg/ha. The results on grain yield from this study are consistent with what others have reported as yield ranges in Bambara groundnuts (Akpalu et al., Reference Akpalu, Sarkodie-Addo and Akpalu2012; Pungulani et al., Reference Pungulani, Kadyampakeni, Nsapato and Kachapila2012; Mohammed, Reference Mohammed2014; Gbaguidi et al., Reference Gbaguidi, Dansi, Dossou-Aminon, Gbemavo, Orobiyi, Sanoussi and Yedomonhan2018; Valombola et al., Reference Valombola, Akundabweni, Awala and Hove2019; Kunene, Reference Kunene2021; Olanrewaju et al., Reference Olanrewaju, Oyatomi, Babalola and Abberton2021).

The 100 seed weight is a useful indicator for yield and this study has also shown some significant positive correlation between the two traits. Others have used 100 seed weight in evaluating morphological features such as seed size (Massawe et al., Reference Massawe, Mwale and Robert2005; Mohammed, Reference Mohammed2014; Unigwe et al., Reference Unigwe, Gerrano, Adebola, Pillay and Monrovia2016) signifying its role in determining grain yield. The positive association between grain yield and other traits suggested that indirect selection for yield trait can be achieved to an extent through selection based on other traits especially in the early stages of selection where multi location trials may not be possible. This also suggests that breeding for any one of these correlated traits will not reduce the desired level of the other and therefore, selection of high yielding genotypes could be based on the other traits such as 100 seed weight and number of pods per plant. These findings auger well with other studies (Ofori, Reference Ofori1996; Aliyu et al., Reference Aliyu, Massawe and Mayes2016; Gbaguidi et al., Reference Gbaguidi, Dansi, Dossou-Aminon, Gbemavo, Orobiyi, Sanoussi and Yedomonhan2018; Valombola et al., Reference Valombola, Akundabweni, Awala and Hove2019; Olanrewaju et al., Reference Olanrewaju, Oyatomi, Babalola and Abberton2021) which had also reported significant correlations in Bambara groundnuts for different traits.

GD is an important enabler for crop improvement. The current study observed higher GD within population (90%) than among populations, suggesting substancial gene flow among the landraces from the different origins, across large geographical distances. Though some literature have reported higher genetic variation among populations than within population, this study is consistent with what Minnaar-Ontong et al. (Reference Minnaar-Ontong, Gerrano and Labuschagne2021) and Uba et al. (Reference Uba, Oselebe, Tesfaye and Abtew2021) observed in their respective studies.

Genotypes in the same cluster reflect genetic relationships and most diverse accessions were out grouped. The understanding of GD within Bambara genotypes will enhance future planning in breeding programmes and provide helpful information in maintaining and monitoring GD required for a robust breeding programme (Eltaher et al., Reference Eltaher, Sallam, Belamkar, Emara, Nower, Salem K, Poland and Baenziger2018).

Clustering analysis showed that some of the Bambara groundnut genotypes from the same region did not cluster accurately based on their origin. According to Uba et al. (Reference Uba, Oselebe, Tesfaye and Abtew2021) this could be due to low genetic differentiation among the populations, which suggests that the genetic background of Bambara groundnut populations does not always correlate with their geographical area of origin. However, it was worth noting that majority of the genotypes from the same origin clustered together which provides evidence in the origin and relationship of Bambara groundnut genotypes. Similar patterns of grouping of landraces according to geographical origin have been reported using morphological markers (Ntundu et al., Reference Ntundu, Shillah, Marandu and Christiansen2006; Aliyu and Massawe, Reference Aliyu and Massawe2013) and molecular markers (Jaccoud et al., Reference Jaccoud, Peng, Feinstein and Kilian2001; Rungnoi et al., Reference Rungnoi, Suwanprasert, Somta and Srinives2012; Odongo et al., Reference Odongo, Oyoo, Wasike, Owuche, Karanja and Korir2015; Uba et al., Reference Uba, Oselebe, Tesfaye and Abtew2021) for collections of Bambara groundnut landraces within countries and among regions. The clustering of some genotypes from the different regions into the same cluster might suggest some relatedness between accessions from different regions which could partly be attributed to seed exchange as previously alluded to by Uba et al. (Reference Uba, Oselebe, Tesfaye and Abtew2021). Population structure analysis of this study suggested that the genotypes with unknown origin and United Kingdom genotypes might have originated from Western Africa region, because they clustered consistently with the subpopulation consisting of genotypes from West Africa.

Other studies (Somta et al., Reference Somta, Chankaew, Rungnoi and Srinives2011; Rungnoi et al., Reference Rungnoi, Suwanprasert, Somta and Srinives2012; Uba et al., Reference Uba, Oselebe, Tesfaye and Abtew2021) had reported similar trend where they observed that genotypes from locations close to each other were more related and had small genetic distances between them. High distant genetic relatedness reported here has important implications in parental selection for the genetic improvement of Bambara. If hybridization is made between distant related genotypes, it may result in high genetic recombination in the progenies that might be useful in enhancing the genetic gain from selection in the crop.

Conclusion

The significant differences observed on various morphological traits among Bambara groundnut genotypes such as grain yield and weight of 100 seeds confirm that there is genetic variation within and between Bambara groundnut genotypes. While some Bambara groundnut genotypes can be identified and selected and grown directly for their high yields, some may be identified, selected and exploited for utilization in Bambara groundnut breeding programmes. DArT SNP marker molecular characterization of Bambara groundnut genotypes from different regions revealed that variation exists among the genotypes. This level of GD could be the basis for developing new Bambara cultivars with desirable characteristics such as high yield while being adapted to diverse environments. The study identified three subpopulations which could be explained by their geographical differentiation, natural selection and regional adaptation history. Based on grain yield, the study identified the following genotypes as top performers across the evaluation seasons TVSu-1321, TVSu-1557, TVSu-1256, TVSu-1733, Mazira-a-mpheta, BC7, Mazira-a-mpheta-1, TVSu1794, TVSu-1448, TVSu-793, TVSu-388, TVSu-1312, TVSu-1797, BC13, TVSu-542, Kayera, BC3 and TVSu-1045. These will need to be further tested under farmers’ field conditions to confirm the high yielding ability and most preferred ones and later on released for production by farmers.

Supplementary material

The supplementary material for this article can be found at https://doi.org/10.1017/S1479262123000503.

Acknowledgements

Research assistants at Bunda College of Agriculture provided technical support for the trials. International Foundation for Science provided the financial support for the study.

Author contributions

Michael Chipeta conceptualized, designed and implemented the research, collected data, analysed phenotyping data and wrote the first draft of the manuscript. Davis Gimode analysed molecular data and revised the draft. Both authors read and approved the final manuscript.

Financial support

This research was supported by the International Foundation for Science, Stockholm, Sweden, through a grant number 1-1-c-6190-1 to Michael Chipeta.

Competing interest

The authors have no competing interests to declare that are relevant to the content of this article.

References

Abu, HB and Buah, JSS (2011) Characterization of Bambara groundnut landraces and their evaluation by farmers in the Upper West Region of Ghana. Journal of Developments in Sustainable Agriculture 6, 6474.Google Scholar
Ahmad, NS, Chai, HH, Basu, S, Sri Redjeki, E, Moreton, J, Mayes, K, Ho, WK, Massawe, F and Mayes, S (2016) Exploring the domestication of Bambara groundnut. In Proceedings of the XXIX International Horticultural Congress on Horticulture: Sustaining Lives, Livelihoods and Landscapes (IHC2014): IV International Symposium on Plant Genetic Resources, Brisbane, Australia, 17–22 August 2014; pp. 183–190.Google Scholar
Akpalu, MM, Sarkodie-Addo, J and Akpalu, SE (2012) Effect of spacing on growth and yield of five Bambara groundnut (Vigna subterranea (L) Verdc.) landraces. Journal of Science and Technology 32, 919.Google Scholar
Aliyu, S and Massawe, FJ (2013) Microsatellites based marker molecular analysis of Ghanaian Bambara groundnut (Vigna subterranea (L.) Verdc.) landraces alongside morphological characterization. Genetic Resources and Crop Evolution 60, 777787.Google Scholar
Aliyu, S, Massawe, F and Mayes, S (2016) Genetic diversity and population structure of Bambara groundnut (Vigna subterranea (L.) Verdc.): synopsis of the past two decades of analysis and implications for crop improvement programmes. Genetic Resources and Crop Evolution 63, 925943.CrossRefGoogle Scholar
Amadou, HI, Bebeli, PJ and Kaltsikes, PJ (2001) Genetic diversity in Bambara groundnut (Vigna subterranea L.) germplasm revealed by RAPD markers. Genome 44, 995999.CrossRefGoogle ScholarPubMed
Azam-Ali, S, Sesay, A, Karikari, S, Massawe, F, Aguilar-Manjarrez, J, Bannayan, M and Hampson, K (2001) Assessing the potential of an underutilized crop – a case study using Bambara groundnut. Experimental Agriculture 37, 433472.CrossRefGoogle Scholar
Basu, S, Roberts, JA, Azam-Ali, SN and Mayes, S (2007) Bambara groundnut. In Kole, C (ed.), Pulses, Sugar and Tuber Crops. Genome Mapping and Molecular Breeding in Plants. Berlin, Heidelberg: Springer, vol. 3, pp. 159173.CrossRefGoogle Scholar
Bates, D, Mächler, M, Bolker, B and Walker, S (2015) Fitting linear mixed-effects models using lme4. Journal of Statistical Software 67, 148.CrossRefGoogle Scholar
Berchie, JN, Dapaah, HA, Agyeman, A, Sarkodie-Addo, J, Addo, JK, Addy, S and Blankson, E (2016) Performance of five Bambara groundnut (Vigna subterranea (L.) Verdc.) landraces in the transition agroecology of Ghana under different Sowing Dates. Agricultural and Food Science Journal of Ghana 9, 718729.Google Scholar
Eltaher, S, Sallam, A, Belamkar, V, Emara, HA, Nower, AA, Salem K, FM, Poland, J and Baenziger, PS (2018) Genetic diversity and population structure of F3:6 nebraska winter wheat genotypes using genotyping-by-sequencing. Frontiers in Genetics 9, 76.CrossRefGoogle Scholar
Gbaguidi, AA, Dansi, A, Dossou-Aminon, I, Gbemavo, DSJC, Orobiyi, A, Sanoussi, F and Yedomonhan, H (2018) Agromorphological diversity of local Bambara groundnut (Vigna subterranea (L.) Verdc.) collected in Benin. Genetic Resources and Crop Evolution 65, 11591171.CrossRefGoogle Scholar
Goli, AE (1995) Characterization and evaluation of IITA's Bambara groundnut collection. In Begemann, JH and Mushonga, J (eds), Proceedings of the Workshop on Conservation and Improvement of Bambara Groundnut (Vigna subtarranea (L.) Verdc.). Zimbabwe: International Plant Genetic Resources Institute (IPGRI), pp. 410.Google Scholar
Granato, IS, Galli, G, de Oliveira Couto, EG, e Souza, MB, Mendonça, LF and Fritsche-Neto, R (2018) snpReady: a tool to assist breeders in genomic analysis. Molecular Breeding 38, 17.CrossRefGoogle Scholar
Gruber, B, Unmack, PJ, Berry, OF and Georges, A (2018) dartR: an R package to facilitate analysis of SNP data generated from reduced representation genome sequencing. Molecular Ecology Resources 18, 691699.CrossRefGoogle Scholar
Hillocks, RJ, Bennett, C and Mponda, OM (2012) Bambara nut: a review of utilisation, market potential and crop improvement. African Crop Science Journal 20, 116.Google Scholar
Intertek-Agritech (2016) Agri-Services CGIAR HTPG PROJECT Sampling Instructions for SNP verification and routine SNP Analysis. Laboratory Manual published by Excellence in Breeding (EiB) platform.Google Scholar
IPGRI/IITA/BAMNET (2000) Descriptors for Bambara groundnut (Vigna subterranea). International Plant Genetic Resources Institute, Rome, Italy; International Institute of Tropical Agriculture, Ibadan, Nigeria; The International Bambara groundnut Network, Germany. ISBN 9290434619.Google Scholar
Jaccoud, D, Peng, K, Feinstein, D and Kilian, A (2001) Diversity Arrays: a solid state technology for sequence information independent genotyping. Nucleic Acids Research 29, 17.CrossRefGoogle ScholarPubMed
Jombart, T, Kamvar, ZN, Collins, C, Lustrik, R, Beugin, M, Knaus, BJ and Jombart, MT (2020) Adegenet: exploratory analysis of genetic and genomic data. The Comprehensive R Archive Network; https://cran.r-project.org/web/packages/adegenet/.Google Scholar
Karikari, SK and Tabona, TT (2004) Constitutive traits and selective indices of Bambara groundnut (Vigna subterranea (L.) Verdc.) landraces for drought tolerance under Botswana conditions. Physics and Chemistry of the Earth 29, 10291034.CrossRefGoogle Scholar
Kilian, A, Wenzl, P, Huttner, E, Carling, J, Xia, L, Blois, H, Caig, V, Heller-Uszynska, K, Jaccoud, D, Hopper, C, Aschenbrenner-Kilian, M, Evers, M, Peng, K, Cayla, C, Hok, P and Uszynski, G (2012) Diversity arrays technology: a generic genome profiling technology on open platforms. Methods in Molecular Biology 888, 6789.CrossRefGoogle ScholarPubMed
Kunene, SP (2021) Evaluation of Bambara groundnut (Vigna subterranea (L) Verdc.) accessions for drought tolerant and yield performance using agro-morphological and physiological traits (PhD thesis). University of KwaZulu-Natal, Pietermaritzburg, South Africa.Google Scholar
, S, Josse, J and Husson, F (2008) FactoMineR: an R package for multivariate analysis. Journal of Statistical Software 25, 118.CrossRefGoogle Scholar
Maechler, M, Rousseeuw, P, Struyf, A, Hubert, M and Hornik, K (2019) cluster: Cluster Analysis Basics and Extensions. R Package Version 2.1.0.Google Scholar
Massawe, FJ, Dickinson, M, Roberts, JA and Azam-Ali, SN (2002) Genetic diversity in Bambara groundnut (Vigna subterranean (L.) Verdc.) landraces revealed by AFLP markers. Genome 45, 11751180.CrossRefGoogle Scholar
Massawe, FJ, Roberts, JA and Davey, MR (2003) Genetic diversity in Bambara groundnut (Vigna subterranea (L.) Verdc.) Landraces assessed by Random Amplified Polymorphic DNA (RAPD) markers. Genetic Resources and Crop Evolution 50, 737741.CrossRefGoogle Scholar
Massawe, FJ, Mwale, SS and Robert, JA (2005) Breeding in Bambara groundnut (Vigna subterranea (L.) Verdc.): strategic considerations. African Journal of Biotechnology 4, 463471.Google Scholar
Massawe, FJ, Mwale, SS, Azam-Ali, SN and Roberts, JA (2007) Towards genetic improvement of Bambara groundnut [Vigna subterranea (L.) Verdc.]. In Ochatt, S and Jain, SM (eds), Breeding of Neglected and Under-Utilized Crops: Spices and Herbs. Science Publishers, p. 468.Google Scholar
Mayes, S, Basu, SM, Molosiwa, O, Redjek, ES, Ahmad, NS, Khan, F, Zehra, S, Noah, S, Mayes, K, Roberts, J, Stadler, F, Massawe, F, Kilian, A and Azam-Ali, S (2013) Molecular analysis of Bambara Groundnut, an underutilized African Legume Crop as Part of the BAMLINK Project—What Lessons Can We Learn? In: Massawe F, Mayes S, Alderson P (eds) 2nd international symposium on underutilized species.CrossRefGoogle Scholar
Minnaar-Ontong, A, Gerrano, AS and Labuschagne, MT (2021) Assessment of genetic diversity and structure of Bambara groundnut [Vigna subterranea (L.) verdc.] landraces in South Africa. Science Reports 11, 7408.CrossRefGoogle ScholarPubMed
Mohammed, MS (2014) Pre-breeding of Bambara groundnut (Vigna subterranea (L.) Verdc.) (PhD thesis). University of KwaZulu-Natal, Pietermaritzburg, South Africa.Google Scholar
Molosiwa, OO, Aliyu, S, Stadler, F, Mayes, K, Massawe, F, Kilian, A and Mayes, S (2015) SSR marker development, genetic diversity and population structure analysis of Bambara groundnut [Vigna subterranea (L.) Verdc.] landraces. Genetics Resources and Crop Evolution 62, 12251243.CrossRefGoogle Scholar
Mubaiwa, J, Fogliano, V, Chidewe, C, Jan Bakker, E and Linnemann, AR (2018) Utilization of Bambara groundnut (Vigna subterranea (L.) Verdc.) for sustainable food and nutrition security in semi-arid regions of Zimbabwe. PLoS One 13, e0204817tt.CrossRefGoogle ScholarPubMed
National Research Council (2006) Lost Crops of Africa: Vegetables, Vol. II, Washington DC: National Academies Press.Google Scholar
Ntundu, WH, Bach, IC, Christiansen, JL and Andersen, SB (2004) Analysis of genetic diversity in bambara groundnut [Vigna subterranea (L.) Verdc.] landraces using amplified fragment length polymorphism (AFLP) markers. African Journal of Biotechnology 3, 220225.Google Scholar
Ntundu, WH, Shillah, SA, Marandu, WYF and Christiansen, JL (2006) Morphological diversity of Bambara groundnut (Vigna subterranea (L) Verdc.) landraces in Tanzania. Genetics Resources and Crop Evolution 53, 367378.CrossRefGoogle Scholar
Odongo, FO, Oyoo, ME, Wasike, V, Owuche, JO, Karanja, L and Korir, P (2015) Genetic diversity of Bambara groundnut (Vigna subterranea (L.) Verdc.) landraces in Kenya using microsatellite markers. African Journal of Biotechnology 14, 283291.Google Scholar
Ofori, I (1996) Correlation and path-coefficient analysis of components of seed yield in Bambara groundnut (Vigna subterranea). Euphytica 91, 103107.CrossRefGoogle Scholar
Olanrewaju, OS, Oyatomi, O, Babalola, OO and Abberton, M (2021) Genetic diversity and environmental influence on growth and Yield parameters of Bambara groundnut. Frontiers in Plant Science 12, 796352.CrossRefGoogle ScholarPubMed
Olukolu, BA, Mayes, S, Stadler, F, Ng, NQ, Fawole, I, Dominique, D, Azam-Ali, SN, Abbott, AG and Kole, C (2012) Genetic diversity in Bambara groundnut (Vigna subterranea (L.) Verdc.) As revealed by phenotypic descriptors and DArT marker analysis. Genetics Resources and Crop Evolution 59, 347358.CrossRefGoogle Scholar
Onwubiko, NC (2021) Genetic diversity in Bambara groundnut (Vigna subterranea (L) Verdc.). Agricultura tropica et subtropica 54, 8996.CrossRefGoogle Scholar
Ouedraogo, M, Ouedraogo, JT, Tignere, JB, Bilma, D, Dabire, CB and Konate, G (2008) Characterization and evaluation of accessions of Bambara groundnut (Vigna subterranea (L.) Verdcourt) from Burkina Faso. Sciences and Nature 5, 191197.CrossRefGoogle Scholar
Peakall, R and Smouse, P (2012) GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research-an update. Bioinformatics (Oxford, England) 28, 25372539.Google ScholarPubMed
Pungulani, L, Kadyampakeni, D, Nsapato, L and Kachapila, M (2012) Selection of high yielding and farmers’ preferred genotypes of Bambara nut (Vigna subterranea (L.) Verdc) in Malawi. American Journal of Plant Science 3, 18021808.CrossRefGoogle Scholar
R Core Team (2021) R: A language and environment for statistical computing. Published online 2020.Google Scholar
Rungnoi, O, Suwanprasert, J, Somta, P and Srinives, P (2012) Molecular genetic diversity of Bambara groundnut (Vigna subterranea L. Verdc.) revealed by RAPD and ISSR marker analysis. SABRAO Journal of Breeding and Genetics 44, 87101.Google Scholar
Shegro, AG, Jansen Van Rensburg, WS and Adebola, PO (2013) Assessment of genetic variability in Bambara groundnut (Vigna subterranea (L.) Verdc.) using morphological quantitate traits. Academia Journal of Agricultural Research 1, 4551.Google Scholar
Somta, P, Chankaew, S, Rungnoi, O and Srinives, P (2011) Genetic diversity of the Bambara groundnut (Vigna subterranea (L.) Verdc.) as assessed by SSR markers. Genome 54, 898910.CrossRefGoogle ScholarPubMed
Stadler, F (2009) Analysis of differential gene expression under water-deficit stress and genetic diversity in Bambara groundnut [Vigna subterranea (L.) Verdc.] using novel high-throughput technologies. Ph.D. thesis, Technischen Universität München, Germany.Google Scholar
Tan, XL, Azam-Ali, S, Goh, EV, Mustafa, M, Chai, HH, Ho, WK, Mayes, S, Mabhaudhi, T, Azam-Ali, S and Massawe, F (2020) Bambara groundnut: an underutilized leguminous crop for global food security and nutrition. Frontiers in Nutrition 7, 116.CrossRefGoogle ScholarPubMed
Uba, CU, Oselebe, HO, Tesfaye, AA and Abtew, WG (2021) Genetic diversity and population structure analysis of Bambara groundnut (Vigna subterrenea L) landraces using DArT SNP markers. PLoS One 16, e0253600.CrossRefGoogle ScholarPubMed
Unigwe, AE, Gerrano, AS, Adebola, P, Pillay, M and Monrovia, L (2016) Morphological variation in selected accessions of Bambara groundnut (Vigna subterranea L. Verdc) in South Africa. Journal of Agricultural Science 8, 6980.CrossRefGoogle Scholar
Valombola, J, Akundabweni, L, Awala, S and Hove, K (2019) Agronomic and morphological diversity of Bambara groundnut (Vigna subterranea (L.) Verdc.) accessions in North-Central Namibia. Welwitschia International Journal of Agricultural Sciences 1, 8899.Google Scholar
Vurayai, R, Emongor, V and Moseki, B (2011) Physiological responses of Bambara groundnut (Vigna subterranea (L) Verdc.) to short periods of water stress during different developmental stages. Asian Journal of Agricultural Science 3, 3743.Google Scholar
Figure 0

Table 1. Mean square values for traits scored in 2019/2020 season

Figure 1

Table 2. Mean square values for traits scored in 2021/2022 season

Figure 2

Figure 1. A scatter plot matrix showing the correlation coefficients among traits and the significance levels scored in 2021/2022 season.

Figure 3

Figure 2. Hierarchical clustering dendrogram analysis; Euclidean distance was used and the associations between groups were done by the Ward method for the morphological traits.

Figure 4

Table 3. Pairwise Nei genetic distances among three subpopulations

Figure 5

Table 4. Analysis of molecular variance (AMOVA) for the three geographical regions of the Bambara groundnut Genotypes

Supplementary material: File

Chipeta and Gimode supplementary material
Download undefined(File)
File 259.6 KB
Supplementary material: File

Chipeta and Gimode supplementary material
Download undefined(File)
File 193.4 KB
Supplementary material: File

Chipeta and Gimode supplementary material
Download undefined(File)
File 278.3 KB
Supplementary material: File

Chipeta and Gimode supplementary material
Download undefined(File)
File 40.8 KB