Mol. Cells 2020; 43(7): 591~599  https://doi.org/10.14348/molcells.2020.0020
Mapping Cellular Coordinates through Advances in Spatial Transcriptomics Technology
Joji Marie Teves 1,2 and Kyoung Jae Won 1,2,*
1Biotech Research and Innovation Centre (BRIC), University of Copenhagen, DK-2200 Copenhagen, Denmark, 2Novo Nordisk Foundation Center for Stem Cell Biology (DanStem), Faculty of Health Sciences, Faculty of Health and Medical Sciences, University of Copenhagen, DK-2200 Copenhagen, Denmark
Received January 14, 2020; Revised April 22, 2020; Accepted May 10, 2020.; Published online June 4, 2020.
© Korean Society for Molecular and Cellular Biology. All rights reserved.

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported License. To view a copy of this license, visit (http://creativecommons.org/licenses/by-nc-sa/3.0/).
ABSTRACT
Complex cell-to-cell communication underlies the basic processes essential for homeostasis in the given tissue architecture. Obtaining quantitative gene-expression of cells in their native context has significantly advanced through single-cell RNA sequencing technologies along with mechanical and enzymatic tissue manipulation. This approach, however, is largely reliant on the physical dissociation of individual cells from the tissue, thus, resulting in a library with unaccounted positional information. To overcome this, positional information can be obtained by integrating imaging and positional barcoding. Collectively, spatial transcriptomics strategies provide tissue architecture-dependent as well as position-dependent cellular functions. This review discusses the current technologies for spatial transcriptomics ranging from the methods combining mechanical dissociation and single-cell RNA sequencing to computational spatial re-mapping.
Keywords: cellular communication, single-cell RNA, spatial transcriptomics, tissue architecture
INTRODUCTION

Cell-to-cell communication is essential to maintain proper tissue homeostasis. Disruption of homeostatic cellular communication underlies many pathologic cellular transformations including cancer (Oktay et al., 2015). Studying the complexity of healthy tissue architecture and abnormal transformations both at the cellular and transcriptional level is important in improving the understanding of key pathways that can be targeted for therapeutic strategies. Recently, single cell RNA-sequencing (scRNAseq) technologies have revolutionized our understanding of gene expression by quantifying the transcriptome of individual cells. Moreover, the development of computational approaches to quantify large transcriptomic data alongside scRNAseq technology has provided transcriptomic information for previously uncharacterized cell types and has made it possible to study their dynamics at a population level (Grun et al., 2015; Patel et al., 2014; Pollen et al., 2014; Proserpio et al., 2016; Shalek et al., 2014; Trapnell et al., 2014).

However, scRNAseq technologies have the intrinsic limitation of losing positional information during tissue dissociation into single cells. Positional information is a critical aspect when studying tissue architecture to account for how physically interacting cells and signal exchanges maintain or alter homeostasis. Indeed, several strategies have been employed to provide positional information relative to transcriptomic data. In particular, tissue microdissection followed by RNA sequencing provides approximated positional information based on microdissected fragments (Combs and Eisen, 2013). Moreover, recent advances in a set of techniques collectively called spatial transcriptomics allow positional information to be identified at a single cell resolution.

By showing both expression and transcriptome at a single cell resolution, spatial transcriptomics can provide information such as tissue architecture and cell interactions. However, obtaining both position and expression information at a single cell resolution is still technically challenging. Historically, positional information has been obtained by physically taking cells positioned in a tissue, followed by expression profiling. Later, cells expressing specific genes were captured in the image to provide the position of cells as well as the gene expression levels (Chen et al., 2018). Also, a variety of barcode-based approaches to point the position of cells in the two-dimensional tissue section have also been introduced (Chen et al., 2018). Currently 10X Genomics (Visium; https://www.10xgenomics.com/products/spatial-gene-expression/) produces kit for spatial transcriptomics based on barcode-based approaches, and CARTANA (http://cartana.se) is based on in situ sequencing. Indeed, spatial transcriptomics could characterize information such as position-dependent gene expression profiles and cell composition, spatio-temporal gene regulation and cellular communication. In this review, we discuss various spatial transcriptomics methodologies and their applications. Table 1 summarizes the information such as cell/tissue type, experimental method and detection sensitivity for the for each technology we discuss.

TISSUE DISSECTION FOLLOWED BY TRANSCRIPTOMIC PROFILING

One of the earliest methods to study spatial transcriptomics was through the use of tissue microdissection followed by RNA sequencing. In Drosophila embryos, this approach has allowed to uncover spatial-specific gene expression patterns using coordinate-bound cryosectioning (Fig. 1A) (Combs and Eisen, 2013). Similarly, using 50 to 100 cryosectioned thin slices from zebrafish embryos, the Tomo-Seq method has provided three-dimensional (3D) spatial expression patterns with the aid of computational reconstruction of the zebrafish tissue architecture (Junker et al., 2014). The microdissected slices were further used to the reconstruction of murine brain, providing a 3D image of gene expression (Okamura-Oho et al., 2012).

Advancements such as laser capture microdissection (LCM) enabled a precise capture of targeted cells, or even single cells, while retaining intact tissue structure (Fig. 1B) (Datta et al., 2015). Subsequently, combining LCM and RNA sequencing was used to resolve spatially bound transcriptomic profiles of rare cell population (Nichterwitz et al., 2016). Comparably, geographical position sequencing (Geo-seq) is technique combining LCM with scRNAseq (Chen et al., 2017; Xue et al., 2019). Moreover, LCM has been used in various applications to provide position-based transcriptional information. For instance, LCM followed by RNAseq in mouse intestinal epithelium revealed the transcriptome of spatially zoned areas along the villus axis, which leads to spatial reconstruction of the tissue from scRNAseq data (Moor et al., 2018). LCM enables accurate separation of a small number of cells while preserving the tissue morphology. However, LCM procedures are labor-intensive and expensive to perform (Chung and Shen, 2015).

FLUORESCENCE IN SITU HYBRIDIZATION (FISH) APPROACHES FOR SPATIAL TRANSCRIPTOMICS

Multiplexed image-based transcriptomics is an emerging technology for spatial detection of RNAs. Particularly, in situ hybridization followed by microscopic analysis of intact tissues has the potential to provide direct information regarding the spatial organization of the cell’s transcriptome (Levsky and Singer, 2003). This approach utilizes fluorescent-labelled nucleotide probes to detect specific mRNAs localized at different spatial coordinates in a cell or tissue. Single molecule FISH (smFISH) is a pioneering method that uses many short oligonucleotide probes to target mRNAs (Raj et al., 2008). smFISH provides the cell-to-cell variation in transcript abundance and even the subcellular localization of a given RNA (Raj et al., 2010). However, the number of RNA species that can be simultaneously measured by smFISH is limited. To increase the number of detected RNA species, combinatorial fluorescence has been developed to visualize 10 or more genes simultaneously (Jakt et al., 2013). smFISH can obtain spatial information of targeted mRNA with high sensitivity and low false positives or false negatives. However, smFISH is limited by technical factors regarding probe design and binding such as relatively low signal intensity and tendencies towards nonspecific probe binding. This limitation was improved through the development of click-amplifying FISH (clampFISH), which utilizes bio-orthogonal click-chemistry to fix a padlock-style probe/target interaction that subsequently enhances the specificity and fluorescence of signal amplification (Rouhanifard et al., 2018).

Another FISH technology development includes a cyclic approach to single-molecule fluorescence in situ hybridization methodology (osmFISH) wherein multiple smFISH rounds are repeatedly applied to increase the number of detected RNA species (Fig. 1C) (Codeluppi et al., 2018). In osmFISH, the number of targets becomes the number of fluorescence channels multiplied by the number of hybridization cycles thus, significantly increasing the number of transcripts quantified in each round. In a recent study, osmFISM profiled 33 genes over 13 imaging rounds in mouse somatosensory cortex (Codeluppi et al., 2018). Similarly, sequential FISH (seqFISH) utilizes sequential labelling of mRNAs using a set of FISH probes designed with a single type of fluorophore, generating a barcode of color labels (Fig. 1C) (Lubeck et al., 2014). For this, mRNA hybridization is performed at each round using the same FISH probes but labeled with a different dye. Therefore, seqFISH can scale the number of targets exponentially for each round. SeqFISH has been applied to detect thousands of transcripts per cells in mouse hippocampus using multiple colors (Shah et al., 2016). However, seqFISH is hampered by the optical density of each targeted transcript, which limits the number of detectable RNAs despite the multiplexing capacity. In an improved version, SeqFISH+ made it possible to probe 10,000 genes by solving the optical density problem using ‘pseudocolours’ to dilute cellular RNAs and enabling spatial reconstructions at high resolution (Eng et al., 2019).

Alternatively, the multiplexed error-robust FISH (MERFISH) technology, instead of barcoding with different color combinations, determines the presence or absence of fluorescence using a two-stop hybridization protocol (Fig. 1C) (Moffitt and Zhuang, 2016). Each RNA specie is encoded with a N-bit binary word from N rounds of hybridization that shows the presence or absence of a color. MERFISH, as well as other multiplexed smFISH-based RNA profiling methods, requires non-overlapping signals from individual RNAs. Improvements in microscopic detection combined with MERFISH has been used to overcome technical drawbacks such as signal detection and background correction (Wang et al., 2018a). Spatial information resolution was increased when MERFISH was combined with high-content imaging of sub-cellular structures to accurately determine the compartmentalization of RNAs (Xia et al., 2019).

HYBRID COMBINATIONS OF IN SITU AND SEQUENCING-BASED TRANSCRIPTOMIC PROFILING

To maximize the features of in situ profiling, while also obtaining multi-omic spatially resolved information from the same biological samples, combined adaptations of in situ profiling and biochemical modifications of sequencing library preparations have been developed. One alternative to single cell extraction from tissues or multiplexed FISH is sequencing of transcripts directly on tissue sections using nucleotide-fluorescent sequencing or enzymatic amplification (Ke et al., 2013; Lee et al., 2014). In particular, the fluorescent in-situ sequencing (FISSEQ) technology utilizes in situ conversion of RNA into cross-linked cDNA amplicons followed by manual sequencing by fluorescent tag detection on a confocal microscope. This results in the detection of thousands of spatially resolved nucleotide sequences of genes from high content-images (Lee et al., 2015). Additionally, RNA sequencing of intact tissue using STARmap (spatially-resolved transcript amplicon readout mapping) provides 3D sequence information by labeling cellular RNA with pairs of DNA probes followed by enzymatic amplification. The resulting amplification produces a DNA nanoball that eliminates background signal caused by mislabeling of single probes (Wang et al., 2018b). CARTANA is a company that employed in situ sequencing for customizable gene panels (Hernandez et al., 2019). CARTANA can profile more than 100 genes in a high throughput pipeline (10 cm2 [1 million cells] per week per microscope).

Barcode based approaches for spatial transcriptomics

Recent strategies using barcodes on native tissue on a slide, namely Spatial transcriptomics and Slide-seq, have been developed for high-resolution spatial resolution (Rodriques et al., 2019; Stahl et al., 2016). This approach dissects a histological section with a grid, where each spot is labelled with barcoded oligonucleotide primers to capture adjacent tissue mRNA. The resulting transcripts are then reverse transcribed to cDNA and positional information is contained within their respective barcodes (Fig. 1D). The Visium Spatial Gene Expression Solution (10× genomics) is based on the barcode-based approach (https://www.10xgenomics.com/products/spatial-gene-expression/). The current resolution of commercial Spatial Transcriptomics is limited to 100 µm, capturing an average of 3 to 30 cells per regular grid. High-definition spatial transcriptomics (HDST) uses 2-μm beads to enhance the resolution (Vickovic et al., 2019) while the Slide-seq technology uses 10-μm beads containing individual position barcode (Rodriques et al., 2019).

Sequencing direct cell contact

Although FISH and barcode-based approaches provide cellular positioning within a tissue, approaches to use cell contact information (thus provide relative positional information) have also been developed. This approach utilizes trasncriptomic information obtained from physically-interacting cells. ProximID is an approach where cells are gently dissociated to retain cells that are in physical contact, and then subsequently processed for RNAseq to obtain transcriptomic information of doublets or triplets of interacting cells (Fig. 1E) (Boisset et al., 2018). In parallel, RNAs from single cells are also sequenced (scRNAseq), which provide the reference information of the doublets and triplets. Without microdissection, ProximID identified that Tac1+ enteroendocrine cell–Lgr5+ stem cells interact in the intestinal crypt (Boisset et al., 2018), thus providing insights to previously uncharacterized interactions that may play a role in stem cell driven regeneration. ProximID collects the interacting cells that survived mild dissociation therefor this method can potentially be biased towards collecting strongly interacting cells. Similarly, paired cell RNAseq (pcRNAseq) sequenced RNA from cellular doublets composed of liver endothelial cells (LECs) and hepatocytes generated from extensive flow cytometric sorting (Halpern et al., 2018). pcRNAseq and parallel scRNAseq allowed the identification of zone-dependent LEC gene expression relative to the hepatocyte gradient (Halpern et al., 2018). These approaches utilize the innate interactions between cells in the tissue. Similarly, PIC-seq sequences RNAs from physically interacting cells and identified the composition of the interacting cells and investigated genes that are differentially expressed when interacting with other cell types (Giladi et al., 2020).

Computational reconstruction of spatial information from scRNAseq data

Tissue architecture reconstruction is an algorithm-based strategy used for reverse identification of a cell’s spatial origin based on gene expression data and a positional reference (Fig. 1F). Tomo-Seq used slices of cryosectioned zebrafish embryo to reconstruct 3D expression patterns (Junker et al., 2014). Similarly, scRNAseq information was used to reconstruct zoned areas along the villus axis of the intestinal epithelium with respect to the location information provided by LCM (Moor et al., 2018). Computational algorithms have been developed to predict the positional origin of cells using a model trained using the in situ hybridization (ISH)-based gene expression data as reference (Achim et al., 2015; Satija et al., 2015). Lobule zonation markers identified using smFISH in the liver has been used to discover zonation-specific genes (Halpern et al., 2017). Additionally, Distmap enabled the reconstruction of Drosophila embryo using scRNAseq and gene expression of 84 marker genes and performed virtual in situ hybridization (Karaiskos et al., 2017). Similar to pseudo-time analysis for scRNAseq data, pseudo-spatial trajectory analysis using scRNAseq from the dissected colon samples identified location-dependent gene expression that controls epithelial-to-mesenchymal transition (EMT) (McFaline-Figueroa et al., 2019). Recently, architecture reconstruction has been achieved without a spatial reference in novoSpaRc under that assumption that gene expression between nearby cells is generally more similar than gene expression between cells separated by a larger distance (Nitzan et al., 2019). Computational strategies in spatial transcriptomics complement the developments in biochemical and mechanical techniques.

Cellular coordinates as a critical component for studying cell-cell interactions

Cellular location is important to understand how cells interact with each other within a tissue. Current technologies in spatial transcriptomics such the ones mentioned in this review have provided vital information on the cellular identity, transcriptional changes, critical gene markers for cell communication and even previously uncharacterized cellular interactions relative to its position in the tissue. Spatial reconstruction allowed the detection of a comprehensive set of liver genes differentially expressed across the lobule (Halpern et al., 2017). The LECs, the interacting partner with hepatocytes, also exhibited differential gene expression across the lobule (Halpern et al., 2018). Moreover, proximity-based reconstruction of hematopoietic stem/progenitor cell’s (HPSC) microanatomy in bone marrow followed by scRNAseq post-transplantation allowed the retrieval of osteolineage cells (OLCs) that were proximal or distal to HPSCs (Silberstein et al., 2016). Through proximity-based profiling, the study identified differentially expressed genes in OLCs affected by nearby HSPC and niche factors regulating HSPC quiescence in vivo (Silberstein et al., 2016). Likewise, proximity-based “clumplet” cell sequencing technique ProximID identified the physical interaction of hormone-producing Tac-enteroendocrine cells with Lgr5 stem cell in the intestinal crypt (Boisset et al., 2018). These results suggest a previously uncharacterized hormone-related secretory pathway involved in homeostatic maintenance of intestinal crypts. Furthermore, spatial transcriptomics could also spatially resolve interactions between foreign cells as demonstrated through scDual-Seq by quantifying host-pathogen interactions (Avital et al., 2017). Indeed, accounting for a cell’s positional information allows an extra component to help resolve the complex interplay of key regulatory pathways and cellular interactions involved in homeostatic tissue maintenance. Advancements in spatial transcriptomics technologies ranging from mechanical cell capture, in situ biochemistry, enhanced microscopy and/or combinations thereof, as well as sequencing innovations along with advanced computational algorithms will further shed light on complex tissue architecture in the context of tissue homeostasis and pathologies.

BEYOND CELLULAR COORDINATES: SPATIAL TRANSCRIPTOMICS ACROSS BIOLOGY

Spatial transcriptomic methods can serve as a powerful tool in capturing cell-to-cell properties at a tissue-level resolution. Using spatial transcriptomics enabled the quantitative elucidation of important cellular dynamics. Besides the cell and the tissues that were used to develop the methodology in Table 1, the number of applications for spatial transcriptomics are increasing. For instance, spatial quantification of neural and mesodermal lineage trajectories from Tomo-Seq was associated with the anterior-posterior axis of gastruloids developmental (van den Brink et al., 2020). Barcode based spatial transcriptomics have been applied to detect spatial crosstalk in Alzheimer’s disease wherein novel cellular interactions were uncovered between microglia and astroglial cells in amyloid-β plaque niches (Chen et al., 2019). Also, these approaches have been applied to study spatial localization of cell population in pancreatic ductal adenocarcinoma (Moncada et al., 2018) and HER2+ breast tumors (Salmén et al., 2018). In situ sequencing has been applied to study spatio-temporal cell fate decisions in murine neural crest (Soldatov et al., 2019), the immunological architecture of the murine tuberculosis granulomas (Carow et al., 2019), mouse brain (Tiklova et al., 2019) and human heart development (Asp et al., 2019) as well as primary colorectal cancers (Mathot et al., 2017). PIC-seq detected preferential interaction between two immune cells (T-cell and dendritic cells) which revealed interaction-specific phenomena important during pathogen stimulation (Giladi et al., 2020). Ligand-partner expression patterns were used to extract cross-interactions between breast cancer and cancer-associated fibroblast (CAF) by ICELLNet (Noël et al., 2020). A similar approach was implemented in NicheNet describing the interactions between tumor and immune cell microenvironment (Browaeys et al., 2020).

DISCUSSION

Studies on cellular identity combined with their native location has been significantly advanced by scRNAseq and spatial transcriptomic technologies. This review highlights the current technological developments that have refined spatial transcriptomics. Pioneering image-based approaches use FISH to visualize cellular localization. However, overlaps between fluorescent signals have been a major limitation in image resolution (Burgess, 2019). This was surmounted by barcoding-based improvements that lead to the development of 2-μm beads for HDST to measure single cell transcriptomic profiles (Vickovic et al., 2019). Although combining scRNAseq and barcoding-based approaches can determine transcriptome-wide profiles, they are still currently limited by their lack of sensitivity towards RNA detection (Kolodziejczyk et al., 2015) thus providing opportunities for further technological refinements.

As cells constantly communicate with other cells within a tissue, the study of cell communication will become more popular using scRNAseq and spatial transcriptomics. Previous approaches to study cell communication relied on the co-expression of ligand-receptor pairs between cell types (Kumar et al., 2018; Skelly et al., 2018). Neighboring cells communicate via complex paracrine signaling networks (Roy and Kornberg, 2015). Studying paracrine signaling using single cell resolution spatial transcriptomics provides an opportunity to enhance our understanding of cell communication.

Subsequently, we can then begin to question spatio-temporal changes of cells during development or disease progression. The majority of current spatial transcriptomics required a section of a tissue and cannot reveal spatio-temporal changes effectively. Attempts such as 4D-seq (unpblished; https://techtransfer.universityofcalifornia.edu/NCD/30312.html) uses two-photon microscopy and DNA-labelling to capture cellular spatio-temporal variation within a tissue prior to scRNA-seq analysis. As technical resolution continually improves, it is possible to foresee studies integrating high-content 4D image data with scRNAseq in order to gain a deeper understanding of spatio-temporal cellular regulation.

ACKNOWLEDGMENTS

We are grateful to Sung-Whan Bae for graphic illustration and Patrick Martin for proofreading of the manuscript. The Novo Nordisk Foundation Center for Stem Cell Biology is supported by a Novo Nordisk Foundation grant (NNF17CC0027852). This work is also supported by Lundbeck Foundation Experiment grant (R324-2019-1649) to K.J.W. and Marie Sklodowska Curie COFUND Fellowship to J.M.T.

AUTHOR CONTRIBUTIONS

Both J.M.T. and K.J.W. wrote the manuscript.

CONFLICT OF INTEREST

The authors have no potential conflicts of interest to disclose.

FIGURES
Fig. 1. Diverse approaches to associate spatial information with transcriptomics. (A) Cryosection provides positional information. (B) LCM provides fine resolution (even to single cell) positional information. (C) Image-based single cell level spatial transcriptomic approaches. osmFISH labels RNA with a number of colors each time for different genes. seqFish uses a combination of colors to mark RNAs. MERFISH labels presence or absence of fluorescence. (D) Spatial transcriptomics uses barcodes to spatially distinguish each spot. (E) RNAseq for interacting cells provides relative spatial information. (F) Spatial reconstruction uses transcriptomic information to reconstruct original spatial information.
TABLES

Table 1

Summary of selected technologies for spatial profiling of cells

Approach Technology Input material Experimental method Quantification method Representative detection sensitivity Detection rangea Reference
LCM-seq LCM-seq Primary mouse brain and spinal cord tissue;Post-mortem human brain Laser capture microdissection NGS data analysis (DESeq2 + GO analysis) ~1,743 to 14,893 genes per 0.1 RPKM SG to HT (Combs and Eisen, 2013)
FISH-based smFISH A549 and CHO cell line;Primary rat hippocampus neurons; Yeast,Caenorhabditis elegans, Drosophila melanogaster Fluorescence imaging + photo-bleaching on fixed cells Probe-based computational identification of mRNA targets ~3 mRNA species per cell SG (Halpern et al., 2017)
MERFISH IMR90 cell line;U2OS cell line Multiplexed fluorescence imaging of target probes on fixed cells (+ clearing) Probe-based encoding + GO analysis ~100 to 1,000 RNA species per cell SG to MT (Salmén et al., 2018)(Lee et al., 2015)
seqFISH+ NIH3T3 cell line;Primary mouse brain tissue Sequential fluorescence of pseudocolor probes Probe-based encoding + scRNA-seq-based spatial localization mapping ~10,000 genes per cell SG to HT (Kumar et al., 2018)
In situ sequencing (ISS) barcode-based FISSEQ HeLa, 293A, COS1, U2OS, iPSC, primary fibroblasts and bipolar neurons cell lines;iPS-derived 3D organoids;Primary mouse embryo and brain tissue;Drosophila embryos Reverse transcript probes + sequence-by-ligation Probe-based calling + 3D image deconvolution ~200 to 400 mRNA per cell; scalable 5X MT (McFaline-Figueroa et al., 2019)
STARmap Primary mouse cortical neuron cells;Primary mouse brain tissue Hydrogel-based isolation of target probes + SEDAL sequencing Probe-based calling + 2D/3D cell segmentation + differential gene expression analysis ~160 to 1,020 genes simultaneously; scalable to ~30,000 cells MT to HT (Moffitt and Zhuang, 2016)
Spatial and single-cell sequencing-based Spatial reconstruction from single-cell transcriptomics (Seurat) Danio rerio embryo tissue Tissue dissociation + strand-specific, scRNA-seq modified from SMART protocol NGS analysis + spatial location inference Spatial reconstruction from 851 single-cell reference HT (Patel et al., 2014)
Spatial transcriptomics Primary mouse olfactory bulbs and brain tissue;Breast cancer biopsy tissue Spatial oligodT barcode array + cDNA synthesis + RNA-seq Transcriptome analysis On surface: 9.6 M unique transcripts per 400 M readsIn solution: 18 M unique transcripts per 290 reads HT (Moncada et al., 2018)
Spatially-interacting cells pcRNAseq (to detect liver zonation) Primary mouse liver tissue; Primary mouse hepatocyte cells Mild dissociation + single-cell and paired cell sorting + MARS-Seq ScRNAseq analysis + landmark gene identification + spatial zonation inference ~70 landmark genes MT to HT (Pollen et al., 2014)
ProximID (to detect interacting cell network) Primary mouse bone marrow and fetal liver tissue; Primary mouse intestinal crypts Mild dissociation + single-cell sorting or pipette dissociation + CEL-Seq ScRNAseq analysis + cell interaction inference + network analysis Example: ~17 to 78 simulated preferential cell interactions SG to HT (Nitzan et al., 2019)
PIC-Seq (to detect physical interactions) Primary mouse spleen tissue Mild dissociation + population-based low cytometric sorting + MARS-seq single-cell RNA sequencing NGS analysis + computational calling of physically interacting cells + differential analysis Example: ~348 differential genes from ~2,389 simulated physically interacting cells HT (Okamura-Oho et al., 2012)

A549, adenocarcinomic human alveolar basal epithelial cells; CHO, epithelial cell line derived from the ovary of the Chinese hamster; IMR90, human foetal lung cells; U2OS, Human Bone Osteosarcoma Epithelial Cells; NIH3T3, mouse embryonic fibroblast cells; HeLa, human cervical cancer cells; COS1, African green monkey kidney fibroblast-like cell; 293A, human embryonic kidney cells; iPSC, induced pluripotent stem cells.

a Detection range: SG, single gene; MT, medium throughput (targeted transcript capture sequencing); HT, high throughput (non-targetted transcriptome-wide sequencing).


REFERENCES
  1. Achim K., Pettit J.B., Saraiva L.R., Gavriouchkina D., Larsson T., Arendt D., and Marioni J.C. (2015). High-throughput spatial mapping of single-cell RNA-seq data to tissue of origin. Nat. Biotechnol. 33, 503-509.
    Pubmed CrossRef
  2. Asp M., Giacomello S., Larsson L., Wu C., Furth D., Qian X., Wardell E., Custodio J., Reimegard J., and Salmen F., et al. (2019). A spatiotemporal organ-wide gene expression and cell atlas of the developing human heart. Cell 179, 1647-1660.e19.
    Pubmed CrossRef
  3. Avital G., Avraham R., Fan A., Hashimshony T., Hung D.T., and Yanai I. (2017). scDual-Seq: mapping the gene regulatory program of Salmonella infection by host and pathogen single-cell RNA-sequencing. Genome Biol. 18, 200.
    Pubmed KoreaMed CrossRef
  4. Boisset J.C., Vivie J., Grun D., Muraro M.J., Lyubimova A., and van Oudenaarden A. (2018). Mapping the physical network of cellular interactions. Nat. Methods 15, 547-553.
    Pubmed CrossRef
  5. Browaeys R., Saelens W., and Saeys Y. (2020). NicheNet: modeling intercellular communication by linking ligands to target genes. Nat. Methods 17, 159-162.
    Pubmed CrossRef
  6. Burgess D.J. (2019). Spatial transcriptomics coming of age. Nat. Rev. Genet. 20, 317.
    Pubmed CrossRef
  7. Carow B., Hauling T., Qian X., Kramnik I., Nilsson M., and Rottenberg M.E. (2019). Spatial and temporal localization of immune transcripts defines hallmarks and diversity in the tuberculosis granuloma. Nat. Commun. 10, 1823.
    Pubmed KoreaMed CrossRef
  8. Chen J., Suo S., Tam P.P., Han J.J., Peng G., and Jing N. (2017). Spatial transcriptomic analysis of cryosectioned tissue samples with Geo-seq. Nat. Protoc. 12, 566-580.
    Pubmed CrossRef
  9. Chen W.T., Lu A., Craessaerts K., Pavie B., Sala Frigerio C., Mancuso R., Qian X., Lalakova J., Kühnemund M., and Voytyuk I., et al. (2019). Spatial and temporal transcriptomics reveal microglia-astroglia crosstalk in the amyloid-β plaque cell niche of Alzheimer's disease. BioRxiv . https://doi.org/10.1101/719930.
    CrossRef
  10. Chen X., Teichmann S.A., and Meyer K.B. (2018). From tissues to cell types and back: single-cell gene expression analysis of tissue architecture. Annu. Rev. Biomed. Data Sci. 1, 29-51.
    CrossRef
  11. Chung S.H. and Shen W. (2015). Laser capture microdissection: from its principle to applications in research on neurodegeneration. Neural Regen. Res. 10, 897-898.
    Pubmed KoreaMed CrossRef
  12. Codeluppi S., Borm L.E., Zeisel A., La Manno G., van Lunteren J.A., Svensson C.I., and Linnarsson S. (2018). Spatial organization of the somatosensory cortex revealed by osmFISH. Nat. Methods 15, 932-935.
    Pubmed CrossRef
  13. Combs P.A. and Eisen M.B. (2013). Sequencing mRNA from cryo-sliced Drosophila embryos to determine genome-wide spatial patterns of gene expression. PLoS One 8, e71820.
    Pubmed KoreaMed CrossRef
  14. Datta S., Malhotra L., Dickerson R., Chaffee S., Sen C.K., and Roy S. (2015). Laser capture microdissection: big data from small samples. Histol. Histopathol. 30, 1255-1269.
    Pubmed KoreaMed CrossRef
  15. Eng C.L., Lawson M., Zhu Q., Dries R., Koulena N., Takei Y., Yun J., Cronin C., Karp C., and Yuan G.C., et al. (2019). Transcriptome-scale super-resolved imaging in tissues by RNA seqFISH. Nature 568, 235-239.
    Pubmed KoreaMed CrossRef
  16. Giladi A., Cohen M., Medaglia C., Baran Y., Li B., Zada M., Bost P., Blecher-Gonen R., Salame T.M., and Mayer J.U., et al. (2020). Dissecting cellular crosstalk by sequencing physically interacting cells. Nat. Biotechnol. 38, 629-637.
    Pubmed CrossRef
  17. Grun D., Lyubimova A., Kester L., Wiebrands K., Basak O., Sasaki N., Clevers H., and van Oudenaarden A. (2015). Single-cell messenger RNA sequencing reveals rare intestinal cell types. Nature 525, 251-255.
    Pubmed CrossRef
  18. Halpern K.B., Shenhav R., Massalha H., Toth B., Egozi A., Massasa E.E., Medgalia C., David E., Giladi A., and Moor A.E., et al. (2018). Paired-cell sequencing enables spatial gene expression mapping of liver endothelial cells. Nat. Biotechnol. 36, 962-970.
    Pubmed KoreaMed CrossRef
  19. Halpern K.B., Shenhav R., Matcovitch-Natan O., Toth B., Lemze D., Golan M., Massasa E.E., Baydatch S., Landen S., and Moor A.E., et al. (2017). Single-cell spatial reconstruction reveals global division of labour in the mammalian liver. Nature 542, 352-356.
    Pubmed KoreaMed CrossRef
  20. Hernandez I., Qian X.Y., Lalakova J., Verheyen T., Hilscher M., and Kuhnemund M. (2019). Mapping brain cell types with CARTANA in situ sequencing on the Nikon Ti2-E microscope. Nat. Methods 16 .
  21. Jakt L.M., Moriwaki S., and Nishikawa S. (2013). A continuum of transcriptional identities visualized by combinatorial fluorescent in situ hybridization. Development 140, 216-225.
    Pubmed CrossRef
  22. Junker J.P., Noel E.S., Guryev V., Peterson K.A., Shah G., Huisken J., McMahon A.P., Berezikov E., Bakkers J., and van Oudenaarden A. (2014). Genome-wide RNA tomography in the zebrafish embryo. Cell 159, 662-675.
    Pubmed CrossRef
  23. Karaiskos N., Wahle P., Alles J., Boltengagen A., Ayoub S., Kipar C., Kocks C., Rajewsky N., and Zinzen R.P. (2017). The Drosophila embryo at single-cell transcriptome resolution. Science 358, 194-199.
    Pubmed CrossRef
  24. Ke R., Mignardi M., Pacureanu A., Svedlund J., Botling J., Wahlby C., and Nilsson M. (2013). In situ sequencing for RNA analysis in preserved tissue and cells. Nat. Methods 10, 857-860.
    Pubmed CrossRef
  25. Kolodziejczyk A.A., Kim J.K., Svensson V., Marioni J.C., and Teichmann S.A. (2015). The technology and biology of single-cell RNA sequencing. Mol. Cell 58, 610-620.
    Pubmed CrossRef
  26. Kumar M.P., Du J., Lagoudas G., Jiao Y., Sawyer A., Drummond D.C., Lauffenburger D.A., and Raue A. (2018). Analysis of single-cell RNA-seq identifies cell-cell communication associated with tumor characteristics. Cell Rep. 25, 1458-1468.e4.
    Pubmed KoreaMed CrossRef
  27. Lee J.H., Daugharthy E.R., Scheiman J., Kalhor R., Ferrante T.C., Terry R., Turczyk B.M., Yang J.L., Lee H.S., and Aach J., et al. (2015). Fluorescent in situ sequencing (FISSEQ) of RNA for gene expression profiling in intact cells and tissues. Nat. Protoc. 10, 442-458.
    Pubmed KoreaMed CrossRef
  28. Lee J.H., Daugharthy E.R., Scheiman J., Kalhor R., Yang J.L., Ferrante T.C., Terry R., Jeanty S.S., Li C., and Amamoto R., et al. (2014). Highly multiplexed subcellular RNA sequencing in situ. Science 343, 1360-1363.
    Pubmed KoreaMed CrossRef
  29. Levsky J.M. and Singer R.H. (2003). Fluorescence in situ hybridization: past, present and future. J. Cell Sci. 116, 2833-2838.
    Pubmed CrossRef
  30. Lubeck E., Coskun A.F., Zhiyentayev T., Ahmad M., and Cai L. (2014). Single-cell in situ RNA profiling by sequential hybridization. Nat. Methods 11, 360-361.
    Pubmed KoreaMed CrossRef
  31. Mathot L., Kundu S., Ljungstrom V., Svedlund J., Moens L., Adlerteg T., Falk-Sorqvist E., Rendo V., Bellomo C., and Mayrhofer M., et al. (2017). Somatic ephrin receptor mutations are associated with metastasis in primary colorectal cancer. Cancer Res. 77, 1730-1740.
    Pubmed CrossRef
  32. McFaline-Figueroa J.L., Hill A.J., Qiu X., Jackson D., Shendure J., and Trapnell C. (2019). A pooled single-cell genetic screen identifies regulatory checkpoints in the continuum of the epithelial-to-mesenchymal transition. Nat. Genet. 51, 1389-1398.
    Pubmed KoreaMed CrossRef
  33. Moffitt J.R. and Zhuang X. (2016). RNA imaging with multiplexed error-robust fluorescence in situ hybridization (MERFISH). Methods Enzymol. 572, 1-49.
    Pubmed KoreaMed CrossRef
  34. Moncada R., Chiodin M., Devlin J.C., Baron M., Hajdu C.H., Simeone D., and Yanai I. (2018). Integrating single-cell RNA-Seq with spatial transcriptomics in pancreatic ductal adenocarcinoma using multimodal intersection analysis. BioRxiv . https://doi.org/10.1101/254375.
    CrossRef
  35. Moor A.E., Harnik Y., Ben-Moshe S., Massasa E.E., Rozenberg M., Eilam R., Bahar Halpern K., and Itzkovitz S. (2018). Spatial reconstruction of single enterocytes uncovers broad zonation along the intestinal villus axis. Cell 175, 1156-1167.e15.
    Pubmed CrossRef
  36. Nichterwitz S., Chen G., Aguila Benitez J., Yilmaz M., Storvall H., Cao M., Sandberg R., Deng Q., and Hedlund E. (2016). Laser capture microscopy coupled with Smart-seq2 for precise spatial transcriptomic profiling. Nat. Commun. 7, 12139.
    Pubmed KoreaMed CrossRef
  37. Nitzan M., Karaiskos N., Friedman N., and Rajewsky N. (2019). Gene expression cartography. Nature 576, 132-137.
    Pubmed CrossRef
  38. Noël F., Massenet-Regad L., Carmi-Levy I., Cappuccio A., Grandclaudon M., Trichot C., Kieffer Y., Mechta-Grigoriou F., and Soumelis V. (2020). ICELLNET: a transcriptome-based framework to dissect intercellular communication. BioRxiv . https://doi.org/10.1101/2020.03.05.976878.
    CrossRef
  39. Okamura-Oho Y., Shimokawa K., Takemoto S., Hirakiyama A., Nakamura S., Tsujimura Y., Nishimura M., Kasukawa T., Masumoto K.H., and Nikaido I., et al. (2012). Transcriptome tomography for brain analysis in the web-accessible anatomical space. PLoS One 7, e45373.
    Pubmed KoreaMed CrossRef
  40. Oktay M.H., Lee Y.F., Harney A., Farrell D., Kuhn N.Z., Morris S.A., Greenspan E., Mohla S., Grodzinski P., and Norton L. (2015). Cell-to-cell communication in cancer: workshop report. NPJ Breast Cancer 1, 15022.
    Pubmed KoreaMed CrossRef
  41. Patel A.P., Tirosh I., Trombetta J.J., Shalek A.K., Gillespie S.M., Wakimoto H., Cahill D.P., Nahed B.V., Curry W.T., and Martuza R.L., et al. (2014). Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma. Science 344, 1396-1401.
    Pubmed KoreaMed CrossRef
  42. Pollen A.A., Nowakowski T.J., Shuga J., Wang X., Leyrat A.A., Lui J.H., Li N., Szpankowski L., Fowler B., and Chen P., et al. (2014). Low-coverage single-cell mRNA sequencing reveals cellular heterogeneity and activated signaling pathways in developing cerebral cortex. Nat. Biotechnol. 32, 1053-1058.
    Pubmed KoreaMed CrossRef
  43. Proserpio V., Piccolo A., Haim-Vilmovsky L., Kar G., Lonnberg T., Svensson V., Pramanik J., Natarajan K.N., Zhai W., and Zhang X., et al. (2016). Single-cell analysis of CD4+ T-cell differentiation reveals three major cell states and progressive acceleration of proliferation. Genome Biol. 17, 103.
    Pubmed KoreaMed CrossRef
  44. Raj A., Rifkin S.A., Andersen E., and van Oudenaarden A. (2010). Variability in gene expression underlies incomplete penetrance. Nature 463, 913-918.
    Pubmed KoreaMed CrossRef
  45. Raj A., van den Bogaard P., Rifkin S.A., van Oudenaarden A., and Tyagi S. (2008). Imaging individual mRNA molecules using multiple singly labeled probes. Nat. Methods 5, 877-879.
    Pubmed KoreaMed CrossRef
  46. Rodriques S.G., Stickels R.R., Goeva A., Martin C.A., Murray E., Vanderburg C.R., Welch J., Chen L.M., Chen F., and Macosko E.Z. (2019). Slide-seq: a scalable technology for measuring genome-wide expression at high spatial resolution. Science 363, 1463-1467.
    Pubmed KoreaMed CrossRef
  47. Rouhanifard S.H., Mellis I.A., Dunagin M., Bayatpour S., Jiang C.L., Dardani I., Symmons O., Emert B., Torre E., and Cote A., et al. (2018). ClampFISH detects individual nucleic acid molecules using click chemistry-based amplification. Nat. Biotechnol . 2018 Nov 12 [Epub]. https://doi.org/10.1038/nbt.4286.
    Pubmed KoreaMed CrossRef
  48. Roy S. and Kornberg T.B. (2015). Paracrine signaling mediated at cell-cell contacts. Bioessays 37, 25-33.
    Pubmed KoreaMed CrossRef
  49. Salmén F., Vickovic S., Larsson L., Stenbeck L., Vallon-Christersson J., Ehinger A., Häkkinen J., Borg Å., Frisén J., and Ståhl P.L., et al. (2018). Multidimensional transcriptomics provides detailed information about immune cell distribution and identity in HER2+ breast tumors. BioRxiv . https://doi.org/10.1101/358937.
    CrossRef
  50. Satija R., Farrell J.A., Gennert D., Schier A.F., and Regev A. (2015). Spatial reconstruction of single-cell gene expression data. Nat. Biotechnol. 33, 495-502.
    Pubmed KoreaMed CrossRef
  51. Shah S., Lubeck E., Zhou W., and Cai L. (2016). In situ transcription profiling of single cells reveals spatial organization of cells in the mouse hippocampus. Neuron 92, 342-357.
    Pubmed KoreaMed CrossRef
  52. Shalek A.K., Satija R., Shuga J., Trombetta J.J., Gennert D., Lu D., Chen P., Gertner R.S., Gaublomme J.T., and Yosef N., et al. (2014). Single-cell RNA-seq reveals dynamic paracrine control of cellular variation. Nature 510, 363-369.
    Pubmed KoreaMed CrossRef
  53. Silberstein L., Goncalves K.A., Kharchenko P.V., Turcotte R., Kfoury Y., Mercier F., Baryawno N., Severe N., Bachand J., and Spencer J.A., et al. (2016). Proximity-based differential single-cell analysis of the niche to identify stem/progenitor cell regulators. Cell Stem Cell 19, 530-543.
    Pubmed KoreaMed CrossRef
  54. Skelly D.A., Squiers G.T., McLellan M.A., Bolisetty M.T., Robson P., Rosenthal N.A., and Pinto A.R. (2018). Single-cell transcriptional profiling reveals cellular diversity and intercommunication in the mouse heart. Cell Rep. 22, 600-610.
    Pubmed CrossRef
  55. Soldatov R., Kaucka M., Kastriti M.E., Petersen J., Chontorotzea T., Englmaier L., Akkuratova N., Yang Y., Haring M., and Dyachuk V., et al. (2019). Spatiotemporal structure of cell fate decisions in murine neural crest. Science 364, eaas9536.
    Pubmed CrossRef
  56. Stahl P.L., Salmen F., Vickovic S., Lundmark A., Navarro J.F., Magnusson J., Giacomello S., Asp M., Westholm J.O., and Huss M., et al. (2016). Visualization and analysis of gene expression in tissue sections by spatial transcriptomics. Science 353, 78-82.
    Pubmed CrossRef
  57. Tiklova K., Bjorklund A.K., Lahti L., Fiorenzano A., Nolbrant S., Gillberg L., Volakakis N., Yokota C., Hilscher M.M., and Hauling T., et al. (2019). Single-cell RNA sequencing reveals midbrain dopamine neuron diversity emerging during mouse brain development. Nat. Commun. 10, 581.
    Pubmed KoreaMed CrossRef
  58. Trapnell C., Cacchiarelli D., Grimsby J., Pokharel P., Li S., Morse M., Lennon N.J., Livak K.J., Mikkelsen T.S., and Rinn J.L. (2014). The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells. Nat. Biotechnol. 32, 381-386.
    Pubmed KoreaMed CrossRef
  59. van den Brink S.C., Alemany A., van Batenburg V., Moris N., Blotenburg M., Vivie J., Baillie-Johnson P., Nichols J., Sonnen K.F., and Martinez Arias A., et al. (2020). Single-cell and spatial transcriptomics reveal somitogenesis in gastruloids. Nature . 2020 Feb 19 [Epub]. https://doi.org/10.1038/s41586-020-2024-3.
    Pubmed CrossRef
  60. Vickovic S., Eraslan G., Salmén F., Klughammer J., Stenbeck L., Äijö T., Bonneau R., Bergenstråhle L., Fernandéz Navarro J., and Gould J., et al. (2019). High-definition spatial transcriptomics for in situ tissue profiling. Nat. Methods 16, 987-990.
    Pubmed KoreaMed CrossRef
  61. Wang G., Moffitt J.R., and Zhuang X. (2018a). Multiplexed imaging of high-density libraries of RNAs with MERFISH and expansion microscopy. Sci. Rep. 8, 4847.
    Pubmed KoreaMed CrossRef
  62. Wang X., Allen W.E., Wright M.A., Sylwestrak E.L., Samusik N., Vesuna S., Evans K., Liu C., Ramakrishnan C., and Liu J., et al. (2018b). Three-dimensional intact-tissue sequencing of single-cell transcriptional states. Science 361. eaat5691.
    Pubmed KoreaMed CrossRef
  63. Xia C., Fan J., Emanuel G., Hao J., and Zhuang X. (2019). Spatial transcriptome profiling by MERFISH reveals subcellular RNA compartmentalization and cell cycle-dependent gene expression. Proc. Natl. Acad. Sci. U. S. A. 116, 19490-19499.
    Pubmed KoreaMed CrossRef
  64. Xue Y., Liu D., Cui G., Ding Y., Ai D., Gao S., Zhang Y., Suo S., Wang X., and Lv P., et al. (2019). A 3D atlas of hematopoietic stem and progenitor cell expansion by multi-dimensional RNA-seq analysis. Cell Rep. 27, 1567-1578.e5.
    Pubmed CrossRef


Current Issue

31 July 2020 Volume 43,
Number 7, pp. 591~670

This Article


Cited By Articles
  • CrossRef (0)

Social Network Service
Services

Indexed in

  • Science Central
  • CrossMark