Application of Genomics to Molecular Breeding of Wheat and Barley

Enabling Molecular Technologies for Trait Improvement in Wheat

Wheat is the major staple food crop and a source of calories for humans worldwide. A steady increase in the wheat production is essential to meet the demands of an ever-increasing global population and to achieve food security. The large size and structurally intricate genome of polyploid wheat had hindered the genomic analysis. However, with the advent of new genomic technologies such as next generation sequencing has led to genome drafts for bread wheat and its progenitors and has paved the way to design new strategies for crop improvement. Here we provide an overview of the advancements made in wheat genomics together with the available "omics approaches" and bioinformatics resources developed for wheat research. Advances in genomic, transcriptomic, and metabolomic technologies are highlighted as options to circumvent existing bottlenecks in the phenotypic and genomic selection and gene transfer. The contemporary reverse genetics approaches, including the novel genome editing techniques to inform targeted manipulation of a single/multiple genes and strategies for generating marker-free transgenic wheat plants, emphasize potential to revolutionize wheat improvement shortly.

A high-density, SNP-based consensus map of tetraploid wheat as a bridge to integrate durum and bread wheat genomics and breeding

Consensus linkage maps are important tools in crop genomics. We have assembled a high-density tetraploid wheat consensus map by integrating 13 data sets from independent biparental populations involving durum wheat cultivars (Triticum turgidum ssp. durum), cultivated emmer (T. turgidum ssp. dicoccum) and their ancestor (wild emmer, T. turgidum ssp. dicoccoides). The consensus map harboured 30 144 markers (including 26 626 SNPs and 791 SSRs) half of which were present in at least two component maps. The final map spanned 2631 cM of all 14 durum wheat chromosomes and, differently from the individual component maps, all markers fell within the 14 linkage groups. Marker density per genetic distance unit peaked at centromeric regions, likely due to a combination of low recombination rate in the centromeric regions and even gene distribution along the chromosomes. Comparisons with bread wheat indicated fewer regions with recombination suppression, making this consensus map valuable for mapping in the A and B genomes of both durum and bread wheat. Sequence similarity analysis allowed us to relate mapped gene-derived SNPs to chromosome-specific transcripts. Dense patterns of homeologous relationships have been established between the A- and B-genome maps and between nonsyntenic homeologous chromosome regions as well, the latter tracing to ancient translocation events. The gene-based homeologous relationships are valuable to infer the map location of homeologs of target loci/QTLs. Because most SNP and SSR markers were previously mapped in bread wheat, this consensus map will facilitate a more effective integration and exploitation of genes and QTL for wheat breeding purposes.

Transcriptome sequencing of two wild barley (Hordeum spontaneum L.) ecotypes differentially adapted to drought stress reveals ecotype-specific transcripts

BACKGROUND

Wild barley is adapted to highly diverse environments throughout its geographical distribution range. Transcriptome sequencing of differentially adapted wild barley ecotypes from contrasting environments contributes to the identification of genes and genetic variation involved in abiotic stress tolerance and adaptation.

RESULTS

Two differentially adapted wild barley ecotypes from desert (B1K2) and Mediterranean (B1K30) environments were analyzed for drought stress response under controlled conditions. The desert ecotype lost more water under both irrigation and drought, but exhibited higher relative water content (RWC) and better water use efficiency (WUE) than the coastal ecotype. We sequenced normalized cDNA libraries from drought-stressed leaves of both ecotypes with the 454 platform to identify drought-related transcripts. Over half million reads per ecotype were de novo assembled into 20,439 putative unique transcripts (PUTs) for B1K2, 21,494 for B1K30 and 28,720 for the joint assembly. Over 50% of PUTs of each ecotype were not shared with the other ecotype. Furthermore, 16% (3,245) of B1K2 and 17% (3,674) of B1K30 transcripts did not show orthologous sequence hits in the other wild barley ecotype and cultivated barley, and are candidates of ecotype-specific transcripts. Over 800 unique transcripts from each ecotype homologous to over 30 different stress-related genes were identified. We extracted 1,017 high quality SNPs that differentiated the two ecotypes. The genetic distance between the desert ecotype and cultivated barley was 1.9-fold higher than between the Mediterranean ecotype and cultivated barley. Moreover, the desert ecotype harbored a larger proportion of non-synonymous SNPs than the Mediterranean ecotype suggesting different demographic histories of these ecotypes.

CONCLUSIONS

The results indicate a strong physiological and genomic differentiation between the desert and Mediterranean wild barley ecotypes and a closer relationship of the Mediterranean to cultivated barley. A significant number of novel transcripts specific to wild barley were identified. The higher SNP density and larger proportion of SNPs with functional effects in the desert ecotype suggest different demographic histories and effects of natural selection in Mediterranean and desert wild barley. The data are a valuable genomic resource for an improved genome annotation, transcriptome studies of drought adaptation and a source of new genetic markers for future barley improvement.

A consensus framework map of durum wheat (Triticum durum Desf.) suitable for linkage disequilibrium analysis and genome-wide association mapping

BACKGROUND

Durum wheat (Triticum durum Desf.) is a tetraploid cereal grown in the medium to low-precipitation areas of the Mediterranean Basin, North America and South-West Asia. Genomics applications in durum wheat have the potential to boost exploitation of genetic resources and to advance understanding of the genetics of important complex traits (e.g. resilience to environmental and biotic stresses). A dense and accurate consensus map specific for T. durum will greatly facilitate genetic mapping, functional genomics and marker-assisted improvement.

RESULTS

High quality genotypic data from six core recombinant inbred line populations were used to obtain a consensus framework map of 598 simple sequence repeats (SSR) and Diversity Array Technology® (DArT) anchor markers (common across populations). Interpolation of unique markers from 14 maps allowed us to position a total of 2,575 markers in a consensus map of 2,463 cM. The T. durum A and B genomes were covered in their near totality based on the reference SSR hexaploid wheat map. The consensus locus order compared to those of the single component maps showed good correspondence, (average Spearman's rank correlation rho ρ value of 0.96). Differences in marker order and local recombination rate were observed between the durum and hexaploid wheat consensus maps. The consensus map was used to carry out a whole-genome search for genetic differentiation signatures and association to heading date in a panel of 183 accessions adapted to the Mediterranean areas. Linkage disequilibrium was found to decay below the r2 threshold=0.3 within 2.20 cM, on average. Strong molecular differentiations among sub-populations were mapped to 87 chromosome regions. A genome-wide association scan for heading date from 27 field trials in the Mediterranean Basin and in Mexico yielded 50 chromosome regions with evidences of association in multiple environments.

CONCLUSIONS

The consensus map presented here was used as a reference for genetic diversity and mapping analyses in T. durum, providing nearly complete genome coverage and even marker density. Markers previously mapped in hexaploid wheat constitute a strong link between the two species. The consensus map provides the basis for high-density single nucleotide polymorphic (SNP) marker implementation in durum wheat.

An assessment of the biotechnological use of hemoglobin modulation in cereals

Non-symbiotic hemoglobin (nsHb) genes are ubiquitous in plants, but their biological functions have mostly been studied in model plant species rather than in crops. nsHb influences cell signaling and metabolism by modulating the levels of nitric oxide (NO). Class 1 nsHb is upregulated under hypoxia and is involved in various biotic and abiotic stress responses. Ectopic overexpression of nsHb in Arabidopsis thaliana accelerates development, whilst targeted overexpression in seeds can increase seed yield. Such observations suggest that manipulating nsHb could be a valid biotechnological target. We studied the effects of overexpression of class 1 nsHb in the monocotyledonous crop plant barley (Hordeum vulgare cv. Golden Promise). nsHb was shown to be involved in NO metabolism in barley, as ectopic overexpression reduced the amount of NO released during hypoxia. Further, as in Arabidopsis, nsHb overexpression compromised basal resistance toward pathogens in barley. However, unlike Arabidopsis, nsHb ectopic overexpression delayed growth and development in barley, and seed specific overexpression reduced seed yield. Thus, nsHb overexpression in barley does not seem to be an efficient strategy for increasing yield in cereal crops. These findings highlight the necessity for using actual crop plants rather than laboratory model plants when assessing the effects of biotechnological approaches to crop improvement.

Glutathione S-transferase (GST) family in barley: identification of members, enzyme activity, and gene expression pattern

Barley (Hordeum vulgare) is one of the most important cereals in many developing countries where drought stress considerably diminishes agricultural production. Glutathione S-transferases (GSTs EC 2.5.1.18) are multifunctional enzymes which play a crucial role in cellular detoxification and oxidative stress tolerance. In this study, 84 GST genes were identified in barley by a comprehensive in silico approach. Sequence alignment and phylogenetic analysis grouped these HvGST proteins in eight classes. The largest numbers of the HvGST genes (50) were included in the Tau class followed by 21 genes in Phi, five in Zeta, two in DHAR, two in EF1G, two in Lambda, and one each in TCHQD and Theta classes. Phylogenetic analysis of the putative GSTs from Arabidopsis, rice, and barley indicated that major functional diversification within the GST family predated the monocot/dicot divergence. However, intra-specious duplication seems to be common. Expression patterns of five GST genes from Phi and Tau classes were investigated in three barley genotypes (Yusof [drought-tolerant], Moroc9-75 [drought-sensitive], and HS1 [wild ecotype]) under control and drought-stressed conditions, during the vegetative stage. All investigated genes were up-regulated significantly under drought stress and/or showed a higher level of transcripts in the tolerant cultivar. Additionally, GST enzyme activity was superior in Yusof and induced in the extreme-drought-treated leaves, while it was not changed in Moroc9-75 under drought conditions. Moreover, the lowest and highest levels of lipid peroxidation were observed in the Yusof and Moroc9-75 cultivars, respectively. Based on the achieved results, detoxification and antioxidant activity of GSTs might be considered an important factor in the drought tolerance of barley genotypes for further investigations.

Towards positional isolation of three quantitative trait loci conferring resistance to powdery mildew in two Spanish barley landraces

Three quantitative trait loci (QTL) conferring broad spectrum resistance to powdery mildew, caused by the fungus Blumeria graminis f. sp. hordei, were previously identified on chromosomes 7HS, 7HL and 6HL in the Spanish barley landrace-derived lines SBCC097 and SBCC145. In the present work, a genome-wide putative linear gene index of barley (Genome Zipper) and the first draft of the physical, genetic and functional sequence of the barley genome were used to go one step further in the shortening and explicit demarcation on the barley genome of these regions conferring resistance to powdery mildew as well as in the identification of candidate genes. First, a comparative analysis of the target regions to the barley Genome Zippers of chromosomes 7H and 6H allowed the development of 25 new gene-based molecular markers, which slightly better delimit the QTL intervals. These new markers provided the framework for anchoring of genetic and physical maps, figuring out the outline of the barley genome at the target regions in SBCC097 and SBCC145. The outermost flanking markers of QTLs on 7HS, 7HL and 6HL defined a physical area of 4 Mb, 3.7 Mb and 3.2 Mb, respectively. In total, 21, 10 and 16 genes on 7HS, 7HL and 6HL, respectively, could be interpreted as potential candidates to explain the resistance to powdery mildew, as they encode proteins of related functions with respect to the known pathogen defense-related processes. The majority of these were annotated as belonging to the NBS-LRR class or protein kinase family.

Impact of defoliation severity on photosynthesis, carbon metabolism and transport gene expression in perennial ryegrass

Defoliation severity affects grass regrowth. The changes to biological processes affecting regrowth induced by severe defoliation are not fully understood, nor have they been investigated at a molecular level in field-grown plants. Field-grown perennial ryegrass (Lolium perenne L.) plants were defoliated to 20, 40 or 60mm during winter. Throughout regrowth, transcript profiles of 17 genes involved in photosynthesis and carbon metabolism or transport were characterised in stubble and lamina tissue. Although defoliation to 20mm reduced residual lamina area and stubble water-soluble carbohydrate reserves compared with plants defoliated to 40 or 60mm, net herbage regrowth was not reduced. Transcript profiles indicated a potential compensatory mechanism that may have facilitated regrowth. At the one-leaf regrowth stage, plants defoliated to 20mm had greater abundance of photosynthesis-related gene transcripts (rca, rbcS1, rbcS2, fba, fbp and fnr) and 20% greater stubble total nitrogen than plants defoliated to 60mm. A greater capacity for photosynthesis in outer leaf sheaths may be one potential mechanism used by severely defoliated plants to compensate for the reduced residual lamina area; however, this premise requires further investigation.

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

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

Detection of segregation distortion loci in triticale (x Triticosecale Wittmack) based on a high-density DArT marker consensus genetic linkage map

Background

Triticale is adapted to a wide range of abiotic stress conditions, is an important high-quality feed stock and produces similar grain yield but more biomass compared to other crops. Modern genomic approaches aimed at enhancing breeding progress in cereals require high-quality genetic linkage maps. Consensus maps are genetic maps that are created by a joint analysis of the data from several segregating populations and different approaches are available for their construction. The phenomenon that alleles at a locus deviate from the Mendelian expectation has been defined as segregation distortion. The study of segregation distortion is of particular interest in doubled haploid (DH) populations due to the selection pressure exerted on the plants during the process of their establishment.

Results

The final consensus map, constructed out of six segregating populations derived from nine parental lines, incorporated 2555 DArT markers mapped to 2602 loci (1929 unique). The map spanned 2309.9 cM with an average number of 123.9 loci per chromosome and an average marker density of one unique locus every 1.2 cM. The R genome showed the highest marker coverage followed by the B genome and the A genome. In general, locus order was well maintained between the consensus linkage map and the component maps. However, we observed several groups of loci for which the colinearity was slightly uneven. Among the 2602 loci mapped on the consensus map, 886 showed distorted segregation in at least one of the individual mapping populations. In several DH populations derived by androgenesis, we found chromosomes (2B, 3B, 1R, 2R, 4R and 7R) containing regions where markers exhibited a distorted segregation pattern. In addition, we observed evidence for segregation distortion between pairs of loci caused either by a predominance of parental or recombinant genotypes.

Conclusions

We have constructed a reliable, high-density DArT marker consensus genetic linkage map as a basis for genomic approaches in triticale research and breeding, for example for multiple-line cross QTL mapping experiments. The results of our study exemplify the tremendous impact of different DH production techniques on allele frequencies and segregation distortion covering whole chromosomes.

Comprehensive sequence analysis of 24,783 barley full-length cDNAs derived from 12 clone libraries

Full-length cDNA (FLcDNA) libraries consisting of 172,000 clones were constructed from a two-row malting barley cultivar (Hordeum vulgare 'Haruna Nijo') under normal and stressed conditions. After sequencing the clones from both ends and clustering the sequences, a total of 24,783 complete sequences were produced. By removing duplicates between these and publicly available sequences, 22,651 representative sequences were obtained: 17,773 were novel barley FLcDNAs, and 1,699 were barley specific. Highly conserved genes were found in the barley FLcDNA sequences for 721 of 881 rice (Oryza sativa) trait genes with 50% or greater identity. These FLcDNA resources from our Haruna Nijo cDNA libraries and the full-length sequences of representative clones will improve our understanding of the biological functions of genes in barley, which is the cereal crop with the fourth highest production in the world, and will provide a powerful tool for annotating the barley genome sequences that will become available in the near future.

Insertion site-based polymorphism markers open new perspectives for genome saturation and marker-assisted selection in wheat

In wheat, the deployment of marker-assisted selection has long been hampered by the lack of markers compatible with high-throughput cost-effective genotyping techniques. Recently, insertion site-based polymorphism (ISBP) markers have appeared as very powerful new tools for genomics and genetic studies in hexaploid wheat. To demonstrate their possible use in wheat breeding programmes, we assessed their potential to meet the five main requirements for utilization in MAS: flexible and high-throughput detection methods, low quantity and quality of DNA required, low cost per assay, tight link to target loci and high level of polymorphism in breeding material. Toward this aim, we developed a programme, IsbpFinder, for the automated design of ISBP markers and adapted three detection methods (melting curve analysis, SNaPshot Multiplex System and Illumina BeadArray technology) for high throughput and flexible detection of ISBP or ISBP-derived SNP markers. We demonstrate that the high level of polymorphism of the ISBPs combined with cost-effective genotyping methods can be used to efficiently saturate genetic maps, discriminate between elite cultivars, and design tightly linked diagnostic markers for virtually all target loci in the wheat genome. All together, our results suggest that ISBP markers have the potential to lead to a breakthrough in wheat marker-assisted selection.

DArT markers: diversity analyses, genomes comparison, mapping and integration with SSR markers in Triticum monococcum

BACKGROUND

Triticum monococcum (2n = 2x = 14) is an ancient diploid wheat with many useful traits and is used as a model for wheat gene discovery. DArT (Diversity Arrays Technology) employs a hybridisation-based approach to type thousands of genomic loci in parallel. DArT markers were developed for T. monococcum to assess genetic diversity, compare relationships with hexaploid genomes, and construct a genetic linkage map integrating DArT and microsatellite markers.

RESULTS

A DArT array, consisting of 2304 hexaploid wheat, 1536 tetraploid wheat, 1536 T. monococcum as well as 1536 T. boeoticum representative genomic clones, was used to fingerprint 16 T. monococcum accessions of diverse geographical origins. In total, 846 polymorphic DArT markers were identified, of which 317 were of T. monococcum origin, 246 of hexaploid, 157 of tetraploid, and 126 of T. boeoticum genomes. The fingerprinting data indicated that the geographic origin of T. monococcum accessions was partially correlated with their genetic variation. DArT markers could also well distinguish the genetic differences amongst a panel of 23 hexaploid wheat and nine T. monococcum genomes. For the first time, 274 DArT markers were integrated with 82 simple sequence repeat (SSR) and two morphological trait loci in a genetic map spanning 1062.72 cM in T. monococcum. Six chromosomes were represented by single linkage groups, and chromosome 4Am was formed by three linkage groups. The DArT and SSR genetic loci tended to form independent clusters along the chromosomes. Segregation distortion was observed for one third of the DArT loci. The Ba (black awn) locus was refined to a 23.2 cM region between the DArT marker locus wPt-2584 and the microsatellite locus Xgwmd33 on 1Am; and the Hl (hairy leaf) locus to a 4.0 cM region between DArT loci 376589 and 469591 on 5Am.

CONCLUSION

DArT is a rapid and efficient approach to develop many new molecular markers for genetic studies in T. monococcum. The constructed genetic linkage map will facilitate localisation and map-based cloning of genes of interest, comparative mapping as well as genome organisation and evolution studies between this ancient diploid species and other crops.

DArT markers: diversity analyses, genomes comparison, mapping and integration with SSR markers in Triticum monococcum

Background

Triticum monococcum (2n = 2x = 14) is an ancient diploid wheat with many useful traits and is used as a model for wheat gene discovery. DArT (Diversity Arrays Technology) employs a hybridisation-based approach to type thousands of genomic loci in parallel. DArT markers were developed for T. monococcum to assess genetic diversity, compare relationships with hexaploid genomes, and construct a genetic linkage map integrating DArT and microsatellite markers.

Results

A DArT array, consisting of 2304 hexaploid wheat, 1536 tetraploid wheat, 1536 T. monococcum as well as 1536 T. boeoticum representative genomic clones, was used to fingerprint 16 T. monococcum accessions of diverse geographical origins. In total, 846 polymorphic DArT markers were identified, of which 317 were of T. monococcum origin, 246 of hexaploid, 157 of tetraploid, and 126 of T. boeoticum genomes. The fingerprinting data indicated that the geographic origin of T. monococcum accessions was partially correlated with their genetic variation. DArT markers could also well distinguish the genetic differences amongst a panel of 23 hexaploid wheat and nine T. monococcum genomes. For the first time, 274 DArT markers were integrated with 82 simple sequence repeat (SSR) and two morphological trait loci in a genetic map spanning 1062.72 cM in T. monococcum. Six chromosomes were represented by single linkage groups, and chromosome 4A^m was formed by three linkage groups. The DArT and SSR genetic loci tended to form independent clusters along the chromosomes. Segregation distortion was observed for one third of the DArT loci. The Ba (black awn) locus was refined to a 23.2 cM region between the DArT marker locus wPt-2584 and the microsatellite locus Xgwmd33 on 1A^m; and the Hl (hairy leaf) locus to a 4.0 cM region between DArT loci 376589 and 469591 on 5A^m.

Conclusion

DArT is a rapid and efficient approach to develop many new molecular markers for genetic studies in T. monococcum. The constructed genetic linkage map will facilitate localisation and map-based cloning of genes of interest, comparative mapping as well as genome organisation and evolution studies between this ancient diploid species and other crops.

High-resolution melting analysis of cDNA-derived PCR amplicons for rapid and cost-effective identification of novel alleles in barley

An original method has been established for the identification of novel alleles of eukaryotic translation initiation factor 4E (eIF4E) gene, which is required for resistance to agronomically important bymoviruses, in barley germplasm. This method involves scanning for sequence variations in cDNA-derived PCR amplicons using High-resolution melting (HRM) followed by direct Sanger sequencing of only those amplicons which were predicted to carry nucleotide changes. HRM is a simple, cost-effective, rapid and high-throughput assay, which so far has only been widely used in clinical pathology for molecular diagnostic of diseases and patient genotyping. Application of HRM allowed significant reduction in the amount of expensive Sanger sequencing required for allele mining in plants. The method described here involved an investigation of total cDNA rather than genomic DNA, thus permitting the analyses of shorter (up to 300-bp) and fewer overlapping amplicons to cover the coding sequence. This strategy further reduced the allele mining costs. The sensitivity and accuracy of HRM for predicting genotypes carrying a wide range of nucleotide polymorphisms in eIF4E approached 100%. Results of the current study are promising and suggest that this method could also potentially be applied to the discovery of superior alleles controlling other important traits in barley as well in other model and crop plant species.

TILLMore, a resource for the discovery of chemically induced mutants in barley

A sodium azide-mutagenized population of barley (cv. 'Morex') was developed and utilized to identify mutants at target genes using the 'targeting induced local lesions in genomes' (TILLING) procedure. Screening for mutations at four agronomically important genes (HvCO1, Rpg1, eIF4E and NR) identified a total of 22 new mutant alleles, equivalent to the extrapolated rate of one mutation every 374 kb. All mutations except one were G/C to A/T transitions and several (approximately 68%) implied a change in protein amino acid sequence and therefore a possible effect on phenotype. The high rate of mutation detected through TILLING is in keeping with the high frequency (32.7%) of variant phenotypes observed amongst the M(3) families. Our results indicate the feasibility of using this resource for both reverse and forward genetics approaches to investigate gene function in barley and related crops.

High-density genetic map of durum wheat x wild emmer wheat based on SSR and DArT markers

A genetic linkage map of tetraploid wheat was constructed based on a cross between durum wheat [Triticum turgidum ssp. durum (Desf.) MacKey] cultivar Langdon and wild emmer wheat [T. turgidum ssp. dicoccoides (Körn.) Thell.] accession G18-16. One hundred and fifty-two single-seed descent derived F(6) recombinant inbred lines (RILs) were analyzed with a total of 690 loci, including 197 microsatellite and 493 DArT markers. Linkage analysis defined 14 linkage groups. Most markers were mapped to the B-genome (60%), with an average of 57 markers per chromosome and the remaining 40% mapped to the A-genome, with an average of 39 markers per chromosome. To construct a stabilized (skeleton) map, markers interfering with map stability were removed. The skeleton map consisted of 307 markers with a total length of 2,317 cM and average distance of 7.5 cM between adjacent markers. The length of individual chromosomes ranged between 112 cM for chromosome 4B to 217 cM for chromosome 3B. A fraction (30.1%) of the markers deviated significantly from the expected Mendelian ratios; clusters of loci showing distorted segregation were found on chromosomes 1A, 1BL, 2BS, 3B, and 4B. DArT markers showed high proportion of clustering, which may be indicative of gene-rich regions. Three hundred and fifty-two new DArT markers were mapped for the first time on the current map. This map provides a useful groundwork for further genetic analyses of important quantitative traits, positional cloning, and marker-assisted selection, as well as for genome comparative genomics and genome organization studies in wheat and other cereals.