Dynamics and Spatial Genomics of the Nascent Transcriptome by Intron seqFISH.

Visualization of the transcriptome and the nuclear organization in situ has been challenging for single-cell analysis. Here, we demonstrate a multiplexed single-molecule in situ method, intron seqFISH, that allows imaging of 10,421 genes at their nascent transcription active sites in single cells, followed by mRNA and lncRNA seqFISH and immunofluorescence. This nascent transcriptome-profiling method can identify different cell types and states with mouse embryonic stem cells and fibroblasts. The nascent sites of RNA synthesis tend to be localized on the surfaces of chromosome territories, and their organization in individual cells is highly variable. Surprisingly, the global nascent transcription oscillated asynchronously in individual cells with a period of 2 hr in mouse embryonic stem cells, as well as in fibroblasts. Together, spatial genomics of the nascent transcriptome by intron seqFISH reveals nuclear organizational principles and fast dynamics in single cells that are otherwise obscured.

Shared retinoic acid responsive enhancers coordinately regulate nascent transcription of Hoxb coding and non-coding RNAs in the developing mouse neural tube.

Signaling pathways regulate the patterns of Hox gene expression that underlie their functions in the specification of axial identity. Little is known about the properties of cis-regulatory elements and underlying transcriptional mechanisms that integrate graded signaling inputs to coordinately control Hox expression. Here, we optimized a single molecule fluorescent in situ hybridization (smFISH) technique with probes spanning introns to evaluate how three shared retinoic acid response element (RARE)-dependent enhancers in the Hoxb cluster regulate patterns of nascent transcription in vivo at the level of single cells in wild-type and mutant embryos. We predominately detect nascent transcription of only a single Hoxb gene in each cell, with no evidence for simultaneous co-transcriptional coupling of all or specific subsets of genes. Single and/or compound RARE mutations indicate that each enhancer differentially impacts global and local patterns of nascent transcription, suggesting that selectivity and competitive interactions between these enhancers is important to robustly maintain the proper levels and patterns of nascent Hoxb transcription. This implies that rapid and dynamic regulatory interactions potentiate transcription of genes through combined inputs from these enhancers in coordinating the retinoic acid response.

A Growing Toolbox to Image Gene Expression in Single Cells: Sensitive Approaches for Demanding Challenges.

The spatiotemporal regulation of gene expression is key to many biological processes. Recent imaging approaches opened exciting perspectives for understanding the intricate mechanisms regulating RNA metabolism, from synthesis to decay. Imaging techniques allow their observation at high spatial and temporal resolution, while keeping cellular morphology and micro-environment intact. Here, we focus on approaches for imaging single RNA molecules in cells, tissues, and embryos. In fixed cells, the rapid development of smFISH multiplexing opens the way to large-scale single-molecule studies, while in live cells, gene expression can be observed in real time in its native context. We highlight the strengths and limitations of these methods, as well as future challenges. We present how they advanced our understanding of gene expression heterogeneity and bursting, as well as the spatiotemporal aspects of splicing, translation, and RNA decay. These insights yield a dynamic and stochastic view of gene expression in single cells.

Spatial Organization of Single mRNPs at Different Stages of the Gene Expression Pathway.

mRNAs form ribonucleoprotein complexes (mRNPs) by association with proteins that are crucial for mRNA metabolism. While the mRNP proteome has been well characterized, little is known about mRNP organization. Using a single-molecule approach, we show that mRNA conformation changes depending on its cellular localization and translational state. Compared to nuclear mRNPs and lncRNPs, association with ribosomes decompacts individual mRNAs, while pharmacologically dissociating ribosomes or sequestering them into stress granules leads to increased compaction. Moreover, translating mRNAs rarely show co-localized 5' and 3' ends, indicating either that mRNAs are not translated in a closed-loop configuration, or that mRNA circularization is transient, suggesting that a stable closed-loop conformation is not a universal state for all translating mRNAs.

Single-Molecule Analysis Reveals Linked Cycles of RSC Chromatin Remodeling and Ace1p Transcription Factor Binding in Yeast.

It is unknown how the dynamic binding of transcription factors (TFs) is molecularly linked to chromatin remodeling and transcription. Using single-molecule tracking (SMT), we show that the chromatin remodeler RSC speeds up the search process of the TF Ace1p for its response elements (REs) at the CUP1 promoter. We quantified smFISH mRNA data using a gene bursting model and demonstrated that RSC regulates transcription bursts of CUP1 only by modulating TF occupancy but does not affect initiation and elongation rates. We show by SMT that RSC binds to activated promoters transiently, and based on MNase-seq data, that RSC does not affect the nucleosomal occupancy at CUP1. Therefore, transient binding of Ace1p and rapid bursts of transcription at CUP1 may be dependent on short repetitive cycles of nucleosome mobilization. This type of regulation reduces the transcriptional noise and ensures a homogeneous response of the cell population to heavy metal stress.

Pairwise biosynthesis of ion channels stabilizes excitability and mitigates arrhythmias.

Coordinated expression of ion channels is crucial for cardiac rhythms, neural signaling, and cell cycle progression. Perturbation of this balance results in many disorders including cardiac arrhythmias. Prior work revealed association of mRNAs encoding cardiac NaV1.5 (SCN5A) and hERG1 (KCNH2), but the functional significance of this association was not established. Here, we provide a more comprehensive picture of KCNH2, SCN5A, CACNA1C, and KCNQ1 transcripts collectively copurifying with nascent hERG1, NaV1.5, CaV1.2, or KCNQ1 channel proteins. Single-molecule fluorescence in situ hybridization (smFISH) combined with immunofluorescence reveals that the channel proteins are synthesized predominantly as heterotypic pairs from discrete molecules of mRNA, not as larger cotranslational complexes. Puromycin disrupted colocalization of mRNA with its encoded protein, as expected, but remarkably also pairwise mRNA association, suggesting that transcript association relies on intact translational machinery or the presence of the nascent protein. Targeted depletion of KCHN2 by specific shRNA resulted in concomitant reduction of all associated mRNAs, with a corresponding reduction in the encoded channel currents. This co-knockdown effect, originally described for KCNH2 and SCN5A, thus appears to be a general phenomenon among transcripts encoding functionally related proteins. In multielectrode array recordings, proarrhythmic behavior arose when IKr was reduced by the selective blocker dofetilide at IC50 concentrations, but not when equivalent reductions were mediated by shRNA, suggesting that co-knockdown mitigates proarrhythmic behavior expected from the selective reduction of a single channel species. We propose that coordinated, cotranslational association of functionally related ion channel mRNAs confers electrical stability by co-regulating complementary ion channels in macromolecular complexes.

LST-1 is a bifunctional regulator that feeds back on Notch-dependent transcription to regulate C. elegans germline stem cells.

Notch signaling regulates stem cells across animal phylogeny. C. elegans Notch signaling activates transcription of two genes, lst-1 and sygl-1, that encode potent regulators of germline stem cells. The LST-1 protein regulates stem cells in two distinct ways: It promotes self-renewal posttranscriptionally and also restricts self-renewal by a poorly understood mechanism. Its self-renewal promoting activity resides in its N-terminal region, while its self-renewal restricting activity resides in its C-terminal region and requires the Zn finger. Here, we report that LST-1 limits self-renewal by down-regulating Notch-dependent transcription. We detect LST-1 in the nucleus, in addition to its previously known cytoplasmic localization. LST-1 lowers nascent transcript levels at both lst-1 and sygl-1 loci but not at let-858, a Notch-independent locus. LST-1 also lowers levels of two key components of the Notch activation complex, the LAG-1 DNA binding protein and Notch intracellular domain (NICD). Genetically, an LST-1 Zn finger mutant increases Notch signaling strength in both gain- and loss-of-function GLP-1/Notch receptor mutants. Biochemically, LST-1 co-immunoprecipitates with LAG-1 from nematode extracts, suggesting a direct effect. LST-1 is thus a bifunctional regulator that coordinates posttranscriptional and transcriptional mechanisms in a single protein. This LST-1 bifunctionality relies on its bipartite protein architecture and is bolstered by generation of two LST-1 isoforms, one specialized for Notch downregulation. A conserved theme from worms to human is the coupling of PUF-mediated RNA repression together with Notch feedback in the same protein.

Notch-dependent DNA cis-regulatory elements and their dose-dependent control of C. elegans stem cell self-renewal.

A long-standing biological question is how DNA cis-regulatory elements shape transcriptional patterns during metazoan development. Reporter constructs, cell culture assays and computational modeling have made major contributions to answering this question, but analysis of elements in their natural context is an important complement. Here, we mutate Notch-dependent LAG-1 binding sites (LBSs) in the endogenous Caenorhabditis elegans sygl-1 gene, which encodes a key stem cell regulator, and analyze the consequences on sygl-1 expression (nascent transcripts, mRNA, protein) and stem cell maintenance. Mutation of one LBS in a three-element cluster approximately halved both expression and stem cell pool size, whereas mutation of two LBSs essentially abolished them. Heterozygous LBS mutant clusters provided intermediate values. Our results lead to two major conclusions. First, both LBS number and configuration impact cluster activity: LBSs act additively in trans and synergistically in cis. Second, the SYGL-1 gradient promotes self-renewal above its functional threshold and triggers differentiation below the threshold. Our approach of coupling CRISPR/Cas9 LBS mutations with effects on both molecular and biological readouts establishes a powerful model for in vivo analyses of DNA cis-regulatory elements.

PRD-2 mediates clock-regulated perinuclear localization of clock gene RNAs within the circadian cycle of Neurospora.

The transcription-translation negative feedback loops underlying animal and fungal circadian clocks are remarkably similar in their molecular regulatory architecture and, although much is understood about their central mechanism, little is known about the spatiotemporal dynamics of the gene products involved. A common feature of these circadian oscillators is a significant temporal delay between rhythmic accumulation of clock messenger RNAs (mRNAs) encoding negative arm proteins, for example, frq in Neurospora and Per1-3 in mammals, and the appearance of the clock protein complexes assembled from the proteins they encode. Here, we report use of single-molecule RNA fluorescence in situ hybridization (smFISH) to show that the fraction of nuclei actively transcribing the clock gene frq changes in a circadian manner, and that these mRNAs cycle in abundance with fewer than five transcripts per nucleus at any time. Spatial point patterning statistics reveal that frq is spatially clustered near nuclei in a time of day-dependent manner and that clustering requires an RNA-binding protein, PRD-2 (PERIOD-2), recently shown also to bind to mRNA encoding another core clock component, casein kinase 1. An intrinsically disordered protein, PRD-2 displays behavior in vivo and in vitro consistent with participation in biomolecular condensates. These data are consistent with a role for phase-separating RNA-binding proteins in spatiotemporally organizing clock mRNAs to facilitate local translation and assembly of clock protein complexes.

Uniform gene expression in embryos is achieved by temporal averaging of transcription noise.

Transcription is often stochastic. This is seemingly incompatible with the importance of gene expression during development. Here we show that during zebrafish embryogenesis, transcription activation is stochastic due to (1) genes acquiring transcriptional competence at different times in different cells, (2) differences in cell cycle stage between cells, and (3) the stochastic nature of transcription. Initially, stochastic transcription causes large cell-to-cell differences in transcript levels. However, variability is reduced by lengthening cell cycles and the accumulation of transcription events in each cell. Temporal averaging might provide a general context in which to understand how embryos deal with stochastic transcription.

FISH-quant v2: a scalable and modular tool for smFISH image analysis.

Regulation of RNA abundance and localization is a key step in gene expression control. Single-molecule RNA fluorescence in situ hybridization (smFISH) is a widely used single-cell-single-molecule imaging technique enabling quantitative studies of gene expression and its regulatory mechanisms. Today, these methods are applicable at a large scale, which in turn come with a need for adequate tools for data analysis and exploration. Here, we present FISH-quant v2, a highly modular tool accessible for both experts and non-experts. Our user-friendly package allows the user to segment nuclei and cells, detect isolated RNAs, decompose dense RNA clusters, quantify RNA localization patterns and visualize these results both at the single-cell level and variations within the cell population. This tool was validated and applied on large-scale smFISH image data sets, revealing diverse subcellular RNA localization patterns and a surprisingly high degree of cell-to-cell heterogeneity.

Quantitative intracellular retention of delivered RNAs through optimized cell fixation and immunostaining.

Detection of nucleic acids within subcellular compartments is key to understanding their function. Determining the intracellular distribution of nucleic acids requires quantitative retention and estimation of their association with different organelles by immunofluorescence microscopy. This is particularly important for the delivery of nucleic acid therapeutics, which depends on endocytic uptake and endosomal escape. However, the current protocols fail to preserve the majority of exogenously delivered nucleic acids in the cytoplasm. To solve this problem, by monitoring Cy5-labeled mRNA delivered to primary human adipocytes via lipid nanoparticles (LNP), we optimized cell fixation, permeabilization, and immunostaining of a number of organelle markers, achieving quantitative retention of mRNA and allowing visualization of levels that escape detection using conventional procedures. The optimized protocol proved effective on exogenously delivered siRNA, miRNA, as well as endogenous miRNA. Our protocol is compatible with RNA probes of single molecule fluorescence in situ hybridization (smFISH) and molecular beacon, thus demonstrating that it is broadly applicable to study a variety of nucleic acids in cultured cells.

Visualizing active viral infection reveals diverse cell fates in synchronized algal bloom demise.

Marine viruses are the most abundant biological entity in the ocean and are considered as major evolutionary drivers of microbial life [C. A. Suttle, Nat. Rev. Microbiol. 5, 801-812 (2007)]. Yet, we lack quantitative approaches to assess their impact on the marine ecosystem. Here, we provide quantification of active viral infection in the bloom forming single-celled phytoplankton Emiliania huxleyi infected by the large virus EhV, using high-throughput single-molecule messenger RNA in situ hybridization (smFISH) of both virus and host transcripts. In natural samples, viral infection reached only 25% of the population despite synchronized bloom demise exposing the coexistence of infected and noninfected subpopulations. We prove that photosynthetically active cells chronically release viral particles through nonlytic infection and that viral-induced cell lysis can occur without viral release, thus challenging major assumptions regarding the life cycle of giant viruses. We could also assess active infection in cell aggregates linking viral infection and carbon export to the deep ocean [C. P. Laber et al., Nat. Microbiol. 3, 537-547 (2018)] and suggest a potential host defense strategy by enrichment of infected cells in sinking aggregates. Our approach can be applied to diverse marine microbial systems, opening a mechanistic dimension to the study of biotic interactions in the ocean.

What Is a Transcriptional Burst?

The idea that gene activity can be discontinuous will not surprise many biologists - many genes are restricted in when and where they can be expressed. Yet during the past 15 years, a collection of observations compiled under the umbrella term 'transcriptional bursting' has received considerable interest. Direct visualization of the dynamics of discontinuous transcription has expanded our understanding of basic transcriptional mechanisms and their regulation and provides a real-time readout of gene activity during the life of a cell. In this review, we try to reconcile the different views of the transcriptional process emerging from studies of bursting, and how this work contextualizes the relative importance of different regulatory inputs to normal dynamic ranges of gene activity.

Reference genes for quantitative Arabidopsis single molecule RNA fluorescence in situ hybridization.

Subcellular mRNA quantities and spatial distributions are fundamental for driving gene regulatory programmes. Single molecule RNA fluorescence in situ hybridization (smFISH) uses fluorescent probes to label individual mRNA molecules, thereby facilitating both localization and quantitative studies. Validated reference mRNAs function as positive controls and are required for calibration. Here we present selection criteria for the first set of Arabidopsis smFISH reference genes. Following sequence and transcript data assessments, four mRNA probe sets were selected for imaging. Transcript counts per cell, correlations with cell size, and corrected fluorescence intensities were all calculated for comparison. In addition to validating reference probe sets, we present sample preparation steps that can retain green fluorescent protein fluorescence, thereby providing a method for simultaneous RNA and protein detection. In summary, our reference gene analyses, modified protocol, and simplified quantification method together provide a firm foundation for future quantitative single molecule RNA studies in Arabidopsis root apical meristem cells.

The circadian oscillator analysed at the single-transcript level.

The circadian clock is an endogenous and self-sustained oscillator that anticipates daily environmental cycles. While rhythmic gene expression of circadian genes is well-described in populations of cells, the single-cell mRNA dynamics of multiple core clock genes remain largely unknown. Here we use single-molecule fluorescence in situ hybridisation (smFISH) at multiple time points to measure pairs of core clock transcripts, Rev-erbα (Nr1d1), Cry1 and Bmal1, in mouse fibroblasts. The mean mRNA level oscillates over 24 h for all three genes, but mRNA numbers show considerable spread between cells. We develop a probabilistic model for multivariate mRNA counts using mixtures of negative binomials, which accounts for transcriptional bursting, circadian time and cell-to-cell heterogeneity, notably in cell size. Decomposing the mRNA variability into distinct noise sources shows that clock time contributes a small fraction of the total variability in mRNA number between cells. Thus, our results highlight the intrinsic biological challenges in estimating circadian phase from single-cell mRNA counts and suggest that circadian phase in single cells is encoded post-transcriptionally.

24-nt phasiRNAs move from tapetal to meiotic cells in maize anthers.

In maize, 24-nt phased, secondary small interfering RNAs (phasiRNAs) are abundant in meiotic stage anthers, but their distribution and functions are not precisely known. Using laser capture microdissection, we analyzed tapetal cells, meiocytes and other somatic cells at several stages of anther development to establish the timing of 24-PHAS precursor transcripts and the 24-nt phasiRNA products. By integrating RNA and small RNA profiling plus single-molecule and small RNA FISH (smFISH or sRNA-FISH) spatial detection, we demonstrate that the tapetum is the primary site of 24-PHAS precursor and Dcl5 transcripts and the resulting 24-nt phasiRNAs. Interestingly, 24-nt phasiRNAs accumulate in all cell types, with the highest levels in meiocytes, followed by tapetum. Our data support the conclusion that 24-nt phasiRNAs are mobile from tapetum to meiocytes and to other somatic cells. We discuss possible roles for 24-nt phasiRNAs in anther cell types.

The spatial and temporal dynamics of nuclear RNAi-targeted retrotransposon transcripts in Caenorhabditis elegans.

Nuclear RNA interference provides a unique approach to the study of RNA-mediated transgenerational epigenetic inheritance. A paradox in the field is that expression of target loci is necessary for the initiation and maintenance of their silencing. How expression and repression are coordinated during animal development is poorly understood. To resolve this gap, we took imaging, deep-sequencing and genetic approaches towards delineating the developmental regulation and subcellular localization of RNA transcripts of two representative endogenous targets, the LTR retrotransposons Cer3 and Cer8. By examining wild-type worms and a collection of mutant strains, we found that the expression and silencing cycle of Cer3 and Cer8 is coupled with embryonic and germline development. Strikingly, endogenous targets exhibit a hallmark of nuclear enrichment of their RNA transcripts. In addition, germline and somatic repressions of Cer3 have different genetic requirements for three heterochromatin enzymes, MET-2, SET-25 and SET-32, in conjunction with the nuclear Argonaute protein HRDE-1. These results provide the first comprehensive cellular and developmental characterization of nuclear RNAi activities throughout the animal reproductive cycle.

Cell Atlas technologies and insights into tissue architecture.

Since Robert Hooke first described the existence of 'cells' in 1665, scientists have sought to identify and further characterise these fundamental units of life. While our understanding of cell location, morphology and function has expanded greatly; our understanding of cell types and states at the molecular level, and how these function within tissue architecture, is still limited. A greater understanding of our cells could revolutionise basic biology and medicine. Atlasing initiatives like the Human Cell Atlas aim to identify all cell types at the molecular level, including their physical locations, and to make this reference data openly available to the scientific community. This is made possible by a recent technology revolution: both in single-cell molecular profiling, particularly single-cell RNA sequencing, and in spatially resolved methods for assessing gene and protein expression. Here, we review available and upcoming atlasing technologies, the biological insights gained to date and the promise of this field for the future.

Imaging mRNA and protein interactions within neurons.

RNA-protein interactions are essential for proper gene expression regulation, particularly in neurons with unique spatial constraints. Currently, these interactions are defined biochemically, but a method is needed to evaluate them quantitatively within morphological context. Colocalization of two-color labels using wide-field microscopy is a method to infer these interactions. However, because of chromatic aberrations in the objective lens, this approach lacks the resolution to determine whether two molecules are physically in contact or simply nearby by chance. Here, we developed a robust super registration methodology that corrected the chromatic aberration across the entire image field to within 10 nm, which is capable of determining whether two molecules are physically interacting or simply in proximity by random chance. We applied this approach to image single-molecule FISH in combination with immunofluorescence (smFISH-IF) and determined whether the association between an mRNA and binding protein(s) within a neuron was significant or accidental. We evaluated several mRNA-binding proteins identified from RNA pulldown assays to determine which of these exhibit bona fide interactions. Surprisingly, many known mRNA-binding proteins did not bind the mRNA in situ, indicating that adventitious interactions are significant using existing technology. This method provides an ability to evaluate two-color registration compatible with the scale of molecular interactions.