Chromatin Biology Impacts Adaptive Evolution of Filamentous Plant Pathogens

Abstract

In their adaptive response, pathogens evolve strategies, often involving secreted effector molecules, to overcome host immunity and support host colonization [3]. [...]it can be anticipated that this coevolutionary arms race leads to highly specific interactions between adapted pathogens and their specific hosts. Additionally, chromatin structure also plays crucial regulatory roles in establishing symbiotic interaction between the fungus Epichloë festucae and its plant host [27]. [...]studying the impact of chromatin biology on genome organization can broaden our knowledge and potentially provide mechanistic understanding of the evolution of two-speed genomes in plant pathogens and, in general, of adaptive genome evolution in plant-fungus interactions.

About the Authors:

Michael F. Seidl

Contributed equally to this work with: Michael F. Seidl, David E. Cook

Affiliation: Laboratory of Phytopathology, Wageningen University, The Netherlands

David E. Cook

Contributed equally to this work with: Michael F. Seidl, David E. Cook

Affiliation: Laboratory of Phytopathology, Wageningen University, The Netherlands

Bart P. H. J. Thomma

* E-mail: [email protected]

Affiliation: Laboratory of Phytopathology, Wageningen University, The NetherlandsCitation: Seidl MF, Cook DE, Thomma BPHJ (2016) Chromatin Biology Impacts Adaptive Evolution of Filamentous Plant Pathogens. PLoS Pathog 12(11): e1005920. doi:10.1371/journal.ppat.1005920

Editor: Cyril Zipfel, THE SAINSBURY LABORATORY, UNITED KINGDOM

Published: November 3, 2016

Copyright: © 2016 Seidl et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Funding: MFS acknowledges the receipt of a VENI grant (863.15.005) of the Research Council Earth and Life Sciences (ALW) of the Netherlands Organization of Scientific Research (NWO). Funding was provided to DEC by the Human Frontiers in Science Program (HFSP) Long-term fellowship (LT000627/2014-L). The work in the laboratory of BPJHT is supported by ALW-NWO. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

An organism's adaptation to changing environments is fueled by its genetic variability, which is established by mechanisms ranging from single-nucleotide polymorphisms to large-scale structural variations, all of which affect chromosomal shape, organization, and gene content [1]. These processes are particularly relevant for pathogens that must respond to continual selection pressure arising from host immune systems that evolved to detect the presence or activity of potential microbial pathogens through a variety of invasion patterns [2]. In their adaptive response, pathogens evolve strategies, often involving secreted effector molecules, to overcome host immunity and support host colonization [3]. Thus, it can be anticipated that this coevolutionary arms race leads to highly specific interactions between adapted pathogens and their specific hosts. Paradoxically, particular pathogens successfully colonize a broad range of hosts, yet how such pathogens cope in arms races with such a diversity of hosts remains unknown.

A structured genome drives adaptive evolution

It has been proposed that filamentous fungal and oomycete plant pathogens often evolved structured genomes with two distinct types of genomic regions: (1) gene-rich regions containing slowly evolving genes that mediate general physiology and (2) gene-poor regions that are dynamic and enriched for repetitive DNA such as transposable elements (TEs) and virulence-related genes that often display signs of accelerated evolution [4,5]. Extensive structural variation often occurs in these fast-evolving regions, leading to translocation, duplication, or loss of genetic material [1,4,5]. The highly dynamic regions can either be embedded within the core chromosomes or be located on separate chromosomes that often display presence/absence variations within a population, known as conditionally dispensable or accessory chromosomes [4,6]. The common occurrence of these bipartite genomes led to the "two-speed" model for pathogen genome evolution [4], suggesting that specific regions form sites of rapid genomic diversification to facilitate coevolution during host interactions [1,4,5].

Transposable elements shape "two-speed" genomes

It is generally observed that the dynamic regions of two-speed genomes are enriched in TEs, yet it remains undemonstrated how they mechanistically contribute to genome variability [1,4,5]. Recently, the contribution of TEs towards the evolution of the two-speed genome in the vascular wilt pathogen Verticillium dahliae was reported [7]. Extensive genome rearrangements are generated by double-strand repair pathways that erroneously utilize stretches of homologous sequence at an unlinked locus, the majority of which occur around TEs simply as a consequence of their abundance and sequence similarity (Fig 1A) [7]. Genomic rearrangements often occur at dynamic regions that are enriched for recent segmental duplications, generating genetic material subject to evolutionary diversification, e.g., by reciprocal gene loss (Fig 1A) [7,8]. Furthermore, these dynamic regions are enriched in in planta induced effectors [8] and evolutionary young and "active" TEs (Fig 1B). These TEs likely contribute to accelerated genomic diversification of dynamic effector regions [7].

Download:

*

PPT

PowerPoint slide

*

PNG

larger image

*

TIFF

original image

Fig 1. Dynamic genomic regions are associated with transposons and with distinct chromatin landscapes.

(A) In V. dahliae strain JR2, repeat-rich dynamic effector regions that evolve by genome rearrangements (indicated by red arrow heads) and by extensive segmental duplications (links between highly similar duplicated regions shown in grey) are located on chromosomes 2, 4, and 5 (red highlights). Repeat density along the chromosomes is indicated by a pink line (summarized as percent nucleotides in genomic windows of 5 kb, with a slide of 0.5 kb). (B) Dynamic genomic regions in V. dahliae are enriched in transcriptionally "active" and evolutionary "young" repetitive elements when compared with the core genome [7]. (C) Different histone modifications can be associated with core (Chr 13) and conditionally dispensable (Chr 14) chromosomes in the wheat pathogen Zymoseptoria tritici (as previously reported [13]). Repeat density along the chromosomes is indicated by a pink line (summarized as percent nucleotide in genomic windows of 5 kb, with a slide of 0.5 kb). Publicly available chromatin immunoprecipitation sequencing (ChIP-seq) samples [13] were mapped to the Z. tritici genome, and enriched regions were subsequently identified using RSEG [28]. DNA associated with euchromatic (H3K4me2, green) and heterochromatic (H3K27me3, orange; H3K9me3, red) marks are indicated, and significantly enriched genomic regions are marked with a solid line. Structural variations (duplications, black; deletions, blue) were identified by CNVnator, using publicly available resequencing data of multiple Z. tritici strains [29]. (D) The number of nucleotides (in Mb) of the Z. tritici genome covered by different histone regions (defined by RSEG) for euchromatin (green) and heterochromatin (orange, red) are shown by bar charts. (E) The number of duplications and deletions overlapping with histone regions (see above) are shown by bar charts.

http://dx.doi.org/10.1371/journal.ppat.1005920.g001

Chromatin biology impacts the adaptive evolution of filamentous plant pathogens

Chromatin, a complex of nucleic acids and proteins, determines the physical shape and organization of DNA within the nucleus [9]. In many eukaryotes, highly repetitive regions are composed of tightly condensed chromatin, referred to as heterochromatin, as opposed to open chromatin or euchromatin. Heterochromatin functions to silence repetitive and neighboring DNA and to repress recombination in many eukaryotic genomes [10]. However, these observations are inconsistent with data from fungal plant pathogens since repeat-rich regions of the two-speed genome are enriched in structural variations [7,11], and genes located in these regions often show highly coordinated expression [4]. To address this conundrum, and to understand the formation, maintenance, and transcriptional regulation of bipartite genomes, it is necessary to analyze genome structure and organization through the study of chromatin biology.

Chromatin and genome organization

Chromatin features such as accessibility, histone modifications, and chromosome conformation contribute to the organization of a genome. Recent modeling of mammalian and yeast genomes suggests that genomic rearrangements can be accurately predicted based on two chromatin features alone; chromatin openness and DNA contact in the nucleus (Fig 2A) [12]. Chromatin studies in most fungi remain scarce, but recent work on the fungal wheat pathogen Z. tritici shows marked differences in histone modifications between core and conditionally dispensable chromosomes (Fig 1C) [13]. Only half of the genome is occupied by histones carrying modifications that are commonly associated with heterochromatic regions, such as histone methylation of lysine residues on histone 3 at position 9 or 27 (H3K9me3 or H3K27me3). Notably, the majority of the detected structural variations colocalize with these heterochromatic regions (Fig 1C-1E). This suggests a further link between chromatin and genome structure, and despite the general thought that heterochromatin suppresses genomic variation, heterochromatic regions appear highly variable in the plant pathogenic fungi analyzed to date (see [6] for additional examples).

Download:

*

PPT

PowerPoint slide

*

PNG

larger image

*

TIFF

original image

Fig 2. Impact of chromatin organization on adaptive (genome) evolution.

(A) Genome rearrangements (red arrows) occur in open chromatin regions (euchromatin, green background; heterochromatin, red background) that are in spatial proximity within the nucleus. Chromosomes are shown as lines and genes as differentially shaped symbols. Spatial organization of the nucleus is highlighted by nuclear regions that are occupied by distinct chromosomes (different colors). (B) Chromatin influences gene expression, as genes located in open chromatin are transcribed (arrows), while genes located in heterochromatin are transcriptionally silent. (C) Interspecific genome hybridization leads to genome restructuring, often accompanied by extensive gene loss, and changes in transcriptional profiles, which can be influenced by differences in chromatin between parental species (genes from two parental lineages are indicated by orange and blue, respectively). (D) Host jumps, specialization towards a specific host, and evasion of host immunity can be influenced by changes in chromatin that translate to alterations in gene expression or DNA content.

http://dx.doi.org/10.1371/journal.ppat.1005920.g002

A possible mechanism to explain the link between chromatin and genome structure is that heterochromatic regions are more prone to DNA breaks during replication that, when repaired, could result in structural variations [1,14]. In Neurospora crassa, heterochromatin is highly enriched for a specific modification of histone 2, [gamma]H2A ([gamma]H2A.X in mammals), which generally marks double-strand DNA breaks [15]. It is not known why this histone modification marks these regions, but if it follows its canonical role, these regions showing elevated rates of DNA breaks would be prone to erroneous replication-based repair. Additionally, N. crassa chromosome conformation maps indicate frequent and strong contact between heterochromatic regions [16]. Taken together, the increased rates of genomic variation at heterochromatic regions may arise from DNA damage during the replication of heterochromatic DNA, followed by error-prone replication-based repair or template switching between contacting loci.

Chromatin and gene expression

The majority of chromatin studies using plant-pathogenic fungi have focused on their transcriptional impacts (Fig 2B) [17-19]. Canonical repressive chromatin marks, such as H3K9me2/3 or H3K27me2/3, are significantly enriched at regions known to harbor genes important for pathogenicity [17,18]. Additionally, the canonical activating mark H3K4me2/3 was shown to play a significant role in regulating gene expression and promoting growth and virulence in Magnaporthe oryzae [19]. Collectively, these studies indicate that deregulating chromatin significantly impacts fungal transcription, growth, and development. To gain further understanding into the regulation of pathogenicity, it is necessary to determine the key enzymes responsible for reading, writing, and erasing these epigenetic modifications.

Chromatin and adaptation

Chromatin modifications can functionally diversify between species. For example, the addition of methyl groups to histone 3 at lysine 36 (H3K36me3) functions in guiding both DNA mismatch repair and messenger RNA (mRNA) splicing machinery in humans [20,21], while in the yeast Saccharomyces cerevisiae, the same mark functions to suppress intragenic transcripts [22]. Thus, although the overall localization of the mark at transcribed genes is conserved, species evolved to utilize the mark differently. Along with functional diversification between species, chromatin architecture can vary within a population [23]. It has been proposed that environmental conditions can result in heritable epigenetic variation and act as a substrate for so-called "Neo-Darwin" selection [24].

Given that chromatin regulation can evolve, is variable within populations, and can be environmentally influenced, it can be anticipated that chromatin affects host-microbe interactions. For example, several documented cases show that plant pathogens may undergo, or are the result of, interspecific genome-hybridization events [25]. These events are accompanied by complex gene expression changes, likely influenced by parental chromatin structure, and by genome reorganization in the hybrids (Fig 2C). These changes in expression and DNA content can lead to alterations in the hybrid's ability to compete for a particular niche and could facilitate a host jump (Fig 2D). Pathogens that are well adapted to a specific host may also undergo chromatin-based regulatory changes that alter their host interaction. For example, expression of the Avr3a effector in Phytophthora sojae can change between generations, allowing for evasion of recognition by the corresponding host immune receptor Rps3a in soybean [26]. The mechanism for this switch has not been reported, but experimental evidence shows an enrichment of small RNA (sRNA) at the locus in silenced versus Avr3a-expressing lines with no genetic mutations reported, suggesting epigenetic regulation (Fig 2D). The prevalence of such an epigenetic evasion strategy in other plant-microbe interactions remains unknown. It is also interesting to speculate that variation in chromatin structure within a pathogen population could influence a broad-host-range pathogen to form strains with increased fitness on a particular host. Conceptually, this could be achieved through alterations in the chromatin-based regulation for the optimal timing, rate, and order of transcriptional events leading to successful infection. That is, alterations in transcription to dampen or augment an infection strategy better suited for one host versus another could influence the evolution of a plant-microbe relationship (Fig 2D). Future experiments to address these possibilities may uncover yet additional layers influencing the evolution of plant-microbe interactions.

Outlook

Chromatin biology can impact filamentous pathogens across spatial and temporal scales, from governing genome organization to controlling the expression of its individual parts. We conceive that complex chromatin structures in pathogenic fungi will influence not only coordinated effector expression in dynamic genomic regions but also structural variations, thereby further linking genome and chromatin structure to genome evolution and adaptation. Additionally, chromatin structure also plays crucial regulatory roles in establishing symbiotic interaction between the fungus Epichloë festucae and its plant host [27]. Thus, studying the impact of chromatin biology on genome organization can broaden our knowledge and potentially provide mechanistic understanding of the evolution of two-speed genomes in plant pathogens and, in general, of adaptive genome evolution in plant-fungus interactions.

References

Seidl MF, Thomma BPHJ. Sex or no sex: evolutionary adaptation occurs regardless. Bioessays. 2014;36: 335-345 doi: 10.1002/bies.201300155. pmid:24531982

Cook DE, Mesarich CH, Thomma BPHJ. Understanding plant immunity as a surveillance system to detect invasion. Annu Rev Phytopathol. 2015;53: 541-563 doi: 10.1146/annurev-phyto-080614-120114. pmid:26047564

Rovenich H, Boshoven JC, Thomma BPHJ. Filamentous pathogen effector functions: of pathogens, hosts and microbiomes. Curr Opin Plant Biol. 2014;20C: 96-103

Raffaele S, Kamoun S. Genome evolution in filamentous plant pathogens: why bigger can be better. Nat Rev Microbiol. 2012;10: 417-430 doi: 10.1038/nrmicro2790. pmid:22565130

Dong S, Raffaele S, Kamoun S. The two-speed genomes of filamentous pathogens: waltz with plants. Curr Opin Genet Dev. 2015;35: 57-65 doi: 10.1016/j.gde.2015.09.001. pmid:26451981

Galazka JM, Freitag M. Variability of chromosome structure in pathogenic fungi-of "ends and odds." Curr Opin Microbiol. 2014;20C: 19-26

Faino L, Seidl MF, Shi-Kunne X, Pauper M, van den Berg GCM, Wittenberg AHJ, et al. Transposons passively and actively contribute to evolution of the two-speed genome of a fungal pathogen. Genome Res. 2016;26: 1091-1100 doi: 10.1101/gr.204974.116. pmid:27325116

de Jonge R, Bolton MD, Kombrink A, van den Berg GCM, Yadeta KA, Thomma BPHJ. Extensive chromosomal reshuffling drives evolution of virulence in an asexual pathogen. Genome Res. 2013;23: 1271-1282 doi: 10.1101/gr.152660.112. pmid:23685541

Friedman N, Rando OJ. Epigenomics and the structure of the living genome. Genome Res. 2015;25: 1482-1490 doi: 10.1101/gr.190165.115. pmid:26430158

Grewal SIS, Jia S. Heterochromatin revisited. Nat Rev Genet. 2007;8: 35-46 doi: 10.1038/nrg2008. pmid:17173056

Raffaele S, Farrer RA, Cano LM, Studholme DJ, MacLean D, Thines M, et al. Genome evolution following host jumps in the Irish potato famine pathogen lineage. Science. 2010;330: 1540-1543 doi: 10.1126/science.1193070. pmid:21148391

Berthelot C, Muffato M, Abecassis J, Roest Crollius H. The 3D organization of chromatin explains evolutionary fragile genomic regions. Cell Rep. 2015;10: 1913-1924 doi: 10.1016/j.celrep.2015.02.046. pmid:25801028

Schotanus K, Soyer JL, Connolly LR, Grandaubert J, Happel P, Smith KM, et al. Histone modifications rather than the novel regional centromeres of Zymoseptoria tritici distinguish core and accessory chromosomes. Epigenetics Chromatin. 2015;8: 41 doi: 10.1186/s13072-015-0033-5. pmid:26430472

Carvalho CMB, Lupski JR. Mechanisms underlying structural variant formation in genomic disorders. Nat Rev Genet. 2016;17: 224-238 doi: 10.1038/nrg.2015.25. pmid:26924765

Sasaki T, Lynch KL, Mueller CV, Friedman S, Freitag M, Lewis ZA. Heterochromatin controls [gamma]H2A localization in Neurospora crassa. Eukaryotic Cell. 2014;13: 990-1000 doi: 10.1128/EC.00117-14. pmid:24879124

Galazka JM, Klocko AD, Uesake M, Honda S, Selker EU, Freitag M. Neurospora chromosomes are organized by blocs of importin alpha-dependent heterochromatin that are largely independent of H3K9me3. Genome Res. 2016; gr.203182.115

Connolly LR, Smith KM, Freitag M. The Fusarium graminearum histone H3 K27 methyltransferase KMT6 regulates development and expression of secondary metabolite gene clusters. PLoS Genet. 2013;9: e1003916 doi: 10.1371/journal.pgen.1003916. pmid:24204317

Soyer JL, Ghalid El M, Glaser N, Ollivier B, Linglin J, Grandaubert J, et al. Epigenetic control of effector gene expression in the plant pathogenic fungus Leptosphaeria maculans. PLoS Genet. 2014;10: e1004227. doi: 10.1371/journal.pgen.1004227. pmid:24603691

Pham KTM, Inoue Y, Vu BV, Nguyen HH, Nakayashiki T, Ikeda K-I, et al. MoSET1 (Histone H3K4 Methyltransferase in Magnaporthe oryzae) regulates global gene expression during infection-related morphogenesis. PLoS Genet. 2015;11: e1005385. doi: 10.1371/journal.pgen.1005385. pmid:26230995

Luco RF, Pan Q, Tominaga K, Blencowe BJ, Pereira-Smith OM, Misteli T. Regulation of alternative splicing by histone modifications. Science. 2010;327: 996-1000 doi: 10.1126/science.1184208. pmid:20133523

Li F, Mao G, Tong D, Huang J, Gu L, Yang W, et al. The histone mark H3K36me3 regulates human DNA mismatch repair through its interaction with MutS[alpha] . Cell. 2013;153: 590-600 doi: 10.1016/j.cell.2013.03.025. pmid:23622243

Carrozza MJ, Li B, Florens L, Suganuma T, Swanson SK, Lee KK, et al. Histone H3 methylation by Set2 directs deacetylation of coding regions by Rpd3S to suppress spurious intragenic transcription. Cell. 2005;123: 581-592 doi: 10.1016/j.cell.2005.10.023. pmid:16286007

Waszak SM, Delaneau O, Gschwind AR, Kilpinen H, Raghav SK, Witwicki RM, et al. Population variation and genetic control of modular chromatin architecture in humans. Cell. 2015;162: 1039-1050 doi: 10.1016/j.cell.2015.08.001. pmid:26300124

Skinner MK. Environmental epigenetics and a unified theory of the molecular aspects of evolution: A Neo-Lamarckian concept that facilitates Neo-Darwinian evolution. Genome Biol Evol. 2015;7: 1296-1302 doi: 10.1093/gbe/evv073. pmid:25917417

Depotter JR, Seidl MF, Wood TA, Thomma BPHJ. Interspecific hybridization impacts host range and pathogenicity of filamentous microbes. Curr Opin Microbiol. 2016;32: 7-13 doi: 10.1016/j.mib.2016.04.005. pmid:27116367

Qutob D, Chapman BP, Gijzen M. Transgenerational gene silencing causes gain of virulence in a plant pathogen. Nat Commun. 2013;4: 1349 doi: 10.1038/ncomms2354. pmid:23322037

Chujo T, Scott B. Histone H3K9 and H3K27 methylation regulates fungal alkaloid biosynthesis in a fungal endophyte-plant symbiosis. Mol Microbiol. 2014;92: 413-434 doi: 10.1111/mmi.12567. pmid:24571357

Song Q, Smith AD. Identifying dispersed epigenomic domains from ChIP-Seq data. Bioinformatics. 2011;27: 870-871 doi: 10.1093/bioinformatics/btr030. pmid:21325299

Croll D, Zala M, McDonald BA. Breakage-fusion-bridge cycles and large insertions contribute to the rapid evolution of accessory chromosomes in a fungal pathogen. PLoS Genet. 2013;9: e1003567 doi: 10.1371/journal.pgen.1003567. pmid:23785303