BACKGROUND
wild relatives of wheat are the tertiary gene pool of bread wheat and contain many favorable genes for wheat improvement, such as disease resistance, drought tolerance, salt tolerance, winter hardiness and adaptability to poor soil .
for the utilization of these elite genes, chromosome engineering has been used to produce small alien introgressions.
irrespective of the enormous genetic variation in wild germplasm and sophisticated techniques available for alien gene transfer, it is still difficult to efficiently identify introgressed chromatin when alien introgression was used in wheat improvement.
cytological approaches including chromosome banding, genomic in situ hybridization  and fluorescence in situ hybridization , have been extensively applied to identify and characterize introgression lines in wheat, but with limitation of low throughput  <cit> .
molecular markers are high efficient for identifying alien chromatin, while the numbers of molecular markers are very low in particular for the tertiary gene pool species of wheat.
hence, there is urgent need to exploit a system for high-throughput method for developing molecular markers to identify alien introgression and translocation lines.
introns are an attractive source of polymorphism for marker development, because insertions/deletions and base substitutions are more common within intron than within exon sequences  <cit> .
intron length polymorphism has been considered as a convenient and reliable source of informative markers with high interspecies transferability  <cit> .
based on the orthologous gene conservation between rice and wheat  <cit> , a set of pcr-based landmark unique gene  primers were developed, and subsequently the plug primers were shown as being suitable to detect polymorphisms among wheat a, b and d genomes  <cit> .
the plug markers not only can identify homology between wheat and alien chromosomes, thus being useful in marker assisted selection , comparative genomics, alien chromosome tracing, taxonomic studies and genotyping .
genome mapping and sequencing in large and polyploid genomes especially wheat remain daunting tasks.
however, the recent technological advances in flow cytometric sorting makes possible the dissecting of large genomes into individual chromosome, reducing sample complexity and enabling analysis at the subgenomic level  <cit> .
flow cytometric chromosome sorting has been successfully applied in bread and durum wheat  <cit> .
recently, flow-sorting of individual chromosomes has been performed in wild relatives of wheat, such as ae.
geniculata, ae.
umbellulata, ae.
comosa, ae.
markgrafii, ae.
triuncialis, ae.cylindrica, t. urartu, ae.
speltoides, ae.
tauschii and ae.
biuncialis .
in combination with the flow sorting and the next-generation sequencing to develop dna sequence-based markers of wild relatives of wheat have already been used to trace alien segments in wheat breeding.
for example, the single-nucleotide polymorphism  markers specific for the short arm of chromosome 5mgs of aegilops geniculata were successfully developed by chromosome sorting and next generation sequencing platform and the snps identified could be employed to accurately detect 5mgs introgressions in common wheat  <cit> .
haynaldia villosa  schur  is the tertiary gene pool of wheat, carrying many important genes for wheat improvement, such as resistant to several wheat diseases, including powdery mildew, eyespot, take-all, wheat spindle streak mosaic virus , and tolerant to drought and cold stresses, good tiller ability, and high grain protein content  <cit> .
previously, we identified wheat yellow mosaic resistance gene wss <dig>  <cit>  located on the short arm of chromosome 4 v .
a set of small fragment translocation lines involving 4vs were obtained using ph1b induction system, by which the wss <dig> was mapped to the terminal region of 4vs  <cit> .
here, we report an efficient method for developing chromosome-specific it markers for the alien 4vs by chromosome sorting and next-generation sequencing, in order to utilize the beneficial genes of wild relatives for wheat breeding.
methods
plant materials
t. durum-h.villosa amphiploid , t. aestivum-h. villosa substitution line ds4v, 4vs ditelosomic addition line, t. aestivum-h. villosa translocation lines t4dl · 4vs, t4vl · 4ds, nau <dig>  nau <dig>  nau <dig>  nau <dig> and nau <dig>  <cit> were developed at the cytogenetics institute, nanjing agricultural university .
h. villosa  was obtained from cambridge botanical garden, uk.
common wheat cv.
chinese spring  maintained at cinau was used as control.
chromosome composition of these materials is shown in additional file 1: figure s <dig>
chromosome sorting and dna sequencing
aqueous suspensions of chromosome 4vs of h. villosa were prepared from synchronized meristem root tip cells following vrana et al.
<cit>  and kubaláková et al.
<cit> .
the chromosomes in suspension were stained with 2 μg/ml dapi  and the 4vs telosomes were sorted using a facsvantage se flow cytometer and sorter .
purity in the sorted fractions was determined after fish with probes for gaa microsatellite and psc <dig>  repeat on chromosomes sorted onto microscope slides.
dna of the sorted chromosome arms was purified and amplified by multiple displacement amplification  as described by Šimková et al.
<cit> .
three independent amplification products were combined to reduce amplification bias.
the amplified dna was purified by ethanol precipitation before sequencing.
about 10 μg of mda-amplified dna was used to create the two shotgun dna-seq libraries of 500-700 bp and 700-1300 bp inserted-size.
the libraries were sequenced in a single lane of illumina hiseq  <dig> platform.
the sequence read data were deposited in the  sequence read archive  and is available under accession number srr <dig>  de novo assembly of the illumina paired-end reads was performed using the software hecate  using different k-mer sizes .
the result of the 45-mer run provided the assembly with the best sequence coverage and n <dig> size, and therefore was determined to generate 4vs scaffolds.
sequence resources for primer design
dna sequences used in this work included h. villosa 4vs assembly, common wheat cv.
chinese spring chromosome arms 4al, 4bs and 4ds assemblies and annotated genes on 4ds  <cit>  and annotated genes in the draft genome sequence assembly of ae.
tauschii  <cit> .
homology and alignment analysis
the flowchart of designing 4vs specific markers is shown in fig.
<dig>  in the first step, we choose a set of genes to calculate exon-exon junction sizes in genomic sequences of homoeolog arms 4al, 4bs, 4ds and 4ds.
the set contained genes from the annotated ae.
tauschii chromosome 4d mapped within 60 mb of pseudomolecule and all annotated genes on chinese spring 4ds.
we removed repeats in the two gene sets by blastn analysis using the cut-off parameters e-value > 1e- <dig>  coverage >80%.
in the second step, all genes were compared with genomic sequences of chinese spring 4al, 4bs, 4ds and h. villosa 4vs using a local blastn program.
all genes matching 4al, 4bs, 4ds and 4vs assemblies and possessing at least one predicted exon-exon junction were selected.
intron sizes of corresponding genes were then calculated and compared against each other.
genes whose intron size in 4vs differed at least 10% from that of 4al, 4bs and 4ds simultaneously in the same homologues allele were chosen for designing the markers.
the primers were designed in the exons flanking the targeted introns.fig.
<dig> schematic representation of the development system of pcr-based it markers specific for the chromosome 4vs of haynaldia villosa and 8% non-denaturing poly-acrylamide gels electrophoresis of pcr product of it marker cinau <dig>  a schematic representation of the development system of pcr-based it markers specific for the chromosome 4vs of haynaldia villosa.
b 8% non-denaturing poly-acrylamide gels electrophoresis of pcr product of it marker cinau <dig>  grey boxes represented the exons and grey lines represented the introns.
red line, light blue line, deep blue line and light green line showed the intron  <dig> region of gene fragment aegta <dig> corresponding to the intron regions of the scaffold of iwgsc_css_4al_scaff_ <dig>  iwgsc_css_4bs_scaff_ <dig>  iwgsc_css_4ds_scaff_ <dig> and hecate_scf <dig>  respectively).
the base number of intron  <dig> in the gene fragment of subgenome aa, bb, dd and vv was 658 bp, 367 bp, 415 bp and 468 bp, respectively
primer design
it primers were designed in exon regions flanking a targeted intron using the online software primer  <dig>  according to the h. villosa sequences, assuming that exon regions and exon-intron structures of orthologous genes are highly conserved among grass species.
the primers designed were expected to amplify 4al, 4bs, 4ds and 4vs genomic dna, simultaneously.
the following criteria were used to design pcr primers for these markers: melting temperature 55–65 °c , primer length ranging between 18– <dig> bases , desired size of amplified fragments estimated was approximately 50 bp more than targeted intron.
all primers were synthesized by shanghai invitrogen biotechnology company ltd. .
dna extraction and pcr
genomic dna was isolated from 2 g fresh leaves of plants at three-leaf stage with sds-phenol-chloroform method according to sharp et al.
<cit>  and devos et al.
<cit>  and purified to eliminate rna, amylase and other unwanted compounds.
the purity and concentration of dna was assessed with microplate reader .
the dna was finally diluted to concentration of 50 ng/μl and stored at −20 °c until use.
pcr amplification was carried out in a 10 μl reaction volume containing 40 ng genomic dna, 2 μm each of the primer pairs,  <dig>  mm each dntps,  <dig>  mm mgcl <dig>  1 × pcr buffer , and  <dig>  u taq dna polymerase with a ptc- <dig> thermal cycler .
the amplification was conducted an initial denaturation at 94 °c for 3 min, followed by 35 cycles at 94 °c for 30 s, annealing of different primers at  <dig>   <dig>  or 60 °c for 50 s at a ramp rate of  <dig>  °c/s, 72 °c for 1 min 10 s, and a final extension at 72 °c for 10 min.
pcr products were resolved in 8% non-denaturing poly-acrylamide gels  and the band patterns were visualized by silver staining  <cit> .
RESULTS
shotguns sequencing of h. villosa chromosome 4vs and its assembly
the analysis of dapi-stained, chromosome suspensions prepared from a wheat-h. villosa 4vs addition line resulted in histograms with five peaks of fluorescence intensity .
the leftmost peak represented telocentric chromosome 4vs, which was well resolved from the chromosome composite peaks i, ii, iii and chromosome 3b peak of the bread wheat .
the flow-sorted telocentric 4vs has more than 89% purity.
dna amplified from flow-sorted t4vs was sequenced by the illumina technology.
in total, we generated a high-quality of  <dig>  gb paired-end reads.
a total length of  <dig>  mb assembled sequences was obtained, comprising  <dig>  scaffolds.
development of pcr-based it markers
a gene set is collected first which might have homologous alleles on 4al, 4bs, 4ds and 4vs, simultaneously.
genes from 4vs are the first choice but poorly annotated.
genes on 4ds either from ae.
tauschii or common wheat, are the second choice as they are presumably most related to 4vs.
ae.
tauschii 4ds was defined as follow.
all annotated genes along the pseudomolecule of ae.
tauschii 4d were aligned by blastn search against the 4vs assembly.
genes within about 60 mb had the highest frequency of blastn hits .
this 60-mb interval is presumed to be ae.
tauschii 4ds which contained  <dig> genes.
a total of  <dig> genes were annotated on 4ds of chinses spring.
if one gene in chinese spring 4ds showed high homologous to that in the ae.
tauschii 4ds by blastn analysis, they are acknowledged as the same .
all genes from ae.
tauschii were reserved, while redundant genes from chinese spring were removed, considering that the average length of genes in ae.
tauschii is longer than that in chinese spring release.
finally, we identified  <dig> genes on chinese spring 4ds were retained.
finally, a gene set totaling  <dig> genes from 4ds were established and aligned to genomic sequences of 4al, 4bs, 4ds and 4vs using blastn.
a total of  <dig> genes have blastn hits on all of four chromosome arms.
out of them,  <dig> genes contained at least one intron and a total of  <dig> introns in the sequences of 4al, 4bs, 4ds and 4vs with an average of  <dig>  introns per gene.
a total of  <dig> introns in 4vs differed by at least 10% as compared to those in wheat 4al, 4bs and 4ds.
they are termed “targeted intron” which meet the criteria for primer design.
these  <dig> introns were assigned to  <dig> genes of which some have two or more introns.
to abide “one genome marker” rule, only one such intron for each gene were selected to develop it marker.
except for  <dig> genes that failed to design primers, a total of  <dig> primer pairs were designed which spanned the targeted intron, using the online software primer 3 v <dig> .
<dig> .
take one gene for designing primers at the targeted intron for an example.
the principle to develop specific it markers of h. villosa 4vs was displayed in fig.
1
.
a gene  fragment from ae.
tauschii was aligned to dna sequences of chromosome arms 4al, 4bs, and 4ds of chinese spring and h. villosa 4vs.
the gene has blastn hits on scaffolds iwgsc_css_4al_scaff_ <dig>  iwgsc_css_4bs_scaff_ <dig>  and iwgsc_css_4ds_scaff_ <dig> on 4al, 4bs and 4ds, and 4vs scaffold hecate_scf <dig>  respectively.
we compared the size of every intron of homologues alleles and found that the length of intron  <dig> differed among a, b, d of wheat and v genome of h. villosa, that they are 685 bp, 367 bp, 415 bp and 468 bp, respectively.
it marker was developed by designing primers on exons  <dig> and  <dig> flanking the intron  <dig> locus.
validation and the efficiency of the it markers development by pcr
to verify the performance of the it marker system, a total of  <dig> it markers were designed using the online software primer 3 v <dig> .
<dig> .
pcr was performed with the designed  <dig> primer sets using genomic dna from the wheat cv.
chinese spring, t. durum-h. villosa amphiploid , h. villosa and t. aestivum-h. villosa translocation line t4dl · 4vs, t. aestivum-h. villosa substitution line ds4v, 4vs ditelosomic addition line and t. aestivum-h. villosa translocation line t4vl · 4ds as template.
only  <dig> out of the  <dig> primer pairs produced no amplicons, while pcr with  <dig> primer pairs resulted in amplicons in all four templates.
in order to assign the  <dig> it markers to chromosome arm 4vs of h. villosa, pcr products obtained in wheat cv.
chinese spring, t. durum-h. villosa amphiploid , h. villosa and t. aestivum-h. villosa translocation line t4dl · 4vs were separated on non-denaturing polyacrylamide gel.
if the amplification of a primer pair generated a distinct pcr product shown as polymorphic band in h. villosa, t. durum-h. villosa amphiploid and t. aestivum-h. villosa translocation line t4dl · 4vs, but not in cv.
chinese spring, the primer pair can be used as chromosome arm 4vs-specific marker.
in total,  <dig> markers were found specific for 4vs .
consequently, the success rate of developing 4vs chromosome arm-specific molecular markers was as high as  <dig> %.
polyacrylamide gel electrophoresis separated pcr products into two, three or four bands.
four bands were obtained with  <dig> primer sets , three bands with the  <dig> primer sets , two bands with the  <dig> primer sets  .
examples of amplification products obtained using six primer pairs are given in fig.
<dig>  the representative  <dig> types it markers corresponding to the sequence information in subgenome aa, bb, dd and vv and the sequences of forward and reverse primers and the size of intron of the representative  <dig> types it markers are given in tables  <dig> and  <dig>  respectively.
the size of the products amplified by these six primer pairs ranged from  <dig> to 1800 bp, and each product was larger than the size predicted from the sequences of wheat and h. villosa intron, suggesting that all pcr products contained introns and parts of exons.fig.
<dig> 8% non-denaturing poly-acrylamide gels electrophoresis of pcr products.
arrows show 4vs-specific fragments
considering the size of pcr products, the four bands amplified by the type i markers originated from chromosomes 4a, 4b, 4d and 4vs, respectively.
these 4vs-specific molecular markers developed are co-dominant and are useful to simultaneously trace the alien 4vs chromosome arm and its wheat homoeologous group.
for the type ii, type iii and type iv markers, according to the size of introns in the subgenomes a, b, d and v, the three bands were amplified from chromosome 4d, 4a and 4 v; 4d, 4b and 4 v; 4a, 4b and 4 v, respectively.
for the type v and vii markers, the two bands detected on polyacrylamide gels consisted of the products derived from 4 vs chromosome arm, and a chinese spring chromosome .
based on pcr product size polymorphisms, a total of  <dig>  it markers could be assigned to subgenomes a, b and d,  <dig> it markers could not be assigned to subgenomes a, b and d. we compared the intron sizes of subgenomes a, b and d and found that the difference in the intron size between homoelogs was small.
discussion
wild relatives of wheat contain a large number of favorable genes for resistance to biotic and abiotic stresses and can be used for wheat improvement.
but now, monitoring alien introgressions in wheat background is difficult because of the shortage of genetic and molecular mapping information on the wild relatives.
dna markers can detect small amounts of alien chromatin that cannot be recognized cytogenetically, thus having high efficiency in the identification of the introgressed fragments.
simple sequence repeat  was used as the molecular markers to screen the individual chromosome of wild relatives of wheat.
however, they were low in transferability to wild relatives and were lack in locus specificity, and therefore had limitation in application to wheat breeding  <cit> .
taken the h. villosa as an example, only  <dig> of the  <dig> wheat microsatellites had high polymorphisms suitable for molecular markers  <cit> .
random amplified polymorphic dnas , although easy to develop, still suffer from poor reproducibility  <cit> .
bin-mapped expressed sequence tags  were also explored as a source of markers, but only a few of them were polymorphic  <cit> .
zhao et al.
<cit>  designed  <dig> primer pairs from wheat est sequences and found that only  <dig>  of the primers amplified specific bands from chromosome 4 v. up to now, there is still fewer molecular markers specific for individual chromosome of wild relatives that can be used, this is partially due to the lack of genome information of the relatives, as reflected by the fact that most of the markers were designed on the basis of genome of wheat, rice, or other species.
flow cytometric chromosome sorting is an efficient method for dramatically simplifying genome analysis by reducing dna sequence complexity  <cit> .
and it combination with next-generation sequencing technology could obtain the sequence composition of single chromosome of wild relatives, thus being highly productive in development and validation of molecular markers  <cit> .
as next-generation sequencing technologies become more economical, making possible the high throughput mining of molecular markers specific for chromosome.
moreover, all wheat chromosomes have been sequenced using flow sorting and next-generation sequencing, which should be useful for the developing molecular markers specific for alien chromosome  <cit> .
tiwari et al.
<cit>  developed a total of  <dig> unique, 5mgs-specific snps of the ae.
geniculata by the combination of chromosome sorting and the next generation sequencing platform, and showed that this approach has high-throughput for the discovery of markers specific for wild relative.
in the present study, we first flow-sorted the chromosome arm 4vs of h. villosa, thereby drastically reducing the amount of work on dna sequencing of this wild relatives.
based on the genome sequence of 4vs, we designed it markers, which are based on the sequence conservation of orthologous genes and therefore have higher transferability between triticeae species.
a total of  <dig> it markers were designed and  <dig>  markers proved to be specific for h. villosa 4vs.
compared to rapd, ssr, and est-pcr markers developed by conventional methods  <cit> , the it markers designing by chromosome sorting and next-generation sequencing is more efficient with high success rate and specificity.
in our lab, zhao et al.
used wheat ph1b mutant to induce translocations of 4vs chromosome fragments to further physically map the wss <dig> to specific chromosome region  <cit> .
however, due to the lower density of markers used in determination of the translocated fragments, a limited resolution of physical maps consisting of  <dig> bins were obtained.
the it markers developed in this study will dramatically increase the density of 4vs physical map or cytological map using these structural aberrants involving 4vs chromosome.
the development of a large number of it markers will be invaluable to trace alien chromatin in a wheat background, for comparative genome mapping, chromosomal evolutionary analysis, and gene introgression during wheat improvement using h. villosa as gene donor.
twenty-one 4vs-specific molecular markers developed in this work are co-dominant and are useful to simultaneously trace the alien 4vs chromosome arm and its wheat homoeologues.
in the present research, we randomly used  <dig> specific it markers for the physical mapping analysis and tracing alien chromatin in wheat background.
using template dna from five wheat-h. villosa 4vs translocation lines  and wheat-h. villosa t4dl · 4vs translocation line, we examined the presence or absence of the  <dig> specific it markers.
thus, the  <dig>   <dig>   <dig>   <dig> and  <dig> it markers were assigned to the region of wheat-h. villosa 4vs translocation lines , respectively.
wheat spindle streak mosaic virus  and wheat yellow mosaic virus  resistance gene wss <dig> were previously mapped to the 4vs arm of h. villosa, using 4 v  substitution and t4dl · 4vs translocation lines  <cit> .
for a more accurate mapping of wss <dig> gene, a ph1b mutant of cv.
chinese spring was used to induce new translocations with smaller 4vs chromosome fragments.
based on the resistance evaluation, gish and molecular marker analysis of the available translocations, the gene conferring the wymv resistance on 4vs were mapped to the distal region of 4vs in the bin of fl  <dig> – <dig>   <cit> .
if the resistance genes are fully explored and used, they would greatly enrich the available resistance germplasm resources for wheat.
wheat-h. villosa 4vs translocation lines  contained wheat yellow mosaic virus  resistance gene, wss <dig>  so these it markers located in the terminal region will be helpful for marker-assisted introgression of the genes of interest, such as wss <dig> gene, into elite cultivars of the common wheat.
miftahudin et al.
<cit>  demonstrated that chromosome 4a have undergone two reciprocal translocations and two inversions events that placed most of the ancestral short arm on the modern long arm .
in our lab, we also used chromowiz to define the 4vs syntenic regions in wheat chromosomes and found that they were enriched on wheat group  <dig> chromosomes 4al.
to develop 4vs-specific it marker, we used dna sequences obtained from flow-sorted chromosome arm 4vs of h. villosa was compared with the sequence assemblies of wheat homoeologous group  <dig> chromosomes  and ae.
tauschii 4ds to identify exon-exon junctions and localize introns on chromosome 4vs of h. villosa.
CONCLUSIONS
in this study, we flow-sorted a ditelosomic addition wheat-h. villosa line to isolate the short arm of the h.villosa 4vs chromosome with ~89% purity identified by fish using a facsvantage se flow cytometer and sorter.
this approach reduced dna sample complexity and permitted the development of markers specific for the short arm of 4vs.
the sequence of h. villosa 4vs was assembled by the software hecate, and then was compared with the sequence assemblies of wheat homoeologous group  <dig> chromosomes and ae.
tauschii 4ds to identify exon-exon junctions and localize introns on chromosome 4vs of h. villosa.
the intron length polymorphisms suitable for designing h. villosa primers were evaluated with criteria, whose intron size of genes in 4vs differed at least 10% from that of 4al, 4bs and 4ds simultaneously in the same homologues allele were chosen for designing the markers.
lastly, we designed a total of  <dig> it markers, among which  <dig>  markers were found specific for 4vs, with the success rate being as high as  <dig> %.
collectively, this approach of combination of chromosome sorting and genome sequencing can be applicable for development of species/genome-specific markers to trace alien chromatin in wheat background, or for comparative genome mapping, chromosomal evolutionary analysis.
additional files
additional file 1: figure s <dig>  the chromosome composition of these materials used in the experiment.
additional file 2: figure s <dig>  histogram of relative dapi fluorescence .
histogram obtained after flow cytometric analysis of mitotic metaphase chromosomes of t. aestivum-h. villosa ditelosomic additional line dt4vs.
peak corresponding to telosomes 4vs  is well discriminated, which facilitated their flow sorting.
sorted chromosome arms were identified after fish with probes for gaa  and psc <dig>   repeat, which results in characteristic banding pattern .
x-axis: relative dapi fluorescence intensity; y-axis: number of particles.
additional file 3: figure s <dig>  the numbers of genes annotated in chinese spring 4ds chromosome and 60-mb of ae.
tauschii 4d.
the venn diagrams displayed the numbers of genes differently annotated in chinese spring 4ds chromosome and 60-mb of ae.
tauschii 4d , and the number of shared conserved genes among the two chromosomes .
additional file 4: table s <dig>  the it markers corresponding to the sequence information in subgenome aa, bb, dd and vv
additional file 5: table s <dig>  the sequences of forward and reverse primers and the size of intron of the it markers
additional file 6: table s <dig>  pcr products of  <dig> it markers.
abbreviations
estsexpressed sequence tags
fishfluorescence in situ hybridization
gishgenomic in situ hybridization
itintron targeting
masmarker assisted selection
plugpcr-based landmark unique gene
rapdrandom amplified polymorphic dnas
snpsingle-nucleotide polymorphism
ssrsimple sequence repeat
wssmvwheat spindle streak mosaic virus
wymvwheat yellow mosaic virus.

 wild relatives of wheat are the tertiary gene pool of bread wheat and contain many favorable genes for wheat improvement, such as disease resistance, drought tolerance, salt tolerance, winter hardiness and adaptability to poor soil. for the utilization of these elite genes, chromosome engineering has been used to produce small alien introgressions. irrespective of the enormous genetic variation in wild germplasm and sophisticated techniques available for alien gene transfer, it is still difficult to efficiently identify introgressed chromatin when alien introgression was used in wheat improvement. cytological approaches including chromosome banding, genomic in situ hybridization and fluorescence in situ hybridization, have been extensively applied to identify and characterize introgression lines in wheat, but with limitation of low throughput. molecular markers are high efficient for identifying alien chromatin, while the numbers of molecular markers are very low in particular for the tertiary gene pool species of wheat. hence, there is urgent need to exploit a system for high-throughput method for developing molecular markers to identify alien introgression and translocation lines. introns are an attractive source of polymorphism for marker development, because insertions/deletions and base substitutions are more common within intron than within exon sequences. intron length polymorphism has been considered as a convenient and reliable source of informative markers with high interspecies transferability. based on the orthologous gene conservation between rice and wheat, a set of pcr-based landmark unique gene primers were developed, and subsequently the plug primers were shown as being suitable to detect polymorphisms among wheat a, b and d genomes. the plug markers not only can identify homology between wheat and alien chromosomes, thus being useful in marker assisted selection, comparative genomics, alien chromosome tracing, taxonomic studies and genotyping. genome mapping and sequencing in large and polyploid genomes especially wheat remain daunting tasks. however, the recent technological advances in flow cytometric sorting makes possible the dissecting of large genomes into individual chromosome, reducing sample complexity and enabling analysis at the subgenomic level. flow cytometric chromosome sorting has been successfully applied in bread and durum wheat. recently, flow-sorting of individual chromosomes has been performed in wild relatives of wheat, such as ae. geniculata, ae. umbellulata, ae. comosa, ae. markgrafii, ae. triuncialis, ae.cylindrica, t. urartu, ae. speltoides, ae. tauschii and ae. biuncialis. in combination with the flow sorting and the next-generation sequencing to develop dna sequence-based markers of wild relatives of wheat have already been used to trace alien segments in wheat breeding. for example, the single-nucleotide polymorphism markers specific for the short arm of chromosome 5mgs of aegilops geniculata were successfully developed by chromosome sorting and next generation sequencing platform and the snps identified could be employed to accurately detect 5mgs introgressions in common wheat. haynaldia villosa schur is the tertiary gene pool of wheat, carrying many important genes for wheat improvement, such as resistant to several wheat diseases, including powdery mildew, eyespot, take-all, wheat spindle streak mosaic virus, and tolerant to drought and cold stresses, good tiller ability, and high grain protein content. previously, we identified wheat yellow mosaic resistance gene wss located on the short arm of chromosome 4 v. a set of small fragment translocation lines involving 4vs were obtained using ph1b induction system, by which the wss was mapped to the terminal region of 4vs. here, we report an efficient method for developing chromosome-specific it markers for the alien 4vs by chromosome sorting and next-generation sequencing, in order to utilize the beneficial genes of wild relatives for wheat breeding. methods plant materials t. durum-h.villosa amphiploid, t. aestivum-h. villosa substitution line ds4v, 4vs ditelosomic addition line, t. aestivum-h. villosa translocation lines t4dl · 4vs, t4vl · 4ds, nau nau nau nau and nau were developed at the cytogenetics institute, nanjing agricultural university. h. villosa was obtained from cambridge botanical garden, uk. common wheat cv. chinese spring maintained at cinau was used as control. chromosome composition of these materials is shown in additional file 1: figure s chromosome sorting and dna sequencing aqueous suspensions of chromosome 4vs of h. villosa were prepared from synchronized meristem root tip cells following vrana et al. and kubaláková et al.. the chromosomes in suspension were stained with 2 μg/ml dapi and the 4vs telosomes were sorted using a facsvantage se flow cytometer and sorter. purity in the sorted fractions was determined after fish with probes for gaa microsatellite and psc repeat on chromosomes sorted onto microscope slides. dna of the sorted chromosome arms was purified and amplified by multiple displacement amplification as described by Šimková et al.. three independent amplification products were combined to reduce amplification bias. the amplified dna was purified by ethanol precipitation before sequencing. about 10 μg of mda-amplified dna was used to create the two shotgun dna-seq libraries of 500-700 bp and 700-1300 bp inserted-size. the libraries were sequenced in a single lane of illumina hiseq platform. the sequence read data were deposited in the sequence read archive and is available under accession number srr de novo assembly of the illumina paired-end reads was performed using the software hecate using different k-mer sizes. the result of the 45-mer run provided the assembly with the best sequence coverage and n size, and therefore was determined to generate 4vs scaffolds. sequence resources for primer design dna sequences used in this work included h. villosa 4vs assembly, common wheat cv. chinese spring chromosome arms 4al, 4bs and 4ds assemblies and annotated genes on 4ds and annotated genes in the draft genome sequence assembly of ae. tauschii. homology and alignment analysis the flowchart of designing 4vs specific markers is shown in fig.  in the first step, we choose a set of genes to calculate exon-exon junction sizes in genomic sequences of homoeolog arms 4al, 4bs, 4ds and 4ds. the set contained genes from the annotated ae. tauschii chromosome 4d mapped within 60 mb of pseudomolecule and all annotated genes on chinese spring 4ds. we removed repeats in the two gene sets by blastn analysis using the cut-off parameters e-value > 1e- coverage >80%. in the second step, all genes were compared with genomic sequences of chinese spring 4al, 4bs, 4ds and h. villosa 4vs using a local blastn program. all genes matching 4al, 4bs, 4ds and 4vs assemblies and possessing at least one predicted exon-exon junction were selected. intron sizes of corresponding genes were then calculated and compared against each other. genes whose intron size in 4vs differed at least 10% from that of 4al, 4bs and 4ds simultaneously in the same homologues allele were chosen for designing the markers. the primers were designed in the exons flanking the targeted introns.fig. schematic representation of the development system of pcr-based it markers specific for the chromosome 4vs of haynaldia villosa and 8% non-denaturing poly-acrylamide gels electrophoresis of pcr product of it marker cinau a schematic representation of the development system of pcr-based it markers specific for the chromosome 4vs of haynaldia villosa. b 8% non-denaturing poly-acrylamide gels electrophoresis of pcr product of it marker cinau grey boxes represented the exons and grey lines represented the introns. red line, light blue line, deep blue line and light green line showed the intron region of gene fragment aegta corresponding to the intron regions of the scaffold of iwgsc_css_4al_scaff_ iwgsc_css_4bs_scaff_ iwgsc_css_4ds_scaff_ and hecate_scf respectively). the base number of intron in the gene fragment of subgenome aa, bb, dd and vv was 658 bp, 367 bp, 415 bp and 468 bp, respectively primer design it primers were designed in exon regions flanking a targeted intron using the online software primer according to the h. villosa sequences, assuming that exon regions and exon-intron structures of orthologous genes are highly conserved among grass species. the primers designed were expected to amplify 4al, 4bs, 4ds and 4vs genomic dna, simultaneously. the following criteria were used to design pcr primers for these markers: melting temperature 55–65 °c, primer length ranging between 18– bases, desired size of amplified fragments estimated was approximately 50 bp more than targeted intron. all primers were synthesized by shanghai invitrogen biotechnology company ltd.. dna extraction and pcr genomic dna was isolated from 2 g fresh leaves of plants at three-leaf stage with sds-phenol-chloroform method according to sharp et al. and devos et al. and purified to eliminate rna, amylase and other unwanted compounds. the purity and concentration of dna was assessed with microplate reader. the dna was finally diluted to concentration of 50 ng/μl and stored at −20 °c until use. pcr amplification was carried out in a 10 μl reaction volume containing 40 ng genomic dna, 2 μm each of the primer pairs,  mm each dntps,  mm mgcl 1 × pcr buffer, and u taq dna polymerase with a ptc- thermal cycler. the amplification was conducted an initial denaturation at 94 °c for 3 min, followed by 35 cycles at 94 °c for 30 s, annealing of different primers at or 60 °c for 50 s at a ramp rate of  °c/s, 72 °c for 1 min 10 s, and a final extension at 72 °c for 10 min. pcr products were resolved in 8% non-denaturing poly-acrylamide gels and the band patterns were visualized by silver staining.  shotguns sequencing of h. villosa chromosome 4vs and its assembly the analysis of dapi-stained, chromosome suspensions prepared from a wheat-h. villosa 4vs addition line resulted in histograms with five peaks of fluorescence intensity. the leftmost peak represented telocentric chromosome 4vs, which was well resolved from the chromosome composite peaks i, ii, iii and chromosome 3b peak of the bread wheat. the flow-sorted telocentric 4vs has more than 89% purity. dna amplified from flow-sorted t4vs was sequenced by the illumina technology. in total, we generated a high-quality of gb paired-end reads. a total length of  mb assembled sequences was obtained, comprising scaffolds. development of pcr-based it markers a gene set is collected first which might have homologous alleles on 4al, 4bs, 4ds and 4vs, simultaneously. genes from 4vs are the first choice but poorly annotated. genes on 4ds either from ae. tauschii or common wheat, are the second choice as they are presumably most related to 4vs. ae. tauschii 4ds was defined as follow. all annotated genes along the pseudomolecule of ae. tauschii 4d were aligned by blastn search against the 4vs assembly. genes within about 60 mb had the highest frequency of blastn hits. this 60-mb interval is presumed to be ae. tauschii 4ds which contained genes. a total of genes were annotated on 4ds of chinses spring. if one gene in chinese spring 4ds showed high homologous to that in the ae. tauschii 4ds by blastn analysis, they are acknowledged as the same. all genes from ae. tauschii were reserved, while redundant genes from chinese spring were removed, considering that the average length of genes in ae. tauschii is longer than that in chinese spring release. finally, we identified genes on chinese spring 4ds were retained. finally, a gene set totaling genes from 4ds were established and aligned to genomic sequences of 4al, 4bs, 4ds and 4vs using blastn. a total of genes have blastn hits on all of four chromosome arms. out of them, genes contained at least one intron and a total of introns in the sequences of 4al, 4bs, 4ds and 4vs with an average of introns per gene. a total of introns in 4vs differed by at least 10% as compared to those in wheat 4al, 4bs and 4ds. they are termed “targeted intron” which meet the criteria for primer design. these introns were assigned to genes of which some have two or more introns. to abide “one genome marker” rule, only one such intron for each gene were selected to develop it marker. except for genes that failed to design primers, a total of primer pairs were designed which spanned the targeted intron, using the online software primer 3 v.. take one gene for designing primers at the targeted intron for an example. the principle to develop specific it markers of h. villosa 4vs was displayed in fig. 1. a gene fragment from ae. tauschii was aligned to dna sequences of chromosome arms 4al, 4bs, and 4ds of chinese spring and h. villosa 4vs. the gene has blastn hits on scaffolds iwgsc_css_4al_scaff_ iwgsc_css_4bs_scaff_ and iwgsc_css_4ds_scaff_ on 4al, 4bs and 4ds, and 4vs scaffold hecate_scf respectively. we compared the size of every intron of homologues alleles and found that the length of intron differed among a, b, d of wheat and v genome of h. villosa, that they are 685 bp, 367 bp, 415 bp and 468 bp, respectively. it marker was developed by designing primers on exons and flanking the intron locus. validation and the efficiency of the it markers development by pcr to verify the performance of the it marker system, a total of it markers were designed using the online software primer 3 v.. pcr was performed with the designed primer sets using genomic dna from the wheat cv. chinese spring, t. durum-h. villosa amphiploid, h. villosa and t. aestivum-h. villosa translocation line t4dl · 4vs, t. aestivum-h. villosa substitution line ds4v, 4vs ditelosomic addition line and t. aestivum-h. villosa translocation line t4vl · 4ds as template. only out of the primer pairs produced no amplicons, while pcr with primer pairs resulted in amplicons in all four templates. in order to assign the it markers to chromosome arm 4vs of h. villosa, pcr products obtained in wheat cv. chinese spring, t. durum-h. villosa amphiploid, h. villosa and t. aestivum-h. villosa translocation line t4dl · 4vs were separated on non-denaturing polyacrylamide gel. if the amplification of a primer pair generated a distinct pcr product shown as polymorphic band in h. villosa, t. durum-h. villosa amphiploid and t. aestivum-h. villosa translocation line t4dl · 4vs, but not in cv. chinese spring, the primer pair can be used as chromosome arm 4vs-specific marker. in total, markers were found specific for 4vs. consequently, the success rate of developing 4vs chromosome arm-specific molecular markers was as high as %. polyacrylamide gel electrophoresis separated pcr products into two, three or four bands. four bands were obtained with primer sets, three bands with the primer sets, two bands with the primer sets. examples of amplification products obtained using six primer pairs are given in fig.  the representative types it markers corresponding to the sequence information in subgenome aa, bb, dd and vv and the sequences of forward and reverse primers and the size of intron of the representative types it markers are given in tables  and respectively. the size of the products amplified by these six primer pairs ranged from to 1800 bp, and each product was larger than the size predicted from the sequences of wheat and h. villosa intron, suggesting that all pcr products contained introns and parts of exons.fig. 8% non-denaturing poly-acrylamide gels electrophoresis of pcr products. arrows show 4vs-specific fragments considering the size of pcr products, the four bands amplified by the type i markers originated from chromosomes 4a, 4b, 4d and 4vs, respectively. these 4vs-specific molecular markers developed are co-dominant and are useful to simultaneously trace the alien 4vs chromosome arm and its wheat homoeologous group. for the type ii, type iii and type iv markers, according to the size of introns in the subgenomes a, b, d and v, the three bands were amplified from chromosome 4d, 4a and 4 v; 4d, 4b and 4 v; 4a, 4b and 4 v, respectively. for the type v and vii markers, the two bands detected on polyacrylamide gels consisted of the products derived from 4 vs chromosome arm, and a chinese spring chromosome. based on pcr product size polymorphisms, a total of it markers could be assigned to subgenomes a, b and d, it markers could not be assigned to subgenomes a, b and d. we compared the intron sizes of subgenomes a, b and d and found that the difference in the intron size between homoelogs was small. discussion wild relatives of wheat contain a large number of favorable genes for resistance to biotic and abiotic stresses and can be used for wheat improvement. but now, monitoring alien introgressions in wheat background is difficult because of the shortage of genetic and molecular mapping information on the wild relatives. dna markers can detect small amounts of alien chromatin that cannot be recognized cytogenetically, thus having high efficiency in the identification of the introgressed fragments. simple sequence repeat was used as the molecular markers to screen the individual chromosome of wild relatives of wheat. however, they were low in transferability to wild relatives and were lack in locus specificity, and therefore had limitation in application to wheat breeding. taken the h. villosa as an example, only of the wheat microsatellites had high polymorphisms suitable for molecular markers. random amplified polymorphic dnas, although easy to develop, still suffer from poor reproducibility. bin-mapped expressed sequence tags were also explored as a source of markers, but only a few of them were polymorphic. zhao et al. designed primer pairs from wheat est sequences and found that only of the primers amplified specific bands from chromosome 4 v. up to now, there is still fewer molecular markers specific for individual chromosome of wild relatives that can be used, this is partially due to the lack of genome information of the relatives, as reflected by the fact that most of the markers were designed on the basis of genome of wheat, rice, or other species. flow cytometric chromosome sorting is an efficient method for dramatically simplifying genome analysis by reducing dna sequence complexity. and it combination with next-generation sequencing technology could obtain the sequence composition of single chromosome of wild relatives, thus being highly productive in development and validation of molecular markers. as next-generation sequencing technologies become more economical, making possible the high throughput mining of molecular markers specific for chromosome. moreover, all wheat chromosomes have been sequenced using flow sorting and next-generation sequencing, which should be useful for the developing molecular markers specific for alien chromosome. tiwari et al. developed a total of unique, 5mgs-specific snps of the ae. geniculata by the combination of chromosome sorting and the next generation sequencing platform, and showed that this approach has high-throughput for the discovery of markers specific for wild relative. in the present study, we first flow-sorted the chromosome arm 4vs of h. villosa, thereby drastically reducing the amount of work on dna sequencing of this wild relatives. based on the genome sequence of 4vs, we designed it markers, which are based on the sequence conservation of orthologous genes and therefore have higher transferability between triticeae species. a total of it markers were designed and markers proved to be specific for h. villosa 4vs. compared to rapd, ssr, and est-pcr markers developed by conventional methods, the it markers designing by chromosome sorting and next-generation sequencing is more efficient with high success rate and specificity. in our lab, zhao et al. used wheat ph1b mutant to induce translocations of 4vs chromosome fragments to further physically map the wss to specific chromosome region. however, due to the lower density of markers used in determination of the translocated fragments, a limited resolution of physical maps consisting of bins were obtained. the it markers developed in this study will dramatically increase the density of 4vs physical map or cytological map using these structural aberrants involving 4vs chromosome. the development of a large number of it markers will be invaluable to trace alien chromatin in a wheat background, for comparative genome mapping, chromosomal evolutionary analysis, and gene introgression during wheat improvement using h. villosa as gene donor. twenty-one 4vs-specific molecular markers developed in this work are co-dominant and are useful to simultaneously trace the alien 4vs chromosome arm and its wheat homoeologues. in the present research, we randomly used specific it markers for the physical mapping analysis and tracing alien chromatin in wheat background. using template dna from five wheat-h. villosa 4vs translocation lines and wheat-h. villosa t4dl · 4vs translocation line, we examined the presence or absence of the specific it markers. thus, the and it markers were assigned to the region of wheat-h. villosa 4vs translocation lines, respectively. wheat spindle streak mosaic virus and wheat yellow mosaic virus resistance gene wss were previously mapped to the 4vs arm of h. villosa, using 4 v substitution and t4dl · 4vs translocation lines. for a more accurate mapping of wss gene, a ph1b mutant of cv. chinese spring was used to induce new translocations with smaller 4vs chromosome fragments. based on the resistance evaluation, gish and molecular marker analysis of the available translocations, the gene conferring the wymv resistance on 4vs were mapped to the distal region of 4vs in the bin of fl –. if the resistance genes are fully explored and used, they would greatly enrich the available resistance germplasm resources for wheat. wheat-h. villosa 4vs translocation lines contained wheat yellow mosaic virus resistance gene, wss so these it markers located in the terminal region will be helpful for marker-assisted introgression of the genes of interest, such as wss gene, into elite cultivars of the common wheat. miftahudin et al. demonstrated that chromosome 4a have undergone two reciprocal translocations and two inversions events that placed most of the ancestral short arm on the modern long arm. in our lab, we also used chromowiz to define the 4vs syntenic regions in wheat chromosomes and found that they were enriched on wheat group chromosomes 4al. to develop 4vs-specific it marker, we used dna sequences obtained from flow-sorted chromosome arm 4vs of h. villosa was compared with the sequence assemblies of wheat homoeologous group chromosomes and ae. tauschii 4ds to identify exon-exon junctions and localize introns on chromosome 4vs of h. villosa.  in this study, we flow-sorted a ditelosomic addition wheat-h. villosa line to isolate the short arm of the h.villosa 4vs chromosome with ~89% purity identified by fish using a facsvantage se flow cytometer and sorter. this approach reduced dna sample complexity and permitted the development of markers specific for the short arm of 4vs. the sequence of h. villosa 4vs was assembled by the software hecate, and then was compared with the sequence assemblies of wheat homoeologous group chromosomes and ae. tauschii 4ds to identify exon-exon junctions and localize introns on chromosome 4vs of h. villosa. the intron length polymorphisms suitable for designing h. villosa primers were evaluated with criteria, whose intron size of genes in 4vs differed at least 10% from that of 4al, 4bs and 4ds simultaneously in the same homologues allele were chosen for designing the markers. lastly, we designed a total of it markers, among which markers were found specific for 4vs, with the success rate being as high as %. collectively, this approach of combination of chromosome sorting and genome sequencing can be applicable for development of species/genome-specific markers to trace alien chromatin in wheat background, or for comparative genome mapping, chromosomal evolutionary analysis. additional files additional file 1: figure s the chromosome composition of these materials used in the experiment. additional file 2: figure s histogram of relative dapi fluorescence. histogram obtained after flow cytometric analysis of mitotic metaphase chromosomes of t. aestivum-h. villosa ditelosomic additional line dt4vs. peak corresponding to telosomes 4vs is well discriminated, which facilitated their flow sorting. sorted chromosome arms were identified after fish with probes for gaa and psc repeat, which results in characteristic banding pattern. x-axis: relative dapi fluorescence intensity; y-axis: number of particles. additional file 3: figure s the numbers of genes annotated in chinese spring 4ds chromosome and 60-mb of ae. tauschii 4d. the venn diagrams displayed the numbers of genes differently annotated in chinese spring 4ds chromosome and 60-mb of ae. tauschii 4d, and the number of shared conserved genes among the two chromosomes. additional file 4: table s the it markers corresponding to the sequence information in subgenome aa, bb, dd and vv additional file 5: table s the sequences of forward and reverse primers and the size of intron of the it markers additional file 6: table s pcr products of it markers. abbreviations estsexpressed sequence tags fishfluorescence in situ hybridization gishgenomic in situ hybridization itintron targeting masmarker assisted selection plugpcr-based landmark unique gene rapdrandom amplified polymorphic dnas snpsingle-nucleotide polymorphism ssrsimple sequence repeat wssmvwheat spindle streak mosaic virus wymvwheat yellow mosaic virus. 
among other, wheat yellow mosaic resistance gene wss and a take-all resistance gene were identified on the short arm of chromosome 4 v of h. villosa. the short arm of chromosome 4vs of h.villosa was flow-sorted using a facsvantage se flow cytometer and sorter, and then sequenced by illumina sequencing. however, monitoring these introgressions has been a daunting task because of inadequate knowledge as to h.villosa genome, as reflected by the lack of specific markers. this study aims to develop 4vs-specific markers by combination of chromosome sorting and next-generation sequencing. the combination of chromosome sorting and next-generation sequencing to develop specific it markers for 4vs of h. villosa has high success rate and specificity, thus being applicable for the development of chromosome-specific markers for alien chromatin in wheat breeding. consequently, we designed a total of intron targeting markers, among which markers were specific for tracing the 4vs chromatin in the wheat background. haynaldia villosa schur is the tertiary gene pool of wheat, and thus a potential resource of genes for wheat improvement. the sequence of h. villosa 4vs was assembled by the software hecate, and then was compared with the sequence assemblies of wheat chromosome arms 4al, 4bs and 4ds and ae. tauschii 4ds, with the objectives of identifying exon-exon junctions and localizing introns on chromosome 4vs of h. villosa. the intron length polymorphisms suitable for designing h. villosa primers were evaluated with criteria.
 haynaldia villosa schur is the tertiary gene pool of wheat, and thus a potential resource of genes for wheat improvement. among other, wheat yellow mosaic resistance gene wss and a take-all resistance gene were identified on the short arm of chromosome 4 v of h. villosa. we had obtained introgressions on 4vs chromosome arm, with the objective of utilizing the target genes. however, monitoring these introgressions has been a daunting task because of inadequate knowledge as to h.villosa genome, as reflected by the lack of specific markers.  this study aims to develop 4vs-specific markers by combination of chromosome sorting and next-generation sequencing. the short arm of chromosome 4vs of h.villosa was flow-sorted using a facsvantage se flow cytometer and sorter, and then sequenced by illumina sequencing. the sequence of h. villosa 4vs was assembled by the software hecate, and then was compared with the sequence assemblies of wheat chromosome arms 4al, 4bs and 4ds and ae. tauschii 4ds, with the objectives of identifying exon-exon junctions and localizing introns on chromosome 4vs of h. villosa. the intron length polymorphisms suitable for designing h. villosa primers were evaluated with criteria. consequently, we designed a total of intron targeting markers, among which markers were specific for tracing the 4vs chromatin in the wheat background.  the combination of chromosome sorting and next-generation sequencing to develop specific it markers for 4vs of h. villosa has high success rate and specificity, thus being applicable for the development of chromosome-specific markers for alien chromatin in wheat breeding. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords triticum aestivumhaynaldia villosamolecular markerintron polymorphism national high technology research program of china2011aa10010201wang xiue the national natural science foundation of china3157165331501305wang haiyan xiao jin the program of introducing talents of discipline to universitiesno.b08025wang xiue the priority academic program development of jiangsu higher education institutions six talent peaks project in jiangsu province the ministry of education, youth and sports of the czech republiclo1204doležel jaroslav ‘948’ project of ministry of agriculture2015-z41wang xiue the fundamental research funds for the central universitieskyz201403wang haiyan the fundamental research funds for the central universitieskj2013003xiao jin the national key research and development program 2016yfd0102001wang xiue issue-copyright-statement© the author 2017 