High-throughput single nucleotide polymorphism genotyping in wheat (Triticum spp.)

Reviewing the essential roles of remote phenotyping, GWAS and explainable AI in practical marker-assisted selection for drought-tolerant winter wheat breeding

Marker-assisted selection (MAS) plays a crucial role in crop breeding improving the speed and precision of conventional breeding programmes by quickly and reliably identifying and selecting plants with desired traits. However, the efficacy of MAS depends on several prerequisites, with precise phenotyping being a key aspect of any plant breeding programme. Recent advancements in high-throughput remote phenotyping, facilitated by unmanned aerial vehicles coupled to machine learning, offer a non-destructive and efficient alternative to traditional, time-consuming, and labour-intensive methods. Furthermore, MAS relies on knowledge of marker-trait associations, commonly obtained through genome-wide association studies (GWAS), to understand complex traits such as drought tolerance, including yield components and phenology. However, GWAS has limitations that artificial intelligence (AI) has been shown to partially overcome. Additionally, AI and its explainable variants, which ensure transparency and interpretability, are increasingly being used as recognised problem-solving tools throughout the breeding process. Given these rapid technological advancements, this review provides an overview of state-of-the-art methods and processes underlying each MAS, from phenotyping, genotyping and association analyses to the integration of explainable AI along the entire workflow. In this context, we specifically address the challenges and importance of breeding winter wheat for greater drought tolerance with stable yields, as regional droughts during critical developmental stages pose a threat to winter wheat production. Finally, we explore the transition from scientific progress to practical implementation and discuss ways to bridge the gap between cutting-edge developments and breeders, expediting MAS-based winter wheat breeding for drought tolerance.

Development of breeder chip for gene detection and molecular-assisted selection by target sequencing in wheat

Wheat is an essential food crop and its high and stable yield is suffering from great challenges due to the limitations of current breeding technology and various stresses. Accelerating molecularly assisted stress-resistance breeding is critical. Through a meta-analysis of published loci in wheat over the last two decades, we selected 60 loci with main breeding objectives, high heritability, and reliable genotyping, such as stress resistance, yield, plant height, and resistance to spike germination. Then, using genotyping by target sequencing (GBTS) technology, we developed a liquid phase chip based on 101 functional or closely linked markers. The genotyping of 42 loci was confirmed in an extensive collection of Chinese wheat cultivars, indicating that the chip can be used in molecular-assisted selection (MAS) for target breeding goals. Besides, we can perform the preliminary parentage analysis with the genotype data. The most significant contribution of this work lies in translating a large number of molecular markers into a viable chip and providing reliable genotypes. Breeders can quickly screen germplasm resources, parental breeding materials, and intermediate materials for the presence of excellent allelic variants using the genotyping data by this chip, which is high throughput, convenient, reliable, and cost-efficient.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01359-3.

Role of NGS and SNP genotyping methods in sugarcane improvement programs

Sugarcane (Saccharum spp.) is one of the most economically significant crops because of its high sucrose content and it is a promising biomass feedstock for biofuel production. Sugarcane genome sequencing and analysis is a difficult task due to its heterozygosity and polyploidy. Long sequence read technologies, PacBio Single-Molecule Real-Time (SMRT) sequencing, the Illumina TruSeq, and the Oxford Nanopore sequencing could solve the problem of genome assembly. On the applications side, next generation sequencing (NGS) technologies played a major role in the discovery of single nucleotide polymorphism (SNP) and the development of low to high throughput genotyping platforms. The two mainstream high throughput genotyping platforms are the SNP microarray and genotyping by sequencing (GBS). This paper reviews the NGS in sugarcane genomics, genotyping methodologies, and the choice of these methods. Array-based SNP genotyping is robust, provides consistent SNPs, and relatively easier downstream data analysis. The GBS method identifies large scale SNPs across the germplasm. A combination of targeted GBS and array-based genotyping methods should be used to increase the accuracy of genomic selection and marker-assisted breeding.

Genome-wide introgression from a bread wheat × Lophopyrum elongatum amphiploid into wheat

We introgressed wheatgrass germplasm from the octoploid amphiploid Triticum aestivum× Lophopyrum elongatum into wheat by manipulating the wheat Ph1 gene and discovered and characterized 130 introgression lines harboring single or, in various combinations, complete and recombined L. elongatum chromosomes. Diploid wheatgrass Lophopyrum elongatum (genomes EE) possesses valuable traits for wheat genetics and breeding. We evaluated several strategies for introgression of this germplasm into wheat. To detect it, we developed and validated multiplexed sets of Sequenom MassARRAY single nucleotide polymorphism (SNP) markers, which differentiated disomic and monosomic L. elongatum chromosomes from wheat chromosomes. We identified 130 introgression lines (ILs), which harbored 108 complete and 89 recombined L. elongatum chromosomes. Of the latter, 59 chromosomes were recombined by one or more crossovers and 30 were involved in centromeric (Robertsonian) translocations or were telocentric. To identify wheat chromosomes substituted for or recombined with L. elongatum chromosomes, we genotyped the ILs with the wheat 90-K Infinium SNP array. We found that most of the wheat 90-K probes correctly detected their targets in the L. elongatum genome and showed that some wheat SNPs are ancient and had originated prior to the divergence of the wheat and L. elongatum lineages. Of the 130 ILs, 52% were homozygous for Ph1 deletion and thus are staged to be recombined further. We failed to detect in the L. elongatum genome the 4/5 reciprocal translocation that has been reported in Thinopyrum bessarabicum and several other Triticeae genomes.

Genetic Diversity and Population Structure in a Vitis spp. Core Collection Investigated by SNP Markers

Single nucleotide polymorphism (SNP) genotyping arrays are powerful tools to measure the level of genetic polymorphism within a population. The coming of next-generation sequencing technologies led to identifying thousands and millions of SNP loci useful in assessing the genetic diversity. The Vitis genotyping array, containing 18k SNP loci, has been developed and used to detect genetic diversity of Vitis vinifera germplasm. So far, this array was not validated on non-vinifera genotypes used as grapevine rootstocks. In this work, a core collection of 70 grapevine rootstocks, composed of individuals belonging to Vitis species not commonly used in the breeding programs, was genotyped using the 18k SNP genotyping array. SNP results were compared to the established SSR (Simple Sequence Repeat) markers in terms of heterozygosity and genetic structure of the core collection. Genotyping array has proved to be a valuable tool for genotyping of grapevine rootstocks, with more than 90% of SNPs successfully amplified. Structure analysis detected a high degree of admixed genotypes, supported by the complex genetic background of non-vinifera germplasm. Moreover, SNPs clearly differentiated non-vinifera and vinifera germplasm. These results represent a first step in studying the genetic diversity of non-conventional breeding material that will be used to select rootstocks with high tolerance to limiting environmental conditions.

SNP markers for early identification of high molecular weight glutenin subunits (HMW-GSs) in bread wheat

A set of eight SNP markers was developed to facilitate the early selection of HMW-GS alleles in breeding programmes. In bread wheat (Triticum aestivum), the high molecular weight glutenin subunits (HMW-GSs) are the most important determinants of technological quality. Known to be very diverse, HMW-GSs are encoded by the tightly linked genes Glu-1-1 and Glu-1-2. Alleles that improve the quality of dough have been identified. Up to now, sodium dodecyl sulphate-polyacrylamide gel electrophoresis (SDS-PAGE) of grain proteins is the most widely used for their identification. To facilitate the early selection of HMW-GS alleles in breeding programmes, we developed DNA-based molecular markers. For each accession of a core collection (n = 364 lines) representative of worldwide bread wheat diversity, HMW-GSs were characterized by both genotyping and SDS-PAGE. Based on electrophoresis, we observed at least 8, 22 and 9 different alleles at the Glu-A1, Glu-B1 and Glu-D1 loci, respectively, including new variants. We designed a set of 17 single-nucleotide polymorphism (SNP) markers that were representative of the most frequent SDS-PAGE alleles at each locus. At Glu-A1 and Glu-D1, two and three marker-based haplotypes, respectively, captured the diversity of the SDS-PAGE alleles rather well. Discrepancies were found mainly for the Glu-B1 locus. However, statistical tests revealed that two markers at each Glu-B1 gene and their corresponding haplotypes were more significantly associated with the rheological properties of the dough than were the relevant SDS-PAGE alleles. To conclude, this study demonstrates that the SNP markers developed provide additional information on HMW-GS diversity. Two markers at Glu-A1, four at Glu-B1 and two at Glu-D1 constitute a useful toolbox for breeding wheat to improve end-use value.

Genetic diversity and linkage disequilibrium using SNP (KASP) and AFLP markers in a worldwide durum wheat (Triticum turgidum L. var durum) collection

The aim of this work was to analyze the genetic diversity and linkage disequilibrium in a collection of 168 durum wheat accessions (Triticum turgidum L. var. durum) of different origins. Our collection was mainly composed of released and unreleased Argentinian germplasm, with additional genotypes from Italy, Chile, France, CIMMYT, Cyprus, USA and WANA region. To this end, the entire collection was characterized with 85 Single Nucleotide Polymorphism (SNP) markers obtained by Kompetitive Allele Specific PCR (KASP), giving a heterozygosity (He) mean value of 0.183 and a coefficient of genetic differentiation (Gst) value of 0.139. A subset of 119 accessions was characterized with six Amplified Fragment Length Polymorphism (AFLP) primer combinations. A total of 181 polymorphic markers (125 AFLP and 56 SNP) amplified across this subset revealed He measures of 0.352 and 0.182, respectively. Of these, 134 were selected to estimate the genome-wide linkage disequilibrium obtaining low significant values (r2 = 0.11) in the subset, indicating its suitability for future genome-wide association studies (GWAS). The structure analysis conducted in the entire collection with SNP detected two subpopulations. However, the structure analysis conducted with AFLP markers in the subset of 119 accessions proved to have greater degree of resolution and detect six subpopulations. The information provided by both marker types was complementary and showed a strong association between old Argentinian and Italian germplasm and a contribution of CIMMYT germplasm to modern Argentinian, Chilean and Cypriot accessions. The influence of Mediterranean germplasm, mainly from Italy, on part of the modern Argentinian cultivars or breeding lines was also clearly evidenced. Although our analysis yields conclusive results and useful information for association mapping studies, further analyses are needed to refine the number of subpopulations present in the germplasm collection analyzed.

Identification of favorable SNP alleles and candidate genes responsible for inflorescence-related traits via GWAS in chrysanthemum

81 SNPs were identified for three inflorescence-related traits, in which 15 were highly favorable. Two dCAPS markers were developed for future MAS breeding, and six candidate genes were predicted. Chrysanthemum is a leading ornamental species worldwide and demonstrates a wealth of morphological variation. Knowledge about the genetic basis of its phenotypic variation for key horticultural traits can contribute to its effective management and genetic improvement. In this study, we conducted a genome-wide association study (GWAS) based on two years of phenotype data and a set of 92,617 single nucleotide polymorphisms (SNPs) using a panel of 107 diverse cut chrysanthemums to dissect the genetic control of three inflorescence-related traits. A total of 81 SNPs were significantly associated with the three inflorescence-related traits (capitulum diameter, number of ray florets and flowering time) in at least one environment, with an individual allele explaining 22.72-38.67% of the phenotypic variation. Fifteen highly favorable alleles were identified for the three target traits by computing the phenotypic effect values for the stable associations detected in 2 year-long trials at each locus. Dosage pyramiding effects of the highly favorable SNP alleles and significant linear correlations between highly favorable allele numbers and corresponding phenotypic performance were observed. Two highly favorable SNP alleles correlating to flowering time and capitulum diameter were converted to derived cleaved amplified polymorphic sequence (dCAPS) markers to facilitate future breeding. Finally, six putative candidate genes were identified that contribute to flowering time and capitulum diameter. These results serve as a foundation for analyzing the genetic mechanisms underlying important horticultural traits and provide valuable insights into molecular marker-assisted selection (MAS) in chrysanthemum breeding programs.

Characterization of a new Pm2 allele associated with broad-spectrum powdery mildew resistance in wheat line Subtil

Wheat powdery mildew is a severe disease affecting yield and quality. Host resistance was proved to be effective and environment-friendly. Wheat line Subtil is an elite germplasm resource resistant to 28 of 30 tested Bgt isolates. Genetic analysis showed that the powdery mildew resistance in Subtil was conferred by a single dominant gene, temporarily designated PmSub. Using bulked segregant analysis, PmSub was mapped to chromosome arm 5DS, and flanked by the markers Bwm16 and Cfd81/Bwm21 at 5.0 and 0.9 cM, respectively. Allelism tests further confirmed PmSub was allelic with documented Pm2 alleles. Then, homologous sequences of Pm2a related sequence was cloned from Subtil and Chinese Spring. It was completely identical to the reported Pm2a sequence, but significantly different from that of Chinese Spring. A marker SWGI067 was developed based on the sequence divergence of homologous sequence in Subtil and Chinese Spring. SWGI067 was closely linked to PmSub, indicating that the gene PmSub itself was different from the cloned Pm2a related sequence. Meanwhile, Subtil produced significantly different reaction pattern compared with other genotypes with Pm genes at or near Pm2 locus. Therefore, PmSub was most likely a new allele of Pm2. PmSub has opportunities for marker-assisted selecting for high-efficiency wheat improvement.

Crop Breeding Chips and Genotyping Platforms: Progress, Challenges, and Perspectives

There is a rapidly rising trend in the development and application of molecular marker assays for gene mapping and discovery in field crops and trees. Thus far, more than 50 SNP arrays and 15 different types of genotyping-by-sequencing (GBS) platforms have been developed in over 25 crop species and perennial trees. However, much less effort has been made on developing ultra-high-throughput and cost-effective genotyping platforms for applied breeding programs. In this review, we discuss the scientific bottlenecks in existing SNP arrays and GBS technologies and the strategies to develop targeted platforms for crop molecular breeding. We propose that future practical breeding platforms should adopt automated genotyping technologies, either array or sequencing based, target functional polymorphisms underpinning economic traits, and provide desirable prediction accuracy for quantitative traits, with universal applications under wide genetic backgrounds in crops. The development of such platforms faces serious challenges at both the technological level due to cost ineffectiveness, and the knowledge level due to large genotype-phenotype gaps in crop plants. It is expected that such genotyping platforms will be achieved in the next ten years in major crops in consideration of (a) rapid development in gene discovery of important traits, (b) deepened understanding of quantitative traits through new analytical models and population designs, (c) integration of multi-layer -omics data leading to identification of genes and pathways responsible for important breeding traits, and (d) improvement in cost effectiveness of large-scale genotyping. Crop breeding chips and genotyping platforms will provide unprecedented opportunities to accelerate the development of cultivars with desired yield potential, quality, and enhanced adaptation to mitigate the effects of climate change.

High-throughput SNP genotyping of modern and wild emmer wheat for yield and root morphology using a combined association and linkage analysis

Durum wheat (Triticum turgidum var. durum Desf.) is a major world crop that is grown primarily in areas of the world that experience periodic drought, and therefore, breeding climate-resilient durum wheat is a priority. High-throughput single nucleotide polymorphism (SNP) genotyping techniques have greatly increased the power of linkage and association mapping analyses for bread wheat, but as yet there is no durum wheat-specific platform available. In this study, we evaluate the new 384HT Wheat Breeders Array for its usefulness in tetraploid wheat breeding by genotyping a breeding population of F₆ hybrids, derived from multiple crosses between T. durum cultivars and wild and cultivated emmer wheat accessions. Using a combined linkage and association mapping approach, we generated a genetic map including 1345 SNP markers, and identified markers linked to 6 QTLs for coleoptile length (2), heading date (1), anthocyanin accumulation (1) and osmotic stress tolerance (2). We also developed a straightforward approach for combining genetic data from multiple families of limited size that will be useful for evaluating and mapping pre-existing breeding material.

Genetic analysis of a novel broad-spectrum powdery mildew resistance gene from the wheat-Agropyron cristatum introgression line Pubing 74

A novel broad-spectrum powdery mildew resistance gene PmPB74 was identified in wheat- Agropyron cristatum introgression line Pubing 74. Development of wheat cultivars with broad-spectrum, durable resistance to powdery mildew has been restricted by lack of superior genetic resources. In this study, a wheat-A. cristatum introgression line Pubing 74, originally selected from a wide cross between the common wheat cultivar Fukuhokomugi (Fukuho) and Agropyron cristatum (L.) Gaertn (2n = 4x = 28; genome PPPP), displayed resistance to powdery mildew at both the seedling and adult stages. The putative alien chromosomal fragment in Pubing 74 was below the detection limit of genomic in situ hybridization (GISH), but evidence for other non-GISH-detectable introgressions was provided by the presence of three STS markers specific to A. cristatum. Genetic analysis indicated that Pubing 74 carried a single dominant gene for powdery mildew resistance, temporarily designated PmPB74. Molecular mapping showed that PmPB74 was located on wheat chromosome arm 5DS, and flanked by markers Xcfd81 and HRM02 at genetic distances of 2.5 and 1.7 cM, respectively. Compared with other lines with powdery mildew resistance gene(s) on wheat chromosome arm 5DS, Pubing 74 was resistant to all 28 Blumeria graminis f. sp tritici (Bgt) isolates from different wheat-producing regions of northern China. Allelism tests indicated that PmPB74 was not allelic to PmPB3558 or Pm2. Our work showed that PmPB74 is a novel gene with broad resistance to powdery mildew, and hence will be helpful in broadening the genetic basis of powdery mildew resistance in wheat.

Characterization of a New Pm2 Allele Conferring Powdery Mildew Resistance in the Wheat Germplasm Line FG-1

Powdery mildew has a negative impact on wheat production. Novel host resistance increases the diversity of resistance genes and helps to control the disease. In this study, wheat line FG-1 imported from France showed a high level of powdery mildew resistance at both the seedling and adult stages. An F2 population and F2:3 families from the cross FG-1 × Mingxian 169 both fit Mendelian ratios for a single dominant resistance gene when tested against multiple avirulent Blumeria tritici f. sp. tritici (Bgt) races. This gene was temporarily designated PmFG. PmFG was mapped on the multi-allelic Pm2 locus of chromosome 5DS using seven SSR, 10 single nucleotide polymorphism (SNP)-derived and two SCAR markers with the flanking markers Xbwm21/Xcfd81/Xscar112 (distal) and Xbwm25 (proximal) at 0.3 and 0.5 cM being the closest. Marker SCAR203 co-segregated with PmFG. Allelism tests between PmFG and documented Pm2 alleles confirmed that PmFG was allelic with Pm2. Line FG-1 produced a significantly different reaction pattern compared to other lines with genes at or near Pm2 when tested against 49 Bgt isolates. The PmFG-linked marker alleles detected by the SNP-derived markers revealed significant variation between FG-1 and other lines with genes at or near Pm2. It was concluded that PmFG is a new allele at the Pm2 locus. Data from seven closely linked markers tested on 31 wheat cultivars indicated opportunities for marker-assisted pyramiding of this gene with other genes for powdery mildew resistance and additional traits.

Amand P, Bai G. 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.

Mapping quantitative trait loci for peroxidase activity and developing gene-specific markers for TaPod-A1 on wheat chromosome 3AL

Three novel QTL for peroxidase activity were mapped, and gene-specific markers for TaPod-A1 were developed and validated using RILs derived from the Doumai/Shi 4185 cross and 281 wheat cultivars. TaPod-A1 is within one of the three QTL. Peroxidase (POD) activity in grain is an important factor determining the color of flour and end-use products of wheat, such as noodles and steamed bread. Mapping QTL for POD activity, characterization of POD genes and development of gene-specific markers are important for molecular marker-assisted selection in wheat breeding. Quantitative trait loci (QTL) for POD activity in common wheat were mapped using a recombinant inbred line (RIL) population derived from a Doumai/Shi 4185 cross grown in four environments and genotyped using the wheat 90 K iSelect assay. Three novel QTL for POD activity, QPod.caas-3AL, QPod.caas-4BS and QPod.caas-5AS, were identified on chromosomes 3AL, 4BS and 5AS, explaining 5.3-21.2% of phenotypic variance across environments. The full-length genomic DNA (gDNA) sequence of a POD gene, designated TaPod-A1, on chromosome 3A was characterized by homolog cloning and PCR verification. Two complementary dominant sequence-tagged site (STS) markers, POD-3A1 and POD-3A2, were developed based on single nucleotide polymorphisms (SNPs) between two alleles at the TaPod-A1 locus, amplifying 291- and 766-bp fragments in cultivars with lower and higher POD activities, respectively. The two gene-specific markers were mapped on chromosome 3AL using a set of Chinese Spring (CS) nulli-tetrasomic lines, and ditelosomic lines 3AL and 3AS. QTL analysis indicated that QPod.caas-3AL co-segregated with the gene-specific markers POD-3A1 and POD-3A2. POD-3A1 and POD-3A2 were verified on 281 wheat cultivars and advanced lines, and showed significant (P < 0.05) associations with POD activities. POD-3A1 and POD-3A2 may be useful as markers for improving color attributes in wheat breeding programs.

High-density genetic linkage map construction and QTL mapping of grain shape and size in the wheat population Yanda1817 × Beinong6

High-density genetic linkage maps are necessary for precisely mapping quantitative trait loci (QTLs) controlling grain shape and size in wheat. By applying the Infinium iSelect 9K SNP assay, we have constructed a high-density genetic linkage map with 269 F ₈ recombinant inbred lines (RILs) developed between a Chinese cornerstone wheat breeding parental line Yanda1817 and a high-yielding line Beinong6. The map contains 2431 SNPs and 128 SSR & EST-SSR markers in a total coverage of 3213.2 cM with an average interval of 1.26 cM per marker. Eighty-eight QTLs for thousand-grain weight (TGW), grain length (GL), grain width (GW) and grain thickness (GT) were detected in nine ecological environments (Beijing, Shijiazhuang and Kaifeng) during five years between 2010–2014 by inclusive composite interval mapping (ICIM) (LOD≥2.5). Among which, 17 QTLs for TGW were mapped on chromosomes 1A, 1B, 2A, 2B, 3A, 3B, 3D, 4A, 4D, 5A, 5B and 6B with phenotypic variations ranging from 2.62% to 12.08%. Four stable QTLs for TGW could be detected in five and seven environments, respectively. Thirty-two QTLs for GL were mapped on chromosomes 1B, 1D, 2A, 2B, 2D, 3B, 3D, 4A, 4B, 4D, 5A, 5B, 6B, 7A and 7B, with phenotypic variations ranging from 2.62% to 44.39%. QGl.cau-2A.2 can be detected in all the environments with the largest phenotypic variations, indicating that it is a major and stable QTL. For GW, 12 QTLs were identified with phenotypic variations range from 3.69% to 12.30%. We found 27 QTLs for GT with phenotypic variations ranged from 2.55% to 36.42%. In particular, QTL QGt.cau-5A.1 with phenotypic variations of 6.82–23.59% was detected in all the nine environments. Moreover, pleiotropic effects were detected for several QTL loci responsible for grain shape and size that could serve as target regions for fine mapping and marker assisted selection in wheat breeding programs.

Identification and development of a functional marker from 6-SFT-A2 associated with grain weight in wheat

As a class of water-soluble, fructose-based oligo- and polysaccharides, fructans are major nonstructural carbohydrates and an important carbon source for grain filling in wheat (Triticum aestivum L.). Four enzymes are involved in fructan synthesis in higher plants, and 6-SFT is a key enzyme in fructan biosynthesis. In this study, thirteen single nucleotide polymorphisms were detected in 6-SFT-A2 in 24 wheat accessions, forming three haplotypes. Two cleaved amplified polymorphic sequence markers developed based on polymorphisms at sites 1870(A-G) and 1951(A-G) distinguished the three haplotypes. 6-SFT-A2 was located on chromosome 4A, between markers P2454.3 and P3465.1 in a doubled haploid (DH) population derived from the cross Hanxuan 10 × Lumai 14. The DH population comprising 150 lines and a historical population consisting of 154 accessions were used in a 6-SFT-A2 marker-trait association analysis. The three haplotypes were significantly associated with thousand-grain weight (TGW) under rainfed conditions. HapIII had a significant positive effect on TGW. There were significant differences between the Hanxuan 10 and Lumai 14 genotypes in both rainfed and irrigated environments. The average TGW of Lumai 14 (HapIII) was higher than that of Hanxuan 10 (HapI). The frequencies of 6-SFT-A2HapIII in cultivars released at different periods showed that it had been strongly positively selected in breeding programs. The preferred HapIII for TGW occurred at higher frequencies in Gansu, Beijing, Shanxi, and Hebei than other regions in northern China.

Quantitative SNP genotyping of polyploids with MassARRAY and other platforms

Accurate genotyping is essential for building genetic maps and performing genome assembly of polyploid species. Recent high-throughput techniques, such as Illumina GoldenGate™ and Sequenom iPLEX MassARRAY®, have made it possible to accurately estimate the relative abundances of different alleles even when the ploidy of the population is unknown. Here we describe the experimental methods for collecting these relative allele intensities and then demonstrate the practical concerns for inferring genotypes using Bayesian inference via the software package SuperMASSA.

Application of next-generation sequencing technology to study genetic diversity and identify unique SNP markers in bread wheat from Kazakhstan

BACKGROUND

New SNP marker platforms offer the opportunity to investigate the relationships between wheat cultivars from different regions and assess the mechanism and processes that have led to adaptation to particular production environments. Wheat breeding has a long history in Kazakhstan and the aim of this study was to explore the relationship between key varieties from Kazakhstan and germplasm from breeding programs for other regions.

RESULTS

The study revealed 5,898 polymorphic markers amongst ten cultivars, of which 2,730 were mapped in the consensus genetic map. Mapped SNP markers were distributed almost equally across the A and B genomes, with between 279 and 484 markers assigned to each chromosome. Marker coverage was approximately 10-fold lower in the D genome. There were 863 SNP markers identified as unique to specific cultivars, and clusters of these markers (regions containing more than three closely mapped unique SNPs) showed specific patterns on the consensus genetic map for each cultivar. Significant intra-varietal genetic polymorphism was identified in three cultivars (Tzelinnaya 3C, Kazakhstanskaya rannespelaya and Kazakhstanskaya 15). Phylogenetic analysis based on inter-varietal polymorphism showed that the very old cultivar Erythrospermum 841 was the most genetically distinct from the other nine cultivars from Kazakhstan, falling in a clade together with the American cultivar Sonora and genotypes from Central and South Asia. The modern cultivar Kazakhstanskaya 19 also fell into a separate clade, together with the American cultivar Thatcher. The remaining eight cultivars shared a single sub-clade but were categorised into four clusters.

CONCLUSION

The accumulated data for SNP marker polymorphisms amongst bread wheat genotypes from Kazakhstan may be used for studying genetic diversity in bread wheat, with potential application for marker-assisted selection and the preparation of a set of genotype-specific markers.

Single-nucleotide polymorphisms and association analysis of drought-resistance gene TaSnRK2.8 in common wheat

TaSnRK2.8, an SnRK2 (sucrose non-fermenting1-related protein kinase 2) member of wheat, confers enhanced multi-stress tolerances in carbohydrate metabolism. In the study, two types of genomic sequences of TaSnRK2.8 were detected in common wheat. Sequencing analysis showed that there was a variation-enriched region, designated TaSnRK2.8-A-C, covering the eighth intron, the ninth exon and the 3'-flanking region of TaSnRK2.8-A, and no divergence occurred in TaSnRK2.8-B. Single nucleotide polymorphisms in the TaSnRK2.8-A-C region were investigated in 165 wheat accessions. Three of 751 sequenced nucleotide sites were polymorphic. Nucleotide diversity (π) in the region was 0.00068. Sliding-window analysis demonstrated that the nucleotide diversity was highest in the 3'-flanking sequence. As predicted, the highly frequent SNP was significantly associated with seedling biomass under normal conditions, plant height, flag leaf width and water-soluble carbohydrate content under drought conditions. Analysis of variance of correlated traits between accessions with the A and G genotypes indicated that the A variant was the more favorable allele associated with significantly increased seedling biomass and water-soluble carbohydrates. Based on the SNP, we developed a functional marker of TaSnRK2.8-A-C, that could be utilized in wheat breeding programs aimed at improving seedling biomass and water-soluble carbohydrates, and consequently to enhance stress resistance in wheat.

Polymorphism of TaSAP1-A1 and its association with agronomic traits in wheat

TaSAP1, a member of the stress association protein (SAP) gene family from wheat, is involved in response to several abiotic stresses, including drought, salt and cold. In this study, TaSAP1-A1, a TaSAP1 member on chromosome 7A, and its flanking sequences were isolated. Polymorphism analysis indicated that the average nucleotide diversity (π) of the whole region was 0.00296. The highest nucleotide diversity occurred in the promoter region (π = 0.00631) and no polymorphism was identified in the coding region. Three markers T7AM5, T7AM2606 and T7AM39 located in the promoter region, were developed from sequence variations (InDel5-1810, SNP-2606 and InDel39-1637). Six haplotypes were identified among 300 accessions based on the three markers. TaSAP1-A1 was located on chromosome 7A using marker T7AM39 and was flanked by markers Xwmc530 and Xbarc174. QTL for yield-related traits, including 1,000-grain weight, number of grains per spike and grain yield, were located in the same region. In marker- and haplotype-trait association analyses, TaSAP1-A1 was significantly associated with 1,000-grain weight, number of grains per spike, spike length, peduncle length and total number of spikelets per spike in multiple environments. These results provide useful information for marker-assisted selection for yield-related traits under well-watered and drought-stressed conditions.

Assignment of SNP allelic configuration in polyploids using competitive allele-specific PCR: application to citrus triploid progeny

BACKGROUND

Polyploidy is a major component of eukaryote evolution. Estimation of allele copy numbers for molecular markers has long been considered a challenge for polyploid species, while this process is essential for most genetic research. With the increasing availability and whole-genome coverage of single nucleotide polymorphism (SNP) markers, it is essential to implement a versatile SNP genotyping method to assign allelic configuration efficiently in polyploids.

SCOPE

This work evaluates the usefulness of the KASPar method, based on competitive allele-specific PCR, for the assignment of SNP allelic configuration. Citrus was chosen as a model because of its economic importance, the ongoing worldwide polyploidy manipulation projects for cultivar and rootstock breeding, and the increasing availability of SNP markers.

CONCLUSIONS

Fifteen SNP markers were successfully designed that produced clear allele signals that were in agreement with previous genotyping results at the diploid level. The analysis of DNA mixes between two haploid lines (Clementine and pummelo) at 13 different ratios revealed a very high correlation (average = 0·9796; s.d. = 0·0094) between the allele ratio and two parameters [θ angle = tan(-1) (y/x) and y' = y/(x + y)] derived from the two normalized allele signals (x and y) provided by KASPar. Separated cluster analysis and analysis of variance (ANOVA) from mixed DNA simulating triploid and tetraploid hybrids provided 99·71 % correct allelic configuration. Moreover, triploid populations arising from 2n gametes and interploid crosses were easily genotyped and provided useful genetic information. This work demonstrates that the KASPar SNP genotyping technique is an efficient way to assign heterozygous allelic configurations within polyploid populations. This method is accurate, simple and cost-effective. Moreover, it may be useful for quantitative studies, such as relative allele-specific expression analysis and bulk segregant analysis.

Genetic diversity revealed by single nucleotide polymorphism markers in a worldwide germplasm collection of durum wheat

Evaluation of genetic diversity and genetic structure in crops has important implications for plant breeding programs and the conservation of genetic resources. Newly developed single nucleotide polymorphism (SNP) markers are effective in detecting genetic diversity. In the present study, a worldwide durum wheat collection consisting of 150 accessions was used. Genetic diversity and genetic structure were investigated using 946 polymorphic SNP markers covering the whole genome of tetraploid wheat. Genetic structure was greatly impacted by multiple factors, such as environmental conditions, breeding methods reflected by release periods of varieties, and gene flows via human activities. A loss of genetic diversity was observed from landraces and old cultivars to the modern cultivars released during periods of the Early Green Revolution, but an increase in cultivars released during the Post Green Revolution. Furthermore, a comparative analysis of genetic diversity among the 10 mega ecogeographical regions indicated that South America, North America, and Europe possessed the richest genetic variability, while the Middle East showed moderate levels of genetic diversity.

Phenotypic diversity and association mapping for fruit quality traits in cultivated tomato and related species

Association mapping has been proposed as an efficient approach to assist in the identification of the molecular basis of agronomical traits in plants. For this purpose, we analyzed the phenotypic and genetic diversity of a large collection of tomato accessions including 44 heirloom and vintage cultivars (Solanum lycopersicum), 127 S. lycopersicum var. cerasiforme (cherry tomato) and 17 Solanum pimpinellifolium accessions. The accessions were genotyped using a SNPlex™ assay of 192 SNPs, among which 121 were informative for subsequent analysis. Linkage disequilibrium (LD) of pairwise loci and population structure were analyzed, and the association analysis between SNP genotypes and ten fruit quality traits was performed using a mixed linear model. High level of LD was found in the collection at the whole genome level. It was lower when considering only the 127 S. lycopersicum var. cerasiforme accessions. Genetic structure analysis showed that the population was structured into two main groups, corresponding to cultivated and wild types and many intermediates. The number of associations detected per trait varied, according to the way the structure was taken into account, with 0-41 associations detected per trait in the whole collection and a maximum of four associations in the S. lycopersicum var. cerasiforme accessions. A total of 40 associations (30 %) were co-localized with previously identified quantitative trait loci. This study thus showed the potential and limits of using association mapping in tomato populations.

Endosperm tolerance of paternal aneuploidy allows radiation hybrid mapping of the wheat D-genome and a measure of γ ray-induced chromosome breaks

Physical mapping and genome sequencing are underway for the ≈17 Gb wheat genome. Physical mapping methods independent of meiotic recombination, such as radiation hybrid (RH) mapping, will aid precise anchoring of BAC contigs in the large regions of suppressed recombination in Triticeae genomes. Reports of endosperm development following pollination with irradiated pollen at dosages that cause embryo abortion prompted us to investigate endosperm as a potential source of RH mapping germplasm. Here, we report a novel approach to construct RH based physical maps of all seven D-genome chromosomes of the hexaploid wheat ‘Chinese Spring’, simultaneously. An 81-member subset of endosperm samples derived from 20-Gy irradiated pollen was genotyped for deletions, and 737 markers were mapped on seven D-genome chromosomes. Analysis of well-defined regions of six chromosomes suggested a map resolution of ∼830 kb could be achieved; this estimate was validated with assays of markers from a sequenced contig. We estimate that the panel contains ∼6,000 deletion bins for D-genome chromosomes and will require ∼18,000 markers for high resolution mapping. Map-based deletion estimates revealed a majority of 1–20 Mb interstitial deletions suggesting mutagenic repair of double-strand breaks in pollen provides a useful resource for RH mapping and map based cloning studies.

A multiparent advanced generation inter-cross population for genetic analysis in wheat

We present the first results from a novel multiparent advanced generation inter-cross (MAGIC) population derived from four elite wheat cultivars. The large size of this MAGIC population (1579 progeny), its diverse genetic composition and high levels of recombination all contribute to its value as a genetic resource. Applications of this resource include interrogation of the wheat genome and the analysis of gene-trait association in agronomically important wheat phenotypes. Here, we report the utilization of a MAGIC population for the first time for linkage map construction. We have constructed a linkage map with 1162 DArT, single nucleotide polymorphism and simple sequence repeat markers distributed across all 21 chromosomes. We benchmark this map against a high-density DArT consensus map created by integrating more than 100 biparental populations. The linkage map forms the basis for further exploration of the genetic architecture within the population, including characterization of linkage disequilibrium, founder contribution and inclusion of an alien introgression into the genetic map. Finally, we demonstrate the application of the resource for quantitative trait loci mapping using the complex traits plant height and hectolitre weight as a proof of principle.

Identification, characterization and interpretation of single-nucleotide sequence variation in allopolyploid crop species

An understanding of nature and extent of nucleotide sequence variation is required for programmes of discovery and characterization of single nucleotide polymorphisms (SNPs), which provide the most versatile class of molecular genetic marker. A majority of higher plant species are polyploids, and allopolyploidy, because of hybrid formation between closely related taxa, is very common. Mutational variation may arise both between allelic (homologous) sequences within individual subgenomes and between homoeologous sequences among subgenomes, in addition to paralogous variation between duplicated gene copies. Successful SNP validation in allopolyploids depends on differentiation of the sequence variation classes. A number of biological factors influence the feasibility of discrimination, including degree of gene family complexity, inbreeding or outbreeding reproductive habit, and the level of knowledge concerning progenitor diploid species. In addition, developments in high-throughput DNA sequencing and associated computational analysis provide general solutions for the genetic analysis of allopolyploids. These issues are explored in the context of experience from a range of allopolyploid species, representing grain (wheat and canola), forage (pasture legumes and grasses), and horticultural (strawberry) crop. Following SNP discovery, detection in routine genotyping applications also presents challenges for allopolyploids. Strategies based on either design of subgenome-specific SNP assays through homoeolocus-targeted polymerase chain reaction (PCR) amplification, or detection of incremental changes in nucleotide variant dosage, are described.

Genome-wide association analysis to identify chromosomal regions determining components of earliness in wheat

The modification of flowering date is considered an important way to escape the current or future climatic constraints that affect wheat crops. A better understanding of its genetic bases would enable a more efficient and rapid modification through breeding. The objective of this study was to identify chromosomal regions associated with earliness in wheat. A 227-wheat core collection chosen to be highly contrasted for earliness was characterized for heading date. Experiments were conducted in controlled conditions and in the field for 3 years to break down earliness in the component traits: photoperiod sensitivity, vernalization requirement and narrow-sense earliness. Whole-genome association mapping was carried out using 760 molecular markers and taking into account the five ancestral group structure. We identified 62 markers individually associated to earliness components corresponding to 33 chromosomal regions. In addition, we identified 15 other significant markers and seven more regions by testing marker pair interactions. Co-localizations were observed with the Ppd-1, Vrn-1 and Rht-1 candidate genes. Using an independent set of lines to validate the model built for heading date, we were able to explain 34% of the variation using the structure and the significant markers. Results were compared with already published data using bi-parental populations giving an insight into the genetic architecture of flowering time in wheat.

SNP markers and their impact on plant breeding

The use of molecular markers has revolutionized the pace and precision of plant genetic analysis which in turn facilitated the implementation of molecular breeding of crops. The last three decades have seen tremendous advances in the evolution of marker systems and the respective detection platforms. Markers based on single nucleotide polymorphisms (SNPs) have rapidly gained the center stage of molecular genetics during the recent years due to their abundance in the genomes and their amenability for high-throughput detection formats and platforms. Computational approaches dominate SNP discovery methods due to the ever-increasing sequence information in public databases; however, complex genomes pose special challenges in the identification of informative SNPs warranting alternative strategies in those crops. Many genotyping platforms and chemistries have become available making the use of SNPs even more attractive and efficient. This paper provides a review of historical and current efforts in the development, validation, and application of SNP markers in QTL/gene discovery and plant breeding by discussing key experimental strategies and cases exemplifying their impact.

A taxonomist's view on genomic authentication

A brief history of taxonomy, for the most part plant oriented, is provided, which demonstrates the use of morphology early on, through the stages when different technologies became available at different times until the present use of genomic tools. Genomic authentication facilitates with greater precision than ever before the identification of an organism or part thereof. In this chapter I made an attempt to stress that, in general, but more so for genomic authentication, the use of the variation inherent in taxa down to the lowest level of the hierarchy of classification needs to be used to achieve a high degree of correct authentication.

Discovery and mapping of a new expressed sequence tag-single nucleotide polymorphism and simple sequence repeat panel for large-scale genetic studies and breeding of Theobroma cacao L

Theobroma cacao is an economically important tree of several tropical countries. Its genetic improvement is essential to provide protection against major diseases and improve chocolate quality. We discovered and mapped new expressed sequence tag-single nucleotide polymorphism (EST-SNP) and simple sequence repeat (SSR) markers and constructed a high-density genetic map. By screening 149 650 ESTs, 5246 SNPs were detected in silico, of which 1536 corresponded to genes with a putative function, while 851 had a clear polymorphic pattern across a collection of genetic resources. In addition, 409 new SSR markers were detected on the Criollo genome. Lastly, 681 new EST-SNPs and 163 new SSRs were added to the pre-existing 418 co-dominant markers to construct a large consensus genetic map. This high-density map and the set of new genetic markers identified in this study are a milestone in cocoa genomics and for marker-assisted breeding. The data are available at http://tropgenedb.cirad.fr.

Transcript-specific, single-nucleotide polymorphism discovery and linkage analysis in hexaploid bread wheat (Triticum aestivum L.)

Food security is a global concern and substantial yield increases in cereal crops are required to feed the growing world population. Wheat is one of the three most important crops for human and livestock feed. However, the complexity of the genome coupled with a decline in genetic diversity within modern elite cultivars has hindered the application of marker-assisted selection (MAS) in breeding programmes. A crucial step in the successful application of MAS in breeding programmes is the development of cheap and easy to use molecular markers, such as single-nucleotide polymorphisms. To mine selected elite wheat germplasm for intervarietal single-nucleotide polymorphisms, we have used expressed sequence tags derived from public sequencing programmes and next-generation sequencing of normalized wheat complementary DNA libraries, in combination with a novel sequence alignment and assembly approach. Here, we describe the development and validation of a panel of 1114 single-nucleotide polymorphisms in hexaploid bread wheat using competitive allele-specific polymerase chain reaction genotyping technology. We report the genotyping results of these markers on 23 wheat varieties, selected to represent a broad cross-section of wheat germplasm including a number of elite UK varieties. Finally, we show that, using relatively simple technology, it is possible to rapidly generate a linkage map containing several hundred single-nucleotide polymorphism markers in the doubled haploid mapping population of Avalon × Cadenza.

Sequence-based marker development in wheat: advances and applications to breeding

In the past two decades, the wheat community has made remarkable progress in developing molecular resources for breeding. A wide variety of molecular tools has been established to accelerate genetic and physical mapping for facilitating the efficient identification of molecular markers linked to genes and QTL of agronomic interest. Already, wheat breeders are benefiting from a wide range of techniques to follow the introgression of the most favorable alleles in elite material and develop improved varieties. Breeders soon will be able to take advantage of new technological developments based on Next Generation Sequencing. In this paper, we review the molecular toolbox available to wheat scientists and breeders for performing fundamental genomic studies and breeding. Special emphasis is given on the production and detection of single nucleotide polymorphisms (SNPs) that should enable a step change in saturating the wheat genome for more efficient genetic studies and for the development of new selection methods. The perspectives offered by the access to an ordered full genome sequence for further marker development and enhanced precision breeding is also discussed. Finally, we discuss the advantages and limitations of marker-assisted selection for supporting wheat improvement.

Genetic diversity analysis of abiotic stress response gene TaSnRK2.7-A in common wheat

Sucrose non-fermenting1-related protein kinase 2 (SnRK2) plays a key role in plant stress signaling transduction pathways. In this study, one copy of TaSnRK2.7, a SnRK2 member of common wheat, was isolated and characterized for nucleotide diversity among 45 wheat accessions with different stress-response features. Most of the accessions were elite wheat cultivars, which had been subject to population bottlenecks and intensive selection during breeding. Nucleotide and haplotype diversity across the entire TaSnRK2.7-A region was 0.00076 and 0.590, respectively, and diversity in non-coding regions was higher than that in coding regions. Sliding-window analysis showed variable levels of nucleotide variation along the entire TaSnRK2.7-A region; the sixth intron and ninth exon represented variation-enriched regions. As predicted, neutrality tests revealed that population bottlenecks or purifying selection had acted on the TaSnRK2.7-A gene, a relatively conserved gene. Furthermore, strong linkage disequilibrium between SNP loci extends across the entire TaSnRK2.7-A region. These findings demonstrate that the TaSnRK2.7-A genomic region has evolved under extensive selection pressure during crop breeding.

SNiPlay: a web-based tool for detection, management and analysis of SNPs. Application to grapevine diversity projects

BACKGROUND

High-throughput re-sequencing, new genotyping technologies and the availability of reference genomes allow the extensive characterization of Single Nucleotide Polymorphisms (SNPs) and insertion/deletion events (indels) in many plant species. The rapidly increasing amount of re-sequencing and genotyping data generated by large-scale genetic diversity projects requires the development of integrated bioinformatics tools able to efficiently manage, analyze, and combine these genetic data with genome structure and external data.

RESULTS

In this context, we developed SNiPlay, a flexible, user-friendly and integrative web-based tool dedicated to polymorphism discovery and analysis. It integrates:1) a pipeline, freely accessible through the internet, combining existing softwares with new tools to detect SNPs and to compute different types of statistical indices and graphical layouts for SNP data. From standard sequence alignments, genotyping data or Sanger sequencing traces given as input, SNiPlay detects SNPs and indels events and outputs submission files for the design of Illumina's SNP chips. Subsequently, it sends sequences and genotyping data into a series of modules in charge of various processes: physical mapping to a reference genome, annotation (genomic position, intron/exon location, synonymous/non-synonymous substitutions), SNP frequency determination in user-defined groups, haplotype reconstruction and network, linkage disequilibrium evaluation, and diversity analysis (Pi, Watterson's Theta, Tajima's D).Furthermore, the pipeline allows the use of external data (such as phenotype, geographic origin, taxa, stratification) to define groups and compare statistical indices.2) a database storing polymorphisms, genotyping data and grapevine sequences released by public and private projects. It allows the user to retrieve SNPs using various filters (such as genomic position, missing data, polymorphism type, allele frequency), to compare SNP patterns between populations, and to export genotyping data or sequences in various formats.

CONCLUSIONS

Our experiments on grapevine genetic projects showed that SNiPlay allows geneticists to rapidly obtain advanced results in several key research areas of plant genetic diversity. Both the management and treatment of large amounts of SNP data are rendered considerably easier for end-users through automation and integration. Current developments are taking into account new advances in high-throughput technologies.SNiPlay is available at: http://sniplay.cirad.fr/.

Association mapping of common bacterial blight resistance QTL in Ontario bean breeding populations

BACKGROUND

Common bacterial blight (CBB), incited by Xanthomonas axonopodis pv. phaseoli (Xap), is a major yield-limiting factor of common bean (Phaseolus vulgaris L.) production around the world. Host resistance is practically the most effective and environmentally-sound approach to control CBB. Unlike conventional QTL discovery strategies, in which bi-parental populations (F2, RIL, or DH) need to be developed, association mapping-based strategies can use plant breeding populations to synchronize QTL discovery and cultivar development.

RESULTS

A population of 469 dry bean lines of different market classes representing plant materials routinely developed in a bean breeding program were used. Of them, 395 lines were evaluated for CBB resistance at 14 and 21 DAI (Days After Inoculation) in the summer of 2009 in an artificially inoculated CBB nursery in south-western Ontario. All lines were genotyped using 132 SNPs (Single Nucleotide Polymorphisms) evenly distributed across the genome. Of the 132 SNPs, 26 SNPs had more than 20% missing data, 12 SNPs were monomorphic, and 17 SNPs had a MAF (Minor Allelic Frequency) of less than 0.20, therefore only 75 SNPs were used for association study, based on one SNP per locus. The best possible population structure was to assign 36% and 64% of the lines into Andean and Mesoamerican subgroups, respectively. Kinship analysis also revealed complex familial relationships among all lines, which corresponds with the known pedigree history. MLM (Mixed Linear Model) analysis, including population structure and kinship, was used to discover marker-trait associations. Eighteen and 22 markers were significantly associated with CBB rating at 14 and 21 DAI, respectively. Fourteen markers were significant for both dates and the markers UBC420, SU91, g321, g471, and g796 were highly significant (p ≤ 0.001). Furthermore, 12 significant SNP markers were co-localized with or close to the CBB-QTLs identified previously in bi-parental QTL mapping studies.

CONCLUSIONS

This study demonstrated that association mapping using a reasonable number of markers, distributed across the genome and with application of plant materials that are routinely developed in a plant breeding program can detect significant QTLs for traits of interest.

SNPs discovery and CAPS marker conversion in soybean

To discover new SNPs and develop an easy assay method in soybean, we compared the high-throughput pyrosequencing ESTs with whole genome sequences in different soybean varieties and identified 3899 SNPs. Transitions were found to be much more frequent than transversions in these SNPs. We found that SNPs were widely distributed in the soybean genome, targeting numerous genes involved in various physiological and biochemical processes influencing important agronomic traits. A set of 16 SNPs were validated in nine soybean varieties, and seven SNPs were converted into CAPS. From functional gene association analysis, the marker CAPS282 on the 3'-UTR of gene Glyma07g03490 was identified as associated with 100-seed weight in soybean. The SNP discovery and CAPS markers conversion system developed in this study is fast and cost effective, and holds great promise for molecular-assisted breeding of soybean.