Global diversity and distribution of three necrotrophic effectors in Phaeosphaeria nodorum and related species

Discovery of a New Starship Transposon Driving the Horizontal Transfer of the ToxA Virulence Gene in Alternaria ventricosa

The virulence gene ToxA has been proposed to be horizontally transferred between three fungal wheat pathogens (Parastagonospora nodorum, Pyrenophora tritici-repentis, and Bipolaris sorokiniana) as part of a conserved ~14 kb ToxhAT transposon. Here, our analysis of 2137 fungal species-representative assemblies revealed that the ToxA gene is an isolate of Alternaria ventricosa and shows a remarkable 99.5% similarity to those found in B. sorokiniana and P. tritici-repentis. Analysis of the regions flanking ToxA within A. ventricosa revealed that it was embedded within a 14 kb genomic element nearly identical to the corresponding ToxhAT regions in B. sorokiniana, P. nodorum, and P. tritici-repentis. Comparative analysis further showed that ToxhAT in A. ventricosa resides within a larger mobile genetic element, which we identified as a member of the Starship transposon superfamily, named Frontier. Our analysis demonstrated that ToxhAT has been independently captured by three distinct Starships-Frontier, Sanctuary, and Horizon-which, despite having minimal sequence similarity outside of ToxhAT, facilitate its mobilization. These findings place Frontier, Sanctuary, and Horizon within a growing class of Starships implicated in the horizontal transfer of adaptive genes among fungal species. Moreover, we identified three distinct HGT events involving ToxA across these four fungal species, reinforcing the hypothesis of a single evolutionary origin for the ToxhAT transposon. These findings underscore the pivotal role of transposon-mediated HGT in the adaptive evolution of eukaryotic pathogens, offering new insights into how transposons facilitate genetic exchange and shape host-pathogen interactions in fungi.

Repeat-induced point mutations driving Parastagonospora nodorum genomic diversity are balanced by selection against non-synonymous mutations

Parastagonospora nodorum is necrotrophic fungal pathogen of wheat with significant genomic resources. Population-level pangenome data for 173 isolates, of which 156 were from Western Australia (WA) and 17 were international, were examined for overall genomic diversity and effector gene content. A heterothallic core population occurred across all regions of WA, with asexually-reproducing clonal clusters in dryer northern regions. High potential for SNP diversity in the form of repeat-induced point mutation (RIP)-like transitions, was observed across the genome, suggesting widespread 'RIP-leakage' from transposon-rich repetitive sequences into non-repetitive regions. The strong potential for RIP-like mutations was balanced by negative selection against non-synonymous SNPs, that was observed within protein-coding regions. Protein isoform profiles of known effector loci (SnToxA, SnTox1, SnTox3, SnTox267, and SnTox5) indicated low-levels of non-synonymous and high-levels of silent RIP-like mutations. Effector predictions identified 186 candidate secreted predicted effector proteins (CSEPs), 69 of which had functional annotations and included confirmed effectors. Pangenome-based effector isoform profiles across WA were distinct from global isolates and were conserved relative to population structure, and may enable new approaches for monitoring crop disease pathotypes.

Evolution, diversity, and function of the disease susceptibility gene Snn1 in wheat

Septoria nodorum blotch (SNB), caused by Parastagonospora nodorum, is a disease of durum and common wheat initiated by the recognition of pathogen-produced necrotrophic effectors (NEs) by specific wheat genes. The wheat gene Snn1 was previously cloned, and it encodes a wall-associated kinase that directly interacts with the NE SnTox1 leading to programmed cell death and ultimately the development of SNB. Here, sequence analysis of Snn1 from 114 accessions including diploid, tetraploid, and hexaploid wheat species revealed that some wheat lines possess two copies of Snn1 (designated Snn1-B1 and Snn1-B2) approximately 120 kb apart. Snn1-B2 evolved relatively recently as a paralog of Snn1-B1, and both genes have undergone diversifying selection. Three point mutations associated with the formation of the first SnTox1-sensitive Snn1-B1 allele from a primitive wild wheat were identified. Four subsequent and independent SNPs, three in Snn1-B1 and one in Snn1-B2, converted the sensitive alleles to insensitive forms. Protein modeling indicated these four mutations could abolish Snn1-SnTox1 compatibility either through destabilization of the Snn1 protein or direct disruption of the protein-protein interaction. A high-throughput marker was developed for the absent allele of Snn1, and it was 100% accurate at predicting SnTox1-insensitive lines in both durum and spring wheat. Results of this study increase our understanding of the evolution, diversity, and function of Snn1-B1 and Snn1-B2 genes and will be useful for marker-assisted elimination of these genes for better host resistance.

A New ToxA Haplotype in the Wheat Fungal Pathogen Bipolaris sorokiniana

The necrotrophic effector ToxA is a well-studied virulence factor produced by several fungal necrotrophs. Initially cloned from the wheat tan spot pathogen Pyrenophora tritici-repentis in 1996, ToxA was found almost a decade later in another fungal pathogen, Parastagonospora nodorum, and its sister species, Parastagonospora pseudonodorum. In 2018, ToxA was detected in a third wheat fungal pathogenic species, Bipolaris sorokiniana, which causes spot blotch disease. However, unlike the case with P. tritici-repentis and P. nodorum, the ToxA in B. sorokiniana has only been investigated in recent years. In this report, five Australian B. sorokiniana isolates were assessed for the presence of ToxA. Four isolates were found to contain ToxA. While one isolate harbored the previously reported ToxA haplotype sequence (ToxA19), three isolates contain a different haplotype, designated herein as ToxA25, which has a nonsynonymous mutation resulting in an amino acid change of glycine to arginine at position 168. Both B. sorokiniana ToxA isoforms, when heterologously expressed in Escherichia coli, exhibited the classic ToxA necrosis-inducing activity on ToxA-sensitive Tsn1 cultivars. Preliminary analysis of the B. sorokiniana isolates in Australian wheat cultivars showed that isolates with ToxA19, ToxA25, or ToxA-deficient displayed various degrees of virulence, with the most aggressive isolates observed for those producing ToxA. Differences in spot blotch disease severity between Tsn1 and tsn1 cultivars were observed; however, this was not limited to the ToxA-producing isolates. The overall results suggest that the virulence of the Australian B. sorokiniana isolates is diverse, with the significance of ToxA-Tsn1 interactions depending on individual isolates.

Evolution of the Toxb Gene in Pyrenophora tritici-repentis and Related Species

Pyrenophora tritici-repentis (tan spot) is a destructive foliar pathogen of wheat with global impact. This ascomycete fungus possesses a highly plastic open pangenome shaped by the gain and loss of effector genes. This study investigated the allelic variations in the chlorosis-encoding gene ToxB across 422 isolates representing all identified pathotypes and worldwide origins. To gain better insights into ToxB evolution, we examined its presence and variability in other Pyrenophora spp. A ToxB haplotype network was constructed, revealing the evolutionary relationships of this gene (20 haplotypes) across four Pyrenophora species. Notably, toxb, the homolog of ToxB, was detected for the first time in the barley pathogen Pyrenophora teres. The ToxB/toxb genes display evidence of selection that is characterized by loss of function, duplication, and diverse mutations. Within the ToxB/toxb open reading frame, 72 mutations were identified, including 14 synonymous, 55 nonsynonymous, and 3 indel mutations. Remarkably, a, ∼5.6-kb Copia-like retrotransposon, named Copia-1_Ptr, was found inserted in the toxb gene of a race 3 isolate. This insert disrupted the ToxB gene's function, a first case of effector gene disruption by a transposable element in P. tritici-repentis. Additionally, a microsatellite with 25 nucleotide repeats (0 to 10) in the upstream region of ToxB suggested a potential mechanism influencing ToxB expression and regulation. Exploring ToxB-like protein distribution in other ascomycetes revealed the presence of ToxB-like proteins in 19 additional species, including the Leotiomycetes class for the first time. The presence/absence pattern of ToxB-like proteins defied species relatedness compared with a phylogenetic tree, suggesting a past horizontal gene transfer event during the evolution of the ToxB gene. [Formula: see text] Copyright © 2024 His Majesty the King in Right of Canada, as represented by the Minister of Agriculture and Agri-Food. This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Evaluation of Stagonospora Nodorum Blotch Severity and Parastagonospora nodorum Population Structure and Genetic Diversity Across Multiple Locations and Wheat Varieties in Virginia

Parastagonospora nodorum is a necrotrophic pathogen that causes Stagonospora nodorum blotch (SNB) in wheat. Wheat varieties grown in Virginia vary in susceptibility to SNB, and the severity of SNB varies across locations and years. However, the impacts of wheat genetic backgrounds and environments on SNB severity and the structure of P. nodorum populations in the region have not been well studied. Thus, a population genetic study was conducted utilizing P. nodorum isolates collected from different wheat varieties and locations in Virginia. A total of 320 isolates were collected at seven locations over 2 years from five wheat varieties. Isolates were genotyped using multilocus simple sequence repeat markers, and necrotrophic effector (NE) and mating type genes were amplified using gene-specific primers. Wheat varieties varied in susceptibility to SNB, but site-specific environmental conditions were the primary drivers of disease severity. Fungal populations were genetically diverse, but no genetic subdivision was observed among locations or varieties. The ratio of the two mating type idiomorphs was not significantly different from 1:1, consistent with the P. nodorum population undergoing sexual reproduction. Three major NE genes were detected within the P. nodorum population, but not with equal frequency. However, NE gene profiles were similar for groups of isolates originating from different varieties, suggesting that wheat genetic backgrounds do not differentially select for NEs. There was no evidence of population structure among P. nodorum populations in Virginia and, thus, no support for wheat genetic backgrounds shaping these populations. Finally, although varieties only exhibited moderate resistance to SNB, current levels of resistance are likely to be durable over time and remain a useful tool for integrated management of SNB in the region. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

The Necrotrophic Pathogen Parastagonospora nodorum Is a Master Manipulator of Wheat Defense

Parastagonospora nodorum is a necrotrophic pathogen of wheat that is particularly destructive in major wheat-growing regions of the United States, northern Europe, Australia, and South America. P. nodorum secretes necrotrophic effectors that target wheat susceptibility genes to induce programmed cell death (PCD), resulting in increased colonization of host tissue and, ultimately, sporulation to complete its pathogenic life cycle. Intensive research over the last two decades has led to the functional characterization of five proteinaceous necrotrophic effectors, SnTox1, SnToxA, SnTox267, SnTox3, and SnTox5, and three wheat susceptibility genes, Tsn1, Snn1, and Snn3D-1. Functional characterization has revealed that these effectors, in addition to inducing PCD, have additional roles in pathogenesis, including chitin binding that results in protection from wheat chitinases, blocking defense response signaling, and facilitating plant colonization. There are still large gaps in our understanding of how this necrotrophic pathogen is successfully manipulating wheat defense to complete its life cycle. This review summarizes our current knowledge, identifies knowledge gaps, and provides a summary of well-developed tools and resources currently available to study the P. nodorum-wheat interaction, which has become a model for necrotrophic specialist interactions. Further functional characterization of the effectors involved in this interaction and work toward a complete understanding of how P. nodorum manipulates wheat defense will provide fundamental knowledge about this and other necrotrophic interactions. Additionally, a broader understanding of this interaction will contribute to the successful management of Septoria nodorum blotch disease on wheat. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Virulence Associations and Global Context of AvrStb6 Genetic Diversity in Iranian Populations of Zymoseptoria tritici

Managing pathogen damage in wheat production is important for sustaining yields. Fungal plant pathogen genomes encode many small secreted proteins acting as effectors that play key roles in the successful colonization of host tissue and triggering host defenses. AvrStb6 is the first described Zymoseptoria tritici avirulence effector, which triggers Stb6-mediated immunity in the wheat host in a gene-for-gene manner. Evasion of major resistance factors such as Stb6 challenges deployment decisions on wheat cultivars. In this study, we analyzed the evolution of the AvrStb6 effector in Iranian isolates of Z. tritici. In total, 78 isolates were isolated and purified from 30 infected wheat specimens collected from the East Azerbaijan and Ardabil provinces of Iran. The pathogenicity of all isolates was evaluated on the susceptible wheat cultivar 'Tajan'. A subset of 40 isolates were also tested for pathogenicity on the resistant cultivar 'Shafir' carrying Stb6. Genetic diversity at the AvrStb6 locus was analyzed for 14 isolates covering the breadth of the observed disease severity. The AvrStb6 sequence variation was high, with virulent isolates carrying highly diverse AvrStb6 haplotypes. In an analysis including more than 1,000 additional AvrStb6 sequences from a global set of isolates, we found that virulent isolates carried AvrStb6 haplotypes either clustering with known virulent haplotypes on different continents or constituting previously unknown haplotypes. Furthermore, we found that AvrStb6 variants from avirulent isolates clustered with known avirulent genotypes from Europe. Our study highlights the relevance of AvrStb6 for Z. tritici virulence and the exceptional global diversity patterns of this effector.

Genome-wide association mapping of resistance to the foliar diseases septoria nodorum blotch and tan spot in a global winter wheat collection

Septoria nodorum blotch (SNB) and tan spot, caused by the necrotrophic fungal pathogens Parastagonospora nodorum and Pyrenophora tritici-repentis, respectively, often occur together as a leaf spotting disease complex on wheat (Triticum aestivum L.). Both pathogens produce necrotrophic effectors (NEs) that contribute to the development of disease. Here, genome-wide association analysis of a diverse panel of 264 winter wheat lines revealed novel loci on chromosomes 5A and 5B associated with sensitivity to the NEs SnTox3 and SnTox5 in addition to the known sensitivity genes for NEs Ptr/SnToxA, SnTox1, SnTox3, and SnTox5. Sensitivity loci for SnTox267 and Ptr ToxB were not detected. Evaluation of the panel with five P. nodorum isolates for SNB development indicated the Snn3-SnTox3 and Tsn1-SnToxA interactions played significant roles in disease development along with additional QTL on chromosomes 2A and 2D, which may correspond to the Snn7-SnTox267 interaction. For tan spot, the Tsc1-Ptr ToxC interaction was associated with disease caused by two isolates, and a novel QTL on chromosome 7D was associated with a third isolate. The Tsn1-ToxA interaction was associated with SNB but not tan spot. Therefore some, but not all, of the previously characterized host gene-NE interactions in these pathosystems play significant roles in disease development in winter wheat. Based on these results, breeders should prioritize the selection of resistance alleles at the Tsc1, Tsn1, Snn3, and Snn7 loci as well as the 2A and 7D QTL to obtain good levels of resistance to SNB and tan spot in winter wheat.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01400-5.

A Revised Nomenclature for ToxA Haplotypes Across Multiple Fungal Species

ToxA is one of the most studied proteinaceous necrotrophic effectors produced by plant pathogens. It has been identified in four pathogens (Pyrenophora tritici-repentis, Parastagonospora nodorum, Parastagonospora pseudonodorum [formerly Parastagonospora avenaria f. sp. tritici], and Bipolaris sorokiniana) causing leaf spot diseases on cereals worldwide. To date, 24 different ToxA haplotypes have been identified. Some P. tritici-repentis and related species also express ToxB, another small protein necrotrophic effector. We present here a revised and standardized nomenclature for these effectors, which could be extended to other poly-haplotypic genes found across multiple species.

Genetics and breeding for resistance against four leaf spot diseases in wheat (Triticum aestivum L.)

In wheat, major yield losses are caused by a variety of diseases including rusts, spike diseases, leaf spot and root diseases. The genetics of resistance against all these diseases have been studied in great detail and utilized for breeding resistant cultivars. The resistance against leaf spot diseases caused by each individual necrotroph/hemi-biotroph involves a complex system involving resistance (R) genes, sensitivity (S) genes, small secreted protein (SSP) genes and quantitative resistance loci (QRLs). This review deals with resistance for the following four-leaf spot diseases: (i) Septoria nodorum blotch (SNB) caused by Parastagonospora nodorum; (ii) Tan spot (TS) caused by Pyrenophora tritici-repentis; (iii) Spot blotch (SB) caused by Bipolaris sorokiniana and (iv) Septoria tritici blotch (STB) caused by Zymoseptoria tritici.

Genetics of resistance to septoria nodorum blotch in wheat

Septoria nodorum blotch (SNB) is a foliar disease of wheat caused by the necrotrophic fungal pathogen Parastagonospora nodorum. Research over the last two decades has shown that the wheat-P. nodorum pathosystem mostly follows an inverse gene-for-gene model. The fungus produces necrotrophic effectors (NEs) that interact with specific host gene products encoded by dominant sensitivity (S) genes. When a compatible interaction occurs, a 'defense response' in the host leads to programmed cell death thereby provided dead/dying cells from which the pathogen, being a necrotroph, can acquire nutrients allowing it to grow and sporulate. To date, nine S gene-NE interactions have been characterized in this pathosystem. Five NE-encoding genes, SnTox1, SnTox3, SnToxA, SnTox5, and SnTox267, have been cloned along with three host S genes, Tsn1, Snn1, and Snn3-D1. Studies have shown that P. nodorum hijacks multiple and diverse host targets to cause disease. SNB resistance is often quantitative in nature because multiple compatible interactions usually occur concomitantly. NE gene expression plays a key role in disease severity, and the effect of each compatible interaction can vary depending on the other existing compatible interactions. Numerous SNB-resistance QTL have been identified in addition to the known S genes, and more research is needed to understand the nature of these resistance loci. Marker-assisted elimination of S genes through conventional breeding practices and disruption of S genes using gene editing techniques are both effective strategies for the development of SNB-resistant wheat cultivars, which will become necessary as the global demand for sustenance grows.

Identification of a Novel ToxA Haplotype of Pyrenophora tritici-repentis from Japan

Pyrenophora tritici-repentis was described first as a pathogen of wheat (tan spot) in Japan in the 1920s, but, since then, no reports on P. tritici-repentis race structure or its effectors in Japan have been published. In this study, 10 single-spore isolates of P. tritici-repentis were collected from bread wheat in Japan. These isolates were evaluated for virulence on four differential wheat genotypes and tested for the presence/absence of the effector-encoding genes, ToxA and ToxB, in multiplex PCR assays. These isolates were identified as ToxA producers, of which eight were designated as race 2 (ToxA producers) and two were classified as race 1 (ToxA and ToxC producers) based on their virulence patterns. Sequence analysis of the ToxA amplicons from these 10 isolates indicated the presence of a novel ToxA haplotype (denoted PtrA2). A comparative sequence analysis and resequencing of ToxA from reference P. tritici-repentis isolates showed that all previously published ToxA haplotypes in P. tritici-repentis were identical, and are hence denoted PtrA1 in this study. A total of 163 PtrToxA sequences from global origins were already deposited in GenBank and were confirmed identical to PtrA1. Sequence variation in PtrA1 and PtrA2 open reading frames were found at three positions: one synonymous mutation at position 412 (C/G) and two nonsynonymous mutations at positions 342 and 362 that alter amino acid sequence. These mutations did not seem to affect the necrosis development on a ToxA-sensitive wheat genotype when rated for symptoms 5 to 7 days after inoculation. This is the first report correctly confirming the presence of an additional novel ToxA haplotype in P. tritici-repentis for which we have predicted its isoform and updated the ToxA haplotype evolutionary network.

Whole-genome sequence analysis of Bipolaris sorokiniana infecting wheat in India and characterization of ToxA gene in different isolates as pathogenicity determinants

Spot blotch disease of wheat caused by Bipolaris sorokiniana Boerma (Sacc.) is an emerging problem in South Asian countries. Whole genome of a highly virulent isolate of B. sorokiniana BS112 (BHU, Uttar Pradesh; Accession no. GCA_004329375.1) was sequenced using a hybrid assembly approach. Secreted proteins, virulence gene(s), pathogenicity-related gene(s) were identified and the role of ToxA gene present in this genome, in the development of disease was recognized. ToxA gene (535 bp) was analyzed and identified in the genome of B. sorokiniana (BS112) which revealed 100% homology with the ToxA gene of Pyrenophora tritici repentis (Accession no. MH017419). Furthermore, ToxA gene was amplified, sequenced and validated in 39 isolates of B. sorokiniana which confirmed the presence of ToxA gene in all the isolates taken for this study. All ToxA sequences were submitted in NCBI database (MN601358-MN601396). As ToxA gene interacts with Tsn1 gene of host, 13 wheat genotypes were evaluated out of which 5 genotypes (38.4%) were found to be Tsn1 positive with more severe necrotic lesions compared to Tsn1-negative wheat genotypes. In vitro expression analysis of ToxA gene in the pathogen B. sorokiniana using qPCR revealed maximum upregulation (14.67 fold) at 1st day after inoculation (DAI) in the medium. Furthermore, in planta expression analysis of ToxA gene in Tsn1-positive and Tsn1-negative genotypes, revealed maximum expression (7.89-fold) in Tsn1-positive genotype, Agra local at 5th DAI compared to Tsn1-negative genotype Chiriya 7 showing minimum expression (0.048-fold) at 5th DAI. In planta ToxA-Tsn1 interaction studies suggested that spot blotch disease is more severe in Tsn1-positive genotypes, which will be helpful in better understanding and management of spot blotch disease of wheat.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03213-3.

A Conserved Hypothetical Gene Is Required but Not Sufficient for Ptr ToxC Production in Pyrenophora tritici-repentis

The fungus Pyrenophora tritici-repentis causes tan spot, an important foliar disease of wheat worldwide. The fungal pathogen produces three necrotrophic effectors, namely Ptr ToxA, Ptr ToxB, and Ptr ToxC to induce necrosis or chlorosis in wheat. Both Ptr ToxA and Ptr ToxB are proteins, and their encoding genes have been cloned. Ptr ToxC was characterized as a low-molecular weight molecule 20 years ago but the one or more genes controlling its production in P. tritici-repentis are unknown. Here, we report the genetic mapping, molecular cloning, and functional analysis of a fungal gene that is required for Ptr ToxC production. The genetic locus controlling the production of Ptr ToxC, termed ToxC, was mapped to a subtelomeric region using segregating biparental populations, genome sequencing, and association analysis. Additional marker analysis further delimited ToxC to a 173-kb region. The predicted genes in the region were examined for presence/absence polymorphism in different races and isolates leading to the identification of a single candidate gene. Functional validation showed that this gene was required but not sufficient for Ptr ToxC production, thus it is designated as ToxC1. ToxC1 encoded a conserved hypothetical protein likely located on the vacuole membrane. The gene was highly expressed during infection, and only one haplotype was identified among 120 isolates sequenced. Our work suggests that Ptr ToxC is not a protein and is likely produced through a cascade of biosynthetic pathway. The identification of ToxC1 is a major step toward revealing the Ptr ToxC biosynthetic pathway and studying its molecular interactions with host factors.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

A triple threat: the Parastagonospora nodorum SnTox267 effector exploits three distinct host genetic factors to cause disease in wheat

Parastagonospora nodorum is a fungal pathogen of wheat. As a necrotrophic specialist, it deploys effector proteins that target dominant host susceptibility genes to elicit programmed cell death (PCD). Here we identify and functionally validate the effector targeting the host susceptibility genes Snn2, Snn6 and Snn7. We utilized whole-genome sequencing, association mapping, gene-disrupted mutants, gain-of-function transformants, virulence assays, bioinformatics and quantitative PCR to characterize these interactions. A single proteinaceous effector, SnTox267, targeted Snn2, Snn6 and Snn7 to trigger PCD. Snn2 and Snn6 functioned cooperatively to trigger PCD in a light-dependent pathway, whereas Snn7-mediated PCD functioned in a light-independent pathway. Isolates harboring 20 SnTox267 protein isoforms quantitatively varied in virulence. The diversity and distribution of isoforms varied between populations, indicating adaptation to local selection pressures. SnTox267 deletion resulted in the upregulation of effector genes SnToxA, SnTox1 and SnTox3. We validated a novel effector operating in an inverse-gene-for-gene manner to target three genetically distinct host susceptibility genes and elicit PCD. The discovery of the complementary gene action of Snn2 and Snn6 indicates their potential function in a guard or decoy model. Additionally, differences in light dependency in the elicited pathways and upregulation of unlinked effectors sheds new light onto a complex fungal necrotroph-host interaction.

Variability in an effector gene promoter of a necrotrophic fungal pathogen dictates epistasis and effector-triggered susceptibility in wheat

The fungus Parastagonospora nodorum uses proteinaceous necrotrophic effectors (NEs) to induce tissue necrosis on wheat leaves during infection, leading to the symptoms of septoria nodorum blotch (SNB). The NEs Tox1 and Tox3 induce necrosis on wheat possessing the dominant susceptibility genes Snn1 and Snn3B1/Snn3D1, respectively. We previously observed that Tox1 is epistatic to the expression of Tox3 and a quantitative trait locus (QTL) on chromosome 2A that contributes to SNB resistance/susceptibility. The expression of Tox1 is significantly higher in the Australian strain SN15 compared to the American strain SN4. Inspection of the Tox1 promoter region revealed a 401 bp promoter genetic element in SN4 positioned 267 bp upstream of the start codon that is absent in SN15, called PE401. Analysis of the world-wide P. nodorum population revealed that a high proportion of Northern Hemisphere isolates possess PE401 whereas the opposite was observed in representative P. nodorum isolates from Australia and South Africa. The presence of PE401 removed the epistatic effect of Tox1 on the contribution of the SNB 2A QTL but not Tox3. PE401 was introduced into the Tox1 promoter regulatory region in SN15 to test for direct regulatory roles. Tox1 expression was markedly reduced in the presence of PE401. This suggests a repressor molecule(s) binds PE401 and inhibits Tox1 transcription. Infection assays also demonstrated that P. nodorum which lacks PE401 is more pathogenic on Snn1 wheat varieties than P. nodorum carrying PE401. An infection competition assay between P. nodorum isogenic strains with and without PE401 indicated that the higher Tox1-expressing strain rescued the reduced virulence of the lower Tox1-expressing strain on Snn1 wheat. Our study demonstrated that Tox1 exhibits both 'selfish' and 'altruistic' characteristics. This offers an insight into a complex NE-NE interaction that is occurring within the P. nodorum population. The importance of PE401 in breeding for SNB resistance in wheat is discussed.

The crystal structure of SnTox3 from the necrotrophic fungus Parastagonospora nodorum reveals a unique effector fold and provides insight into Snn3 recognition and pro-domain protease processing of fungal effectors

Plant pathogens cause disease through secreted effector proteins, which act to promote infection. Typically, the sequences of effectors provide little functional information and further targeted experimentation is required. Here, we utilized a structure/function approach to study SnTox3, an effector from the necrotrophic fungal pathogen Parastagonospora nodorum, which causes cell death in wheat-lines carrying the sensitivity gene Snn3. We developed a workflow for the production of SnTox3 in a heterologous host that enabled crystal structure determination and functional studies. We show this approach can be successfully applied to study effectors from other pathogenic fungi. The β-barrel fold of SnTox3 is a novel fold among fungal effectors. Structure-guided mutagenesis enabled the identification of residues required for Snn3 recognition. SnTox3 is a pre-pro-protein, and the pro-domain of SnTox3 can be cleaved in vitro by the protease Kex2. Complementing this, an in silico study uncovered the prevalence of a conserved motif (LxxR) in an expanded set of putative pro-domain-containing fungal effectors, some of which can be cleaved by Kex2 in vitro. Our in vitro and in silico study suggests that Kex2-processed pro-domain (designated here as K2PP) effectors are common in fungi and this may have broad implications for the approaches used to study their functions.

Remarkable recent changes in the genetic diversity of the avirulence gene AvrStb6 in global populations of the wheat pathogen Zymoseptoria tritici

Septoria tritici blotch (STB), caused by the fungus Zymoseptoria tritici, is one of the most economically important diseases of wheat. Recently, both factors of a gene-for-gene interaction between Z. tritici and wheat, the wheat receptor-like kinase Stb6 and the Z. tritici secreted effector protein AvrStb6, have been identified. Previous analyses revealed a high diversity of AvrStb6 haplotypes present in earlier Z. tritici isolate collections, with up to c.18% of analysed isolates possessing the avirulence isoform of AvrStb6 identical to that originally identified in the reference isolate IPO323. With Stb6 present in many commercial wheat cultivars globally, we aimed to assess potential changes in AvrStb6 genetic diversity and the incidence of haplotypes allowing evasion of Stb6-mediated resistance in more recent Z. tritici populations. Here we show, using targeted resequencing of AvrStb6, that this gene is universally present in field isolates sampled from major wheat-growing regions of the world in 2013-2017. However, in contrast to the data from previous AvrStb6 population studies, we report a complete absence of the originally described avirulence isoform of AvrStb6 amongst modern Z. tritici isolates. Moreover, a remarkably small number of haplotypes, each encoding AvrStb6 protein isoforms conditioning virulence on Stb6-containing wheat, were found to predominate among modern Z. tritici isolates. A single virulence isoform of AvrStb6 was found to be particularly abundant throughout the global population. These findings indicate that, despite the ability of Z. tritici to sexually reproduce on resistant hosts, AvrStb6 avirulence haplotypes tend to be eliminated in subsequent populations.

Characterization of Effector-Target Interactions in Necrotrophic Pathosystems Reveals Trends and Variation in Host Manipulation

Great strides have been made in defining the details of the plant defense response involving biotrophic fungal and bacterial pathogens. The groundwork for the current model was laid by H.H. Flor and others who defined the gene-for-gene hypothesis, which is now known to involve effector-triggered immunity (ETI). PAMP-triggered immunity (PTI) is also a highly effective response to most pathogens because of the recognition of common pathogen molecules by pattern recognition receptors. In this article, we consider the three pathogens that make up the foliar disease complex of wheat, Zymoseptoria tritici, Pyrenophora tritici-repentis, and Parastagonospora nodorum, to review the means by which necrotrophic pathogens circumvent, or outright hijack, the ETI and PTI pathways to cause disease.

Population genomics of transposable element activation in the highly repressive genome of an agricultural pathogen

The activity of transposable elements (TEs) can be an important driver of genetic diversity with TE-mediated mutations having a wide range of fitness consequences. To avoid deleterious effects of TE activity, some fungi have evolved highly sophisticated genomic defences to reduce TE proliferation across the genome. Repeat-induced point mutation (RIP) is a fungal-specific TE defence mechanism efficiently targeting duplicated sequences. The rapid accumulation of RIPs is expected to deactivate TEs over the course of a few generations. The evolutionary dynamics of TEs at the population level in a species with highly repressive genome defences is poorly understood. Here, we analyse 366 whole-genome sequences of Parastagonospora nodorum, a fungal pathogen of wheat with efficient RIP. A global population genomics analysis revealed high levels of genetic diversity and signs of frequent sexual recombination. Contrary to expectations for a species with RIP, we identified recent TE activity in multiple populations. The TE composition and copy numbers showed little divergence among global populations regardless of the demographic history. Miniature inverted-repeat transposable elements (MITEs) and terminal repeat retrotransposons in miniature (TRIMs) were largely underlying recent intra-species TE expansions. We inferred RIP footprints in individual TE families and found that recently active, high-copy TEs have possibly evaded genomic defences. We find no evidence that recent positive selection acted on TE-mediated mutations rather that purifying selection maintained new TE insertions at low insertion frequencies in populations. Our findings highlight the complex evolutionary equilibria established by the joint action of TE activity, selection and genomic repression.

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.

Parastagonospora nodorum and Related Species in Western Canada: Genetic Variability and Effector Genes

Parastagonospora nodorum is an important fungal pathogen that causes Septoria nodorum blotch (SNB) in wheat. This pathogen produces several necrotrophic effectors that act as virulence factors; three have been cloned, SnToxA, SnTox1, and SnTox3. In this study, P. nodorum and its sister species P. avenaria f. tritici (Pat1) were isolated from wheat node and grain samples collected from distanced sites in western Canada during 2018. The presence of effector genes and associated haplotypes were determined by PCR and sequence analysis. An internal transcribed spacer-restriction fragment length polymorphism test was developed to distinguish between leaf spotting pathogens (P. nodorum, Pat1, Pyrenophora tritici-repentis, and Bipolaris sorokiniana). P. nodorum was mainly recovered from wheat nodes and to a lesser extent from the grains, while Pat1 was exclusively isolated from grain samples. The effector genes were present in almost all P. nodorum isolates, with the ToxA haplotype 5 (H5) being most prevalent, while a novel ToxA haplotype (denoted here H21) is reported for the first time. In Pat1, only combinations of SnTox1 and SnTox3 genes were present. A ToxA haplotype network was also constructed to assess the evolutionary relationship among globally found haplotypes to date. Finally, cultivars representing wheat development in Canada for the last century were tested for sensitivity to Sn-effectors and to the presence of Tsn1, the ToxA sensitivity gene. Of tested cultivars, 32.9 and 56.9% were sensitive to SnTox1 and SnTox3, respectively, and Tsn1 was present in 59% of the cultivars. In conclusion, P. nodorum and Pat1 were prevalent wheat pathogens in Canada with a potential tissue-specific colonization capacity, while producing necrotrophic effectors to which wheat is sensitive.

The Genetic Architecture of Emerging Fungicide Resistance in Populations of a Global Wheat Pathogen

Containing fungal diseases often depends on the application of fungicidal compounds. Fungicides can rapidly lose effectiveness due to the rise of resistant individuals in populations. However, the lack of knowledge about resistance mutations beyond known target genes challenges investigations into pathways to resistance. We used whole-genome sequencing data and association mapping to reveal the multilocus genetic architecture of fungicide resistance in a global panel of 159 isolates of Parastagonospora nodorum, an important fungal pathogen of wheat. We found significant differences in azole resistance among global field populations. The populations evolved distinctive combinations of resistance alleles which can interact when co-occurring in the same genetic background. We identified 34 significantly associated single nucleotide polymorphisms located in close proximity to genes associated with fungicide resistance in other fungi, including a major facilitator superfamily transporter. Using fungal colony growth rates and melanin production at different temperatures as fitness proxies, we found no evidence that resistance was constrained by genetic trade-offs. Our study demonstrates how genome-wide association studies of a global collection of pathogen strains can recapitulate the emergence of fungicide resistance. The distinct complement of resistance mutations found among populations illustrates how the evolutionary trajectory of fungicide adaptation can be complex and challenging to predict.

Genetic Structure of the Norwegian Parastagonospora nodorum Population

The necrotrophic fungal pathogen Parastagonospora nodorum causes Septoria nodorum blotch (SNB), which is one of the dominating leaf blotch diseases of wheat in Norway. A total of 165 P. nodorum isolates were collected from three wheat growing regions in Norway from 2015 to 2017. These isolates, as well as nine isolates from other countries, were analyzed for genetic variation using 20 simple sequence repeat (SSR) markers. Genetic analysis of the isolate collection indicated that the P. nodorum pathogen population infecting Norwegian spring and winter wheat underwent regular sexual reproduction and exhibited a high level of genetic diversity, with no genetic subdivisions between sampled locations, years or host cultivars. A high frequency of the presence of necrotrophic effector (NE) gene SnToxA was found in Norwegian P. nodorum isolates compared to other parts of Europe, and we hypothesize that the SnToxA gene is the major virulence factor among the three known P. nodorum NE genes (SnToxA, SnTox1, and SnTox3) in the Norwegian pathogen population. While the importance of SNB has declined in much of Europe, Norway has remained as a P. nodorum hotspot, likely due at least in part to local adaptation of the pathogen population to ToxA sensitive Norwegian spring wheat cultivars.

Natural selection drives population divergence for local adaptation in a wheat pathogen

Evolution favors the emergence of locally-adapted optimum phenotypes that are likely to differ across a wide array of environmental conditions. The emergence of favorable adaptive characteristics is accelerated in agricultural pathogens due to the unique properties of agro-ecosystems. We performed a Q_ST - F_ST comparison using 164 strains of Parastagonospora nodorum sampled from eight global field populations to disentangle the predominant evolutionary forces driving population divergence in a wheat pathogen. We used digital image analysis to obtain quantitative measurements of growth rate and melanization at different temperatures and under different fungicide concentrations in a common garden experiment. F_ST measures were based on complete genome sequences obtained for all 164 isolates. Our analyses indicated that all measured traits were under selection. Growth rates at 18 °C and 24 °C were under stabilizing selection (Q_ST < F_ST), while diversifying selection (Q_ST > F_ST) was the predominant evolutionary force affecting growth under fungicide and high temperature stress. Stabilizing selection (Q_ST < F_ST) was the predominant force affecting melanization across the different environments. Melanin production increased at 30 °C but was negatively correlated with higher growth rates, consistent with a trade-off under heat stress. Our results demonstrate that global populations of P. nodorum possess significant evolutionary potential to adapt to changing local conditions, including warmer temperatures and applications of fungicides.

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.

Role of Effector-Sensitivity Gene Interactions and Durability of Quantitative Resistance to Septoria Nodorum Blotch in Eastern U.S. Wheat

Important advances have been made in understanding the relationship of necrotrophic effectors (NE) and host sensitivity (Snn) genes in the Parastagonospora nodorum-wheat pathosystem. Yet much remains to be learned about the role of these interactions in determining wheat resistance levels in the field, and there is mixed evidence on whether breeding programs have selected against Snn genes due to their role in conferring susceptibility. SNB occurs ubiquitously in the U.S. Atlantic seaboard, and the environment is especially well suited to field studies of resistance to natural P. nodorum populations, as there are no other important wheat leaf blights. Insights into the nature of SNB resistance have been gleaned from multi-year data on phenotypes and markers in cultivars representative of the region's germplasm. In this perspective article, we review the evidence that in this eastern region of the U.S., wheat cultivars have durable quantitative SNB resistance and Snn-NE interactions are of limited importance. This conclusion is discussed in light of the relevant available information from other parts of the world.

Pyrenophora tritici-repentis in Tunisia: Race Structure and Effector Genes

Tan spot is a destructive foliar wheat disease worldwide and caused by the ascomycete fungus Pyrenophora tritici-repentis (Ptr); it has become more frequent in Tunisia over the last decade. In this study, the virulence of 73 single-spore isolates, collected from durum and bread wheat fields during 2017-2018 growing season, was evaluated on four differential wheat genotypes. This was followed by polymerase chain reaction tests with specific primers for the effector genes ToxA, ToxB, and toxb (ToxB-homolog). Sequence analysis to validate the identity of the amplified genes was followed, and ToxA amplicons from a subset of 22 isolates were analyzed to determine its haplotype identity. Ptr isolates from Tunisia were grouped in races 2, 4, 5, and 7, and 44% of the tested isolates did not fit under any known race, and were denoted here as atypical. These atypical isolates induced the same symptoms as race 7 isolates, extensive necrosis, and chlorosis on susceptible genotypes, but lacked the ToxA gene. ToxA is the only identified necrosis-inducing effector in Ptr, and was amplified in 51% of tested isolates, and shared identical sequence to previously identified haplotype (H15). ToxB and its homolog toxb were present in 97% and 93% of tested isolates, respectively. Ptr in Tunisia lacked Ptr ToxC activity, and none of the tested isolates induced the specific symptoms of that effector. Race 7 and the atypical isolates dominated the Tunisian Ptr population, while races 2, 4, and 5 were found at low percentages. In conclusion, ToxB and its homolog were the most dominant genes in Ptr from Tunisia, and the majority of the isolates induced necrosis and chlorosis on Ptr ToxA and Ptr ToxB susceptible wheat genotypes. However, only about half of that necrosis can be attributed to ToxA presence, this result necessitates further research to investigate the prevalence of additional necrotic effector(s). Terminology: in this paper, Pyrenophora tritici-repentis abbreviated as Ptr, the effectors are referred to by Ptr ToxA, Ptr ToxB and Ptr ToxC, and the genes coding for them are written in italic as ToxA, ToxB, and ToxC, respectively.

Local adaptation drives the diversification of effectors in the fungal wheat pathogen Parastagonospora nodorum in the United States

Filamentous fungi rapidly evolve in response to environmental selection pressures in part due to their genomic plasticity. Parastagonospora nodorum, a fungal pathogen of wheat and causal agent of septoria nodorum blotch, responds to selection pressure exerted by its host, influencing the gain, loss, or functional diversification of virulence determinants, known as effector genes. Whole genome resequencing of 197 P. nodorum isolates collected from spring, durum, and winter wheat production regions of the United States enabled the examination of effector diversity and genomic regions under selection specific to geographically discrete populations. 1,026,859 SNPs/InDels were used to identify novel loci, as well as SnToxA and SnTox3 as factors in disease. Genes displaying presence/absence variation, predicted effector genes, and genes localized on an accessory chromosome had significantly higher pN/pS ratios, indicating a higher rate of sequence evolution. Population structure analyses indicated two P. nodorum populations corresponding to the Upper Midwest (Population 1) and Southern/Eastern United States (Population 2). Prevalence of SnToxA varied greatly between the two populations which correlated with presence of the host sensitivity gene Tsn1 in the most prevalent cultivars in the corresponding regions. Additionally, 12 and 5 candidate effector genes were observed to be under diversifying selection among isolates from Population 1 and 2, respectively, but under purifying selection or neutrally evolving in the opposite population. Selective sweep analysis revealed 10 and 19 regions that had recently undergone positive selection in Population 1 and 2, respectively, involving 92 genes in total. When comparing genes with and without presence/absence variation, those genes exhibiting this variation were significantly closer to transposable elements. Taken together, these results indicate that P. nodorum is rapidly adapting to distinct selection pressures unique to spring and winter wheat production regions by rapid adaptive evolution and various routes of genomic diversification, potentially facilitated through transposable element activity.

Transposon-Mediated Horizontal Transfer of the Host-Specific Virulence Protein ToxA between Three Fungal Wheat Pathogens

This work dissects the tripartite horizontal transfer of ToxA, a gene that has a direct negative impact on global wheat yields. Defining the extent of horizontally transferred DNA is important because it can provide clues to the mechanisms that facilitate HGT. Our analysis of ToxA and its surrounding 14 kb suggests that this gene was horizontally transferred in two independent events, with one event likely facilitated by a type II DNA transposon. These horizontal transfer events are now in various processes of decay in each species due to the repeated insertion of new transposons and subsequent rounds of targeted mutation by a fungal genome defense mechanism known as repeat induced point mutation. This work highlights the role that HGT plays in the evolution of host adaptation in eukaryotic pathogens. It also increases the growing body of evidence indicating that transposons facilitate adaptive HGT events between fungi present in similar environments and hosts.

Most known examples of horizontal gene transfer (HGT) between eukaryotes are ancient. These events are identified primarily using phylogenetic methods on coding regions alone. Only rarely are there examples of HGT where noncoding DNA is also reported. The gene encoding the wheat virulence protein ToxA and the surrounding 14 kb is one of these rare examples. ToxA has been horizontally transferred between three fungal wheat pathogens (Parastagonospora nodorum, Pyrenophora tritici-repentis, and Bipolaris sorokiniana) as part of a conserved ∼14 kb element which contains coding and noncoding regions. Here we used long-read sequencing to define the extent of HGT between these three fungal species. Construction of near-chromosomal-level assemblies enabled identification of terminal inverted repeats on either end of the 14 kb region, typical of a type II DNA transposon. This is the first description of ToxA with complete transposon features, which we call ToxhAT. In all three species, ToxhAT resides in a large (140-to-250 kb) transposon-rich genomic island which is absent in isolates that do not carry the gene (annotated here as toxa⁻). We demonstrate that the horizontal transfer of ToxhAT between P. tritici-repentis and P. nodorum occurred as part of a large (∼80 kb) HGT which is now undergoing extensive decay. In B. sorokiniana, in contrast, ToxhAT and its resident genomic island are mobile within the genome. Together, these data provide insight into the noncoding regions that facilitate HGT between eukaryotes and into the genomic processes which mask the extent of HGT between these species.

Necrotrophic Exploitation and Subversion of Plant Defense: A Lifestyle or Just a Phase, and Implications in Breeding Resistance

Breeding disease-resistant plants is a critical, environmentally friendly component of any strategy to sustainably feed and clothe the 9.8 billion people expected to live on Earth by 2050. Here, I review current literature detailing plant defense responses as they relate to diverse biological outcomes; disease resistance, susceptibility, and establishment of mutualistic plant-microbial relationships. Of particular interest is the degree to which these outcomes are a function of plant-associated microorganisms' lifestyles; biotrophic, hemibiotrophic, necrotrophic, or mutualistic. For the sake of brevity, necrotrophic pathogens and the necrotrophic phase of pathogenicity are emphasized in this review, with special attention given to the host-specific pathogens that exploit defense. Defense responses related to generalist necrotrophs and mutualists are discussed in the context of excellent reviews by others. In addition, host evolutionary trade-offs of disease resistance with other desirable traits are considered in the context of breeding for durable disease resistance.

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.

Genome-wide evidence for divergent selection between populations of a major agricultural pathogen

The genetic and environmental homogeneity in agricultural ecosystems is thought to impose strong and uniform selection pressures. However, the impact of this selection on plant pathogen genomes remains largely unknown. We aimed to identify the proportion of the genome and the specific gene functions under positive selection in populations of the fungal wheat pathogen Zymoseptoria tritici. First, we performed genome scans in four field populations that were sampled from different continents and on distinct wheat cultivars to test which genomic regions are under recent selection. Based on extended haplotype homozygosity and composite likelihood ratio tests, we identified 384 and 81 selective sweeps affecting 4% and 0.5% of the 35 Mb core genome, respectively. We found differences both in the number and the position of selective sweeps across the genome between populations. Using a XtX‐based outlier detection approach, we identified 51 extremely divergent genomic regions between the allopatric populations, suggesting that divergent selection led to locally adapted pathogen populations. We performed an outlier detection analysis between two sympatric populations infecting two different wheat cultivars to identify evidence for host‐driven selection. Selective sweep regions harboured genes that are likely to play a role in successfully establishing host infections. We also identified secondary metabolite gene clusters and an enrichment in genes encoding transporter and protein localization functions. The latter gene functions mediate responses to environmental stress, including interactions with the host. The distinct gene functions under selection indicate that both local host genotypes and abiotic factors contributed to local adaptation.

Just the surface: advances in the discovery and characterization of necrotrophic wheat effectors

For many years pathogens of wheat have remained poorly understood. Hindered by an inaccessible host and the obligate nature of many of the pathogens, our understanding of these interactions has been limited compared to other more amenable pathosystems. However, breakthroughs over recent years have shed new light on diseases of wheat, particularly those caused by the genetically tractable necrotrophic pathogens. We now understand that many of the necrotrophic fungal pathogens do interact with wheat in a strict gene-for-gene relationship, and that pathogen and host partners in these interactions have now been identified. This improved understanding of necrotrophic effector biology has fundamentally changed the way we consider these important wheat diseases.

The discovery of the virulence gene ToxA in the wheat and barley pathogen Bipolaris sorokiniana

Bipolaris sorokiniana is the causal agent of multiple diseases on wheat and barley and is the primary constraint to cereal production throughout South Asia. Despite its significance, the molecular basis of disease is poorly understood. To address this, the genomes of three Australian isolates of B. sorokiniana were sequenced and screened for known pathogenicity genes. Sequence analysis revealed that the isolate BRIP10943 harboured the ToxA gene, which has been associated previously with disease in the wheat pathogens Parastagonospora nodorum and Pyrenophora tritici-repentis. Analysis of the regions flanking ToxA within B. sorokiniana revealed that it was embedded within a 12-kb genomic element nearly identical to the corresponding regions in P. nodorum and P. tritici-repentis. A screen of 35 Australian B. sorokiniana isolates confirmed that ToxA was present in 12 isolates. Sequencing of the ToxA genes within these isolates revealed two haplotypes, which differed by a single non-synonymous nucleotide substitution. Pathogenicity assays showed that a B. sorokiniana isolate harbouring ToxA was more virulent on wheat lines that contained the sensitivity gene when compared with a non-ToxA isolate. This work demonstrates that proteins that confer host-specific virulence can be horizontally acquired across multiple species. This acquisition can dramatically increase the virulence of pathogenic strains on susceptible cultivars, which, in an agricultural setting, can have devastating economic and social impacts.

Evolutionary analyses of the avirulence effector AvrStb6 in global populations of Zymoseptoria tritici identify candidate amino acids involved in recognition

We analysed the population genetic diversity of AvrStb6, the first avirulence gene cloned from the wheat pathogen Zymoseptoria tritici, using 142 Z. tritici strains sampled from four wheat fields growing on three continents. Although AvrStb6 was located in a recombination hotspot, it was found in every strain, with 71 polymorphic sites that produced 41 distinct DNA haplotypes encoding 30 AvrStb6 protein isoforms. An AvrStb6 homologue was found in the closest known relative, Z. pseudotritici, but not in three other closely related Zymoseptoria species, indicating that this gene has emerged in Zymoseptoria quite recently. Two AvrStb6 homologues with nucleotide similarities greater than 70% were identified on chromosome 10 in all Z. tritici isolates, suggesting that AvrStb6 belongs to a multigene family of candidate effectors that has expanded recently through gene duplication. The AvrStb6 sequences exhibited strong evidence for non-neutral evolution, including a large number of non-synonymous mutations, with significant positive diversifying selection operating on nine of the 82 codons. It appears that balancing selection is operating across the entire gene in natural field populations. There was also evidence for co-evolving codons within the gene that may reflect compensatory mutations associated with the evasion of recognition by Stb6. Intragenic recombination also appears to have affected the diversity of AvrStb6.

Mapping of SnTox3-Snn3 as a major determinant of field susceptibility to Septoria nodorum leaf blotch in the SHA3/CBRD × Naxos population

The effect of the SnTox3-Snn3 interaction was documented for the first time under natural infection at the adult plant stage in the field. Co-segregating SNP markers were identified. Parastagonospora nodorum is a necrotrophic pathogen of wheat, causing Septoria nodorum blotch (SNB) affecting both the leaf and glume. P. nodorum is the major leaf blotch pathogen on spring wheat in Norway. Resistance to the disease is quantitative, but several host-specific interactions between necrotrophic effectors (NEs) and host sensitivity (Snn) genes have been identified, playing a major role at the seedling stage. However, the effect of these interactions in the field under natural infection has not been investigated. In the present study, we saturated the genetic map of the recombinant inbred (RI) population SHA3/CBRD × Naxos using the Illumina 90 K SNP chip. The population had previously been evaluated for segregation of SNB susceptibility in field trials. Here, we infiltrated the population with the purified NEs SnToxA, SnTox1 and SnTox3, and mapped the Snn3 locus on 5BS based on sensitivity segregation and SNP marker data. We also conducted inoculation and culture filtrate (CF) infiltration experiments on the population with four selected P. nodorum isolates from Norway and North America. Remapping of quantitative trait loci (QTL) for field resistance showed that the SnTox3-Snn3 interaction could explain 24% of the phenotypic variation in the field, and more than 51% of the variation in seedling inoculations. To our knowledge, this is the first time the effect of this interaction has been documented at the adult plant stage under natural infection in the field.

Rapid emergence of pathogens in agro-ecosystems: global threats to agricultural sustainability and food security

Agricultural ecosystems are composed of genetically depauperate populations of crop plants grown at a high density and over large spatial scales, with the regional composition of crop species changing little from year to year. These environments are highly conducive for the emergence and dissemination of pathogens. The uniform host populations facilitate the specialization of pathogens to particular crop cultivars and allow the build-up of large population sizes. Population genetic and genomic studies have shed light on the evolutionary mechanisms underlying speciation processes, adaptive evolution and long-distance dispersal of highly damaging pathogens in agro-ecosystems. These studies document the speed with which pathogens evolve to overcome crop resistance genes and pesticides. They also show that crop pathogens can be disseminated very quickly across and among continents through human activities. In this review, we discuss how the peculiar architecture of agro-ecosystems facilitates pathogen emergence, evolution and dispersal. We present four example pathosystems that illustrate both pathogen specialization and pathogen speciation, including different time frames for emergence and different mechanisms underlying the emergence process. Lastly, we argue for a re-design of agro-ecosystems that embraces the concept of dynamic diversity to improve their resilience to pathogens. This article is part of the themed issue 'Tackling emerging fungal threats to animal health, food security and ecosystem resilience'.

Validation of Genome-Wide Association Studies as a Tool to Identify Virulence Factors in Parastagonospora nodorum

Parastagonospora nodorum is a necrotrophic fungal pathogen causing Septoria nodorum blotch on wheat. We have identified nine necrotrophic effector-host dominant sensitivity gene interactions, and we have cloned three of the necrotrophic effector genes, including SnToxA, SnTox1, and SnTox3. Because sexual populations of P. nodorum are difficult to develop under lab conditions, genome-wide association study (GWAS) is the best population genomic approach to identify genomic regions associated with traits using natural populations. In this article, we used a global collection of 191 P. nodorum isolates from which we identified 2,983 single-nucleotide polymorphism (SNP) markers and gene markers for SnToxA and SnTox3 to evaluate the power of GWAS on two popular wheat breeding lines that were sensitive to SnToxA and SnTox3. Strong marker trait associations (MTA) with P. nodorum virulence that mapped to SnTox3 and SnToxA were first identified using the marker set described above. A novel locus in the P. nodorum genome associated with virulence was also identified as a result of this analysis. To evaluate whether a sufficient level of marker saturation was available, we designed a set of primers every 1 kb in the genomic regions containing SnToxA and SnTox3. Polymerase chain reaction amplification was performed across the 191 isolates and the presence/absence polymorphism was scored and used as the genotype. The marker proximity necessary to identify MTA flanking SnToxA and SnTox3 ranged from 4 to 5 and 1 to 7 kb, respectively. Similar analysis was performed on the novel locus. Using a 45% missing data threshold, two more SNP were identified spanning a 4.6-kb genomic region at the novel locus. These results showed that the rate of linkage disequilibrium (LD) decay in P. nodorum and, likely, other fungi is high compared with plants and animals. The fast LD decay in P. nodorum is an advantage only if sufficient marker density is attained. Based on our results with the SnToxA and SnTox3 regions, markers are needed every 9 or 8 kb, respectively, or in every gene, to guarantee that genes associated with a quantitative trait such as virulence are not missed.

Population Genomics of Fungal and Oomycete Pathogens

We are entering a new era in plant pathology in which whole-genome sequences of many individuals of a pathogen species are becoming readily available. Population genomics aims to discover genetic mechanisms underlying phenotypes associated with adaptive traits such as pathogenicity, virulence, fungicide resistance, and host specialization, as genome sequences or large numbers of single nucleotide polymorphisms become readily available from multiple individuals of the same species. This emerging field encompasses detailed genetic analyses of natural populations, comparative genomic analyses of closely related species, identification of genes under selection, and linkage analyses involving association studies in natural populations or segregating populations resulting from crosses. The era of pathogen population genomics will provide new opportunities and challenges, requiring new computational and analytical tools. This review focuses on conceptual and methodological issues as well as the approaches to answering questions in population genomics. The major steps start with defining relevant biological and evolutionary questions, followed by sampling, genotyping, and phenotyping, and ending in analytical methods and interpretations. We provide examples of recent applications of population genomics to fungal and oomycete plant pathogens.

Plant-Pathogen Effectors: Cellular Probes Interfering with Plant Defenses in Spatial and Temporal Manners

Plants possess large arsenals of immune receptors capable of recognizing all pathogen classes. To cause disease, pathogenic organisms must be able to overcome physical barriers, suppress or evade immune perception, and derive nutrients from host tissues. Consequently, to facilitate some of these processes, pathogens secrete effector proteins that promote colonization. This review covers recent advances in the field of effector biology, focusing on conserved cellular processes targeted by effectors from diverse pathogens. The ability of effectors to facilitate pathogen entry into the host interior, suppress plant immune perception, and alter host physiology for pathogen benefit is discussed. Pathogens also deploy effectors in a spatial and temporal manner, depending on infection stage. Recent advances have also enhanced our understanding of effectors acting in specific plant organs and tissues. Effectors are excellent cellular probes that facilitate insight into biological processes as well as key points of vulnerability in plant immune signaling networks.

SnTox1, a Parastagonospora nodorum necrotrophic effector, is a dual-function protein that facilitates infection while protecting from wheat-produced chitinases

SnTox1 induces programmed cell death and the up-regulation of pathogenesis-related genes including chitinases. Additionally, SnTox1 has structural homology to several plant chitin-binding proteins. Therefore, we evaluated SnTox1 for chitin binding and localization. We transformed an avirulent strain of Parastagonospora nodorum as well as three nonpathogens of wheat (Triticum aestivum), including a necrotrophic pathogen of barley, a hemibiotrophic pathogen of sugar beet and a saprotroph, to evaluate the role of SnTox1 in infection and in protection from wheat chitinases. SnTox1 bound chitin and an SnTox1-green fluorescent fusion protein localized to the mycelial cell wall. Purified SnTox1 induced necrosis in the absence of the pathogen when sprayed on the leaf surface and appeared to remain on the leaf surface while inducing both epidermal and mesophyll cell death. SnTox1 protected the different fungi from chitinase degradation. SnTox1 was sufficient to change the host range of a necrotrophic pathogen but not a hemibiotroph or saprotroph. Collectively, this work shows that SnTox1 probably interacts with a receptor on the outside of the cell to induce cell death to acquire nutrients, but SnTox1 accomplishes a second role in that it protects against one aspect of the defense response, namely the effects of wheat chitinases.

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.

Utilizing Gene Tree Variation to Identify Candidate Effector Genes in Zymoseptoria tritici

Zymoseptoria tritici is a host-specific, necrotrophic pathogen of wheat. Infection by Z. tritici is characterized by its extended latent period, which typically lasts 2 wks, and is followed by extensive host cell death, and rapid proliferation of fungal biomass. This work characterizes the level of genomic variation in 13 isolates, for which we have measured virulence on 11 wheat cultivars with differential resistance genes. Between the reference isolate, IPO323, and the 13 Australian isolates we identified over 800,000 single nucleotide polymorphisms, of which ∼10% had an effect on the coding regions of the genome. Furthermore, we identified over 1700 probable presence/absence polymorphisms in genes across the Australian isolates using de novo assembly. Finally, we developed a gene tree sorting method that quickly identifies groups of isolates within a single gene alignment whose sequence haplotypes correspond with virulence scores on a single wheat cultivar. Using this method, we have identified < 100 candidate effector genes whose gene sequence correlates with virulence toward a wheat cultivar carrying a major resistance gene.

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

The necrotrophic fungus Parastagonospora nodorum is an important pathogen of one of the world's most economically important cereal crops, wheat (Triticum aestivum L.). P. nodorum produces necrotrophic protein effectors that mediate host cell death, providing nutrients for continuation of the infection process. The recent discovery of pathogen effectors has revolutionized disease resistance breeding for necrotrophic diseases in crop species, allowing often complex genetic resistance mechanisms to be broken down into constituent parts. To date, three effectors have been identified in P. nodorum. Here we use the effector, SnTox1, to screen 642 progeny from an eight-parent multiparent advanced generation inter-cross (i.e., MAGIC) population, genotyped with a 90,000-feature single-nucleotide polymorphism array. The MAGIC founders showed a range of sensitivity to SnTox1, with transgressive segregation evident in the progeny. SnTox1 sensitivity showed high heritability, with quantitative trait locus analyses fine-mapping the Snn1 locus to the short arm of chromosome 1B. In addition, a previously undescribed SnTox1 sensitivity locus was identified on the long arm of chromosome 5A, termed here QSnn.niab-5A.1. The peak single-nucleotide polymorphism for the Snn1 locus was converted to the KASP genotyping platform, providing breeders and researchers a simple and cheap diagnostic marker for allelic state at Snn1.

A ToxA-like protein from Cochliobolus heterostrophus induces light-dependent leaf necrosis and acts as a virulence factor with host selectivity on maize

ToxA, the first discovered fungal proteinaceous host-selective toxin (HST), was originally identified in 1989 from the tan spot fungus Pyrenophora tritici-repentis (Ptr). About 25years later, a homolog was identified in the leaf/glume blotch fungus Stagonospora nodorum (Parastagonospora nodorum), also a pathogen of wheat. Here we report the identification and function of a ToxA-like protein from the maize pathogen Cochliobolus heterostrophus (Ch) that possesses necrosis-inducing activity specifically against maize. ChToxA is encoded by a 535-bp open reading frame featuring a ToxA-specific intron with unusual splicing sites (5'-ATAAGT…TAC-3') at conserved positions relative to PtrToxA. The protein shows 64% similarity to PtrToxA and is predicted to adopt a similar three-dimensional structure, although lacking the arginyl-glycyl-aspartic acid (RGD) motif reported to be required for internalization into sensitive wheat mesophyll cells. Reverse-transcriptase PCR revealed that the ChTOXA gene expression is up-regulated in planta, relative to axenic culture. Plant assays indicated that the recombinant ChToxA protein induces light-dependent leaf necrosis in a host-selective manner on maize inbred lines. Gene deletion experiments confirmed that ChtoxA mutants are reduced in virulence on specific ChToxA-sensitive maize lines, relative to virulence caused by wild-type strains. Database searches identified potential ChToxA homologues in other plant-pathogenic ascomycetes. Sequence and phylogenetic analyses revealed that the corresponding ToxA-like proteins include one member recently shown to be associated with formation of penetration hypha. These results provide the first evidence that C. heterostrophus is capable of producing proteinaceous HSTs as virulence factors in addition to well-known secondary metabolite-type toxins produced biosynthetically by polyketide synthase megaenzymes. Further studies on ChToxA may provide new insights into effector evolution in host-pathogen interactions.

Blocking the QB-binding site of photosystem II by tenuazonic acid, a non-host-specific toxin of Alternaria alternata, activates singlet oxygen-mediated and EXECUTER-dependent signalling in Arabidopsis

Necrotrophic fungal pathogens produce toxic compounds that induce cell death in infected plants. Often, the primary targets of these toxins and the way a plant responds to them are not known. In the present work, the effect of tenuazonic acid (TeA), a non-host-specific toxin of Alternaria alternata, on Arabidopsis thaliana has been analysed. TeA blocks the QB -binding site at the acceptor side of photosystem II (PSII). As a result, charge recombination at the reaction centre (RC) of PSII is expected to enhance the formation of the excited triplet state of the RC chlorophyll that promotes generation of singlet oxygen ((1)O₂). (1)O₂ activates a signalling pathway that depends on the two EXECUTER (EX) proteins EX1 and EX2 and triggers a programmed cell death response. In seedlings treated with TeA at half-inhibition concentration (1)O₂-mediated and EX-dependent signalling is activated as indicated by the rapid and transient up-regulation of (1)O₂-responsive genes in wild type, and its suppression in ex1/ex2 mutants. Lesion formation occurs when seedlings are exposed to higher concentrations of TeA for a longer period of time. Under these conditions, the programmed cell death response triggered by (1)O₂-mediated and EX-dependent signalling is superimposed by other events that also contribute to lesion formation.