Destabilization of chromosome structure by histone H3 lysine 27 methylation

From Lesions to Lessons: Two Decades of Filamentous Plant Pathogen Genomics

Many filamentous microorganisms, such as fungi and oomycetes, have evolved the ability to colonize plants and cause devastating crop diseases. Coevolutionary conflicts with their hosts have shaped the genomes of these plant pathogens. Over the past 20 years, genomics and genomics-enabled technologies have revealed remarkable diversity in genome size, architecture, and gene regulatory mechanisms. Technical and conceptual advances continue to provide novel insights into evolutionary dynamics, diversification of distinct genomic compartments, and facilitated molecular disease diagnostics. In this review, we discuss how genomics has advanced our understanding of genome organization and plant-pathogen coevolution and provide a perspective on future developments in the field. [Formula: see text] Copyright © 2025 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

The Molecular Dialogue Between Zymoseptoria tritici and Wheat

Zymoseptoria tritici is a highly damaging pathogen that causes high wheat yield losses in temperate climates. Z. tritici emerged during the domestication of wheat in the Fertile Crescent and has been extensively used as a model system for population genetic and genomic studies. New genetic tools and resources have provided a better understanding of the molecular components involved in the wheat-Z. tritici interaction, which is highlighted by the cloning of three wheat resistance genes and four Z. tritici avirulence genes. Despite the considerable progress made in the last few years, the mechanisms that mediate Z. tritici colonization remain largely unknown. In this review, we summarize the latest advances in understanding the molecular components mediating wheat-Z. tritici interactions, and we discuss future research lines to close current knowledge gaps. [Formula: see text] Copyright © 2025 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Comparative Analyses of Compatible and Incompatible Host-Pathogen Interactions Provide Insight into Divergent Host Specialization of Closely Related Pathogens

Host-pathogen co-evolutionary dynamics drive constant changes in plant pathogens to thrive in their plant host. Factors that determine host specificity are diverse and range from molecular and morphological strategies to metabolic and reproductive adaptations. We applied an experimental approach and conducted comparative microscopy, transcriptome analyses, and functional analyses of selected pathogen traits to identify determinants of host specificity in an important wheat pathogen. We included three closely related fungal pathogens, Zymoseptoria tritici, Z. pseudotritici, and Z. ardabiliae, that establish compatible and incompatible interactions with wheat. Although infections of the incompatible species induce plant defenses during invasion of stomatal openings, we found a conserved early-infection program among the three species whereby only 9.2% of the 8,885 orthologous genes are significantly differentially expressed during initial infection. The genes upregulated in Z. tritici likely reflect specialization to wheat, whereas upregulated genes in the incompatible interaction may reflect processes to counteract cellular stress associated with plant defenses. We selected nine candidate genes encoding putative effectors and host-specificity determinants in Z. tritici and deleted these to study their functional relevance. Despite the particular expression patterns of the nine genes, only two mutants were impaired in virulence. We further expressed the Z. tritici proteins in Nicotiana benthamiana to investigate protein function and assess cell death reaction. Hereby, we identify three effectors with cell-death-inducing properties. From the functional analyses, we conclude that the successful infection of Z. tritici in wheat relies on an extensive redundancy of virulence determinants. [Formula: see text] Copyright © 2025 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

KDM6A-SND1 interaction maintains genomic stability by protecting the nascent DNA and contributes to cancer chemoresistance

Genomic instability is one of the hallmarks of cancer. While loss of histone demethylase KDM6A increases the risk of tumorigenesis, its specific role in maintaining genomic stability remains poorly understood. Here, we propose a mechanism in which KDM6A maintains genomic stability independently on its demethylase activity. This occurs through its interaction with SND1, resulting in the establishment of a protective chromatin state that prevents replication fork collapse by recruiting of RPA and Ku70 to nascent DNA strand. Notably, KDM6A-SND1 interaction is up-regulated by KDM6A SUMOylation, while KDM6AK90A mutation almost abolish the interaction. Loss of KDM6A or SND1 leads to increased enrichment of H3K9ac and H4K8ac but attenuates the enrichment of Ku70 and H3K4me3 at nascent DNA strand. This subsequently results in enhanced cellular sensitivity to genotoxins and genomic instability. Consistent with these findings, knockdown of KDM6A and SND1 in esophageal squamous cell carcinoma (ESCC) cells increases genotoxin sensitivity. Intriguingly, KDM6A H101D & P110S, N1156T and D1216N mutations identified in ESCC patients promote genotoxin resistance via increased SND1 association. Our finding provides novel insights into the pivotal role of KDM6A-SND1 in genomic stability and chemoresistance, implying that targeting KDM6A and/or its interaction with SND1 may be a promising strategy to overcome the chemoresistance.

Segmental duplications drive the evolution of accessory regions in a major crop pathogen

Many pathogens evolved compartmentalized genomes with conserved core and variable accessory regions (ARs) that carry effector genes mediating virulence. The fungal plant pathogen Fusarium oxysporum has such ARs, often spanning entire chromosomes. The presence of specific ARs influences the host range, and horizontal transfer of ARs can modify the pathogenicity of the receiving strain. However, how these ARs evolve in strains that infect the same host remains largely unknown. We defined the pan-genome of 69 diverse F. oxysporum strains that cause Fusarium wilt of banana, a significant constraint to global banana production, and analyzed the diversity and evolution of the ARs. Accessory regions in F. oxysporum strains infecting the same banana cultivar are highly diverse, and we could not identify any shared genomic regions and in planta-induced effectors. We demonstrate that segmental duplications drive the evolution of ARs. Furthermore, we show that recent segmental duplications specifically in accessory chromosomes cause the expansion of ARs in F. oxysporum. Taken together, we conclude that extensive recent duplications drive the evolution of ARs in F. oxysporum, which contribute to the evolution of virulence.

Regulation of effector gene expression as concerted waves in Leptosphaeria maculans: a two-player game

Effector genes, encoding molecules involved in disease establishment, are concertedly expressed throughout the lifecycle of plant-pathogenic fungi. However, little is known about how effector gene expression is regulated. Since many effector genes are located in repeat-rich regions, the role of chromatin remodeling in their regulation was recently investigated, notably establishing that the repressive histone modification H3K9me3, deposited by KMT1, was involved in several fungal species including Leptosphaeria maculans. Nevertheless, previous data suggest that a second regulatory layer, probably involving a specific transcription factor (TF), might be required. In L. maculans, a Dothideomycete causing stem canker of oilseed rape, we identified the ortholog of Pf2, a TF belonging to the Zn2Cys6 fungal-specific family, and described as essential for pathogenicity and effector gene expression. We investigated its role together with KMT1, by inactivating and over-expressing LmPf2 in a wild-type strain and a ∆kmt1 mutant. Functional analyses of the corresponding transformants highlighted an essential role of LmPf2 in the establishment of pathogenesis and we found a major effect of LmPf2 on the induction of effector gene expression once KMT1 repression is lifted. Our results show, for the first time, a dual control of effector gene expression.

Implications of the three-dimensional chromatin organization for genome evolution in a fungal plant pathogen

The spatial organization of eukaryotic genomes is linked to their biological functions, although it is not clear how this impacts the overall evolution of a genome. Here, we uncover the three-dimensional (3D) genome organization of the phytopathogen Verticillium dahliae, known to possess distinct genomic regions, designated adaptive genomic regions (AGRs), enriched in transposable elements and genes that mediate host infection. Short-range DNA interactions form clear topologically associating domains (TADs) with gene-rich boundaries that show reduced levels of gene expression and reduced genomic variation. Intriguingly, TADs are less clearly insulated in AGRs than in the core genome. At a global scale, the genome contains bipartite long-range interactions, particularly enriched for AGRs and more generally containing segmental duplications. Notably, the patterns observed for V. dahliae are also present in other Verticillium species. Thus, our analysis links 3D genome organization to evolutionary features conserved throughout the Verticillium genus.

Recent reactivation of a pathogenicity-associated transposable element is associated with major chromosomal rearrangements in a fungal wheat pathogen

Transposable elements (TEs) are key drivers of genomic variation contributing to recent adaptation in most species. Yet, the evolutionary origins and insertion dynamics within species remain poorly understood. We recapitulate the spread of the pathogenicity-associated Styx element across five species that last diverged ∼11 000 years ago. We show that the element likely originated in the Zymoseptoria fungal pathogen genus and underwent multiple independent reactivation events. Using a global 900-genome panel of the wheat pathogen Zymoseptoria tritici, we assess Styx copy number variation and identify renewed transposition activity in Oceania and South America. We show that the element can mobilize to create additional Styx copies in a four-generation pedigree. Importantly, we find that new copies of the element are not affected by genomic defenses suggesting minimal control against the element. Styx copies are preferentially located in recombination breakpoints and likely triggered multiple types of large chromosomal rearrangements. Taken together, we establish the origin, diversification and reactivation of a highly active TE with likely major consequences for chromosomal integrity and the expression of disease.

H3K27me3 is vital for fungal development and secondary metabolite gene silencing, and substitutes for the loss of H3K9me3 in the plant pathogen Fusarium proliferatum

Facultative heterochromatin marked by histone H3 lysine 27 trimethylation (H3K27me3) is an important regulatory layer involved in secondary metabolite (SM) gene silencing and crucial for fungal development in the genus Fusarium. While this histone mark is essential in some (e.g., the rice pathogen Fusarium fujikuroi), it appears dispensable in other fusaria. Here, we show that deletion of FpKMT6 is detrimental but not lethal in the plant pathogen Fusarium proliferatum, a member of the Fusarium fujikuroi species complex (FFSC). Loss of FpKmt6 results in aberrant growth, and expression of a large set of previously H3K27me3-silenced genes is accompanied by increased H3K27 acetylation (H3K27ac) and an altered H3K36me3 pattern. Next, H3K9me3 patterns are affected in Δfpkmt6, indicating crosstalk between both heterochromatic marks that became even more obvious in a strain deleted for FpKMT1 encoding the H3K9-specific histone methyltransferase. In Δfpkmt1, all H3K9me3 marks present in the wild-type strain are replaced by H3K27me3, a finding that may explain the subtle phenotype of the Δfpkmt1 strain which stands in marked contrast to other filamentous fungi. A large proportion of SM-encoding genes is allocated with H3K27me3 in the wild-type strain and loss of H3K27me3 results in elevated expression of 49% of them. Interestingly, genes involved in the biosynthesis of the phytohormones gibberellins (GA) are among the most upregulated genes in Δfpkmt6. Although several FFSC members harbor GA biosynthetic genes, its production is largely restricted to F. fujikuroi, possibly outlining the distinct lifestyles of these notorious plant pathogens. We show that H3K27me3 is involved in GA gene silencing in F. proliferatum and at least one additional FFSC member, and thus, may serve as a regulatory layer for gene silencing under non-favoring conditions.

H4K20me3 is important for Ash1-mediated H3K36me3 and transcriptional silencing in facultative heterochromatin in a fungal pathogen

Facultative heterochromatin controls development and differentiation in many eukaryotes. In metazoans, plants, and many filamentous fungi, facultative heterochromatin is characterized by transcriptional repression and enrichment with nucleosomes that are trimethylated at histone H3 lysine 27 (H3K27me3). While loss of H3K27me3 results in derepression of transcriptional gene silencing in many species, additional up- and downstream layers of regulation are necessary to mediate control of transcription in chromosome regions enriched with H3K27me3. Here, we investigated the effects of one histone mark on histone H4, namely H4K20me3, in the fungus Zymoseptoria tritici, a globally important pathogen of wheat. Deletion of kmt5, the gene encoding the sole methyltransferase responsible for H4K20 methylation, resulted in global derepression of transcription, especially in regions of facultative heterochromatin. Derepression in the absence of H4K20me3 not only affected known genes but also a large number of novel, previously undetected transcripts generated from regions of facultative heterochromatin on accessory chromosomes. Transcriptional activation in kmt5 deletion strains was accompanied by a complete loss of Ash1-mediated H3K36me3 and chromatin reorganization affecting H3K27me3 and H3K4me2 distribution in regions of facultative heterochromatin. Strains with H4K20L, M or Q mutations in the single histone H4 gene of Z. tritici recapitulated these chromatin changes, suggesting that H4K20me3 is important for Ash1-mediated H3K36me3. The ∆kmt5 mutants we obtained were more sensitive to genotoxic stressors than wild type and both, ∆kmt5 and ∆ash1, showed greatly increased rates of accessory chromosome loss. Taken together, our results provide insights into an unsuspected mechanism involved in the assembly and maintenance of facultative heterochromatin.

Epigenetic regulation of nuclear processes in fungal plant pathogens

Through the association of protein complexes to DNA, the eukaryotic nuclear genome is broadly organized into open euchromatin that is accessible for enzymes acting on DNA and condensed heterochromatin that is inaccessible. Chemical and physical alterations to chromatin may impact its organization and functionality and are therefore important regulators of nuclear processes. Studies in various fungal plant pathogens have uncovered an association between chromatin organization and expression of in planta-induced genes that are important for pathogenicity. This review discusses chromatin-based regulation mechanisms as determined in the fungal plant pathogen Verticillium dahliae and relates the importance of epigenetic transcriptional regulation and other nuclear processes more broadly in fungal plant pathogens.

Nuclear genome organization in fungi: from gene folding to Rabl chromosomes

Comparative genomics has recently provided unprecedented insights into the biology and evolution of the fungal lineage. In the postgenomics era, a major research interest focuses now on detailing the functions of fungal genomes, i.e. how genomic information manifests into complex phenotypes. Emerging evidence across diverse eukaryotes has revealed that the organization of DNA within the nucleus is critically important. Here, we discuss the current knowledge on the fungal genome organization, from the association of chromosomes within the nucleus to topological structures at individual genes and the genetic factors required for this hierarchical organization. Chromosome conformation capture followed by high-throughput sequencing (Hi-C) has elucidated how fungal genomes are globally organized in Rabl configuration, in which centromere or telomere bundles are associated with opposite faces of the nuclear envelope. Further, fungal genomes are regionally organized into topologically associated domain-like (TAD-like) chromatin structures. We discuss how chromatin organization impacts the proper function of DNA-templated processes across the fungal genome. Nevertheless, this view is limited to a few fungal taxa given the paucity of fungal Hi-C experiments. We advocate for exploring genome organization across diverse fungal lineages to ensure the future understanding of the impact of nuclear organization on fungal genome function.

Repeat-Induced Point Mutation and Gene Conversion Coinciding with Heterochromatin Shape the Genome of a Plant-Pathogenic Fungus

Meiosis is associated with genetic changes in the genome-via recombination, gene conversion, and mutations. The occurrence of gene conversion and mutations during meiosis may further be influenced by the chromatin conformation, similar to the effect of the chromatin conformation on the mitotic mutation rate. To date, however, the exact distribution and type of meiosis-associated changes and the role of the chromatin conformation in this context are largely unexplored. Here, we determine recombination, gene conversion, and de novo mutations using whole-genome sequencing of all meiotic products of 23 individual meioses in Zymoseptoria tritici, an important pathogen of wheat. We confirm a high genome-wide recombination rate of 65 centimorgan (cM)/Mb and see higher recombination rates on the accessory compared to core chromosomes. A substantial fraction of 0.16% of all polymorphic markers was affected by gene conversions, showing a weak GC-bias and occurring at higher frequency in regions of constitutive heterochromatin, indicated by the histone modification H3K9me3. The de novo mutation rate associated with meiosis was approximately three orders of magnitude higher than the corresponding mitotic mutation rate. Importantly, repeat-induced point mutation (RIP), a fungal defense mechanism against duplicated sequences, is active in Z. tritici and responsible for the majority of these de novo meiotic mutations. Our results indicate that the genetic changes associated with meiosis are a major source of variability in the genome of an important plant pathogen and shape its evolutionary trajectory. IMPORTANCE The impact of meiosis on the genome composition via gene conversion and mutations is mostly poorly understood, in particular, for non-model species. Here, we sequenced all four meiotic products for 23 individual meioses and determined the genetic changes caused by meiosis for the important fungal wheat pathogen Zymoseptoria tritici. We found a high rate of gene conversions and an effect of the chromatin conformation on gene conversion rates. Higher conversion rates were found in regions enriched with the H3K9me3-a mark for constitutive heterochromatin. Most importantly, meiosis was associated with a much higher frequency of de novo mutations than mitosis; 78% of the meiotic mutations were caused by repeat-induced point mutations-a fungal defense mechanism against duplicated sequences. In conclusion, the genetic changes associated with meiosis are therefore a major factor shaping the genome of this fungal pathogen.

Biosynthesis of the Isocoumarin Derivatives Fusamarins is Mediated by the PKS8 Gene Cluster in Fusarium

Fusarium mangiferae causes the mango malformation disease (MMD) on young mango trees and seedlings resulting in economically significant crop losses. In addition, F. mangiferae produces a vast array of secondary metabolites (SMs), including mycotoxins that may contaminate the harvest. Their production is tightly regulated at the transcriptional level. Here, we show that lack of the H3 K9-specific histone methyltransferase, FmKmt1, influences the expression of the F. mangiferae polyketide synthase (PKS) 8 (FmPKS8), a so far cryptic PKS. By a combination of reverse genetics, untargeted metabolomics, bioinformatics and chemical analyses including structural elucidation, we determined the FmPKS8 biosynthetic gene cluster (BGC) and linked its activity to the production of fusamarins (FMN), which can be structurally classified as dihydroisocoumarins. Functional characterization of the four FMN cluster genes shed light on the biosynthetic pathway. Cytotoxicity assays revealed moderate toxicities with IC50 values between 1 and 50 μM depending on the compound.

Potential and limits of (mal)adaptive mutation rate plasticity in plants

Genetic mutations provide the heritable material for plant adaptation to their environments. At the same time, the environment can affect the mutation rate across plant genomes. However, the extent to which environmental plasticity in mutation rates can facilitate or hinder adaptation remains a longstanding and unresolved question. Emerging discoveries of mechanisms affecting mutation rate variability provide opportunities to consider this question in a new light. Links between chromatin states, transposable elements, and DNA repair suggest cases of adaptive mutation rate plasticity could occur. Yet, numerous evolutionary and biological forces are expected to limit the impact of any such mutation rate plasticity on adaptive evolution. Persistent uncertainty about the significance of mutation rate plasticity on adaptation motivates new experimental and theoretical research relevant to understanding plant responses in changing environments.

A thousand-genome panel retraces the global spread and adaptation of a major fungal crop pathogen

Human activity impacts the evolutionary trajectories of many species worldwide. Global trade of agricultural goods contributes to the dispersal of pathogens reshaping their genetic makeup and providing opportunities for virulence gains. Understanding how pathogens surmount control strategies and cope with new climates is crucial to predicting the future impact of crop pathogens. Here, we address this by assembling a global thousand-genome panel of Zymoseptoria tritici, a major fungal pathogen of wheat reported in all production areas worldwide. We identify the global invasion routes and ongoing genetic exchange of the pathogen among wheat-growing regions. We find that the global expansion was accompanied by increased activity of transposable elements and weakened genomic defenses. Finally, we find significant standing variation for adaptation to new climates encountered during the global spread. Our work shows how large population genomic panels enable deep insights into the evolutionary trajectory of a major crop pathogen.

When ferroptosis meets pathogenic infections

Apoptosis, necrosis, or autophagy are diverse types of regulated cell death (RCD), recognized as the strategies that host cells use to defend against pathogens such as viruses, bacteria, or fungi. Pathogens can induce or block different types of host cell RCD, promoting propagation or evading host immune surveillance. Ferroptosis is a newly identified RCD. Evidence has demonstrated how pathogens regulate ferroptosis to promote their replication, dissemination, and pathogenesis. However, the interaction between ferroptosis and pathogenic infections still needs to be completely elucidated. This review summarizes the advances in the interaction between pathogenic infections and host ferroptotic processes, focusing on the underlying mechanisms of how pathogens exploit ferroptosis, and discussing possible therapeutic measures against pathogen-associated diseases in a ferroptosis-dependent manner.

Heterochromatin in the fungal plant pathogen, Zymoseptoria tritici: Control of transposable elements, genome plasticity and virulence

Heterochromatin is a repressive chromatin state that plays key roles in the functional organisation of eukaryotic genomes. In fungal plant pathogens, effector genes that are required for host colonization tend to be associated with heterochromatic regions of the genome that are enriched with transposable elements. It has been proposed that the heterochromatin environment silences effector genes in the absence of host and dynamic chromatin remodelling facilitates their expression during infection. Here we discuss this model in the context of the key wheat pathogen, Zymoseptoria tritici. We cover progress in understanding the deposition and recognition of heterochromatic histone post translational modifications in Z. tritici and the role that heterochromatin plays in control of genome plasticity and virulence.

The contribution of DNA repair pathways to genome editing and evolution in filamentous pathogens

DNA double-strand breaks require repair or risk corrupting the language of life. To ensure genome integrity and viability, multiple DNA double-strand break repair pathways function in eukaryotes. Two such repair pathways, canonical non-homologous end joining and homologous recombination, have been extensively studied, while other pathways such as microhomology-mediated end joint and single-strand annealing, once thought to serve as back-ups, now appear to play a fundamental role in DNA repair. Here, we review the molecular details and hierarchy of these four DNA repair pathways, and where possible, a comparison for what is known between animal and fungal models. We address the factors contributing to break repair pathway choice, and aim to explore our understanding and knowledge gaps regarding mechanisms and regulation in filamentous pathogens. We additionally discuss how DNA double-strand break repair pathways influence genome engineering results, including unexpected mutation outcomes. Finally, we review the concept of biased genome evolution in filamentous pathogens, and provide a model, termed Biased Variation, that links DNA double-strand break repair pathways with properties of genome evolution. Despite our extensive knowledge for this universal process, there remain many unanswered questions, for which the answers may improve genome engineering and our understanding of genome evolution.

Proteasome-dependent truncation of the negative heterochromatin regulator Epe1 mediates antifungal resistance

Epe1 histone demethylase restricts H3K9-methylation-dependent heterochromatin, preventing it from spreading over, and silencing, gene-containing regions in fission yeast. External stress induces an adaptive response allowing heterochromatin island formation that confers resistance on surviving wild-type lineages. Here we investigate the mechanism by which Epe1 is regulated in response to stress. Exposure to caffeine or antifungals results in Epe1 ubiquitylation and proteasome-dependent removal of the N-terminal 150 residues from Epe1, generating truncated Epe1 (tEpe1) which accumulates in the cytoplasm. Constitutive tEpe1 expression increases H3K9 methylation over several chromosomal regions, reducing expression of underlying genes and enhancing resistance. Reciprocally, constitutive non-cleavable Epe1 expression decreases resistance. tEpe1-mediated resistance requires a functional JmjC demethylase domain. Moreover, caffeine-induced Epe1-to-tEpe1 cleavage is dependent on an intact cell integrity MAP kinase stress signaling pathway, mutations in which alter resistance. Thus, environmental changes elicit a mechanism that curtails the function of this key epigenetic modifier, allowing heterochromatin to reprogram gene expression, thereby bestowing resistance to some cells within a population. H3K9me-heterochromatin components are conserved in human and crop-plant fungal pathogens for which a limited number of antifungals exist. Our findings reveal how transient heterochromatin-dependent antifungal resistant epimutations develop and thus inform on how they might be countered.

Non-Mendelian transmission of accessory chromosomes in fungi

Non-Mendelian transmission has been reported for various genetic elements, ranging from small transposons to entire chromosomes. One prime example of such a transmission pattern are B chromosomes in plants and animals. Accessory chromosomes in fungi are similar to B chromosomes in showing presence/absence polymorphism and being non-essential. How these chromosomes are transmitted during meiosis is however poorly understood—despite their often high impact on the fitness of the host. For several fungal organisms, a non-Mendelian transmission or a mechanistically unique meiotic drive of accessory chromosomes have been reported. In this review, we provide an overview of the possible mechanisms that can cause the non-Mendelian transmission or meiotic drives of fungal accessory chromosomes. We compare processes responsible for the non-Mendelian transmission of accessory chromosomes for different fungal eukaryotes and discuss the structural traits of fungal accessory chromosomes affecting their meiotic transmission. We conclude that research on fungal accessory chromosomes, due to their small size, ease of sequencing, and epigenetic profiling, can complement the study of B chromosomes in deciphering factors that influence and regulate the non-Mendelian transmission of entire chromosomes.

KDM6B interacts with TFDP1 to activate P53 signalling in regulating mouse palatogenesis

Epigenetic regulation plays extensive roles in diseases and development. Disruption of epigenetic regulation not only increases the risk of cancer, but can also cause various developmental defects. However, the question of how epigenetic changes lead to tissue-specific responses during neural crest fate determination and differentiation remains understudied. Using palatogenesis as a model, we reveal the functional significance of Kdm6b, an H3K27me3 demethylase, in regulating mouse embryonic development. Our study shows that Kdm6b plays an essential role in cranial neural crest development, and loss of Kdm6b disturbs P53 pathway-mediated activity, leading to complete cleft palate along with cell proliferation and differentiation defects in mice. Furthermore, activity of H3K27me3 on the promoter of Trp53 is antagonistically controlled by Kdm6b, and Ezh2 in cranial neural crest cells. More importantly, without Kdm6b, the transcription factor TFDP1, which normally binds to the promoter of Trp53, cannot activate Trp53 expression in palatal mesenchymal cells. Furthermore, the function of Kdm6b in activating Trp53 in these cells cannot be compensated for by the closely related histone demethylase Kdm6a. Collectively, our results highlight the important role of the epigenetic regulator KDM6B and how it specifically interacts with TFDP1 to achieve its functional specificity in regulating Trp53 expression, and further provide mechanistic insights into the epigenetic regulatory network during organogenesis.

Differential regulation and production of secondary metabolites among isolates of the fungal wheat pathogen Zymoseptoria tritici

The genome of the wheat pathogenic fungus, Zymoseptoria tritici, represents extensive presence-absence variation in gene content. Here, we addressed variation in biosynthetic gene clusters (BGCs) content and biochemical profiles among three isolates. We analysed secondary metabolite properties based on genome, transcriptome and metabolome data. The isolates represent highly distinct genome architecture, but harbor similar repertoire of BGCs. Expression profiles for most BGCs show comparable patterns of regulation among the isolates, suggesting a conserved "biochemical infection program". For all three isolates, we observed a strong up-regulation of a putative abscisic acid (ABA) gene cluster during biotrophic host colonization, indicating that Z. tritici potentially interfere with host defenses by the biosynthesis of this phytohormone. Further, during in vitro growth the isolates show similar metabolomes congruent with the predicted BGC content. We assessed if secondary metabolite production is regulated by histone methylation using a mutant impaired in formation of facultative heterochromatin (H3K27me3). In contrast to other ascomycete fungi, chromatin modifications play a less prominent role in regulation of secondary metabolites. In summary, we show that Z. tritici has a conserved program of secondary metabolite production contrasting the immense variation in effector expression, some of these metabolites might play a key role during host colonization. Importance Zymoseptoria tritici is one of the most devastating pathogens of wheat. So far the molecular determinants of virulence and their regulation are poorly understood. Previous studies have focused on proteinasous virulence factors and their extensive diversity. In this study, we focus on secondary metabolites produced by Z. tritici. Using a comparative framework, we here characterize core and non-core metabolites produced by Z. tritici by combining genome, transcriptome and metabolome datasets. Our findings indicate highly conserved biochemical profiles contrasting genetic and phenotypic diversity of the field isolates investigated here. This discovery has relevance for future crop protection strategies.

Polycomb Repressive Complex 2 in Eukaryotes-An Evolutionary Perspective

Polycomb repressive complex 2 (PRC2) represents a group of evolutionarily conserved multi-subunit complexes that repress gene transcription by introducing trimethylation of lysine 27 on histone 3 (H3K27me3). PRC2 activity is of key importance for cell identity specification and developmental phase transitions in animals and plants. The composition, biochemistry, and developmental function of PRC2 in animal and flowering plant model species are relatively well described. Recent evidence demonstrates the presence of PRC2 complexes in various eukaryotic supergroups, suggesting conservation of the complex and its function. Here, we provide an overview of the current understanding of PRC2-mediated repression in different representatives of eukaryotic supergroups with a focus on the green lineage. By comparison of PRC2 in different eukaryotes, we highlight the possible common and diverged features suggesting evolutionary implications and outline emerging questions and directions for future research of polycomb repression and its evolution.

Not all Is SET for Methylation: Evolution of Eukaryotic Protein Methyltransferases

Dynamic posttranslational modifications to canonical histones that constitute the nucleosome (H2A, H2B, H3, and H4) control all aspects of enzymatic transactions with DNA. Histone methylation has been studied heavily for the past 20 years, and our mechanistic understanding of the control and function of individual methylation events on specific histone arginine and lysine residues has been greatly improved over the past decade, driven by excellent new tools and methods. Here, we will summarize what is known about the distribution and some of the functions of protein methyltransferases from all major eukaryotic supergroups. The main conclusion is that protein, and specifically histone, methylation is an ancient process. Many taxa in all supergroups have lost some subfamilies of both protein arginine methyltransferases (PRMT) and the heavily studied SET domain lysine methyltransferases (KMT). Over time, novel subfamilies, especially of SET domain proteins, arose. We use the interactions between H3K27 and H3K36 methylation as one example for the complex circuitry of histone modifications that make up the "histone code," and we discuss one recent example (Paramecium Ezl1) for how extant enzymes that may resemble more ancient SET domain KMTs are able to modify two lysine residues that have divergent functions in plants, fungi, and animals. Complexity of SET domain KMT function in the well-studied plant and animal lineages arose not only by gene duplication but also acquisition of novel DNA- and histone-binding domains in certain subfamilies.

Chromatin profiling reveals heterogeneity in clinical isolates of the human pathogen Aspergillus fumigatus

Invasive Pulmonary Aspergillosis, which is caused by the filamentous fungus Aspergillus fumigatus, is a life-threatening infection for immunosuppressed patients. Chromatin structure regulation is important for genome stability maintenance and has the potential to drive genome rearrangements and affect virulence and pathogenesis of pathogens. Here, we performed the first A. fumigatus global chromatin profiling of two histone modifications, H3K4me3 and H3K9me3, focusing on the two most investigated A. fumigatus clinical isolates, Af293 and CEA17. In eukaryotes, H3K4me3 is associated with active transcription, while H3K9me3 often marks silent genes, DNA repeats, and transposons. We found that H3K4me3 deposition is similar between the two isolates, while H3K9me3 is more variable and does not always represent transcriptional silencing. Our work uncovered striking differences in the number, locations, and expression of transposable elements between Af293 and CEA17, and the differences are correlated with H3K9me3 modifications and higher genomic variations among strains of Af293 background. Moreover, we further showed that the Af293 strains from different laboratories actually differ in their genome contents and found a frequently lost region in chromosome VIII. For one such Af293 variant, we identified the chromosomal changes and demonstrated their impacts on its secondary metabolites production, growth and virulence. Overall, our findings not only emphasize the influence of genome heterogeneity on A. fumigatus fitness, but also caution about unnoticed chromosomal variations among common laboratory strains.

UvKmt6-mediated H3K27 trimethylation is required for development, pathogenicity, and stress response in Ustilaginoidea virens

Polycomb repressive complex 2 (PRC2) is responsible for the trimethylation of lysine 27 of histone H3 (H3K27me3)-mediated transcriptional silencing. At present, its biological roles in the devastating rice pathogenic fungus Ustilaginoidea virens remain unclear. In this study, we analyzed the function of a putative PRC2 catalytic subunit UvKmt6. The results showed that disruption of UvKMT6 resulted in reduced growth, conidiation and pathogenicity in U. virens. Furthermore, UvKmt6 is essential for establishment of H3K27me3 modification, which covers 321 genes in the genome. Deletion of UvKMT6 led to transcriptional derepression of 629 genes, 140 of which were occupied with H3K27me3 modification. Consistent with RNA-seq and ChIP-seq analysis, UvKmt6 was further confirmed to participate in the transcriptional repression of genes encoding effectors and genes associated with secondary metabolites production, such as PKSs, NRPSs and Cytochrome P450s. Notably, we found that UvKmt6 is involved in transcriptional repression of oxidative, osmotic, cell wall and nutrient starvation stresses response-related genes. From the perspective of gene expression and phenotype, in addition to the relatively conservative role in fungal development, virulence and production of secondary metabolites, we further reported that UvKmt6-mdediated H3K27me3 plays a critical role in the response to various stresses in U. virens.

Histone H3K27 methylation perturbs transcriptional robustness and underpins dispensability of highly conserved genes in fungi

Epigenetic modifications are key regulators of gene expression and underpin genome integrity. Yet, how epigenetic changes affect the evolution and transcriptional robustness of genes remains largely unknown. Here, we show how the repressive histone mark H3K27me3 underpins the trajectory of highly conserved genes in fungi. We first performed transcriptomic profiling on closely related species of the plant pathogen Fusarium graminearum species complex. We determined transcriptional responsiveness of genes across environmental conditions to determine expression robustness. To infer evolutionary conservation, we used a framework of 23 species across the Fusarium genus including three species covered with histone methylation data. Gene expression variation is negatively correlated with gene conservation confirming that highly conserved genes show higher expression robustness. In contrast, genes marked by H3K27me3 do not show such associations. Furthermore, highly conserved genes marked by H3K27me3 encode smaller proteins, exhibit weaker codon usage bias, higher levels of hydrophobicity, show lower intrinsically disordered regions, and are enriched for functions related to regulation and membrane transport. The evolutionary age of conserved genes with H3K27me3 histone marks falls typically within the origins of the Fusarium genus. We show that highly conserved genes marked by H3K27me3 are more likely to be dispensable for survival during host infection. Lastly, we show that conserved genes exposed to repressive H3K27me3 marks across distantly related Fusarium fungi are associated with transcriptional perturbation at the microevolutionary scale. In conclusion, we show how repressive histone marks are entangled in the evolutionary fate of highly conserved genes across evolutionary timescales.

Epigenetic modifications affect the rate of spontaneous mutations in a pathogenic fungus

Mutations are the source of genetic variation and the substrate for evolution. Genome-wide mutation rates appear to be affected by selection and are probably adaptive. Mutation rates are also known to vary along genomes, possibly in response to epigenetic modifications, but causality is only assumed. In this study we determine the direct impact of epigenetic modifications and temperature stress on mitotic mutation rates in a fungal pathogen using a mutation accumulation approach. Deletion mutants lacking epigenetic modifications confirm that histone mark H3K27me3 increases whereas H3K9me3 decreases the mutation rate. Furthermore, cytosine methylation in transposable elements (TE) increases the mutation rate 15-fold resulting in significantly less TE mobilization. Also accessory chromosomes have significantly higher mutation rates. Finally, we find that temperature stress substantially elevates the mutation rate. Taken together, we find that epigenetic modifications and environmental conditions modify the rate and the location of spontaneous mutations in the genome and alter its evolutionary trajectory.

Role of Polycomb in the control of transposable elements

It is generally considered that Polycomb Repressive Complex (PRC)2 deposits the histone mark H3K27me3 on silent protein-coding genes, while transposable elements are repressed by DNA and/or H3K9 methylation. Yet, there is increasing evidence that PRC2 also targets and even silences transposable elements in representatives of several distantly related eukaryotic lineages. In plants and animals, H3K27me3 is present on transposable elements in mutants and specific cell types devoid of DNA methylation. In this Opinion, we summarize the experimental evidence for this phenomenon across the eukaryotic kingdom, and discuss its functional and evolutionary significance. We hypothesize that an ancestral role of Polycomb group (PcG) proteins was to silence transposable elements.

Dissecting the Biology of the Fungal Wheat Pathogen Zymoseptoria tritici: A Laboratory Workflow

The fungus Zymoseptoria tritici is one of the most devastating pathogens of wheat. Aside from its importance as a disease-causing agent, this species has emerged as a powerful model system for evolutionary genetic studies of crop-infecting fungal pathogens. Z. tritici exhibits exceptionally high levels of genetic and phenotypic diversity as well as morphological plasticity, which can make experimental studies and comparability of results obtained in different laboratories, e.g., from infection assays, challenging. Therefore, standardized experimental methods are crucial for research on Z. tritici biology and the interaction of this fungus with its wheat host. Here, we describe a suite of well-tested and optimized protocols ranging from isolation of Z. tritici field specimens to analyses of virulence assays under controlled conditions. Several biological and technical aspects of working with Z. tritici under laboratory conditions are considered and carefully described in each protocol. © 2020 The Authors. Basic Protocol 1: Purification of Z. tritici field isolates from leaf material Basic Protocol 2: Molecular identification of Z. tritici isolates Support Protocol 1: Rapid extraction of Z. tritici genomic DNA Support Protocol 2: Extraction of high-quality Z. tritici genomic DNA Basic Protocol 3: In vitro culture and long-term storage of Z. tritici isolates Basic Protocol 4: Analysis of Z. tritici virulence in wheat Support Protocol 3: Preparation of Z. tritici inoculum.

Biosynthesis of Fusapyrone Depends on the H3K9 Methyltransferase, FmKmt1, in Fusarium mangiferae

The phytopathogenic fungus Fusarium mangiferae belongs to the Fusarium fujikuroi species complex (FFSC). Members of this group cause a wide spectrum of devastating diseases on diverse agricultural crops. F. mangiferae is the causal agent of the mango malformation disease (MMD) and as such detrimental for agriculture in the southern hemisphere. During plant infection, the fungus produces a plethora of bioactive secondary metabolites (SMs), which most often lead to severe adverse defects on plants health. Changes in chromatin structure achieved by posttranslational modifications (PTM) of histones play a key role in regulation of fungal SM biosynthesis. Posttranslational tri-methylation of histone 3 lysine 9 (H3K9me3) is considered a hallmark of heterochromatin and established by the SET-domain protein Kmt1. Here, we show that FmKmt1 is involved in H3K9me3 in F. mangiferae. Loss of FmKmt1 only slightly though significantly affected fungal hyphal growth and stress response and is required for wild type-like conidiation. While FmKmt1 is largely dispensable for the biosynthesis of most known SMs, removal of FmKMT1 resulted in an almost complete loss of fusapyrone and deoxyfusapyrone, γ-pyrones previously only known from Fusarium semitectum. Here, we identified the polyketide synthase (PKS) FmPKS40 to be involved in fusapyrone biosynthesis, delineate putative cluster borders by co-expression studies and provide insights into its regulation.

Transposable Elements Contribute to Genome Dynamics and Gene Expression Variation in the Fungal Plant Pathogen Verticillium dahliae

Transposable elements (TEs) are a major source of genetic and regulatory variation in their host genome and are consequently thought to play important roles in evolution. Many fungal and oomycete plant pathogens have evolved dynamic and TE-rich genomic regions containing genes that are implicated in host colonization and adaptation. TEs embedded in these regions have typically been thought to accelerate the evolution of these genomic compartments, but little is known about their dynamics in strains that harbor them. Here, we used whole-genome sequencing data of 42 strains of the fungal plant pathogen Verticillium dahliae to systematically identify polymorphic TEs that may be implicated in genomic as well as in gene expression variation. We identified 2,523 TE polymorphisms and characterize a subset of 8% of the TEs as polymorphic elements that are evolutionary younger, less methylated, and more highly expressed when compared with the remaining 92% of the total TE complement. As expected, the polyrmorphic TEs are enriched in the adaptive genomic regions. Besides, we observed an association of polymorphic TEs with pathogenicity-related genes that localize nearby and that display high expression levels. Collectively, our analyses demonstrate that TE dynamics in V. dahliae contributes to genomic variation, correlates with expression of pathogenicity-related genes, and potentially impacts the evolution of adaptive genomic regions.

Ago1 Affects the Virulence of the Fungal Plant Pathogen Zymoseptoria tritici

In host-pathogen interactions RNA interference (RNAi) has emerged as a pivotal mechanism to modify both, the immune responses of the host as well as the pathogenicity and virulence of the pathogen. In addition, in some fungi RNAi is also known to affect chromosome biology via its effect on chromatin conformation. Previous studies reported no effect of the RNAi machinery on the virulence of the fungal plant pathogen Zymoseptoria tritici however the role of RNAi is still poorly understood in this species. Herein, we elucidate whether the RNAi machinery is conserved within the genus Zymoseptoria. Moreover, we conduct functional analyses of Argonaute and Dicer-like proteins and test if the RNAi machinery affects chromosome stability. We show that the RNAi machinery is conserved among closely related Zymoseptoria species while an exceptional pattern of allelic diversity was possibly caused by introgression. The deletion of Ago1 reduced the ability of the fungus to produce asexual propagules in planta in a quantitative matter. Chromosome stability of the accessory chromosome of Z. tritici was not prominently affected by the RNAi machinery. These results indicate, in contrast to previous finding, a role of the RNAi pathway during host infection, but not in the stability of accessory chromosomes in Z. tritici.

Machine-learning predicts genomic determinants of meiosis-driven structural variation in a eukaryotic pathogen

Species harbor extensive structural variation underpinning recent adaptive evolution. However, the causality between genomic features and the induction of new rearrangements is poorly established. Here, we analyze a global set of telomere-to-telomere genome assemblies of a fungal pathogen of wheat to establish a nucleotide-level map of structural variation. We show that the recent emergence of pesticide resistance has been disproportionally driven by rearrangements. We use machine learning to train a model on structural variation events based on 30 chromosomal sequence features. We show that base composition and gene density are the major determinants of structural variation. Retrotransposons explain most inversion, indel and duplication events. We apply our model to Arabidopsis thaliana and show that our approach extends to more complex genomes. Finally, we analyze complete genomes of haploid offspring in a four-generation pedigree. Meiotic crossover locations are enriched for new rearrangements consistent with crossovers being mutational hotspots. The model trained on species-wide structural variation accurately predicts the position of >74% of newly generated variants along the pedigree. The predictive power highlights causality between specific sequence features and the induction of chromosomal rearrangements. Our work demonstrates that training sequence-derived models can accurately identify regions of intrinsic DNA instability in eukaryotic genomes.

Genome-wide mapping of histone modifications during axenic growth in two species of Leptosphaeria maculans showing contrasting genomic organization

Leptosphaeria maculans 'brassicae' (Lmb) and Leptosphaeria maculans 'lepidii' (Lml) are closely related phytopathogenic species that exhibit a large macrosynteny but contrasting genome structure. Lmb has more than 30% of repeats clustered in large repeat-rich regions, while the Lml genome has only a small amount of evenly distributed repeats. Repeat-rich regions of Lmb are enriched in effector genes, expressed during plant infection. The distinct genome structures of Lmb and Lml provide an excellent model for comparing the organization of pathogenicity genes in relation to the chromatin landscape in two closely related phytopathogenic fungi. Here, we performed chromatin immunoprecipitation (ChIP) during axenic culture, targeting histone modifications typical for heterochromatin or euchromatin, combined with transcriptomic analysis to analyze the influence of chromatin organization on gene expression. In both species, we found that facultative heterochromatin is enriched with genes lacking functional annotation, including numerous effector and species-specific genes. Notably, orthologous genes located in H3K27me3 domains are enriched with effector genes. Compared to other fungal species, including Lml, Lmb is distinct in having large H3K9me3 domains associated with repeat-rich regions that contain numerous species-specific effector genes. Discovery of these two distinctive heterochromatin landscapes now raises questions about their involvement in the regulation of pathogenicity, the dynamics of these domains during plant infection and the selective advantage to the fungus to host effector genes in H3K9me3 or H3K27me3 domains.

Loss of EZH2-like or SU(VAR)3-9-like proteins causes simultaneous perturbations in H3K27 and H3K9 tri-methylation and associated developmental defects in the fungus Podospora anserina

BACKGROUND

Selective gene silencing is key to development. It is generally accepted that H3K27me3-enriched heterochromatin maintains transcriptional repression established during early development and regulates cell fate. Conversely, H3K9me3-enriched heterochromatin prevents differentiation but constitutes protection against transposable elements. We exploited the fungus Podospora anserina, a valuable alternative to higher eukaryote models, to question the biological relevance and functional interplay of these two distinct heterochromatin conformations.

RESULTS

We established genome-wide patterns of H3K27me3 and H3K9me3 modifications, and found these marks mutually exclusive within gene-rich regions but not within repeats. We generated the corresponding histone methyltransferase null mutants and showed an interdependence of H3K9me3 and H3K27me3 marks. Indeed, removal of the PaKmt6 EZH2-like enzyme resulted not only in loss of H3K27me3 but also in significant H3K9me3 reduction. Similarly, removal of PaKmt1 SU(VAR)3-9-like enzyme caused loss of H3K9me3 and substantial decrease of H3K27me3. Removal of the H3K9me binding protein PaHP1 provided further support to the notion that each type of heterochromatin requires the presence of the other. We also established that P. anserina developmental programs require H3K27me3-mediated silencing, since loss of the PaKmt6 EZH2-like enzyme caused severe defects in most aspects of the life cycle including growth, differentiation processes and sexual reproduction, whereas loss of the PaKmt1 SU(VAR)3-9-like enzyme resulted only in marginal defects, similar to loss of PaHP1.

CONCLUSIONS

Our findings support a conserved function of the PRC2 complex in fungal development. However, we uncovered an intriguing evolutionary fluidity in the repressive histone deposition machinery, which challenges canonical definitions of constitutive and facultative heterochromatin.

Three putative DNA methyltransferases of Verticillium dahliae differentially contribute to DNA methylation that is dispensable for growth, development and virulence

BACKGROUND

DNA methylation is an important epigenetic control mechanism that in many fungi is restricted to genomic regions containing transposable elements (TEs). Two DNA methyltransferases, Dim2 and Dnmt5, are known to perform methylation at cytosines in fungi. While most ascomycete fungi encode both Dim2 and Dnmt5, only few functional studies have been performed in species containing both.

METHODS

In this study, we report functional analysis of both Dim2 and Dnmt5 in the plant pathogenic fungus Verticillium dahliae.

RESULTS

Our results show that Dim2, but not Dnmt5 or the putative sexual-cycle-related DNA methyltransferase Rid, is responsible for the majority of DNA methylation under the tested conditions. Single or double DNA methyltransferase mutants did not show altered development, virulence, or transcription of genes or TEs. In contrast, Hp1 and Dim5 mutants that are impacted in chromatin-associated processes upstream of DNA methylation are severely affected in development and virulence and display transcriptional reprogramming in specific hypervariable genomic regions (so-called adaptive genomic regions) that contain genes associated with host colonization. As these adaptive genomic regions are largely devoid of DNA methylation and of Hp1- and Dim5-associated heterochromatin, the differential transcription is likely caused by pleiotropic effects rather than by differential DNA methylation.

CONCLUSION

Overall, our study suggests that Dim2 is the main DNA methyltransferase in V. dahliae and, in conjunction with work on other fungi, is likely the main active DNMT in ascomycetes, irrespective of Dnmt5 presence. We speculate that Dnmt5 and Rid act under specific, presently enigmatic, conditions or, alternatively, act in DNA-associated processes other than DNA methylation.

Recent loss of the Dim2 DNA methyltransferase decreases mutation rate in repeats and changes evolutionary trajectory in a fungal pathogen

DNA methylation is found throughout all domains of life, yet the extent and function of DNA methylation differ among eukaryotes. Strains of the plant pathogenic fungus Zymoseptoria tritici appeared to lack cytosine DNA methylation (5mC) because gene amplification followed by Repeat-Induced Point mutation (RIP) resulted in the inactivation of the dim2 DNA methyltransferase gene. 5mC is, however, present in closely related sister species. We demonstrate that inactivation of dim2 occurred recently as some Z. tritici isolates carry a functional dim2 gene. Moreover, we show that dim2 inactivation occurred by a different path than previously hypothesized. We mapped the genome-wide distribution of 5mC in strains with or without functional dim2 alleles. Presence of functional dim2 correlates with high levels of 5mC in transposable elements (TEs), suggesting a role in genome defense. We identified low levels of 5mC in strains carrying non-functional dim2 alleles, suggesting that 5mC is maintained over time, presumably by an active Dnmt5 DNA methyltransferase. Integration of a functional dim2 allele in strains with mutated dim2 restored normal 5mC levels, demonstrating de novo cytosine methylation activity of Dim2. To assess the importance of 5mC for genome evolution, we performed an evolution experiment, comparing genomes of strains with high levels of 5mC to genomes of strains lacking functional dim2. We found that presence of a functional dim2 allele alters nucleotide composition by promoting C to T transitions (C→T) specifically at CpA (CA) sites during mitosis, likely contributing to TE inactivation. Our results show that 5mC density at TEs is a polymorphic trait in Z. tritici populations that can impact genome evolution.

Local Rather than Global H3K27me3 Dynamics Are Associated with Differential Gene Expression in Verticillium dahliae

Differential growth conditions typically trigger global transcriptional responses in filamentous fungi. Such fungal responses to environmental cues involve epigenetic regulation, including chemical histone modifications. It has been proposed that conditionally expressed genes, such as those that encode secondary metabolites but also effectors in pathogenic species, are often associated with a specific histone modification, lysine27 methylation of H3 (H3K27me3). However, thus far, no analyses on the global H3K27me3 profiles have been reported under differential growth conditions in order to assess if H3K27me3 dynamics govern differential transcription. Using chromatin immunoprecipitation sequencing (ChIP-seq) and RNA sequencing data from the plant-pathogenic fungus Verticillium dahliae grown in three in vitro cultivation media, we now show that a substantial number of the identified H3K27me3 domains globally display stable profiles among these growth conditions. However, we observe local quantitative differences in H3K27me3 ChIP-seq signals that are associated with a subset of differentially transcribed genes between media. Comparing the in vitro results to expression during plant infection suggests that in planta-induced genes may require chromatin remodeling to achieve expression. Overall, our results demonstrate that some loci display H3K27me3 dynamics associated with concomitant transcriptional variation, but many differentially expressed genes are associated with stable H3K27me3 domains. Thus, we conclude that while H3K27me3 is required for transcriptional repression, it does not appear that transcriptional activation requires the global erasure of H3K27me3. We propose that the H3K27me3 domains that do not undergo dynamic methylation may contribute to transcription through other mechanisms or may serve additional genomic regulatory functions. IMPORTANCE In many organisms, including filamentous fungi, epigenetic mechanisms that involve chemical and physical modifications of DNA without changing the genetic sequence have been implicated in transcriptional responses upon developmental or environmental cues. In fungi, facultative heterochromatin that can decondense to allow transcription in response to developmental changes or environmental stimuli is characterized by the trimethylation of lysine 27 on histone H3 (H3K27me3), and H3K27me3 has been implicated in transcriptional regulation, although the precise mechanisms and functions remain enigmatic. Based on ChIP and RNA sequencing data, we show for the soilborne broad-host-range vascular wilt plant-pathogenic fungus Verticillium dahliae that although some loci display H3K27me3 dynamics that can contribute to transcriptional variation, other loci do not show such a dependence. Thus, although we recognize that H3K27me3 is required for transcriptional repression, we also conclude that this mark is not a conditionally responsive global regulator of differential transcription upon responses to environmental cues.

A unique chromatin profile defines adaptive genomic regions in a fungal plant pathogen

Genomes store information at scales beyond the linear nucleotide sequence, which impacts genome function at the level of an individual, while influences on populations and long-term genome function remains unclear. Here, we addressed how physical and chemical DNA characteristics influence genome evolution in the plant pathogenic fungus Verticillium dahliae. We identified incomplete DNA methylation of repetitive elements, associated with specific genomic compartments originally defined as Lineage-Specific (LS) regions that contain genes involved in host adaptation. Further chromatin characterization revealed associations with features such as H3 Lys-27 methylated histones (H3K27me3) and accessible DNA. Machine learning trained on chromatin data identified twice as much LS DNA as previously recognized, which was validated through orthogonal analysis, and we propose to refer to this DNA as adaptive genomic regions. Our results provide evidence that specific chromatin profiles define adaptive genomic regions, and highlight how different epigenetic factors contribute to the organization of these regions.

Partial pathogenicity chromosomes in Fusarium oxysporum are sufficient to cause disease and can be horizontally transferred

In Fusarium oxysporum f.sp. lycopersici, all effector genes reported so far - also called SIX genes - are located on a single accessory chromosome which is required for pathogenicity and can also be horizontally transferred to another strain. To narrow down the minimal region required for virulence, we selected partial pathogenicity chromosome deletion strains by fluorescence-assisted cell sorting of a strain in which the two arms of the pathogenicity chromosome were labelled with GFP and RFP respectively. By testing the virulence of these deletion mutants, we show that the complete long arm and part of the short arm of the pathogenicity chromosome are not required for virulence. In addition, we demonstrate that smaller versions of the pathogenicity chromosome can also be transferred to a non-pathogenic strain and they are sufficient to turn the non-pathogen into a pathogen. Surprisingly, originally non-pathogenic strains that had received a smaller version of the pathogenicity chromosome were much more aggressive than recipients with a complete pathogenicity chromosome. Whole genome sequencing analysis revealed that partial deletions of the pathogenicity chromosome occurred mainly close to repeats, and that spontaneous duplication of sequences in accessory regions is frequent both in chromosome deletion strains and in horizontal transfer strains.

Chromatin Dynamics Contribute to the Spatiotemporal Expression Pattern of Virulence Genes in a Fungal Plant Pathogen

Dynamic changes in transcription profiles are key for the success of pathogens in colonizing their hosts. In many pathogens, genes associated with virulence, such as effector genes, are located in regions of the genome that are rich in transposable elements and heterochromatin. The contribution of chromatin modifications to gene expression in pathogens remains largely unknown. Using a combination of a reporter gene-based approach and chromatin immunoprecipitation, we show that the heterochromatic environment of effector genes in the fungal plant pathogen Zymoseptoria tritici is a key regulator of their specific spatiotemporal expression patterns. Enrichment in trimethylated lysine 27 of histone H3 dictates the repression of effector genes in the absence of the host. Chromatin decondensation during host colonization, featuring a reduction in this repressive modification, indicates a major role for epigenetics in effector gene induction. Our results illustrate that chromatin modifications triggered during host colonization determine the specific expression profile of effector genes at the cellular level and, hence, provide new insights into the regulation of virulence in fungal plant pathogens.IMPORTANCE Fungal plant pathogens possess a large repertoire of genes encoding putative effectors, which are crucial for infection. Many of these genes are expressed at low levels in the absence of the host but are strongly induced at specific stages of the infection. The mechanisms underlying this transcriptional reprogramming remain largely unknown. We investigated the role of the genomic environment and associated chromatin modifications of effector genes in controlling their expression pattern in the fungal wheat pathogen Zymoseptoria tritici Depending on their genomic location, effector genes are epigenetically repressed in the absence of the host and during the initial stages of infection. Derepression of effector genes occurs mainly during and after penetration of plant leaves and is associated with changes in histone modifications. Our work demonstrates the role of chromatin in shaping the expression of virulence components and, thereby, the interaction between fungal pathogens and their plant hosts.

Accessory Chromosomes in Fusarium oxysporum

Most genomes within the species complex of Fusarium oxysporum are organized into two compartments: the core chromosomes (CCs) and accessory chromosomes (ACs). As opposed to CCs, which are conserved and vertically transmitted to carry out essential housekeeping functions, lineage- or strain-specific ACs are believed to be initially horizontally acquired through unclear mechanisms. These two genomic compartments are different in terms of gene density, the distribution of transposable elements, and epigenetic markers. Although common in eukaryotes, the functional importance of ACs is uniquely emphasized among fungal species, specifically in relationship to fungal pathogenicity and their adaptation to diverse hosts. With a focus on the cross-kingdom fungal pathogen F. oxysporum, this review provides a summary of the differences between CCs and ACs based on current knowledge of gene functions, genome structures, and epigenetic signatures, and explores the transcriptional crosstalk between the core and accessory genomes.

Genome compartmentalization predates species divergence in the plant pathogen genus Zymoseptoria

Background

Antagonistic co-evolution can drive rapid adaptation in pathogens and shape genome architecture. Comparative genome analyses of several fungal pathogens revealed highly variable genomes, for many species characterized by specific repeat-rich genome compartments with exceptionally high sequence variability. Dynamic genome structure may enable fast adaptation to host genetics. The wheat pathogen Zymoseptoria tritici with its highly variable genome, has emerged as a model organism to study genome evolution of plant pathogens. Here, we compared genomes of Z. tritici isolates and of sister species infecting wild grasses to address the evolution of genome composition and structure.

Results

Using long-read technology, we sequenced and assembled genomes of Z. ardabiliae, Z. brevis, Z. pseudotritici and Z. passerinii, together with two isolates of Z. tritici. We report a high extent of genome collinearity among Zymoseptoria species and high conservation of genomic, transcriptomic and epigenomic signatures of compartmentalization. We identify high gene content variability both within and between species. In addition, such variability is mainly limited to the accessory chromosomes and accessory compartments. Despite strong host specificity and non-overlapping host-range between species, predicted effectors are mainly shared among Zymoseptoria species, yet exhibiting a high level of presence-absence polymorphism within Z. tritici. Using in planta transcriptomic data from Z. tritici, we suggest different roles for the shared orthologs and for the accessory genes during infection of their hosts.

Conclusion

Despite previous reports of high genomic plasticity in Z. tritici, we describe here a high level of conservation in genomic, epigenomic and transcriptomic composition and structure across the genus Zymoseptoria. The compartmentalized genome allows the maintenance of a functional core genome co-occurring with a highly variable accessory genome.

Advances in Therapeutic Targeting of Cancer Stem Cells within the Tumor Microenvironment: An Updated Review

Despite great strides being achieved in improving cancer patients' outcomes through better therapies and combinatorial treatment, several hurdles still remain due to therapy resistance, cancer recurrence and metastasis. Drug resistance culminating in relapse continues to be associated with fatal disease. The cancer stem cell theory posits that tumors are driven by specialized cancer cells called cancer stem cells (CSCs). CSCs are a subpopulation of cancer cells known to be resistant to therapy and cause metastasis. Whilst the debate on whether CSCs are the origins of the primary tumor rages on, CSCs have been further characterized in many cancers with data illustrating that CSCs display great abilities to self-renew, resist therapies due to enhanced epithelial to mesenchymal (EMT) properties, enhanced expression of ATP-binding cassette (ABC) membrane transporters, activation of several survival signaling pathways and increased immune evasion as well as DNA repair mechanisms. CSCs also display great heterogeneity with the consequential lack of specific CSC markers presenting a great challenge to their targeting. In this updated review we revisit CSCs within the tumor microenvironment (TME) and present novel treatment strategies targeting CSCs. These promising strategies include targeting CSCs-specific properties using small molecule inhibitors, immunotherapy, microRNA mediated inhibitors, epigenetic methods as well as targeting CSC niche-microenvironmental factors and differentiation. Lastly, we present recent clinical trials undertaken to try to turn the tide against cancer by targeting CSC-associated drug resistance and metastasis.

A chromosome-level genome assembly provides new insights into paternal genome elimination in the cotton mealybug Phenacoccus solenopsis

Mealybugs (Hemiptera: Pseudococcidae) are economically important agricultural pests with several compelling biological phenomena including paternal genome elimination (PGE). However, limited high-quality genome assemblies of mealybugs hinder a full understanding of this striking and unusual biological phenomenon. Here, we generated a chromosome-level genome assembly of cotton mealybug, Phenacoccus solenopsis, by combining Illumina short reads, PacBio long reads and Hi-C scaffolding. The assembled genome was 292.54 Mb with a contig N50 of 489.8 kb and a scaffold N50 of 49.0 Mb. Hi-C scaffolding assigned 84.42% of the bases to five chromosomes. A total of 110.75 Mb (37.9%) repeat sequences and 11,880 protein-coding genes were predicted. The completeness of the genome assembly was estimated to be 95.5% based on BUSCO genes. In addition, 27,086 (95.3%) full-length PacBio transcripts were uniquely mapped to the assembled scaffolds, suggesting the high quality of the genome assembly. We showed that cotton mealybugs lack differentiated sex chromosomes by analysing genome resequencing data of males and females. DAPI staining confirmed that one chromosome set in males becomes heterochromatin at an early embryo stage. Chromatin immunoprecipitation assays with sequencing analysis demonstrated that the epigenetic modifications H3K9me3 and H3K27me3 are distributed across the whole genome in males, suggesting that these two modifications might be involved in maintaining heterochromatin status. Both markers were more likely to be distributed in repeat regions, while H3K27me3 had higher overall enrichment. Our results provide a valuable genomic resource and shed new light on the genomic and epigenetic basis of PGE in cotton mealybugs.

Imprinting aberrations of SNRPN, ZAC1 and INPP5F genes involved in the pathogenesis of congenital heart disease with extracardiac malformations

Congenital heart disease (CHD) with extracardiac malformations (EM) is the most common multiple malformation, resulting from the interaction between genetic abnormalities and environmental factors. Most studies have attributed the causes of CHD with EM to chromosomal abnormalities. However, multi‐system dysplasia is usually caused by both genetic mutations and epigenetic dysregulation. The epigenetic mechanisms underlying the pathogenesis of CHD with EM remain unclear. In this study, we investigated the mechanisms of imprinting alterations, including those of the Small nuclear ribonucleoprotein polypeptide N (SNRPN), PLAG1 like zinc finger 1 (ZAC1) and inositol polyphosphate‐5‐phosphatase F (INPP5F) genes, in the pathogenesis of CHD with EM. The methylation levels of SNRPN, ZAC1, and INPP5F genes were analysed by the MassARRAY platform in 24 children with CHD with EM and 20 healthy controls. The expression levels of these genes were detected by real‐time polymerase chain reaction (PCR). The correlation between methylation regulation and gene expression was confirmed using 5‐azacytidine (5‐Aza) treated cells. The methylation levels of SNRPN and ZAC1 genes were significantly increased in CHD with EM, while that of INPP5F was decreased. The methylation alterations of these genes were negatively correlated with expression. Risk analysis showed that abnormal hypermethylation of SNRPN and ZAC1 resulted in 5.545 and 7.438 times higher risks of CHD with EM, respectively, and the abnormal hypomethylation of INPP5F was 8.38 times higher than that of the control group. We concluded that abnormally high methylation levels of SNRPN and ZAC1 and decreased levels of INPP5F imply an increased risk of CHD with EM by altering their gene functions. This study provides evidence of imprinted regulation in the pathogenesis of multiple malformations.

Polycomb Repression without Bristles: Facultative Heterochromatin and Genome Stability in Fungi

Genome integrity is essential to maintain cellular function and viability. Consequently, genome instability is frequently associated with dysfunction in cells and associated with plant, animal, and human diseases. One consequence of relaxed genome maintenance that may be less appreciated is an increased potential for rapid adaptation to changing environments in all organisms. Here, we discuss evidence for the control and function of facultative heterochromatin, which is delineated by methylation of histone H3 lysine 27 (H3K27me) in many fungi. Aside from its relatively well understood role in transcriptional repression, accumulating evidence suggests that H3K27 methylation has an important role in controlling the balance between maintenance and generation of novelty in fungal genomes. We present a working model for a minimal repressive network mediated by H3K27 methylation in fungi and outline challenges for future research.