Abstract
Direct reprogramming using defined sets of transcription factors (TFs) is a recent strategy for generating induced hepatocytes (iHeps) from fibroblasts for use in regenerative medicine and drug development. Comprehensive studies detailing the regulatory role of TFs during this reprogramming process could help increase its efficiency. This study aimed to find the TFs with the greatest influences on the generation of iHeps from fibroblasts, and to further understand their roles in the regulation of the gene expression program. Here, we used systems biology approaches to analyze high quality expression data sets in combination with TF-binding sites data and protein-protein interactions data during the direct reprogramming of fibroblasts to iHeps. Our results revealed two main patterns for differentially expressed genes (DEGs): up-regulated genes were categorized as hepatic-specific pattern, and down-regulated genes were categorized as mesoderm- and fibroblast-specific pattern. Interestingly, hepatic-specific genes co-expressed and were regulated by hepatic-specific TFs, specifically Hnf4a and Foxa2. Conversely, the mesoderm- and fibroblast-specific pattern was mainly silenced by polycomb repressive complex 2 (PRC2) members, including Suz12, Mtf2, Ezh2, and Jarid2. Independent analysis of both the gene and core regulatory network of DE-TFs showed significant roles for Hnf4a, Foxa2, and PRC2 members in the regulation of the gene expression program and in biological processes during the direct conversion process. Altogether, using systems biology approaches, we clarified the role of Hnf4a and Foxa2 as hepatic-specific TFs, and for the first time, introduced the PRC2 complex as the main regulator that favors the direct reprogramming process in cooperation with hepatic-specific factors.