Genetic architecture of male floral traits required for hybrid wheat breeding

Wheat improvement through advances in single nucleotide polymorphism (SNP) detection and genotyping with a special emphasis on rust resistance

KEY MESSAGE

Single nucleotide polymorphism (SNP) markers in wheat and their prospects in breeding with special reference to rust resistance. Single nucleotide polymorphism (SNP)-based markers are increasingly gaining momentum for screening and utilizing vital agronomic traits in wheat. To date, more than 260 million SNPs have been detected in modern cultivars and landraces of wheat. This rapid SNP discovery was made possible through the release of near-complete reference and pan-genome assemblies of wheat and its wild relatives, coupled with whole genome sequencing (WGS) of thousands of wheat accessions. Further, genotyping customized SNP sites were facilitated by a series of arrays (9 to 820Ks), a cost effective substitute WGS. Lately, germplasm-specific SNP arrays have been introduced to characterize novel traits and detect closely linked SNPs for marker-assisted breeding. Subsequently, the kompetitive allele-specific PCR (KASP) assay was introduced for rapid and large-scale screening of specific SNP markers. Moreover, with the advances and reduction in sequencing costs, ample opportunities arise for generating SNPs artificially through mutations and in combination with next-generation sequencing and comparative genomic analyses. In this review, we provide historical developments and prospects of SNP markers in wheat breeding with special reference to rust resistance where over 50 genetic loci have been characterized through SNP markers. Rust resistance is one of the most essential traits for wheat breeding as new strains of the Puccinia fungus, responsible for rust diseases, evolve frequently and globally.

Beyond pathways: Accelerated flavonoids candidate identification and novel exploration of enzymatic properties using combined mapping populations of wheat

Although forward-genetics-metabolomics methods such as mGWAS and mQTL have proven effective in providing myriad loci affecting metabolite contents, they are somehow constrained by their respective constitutional flaws such as the hidden population structure for GWAS and insufficient recombinant rate for QTL. Here, the combination of mGWAS and mQTL was performed, conveying an improved statistical power to investigate the flavonoid pathways in common wheat. A total of 941 and 289 loci were, respectively, generated from mGWAS and mQTL, within which 13 of them were co-mapped using both approaches. Subsequently, the mGWAS or mQTL outputs alone and their combination were, respectively, utilized to delineate the metabolic routes. Using this approach, we identified two MYB transcription factor encoding genes and five structural genes, and the flavonoid pathway in wheat was accordingly updated. Moreover, we have discovered some rare-activity-exhibiting flavonoid glycosyl- and methyl-transferases, which may possess unique biological significance, and harnessing these novel catalytic capabilities provides potentially new breeding directions. Collectively, we propose our survey illustrates that the forward-genetics-metabolomics approaches including multiple populations with high density markers could be more frequently applied for delineating metabolic pathways in common wheat, which will ultimately contribute to metabolomics-assisted wheat crop improvement.

Green Revolution dwarfing Rht genes negatively affected wheat floral traits related to cross-pollination efficiency

Hybrid breeding is a promising strategy to quickly improve wheat yield and stability. Due to the usefulness of the Rht 'Green Revolution' dwarfing alleles, it is important to gain a better understanding of their impact on traits related to hybrid development. Traits associated with cross-pollination efficiency were studied using Near Isogenic Lines carrying the different sets of alleles in Rht genes: Rht1 (semi-dwarf), Rht2 (semi-dwarf), Rht1 + 2 (dwarf), Rht3 (extreme dwarf), Rht2 + 3 (extreme dwarf), and rht (tall) during four growing seasons. Results showed that the extreme dwarfing alleles Rht2 + 3, Rht3, and Rht1 + 2 presented the greatest effects in all the traits analyzed. Plant height showed reductions up to 64% (Rht2 + 3) compared to rht. Decreases up to 20.2% in anther length and 33% in filament length (Rht2 + 3) were observed. Anthers extrusion decreased from 40% (rht) to 20% (Rht1 and Rht2), 11% (Rht3), 8.3% (Rht1 + 2), and 6.5% (Rht2 + 3). Positive correlations were detected between plant height and anther extrusion, anther, and anther filament lengths, suggesting the negative effect of dwarfing alleles. Moreover, the magnitude of these negative impacts depends on the combination of the alleles: Rht2 + 3 > Rht3/Rht1 + 2 > Rht2/Rht1 > rht (tall). Reductions were consistent across genotypes and environments with interactions due to magnitude effects. Our results indicate that Rht alleles are involved in multiple traits of interest for hybrid wheat production and the need to select alternative sources for reduced height/lodging resistance for hybrid breeding programs.

Genome-Wide Association Study to Identify Marker-Trait Associations for Seed Color in Colored Wheat (Triticum aestivum L.)

This study conducted phenotypic evaluations on a wheat F3 population derived from 155 F2 plants. Traits related to seed color, including chlorophyll a, chlorophyll b, carotenoid, anthocyanin, L*, a*, and b*, were assessed, revealing highly significant correlations among various traits. Genotyping using 81,587 SNP markers resulted in 3969 high-quality markers, revealing a genome-wide distribution with varying densities across chromosomes. A genome-wide association study using fixed and random model circulating probability unification (FarmCPU) and Bayesian-information and linkage-disequilibrium iteratively nested keyway (BLINK) identified 11 significant marker-trait associations (MTAs) associated with L*, a*, and b*, and chromosomal distribution patterns revealed predominant locations on chromosomes 2A, 2B, and 4B. A comprehensive annotation uncovered 69 genes within the genomic vicinity of each MTA, providing potential functional insights. Gene expression analysis during seed development identified greater than 2-fold increases or decreases in expression in colored wheat for 16 of 69 genes. Among these, eight genes, including transcription factors and genes related to flavonoid and ubiquitination pathways, exhibited distinct expression patterns during seed development, providing further approaches for exploring seed coloration. This comprehensive exploration expands our understanding of the genetic basis of seed color and paves the way for informed discussions on the molecular intricacies contributing to this phenotypic trait.

ZmRPN1 confers quantitative variation in pollen number and boosts hybrid seed production in maize

The number of pollen grains is a critical determinant of reproductive success in seed plants and varies among species and individuals. However, in contrast with many mutant-screening studies relevant to anther and pollen development, the natural genetic basis for variations in pollen number remains largely unexplored. To address this issue, we carried out a genome-wide association study in maize, ultimately revealing that a large presence/absence variation in the promoter region of ZmRPN1 alters its expression level and thereby contributes to pollen number variation. Molecular analyses showed that ZmRPN1 interacts with ZmMSP1, which is known as a germline cell number regulator, and facilitates ZmMSP1 localization to the plasma membrane. Importantly, ZmRPN1 dysfunction resulted in a substantial increase in pollen number, consequently boosting seed production by increasing female-male planting ratio. Together, our findings uncover a key gene controlling pollen number, and therefore, modulation of ZmRPN1 expression could be efficiently used to develop elite pollinators for modern hybrid maize breeding.

Population Structure of Modern Winter Wheat Accessions from Central Asia

Despite the importance of winter wheat in Central Asian countries, there are limited reports describing their diversity within this region. In this study, the population structures of 115 modern winter wheat cultivars from four Central Asian countries were compared to germplasms from six other geographic origins using 10,746 polymorphic single-nucleotide polymorphism (SNP) markers. After applying the STRUCTURE package, we found that in terms of the most optimal K steps, samples from Kazakhstan and Kyrgyzstan were grouped together with samples from Russia, while samples from Tajikistan and Uzbekistan were grouped with samples from Afghanistan. The mean value of Nei's genetic diversity index for the germplasm from four groups from Central Asia was 0.261, which is comparable to that of the six other groups studied: Europe, Australia, the USA, Afghanistan, Turkey, and Russia. The Principal Coordinate Analysis (PCoA) showed that samples from Kyrgyzstan, Tajikistan, and Uzbekistan were close to samples from Turkey, while Kazakh accessions were located near samples from Russia. The evaluation of 10,746 SNPs in Central Asian wheat suggested that 1006 markers had opposing allele frequencies. Further assessment of the physical positions of these 1006 SNPs in the Wheat Ensembl database indicated that most of these markers are constituents of genes associated with plant stress tolerance and adaptability. Therefore, the SNP markers identified can be effectively used in regional winter wheat breeding projects for facilitating plant adaptation and stress resistance.

Identification of major QTLs for yield-related traits with improved genetic map in wheat

Introduction

Identification of stable major quantitative trait loci (QTLs) for yield-related traits is important for yield potential improvement in wheat breeding.

Methods

In the present study, we genotyped a recombinant inbred line (RIL) population using the Wheat 660K SNP array and constructed a high-density genetic map. The genetic map showed high collinearity with the wheat genome assembly. Fourteen yield-related traits were evaluated in six environments for QTL analysis.

Results and Discussion

A total of 12 environmentally stable QTLs were identified in at least three environments, explaining up to 34.7% of the phenotypic variation. Of these, QTkw-1B.2 for thousand kernel weight (TKW), QPh-2D.1 (QSl-2D.2/QScn-2D.1) for plant height (PH), spike length (SL) and spikelet compactness (SCN), QPh-4B.1 for PH, and QTss-7A.3 for total spikelet number per spike (TSS) were detected in at least five environments. A set of Kompetitive Allele Specific PCR (KASP) markers were converted based on the above QTLs and used to genotype a diversity panel comprising of 190 wheat accessions across four growing seasons. QPh-2D.1 (QSl-2D.2/QScn-2D.1), QPh-4B.1 and QTss-7A.3 were successfully validated. Compared with previous studies, QTkw-1B.2 and QPh-4B.1 should be novel QTLs. These results provided a solid foundation for further positional cloning and marker-assisted selection of the targeted QTLs in wheat breeding programs.

Genome-Wide Association Analysis for Hybrid Breeding in Wheat

Disclosure of markers that are significantly associated with plant traits can help develop new varieties with desirable properties. This study determined the genome-wide associations based on DArTseq markers for six agronomic traits assessed in eight environments for wheat. Moreover, the association study for heterosis and analysis of the effects of markers grouped by linkage disequilibrium were performed based on mean values over all experiments. All results were validated using data from post-registration trials. GWAS revealed 1273 single nucleotide polymorphisms with biologically significant effects. Most polymorphisms were predicted to be modifiers of protein translation, with only two having a more pronounced effect. Markers significantly associated with the considered set of features were clustered within chromosomes based on linkage disequilibrium in 327 LD blocks. A GWAS for heterosis revealed 1261 markers with significant effects.

Domestication of newly evolved hexaploid wheat—A journey of wild grass to cultivated wheat

Domestication of wheat started with the dawn of human civilization. Since then, improvement in various traits including resistance to diseases, insect pests, saline and drought stresses, grain yield, and quality were improved through selections by early farmers and then planned hybridization after the discovery of Mendel’s laws. In the 1950s, genetic variability was created using mutagens followed by the selection of superior mutants. Over the last 3 decades, research was focused on developing superior hybrids, initiating marker-assisted selection and targeted breeding, and developing genetically modified wheat to improve the grain yield, tolerance to drought, salinity, terminal heat and herbicide, and nutritive quality. Acceptability of genetically modified wheat by the end-user remained a major hurdle in releasing into the environment. Since the beginning of the 21st century, changing environmental conditions proved detrimental to achieving sustainability in wheat production particularly in developing countries. It is suggested that high-tech phenotyping assays and genomic procedures together with speed breeding procedures will be instrumental in achieving food security beyond 2050.

Overexpression of the WAPO-A1 gene increases the number of spikelets per spike in bread wheat

Two homoeologous QTLs for number of spikelets per spike (SPS) were mapped on chromosomes 7AL and 7BL using two wheat MAGIC populations. Sets of lines contrasting for the QTL on 7AL were developed which allowed for the validation and fine mapping of the 7AL QTL and for the identification of a previously described candidate gene, WHEAT ORTHOLOG OF APO1 (WAPO1). Using transgenic overexpression in both a low and a high SPS line, we provide a functional validation for the role of this gene in determining SPS also in hexaploid wheat. We show that the expression levels of this gene positively correlate with SPS in multiple MAGIC founder lines under field conditions as well as in transgenic lines grown in the greenhouse. This work highlights the potential use of WAPO1 in hexaploid wheat for further yield increases. The impact of WAPO1 and SPS on yield depends on other genetic and environmental factors, hence, will require a finely balanced expression level to avoid the development of detrimental pleiotropic phenotypes.

Genome-Wide Association Study of Leaf Rust and Stem Rust Seedling and Adult Resistances in Tetraploid Wheat Accessions Harvested in Kazakhstan

Leaf rust (LR) and stem rust (SR) are diseases increasingly impacting wheat production worldwide. Fungal pathogens producing rust diseases in wheat may cause yield losses of up to 50−60%. One of the most effective methods for preventing such losses is the development of resistant cultivars with high yield potential. This goal can be achieved through complex breeding studies, including the identification of key genetic factors controlling rust disease resistance. The objective of this study was to identify sources of tetraploid wheat resistance to LR and SR races, both at the seedling growth stage in the greenhouse and at the adult plant stage in field experiments, under the conditions of the North Kazakhstan region. A panel consisting of 193 tetraploid wheat accessions was used in a genome-wide association study (GWAS) for the identification of quantitative trait loci (QTLs) associated with LR and SR resistance, using 16,425 polymorphic single-nucleotide polymorphism (SNP) markers in the seedling and adult stages of plant development. The investigated panel consisted of seven tetraploid subspecies (Triticum turgidum ssp. durum, ssp. turanicum, ssp. turgidum, ssp. polonicum, ssp. carthlicum, ssp. dicoccum, and ssp. dicoccoides). The GWAS, based on the phenotypic evaluation of the tetraploid collection’s reaction to the two rust species at the seedling (in the greenhouse) and adult (in the field) stages, revealed 38 QTLs (p < 0.001), comprising 17 for LR resistance and 21 for SR resistance. Ten QTLs were associated with the reaction to LR at the seedling stage, while six QTLs were at the adult plant stage and one QTL was at both the seedling and adult stages. Eleven QTLs were associated with SR response at the seedling stage, while nine QTLs were at the adult plant stage and one QTL was at both the seedling and adult stages. A comparison of these results with previous LR and SR studies indicated that 11 of the 38 QTLs are presumably novel loci. The QTLs identified in this work can potentially be used for marker-assisted selection of tetraploid and hexaploid wheat for the breeding of new LR- and SR-resistant cultivars.

The effect of the Rht1 haplotype on Fusarium head blight resistance in relation to type and level of background resistance and in combination with Fhb1 and Qfhs.ifa-5A

The effect of the Rht1-genes on FHB resistance depends on anther extrusion and level of background resistance. Qfhs.ifa-5A increases resistance and anther extrusion as efficiently as semi-dwarfing alleles decrease it. The semi-dwarfing reduced height alleles Rht-D1b and Rht-B1b have been deployed in modern wheat cultivars throughout the world, but they increase susceptibility to Fusarium head blight (FHB). Here, we investigated the impact of the Rht1 genes on anther retention (AR) in relation to FHB resistance using four different sets of near-isogenic lines (NILs) with contrasting levels and types of background FHB resistance. NILs were evaluated for FHB severity, plant height and AR in three greenhouse and three field trials using artificial spray inoculation. Rht-B1b and Rht-D1b alleles increased AR and FHB susceptibility in all genetic backgrounds. The magnitude of the effects differed between NIL groups. Increased FHB susceptibility largely followed increased AR. Differences in FHB susceptibility between tall and dwarf haplotypes were largest in the NIL group with the highest changes in AR. In the most resistant NIL group, dwarfed lines had only slightly higher AR than tall lines and maintained good resistance, while both tall and dwarf lines had high levels of retained anthers in the most susceptible NIL group. We further investigated the effect of the major Fusarium resistance QTL Fhb1 and Qfhs.ifa-5A in combination with the Rht1 genes. Qfhs.ifa-5A enhanced anther extrusion in tall as well as semi-dwarf haplotypes, whereas Fhb1 did not affect AR. Qfhs.ifa-5A supported FHB resistance more efficiently than Fhb1 in lines that were more responsive to AR, while both Fhb1 and Qfhs.ifa-5A were equally efficient in NILs that had high background resistance and low response to AR.

Genetic Incorporation of the Favorable Alleles for Three Genes Associated With Spikelet Development in Wheat

The number of spikelets per spike is an important trait that directly affects grain yield in wheat. Three quantitative trait loci (QTLs) associated with spikelet nodes per spike (SNS) were mapped in a population of recombinant inbred lines generated from a cross between two advanced breeding lines of winter wheat based on the phenotypic variation evaluated over six locations/years. Two of the three QTLs are QSns.sxau-2A at the WHEATFRIZZY PANICLE (WFZP) loci and QSns.sxau-7A at the WHEAT ORTHOLOG OF APO1 (WAPO1) loci. The WFZP-A1b allele with a 14-bp deletion at QSns.sxau-2A was associated with increased spikelets per spike. WAPO-A1e, as a novel allele at WAPO1, were regulated at the transcript level that was associated with the SNS trait. The third SNS QTL, QSns.sxau-7D on chromosome 7D, was not associated with homoeologous WAPO-D1 or any other genes known to regulate SNS. The favorable alleles for each of WZFP-A1, WAPO-A1, and QSns.sxau-7D are identified and incorporated to increase up to 3.4 spikelets per spike in the RIL lines. Molecular markers for the alleles were developed. This study has advanced our understanding of the genetic basis of natural variation in spikelet development in wheat.

Induction of Semi-Dwarf Trait to a Three Pistil Tall Mutant Through Breeding With Increased Grain Numbers in Wheat

Multi-ovary wheat (three pistil) is a unique germplasm for the seed production of hybrid wheat. The purpose of the present study was to transfer the multi-ovary trait to semi-dwarf plants to increase the production of grains in wheat crops. Therefore, tall, semi-dwarf, and dwarf plants were crossed with plants with the three-pistil trait. A three-pistil tall plant was used as the female parent, while tall (Synthetic hexaploid), semi-dwarf, and dwarf plants were used as male parents. F1 and F2 progenies with parents were planted in 2015-16 using RCBD. The outcome of the crosses showed that multi-ovary tall plants gave significant difference for all five traits (days to maturity, plant height, number of seeds per spike, grain weight per spike, and grain yield per unit area) in both generations. The greatest number of grains per spike and grain yield per unit area were obtained from the cross of three-pistil tall and dwarf parent (P1/P6) in the F1 and F2 generations. The cross also resulted in a significant reduction in height (96 cm). Further heterosis studies conducted with crosses between three-pistil tall and dwarf parent (P1/P6) showed the greatest heterosis and heterobeltiosis for the number of grains per spike (60.0 and 26.19%, respectively) and grain yield per m² (27.68 and 2.85%, respectively). In the case of grain weight per spike, the heterosis value was also positive and significant (17.7). Meanwhile, for other traits, their values were negative for heterosis and heterobeltiosis. High numbers of grains and grain weight were found to be associated with positive heterobeltiosis and in turn the grain yield per m², but plant height and maturity had negative affirmation with heterobeltiosis. Heterosis studies also indicated the dominant gene action for the three-pistil trait. Thus, the study clearly signified that grain yield can be increased by using the multi-ovary genotype with the semi-dwarf height. This new germplasm will be helpful for breeders to increase the production of wheat crops in the southern climate of Pakistan.

Hybrid Seed Set in Relation with Male Floral Traits, Estimation of Heterosis and Combining Abilities for Yield and Its Components in Wheat (Triticum aestivum L.)

Breeding hybrids with maximum heterosis requires efficient cross-pollination and an improved male sterility system. Renewed efforts have been made to dissect the phenotypic variation and genetic basis of hybrid floral traits, although the potential of tailoring the appropriate flower design on seed setting is less known. To this end, elite wheat genotypes were crossed using a chemical hybridizing agent at different doses. A total of 23 hybrids were developed from a partial diallel design; and planted in an alpha lattice design with their parents at two locations in Morocco, for two years, to evaluate for yield components, heterosis and combining abilities. The 13.5 L ha^-1 dose induced a maximum level of sterility (95%) and seed set showed large phenotypic variation and high heritability. In parallel, seed set showed tight correlation with pollen mass (0.97), visual anther extrusion (0.94) and pollen shedding (0.91) (p < 0.001), allowing direct selection of the associated traits. Using the combined data, mid-parent heterosis ranges were -7.64-14.55% for biomass (BM), -8.34-12.51% for thousand kernel weight (TKW) and -5.29-26.65% for grain yield (YLD); while best-parent heterosis showed ranges of -11.18-7.20%, -11.35-11.26% and -8.27-24.04% for BM, TKW and YLD, respectively. The magnitude of general combining ability (GCA) variance was greater than the specific combining ability (SCA) variance suggesting a greater additive gene action for BM, TKW and YLD. The favorable GCA estimates showed a simple method to predict additive effects contributing to high heterosis and thus could be an effective approach for the selection of promising parents in early generations.

Pollen Number and Ribosome Gene Expression Altered in a Genome-Editing Mutant of REDUCED POLLEN NUMBER1 Gene

The number of pollen grains varies within and between species. However, little is known about the molecular basis of this quantitative trait, in contrast with the many studies available on cell differentiation in the stamen. Recently, the first gene responsible for pollen number variation, REDUCED POLLEN NUMBER1 (RDP1), was isolated by genome-wide association studies of Arabidopsis thaliana and exhibited the signature of natural selection. This gene encodes a homolog of yeast Mrt4 (mRNA turnover4), which is an assembly factor of the large ribosomal subunit. However, no further data were available to link ribosome function to pollen development. Here, we characterized the RDP1 gene using the standard A. thaliana accession Col-0. The frameshift mutant, rdp1-3 generated by CRISPR/Cas9 revealed the pleiotropic effect of RDP1 in flowering, thus demonstrating that this gene is required for a broad range of processes other than pollen development. We found that the natural Col-0 allele conferred a reduced pollen number against the Bor-4 allele, as assessed using the quantitative complementation test, which is more sensitive than transgenic experiments. Together with a historical recombination event in Col-0, which was identified by sequence alignment, these results suggest that the coding sequence of RDP1 is the candidate region responsible for the natural phenotypic variation. To elucidate the biological processes in which RDP1 is involved, we conducted a transcriptome analysis. We found that genes responsible for ribosomal large subunit assembly/biogenesis were enriched among the differentially regulated genes, which supported the hypothesis that ribosome biogenesis is disturbed in the rdp1-3 mutant. Among the pollen-development genes, three key genes encoding basic helix-loop-helix (bHLH) transcription factors (ABORTED MICROSPORES (AMS), bHLH010, and bHLH089), as well as direct downstream genes of AMS, were downregulated in the rdp1-3 mutant. In summary, our results suggest a specialized function of ribosomes in pollen development through RDP1, which harbors natural variants under selection.

WAPO-A1 is the causal gene of the 7AL QTL for spikelet number per spike in wheat

Improving our understanding of the genes regulating grain yield can contribute to the development of more productive wheat varieties. Previously, a highly significant QTL affecting spikelet number per spike (SNS), grain number per spike (GNS) and grain yield was detected on chromosome arm 7AL in multiple genome-wide association studies. Using a high-resolution genetic map, we established that the A-genome homeolog of WHEAT ORTHOLOG OF APO1 (WAPO-A1) was a leading candidate gene for this QTL. Using mutants and transgenic plants, we demonstrate in this study that WAPO-A1 is the causal gene underpinning this QTL. Loss-of-function mutants wapo-A1 and wapo-B1 showed reduced SNS in tetraploid wheat, and the effect was exacerbated in wapo1 combining both mutations. By contrast, spikes of transgenic wheat plants carrying extra copies of WAPO-A1 driven by its native promoter had higher SNS, a more compact spike apical region and a smaller terminal spikelet than the wild type. Taken together, these results indicate that WAPO1 affects SNS by regulating the timing of terminal spikelet formation. Both transgenic and wapo1 mutant plants showed a wide range of floral abnormalities, indicating additional roles of WAPO1 on wheat floral development. Previously, we found three widespread haplotypes in the QTL region (H1, H2 and H3), each associated with particular WAPO-A1 alleles. Results from this and our previous study show that the WAPO-A1 allele in the H1 haplotype (115-bp deletion in the promoter) is expressed at significantly lower levels in the developing spikes than the alleles in the H2 and H3 haplotypes, resulting in reduced SNS. Field experiments also showed that the H2 haplotype is associated with the strongest effects in increasing SNS and GNS (H2>H3>H1). The H2 haplotype is already present in most modern common wheat varieties but is rare in durum wheat, where it might be particularly useful to improve grain yield.

Comprehensive evaluation of mapping complex traits in wheat using genome-wide association studies

Genome-wide association studies (GWAS) are effectively applied to detect the marker trait associations (MTAs) using whole genome-wide variants for complex quantitative traits in different crop species. GWAS has been applied in wheat for different quality, biotic and abiotic stresses, and agronomic and yield-related traits. Predictions for marker-trait associations are controlled with the development of better statistical models taking population structure and familial relatedness into account. In this review, we have provided a detailed overview of the importance of association mapping, population design, high-throughput genotyping and phenotyping platforms, advancements in statistical models and multiple threshold comparisons, and recent GWA studies conducted in wheat. The information about MTAs utilized for gene characterization and adopted in breeding programs is also provided. In the literature that we surveyed, as many as 86,122 wheat lines have been studied under various GWA studies reporting 46,940 loci. However, further utilization of these is largely limited. The future breakthroughs in area of genomic selection, multi-omics-based approaches, machine, and deep learning models in wheat breeding after exploring the complex genetic structure with the GWAS are also discussed. This is a most comprehensive study of a large number of reports on wheat GWAS and gives a comparison and timeline of technological developments in this area. This will be useful to new researchers or groups who wish to invest in GWAS.

Genome-Wide Association Analysis of Stable Stripe Rust Resistance Loci in a Chinese Wheat Landrace Panel Using the 660K SNP Array

Stripe rust (caused by Puccinia striiformis f. sp. tritici) is one of the most severe diseases affecting wheat production. The disease is best controlled by developing and growing resistant cultivars. Chinese wheat (Triticum aestivum) landraces have excellent resistance to stripe rust. The objectives of this study were to identify wheat landraces with stable resistance and map quantitative trait loci (QTL) for resistance to stripe rust from 271 Chinese wheat landraces using a genome-wide association study (GWAS) approach. The landraces were phenotyped for stripe rust responses at the seedling stage with two predominant Chinese races of P. striiformis f. sp. tritici in a greenhouse and the adult-plant stage in four field environments and genotyped using the 660K wheat single-nucleotide polymorphism (SNP) array. Thirteen landraces with stable resistance were identified, and 17 QTL, including eight associated to all-stage resistance and nine to adult-plant resistance, were mapped on chromosomes 1A, 1B, 2A, 2D, 3A, 3B, 5A, 5B, 6D, and 7A. These QTL explained 6.06-16.46% of the phenotypic variation. Five of the QTL, QYrCL.sicau-3AL, QYrCL.sicau-3B.4, QYrCL.sicau-3B.5, QYrCL.sicau-5AL.1 and QYrCL.sicau-7AL, were likely new. Five Kompetitive allele specific PCR (KASP) markers for four of the QTL were converted from the significant SNP markers. The identified wheat landraces with stable resistance to stripe rust, significant QTL, and KASP markers should be useful for breeding wheat cultivars with durable resistance to stripe rust.

Training set design in genomic prediction with multiple biparental families

Genomic selection is a powerful tool to reduce the cycle length and enhance the genetic gain of complex traits in plant breeding. However, questions remain about the optimum design and composition of the training set. In this study, we used 944 soybean [Glycine max (L.) Merr.] recombinant inbred lines from eight families derived through a partial-diallel mating design among five parental lines. The cross-validated prediction accuracies for the six traits seed yield, 1,000-seed weight, protein yield, plant height, protein content, and oil content were high, ranging from 0.79 to 0.87. We investigated among-family predictions, making use of the special mating design with different degrees of relatedness among families. Generally, the prediction accuracy decreased from full-sibs to half-sib families to unrelated families. However, half-sib and unrelated families also showed substantial variation in their prediction accuracy for a given family, which appeared to be caused at least in part by the shared segregation of quantitative trait loci in both the training and prediction sets. Combining several half-sib families in composite training sets generally led to an increase in the prediction accuracy compared with the best family alone. The prediction accuracy increased with the size of the training set, but for comparable prediction accuracy, substantially more half-sibs were required than full-sibs. Collectively, our results highlight the potential of genomic selection for soybean breeding and, in a broader context, illustrate the importance of the targeted design of the training set.

Identification of genetic loci associated with major agronomic traits of wheat (Triticum aestivum L.) based on genome-wide association analysis

BACKGROUND

Bread wheat (Triticum aestivum L.) is one of the most widely consumed cereal crops, but its complex genome makes it difficult to investigate the genetic effect on important agronomic traits. Genome-wide association (GWA) analysis is a useful method to identify genetic loci controlling complex phenotypic traits. With the RNA-sequencing based gene expression analysis, putative candidate genes governing important agronomic trait can be suggested and also molecular markers can be developed.

RESULTS

We observed major quantitative agronomic traits of wheat; the winter survival rate (WSR), days to heading (DTH), days to maturity (DTM), stem length (SL), spike length (SPL), awn length (AL), liter weight (LW), thousand kernel weight (TKW), and the number of seeds per spike (SPS), of 287 wheat accessions from diverse country origins. A significant correlation was observed between the observed traits, and the wheat genotypes were divided into three subpopulations according to the population structure analysis. The best linear unbiased prediction (BLUP) values of the genotypic effect for each trait under different environments were predicted, and these were used for GWA analysis based on a mixed linear model (MLM). A total of 254 highly significant marker-trait associations (MTAs) were identified, and 28 candidate genes closely located to the significant markers were predicted by searching the wheat reference genome and RNAseq data. Further, it was shown that the phenotypic traits were significantly affected by the accumulation of favorable or unfavorable alleles.

CONCLUSIONS

From this study, newly identified MTA and putative agronomically useful genes will help to study molecular mechanism of each phenotypic trait. Further, the agronomically favorable alleles found in this study can be used to develop wheats with superior agronomic traits.

Pre-emptive Breeding Against Karnal Bunt Infection in Common Wheat: Combining Genomic and Agronomic Information to Identify Suitable Parents

Wheat (Triticum aestivum L.) is the most widely grown cereal crop in the world and is staple food to half the world's population. The current world population is expected to reach 9.8 billion people by 2050, but food production is not expected to keep pace with demand in developing countries. Significant opportunities exist for traditional grain exporters to produce and export greater amounts of wheat to fill the gap. Karnal bunt, however, is a major threat, due to its use as a non-tariff trade barrier by several wheat-importing countries. The cultivation of resistant varieties remains the most cost-effective approach to manage the disease, but in countries that are free of the disease, genetic improvement is difficult due to quarantine restrictions. Here we report a study on pre-emptive breeding designed to identify linked molecular markers, evaluate the prospects of genomic selection as a tool, and prioritise wheat genotypes suitable for use as parents. In a genome-wide association (GWAS) study, we identified six DArTseq markers significantly linked to Karnal bunt resistance, which explained between 7.6 and 29.5% of the observed phenotypic variation. The accuracy of genomic prediction was estimated to vary between 0.53 and 0.56, depending on whether it is based solely on the identified Quantitative trait loci (QTL) markers or the use of genome-wide markers. As genotypes used as parents would be required to possess good yield and phenology, further research was conducted to assess the agronomic value of Karnal bunt resistant germplasm from the International Maize and Wheat Improvement Center (CIMMYT). We identified an ideal genotype, ZVS13_385, which possessed similar agronomic attributes to the highly successful Australian wheat variety, Mace. It is phenotypically resistant to Karnal bunt infection (<1% infection) and carried all the favourable alleles detected for resistance in this study. The identification of a genotype combining Karnal bunt resistance with adaptive agronomic traits overcomes the concerns of breeders regarding yield penalty in the absence of the disease.

Genome-wide association study of yield components in spring wheat collection harvested under two water regimes in Northern Kazakhstan

Background

Bread wheat is the most important cereal in Kazakhstan, where it is grown on over 12 million hectares. One of the major constraints affecting wheat grain yield is drought due to the limited water supply. Hence, the development of drought-resistant cultivars is critical for ensuring food security in this country. Therefore, identifying quantitative trait loci (QTLs) associated with drought tolerance as an essential step in modern breeding activities, which rely on a marker-assisted selection approach.

Methods

A collection of 179 spring wheat accessions was tested under irrigated and rainfed conditions in Northern Kazakhstan over three years (2018, 2019, and 2020), during which data was collected on nine traits: heading date (HD), seed maturity date (SMD), plant height (PH), peduncle length (PL), number of productive spikes (NPS), spike length (SL), number of kernels per spike (NKS), thousand kernel weight (TKW), and kernels yield per m² (YM2). The collection was genotyped using a 20,000 (20K) Illumina iSelect SNP array, and 8,662 polymorphic SNP markers were selected for a genome-wide association study (GWAS) to identify QTLs for targeted agronomic traits.

Results

Out of the total of 237 discovered QTLs, 50 were identified as being stable QTLs for irrigated and rainfed conditions in the Akmola region, Northern Kazakhstan; the identified QTLs were associated with all the studied traits except PH. The results indicate that nine QTLs for HD and 11 QTLs for SMD are presumably novel genetic factors identified in the irrigated and rainfed conditions of Northern Kazakhstan. The identified SNP markers of the QTLs for targeted traits in rainfed conditions can be applied to develop new competitive spring wheat cultivars in arid zones using a marker-assisted selection approach.

Exploring the genic resources underlying metabolites through mGWAS and mQTL in wheat: From large-scale gene identification and pathway elucidation to crop improvement

Common wheat (Triticum aestivum L.) is a leading cereal crop, but has lagged behind with respect to the interpretation of the molecular mechanisms of phenotypes compared with other major cereal crops such as rice and maize. The recently available genome sequence of wheat affords the pre-requisite information for efficiently exploiting the potential molecular resources for decoding the genetic architecture of complex traits and identifying valuable breeding targets. Meanwhile, the successful application of metabolomics as an emergent large-scale profiling methodology in several species has demonstrated this approach to be accessible for reaching the above goals. One such productive avenue is combining metabolomics approaches with genetic designs. However, this trial is not as widespread as that for sequencing technologies, especially when the acquisition, understanding, and application of metabolic approaches in wheat populations remain more difficult and even arguably underutilized. In this review, we briefly introduce the techniques used in the acquisition of metabolomics data and their utility in large-scale identification of functional candidate genes. Considerable progress has been made in delivering improved varieties, suggesting that the inclusion of information concerning these metabolites and genes and metabolic pathways enables a more explicit understanding of phenotypic traits and, as such, this procedure could serve as an -omics-informed roadmap for executing similar improvement strategies in wheat and other species.

Genome-Wide Association and Prediction of Male and Female Floral Hybrid Potential Traits in Elite Spring Bread Wheat Genotypes

Hybrid wheat breeding is one of the most promising technologies for further sustainable yield increases. However, the cleistogamous nature of wheat displays a major bottleneck for a successful hybrid breeding program. Thus, an optimized breeding strategy by developing appropriate parental lines with favorable floral trait combinations is the best way to enhance the outcrossing ability. This study, therefore, aimed to dissect the genetic basis of various floral traits using genome-wide association study (GWAS) and to assess the potential of genome-wide prediction (GP) for anther extrusion (AE), visual anther extrusion (VAE), pollen mass (PM), pollen shedding (PSH), pollen viability (PV), anther length (AL), openness of the flower (OPF), duration of floret opening (DFO) and stigma length. To this end, we employed 196 ICARDA spring bread wheat lines evaluated for three years and genotyped with 10,477 polymorphic SNP. In total, 70 significant markers were identified associated to the various assessed traits at FDR ≤ 0.05 contributing a minor to large proportion of the phenotypic variance (8–26.9%), affecting the traits either positively or negatively. GWAS revealed multi-marker-based associations among AE, VAE, PM, OPF and DFO, most likely linked markers, suggesting a potential genomic region controlling the genetic association of these complex traits. Of these markers, Kukri_rep_c103359_233 and wsnp_Ex_rep_c107911_91350930 deserve particular attention. The consistently significant markers with large effect could be useful for marker-assisted selection. Genomic selection revealed medium to high prediction accuracy ranging between 52% and 92% for the assessed traits with the least and maximum value observed for stigma length and visual anther extrusion, respectively. This indicates the feasibility to implement genomic selection to predict the performance of hybrid floral traits with high reliability.

Factors Affecting the Number of Pollen Grains per Male Strobilus in Japanese Cedar (Cryptomeria japonica)

Japanese cedar (Cryptomeria japonica) is the most important timber species in Japan; however, its pollen is the primary cause of pollinosis in Japan. The total number of pollen grains produced by a single tree is determined by the number of male strobili (male flowers) and the number of pollen grains per male strobilus. While the number of male strobili is a visible and well-investigated trait, little is known about the number of pollen grains per male strobilus. We hypothesized that genetic and environmental factors affect the pollen number per male strobilus and explored the factors that affect pollen production and genetic variation among clones. We counted pollen numbers of 523 male strobili from 26 clones using a cell counter method that we recently developed. Piecewise Structural Equation Modeling (pSEM) revealed that the pollen number is mostly affected by genetic variation, male strobilus weight, and pollen size. Although we collected samples from locations with different environmental conditions, statistical modeling succeeded in predicting pollen numbers for different clones sampled from branches facing different directions. Comparison of predicted pollen numbers revealed that they varied >3-fold among the 26 clones. The determination of the factors affecting pollen number and a precise evaluation of genetic variation will contribute to breeding strategies to counter pollinosis. Furthermore, the combination of our efficient counting method and statistical modeling will provide a powerful tool not only for Japanese cedar but also for other plant species.

Spikelet movements, anther extrusion and pollen production in wheat cultivars with contrasting tendencies to cleistogamy

BACKGROUND

Cleistogamic flowers are a main barrier in pollen dispersal for cross-pollination necessary in wheat hybrid breeding. The aim of our study was to gain new knowledge on the biology of wheat flowering, in particular on the differences between the cleisto- and chasmogamic forms which has certainly cognitive significance, but it can also be used in practice when seeking a female and male ideotypes for cross hybridization.

RESULTS

We characterized the most significant features defining the flowering specificity in two wheat cultivars with contrasting tendency to cleistogamy: Piko (chasmogamous) and Dacanto (cleistogamous). In the field observations we assessed diurnal pattern of anther extrusion and anther extrusion capacity. For the first time we adapted the time lapse method for measuring kinetics of the spikelet movement and 3-D image correlation technique for the non-invasive measurements of potential deformations of the spikelet lemmas. We found that the two cultivars differ in the potential of pollen dispersion for-cross-pollination and in the spikelet kinetics. We also described some anatomical traits that can have potential functional role in floret opening. None of the cultivars showed any symptoms of lemma surface deformation.

CONCLUSIONS

The cleistogamic and chasmogamic wheat cultivars differ significantly in the potential for pollen dispersion for cross-pollination, which is mainly related to anther extrusion capacity. Although none of these features differentiated the cultivars clearly, we assume, based on spikelet kinetics and the lack of lemmas surface deformation, that the water transport and turgor of cells is essential for the floret opening and anther extrusion in wheat. The search for parental ideotype should be supported by marker assisted selection, e.g. based of polymorphisms in genes related to aquaporin biosynthesis.

Identification of quantitative trait loci for Fusarium head blight (FHB) resistance in the cross between wheat landrace N553 and elite cultivar Yangmai 13

Fusarium head blight (FHB) of wheat poses a serious threat to food security in the Yellow-Huai River Valley Winter Wheat Region (YHW) of China. Discovery of new resistant quantitative trait loci (QTLs) or genes and application of them to highly susceptible varieties in the YHW are of great significance for ensuring the grain yield. Here, 160 recombinant inbred lines (RILs) from the cross between N553 (resistant) and Yangmai 13 (moderately susceptible) were used to evaluate FHB resistance by point inoculation, spray inoculation, and natural infection. A high-density genetic map was constructed by using a 15K SNP array and 128 polymorphism SSR markers. A total of 1452 polymorphic markers were identified, which formed 21 linkage groups and covered a total of 3555.1 cM in length. Two and four QTLs respectively related to type I and type II resistance were detected, among which QFhb-hnau.3BS.1 and QFhb-hnau.2DL were stably identified in most environments in Yangzhou and Zhengzhou, whereas QFhbn-hnau.5AL was only identified under natural infection in Jianyang. Based on the physical position (IWGSC RefSeq v1.0), QFhb-hnau.3BS.1 from the landrace N553 is likely to be Fhb1, while QFhb-hnau.2DL from Yangmai 13 may be a novel QTL. Significantly higher FHB resistance was observed in the lines with both QFhb-hnau.3BS.1 and QFhb-hnau.2DL, indicating that these two QTLs have apparent additive effects, and the RILs harboring both the two QTLs may have great application potential for the improvement of FHB resistance in wheat breeding.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01220-5.

Genetic factors associated with favourable pollinator traits in the wheat cultivar Piko

Hybrid breeding in wheat has the potential to boost yields. An efficient hybrid seed production system requires elite pollinators; however, such germplasm is limited among modern cultivars. Piko, a winter wheat (Triticum aestivum L.) cultivar, has been identified as a superior pollinator and has been used in Europe. Piko has favourable pollinator traits for anther extrusion, anther length, pollen mass and hybrid seed set. However, the genetic factors responsible for Piko's favourable traits are largely unknown. Here, we report on the genetic analysis of a Piko-derived F2 mapping population. We confirmed that Piko's Rht-D1a allele for tall stature is associated with large anthers and high anther extrusion. However, Rht-D1 was not found to be associated with anther filament length, confirmed by near isogenic lines. Piko's photoperiod sensitive Ppd-B1b allele shows an association with increased spike length, more spikelets and spike architecture traits, while the insensitive Ppd-B1a allele is linked with high anther extrusion and larger anthers. We identified an anther extrusion quantitative trait locus (QTL) on chromosome 6A that showed significantly biased transmission of the favourable Piko allele amongst F2 progenies. The Piko allele is completely absent in the distal 6AS region and the central 6A region revealed a significantly lower ratio (<8%) of F2 with homozygous Piko alleles. Our study provided further evidence for the effects of Rht-D1 and Ppd-B1 loci on multiple pollinator traits and a novel anther extrusion QTL that exhibits segregation distortion.

From Genetic Maps to QTL Cloning: An Overview for Durum Wheat

Durum wheat is one of the most important cultivated cereal crops, providing nutrients to humans and domestic animals. Durum breeding programs prioritize the improvement of its main agronomic traits; however, the majority of these traits involve complex characteristics with a quantitative inheritance (quantitative trait loci, QTL). This can be solved with the use of genetic maps, new molecular markers, phenotyping data of segregating populations, and increased accessibility to sequences from next-generation sequencing (NGS) technologies. This allows for high-density genetic maps to be developed for localizing candidate loci within a few Kb in a complex genome, such as durum wheat. Here, we review the identified QTL, fine mapping, and cloning of QTL or candidate genes involved in the main traits regarding the quality and biotic and abiotic stresses of durum wheat. The current knowledge on the used molecular markers, sequence data, and how they changed the development of genetic maps and the characterization of QTL is summarized. A deeper understanding of the trait architecture useful in accelerating durum wheat breeding programs is envisioned.

Intercontinental trials reveal stable QTL for Northern corn leaf blight resistance in Europe and in Brazil

KEY MESSAGE

NCLB is the most devastating leaf disease in European maize, and the introduction of Brazilian resistance donors can efficiently increase the resistance levels of European maize germplasm. Northern corn leaf blight (NCLB) is one of the most devastating leaf pathogens in maize (Zea mays L.). Maize cultivars need to be equipped with broad and stable NCLB resistance to cope with production intensification and climate change. Brazilian germplasm is a great source to increase low NCLB resistance levels in European materials, but little is known about their effect in European environments. To investigate the usefulness of Brazilian germplasm as NCLB resistance donors, we conducted multi-parent QTL mapping, evaluated the potential of marker-assisted selection as well as genome-wide selection of 742 F1-derived DH lines. The line per se performance was evaluated in one location in Brazil and six location-by-year combinations (= environments) in Europe, while testcrosses were assessed in two locations in Brazil and further 10 environments in Europe. Jointly, we identified 17 QTL for NCLB resistance explaining 3.57-30.98% of the genotypic variance each. Two of these QTL were detected in both Brazilian and European environments indicating the stability of these QTL in contrasting ecosystems. We observed moderate to high genomic prediction accuracies between 0.58 and 0.83 depending on population and continent. Collectively, our study illustrates the potential use of tropical resistance sources to increase NCLB resistance level in applied European maize breeding programs.

Genome-Wide Association Mapping and Genomic Prediction of Anther Extrusion in CIMMYT Hybrid Wheat Breeding Program via Modeling Pedigree, Genomic Relationship, and Interaction With the Environment

Anther extrusion (AE) is the most important male floral trait for hybrid wheat seed production. AE is a complex quantitative trait that is difficult to phenotype reliably in field experiments not only due to high genotype-by-environment effects but also due to the short expression window in the field condition. In this study, we conducted a genome-wide association scan (GWAS) and explored the possibility of applying genomic prediction (GP) for AE in the CIMMYT hybrid wheat breeding program. An elite set of male lines (n = 603) were phenotype for anther count (AC) and anther visual score (VS) across three field experiments in 2017–2019 and genotyped with the 20K Infinitum is elect SNP array. GWAS produced five marker trait associations with small effects. For GP, the main effects of lines (L), environment (E), genomic (G) and pedigree relationships (A), and their interaction effects with environments were used to develop seven statistical models of incremental complexity. The base model used only L and E, whereas the most complex model included L, E, G, A, and G × E and A × E. These models were evaluated in three cross-validation scenarios (CV0, CV1, and CV2). In cross-validation CV0, data from two environments were used to predict an untested environment; in random cross-validation CV1, the test set was never evaluated in any environment; and in CV2, the genotypes in the test set were evaluated in only a subset of environments. The prediction accuracies ranged from −0.03 to 0.74 for AC and −0.01 to 0.54 for VS across different models and CV schemes. For both traits, the highest prediction accuracies with low variance were observed in CV2, and inclusion of the interaction effects increased prediction accuracy for AC only. In CV0, the prediction accuracy was 0.73 and 0.45 for AC and VS, respectively, indicating the high reliability of across environment prediction. Genomic prediction appears to be a very reliable tool for AE in hybrid wheat breeding. Moreover, high prediction accuracy in CV0 demonstrates the possibility of implementing genomic selection across breeding cycles in related germplasm, aiding the rapid breeding cycle.

Hybrid breeding in wheat: how shaping floral biology can offer new perspectives

Hybrid breeding in wheat (Triticum aestivum L.) has the potential to deliver major yield increases. This is a requisite to guarantee food security for increasing population demands and to counterbalance the effects of extreme environmental conditions. Successful hybrid breeding in wheat relies on forced outcrossing while preventing self-pollination. To achieve this, research has been directed towards identifying and improving fertility control systems. To maximise cross-pollination and seed set, however, fertility control systems need to be complemented by breeding phenotypically distinct male and female lines. This review summarises existing and novel male sterility systems for wheat hybridisation. We also consider the genetic resources that can be used to alter wheat's floral development and spike morphology, with a focus on the genetic variation already available. Exploiting these resources can lead to enhanced outcrossing, a key requirement in the progress towards hybrid wheat breeding.

Pollen Grain Counting Using a Cell Counter

The number of pollen grains is a critical part of the reproductive strategies in plants and varies greatly between and within species. In agriculture, pollen viability is important for crop breeding. It is a laborious work to count pollen tubes using a counting chamber under a microscope. Here, we present a method of counting the number of pollen grains using a cell counter. In this method, the counting step is shortened to 3 min per flower, which, in our setting, is more than five times faster than the counting chamber method. This technique is applicable to species with a lower and higher number of pollen grains, as it can count particles in a wide range, from 0 to 20,000 particles, in one measurement. The cell counter also estimates the size of the particles together with the number. Because aborted pollen shows abnormal membrane characteristics and/or a distorted or smaller shape, a cell counter can quantify the number of normal and aborted pollen separately. We explain how to count the number of pollen grains and measure pollen size in Arabidopsis thaliana, Arabidopsis kamchatica, and wheat (Triticum aestivum).

Adaptive reduction of male gamete number in the selfing plant Arabidopsis thaliana

The number of male gametes is critical for reproductive success and varies between and within species. The evolutionary reduction of the number of pollen grains encompassing the male gametes is widespread in selfing plants. Here, we employ genome-wide association study (GWAS) to identify underlying loci and to assess the molecular signatures of selection on pollen number-associated loci in the predominantly selfing plant Arabidopsis thaliana. Regions of strong association with pollen number are enriched for signatures of selection, indicating polygenic selection. We isolate the gene REDUCED POLLEN NUMBER1 (RDP1) at the locus with the strongest association. We validate its effect using a quantitative complementation test with CRISPR/Cas9-generated null mutants in nonstandard wild accessions. In contrast to pleiotropic null mutants, only pollen numbers are significantly affected by natural allelic variants. These data support theoretical predictions that reduced investment in male gametes is advantageous in predominantly selfing species.

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.

Negative dominance and dominance-by-dominance epistatic effects reduce grain-yield heterosis in wide crosses in wheat

The genetics underlying heterosis, the difference in performance of crosses compared with midparents, is hypothesized to vary with relatedness between parents. We established a unique germplasm comprising three hybrid wheat sets differing in the degree of divergence between parents and devised a genetic distance measure giving weight to heterotic loci. Heterosis increased steadily with heterotic genetic distance for all 1903 hybrids. Midparent heterosis, however, was significantly lower in the hybrids including crosses between elite and exotic lines than in crosses among elite lines. The analysis of the genetic architecture of heterosis revealed this to be caused by a higher portion of negative dominance and dominance-by-dominance epistatic effects. Collectively, these results expand our understanding of heterosis in crops, an important pillar toward global food security.

The Wheat 660K SNP array demonstrates great potential for marker-assisted selection in polyploid wheat

The rapid development and application of molecular marker assays have facilitated genomic selection and genome-wide linkage and association studies in wheat breeding. Although PCR-based markers (e.g. simple sequence repeats and functional markers) and genotyping by sequencing have contributed greatly to gene discovery and marker-assisted selection, the release of a more accurate and complete bread wheat reference genome has resulted in the design of single-nucleotide polymorphism (SNP) arrays based on different densities or application targets. Here, we evaluated seven types of wheat SNP arrays in terms of their SNP number, distribution, density, associated genes, heterozygosity and application. The results suggested that the Wheat 660K SNP array contained the highest percentage (99.05%) of genome-specific SNPs with reliable physical positions. SNP density analysis indicated that the SNPs were almost evenly distributed across the whole genome. In addition, 229 266 SNPs in the Wheat 660K SNP array were located in 66 834 annotated gene or promoter intervals. The annotated genes revealed by the Wheat 660K SNP array almost covered all genes revealed by the Wheat 35K (97.44%), 55K (99.73%), 90K (86.9%) and 820K (85.3%) SNP arrays. Therefore, the Wheat 660K SNP array could act as a substitute for other 6 arrays and shows promise for a wide range of possible applications. In summary, the Wheat 660K SNP array is reliable and cost-effective and may be the best choice for targeted genotyping and marker-assisted selection in wheat genetic improvement.

Genetic mapping and genomic selection for maize stalk strength

BACKGROUND

Maize is one of the most important staple crops and is widely grown throughout the world. Stalk lodging can cause enormous yield losses in maize production. However, rind penetrometer resistance (RPR), which is recognized as a reliable measurement to evaluate stalk strength, has been shown to be efficient and useful for improving stalk lodging-resistance. Linkage mapping is an acknowledged approach for exploring the genetic architecture of target traits. In addition, genomic selection (GS) using whole genome markers enhances selection efficiency for genetically complex traits. In the present study, two recombinant inbred line (RIL) populations were utilized to dissect the genetic basis of RPR, which was evaluated in seven growth stages.

RESULTS

The optimal stages to measure stalk strength are the silking phase and stages after silking. A total of 66 and 45 quantitative trait loci (QTL) were identified in each RIL population. Several potential candidate genes were predicted according to the maize gene annotation database and were closely associated with the biosynthesis of cell wall components. Moreover, analysis of gene ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway further indicated that genes related to cell wall formation were involved in the determination of RPR. In addition, a multivariate model of genomic selection efficiently improved the prediction accuracy relative to a univariate model and a model considering RPR-relevant loci as fixed effects.

CONCLUSIONS

The genetic architecture of RPR is highly genetically complex. Multiple minor effect QTL are jointly involved in controlling phenotypic variation in RPR. Several pleiotropic QTL identified in multiple stages may contain reliable genes and can be used to develop functional markers for improving the selection efficiency of stalk strength. The application of genomic selection to RPR may be a promising approach to accelerate breeding process for improving stalk strength and enhancing lodging-resistance.

Refining the genetic architecture of flag leaf glaucousness in wheat

The cuticle is the plant's barrier against abiotic and biotic stresses, and the deposition of epicuticular wax crystals results in the scattering of light, an effect termed glaucousness. Here, we dissect the genetic architecture of flag leaf glaucousness in wheat toward a future targeted design of the cuticle. The cuticle serves as a barrier that protects plants against abiotic and biotic stresses. Differences in cuticle composition can be detected by the scattering of light on epicuticular wax crystals, which causes a phenotype termed glaucousness. In this study, we dissected the genetic architecture of flag leaf glaucousness in a panel of 1106 wheat cultivars of global origin. We observed a large genotypic variation, but the geographic pattern suggests that other wax layer characteristics besides glaucousness may be important in conferring tolerance to abiotic stresses such as heat and drought. Genome-wide association mapping identified two major quantitative trait loci (QTL) on chromosomes 3A and 2B. The latter corresponds to the W1 locus, but further characterization revealed that it is likely to contain additional QTL. The same holds true for the major QTL on 3A, which was also found to show an epistatic interaction with another locus located a few centiMorgan distal to it. Genome-wide prediction and the identification of a few additional putative QTL revealed that small-effect QTL also contribute to the trait. Collectively, our results illustrate the complexity of the genetic control of flag leaf glaucousness, with additive effects and epistasis, and lay the foundation for the cloning of the underlying genes toward a more targeted design of the cuticle by plant breeding.

Genome-Wide Association Mapping and Genomic Prediction of Anther Extrusion in CIMMYT Hybrid Wheat Breeding Program via Modeling Pedigree, Genomic Relationship, and Interaction With the Environment

Anther extrusion (AE) is the most important male floral trait for hybrid wheat seed production. AE is a complex quantitative trait that is difficult to phenotype reliably in field experiments not only due to high genotype-by-environment effects but also due to the short expression window in the field condition. In this study, we conducted a genome-wide association scan (GWAS) and explored the possibility of applying genomic prediction (GP) for AE in the CIMMYT hybrid wheat breeding program. An elite set of male lines (n = 603) were phenotype for anther count (AC) and anther visual score (VS) across three field experiments in 2017-2019 and genotyped with the 20K Infinitum is elect SNP array. GWAS produced five marker trait associations with small effects. For GP, the main effects of lines (L), environment (E), genomic (G) and pedigree relationships (A), and their interaction effects with environments were used to develop seven statistical models of incremental complexity. The base model used only L and E, whereas the most complex model included L, E, G, A, and G × E and A × E. These models were evaluated in three cross-validation scenarios (CV0, CV1, and CV2). In cross-validation CV0, data from two environments were used to predict an untested environment; in random cross-validation CV1, the test set was never evaluated in any environment; and in CV2, the genotypes in the test set were evaluated in only a subset of environments. The prediction accuracies ranged from -0.03 to 0.74 for AC and -0.01 to 0.54 for VS across different models and CV schemes. For both traits, the highest prediction accuracies with low variance were observed in CV2, and inclusion of the interaction effects increased prediction accuracy for AC only. In CV0, the prediction accuracy was 0.73 and 0.45 for AC and VS, respectively, indicating the high reliability of across environment prediction. Genomic prediction appears to be a very reliable tool for AE in hybrid wheat breeding. Moreover, high prediction accuracy in CV0 demonstrates the possibility of implementing genomic selection across breeding cycles in related germplasm, aiding the rapid breeding cycle.

Identification of consistent QTL with large effect on anther extrusion in doubled haploid populations developed from spring wheat accessions in German Federal ex situ Genebank

Novel large-effect consistent QTL for anther extrusion (AE) to improve cross-pollination were mapped in doubled haploid populations derived from IPK gene bank spring wheat accessions. TaAP2-D, an ortholog of Cleistogamy1 in barley, is a likely candidate gene for AE in wheat. To establish a robust hybrid wheat breeding system, male lines harboring alleles that promote outcrossing should be developed. In this study, we developed two doubled haploid (DH) populations of hexaploid spring wheat (Triticum aestivum L.) by crossing accessions taken from IPK gene bank. In both populations, the phenotypic data of anther extrusion (AE) based on three years of field trials showed a wide variation and approximated a normal distribution. Both populations were genotyped with a 15 k Infinium single nucleotide polymorphism (SNP) array resulting in 3567 and 3457 polymorphic SNP markers for DH population-1 and DH population-2, respectively. Composite interval mapping identified quantitative trait loci (QTL) on chromosomes 1D, 2D, 4A, 4B, 5A, 5D, 6A, and 6B; with consistent QTL (that are identified in all the years) on chromosome 4A in DH population-1, and on chromosomes 2D and 6B in DH population-2. The consistent QTL explained 17.2%, 32.9%, and 12.3% of the phenotypic variances, respectively. Genic scan of the chromosome 2D-QTL showed that the wheat gene TaAP2-D, an ortholog of Cleistogamy1 which promotes AE via swelling of the lodicules in barley, lies within the QTL region. A diagnostic marker was developed for TaAP2-D that showed co-segregation with the AE phenotype. This study shows the use of gene bank diversity reservoir to find alleles which are otherwise difficult to detect in elite populations. The identification of large-effect consistent QTL for AE is expected to help form efficient male parental lines suitable for hybrid wheat seed production and serve as a source for map-based cloning.

Identification of a candidate gene for a QTL for spikelet number per spike on wheat chromosome arm 7AL by high-resolution genetic mapping

A high-resolution genetic map combined with haplotype analyses identified a wheat ortholog of rice gene APO1 as the best candidate gene for a 7AL locus affecting spikelet number per spike. A better understanding of the genes controlling differences in wheat grain yield components can accelerate the improvements required to satisfy future food demands. In this study, we identified a promising candidate gene underlying a quantitative trait locus (QTL) on wheat chromosome arm 7AL regulating spikelet number per spike (SNS). We used large heterogeneous inbred families ( > 10,000 plants) from two crosses to map the 7AL QTL to an 87-kb region (674,019,191-674,106,327 bp, RefSeq v1.0) containing two complete and two partial genes. In this region, we found three major haplotypes that were designated as H1, H2 and H3. The H2 haplotype contributed the high-SNS allele in both H1 × H2 and H2 × H3 segregating populations. The ancestral H3 haplotype is frequent in wild emmer (48%) but rare (~ 1%) in cultivated wheats. By contrast, the H1 and H2 haplotypes became predominant in modern cultivated durum and common wheat, respectively. Among the four candidate genes, only TraesCS7A02G481600 showed a non-synonymous polymorphism that differentiated H2 from the other two haplotypes. This gene, designated here as WHEAT ORTHOLOG OF APO1 (WAPO1), is an ortholog of the rice gene ABERRANT PANICLE ORGANIZATION 1 (APO1), which affects spikelet number. Taken together, the high-resolution genetic map, the association between polymorphisms in the different mapping populations with differences in SNS, and the known role of orthologous genes in other grass species suggest that WAPO-A1 is the most likely candidate gene for the 7AL SNS QTL among the four genes identified in the candidate gene region.

Hybrid wheat: past, present and future

The review outlines past failures, present status and future prospects of hybrid wheat, and includes information on CMS/CHA/transgenic approaches for male sterility, heterotic groups and cost-effective hybrid seed production. Hybrid varieties give increased yield and improved grain quality in both cross- and self-pollinated crops. However, hybrid varieties in self-pollinated crops (particularly cereals) have not been very successful, except for hybrid rice in China. In case of hybrid wheat, despite the earlier failures, renewed efforts in recent years have been made and hybrid varieties with desirable attributes have been produced and marketed in some European countries. This review builds upon previous reviews, with a new outlook and improved knowledge base, not covered in earlier reviews. New technologies have been described, which include the Hordeum chilense-based CMS-fertility restorer system, chromosomal XYZ-4E-ms system and the following transgenic technologies: (1) conditional male sterility involving use of tapetum-specific expression of a gene that converts a pro-toxin into a phytotoxin causing male sterility; (2) barnase-barstar SeedLink system of Bayer CropScience; (3) split-barnase system that obviates the need of a barstar-based male restorer line; and (4) seed production technology of DuPont-Pioneer that makes use of transgenes in production of male-sterile lines, but gives hybrid seed with no transgenes. This review also includes a brief account of studies for discovery of heterotic QTL, genomic prediction of hybrid vigour and the development of heterotic groups/patterns and their importance in hybrid wheat production. The problem of high cost of hybrid seed due to required high seed rate in wheat relative to hybrid rice has also been addressed. The review concludes with a brief account of the current efforts and future possibilities in making hybrid wheat a commercial success.

High-resolution mapping of rachis nodes per rachis, a critical determinant of grain yield components in wheat

Exploring large genomic data sets based on the latest reference genome assembly identifies the rice ortholog APO1 as a key candidate gene for number of rachis nodes per spike in wheat. Increasing grain yield in wheat is a key breeding objective worldwide. Several component traits contribute to grain yield with spike attributes being among the most important. In this study, we performed a genome-wide association analysis for 12 grain yield and component traits measured in field trials with contrasting agrochemical input levels in a panel of 220 hexaploid winter wheats. A highly significant, environmentally consistent QTL was detected for number of rachis nodes per rachis (NRN) on chromosome 7AL. The five most significant SNPs formed a strong linkage disequilibrium (LD) block and tagged a 2.23 Mb region. Using pairwise LD for exome SNPs located across this interval in a large worldwide hexaploid wheat collection, we reduced the genomic region for NRN to a 258 Kb interval containing four of the original SNP and six high-confidence genes. The ortholog of one (TraesCS7A01G481600) of these genes in rice was ABBERANT PANICLE ORGANIZATION1 (APO1), which is known to have significant effects on panicle attributes. The APO1 ortholog was the best candidate for NRN and was associated with a 115 bp promoter deletion and two amino acid (C47F and D384 N) changes. Using a large worldwide collection of tetraploid and hexaploid wheat, we found 12 haplotypes for the NRN QTL and evidence for positive enrichment of two haplotypes in modern germplasm. Comparison of five QTL haplotypes in Australian yield trials revealed their relative, context-dependent contribution to grain yield. Our study provides diagnostic SNPs and value propositions to support deployment of the NRN trait in wheat breeding.

Fine-mapping of the Fusarium head blight resistance QTL Qfhs.ifa-5A identifies two resistance QTL associated with anther extrusion

Fine-mapping separated Qfhs.ifa-5A into a major QTL mapping across the centromere and a minor effect QTL positioned at the distal half of 5AS. Both increase Fusarium resistance and anther extrusion. The Fusarium head blight (FHB) resistance QTL Qfhs.ifa-5A resides in the low-recombinogenic pericentromeric region of chromosome 5A making fine-mapping particularly arduous. Qfhs.ifa-5A primarily contributes resistance to fungal entry with the favorable allele descending from the highly Fusarium resistant cultivar Sumai-3. Fine-mapping a near-isogenic recombinant inbred line population partitioned the Qfhs.ifa-5A interval into 12 bins. Near-isogenic lines recombining at the interval were phenotyped for FHB severity, anther retention and plant height. Composite interval mapping separated the initially single QTL into two QTL. The major effect QTL Qfhs.ifa-5Ac mapped across the centromere and the smaller effect QTL Qfhs.ifa-5AS mapped to the distal half of 5AS. Although Qfhs.ifa-5Ac and Qfhs.ifa-5AS intervals were as small as 0.1 and 0.2 cM, their corresponding physical distances were large, comprising 44.1 Mbp and 49.2 Mbp, respectively. Sumai-3 alleles at either QTL improved FHB resistance and increased anther extrusion suggesting a pleiotropic effect of anthers on resistance. This hypothesis was supported by greenhouse experiments using the susceptible cultivar Remus and its resistant near-isogenic line NIL3 carrying the entire Qfhs.ifa-5A segment. By manually removing anthers prior to spray inoculation both, Remus and NIL3 became almost equally resistant in the early phase of the disease development and were significantly less diseased than variants without anther manipulation. At late time points the positive effect of the anther removal became smaller for Remus and disappeared completely for NIL3. Results affirm that absence of anthers enhanced resistance to initial infection but did not protect plants from fungal spreading within spikes.

Effects of Rht-B1 and Ppd-D1 loci on pollinator traits in wheat

Elite wheat pollinators are critical for successful hybrid breeding. We identified Rht-B1 and Ppd-D1 loci affecting multiple pollinator traits and therefore represent major targets for improving hybrid seed production. Hybrid breeding has a great potential to significantly boost wheat yields. Ideal male pollinators would be taller in stature, contain many spikelets well-spaced along the spike and exhibit high extrusion of large anthers. Most importantly, flowering time would match with that of the female parent. Available genetic resources for developing an elite wheat pollinator are limited, and the genetic basis for many of these traits is largely unknown. Here, we report on the genetic analysis of pollinator traits using biparental mapping populations. We identified two anther extrusion QTLs of medium effect, one on chromosome 1BL and the other on 4BS coinciding with the semi-dwarfing Rht-B1 locus. The effect of Rht-B1 alleles on anther extrusion is genotype dependent, while tall plant Rht-B1a allele is consistently associated with large anthers. Multiple QTLs were identified at the Ppd-D1 locus for anther length, spikelet number and spike length, with the photoperiod-sensitive Ppd-D1b allele associated with favourable pollinator traits in the populations studied. We also demonstrated that homeoloci, Rht-D1 and Ppd-B1, influence anther length among other traits. These results suggest that combinations of Rht-B1 and Ppd-D1 alleles control multiple pollinator traits and should be major targets of hybrid wheat breeding programs.