De Novo Assembly and Phasing of Dikaryotic Genomes from Two Isolates of Puccinia coronata f. sp. avenae, the Causal Agent of Oat Crown Rust

Intrinsically Disordered Kiwellin Protein-Like Effectors Target Plant Chloroplasts and are Extensively Present in Rust Fungi

The effector proteins produced by plant pathogens are one of the essential components of host-pathogen interaction. Despite being important, most of the effector proteins remain unexplored due to the diversity in their primary sequence generated by the high selection pressure of the host immune system. However to maintain the primary function in the infection process, these effectors may tend to maintain their native protein fold to perform the corresponding biological function. In the present study, unannotated candidate secretory effector proteins of sixteen major plant fungal pathogens were analyzed to find the conserved known protein folds using homology, ab initio, and Alpha Fold/Rosetta Fold protein dimensional (3D) structure approaches. Several unannotated candidate effector proteins were found to match various known conserved protein families potentially involved in host defense manipulation in different plant pathogens. Surprisingly a large number of plant Kiwellin proteins fold like secretory proteins (> 100) were found in studied rust fungal pathogens. Many of them were predicted as potential effector proteins. Furthermore, template independent modelling using Alpha Fold/Rosetta Fold analysis and structural comparison of these candidates also predicted them to match with plant Kiwellin proteins. We also found plant Kiwellin matching proteins outside rusts including several non-pathogenic fungi suggesting the broad function of these proteins. One of the highest confidently modeled Kiwellin matching candidates effectors, Pstr_13960 (97.8%), from the Indian P. striiformis race Yr9 was characterized using overexpression, localization, and deletion studies in Nicotiana benthamiana. The Pstr_13960 suppressed the BAX-induced cell death and localized in the chloroplast. Furthermore, the expression of the Kiwellin matching region (Pst_13960_kiwi) alone suppressed the BAX-induced cell death in N. benthamiana despite the change of location to the cytoplasm and nucleus, suggesting the novel function of the Kiwellin core fold in rust fungi. Molecular docking showed that Pstr_13960 can interact with plant Chorismate mutases (CMs) using three loops conserved in plant and rust Kiwellins. Further analysis of Pstr_13960 showed to contain Intrinsically disordered regions (IDRs) in place of the N-terminal β1/β2 region found in plant Kiwellins suggesting the evolution of rust Kiwellins-like effectors (KLEs). Overall, this study reports the presence of a Kiwellin protein-like fold containing a novel effector protein family in rust fungi depicting a classical example of the evolution of effectors at the structure level as Kiwellin effectors show very low significant similarity to plant Kiwellin at the sequence level.

Genome biology and evolution of mating-type loci in four cereal rust fungi

Permanent heterozygous loci, such as sex- or mating-compatibility regions, often display suppression of recombination and signals of genomic degeneration. In Basidiomycota, two distinct loci confer mating compatibility. These loci encode homeodomain (HD) transcription factors and pheromone receptor (Pra)-ligand allele pairs. To date, an analysis of genome level mating-type (MAT) loci is lacking for obligate biotrophic basidiomycetes in the Pucciniales, an order containing serious agricultural plant pathogens. Here, we focus on four species of Puccinia that infect oat and wheat, including P. coronata f. sp. avenae, P. graminis f. sp. tritici, P. triticina and P. striiformis f. sp. tritici. MAT loci are located on two separate chromosomes supporting previous hypotheses of a tetrapolar mating compatibility system in the Pucciniales. The HD genes are multiallelic in all four species while the PR locus appears biallelic, except for P. graminis f. sp. tritici, which potentially has multiple alleles. HD loci are largely conserved in their macrosynteny, both within and between species, without strong signals of recombination suppression. Regions proximal to the PR locus, however, displayed signs of recombination suppression and genomic degeneration in the three species with a biallelic PR locus. Our observations support a link between recombination suppression, genomic degeneration, and allele diversity of MAT loci that is consistent with recent mathematical modelling and simulations. Finally, we confirm that MAT genes are expressed during the asexual infection cycle, and we propose that this may support regulating nuclear maintenance and pairing during infection and spore formation. Our study provides insights into the evolution of MAT loci of key pathogenic Puccinia species. Understanding mating compatibility can help predict possible combinations of nuclear pairs, generated by sexual reproduction or somatic recombination, and the potential evolution of new virulent isolates of these important plant pathogens.

Genome-Enabled Analysis of Population Dynamics and Virulence-Associated Loci in the Oat Crown Rust Fungus Puccinia coronata f. sp. avenae

Puccinia coronata f. sp. avenae (Pca) is an important fungal pathogen causing crown rust that impacts oat production worldwide. Genetic resistance for crop protection against Pca is often overcome by the rapid virulence evolution of the pathogen. This study investigated the factors shaping adaptive evolution of Pca using pathogen populations from distinct geographic regions within the United States and South Africa. Phenotypic and genome-wide sequencing data of these diverse Pca collections, including 217 isolates, uncovered phylogenetic relationships and established distinct genetic composition between populations from northern and southern regions from the United States and South Africa. The population dynamics of Pca involve a bidirectional movement of inoculum between northern and southern regions of the United States and contributions from clonality and sexuality. The population from South Africa is solely clonal. A genome-wide association study (GWAS) employing a haplotype-resolved Pca reference genome was used to define 11 virulence-associated loci corresponding to 25 oat differential lines. These regions were screened to determine candidate Avr effector genes. Overall, the GWAS results allowed us to identify the underlying genetic factors controlling pathogen recognition in an oat differential set used in the United States to assign pathogen races (pathotypes). Key GWAS findings support complex genetic interactions in several oat lines, suggesting allelism among resistance genes or redundancy of genes included in the differential set, multiple resistance genes recognizing genetically linked Avr effector genes, or potentially epistatic relationships. A careful evaluation of the composition of the oat differential set accompanied by the development or implementation of molecular markers is recommended. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Soybean-Phakopsora pachyrhizi interactions: towards the development of next-generation disease-resistant plants

Soybean rust (SBR), caused by the obligate biotrophic fungus Phakopsora pachyrhizi, is a devastating foliar disease threatening soybean production. To date, no commercial cultivars conferring durable resistance to SBR are available. The development of long-lasting SBR resistance has been hindered by the lack of understanding of this complex pathosystem, encompassing challenges posed by intricate genetic structures in both the host and pathogen, leading to a gap in the knowledge of gene-for-gene interactions between soybean and P. pachyrhizi. In this review, we focus on recent advancements and emerging technologies that can be used to improve our understanding of the P. pachyrhizi-soybean molecular interactions. We further explore approaches used to combat SBR, including conventional breeding, transgenic approaches and RNA interference, and how advances in our understanding of plant immune networks, the availability of new molecular tools, and the recent sequencing of the P. pachyrhizi genome could be used to aid in the development of better genetic resistance against SBR. Lastly, we discuss the research gaps of this pathosystem and how new technologies can be used to shed light on these questions and to develop durable next-generation SBR-resistant soybean plants.

Plant and pathogen genomics: essential approaches for stem rust resistance gene stacks in wheat

The deployment of disease resistance genes is currently the most economical and environmentally sustainable method of crop protection. However, disease resistance genes can rapidly break down because of constant pathogen evolution, particularly when they are deployed singularly. Polygenic resistance is, therefore, considered the most durable, but combining and maintaining these genes by breeding is a laborious process as effective genes are usually unlinked. The deployment of polygenic resistance with single-locus inheritance is a promising innovation that overcomes these difficulties while enhancing resistance durability. Because of major advances in genomic technologies, increasing numbers of plant resistance genes have been cloned, enabling the development of resistance transgene stacks (RTGSs) that encode multiple genes all located at a single genetic locus. Gene stacks encoding five stem rust resistance genes have now been developed in transgenic wheat and offer both breeding simplicity and potential resistance durability. The development of similar genomic resources in phytopathogens has advanced effector gene isolation and, in some instances, enabled functional validation of individual resistance genes in RTGS. Here, the wheat stem rust pathosystem is used as an illustrative example of how host and pathogen genomic advances have been instrumental in the development of RTGS, which is a strategy applicable to many other agricultural crop species.

Genetic Analysis and Physical Mapping of Oat Adult Plant Resistance Loci Against Puccinia coronata f. sp. avenae

Six quantitative trait loci (QTLs) for adult plant resistance against oat crown rust (Puccinia coronata f. sp. avenae) were identified from mapping three recombinant inbred populations. Using genotyping-by-sequencing with markers called against the OT3098 v1 reference genome, the QTLs were mapped on six different chromosomes: Chr1D, Chr4D, Chr5A, Chr5D, Chr7A, and Chr7C. Composite interval mapping with marker cofactor selection showed that the phenotypic variance explained by all identified QTLs for coefficient of infection range from 12.2 to 46.9%, whereas heritability estimates ranged from 0.11 to 0.38. The significant regions were narrowed down to intervals of 3.9 to 25 cM, equivalent to physical distances of 11 to 133 Mb. At least two flanking single-nucleotide polymorphism markers were identified within 10 cM of each QTL that could be used in marker-assisted introgression, pyramiding, and selection. The additive effects of the QTLs in each population were determined using single-nucleotide polymorphism haplotype data, which showed a significantly lower coefficient of infection in lines homozygous for the resistant alleles. Analysis of pairwise linkage disequilibrium also revealed high correlation of markers and presence of linkage blocks in the significant regions. To further facilitate marker-assisted breeding, polymerase chain reaction allelic competitive extension (PACE) markers for the adult plant resistance loci were developed. Putative candidate genes were also identified in each of the significant regions, which include resistance gene analogs that encode for kinases, ligases, and predicted receptors of avirulence proteins from pathogens.

Major proliferation of transposable elements shaped the genome of the soybean rust pathogen Phakopsora pachyrhizi

With >7000 species the order of rust fungi has a disproportionately large impact on agriculture, horticulture, forestry and foreign ecosystems. The infectious spores are typically dikaryotic, a feature unique to fungi in which two haploid nuclei reside in the same cell. A key example is Phakopsora pachyrhizi, the causal agent of Asian soybean rust disease, one of the world's most economically damaging agricultural diseases. Despite P. pachyrhizi's impact, the exceptional size and complexity of its genome prevented generation of an accurate genome assembly. Here, we sequence three independent P. pachyrhizi genomes and uncover a genome up to 1.25 Gb comprising two haplotypes with a transposable element (TE) content of ~93%. We study the incursion and dominant impact of these TEs on the genome and show how they have a key impact on various processes such as host range adaptation, stress responses and genetic plasticity.

Exogenous double-stranded RNA inhibits the infection physiology of rust fungi to reduce symptoms in planta

Rust fungi (Pucciniales) are a diverse group of plant pathogens in natural and agricultural systems. They pose ongoing threats to the diversity of native flora and cause annual crop yield losses. Agricultural rusts are predominantly managed with fungicides and breeding for resistance, but new control strategies are needed on non-agricultural plants and in fragile ecosystems. RNA interference (RNAi) induced by exogenous double-stranded RNA (dsRNA) has promise as a sustainable approach for managing plant-pathogenic fungi, including rust fungi. We investigated the mechanisms and impact of exogenous dsRNA on rust fungi through in vitro and whole-plant assays using two species as models, Austropuccinia psidii (the cause of myrtle rust) and Coleosporium plumeriae (the cause of frangipani rust). In vitro, dsRNA either associates externally or is internalized by urediniospores during the early stages of germination. The impact of dsRNA on rust infection architecture was examined on artificial leaf surfaces. dsRNA targeting predicted essential genes significantly reduced germination and inhibited development of infection structures, namely appressoria and penetration pegs. Exogenous dsRNA sprayed onto 1-year-old trees significantly reduced myrtle rust symptoms. Furthermore, we used comparative genomics to assess the wide-scale amenability of dsRNA to control rust fungi. We sequenced genomes of six species of rust fungi, including three new families (Araucariomyceaceae, Phragmidiaceae, and Skierkaceae) and identified key genes of the RNAi pathway across 15 species in eight families of Pucciniales. Together, these findings indicate that dsRNA targeting essential genes has potential for broad-use management of rust fungi across natural and agricultural systems.

Gapless Genome Assembly of Puccinia triticina Provides Insights into Chromosome Evolution in Pucciniales

Chromosome evolution drives species evolution, speciation, and adaptive radiation. Accurate genome assembly is crucial to understanding chromosome evolution of species, such as dikaryotic fungi. Rust fungi (Pucciniales) in dikaryons represent the largest group of plant pathogens, but the evolutionary process of adaptive radiation in Pucciniales remains poorly understood. Here, we report a gapless genome for the wheat leaf rust fungus Puccinia triticina determined using PacBio high-fidelity (HiFi) sequencing. This gapless assembly contains two sets of chromosomes, showing that one contig represents one chromosome. Comparisons of homologous chromosomes between the phased haplotypes revealed that highly frequent small-scale sequence divergence shapes haplotypic variation. Genome analyses of Puccinia triticina along with other rusts revealed that recent transposable element bursts and extensive segmental gene duplications synergistically highlight the evolution of chromosome structures. Comparative analysis of chromosomes indicated that frequent chromosomal rearrangements may act as a major contributor to rapid radiation of Pucciniales. This study presents the first gapless, phased assembly for a dikaryotic rust fungus and provides insights into adaptive evolution and species radiation in Pucciniales. IMPORTANCE Rust fungi (Pucciniales) are the largest group of plant pathogens. Adaptive radiation is a predominant feature in Pucciniales evolution. Chromosome evolution plays an important role in adaptive evolution. Accurate chromosome-scale assembly is required to understand the role of chromosome evolution in Pucciniales. We took advantage of HiFi sequencing to construct a gapless, phased genome for Puccinia triticina. Further analyses revealed that the evolution of chromosome structures in rust lineage is shaped by the combination of transposable element bursts and segmental gene duplications. Chromosome comparisons of Puccinia triticina and other rusts suggested that frequent chromosomal arrangements may make remarkable contributions to high species diversity of rust fungi. Our results present the first gapless genome for Pucciniales and shed light on the feature of chromosome evolution in Pucciniales.

Whole Genome Sequencing and Comparative Genomics of Indian Isolates of Wheat Spot Blotch Pathogen Bipolaris sorokiniana Reveals Expansion of Pathogenicity Gene Clusters

Spot blotch is a highly destructive disease in wheat caused by the fungal pathogen Bipolaris sorokiniana (teleomorph, Cochliobolus sativus). It is prevalent in warm and humid areas, including Africa, Asia, Latin America, and the USA. In the present study, twelve isolates of B. sorokiniana were collected from wheat fields in three different geographical locations in India. The pathogenicity of seven sporulating isolates was assessed on 'DDK 1025', a spot blotch-susceptible wheat variety under greenhouse conditions. The isolate 'D2' illustrated the highest virulence, followed by 'SI' and 'BS52'. These three isolates were sequenced using the Illumina HiSeq1000 platform. The estimated genome sizes of the isolates BS52, D2, and SI were 35.19 MB, 39.32 MB, and 32.76 MB, with GC contents of 48.48%, 50.43%, and 49.42%, respectively. The numbers of pathogenicity genes identified in BS52, D2, and SI isolates were 2015, 2476, and 2018, respectively. Notably, the isolate D2 exhibited a relatively larger genome with expanded arsenals of Biosynthetic Gene Clusters (BGCs), CAZymes, secretome, and pathogenicity genes, which could have contributed to its higher virulence among the tested isolates. This study provides the first comparative genome analysis of the Indian isolates of B. sorokiniana using whole genome sequencing.

Population genomics of Puccinia graminis f.sp. tritici highlights the role of admixture in the origin of virulent wheat rust races

Puccinia graminis f.sp. tritici (Pgt) causes stem rust disease in wheat that can result in severe yield losses. The factors driving the evolution of its virulence and adaptation remain poorly characterized. We utilize long-read sequencing to develop a haplotype-resolved genome assembly of a U.S. isolate of Pgt. Using Pgt haplotypes as a reference, we characterize the structural variants (SVs) and single nucleotide polymorphisms in a diverse panel of isolates. SVs impact the repertoire of predicted effectors, secreted proteins involved in host-pathogen interaction, and show evidence of purifying selection. By analyzing global and local genomic ancestry we demonstrate that the origin of 8 out of 12 Pgt clades is linked with either somatic hybridization or sexual recombination between the diverged donor populations. Our study shows that SVs and admixture events appear to play an important role in broadening Pgt virulence and the origin of highly virulent races, creating a resource for studying the evolution of Pgt virulence and preventing future epidemic outbreaks.

Current Status and Future Perspectives of Genomics Research in the Rust Fungi

Rust fungi in Pucciniales have caused destructive plant epidemics, have become more aggressive with new virulence, rapidly adapt to new environments, and continually threaten global agriculture. With the rapid advancement of genome sequencing technologies and data analysis tools, genomics research on many of the devastating rust fungi has generated unprecedented insights into various aspects of rust biology. In this review, we first present a summary of the main findings in the genomics of rust fungi related to variations in genome size and gene composition between and within species. Then we show how the genomics of rust fungi has promoted our understanding of the pathogen virulence and population dynamics. Even with great progress, many questions still need to be answered. Therefore, we introduce important perspectives with emphasis on the genome evolution and host adaptation of rust fungi. We believe that the comparative genomics and population genomics of rust fungi will provide a further understanding of the rapid evolution of virulence and will contribute to monitoring the population dynamics for disease management.

A chromosome-level, fully phased genome assembly of the oat crown rust fungus Puccinia coronata f. sp. avenae: a resource to enable comparative genomics in the cereal rusts

Advances in sequencing technologies as well as development of algorithms and workflows have made it possible to generate fully phased genome references for organisms with nonhaploid genomes such as dikaryotic rust fungi. To enable discovery of pathogen effectors and further our understanding of virulence evolution, we generated a chromosome-scale assembly for each of the 2 nuclear genomes of the oat crown rust pathogen, Puccinia coronata f. sp. avenae (Pca). This resource complements 2 previously released partially phased genome references of Pca, which display virulence traits absent in the isolate of historic race 203 (isolate Pca203) which was selected for this genome project. A fully phased, chromosome-level reference for Pca203 was generated using PacBio reads and Hi-C data and a recently developed pipeline named NuclearPhaser for phase assignment of contigs and phase switch correction. With 18 chromosomes in each haplotype and a total size of 208.10 Mbp, Pca203 has the same number of chromosomes as other cereal rust fungi such as Puccinia graminis f. sp. tritici and Puccinia triticina, the causal agents of wheat stem rust and wheat leaf rust, respectively. The Pca203 reference marks the third fully phased chromosome-level assembly of a cereal rust to date. Here, we demonstrate that the chromosomes of these 3 Puccinia species are syntenous and that chromosomal size variations are primarily due to differences in repeat element content.

Infection Strategies and Pathogenicity of Biotrophic Plant Fungal Pathogens

Biotrophic plant pathogenic fungi are widely distributed and are among the most damaging pathogenic organisms of agriculturally important crops responsible for significant losses in quality and yield. However, the pathogenesis of obligate parasitic pathogenic microorganisms is still under investigation because they cannot reproduce and complete their life cycle on an artificial medium. The successful lifestyle of biotrophic fungal pathogens depends on their ability to secrete effector proteins to manipulate or evade plant defense response. By integrating genomics, transcriptomics, and effectoromics, insights into how the adaptation of biotrophic plant fungal pathogens adapt to their host populations can be gained. Efficient tools to decipher the precise molecular mechanisms of rust–plant interactions, and standardized routines in genomics and functional pipelines have been established and will pave the way for comparative studies. Deciphering fungal pathogenesis not only allows us to better understand how fungal pathogens infect host plants but also provides valuable information for plant diseases control, including new strategies to prevent, delay, or inhibit fungal development. Our review provides a comprehensive overview of the efforts that have been made to decipher the effector proteins of biotrophic fungal pathogens and demonstrates how rapidly research in the field of obligate biotrophy has progressed.

Role of Sexual Reproduction in the Evolution of the Wheat Stripe Rust Fungus Races in China

Experimental and population genetic approaches have reshaped our view of how fungal pathogens reproduce, with consequences for our understanding of fungal invasions. Puccinia striiformis f. sp. tritici, the causal agent of stripe rust, poses a severe threat to wheat production worldwide. The sexual stage of P. striiformis f. sp. tritici was discovered >10 years ago, but how it affects the evolution of the pathogen, especially the emergence of the new virulent races, remains largely unknown. Here, using population genetic analyses, we demonstrate that sexual reproduction plays an important role in the evolution of P. striiformis f. sp. tritici races in China, specifically the newly emerged and devastating race virulent to resistance gene Yr26, which is widely used in China and exerts strong selective pressure on the pathogen population. Association analysis identified six genes encoding secreted proteins as candidates for virulence on wheat cultivars carrying the Yr26 resistance gene. Our results highlight the important role of sexual reproduction and selection exerted by hosts in the emergence of new virulent races in China.

Draft Genome Sequence Resource for the Orange Rust Pathogen of Sugarcane Puccinia kuehnii

Puccinia kuehnii is an obligate biotrophic fungal pathogen that causes orange rust of sugarcane, which is prevalent in many countries around the globe. In the United States, orange rust was first detected in sugarcane in Florida in 2007 and poses a persistent and economically damaging threat to the sugarcane industry in this region. Here, we generated the first genome assemblies for two isolates of P. kuehnii (1040 and 2143) collected in Florida in 2017 from two sugarcane cultivars, CL85-1040 and CP89-2143, respectively. These two rust genome resources will be of immense value for future genomic studies, particularly further exploration of the predicted secretomes that may help define key pathogenicity determinants for this economically important pathogen.

Multiple Antioxidative and Bioactive Molecules of Oats (Avena sativa L.) in Human Health

Oats (Avena sativa L.) are rich in protein, fiber, calcium, vitamins (B, C, E, and K), amino acids, and antioxidants (beta-carotene, polyphenols, chlorophyll, and flavonoids). β-glucan and avenanthramides improve the immune system, eliminate harmful substances from the body, reduce blood cholesterol, and help with dietary weight loss by enhancing the lipid profile and breaking down fat in the body. β-glucan regulates insulin secretion, preventing diabetes. Progladins also lower cholesterol levels, suppress the accumulation of triglycerides, reduce blood sugar levels, suppress inflammation, and improve skin health. Saponin-based avanacosidase and functional substances of flavone glycoside improve the immune function, control inflammation, and prevent infiltration in the skin. Moreover, lignin and phytoestrogen prevent hormone-related cancer and improve the quality of life of postmenopausal women. Sprouted oats are rich in saponarin in detoxifying the liver. The literatures have been reviewed and the recent concepts and prospects have been summarized with figures and tables. This review discusses recent trends in research on the functionality of oats rather than their nutritional value with individual immunity for self-medication. The oat and its acting components have been revisited for the future prospect and development of human healthy and functional sources.

A Chromosome-Scale Assembly of the Wheat Leaf Rust Pathogen Puccinia triticina Provides Insights Into Structural Variations and Genetic Relationships With Haplotype Resolution

Despite the global economic importance of the wheat leaf rust pathogen Puccinia triticina (Pt), genomic resources for Pt are limited and chromosome-level assemblies of Pt are lacking. Here, we present a complete haplotype-resolved genome assembly at a chromosome-scale for Pt using the Australian pathotype 64-(6),(7),(10),11 (Pt64; North American race LBBQB) built upon the newly developed technologies of PacBio and Hi-C sequencing. PacBio reads with ∼200-fold coverage (29.8 Gb data) were assembled by Falcon and Falcon-unzip and subsequently scaffolded with Hi-C data using Falcon-phase and Proximo. This approach allowed us to construct 18 chromosome pseudomolecules ranging from 3.5 to 12.3 Mb in size for each haplotype of the dikaryotic genome of Pt64. Each haplotype had a total length of ∼147 Mb, scaffold N₅₀ of ∼9.4 Mb, and was ∼93% complete for BUSCOs. Each haplotype had ∼29,800 predicted genes, of which ∼2,000 were predicted as secreted proteins (SPs). The investigation of structural variants (SVs) between haplotypes A and B revealed that 10% of the total genome was spanned by SVs, highlighting variations previously undetected by short-read based assemblies. For the first time, the mating type (MAT) genes on each haplotype of Pt64 were identified, which showed that MAT loci a and b are located on two chromosomes (chromosomes 7 and 14), representing a tetrapolar type. Furthermore, the Pt64 assembly enabled haplotype-based evolutionary analyses for 21 Australian Pt isolates, which highlighted the importance of a haplotype resolved reference when inferring genetic relationships using whole genome SNPs. This Pt64 assembly at chromosome-scale with full phase information provides an invaluable resource for genomic and evolutionary research, which will accelerate the understanding of molecular mechanisms underlying Pt-wheat interactions and facilitate the development of durable resistance to leaf rust in wheat and sustainable control of rust disease.

Genomic Characteristics and Comparative Genomics Analysis of Two Chinese Corynespora cassiicola Strains Causing Corynespora Leaf Fall (CLF) Disease

Rubber tree Corynespora leaf fall (CLF) disease, caused by the fungus Corynespora cassiicola, is one of the most damaging diseases in rubber tree plantations in Asia and Africa, and this disease also threatens rubber nurseries and young rubber plantations in China. C. cassiicola isolates display high genetic diversity, and virulence profiles vary significantly depending on cultivar. Although one phytotoxin (cassicolin) has been identified, it cannot fully explain the diversity in pathogenicity between C. cassiicola species, and some virulent C. cassiicola strains do not contain the cassiicolin gene. In the present study, we report high-quality gapless genome sequences, obtained using short-read sequencing and single-molecule long-read sequencing, of two Chinese C. cassiicola virulent strains. Comparative genomics of gene families in these two stains and a virulent CPP strain from the Philippines showed that all three strains experienced different selective pressures, and metabolism-related gene families vary between the strains. Secreted protein analysis indicated that the quantities of secreted cell wall-degrading enzymes were correlated with pathogenesis, and the most aggressive CCP strain (cassiicolin toxin type 1) encoded 27.34% and 39.74% more secreted carbohydrate-active enzymes (CAZymes) than Chinese strains YN49 and CC01, respectively, both of which can only infect rubber tree saplings. The results of antiSMASH analysis showed that all three strains encode ~60 secondary metabolite biosynthesis gene clusters (SM BGCs). Phylogenomic and domain structure analyses of core synthesis genes, together with synteny analysis of polyketide synthase (PKS) and non-ribosomal peptide synthetase (NRPS) gene clusters, revealed diversity in the distribution of SM BGCs between strains, as well as SM polymorphisms, which may play an important role in pathogenic progress. The results expand our understanding of the C. cassiicola genome. Further comparative genomic analysis indicates that secreted CAZymes and SMs may influence pathogenicity in rubber tree plantations. The findings facilitate future exploration of the molecular pathogenic mechanism of C. cassiicola.

Population Genomics of Fungal Plant Pathogens and the Analyses of Rapidly Evolving Genome Compartments

Genome sequencing of fungal pathogens have documented extensive variation in genome structure and composition between species and in many cases between individuals of the same species. This type of genomic variation can be adaptive for pathogens to rapidly evolve new virulence phenotypes. Analyses of genome-wide variation in fungal pathogen genomes rely on high quality assemblies and methods to detect and quantify structural variation. Population genomic studies in fungi have addressed the underlying mechanisms whereby structural variation can be rapidly generated. Transposable elements, high mutation and recombination rates as well as incorrect chromosome segregation during mitosis and meiosis contribute to extensive variation observed in many species. We here summarize key findings in the field of fungal pathogen genomics and we discuss methods to detect and characterize structural variants including an alignment-based pipeline to study variation in population genomic data.

Increased virulence of Puccinia coronata f. sp.avenae populations through allele frequency changes at multiple putative Avr loci

Pathogen populations are expected to evolve virulence traits in response to resistance deployed in agricultural settings. However, few temporal datasets have been available to characterize this process at the population level. Here, we examined two temporally separated populations of Puccinia coronata f. sp. avenae (Pca), which causes crown rust disease in oat (Avena sativa) sampled from 1990 to 2015. We show that a substantial increase in virulence occurred from 1990 to 2015 and this was associated with a genetic differentiation between populations detected by genome-wide sequencing. We found strong evidence for genetic recombination in these populations, showing the importance of the alternate host in generating genotypic variation through sexual reproduction. However, asexual expansion of some clonal lineages was also observed within years. Genome-wide association analysis identified seven Avr loci associated with virulence towards fifteen Pc resistance genes in oat and suggests that some groups of Pc genes recognize the same pathogen effectors. The temporal shift in virulence patterns in the Pca populations between 1990 and 2015 is associated with changes in allele frequency in these genomic regions. Nucleotide diversity patterns at a single Avr locus corresponding to Pc38, Pc39, Pc55, Pc63, Pc70, and Pc71 showed evidence of a selective sweep associated with the shift to virulence towards these resistance genes in all 2015 collected isolates.

A simple method for high molecular-weight genomic DNA extraction suitable for long-read sequencing from spores of an obligate biotroph oomycete

Long-read sequencing technologies are having a major impact on our approaches to studying non-model organisms and microbial communities. By significantly reducing the cost and facilitating the genome assembly pipelines, any laboratory can now develop its own genomics program regardless of the complexity of the genome studied. The most crucial current challenge is to develop efficient protocols for extracting genomic DNA (gDNA) with high quality and integrity adapted to the organism of interest. This can be particularly complex for obligate pathogens that must maintain intimate interactions inside infected host tissues. Here we propose a simple and cost-effective method for high molecular weight gDNA extraction from spores of Plasmopara halstedii, an obligate biotroph oomycete pathogen responsible for downy mildew in sunflower. We optimized the yield, the quality and the integrity of the extracted gDNA by fine-tuning three critical parameters, the grinding, the lysis temperature and the lysis duration. We obtained gDNA with a fragment size distribution reaching a peak ranging from 79 to 145 kb. More than half of the extracted gDNA consisted of DNA fragments larger than 42 kb, with 23% of fragments larger than 100 kb. We then demonstrated the relevance of this protocol for long-read sequencing using PacBio RSII technology. With this protocol, we were able to obtain a mean read length of 9.3 kb, a max read length of 71 kb and an N50 of 13.3 kb. The development of such DNA extraction protocols is an essential prerequisite for fully exploiting technologies requiring high molecular weight gDNA (e.g. long-read sequencing or optical mapping). These technological advances will help generate data to answer questions such as the role of newly duplicated gene clusters, repeated regions, genomic structural variations or to define number of chromosomes that still remains undefined in many species of pathogenic fungi and oomycetes.

Effector Biology of Biotrophic Plant Fungal Pathogens: Current Advances and Future Prospects

The interaction of fungal pathogens with their host requires a novel invading mechanism and the presence of various virulence-associated components responsible for promoting the infection. The small secretory proteins, explicitly known as effector proteins, are one of the prime mechanisms of host manipulation utilized by the pathogen to disarm the host. Several effector proteins are known to translocate from fungus to the plant cell for host manipulation. Many fungal effectors have been identified using genomic, transcriptomic, and bioinformatics approaches. Most of the effector proteins are devoid of any conserved signatures, and their prediction based on sequence homology is very challenging, therefore by combining the sequence consensus based upon machine learning features, multiple tools have also been developed for predicting apoplastic and cytoplasmic effectors. Various post-genomics approaches like transcriptomics of virulent isolates have also been utilized for identifying active consortia of effectors. Significant progress has been made in understanding biotrophic effectors; however, most of it is underway due to their complex interaction with host and complicated recognition and signaling networks. This review discusses advances, and challenges in effector identification and highlighted various features of the potential effector proteins and approaches for understanding their genetics and strategies for regulation.

Evolution of virulence in rust fungi - multiple solutions to one problem

Rust fungi are major pathogens that negatively affect crops and ecosystems. Recent rust disease epidemics driven by the emergence of strains with novel virulence profiles demand a better understanding of the evolutionary mechanisms of these organisms. Here, we review research advances in genome-scale analysis coupled with functional validation of effector candidate genes that have been instrumental to elucidate processes that contribute to changes in virulence phenotypes. We highlight how haplotype-phased genome references have paved the road to link these processes to the reproductive phases of rust fungi and have provided evidence for somatic exchange between strains as an important mechanism for generating diversity in asexual populations. With increasing data availability, we envision the future development of molecular virulence diagnostic tools.

Long-Read-Based de novo Genome Assembly and Comparative Genomics of the Wheat Leaf Rust Pathogen Puccinia triticina Identifies Candidates for Three Avirulence Genes

Leaf rust, caused by Puccinia triticina (Pt), is one of the most devastating diseases of wheat, affecting production in nearly all wheat-growing regions worldwide. Despite its economic importance, genomic resources for Pt are very limited. In the present study, we have used long-read sequencing (LRS) and the pipeline of FALCON and FALCON-Unzip (v4.1.0) to carry out the first LRS-based de novo genome assembly for Pt. Using 22.4-Gb data with an average read length of 11.6 kb and average coverage of 150-fold, we generated a genome assembly for Pt104 [strain 104-2,3,(6),(7),11; isolate S423], considered to be the founding isolate of a clonal lineage of Pt in Australia. The Pt104 genome contains 162 contigs with a total length of 140.5 Mb and N50 of 2 Mb, with the associated haplotigs providing haplotype information for 91% of the genome. This represents the best quality of Pt genome assembly to date, which reduces the contig number by 91-fold and improves the N50 by 4-fold as compared to the previous Pt race1 assembly. An annotation pipeline that combined multiple lines of evidence including the transcriptome assemblies derived from RNA-Seq, previously identified expressed sequence tags and Pt race 1 protein sequences predicted 29,043 genes for Pt104 genome. Based on the presence of a signal peptide, no transmembrane segment, and no target location to mitochondria, 2,178 genes were identified as secreted proteins (SPs). Whole-genome sequencing (Illumina paired-end) was performed for Pt104 and six additional strains with differential virulence profile on the wheat leaf rust resistance genes Lr26, Lr2a, and Lr3ka. To identify candidates for the corresponding avirulence genes AvrLr26, AvrLr2a, and AvrLr3ka, genetic variation within each strain was first identified by mapping to the Pt104 genome. Variants within predicted SP genes between the strains were then correlated to the virulence profiles, identifying 38, 31, and 37 candidates for AvrLr26, AvrLr2a, and AvrLr3ka, respectively. The identification of these candidate genes lays a good foundation for future studies on isolating these avirulence genes, investigating the molecular mechanisms underlying host-pathogen interactions, and the development of new diagnostic tools for pathogen monitoring.

How Target-Sequence Enrichment and Sequencing (TEnSeq) Pipelines Have Catalyzed Resistance Gene Cloning in the Wheat-Rust Pathosystem

The wheat-rust pathosystem has been well-studied among host-pathogen interactions since last century due to its economic importance. Intensified efforts toward cloning of wheat rust resistance genes commenced in the late 1990s with the first successful isolation published in 2003. Currently, a total of 24 genes have been cloned from wheat that provides resistance to stem rust, leaf rust, and stripe rust. Among them, more than half (15) were cloned over the last 4 years. This rapid cloning of resistance genes from wheat can be largely credited to the development of approaches for reducing the genome complexity as 10 out of the 15 genes cloned recently were achieved by approaches that are summarized as TEnSeq (Target-sequence Enrichment and Sequencing) pipelines in this review. The growing repertoire of cloned rust resistance genes provides new tools to support deployment strategies aimed at achieving durable resistance. This will be supported by the identification of genetic variation in corresponding Avr genes from rust pathogens, which has recently begun. Although developed with wheat resistance genes as the primary targets, TEnSeq approaches are also applicable to other classes of genes as well as for other crops with complex genomes.

Distinct Life Histories Impact Dikaryotic Genome Evolution in the Rust Fungus Puccinia striiformis Causing Stripe Rust in Wheat

Stripe rust of wheat, caused by the obligate biotrophic fungus Puccinia striiformis f.sp. tritici, is a major threat to wheat production worldwide with an estimated yearly loss of US $1 billion. The recent advances in long-read sequencing technologies and tailored-assembly algorithms enabled us to disentangle the two haploid genomes of Pst. This provides us with haplotype-specific information at a whole-genome level. Exploiting this novel information, we perform whole-genome comparative genomics of two P. striiformis f.sp. tritici isolates with contrasting life histories. We compare one isolate of the old European lineage (PstS0), which has been asexual for over 50 years, and a Warrior isolate (PstS7 lineage) from a novel incursion into Europe in 2011 from a sexual population in the Himalayan region. This comparison provides evidence that long-term asexual evolution leads to genome expansion, accumulation of transposable elements, and increased heterozygosity at the single nucleotide, structural, and allele levels. At the whole-genome level, candidate effectors are not compartmentalized and do not exhibit reduced levels of synteny. Yet we were able to identify two subsets of candidate effector populations. About 70% of candidate effectors are invariant between the two isolates, whereas 30% are hypervariable. The latter might be involved in host adaptation on wheat and explain the different phenotypes of the two isolates. Overall, this detailed comparative analysis of two haplotype-aware assemblies of P. striiformis f.sp. tritici is the first step in understanding the evolution of dikaryotic rust fungi at a whole-genome level.

Distinct Life Histories Impact Dikaryotic Genome Evolution in the Rust Fungus Puccinia striiformis Causing Stripe Rust in Wheat

Stripe rust of wheat, caused by the obligate biotrophic fungus Puccinia striiformis f.sp. tritici, is a major threat to wheat production worldwide with an estimated yearly loss of US $1 billion. The recent advances in long-read sequencing technologies and tailored-assembly algorithms enabled us to disentangle the two haploid genomes of Pst. This provides us with haplotype-specific information at a whole-genome level. Exploiting this novel information, we perform whole-genome comparative genomics of two P. striiformis f.sp. tritici isolates with contrasting life histories. We compare one isolate of the old European lineage (PstS0), which has been asexual for over 50 years, and a Warrior isolate (PstS7 lineage) from a novel incursion into Europe in 2011 from a sexual population in the Himalayan region. This comparison provides evidence that long-term asexual evolution leads to genome expansion, accumulation of transposable elements, and increased heterozygosity at the single nucleotide, structural, and allele levels. At the whole-genome level, candidate effectors are not compartmentalized and do not exhibit reduced levels of synteny. Yet we were able to identify two subsets of candidate effector populations. About 70% of candidate effectors are invariant between the two isolates, whereas 30% are hypervariable. The latter might be involved in host adaptation on wheat and explain the different phenotypes of the two isolates. Overall, this detailed comparative analysis of two haplotype-aware assemblies of P. striiformis f.sp. tritici is the first step in understanding the evolution of dikaryotic rust fungi at a whole-genome level.

Progress on Molecular Genetics and Manipulation of Rust Fungi

Among the thousands of rust species described, many are known for their devastating effects on their hosts, which include major agriculture crops and trees. Hence, for over a century, these basidiomycete pathogenic fungi have been researched and experimented with. However, due to their biotrophic nature, they are challenging organisms to work with and, needing their hosts for propagation, represent pathosystems that are not easily experimentally accessible. Indeed, efforts to perform genetics have been few and far apart for the rust fungi, though one study performed in the 1940s was famously instrumental in formulating the gene-for-gene hypothesis describing pathogen-host interactions. By taking full advantage of the molecular genetic tools developed in the 1980s, research on many plant pathogenic microbes thrived, yet similar work on the rusts remained very challenging though not without some successes. However, the genomics era brought real breakthrough research for the biotrophic fungi and with innovative experimentation and the use of heterologous systems, molecular genetic analyses over the last 2 decades have significantly advanced our insight into the function of many rust fungus genes and their role in the interaction with their hosts. This has allowed optimizing efforts for resistance breeding and the design and testing of various novel strategies to reduce the devastating diseases they cause.

Lack of gene flow between Phytophthora infestans populations of two neighboring countries with the largest potato production

Gene flow is an important evolutionary force that enables adaptive responses of plant pathogens in response to changes in the environment and plant disease management strategies. In this study, we made a direct inference concerning gene flow in the Irish famine pathogen Phytophthora infestans between two of its hosts (potato and tomato) as well as between China and India. This was done by comparing sequence characteristics of the eukaryotic translation elongation factor 1 alpha (eEF-1α) gene, generated from 245 P. infestans isolates sampled from two countries and hosts. Consistent with previous results, we found that eEF-1α gene was highly conserved and point mutation was the only mechanism generating any sequence variation. Higher genetic variation was found in the eEF-1α sequences in the P. infestans populations sampled from tomato compared to those sampled from potato. We also found the P. infestans population from India displayed a higher genetic variation in the eEF-1α sequences compared to China. No gene flow was detected between the pathogen populations from the two countries, which is possibly attributed to the geographic barrier caused by Himalaya Plateau and the minimum cross-border trade of potato and tomato products. The implications of these results for a sustainable management of late blight diseases are discussed.

Flax rust infection transcriptomics reveals a transcriptional profile that may be indicative for rust Avr genes

Secreted effectors of fungal pathogens are essential elements for disease development. However, lack of sequence conservation among identified effectors has long been a problem for predicting effector complements in fungi. Here we have explored the expression characteristics of avirulence (Avr) genes and candidate effectors of the flax rust fungus, Melampsora lini. We performed transcriptome sequencing and real-time quantitative PCR (qPCR) on RNA extracted from ungerminated spores, germinated spores, isolated haustoria and flax seedlings inoculated with M. lini isolate CH5 during plant infection. Genes encoding two categories of M. lini proteins, namely Avr proteins and plant cell wall degrading enzymes (CWDEs), were investigated in detail. Analysis of the expression profiles of 623 genes encoding predicted secreted proteins in the M. lini transcriptome shows that the six known Avr genes (i.e. AvrM (avrM), AvrM14, AvrL2, AvrL567, AvrP123 (AvrP) and AvrP4) fall within a group of 64 similarly expressed genes that are induced in planta and show a peak of expression early in infection with a subsequent decline towards sporulation. Other genes within this group include two paralogues of AvrL2, an AvrL567 virulence allele, and a number of genes encoding putative effector proteins. By contrast, M. lini genes encoding CWDEs fall into different expression clusters with their distribution often unrelated to their catalytic activity or substrate targets. These results suggest that synthesis of M. lini Avr proteins may be regulated in a coordinated fashion and that the expression profiling-based analysis has significant predictive power for the identification of candidate Avr genes.

Emergence of the Ug99 lineage of the wheat stem rust pathogen through somatic hybridisation

Parasexuality contributes to diversity and adaptive evolution of haploid (monokaryotic) fungi. However, non-sexual genetic exchange mechanisms are not defined in dikaryotic fungi (containing two distinct haploid nuclei). Newly emerged strains of the wheat stem rust pathogen, Puccinia graminis f. sp. tritici (Pgt), such as Ug99, are a major threat to global food security. Here, we provide genomics-based evidence supporting that Ug99 arose by somatic hybridisation and nuclear exchange between dikaryons. Fully haplotype-resolved genome assembly and DNA proximity analysis reveal that Ug99 shares one haploid nucleus genotype with a much older African lineage of Pgt, with no recombination or chromosome reassortment. These findings indicate that nuclear exchange between dikaryotes can generate genetic diversity and facilitate the emergence of new lineages in asexual fungal populations.

De Novo Genome Assembly and Comparative Genomics of the Barley Leaf Rust Pathogen Puccinia hordei Identifies Candidates for Three Avirulence Genes

Puccinia hordei (Ph) is a damaging pathogen of barley throughout the world. Despite its importance, almost nothing is known about the genomics of this pathogen, and a reference genome is lacking. In this study, the first reference genome was assembled for an Australian isolate of Ph ("Ph560") using long-read SMRT sequencing. A total of 838 contigs were assembled, with a total size of 207 Mbp. This included both haplotype collapsed and separated regions, consistent with an estimated haploid genome size of about 150Mbp. An annotation pipeline that combined RNA-Seq of Ph-infected host tissues and homology to proteins from four other Puccinia species predicted 25,543 gene models of which 1,450 genes were classified as encoding secreted proteins based on the prediction of a signal peptide and no transmembrane domain. Genome resequencing using short-read technology was conducted for four additional Australian strains, Ph612, Ph626, Ph608 and Ph584, which are considered to be simple mutational derivatives of Ph560 with added virulence to one or two of three barley leaf rust resistance genes (viz. Rph3, Rph13 and Rph19). To identify candidate genes for the corresponding avirulence genes AvrRph3, AvrRph13 and AvrRph19, genetic variation in predicted secreted protein genes between the strains was correlated to the virulence profiles of each, identifying 35, 29 and 46 candidates for AvrRph13, AvrRph3 and AvrRph19, respectively. The identification of these candidate genes provides a strong foundation for future efforts to isolate these three avirulence genes, investigate their biological properties, and develop new diagnostic tests for monitoring pathogen virulence.

Chromium sequencing: the doors open for genomics of obligate plant pathogens

It is challenging to sequence and assemble genomes of obligate plant pathogens and microorganisms because of limited amounts of DNA, comparatively large genomes and high numbers of repeat regions. We sequenced the 1.2 gigabase genome of an obligate rust fungus, Austropuccinia psidii, the cause of rust on Myrtaceae, with a Chromium 10X library. This technology has mostly been applied for single-cell sequencing in immunological studies of mammals. We compared scaffolds of a genome assembled from the Chromium library with one assembled from combined paired-end and mate-pair libraries, sequenced with Illumina HiSeq. Chromium 10X provided a superior assembly, in terms of number of scaffolds, N50 and number of genes recovered. It required less DNA than other methods and was sequenced and assembled at a lower cost. Chromium sequencing could provide a solution to sequence and assemble genomes of obligate plant pathogens where the amount of available DNA is a limiting factor.

Virulence Structure and Diversity of Puccinia coronata f. sp. avenae P. Syd. & Syd. in Poland During 2013 to 2015

The crown rust fungus Puccinia coronata f. sp. avenae P. Syd. & Syd. (Pca) attacks cultivated oat and its wild relatives, causing significant losses to the crop worldwide. Although understanding the origin and dynamics of the pathogen's diversity is critical to developing methods for its control, there are little relevant data on Pca virulence diversity in Europe, the global center of oat production. The goal of this study was to analyze the diversity of Pca populations in Poland in 2013 to 2015 based on their ability to overcome currently available host resistance Pc genes. Pca isolate virulence was evaluated on a panel of lines containing 26 major resistance genes of oat. The isolates were able to overcome from 1 to 16 resistance genes each, with most isolates being virulent on five to seven lines. In all years, a very high level of crown rust pathotype diversity was observed, with Simpson and Evenness indices of 0.99. In total, 156 different pathotypes were detected, with no prevalent pathotype in any of the 3 years analyzed. The results showed that the virulence level of P. coronata isolates was relatively low for each year studied (21% on average), most likely owing to the low take up of Pc genes in Polish oat cultivars, meaning that many sources of resistance are still effective against Pca races occurring in Poland. The long-range dispersal of Puccinia spores supported by the availability of wild, weedy, and cultivated Avena species makes it likely that the virulence profile seen in Poland is representative of much of central Europe and beyond.

A chromosome-scale genome assembly reveals a highly dynamic effector repertoire of wheat powdery mildew

Blumeria graminis f. sp. tritici (B.g. tritici) is the causal agent of the wheat powdery mildew disease. The highly fragmented B.g. tritici genome available so far has prevented a systematic analysis of effector genes that are known to be involved in host adaptation. To study the diversity and evolution of effector genes we produced a chromosome-scale assembly of the B.g. tritici genome. The genome assembly and annotation was achieved by combining long-read sequencing with high-density genetic mapping, bacterial artificial chromosome fingerprinting and transcriptomics. We found that the 166.6 Mb B.g. tritici genome encodes 844 candidate effector genes, over 40% more than previously reported. Candidate effector genes have characteristic local genomic organization such as gene clustering and enrichment for recombination-active regions and certain transposable element families. A large group of 412 candidate effector genes shows high plasticity in terms of copy number variation in a global set of 36 isolates and of transcription levels. Our data suggest that copy number variation and transcriptional flexibility are the main drivers for adaptation in B.g. tritici. The high repeat content may play a role in providing a genomic environment that allows rapid evolution of effector genes with selection as the driving force.

Structural and functional insights into the modulation of the activity of a flax cytokinin oxidase by flax rust effector AvrL567-A

During infection, plant pathogens secrete effector proteins to facilitate colonization. In comparison with our knowledge of bacterial effectors, the current understanding of how fungal effectors function is limited. In this study, we show that the effector AvrL567-A from the flax rust fungus Melampsora lini interacts with a flax cytosolic cytokinin oxidase, LuCKX1.1, using both yeast two-hybrid and in planta bimolecular fluorescence assays. Purified LuCKX1.1 protein shows catalytic activity against both N6-(Δ2-isopentenyl)-adenine (2iP) and trans-zeatin (tZ) substrates. Incubation of LuCKX1.1 with AvrL567-A results in increased catalytic activity against both substrates. The crystal structure of LuCKX1.1 and docking studies with AvrL567-A indicate that the AvrL567 binding site involves a flexible surface-exposed region that surrounds the cytokinin substrate access site, which may explain its effect in modulating LuCKX1.1 activity. Expression of AvrL567-A in transgenic flax plants gave rise to an epinastic leaf phenotype consistent with hormonal effects, although no difference in overall cytokinin levels was observed. We propose that, during infection, plant pathogens may differentially modify the levels of extracellular and intracellular cytokinins.

Detection of Race-Specific Resistance Against Puccinia coronata f. sp. avenae in Brachypodium Species

Oat crown rust caused by Puccinia coronata f. sp. avenae is the most destructive foliar disease of cultivated oat. Characterization of genetic factors controlling resistance responses to Puccinia coronata f. sp. avenae in nonhost species could provide new resources for developing disease protection strategies in oat. We examined symptom development and fungal colonization levels of a collection of Brachypodium distachyon and B. hybridum accessions infected with three North American P. coronata f. sp. avenae isolates. Our results demonstrated that colonization phenotypes are dependent on both host and pathogen genotypes, indicating a role for race-specific responses in these interactions. These responses were independent of the accumulation of reactive oxygen species. Expression analysis of several defense-related genes suggested that salicylic acid and ethylene-mediated signaling but not jasmonic acid are components of resistance reaction to P. coronata f. sp. avenae. Our findings provide the basis to conduct a genetic inheritance study to examine whether effector-triggered immunity contributes to nonhost resistance to P. coronata f. sp. avenae in Brachypodium spp.