High throughput SNP discovery and genotyping in hexaploid wheat.

Because of their abundance and their amenability to high-throughput genotyping techniques, Single Nucleotide Polymorphisms (SNPs) are powerful tools for efficient genetics and genomics studies, including characterization of genetic resources, genome-wide association studies and genomic selection. In wheat, most of the previous SNP discovery initiatives targeted the coding fraction, leaving almost 98% of the wheat genome largely unexploited. Here we report on the use of whole-genome resequencing data from eight wheat lines to mine for SNPs in the genic, the repetitive and non-repetitive intergenic fractions of the wheat genome. Eventually, we identified 3.3 million SNPs, 49% being located on the B-genome, 41% on the A-genome and 10% on the D-genome. We also describe the development of the TaBW280K high-throughput genotyping array containing 280,226 SNPs. Performance of this chip was examined by genotyping a set of 96 wheat accessions representing the worldwide diversity. Sixty-nine percent of the SNPs can be efficiently scored, half of them showing a diploid-like clustering. The TaBW280K was proven to be a very efficient tool for diversity analyses, as well as for breeding as it can discriminate between closely related elite varieties. Finally, the TaBW280K array was used to genotype a population derived from a cross between Chinese Spring and Renan, leading to the construction a dense genetic map comprising 83,721 markers. The results described here will provide the wheat community with powerful tools for both basic and applied research.

Small-scale gene duplications played a major role in the recent evolution of wheat chromosome 3B.

BACKGROUND: Bread wheat is not only an important crop, but its large (17 Gb), highly repetitive, and hexaploid genome makes it a good model to study the organization and evolution of complex genomes. Recently, we produced a high quality reference sequence of wheat chromosome 3B (774 Mb), which provides an excellent opportunity to study the evolutionary dynamics of a large and polyploid genome, specifically the impact of single gene duplications. RESULTS: We find that 27 % of the 3B predicted genes are non-syntenic with the orthologous chromosomes of Brachypodium distachyon, Oryza sativa, and Sorghum bicolor, whereas, by applying the same criteria, non-syntenic genes represent on average only 10 % of the predicted genes in these three model grasses. These non-syntenic genes on 3B have high sequence similarity to at least one other gene in the wheat genome, indicating that hexaploid wheat has undergone massive small-scale interchromosomal gene duplications compared to other grasses. Insertions of non-syntenic genes occurred at a similar rate along the chromosome, but these genes tend to be retained at a higher frequency in the distal, recombinogenic regions. The ratio of non-synonymous to synonymous substitution rates showed a more relaxed selection pressure for non-syntenic genes compared to syntenic genes, and gene ontology analysis indicated that non-syntenic genes may be enriched in functions involved in disease resistance. CONCLUSION: Our results highlight the major impact of single gene duplications on the wheat gene complement and confirm the accelerated evolution of the Triticeae lineage among grasses.

Organization and evolution of transposable elements along the bread wheat chromosome 3B.

BACKGROUND: The 17 Gb bread wheat genome has massively expanded through the proliferation of transposable elements (TEs) and two recent rounds of polyploidization. The assembly of a 774 Mb reference sequence of wheat chromosome 3B provided us with the opportunity to explore the impact of TEs on the complex wheat genome structure and evolution at a resolution and scale not reached so far. RESULTS: We develop an automated workflow, CLARI-TE, for TE modeling in complex genomes. We delineate precisely 56,488 intact and 196,391 fragmented TEs along the 3B pseudomolecule, accounting for 85% of the sequence, and reconstruct 30,199 nested insertions. TEs have been mostly silent for the last one million years, and the 3B chromosome has been shaped by a succession of bursts that occurred between 1 to 3 million years ago. Accelerated TE elimination in the high-recombination distal regions is a driving force towards chromosome partitioning. CACTAs overrepresented in the high-recombination distal regions are significantly associated with recently duplicated genes. In addition, we identify 140 CACTA-mediated gene capture events with 17 genes potentially created by exon shuffling and show that 19 captured genes are transcribed and under selection pressure, suggesting the important role of CACTAs in the recent wheat adaptation. CONCLUSION: Accurate TE modeling uncovers the dynamics of TEs in a highly complex and polyploid genome. It provides novel insights into chromosome partitioning and highlights the role of CACTA transposons in the high level of gene duplication in wheat.

QTug.sau-3B is a major quantitative trait locus for wheat hexaploidization.

Meiotic nonreduction resulting in unreduced gametes is thought to be the predominant mechanism underlying allopolyploid formation in plants. Until now, however, its genetic base was largely unknown. The allohexaploid crop common wheat (Triticum aestivum L.), which originated from hybrids of T. turgidum L. with Aegilops tauschii Cosson, provides a model to address this issue. Our observations of meiosis in pollen mother cells from T. turgidum×Ae. tauschii hybrids indicated that first division restitution, which exhibited prolonged cell division during meiosis I, was responsible for unreduced gamete formation. A major quantitative trait locus (QTL) for this trait, named QTug.sau-3B, was detected on chromosome 3B in two T. turgidum×Ae. tauschii haploid populations. This QTL is situated between markers Xgwm285 and Xcfp1012 and covered a genetic distance of 1 cM in one population. QTug.sau-3B is a haploid-dependent QTL because it was not detected in doubled haploid populations. Comparative genome analysis indicated that this QTL was close to Ttam-3B, a collinear homolog of tam in wheat. Although the relationship between QTug.sau-3B and Ttam requires further study, high frequencies of unreduced gametes may be related to reduced expression of Ttam in wheat.

Structural and functional partitioning of bread wheat chromosome 3B.

We produced a reference sequence of the 1-gigabase chromosome 3B of hexaploid bread wheat. By sequencing 8452 bacterial artificial chromosomes in pools, we assembled a sequence of 774 megabases carrying 5326 protein-coding genes, 1938 pseudogenes, and 85% of transposable elements. The distribution of structural and functional features along the chromosome revealed partitioning correlated with meiotic recombination. Comparative analyses indicated high wheat-specific inter- and intrachromosomal gene duplication activities that are potential sources of variability for adaption. In addition to providing a better understanding of the organization, function, and evolution of a large and polyploid genome, the availability of a high-quality sequence anchored to genetic maps will accelerate the identification of genes underlying important agronomic traits.

Meiotic gene evolution: can you teach a new dog new tricks?

Meiosis, the basis of sex, evolved through iterative gene duplications. To understand whether subsequent duplications have further enriched the core meiotic "tool-kit," we investigated the fate of meiotic gene duplicates following whole genome duplication (WGD), a common occurrence in eukaryotes. We show that meiotic genes return to a single copy more rapidly than genome-wide average in angiosperms, one of the lineages in which WGD is most vividly exemplified. The rate at which duplicates are lost decreases through time, a tendency that is also observed genome-wide and may thus prove to be a general trend post-WGD. The sharpest decline is observed for the subset of genes mediating meiotic recombination; however, we found no evidence that the presence of these duplicates is counterselected in two recent polyploid crops selected for fertility. We therefore propose that their loss is passive, highlighting how quickly WGDs are resolved in the absence of selective duplicate retention.

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

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

The physical map of wheat chromosome 1BS provides insights into its gene space organization and evolution.

BACKGROUND: The wheat genome sequence is an essential tool for advanced genomic research and improvements. The generation of a high-quality wheat genome sequence is challenging due to its complex 17 Gb polyploid genome. To overcome these difficulties, sequencing through the construction of BAC-based physical maps of individual chromosomes is employed by the wheat genomics community. Here, we present the construction of the first comprehensive physical map of chromosome 1BS, and illustrate its unique gene space organization and evolution. RESULTS: Fingerprinted BAC clones were assembled into 57 long scaffolds, anchored and ordered with 2,438 markers, covering 83% of chromosome 1BS. The BAC-based chromosome 1BS physical map and gene order of the orthologous regions of model grass species were consistent, providing strong support for the reliability of the chromosome 1BS assembly. The gene space for chromosome 1BS spans the entire length of the chromosome arm, with 76% of the genes organized in small gene islands, accompanied by a two-fold increase in gene density from the centromere to the telomere. CONCLUSIONS: This study provides new evidence on common and chromosome-specific features in the organization and evolution of the wheat genome, including a non-uniform distribution of gene density along the centromere-telomere axis, abundance of non-syntenic genes, the degree of colinearity with other grass genomes and a non-uniform size expansion along the centromere-telomere axis compared with other model cereal genomes. The high-quality physical map constructed in this study provides a solid basis for the assembly of a reference sequence of chromosome 1BS and for breeding applications.

Sequence-based analysis of translocations and inversions in bread wheat (Triticum aestivum L.).

Structural changes of chromosomes are a primary mechanism of genome rearrangement over the course of evolution and detailed knowledge of such changes in a given species and its close relatives should increase the efficiency and precision of chromosome engineering in crop improvement. We have identified sequences bordering each of the main translocation and inversion breakpoints on chromosomes 4A, 5A and 7B of the modern bread wheat genome. The locations of these breakpoints allow, for the first time, a detailed description of the evolutionary origins of these chromosomes at the gene level. Results from this study also demonstrate that, although the strategy of exploiting sorted chromosome arms has dramatically simplified the efforts of wheat genome sequencing, simultaneous analysis of sequences from homoeologous and non-homoeologous chromosomes is essential in understanding the origins of DNA sequences in polyploid species.

Wheat syntenome unveils new evidences of contrasted evolutionary plasticity between paleo- and neoduplicated subgenomes.

Bread wheat derives from a grass ancestor structured in seven protochromosomes followed by a paleotetraploidization to reach a 12 chromosomes intermediate and a neohexaploidization (involving subgenomes A, B and D) event that finally shaped the 21 modern chromosomes. Insights into wheat syntenome in sequencing conserved orthologous set (COS) genes unravelled differences in genomic structure (such as gene conservation and diversity) and genetical landscape (such as recombination pattern) between ancestral as well as recent duplicated blocks. Contrasted evolutionary plasticity is observed where the B subgenome appears more sensitive (i.e. plastic) in contrast to A as dominant (i.e. stable) in response to the neotetraploidization and D subgenome as supra-dominant (i.e. pivotal) in response to the neohexaploidization event. Finally, the wheat syntenome, delivered through a public web interface PlantSyntenyViewer at http://urgi.versailles.inra.fr/synteny-wheat, can be considered as a guide for accelerated dissection of major agronomical traits in wheat.

Wheat centromeric retrotransposons: the new ones take a major role in centromeric structure.

The physical map of the hexaploid wheat chromosome 3B was screened using centromeric DNA probes. A 1.1-Mb region showing the highest number of positive bacterial artificial chromosome (BAC) clones was fully sequenced and annotated, revealing that 96% of the DNA consisted of transposable elements, mainly long terminal repeat (LTR) retrotransposons (88%). Estimation of the insertion time of the transposable elements revealed that CRW (also called Cereba) and Quinta are the youngest elements at the centromeres of common wheat (Triticum spp.) and its diploid ancestors, with Quinta being younger than CRW in both diploid and hexaploid wheats. Chromatin immunoprecipitation experiments revealed that both CRW and Quinta families are targeted by the centromere-specific histone H3 variant CENH3. Immuno colocalization of retroelements and CENH3 antibody indicated that a higher proportion of Quinta than CRWs was associated with CENH3, although CRWs were more abundant. Long arrays of satellite repeats were also identified in the wheat centromere regions, but they lost the ability to bind with CENH3. In addition to transposons, two functional genes and one pseudogene were identified. The gene density in the centromere appeared to be between three and four times lower than the average gene density of chromosome 3B. Comparisons with related grasses also indicated a loss of microcollinearity in this region. Finally, comparison of centromeric sequences of Aegilops tauschii (DD), Triticum boeoticum (AA) and hexaploid wheat revealed that the centromeres in both the polyploids and diploids are still undergoing dynamic changes, and that the new CRWs and Quintas may have undertaken the core role in kinetochore formation.

Intraspecific sequence comparisons reveal similar rates of non-collinear gene insertion in the B and D genomes of bread wheat.

BACKGROUND: Polyploidization is considered one of the main mechanisms of plant genome evolution. The presence of multiple copies of the same gene reduces selection pressure and permits sub-functionalization and neo-functionalization leading to plant diversification, adaptation and speciation. In bread wheat, polyploidization and the prevalence of transposable elements resulted in massive gene duplication and movement. As a result, the number of genes which are non-collinear to genomes of related species seems markedly increased in wheat. RESULTS: We used new-generation sequencing (NGS) to generate sequence of a Mb-sized region from wheat chromosome arm 3DS. Sequence assembly of 24 BAC clones resulted in two scaffolds of 1,264,820 and 333,768 bases. The sequence was annotated and compared to the homoeologous region on wheat chromosome 3B and orthologous loci of Brachypodium distachyon and rice. Among 39 coding sequences in the 3DS scaffolds, 32 have a homoeolog on chromosome 3B. In contrast, only fifteen and fourteen orthologs were identified in the corresponding regions in rice and Brachypodium, respectively. Interestingly, five pseudogenes were identified among the non-collinear coding sequences at the 3B locus, while none was found at the 3DS locus. CONCLUSION: Direct comparison of two Mb-sized regions of the B and D genomes of bread wheat revealed similar rates of non-collinear gene insertion in both genomes with a majority of gene duplications occurring before their divergence. Relatively low proportion of pseudogenes was identified among non-collinear coding sequences. Our data suggest that the pseudogenes did not originate from insertion of non-functional copies, but were formed later during the evolution of hexaploid wheat. Some evidence was found for gene erosion along the B genome locus.

DNA repair and crossing over favor similar chromosome regions as discovered in radiation hybrid of Triticum.

BACKGROUND: The uneven distribution of recombination across the length of chromosomes results in inaccurate estimates of genetic to physical distances. In wheat (Triticum aestivum L.) chromosome 3B, it has been estimated that 90% of the cross over events occur in distal sub-telomeric regions representing 40% of the chromosome. Radiation hybrid (RH) mapping which does not rely on recombination is a strategy to map genomes and has been widely employed in animal species and more recently in some plants. RH maps have been proposed to provide i) higher and ii) more uniform resolution than genetic maps, and iii) to be independent of the distribution patterns observed for meiotic recombination. An in vivo RH panel was generated for mapping chromosome 3B of wheat in an attempt to provide a complete scaffold for this ~1 Gb segment of the genome and compare the resolution to previous genetic maps. RESULTS: A high density RH map with 541 marker loci anchored to chromosome 3B spanning a total distance of 1871.9 cR was generated. Detailed comparisons with a genetic map of similar quality confirmed that i) the overall resolution of the RH map was 10.5 fold higher and ii) six fold more uniform. A significant interaction (r = 0.879 at p = 0.01) was observed between the DNA repair mechanism and the distribution of crossing-over events. This observation could be explained by accepting the possibility that the DNA repair mechanism in somatic cells is affected by the chromatin state in a way similar to the effect that chromatin state has on recombination frequencies in gametic cells. CONCLUSIONS: The RH data presented here support for the first time in vivo the hypothesis of non-casual interaction between recombination hot-spots and DNA repair. Further, two major hypotheses are presented on how chromatin compactness could affect the DNA repair mechanism. Since the initial RH application 37 years ago, we were able to show for the first time that the iii) third hypothesis of RH mapping might not be entirely correct.

A 3,000-loci transcription map of chromosome 3B unravels the structural and functional features of gene islands in hexaploid wheat.

To improve our understanding of the organization and regulation of the wheat (Triticum aestivum) gene space, we established a transcription map of a wheat chromosome (3B) by hybridizing a newly developed wheat expression microarray with bacterial artificial chromosome pools from a new version of the 3B physical map as well as with cDNA probes derived from 15 RNA samples. Mapping data for almost 3,000 genes showed that the gene space spans the whole chromosome 3B with a 2-fold increase of gene density toward the telomeres due to an increase in the number of genes in islands. Comparative analyses with rice (Oryza sativa) and Brachypodium distachyon revealed that these gene islands are composed mainly of genes likely originating from interchromosomal gene duplications. Gene Ontology and expression profile analyses for the 3,000 genes located along the chromosome revealed that the gene islands are enriched significantly in genes sharing the same function or expression profile, thereby suggesting that genes in islands acquired shared regulation during evolution. Only a small fraction of these clusters of cofunctional and coexpressed genes was conserved with rice and B. distachyon, indicating a recent origin. Finally, genes with the same expression profiles in remote islands (coregulation islands) were identified suggesting long-distance regulation of gene expression along the chromosomes in wheat.

Palaeogenomics in cereals: modeling of ancestors for modern species improvement.

During the last decade, technological improvements led to the development of large sets of plant genomic resources permitting the emergence of high-resolution comparative genomic studies. Synteny-based identification of seven shared duplications in cereals led to the modeling of a common ancestral genome structure of 33.6 Mb structured in five protochromosomes containing 9138 protogenes and provided new insights into the evolution of cereal genomes from their extinct ancestors. Recent palaeogenomic data indicate that whole genome duplications were a driving force in the evolutionary success of cereals over the last 50 to 70 millions years. Finally, detailed synteny and duplication relationships led to an improved representation of cereal genomes in concentric circles, thus providing a new reference tool for improved gene annotation and cross-genome markers development.

Fat element-a new marker for chromosome and genome analysis in the Triticeae.

Chromosomal distribution of the Fat element that was isolated from bacterial artificial chromosome (BAC) end sequences of wheat chromosome 3B was studied in 45 species representing eight genera of Poaceae (Aegilops, Triticum, Agropyron, Elymus, Secale, Hordeum, Avena and Triticale) using fluorescence in situ hybridisation (FISH). The Fat sequence was not present in oats and in two barley species, Hordeum vulgare and Hordeum spontaneum, that we investigated. Only very low amounts of the Fat element were detected on the chromosomes of two other barley species, Hordeum geniculatum and Hordeum chilense, with different genome compositions. The chromosomes of other cereal species exhibited distinct hybridisation patterns with the Fat probe, and labelling intensity varied significantly depending on the species or genome. The highest amount of hybridisation was detected on chromosomes of the D genome of Aegilops and Triticum and on chromosomes of the S genome of Agropyron. Despite the bioinformatics analysis of several BAC clones that revealed the tandem organisation of the Fat element, hybridisation with the Fat probe produces uneven, diffuse signals in the proximal regions of chromosomes. In some of the genomes we investigated, however, it also forms distinct, sharp clusters in chromosome-specific positions, and the brightest fluorescence was always observed on group 4 chromosomes. Thus, the Fat element represents a new family of Triticeae-specific, highly repeated DNA elements with a clustered-dispersed distribution pattern. These elements may have first emerged in cereal genomes at the time of divergence of the genus Hordeum from the last common ancestor. During subsequent evolution, the amount and chromosomal distribution of the Fat element changed due to amplification, elimination and re-distribution of this sequence. Because the labelling patterns that we detected were highly specific, the Fat element can be used as an accessory probe in FISH analysis for chromosome identification and investigation of evolutionary processes at the chromosomal level.

Palaeogenomics of plants: synteny-based modelling of extinct ancestors.

In the past ten years, international initiatives have led to the development of large sets of genomic resources that allow comparative genomic studies between plant genomes at a high level of resolution. Comparison of map-based genomic sequences revealed shared intra-genomic duplications, providing new insights into the evolution of flowering plant genomes from common ancestors. Plant genomes can be presented as concentric circles, providing a new reference for plant chromosome evolutionary relationships and an efficient tool for gene annotation and cross-genome markers development. Recent palaeogenomic data demonstrate that whole-genome duplications have provided a motor for the evolutionary success of flowering plants over the last 50-70 million years.

A highly conserved gene island of three genes on chromosome 3B of hexaploid wheat: diverse gene function and genomic structure maintained in a tightly linked block.

BACKGROUND: The complexity of the wheat genome has resulted from waves of retrotransposable element insertions. Gene deletions and disruptions generated by the fast replacement of repetitive elements in wheat have resulted in disruption of colinearity at a micro (sub-megabase) level among the cereals. In view of genomic changes that are possible within a given time span, conservation of genes between species tends to imply an important functional or regional constraint that does not permit a change in genomic structure. The ctg1034 contig completed in this paper was initially studied because it was assigned to the Sr2 resistance locus region, but detailed mapping studies subsequently assigned it to the long arm of 3B and revealed its unusual features. RESULTS: BAC shotgun sequencing of the hexaploid wheat (Triticum aestivum cv. Chinese Spring) genome has been used to assemble a group of 15 wheat BACs from the chromosome 3B physical map FPC contig ctg1034 into a 783,553 bp genomic sequence. This ctg1034 sequence was annotated for biological features such as genes and transposable elements. A three-gene island was identified among >80% repetitive DNA sequence. Using bioinformatics analysis there were no observable similarity in their gene functions. The ctg1034 gene island also displayed complete conservation of gene order and orientation with syntenic gene islands found in publicly available genome sequences of Brachypodium distachyon, Oryza sativa, Sorghum bicolor and Zea mays, even though the intergenic space and introns were divergent. CONCLUSION: We propose that ctg1034 is located within the heterochromatic C-band region of deletion bin 3BL7 based on the identification of heterochromatic tandem repeats and presence of significant matches to chromodomain-containing gypsy LTR retrotransposable elements. We also speculate that this location, among other highly repetitive sequences, may account for the relative stability in gene order and orientation within the gene island.Sequence data from this article have been deposited with the GenBank Data Libraries under accession no. GQ422824.

Two different CC-NBS-LRR genes are required for Lr10-mediated leaf rust resistance in tetraploid and hexaploid wheat.

Comparative study of disease resistance genes in crop plants and their relatives provides insight on resistance gene function, evolution and diversity. Here, we studied the allelic diversity of the Lr10 leaf rust resistance gene, a CC-NBS-LRR coding gene originally isolated from hexaploid wheat, in 20 diploid and tetraploid wheat lines. Besides a gene in the tetraploid wheat variety 'Altar' that is identical to the hexaploid wheat Lr10, two additional, functional resistance alleles showing sequence diversity were identified by virus-induced gene silencing in tetraploid wheat lines. In contrast to most described NBS-LRR proteins, the N-terminal CC domain of LR10 was found to be under strong diversifying selection. A second NBS-LRR gene at the Lr10 locus, RGA2, was shown through silencing to be essential for Lr10 function. Interestingly, RGA2 showed much less sequence diversity than Lr10. These data demonstrate allelic diversity of functional genes at the Lr10 locus in tetraploid wheat, and these new genes can now be analyzed for agronomic relevance. Lr10-based resistance is highly unusual both in its dependence on two, only distantly, related CC-NBS-LRR proteins, as well as in the pattern of diversifying selection in the N-terminal domain. This indicates a new and complex molecular mechanism of pathogen detection and signal transduction.

Fine mapping and marker development for the crossability gene SKr on chromosome 5BS of hexaploid wheat (Triticum aestivum L.).

Most elite wheat varieties cannot be crossed with related species thereby restricting greatly the germplasm that can be used for alien introgression in breeding programs. Inhibition to crossability is controlled genetically and a number of QTL have been identified to date, including the major gene Kr1 on 5BL and SKr, a strong QTL affecting crossability between wheat and rye on chromosome 5BS. In this study, we used a recombinant SSD population originating from a cross between the poorly crossable cultivar Courtot (Ct) and the crossable line MP98 to characterize the major dominant effect of SKr and map the gene at the distal end of the chromosome near the 5B homeologous GSP locus. Colinearity with barley and rice was used to saturate the SKr region with new markers and establish orthologous relationships with a 54-kb region on rice chromosome 12. In total, five markers were mapped within a genetic interval of 0.3 cM and 400 kb of BAC contigs were established on both sides of the gene to lay the foundation for map-based cloning of SKr. Two SSR markers completely linked to SKr were used to evaluate a collection of crossable wheat progenies originating from primary triticale breeding programs. The results confirm the major effect of SKr on crossability and the usefulness of the two markers for the efficient introgression of crossability in elite wheat varieties.