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

Deciphering spike architecture formation towards yield improvement in wheat

Wheat is the most widely grown crop globally, providing 20% of the daily consumed calories and protein content around the world. With the growing global population and frequent occurrence of extreme weather caused by climate change, ensuring adequate wheat production is essential for food security. The architecture of the inflorescence plays a crucial role in determining the grain number and size, which is a key trait for improving yield. Recent advances in wheat genomics and gene cloning techniques have improved our understanding of wheat spike development and its applications in breeding practices. Here, we summarize the genetic regulation network governing wheat spike formation, the strategies used for identifying and studying the key factors affecting spike architecture, and the progress made in breeding applications. Additionally, we highlight future directions that will aid in the regulatory mechanistic study of wheat spike determination and targeted breeding for grain yield improvement.

Precisely mapping a major QTL for grain weight on chromosome 5B of the founder parent Chuanmai42 in the wheat-growing region of southwestern China

KEY MESSAGE

QTgw.saas-5B was validated as a major thousand-grain weight-related QTL in a founder parent used for wheat breeding and then precisely mapped to a 0.6 cM interval. Increasing the thousand-grain weight (TGW) is considered to be one of the most important ways to improve yield, which is a core objective among wheat breeders. Chuanmai42, which is a wheat cultivar with high TGW and a high and stable yield, is a parent of more than 30 new varieties grown in southwestern China. In this study, a Chuanmai42-derived recombinant inbred line (RIL) population was used to dissect the genetic basis of TGW. A major QTL (QTgw.saas-5B) mapped to the Xgwm213-Xgwm540 interval on chromosome 5B of Chuanmai42 explained up to 20% of the phenotypic variation. Using 71 recombinants with a recombination in the QTgw.saas-5B interval identified from a secondary RIL population comprising 1818 lines constructed by crossing the QTgw.saas-5B near-isogenic line with the recurrent parent Chuannong16, QTgw.saas-5B was delimited to a 0.6 cM interval, corresponding to a 21.83 Mb physical interval in the Chinese Spring genome. These findings provide the foundation for QTgw.saas-5B cloning and its use in molecular marker-assisted breeding.

Metabolome dynamics during wheat domestication

One of the most important crops worldwide is wheat. Wheat domestication took place about 10,000 years ago. Not only that its wild progenitors have been discovered and phenotypically characterized, but their genomes were also sequenced and compared to modern wheat. While comparative genomics is essential to track genes that contribute to improvement in crop yield, comparative analyses of functional biological end-products, such as metabolites, are still lacking. With the advent of rigorous mass-spectrometry technologies, it is now possible to address that problem on a big-data scale. In attempt to reveal classes of metabolites, which are associated with wheat domestication, we analyzed the metabolomes of wheat kernel samples from various wheat lines. These wheat lines represented subspecies of tetraploid wheat along primary and secondary domestications, including wild emmer, domesticated emmer, landraces durum, and modern durum. We detected that the groups of plant metabolites such as plant-defense metabolites, antioxidants and plant hormones underwent significant changes during wheat domestication. Our data suggest that these metabolites may have contributed to the improvement in the agricultural fitness of wheat. Closer evaluation of specific metabolic pathways may result in the future in genetically-engineered high-yield crops.

Genotyping-by-Sequencing Based Molecular Genetic Diversity of Pakistani Bread Wheat (Triticum aestivum L.) Accessions

Spring wheat (Triticum aestivum L.) is one of the most imperative staple food crops, with an annual production of 765 million tons globally to feed ∼40% world population. Genetic diversity in available germplasm is crucial for sustainable wheat improvement to ensure global food security. A diversity panel of 184 Pakistani wheat accessions was genotyped using 123,596 high-quality single nucleotide polymorphism (SNP) markers generated by genotyping-by-sequencing with 42% of the SNPs mapped on B, 36% on A, and 22% on D sub-genomes of wheat. Chromosome 2B contains the most SNPs (9,126), whereas 4D has the least (2,660) markers. The mean polymorphic information content, genetic diversity, and major allele frequency of the population were 0.157, 0.1844, and 0.87, respectively. Analysis of molecular variance revealed a higher genetic diversity (80%) within the sub-population than among the sub-populations (20%). The genome-wide linkage disequilibrium was 0.34 Mbp for the whole wheat genome. Among the three subgenomes, A has the highest LD decay value (0.29 Mbp), followed by B (0.2 Mbp) and D (0.07 Mbp) genomes, respectively. The results of population structure, principal coordinate analysis, phylogenetic tree, and kinship analysis also divided the whole population into three clusters comprising 31, 33, and 120 accessions in group 1, group 2, and group 3, respectively. All groups were dominated by the local wheat accessions. Estimation of genetic diversity will be a baseline for the selection of breeding parents for mutations and the genome-wide association and marker-assisted selection studies.

Identification and Validation of Quantitative Trait Loci for Wheat Dwarf Virus Resistance in Wheat (Triticum spp.)

Wheat dwarf virus (WDV) is transmitted by the leafhopper Psammotettix alienus. As a major pathogen in wheat and other cereals, WDV causes high yield losses in many European countries. Due to climate change, insect-transmitted viruses will become more important and the restrictions in the use of insecticides efficient against P. alienus renders growing of WDV resistant/tolerant varieties the only effective strategy to control WDV. So far, there is little information about the possible sources of resistance and no known information about the genome regions responsible for the resistance. In a screening for WDV resistance using artificial inoculation in gauze houses, a panel of 500 wheat accessions including cultivars, gene bank accessions, and wild relatives of wheat was phenotyped for virus titer, infection rate, as well as plant height and yield parameters relative to healthy controls of the same genotype. Additionally, 85 T. aestivum-Ae. tauschii intogression lines were tested for WDV resistance in the greenhouse. A subset of 250 hexaploid wheat accessions was genotyped with the 15k iSelect SNP Chip. By genome-wide association study (GWAS), the quantitative trait loci (QTL) for partial WDV resistance were identified. Within these studies, one cultivar was identified showing an average infection rate of only 5.7%. By analyzing single seed descent (SSD) and doubled haploid (DH) populations comprising 153 and 314 individuals for WDV resistance and by genotyping these with the 25k iSelect SNP Chip, QTL for yield per plant, thousand-grain weight, and relative virus titer were validated on chromosomes 1B, 2B, 3B, 4B, 4A, 5A, 6A, and 7A. These results will be the basis for marker-assisted selection for WDV resistance to replacing the laborious, time-consuming, and technically challenging phenotyping with WDV bearing leafhoppers.

Chromosome-specific KASP markers for detecting Amblyopyrum muticum segments in wheat introgression lines

Many wild-relative species are being used in prebreeding programs to increase the genetic diversity of wheat (Triticum aestivum L.). Genotyping tools such as single nucleotide polymorphism (SNP)-based arrays and molecular markers have been widely used to characterize wheat-wild relative introgression lines. However, due to the polyploid nature of the recipient wheat genome, it is difficult to develop SNP-based Kompetitive allele-specific polymerase chain reaction (KASP) markers that are codominant to track the introgressions from the wild species. Previous attempts to develop KASP markers have involved both exome- and polymerase chain reaction (PCR)-amplicon-based sequencing of the wild species. But chromosome-specific KASP assays have been hindered by homoeologous SNPs within the wheat genome. This study involved whole genome sequencing of the diploid wheat wild relative Amblyopyrum muticum (Boiss.) Eig and development of a de novo SNP discovery pipeline that generated ∼38,000 SNPs in unique wheat genome sequences. New assays were designed to increase the density of Am. muticum polymorphic KASP markers. With a goal of one marker per 60 Mbp, 335 new KASP assays were validated as diagnostic for Am. muticum in a wheat background. Together with assays validated in previous studies, 498 well distributed chromosome-specific markers were used to recharacterize previously genotyped wheat-Am. muticum doubled haploid (DH) introgression lines. The chromosome-specific nature of the KASP markers allowed clarification of which wheat chromosomes were involved with recombination events or substituted with Am. muticum chromosomes and the higher density of markers allowed detection of new small introgressions in these DH lines.

CerealsDB: A Whistle-Stop Tour of an Open Access SNP Resource

The CerealsDB website, created by members of the Functional Genomics Group at the University of Bristol, provides access to a database containing SNP and genotyping data for hexaploid wheat and, to a lesser extent, its progenitors and several of its relatives. The site is principally aimed at plant breeders and research scientists who wish to obtain information regarding SNP markers; for example, obtain primers used for their identification or the sequences upon which they are based. The database underpinning the website contains circa one million putative varietal SNPs of which several hundreds of thousands have been experimentally validated on a range of common genotyping platforms. For each SNP marker, the site also hosts the allelic scores for thousands of elite wheat varieties, landrace cultivars, and wheat relatives. Tools are available to help negotiate and visualize the datasets. The website has been designed to be simple and straightforward to use and is completely open access.

Wheat in vivo RNA structure landscape reveals a prevalent role of RNA structure in modulating translational subgenome expression asymmetry

BACKGROUND

Polyploidy, especially allopolyploidy, which entails merging divergent genomes via hybridization and whole-genome duplication (WGD), is a major route to speciation in plants. The duplication among the parental genomes (subgenomes) often leads to one subgenome becoming dominant over the other(s), resulting in subgenome asymmetry in gene content and expression. Polyploid wheats are allopolyploids with most genes present in two (tetraploid) or three (hexaploid) functional copies, which commonly show subgenome expression asymmetry. It is unknown whether a similar subgenome asymmetry exists during translation. We aim to address this key biological question and explore the major contributing factors to subgenome translation asymmetry.

RESULTS

Here, we obtain the first tetraploid wheat translatome and reveal that subgenome expression asymmetry exists at the translational level. We further perform in vivo RNA structure profiling to obtain the wheat RNA structure landscape and find that mRNA structure has a strong impact on translation, independent of GC content. We discover a previously uncharacterized contribution of RNA structure in subgenome translation asymmetry. We identify 3564 single-nucleotide variations (SNVs) across the transcriptomes between the two tetraploid wheat subgenomes, which induce large RNA structure disparities. These SNVs are highly conserved within durum wheat cultivars but are divergent in both domesticated and wild emmer wheat.

CONCLUSIONS

We successfully determine both the translatome and in vivo RNA structurome in tetraploid wheat. We reveal that RNA structure serves as an important modulator of translational subgenome expression asymmetry in polyploids. Our work provides a new perspective for molecular breeding of major polyploid crops.

Quantitative trait loci for yellow rust resistance in spring wheat doubled haploid populations developed from the German Federal ex situ genebank genetic resources

Novel resistance sources to the pathogen Puccinia striiformis f. sp. tritici, which causes yellow rust (stripe rust), a widespread devastating foliar disease in wheat (Triticum aestivum L.), are in demand. Here, we tested two doubled haploid (DH) spring wheat populations derived from the genetic resources for resistance to yellow rust in field trials in Germany and Egypt. Additionally, we performed tests for all-stage resistance (seedling resistance). We performed linkage mapping based on 15k Infinium SNP chip genotyping data that resulted in 3,567 and 3,457 polymorphic markers for DH Population 1 (103 genotypes) and DH Population 2 (148 genotypes), respectively. In DH Population 1, we identified a major and consistent quantitative trait locus (QTL) on chromosome 1B that explained up to 28 and 39% of the phenotypic variation in the field and seedling tests, respectively. The favorable allele was contributed by the line 'TRI-5645', a landrace from Iran, and is most probably the yellow rust resistance (Yr) gene Yr10. In DH Population 2, the favorable allele of a major QTL on chromosome 6B was contributed by the line 'TRI-5310', representing the variety 'Eureke' from France. This QTL was mainly effective in the German environments and explained up to 36% of the phenotypic variation. In Egypt, however, only a moderate resistance QTL was identified in the field tests and no resistance QTL was observed in the seedling tests. Our results demonstrate the usefulness of genetic resources to identify novel sources of resistance to yellow rust, including the "Warrior" race PstS10.

Identification of a major and stable QTL on chromosome 5A confers spike length in wheat (Triticum aestivum L.)

Spike length (SL) is the key determinant of plant architecture and yield potential. In this study, 193 recombinant inbred lines (RILs) derived from a cross between 13F10 and Chuanmai 42 (CM42) were evaluated for spike length in six environments. Sixty RILs consisting of 30 high and 30 low SLs were genotyped using the bulked segregant analysis exome sequencing (BSE-Seq) analysis for preliminary quantitative trait locus (QTL) mapping. A 6.69 Mb (518.43-525.12 Mb) region on chromosome 5AL was found to have a significant effect on the SL trait. Fifteen competitive allele-specific PCR (KASP) markers were successfully converted from the single nucleotide polymorphisms (SNPs) in the SL target region. Combined with four novel simple sequence repeat (SSR) markers, a genetic linkage map spanning 21.159 cM was constructed. The mapping result confirmed the identity of a major and stable QTL named QSl.cib-5A in the targeted region that explained 7.88-26.60% of the phenotypic variation in SL. QSl.cib-5A was narrowed to a region of 4.84 cM interval corresponding to a 4.67 Mb (516.60-521.27 Mb) physical region in the Chinese Spring RefSeq v2.0 containing 17 high-confidence genes with 25 transcripts. In addition, this QTL exhibited pleiotropic effects on spikelet density (SD), with the phenotypic variances proportion ranging from 11.34 to 19.92%. This study provides a foundational step for cloning the QSl.cib-5A, which is involved in the regulation of spike morphology in common wheat.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01249-6.

Genetic Diversity and Population Structure Analysis of Triticum aestivum L. Landrace Panel from Afghanistan

Landraces are a potential source of genetic diversity and provide useful genetic resources to cope with the current and future challenges in crop breeding. Afghanistan is located close to the centre of origin of hexaploid wheat. Therefore, understanding the population structure and genetic diversity of Afghan wheat landraces is of enormous importance in breeding programmes for the development of high-yielding cultivars as well as broadening the genetic base of bread wheat. Here, a panel of 363 bread wheat landraces collected from seven north and north-eastern provinces of Afghanistan were evaluated for population structure and genetic diversity using single nucleotide polymorphic markers (SNPs). The genotyping-by-sequencing of studied landraces after quality control provided 4897 high-quality SNPs distributed across the genomes A (33.75%), B (38.73%), and D (27.50%). The population structure analysis was carried out by two methods using model-based STRUCTURE analysis and cluster-based discriminant analysis of principal components (DAPC). The analysis of molecular variance showed a higher proportion of variation within the sub-populations compared with the variation observed as a whole between sub-populations. STRUCTURE and DAPC analysis grouped the majority of the landraces from Badakhshan and Takhar together in one cluster and the landraces from Baghlan and Kunduz in a second cluster, which is in accordance with the micro-climatic conditions prevalent within the north-eastern agro-ecological zone. Genetic distance analysis was also studied to identify differences among the Afghan regions; the strongest correlation was observed for the Badakhshan and Takhar (0.003), whereas Samangan and Konarha (0.399) showed the highest genetic distance. The population structure and genetic diversity analysis highlighted the complex genetic variation present in the landraces which were highly correlated to the geographic origin and micro-climatic conditions within the agro-climatic zones of the landraces. The higher proportions of admixture could be attributed to historical unsupervised exchanges of seeds between the farmers of the central and north-eastern provinces of Afghanistan. The results of this study will provide useful information for genetic improvement in wheat and is essential for association mapping and genomic prediction studies to identify novel sources for resistance to abiotic and biotic stresses.

QTL Mapping and Validation for Kernel Area and Circumference in Common Wheat via High-Density SNP-Based Genotyping

As an important component, 1,000 kernel weight (TKW) plays a significant role in the formation of yield traits of wheat. Kernel size is significantly positively correlated to TKW. Although numerous loci for kernel size in wheat have been reported, our knowledge on loci for kernel area (KA) and kernel circumference (KC) remains limited. In the present study, a recombinant inbred lines (RIL) population containing 371 lines genotyped using the Wheat55K SNP array was used to map quantitative trait loci (QTLs) controlling the KA and KC in multiple environments. A total of 54 and 44 QTLs were mapped by using the biparental population or multienvironment trial module of the inclusive composite interval mapping method, respectively. Twenty-two QTLs were considered major QTLs. BLAST analysis showed that major and stable QTLs QKc.sau-6A.1 (23.12-31.64 cM on 6A) for KC and QKa.sau-6A.2 (66.00-66.57 cM on 6A) for KA were likely novel QTLs, which explained 22.25 and 20.34% of the phenotypic variation on average in the 3 year experiments, respectively. Two Kompetitive allele-specific PCR (KASP) markers, KASP-AX-109894590 and KASP-AX-109380327, were developed and tightly linked to QKc.sau-6A.1 and QKa.sau-6A.2, respectively, and the genetic effects of the different genotypes in the RIL population were successfully confirmed. Furthermore, in the interval where QKa.sau-6A.2 was located on Chinese Spring and T. Turgidum ssp. dicoccoides reference genomes, only 11 genes were found. In addition, digenic epistatic QTLs also showed a significant influence on KC and KA. Altogether, the results revealed the genetic basis of KA and KC and will be useful for the marker-assisted selection of lines with different kernel sizes, laying the foundation for the fine mapping and cloning of the gene(s) underlying the stable QTLs detected in this study.

Genome-Wide Association Mapping of Seedling Drought Tolerance in Winter Wheat

In the southern Great Plains of the United States, winter wheat grown for dual-purpose is often planted early, which puts it at risk for drought stress at the seedling stage in the autumn. To map quantitative trait loci (QTL) associated with seedling drought tolerance, a genome-wide association study (GWAS) was performed on a hard winter wheat association mapping panel. Two sets of plants were planted in the greenhouse initially under well-watered conditions. At the five-leaf stage, one set continued to receive the optimum amount of water, whereas watering was withdrawn from the other set (drought stress treatment) for 14 days to mimic drought stress. Large phenotypic variation was observed in leaf chlorophyll content, leaf chlorophyll fluorescence, shoot length, number of leaves per seedling, and seedling recovery. A mixed linear model analysis detected multiple significant QTL associated with seedling drought tolerance-related traits on chromosomes 1B, 2A, 2B, 2D, 3A, 3B, 3D, 4B, 5A, 5B, 6B, and 7B. Among those, 12 stable QTL responding to drought stress for various traits were identified. Shoot length and leaf chlorophyll fluorescence were good indicators in responding to drought stress because most of the drought responding QTL detected using means of these two traits were also detected in at least two experimental repeats. These stable QTL are more valuable for use in marker-assisted selection during wheat breeding. Moreover, different traits were mapped on several common chromosomes, such as 1B, 2B, 3B, and 6B, and two QTL clusters associated with three or more traits were located at 107-130 and 80-83 cM on chromosomes 2B and 6B, respectively. Furthermore, some QTL detected in this study co-localized with previously reported QTL for root and shoot traits at the seedling stage and canopy temperature at the grain-filling stage of wheat. In addition, several of the mapped chromosomes were also associated with drought tolerance during the flowering or grain-filling stage in wheat. Some significant single-nucleotide polymorphisms (SNPs) were aligned to candidate genes playing roles in plant abiotic stress responses. The SNP markers identified in this study will be further validated and used for marker-assisted breeding of seedling drought tolerance during dual-purpose wheat breeding.

An extended KASP-SNP resource for molecular breeding in Chinese cabbage(Brassica rapa L. ssp. pekinensis)

Kompetitive allele-specific PCR (KASP) is a cost-effective single-step SNP genotyping technology, With an objective to enhance the marker repertoire and develop high efficient KASP-SNP markers in Chinese cabbage, we re-sequenced four Chinese cabbage doubled haploid (DH) lines, Y177-47, Y635-10, Y510-1 and Y510-9, and generated a total of more than 38.5 billion clean base pairs. A total of 827,720 SNP loci were identified with an estimated density of 3,217 SNPs/Mb. Further, a total of 387,354 SNPs with at least 30 bp to the next most adjacent SNPs on either side were selected as resource for KASP markers. From this resource, 258 (96.27%) of 268 SNP loci were successfully transformed into KASP-SNP markers using a Roche LightCycler 480-II instrument. Among these markers, 221 (85.66%) were co-dominant markers, 220 (85.27%) were non-synonymous SNPs, and 257 (99.6%) were newly developed markers. In addition, 53 markers were applied for genotyping of 34 Brassica rapa accessions. Cluster analysis separated these 34 accessions into three clusters based on heading types. The millions of SNP loci, a large set of resource for KASP markers, as well as the newly developed KASP markers in this study may facilitate further genetic and molecular breeding studies in Brassica rapa.

Detection of genomic regions associated with tiller number in Iranian bread wheat under different water regimes using genome-wide association study

Two of the important traits for wheat yield are tiller and fertile tiller number, both of which have been thought to increase cereal yield in favorable and unfavorable environments. A total of 6,349 single nucleotide polymorphism (SNP) markers from the 15 K wheat Infinium array were employed for genome-wide association study (GWAS) of tillering number traits, generating a physical distance of 14,041.6 Mb based on the IWGSC wheat genome sequence. GWAS analysis using Fixed and random model Circulating Probability Unification (FarmCPU) identified a total of 47 significant marker-trait associations (MTAs) for total tiller number (TTN) and fertile tiller number (FTN) in Iranian bread wheat under different water regimes. After applying a 5% false discovery rate (FDR) threshold, a total of 13 and 11 MTAs distributed on 10 chromosomes were found to be significantly associated with TTN and FTN, respectively. Linked single nucleotide polymorphisms for IWB39005 (2A) and IWB44377 (7A) were highly significantly associated (FDR < 0.01) with TTN and FTN traits. Moreover, to validate GWAS results, meta-analysis was performed and 30 meta-QTL regions were identified on 11 chromosomes. The integration of GWAS and meta-QTLs revealed that tillering trait in wheat is a complex trait which is conditioned by the combined effects of minor changes in multiple genes. The information provided by this study can enrich the currently available candidate genes and genetic resources pools, offering evidence for subsequent analysis of genetic adaptation of wheat to different climatic conditions of Iran and other countries.

Overcoming polyploidy pitfalls: a user guide for effective SNP conversion into KASP markers in wheat

Conversion of SNP chip assays into locus-specific KASP markers requires adapted strategies in polyploid species with high genome homeology. Procedures are exemplified by QTL-associated SNPs in hexaploid wheat. Kompetitive allele-specific PCR (KASP) markers are commonly used in marker-assisted commercial plant breeding due to their cost-effectiveness and throughput for high sample volumes. However, conversion of trait-linked SNP markers from array-based SNP detection technologies into KASP markers is particularly challenging in polyploid crop species, due to the presence of highly similar homeologous and paralogous genome sequences. We evaluated strategies and identified key requirements for successful conversion of Illumina Infinium assays from the wheat 90 K SNP array into robust locus-specific KASP markers. Numerous examples showed that commonly used software for semiautomated KASP primer design frequently fails to achieve locus-specificity of KASP assays in wheat. Instead, alignment of SNP probes with multiple reference genomes and Sanger sequencing of relevant genotypes, followed by visual KASP primer placement, was critical for locus-specificity. To identify KASP assays resulting in false calling of heterozygous individuals, validation of KASP assays using extended reference genotype sets including heterozygous genotypes is strongly advised for polyploid crop species. Applying this strategy, we developed highly reproducible, stable KASP assays that are predictive for root biomass QTL haplotypes from highly homoeologous wheat chromosome regions. Due to their locus-specificity, these assays predicted root biomass considerably better than the original trait-associated markers from the Illumina array.

Genome-Wide Association Studies in Diverse Spring Wheat Panel for Stripe, Stem, and Leaf Rust Resistance

Among several important wheat foliar diseases, Stripe rust (YR), Leaf rust (LR), and Stem rust (SR) have always been an issue of concern to the farmers and wheat breeders. Evolution of virulent pathotypes of these rusts has posed frequent threats to an epidemic. Pyramiding rust-resistant genes are the most economical and environment-friendly approach in postponing this inevitable threat. To achieve durable long term resistance against the three rusts, an attempt in this study was made searching for novel sources of resistant alleles in a panel of 483 spring wheat genotypes. This is a unique and comprehensive study where evaluation of a diverse panel comprising wheat germplasm from various categories and adapted to different wheat agro-climatic zones was challenged with 18 pathotypes of the three rusts with simultaneous screening in field conditions. The panel was genotyped using 35K SNP array and evaluated for each rust at two locations for two consecutive crop seasons. High heritability estimates of disease response were observed between environments for each rust type. A significant effect of population structure in the panel was visible in the disease response. Using a compressed mixed linear model approach, 25 genomic regions were found associated with resistance for at least two rusts. Out of these, seven were associated with all the three rusts on chromosome groups 1 and 6 along with 2B. For resistance against YR, LR, and SR, there were 16, 18, and 27 QTL (quantitative trait loci) identified respectively, associated at least in two out of four environments. Several of these regions got annotated with resistance associated genes viz. NB-LRR, E3-ubiquitin protein ligase, ABC transporter protein, etc. Alien introgressed (on 1B and 3D) and pleiotropic (on 7D) resistance genes were captured in seedling and adult plant disease responses, respectively. The present study demonstrates the use of genome-wide association for identification of a large number of favorable alleles for leaf, stripe, and stem rust resistance for broadening the genetic base. Quick conversion of these QTL into user-friendly markers will accelerate the deployment of these resistance loci in wheat breeding programs.

A framework for gene mapping in wheat demonstrated using the Yr7 yellow rust resistance gene

We used three approaches to map the yellow rust resistance gene Yr7 and identify associated SNPs in wheat. First, we used a traditional QTL mapping approach using a double haploid (DH) population and mapped Yr7 to a low-recombination region of chromosome 2B. To fine map the QTL, we then used an association mapping panel. Both populations were SNP array genotyped allowing alignment of QTL and genome-wide association scans based on common segregating SNPs. Analysis of the association panel spanning the QTL interval, narrowed the interval down to a single haplotype block. Finally, we used mapping-by-sequencing of resistant and susceptible DH bulks to identify a candidate gene in the interval showing high homology to a previously suggested Yr7 candidate and to populate the Yr7 interval with a higher density of polymorphisms. We highlight the power of combining mapping-by-sequencing, delivering a complete list of gene-based segregating polymorphisms in the interval with the high recombination, low LD precision of the association mapping panel. Our mapping-by-sequencing methodology is applicable to any trait and our results validate the approach in wheat, where with a near complete reference genome sequence, we are able to define a small interval containing the causative gene.

Exome sequencing of bulked segregants identified a novel TaMKK3-A allele linked to the wheat ERA8 ABA-hypersensitive germination phenotype

Using bulked segregant analysis of exome sequence, we fine-mapped the ABA-hypersensitive mutant ERA8 in a wheat backcross population to the TaMKK3-A locus of chromosome 4A. Preharvest sprouting (PHS) is the germination of mature grain on the mother plant when it rains before harvest. The ENHANCED RESPONSE TO ABA8 (ERA8) mutant increases seed dormancy and, consequently, PHS tolerance in soft white wheat 'Zak.' ERA8 was mapped to chromosome 4A in a Zak/'ZakERA8' backcross population using bulked segregant analysis of exome sequenced DNA (BSA-exome-seq). ERA8 was fine-mapped relative to mutagen-induced SNPs to a 4.6 Mb region containing 70 genes. In the backcross population, the ERA8 ABA-hypersensitive phenotype was strongly linked to a missense mutation in TaMKK3-A-G1093A (LOD 16.5), a gene associated with natural PHS tolerance in barley and wheat. The map position of ERA8 was confirmed in an 'Otis'/ZakERA8 but not in a 'Louise'/ZakERA8 mapping population. This is likely because Otis carries the same natural PHS susceptible MKK3-A-A660^S allele as Zak, whereas Louise carries the PHS-tolerant MKK3-A-C660^R allele. Thus, the variation for grain dormancy and PHS tolerance in the Louise/ZakERA8 population likely resulted from segregation of other loci rather than segregation for PHS tolerance at the MKK3 locus. This inadvertent complementation test suggests that the MKK3-A-G1093A mutation causes the ERA8 phenotype. Moreover, MKK3 was a known ABA signaling gene in the 70-gene 4.6 Mb ERA8 interval. None of these 70 genes showed the differential regulation in wild-type Zak versus ERA8 expected of a promoter mutation. Thus, the working model is that the ERA8 phenotype results from the MKK3-A-G1093A mutation.

Rapid identification of homozygosity and site of wild relative introgressions in wheat through chromosome-specific KASP genotyping assays

For future food security, it is important that wheat, one of the most widely consumed crops in the world, can survive the threat of abiotic and biotic stresses. New genetic variation is currently being introduced into wheat through introgressions from its wild relatives. For trait discovery, it is necessary that each introgression is homozygous and hence stable. Breeding programmes rely on efficient genotyping platforms for marker-assisted selection (MAS). Recently, single nucleotide polymorphism (SNP)-based markers have been made available on high-throughput Axiom^® SNP genotyping arrays. However, these arrays are inflexible in their design and sample numbers, making their use unsuitable for long-term MAS. SNPs can potentially be converted into Kompetitive allele-specific PCR (KASP™) assays that are comparatively cost-effective and efficient for low-density genotyping of introgression lines. However, due to the polyploid nature of wheat, KASP assays for homoeologous SNPs can have difficulty in distinguishing between heterozygous and homozygous hybrid lines in a backcross population. To identify co-dominant SNPs, that can differentiate between heterozygotes and homozygotes, we PCR-amplified and sequenced genomic DNA from potential single-copy regions of the wheat genome and compared them to orthologous copies from different wild relatives. A panel of 620 chromosome-specific KASP assays have been developed that allow rapid detection of wild relative segments and provide information on their homozygosity and site of introgression in the wheat genome. A set of 90 chromosome-nonspecific assays was also produced that can be used for genotyping introgression lines. These multipurpose KASP assays represent a powerful tool for wheat breeders worldwide.

QTL analysis and fine mapping of a QTL for yield-related traits in wheat grown in dry and hot environments

Genetic control of grain yield and phenology was examined in the Excalibur/Kukri doubled haploid mapping population grown in 32 field experiments across the climatic zones of southern Australia, India and north-western Mexico where the wheat crop experiences drought and heat stress. A total of 128 QTL were identified for four traits: grain yield, thousand grain weight (TGW), days to heading and grain filling duration. These QTL included 24 QTL for yield and 27 for TGW, showing significant interactions with the environment (Q * E). We also identified 14 QTL with a significant, small main effects on yield across environments. The study focussed on a region of chromosome 1B where two main effect QTL were found for yield and TGW without the confounding effect of phenology. Excalibur was the source of favourable alleles: QYld.aww-1B.2 with a peak at 149.5-150.1 cM and QTgw.aww-1B at 168.5-171.4 cM. We developed near isogenic lines (NIL) for the interval including QYld.aww-1B.2 and QTgw.aww-1B and evaluated them under semi-controlled conditions. Significant differences in four pairs of NIL were observed for grain yield but not for TGW, confirming a positive effect of the Excalibur allele for QYld.aww-1B.2. The interval containing QYld.aww-1B.2 was narrowed down to 2.9 cM which corresponded to a 2.2 Mbp genomic region on the chromosome 1B genomic reference sequence of cv. Chinese Spring and contained 39 predicted genes.

Genotyping-by-Sequencing Derived Single Nucleotide Polymorphisms Provide the First Well-Resolved Phylogeny for the Genus Triticum (Poaceae)

Wheat (Triticum spp.) has been an important staple food crop for mankind since the beginning of agriculture. The genus Triticum L. is composed of diploid, tetraploid, and hexaploid species, majority of which have not yet been discriminated clearly, and hence their phylogeny and classification remain unresolved. Genotyping-by-sequencing (GBS) is an easy and affordable method that allows us to generate genome-wide single nucleotide polymorphism (SNP) markers. In this study, we used GBS to obtain SNPs covering all seven chromosomes from 283 accessions of Triticum-related genera. After filtering low-quality and redundant SNPs based on haplotype information, the GBS assay provided 14,188 high-quality SNPs that were distributed across the A (71%), B (26%), and D (2.4%) genomes. Cluster analysis and discriminant analysis of principal components (DAPC) allowed us to distinguish six distinct groups that matched well with Triticum species complexity. We constructed a Bayesian phylogenetic tree using 14,188 SNPs, in which 17 Triticum species and subspecies were discriminated. Dendrogram analysis revealed that the polyploid wheat species could be divided into groups according to the presence of A, B, D, and G genomes with strong nodal support and provided new insight into the evolution of spelt wheat. A total of 2,692 species-specific SNPs were identified to discriminate the common (T. aestivum) and durum (T. turgidum) wheat cultivar and landraces. In principal component analysis grouping, the two wheat species formed individual clusters and the SNPs were able to distinguish up to nine groups of 10 subspecies. This study demonstrated that GBS-derived SNPs could be used efficiently in genebank management to classify Triticum species and subspecies that are very difficult to distinguish by their morphological characters.

Characterization of genetic diversity and population structure in wheat using array based SNP markers

Genetic diversity is crucial for successful adaptation and sustained improvement in crops. India is bestowed with diverse agro-climatic conditions which makes it rich in wheat germplasm adapted to various niches. Germplasm repository consists of local landraces, trait specific genetic stocks including introgressions from wild relatives, exotic collections, released varieties, and improved germplasm. Characterization of genetic diversity is done using morpho-physiological characters as well as by analyzing variations at DNA level. However, there are not many reports on array based high throughput SNP markers having characteristics of genome wide coverage employed in Indian spring wheat germplasm. Amongst wheat SNP arrays, 35K Axiom Wheat Breeder's Array has the highest SNP polymorphism efficiency suitable for genetic mapping and genetic diversity characterization. Therefore, genotyping was done using 35K in 483 wheat genotypes resulting in 14,650 quality filtered SNPs, that were distributed across the B (~ 50%), A (~ 39%), and D (~ 10%) genomes. The total genetic distance coverage was 4477.85 cM with 3.27 SNP/cM and 0.49 cM/SNP as average marker density and average inter-marker distance, respectively. The PIC ranged from 0.09 to 0.38 with an average of 0.29 across genomes. Population structure and Principal Coordinate Analysis resulted in two subpopulations (SP1 and SP2). The analysis of molecular variance revealed the genetic variation of 2% among and 98% within subpopulations indicating high gene flow between SP1 and SP2. The subpopulation SP2 showed high level of genetic diversity based on genetic diversity indices viz. Shannon's information index (I) = 0.648, expected heterozygosity (He) = 0.456 and unbiased expected heterozygosity (uHe) = 0.456. To the best of our knowledge, this study is the first to include the largest set of Indian wheat genotypes studied exclusively for genetic diversity. These findings may serve as a potential source for the identification of uncharacterized QTL/gene using genome wide association studies and marker assisted selection in wheat breeding programs.

Novel Informatic Tools to Support Functional Annotation of the Durum Wheat Genome

Seed mutagenesis is one strategy to create a population with thousands of useful mutations for the direct selection of desirable traits, to introduce diversity into varietal improvement programs, or to generate a mutant collection to support gene functional analysis. However, phenotyping such large collections, where each individual may carry many mutations, is a bottleneck for downstream analysis. Targeting Induced Local Lesions in Genomes (TILLinG), when coupled with next-generation sequencing allows high-throughput mutation discovery and selection by genotyping. We mutagenized an advanced durum breeding line, UAD0951096_F2:5 and performed short-read (2x125 bp) Illumina sequencing of the exome of 100 lines using an available exome capture platform. To improve variant calling, we generated a consolidated exome reference using the recently available genome sequences of the cultivars Svevo and Kronos to facilitate the alignment of reads from the UAD0951096_F2:5 derived mutants. The resulting exome reference was 484.4 Mbp. We also developed a user-friendly, searchable database and bioinformatic analysis pipeline that allowed us to predict zygosity of the mutations discovered and extracts flanking sequences for rapid marker development. Here, we present these tools with the aim of allowing researchers fast and accurate downstream selection of mutations discovered by TILLinG by sequencing to support functional annotation of the durum wheat genome.

Genetic Contribution of Synthetic Hexaploid Wheat to CIMMYT's Spring Bread Wheat Breeding Germplasm

Synthetic hexaploid (SH) wheat (AABBD'D') is developed by artificially generating a fertile hybrid between tetraploid durum wheat (Triticum turgidum, AABB) and diploid wild goat grass (Aegilops tauschii, D'D'). Over three decades, the International Maize and Wheat Improvement Center (CIMMYT) has developed and utilized SH wheat to bridge gene transfer from Ae. tauschii and durum wheat to hexaploid bread wheat. This is a unique example of success utilizing wild relatives in mainstream breeding at large scale worldwide. Our study aimed to determine the genetic contribution of SH wheat to CIMMYT's global spring bread wheat breeding program. We estimated the theoretical and empirical contribution of D' to synthetic derivative lines using the ancestral pedigree and marker information using over 1,600 advanced lines and their parents. The average marker-estimated D' contribution was 17.5% with difference in genome segments suggesting application of differential selection pressure. The pedigree-based contribution was correlated with marker-based estimates without providing chromosome segment specific variation. Results from international yield trials showed that 20% of the lines were synthetic derived with an average D' contribution of 15.6%. Our results underline the importance of SH wheat in maintaining and enhancing genetic diversity and genetic gain over years and is important for development of a more targeted introgression strategy. The study provides retrospective view into development and utilization of SH in the CIMMYT Global Wheat Program.

Genetic diversity and population structure analysis of synthetic and bread wheat accessions in Western Siberia

Recurrent selection and intercrossing between best of the best parents in each generation of breeding cycle resulted in a narrower genetic diversity in elite wheat (Triticum aestivum L.) germplasm. Therefore, we investigated diverse source of 143 synthetic and bread wheat accessions for identifying potentially rich genetic resources for improving the genetic diversity in wheat. This study identified 47,526 genotyping-by-sequencing-derived SNP markers that were nearly evenly distributed across three genomes of wheat. The population structure analysis identified three distinct clusters (Japan synthetics, CIMMYT synthetics, and bread wheat) of wheat genotypes on the basis of type and geographical origin of wheat accessions. Population differentiation using analysis of molecular variance indicated 21% of the total genetic variance among subgroups and the remainder within subgroups. This study also identified that the Japan synthetic group was the most divergent group compared with other subgroups. The genetic diversity comparisons between synthetic and bread wheat lines showed that the gene diversity of synthetic wheat was 33% higher than bread wheat accessions, indicating the potential use of these lines for broadening the genetic diversity of modern wheat cultivars. The results from this study will be helpful in further understanding genomic features of wheat and facilitate their use in wheat breeding programs.

Development of EST-based SSR and SNP markers in Gastrodia elata (herbal medicine) by sequencing, de novo assembly and annotation of the transcriptome

Tianma (Gastrodia elata Blume) has unique biological characteristics and high medicinal value. The wild resource of G. elata is being overutilized and should be conserved as it is already included in the list of endangered species in China. The population size of cultivated G. elata is small because of domestication bottleneck. Therefore, it is of utmost importance to evolve high-quality varieties and conserve wild resources of G. elata. In this study, we sequenced tuber transcriptomes of three major cultivated sub-species of Gastrodia elata, namely G. elata BI. f. elata, G. elata Bl. f. glauca S. Chow, and G. elata Bl. f. Viridis, and obtained about 7.8G clean data. The assembled high-quality reads of three sub-species were clustered into 56,884 unigenes. Of these, 31,224 (54.89%), 25,733 (45.24%), 22,629 (39.78%), and 11,856 (20.84%) unigenes were annotated by Nr, Swiss-Port, Eukaryotic Ortholog Groups (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, respectively. Here, a total of 3766 EST-SSRs and 128,921 SNPs were identified from the unigenes. The results not only offer huge number of genes that were responsible for the growth, development, and metabolism of bioactive components, but also a large number of molecular markers were detected for future studies on the conservation genetics and molecular breeding of G. elata.

Genome-wide Mapping for Stripe Rust Resistance Loci in Common Wheat Cultivar Qinnong 142

Stripe rust caused by Puccinia striiformis f. sp. tritici threatens worldwide wheat production. Growing resistant cultivars is the best way to control this disease. Chinese wheat cultivar Qinnong 142 (QN142) has a high level of adult-plant resistance to stripe rust. To identify quantitative trait loci (QTLs) related to stripe rust resistance, we developed a recombinant inbred line (RIL) population from a cross between QN142 and susceptible cultivar Avocet S. The parents and 165 F₆ RILs were evaluated in terms of their stripe rust infection type and disease severity in replicated field tests with six site-year environments. The parents and RILs were genotyped with single-nucleotide polymorphism (SNP) markers. Four stable QTLs were identified in QN142 and mapped to chromosome arms 1BL, 2AL, 2BL, and 6BS. The 1BL QTL was probably the known resistance gene Yr29, the 2BL QTL was in a resistance gene-rich region, and the 2AL and 6BS QTLs might be new. Kompetitive allele specific polymerase chain reaction markers developed from the SNP markers flanking these QTLs were highly polymorphic in a panel of 150 wheat cultivars and breeding lines. These markers could be used in marker-assisted selection for incorporating the stripe rust resistance QTL into new wheat cultivars.

Contribution of recent technological advances to future resistance breeding

The development of durable host resistance strategies to control crop diseases is a primary need for sustainable agricultural production in the future. This article highlights the potential of recent progress in the understanding of host resistance for future cereal breeding. Much of the novel work is based on advancements in large-scale sequencing and genomics, rapid gene isolation techniques and high-throughput molecular marker technologies. Moreover, emerging applications on the pathogen side like effector identification or field pathogenomics are discussed. The combination of knowledge from both sides of cereal pathosystems will result in new approaches for resistance breeding. We describe future applications and innovative strategies to implement effective and durable strategies to combat diseases of major cereal crops while reducing pesticide dependency.

Progress in plant genome sequencing: research directions

Since the first plant genome of Arabidopsis thaliana has been sequenced and published, genome sequencing technologies have undergone significant changes. New algorithms, sequencing technologies and bioinformatic approaches were adopted to obtain genome, transcriptome and exome sequences for model and crop species, which have permitted deep inferences into plant biology. As a result of an improved genome assembly and analysis methods, genome sequencing costs plummeted and the number of high-quality plant genome sequences is constantly growing. Consequently, more than 300 plant genome sequences have been published over the past twenty years. Although many of the published genomes are considered incomplete, they proved to be a valuable tool for identifying genes involved in the formation of economically valuable plant traits, for marker-assisted and genomic selection and for comparative analysis of plant genomes in order to determine the basic patterns of origin of various plant species. Since a high coverage and resolution of a genome sequence is not enough to detect all changes in complex samples, targeted sequencing, which consists in the isolation and sequencing of a specific region of the genome, has begun to develop. Targeted sequencing has a higher detection power (the ability to identify new differences/variants) and resolution (up to one basis). In addition, exome sequencing (the method of sequencing only protein-coding genes regions) is actively developed, which allows for the sequencing of non-expressed alleles and genes that cannot be found with RNA-seq. In this review, an analysis of sequencing technologies development and the construction of “reference” genomes of plants is performed. A comparison of the methods of targeted sequencing based on the use of the reference DNA sequence is accomplished.

Development of SNP, KASP, and SSR Markers by BSR-Seq Technology for Saturation of Genetic Linkage Map and Efficient Detection of Wheat Powdery Mildew Resistance Gene Pm61

The gene Pm61 that confers powdery mildew resistance has been previously identified on chromosome arm 4AL in Chinese wheat landrace Xuxusanyuehuang (XXSYH). To facilitate the use of Pm61 in breeding practices, the bulked segregant analysis-RNA-Seq (BSR-Seq) analysis, in combination with the information on the Chinese Spring reference genome sequence, was performed in the F_2:3 mapping population of XXSYH × Zhongzuo 9504. Two single nucleotide polymorphism (SNP), two Kompetitive Allele Specific PCR (KASP), and six simple sequence repeat (SSR) markers, together with previously identified polymorphic markers, saturated the genetic linkage map for Pm61, especially in the proximal side of the target gene that was short of gene-linked markers. In the newly established genetic linkage map, Pm61 was located in a 0.71 cM genetic interval and can be detected in a high throughput scale by the KASP markers Xicsk8 and Xicsk13 or by the standard PCR-based markers Xicscx497 and Xicsx538. The newly saturated genetic linkage map will be useful in molecular marker assisted-selection of Pm61 in breeding for disease resistant cultivar and in its map-based cloning.

Developing a High-Throughput SNP-Based Marker System to Facilitate the Introgression of Traits From Aegilops Species Into Bread Wheat (Triticum aestivum)

The genus Aegilops contains a diverse collection of wild species exhibiting variation in geographical distribution, ecological adaptation, ploidy and genome organization. Aegilops is the most closely related genus to Triticum which includes cultivated wheat, a globally important crop that has a limited gene pool for modern breeding. Aegilops species are a potential future resource for wheat breeding for traits, such as adaptation to different ecological conditions and pest and disease resistance. This study describes the development and application of the first high-throughput genotyping platform specifically designed for screening wheat relative species. The platform was used to screen multiple accessions representing all species in the genus Aegilops. Firstly, the data was demonstrated to be useful for screening diversity and examining relationships within and between Aegilops species. Secondly, markers able to characterize and track introgressions from Aegilops species in hexaploid wheat were identified and validated using two different approaches.

Imputation accuracy of wheat genotyping-by-sequencing (GBS) data using barley and wheat genome references

Genotyping-by-sequencing (GBS) provides high SNP coverage and has recently emerged as a popular technology for genetic and breeding applications in bread wheat (Triticum aestivum L.) and many other plant species. Although GBS can discover millions of SNPs, a high rate of missing data is a major concern for many applications. Accurate imputation of those missing data can significantly improve the utility of GBS data. This study compared imputation accuracies among four genome references including three wheat references (Chinese Spring survey sequence, W7984, and IWGSC RefSeq v1.0) and one barley reference genome by comparing imputed data derived from low-depth sequencing to actual data from high-depth sequencing. After imputation, the average number of imputed data points was the highest in the B genome (~48.99%). The D genome had the lowest imputed data points (~15.02%) but the highest imputation accuracy. Among the four reference genomes, IWGSC RefSeq v1.0 reference provided the most imputed data points, but the lowest imputation accuracy for the SNPs with < 10% minor allele frequency (MAF). The W7984 reference, however, provided the highest imputation accuracy for the SNPs with < 10% MAF.

Mapping of novel salt tolerance QTL in an Excalibur × Kukri doubled haploid wheat population

KEY MESSAGE

Novel QTL for salinity tolerance traits have been detected using non-destructive and destructive phenotyping in bread wheat and were shown to be linked to improvements in yield in saline fields. Soil salinity is a major limitation to cereal production. Breeding new salt-tolerant cultivars has the potential to improve cereal crop yields. In this study, a doubled haploid bread wheat mapping population, derived from the bi-parental cross of Excalibur × Kukri, was grown in a glasshouse under control and salinity treatments and evaluated using high-throughput non-destructive imaging technology. Quantitative trait locus (QTL) analysis of this population detected multiple QTL under salt and control treatments. Of these, six QTL were detected in the salt treatment including one for maintenance of shoot growth under salinity (QG(1-5).asl-7A), one for leaf Na+ exclusion (QNa.asl-7A) and four for leaf K+ accumulation (QK.asl-2B.1, QK.asl-2B.2, QK.asl-5A and QK:Na.asl-6A). The beneficial allele for QG(1-5).asl-7A (the maintenance of shoot growth under salinity) was present in six out of 44 mainly Australian bread and durum wheat cultivars. The effect of each QTL allele on grain yield was tested in a range of salinity concentrations at three field sites across 2 years. In six out of nine field trials with different levels of salinity stress, lines with alleles for Na+ exclusion and/or K+ maintenance at three QTL (QNa.asl-7A, QK.asl-2B.2 and QK:Na.asl-6A) excluded more Na+ or accumulated more K+ compared to lines without these alleles. Importantly, the QK.asl-2B.2 allele for higher K+ accumulation was found to be associated with higher grain yield at all field sites. Several alleles at other QTL were associated with higher grain yields at selected field sites.

Genomic approaches for studying crop evolution

Understanding how crop plants evolved from their wild relatives and spread around the world can inform about the origins of agriculture. Here, we review how the rapid development of genomic resources and tools has made it possible to conduct genetic mapping and population genetic studies to unravel the molecular underpinnings of domestication and crop evolution in diverse crop species. We propose three future avenues for the study of crop evolution: establishment of high-quality reference genomes for crops and their wild relatives; genomic characterization of germplasm collections; and the adoption of novel methodologies such as archaeogenetics, epigenomics, and genome editing.

Identification of induced mutations in hexaploid wheat genome using exome capture assay

Wheat is a staple food crop of many countries. Improving resilience to biotic and abiotic stresses remain key breeding targets. Among these, rust diseases are the most detrimental in terms of depressing wheat production. In the present study, chemical mutagenesis was used to induce mutations in the wheat variety NN-Gandum-1. This cultivar is moderately resistant to leaf and yellow rust. The aim of mutagenesis was to improve resistance to the disease as well as to study function of genes conferring resistance to the disease. In the present investigation, a 0.8% EMS dose was found optimum for supporting 45-55% germination of NN-Gandum-1. A total of 3,634 M2 fertile plants were produced from each of the M1 plant. Out of these, 33 (0.91%) and 20 plants (0.55%) showed absolute resistance to leaf and yellow rust, respectively. While 126 (3.46%) and 127 plants (3.49%) exhibited high susceptibility to the leaf and yellow rust, respectively. In the M4 generation, a total of 11 M4 lines (nine absolute resistant and two highly susceptible) and one wild type were selected for NGS-based exome capture assay. A total of 104,779 SNPs were identified that were randomly distributed throughout the wheat sub genomes (A, B and D). Induced mutations in intronic sequences predominated. The highest total number of SNPs detected in this assay were mapped to chr.2B (14,273 SNPs), which contains the highest number of targeted base pairs in the assay. The average mutation density across all regions interrogated was estimated to be one mutation per 20.91 Mb. The highest mutation frequency was found in chr.2D (1/11.7 kb) and the lowest in chr.7D (1/353.4 kb). Out of the detected mutations, 101 SNPs were filtered using analysis criteria aimed to enrich for mutations that may affect gene function. Out of these, one putative SNP detected in Lr21 were selected for further analysis. The SNP identified in chimeric allele (Lr21) of a resistant mutant (N1-252) was located in a NBS domain of chr.1BS at 3.4 Mb position. Through computational analysis, it was demonstrated that this identified SNP causes a substitution of glutamic acid with alanine, resulting in a predicted altered protein structure. This mutation, therefore, is a candidate for contributing to the resistance phenotype in the mutant line. Based on this work, we conclude that the wheat mutant resource developed is useful as a source of novel genetic variation for forward-genetic screens and also as a useful tool for gaining insights into the important biological circuits of different traits of complex genomes like wheat.

Unlocking the novel genetic diversity and population structure of synthetic Hexaploid wheat

Background

Synthetic hexaploid wheat (SHW) is a reconstitution of hexaploid wheat from its progenitors (Triticum turgidum ssp. durum L.; AABB x Aegilops tauschii Coss.; DD) and has novel sources of genetic diversity for broadening the genetic base of elite bread wheat (BW) germplasm (T. aestivum L). Understanding the diversity and population structure of SHWs will facilitate their use in wheat breeding programs. Our objectives were to understand the genetic diversity and population structure of SHWs and compare the genetic diversity of SHWs with elite BW cultivars and demonstrate the potential of SHWs to broaden the genetic base of modern wheat germplasm.

Results

The genotyping-by-sequencing of SHW provided 35,939 high-quality single nucleotide polymorphisms (SNPs) that were distributed across the A (33%), B (36%), and D (31%) genomes. The percentage of SNPs on the D genome was nearly same as the other two genomes, unlike in BW cultivars where the D genome polymorphism is generally much lower than the A and B genomes. This indicates the presence of high variation in the D genome in the SHWs. The D genome gene diversity of SHWs was 88.2% higher than that found in a sample of elite BW cultivars. Population structure analysis revealed that SHWs could be separated into two subgroups, mainly differentiated by geographical location of durum parents and growth habit of the crop (spring and winter type). Further population structure analysis of durum and Ae. parents separately identified two subgroups, mainly based on type of parents used. Although Ae. tauschii parents were divided into two sub-species: Ae. tauschii ssp. tauschii and ssp. strangulate, they were not clearly distinguished in the diversity analysis outcome. Population differentiation between SHWs (Spring_SHW and Winter_SHW) samples using analysis of molecular variance indicated 17.43% of genetic variance between populations and the remainder within populations.

Conclusions

SHWs were diverse and had a clearly distinguished population structure identified through GBS-derived SNPs. The results of this study will provide valuable information for wheat genetic improvement through inclusion of novel genetic variation and is a prerequisite for association mapping and genomic selection to unravel economically important marker-trait associations and for cultivar development.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4969-2) contains supplementary material, which is available to authorized users.

Genomic Selection in Preliminary Yield Trials in a Winter Wheat Breeding Program

Genomic prediction (GP) is now routinely performed in crop plants to predict unobserved phenotypes. The use of predicted phenotypes to make selections is an active area of research. Here, we evaluate GP for predicting grain yield and compare genomic and phenotypic selection by tracking lines advanced. We examined four independent nurseries of F3:6 and F3:7 lines trialed at 6 to 10 locations each year. Yield was analyzed using mixed models that accounted for experimental design and spatial variations. Genotype-by-sequencing provided nearly 27,000 high-quality SNPs. Average genomic predictive ability, estimated for each year by randomly masking lines as missing in steps of 10% from 10 to 90%, and using the remaining lines from the same year as well as lines from other years in a training set, ranged from 0.23 to 0.55. The predictive ability estimated for a new year using the other years ranged from 0.17 to 0.28. Further, we tracked lines advanced based on phenotype from each of the four F3:6 nurseries. Lines with both above average genomic estimated breeding value (GEBV) and phenotypic value (BLUP) were retained for more years compared to lines with either above average GEBV or BLUP alone. The number of lines selected for advancement was substantially greater when predictions were made with 50% of the lines from the testing year added to the training set. Hence, evaluation of only 50% of the lines yearly seems possible. This study provides insights to assess and integrate genomic selection in breeding programs of autogamous crops.

SNP-based pool genotyping and haplotype analysis accelerate fine-mapping of the wheat genomic region containing stripe rust resistance gene Yr26

NGS-assisted super pooling emerging as powerful tool to accelerate gene mapping and haplotype association analysis within target region uncovering specific linkage SNPs or alleles for marker-assisted gene pyramiding. Conventional gene mapping methods to identify genes associated with important agronomic traits require significant amounts of financial support and time. Here, a single nucleotide polymorphism (SNP)-based mapping approach, RNA-Seq and SNP array assisted super pooling analysis, was used for rapid mining of a candidate genomic region for stripe rust resistance gene Yr26 that has been widely used in wheat breeding programs in China. Large DNA and RNA super-pools were genotyped by Wheat SNP Array and sequenced by Illumina HiSeq, respectively. Hundreds of thousands of SNPs were identified and then filtered by multiple filtering criteria. Among selected SNPs, over 900 were found within an overlapping interval of less than 30 Mb as the Yr26 candidate genomic region in the centromeric region of chromosome arm 1BL. The 235 chromosome-specific SNPs were converted into KASP assays to validate the Yr26 interval in different genetic populations. Using a high-resolution mapping population (> 30,000 gametes), we confined Yr26 to a 0.003-cM interval. The Yr26 target region was anchored to the common wheat IWGSC RefSeq v1.0 and wild emmer WEWSeq v.1.0 sequences, from which 488 and 454 kb fragments were obtained. Several candidate genes were identified in the target genomic region, but there was no typical resistance gene in either genome region. Haplotype analysis identified specific SNPs linked to Yr26 and developed robust and breeder-friendly KASP markers. This integration strategy can be applied to accelerate generating many markers closely linked to target genes/QTL for a trait of interest in wheat and other polyploid species.

Conversion of array-based single nucleotide polymorphic markers for use in targeted genotyping by sequencing in hexaploid wheat (Triticum aestivum)

Wheat breeders and academics alike use single nucleotide polymorphisms (SNPs) as molecular markers to characterize regions of interest within the hexaploid wheat genome. A number of SNP-based genotyping platforms are available, and their utility depends upon factors such as the available technologies, number of data points required, budgets and the technical expertise required. Unfortunately, markers can rarely be exchanged between existing and newly developed platforms, meaning that previously generated data cannot be compared, or combined, with more recently generated data sets. We predict that genotyping by sequencing will become the predominant genotyping technology within the next 5-10 years. With this in mind, to ensure that data generated from current genotyping platforms continues to be of use, we have designed and utilized SNP-based capture probes from several thousand existing and publicly available probes from Axiom® and KASP™ genotyping platforms. We have validated our capture probes in a targeted genotyping by sequencing protocol using 31 previously genotyped UK elite hexaploid wheat accessions. Data comparisons between targeted genotyping by sequencing, Axiom® array genotyping and KASP™ genotyping assays, identified a set of 3256 probes which reliably bring together targeted genotyping by sequencing data with the previously available marker data set. As such, these probes are likely to be of considerable value to the wheat community. The probe details, full probe sequences and a custom built analysis pipeline may be freely downloaded from the CerealsDB website (http://www.cerealsdb.uk.net/cerealgenomics/CerealsDB/sequence_capture.php).

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.

Combining Single Nucleotide Polymorphism Genotyping Array with Bulked Segregant Analysis to Map a Gene Controlling Adult Plant Resistance to Stripe Rust in Wheat Line 03031-1-5 H62

Stripe rust, caused by Puccinia striiformis f. sp. tritici, is one of the most devastating diseases of wheat worldwide. Growing resistant cultivars is considered the best approach to manage this disease. In order to identify the resistance gene(s) in wheat line 03031-1-5 H62, which displayed high resistance to stripe rust at adult plant stage, a cross was made between 03031-1-5 H62 and susceptible cultivar Avocet S. The mapping population was tested with Chinese P. striiformis f. sp. tritici race CYR32 through artificial inoculation in a field in Yangling, Shaanxi Province and under natural infection in Tianshui, Gansu Province. The segregation ratios indicated that the resistance was conferred by a single dominant gene, temporarily designated as YrH62. A combination of bulked segregant analysis (BSA) with wheat 90K single nucleotide polymorphism (SNP) array was used to identify molecular markers linked to YrH62. A total of 376 polymorphic SNP loci identified from the BSA analysis were located on chromosome 1B, from which 35 kompetitive allele-specific PCR (KASP) markers selected together with 84 simple sequence repeat (SSR) markers on 1B were used to screen polymorphism and a chromosome region associated with rust resistance was identified. To saturate the chromosomal region covering the YrH62 locus, a 660K SNP array was used to identify more SNP markers. To develop tightly linked markers for marker-assisted selection of YrH62 in wheat breeding, 18 SNPs were converted into KASP markers. A final linkage map consisting of 15 KASP and 3 SSR markers was constructed with KASP markers AX-109352427 and AX-109862469 flanking the YrH62 locus in a 1.0 cM interval. YrH62 explained 63.8 and 69.3% of the phenotypic variation for disease severity and infection type, respectively. YrH62 was located near the centromeric region of chromosome 1BS based on the positions of the SSR markers in 1B deletion bins. Based on the origin, responses to P. striiformis f. sp. tritici races, and marker distances, YrH62 is likely different from the other reported stripe rust resistance genes/quantitative trait loci on 1B. The gene and tightly linked KASP markers will be useful for breeding wheat cultivars with resistance to stripe rust.

High-density SNP mapping reveals closely linked QTL for resistance to Stagonospora nodorum blotch (SNB) in flag leaf and glume of hexaploid wheat

The genetic control of adult plant resistance to Stagonospora nodorum blotch (SNB) is complex, consisting of genes with minor effects interacting in an additive manner. Earlier studies detected quantitative trait loci (QTL) for flag leaf resistance in successive years on chromosomes 1B, 2A, 2D, and 5B using SSR- and DArT-based genetic maps of progeny from the crosses EGA Blanco/Millewa, 6HRWSN125/WAWHT2074, and P92201D5/P91193D1. Similarly, QTL for glume resistance detected in successive years and multiple environments were identified on chromosomes 2D and 4B from genetic maps of P92201D5/P91193D1 and 6HRWSN125/WAWHT2074, respectively. The SSR- and DArT-based genetic maps had an average distance of 6.5, 7.8, and 9.7 cM between marker loci for populations EGA/Millewa, P92201D5/P91193D1, and 6HRWSN125/WAWHT2074, respectively. This study used single nucleotide polymorphism (SNP) markers from the iSelect Infinium 90K genotyping array to fine-map genomic regions harbouring QTL for flag leaf and glume SNB resistance, reducing the average distance between markers to 2.9, 3.3, and 3.4 cM for populations P92201D5/P91193D1, EGA/Millewa, and 6HRWSN125/WAWHT2074, respectively. Increasing the marker density of the genetic maps with SNPs did not identify any new QTL for SNB resistance but discriminated previously identified co-located QTL into separate but closely linked QTL.

Genotyping-by-Sequencing (GBS) Revealed Molecular Genetic Diversity of Iranian Wheat Landraces and Cultivars

Background: Genetic diversity is an essential resource for breeders to improve new cultivars with desirable characteristics. Recently, genotyping-by-sequencing (GBS), a next-generation sequencing (NGS) technology that can simplify complex genomes, has now be used as a high-throughput and cost-effective molecular tool for routine breeding and screening in many crop species, including the species with a large genome. Results: We genotyped a diversity panel of 369 Iranian hexaploid wheat accessions including 270 landraces collected between 1931 and 1968 in different climate zones and 99 cultivars released between 1942 to 2014 using 16,506 GBS-based single nucleotide polymorphism (GBS-SNP) markers. The B genome had the highest number of mapped SNPs while the D genome had the lowest on both the Chinese Spring and W7984 references. Structure and cluster analyses divided the panel into three groups with two landrace groups and one cultivar group, suggesting a high differentiation between landraces and cultivars and between landraces. The cultivar group can be further divided into four subgroups with one subgroup was mostly derived from Iranian ancestor(s). Similarly, landrace groups can be further divided based on years of collection and climate zones where the accessions were collected. Molecular analysis of variance indicated that the genetic variation was larger between groups than within group. Conclusion: Obvious genetic diversity in Iranian wheat was revealed by analysis of GBS-SNPs and thus breeders can select genetically distant parents for crossing in breeding. The diverse Iranian landraces provide rich genetic sources of tolerance to biotic and abiotic stresses, and they can be useful resources for the improvement of wheat production in Iran and other countries.

Crop Breeding Chips and Genotyping Platforms: Progress, Challenges, and Perspectives

There is a rapidly rising trend in the development and application of molecular marker assays for gene mapping and discovery in field crops and trees. Thus far, more than 50 SNP arrays and 15 different types of genotyping-by-sequencing (GBS) platforms have been developed in over 25 crop species and perennial trees. However, much less effort has been made on developing ultra-high-throughput and cost-effective genotyping platforms for applied breeding programs. In this review, we discuss the scientific bottlenecks in existing SNP arrays and GBS technologies and the strategies to develop targeted platforms for crop molecular breeding. We propose that future practical breeding platforms should adopt automated genotyping technologies, either array or sequencing based, target functional polymorphisms underpinning economic traits, and provide desirable prediction accuracy for quantitative traits, with universal applications under wide genetic backgrounds in crops. The development of such platforms faces serious challenges at both the technological level due to cost ineffectiveness, and the knowledge level due to large genotype-phenotype gaps in crop plants. It is expected that such genotyping platforms will be achieved in the next ten years in major crops in consideration of (a) rapid development in gene discovery of important traits, (b) deepened understanding of quantitative traits through new analytical models and population designs, (c) integration of multi-layer -omics data leading to identification of genes and pathways responsible for important breeding traits, and (d) improvement in cost effectiveness of large-scale genotyping. Crop breeding chips and genotyping platforms will provide unprecedented opportunities to accelerate the development of cultivars with desired yield potential, quality, and enhanced adaptation to mitigate the effects of climate change.

Toward Genomics-Based Breeding in C3 Cool-Season Perennial Grasses

Most important food and feed crops in the world belong to the C3 grass family. The future of food security is highly reliant on achieving genetic gains of those grasses. Conventional breeding methods have already reached a plateau for improving major crops. Genomics tools and resources have opened an avenue to explore genome-wide variability and make use of the variation for enhancing genetic gains in breeding programs. Major C3 annual cereal breeding programs are well equipped with genomic tools; however, genomic research of C3 cool-season perennial grasses is lagging behind. In this review, we discuss the currently available genomics tools and approaches useful for C3 cool-season perennial grass breeding. Along with a general review, we emphasize the discussion focusing on forage grasses that were considered orphan and have little or no genetic information available. Transcriptome sequencing and genotype-by-sequencing technology for genome-wide marker detection using next-generation sequencing (NGS) are very promising as genomics tools. Most C3 cool-season perennial grass members have no prior genetic information; thus NGS technology will enhance collinear study with other C3 model grasses like Brachypodium and rice. Transcriptomics data can be used for identification of functional genes and molecular markers, i.e., polymorphism markers and simple sequence repeats (SSRs). Genome-wide association study with NGS-based markers will facilitate marker identification for marker-assisted selection. With limited genetic information, genomic selection holds great promise to breeders for attaining maximum genetic gain of the cool-season C3 perennial grasses. Application of all these tools can ensure better genetic gains, reduce length of selection cycles, and facilitate cultivar development to meet the future demand for food and fodder.

High-throughput SNP genotyping of modern and wild emmer wheat for yield and root morphology using a combined association and linkage analysis

Durum wheat (Triticum turgidum var. durum Desf.) is a major world crop that is grown primarily in areas of the world that experience periodic drought, and therefore, breeding climate-resilient durum wheat is a priority. High-throughput single nucleotide polymorphism (SNP) genotyping techniques have greatly increased the power of linkage and association mapping analyses for bread wheat, but as yet there is no durum wheat-specific platform available. In this study, we evaluate the new 384HT Wheat Breeders Array for its usefulness in tetraploid wheat breeding by genotyping a breeding population of F₆ hybrids, derived from multiple crosses between T. durum cultivars and wild and cultivated emmer wheat accessions. Using a combined linkage and association mapping approach, we generated a genetic map including 1345 SNP markers, and identified markers linked to 6 QTLs for coleoptile length (2), heading date (1), anthocyanin accumulation (1) and osmotic stress tolerance (2). We also developed a straightforward approach for combining genetic data from multiple families of limited size that will be useful for evaluating and mapping pre-existing breeding material.