Systems Biology in Heterogenous Tissues: Integrating Multiple *Omics Datasets to Understand Hematopoietic Differentiation

Abstract

Motivation: Systems biology integrates expression, methylation, transcription factor binding and histone modification profiles with other physiological characteristics of a specific organ. Repositories that provide the required data, like ENCODE, generally work on a high level and do not take the heterogeneity of cell types within an organ into consideration. The hematopoietic system allows the characterization and study of each cell type involved in the generation of blood cells from bone marrow stem cells and thus provides a good foundation for systems biology studies. Here we compare RNA expression, DNA methylation, chromatin accessibility, DNA binding proteins and histone modification profiles in seven different hematopoietic populations using a Bayesian non-parametric hierarchical latent-class mixed-effect model known as IDEAS to characterize epigenetic changes associated with hematopoietic differentiation. Unlike other existing approaches IDEAS considers various cell types of a biological systems in concert instead of disjointly. Results: Using the VISION database and the IDEAS toolkit we provide insights into the transcriptional, epigenetic and regulatory programs governing the hematopoietic differentiation process. The characterization of the different hematopoietic components and their interactions provide the foundations for a systems biology model of hematopoiesis. Previous hematopoietic epigenome segmentation studies have focused on histone modifications, chromatin accessibility and DNA binding protein profiles. DNA methylation has been shown to vary markedly in hematopoietic populations. Inclusion of DNA methylation in these segmentation studies increased the original 36-state model of regulatory interactions to 41 states. These new DNA methylation-related states were associated with repressive marks, active RNA transcription, and a novel state regulated by DNA methylation alone. Imputing epigenetic models on inputs systematically perturbed for hematopoietic populations resulted in models of varying degrees of overlap, which were quantified and set in context with underlying biological processes. Conclusion: Our data show that methylation has a strong impact on functional genomic modeling and can be used to discern cell type specific epigenetic regulatory behavior by leveraging imputation for missing cell type data.