- Open Access
MAPK-triggered chromatin reprogramming by histone deacetylase in plant innate immunity
Genome Biologyvolume 18, Article number: 131 (2017)
Microbial-associated molecular patterns activate several MAP kinases, which are major regulators of the innate immune response in Arabidopsis thaliana that induce large-scale changes in gene expression. Here, we determine whether microbial-associated molecular pattern-triggered gene expression involves modifications at the chromatin level.
Histone acetylation and deacetylation are major regulators of microbial-associated molecular pattern-triggered gene expression and implicate the histone deacetylase HD2B in the reprogramming of defence gene expression and innate immunity. The MAP kinase MPK3 directly interacts with and phosphorylates HD2B, thereby regulating the intra-nuclear compartmentalization and function of the histone deacetylase.
By studying a number of gene loci that undergo microbial-associated molecular pattern-dependent activation or repression, our data reveal a mechanistic model for how protein kinase signaling directly impacts chromatin reprogramming in plant defense.
Due to their sessile nature, plants have developed sophisticated ways to respond and adapt to a variety of external stress factors that would otherwise compromise proper development, reproductive success, and ultimately survival. Numerous cellular proteins interact and communicate in response to extracellular stimuli and, through multiple signaling networks, transmit signals to the nucleus for reprogramming chromosomal gene expression. These dynamic regulatory mechanisms contribute to the capacity of plants to adapt to the onslaught of both biotic and abiotic challenges. Indeed, the success of photosynthetic eukaryotes is influenced by the adaptive dynamics of chromatin regulatory mechanisms, like histone modifications, which are rapid and reversible.
The choice of gene expression ultimately determines the fate of cells, forming the basis of biological diversity. The regulation of gene expression is closely coupled to chromatin structure and its modifications, which determine the accessibility of many regulatory proteins and non-coding RNAs (ncRNAs) to the DNA, adding a further layer of complexity to the genetic information encoded by the DNA sequence [1,2,3]. Chromatin is a tightly contained higher order structure that compacts genomic DNA to fit within the nucleus. The fundamental unit of chromatin is the nucleosome, which is composed of DNA that is wrapped around an octamer of histone proteins. Chromatin structure is modulated by a variety of mechanisms including DNA methylation catalyzed by DNA cytosine methyltransferases, histone post-translational modifications, such as acetylation and methylation, catalyzed by a wide range of enzymes specific for each modification, alterations in histone-DNA interactions that facilitate nucleosome sliding and are catalyzed by chromatin remodeling complexes, histone variants, and long and small ncRNAs that can act directly on chromatin and induce RNA-dependent DNA methylation (RdDM) [4,5,6].
Many developmental and environmental cues induce changes in chromatin structure. Plants sense pathogens through the perception of MAMPs, which induce signaling cascades to activate transcription factors and invoke chromatin regulatory mechanisms to reorganize the chromatin structure and, ultimately, provoke the changes in gene expression necessary for plant defense. Thus, in response to different stimuli, a single eukaryotic genome (DNA sequence) can give rise to distinct epigenomes.
Here, we examine the role of chromatin remodeling in Arabidopsis thaliana upon challenge with a synthetically produced 22 amino-acid long flagellin peptide (flg22) that mimics the response to bacterial pathogens. Flg22 is recognized in Arabidopsis by the plasma membrane leucine-rich repeat-receptor kinase (LRR-RK) FLS2 and activates two MAPK signaling pathways that initiate an array of defense responses including the production of several hormones, reactive oxygen species, and the induction of a large set of defense genes, processes generally referred to as MAMP-triggered immunity (MTI). These two main cascade branches involve the successive recruitment of MAP kinase kinase kinases (MAP3Ks), which phosphorylate and activate MAP kinase kinases (MAP2Ks) that phosphorylate and activate the MAP kinases MPK3, 4, 6, and 11 [7,8,9]. Ultimately, these MAPKs phosphorylate and thereby regulate protein factors responsible for the increased or decreased expression of specific gene sets with the goal to counteract the pathogen assault. Although specific transcription factors have been identified as targets for MPK3 , MPK4 , and MPK6 , much of the protein machinery orchestrating gene regulation in response to flagellin has not been identified.
Here, we report a role for the MAP kinase MPK3 in chromatin modulation and dissect the defense mechanism during the response to flagellin. We show that MPK3 interacts directly with and phosphorylates the histone deacetylase HD2B, which has been shown to deacetylate the lysine 9 residue of histone 3 (H3K9), a modification generally linked to the compaction of chromatin. We also show that, upon flagellin perception, HD2B is relocalized from the nucleolus to the nuclear compartment, leading to global genome-wide shifts in the H3K9 acetylation landscape. We identified flg22-regulated defense genes targeted by the MPK3-HD2B regulatory module and show that HD2B is directly implicated in bacterial defense in plants. Our results mechanistically define how a MAMP-activated MAP kinase regulates global changes in the chromatin landscape.
MPK3 interacts with and phosphorylates HD2B in vivo
Our previous phosphoproteomic approach aimed at identifying MAPK substrates revealed 303 in vivo phosphorylation sites in proteins isolated from Arabidopsis root cells . Among these sites, 91 matched the proline-directed motifs pS/pT-P that commonly serve as phosphorylation sites for MAPKs. Because we were interested to identify global regulators of gene expression that play a role in pathogen defense, we selected the histone deacetylase (HDAC) HD2B, which was phosphorylated at the amino acid positions T249 and S266, matching the S/T-P consensus motif for MAPK substrates (Additional file 1: Figures S1–S5).
We first assayed whether HD2B was a substrate of the canonical MAPK pathways activated by pathogen recognition. To this end, we immuno-purified endogenous flg22-activated MPK3, MPK4, and MPK6 from a root cell culture extract and tested whether the MAPKs phosphorylated HD2B in vitro. Using MBP and GST as positive and negative controls, respectively, MPK3 showed the highest phosphorylation activity toward GST-tagged HD2B among the three tested MPKs (Fig. 1a), whereas MPK4 and MPK6 showed a weaker substrate preference to HD2B when compared with the artificial substrate MBP (Fig. 1a). To control our MAPK antibody specificity, we performed kinase assays with protein extracts from mpk3, mpk4, mpk6, and mpk7 mutants and their respective wild-type (WT) lines after immunoprecipitation with the corresponding antibodies (Additional file 1: Figure S6). MPK7 antibody was used as a negative control as this kinase is not activated upon flg22 treatment. The results indicated that the antibodies specifically immunoprecipitated the respective kinases (Additional file 1: Figure S6). To confirm the substrate specificity of the phosphorylation of HD2B by MPK3, we also tested several additional substrates that were identified in the phosphoproteome screen. As shown in Additional file 1: Figure S7, PI-4Kß1 (At5g64070) was specifically phosphorylated by MPK3 and MPK6, SCF (At5g13300) was not phosphorylated by any of the three immune MAPKs, MSL9 (At5g19520) and GOS12 (At2g45200) were preferentially phosphorylated by MPK6. To determine whether MPK3 phosphorylated the previously identified HD2B in vivo sites T249 and S266 , the T249 amino acid residue was mutated individually to glutamate or in combination with S266 to aspartate residues (HD2B-T249E and HD2B-T249E/S266D, respectively; denominated further as HD2B-ED), and these mutant HD2B proteins together with WT HD2B were tested for their ability to serve as MPK3 substrates (Fig. 1b). In comparison to WT HD2B, HD2B-T249E showed reduced and HD2B-T249E/S266D complete lack of phosphorylation by MPK3, indicating that the in vivo phosphorylation sites T249 and S266 can be targeted by MPK3 in vitro (Fig. 1b).
To better understand where the MPK3 and HD2B interaction takes place, we investigated the subcellular localization of the HD2B and MPK3 proteins. A GFP-HD2B fusion protein was localized in the nucleolus [14, 15], whereas MPK3 was observed in both cytoplasmic and nuclear cell compartments (Fig. 1c and Additional file 1: Figure S8). Using indirect immunofluorescence staining, we found that MPK3 and HD2B were co-localized in the nucleoli of the same plant tissue (Fig. 1d). To determine if MPK3 and HD2B can interact in vivo, we performed a bimolecular fluorescence complementation assay (BiFC). To this end, the MPK3 and HD2B complementary DNAs (cDNAs) were inserted into binary vectors, containing the split yellow fluorescent protein (YFP) N-terminal fragment (YFPn) and the C-terminal fragment (YFPc), respectively. When YFPc-HD2B and YFPn-HD2B were expressed together in transiently transformed Nicotiana benthamiana leaves, YFP fluorescence was observed inside nuclei revealing an in planta interaction between the two proteins (Fig. 1e, upper panels). No fluorescence was observed when each tagged protein was expressed separately with the corresponding empty vector (Fig. 1e, middle and bottom panels). Because both MPK4 and MPK6 showed a capacity to phosphorylate HD2B, but with a weaker preference than MPK3, we tested their interaction with HD2B by BiFC. As expected we found that MPK4 and MPK6 can also interact with HD2B but with a lower efficiency (Additional file 1: Figure S9). Due to the fact that MPK3 showed the highest phosphorylation activity we decided to focus our study on this specific MAPK. To confirm the MPK3/HD2B interaction in vivo we performed co-immunoprecipitation experiments. First, protein extracts obtained from protoplasts transiently expressing HD2B-c-Myc and MPK3-HA tagged proteins were immunoprecipitated with an anti-c-Myc antibody and analyzed by immunoblotting with anti-c-Myc or anti-HA antibodies. We observed that HD2B interacts with MPK3 (Additional file 1: Figure S10A) consistent with our in vitro kinase assays (Fig. 1a). Second, co-IP assays performed with transgenic plants expressing both GFP-HD2B and MPK3-c-Myc tagged proteins (Additional file 1: Figure S10B) confirmed our BiFC assays (Fig. 1e) that HD2B interacts with MPK3 in vivo. Taken together, our results showed that MPK3 interacts with and phosphorylates HD2B.
The MPK3-HD2B module controls the transcription of biotic stress response genes by modulating H3K9ac levels
To gain insight into the role of HD2B and its interaction with MPK3 in gene regulation and plant defense, we examined the gene expression consequences of HD2B and MPK3 loss-of-function mutants in resting conditions and in flagellin-challenged plants. To this end, we performed RNA-sequencing (RNA-seq) analyses on hd2b and mpk3 mutants and WT plants (Additional files 2 and 3: Tables SI and SII) either mock-treated for 30 min or treated for the same period of time with 1 μM of flg22. First, when compared to WT, we observed that 1714 genes were deregulated in hd2b with a fold change > 2 and p value < 0.05. From these 1714 genes, 74% were upregulated and 26% downregulated (Fig. 2a). This pattern is consistent with the role of HD2B as a repressor. A functional annotation of hd2b upregulated genes under control conditions showed a clear enrichment in genes involved in transcription (Fig. 2a and Additional file 1: Figure S11). By contrast, an annotation of downregulated genes showed enrichment in genes involved in transport activity (Additional file 1: Figure S12A and S12B). Because HD2B is a histone deacetylase and as such a putative repressor of gene expression, we focused our subsequent analyses on the set of genes that was upregulated in the hd2b mutant. Comparative analysis between the hd2b and mpk3 transcriptomes revealed an overlap among the identity of the deregulated genes in the two mutants. Indeed, a total of 414 genes were commonly upregulated in the two mutants (Fig. 2b) representing 32% of the mpk3 upregulated genes. This result was not surprising because MPK3 likely has many substrates in addition to HD2B that contribute to the gene expression modifications observed in the transcriptome of our mutants. Reciprocally, HD2B likely is involved in diverse cellular processes that include, but are not limited to, pathogen defense. Gene ontology (GO) analysis of this class of genes revealed a significant enrichment in RNA biosynthetic processes and transcription (Fig. 2c).
To put these responses in the context of defense responses, the transcriptome results of untreated hd2b and mpk3 mutants were compared with data obtained on WT plants after flagellin treatment, which allowed us to distinguish 150 genes that were upregulated in the mpk3 and hd2b mutants and upon flagellin treatment (Fig. 2d), suggesting that the HD2B-MPK3 module normally represses the expression of these genes but that the inhibitory effect is alleviated during pathogen defense.
To further dissect the mechanisms by which MPK3 and HD2B regulate these genes, and based on the changes in their expression levels from the RNA-seq data, we selected four genes, AT2G31290, AT5G03350, AT3G11380, and AT2G28180, for further analysis. Quantitative reverse transcription polymerase chain reaction (RT-qPCR) confirmed the RNA-seq results indicating an over-expression of AT2G31290, AT5G03350, AT3G11380, and AT2G28180 in mpk3 and hd2b mutants (Fig. 2e).
Because HD2B is a histone deacetylase, we hypothesized that HD2B could be directly interacting with the gene loci to repress their transcription. To test if HD2B directly binds the four regulated genes, we performed chromatin immunoprecipitation (ChIP)-qPCR on seedlings expressing a GFP-HD2B fusion protein under the control of the endogenous HD2B promoter (pHD2B::GFP-HD2B). As a negative control, we included the AT1G07690 gene that is upregulated in hd2b mutants (Additional file 1: Figure S13A) but was not identified as a HD2B target by ChIP-seq. A marked enrichment of HD2B binding on each of the four candidate genes but not on AT1G07690 was found in the absence of flagellin treatment (Fig. 2f and Additional file 1: Figure S13B). Therefore, the HD2B protein is associated with chromatin and binds to these four genes, suggesting that HD2B directly regulates their transcription.
Acetylation and deacetylation are dynamic processes that can be rapidly reprogrammed, depending on the signals received by the cell. To determine the level of H3K9 acetylation, we performed ChIP-qPCR assays on these representative genes. Interestingly, we found hyper-acetylation on these four genes in both mpk3 and hd2b mutants when compared with WT (Fig. 2g). Altogether, our results indicate that the HD2B-MPK3 module is required to regulate the basal expression of a subset of genes. In hd2b and mpk3 mutants, the absence of deacetylation results in the constitutive induction of these genes (Fig. 2e). Thus, in addition to responding to pathogen stress, MPK3 activity appears to be required under normal conditions to promote HD2B-directed histone deacetylation to repress unwanted expression of these genes (Figs. 2e–g).
HD2B is involved in pathogen defense
Because the MPK3-HD2B module regulates pathogen response genes, we investigated the direct relevance of HD2B during biotic stress. In Arabidopsis, MPK3 has been shown to negatively regulate the basal activity of defense genes in the absence of MAMPs, but to be required for full expression of defense-related genes upon pathogen challenge . Because MPK3 interacts and phosphorylates HD2B, we hypothesized that HD2B might also be implicated in defense against pathogens. To test this hypothesis, hd2b mutants were challenged with the non-pathogenic bacterial strain Pseudomonas syringae pv. tomato hrcC- that carries a mutation in the type 3 secretion apparatus and is hence deficient in effector deployment. Mutants showed increased sensitivity compared with WT plants (Fig. 3a). As a control, we also tested the hd2b mutant complemented by expression of HD2B under its own promoter (pHD2B::GFP-HD2B). hd2b complemented with pHD2B::GFP-HD2B showed infection levels of Pst hrcC- that were similar to WT col-0 plants. We also compared the resistance of the hd2b mutant and HD2B over-expressing lines with WT Col-0 plants upon infection to virulent Pseudomonas syringae pv. tomato PstDC3000. hd2b mutants were also more susceptible than WT, whereas plants over-expressing HD2B under the constitutive 35S cauliflower mosaic virus promoter (35S::HD2B) were more resistant to Pst DC3000 (Additional file 1: Figure S14). These results show that HD2B is implicated in the defense against pathogens probably by regulating defense genes. We hypothesized that HD2B is not involved in MAMP signaling upstream of the MAP kinase pathway. To confirm this, we first analyzed the MPKs activation after flg22 treatment in both WT and hd2b mutant and no changes were observed (Additional file 1: Figure S15). This result suggested that the signaling cascade is not affected in the hd2b mutant. We next asked if flagellin sensitivity is affected in this mutant. We did not observe any difference in growth inhibition (Additional file 1: Figure S16), suggesting that HD2B is not involved in MAMP-induced growth inhibition which is MAPK-independent.
The HD2B phosphorylation status determines its sub-nuclear localization and choice of target genes
Having established that HD2B contributes to pathogen defense and that it can be phosphorylated by MPK3, we investigated in more detail the consequences of MPK3-dependent HD2B phosphorylation in response to flagellin. MPK3-mediated phosphorylation could affect HD2B activity in several ways, such as modulating its catalytic activity, its stability, or cellular localization. To discern if MPK3-mediated phosphorylation affects any of these parameters, we first investigated the subcellular localization of HD2B after flagellin-induced MPK3 activation. In pHD2B::GFP-HD2B leaves, GFP-HD2B was observed to be relocated from the nucleolus to the nucleoplasm (diffuse and speckle pattern) uniquely after flg22 recognition but not after mock treatment (Fig. 3b). This relocalization was confirmed also to occur in protoplasts (Additional file 1: Figure S17), indicating that HD2B moves from the nucleolus to other sites in the nuclear compartment in response to flg22 treatment. To examine the role of HD2B phosphorylation by MPK3 on the sub-nuclear localization of HD2B, we analyzed the localization of WT and mutated HD2B. Consistent with phosphorylation playing a role in HD2B localization, YFP-HD2B-ED (HD2B-T249E/S266D), a mutated form of HD2B that mimics a constitutively phosphorylated version of the HD2B protein, accumulated in the whole nucleus of Arabidopsis protoplasts even in the absence of flagellin treatment (Additional file 1: Figure S18). We confirmed this nuclear relocalization in stably transformed Arabidopsis lines expressing pHD2B::GFP-HD2B-ED (Fig. 3c and Additional file 1: Figure S19) suggesting that flg22-triggered MPK3 phosphorylation of HD2B induces the nuclear relocalization of HD2B. To test the role of HD2B phosphorylation in vivo, we introduced a non-phosphorylatable HD2B version HD2B-AA (HD2B-T249A/S266A) in the hd2b mutant. This mutated version of the protein accumulated only in the nucleolus (Fig. 3c). We next tested the sensitivity of HD2B-AA complemented hd2b lines to P. syringae hrcC-. These plants displayed enhanced sensitivity to the pathogen when compared to plants expressing WT HD2B (Fig. 3a), further confirming that phosphorylation of HD2B is required for proper activation of defense mechanisms. Consistently, the HD2B-AA protein fails to restore a normal expression level to HD2B target genes that are constitutively activated in the hd2b mutant (Additional file 1: Figure S20).
To determine if the change in HD2B localization has an impact on the genes that this HDAC targets, we performed ChIP-seq in pHD2B::GFP-HD2B seedlings after a 30-min flagellin or mock treatment. First, we controlled the HD2B enrichment after ChIP-seq (Fig. 4a) and then performed a peak detection using MACS2. In control mock conditions, HD2B bound 5460 genomic loci, whereas it bound 8149 loci in flg22-treated plants (Additional file 2: Table SIV), indicating that flg22 treatment induced a recruitment of HD2B on chromatin. These differences could be attributed to a relocalization of the protein rather than to modifications of its accumulation since the expression of HD2B remains constant upon flagellin treatment (Additional file 1: Figure S21). 4431 HD2B target sites overlapped between mock- and flg22-treated plants, revealing a redistribution of HD2B on chromatin.
In mock-treated and flg22-treated conditions, most HD2B peaks were found approximately 300 base pairs (bp) downstream of the transcription start site (TSS) (Fig. 4c and d), and most target genes (84% in mock-treated and 78% in flg22-treated plants) were protein-coding genes (Fig. 4b). Several of these ChIP-seq identified targets were confirmed by ChIP-qPCR in this study (Fig. 2f). A GO analysis of the protein-coding HD2B-targeted genes revealed a significant enrichment for genes involved in defense response under mock conditions, whereas HD2B targeted genes after flg22 treatment were mainly involved in plastid organization and interestingly these genes are downregulated after a pathogen attack (Fig. 4e and Additional file 1: Figure S22).
To correlate HD2B binding to changes in gene expression, we compared the ChIP-seq results obtained for HD2B binding sites with the transcriptome of the hd2b mutant. In total, 277 of HD2B direct targets were upregulated in the hd2b mutant, whereas only 79 were downregulated (Additional file 1: Figure S23). This low level of correlation could be explained by the redundancy in the HD2 protein family. Together, these results indicate that after pathogen recognition, MPK3 modulates HD2B localization and activity to promote its dynamic redistribution on a new set of target genes.
Modulation of H3K9ac levels is a hallmark of the flagellin response
Because HD2B targets are a highly varied set of genes before and after flg22 treatment, it is possible to assume that changes in histone acetylation are an important component of the cellular response to pathogens. To analyze the contribution of histone acetylation to the flg22-induced immune response, we analyzed the H3K9ac landscape before and after 30 min of flg22 treatment using a ChIP-seq approach. We identified about 15,000 peaks both in control and flg22-treated plants (Additional file 3: Table SIII). Most of these peaks were observed at core promoters or protein-coding genes where they predominantly localized to the first nucleosome after the TSS (Additional file 1: Figure S24A–E). According to the agriGO toolkit, genes harboring the H3K9ac mark were associated with developmental processes, responses to abiotic and biotic stimulus, signaling cascades and responses to hormone stimuli, both in control conditions and after flg22 treatment (Additional file 1: Figure S24F).
We next compared the two ChIP-seq datasets using seqMINER and diffReps, which allows qualitative and quantitative comparisons to be made between a reference set of genomic positions and multiple ChIP-seq datasets . Using this approach, we observed three classes of genes. The first class, consisting of 731 genes (Additional file 3: Table SIII), was hyper-acetylated after flg22 treatment (Fig. 5a and b). GO analysis of this class revealed a significant enrichment in genes involved in defense and immune responses and phosphorylation (class I, Fig. 5c and Additional file 1: Figure S26A). The second class of 13,159 genes did not display any obvious changes in H3K9ac marks (class II, Fig. 5a). The third class totaling 787 genes was hypo-acetylated after flg22 treatment (class III, Fig. 5a and e and Additional file 3: Table SIII) and mainly consisted of genes involved in chloroplast and plastid organization and in metabolic processes (Fig. 5f and Additional file 1: Figure S26B). Gene expression analysis revealed that most hyper-acetylated genes were induced (Fig. 5d) and most hypo-acetylated genes were repressed by diverse biotic stresses (Fig. 5g), providing strong evidence for a positive correlation between the acetylation and the expression level of biotic stress responsive genes. Consistently, comparing these results with the transcriptome of WT plants treated with flg22 revealed that about 98% of hyper-acetylated genes were induced by flg22 treatment, whereas only 2% were repressed, and about 64% of hypo-acetylated genes were repressed by flg22 treatment, whereas only 36% were induced (Additional file 1: Figure S27), clearly showing that H3K9ac is associated with gene activation. Our results indicated that flg22-induced genes were regulated mainly by de novo acetylation, whereas flg22-repressed genes were regulated by dynamic deacetylation, highlighting the role of HDAC proteins in the control of flg22-regulated genes.
HD2B is a major contributor of flg22-induced changes in histone acetylation
Because histone acetylation appears to be a major factor controlling the transcriptional reprogramming induced by pathogen recognition, we investigated whether HD2B could be a central contributor to this process. Among the 787 genes hypo-acetylated after flg22 treatment, 60% were HD2B targets, indicating that HD2B plays an important role in the transcriptional reprogramming induced by pathogens (Fig. 6a). As expected, among the 731 hyper-acetylated genes only 14% were HD2B targets (Fig. 6a).
To clarify the mechanism by which the MPK3-HD2B module acts on genes after flg22 treatment, we selected four HD2B target genes that were hypo-acetylated after flg22 treatment. We first confirmed their acetylation status (Fig. 6b) and then verified that they were bound by HD2B after flg22 treatment (Fig. 6c). Second, we analyzed their expression and acetylation levels before and after flg22 treatment in WT plants and compared them with the levels in mpk3 and hd2b mutants (Fig. 6d). As expected, in WT control plants, hypo-acetylated genes were repressed and were further hypo-acetylated after flg22 treatment, but neither gene repression nor a change in acetylation occurred in the mpk3 and hd2b mutants, confirming the role of the MPK3-HD2B module in the repression and acetylation dynamics of these genes (Fig. 6b–d). Furthermore, to functionally prove that HD2B phosphorylation plays a major role in this process, we analyzed the expression levels of these specific four genes by qRT-PCR in our HD2B-AA line. We observed that this repression is not restored in the HD2B-AA line, providing further evidence for the functional relevance of HD2B phosphorylation (Additional file 1: Figure S28).
Plants are sessile organisms that constantly have to respond to changes in environmental conditions such as biotic or abiotic stress. In the case of biotic stresses, plants have adapted the capacity to recognize pathogens through MAMPs via specific receptors. The multi-faceted responses of plants to pathogens have been examined extensively and these studies have revealed that receptor-mediated pathogen recognition triggers MAP kinase signaling cascades, which, through the activation and repression of large gene sets, ultimately result in the establishment of MAMP-triggered immunity.
Previous studies have identified several transcription factors targeted by MAP kinases. For example, the bZIP (basic leucine zipper) transcription factor VIP1 (VirE1-INTERACTING PROTEIN 1) is specifically phosphorylated by MPK3 upon flg22 treatment and activates a number of defense-related genes such as PR1 (PATHOGENESIS-RELATED 1) [10, 18]. However, extensive changes in the transcriptional program of cells likely rely not only on the activation of specific transcription factors but also on chromatin modifications that can dynamically and reversibly regulate gene expression. Consistently, we show here that flagellin-activated MAP kinases interact with and phosphorylate the histone deacetylase HD2B to control the expression level of biotic stress-regulated genes through modulation of the H3K9ac histone mark.
To date, a handful of histone modifiers have been implicated in plant innate immunity. For example, salicylic acid (SA) signaling plays an essential role in plant pathogen resistance and is controlled partially by the HDAC SIRTUIN2 (SRT2), which represses the expression of several SA biosynthetic genes such as PAD4 and SID2 . Consistently, srt2 mutant Arabidopsis plants were reported to be more resistant to pathogen infection than WT control plants, whereas an SRT2 over-expressing line was more susceptible. In addition, it was reported in Arabidopsis that mutations in the HDAC HDA19 result in enhanced basal expression of several biotic responsive genes  and improve tolerance to P. syringae , although contradictory results have been described previously . Moreover, the rice HDAC HDT701 negatively regulates innate immunity by directly binding and modulating the histone H4 acetylation levels of PRR and defense-related genes .
Here, we measured the resistance of hd2b mutants to the both virulent and non-pathogenic Pseudomonas syringae strains and concluded that hd2b mutants were more susceptible than WT control plants, whereas HD2B over-expressing plants were more resistant. Hence, different HDACs can lead to similar or contrasting outcomes for plant immunity likely depending on the genes they target; our results and those reported by others clearly highlight the crucial contribution of chromatin regulation to pathogen defense.
Although several studies provide evidence for the involvement of histone modifiers in pathogen response, the way in which their activity is connected to biotic stress signaling cascades has not been unraveled. In this work, we connected MAPK-induced signal transduction to HDAC-dependent histone modifications. Indeed, our results show that H2DB is a direct target of MAPKs and that HD2B regulates a large number of genes involved in pathogen defense. Although our study is the first to describe a general histone modifying protein as a direct target of a MAPK in plants, such mechanisms seem to be highly conserved in eukaryotes. Indeed, in animals, MAPKs drive histone modifications and direct chromatin remodeling. For example, phosphorylation of the histone acetyl-transferase p300 by the MAP kinase ERK2 promotes its localized histone binding . Likewise, ERK1 and ERK2 modulate the assembly of chromatin remodeling complexes and, thereby, control the expression of vitamin D-responsive genes .
The Arabidopsis genome codes for 18 HDACs, which fall into four groups [26, 27]: class I and II HDACs correspond to yeast RPD3 and HDA1, respectively, class III enzymes are NAD-dependent HDACs related to yeast SIR2, and class IV HDACs are plant-specific and consist of four members called HD2A-D . The four class IV HDACs all contain a conserved N-terminal catalytic and a central acidic domain and a divergent C-terminal region that suggest specific functioning of the different proteins . In contrast to the other three Arabidopsis HD2s, HD2B does not contain a Zn-finger motif in its C-terminus. However, with the exception of Arabidopsis HD2A, the two MAPK-targeted C-terminally located HD2B in vivo phosphorylation sites are highly conserved in other species, including the equivalent rice and maize HD2s . Our data indicate that phosphorylation of these sites is important for the molecular function of HD2B and that the phosphorylation status of the two residues determines the intra-nuclear localization of the HD2B. Possibly, a nucleolar protein could anchor unphosphorylated HD2B to this sub-nuclear compartment and MAPK-directed phosphorylation of HD2B could release HD2B to enter the remaining nuclear compartment and associate with novel chromosomal loci. In this way, the nucleolus could serve as a stocking center to prevent HD2B from potentially associating with certain labile nuclear sites and could provide a protein synthesis-independent mechanism to make HD2B rapidly available for gene regulation upon signaling.
Our results show that MPK3 likely modulates HD2B activity by altering its nuclear distribution. Indeed, flagellin treatment induces the relocalization of HD2B in an MPK3-dependent manner, and this effect can be obtained by expressing a phosphomimetic version of HD2B. Consistently, control of the intracellular location has been described in yeast and mammals as an important regulatory mechanism of HDAC activity [28, 29]. As expected from their ability to deacetylate histones, HDACs are found generally in the nucleus of most organisms. Nevertheless, in mammals, the Class IIa HDACs were found both in the cytoplasm and in the nucleus. Indeed, specific phosphorylation of HDAC7 induces its nuclear exclusion, rendering it unable to impact transcription . In the case of HD2B, phosphorylation leads to a sub-nuclear relocalization from the nucleolus to the nucleoplasm in specific chromatin regions. In agreement with the flg22-induced relocalization of HD2B, our ChIP-seq analyses revealed a drastic shift in the HD2B chromatin targeted sites upon flagellin treatment. In mammals, it has been reported that the phosphorylation of HDAC2 leads to its redistribution on chromatin and increased recruitment to promoters [31, 32], suggesting that HDAC phosphorylation could be a mechanism conserved from mammals to plants to drive their redistribution in response to a stimulus.
HD2B is a HDAC and is predicted to function as a repressor of gene expression. As expected, a proportion of HD2B target genes were constitutively upregulated in hd2b mutants (Additional file 4). Furthermore, 64% of the genes that are hypo-acetylated after flagellin treatment are direct HD2B targets, indicating that HD2B contributes substantially to MAMP-triggered transcriptional reprogramming (Additional file 5). The study of pathogen-induced histone modifications is a fledgling field of research. In mammals, although four bacteria have been reported to modulate histone acetylation levels upon infection, the underlying mechanisms are not known [33,34,35,36]. In plants, modulation of histone acetylation levels in response to biotic stress has also been described, but again the underlying mechanisms remained unclear [19, 21, 23, 37, 38]. Our data show that genome-wide modulation of H3K9ac levels is an early chromatin response to a biotic stress. Histone acetylation is a labile chromatin mark [39, 40] and, as such, is a fast and reversible post-translational process that can allow plants to rapidly modulate gene expression. Thus, this mark is ideally suited for biological programs that require an immediate response to fluctuating environmental conditions. In line with this characteristic, we observed a strong correlation between upregulation of gene expression and hyper-acetylation, whereas downregulation correlated with hypo-acetylation. Altogether, our data suggest that modulation of the genome-wide H3K9ac landscape is a hallmark of the flg22 response and that HD2B is a key player in this process. However, our studies do not exclude the involvement of additional HDACs in this process.
Finally, our detailed analysis of the expression levels, HD2B binding, and histone acetylation of HD2B target genes under control and flg22 treatment conditions revealed that HD2B has several roles. We have shown that HD2B regulates the basal expression level of a subset of genes in the absence of pathogen challenge, probably by acting as a counter-balance to the activity of HATs. Indeed, expression of these HD2B-targeted genes in unchallenged plants is constitutively high in hd2b mutants. These genes can be distributed into two classes after flagellin treatment: class I genes are repressed by flagellin due to increased HD2B recruitment on these sites, whereas class II genes are induced by flagellin due to eviction of HD2B from these loci.
In summary, our work defines the first example of a MAP kinase-regulated chromatin mechanism and details how MAMP-triggered MAP kinase signaling regulates global changes in the chromatin landscape. As such, it sets the stage for other large-scale studies examining the contribution of protein kinase-mediated chromatin regulation in plants and mammals.
Plant material, growth conditions, and treatments
T-DNA insertion lines hd2b (At5g22650) from SAIL collection  Sail_1247_A02 (Additional file 1: Figure S29) and mpk3 SALK_151594  were obtained from the Nottingham Arabidopsis Stock Centre (NASC).To produce transgenic plants expressing a 35S::GFP:MPK3 construct, the coding region of MPK3 (At3g45640) was amplified from cDNA of Col-0 by PCR using the following primers: forward, 5′gc gga tcc atg aac acc ggc ggt ggc3′, reverse 5′gc act agt cta acc gta tgt tgg ctt gag3′. The restriction enzyme sites, set in bold, that had been added to the primers were utilized to ligate the open reading frame to BamHI/SpeI sites of pCAT-GFP, resulting in p35S::GFP-MPK3. For the generation of stably transformed A. thaliana plants, the expression cassette was excised by Sse8387I and cloned into the PstI site of the binary vector pCB302 . This vector was transformed to A. tumefaciens strain GV3101 by electroporation. A. thaliana Col-0 plants were transformed by floral dip method according to Clough and Bent . Transformed plants were selected on BASTA.
Both pHD2B::GFP-HD2B and 35S::GFP-HD2B constructs were made by using Gateway technology (Invitrogen). The CDS of HD2B was PCR amplified from cDNA with forward primer 5′-caccATGGAGTTCTGGGGAGTTGC-3′ and reverse primer 5′-AGCTCTACCCTTTCCCTTGC-3′ using Phusion High Fidelity DNA polymerase (New England BioLabs). The putative promoter (1.055 kb upstream of the start codon) was PCR amplified from genomic DNA with forward primer 5′-caccGTTTTGGATCTGCAGACAAGG-3′ and reverse primer 5′-TGTTGTTGAACGAGGAAGAGAG-3′. Both PCR fragments were cloned into pENTR/D-TOPO (Invitrogen). The promoter was then recloned from the pENTR/D-TOPO into a pENTR4-1 vector (Invitrogen) in front of a GFP using the NotI and AscI restriction sites. The pENTR4-1 containing the HD2B promoter and GFP, the pENTR-D-TOPO-AtHD2B vector and a pENTR2-3 vector containing a CaMV 35S terminator were recombined by a multisite gateway reaction (Invitrogen) into the binary destination pBnRGW vector. This is a modified vector based on pKGW  in which the kanamycin resistance was replaced with basta resistance and the NAP::DsRed expression cassette from pFluar 101  was introduced for easy selection of red fluorescent transformed seeds. To create 35S::GFP-HD2B construct the same strategy was used but pENTR4-1 contained CaMV35S promoter. A. tumefaciens strain C58 mediated transformation for both constructs was performed as described by Bechtold and Pelletier . Seed-specific expression of red fluorescent protein DsRed permitted the identification of mature transformed seeds by fluorescence stereo microscopy (Leica).
For flagellin experiments, seeds were surface-sterilized by treatment with bayrochlore and then soiled in sterile half-strength MS liquid medium, placed for 2–4 days at 4 °C to obtain homogeneous germination, and plants were grown in chambers at 20 °C in long-day (16 h of light) conditions. After 14 days, flg22 peptide was added in the medium to a final concentration of 1 μm when necessary.
Infection assays in Arabidopsis rosette leaves with Pseudomonas syringae pv. tomato PstDC3000 were performed as described previously  on leaves of 4–5-week-old plants grown on soil in environmental chambers at 22 °C under short-day conditions (8 h of light).
For Pseudomonas pathogen assays in Arabidopsis seedlings, 0.5 MS plates containing 14-day-old Arabidopsis plantlets were flooded with a bacterial suspension of the hrcC- mutant (defective in type III secretion system) of Pseudomonas psyringae pv. tomato DC3000 for 3 min. The bacteria were grown ON at 30 °C in LB broth + Rif and resuspended in 10 mM MgCl2 after several washings for the removal of the media and the antibiotic. The final concentration of the suspension was adjusted to an OD600 = 0.1 and supplemented with 0.025% Silwet L-77. After discarding the bacterial suspension, the plates were sealed with micropore and put back in the incubation room. Day 0 samples were taken 2 h post-inoculation and day 2 samples after 48 h. For the sampling process, the rosettes of four plantlets were cut and weighed together, registering each time the fresh weight in milligrams. Afterwards, the rosette surface was sterilized by a 5-s wash in 70% ethanol, followed by two washes in sterile water. The samples were ground in order to release the bacteria and serial dilutions performed. A total of 10 μL of each dilution was plated on plates with LB + Rif and incubated for two days at 30 °C. The growth of the bacterial population inside the plantlets was determined by calculating the total amount of bacterial CFU per mg of fresh weight. Each experiment was performed in quintuplicate, analyzing eight samples per accession on each experiment.
Full-length HD2B (At5g22650) was amplified and cloned into pGEX4T-1. For mutagenesis of HD2B, PCR was performed on plasmids with Pfu Ultra. PCR mixtures were digested for 2 h with DpnI and transformed into Escherichia coli. Clones were sequenced and each mutation was transformed into E. coli BL21. GST protein expression was induced for 4 h at 37 °C with 1 mM IPTG and purified according to the manufacturer’s protocol. For producing HD2B T249E and S266D, the following primers were used: GGAGGACACACCGCCGAACCACACCCAGCT, AGCTGGGTGTGGTTCGGCGGTGTGTCCTCC, GTGAATGCTAACCAGGACCCCAAGTCTGGA, and TCCAGACTTGGGGTCCTGGTTAGCATTCAC, respectively.
Radioactive kinase assays were performed with MPK3, 4, and 6 immunoprecipitated from extracts of A. thaliana cell cultures that were treated for 10 min with 1 μM flg22. Non-radioactive kinase assays for phosphosite mapping were performed the same for 45 min at 30 °C.
Protein localization in protoplasts
Full-length HD2B and HD2B-DD mutated forms were cloned in a pGreen plasmid behind a YFP gene driven by a 35S promoter. Arabidopsis protoplasts were prepared from a suspension culture as described . One day after transformation, flg22 peptide was added in the medium to a final concentration of 2 μM for 15 min when necessary and localization was checked by fluorescence microscopy.
Leaves of plants stably transformed with p35S::GFP-HD2B construct were treated for 30 min with or without 1 μM flg22 and fixed in PFA 4% in PHEM (PIPES 60 mM; HEPES 25 mM; EGTA 10 mM; MgCl2 2 mM pH 6.9) for 1 h at room temperature (apply vacuum 20 min to facilitate uptake of the fixation solution). Seedlings were washed 5 min in PHEM and 5 min in phosphate buffered solution (PBS) pH6.9 and chopped on a Petri dish in PBS supplemented with 0.1% Triton X-100 (w/v). The mixture was filtered (50 μm) and centrifuged for 10 min at 2000 g. The supernatant was carefully removed and the pellet washed once with PBS, gently resuspended in 20 μL PBS, and a drop was placed on a poly-lysine slide and air dried. Slides were rehydrated with PBS and permealized twice by 10 min incubation in PBST (PBS, 0.1% Tween-20 v/v). Slides were placed in a moist chamber and incubated overnight at 4 °C with primary antibody anti-GFP (Clontech, ref. 632592) in PBST supplemented with BSA (3% w/v). Slides were washed 5 × 10 min in PBST (at RT) and incubated for 1 h at room temperature in the dark with the secondary antibody (A11037 Invitrogen, Alexa Fluor 594 goat anti-rabbit) diluted (1/400 v/v) in PBST, 3% BSA. Slides were washed 5 × 10 min in PBST and then mounted with a drop of Vectashield with DAPI and observed as described for pollen mitosis analysis with the suitable cube fluorescence filters (BP340-380, DS 400, BP 450-490 for DAPI) (BP570-590, DS595, BP 605-655 for A594).
For whole mount immunofluorescence labeling, six-day-old pHD2B::GFP-HD2B seedlings were labeled according to the method of Sauer et al.  with small modifications: seedlings were fixed in 1.5% paraformaldehyde and 0.5% glutaraldehyde in 0.5 MTSB buffer (50 mM K-PIPES, 5 mM MgSO4-7H20, 5 mM EGTA) at pH 6.8; cell-wall digestion enzyme mixture contained 1% meicelase, 1% cellulase, and 1% macerozyne in PBS; samples were incubated with rabbit polyclonal anti-AtMPK3 (SIGMA-ALDRICH) primary antibody diluted 1:350 in PBS supplemented with 2% BSA at 4 °C overnight, and subsequently with secondary antibody Alexa-Fluor 546 goat anti-rabbit IgGs (H + L) (Invitrogen) diluted 1:500 in PBS containing 2% BSA for 3 h (1.5 h at 37 °C and 1.5 h at room temperature). Finally, samples were counterstained with DAPI and mounted in one drop of 0.1% [w/v] para-phenylenediamine prepared in 90% [v/v] glycerol in PBS. Microscopic analysis of immunolabeled samples was performed using a Zeiss LSM710 (Carl Zeiss Jena, Germany). DAPI was excited at 405 nm and emission was detected between 410 and 476 nm. GFP was excited at 488 nm and emission was detected between 500 and 535 nm. Alexa 546-conjugated antibody was excited at 561 nm and fluorescence was detected between 566 and 591 nm. The post-processing of images was done using ZEN 2010 software, Photoshop 6.0/CS, and Microsoft PowerPoint.
LR recombinations of HD2B and MPK3 coding sequences in pENTR3C were done with split-YFP destination vectors pBiFC3 and pBiFC2 that allow N-terminal fusion with the C- and N-terminal YFP moieties, respectively . Recombined vectors were transformed into the Agrobacterium C58C1 strain. Nicotiana benthamiana leaves were agro-infiltrated as previously described . After three days, the YFP fluorescence was visualized using a confocal laser scanning microscope (Leica, Germany).
For microscopy, six-day-old plants stably transformed with p35S::GFP-MPK3 construct were transferred to microchambers mounted between microscopic slides and coverslips with one Parafilm layer as a spacer. The chambers were filled with liquid half-strength MS medium. Microscopic analysis was performed using a Zeiss LSM 710 confocal laser scanning microscope (Carl Zeiss Jena, Germany). All images were acquired using a 40X objective lens (NA 1.42). GFP was excited at 488 nm and detected between 500 and 535 nm. Post-processing of images was done with the aid of Zeiss ZEN software (Ver.2010b) and Microsoft PowerPoint applications.
Two grams of 16-day-old Col-0 seedlings were ground in a mortar with liquid nitrogen and were resuspended in buffer A containing 2.5% Ficoll type 400 (F-4375, Sigma-Aldrich), 5% Dextran (D1662, Sigma-Aldrich), 0.4 M sucrose, 25 mM Tris-HCl pH 7.5, 10 mM MgCl2, 5 mM DTT (D0632, Sigma-Aldrich), protease inhibitors (Complete cocktail, Roche), and phosphatase inhibitors (1 mM NaF, 0.5 mM Na3VO4, 15 mM B-glycerophosphate, 15 mM 4-nitrophenyl phosphate, Sigma-Aldrich chemicals). After a few minutes’ incubation on ice, the samples were filtrated through one layer of 62-μm nylon mesh by centrifugating for 3 min at 212 g at 4 °C. Triton X-100 was then added to a final concentration of 0.5% and the samples were gently mixed and incubated on ice for 15 min. An aliquot was kept at this step, hereafter referred to as Input. The samples were then centrifuged for 5 min at 1500 g at 4 °C and an aliquot of the supernatants was kept at this step, hereafter referred to as Cytoplasmic fraction. The pellets were gently resuspended in buffer B (buffer A and 0.1% Triton X-100) and centrifuged for 5 min at 1500 g at 4 °C. The pellets were then resuspended in buffer A and centrifuged for 5 min at 2000 g at 4 °C. The pellets were finally resuspended in SDS-sample buffer; they are hereafter referred to as nuclear fractions. SDS-sample buffer was also added to the Input and Cytoplasmic fractions and all three kinds of aliquots were heated at 95 °C for 10 min for complete denaturation. Protein samples were resolved by SDS-PAGE at a constant amperage of 15 mA per gel and transferred onto methanol-activated PVDF membranes (GE Healthcare) for 1 h at a constant voltage of 100 V. Blots were blocked with 5% non-fat dry milk in 1× TBST for 1 h and probed with primary antibodies overnight at 4 °C: anti-H3 (Abcam ref. Ab1791, diluted 1:10,000), anti-PEPC (Tebu-bio ref. 100-4163, diluted 1:15,000), and anti-MPK3 (described in , diluted 1:4000). The membranes were washed four times with 1× TBST. Goat anti-rabbit antibodies, at a dilution of 1:20,000 in 5% non-fat dry milk in 1× TBST, conjugated to horseradish peroxidase (A6154, Sigma-Aldrich), were used as secondary antibodies for 1 h at room temperature. The membranes were washed again four times with 1× TBST and the antigen-antibody interaction was detected with enhanced chemiluminescence reagent (ECL Prime, GE Healthcare) using a GeneGnome imaging system (Syngene).
For co-immunoprecipitation assays in protoplasts, protein extracts were prepared from protoplasts as described . For co-immunoprecipitation, samples were extracted in 50 mM Tris-HCl, pH 7.8, 150 mM NaCl, 1 mM EDTA, 0.1% Nonidet P-40 (NP-40), and proteinase inhibitor mix (Complete EDTA free, Roche). Three independent protoplast transformations of 100 μL each were pulled and used for protein extraction. Protein extracts were precleared with 15 mL of protein A–Sepharose beads for 2 h at room temperature, then immunoprecipitated overnight in the presence of antibodies with 25 μL of beads. Samples were washed three times with wash buffer (50 mM Tris-HCl pH 7.8, 150 mM NaCl, 5 mM EGTA, 5 mM EDTA, and 0.1% Tween 20) and subjected to immunoblotting.
For co-immunoprecipitation experiments in planta, transgenic lines expressing p35S::GFP-HD2B or p35S::MPK3-c-Myc constructs were crossed together to obtain plants expressing both constructs. Plant leaves were ground to a fine powder in liquid nitrogen. Cells were homogenized and lysed in Chris buffer containing 50 mM Tris-HCl pH 8.0, 0.5% NP-40, 200 mM NaCl, 0.1 mM EDTA, 10% glycerol, 10 mM NEM, and protease inhibitor mix (Complete EDTA-free, Roche). GFP-HD2B and MPK3-c-Myc were then immunoprecipitated using a polyclonal anti-GFP antibody (Clonetech) and a polyclonal anti-c-Myc antibody (Sigma), respectively. Immune complexes were collected by incubation for 2 h at 4 °C with Protein A/G Ultralink Resin (Thermo Scientific) and washed six times in lysis buffer. Immunoprecipitated proteins were detected by immuno-blot using a polyclonal anti-c-Myc antibody (Sigma).
About 200 mg of 14-day-old seedlings were ground with liquid nitrogen and resuspended in 500 μL of a buffer containing 50 mM Tris-HCl pH 7.5, 150 mM NaCl, 0.1% NP40, 5 mM EGTA, 0.1 mM DTT (Sigma-Aldrich chemicals), protease inhibitors (Complete cocktail, Roche), and phosphatase inhibitors (1 mM NaF, 0.5 mM Na3VO4, 15 mM β-glycerophosphate, 15 mM 4-nitrophenyl phosphate, Sigma-Aldrich chemicals). The suspension was centrifuged at 20,000 g for 15 min at 4 °C, and the supernatant was collected. Protein quantification was carried out by Bradford method (Protein Assay Kit, Thermo Fisher Scientific) and the normalized protein amounts of all samples were denatured with SDS-sample buffer by boiling them at 95 °C for 10 min. Protein samples were resolved by SDS-PAGE and transferred onto PVDF membranes (Bio-Rad). Blots were blocked with 5% BSA in 1× TBST and then probed with Phospho-p44/42 MAPK (Erk1/2) (Thr202/Tyr204) (D13.14.4E) XP rabbit monoclonal antibody (#4370, Cell Signaling), hereafter referred to as anti-pTpY antibody. Goat anti-rabbit antibodies conjugated to horseradish peroxidase were used as secondary antibodies (A6154, Sigma-Aldrich). The antigen − antibody interaction was detected with chemiluminescent reagents (Clarity ECL substrate, Bio-Rad) using an imaging system (ChemiDoc MP System, Bio-Rad). Coomassie blue staining of blots was then carried out for protein visualization.
Gene expression analysis
Total RNA was extracted from seedlings with the RNeasy MiniPrep kit (Qiagen), according to the manufacturer’s instructions. First strand cDNA was synthesized from 2 μg of total RNA using Improm-II reverse transcriptase (A3802, Promega) according to the manufacturer’s instructions. of the synthesized cDNA, 1/25th was mixed with 100 nM each of primer and LightCycler® 480 Sybr Green I master mix (Roche Applied Science) for qPCR analysis. Products were amplified and fluorescent signals acquired with a LightCycler® 480 detection system. The specificity of amplification products was determined by melting curves. UBQ10 was used as internal control for signals normalization. Exor4 relative quantification software (Roche Applied Science) automatically calculated relative expression level of the selected genes with algorithms based on ΔΔCt method. Data were from duplicates of at least two biological replicates. The sequences of primers can be found in Additional file 6: Table SV.
Transcriptomic analyses were performed RNA-seq at the URGV plateform (Evry, France). Three biological replicates were run for each condition.
For GO analyses, the GO Analysis Toolkit and Database for Agriculture Community  was used.
ChIP assays were performed on 14-day-old in vitro seedlings using anti-GFP (Santa Cruz), IgG control (Millipore), or anti-H3K9ac (Millipore) antibodies, using a procedure adapted from Gendrel et al. . Briefly, after plant material fixation in 1% (v/v) formaldehyde, tissues were homogenized, nuclei isolated, and lysed. Cross-linked chromatin was sonicated using a water bath Bioruptor UCD-200 (Diagenode, Liège, Belgium) (30 s on/off pulses, at high intensity for 60 min). Protein/DNA complexes were immunoprecipitated with antibodies, overnight at 4 °C with gentle shaking, and incubated for 1 h at 4 °C with 50 μL of Dynabeads Protein A (Invitrogen, Ref. 100-02D). Immunoprecipitated DNA was then recovered using the IPure kit (Diagenode, Liège, Belgium) and analyzed by RT-qPCR. An aliquot of untreated sonicated chromatin was processed in parallel and used as the total input DNA control.
For ChIP-qPCR experiments, fold enrichment of targets in ChIPed DNA relative to input was calculated from an average of three replicate qPCR reactions. The sequences of primers can be found in Additional file 6: Table SV. Positions of the amplified regions on the different loci are indicated in Additional file 1: Figure S30.
After immunoprecipitation of the chromatin, ChIP-seq libraries were generated and sequenced. Alignment was performed using Bowtie  v0.12.7 on Arabidopsis thaliana genome TAIR10. Default Bowtie parameters were used except for: -best -strata (used to get the best mapping position with the minimum of mismatches). Peak calling was performed with MACS (http://liulab.dfci.harvard.edu/MACS/) . Gene annotation and peak distribution relative to annotated Arabidopsis transcription start site was performed with GPAT .
Global clustering of the H3K9ac ChIP-seq data and quantitative comparisons were performed using the seqMINER program (http://bips.u-strasbg.fr/seqminer/) . As reference coordinates, we used the MACS determined peaks for H3K9ac. Tag densities from each ChIP-seq dataset were collected in a window of 1 kb around the reference peak. The collected values were subjected to k-means clustering coupled to linear-based normalization. The normalization procedure reduces bias in the clustering due to inherent differences between ChIP-seq experiments.
RNA-seq and ChIP-seq data from this article were deposited at GEO Project GSE99936.
Kouzarides T. Chromatin modifications and their function. Cell. 2007;128:693–705.
Mercer TR, Mattick JS. Understanding the regulatory and transcriptional complexity of the genome through structure. Genome Res. 2013;23(7):1081–108.
Saze H, Tsugane K, Kanno T, Nishimura T. DNA methylation in plants: relationship to small RNAs and histone modifications, and functions in transposon inactivation. Plant Cell Physiol. 2012;53:766–84.
Bannister AJ, Kouzarides T. Regulation of chromatin by histone modifications. Cell Res. 2011;21(3):381–95.
Bender J. RNA-directed DNA, methylation: getting a grip on mechanism. Curr Biol. 2012;22(10):R400–1.
Keller C, Kulasegaran-Shylini R, Shimada Y, Hotz HR, Buhler M. Noncoding RNAs prevent spreading of a repressive histone mark. Nat Struct Mol Biol. 2013;20(8):994–1000.
Asai T, Tena G, Plotnikova J, Willmann MR, Chiu WL, Gomez-Gomez L, et al. MAP kinase signalling cascade in Arabidopsis innate immunity. Nature. 2002;415:977–83.
Gao M, Liu J, Bi D, Zhang Z, Cheng F, Chen S, et al. MEKK1, MKK1/MKK2 and MPK4 function together in a mitogen-activated protein kinase cascade to regulate innate immunity in plants. Cell Res. 2008;18(12):1190–8.
Bethke G, Pecher P, Eschen-Lippold L, Tsuda K, Katagiri F, Glazebrook J, et al. Activation of the Arabidopsis thaliana mitogen-activated protein kinase MPK11 by the flagellin-derived elicitor peptide, flg22. Mol Plant Microbe Interact. 2012;25:471–80.
Djamei A, Pitzschke A, Nakagami H, Rajh I, Hirt H. Trojan horse strategy in Agrobacterium transformation: abusing MAPK defense signaling. Science. 2007;318:453–6.
Andreasson E, Jenkins T, Brodersen P, Thorgrimsen S, Petersen NH, Zhu S, et al. The MAP kinase substrate MKS1 is a regulator of plant defense responses. EMBO J. 2005;24:2579–89.
Bethke G, Unthan T, Uhrig JF, Poschl Y, Gust AA, Scheel D, et al. Flg22 regulates the release of an ethylene response factor substrate from MAP kinase 6 in Arabidopsis thaliana via ethylene signaling. Proc Natl Acad Sci U S A. 2009;106:8067–72.
de la Fuente van Bentem S, Anrather D, Dohnal I, Roitinger E, Csaszar E, Joore J, et al. Site-specific phosphorylation profiling of Arabidopsis proteins by mass spectrometry and peptide chip analysis. J Proteome Res. 2008;7(6):2458–70.
Zhou C, Labbe H, Sridha S, Wang L, Tian L, Latoszek-Green M, et al. Expression and function of HD2-type histone deacetylases in Arabidopsis development. Plant J. 2004;38:715–24.
Pendle AF, Clark GP, Boon R, Lewandowska D, Lam YW, Andersen J, et al. Proteomic analysis of the Arabidopsis nucleolus suggests novel nucleolar functions. Mol Biol Cell. 2005;16(1):260–9.
Frei dit Frey N, Garcia AV, Bigeard J, Zaag R, Bueso E, Garmier M, et al. Functional analysis of Arabidopsis immune-related MAPKs uncovers a role for MPK3 as negative regulator of inducible defences. Genome Biol. 2014;15(6):R87.
Ye T, Krebs AR, Choukrallah MA, Keime C, Plewniak F, Davidson I, et al. seqMINER: an integrated ChIP-seq data interpretation platform. Nucleic Acids Res. 2011;39(6):e35.
Pitzschke A, Djamei A, Teige M, Hirt H. VIP1 response elements mediate mitogen-activated protein kinase 3-induced stress gene expression. Proc Natl Acad Sci U S A. 2009;106(43):18414–9.
Wang C, Gao F, Wu J, Dai J, Wei C, Li Y. Arabidopsis putative deacetylase AtSRT2 regulates basal defense by suppressing PAD4, EDS5 and SID2 expression. Plant Cell Physiol. 2010;51(8):1291–9.
Tian L, Fong MP, Wang JJ, Wei NE, Jiang H, Doerge RW, et al. Reversible histone acetylation and deacetylation mediate genome-wide, promoter-dependent and locus-specific changes in gene expression during plant development. Genetics. 2005;169(1):337–45.
Choi SM, Song HR, Han SK, Han M, Kim CY, Park J, et al. HDA19 is required for the repression of salicylic acid biosynthesis and salicylic acid-mediated defense responses in Arabidopsis. Plant J. 2012;71(1):135–46.
Kim KC, Lai Z, Fan B, Chen Z. Arabidopsis WRKY38 and WRKY62 transcription factors interact with histone deacetylase 19 in basal defense. Plant Cell. 2008;20(9):2357–71.
Ding B, Bellizzi Mdel R, Ning Y, Meyers BC, Wang GL. HDT701, a histone H4 deacetylase, negatively regulates plant innate immunity by modulating histone H4 acetylation of defense-related genes in rice. Plant Cell. 2012;24(9):3783–94.
Chen YJ, Wang YN, Chang WC. ERK2-mediated C-terminal serine phosphorylation of p300 is vital to the regulation of epidermal growth factor-induced keratin 16 gene expression. J Biol Chem. 2007;282:27215–28.
Oya H, Yokoyama A, Yamaoka I, Fujiki R, Yonezawa M, Youn MY, et al. Phosphorylation of Williams syndrome transcription factor by MAPK induces a switching between two distinct chromatin remodeling complexes. J Biol Chem. 2009;284(47):32472–82.
Dangl M, Brosch G, Haas H, Loidl P, Lusser A. Comparative analysis of HD2 type histone deacetylases in higher plants. Planta. 2001;213:280–5.
Wu K, Tian L, Zhou C, Brown D, Miki B. Repression of gene expression by Arabidopsis HD2 histone deacetylases. Plant J. 2003;34:241–7.
Bjerling P, Silverstein RA, Thon G, Caudy A, Grewal S, Ekwall K. Functional divergence between histone deacetylases in fission yeast by distinct cellular localization and in vivo specificity. Mol Cell Biol. 2002;22:2170–81.
Segre CV, Chiocca S. Regulating the regulators: the post-translational code of class I HDAC1 and HDAC2. J Biomed Biotechnol. 2011;2011:690848.
Dequiedt F, Martin M, Von Blume J, Vertommen D, Lecomte E, Mari N, et al. New role for hPar-1 kinases EMK and C-TAK1 in regulating localization and activity of class IIa histone deacetylases. Mol Cell Biol. 2006;26(19):7086–102.
Sun JM, Chen HY, Moniwa M, Litchfield DW, Seto E, Davie JR. The transcriptional repressor Sp3 is associated with CK2-phosphorylated histone deacetylase 2. J Biol Chem. 2002;277(39):35783–6.
Sun JM, Chen HY, Davie JR. Differential distribution of unmodified and phosphorylated histone deacetylase 2 in chromatin. J Biol Chem. 2007;282(45):33227–36.
Wang Y, Curry HM, Zwilling BS, Lafuse WP. Mycobacteria inhibition of IFN-gamma induced HLA-DR gene expression by up-regulating histone deacetylation at the promoter region in human THP-1 monocytic cells. J Immunol. 2005;174(9):5687–94.
Basu S, Pathak S, Pathak SK, Bhattacharyya A, Banerjee A, Kundu M, et al. Mycobacterium avium-induced matrix metalloproteinase-9 expression occurs in a cyclooxygenase-2-dependent manner and involves phosphorylation- and acetylation-dependent chromatin modification. Cell Microbiol. 2007;9(12):2804–16.
Garcia-Garcia JC, Barat NC, Trembley SJ, Dumler JS. Epigenetic silencing of host cell defense genes enhances intracellular survival of the rickettsial pathogen Anaplasma phagocytophilum. PLoS Pathog. 2009;5(6):e1000488.
Eskandarian HA, Impens F, Nahori MA, Soubigou G, Coppee JY, Cossart P, et al. A role for SIRT2-dependent histone H3K18 deacetylation in bacterial infection. Science. 2013;341(6145):1238858.
Kim JM, To TK, Nishioka T, Seki M. Chromatin regulation functions in plant abiotic stress responses. Plant Cell Environ. 2010;33:604–11.
Bourque S, Dutartre A, Hammoudi V, Blanc S, Dahan J, Jeandroz S, et al. Type-2 histone deacetylases as new regulators of elicitor-induced cell death in plants. New Phytol. 2011;192:127–39.
Waterborg JH. Dynamics of histone acetylation in Saccharomyces cerevisiae. Biochemistry. 2001;40:2599–605.
Scott I. Regulation of cellular homoeostasis by reversible lysine acetylation. Essays Biochem. 2012;52:13–22.
Sessions A, Burke E, Presting G, Aux G, McElver J, Patton D, et al. A high-throughput Arabidopsis reverse genetics system. Plant Cell. 2002;14:2985–94.
Wang H, Ngwenyama N, Liu Y, Walker JC, Zhang S. Stomatal development and patterning are regulated by environmentally responsive mitogen-activated protein kinases in Arabidopsis. Plant Cell. 2007;19:63–73.
Xiang C, Han P, Lutziger I, Wang K, Oliver DJ. A mini binary vector series for plant transformation. Plant Mol Biol. 1999;40:711–7.
Clough SJ, Bent AF. Floral dip: a simplified method for Agrobacterium-mediated transformation of Arabidopsis thaliana. Plant J. 1998;16:735–43.
Karimi M, Inze D, Depicker A. GATEWAY vectors for Agrobacterium-mediated plant transformation. Trends Plant Sci. 2002;7(5):193–5.
Stuitje AR, Verbree EC, van der Linden KH, Mietkiewska EM, Nap JP, Kneppers TJ. Seed-expressed fluorescent proteins as versatile tools for easy (co)transformation and high-throughput functional genomics in Arabidopsis. Plant Biotechnol J. 2003;1:301–9.
Bechtold N, Pelletier G. In planta Agrobacterium-mediated transformation of adult Arabidopsis thaliana plants by vacuum infiltration. Methods Mol Biol. 1998;82:259–66.
Kemmerling B, Schwedt A, Rodriguez P, Mazzotta S, Frank M, Qamar SA, et al. The BRI1-associated kinase 1, BAK1, has a brassinolide-independent role in plant cell-death control. Curr Biol. 2007;17(13):1116–22.
Cardinale F, Jonak C, Ligterink W, Niehaus K, Boller T, Hirt H. Differential activation of four specific MAPK pathways by distinct elicitors. J Biol Chem. 2000;275(47):36734–40.
Sauer M, Paciorek T, Benkova E, Friml J. Immunocytochemical techniques for whole-mount in situ protein localization in plants. Nat Protoc. 2006;1(1):98–103.
Azimzadeh J, Nacry P, Christodoulidou A, Drevensek S, Camilleri C, Amiour N, Parcy F, Pastuglia M, Bouchez D. Arabidopsis TONNEAU1 Proteins Are Essential for Preprophase Band Formation and Interact with Centrin. THE PLANT CELL ONLINE. 2008;20(8):2146–2159.
Colcombet J, Lopez-Obando M, Heurtevin L, Bernard C, Martin K, Berthomé R, Lurin C. Systematic study of subcellular localization of Arabidopsis PPR proteins confirms a massive targeting to organelles. RNA Biology. 2014;10(9):1557–1575.
Nakagami H, Soukupová H, Schikora A, Zárský V, Hirt H. A Mitogenactivated Protein Kinase Kinase Kinase Mediates Reactive Oxygen Species Homeostasis in . Journal of Biological Chemistry. 2006;281(50):38697–38704.
Cardinale F, Meskiene I, Ouaked F, Hirt H. Convergence and divergence of stress-induced mitogen-activated protein kinase signaling pathways at the level of two distinct mitogen-activated protein kinase kinases. Plant Cell. 2002;14:703–11.
Du Z, Zhou X, Ling Y, Zhang Z, Su Z. agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 2010;8:W64–70.
Gendrel AV, Lippman Z, Martienssen R, Colot V. Profiling histone modification patterns in plants using genomic tiling microarrays. Nat Methods. 2005;2:213–8.
Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25.
Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008;9:R137.
Krebs A, Frontini M, Tora L, GPAT: Retrieval of genomic annotation from large genomic position datasets. BMC Bioinformatics. 2008;9(1):533.
Hruz T, Laule O, Szabo G, Wessendorp F, Bleuler S, Oertie L, et al. Genevestigator v3: a reference expression database for the meta-analysis of transcriptomes. Adv Bioinf. 2008;2008:420747.
Funding for the project was provided to HH and MB by a grant of the French National Agency of Research (LabEx Saclay Plant Sciences-SPS ANR-10-LABX-0040-SPS) as well as by grant CRG4 of the King Abdullah University of Science and Technology. Moreover, OS was supported by the Czech Science Foundation GA ČR grant no. GA16-22044S. No ethics approval was required for this project.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplemental figures. Figure S1–S5. HD2B phosphorylation sites. Figure S6. MPK3, MPK4, MPK6, and MPK7 antibody specificity. Figure S7. MPK3, MPK4, and MPK6 substrate specificity. Figure S8. Subcellular localization of MPK3. Figure S9. BiFC analysis of HD2B interaction with MPK3, 4, and 6 in epidermal cells of Agrobacterium-infiltrated Nicotiana benthamiana. Figure S10. Co-imunoprecipitation of HD2B with MPK3. Figure S11. GO analysis of transcription factors up-regulated in hd2b mutant. Figure S12. GO analysis of downregulated genes in hd2b mutant. Figure S13. Validation of transcriptomic data by RT-qPCR. Figure S14. Susceptibility of hd2b mutant and HD2B overexpressing lines. Figure S15. Flg22-induced activation of MPK3, MPK4, and MPK6 is similar in Col-0 and hd2b mutant. Figure S16. Flg22 growth inhibition assay of hd2b mutant. Figure S17. Flagellin-induced relocalization of HD2B. Figure S18. Nuleolar HD2B and nuclear HD2B-ED localization. Figure S19. Nuclear HD2B localization. Figure S20. Gene expression analyses of hd2b and mpk3 over-expressed genes in the HD2B-AA line. Figure S21. Flagellin-independent expression of HD2B. Figure S22. Specific HD2B targets in mock conditions are involved in defense response. Figure S23. Correlation between HD2B binding and gene expression in hd2b mutants in mock and flg22 conditions. Figure S24. Characterization of H3K9 acetylated regions in mock- and flg22-treated seedlings. Figure S25. Comparison between HD2B and H3K9ac peak positions relative to TSS. Figure S26. GO analysis of H3K9-hyper and H3K9-hypo acetylated genes after flg22 treatment. Figure S27. H3K9 acetylation levels are directly correlated with expression levels after flg22 treatment. Figure S28. Gene expression analyses of hd2b and mpk3 flg22 mis-regulated genes in the HD2B-AA line. Figure S29. Characterization of hd2b mutant line. Figure S30. Positions of the regions analyzed by ChIP-qPCR experiments on the different loci. (PDF 1569 kb)
List of HD2B target genes in mock and flg22 conditions. (XLSX 167 kb)
List of H3K9 hyper- and hypo-acetylated genes after flg22 treatment. (XLSX 60 kb)
List of deregulated genes in hd2b and mpk3 mutants. (XLSX 51 kb)
List of deregulated genes after flg22 treatment in Wt. (XLSX 52 kb)
List of primers used in this study. (XLSX 13 kb)