SNP discovery and genetic mapping using genotyping by sequencing of whole genome genomic DNA from a pea RIL population

Abstract

BackgroundProgress in genetics and breeding in pea still suffers from the limited availability of molecular resources. SNP markers that can be identified through affordable sequencing processes, without the need for prior genome reduction or a reference genome to assemble sequencing data would allow the discovery and genetic mapping of thousands of molecular markers. Such an approach could significantly speed up genetic studies and marker assisted breeding for non-model species.ResultsA total of 419,024 SNPs were discovered using HiSeq whole genome sequencing of four pea lines, followed by direct identification of SNP markers without assembly using the discoSnp tool. Subsequent filtering led to the identification of 131,850 highly designable SNPs, polymorphic between at least two of the four pea lines.A subset of 64,754 SNPs was called and genotyped by short read sequencing on a subpopulation of 48 RILs from the cross ‘Baccara’ x ‘PI180693’. This data was used to construct a WGGBS-derived pea genetic map comprising 64,263 markers. This map is collinear with previous pea consensus maps and therefore with the Medicago truncatula genome. Sequencing of four additional pea lines showed that 33 % to 64 % of the mapped SNPs, depending on the pairs of lines considered, are polymorphic and can therefore be useful in other crosses.The subsequent genotyping of a subset of 1000 SNPs, chosen for their mapping positions using a KASP™ assay, showed that almost all generated SNPs are highly designable and that most (95 %) deliver highly qualitative genotyping results. Using rather low sequencing coverages in SNP discovery and in SNP inferring did not hinder the identification of hundreds of thousands of high quality SNPs.ConclusionsThe development and optimization of appropriate tools in SNP discovery and genetic mapping have allowed us to make available a massive new genomic resource in pea. It will be useful for both fine mapping within chosen QTL confidence intervals and marker assisted breeding for important traits in pea improvement.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-2447-2) contains supplementary material, which is available to authorized users.

Large-scale SNP discovery and construction of a high-density genetic map of Colossoma macropomum through genotyping-by-sequencing

“…Genotype-by-sequencing (GBS) is one next-generation sequencing technique, and is based on the reduction of genome complexity using restriction enzymes23. GBS is characterized as a simple, quick, specific, reproducible technique23, and has been extensively used to identify a large number of SNPs in various model and non-model species25262728.…”

NGS-Based Genotyping, High-Throughput Phenotyping and Genome-Wide Association Studies Laid the Foundations for Next-Generation Breeding in Horticultural Crops

“…GBS was originally applied in maize to identify SNP markers in a RIL population [65]. Since then, several studies have been reported in literature for large-scale SNP discovery and genotyping of RILs [68,69], bi-parental mapping population segregating for important agronomical traits [70,71]; ethylmethane sulfonate-(EMS) induced mutant populations [72,73]; and unrelated individuals in small or large size populations [74,75]. All of these efforts, most of which were carried out on vegetable crop species, aimed to have available a high number of SNP data points for concomitant or future genome-wide association studies.…”

Transcriptomic Studies in Non-Model Plants: Case of Pisum sativum L. and Medicago lupulina L.

“…RAD-Seq yields fragments distributed randomly over a genome and is suitable for discovering indels (insertion-deletion polymorphisms), SNVs (single nucleotide variations) and microsatellites simple sequence repeats (SSR). Using RAD-Seq, Boutet et al [57] discovered a total of 419,024 SNVs between at least two of the four pea lines analysed in their work. Pea genetic map constructed by genotyping a subset of 64,754 SNVs on a subpopulation of 48 RILs (recombinant inbred lines) was collinear with previous pea consensus maps and therefore with the M. truncatula genome.…”