Discovery and development of exome-based, co-dominant single nucleotide polymorphism markers in hexaploid wheat (Triticum aestivum L.)

Saturation Mapping of a Major Effect QTL for Stripe Rust Resistance on Wheat Chromosome 2B in Cultivar Napo 63 Using SNP Genotyping Arrays

Stripe rust or yellow rust (YR), caused by Puccinia striiformis f. sp. tritici (Pst), is one of the most important diseases of wheat (Triticum aestivum L.). Widespread deployment of resistant cultivars is the best means of achieving durable disease control. The red grain, spring wheat cultivar Napo 63 produced by CIMMYT in the 1960s shows a high level of adult-plant resistance to stripe rust in the field. To elucidate the genetic basis of resistance in this cultivar we evaluated 224 F2:3 lines and 175 F2:6 recombinant inbred lines (RILs) derived from a cross between Napo 63 and the Pst-susceptible line Avocet S. The maximum disease severity (MDS) data of F2:3 lines and the relative area under the disease progress curve (rAUDPC) data of RILs were collected during the 2014-2015 and 2015-2016 wheat growing seasons, respectively. Combined bulked segregant analysis and 90K single nucleotide polymorphism (SNP) arrays placed 275 of 511 polymorphic SNPs on chromosome 2B. Sixty four KASP markers selected from the 275 SNPs and 76 SSR markers on 2B were used to identify a chromosome region associated with rust response. A major effect QTL, named Qyrnap.nwafu-2BS, was identified by inclusive composite interval mapping and was preliminarily mapped to a 5.46 cM interval flanked by KASP markers 90K-AN34 and 90K-AN36 in chromosome 2BS. Fourteen KASP markers more closely linked to the locus were developed following a 660K SNP array analysis. The QTL region was finally narrowed to a 0.9 cM interval flanked by KASP markers 660K-AN21 and 660K-AN57 in bin region 2BS-1-0.53. The resistance of Napo 63 was stable across all environments, and as a QTL, explained an average 66.1% of the phenotypic variance in MDS of F2:3 lines and 55.7% of the phenotypic variance in rAUDPC of F5:6 RILs. The short genetic interval and flanking KASP markers developed in the study will facilitate marker-assisted selection, gene pyramiding, and eventual positional cloning of Qyrnap.nwafu-2BS.

Characterization of a Wheat Breeders' Array suitable for high-throughput SNP genotyping of global accessions of hexaploid bread wheat (Triticum aestivum)

Targeted selection and inbreeding have resulted in a lack of genetic diversity in elite hexaploid bread wheat accessions. Reduced diversity can be a limiting factor in the breeding of high yielding varieties and crucially can mean reduced resilience in the face of changing climate and resource pressures. Recent technological advances have enabled the development of molecular markers for use in the assessment and utilization of genetic diversity in hexaploid wheat. Starting with a large collection of 819 571 previously characterized wheat markers, here we describe the identification of 35 143 single nucleotide polymorphism-based markers, which are highly suited to the genotyping of elite hexaploid wheat accessions. To assess their suitability, the markers have been validated using a commercial high-density Affymetrix Axiom^® genotyping array (the Wheat Breeders' Array), in a high-throughput 384 microplate configuration, to characterize a diverse global collection of wheat accessions including landraces and elite lines derived from commercial breeding communities. We demonstrate that the Wheat Breeders' Array is also suitable for generating high-density genetic maps of previously uncharacterized populations and for characterizing novel genetic diversity produced by mutagenesis. To facilitate the use of the array by the wheat community, the markers, the associated sequence and the genotype information have been made available through the interactive web site 'CerealsDB'.

Genotyping-by-sequencing approaches to characterize crop genomes: choosing the right tool for the right application

In the last decade, the revolution in sequencing technologies has deeply impacted crop genotyping practice. New methods allowing rapid, high-throughput genotyping of entire crop populations have proliferated and opened the door to wider use of molecular tools in plant breeding. These new genotyping-by-sequencing (GBS) methods include over a dozen reduced-representation sequencing (RRS) approaches and at least four whole-genome resequencing (WGR) approaches. The diversity of methods available, each often producing different types of data at different cost, can make selection of the best-suited method seem a daunting task. We review the most common genotyping methods used today and compare their suitability for linkage mapping, genomewide association studies (GWAS), marker-assisted and genomic selection and genome assembly and improvement in crops with various genome sizes and complexity. Furthermore, we give an outline of bioinformatics tools for analysis of genotyping data. WGR is well suited to genotyping biparental cross populations with complex, small- to moderate-sized genomes and provides the lowest cost per marker data point. RRS approaches differ in their suitability for various tasks, but demonstrate similar costs per marker data point. These approaches are generally better suited for de novo applications and more cost-effective when genotyping populations with large genomes or high heterozygosity. We expect that although RRS approaches will remain the most cost-effective for some time, WGR will become more widespread for crop genotyping as sequencing costs continue to decrease.

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.

From Genomes to GENE-omes: Exome Sequencing Concept and Applications in Crop Improvement

Exome sequencing represents targeted capture and sequencing of 1-2% of 'high-value genomic regions' (subset of the genome) which are enriched for functional variants and harbors low level of repetitive regions. We discuss here an overview of exome sequencing, ways to approach plant exomes, and advantages and applicability of this powerful approach in deciphering functional regions of genomes. Though initially this approach was developed as an alternative to whole genome sequencing (WGS), but the multitude of benefits conferred by sequence capture via hybridization approaches created a niche for itself to solve many of biological riddles, particularly for resolving phylogenetic distances. The technique has also proved to be successful in understanding the basis of natural and induced molecular variation, marker development and developing genomic resources for complex, wild and non-model species, which are still intractable for WGS efforts. Thus, with profound applications of this powerful sequencing strategy, near future is expected to witness a collective expansion of both techniques, i.e., sequence capture via hybridization for evolutionary and ecological research and WGS approaches for its universal accessibility.

Combined Genomic and Genetic Data Integration of Major Agronomical Traits in Bread Wheat (Triticum aestivum L.)

The high resolution integration of bread wheat genetic and genomic resources accumulated during the last decades offers the opportunity to unveil candidate genes driving major agronomical traits to an unprecedented scale. We combined 27 public quantitative genetic studies and four genetic maps to deliver an exhaustive consensus map consisting of 140,315 molecular markers hosting 221, 73, and 82 Quantitative Trait Loci (QTL) for respectively yield, baking quality, and grain protein content (GPC) related traits. Projection of the consensus genetic map and associated QTLs onto the wheat syntenome made of 99,386 genes ordered on the 21 chromosomes delivered a complete and non-redundant repertoire of 18, 8, 6 metaQTLs for respectively yield, baking quality and GPC, altogether associated to 15,772 genes (delivering 28,630 SNP-based makers) including 37 major candidates. Overall, this study illustrates a translational research approach in transferring information gained from grass relatives to dissect the genomic regions hosting major loci governing key agronomical traits in bread wheat, their flanking markers and associated candidate genes to be now considered as a key resource for breeding programs.

Building strong relationships between conservation genetics and primary industry leads to mutually beneficial genomic advances

Several reviews in the past decade have heralded the benefits of embracing high-throughput sequencing technologies to inform conservation policy and the management of threatened species, but few have offered practical advice on how to expedite the transition from conservation genetics to conservation genomics. Here, we argue that an effective and efficient way to navigate this transition is to capitalize on emerging synergies between conservation genetics and primary industry (e.g., agriculture, fisheries, forestry and horticulture). Here, we demonstrate how building strong relationships between conservation geneticists and primary industry scientists is leading to mutually-beneficial outcomes for both disciplines. Based on our collective experience as collaborative New Zealand-based scientists, we also provide insight for forging these cross-sector relationships.

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

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

CerealsDB 3.0: expansion of resources and data integration

BACKGROUND

The increase in human populations around the world has put pressure on resources, and as a consequence food security has become an important challenge for the 21st century. Wheat (Triticum aestivum) is one of the most important crops in human and livestock diets, and the development of wheat varieties that produce higher yields, combined with increased resistance to pests and resilience to changes in climate, has meant that wheat breeding has become an important focus of scientific research. In an attempt to facilitate these improvements in wheat, plant breeders have employed molecular tools to help them identify genes for important agronomic traits that can be bred into new varieties. Modern molecular techniques have ensured that the rapid and inexpensive characterisation of SNP markers and their validation with modern genotyping methods has produced a valuable resource that can be used in marker assisted selection. CerealsDB was created as a means of quickly disseminating this information to breeders and researchers around the globe.

DESCRIPTION

CerealsDB version 3.0 is an online resource that contains a wide range of genomic datasets for wheat that will assist plant breeders and scientists to select the most appropriate markers for use in marker assisted selection. CerealsDB includes a database which currently contains in excess of a million putative varietal SNPs, of which several hundreds of thousands have been experimentally validated. In addition, CerealsDB also contains new data on functional SNPs predicted to have a major effect on protein function and we have constructed a web service to encourage data integration and high-throughput programmatic access.

CONCLUSION

CerealsDB is an open access website that hosts information on SNPs that are considered useful for both plant breeders and research scientists. The recent inclusion of web services designed to federate genomic data resources allows the information on CerealsDB to be more fully integrated with the WheatIS network and other biological databases.

Sequence capture by hybridization to explore modern and ancient genomic diversity in model and nonmodel organisms

The recent expansion of next-generation sequencing has significantly improved biological research. Nevertheless, deep exploration of genomes or metagenomic samples remains difficult because of the sequencing depth and the associated costs required. Therefore, different partitioning strategies have been developed to sequence informative subsets of studied genomes. Among these strategies, hybridization capture has proven to be an innovative and efficient tool for targeting and enriching specific biomarkers in complex DNA mixtures. It has been successfully applied in numerous areas of biology, such as exome resequencing for the identification of mutations underlying Mendelian or complex diseases and cancers, and its usefulness has been demonstrated in the agronomic field through the linking of genetic variants to agricultural phenotypic traits of interest. Moreover, hybridization capture has provided access to underexplored, but relevant fractions of genomes through its ability to enrich defined targets and their flanking regions. Finally, on the basis of restricted genomic information, this method has also allowed the expansion of knowledge of nonreference species and ancient genomes and provided a better understanding of metagenomic samples. In this review, we present the major advances and discoveries permitted by hybridization capture and highlight the potency of this approach in all areas of biology.

High-density SNP genotyping array for hexaploid wheat and its secondary and tertiary gene pool

In wheat, a lack of genetic diversity between breeding lines has been recognized as a significant block to future yield increases. Species belonging to bread wheat's secondary and tertiary gene pools harbour a much greater level of genetic variability, and are an important source of genes to broaden its genetic base. Introgression of novel genes from progenitors and related species has been widely employed to improve the agronomic characteristics of hexaploid wheat, but this approach has been hampered by a lack of markers that can be used to track introduced chromosome segments. Here, we describe the identification of a large number of single nucleotide polymorphisms that can be used to genotype hexaploid wheat and to identify and track introgressions from a variety of sources. We have validated these markers using an ultra-high-density Axiom(®) genotyping array to characterize a range of diploid, tetraploid and hexaploid wheat accessions and wheat relatives. To facilitate the use of these, both the markers and the associated sequence and genotype information have been made available through an interactive web site.

Exploiting the Repetitive Fraction of the Wheat Genome for High-Throughput Single-Nucleotide Polymorphism Discovery and Genotyping

Transposable elements (TEs) account for more than 80% of the wheat genome. Although they represent a major obstacle for genomic studies, TEs are also a source of polymorphism and consequently of molecular markers such as insertion site-based polymorphism (ISBP) markers. Insertion site-based polymorphisms have been found to be a great source of genome-specific single-nucleotide polymorphism (SNPs) in the hexaploid wheat ( L.) genome. Here, we report on the development of a high-throughput SNP discovery approach based on sequence capture of ISBP markers. By applying this approach to the reference sequence of chromosome 3B from hexaploid wheat, we designed 39,077 SNPs that are evenly distributed along the chromosome. We demonstrate that these SNPs can be efficiently scored with the KASPar (Kompetitive allele-specific polymerase chain reaction) genotyping technology. Finally, through genetic diversity and genome-wide association studies, we also demonstrate that ISBP-derived SNPs can be used in marker-assisted breeding programs.

Genomic-assisted haplotype analysis and the development of high-throughput SNP markers for salinity tolerance in soybean

Soil salinity is a limiting factor of crop yield. The soybean is sensitive to soil salinity, and a dominant gene, Glyma03g32900 is primarily responsible for salt-tolerance. The identification of high throughput and robust markers as well as the deployment of salt-tolerant cultivars are effective approaches to minimize yield loss under saline conditions. We utilized high quality (15x) whole-genome resequencing (WGRS) on 106 diverse soybean lines and identified three major structural variants and allelic variation in the promoter and genic regions of the GmCHX1 gene. The discovery of single nucleotide polymorphisms (SNPs) associated with structural variants facilitated the design of six KASPar assays. Additionally, haplotype analysis and pedigree tracking of 93 U.S. ancestral lines were performed using publically available WGRS datasets. Identified SNP markers were validated, and a strong correlation was observed between the genotype and salt treatment phenotype (leaf scorch, chlorophyll content and Na(+) accumulation) using a panel of 104 soybean lines and, an interspecific bi-parental population (F8) from PI483463 x Hutcheson. These markers precisely identified salt-tolerant/sensitive genotypes (>91%), and different structural-variants (>98%). These SNP assays, supported by accurate phenotyping, haplotype analyses and pedigree tracking information, will accelerate marker-assisted selection programs to enhance the development of salt-tolerant soybean cultivars.

Genomics as the key to unlocking the polyploid potential of wheat

Polyploidy has played a central role in plant genome evolution and in the formation of new species such as tetraploid pasta wheat and hexaploid bread wheat. Until recently, the high sequence conservation between homoeologous genes, together with the large genome size of polyploid wheat, had hindered genomic analyses in this important crop species. In the past 5 yr, however, the advent of next-generation sequencing has radically changed the wheat genomics landscape. Here, we review a series of advances in genomic resources and tools for functional genomics that are shifting the paradigm of what is possible in wheat molecular genetics and breeding. We discuss how understanding the relationship between homoeologues can inform approaches to modulate the response of quantitative traits in polyploid wheat; we also argue that functional redundancy has 'locked up' a wide range of phenotypic variation in wheat. We explore how genomics provides key tools to inform targeted manipulation of multiple homoeologues, thereby allowing researchers and plant breeders to unlock the full polyploid potential of wheat.

Application of Population Sequencing (POPSEQ) for Ordering and Imputing Genotyping-by-Sequencing Markers in Hexaploid Wheat

The advancement of next-generation sequencing technologies in conjunction with new bioinformatics tools enabled fine-tuning of sequence-based, high-resolution mapping strategies for complex genomes. Although genotyping-by-sequencing (GBS) provides a large number of markers, its application for association mapping and genomics-assisted breeding is limited by a large proportion of missing data per marker. For species with a reference genomic sequence, markers can be ordered on the physical map. However, in the absence of reference marker order, the use and imputation of GBS markers is challenging. Here, we demonstrate how the population sequencing (POPSEQ) approach can be used to provide marker context for GBS in wheat. The utility of a POPSEQ-based genetic map as a reference map to create genetically ordered markers on a chromosome for hexaploid wheat was validated by constructing an independent de novo linkage map of GBS markers from a Synthetic W7984 × Opata M85 recombinant inbred line (SynOpRIL) population. The results indicated that there is strong agreement between the independent de novo linkage map and the POPSEQ mapping approach in mapping and ordering GBS markers for hexaploid wheat. After ordering, a large number of GBS markers were imputed, thus providing a high-quality reference map that can be used for QTL mapping for different traits. The POPSEQ-based reference map and whole-genome sequence assemblies are valuable resources that can be used to order GBS markers and enable the application of highly accurate imputation methods to leverage the application GBS markers in wheat.

Physical Mapping of Bread Wheat Chromosome 5A: An Integrated Approach

The huge size, redundancy, and highly repetitive nature of the bread wheat [Triticum aestivum (L.)] genome, makes it among the most difficult species to be sequenced. To overcome these limitations, a strategy based on the separation of individual chromosomes or chromosome arms and the subsequent production of physical maps was established within the frame of the International Wheat Genome Sequence Consortium (IWGSC). A total of 95,812 bacterial artificial chromosome (BAC) clones of short-arm chromosome 5A (5AS) and long-arm chromosome 5A (5AL) arm-specific BAC libraries were fingerprinted and assembled into contigs by complementary analytical approaches based on the FingerPrinted Contig (FPC) and Linear Topological Contig (LTC) tools. Combined anchoring approaches based on polymerase chain reaction (PCR) marker screening, microarray, and sequence homology searches applied to several genomic tools (i.e., genetic maps, deletion bin map, neighbor maps, BAC end sequences (BESs), genome zipper, and chromosome survey sequences) allowed the development of a high-quality physical map with an anchored physical coverage of 75% for 5AS and 53% for 5AL with high portions (64 and 48%, respectively) of contigs ordered along the chromosome. In the genome of grasses, Brachypodium [Brachypodium distachyon (L.) Beauv.], rice (Oryza sativa L.), and sorghum [Sorghum bicolor (L.) Moench] homologs of genes on wheat chromosome 5A were separated into syntenic blocks on different chromosomes as a result of translocations and inversions during evolution. The physical map presented represents an essential resource for fine genetic mapping and map-based cloning of agronomically relevant traits and a reference for the 5A sequencing projects.

Mapping a Type 1 FHB resistance on chromosome 4AS of Triticum macha and deployment in combination with two Type 2 resistances

Markers closely flanking a Type 1 FHB resistance have been produced and the potential of combining this with Type 2 resistances to improve control of FHB has been demonstrated. Two categories of resistance to Fusarium head blight (FHB) in wheat are generally recognised: resistance to initial infection (Type 1) and resistance to spread within the head (Type 2). While numerous sources of Type 2 resistance have been reported, relatively fewer Type 1 resistances have been characterised. Previous study identified a Type 1 FHB resistance (QFhs.jic-4AS) on chromosome 4A in Triticum macha. Little is known about the effect of combining Type 1 and Type 2 resistances on overall FHB symptoms or accumulation of the mycotoxin deoxynivalenol (DON). QFhs.jic-4AS was combined independently with two Type 2 FHB resistances (Fhb1 and one associated with the 1BL/1RS translocation). While combining Type 1 and Type 2 resistances generally reduced visual symptom development, the effect on DON accumulation was marginal. A lack of polymorphic markers and a limited number of recombinants had originally prevented accurate mapping of the QFhs.jic-4AS resistance. Using an array of recently produced markers in combination with new populations, the position of QFhs.jic-4AS has been determined to allow this resistance to be followed in breeding programmes.

The physical map of wheat chromosome 5DS revealed gene duplications and small rearrangements

BACKGROUND

The substantially large bread wheat genome, organized into highly similar three sub-genomes, renders genomic research challenging. The construction of BAC-based physical maps of individual chromosomes reduces the complexity of this allohexaploid genome, enables elucidation of gene space and evolutionary relationships, provides tools for map-based cloning, and serves as a framework for reference sequencing efforts. In this study, we constructed the first comprehensive physical map of wheat chromosome arm 5DS, thereby exploring its gene space organization and evolution.

RESULTS

The physical map of 5DS was comprised of 164 contigs, of which 45 were organized into 21 supercontigs, covering 176 Mb with an N50 value of 2,173 kb. Fifty-eight of the contigs were larger than 1 Mb, with the largest contig spanning 6,649 kb. A total of 1,864 molecular markers were assigned to the map at a density of 10.5 markers/Mb, anchoring 100 of the 120 contigs (>5 clones) that constitute ~95 % of the cumulative length of the map. Ordering of 80 contigs along the deletion bins of chromosome arm 5DS revealed small-scale breaks in syntenic blocks. Analysis of the gene space of 5DS suggested an increasing gradient of genes organized in islands towards the telomere, with the highest gene density of 5.17 genes/Mb in the 0.67-0.78 deletion bin, 1.4 to 1.6 times that of all other bins.

CONCLUSIONS

Here, we provide a chromosome-specific view into the organization and evolution of the D genome of bread wheat, in comparison to one of its ancestors, revealing recent genome rearrangements. The high-quality physical map constructed in this study paves the way for the assembly of a reference sequence, from which breeding efforts will greatly benefit.

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.

RNA-Seq bulked segregant analysis enables the identification of high-resolution genetic markers for breeding in hexaploid wheat

The identification of genetic markers linked to genes of agronomic importance is a major aim of crop research and breeding programmes. Here, we identify markers for Yr15, a major disease resistance gene for wheat yellow rust, using a segregating F2 population. After phenotyping, we implemented RNA sequencing (RNA-Seq) of bulked pools to identify single-nucleotide polymorphisms (SNP) associated with Yr15. Over 27 000 genes with SNPs were identified between the parents, and then classified based on the results from the sequenced bulks. We calculated the bulk frequency ratio (BFR) of SNPs between resistant and susceptible bulks, selecting those showing sixfold enrichment/depletion in the corresponding bulks (BFR > 6). Using additional filtering criteria, we reduced the number of genes with a putative SNP to 175. The 35 SNPs with the highest BFR values were converted into genome-specific KASP assays using an automated bioinformatics pipeline (PolyMarker) which circumvents the limitations associated with the polyploid wheat genome. Twenty-eight assays were polymorphic of which 22 (63%) mapped in the same linkage group as Yr15. Using these markers, we mapped Yr15 to a 0.77-cM interval. The three most closely linked SNPs were tested across varieties and breeding lines representing UK elite germplasm. Two flanking markers were diagnostic in over 99% of lines tested, thus providing a reliable haplotype for marker-assisted selection in these breeding programmes. Our results demonstrate that the proposed methodology can be applied in polyploid F2 populations to generate high-resolution genetic maps across target intervals.

From genome to gene: a new epoch for wheat research?

Genetic research for bread wheat (Triticum aestivum), a staple crop around the world, has been impeded by its complex large hexaploid genome that contains a high proportion of repetitive DNA. Recent advances in sequencing technology have now overcome these challenges and led to genome drafts for bread wheat and its progenitors as well as high-resolution transcriptomes. However, the exploitation of these data for identifying agronomically important genes in wheat is lagging behind. We review recent wheat genome sequencing achievements and focus on four aspects of strategies and future hotspots for wheat improvement: positional cloning, 'omics approaches, combining forward and reverse genetics, and epigenetics.

Single Nucleotide Polymorphism Identification in Polyploids: A Review, Example, and Recommendations

Understanding the relationship between genotype and phenotype is a major biological question and being able to predict phenotypes based on molecular genotypes is integral to molecular breeding. Whole-genome duplications have shaped the history of all flowering plants and present challenges to elucidating the relationship between genotype and phenotype, especially in neopolyploid species. Although single nucleotide polymorphisms (SNPs) have become popular tools for genetic mapping, discovery and application of SNPs in polyploids has been difficult. Here, we summarize common experimental approaches to SNP calling, highlighting recent polyploid successes. To examine the impact of software choice on these analyses, we called SNPs among five peanut genotypes using different alignment programs (BWA-mem and Bowtie 2) and variant callers (SAMtools, GATK, and Freebayes). Alignments produced by Bowtie 2 and BWA-mem and analyzed in SAMtools shared 24.5% concordant SNPs, and SAMtools, GATK, and Freebayes shared 1.4% concordant SNPs. A subsequent analysis of simulated Brassica napus chromosome 1A and 1C genotypes demonstrated that, of the three software programs, SAMtools performed with the highest sensitivity and specificity on Bowtie 2 alignments. These results, however, are likely to vary among species, and we therefore propose a series of best practices for SNP calling in polyploids.

Deregulation of Plant Cell Death Through Disruption of Chloroplast Functionality Affects Asexual Sporulation of Zymoseptoria tritici on Wheat

Chloroplasts have a critical role in plant defense as sites for the biosynthesis of the signaling compounds salicylic acid (SA), jasmonic acid (JA), and nitric oxide (NO) and as major sites of reactive oxygen species production. Chloroplasts, therefore, regarded as important players in the induction and regulation of programmed cell death (PCD) in response to abiotic stresses and pathogen attack. The predominantly foliar pathogen of wheat Zymoseptoria tritici is proposed to exploit the plant PCD, which is associated with the transition in the fungus to the necrotrophic phase of infection. In this study virus-induced gene silencing was used to silence two key genes in carotenoid and chlorophyll biosynthesis, phytoene desaturase (PDS) and Mg-chelatase H subunit (ChlH). The chlorophyll-deficient, PDS- and ChlH-silenced leaves of susceptible plants underwent more rapid pathogen-induced PCD but were significantly less able to support the subsequent asexual sporulation of Z. tritici. Conversely, major gene (Stb6)-mediated resistance to Z. tritici was partially compromised in PDS- and ChlH-silenced leaves. Chlorophyll-deficient wheat ears also displayed increased Z. tritici disease lesion formation accompanied by increased asexual sporulation. These data highlight the importance of chloroplast functionality and its interaction with regulated plant cell death in mediating different genotype and tissue-specific interactions between Z. tritici and wheat.

Identification and characterization of more than 4 million intervarietal SNPs across the group 7 chromosomes of bread wheat

Despite being a major international crop, our understanding of the wheat genome is relatively poor due to its large size and complexity. To gain a greater understanding of wheat genome diversity, we have identified single nucleotide polymorphisms between 16 Australian bread wheat varieties. Whole-genome shotgun Illumina paired read sequence data were mapped to the draft assemblies of chromosomes 7A, 7B and 7D to identify more than 4 million intervarietal SNPs. SNP density varied between the three genomes, with much greater density observed on the A and B genomes than the D genome. This variation may be a result of substantial gene flow from the tetraploid Triticum turgidum, which possesses A and B genomes, during early co-cultivation of tetraploid and hexaploid wheat. In addition, we examined SNP density variation along the chromosome syntenic builds and identified genes in low-density regions which may have been selected during domestication and breeding. This study highlights the impact of evolution and breeding on the bread wheat genome and provides a substantial resource for trait association and crop improvement. All SNP data are publically available on a generic genome browser GBrowse at www.wheatgenome.info.

Applying association mapping and genomic selection to the dissection of key traits in elite European wheat

We show the application of association mapping and genomic selection for key breeding targets using a large panel of elite winter wheat varieties and a large volume of agronomic data. The heightening urgency to increase wheat production in line with the needs of a growing population, and in the face of climatic uncertainty, mean new approaches, including association mapping (AM) and genomic selection (GS) need to be validated and applied in wheat breeding. Key adaptive responses are the cornerstone of regional breeding. There is evidence that new ideotypes for long-standing traits such as flowering time may be required. In order to detect targets for future marker-assisted improvement and validate the practical application of GS for wheat breeding we genotyped 376 elite wheat varieties with 3,046 DArT, single nucleotide polymorphism and gene markers and measured seven traits in replicated yield trials over 2 years in France, Germany and the UK. The scale of the phenotyping exceeds the breadth of previous AM and GS studies in these key economic wheat production regions of Northern Europe. Mixed-linear modelling (MLM) detected significant marker-trait associations across and within regions. Genomic prediction using elastic net gave low to high prediction accuracies depending on the trait, and could be experimentally increased by modifying the constituents of the training population (TP). We also tested the use of differentially penalised regression to integrate candidate gene and genome-wide markers to predict traits, demonstrating the validity and simplicity of this approach. Overall, our results suggest that whilst AM offers potential for application in both research and breeding, GS represents an exciting opportunity to select key traits, and that optimisation of the TP is crucial to its successful implementation.

Association mapping and validation of QTLs for flour yield in the soft winter wheat variety Kitahonami

The winter wheat variety Kitahonami shows a superior flour yield in comparison to other Japanese soft wheat varieties. To map the quantitative trait loci (QTL) associated with this trait, association mapping was performed using a panel of lines from Kitahonami's pedigree, along with leading Japanese varieties and advanced breeding lines. Using a mixed linear model corrected for kernel types and familial relatedness, 62 marker-trait associations for flour yield were identified and classified into 21 QTLs. In eighteen of these, Kitahonami alleles showed positive effects. Pedigree analysis demonstrated that a continuous pyramiding of QTLs had occurred throughout the breeding history of Kitahonami. Linkage analyses using three sets of doubled haploid populations from crosses in which Kitahonami was used as a parent were performed, leading to the validation of five of the eight QTLs tested. Among these, QTLs on chromosomes 3B and 7A showed highly significant and consistent effects across the three populations. This study shows that pedigree-based association mapping using breeding materials can be a useful method for QTL identification at the early stages of breeding programs.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Genome-wide association mapping for seedling and adult plant resistance to stripe rust in synthetic hexaploid wheat

Use of genetic diversity from related wild and domesticated species has made a significant contribution to improving wheat productivity. Synthetic hexaploid wheats (SHWs) exhibit natural genetic variation for resistance and/or tolerance to biotic and abiotic stresses. Stripe rust caused by (Puccinia striiformis f. sp. tritici; Pst), is an important disease of wheat worldwide. To characterise loci conferring resistance to stripe rust in SHWs, we conducted a genome-wide association study (GWAS) with a panel of 181 SHWs using the wheat 9 K SNP iSelect array. The SHWs were evaluated for their response to the prevailing races of Pst at the seedling and adult plant stages, the latter in replicated field trials at two sites in Ethiopia in 2011. About 28% of the SHWs exhibited immunity at the seedling stage while 56% and 83% were resistant to Pst at the adult plant stage at Meraro and Arsi Robe, respectively. A total of 27 SNPs in nine genomic regions (1 BS, 2 AS, 2 BL, 3 BL, 3 DL, 5A, 5 BL, 6DS and 7A) were linked with resistance to Pst at the seedling stage, while 38 SNPs on 18 genomic regions were associated with resistance at the adult plant stage. Six genomic regions were commonly detected at both locations using a mixed linear model corrected for population structure, kinship relatedness and adjusted for false discovery rate (FDR). The loci on chromosome regions 1 AS, 3 DL, 6 DS and 7 AL appeared to be novel QTL; our results confirm that resynthesized wheat involving its progenitor species is a rich source of new stripe (yellow) rust resistance that may be useful in choosing SHWs and incorporating diverse yellow rust (YR) resistance loci into locally adapted wheat cultivars.

A chromosome-based draft sequence of the hexaploid bread wheat (Triticum aestivum) genome

An ordered draft sequence of the 17-gigabase hexaploid bread wheat (Triticum aestivum) genome has been produced by sequencing isolated chromosome arms. We have annotated 124,201 gene loci distributed nearly evenly across the homeologous chromosomes and subgenomes. Comparative gene analysis of wheat subgenomes and extant diploid and tetraploid wheat relatives showed that high sequence similarity and structural conservation are retained, with limited gene loss, after polyploidization. However, across the genomes there was evidence of dynamic gene gain, loss, and duplication since the divergence of the wheat lineages. A high degree of transcriptional autonomy and no global dominance was found for the subgenomes. These insights into the genome biology of a polyploid crop provide a springboard for faster gene isolation, rapid genetic marker development, and precise breeding to meet the needs of increasing food demand worldwide.

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.

Transcriptome and complexity-reduced, DNA-based identification of intraspecies single-nucleotide polymorphisms in the polyploid Gossypium hirsutum L

Varietal single nucleotide polymorphisms (SNPs) are the differences within one of the two subgenomes between different tetraploid cotton varieties and have not been practically used in cotton genetics and breeding because they are difficult to identify due to low genetic diversity and very high sequence identity between homeologous genes in cotton. We have used transcriptome and restriction site-associated DNA sequencing to identify varietal SNPs among 18 G. hirsutum varieties based on the rationale that varietal SNPs can be more confidently called when flanked by subgenome-specific SNPs. Using transcriptome data, we successfully identified 37,413 varietal SNPs and, of these, 22,121 did not have an additional varietal SNP within their 20-bp flanking regions so can be used in most SNP genotyping assays. From restriction site-associated DNA sequencing data, we identified an additional 3090 varietal SNPs between two of the varieties. Of the 1583 successful SNP assays achieved using different genotyping platforms, 1363 were verified. Many of the SNPs behaved as dominant markers because of coamplification from homeologous loci, but the number of SNPs acting as codominant markers increased when one or more subgenome-specific SNP(s) were incorporated in their assay primers, giving them greater utility for breeding applications. A G. hirsutum genetic map with 1244 SNP markers was constructed covering 5557.42 centiMorgan and used to map qualitative and quantitative traits. This collection of G. hirsutum varietal SNPs complements existing intra-specific SNPs and provides the cotton community with a valuable marker resource applicable to genetic analyses and breeding programs.

Identification and independent validation of a stable yield and thousand grain weight QTL on chromosome 6A of hexaploid wheat (Triticum aestivum L.)

BACKGROUND

Grain yield in wheat is a polygenic trait that is influenced by environmental and genetic interactions at all stages of the plant's growth. Yield is usually broken down into three components; number of spikes per area, grain number per spike, and grain weight (TGW). In polyploid wheat, studies have identified quantitative trait loci (QTL) which affect TGW, yet few have been validated and fine-mapped using independent germplasm, thereby having limited impact in breeding.

RESULTS

In this study we identified a major QTL for TGW, yield and green canopy duration on wheat chromosome 6A of the Spark x Rialto population, across 12 North European environments. Using independent germplasm in the form of BC2 and BC4 near isogenic lines (NILs), we validated the three QTL effects across environments. In four of the five experiments the Rialto 6A introgression gave significant improvements in yield (5.5%) and TGW (5.1%), with morphometric measurements showing that the increased grain weight was a result of wider grains. The extended green canopy duration associated with the high yielding/TGW Rialto allele was comprised of two independent effects; earlier flowering and delayed final maturity, and was expressed stably across the five environments. The wheat homologue (TaGW2) of a rice gene associated with increased TGW and grain width was mapped within the QTL interval. However, no polymorphisms were identified in the coding sequence between the parents.

CONCLUSION

The discovery and validation through near-isogenic lines of robust QTL which affect yield, green canopy duration, thousand grain weight, and grain width on chromosome 6A of hexaploid wheat provide an important first step to advance our understanding of the genetic mechanisms regulating the complex processes governing grain size and yield in polyploid wheat.

Genetic and molecular bases of yield-associated traits: a translational biology approach between rice and wheat

Transferring the knowledge bases between related species may assist in enlarging the yield potential of crop plants. Being cereals, rice and wheat share a high level of gene conservation; however, they differ at metabolic levels as a part of the environmental adaptation resulting in different yield capacities. This review focuses on the current understanding of genetic and molecular regulation of yield-associated traits in both crop species, highlights the similarities and differences and presents the putative knowledge gaps. We focus on the traits associated with phenology, photosynthesis, and assimilate partitioning and lodging resistance; the most important drivers of yield potential. Currently, there are large knowledge gaps in the genetic and molecular control of such major biological processes that can be filled in a translational biology approach in transferring genomics and genetics informations between rice and wheat.

Biofortification of wheat grain with iron and zinc: integrating novel genomic resources and knowledge from model crops

Wheat, like many other staple cereals, contains low levels of the essential micronutrients iron and zinc. Up to two billion people worldwide suffer from iron and zinc deficiencies, particularly in regions with predominantly cereal-based diets. Although wheat flour is commonly fortified during processing, an attractive and more sustainable solution is biofortification, which requires developing new varieties of wheat with inherently higher iron and zinc content in their grains. Until now most studies aimed at increasing iron and zinc content in wheat grains have focused on discovering natural variation in progenitor or related species. However, recent developments in genomics and transformation have led to a step change in targeted research on wheat at a molecular level. We discuss promising approaches to improve iron and zinc content in wheat using knowledge gained in model grasses. We explore how the latest resources developed in wheat, including sequenced genomes and mutant populations, can be exploited for biofortification. We also highlight the key research and practical challenges that remain in improving iron and zinc content in wheat.

Genome-wide marker development for the wheat D genome based on single nucleotide polymorphisms identified from transcripts in the wild wheat progenitor Aegilops tauschii

13,347 high-confidence SNPs were discovered through transcriptome sequencing of Aegilops tauschii, which are useful for genomic analysis and molecular breeding of hexaploid wheat. In organisms with large and complex genomes, such as wheat, RNA-seq analysis is cost-effective for discovery of genome-wide single nucleotide polymorphisms (SNPs). In this study, deep sequencing of the spike transcriptome from two Aegilops tauschii accessions representing two major lineages led to the discovery of 13,347 high-confidence (HC) SNPs in 4,872 contigs. After removing redundant SNPs detected in the leaf transcriptome from the same accessions in an earlier study, 10,589 new SNPs were discovered. In total, 5,642 out of 5,808 contigs with HC SNPs were assigned to the Ae. tauschii draft genome sequence. On average, 732 HC polymorphic contigs were mapped in silico to each Ae. tauschii chromosome. Based on the polymorphic data, we developed markers to target the short arm of chromosome 2D and validated the polymorphisms using 20 Ae. tauschii accessions. Of the 29 polymorphic markers, 28 were successfully mapped to 2DS in the diploid F2 population of Ae. tauschii. Among ten hexaploid wheat lines, which included wheat synthetics and common wheat cultivars, 25 of the 43 markers were polymorphic. In the hexaploid F2 population between a common wheat cultivar and a synthetic wheat line, 23 of the 25 polymorphic markers between the parents were available for genotyping of the F2 plants and 22 markers mapped to chromosome 2DS. These results indicate that molecular markers that developed from polymorphisms between two distinct lineages of Ae. tauschii might be useful for analysis not only of the diploid, but also of the hexaploid wheat genome.

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.

Strategy for exploiting exotic germplasm using genetic, morphological, and environmental diversity: the Aegilops tauschii Coss. example

Hexaploid bread wheat evolved from a rare hybridisation, which resulted in a loss of genetic diversity in the wheat D-genome with respect to the ancestral donor, Aegilops tauschii. Novel genetic variation can be introduced into modern wheat by recreating the above hybridisation; however, the information associated with the Ae. tauschii accessions in germplasm collections is limited, making rational selection of accessions into a re-synthesis programme difficult. We describe methodologies to identify novel diversity from Ae. tauschii accessions that combines Bayesian analysis of genotypic data, sub-species diversity and geographic information that summarises variation in climate and habitat at the collection point for each accession. Comparisons were made between diversity discovered amongst a panel of Ae. tauschii accessions, bread wheat varieties and lines from the CIMMYT synthetic hexaploid wheat programme. The selection of Ae. tauschii accessions based on differing approaches had significant effect on diversity within each set. Our results suggest that a strategy that combines several criteria will be most effective in maximising the sampled variation across multiple parameters. The analysis of multiple layers of variation in ex situ Ae. tauschii collections allows for an informed and rational approach to the inclusion of wild relatives into crop breeding programmes.