During the reprogramming of mouse embryonic fibroblasts (MEFs) to induced pluripotent stem cells, the activation of pluripotency genes such as NANOG occurs after the mesenchymal to epithelial transition. Here we report that both adult stem cells (neural stem cells) and differentiated cells (astrocytes) of the neural lineage can activate NANOG in the absence of cadherin expression during reprogramming. Gene expression analysis revealed that only the NANOG+E-cadherin+ populations expressed stabilization markers, had upregulated several cell cycle genes; and were transgene independent. Inhibition of DOT1L activity enhanced both the numbers of NANOG+ and NANOG+E-cadherin+ colonies in neural stem cells. Expressing SOX2 in MEFs prior to reprogramming did not alter the ratio of NANOG colonies that express E-cadherin. Taken together these results provide a unique pathway for reprogramming taken by cells of the neural lineage.

Long range regulatory interactions among distal enhancers and target genes are important for tissue-specific gene expression. Genome-scale identification of these interactions in a cell line-specific manner, especially using the fewest possible datasets, is a significant challenge. We develop a novel computational approach, Regulatory Interaction Prediction for Promoters and Long-range Enhancers (RIPPLE), that integrates published Chromosome Conformation Capture (3C) data sets with a minimal set of regulatory genomic data sets to predict enhancer-promoter interactions in a cell line-specific manner. Our results suggest that CTCF, RAD21, a general transcription factor (TBP) and activating chromatin marks are important determinants of enhancer-promoter interactions. To predict interactions in a new cell line and to generate genome-wide interaction maps, we develop an ensemble version of RIPPLE and apply it to generate interactions in five human cell lines. Computational validation of these predictions using existing ChIA-PET and Hi-C data sets showed that RIPPLE accurately predicts interactions among enhancers and promoters. Enhancer-promoter interactions tend to be organized into subnetworks representing coordinately regulated sets of genes that are enriched for specific biological processes and cis-regulatory elements. Overall, our work provides a systematic approach to predict and interpret enhancer-promoter interactions in a genome-wide cell-type specific manner using a few experimentally tractable measurements.

An open chromatin largely devoid of heterochromatin is a hallmark of stem cells. It remains unknown whether an open chromatin is necessary for the differentiation potential of stem cells, and which molecules are needed to maintain open chromatin. Here we show that the chromatin remodelling factor Chd1 is required to maintain the open chromatin of pluripotent mouse embryonic stem cells. Chd1 is a euchromatin protein that associates with the promoters of active genes, and downregulation of Chd1 leads to accumulation of heterochromatin. Chd1-deficient embryonic stem cells are no longer pluripotent, because they are incapable of giving rise to primitive endoderm and have a high propensity for neural differentiation. Furthermore, Chd1 is required for efficient reprogramming of fibroblasts to the pluripotent stem cell state. Our results indicate that Chd1 is essential for open chromatin and pluripotency of embryonic stem cells, and for somatic cell reprogramming to the pluripotent state.

Elucidating the mechanism of reprogramming is confounded by heterogeneity due to the low efficiency and differential kinetics of obtaining induced pluripotent stem cells (iPSCs) from somatic cells. Therefore, we increased the efficiency with a combination of epigenomic modifiers and signaling molecules and profiled the transcriptomes of individual reprogramming cells. Contrary to the established temporal order, somatic gene inactivation and upregulation of cell cycle, epithelial, and early pluripotency genes can be triggered independently such that any combination of these events can occur in single cells. Sustained co-expression of Epcam, Nanog, and Sox2 with other genes is required to progress toward iPSCs. Ehf, Phlda2, and translation initiation factor Eif4a1 play functional roles in robust iPSC generation. Using regulatory network analysis, we identify a critical role for signaling inhibition by 2i in repressing somatic expression and synergy between the epigenomic modifiers ascorbic acid and a Dot1L inhibitor for pluripotency gene activation.

Inhibiting the histone 3 lysine 79 (H3K79) methyltransferase, disruptor of telomeric silencing 1-like (DOT1L), increases the efficiency of reprogramming somatic cells to induced pluripotent stem cells (iPSCs). Here, we find that, despite the enrichment of H3K79 methylation on thousands of actively transcribed genes in somatic cells, DOT1L inhibition (DOT1Li) does not immediately cause the shutdown of the somatic transcriptional profile to enable transition to pluripotency. Contrary to the prevalent view, DOT1Li promotes iPSC generation beyond the mesenchymal to epithelial transition and even from already epithelial cell types. DOT1Li is most potent at the midpoint of reprogramming in part by repressing Nfix that persists at late stages of reprogramming. Importantly, regulation of single genes cannot substitute for DOT1Li, demonstrating that H3K79 methylation has pleiotropic effects in maintaining cell identity.

Changes in transcriptional regulatory networks can significantly alter cell fate. To gain insight into transcriptional dynamics, several studies have profiled bulk multi-omic data sets with parallel transcriptomic and epigenomic measurements at different stages of a developmental process. However, integrating these data to infer cell type-specific regulatory networks is a major challenge. We present dynamic regulatory module networks (DRMNs), a novel approach to infer cell type-specific cis-regulatory networks and their dynamics. DRMN integrates expression, chromatin state, and accessibility to predict cis-regulators of context-specific expression, where context can be cell type, developmental stage, or time point, and uses multitask learning to capture network dynamics across linearly and hierarchically related contexts. We applied DRMNs to study regulatory network dynamics in three developmental processes, each showing different temporal relationships and measuring a different combination of regulatory genomic data sets: cellular reprogramming, liver dedifferentiation, and forward differentiation. DRMN identified known and novel regulators driving cell type-specific expression patterns, showing its broad applicability to examine dynamics of gene regulatory networks from linearly and hierarchically related multi-omic data sets.

Reprogramming of somatic cells to induced pluripotent stem cells (iPSCs) represents a profound change in cell fate. Here, we show that combining ascorbic acid (AA) and 2i (MAP kinase and GSK inhibitors) increases the efficiency of reprogramming from fibroblasts and synergistically enhances conversion of partially reprogrammed intermediates to the iPSC state. AA and 2i induce differential transcriptional responses, each leading to the activation of specific pluripotency loci. A unique cohort of pluripotency genes including Esrrb require both stimuli for activation. Temporally, AA-dependent histone demethylase effects are important early, whereas Tet enzyme effects are required throughout the conversion. 2i function could partially be replaced by depletion of components of the epidermal growth factor (EGF) and insulin growth factor pathways, indicating that they act as barriers to reprogramming. Accordingly, reduction in the levels of the EGF receptor gene contributes to the activation of Esrrb. These results provide insight into the rewiring of the pluripotency network at the late stage of reprogramming.

The heterochromatin protein 1 (HP1) family is involved in various functions with maintenance of chromatin structure. During murine somatic cell reprogramming, we find that early depletion of HP1γ reduces the generation of induced pluripotent stem cells, while late depletion enhances the process, with a concomitant change from a centromeric to nucleoplasmic localization and elongation-associated histone H3.3 enrichment. Depletion of heterochromatin anchoring protein SENP7 increased reprogramming efficiency to a similar extent as HP1γ, indicating the importance of HP1γ release from chromatin for pluripotency acquisition. HP1γ interacted with OCT4 and DPPA4 in HP1α and HP1β knockouts and in H3K9 methylation depleted H3K9M embryonic stem cell (ESC) lines. HP1α and HP1γ complexes in ESCs differed in association with histones, the histone chaperone CAF1 complex, and specific components of chromatin-modifying complexes such as DPY30, implying distinct functional contributions. Taken together, our results reveal the complex contribution of the HP1 proteins to pluripotency.

S-adenosylmethionine (SAM) is the methyl-donor substrate for DNA and histone methyltransferases that regulate epigenetic states and subsequent gene expression. This metabolism-epigenome link sensitizes chromatin methylation to altered SAM abundance, yet the mechanisms that allow organisms to adapt and protect epigenetic information during life-experienced fluctuations in SAM availability are unknown. We identified a robust response to SAM depletion that is highlighted by preferential cytoplasmic and nuclear mono-methylation of H3 Lys 9 (H3K9) at the expense of broad losses in histone di- and tri-methylation. Under SAM-depleted conditions, H3K9 mono-methylation preserves heterochromatin stability and supports global epigenetic persistence upon metabolic recovery. This unique chromatin response was robust across the mouse lifespan and correlated with improved metabolic health, supporting a significant role for epigenetic adaptation to SAM depletion in vivo. Together, these studies provide evidence for an adaptive response that enables epigenetic persistence to metabolic stress.

Changes in chromatin state play important roles in cell fate transitions. Current computational approaches to analyze chromatin modifications across multiple cell types do not model how the cell types are related on a lineage or over time. To overcome this limitation, we developed a method called Chromatin Module INference on Trees (CMINT), a probabilistic clustering approach to systematically capture chromatin state dynamics across multiple cell types. Compared to existing approaches, CMINT can handle complex lineage topologies, capture higher quality clusters, and reliably detect chromatin transitions between cell types. We applied CMINT to gain novel insights in two complex processes: reprogramming to induced pluripotent stem cells (iPSCs) and hematopoiesis. In reprogramming, chromatin changes could occur without large gene expression changes, different combinations of activating marks were associated with specific reprogramming factors, there was an order of acquisition of chromatin marks at pluripotency loci, and multivalent states (comprising previously undetermined combinations of activating and repressive histone modifications) were enriched for CTCF. In the hematopoietic system, we defined critical decision points in the lineage tree, identified regulatory elements that were enriched in cell-type-specific regions, and found that the underlying chromatin state was achieved by specific erasure of preexisting chromatin marks in the precursor cell or by de novo assembly. Our method provides a systematic approach to model the dynamics of chromatin state to provide novel insights into the relationships among cell types in diverse cell-fate specification processes.

Immune-mediated destruction of insulin-producing β cells causes type 1 diabetes (T1D). However, how β cells participate in their own destruction during the disease process is poorly understood. Here, we report that modulating the unfolded protein response (UPR) in β cells of non-obese diabetic (NOD) mice by deleting the UPR sensor IRE1α prior to insulitis induced a transient dedifferentiation of β cells, resulting in substantially reduced islet immune cell infiltration and β cell apoptosis. Single-cell and whole-islet transcriptomics analyses of immature β cells revealed remarkably diminished expression of β cell autoantigens and MHC class I components, and upregulation of immune inhibitory markers. IRE1α-deficient mice exhibited significantly fewer cytotoxic CD8+ T cells in their pancreata, and adoptive transfer of their total T cells did not induce diabetes in Rag1-/- mice. Our results indicate that inducing β cell dedifferentiation, prior to insulitis, allows these cells to escape immune-mediated destruction and may be used as a novel preventive strategy for T1D in high-risk individuals.

Cell type-specific gene expression patterns are outputs of transcriptional gene regulatory networks (GRNs) that connect transcription factors and signaling proteins to target genes. Single-cell technologies such as single cell RNA-sequencing (scRNA-seq) and single cell Assay for Transposase-Accessible Chromatin using sequencing (scATAC-seq), can examine cell-type specific gene regulation at unprecedented detail. However, current approaches to infer cell type-specific GRNs are limited in their ability to integrate scRNA-seq and scATAC-seq measurements and to model network dynamics on a cell lineage. To address this challenge, we have developed single-cell Multi-Task Network Inference (scMTNI), a multi-task learning framework to infer the GRN for each cell type on a lineage from scRNA-seq and scATAC-seq data. Using simulated and real datasets, we show that scMTNI is a broadly applicable framework for linear and branching lineages that accurately infers GRN dynamics and identifies key regulators of fate transitions for diverse processes such as cellular reprogramming and differentiation.

The heterochromatin protein 1 (HP1) family members are canonical effectors and propagators of gene repression mediated by histone H3 lysine 9 (H3K9) methylation. HP1γ exhibits an increased interaction with active transcription elongation-associated factors in embryonic stem cells (ESCs) compared to somatic cells. However, whether this association has a functional consequence remains elusive. Here we find that genic HP1γ colocalizes and enhances enrichment of transcription elongation-associated H3K36me3 rather than H3K9me3. Unexpectedly, sustained H3K36me3 deposition is dependent on HP1γ. HP1γ-deleted ESCs display reduced H3K36me3 enrichment, concomitant with decreased expression at shared genes which function to maintain cellular homeostasis. Both the H3K9me3-binding chromodomain and histone binding ability of HP1γ are dispensable for maintaining H3K36me3 levels. Instead, the chromoshadow together with the hinge domain of HP1γ that confer protein and nucleic acid-binding ability are sufficient because they retain the ability to interact with NSD1, an H3K36 methyltransferase. HP1γ-deleted ESCs have a slower self-renewal rate and an impaired ability to differentiate towards cardiac mesoderm. Our findings reveal a requirement for HP1γ in faithful establishment of transcription elongation in ESCs, which regulates pluripotency.

Histone modifying enzymes play a central role in maintaining cell identity by establishing a conducive chromatin environment for lineage specific transcription factor activity. Pluripotent embryonic stem cell (ESC) identity is characterized by a lower abundance of gene repression associated histone modifications that enables rapid response to differentiation cues. The KDM3 family of histone demethylases removes the repressive histone H3 lysine 9 dimethylation (H3K9me2). Here we uncover a surprising role for the KDM3 proteins in the maintenance of the pluripotent state through post-transcriptional regulation. We find that KDM3A and KDM3B interact with RNA processing factors such as EFTUD2 and PRMT5. Acute selective degradation of the endogenous KDM3A and KDM3B proteins resulted in altered splicing independent of H3K9me2 status or catalytic activity. These splicing changes partially resemble the splicing pattern of the more blastocyst-like ground state of pluripotency and occurred in important chromatin and transcription factors such as Dnmt3b, Tbx3 and Tcf12. Our findings reveal non-canonical roles of histone demethylating enzymes in splicing to regulate cell identity.

Leucine-rich repeat containing 10 (LRRC10) is a cardiomyocyte-specific protein, but its role in cardiac biology is little understood. Recently Lrrc10 was identified as required for endogenous cardiac regeneration in zebrafish; however, whether LRRC10 plays a role in mammalian heart regeneration remains unclear. In this study, we demonstrate that Lrrc10-/- knockout mice exhibit a loss of the neonatal mouse regenerative response, marked by reduced cardiomyocyte cytokinesis and increased cardiomyocyte binucleation. Interestingly, LRRC10 deletion disrupts the regenerative transcriptional landscape of the regenerating neonatal mouse heart. Remarkably, cardiac overexpression of LRRC10 restores cardiomyocyte cytokinesis, increases cardiomyocyte mononucleation, and the cardiac regenerative capacity of Lrrc10-/- mice. Our results are consistent with a model in which LRRC10 is required for cardiomyocyte cytokinesis as well as regulation of the transcriptional landscape during mammalian heart regeneration.

Reprogramming of somatic cells into induced pluripotent stem cells (iPSCs) involves a marked reorganization of chromatin. To identify post-translational histone modifications that change in global abundance during this process, we have applied a quantitative mass-spectrometry-based approach. We found that iPSCs, compared with both the starting fibroblasts and a late reprogramming intermediate (pre-iPSCs), are enriched for histone modifications associated with active chromatin, and depleted for marks of transcriptional elongation and a subset of repressive modifications including H3K9me2/me3. Dissecting the contribution of H3K9 methylation to reprogramming, we show that the H3K9 methyltransferases Ehmt1, Ehmt2 and Setdb1 regulate global H3K9me2/me3 levels and that their depletion increases iPSC formation from both fibroblasts and pre-iPSCs. Similarly, we find that inhibition of heterochromatin protein-1γ (Cbx3), a protein known to recognize H3K9 methylation, enhances reprogramming. Genome-wide location analysis revealed that Cbx3 predominantly binds active genes in both pre-iPSCs and pluripotent cells but with a strikingly different distribution: in pre-iPSCs, but not in embryonic stem cells, Cbx3 associates with active transcriptional start sites, suggesting a developmentally regulated role for Cbx3 in transcriptional activation. Despite largely non-overlapping functions and the predominant association of Cbx3 with active transcription, the H3K9 methyltransferases and Cbx3 both inhibit reprogramming by repressing the pluripotency factor Nanog. Together, our findings demonstrate that Cbx3 and H3K9 methylation restrict late reprogramming events, and suggest that a marked change in global chromatin character constitutes an epigenetic roadblock for reprogramming.

Much has been learned about the mechanisms of action of pluripotency factors Oct4 and Sox2. However, as with other regulators of cell identity, little is known about the impact of disrupting their binding motifs in a native environment or the characteristics of genes they regulate. By quantitatively examining dynamic ranges of gene expression instead of focusing on conventional measures of differential expression, we found that Oct4 and Sox2 enhancer binding is strongly enriched near genes subject to large dynamic ranges of expression among cell types, with binding sites near these genes usually within superenhancers. Mutagenesis of representative Oct4:Sox2 motifs near such active, dynamically regulated genes revealed critical roles in transcriptional activation during reprogramming, with more limited roles in transcriptional maintenance in the pluripotent state. Furthermore, representative motifs near silent genes were critical for establishing but not maintaining the fully silent state, while genes whose transcript levels varied by smaller magnitudes among cell types were unaffected by nearby Oct4:Sox2 motifs. These results suggest that Oct4 and Sox2 directly establish both active and silent transcriptional states in pluripotent cells at a large number of genes subject to dynamic regulation during mammalian development, but are less important than expected for maintaining transcriptional states.