YY1 Is a Structural Regulator of Enhancer-Promoter Loops.

There is considerable evidence that chromosome structure plays important roles in gene control, but we have limited understanding of the proteins that contribute to structural interactions between gene promoters and their enhancer elements. Large DNA loops that encompass genes and their regulatory elements depend on CTCF-CTCF interactions, but most enhancer-promoter interactions do not employ this structural protein. Here, we show that the ubiquitously expressed transcription factor Yin Yang 1 (YY1) contributes to enhancer-promoter structural interactions in a manner analogous to DNA interactions mediated by CTCF. YY1 binds to active enhancers and promoter-proximal elements and forms dimers that facilitate the interaction of these DNA elements. Deletion of YY1 binding sites or depletion of YY1 protein disrupts enhancer-promoter looping and gene expression. We propose that YY1-mediated enhancer-promoter interactions are a general feature of mammalian gene control.

Super-enhancers in the control of cell identity and disease.

Super-enhancers are large clusters of transcriptional enhancers that drive expression of genes that define cell identity. Improved understanding of the roles that super-enhancers play in biology would be afforded by knowing the constellation of factors that constitute these domains and by identifying super-enhancers across the spectrum of human cell types. We describe here the population of transcription factors, cofactors, chromatin regulators, and transcription apparatus occupying super-enhancers in embryonic stem cells and evidence that super-enhancers are highly transcribed. We produce a catalog of super-enhancers in a broad range of human cell types and find that super-enhancers associate with genes that control and define the biology of these cells. Interestingly, disease-associated variation is especially enriched in the super-enhancers of disease-relevant cell types. Furthermore, we find that cancer cells generate super-enhancers at oncogenes and other genes important in tumor pathogenesis. Thus, super-enhancers play key roles in human cell identity in health and in disease.

Revisiting global gene expression analysis.

Gene expression analysis is a widely used and powerful method for investigating the transcriptional behavior of biological systems, for classifying cell states in disease, and for many other purposes. Recent studies indicate that common assumptions currently embedded in experimental and analytical practices can lead to misinterpretation of global gene expression data. We discuss these assumptions and describe solutions that should minimize erroneous interpretation of gene expression data from multiple analysis platforms.

JMnorm: a novel joint multi-feature normalization method for integrative and comparative epigenomics.

Combinatorial patterns of epigenetic features reflect transcriptional states and functions of genomic regions. While many epigenetic features have correlated relationships, most existing data normalization approaches analyze each feature independently. Such strategies may distort relationships between functionally correlated epigenetic features and hinder biological interpretation. We present a novel approach named JMnorm that simultaneously normalizes multiple epigenetic features across cell types, species, and experimental conditions by leveraging information from partially correlated epigenetic features. We demonstrate that JMnorm-normalized data can better preserve cross-epigenetic-feature correlations across different cell types and enhance consistency between biological replicates than data normalized by other methods. Additionally, we show that JMnorm-normalized data can consistently improve the performance of various downstream analyses, which include candidate cis-regulatory element clustering, cross-cell-type gene expression prediction, detection of transcription factor binding and changes upon perturbations. These findings suggest that JMnorm effectively minimizes technical noise while preserving true biologically significant relationships between epigenetic datasets. We anticipate that JMnorm will enhance integrative and comparative epigenomics.

The rate of protein synthesis in hematopoietic stem cells is limited partly by 4E-BPs.

Adult stem cells must limit their rate of protein synthesis, but the underlying mechanisms remain largely unexplored. Differences in protein synthesis among hematopoietic stem cells (HSCs) and progenitor cells did not correlate with differences in proteasome activity, total RNA content, mRNA content, or cell division rate. However, adult HSCs had more hypophosphorylated eukaryotic initiation factor 4E-binding protein 1 (4E-BP1) and 4E-BP2 as compared with most other hematopoietic progenitors. Deficiency for 4E-BP1 and 4E-BP2 significantly increased global protein synthesis in HSCs, but not in other hematopoietic progenitors, and impaired their reconstituting activity, identifying a mechanism that promotes HSC maintenance by attenuating protein synthesis.

Divergent transcription of long noncoding RNA/mRNA gene pairs in embryonic stem cells.

Many long noncoding RNA (lncRNA) species have been identified in mammalian cells, but the genomic origin and regulation of these molecules in individual cell types is poorly understood. We have generated catalogs of lncRNA species expressed in human and murine embryonic stem cells and mapped their genomic origin. A surprisingly large fraction of these transcripts (>60%) originate from divergent transcription at promoters of active protein-coding genes. The divergently transcribed lncRNA/mRNA gene pairs exhibit coordinated changes in transcription when embryonic stem cells are differentiated into endoderm. Our results reveal that transcription of most lncRNA genes is coordinated with transcription of protein-coding genes.

Transcriptional characterization of iPSC-derived microglia as a model for therapeutic development in neurodegeneration.

Microglia are the resident immune cells in the brain that play a key role in driving neuroinflammation, a hallmark of neurodegenerative disorders. Inducible microglia-like cells have been developed as an in vitro platform for molecular and therapeutic hypothesis generation and testing. However, there has been no systematic assessment of similarity of these cells to primary human microglia along with their responsiveness to external cues expected of primary cells in the brain. In this study, we performed transcriptional characterization of commercially available human inducible pluripotent stem cell (iPSC)-derived microglia-like (iMGL) cells by bulk and single cell RNA sequencing to assess their similarity with primary human microglia. To evaluate their stimulation responsiveness, iMGL cells were treated with Liver X Receptor (LXR) pathway agonists and their transcriptional responses characterized by bulk and single cell RNA sequencing. Bulk transcriptome analyses demonstrate that iMGL cells have a similar overall expression profile to freshly isolated human primary microglia and express many key microglial transcription factors and functional and disease-associated genes. Notably, at the single-cell level, iMGL cells exhibit distinct transcriptional subpopulations, representing both homeostatic and activated states present in normal and diseased primary microglia. Treatment of iMGL cells with LXR pathway agonists induces robust transcriptional changes in lipid metabolism and cell cycle at the bulk level. At the single cell level, we observe heterogeneity in responses between cell subpopulations in homeostatic and activated states and deconvolute bulk expression changes into their corresponding single cell states. In summary, our results demonstrate that iMGL cells exhibit a complex transcriptional profile and responsiveness, reminiscent of in vivo microglia, and thus represent a promising model system for therapeutic development in neurodegeneration.

Discovery and Targeting of the Signaling Controls of PNPLA3 to Effectively Reduce Transcription, Expression, and Function in Pre-Clinical NAFLD/NASH Settings.

Non-alcoholic fatty liver disease (NAFLD) and non-alcoholic steatohepatitis (NASH) are emerging worldwide epidemics, projected to become the leading cause of liver transplants. The strongest genetic risk factor for NAFLD/NASH susceptibility and progression is a single-nucleotide polymorphism (SNP) in the patatin-like phospholipase domain-containing 3 gene (PNPLA3), rs738409, encoding the missense mutation I148M. This aminoacidic substitution interferes with the normal remodeling of lipid droplets in hepatocytes. It is also thought to play a key role in promoting liver fibrosis by inhibiting the release of retinol from hepatic stellate cells. Reducing PNPLA3 levels in individuals homozygous for 148M may be an effective treatment for the entire spectrum of NAFLD, based on gene dosage analysis in the human population, as well as the protective effect of another naturally occurring SNP (rs2294918) in PNPLA3 which, when co-inherited, reduces PNPLA3 mRNA levels to 50% and counteracts disease risk. By screening a clinical compound library targeting specific signaling pathways active in primary human hepatocytes, we identified momelotinib, a drug evaluated in clinical trials to treat myelofibrosis, as a potent down-regulator of PNPLA3 expression, across all genotypes. We found that momelotinib treatment yielded >80% reduction in PNPLA3 mRNA in human primary hepatocytes and stellate cells, as well as in vivo via acute and chronic treatment of WT mice. Using a human multilineage 3D spheroid model of NASH homozygous for the PNPLA3 mutant protein, we additionally show that it decreases PNPLA3 mRNA as well as intracellular lipid content. Furthermore, we show that the effects on PNPLA3 coincide with changes in chromatin accessibility within regulatory regions of the PNPLA3 locus, consistent with inhibition occurring at the level of transcription. In addition to its primary reported targets, the JAK kinases, momelotinib inhibits several non-JAK kinases, including ACVR1. Using a combination of targeted siRNA knockdowns and signaling pathway perturbations, we show that momelotinib reduces the expression of the PNPLA3 gene largely through the inhibition of BMP signaling rather than the JAK/STAT pathway. Overall, our work identified momelotinib as a potential NASH therapeutic and uncovered previously unrecognized connections between signaling pathways and PNPLA3. These pathways may be exploited by drug modalities to "tune down" the level of gene expression, and therefore offer a potential therapeutic benefit to a high at-risk subset of NAFLD/NASH patients.

DIGIT Is a Conserved Long Noncoding RNA that Regulates GSC Expression to Control Definitive Endoderm Differentiation of Embryonic Stem Cells.

Long noncoding RNAs (lncRNAs) exhibit diverse functions, including regulation of development. Here, we combine genome-wide mapping of SMAD3 occupancy with expression analysis to identify lncRNAs induced by activin signaling during endoderm differentiation of human embryonic stem cells (hESCs). We find that DIGIT is divergent to Goosecoid (GSC) and expressed during endoderm differentiation. Deletion of the SMAD3-occupied enhancer proximal to DIGIT inhibits DIGIT and GSC expression and definitive endoderm differentiation. Disruption of the gene encoding DIGIT and depletion of the DIGIT transcript reveal that DIGIT is required for definitive endoderm differentiation. In addition, we identify the mouse ortholog of DIGIT and show that it is expressed during development and promotes definitive endoderm differentiation of mouse ESCs. DIGIT regulates GSC in trans, and activation of endogenous GSC expression is sufficient to rescue definitive endoderm differentiation in DIGIT-deficient hESCs. Our study defines DIGIT as a conserved noncoding developmental regulator of definitive endoderm.

Activation of proto-oncogenes by disruption of chromosome neighborhoods.

Oncogenes are activated through well-known chromosomal alterations such as gene fusion, translocation, and focal amplification. In light of recent evidence that the control of key genes depends on chromosome structures called insulated neighborhoods, we investigated whether proto-oncogenes occur within these structures and whether oncogene activation can occur via disruption of insulated neighborhood boundaries in cancer cells. We mapped insulated neighborhoods in T cell acute lymphoblastic leukemia (T-ALL) and found that tumor cell genomes contain recurrent microdeletions that eliminate the boundary sites of insulated neighborhoods containing prominent T-ALL proto-oncogenes. Perturbation of such boundaries in nonmalignant cells was sufficient to activate proto-oncogenes. Mutations affecting chromosome neighborhood boundaries were found in many types of cancer. Thus, oncogene activation can occur via genetic alterations that disrupt insulated neighborhoods in malignant cells.

Transcription factor trapping by RNA in gene regulatory elements.

Transcription factors (TFs) bind specific sequences in promoter-proximal and -distal DNA elements to regulate gene transcription. RNA is transcribed from both of these DNA elements, and some DNA binding TFs bind RNA. Hence, RNA transcribed from regulatory elements may contribute to stable TF occupancy at these sites. We show that the ubiquitously expressed TF Yin-Yang 1 (YY1) binds to both gene regulatory elements and their associated RNA species across the entire genome. Reduced transcription of regulatory elements diminishes YY1 occupancy, whereas artificial tethering of RNA enhances YY1 occupancy at these elements. We propose that RNA makes a modest but important contribution to the maintenance of certain TFs at gene regulatory elements and suggest that transcription of regulatory elements produces a positive-feedback loop that contributes to the stability of gene expression programs.

Genome-wide localization of small molecules.

A vast number of small-molecule ligands, including therapeutic drugs under development and in clinical use, elicit their effects by binding specific proteins associated with the genome. An ability to map the direct interactions of a chemical entity with chromatin genome-wide could provide important insights into chemical perturbation of cellular function. Here we describe a method that couples ligand-affinity capture and massively parallel DNA sequencing (Chem-seq) to identify the sites bound by small chemical molecules throughout the human genome. We show how Chem-seq can be combined with ChIP-seq to gain unique insights into the interaction of drugs with their target proteins throughout the genome of tumor cells. These methods will be broadly useful to enhance understanding of therapeutic action and to characterize the specificity of chemical entities that interact with DNA or genome-associated proteins.

A distinct small RNA pathway silences selfish genetic elements in the germline.

In the Drosophila germline, repeat-associated small interfering RNAs (rasiRNAs) ensure genomic stability by silencing endogenous selfish genetic elements such as retrotransposons and repetitive sequences. Whereas small interfering RNAs (siRNAs) derive from both the sense and antisense strands of their double-stranded RNA precursors, rasiRNAs arise mainly from the antisense strand. rasiRNA production appears not to require Dicer-1, which makes microRNAs (miRNAs), or Dicer-2, which makes siRNAs, and rasiRNAs lack the 2',3' hydroxy termini characteristic of animal siRNA and miRNA. Unlike siRNAs and miRNAs, rasiRNAs function through the Piwi, rather than the Ago, Argonaute protein subfamily. Our data suggest that rasiRNAs protect the fly germline through a silencing mechanism distinct from both the miRNA and RNA interference pathways.

A single Argonaute protein mediates both transcriptional and posttranscriptional silencing in Schizosaccharomyces pombe.

The Schizosaccharomyces pombe genome encodes only one of each of the three major classes of proteins implicated in RNA silencing: Dicer (Dcr1), RNA-dependent RNA polymerase (RdRP; Rdp1), and Argonaute (Ago1). These three proteins are required for silencing at centromeres and for the initiation of transcriptionally silent heterochromatin at the mating-type locus. Here, we show that the introduction of a double-stranded RNA (dsRNA) hairpin corresponding to a green fluorescent protein (GFP) transgene triggers classical RNA interference (RNAi) in S. pombe. That is, GFP silencing triggered by dsRNA reflects a change in the steady-state concentration of GFP mRNA, but not in the rate of GFP transcription. RNAi in S. pombe requires dcr1, rdp1, and ago1, but does not require chp1, tas3, or swi6, genes required for transcriptional silencing. Thus, the RNAi machinery in S. pombe can direct both transcriptional and posttranscriptional silencing using a single Dicer, RdRP, and Argonaute protein. Our findings suggest that these three proteins fulfill a common biochemical function in distinct siRNA-directed silencing pathways.