An evolutionarily conserved histone modification H3K37ac activates gene transcription in response to salt stress in rice

Abstract

Histone acetylation is a well-known chromatin mark that is associated with gene transcriptional activation (Y. Chen et al., 2024). In plants, histone acetylation modification mainly occurs at the N-terminal lysine (Lys) residues of histone H3 (K9, K14, K18, K23, K27, K36) and H4 (K5, K8, K12, K16, K20) (Waterborg, 2011). The acetylation status of specific Lys residues on histone proteins is reversibly modulated by histone acetyltransferases (HATs) and histone deacetylases (HDACs). These enzymes play essential roles in plant growth, development, and responses to abiotic stresses by regulating gene transcription. Many HDACs have been reported to be involved in modulating plants' tolerance to abiotic stresses (Ueda & Seki, 2020; Cui et al., 2023). For instance, HDAC, HDA15, interacts with MYB96 to repress negative regulators of early abscisic acid signaling and improve drought resistance in Arabidopsis (Lee & Seo, 2019). Similarly, the rice HDA706 was found to enhance salt tolerance by deacetylating the H4K5/K8ac at the promoter region of PP2C49 (Liu et al., 2023). By contrast, the Arabidopsis HDA9, responsible for H3K9ac deacetylation, was reported to negatively regulate salt and drought stress responses (Zheng et al., 2016). These studies highlight the pivotal roles of histone acetylation in modulating the transcriptional programs underlying plant abiotic stresses tolerance. However, our understanding of the full repertoire of histone Lys acetylation sites and their specific functions remains incomplete. Herein, we identified an unexplored histone Lys acetylation in rice, Lys-37 of histone H3 (H3K37ac); we found this acetylation modification is highly conserved across plant species. Notably, the level of H3K37 acetylation in rice is strongly induced by salinity stress. Furthermore, we have identified a HDAC, HDA705, that is responsible for removing the acetylation mark from H3K37, and mutation of HDA705 significantly improved the salt tolerance of rice plants. The ‘Zhonghua 11’ (ZH11) rice variety (Oryza sativa ssp. japonica) was used in this study. The hda704, hda705, hda706, hda710, hda712, and srt2 mutants (ZH11 background) used in this study were previously characterized (Xu et al., 2021, 2025). For salt stress treatment, both wild-type (WT) and mutant plants (hda705) were first hydroponically grown in liquid Yoshida's solution (Z. Chen et al., 2024). Fourteen-day-old seedlings, cultivated under 25°C with a 14 h : 10 h, light : dark cycle, were then subjected to salt stress by adding 180 mM NaCl to the nutrient solution for 3 d. Control seedlings were grown under the same temperature and light conditions but without salt stress. Following the stress treatments, plants were returned to normal conditions for 7 d, after which the survival rate was determined by calculating the percentage of viable plants (live plants/total plants × 100%). The survival rates were assessed from three groups, with each group containing 40 plants (40 replicates). For osmotic stress, 14-d-old seedlings were exposed to 20% PEG6000 in the nutrient solution for 6 h, after which leaf samples were collected for histone extraction. Heat stress was applied by incubating 14-d-old seedlings at 42°C for 6 h in a growth chamber, followed by leaf collection for histone extraction. Cold stress treatments involved incubating 14-d-old seedlings at 4°C for 6 h, with leaf samples harvested afterward for histone extraction. The histone H3K37 acetylation (H3K37ac) was identified from a global rice acetyl-proteome dataset generated by liquid chromatography tandem mass spectrometry (LC-MS/MS; Ma et al., 2023). Briefly, proteins extracted from rice panicles were digested with trypsin after reduction and alkylation. The resulting peptides were incubated with anti-acetyl-lysine beads (cat no.: PTM-104; PTM Bio, CN) for immunoaffinity enrichment. The bound acetylated peptides were eluted, vacuum-dried, and desalted using C18 ZipTips (Millipore, Billerica, MA, USA) before LC-MS/MS analysis. For LC-MS/MS analysis, the enriched acetyl-peptides were analyzed by online LC-MS/MS using an EASY-nLC 1200 UPLC system (Thermo Fisher Scientific, Waltham, MA, USA) coupled to a Q Exactive™ Plus mass spectrometer (Thermo Fisher Scientific). Peptides were separated on a reversed-phase column (15 cm and 75 μm i.d.) with a 60-min linear gradient from 6% to 35% solvent B (0.1% formic acid in 98% acetonitrile), followed by a step to 80% B, at a flow rate of 400 nL min−1. MS analysis was performed with a nano-electrospray ionization source at 2.0 kV. Full-scan MS spectra (m/z 350–1800) were acquired at a resolution of 70 000, followed by data-dependent MS/MS scans of the top 20 precursors (NCE: 28; resolution: 17500). Dynamic exclusion was set to 15 s. The raw MS/MS data were processed with MaxQuant (v.1.6.6.0) against the Oryza sativa protein database (http://rice.hzau.edu.cn/rice_rs2/). Search parameters included: trypsin/P as the cleavage enzyme, allowing up to two miscleavages; peptides with a minimum length of seven amino acids and a maximum of five modifications; precursor mass tolerance of 5 ppm (main search) and fragment mass tolerance of 0.04 Da; fixed modification of carbamidomethylation (C); variable modifications of oxidation (M), acetylation (K), and N-terminal acetylation. The false discovery rate (FDR) threshold for proteins, peptides, and modified sites was set to 40. Total histone-enriched fractions were extracted from Arabidopsis, tobacco, maize, and rice seedlings using an EpiQuik Total Histone Extraction Kit (OP-0006-100; Epigentek, Farmingdale, NY, USA). Before immunoblotting, half of the volume of 2 × SDS loading buffer (120 mM Tris–HCl, pH 6.8, 20% glycerol, 4% SDS, 0.04% bromophenol blue, and 10% β-mercaptoethanol) was added to the supernatant of extracted. The mixture was thoroughly combined and heated at 95°C for 10 min. The histone proteins were then subjected to immunoblot analysis using anti-H3 (1 : 1000, ab1791; Abcam, Cambridge, UK) and anti-H3K37ac (1 : 1000, PTM-128; PTM, CN) antibodies. The specificity of the H3K37ac antibody was examined by dot blot and peptide competition assays. For the dot blot assay, 1, 5, or 10 ng of each synthetic peptide (H3K37ac, H3K36ac, and H3K27ac) was spotted onto a nitrocellulose membrane and air-dried. The membrane was blocked with 5% nonfat milk in Tris Buffered Saline with Tween-20 (TBST) for 1 h and then incubated with the H3K37ac antibody (1 : 1000 dilution) overnight at 4°C. Signals were detected using an horseradish peroxide (HRP)-conjugated secondary antibody and chemiluminescence. For the peptide competition assay, the H3K37ac antibody (1 μg) was preincubated with synthetic H3K37ac peptide (0, 1, 2, 5, or 10 μg) in 10 ml of TBST at room temperature for 1 h before being used to probe histones separated by 15% Sodium dodecyl sulfate-Polyacrylamide gel electrophoresis (SDS-PAGE) and transferred onto polyvinylidene fluoride membranes. The membranes were subsequently incubated with HRP-conjugated secondary antibodies, and signals were visualized using ECL detection reagents. The full-length coding sequences of SRT2 and HDA705 were cloned into the PGEX-6P-1 expression vector for protein production. These plasmids were used to express the respective proteins in Escherichia coli BL21 (DE3) cells. The GST-tagged proteins GST, SRT2-GST, and HDA705-GST were then purified using GST beads (17–5132-01; GE Healthcare, Chicago, IL, USA). Nuclear proteins were extracted from WT rice leaves and used for the in vitro de-Kac activity assay. Purified GST-tagged HDACs were incubated with the extracted nuclear proteins in a reaction buffer containing 50 mmol l−1 Tris–HCl (pH 8.5), 2.7 mmol l−1 KCl, 1 mmol l−1 MgCl2, 137 mmol l−1 NaCl, and 1 mmol l−1 dithiothreitol at 37°C for 4 h. The reaction products were analyzed by immunoblotting with anti-H3K37ac antibodies (1 : 1000, PTM-128; PTM, CN). For the HDAC activity assay, constructs expressing HDA705-3 × FLAG, SRT2-3 × FLAG, H3-6 × His, and H3K37A-6 × His were introduced into Nicotiana benthamiana leaf epidermal cells. After 48 h of infiltration, the transformed leaves were harvested and frozen in liquid nitrogen for tissue disruption. Total proteins were extracted using a lysis buffer (10 mM Tris–HCl, pH 7.5; 150 mM NaCl; 0.5 mM EDTA; 0.5% Nonidet™ P40 Substitute), purified with Ni-NTA Magnetic Agarose Beads (36113; QIAGEN, Hilden, Germany), and analyzed by western blotting. One microgram of total RNA was reverse-transcribed in a 20 μl reaction volume using DNase and reverse transcriptase (R233-01; Vazyme, Nanjing, China) according to the manufacturer's instructions to generate cDNA. RT-qPCR was performed on an ABI PRISM 7500 device. Three biological replicates were conducted for each gene. The rice ACTIN1 gene was used as the internal control. Chromatin immunoprecipitation (ChIP) protocol followed a previously established method (Xu et al., 2023). Two grams of 14-d-old rice seedlings were cross-linked by 1% (v/v) formaldehyde and used for chromatin extraction. Chromatin was fragmented to c. 200 bp by sonication and then incubated with antibody-coated beads (10001D; Invitrogen/Life Technologies, Carlsbad, CA, USA) overnight. After extensive washing, immunoprecipitated chromatin was de-cross-linked and retrieved, and un-immunoprecipitated chromatin was used as input. The precipitated and input DNA samples were analyzed by real-time PCR with gene-specific primers (Supporting Information Table S1). The raw reads of RNA-seq were processed to remove low-quality sequence using FastP (v.0.232). Clean reads were then aligned to the rice genome sequence (MSU v.7.0) using the Subjunc aligner with default settings. Differential gene expression analysis was carried out using DESeq2 (v.1.36.0), with genes exhibiting a fold change (FC) > 2 and a P-adjusted value of < were under salt stress (Fig. Table S3). analysis of genes verified the CUT&Tag results (Fig. GO analysis revealed that the genes with H3K37ac were mainly enriched in regulation of gene response to stimulus, and response (Fig. To H3K37ac modification is involved in modulating salt stress we analyzed the of rice seedlings under salt stress same as for with three biological replicates for both and In we identified genes < Table and genes < Table By CUT&Tag we found that genes H3K37ac enrichment at their salinity genes (Fig. analysis of and CUT&Tag data revealed that genes both transcription levels and H3K37ac modification levels under salt stress (Fig. Table GO enrichment analysis revealed that these genes mainly involved in response to stimulus, response to stress (Fig. genes previously such as and et al., et al., To we performed real-time and to the transcription and H3K37ac modification levels of five and with the analysis, the expression levels and H3K37ac modification levels of the five genes were significantly under salt stress (Fig. Collectively, these results that H3K37ac modification is involved in modulating salt stress response in The correlation between H3K37ac levels and the expression of genes that this novel chromatin mark an role in the transcriptional regulation of the stress To the enzymes responsible for H3K37ac we analyzed the modification levels in different rice HDACs assays showed that H3K37ac level was in mutant and was in HDAC mutant srt2 (Fig. To their we histone H3 with SRT2 or HDA705 in tobacco purified the and performed western with the H3K37ac antibody. of the enzymes significantly H3K37ac their to histone H3 at the K37 (Fig. Furthermore, in vitro HDAC activity assays that HDA705 activity (Fig. To investigate whether HDA705 salt stress responses H3K37ac regulation in rice, WT plants and two mutant and were subjected to salt stress The results that the mutants with a higher survival rate the WT after from salt stress mM NaCl; and revealed that the expression levels and H3K37ac levels of genes in mutants were significantly to WT (Fig. that the of H3K37ac on the genes their expression and salt tolerance in mutant plants. these data that HDA705 is the for H3K37ac and that this a role in modulating rice salt tolerance. In this study identified H3K37ac as a previously histone acetylation mark that is associated with transcriptional activation and to salt stress in We that the HDAC HDA705 is responsible for the of and that the mutant tolerance to These into the potential underlying plant responses to salt stress. our data indicated that salinity stress enrichment of H3K37ac in the rice genome, remains whether this histone modification the of such as or of be to whether this modification gene Furthermore, we established that the HDAC, HDA705, H3K37ac and that HDA705 mutants salt a role for HDA705 in salt tolerance H3K37ac However, HDA705 been reported to other histone lysine sites and various histone modifications (Xu et al., 2021, the that these to salt as be essential to determine whether HDA705 H3K37ac other histone modifications and whether its activity or specificity is by specific or these be for understanding the of regulation in plant Collectively, our study a for understanding the regulation of in plants and a for the roles of histone acetylation in stress tolerance. study was by for the and the for of the of and the of the of the for from and the and of the of and the and performed most of the and performed the data and to the and and the have and the and to this The CUT&Tag data of H3K37ac in this study were to the database at the under the no.: The RNA-seq data used in this study were from with the no.: acid sequence alignment of histone protein H3 from different species. of H3K37ac antibodies by dot competition with H3K37ac of histone H3K37ac activity of HDA705 and SRT2 in tobacco cells. of H3K37ac data in The of H3K37ac modifications in rice GO pathway enrichment analysis of the H3K37ac-marked genes in wild-type rice plants. of the expression levels of the H3K37ac modified genes in of and regions in analysis of H3K37ac modification levels of genes in and plants. showing H3K37ac enrichment the of genes induced by salt and genes in rice plants. Table used in this study. Table H3K37ac modification detected in Table Dynamic of H3K37ac in salt plants. Table Differentially expressed genes in salt plants. Table with higher expression and H3K37ac in salt plants. is responsible for the or of Information by the be to the The is responsible for the or of by the be to the for the The remains with to in and in