Influence of histone modifications on mRNA abundance and structure

Abstract

Histones are frequently decorated with covalent modifications. These histone modifications are thought to be involved in various chromatin-dependent processes including transcription and splicing. To elucidate the relationship between histone modifications and these two processes, we derived models to predict the expression level of genes and the structure of transcribed mRNA from histone modification levels. We found that histone modification levels and gene expression are very well correlated. Moreover, we show that only a small number of histone modifications are necessary to accurately predict gene expression. We show that different sets of histone modifications are necessary to predict gene expression driven by high CpG content promoters (HCPs) or low CpG content promoters (LCPs). Quantitative models involving H3K4me3 and H3K79me1 are the most predictive of the expression levels in LCPs, whereas HCPs require H3K27ac and H4K20me1. We propose a preliminary “Histone Code of Transcription”, where H3K4me3 is involved in RNA polymerase II (PolII) recruitment and/or initiation, the combinatorial action of H3K27ac and H4K20me1 leads to the transition to elongation, and finally H3K79me1 and H4K20me1 signal the transition to an elongating PolII. The preliminary “Histone Code of Transcription” awaits confirmation by further experimental studies. We furthermore show that the connections between histone modifications and gene expression seem to be general, as we were able to predict gene expression levels of one cell type using a model trained on another one. We propose that our model could be further improved by including information about different mechanisms of regulation of mRNA stability and degradation, giving rise to more accurate predictions of gene expression levels. Using logistic models, we showed that levels of histone modifications, nucleosomes and PolII are predictive of the splicing outcomes of alternative exons, and that this result is not a consequence of experimental artifacts. Furthermore, we identified four histone modifications, namely H3K27me1, H3K36ac, H3K36me3 and H3K79me2, which consistently have a significant contribution to prediction accuracy of our models. This finding implies that they could be directly related to the splicing process, in agreement with recent analyses of the relationship of chromatin structure and the splicing process. We also established that histone modifications convey information about alternative splicing different from the one encoded in the DNA sequence of exons and surrounding regions, suggesting a possible interplay between these two mechanisms of splicing regulation. Finally, we confirmed the existence of functional coupling of transcription and splicing, by studying the dependence of structure of transcripts on their expression levels. The exact mechanisms behind these observations will have to be studied further.

Histonproteinen liegen häufig chemisch modifiziert vor. Diese Modifikationen sind an vielen chromatinabhängigen Prozessen beteiligt. Diese Arbeit untersucht den Zusammenhang zwischen Histonmodifikationen und Transkription bzw. Splicing anhand von Modellen, die den Expressionslevel von Genen bzw. die Struktur von mRNAs mit Hilfe der Histonmodifikationen vorhersagen. Unsere Ergebnisse zeigen, dass die Häufigkeit von Histonmodifikationen am Promotor und der Genexpressionslevel stark miteinander korreliert sind. Die Vorhersage der Genexpression hängt dabei nur von wenigen Histonmodifikationen ab. Dabei haben wir unterschiedliche Gruppen von Modifikationen identifiziert, die für eine gute Vorhersagequalität in Promotoren mit hohem bzw. niedrigen CpG Gehalt notwendig sind. Quantitative Modelle, die die Information von H4K4me3 und H3K79me1 beinhalten, haben die beste Qualität in Promotoren mit niedrigem CpG Gehalt, während Modelle, die H3K27ac und H4K20me1 verwenden, am Besten sind für Promotoren mit hohem CpG Gehalt. Basierend auf diesen Ergebnissen schlagen wir einen vorläufigen „Histoncode für die Transkription“ vor, in dem H3K4me3 an der Rekrutierung und/oder Initiation von RNA Polymerase II (Pol II) beteiligt ist, H3K27ac und H4K20me1 den Übergang zur Elongation ermöglicht und H3K79me1 und H4K20me1 den erfolgreichen Übergang anzeigt. Dieser vorläufige „Histoncode für die Transkription“ muss in zukünftigen experimentellen Studien kritisch überprüft werden. Die gefundenen Zusammenhänge zwischen Histonmodifikationen und Genexpression sind von allgemeiner Natur, da es möglich war, die Genexpression von Zellen mit Modellen vorherzusagen, die in einem anderen Zelltyp erstellt worden sind. Unsere Ergebnisse zeigen weiterhin, dass unsere Modelle durch die Berücksichtigung von Prozessen, die die mRNA Stabilität beeinflussen, weiter verbessert werden können. Unsere Untersuchungen zeigen, dass die Häufigkeit von Histonmodifikationen, Nukleosomen und Pol II verwendet werden können, um alternatives Splicing vorherzusagen. Wir haben vier Histonmodifikationen (H3K27me1, H3K36ac, H3K36me3 und H3K79me2) identifiziert, die signifikant zur Vorhersagequalität unserer Modelle beitragen. Dieses Ergebnis legt nahe, dass diese Modifikationen im direkten Zusammenhang mit dem Splicingprozess stehen könnten. Histonmodifikationen tragen Information über alternatives Splicing, die z.T. komplementär zu den Sequenzinformationen in Exons und den umgebenden Regionen sind, was auf ein Zusammenspiel von Histonmodifikations- und Sequenzabhängigen Prozessen in der Regulation von Splicing hindeutet. Unsere Ergbnisse zeigen eine Abhängigkeit zwischen der Struktur von Transkripten und deren Expressionslevel, was eine funkti onelle Kopplung zwischen Transkription und Splicing bestätigt. Die Mechanismen, die diesen Beobachtungen zu Grunde liegen, müssen in Zukunft weiter untersucht werden.