Search

An inducible gene expression program is a hallmark of the host inflammatory response. Recently, long intergenic non-coding RNAs (lincRNAs) have been shown to regulate the magnitude, duration, and resolution of these responses. Among these is lincRNA-Cox2, a dynamically regulated gene that broadly controls immune gene expression. To evaluate the in vivo functions of this lincRNA, we characterized multiple models of lincRNA-Cox2-deficient mice. LincRNA-Cox2-deficient macrophages and murine tissues had altered expression of inflammatory genes. Transcriptomic studies from various tissues revealed that deletion of the lincRNA-Cox2 locus also strongly impaired the basal and inducible expression of the neighboring gene prostaglandin-endoperoxide synthase (Ptgs2), encoding cyclooxygenase-2, a key enzyme in the prostaglandin biosynthesis pathway. By utilizing different genetic manipulations in vitro and in vivo, we found that lincRNA-Cox2 functions through an enhancer RNA mechanism to regulate Ptgs2. More importantly, lincRNA-Cox2 also functions in trans, independently of Ptgs2, to regulate critical innate immune genes in vivo., (Published by Elsevier Inc.)

Cellular reprogramming through manipulation of defined factors holds great promise for large-scale production of cell types needed for use in therapy and for revealing principles of gene regulation. However, most reprogramming systems are inefficient, converting only a fraction of cells to the desired state. Here, we analyze MYOD-mediated reprogramming of human fibroblasts to myotubes, a well-characterized model system for direct conversion by defined factors, at pseudotemporal resolution using single-cell RNA-seq. To expose barriers to efficient conversion, we introduce a novel analytic technique, trajectory alignment, which enables quantitative comparison of gene expression kinetics across two biological processes. Reprogrammed cells navigate a trajectory with branch points that correspond to two alternative decision points, with cells that select incorrect branches terminating at aberrant or incomplete reprogramming outcomes. Analysis of these branch points revealed insulin and BMP signaling as crucial molecular determinants of reprogramming. Single-cell trajectory alignment enables rigorous quantitative comparisons between biological trajectories found in diverse processes in development, reprogramming, and other contexts., (Copyright © 2018 The Authors. Published by Elsevier Inc. All rights reserved.)

Chromosomes occupy distinct interphase territories in the three-dimensional nucleus. However, how these chromosome territories are arranged relative to one another is poorly understood. Here, we investigated the inter- chromosomal interactions between chromosomes 2q, 12, and 17 in human mesenchymal stem cells (MSCs) and MSC-derived cell types by DNA-FISH We compared our findings in normal karyotypes with a three-generation family harboring a 2q37-deletion syndrome, featuring a heterozygous partial deletion of histone deacetylase 4 ( HDAC4 ) on chr2q37. In normal karyotypes, we detected stable, recurring arrangements and interactions between the three chromosomal territories with a tissue-specific interaction bias at certain loci. These inter- chromosomal interactions were confirmed by Hi-C. Interestingly, the disease-related HDAC4 deletion resulted in displaced inter- chromosomal arrangements and altered interactions between the deletion-affected chromosome 2 and chromosome 12 and/or 17 in 2q37-deletion syndrome patients. Our findings provide evidence for a direct link between a structural chromosomal aberration and altered interphase architecture that results in a nuclear configuration, supporting a possible molecular pathogenesis., (© 2018 The Authors.)

Resistance to chemotherapy plays a significant role in cancer mortality. To identify genetic units affecting sensitivity to cytarabine, the mainstay of treatment for acute myeloid leukemia (AML), we developed a comprehensive and integrated genome-wide platform based on a dual protein-coding and non-coding integrated CRISPRa screening (DICaS). Putative resistance genes were initially identified using pharmacogenetic data from 760 human pan-cancer cell lines. Subsequently, genome scale functional characterization of both coding and long non-coding RNA (lncRNA) genes by CRISPR activation was performed. For lncRNA functional assessment, we developed a CRISPR activation of lncRNA (CaLR) strategy, targeting 14,701 lncRNA genes. Computational and functional analysis identified novel cell-cycle, survival/apoptosis, and cancer signaling genes. Furthermore, transcriptional activation of the GAS6-AS2 lncRNA, identified in our analysis, leads to hyperactivation of the GAS6/TAM pathway, a resistance mechanism in multiple cancers including AML. Thus, DICaS represents a novel and powerful approach to identify integrated coding and non-coding pathways of therapeutic relevance., (Copyright © 2018 Elsevier Inc. All rights reserved.)

The binding of the transcriptional regulator CTCF to the genome has been implicated in the formation of topologically associated domains (TADs). However, the general mechanisms of folding the genome into TADs are not fully understood. Here we test the effects of deleting a CTCF-rich locus on TAD boundary formation. Using genome-wide chromosome conformation capture (Hi-C), we focus on one TAD boundary on chromosome X harboring ~ 15 CTCF binding sites and located at the long non-coding RNA (lncRNA) locus Firre. Specifically, this TAD boundary is invariant across evolution, tissues, and temporal dynamics of X-chromosome inactivation. We demonstrate that neither the deletion of this locus nor the ectopic insertion of Firre cDNA or its ectopic expression are sufficient to alter TADs in a sex-specific or allele-specific manner. In contrast, Firre's deletion disrupts the chromatin super-loop formation of the inactive X-chromosome. Collectively, our findings suggest that apart from CTCF binding, additional mechanisms may play roles in establishing TAD boundary formation.

While gene expression dynamics have been extensively cataloged during hematopoietic differentiation in the adult, less is known about transcriptome diversity of human hematopoietic stem cells (HSCs) during development. To characterize transcriptional and post-transcriptional changes in HSCs during development, we leveraged high-throughput genomic approaches to profile miRNAs, lincRNAs, and mRNAs. Our findings indicate that HSCs manifest distinct alternative splicing patterns in key hematopoietic regulators. Detailed analysis of the splicing dynamics and function of one such regulator, HMGA2, identified an alternative isoform that escapes miRNA-mediated targeting. We further identified the splicing kinase CLK3 that, by regulating HMGA2 splicing, preserves HMGA2 function in the setting of an increase in let-7 miRNA levels, delineating how CLK3 and HMGA2 form a functional axis that influences HSC properties during development. Collectively, our study highlights molecular mechanisms by which alternative splicing and miRNA-mediated post-transcriptional regulation impact the molecular identity and stage-specific developmental features of human HSCs., (Copyright © 2018 Elsevier Inc. All rights reserved.)

In the post-genomic era, thousands of putative noncoding regulatory regions have been identified, such as enhancers, promoters, long noncoding RNAs (lncRNAs), and a cadre of small peptides. These ever-growing catalogs require high-throughput assays to test their functionality at scale. Massively parallel reporter assays have greatly enhanced the understanding of noncoding DNA elements en masse Here, we present a massively parallel RNA assay (MPRNA) that can assay 10,000 or more RNA segments for RNA-based functionality. We applied MPRNA to identify RNA-based nuclear localization domains harbored in lncRNAs. We examined a pool of 11,969 oligos densely tiling 38 human lncRNAs that were fused to a cytosolic transcript. After cell fractionation and barcode sequencing, we identified 109 unique RNA regions that significantly enriched this cytosolic transcript in the nucleus including a cytosine-rich motif. These nuclear enrichment sequences are highly conserved and over-represented in global nuclear fractionation sequencing. Importantly, many of these regions were independently validated by single-molecule RNA fluorescence in situ hybridization. Overall, we demonstrate the utility of MPRNA for future investigation of RNA-based functionalities., (© 2018 The Authors. Published under the terms of the CC BY 4.0 license.)

Imaging and chromatin capture techniques have provided important insights into our understanding of nuclear organization. A limitation of these techniques is the inability to resolve allele-specific spatiotemporal properties of genomic loci in living cells. Here, we describe an allele-specific CRISPR live-cell DNA imaging technique (SNP-CLING) to provide the first comprehensive insights into allelic positioning across space and time in mouse embryonic stem cells and fibroblasts. With 3D imaging, we studied alleles on different chromosomes in relation to one another and relative to nuclear substructures such as the nucleolus. We find that alleles maintain similar positions relative to each other and the nucleolus; however, loci occupy unique positions. To monitor spatiotemporal dynamics by SNP-CLING, we performed 4D imaging and determined that alleles are either stably positioned or fluctuating during cell state transitions, such as apoptosis. SNP-CLING is a universally applicable technique that enables the dissection of allele-specific spatiotemporal genome organization in live cells.

Transcription factors (TFs) direct developmental transitions by binding to target DNA sequences, influencing gene expression and establishing complex gene-regultory networks. To systematically determine the molecular components that enable or constrain TF activity, we investigated the genomic occupancy of FOXA2, GATA4 and OCT4 in several cell types. Despite their classification as pioneer factors, all three TFs exhibit cell-type-specific binding, even when supraphysiologically and ectopically expressed. However, FOXA2 and GATA4 can be distinguished by low enrichment at loci that are highly occupied by these factors in alternative cell types. We find that expression of additional cofactors increases enrichment at a subset of these sites. Finally, FOXA2 occupancy and changes to DNA accessibility can occur in G 1 -arrested cells, but subsequent loss of DNA methylation requires DNA replication.

The spatial organization of RNA within cells is a crucial factor influencing a wide range of biological functions throughout all kingdoms of life. However, a general understanding of RNA localization has been hindered by a lack of simple, high-throughput methods for mapping the transcriptomes of subcellular compartments. Here, we develop such a method, termed APEX-RIP, which combines peroxidase-catalyzed, spatially restricted in situ protein biotinylation with RNA-protein chemical crosslinking. We demonstrate that, using a single protocol, APEX-RIP can isolate RNAs from a variety of subcellular compartments, including the mitochondrial matrix, nucleus, cytosol, and endoplasmic reticulum (ER), with specificity and sensitivity that rival or exceed those of conventional approaches. We further identify candidate RNAs localized to mitochondria-ER junctions and nuclear lamina, two compartments that are recalcitrant to classical biochemical purification. Since APEX-RIP is simple, versatile, and does not require special instrumentation, we envision its broad application in a variety of biological contexts.

The tumor suppressor p53 is a well-characterized transcription factor that can bind gene promoters and regulate target gene transcription in response to DNA damage. Recent studies, however, have revealed that p53 binding events occur predominantly within regulatory enhancer elements. The effect of p53 binding on enhancer function has not been systematically evaluated. Here, we perform a genome-scale analysis of enhancer activity from p53-bound sequences using a series of massively parallel reporter assays (MPRAs) coupled with the assay for transposase-accessible chromatin (ATAC-Seq). We find that the majority of sequences examined display p53-dependent enhancer activity during the DNA damage response. Furthermore, we observe that p53 is bound to enhancer elements in healthy fibroblasts and poised for rapid activation in response to DNA damage. Surprisingly, our analyses revealed that most p53-bound enhancers are located within regions of inaccessible chromatin. A large subset of these enhancers become accessible following DNA damage indicating that p53 regulates their activity, in part, by modulating chromatin accessibility. The recognition and activation of enhancer elements located within inaccessible chromatin may contribute to the ability of the p53 network to function across the diverse chromatin landscapes of different tissues and cell types., (© The Author(s) 2017. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Commitment to the innate lymphoid cell (ILC) lineage is determined by Id2, a transcriptional regulator that antagonizes T and B cell-specific gene expression programs. Yet how Id2 expression is regulated in each ILC subset remains poorly understood. We identified a cis-regulatory element demarcated by a long non-coding RNA (lncRNA) that controls the function and lineage identity of group 1 ILCs, while being dispensable for early ILC development and homeostasis of ILC2s and ILC3s. The locus encoding this lncRNA, which we termed Rroid, directly interacted with the promoter of its neighboring gene, Id2, in group 1 ILCs. Moreover, the Rroid locus, but not the lncRNA itself, controlled the identity and function of ILC1s by promoting chromatin accessibility and deposition of STAT5 at the promoter of Id2 in response to interleukin (IL)-15. Thus, non-coding elements responsive to extracellular cues unique to each ILC subset represent a key regulatory layer for controlling the identity and function of ILCs., (Copyright © 2017 Elsevier Inc. All rights reserved.)

A substantial fraction of the genome is transcribed in a cell-type-specific manner, producing long non-coding RNAs (lncRNAs), rather than protein-coding transcripts. Here, we systematically characterize transcriptional dynamics during hematopoiesis and in hematological malignancies. Our analysis of annotated and de novo assembled lncRNAs showed many are regulated during differentiation and mis-regulated in disease. We assessed lncRNA function via an in vivo RNAi screen in a model of acute myeloid leukemia. This identified several lncRNAs essential for leukemia maintenance, and found that a number act by promoting leukemia stem cell signatures. Leukemia blasts show a myeloid differentiation phenotype when these lncRNAs were depleted, and our data indicates that this effect is mediated via effects on the MYC oncogene. Bone marrow reconstitutions showed that a lncRNA expressed across all progenitors was required for the myeloid lineage, whereas the other leukemia-induced lncRNAs were dispensable in the normal setting.

While long intergenic noncoding RNAs (lincRNAs) and mRNAs share similar biogenesis pathways, these transcript classes differ in many regards. LincRNAs are less evolutionarily conserved, less abundant, and more tissue-specific, suggesting that their pre- and post-transcriptional regulation is different from that of mRNAs. Here, we perform an in-depth characterization of the features that contribute to lincRNA regulation in multiple human cell lines. We find that lincRNA promoters are depleted of transcription factor (TF) binding sites, yet enriched for some specific factors such as GATA and FOS relative to mRNA promoters. Surprisingly, we find that H3K9me3-a histone modification typically associated with transcriptional repression-is more enriched at the promoters of active lincRNA loci than at those of active mRNAs. Moreover, H3K9me3-marked lincRNA genes are more tissue-specific. The most discriminant differences between lincRNAs and mRNAs involve splicing. LincRNAs are less efficiently spliced, which cannot be explained by differences in U1 binding or the density of exonic splicing enhancers but may be partially attributed to lower U2AF65 binding and weaker splicing-related motifs. Conversely, the stability of lincRNAs and mRNAs is similar, differing only with regard to the location of stabilizing protein binding sites. Finally, we find that certain transcriptional properties are correlated with higher evolutionary conservation in both DNA and RNA motifs and are enriched in lincRNAs that have been functionally characterized., (© 2017 Melé et al. Published by Cold Spring Harbor Laboratory Press.)

Long non-coding RNAs (lncRNAs) have been implicated in normal cellular homeostasis as well as pathophysiological conditions, including cancer. Here we performed global gene expression profiling of mammary epithelial cells transformed by oncogenic v-Src, and identified a large subset of uncharacterized lncRNAs potentially involved in breast cancer development. Specifically, our analysis revealed a novel lncRNA, LINC00520 that is upregulated upon ectopic expression of oncogenic v-Src, in a manner that is dependent on the transcription factor STAT3. Similarly, LINC00520 is also increased in mammary epithelial cells transformed by oncogenic PI3K and its expression is decreased upon knockdown of mutant PIK3CA. Additional expression profiling highlight that LINC00520 is elevated in a subset of human breast carcinomas, with preferential enrichment in the basal-like molecular subtype. ShRNA-mediated depletion of LINC00520 results in decreased cell migration and loss of invasive structures in 3D. RNA sequencing analysis uncovers several genes that are differentially expressed upon ectopic expression of LINC00520, a significant subset of which are also induced in v-Src-transformed MCF10A cells. Together, these findings characterize LINC00520 as a lncRNA that is regulated by oncogenic Src, PIK3CA and STAT3, and which may contribute to the molecular etiology of breast cancer.

Neutrophils, eosinophils and 'classical' monocytes collectively account for about 70% of human blood leukocytes and are among the shortest-lived cells in the body. Precise regulation of the lifespan of these myeloid cells is critical to maintain protective immune responses and minimize the deleterious consequences of prolonged inflammation. However, how the lifespan of these cells is strictly controlled remains largely unknown. Here we identify a long non-coding RNA that we termed Morrbid, which tightly controls the survival of neutrophils, eosinophils and classical monocytes in response to pro-survival cytokines in mice. To control the lifespan of these cells, Morrbid regulates the transcription of the neighbouring pro-apoptotic gene, Bcl2l11 (also known as Bim), by promoting the enrichment of the PRC2 complex at the Bcl2l11 promoter to maintain this gene in a poised state. Notably, Morrbid regulates this process in cis, enabling allele-specific control of Bcl2l11 transcription. Thus, in these highly inflammatory cells, changes in Morrbid levels provide a locus-specific regulatory mechanism that allows rapid control of apoptosis in response to extracellular pro-survival signals. As MORRBID is present in humans and dysregulated in individuals with hypereosinophilic syndrome, this long non-coding RNA may represent a potential therapeutic target for inflammatory disorders characterized by aberrant short-lived myeloid cell lifespan.

Background: Ebola virus is the causative agent of a severe syndrome in humans with a fatality rate that can approach 90 %. During infection, the host immune response is thought to become dysregulated, but the mechanisms through which this happens are not entirely understood. In this study, we analyze RNA sequencing data to determine the host response to Ebola virus infection in circulating immune cells., Results: Approximately half of the 100 genes with the strongest early increases in expression were interferon-stimulated genes, such as ISG15, OAS1, IFIT2, HERC5, MX1 and DHX58. Other highly upregulated genes included cytokines CXCL11, CCL7, IL2RA, IL2R1, IL15RA, and CSF2RB, which have not been previously reported to change during Ebola virus infection. Comparing this response in two different models of exposure (intramuscular and aerosol) revealed a similar signature of infection. The strong innate response in the aerosol model was seen not only in circulating cells, but also in primary and secondary target tissues. Conversely, the innate immune response of vaccinated macaques was almost non-existent. This suggests that the innate response is a major aspect of the cellular response to Ebola virus infection in multiple tissues., Conclusions: Ebola virus causes a severe infection in humans that is associated with high mortality. The host immune response to virus infection is thought to be an important aspect leading to severe pathology, but the components of this overactive response are not well characterized. Here, we analyzed how circulating immune cells respond to the virus and found that there is a strong innate response dependent on active virus replication. This finding is in stark contrast to in vitro evidence showing a suppression of innate immune signaling, and it suggests that the strong innate response we observe in infected animals may be an important contributor to pathogenesis.

The Linc-p21 locus, encoding a long non-coding RNA, plays an important role in p53 signaling, cell-cycle regulation, and tumor suppression. However, despite extensive study, confusion exists regarding its mechanism of action: is activity driven by the transcript acting in trans, in cis, or by an underlying functional enhancer? Here, using a knockout mouse model and a massively parallel enhancer assay, we delineate the functional elements at this locus. We observe that, even in tissues with no detectable Linc-p21 transcript, deletion of the locus significantly affects local gene expression, including of the cell-cycle regulator Cdkn1a. To characterize this RNA-independent regulatory effect, we systematically interrogated the underlying DNA sequence for enhancer activity at nucleotide resolution and confirmed the existence of multiple enhancer elements. Together, these data suggest that, in vivo, the cis-regulatory effects mediated by Linc-p21, in the presence or absence of transcription, are due to DNA enhancer elements., (Copyright © 2016 The Authors. Published by Elsevier Inc. All rights reserved.)

The complex language of eukaryotic gene expression remains incompletely understood. Despite the importance suggested by many noncoding variants statistically associated with human disease, nearly all such variants have unknown mechanisms. Here, we address this challenge using an approach based on a recent machine learning advance-deep convolutional neural networks (CNNs). We introduce the open source package Basset to apply CNNs to learn the functional activity of DNA sequences from genomics data. We trained Basset on a compendium of accessible genomic sites mapped in 164 cell types by DNase-seq, and demonstrate greater predictive accuracy than previous methods. Basset predictions for the change in accessibility between variant alleles were far greater for Genome-wide association study (GWAS) SNPs that are likely to be causal relative to nearby SNPs in linkage disequilibrium with them. With Basset, a researcher can perform a single sequencing assay in their cell type of interest and simultaneously learn that cell's chromatin accessibility code and annotate every mutation in the genome with its influence on present accessibility and latent potential for accessibility. Thus, Basset offers a powerful computational approach to annotate and interpret the noncoding genome., (© 2016 Kelley et al.; Published by Cold Spring Harbor Laboratory Press.)

There is growing evidence that transcription and nuclear organization are tightly linked. Yet, whether transcription of thousands of long noncoding RNAs (lncRNAs) could play a role in this packaging process remains elusive. Although some lncRNAs have been found to have clear roles in nuclear architecture (e.g., FIRRE, NEAT1, XIST, and others), the vast majority remain poorly understood. In this Perspective, we highlight how the act of transcription can affect nuclear architecture. We synthesize several recent findings into a proposed model where the transcription of lncRNAs can serve as guide-posts for shaping genome organization. This model is similar to the game "cat's cradle," where the shape of a string is successively changed by opening up new sites for finger placement. Analogously, transcription of lncRNAs could serve as "grip holds" for nuclear proteins to pull the genome into new positions. This model could explain general lncRNA properties such as low abundance and tissue specificity. Overall, we propose a general framework for how the act of lncRNA transcription could play a role in organizing the 3D genome., (Copyright © 2016 Elsevier Inc. All rights reserved.)

HLA-G, a nonclassical HLA molecule uniquely expressed in the placenta, is a central component of fetus-induced immune tolerance during pregnancy. The tissue-specific expression of HLA-G, however, remains poorly understood. Here, systematic interrogation of the HLA-G locus using massively parallel reporter assay (MPRA) uncovered a previously unidentified cis-regulatory element 12 kb upstream of HLA-G with enhancer activity, Enhancer L Strikingly, clustered regularly-interspaced short palindromic repeats (CRISPR)/Cas9-mediated deletion of this enhancer resulted in ablation of HLA-G expression in JEG3 cells and in primary human trophoblasts isolated from placenta. RNA-seq analysis demonstrated that Enhancer L specifically controls HLA-G expression. Moreover, DNase-seq and chromatin conformation capture (3C) defined Enhancer L as a cell type-specific enhancer that loops into the HLA-G promoter. Interestingly, MPRA-based saturation mutagenesis of Enhancer L identified motifs for transcription factors of the CEBP and GATA families essential for placentation. These factors associate with Enhancer L and regulate HLA-G expression. Our findings identify long-range chromatin looping mediated by core trophoblast transcription factors as the mechanism controlling tissue-specific HLA-G expression at the maternal-fetal interface. More broadly, these results establish the combination of MPRA and CRISPR/Cas9 deletion as a powerful strategy to investigate human immune gene regulation.

More than half the human and mouse genomes are comprised of repetitive sequences, such as transposable elements (TEs), which have been implicated in many biological processes. In contrast, much less is known about other repeats, such as local repeats that occur in multiple instances within a given locus in the genome but not elsewhere. Here, we systematically characterize local repeats in the genomic locus of the Firre long noncoding RNA (lncRNA). We find a conserved function for the RRD repeat as a ribonucleic nuclear retention signal that is sufficient to retain an otherwise cytoplasmic mRNA in the nucleus. We also identified a repeat, termed R0, that can function as a DNA enhancer element within the intronic sequences of Firre. Collectively, our data suggest that local repeats can have diverse functionalities and molecular modalities in the Firre locus and perhaps more globally in other lncRNAs.

Mounting evidence suggests that long noncoding RNAs (lncRNAs) can function as microRNA sponges and compete for microRNA binding to protein-coding transcripts. However, the prevalence, functional significance and targets of lncRNA-mediated sponge regulation of cancer are mostly unknown. Here we identify a lncRNA-mediated sponge regulatory network that affects the expression of many protein-coding prostate cancer driver genes, by integrating analysis of sequence features and gene expression profiles of both lncRNAs and protein-coding genes in tumours. We confirm the tumour-suppressive function of two lncRNAs (TUG1 and CTB-89H12.4) and their regulation of PTEN expression in prostate cancer. Surprisingly, one of the two lncRNAs, TUG1, was previously known for its function in polycomb repressive complex 2 (PRC2)-mediated transcriptional regulation, suggesting its sub-cellular localization-dependent function. Our findings not only suggest an important role of lncRNA-mediated sponge regulation in cancer, but also underscore the critical influence of cytoplasmic localization on the efficacy of a sponge lncRNA.

Background: Recent evidence suggests that RNA interaction can regulate the activity and localization of chromatin-associated proteins. However, it is unknown if these observations are specialized instances for a few key RNAs and chromatin factors in specific contexts, or a general mechanism underlying the establishment of chromatin state and regulation of gene expression., Results: Here, we perform formaldehyde RNA immunoprecipitation (fRIP-Seq) to survey the RNA associated with a panel of 24 chromatin regulators and traditional RNA binding proteins. For each protein that reproducibly bound measurable quantities of bulk RNA (90% of the panel), we detect enrichment for hundreds to thousands of both noncoding and mRNA transcripts., Conclusion: For each protein, we find that the enriched sets of RNAs share distinct biochemical, functional, and chromatin properties. Thus, these data provide evidence for widespread specific and relevant RNA association across diverse classes of chromatin-modifying complexes.

Only relatively recently has it become clear that mammalian genomes encode tens of thousands of long non-coding RNAs (lncRNAs). A striking 40% of these are expressed specifically in the brain, where they show precisely regulated temporal and spatial expression patterns. This begs the question, what is the functional role of these many lncRNA transcripts in the brain? Here we canvass a growing number of mechanistic studies that have elucidated central roles for lncRNAs in the regulation of nervous system development and function. We also survey studies indicating that neurological and psychiatric disorders may ensue when these mechanisms break down. Finally, we synthesize these insights with evidence from comparative genomics to argue that lncRNAs may have played important roles in brain evolution, by virtue of their abundant sequence innovation in mammals and plausible mechanistic connections to the adaptive processes that occurred recently in the primate and human lineages., (Copyright © 2015 Elsevier Inc. All rights reserved.)

Spinal muscular atrophy (SMA) is caused by mutations in the SMN1 gene. Because this gene is expressed ubiquitously, it remains poorly understood why motor neurons (MNs) are one of the most affected cell types. To address this question, we carried out RNA sequencing studies using fixed, antibody-labeled, and purified MNs produced from control and SMA patient-derived induced pluripotent stem cells (iPSCs). We found SMA-specific changes in MNs, including hyper-activation of the ER stress pathway. Functional studies demonstrated that inhibition of ER stress improves MN survival in vitro even in MNs expressing low SMN. In SMA mice, systemic delivery of an ER stress inhibitor that crosses the blood-brain barrier led to the preservation of spinal cord MNs. Therefore, our study implies that selective activation of ER stress underlies MN death in SMA. Moreover, the approach we have taken would be broadly applicable to the study of disease-prone human cells in heterogeneous cultures., (Copyright © 2015 Elsevier Inc. All rights reserved.)

The equivalence of human induced pluripotent stem cells (hiPSCs) and human embryonic stem cells (hESCs) remains controversial. Here we use genetically matched hESC and hiPSC lines to assess the contribution of cellular origin (hESC vs. hiPSC), the Sendai virus (SeV) reprogramming method and genetic background to transcriptional and DNA methylation patterns while controlling for cell line clonality and sex. We find that transcriptional and epigenetic variation originating from genetic background dominates over variation due to cellular origin or SeV infection. Moreover, the 49 differentially expressed genes we detect between genetically matched hESCs and hiPSCs neither predict functional outcome nor distinguish an independently derived, larger set of unmatched hESC and hiPSC lines. We conclude that hESCs and hiPSCs are molecularly and functionally equivalent and cannot be distinguished by a consistent gene expression signature. Our data further imply that genetic background variation is a major confounding factor for transcriptional and epigenetic comparisons of pluripotent cell lines, explaining some of the previously observed differences between genetically unmatched hESCs and hiPSCs.

The regulatory potential of RNA has never ceased to amaze: from RNA catalysis, to RNA-mediated splicing, to RNA-based silencing of an entire chromosome during dosage compensation. More recently, thousands of long noncoding RNA (lncRNA) transcripts have been identified, the majority with unknown function. Thus, it is tempting to think that these lncRNAs represent a cadre of new factors that function through ribonucleic mechanisms. Some evidence points to several lncRNAs with tantalizing physiological contributions and thought-provoking molecular modalities. However, dissecting the RNA biology of lncRNAs has been difficult, and distinguishing the independent contributions of functional RNAs from underlying DNA elements, or the local act of transcription, is challenging. Here, we aim to survey the existing literature and highlight future approaches that will be needed to link the RNA-based biology and mechanisms of lncRNAs in vitro and in vivo., (© 2015 Goff and Rinn; Published by Cold Spring Harbor Laboratory Press.)

Induced pluripotency is a promising avenue for disease modeling and therapy, but the molecular principles underlying this process, particularly in human cells, remain poorly understood due to donor-to-donor variability and intercellular heterogeneity. Here, we constructed and characterized a clonal, inducible human reprogramming system that provides a reliable source of cells at any stage of the process. This system enabled integrative transcriptional and epigenomic analysis across the human reprogramming timeline at high resolution. We observed distinct waves of gene network activation, including the ordered re-activation of broad developmental regulators followed by early embryonic patterning genes and culminating in the emergence of a signature reminiscent of pre-implantation stages. Moreover, complementary functional analyses allowed us to identify and validate novel regulators of the reprogramming process. Altogether, this study sheds light on the molecular underpinnings of induced pluripotency in human cells and provides a robust cell platform for further studies. PAPERCLIP., (Copyright © 2015 Elsevier Inc. All rights reserved.)

Noncoding RNAs play diverse roles throughout biology and exhibit broad functional capacity. To investigate and harness these capabilities, we developed clustered regularly interspaced short palindromic repeats (CRISPR)-Display (CRISP-Disp), a targeted localization method that uses Cas9 to deploy large RNA cargos to DNA loci. We demonstrate that functional RNA domains up to at least 4.8 kb long can be inserted in CRISPR guide RNA at multiple points, allowing the construction of Cas9 complexes with protein-binding cassettes, artificial aptamers, pools of random sequences and natural long noncoding RNAs. A unique feature of CRISP-Disp is the multiplexing of distinct functions at multiple targets, limited only by the availability of functional RNA motifs. We anticipate the use of CRISP-Disp for ectopically targeting functional RNAs and ribonucleoprotein (RNP) complexes to genomic loci.

Background: The human malaria parasite Plasmodium falciparum has a complex and multi-stage life cycle that requires extensive and precise gene regulation to allow invasion and hijacking of host cells, transmission, and immune escape. To date, the regulatory elements orchestrating these critical parasite processes remain largely unknown. Yet it is becoming increasingly clear that long non-coding RNAs (lncRNAs) could represent a missing regulatory layer across a broad range of organisms., Results: To investigate the regulatory capacity of lncRNA in P. falciparum, we harvested fifteen samples from two time-courses. Our sample set profiled 56 h of P. falciparum blood stage development. We then developed and validated strand-specific, non-polyA-selected RNA sequencing methods, and pursued the first assembly of P. falciparum strand-specific transcript structures from RNA sequencing data. This approach enabled the annotation of over one thousand lncRNA transcript models and their comprehensive global analysis: coding prediction, periodicity, stage-specificity, correlation, GC content, length, location relative to annotated transcripts, and splicing. We validated the complete splicing structure of three lncRNAs with compelling properties. Non-polyA-selected deep sequencing also enabled the prediction of hundreds of intriguing P. falciparum circular RNAs, six of which we validated experimentally., Conclusions: We found that a subset of lncRNAs, including all subtelomeric lncRNAs, strongly peaked in expression during invasion. By contrast, antisense transcript levels significantly dropped during invasion. As compared to neighboring mRNAs, the expression of antisense-sense pairs was significantly anti-correlated during blood stage development, indicating transcriptional interference. We also validated that P. falciparum produces circRNAs, which is notable given the lack of RNA interference in the organism, and discovered that a highly expressed, five-exon antisense RNA is poised to regulate P. falciparum gametocyte development 1 (PfGDV1), a gene required for early sexual commitment events.

Long noncoding RNAs (lncRNAs) have been implicated in numerous cellular processes including brain development. However, the in vivo expression dynamics and molecular pathways regulated by these loci are not well understood. Here, we leveraged a cohort of 13 lncRNAnull mutant mouse models to investigate the spatiotemporal expression of lncRNAs in the developing and adult brain and the transcriptome alterations resulting from the loss of these lncRNA loci. We show that several lncRNAs are differentially expressed both in time and space, with some presenting highly restricted expression in only selected brain regions. We further demonstrate altered regulation of genes for a large variety of cellular pathways and processes upon deletion of the lncRNA loci. Finally, we found that 4 of the 13 lncRNAs significantly affect the expression of several neighboring proteincoding genes in a cis-like manner. By providing insight into the endogenous expression patterns and the transcriptional perturbations caused by deletion of the lncRNA locus in the developing and postnatal mammalian brain, these data provide a resource to facilitate future examination of the specific functional relevance of these genes in neural development, brain function, and disease.

The tumor suppressor p53 has been studied extensively as a direct transcriptional activator of protein-coding genes. Recent studies, however, have shed light on novel regulatory functions of p53 within noncoding regions of the genome. Here, we use a systematic approach that integrates transcriptome-wide expression analysis, genome-wide p53 binding profiles and chromatin state maps to characterize the global regulatory roles of p53 in response to DNA damage. Notably, our approach identified conserved features of the p53 network in both human and mouse primary fibroblast models. In addition to known p53 targets, we identify many previously unappreciated mRNAs and long noncoding RNAs that are regulated by p53. Moreover, we find that p53 binding occurs predominantly within enhancers in both human and mouse model systems. The ability to modulate enhancer activity offers an additional layer of complexity to the p53 network and greatly expands the diversity of genomic elements directly regulated by p53., (© The Author(s) 2015. Published by Oxford University Press on behalf of Nucleic Acids Research.)

DNA methylation is a key epigenetic modification involved in regulating gene expression and maintaining genomic integrity. Here we inactivated all three catalytically active DNA methyltransferases (DNMTs) in human embryonic stem cells (ESCs) using CRISPR/Cas9 genome editing to further investigate the roles and genomic targets of these enzymes. Disruption of DNMT3A or DNMT3B individually as well as of both enzymes in tandem results in viable, pluripotent cell lines with distinct effects on the DNA methylation landscape, as assessed by whole-genome bisulfite sequencing. Surprisingly, in contrast to findings in mouse, deletion of DNMT1 resulted in rapid cell death in human ESCs. To overcome this immediate lethality, we generated a doxycycline-responsive tTA-DNMT1* rescue line and readily obtained homozygous DNMT1-mutant lines. However, doxycycline-mediated repression of exogenous DNMT1* initiates rapid, global loss of DNA methylation, followed by extensive cell death. Our data provide a comprehensive characterization of DNMT-mutant ESCs, including single-base genome-wide maps of the targets of these enzymes.

Background: Long non-coding RNAs (lncRNAs) have been implicated in diverse biological processes. In contrast to extensive genomic annotation of lncRNA transcripts, far fewer have been characterized for subcellular localization and cell-to-cell variability. Addressing this requires systematic, direct visualization of lncRNAs in single cells at single-molecule resolution., Results: We use single-molecule RNA-FISH to systematically quantify and categorize the subcellular localization patterns of a representative set of 61 lncRNAs in three different cell types. Our survey yields high-resolution quantification and stringent validation of the number and spatial positions of these lncRNA, with an mRNA set for comparison. Using this highly quantitative image-based dataset, we observe a variety of subcellular localization patterns, ranging from bright sub-nuclear foci to almost exclusively cytoplasmic localization. We also find that the low abundance of lncRNAs observed from cell population measurements cannot be explained by high expression in a small subset of 'jackpot' cells. Additionally, nuclear lncRNA foci dissolve during mitosis and become widely dispersed, suggesting these lncRNAs are not mitotic bookmarking factors. Moreover, we see that divergently transcribed lncRNAs do not always correlate with their cognate mRNA, nor do they have a characteristic localization pattern., Conclusions: Our systematic, high-resolution survey of lncRNA localization reveals aspects of lncRNAs that are similar to mRNAs, such as cell-to-cell variability, but also several distinct properties. These characteristics may correspond to particular functional roles. Our study also provides a quantitative description of lncRNAs at the single-cell level and a universally applicable framework for future study and validation of lncRNAs.

Neuronal development requires a complex choreography of transcriptional decisions to obtain specific cellular identities. Realizing the ultimate goal of identifying genome-wide signatures that define and drive specific neuronal fates has been hampered by enormous complexity in both time and space during development. Here, we have paired high-throughput purification of pyramidal neuron subclasses with deep profiling of spatiotemporal transcriptional dynamics during corticogenesis to resolve lineage choice decisions. We identified numerous features ranging from spatial and temporal usage of alternative mRNA isoforms and promoters to a host of mRNA genes modulated during fate specification. Notably, we uncovered numerous long noncoding RNAs with restricted temporal and cell-type-specific expression. To facilitate future exploration, we provide an interactive online database to enable multidimensional data mining and dissemination. This multifaceted study generates a powerful resource and informs understanding of the transcriptional regulation underlying pyramidal neuron diversity in the neocortex., (Copyright © 2015 Elsevier Inc. All rights reserved.)

Recent advances in RNA-sequencing technologies have led to the discovery of thousands of previously unannotated noncoding transcripts, including many long noncoding RNAs (lncRNAs) whose functions remain largely unknown. Here we discuss considerations and best practices in lncRNA identification and annotation, which we hope will foster functional and mechanistic exploration.

Background: Transposable elements (TEs) have significantly influenced the evolution of transcriptional regulatory networks in the human genome. Post-transcriptional regulation of human genes by TE-derived sequences has been observed in specific contexts, but has yet to be systematically and comprehensively investigated. Here, we study a collection of 75 CLIP-Seq experiments mapping the RNA binding sites for a diverse set of 51 human proteins to explore the role of TEs in post-transcriptional regulation of human mRNAs and lncRNAs via RNA-protein interactions., Results: We detect widespread interactions between RNA binding proteins (RBPs) and many families of TE-derived sequence in the CLIP-Seq data. Further, alignment coverage peaks on specific positions of the TE consensus sequences, illuminating a diversity of TE-specific RBP binding motifs. Evidence of binding and conservation of these motifs in the nonrepetitive transcriptome suggests that TEs have generally appropriated existing sequence preferences of the RBPs. Depletion assays for numerous RBPs show that TE-derived binding sites affect transcript abundance and splicing similarly to nonrepetitive sites. However, in a few cases the effect of RBP binding depends on the specific TE family bound; for example, the ubiquitously expressed RBP HuR confers transcript stability unless bound to an Alu element., Conclusions: Our meta-analysis suggests a widespread role for TEs in shaping RNA-protein regulatory networks in the human genome.

The neocortex contains an unparalleled diversity of neuronal subtypes, each defined by distinct traits that are developmentally acquired under the control of subtype-specific and pan-neuronal genes. The regulatory logic that orchestrates the expression of these unique combinations of genes is unknown for any class of cortical neuron. Here, we report that Fezf2 is a selector gene able to regulate the expression of gene sets that collectively define mouse corticospinal motor neurons (CSMN). We find that Fezf2 directly induces the glutamatergic identity of CSMN via activation of Vglut1 (Slc17a7) and inhibits a GABAergic fate by repressing transcription of Gad1. In addition, we identify the axon guidance receptor EphB1 as a target of Fezf2 necessary to execute the ipsilateral extension of the corticospinal tract. Our data indicate that co-regulated expression of neuron subtype-specific and pan-neuronal gene batteries by a single transcription factor is one component of the regulatory logic responsible for the establishment of CSMN identity.

Numerous studies over the past decade have identified increasing numbers of long noncoding RNAs (lncRNAs) across many organisms. Research since has shown that lncRNAs constitute an important layer of genome regulation in diverse biological processes and disease. Here, we discuss the common emerging theme of lncRNAs interfacing with epigenetic machinery. This, in turn, modulates the activity and localization of the epigenetic machinery during cell fate specification., (Copyright © 2014 Cold Spring Harbor Laboratory Press; all rights reserved.)

Long noncoding RNAs (lncRNAs) are emerging as important functional components in the establishment of long-range chromosomal interactions. In a recent paper published in Cell Research, Xiang et al. provide mechanistic insight into this phenomenon by characterizing the role of CCAT1-L, a colorectal cancer-specific lncRNA, in intra-chromosome looping between the MYC gene promoter and distal upstream enhancer elements that regulate MYC transcription.

Defining the transcriptional dynamics of a temporal process such as cell differentiation is challenging owing to the high variability in gene expression between individual cells. Time-series gene expression analyses of bulk cells have difficulty distinguishing early and late phases of a transcriptional cascade or identifying rare subpopulations of cells, and single-cell proteomic methods rely on a priori knowledge of key distinguishing markers. Here we describe Monocle, an unsupervised algorithm that increases the temporal resolution of transcriptome dynamics using single-cell RNA-Seq data collected at multiple time points. Applied to the differentiation of primary human myoblasts, Monocle revealed switch-like changes in expression of key regulatory factors, sequential waves of gene regulation, and expression of regulators that were not known to act in differentiation. We validated some of these predicted regulators in a loss-of function screen. Monocle can in principle be used to recover single-cell gene expression kinetics from a wide array of cellular processes, including differentiation, proliferation and oncogenic transformation.

RNA, including long noncoding RNA (lncRNA), is known to be an abundant and important structural component of the nuclear matrix. However, the molecular identities, functional roles and localization dynamics of lncRNAs that influence nuclear architecture remain poorly understood. Here, we describe one lncRNA, Firre, that interacts with the nuclear-matrix factor hnRNPU through a 156-bp repeating sequence and localizes across an ~5-Mb domain on the X chromosome. We further observed Firre localization across five distinct trans-chromosomal loci, which reside in spatial proximity to the Firre genomic locus on the X chromosome. Both genetic deletion of the Firre locus and knockdown of hnRNPU resulted in loss of colocalization of these trans-chromosomal interacting loci. Thus, our data suggest a model in which lncRNAs such as Firre can interface with and modulate nuclear architecture across chromosomes.

Although numerous approaches have been developed to map RNA-binding sites of individual RNA-binding proteins (RBPs), few methods exist that allow assessment of global RBP-RNA interactions. Here, we describe PIP-seq, a universal, high-throughput, ribonuclease-mediated protein footprint sequencing approach that reveals RNA-protein interaction sites throughout a transcriptome of interest. We apply PIP-seq to the HeLa transcriptome and compare binding sites found using different cross-linkers and ribonucleases. From this analysis, we identify numerous putative RBP-binding motifs, reveal novel insights into co-binding by RBPs, and uncover a significant enrichment for disease-associated polymorphisms within RBP interaction sites.

Many studies are uncovering functional roles for long noncoding RNAs (lncRNAs), yet few have been tested for in vivo relevance through genetic ablation in animal models. To investigate the functional relevance of lncRNAs in various physiological conditions, we have developed a collection of 18 lncRNA knockout strains in which the locus is maintained transcriptionally active. Initial characterization revealed peri- and postnatal lethal phenotypes in three mutant strains (Fendrr, Peril, and Mdgt), the latter two exhibiting incomplete penetrance and growth defects in survivors. We also report growth defects for two additional mutant strains (linc-Brn1b and linc-Pint). Further analysis revealed defects in lung, gastrointestinal tract, and heart in Fendrr(-/-) neonates, whereas linc-Brn1b(-/-) mutants displayed distinct abnormalities in the generation of upper layer II-IV neurons in the neocortex. This study demonstrates that lncRNAs play critical roles in vivo and provides a framework and impetus for future larger-scale functional investigation into the roles of lncRNA molecules. DOI: http://dx.doi.org/10.7554/eLife.01749.001.

DNA methylation was first described almost a century ago; however, the rules governing its establishment and maintenance remain elusive. Here we present data demonstrating that active transcription regulates levels of genomic methylation. We identify a novel RNA arising from the CEBPA gene locus that is critical in regulating the local DNA methylation profile. This RNA binds to DNMT1 and prevents CEBPA gene locus methylation. Deep sequencing of transcripts associated with DNMT1 combined with genome-scale methylation and expression profiling extend the generality of this finding to numerous gene loci. Collectively, these results delineate the nature of DNMT1-RNA interactions and suggest strategies for gene-selective demethylation of therapeutic targets in human diseases.