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. Parastagonospora nodorum is an economically important pathogen of wheat, employing proteinaceous effectors to cause disease. Recognition of effectors by host susceptibility genes often leads to the elicitation of programmed cell death. However, little is known on the correlation between effector diversity and the spatial distribution of host resistance/susceptibility or the genomic mechanisms of diversification. This research presents the genome resequencing of 197 P. nodorum isolates collected from spring, winter, and durum wheat production regions of the United States, enabling the investigation of genome dynamics and evolution. Results illustrate local adaptation to host resistance or susceptibility, as evidenced by population-specific evolution of predicted effector genes and positively selected selective sweeps. Predicted effector genes, genes exhibiting presence/absence variation, and genes residing on an accessory chromosome, were found to be diversifying more rapidly. Additionally, transposable elements were predicted to play a role in the maintenance or elimination of genes. A GWAS approach identified the previously reported SnToxA and SnTox3 as well as novel virulence candidates as major elicitors of disease. These results highlight the flexibility of the P. nodorum genome in response to population-specific selection pressures and illustrate the utility of whole genome resequencing for the identification of putative virulence mechanisms.

The rpg4 gene confers recessive resistance to several races of wheat stem rust (Puccinia graminis f. sp. tritici) and Rpg5 provides dominant resistance against isolates of the rye stem rust (P. graminis f. sp. secalis) in barley. The rpg4 and Rpg5 genes are tightly linked on chromosome 5H, and positional cloning using high-resolution populations clearly separated the genes, unambiguously identifying Rpg5; however, the identity of rpg4 remained unclear. High-resolution genotyping of critical recombinants at the rpg4/Rpg5 locus, designated here as rpg4-mediated resistance locus (RMRL) delimited two distinct yet tightly linked loci required for resistance, designated as RMRL1 and RMRL2. Utilizing virus-induced gene silencing, each gene at RMRL1, i.e., HvRga1 (a nucleotide-binding site leucine-rich repeat [NBS-LRR] domain gene), Rpg5 (an NBS-LRR-protein kinase domain gene), and HvAdf3 (an actin depolymerizing factor-like gene), was individually silenced followed by inoculation with P. graminis f. sp. tritici race QCCJ. Silencing each gene changed the reaction type from incompatible to compatible, indicating that all three genes are required for rpg4-mediated resistance. This stem rust resistance mechanism in barley follows the emerging theme of unrelated pairs of genetically linked NBS-LRR genes required for specific pathogen recognition and resistance. It also appears that actin cytoskeleton dynamics may play an important role in determining resistance against several races of stem rust in barley.

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.

Summary: Parastagonospora nodorum is an economically important necrotrophic fungal pathogen of wheat. Parastagonospora nodorum secretes necrotrophic effectors that target wheat susceptibility genes to induce programmed cell death (PCD). In this study, we cloned and functionally validated SnTox5 and characterized its role in pathogenesis. We used whole genome sequencing, genome‐wide association study (GWAS) mapping, CRISPR‐Cas9‐based gene disruption, gain‐of‐function transformation, quantitative trait locus (QTL) analysis, haplotype and isoform analysis, protein modeling, quantitative PCR, and laser confocal microscopy to validate SnTox5 and functionally characterize SnTox5. SnTox5 is a mature 16.26 kDa protein with high structural similarity to SnTox3. Wild‐type and mutant P. nodorum strains and wheat genotypes of SnTox5 and Snn5, respectively, were used to show that SnTox5 not only targets Snn5 to induce PCD but also facilitates the colonization of the mesophyll layer even in the absence of Snn5. Here we show that SnTox5 facilitates the efficient colonization of the mesophyll tissue and elicits PCD specific to host lines carrying Snn5. The homology to SnTox3 and the ability of SnTox5 to facilitate the colonizing of the mesophyll also suggest a role in the suppression of host defense before PCD induction.

Summary: 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.