Melania Figueroa

Comparative genomics of a plant-pathogenic fungus, Pyrenophora tritici-repentis, reveals transduplication and the impact of repeat elements on pathogenicity and population divergence.

Pyrenophora tritici-repentis is a necrotrophic fungus causal to the disease tan spot of wheat, whose contribution to crop loss has increased significantly during the last few decades. Pathogenicity by this fungus is attributed to the production of host-selective toxins (HST), which are recognized by their host in a genotype-specific manner. To better understand the mechanisms that have led to the increase in disease incidence related to this pathogen, we sequenced the genomes of three P. tritici-repentis isolates. A pathogenic isolate that produces two known HSTs was used to assemble a reference nuclear genome of approximately 40 Mb composed of 11 chromosomes that encode 12,141 predicted genes. Comparison of the reference genome with those of a pathogenic isolate that produces a third HST, and a nonpathogenic isolate, showed the nonpathogen genome to be more diverged than those of the two pathogens. Examination of gene-coding regions has provided candidate pathogen-specific proteins and revealed gene families that may play a role in a necrotrophic lifestyle. Analysis of transposable elements suggests that their presence in the genome of pathogenic isolates contributes to the creation of novel genes, effector diversification, possible horizontal gene transfer events, identified copy number variation, and the first example of transduplication by DNA transposable elements in fungi. Overall, comparative analysis of these genomes provides evidence that pathogenicity in this species arose through an influx of transposable elements, which created a genetically flexible landscape that can easily respond to environmental changes.

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.

Widespread Polycistronic Transcripts in Fungi Revealed by Single-Molecule mRNA Sequencing.

Genes in prokaryotic genomes are often arranged into clusters and co-transcribed into polycistronic RNAs. Isolated examples of polycistronic RNAs were also reported in some higher eukaryotes but their presence was generally considered rare. Here we developed a long-read sequencing strategy to identify polycistronic transcripts in several mushroom forming fungal species including Plicaturopsis crispa, Phanerochaete chrysosporium, Trametes versicolor, and Gloeophyllum trabeum. We found genome-wide prevalence of polycistronic transcription in these Agaricomycetes, involving up to 8% of the transcribed genes. Unlike polycistronic mRNAs in prokaryotes, these co-transcribed genes are also independently transcribed. We show that polycistronic transcription may interfere with expression of the downstream tandem gene. Further comparative genomic analysis indicates that polycistronic transcription is conserved among a wide range of mushroom forming fungi. In summary, our study revealed, for the first time, the genome prevalence of polycistronic transcription in a phylogenetic range of higher fungi. Furthermore, we systematically show that our long-read sequencing approach and combined bioinformatics pipeline is a generic powerful tool for precise characterization of complex transcriptomes that enables identification of mRNA isoforms not recovered via short-read assembly.

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.

Localizing gene regulation reveals a staggered wood decay mechanism for the brown rot fungus Postia placenta

Significance Wood-decomposing fungi are key players in the carbon cycle and are models for making energy from lignocellulose, sustainably. Our study focuses on brown rot fungi that selectively remove carbohydrates, leaving most lignin behind. These fungi often decompose wood faster than their lignin-degrading white rot ancestors, despite losses in genes involved in plant cell wall hydrolysis. To explain brown rot, many have implicated reactive oxygen species (ROS) in facilitating hydrolysis, with microenvironmental gradients partitioning ROS from enzymes. By spatially colocalizing gene expression and enzyme activities as Postia placenta colonizes wood, we provide evidence of an oxidative-hydrolytic two-step mechanism controlled by differential expression, not microenvironments, and we highlight 549 genes (∼4% of the genome) that are upregulated during this unique pretreatment. Wood-degrading brown rot fungi are essential recyclers of plant biomass in forest ecosystems. Their efficient cellulolytic systems, which have potential biotechnological applications, apparently depend on a combination of two mechanisms: lignocellulose oxidation (LOX) by reactive oxygen species (ROS) and polysaccharide hydrolysis by a limited set of glycoside hydrolases (GHs). Given that ROS are strongly oxidizing and nonselective, these two steps are likely segregated. A common hypothesis has been that brown rot fungi use a concentration gradient of chelated metal ions to confine ROS generation inside wood cell walls before enzymes can infiltrate. We examined an alternative: that LOX components involved in ROS production are differentially expressed by brown rot fungi ahead of GH components. We used spatial mapping to resolve a temporal sequence in Postia placenta, sectioning thin wood wafers colonized directionally. Among sections, we measured gene expression by whole-transcriptome shotgun sequencing (RNA-seq) and assayed relevant enzyme activities. We found a marked pattern of LOX up-regulation in a narrow (5-mm, 48-h) zone at the hyphal front, which included many genes likely involved in ROS generation. Up-regulation of GH5 endoglucanases and many other GHs clearly occurred later, behind the hyphal front, with the notable exceptions of two likely expansins and a GH28 pectinase. Our results support a staggered mechanism for brown rot that is controlled by differential expression rather than microenvironmental gradients. This mechanism likely results in an oxidative pretreatment of lignocellulose, possibly facilitated by expansin- and pectinase-assisted cell wall swelling, before cellulases and hemicellulases are deployed for polysaccharide depolymerization.

Host-Selective Toxins of Pyrenophora tritici-repentis Induce Common Responses Associated with Host Susceptibility

Pyrenophora tritici-repentis (Ptr), a necrotrophic fungus and the causal agent of tan spot of wheat, produces one or a combination of host-selective toxins (HSTs) necessary for disease development. The two most studied toxins produced by Ptr, Ptr ToxA (ToxA) and Ptr ToxB (ToxB), are proteins that cause necrotic or chlorotic symptoms respectively. Investigation of host responses induced by HSTs provides better insight into the nature of the host susceptibility. Microarray analysis of ToxA has provided evidence that it can elicit responses similar to those associated with defense. In order to evaluate whether there are consistent host responses associated with susceptibility, a similar analysis of ToxB-induced changes in the same sensitive cultivar was conducted. Comparative analysis of ToxA- and ToxB-induced transcriptional changes showed that similar groups of genes encoding WRKY transcription factors, RLKs, PRs, components of the phenylpropanoid and jasmonic acid pathways are activated. ROS accumulation and photosystem dysfunction proved to be common mechanism-of-action for these toxins. Despite similarities in defense responses, transcriptional and biochemical responses as well as symptom development occur more rapidly for ToxA compared to ToxB, which could be explained by differences in perception as well as by differences in activation of a specific process, for example, ethylene biosynthesis in ToxA treatment. Results of this study suggest that perception of HSTs will result in activation of defense responses as part of a susceptible interaction and further supports the hypothesis that necrotrophic fungi exploit defense responses in order to induce cell death.

Solution NMR Structures of Pyrenophora tritici-repentis ToxB and Its Inactive Homolog Reveal Potential Determinants of Toxin Activity

Background: ToxB is a proteinaceous toxin but its homolog toxb has no toxic activity. Results: Both adopt a β-sandwich fold stabilized by two disulfide bonds but differ in the dynamics of one sandwich half. Conclusion: Toxicity is correlated with decreased compactness, increased flexibility, and polymorphism in an active site loop. Significance: ToxB activity depends on interplay between internal dynamics and interactions with putative targets. Pyrenophora tritici-repentis Ptr ToxB (ToxB) is a proteinaceous host-selective toxin produced by Pyrenophora tritici-repentis (P. tritici-repentis), a plant pathogenic fungus that causes the disease tan spot of wheat. One feature that distinguishes ToxB from other host-selective toxins is that it has naturally occurring homologs in non-pathogenic P. tritici-repentis isolates that lack toxic activity. There are no high-resolution structures for any of the ToxB homologs, or for any protein with >30% sequence identity, and therefore what underlies activity remains an open question. Here, we present the NMR structures of ToxB and its inactive homolog Ptr toxb. Both proteins adopt a β-sandwich fold comprising three strands in each half that are bridged together by two disulfide bonds. The inactive toxb, however, shows higher flexibility localized to the sequence-divergent β-sandwich half. The absence of toxic activity is attributed to a more open structure in the vicinity of one disulfide bond, higher flexibility, and residue differences in an exposed loop that likely impacts interaction with putative targets. We propose that activity is regulated by perturbations in a putative active site loop and changes in dynamics distant from the site of activity. Interestingly, the new structures identify AvrPiz-t, a secreted avirulence protein produced by the rice blast fungus, as a structural homolog to ToxB. This homology suggests that fungal proteins involved in either disease susceptibility such as ToxB or resistance such as AvrPiz-t may have a common evolutionary origin.

Extensive gene content variation in the Brachypodium distachyon pan-genome correlates with population structure

While prokaryotic pan-genomes have been shown to contain many more genes than any individual organism, the prevalence and functional significance of differentially present genes in eukaryotes remains poorly understood. Whole-genome de novo assembly and annotation of 54 lines of the grass Brachypodium distachyon yield a pan-genome containing nearly twice the number of genes found in any individual genome. Genes present in all lines are enriched for essential biological functions, while genes present in only some lines are enriched for conditionally beneficial functions (e.g., defense and development), display faster evolutionary rates, lie closer to transposable elements and are less likely to be syntenic with orthologous genes in other grasses. Our data suggest that differentially present genes contribute substantially to phenotypic variation within a eukaryote species, these genes have a major influence in population genetics, and transposable elements play a key role in pan-genome evolution.The role of differential gene content in the evolution and function of eukaryotic genomes remains poorly explored. Here the authors assemble and annotate the Brachypodium distachyon pan-genome consisting of 54 diverse lines and reveal the differential present genes as a major driver of phenotypic variation.

Changing the Game: Using Integrative Genomics to Probe Virulence Mechanisms of the Stem Rust Pathogen Puccinia graminis f. sp. tritici

The recent resurgence of wheat stem rust caused by new virulent races of Puccinia graminis f. sp. tritici (Pgt) poses a threat to food security. These concerns have catalyzed an extensive global effort toward controlling this disease. Substantial research and breeding programs target the identification and introduction of new stem rust resistance (Sr) genes in cultivars for genetic protection against the disease. Such resistance genes typically encode immune receptor proteins that recognize specific components of the pathogen, known as avirulence (Avr) proteins. A significant drawback to deploying cultivars with single Sr genes is that they are often overcome by evolution of the pathogen to escape recognition through alterations in Avr genes. Thus, a key element in achieving durable rust control is the deployment of multiple effective Sr genes in combination, either through conventional breeding or transgenic approaches, to minimize the risk of resistance breakdown. In this situation, evolution of pathogen virulence would require changes in multiple Avr genes in order to bypass recognition. However, choosing the optimal Sr gene combinations to deploy is a challenge that requires detailed knowledge of the pathogen Avr genes with which they interact and the virulence phenotypes of Pgt existing in nature. Identifying specific Avr genes from Pgt will provide screening tools to enhance pathogen virulence monitoring, assess heterozygosity and propensity for mutation in pathogen populations, and confirm individual Sr gene functions in crop varieties carrying multiple effective resistance genes. Toward this goal, much progress has been made in assembling a high quality reference genome sequence for Pgt, as well as a Pan-genome encompassing variation between multiple field isolates with diverse virulence spectra. In turn this has allowed prediction of Pgt effector gene candidates based on known features of Avr genes in other plant pathogens, including the related flax rust fungus. Upregulation of gene expression in haustoria and evidence for diversifying selection are two useful parameters to identify candidate Avr genes. Recently, we have also applied machine learning approaches to agnostically predict candidate effectors. Here, we review progress in stem rust pathogenomics and approaches currently underway to identify Avr genes recognized by wheat Sr genes.

An overview of genetic rust resistance: From broad to specific mechanisms

1 Commonwealth Scientific and Industrial Research Organisation (CSIRO), Agriculture and Food, Canberra, Australian Capital Territory, Australia, 2 Research School of Biology, The Australian National University, Canberra Australian Capital Territory, Australia, 3 Department of Plant Pathology and The Stakman-Borlaug Center for Sustainable Plant Health, University of Minnesota, St. Paul, Minnesota, United States of America