Natural genomic variation in rice blast genomes is associated with specific heterochromatinmodifications

Genome organization in eukaryotes exhibits non-random patterns tied to transcription, replication, and chromatin. However, the driving forces across these processes, and their impacts on genome evolution remain unclear. To address this, we analyzed sequence data from 86 Magnaporthe oryzae isolates, a globally important plant pathogenic fungus, alongside new reference genomes, to investigate DNA sequence variation and the epigenome. Histone modification profiles were obtained through genome-wide chromatin immunoprecipitation-sequencing of the four new reference strains, which revealed that repressive histone marks (H3K27me3, H3K9me3) were associated with SNP and INDEL frequency. Densely grouped SNPs were found to reside in heterochromatin, often outside transposable elements, highlighting the link between heterochromatin and DNA variation. Even when controlling for selection, silent SNP frequency/kb was higher in H3K27me3-associated genes. Effector genes, key to pathogenicity, also displayed this trend. Comparing the reference strains, euchromatic regions were often syntenic, while heterochromatic regions trended towards non-syntenic. Heterochromatin emerged as a major factor associated with diverse DNA variations in M. oryzae populations, even when accounting for selective pressure. This underscores heterochromatins pivotal role in shaping genetic diversity in these mainly asexually reproducing fungi.