Emergence and spread of the barley net blotch pathogen coincided with crop domestication and cultivation history

Fungal pathogens cause devastating disease in crops. Understanding the evolutionary origin of pathogens is essential to the prediction of future disease emergence and the potential of pathogens to disperse. The fungus Pyrenophora teres f. teres causes net form net blotch (NFNB), an economically significant disease of barley. In this study, we have used 104 P. teres f. teres genomes from four continents to explore the population structure and demographic history of the fungal pathogen. We showed that P. teres f. teres is structured into populations that tend to be geographically restricted to different regions. Using Multiple Sequentially Markovian Coalescent and machine learning approaches we demonstrated that the demographic history of the pathogen correlates with the history of barley, highlighting the importance of human migration and trade in spreading the pathogen. Exploring signatures of natural selection, we identified several population-specific selective sweeps that colocalized with genomic regions enriched in putative virulence genes, and loci previously identified as determinants of virulence specificities by quantitative trait locus analyses. This reflects rapid adaptation to local hosts and environmental conditions of P. teres f. teres as it spread with barley. Our research highlights how human activities can contribute to the spread of pathogens that significantly impact the productivity of field crops.

A Moroccan Pyrenophora teres f. teres Population Defeats Rpt5, the Broadly Effective Resistance on Barley Chromosome 6H

Net form net blotch (NFNB), caused by Pyrenophora teres f. teres, is an important barley disease. The centromeric region of barley chromosome 6H has often been associated with resistance or susceptibility to NFNB, including the broadly effective dominant resistance gene Rpt5 derived from barley line CIho 5791. We characterized a population of Moroccan P. teres f. teres isolates that had overcome Rpt5 resistance and identified quantitative trait loci (QTL) that were effective against these isolates. Eight Moroccan P. teres f. teres isolates were phenotyped on barley lines CIho 5791 and Tifang. Six isolates were virulent on CIho 5791, and two were avirulent. A CIho 5791 × Tifang recombinant inbred line (RIL) population was phenotyped with all eight isolates and confirmed the defeat of the 6H resistance locus formerly mapped as Rpt5 in barley line CI9819. A major QTL on chromosome 3H with the resistance allele derived from Tifang, as well as minor QTL, was identified and provided resistance against these isolates. F2 segregation ratios supported dominant inheritance for both the 3H and 6H resistance. Furthermore, inoculation of progeny isolates derived from a cross of P. teres f. teres isolates 0-1 (virulent on Tifang/avirulent on CIho 5791) and MorSM 40-3 (avirulent on Tifang/virulent on CIho 5791) onto the RIL and F2 populations determined that recombination between isolates can generate novel genotypes that overcome both resistance genes. Markers linked to the QTL identified in this study can be used to incorporate both resistance loci into elite barley cultivars for durable resistance.

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.

Association Mapping of Resistance to Tan Spot in the Global Durum Panel

Tan spot, caused by the necrotrophic fungal pathogen Pyrenophora tritici-repentis (Ptr), is an important disease of durum and common wheat worldwide. Compared with common wheat, less is known about the genetics and molecular basis of tan spot resistance in durum wheat. We evaluated 510 durum lines from the Global Durum Wheat Panel (GDP) for sensitivity to the necrotrophic effectors (NEs) Ptr ToxA and Ptr ToxB and for reaction to Ptr isolates representing races 1 to 5. Overall, susceptible durum lines were most prevalent in South Asia, the Middle East, and North Africa. Genome-wide association analysis showed that the resistance locus Tsr7 was significantly associated with tan spot caused by races 2 and 3, but not races 1, 4, or 5. The NE sensitivity genes Tsc1 and Tsc2 were associated with susceptibility to Ptr ToxC- and Ptr ToxB-producing isolates, respectively, but Tsn1 was not associated with tan spot caused by Ptr ToxA-producing isolates, which further validates that the Tsn1-Ptr ToxA interaction does not play a significant role in tan spot development in durum. A unique locus on chromosome arm 2AS was associated with tan spot caused by race 4, a race once considered avirulent. A novel trait characterized by expanding chlorosis leading to increased disease severity caused by the Ptr ToxB-producing race 5 isolate DW5 was identified, and this trait was governed by a locus on chromosome 5B. We recommend that durum breeders select resistance alleles at the Tsr7, Tsc1, Tsc2, and the chromosome 2AS loci to obtain broad resistance to tan spot.

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.

Genetic and physical localization of a major susceptibility gene to Pyrenophora teres f. maculata in barley

KEY MESSAGE

Genetic characterization of a major spot form net blotch susceptibility locus to using linkage mapping to identify a candidate gene and user-friendly markers in barley. Spot form net blotch (SFNB), caused by the necrotrophic fungal pathogen Pyrenophora teres f. maculata (Ptm), is an economically important foliar diseases in barley. Although various resistance loci have been identified, breeding for SFNB-resistant varieties has been hampered due to the complex virulence profile of Ptm populations. One resistance locus in the host may be effective against one specific isolate, but it may confer susceptibility to other isolates. A major susceptibility QTL on chromosome 7H, named Sptm1, was consistently identified in many studies. In the present study, we conduct fine mapping to localize Sptm1 with high resolution. A segregating population was developed from selected F₂ progenies of the cross Tradition (S) × PI 67381 (R), in which the disease phenotype was determined by the Sptm1 locus alone. Disease phenotypes of critical recombinants were confirmed in the following two consecutive generations. Genetic mapping anchored the Sptm1 gene to an ⁓400 kb region on chromosome 7H. Gene prediction and annotation identified six protein-coding genes in the delimited Sptm1 region, and the gene encoding a putative cold-responsive protein kinase was selected as a strong candidate. Therefore, providing fine localization and candidate of Sptm1 for functional validation, our study will facilitate the understanding of susceptibility mechanism underlying the barley-Ptm interaction and offers a potential target for gene editing to develop valuable materials with broad-spectrum resistance to SFNB.

Hybrid inferiority and genetic incompatibilities drive divergence of fungal pathogens infecting the same host

Agro-ecosystems provide environments that are conducive for rapid evolution and dispersal of plant pathogens. Previous studies have demonstrated that hybridization of crop pathogens can give rise to new lineages with altered virulence profiles. Currently, little is known about either the genetics of fungal pathogen hybridization or the mechanisms that may prevent hybridization between related species. The fungus Pyrenophora teres is a global pathogen of barley. P. teres that exists as two distinct lineages P. teres f. teres and P. teres f. maculata (Ptt and Ptm, respectively), which both infect barley but produce very distinct lesions and rarely interbreed. Interestingly, Ptt and Ptm can, by experimental mating, produce viable progenies. Here, we addressed the underlying genetics of reproductive barriers of P. teres. We hypothesize that Ptt and Ptm diverged in the past, possibly by adapting to distinct hosts, and only more recently colonized the same host in agricultural fields. Using experimental mating and in planta phenotyping in barley cultivars susceptible to both P. teres forms, we demonstrate that hybrids produce mixed infection phenotypes but overall show inferior pathogenic fitness relative to the pure parents. Based on analyses of 104 hybrid genomes, we identify signatures of negative epistasis between parental alleles at distinct loci (Dobzhansky-Müller incompatibilities). Most DMI regions are not involved in virulence but certain genes are predicted or known to play a role in virulence. These results potentially suggest that divergent niche adaptation - albeit in the same host plant - contributes to speciation in P. teres.

Optimization of highly efficient exogenous-DNA-free Cas9-ribonucleoprotein mediated gene editing in disease susceptibility loci in wheat (Triticum aestivum L.)

The advancement of precision engineering for crop trait improvement is important in the face of rapid population growth, climate change, and disease. To this end, targeted double-stranded break technology using RNA-guided Cas9 has been adopted widely for genome editing in plants. Agrobacterium or particle bombardment-based delivery of plasmids encoding Cas9 and guide RNA (gRNA) is common, but requires optimization of expression and often results in random integration of plasmid DNA into the plant genome. Recent advances have described gene editing by the delivery of Cas9 and gRNA as pre-assembled ribonucleoproteins (RNPs) into various plant tissues, but with moderate efficiency in resulting regenerated plants. In this report we describe significant improvements to Cas9-RNP mediated gene editing in wheat. We demonstrate that Cas9-RNP assays in protoplasts are a fast and effective tool for rational selection of optimal gRNAs for gene editing in regenerable immature embryos (IEs), and that high temperature treatment enhances gene editing rates in both tissue types. We also show that Cas9-mediated editing persists for at least 14 days in gold particle bombarded wheat IEs. The regenerated edited wheat plants in this work are recovered at high rates in the absence of exogenous DNA and selection. With this method, we produce knockouts of a set of three homoeologous genes and two pathogenic effector susceptibility genes, engineering insensitivity to corresponding necrotrophic effectors produced by Parastagonospora nodorum. The establishment of highly efficient, exogenous DNA-free gene editing technology holds promise for accelerated trait diversity production in an expansive array of crops.

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.

Host and pathogen genetics reveal an inverse gene-for-gene association in the P. teres f. maculata-barley pathosystem

Pathogen and host genetics were used to uncover an inverse gene-for-gene interaction where virulence genes from the pathogen Pyrenophora teres f. maculata target barley susceptibility genes, resulting in disease. Although models have been proposed to broadly explain how plants and pathogens interact and coevolve, each interaction evolves independently, resulting in various scenarios of host manipulation and plant defense. Spot form net blotch is a foliar disease of barley caused by Pyrenophora teres f. maculata. We developed a barley population (Hockett × PI 67381) segregating for resistance to a diverse set of P. teres f. maculata isolates. Quantitative trait locus analysis identified major loci on barley chromosomes (Chr) 2H and 7H associated with resistance/susceptibility. Subsequently, we used avirulent and virulent P. teres f. maculata isolates to develop a pathogen population, identifying two major virulence loci located on Chr1 and Chr2. To further characterize this host-pathogen interaction, progeny from the pathogen population harboring virulence alleles at either the Chr1 or Chr2 locus was phenotyped on the Hockett × PI 67381 population. Progeny harboring only the Chr1 virulence allele lost the barley Chr7H association but maintained the 2H association. Conversely, isolates harboring only the Chr2 virulence allele lost the barley Chr2H association but maintained the 7H association. Hockett × PI 67381 F₂ individuals showed susceptible/resistant ratios not significantly different than 15:1 and results from F₂ inoculations using the single virulence genotypes were not significantly different from a 3:1 (S:R) ratio, indicating two dominant susceptibility genes. Collectively, this work shows that P. teres f. maculata virulence alleles at the Chr1 and Chr2 loci are targeting the barley 2H and 7H susceptibility alleles in an inverse gene-for-gene manner to facilitate colonization.

Haplotype variants of Sr46 in Aegilops tauschii, the diploid D genome progenitor of wheat

Stem rust resistance genes, SrRL5271 and Sr672.1 as well as SrCPI110651, from Aegilops tauschii, the diploid D genome progenitor of wheat, are sequence variants of Sr46 differing by 1-2 nucleotides leading to non-synonymous amino acid substitutions. The Aegilops tauschii (wheat D-genome progenitor) accessions RL 5271 and CPI110672 were identified as resistant to multiple races (including the Ug99) of the wheat stem rust pathogen Puccinia graminis f. sp. tritici (Pgt). This study was conducted to identify the stem rust resistance (Sr) gene(s) in both accessions. Genetic analysis of the resistance in RL 5271 identified a single dominant allele (SrRL5271) controlling resistance, whereas resistance segregated at two loci (SR672.1 and SR672.2) for a cross of CPI110672. Bulked segregant analysis placed SrRL5271 and Sr672.1 in a region on chromosome arm 2DS that encodes Sr46. Molecular marker screening, mapping and genomic sequence analysis demonstrated SrRL5271 and Sr672.1 are alleles of Sr46. The amino acid sequence of SrRL5271 and Sr672.1 is identical but differs from Sr46 (hereafter referred to as Sr46_h1 by following the gene nomenclature in wheat) by a single amino acid (N763K) and is thus designated Sr46_h2. Screening of a panel of Ae. tauschii accessions identified an additional allelic variant that differed from Sr46_h2 by a different amino acid (A648V) and was designated Sr46_h3. By contrast, the protein encoded by the susceptible allele of Ae. tauschii accession AL8/78 differed from these resistance proteins by 54 amino acid substitutions (94% nucleotide sequence gene identity). Cloning and complementation tests of the three resistance haplotypes confirmed their resistance to Pgt race 98-1,2,3,5,6 and partial resistance to Pgt race TTRTF in bread wheat. The three Sr46 haplotypes, with no virulent races detected yet, represent a valuable source for improving stem resistance in wheat.

Association mapping reveals a reciprocal virulence/avirulence locus within diverse US Pyrenophora teres f. maculata isolates

Background

Spot form net blotch (SFNB) caused by the necrotrophic fungal pathogen Pyrenophora teres f. maculata (Ptm) is an economically important disease of barley that also infects wheat. Using genetic analysis to characterize loci in Ptm genomes associated with virulence or avirulence is an important step to identify pathogen effectors that determine compatible (virulent) or incompatible (avirulent) interactions with cereal hosts. Association mapping (AM) is a powerful tool for detecting virulence loci utilizing phenotyping and genotyping data generated for natural populations of plant pathogenic fungi.

Results

Restriction-site associated DNA genotyping-by-sequencing (RAD-GBS) was used to generate 4,836 single nucleotide polymorphism (SNP) markers for a natural population of 103 Ptm isolates collected from Idaho, Montana and North Dakota. Association mapping analyses were performed utilizing the genotyping and infection type data generated for each isolate when challenged on barley seedlings of thirty SFNB differential barley lines. A total of 39 marker trait associations (MTAs) were detected across the 20 barley lines corresponding to 30 quantitative trait loci (QTL); 26 novel QTL and four that were previously mapped in Ptm biparental populations. These results using diverse US isolates and barley lines showed numerous barley-Ptm genetic interactions with seven of the 30 Ptm virulence/avirulence loci falling on chromosome 3, suggesting that it is a reservoir of diverse virulence effectors. One of the loci exhibited reciprocal virulence/avirulence with one haplotype predominantly present in isolates collected from Idaho increasing virulence on barley line MXB468 and the alternative haplotype predominantly present in isolates collected from North Dakota and Montana increasing virulence on barley line CI9819.

Conclusions

Association mapping provided novel insight into the host pathogen genetic interactions occurring in the barley-Ptm pathosystem. The analysis suggests that chromosome 3 of Ptm serves as an effector reservoir in concordance with previous reports for Pyrenophora teres f. teres, the causal agent of the closely related disease net form net blotch. Additionally, these analyses identified the first reported case of a reciprocal pathogen virulence locus. However, further investigation of the pathosystem is required to determine if multiple genes or alleles of the same gene are responsible for this genetic phenomenon.

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.

The Parastagonospora nodorum necrotrophic effector SnTox5 targets the wheat gene Snn5 and facilitates entry into the leaf mesophyll

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.

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.

Mutations in a barley cytochrome P450 gene enhances pathogen induced programmed cell death and cutin layer instability

Disease lesion mimic mutants (DLMMs) are characterized by the spontaneous development of necrotic spots with various phenotypes designated as necrotic (nec) mutants in barley. The nec mutants were traditionally considered to have aberrant regulation of programmed cell death (PCD) pathways, which have roles in plant immunity and development. Most barley nec3 mutants express cream to orange necrotic lesions contrasting them from typical spontaneous DLMMs that develop dark pigmented lesions indicative of serotonin/phenolics deposition. Barley nec3 mutants grown under sterile conditions did not exhibit necrotic phenotypes until inoculated with adapted pathogens, suggesting that they are not typical DLMMs. The F₂ progeny of a cross between nec3-γ1 and variety Quest segregated as a single recessive susceptibility gene post-inoculation with Bipolaris sorokiniana, the causal agent of the disease spot blotch. Nec3 was genetically delimited to 0.14 cM representing 16.5 megabases of physical sequence containing 149 annotated high confidence genes. RNAseq and comparative analysis of the wild type and five independent nec3 mutants identified a single candidate cytochrome P450 gene (HORVU.MOREX.r2.6HG0460850) that was validated as nec3 by independent mutations that result in predicted nonfunctional proteins. Histology studies determined that nec3 mutants had an unstable cutin layer that disrupted normal Bipolaris sorokiniana germ tube development.

At the site of pathogen infection, plant defense mechanisms rely on controlled programmed cell death (PCD) to sequester biotrophic pathogens that require living cells to extract nutrients from the host. However, these defense mechanisms are hijacked by necrotrophic plant pathogens that purposefully induce PCD to feed on the dead cells, thus facilitating further disease development. Thus, understanding PCD responses is important for resistance to both classes of pathogens. We characterized five independent disease lesion mimic mutants of barley designated necrotic 3 (nec3) that show aberrant regulation of PCD responses upon pathogen challenge. A cytochrome P450 gene was identified as Nec3 encoding a Tryptamine 5-Hydroxylase that functions as a terminal serotonin biosynthetic enzyme in the Tryptophan pathway of plants. We posit that nec3 mutants have disrupted serotonin biosynthesis resulting in expanded PCD, necrotrophic pathogen susceptibility and cutin layer instability. The nec3 mutants show expanded PCD and disease susceptibility of pathogen-induced necrotic lesions, suggesting a role of serotonin to sequester PCD and suppress pathogen colonization. The identification of Nec3 will facilitate functional analysis to elucidate the role that serotonin plays in the elicitation or suppression of PCD immunity responses to diverse pathogens and the effects it has on cutin layer biosynthesis.

Genetic mapping of host resistance to the Pyrenophora teres f. maculata isolate 13IM8.3

Spot form net blotch (SFNB), caused by the necrotrophic fungal pathogen Pyrenophora teres f. maculata (Ptm), is a foliar disease of barley that results in significant yield losses in major growing regions worldwide. Understanding the host-parasite interactions between pathogen virulence/avirulence genes and the corresponding host susceptibility/resistance genes is important for the deployment of genetic resistance against SFNB. Two recombinant inbred mapping populations were developed to characterize genetic resistance/susceptibility to the Ptm isolate 13IM8.3, which was collected from Idaho (ID). An Illumina Infinium array was used to produce a genome wide marker set. Quantitative trait loci (QTL) analysis identified ten significant resistance/susceptibility loci, with two of the QTL being common to both populations. One of the QTL on 5H appears to be novel, while the remaining loci have been reported previously. Single nucleotide polymorphisms (SNPs) closely linked to or delimiting the significant QTL have been converted to user-friendly markers. Loci and associated molecular markers identified in this study will be useful in genetic mapping and deployment of the genetic resistance to SFNB in barley.

Genome-Wide Association and Selective Sweep Studies Reveal the Complex Genetic Architecture of DMI Fungicide Resistance in Cercospora beticola

The rapid and widespread evolution of fungicide resistance remains a challenge for crop disease management. The demethylation inhibitor (DMI) class of fungicides is a widely used chemistry for managing disease, but there has been a gradual decline in efficacy in many crop pathosystems. Reliance on DMI fungicides has increased resistance in populations of the plant pathogenic fungus Cercospora beticola worldwide. To better understand the genetic and evolutionary basis for DMI resistance in C. beticola, a genome-wide association study (GWAS) and selective sweep analysis were conducted for the first time in this species. We performed whole-genome resequencing of 190 C. beticola isolates infecting sugar beet (Beta vulgaris ssp. vulgaris). All isolates were phenotyped for sensitivity to the DMI tetraconazole. Intragenic markers on chromosomes 1, 4, and 9 were significantly associated with DMI fungicide resistance, including a polyketide synthase gene and the gene encoding the DMI target CbCYP51. Haplotype analysis of CbCYP51 identified a synonymous mutation (E170) and nonsynonymous mutations (L144F, I387M, and Y464S) associated with DMI resistance. Genome-wide scans of selection showed that several of the GWAS mutations for fungicide resistance resided in regions that have recently undergone a selective sweep. Using radial plate growth on selected media as a fitness proxy, we did not find a trade-off associated with DMI fungicide resistance. Taken together, we show that population genomic data from a crop pathogen can allow the identification of mutations conferring fungicide resistance and inform about their origins in the pathogen population.

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.

Function and evolution of allelic variations of Sr13 conferring resistance to stem rust in tetraploid wheat (Triticum turgidum L.)

The resistance gene Sr13 is one of the most important genes in durum wheat for controlling stem rust caused by Puccinia graminis f. sp. tritici (Pgt). The Sr13 functional gene CNL13 has haplotypes R1, R2 and R3. The R1/R3 and R2 haplotypes were originally designated as alleles Sr13a and Sr13b, respectively. To detect additional Sr13 alleles, we developed Kompetitive allele specific PCR (KASP™) marker KASPSr13 and four semi-thermal asymmetric reverse PCR markers, rwgsnp37-rwgsnp40, based on the CNL13 sequence. These markers were shown to detect R1, R2 and R3 haplotypes in a panel of diverse tetraploid wheat accessions. We also observed the presence of Sr13 in durum line CAT-A1, although it lacked any of the known haplotypes. Sequence analysis revealed that CNL13 of CAT-A1 differed from the susceptible haplotype S1 by a single nucleotide (C2200T) in the leucine-rich repeat region and differed from the other three R haplotypes by one or two additional nucleotides, confirming that CAT-A1 carries a new (R4) haplotype. Stem rust tests on the monogenic, transgenic and mutant lines showed that R1 differed from R3 in its susceptibility to races TCMJC and THTSC, whereas R4 differed from all other haplotypes for susceptibility to TTKSK, TPPKC and TCCJC. Based on these differences, we designate the R1, R3 and R4 haplotypes as alleles Sr13a, Sr13c and Sr13d, respectively. This study indicates that Sr13d may be the primitive functional allele originating from the S1 haplotype via a point mutation, with the other three R alleles probably being derived from Sr13d through one or two additional point mutations.

A genome-wide genetic linkage map and reference quality genome sequence for a new race in the wheat pathogen Pyrenophora tritici-repentis

Pyrenophora tritici-repentis is an ascomycete fungus that causes tan spot of wheat. The disease has a worldwide distribution and can cause significant yield and quality losses in wheat production. The fungal pathogen is homothallic in nature, which means it can undergo sexual reproduction by selfing to produce pseudothecia on wheat stubble for seasonal survival. Since homothallism precludes the development of bi-parental fungal populations, no genetic linkage map has been developed for P. tritici-repentis for mapping and map-based cloning of fungal virulence genes. In this work, we created two heterothallic strains by deleting one of the mating type genes in each of two parental isolates 86-124 (race 2) and AR CrossB10 (a new race) and developed a bi-parental fungal population between them. The draft genome sequences of the two parental isolates were aligned to the Pt-1C-BFP reference sequence to mine single nucleotide polymorphisms (SNPs). A total of 225 SNP markers were developed for genotyping the entire population. Additionally, 75 simple sequence repeat, and two gene markers were also developed and used in the genotyping. The resulting linkage map consisted of 13 linkage groups spanning 5,075.83 cM in genetic distance. Because the parental isolate AR CrossB10 is a new race and produces Ptr ToxC, it was sequenced using long-read sequencing platforms and de novo assembled into contigs. The majority of the contigs were further anchored into chromosomes with the aid of the linkage maps. The whole genome comparison of AR CrossB10 to the reference genome of M4 revealed a few chromosomal rearrangements. The genetic linkage map and the new AR CrossB10 genome sequence are valuable tools for gene cloning in P. tritici-repentis.

A protein kinase-major sperm protein gene hijacked by a necrotrophic fungal pathogen triggers disease susceptibility in wheat

Septoria nodorum blotch (SNB), a disease caused by the necrotrophic fungal pathogen Parastagonospora nodorum, is a threat to wheat (Triticum aestivum) production worldwide. Multiple inverse gene-for-gene interactions involving the recognition of necrotrophic effectors (NEs) by wheat sensitivity genes play major roles in causing SNB. One interaction involves the wheat gene Snn3 and the P. nodorum NE SnTox3. Here, we used a map-based strategy to clone the Snn3-D1 gene from Aegilops tauschii, the D-genome progenitor of common wheat. Snn3-D1 contained protein kinase and major sperm protein domains, both of which were essential for function as confirmed by mutagenesis. As opposed to other characterized interactions in this pathosystem, a compatible Snn3-D1-SnTox3 interaction was light-independent, and Snn3-D1 transcriptional expression was downregulated by light and upregulated by darkness. Snn3-D1 likely emerged in Ae. tauschii due to an approximately 218-kb insertion that occurred along the west bank of the Caspian Sea. The identification of this new class of NE sensitivity genes combined with the previously cloned sensitivity genes demonstrates that P. nodorum can take advantage of diverse host targets to trigger SNB susceptibility in wheat.

Identification and characterization of Cercospora beticola necrosis-inducing effector CbNip1

Cercospora beticola is a hemibiotrophic fungus that causes cercospora leaf spot disease of sugar beet (Beta vulgaris). After an initial symptomless biotrophic phase of colonization, necrotic lesions appear on host leaves as the fungus switches to a necrotrophic lifestyle. The phytotoxic secondary metabolite cercosporin has been shown to facilitate fungal virulence for several Cercospora spp. However, because cercosporin production and subsequent cercosporin-initiated formation of reactive oxygen species is light-dependent, cell death evocation by this toxin is only fully ensured during a period of light. Here, we report the discovery of the effector protein CbNip1 secreted by C. beticola that causes enhanced necrosis in the absence of light and, therefore, may complement light-dependent necrosis formation by cercosporin. Infiltration of CbNip1 protein into sugar beet leaves revealed that darkness is essential for full CbNip1-triggered necrosis, as light exposure delayed CbNip1-triggered host cell death. Gene expression analysis during host infection shows that CbNip1 expression is correlated with symptom development in planta. Targeted gene replacement of CbNip1 leads to a significant reduction in virulence, indicating the importance of CbNip1 during colonization. Analysis of 89 C. beticola genomes revealed that CbNip1 resides in a region that recently underwent a selective sweep, suggesting selection pressure exists to maintain a beneficial variant of the gene. Taken together, CbNip1 is a crucial effector during the C. beticola-sugar beet disease process.

The Barley HvWRKY6 Transcription Factor Is Required for Resistance Against Pyrenophora teres f. teres

Barley is an important cereal crop worldwide because of its use in the brewing and distilling industry. However, adequate supplies of quality malting barley are threatened by global climate change due to drought in some regions and excess precipitation in others, which facilitates epidemics caused by fungal pathogens. The disease net form net blotch caused by the necrotrophic fungal pathogen Pyrenophora teres f. teres (Ptt) has emerged as a global threat to barley production and diverse populations of Ptt have shown a capacity to overcome deployed genetic resistances. The barley line CI5791 exhibits remarkably effective resistance to diverse Ptt isolates from around the world that maps to two major QTL on chromosomes 3H and 6H. To identify genes involved in this effective resistance, CI5791 seed were γ-irradiated and two mutants, designated CI5791-γ3 and CI5791-γ8, with compromised Ptt resistance were identified from an M₂ population. Phenotyping of CI5791-γ3 and -γ8 × Heartland F₂ populations showed three resistant to one susceptible segregation ratios and CI5791-γ3 × -γ8 F₁ individuals were susceptible, thus these independent mutants are in a single allelic gene. Thirty-four homozygous mutant (susceptible) CI5791-γ3 × Heartland F₂ individuals, representing 68 recombinant gametes, were genotyped via PCR genotype by sequencing. The data were used for single marker regression mapping placing the mutation on chromosome 3H within an approximate 75 cM interval encompassing the 3H CI5791 resistance QTL. Sequencing of the mutants and wild-type (WT) CI5791 genomic DNA following exome capture identified independent mutations of the HvWRKY6 transcription factor located on chromosome 3H at ∼50.7 cM, within the genetically delimited region. Post transcriptional gene silencing of HvWRKY6 in barley line CI5791 resulted in Ptt susceptibility, confirming that it functions in NFNB resistance, validating it as the gene underlying the mutant phenotypes. Allele analysis and transcript regulation of HvWRKY6 from resistant and susceptible lines revealed sequence identity and upregulation upon pathogen challenge in all genotypes analyzed, suggesting a conserved transcription factor is involved in the defense against the necrotrophic pathogen. We hypothesize that HvWRKY6 functions as a conserved signaling component of defense mechanisms that restricts Ptt growth in barley.

Four Reference Quality Genome Assemblies of Pyrenophora teres f. maculata: A Resource for Studying the Barley Spot Form Net Blotch Interaction

Pyrenophora teres is the causal agent of net blotch, the most devastating foliar disease of barley. In nature, net blotch is seen in two forms, net form net blotch, caused by P. teres f. teres, and spot form net blotch, caused by P. teres f. maculata. To date, 11 P. teres f. teres genomes have been sequenced and deposited in publicly available repositories, but only one P. teres f. maculata genome has been publicly deposited. Here, we present four additional reference-quality full-genome sequences of P. teres f. maculata isolates with good geographical and phenotypic diversity, with accompanying RNA sequencing-based genome annotations. These additional P. teres f. maculata genomes will aid in the understanding of the genomic complexities of this important barley pathogen.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

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.

Plant genes hijacked by necrotrophic fungal pathogens

Plant fungal pathogens can be classified according to their lifestyles. Biotrophs feed on living tissue and constitute an economically significant group of pathogens historically. Necrotrophs, which feed on dead tissue, have become economically significant over recent decades, especially those of the Dothideomycetes, which produce necrotrophic effectors (NEs) to modulate the host response. Some of these pathogens interact with their hosts in an inverse gene-for-gene manner, where NEs are recognized by specific dominant genes in the host leading to host-mediated programmed cell death allowing the pathogen to cause disease. Whereas the NE genes tend to be unique, several of the plant 'susceptibility' genes belong to the nucleotide-binding leucine-rich repeat class of disease 'resistance' genes, and one is a wall-associated kinase. These susceptible interactions exhibit hallmarks of defense responses to biotrophic pathogens. Therefore, there is now accumulating evidence that many necrotrophic specialists hijack the resistance mechanisms that are effective against biotrophic pathogens.

Research advances in the Pyrenophora teres-barley interaction

Pyrenophora teres f. teres and P. teres f. maculata are significant pathogens that cause net blotch of barley. An increased number of loci involved in P. teres resistance or susceptibility responses of barley as well as interacting P. teres virulence effector loci have recently been identified through biparental and association mapping studies of both the pathogen and host. Characterization of the resistance/susceptibility loci in the host and the interacting effector loci in the pathogen will provide a path for targeted gene validation for better-informed release of resistant barley cultivars. This review assembles concise consensus maps for all loci published for both the host and pathogen, providing a useful resource for the community to be used in pathogen characterization and barley breeding for resistance to both forms of P. teres.

A Comparative Genomic Analysis of the Barley Pathogen Pyrenophora teres f. teres Identifies Subtelomeric Regions as Drivers of Virulence

Pyrenophora teres f. teres causes net form net blotch of barley and is an economically important pathogen throughout the world. However, P. teres f. teres is lacking in the genomic resources necessary to characterize the mechanisms of virulence. Recently a high-quality reference genome was generated for P. teres f. teres isolate 0-1. Here, we present the reference quality sequence and annotation of four new isolates and we use the five available P. teres f. teres genomes for an in-depth comparison, resulting in the generation of hypotheses pertaining to the potential mechanisms and evolution of virulence. Comparative analyses were performed between all five P. teres f. teres genomes, examining genomic organization, structural variations, and core and accessory genomic content, specifically focusing on the genomic characterization of known virulence loci and the localization of genes predicted to encode secreted and effector proteins. We showed that 14 of 15 currently published virulence quantitative trait loci (QTL) span accessory genomic regions, consistent with these accessory regions being important drivers of host adaptation. Additionally, these accessory genomic regions were frequently found in subtelomeric regions of chromosomes, with 10 of the 14 accessory region QTL localizing to subtelomeric regions. Comparative analysis of the subtelomeric regions of P. teres f. teres chromosomes revealed translocation events in which homology was detected between nonhomologous chromosomes at a significantly higher rate than the rest of the genome. These results indicate that the subtelomeric accessory genomic compartments not only harbor most of the known virulence loci but, also, that these regions have the capacity to rapidly evolve.

Mapping and characterization of two stem rust resistance genes derived from cultivated emmer wheat accession PI 193883

Two stem rust resistance genes identified on chromosome arms 2BL and 6AL of the cultivated emmer wheat accession PI 193883 can be used for protecting modern varieties against Ug99 strains. The wheat research community consistently strives to identify new genes that confer resistance to stem rust caused by the fungal pathogen Puccinia graminis f. sp. tritici Eriks & E. Henn (Pgt). In the current study, our objective was to identify and genetically characterize the stem rust resistance derived from the cultivated emmer accession PI 193883. A recombinant inbred line population developed from a cross between the susceptible durum wheat line Rusty and PI 193883 was genotyped and evaluated for reaction to Pgt races TTKSK, TRTTF, and TMLKC. Two QTLs conferring resistance were identified on chromosome arms 2BL (QSr.fcu-2B) and 6AL (QSr.fcu-6A). The stem rust resistance gene (Sr883-2B) underlying QSr.fcu-2B was recessive, and based on its physical location it is located proximal to the Sr9 region. QSr.fcu-6A was located in the Sr13 region, but PI 193883 is known to carry the susceptible haplotype S4 for Sr13, indicating that the gene underlying QSr.fcu-6A (Sr883-6A) is likely a new allele of Sr13 or a gene residing close to Sr13. Three IWGSC scaffold-based simple sequence repeat (SSR) and two SNP-based semi-thermal asymmetric reverse PCR (STARP) markers were developed for the Sr883-2B region, and one STARP marker was developed for Sr883-6A. Sr883-2B was epistatic to Sr883-6A for reaction to TTKSK and TRTTF, and the two genes had additive effects for TMLKC. These two genes and the markers developed in this research provide additional resources and tools for the improvement in stem rust resistance in durum and common wheat breeding programs.

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.

Mapping of barley susceptibility/resistance QTL against spot form net blotch caused by Pyrenophora teres f. maculata using RIL populations

Spot form net blotch (SFNB) caused by the necrotrophic fungal pathogen Pyrenophora teres f. maculata (Ptm) is an important disease of barley worldwide including the major barley production regions of North America. To characterize SFNB resistance/susceptibility quantitative trait loci (QTL), three recombinant inbred line (RIL) populations were developed from crosses between the malting barley cultivars, Tradition (six row) and Pinnacle (two row), and the two world barley core collection lines, PI67381 and PI84314. Tradition and Pinnacle were susceptible to many North American Ptm isolates, while PI67381 and PI84314 carry resistances to diverse Ptm isolates from across the globe. The RIL populations, Tradition/PI67381, Pinnacle/PI67381, and Pinnacle/PI84314 were genotyped using polymerase chain reaction-mediated genotype-by-sequencing single nucleotide polymorphism marker panels and phenotyped at the seedling stage with six geographically distinct Ptm isolates: FGOB10Ptm-1 (North Dakota, USA), Pin-A14 (Montana, USA), Cel-A17 (Montana, USA), SG1 (Australia), NZKF2 (New Zealand) and DEN2.6 (Denmark). The goal was to determine if the susceptible elite lines contained common susceptibility genes/QTL or if the resistant lines had common resistant genes/QTL effective against diverse Ptm isolates. The QTL analyses identified a total of 12 resistance and/or susceptibility loci on chromosomes 2H, 3H, 4H, 6H, and 7H of which three had not been previously reported. Common major QTL were detected on chromosome 2H (R² = 14-40%) and 7H (R² = 24-80%) in all three RIL populations, suggesting underlying genes with broad resistance specificity. The major 7H QTL was shown to be a dominant susceptibility gene in both susceptible malting barley varieties.

Genetics of Variable Disease Expression Conferred by Inverse Gene-For-Gene Interactions in the Wheat-Parastagonospora nodorum Pathosystem

The wheat-Parastagonospora nodorum pathosystem involves the recognition of pathogen-secreted necrotrophic effectors (NEs) by corresponding wheat NE sensitivity genes. This inverse gene-for-gene recognition leads to necrotrophic effector-triggered susceptibility and ultimately septoria nodorum blotch disease. Here, we used multiple pathogen isolates to individually evaluate the effects of the host gene-NE interactions Tan spot necrosis1-Stagonospora nodorum ToxinA (Tsn1-SnToxA), Stagonospora nodorum necrosis1-Stagonospora nodorum Toxin1 (Snn1-SnTox1), and Stagonospora nodorum necrosis3-B genome homeolog1-Stagonospora nodorum Toxin3 (Snn3-B1-SnTox3), alone and in various combinations, to determine the relative importance of these interactions in causing disease. Genetic analysis of a recombinant inbred wheat population inoculated separately with three P. nodorum isolates, all of which produce all three NEs, indicated that the Tsn1-SnToxA and Snn3-B1-SnTox3 interactions contributed to disease caused by all four isolates, but their effects varied and ranged from epistatic to additive. The Snn1-SnTox1 interaction was associated with increased disease for one isolate, but for other isolates, there was evidence that this interaction inhibited the expression of other host gene-NE interactions. RNA sequencing analysis in planta showed that SnTox1 was differentially expressed between these three isolates after infection. Further analysis of NE gene-knockout isolates showed that the effect of some interactions could be masked or inhibited by other compatible interactions, and the regulation of this occurs at the level of NE gene transcription. Collectively, these results show that the inverse gene-for-gene interactions leading to necrotrophic effector-triggered susceptibility in the wheat-P. nodorum pathosystem vary in their effects depending on the genetic backgrounds of the pathogen and host, and interplay among the interactions is complex and intricately regulated.

Genetic Diversity and Resistance to Fusarium Head Blight in Synthetic Hexaploid Wheat Derived From Aegilops tauschii and Diverse Triticum turgidum Subspecies

Synthetic hexaploid wheat (SHW) can serve as a bridge for the transfer of useful genes from Aegilops tauschii and tetraploid wheat (Triticum turgidum) into common wheat (T. aestivum). The objective of this study was to evaluate 149 SHW lines and their 74 tetraploid parents for their genetic diversity, breeding values and inter-genomic interactions for resistance to Fusarium head blight (FHB). The genetic diversity analysis was performed based on the population structure established using 4,674 and 3,330 polymorphic SNP markers among the SHW lines and tetraploid parents, respectively. The results showed that all T. carthlicum and most T. dicoccum accessions formed different clusters and subpopulations, respectively, whereas all the T. durum, T. polonicum, T. turgidum, and T. turanicum accessions were clustered together, suggesting that T. durum was more closely related to T. polonicum, T. turgidum, and T. turanicum than to T. dicoccum. The genetic diversity of the SHW lines mainly reflected that of the tetraploid parents. The SHW lines and their tetraploid parents were evaluated for reactions to FHB in two greenhouse seasons and at two field nurseries for 2 years. As expected, most of the SHW lines were more resistant than their tetraploid parents in all environments. The FHB severities of the SHW lines varied greatly depending on the Ae. tauschii and tetraploid genotypes involved. Most of the SHW lines with a high level of FHB resistance were generally derived from the tetraploid accessions with a high level of FHB resistance. Among the 149 SHW lines, 140 were developed by using three Ae. tauschii accessions CIae 26, PI 268210, and RL 5286. These SHW lines showed FHB severities reduced by 21.7%, 17.3%, and 11.5%, respectively, with an average reduction of 18.3%, as compared to the tetraploid parents, suggesting that the D genome may play a major role in reducing disease severity in the SHW lines. Thirteen SHW lines consistently showed a high level of FHB resistance compared to the resistant check, Sumai 3, in each environment. These SHW lines will be useful for the development of FHB-resistant wheat germplasm and populations for discovery of novel FHB resistance genes.

ToxA Is Present in the U.S. Bipolaris sorokiniana Population and Is a Significant Virulence Factor on Wheat Harboring Tsn1

ToxA, a necrotrophic effector originally identified from the tan spot fungus Pyrenophora tritici-repentis in 1987, was subsequently identified from Parastagonospora nodorum in 2006. More recently, the ToxA gene was identified in the spot blotch fungus Bipolaris sorokiniana in Australia. Here we show that the ToxA gene is also present in the B. sorokiniana population in the winter wheat region of southcentral Texas. Leaves from 'Duster' wheat showing strong necrotic lesions were collected in Castroville, TX. Fifteen single-spore isolates were collected from separate lesions, and 13 of them harbored the BsToxA gene and secreted ToxA in culture based on sensitivity of BG261, the differential line containing the dominant ToxA sensitivity gene, Tsn1. Four isolates harboring BsToxA and one deficient in BsToxA were used to infiltrate two wheat lines harboring Tsn1 as well as their corresponding tsn1 mutant lines. Culture filtrates of the isolate lacking BsToxA did not induce necrosis on any of the lines. Culture filtrates of the four BsToxA-containing isolates induced necrosis on the wild type (Tsn1) lines but not on the corresponding tsn1 mutant lines. Sensitivity to these culture filtrates also mapped to the previously identified location for Tsn1 in the winter wheat mapping population Arina × Forno. Inoculation of one of these ToxA-producing isolates on the same population showed that the Tsn1 locus accounted for 24.4% of the disease variation. All 13 isolates harbored the same BsToxA nucleotide sequence, which was identical to one of the two haplotypes previously identified in Australia. Sensitivity to ToxA is prevalent in popular hard winter wheat cultivars in the central and southcentral winter wheat regions of the United States, showing the potential of a selective advantage for B. sorokiniana isolates that harbor the ToxA gene.

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.

Transposable Element Genomic Fissuring in Pyrenophora teres Is Associated With Genome Expansion and Dynamics of Host-Pathogen Genetic Interactions

Pyrenophora teres, P. teres f. teres (PTT) and P. teres f. maculata (PTM) cause significant diseases in barley, but little is known about the large-scale genomic differences that may distinguish the two forms. Comprehensive genome assemblies were constructed from long DNA reads, optical and genetic maps. As repeat masking in fungal genomes influences the final gene annotations, an accurate and reproducible pipeline was developed to ensure comparability between isolates. The genomes of the two forms are highly collinear, each composed of 12 chromosomes. Genome evolution in P. teres is characterized by genome fissuring through the insertion and expansion of transposable elements (TEs), a process that isolates blocks of genic sequence. The phenomenon is particularly pronounced in PTT, which has a larger, more repetitive genome than PTM and more recent transposon activity measured by the frequency and size of genome fissures. PTT has a longer cultivated host association and, notably, a greater range of host-pathogen genetic interactions compared to other Pyrenophora spp., a property which associates better with genome size than pathogen lifestyle. The two forms possess similar complements of TE families with Tc1/Mariner and LINE-like Tad-1 elements more abundant in PTT. Tad-1 was only detectable as vestigial fragments in PTM and, within the forms, differences in genome sizes and the presence and absence of several TE families indicated recent lineage invasions. Gene differences between P. teres forms are mainly associated with gene-sparse regions near or within TE-rich regions, with many genes possessing characteristics of fungal effectors. Instances of gene interruption by transposons resulting in pseudogenization were detected in PTT. In addition, both forms have a large complement of secondary metabolite gene clusters indicating significant capacity to produce an array of different molecules. This study provides genomic resources for functional genetics to help dissect factors underlying the host-pathogen interactions.

Reference Quality Genome Assemblies of Three Parastagonospora nodorum Isolates Differing in Virulence on Wheat

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.

Characterizing the Pyrenophora teres f. maculata-Barley Interaction Using Pathogen Genetics

Pyrenophora teres f. maculata is the cause of the foliar disease spot form net blotch (SFNB) on barley. To evaluate pathogen genetics underlying the P. teres f. maculata-barley interaction, we developed a 105-progeny population by crossing two globally diverse isolates, one from North Dakota and the other from Western Australia. Progeny were phenotyped on a set of four barley genotypes showing a differential reaction to the parental isolates, then genotyped using a restriction site-associated-genotype-by-sequencing (RAD-GBS) approach. Genetic maps were developed for use in quantitative trait locus (QTL) analysis to identify virulence-associated QTL. Six QTL were identified on five different linkage groups and individually accounted for 20-37% of the disease variation, with the number of significant QTL ranging from two to four for the barley genotypes evaluated. The data presented demonstrate the complexity of virulence involved in the P. teres f. maculata-barley pathosystem and begins to lay the foundation for understanding this important interaction.

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.

Association mapping utilizing diverse barley lines reveals net form net blotch seedling resistance/susceptibility loci

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.

Genetic analysis of net form net blotch resistance in barley lines CIho 5791 and Tifang against a global collection of P. teres f. teres isolates

A CIho 5791 × Tifang recombinant inbred mapping population was developed and used to identify major dominant resistance genes on barley chromosomes 6H and 3H in CI5791 and on 3H in Tifang. The barley line CIho 5791 confers high levels of resistance to Pyrenophora teres f. teres, causal agent of net form net blotch (NFNB), with few documented isolates overcoming this resistance. Tifang barley also harbors resistance to P. teres f. teres which was previously shown to localize to barley chromosome 3H. A CIho 5791 × Tifang F₆ recombinant inbred line (RIL) population was developed using single seed descent. The Illumina iSelect SNP platform was used to identify 2562 single nucleotide polymorphism (SNP) markers across the barley genome, resulting in seven linkage maps, one for each barley chromosome. The CIho 5791 × Tifang RIL population was evaluated for NFNB resistance using nine P. teres f. teres isolates collected globally. Tifang was resistant to four of the isolates tested whereas CIho 5791 was highly resistant to all nine isolates. QTL analysis indicated that the CIho 5791 resistance mapped to chromosome 6H whereas the Tifang resistance mapped to chromosome 3H. Additionally, CIho 5791 also harbored resistance to two Japanese isolates that mapped to a 3H region similar to that of Tifang. SNP markers and RILs harboring both 3H and 6H resistance will be useful in resistance breeding against NFNB.

Nicotiana benthamiana as a nonhost of Zymoseptoria tritici

This article is a Commentary on Kettles et al., 213: 338–350.

New Insights into the Roles of Host Gene-Necrotrophic Effector Interactions in Governing Susceptibility of Durum Wheat to Tan Spot and Septoria nodorum Blotch

Tan spot and Septoria nodorum blotch (SNB) are important diseases of wheat caused by the necrotrophic fungi Pyrenophora tritici-repentis and Parastagonospora nodorum, respectively. The P. tritici-repentis necrotrophic effector (NE) Ptr ToxB causes tan spot when recognized by the Tsc2 gene. The NE ToxA is produced by both pathogens and has been associated with the development of both tan spot and SNB when recognized by the wheat Tsn1 gene. Most work to study these interactions has been conducted in common wheat, but little has been done in durum wheat. Here, quantitative trait loci (QTL) analysis of a segregating biparental population indicated that the Tsc2-Ptr ToxB interaction plays a prominent role in the development of tan spot in durum. However, analysis of two biparental populations indicated that the Tsn1-ToxA interaction was not associated with the development of tan spot, but was strongly associated with the development of SNB. Pa. nodorum expressed ToxA at high levels in infected Tsn1 plants, whereas ToxA expression in P. tritici-repentis was barely detectable, suggesting that the differences in disease levels associated with the Tsn1-ToxA interaction were due to differences in pathogen expression of ToxA These and previous results together indicate that: (1) the effects of Tsn1-ToxA on tan spot in common wheat can range from nonsignificant to highly significant depending on the host genetic background; (2) Tsn1-ToxA is not a significant factor for tan spot development in durum wheat; and (3) Tsn1-ToxA plays a major role in SNB development in both common and durum wheat. Durum and common wheat breeders alike should strive to remove both Tsc2 and Tsn1 from their materials to achieve disease resistance.

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.

The hijacking of a receptor kinase-driven pathway by a wheat fungal pathogen leads to disease

Necrotrophic pathogens live and feed on dying tissue, but their interactions with plants are not well understood compared to biotrophic pathogens. The wheat Snn1 gene confers susceptibility to strains of the necrotrophic pathogen Parastagonospora nodorum that produce the SnTox1 protein. We report the positional cloning of Snn1, a member of the wall-associated kinase class of receptors, which are known to drive pathways for biotrophic pathogen resistance. Recognition of SnTox1 by Snn1 activates programmed cell death, which allows this necrotroph to gain nutrients and sporulate. These results demonstrate that necrotrophic pathogens such as P. nodorum hijack host molecular pathways that are typically involved in resistance to biotrophic pathogens, revealing the complex nature of susceptibility and resistance in necrotrophic and biotrophic pathogen interactions with plants.