Diverse Cellular Contexts Research Articles

A heat shock–responsive lncRNA Heat acts as a HSF1-directed transcriptional brake via m6A modification

Long noncoding RNAs (lncRNAs) are key regulators of gene expression in diverse cellular contexts and biological processes. Given the surprising range of shapes and sizes, how distinct lncRNAs achieve functional specificity remains incompletely understood. Here, we identified a heat shock-inducible lncRNA, Heat, in mouse cells that acts as a transcriptional brake to restrain stress gene expression. Functional characterization reveals that Heat directly binds to heat shock transcription factor 1 (HSF1), thereby targeting stress genes in a trans-acting manner. Intriguingly, Heat is heavily methylated in the form of m6A. Although dispensable for HSF1 binding, Heat methylation is required for silencing stress genes to attenuate heat shock response. Consistently, m6A depletion results in prolonged activation of stress genes. Furthermore, Heat mediates these effects via the nuclear m6A reader YTHDC1, forming a transcriptional silencing complex for stress genes. Our study reveals a crucial role of nuclear epitranscriptome in the transcriptional regulation of heat shock response.

Modeling the Evolutionary Architectures of Transcribed Human Enhancer Sequences Reveals Distinct Origins, Functions, and Associations with Human Trait Variation.

Despite the importance of gene regulatory enhancers in human biology and evolution, we lack a comprehensive model of enhancer evolution and function. This substantially limits our understanding of the genetic basis of species divergence and our ability to interpret the effects of noncoding variants on human traits.To explore enhancer sequence evolution and its relationship to regulatory function, we traced the evolutionary origins of transcribed human enhancer sequences with activity across diverse tissues and cellular contexts from the FANTOM5 consortium. The transcribed enhancers are enriched for sequences of a single evolutionary age (“simple” evolutionary architectures) compared with enhancers that are composites of sequences of multiple evolutionary ages (“complex” evolutionary architectures), likely indicating constraint against genomic rearrangements. Complex enhancers are older, more pleiotropic, and more active across species than simple enhancers. Genetic variants within complex enhancers are also less likely to associate with human traits and biochemical activity. Transposable-element-derived sequences (TEDS) have made diverse contributions to enhancers of both architectures; the majority of TEDS are found in enhancers with simple architectures, while a minority have remodeled older sequences to create complex architectures. Finally, we compare the evolutionary architectures of transcribed enhancers with histone-mark-defined enhancers.Our results reveal that most human transcribed enhancers are ancient sequences of a single age, and thus the evolution of most human enhancers was not driven by increases in evolutionary complexity over time. Our analyses further suggest that considering enhancer evolutionary histories provides context that can aid interpretation of the effects of variants on enhancer function. Based on these results, we propose a framework for analyzing enhancer evolutionary architecture.

DNA nuclear targeting sequences for enhanced non-viral gene transfer: An in vitro and in vivo study

An important bottleneck for non-viral gene transfer commonly relates to translocation of nucleic acids into the nuclear compartment of target cells. So-called 3NFs are optimized short nucleotide sequences able to interact with the transcription factor nuclear factor κB (NF-κB), which can enhance the nuclear import of plasmid DNA (pDNA) carrying such motifs. In this work, we first designed a consistent set of six pDNAs featuring a common backbone and only varying in their 3NF sequences. These constructions were then transfected under various experimental settings. In vitro, cationic polymer-assisted pDNA delivery in five human-derived cell lines showed the potential advantage of 3NF carrying pDNA in diverse cellular contexts. In vivo, naked pDNAs were hydrodynamically delivered to muscle hindlimbs in healthy mice; this direct accurate comparative (in the absence of any gene carrier) revealed modest but consistent trends in favor of the pDNAs equipped with 3NF. In summary, the results reported emphasize the implications of various parameters on NF-κB-mediated pDNA nuclear import; under specific conditions, 3NF can provide modest to substantial advantages for pDNA gene transfer, in vitro as well as in vivo. This study thus further underscores the potential of optimized nuclear import for more efficient non-viral gene transfer applications.

The journey of Ca2+ through the cell - pulsing through the network of ER membrane contact sites.

Calcium is the third most abundant metal on earth, and the fundaments of its homeostasis date back to pre-eukaryotic life forms. In higher organisms, Ca2+ serves as a cofactor for a wide array of (enzymatic) interactions in diverse cellular contexts and constitutes the most important signaling entity in excitable cells. To enable responsive behavior, cytosolic Ca2+ concentrations are kept low through sequestration into organellar stores, particularly the endoplasmic reticulum (ER), but also mitochondria and lysosomes. Specific triggers are then used to instigate a local release of Ca2+ on demand. Here, communication between organelles comes into play, which is accomplished through intimate yet dynamic contacts, termed membrane contact sites (MCSs). The field of MCS biology in relation to cellular Ca2+ homeostasis has exploded in recent years. Taking advantage of this new wealth of knowledge, in this Review, we invite the reader on a journey of Ca2+ flux through the ER and its associated MCSs. New mechanistic insights and technological advances inform the narrative on Ca2+ acquisition and mobilization at these sites of communication between organelles, and guide the discussion of their consequences for cellular physiology.

A functional variant on 20q13.33 related to glioma risk alters enhancer activity and modulates expression of multiple genes.

Genome‐wide association studies (GWAS) have identified single‐nucleotide polymorphisms (SNPs) associated with glioma risk on 20q13.33, but the biological mechanisms underlying this association are unknown. We tested the hypothesis that a functional SNP on 20q13.33 impacted the activity of an enhancer, leading to an altered expression of nearby genes. To identify candidate functional SNPs, we identified all SNPs in linkage disequilibrium with the risk‐associated SNP rs2297440 that mapped to putative enhancers. Putative enhancers containing candidate functional SNPs were tested for allele‐specific effects in luciferase enhancer activity assays against glioblastoma multiforme (GBM) cell lines. An enhancer containing SNP rs3761124 exhibited allele‐specific effects on activity. Deletion of this enhancer by CRISPR‐Cas9 editing in GBM cell lines correlated with an altered expression of multiple genes, including STMN3, RTEL1, RTEL1‐TNFRSF6B, GMEB2, and SRMS. Expression quantitative trait loci (eQTL) analyses using nondiseased brain samples, isocitrate dehydrogenase 1 (IDH1) wild‐type glioma, and neurodevelopmental tissues showed STMN3 to be a consistent significant eQTL with rs3761124. RTEL1 and GMEB2 were also significant eQTLs in the context of early CNS development and/or in IDH1 wild‐type glioma. We provide evidence that rs3761124 is a functional variant on 20q13.33 related to glioma/GBM risk that modulates the expression of STMN3 and potentially other genes across diverse cellular contexts.

Mutation of NEKL-4/NEK10 and TTLL genes suppress neuronal ciliary degeneration caused by loss of CCPP-1 deglutamylase function.

Ciliary microtubules are subject to post-translational modifications that act as a “Tubulin Code” to regulate motor traffic, binding proteins and stability. In humans, loss of CCP1, a cytosolic carboxypeptidase and tubulin deglutamylating enzyme, causes infantile-onset neurodegeneration. In C. elegans, mutations in ccpp-1, the homolog of CCP1, result in progressive degeneration of neuronal cilia and loss of neuronal function. To identify genes that regulate microtubule glutamylation and ciliary integrity, we performed a forward genetic screen for suppressors of ciliary degeneration in ccpp-1 mutants. We isolated the ttll-5(my38) suppressor, a mutation in a tubulin tyrosine ligase-like glutamylase gene. We show that mutation in the ttll-4, ttll-5, or ttll-11 gene suppressed the hyperglutamylation-induced loss of ciliary dye filling and kinesin-2 mislocalization in ccpp-1 cilia. We also identified the nekl-4(my31) suppressor, an allele affecting the NIMA (Never in Mitosis A)-related kinase NEKL-4/NEK10. In humans, NEK10 mutation causes bronchiectasis, an airway and mucociliary transport disorder caused by defective motile cilia. C. elegans NEKL-4 localizes to the ciliary base but does not localize to cilia, suggesting an indirect role in ciliary processes. This work defines a pathway in which glutamylation, a component of the Tubulin Code, is written by TTLL-4, TTLL-5, and TTLL-11; is erased by CCPP-1; is read by ciliary kinesins; and its downstream effects are modulated by NEKL-4 activity. Identification of regulators of microtubule glutamylation in diverse cellular contexts is important to the development of effective therapies for disorders characterized by changes in microtubule glutamylation. By identifying C. elegans genes important for neuronal and ciliary stability, our work may inform research into the roles of the tubulin code in human ciliopathies and neurodegenerative diseases.

Sequence-Independent Self-Assembly of Germ Granule mRNAs into Homotypic Clusters

mRNAs enriched in membraneless condensates provide functional compartmentalization within cells. The mechanisms that recruit transcripts to condensates are under intense study; however, how mRNAs organize once they reach a granule remains poorly understood. Here, we report on a self-sorting mechanism by which multiple mRNAs derived from the same gene assemble into discrete homotypic clusters. We demonstrate that invivo mRNA localization to granules and self-assembly within granules are governed by different mRNA features: localization is encoded by specific RNA regions, whereas self-assembly involves the entire mRNA, does not involve sequence-specific, ordered intermolecular RNA:RNA interactions, and is thus RNA sequence independent. We propose that the ability of mRNAs to self-sort into homotypic assemblies is an inherent property of an messenger ribonucleoprotein (mRNP) that is augmented under conditions that increase RNA concentration, such as upon enrichment in RNA-protein granules, a process that appears conserved in diverse cellular contexts and organisms.

USP7-dependent control of myogenin stability is required for terminal differentiation in skeletal muscle progenitors.

Satellite cells (SCs) are myogenic progenitors responsible for skeletal muscle regeneration and maintenance. Upon activation, SCs enter a phase of robust proliferation followed by terminal differentiation. Underlying this myogenic progression, the sequential expression of muscle regulatory transcription factors (MRFs) and the downregulation of transcription factor paired box gene 7 (Pax7) are key steps regulating SC fate. In addition to transcriptional regulation, post-translational control of Pax7 and the MRFs provides another layer of spatiotemporal control to the myogenic process. In this context, previous work showed that Pax7 is ubiquitinated by the E3 ligase neural precursor cell-expressed developmentally downregulated protein 4 and interacts with several proteins related to the ubiquitin-proteasome system, including the deubiquitinase ubiquitin-specific protease 7 (USP7). Although USP7 functions in diverse cellular contexts, its role(s) during myogenesis remains poorly explored. Here, we show that USP7 is transiently expressed in adult muscle progenitors, correlating with the onset of myogenin expression, while it is downregulated in newly formed myotubes/myofibers. Acute inhibition of USP7 activity upon muscle injury results in persistent expression of early regeneration markers and a significant reduction in the diameter of regenerating myofibers. At the molecular level, USP7 downregulation or pharmacological inhibition impairs muscle differentiation by affecting myogenin stability. Co-immunoprecipitation and invitro activity assays indicate that myogenin is a novel USP7 target for deubiquitination. These results suggest that USP7 regulates SC myogenic progression by enhancing myogenin stability.

Sequence Independent Self-Assembly of Germ Granule mRNAs into Homotypic Clusters

mRNAs enriched in membraneless condensates provide functional compartmentalization within cells. How mRNAs sort and organize within these condensates is poorly understood. Here, we report on a self-sorting mechanism by which multiple mRNAs derived from the same gene assemble into homotypic mRNA clusters. In Drosophila germ granules, clusters formed by mRNAs originating from different genes are de-mixed from each other within the same granule and located at distinct positions within the homogeneously-distributed protein environment of the granule. The specificity for self-sorting is not dependent on any particular RNA sequence and the intermolecular RNA:RNA interactions among clustered mRNAs are transient. We propose that the ability of mRNAs to self-sort into homotypic assemblies is their inherent property that is augmented under conditions that increase local mRNA concentration, such as upon enrichment in RNA-protein condensates, a process that appears conserved in diverse cellular contexts and organisms.

Integrating regulatory DNA sequence and gene expression to predict genome-wide chromatin accessibility across cellular contexts.

MotivationGenome-wide profiles of chromatin accessibility and gene expression in diverse cellular contexts are critical to decipher the dynamics of transcriptional regulation. Recently, convolutional neural networks have been used to learn predictive cis-regulatory DNA sequence models of context-specific chromatin accessibility landscapes. However, these context-specific regulatory sequence models cannot generalize predictions across cell types.ResultsWe introduce multi-modal, residual neural network architectures that integrate cis-regulatory sequence and context-specific expression of trans-regulators to predict genome-wide chromatin accessibility profiles across cellular contexts. We show that the average accessibility of a genomic region across training contexts can be a surprisingly powerful predictor. We leverage this feature and employ novel strategies for training models to enhance genome-wide prediction of shared and context-specific chromatin accessible sites across cell types. We interpret the models to reveal insights into cis- and trans-regulation of chromatin dynamics across 123 diverse cellular contexts.Availability and implementationThe code is available at https://github.com/kundajelab/ChromDragoNN.Supplementary information Supplementary data are available at Bioinformatics online.

A Machine Learning Approach for the Discovery of Ligand-Specific Functional Mechanisms of GPCRs.

G protein-coupled receptors (GPCRs) play a key role in many cellular signaling mechanisms, and must select among multiple coupling possibilities in a ligand-specific manner in order to carry out a myriad of functions in diverse cellular contexts. Much has been learned about the molecular mechanisms of ligand-GPCR complexes from Molecular Dynamics (MD) simulations. However, to explore ligand-specific differences in the response of a GPCR to diverse ligands, as is required to understand ligand bias and functional selectivity, necessitates creating very large amounts of data from the needed large-scale simulations. This becomes a Big Data problem for the high dimensionality analysis of the accumulated trajectories. Here we describe a new machine learning (ML) approach to the problem that is based on transforming the analysis of GPCR function-related, ligand-specific differences encoded in the MD simulation trajectories into a representation recognizable by state-of-the-art deep learning object recognition technology. We illustrate this method by applying it to recognize the pharmacological classification of ligands bound to the 5-HT2A and D2 subtypes of class-A GPCRs from the serotonin and dopamine families. The ML-based approach is shown to perform the classification task with high accuracy, and we identify the molecular determinants of the classifications in the context of GPCR structure and function. This study builds a framework for the efficient computational analysis of MD Big Data collected for the purpose of understanding ligand-specific GPCR activity.

LncRedibly versatile: biochemical and biological functions of long noncoding RNAs.

Long noncoding RNAs (lncRNAs) are transcripts that do not code for proteins, but nevertheless exert regulatory effects on various biochemical pathways, in part via interactions with proteins, DNA, and other RNAs. LncRNAs are thought to regulate transcription and other biological processes by acting, for example, as guides that target proteins to chromatin, scaffolds that facilitate protein-protein interactions and complex formation, and orchestrators of phase-separated compartments. The study of lncRNAs has reached an exciting time, as recent advances in experimental and computational methods allow for genome-wide interrogation of biochemical and biological mechanisms of these enigmatic transcripts. A better appreciation for the biochemical versatility of lncRNAs has allowed us to begin closing gaps in our knowledge of how they act in diverse cellular and organismal contexts, including development and disease.

Synchronous Transcription of Cytoskeleton-Associated Genes is Critical to Cotton Fiber Elongation

Cytoskeleton-associated genes’ expression in the modern cotton fibers has selectively been evolved under domestication. Since the magnitude of transcriptional dynamicity of cytoskeletal genes in fibers has not yet been thoroughly explored, comprehensive expression analyses of the domestication-driven profilins (GhPRFs) were performed by assaying 12 tissues from various developmental stages of elite cotton genotypes. Results from 3240 successful profilin (GhPRF1–GhPRF6) expression assays revealed enormous spatial transcriptional biases with their predominant expression in the elongating fibers. Moreover, fiber developmental stages have more pronounced impacts on differential gene expression than the genetic backgrounds. Transcriptional comparisons of GhPRFs-associated key cytoskeletal genes GhAnnex and GhACT show perfect (r = 1.0) to strong positive (r > 0.95) correlation in the developing fibers. RNAi-mediated constitutive suppression of domestication-driven GhPRF1 gene resulted in significant expression modulation of GhAnnex and GhACT genes in floral/fiber tissues and consequently sustained the cell wall development. These results highlight that temporal expression of GhPRF–GhAnnex–GhACT triad is exceptionally coordinated which illustrates fascinating crosstalk in diverse cellular contexts and cell types. Temporal fold-expression variations of other cytoskeletal-‘structural’ genes Tub–Exp–FLA triad were transversely synchronous to the GhPRF–GhAnnex–GhACT genic triad. However, the CMW-‘modifying’ genes GalT1, ABP1, and CesA were unconditionally divergent from the synchronous temporal dynamics of genic triad investigated beforehand. To our knowledge, this study represents the first in-depth description of synchronous transcriptional dynamicity of cytoskeletal-‘structural’ genes in the developing fibers. These results present a novel viewpoint on the temporal expression evolution of cytoskeletal-‘structural’ genes, augmenting our understanding of the evolutionary genetics of cytoskeletal genes in fiber development.

Epigenetics in Insects: Genome Regulation and the Generation of Phenotypic Diversity.

Epigenetic inheritance is fundamentally important to cellular differentiation and developmental plasticity. In this review, we provide an introduction to the field of molecular epigenetics in insects. Epigenetic information is passed across cell divisions through the methylation of DNA, the modification of histone proteins, and the activity of noncoding RNAs. Much of our knowledge of insect epigenetics has been gleaned from a few model species. However, more studies of epigenetic information in traditionally nonmodel taxa will help advance our understanding of the developmental and evolutionary significance of epigenetic inheritance in insects. To this end, we also provide a brief overview of techniques for profiling and perturbing individual facets of the epigenome. Doing so in diverse cellular, developmental, and taxonomic contexts will collectively help shed new light on how genome regulation results in the generation of diversity in insect form and function.

A Hearty Dose of Noncoding RNAs: The Imprinted DLK1-DIO3 Locus in Cardiac Development and Disease.

The imprinted Dlk1-Dio3 genomic region harbors a noncoding RNA cluster encoding over fifty microRNAs (miRNAs), three long noncoding RNAs (lncRNAs), and a small nucleolar RNA (snoRNA) gene array. These distinct noncoding RNAs (ncRNAs) are thought to arise from a single polycistronic transcript that is subsequently processed into individual ncRNAs, each with important roles in diverse cellular contexts. Considering these ncRNAs are derived from a polycistron, it is possible that some coordinately regulate discrete biological processes in the heart. Here, we provide a comprehensive summary of Dlk1-Dio3 miRNAs and lncRNAs, as they are currently understood in the cellular and organ-level context of the cardiovascular system. Highlighted are expression profiles, mechanistic contributions, and functional roles of these ncRNAs in heart development and disease. Notably, a number of these ncRNAs are implicated in processes often perturbed in heart disease, including proliferation, differentiation, cell death, and fibrosis. However, most literature falls short of characterizing precise mechanisms for many of these ncRNAs, warranting further investigation. Taken together, the Dlk1-Dio3 locus represents a largely unexplored noncoding regulator of cardiac homeostasis, harboring numerous ncRNAs that may serve as therapeutic targets for cardiovascular disease.

Abstract 954: Predicting synergistic drug combinations and resistance mechanisms from genomic features and single-agent response profiles

Abstract Drug combinations promise to improve clinical responses and/or forestall drug resistance. To capitalize on this promise, we need to know which drugs to combine, and which patients to give them to based on the genetic or pathological features of their disease. However, progress towards this goal has been hindered by the infeasibility of performing comprehensive drug-combination studies across thousands of cellular contexts. We hypothesized that the basal gene-transcription state of cancer cell lines, in concert with the cell-viability profiles of single-agent small molecules, might be leveraged to nominate specific synergistic drug combinations and identify mechanisms of drug resistance, eliminating the need to test all possible drug/drug combinations across cellular models. Specifically, we predicted that inhibiting the protein product of transcripts associated with drug resistance to a given small molecule might induce drug synergy. To test this notion, we analyzed nearly 400,000 drug-sensitivity profiles in >800 cancer cell lines to identify candidate compound-gene pairs. We identified over 100 examples where outlier expression of a single transcript was correlated with resistance to a small molecule. Of these gene/drug pairs, 9 genes represented imminently druggable targets, including established clinically-relevant relationships between the alkylating agent temozolomide and MGMT expression, and between a subset of chemotherapeutics including paclitaxel and the efflux pump ABCB1. Inhibition of candidate “co-targets”, which included 3 previously characterized relationships and 6 novel relationships, resulted in cell-line-specific synergistic cell killing across multiple cell-line models. For validated compound-gene pairs, exogenous expression of the “co-target” was sufficient to confer resistance. For example, we found that high expression of MGLL, encoding monoglyceride lipase, was uniquely associated with lack of response to the histone lysine demethylase inhibitor GSK-J4. Endogenous or exogenous MGLL expression conferred resistance to GSK-J4, while MGLL-proficient cell lines could be sensitized to GSK-J4 up to 50-fold by co-treatment with an irreversible MGLL inhibitor. These initial studies highlight the potential of integrating basal gene expression features with small-molecule response to nominate rational candidates for drug combinations. As public repositories of single agent response data from diverse cellular contexts continue to expand, so too will our repertoire of therapeutic combinations. Moreover, this approach permits the parallel identification of genomic features that indicate which patient populations are most likely to benefit from such combinations. Citation Format: Matthew G. Rees, Lisa Brenan, Patrick Duggan, Cory M. Johannessen. Predicting synergistic drug combinations and resistance mechanisms from genomic features and single-agent response profiles [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 954.

Why does biopolymer condensation matter?

Biomolecular condensation that occurs through phase separation is an emerging topic in cell biology. Ongoing research will continue to provide novel mechanistic insights into the involvement of biological condensates in a variety of normal and pathological cellular processes. In particular, a deeper understanding of the physical basis for biomolecular phase separation and of the properties of biopolymer condensates should emerge in the near future. Accordingly, connections between these processes and the organization of matter in living organisms will undergo fresh scrutiny. Biomolecules can phase separate and form condensates that have roles in diverse cellular processes and contexts. Michnick and Bergeron-Sandoval comment on this rapidly progressing field and envisage that the study of biological phase separation will bring new understanding of cell and developmental biology.

RNA localization regulates diverse and dynamic cellular processes.

At the nexus of specialized cellular responses are localized enrichments of protein activity. The localization of messenger RNA (mRNA) coupled with translational control often plays a crucial role in the generation of protein concentrations at defined subcellular domains. Although mRNA localization is classically associated with large specialized cells, such as neurons and embryos, RNA localization is a highly conserved paradigm of post-transcriptional regulation observed in diverse cellular contexts. Functions of localized mRNAs extend far beyond the well-studied examples of neuronal polarization and developmental patterning. Since the initial discovery of the intracellular localization of cytoskeletal mRNAs within migrating cells, hundreds of mRNAs are now known to be enriched at specific organelles where they contribute to cell function. In this short review, we discuss basic principles regulating RNA localization and consider the contribution of localized mRNA to several essential cellular behaviors. We consider RNA localization as a mechanism with widespread implications for cellular function.

Nuclear Import Receptor Inhibits Phase Separation of FUS through Binding to Multiple Sites

Liquid-liquid phase separation (LLPS) is believed to underlie formation of biomolecular condensates, cellular compartments that concentrate macromolecules without surrounding membranes. Physical mechanisms that control condensate formation/dissolution are poorly understood. The RNA-binding protein fused in sarcoma (FUS) undergoes LLPS invitro and associates with condensates in cells. We show that the importin karyopherin-β2/transportin-1 inhibits LLPS of FUS. This activity depends on tight binding of karyopherin-β2 to the C-terminal proline-tyrosine nuclear localization signal (PY-NLS) of FUS. Nuclear magnetic resonance (NMR) analyses reveal weak interactions of karyopherin-β2 with sequence elements and structural domains distributed throughout the entirety of FUS. Biochemical analyses demonstrate that most of these same regions also contribute to LLPS of FUS. The data lead to a model where high-affinity binding of karyopherin-β2 to the FUS PY-NLS tethers the proteins together, allowing multiple, distributed weak intermolecular contacts to disrupt FUS self-association, blocking LLPS. Karyopherin-β2 may act analogously to control condensates in diverse cellular contexts.

Dynamics of promoter bivalency and RNAP II pausing in mouse stem and differentiated cells

BackgroundMammalian embryonic stem cells display a unique epigenetic and transcriptional state to facilitate pluripotency by maintaining lineage-specification genes in a poised state. Two epigenetic and transcription processes involved in maintaining poised state are bivalent chromatin, characterized by the simultaneous presence of activating and repressive histone methylation marks, and RNA polymerase II (RNAPII) promoter proximal pausing. However, the dynamics of histone modifications and RNAPII at promoters in diverse cellular contexts remains underexplored.ResultsWe collected genome wide data for bivalent chromatin marks H3K4me3 and H3K27me3, and RNAPII (8WG16) occupancy together with expression profiling in eight different cell types, including ESCs, in mouse. The epigenetic and transcription profiles at promoters grouped in over thirty clusters with distinct functional identities and transcription control.ConclusionThe clustering analysis identified distinct bivalent clusters where genes in one cluster retained bivalency across cell types while in the other were mostly cell type specific, but neither showed a high RNAPII pausing. We noted that RNAPII pausing is more associated with active genes than bivalent genes in a cell type, and was globally reduced in differentiated cell types compared to multipotent.