Elevated CO2 Can Worsen Fusarium Head Blight Disease Severity in Wheat but the Fhb1 QTL Provides Reliable Disease Resistance.

Fusarium head blight (FHB) is a destructive fungal disease of wheat that causes significant economic loss due to lower yields and the contamination of grain with fungal toxins (mycotoxins), particularly deoxynivalenol (DON). FHB disease spread and mycotoxin contamination has been shown to worsen at elevated CO2, therefore, it is important to identify climate-resilient FHB resistance. This work evaluates whether wheat with the Fhb1 quantitative trait locus (QTL), the most widely deployed FHB resistance locus in wheat breeding programs, provides reliable disease resistance at elevated CO2. Near-isogenic wheat lines (NILs) derived from either a highly FHB susceptible or a more FHB resistant genetic background, with or without the Fhb1 QTL, were grown in growth chambers at ambient (400 ppm) and elevated (1000 ppm) CO2 conditions. Wheat was inoculated with Fusarium graminearum and evaluated for FHB severity. At elevated CO2, the NILs derived from more FHB-resistant wheat had increased disease spread, greater pathogen biomass and mycotoxin contamination, and lower rates of DON detoxification; this was not observed in wheat from a FHB susceptible genetic background. The Fhb1 QTL was not associated with increased disease severity in wheat grown at elevated CO2 and provided reliable disease resistance.

Fhb1 disease resistance QTL does not exacerbate wheat grain protein loss at elevated CO2.

Fusarium head blight, a devastating cereal crop disease, can cause significant yield losses and contaminate grain with hazardous fungal toxins. Concerningly, recent evidence indicates that substantial grain protein content loss is likely to occur in wheat that is moderately resistant to head blight when it is grown at elevated CO2. Although wheat breeders in North America utilize a number of resistance sources and genes to reduce pathogen damage, the Fhb1 gene is widely deployed. To determine whether Fhb1 is associated with the protein content loss at elevated CO2, twelve near-isogenic spring wheat lines from either a susceptible or moderately susceptible genetic background, and with, or without the Fhb1 QTL, were grown at ambient and elevated CO2 conditions. The near-isogenic lines were evaluated for differences in physiology, productivity, and grain protein content. Our results showed that the Fhb1 QTL did not have any significant effect on plant growth, development, yield, or grain protein content at ambient or elevated CO2. Therefore, other factors in the moderately susceptible wheat genetic background are likely responsible for the more severe grain protein loss at elevated CO2.

A Homolog of the Arabidopsis TIME FOR COFFEE Gene Is Involved in Nonhost Resistance to Wheat Stem Rust in Brachypodium distachyon.

Plants resist infection by pathogens using both preexisting barriers and inducible defense responses. Inducible responses are governed in a complex manner by various hormone signaling pathways. The relative contribution of hormone signaling pathways to nonhost resistance to pathogens is not well understood. In this study, we examined the molecular basis of disrupted nonhost resistance to the fungal species Puccinia graminis, which causes stem rust of wheat, in an induced mutant of the model grass Brachypodium distachyon. Through bioinformatic analysis, a 1-bp deletion in the mutant genotype was identified that introduces a premature stop codon in the gene Bradi1g24100, which is a homolog of the Arabidopsis thaliana gene TIME FOR COFFEE (TIC). In Arabidopsis, TIC is central to the regulation of the circadian clock and plays a crucial role in jasmonate signaling by attenuating levels of the transcription factor protein MYC2, and its mutational disruption results in enhanced susceptibility to the hemibiotroph Pseudomonas syringae. Our similar finding for an obligate biotroph suggests that the biochemical role of TIC in mediating disease resistance to biotrophs is conserved in grasses, and that the correct modulation of jasmonate signaling during infection by Puccinia graminis may be essential for nonhost resistance to wheat stem rust in B. distachyon.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Heritable temporal gene expression patterns correlate with metabolomic seed content in developing hexaploid oat seed.

Oat ranks sixth in world cereal production and has a higher content of health-promoting compounds compared with other cereals. However, there is neither a robust oat reference genome nor transcriptome. Using deeply sequenced full-length mRNA libraries of oat cultivar Ogle-C, a de novo high-quality and comprehensive oat seed transcriptome was assembled. With this reference transcriptome and QuantSeq 3' mRNA sequencing, gene expression was quantified during seed development from 22 diverse lines across six time points. Transcript expression showed higher correlations between adjacent time points. Based on differentially expressed genes, we identified 22 major temporal co-expression (TCoE) patterns of gene expression and revealed enriched gene ontology biological processes. Within each TCoE set, highly correlated transcripts, putatively commonly affected by genetic background, were clustered and termed genetic co-expression (GCoE) sets. Seventeen of the 22 TCoE sets had GCoE sets with median heritabilities higher than 0.50, and these heritability estimates were much higher than that estimated from permutation analysis, with no divergence observed in cluster sizes between permutation and non-permutation analyses. Linear regression between 634 metabolites from mature seeds and the PC1 score of each of the GCoE sets showed significantly lower p-values than permutation analysis. Temporal expression patterns of oat avenanthramides and lipid biosynthetic genes were concordant with previous studies of avenanthramide biosynthetic enzyme activity and lipid accumulation. This study expands our understanding of physiological processes that occur during oat seed maturation and provides plant breeders the means to change oat seed composition through targeted manipulation of key pathways.

Genetic Variation for Seed Metabolite Levels in Brachypodium distachyon.

Metabolite composition and concentrations in seed grains are important traits of cereals. To identify the variation in the seed metabolotypes of a model grass, namely Brachypodium distachyon, we applied a widely targeted metabolome analysis to forty inbred lines of B. distachyon and examined the accumulation patterns of 183 compounds in the seeds. By comparing the metabolotypes with the population structure of these lines, we found signature metabolites that represent different accumulation patterns for each of the three B. distachyon subpopulations. Moreover, we found that thirty-seven metabolites exhibited significant differences in their accumulation between the lines Bd21 and Bd3-1. Using a recombinant inbred line (RIL) population from a cross between Bd3-1 and Bd21, we identified the quantitative trait loci (QTLs) linked with this variation in the accumulation of thirteen metabolites. Our metabolite QTL analysis illustrated that different genetic factors may presumably regulate the accumulation of 4-pyridoxate and pyridoxamine in vitamin B6 metabolism. Moreover, we found two QTLs on chromosomes 1 and 4 that affect the accumulation of an anthocyanin, chrysanthemin. These QTLs genetically interacted to regulate the accumulation of this compound. This study demonstrates the potential for metabolite QTL mapping in B. distachyon and provides new insights into the genetic dissection of metabolomic traits in temperate grasses.

Genomic Dissection of Nonhost Resistance to Wheat Stem Rust in Brachypodium distachyon.

The emergence of new races of Puccinia graminis f. sp. tritici, the causal pathogen of wheat stem rust, has spurred interest in developing durable resistance to this disease in wheat. Nonhost resistance holds promise to help control this and other diseases because it is durable against nonadapted pathogens. However, the genetic and molecular basis of nonhost resistance to wheat stem rust is poorly understood. In this study, the model grass Brachypodium distachyon, a nonhost of P. graminis f. sp. tritici, was used to genetically dissect nonhost resistance to wheat stem rust. A recombinant inbred line (RIL) population segregating for response to wheat stem rust was evaluated for resistance. Evaluation of genome-wide cumulative single nucleotide polymorphism allele frequency differences between contrasting pools of resistant and susceptible RILs followed by molecular marker analysis identified six quantitative trait loci (QTL) that cumulatively explained 72.5% of the variation in stem rust resistance. Two of the QTLs explained 31.7% of the variation, and their interaction explained another 4.6%. Thus, nonhost resistance to wheat stem rust in B. distachyon is genetically complex, with both major and minor QTLs acting additively and, in some cases, interacting. These findings will guide future research to identify genes essential to nonhost resistance to wheat stem rust.

Quantitative Trait Loci Associated with Drought Tolerance in Brachypodium distachyon.

The temperate wild grass Brachypodium distachyon (Brachypodium) serves as model system for studying turf and forage grasses. Brachypodium collections show diverse responses to drought stress, but little is known about the genetic mechanisms of drought tolerance of this species. The objective of this study was to identify quantitative trait loci (QTLs) associated with drought tolerance traits in Brachypodium. We assessed leaf fresh weight (LFW), leaf dry weight (LDW), leaf water content (LWC), leaf wilting (WT), and chlorophyll fluorescence (Fv/Fm) under well-watered and drought conditions on a recombinant inbred line (RIL) population from two parents (Bd3-1 and Bd1-1) known to differ in their drought adaptation. A linkage map of the RIL population was constructed using 467 single nucleotide polymorphism (SNP) markers obtained from genotyping-by-sequencing. The Bd3-1/Bd1-1 map spanned 1,618 cM and had an average distance of 3.5 cM between adjacent single nucleotide polymorphisms (SNPs). Twenty-six QTLs were identified in chromosome 1, 2, and 3 in two experiments, with 14 of the QTLs under well-watered conditions and 12 QTLs under drought stress. In Experiment 1, a QTL located on chromosome 2 with a peak at 182 cM appeared to simultaneously control WT, LWC, and Fv/Fm under drought stress, accounting for 11-18.7% of the phenotypic variation. Allelic diversity of candidate genes DREB2B, MYB, and SPK, which reside in one multi-QTL region, may play a role in the natural variation in whole plant drought tolerance in Brachypodium. Co-localization of QTLs for multiple drought-related traits suggest that the gene(s) involved are important regulators of drought tolerance in Brachypodium.

De Novo Transcriptome Assembly in Polyploid Species.

In the absence of a reference genome, the ultimate goal of a de novo transcriptome assembly is to accurately and comprehensively reconstruct the set of messenger RNA transcripts represented in the sample. Non-reference assembly of the transcriptome of polyploid species poses a particular challenge because of the presence of homeologs that are difficult to disentangle at the sequence level. This is especially true for hexaploid oats, which have three highly similar subgenomes, two of which are thought to be nearly identical. Under these circumstances, most software packages and established pipelines encounter difficulties in rendering an accurate transcriptome because they are typically developed, refined, and tested for diploid organisms. We present a protocol for transcriptome assembly in oats that can be extended both to other polyploids and species with highly duplicated genomes.

Genetic Architecture of Flowering-Time Variation in Brachypodium distachyon.

The transition to reproductive development is a crucial step in the plant life cycle, and the timing of this transition is an important factor in crop yields. Here, we report new insights into the genetic control of natural variation in flowering time in Brachypodium distachyon, a nondomesticated pooid grass closely related to cereals such as wheat (Triticum spp.) and barley (Hordeum vulgare L.). A recombinant inbred line population derived from a cross between the rapid-flowering accession Bd21 and the delayed-flowering accession Bd1-1 were grown in a variety of environmental conditions to enable exploration of the genetic architecture of flowering time. A genotyping-by-sequencing approach was used to develop SNP markers for genetic map construction, and quantitative trait loci (QTLs) that control differences in flowering time were identified. Many of the flowering-time QTLs are detected across a range of photoperiod and vernalization conditions, suggesting that the genetic control of flowering within this population is robust. The two major QTLs identified in undomesticated B. distachyon colocalize with VERNALIZATION1/PHYTOCHROME C and VERNALIZATION2, loci identified as flowering regulators in the domesticated crops wheat and barley. This suggests that variation in flowering time is controlled in part by a set of genes broadly conserved within pooid grasses.

Quantitative trait loci associated with natural diversity in water-use efficiency and response to soil drying in Brachypodium distachyon.

All plants must optimize their growth with finite resources. Water use efficiency (WUE) measures the relationship between biomass acquisition and transpired water. In the present study, we performed two experiments to understand the genetic basis of WUE and other parameters of plant-water interaction under control and water-limited conditions. Our study used two inbred natural accessions of Brachypodium distachyon, a model grass species with close phylogenetic affinity to temperate forage and cereal crops. First, we identify the soil water content which causes a reduction in leaf relative water content and an increase in WUE. Second, we present results from a large phenotyping experiment utilizing a recombinant inbred line mapping population derived from these same two natural accessions. We identify QTLs associated with environmentally-insensitive genetic variation in WUE, including a pair of epistatically interacting loci. We also identify QTLs associated with constitutive differences in biomass and a QTL describing an environmentally-sensitive difference in leaf carbon content. Finally, we present a new linkage map for this mapping population based on new SNP markers as well as updated genomic positions for previously described markers. Our studies provide an initial characterization of plant-water relations in B. distachyon and identify candidate genomic regions involved in WUE.

Cell Wall Composition and Biomass Recalcitrance Differences Within a Genotypically Diverse Set of Brachypodium distachyon Inbred Lines.

Brachypodium distachyon (Brachypodium) has emerged as a useful model system for studying traits unique to graminaceous species including bioenergy crop grasses owing to its amenability to laboratory experimentation and the availability of extensive genetic and germplasm resources. Considerable natural variation has been uncovered for a variety of traits including flowering time, vernalization responsiveness, and above-ground growth characteristics. However, cell wall composition differences remain underexplored. Therefore, we assessed cell wall-related traits relevant to biomass conversion to biofuels in seven Brachypodium inbred lines that were chosen based on their high level of genotypic diversity as well as available genome sequences and recombinant inbred line (RIL) populations. Senesced stems plus leaf sheaths from these lines exhibited significant differences in acetyl bromide soluble lignin (ABSL), cell wall polysaccharide-derived sugars, hydroxycinnamates content, and syringyl:guaiacyl:p-hydroxyphenyl (S:G:H) lignin ratios. Free glucose, sucrose, and starch content also differed significantly in senesced stems, as did the amounts of sugars released from cell wall polysaccharides (digestibility) upon exposure to a panel of thermochemical pretreatments followed by hydrolytic enzymatic digestion. Correlations were identified between inbred line lignin compositions and plant growth characteristics such as biomass accumulation and heading date (HD), and between amounts of cell wall polysaccharides and biomass digestibility. Finally, stem cell wall p-coumarate and ferulate contents and free-sugars content changed significantly with increased duration of vernalization for some inbred lines. Taken together, these results show that Brachypodium displays substantial phenotypic variation with respect to cell wall composition and biomass digestibility, with some compositional differences correlating with growth characteristics. Moreover, besides influencing HD and biomass accumulation, vernalization was found to affect cell wall composition and free sugars accumulation in some Brachypodium inbred lines, suggesting genetic differences in how vernalization affects carbon flux to polysaccharides. The availability of related RIL populations will allow for the genetic and molecular dissection of this natural variation, the knowledge of which may inform ways to genetically improve bioenergy crop grasses.

Subgenome-specific assembly of vitamin E biosynthesis genes and expression patterns during seed development provide insight into the evolution of oat genome.

Vitamin E is essential for humans and thus must be a component of a healthy diet. Among the cereal grains, hexaploid oats (Avena sativa L.) have high vitamin E content. To date, no gene sequences in the vitamin E biosynthesis pathway have been reported for oats. Using deep sequencing and orthology-guided assembly, coding sequences of genes for each step in vitamin E synthesis in oats were reconstructed, including resolution of the sequences of homeologs. Three homeologs, presumably representing each of the three oat subgenomes, were identified for the main steps of the pathway. Partial sequences, likely representing pseudogenes, were recovered in some instances as well. Pairwise comparisons among homeologs revealed that two of the three putative subgenome-specific homeologs are almost identical for each gene. Synonymous substitution rates indicate the time of divergence of the two more similar subgenomes from the distinct one at 7.9-8.7 MYA, and a divergence between the similar subgenomes from a common ancestor 1.1 MYA. A new proposed evolutionary model for hexaploid oat formation is discussed. Homeolog-specific gene expression was quantified during oat seed development and compared with vitamin E accumulation. Homeolog expression largely appears to be similar for most of genes; however, for some genes, homoeolog-specific transcriptional bias was observed. The expression of HPPD, as well as certain homoeologs of VTE2 and VTE4, is highly correlated with seed vitamin E accumulation. Our findings expand our understanding of oat genome evolution and will assist efforts to modify vitamin E content and composition in oats.

A spontaneous segmental deletion from chromosome arm 3DL enhances Fusarium head blight resistance in wheat.

Much effort has been directed at identifying sources of resistance to Fusarium head blight (FHB) in wheat. We sought to identify molecular markers for what we hypothesized was a new major FHB resistance locus originating from the wheat cultivar 'Freedom' and introgressed into the susceptible wheat cultivar 'USU-Apogee'. An F2:3 mapping population from a cross between Apogee and A30, its BC4 near-isoline exhibiting improved FHB resistance, was evaluated for resistance. The distribution of FHB resistance in the population approximated a 1:3 moderately resistant : moderately susceptible + susceptible ratio. Separate disease evaluations established that A30 accumulated less deoxynivalenol and yielded a greater proportion of sound grain than Apogee. Molecular mapping revealed that the FHB resistance of A30 is associated with molecular markers on chromosome arm 3DL that exhibit a null phenotype in A30 but are present in both Apogee and Freedom, indicating a spontaneous deletion occurred during the development of A30. Aneuploid analysis revealed that the size of the deleted segment is approximately 19% of the arm's length. Our results suggest that the deleted interval of chromosome arm 3DL in Apogee may harbor FHB susceptibility genes that promote disease spread in infected spikes, and that their elimination increases FHB resistance in a novel manner.

Brachypodium distachyon-Cochliobolus sativus Pathosystem is a New Model for Studying Plant-Fungal Interactions in Cereal Crops.

Cochliobolus sativus (anamorph: Bipolaris sorokiniana) causes spot blotch, common root rot, and kernel blight or black point in barley and wheat. However, little is known about the molecular mechanisms underlying the pathogenicity of C. sativus or the molecular basis of resistance and susceptibility in the hosts. This study aims to establish the model grass Brachypodium distachyon as a new model for studying plant-fungus interactions in cereal crops. Six B. distachyon lines were inoculated with five C. sativus isolates. The results indicated that all six B. distachyon lines were infected by the C. sativus isolates, with their levels of resistance varying depending on the fungal isolates used. Responses ranging from hypersensitive response-mediated resistance to complete susceptibility were observed in a large collection of B. distachyon (2n=2x=10) and B. hybridum (2n=4x=30) accessions inoculated with four of the C. sativus isolates. Evaluation of an F2 population derived from the cross between two of the B. distachyon lines, Bd1-1 and Bd3-1, with isolate Cs07-47-1 showed quantitative and transgressive segregation for resistance to C. sativus, suggesting that the resistance may be governed by quantitative trait loci from both parents. The availability of whole-genome sequences of both the host (B. distachyon) and the pathogen (C. sativus) makes this pathosystem an attractive model for studying this important disease of cereal crops.

Genome diversity in Brachypodium distachyon: deep sequencing of highly diverse inbred lines.

Brachypodium distachyon is small annual grass that has been adopted as a model for the grasses. Its small genome, high-quality reference genome, large germplasm collection, and selfing nature make it an excellent subject for studies of natural variation. We sequenced six divergent lines to identify a comprehensive set of polymorphisms and analyze their distribution and concordance with gene expression. Multiple methods and controls were utilized to identify polymorphisms and validate their quality. mRNA-Seq experiments under control and simulated drought-stress conditions, identified 300 genes with a genotype-dependent treatment response. We showed that large-scale sequence variants had extremely high concordance with altered expression of hundreds of genes, including many with genotype-dependent treatment responses. We generated a deep mRNA-Seq dataset for the most divergent line and created a de novo transcriptome assembly. This led to the discovery of >2400 previously unannotated transcripts and hundreds of genes not present in the reference genome. We built a public database for visualization and investigation of sequence variants among these widely used inbred lines.

Update on the genomics and basic biology of Brachypodium: International Brachypodium Initiative (IBI).

The scientific presentations at the First International Brachypodium Conference (abstracts available at http://www.brachy2013.unimore.it) are evidence of the widespread adoption of Brachypodium distachyon as a model system. Furthermore, the wide range of topics presented (genome evolution, roots, abiotic and biotic stress, comparative genomics, natural diversity, and cell walls) demonstrates that the Brachypodium research community has achieved a critical mass of tools and has transitioned from resource development to addressing biological questions, particularly those unique to grasses.

Molecular and comparative mapping of genes governing spike compactness from wild emmer wheat.

The development and morphology of the wheat spike is important because the spike is where reproduction occurs and it holds the grains until harvest. Therefore, genes that influence spike morphology are of interest from both theoretical and practical stand points. When substituted for the native chromosome 2A in the tetraploid Langdon (LDN) durum wheat background, the Triticum turgidum ssp. dicoccoides chromosome 2A from accession IsraelA confers a short, compact spike with fewer spikelets per spike compared to LDN. Molecular mapping and quantitative trait loci (QTL) analysis of these traits in a homozygous recombinant population derived from LDN × the chromosome 2A substitution line (LDNIsA-2A) indicated that the number of spikelets per spike and spike length were controlled by linked, but different, loci on the long arm of 2A. A QTL explaining most of the variation for spike compactness coincided with the QTL for spike length. Comparative mapping indicated that the QTL for number of spikelets per spike overlapped with a previously mapped QTL for Fusarium head blight susceptibility. The genes governing spike length and compactness were not orthologous to either sog or C, genes known to confer compact spikes in diploid and hexaploid wheat, respectively. Mapping and sequence analysis indicated that the gene governing spike length and compactness derived from wild emmer could be an ortholog of the barley Cly1/Zeo gene, which research indicates is an AP2-like gene pleiotropically affecting cleistogamy, flowering time, and rachis internode length. This work provides researchers with knowledge of new genetic loci and associated markers that may be useful for manipulating spike morphology in durum wheat.

Parallel analysis of RNA ends enhances global investigation of microRNAs and target RNAs of Brachypodium distachyon.

BACKGROUND: The wild grass Brachypodium distachyon has emerged as a model system for temperate grasses and biofuel plants. However, the global analysis of miRNAs, molecules known to be key for eukaryotic gene regulation, has been limited in B. distachyon to studies examining a few samples or that rely on computational predictions. Similarly an in-depth global analysis of miRNA-mediated target cleavage using parallel analysis of RNA ends (PARE) data is lacking in B. distachyon. RESULTS: B. distachyon small RNAs were cloned and deeply sequenced from 17 libraries that represent different tissues and stresses. Using a computational pipeline, we identified 116 miRNAs including not only conserved miRNAs that have not been reported in B. distachyon, but also non-conserved miRNAs that were not found in other plants. To investigate miRNA-mediated cleavage function, four PARE libraries were constructed from key tissues and sequenced to a total depth of approximately 70 million sequences. The roughly 5 million distinct genome-matched sequences that resulted represent an extensive dataset for analyzing small RNA-guided cleavage events. Analysis of the PARE and miRNA data provided experimental evidence for miRNA-mediated cleavage of 264 sites in predicted miRNA targets. In addition, PARE analysis revealed that differentially expressed miRNAs in the same family guide specific target RNA cleavage in a correspondingly tissue-preferential manner. CONCLUSIONS: B. distachyon miRNAs and target RNAs were experimentally identified and analyzed. Knowledge gained from this study should provide insights into the roles of miRNAs and the regulation of their targets in B. distachyon and related plants.

Analysis and annotation of the hexaploid oat seed transcriptome.

BACKGROUND: Next generation sequencing provides new opportunities to explore transcriptomes. However, challenges remain for accurate differentiation of homoeoalleles and paralogs, particularly in polyploid organisms with no supporting genome sequence. In this study, RNA-Seq was employed to generate and characterize the first gene expression atlas for hexaploid oat. RESULTS: The software packages Trinity and Oases were used to produce a transcript assembly from nearly 134 million 100-bp paired-end reads from developing oat seeds. Based on the quality-parameters employed, Oases assemblies were superior. The Oases 67-kmer assembly, denoted dnOST (de novo Oat Seed Transcriptome), is over 55 million nucleotides in length and the average transcript length is 1,043 nucleotides. The 74.8× sequencing depth was adequate to differentiate a large proportion of putative homoeoalleles and paralogs. To assess the robustness of dnOST, we successfully identified gene transcripts associated with the biosynthetic pathways of three compounds with health-promoting properties (avenanthramides, tocols, ß-glucans), and quantified their expression. CONCLUSIONS: To our knowledge, this study provides the first direct performance comparison between two major assemblers in a polyploid organism. The workflow we developed provides a useful guide for comparable analyses in other organisms. The transcript assembly developed here is a major advance. It expands the number of oat ESTs 3-fold, and constitutes the first comprehensive transcriptome study in oat. This resource will be a useful new tool both for analysis of genes relevant to nutritional enhancement of oat, and for improvement of this crop in general.

Infection of Brachypodium distachyon by formae speciales of Puccinia graminis: early infection events and host-pathogen incompatibility.

Puccinia graminis causes stem rust, a serious disease of cereals and forage grasses. Important formae speciales of P. graminis and their typical hosts are P. graminis f. sp. tritici (Pg-tr) in wheat and barley, P. graminis f. sp. lolii (Pg-lo) in perennial ryegrass and tall fescue, and P. graminis f. sp. phlei-pratensis (Pg-pp) in timothy grass. Brachypodium distachyon is an emerging genetic model to study fungal disease resistance in cereals and temperate grasses. We characterized the P. graminis-Brachypodium pathosystem to evaluate its potential for investigating incompatibility and non-host resistance to P. graminis. Inoculation of eight Brachypodium inbred lines with Pg-tr, Pg-lo or Pg-pp resulted in sporulating lesions later accompanied by necrosis. Histological analysis of early infection events in one Brachypodium inbred line (Bd1-1) indicated that Pg-lo and Pg-pp were markedly more efficient than Pg-tr at establishing a biotrophic interaction. Formation of appressoria was completed (60-70% of germinated spores) by 12 h post-inoculation (hpi) under dark and wet conditions, and after 4 h of subsequent light exposure fungal penetration structures (penetration peg, substomatal vesicle and primary infection hyphae) had developed. Brachypodium Bd1-1 exhibited pre-haustorial resistance to Pg-tr, i.e. infection usually stopped at appressorial formation. By 68 hpi, only 0.3% and 0.7% of the Pg-tr urediniospores developed haustoria and colonies, respectively. In contrast, development of advanced infection structures by Pg-lo and Pg-pp was significantly more common; however, Brachypodium displayed post-haustorial resistance to these isolates. By 68 hpi the percentage of urediniospores that only develop a haustorium mother cell or haustorium in Pg-lo and Pg-pp reached 8% and 5%, respectively. The formation of colonies reached 14% and 13%, respectively. We conclude that Brachypodium is an apt grass model to study the molecular and genetic components of incompatiblity and non-host resistance to P. graminis.