Flowchart illustrating the experimental pipeline presented in this study. ChIP/DNase-seq data were used to construct regression models of mRNA transcript abundance for a set of genes. The prediction accuracy of each model was evaluated relative to RNA-sequencing data. By constructing groups of genes categorised by biological process and applying the above methodology, it was possible to identify heterogeneity in the relative predictive power of TFs and HMs. These groups were later analysed for enrichment for housekeeping genes.