Genome-wide Association Mapping and Prediction of Adult Stage Septoria tritici Blotch Infection in European Winter Wheat via High-Density Marker Arrays

Haplotype-based association mapping of genomic regions associated with Zymoseptoria tritici resistance using 217 diverse wheat genotypes

BACKGROUND

Septoria tritici blotch (STB) is considered to be one of the most destructive foliar wheat diseases and is caused by Zymoseptoria tritici. The yield losses are severe and in Northwestern Europe can reach up to 50%. The efficacy of fungicides is diminishing due to changes in the genetic structure of the pathogen. Therefore, resistance breeding is the most effective strategy of disease management. Recently, genome-wide association studies (GWAS) have become more popular due to their robustness in dissecting complex traits, including STB resistance in wheat. This was made possible by the use of large mapping populations and new sequencing technologies. High-resolution mapping benefits from historical recombination and greater allele numbers in GWAS.

RESULTS

In our study, 217 wheat genotypes of diverse origin were phenotyped against five Z. tritici isolates (IPO323, IPO88004, IPO92004, IPO86036 and St1-03) and genotyped on the DArTseq platform. In polytunnel tests two disease parameters were evaluated: the percentage of leaf area covered by necrotic lesions (NEC) and the percentage of leaf area covered by lesions bearing pycnidia (PYC). The disease escape parameters heading date (Hd) and plant height (Ht) were also measured. Pearson's correlation showed a positive effect between disease parameters, providing additional information. The Structure analysis indicated four subpopulations which included from 28 (subpopulation 2) to 79 genotypes (subpopulation 3). All of the subpopulations showed a relatively high degree of admixture, which ranged from 60% of genotypes with less than 80% of proportions of the genome attributed to assigned subpopulation for group 2 to 85% for group 4. Haplotype-based GWAS analysis allowed us to identify 27 haploblocks (HBs) significantly associated with analysed traits with a p-value above the genome-wide significance threshold (5%, which was -log10(p) > 3.64) and spread across the wheat genome. The explained phenotypic variation of identified significant HBs ranged from 0.2% to 21.5%. The results of the analysis showed that four haplotypes (HTs) associated with disease parameters cause a reduction in the level of leaf coverage by necrosis and pycnidia, namely: Chr3A_HB98_HT2, Chr5B_HB47_HT1, Chr7B_HB36_HT1 and Chr5D_HB10_HT3.

CONCLUSIONS

GWAS analysis enabled us to identify four significant chromosomal regions associated with a reduction in STB disease parameters. The list of valuable HBs and wheat varieties possessing them provides promising material for further molecular analysis of resistance loci and development of breeding programmes.

Genomic prediction and genome-wide association study using combined genotypic data from different genotyping systems: application to apple fruit quality traits

With advances in next-generation sequencing technologies, various marker genotyping systems have been developed for genomics-based approaches such as genomic selection (GS) and genome-wide association study (GWAS). As new genotyping platforms are developed, data from different genotyping platforms must be combined. However, the potential use of combined data for GS and GWAS has not yet been clarified. In this study, the accuracy of genomic prediction (GP) and the detection power of GWAS increased for most fruit quality traits of apples when using combined data from different genotyping systems, Illumina Infinium single-nucleotide polymorphism array and genotyping by random amplicon sequencing-direct (GRAS-Di) systems. In addition, the GP model, which considered the inbreeding effect, further improved the accuracy of the seven fruit traits. Runs of homozygosity (ROH) islands overlapped with the significantly associated regions detected by the GWAS for several fruit traits. Breeders may have exploited these regions to select promising apples by breeders, increasing homozygosity. These results suggest that combining genotypic data from different genotyping platforms benefits the GS and GWAS of fruit quality traits in apples. Information on inbreeding could be beneficial for improving the accuracy of GS for fruit traits of apples; however, further analysis is required to elucidate the relationship between the fruit traits and inbreeding depression (e.g. decreased vigor).

Association mapping of septoria tritici blotch resistance in bread wheat in Bale and Arsi highlands, Ethiopia

Septoria tritici blotch (STB) caused by the fungal pathogen Zymoseptoria tritici, anamorph Septoria tritici Rob. ex Desm., is an important wheat pathogen worldwide, reported to be major wheat production threating factor, posing considerable yield loss every year. Developing resistant cultivars is an efficient, economical, environmentally friendly and simple approach for managing STB. This study was carried out to evaluate spring bread wheat lines for their reaction to STB disease under field conditions; to associate phenotypic and genotypic data for identification of STB disease resistance; and to identify genomic region(s) associated with resistance to STB in spring bread wheat lines. Two hundred forty (240) spring bread wheat lines were evaluated under field conditions in non-replicated trials, using an augmented design. The trials were conducted at three locations (Kulumsa Agricultural Research Center, Madda Walabu University Research Site and Sinana Agricultural Research Center) in 2017 main cropping season (July to December). Out of these 240 wheat lines, 123 of them were genotyped with 10263 single nucleotide polymorphism (SNPs) markers and population structure and association mapping analysis was done. The wheat lines showed significant variations in percentage disease severity and area under the disease progress curve at all the three locations they were evaluated. The wheat lines were classified as resistant, moderately resistant, moderately susceptible and susceptible based on the percentage disease severity scored. Five wheat lines were found to be resistant to STB in all the three locations and are recommended for direct release by the national program and parentage purposes in wheat breeding programs. The 123 wheat lines were clustered into 3 subpopulations in which the first cluster contained 99 wheat lines; the second 17 and the last one 7. Among the polymorphic 8127 SNPs markers, 26 markers on chromosomes 7B, 1D, 3A, 2B, 6B and 3D were found to be significantly (P < 0.001) associated with STB resistance so that they can be utilized for marker assisted selection and gene pyramiding in resistance breeding programs.

Genome-Wide Association Study Reveals Novel Powdery Mildew Resistance Loci in Bread Wheat

Powdery mildew (PM), caused by the fungal pathogen Blumeria graminis f. sp. tritici (Bgt), significantly threatens global bread wheat production. Although the use of resistant cultivars is an effective strategy for managing PM, currently available wheat cultivars lack sufficient levels of resistance. To tackle this challenge, we conducted a comprehensive genome-wide association study (GWAS) using a diverse panel of 286 bread wheat genotypes. Over three consecutive years (2020-2021, 2021-2022, and 2022-2023), these genotypes were extensively evaluated for PM severity under field conditions following inoculation with virulent Bgt isolates. The panel was previously genotyped using the Illumina 90K Infinium iSelect assay to obtain genome-wide single-nucleotide polymorphism (SNP) marker coverage. By applying FarmCPU, a multilocus mixed model, we identified a total of 113 marker-trait associations (MTAs) located on chromosomes 1A, 1B, 2B, 3A, 3B, 4A, 4B, 5A, 5B, 6B, 7A, and 7B at a significance level of p ≤ 0.001. Notably, four novel MTAs on chromosome 6B were consistently detected in 2020-2021 and 2021-2022. Furthermore, within the confidence intervals of the identified SNPs, we identified 96 candidate genes belonging to different proteins including 12 disease resistance/host-pathogen interaction-related protein families. Among these, protein kinases, leucine-rich repeats, and zinc finger proteins were of particular interest due to their potential roles in PM resistance. These identified loci can serve as targets for breeding programs aimed at developing disease-resistant wheat cultivars.

Genetic control and prospects of predictive breeding for European winter wheat's Zeleny sedimentation values and Hagberg-Perten falling number

KEY MESSAGE

Sedimentation values and falling number in the last decades have helped maintain high baking quality despite rigorous selection for grain yield in wheat. Allelic combinations of major loci sustained the bread-making quality while improving grain yield. Glu-D1, Pinb-D1, and non-gluten proteins are associated with sedimentation values and falling number in European wheat. Zeleny sedimentation values (ZSV) and Hagberg-Perten falling number (HFN) are among the most important parameters that help determine the baking quality classes of wheat and, thus, influence the monetary benefits for growers. We used a published data set of 372 European wheat varieties evaluated in replicated field trials in multiple environments. ZSV and HFN traits hold a wide and significant genotypic variation and high broad-sense heritability. The genetic correlations revealed positive and significant associations of ZSV and HFN with each other, grain protein content (GPC) and grain hardness; however, they were all significantly negatively correlated with grain yield. Besides, GPC appeared to be the major predictor for ZSV and HFN. Our genome-wide association analyses based on high-quality SSR, SNP, and candidate gene markers revealed a strong quantitative genetic nature of ZSV and HFN by explaining their total genotypic variance as 41.49% and 38.06%, respectively. The association of known Glutenin (Glu-1) and Puroindoline (Pin-1) with ZSV provided positive analytic proof of our studies. We report novel candidate loci associated with globulins and albumins-the non-gluten monomeric proteins in wheat. In addition, predictive breeding analyses for ZSV and HFN suggest using genomic selection in the early stages of breeding programs with an average prediction accuracy of 81 and 59%, respectively.

Exploring GWAS and genomic prediction to improve Septoria tritici blotch resistance in wheat

Septoria tritici blotch (STB) is a destructive foliar diseases threatening wheat grain yield. Wheat breeding for STB disease resistance has been identified as the most sustainable and environment-friendly approach. In this work, a panel of 316 winter wheat breeding lines from a commercial breeding program were evaluated for STB resistance at the seedling stage under controlled conditions followed by genome-wide association study (GWAS) and genomic prediction (GP). The study revealed a significant genotypic variation for STB seedling resistance, while disease severity scores exhibited a normal frequency distribution. Moreover, we calculated a broad-sense heritability of 0.62 for the trait. Nine single- and multi-locus GWAS models identified 24 marker-trait associations grouped into 20 quantitative trait loci (QTLs) for STB seedling-stage resistance. The seven QTLs located on chromosomes 1B, 2A, 2B, 5B (two), 7A, and 7D are reported for the first time and could potentially be novel. The GP cross-validation analysis in the RR-BLUP model estimated the genomic-estimated breeding values (GEBVs) of STB resistance with a prediction accuracy of 0.49. Meanwhile, the GWAS assisted wRR-BLUP model improved the accuracy to 0.58. The identified QTLs can be used for marker-assisted backcrossing against STB in winter wheat. Moreover, the higher prediction accuracy recorded from the GWAS-assisted GP analysis implies its power to successfully select superior candidate lines based on their GEBVs for STB resistance.

Biology, taxonomy, genetics, and management of Zymoseptoria tritici: the causal agent of wheat leaf blotch

Septoria tritici blotch or Septoria leaf blotch has been used for long time, but leaf blotch is a correct disease name. Moreover, Lb resistant gene is the correct name, but, not Stb gene. It has sexual and asexual parts on the mycelia, known as heterothallic fungi. Its pathogenic diversity ranged from 40% to 93% and has produced a wide variety of AvrLb6 haplotypes. M. graminicola has a plasmogamy and karyogamy sexual process. The pathogen can use macropycnidiospores, micropycnidiospores, and pycnidia vegetative growths for infection and overwintering. Synthetic M3, Kavkaz-K4500, Synthetic 6×, and TE9111 wheat genotypes have horizontal resistance. Avirulence (Avr) genes in Z. tritici and their matching wheat (R) genes indicate gene for gene mechanisms of resistance. Twenty-two R genes (vertical resistance) have been identified. In both horizontal and vertical resistance, different Lb genes have been broken down due to new Z.tritici virulent gene and currently Lb19 resistant gene is being recommended. Mixing of resistant and susceptible cultivars is also the most effective management strategy. Moreover, different cultural practices and biological control have been proposed. Lastly, different fungicides are also available. However, in developing countries cultivar mixture, isolates diversity, biological control, and epidemic studies have been greatly missed.

Multi-omics assisted breeding for biotic stress resistance in soybean

Biotic stress is a critical factor limiting soybean growth and development. Soybean responses to biotic stresses such as insects, nematodes, fungal, bacterial, and viral pathogens are governed by complex regulatory and defense mechanisms. Next-generation sequencing has availed research techniques and strategies in genomics and post-genomics. This review summarizes the available information on marker resources, quantitative trait loci, and marker-trait associations involved in regulating biotic stress responses in soybean. We discuss the differential expression of related genes and proteins reported in different transcriptomics and proteomics studies and the role of signaling pathways and metabolites reported in metabolomic studies. Recent advances in omics technologies offer opportunities to reshape and improve biotic stress resistance in soybean by altering gene regulation and/or other regulatory networks. We suggest using 'integrated omics' to precisely understand how soybean responds to different biotic stresses. We also discuss the potential challenges of integrating multi-omics for the functional analysis of genes and their regulatory networks and the development of biotic stress-resistant cultivars. This review will help direct soybean breeding programs to develop resistance against different biotic stresses.

Multi-stage resistance to Zymoseptoria tritici revealed by GWAS in an Australian bread wheat diversity panel

Septoria tritici blotch (STB) has been ranked the third most important wheat disease in the world, threatening a large area of wheat production. Although major genes play an important role in the protection against Zymoseptoria tritici infection, the lifespan of their resistance unfortunately is very short in modern wheat production systems. Combinations of quantitative resistance with minor effects, therefore, are believed to have prolonged and more durable resistance to Z. tritici. In this study, new quantitative trait loci (QTLs) were identified that are responsible for seedling-stage resistance and adult-plant stage resistance (APR). More importantly was the characterisation of a previously unidentified QTL that can provide resistance during different stages of plant growth or multi-stage resistance (MSR). At the seedling stage, we discovered a new isolate-specific QTL, QSt.wai.1A.1. At the adult-plant stage, the new QTL QStb.wai.6A.2 provided stable and consistent APR in multiple sites and years, while the QTL QStb.wai.7A.2 was highlighted to have MSR. The stacking of multiple favourable MSR alleles was found to improve resistance to Z. tritici by up to 40%.

QTL mapping and genomic prediction of resistance to wheat head blight caused by Fusarium verticillioides

Fusarium head blight (FHB), is one of the destructive fugue diseases of wheat worldwide caused by the Fusarium verticillioides (F.v). In this study, a population consisting of 262 recombinant inbred lines (RILs) derived from Zhongmai 578 and Jimai 22 was used to map Quantitative Trait Locus (QTL) for FHB resistance, with the genotype data using the wheat 50 K single nucleotide polymorphism (SNP) array. The percentage of symptomatic spikelet (PSS) and the weighted average of PSS (PSSW) were collected for each RIL to represent their resistance to wheat head blight caused by F.v. In total, 22 QTL associated with FHB resistance were identified on chromosomes 1D, 2B, 3B, 4A, 5D, 7A, 7B, and 7D, respectively, from which 10 and 12 QTL were detected from PSS and PSSW respectively, explaining 3.82%-10.57% of the phenotypic variances using the inclusive composite interval mapping method. One novel QTL, Qfhb. haust-4A.1, was identified, explaining 10.56% of the phenotypic variation. One stable QTL, Qfhb. haust-1D.1 was detected on chromosome 1D across multiple environments explaining 4.39%-5.70% of the phenotypic variation. Forty-seven candidate genes related to disease resistance were found in the interval of Qfhb. haust-1D.1 and Qfhb. haust-4A.1. Genomic prediction accuracies were estimated from the five-fold cross-validation scheme ranging from 0.34 to 0.40 for PSS, and from 0.34 to 0.39 for PSSW in in-vivo inoculation treatment. This study provided new insight into the genetic analysis of resistance to wheat head blight caused by F.v, and genomic selection (GS) as a potential approach for improving the resistance of wheat head blight.

Gabi wheat a panel of European elite lines as central stock for wheat genetic research

In plant sciences, curation and availability of interoperable phenotypic and genomic data is still in its infancy and represents an obstacle to rapid scientific discoveries in this field. To that end, supplementing the efforts being made to generate open access wheat genome, pan wheat genome and other bioinformatic resources, we present the GABI-WHEAT panel of elite European cultivars comprising 358 winter and 14 summer wheat varieties released between 1975 to 2007. The panel has been genotyped with SNP arrays of increasing density to investigate several important agronomic, quality and disease resistance traits. The robustness of investigated traits and interoperability of genomic and phenotypic data was assessed in the current publication with the aim to transform this panel into a public data resource for future genetic research in wheat. Consecutively, the phenotypic data was formatted to comply with FAIR principles and linked to online databases to substantiate panel origin information and quality. Thus, we were able to make a valuable resource available for plant science in a sustainable way.

Measurement(s)	agronomic, quality and disease traits
Technology Type(s)	manual measurement in the field
Sample Characteristic - Organism	Triticum aestivum L.

Genetic analysis of novel resistance sources and genome-wide association mapping identified novel QTLs for resistance to Zymoseptoria tritici, the causal agent of septoria tritici blotch in wheat

Septoria tritici blotch (STB) caused by Zymoseptoria tritici is one of the most important foliar diseases of wheat causing significant yield losses worldwide. In this study, a panel of bread wheat genotypes comprised 185 globally diverse genotypes were tested against 10 Z. tritici isolates at the seedling stage. Genome-wide association study (GWAS) using high-throughput DArTseq markers was performed and further gene expression analysis of significant markers trait association (MTAs) associated with resistance to STB was analyzed. Disease severity level showed significant differences among wheat genotypes for resistance to different Z. tritici isolates. We found novel landrace genotypes that showed highly resistance spectra to all tested isolates. GWAS analysis resulted in 19 quantitative trait loci (QTLs) for resistance to STB that were located on 14 chromosomes. Overall, 14 QTLs were overlapped with previously known QTLs or resistance genes, as well as five potentially novel QTLs on chromosomes 1A, 4A, 5B, 5D, and 6D. Identified novel resistance sources and also novel QTLs for resistance to different Z. tritici isolates can be used for gene pyramiding and development of durable resistance cultivars in future wheat breeding programs.

Meta-analysis reveals consensus genomic regions associated with multiple disease resistance in wheat (Triticum aestivum L.)

In wheat, meta-QTLs (MQTLs) and candidate genes (CGs) were identified for multiple disease resistance (MDR). For this purpose, information was collected from 58 studies for mapping QTLs for resistance to one or more of the five diseases. As many as 493 QTLs were available from these studies, which were distributed in five diseases as follows: septoria tritici blotch (STB) 126 QTLs; septoria nodorum blotch (SNB), 103 QTLs; fusarium head blight (FHB), 184 QTLs; karnal bunt (KB), 66 QTLs; and loose smut (LS), 14 QTLs. Of these 493 QTLs, only 291 QTLs could be projected onto a consensus genetic map, giving 63 MQTLs. The CI of the MQTLs ranged from 0.04 to 15.31 cM with an average of 3.09 cM per MQTL. This is a ~ 4.39 fold reduction from the CI of QTLs, which ranged from 0 to 197.6 cM, with a mean of 13.57 cM. Of 63 MQTLs, 60 were anchored to the reference physical map of wheat (the physical interval of these MQTLs ranged from 0.30 to 726.01 Mb with an average of 74.09 Mb). Thirty-eight (38) of these MQTLs were verified using marker-trait associations (MTAs) derived from genome-wide association studies. As many as 874 CGs were also identified which were further investigated for differential expression using data from five transcriptome studies, resulting in 194 differentially expressed candidate genes (DECGs). Among the DECGs, 85 genes had functions previously reported to be associated with disease resistance. These results should prove useful for fine mapping and cloning of MDR genes and marker-assisted breeding.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01282-z.

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.

Phenotypic Evaluation and Genetic Analysis of Seedling Emergence in a Global Collection of Wheat Genotypes (Triticum aestivum L.) Under Limited Water Availability

The challenge in establishing an early-sown wheat crop in southern Australia is the need for consistently high seedling emergence when sowing deep in subsoil moisture (>10 cm) or into dry top-soil (4 cm). However, the latter is strongly reliant on a minimum soil water availability to ensure successful seedling emergence. This study aimed to: (1) evaluate 233 Australian and selected international wheat genotypes for consistently high seedling emergence under limited soil water availability when sown in 4 cm of top-soil in field and glasshouse (GH) studies; (2) ascertain genetic loci associated with phenotypic variation using a genome-wide association study (GWAS); and (3) compare across loci for traits controlling coleoptile characteristics, germination, dormancy, and pre-harvest sprouting. Despite significant (P < 0.001) environment and genotype-by-environment interactions within and between field and GH experiments, eight genotypes that included five cultivars, two landraces, and one inbred line had consistently high seedling emergence (mean value > 85%) across nine environments. Moreover, 21 environment-specific quantitative trait loci (QTL) were detected in GWAS analysis on chromosomes 1B, 1D, 2B, 3A, 3B, 4A, 4B, 5B, 5D, and 7D, indicating complex genetic inheritance controlling seedling emergence. We aligned QTL for known traits and individual genes onto the reference genome of wheat and identified 16 QTL for seedling emergence in linkage disequilibrium with coleoptile length, width, and cross-sectional area, pre-harvest sprouting and dormancy, germination, seed longevity, and anthocyanin development. Therefore, it appears that seedling emergence is controlled by multifaceted networks of interrelated genes and traits regulated by different environmental cues.

Genome-Wide Association Study Reveals Novel Genetic Loci for Quantitative Resistance to Septoria Tritici Blotch in Wheat (Triticum aestivum L.)

Septoria tritici blotch, caused by the fungus Zymoseptoria titici, poses serious and persistent challenges to wheat cultivation in Ethiopia and worldwide. Deploying resistant cultivars is a major component of controlling septoria tritici blotch (STB). Thus, the objective of this study was to elucidate the genomic architecture of STB resistance in an association panel of 178 bread wheat genotypes. The association panel was phenotyped for STB resistance, phenology, yield, and yield-related traits in three locations for 2 years. The panel was also genotyped for single nucleotide polymorphism (SNP) markers using the genotyping-by-sequencing (GBS) method, and a total of 7,776 polymorphic SNPs were used in the subsequent analyses. Marker-trait associations were also computed using a genome association and prediction integrated tool (GAPIT). The study then found that the broad-sense heritability for STB resistance ranged from 0.58 to 0.97 and 0.72 to 0.81 at the individual and across-environment levels, respectively, indicating the presence of STB resistance alleles in the association panel. Population structure and principal component analyses detected two sub-groups with greater degrees of admixture. A linkage disequilibrium (LD) analysis in 338,125 marker pairs also detected the existence of significant (p ≤ 0.01) linkage in 27.6% of the marker pairs. Specifically, in all chromosomes, the LD between SNPs declined within 2.26-105.62 Mbp, with an overall mean of 31.44 Mbp. Furthermore, the association analysis identified 53 loci that were significantly (false discovery rate, FDR, <0.05) associated with STB resistance, further pointing to 33 putative quantitative trait loci (QTLs). Most of these shared similar chromosomes with already published Septoria resistance genes, which were distributed across chromosomes 1B, 1D, 2A, 2B, 2D, 3A,3 B, 3D, 4A, 5A, 5B, 6A, 7A, 7B, and 7D. However, five of the putative QTLs identified on chromosomes 1A, 5D, and 6B appeared to be novel. Dissecting the detected loci on IWGSC RefSeq Annotation v2.1 revealed the existence of disease resistance-associated genes in the identified QTL regions that are involved in plant defense responses. These putative QTLs explained 2.7-13.2% of the total phenotypic variation. Seven of the QTLs (R ² = 2.7-10.8%) for STB resistance also co-localized with marker-trait associations (MTAs) for agronomic traits. Overall, this analysis reported on putative QTLs for adult plant resistance to STB and some important agronomic traits. The reported and novel QTLs have been identified previously, indicating the potential to improve STB resistance by pyramiding QTLs by marker-assisted selection.

Evaluation of Septoria Nodorum Blotch (SNB) Resistance in Glumes of Wheat (Triticum aestivum L.) and the Genetic Relationship With Foliar Disease Response

Septoria nodorum blotch (SNB) is a necrotrophic disease of wheat prominent in some parts of the world, including Western Australia (WA) causing significant losses in grain yield. The genetic mechanisms for resistance are complex involving multiple quantitative trait loci. In order to decipher comparable or independent regulation, this study identified the genetic control for glume compared to foliar resistance across four environments in WA against 37 different isolates. High proportion of the phenotypic variation across environments was contributed by genotype (84.0% for glume response and 82.7% for foliar response) with genotype-by-environment interactions accounting for a proportion of the variation for both glume and foliar response (14.7 and 16.2%, respectively). Despite high phenotypic correlation across environments, most of the eight and 14 QTL detected for glume and foliar resistance using genome wide association analysis (GWAS), respectively, were identified as environment-specific. QTL for glume and foliar resistance neither co-located nor were in LD in any particular environment indicating autonomous genetic mechanisms control SNB response in adult plants, regulated by independent biological mechanisms and influenced by significant genotype-by- environment interactions. Known Snn and Tsn loci and QTL were compared with 22 environment-specific QTL. None of the eight QTL for glume or the 14 for foliar response were co-located or in linkage disequilibrium with Snn and only one foliar QTL was in LD with Tsn loci on the physical map. Therefore, glume and foliar response to SNB in wheat is regulated by multiple environment-specific loci which function independently, with limited influence of known NE-Snn interactions for disease progression in Western Australian environments. Breeding for stable resistance would consequently rely on recurrent phenotypic selection to capture and retain favorable alleles for both glume and foliar resistance relevant to a particular environment.

Identification of Resistance Sources and Genome-Wide Association Mapping of Septoria Tritici Blotch Resistance in Spring Bread Wheat Germplasm of ICARDA

Septoria tritici blotch (STB) of wheat, caused by the ascomycete Zymoseptoria tritici (formerly Mycosphaerella graminicola), is one of the most important foliar diseases of wheat. In Morocco, STB is a devastating disease in temperate wheat-growing regions, and the yield losses can exceed up to 50% under favorable conditions. The aims of this study were to identify sources of resistance to STB in Septoria Association Mapping Panel (SAMP), which is composed of 377 advanced breeding lines (ABLs) from spring bread wheat breeding program of ICARDA, and to identify loci associated with resistance to STB at seedling (SRT) as well as at the adult plant (APS) stages using genome-wide association mapping (GWAM). Seedling resistance was evaluated under controlled conditions with two virulent isolates of STB (SAT-2 and 71-R3) from Morocco, whereas adult plant resistance was assessed at two hot spot locations in Morocco (Sidi Allal Tazi, Marchouch) under artificial inoculation with a mixture of STB isolates. At seedling stage, 45 and 32 ABLs were found to be resistant to 71-R3 and SAT-2 isolates of STB, respectively. At adult plant stage, 50 ABLs were found to be resistant at hot spot locations in Morocco. Furthermore, 10 genotypes showed resistance in both locations during two cropping seasons. GWAM was conducted with 9,988 SNP markers using phenotypic data for seedling and the adult plant stage. MLM model was employed in TASSEL 5 (v 5.2.53) using principal component analysis and Kinship Matrix as covariates. The GWAM analysis indicated 14 quantitative trait loci (QTL) at the seedling stage (8 for isolate SAT-2 and 6 for isolate 71-R3), while 23 QTL were detected at the adult plant stage resistance (4 at MCH-17, 16 at SAT-17, and 3 at SAT-18). SRT QTL explained together 33.3% of the phenotypic variance for seedling resistance to STB isolate SAT-2 and 28.3% for 71-R3, respectively. QTL for adult plant stage resistance explained together 13.1, 68.6, and 11.9% of the phenotypic variance for MCH-17, SAT-17, and SAT-18, respectively. Identification of STB-resistant spring bread wheat germplasm in combination with QTL detected both at SRT and APS stage will serve as an important resource in STB resistance breeding efforts.

Genome-Wide Association Analysis and Genomic Prediction for Adult-Plant Resistance to Septoria Tritici Blotch and Powdery Mildew in Winter Wheat

Septoria tritici blotch (STB) caused by the fungal pathogen Zymoseptoria tritici and powdery mildew (PM) caused by Blumeria graminis f.sp tritici (Bgt) are among the forefront foliar diseases of wheat that lead to a significant loss of grain yield and quality. Resistance breeding aimed at developing varieties with inherent resistance to STB and PM diseases has been the most sustainable and environment-friendly approach. In this study, 175 winter wheat landraces and historical cultivars originated from the Nordic region were evaluated for adult-plant resistance (APR) to STB and PM in Denmark, Estonia, Lithuania, and Sweden. Genome-wide association study (GWAS) and genomic prediction (GP) were performed based on the adult-plant response to STB and PM in field conditions using 7,401 single-nucleotide polymorphism (SNP) markers generated by 20K SNP chip. Genotype-by-environment interaction was significant for both disease scores. GWAS detected stable and environment-specific quantitative trait locis (QTLs) on chromosomes 1A, 1B, 1D, 2B, 3B, 4A, 5A, 6A, and 6B for STB and 2A, 2D, 3A, 4B, 5A, 6B, 7A, and 7B for PM adult-plant disease resistance. GP accuracy was improved when assisted with QTL from GWAS as a fixed effect. The GWAS-assisted GP accuracy ranged within 0.53–0.75 and 0.36–0.83 for STB and PM, respectively, across the tested environments. This study highlights that landraces and historical cultivars are a valuable source of APR to STB and PM. Such germplasm could be used to identify and introgress novel resistance genes to modern breeding lines.

Genome-Wide Association Studies and Prediction of Tan Spot (Pyrenophora tritici-repentis) Infection in European Winter Wheat via Different Marker Platforms

Tan spot, caused by the fungus Pyrenophoratritici-repentis (Ptr), is a severe foliar disease of wheat (Triticumaestivum L.). Improving genetic resistance is a durable strategy to reduce Ptr-related losses. Here, we dissected Ptr-infection’s genetic basis in 372 European wheat varieties via simple sequence repeats (SSRs) plus 35k and 90k single nucleotide polymorphism (SNP) marker platforms. In our phenotypic data analyses, Ptr infection showed a significant genotypic variance and a significant negative correlation with plant height. Genome-wide association studies revealed a highly quantitative nature of Ptr infection and identified two quantitative trait loci (QTL), viz., QTs.ipk-7A and QTs.ipk-7B, which imparted 21.23 and 5.84% of the genotypic variance, respectively. Besides, the Rht-D1 gene showed a strong allelic influence on the infection scores. Due to the complex genetic nature of the Ptr infection, the potential of genome-wide prediction (GP) was assessed via three different genetic models on individual and combined marker platforms. The GP results indicated that the marker density and marker platforms do not considerably impact prediction accuracy (~40–42%) and that higher-order epistatic interactions may not be highly pervasive. Our results provide a further understanding of Ptr-infection’s genetic nature, serve as a resource for marker-assisted breeding, and highlight the potential of genome-wide selection for improved Ptr resistance.

Comparison of Genomic Prediction Methods for Yellow, Stem, and Leaf Rust Resistance in Wheat Landraces from Afghanistan

Wheat rust diseases, including yellow rust (Yr; also known as stripe rust) caused by Puccinia striiformis Westend. f. sp. tritici, leaf rust (Lr) caused by Puccinia triticina Eriks. and stem rust (Sr) caused by Puccinia graminis Pres f. sp. tritici are major threats to wheat production all around the globe. Durable resistance to wheat rust diseases can be achieved through genomic-assisted prediction of resistant accessions to increase genetic gain per unit time. Genomic prediction (GP) is a promising technology that uses genomic markers to estimate genomic-assisted breeding values (GBEVs) for selecting resistant plant genotypes and accumulating favorable alleles for adult plant resistance (APR) to wheat rust diseases. To evaluate GP we compared the predictive ability of nine different parametric, semi-parametric and Bayesian models including Genomic Unbiased Linear Prediction (GBLUP), Ridge Regression (RR), Least Absolute Shrinkage and Selection Operator (LASSO), Elastic Net (EN), Bayesian Ridge Regression (BRR), Bayesian A (BA), Bayesian B (BB), Bayesian C (BC) and Reproducing Kernel Hilbert Spacing model (RKHS) to estimate GEBV’s for APR to yellow, leaf and stem rust of wheat in a panel of 363 bread wheat landraces of Afghanistan origin. Based on five-fold cross validation the mean predictive abilities were 0.33, 0.30, 0.38, and 0.33 for Yr (2016), Yr (2017), Lr, and Sr, respectively. No single model outperformed the rest of the models for all traits. LASSO and EN showed the lowest predictive ability in four of the five traits. GBLUP and RR gave similar predictive abilities, whereas Bayesian models were not significantly different from each other as well. We also investigated the effect of the number of genotypes and the markers used in the analysis on the predictive ability of the GP model. The predictive ability was highest with 1000 markers and there was a linear trend in the predictive ability and the size of the training population. The results of the study are encouraging, confirming the feasibility of GP to be effectively applied in breeding programs for resistance to all three wheat rust diseases.

Genomics-Assisted Breeding for Quantitative Disease Resistances in Small-Grain Cereals and Maize

Generating genomics-driven knowledge opens a way to accelerate the resistance breeding process by family or population mapping and genomic selection. Important prerequisites are large populations that are genomically analyzed by medium- to high-density marker arrays and extensive phenotyping across locations and years of the same populations. The latter is important to train a genomic model that is used to predict genomic estimated breeding values of phenotypically untested genotypes. After reviewing the specific features of quantitative resistances and the basic genomic techniques, the possibilities for genomics-assisted breeding are evaluated for six pathosystems with hemi-biotrophic fungi: Small-grain cereals/Fusarium head blight (FHB), wheat/Septoria tritici blotch (STB) and Septoria nodorum blotch (SNB), maize/Gibberella ear rot (GER) and Fusarium ear rot (FER), maize/Northern corn leaf blight (NCLB). Typically, all quantitative disease resistances are caused by hundreds of QTL scattered across the whole genome, but often available in hotspots as exemplified for NCLB resistance in maize. Because all crops are suffering from many diseases, multi-disease resistance (MDR) is an attractive aim that can be selected by specific MDR QTL. Finally, the integration of genomic data in the breeding process for introgression of genetic resources and for the improvement within elite materials is discussed.

Genome-wide association analysis of stem water-soluble carbohydrate content in bread wheat

GWAS identified 36 potentially new loci for wheat stem water-soluble carbohydrate (WSC) contents and 13 pleiotropic loci affecting WSC and thousand-kernel weight. Five KASP markers were developed and validated. Water-soluble carbohydrates (WSC) reserved in stems contribute significantly to grain yield (GY) in wheat. However, knowledge of the genetic architecture underlying stem WSC content (SWSCC) is limited. In the present study, 166 diverse wheat accessions from the Yellow and Huai Valleys Winter Wheat Zone of China and five other countries were grown in four well-watered environments. SWSCC at 10 days post-anthesis (10DPA), 20DPA and 30DPA, referred as WSC10, WSC20 and WSC30, respectively, and thousand-kernel weight (TKW) were assessed. Correlation analysis showed that TKW was significantly and positively correlated with WSC10 and WSC20. Genome-wide association study was performed on SWSCC and TKW with 373,106 markers from the wheat 660 K and 90 K SNP arrays. Totally, 62 stable loci were detected for SWSCC, with 36, 24 and 19 loci for WSC10, WSC20 and WSC30, respectively; among these, 36 are potentially new, 16 affected SWSCC at two or three time-points, and 13 showed pleiotropic effects on both SWSCC and TKW. Linear regression showed clear cumulative effects of favorable alleles for increasing SWSCC and TKW. Genetic gain analyses indicated that pyramiding favorable alleles of SWSCC had simultaneously improved TKW. Kompetitive allele-specific PCR markers for five pleiotropic loci associated with both SWSCC and TKW were developed and validated. This study provided a genome-wide landscape of the genetic architecture of SWSCC, gave a perspective for understanding the relationship between WSC and GY and explored the theoretical basis for co-improvement of WSC and GY. It also provided valuable loci and markers for future breeding.

Appraising the Genetic Architecture of Kernel Traits in Hexaploid Wheat Using GWAS

Kernel morphology is one of the major yield traits of wheat, the genetic architecture of which is always important in crop breeding. In this study, we performed a genome-wide association study (GWAS) to appraise the genetic architecture of the kernel traits of 319 wheat accessions using 22,905 single nucleotide polymorphism (SNP) markers from a wheat 90K SNP array. As a result, 111 and 104 significant SNPs for Kernel traits were detected using four multi-locus GWAS models (mrMLM, FASTmrMLM, FASTmrEMMA, and pLARmEB) and three single-locus models (FarmCPU, MLM, and MLMM), respectively. Among the 111 SNPs detected by the multi-locus models, 24 SNPs were simultaneously detected across multiple models, including seven for kernel length, six for kernel width, six for kernels per spike, and five for thousand kernel weight. Interestingly, the five most stable SNPs (RAC875_29540_391, Kukri_07961_503, tplb0034e07_1581, BS00074341_51, and BobWhite_049_3064) were simultaneously detected by at least three multi-locus models. Integrating these newly developed multi-locus GWAS models to unravel the genetic architecture of kernel traits, the mrMLM approach detected the maximum number of SNPs. Furthermore, a total of 41 putative candidate genes were predicted to likely be involved in the genetic architecture underlining kernel traits. These findings can facilitate a better understanding of the complex genetic mechanisms of kernel traits and may lead to the genetic improvement of grain yield in wheat.

Genetic Analysis Using a Multi-Parent Wheat Population Identifies Novel Sources of Septoria Tritici Blotch Resistance

Zymoseptoria tritici is the causative fungal pathogen of septoria tritici blotch (STB) disease of wheat (Triticum aestivum L.) that continuously threatens wheat crops in Ireland and throughout Europe. Under favorable conditions, STB can cause up to 50% yield losses if left untreated. STB is commonly controlled with fungicides; however, a combination of Z. tritici populations developing fungicide resistance and increased restrictions on fungicide use in the EU has led to farmers relying on fewer active substances. Consequently, this serves to drive the emergence of Z. tritici resistance against the remaining chemistries. In response, the use of resistant wheat varieties provides a more sustainable disease management strategy. However, the number of varieties offering an adequate level of resistance against STB is limited. Therefore, new sources of resistance or improved stacking of existing resistance loci are needed to develop varieties with superior agronomic performance. Here, we identified quantitative trait loci (QTL) for STB resistance in the eight-founder "NIAB Elite MAGIC" winter wheat population. The population was screened for STB response in the field under natural infection for three seasons from 2016 to 2018. Twenty-five QTL associated with STB resistance were identified in total. QTL either co-located with previously reported QTL or represent new loci underpinning STB resistance. The genomic regions identified and the linked genetic markers serve as useful resources for STB resistance breeding, supporting rapid selection of favorable alleles for the breeding of new wheat cultivars with improved STB resistance.

Prospects of GWAS and predictive breeding for European winter wheat's grain protein content, grain starch content, and grain hardness

Grain quality traits determine the classification of registered wheat (Triticum aestivum L.) varieties. Although environmental factors and crop management practices exert a considerable influence on wheat quality traits, a significant proportion of the variance is attributed to the genetic factors. To identify the underlying genetic factors of wheat quality parameters viz., grain protein content (GPC), grain starch content (GSC), and grain hardness (GH), we evaluated 372 diverse European wheat varieties in replicated field trials in up to eight environments. We observed that all of the investigated traits hold a wide and significant genetic variation, and a significant negative correlation exists between GPC and GSC plus grain yield. Our association analyses based on 26,694 high-quality single nucleotide polymorphic markers revealed a strong quantitative genetic nature of GPC and GSC with associations on groups 2, 3, and 6 chromosomes. The identification of known Puroindoline-b gene for GH provided a positive analytic proof for our studies. We report that a locus QGpc.ipk-6A controls both GPC and GSC with opposite allelic effects. Based on wheat's reference and pan-genome sequences, the physical characterization of two loci viz., QGpc.ipk-2B and QGpc.ipk-6A facilitated the identification of the candidate genes for GPC. Furthermore, by exploiting additive and epistatic interactions of loci, we evaluated the prospects of predictive breeding for the investigated traits that suggested its efficient use in the breeding programs.

Multi-Location Evaluation of Global Wheat Lines Reveal Multiple QTL for Adult Plant Resistance to Septoria Nodorum Blotch (SNB) Detected in Specific Environments and in Response to Different Isolates

The slow rate of genetic gain for improving resistance to Septoria nodorum blotch (SNB) is due to the inherent complex interactions between host, isolates, and environments. Breeding for improved SNB resistance requires evaluation and selection of wheat genotypes consistently expressing low SNB response in different target production environments. The study focused on evaluating 232 genotypes from global origins for resistance to SNB in the flag leaf expressed in different Western Australian environments. The aim was to identify resistant donor germplasm against historical and contemporary pathogen isolates and enhance our knowledge of the genetic basis of genotype-by-environment interactions for SNB response. Australian wheat varieties, inbred lines from Centro Internacional de Mejoramiento de Maiz y Trigo (CIMMYT), and International Center for Agricultural Research in the Dry Areas (ICARDA), and landraces from discrete regions of the world showed low to moderate phenotypic correlation for disease response amongst genotypes when evaluated with historical and contemporary isolates at two locations across 3 years in Western Australia (WA). Significant (P < 0.001) genotype-by-environment interactions were detected regardless of same or different isolates used as an inoculum source. Joint regression analysis identified 19 genotypes that consistently expressed low disease severity under infection with different isolates in multi-locations. The CIMMYT inbred lines, 30ZJN09 and ZJN12 Qno25, were particularly pertinent as they had low SNB response and highest trait stability at two locations across 3 years. Genome wide association studies detected 20 QTL associated with SNB resistance on chromosomes 1A, 1B, 4B, 5A, 5B, 6A, 7A, 7B, and 7D. QTL on chromosomes 1B and 5B were previously reported in similar genomic regions. Multiple QTL were identified on 1B, 5B, 6A, and 5A and detected in response to SNB infection against different isolates and specific environments. Known SnTox-Snn interactions were either not evident or variable across WA environments and SNB response may involve other multiple complex biological mechanisms.

GWAS-Assisted Genomic Prediction to Predict Resistance to Septoria Tritici Blotch in Nordic Winter Wheat at Seedling Stage

Septoria tritici blotch (STB) disease caused by Zymoseptoria tritici is one of the most damaging diseases of wheat causing significant yield losses worldwide. Identification and employment of resistant germplasm is the most cost-effective method to control STB. In this study, we characterized seedling stage resistance to STB in 175 winter wheat landraces and old cultivars of Nordic origin. The study revealed significant (p < 0.05) phenotypic differences in STB severity in the germplasm. Genome-wide association analysis (GWAS) using five different algorithms identified ten significant markers on five chromosomes. Six markers were localized within a region of 2 cM that contained seven candidate genes on chromosome 1B. Genomic prediction (GP) analysis resulted in a model with an accuracy of 0.47. To further improve the prediction efficiency, significant markers identified by GWAS were included as fixed effects in the GP model. Depending on the number of fixed effect markers, the prediction accuracy improved from 0.47 (without fixed effects) to 0.62 (all non-redundant GWAS markers as fixed effects), respectively. The resistant genotypes and single-nucleotide polymorphism (SNP) markers identified in the present study will serve as a valuable resource for future breeding for STB resistance in wheat. The results also highlight the benefits of integrating GWAS with GP to further improve the accuracy of GP.

TaAPO-A1, an ortholog of rice ABERRANT PANICLE ORGANIZATION 1, is associated with total spikelet number per spike in elite European hexaploid winter wheat (Triticum aestivum L.) varieties

We dissected the genetic basis of total spikelet number (TSN) along with other traits, viz. spike length (SL) and flowering time (FT) in a panel of 518 elite European winter wheat varieties. Genome-wide association studies (GWAS) based on 39,908 SNP markers revealed highly significant quantitative trait loci (QTL) for TSN on chromosomes 2D, 7A, and 7B, for SL on 5A, and FT on 2D, with 2D-QTL being the functional marker for the gene Ppd-D1. The physical region of the 7A-QTL for TSN revealed the presence of a wheat ortholog (TaAPO-A1) to APO1–a rice gene that positively controls the spikelet number on the panicles. Interspecific analyses of the TaAPO-A1 orthologs showed that it is a highly conserved gene important for floral development and present in a wide range of terrestrial plants. Intraspecific studies of the TaAPO-A1 across wheat genotypes revealed a polymorphism in the conserved F-box domain, defining two haplotypes. A KASP marker developed on the polymorphic site showed a highly significant association of TaAPO-A1 with TSN, explaining 23.2% of the total genotypic variance. Also, the TaAPO-A1 alleles showed weak but significant differences for SL and grain yield. Our results demonstrate the importance of wheat sequence resources to identify candidate genes for important traits based on genetic analyses.

Precision Phenotyping Reveals Novel Loci for Quantitative Resistance to Septoria Tritici Blotch

Accurate, high-throughput phenotyping for quantitative traits is a limiting factor for progress in plant breeding. We developed an automated image analysis to measure quantitative resistance to septoria tritici blotch (STB), a globally important wheat disease, enabling identification of small chromosome intervals containing plausible candidate genes for STB resistance. 335 winter wheat cultivars were included in a replicated field experiment that experienced natural epidemic development by a highly diverse but fungicide-resistant pathogen population. More than 5.4 million automatically generated phenotypes were associated with 13,648 SNP markers to perform the GWAS. We identified 26 chromosome intervals explaining 1.9-10.6% of the variance associated with four independent resistance traits. Sixteen of the intervals overlapped with known STB resistance intervals, suggesting that our phenotyping approach can identify simultaneously (i.e., in a single experiment) many previously defined STB resistance intervals. Seventeen of the intervals were less than 5 Mbp in size and encoded only 173 genes, including many genes associated with disease resistance. Five intervals contained four or fewer genes, providing high priority targets for functional validation. Ten chromosome intervals were not previously associated with STB resistance, perhaps representing resistance to pathogen strains that had not been tested in earlier experiments. The SNP markers associated with these chromosome intervals can be used to recombine different forms of quantitative STB resistance that are likely to be more durable than pyramids of major resistance genes. Our experiment illustrates how high-throughput automated phenotyping can accelerate breeding for quantitative disease resistance.