Development and molecular characterization of 55 novel polymorphic cDNA-SSR markers in faba bean (Vicia faba L.) using 454 pyrosequencing

Genome-wide discovery of di-nucleotide SSR markers based on whole genome re-sequencing data of Cicer arietinum L. and Cicer reticulatum Ladiz

Simple sequence repeats (SSRs) are valuable genetic markers due to their co-dominant inheritance, multi-allelic and reproducible nature. They have been largely used for exploiting genetic architecture of plant germplasms, phylogenetic analysis, and mapping studies. Among the SSRs, di-nucleotide repeats are the most frequent of the simple repeats distributed throughout the plant genomes. In present study, we aimed to discover and develop di-nucleotide SSR markers by using the whole genome re-sequencing (WGRS) data from Cicer arietinum L. and C. reticulatum Ladiz. A total of 35,329 InDels were obtained in C. arietinum, whereas 44,331 InDels in C. reticulatum. 3387 InDels with 2 bp length were detected in C. arietinum, there were 4704 in C. reticulatum. Among 8091 InDels, 58 di-nucleotide regions that were polymorphic between two species were selected and used for validation. We tested primers for evaluation of genetic diversity in 30 chickpea genotypes including C. arietinum, C. reticulatum, C. echinospermum P.H. Davis, C. anatolicum Alef., C. canariense A. Santos & G.P. Lewis, C. microphyllum Benth., C. multijugum Maesen, C. oxyodon Boiss. & Hohen. and C. songaricum Steph ex DC. A total of 244 alleles were obtained for 58 SSR markers giving an average of 2.36 alleles per locus. The observed heterozygosity was 0.08 while the expected heterozygosity was 0.345. Polymorphism information content was found to be 0.73 across all loci. Phylogenetic tree and principal coordinate analysis clearly divided the accessions into four groups. The SSR markers were also evaluated in 30 genotypes of a RIL population obtained from an interspecific cross between C. arietinum and C. reticulatum. Chi-square (χ²) test revealed an expected 1:1 segregation ratio in the population. These results demonstrated the success of SSR identification and marker development for chickpea with the use of WGRS data. The newly developed 58 SSR markers are expected to be useful for chickpea breeders.

VfODB: a comprehensive database of ESTs, EST-SSRs, mtSSRs, microRNA-target markers and genetic maps in Vicia faba

Faba bean (Vicia faba) is an essential food and fodder legume crop worldwide due to its high content of proteins and fibres. Molecular markers tools represent an invaluable tool for faba bean breeders towards rapid crop improvement. Although there have historically been few V. faba genome resources available, several transcriptomes and mitochondrial genome sequence data have been released. These data in addition to previously developed genetic linkage maps represent a great resource for developing functional markers and maps that can accelerate the faba bean breeding programmes. Here, we present the Vicia faba Omics database (VfODB) as a comprehensive database integrating germplasm information, expressed sequence tags (ESTs), expressed sequence tags-simple sequence repeats (EST-SSRs), and mitochondrial-simple sequence repeats (mtSSRs), microRNA-target markers and genetic maps in faba bean. In addition, KEGG pathway-based markers and functional maps are integrated as a novel class of annotation-based markers/maps. Collectively, we developed 31 536 EST markers, 9071 EST-SSR markers and 3023 microRNA-target markers based on V. faba RefTrans V2 mining. By mapping 7940 EST and 2282 EST-SSR markers against the KEGG pathways database we successfully developed 107 functional maps. Also, 40 mtSSR markers were developed based on mitochondrial genome mining. On the data curation level, we retrieved 3461 markers representing 12 types of markers (CAPS, EST, EST-SSR, Gene marker, INDEL, Isozyme, ISSR, RAPD, SCAR, RGA, SNP and SSR), which mapped across 18 V. faba genetic linkage maps. VfODB provides two user-friendly tools to identify, classify SSR motifs and in silico amplify their targets. VfODB can serve as a powerful database and helpful platform for faba bean research community as well as breeders interested in Genomics-Assisted Breeding.

Mining and Development of Novel SSR Markers Using Next Generation Sequencing (NGS) Data in Plants

Microsatellites, or simple sequence repeats (SSRs), are one of the most informative and multi-purpose genetic markers exploited in plant functional genomics. However, the discovery of SSRs and development using traditional methods are laborious, time-consuming, and costly. Recently, the availability of high-throughput sequencing technologies has enabled researchers to identify a substantial number of microsatellites at less cost and effort than traditional approaches. Illumina is a noteworthy transcriptome sequencing technology that is currently used in SSR marker development. Although 454 pyrosequencing datasets can be used for SSR development, this type of sequencing is no longer supported. This review aims to present an overview of the next generation sequencing, with a focus on the efficient use of de novo transcriptome sequencing (RNA-Seq) and related tools for mining and development of microsatellites in plants.

Genetic diversity of Xanthoceras sorbifolium bunge germplasm using morphological traits and microsatellite molecular markers

Xanthoceras sorbifolium Bunge has great potential for producing biodiesel. In order to select and evaluate appropriate germplasm to produce biodiesel, we analyzed the genetic diversity of Xanthoceras sorbifolium Bunge germplasm based on morphological traits and simple sequence repeats (SSRs) in this study. Fifty-six germplasm samples were evaluated using nine morphological traits and 23 SSR loci. Significant differences among germplasms were observed in eight morphological characters. The SSR markers analysis showed high genetic diversity among the germplasms. All SSRs had polymorphisms, and we detected 77 alleles in total. The number of alleles at each locus ranged from two to six, averaging 3.35 per marker. The polymorphic information content values ranged from 0.36 to 0.61, averaging 0.49. Expected heterozygosity, observed heterozygosity, and Shannon's information index calculations detected large genetic variations among germplasms. The high average number of alleles per locus and the allelic diversity observed in the set of genotypes analyzed indicated that the genetic base of this species is relatively wide. Thus, microsatellite markers can be used to efficiently distinguish Xanthoceras sorbifolium Bunge germplasms and assess their genetic diversity. Hundred-grain weight and lateral diameter were positively correlated with monounsaturated fatty acids and depended on genotype. These results suggest that seeds with higher hundred-grain weight and lateral diameter could be more suitable to produce biodiesel. Our data will lay a foundation for selecting appropriate germplasm to produce biodiesel based on seed phenotype and will contribute to the conservation and management of this important plant genetic resource.

Relationship between geographical origin, seed size and genetic diversity in faba bean (Vicia faba L.) as revealed by SSR markers

Faba bean (Vicia faba L.) is an important legume species because of its high protein and starch content. Broad bean can be grown in different climatic conditions and is an ideal rotation crop because of the nitrogen fixing bacteria in its roots. In this work, 255 faba bean germplasm accessions were characterized using 32 SSR primers which yielded 302 polymorphic fragments. According to the results, faba bean individuals were divided into two main groups based on the neighbor-joining algorithm (r = 0.91) with some clustering based on geographical origin as well as seed size. Population structure was also determined and agreed with the dendrogram analysis in splitting the accessions into two subpopulations. Analysis of molecular variance (AMOVA) revealed high levels of within population genetic variation. Genetic similarity and geographical proximity were related with separation of European accessions from African and Asian ones. Interestingly, there was no significant difference between landrace (38%) and cultivar (40%) diversity indicating that genetic variability has not yet been lost due to breeding. A total of 44 genetically well-characterized faba bean individuals were selected for a core collection to be further examined for yield and nutritional traits.

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.

Genetic Diversity and Population Structure in Vicia faba L. Landraces and Wild Related Species Assessed by Nuclear SSRs

Faba bean (Vicia faba L.) is a facultative cross-pollinating legume crop with a great importance for food and feed due to its high protein content as well as the important role in soil fertility and nitrogen fixation. In this work we evaluated genetic diversity and population structure of faba bean accessions from the Western Mediterranean basin and wild related species. For that purpose we screened 53 V. faba, 2 V. johannis and 7 V. narbonensis accessions from Portugal, Spain and Morocco with 28 faba bean Single Sequence Repeats (SSR). SSR genotyping showed that the number of alleles detected per locus for the polymorphic markers ranged between 2 and 10, with Polymorphic Information Content (PIC) values between 0.662 and 0.071, and heterozygosity (H_O) between 0–0.467. Heterozygosity and inbreeding coefficient levels indicate a higher level of inbreeding in wild related species than in cultivated Vicia. The analysis of molecular variance (AMOVA) showed a superior genetic diversity within accessions than between accessions even from distant regions. These results are in accordance to population structure analysis showing that individuals from the same accession can be genetically more similar to individuals from far away accessions, than from individuals from the same accession. In all three levels of analysis (whole panel of cultivated and wild accessions, cultivated faba bean accessions and Portuguese accessions) no population structure was observed based on geography or climatic factors. Differences between V. narbonensis and V. johannis are undetectable although these wild taxa are clearly distinct from V. faba accessions. Thus, a limited gene flow occurred between cultivated accessions and wild relatives. Contrastingly, the lack of population structure seems to indicate a high degree of gene flow between V. faba accessions, possibly explained by the partially allogamous habit in association with frequent seed exchange/introduction.

Identification and Verification of QTL Associated with Frost Tolerance Using Linkage Mapping and GWAS in Winter Faba Bean

Frost stress is one of the abiotic stresses that causes a significant reduction in winter faba bean yield in Europe. The main objective of this work is to genetically improve frost tolerance in winter faba bean by identifying and validating QTL associated with frost tolerance to be used in marker-assisted selection (MAS). Two different genetic backgrounds were used: a biparental population (BPP) consisting of 101 inbred lines, and 189 genotypes from single seed descent (SSD) from the Gottingen Winter bean Population (GWBP). All experiments were conducted in a frost growth chamber under controlled conditions. Both populations were genotyped using the same set of 189 SNP markers. Visual scoring for frost stress symptoms was used to define frost tolerance in both populations. In addition, leaf fatty acid composition (FAC) and proline content were analyzed in BPP as physiological traits. QTL mapping (for BPP) and genome wide association studies (for GWBP) were performed to detect QTL associated with frost tolerance. High genetic variation between genotypes, and repeatability estimates, were found for all traits. QTL mapping and GWAS identified new putative QTL associated with promising frost tolerance and related traits. A set of 54 SNP markers common in both genetic backgrounds showed a high genetic diversity with polymorphic information content (PIC) ranging from 0.31 to 0.37 and gene diversity ranging from 0.39 to 0.50. This indicates that these markers may be polymorphic for many faba bean populations. Five SNP markers showed a significant marker-trait association with frost tolerance and related traits in both populations. Moreover, synteny analysis between Medicago truncatula (a model legume) and faba bean genomes was performed to identify candidate genes for these markers. Collinearity was evaluated between the faba bean genetic map constructed in this study and the faba bean consensus map, resulting in identifying possible genomic regions in faba bean which may control frost tolerance genes. The two genetic backgrounds were useful in detecting new variation for improving frost tolerance in winter faba bean. Of the five validated SNP markers, one (VF_Mt3g086600) was found to be associated with frost tolerance and FAC in both populations. This marker was also associated with winter hardiness and high yield in earlier studies. This marker is located in a gene of unknown function.

De Novo Sequencing and Analysis of the Safflower Transcriptome to Discover Putative Genes Associated with Safflor Yellow in Carthamus tinctorius L

Safflower (Carthamus tinctorius L.), an important traditional Chinese medicine, is cultured widely for its pharmacological effects, but little is known regarding the genes related to the metabolic regulation of the safflower’s yellow pigment. To investigate genes related to safflor yellow biosynthesis, 454 pyrosequencing of flower RNA at different developmental stages was performed, generating large databases.In this study, we analyzed 454 sequencing data from different flowering stages in safflower. In total, 1,151,324 raw reads and 1,140,594 clean reads were produced, which were assembled into 51,591 unigenes with an average length of 679 bp and a maximum length of 5109 bp. Among the unigenes, 40,139 were in the early group, 39,768 were obtained from the full group and 28,316 were detected in both samples. With the threshold of “log2 ratio ≥ 1”, there were 34,464 differentially expressed genes, of which 18,043 were up-regulated and 16,421 were down-regulated in the early flower library. Based on the annotations of the unigenes, 281 pathways were predicted. We selected 12 putative genes and analyzed their expression levels using quantitative real time-PCR. The results were consistent with the 454 sequencing results. In addition, the expression of chalcone synthase, chalcone isomerase and anthocyanidin synthase, which are involved in safflor yellow biosynthesis and safflower yellow pigment (SYP) content, were analyzed in different flowering periods, indicating that their expression levels were related to SYP synthesis. Moreover, to further confirm the results of the 454 pyrosequencing, full-length cDNA of chalcone isomerase (CHI) and anthocyanidin synthase (ANS) were cloned from safflower petal by RACE (Rapid-amplification of cDNA ends) method according to fragment of the transcriptome.

Isolation and characterization of microsatellite markers for Jasminum sambac (Oleaceae) using Illumina shotgun sequencing

PREMISE OF THE STUDY

Microsatellite markers of Jasminum sambac (Oleaceae) were isolated to investigate wild germplasm resources and provide markers for breeding.

METHODS AND RESULTS

Illumina sequencing was used to isolate microsatellite markers from the transcriptome of J. sambac. A total of 1322 microsatellites were identified from 49,772 assembled unigenes. One hundred primer pairs were randomly selected to verify primer amplification efficiency. Out of these tested primer pairs, 31 were successfully amplified: 18 primer pairs yielded a single allele, seven exhibited fixed heterozygosity with two alleles, and only six displayed polymorphisms.

CONCLUSIONS

This study obtained the first set of microsatellite markers for J. sambac, which will be helpful for the assessment of wild germplasm resources and the development of molecular marker-assisted breeding.

Association mapping for frost tolerance using multi-parent advanced generation inter-cross (MAGIC) population in faba bean (Vicia faba L.)

A multi-parent advanced generation inter-cross (MAGIC) derived from 11 founder lines in faba bean was used in this study to identify quantitative trait loci (QTL) for frost tolerance traits using the association mapping method with 156 SNP markers. This MAGIC population consists of a set of 189 genotypes from the Göttingen Winter Bean Population. The association panel was tested in two different experiments, i.e. a frost and a hardening experiment. Six morphological traits, leaf fatty acid composition, relative water content in shoots were scored in this study. The genotypes presented a large genetic variation for all traits that were highly heritable after frost and after hardening. High phenotypic significant correlations were established between traits. The principal coordinates analysis resulted in no clear structure in the current population. Association mapping was performed using a general linear model and mixed linear model with kinship. A False discovery rate of 0.20 (and 0.05) was used to test the significance of marker-trait association. As a result, many putative QTLs for 13 morphological and physiological traits were detected using both models. The results reveal that QTL mapping by association analysis is a powerful method of detecting the alleles associated with frost tolerance in the winter faba bean which can be used in accelerating breeding programs.

Cross-amplification of Vicia sativa subsp. sativa microsatellites across 22 other Vicia species

The temperate and herbaceous genus Vicia L. is a member of the legume tribe Fabeae of the subfamily Papilionoideae. The genus Vicia comprises 166 annual or perennial species distributed mainly in Europe, Asia, and North America, but also extending to the temperate regions of South America and tropical Africa. The use of simple sequence repeat (SSR) markers for Vicia species has not been investigated as extensively as for other crop species. In this study, we assessed the potential for cross-species amplification of cDNA microsatellite markers developed from common vetch (Vicia sativa subsp. sativa). For cross-species amplification of the SSRs, amplification was carried out with genomic DNA isolated from two to eight accessions of 22 different Vicia species. For individual species or subspecies, the transferability rates ranged from 33% for V. ervilia to 82% for V. sativa subsp. nigra with an average rate of 52.0%. Because the rate of successful SSR marker amplification generally correlates with genetic distance, these SSR markers are potentially useful for analyzing genetic relationships between or within Vicia species.

Hemolytic crisis in a G6PD-deficient infant after ingestion of pumpkin

A 8 month-old infant presented with acute onset of severe jaundice, anemia requiring transfusion and Glucose-6-Phosphate Dehydrogenase deficiency. The infant did not take drugs, he did not consume fava beans, but fava beans DNA was found on pumpkin he consumed the day before jaundice onset. This is the first case of hemolysis triggered by ingestion of food cross-contaminated with fava beans.

Exploiting Illumina sequencing for the development of 95 novel polymorphic EST-SSR markers in common vetch (Vicia sativa subsp. sativa)

The common vetch (Vicia sativa subsp. sativa), a self-pollinating and diploid species, is one of the most important annual legumes in the world due to its short growth period, high nutritional value, and multiple usages as hay, grain, silage, and green manure. The available simple sequence repeat (SSR) markers for common vetch, however, are insufficient to meet the developing demand for genetic and molecular research on this important species. Here, we aimed to develop and characterise several polymorphic EST-SSR markers from the vetch Illumina transcriptome. A total number of 1,071 potential EST-SSR markers were identified from 1025 unigenes whose lengths were greater than 1,000 bp, and 450 primer pairs were then designed and synthesized. Finally, 95 polymorphic primer pairs were developed for the 10 common vetch accessions, which included 50 individuals. Among the 95 EST-SSR markers, the number of alleles ranged from three to 13, and the polymorphism information content values ranged from 0.09 to 0.98. The observed heterozygosity values ranged from 0.00 to 1.00, and the expected heterozygosity values ranged from 0.11 to 0.98. These 95 EST-SSR markers developed from the vetch Illumina transcriptome could greatly promote the development of genetic and molecular breeding studies pertaining to in this species.

Isolation and characterization of microsatellite markers for Cotinus coggygria Scop. (Anacardiaceae) by 454 pyrosequencing

Cotinus coggygria Scop. (Anacardiaceae) is a deciduous shrub or small tree that is native to a large area covering from southern Europe, east across central Asia, and the Himalayas in northern China. Shotgun 454 pyrosequencing was used to develop microsatellite markers from the genome of C. coggygria. In this study, 349 microsatellite loci were identified from 40,074 individual sequence reads produced by one-sixteenth run, and primer pairs were designed for these loci. To test the primer amplification efficiency, 50 microsatellite primer pairs were tested across 12 individuals from two C. coggygria populations (Wuzhi Mountain: 36°30'N, 113°39'E; Tianlong Mountain: 37°42'N, 112°26'E). Among the 50 tested primer pairs, eight were found to be polymorphic. The average allele number of the microsatellites was 3.5 per locus, with a range from two to five. The inbreeding coefficient ranged from -0.478 to 0.222. The observed and expected heterozygosities varied from 0.167 to 0.750 and from 0.163 to 0.743, respectively. This set of markers is potentially useful for assessing the genetic diversity, as well as for understanding the population structure and phylogeographical and landscape genetic patterns, of C. coggygria.

SNP discovery and high-density genetic mapping in faba bean (Vicia faba L.) permits identification of QTLs for ascochyta blight resistance

Ascochyta blight, caused by the fungus Ascochyta fabae Speg., is a common and destructive disease of faba bean (Vicia faba L.) on a global basis. Yield losses vary from typical values of 35-40% to 90% under specific environmental conditions. Several sources of resistance have been identified and used in breeding programs. However, introgression of the resistance gene determinants into commercial cultivars as a gene pyramiding approach is reliant on selection of closely linked genetic markers. A total of 14,552 base variants were identified from a faba bean expressed sequence tag (EST) database, and were further quality assessed to obtain a set of 822 high-quality single nucleotide polymorphisms (SNPs). Sub-sets of 336 EST-derived simple sequence repeats (SSRs) and 768 SNPs were further used for high-density genetic mapping of a biparental faba bean mapping population (Icarus×Ascot) that segregates for resistance to ascochyta blight. The linkage map spanned a total length of 1216.8 cM with 12 linkage groups (LGs) and an average marker interval distance of 2.3 cM. Comparison of map structure to the genomes of closely related legume species revealed a high degree of conserved macrosynteny, as well as some rearrangements. Based on glasshouse evaluation of ascochyta blight resistance performed over two years, four genomic regions controlling resistance were identified on Chr-II, Chr-VI and two regions on Chr-I.A. Of these, one (QTL-3) may be identical with quantitative trait loci (QTLs) identified in prior studies, while the others (QTL-1, QTL-2 and QTL-4) may be novel. Markers in close linkage to ascochyta blight resistance genes identified in this study can be further validated and effectively implemented in faba bean breeding programs.

Use of EST-SSR markers for evaluating genetic diversity and fingerprinting celery (Apium graveolens L.) cultivars

Celery (Apium graveolens L.) is one of the most economically important vegetables worldwide, but genetic and genomic resources supporting celery molecular breeding are quite limited, thus few studies on celery have been conducted so far. In this study we made use of simple sequence repeat (SSR) markers generated from previous celery transcriptome sequencing and attempted to detect the genetic diversity and relationships of commonly used celery accessions and explore the efficiency of the primers used for cultivars identification. Analysis of molecular variance (AMOVA) of Apium graveolens L. var. dulce showed that approximately 43% of genetic diversity was within accessions, 45% among accessions, and 22% among horticultural types. The neighbor-joining tree generated by unweighted pair group method with arithmetic mean (UPGMA), and population structure analysis, as well as principal components analysis (PCA), separated the cultivars into clusters corresponding to the geographical areas where they originated. Genetic distance analysis suggested that genetic variation within Apium graveolens was quite limited. Genotypic diversity showed any combinations of 55 genic SSRs were able to distinguish the genotypes of all 30 accessions.