Since Iran is located in the Fertile Crescent where wheat was initially domesticated more than 12,000 years ago (Blixt, 2009), it is suggested that wheats and their pathogens are in a coevolutionary arms race (Stukenbrock et al., 2007). Our understanding of the interactions between wheat and P. nodorum in natural populations derived from Iran remains limited to date. Hence, the objective of this study was to assess the virulence profile of 33 isolates of P. nodorum on a diverse panel of wheat genotypes, and to analyze the molecular patterns of NEs and S genes in P. nodorum and wheat genotypes, respectively. In a previous study, the genetic structure of P. nodorum was assessed using microsatellite markers, revealing that Iranian isolates of P. nodorum exhibited the highest genetic diversity among 693 isolates collected from five continents worldwide. These findings provided supporting evidence for the hypothesis that P. nodorum has undergone co-evolution with its host in the Fertile Crescent (Ghaderi et al., 2020). This is in line with our current findings, as we have demonstrated that most Iranian isolates exhibited higher DS values compared to the reference isolate SN15. Fars (Pn-FN) isolates had the highest DS value, while Golestan (Pn-GAl) isolates had the lowest DS value compared to other isolates derived from other provinces. Previus study also showed that Fars and Golestan provinces had the highest and lowest genetic diversity, respectively among the Iranian population of P. nodorum (Ghaderi et al., 2020). In addition to their low genetic diversity and virulence, the Golestan isolates exhibited distinct differentiation from other populations. This can be attributed to the geographical separation of Golestan, which is over 800 km distant from the other three collection sites. In addition, physical barriers, such as mountains and deserts, act as differentiating factors between this population and the other populations examined in this study. In contrast to the Golestan population, there was no significant difference in the mean DS value between the Fars and Kohgiluyeh Boyer-Ahmad populations. This lack of difference could potentially be attributed to the presence of gene flow between these two populations. The findings of a previous population structure study also supported the presence of gene flow between the Fars and Kohgiluyeh Boyer-Ahmad populations (Ghaderi et al., 2020).
Our results may also suggest that high genetic diversity is positively correlated with increased virulence. It is reported that a diverse gene pool accelerates the generation of appropriate recombinants compensating for the lack of virulence (Phan et al., 2020). With higher genetic diversity there is a higher chance of selection for highly virulent isolates. In contrast, populations with low genetic diversity tend to be vulnerable to the introduction of new cultivars as they are unable to quickly recombine (Phan et al., 2020). Active sexual reproduction during the growing season enables P. nodorum to evolve quickly leading to the appearance of new virulent pathotypes occupying new ecological niches. Hence, the previously observed high genetic diversity in Fars province (Ghaderi et al., 2020), along with the increased virulence of isolates found in this study, suggests that the Fars population likely engages in more active sexual reproduction compared to the Golestan isolates. However, further comprehensive investigations are required in the future to address this aspect in greater detail.
Pathogenicity assay showed that all Iranian wheat cultivars fall into highly susceptible (47.37%), susceptible (31.58%) and moderately susceptible (21.06%) cultivars to SNB. SnTox1-Snn1 was the most abundant NE-S gene frequency (100%) in Iranian P. nodorum isolates and Iranian wheat cultivars. Contrary to the most NEs, the C-terminal chitin-binding region of NE SnTox1 acts as a PAMP-type stimulator and interacts with the cell surface pattern recognition receptor Snn1 (Liu et al., 2012) triggering PTI (rather than ETI) response in the host, which is one of the general defense mechanisms (Liu et al., 2016; Faris & Friesen, 2020). Therefore, considering the dual function of SnTox1 both as a necrotrophic effector stimulating the PTI in the plant and its role in protecting the fungus against wheat chitinases (Liu et al., 2016), the high-frequency of SnTox1 in the Iranian P. nodorum population is logical, which is comparable with its relatively high frequency (85%) in a large global P. nodorum populations as determined previously (Liu et al., 2012; Faris & Friesen, 2020). Considering the protective role of SnTox1 in disease development, which operates independently of its function in inducing cell death (Liu et al., 2016; McDonald et al., 2022), it was expected to observe a high frequency (100%) of SnTox1 presence across all disease severity (DS) categories (Table 6).
Table 6
The proportion of SNB susceptibility (S) and necrotrophic effector (NE) genes under different disease severity (DS) classification
S /NE gene | DS > 41 | 26 ≤ DS ≥ 40 | 16 ≤ DS ≥ 25 | 11 ≤ DS ≥ 15 | 6 ≤ DS ≥ 10 | 1 ≤ DS ≥ 5 |
Tsn1 | 70 | 50 | 50 | 0 | 0 | 0 |
Snn1 | 90 | 80 | 75 | 100 | 100 | 0 |
Snn3 | 100 | 70 | 75 | 100 | 57 | 0 |
SnToxA | 100 | 69 | 50 | 50 | - | - |
SnTox1 | 100 | 100 | 100 | 100 | - | - |
SnTox3 | 100 | 69 | 67 | 75 | - | - |
SnTox5 | 100 | 94 | 100 | 100 | - | - |
S gene SnToxA is the main determinant in causing SNB, tan spot, and spot blotch diseases in plant lines harboring Tsn1. In the absence of Tsn1, the severity of the disease reduces significantly even in the presence of other S genes (Peters Haugrud et al., 2022; McCombe et al., 2022; Liu et al., 2012) which is consistent with the results here. S gene Tsn1 was present in 70% of highly susceptible cultivars. Tsn1 frequency decreased to 50% of susceptible and moderately susceptible cultivars and was completely absent in moderately resistant, resistant and highly resistant cultivars (Table 6).
The mean DS (36.87%) of Iranian wheat cultivars (regarded as susceptible) was correlated with the presence of all S genes Tsn1, Snn1 and Snn3 (84.21%). Additionally, the mean disease severity (DS) of world cultivars (14.81%), categorized as moderately resistant, corresponded to the low proportion (58.33%) of having all three S genes. Those cultivars possessing all three S genes Tsn1, Snn1 and Snn3 were all categorized as highly susceptible, susceptible and moderately susceptible. Also, the frequency of S genes in highly susceptible cultivars (86.66%) was much higher than the frequency of these genes in resistant cultivars (52%). Furthermore, an increased number of NE gene SnToxA was positively correlated to virulence. In DS > 41, the frequency of SnToxA was 100% decreasing gradually to 50% of P. nodorum samples in 16 ≤ DS ≥ 25 and 11 ≤ DS ≥ 15 (Table 6). Despite the absence of certain NEs in some isolates of P. nodorum, they exhibited high virulence, indicating the potential existence of other unidentified NEs and additional disease determinants. This finding further confirms the complex nature of the P. nodorum-wheat interaction (Oliver et al., 2009). It was suggested that susceptible cultivars carrying multiple S genes could rapidly select pathogen isolates having the corresponding NE genes (McDonald et al., 2013). Therefore, the frequency of the NEs in the pathogen population could be directly related to the frequency of the S genes in its host. In the absence of host S genes, the NE will be neutral and subjected to genetic drift (McDonald et al., 2013). Although the overall higher occurrence of both NEs and S genes in Iranian P. nodorum isolates and wheat genotypes supports this hypothesis, drawing a definitive conclusion requires monitoring a large population size of P. nodorum and host genotypes over multiple years.
Wheat cultivars were categorized into six groups (highly resistant, resistant, moderately resistant, moderately susceptible and highly susceptible) based on their variable reaction towards P. nodorum isolates. Most previous studies rated wheat plants for their reaction to SNB by a qualitative 0–5 scale (Liu et al., 2004; Károlyiné Cséplő et al., 2013; Singh et al., 2006; Feng et al., 2004). The quantitative scale implemented in this study can be compared to the previous qualitative scale, where values 0, 1, 2, 3, 4, and 5 correspond to immune, highly resistant, resistant, moderately susceptible, susceptible, and highly susceptible, respectively (Singh et al., 2006; Feng et al., 2004). Therefore, based on the quantitative scale defined here, Iranian wheat cultivars with a mean DS value of 36.87% were categorized as susceptible while commercial world cultivars with a mean DS value of 14.81% were classified as moderately resistant cultivars to SNB. The findings regarding the response of differential and world commercial cultivars to SNB were in line with qualitative assays conducted by previous researchers where BR34 was classified as highly resistant, Erik, M-3, Vernapolis, and Arina as resistant, ND 495 and Katepwa as moderately susceptible, 4B-160 and Coulter as susceptible, and Grandin as highly susceptible cultivar to SNB (Liu et al., 2004; Singh et al., 2006; Károlyiné Cséplő et al., 2013; Tommasini et al., 2007). Hence, the utilization of the quantitative scale as a substitutional assay for assessing the response of wheat cultivars to P. nodorum provides a more precise and reliable approach, which can serve as a valuable parameter in breeding programs.
NE infiltration represents one of the simplest and quickest strategies for evaluating the sensitivity reaction of host plants. Response to toxins occurs quickly, which provides accurate and reproducible results by minimizing the environmental effects on symptom development facilitating the early detection of resistant and susceptible cultivars (Sharma, 2012; Waters et al., 2011). As an example, the sensitivity of 60 bread wheat cultivars from Western Australian was assessed through NE SnTox3 infiltration, which served as an expedited approach for breeding purposes (Waters et al., 2011). The screening of wheat cultivars using purified toxins allows breeders to promptly identify and eliminate susceptible genotypes in the early stages of breeding cycles (Oliver et al., 2009). As mentioned earlier, the interaction between SnToxA and Tsn1 is an important factor in the development of SNB. Therefore, in order to achieve a high level of resistance against both SNB and tan spot, breeders should prioritize the selection of genotypes that lack the Tsn1 gene (Peters Haugrud et al., 2023). Considering the coexistence of SNB and tan spot, two major diseases in wheat (Phan et al., 2020), the application of NE SnToxA holds potential for marker-assisted selection in screening wheat cultivars against SNB, tan spot, and spot blotch diseases. The developed SnToxA cloning and infiltration assay demonstrated the capability to identify Tsn1-resistant cultivars within a 72-hour timeframe, enabling early detection and selection.
Genetic analysis of wheat lines revealed that the resistance to SNB can be inherited and controlled in a quantitative trait (Liu et al., 2004; Friesen et al., 2008). Most quantitative trait loci have minor effects posing a challenge for their utilization in breeding programs (Solomon et al., 2006; Liu et al., 2004). Hence, adopting new screening strategies developed based on the mechanisms of interaction between P. nodorum and wheat is essential. These methods, including NE infiltration assays, glasshouse assays, and PCR-based disease susceptibility screening (Waters et al., 2011), effectively contribute to the selection of non-susceptible wheat cultivars, thereby enhancing our capability to improve wheat breeding programs. Identifying the presence of NE genes alerts breeders to avoid using corresponding sensitivity genes in the host. In this study, we have demonstrated the effectiveness of NE infiltration coupled with multiplex PCR as a rapid and efficient method for detecting the S gene in the host and NE genes in P. nodorum. This offers valuable directions and tools that can significantly enhance breeding programs aimed at developing wheat cultivars with increased resistance against SNB.