A unidirectional histone code in bidirectional promoters across cell types

Recent studies indicate that thousands of genes are expressed from bidirectional promoters (BDPs). Gene regulation at BDPs is poorly understood, in particular how the cell is able to regulate them differently. Here we investigate the effect of histone Modifications in BDP regulation. In this study, we model the gene activity using different gene expression assays, such as RNA-Seq, GRO-cap, and CAGE around the BDP transcription start sites (TSSs) using different histone modifications in various cell types. We develop a new statistical approach that links histone modifications to gene expression at BDPs. It improves over previous methods, because it is able to capture spatial dependencies of histone modifications along a promoter and leads to more interpretable results. We predict a general histone code that is independent of transcript orientation, cell type, and promoter configuration. The histone code at BDPs reveals that promoters are regulated unidirectionally, such that the majority of histone marks associated with gene expression occur downstream of the gene's TSS. By contrasting associations of histone marks with steady-state levels of capped RNAs in CAGE and nascent capped RNAs in GRO-cap, we show which histone marks have a preference for initiating polymerases and actively elongating polymerases. Using single-cell data we show that the bidirectional histone signal of activating marks is often due averaging of a heterogeneous cell population. Our results have implications for other experiments where the relationship between histone modification and gene initiation or gene elongation is investigated on a genome-wide scale.EGA study EGAS00001001656