Dynamic changes in genome-wide histone H3 lysine 4 methylation patterns in response to dehydration stress in Arabidopsis thaliana
© van Dijk et al; licensee BioMed Central Ltd. 2010
Received: 29 June 2010
Accepted: 5 November 2010
Published: 5 November 2010
The molecular mechanisms of genome reprogramming during transcriptional responses to stress are associated with specific chromatin modifications. Available data, however, describe histone modifications only at individual plant genes induced by stress. We have no knowledge of chromatin modifications taking place at genes whose transcription has been down-regulated or on the genome-wide chromatin modification patterns that occur during the plant's response to dehydration stress.
Using chromatin immunoprecipitation and deep sequencing (ChIP-Seq) we established the whole-genome distribution patterns of histone H3 lysine 4 mono-, di-, and tri-methylation (H3K4me1, H3K4me2, and H3K4me3, respectively) in Arabidopsis thaliana during watered and dehydration stress conditions. In contrast to the relatively even distribution of H3 throughout the genome, the H3K4me1, H3K4me2, and H3K4me3 marks are predominantly located on genes. About 90% of annotated genes carry one or more of the H3K4 methylation marks. The H3K4me1 and H3K4me2 marks are more widely distributed (80% and 84%, respectively) than the H3K4me3 marks (62%), but the H3K4me2 and H3K4me1 levels changed only modestly during dehydration stress. By contrast, the H3K4me3 abundance changed robustly when transcripts levels from responding genes increased or decreased. In contrast to the prominent H3K4me3 peaks present at the 5'-ends of most transcribed genes, genes inducible by dehydration and ABA displayed atypically broader H3K4me3 distribution profiles that were present before and after the stress.
A higher number (90%) of annotated Arabidopsis genes carry one or more types of H3K4me marks than previously reported. During the response to dehydration stress the changes in H3K4me1, H3K4me2, and H3K4me3 patterns show different dynamics and specific patterns at up-regulated, down-regulated, and unaffected genes. The different behavior of each methylation mark during the response process illustrates that they have distinct roles in the transcriptional response of implicated genes. The broad H3K4me3 distribution profiles on nucleosomes of stress-induced genes uncovered a specific chromatin pattern associated with many of the genes involved in the dehydration stress response.
Plants respond to external signals by activating signaling pathways that rapidly alter physiological reactions and transcription rates of responding genes. The molecular mechanisms regulating gene expression are coordinated at the genome level where the accessibility of DNA sequences is determined by the structure of chromatin. Chemical modifications of the histone amino-terminal tails are implicated in facilitating nucleosome remodeling and repositioning, and in the recruitment of specific transcription factors. Changes in DNA methylation and/or histone modifications are associated with altered gene expression under various stresses in plant and cell culture systems [1–5]. Changes in nucleosome occupancy and in the levels of histone H3 tri-methylation (lysine 4) or acetylation (lysine 9, 14 or 27) during dehydration stress occurred at four inducible Arabidopsis genes . Dehydration stress induces biosynthesis of abscisic acid (ABA), which as a second messenger activates signaling cascades that trigger stomatal closure and altered gene expression via chromatin remodeling. In turn, ABA regulation is achieved through genetic and chromatin modification mechanisms [6, 7]. Thus, chromatin structure and chromatin modifications are emerging as critical factors in plants' responses to environmental cues.
It is important to note, however, that currently available data describe histone modifications only at individual stress-induced plant genes. We also have limited knowledge of chromatin modifications taking place at genes whose transcription has been down-regulated. Moreover, there are no data on global modification trends occurring during the reprogramming of the entire genome in response to stress.
Here, we provide a genome-wide view and analysis of the histone H3 lysine 4 mono-, di-, and tri-methylation (H3K4me1, H3K4me2, H3K4me3, respectively) patterns in chromatin isolated from Arabidopsis rosette leaves before and after dehydration stress. Genome-wide transcript patterns in watered and dehydration stressed plants were compared with changes in the H3K4me1, H3K4me2, and H3K4me3 levels of nucleosomes associated with responding genes. Using chromatin immunoprecipitation (ChIP) with H3K4 methylation specific antibodies and genome-wide sequencing (ChIP-Seq), we revealed different dynamics and different magnitudes of the changes in H3K4me1, H3K4me2, and H3K4me3 profiles taking place upon dehydration stress. We demonstrate specific patterns of the H3K4me1, H3K4me2, and H3K4me3 distributions at up-regulated, down-regulated, and unaffected genes during the stress.
Plant materials and growth conditions
Col-0 plants were grown in 12 h light for 4 weeks in pots with soil covered with miracloth to prevent soil contamination of leaf tissues at harvest. Control plants continued to receive water while water was withheld from water-deficit treated plants until a relative water content (RWC) of 65% was reached in 4 to 6 days. Vegetative rosettes were harvested and frozen for RNA isolation or immediately fixed in formaldehyde for isolation of chromatin as described below.
Isolation and immunoprecipitation of chromatin
Watered or water-deficit treated vegetative tissues were fixed in 1% formaldehyde for 10 min after vacuum infiltration of the tissues in phosphate buffered saline (PBS) with 0.007% silwet. Fixation was stopped by quenching with 125 mM glycine. Next, the tissue was rinsed twice in ice cold water, blotted dry and flash frozen in liquid nitrogen. Chromatin was isolated and resuspended in 300 μl nuclei lysis buffer as described by  and sonicated for four 15 sec pulses using a Branson sonifier 450 at an output between 10 - 15%. The sheared chromatin was used in chromatin immunoprecipitation assays as described  with some modifications. In brief, 100 μl of sheared chromatin was diluted in 900 μl ChIP dilution buffer and precleared, using 15 μl of protein A Dynabeads (Invitrogen, Carlsbad, CA). Protein A beads were collected using a magnet, and the precleared chromatin was incubated overnight at 4 °C with antibodies against either H3 (6 μL, Abcam ab1791, Cambridge, MA), mono-methylated H3K4 (10 μL, Abcam ab8895, Cambridge, MA), dimethylated H3K4 (6 μL, Upstate 07-212, Millipore, Billerica, MA), or trimethylated H3K4 (6 μL, Abcam ab8580, Cambridge, MA). Western blotting of the H3K4 methylation antibodies against mono-, di- or tri-methylated H3K4 or monomethylated H3K9 synthetic peptides was performed and confirmed their expected specificity (Additional file 1, Figure S1). A sample without antibody addition was also used as a control. Next the antibody/chromatin complex was precipitated by incubating the mix with 40 μl of protein A beads for 1 h followed by collection of the beads using a magnet. The beads were washed, chromatin eluted, and DNA de-crosslinked using procedures as described . DNA was recovered by phenol/chloroform extraction and ethanol precipitation in the presence of Novagen pellet paint (EMD Chemicals, Gibbstown, NJ) and resuspended in 30 μL of TE buffer. Purified DNA was end modified, ligated to amplification primers, amplified, and sequenced according to the manufacturer's protocols (Illumina, San Diego, CA).
RNA isolation and labeling
Total RNA was extracted from the frozen rosette tissues using grinding and TRIzol (Invitrogen, Carlsbad, CA) reagent. Three independent replicates of watered and dehydration-stressed tissues were used. Fifteen μg of total RNA from each sample was used to produce labeled cRNA for hybridization to Affymetrix Arabidopsis ATH microarrays using the manufacturer's protocols (Affymetrix, Santa Clara, CA).
Microarray data analysis
Affymetrix® GeneChip Operating Software (GCOS) generated the image and probe set signal values as CEL files, which were loaded into Rosetta Resolver® (version 6.0, build 126.96.36.1991; Rosetta Inpharmatics LLC, Seattle, WA). The water deficit and control condition replicate arrays were combined into their respective intensity profiles, and subsequently used to generate ratio experiments in Rosetta Resolver and analyzed using the software's background correction, normalization, and proprietary ratio error model . The gene expression arrays have been submitted to the NCBI Gene Expression Omnibus under GSE11538 (part of SuperSeries record GSE11658).
Computational Methods for the analysis of sequence data
The primary analysis of the Illumina Genome Analyzer output was performed by the manufacturer's ELAND program. ELAND mapped sequencing reads of 35 bases uniquely to the Arabidopsis genome with no more than 2 mismatches. As a quality control, we used BLAST  to re-map all reads to the Arabidopsis genome  with a word length of 9 bases and low-complexity filtering turned off. This analysis revealed a small number of further non-unique sequencing reads. From the remaining reads that uniquely map to the genome with no more than 2 mismatches, the coverage of each genomic position by the sequencing reads was computed. The DNA sequencing reads for total genomic DNA from Arabidopsis thaliana cultivar Col-0 were obtained and mapped to the Arabidopsis thaliana Col-0 genome, which is used as the control data set. To check for biases in Illumina sequencing, the distribution of sequencing coverage along the genome was calculated. Statistical analyses were performed in MATLAB (MathWorks, Nantucket, MA). Tests for the statistical significance in the difference in the coverage profiles of H3K4me3 along the normalized gene length between the sets of all expressed and ABA induced genes under both watered and soil water deficit conditions were performed by the non-parametric Wilcoxon rank sum test on the total (peak-normalized) coverage over the normalized gene length for individual genes.
Coverage peaks were identified as maximal scoring segments using a modified dynamic programming algorithm, analogous to pairwise local alignments . Background noise was estimated as the mean coverage of the intergenic regions (ranging from 1.10 to 3.25 depending on the sequencing run) and was subtracted from the raw coverage values. We extended peaks from each genomic position that had an above background coverage so as to reach the maximum of the cumulative coverage, analogous to the cumulative score of a BLAST High-Scoring Segment Pair. During peak extension, we allowed for local drops in the peaks provided that they were compensated by subsequent above-background regions.
Position-wise coverage data for the three H3K4 methylation levels for watered and dehydration-stressed conditions, peak calls, our gene expression results, exons, UTRs, alternative transcription start sites are displayed in the implementation of the UC Santa Cruz Genome Browser  at the Bioinformatics Core Facility at the University of Nebraska-Lincoln. Exons, transcription start sites and untranslated regions were downloaded from the TAIR ftp://ftp.Arabidopsis.org/home/tair/Genes/. The results can be accessed at: http://bioinformatics.unl.edu/h3k4me/cgi-bin/hgGateway. The ChIP-Seq data has been deposited at the NCBI Gene Expression Omnibus under series record GSE11657 (part of SuperSeries record GSE11658).
Changes in gene expression in plants exposed to soil water deficit conditions
Plants were watered or subjected to dehydration stress by withholding water until a relative water content (RWC) of 65% was reached. A RWC of 65% causes leaf wilting and is considered a moderate level of dehydration stress . The mRNAs from rosette leaves of non-stressed (control) and dehydration-stressed Arabidopsis plants were analyzed by microarrays (Additional file 2, Table S1). Further analysis indicated 640 genes increased, and 525 genes decreased, their transcript levels by 4-fold or more (p ≤ 0.001) as result of the dehydration stress (Additional file 3, Table S2). Analysis of the most highly represented gene categories (using the Gene Ontology classification) indicated that the majority of the genes with altered expression was involved in various stress-responses (Additional file 4, Table S3). For example, within the list of genes responding to dehydration-stress, many genes have been also recognized as components of pathways stimulated by cold, temperature, light, and radiation (~750 genes), by hormones, endogenous, and chemical stimuli (~300 genes), as well as genes encoding transcription factors and regulators (131 genes). Genes for metabolic enzymes, and for oxidoreductase functions, were also affected (Additional file 4, Table S3). Among the genes up-regulated by dehydration stress, well-known genes implicated in this response in other studies were present, including genes ABF3 (ABA responsive element binding factor 3, AT4G34000, ), CBF4 (C-repeat binding factor 4, AT5G51990 ), RD26 (Responsive to desiccation 26, AT4G27410, ), RD29A (AT5G52310), RD29B (AT5G52300), ATHB7 (Arabidopsis thaliana homeobox 7, AT2G46680) and ATHB12 (AT3G61890; ).
Genome-wide distributions of H3 and mono-, di-, and tri-methylated H3K4 in watered and in dehydration stressed plants
Percentage of genes with H3K4 methylation regions1
Treatment and type of H3K4 methylation
Number of genes
Percentage of genes
Genes with one or more types of H3K4 methylation
Genes lacking H3K4 methylation
Genome-wide distribution of histone H3, H3K4me1, me2 and me3 at gene sequences: correlation with transcription activity
The profiles of mono-, di-, and tri-methylated H3K4 displayed specific patterns associated with transcript levels. H3K4me1 levels remained fairly constant throughout genes (Figure 2B). One pattern of interest was that the 5'-regions of the genes with the lowest expression levels had higher H3K4me1 levels than more highly expressed genes (Figure 2B). The profiles of H3K4me2 distribution had no correlation with transcript levels, and were relatively constant throughout the transcribed regions (Figure 2C). As such, H3K4me2 marked 84% of genes (Table 1) irrespective of the transcript levels produced from these genes. Both H3K4me1 and H3K4me2 were fairly constant across the middle and end of the genes, and accurately marked the ends of the genes.
In contrast, the H3K4me3 distribution showed a clear peak around 300 bp downstream of the TSS. However, the magnitudes of signal values for the five expression-based quintiles were in reverse to those displayed by H3K4me1: the highest H3K4me3 peaks were for genes with the highest transcript levels. The H3K4me3 signal at the 5'-ends of the least actively transcribed genes displayed the lowest coverage values (Figure 2D). Characteristically, the peak of H3K4me3 distribution tapers off towards the 3' end of the genes and reaches intergenic levels at the 3'-gene ends.
H3K4 methylation changes when transcript levels change
The most significant changes were displayed by the H3K4 tri-methylation profiles. The median and the majority of the genes within one standard deviation of the median increased for all up-regulated quintiles and decreased for all down regulated quintiles (Figure 3). In addition, the largest increase in H3K4me3 levels was found on nucleosomes of the genes in the highest expression quintile, indicating that the magnitude of the H3K4me3 increase was positively associated with increased transcript levels. Similarly, the down-regulated quintile with the largest decrease in transcript levels had the largest decrease in H3K4me3 levels (Figure 3).
As a control for our method of analysis, we analyzed genes that did not change transcript levels and found that these genes did not show significant variation in their H3K4 methylation levels. Correlation coefficients of 0.95, 0.88, and 0.93 were observed between the watered and dehydration-stressed samples for the H3K4 mono-, di-, or tri-methylation comparisons of genes that did not change transcript levels, respectively. These high correlation coefficients indicated that there was very little change in H3K4 methylation in this set of genes (Additional file 8, Figure S2).
H3K4 methylation profiles of selected genes during watered or dehydration stressed conditions
Similar tendencies were demonstrated by other dehydration stress responsive genes as well. Six representative chromosomal regions, representing genes that were induced, repressed, or unaffected by dehydration stress conditions, are shown in a UCSC Genome Browser view (Additional file 9, Figure S3). For induced genes, the H3K4me3 levels increased during dehydration stress as their expression increased, with a corresponding decrease in H3K4me1, while H3K4me2 levels showed little change (Additional file 9, Figure S3A, S3B). The two dehydration-stress repressed genes showed decreased H3K4me3 levels, and a slight increase in H3K4me1 (Additional file 9, Figure S3C and additional file 10, Figure S3D). By comparison, genes without changes in gene expression did not display significant changes in any of the H3K4 methylation-forms (Additional file 10, Figure S3E, S3F). Collectively, the changes, or lack of changes, in the H3K4 methylation patterns displayed at individual gene loci were consistent with the overall genomic patterns as described above.
Dehydration stress responsive genes show broader H3K4me3 distribution profiles
Here we report the first whole genome map of the histone H3K4 methylation patterns of chromatin during a plant's response to dehydration stress. Using the ChIP-Seq approach, we found that the distribution of mono-, di-, and tri-methylation of H3K4 in Arabidopsis were predominantly located on genes, in general agreement with prior studies . However, in watered plants we identified 21,875 genes carrying H3K4me1 versus 15,475 genes found by the ChIP-chip method; 22,708 genes with the H3K4me2 mark versus 12,781, and 16,769 H3K4me3-tagged genes versus 15,894 genes reported in . Furthermore, our results indicated that only 10% of expressed genes did not contain at least one H3K4me type, while ~1/3 of the genes was found to be free of H3K4 methylation by the ChIP-chip approach. These differences might reflect different sensitivities of the two methods and/or different tissue distributions in the whole seedlings versus rosette leaves analyzed (this study). The gene-associated location of H3K4me modifications did not change during the dehydration stress-responding process. Within this gene-centric framework, however, the abundances of the individual forms of H3K4 methylation displayed specific profile changes associated with altered transcription from these genes.
H3 distribution and nucleosome free regions
The genome-wide distribution of histone H3 was used to gain insight to whether a region lacked nucleosomes or lacked H3K4 methylation of those nucleosomes. The H3 distribution showed a fairly continuous presence across genic and intergenic regions. By contrast, the H3K4 methylated forms were practically absent from intergenic regions. In agreement with observations made in yeast  and Arabidopsis , we also observed an H3 enrichment in gene regions relative to intergenic regions. Furthermore, in agreement with other genomic analyses in Arabidopsis [24–26], our analysis of the H3 profiles in genes grouped by transcript levels, revealed regions with low H3 presence, "nucleosome free regions" (NFRs), identified immediately upstream of the TSS for more highly expressed genes. NFRs have also been observed in yeast [27–29], Drosophila , and human cells . Interestingly, nucleosome density has also been observed to be different in exon and intron gene regions . It has been suggested that the nucleosome occupancy is influenced by the DNA sequence, which might facilitate nucleosome removal upon transcription . Our finding that NFRs are associated with transcriptionally active genes is consistent with the idea of nucleosome displacement by transcription factors  or complexes involved in transcription . Similarly, transcription from the Arabidopsis AP1, a gene required for the initiation of flowering, requires the reprogramming of the AP1 locus to an actively transcribed state that is accompanied by the removal of a nucleosome from the transcription start site and its dynamic re-positioning in a developmentally regulated process .
H3K4me1 distribution changes
The majority of genes with increased transcript levels showed modest reductions in H3K4me1 levels. Conversely, the most strongly down-regulated genes showed an increase in H3K4me1 marks, while genes with smaller reductions in transcript levels did not detectably change their average H3K4me1 levels. A negative association of H3K4me1 levels with increased transcription might be due to conversion to the more abundant H3K4me2 and H3K4me3 forms at higher transcription rates. Alternatively, higher H3K4me1 levels, particularly near the TSS, might repress transcription. A repressive role has been reported in Chlamydomonas where H3K4 mono-methylation was found associated with repressed transgenes and transposable elements. The loss of this mark in Chlamydomonas mutant Mut11 resulted in expression of these genes . The possible role of H3K4me1 in gene repression might be unique to the plant kingdom as H3K4me1 levels are associated with enhancers in animal cells .
H3K4me2 distribution changes
In our studies, 84% of Arabidopsis genes carried H3K4me2 marks, and on average, these marks were distributed relatively uniformly along gene sequences. In agreement with an earlier genomic study , H3K4me2 levels were not associated with transcript levels in watered plants. In response to dehydration stress, H3K4me2 levels were only slightly reduced when transcript levels increased, and showed little change when transcript levels were reduced. Previous ChIP-PCR analysis of individual genes in Arabidopsis found that H3K4me2 marked 5' boundaries of the tested genes regardless of expression levels . These observations suggest H3K4me2 marks do not correlate with transcript levels and may play a role in gene regulation different from the roles of mono-, and tri-methylated nucleosomes.
H3K4me3 distribution changes
H3K4me3 modification was predominantly associated with actively transcribed genes. Our observation that higher H3K4me3 levels were associated with higher transcript levels was in agreement with studies in yeast , rice , Arabidopsis [22, 24], and human cells [40, 41]. This tendency was dynamic, as we found up-regulated genes had increased H3K4me3 levels, and down-regulated genes had reduced H3K4me3 levels. These findings are in agreement with numerous results from individual plant genes, where H3K4me3 was found to increase with increased transcription at dehydration induced promoters , or the phaseolin , alcohol dehydrogenase 1 (ADH1) and pyruvate decarboxylase 1 (PDC1) genes . The magnitude of the change, however, only weakly correlated with the magnitude of the change in transcript levels as shown by the wide range of H3K4me3 values in each expression quintile (Figure 3).
The genome-wide H3K4me3 distribution in Arabidopsis displays a peak at +300 bp relative to the TSS (Figure 2D, ), similar to the location observed in human cells , but slightly upstream of the peak location in rice observed at about +500 bp . A surprising deviation from this pattern was the broad H3K4me3 distribution profiles present on many dehydration/ABA inducible genes. Interestingly, this broad pattern was present on many of these genes prior to their induction, although the lower H3K4me3 levels made this pattern less striking than in the induced state. For example, both RD29A and RD29B had low levels of H3K4me3 in a broad distribution pattern that was similar to the broad pattern of the higher levels of H3K4me3 observed during dehydration stress (Figure 4). Although the occurrence of this pattern was statistically significant, its functional significance is unclear.
Although there is a strong association of H3K4me3 with transcriptionally active genes, the biochemical roles of H3K4me3 are emerging as occurring predominantly after transcription initiation [43–45]. H3K4me3 marks are recognized by chromatin remodeling factors facilitating transcription by altering the structure, composition, and positioning of nucleosomes, by components of the spliceosome, and by proteins involved in mRNA capping and stability . Recognition and binding to H3K4me3 has been traced to plant homeodomains or chromodomains present in these proteins . Additionally, in support of a post-transcriptional function, the yeast SET1 H3K4 methyltransferase responsible for H3K4me3 modification is recruited to the +300 bp region of genes after transcription initiation in yeast . Elevated H3K4me3 levels occurring on genes with increased transcript production might effectively recruit more chromosomal and RNA processing proteins, thereby facilitating the increased transcriptional activity on these induced genes.
We thank Dr. Richard Mott for communicating unpublished data. This work was supported by National Science Foundation grant EPS-0701892.
- Steward N, Ito M, Yamaguchi Y, Koizumi N, Sano H: Periodic DNA methylation in maize nucleosomes and demethylation by environmental stress. J Biol Chem. 2002, 277: 37741-37746. 10.1074/jbc.M204050200.PubMedView ArticleGoogle Scholar
- Manzanero S, Rutten T, Kotseruba V, Houben A: Alterations in the distribution of histone H3 phosphorylation in mitotic plant chromosomes in response to cold treatment and the protein phosphatase inhibitor cantharidin. Chromosome Res. 2002, 10: 467-476. 10.1023/A:1020940313841.PubMedView ArticleGoogle Scholar
- Chua YL, Watson LA, Gray JC: The transcriptional enhancer of the pea plastocyanin gene associates with the nuclear matrix and regulates gene expression through histone acetylation. Plant Cell. 2003, 15: 1468-1479. 10.1105/tpc.011825.PubMedPubMed CentralView ArticleGoogle Scholar
- Tsuji H, Saika H, Tsutsumi N, Hirai A, Nakazono M: Dynamic and reversible changes in histone H3-Lys4 methylation and H3 acetylation occurring at submergence-inducible genes in rice. Plant Cell Physiol. 2006, 47: 995-1003. 10.1093/pcp/pcj072.PubMedView ArticleGoogle Scholar
- Sokol A, Kwiatkowska A, Jerzmanowski A, Prymakowska-Bosak M: Up-regulation of stress-inducible genes in tobacco and Arabidopsis cells in response to abiotic stresses and ABA treatment correlates with dynamic changes in histone H3 and H4 modifications. Planta. 2007, 227: 245-254. 10.1007/s00425-007-0612-1.PubMedView ArticleGoogle Scholar
- Kim JM, To TK, Ishida J, Morosawa T, Kawashima M, Matsui A, Toyoda T, Kimura H, Shinozaki K, Seki M: Alterations of lysine modifications on the histone H3 N-tail under drought stress conditions in Arabidopsis thaliana. Plant Cell Physiol. 2008, 49: 1580-1588. 10.1093/pcp/pcn133.PubMedView ArticleGoogle Scholar
- Chinnusamy V, Gong Z, Zhu JK: Abscisic acid-mediated epigenetic processes in plant development and stress responses. J Integr Plant Biol. 2008, 50: 1187-1195. 10.1111/j.1744-7909.2008.00727.x.PubMedPubMed CentralView ArticleGoogle Scholar
- Bowler C, Benvenuto G, Laflamme P, Molino D, Probst AV, Tariq M, Paszkowski J: Chromatin techniques for plant cells. Plant J. 2004, 39: 776-789. 10.1111/j.1365-313X.2004.02169.x.PubMedView ArticleGoogle Scholar
- Weng L, Dai H, Zhan Y, He Y, Stepaniants SB, Bassett DE: Rosetta error model for gene expression analysis. Bioinformatics. 2006, 22: 1111-1121. 10.1093/bioinformatics/btl045.PubMedView ArticleGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMedPubMed CentralView ArticleGoogle Scholar
- Swarbreck D, Wilks C, Lamesch P, Berardini TZ, Garcia-Hernandez M, Foerster H, Li D, Meyer T, Muller R, Ploetz L, Radenbaugh A, Singh S, Swing V, Tissier C, Zhang P, Huala E: The Arabidopsis Information Resource (TAIR): gene structure and function annotation. Nucleic Acids Res. 2008, 36: D1009-1014. 10.1093/nar/gkm965.PubMedPubMed CentralView ArticleGoogle Scholar
- Smith TF, Waterman MS: Identification of common molecular subsequences. J Mol Biol. 1981, 147: 195-197. 10.1016/0022-2836(81)90087-5.PubMedView ArticleGoogle Scholar
- Kuhn RM, Karolchik D, Zweig AS, Trumbower H, Thomas DJ, Thakkapallayil A, Sugnet CW, Stanke M, Smith KE, Siepel A, Rosenbloom KR, Rhead B, Raney BJ, Pohl A, Pedersen JS, Hsu F, Hinrichs AS, Harte RA, Diekhans M, Clawson H, Bejerano G, Barber GP, Baertsch R, Haussler D, Kent WJ: The UCSC genome browser database: update 2007. Nucleic Acids Research. 2007, 35: D668-673. 10.1093/nar/gkl928.PubMedPubMed CentralView ArticleGoogle Scholar
- Kawaguchi R, Girke T, Bray EA, Bailey-Serres J: Differential mRNA translation contributes to gene regulation under non-stress and dehydration stress conditions in Arabidopsis thaliana. Plant J. 2004, 38: 823-839. 10.1111/j.1365-313X.2004.02090.x.PubMedView ArticleGoogle Scholar
- Kang JY, Choi HI, Im MY, Kim SY: Arabidopsis basic leucine zipper proteins that mediate stress-responsive abscisic acid signaling. Plant Cell. 2002, 14: 343-357. 10.1105/tpc.010362.PubMedPubMed CentralView ArticleGoogle Scholar
- Haake V, Cook D, Riechmann JL, Pineda O, Thomashow MF, Zhang JZ: Transcription factor CBF4 is a regulator of drought adaptation in Arabidopsis. Plant Physiol. 2002, 130: 639-648. 10.1104/pp.006478.PubMedPubMed CentralView ArticleGoogle Scholar
- Fujita M, Fujita Y, Maruyama K, Seki M, Hiratsu K, Ohme-Takagi M, Tran LS, Yamaguchi-Shinozaki K, Shinozaki K: A dehydration-induced NAC protein, RD26, is involved in a novel ABA-dependent stress-signaling pathway. Plant J. 2004, 39: 863-876. 10.1111/j.1365-313X.2004.02171.x.PubMedView ArticleGoogle Scholar
- Olsson AS, Engstrom P, Soderman E: The homeobox genes ATHB12 and ATHB7 encode potential regulators of growth in response to water deficit in Arabidopsis. Plant Mol Biol. 2004, 55: 663-677. 10.1007/s11103-004-1581-4.PubMedView ArticleGoogle Scholar
- Karolchik D, Bejerano G, Hinrichs AS, Kuhn RM, Miller W, Rosenbloom KR, Zweig AS, Haussler D, Kent WJ: Comparative genomic analysis using the UCSC genome browser. Methods Mol Biol. 2007, 395: 17-34.PubMedView ArticleGoogle Scholar
- Whitehouse I, Rando OJ, Delrow J, Tsukiyama T: Chromatin remodelling at promoters suppresses antisense transcription. Nature. 2007, 450: 1031-1035. 10.1038/nature06391.PubMedView ArticleGoogle Scholar
- Huang D, Wu W, Abrams SR, Cutler AJ: The relationship of drought-related gene expression in Arabidopsis thaliana to hormonal and environmental factors. J Exp Bot. 2008, 59: 2991-3007. 10.1093/jxb/ern155.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang X, Bernatavichute YV, Cokus S, Pellegrini M, Jacobsen SE: Genome-wide analysis of mono-, di- and trimethylation of histone H3 lysine 4 in Arabidopsis thaliana. Genome Biol. 2009, 10: R62-10.1186/gb-2009-10-6-r62.PubMedPubMed CentralView ArticleGoogle Scholar
- Lee CK, Shibata Y, Rao B, Strahl BD, Lieb JD: Evidence for nucleosome depletion at active regulatory regions genome-wide. Nat Genet. 2004, 36: 900-905. 10.1038/ng1400.PubMedView ArticleGoogle Scholar
- Oh S, Park S, van Nocker S: Genic and global functions for Paf1C in chromatin modification and gene expression in Arabidopsis. PLoS Genet. 2008, 4: e1000077-10.1371/journal.pgen.1000077.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang X, Clarenz O, Cokus S, Bernatavichute YV, Pellegrini M, Goodrich J, Jacobsen SE: Whole-genome analysis of histone H3 lysine 27 trimethylation in Arabidopsis. PLoS Biol. 2007, 5: e129-10.1371/journal.pbio.0050129.PubMedPubMed CentralView ArticleGoogle Scholar
- Kodama Y, Nagaya S, Shinmyo A, Kato K: Mapping and characterization of DNase I hypersensitive sites in Arabidopsis chromatin. Plant Cell Physiol. 2007, 48: 459-470. 10.1093/pcp/pcm017.PubMedView ArticleGoogle Scholar
- Yuan GC, Liu YJ, Dion MF, Slack MD, Wu LF, Altschuler SJ, Rando OJ: Genome-scale identification of nucleosome positions in S. cerevisiae. Science. 2005, 309: 626-630. 10.1126/science.1112178.PubMedView ArticleGoogle Scholar
- Lam FH, Steger DJ, O'Shea EK: Chromatin decouples promoter threshold from dynamic range. Nature. 2008, 453: 246-250. 10.1038/nature06867.PubMedPubMed CentralView ArticleGoogle Scholar
- Weiner A, Hughes A, Yassour M, Rando OJ, Friedman N: High-resolution nucleosome mapping reveals transcription-dependent promoter packaging. Genome Res. 2010, 20: 90-100. 10.1101/gr.098509.109.PubMedPubMed CentralView ArticleGoogle Scholar
- Mavrich TN, Jiang C, Ioshikhes IP, Li X, Venters BJ, Zanton SJ, Tomsho LP, Qi J, Glaser RL, Schuster SC, Gilmour DS, Albert I, Pugh BF: Nucleosome organization in the Drosophila genome. Nature. 2008, 453: 358-362. 10.1038/nature06929.PubMedPubMed CentralView ArticleGoogle Scholar
- Schones DE, Cui K, Cuddapah S, Roh TY, Barski A, Wang Z, Wei G, Zhao K: Dynamic regulation of nucleosome positioning in the human genome. Cell. 2008, 132: 887-898. 10.1016/j.cell.2008.02.022.PubMedView ArticleGoogle Scholar
- Chodavarapu RK, Feng S, Bernatavichute YV, Chen PY, Stroud H, Yu Y, Hetzel JA, Kuo F, Kim J, Cokus SJ, Casero D, Bernal M, Huijser P, Clark AT, Kramer U, Merchant SS, Zhang X, Jacobsen SE, Pellegrini M: Relationship between nucleosome positioning and DNA methylation. Nature. 2010, 466: 388-392. 10.1038/nature09147.PubMedPubMed CentralView ArticleGoogle Scholar
- Segal E, Fondufe-Mittendorf Y, Chen L, Thastrom A, Field Y, Moore IK, Wang JP, Widom J: A genomic code for nucleosome positioning. Nature. 2006, 442: 772-778. 10.1038/nature04979.PubMedPubMed CentralView ArticleGoogle Scholar
- Saleh A, Alvarez-Venegas R, Avramova Z: Dynamic and stable histone H3 methylation patterns at the Arabidopsis FLC and AP1 loci. Gene. 2008, 423: 43-47. 10.1016/j.gene.2008.06.022.PubMedView ArticleGoogle Scholar
- van Dijk K, Marley KE, Jeong BR, Xu J, Hesson J, Cerny RL, Waterborg JH, Cerutti H: Monomethyl histone H3 lysine 4 as an epigenetic mark for silenced euchromatin in Chlamydomonas. Plant Cell. 2005, 17: 2439-2453. 10.1105/tpc.105.034165.PubMedPubMed CentralView ArticleGoogle Scholar
- Heintzman ND, Stuart RK, Hon G, Fu Y, Ching CW, Hawkins RD, Barrera LO, Van Calcar S, Qu C, Ching KA, Wang W, Weng Z, Green RD, Crawford GE, Ren B: Distinct and predictive chromatin signatures of transcriptional promoters and enhancers in the human genome. Nat Genet. 2007, 39: 311-318. 10.1038/ng1966.PubMedView ArticleGoogle Scholar
- Alvarez-Venegas R, Avramova Z: Methylation patterns of histone H3 Lys 4, Lys 9 and Lys 27 in transcriptionally active and inactive Arabidopsis genes and in atx1 mutants. Nucleic Acids Res. 2005, 33: 5199-5207. 10.1093/nar/gki830.PubMedPubMed CentralView ArticleGoogle Scholar
- Pokholok DK, Harbison CT, Levine S, Cole M, Hannett NM, Lee TI, Bell GW, Walker K, Rolfe PA, Herbolsheimer E, Zeitlinger J, Lewitter F, Gifford DK, Young RA: Genome-wide map of nucleosome acetylation and methylation in yeast. Cell. 2005, 122: 517-527. 10.1016/j.cell.2005.06.026.PubMedView ArticleGoogle Scholar
- Li X, Wang X, He K, Ma Y, Su N, He H, Stolc V, Tongprasit W, Jin W, Jiang J, Terzaghi W, Li S, Deng XW: High-resolution mapping of epigenetic modifications of the rice genome uncovers interplay between dna methylation, histone methylation, and gene expression. Plant Cell. 2008, 20: 259-276. 10.1105/tpc.107.056879.PubMedPubMed CentralView ArticleGoogle Scholar
- Barski A, Cuddapah S, Cui K, Roh TY, Schones DE, Wang Z, Wei G, Chepelev I, Zhao K: High-resolution profiling of histone methylations in the human genome. Cell. 2007, 129: 823-837. 10.1016/j.cell.2007.05.009.PubMedView ArticleGoogle Scholar
- Guenther MG, Levine SS, Boyer LA, Jaenisch R, Young RA: A chromatin landmark and transcription initiation at most promoters in human cells. Cell. 2007, 130: 77-88. 10.1016/j.cell.2007.05.042.PubMedPubMed CentralView ArticleGoogle Scholar
- Ng DW, Chandrasekharan MB, Hall TC: Ordered Histone Modifications Are Associated with Transcriptional Poising and Activation of the phaseolin Promoter. Plant Cell. 2006, 18: 119-132. 10.1105/tpc.105.037010.PubMedPubMed CentralView ArticleGoogle Scholar
- Kouzarides T: Chromatin modifications and their function. Cell. 2007, 128: 693-705. 10.1016/j.cell.2007.02.005.PubMedView ArticleGoogle Scholar
- Shilatifard A: Molecular implementation and physiological roles for histone H3 lysine 4 (H3K4) methylation. Curr Opin Cell Biol. 2008, 20: 341-348. 10.1016/j.ceb.2008.03.019.PubMedPubMed CentralView ArticleGoogle Scholar
- Shukla A, Chaurasia P, Bhaumik SR: Histone methylation and ubiquitination with their cross-talk and roles in gene expression and stability. Cell Mol Life Sci. 2009, 66: 1419-1433. 10.1007/s00018-008-8605-1.PubMedView ArticleGoogle Scholar
- Ansari KI, Mandal SS: Mixed lineage leukemia: roles in gene expression, hormone signaling and mRNA processing. Febs J. 2010, 277: 1790-1804. 10.1111/j.1742-4658.2010.07606.x.PubMedView ArticleGoogle Scholar
- Ruthenburg AJ, Li H, Patel DJ, Allis CD: Multivalent engagement of chromatin modifications by linked binding modules. Nat Rev Mol Cell Biol. 2007, 8: 983-994. 10.1038/nrm2298.PubMedPubMed CentralView ArticleGoogle Scholar
- Ng HH, Robert F, Young RA, Struhl K: Targeted recruitment of Set1 histone methylase by elongating Pol II provides a localized mark and memory of recent transcriptional activity. Mol Cell. 2003, 11: 709-719. 10.1016/S1097-2765(03)00092-3.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.