Oct4:Sox2 binding is essential for establishing but not maintaining active and silent states of dynamically regulated genes in pluripotent cells.

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.

Methyl-Metabolite Depletion Elicits Adaptive Responses to Support Heterochromatin Stability and Epigenetic Persistence.

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.

Dynamic regulatory module networks for inference of cell type-specific transcriptional networks.

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.

Role of the murine reprogramming factors in the induction of pluripotency.

Induced pluripotent stem (iPS) cells can be obtained from fibroblasts upon expression of Oct4, Sox2, Klf4, and c-Myc. To understand how these factors induce pluripotency, we carried out genome-wide analyses of their promoter binding and expression in iPS and partially reprogrammed cells. We find that target genes of the four factors strongly overlap in iPS and embryonic stem (ES) cells. In partially reprogrammed cells, many genes co-occupied by c-Myc and any of the other three factors already show an ES cell-like binding and expression pattern. In contrast, genes that are specifically co-bound by Oct4, Sox2, and Klf4 in ES cells and encode pluripotency regulators severely lack binding and transcriptional activation. Among the four factors, c-Myc promotes the most ES cell-like transcription pattern when expressed individually in fibroblasts. These data uncover temporal and separable contributions of the four factors during the reprogramming process and indicate that ectopic c-Myc predominantly acts before pluripotency regulators are activated.

Coordinated removal of repressive epigenetic modifications during induced reversal of cell identity.

Epigenetic modifications operate in concert to maintain cell identity, yet how these interconnected networks suppress alternative cell fates remains unknown. Here, we uncover a link between the removal of repressive histone H3K9 methylation and DNA methylation during the reprogramming of somatic cells to pluripotency. The H3K9me2 demethylase, Kdm3b, transcriptionally controls DNA hydroxymethylase Tet1 expression. Unexpectedly, in the absence of Kdm3b, loci that must be DNA demethylated are trapped in an intermediate hydroxymethylated (5hmC) state and do not resolve to unmethylated cytosine. Ectopic 5hmC trapping precludes the chromatin association of master pluripotency factor, POU5F1, and pluripotent gene activation. Increased Tet1 expression is important for the later intermediates of the reprogramming process. Taken together, coordinated removal of distinct chromatin modifications appears to be an important mechanism for altering cell identity.

HP1γ regulates H3K36 methylation and pluripotency in embryonic stem cells.

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.

The role of α-ketoglutarate-dependent proteins in pluripotency acquisition and maintenance.

α-Ketoglutarate is an important metabolic intermediate that acts as a cofactor for several chromatin-modifying enzymes, including histone demethylases and the Tet family of enzymes that are involved in DNA demethylation. In this review, we focus on the function and genomic localization of these α-ketoglutarate-dependent enzymes in the maintenance of pluripotency during cellular reprogramming to induced pluripotent stem cells and in disruption of pluripotency during in vitro differentiation. The enzymatic function of many of these α-ketoglutarate-dependent proteins is required for pluripotency acquisition and maintenance. A better understanding of their specific function will be essential in furthering our knowledge of pluripotency.

Chromatin module inference on cellular trajectories identifies key transition points and poised epigenetic states in diverse developmental processes.

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.

Mediator coordinates PIC assembly with recruitment of CHD1.

Murine Chd1 (chromodomain helicase DNA-binding protein 1), a chromodomain-containing chromatin remodeling protein, is necessary for embryonic stem (ES) cell pluripotency. Chd1 binds to nucleosomes trimethylated at histone 3 Lys 4 (H3K4me3) near the beginning of active genes but not to bivalent domains also containing H3K27me3. To address the mechanism of this specificity, we reproduced H3K4me3- and CHD1-stimulated gene activation in HeLa extracts. Multidimensional protein identification technology (MuDPIT) and immunoblot analyses of purified preinitiation complexes (PICs) revealed the recruitment of CHD1 to naive chromatin but enhancement on H3K4me3 chromatin. Studies in depleted extracts showed that the Mediator coactivator complex, which controls PIC assembly, is also necessary for CHD1 recruitment. MuDPIT analyses of CHD1-associated proteins support the recruitment data and reveal numerous components of the PIC, including Mediator. In vivo, CHD1 and Mediator are recruited to an inducible gene, and genome-wide binding of the two proteins correlates well with active gene transcription in mouse ES cells. Finally, coimmunoprecipitation of CHD1 and Mediator from cell extracts can be ablated by shRNA knockdown of a specific Mediator subunit. Our data support a model in which the Mediator coordinates PIC assembly along with the recruitment of CHD1. The combined action of the PIC and H3K4me3 provides specificity in targeting CHD1 to active genes.

Elucidating Protein-DNA Interactions in Human Alphoid Chromatin via Hybridization Capture and Mass Spectrometry.

The centromere is the chromosomal locus where the kinetochore forms and is critical for ensuring proper segregation of sister chromatids during cell division. A substantial amount of effort has been devoted to understanding the characteristic features and roles of the centromere, yet some fundamental aspects of the centromere, such as the complete list of elements that define it, remain obscure. It is well-known that human centromeres include a highly repetitive class of DNA known as alpha satellite, or alphoid, DNA. We present here the first DNA-centric examination of human protein-alpha satellite interactions, employing an approach known as HyCCAPP (hybridization capture of chromatin-associated proteins for proteomics) to identify the protein components of alphoid chromatin in a human cell line. Using HyCCAPP, cross-linked alpha satellite chromatin was isolated from cell lysate, and captured proteins were analyzed via mass spectrometry. After being compared to proteins identified in control pulldown experiments, 90 proteins were identified as enriched at alphoid DNA. This list included many known centromere-binding proteins in addition to multiple novel alpha satellite-binding proteins, such as LRIF1, a heterochromatin-associated protein. The ability of HyCCAPP to reveal both known as well as novel alphoid DNA-interacting proteins highlights the validity and utility of this approach.

Histone h3 lysine 56 acetylation is linked to the core transcriptional network in human embryonic stem cells.

Lysine 56 acetylation in the helical core of histone H3 opens yeast chromatin and enables histone gene transcription, DNA replication, and DNA repair and prevents epigenetic silencing. While K56Ac is globally abundant in yeast and flies, its presence has been uncertain in mammals. We show here using mass spectrometry and genome-wide analyses that K56Ac is present in human embryonic stem cells (hESCs), overlapping strongly at active and inactive promoters with the binding of the key regulators of pluripotency, NANOG, SOX2, and OCT4. This includes also the canonical histone gene promoters and those for the hESC-specific microRNAs. K56Ac then relocates to developmental genes upon cellular differentiation. Thus the K56Ac state more accurately reflects the epigenetic differences between hESCs and somatic cells than other active histone marks such as H3 K4 trimethylation and K9 acetylation. These results suggest that K56Ac is involved in the human core transcriptional network of pluripotency.

A predictive modeling approach for cell line-specific long-range regulatory interactions.

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.

Chd1 regulates open chromatin and pluripotency of embryonic stem cells.

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.

KDM3A and KDM3B Maintain Naïve Pluripotency Through the Regulation of Alternative Splicing.

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.

DOT1L interaction partner AF10 controls patterning of H3K79 methylation and RNA polymerase II to maintain cell identity.

Histone 3 lysine 79 methylation (H3K79me) is enriched on gene bodies proportional to gene expression levels and serves as a strong barrier for the reprogramming of somatic cells to induced pluripotent stem cells (iPSCs). DOT1L is the sole histone methyltransferase that deposits all three orders-mono (me1), di (me2), and tri (me3) methylation-at H3K79. Here, we leverage genetic and chemical approaches to parse the specific functions of orders of H3K79me in maintaining cell identity. DOT1L interacts with AF10 (Mllt10), which recognizes unmodified H3K27 and boosts H3K79me2/3 methylation. AF10 deletion evicts H3K79me2/3 and reorganizes H3K79me1 to the transcription start site to facilitate iPSC formation in the absence of steady-state transcriptional changes. Instead, AF10 loss redistributes RNA polymerase II to a uniquely pluripotent pattern at highly expressed, rapidly transcribed housekeeping genes. Taken together, we reveal a specific mechanism for H3K79me2/3 located at the gene body in reinforcing cell identity.

DOT1L is a barrier to histone acetylation during reprogramming to pluripotency.

Embryonic stem cells (ESCs) have transcriptionally permissive chromatin enriched for gene activation-associated histone modifications. A striking exception is DOT1L-mediated H3K79 dimethylation (H3K79me2) that is considered a positive regulator of transcription. We find that ESCs are depleted for H3K79me2 at shared locations of enrichment with somatic cells, which are highly and ubiquitously expressed housekeeping genes, and have lower RNA polymerase II (RNAPII) at the transcription start site (TSS) despite greater nascent transcription. Inhibiting DOT1L increases the efficiency of reprogramming of somatic to induced pluripotent stem cells, enables an ESC-like RNAPII pattern at the TSS, and functionally compensates for enforced RNAPII pausing. DOT1L inhibition increases H3K27 methylation and RNAPII elongation-enhancing histone acetylation without changing the expression of the causal histone-modifying enzymes. Only the maintenance of elevated histone acetylation is essential for enhanced reprogramming and occurs at loci that are depleted for H3K79me2. Thus, DOT1L inhibition promotes the hyperacetylation and hypertranscription pluripotent properties.

Inference of cell type-specific gene regulatory networks on cell lineages from single cell omic datasets.

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.

LRRC10 regulates mammalian cardiomyocyte cell cycle during heart regeneration.

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.

Connecting the DOTs on Cell Identity.

DOT1-Like (DOT1L) is the sole methyltransferase of histone H3K79, a modification enriched mainly on the bodies of actively transcribing genes. DOT1L has been extensively studied in leukemia were some of the most frequent onco-fusion proteins contain portions of DOT1L associated factors that mislocalize H3K79 methylation and drive oncogenesis. However, the role of DOT1L in non-transformed, developmental contexts is less clear. Here we assess the known functional roles of DOT1L both in vitro cell culture and in vivo models of mammalian development. DOT1L is evicted during the 2-cell stage when cells are totipotent and massive epigenetic and transcriptional alterations occur. Embryonic stem cell lines that are derived from the blastocyst tolerate the loss of DOT1L, while the reduction of DOT1L protein levels or its catalytic activity greatly enhances somatic cell reprogramming to induced pluripotent stem cells. DOT1L knockout mice are embryonically lethal when organogenesis commences. We catalog the rapidly increasing studies of total and lineage specific knockout model systems that show that DOT1L is broadly required for differentiation. Reduced DOT1L activity is concomitant with increased developmental potential. Contrary to what would be expected of a modification that is associated with active transcription, loss of DOT1L activity results in more upregulated than downregulated genes. DOT1L also participates in various epigenetic networks that are both cell type and developmental stage specific. Taken together, the functions of DOT1L during development are pleiotropic and involve gene regulation at the locus specific and global levels.