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.

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.

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.

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.

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.

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.

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.

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.