We recently showed that the mammalian genome encodes >1,000 large intergenic noncoding (linc)RNAs that are clearly conserved across mammals and, thus, functional. Gene expression patterns have implicated these lincRNAs in diverse biological processes, including cell-cycle regulation, immune surveillance, and embryonic stem cell pluripotency. However, the mechanism by which these lincRNAs function is unknown. Here, we expand the catalog of human lincRNAs to Ϸ3,300 by analyzing chromatin-state maps of various human cell types. Inspired by the observation that the well-characterized lincRNA HOTAIR binds the polycomb repressive complex (PRC)2, we tested whether many lincRNAs are physically associated with PRC2. Remarkably, we observe that Ϸ20% of lincRNAs expressed in various cell types are bound by PRC2, and that additional lincRNAs are bound by other chromatin-modifying complexes. Also, we show that siRNAmediated depletion of certain lincRNAs associated with PRC2 leads to changes in gene expression, and that the up-regulated genes are enriched for those normally silenced by PRC2. We propose a model in which some lincRNAs guide chromatin-modifying complexes to specific genomic loci to regulate gene expression.histone modifications ͉ epigenetic regulation ͉ polycomb

Gene expression is a fundamentally stochastic process, with randomness in transcription and translation leading to significant cell-to-cell variations in mRNA and protein levels. This variation appears in organisms ranging from microbes to metazoans and its characteristics depend both on the biophysical parameters governing gene expression and on gene network structure. Stochastic gene expression can have important consequences for cellular function, being beneficial in some contexts and harmful in others. These situations include stress response, pathogenesis, metabolism, development, cell cycle, circadian rhythms and aging.

Individual cells in genetically homogeneous populations have been found to express different numbers of molecules of specific proteins. We investigated the origins of these variations in mammalian cells by counting individual molecules of mRNA produced from a reporter gene that was stably integrated into the cell's genome. We found that there are massive variations in the number of mRNA molecules present in each cell. These variations occur because mRNAs are synthesized in short but intense bursts of transcription beginning when the gene transitions from an inactive to an active state and ending when they transition back to the inactive state. We show that these transitions are intrinsically random and not due to global, extrinsic factors such as the levels of transcriptional activators. Moreover, the gene activation causes burst-like expression of all genes within a wider genomic locus. We further found that bursts are also exhibited in the synthesis of natural genes. The bursts of mRNA expression can be buffered at the protein level by slow protein degradation rates. A stochastic model of gene activation and inactivation was developed to explain the statistical properties of the bursts. The model showed that increasing the level of transcription factors increases the average size of the bursts rather than their frequency. These results demonstrate that gene expression in mammalian cells is subject to large, intrinsically random fluctuations and raise questions about how cells are able to function in the face of such noise.

We describe a method for imaging individual mRNA molecules in fixed cells by probing each mRNA species with 48 or more short, singly labeled oligonucleotide probes. This makes each mRNA molecule visible as a computationally identifiable fluorescent spot via fluorescence microscopy. We demonstrate simultaneous detection of three mRNA species in single cells and mRNA detection in yeast, nematodes, fruit fly wing discs, mammalian cell lines and neurons.

Therapies targeting signaling molecules mutated in cancers can often have striking short-term effects, but the emergence of resistant cancer cells is a major barrier to full cures 1,2 . Resistance can result from a secondary mutations 3,4 , but other times there is no clear genetic cause, raising the possibility of non-genetic rare cell variability [5][6][7][8][9][10][11] . Here, we show that melanoma cells can display profound transcriptional variability at the single cell level that predicts which cells will ultimately resist drug treatment. This variability involves infrequent, semi-coordinated transcription of a number of resistance markers at high levels in a very small percentage of cells. The addition of Users may view, print, copy, and download text and data-mine the content in such documents, for the purposes of academic research, subject always to the full Conditions of use: http://www.nature.com/authors/editorial_policies/license.html#terms Author contributions: SMS, AR designed the study. SMS performed all experiments and analysis except: MD, ST assisted with fluctuation analysis and RNA-sequencing; EAT, BE performed NGFR and AXL sort experiments; CK, MB, KS performed PDX experiments; PB, MH provided cell lines; MX performed WM989-A6 characterization; EE developed iterative RNA FISH protocol; INA, KN performed DNA sequencing. MH provided guidance. SMS, AR wrote the paper. Author information:AR receives consulting income and AR and SMS receive royalties related to Stellaris™ RNA FISH probes.

Several of the thousands of human long non-coding RNAs (lncRNAs) have been functionally characterized1–4; however, potential roles for lncRNAs in somatic tissue differentiation remain poorly understood. Here we show that a 3.7-kilobase lncRNA, terminal differentiation-induced ncRNA (TINCR), controls human epidermal differentiation by a post-transcriptional mechanism. TINCR is required for high messenger RNA abundance of key differentiation genes, many of which are mutated in human skin diseases, including FLG, LOR, ALOXE3, ALOX12B, ABCA12, CASP14 and ELOVL3. TINCR-deficient epidermis lacked terminal differentiation ultrastructure, including keratohyalin granules and intact lamellar bodies. Genome-scale RNA interactome analysis revealed that TINCR interacts with a range of differentiation mRNAs. TINCR–mRNA interaction occurs through a 25-nucleotide ‘TINCR box’ motif that is strongly enriched in interacting mRNAs and required for TINCR binding. A high-throughput screen to analyse TINCR binding capacity to approximately 9,400 human recombinant proteins revealed direct binding of TINCR RNA to the staufen1 (STAU1) protein. STAU1-deficient tissue recapitulated the impaired differentiation seen with TINCR depletion. Loss of UPF1 and UPF2, both of which are required for STAU1-mediated RNA decay, however, did not have differentiation effects. Instead, the TINCR–STAU1 complex seems to mediate stabilization of differentiation mRNAs, such as KRT80. These data identify TINCR as a key lncRNA required for somatic tissue differentiation, which occurs through lncRNA binding to differentiation mRNAs to ensure their expression.

The phenotypic differences between individual organisms can often be ascribed to underlying genetic and environmental variation. However, even genetically identical organisms in homogenous environments vary, suggesting that randomness in developmental processes such as gene expression may also generate diversity. In order to examine the consequences of gene expression variability in multicellular organisms, we studied intestinal specification in the roundworm Caenorhabditis elegans in which wild-type cell fate is invariant and controlled by a small transcriptional network. Mutations in elements of this network can have indeterminate effects: some mutant embryos fail to develop intestinal cells, while others produce intestinal precursors. By counting transcripts of the genes in this network in individual embryos, we show that the expression of an otherwise redundant gene becomes highly variable in the mutants and that this variation is thresholded to produce an ON/OFF expression pattern of the master regulatory gene of intestinal differentiation. Our results demonstrate that mutations in developmental networks can expose otherwise buffered stochastic variability in gene expression, leading to pronounced phenotypic variation.

Random cell-to-cell variations in gene expression within an isogenic population can lead to transitions between alternative states of gene expression. Little is known about how these variations (noise) in natural systems affect such transitions. In Bacillus subtilis, noise in ComK, the protein that regulates competence for DNA uptake, is thought to cause cells to transition to the competent state in which genes encoding DNA uptake proteins are expressed. We demonstrate that noise in comK expression selects cells for competence and that experimental reduction of this noise decreases the number of competent cells. We also show that transitions are limited temporally by a reduction in comK transcription. These results illustrate how such stochastic transitions are regulated in a natural system and suggest that noise characteristics are subject to evolutionary forces.Variability in gene expression within a population of genetically identical cells enables those cells to maintain a diversity of phenotypes, potentially enhancing fitness (1, 2). When the underlying gene network contains regulatory positive feedback loops, individual cells can exist in different states; some cells may, for example, live in the "off" expression state of a particular gene, whereas others are in the "on" expression state (this is an example of bistable gene expression). These stochastic fluctuations in gene expression, commonly referred to as noise, have been proposed to cause transitions between these states (3-7). We apply recently developed theories of noise (8, 9) to examine how noise influences these transitions in a natural system. An example of bistable expression with associated stochastic transitions (10-16) involves the ability of the soil bacterium Bacillus subtilis to develop "competence" for DNA uptake as it enters stationary growth phase, potentially allowing bacteria to increase their fitness by incorporating new genetic material. The genes needed for competence are transcribed only in the presence of ComK, the master regulator of competence. comK expression is subject to positive autoregulation effected by the cooperative binding of ComK to its own promoter ( Fig. 1A) (21), effectively lowering the rate of ComK degradation and allowing random fluctuations in the level of ComK to occasionally cause transitions to the competent state. Cells continue to randomly transition to competence for 2 hours, by which time (T 2 ) transitions have ceased to occur (16) and the 15% of the cells that have become competent remain so until diluted into fresh growth medium ( Fig. 1C and movie S1). In this report, we ask why cells only transition to competence for a limited duration of time and investigate the source of the fluctuations that actuate the ComK feedback loop in a minority of cells.To understand why cells only transition to the competent state for ~2 hours during stationary phase, we examined the dynamics of comK expression in noncompetent cells. Because the level of ComK in noncompetent cells is very low, we used fluoresce...