Parastagonospora nodorum, the causal agent of Septoria nodorum blotch in wheat, has emerged as a model necrotrophic fungal organism for the study of host–microbe interactions. To date, three necrotrophic effectors have been identified and characterized from this pathogen, including SnToxA, SnTox1, and SnTox3. Necrotrophic effector identification was greatly aided by the development of a draft genome of Australian isolate SN15 via Sanger sequencing, yet it remained largely fragmented. This research presents the development of nearly finished genomes of P. nodorum isolates Sn4, Sn2000, and Sn79-1087 using long-read sequencing technology. RNAseq analysis of isolate Sn4, consisting of eight time points covering various developmental and infection stages, mediated the annotation of 13,379 genes. Analysis of these genomes revealed large-scale polymorphism between the three isolates, including the complete absence of contig 23 from isolate Sn79-1087, and a region of genome expansion on contig 10 in isolates Sn4 and Sn2000. Additionally, these genomes exhibit the hallmark characteristics of a “two-speed” genome, being partitioned into two distinct GC-equilibrated and AT-rich compartments. Interestingly, isolate Sn79-1087 contains a lower proportion of AT-rich segments, indicating a potential lack of evolutionary hotspots. These newly sequenced genomes, consisting of telomere-to-telomere assemblies of nearly all 23 P. nodorum chromosomes, provide a robust foundation for the further examination of effector biology and genome evolution.

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.

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.

A diverse collection of barley lines was phenotyped with three North American Pyrenophora teres f. teres isolates and association analyses detected 78 significant marker-trait associations at 16 genomic loci. Pyrenophora teres f. teres is a necrotrophic fungal pathogen and the causal agent of the economically important foliar disease net form net blotch (NFNB) of barley. The deployment of effective and durable resistance against P. teres f. teres has been hindered by the complexity of quantitative resistance and susceptibility. Several bi-parental mapping populations have been used to identify QTL associated with NFNB disease on all seven barley chromosomes. Here, we report the first genome-wide association study (GWAS) to detect marker-trait associations for resistance or susceptibility to P. teres f. teres. Geographically diverse barley genotypes from a world barley core collection (957) were genotyped with the Illumina barley iSelect chip and phenotyped with three P. teres f. teres isolates collected in two geographical regions of the USA (15A, 6A and LDNH04Ptt19). The best of nine regression models tested were identified for each isolate and used for association analysis resulting in the identification of 78 significant marker-trait associations (MTA; -log10p value >3.0). The MTA identified corresponded to 16 unique genomic loci as determined by analysis of local linkage disequilibrium between markers that did not meet a correlation threshold of R ≥ 0.1, indicating that the markers represented distinct loci. Five loci identified represent novel QTL and were designated QRptts-3HL, QRptts-4HS, QRptts-5HL.1, QRptts-5HL.2, and QRptts-7HL.1. In addition, 55 of the barley lines examined exhibited a high level of resistance to all three isolates and the SNP markers identified will provide useful genetic resources for barley breeding programs.

Net form net blotch, caused by the necrotrophic fungal pathogen Pyrenophora teres f. teres, is a destructive foliar disease of barley with the potential to cause significant yield loss in major production regions throughout the world. The complexity of the host–parasite genetic interactions in this pathosystem hinders the deployment of effective resistance in barley cultivars, warranting a deeper understanding of the interactions. Here, we report on the high-resolution mapping of the dominant susceptibility locus near the centromere of chromosome 6H in the barley cultivars Rika and Kombar, which are putatively targeted by necrotrophic effectors from P. teres f. teres isolates 6A and 15A, respectively. Utilization of progeny isolates derived from a cross of P. teres f. teres isolates 6A × 15A harboring single major virulence loci (VK1, VK2, and VR2) allowed for the Mendelization of single inverse gene-for-gene interactions in a high-resolution population consisting of 2976 Rika × Kombar recombinant gametes. Brachypodium distachyon synteny was exploited to develop and saturate the susceptibility region with markers, delimiting it to ∼0.24 cM and a partial physical map was constructed. This genetic and physical characterization further resolved the dominant susceptibility locus, designated Spt1 (susceptibility to P. teres f. teres). The high-resolution mapping and cosegregation of the Spt1.R and Spt1.K gene/s indicates tightly linked genes in repulsion or alleles possibly targeted by different necrotrophic effectors. Newly developed barley genomic resources greatly enhance the efficiency of positional cloning efforts in barley, as demonstrated by the Spt1 fine mapping and physical contig identification reported here.