Assessment of genetic variation within a global collection of lentil (Lens culinaris Medik.) cultivars and landraces using SNP markers

Background

Lentil is a self-pollinated annual diploid (2n = 2× = 14) crop with a restricted history of genetic improvement through breeding, particularly when compared to cereal crops. This limited breeding has probably contributed to the narrow genetic base of local cultivars, and a corresponding potential to continue yield increases and stability. Therefore, knowledge of genetic variation and relationships between populations is important for understanding of available genetic variability and its potential for use in breeding programs. Single nucleotide polymorphism (SNP) markers provide a method for rapid automated genotyping and subsequent data analysis over large numbers of samples, allowing assessment of genetic relationships between genotypes.

Results

In order to investigate levels of genetic diversity within lentil germplasm, 505 cultivars and landraces were genotyped with 384 genome-wide distributed SNP markers, of which 266 (69.2%) obtained successful amplification and detected polymorphisms. Gene diversity and PIC values varied between 0.108-0.5 and 0.102-0.375, with averages of 0.419 and 0.328, respectively. On the basis of clarity and interest to lentil breeders, the genetic structure of the germplasm collection was analysed separately for cultivars and landraces. A neighbour-joining (NJ) dendrogram was constructed for commercial cultivars, in which lentil cultivars were sorted into three major groups (G-I, G-II and G-III). These results were further supported by principal coordinate analysis (PCoA) and STRUCTURE, from which three clear clusters were defined based on differences in geographical location. In the case of landraces, a weak correlation between geographical origin and genetic relationships was observed. The landraces from the Mediterranean region, predominantly Greece and Turkey, revealed very high levels of genetic diversity.

Conclusions

Lentil cultivars revealed clear clustering based on geographical origin, but much more limited correlation between geographic origin and genetic diversity was observed for landraces. These results suggest that selection of divergent parental genotypes for breeding should be made actively on the basis of systematic assessment of genetic distance between genotypes, rather than passively based on geographical distance.

Electronic supplementary material

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

EST-SNP discovery and dense genetic mapping in lentil (Lens culinaris Medik.) enable candidate gene selection for boron tolerance

Large-scale SNP discovery and dense genetic mapping in a lentil intraspecific cross permitted identification of a single chromosomal region controlling tolerance to boron toxicity, an important breeding objective. Lentil (Lens culinaris Medik.) is a highly nutritious food legume crop that is cultivated world-wide. Until recently, lentil has been considered a genomic 'orphan' crop, limiting the feasibility of marker-assisted selection strategies in breeding programs. The present study reports on the identification of single-nucleotide polymorphisms (SNPs) from transcriptome sequencing data, utilisation of expressed sequence tag (EST)-derived simple sequence repeat (SSR) and SNP markers for construction of a gene-based genetic linkage map, and identification of markers in close linkage to major QTLs for tolerance to boron (B) toxicity. A total of 2,956 high-quality SNP markers were identified from a lentil EST database. Sub-sets of 546 SSRs and 768 SNPs were further used for genetic mapping of an intraspecific mapping population (Cassab × ILL2024) that exhibits segregation for B tolerance. Comparative analysis of the lentil linkage map with the sequenced genomes of Medicago truncatula Gaertn., soybean (Glycine max [L.] Merr.) and Lotus japonicus L. indicated blocks of conserved macrosynteny, as well as a number of rearrangements. A single genomic region was found to be associated with variation for B tolerance in lentil, based on evaluation performed over 2 years. Comparison of flanking markers to genome sequences of model species (M. truncatula, soybean and Arabidopsis thaliana) identified candidate genes that are functionally associated with B tolerance, and could potentially be used for diagnostic marker development in lentil.

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.

SNP marker discovery, linkage map construction and identification of QTLs for enhanced salinity tolerance in field pea (Pisum sativum L.)

BACKGROUND

Field pea (Pisum sativum L.) is a self-pollinating, diploid, cool-season food legume. Crop production is constrained by multiple biotic and abiotic stress factors, including salinity, that cause reduced growth and yield. Recent advances in genomics have permitted the development of low-cost high-throughput genotyping systems, allowing the construction of saturated genetic linkage maps for identification of quantitative trait loci (QTLs) associated with traits of interest. Genetic markers in close linkage with the relevant genomic regions may then be implemented in varietal improvement programs.

RESULTS

In this study, single nucleotide polymorphism (SNP) markers associated with expressed sequence tags (ESTs) were developed and used to generate comprehensive linkage maps for field pea. From a set of 36,188 variant nucleotide positions detected through in silico analysis, 768 were selected for genotyping of a recombinant inbred line (RIL) population. A total of 705 SNPs (91.7%) successfully detected segregating polymorphisms. In addition to SNPs, genomic and EST-derived simple sequence repeats (SSRs) were assigned to the genetic map in order to obtain an evenly distributed genome-wide coverage. Sequences associated with the mapped molecular markers were used for comparative genomic analysis with other legume species. Higher levels of conserved synteny were observed with the genomes of Medicago truncatula Gaertn. and chickpea (Cicer arietinum L.) than with soybean (Glycine max [L.] Merr.), Lotus japonicus L. and pigeon pea (Cajanus cajan [L.] Millsp.). Parents and RIL progeny were screened at the seedling growth stage for responses to salinity stress, imposed by addition of NaCl in the watering solution at a concentration of 18 dS m-1. Salinity-induced symptoms showed normal distribution, and the severity of the symptoms increased over time. QTLs for salinity tolerance were identified on linkage groups Ps III and VII, with flanking SNP markers suitable for selection of resistant cultivars. Comparison of sequences underpinning these SNP markers to the M. truncatula genome defined genomic regions containing candidate genes associated with saline stress tolerance.

CONCLUSION

The SNP assays and associated genetic linkage maps developed in this study permitted identification of salinity tolerance QTLs and candidate genes. This constitutes an important set of tools for marker-assisted selection (MAS) programs aimed at performance enhancement of field pea cultivars.

Improving yield potential in crops under elevated CO(2): Integrating the photosynthetic and nitrogen utilization efficiencies

Increasing crop productivity to meet burgeoning human food demand is challenging under changing environmental conditions. Since industrial revolution atmospheric CO(2) levels have linearly increased. Developing crop varieties with increased utilization of CO(2) for photosynthesis is an urgent requirement to cope with the irreversible rise of atmospheric CO(2) and achieve higher food production. The primary effects of elevated CO(2) levels in most crop plants, particularly C(3) plants, include increased biomass accumulation, although initial stimulation of net photosynthesis rate is only temporal and plants fail to sustain the maximal stimulation, a phenomenon known as photosynthesis acclimation. Despite this acclimation, grain yield is known to marginally increase under elevated CO(2). The yield potential of C(3) crops is limited by their capacity to exploit sufficient carbon. The "C fertilization" through elevated CO(2) levels could potentially be used for substantial yield increase. Rubisco is the rate-limiting enzyme in photosynthesis and its activity is largely affected by atmospheric CO(2) and nitrogen availability. In addition, maintenance of the C/N ratio is pivotal for various growth and development processes in plants governing yield and seed quality. For maximizing the benefits of elevated CO(2), raising plant nitrogen pools will be necessary as part of maintaining an optimal C/N balance. In this review, we discuss potential causes for the stagnation in yield increases under elevated CO(2) levels and explore possibilities to overcome this limitation by improved photosynthetic capacity and enhanced nitrogen use efficiency. Opportunities of engineering nitrogen uptake, assimilatory, and responsive genes are also discussed that could ensure optimal nitrogen allocation toward expanding source and sink tissues. This might avert photosynthetic acclimation partially or completely and drive for improved crop production under elevated CO(2) levels.

Transcriptome sequencing of field pea and faba bean for discovery and validation of SSR genetic markers

BACKGROUND

Field pea (Pisum sativum L.) and faba bean (Vicia faba L.) are cool-season grain legume species that provide rich sources of food for humans and fodder for livestock. To date, both species have been relative 'genomic orphans' due to limited availability of genetic and genomic information. A significant enrichment of genomic resources is consequently required in order to understand the genetic architecture of important agronomic traits, and to support germplasm enhancement, genetic diversity, population structure and demographic studies.

RESULTS

cDNA samples obtained from various tissue types of specific field pea and faba bean genotypes were sequenced using 454 Roche GS FLX Titanium technology. A total of 720,324 and 304,680 reads for field pea and faba bean, respectively, were de novo assembled to generate sets of 70,682 and 60,440 unigenes. Consensus sequences were compared against the genome of the model legume species Medicago truncatula Gaertn., as well as that of the more distantly related, but better-characterised genome of Arabidopsis thaliana L.. In comparison to M. truncatula coding sequences, 11,737 and 10,179 unique hits were obtained from field pea and faba bean. Totals of 22,057 field pea and 18,052 faba bean unigenes were subsequently annotated from GenBank. Comparison to the genome of soybean (Glycine max L.) resulted in 19,451 unique hits for field pea and 16,497 unique hits for faba bean, corresponding to c. 35% and 30% of the known gene space, respectively. Simple sequence repeat (SSR)-containing expressed sequence tags (ESTs) were identified from consensus sequences, and totals of 2,397 and 802 primer pairs were designed for field pea and faba bean. Subsets of 96 EST-SSR markers were screened for validation across modest panels of field pea and faba bean cultivars, as well as related non-domesticated species. For field pea, 86 primer pairs successfully obtained amplification products from one or more template genotypes, of which 59% revealed polymorphism between 6 genotypes. In the case of faba bean, 81 primer pairs displayed successful amplification, of which 48% detected polymorphism.

CONCLUSIONS

The generation of EST datasets for field pea and faba bean has permitted effective unigene identification and functional sequence annotation. EST-SSR loci were detected at incidences of 14-17%, permitting design of comprehensive sets of primer pairs. The subsets from these primer pairs proved highly useful for polymorphism detection within Pisum and Vicia germplasm.

Transcriptome sequencing of lentil based on second-generation technology permits large-scale unigene assembly and SSR marker discovery

Background

Lentil (Lens culinaris Medik.) is a cool-season grain legume which provides a rich source of protein for human consumption. In terms of genomic resources, lentil is relatively underdeveloped, in comparison to other Fabaceae species, with limited available data. There is hence a significant need to enhance such resources in order to identify novel genes and alleles for molecular breeding to increase crop productivity and quality.

Results

Tissue-specific cDNA samples from six distinct lentil genotypes were sequenced using Roche 454 GS-FLX Titanium technology, generating c. 1.38 × 10⁶ expressed sequence tags (ESTs). De novo assembly generated a total of 15,354 contigs and 68,715 singletons. The complete unigene set was sequence-analysed against genome drafts of the model legume species Medicago truncatula and Arabidopsis thaliana to identify 12,639, and 7,476 unique matches, respectively. When compared to the genome of Glycine max, a total of 20,419 unique hits were observed corresponding to c. 31% of the known gene space. A total of 25,592 lentil unigenes were subsequently annoated from GenBank. Simple sequence repeat (SSR)-containing ESTs were identified from consensus sequences and a total of 2,393 primer pairs were designed. A subset of 192 EST-SSR markers was screened for validation across a panel 12 cultivated lentil genotypes and one wild relative species. A total of 166 primer pairs obtained successful amplification, of which 47.5% detected genetic polymorphism.

Conclusions

A substantial collection of ESTs has been developed from sequence analysis of lentil genotypes using second-generation technology, permitting unigene definition across a broad range of functional categories. As well as providing resources for functional genomics studies, the unigene set has permitted significant enhancement of the number of publicly-available molecular genetic markers as tools for improvement of this species.

Genotyping elite genotypes within the Australian lentil breeding program with lentil-specific sequenced tagged microsatellite site (STMS) markers

Lentil (Lens culinaris ssp. culinaris) is consumed in many countries as a rich source of protein in largely vegetarian diets. Australia grows lentil as a cash crop in rotation with cereal and produces predominantly red lentils that are exported throughout the world, particularly to countries in South Asia and the Middle East. Differentiation of varieties is important when exporting products to such markets, maintaining variety purity during seed production and in the collection of end-point royalties. Lentil-specific and fluorescent sequenced tagged microsatellite markers (STMS) markers were used to construct a DNA fingerprint database for 10 Lens culinaris ssp. culinaris genotypes (Northfield, Digger, ILL7537, Nugget, Indianhead, ILL2024, ILL6788, Palouse, Nipper and Boomer) that represent major new cultivars and key breeding lines within the Australian breeding program. All 10 lentil genotypes were distinguished using the assessed STMS loci. Unique alleles were observed for several lines, including Boomer and Nipper, varieties recently released in Australia. This database will play an important role in seed typing for commercial export certification and the commercial management of cultivars.

Molecular genetic diversity and linked resistance to ascochyta blight in Australian chickpea breeding materials and their wild relatives

Simple sequence-repeat (SSR) and sequence characterised amplified regions (SCARs) have been used to characterise the genetic diversity of chickpea germplasm. A set of 48 genotypes comprising cultigen, landraces, and wild relatives important for breeding purposes was used to determine the genetic similarity between genotypes and to assess the association between ascochyta blight (AB) and SCAR phenotypes. The 21 SSR markers amplified a total of 370 alleles, with an average of ~17 alleles per SSR locus among the 48 genotypes. Polymorphic information content (PIC) values ranged from 0.37 for the XGA13 locus to 0.93 for the XGA106. Principal coordinate analysis (PCO) of genetic similarity (GS) estimates revealed a clear differentiation of the chickpea genotypes into 5 groups, which were generally consistent with available pedigree information. Comparison of SCAR and AB phenotypes enabled us to tag the common source(s) of AB resistance in the breeding collection. Based on the SCAR phenotypes, it was evident that the studied chickpea genotypes, including worldwide-known AB-resistant lines (ICC12004, ILC72, ILC3279), carry at least one common source of resistance to AB. Since SSR markers are polymerase chain reaction (PCR)-based markers, highly polymorphic, and amenable to high-throughput technologies, they are therefore well suited for studies of genetic diversity and cultivar identification in chickpea. The broad level of genetic diversity detected in the chickpea germplasm should be useful for selective breeding for specific traits such as AB, backcrossing, and in enhancing the genetic base of breeding programs.

Extensive macrosynteny between Medicago truncatula and Lens culinaris ssp. culinaris

The first predominantly gene-based genetic linkage map of lentil (Lens culinaris ssp. culinaris) was constructed using an F5 population developed from a cross between the cultivars Digger (ILL5722) and Northfield (ILL5588) using 79 intron-targeted amplified polymorphic (ITAP) and 18 genomic simple sequence repeat (SSR) markers. Linkage analysis revealed seven linkage groups (LGs) comprised of 5-25 markers that varied in length from 80.2 to 274.6 cM. The genome map spanned a total length of 928.4 cM. Clear evidence of a simple and direct macrosyntenic relationship between lentil and Medicago truncatula was observed. Sixty-six out of the 71 gene-based markers, which were previously assigned to M. truncatula genetic and physical maps, were found in regions syntenic between the Lens c. ssp. culinaris and M. truncatula genomes. However, there was evidence of moderate chromosomal rearrangements which may account for the difference in chromosome numbers between these two legume species. Eighteen common SSR markers were used to connect the current map with the most comprehensive and recent map that exists for lentil, providing the syntenic context of four important domestication traits. The composite map presented, anchored with orthologous markers mapped in M. truncatula, provides a strong foundation for the future use of genomic and genetic information in lentil genetic analysis and breeding.