Alterations of lysine modifications on the histone H3 N-tail under drought stress conditions in Arabidopsis thaliana

Post-translational modification of histone N-tails affects eukaryotic gene activity. In Arabidopsis, the histone modification level correlates with gene activation and repression in vernalization and flowering processes, but there is little information on changes in histone modification status and nucleosome structure under abiotic stresses. We determined the temporal and spatial changes in nucleosome occupancy and levels of H3K4me3, H3K9ac, H3K14ac, H3K23ac and H3K27ac in the histone H3 N-tail on the regions of four Arabidopsis drought stress-inducible genes, RD29A, RD29B, RD20 and At2g20880 [corrected], under drought stress conditions by chromatin immunoprecipitation analysis. We found two types of regulatory mechanisms of nucleosome occupancy function in the drought stress response. For RD29A and RD29B genes, nucleosome occupancy of promoter regions is low compared with that of coding regions, and no notable nucleosome loss occurs under drought stress. In contrast, nucleosome density is gradually decreased in response to drought stress on RD20 and At2g20880 [corrected] genes. Enrichments of H3K4me3 and H3K9ac correlate with gene activation in response to drought stress in all four genes. Interestingly, establishment of H3K4me3 occurs after accumulation of RNAPII on the coding regions of RD29A and At2g20880 [corrected]. Enrichment of H3K23ac and H3K27ac occurs in response to drought stress on the coding regions of RD29B, RD20 and At2g20880 [corrected], but not on the coding region of At2g20880 [corrected]. Our results indicate that histone modifications on the H3 N-tail are altered with gene activation on the coding regions of drought stress-responsive genes under drought stress conditions and that several patterns of nucleosome changes function in the drought stress response.     

Hdac3, Setdb1, and Kap1 mark H3K9me3/H3K14ac bivalent regions in young and aged liver

Post‐translational modifications of histone tails play a crucial role in gene regulation. Here, we performed chromatin profiling by quantitative targeted mass spectrometry to assess all possible modifications of the core histones. We identified a bivalent combination, a dually marked H3K9me3/H3K14ac modification in the liver, that is significantly decreased in old hepatocytes. Subsequent sequential ChIP‐Seq identified dually marked single nucleosome regions, with reduced number of sites and decreased signal in old livers, confirming mass spectrometry results. We detected H3K9me3 and H3K14ac bulk ChIP‐Seq signal in reChIP nucleosome regions, suggesting a correlation between H3K9me3/H3K14ac bulk bivalent genomic regions and dually marked single nucleosomes. Histone H3K9 deacetylase Hdac3, as well as H3K9 methyltransferase Setdb1, found in complex Kap1, occupied both bulk and single nucleosome bivalent regions in both young and old livers, correlating to presence of H3K9me3. Expression of genes associated with bivalent regions in young liver, including those regulating cholesterol secretion and triglyceride synthesis, is upregulated in old liver once the bivalency is lost. Hence, H3K9me3/H3K14ac dually marked regions define a poised inactive state that is resolved with loss of one or both of the chromatin marks, which subsequently leads to change in gene expression.     

Analysis of biological features associated with meiotic recombination hot and cold spots in Saccharomyces cerevisiae

Meiotic recombination is not distributed uniformly throughout the genome. There are regions of high and low recombination rates called hot and cold spots, respectively. The recombination rate parallels the frequency of DNA double-strand breaks (DSBs) that initiate meiotic recombination. The aim is to identify biological features associated with DSB frequency. We constructed vectors representing various chromatin and sequence-based features for 1179 DSB hot spots and 1028 DSB cold spots. Using a feature selection approach, we have identified five features that distinguish hot from cold spots in Saccharomyces cerevisiae with high accuracy, namely the histone marks H3K4me3, H3K14ac, H3K36me3, and H3K79me3; and GC content. Previous studies have associated H3K4me3, H3K36me3, and GC content with areas of mitotic recombination. H3K14ac and H3K79me3 are novel predictions and thus represent good candidates for further experimental study. We also show nucleosome occupancy maps produced using next generation sequencing exhibit a bias at DSB hot spots and this bias is strong enough to obscure biologically relevant information. A computational approach using feature selection can productively be used to identify promising biological associations. H3K14ac and H3K79me3 are novel predictions of chromatin marks associated with meiotic DSBs. Next generation sequencing can exhibit a bias that is strong enough to lead to incorrect conclusions. Care must be taken when interpreting high throughput sequencing data where systematic biases have been documented.     

Genome-Wide Identification of Chromatin Transitional Regions Reveals Diverse Mechanisms Defining the Boundary of Facultative Heterochromatin

Due to the self-propagating nature of the heterochromatic modification H3K27me3, chromatin barrier activities are required to demarcate the boundary and prevent it from encroaching into euchromatic regions. Studies in Drosophila and vertebrate systems have revealed several important chromatin barrier elements and their respective binding factors. However, epigenomic data indicate that the binding of these factors are not exclusive to chromatin boundaries. To gain a comprehensive understanding of facultative heterochromatin boundaries, we developed a two-tiered method to identify the Chromatin Transitional Region (CTR), i.e. the nucleosomal region that shows the greatest transition rate of the H3K27me3 modification as revealed by ChIP-Seq. This approach was applied to identify CTRs in Drosophila S2 cells and human HeLa cells. Although many insulator proteins have been characterized in Drosophila, less than half of the CTRs in S2 cells are associated with known insulator proteins, indicating unknown mechanisms remain to be characterized. Our analysis also revealed that the peak binding of insulator proteins are usually 1–2 nucleosomes away from the CTR. Comparison of CTR-associated insulator protein binding sites vs. those in heterochromatic region revealed that boundary-associated binding sites are distinctively flanked by nucleosome destabilizing sequences, which correlates with significant decreased nucleosome density and increased binding intensities of co-factors. Interestingly, several subgroups of boundaries have enhanced H3.3 incorporation but reduced nucleosome turnover rate. Our genome-wide study reveals that diverse mechanisms are employed to define the boundaries of facultative heterochromatin. In both Drosophila and mammalian systems, only a small fraction of insulator protein binding sites co-localize with H3K27me3 boundaries. However, boundary-associated insulator binding sites are distinctively flanked by nucleosome destabilizing sequences, which correlates with significantly decreased nucleosome density and increased binding of co-factors.