Histone modification profiles characterize function-specific gene regulation

Whole-genome landscape of H3K4me3, H3K36me3 and H3K9ac and their association with gene expression during Paulownia witches' broom disease infection and recovery processes

Histone methylation and acetylation participate in the modulation of gene expression. Here, chromatin immunoprecipitation sequencing (ChIP-Seq) was used to determine genome-wide patterns of three histone modifications, H3K4me3, H3K36me3, and H3K9ac (associated with actively expressed genes) and their associations with gene expression in Paulownia fortunei following phytoplasma infection and recovery from Paulownia witches' broom (PaWB) disease after methyl methane sulfonate treatment. The three histone marks were preferentially deposited in genic regions, especially downstream of transcription start sites, and were highly concurrent with gene expression. Genes with all three histone marks exhibited the highest expression levels. Based on the comparison scheme, we detected 365, 2244, and 752 PaWB-associated genes with H3K4me3, H3K36me3, and H3K9ac marks, separately. KEGG pathway analysis showed that these genes were involved in plant-pathogen interaction, plant hormone signal transduction, and starch and sucrose metabolism. A small proportion of differentially modified genes showed changes in expression in response to phytoplasma infection, including genes involved in calcium ion signal transduction, abscisic acid signal transduction, and ethylene biosynthesis. This comprehensive analysis of genome-wide histone modifications and gene expression in Paulownia following phytoplasma infection provides new insights into the epigenetic responses to phytoplasma infection and will be useful for further studies on epigenetic regulation mechanisms in plants under biotic stress.

Current evidence for a role of epigenetic mechanisms in response to ionizing radiation in an ecotoxicological context

The issue of potential long-term or hereditary effects for both humans and wildlife exposed to low doses (or dose rates) of ionising radiation is a major concern. Chronic exposure to ionising radiation, defined as an exposure over a large fraction of the organism's lifespan or even over several generations, can possibly have consequences in the progeny. Recent work has begun to show that epigenetics plays an important role in adaptation of organisms challenged to environmental stimulae. Changes to so-called epigenetic marks such as histone modifications, DNA methylation and non-coding RNAs result in altered transcriptomes and proteomes, without directly changing the DNA sequence. Moreover, some of these environmentally-induced epigenetic changes tend to persist over generations, and thus, epigenetic modifications are regarded as the conduits for environmental influence on the genome. Here, we review the current knowledge of possible involvement of epigenetics in the cascade of responses resulting from environmental exposure to ionising radiation. In addition, from a comparison of lab and field obtained data, we investigate evidence on radiation-induced changes in the epigenome and in particular the total or locus specific levels of DNA methylation. The challenges for future research and possible use of changes as an early warning (biomarker) of radiosensitivity and individual exposure is discussed. Such a biomarker could be used to detect and better understand the mechanisms of toxic action and inter/intra-species susceptibility to radiation within an environmental risk assessment and management context.

Genetic Contexts Characterize Dynamic Histone Modification Patterns Among Cell Types

Histone modifications play critical roles in mammalian development, regulating chromatin structure and gene expression. Dynamic histone modifications among cell types have been shown to associate with changes in mammalian development. However, how to quantitatively measure the histone modification alterations and how histone modifications vary across cell types under different genetic contexts remain largely unexplored and whether these changes are related to the primary DNA sequence remains limited. Here, we employed an entropy-based method to measure histone modification alterations in six definite genomic regions across five cell types and identified lineage-specific histone modification genes. We observed that histone modification alterations prefer to enrich in 5'-UTR exons, and also in 3'-UTR exons and its downstream. Then we built a model to predict the histone modification patterns from the primary DNA sequence. We found that the frequencies of k-mer sequence compositions are predictive of histone modification patterns, suggesting that the primary DNA sequence correlated with the histone modification alterations among cell types. Additionally, the lineage-specific histone modification genes display a higher conservation and lower GC-content. Together, we performed a systematic analysis for histone modification alterations and demonstrated how to identify genomic region-specific elements of epigenetic and genetic regulation and histone modification patterns across different cell types.

Role of RbBP5 and H3K4me3 in the vicinity of Snail transcription start site during epithelial-mesenchymal transition in prostate cancer cell

EMT (epithelial-mesenchymal transition) occurs in a wide range of tumor types, and has been shown to be crucial for metastasis. Epigenetic modifications of histones contribute to chromatin structure and result in the alterations in gene expression. Tri-methylation of histone H3 lysine 4 (H3K4me3) is associated with the promoters of actively transcribed genes and can serve as a transcriptional on/off switch. RbBP5 is a component of the COMPASS/ -like complex, which catalyzes H3K4me3 formation. In this study, we found that in the process of TGF-Beta1 induced EMT in the prostate cancer cell line DU145, H3K4me3 enrichment and RbBP5 binding increased in the vicinity of Snail (SNAI1) transcription start site. Knocking-down of RbBP5 notably decreased Snail expression and EMT. Recruitment of RbBP5 and formation of H3K4me3 at Snail TSS during EMT depend on binding of SMAD2/3 and CBP at Snail TSS. This study links the SMAD2/3 signal with Snail transcription via a histone modification - H3K4me3. Furthermore, our research also demonstrates that RbBP5 and even WRAD may be a promising therapeutic candidates in treating prostate cancer metastasis, and that DU145 cells maintain their incomplete mesenchymal state in an auto/paracrine manner.

Curcumin Derivative Epigenetically Reactivates Nrf2 Antioxidative Stress Signaling in Mouse Prostate Cancer TRAMP C1 Cells

The carcinogenesis of prostate cancer (PCa) in TRAMP model is highly correlated with hypermethylation in the promoter region of Nrf2 and the accompanying reduced transcription of Nrf2 and its regulated detoxifying genes. We aimed to investigate the effects of (3E,5E)-3,5-bis-(3,4,5-trimethoxybenzylidene)-tetrahydro-thiopyran-4-one (F10) and (3E,5E)-3,5-bis-(3,4,5-trimethoxy-benzylidene)-tetrahydropyran-4-one (E10), two synthetic curcumin derivatives, on restoring Nrf2 activity in TRAMP C1 cells. HepG2-C8 cells transfected with an antioxidant-response element (ARE)-luciferase vector were treated with F10, E10, curcumin, and sulforaphane (SFN) to compare their effects on Nrf2-ARE pathways. We performed real-time quantitative PCR and Western blotting to investigate the effects of F10 and E10 on Nrf2, correlated phase II detoxification genes. We also measured expression and activity of DNMTand HDAC enzymes. Enrichment of H3K27me3 on the promoter region of Nrf2 was explored with a chromatin immunoprecipitation (ChIP) assay. Methylation of the CpG region in Nrf2 promoter was doubly examined by bisulfite genomic sequencing (BGS) and methylation DNA immunoprecipitation (MeDIP). Compared with curcumin and SFN, F10 is more potent in activating Nrf2-ARE pathways. Both F10 and E10 enhanced level of Nrf2 and the correlated phase II detoxifying genes. BGS and MeDIP assays indicated that F10 but not E10 hypomethylated the Nrf2 promoter. F10 also downregulated the protein level of DNMT1, DNMT3a, DNMT3b, HDAC1, HDAC4, and HDAC7 and the activity of DNMTs and HDACs. F10 but not E10 effectively reduced the accumulation of H3k27me3 on the promoter of Nrf2. F10 and E10 can activate the Nrf2-ARE pathway and increase the level of Nrf2 and correlated phase II detoxification genes. The reactivation effect on Nrf2 by F10 in TRAMP C1 may come from demethylation, decrease of HDACs, and inhibition of H3k27me3 accumulation.

In silico prediction and characterization of protein post-translational modifications

This review outlines the computational approaches and procedures for predicting post translational modification (PTM)-induced changes in protein conformation and their influence on protein function(s), the latter being assessed as differential affinity in interaction with either low (ligands for receptors or transporters, substrates for enzymes) or high molecular mass molecules (proteins or nucleic acids in supramolecular assemblies). The scope for an in silico approach is discussed against a summary of the in vitro evidence on the structural and functional outcome of protein PTM.

Association analysis between the distributions of histone modifications and gene expression in the human embryonic stem cell

It is well known that histone modifications are associated with gene expression. In order to further study this relationship, 16 kinds of Chip-seq histone modification data and mRNA-seq data of the human embryonic stem cell H1 are chosen. The distributions of histone modifications in the regions flanking transcription start sites (TSSs) for highly expressed and lowly expressed genes are computed, respectively. And four types of distributions of histone modifications in regions flanking TSSs and the spatial patterning of the correlations between histone modifications and gene expression are detected. Our results suggest that the correlations between the regions overlapped by peaks are higher than the non-overlapped ones for each histone modification. In addition, to obtain the effect of the cooperative action of histone modification on gene expression, five histone modification clusters are found in highly expressed and lowly expressed genes, histone modification and gene expression interaction network is constructed. To further explore which region is the main target region for the specific histone modification, the human genes are divided into five functional regions. The results indicate that histone modifications are mostly located in the promoters of highly expressed genes versus the exons of lowly expressed genes, and exons have a smaller range of normalized tag counts than other gene elements in the two groups of genes. Finally, the type specificity and regional bias of histone modifications for 11 key transcription factor genes regulating the stem cell renewal are analyzed.

MRHMMs: multivariate regression hidden Markov models and the variantS

SUMMARY

Hidden Markov models (HMMs) are flexible and widely used in scientific studies. Particularly in genomics and genetics, there are multiple distinct regimes in the genome within each of which the relationships among multivariate features are distinct. Examples include differential gene regulation depending on gene functions and experimental conditions, and varying combinatorial patterns of multiple transcription factors. We developed a software package called MRHMMs (Multivariate Regression Hidden Markov Models and the variantS) that accommodates a variety of HMMs that can be flexibly applied to many biological studies and beyond. MRHMMs supplements existing HMM software packages in two aspects. First, MRHMMs provides a diverse set of emission probability structures, including mixture of multivariate normal distributions and (logistic) regression models. Second, MRHMMs is computationally efficient for analyzing large data-sets generated in current genome-wide studies. Especially, the software is written in C for the speed advantage and further amenable to implement alternative models to meet users' own purposes.

AVAILABILITY AND IMPLEMENTATION

http://sourceforge.net/projects/mrhmms/

Regression hidden Markov modeling reveals heterogeneous gene expression regulation: a case study in mouse embryonic stem cells

BACKGROUND

Studies have shown the strong association between histone modification levels and gene expression levels. The detailed relationships between the two can vary substantially due to differential regulation, and hence a simple regression model may not be adequate. We apply a regression hidden Markov model (regHMM) to further investigate the potential multiple relationships between genes and histone methylation levels in mouse embryonic stem cells.

RESULTS

Seven histone methylation levels are used in the study. Averaged histone modifications over non-overlapping 200 bp windows on the range transcription starting site (TSS) ± 1 Kb are used as predictors, and in total 70 explanatory variables are generated. Based on regHMM results, genes segregated into two groups, referred to as State 1 and State 2, have distinct association strengths. Genes in State 1 are better explained by histone methylation levels with R(2)=.72 while those in State 2 have weaker association strength with R(2)=.38. The regression coefficients in the two states are not very different in magnitude except in the intercept,.25 and 1.15 for State 1 and State 2, respectively. We found specific GO categories that may be attributed to the different relationships. The GO categories more frequently observed in State 2 match those of housekeeping genes, such as cytoplasm, nucleus, and protein binding. In addition, the housekeeping gene expression levels are significantly less explained by histone methylation in mouse embryonic stem cells, which is consistent with the constitutive expression patterns that would be expected.

CONCLUSION

Gene expression levels are not universally affected by histone methylation levels, and the relationships between the two differ by the gene functions. The expression levels of the genes that perform the most common housekeeping genes' GO categories are less strongly associated with histone methylation levels. We suspect that additional biological factors may also be strongly associated with the gene expression levels in State 2. We discover that the effect of the presence of CpG island in TSS ± 1 Kb is larger in State 2.

LinkNMF: identification of histone modification modules in the human genome using nonnegative matrix factorization

Histone modifications are ubiquitous processes involved in various cellular mechanisms. Systemic analysis of multiple chromatin modifications has been used to characterize various chromatin states associated with functional DNA elements, gene expression, and specific biological functions. However, identification of modular modification patterns is still required to understand the functional associations between histone modification patterns and specific chromatin/DNA binding factors. To recognize modular modification patterns, we developed a novel algorithm that combines nonnegative matrix factorization (NMF) and a clique-detection algorithm. We applied it, called LinkNMF, to generate a comprehensive modification map in human CD4+ T cell promoter regions. Initially, we identified 11 modules not recognized by conventional approaches. The modules were grouped into two major classes: gene activation and repression. We found that genes targeted by each module were enriched with distinguishable biological functions, suggesting that each modular pattern plays a unique functional role. To explain the formation of modular patterns, we investigated the module-specific binding patterns of chromatin regulators. Application of LinkNMF to histone modification maps of diverse cells and developmental stages will be helpful for understanding how histone modifications regulate gene expression. The algorithm is available on our website at biodb.kaist.ac.kr/LinkNMF.