Next Article in Journal
New Virus Variant Detection Based on the Optimal Natural Metric
Previous Article in Journal
Pathway-Based Mendelian Randomization for Pre-Infection IL-6 Levels Highlights Its Role in Coronavirus Disease
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Screening of CIMMYT and South Asian Bread Wheat Germplasm Reveals Marker–Trait Associations for Seedling Resistance to Septoria Nodorum Blotch

by
Rupsanatan Mandal
1,2,
Xinyao He
3,
Gyanendra Singh
4,
Muhammad Rezaul Kabir
5,
Arun Kumar Joshi
6,7 and
Pawan Kumar Singh
3,*
1
Visiting Scientist, International Maize and Wheat Improvement Center (CIMMYT), Texcoco 56237, Mexico
2
Department of Genetics and Plant Breeding, Uttar Banga Krishi Viswavidyalaya, Cooch Behar 736165, India
3
International Maize and Wheat Improvement Centre, Texcoco 56237, Mexico
4
ICAR-Indian Institute of Wheat and Barley Research, Karnal 132001, India
5
Bangladesh Wheat and Maize Research Institute, Dinajpur 5200, Bangladesh
6
International Maize and Wheat Improvement Center (CIMMYT)-India Office, New Delhi 110012, India
7
Borlaug Institute for South Asia, New Delhi 110012, India
*
Author to whom correspondence should be addressed.
Genes 2024, 15(7), 890; https://doi.org/10.3390/genes15070890
Submission received: 14 June 2024 / Revised: 2 July 2024 / Accepted: 5 July 2024 / Published: 7 July 2024
(This article belongs to the Special Issue Genetics of Disease Resistance in Wheat)

Abstract

:
Wheat (Triticum aestivum L.) production is adversely impacted by Septoria nodorum blotch (SNB), a fungal disease caused by Parastagonospora nodorum. Wheat breeders are constantly up against this biotic challenge as they try to create resistant cultivars. The genome-wide association study (GWAS) has become an efficient tool for identifying molecular markers linked with SNB resistance. This technique is used to acquire an understanding of the genetic basis of resistance and to facilitate marker-assisted selection. In the current study, a total of 174 bread wheat accessions from South Asia and CIMMYT were assessed for SNB reactions at the seedling stage in three greenhouse experiments at CIMMYT, Mexico. The results indicated that 129 genotypes were resistant to SNB, 39 were moderately resistant, and only 6 were moderately susceptible. The Genotyping Illumina Infinium 15K Bead Chip was used, and 11,184 SNP markers were utilized to identify marker–trait associations (MTAs) after filtering. Multiple tests confirmed the existence of significant MTAs on chromosomes 5B, 5A, and 3D, and the ones at Tsn1 on 5B were the most stable and conferred the highest phenotypic variation. The resistant genotypes identified in this study could be cultivated in South Asian countries as a preventative measure against the spread of SNB. This work also identified molecular markers of SNB resistance that could be used in future wheat breeding projects.

1. Introduction

Septoria nodorum blotch (SNB) is a significant foliar fungal disease caused by P. nodorum that affects wheat production globally [1]. This fungus prefers the moist and warm areas of the world [2,3], and recently SNB has been reported in India [4,5]. This disease leads to substantial reductions in yield losses and creates an obstacle for wheat breeders striving to develop resistant genotypes [6]. The pathogen mostly infects leaves, resulting in necrotic lesions that reduce their photosynthetic capability and subsequently affect grain yield [7]. The pathogen P. nodorum has the ability to generate numerous virulence factors, including SnToxA, SnTox1, and SnTox3, which facilitate its ability to invade wheat plants [8,9]. These factors consist of enzymes that degrade the cell walls of plants, toxins that harm plant cells, and other chemicals that promote fungal growth [1,10]. The genetic heterogeneity of the pathogen enables it to bypass plant resistance mechanisms and adapt to different wheat cultivars [11,12,13,14], while wheat plants possess an intrinsic genetic capability for resistance against SNB [15,16]. Many investigations of SNB resistance, which is often regulated by several genes with minor individual effects, have shown that it is known as quantitative resistance [17,18,19,20,21]. Through the use of advanced genotyping technologies, researchers have conducted an in-depth study of the genetic factors that contribute to SNB resistance. Genome-wide association studies (GWASs) have proved their efficacy in identifying genetic markers associated with different phenotypes, such as SNB resistance, which offers useful insights into the genetic basis of resistance and facilitates marker-assisted selection [22,23,24,25,26,27,28]. Several GWASs conducted on wheat have identified multiple quantitative trait loci (QTL) that are associated with resistance to SNB. These QTLs are dispersed throughout various chromosomes, including 1A, 2A, 3A, 4A, 5A, 6A, 7A, 1B, 3B, 4B, 5B, 6B, 2D, 5D, and 7D [20,21,29,30,31]. GWAS has provided important insights, but issues including genetic variability, population structure, and linkage disequilibrium need to be properly addressed [20,21,22,23,24,25]. The molecular mechanism underlying SNB resistance in wheat may be better understood by combining GWAS data with functional genomics and transcriptomics information. The identification of novel loci and markers linked with SNB resistance genes through GWASs holds great promise for accelerating the development of resistant varieties of wheat, contributing to sustainable wheat production. The aims of the study were to determine the genotypes that are resistant to SNB and to perform a thorough GWAS to find the specific locations that are linked with SNB resistance utilizing 174 wheat accessions from India, Nepal, Bangladesh, and CIMMYT-Mexico.

2. Materials and Methods

2.1. Plant Materials

One hundred seventy-four spring wheat genotypes were used for the current investigation. Among the accessions, 97 originated from CIMMYT, Mexico, 30 from India, 28 from Nepal, and the remaining 19 from Bangladesh (Supplementary Table S1), which are elite breeding lines and modern varieties in the respective countries.

2.2. Disease Screening

Three separate tests were implemented in a greenhouse at CIMMYT-Mexico to evaluate disease severity in wheat seedlings. The experiments were carried out using an RCBD (randomized complete block design), with three replications. The experiments were conducted at temperatures of 22 °C during the day and 18 °C at night, with a 16 h photoperiod. The tests utilized plastic containers as experimental units, with four plants for each entry. The resistant control genotype used was “Erik”, while the susceptible control genotype used was “Glenlea”(Supplementary Table S3). The inoculum was made from the Mexican P. nodorum strain MexSn4 (CIMMYT ID: CIMFU-463), which was identified as a ToxA producer based on inoculation tests, filtering experiments, and the ToxA marker assay. Conidiospores were prepared from V8 Juice–PDA medium and were adjusted to 1 × 107 spores mL−1 for inoculation, following the protocols described by Hu et al. [32]. Artificial inoculation was carried out when the second leaf had fully grown, at approximately 14 DAP (days after planting). This was carried out by using a hand sprayer to apply the inoculum until it covered the entire plant, with an approximate volume of 0.5 mL per plant. After the leaves had dried, they were moved to a humid chamber with relative humidity of 100% and a temperature of 20 °C to promote infection. After a period of 24 h, the plants were placed back onto the greenhouse bench. The disease scores were evaluated using a linear scale ranging from 1 to 5, where 1 represented the least infected and 5 represented the most infected [32,33,34].

2.3. Statistical Analysis

Analysis of variance (ANOVA) was conducted using R software version 4.3.1 with the “metan” package. Additionally, stability and GGE biplot analyses were performed using the “plant breeding” and “GGEBiplotGUI” packages [35,36].

2.4. Genotyping

The panel was genotyped at Trait Genetics GmbH, Gatersleben, Germany, using the Illumina Infinium 15 K Bead Chip, which produced a total of 16,028 SNPs. The genotypic data were filtered to exclude SNP markers with minor allele frequencies less than 0.05, unknown chromosome positions, or more than 10% missing SNP data, resulting in 11,184 SNP markers that were included in the subsequent study. Genomic locations of the SNP markers were determined either from the consensus map [37] or via BLAST searches against Chinese Spring reference genome v.1.0 using Phytozome v13 [38].

2.5. Population Structure and Principal Coordinate Analysis

The number of subpopulations among the 174 wheat genotypes was estimated using STRUCTURE software v2.3.4, which is a model-based Bayesian cluster analysis tool [39]. The number of presumed groups was set between k = 1 and 10, and the admixture model was employed. Each of the five independent replications of the analysis involved 50,000 burn-in repeats and 10,000 MCMC (Markov Chain Monte Carlo) iterations. The statistical measure ΔK, which is derived from the rate of change in the log-likelihood of the provided data, was employed with Structure Harvester to ascertain the optimal number of clusters [40]. Structure Plot v2.0 was used to create the structure bar plot with the most clusters possible [41]. PCoA was analyzed using the GenAIEx 6.5 software [42].

2.6. Kinship and Linkage Disequilibrium Analysis

TASSEL v5.069 was used to analyze all 11,184 SNP markers to generate a kinship matrix as well as clusters among the individual genotypes. A kinship matrix heat map was generated using the R program v4.3.1 [35]. Linkage disequilibrium (LD) of the panel was examined using TASSEL v5.069 [43] and the 11,184 SNP markers. By graphing the R2 values against the physical distance (Mb) between the markers using R software, the LD decay distances were obtained for the entire genome and individual chromosomes.

2.7. GWAS Analysis for SNB

GLM, MLM, and FarmCPU were employed to identify marker–trait associations. The FarmCPU model was carried out using the R software package rMVP [44], whereas the PCA-based GLM model and the MLM model, which consider both Kinship (K matrix) and population structure (Q value), were executed in Tassel v5.069 [43]. The three experiments were each subjected to a separate GWAS and then examined using the grand mean data across experiments. Significant markers in all models were declared with a threshold of p < 0.001.

3. Results

3.1. Phenotypic Evaluation

The distribution of disease scores in the panel throughout the three experiments was depicted in a box plot diagram (Figure 1A). The examined genotypes showed a skewed distribution towards the resistance side, and around 74% of accessions exhibited resistant reactions (disease scores between 1.0 and 2.0), 39 genotypes showed moderately resistant reactions (2.1–3.0), and only 6 genotypes were moderately susceptible to susceptible, with disease scores higher than 3.1 (Figure 1B, Supplementary Table S1).
Two genotypes (CIM50 and CIM82), which were the most resistant, did not exhibit any infection symptoms in all three experiments, and many other genotypes showed a very low degree of infection symptoms. The 174 wheat genotypes demonstrated a grand mean SNB score of 1.71, and significant correlations were found among individual experiments, ranging from 0.57 (Experiment 1 vs. Experiment 3) to 0.66 (Experiment 2 vs. Experiment 3). The Bartlett test revealed homogeneity of variances among experiments, with p < 0.001. ANOVA revealed significant effects in “Genotype” and “Experiment”, as well as their interaction (Table 1). High broad-sense heritability (h2bs) estimates were obtained for the three experiments, ranging from 91% (Experiment 3) to 95% (Experiment 1).
Based on the SNB scores, the three experiments were divided into two sectors with various resistant and susceptible genotypes in the “which-won-where” view of the GGE biplot (Figure 2). The first two principal components accounted for 90.24% of the whole G+GE (Genotype + Genotype × Experiment) variation.
The three experiments had distinct contributions to the observed difference in SNB resistance, i.e., the first experiment (E1) was more different from the other two experiments (E2 and E3) and made significant contributions to the genotype-by-environment interaction (G × E). Additionally, E1 was situated farther from the origin than E2 and E3, indicating a higher phenotypic variation among genotypes (Figure 2). The tested genotypes displayed varying principal component scores, implying variable genotype-by-environment performance. The lines NPL-12, NPL-17, CIM-22, CIM-51, CIM-83, CIM-87, CIM-97, NPL-8, NPL-28, CIM-39, and NPL-6 demonstrated distinct adaptations for specific environments, as indicated by their performance in experiments E1, E2, and E3. On the contrary, genotypes like CIM-50 and CIM-82 had the lowest principal component values, indicating a reduced genotype-by-environment interaction and increased stability (Figure 2).

3.2. Population Structure and Principal Coordinate Analysis

Based on the ΔK statistic, population structure analysis identified two subpopulations among the 174 genotypes (Figure 3A,B), with 112 and 62 members in subpopulations SP1 and SP2, respectively (Supplementary Table S1). The majority of CIMMYT genotypes (87 out of 97) were found in SP1, whereas most genotypes from Bangladesh (17 out of 19) and those from Nepal (20 out of 28) were found in SP2. Indian accessions did not show obvious differences in the two subpopulations (16 in SP1 and 14 in SP2). Genotypes in SP1 generally exhibited better resistance than those from SP2, according to the mean SNB scores of 1.68 for SP1 and 1.76 for SP2. PCoA analysis also divided the panel into two groups, which were very similar to the two subpopulations identified in the structure analysis (Figure 3C).

3.3. Kinship and Linkage Disequilibrium Analysis

Using the 11,184 SNPs, kinship analysis virtually divided the 174 wheat genotypes into two groups, which were very similar to those in the STRUCTURE and PCoA analyses: CIMMYT vs. non-CIMMYT genotypes (as shown in Figure 4). As for the LD analysis, the critical R2 value was estimated at 2.24 Mb for the entire genome, and at 2.23 Mb, 4.69 Mb, and 1.78 Mb for the A, B, and D genomes, respectively (Supplementary Figure S1).

3.4. Marker Trait Association (MTA) for SNB

Of the three GWAS algorithms tested in this study, namely GLM, MLM, and FarmCPU, MLM fit the data the best (Supplementary Figure S2) and was thus chosen for subsequent analysis. A total of 11 significant MTAs for SNB resistance were identified (Figure 5), and the SNPs Tdurum_contig25513_123, tplb0027f13_1493, and Tdurum_contig31131_198 on chromosome 5B were the most stable and were detected in two individual experiments. The remaining seven MTAs, which were significant in only one of the experiments and the mean data, were distributed on chromosomes 3D, 5A, and 5B (Figure 5). The phenotypic variation explained by these significant markers ranged from 7.35 to 15.66%, with the highest value from tplb0027f13_1493, followed by 13.30% from IACX9261, and both of these markers are on the long arm of chromosome 5B. There are additional MTAs on this chromosome arm that were significant in this study, and four of them, tplb0027f13_1493, Tdurum_contig25513_123, Tdurum_contig31131_198T, and IACX9261, exhibited significant phenotypic effects (Table 2). The genotypes in this panel that exhibited high resistance had a resistance allele for all four markers, as anticipated.

3.5. Candidate Genes for the Significant MTAs

Nucleotide sequences of the significant SNP markers were used in BLAST searches against the Chinese Spring reference genome v2.2 available using the Phytozome v13 online tool to identify putative genes associated with disease resistance mechanisms in plants. We employed SNPs for the 11 significant MTAs to explore the putative genes. Each MTA was analyzed using a 2 Mb window to find the putative genes, resulting in a total of eight genes in the genomic intervals (Supplementary Table S2), all of which have known functions relevant to disease resistance.

4. Discussion

The potential yield loss caused by SNB is estimated to be up to 16% worldwide, which warrants great efforts regarding resistance breeding and disease management for this disease [45]. Exploring resistant sources in the current wheat panel and employing them in wheat breeding could be a useful approach to managing diseases. In this study, we conducted greenhouse screening for seedling SNB resistance, which avoided the influence of other naturally occurring foliar diseases that mimic SNB symptoms, like spot blotch and tan spot; in addition, greenhouse experiments provide the optimized light and humidity requirements for SNB infection that are often not available for field experiments [32,33,34]. As reported by Thapa et al. [46], greenhouse screening against SNB was effective for the identification of resistant genotypes as well as for the GWAS study.
Using the same panel, Phuke et al. [47] reported high levels of resistance of a CIMMYT germplasm to tan spot, and the present study has drawn a similar conclusion for resistance to SNB (around 74% of genotypes were resistant). However, in a previous study on greenhouse SNB screening of CIMMYT germplasms, only 51 out of 385 lines tested were classified as resistant [48], indicating a significant increase in the frequency of SNB-resistant genotypes in recent CIMMYT germplasms. A possible explanation for this could be the continuous selection for SNB resistance, and the consequent decreased frequency of Tsn1, a susceptible gene for multiple wheat diseases, e.g., SNB, tan spot, and spot blotch, which will be further discussed later.
A moderate population structure was found in the current study, which could be ascribed to the fact that CIMMYT wheat germplasms have been widely distributed and utilized globally [47,49], and the frequent and historical germplasm exchange between CIMMYT and the South Asian countries resulted in the reduced population structure observed in this study. LD decay for the A and B subgenomes was smaller than for other spring wheat populations [50,51], and LD for subgenome D was significantly lower than for subgenomes A and B.
In our study, the genomic region on the long arm of chromosome 5B represented by the markers Tdurum_contig25513_123, tplb0027f13_1493, and IACX9261 was the most stable, which must represent the Tsn1 gene that plays a substantial role in SNB resistance [28]. Similar conclusions were made by Al Tameemi et al. [38] using 274 US wheat accessions and by Francki et al. [52] using 232 global wheat accessions. In our study, three MTAs on the short arm of chromosome 5B (BS00091519_51, GENE-3324_338, and Excalibur_c47452_183) were significantly associated with SNB disease, which, based on their physical positions, must represent the Snn3-B1 locus, which has a significant impact on SNB resistance, according to Downie et al. [26]. It is noteworthy that this genomic region was found to be significant in one of our recent studies on SNB resistance for a panel of 296 Indian genotypes [53]. The significant MTAs identified on chromosomes 3D (AX-94978939) and 5B (Tdurum_contig31131_198) do not match with previously identified QTL; therefore, these may be novel genomic regions for SNB resistance, but further validation is needed before they are utilized in breeding.
Resistance breeding in wheat is a tedious task, mainly due to the quantitative nature of resistance, which is true for leaf spotting diseases including SNB, as well as for adult plant resistance against rusts and powdery mildew [54]. This, together with other important quantitative traits like yield potential and abiotic stress resistance, often slows progress regarding the genetic gain of wheat breeding. Nevertheless, the existence of multi-disease resistance genes/loci provides breeders an opportunity to incorporate a single gene/locus for multiple traits, like Lr34/Yr18/Sr57/Pm38/Sb1 and Lr46/Yr29/Sr58/Pm39, which should be fully utilized by breeders. Tsn1 is one such gene, which was initially associated with resistance against SNB and tan spot, and later spot blotch, too [53]. For the time being, the former two diseases are not important in South Asia, but spot blotch has been a major disease in the region and accounts for significant annual yield losses [55]. Considering this, it is strongly recommended that wheat breeders in South Asia eliminate this susceptible gene from their breeding materials to first increase spot blotch resistance and then be prepared for possible future epidemics of SNB and tan spot. If an SNB epidemic occurs in the region in the future, which is not unlikely considering the changing climate and the wide adoption of conservation agricultural practices in the region, additional SNB resistance genes/loci could be incorporated, among which are the resistance alleles for the markers identified in the current study, especially Snn3-B1 on chromosome 5B, which has shown importance in multiple studies.

5. Conclusions

In conclusion, a panel of 174 wheat genotypes from CIMMYT and South Asian countries was evaluated in this study for seedling resistance against SNB. About 74% of genotypes exhibited resistant reactions, which could be ascribed to the low frequency of the susceptibility gene Tsn1. Using the MLM algorithm for GWAS, 11 MTAs were identified, representing the known susceptibility genes Tsn1 and Snn3-B1, as well as two potentially novel MTAs on chromosomes 3D and 5B that confer SNB resistance. The identified MTAs could be useful for marker-assisted selection in wheat breeding using the resistant genotypes identified in this study as donors.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/genes15070890/s1. Supplementary Table S1: List of genotypes with their phenotypic values in three experiments, and Q values from the structure analysis; Supplementary Table S2: Transcript id and functional annotation of MTAs; Supplementary Table S3: SNB disease scores of check genotypes for three experiments and their pooled; Supplementary Figure S1: LD decay plots for whole genomes and individual chromosome; Supplementary Figure S2: QQ plots for GWAS using GLM, MLM and FarmCPU models.

Author Contributions

P.K.S. designed the research and supervised the experiments. G.S., A.K.J. and M.R.K. provided plant materials. X.H. contributed to the genotyping work, and P.K.S. contributed to the phenotyping work. R.M. and X.H. analyzed the data. R.M. wrote the first draft of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

Financial support from the Indian Council of Agriculture Research (ICAR), India, the Accelerating Genetic Gain (AGG) in Maize and Wheat Project (grant no. INV-003439) funded by the Bill and Melinda Gates Foundation (BMGF), the Foreign and Commonwealth Development Office (FCDO), and the Foundation for Food and Agriculture Research (FFAR), USAID, and One CGIAR Initiatives for conducting this research is gratefully acknowledged.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding author.

Acknowledgments

The first author is grateful for financial support from the ICAR-NAHEB International Travel Grant. Technical support from Nerida Lozano, Elizabeth Blanco Venegas, Mayela Flores Romero, and Mónica Preciado from Wheat Pathology, CIMMYT, is highly acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ANOVAAnalysis of Variance
FarmCPUFixed and random model Circulating Probability Unification
GWASGenome-Wide Association Study
GLMGeneral Linear Model
MLMMultiple Linear Model
MTAMarker Traits Association
QTLQuantitative Trait Loci
rMVPMemory-Efficient, Visualize-Enhanced, Parallel-Accelerated GWAS Tool
SNBSeptoria Nodorum Blotch
SNPSingle Nucleotide Polymorphism

References

  1. Ruud, A.K.; Lillemo, M. Diseases affecting wheat: Septoria nodorum blotch. In Burleigh Dodds Series in Agricultural Science; Burleigh Dodds Science Publishing Limited: Cambridge, UK, 2018; pp. 109–144. [Google Scholar] [CrossRef]
  2. Bhathal, J.S.; Loughman, R.; Speijers, J. Yield reduction in wheat in relation to leaf disease from yellow (tan) spot and Septoria nodorum blotch. Eur. J. Plant Pathol. 2003, 109, 435–443. [Google Scholar] [CrossRef]
  3. Juliana, P.; Singh, R.P.; Singh, P.K.; Crossa, J.; Rutkoski, J.E.; Poland, J.A.; Bergstrom, G.C.; Sorrells, M.E. Comparison of models and whole-genome profiling approaches for genomic-enabled prediction of Septoria tritici blotch, Stagonospora nodorum blotch, and tan spot resistance in wheat. Plant Genome 2017, 10, plantgenome2016-08. [Google Scholar] [CrossRef] [PubMed]
  4. Katoch, S.; Rana, S.K.; Sharma, P.N. Application of PCR based diagnostics in the exploration of Parastagonospora nodorum prevalence in wheat growing regions of Himachal Pradesh. J. Plant Biochem. Biotechnol. 2019, 28, 169–175. [Google Scholar] [CrossRef]
  5. Navathe, S.; He, X.; Kamble, U.; Kumar, M.; Patial, M.; Singh, G.; Singh, G.P.; Joshi, A.K.; Singh, P.K. Assessment of Indian wheat germplasm for Septoria nodorum blotch and tan spot reveals new QTLs conferring resistance along with recessive alleles of Tsn1 and Snn3. Front. Plant Sci. 2023, 14, 1223959. [Google Scholar] [CrossRef] [PubMed]
  6. Ficke, A.; Cowger, G.C.; Bergstrom, G.C.; Guro, B. Understanding yield loss and pathogen biology to improve disease management: Septoria Nodorum Blotch—A case study in wheat. Plant Dis. 2017, 102, 696–707. [Google Scholar] [CrossRef] [PubMed]
  7. Blixt, E.; Olson, Å.; Lindahl, B.; Djurle, A.; Yuen, J. Spatiotemporal variation in the fungal community associated with wheat leaves showing symptoms similar to Stagonospora nodorum blotch. Eur. J. Plant Pathol. 2010, 126, 373–386. [Google Scholar] [CrossRef]
  8. Friesen, T.L.; Holmes, D.J.; Bowden, R.L.; Faris, J.D. ToxA is present in the U.S. Bipolaris sorokiniana population and is a significant virulence factor on wheat harboring Tsn1. Plant Dis. 2018, 102, 2446–2452. [Google Scholar] [CrossRef] [PubMed]
  9. McDonald, M.C.; Ahren, D.; Simpfendorfer, S.; Milgate, A.; Solomon, P.S. The discovery of the virulence gene ToxA in the wheat and barley pathogen Bipolaris sorokiniana. Mol. Plant Pathol. 2018, 19, 432–439. [Google Scholar] [CrossRef] [PubMed]
  10. Peters Haugrud, A.R.; Zhang, Z.; Richards, J.K.; Friesen, T.L.; Faris, J.D. Genetics of variable disease expression conferred by inverse gene-for-gene interactions in the wheat-Parastagonospora nodorum pathosystem. Plant Physiol. 2019, 180, 420–434. [Google Scholar] [CrossRef]
  11. Murphy, N.E.; Loughman, R.; Appels, R.; Lagudah, E.S.; Jones, M.G.K. Genetic variability in a collection of Stagonospora nodorum isolates from Western Australia. Aust. J. Agric. Res. 2000, 51, 679–684. [Google Scholar] [CrossRef]
  12. McDonald, B.A.; Linde, C. Pathogen population genetics, evolutionary potential, and durable resistance. Annu. Rev. Phytopathol. 2002, 40, 349–379 . [Google Scholar] [CrossRef] [PubMed]
  13. Stukenbrock, E.H.; Banke, S.; McDonald, B.A. Global migration patterns in the fungal wheat pathogen Phaeosphaeria nodorum. Mol. Ecol. 2006, 15, 2895–2904. [Google Scholar] [CrossRef]
  14. McDonald, M.C.; Razavi, M.; Friesen, T.L.; Brunner, P.C.; McDonald, B.A. Phylogenetic and population genetic analyses of Phaeosphaeria nodorum and its pathosystem parallels that of the wheat–tan spot system. Genome 2012, 49, 1265–1273. [Google Scholar] [CrossRef]
  15. Bostwick, D.E.; Ohm, H.W.; Shaner, G. Inheritance of Septoria-glume blotch resistance in wheat. Crop Sci. 1993, 33, 439–443. [Google Scholar] [CrossRef]
  16. Wicki, W.; Winzeler, M.; Schmid, J.E.; Stamp, P.; Messmer, M. Inheritance of resistance to leaf and glume blotch caused by Septoria nodorum Berk. in winter wheat. Theor. Appl. Genet. 1999, 99, 1265–1272. [Google Scholar] [CrossRef]
  17. Adhikari, T.B.; Jackson, E.W.; Gurung, S.; Hansen, J.M.; Bonman, J.M. Association mapping of quantitative resistance to Phaeosphaeria nodorum in spring wheat landraces from the USDA National Small Grains Collection. Phytopathology 2011, 101, 1301–1310. [Google Scholar] [CrossRef]
  18. Gurung, S.; Mamidi, S.; Bonman, J.M.; Xiong, M.; Brown-Guedira, G.; Adhikari, T.B. Genome-wide association study reveals novel quantitative trait loci associated with resistance to multiple leaf spot diseases of spring wheat. PLoS ONE 2014, 9, e108179. [Google Scholar] [CrossRef]
  19. Jighly, A.; Alagu, M.; Makdis, F.; Singh, M.; Singh, S.; Emebiri, L.C.; Ogbonnaya, F.C. Genomic regions conferring resistance to multiple fungal pathogens in synthetic hexaploid wheat. Mol. Breed. 2016, 36, 1–9. [Google Scholar] [CrossRef]
  20. Liu, Z.; El-Basyoni, I.; Kariyawasam, G.; Zhang, G.; Fritz, A.; Hansen, J.; Marais, F.; Friskop, A.; Chao, S.; Akhunov, E.; et al. Evaluation and association mapping of resistance to tan spot and Stagonospora nodorum blotch in adapted winter wheat germplasm. Plant Dis. 2015, 99, 1333–1341. [Google Scholar] [CrossRef]
  21. Ruud, A.K.; Dieseth, J.A.; Ficke, A.; Furuki, E.; Phan, H.T.T.; Oliver, R.P.; Tan, K.C.; Lillemo, M. Genome-wide association mapping of resistance to Septoria nodorum leaf blotch in a Nordic spring wheat collection. Plant Genome 2019, 12, 180105. [Google Scholar] [CrossRef]
  22. Breseghello, F.; Sorrells, M.E. Association analysis as a strategy for improvement of quantitative traits in plants. Crop Sci. 2006, 46, 1323–1330. [Google Scholar] [CrossRef]
  23. Gupta, P.K.; Kulwal, P.L.; Jaiswal, V. Association mapping in crop plants: Opportunities and challenges. Adv. Genet. 2014, 85, 109–147. [Google Scholar] [PubMed]
  24. Cockram, J.; Scuderi, A.; Barber, T.; Furuki, E.; Gardner, K.A.; Gosman, N.; Kowalczyk, R.; Phan, H.P.; Rose, G.A.; Tan, K.C.; et al. Fine-mapping the wheat Snn1 locus conferring sensitivity to the Parastagonospora nodorum necrotrophic effector SnTox1 using an eight founder multiparent advanced generation inter-cross population. G3 (Bethesda) 2015, 5, 2257–2266. [Google Scholar] [CrossRef] [PubMed]
  25. Bernardo, R. Bandwagons I, too, have known. Theor. Appl. Genet. 2016, 129, 2323–2332. [Google Scholar] [CrossRef] [PubMed]
  26. Downie, R.C.; Bouvet, L.; Furuki, E.; Gosman, N.; Gardner, K.A.; Mackay, I.J.; Campos Mantello, C.; Mellers, G.; Phan, H.T.; Rose, G.A.; et al. Assessing European wheat sensitivities to Parastagonospora nodorum necrotrophic effectors and fine-mapping the Snn3-B1 locus conferring sensitivity to the effector SnTox3. Front. Plant Sci. 2018, 9, 881. [Google Scholar] [CrossRef]
  27. Lin, M.; Corsi, B.; Ficke, A.; Tan, K.C.; Cockram, J.; Lillemo, M. Genetic mapping using a wheat multi-founder population reveals a locus on chromosome 2A controlling resistance to both leaf and glume blotch caused by the necrotrophic fungal pathogen Parastagonospora nodorum. Theor. Appl. Genet. 2020, 133, 785–808. [Google Scholar] [CrossRef] [PubMed]
  28. Lin, M.; Stadlmeier, M.; Mohler, V.; Tan, K.C.; Ficke, A.; Cockram, J.; Lillemo, M. Identification and cross-validation of genetic loci conferring resistance to Septoria nodorum blotch using a German multi-founder winter wheat population. Theor. Appl. Genet. 2021, 134, 125–142. [Google Scholar] [CrossRef] [PubMed]
  29. Tommasini, L.; Schnurbusch, T.; Fossati, D.; Mascher, F.; Keller, B. Association mapping of Stagonospora nodorum blotch resistance in modern European winter wheat varieties. Theor. Appl. Genet. 2007, 115, 697–708. [Google Scholar] [CrossRef]
  30. Phan, H.T.; Rybak, K.; Furuki, E.; Breen, S.; Solomon, P.S.; Oliver, R.P.; Tan, K.C. Differential effector gene expression underpins epistasis in a plant fungal disease. Plant J. 2016, 87, 343–354. [Google Scholar] [CrossRef]
  31. Halder, J.; Zhang, J.; Ali, S.; Sidhu, J.S.; Gill, H.S.; Talukder, S.K.; Kleinjan, J.; Turnipseed, B.; Sehgal, S.K. Mining and genomic characterization of resistance to tan spot, Stagonospora nodorum blotch (SNB), and Fusarium head blight in Watkins core collection of wheat landraces. BMC Plant Biol. 2019, 19, 480. [Google Scholar] [CrossRef]
  32. Hu, W.; He, X.; Dreisigacker, S.; Sansaloni, C.P.; Juliana, P.; Singh, P.K. A wheat chromosome 5AL region confers seedling resistance to both tan spot and Septoria nodorum blotch in two mapping populations. Crop J. 2019, 7, 809–818. [Google Scholar] [CrossRef]
  33. Feng, J.; Ma, H.; Hughes, G.R. Genetics of resistance to Stagonospora nodorum blotch of hexaploid wheat. Crop Sci. 2004, 44, 2043–2048. [Google Scholar] [CrossRef]
  34. Singh, P.K.; Mergoum, M.; Adhikari, T.B.; Kianian, S.F.; Elias, E.M. Chromosomal location of genes for seedling resistance to tan spot and Stagonospora nodorum blotch in tetraploid wheat. Euphytica 2007, 155, 27–34. [Google Scholar] [CrossRef]
  35. Kumar, V.; Kumar, S.; Kharub, A.S.; Verma, R.P.S.; Singh, G.P. GGE Biplot Model for Environmental Delineation and Identification of Resistant Sources for Stripe Rust and Spot Blotch in India; Rabat, Morocco. 2017. Avaliable online: https://repo.mel.cgiar.org/handle/20.500.11766/8026 (accessed on 31 December 2023).
  36. Chattopadhyay, N.; Mandal, R.; Roy, A.; Bhattacharya, P.M.; Chowdhury, A.K. Correction to: Assessment of wheat genotypes based on genotype-by-environment interaction for durable resistance to spot blotch disease in hot spot. Cereal Res. Commun. 2022, 50, 103. [Google Scholar] [CrossRef]
  37. Wang, S.; Wong, D.; Forrest, K.; Allen, A.; Chao, S.; Huang, B.E.; Maccaferri, M.; Salvi, S.; Milner, S.G.; Cattivelli, L.; et al. Characterization of polyploid wheat genomic diversity using a high-density 90 000 single nucleotide polymorphism array. Plant Biotechnol. J. 2014, 12, 787–796. [Google Scholar] [CrossRef] [PubMed]
  38. AlTameemi, R.; Gill, H.S.; Ali, S.; Ayana, G.; Halder, J.; Sidhu, J.S.; Gill, U.S.; Turnipseed, B.; Hernandez, J.L.G.; Sehgal, S.K. Genome-wide association analysis permits characterization of Stagonospora nodorum blotch (SNB) resistance in hard winter wheat. Sci. Rep. 2021, 11, 12570. [Google Scholar] [CrossRef] [PubMed]
  39. Pritchard, J.K.; Stephens, M.; Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 2000, 155, 945–959. [Google Scholar] [CrossRef]
  40. Earl, D.A.; Vonholdt, B. Structure Harvester: A website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv. Genet. Resour. 2012, 4, 359–361. [Google Scholar] [CrossRef]
  41. Ramasamy, R.K.; Ramasamy, S.; Bindroo, B.B.; Naik, V.G. STRUCTURE PLOT: A program for drawing elegant STRUCTURE bar plots in user friendly interface. Springerplus 2014, 3, 1–3. [Google Scholar] [CrossRef]
  42. Peakall, R.; Smouse, P.E. GenAlEx 6.5: Genetic analysis in Excel. Population genetic software for teaching and research—An update. Bioinformatics 2012, 28, 2537–2539. [Google Scholar] [CrossRef]
  43. Bradbury, P.; Zhang, Z.; Kroon, D.; Casstevens, T.; Ramdoss, Y.; Buckler, E. TASSEL: Software for association mapping of complex traits in diverse samples. Bioinformatics 2007, 23, 2633–2635. [Google Scholar] [CrossRef] [PubMed]
  44. Yin, L.; Zhang, H.; Tang, Z.; Xu, J.; Yin, D.; Zhang, Z.; Yuan, X.; Zhu, M.; Zhao, S.; Li, X.; et al. rMVP: A Memory-efficient, visualization-enhanced, and parallel-accelerated tool for genome-wide association study. Genom. Proteom. Bioinform. 2021, 19, 619–628. [Google Scholar] [CrossRef] [PubMed]
  45. Ficke, A.; Dieseth, J.A.; Kim, M.O.; Lillemo, M. Bladsjukdommer i norsk hvete. Forekomst, betydning of tiltak. In Jord- og Plantekultur 2018/NIBIO BOK; Norwegian Institute of Bioeconomy Research: Ås, Norway, 2018; pp. 108–112. [Google Scholar]
  46. Thapa, R.; Brown-Guedira, G.; Ohm, H.W.; Mateos-Hernandez, M.; Wise, K.A.; Goodwin, S.B. Determining the order of resistance genes against Stagonospora nodorum blotch, Fusarium head blight and stem rust on wheat chromosome arm 3BS. BMC Res Notes 2016, 9, 58. [Google Scholar] [CrossRef] [PubMed]
  47. Phuke, R.M.; He, X.; Juliana, P.; Bishnoi, S.K.; Singh, G.P.; Kabir, M.R.; Roy, K.K.; Joshi, A.K.; Singh, R.P.; Singh, P.K. Association mapping of seedling resistance to tan spot (Pyrenophora tritici-repentis race 1) in CIMMYT and South Asian wheat germplasm. Front Plant Sci. 2020, 11, 1309. [Google Scholar] [CrossRef] [PubMed]
  48. Singh, P.K.; Duveiller, E.; Singh, R.P. Evaluation of CIMMYT germplasm for resistance to leaf spotting diseases of wheat Czech. J. Genet. Plant Breed. 2011, 47, S102–S108. [Google Scholar] [CrossRef]
  49. Ortiz, R.; Sayre, K.D.; Govaerts, B.; Gupta, R.; Subbarao, G.V.; Ban, T.; Hodson, D.; Dixon, J.M.; Ortiz-Monasterio, J.I.; Reynolds, M. Climate change: Can wheat beat the heat? Agric. Ecosyst. Environ. 2008, 126, 46–58. [Google Scholar] [CrossRef]
  50. Chao, S.; Dubcovsky, J.; Dvorak, J.; Luo, M.-C.; Baenziger, S.P.; Matnyazov, R.; Clark, D.R.; Talbert, L.E.; Anderson, J.A.; Dreisigacker, S.; et al. Population- and genome-specific patterns of linkage disequilibrium and SNP variation in spring and winter wheat (Triticum aestivum L.). BMC Genom. 2010, 11, 727. [Google Scholar] [CrossRef] [PubMed]
  51. Bajgain, P.; Rouse, M.N.; Bulli, P.; Bhavani, S.; Gordon, T.; Wanyera, R.; Njau, P.N.; Legesse, W.; Anderson, J.A.; Pumphrey, M.O. Association mapping of North American spring wheat breeding germplasm reveals loci conferring resistance to Ug99 and other African stem rust races. BMC Plant Biol. 2015, 15, 249. [Google Scholar] [CrossRef]
  52. Francki, M.G.; Walker, E.; McMullan, C.J.; Morris, W.G. Evaluation of Septoria nodorum blotch (SNB) resistance in glumes of wheat (Triticum aestivum L.) and the genetic relationship with foliar disease response. Front Genet. 2021, 12, 681768. [Google Scholar] [CrossRef]
  53. Navathe, S.; Yadav, P.S.; Chand, R.; Mishra, V.K.; Vasistha, N.K.; Meher, P.K.; Joshi, A.K.; Gupta, P.K. ToxA-Tsn1 interaction for spot blotch susceptibility in Indian wheat: An example of inverse gene-for-gene relationship. Plant Dis. 2020, 104, 71–81. [Google Scholar] [CrossRef]
  54. Singh, R.P.; Singh, P.K.; Rutkoski, J.; Hodson, D.P.; He, X.; Jørgensen, L.N.; Hovmøller, M.S.; Huerta-Espino, J. Disease impact on wheat yield potential and prospects of genetic control. Annu. Rev. Phytopathol. 2016, 54, 303–322. [Google Scholar] [CrossRef] [PubMed]
  55. Roy, C.; He, X.; Gahtyari, N.C.; Mahapatra, S.; Singh, P.K. Managing spot blotch disease in wheat: Conventional to molecular aspects. Front Plant Sci. 2023, 14, 1098648. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Phenotypic reactions of the tested genotypes. (A) Boxplot distribution of SNB scores for the 174 wheat genotypes in individual experiments. (B) Histogram distribution of mean disease scores over the three experiments.
Figure 1. Phenotypic reactions of the tested genotypes. (A) Boxplot distribution of SNB scores for the 174 wheat genotypes in individual experiments. (B) Histogram distribution of mean disease scores over the three experiments.
Genes 15 00890 g001
Figure 2. The “which-won-where” view (A) and GGE biplot (B) based on the G × E interaction of SNB disease scores in three different experiments.
Figure 2. The “which-won-where” view (A) and GGE biplot (B) based on the G × E interaction of SNB disease scores in three different experiments.
Genes 15 00890 g002
Figure 3. Population structure (A,B) and PCoA plot (C) of the 174 genotypes.
Figure 3. Population structure (A,B) and PCoA plot (C) of the 174 genotypes.
Genes 15 00890 g003
Figure 4. Heatmap and dendrogram of kinship matrix estimated using Van Randen algorithm based on 11,184 SNP markers and 174 wheat genotypes. Yellow colors represent the haploblocks.
Figure 4. Heatmap and dendrogram of kinship matrix estimated using Van Randen algorithm based on 11,184 SNP markers and 174 wheat genotypes. Yellow colors represent the haploblocks.
Genes 15 00890 g004
Figure 5. Manhattan plots for MTAs detected from individual experiments such as E1 (A), E2 (B), E3 (C) and the pooled data (D). The red horizontal line indicates the significance threshold (p = 0.001) used in the present study.
Figure 5. Manhattan plots for MTAs detected from individual experiments such as E1 (A), E2 (B), E3 (C) and the pooled data (D). The red horizontal line indicates the significance threshold (p = 0.001) used in the present study.
Genes 15 00890 g005
Table 1. Analysis of variance (ANOVA) for SNB resistance in the GWAS panel.
Table 1. Analysis of variance (ANOVA) for SNB resistance in the GWAS panel.
SourceDfMSp-Value
Experiment218.38816p < 0.001
Replication60.128591p > 0.05
Genotype1733.042647p < 0.001
Genotype × Experiment3460.553166p < 0.001
Residual10380.092841
Table 2. Markers significantly associated with seedling SNB resistance using the MLM algorithm.
Table 2. Markers significantly associated with seedling SNB resistance using the MLM algorithm.
Exp.MarkerChrPositionp-ValueR2 (%)
Exp2AX-946776405B6,838,1682 × 10−48.96
MeanAX-946776405B6,838,1683 × 10−48.48
Exp2AX-948645775B6,974,7541 × 10−49.68
MeanAX-948645775B6,974,7543 × 10−48.63
MeanAX-949789393D52,806,8486 × 10−47.65
Exp2AX-949789393D52,806,8487 × 10−47.35
Exp2BS00091519_515B6,648,5172 × 10−48.80
MeanBS00091519_515B6,648,5173 × 10−48.44
Exp2Excalibur_c47452_1835B6,654,1163 × 10−48.22
MeanExcalibur_c47452_1835B6,654,1165 × 10−47.83
Exp2GENE-3324_3385B3,521,2255 × 10−47.76
MeanGENE-3324_3385B3,521,2255 × 10−47.40
Exp2IACX92615B546,703,9366 × 10−413.30
MeanIACX92615B546,703,9366 × 10−410.31
Exp2Tdurum_contig25513_1235B565,753,5032 × 10−615.23
MeanTdurum_contig25513_1235B565,753,5031 × 10−512.93
Exp1Tdurum_contig25513_1235B565,753,5033 × 10−48.59
Exp3Tdurum_contig31131_1985B417,443,5451 × 10−411.34
MeanTdurum_contig31131_1985B417,443,5453 × 10−49.49
Exp2Tdurum_contig31131_1985B417,443,5459 × 10−47.86
Exp2tplb0027f13_14935B565,754,2151 × 10−615.67
Meantplb0027f13_14935B565,754,2152 × 10−511.77
Exp1tplb0027f13_14935B565,754,2152 × 10−47.68
Exp3wsnp_Ex_c210_4116045A584,054,4618 × 10−47.50
Exp3wsnp_Ex_c29304_383554345A584,451,5111 × 10−37.58
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mandal, R.; He, X.; Singh, G.; Kabir, M.R.; Joshi, A.K.; Singh, P.K. Screening of CIMMYT and South Asian Bread Wheat Germplasm Reveals Marker–Trait Associations for Seedling Resistance to Septoria Nodorum Blotch. Genes 2024, 15, 890. https://doi.org/10.3390/genes15070890

AMA Style

Mandal R, He X, Singh G, Kabir MR, Joshi AK, Singh PK. Screening of CIMMYT and South Asian Bread Wheat Germplasm Reveals Marker–Trait Associations for Seedling Resistance to Septoria Nodorum Blotch. Genes. 2024; 15(7):890. https://doi.org/10.3390/genes15070890

Chicago/Turabian Style

Mandal, Rupsanatan, Xinyao He, Gyanendra Singh, Muhammad Rezaul Kabir, Arun Kumar Joshi, and Pawan Kumar Singh. 2024. "Screening of CIMMYT and South Asian Bread Wheat Germplasm Reveals Marker–Trait Associations for Seedling Resistance to Septoria Nodorum Blotch" Genes 15, no. 7: 890. https://doi.org/10.3390/genes15070890

APA Style

Mandal, R., He, X., Singh, G., Kabir, M. R., Joshi, A. K., & Singh, P. K. (2024). Screening of CIMMYT and South Asian Bread Wheat Germplasm Reveals Marker–Trait Associations for Seedling Resistance to Septoria Nodorum Blotch. Genes, 15(7), 890. https://doi.org/10.3390/genes15070890

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop