Long Noncoding RNAs: Molecular Modalities to Organismal Functions.

We have known for decades that long noncoding RNAs (lncRNAs) can play essential functions across most forms of life. The maintenance of chromosome length requires an lncRNA (e.g., hTERC) and two lncRNAs in the ribosome that are required for protein synthesis. Thus, lncRNAs can represent powerful RNA machines. More recently, it has become clear that mammalian genomes encode thousands more lncRNAs. Thus, we raise the question: Which, if any, of these lncRNAs could also represent RNA-based machines? Here we synthesize studies that are beginning to address this question by investigating fundamental properties of lncRNA genes, revealing new insights into the RNA structure-function relationship, determining cis- and trans-acting lncRNAs in vivo, and generating new developments in high-throughput screening used to identify functional lncRNAs. Overall, these findings provide a context toward understanding the molecular grammar underlying lncRNA biology.

Long non-coding RNAs: definitions, functions, challenges and recommendations.

Genes specifying long non-coding RNAs (lncRNAs) occupy a large fraction of the genomes of complex organisms. The term 'lncRNAs' encompasses RNA polymerase I (Pol I), Pol II and Pol III transcribed RNAs, and RNAs from processed introns. The various functions of lncRNAs and their many isoforms and interleaved relationships with other genes make lncRNA classification and annotation difficult. Most lncRNAs evolve more rapidly than protein-coding sequences, are cell type specific and regulate many aspects of cell differentiation and development and other physiological processes. Many lncRNAs associate with chromatin-modifying complexes, are transcribed from enhancers and nucleate phase separation of nuclear condensates and domains, indicating an intimate link between lncRNA expression and the spatial control of gene expression during development. lncRNAs also have important roles in the cytoplasm and beyond, including in the regulation of translation, metabolism and signalling. lncRNAs often have a modular structure and are rich in repeats, which are increasingly being shown to be relevant to their function. In this Consensus Statement, we address the definition and nomenclature of lncRNAs and their conservation, expression, phenotypic visibility, structure and functions. We also discuss research challenges and provide recommendations to advance the understanding of the roles of lncRNAs in development, cell biology and disease.

Single-Cell Profiling of Ebola Virus Disease In Vivo Reveals Viral and Host Dynamics.

Ebola virus (EBOV) causes epidemics with high mortality yet remains understudied due to the challenge of experimentation in high-containment and outbreak settings. Here, we used single-cell transcriptomics and CyTOF-based single-cell protein quantification to characterize peripheral immune cells during EBOV infection in rhesus monkeys. We obtained 100,000 transcriptomes and 15,000,000 protein profiles, finding that immature, proliferative monocyte-lineage cells with reduced antigen-presentation capacity replace conventional monocyte subsets, while lymphocytes upregulate apoptosis genes and decline in abundance. By quantifying intracellular viral RNA, we identify molecular determinants of tropism among circulating immune cells and examine temporal dynamics in viral and host gene expression. Within infected cells, EBOV downregulates STAT1 mRNA and interferon signaling, and it upregulates putative pro-viral genes (e.g., DYNLL1 and HSPA5), nominating pathways the virus manipulates for its replication. This study sheds light on EBOV tropism, replication dynamics, and elicited immune response and provides a framework for characterizing host-virus interactions under maximum containment.

An Integrated Genome-wide CRISPRa Approach to Functionalize lncRNAs in Drug Resistance.

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.

Integrative Analyses of Human Reprogramming Reveal Dynamic Nature of Induced Pluripotency.

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.

From genotype to phenotype: genetics of mammalian long non-coding RNAs in vivo.

Genome-wide sequencing has led to the discovery of thousands of long non-coding RNA (lncRNA) loci in the human genome, but evidence of functional significance has remained controversial for many lncRNAs. Genetically engineered model organisms are considered the gold standard for linking genotype to phenotype. Recent advances in CRISPR-Cas genome editing have led to a rapid increase in the use of mouse models to more readily survey lncRNAs for functional significance. Here, we review strategies to investigate the physiological relevance of lncRNA loci by highlighting studies that have used genetic mouse models to reveal key in vivo roles for lncRNAs, from fertility to brain development. We illustrate how an investigative approach, starting with whole-gene deletion followed by transcription termination and/or transgene rescue strategies, can provide definitive evidence for the in vivo function of mammalian lncRNAs.

Identifying recent adaptations in large-scale genomic data.

Although several hundred regions of the human genome harbor signals of positive natural selection, few of the relevant adaptive traits and variants have been elucidated. Using full-genome sequence variation from the 1000 Genomes (1000G) Project and the composite of multiple signals (CMS) test, we investigated 412 candidate signals and leveraged functional annotation, protein structure modeling, epigenetics, and association studies to identify and extensively annotate candidate causal variants. The resulting catalog provides a tractable list for experimental follow-up; it includes 35 high-scoring nonsynonymous variants, 59 variants associated with expression levels of a nearby coding gene or lincRNA, and numerous variants associated with susceptibility to infectious disease and other phenotypes. We experimentally characterized one candidate nonsynonymous variant in Toll-like receptor 5 (TLR5) and show that it leads to altered NF-κB signaling in response to bacterial flagellin. PAPERFLICK:

Transcriptional and epigenetic dynamics during specification of human embryonic stem cells.

Differentiation of human embryonic stem cells (hESCs) provides a unique opportunity to study the regulatory mechanisms that facilitate cellular transitions in a human context. To that end, we performed comprehensive transcriptional and epigenetic profiling of populations derived through directed differentiation of hESCs representing each of the three embryonic germ layers. Integration of whole-genome bisulfite sequencing, chromatin immunoprecipitation sequencing, and RNA sequencing reveals unique events associated with specification toward each lineage. Lineage-specific dynamic alterations in DNA methylation and H3K4me1 are evident at putative distal regulatory elements that are frequently bound by pluripotency factors in the undifferentiated hESCs. In addition, we identified germ-layer-specific H3K27me3 enrichment at sites exhibiting high DNA methylation in the undifferentiated state. A better understanding of these initial specification events will facilitate identification of deficiencies in current approaches, leading to more faithful differentiation strategies as well as providing insights into the rewiring of human regulatory programs during cellular transitions.

Genome regulation by long noncoding RNAs.

The central dogma of gene expression is that DNA is transcribed into messenger RNAs, which in turn serve as the template for protein synthesis. The discovery of extensive transcription of large RNA transcripts that do not code for proteins, termed long noncoding RNAs (lncRNAs), provides an important new perspective on the centrality of RNA in gene regulation. Here, we discuss genome-scale strategies to discover and characterize lncRNAs. An emerging theme from multiple model systems is that lncRNAs form extensive networks of ribonucleoprotein (RNP) complexes with numerous chromatin regulators and then target these enzymatic activities to appropriate locations in the genome. Consistent with this notion, lncRNAs can function as modular scaffolds to specify higher-order organization in RNP complexes and in chromatin states. The importance of these modes of regulation is underscored by the newly recognized roles of long RNAs for proper gene control across all kingdoms of life.

Group 1 Innate Lymphoid Cell Lineage Identity Is Determined by a cis-Regulatory Element Marked by a Long Non-coding RNA.

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.

Expanded encyclopaedias of DNA elements in the human and mouse genomes.

The human and mouse genomes contain instructions that specify RNAs and proteins and govern the timing, magnitude, and cellular context of their production. To better delineate these elements, phase III of the Encyclopedia of DNA Elements (ENCODE) Project has expanded analysis of the cell and tissue repertoires of RNA transcription, chromatin structure and modification, DNA methylation, chromatin looping, and occupancy by transcription factors and RNA-binding proteins. Here we summarize these efforts, which have produced 5,992 new experimental datasets, including systematic determinations across mouse fetal development. All data are available through the ENCODE data portal (https://www.encodeproject.org), including phase II ENCODE1 and Roadmap Epigenomics2 data. We have developed a registry of 926,535 human and 339,815 mouse candidate cis-regulatory elements, covering 7.9 and 3.4% of their respective genomes, by integrating selected datatypes associated with gene regulation, and constructed a web-based server (SCREEN; http://screen.encodeproject.org) to provide flexible, user-defined access to this resource. Collectively, the ENCODE data and registry provide an expansive resource for the scientific community to build a better understanding of the organization and function of the human and mouse genomes.

Inter-chromosomal Contact Properties in Live-Cell Imaging and in Hi-C.

Imaging (fluorescence in situ hybridization [FISH]) and genome-wide chromosome conformation capture (Hi-C) are two major approaches to the study of higher-order genome organization in the nucleus. Intra-chromosomal and inter-chromosomal interactions (referred to as non-homologous chromosomal contacts [NHCCs]) have been observed by several FISH-based studies, but locus-specific NHCCs have not been detected by Hi-C. Due to crosslinking, neither of these approaches assesses spatiotemporal properties. Toward resolving the discrepancies between imaging and Hi-C, we sought to understand the spatiotemporal properties of NHCCs in living cells by CRISPR/Cas9 live-cell imaging (CLING). In mammalian cells, we find that NHCCs are stable and occur as frequently as intra-chromosomal interactions, but NHCCs occur at farther spatial distance that could explain their lack of detection in Hi-C. By revealing the spatiotemporal properties in living cells, our study provides fundamental insights into the biology of NHCCs.

Massively parallel dissection of RNA in RNA-protein interactions in vivo.

Many of the biological functions performed by RNA are mediated by RNA-binding proteins (RBPs), and understanding the molecular basis of these interactions is fundamental to biology. Here, we present massively parallel RNA assay combined with immunoprecipitation (MPRNA-IP) for in vivo high-throughput dissection of RNA-protein interactions and describe statistical models for identifying RNA domains and parsing the structural contributions of RNA. By using custom pools of tens of thousands of RNA sequences containing systematically designed truncations and mutations, MPRNA-IP is able to identify RNA domains, sequences, and secondary structures necessary and sufficient for protein binding in a single experiment. We show that this approach is successful for multiple RNAs of interest, including the long noncoding RNA NORAD, bacteriophage MS2 RNA, and human telomerase RNA, and we use it to interrogate the hitherto unknown sequence or structural RNA-binding preferences of the DNA-looping factor CTCF. By integrating systematic mutation analysis with crosslinking immunoprecipitation, MPRNA-IP provides a novel high-throughput way to elucidate RNA-based mechanisms behind RNA-protein interactions in vivo.

DNMT1 inhibition by pUG-fold quadruplex RNA.

Aberrant DNA methylation is one of the earliest hallmarks of cancer. DNMT1 is responsible for methylating newly replicated DNA, but the precise regulation of DNMT1 to ensure faithful DNA methylation remains poorly understood. A link between RNA and chromatin-associated proteins has recently emerged, and several studies have shown that DNMT1 can be regulated by a variety of RNAs. In this study, we have confirmed that human DNMT1 indeed interacts with multiple RNAs, including its own nuclear mRNA. Unexpectedly, we found that DNMT1 exhibits a strong and specific affinity for GU-rich RNAs that form a pUG-fold, a noncanonical G-quadruplex. We find that pUG-fold-capable RNAs inhibit DNMT1 activity by inhibiting binding of hemimethylated DNA, and we additionally provide evidence for multiple RNA binding modes with DNMT1. Together, our data indicate that a human chromatin-associated protein binds to and is regulated by pUG-fold RNA.

A large intergenic noncoding RNA induced by p53 mediates global gene repression in the p53 response.

Recently, more than 1000 large intergenic noncoding RNAs (lincRNAs) have been reported. These RNAs are evolutionarily conserved in mammalian genomes and thus presumably function in diverse biological processes. Here, we report the identification of lincRNAs that are regulated by p53. One of these lincRNAs (lincRNA-p21) serves as a repressor in p53-dependent transcriptional responses. Inhibition of lincRNA-p21 affects the expression of hundreds of gene targets enriched for genes normally repressed by p53. The observed transcriptional repression by lincRNA-p21 is mediated through the physical association with hnRNP-K. This interaction is required for proper genomic localization of hnRNP-K at repressed genes and regulation of p53 mediates apoptosis. We propose a model whereby transcription factors activate lincRNAs that serve as key repressors by physically associating with repressive complexes and modulate their localization to sets of previously active genes.

Neat1 is a p53-inducible lincRNA essential for transformation suppression.

The p53 gene is mutated in over half of all cancers, reflecting its critical role as a tumor suppressor. Although p53 is a transcriptional activator that induces myriad target genes, those p53-inducible genes most critical for tumor suppression remain elusive. Here, we leveraged p53 ChIP-seq (chromatin immunoprecipitation [ChIP] combined with high-throughput sequencing) and RNA-seq (RNA sequencing) data sets to identify new p53 target genes, focusing on the noncoding genome. We identify Neat1, a noncoding RNA (ncRNA) constituent of paraspeckles, as a p53 target gene broadly induced by mouse and human p53 in different cell types and by diverse stress signals. Using fibroblasts derived from Neat1-/- mice, we examined the functional role of Neat1 in the p53 pathway. We found that Neat1 is dispensable for cell cycle arrest and apoptosis in response to genotoxic stress. In sharp contrast, Neat1 plays a crucial role in suppressing transformation in response to oncogenic signals. Neat1 deficiency enhances transformation in oncogene-expressing fibroblasts and promotes the development of premalignant pancreatic intraepithelial neoplasias (PanINs) and cystic lesions in KrasG12D-expressing mice. Neat1 loss provokes global changes in gene expression, suggesting a mechanism by which its deficiency promotes neoplasia. Collectively, these findings identify Neat1 as a p53-regulated large intergenic ncRNA (lincRNA) with a key role in suppressing transformation and cancer initiation, providing fundamental new insight into p53-mediated tumor suppression.

Evolutionary divergence of Firre localization and expression.

Long noncoding RNAs (lncRNAs) are rapidly evolving and thus typically poorly conserved in their sequences. How these sequence differences affect the characteristics and potential functions of lncRNAs with shared synteny remains unclear. Here we show that the syntenically conserved lncRNA Firre displays distinct expression and localization patterns in human and mouse. Single molecule RNA FISH reveals that in a range of cell lines, mouse Firre (mFirre) is predominantly nuclear, while human FIRRE (hFIRRE) is distributed between the cytoplasm and nucleus. This localization pattern is maintained in human/mouse hybrid cells expressing both human and mouse Firre, implying that the localization of the lncRNA is species autonomous. We find that the majority of hFIRRE transcripts in the cytoplasm are comprised of isoforms that are enriched in RRD repeats. We furthermore determine that in various tissues, mFirre is more highly expressed than its human counterpart. Our data illustrate that the rapid evolution of syntenic lncRNAs can lead to variations in lncRNA localization and abundance, which in turn may result in disparate lncRNA functions even in closely related species.

"Cat's Cradling" the 3D Genome by the Act of LncRNA Transcription.

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.

An in vivo screen of noncoding loci reveals that Daedalus is a gatekeeper of an Ikaros-dependent checkpoint during haematopoiesis.

Haematopoiesis relies on tightly controlled gene expression patterns as development proceeds through a series of progenitors. While the regulation of hematopoietic development has been well studied, the role of noncoding elements in this critical process is a developing field. In particular, the discovery of new regulators of lymphopoiesis could have important implications for our understanding of the adaptive immune system and disease. Here we elucidate how a noncoding element is capable of regulating a broadly expressed transcription factor, Ikaros, in a lymphoid lineage-specific manner, such that it imbues Ikaros with the ability to specify the lymphoid lineage over alternate fates. Deletion of the Daedalus locus, which is proximal to Ikaros, led to a severe reduction in early lymphoid progenitors, exerting control over the earliest fate decisions during lymphoid lineage commitment. Daedalus locus deletion led to alterations in Ikaros isoform expression and a significant reduction in Ikaros protein. The Daedalus locus may function through direct DNA interaction as Hi-C analysis demonstrated an interaction between the two loci. Finally, we identify an Ikaros-regulated erythroid-lymphoid checkpoint that is governed by Daedalus in a lymphoid-lineage-specific manner. Daedalus appears to act as a gatekeeper of Ikaros's broad lineage-specifying functions, selectively stabilizing Ikaros activity in the lymphoid lineage and permitting diversion to the erythroid fate in its absence. These findings represent a key illustration of how a transcription factor with broad lineage expression must work in concert with noncoding elements to orchestrate hematopoietic lineage commitment.