Development of 65 novel polymorphic cDNA-SSR markers in common vetch (Vicia sativa subsp. sativa) using next generation sequencing

Designing Novel Strategies for Improving Old Legumes: An Overview from Common Vetch

Common vetch (Vicia sativa L.) is a grain legume used in animal feeding, rich in protein content, fatty acid, and mineral composition that makes for a very adequate component to enrich feedstuff. In addition, relevant pharmacological properties have been reported in humans. The common vetch, similar to other legumes, can fix atmospheric nitrogen, a crucial feature for sustainable agricultural systems. These properties enhance the use of vetch as a cover crop and its sowing in intercropping systems. Moreover, several studies have recently pointed out the potential of vetch in the phytoremediation of contaminated soils. These characteristics make vetch a relevant crop, which different potential improvements target. Varieties with different yields, flowering times, shattering resistance, nutritional composition, rhizobacteria associations, drought tolerance, nitrogen fixation capacity, and other agronomic-relevant traits have been identified when different vetch accessions are compared. Recently, the analysis of genomic and transcriptomic data has allowed the development of different molecular markers to be used for assisted breeding purposes, promoting crop improvement. Here, we review the potential of using the variability of V. sativa genetic resources and new biotechnological and molecular tools for selecting varieties with improved traits to be used in sustainable agriculture systems.

Screening and identification of multiple abiotic stress responsive candidate genes based on hybrid-sequencing in Vicia sativa

Common vetch is an important leguminous forage for both livestock fodder and green manure and has a tremendous latent capacity in a sustainable agroecosystem. In the present study, a comprehensive transcriptome analysis of the aboveground leaves and underground roots of common vetch under multiple abiotic stress treatments, including NaCl, drought, cold, and cold drought, was performed using hybrid-sequencing technology, i. e. single-molecule real-time sequencing technology (SMRT) and supplemented by next-generation sequencing (NGS) technology. A total of 485,038 reads of insert (ROIs) with a mean length of 2606 bp and 228,261 full-length nonchimeric (FLNC) reads were generated. After deduplication, 39,709 transcripts were generated. Of these transcripts, we identified 1059 alternative splicing (AS) events, 17,227 simple sequence repeats (SSRs), and 1647 putative transcription factors (TFs). Furthermore, 640 candidates long noncoding RNAs (lncRNAs) and 28,256 complete coding sequences (CDSs) were identified. In gene annotation analyses, a total of 38,826 transcripts (97.78%) were annotated in eight public databases. Finally, seven multiple abiotic stress-responsive candidate genes were obtained through gene expression, annotation information, and protein-protein interaction (PPI) networks. Our research not only enriched the structural information of FL transcripts in common vetch, but also provided useful information for exploring the molecular mechanism of multiple abiotic stress tolerance between aboveground and underground tissues in common vetch and related legumes.

Genome-Wide Identification and Development of LTR Retrotransposon-Based Molecular Markers for the Melilotus Genus

Melilotus is an important genus of legumes with industrial and medicinal value, partly due to the production of coumarin. To explore the genetic diversity and population structure of Melilotus, 40 accessions were analyzed using long terminal repeat (LTR) retrotransposon-based markers. A total of 585,894,349 bp of LTR retrotransposon sequences, accounting for 55.28% of the Melilotus genome, were identified using bioinformatics tools. A total of 181,040 LTR retrotransposons were identified and classified as Gypsy, Copia, or another type. A total of 350 pairs of primers were designed for assessing polymorphisms in 15 Melilotus albus accessions. Overall, 47 polymorphic primer pairs were screened for their availability and transferability in 18 Melilotus species. All the primer pairs were transferable, and 292 alleles were detected at 47 LTR retrotransposon loci. The average polymorphism information content (PIC) value was 0.66, which indicated that these markers were highly informative. Based on unweighted pair group method with arithmetic mean (UPGMA) dendrogram cluster analysis, the 18 Melilotus species were classified into three clusters. This study provides important data for future breeding programs and for implementing genetic improvements in the Melilotus genus.

Common Vetch, Valuable Germplasm for Resilient Agriculture: Genetic Characterization and Spanish Core Collection Development

Common vetch (Vicia sativa L.) is a legume used for animal feed because of its high protein content and great capacity for nitrogen fixation, making this crop relevant in sustainable agriculture. The Spanish vetch collection, conserved at the Spanish Plant Genetic Resources Center (CRF), is one of the largest collections of this species worldwide, including landraces, wild relatives mainly collected in Spain, and commercial cultivars, but also accessions of international origin. The analysis of the genetic diversity of this material, whose genome has not been sequenced yet, and the assembly of a representative collection could play a pivotal role in conserving and exploiting these genetic resources in breeding programs mainly in those focused on consequences and demands of climate change. In this work, a set of 14 simple sequence repeat (SSR) reference alleles for genetic diversity analysis of the CRF vetch collection has been developed, used for genotyping more than 545 common vetch accessions from all over the world and validated. All the tested markers were polymorphic for the analyzed accessions. Overall, at least 86 different loci were identified with 2-11 alleles per locus with an average of 6.1 alleles per locus. Also, the analyses of the generated SSR database support that most of these SSR markers are transferable across closely related species of Vicia genus. Analysis of molecular variance revealed that wild relatives have a higher genetic diversity than landraces. However, cultivars have similar diversity than landraces, indicating that genetic variability has been barely lost due to the breeding of this legume. Low differences of genetic variations between Spanish and non-Spanish accessions have been observed, suggesting a high degree of diversity within Spanish genotypes, which provide 95% of the total genetic variation, so we have focused our efforts on characterizing genotypes of Spanish origin that were further studied using storage protein profiles. Based on SSR, seed protein profiles, and agromorphological and passport data, a vetch core collection (VCC) containing 47 V. sativa accessions of Spanish origin has been established. In this collection, the characterization has been expanded using ISSR markers, and it has been reevaluated with new agromorphological data, including drought tolerance characters. This VCC presents a minimum loss of genetic diversity concerning the total collection and constitutes an invaluable material that can be used in future breeding programs for direct use in a resilient agricultural system.

Diversity of Naturalized Hairy Vetch (Vicia villosa Roth) Populations in Central Argentina as a Source of Potential Adaptive Traits for Breeding

Hairy vetch (Vicia villosa ssp. villosa Roth) is native of Europe and Western Asia and it is the second most cultivated vetch worldwide. Hairy vetch is used as forage species in semiarid environments and as a legume cover crop in sub-humid and humid regions. Being an incompletely domesticated species, hairy vetch can form spontaneous populations in a new environment. These populations might contain novel and adaptive traits valuable for breeding. Niche occupancy based on geographic occurrence and environmental data of naturalized populations in central Argentina showed that these populations were distributed mainly on disturbed areas with coarse soil texture and alkaline-type soils. Low rainfall and warm temperatures during pre- and post-seed dispersal explained the potential distribution under sub-humid and semiarid conditions from Pampa and Espinal ecoregions. Conversely, local adaptation along environmental gradients did not drive the divergence among recently established Argentinian (AR) populations. The highest genetic diversity revealed by microsatellite analysis was observed within accessions (72%) while no clear separation was detected between AR and European (EU) genotypes, although naturalized AR populations showed strong differentiation with the wild EU accessions. Common garden experiments were conducted in 2014-16 in order to evaluate populations' germination, flowering, and biomass traits. European cultivars were characterized by low physical seed dormancy (PY), while naturalized AR accessions showed higher winter biomass production. Detected variation in the quantitative assessment of populations could be useful for selection in breeding for traits that convey favorable functions within specific contexts.

De Novo Assembly and Discovery of Genes That Involved in Drought Tolerance in the Common Vetch

The common vetch (Vicia sativa) is often used as feed for livestock because of its high nutritional value. However, drought stress reduces forage production through plant damage. Here, we studied the transcriptional profiles of common vetch exposed to drought in order to understand the molecular mechanisms of drought tolerance in this species. The genome of the common vetch has not been sequenced, therefore we used Illumina sequencing to generate de novo transcriptomes. Nearly 500 million clean reads were used to generate 174,636 transcripts, including 122,299 unigenes. In addition, 5313 transcription factors were identified and these transcription factors were classified into 79 different gene families. We also identified 11,181 SSR loci from di- to hexa-nucleotides whose repeat number was greater than five. On the basis of differentially expressed genes, Gene Ontology analysis identified many drought-relevant categories, including “oxidation-reduction process”, “lipid metabolic process” and “oxidoreductase activity”. In addition to these, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis identified pathways, such as “Plant hormone signal transduction”, “Glycolysis/Gluconeogenesis” and “Phenylpropanoid biosynthesis”, as differentially expressed in the plants exposed to drought. The expression results in this study will be useful for further extending our knowledge on the drought tolerance of common vetch.

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.

Cross-species transferability of EST-SSR markers developed from the transcriptome of Melilotus and their application to population genetics research

Melilotus is one of the most important legume forages, but the lack of molecular markers has limited the development and utilization of Melilotus germplasm resources. In the present study, 151 M clean reads were generated from various genotypes of Melilotus albus using Illumina sequencing. A total of 19,263 potential EST-SSRs were identified from 104,358 unigene sequences. Moreover, 18,182 primer pairs were successfully designed, and 550 primer pairs were selected using criteria of base repeat type, fragment length and annealing temperature. In addition, 550 primer pairs were screened by using PCR amplification products and used to assess polymorphisms in 15 M. albus accessions. A total of 114 primer pairs were detected as being highly polymorphic, and the average polymorphism information content (PIC) value was 0.79. Furthermore, those 114 polymorphic primer pairs were used to evaluate the transferability to 18 species of the genus Melilotus, and 70 EST-SSR markers were found to be transferable among the 18 Melilotus species. According to the UPGMA dendrogram and STRUCTURE analysis, the 18 Melilotus species were classified into three clusters. This study offers a valuable resource for the genetic diversity and molecular assisted breeding of germplasm resources in the genus Melilotus.

Determination of neurotoxic agents as markers of common vetch adulteration in lentil by LC-MS/MS

Vicia sativa (common vetch), which contains β-cyanoalanine (BCA) and γ-glutamyl β-cyanoalanine (γGBCA) is used for adulteration of lentil and have neuro toxic effects on people and animals. Therefore the determination of BCA and γGBCA analytes in accurate and reliable manner has high importance for human and animal health. LC-IDMS/MS method has been developed and validated for quantitative analyses of β-cyanoalanine and ¹³C labeled BCA is used as an internal standard in γGBCA analysis. In order to validate the method, linearity, recovery, precision (repeatability), intermediate precision, limit of detection, and limit of quantification parameters were investigated. The correlation coefficient was found to be greater than 0.99 for both analytes. The recoveriess were determined as 95.8% for BCA and 97.4% for γGBCA. The relative expanded measurement uncertainties of β-cyanoalanine and γ-glutamyl-β-cyanoalanine were obtained as 4.6% and 5.8%, respectively with the coverage factor, k, is 2 and at 95% confidence level.

Optimizing Sample Size to Assess the Genetic Diversity in Common Vetch (Vicia sativa L.) Populations Using Start Codon Targeted (SCoT) Markers

Common vetch (Vicia sativa subsp. sativa L.) is a self-pollinating annual forage legume with worldwide importance. Here, we investigate the optimal number of individuals that may represent the genetic diversity of a single population, using Start Codon Targeted (SCoT) markers. Two cultivated varieties and two wild accessions were evaluated using five SCoT primers, also testing different sampling sizes: 1, 2, 3, 5, 8, 10, 20, 30, 40, 50, and 60 individuals. The results showed that the number of alleles and the Polymorphism Information Content (PIC) were different among the four accessions. Cluster analysis by Unweighted Pair Group Method with Arithmetic Mean (UPGMA) and STRUCTURE placed the 240 individuals into four distinct clusters. The Expected Heterozygosity (H_E) and PIC increased along with an increase in sampling size from 1 to 10 plants but did not change significantly when the sample sizes exceeded 10 individuals. At least 90% of the genetic variation in the four germplasms was represented when the sample size was 10. Finally, we concluded that 10 individuals could effectively represent the genetic diversity of one vetch population based on the SCoT markers. This study provides theoretical support for genetic diversity, cultivar identification, evolution, and marker-assisted selection breeding in common vetch.

Transcriptome Analyses Reveal Candidate Pod Shattering-Associated Genes Involved in the Pod Ventral Sutures of Common Vetch (Vicia sativa L.)

The seed dispersion caused by pod shattering is a form of propagation used by many wild species. Loss of seeds from pod shattering is frequent in the common vetch (Vicia sativa L.), an important self-pollinating annual forage legume. However, pod shattering is one of the most important defects that limits the reproduction of the vetch in the field and the usage as a leguminous forage crop. To better understand the vetch pod shattering mechanism, we used high-throughput RNA sequencing to assess the global changes in the transcriptomes of the pod ventral sutures of shattering-susceptible and shattering-resistant vetch accessions screened from 541 vetch germplasms. A total of 1,285 significantly differentially expressed unigenes (DEGs) were detected, including 575 up-regulated unigenes and 710 down-regulated unigenes. Analyses of Gene Ontology and KEGG metabolic enrichment pathways of 1,285 DEGs indicated that 22 DEGs encoding cell wall modifications and hydrolases associated with pod shattering were highly expressed in shattering-susceptible accessions. These genes were mainly enriched in "hydrolase activity," "cytoplasm," and "carbohydrate metabolic process" systems. These cell wall modifications and hydrolases genes included β-glucosidase and endo-polygalacturonase, which work together to break down the glycosidic bonds of pectin and cellulose, and to promote the dissolution and disappearance of the cell wall in the ventral suture of the pod and make the pod more susceptible to shattering. We demonstrated the differences in gene transcription levels between the shattering-susceptible and shattering-resistant vetch accessions for the first time and our results provided valuable information for the identifying and characterizing of pod shattering regulation networks in vetch. This information may facilitate the future identification of pod shattering-related genes and their underlying molecular mechanisms in the common vetch.

Development and cross-species transferability of EST-SSR markers in Siberian wildrye (Elymus sibiricus L.) using Illumina sequencing

Siberian wildrye (Elymus sibiricus L.) is a perennial, self-fertilizing grass that plays an important role in animal husbandry and environmental sustenance. However, the transcriptomic and genomic information on this species is very limited, which hinders genetic and breeding studies. In the present study, 76,686,804 clean reads were generated from 11 different tissue samples of E. sibiricus by Illumina paired-end sequencing, and the reads were deposited into the NCBI SRA database (SRX574376). A total of 8,769 EST-SSRs were identified from 94,458 unigene sequences, which were obtained by de novo assembly. Moreover, 1,078 primer pairs were successfully designed, and 500 pairs were randomly selected to assess polymorphisms in 15 E. sibiricus accessions. A total of 112 primer pairs were polymorphic, and the polymorphism information content (PIC) values ranged from 0.39 to 0.81, indicating a high level of informativeness. Furthermore, these 112 polymorphic primer pairs were used to evaluate the transferability to 13 other related species, and 55 EST-SSR markers were found to be polymorphic among these 13 Elymus species. This study collected the global sequence data for E. sibiricus, and the newly developed markers will prove valuable in facilitating genetic diversity in E. sibiricus and related Elymus species.

Cross-Species, Amplifiable EST-SSR Markers for Amentotaxus Species Obtained by Next-Generation Sequencing

Amentotaxus, a genus of Taxaceae, is an ancient lineage with six relic and endangered species. Four Amentotaxus species, namely A. argotaenia, A. formosana, A. yunnanensis, and A. poilanei, are considered a species complex because of their morphological similarities. Small populations of these species are allopatrically distributed in Asian forests. However, only a few codominant markers have been developed and applied to study population genetic structure of these endangered species. In this study, we developed and characterized polymorphic expressed sequence tag-simple sequence repeats (EST-SSRs) from the transcriptome of A. formosana. We identified 4955 putative EST-SSRs from 68,281 unigenes as potential molecular markers. Twenty-six EST-SSRs were selected for estimating polymorphism and transferability among Amentotaxus species, of which 23 EST-SSRs were polymorphic within Amentotaxus species. Among these, the number of alleles ranged from 1-4, the polymorphism information content ranged from 0.000-0.692, and the observed and expected heterozygosity were 0.000-1.000 and 0.080-0.740, respectively. Population genetic structure analyses confirmed that A. argotaenia and A. formosana were separate species and A. yunnanensis and A. poilanei were the same species. These novel EST-SSRs can facilitate further population genetic structure research of Amentotaxus species.

Transcriptome Analysis of Two Vicia sativa Subspecies: Mining Molecular Markers to Enhance Genomic Resources for Vetch Improvement

The vetch (Vicia sativa) is one of the most important annual forage legumes globally due to its multiple uses and high nutritional content. Despite these agronomical benefits, many drawbacks, including cyano-alanine toxin, has reduced the agronomic value of vetch varieties. Here, we used 454 technology to sequence the two V. sativa subspecies (ssp. sativa and ssp. nigra) to enrich functional information and genetic marker resources for the vetch research community. A total of 86,532 and 47,103 reads produced 35,202 and 18,808 unigenes with average lengths of 735 and 601 bp for V. sativa sativa and V. sativa nigra, respectively. Gene Ontology annotations and the cluster of orthologous gene classes were used to annotate the function of the Vicia transcriptomes. The Vicia transcriptome sequences were then mined for simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers. About 13% and 3% of the Vicia unigenes contained the putative SSR and SNP sequences, respectively. Among those SSRs, 100 were chosen for the validation and the polymorphism test using the Vicia germplasm set. Thus, our approach takes advantage of the utility of transcriptomic data to expedite a vetch breeding program.

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.

Characterization and development of EST-SSR markers derived from transcriptome of yellow catfish

Yellow catfish (Pelteobagrus fulvidraco) is one of the most important freshwater fish due to its delicious flesh and high nutritional value. However, lack of sufficient simple sequence repeat (SSR) markers has hampered the progress of genetic selection breeding and molecular research for yellow catfish. To this end, we aimed to develop and characterize polymorphic expressed sequence tag (EST)–SSRs from the 454 pyrosequencing transcriptome of yellow catfish. Totally, 82,794 potential EST-SSR markers were identified and distributed in the coding and non-coding regions. Di-nucleotide (53,933) is the most abundant motif type, and AC/GT, AAT/ATT, AAAT/ATTT are respective the most frequent di-, tri-, tetra-nucleotide repeats. We designed primer pairs for all of the identified EST-SSRs and randomly selected 300 of these pairs for further validation. Finally, 263 primer pairs were successfully amplified and 57 primer pairs were found to be consistently polymorphic when four populations of 48 individuals were tested. The number of alleles for the 57 loci ranged from 2 to 17, with an average of 8.23. The observed heterozygosity (H_O), expected heterozygosity (H_E), polymorphism information content (PIC) and fixation index (FIS) values ranged from 0.04 to 1.00, 0.12 to 0.92, 0.12 to 0.91 and −0.83 to 0.93, respectively. These EST-SSR markers generated in this study could greatly facilitate future studies of genetic diversity and molecular breeding in yellow catfish.

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.

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.

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.