GWAS analysis reveals distinct pathogenicity profiles of Australian Parastagonospora nodorum isolates and identification of marker-trait-associations to septoria nodorum blotch.

The fungus Parastagonospora nodorum is the causal agent of septoria nodorum leaf blotch (SNB) and glume blotch which are common in many wheat growing regions in the world. The disease is complex and could be explained by multiple interactions between necrotrophic effectors secreted by the pathogen and matching susceptibility genes in wheat. An Australian P. nodorum population was clustered into five groups with contrasting properties. This study was set to identify their pathogenicity profiles using a diverse wheat panel of 134 accessions which are insensitive to SnToxA and SnTox1 in both in vitro and in vivo conditions. SNB seedling resistance/susceptibility to five representative isolates from the five clusters, responses to crude culture-filtrates (CFs) of three isolates and sensitivity to SnTox3 semi-purified effector together with 11,455 SNP markers have been used for linkage disequilibrium (LD) and association analyses. While quantitative trait loci (QTL) on 1D, 2A, 2B, 4B, 5B, 6A, 6B, 7A, 7D chromosomes were consistently detected across isolates and conditions, distinct patterns and isolate specific QTL were also observed among these isolates. In this study, SnTox3-Snn3-B1 interaction for the first time in Australia and SnTox3-Snn3-D1 interaction for the first time in bread wheat were found active using wild-type isolates. These findings could be due to new SnTox3 haplotype/isoform and exotic CIMMYT/ICARDA and Vavilov germplasm used, respectively. This study could provide useful information for dissecting novel and different SNB disease components, helping to prioritise research targets and contributing valuable information on genetic loci/markers for marker-assisted selection in SNB resistance wheat breeding programme.

A specific fungal transcription factor controls effector gene expression and orchestrates the establishment of the necrotrophic pathogen lifestyle on wheat.

The fungus Parastagonospora nodorum infects wheat through the use of necrotrophic effector (NE) proteins that cause host-specific tissue necrosis. The Zn2Cys6 transcription factor PnPf2 positively regulates NE gene expression and is required for virulence on wheat. Little is known about other downstream targets of PnPf2. We compared the transcriptomes of the P. nodorum wildtype and a strain deleted in PnPf2 (pf2-69) during in vitro growth and host infection to further elucidate targets of PnPf2 signalling. Gene ontology enrichment analysis of the differentially expressed (DE) genes revealed that genes associated with plant cell wall degradation and proteolysis were enriched in down-regulated DE gene sets in pf2-69 compared to SN15. In contrast, genes associated with redox control, nutrient and ion transport were up-regulated in the mutant. Further analysis of the DE gene set revealed that PnPf2 positively regulates twelve genes that encode effector-like proteins. Two of these genes encode proteins with homology to previously characterised effectors in other fungal phytopathogens. In addition to modulating effector gene expression, PnPf2 may play a broader role in the establishment of a necrotrophic lifestyle by orchestrating the expression of genes associated with plant cell wall degradation and nutrient assimilation.

Low Amplitude Boom-and-Bust Cycles Define the Septoria Nodorum Blotch Interaction.

INTRODUCTION: Septoria nodorum blotch (SNB) is a complex fungal disease of wheat caused by the Dothideomycete fungal pathogen Parastagonospora nodorum. The fungus infects through the use of necrotrophic effectors (NEs) that cause necrosis on hosts carrying matching dominant susceptibility genes. The Western Australia (WA) wheatbelt is a SNB "hot spot" and experiences significant under favorable conditions. Consequently, SNB has been a major target for breeders in WA for many years. MATERIALS AND METHODS: In this study, we assembled a panel of 155 WA P. nodorum isolates collected over a 44-year period and compared them to 23 isolates from France and the USA using 28 SSR loci. RESULTS: The WA P. nodorum population was clustered into five groups with contrasting properties. 80% of the studied isolates were assigned to two core groups found throughout the collection location and time. The other three non-core groups that encompassed transient and emergent populations were found in restricted locations and time. Changes in group genotypes occurred during periods that coincided with the mass adoption of a single or a small group of widely planted wheat cultivars. When introduced, these cultivars had high scores for SNB resistance. However, the field resistance of these new cultivars often declined over subsequent seasons prompting their replacement with new, more resistant varieties. Pathogenicity assays showed that newly emerged isolates non-core are more pathogenic than old isolates. It is likely that the non-core groups were repeatedly selected for increased virulence on the contemporary popular cultivars. DISCUSSION: The low level of genetic diversity within the non-core groups, difference in virulence, low abundance, and restriction to limited locations suggest that these populations more vulnerable to a population crash when the cultivar was replaced by one that was genetically different and more resistant. We characterize the observed pattern as a low-amplitude boom-and-bust cycle in contrast with the classical high amplitude boom-and-bust cycles seen for biotrophic pathogens where the contrast between resistance and susceptibility is typically much greater. Implications of the results are discussed relating to breeding strategies for more sustainable SNB resistance and more generally for pathogens with NEs.

Pan-Parastagonospora Comparative Genome Analysis-Effector Prediction and Genome Evolution.

We report a fungal pan-genome study involving Parastagonospora spp., including 21 isolates of the wheat (Triticum aestivum) pathogen Parastagonospora nodorum, 10 of the grass-infecting Parastagonospora avenae, and 2 of a closely related undefined sister species. We observed substantial variation in the distribution of polymorphisms across the pan-genome, including repeat-induced point mutations, diversifying selection and gene gains and losses. We also discovered chromosome-scale inter and intraspecific presence/absence variation of some sequences, suggesting the occurrence of one or more accessory chromosomes or regions that may play a role in host-pathogen interactions. The presence of known pathogenicity effector loci SnToxA, SnTox1, and SnTox3 varied substantially among isolates. Three P. nodorum isolates lacked functional versions for all three loci, whereas three P. avenae isolates carried one or both of the SnTox1 and SnTox3 genes, indicating previously unrecognized potential for discovering additional effectors in the P. nodorum-wheat pathosystem. We utilized the pan-genomic comparative analysis to improve the prediction of pathogenicity effector candidates, recovering the three confirmed effectors among our top-ranked candidates. We propose applying this pan-genomic approach to identify the effector repertoire involved in other host-microbe interactions involving necrotrophic pathogens in the Pezizomycotina.

Novel sources of resistance to Septoria nodorum blotch in the Vavilov wheat collection identified by genome-wide association studies.

KEY MESSAGE: The fungus Parastagonospora nodorum causes Septoria nodorum blotch (SNB) of wheat. A genetically diverse wheat panel was used to dissect the complexity of SNB and identify novel sources of resistance. The fungus Parastagonospora nodorum is the causal agent of Septoria nodorum blotch (SNB) of wheat. The pathosystem is mediated by multiple fungal necrotrophic effector-host sensitivity gene interactions that include SnToxA-Tsn1, SnTox1-Snn1, and SnTox3-Snn3. A P. nodorum strain lacking SnToxA, SnTox1, and SnTox3 (toxa13) retained wild-type-like ability to infect some modern wheat cultivars, suggesting evidence of other effector-mediated susceptibility gene interactions or the lack of host resistance genes. To identify genomic regions harbouring such loci, we examined a panel of 295 historic wheat accessions from the N. I. Vavilov Institute of Plant Genetic Resources in Russia, which is comprised of genetically diverse landraces and breeding lines registered from 1920 to 1990. The wheat panel was subjected to effector bioassays, infection with P. nodorum wild type (SN15) and toxa13. In general, SN15 was more virulent than toxa13. Insensitivity to all three effectors contributed significantly to resistance against SN15, but not toxa13. Genome-wide association studies using phenotypes from SN15 infection detected quantitative trait loci (QTL) on chromosomes 1BS (Snn1), 2DS, 5AS, 5BS (Snn3), 3AL, 4AL, 4BS, and 7AS. For toxa13 infection, a QTL was detected on 5AS (similar to SN15), plus two additional QTL on 2DL and 7DL. Analysis of resistance phenotypes indicated that plant breeders may have inadvertently selected for effector insensitivity from 1940 onwards. We identify accessions that can be used to develop bi-parental mapping populations to characterise resistance-associated alleles for subsequent introgression into modern bread wheat to minimise the impact of SNB.

A functionally conserved Zn2 Cys6 binuclear cluster transcription factor class regulates necrotrophic effector gene expression and host-specific virulence of two major Pleosporales fungal pathogens of wheat.

The fungus Parastagonospora nodorum is the causal agent of Septoria nodorum blotch of wheat (Triticum aestivum). The interaction is mediated by multiple fungal necrotrophic effector-dominant host sensitivity gene interactions. The three best-characterized effector-sensitivity gene systems are SnToxA-Tsn1, SnTox1-Snn1 and SnTox3-Snn3. These effector genes are highly expressed during early infection, but expression decreases as the infection progresses to tissue necrosis and sporulation. However, the mechanism of regulation is unknown. We have identified and functionally characterized a gene, referred to as PnPf2, which encodes a putative zinc finger transcription factor. PnPf2 deletion resulted in the down-regulation of SnToxA and SnTox3 expression. Virulence on Tsn1 and Snn3 wheat cultivars was strongly reduced. The SnTox1-Snn1 interaction remained unaffected. Furthermore, we have also identified and deleted an orthologous PtrPf2 from the tan spot fungus Pyrenophora tritici-repentis which possesses a near-identical ToxA that was acquired from P. nodorum via horizontal gene transfer. PtrPf2 deletion also resulted in the down-regulation of PtrToxA expression and a near-complete loss of virulence on Tsn1 wheat. We have demonstrated, for the first time, evidence for a functionally conserved signalling component that plays a role in the regulation of a common/horizontally transferred effector found in two major fungal pathogens of wheat.

Differential effector gene expression underpins epistasis in a plant fungal disease.

Fungal effector-host sensitivity gene interactions play a key role in determining the outcome of septoria nodorum blotch disease (SNB) caused by Parastagonospora nodorum on wheat. The pathosystem is complex and mediated by interaction of multiple fungal necrotrophic effector-host sensitivity gene systems. Three effector sensitivity gene systems are well characterized in this pathosystem; SnToxA-Tsn1, SnTox1-Snn1 and SnTox3-Snn3. We tested a wheat mapping population that segregated for Snn1 and Snn3 with SN15, an aggressive P. nodorum isolate that produces SnToxA, SnTox1 and SnTox3, to study the inheritance of sensitivity to SnTox1 and SnTox3 and disease susceptibility. Interval quantitative trait locus (QTL) mapping showed that the SnTox1-Snn1 interaction was paramount in SNB development on both seedlings and adult plants. No effect of the SnTox3-Snn3 interaction was observed under SN15 infection. The SnTox3-Snn3 interaction was however, detected in a strain of SN15 in which SnTox1 had been deleted (tox1-6). Gene expression analysis indicates increased SnTox3 expression in tox1-6 compared with SN15. This indicates that the failure to detect the SnTox3-Snn3 interaction in SN15 is due - at least in part - to suppressed expression of SnTox3 mediated by SnTox1. Furthermore, infection of the mapping population with a strain deleted in SnToxA, SnTox1 and SnTox3 (toxa13) unmasked a significant SNB QTL on 2DS where the SnTox2 effector sensitivity gene, Snn2, is located. This QTL was not observed in SN15 and tox1-6 infections and thus suggesting that SnToxA and/or SnTox3 were epistatic. Additional QTLs responding to SNB and effectors sensitivity were detected on 2AS1 and 3AL.

Dissecting the role of histidine kinase and HOG1 mitogen-activated protein kinase signalling in stress tolerance and pathogenicity of Parastagonospora nodorum on wheat.

The HOG1 mitogen-activated protein kinase (MAPK) pathway is activated through two-component histidine kinase (HK) signalling. This pathway was first characterized in the budding yeast Saccharomyces cerevisiae as a regulator of osmotolerance. The fungus Parastagonospora nodorum is the causal agent of septoria nodorum blotch of wheat. This pathogen uses host-specific effectors in tandem with general pathogenicity mechanisms to carry out its infection process. Genes showing strong sequence homology to S. cerevisiae HOG1 signalling pathway genes have been identified in the genome of P. nodorum. In this study, we examined the role of the pathway in the virulence of P. nodorum on wheat by disrupting putative pathway component genes: HOG1 (SNOG_13296) MAPK and NIK1 (SNOG_11631) hybrid HK. Mutants deleted in NIK1 and HOG1 were insensitive to dicarboximide and phenylpyrrole fungicides, but not a fungicide that targets ergosterol biosynthesis. Furthermore, both Δnik1 and Δhog1 mutants showed increased sensitivity to hyperosmotic stress. However, HOG1, but not NIK1, is required for tolerance to elevated temperatures. HOG1 deletion conferred increased tolerance to 6-methoxy-2-benzoxazolinone, a cereal phytoalexin. This suggests that the HOG1 signalling pathway is not exclusively associated with NIK1. Both Δnik1 and Δhog1 mutants retained the ability to infect and cause necrotic lesions on wheat. However, we observed that the Δhog1 mutation resulted in reduced production of pycnidia, asexual fruiting bodies that facilitate spore dispersal during late infection. Our study demonstrated the overlapping and distinct roles of a HOG1 MAPK and two-component HK signalling in P. nodorum growth and pathogenicity.

Functional redundancy of necrotrophic effectors - consequences for exploitation for breeding.

Necrotrophic diseases of wheat cause major losses in most wheat growing areas of world. Tan spot (caused by Pyrenophora tritici-repentis) and septoria nodorum blotch (SNB; Parastagonospora nodorum) have been shown to reduce yields by 10-20% across entire agri-ecological zones despite the application of fungicides and a heavy focus over the last 30 years on resistance breeding. Efforts by breeders to improve the resistance of cultivars has been compromised by the universal finding that resistance was quantitative and governed by multiple quantitative trait loci (QTL). Most QTL had a limited effect that was hard to measure precisely and varied significantly from site to site and season to season. The discovery of necrotrophic effectors has given breeding for disease resistance new methods and tools. In the case of tan spot in West Australia, a single effector, PtrToxA and its recogniser gene Tsn1, has a dominating impact in disease resistance. The delivery of ToxA to breeders has had a major impact on cultivar choice and breeding strategies. For P. nodorum, three effectors - SnToxA, SnTox1, and SnTox3 - have been well characterized. Unlike tan spot, no one effector has a dominating role. Genetic analysis of various mapping populations and pathogen isolates has shown that different effectors have varying impact and that epistatic interactions also occur. As a result of these factors the deployment of these effectors for SNB resistance breeding is more complex. We have deleted the three effectors in a strain of P. nodorum and measured effector activity and disease potential of the triple knockout mutant. The culture filtrate causes necrosis in several cultivars and the strain causes disease, albeit the overall levels are less than in the wild type. Modeling of the field disease resistance scores of cultivars from their reactions to the microbially expressed effectors SnToxA, SnTox1, and SnTox3 is significantly improved by including the response to the triple knockout mutant culture filtrate. This indicates that one or more further effectors are secreted into the culture filtrate. We conclude that the in vitro-secreted necrotrophic effectors explain a very large part of the disease response of wheat germplasm and that this method of resistance breeding promises to further reduce the impact of these globally significant diseases.

Quantitative variation in effector activity of ToxA isoforms from Stagonospora nodorum and Pyrenophora tritici-repentis.

ToxA is a proteinaceous necrotrophic effector produced by Stagonospora nodorum and Pyrenophora tritici-repentis. In this study, all eight mature isoforms of the ToxA protein were purified and compared. Circular dichroism spectra indicated that all isoforms were structurally intact and had indistinguishable secondary structural features. ToxA isoforms were infiltrated into wheat lines that carry the sensitivity gene Tsn1. It was observed that different wheat lines carrying identical Tsn1 alleles varied in sensitivity to ToxA. All ToxA isoforms induced necrosis when introduced into any Tsn1 wheat line but we observed quantitative variation in effector activity, with the least-active version found in isolates of P. tritici-repentis. Pathogen sporulation increased with higher doses of ToxA. The isoforms that induced the most rapid necrosis also induced the most sporulation, indicating that pathogen fitness is affected by differences in ToxA activity. We show that differences in toxin activity encoded by a single gene can contribute to the quantitative inheritance of necrotrophic virulence. Our findings support the hypothesis that the variation at ToxA results from selection that favors increased toxin activity.

Trehalose biosynthesis is involved in sporulation of Stagonospora nodorum.

Stagonospora nodorum is a necrotrophic fungal pathogen that is the causal agent of leaf and glume blotch on wheat. S. nodorum is a polycyclic pathogen, whereby rain-splashed pycnidiospores attach to and colonise wheat tissue and subsequently sporulate again within 2-3weeks. As several cycles of infection are needed for a damaging infection, asexual sporulation is a critical phase of its infection cycle. A non-targeted metabolomics screen for sporulation-associated metabolites identified that trehalose accumulated significantly in concert with asexual sporulation both in vitro and in planta. A reverse-genetics approach was used to investigate the role of trehalose in asexual sporulation. Trehalose biosynthesis was disrupted by deletion of the gene Tps1, encoding a trehalose 6-phosphate synthase, resulting in almost total loss of trehalose during in vitro growth and in planta. In addition, lesion development and pycnidia formation were also significantly reduced in tps1 mutants. Reintroduction of the Tps1 gene restored trehalose biosynthesis, pathogenicity and sporulation to wild-type levels. Microscopic examination of tps1 infected wheat leaves showed that pycnidial formation often halted at an early stage of development. Further examination of the tps1 phenotype revealed that tps1 pycnidiospores exhibited a reduced germination rate while under heat stress, and tps1 mutants had a reduced growth rate while under oxidative stress. This study confirms a link between trehalose biosynthesis and pathogen fitness in S.nodorum.

A metabolomic approach to dissecting osmotic stress in the wheat pathogen Stagonospora nodorum.

A non-targeted metabolomics approach was used to identify significant changes in metabolism upon exposure of the wheat pathogen Stagonospora nodorum to 0.5M NaCl. The polyol arabitol, and to a lesser extent glycerol, was found to accumulate in response to the osmotic stress treatment. Amino acid synthesis was strongly down-regulated whilst mannitol levels were unaffected. A reverse genetic approach was undertaken to dissect the role of arabitol metabolism during salt stress. Strains of S. nodorum lacking a gene encoding an l-arabitol dehydrogenase (abd1), a xylitol dehydrogenase (xdh1) and a double-mutant lacking both genes (abd1xdh1) were exposed to salt and the intracellular metabolites analysed. Arabitol levels were significantly up-regulated upon salt stress in the xdh1 strains but were significantly lower than the wild-type. Arabitol was not significantly different in either the abd1 or the abd1xdh1 strains during osmotic stress but the concentration of glycerol was significantly higher indicating a compensatory mechanism in operation. Genome sequence analysis identified a second possible enzyme capable of synthesizing arabitol explaining the basal level of arabitol present in the abd1xdh1 strains. This study identified that arabitol is the primary compatible solute in S. nodorum but in-built levels of redundancy are present allowing the fungus to tolerate osmotic stress.

Investigating the role of calcium/calmodulin-dependent protein kinases in Stagonospora nodorum.

Three genes encoding different Ca2+/calmodulin-dependent protein kinases have been characterized in the wheat phytopathogenic fungus Stagonospora nodorum. The kinases were identified from the S. nodorum genome sequence on the basis of sequence homology to known Ca2+/calmodulin-dependent protein kinases. Expression analysis determined that each of the kinases was expressed during growth in vitro and also during infection. The onset of sporulation triggered increased transcript levels of each of the kinases, particularly CpkA where an 11-fold increase in expression was observed during sporulation in planta. The role of the kinases was further determined via a reverse genetics approach. The disruption of CpkA affected vegetative growth in vitro and also sporulation. The cpkA strains produced 20-fold less spores on complex media and were unable to sporulate on defined minimal media. Infection assays showed that CpkA was not required for lesion development but was essential for sporulation at the completion of the infection cycle. Microscopic analysis revealed that the disruption of CpkA resulted in Stagonospora nodorum being unable to differentiate the mycelial knot into immature pycnidia during sporulation. A metabolite analysis of infected leaves during sporulation excluded the possible involvement of mannitol, a compound previously shown to be involved in the sporulation of Stagonospora nodorum. The disruption of CpkB did not effect growth in vitro or pathogenicity. Stagonospora nodorum strains lacking CpkC appeared unaffected during growth in planta but showed delayed lesion development and sporulation during infection.