Characterization of microsatellites and gene contents from genome shotgun sequences of mungbean (Vigna radiata (L.) Wilczek)

Genetics, genomics, and breeding of black gram [Vigna mungo (L.) Hepper]

Black gram [Vigna mungo (L.) Hepper] is a highly nutritious grain legume crop, mainly grown in South and Southeast Asia, with the largest area in India, where the crop is challenged by several biotic and abiotic stresses leading to significant yield losses. Improving genetic gains to increase on-farm yields is the primary goal of black gram breeding programs. This could be achieved by developing varieties resistant to major diseases like mungbean yellow mosaic disease, urdbean leaf crinkle virus, Cercospora leaf spot, anthracnose, powdery mildew, and insect pests such as whitefly, cowpea aphids, thrips, stem flies, and bruchids. Along with increasing on-farm yields, incorporating market-preferred traits ensures the adoption of improved varieties. Black gram breeding programs rely upon a limited number of parental lines, leading to a narrow genetic base of the developed varieties. For accelerating genetic gain, there is an urgent need to include more diverse genetic material for improving traits for better adaptability and stress resistance in breeding populations. The present review summarizes the importance of black gram, the major biotic and abiotic stresses, available genetic and genomic resources, major traits for potential crop improvement, their inheritance, and the breeding approaches being used in black gram for the development of new varieties.

Insights into morphological and physio-biochemical adaptive responses in mungbean (Vigna radiata L.) under heat stress

Mungbean (Vigna radiata L. Wilczek) is an important food legume crop which contributes significantly to nutritional and food security of South and Southeast Asia. The crop thrives in hot and humid weather conditions, with an optimal temperature range of 28°-35°C, and is mainly cultivated under rainfed environments. However, the rising global temperature has posed a serious threat to mungbean cultivation. Optimal temperature is a vital factor in cellular processes, and every crop species has evolved with its specific temperature tolerance ability. Moreover, variation within a crop species is inevitable, given the diverse environmental conditions under which it has evolved. For instance, various mungbean germplasm can grow and produce seeds in extreme ambient temperatures as low as 20°C or as high as 45°C. This range of variation in mungbean germplasm for heat tolerance plays a crucial role in developing heat tolerant and high yielding mungbean cultivars. However, heat tolerance is a complex mechanism which is extensively discussed in this manuscript; and at the same time individual genotypes have evolved with various ways of heat stress tolerance. Therefore, to enhance understanding towards such variability in mungbean germplasm, we studied morphological, anatomical, physiological, and biochemical traits which are responsive to heat stress in plants with more relevance to mungbean. Understanding heat stress tolerance attributing traits will help in identification of corresponding regulatory networks and associated genes, which will further help in devising suitable strategies to enhance heat tolerance in mungbean. The major pathways responsible for heat stress tolerance in plants are also discussed.

Small cardamom genome: development and utilization of microsatellite markers from a draft genome sequence of Elettaria cardamomum Maton

Small cardamom (Elettaria cardamomum Maton), the queen of spices, is the third most expensive spice in the world after saffron and vanilla, valued highly for its aroma and taste. This perennial herbaceous plant is a native of coastal parts of Southern India and displays a significant amount of morphological diversity. Its genetic potential has not been exploited due to lack of genomic resources limiting our understanding of the genome and important metabolic pathways which give it the economic advantage in the spice industry. Here, we report upon the de novo assembled, draft whole genome sequence of cardamom variety, Njallani Green Gold. We used a hybrid assembly strategy using the reads from the Oxford Nanopore, Illumina and 10x Genomics GemCode sequencing chemistries. The assembled genome length was 1.06 Gb (gigabases) which is close to the estimated genome size of cardamom. More than 75% of the genome was captured in 8000 scaffolds with a N50 of 0.15 Mb. The genome appears to have a high repeat content and 68055 gene models were predicted. The genome is close to Musa species and displays an expansion and contraction in different gene families. The draft assembly was used for in silico mining of simple sequence repeats (SSRs). A total of 2,50,571 SSRs were identified of which 2,18,270 were perfect SSRs and 32,301 were compound SSRs. Among the perfect SSRs, trinucleotides were most abundant (1,25,329) and hexanucleotide repeats appear least (2,380). From the 2,50,571 SSRs mined, 2,27,808 primer pairs were designed based on flanking sequence information. Wet lab validation was performed for 246 SSR loci and based on their amplification profiles, 60 SSR markers were used for diversity analysis of a set of 60 diverse cardamom accessions. The average number of alleles detected per locus were 14.57 with a minimum of 4 and maximum of 30 alleles. Population structure analysis revealed the presence of high degree of admixtures which could primarily be due to cross-pollination prevalent in this species. The SSR markers identified would help in the development of gene or trait-linked markers which can be subsequently used for marker-assisted breeding for crop improvement in cardamom. The information on utilization of the SSR loci for generation of markers has been developed into a public database, 'cardamomSSRdb' that is freely available for use by the cardamom community.

Genome-wide characterization and development of SSR markers for genetic diversity analysis in northwestern Himalayas Walnut (Juglans regia L.)

In the present study, we designed and validated genome-wide polymorphic SSR markers (110 SSRs) by mining the walnut genome. A total of 198,924 SSR loci were identified. Among these, successful primers were designed for 162,594 (81.73%) SSR loci. Dinucleotides were the most predominant accounting for 88.40% (175,075) of total SSRs. The SSR frequency was 377.312 SSR/Mb and it showed a decreasing trend from dinucleotide to octanucleotide motifs. We identified 20 highly polymorphic SSR markers and used them to genotype 72 walnut accessions. Over all, we obtained 118 alleles that ranged from 2 to 12 with an average value of 5.9. The higher SSR PIC values indicate their robustness in discriminating walnut genotypes. Heat map, PCA, and population structure categorized 72 walnut genotypes into 2 distinct clusters. The genetic variation within population was higher than among population as inferred by analysis of molecular variance (AMOVA). For walnut improvement, it is necessary to have a large repository of SSRs with high discriminative power. The present study reports 150,000 SSRs, which is the largest SSR repository for this important nut crop. Scientific communities may use this repository for walnut improvement such as QTL mapping, genetic studies, linkage map construction, and marker-assisted selection.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03563-6.

Vignette of Vigna domestication: From archives to genomics

The genus Vigna comprises fast-growing, diploid legumes, cultivated in tropical and subtropical parts of the world. It comprises more than 200 species among which Vigna angularis, Vigna radiata, Vigna mungo, Vigna aconitifolia, Vigna umbellata, Vigna unguiculata, and Vigna vexillata are of enormous agronomic importance. Human selection along with natural variability within these species encompasses a vital source for developing new varieties. The present review convokes the early domestication history of Vigna species based on archeological pieces of evidence and domestication-related traits (DRTs) together with genetics of domestication. Traces of early domestication of Vigna have been evidenced to spread across several temperate and tropical regions of Africa, Eastern Asia, and few parts of Europe. Several DRTs of Vigna species, such as pod shattering, pod and seed size, dormancy, seed coat, seed color, maturity, and pod dehiscence, can clearly differentiate wild species from their domesticates. With the advancement in next-generation high-throughput sequencing techniques, exploration of genetic variability using recently released reference genomes along with de novo sequencing of Vigna species have provided a framework to perform genome-wide association and functional studies to figure out different genes related to DRTs. In this review, genes and quantitative trait loci (QTLs) related to DRTs of different Vigna species have also been summarized. Information provided in this review will enhance the in-depth understanding of the selective pressures that causes crop domestication along with nature of evolutionary selection made in unexplored Vigna species. Furthermore, correlated archeological and domestication-related genetic evidence will facilitate Vigna species to be considered as suitable model plants.

Thirty Years of Mungbean Genome Research: Where Do We Stand and What Have We Learned?

Mungbean is a socioeconomically important legume crop in Asia that is currently in high demand by consumers and industries both as dried beans and in plant-based protein foods. Marker-assisted and genomics-assisted breeding are promising approaches to efficiently and rapidly develop new cultivars with improved yield, quality, and resistance to biotic and abiotic stresses. Although mungbean was at the forefront of research at the dawn of the plant genomics era 30 years ago, the crop is a "slow runner" in genome research due to limited genomic resources, especially DNA markers. Significant progress in mungbean genome research was achieved only within the last 10 years, notably after the release of the VC1973A draft reference genome constructed using next-generation sequencing technology, which enabled fast and efficient DNA marker development, gene mapping, and identification of candidate genes for complex traits. Resistance to biotic stresses has dominated mungbean genome research to date; however, research is on the rise. In this study, we provide an overview of the past progress and current status of mungbean genomics research. We also discuss and evaluate some research results to provide a better understanding of mungbean genomics.

Association Mapping for Yield Attributing Traits and Yellow Mosaic Disease Resistance in Mung Bean [Vigna radiata (L.) Wilczek]

Mung bean [Vigna radiata (L.) Wilczek] is an important short-duration grain legume widely known for its nutritional, soil ameliorative, and cropping system intensification properties. This study aims at evaluating genetic diversity among mung bean genotypes and detecting genomic regions associated with various yield attributing traits and yellow mosaic disease (YMD) resistance by association mapping. A panel of 80 cultivars and advanced breeding lines was evaluated for 10 yield-related and YMD resistance traits during kharif (monsoon) and summer seasons of 2018-2019 and 2019-2020. A total of 164 genome-wide simple sequence repeat (SSR) markers were initially screened, out of which 89 were found polymorphic which generated 317 polymorphic alleles with an average of 3.56 alleles per SSR locus. The number of alleles at each locus varied from 2 to 7. The population genetic structure analysis grouped different genotypes in three major clusters and three genetically distinct subpopulations (SPs) (i.e., SP-1, SP-2, and SP-3) with one admixture subpopulation (SP-4). Both cluster and population genetic structure analysis categorized the advanced mung bean genotypes in a single group/SP and the released varieties in other groups/SPs, suggesting that the studied genotypes may have common ancestral history at some level. The population genetic structure was also in agreement with the genetic diversity analysis. The estimate of the average degree of linkage disequilibrium (LD) present at the genome level in 80 mung bean genotypes unveiled significant LD blocks. Over the four seasons, 10 marker-trait associations were observed significant for YMD and four seed yield (SY)-related traits viz., days to flowering, days to maturity, plant height, and number of pods per plant using the mixed linear model (MLM) method. These associations may be useful for marker-assisted mung bean yield improvement programs and YMD resistance.

De novo assembly, transcriptome characterization and marker discovery in Indian major carp, Labeo rohita through pyrosequencing

Labeo rohita, one of the Indian major carps, is the most popular culture species in Indian subcontinent due to its consumer preference and delicacy. A selective breeding program for harvest body weight has resulted in an average genetic gain of 17% per generation. Transcriptome resource for this species is scanty. Here, we have characterized the liver and muscle transcriptomes of rohu using Roche 454 GS-FLX next generation sequencing platform. In total, 1.2 million reads were generated, de novo assembly and clustering resulted in 4171 transcripts. Out of these, 4171 had significant blast hit against NCBI nr database, and 2130 transcripts were successfully annotated. In total, 289 SSRs were identified with an identification rate of 5.8%, and dinucleotide repeat motifs were observed to be the most abundant SSRs. Further, 2231 putative SNPs were identified with high confidence. Validation of eight putative SNPs using Sanger sequencing resulted in 100% true SNPs. Significant allelic imbalance of M1, M4 and M5 loci between growth selected and control individual were observed. Furthermore, 13 transcription factors were identified in the present study belonging to six different transcription factor families. The present study demonstrated the utility of RNAseq to develop genomics resources in non-model fish species, and the marker resources developed would support the genetic improvement program of this species.

Genetic and phenotypic diversity of selected Kenyan mung bean (Vigna radiata L. Wilckzek) genotypes

BACKGROUND

Mung bean is a pulse crop principally grown in the tropic and subtropic parts of the world for its nutrient-rich seeds. Seven mung beans accessions from Eastern Kenya were evaluated using thirteen phenotypic traits. In addition, 10 SSR markers were used to determine their genetic diversity and population structure. This aimed at enhancing germplasm utilization for subsequent mung bean breeding programs.

RESULTS

Analysis of variance for most of the phenology traits showed significant variation, with the yield traits recording the highest. The first three principal components (PC) explained 83.4% of the overall phenotypic variation, with the highest (PC1) being due to variation of majority of the traits studied such as pod length, plant height, and seeds per pod. The dendogram revealed that the improved genotypes had common ancestry with the local landraces. The seven mung beans were also genotyped using 10 microsatellite markers, eight of which showed clear and consistent amplification profiles with scorable polymorphisms in all the studied genotypes. Genetic diversity, allele number, and polymorphic information content (PIC) were determined using powermarker (version 3.25) and phylogenetic tree constructed using DARWIN version 6.0.12. Analysis of molecular variance (AMOVA) was calculated using GenALEx version 6.5. A total of 23 alleles were detected from the seven genotypes on all the chromosomes studied with an average of 2.875 across the loci. The PIC values ranged from 0.1224 (CEDG056) to 0.5918 (CEDG092) with a mean of 0.3724. Among the markers, CEDG092 was highly informative while the rest were reasonably informative except CEDG056, which was less informative. Gene diversity ranged from 0.1836 (CEDG050) to 0.5102 (CDED088) with an average of 0.3534. The Jaccards dissimilarity matrix indicated that genotypes VC614850 and N26 had the highest level of dissimilarity while VC637245 and N26 had lowest dissimilarity index. The phylogenetic tree grouped the genotypes into three clusters as revealed by population structure analysis (K = 3), with cluster III having one unique genotype (VC6137B) only. AMOVA indicated that the highest variation (99%) was between individual genotype. In addition, marker traits association analysis revealed 18 significant associations (P < 0.05).

CONCLUSION

These findings indicate sufficient variation among the studied genotypes that can be considered for germplasm breeding programs.

Development, characterization, functional annotation and validation of genomic and genic-SSR markers using de novo next generation sequencing in Melia dubia Cav

Melia dubia Cav. (Meliaceae), a fast-growing tropical tree finds use in plywood, pulp and high-value solid wood products. To increase its productivity, we must essentially capture genetic diversity and identify genotypes with superior wood properties. This study aimed to develop novel microsatellite markers from genomic data and validate the markers in M. dubia. Direct Seq-to-SSR approach was adopted and using an in-house Perl script, 426,390 SSR markers identified. For validation, selected 151 markers, of which 50 were genomic markers chosen randomly, and 101 were genic markers identified through BLAST2GO. Amplification was observed in all loci, and 81.4% generated high-quality, reproducible amplicons of the expected size. Out of 50 genomic markers, we used ten highly polymorphic markers to assess genetic diversity among 75 genotypes from three populations. One hundred fourteen alleles were recorded, with a moderate level of diversity and a positive fixation index. Twenty-nine genic markers representing 13 enzymes showing polymorphism for wood stiffness were selected for diversity assessment of 24 genotypes (12 genotypes each with high and low-stress wave velocity). The product size ranged from 87 to 279, covering the majority of the genome. Cluster and structure analysis segregated ~ 80% of the genotypes based on the trait. This is the first report of the development of genic markers from a genomic survey and has proved efficient in differentiating genotypes based on the trait. The markers developed in this study will be useful for genetic mapping, diversity estimation, marker-assisted selection for desired traits and breeding for wood traits in M. dubia.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02858-w.

Genome Wide Characterization, Comparative and Genetic Diversity Analysis of Simple Sequence Repeats in Cucurbita Species

Simple sequence repeats (SSRs) are widely used in mapping constructions and comparative and genetic diversity analyses. Here, 103,056 SSR loci were found in Cucurbita species by in silico PCR. In general, the frequency of these SSRs decreased with the increase in the motif length, and di-nucleotide motifs were the most common type. For the same repeat types, the SSR frequency decreased sharply with the increase in the repeat number. The majority of the SSR loci were suitable for marker development (84.75% in Cucurbita moschata, 94.53% in Cucurbita maxima, and 95.09% in Cucurbita pepo). Using these markers, the cross-species transferable SSR markers between C. pepo and other Cucurbitaceae species were developed, and the complicated mosaic relationships among them were analyzed. Especially, the main syntenic relationships between C. pepo and C. moschata or C. maxima indicated that the chromosomes in the Cucurbita genomes were highly conserved during evolution. Furthermore, 66 core SSR markers were selected to measure the genetic diversity in 61 C. pepo germplasms, and they were divided into two groups by structure and unweighted pair group method with arithmetic analysis. These results will promote the utilization of SSRs in basic and applied research of Cucurbita species.

Genetic diversity and structure of landrace of lablab (Lablab purpureus (L.) Sweet) cultivars in Thailand revealed by SSR markers

Lablab (Lablab purpureus (L.) Sweet) is a legume crop widely cultivated in tropical and subtropical regions of Africa and Asia. In this study, we assessed genetic diversity and population structure of 299 individuals of subspecies purpureus and bengalensis of lablab from Thailand using 13 simple sequence repeat (SSR) markers. The SSR markers detected only 34 alleles in total with a mean of 2.6 alleles per locus. Overall gene diversity was 0.360. Gene diversity (H _E) and allelic richness (A _R) in different geographic regions was comparable. Similarly, both H _E and A _R between subspecies purpureus and bengalensis were similar. STRUCTURE and neighbor-joining (NJ) analyses revealed that the 299 individuals were clustered into two major groups. In contrast, principal coordinate analysis (PCoA) revealed admixture of the lablab germplasm. STRUCTURE, NJ and PCoA analyses also revealed that the subspecies purpureus and bengalensis are not genetically differentiated. Although the number of individuals from the west of Thailand was small and all of them were collected from the same province, they possessed comparable gene diversity with those from the other geographic regions. These results demonstrated that there is moderately low genetic diversity of lablab in Thailand and the west of the country possesses high diversity of lablab.

Genome-Wide SNP Identification and Association Mapping for Seed Mineral Concentration in Mung Bean (Vigna radiata L.)

Mung bean (Vigna radiata L.) quality is dependent on seed chemical composition, which in turn determines the benefits of its consumption for human health and nutrition. While mung bean is rich in a range of nutritional components, such as protein, carbohydrates and vitamins, it remains less well studied than other legume crops in terms of micronutrients. In addition, mung bean genomics and genetic resources are relatively sparse. The objectives of this research were three-fold, namely: to develop a genome-wide marker system for mung bean based on genotyping by sequencing (GBS), to evaluate diversity of mung beans available to breeders in the United States and finally, to perform a genome-wide association study (GWAS) for nutrient concentrations based on a seven mineral analysis using inductively coupled plasma (ICP) spectroscopy. All parts of our research were performed with 95 cultivated mung bean genotypes chosen from the USDA core collection representing accessions from 13 countries. Overall, we identified a total of 6,486 high quality single nucleotide polymorphisms (SNPs) from the GBS dataset and found 43 marker × trait associations (MTAs) with calcium, iron, potassium, manganese, phosphorous, sulfur or zinc concentrations in mung bean grain produced in either of two consecutive years' field experiments. The MTAs were scattered across 35 genomic regions explaining on average 22% of the variation for each seed nutrient in each year. Most of the gene regions provided valuable candidate loci to use in future breeding of new varieties of mung bean and further the understanding of genetic control of nutritional properties in the crop. Other SNPs identified in this study will serve as important resources to enable marker-assisted selection (MAS) for nutritional improvement in mung bean and to analyze cultivars of mung bean.

J, A. S, Bhat KV, Gaikwad AB. Genome wide identification and characterization of microsatellite markers in black pepper (Piper nigrum): A valuable resource for boosting genomics applications

Black pepper is one of the most valued and widely used spices in the world and dominates multi-billion dollar global spices trade. India is amongst the major producers, consumers and exporters of black pepper. In spite of its commercial and cultural importance, black pepper has received meagre attention in terms of generation of genomic resources. Availability of markers distributed throughout the genome would facilitate and accelerate genetic studies, QTL identification, genetic enhancement and crop improvement in black pepper. In this perspective, the sequence information from the recently sequenced black pepper (Piper nigrum) genome has been used for identification and characterisation of Simple Sequence Repeats (SSRs). Total 69,126 SSRs were identified from assembled genomic sequence of P. nigrum. The SSR frequency was 158 per MB making it, one SSR for every 6.3 kb in the assembled genome. Among the different types of microsatellite repeat motifs, dinucleotides were the most abundant (48.6%), followed by trinucleotide (23.7%) and compound repeats (20.62%). A set of 85 SSRs were used for validation, of which 74 produced amplification products of expected size. Genetic diversity of 30 black pepper accessions using 50 SSRs revealed four distinct clusters. Further, the cross species transferability of the SSRs was checked in nine other Piper species. Out of 50 SSRs used, 19 and 31 SSRs were amplified in nine and seven species, respectively. Thus the identified SSRs may have application in other species of the genus Piper where genome sequence is not available yet. Present study reports the first NGS based genomic SSRs in black pepper and thus constitute a valuable resource for a whole fleet of applications in genetics and plant breeding studies such as genetic map construction, QTL identification, map-based gene cloning, marker-assisted selection and evolutionary studies in Piper nigrum and related species.

Interspecific delimitation and relationships among four Ostrya species based on plastomes

BACKGROUND

The genus Ostrya (Betulaceae) contains eight species and four of them are distributed in China. However, studies based on limited informative sites of several chloroplast markers failed to resolve interspecific delimitation and relationships among the four Chinese species. In this study, we aimed to use the whole chloroplast genomes to address these two issues.

RESULTS

We assembled and annotated 33 complete chloroplast genomes (plastomes) of the four Chinese species, representing 17 populations across most of their geographical distributions. Each species contained samples of several individuals that cover most of geographic distributions of the species. All plastomes are highly conserved in genome structure and gene order, with a total length of 158-159 kb and 122 genes. Phylogenetic analyses of whole plastomes, non-coding regions and protein-coding genes produced almost the same topological relationships. In contrast to the well-delimitated species boundary inferred from the nuclear ITS sequence variations, three of the four species are non-monophyletic in the plastome trees, which is consistent with previous studies based on a few chloroplast markers.

CONCLUSIONS

The high incongruence between the ITS and plastome trees may suggest the widespread occurrences of hybrid introgression and incomplete lineage sorting during the divergence of these species. In addition, the plastomes with more informative sites compared with a few chloroplast markers still failed to resolve the phylogenetic relationships of the four species, and further studies involving population genomic data may be needed to better understand their evolutionary histories.

Construction of genetic linkage map and genome dissection of domestication-related traits of moth bean (Vigna aconitifolia), a legume crop of arid areas

The moth bean (Vigna aconitifolia), possibly the most primitive crop of the genus Vigna, is a highly drought- and heat-resistant legume grown in arid areas. Moth bean domestication involved phenotypic changes, including reduction of seed dormancy and pod shattering, increased organ size, and earlier flowering and maturity. However, the genetics of the domestication process in moth bean is not known. In this study, we constructed a genetic linkage map for moth bean and used the map to identify quantitative trait loci (QTL) for domestication-related traits of an F₂ population of 188 individuals produced from a cross of wild moth bean (TN67) and cultivated moth bean (ICPMO056). The genetic linkage map comprised 11 linkage groups (LG) of 172 simple sequence repeat markers and spanned a total length of 1016.8 centiMorgan (cM), with an average marker distance of 7.34 cM. A comparative genome analysis showed high genome synteny between moth bean and mungbean (Vigna radiata), adzuki bean (Vigna angularis), rice bean (Vigna umbellata), and yardlong bean (Vigna unguiculata). In total, 50 QTLs and 3 genes associated with 20 domestication-related traits were identified. Most of the QTLs belonged to five LGs (1, 2, 4, 7, and 10). Key traits related to domestication such as seed dormancy and pod shattering were controlled by large-effect QTLs (PVE > 20%) with one or two minor QTLs, whereas all other traits were controlled by one-seven minor QTLs, apart from seed weight, which was controlled by one major and seven minor QTLs. These results suggest that a small number of mutations with large phenotypic effects have contributed to the domestication of the moth bean. Comparative analysis of QTLs with related Vigna crops revealed that there are several domestication-related large-effect QTLs that had not been used in moth bean domestication. This study provides a basic genetic map and identified genome regions associated with domestication-related traits, which will be useful for the genetic improvement of the moth bean and related Vigna species.

VigSatDB: genome-wide microsatellite DNA marker database of three species of Vigna for germplasm characterization and improvement

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.

QTL analysis of domestication syndrome in zombi pea (Vigna vexillata), an underutilized legume crop

Zombi pea (Vigna vexillata (L.) A. Rich) is an underutilized crop belonging to the genus Vigna. Two domesticated forms of zombi pea are cultivated as crop plants; seed and tuber forms. The cultivated seed form is present in Africa, while the cultivated tuber form is present in a very limited part of Asia. Genetics of domestication have been investigated in most of cultivated Vigna crops by means of quantitative trait locus (QTL) mapping. In this study, we investigated genetics of domestication in zombi pea by QTL analysis using an F2 population of 139 plants derived from a cross between cultivated tuber form of V. vexillata (JP235863) and wild V. vexillata (AusTRCF66514). A linkage map with 11 linkage groups (LGs) was constructed from this F2 population using 145 SSR, 117 RAD-seq and 2 morphological markers. Many highly segregation distorted markers were found on LGs 5, 6, 7, 8, 10 and 11. Most of the distorted markers were clustered together and all the markers on LG8 were highly distorted markers. Comparing this V. vexillata linkage map with linkage maps of other four Vigna species demonstrated several genome rearrangements in V. vexillata. QTL analysis for 22 domestication-related traits was investigated by inclusive composite interval mapping in which 37 QTLs were identified for 18 traits; no QTL was detected for 4 traits. Number of QTLs detected in each trait ranged from 1 to 5 with an average of only 2.3. Five QTLs for tuber width and three QTLs for tuber weight. Interestingly, 2 QTLs each for tuber width and tuber weight detected on LG2 and LG4 were located at similar position and wild allele increased tuber width and weight. This indicated wild germplasm having small tuber have potential to increase yield of large tuber cultivated type. Large-effect QTLs (PVE > 20%) were on LG4 (pod length), LG5 (leaf size and seed thickness), and LG7 (for seed-related traits). Comparison of domestication-related QTLs of the zombi pea with those of cowpea (Vigna unguiculata), azuki bean (Vigna angularis), mungbean (Vigna radiata) and rice bean (Vigna umbellata) revealed that there was conservation of some QTLs for seed size, pod size and leaf size between zombi pea and cowpea and that QTLs associated with seed size (weight, length, width and thickness) in each species were clustered on same linkage.

The Linum usitatissimum L. plastome reveals atypical structural evolution, new editing sites, and the phylogenetic position of Linaceae within Malpighiales

The plastome of Linum usitatissimum was completely sequenced allowing analyses of evolution of genome structure, RNA editing sites, molecular markers, and indicating the position of Linaceae within Malpighiales. Flax (Linum usitatissimum L.) is an economically important crop used as food, feed, and industrial feedstock. It belongs to the Linaceae family, which is noted by high morphological and ecological diversity. Here, we reported the complete sequence of flax plastome, the first species within Linaceae family to have the plastome sequenced, assembled and characterized in detail. The plastome of flax is a circular DNA molecule of 156,721 bp with a typical quadripartite structure including two IRs of 31,990 bp separating the LSC of 81,767 bp and the SSC of 10,974 bp. It shows two expansion events from IRB to LSC and from IRB to SSC, and a contraction event in the IRA-LSC junction, which changed significantly the size and the gene content of LSC, SSC and IRs. We identified 109 unique genes and 2 pseudogenes (rpl23 and ndhF). The plastome lost the conserved introns of clpP gene and the complete sequence of rps16 gene. The clpP, ycf1, and ycf2 genes show high nucleotide and aminoacid divergence, but they still possibly retain the functionality. Moreover, we also identified 176 SSRs, 20 tandem repeats, and 39 dispersed repeats. We predicted in 18 genes a total of 53 RNA editing sites of which 32 were not found before in other species. The phylogenetic inference based on 63 plastid protein-coding genes of 38 taxa supports three major clades within Malpighiales order. One of these clades has flax (Linaceae) sister to Chrysobalanaceae family, differing from earlier studies that included Linaceae into the euphorbioid clade.

Quantitative trait locus mapping under irrigated and drought treatments based on a novel genetic linkage map in mungbean (Vigna radiata L.)

A novel genetic linkage map was constructed using SSR markers and stable QTLs were identified for six drought tolerance related-traits using single-environment analysis under irrigation and drought treatments. Mungbean (Vigna radiata L.) is one of the most important leguminous food crops. However, mungbean production is seriously constrained by drought. Isolation of drought-responsive genetic elements and marker-assisted selection breeding will benefit from the detection of quantitative trait locus (QTLs) for traits related to drought tolerance. In this study, we developed a full-coverage genetic linkage map based on simple sequence repeat (SSR) markers using a recombinant inbred line (RIL) population derived from an intra-specific cross between two drought-resistant varieties. This novel map was anchored with 313 markers. The total map length was 1010.18 cM across 11 linkage groups, covering the entire genome of mungbean with a saturation of one marker every 3.23 cM. We subsequently detected 58 QTLs for plant height (PH), maximum leaf area (MLA), biomass (BM), relative water content, days to first flowering, and seed yield (Yield) and 5 for the drought tolerance index of 3 traits in irrigated and drought environments at 2 locations. Thirty-eight of these QTLs were consistently detected two or more times at similar linkage positions. Notably, qPH5A and qMLA2A were consistently identified in marker intervals from GMES5773 to MUS128 in LG05 and from Mchr11-34 to the HAAS_VR_1812 region in LG02 in four environments, contributing 6.40-20.06% and 6.97-7.94% of the observed phenotypic variation, respectively. None of these QTLs shared loci with previously identified drought-related loci from mungbean. The results of these analyses might facilitate the isolation of drought-related genes and help to clarify the mechanism of drought tolerance in mungbean.

Multiple organ gigantism caused by mutation in VmPPD gene in blackgram (Vigna mungo)

Seed size is one of the most important traits in leguminous crops. We obtained a recessive mutant of blackgram that had greatly enlarged leaves, stems and seeds. The mutant produced 100% bigger leaves, 50% more biomass and 70% larger seeds though it produced 40% less number of seeds. We designated the mutant as multiple-organ-gigantism (mog) and found the mog phenotype was due to increase in cell numbers but not in cell size. We also found the mog mutant showed a rippled leaf (rl) phenotype, which was probably caused by a pleiotropic effect of the mutation. We performed a map-based cloning and successfully identified an 8 bp deletion in the coding sequence of VmPPD gene, an orthologue of Arabidopsis PEAPOD (PPD) that regulates arrest of cell divisions in meristematic cells. We found no other mutations in the neighboring genes between the mutant and the wild type. We also knocked down GmPPD genes and reproduced both the mog and rl phenotypes in soybean. Controlling PPD genes to produce the mog phenotype is highly valuable for breeding since larger seed size could directly increase the commercial values of grain legumes.

Genetic diversity and structure of the zombi pea (Vigna vexillata (L.) A. Rich) gene pool based on SSR marker analysis

Zombi pea (Vigna vexillata (L.) A. Rich) is an underutilized legume species and a useful gene source for resistance to biotic and abiotic stresses, although there is little understanding on its genetic diversity and structure. In this study, 422 (408 wild and 14 cultivated) accessions of zombi pea from diverse origins (201 from Africa, 126 from America, 85 from Australia, 5 from Asia and 5 from unknown origin) were analyzed with 20 simple sequence repeat (SSR) markers to determine its genetic diversity and genetic structure. The SSR markers detected 273 alleles in total with a mean of 13.6 alleles per locus. Polymorphism information content values of the markers varied from 0.58 to 0.90 with an average of 0.76. Overall gene diversity was 0.715. Gene diversity and average allelic richness was highest in Africa (0.749 and 8.08, respectively) and lowest in America (0.435 and 4.10, respectively). Nei's genetic distance analysis revealed that the highest distance was between wild Australia and cultivated Africa (0.559), followed by wild West Africa and wild Australia (0.415). STRUCTURE, neighbor-joining (NJ), and principal coordinate analyses consistently showed that these zombi pea accessions were clustered into three major groups, viz. America, Africa and Asia, and Australia. NJ tree also suggested that American and Australian accessions are originated from East African zombi peas, and that the cultivated accessions from Africa and Asia were genetically distinct, while those from America were clustered with some cultivated accessions from Africa. These results suggest that Africa is the center of origin and diversity of zombi pea, and that domestication of this pea took place more than once in different regions.

Genome wide characterization of simple sequence repeats in watermelon genome and their application in comparative mapping and genetic diversity analysis

BACKGROUND

Microsatellite markers are one of the most informative and versatile DNA-based markers used in plant genetic research, but their development has traditionally been difficult and costly. The whole genome sequencing with next-generation sequencing (NGS) technologies provides large amounts of sequence data to develop numerous microsatellite markers at whole genome scale. SSR markers have great advantage in cross-species comparisons and allow investigation of karyotype and genome evolution through highly efficient computation approaches such as in silico PCR. Here we described genome wide development and characterization of SSR markers in the watermelon (Citrullus lanatus) genome, which were then use in comparative analysis with two other important crop species in the Cucurbitaceae family: cucumber (Cucumis sativus L.) and melon (Cucumis melo L.). We further applied these markers in evaluating the genetic diversity and population structure in watermelon germplasm collections.

RESULTS

A total of 39,523 microsatellite loci were identified from the watermelon draft genome with an overall density of 111 SSRs/Mbp, and 32,869 SSR primers were designed with suitable flanking sequences. The dinucleotide SSRs were the most common type representing 34.09 % of the total SSR loci and the AT-rich motifs were the most abundant in all nucleotide repeat types. In silico PCR analysis identified 832 and 925 SSR markers with each having a single amplicon in the cucumber and melon draft genome, respectively. Comparative analysis with these cross-species SSR markers revealed complicated mosaic patterns of syntenic blocks among the genomes of three species. In addition, genetic diversity analysis of 134 watermelon accessions with 32 highly informative SSR loci placed these lines into two groups with all accessions of C.lanatus var. citorides and three accessions of C. colocynthis clustered in one group and all accessions of C. lanatus var. lanatus and the remaining accessions of C. colocynthis clustered in another group. Furthermore, structure analysis was consistent with the dendrogram indicating the 134 watermelon accessions were classified into two populations.

CONCLUSION

The large number of genome wide SSR markers developed herein from the watermelon genome provides a valuable resource for genetic map construction, QTL exploration, map-based gene cloning and marker-assisted selection in watermelon which has a very narrow genetic base and extremely low polymorphism among cultivated lines. Furthermore, the cross-species transferable SSR markers identified herein should also have practical uses in many applications in species of Cucurbitaceae family whose whole genome sequences are not yet available.

Identification and Characterization of Microsatellite Loci in Maqui (Aristotelia chilensis [Molina] Stunz) Using Next-Generation Sequencing (NGS)

Maqui (Aristotelia chilensis [Molina] Stunz) is a small dioecious tree native to South America with edible fruit characterized by very high antioxidant capacity and anthocyanin content. To preserve maqui as a genetic resource it is essential to study its genetic diversity. However, the complete genome is unknown and only a few gene sequences are available in databases. Simple sequence repeats (SSR) markers, which are neutral, co-dominant, reproducible and highly variable, are desirable to support genetic studies in maqui populations. By means of identification and characterization of microsatellite loci from a maqui genotype, using 454 sequencing technology, we develop a set of SSR for this species. Obtaining a total of 165,043 shotgun genome sequences, with an average read length of 387 bases, we covered 64 Mb of the maqui genome. Reads were assembled into 4,832 contigs, while 98,546 reads remained as singletons, generating a total of 103,378 consensus genomic sequences. A total of 24,494 SSR maqui markers were identified. Of them, 15,950 SSR maqui markers were classified as perfects. The most common SSR motifs were dinucleotide (31%), followed by tetranucleotide (26%) and trinucleotide motifs (24%). The motif AG/CT (28.4%) was the most abundant, while the motif AC (89 bp) was the largest. Eleven polymorphic SSRs were selected and used to analyze a population of 40 maqui genotypes. Polymorphism information content (PIC) ranged from 0.117 to 0.82, with an average of 0.58. Non-significant groups were observed in the maqui population, showing a panmictic genetic structure. In addition, we also predicted 11150 putative genes and 3 microRNAs (miRNAs) in maqui sequences. This results, including partial sequences of genes, some miRNAs and SSR markers from high throughput next generation sequencing (NGS) of maqui genomic DNA, constitute the first platform to undertake genetic and molecular studies of this important species.

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

Mungbean (Vigna radiata L.) is a cleistogamous plant in which flowers are pollinated before they open, which prevents yield improvements through heterosis. We previously generated a chasmogamous mutant (CM) mungbean in which open flowers are pollinated. In this study, we developed insertion/deletion (indel) markers based on the transcriptome differences between CM and Sulu-1 (i.e., normal flowering) plants. An F2 population derived from a cross between CM and Sulu-1 was used for gene mapping. Segregation analyses revealed that a single recessive gene regulates the production of chasmogamous flowers. Using newly developed indel and simple sequence repeat markers, the cha gene responsible for the chasmogamous flower trait was mapped to a 277.1-kb segment on chromosome 6. Twelve candidate genes were detected in this segment, including Vradi06g12650, which encodes a YUCCA family protein associated with floral development. A single base pair deletion producing a frame-shift mutation and a premature stop codon in Vradi06g12650 was detected only in CM plants. This suggested that Vradi06g12650 is a cha candidate gene. Our results provide important information for the molecular breeding of chasmogamous mungbean lines, which may serve as new genetic resources for hybrid cultivar development.

High-Throughput Development of SSR Markers from Pea (Pisum sativum L.) Based on Next Generation Sequencing of a Purified Chinese Commercial Variety

Pea (Pisum sativum L.) is an important food legume globally, and is the plant species that J.G. Mendel used to lay the foundation of modern genetics. However, genomics resources of pea are limited comparing to other crop species. Application of marker assisted selection (MAS) in pea breeding has lagged behind many other crops. Development of a large number of novel and reliable SSR (simple sequence repeat) or microsatellite markers will help both basic and applied genomics research of this crop. The Illumina HiSeq 2500 System was used to uncover 8,899 putative SSR containing sequences, and 3,275 non-redundant primers were designed to amplify these SSRs. Among the 1,644 SSRs that were randomly selected for primer validation, 841 yielded reliable amplifications of detectable polymorphisms among 24 genotypes of cultivated pea (Pisum sativum L.) and wild relatives (P. fulvum Sm.) originated from diverse geographical locations. The dataset indicated that the allele number per locus ranged from 2 to 10, and that the polymorphism information content (PIC) ranged from 0.08 to 0.82 with an average of 0.38. These 1,644 novel SSR markers were also tested for polymorphism between genotypes G0003973 and G0005527. Finally, 33 polymorphic SSR markers were anchored on the genetic linkage map of G0003973 × G0005527 F2 population.

Development of Genomic Microsatellite Markers in Carthamus tinctorius L. (Safflower) Using Next Generation Sequencing and Assessment of Their Cross-Species Transferability and Utility for Diversity Analysis

BACKGROUND

Safflower (Carthamus tinctorius L.), an Asteraceae member, yields high quality edible oil rich in unsaturated fatty acids and is resilient to dry conditions. The crop holds tremendous potential for improvement through concerted molecular breeding programs due to the availability of significant genetic and phenotypic diversity. Genomic resources that could facilitate such breeding programs remain largely underdeveloped in the crop. The present study was initiated to develop a large set of novel microsatellite markers for safflower using next generation sequencing.

PRINCIPAL FINDINGS

Low throughput genome sequencing of safflower was performed using Illumina paired end technology providing ~3.5X coverage of the genome. Analysis of sequencing data allowed identification of 23,067 regions harboring perfect microsatellite loci. The safflower genome was found to be rich in dinucleotide repeats followed by tri-, tetra-, penta- and hexa-nucleotides. Primer pairs were designed for 5,716 novel microsatellite sequences with repeat length ≥ 20 bases and optimal flanking regions. A subset of 325 microsatellite loci was tested for amplification, of which 294 loci produced robust amplification. The validated primers were used for assessment of 23 safflower accessions belonging to diverse agro-climatic zones of the world leading to identification of 93 polymorphic primers (31.6%). The numbers of observed alleles at each locus ranged from two to four and mean polymorphism information content was found to be 0.3075. The polymorphic primers were tested for cross-species transferability on nine wild relatives of cultivated safflower. All primers except one showed amplification in at least two wild species while 25 primers amplified across all the nine species. The UPGMA dendrogram clustered C. tinctorius accessions and wild species separately into two major groups. The proposed progenitor species of safflower, C. oxyacantha and C. palaestinus were genetically closer to cultivated safflower and formed a distinct cluster. The cluster analysis also distinguished diploid and tetraploid wild species of safflower.

CONCLUSION

Next generation sequencing of safflower genome generated a large set of microsatellite markers. The novel markers developed in this study will add to the existing repertoire of markers and can be used for diversity analysis, synteny studies, construction of linkage maps and marker-assisted selection.

De novo Assembly, Characterization of Immature Seed Transcriptome and Development of Genic-SSR Markers in Black Gram [Vigna mungo (L.) Hepper]

Black gram [V. mungo (L.) Hepper] is an important legume crop extensively grown in south and south-east Asia, where it is a major source of dietary protein for its predominantly vegetarian population. However, lack of genomic information and markers has become a limitation for genetic improvement of this crop. Here, we report the transcriptome sequencing of the immature seeds of black gram cv. TU94-2, by Illumina paired end sequencing technology to generate transcriptome sequences for gene discovery and genic-SSR marker development. A total of 17.2 million paired-end reads were generated and 48,291 transcript contigs (TCS) were assembled with an average length of 443 bp. Based on sequence similarity search, 33,766 TCS showed significant similarity to known proteins. Among these, only 29,564 TCS were annotated with gene ontology (GO) functional categories. A total number of 138 unique KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways were identified, of which majority of TCS are grouped into purine metabolism (678) followed by pyrimidine metabolism (263). A total of 48,291 TCS were searched for SSRs and 1,840 SSRs were identified in 1,572 TCS with an average frequency of one SSR per 11.9 kb. The tri-nucleotide repeats were most abundant (35%) followed by di-nucleotide repeats (32%). PCR primer pairs were successfully designed for 933 SSR loci. Sequences analyses indicate that about 64.4% and 35.6% of the SSR motifs were present in the coding sequences (CDS) and untranslated regions (UTRs) respectively. Tri-nucleotide repeats (57.3%) were preferentially present in the CDS. The rate of successful amplification and polymorphism were investigated using selected primers among 18 black gram accessions. Genic-SSR markers developed from the Illumina paired end sequencing of black gram immature seed transcriptome will provide a valuable resource for genetic diversity, evolution, linkage mapping, comparative genomics and marker-assisted selection in black gram.

The AVRDC - The World Vegetable Center mungbean (Vigna radiata) core and mini core collections

Background

Large ex situ germplasm collections generally harbor a wide range of crop diversity. AVRDC – The World Vegetable Center is holding in trust the world’s second largest mungbean (Vigna radiata) germplasm collection with more than 6,700 accessions. Screening large collections for traits of interest is laborious and expensive. To enhance the access of breeders to the diversity of the crop, mungbean core and mini core collections have been established.

Results

The core collection of 1,481 entries has been built by random selection of 20% of the accessions after geographical stratification and subsequent cluster analysis of eight phenotypic descriptors in the whole collection. Summary statistics, especially the low differences of means, equal variance of the traits in both the whole and core collection and the visual inspection of quantile-quantile plots comparing the variation of phenotypic traits present in both collections indicated that the core collection well represented the pattern of diversity of the whole collection. The core collection was genotyped with 20 simple sequence repeat markers and a mini core set of 289 accessions was selected, which depicted the allele and genotype diversity of the core collection.

Conclusions

The mungbean core and mini core collections plus their phenotypic and genotypic data are available for distribution to breeders. It is expected that these collections will enhance the access to biodiverse mungbean germplasm for breeding.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1556-7) contains supplementary material, which is available to authorized users.

Assessment of genetic diversity and population structure of mung bean (Vigna radiata) germplasm using EST-based and genomic SSR markers

Mung bean is an important legume crop in tropical and subtropical countries of Asia and has high nutritional and economic value. However the genetic diversity of mung bean is poorly characterized. In this study, our goal was to develop and use microsatellite simple sequence repeat (SSR) markers for germplasm evaluation. In total, 500 novel expression sequence tag EST-based SSRs (eSSRs) and genomic SSRs (gSSRs) were developed from mung bean transcriptome and genome sequences. Of these, only 58 were useful for diversity evaluation in a panel of 157 cultivated and wild mung bean accessions from different collection sites in East Asia. A total of 2.66 alleles were detected on average per locus which shows that polymorphism is generally low for the species. The average polymorphic information content (PIC) of gSSRs was higher than eSSRs and most of the polymorphic gSSRs were composed of di- and tri-nucleotide repeats (52.4% and 38.1% of all loci, respectively). The genotypes were differentiated into nine subgroups by cluster analysis, and the wild mung bean accessions separated well from the cultivated accessions. Analysis of molecular variance indicated that 22% of variance was observed among populations and 78% was due to differences within populations. Clustering, population structure analyses showed that non-Chinese cultivated and wild mung bean accessions were separated from Chinese accessions, but no geographical distinctions existed between genotypes collected in China. Interestingly, the average PIC value of cultivated mung bean (0.36) was higher than that of wild mung bean (0.25) showing that further collecting and wide crosses are necessary for mung bean improvement.

Genetic diversity of the black gram [Vigna mungo (L.) Hepper] gene pool as revealed by SSR markers

In this study, 520 cultivated and 14 wild accessions of black gram (Vigna mungo (L.) Hepper) were assessed for diversity using 22 SSR markers. Totally, 199 alleles were detected with a mean of 9.05 alleles per locus. Wild black gram showed higher gene diversity than cultivated black gram. Gene diversity of cultivated accessions among regions was comparable, while allelic richness of South Asia was higher than that of other regions. 78.67% of the wild gene diversity presented in cultivated accessions, indicating that the domestication bottleneck effect in black gram is relatively low. Genetic distance analysis revealed that cultivated black gram was more closely related to wild black gram from South Asia than that from Southeast Asia. STRUCTURE, principal coordinate and neighbor-joining analyses consistently revealed that 534 black gram accessions were grouped into three major subpopulations. The analyses also revealed that cultivated black gram from South Asia was genetically distinct from that from West Asia. Comparison by SSR analysis with other closely related Vigna species, including mungbean, azuki bean, and rice bean, revealed that level of gene diversity of black gram is comparable to that of mungbean and rice bean but lower than that of azuki bean.

Transcriptome sequencing and marker development for four underutilized legumes

•

PREMISE OF THE STUDY

Combating threats to food and nutrition security in the context of climate change and global population increase is one of the highest priorities of major international organizations. Hundreds of species are grown on a small scale in some of the most drought/flood-prone regions of the world and as such may harbor some of the most environmentally tolerant crops (and alleles). •

METHODS AND RESULTS

In this study, transcriptomes were sequenced, assembled, and annotated for four underutilized legume crops. Microsatellite markers were identified in each species, as well as a conserved orthologous set of markers for cross-family phylogenetics and comparative mapping, which were ground-truthed on a panel of diverse legume germplasm. •

CONCLUSIONS

An understanding of these underutilized legumes will inform crop selection and breeding by allowing the investigation of genetic variation and the genetic basis of adaptive traits to be established.

Transcriptome sequencing and marker development for four underutilized legumes

•

PREMISE OF THE STUDY

Combating threats to food and nutrition security in the context of climate change and global population increase is one of the highest priorities of major international organizations. Hundreds of species are grown on a small scale in some of the most drought/flood-prone regions of the world and as such may harbor some of the most environmentally tolerant crops (and alleles). •

METHODS AND RESULTS

In this study, transcriptomes were sequenced, assembled, and annotated for four underutilized legume crops. Microsatellite markers were identified in each species, as well as a conserved orthologous set of markers for cross-family phylogenetics and comparative mapping, which were ground-truthed on a panel of diverse legume germplasm. •

CONCLUSIONS

An understanding of these underutilized legumes will inform crop selection and breeding by allowing the investigation of genetic variation and the genetic basis of adaptive traits to be established.

Genomic resources in mungbean for future breeding programs

Among the legume family, mungbean (Vigna radiata) has become one of the important crops in Asia, showing a steady increase in global production. It provides a good source of protein and contains most notably folate and iron. Beyond the nutritional value of mungbean, certain features make it a well-suited model organism among legume plants because of its small genome size, short life-cycle, self-pollinating, and close genetic relationship to other legumes. In the past, there have been several efforts to develop molecular markers and linkage maps associated with agronomic traits for the genetic improvement of mungbean and, ultimately, breeding for cultivar development to increase the average yields of mungbean. The recent release of a reference genome of the cultivated mungbean (V. radiata var. radiata VC1973A) and an additional de novo sequencing of a wild relative mungbean (V. radiata var. sublobata) has provided a framework for mungbean genetic and genome research, that can further be used for genome-wide association and functional studies to identify genes related to specific agronomic traits. Moreover, the diverse gene pool of wild mungbean comprises valuable genetic resources of beneficial genes that may be helpful in widening the genetic diversity of cultivated mungbean. This review paper covers the research progress on molecular and genomics approaches and the current status of breeding programs that have developed to move toward the ultimate goal of mungbean improvement.

Genomics and molecular breeding in lesser explored pulse crops: current trends and future opportunities

Pulses are multipurpose crops for providing income, employment and food security in the underprivileged regions, notably the FAO-defined low-income food-deficit countries. Owing to their intrinsic ability to endure environmental adversities and the least input/management requirements, these crops remain central to subsistence farming. Given their pivotal role in rain-fed agriculture, substantial research has been invested to boost the productivity of these pulse crops. To this end, genomic tools and technologies have appeared as the compelling supplement to the conventional breeding. However, the progress in minor pulse crops including dry beans (Vigna spp.), lupins, lablab, lathyrus and vetches has remained unsatisfactory, hence these crops are often labeled as low profile or lesser researched. Nevertheless, recent scientific and technological breakthroughs particularly the next generation sequencing (NGS) are radically transforming the scenario of genomics and molecular breeding in these minor crops. NGS techniques have allowed de novo assembly of whole genomes in these orphan crops. Moreover, the availability of a reference genome sequence would promote re-sequencing of diverse genotypes to unlock allelic diversity at a genome-wide scale. In parallel, NGS has offered high-resolution genetic maps or more precisely, a robust genetic framework to implement whole-genome strategies for crop improvement. As has already been demonstrated in lupin, sequencing-based genotyping of the representative sample provided access to a number of functionally-relevant markers that could be deployed straight away in crop breeding programs. This article attempts to outline the recent progress made in genomics of these lesser explored pulse crops, and examines the prospects of genomics assisted integrated breeding to enhance and stabilize crop yields.

Detection of genome donor species of neglected tetraploid crop Vigna reflexo-pilosa (créole bean), and genetic structure of diploid species based on newly developed EST-SSR markers from azuki bean (Vigna angularis)

Vigna reflexo-pilosa, which includes a neglected crop, is the only one tetraploid species in genus Vigna. The ancestral species that make up this allotetraploid species have not conclusively been identified, although previous studies suggested that a donor genome of V. reflexo-pilosa is V. trinervia. In this study, 1,429 azuki bean EST-SSR markers were developed of which 38 EST-SSR primer pairs that amplified one product in diploid species and two discrete products in tetraploid species were selected to analyze 268 accessions from eight taxa of seven Asian Vigna species including V. reflexo-pilosa var. glabra, V. reflexo-pilosa var. reflexo-pilosa, V. exilis, V. hirtella, V. minima, V. radiata var. sublobata, V. tenuicaulis and V. trinervia to identify genome donor of V. reflexo-pilosa. Since both diploid and tetraploid species were analyzed and each SSR primer pair detected two loci in the tetraploid species, we separated genomes of the tetraploid species into two different diploid types, viz. A and B. In total, 445 alleles were detected by 38 EST-SSR markers. The highest gene diversity was observed in V. hirtella. By assigning the discrete PCR products of V. reflexo-pilosa into two distinguished genomes, we were able to identify the two genome donor parents of créole bean. Phylogenetic and principal coordinate analyses suggested that V. hirtella is a species complex and may be composed of at least three distinct taxa. Both analyses also clearly demonstrated that V. trinervia and one taxon of V. hirtella are the genome donors of V. reflexo-pilosa. Gene diversity indicates that the evolution rate of EST-SSRs on genome B of créole bean might be faster than that on genome A. Species relationship among the Vigna species in relation to genetic data, morphology and geographical distribution are presented.

Expanding the repertoire of microsatellite markers for polymorphism studies in Indian accessions of mung bean (Vigna radiata L. Wilczek)

Limited availability of validated, polymorphic microsatellite markers in mung bean (Vigna radiata), an important food legume of India, has been a major hurdle towards its improvement and higher yield. The present study was undertaken in order to develop a new set of microsatellite markers and utilize them for the analysis of genetic diversity within mung bean accessions from India. A GA/CT enriched library was constructed from V. radiata which resulted in 1,250 putative recombinant clones of which 850 were sequenced. SSR motifs were identified and their flanking sequences were utilized to design 328 SSR primer pairs. Of these, 48 SSR markers were employed for assessing genetic diversity among 76 mung bean accessions from various geographical locations in India. Two hundred and thirty four alleles with an average of 4.85 alleles per locus were detected at 48 loci. The polymorphic information content (PIC) per locus varied from 0.1 to 0.88 (average: 0.49 per locus). The observed and expected heterozygosities ranged from 0.40 to 0.95 and 0.40 to 0.81 respectively. Based on Jaccard's similarity matrix, a dendrogram was constructed using the unweighted pair-group method with arithmetic averages (UPGMA) analysis which revealed that one accession from Bundi, Rajasthan was clustered out separately while remaining accessions were grouped into two major clusters. The markers generated in this study will help in expanding the repertoire of the available SSR markers thereby facilitating analysis of genetic diversity, molecular mapping and ultimately broadening the scope for genetic improvement of this legume.

Large-scale microsatellite development in grasspea (Lathyrus sativus L.), an orphan legume of the arid areas

BACKGROUND

Grasspea (Lathyrus sativus L., 2n = 14), a member of the family Leguminosae, holds great agronomic potential as grain and forage legume crop in the arid areas for its superb resilience to abiotic stresses such as drought, flood and salinity. The crop could not make much progress through conventional breeding in the past, and there are hardly any detailed molecular biology studies due to paucity of reliable molecular markers representative of the entire genome.

RESULTS

Using the 454 FLX Titanium pyrosequencing technique, 651,827 simple sequence repeat (SSR) loci were identified and 50,144 nonredundant primer pairs were successfully designed, of which 288 were randomly selected for validation among 23 L. sativus and one L. cicera accessions of diverse provenance. 74 were polymorphic, 70 monomorphic, and 144 with no PCR product. The number of observed alleles ranged from two to five, the observed heterozygosity from 0 to 0.9545, and Shannon's information index ranged from 0.1013 to 1.0980, respectively. The dendrogram constructed by using unweighted pair group method with arithmetic mean (UPGMA) based on Nei's genetic distance, showed obvious distinctions and understandable relationships among the 24 accessions.

CONCLUSIONS

The large number of SSR primer pairs developed in this study would make a significant contribution to genomics enabled improvement of grasspea.

QTL mapping for salt tolerance and domestication-related traits in Vigna marina subsp. oblonga, a halophytic species

QTL mapping in F2 population [V. luteola × V. marina subsp. oblonga] revealed that the salt tolerance in V. marina subsp. oblonga is controlled by a single major QTL. The habitats of beach cowpea (Vigna marina) are sandy beaches in tropical and subtropical regions. As a species that grows closest to the sea, it has potential to be a gene source for breeding salt-tolerant crops. We reported here for the first time, quantitative trait loci (QTLs) mapping for salt tolerance in V. marina. A genetic linkage map was constructed from an F2 population of 120 plants derived from an interspecific cross between V. luteola and V. marina subsp. oblonga. The map comprised 150 SSR markers. The markers were clustered into 11 linkage groups spanning 777.6 cM in length with a mean distance between the adjacent markers of 5.59 cM. The F2:3 population was evaluated for salt tolerance under hydroponic conditions at the seedling and developmental stages. Segregation analysis indicated that salt tolerance in V. marina is controlled by a few genes. Multiple interval mapping consistently identified one major QTL which can explain about 50% of phenotypic variance. The flanking markers may facilitate transfer of the salt tolerance allele from V. marina subsp. oblonga into related Vigna crops. The QTL for domestication-related traits from V. marina are also discussed.

Gene discovery and molecular marker development, based on high-throughput transcript sequencing of Paspalum dilatatum Poir

Background

Paspalum dilatatum Poir. (common name dallisgrass) is a native grass species of South America, with special relevance to dairy and red meat production. P. dilatatum exhibits higher forage quality than other C4 forage grasses and is tolerant to frost and water stress. This species is predominantly cultivated in an apomictic monoculture, with an inherent high risk that biotic and abiotic stresses could potentially devastate productivity. Therefore, advanced breeding strategies that characterise and use available genetic diversity, or assess germplasm collections effectively are required to deliver advanced cultivars for production systems. However, there are limited genomic resources available for this forage grass species.

Results

Transcriptome sequencing using second-generation sequencing platforms has been employed using pooled RNA from different tissues (stems, roots, leaves and inflorescences) at the final reproductive stage of P. dilatatum cultivar Primo. A total of 324,695 sequence reads were obtained, corresponding to c. 102 Mbp. The sequences were assembled, generating 20,169 contigs of a combined length of 9,336,138 nucleotides. The contigs were BLAST analysed against the fully sequenced grass species of Oryza sativa subsp. japonica, Brachypodium distachyon, the closely related Sorghum bicolor and foxtail millet (Setaria italica) genomes as well as against the UniRef 90 protein database allowing a comprehensive gene ontology analysis to be performed. The contigs generated from the transcript sequencing were also analysed for the presence of simple sequence repeats (SSRs). A total of 2,339 SSR motifs were identified within 1,989 contigs and corresponding primer pairs were designed. Empirical validation of a cohort of 96 SSRs was performed, with 34% being polymorphic between sexual and apomictic biotypes.

Conclusions

The development of genetic and genomic resources for P. dilatatum will contribute to gene discovery and expression studies. Association of gene function with agronomic traits will significantly enable molecular breeding and advance germplasm enhancement.

DNA fingerprinting in botany: past, present, future

Almost three decades ago Alec Jeffreys published his seminal Nature papers on the use of minisatellite probes for DNA fingerprinting of humans (Jeffreys and colleagues Nature 1985, 314:67-73 and Nature 1985, 316:76-79). The new technology was soon adopted for many other organisms including plants, and when Hilde Nybom, Kurt Weising and Alec Jeffreys first met at the very First International Conference on DNA Fingerprinting in Berne, Switzerland, in 1990, everybody was enthusiastic about the novel method that allowed us for the first time to discriminate between humans, animals, plants and fungi on the individual level using DNA markers. A newsletter coined "Fingerprint News" was launched, T-shirts were sold, and the proceedings of the Berne conference filled a first book on "DNA fingerprinting: approaches and applications". Four more conferences were about to follow, one on each continent, and Alec Jeffreys of course was invited to all of them. Since these early days, methodologies have undergone a rapid evolution and diversification. A multitude of techniques have been developed, optimized, and eventually abandoned when novel and more efficient and/or more reliable methods appeared. Despite some overlap between the lifetimes of the different technologies, three phases can be defined that coincide with major technological advances. Whereas the first phase of DNA fingerprinting ("the past") was dominated by restriction fragment analysis in conjunction with Southern blot hybridization, the advent of the PCR in the late 1980s gave way to the development of PCR-based single- or multi-locus profiling techniques in the second phase. Given that many routine applications of plant DNA fingerprinting still rely on PCR-based markers, we here refer to these methods as "DNA fingerprinting in the present", and include numerous examples in the present review. The beginning of the third phase actually dates back to 2005, when several novel, highly parallel DNA sequencing strategies were developed that increased the throughput over current Sanger sequencing technology 1000-fold and more. High-speed DNA sequencing was soon also exploited for DNA fingerprinting in plants, either in terms of facilitated marker development, or directly in the sense of "genotyping-by-sequencing". Whereas these novel approaches are applied at an ever increasing rate also in non-model species, they are still far from routine, and we therefore treat them here as "DNA fingerprinting in the future".

De novo transcriptome assembly and novel microsatellite marker information in Capsicum annuum varieties Saengryeg 211 and Saengryeg 213

BACKGROUND

Pepper, Capsicum annuum L., Solanaceae, is a major staple economically important vegetable crop worldwide. Limited functional genomics resources and whole genome association studies could be substantially improved through the application of molecular approach for the characterization of gene content and identification of molecular markers. The massive parallel pyrosequencing of two pepper varieties, the highly pungent, Saengryeg 211, and the non-pungent, Saengryeg 213, including de novo transcriptome assembly, functional annotation, and in silico discovery of potential molecular markers is described. We performed 454 GS-FLX Titanium sequencing of polyA-selected and normalized cDNA libraries generated from a single pool of transcripts obtained from mature fruits of two pepper varieties.

RESULTS

A single 454 pyrosequencing run generated 361,671 and 274,269 reads totaling 164.49 and 124.60 Mb of sequence data (average read length of 454 nucleotides), which assembled into 23,821 and 17,813 isotigs and 18,147 and 15,129 singletons for both varieties, respectively. These reads were organized into 20,352 and 15,781 'isogroups' for both varieties. Assembled sequences were functionally annotated based on homology to genes in multiple public databases and assigned with Gene Ontology (GO) terms. Sequence variants analyses identified a total of 3,766 and 2,431 potential (Simple Sequence Repeat) SSR motifs for microsatellite analysis for both varieties, where trinucleotide was the most common repeat unit (84%), followed by di (9.9%), hexa (4.1%) and pentanucleotide repeats (2.1%). GAA repeat (8.6%) was the most frequent repeat motif, followed by TGG (7.2%), TTC (6.5%), and CAG (6.2%).

CONCLUSIONS

High-throughput transcriptome assembly, annotation and large scale of SSR marker discovery has been achieved using next generation sequencing (NGS) of two pepper varieties. These valuable informations for functional genomics resource shall help to further improve the pepper breeding efforts with respect to genetic linkage maps, QTL mapping and marker-assisted trait selection.

Genome-wide SNP discovery in mungbean by Illumina HiSeq

Mungbean [Vigna radiata (L.) Wilczek], a self-pollinated diploid plant with 2n = 22 chromosomes, is an important legume crop with a high-quality amino acid profile. Sequence variation at the whole-genome level was examined by comparing two mungbean cultivars, Sunhwanokdu and Gyeonggijaerae 5, using Illumina HiSeq sequencing data. More than 40 billion bp from both mungbean cultivars were sequenced to a depth of 72×. After de novo assembly of Sunhwanokdu contigs by ABySS 1.3.2 (N50 = 9,958 bp), those longer than 10 kb were aligned with Gyeonggijaerae 5 reads using the Burrows-Wheeler Aligner. SAMTools was used for retrieving single nucleotide polymorphisms (SNPs) between Sunhwanokdu and Gyeonggijaerae 5, defining the lowest and highest depths as 5 and 100, respectively, and the sequence quality as 100. Of the 305,504 single-base changes identified, 40,503 SNPs were considered heterozygous in Gyeonggijaerae 5. Among the remaining 265,001 SNPs, 65.9 % (174,579 cases) were transitions and 34.1 % (90,422 cases) were transversions. For SNP validation, a total of 42 SNPs were chosen among Sunhwanokdu contigs longer than 10 kb and sharing at least 80 % sequence identity with common bean expressed sequence tags as determined with est2genome. Using seven mungbean cultivars from various origins in addition to Sunhwanokdu and Gyeonggijaerae 5, most of the SNPs identified by bioinformatics tools were confirmed by Sanger sequencing. These genome-wide SNP markers could enrich the current molecular resources and might be of value for the construction of a mungbean genetic map and the investigation of genetic diversity.

Revisiting an important component of plant genomes: microsatellites

Microsatellites are some of the most highly variable repetitive DNA tracts in genomes. Few studies focus on whether the characteristic instability of microsatellites is linked to phenotypic effects in plants. We summarise recent data to investigate how microsatellite variations affect gene expression and hence phenotype. We discuss how the basic characteristics of microsatellites may contribute to phenotypic effects. In summary, microsatellites in plants are universal and highly mutable, they coexist and coevolve with transposable elements, and are under selective pressure. The number of motif nucleotides, the type of motif and transposon activity all contribute to the nonrandom generation and decay of microsatellites, and to conservation and distribution biases. Although microsatellites are generated by accident, they mature through responses to environmental change before final decay. This process is mediated by organism adjustment mechanisms, which maintain a balance between birth versus death and growth versus decay in microsatellites. Close relationships also exist between the physical structure, variation and functionality of microsatellites: in most plant species, sequences containing microsatellites are associated with catalytic activity and binding functions, are expressed in the membrane and organelles, and participate in the developmental and metabolic processes. Microsatellites contribute to genome structure and functional plasticity, and may be considered to promote species evolution in plants in response to environmental changes. In conclusion, the generation, loss, functionality and evolution of microsatellites can be related to plant gene expression and functional alterations. The effect of microsatellites on phenotypic variation may be as significant in plants as it is in animals.

Development and validation of microsatellite markers for Brachiaria ruziziensis obtained by partial genome assembly of Illumina single-end reads

Background

Brachiaria ruziziensis is one of the most important forage species planted in the tropics. The application of genomic tools to aid the selection of superior genotypes can provide support to B. ruziziensis breeding programs. However, there is a complete lack of information about the B. ruziziensis genome. Also, the availability of genomic tools, such as molecular markers, to support B. ruziziensis breeding programs is rather limited. Recently, next-generation sequencing technologies have been applied to generate sequence data for the identification of microsatellite regions and primer design. In this study, we present a first validated set of SSR markers for Brachiaria ruziziensis, selected from a de novo partial genome assembly of single-end Illumina reads.

Results

A total of 85,567 perfect microsatellite loci were detected in contigs with a minimum 10X coverage. We selected a set of 500 microsatellite loci identified in contigs with minimum 100X coverage for primer design and synthesis, and tested a subset of 269 primer pairs, 198 of which were polymorphic on 11 representative B. ruziziensis accessions. Descriptive statistics for these primer pairs are presented, as well as estimates of marker transferability to other relevant brachiaria species. Finally, a set of 11 multiplex panels containing the 30 most informative markers was validated and proposed for B. ruziziensis genetic analysis.

Conclusions

We show that the detection and development of microsatellite markers from genome assembled Illumina single-end DNA sequences is highly efficient. The developed markers are readily suitable for genetic analysis and marker assisted selection of Brachiaria ruziziensis. The use of this approach for microsatellite marker development is promising for species with limited genomic information, whose breeding programs would benefit from the use of genomic tools. To our knowledge, this is the first set of microsatellite markers developed for this important species.