Genetic Diversity Analysis of Faba Bean (Vicia faba L.) Based on EST-SSR Markers

Molecular genetic diversity and linkage disequilibrium structure of the Egyptian faba bean using Single Primer Enrichment Technology (SPET)

Faba bean is an important legume crop. The genetic diversity among faba bean genotypes is very important for the genetic improvement of target traits. A set of 128 fab bean genotypes that are originally from Egypt were used in this study to investigate the genetic diversity and population structure. The 128 genotypes were genotyped using the Single Primer Enrichment Technology (SPET) by which a set of 6759 SNP markers were generated after filtration. The SNP markers were distributed on all chromosomes with a range extending from 822 (Chr. 6) to 1872 (Chr.1). The SNP markers had wide ranges of polymorphic information content (PIC), gene diversity (GD), and minor allele frequency. The analysis of population structure divided the Egyptian faba bean population into five subpopulations. Considerable genetic distance was found among all genotypes, ranging from 0.1 to 0.4. The highly divergent genotype was highlighted in this study and the genetic distance among genotypes ranged from 0.1 and 0.6. Moreover, the structure of linkage disequilibrium was studied, and the analysis revealed a low level of LD in the Egyptian faba bean population. A slow LD decay at the genomic and chromosomal levels was observed. Interestingly, the distribution of haplotype blocks was presented in each chromosome and the number of haplotype block ranged from 65 (Chr. 4) to 156 (Chr. 1). Migration and genetic drift are the main reasons for the low LD in the Egyptian faba bean population. The results of this study shed light on the possibility of the genetic improvement of faba bean crop in Egypt and conducting genetic association analyses to identify candidate genes associated with target traits (e.g. protein content, grain yield, etc.) in this panel.

A detailed study on genetic diversity, antioxidant machinery, and expression profile of drought-responsive genes in rice genotypes exposed to artificial osmotic stress

Seasonal variations in rainfall patterns, particularly during sowing, early growing season, and flowering, drastically affect rice production in northeastern India. However, sensitivity to drought stress is genotype-specific. Since 80% of the land in this region is used for rice production, it is crucial to understand how they have adapted to water stress. This study evaluated 112 rice genotypes grown in NE India for seed germination percentage and seedling development under PEG-mediated drought stress. Among the rice genotype, Sahbhagi dhan, RCPL-1-82, Bhalum-3 and RCPL-1-128 showed drought-tolerant traits, while Ketaki Joha, Chakhao, Chandan, RCPL-1-185 and IR-64 were the most drought-sensitive rice genotypes. Drought-tolerant rice also showed significantly higher seed germination potential, proline content, antioxidant activity and expression of drought-responsive genes than drought-sensitive rice genotypes. A similar expression pattern of genes was also observed in the rice genotype treated with a 50% water deficit in pot culture. In addition, drought stress reduced the pollen fertility and yield per plant in sensitive rice genotypes. Molecular markers associated with drought stress were also used to characterize genetic diversity among the rice genotypes studied.

AFLP-Based Analysis of Variation and Population Structure in Mutagenesis Induced Faba Bean

Genetic diversity enrichment is urgently necessary to develop climate-resilient faba bean cultivars. The present study aimed to measure the enrichment of genetic diversity and changes in the population structure of faba bean, following induced mutagenesis. 120 samples, including 116 M2 mutant plants, generated by exposing the ILB4347 accession to four mutagen treatments (25 and 50 Gray gamma radiation and 0.01%, and 0.05% diethyl sulfate) and four reference genotypes were characterized using 11 amplified fragment length polymorphism (AFLP) primer combinations. The AFLP markers generated 1687 polymorphic alleles, including 756 alleles (45%) that were detected infrequently (f ≤ 0.1). The total allele count of the mutant plants ranged from 117 to 545. We observed a wide range of banding patterns and counts among the mutant plants, showing the high genetic diversity induced by mutation. Mutations also changed the population structure, by altering 31.78% of the total membership coefficient (Q). Although mutations changed the population structure, Nei’s genetic distance showed that the mutant population remained closely related to its control parent. This is the first report examining genetic diversity and population changes in faba bean mutant populations and, thus, could facilitate the application of induced mutagenesis during faba bean breeding.

Characterization of drought stress-responsive root transcriptome of faba bean (Vicia faba L.) using RNA sequencing

Drought and salinity are the major factors that limit the faba bean (Vicia faba L.) production worldwide. The aim of this study is to identify the water stress differentially expressed genes (DEGs) through the root transcriptome analyses of the drought-tolerant Hassawi 2 genotype at vegetative and flowering stages. A total of 624.8 M high-quality Illumina reads were generated and assembled into 198,155 all-unigenes with a mean length of 738 bp and an N50 length of 1347 bp. Among all-unigenes, 78,262 were assigned to non-redundant (Nr), 66,254 to nucleotide (Nt), 54,034 to KEGG, and 43,913 to gene ontology (GO) annotations. A total of 36,834 and 35,510 unigenes were differentially expressed at the vegetative and flowering stages of Hassawi 2 under drought stress, respectively. The majority of unigenes were down-regulated at both developmental stages. However, the number of genes up-regulated (15,366) at the flowering stage exceeded the number of those up-regulated (14,097) at the vegetative stage, and the number of genes down-regulated (20,144) at the flowering stage was smaller than the number of those down-regulated (22,737) at the vegetative stage. The drought stress-responsive differentially expressed unigenes coded for various regulatory proteins, including protein kinases and phosphatases, transcription factors and plant hormones and functional proteins including enzymes for osmoprotectant, detoxification and transporters were differentially expressed, most of which were largely up-regulated. Moreover, a substantial proportion of the DEGs identified in this study were novel, most exhibited a significant change in their expression levels under water stress, making them an unexploited resource that might control specific responses to drought stress in the faba bean. Finally, qRT-PCR results were found almost consistent with the results of next-generation sequencing. Our data will help in understanding the drought tolerance mechanisms in plants and will provide resources for functional genomics.

Molecular genetic diversity and population structure of Ethiopian white lupin landraces: Implications for breeding and conservation

White lupin is one of the four economically important species of the Lupinus genus and is an important grain legume in the Ethiopian farming system. However, there has been limited research effort to characterize the Ethiopian white lupin landraces. Fifteen polymorphic simple sequence repeat (SSR) markers were used to assess the genetic diversity and population structure of 212 Ethiopian white lupin (Lupinus albus) landraces and two genotypes from different species (Lupinus angustifolius and Lupinus mutabilis) were used as out-group. The SSR markers revealed 108 different alleles, 98 of them from 212 landraces and 10 from out-group genotypes, with an average of 6.5 alleles per locus. The average gene diversity was 0.31. Twenty eight landraces harbored one or more private alleles from the total of 28 private alleles identified in the 212 white lupin accessions. Seventy-seven rare alleles with a frequency of less than 5% were identified and accounted for 78.6% of the total alleles detected. Analysis of molecular variance (AMOVA) showed that 92% of allelic diversity was attributed to individual accessions within populations while only 8% was distributed among populations. At 70% similarity level, the UPGMA dendrogram resulted in the formation of 13 clusters comprised of 2 to 136 landraces, with the out-group genotypes and five landraces remaining distinct and ungrouped. Population differentiation and genetic distance were relatively high between Gondar and Ethiopian white lupin populations collected by Australians. A model-based population structure analysis divided the white lupin landraces into two populations. All Ethiopian white lupin landrace populations, except most of the landraces collected by Australians (77%) and about 44% from Awi, were grouped together with significant admixtures. The study also suggested that 34 accessions, as core collections, were sufficient to retain 100% of SSR diversity. These accessions (core G-34) represent 16% of the whole 212 Ethiopian white lupin accessions and populations from West Gojam, Awi and Australian collections contributed more accessions to the core collection.

De novo transcriptomic analysis of cowpea (Vigna unguiculata L. Walp.) for genic SSR marker development

Background

Cowpea [Vigna unguiculata (L.) Walp.] is one of the most important legumes in tropical and semi-arid regions. However, there is relatively little genomic information available for genetic research on and breeding of cowpea. The objectives of this study were to analyse the cowpea transcriptome and develop genic molecular markers for future genetic studies of this genus.

Results

Approximately 54 million high-quality cDNA sequence reads were obtained from cowpea based on Illumina paired-end sequencing technology and were de novo assembled to generate 47,899 unigenes with an N50 length of 1534 bp. Sequence similarity analysis revealed 36,289 unigenes (75.8%) with significant similarity to known proteins in the non-redundant (Nr) protein database, 23,471 unigenes (49.0%) with BLAST hits in the Swiss-Prot database, and 20,654 unigenes (43.1%) with high similarity in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Further analysis identified 5560 simple sequence repeats (SSRs) as potential genic molecular markers. Validating a random set of 500 SSR markers yielded 54 polymorphic markers among 32 cowpea accessions.

Conclusions

This transcriptomic analysis of cowpea provided a valuable set of genomic data for characterizing genes with important agronomic traits in Vigna unguiculata and a new set of genic SSR markers for further genetic studies and breeding in cowpea and related Vigna species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0531-5) contains supplementary material, which is available to authorized users.

Genetic variability in Tunisian populations of faba bean (Vicia faba L. var. major) assessed by morphological and SSR markers

The genetic diversity of 21 faba bean populations was examined using morphological and molecular markers. DNA was extracted from 189 individuals and 8 microsatellite markers were genotyped individually in these 21 populations. A total of 53 alleles were obtained in all populations, with an average of 6.62 alleles per locus. The expected and observed heterozygosity was 0.38 and 0.62 respectively. The average polymorphism index content of SSR markers was 0.61, ranging from 0.31 to 0.81. The unweighted pair group method with arithmetic mean dendrogram clustered all the populations into two groups, each for them subdivided into 3 sub-groups according to geographical origin. Morphological variation showed that the populations were not grouped according to their geographical origin. Therefore, patterns of differentiation of morphological traits did not coincide with molecular differentiation, indicating that morphological variation does not reflect genetic subdivision in studied faba bean populations. Analysis of molecular variance revealed high levels of genetic variation (83%) within population and provides a good base for designing genetic improvement programs. The result of Principal Component Analysis (PCA) revealed that three dimensional principal components (PC1, PC2 and PC3) contributed 40.56% of the total variability and accounted with values of 20.64, 11.22 and 8.70%, respectively. Cluster analysis based on PCA indicated three separate groups of populations. The genetic relationships found between the 21 populations samples were the same in both the PCA and STRUCTURE analysis which support the results observed. These data may serve as a foundation for the development of faba bean breeding programs.

Advances in Faba Bean Genetics and Genomics

Vicia faba L, is a globally important grain legume whose main centers of diversity are the Fertile Crescent and Mediterranean basin. Because of its small number (six) of exceptionally large and easily observed chromosomes it became a model species for plant cytogenetics the 70s and 80s. It is somewhat ironic therefore, that the emergence of more genomically tractable model plant species such as Arabidopsis and Medicago coincided with a marked decline in genome research on the formerly favored plant cytogenetic model. Thus, as ever higher density molecular marker coverage and dense genetic and even complete genome sequence maps of key crop and model species emerged through the 1990s and early 2000s, genetic and genome knowledge of Vicia faba lagged far behind other grain legumes such as soybean, common bean and pea. However, cheap sequencing technologies have stimulated the production of deep transcriptome coverage from several tissue types and numerous distinct cultivars in recent years. This has permitted the reconstruction of the faba bean meta-transcriptome and has fueled development of extensive sets of Simple Sequence Repeat and Single Nucleotide Polymorphism (SNP) markers. Genetics of faba bean stretches back to the 1930s, but it was not until 1993 that DNA markers were used to construct genetic maps. A series of Random Amplified Polymorphic DNA-based genetic studies mainly targeted at quantitative loci underlying resistance to a series of biotic and abiotic stresses were conducted during the 1990's and early 2000s. More recently, SNP-based genetic maps have permitted chromosome intervals of interest to be aligned to collinear segments of sequenced legume genomes such as the model legume Medicago truncatula, which in turn opens up the possibility for hypotheses on gene content, order and function to be translated from model to crop. Some examples of where knowledge of gene content and function have already been productively exploited are discussed. The bottleneck in associating genes and their functions has therefore moved from locating gene candidates to validating their function and the last part of this review covers mutagenesis and genetic transformation, two complementary routes to validating gene function and unlocking novel trait variation for the improvement of this important grain legume.

Development of Gene-Based SSR Markers in Rice Bean (Vigna umbellata L.) Based on Transcriptome Data

Rice bean (Vigna umbellata (Thunb.) Ohwi & Ohashi) is a warm season annual legume mainly grown in East Asia. Only scarce genomic resources are currently available for this legume crop species and no simple sequence repeat (SSR) markers have been specifically developed for rice bean yet. In this study, approximately 26 million high quality cDNA sequence reads were obtained from rice bean using Illumina paired-end sequencing technology and assembled into 71,929 unigenes with an average length of 986 bp. Of these unigenes, 38,840 (33.2%) showed significant similarity to proteins in the NCBI non-redundant protein and nucleotide sequence databases. Furthermore, 30,170 (76.3%) could be classified into gene ontology categories, 25,451 (64.4%) into Swiss-Prot categories and 21,982 (55.6%) into KOG database categories (E-value < 1.0E-5). A total of 9,301 (23.5%) were mapped onto 118 pathways using the Kyoto Encyclopedia of Genes and Genome (KEGG) pathway database. A total of 3,011 genic SSRs were identified as potential molecular markers. AG/CT (30.3%), AAG/CTT (8.1%) and AGAA/TTCT (20.0%) are the three main repeat motifs. A total of 300 SSR loci were randomly selected for validation by using PCR amplification. Of these loci, 23 primer pairs were polymorphic among 32 rice bean accessions. A UPGMA dendrogram revealed three major clusters among 32 rice bean accessions. The large number of SSR-containing sequences and genic SSRs in this study will be valuable for the construction of high-resolution genetic linkage maps, association or comparative mapping and genetic analyses of various Vigna species.

Development and Validation of EST-SSR Markers from the Transcriptome of Adzuki Bean (Vigna angularis)

The adzuki bean (Vigna angularis (Ohwi) Ohwi and Ohashi) is an important grain legume of Asia. It is cultivated mainly in China, Japan and Korea. Despite its importance, few genomic resources are available for molecular genetic research of adzuki bean. In this study, we developed EST-SSR markers for the adzuki bean through next-generation sequencing. More than 112 million high-quality cDNA sequence reads were obtained from adzuki bean using Illumina paired-end sequencing technology, and the sequences were de novo assembled into 65,950 unigenes. The average length of the unigenes was 1,213 bp. Among the unigenes, 14,547 sequences contained a unique simple sequence repeat (SSR) and 3,350 sequences contained more than one SSR. A total of 7,947 EST-SSRs were identified as potential molecular markers, with mono-nucleotide A/T repeats (99.0%) as the most abundant motif class, followed by AG/CT (68.4%), AAG/CTT (30.0%), AAAG/CTTT (26.2%), AAAAG/CTTTT (16.1%), and AACGGG/CCCGTT (6.0%). A total of 500 SSR markers were randomly selected for validation, of which 296 markers produced reproducible amplicons with 38 polymorphic markers among the 32 adzuki bean genotypes selected from diverse geographical locations across China. The large number of SSR-containing sequences and EST-SSR markers will be valuable for genetic analysis of the adzuki bean and related Vigna species.

Assessment of genetic diversity among faba bean genotypes using agro-morphological and molecular markers

Forty faba bean (Vicia faba L.) genotypes were evaluated for their agro-morphological performance and molecular diversity under Central Region of Saudi Arabia conditions during 2010–11 and 2011–12 seasons. Field performance results showed that faba genotypes exhibited a significant amount of variation for their agro-morphological studied parameters. Giza40 recorded the tallest genotype (139.5 cm), highest number of seeds per plants (100.8), and the highest seed yield per plant (70.8 g). The best performing genotypes were Giza40, FLIP03-014FB, Gazira1 and Goff1. Genetic variability among genotypes was determined using Sequence Related Amplified Polymorphism (SRAP) and Amplified Fragment Length Polymorphism (AFLP) markers. A total of 183 amplified fragments (alleles) and 1758 polymorphic fragments (bands) in SRAP and 202 alleles and 716 bands in AFLP were obtained using six SRAP and four AFLP primer combinations respectively. Polymorphism information content (PIC) values for AFLP and SRAP markers were higher than 0.8, indicating the existence of a considerable amount of genetic diversity among faba tested genotypes. The UPGMA based clustering of faba genotypes was largely based on origin and/or genetic background. Result of cluster analysis based on SRAP showed weak and not significant correlation while, it was highly significant based on AFLP analysis with agro-morphological characters (r = 0.01, p > 0.54 and r = 0.26, p < 0.004 respectively). Combined SRAP and AFLP markers proved to be significantly useful for genetic diversity assessment at molecular level. They exhibited high discrimination power, and were able to distinguish the faba bean genotypes with high efficiency and accuracy levels.

Determination of the genetic diversity of vegetable soybean [Glycine max (L.) Merr.] using EST-SSR markers

The development of expressed sequence tag-derived simple sequence repeats (EST-SSRs) provided a useful tool for investigating plant genetic diversity. In the present study, 22 polymorphic EST-SSRs from grain soybean were identified and used to assess the genetic diversity in 48 vegetable soybean accessions. Among the 22 EST-SSR loci, tri-nucleotides were the most abundant repeats, accounting for 50.00% of the total motifs. GAA was the most common motif among tri-nucleotide repeats, with a frequency of 18.18%. Polymorphic analysis identified a total of 71 alleles, with an average of 3.23 per locus. The polymorphism information content (PIC) values ranged from 0.144 to 0.630, with a mean of 0.386. Observed heterozygosity (Ho) values varied from 0.0196 to 1.0000, with an average of 0.6092, while the expected heterozygosity (He) values ranged from 0.1502 to 0.6840, with a mean value of 0.4616. Principal coordinate analysis and phylogenetic tree analysis indicated that the accessions could be assigned to different groups based to a large extent on their geographic distribution, and most accessions from China were clustered into the same groups. These results suggest that Chinese vegetable soybean accessions have a narrow genetic base. The results of this study indicate that EST-SSRs from grain soybean have high transferability to vegetable soybean, and that these new markers would be helpful in taxonomy, molecular breeding, and comparative mapping studies of vegetable soybean in the future.

Molecular diversity assessment using Sequence Related Amplified Polymorphism (SRAP) Markers in Vicia faba L

Sequence-related amplified polymorphism (SRAP) markers were used to assess the genetic diversity and relationship among 58 faba bean (Vicia faba L.) genotypes. Fourteen SRAP primer combinations amplified a total of 1036 differently sized well-resolved peaks (fragments), of which all were polymorphic with a 0.96 PIC value and discriminated all of the 58 faba bean genotypes. An average pairwise similarity of 21% was revealed among the genotypes ranging from 2% to 65%. At a similarity of 28%, UPGMA clustered the genotypes into three main groups comprising 78% of the genotypes. The local landraces and most of the Egyptian genotypes in addition to the Sudan genotypes were grouped in the first main cluster. The advanced breeding lines were scattered in the second and third main clusters with breeding lines from the ICARDA and genotypes introduced from Egypt. At a similarity of 47%, all the genotypes formed separated clusters with the exceptions of Hassawi 1 and Hassawi 2. Group analysis of the genotypes according to their geographic origin and type showed that the landraces were grouped according to their origin, while others were grouped according to their seed type. To our knowledge, this is the first application of SRAP markers for the assessment of genetic diversity in faba bean. Such information will be useful to determine optimal breeding strategies to allow continued progress in faba bean breeding.