Charting the cellular biogeography in colitis reveals fibroblast trajectories and coordinated spatial remodeling.

Gut inflammation involves contributions from immune and non-immune cells, whose interactions are shaped by the spatial organization of the healthy gut and its remodeling during inflammation. The crosstalk between fibroblasts and immune cells is an important axis in this process, but our understanding has been challenged by incomplete cell-type definition and biogeography. To address this challenge, we used multiplexed error-robust fluorescence in situ hybridization (MERFISH) to profile the expression of 940 genes in 1.35 million cells imaged across the onset and recovery from a mouse colitis model. We identified diverse cell populations, charted their spatial organization, and revealed their polarization or recruitment in inflammation. We found a staged progression of inflammation-associated tissue neighborhoods defined, in part, by multiple inflammation-associated fibroblasts, with unique expression profiles, spatial localization, cell-cell interactions, and healthy fibroblast origins. Similar signatures in ulcerative colitis suggest conserved human processes. Broadly, we provide a framework for understanding inflammation-induced remodeling in the gut and other tissues.

Omics goes spatial epigenomics.

Spatial omics techniques generate spatially resolved, comprehensive data about molecules that define the identity and function of cells in tissues. Epigenetic multiplexing approaches such as Multiplexed Error-robust FISH (MERFISH), introduced by Lu et al.1 in this issue of Cell, now allows researchers to study the epigenomic regulation of gene expression in a tissue-region specific manner.

Genome-Scale Imaging of the 3D Organization and Transcriptional Activity of Chromatin.

The 3D organization of chromatin regulates many genome functions. Our understanding of 3D genome organization requires tools to directly visualize chromatin conformation in its native context. Here we report an imaging technology for visualizing chromatin organization across multiple scales in single cells with high genomic throughput. First we demonstrate multiplexed imaging of hundreds of genomic loci by sequential hybridization, which allows high-resolution conformation tracing of whole chromosomes. Next we report a multiplexed error-robust fluorescence in situ hybridization (MERFISH)-based method for genome-scale chromatin tracing and demonstrate simultaneous imaging of more than 1,000 genomic loci and nascent transcripts of more than 1,000 genes together with landmark nuclear structures. Using this technology, we characterize chromatin domains, compartments, and trans-chromosomal interactions and their relationship to transcription in single cells. We envision broad application of this high-throughput, multi-scale, and multi-modal imaging technology, which provides an integrated view of chromatin organization in its native structural and functional context.

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.

Identification of Gene Positioning Factors Using High-Throughput Imaging Mapping.

Genomes are arranged non-randomly in the 3D space of the cell nucleus. Here, we have developed HIPMap, a high-precision, high-throughput, automated fluorescent in situ hybridization imaging pipeline, for mapping of the spatial location of genome regions at large scale. High-throughput imaging position mapping (HIPMap) enabled an unbiased siRNA screen for factors involved in genome organization in human cells. We identify 50 cellular factors required for proper positioning of a set of functionally diverse genomic loci. Positioning factors include chromatin remodelers, histone modifiers, and nuclear envelope and pore proteins. Components of the replication and post-replication chromatin re-assembly machinery are prominently represented among positioning factors, and timely progression of cells through replication, but not mitosis, is required for correct gene positioning. Our results establish a method for the large-scale mapping of genome locations and have led to the identification of a compendium of cellular factors involved in spatial genome organization.

Molecularly defined and spatially resolved cell atlas of the whole mouse brain.

In mammalian brains, millions to billions of cells form complex interaction networks to enable a wide range of functions. The enormous diversity and intricate organization of cells have impeded our understanding of the molecular and cellular basis of brain function. Recent advances in spatially resolved single-cell transcriptomics have enabled systematic mapping of the spatial organization of molecularly defined cell types in complex tissues1-3, including several brain regions (for example, refs. 1-11). However, a comprehensive cell atlas of the whole brain is still missing. Here we imaged a panel of more than 1,100 genes in approximately 10 million cells across the entire adult mouse brains using multiplexed error-robust fluorescence in situ hybridization12 and performed spatially resolved, single-cell expression profiling at the whole-transcriptome scale by integrating multiplexed error-robust fluorescence in situ hybridization and single-cell RNA sequencing data. Using this approach, we generated a comprehensive cell atlas of more than 5,000 transcriptionally distinct cell clusters, belonging to more than 300 major cell types, in the whole mouse brain with high molecular and spatial resolution. Registration of this atlas to the mouse brain common coordinate framework allowed systematic quantifications of the cell-type composition and organization in individual brain regions. We further identified spatial modules characterized by distinct cell-type compositions and spatial gradients featuring gradual changes of cells. Finally, this high-resolution spatial map of cells, each with a transcriptome-wide expression profile, allowed us to infer cell-type-specific interactions between hundreds of cell-type pairs and predict molecular (ligand-receptor) basis and functional implications of these cell-cell interactions. These results provide rich insights into the molecular and cellular architecture of the brain and a foundation for functional investigations of neural circuits and their dysfunction in health and disease.

Deconwolf enables high-performance deconvolution of widefield fluorescence microscopy images.

Microscopy-based spatially resolved omic methods are transforming the life sciences. However, these methods rely on high numerical aperture objectives and cannot resolve crowded molecular targets, limiting the amount of extractable biological information. To overcome these limitations, here we develop Deconwolf, an open-source, user-friendly software for high-performance deconvolution of widefield fluorescence microscopy images, which efficiently runs on laptop computers. Deconwolf enables accurate quantification of crowded diffraction limited fluorescence dots in DNA and RNA fluorescence in situ hybridization images and allows robust detection of individual transcripts in tissue sections imaged with ×20 air objectives. Deconvolution of in situ spatial transcriptomics images with Deconwolf increased the number of transcripts identified more than threefold, while the application of Deconwolf to images obtained by fluorescence in situ sequencing of barcoded Oligopaint probes drastically improved chromosome tracing. Deconwolf greatly facilitates the use of deconvolution in many bioimaging applications.

Optimality in the development of intestinal crypts.

Intestinal crypts in mammals are comprised of long-lived stem cells and shorter-lived progenies. These two populations are maintained in specific proportions during adult life. Here, we investigate the design principles governing the dynamics of these proportions during crypt morphogenesis. Using optimal control theory, we show that a proliferation strategy known as a "bang-bang" control minimizes the time to obtain a mature crypt. This strategy consists of a surge of symmetric stem cell divisions, establishing the entire stem cell pool first, followed by a sharp transition to strictly asymmetric stem cell divisions, producing nonstem cells with a delay. We validate these predictions using lineage tracing and single-molecule fluorescence in situ hybridization of intestinal crypts in infant mice, uncovering small crypts that are entirely composed of Lgr5-labeled stem cells, which become a minority as crypts continue to grow. Our approach can be used to uncover similar design principles in other developmental systems.

Imaging the Architecture of Granulomas Induced by Mycobacterium tuberculosis Infection with Single-molecule Fluorescence In Situ Hybridization.

Granulomas are an important hallmark of Mycobacterium tuberculosis infection. They are organized and dynamic structures created when immune cells assemble around the sites of infection in the lungs that locally restrict M. tuberculosis growth and the host's inflammatory responses. The cellular architecture of granulomas is traditionally studied by immunofluorescence labeling of surface markers on the host cells. However, very few Abs are available for model animals used in tuberculosis research, such as nonhuman primates and rabbits, and secreted immunological markers such as cytokines cannot be imaged in situ using Abs. Furthermore, traditional phenotypic surface markers do not provide sufficient resolution for the detection of the many subtypes and differentiation states of immune cells. Using single-molecule fluorescence in situ hybridization (smFISH) and its derivatives, amplified smFISH and iterative smFISH, we developed a platform for imaging mRNAs encoding immune markers in rabbit and macaque tuberculosis granulomas. Multiplexed imaging for several mRNA and protein markers was followed by quantitative measurement of the expression of these markers in single cells. An analysis of the combinatorial expressions of these markers allowed us to classify the cells into several subtypes, and to chart their densities within granulomas. For one mRNA target, hypoxia-inducible factor-1α, we imaged its mRNA and protein in the same cells, demonstrating the specificity of the probes. This method paves the way for defining granular differentiation states and cell subtypes from transcriptomic data, identifying key mRNA markers for these cell subtypes, and then locating the cells in the spatial context of granulomas.

Highly multiplexed spatial mapping of microbial communities.

Mapping the complex biogeography of microbial communities in situ with high taxonomic and spatial resolution poses a major challenge because of the high density1 and rich diversity2 of species in environmental microbiomes and the limitations of optical imaging technology3-6. Here we introduce high-phylogenetic-resolution microbiome mapping by fluorescence in situ hybridization (HiPR-FISH), a versatile technology that uses binary encoding, spectral imaging and decoding based on machine learning to create micrometre-scale maps of the locations and identities of hundreds of microbial species in complex communities. We show that 10-bit HiPR-FISH can distinguish between 1,023 isolates of Escherichia coli, each fluorescently labelled with a unique binary barcode. HiPR-FISH, in conjunction with custom algorithms for automated probe design and analysis of single-cell images, reveals the disruption of spatial networks in the mouse gut microbiome in response to treatment with antibiotics, and the longitudinal stability of spatial architectures in the human oral plaque microbiome. Combined with super-resolution imaging, HiPR-FISH shows the diverse strategies of ribosome organization that are exhibited by taxa in the human oral microbiome. HiPR-FISH provides a framework for analysing the spatial ecology of environmental microbial communities at single-cell resolution.

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.

Polymer Simulations of Heteromorphic Chromatin Predict the 3D Folding of Complex Genomic Loci.

Chromatin folded into 3D macromolecular structures is often analyzed by chromosome conformation capture (3C) and fluorescence in situ hybridization (FISH) techniques, but these frequently provide contradictory results. Chromatin can be modeled as a simple polymer composed of a connected chain of units. By embedding data for epigenetic marks (H3K27ac), chromatin accessibility (assay for transposase-accessible chromatin using sequencing [ATAC-seq]), and structural anchors (CCCTC-binding factor [CTCF]), we developed a highly predictive heteromorphic polymer (HiP-HoP) model, where the chromatin fiber varied along its length; combined with diffusing protein bridges and loop extrusion, this model predicted the 3D organization of genomic loci at a population and single-cell level. The model was validated at several gene loci, including the complex Pax6 gene, and was able to determine locus conformations across cell types with varying levels of transcriptional activity and explain different mechanisms of enhancer use. Minimal a priori knowledge of epigenetic marks is sufficient to recapitulate complex genomic loci in 3D and enable predictions of chromatin folding paths.

Lamins Organize the Global Three-Dimensional Genome from the Nuclear Periphery.

Lamins are structural components of the nuclear lamina (NL) that regulate genome organization and gene expression, but the mechanism remains unclear. Using Hi-C, we show that lamins maintain proper interactions among the topologically associated chromatin domains (TADs) but not their overall architecture. Combining Hi-C with fluorescence in situ hybridization (FISH) and analyses of lamina-associated domains (LADs), we reveal that lamin loss causes expansion or detachment of specific LADs in mouse ESCs. The detached LADs disrupt 3D interactions of both LADs and interior chromatin. 4C and epigenome analyses further demonstrate that lamins maintain the active and repressive chromatin domains among different TADs. By combining these studies with transcriptome analyses, we found a significant correlation between transcription changes and the interaction changes of active and inactive chromatin domains These findings provide a foundation to further study how the nuclear periphery impacts genome organization and transcription in development and NL-associated diseases.

High-throughput single telomere analysis using DNA microarray and fluorescent in situ hybridization.

The human telomere system is highly dynamic. Both short and long leucocyte average telomere lengths (aTL) are associated with an increased risk of cancer and early death, illustrating the complex relationship between TL and human health and the importance of assessing TL distributions with single TL analysis. A DNA microarray and telomere fluorescent in situ hybridization (DNA-array-FISH) approach was developed to measure the base-pair (bp) lengths of single telomeres. On average 32000 telomeres were measured per DNA sample with one microarray chip assaying 96 test DNA samples. Various telomere parameters, i.e. aTL and the frequency of short/long telomeres, were computed to delineate TL distribution. The intra-assay and inter-assay coefficient of variations of aTL ranged from 1.37% to 3.98%. The correlation coefficient (r) of aTL in repeated measurements ranged from 0.91 to 1.00, demonstrating high measurement precision. aTLs measured by DNA-array-FISH predicted aTLs measured by terminal restriction fragment (TRF) analysis with r ranging 0.87-0.99. A new accurate and high-throughput method has been developed to measure the bp lengths of single telomeres. The large number of single TL data provides an opportunity for an in-depth analysis of telomere dynamics and the complex relationship between telomere and age-related diseases.

Imaging and quantification of human and viral circular RNAs.

We present a robust approach for cellular detection, imaging, localization, and quantification of human and viral encoded circular RNAs (circRNA) using amplified fluorescence in situ hybridization (ampFISH). In this procedure, a pair of hairpin probes bind next to each other at contiguous stretches of sequence and then undergo a conformational reorganization which initiates a target-dependent hybridization chain reaction (HCR) resulting in deposition of an amplified fluorescent signal at the site. By harnessing the capabilities of both ampFISH and single-molecule FISH (smFISH), we selectively identified and imaged circular RNAs and their linear counterparts derived from the human genome, SARS-CoV-2 (an RNA virus), and human cytomegalovirus (HCMV, a DNA virus). Computational image processing facilitated accurate quantification of circular RNA molecules in individual cells. The specificity of ampFISH for circular RNA detection was confirmed through an in situ RNase R treatment that selectively degrades linear RNAs without impacting circular RNAs. The effectiveness of circular RNA detection was further validated by using ampFISH probes with mismatches and probe pairs that do not bind to the continuous sequence in their target RNAs but instead bind at segregated sites. An additional specificity test involved probes against the negative strands of the circular RNA sequence, absent in the cell. Importantly, our technique allows simultaneous detection of circular RNAs and their linear counterparts within the same cell with single molecule sensitivity, enabling explorations of circular RNA biogenesis, subcellular localization, and functions.

Spatially resolved whole transcriptome profiling in human and mouse tissue using Digital Spatial Profiling.

Emerging spatial profiling technology has enabled high-plex molecular profiling in biological tissues, preserving the spatial and morphological context of gene expression. Here, we describe expanding the chemistry for the Digital Spatial Profiling platform to quantify whole transcriptomes in human and mouse tissues using a wide range of spatial profiling strategies and sample types. We designed multiplexed in situ hybridization probes targeting the protein-coding genes of the human and mouse transcriptomes, referred to as the human or mouse Whole Transcriptome Atlas (WTA). Human and mouse WTAs were validated in cell lines for concordance with orthogonal gene expression profiling methods in regions ranging from ∼10-500 cells. By benchmarking against bulk RNA-seq and fluorescence in situ hybridization, we show robust transcript detection down to ∼100 transcripts per region. To assess the performance of WTA across tissue and sample types, we applied WTA to biological questions in cancer, molecular pathology, and developmental biology. Spatial profiling with WTA detected expected gene expression differences between tumor and tumor microenvironment, identified disease-specific gene expression heterogeneity in histological structures of the human kidney, and comprehensively mapped transcriptional programs in anatomical substructures of nine organs in the developing mouse embryo. Digital Spatial Profiling technology with the WTA assays provides a flexible method for spatial whole transcriptome profiling applicable to diverse tissue types and biological contexts.

Improved enzymatic labeling of fluorescent in situ hybridization probes applied to the visualization of retained introns in cells.

Fluorescence in situ hybridization (FISH) is a widely used tool for quantifying gene expression and determining the location of RNA molecules in cells. We present an improved method for FISH probe production that yields high-purity probes with a wide range of fluorophores using standard laboratory equipment at low cost. The method modifies an earlier protocol that uses terminal deoxynucleotidyl transferase to add fluorescently labeled nucleotides to synthetic deoxyoligonucleotides. In our protocol, amino-11-ddUTP is joined to an oligonucleotide pool prior to its conjugation to a fluorescent dye, thereby generating pools of probes ready for a variety of modifications. This order of reaction steps allows for high labeling efficiencies regardless of the GC content or terminal base of the oligonucleotides. The degree of labeling (DOL) for spectrally distinct fluorophores (Quasar, ATTO, and Alexa dyes) was mostly >90%, comparable with commercial probes. The ease and low cost of production allowed the generation of probe sets targeting a wide variety of RNA molecules. Using these probes, FISH assays in C2C12 cells showed the expected subcellular localization of mRNAs and pre-mRNAs for Polr2a (RNA polymerase II subunit 2a) and Gapdh, and of the long noncoding RNAs Malat1 and Neat1 Developing FISH probe sets for several transcripts containing retained introns, we found that retained introns in the Gabbr1 and Noc2l transcripts are present in subnuclear foci separate from their sites of synthesis and partially coincident with nuclear speckles. This labeling protocol should have many applications in RNA biology.

WISH-BONE: Whole-mount in situ histology, to label osteocyte mRNA and protein in 3D adult mouse bones.

Bone is a three-dimensional (3D) highly dynamic tissue under constant remodeling. Commonly used tools to investigate bone biology require sample digestion for biomolecule extraction or provide only two-dimensional (2D) spatial information. There is a need for 3D tools to investigate spatially preserved biomarker expression in osteocytes. In this work, we present a new method, WISH-BONE, to label osteocyte messenger RNA (mRNA) and protein in whole-mount mouse bone. For mRNA labeling, we used hybridization chain reaction-fluorescence in situ hybridization (HCR-FISH) to label genes of interest in osteocytes. For protein labeling, samples were preserved using an epoxy-based solution that protects tissue structure and biomolecular components. Then an enzymatic matrix permeabilization step was performed to enable antibody penetration. Immunostaining was used to label various proteins involved in bone homeostasis. We also demonstrate the use of customized fluorescent nanobodies to target and label proteins in the cortical bone (CB). However, the relatively dim signal observed from nanobodies' staining limited detection. mRNA and protein labeling were performed in separate samples. In this study, we share protocols, highlight opportunities, and identify the challenges of this novel 3D labeling method. They are the first protocols for whole-mount osteocyte 3D labeling of mRNA and protein in mature mouse bones. WISH-BONE will allow the investigation of molecular signaling in bone cells in their 3D environment and could be applied to various bone-related fields of research.

Transcriptome-scale super-resolved imaging in tissues by RNA seqFISH.

Imaging the transcriptome in situ with high accuracy has been a major challenge in single-cell biology, which is particularly hindered by the limits of optical resolution and the density of transcripts in single cells1-5. Here we demonstrate an evolution of sequential fluorescence in situ hybridization (seqFISH+). We show that seqFISH+ can image mRNAs for 10,000 genes in single cells-with high accuracy and sub-diffraction-limit resolution-in the cortex, subventricular zone and olfactory bulb of mouse brain, using a standard confocal microscope. The transcriptome-level profiling of seqFISH+ allows unbiased identification of cell classes and their spatial organization in tissues. In addition, seqFISH+ reveals subcellular mRNA localization patterns in cells and ligand-receptor pairs across neighbouring cells. This technology demonstrates the ability to generate spatial cell atlases and to perform discovery-driven studies of biological processes in situ.

In situ visualization of m6A sites in cellular mRNAs.

N 6-methyladenosine (m6A) is an abundant RNA modification which plays critical roles in RNA function and cellular physiology. However, our understanding of how m6A is spatially regulated remains limited due to a lack of methods for visualizing methylated transcripts of interest in cells. Here, we develop DART-FISH, a method for in situ visualization of specific m6A sites in target RNAs which enables simultaneous detection of both m6A-modified and unmodified transcript copies. We demonstrate the ability of DART-FISH to visualize m6A in a variety of mRNAs across diverse cell types and to provide information on the location and stoichiometry of m6A sites at single-cell resolution. Finally, we use DART-FISH to reveal that m6A is not sufficient for mRNA localization to stress granules during oxidative stress. This technique provides a powerful tool for examining m6A-modified transcript dynamics and investigating methylated RNA localization in individual cells.