Abstract

Genus Vigna represented by more than 100 species is a source of nutritious edible seeds and sprouts that are rich sources of protein and dietary supplements. It is further valuable because of therapeutic attributes due to its antioxidant and anti-diabetic properties. A highly diverse and an extremely ecological niche of different species can be valuable genomic resources for productivity enhancement. It is one of the most underutilized crops for food security and animal feeds. In spite of huge species diversity, only three species of Vigna have been sequenced; thus, there is a need for molecular markers for the remaining species. Computational approach of microsatellite marker discovery along with evaluation of polymorphism utilizing available genomic data of different genotypes can be a quick and an economical approach for genomic resource development. Cross-species transferability by e-PCR over available genomes can further prioritize the potential SSR markers, which could be used for genetic diversity and population differentiation of the remaining species saving cost and time. We present VigSatDB—the world’s first comprehensive microsatellite database of genus Vigna, containing >875 K putative microsatellite markers with 772 354 simple and 103 865 compound markers mined from six genome assemblies of three Vigna species, namely, Vigna radiata (Mung bean), Vigna angularis (Adzuki bean) and Vigna unguiculata (Cowpea). It also contains 1976 validated published markers. Markers can be selected on the basis of chromosomes/location specificity, and primers can be generated using Primer3core tool integrated at backend. Efficacy of VigSatDB for microsatellite loci genotyping has been evaluated by 15 markers over a panel of 10 diverse genotype of V. radiata. Our web genomic resources can be used in diversity analysis, population and varietal differentiation, discovery of quantitative trait loci/genes, marker-assisted varietal improvement in endeavor of Vigna crop productivity and management.

Introduction

Genus Vigna is one of the most important Fabaceae family crops, which is mostly found in tropical regions and warm temperate region of Asia and Africa (1), but some of its species are also found in America (2). Myanmar, India and Brazil are the top three dry bean-producing countries in the world with 5.2, 3.9 and 2.6 MT, respectively (3).

Genus Vigna has six subgenera, namely, Ceratotropis, Vigna, Lasiosporon, Plectotropis, Haydonia and Sigmoidotropis (4). Genus Vigna having more than 100 species is one of the most underutilized leguminous crop for human food and animal feeds, as only 10 species are domesticated (5). Mungbean (Vigna radiata), cowpea (Vigna unguiculata), adzuki bean (Vigna angularis) and urdbean (Vigna mungo) are the major domesticated species. The vast genetic diversity in wild species of Vigna is its best attractive feature, and these wild species are highly tolerant to biotic and abiotic stress such as pests, diseases, drought, alkaline soil and high salinity. Since more than 90 species of this genus are well adaptive to very extreme ecological niches like marine beach, deserts and swamps, they could be a potential source of genomic resource for productivity enhancement (6).

Genus Vigna are known for its edible seeds and sprouts, which are highly nutritional, rich sources of proteins, relatively high in iron and high in vitamins, minerals and folate contents. This leguminous crop is known for its antioxidant properties and for having therapeutic attributes for many diseases like cancer, diabetes, Alzheimer’s diseases, arthritis and cardiovascular diseases (7). This crop also plays a role in improving soil fertility and texture and fixing atmospheric nitrogen (8). Intercrop rotation with Vigna is well known to improve yield production of the other crops (9).

Till now, genomes of three Vigna species viz. mungbean (V. radiata), adzuki bean (V. angularis var. angularis) and V. unguiculata (Cowpea) were sequenced. In genus Vigna, mostly are diploid (2n = 2x = 22), while in Vigna reflexo-pilosa, mostly are tetraploid (2n = 4x = 44). The genome size of all the Vigna species may vary from species to species ranging from 416 to 1394 Mb (9).

Simple sequence repeats (SSRs), also known as microsatellites, are the co-dominant markers and plays an important role in biological functions. It is present in the stretch of DNA sequences that contain 1–6 repeating units of nucleotides (10). These repeating units are present in both coding and non-coding regions of DNA, but it is more abundant in intron region (11). SSR markers play an extensive role in marker-assisted selection (MAS), association studies, species identification, fingerprinting, diversity evaluation, linkage mapping and gene mapping for crop improvement due to its co-dominant inheritance, multi-allelic nature, robust amplification and reproducibility nature. SSR became the most powerful and important technique for plant genetic studies (12, 13). SSR markers are highly polymorphic as it depends on the number of tandem repeats. There are high chances of dissimilarity in number of repeats between species or individuals, thus this leads to the chance of getting high polymorphism between two respective individuals (14).

Although reports are available with a limited number of SSR markers, no bulk discovery of such markers is reported by genome-wide SSR mining. More markers are required for MAS, genetic diversity, genetic linkage maps, association or comparative mapping and qualitative and quantitative traits of Genus Vigna (14–17).

Traditional methods of microsatellite discovery using genomic libraries are compromised with a number of markers besides time and cost. To overcome this problem, in silico approach can be used. In silico method has the advantage of predicting target-specific regions in genome, which can be more effective in developing molecular markers required for quantitative trait loci (QTL) and linkage mapping (18, 19). Present linkage map needs more molecular markers to increase the map density for fine mapping (20).

As per TRIPS (trade-related aspects of the intellectual property rights) agreement and intellectual property rights, plant breeders have their ownership of variety on the basis of distinctiveness, uniformity and stability (DUS) features, which can be used for assigning of new variety status as well as resolving legal varietal disputes (21). Microsatellite markers have been successfully used in supplementing DUS feature in variety identification for crops like wheat, barley, soybean, rice and potato (22). Such approach can also be used for varietal differentiation of different Vigna species. Since Vigna genus has more than 100 species and limited genomic resources are available from three species, extensive molecular mining of SSR markers and evaluation of their polymorphism with cross-species transferability can be a more pragmatic approach to cater the need of markers in hitherto untouched species.

Since flanking region of microsatellite markers are well conserved across closely related species, heterologous PCR primers of such locus can cater the immediate need of molecular markers most economically (23). Molecular markers developed from `focal species’ have been successfully used in `non-focal’ species for population differentiation and diversity analysis (24).

Molecular markers, especially SSRs, play an important role for mapping of genes, genetic analysis and QTL mapping. Many SSR markers have been developed for QTL mapping and linkage maps associated with important traits like resistance and yield to abiotic and biotic stress of various Vigna species like mungbean (14, 25, 8, 26) and cowpea (27).

Although a mega plant genome SSR database (PMDBase) contains 110 plant species that includes two Vigna species, namely, V. radiata (cv. VC1973A) and V. angularis (cv. Jingnong 6), it has several limitations like lack of automation for SSR mining, designing of primers, input query sequence limited to 30 Kb (11); thus, there is a need to develop user friendly SSR database. If all available genomic data of genus Vigna are used to develop such database, it can be a valuable genomic resource for several species of genus Vigna.

Present work aims at genome-wide mining of SSR markers in three species of genus Vigna, namely, V. radiata (mung bean), V. angularis (adzuki bean) and V. unguiculata (cowpea). It further aims at development of a web genomic resource having option for chromosome-wise SSR mining and primer designing for genotyping along with e-PCR-based polymorphism discovery. It also aims to provide annotated genic region SSR to be used as functional domain markers (FDMs).

Materials and methods

Data source

Genomic data of V. radiata (mung bean), V. angularis (adzuki bean) and V. unguiculata (cowpea) were retrieved from NCBI (Table 1). Apart from full genome assembly data, we also used other genomic data of Vigna sp. available in the public domain. Vigna radiata contains genomic data of two cultivars viz., VC1973A and RIL59 (bruchid-resistant recombinant inbred line) (9, 28), whereas we found three assemblies of V. angularis cultivars such as Jingnong 6, Kyungwonpat and IT213134 (29, 30, 31). Vigna unguiculata has only one assembly of IT97K-499-35 cultivar (32). Out of six assemblies, three were sequenced chromosome-wise, i.e. VC1973A, Jingnong 6 and Kyungwonpat, and the rest are all scaffold- and contig-wise. Extensive literature survey was made to obtain validated SSR sequences in different species of genus Vigna to include in the web resource.

Table 1

Genomic data of different varieties of Vigna species used in the study

SpeciesCultivarsAccession no.Assembly levelSize (Mb)GC%
V. radiataVC1973AGCA_000741045.2Chromosome463.63834.40
V. radiataRIL59GCA_001584445.1Scaffold454.90732.00
V. angularisJingnong 6GCA_001190045.1Chromosome467.30134.62
V. angularisKyungwonpatGCA_001723775.1Chromosome444.43931.77
V. angularisIT213134GCA_000465365.1Scaffold291.82433.10
V. unguiculataIT97K-499-35GCA_001687525.1Scaffold695.04613.6
SpeciesCultivarsAccession no.Assembly levelSize (Mb)GC%
V. radiataVC1973AGCA_000741045.2Chromosome463.63834.40
V. radiataRIL59GCA_001584445.1Scaffold454.90732.00
V. angularisJingnong 6GCA_001190045.1Chromosome467.30134.62
V. angularisKyungwonpatGCA_001723775.1Chromosome444.43931.77
V. angularisIT213134GCA_000465365.1Scaffold291.82433.10
V. unguiculataIT97K-499-35GCA_001687525.1Scaffold695.04613.6
Table 1

Genomic data of different varieties of Vigna species used in the study

SpeciesCultivarsAccession no.Assembly levelSize (Mb)GC%
V. radiataVC1973AGCA_000741045.2Chromosome463.63834.40
V. radiataRIL59GCA_001584445.1Scaffold454.90732.00
V. angularisJingnong 6GCA_001190045.1Chromosome467.30134.62
V. angularisKyungwonpatGCA_001723775.1Chromosome444.43931.77
V. angularisIT213134GCA_000465365.1Scaffold291.82433.10
V. unguiculataIT97K-499-35GCA_001687525.1Scaffold695.04613.6
SpeciesCultivarsAccession no.Assembly levelSize (Mb)GC%
V. radiataVC1973AGCA_000741045.2Chromosome463.63834.40
V. radiataRIL59GCA_001584445.1Scaffold454.90732.00
V. angularisJingnong 6GCA_001190045.1Chromosome467.30134.62
V. angularisKyungwonpatGCA_001723775.1Chromosome444.43931.77
V. angularisIT213134GCA_000465365.1Scaffold291.82433.10
V. unguiculataIT97K-499-35GCA_001687525.1Scaffold695.04613.6

In silico SSR mining

SSRs were mined using perl script of MISA (MIcroSAtellite identification tool) with default parameters such as 10 repeating units for mono, 6 repeating units for di and 5 repeating units for tri, tetra, penta and hexa (33). Markers were obtained with their respective descriptive information such as marker type, repeat numbers, markers size, guanine-cytosine (GC) content, start and end position. Primer3core executable was integrated in this web genomic resource for the generation of primer pairs from flanking sequence of each locus. Primer generation parameters were as follows: melting temperature, 55–65°C; GC content, 40–70%; primer size, 18–27 bp; length and product size, 150–280 bp (34).

Functional annotation of genic region SSR markers

Available two species of Vigna, namely V. radiata (cv. VC1973A) and V. angularis (cv. Jingnong 6) that are fully annotated for gene function were incorporated in genome browser in the database. This uses general feature format coordinates of each gene file over coordinates of mined SSRs and displays gene ID along with structural and functional details.

In silico discovery of polymorphism and cross-species transferability

In-house perl script was used for prediction of in silico polymorphic SSR within the two species, namely V. radiata (VC1973A and RIL59) and V. angularis (Jingnong 6, Kyungwonpat and IT213134). Cross-species transferability was evaluated over genomic DNA sequence of V. angularis and V. unguiculata by e-PCR of microsatellite loci obtained from V. radiata cv. VC1973A. For this evaluation, loci were selected based on the following criteria: simple di-nucleotide repeat as they are expected to have more polymorphism due to higher probability of slippage in DNA replication (35) and loci having unique sequence of at least 100 bp in its flanking regions to avoid error in primer designing.

Work flow of markers and polymorphism discovery

In order to mine the microsatellite repeat loci, provision has been made for pattern identification using perl script implemented in MISA tool (33). It mines SSR loci from specific genome assembly of selected species/variety by user. In order to get SSR primers for genotyping, Primer3core (34) executable is integrated at the backend. User can select target of mining as per their choice of chromosome-wise or contig-wise. In-house perl script was used to extract flanking length of 500 bp upstream and 500 downstream of the SSR loci in targeted sequence file.

The `Polymorphism’ tab is based on the allelic length polymorphism identification algorithms by counting the amplicon size (in basepairs). This search is based on two modules for primer selection, namely, `self-designed primers’ and `External primers’. The `self-designed primers’ search generates new primers. VigSatDb has the references genome of six assembles with user-friendly option for uploading of genotype file for polymorphic discovery. Both, forward and reverse primers along with marker are located in the genotype file. The markers with different product size in the reference as well as genotype file are considered to be polymorphic markers. The second module, `External primers’ is based on primers that are already published. An in-house perl script has been written to detect the polymorphic markers by searching existing primers of loci in the flanking region of reference genome file. Then the e-PCR product sizes are computed and displayed as polymorphic loci. Using R code, provision of simulated gel has been made. Moreover, all the results can be achieved as HTML files, downloaded directly or through mail.

Database development

Vigna microsatellite database (VigSatDB) a `three tier architecture’ relational database has been developed having client tier, middle tier and database tier. Data of predicted markers and experimentally validated published markers were stored in the backend in MYSQL database and launched on apache server. The web server contains seven tabs viz. Home, Microsatellite, Tools, Polymorphism, Statistics, Tutorial and Team. The database is further linked to web interface using open source PHP scripting language. Primer3core executable was incorporated into the database for primer generation. It gives the output having three sets of primers with their respective PCR product size, GC content, melting point and start position (Figure 1). Further, MISA scripts, Primer3core, NCBI local Blast (36) and The Generic Genome Browser (37) tools were integrated into the database for GUI-based web server. Moreover, in-house perl scripts for polymorphism search were also incorporated in the database for finding self design primers and external primers to evaluate published primers.

A schematic diagram of SSR marker discovery and primer generation.
Figure 1

A schematic diagram of SSR marker discovery and primer generation.

Revalidation of microsatellite primers generated by VigSatDB

In order to check the efficacy of using VigSatDb for primer generation, in vitro validations were done. Fifteen primers were generated for microsatellite genotyping by PCR using a panel of ten genotypes of mungbean (V. radiata), namely NMK1512, SKM 1508, GSM1715, VIRAT, GAM5, GM-4, K-857, GJM1712, MEHA and VMS-15-1. These varieties were collected from Main Pulse Research Station, Junagadh Agricultural University, Junagadh, Gujarat, India. Young leaves were used to isolate DNA by CTAB method using protocol as described by Doyle and Doyle (38). The DNA concentration and purity of all genotypes were estimated spectrophotometrically using Picodrop (Picodrop Ltd Cambridge, UK).

PCR-based genotyping of microsatellite loci

PCR reaction for all the microsatellite loci were carried out in a 20-μl PCR reaction mixture containing 4 μl 2.5 mM dNTPs, 1 μl of 20 ng DNA in 2.5 μl 10× Taq polymerase assay buffer, 1.5 μl MgCl2 and 1.5 U Taq polymerase enzyme (Sigma Aldrich, Germany) along with 1 μl (1 μM) each of SSR forward and reverse primers. SSR-PCR product was further analyzed on a capillary electrophoresis (QIAxcel Advanced System, QIAGEN, Germany) using QIAxcel DNA Screening Kit and the images were analyzed using QIAxcel ScreenGel Software.

Results

In silico SSR mining and analysis of its abundance

SSR loci were mined successfully using MISA. SSR loci details were also mined from the public domain literature. Web genomic resource was developed with provision to generate primers for genotyping of loci. A total of 876 219 putative microsatellite markers were mined having 772 354 simple and 103 865 compound repeats (Table 2). Genotype-wise motif types were obtained (Table 3).

Table 2

Distribution of SSR types in Vigna species

SpeciesTotalSimpleCompound
V. radiata cv. VC1973A200 178174 11326 065
V. radiata cv.RIL59200 642173 40527 237
V. angularis cv. Jingnong6147 634130 52117 113
V. angularis cv. Kyungwonpat143 109126 60116 508
V. angularis cv. IT21313492 46782 40610 061
V. unguiculata cv. IT97K-499-3592 18985 3086881
Total Markers876 219772 354103 865
SpeciesTotalSimpleCompound
V. radiata cv. VC1973A200 178174 11326 065
V. radiata cv.RIL59200 642173 40527 237
V. angularis cv. Jingnong6147 634130 52117 113
V. angularis cv. Kyungwonpat143 109126 60116 508
V. angularis cv. IT21313492 46782 40610 061
V. unguiculata cv. IT97K-499-3592 18985 3086881
Total Markers876 219772 354103 865
Table 2

Distribution of SSR types in Vigna species

SpeciesTotalSimpleCompound
V. radiata cv. VC1973A200 178174 11326 065
V. radiata cv.RIL59200 642173 40527 237
V. angularis cv. Jingnong6147 634130 52117 113
V. angularis cv. Kyungwonpat143 109126 60116 508
V. angularis cv. IT21313492 46782 40610 061
V. unguiculata cv. IT97K-499-3592 18985 3086881
Total Markers876 219772 354103 865
SpeciesTotalSimpleCompound
V. radiata cv. VC1973A200 178174 11326 065
V. radiata cv.RIL59200 642173 40527 237
V. angularis cv. Jingnong6147 634130 52117 113
V. angularis cv. Kyungwonpat143 109126 60116 508
V. angularis cv. IT21313492 46782 40610 061
V. unguiculata cv. IT97K-499-3592 18985 3086881
Total Markers876 219772 354103 865
Table 3

Distribution of motif type found in each cultivar

Unit sizeVC1973ARIL59Jingnong6KyungwonpatIT213134IT97K-499-35
Mono127 648127 84586 20478 01654 25659 935
Di51 74250 82542 70348 68525 64621 813
Tri17 87218 95616 53414 56011 3309260
Tetra1827182212121200867872
Penta9121003734414231187
Hexa177191247234137122
Unit sizeVC1973ARIL59Jingnong6KyungwonpatIT213134IT97K-499-35
Mono127 648127 84586 20478 01654 25659 935
Di51 74250 82542 70348 68525 64621 813
Tri17 87218 95616 53414 56011 3309260
Tetra1827182212121200867872
Penta9121003734414231187
Hexa177191247234137122
Table 3

Distribution of motif type found in each cultivar

Unit sizeVC1973ARIL59Jingnong6KyungwonpatIT213134IT97K-499-35
Mono127 648127 84586 20478 01654 25659 935
Di51 74250 82542 70348 68525 64621 813
Tri17 87218 95616 53414 56011 3309260
Tetra1827182212121200867872
Penta9121003734414231187
Hexa177191247234137122
Unit sizeVC1973ARIL59Jingnong6KyungwonpatIT213134IT97K-499-35
Mono127 648127 84586 20478 01654 25659 935
Di51 74250 82542 70348 68525 64621 813
Tri17 87218 95616 53414 56011 3309260
Tetra1827182212121200867872
Penta9121003734414231187
Hexa177191247234137122

The genotype having chromosome-wise whole-genome data of all three species, viz., V. radiata (VC1973A), V. angularis (Jingnong 6) and V. angularis (Kyungwonpat), was used for the initial mining of SSR loci. Species and genotype-wise variation in abundance of SSR loci are shown in Figure 2.

Chromosome-wise abundance of SSRs in V. radiata (cv. VC1973A), V. angularis (cv. Jingnong 6) and V. angularis (Kyungwonpat) varieties.
Figure 2

Chromosome-wise abundance of SSRs in V. radiata (cv. VC1973A), V. angularis (cv. Jingnong 6) and V. angularis (Kyungwonpat) varieties.

In cv. VC1973A, maximum repeats were found on chromosome 7 (26 280), followed by chromosome 8 (20 521) and chromosome 5 (18 212). Chromosome 3 had minimum number of repeats (6579) having maximum SSR density per Megabase (Mb) (573.2). The least SSR density was found on chromosome 6 (411.3).

Although Jingnong 6 and Kyungwonpat cultivars were from the same species, at chromosome level substantial variation in SSR abundance was observed. In Jingnong 6, maximum markers were found on chromosome 1 (14 964), while minimum markers were observed on chromosome 6 (7391). In the case of cv. Kyungwonpat, chromosomes 6 and 10 were found to be with maximum (13 887) and minimum (3102) markers, respectively (Table 4). In the case of Jingnong 6, maximum and minimum SSR densities were found on chromosome 6 (381.8) and chromosome 8 (291.5), respectively. In the case of Kyungwonpat cultivar, maximum and minimum SSR densities were found on chromosome 5 (431.8) and 2 (317.9), respectively. In this study, mono-nucleotides and di-nucleotide repeats were found to be highly abundant, while penta-nucleotide and hexa-nucleotide repeats shared less proportion of SSR. AT and AG repeats were found to have maximum occurrence among di-nucleotide repeats in all the six genotypes, while AAT/ATT repeats were found to be most frequent among tri-nucleotide repeats.

Table 4

Chromosome-wise distribution of SSRs in Vigna cultivars VC1973A, Jingnong6 and Kyungwonpat

V. radiataV. angularisV. angularis
VC1973ADensity per
Megabase
Jingnong 6Density per
Megabase
KyungwonpatDensity per
Megabase
Chr 115 207416.614 964354.310 229339.1
Chr 211 317446.211 948335.38989317.9
Chr 37424573.213 679331.98383336.8
Chr 49522457.510 632306.69097332.3
Chr 518 212489.89411306.64771431.8
Chr 615 399411.37391381.813 887370.0
Chr 726 280472.612 016304.93976422.0
Chr 820 521448.78976291.54156321.7
Chr 99144435.210 6053345007377.6
Chr 109966474.612 112299.63102401.8
Chr 1197684957808293.610 407421.5
V. radiataV. angularisV. angularis
VC1973ADensity per
Megabase
Jingnong 6Density per
Megabase
KyungwonpatDensity per
Megabase
Chr 115 207416.614 964354.310 229339.1
Chr 211 317446.211 948335.38989317.9
Chr 37424573.213 679331.98383336.8
Chr 49522457.510 632306.69097332.3
Chr 518 212489.89411306.64771431.8
Chr 615 399411.37391381.813 887370.0
Chr 726 280472.612 016304.93976422.0
Chr 820 521448.78976291.54156321.7
Chr 99144435.210 6053345007377.6
Chr 109966474.612 112299.63102401.8
Chr 1197684957808293.610 407421.5
Table 4

Chromosome-wise distribution of SSRs in Vigna cultivars VC1973A, Jingnong6 and Kyungwonpat

V. radiataV. angularisV. angularis
VC1973ADensity per
Megabase
Jingnong 6Density per
Megabase
KyungwonpatDensity per
Megabase
Chr 115 207416.614 964354.310 229339.1
Chr 211 317446.211 948335.38989317.9
Chr 37424573.213 679331.98383336.8
Chr 49522457.510 632306.69097332.3
Chr 518 212489.89411306.64771431.8
Chr 615 399411.37391381.813 887370.0
Chr 726 280472.612 016304.93976422.0
Chr 820 521448.78976291.54156321.7
Chr 99144435.210 6053345007377.6
Chr 109966474.612 112299.63102401.8
Chr 1197684957808293.610 407421.5
V. radiataV. angularisV. angularis
VC1973ADensity per
Megabase
Jingnong 6Density per
Megabase
KyungwonpatDensity per
Megabase
Chr 115 207416.614 964354.310 229339.1
Chr 211 317446.211 948335.38989317.9
Chr 37424573.213 679331.98383336.8
Chr 49522457.510 632306.69097332.3
Chr 518 212489.89411306.64771431.8
Chr 615 399411.37391381.813 887370.0
Chr 726 280472.612 016304.93976422.0
Chr 820 521448.78976291.54156321.7
Chr 99144435.210 6053345007377.6
Chr 109966474.612 112299.63102401.8
Chr 1197684957808293.610 407421.5

In silico discovery of polymorphism and cross-species transferability

Species-wise data on SSR loci polymorphism was generated successfully by e-PCR using genomic data of respective varieties. A total of 164 polymorphic loci were obtained between genotypes VC1973A and RIL59 of V. radiata. Similarly, in species V. angularis, the two sets of comparison between genotypes Jingnong 6 and Kyungwonpat and Jingnong 6 and IT213134 revealed 243 and 195 polymorphic loci, respectively (Supplementary Table 1).

Cross-species transferability of SSR loci was evaluated using V. radiata cv. VC1973A as focal species over two non-focal species, V. angularis and V. unguiculata. Maximum transferability across species was found in non-focal species V. angularis (2202) than V. unguiculata (2373). In order to obtain more SSR loci having cross-species transferability, the second species (V. angularis) was treated as focal species against non-focal species (V. unguiculata), which revealed 5882 loci (Supplementary Table 2). Since flanking region of SSR loci are evolutionary well conserved, computational prediction of such loci have practical significance. Such SSR loci can cater the need of molecular markers even in the absence of whole-genome sequence data. Present finding of cross-species transferable loci can be prioritized for use in diversity and population differentiation analysis of other species of genus Vigna, which are without whole-genome sequence data. Since more than 100 species of genus Vigna are available but genomic sequence for marker development is confined to just three species, present finding is of immense values.

Revalidation of microsatellite primers generated by VigSatDB

PCR genotyping of 15 microsatellite loci in a panel 10 genotype was successfully obtained, demonstrating the efficacy of VigSatDB. Out of the 15 microsatellite loci, successful amplification was obtained for 14 loci. Figure 3 represents one of the representative gels showing amplification of microsatellite loci along with polymorphism. One locus could not be amplified, which may be due to PCR standardization method (Supplementary Table 3).

Microsatellite polymorphism detection in 10 genotypes of V. radiata in QIAxcel system.
Figure 3

Microsatellite polymorphism detection in 10 genotypes of V. radiata in QIAxcel system.

Interface of VigSatDB for various searches.
Figure 4

Interface of VigSatDB for various searches.

Web genomic resources: VigSatDB

VigSatDB has been successfully developed using three-tier architecture. It contains a total of 876 219 SSR loci mined from whole-genome sequence data of three Vigna species represented by six different varieties. It also contains 1976 published SSR loci, which are already validated. This database has seven tabs viz. Home, Microsatellite, Tools, Polymorphism, Statistics, Tutorial and Team. The `Microsatellite’ tab provides two search options viz., in silico mined SSRs and experimentally validated SSRs. In silico mined SSR tab is used for mining of simple and compound SSRs of chromosome-wise or scaffold-wise genome assemblies along with motif type, motif kind, GC content, percentage, copy number and start position. Primer3core executable integrated at the backend can generate locus-specific primers of desired loci displaying their respective values of annealing temperature, product size and GC content. In order to use genic region SSRs as putative FDM, user has the choice to mine genic region SSRs selectively along with gene annotation. Using genome browser, user can get the details of gene name, accession number, descriptions, positions, strands, etc.

Experimentally validated SSR search tab exhibits the pre-existing validated markers from the literature. This can be searched by polymorphic, motif type, repeat motif and repeat kind. The results displayed contains details of its species, accessions, motifs, type, repeat kind, forward and reverse primers, annealing temperature along with reference of literature. Under the `Tool’ tab, provision has been made for SSR search and Blast search (Figure 4). The `Polymorphism’ tab provides two tabs viz., `Self Design Primer’ (for mining of SSRs along with primer generation) and `External Primers’ (for evaluation of published primers). The `Self Design Primer’ tab generates the new primers, and users can search polymorphic markers against six genomes. For this search, all six reference genomes of Vigna were already incorporated in the database, and users have to select the respective option and the input genotype file is in fasta format. For the best possible results, users can upload input genotype in chromosome-wise separately. The outputs of polymorphic markers are displayed on the screen or get the automated result mails. The option of external primer evaluation can be used advantageously to compare the published data with newly generated SSR data of varieties under investigation (Figure 5).

(A) In silico mining of SSRs and primer designing; (B)e-PCR evaluation of published primers using specific variety of Vigna species.
Figure 5

(A) In silico mining of SSRs and primer designing; (B)e-PCR evaluation of published primers using specific variety of Vigna species.

Discussion

This first comprehensive microsatellite database of genus Vigna represented by three species co1tains 876 219 in silico predicted simple and compound markers and 1976 experimentally validated markers. This work clearly demonstrates that by in silico approach, mining of SSR can be done more efficiently (39). Since these markers are ubiquitously distributed over entire chromosomes, it can be a better representative in form of a molecular marker for genome variability analysis (40). Such approach has the advantage of location specificity over chromosome, thus it can be used as molecular marker of specific genes.

Flanking region SSR markers of targeted genes can be used in germplasm improvement program by MAS or introgression (41). Higher abundance of SSRs with respect to compound repeats was observed in all the genotypes. This is due to higher mutation rate of SSRs over compound repeat (42). All the six genotypes showed high abundance of mono-repeat, which might be due to inherent limitation of the Next Generation Sequencing (NGS) chemistry used in data generation (43). Similar higher abundance of di-nucleotide is also reported in other crops (44, 45). SSR relative density across all three species and their genotypes was found in 291–573 markers per Megabase, which is within the range of other reported crops like rice (370–490 markers per Megabase) and sugarbeet (341–671 markers per Megabase) (46).

A total of 200 178, 147 634 and 92 189 SSRs were mined from representative genotypes of three Vigna species such as V. radiata cv. VC1973A, V. angularis cv. Jingnong6 and V. unguiculata cv. IT97K-499-35, respectively. In species V. radiata cv. VC1973A, mono-nucleotides were the highest (127 648) and hexa-nucleotides were the lowest (177) in number. In V. angularis cv. Jingnong6, maximum numbers of SSRs were mono-type (86 204), while minimum number of SSR markers belonged to hexa-type (247). In V. unguiculata, mono-nucleotides were highest (59 935), while hexa-nucleotides were lowest (122) in number. In all the three species, mono-nucleotide repeats were most abundant followed by di- and tri-nucleotide repeats (Table 3).

Present finding of more than 875 000 putative SSR markers can be used for linkage mapping, identification of QTL/candidate genes and hybrid cultivar development. Similar use has already been reported in Vigna crop species for very limited few traits. For example, discovery of cha gene in Vigna is responsible for chasmogamous flower trait. This is highly relevant for molecular breeding to produce high yield in hybrid cultivar of Vigna (47). Similarly, Vigna SSR has also been used for PCR-based confirmation of successful production of interspecific recombinants of two species namely, mungbean (V. radiata) and mashbean (V. mungo) (48). SSR has been successfully used to identify tightly linked yellow mosaic virus resistance gene (49).

Variety characterization is based on phenotypic observation, which is very challenging to differentiate varieties with similar morphological characteristics. To overcome this traditional method, SSRs play an important role in variety identification, variety development, trait improvement, mapping, genetic studies and molecular breeding applications, product traceability, phylogenetic and taxonomic comparison (11, 50, 51). Very recently, 38 SSR markers of black gram (V. mungo) has been used to differentiate more than 30 varieties with high precision by rapid genotyping in DUS testing program of germplasm management (52). Similarly, 44 SSR markers have been used to differentiate five varieties of mungbean (V. radiata) (53). Using five set of SSR markers, 42 varieties (6 HYVs and 36 landraces) have been differentiated in mungbean (V. radiata) (22). Putative SSR markers present in VigSatDB can be a valuable genomic resource for similar varietal differentiation as ready to use primers for rapid genotyping are available along with in vitro validated SSR markers. Similar approach for varietal differentiation has already been reported in many other crops like sugarcane (54), sesame (55), capsicum (56), tea (57), barley (58) and eggplant (59).

Present set of markers can further be used for trait improvement against biotic and abiotic stresses. For example, in halophytic Vigna marina, SSR has been used to discover one major QTL controlling salt tolerance, and such QTL can be transferred in breeding program using flanking region SSR markers having allele size difference between donor and recipient varieties (60). These markers can be used for high-density linkage mapping and gene discovery required for specific trait improvement.

A large number of whole-genome assemblies of different plant species in the public domain are available, which provides the opportunity to study cross-species transferability in closely related species. This can be helpful in cloning candidate genes from different species and orthologous loci among different species (61). Also, in silico prediction of polymorphic markers can be useful in phylogenetic as well as diversity studies. VigSatDB putative markers can be used for germplasm improvement. Similar use in Vigna species has been reported with very limited markers (14). PMDBase is a comprehensive microsatellite database of 110 plant species, but has limitations that have been illustrated in Table 5. Our database can supplement the need of genomic resource for evolution and diversity analysis as reported in black gram (17).

Table 5

Comparison of PMDBase and VigSatDB databases

AttributesPMDBaseVigSatDB
Number of Vigna species covered23
Polymorphism prediction+
Chromosome-wise search option+
Motif-type search option+
Repeat motif+
Repeat kind+
Genic and non-genic marker search by users’ choice+
GC% content search option+
Chromosomal location search option+
Copy no. search option+
Non-nuclear (mitochondrion and chloroplast)+
Flexibility in primer designing+
ePCR option for evaluation of published primers+
AttributesPMDBaseVigSatDB
Number of Vigna species covered23
Polymorphism prediction+
Chromosome-wise search option+
Motif-type search option+
Repeat motif+
Repeat kind+
Genic and non-genic marker search by users’ choice+
GC% content search option+
Chromosomal location search option+
Copy no. search option+
Non-nuclear (mitochondrion and chloroplast)+
Flexibility in primer designing+
ePCR option for evaluation of published primers+
Table 5

Comparison of PMDBase and VigSatDB databases

AttributesPMDBaseVigSatDB
Number of Vigna species covered23
Polymorphism prediction+
Chromosome-wise search option+
Motif-type search option+
Repeat motif+
Repeat kind+
Genic and non-genic marker search by users’ choice+
GC% content search option+
Chromosomal location search option+
Copy no. search option+
Non-nuclear (mitochondrion and chloroplast)+
Flexibility in primer designing+
ePCR option for evaluation of published primers+
AttributesPMDBaseVigSatDB
Number of Vigna species covered23
Polymorphism prediction+
Chromosome-wise search option+
Motif-type search option+
Repeat motif+
Repeat kind+
Genic and non-genic marker search by users’ choice+
GC% content search option+
Chromosomal location search option+
Copy no. search option+
Non-nuclear (mitochondrion and chloroplast)+
Flexibility in primer designing+
ePCR option for evaluation of published primers+

Since genomes of most of Vigna sp. except V. radiata, V. unguiculata and V. angularis are yet to be sequenced, VigSatDB can be used as a genomic resource to cater the need of molecular markers required for QTL, linkage and genome mapping as well as germplasm improvement programs. Our finding of polymorphic SSRs and cross-species highly conserved SSRs can be used to characterize >90% species of genus Vigna. Such high transferability and their use in breeding, genetic, diversity and genomic studies are reported in crops like carrot species of Apiaceae family (62).

Although there are several known limitations in the utilization of cross-species SSR or heterologous SSR in diversity analysis (63), it is interesting to note in genus Vigna that the extent of SSR conservation is so high that even across genera (Glycine max and Phaseolus vulgaris) such conservation is reported in tune of more than one-fifth. Cross-species transferability of SSR loci in genus Vigna (V. angularis) is reported in tune of more than two-thirds in species like V. umbellata, V. mungo, V. aconitifolia and V. radiata (64). This extent of conservation has been in tune of more than three-fourths in V. aconitifolia and V. reflexo-pilosa (65). Higher rate of cross-species transferability of SSR loci in Vigna is due to their higher homology at genomic level; thus, they are much more useful in comparative genomics and genome mapping (8). Very safely it can be concluded that cross-species transferability of SSR loci is likely to be much more successful especially in crop like Vigna.

Conclusion

We report here first comprehensive web genomic resource of genus Vigna covering three of its commercially important species. VigSatDB contains a total of 876 219 putative microsatellite DNA markers, which also includes 1976 validated published markers. A total of 10 457 polymorphic markers obtained by e-PCR across genotypes can be used for germplasm differentiation and diversity analysis saving in vitro cost and time required for in vitro polymorphic discovery. Our finding on cross species transferability of microsatellite locus among three different species of Vigna can be used to cater the need of molecular markers especially for more than 100 species of genus Vigna where there is no whole-genome sequence data even today. This genomic resource can be of immense use for the global community. It can be used for chromosome-wise microsatellite locus mining and primer designing for non-genic and genic FDM-SSR for rapid genotyping. It can also be used to accelerate polymorphism discovery by e-PCR most economically especially needed in future re-sequencing projects. It can be used not only for knowledge discovery research like QTL and gene mapping but also for marker-assisted breeding in Vigna germplasm improvement and management.

Author’s contribution

P.K.Y., D.K., M.A.I., A.R., S.J. and R.S.J. conceived the theme of the study. R.S.J., M.A.I. and S.J. did the computational analysis of generated data. R.S.J., U.B.A. and V.A. contributed to database development. S.B.B. and R.S.T. did in vitro validation of markers. R.S.J., M.A.I. and S.J. drafted the manuscript. D.K., S.J., P.K.Y., M.A.I. and A.R. edited the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We are thankful to Director, ICAR-IASRI, Indian Council of Agricultural Research, Ministry of Agriculture and Farmers' Welfare, Government of India for providing financial and infrastructural support and ASHOKA facility at ICAR-IASRI, New Delhi, for carrying out this research. We extend our sincere thanks to the Department of Computational Biology and Bioinformatics, Sam Higginbottom University of Agriculture, Technology and Sciences, India, for providing facilities and support.

Conflict of interest. None declared.

Database URL:http://webtom.cabgrid.res.in/vigna_ssr/

References

1.

Datta
,
S.
and
Gupta
,
S.
(
2009
)
Use of genomic resources in improvement of Vigna species
.
J. Food. Legum.
,
22
,
1
10
.

2.

Takahashi
,
Y.
,
Somta
,
P.
,
Muto
,
C.
et al.  (
2016
)
Novel genetic resources in the genus Vigna unveiled from gene bank accessions
.
PLoS One
,
11
, e0147568.

3.

Food and Agriculture Organization of the United Nations
(
2016
) FAO Statistical Yearbook: Rome World Food and Agriculture. http://www.fao.org/faostat/en/.

4.

Javadi
,
F.
,
Tun
,
Y.T.
,
Kawase
,
M.
et al.  (
2011
)
Molecular phylogeny of the subgenus Ceratotropis (genus Vigna, Leguminosae) reveals three eco-geographical groups and late Pliocene–Pleistocene diversification: evidence from four plastid DNA region sequences
.
Ann. Bot.
,
108
,
367
380
.

5.

Harouna
,
D.V.
,
Venkataramana
,
P.B.
,
Ndakidemi
,
P.A.
et al.  (
2018
)
Under-exploited wild Vigna species potentials in human and animal nutrition: a review
.
Glob. Food Sec.
,
18
,
1
11
.

6.

Takahashi
,
Y.
,
Iseki
,
K.
,
Kitazawa
,
K.
et al.  (
2015
)
A homoploid hybrid between wild Vigna species found in a limestone karst
.
Front. Plant Sci.
,
6
,
1050
.

7.

Ganesan
,
K.
and
Xu
,
B.
(
2017
)
A critical review on phytochemical profile and health promoting effects of mung bean (Vigna radiata)
.
Food Science and Human Wellness
,
7
,
11
33
.

8.

Tangphatsornruang
,
S.
,
Somta
,
P.
,
Uthaipaisanwong
,
P.
et al.  (
2009
)
Characterization of microsatellites and gene contents from genome shotgun sequences of mungbean (Vigna radiata (L.) Wilczek)
.
BMC Plant Biol.
,
9
,
137
.

9.

Kang
,
Y.J.
,
Kim
,
S.K.
,
Kim
,
M.Y.
et al.  (
2014
)
Genome sequence of mungbean and insights into evolution within Vigna species
.
Nat. Commun.
,
5
,
ncomms6443
.

10.

Kapil
,
A.
,
Rai
,
P.K.
and
Shanker
,
A.
(
2014
)
ChloroSSRdb: a repository of perfect and imperfect chloroplastic simple sequence repeats (cpSSRs) of green plants
.
Database (Oxford)
,
bau107
,
doi:
.

11.

Yu
,
J.
,
Dossa
,
K.
,
Wang
,
L.
et al.  (
2016
)
PMDBase: a database for studying microsatellite DNA and marker development in plants
.
Nucleic Acids Res.
,
45
,
D1046
D1053
.

12.

Wang
,
L.
,
Do Kim
,
K.
,
Gao
,
D.
et al.  (
2016
)
Analysis of simple sequence repeats in rice bean (Vigna umbellata) using an SSR-enriched library
.
Crop J.
,
4
,
40
47
.

13.

Wang
,
L.X.
,
Elbaidouri
,
M.
,
Abernathy
,
B.
et al.  (
2015
)
Distribution and analysis of SSR in mung bean (Vigna radiata L.) genome based on an SSR-enriched library
.
Mol. Breed.
,
35
,
25
.

14.

Chen
,
H.
,
Liu
,
L.
,
Wang
,
L.
et al.  (
2015
)
Development and validation of EST-SSR markers from the transcriptome of adzuki bean (Vigna angularis)
.
PLoS One
,
10
, e0131939.

15.

Gupta
,
S.K.
,
Bansal
,
R.
and
Gopalakrishna
,
T.
(
2014
)
Development and characterization of genic SSR markers for mungbean (Vigna radiata (L.) Wilczek)
.
Euphytica
,
195
,
245
258
.

16.

Chen
,
H.
,
Wang
,
L.
,
Wang
,
S.
et al.  (
2015
)
Transcriptome sequencing of mung bean (Vigna radiate L.) genes and the identification of EST-SSR markers
.
PLoS One
,
10
, e0120273.

17.

Souframanien
,
J.
and
Reddy
,
K.S.
(
2015
)
De novo assembly, characterization of immature seed transcriptome and development of genic-SSR markers in black gram [Vigna mungo (L.) Hepper]
.
PLoS One
,
10
, e0128748.

18.

Sharma
,
P.C.
,
Grover
,
A.
and
Kahl
,
G.
(
2007
)
Mining microsatellites in eukaryotic genomes
.
Trends Biotechnol.
,
25
,
490
498
.

19.

Gupta
,
S.K.
,
Souframanien
,
J.
and
Gopalakrishna
,
T.
(
2008
)
Construction of a genetic linkage map of black gram, Vigna mungo (L.) Hepper, based on molecular markers and comparative studies
.
Genome
,
51
,
628
637
.

20.

Marubodee
,
R.
,
Ogiso-Tanaka
,
E.
,
Isemura
,
T.
et al.  (
2015
)
Construction of an SSR and RAD-marker based molecular linkage map of Vigna vexillata (L.) a. Rich
.
PLoS One
,
10
, e0138942.

21.

Kwon
,
Y.S.
,
Lee
,
J.M.
,
Yi
,
G.B.
et al.  (
2005
)
Use of SSR markers to complement tests of distinctiveness, uniformity, and stability (DUS) of pepper (Capsicum annuum L.) varieties
.
Mol. Cells
,
19
,
428
435
.

22.

Molla
,
M.R.
,
Ahmed
,
I.
,
Rohman
,
M.M.
et al.  (
2016
)
Genetic diversity analysis and DNA fingerprinting of Mungbean (Vigna radiata L.) genotypes using SSR markers
.
J. Plant Sci.
,
6
,
153
164
.

23.

Moore
,
S.S.
,
Sargeant
,
L.L.
,
King
,
T.J.
et al.  (
1991
)
The conservation of dinucleotide microsatellites among mammalian genomes allows the use of heterologous PCR primer pairs in closely related species
.
Genomics
,
10
,
654
660
.

24.

Gao
,
L.Z.
,
Zhang
,
C.H.
and
Jia
,
J.Z.
(
2005
)
Cross-species transferability of rice microsatellites in its wild relatives and the potential for conservation genetic studies
.
Genet. Resour. Crop Evol.
,
52
,
931
940
.

25.

Kitsanachandee
,
R.
,
Somta
,
P.
,
Chatchawankanphanich
,
O.
et al.  (
2013
)
Detection of quantitative trait loci for mungbean yellow mosaic India virus (MYMIV) resistance in mungbean (Vigna radiata (L.) Wilczek) in India and Pakistan
.
Breed Sci.
,
63
,
367
373
.

26.

Chen
,
H.
,
Chen
,
X.
,
Tian
,
J.
et al.  (
2016
)
Development of gene-based SSR markers in rice bean (Vigna umbellata L.) based on transcriptome data
.
PLoS One
,
11
, e0151040.

27.

Chen
,
H.
,
Wang
,
L.
,
Liu
,
X.
et al.  (
2017
)
De novo transcriptomic analysis of cowpea (Vigna unguiculata L. Walp.) for genic SSR marker development
.
BMC Genet.
,
18
,
65
.

28.

Liu
,
M.S.
,
Kuo
,
T.C.Y.
,
Ko
,
C.Y.
et al.  (
2016
)
Genomic and transcriptomic comparison of nucleotide variations for insights into bruchid resistance of mungbean (Vigna radiata [L.] R. Wilczek)
.
BMC Plant Biol.
,
16
,
46
.

29.

Yang
,
K.
,
Tian
,
Z.
,
Chen
,
C.
et al.  (
2015
)
Genome sequencing of adzuki bean (Vigna angularis) provides insight into high starch and low fat accumulation and domestication
.
Proc. Natl. Acad. Sci. U. S. A.
,
112
,
13213
13218
.

30.

Kang
,
Y.J.
,
Satyawan
,
D.
,
Shim
,
S.
et al.  (
2015
)
Draft genome sequence of adzuki bean, Vigna angularis
.
Sci. Rep.
,
5
,
8069
.

31.

Lestari
,
P.
,
Kang
,
Y.J.
,
Han
,
K.S.
et al.  (
2014
)
Genome-wide single nucleotide polymorphism discovery and validation in adzuki bean
.
Mol Breed.
,
33
,
497
501
.

32.

Muñoz-Amatriaín
,
M.
,
Mirebrahim
,
H.
,
Xu
,
P.
et al.  (
2017
)
Genome resources for climate-resilient cowpea, an essential crop for food security
.
Plant J.
,
89
,
1042
1054
.

33.

Thiel
,
T.
,
Michalek
,
W.
,
Varshney
,
R.
et al.  (
2003
)
Exploiting EST databases for the development and characterization of gene-derived SSR-markers in barley (Hordeum vulgare L.)
.
Theor. Appl. Genet.
,
106
,
411
422
.

34.

Untergasser
,
A.
,
Cutcutache
,
I.
,
Koressaar
,
T.
et al.  (
2012
)
Primer3—new capabilities and interfaces
.
Nucleic Acids Res.
,
40
,
e115
e115
.

35.

Schlötterer
,
C.
and
Tautz
,
D.
(
1992
)
Slippage synthesis of simple sequence DNA
.
Nucleic Acids Res.
,
20
,
211
215
.

36.

Altschul
,
S.F.
,
Gish
,
W.
,
Miller
,
W.
et al.  (
1990
)
Basic local alignment search tool
.
J. Mol. Biol.
,
215
,
403
410
.

37.

Stein
,
L.D.
,
Mungall
,
C.
,
Shu
,
S.
et al.  (
2002
)
The generic genome browser: a building block for a model organism system database
.
Genome Res.
,
12
,
1599
1610
.

38.

Doyle
,
J.J.
and
Doyle
,
J.L.
(
1990
)
Isolation of plant DNA from fresh tissue
.
Focus
,
12
,
13
15
.

39.

Xiao
,
Y.
,
Xia
,
W.
,
Ma
,
J.
et al.  (
2016
)
Genome-wide identification and transferability of microsatellite markers between Palmae species
.
Front. Plant Sci.
,
7
,
1578
.

40.

Guo
,
W.J.
,
Ling
,
J.
and
Li
,
P.
(
2009
)
Consensus features of microsatellite distribution: microsatellite contents are universally correlated with recombination rates and are preferentially depressed by centromeres in multicellular eukaryotic genomes
.
Genomics
,
93
,
323
331
.

41.

Collard
,
B.C.
and
Mackill
,
D.J.
(
2008
)
Marker-assisted selection: an approach for precision plant breeding in the twenty-first century
.
Philos. Trans. R. Soc. Lond. B Biol. Sci.
,
363
,
557
572
.

42.

Ngai
,
M.Y.
and
Saitou
,
N.
(
2016
)
The effect of perfection status on mutation rates of microsatellites in primates
.
Anthropol. Sci.
,
124
,
85
92
.

43.

Haseneyer
,
G.
,
Schmutzer
,
T.
and
Seidel
,
M.
(
2011
)
From RNA-seq to large-scale genotyping-genomics resources for rye (Secale cereale L.)
.
BMC Plant Biol.
,
11
,
131
.

44.

Kariin
,
S.
and
Burge
,
C.
(
1995
)
Dinucleotide relative abundance extremes: a genomic signature
.
Trends Genet.
,
11
,
283
290
.

45.

Shioiri
,
C.
and
Takahata
,
N.
(
2001
)
Skew of mononucleotide frequencies, relative abundance of dinucleotides, and DNA strand asymmetry
.
J. Mol. Evol.
,
53
,
364
376
.

46.

Iquebal
,
M.A.
,
Jaiswal
,
S.
,
Angadi
,
U.B.
et al.  (
2015
)
SBMDb: first whole genome putative microsatellite DNA marker database of sugarbeet for bioenergy and industrial applications
.
Database (Oxford)
,
bav111
,
doi:
.

47.

Chen
,
J.
,
Somta
,
P.
,
Chen
,
X.
et al.  (
2016
)
Gene mapping of a mutant mungbean (Vigna radiata L.) using new molecular markers suggests a gene encoding a yuc4-like protein regulates the chasmogamous flower trait
.
Front Plant Sci.
,
7
,
830
.

48.

Abbas
,
G.
,
Hameed
,
A.
and
Rizwan
,
M.
(
2015
)
Genetic confirmation of mungbean (Vigna radiata) and mashbean (Vigna mungo) interspecific recombinants using molecular markers
.
Front Plant Sci.
,
6
,
1107
.

49.

Rambabu
,
E.
,
Anuradha
,
C.
,
Sridhar
,
V.
et al.  (
2018
)
Identification of molecular markers linked to yellow mosaic virus resistance in blackgram (Vigna mungo (L.) Hepper)
.
Int. J. Curr. Microbiol. App. Sci.
,
7
,
3810
3817
.

50.

Iquebal
,
M.A.
,
Sarika
,
Arora
,
V.
et al.  (
2013
)
First whole genome based microsatellite DNA marker database of tomato for mapping and variety identification
.
BMC Plant Biol.
,
13
,
197
.

51.

Zietkiewicz
,
E.
,
Rafalski
,
A.
and
Labuda
,
D.
(
1994
)
Genome fingerprinting by simple sequence repeat (SSR)-anchored polymerase chain reaction amplification
.
Genomics
,
20
,
176
183
.

52.

Prathyusha
,
V.B.
,
Rao
,
Y.K.
,
Rao
,
G.P.
et al.  (
2017
)
SSR polymorphism and distinctiveness, uniformity and stability (DUS) features in selected genotypes of blackgram (Vigna mungo L. Hepper)
.
J. Food Legum.
,
30
,
1
8
.

53.

Kaur
,
R.
,
Toor
,
A.K.
,
Bassi
,
G.
et al.  (
2017
)
Characterization of Mungbean (Vigna radiata L. Wilczek) varieties using morphological and molecular descriptors
.
Int. J. Curr. Microbiol. App. Sci.
,
6
,
1609
1618
.

54.

Manigbas
,
N.L.
and
Villegas
,
L.C.
(
2004
)
Microsatellite markers in hybridity tests to identify true hybrids of sugarcane
.
Philipp. J. Crop Sci.
,
29
,
23
32
.

55.

Bhattacharyya
,
U.
,
Pandey
,
S.K.
and
Dasgupta
,
T.
(
2014
)
Identification of EST-SSRs and FDM in sesame (Sesamum indicum L.) through data mining
.
Sch. J. Agric. Sci.
,
4
,
60
69
.

56.

Shirasawa
,
K.
,
Ishii
,
K.
,
Kim
,
C.
et al.  (
2013
)
Development of capsicum EST–SSR markers for species identification and in silico mapping onto the tomato genome sequence
.
Mol. Breed.
,
31
,
101
110
.

57.

Sahu
,
J.
,
Sarmah
,
R.
,
Dehury
,
B.
et al.  (
2012
)
Mining for SSRs and FDMs from expressed sequence tags of Camellia sinensis
.
Bioinformation
,
8
,
260
.

58.

Karakousis
,
A.
,
Barr
,
A.R.
,
Chalmers
,
K.J.
et al.  (
2003
)
Potential of SSR markers for plant breeding and variety identification in Australian barley germplasm
.
Aust. J. Agric. Res.
,
54
,
1197
1210
.

59.

Stàgel
,
A.
,
Portis
,
E.
,
Toppino
,
L.
et al.  (
2008
)
Gene-based microsatellite development for mapping and phylogeny studies in eggplant
.
BMC Genomics
,
9
,
357
.

60.

Chankaew
,
S.
,
Isemura
,
T.
,
Naito
,
K.
et al.  (
2014
)
QTL mapping for salt tolerance and domestication-related traits in Vigna marina subsp. oblonga, a halophytic species
.
Theor. Appl. Genet.
,
127
,
691
702
.

61.

Biswas
,
M.K.
,
Liu
,
Y.
,
Li
,
C.
et al.  (
2015
)
Genome-wide computational analysis of musa microsatellites: classification, cross-taxon transferability, functional annotation, association with transposons & miRNAs, and genetic marker potential
.
PLoS One
,
10
, e0131312.

62.

Cavagnaro
,
P.F.
,
Chung
,
S.M.
,
Manin
,
S.
et al.  (
2011
)
Microsatellite isolation and marker development in carrot-genomic distribution, linkage mapping, genetic diversity analysis and marker transferability across Apiaceae
.
BMC Genomics
,
12
,
386
.

63.

Barbara
,
T.
,
Palma-Silva
,
C.
,
Paggi
,
G.M.
et al.  (
2007
)
Cross-species transfer of nuclear microsatellite markers: potential and limitations
.
Mol. Ecol.
,
16
,
3759
3767
.

64.

Chaitieng
,
B.
,
Kaga
,
A.
,
Tomooka
,
N.
et al.  (
2006
)
Development of a black gram [Vigna mungo (L.) Hepper] linkage map and its comparison with an azuki bean [Vigna angularis (Willd.) Ohwi and Ohashi] linkage map
.
Theor. Appl. Genet.
,
113
,
1261
1269
.

65.

Somta
,
P.
,
Seehalak
,
W.
and
Srinives
,
P.
(
2009
)
Development, characterization and cross-species amplification of mungbean (Vigna radiata) genic microsatellite markers
.
Conserv. Genet.
,
10
,
1939
.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.