Abstract
Direct conversion of cardiac fibroblast into induced cardiomyocytes (iCMs) by forced expression of cardiac transcription factors, such as Mef2c, Gata4, and Tbx5 (MGT), holds great promise for regenerative medicine. The process of cardiac reprogramming consists of waves of transcriptome remodelling events. However, how this transcriptome remodelling is driven by the upstream chromatin landscape alteration is still unclear. In this study, we performed single-cell ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) on early reprogramming iCMs given the known epigenetic changes as early as day 3. This approach unveiled networks of transcription factors (TFs) involved in the early shift of chromatin accessibility during cardiac reprogramming. Combining our analysis with functional assays, we identified Smad3 to be a bimodal TF in cardiac reprogramming, a barrier in the initiation of reprogramming and a facilitator during the intermediate stage of reprogramming. Moreover, integrative analysis of scATAC-seq with scRNA-seq data led to the identification of active TFs important for iCM conversion. Finally, we discovered a global rewiring of cis-regulatory interactions of cardiac genes along the reprogramming trajectory. Collectively, our scATAC-seq study and the integrative analysis with scRNA-seq data provided valuable resources to understand the epigenomic heterogeneity and its alteration in relation to transcription changes during early stage of cardiac reprogramming.