Many filamentous plant pathogens exhibit high levels of genomic variability, yet the impact of this variation on host–pathogen interactions is largely unknown. We have addressed host specialization in the wheat pathogen Zymoseptoria tritici. Our study builds on comparative analyses of infection and gene expression phenotypes of three isolates and reveals the extent to which genomic variation translates into phenotypic variation. The isolates exhibit genetic and genomic variation but are similarly virulent. By combining confocal microscopy, disease monitoring, staining of ROS, and comparative transcriptome analyses, we conducted a detailed comparison of the infection processes of these isolates in a susceptible wheat cultivar. We characterized four core infection stages: establishment, biotrophic growth, lifestyle transition, and necrotrophic growth and asexual reproduction that are shared by the three isolates. However, we demonstrate differentiated temporal and spatial infection development and significant differences in the expression profiles of the three isolates during the infection stages. More than 20% of the genes were differentially expressed and these genes were located significantly closer to transposable elements, suggesting an impact of epigenetic regulation. Further, differentially expressed genes were enriched in effector candidates suggesting that isolate‐specific strategies for manipulating host defenses are present in Z. tritici. We demonstrate that individuals of a host‐specialized pathogen have highly differentiated infection programs characterized by flexible infection development and functional redundancy. This illustrates how high genetic diversity in pathogen populations results in highly differentiated infection phenotypes, which fact needs to be acknowledged to understand host–pathogen interactions and pathogen evolution.

Yield losses caused by fungal pathogens represent a major threat to global food production. One of the most devastating fungal wheat pathogens is Zymoseptoria tritici. Despite the importance of this fungus, the underlying mechanisms of plant-pathogen interactions are poorly understood. Here we present a conceptual framework based on coinfection assays, comparative metabolomics, and microbiome profiling to study the interaction of Z. tritici in susceptible and resistant wheat. We demonstrate that Z. tritici suppresses the production of immune-related metabolites in a susceptible cultivar. Remarkably, this fungus-induced immune suppression spreads within the leaf and even to other leaves, a phenomenon that we term "systemic induced susceptibility". Using a comparative metabolomics approach, we identify defense-related biosynthetic pathways that are suppressed and induced in susceptible and resistant cultivars, respectively. We show that these fungus-induced changes correlate with changes in the wheat leaf microbiome. Our findings suggest that immune suppression by this hemibiotrophic pathogen impacts specialized plant metabolism, alters its associated microbial communities, and renders wheat vulnerable to further infections.

Chromosome and genome stability are important for normal cell function as instability often correlates with disease and dysfunction of DNA repair mechanisms. Many organisms maintain supernumerary or accessory chromosomes that deviate from standard chromosomes. The pathogenic fungus Zymoseptoria tritici has as many as eight accessory chromosomes, which are highly unstable during meiosis and mitosis, transcriptionally repressed, show enrichment of repetitive elements, and enrichment with heterochromatic histone methylation marks, e.g., trimethylation of H3 lysine 9 or lysine 27 (H3K9me3, H3K27me3). To elucidate the role of heterochromatin on genome stability in Z . tritici , we deleted the genes encoding the methyltransferases responsible for H3K9me3 and H3K27me3, kmt1 and kmt6 , respectively, and generated a double mutant. We combined experimental evolution and genomic analyses to determine the impact of these deletions on chromosome and genome stability, both in vitro and in planta . We used whole genome sequencing, ChIP-seq, and RNA-seq to compare changes in genome and chromatin structure, and differences in gene expression between mutant and wildtype strains. Analyses of genome and ChIP-seq data in H3K9me3-deficient strains revealed dramatic chromatin reorganization, where H3K27me3 is mostly relocalized into regions that are enriched with H3K9me3 in wild type. Many genome rearrangements and formation of new chromosomes were found in the absence of H3K9me3, accompanied by activation of transposable elements. In stark contrast, loss of H3K27me3 actually increased the stability of accessory chromosomes under normal growth conditions in vitro , even without large scale changes in gene activity. We conclude that H3K9me3 is important for the maintenance of genome stability because it disallows H3K27me3 in regions considered constitutive heterochromatin. In this system, H3K27me3 reduces the overall stability of accessory chromosomes, generating a “metastable” state for these quasi-essential regions of the genome.