- Open Access
Dynamics and functional interplay of histone lysine butyrylation, crotonylation, and acetylation in rice under starvation and submergence
- Yue Lu†1,
- Qiutao Xu†1,
- Yuan Liu1,
- Yue Yu1,
- Zhong-Yi Cheng3,
- Yu Zhao1 and
- Dao-Xiu Zhou1, 2Email authorView ORCID ID profile
© The Author(s). 2018
- Received: 25 March 2018
- Accepted: 10 September 2018
- Published: 25 September 2018
Histone lysine acylations by short-chain fatty acids are distinct from the widely studied histone lysine acetylation in chromatin, although both modifications are regulated by primary metabolism in mammalian cells. It remains unknown whether and how histone acylation and acetylation interact to regulate gene expression in plants that have distinct regulatory pathways of primary metabolism.
We identify 4 lysine butyrylation (Kbu) sites (H3K14, H4K12, H2BK42, and H2BK134) and 45 crotonylation (Kcr) sites on rice histones by mass spectrometry. Comparative analysis of genome-wide Kbu and Kcr and H3K9ac in combination with RNA sequencing reveals 25,306 genes marked by Kbu and Kcr in rice and more than 95% of H3K9ac-marked genes are marked by both. Kbu and Kcr are enriched at the 5′ region of expressed genes. In rice under starvation and submergence, Kbu and Kcr appear to be less dynamic and display changes in different sets of genes compared to H3K9ac. Furthermore, Kbu seems to preferentially poise gene activation by external stresses, rather than internal circadian rhythm which has been shown to be tightly associated with H3K9ac. In addition, we show that rice sirtuin histone deacetylase (SRT2) is involved in the removal of Kcr.
Kbu, Kcr, and H3K9ac redundantly mark a large number of active genes but display different responses to external and internal signals. Thus, the proportion of rice histone lysine acetylation and acylation is dynamically regulated by environmental and metabolic cues, which may represent an epigenetic mechanism to fine-tune gene expression for plant adaptation.
Post-translational modifications of histones affect chromatin state and gene expression. For example, histone lysine acetylation on the ε-amine group not only neutralizes the positive charge of the amine group, enhances the hydrophobicity, and increases the size of the lysine side chain, but also provides platforms for binding by proteins involved in chromatin and gene regulations. Histone acetylation/deacetylation is an essential process in the epigenetic regulation of diverse biological processes, including environmental stress responses in plants [1–3]. The level of histone acetylation is determined by the activity of both histone acetyl-transferases (HAT) and deacetylases (HDAC). Histone acetylation can be regulated by primary metabolism, as key metabolites such as acetyl-CoA and NAD+ are involved in protein acetylation and deacetylation processes, and their cellular levels may regulate the activity of HATs and HDACs [4, 5]. Plant metabolism is largely regulated to adapt to the fluctuating environment. Thus, the collaborative program between epigenetic dynamics and metabolism is further interconnected with environmental cues, which may optimize plant adaptive responses and growth .
Recently, it is shown that histones can be acylated with short-chain fatty acids, including propionylation, butyrylation, crotonylation, 2-hydroxyisobutyrylation, succinylation, malonylation, glutarylation, and β-hydroxybutyrylation in animal cells [6–8]. These modifications are similar to lysine acetylation but are distinct in hydrocarbon chain length and hydrophobicity or charge. Studies in animal cells suggest that histone lysine acylations affect gene expression and may be functionally different from histone lysine acetylation . Histone lysine acylations are shown to be regulated by metabolism. For instance, intracellular crotonyl-CoA stimulates transcription through histone crotonylation . In addition, β-hydroxybutyrylation is induced significantly during prolonged fasting in mouse liver and is associated with gene upregulation in starvation-responsive metabolism pathways .
It has been reported that histone acetyltransferases (HATs) and deacetylases (HDACs) could be involved in histone acylation dynamics . For instance, mammalian CBP/p300, a histone acetyltransferase, catalyzes both histone acetyl and non-acetyl acylation [6, 10, 13, 14]. All seven human sirtuins (SIRT1-7) that were originally annotated as histone deacetylases based on sequence homology to the yeast sir2 histone deacetylase, have been described as having diverse deacylase activities, albeit weaker than their deacetylase activities [9, 13, 15]. Other studies have shown that Zn-finger domain-containing HDACs are also involved in decrotonylation of histone lysine residues .
Histone acylation has not been studied yet in plants. Being photosynthetic organisms, plants are autotrophic and there are significant differences in primary and energy metabolic and signaling pathways. For instance, glycolysis is induced under both biotic and abiotic stresses, and tricarboxylic acid (TCA) cycle is particularly high during light to dark transition [16, 17]. It remains unknown how histone acylation and acetylation interact to regulate gene expression in plants that have distinct regulatory pathways of primary metabolism.
In this work, we first analyzed rice histone butyrylation (Kbu) and crotonylation (Kcr) sites by mass spectrometry, then compared genome-wide distribution of Kbu and Kcr and lysine acetylation (H3K9ac) together with transcriptomic analysis in rice. To study whether altered metabolic activity affected histone acylations and acetylation, we treated rice plants by continuous darkness (starvation) and submergence and investigated genome-wide changes of Kbu, Kcr, and H3K9ac and gene expression. Finally, we tested histone deacylation activities of several rice HDACs. Our results indicate that Kbu and Kcr are robust histone modification marks that regulate gene expression in rice. Rice histone acylations and acetylation are dynamically regulated by environmental cues and are likely to be removed by different mechanisms.
Identification of rice histone lysine butyrylation (Kbu) and crotonylation (Kcr) sites
Comparative analysis of rice genome-wide histone Kbu, Kcr, and H3K9ac
To study potential roles of histone acylations in gene expression, we performed RNA-seq analysis of the same plant materials (Additional file 1: Table S2). To investigate functional relationship between Kbu, Kcr, and H3K9ac in gene expression, we analyzed expression levels of genes modified by different combinations of the marks (Fig. 3d). Genes marked by H3K9ac showed higher expression than those marked only by Kbu or Kcr. Genes marked by a combination of two marks showed higher expression than the mono-marked genes. Genes with all three marks displayed the highest expression levels. In genes with the three marks, those with higher portion of H3K9ac relative to Kbu or Kcr displayed higher expression levels (Fig. 3e). However, the total acylation level of genes with the three marks was not higher than those with only two marks (Fig. 3f). The analysis suggests that H3K9ac is predominant in gene activity and that Kbu and Kcr may contribute to H3K9ac-regulated active chromatin state. Genes marked by Kbu or Kcr only (tags per million mapped tags per kb, TPM > 50) showed low expression levels (FPKM < 0.1) and were enriched in stress and abiotic stimulus responses. Many of them are highly induced under stress (Additional file 1: Figure S3), hypothesizing that Kbu and Kcr may play a role in poising genes for stimulus-induced activation.
Kbu, Kcr, and H3K9ac changes in rice under starvation and submergence
The substrates for histone acetylation/acylation, such as acetyl-CoA, butyryl-CoA, and crotonyl-CoA, are important cellular primary metabolites. Stress conditions such as prolonged darkness (starvation) or submergence (switching to anaerobic respiration) affect greatly plant primary metabolism. To study whether starvation and submergence altered Kbu/Kcr/H3K9ac, we treated rice seedlings in darkness for 48 h or under submergence for 6 h. qRT-PCR analysis of starvation and submergence-induced marker genes confirmed the effectiveness of the treatments (Additional file 1: Figure S4A). The materials were harvested at 6 pm and used for ChIP-seq with anti-Kbu, anti-Kcr, and anti-H3K9ac and RNA-seq.
In parallel, RNA-seq analysis revealed that starvation and submergence resulted in differential expression (FDR < 0.05, fold change > 4) of 1967 (1045 upregulated, 922 downregulated) and 1305 (487 upregulated, 818 downregulated) genes, respectively (Additional file 1: Table S2; Figure S4B). GO analysis indicated that the differentially expressed genes induced by starvation are enriched in energy metabolic processes and in stress and stimulus response pathways (Additional file 1: Table S3). Submergence-affected genes were enriched in nucleic acid metabolism and DNA-binding processes in addition to stress and stimulus response pathways (Additional file 1: Table S4). Some of the enriched gene categories were also enriched for H3K9ac, Kbu, and Kcr changes (Additional file 1: Table S3, Table S4). Remarkably, transcription factor genes were found to be enriched among those with increased or decreased H3K9ac under starvation and submergence and with decreased Kbu and Kcr under submergence (Additional file 1: Table S3, Table S4).
Kbu may poise external stress-induced gene expression
To check whether Kbu and Kcr played a role to poise genes for stimulus-induced expression, we analyzed upregulated genes under starvation and submergence and found that the genes marked by Kbu alone or by Kbu and Kcr were significant enriched compared to the genome-wide averages (Fig. 5c). By contrast, no enrichment was observed for genes concurrently modified by the three marks or by Kcr alone. The analysis supported the hypothesis that Kbu might have a function to poise genes for environmental stress-induced activation.
To study whether Kbu had a function to poise gene activation by internal signals, we identified genes that displayed circadian expression in the rice seedlings at four diurnal time points by RNA-seq (Additional file 1: Table S2) and also analyzed H3K9ac by ChIP-seq of the samples (Additional file 1: Table S1). Genes with significant expression changes (FDR < 0.05, fold change > 2) in more than two of the six comparisons were considered to have diurnal rhythm (Additional file 1: Figure S6A). By these criteria, 2867 diurnal cycling genes (about 13% of expressed genes) were found, and their expression dynamics was tightly correlated with that of H3K9ac levels (Additional file 1: Figure S6B and S6C). Genes marked concurrently by the three marks were enriched for the diurnal expression, whereas those marked by Kbu and Kcr were not enriched (Fig. 5d). Thus, Kbu appeared to be preferentially involved in poising genes for stress-induced activation, whereas H3K9ac is associated with internal circadian gene expression.
SRT2 is involved in erasure of histone Kcr in rice
Lysine butyrylation (Kbu) and crotonylation (Kcr) result in extended hydrocarbon chains that increase the hydrophobicity and bulk of the modified lysine residues in histones exceeding that of lysine acetylation. Whether and how these epigenetic modifications affect gene expression and how they are regulated remains unknown in plants. In this work, we identified a large number of rice histone Kbu and Kcr sites in rice. Although many Kbu and Kcr sites are conserved in mammalian histones, several ones appear to be specific in plants. For instance, H3K14bu seems to exist only in plant and yeast. In addition, we showed that rice Kbu and Kcr are peaked at the TSS regions which is different from human or mouse histone Kbu and Kcr that have also a peak in enhancers at the 5′ region in addition to the TSS . This may be due to architectural difference of cis-regulatory elements between plant and mammalian genes. We showed that histone Kbu and Kcr mark a large number of genes in rice. The majority of Kbu- and Kcr-marked genes are also marked by H3K9ac and are transcriptionally active, which seems to be distinct from results obtained in mammalian cells, in which it has been shown that histone acetylation dynamically competes with butyrylation of highly active gene promoters. The high degree of overlapping between H3K9ac and Kbu/Kcr in the rice genome supports the model in which genomic loci are distinguished not by the presence or absence of a particular acylation but probably by the proportional mixture of distinct acylations. The differential changes of H3K9ac and Kbu/Kcr induced by starvation and submergence suggest that environmental and metabolic cues control proportion of histone acylations in plants, which controls the expression of responsive genes. Several recent studies in animal cells have demonstrated that alterations in the proportional mixture of histone acylations have functional consequences and correlate with distinct physiological states [10, 11, 25–27]. Higher correlations between H3K9ac and expression changes compared to Kbu and Kcr observed in this study suggest that histone acylations may provide a platform for H3K9ac during gene activation in plants. In line with this hypothesis is data showing that Kbu may poise for stress-induced gene activation. The observations that starvation and submergence induced changes of H3K9ac and Kbu/Kcr from different sets of genes and that relatively few genes showed concurrent changes between expression and Kbu or Kcr suggest that the two histone acylation marks and H3K9ac have non-redundant functions in these contexts. Therefore, relative levels of histone acylation and acetylation regulated by environmental and metabolic changes may be a mechanism that fine-tunes regulation of gene expression for plant adaptation.
It has been shown that histone acetyltransferases, such as CBP/p300, catalyze both histone acetyl and non-acetyl acylation [6, 10, 13, 14] and that the different acyl-CoAs compete for binding to the enzyme . Acetyl-CoA produced from cytosolic/nucleic citrate by ATP citrate lyase (ACL) is the major source for histone acetylation. Depletion of acetyl-CoA by knockdown of ACL reduces histone lysine acetylation in both animal and plant cells [28, 29]. Reduction of the cytoplasmic and nuclear pools of acetyl-CoA led to an increase in CBP/p300-catalyzed Kcr in animal cells and the increase in Kcr could be reduced by replenishing the acetyl-CoA pools . Thus, under conditions in which acetyl-CoA is reduced, other acyl-CoA forms will be used more often to modify histones. The relative stability of Kbu and Kcr in rice under starvation and submergence may be due to the possibility that under these conditions tricarboxylic acid cycle (TCA) activity is low resulting in lower cytosolic/nuclear acetyl-CoA, which may be in favor of stabilizing histone Kbu/Kcr acylation. Furthermore, the differential dynamism between H3K9ac and Kbu/Kcr during starvation and submergence may be also related to difference in their erasure mechanism, as only the mutation of ossrt2 had an effect on Kcr, while the three tested HDAC (SRT1, SRT2, and HDA705) all resulted in increases of H3K9ac (Fig. 6). This is in line with the suggestion that non-acetyl histone acylation is mainly removed by sirtuin HDACs, while lysine acetylation can be reversed by both sirtuin and Zn-containing HDAC , although mammalian HDAC1 is shown to have a decrotonylase activity  . However, it is not excluded that anti-Kbu and anti-Kcr antibodies that detect the modifications at all histone sites may not readily detect changes at specific histone lysine sites. A mixture of signals (Kcr or Kbu on different lysines) may change less compared to a single modification (H3K9ac). In addition, Kbu/Kcr/H3K9ac may be regulated at different time points during the stresses. It is also not excluded other factors related to stresses may be involved in the differential changes of histone acylations.
We identified many plant histone Kbu and Kcr sites, which mark a large number of genes in rice. Kbu and Kcr sites are enriched at 5′ region of active genes and most of the H3K9ac-marked genes are co-marked by Kbu and Kcr. However, histone Kbu, Kcr, and H3K9ac display different dynamics in responding to environmental and metabolic cues and have non-redundant functions in gene activation in rice under starvation and submergence. Histone acylation (i.e., Kbu) are likely to poise genes for stress induction and provide a platform for H3K9ac during gene activation in plant. Our data indicate that the proportion of rice histone lysine acetylation and acylation is dynamically regulated by environmental and metabolic signals, which may serve as a mechanism to fine-tune epigenetic control of plant adaptation to environmental variation.
Plant material and treatments
Rice (Oryza sativa) wild type and mutant or RNAi plants were germinated and grown on hormone-free, one-half-strength Murashige and Skoog medium under 16/8 h of light/dark at 30 °C/25 °C. For starvation treatment, 10-day-old wild type DongJin (DJ) seedlings were moved into dark at 6:00 pm and cultured for 48 h. For submergence treatment, 12-day-old wild type seedlings were submerged in distilled water at noon (12:00 am) for 6 h. Seedling leaves of all samples were harvested at 6:00 pm for chromatin and total RNA extraction.
Histone extraction and immunoblot analysis
Histone-enriched fractions were extracted from seedlings of Arabidopsis, rice, maize, and tobacco plants using EpiQuik Total Histone Extraction Kit (Epigentek USA, OP-0006-100). The histone-enriched fractions were used for immunoblot analysis; antibodies used in this study were anti-H3 (ab1791; Abcam), anti-H3K9ace (07-352; Millipore), anti-butyryllysine (PTM-301; PTM Biolabs), and anti-crotonyllysine (PTM-501; PTM Biolabs). Immunoblotting results were quantified using ImageJ (v1.6.0_24).
Vector construction for the CRISPR/Cas9 system
The single-guide RNA (sgRNA) design for the CRISPR/Cas9 system and plasmid construction was conducted as previously described . A web application tool, CRISPR-P (http://cbi.hzau.edu.cn/crispr), was used to select sgRNAs targeting exons of the target genes. The Cas9 destination vector was driven by the maize ubiquitin promoter for expression in rice, and sgRNA expression was driven by the pol III type promoter of U3 sgRNA. Gibson Assembly Cloning method  was used to mobilize sgRNA, after which the binary T-DNA vectors for co-expression of Cas9 and sgRNA were transformed into rice calli. Genomic sequences of OsSRT2 (LOC_Os12g07950) and OsHDA705 (LOC_Os08g25570) were obtained from Rice Genome Annotation Project (http://rice.plantbiology.msu.edu/). Mutation sites of OsSRT2 and OsHDA705 are shown in Additional file 1: Figure S7.
LC–MS/MS and data analysis
Rice seedlings were grinded by liquid nitrogen; core histones were extracted by H2SO4 and digested by trypsin. Lysine crotonylation (Kcr) and butyrylation (Kbu) peptides were enriched by pre-washed antibody beads (PTM Biolabs, Hangzhou). The eluted peptides were cleaned with C18 ZipTips (Millipore) according to the manufacturer’s instructions, followed by analysis with LC–MS/MS. The resulting MS/MS data was processed using MaxQuant with integrated Andromeda search engine (v.184.108.40.206). Tandem mass spectra were searched against UniProt_Oryza sativa database (https://www.uniprot.org/) concatenated with reverse decoy database. False discovery rate thresholds for protein, peptide, and modification site were specified as 1%; minimum peptide length was set at 7; and the site localization probability was set as > 0.75.
Two grams of rice seedling leaves was cross-linked by 1% formaldehyde and used for chromatin extraction. After sonication, chromatin fragments were incubated with antibody (anti-H3K9ac, anti-Kbu, and anti-Kcr)-coated beads (Invitrogen/Life Technologies; 10001D) overnight. Specificity of anti-Kbu and anti-Kcr was tested by dot blots (Additional file 1: Figure S8). After extensive washing, immunoprecipitated chromatin was de-cross-linked and retrieved for qPCR or sequencing. Anti-H3K9ace (07-352; Millipore), anti-butyryllysine (PTM-301; PTM Biolabs), and anti-crotonyllysine (PTM-501; PTM Biolabs) antibodies were used.
ChIP-seq and data analysis
DNA from chromatin immunoprecipitation was used to construct sequencing libraries following the protocol provided by the Illumina TruSeqChIP Sample Prep Set A and sequenced on Illumina HiSeq 2000 with PE 150 method.
Trimmomatic (version 0.32) was used to filter out low-quality reads. Clean reads were mapped to the rice genome (MSU7.0) by Bowtie2 (version 2.2.5), allowing up to two mismatches. Samtools (version 0.1.19) was used to remove potential PCR duplicates, and MACS software (version 1.4.2)  was used to call histone modification peaks by default parameters (bandwidth, 300 bp; model fold, 10, 30; p value, 1.00e−5). Wig files produced by MACS software were used for data visualization by IGV (version 2.3.88). Genes that have modification peaks in gene body or 2 kb upstream in both replicates were defined as modified genes. Tags per kilobase of gene length per million mapped reads (TPM) was used in analyzing gene modification level; genes with TPM fold change > 1.5 in both replicates were considered as differentially modified genes.
In analyzing the relationship between different histone modifications, the rice genome were divided into 1-kb bins and the TPM values in each bin were used to calculate correlation coefficient; data were visualized by R package corrplot (v0.77). In drawing curves of modification distribution over genes, gene body and ± 2-kb regions were divided into 50-bp bins and the average TPM in each bin were plotted. Proportion of each modification was calculated by TPM values in gene body divided by the total acylation TPM value (K9ac + Kbu + Kcr) in each gene.
RT-qPCR and ChIP-qPCR analysis
Total RNA was isolated using TRIzol reagent (Invitrogen). Four micrograms of total RNA was used to synthesize complementary DNA with reverse transcription kit (Promega). Real-time PCR was performed using SYBR Premix ExTaq (TaKaRa) on ABI 7500 real-time PCR system. Levels of rice ACTIN or input (chromatin sample without immunoprecipitation) were used for normalization in RT-qPCR or ChIP-qPCR, respectively. Primers used in RT-qPCR and ChIP-qPCR are listed in Additional file 1: Table S5.
RNA-seq and data analysis
RNA samples were isolated using TRIzol reagent (Invitrogen). The RNA-seq libraries were prepared using the Illumina TruSeq RNA Sample Preparation Kit and sequenced on Illumina HiSeq 2000 with PE150 method.
RNA-seq data were filtered by Trimmomatic (version 0.33) to remove contaminations and low-quality reads. Tophat (version 2.0.14) and Cufflinks (version 2.2.1)  were used to map clean reads to rice genome (MSU 7.0) and calculate differentially expressed genes. Genes with adjusted p value < 0.05 and fold change > 4 in starvation or submergence treatments were considered differentially expressed.
In vitro lysine decrotonylation assays
For OsSRT2 point mutations, forward and reverse primers covering the mutation sites (S123Y and H215Y) were designed (Additional file 1: Table S6) and the fusion PCR was adopted to obtain the mutated OsSRT2 coding sequences, which were then ligated to the protein expression vector PGEX4T-1. OsSRT2-GST, mutants of OsSRT2-GST, and GST proteins were expressed in E. coli and purified using Glutathione Sepharose 4B beads (GE Healthcare 17-0756101). The OsSRT2 full-length cDNA was cloned into pCambia 1301 vector between the 35S promoter and GFP sequences. The resulting vector 35S-SRT2-GFP and empty vector (35S-GFP) were transfected into tobacco leaf cells. After 36 h infection, the leaves were harvested. About 2-g samples were ground into fine powder in liquid nitrogen. Nuclear proteins were extracted as described previously  and suspended in 500 μL reaction buffer.
Reactions were performed as described  in a final volume of 200 μL per well in a 96-well microplate. Briefly, 1 μL of Boc-Lys(crotonyl)-AMC stocking solution (10 mM) and 50 μg purified protein or 80 μL nuclear lysis were added to the reaction buffer (25 mM Tris-HCl, 130 mM NaCl, 3.0 mM KCl, 1 mM MgCl2, 0.1% PEG8000, 500 μM NAD+, pH = 8.0) and incubated at 30 °C for 3 h. The reaction was then stopped by adding 25 μL of trypsin solution (25 mM Tris-HCl, 130 mM NaCl, 3.0 mM KCl, 1 mM MgCl, 30% isopropanol, 0.01 mg/mL trypsin, pH = 8.0). The resulting solvent was mixed and incubated at 30 °C for another 1 h. The fluorescence was analyzed by a fluorescence plate reader with excitation and emission wavelengths at 355 nm and 460 nm, respectively. Boc-Lys(crotonyl)-AMC were gifted by Prof. Yingming Zhao’s lab.
Gene ontology analysis
Gene ontology classification provided in Rice Genome Annotation Project (http://rice.plantbiology.msu.edu/) was used to assign genes to a hierarchical biological process following the criteria of the Gene Ontology Consortium databases (http://www.geneontology.org/external2go/tigr2go). The p value of a particular GO item was calculated with Pearson’s chi-squared test and corrected by FDR, with an FDR cutoff of 0.05 as the significance threshold.
We thank Prof. Yingming Zhao for kindly providing Crotonyl-AMC for enzyme activity tests.
This work was supported by grants from the National Natural Science Foundation of China ; the National Key Research and Development Program of China [2016YFD0100802]; Huazhong Agricultural University Scientific & Technological Self-innovation Foundation [Program No. 2016RC003], and the Fundamental Research Funds for the Central Universities [2662015PY228].
Availability of data and materials
The complete mass spectrometry data produced by this research can be found in ProteomeXchange (http://www.proteomexchange.org/) with the accession number PXD009091 . The high throughput sequencing data generated in this study are deposited into GEO databases with the accession number GSE111308 . Publicly available datasets used are as follows: H3K4me3 and H3K27me3 ChIP-seq data were from GSE19602 , H3K4me2, and H4K12ac ChIP-seq data from GSE26734 , H3K9me2 ChIP-seq data from GSE81436 , drought-treated RNA-seq data from GSE81462 , salt-treated RNA-seq data from GSE101734 , and cold-treated RNA-seq data from GSE67373 .
The review history is available as Additional file 5.
YL and QX performed most of the experiments and analyzed the data. YL, YY, and YZ contributed to the experimentation. ZYC was in charge of MS analysis and provided reagents. DXZ supervised the project, analyzed the data, and wrote the paper with inputs from YL and QX. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Servet C. Conde e Silva N, Zhou DX: histone acetyltransferase AtGCN5/HAG1 is a versatile regulator of developmental and inducible gene expression in Arabidopsis. Mol Plant. 2010;3:670–7.View ArticleGoogle Scholar
- Shen Y, Wei W, Zhou DX. Histone acetylation enzymes coordinate metabolism and gene expression. Trends Plant Sci. 2015;20:614–21.View ArticleGoogle Scholar
- Luo M, Cheng K, Xu Y, Yang S, Wu K. Plant responses to abiotic stress regulated by histone deacetylases. Front Plant Sci. 2017;8:2147.View ArticleGoogle Scholar
- Ost A, Pospisilik JA. Epigenetic modulation of metabolic decisions. Curr Opin Cell Biol. 2015;33:88–94.View ArticleGoogle Scholar
- Shen Y, Issakidis-Bourguet E, Zhou DX. Perspectives on the interactions between metabolism, redox, and epigenetics in plants. J Exp Bot. 2016;67:5291–300.View ArticleGoogle Scholar
- Chen Y, Sprung R, Tang Y, Ball H, Sangras B, Kim SC, Falck JR, Peng J, Gu W, Zhao Y. Lysine propionylation and butyrylation are novel post-translational modifications in histones. Mol Cell Proteomics. 2007;6:812–9.View ArticleGoogle Scholar
- Xie Z, Dai J, Dai L, Tan M, Cheng Z, Wu Y, Boeke JD, Zhao Y. Lysine succinylation and lysine malonylation in histones. Mol Cell Proteomics. 2012;11:100–7.View ArticleGoogle Scholar
- Tan M, Luo H, Lee S, Jin F, Yang JS, Montellier E, Buchou T, Cheng Z, Rousseaux S, Rajagopal N, et al. Identification of 67 histone marks and histone lysine crotonylation as a new type of histone modification. Cell. 2011;146:1016–28.View ArticleGoogle Scholar
- Sabari BR, Zhang D, Allis CD, Zhao Y. Metabolic regulation of gene expression through histone acylations. Nat Rev Mol Cell Biol. 2017;18:90–101.View ArticleGoogle Scholar
- Sabari BR, Tang Z, Huang H, Yong-Gonzalez V, Molina H, Kong HE, Dai L, Shimada M, Cross JR, Zhao Y, et al. Intracellular crotonyl-CoA stimulates transcription through p300-catalyzed histone crotonylation. Mol Cell. 2015;58:203–15.View ArticleGoogle Scholar
- Xie Z, Zhang D, Chung D, Tang Z, Huang H, Dai L, Qi S, Li J, Colak G, Chen Y, et al. Metabolic regulation of gene expression by histone lysine beta-hydroxybutyrylation. Mol Cell. 2016;62:194–206.View ArticleGoogle Scholar
- Wei W, Liu X, Chen J, Gao S, Lu L, Zhang H, Ding G, Wang Z, Chen Z, Shi T, et al. Class I histone deacetylases are major histone decrotonylases: evidence for critical and broad function of histone crotonylation in transcription. Cell Res. 2017;27:898–915.View ArticleGoogle Scholar
- Ogryzko VV, Schiltz RL, Russanova V, Howard BH, Nakatani Y. The transcriptional coactivators p300 and CBP are histone acetyltransferases. Cell. 1996;87:953–9.View ArticleGoogle Scholar
- Kaczmarska Z, Ortega E, Goudarzi A, Huang H, Kim S, Marquez JA, Zhao Y, Khochbin S, Panne D. Structure of p300 in complex with acyl-CoA variants. Nat Chem Biol. 2017;13:21–9.View ArticleGoogle Scholar
- Cheng Z, Tang Y, Chen Y, Kim S, Liu H, Li SS, Gu W, Zhao Y. Molecular characterization of propionyllysines in non-histone proteins. Mol Cell Proteomics. 2009;8:45–52.View ArticleGoogle Scholar
- Tcherkez G, Boex-Fontvieille E, Mahe A, Hodges M. Respiratory carbon fluxes in leaves. Curr Opin Plant Biol. 2012;15:308–14.View ArticleGoogle Scholar
- Lee CP, Eubel H, Millar AH. Diurnal changes in mitochondrial function reveal daily optimization of light and dark respiratory metabolism in Arabidopsis. Mol Cell Proteomics. 2010;9:2125–39.View ArticleGoogle Scholar
- Zhong X, Zhang H, Zhao Y, Sun Q, Hu Y, Peng H, Zhou DX. The rice NAD(+)-dependent histone deacetylase OsSRT1 targets preferentially to stress- and metabolism-related genes and transposable elements. PLoS One. 2013;8:e66807.View ArticleGoogle Scholar
- Pena-Castro JM, van Zanten M, Lee SC, Patel MR, Voesenek LA, Fukao T, Bailey-Serres J. Expression of rice SUB1A and SUB1C transcription factors in Arabidopsis uncovers flowering inhibition as a submergence tolerance mechanism. Plant J. 2011;67:434–46.View ArticleGoogle Scholar
- Bao X, Wang Y, Li X, Li XM, Liu Z, Yang T, Wong CF, Zhang J, Hao Q, Li XD. Identification of ‘erasers’ for lysine crotonylated histone marks using a chemical proteomics approach. Elife. 2014;3:e02999.Google Scholar
- Huang L, Sun Q, Qin F, Li C, Zhao Y, Zhou DX. Down-regulation of a SILENT INFORMATION REGULATOR2-related histone deacetylase gene, OsSRT1, induces DNA fragmentation and cell death in rice. Plant Physiol. 2007;144:1508–19.View ArticleGoogle Scholar
- Zhang H, Zhao Y, Zhou DX. Rice NAD+-dependent histone deacetylase OsSRT1 represses glycolysis and regulates the moonlighting function of GAPDH as a transcriptional activator of glycolytic genes. Nucleic Acids Res. 2017;45:12241–55.View ArticleGoogle Scholar
- Hu Y, Qin F, Huang L, Sun Q, Li C, Zhao Y, Zhou DX. Rice histone deacetylase genes display specific expression patterns and developmental functions. Biochem Biophys Res Commun. 2009;388:266–71.View ArticleGoogle Scholar
- Zhao J, Li M, Gu D, Liu X, Zhang J, Wu K, Zhang X, Teixeira da Silva JA, Duan J. Involvement of rice histone deacetylase HDA705 in seed germination and in response to ABA and abiotic stresses. Biochem Biophys Res Commun. 2016;470:439–44.View ArticleGoogle Scholar
- Kebede AF, Nieborak A, Shahidian LZ, Le Gras S, Richter F, Gomez DA, Baltissen MP, Meszaros G, Magliarelli HF, Taudt A, et al. Histone propionylation is a mark of active chromatin. Nat Struct Mol Biol. 2017;24:1048–56.View ArticleGoogle Scholar
- Goudarzi A, Zhang D, Huang H, Barral S, Kwon OK, Qi S, Tang Z, Buchou T, Vitte AL, He T, et al. Dynamic competing histone H4 K5K8 acetylation and butyrylation are hallmarks of highly active gene promoters. Mol Cell. 2016;62:169–80.View ArticleGoogle Scholar
- Li Y, Sabari BR, Panchenko T, Wen H, Zhao D, Guan H, Wan L, Huang H, Tang Z, Zhao Y, et al. Molecular coupling of histone crotonylation and active transcription by AF9 YEATS domain. Mol Cell. 2016;62:181–93.View ArticleGoogle Scholar
- Wellen KE, Hatzivassiliou G, Sachdeva UM, Bui TV, Cross JR, Thompson CB. ATP-citrate lyase links cellular metabolism to histone acetylation. Science. 2009;324:1076–80.View ArticleGoogle Scholar
- Chen C, Li C, Wang Y, Renaud J, Tian G, Kambhampati S, Saatian B, Nguyen V, Hannoufa A, Marsolais F, et al. Cytosolic acetyl-CoA promotes histone acetylation predominantly at H3K27 in Arabidopsis. Nat Plants. 2017;3:814–24.View ArticleGoogle Scholar
- He Y, Zhang T, Yang N, Xu M, Yan L, Wang L, Wang R, Zhao Y. Self-cleaving ribozymes enable the production of guide RNAs from unlimited choices of promoters for CRISPR/Cas9 mediated genome editing. J Genet Genomics. 2017;44:469–72.View ArticleGoogle Scholar
- Gibson DG, Young L, Chuang RY, Venter JC, Hutchison CA 3rd, Smith HO. Enzymatic assembly of DNA molecules up to several hundred kilobases. Nat Methods. 2009;6:343–5.View ArticleGoogle Scholar
- Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nusbaum C, Myers RM, Brown M, Li W, Liu XS. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008;9:R137.View ArticleGoogle Scholar
- Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and cufflinks. Nat Protoc. 2012;7:562–78.View ArticleGoogle Scholar
- Lu Y, Xu QT, Liu Y, Yu Y, Chen ZY, Zhao Y, Zhou DX: Histone lysine acetylation, butyrylation, and crotonylation dynamics and functional interplay in rice under starvation and submergence. EBI.PRIDE. PXD009091. https://www.ebi.ac.uk/pride/archive/projects/PXD009091
- Lu Y, Xu QT: Dynamics and functional interplay of histone lysine butyrylation, crotonylation, and acetylation in rice under starvation and submergence. NCBI GEO. GSE111308. https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE111308
- He G, Zhu X, Elling AA, Chen L, Wang X, Guo L, et al: Global epigenetic and transcriptional trends among two rice subspecies and their reciprocal hybrids. NCBI GEO. GSE19602. https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE19602
- Zhang WL, Wu YF, Schnable, JC, Zeng, ZX, Freeling, M, Crawford GE, Jiang JM: High-resolution mapping of open chromatin in the rice genome. NCBI GEO. GSE26734. https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE26734
- Tan F, Zhou C, Zhou QW, Zhou SL, Yang WJ, Zhao Y, Li GL, Zhou DX: Analysis of chromatin regulators reveals specific features of rice DNA methylation pathways. NCBI GEO. GSE81436. https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE81436
- Zong W, Tang N, Yang J, Peng L, Ma SQ, Xu Y, Li GY, Xiong LZ: Feedback regulation of ABA signaling and biosynthesis by OsbZIP23 that targets diverse drought resistance related genes in rice. NCBI GEO. GSE81462. https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE81462
- Wang SS, Cao M, Ma X, Chen WK, Zhao J, Sun CQ, Tan LB, Liu FX:RNA-seq data of two half-leaf seedlings under salt treatment and control condition. NCBI GEO. GSE101734.https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE101734.
- Zhang JY, Luo W, Zhao Y, Xu YY, Song SH, Chong K: A global profiling of gene expression in chilling stress in rice. NCBI GEO. GSE67373. https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE67373