The Wheat 660K SNP array demonstrates great potential for marker‐assisted selection in polyploid wheat

The rapid development and application of molecular marker assays have facilitated genomic selection and genome‐wide linkage and association studies in wheat breeding. Although PCR‐based markers (e.g. simple sequence repeats and functional markers) and genotyping by sequencing have contributed greatly to gene discovery and marker‐assisted selection, the release of a more accurate and complete bread wheat reference genome has resulted in the design of single‐nucleotide polymorphism (SNP) arrays based on different densities or application targets. Here, we evaluated seven types of wheat SNP arrays in terms of their SNP number, distribution, density, associated genes, heterozygosity and application. The results suggested that the Wheat 660K SNP array contained the highest percentage (99.05%) of genome‐specific SNPs with reliable physical positions. SNP density analysis indicated that the SNPs were almost evenly distributed across the whole genome. In addition, 229 266 SNPs in the Wheat 660K SNP array were located in 66 834 annotated gene or promoter intervals. The annotated genes revealed by the Wheat 660K SNP array almost covered all genes revealed by the Wheat 35K (97.44%), 55K (99.73%), 90K (86.9%) and 820K (85.3%) SNP arrays. Therefore, the Wheat 660K SNP array could act as a substitute for other 6 arrays and shows promise for a wide range of possible applications. In summary, the Wheat 660K SNP array is reliable and cost‐effective and may be the best choice for targeted genotyping and marker‐assisted selection in wheat genetic improvement.     

Using Next Generation Sequencing for Multiplexed Trait-Linked Markers in Wheat

With the advent of next generation sequencing (NGS) technologies, single nucleotide polymorphisms (SNPs) have become the major type of marker for genotyping in many crops. However, the availability of SNP markers for important traits of bread wheat ( Triticum aestivum L.) that can be effectively used in marker-assisted selection (MAS) is still limited and SNP assays for MAS are usually uniplex. A shift from uniplex to multiplex assays will allow the simultaneous analysis of multiple markers and increase MAS efficiency. We designed 33 locus-specific markers from SNP or indel-based marker sequences that linked to 20 different quantitative trait loci (QTL) or genes of agronomic importance in wheat and analyzed the amplicon sequences using an Ion Torrent Proton Sequencer and a custom allele detection pipeline to determine the genotypes of 24 selected germplasm accessions. Among the 33 markers, 27 were successfully multiplexed and 23 had 100% SNP call rates. Results from analysis of "kompetitive allele-specific PCR" (KASP) and sequence tagged site (STS) markers developed from the same loci fully verified the genotype calls of 23 markers. The NGS-based multiplexed assay developed in this study is suitable for rapid and high-throughput screening of SNPs and some indel-based markers in wheat.     

Single nucleotide polymorphism genotyping in polyploid wheat with the Illumina GoldenGate assay

Single nucleotide polymorphisms (SNPs) are indispensable in such applications as association mapping and construction of high-density genetic maps. These applications usually require genotyping of thousands of SNPs in a large number of individuals. Although a number of SNP genotyping assays are available, most of them are designed for SNP genotyping in diploid individuals. Here, we demonstrate that the Illumina GoldenGate assay could be used for SNP genotyping of homozygous tetraploid and hexaploid wheat lines. Genotyping reactions could be carried out directly on genomic DNA without the necessity of preliminary PCR amplification. A total of 53 tetraploid and 38 hexaploid homozygous wheat lines were genotyped at 96 SNP loci. The genotyping error rate estimated after removal of low-quality data was 0 and 1% for tetraploid and hexaploid wheat, respectively. Developed SNP genotyping assays were shown to be useful for genotyping wheat cultivars. This study demonstrated that the GoldenGate assay is a very efficient tool for high-throughput genotyping of polyploid wheat, opening new possibilities for the analysis of genetic variation in wheat and dissection of genetic basis of complex traits using association mapping approach. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Rapid identification of an adult plant stripe rust resistance gene in hexaploid wheat by high-throughput SNP array genotyping of pooled extremes

Key message

High-throughput SNP array analysis of pooled extreme phenotypes in a segregating population by KASP marker genotyping permitted rapid, cost-effective location of a stripe rust resistance QTL in wheat.

Abstract

German wheat cultivar “Friedrichswerther” has exhibited high levels of adult plant resistance (APR) to stripe rust in field environments for many years. F_2:3 lines and F₆ recombinant inbred line (RILs) populations derived from a cross between Friedrichswerther and susceptible landrace Mingxian 169 were evaluated in the field in 2013, 2016 and 2017. Illumina 90K iSelect SNP arrays were used to genotype bulked extreme pools and parents; 286 of 1135 polymorphic SNPs were identified on chromosome 6B. Kompetitive Allele-Specific PCR (KASP) markers were used to verify the chromosome region associated with the resistance locus. A linkage map was constructed with 18 KASP-SNP markers, and a major effect QTL was identified within a 1.4 cM interval flanked by KASP markers IWB71602 and IWB55937 in the region 6BL3-0-0.36. The QTL, named QYr.nwafu-6BL, was stable across environments, and explained average 54.4 and 47.8% of the total phenotypic variation in F_2:3 lines and F₆ RILs, respectively. On the basis of marker genotypes, pedigree analysis and relative genetic distance QYr.nwafu-6BL is likely to be a new APR QTL. Combined high-throughput SNP array genotyping of pooled extremes and validation by KASP assays lowers sequencing costs compared to genome-wide association studies with SNP arrays, and more importantly, permits rapid isolation of major effect QTL in hexaploid wheat as well as improving accuracy of mapping in the QTL region. QYr.nwafu-6BL with flanking KASP markers developed and verified in a subset of 236 diverse lines can be used in marker-assisted selection to improve stripe rust resistance in breeding programs.

     

Accelerating public sector rice breeding with high-density KASP markers derived from whole genome sequencing of <Emphasis Type="Italic">indica</Emphasis> rice

Few public sector rice breeders have the capacity to use NGS-derived markers in their breeding programmes despite rapidly expanding repositories of rice genome sequence data. They rely on >?18,000 mapped microsatellites (SSRs) for marker-assisted selection (MAS) using gel analysis. Lack of knowledge about target SNP and InDel variant loci has hampered the uptake by many breeders of Kompetitive allele-specific PCR (KASP), a proprietary technology of LGC genomics that can distinguish alleles at variant loci. KASP is a cost-effective single-step genotyping technology, cheaper than SSRs and more flexible than genotyping by sequencing (GBS) or array-based genotyping when used in selection programmes. Before this study, there were 2015 rice KASP marker loci in the public domain, mainly identified by array-based screening, leaving large proportions of the rice genome with no KASP coverage. Here we have addressed the urgent need for a wide choice of appropriate rice KASP assays and demonstrated that NGS can detect many more KASP to give full genome coverage. Through re-sequencing of nine indica rice breeding lines or released varieties, this study has identified 2.5 million variant sites. Stringent filtering of variants generated 1.3 million potential KASP assay designs, including 92,500 potential functional markers. This strategy delivers a 650-fold increase in potential selectable KASP markers at a density of 3.1 per 1 kb in the indica crosses analysed and 377,178 polymorphic KASP design sites on average per cross. This knowledge is available to breeders and has been utilised to improve the efficiency of public sector breeding in Nepal, enabling identification of polymorphic KASP at any region or quantitative trait loci in relevant crosses. Validation of 39 new KASP was carried out by genotyping progeny from a range of crosses to show that they detected segregating alleles. The new KASP have replaced SSRs to aid trait selection during marker-assisted backcrossing in these crosses, where target traits include rice blast and BLB resistance loci. Furthermore, we provide the software for plant breeders to generate KASP designs from their own datasets.     

水稻单核苷酸多态性及其应用现状

单核苷酸多态性（single nucleotide polymorphisms, SNPs）在水稻中数量多,分布密度高,遗传稳定性高。水稻SNPs的发现方法主要有对样本DNA的PCR产物直接测序、从SSR区段检测SNPs和从基因组序列直接搜索等。目前已有多种基因分型技术运用到了水稻SNPs检测,SNPs检测的高度自动化使水稻SNPs基因分型非常方便。单核苷酸多态性在水稻遗传图谱的构建、基因克隆和功能基因组学研究、标记辅助选择育种、遗传资源分类及物种进化等方面的应用具有巨大潜力。     

Introgression lines of Solanum sitiens,a wild nightshade of the Atacama Desert,in the genome of cultivated tomato

The wild tomato relative Solanum sitiens is a xerophyte endemic to the Atacama Desert of Chile and a potential source of genes for tolerance to drought, salinity and low‐temperature stresses. However, until recently, strong breeding barriers prevented its hybridization and introgression with cultivated tomato, Solanum lycopersicum L. We overcame these barriers using embryo rescue, bridging lines and allopolyploid hybrids, and synthesized a library of introgression lines (ILs) that captures the genome of S. sitiens in the background of cultivated tomato. The IL library consists of 56 overlapping introgressions that together represent about 93% of the S. sitiens genome: 65% in homozygous and 28% in heterozygous (segregating) ILs. The breakpoints of each segment and the gaps in genome coverage were mapped by single nucleotide polymorphism (SNP) genotyping using the SolCAP SNP array. Marker‐assisted selection was used to backcross selected introgressions into tomato, to recover a uniform genetic background, to isolate recombinant sub‐lines with shorter introgressions and to select homozygous genotypes. Each IL contains a single S. sitiens chromosome segment, defined by markers, in the genetic background of cv. NC 84173, a fresh market inbred line. Large differences were observed between the lines for both qualitative and quantitative morphological traits, suggesting that the ILs contain highly divergent allelic variation. Several loci contributing to unilateral incompatibility or hybrid necrosis were mapped with the lines. This IL population will facilitate studies of the S. sitiens genome and expands the range of genetic variation available for tomato breeding and research.     

Single nucleotide polymorphism genotyping using Kompetitive Allele Specific PCR (KASP): overview of the technology and its application in crop improvement

Single nucleotide polymorphism (SNP) data can be obtained using one of the numerous uniplex or multiplex SNP genotyping platforms that combine a variety of chemistries, detection methods, and reaction formats. Kompetitive Allele Specific PCR (KASP) is one of the uniplex SNP genotyping platforms, and has evolved to be a global benchmark technology. However, there are no publications relating either to the technology itself or to its application in crop improvement programs. In this review, we provide an overview of the different aspects of the KASP genotyping platform, discuss its application in crop improvement, and compare it with the chip-based Illumina GoldenGate platform. The International Maize and Wheat Improvement Center routinely uses KASP, generating in excess of a million data points annually for crop improvement purposes. We found that (1) 81 % of the SNPs used in a custom-designed GoldenGate assay were transferable to KASP; (2) using KASP, negative controls (no template) consistently clustered together and rarely produced signals exceeding the threshold values for allele calling, in contrast to the situation observed using GoldenGate assays; (3) KASP’s average genotyping error in positive control DNA samples was 0.7–1.6 %, which is lower than that observed using GoldenGate (2.0–2.4 %); (4) KASP genotyping costs for marker-assisted recurrent selection were 7.9–46.1 % cheaper than those of the BeadXpress and GoldenGate platforms; and (5) KASP offers cost-effective and scalable flexibility in applications that require small to moderate numbers of markers, such as quality control analysis, quantitative trait loci (QTL) mapping in bi-parental populations, marker-assisted recurrent selection, marker-assisted backcrossing, and QTL fine mapping.     

Evaluation of genotyping concordance for commercial bovine SNP arrays using quality‐assurance samples

SNP arrays are widely used in genetic research and agricultural genomics applications, and the quality of SNP genotyping data is of paramount importance. In the present study, SNP genotyping concordance and discordance were evaluated for commercial bovine SNP arrays based on two types of quality assurance (QA) samples provided by Neogen GeneSeek. The genotyping discordance rates (GDRs) between chips were on average between 0.06% and 0.37% based on the QA type I data and between 0.05% and 0.15% based on the QA type II data. The average genotyping error rate (GER) pertaining to single SNP chips, based on the QA type II data, varied between 0.02% and 0.08% per SNP and between 0.01% and 0.06% per sample. These results indicate that genotyping concordance rate was high (i.e. from 99.63% to 99.99%). Nevertheless, mitochondrial and Y chromosome SNPs had considerably elevated GDRs and GERs compared to the SNPs on the 29 autosomes and X chromosome. The majority of genotyping errors resulted from single allotyping errors, which also included the opposite instances for allele ‘dropout’ (i.e. from AB to AA or BB). Simultaneous allotyping errors on both alleles (e.g. mistaking AA for BB or vice versa) were relatively rare. Finally, a list of SNPs with a GER greater than 1% is provided. Interpretation of association effects of these SNPs, for example in genome‐wide association studies, needs to be taken with caution. The genotyping concordance information needs to be considered in the optimal design of future bovine SNP arrays.     

A new multiplex SNP genotyping assay for detecting hybridization and introgression between the M and S molecular forms of Anopheles gambiae

The M and S forms of Anopheles gambiae have been the subject of intense study, but are morphologically indistinguishable and can only be identified using molecular techniques. PCR‐based assays to distinguish the two forms have been designed and applied widely. However, the application of these assays towards identifying hybrids between the two forms, and backcrossed hybrids in particular, has been problematic as the currently available diagnostic assays are based on single locus and/or are located within a multicopy gene. Here, we present an alternative genotyping method for detecting hybridization and introgression between M and S molecular forms based on a multilocus panel of single‐nucleotide polymorphisms (SNPs) fixed between the M and S forms. The panel of SNPs employed is located in so‐called islands of divergence leading us to describe this method as the ‘Divergence Island SNP’ (DIS) assay. We show this multilocus SNP genotyping approach can robustly and accurately detect F1 hybrids as well as backcrossed individuals.     

RNA‐Seq bulked segregant analysis enables the identification of high‐resolution genetic markers for breeding in hexaploid wheat

The identification of genetic markers linked to genes of agronomic importance is a major aim of crop research and breeding programmes. Here, we identify markers for Yr15, a major disease resistance gene for wheat yellow rust, using a segregating F₂ population. After phenotyping, we implemented RNA sequencing (RNA‐Seq) of bulked pools to identify single‐nucleotide polymorphisms (SNP) associated with Yr15. Over 27 000 genes with SNPs were identified between the parents, and then classified based on the results from the sequenced bulks. We calculated the bulk frequency ratio (BFR) of SNPs between resistant and susceptible bulks, selecting those showing sixfold enrichment/depletion in the corresponding bulks (BFR > 6). Using additional filtering criteria, we reduced the number of genes with a putative SNP to 175. The 35 SNPs with the highest BFR values were converted into genome‐specific KASP assays using an automated bioinformatics pipeline (PolyMarker) which circumvents the limitations associated with the polyploid wheat genome. Twenty‐eight assays were polymorphic of which 22 (63%) mapped in the same linkage group as Yr15. Using these markers, we mapped Yr15 to a 0.77‐cM interval. The three most closely linked SNPs were tested across varieties and breeding lines representing UK elite germplasm. Two flanking markers were diagnostic in over 99% of lines tested, thus providing a reliable haplotype for marker‐assisted selection in these breeding programmes. Our results demonstrate that the proposed methodology can be applied in polyploid F₂ populations to generate high‐resolution genetic maps across target intervals.     

Exploring the tertiary gene pool of bread wheat: sequence assembly and analysis of chromosome 5Mg of Aegilops geniculata

Next‐generation sequencing (NGS) provides a powerful tool for the discovery of important genes and alleles in crop plants and their wild relatives. Despite great advances in NGS technologies, whole‐genome shotgun sequencing is cost‐prohibitive for species with complex genomes. An attractive option is to reduce genome complexity to a single chromosome prior to sequencing. This work describes a strategy for studying the genomes of distant wild relatives of wheat by isolating single chromosomes from addition or substitution lines, followed by chromosome sorting using flow cytometry and sequencing of chromosomal DNA by NGS technology. We flow‐sorted chromosome 5M^g from a wheat/Aegilops geniculata disomic substitution line [DS5M^g (5D)] and sequenced it using an Illumina HiSeq 2000 system at approximately 50 × coverage. Paired‐end sequences were assembled and used for structural and functional annotation. A total of 4236 genes were annotated on 5M^g, in close agreement with the predicted number of genes on wheat chromosome 5D (4286). Single‐gene FISH indicated no major chromosomal rearrangements between chromosomes 5M^g and 5D. Comparing chromosome 5M^g with model grass genomes identified synteny blocks in Brachypodium distachyon, rice (Oryza sativa), sorghum (Sorghum bicolor) and barley (Hordeum vulgare). Chromosome 5M^g‐specific SNPs and cytogenetic probe‐based resources were developed and validated. Deletion bin‐mapped and ordered 5M^g SNP markers will be useful to track 5M‐specific introgressions and translocations. This study provides a detailed sequence‐based analysis of the composition of a chromosome from a distant wild relative of bread wheat, and opens up opportunities to develop genomic resources for wild germplasm to facilitate crop improvement.     

Development of a set of Triticum aestivum-Aegilops tauschii introgression lines

New wheat introgression lines were obtained which contain different segments of individual chromosomes of Aegilops tauschii in the Triticum aestivum cv. 'Chinese Spring' background. The introgression lines were developed to examine various subsets of alleles from the wild grass in the genetic background of common wheat. As starting point substitution lines of 'Chinese Spring' in which single chromosomes of the D genome had been replaced by homologous chromosomes of a synthetic wheat were used. Synthetic wheat had been obtained earlier from a cross between the tetraploid emmer (genomes AABB) and wild grass Aegilops tauschii (genome DD). The seven wheat chromosome substitution lines carrying different chromosomes of Ae. tauschii were crossed twice to T. aestivum cv. 'Chinese Spring' and 259 BC1-progeny plants were analysed. Phenotypic evaluation was carried out for different traits such as plant height, spikelet number, peduncle length, flowering time, spike length, tiller number, grain weight per ear, fertility and thousand kernel weight. Genotypic analysis was performed using a set of 65 microsatellite markers previously mapped on the chromosomes of the D genome of wheat. During this analysis recombinant lines carrying different segments of Ae. tauschii chromosomes were detected. Plants containing small introgressions of the alien genetic material were selfed to get homozygous lines and plants carrying large pieces of the donor chromosome were backcrossed again to get smaller introgressions. Further microsatellite analysis of selected BC1F2-progeny plants resulted in detection of a first set of 36 homozygous lines carrying different pieces of Ae. tauschii genome.     

Characterization of polyploid wheat genomic diversity using a high‐density 90 000 single nucleotide polymorphism array

High‐density single nucleotide polymorphism (SNP) genotyping arrays are a powerful tool for studying genomic patterns of diversity, inferring ancestral relationships between individuals in populations and studying marker–trait associations in mapping experiments. We developed a genotyping array including about 90 000 gene‐associated SNPs and used it to characterize genetic variation in allohexaploid and allotetraploid wheat populations. The array includes a significant fraction of common genome‐wide distributed SNPs that are represented in populations of diverse geographical origin. We used density‐based spatial clustering algorithms to enable high‐throughput genotype calling in complex data sets obtained for polyploid wheat. We show that these model‐free clustering algorithms provide accurate genotype calling in the presence of multiple clusters including clusters with low signal intensity resulting from significant sequence divergence at the target SNP site or gene deletions. Assays that detect low‐intensity clusters can provide insight into the distribution of presence–absence variation (PAV) in wheat populations. A total of 46 977 SNPs from the wheat 90K array were genetically mapped using a combination of eight mapping populations. The developed array and cluster identification algorithms provide an opportunity to infer detailed haplotype structure in polyploid wheat and will serve as an invaluable resource for diversity studies and investigating the genetic basis of trait variation in wheat.     

Construction of Agropyron Gaertn. genetic linkage maps using a wheat 660K SNP array reveals a homoeologous relationship with the wheat genome

Agropyron Gaertn. (P genome) is a wild relative of wheat that harbours many genetic variations that could be used to increase the genetic diversity of wheat. To agronomically transfer important genes from the P genome to a wheat chromosome by induced homoeologous pairing and recombination, it is necessary to determine the chromosomal relationships between Agropyron and wheat. Here, we report using the wheat 660K single nucleotide polymorphism (SNP) array to genotype a segregating Agropyron F₁ population derived from an interspecific cross between two cross‐pollinated diploid collections ‘Z1842’ [A. cristatum (L.) Beauv.] (male parent) and ‘Z2098’ [A. mongolicum Keng] (female parent) and 35 wheat–A. cristatum addition/substitution lines. Genetic linkage maps were constructed using 913 SNP markers distributed among seven linkage groups spanning 839.7 cM. The average distance between adjacent markers was 1.8 cM. The maps identified the homoeologous relationship between the P genome and wheat and revealed that the P and wheat genomes are collinear and relatively conserved. In addition, obvious rearrangements and introgression spread were observed throughout the P genome compared with the wheat genome. Combined with genotyping data, the complete set of wheat–A. cristatum addition/substitution lines was characterized according to their homoeologous relationships. In this study, the homoeologous relationship between the P genome and wheat was identified using genetic linkage maps, and the detection mean for wheat–A. cristatum introgressions might significantly accelerate the introgression of genetic variation from Agropyron into wheat for exploitation in wheat improvement programmes.     

Conversion of array‐based single nucleotide polymorphic markers for use in targeted genotyping by sequencing in hexaploid wheat (Triticum aestivum)

Wheat breeders and academics alike use single nucleotide polymorphisms (SNP s) as molecular markers to characterize regions of interest within the hexaploid wheat genome. A number of SNP ‐based genotyping platforms are available, and their utility depends upon factors such as the available technologies, number of data points required, budgets and the technical expertise required. Unfortunately, markers can rarely be exchanged between existing and newly developed platforms, meaning that previously generated data cannot be compared, or combined, with more recently generated data sets. We predict that genotyping by sequencing will become the predominant genotyping technology within the next 5–10 years. With this in mind, to ensure that data generated from current genotyping platforms continues to be of use, we have designed and utilized SNP ‐based capture probes from several thousand existing and publicly available probes from Axiom® and KASP ? genotyping platforms. We have validated our capture probes in a targeted genotyping by sequencing protocol using 31 previously genotyped UK elite hexaploid wheat accessions. Data comparisons between targeted genotyping by sequencing, Axiom® array genotyping and KASP ? genotyping assays, identified a set of 3256 probes which reliably bring together targeted genotyping by sequencing data with the previously available marker data set. As such, these probes are likely to be of considerable value to the wheat community. The probe details, full probe sequences and a custom built analysis pipeline may be freely downloaded from the CerealsDB website (http://www.cerealsdb.uk.net/cerealgenomics/CerealsDB /sequence_capture.php).     

Meiotic homoeologous recombination-based mapping of wheat chromosome 2B and its homoeologues in <Emphasis Type="Italic">Aegilops speltoides</Emphasis> and <Emphasis Type="Italic">Thinopyrum elongatum</Emphasis>

Key message

We physically dissected and mapped wheat chromosome 2B and its homoeologues in Aegilops speltoides and Thinopyrum elongatum based on meiotic homoeologous recombination, providing a unique physical framework for genome studies.

Abstract

Common wheat has a large and complex genome with narrow genetic diversity and various degrees of recombination between the A, B, and D subgenomes. This has limited the homologous recombination-based genome studies in wheat. Here, we exploited meiotic homoeologous recombination for molecular mapping of wheat chromosome 2B and its homoeologue 2S from Aegilops speltoides and 2E from Thinopyrum elongatum. The 2B–2S and 2B–2E recombination was induced by the ph1b mutant, and recovered using molecular markers and fluorescent genomic in situ hybridization (FGISH). A total of 112 2B–2S and 87 2B–2E recombinants involving different chromosome regions were developed and physically delineated by FGISH. The 2B–2S and 2B–2E recombination hotspots mapped to the subterminal regions on both arms. Recombination hotspots with the highest recombination rates mapped to the short arms. Eighty-three 2B–2S and 67 2B–2E recombinants were genotyped using the wheat 90 K SNP arrays. Based on the genotyping results and FGISH patterns of the recombinants, chromosomes 2B, 2S, and 2E were partitioned into 93, 66, and 46 bins, respectively. In total, 1037 SNPs physically mapped onto distinct bins of these three homoeologous chromosomes. A homoeologous recombination-based bin map was constructed for chromosome 2B, providing a unique physical framework for genome studies in wheat and its relatives. Meiotic homoeologous recombination also facilitates gene introgression to diversify the wheat genome for germplasm development. Therefore, homoeologous recombination-based studies enhance understanding of the wheat genome and its homoeologous counterparts from wild grasses, and expand the genetic variability of the wheat genome.

     

In silico single nucleotide polymorphism discovery and application to marker-assisted selection in soybean

The general approach to discovering single nucleotide polymorphisms (SNPs) requires locus-specific PCR amplification. To enhance the efficiency of SNP discovery in soybean, we used in silico analysis prior to re-sequencing as it is both rapid and inexpensive. In silico analysis was performed to detect putative SNPs in expressed sequence tag (EST) contigs assembled using publicly available ESTs from 18 different soybean genotypes. SNP validation by direct sequencing of six soybean cultivars and a wild soybean genotype was performed with PCR primers designed from EST contigs aligned with at least 5 out of 18 soybean genotypes. The efficiency of SNP discovery among the confirmation genotypes was 81.2%. Furthermore, the efficiency of SNP discovery between Pureunkong and Jinpumkong 2 genotypes was 47.4%, a great improvement on our previous finding based on direct sequencing (22.3%). Using SNPs between Pureunkong and Jinpumkong 2 in EST contigs, which were linked to target traits, we were able to genotype 90 recombinant inbred lines by high-resolution melting (HRM) analysis. These SNPs were mapped onto the expected locations near quantitative trait loci for water-logging tolerance and seed pectin concentration. Thus, our protocol for HRM analysis can be applied successfully not only to genetic diversity studies, but also to marker-assisted selection (MAS). Our study suggests that a combination of in silico analysis and HRM can reduce the cost and labor involved in developing SNP markers and genotyping SNPs. The markers developed in this study can also easily be applied to MAS if the markers are associated with the target traits.     

Development,validation and genetic analysis of a large soybean SNP genotyping array

Cultivated soybean (Glycine max) suffers from a narrow germplasm relative to other crop species, probably because of under‐use of wild soybean (Glycine soja) as a breeding resource. Use of a single nucleotide polymorphism (SNP) genotyping array is a promising method for dissecting cultivated and wild germplasms to identify important adaptive genes through high‐density genetic mapping and genome‐wide association studies. Here we describe a large soybean SNP array for use in diversity analyses, linkage mapping and genome‐wide association analyses. More than four million high‐quality SNPs identified from high‐depth genome re‐sequencing of 16 soybean accessions and low‐depth genome re‐sequencing of 31 soybean accessions were used to select 180 961 SNPs for creation of the Axiom^® SoyaSNP array. Validation analysis for a set of 222 diverse soybean lines showed that 170 223 markers were of good quality for genotyping. Phylogenetic and allele frequency analyses of the validation set data indicated that accessions showing an intermediate morphology between cultivated and wild soybeans collected in Korea were natural hybrids. More than 90 unanchored scaffolds in the current soybean reference sequence were assigned to chromosomes using this array. Finally, dense average spacing and preferential distribution of the SNPs in gene‐rich chromosomal regions suggest that this array may be suitable for genome‐wide association studies of soybean germplasm. Taken together, these results suggest that use of this array may be a powerful method for soybean genetic analyses relating to many aspects of soybean breeding.