Whole-genome landscape of H3K4me3, H3K36me3 and H3K9ac and their association with gene expression during Paulownia witches' broom disease infection and recovery processes

Histone methylation and acetylation participate in the modulation of gene expression. Here, chromatin immunoprecipitation sequencing (ChIP-Seq) was used to determine genome-wide patterns of three histone modifications, H3K4me3, H3K36me3, and H3K9ac (associated with actively expressed genes) and their associations with gene expression inPaulownia fortuneifollowing phytoplasma infection and recovery from Paulownia witches' broom (PaWB) disease after methyl methane sulfonate treatment. The three histone marks were preferentially deposited in genic regions, especially downstream of transcription start sites, and were highly concurrent with gene expression. Genes with all three histone marks exhibited the highest expression levels. Based on the comparison scheme, we detected 365, 2244, and 752 PaWB-associated genes with H3K4me3, H3K36me3, and H3K9ac marks, separately. KEGG pathway analysis showed that these genes were involved in plant-pathogen interaction, plant hormone signal transduction, and starch and sucrose metabolism. A small proportion of differentially modified genes showed changes in expression in response to phytoplasma infection, including genes involved in calcium ion signal transduction, abscisic acid signal transduction, and ethylene biosynthesis. This comprehensive analysis of genome-wide histone modifications and gene expression in Paulownia following phytoplasma infection provides new insights into the epigenetic responses to phytoplasma infection and will be useful for further studies on epigenetic regulation mechanisms in plants under biotic stress.