Insights into the genetic diversity of an underutilized Indian legume, Vigna stipulacea (Lam.) Kuntz., using morphological traits and microsatellite markers

Vigna stipulacea (Lam.) Kuntz., commonly known as Minni payaru is an underutilized legume species and has a great potential to be utilized as food crop. To evaluate and select the best germplasm to be harnessed in the breeding programme, we assessed the genetic diversity of V. stipulacea (94 accessions) conserved in the Indian National Genebank, based on morphological traits and microsatellite markers. Significant variation was recorded for the morphological traits studied. Euclidean distance using UPGMA method grouped all accessions into two major clusters. Accessions were identified for key agronomic traits such as, early flowering (IC331436, IC251436, IC331437); long peduncle length (IC553518, IC550531, IC553557, IC553540, IC550532, IC553564); and more number of seeds per pod (IC553529, IC622865, IC622867, IC553528). To analyse the genetic diversity among the germplasm 33 SSR primers were used anda total of 116 alleles were detected. The number of alleles varied from two to seven, with an average of 3.52 per loci. The polymorphic information content values varied from 0.20 to 0.74, with a mean of 0.40. The high number of alleles per locus and the allelic diversity in the studied germplasm indicated a relatively wider genetic base of V. stipulacea. Phylogenetic analysis clustered accessions into seven clades. Population structure analysis grouped them into five genetic groups, which were partly supported by PCoA and phylogenetic tree. Besides, PCoA and AMOVA also decoded high genetic diversity among the V. stipulacea accessions. Thus, morphological and microsatellite markers distinguished V. stipulacea accessions and assessed their genetic diversity efficiently. The identified promising accessions can be utilized in Vigna improvement programme through introgression breeding and/or can be used for domestication and enhanced utilization of V. stipulacea.


Introduction
The Vigna is an important genus among all pulse crops, containing more than 200 domesticated and wild species [1]. Among these, the domesticated species, namely, green gram or mungbean [ umbellata (Thunb.) Ohwi and Ohashi] are of commercial importance and contribute significantly to the pulse production [1,2]. While these crops have a positive impact on food and nutritional security worldwide and also on environmental sustainability, their productivity and production remain low owing to several biotic and abiotic constraints [3,4]. Several Vigna species have the status of under-exploited, undomesticated or semi-domesticated species or simply 'wild Vignas' [5][6][7][8]. Collection, evaluation and investigation of hitherto wild and exotic species from different agro-climatic regions to domesticate the undomesticated ones remains an important method, which could help mitigate the food and nutritional insecurity concerns [9].
Among various wild and useful Vigna species, Vigna stipulacea (Lam) Kuntz., is an important one with tremendous potential to be utilized in the Vigna improvement programme. The genome size of V. stipulacea is~445.1 Mbp [10], and it belongs to the secondary gene pool of the cultivated green gram and black gram. This species is a reservoir of various useful traits viz., faster growth, shorter duration and broad resistance to various pests and diseases and can be targeted for its domestication as a new crop and while the desirable traits can be transferred to the cultivated Vigna species [10]. Though V. stipulacea has been mistaken with V. trilobata since a long, it is taxonomically distinct [11]. It is being grown traditionally in lowland paddy fields before paddy cultivation in South India and locally known as 'Minni payaru' or "Siru Payaru" [12,13]. This crop is primarily used for fodder and manure production [14,15] and also being utilized as an ingredient for the preparation of vegetable stew (sambhar), condiment and dosa [12]. This species has been recognized as a source of disease and pest resistance [16][17][18][19]. Despite its economic importance, the information on the extent of variability on the V. stipulacea germplasm remains veiled. The National Gene Bank (NGB) of India housed at ICAR-National Bureau of Plant Genetic Resources (ICAR-NBPGR), New Delhi has conserved germplasm of V. stipulacea. This gene bank collection needed to be evaluated using the agro-morphological traits and efficient molecular markers like microsatellite markers (Simple Sequence Repeats) to assess the extent of diversity and identify promising germplasm lines. Information on the extent of genetic diversity of a panel of genotype and its population genetic structure provides useful information for deploying marker assisted breeding for genetic improvement [1,20]. To use wild and underutilized species in crop improvement programmes through introgression breeding, precise information on their genetic architecture and population genetic structure and relationship with other Vigna species are a deciding factor [21]. The objectives of the present investigation were to evaluate the genetic diversity and study the population genetic structure of V. stipulcaea germplasm using agro-morphological and microsatellite markers which will provide a better insight into the extent of relatedness and diversity in the studied accessions. The identified germplasm as a result of the above efforts can be deployed for transferring traits of interest in cultivated background from the untapped wild genetic resources of V. stipulacea. This could help to shape the Vigna improvement progamme better and contribute significantly to ensuring global food and nutrition security.

Plant materials
A diverse panel of 94 accessions (acc.) of V. stipulacea conserved in the NGB housed at ICAR-NBPGR, New Delhi was used in the present study. These accessions belonged to different diversity-rich agro-climatic regions of India and were collected from different states, including Andhra Pradesh (61 acc.), Tamil Nadu (14 acc.), Chhattisgarh (5 acc.), Madhya Pradesh (5 acc.), Odisha (3acc.) and Gujarat, Karnataka and Kerala (2 acc. each). The detailed passport data of all the accessions are presented in the Table 1.
All the accessions were grown at ICAR-Indian Institute of Pulses Research (ICAR-IIPR), Kanpur, Uttar Pradesh (26˚27'N latitude, 80˚14'E longitude, 452.4 meter above mean sea level (AMSL) for two consecutive years, i.e. during the rainy season (Kharif) 2018 and 2019. The long term mean annual rainfall of Kanpur is 820 mm, and >80% of it is received during the southwest monsoon season (July-September). All accessions were raised in customized cement pots using Augmented Block design (ABD)in the wide hybridization garden facility at ICAR-IIPR, Kanpur (S1 Fig) Each pot measured 1 m in diameter and was about 1 m high, and was filled with a potting mixture of sand: soil: FYM in 1:1:1 ratio. Further, each pot has a provision of drainage at the bottom and the sides so that water stagnation does not occur during excess precipitation. Before sowing, mechanical scarification of seeds was done in all accessions to ensure maximum germination [1,22].

Recording of data
The observations for 33 agro-morphological traits were recorded using descriptors of Bioversity International and ICAR-NBPGR [23].
The following 15 quantitative traits were recorded according to descriptors: 12. Number of seeds per pod: recorded after harvesting by manually counting the seeds per pod and averaged for 10 pods per accession. 13. Pod length (cm): recorded for ten random pods and averaged in each accession in cm.
14. Seed length (mm): measured using Vernier calliper and expressed in mm. 15. Seed width (mm): measured using Vernier calliper and expressed in mm.

Qualitative traits
Eighteen qualitative traits included, seed germination habit, hypocotyl colour, attachment of primary leaves, early plant vigour, stipule shape, leaf colour, branch pigmentation, petiole colour, leafiness, leaf pubescence, plant growth habit, attachment of pod to peduncle, pod pubescence, pod colour, constriction between seeds in the pod, seed shape, luster on seed surface and mottling on seed surface. Data on qualitative traits were recorded as per the descriptors and the details are presented in the S1 Table.

Extraction of DNA
For DNA extraction Cetyl trimethylammonium bromide (CTAB) method was followed [24]. DNA quantification were completed with 1% agarose gel electrophoresis and confirmed by using spectrophotometer (Nanodrop™, Thermo-Fisher, USA) and final DNA concentration was adjusted to 20 ng/μL.

SSR markers and Polymerase Chain Reaction (PCR) amplification
Initially, 100 microsatellite markers from different Vigna backgrounds viz., cowpea [25], adzuki bean [26], mungbean [27][28][29][30][31] and common bean [32,33] (listed in supplementary data, S2 Table) [31] with minor modifications. The details of standardized program for amplification were as follows: initial denaturation was carried out for 6 min at 94˚C, followed by 35 cycles of denaturation for 1 min at 94˚C, annealing for 1 min at 48-52˚C temperature (standardization of annealing temperature was carried out by gradient PCR for each SSR marker), extension for 1 min at 72˚C and final extension for 10 min at 72˚C. The amplified product was separated on 3% metaphor agarose (Lonza, USA) gel for 3h at 110 V and gel images were captured by using gel documentation system (Alpha Imager1, USA).

Data scoring and statistical analyses
The data on morphological traits of both the years were pooled to calculate mean, standard error, standard deviation and coefficient of variation (%). These accessions were classified using hierarchical cluster analysis, Euclidean distances were estimated and a dendrogram was constructed following UPGMA tree method. PCA was used to assess the phenotypic diversity. Each principal component was computed by using a linear combination of an Eigen vector of the correlation matrix with a variable. The three-dimensional plot was also developed to depict the view of accession scores. Both hierarchical clustering and PCA were analyzed using NTSYS pc-2.11x software [34]. Among the 94 V. stipulacea accessions, amplified products of 33 SSR loci were scored manually comparing with DNA ladder 100 bp (G-Biosciences, USA). The genetic diversity indices, viz., minor allele frequency (Maf), number of alleles (An), observed heterozygosity (Ho), genetic diversity (GD) or expected heterozygosity and polymorphism information content (PIC) of each SSR locus were calculated using PowerMarker v3.25 [35]. In addition, genetic distances across the V. stipulacea accessions and Neighbour-joining (NJ) tree were constructed using Power Marker v3.25 [35]. Population structure for V. stipulacea accessions was determined using the Bayesian model-based programme, software STRUCTURE 2.3.4 [36]. The number of clusters (K) was ranged from two to ten and three replications were run for analysis of each K value. The STRUCTURE was run with 100,000 generations of each burn-in and a Markov chain Monte Carlo. LnP(D) was obtained for each K and ΔK values are derived after plotting LnP(D) [37]. Final number of the populations of the studied accessions was obtained by using the STRUCTURE harvester. The data sets were subjected to Principle Coordinate Analysis (PCoA) and analysis for molecular variation (AMOVA) using software GenAlEx V6.5 [38].

Quantitative traits exploration
Descriptive statistics of 15 quantitative traits indicated that significant variation was present among all accessions (  Fig)). Among different morphological traits, the stipule length varied from 1.10 to 2.63 mm with a mean of 1.70 mm. Accession IC553553 collected from Andhra Pradesh recorded the highest stipule length. Accessions IC251435, IC331436, IC331437 (belonging to Odisha) and IC521245, IC37804 (Tamil Nadu) had lowest stipule length. In our study, a high coefficient of variation (30.22%) was observed for peduncle length. The longest peduncle length (63.0 cm) was recorded in the accession IC553564. Most of the accessions from Andhra Pradesh had the peduncle length more than the mean value (33.55 cm). Length of the longest branch showed significant variation in all the accessions

Hierarchical clustering and PCA
Dendrogram based on the Euclidean distances grouped 94 accessions of V. stipulacea into two major clusters. Both the clusters were subsequently subdivided into two sub-clusters each, represented as Ia, Ib, IIa and IIb (Fig 2). The accessions belonging to Chhattisgarh were grouped into sub-cluster Ia whereas, accessions from Odisha were grouped in sub-cluster Ib, which comprised of early maturing accessions. Accessions in sub-clusters IIa and IIb did not group according to their geographical origin. The first three principal components summarized a realistic summary of all data explaining 56.26% variance in 94 accessions of V. stipulacea. Most important was the first principal

PLOS ONE
component (PC1) accounting for 27.18% of the total variation. PC2 and PC3 accounted for 17.68% and 11.40% variance, respectively. PC 1 explained variation in stipule length and terminal leaf length, PC2 in days to initial flowering, days to first podding and branch length while PC3 in pod length and number of seeds per pod (Table 3).

Qualitative traits exploration
Qualitative traits were recorded at the appropriate stages of the plant growth when these had full expression. Accordingly, these traits were recorded at the germination, vegetative, podding and maturity stage as follows: Germination stage. The seed germination habit, hypocotyl colour, attachment of leaf at the two-leaf stage (primary leaf) were observed at the time of germination. Most of the accessions showed an intermediate type of germination and greenish-purple coloured hypocotyls. Petiolate type of primary leaf attachment was observed in all accessions (Fig 3).  Vegetative stage. The traits studied at the vegetative stage were early plant vigour, stipule shape, leaf colour, branch pigmentation, petiole colour, leafiness, leaf pubescence and plant growth habit. Fifty per cent of the accessions had good plant vigour 30 days after sowing (DAS). Stipule shape of all accessions was lanceolate. Green and dark green coloured leaves were observed. Further, the majority of the accessions (77 accessions) had green coloured leaves. Greenish purple pigmentation on the stem was recorded while greenish purple coloured petiole was observed on 84 accessions. Significant variation was observed for the leafiness traits among different accessions. A higher proportion (43 accessions) was of sparse leafy types, while the intermediate (25) and abundant (26) leafy type accessions were almost in equal numbers. Accessions with high biomass and foliage had green fodder potential. There was predominance of accessions with glabrous type of leaves, stem and petiole. As far as plant type is concerned, most of the accessions under study were spreading (41) and semi-erect (42) while, 11 accessions were erect type (Fig 4).

PLOS ONE
Pod formation and maturity stage. At pod formation and maturity stage, attachment of pod to peduncle, pod pubescence, constriction of pod between seeds, pod colour, seed colour, seed shape, lustre on seed surface and mottling on seed surface were recorded for all accessions. Pods of all accessions were pendently attached with the peduncle. There was a predominance of accessions with sparse hairs on pods whereas, for pod colour, black was more prevalent (72 accessions). All accessions had mottled, lustrous seed surface with oval seed shape (Fig 5).
Germplasm characterization on the basis of morphological observations. Based on two years' evaluation, promising accessions were selected for the traits analyzed ( Table 4). The accessions which started flowering early (less than 30 days) were IC251436, IC331436 IC331437. More than 50 cm. peduncle length was possessed by the accessions IC553518, IC550531, IC553557,

SSR characterization
Genetic diversity analysis using 33 SSRs revealed significant variation present in the V. stipulacea germplasm ( Table 5

Phylogenetic studies
The Neighbor-joining (N-J) tree genetically differentiated 94 accessions of V. stipulacea into two major clusters (Fig 6). The N-J tree was further simplified into seven clades for its better explanation and understanding. Clade 1, 2, 3 and 5 contained accessions belonging to Andhra Pradesh while Clade 4 possessed accessions of Andhra Pradesh along with the accessions of Madhya Pradesh. Ten accessions of Tamil Nadu and 4 of Chhattisgarh were grouped in the Clade 6. Clade 7 had the mixture of accessions from different geographical regions viz., Tamil Nadu, Gujarat, Kerala, Karnataka, Odisha and Chhattisgarh.

Population structure analysis
The 94 accessions of V. stipulacea were grouped into five population genetic groups recording the most probable number of K (S3 Fig) and grouping of each individual was illustrated using bar plot diagram. Of the 94 accessions, 80 (85.10%) shared their ancestry with one of the five Table 4

PLOS ONE
The allele-frequency divergence of population I with population II, III, IV and V was 0.2093, 0.2184, 0.2778 and 0.2182, respectively. The allele-frequency divergence of population II with population III, IV and V was 0.3366, 0.3447 and 0.3399, respectively. Subsequently population III had 0.1651 and 0.1056 allele-frequency divergence with population IV and V, respectively. Population IV and V had 0.1976 allele-frequency divergence (Table 6).

AMOVA and PCoA
AMOVA of five V. stiupulacea populations (as per the model-based population structure analysis) described 39% variations among the populations; 55% variations among individuals and 06% within individuals (S5 Table).
Besides, principal coordinate analysis (PCoA) showed high genetic variation among 94 accessions of V. stipulacea. Cumulative variation of three axes of PCoA was 40.01%, wherein first axes explained 20.04%; second 10.58% and third 9.39% of genetic variation. In the PCoA,  Table 1.

Discussion
Wild and underutilized germplasm offers novel sources of diversity hitherto not found in the domesticated species and consequently offers further opportunities of selection for agronomic traits [39]. Assessment of genetic diversity generates knowledge which is valuable for managing and conserving genetic resources [40]. In addition, it helps in identifying contrasting parental materials which can be utilized in breeding programmes. This genetic diversity can also be deployed to develop inbred lines with required genes and/or gene complexes introgressed from diverse and underutilized germplasm. The information of genetic diversity of germplasm has been extensively recognized in the selection of genotypes for yield and yield attributing characters, resistance to biotic and abiotic stresses as well as for nutritional improvement of crop species [41]. In the same line of thought, we carried out genetic diversity analysis of potential legume crop V. stipulacea based on morphological traits and microsatellite markers intending to identify and select diverse genotypes for use in further breeding programs. Phenotypic evaluation of the genetic resources based upon morphological descriptors remains the first step for the evaluation, description and categorization of germplasm collections to augment their utilization in crop breeding program [23,42]. The morphological  11]. All such studies suggested that wide range of exploitable variation is available in wild Vigna species. V. stipulacea was proposed as a candidate for neodomestication for disease and pest resistance by Tomooka et al. (2014) [5], whereas resistance to diseases and pests this species was also reported in several other reports (Reviewed by Pratap et al., 2021) [43]. Nutritional potential of V. stipulacea was unveilded by Gore et al., (2021b) [44]. Pandiyan et al., (2010) [45] reported that V. stipulacea is very closely related and crossable to V.radiata.
In the present study, stipule length varied from 1.10 to 2.63 mm with a mean of 1.70 mm. Among different morphological traits, the stipule length has been reported as one of the key characters for identification of V. stipulacea [46,47]. Long peduncles are advantageous in getting the pods above the canopy. This helps in reducing the damage of pods by the pod borer and other insects and is also helpful at the time of harvesting [48]. This trait may be transferred to cultivated Vigna for facilitating mechanical harvesting, especially in mungbean,which is highly desirable keeping in view the increasing cost of labour and drudgery involved in manual picking of mature pods. In our study, a high coefficient of variation (30.22%) was observed for peduncle length. The longest peduncle length (63.0 cm) was recorded in the accession IC553564. Accessions with longer peduncles can be utilized for Vigna improvement programme and also for fibre production. Long peduncle of cowpea cultivar (V. unguiculata cv. textilis) was reported to be used for fibre in Africa [49]. Among the quantitative traits, days to initial flowering and maturity were the most important traits of observation. Crop duration has a direct impact of fitting a cultivar in a cropping system, and early duration varieties of mungbean have been developed which can fit well in short-season windows, thereby offering the opportunities of horizontal increase in area [39]. Days to initial flowering significantly varied among accessions with the mean value of 37 days. Likewise, days to initial maturity also recorded significant variation and ranged between 49 to 76 days. Among all accessions, IC331436 matured at the earliest as it matured within 49 days after sowing. Other accessions namely, IC251436, IC331437, IC622861, IC276983, IC622865, IC622867, IC622860, IC331454, IC331457 matured in <60 days. Therefore, these accessions may be utilized as donors for developing early maturing genotypes in green gram and black gram. Shorter duration cultivars can help in avoiding the terminal heat stress, especially in spring crop and also the terminal rains at the time of harvest during the rainy season in India [50]. Therefore, accessions with early maturity can be tremendously useful in expanding area in the northern and central part of India during the summer season. The wide variability for pod length was also observed, pod of V. stipulacea is harder and more resistant to the attack by stink bug as compared to the green gram and black gram [14]. Accessions IC553529, IC622865 and IC622867 recorded 15 number of seeds per pod, followed by IC553528 (14 seeds) and therefore, these were observed as the most promising for hybridization programmes towards increasing number of seeds per pods and ultimately seed yield as well as for selecting for direct use as a crop in domestication. Several earlier studies have reported significant positive correlation between number of seeds/ pod and grain yield in mungbean [51], blackgram [52] and cowpea [53]. In V. sublobata, V. radiata and V. trilobata, up to 13 number of seeds per pod were recorded [23]. In the cultivated Vigna species like mungbean and urdbean usually 9-11 seeds per pod are observed and increasing 2-3 seeds per pod may increase the seed yield by upto 10 per cent. Therefore, the accessions with more number of seed per pod can be utilized in improving the cultivated Vigna species. These results of correlation analysis are in consonance with those of Bisht et al., 2005 [23] on studies of different wild Vigna species. An intermediate type of germination was most common in V. stipulacea which is the stage between hypogeal and epigeal germination [54]. In most of the accessions cotyledons were above the ground surface but hypocotyl was underground. These observations were similar to reports of Dixit, 2014 [46], Dixit et al., 2011 [55]. Primary leaf attachment is petiolate type, Tomooka et al., 2002 [56] suggested that this trait can be useful for taxonomic delineation among closely related species.
Spreading and semi-erect type of plant growth habit was most prevalent among accessions. Earlier studies also reported spreading type plant growth habit of V. stipulcaea [47,56]. Erect type of plant is important for mechanical harvesting as well for use as fodder crops. In our study, 11 accessions were found to be of erect plant habit that can be utilized in the Vigna improvement programme. There was a predominance of accessions with sparse hairs on pods. Pubescence plays a vital role in protecting the crop from different biotic and abiotic stresses along with its adaptation mechanism. Pod pubescence as unique character in Vigna unguiculata (L.) Walp was reported by Tripathi et al., 2021 [57]. Among the pods black colour was more prevalent at maturity (72 accessions). Pod colour acts as a morphological marker and may be deployed in quality seed production programmes to monitor the mixture of other varieties. All accessions had mottled, lustrous seed surface.
Assessment of genetic resources of V. stipulacea for identification of useful germplasm is essential for enhancing their utilization by the plant breeders towards genetic improvement of different Vigna crops and also for genebank managers to refine genetic resource management strategy. This study led to the identification of desirable germplasm lines for specific traits like long peduncle, green fodder potential, early maturity and more number of seeds per pod. More interestingly, accessions with multiple traits were also identified which included, for example, green fodder potential and more number of seeds per pod (IC553528 and IC553529); long peduncle length and green fodder potential (IC550531); early flowering and early maturity (IC331436, IC251436, IC331437) and early maturity and more number of seeds/pods (IC622865, IC622867). All the identified trait-specific germplasm was validated again in Kharif 2020 at three locations viz., IIPR, Kanpur, NBPGR, main campus, New Delhi, and NBPGR, Issapur Farm. Therefore, the promising accessions identified through this study can be used in Vigna breeding programs to select high yielding varieties with enhanced resistance to biotic stresses and having green fodder potential. Simultaneously, the potential of outstanding accessions can also be exploited by their direct introduction in the farmers' fields for their multifarious uses. Nonetheless, keeping in view that the study involves an evaluation of an Indian germplasm, there is a further need of comprehensive evaluation of a large number of global germplasm lines to identify trait-specific potential germplasm adapted in India and elsewhere.
Among various molecular markers, microsatellite markers (SSRs) were used by many earlier workers and recognized as highly informative and efficient for studying genetic diversity due to their co-dominant nature, high reproducibility, high rate of cross-transferability, easiness in scoring [58]. In the past, Pratap [59] by studying 15 SSR markers on 65 mungbean genotypes for genetic diversity analysis. Seventy genotypes of the Vigna genus for RAPD, URP and SSRs were studied and reported that SSRs detected 211 polymorphic bands in the Vigna species [29]. More number of alleles per locus is one of the important criteria for selection of the marker to be utilized in the assessment of genetic variability [60]. Primers with higher PIC value and multiple-locus amplification have been recommended for germplasm characterization, evaluation, phylogenetic studies and breeding applications. In the current investigation, the PIC value of primers indicated greater diversity within and between studied germplasm and also revealed that selected primers are very informative and can be utilized in the further studies on Vigna genotypes.
The N-J tree clustered 94 accessions of V. stipulacea into seven clades which suggested a high degree of their genetic divergence. The SSRs amply discriminated the accessions under study. Strikingly, the overall clustering of the V. stipulacea accessions was not as per their geographical origins, as many of the accessions from the same geographical origin (for example accession from Tamil Nadu) did not group together. Earlier, Pratap et al. (2015) [1] also reported that accessions of V. trilobata (a close relative of V. stipulacea) grouped in different clusters despite of same geographic origin and ascribed it to their collection across the agroecological regions and also that some of them might have experienced natural hybridization in due course of time which explains their admixture status. Studies on 692 accessions of mungbean by Gwag et al., (2010) [61] collected from 27 countries reported that due to germplasm exchange between countries, genetic associations in accession under study were not associated with geographical origin. Genetic diversity not always follows the geographical origin of germplasm which may be due to extensive germplasm exchange from different geographical regions [62]. The 3-D PCA grouping patterns were also in general agreement with the grouping described in the dendrogram. Similarly, Pratap et al., (2015); Dikshit et al., (2007) [1,60] while studying the Vigna species also reported that the results of the cluster analysis were comparable to the PCA. Population structure analysis based on the model is one of the excellent tools for the assessment of the genetic diversity in any plant species [63], which could be magnificently applied in conservation and utilization of the collected germplasm [64]. In the recent past, population structure analysis was extensively utilized for genetic diversity analysis of the Vigna germplasm. For example, 53 accessions of Vigna species were genetically grouped into five groups using the model-based population structure analysis [1]. Similarly, AMOVA also described high genetic variability among the populations, among and within the individual of the populations. Besides, PCoA partly revalidated this genetic grouping of the V. stipulacea accessions. The potential for utilization of different wild Vigna species has been described by Takahashi et al., (2016) [65] based on estimation of phylogenetic positions by using rDNA-ITS (nuclear) and atpB-rbcL (chloroplast) spacer regions.
Our result revealed significant variations in the morphological traits and microsatellite DNA polymorphisms in V. stipulacea accessions. The high allelic diversity in the analyzed germplasm lines indicated that the genetic spectrum was relatively wide. V. stipulacea has a great potential which can be exploited for improving other domesticated Vigna species by broadening their genetic base which can be useful to fight with biotic and abiotic stresses, especially in the climate change scenario. The study identified superior accessions for several traits which can be deployed in hybridization programme for improvement of mungbean and black gram.
Supporting information S1