BACKGROUND
progress in plant breeding focuses on the rapid development of new cultivars with improved attributes.
molecular markers allow breeders to characterize specific lines without the need for laborious and time-consuming phenotyping.
marker-assisted selection  is used in plant breeding to identify the allele present at a specific locus, allowing the breeder to select based on genotype  <cit> .
mas has been used for plant breeding in many crops to identify specific individuals with known genes of interest  <cit> , primarily to target large-effect, single targets  <cit> .
since each locus is generally genotyped independently, breeders tend to consider per data point costs when utilizing mas within breeding programs.
contemporary marker technologies for assaying single targets that are often used with mas include kasp, targeted amplicon sequencing, and snp arrays.
kasp  is a uniplex, fluorescence-based single nucleotide genotyping technology that utilizes allele-specific oligo extension  <cit> .
kasp markers have been used for breeding, qtl mapping, and are the main genotyping platform for the generation challenge program at cimmyt  <cit> .
the arrival of inexpensive sequencing has led to the development of economical sequence-based genotyping approaches.
targeted amplicon sequencing  amplifies known gene targets and attaches a barcode in a second pcr reaction for multiplexing  <cit> .
samples are pooled, sequenced, and analyzed by parsing the sample-specific barcode and then identifying a priori or newly discovered variants  <cit> .
using a targeted amplicon approach, bybee et al.
<cit>  specifically looked at genes useful for phylogenetic analysis.
tas was further extended to a single pcr reaction that utilized linker sequences which allowed common target primers and a single set of barcoded primers to be utilized across distinct samples and loci  <cit> .
complementary to assaying single loci for mas, whole-genome profiling can be utilized for genomic section, qtl mapping, and diversity analysis  <cit> .
whole-genome profiling approaches focus on assaying large numbers of markers while reducing the per sample cost  <cit> .
two common whole-genome profiling methods are snp arrays and genotyping-by-sequencing .
snp arrays are comprised of a large number of known polymorphisms that allow an individual to be genotyped at all sites simultaneously which reduces the overall cost per data point  <cit> .
snp arrays have been used across a range of species to characterize diversity  <cit>  and for association mapping  <cit> .
snp arrays tend to be robust marker platforms but can have limitations, including the inability to target loci that were not included during the array development  and a relatively high per-sample cost.
gbs is a reduced representation whole-genome profiling strategy that leverages rapidly dropping sequencing cost and increasing output.
multiplexing samples with dna barcodes greatly reduces the per sample cost  <cit> .
gbs is one of several reduced representation marker platforms to take advantage of second-generation sequencing platforms that produce enormous amounts of sequence  <cit> .
however, since many samples are sequenced together to minimize cost, the reduced sequencing coverage per sample often results in higher levels of missing data.
since sequencing is only targeted to regions flanking restriction sites, gbs is unable to directly ascertain specific loci, leading to considerable informatics challenges when used in mas.
spiked genotyping-by-sequencing  takes advantage of the abundant sequencing output by combining reduced representation gbs libraries with multiple, targeted amplicons.
sgbs assesses known alleles via targeted amplicon sequencing and individual genotypes are determined by allele frequency counts.
multiple loci can be assayed concurrently since genotyping relies on the direct sequence output.
a similar approach to sgbs was developed by wells et al.
<cit>  that utilizes sequencing-based variant detection by barcoding amplicons.
sgbs is more economical since it uses only a small fraction of available sequencing capacity, the majority of which is simultaneously being used to generate independent, whole-genome profiles.
by combining both approaches, breeders and geneticists are able to employ multi-faceted selection strategies and marker assays with nominal resource expenditure.
to evaluate this approach, we performed sgbs on a winter wheat panel that was screened for six converted kasp markers, four known polymorphisms in the leaf rust resistance gene lr <dig>  and one newly-designed marker targeting a known deletion in lr <dig>
methods
plant material
a panel of  <dig> diverse, advanced wheat lines  was assembled and dna was extracted from seedling leaf tissue using a biosprint  <dig> dna plant kit .
dna was quantified in plates using picogreen and concentrations were normalized to  <dig> ng/μl.
markers
eleven single nucleotide markers were tested for the sgbs approach.
six of the markers were converted from a set of the kasp core markers: bs <dig>  bs <dig>  bs <dig>  bs <dig>  bs <dig>  and bs <dig>  <cit> .
four of the markers were developed from previously designed lr <dig> kasp markers: lr34exon11kasp, lr34exon12kasp, lr34intron4kasp, and lr34exon22kasp  <cit> .
the ‘lr34exon11’ marker from lagudah et al.
<cit>  was also adapted for sgbs, by targeting a  <dig> bp insertion in exon  <dig>  indicative of a non-functional allele .
all primer and allele sequences are provided in additional file 2: table s <dig>  two of the markers from the kasp core collection did not amplify  and were not included in the subsequent analysis.
primer design
primers were designed to amplify the full sequencing construct in a single pcr reaction .
a set of  <dig> unique barcoded primers was developed for multiplexing and to differentiate spiked amplicons from gbs reads .
each barcode primer contains a sequencing platform forward priming site, a unique 10-base barcode, and a m <dig> tail sequence .
these were combined with allele-specific primers that also included the m <dig> tail sequence on the forward primer  <cit> .
the allele-specific reverse primer includes both the flanking sequence reverse primer and the sequencer-specific reverse priming site.
incorporating the m <dig> tail design on both the barcoded primer and allele-specific primer enables the utilization of the same set of barcode oligos for any target sequence, amortizing the cost of oligo synthesis for barcodes across many samples.
the alternative of making barcoded allele specific primers for each target locus would be cost-prohibitive.figure  <dig>
primer and amplicon construction.
the first round of pcr uses a forward primer containing the m <dig> sequence to amplify the target region.
the second round of pcr extends from the m <dig> tail and incorporates a unique barcode, leading to a final product containing the sequencer primers, barcode, m <dig> sequence, and polymorphic target.
kasp markers were converted to primers for sgbs by removing the selective base on the end of each forward primer, effectively creating a single, common forward primer for each locus rather than the two allele-specific primers used for kasp genotyping.
integrating the respective m <dig> and reverse ion torrent sequences on the primer pair made the kasp primer sequences compatible with sgbs.
allele-specific amplification
in a  <dig> well plate,  <dig> ng of dna was combined with  <dig> pmol of m <dig> barcode primer .
a master mix consisting of buffer ,  <dig>  μl mgcl <dig> at  <dig> mm ,  <dig>  μl dntp mix at  <dig>  mm for each nucleotide ,  <dig>  pmol forward-tailed primer ,  <dig> pmol reverse primer ,  <dig>  u taq polymerase, and  <dig>  μl h2o were combined with the dna for a total volume of  <dig> μl for each reaction.
plates were pcr-amplified for  <dig> cycles consisting of 95c , 57c , and 72c .
all samples in the plates were pooled and added to the quantified gbs libraries.
library construction and sequencing
two gbs libraries were prepared for ion torrent™  sequencing following the protocol from mascher et al.
<cit> .
libraries were size-selected on a 2% agarose gel between  <dig> and  <dig> bp, quantified using quant-ittm picogreen® , and normalized to  <dig> nm.
after pooling, the amplicon libraries were quantified using picogreen and normalized to  <dig>  nm.
five μl of the pooled amplicons were added to  <dig> μl of each gbs library for a final concentration of 1% .
the libraries were prepared using the ion pi™ template ot <dig>  <dig> kit  and then sequenced on an ion proton™ system using the ion pi™ chip kit v <dig>  the full protocol for library preparation is provided in additional file  <dig> figure  <dig>
library construction flow chart.
gbs libraries are created following standard protocols.
each spiked library amplifies a single target locus.
spiked libraries are pooled, combined with gbs libraries, and sequenced.
sequence data for the amplicon library is parsed using the m <dig> and unique barcode sequence.
data processing
a tassel pipeline designed for illumina sequence data was modified to identify snps from the gbs tags  <cit> .
specifically, tassel was modified to process ion torrent sequencing sites and with variable length sequence reads.
snp genotypes were called according to the approach of poland et al.
<cit>  using a population-based filter.
a tassel-based custom pipeline was written to determine the allele counts at each amplified locus by identifying the presence of both the m <dig> sequence and the target snp alleles.
reads with the m <dig> tail sequence were parsed by barcode and the number of reads at each allele for a given locus was counted by exact matching to one of the target sequences.
genotype calling for allele-specific amplicons
lines with less than 10x read coverage were not included when clustering and calling genotypes.
genotypes were called using k-means clustering and dbscan clustering, both performed in r  <cit> .
for k-means, the relative proportion of reads for each allele were plotted to determine the appropriate number of clusters to use for this input parameter.
dbscan relies on reachability distance to determine the appropriate number of clusters  <cit> .
varying reachability distances were empirically tested to ascertain an appropriate value.
observationally, a reachability distance of  <dig>  ideally grouped all but one locus.
for bs <dig>  the optimal reachability distance was  <dig> .
RESULTS
to test the approach of spiked gbs, we assayed a panel of diverse wheat lines using gbs to create a whole-genome profile and sgbs to target  <dig> known polymorphic sites.
dna was extracted and normalized and gbs libraries were constructed for the ion proton sequencing platform.
the two sequenced gbs libraries contained  <dig> m and  <dig> m reads with a respective mean read length of  <dig> bp and  <dig> bp.
consistent with previous experience with unspiked gbs libraries,  <dig> % and  <dig> % of reads contained a good gbs barcode and a barcode plus enzyme cut site, respectively.
internal alignment-based discovery resulted in the identification of  <dig>  snps with less than 20% missing data.
this is also consistent with previous unspiked gbs libraries  <cit> .
as a proportion of total sequencing output, the spiked amplicons constituted  <dig> % and  <dig> % of each library as determined by a count of m <dig> sequences.
amplicon libraries were individually analyzed to avoid bias due read number differences.
for each locus, the allelic state of each line was determined by counting the number of reads containing both the sample-specific barcode and a given allele.
genotypes were called using k-means clustering in r and dbscan clustering using the fpc package in r  <cit> .
relative read frequency was used to group individuals into one of three classes: a, b, or heterozygous.
k-means requires a parameter specifying the number of expected clusters while dbscan requires the reachability distance  <cit> .
both of these values require individual curation for loci to ensure two  or three  clusters are correctly called.
generally, there were few differences in the results from either method.
for single-copy loci, both methods performed equally and homozygotes and heterozygotes were easily identifiable .
loci with non-zero axis clusters were also easily identified with both methods.
clusters arising from multi-copy loci were often distinct enough to confidently postulate the genotype allelic state .
overall, the level of concordance between the two clustering algorithms was high with  <dig> % of the genotype calls the same between the two methods .
the majority of discordance was due to k-means requiring that all genotypes be classified whereas dbscan did not classify individuals outside of the main clusters.
the dbscan algorithm is therefore likely of more use in polyploid species where a heterozygote may not be as readily identified .
ignoring the individuals that dbscan did not classify, there was 100% agreement between the two methods.figure  <dig>
k-means clustering and dbscan clustering for lr34exon <dig> and bs <dig>  k-means clustering and dbscan clustering were used to cluster genotypes for each individual on relative read frequency of the two snp alleles.
genotypes called within the same group are denoted by color.
unfilled symbols indicate samples that were not classified by the algorithms.
k-means and  dbscan clustering of lr34exon <dig>  lr34exon <dig> locus is a single-copy locus and the two genotypes are easily distinguished by either clustering algorithm.
heterozygotes are characterized by an equal proportion of both alleles.
k-means and  dbscan clustering of bs <dig>  this primer pair presumably amplifies multiple loci in the polyploid wheat genome that can still be distinguished based on relative read frequency.
the three genotypic classes for individual lines are likely aaaaaa, aabbbb, and aaaabb.
the bbbbbb group does not appear to be present, since a null a genotype should fall on the vertical axis with zero reads counts of allele a. dbscan did not classify the unfilled individual, which is potentially heterozygous at one of the loci .
robust conversion of snp markers between different platforms is important for future genotyping applications, but success can vary considerably  <cit> .
in this study, we observed a good level of conversion from the kasp markers.
two attempted primer sets did not result in amplifying the target sequence and further efforts to optimize conditions for these primer sets were not attempted.
for markers that successfully amplified, the average call rate was  <dig> %.
several markers from the kasp core set resulted in non-zero axis read count clusters, likely due to the existence of homologous copies of the target locus.
the percentage of alleles called for each locus and average coverage are reported in table  <dig> table  <dig>
reads/marker
marker
call rate
avg.
depth
marker name, total call rate, and average read depth.
CONCLUSIONS
with sgbs, we have developed a low-cost, flexible platform for whole-genome profiling and targeted, single-locus genotyping.
the open architecture of primer design for the spiked amplicons enables simple inclusion of new or different target loci.
utilizing a unique set of barcodes combined with locus-specific m <dig> tail primers enabled sequencing of amplified targets in parallel with gbs libraries.
while gbs provides a very low-cost approach for whole-genome profiling, it relies on reproducibly sequencing between restriction sites and cannot target a priori selected loci.
targeted amplicons fill this gap by allowing specific loci to be characterized.
however, with the enormous sequencing output from current sequencing platforms, generating a sufficient number of amplicons across an appropriate number of samples to avoid unreasonable sequencing depth and cost is prohibitive.
to minimize cost, we utilized a small fraction of the sequencing run  while generating more than sufficient coverage across all target loci.
any reasonable number of amplicons could likely be combined with a gbs run.
as with any sequencing approach, increasing the number of samples  decreases coverage.
as sequencing output continues to increase, further ‘excess’ capacity can be leveraged in this way.
however, targeted amplicon numbers beyond 10– <dig> are likely to be impractical relative to a fully designed array or whole-genome characterization .
routine implementation of genotyping approaches in large genetic and breeding applications requires simple and robust laboratory pipelines.
in concert with gbs library development, sgbs target amplification is a streamlined procedure affording routine, high-throughput implementation.
the amplicon libraries are generated through a single pcr reaction, collectively normalized, and pooled with a gbs library.
though not attempted here, multiplex pcr reactions for the allele-specific amplification would further simplify the overall protocol.
sgbs was designed for mas in breeding but is also broadly applicable for a large number of other molecular genetics purposes.
many approaches ranging from diversity studies  <cit>  to genetic and association mapping  <cit>  and genomic selection  <cit>  have successfully applied gbs, but the number of genetic markers generated by gbs often exceeds what is needed for genetic studies, such as fine mapping or tilling.
fine mapping for map-based cloning generally requires screening a very large population with flanking markers for the gene of interest.
while gbs is not a suitable marker platform for fine mapping, utilizing the spiked portion of sgbs for these studies would be ideal.
likewise, the targeted amplicons of sgbs could also be used to screen for novel mutations in tilling or eco-tilling populations.
though a priori snps were targeted in the present study, the direct sequencing of targets also enables de novo discovery of novel mutations as in a tilling study.
for plant breeding, sgbs will enable breeders to genotype large collections of germplasm for specific markers by taking advantage of the massive data output of current sequencing platforms.
large numbers of markers are required for genomic selection, but plant breeders are also interested in characterizing important disease or physiological loci in breeding populations.
sgbs provides a low-cost, scalable approach for both requirements and will serve as an important tool as plant breeding continues its use of molecular markers.
since sgbs amplicons are independent of gbs libraries, breeders can generate a whole-genome profile for advanced breeding material while also applying marker-assisted selection to earlier generations.
importantly, the only realized cost for target genotyping using sgbs is a single pcr reaction.
the ability to quickly identify lines containing specific alleles will enhance the capacity and speed of superior cultivar generation in breeding programs.
plant breeding is inherently an exercise in producing and analyzing large amounts of data to discover improved rare and novel variants.
future advancements in plant breeding will fundamentally rely on new technologies being implemented that allow breeders to progress through this process with the most efficient utilization of resources and least disruption to current workflow.
plant breeding programs have historically depended on single-marker germplasm characterization and are beginning to take advantage of whole-genome profiles for genomic selection.
sgbs combines both approaches, eliminating the current necessity of two distinct platforms while leveraging continual advancements in sequencing technology.
this efficient strategy will allow breeders to increase the amount of germplasm and number of loci that are assayed with few changes to workflow and limited expenditure of resources.
developments like sgbs that enable genomics-assisted breeding are crucial to ensuring progress in developing improved plant varieties in the effort to eliminate hunger and poverty across the world.
additional files
additional file 1: table s <dig>  wheat varieties used in this analysis.
additional file 2: table s <dig>  loci, target alleles, and primer sequences used for sgbs.
additional file 3: table s <dig>  barcode sequences and forward oligo sequences.
additional file 4:
full protocol for spiked genotyping-by-sequencing.
abbreviations
masmarker-assisted selection
gsgenomic selection
kaspkompetitive allele specific pcr
tastargeted amplicon sequencing
gbsgenotyping-by-sequencing
competing interests
the authors declare that they have no competing interests.
authors’ contributions
jp and rb conceived and planned the experiments.
sy carried out the experiments.
tr and jp performed data analysis.
tr and jp wrote the manuscript.
all authors read and approved the final manuscript.

 progress in plant breeding focuses on the rapid development of new cultivars with improved attributes. molecular markers allow breeders to characterize specific lines without the need for laborious and time-consuming phenotyping. marker-assisted selection is used in plant breeding to identify the allele present at a specific locus, allowing the breeder to select based on genotype. mas has been used for plant breeding in many crops to identify specific individuals with known genes of interest, primarily to target large-effect, single targets. since each locus is generally genotyped independently, breeders tend to consider per data point costs when utilizing mas within breeding programs. contemporary marker technologies for assaying single targets that are often used with mas include kasp, targeted amplicon sequencing, and snp arrays. kasp is a uniplex, fluorescence-based single nucleotide genotyping technology that utilizes allele-specific oligo extension. kasp markers have been used for breeding, qtl mapping, and are the main genotyping platform for the generation challenge program at cimmyt. the arrival of inexpensive sequencing has led to the development of economical sequence-based genotyping approaches. targeted amplicon sequencing amplifies known gene targets and attaches a barcode in a second pcr reaction for multiplexing. samples are pooled, sequenced, and analyzed by parsing the sample-specific barcode and then identifying a priori or newly discovered variants. using a targeted amplicon approach, bybee et al. specifically looked at genes useful for phylogenetic analysis. tas was further extended to a single pcr reaction that utilized linker sequences which allowed common target primers and a single set of barcoded primers to be utilized across distinct samples and loci. complementary to assaying single loci for mas, whole-genome profiling can be utilized for genomic section, qtl mapping, and diversity analysis. whole-genome profiling approaches focus on assaying large numbers of markers while reducing the per sample cost. two common whole-genome profiling methods are snp arrays and genotyping-by-sequencing. snp arrays are comprised of a large number of known polymorphisms that allow an individual to be genotyped at all sites simultaneously which reduces the overall cost per data point. snp arrays have been used across a range of species to characterize diversity and for association mapping. snp arrays tend to be robust marker platforms but can have limitations, including the inability to target loci that were not included during the array development and a relatively high per-sample cost. gbs is a reduced representation whole-genome profiling strategy that leverages rapidly dropping sequencing cost and increasing output. multiplexing samples with dna barcodes greatly reduces the per sample cost. gbs is one of several reduced representation marker platforms to take advantage of second-generation sequencing platforms that produce enormous amounts of sequence. however, since many samples are sequenced together to minimize cost, the reduced sequencing coverage per sample often results in higher levels of missing data. since sequencing is only targeted to regions flanking restriction sites, gbs is unable to directly ascertain specific loci, leading to considerable informatics challenges when used in mas. spiked genotyping-by-sequencing takes advantage of the abundant sequencing output by combining reduced representation gbs libraries with multiple, targeted amplicons. sgbs assesses known alleles via targeted amplicon sequencing and individual genotypes are determined by allele frequency counts. multiple loci can be assayed concurrently since genotyping relies on the direct sequence output. a similar approach to sgbs was developed by wells et al. that utilizes sequencing-based variant detection by barcoding amplicons. sgbs is more economical since it uses only a small fraction of available sequencing capacity, the majority of which is simultaneously being used to generate independent, whole-genome profiles. by combining both approaches, breeders and geneticists are able to employ multi-faceted selection strategies and marker assays with nominal resource expenditure. to evaluate this approach, we performed sgbs on a winter wheat panel that was screened for six converted kasp markers, four known polymorphisms in the leaf rust resistance gene lr and one newly-designed marker targeting a known deletion in lr methods plant material a panel of diverse, advanced wheat lines was assembled and dna was extracted from seedling leaf tissue using a biosprint dna plant kit. dna was quantified in plates using picogreen and concentrations were normalized to ng/μl. markers eleven single nucleotide markers were tested for the sgbs approach. six of the markers were converted from a set of the kasp core markers: bs bs bs bs bs and bs. four of the markers were developed from previously designed lr kasp markers: lr34exon11kasp, lr34exon12kasp, lr34intron4kasp, and lr34exon22kasp. the ‘lr34exon11’ marker from lagudah et al. was also adapted for sgbs, by targeting a bp insertion in exon indicative of a non-functional allele. all primer and allele sequences are provided in additional file 2: table s two of the markers from the kasp core collection did not amplify and were not included in the subsequent analysis. primer design primers were designed to amplify the full sequencing construct in a single pcr reaction. a set of unique barcoded primers was developed for multiplexing and to differentiate spiked amplicons from gbs reads. each barcode primer contains a sequencing platform forward priming site, a unique 10-base barcode, and a m tail sequence. these were combined with allele-specific primers that also included the m tail sequence on the forward primer. the allele-specific reverse primer includes both the flanking sequence reverse primer and the sequencer-specific reverse priming site. incorporating the m tail design on both the barcoded primer and allele-specific primer enables the utilization of the same set of barcode oligos for any target sequence, amortizing the cost of oligo synthesis for barcodes across many samples. the alternative of making barcoded allele specific primers for each target locus would be cost-prohibitive.figure primer and amplicon construction. the first round of pcr uses a forward primer containing the m sequence to amplify the target region. the second round of pcr extends from the m tail and incorporates a unique barcode, leading to a final product containing the sequencer primers, barcode, m sequence, and polymorphic target. kasp markers were converted to primers for sgbs by removing the selective base on the end of each forward primer, effectively creating a single, common forward primer for each locus rather than the two allele-specific primers used for kasp genotyping. integrating the respective m and reverse ion torrent sequences on the primer pair made the kasp primer sequences compatible with sgbs. allele-specific amplification in a well plate, ng of dna was combined with pmol of m barcode primer. a master mix consisting of buffer, μl mgcl at mm, μl dntp mix at mm for each nucleotide, pmol forward-tailed primer, pmol reverse primer, u taq polymerase, and μl h2o were combined with the dna for a total volume of μl for each reaction. plates were pcr-amplified for cycles consisting of 95c, 57c, and 72c. all samples in the plates were pooled and added to the quantified gbs libraries. library construction and sequencing two gbs libraries were prepared for ion torrent™ sequencing following the protocol from mascher et al.. libraries were size-selected on a 2% agarose gel between and bp, quantified using quant-ittm picogreen®, and normalized to nm. after pooling, the amplicon libraries were quantified using picogreen and normalized to nm. five μl of the pooled amplicons were added to μl of each gbs library for a final concentration of 1%. the libraries were prepared using the ion pi™ template ot kit and then sequenced on an ion proton™ system using the ion pi™ chip kit v the full protocol for library preparation is provided in additional file figure library construction flow chart. gbs libraries are created following standard protocols. each spiked library amplifies a single target locus. spiked libraries are pooled, combined with gbs libraries, and sequenced. sequence data for the amplicon library is parsed using the m and unique barcode sequence. data processing a tassel pipeline designed for illumina sequence data was modified to identify snps from the gbs tags. specifically, tassel was modified to process ion torrent sequencing sites and with variable length sequence reads. snp genotypes were called according to the approach of poland et al. using a population-based filter. a tassel-based custom pipeline was written to determine the allele counts at each amplified locus by identifying the presence of both the m sequence and the target snp alleles. reads with the m tail sequence were parsed by barcode and the number of reads at each allele for a given locus was counted by exact matching to one of the target sequences. genotype calling for allele-specific amplicons lines with less than 10x read coverage were not included when clustering and calling genotypes. genotypes were called using k-means clustering and dbscan clustering, both performed in r. for k-means, the relative proportion of reads for each allele were plotted to determine the appropriate number of clusters to use for this input parameter. dbscan relies on reachability distance to determine the appropriate number of clusters. varying reachability distances were empirically tested to ascertain an appropriate value. observationally, a reachability distance of ideally grouped all but one locus. for bs the optimal reachability distance was.  to test the approach of spiked gbs, we assayed a panel of diverse wheat lines using gbs to create a whole-genome profile and sgbs to target known polymorphic sites. dna was extracted and normalized and gbs libraries were constructed for the ion proton sequencing platform. the two sequenced gbs libraries contained m and m reads with a respective mean read length of bp and bp. consistent with previous experience with unspiked gbs libraries, % and % of reads contained a good gbs barcode and a barcode plus enzyme cut site, respectively. internal alignment-based discovery resulted in the identification of snps with less than 20% missing data. this is also consistent with previous unspiked gbs libraries. as a proportion of total sequencing output, the spiked amplicons constituted % and % of each library as determined by a count of m sequences. amplicon libraries were individually analyzed to avoid bias due read number differences. for each locus, the allelic state of each line was determined by counting the number of reads containing both the sample-specific barcode and a given allele. genotypes were called using k-means clustering in r and dbscan clustering using the fpc package in r. relative read frequency was used to group individuals into one of three classes: a, b, or heterozygous. k-means requires a parameter specifying the number of expected clusters while dbscan requires the reachability distance. both of these values require individual curation for loci to ensure two or three clusters are correctly called. generally, there were few differences in the results from either method. for single-copy loci, both methods performed equally and homozygotes and heterozygotes were easily identifiable. loci with non-zero axis clusters were also easily identified with both methods. clusters arising from multi-copy loci were often distinct enough to confidently postulate the genotype allelic state. overall, the level of concordance between the two clustering algorithms was high with % of the genotype calls the same between the two methods. the majority of discordance was due to k-means requiring that all genotypes be classified whereas dbscan did not classify individuals outside of the main clusters. the dbscan algorithm is therefore likely of more use in polyploid species where a heterozygote may not be as readily identified. ignoring the individuals that dbscan did not classify, there was 100% agreement between the two methods.figure k-means clustering and dbscan clustering for lr34exon and bs k-means clustering and dbscan clustering were used to cluster genotypes for each individual on relative read frequency of the two snp alleles. genotypes called within the same group are denoted by color. unfilled symbols indicate samples that were not classified by the algorithms. k-means and dbscan clustering of lr34exon lr34exon locus is a single-copy locus and the two genotypes are easily distinguished by either clustering algorithm. heterozygotes are characterized by an equal proportion of both alleles. k-means and dbscan clustering of bs this primer pair presumably amplifies multiple loci in the polyploid wheat genome that can still be distinguished based on relative read frequency. the three genotypic classes for individual lines are likely aaaaaa, aabbbb, and aaaabb. the bbbbbb group does not appear to be present, since a null a genotype should fall on the vertical axis with zero reads counts of allele a. dbscan did not classify the unfilled individual, which is potentially heterozygous at one of the loci. robust conversion of snp markers between different platforms is important for future genotyping applications, but success can vary considerably. in this study, we observed a good level of conversion from the kasp markers. two attempted primer sets did not result in amplifying the target sequence and further efforts to optimize conditions for these primer sets were not attempted. for markers that successfully amplified, the average call rate was %. several markers from the kasp core set resulted in non-zero axis read count clusters, likely due to the existence of homologous copies of the target locus. the percentage of alleles called for each locus and average coverage are reported in table  table reads/marker marker call rate avg. depth marker name, total call rate, and average read depth.  with sgbs, we have developed a low-cost, flexible platform for whole-genome profiling and targeted, single-locus genotyping. the open architecture of primer design for the spiked amplicons enables simple inclusion of new or different target loci. utilizing a unique set of barcodes combined with locus-specific m tail primers enabled sequencing of amplified targets in parallel with gbs libraries. while gbs provides a very low-cost approach for whole-genome profiling, it relies on reproducibly sequencing between restriction sites and cannot target a priori selected loci. targeted amplicons fill this gap by allowing specific loci to be characterized. however, with the enormous sequencing output from current sequencing platforms, generating a sufficient number of amplicons across an appropriate number of samples to avoid unreasonable sequencing depth and cost is prohibitive. to minimize cost, we utilized a small fraction of the sequencing run while generating more than sufficient coverage across all target loci. any reasonable number of amplicons could likely be combined with a gbs run. as with any sequencing approach, increasing the number of samples decreases coverage. as sequencing output continues to increase, further ‘excess’ capacity can be leveraged in this way. however, targeted amplicon numbers beyond 10– are likely to be impractical relative to a fully designed array or whole-genome characterization. routine implementation of genotyping approaches in large genetic and breeding applications requires simple and robust laboratory pipelines. in concert with gbs library development, sgbs target amplification is a streamlined procedure affording routine, high-throughput implementation. the amplicon libraries are generated through a single pcr reaction, collectively normalized, and pooled with a gbs library. though not attempted here, multiplex pcr reactions for the allele-specific amplification would further simplify the overall protocol. sgbs was designed for mas in breeding but is also broadly applicable for a large number of other molecular genetics purposes. many approaches ranging from diversity studies to genetic and association mapping and genomic selection have successfully applied gbs, but the number of genetic markers generated by gbs often exceeds what is needed for genetic studies, such as fine mapping or tilling. fine mapping for map-based cloning generally requires screening a very large population with flanking markers for the gene of interest. while gbs is not a suitable marker platform for fine mapping, utilizing the spiked portion of sgbs for these studies would be ideal. likewise, the targeted amplicons of sgbs could also be used to screen for novel mutations in tilling or eco-tilling populations. though a priori snps were targeted in the present study, the direct sequencing of targets also enables de novo discovery of novel mutations as in a tilling study. for plant breeding, sgbs will enable breeders to genotype large collections of germplasm for specific markers by taking advantage of the massive data output of current sequencing platforms. large numbers of markers are required for genomic selection, but plant breeders are also interested in characterizing important disease or physiological loci in breeding populations. sgbs provides a low-cost, scalable approach for both requirements and will serve as an important tool as plant breeding continues its use of molecular markers. since sgbs amplicons are independent of gbs libraries, breeders can generate a whole-genome profile for advanced breeding material while also applying marker-assisted selection to earlier generations. importantly, the only realized cost for target genotyping using sgbs is a single pcr reaction. the ability to quickly identify lines containing specific alleles will enhance the capacity and speed of superior cultivar generation in breeding programs. plant breeding is inherently an exercise in producing and analyzing large amounts of data to discover improved rare and novel variants. future advancements in plant breeding will fundamentally rely on new technologies being implemented that allow breeders to progress through this process with the most efficient utilization of resources and least disruption to current workflow. plant breeding programs have historically depended on single-marker germplasm characterization and are beginning to take advantage of whole-genome profiles for genomic selection. sgbs combines both approaches, eliminating the current necessity of two distinct platforms while leveraging continual advancements in sequencing technology. this efficient strategy will allow breeders to increase the amount of germplasm and number of loci that are assayed with few changes to workflow and limited expenditure of resources. developments like sgbs that enable genomics-assisted breeding are crucial to ensuring progress in developing improved plant varieties in the effort to eliminate hunger and poverty across the world. additional files additional file 1: table s wheat varieties used in this analysis. additional file 2: table s loci, target alleles, and primer sequences used for sgbs. additional file 3: table s barcode sequences and forward oligo sequences. additional file 4: full protocol for spiked genotyping-by-sequencing. abbreviations masmarker-assisted selection gsgenomic selection kaspkompetitive allele specific pcr tastargeted amplicon sequencing gbsgenotyping-by-sequencing competing interests the authors declare that they have no competing interests. authors’ contributions jp and rb conceived and planned the experiments. sy carried out the experiments. tr and jp performed data analysis. tr and jp wrote the manuscript. all authors read and approved the final manuscript. 
specifically in breeding applications, the sgbs approach will allow breeders to obtain a whole-genome profile of important individuals while simultaneously targeting specific genes for a range of selection strategies across the breeding program. this open platform allows multiple, single-target loci to be assayed while simultaneously generating a whole-genome profile. here, we report the application of sgbs on a winter wheat panel that was screened for converted kasp markers and newly-designed markers targeting known polymorphisms in the leaf rust resistance gene lr to minimize the cost of targeted assays, we utilize a small percent of sequencing capacity available in runs of gbs libraries to “spike” amplified targets of a priori alleles tagged with a different set of unique barcodes.
 in plant breeding, there are two primary applications for dna markers in selection: 1) selection of known genes using a single marker assay; and 2) whole-genome profiling and prediction. typically, marker platforms have addressed only one of these objectives.  we have developed spiked genotyping-by-sequencing, which combines targeted amplicon sequencing with reduced representation genotyping-by-sequencing. to minimize the cost of targeted assays, we utilize a small percent of sequencing capacity available in runs of gbs libraries to “spike” amplified targets of a priori alleles tagged with a different set of unique barcodes. this open platform allows multiple, single-target loci to be assayed while simultaneously generating a whole-genome profile. this dual-genotyping approach allows different sets of samples to be evaluated for single markers or whole genome-profiling. here, we report the application of sgbs on a winter wheat panel that was screened for converted kasp markers and newly-designed markers targeting known polymorphisms in the leaf rust resistance gene lr  the flexibility and low-cost of sgbs will enable a range of applications across genetics research. specifically in breeding applications, the sgbs approach will allow breeders to obtain a whole-genome profile of important individuals while simultaneously targeting specific genes for a range of selection strategies across the breeding program. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords plant breedingwheatmarker platformgenotyping-by-sequencingissue-copyright-statement© the author 2015 