An Unbiased Approach to Identifying Cellular Reprogramming-Inducible Enhancers [RNA-seq]

Cellular reprogramming of somatic cells towards induced Pluripotent Stem Cells (iPSCs) is achieved through the over-expression of the transcription factors Oct4, Sox2, Klf4 and c-Myc (OSKM). Due to the high complexity of the process, the vast number of the involved molecules and its stochastic and inefficient nature, the molecular background of reprogramming is only partially revealed. The aim of this study is the identification of genomic regulatory elements which are induced during cellular reprogramming, the examination of their action and their association with the expression of important genes. First, we constructed a combined ChIP-seq OSKM dataset using experiments on MEFs (Mouse Embryonic Fibroblasts) from multiple published studies, in order to examine the OSKM binding upon the genome, irrespective of minor discrepancies originating from the different reprogramming systems and protocols used. We also performed detailed RNA-seq and ATAC-seq experiments in order to monitor the functional results of the OSKM genomic binding. Initially, we focused on the ESC-related sites (Embryonic Stem Cells, natural equivalent of iPSCs) and found out that almost 1/3 of them are bound by OSKM as early as at day 1, or since MEFs. Our analysis highlighted the significant role of Myc in reprogramming, which is often overlooked. We saw that, along with Klf4, they “mark” positions where Oct4 and/or Sox2 will bind after reprogramming initiation. Furthermore, we discovered that various genes are regulated by multiple genomic elements, where OSKM bind at different timepoints in reprogramming. For example, many genes involved in the MET pathway (Mesenchymal to Epithelial Transition) are located near both ESC-, MEF- and Transient-OSKM sites, resulting in stably high expression levels which temporarily increase during the intermediate timepoints. Finally, we observed that OSKM often bind on the chromatin in a dynamic fashion. This “on-and-off” association of OSKM with the genome is related to the need of chromatin reorganization to achieve pluripotency. It is logical to assume that part of the inherent stochasticity of cellular reprogramming is based on the low percentage of cells that will achieve to reproduce this dynamic OSKM binding to regulate efficiently their transcriptional programs and to finally transform to iPSCs. Overall design: We examined the transcriptome of MEFs undergoing reprogramming at different time-points during the process (RNA-seq).