Omics technologies, which provide a comprehensive, global assessment of a set of molecules,1 offer a novel approach for the study of biological activity in life-sciences research. In recent years, the omics field has advanced considerably due to developments in the high-throughput analysis of biological molecules. Transcriptomic technologies in particular have provided rich insights into cellular characteristics, and measure the complete set of RNA transcripts produced by the genome under specific circumstances. Indeed, characterizing the transcriptome of individual cells is fundamental to improve our understanding of complex biological systems. Over the last decade, the rapid development of single-cell sequencing technology has greatly accelerated biomedical research, allowing scientists to overcome the major challenge of heterogeneity within biological samples and leading to the publication of a series of single-cell transcriptome atlases.2–4
However, most single-cell RNA sequencing (scRNA-seq) approaches involve isolating cells from their original position, meaning that spatial information is lost during transcriptome profiling. This omission makes it impossible to examine how a given cellular state interacts with neighboring cells or the surrounding extracellular matrix.2,5,6 To make up for this deficit, a new technology known as spatially resolved transcriptomics (SRT) aims to determine gene expression profiles while preserving information about the spatial context of the tissue. This technique is of great significance, since identifying the specific spatial location and organization of different cell types in multicellular tissues or organs is the foundation for research into cellular biological functions. Transcriptomics that is accompanied by spatial information will undoubtedly enable us to better understand the organization of cells and tissues, and how this organization influences biological function.7 The application of SRT has also opened up new areas of spatially resolved omics, including genomic, transcriptomic, proteomic and potentially other omic data with retained positional information. The ability of SRT to change the way we understand complex tissues earned it the title of Method of the Year 2020 in Nature Methods.8 In this review, we will first summarize and compare the major methods of SRT (Table 1), highlight some current applications of SRT, and finally discuss our perspectives on the potential future applications of this method.
Table 1 -
Major SRT methods.
|Major methods of spatially resolved transcriptomics
||In situ hybridization via complementary base pairing
Without amplification bias
||Optical SRM and combinatorial labeling
||Improved throughput of smFISH
||High-resolution imaging of RNA structure and location
||Specific materials needed
||Sequential barcoding scheme
||Dramatically improved throughput
||Large number of designed probes needed
||Error-robust encoding schemes achieved by introducing a certain Hamming distance
||Further improved throughput
Error detection and correction
|Large number of designed probes needed
||Cross-correlating the images from 2 of out of a total of 3 rounds of hybridization
||Decode high copy number RNAs
||Large number of designed probes needed
||Using secondary probes to dilute target transcripts
||Diluting the density of mRNA
||Large number of designed probes needed
|Sequencing-based methods that do not validate cell types
||Cryosectioning of tissue and performing RNA-seq on individual thin sections
||Genome-wide RNA tomography of tissues or organs
||Often results in 1D data along one axis
Often only useful for elongated samples
|DNA-barcoded beads surface
||Near-cellular spatial resolution
|No tissue imaging
||Labeling nuclei using unmodified DNA oligos on glass slides
Large capture area
|Imaging restricted to nuclei
||DNA nanoball (DNB) barcoded solid-phase RNA capture
||Nano-level resolution (Subcellular)
Large capture area
|Imaging restricted to nuclei
|Sequencing-based methods that validate cell types
||Applying laser energy to remove the cells of interest under microscopy
||Speed, precision, and versatility
Suitable for a small number of cells
Damages adjacent cells
||Combination of LCM and scRNA-seq
||Speed and precision
Suitable for a small number of cells
Damages adjacent cells
||Cells sorted to scRNA-seq after
fluorescent markers photoactivated by 2-photon laser scanning microscopy
|Single-cell level RNA-seq
Reliant on transgenic mice
Not suitable for clinical samples
|Barcoded solid-phase RNA capture
||Combining histological staining
Easy to perform
Requires no special instruments
Available commercial kit
|Low spatial resolution
Limited capture area
||Barcoded solid-phase RNA capture
||Improved spatial resolution
||Limited capture area
||Microfluidic-based method to deliver barcodes to the surface of a tissue slide
||Co-mapping of mRNAs and proteins
||Limited capture area
Spatial resolution needs further improvement
||Barcoded solid-phase RNA capture
||Nano-level resolution (Subcellular)
Large capture area
|Needs more research for demonstration
|Sequencing-based methods in vivo
||Photoactivatable mRNA capture in live cells and tissues
||Compatible with live and intact tissue, mRNA captured with precise spatial resolution, with little bias from RNA contamination or experimentally-related injury
|In situ sequencing methods
||Based on padlock probes and RCA, fixing the cDNA fragments to the cellular protein matrix by amine cross-linker
||Preserve the tissue architecture for RNA localization studies, localize subcellular RNA transcriptome-wide
||Specific instruments required
Random priming is inefficient and may introduce bias
Large number of designed probes needed
||Based on padlock probes and RCA, integrating hydrogel-tissue chemistry, targeted signal amplification, and In situ SEDAL sequencing
||3D volumes of intact tissue, no intrinsic limit to the number of genes
||Large number of designed probes needed
Workflow is complicated
|Methods for FFPE tissue
||Combination of LCM and Smart-3SEQ
||Suitable for small samples or those whose RNA is degraded
||Needs more research for demonstration
||Probes capture RNA after tissue decrosslinked
||For clinical samples from biobanks, immunological staining can be performed
||Needs more research for demonstration
corrFISH = correlation FISH, FFPE = formalin-fixed paraffin-embedded, FISH = fluorescence in situ hybridization, HDST = high-definition spatial transcriptomics, LCM = laser-capture microdissection, MERFISH = multiplexed error-robust FISH, scRNA-seq = single-cell RNA-sequencing, smFISH = single-molecule FISH, seqFISH = sequential barcoded FISH, SRM = super-resolution microscopy, ST = spatial transcriptomics, TIVA= transcriptome in vivo analysis, DBiT-seq =deterministic barcoding in tissue for spatial omics sequencing, ExFISH=expansion microscopy of fluorescence in situ hybridization, FISSEQ=fluorescent in situ RNA sequencing, Geo-seq=geographical position sequencing, RCA= rolling circle amplification..
2. METHODS AND ADVANCES
The core, underlying principle of SRT, involves the acquisition of gene expression profiles while retaining spatial information. Broadly speaking, SRT methods can be classified into 3 groups: imaging-based methods, sequencing-based methods, and in situ sequencing methods. Imaging-based methods are carried out using a microscope, and the transcriptomes are read out through fluorescence in situ hybridization (FISH). Sequencing-based methods are performed using next-generation sequencing technology, and RNA with spatial information is captured and subsequently sequenced. In situ sequencing methods are the combination of imaging- and sequencing-based methods. These methods have their intrinsic strengths and limitations, with differences in target coverage, resolution, and throughput. Like all methods, they must be appropriately matched to the biological question at hand.
2.1. Imaging-based methods
Spatial information about protein expression can be obtained from traditional immunohistochemistry (IHC) approaches, which are mediated by the specific binding reaction between antigens and antibodies. In a variation of the IHC method, spatial information about gene expression can be obtained by identifying the mRNA in cells using in situ hybridization (ISH), which involves a hybridization reaction between labeled nucleotide probes and complementary targeting RNA molecules. The single-molecule FISH (smFISH) method is the gold standard for detecting transcripts, and efforts to develop multiplex FISH at the transcriptome level are underway. This would allow FISH to answer more complex biological questions.
2.1.1. smFISH: the gold standard for transcript detection
smFISH combines quantitative FISH with digital imaging microscopy to sensitively detect single RNA molecules using probes. To amplify the signals, probes are often labeled with multi-fluorochromes per molecule.9 However, the disadvantages, which include difficulty in acquiring heavily labeled oligonucleotides, the easy loss of coupling fluorophores and self-quenching, highly variable signals associated with multi-fluorophores limit the system to a narrow range of applications. To overcome these drawbacks, singly labeled probes, which are easier to generate and purify, were introduced to accurately detect individual mRNA molecules with uniform signals10,11 (Fig. 1A).
smFISH can not only visualize the transcripts in isolated cells or in situ tissues at a high resolution and with high hybridization efficiency, but also provides information about the distribution of subcellular RNA within cells. This feature means that extra analyses concerning subcellular mRNA localization and its functional stages can be achieved, as mRNA localization is currently believed to contribute to the localized regulation of gene expression.12–15 In addition, smFISH can precisely measure the copy numbers of specific RNAs without amplification bias and natural fluctuations in gene expression can be quantitatively measured, permitting the elucidation of the mechanisms that regulate gene expression fluctuations and their role in a variety of biological processes. These features make smFISH a robust tool, and it is the standard method used for detecting transcripts.
2.1.2. Multiplexed FISH: improving target coverage
Many biological questions require transcription information about multiple genes, so efforts have been made to address the challenge of the limited number of RNA species that be simultaneously detected in cells using FISH. Super-resolution microscopy (SRM) can resolve a large number of mRNAs in single cells by labeling mRNAs with unique combinations of fluorophores. Although this expanded the scalability of FISH from a few genes to about 30, it was hard to make further progress due to the limited number of fluorophores.16 Later, a sequential barcoding scheme was introduced to multiplex different mRNAs, called sequential barcoded FISH (seqFISH).17 Here, the mRNA is barcoded by sequential rounds of hybridization, imaging, and probe stripping. The number of RNA species scales as FN, where F is the number of fluorophores and N is the number of rounds of hybridization. Thus, the scalability of seqFISH can be dramatically improved by increasing the number of rounds of hybridization, despite the limited number of fluorophores. In this manner, the entire human transcriptome can be covered, including both mRNAs and non-coding RNAs (Fig. 1B). However, N rounds of hybridizations bring not only exponentially increased throughput, but also exponentially increased detection error rates. When N increases, the error rates increase rapidly, and the majority of RNA molecules are likely to be misidentified after 16 rounds of hybridization.18 Therefore, multiplexed error-robust FISH (MERFISH) method was designed. MERFISH is based on combinatorial labeling and sequential imaging, and uses error-robust encoding schemes to correct labeling and detection errors. The method introduces a certain Hamming distance to encode RNAs, and the minimum Hamming distance is 4 (HD4 code), meaning that at least 4 bits must be read incorrectly to change one codeword into another. As a result, every single-bit error produces a word that is uniquely close to a single codeword, allowing such errors to be detected and corrected. In addition, MERFISH improves the multiplexity of transcripts that can be obtained in single cells to 100 to 1000 through the use of binary sequential barcodes19,20 (Fig. 1C).
2.1.3. Optical crowding: the major challenge for transcriptome-level profiling
Although seqFISH makes scaling to the level of the genome theoretically possible, transcriptome-level profiling in cells is hindered by optical crowding due to the high density of global transcripts, even when combined with SRM. Some progress has been made to improve the resolution of microscopy, with expansion microscopy (ExM) using a swellable polyelectrolyte gel that can expand physically to magnify the sample. RNAs are covalently attached to the gel, a small molecule linker, and FISH (ExFISH) is performed in cultured cells and intact tissues with high yield and specificity.21 Another approach called correlation FISH (corrFISH) was used in seqFISH to read out individual RNA species and barcodes. The principle is simple: RNA species are decoded such that each RNA species appears in only 2 out of a total of 3 rounds of hybridization. By cross-correlating the images from the 2 rounds of hybridization, only the specific RNA species will generate a positive correlation. Although some algorithms have been adapted due to the need for a high degree of temporal resolution, corrFISH still represents a robust method to decode high copy number RNAs in highly multiplexed seqFISH experiments, using conventional fluorescence microscopy.22
Another way to address the issue of molecular crowding is to dilute the density of the mRNA. Thus, seqFISH+ labels the secondary probes with fluorophores rather than the primary probes, which provide target sites for the secondary probes. The primary probes divide the transcript set into 20 subsets, and the secondary probes, which are labeled with three colors, ultimately make it possible to visualize one-sixtieth of the transcripts per image in each round of hybridization (Fig. 1D). seqFISH+ achieves transcriptome-level profiling not only for cultured cells but also in tissues, so that unbiased identification of cells alongside their spatial organization can be realized.23
2.2. Sequencing-based methods
With the development of genetic techniques, next-generation sequencing technology has become a powerful tool for transcriptome analysis. scRNA-seq can precisely characterize cell types and states at the molecular level, without bias, based on global gene expression profiles.24 Generally speaking, sequencing-based methods are performed to capture or barcode the RNA before reverse transcription to ensure each transcript can be mapped to its original spatial spot, or are used to gain gene expression information from cells in a certain spatial context. RNA-seq can be conducted either ex situ or in situ.
2.2.1. Sequencing-based methods that do not validate cell types
Some methods do not provide the spatial resolution of microscopy-based techniques as they do not predict or validate cell types using IHC or ISH. Examples of such methods include tomo-seq and Slide-seq.25,26 Tomo-seq is a serial microtomy-based sequencing method, which cryosections the tissue of interest and performs RNA-seq on individual thin sections (Fig. 2A). RNA extraction is performed for each individual section and the RNA is barcoded with section-specific primers. This approach can provide genome-wide RNA tomography of tissues or organs, but it is often only possible to cryosection tissue in 1 direction, resulting in 1-dimensional data along 1 axis of the body. For this reason, tomo-seq is best suited for elongated samples. Cells with complex organizations in the same individual sections cannot be analyzed well, and often several identical samples are required to perform cryosectioning in different directions along the main body axes to achieve 2- and 3-dimensional resolution.27
Unlike tomo-seq, Slide-seq is a method for transferring RNA from tissue sections onto a surface covered in DNA-barcoded beads with known positions, allowing the locations of the RNA to be inferred by sequencing. It provides a scalable method for obtaining spatially resolved gene expression data at resolutions comparable to the sizes of individual cells.26 Slide-seqV2 is an optimized version of Slide-seq with improvements in library generation, bead synthesis, and array indexing to reach an RNA capture efficiency ~10-fold greater than Slide-seq.28 A new method, called sci-Space, adapts the low-cost sci-Plex procedure for labeling or “hashing” cell nuclei using unmodified DNA oligos to glass slides. The hashed oligos spotted onto the glass slides are transferred onto permeabilized nuclei (not intact cells) in freshly frozen tissue sections. As the nuclei are stained with 4',6-diamidino-2-phenylindole (DAPI) and the spatial coordinates are consequently imaged during oligo transfer, sci-Space can retain a single-cell resolution.29 Obviously, the main limitation of all these methods is that the cell types are purely identified by RNA-seq datasets, which lacks the cell-validation process that microscopy-based techniques provide.
2.2.2. Sequencing-based methods that validate cell types
Although single-cell RNA-seq has the power to define cell types or populations based on gene expression programs, microscopy-based techniques such as IHC or ISH—which validate candidate genes with high cell type specificity—remain the gold standard to define cell type. Thus, a series of sequencing-based methods combined with imaging have been developed to aid the study of solid cell types.
Laser-capture microdissection (LCM) is an optional method to obtain targeted cell subgroups or even single cells quickly and precisely under the microscope. LCM directly harvests cells of interest or isolates specific cells by cutting away unwanted cells via infrared capture systems or ultraviolet cutting systems30 (Fig. 2B). Histological staining or rapid antibody staining are performed on the tissues before LCM so that information about the types of cells captured is validated. The key advantages of LCM are that the positional information of cells is maintained and tissues are not dissected; in this way, subsequent analysis by genomics, transcriptomics, or proteomics can be performed on the LCM-captured cells.31–33 The extracted nucleic acids can be amplified when the captured regions of interest are small, overcoming the requirement for relatively large numbers of cells. This advantage makes it possible to use LCM at the single-cell level, meaning that it is especially useful for scarce tissues and rare cell types. Geo-seq is the combination of LCM and scRNA-seq technology.33 It is also practical to combine LCM with Smart-seq2, a high-throughput scRNA-seq technology with improved sensitivity, accuracy, and full-length coverage across transcripts, making LCM plus Smart-seq2 a promising prospect for SRT.31,34,35 However, the main obstacles to the wide application of LCM SRT are its low throughput and potential RNA contamination from other cells, as well as the inevitable damage to adjacent tissue. This means that it is difficult to study the interplay between individual cells and the microenvironment.
One method with similarities to LCM coupled with Smart-seq2 is NICHE-seq, which couples with scRNA-seq. This technique combines photoactivatable fluorescent markers, 2-photon laser scanning microscopy and flow cytometry-based fluorescence-activated cell sorting (FACS), and is coupled to massively parallel scRNA-seq (Fig. 2C). NICHE-seq requires transgenic mice that ubiquitously express a photoactivatable green fluorescent protein (PA-GFP). Two-photon irradiation is used to activate PA-GFP and in situ labeling. After tissue dissociation, GFP-labeled cells can be sorted to perform scRNA-seq; thus, this method combines information about the cell’s transcriptional state and spatial information.36 As this method requires transgenic mice expressing GFP, however, it is not suitable for clinical samples.
Another method based on barcoded glass slides, called spatial transcriptomics (ST), also validates cell types by histological staining and imaging.37 RNA is spatially captured from tissue sections and barcoded via a microarray on a glass slide. Several primers with unique barcodes are placed on the microscopic glass slide to create a microarray, and a section of tissue is placed on top of the microarray surface. This tissue is then fixed, then stained and imaged, before being permeabilized to release RNA within the cells. In this way, transcripts from the tissue are allowed to meet with the immobilized cDNA synthesis primers in a reverse transcription reaction. Consequently, the cDNA library is sequenced and the data are visualized together with a high-resolution histological image of the tissue section. Each slide consists of about 1000 spots, and ST is both easy to perform and requires no special instruments. The method has been acquired and optimized by 10x Genomics and a commercial kit is available (Fig. 2D).
DBiT-seq is a microfluidic-based method that delivers barcodes to the surface of a tissue slide to allow for spatial omic sequencing. This method enables the co-mapping of mRNAs and proteins on a slide containing formaldehyde-fixed tissue38 (Fig. 3A). Immunofluorescence can be performed on the same tissue slide, so that proteomics can be integrated with transcriptomics. The capacity for integration is of great significance for multi-omics research.
2.2.3. Sequencing-based methods for use in vivo
Current methods that capture RNA, such as LCM, cannot isolate mRNA from individual in vivo cells without damaging adjacent tissue, making it difficult to assess the influence of the microenvironment on the transcriptome. To overcome this problem, a method was engineered to extract mRNA from live cells or tissues at the spatial resolution of a single cell. This method is called transcriptome in vivo analysis (TIVA); it uses a photoactivatable biotinylated tag, known as a TIVA-tag, that can penetrate the cell membrane by virtue of a cell-penetrating peptide. Once inside the cell, photoactivation is used to cleave the linkers between the poly(U) sequence and two poly(A) stretches to enable the poly(U) sequence to bind the poly-A tail of mRNAs (Fig. 3B). With this method, mRNA can be noninvasively captured within live cells and intact tissues, so it provides a useful tool to explore the transcriptomes of single cells in the context of their natural microenvironment. However, the method is low throughput with regard to the numbers of cells that can be analyzed.
2.2.4. Enhanced resolution and a larger area permit the acquisition of more detail
Solid RNA-barcoding is favorable in terms of its untargeted features and high coverage compared with imaging-based methods. Efforts have been made to improve the spatial resolution and capture area of this method, as a higher resolution means the precise location of the transcripts within the tissues, cells or even subcellular organelles can be determined.
ST can achieve a resolution of only 100 μm (which is not at the single-cell level), which means that the 100 μm spatial spots, with a center-to-center distance of 200 μm, typically cover 5 to 100 cells each, depending on tissue type and region.37 Visium technology from 10X Genomics has improved the ST method by reducing the diameter of the capture area from 100 to 55 μm, increasing the capture area to 42.25 mm2 (6.5 mm × 6.5 mm), and decreasing workflow duration (Fig. 2D). Nevertheless, the spatial resolution of 10X Genomics’ Visium is still insufficient to resolve single cells in most samples.
Slide-seq and DBiT-seq both have improved spatial resolution to 10 μm,26,38 but this is still a near-cellular resolution and cannot be described as single-cell resolution. High-definition spatial transcriptomics (HDST), an optimized version of ST, was designed to improve ST resolution to the single-cell level. In this method, hundreds of thousands of barcodes are randomly deposited onto a slide at 2 μm resolution, and their positions are decoded by sequential hybridization and error-correcting strategy.39 One method, called spatiotemporal enhanced resolution omics-sequencing (Stereo-Seq), involves the deposition of a DNA nanoball containing random barcoded sequences to form the chip. Each spot is approximately 220 nm in diameter, with a center-to-center distance of 500 nm.40 Another method, known as Seq-Scope, has a center-to-center resolution of 0.5 to 0.8 μm (0.6 μm on average).41 The nano-level distances of these methods are far less than the diameter of single cells, giving them the high resolution needed to determine the subcellular location of transcripts.
2.3. In situ sequencing methods
The combination of imaging and sequencing processes typically involves in situ sequencing methods. In situ sequencing methods sequence cDNA amplicons based on padlock probes and rolling circle amplification.42,43 The 2 ends of the padlock probe can be ligated to form a circle-like structure upon hybridization, which then can be amplified using an isothermal DNA polymerase. Amplification causes copies of the padlock probe sequence to accumulate and form cDNA nanoballs. An amplified signal can be obtained by visualizing the cDNA nanoballs using fluorescently labeled detection oligonucleotides, which are complementary to the padlock probe.
FISSEQ, a highly multiplexed subcellular in situ RNA sequencing method, uses random hexamer-primed reverse transcription in the presence of aminoallyl dUTP. Aminoallyl dUTP then fixes the cDNA fragments to the cellular protein matrix by a non-reversible amine cross-linker.44 The templates are circularized after degrading the RNA, and subsequently amplified using rolling circle amplification primers complementary to the adapter sequence introduced by RT. The amplicons in the cells are then ready for sequencing-by-ligation (SOLiD sequencing) at room temperature, and imaging on a confocal microscope. To control the signal density of FISSEQ, a partition sequencing strategy was devised using pre-extended sequencing primers.44,45 With 1 base extended, one can randomly sample amplicons at 1/4th of the original density, and with 2 bases extended sampling at 1/16th of the density is possible. One notable highlight of FISSEQ is that it preserves the tissue architecture for RNA localization studies, and localizes subcellular RNA transcriptome-wide. However, the random priming component is inefficient and may introduce bias via uneven binding efficiency, which is due to differences in base composition.
While SRT can be performed on tissues or sections on a slide, applying in situ sequencing methods to 3-dimensional volumes of intact tissue is not easy. To achieve this goal, a method known as spatially resolved transcript amplicon readout mapping (STARmap) integrates hydrogel-tissue chemistry, targeted signal amplification and in situ sequencing.46 The method bypasses the reverse transcription step and uses the SNAIL approach instead, which achieves specific amplification of nucleic acids via intramolecular ligation using a pair of custom primer and padlock probes. Only when both probes hybridize to the same RNA molecule can the padlock probe be circularized and rolling-circle-amplified to generate a DNA nanoball. The cDNA amplicons are modified and copolymerized in situ to embed within a hydrogel network. In situ SEDAL sequencing was devised specifically for STARmap, to sequence 5-base barcoded gene-specific identifier segments in SNAIL probes over 6 rounds of hybridization (Fig. 3C). The highlight of STARmap lies in its ability to be used in 3-dimensional in situ transcriptomics, which is undoubtedly a state-of-the-art SRT technique and serves to deepen our understanding of anatomy and transcriptomics. In parallel, there is no intrinsic limit to the number of genes or RNA species that can be simultaneously and quantitatively accessed with STARmap, since it can be adapted to longer sequencing lengths or higher gene numbers.
In conclusion, ISH is highly sensitive, and can detect transcripts at a subcellular resolution. Many powerful strategies have been designed to improve the target coverage of ISH to the transcriptome level and tackle the problem of optical crowding. However, such methods are typically limited by long image acquisition times and a complex probe-designing process. While next-generation sequencing technology enables high-throughput transcriptomic profiling, many sequencing-based methods are much easier to perform once suitable experimental platforms have been established. As spatially barcoded beads or spots can be arrayed more and more densely, some solid-barcoded methods can reach a cellular or subcellular resolution, and spatial information is more precisely retained. To some extent, in situ sequencing methods are a combination of imaging- and sequencing-based methods, sequencing gene identifiers in the designed primers via FISH. Although in situ sequencing methods are high throughput, the complexity of designing targeted probes cannot be ignored (Table 1).
Most of the in situ-based methods described above rely on the use of fresh tissues. However, fresh clinical material is often difficult to obtain for research purposes. To remove the barriers to clinical research, some methods have been designed to use formalin-fixed paraffin-embedded (FFPE) tissue, in which the RNA is degraded. Smart-3SEQ enables large gene expression profiling experiments to be conducted on even small amounts of total RNA, and effectively characterizes small samples extracted by LCM from FFPE tissue.47 In addition, 10X Genomics markets a commercially available system, called Visium FFPE, that enables SRT to be performed on samples from a biobank. Some studies have already successfully used 10x Visium FFPE on FFPE tissue.48
3. DATA ANALYSIS IN SRT
The additional dimension of spatial information brings not only a novel perspective on the transcriptome, but also significant challenges for data analysis due to increased data volume and complexity. Generally speaking, strategies for SRT data analysis include computational approaches specifically designed for SRT and integration methods (with bulk or single-cell RNA-seq data). Different analytical approaches have been developed for localized gene expression pattern identification, spatial decomposition, gene imputation, and cell–cell communication.49
Gene expression and spatial location are both of great relevance to biological functions. By analyzing the relationship between gene expression and spatial location, some approaches, including Trendsceek,50 SpatialDE,51 SPARK,52 SPARK-X,53 sepal,54 SpaGCN,55 and GLISS, are designed to identify localized gene expression patterns and spatially variable genes (SVGs) based on different algorithms. Another way to discover SVGs is spatial clustering. Unlike standard clustering methods for scRNA-seq, spatial clustering is for the capture locations. Methods like stLearn, MULTILAYER,56 BayesSpace,57 SC-MEB,58 and STAGATE can be applied for spatial clustering. As most of the sequencing-based methods do not guarantee that a single capture spot contains RNA from only 1 cell, a capture location is often a mixture of multiple cell types. Thus, some spatial decomposition algorithms, such as NNLS, spatialDWLS,59 SPOTlight,60 RCTD,61 DSTG, Tangram, and Cell2location attempt to infer the proportions of cell types in each spot, or score a spot for how strongly it corresponds to a single cellular subtype based on bulk or single-cell RNA-seq data. However, certain SRT methods, particularly those that are imaging-based, cannot provide deep transcriptomic information due to limited gene coverage. Gene imputation provides a chance to impute the missing genes from scRNA-seq data in order to improve the quality of SRT data. Methods that include gene imputation functions include Tangram,62 gimVI, Harmony,63 LIGER,64 Seurat,65 SpaGE,66 and stPlus.67 In addition, the reconstruction of spatial information for scRNA-seq data can be achieved by mapping, which leverages a small set of landmark genes to map the single cell back to its context. Methods that incorporate a mapping function include Harmony, LIGER, Seurat, SpaGE, DEEPsc,68 DistMap, SpaOTsc,69 novoSpaRc,70 and CSOmap.71 In addition, cell–cell interactions can be analyzed using some standard algorithms, based on the ligand–receptor interaction pairs from scRNA-seq data and a database of known ligand–receptor interactions. However, the loss of spatial information in scRNA-seq data may lead to inaccurate detection of ligand–receptor interactions, as cellular cross-talk is often spatially restricted. Therefore, SRT data might potentially evaluate the reliability of the ligand–receptor interactions computed from scRNA-seq. In this manner, integrating scRNA-seq with SRT data may accurately reveal specific cell subpopulations and their interactions. Approaches to study cell–cell interactions include SVCA,72 GCNG, NCEM, MISTy, stLearn, Squidpy, novoSpaRc, SpaOTsc, and DEEPsc. Among these methods, comprehensive benchmark studies are needed to help users select methods that best fit their data and hypotheses. According to a study benchmarking 16 integration methods, Tangram, gimVI, and SpaGE outperformed other integration methods when predicting the spatial distribution of RNA transcripts, while Cell2location,73 SpatialDWLS, and RCTD were the top-performing methods for the cell type deconvolution of spots.74
In conclusion, SRT and scRNA-seq data complement each other, since SRT methods cannot always maintain single-cell resolution with the same depth and whole transcriptome coverage of scRNA-seq. The integration of scRNA-seq and SRT data could improve cell type annotation, cell clustering, spatial decomposition, gene imputation, and spatial location reconstruction. A scRNA-seq database could form a solid foundation for SRT data analysis. However, the proper handling of spatial information obtained from SRT data remains a challenge in terms of data processing, and the establishment of proper protocols will be an essential step to reach the full potential of these methods.
4. APPLICATIONS FOR SRT
At present, SRT technologies are widely utilized in basic and clinical research, especially in neuroscience, cancer biology, and developmental biology. However, there are likely further applications for SRT that are yet to be fully utilized, including in the study of hematopoeisis. We reviewed the potential applications in this field in Figure 4.
4.1. Neuroscience: unraveling the organization of brain cell types and functions
The diversity of brain cell types and the complex organization of the brain have long been major challenges to our understanding of brain functions during physiological and pathological conditions.75 A census of brain cell types and their functions will pave the way for systematic research and, potentially, even novel clinical treatments. SRT has advantages over scRNA-seq due to the retained anatomical and spatial contexts, and various spatially resolved transcriptomic methods have been applied to map cell types in neuroscience. For example, an atlas of the mouse somatosensory cortex was built using osmFISH, a semiautomated version of smFISH. Based on scRNA-seq data, the study identified 31 clusters of cell types and 33 cell type-specific marker genes.76 Although multiplexing in osmFISH is relatively low, its sensitivity is higher than scRNA-seq and can provide single-cell-level spatial resolution for gene detection. Multiplexed FISH techniques, such as MERFISH and seqFISH, are applied when more genes need to be analyzed. The multiplexing of 250 genes using seqFISH identified distinct subregions of the hippocampus that contained different combinations of cell types.18 In addition, MERFISH combined with measurements of immediate early gene expression identified ~70 neuronal populations in the mouse hypothalamic preoptic area, and defined discrete cell populations activated by specific social behaviors.77 STARmap has also been used to reconstruct the 3D locations of cells in the murine primary visual cortex and medial prefrontal cortex; in this particular study, >1000 genes were simultaneously mapped in mouse brain sections at a single-cell resolution.46 Sequencing-based methods such as Slide-seq and 10X Genomics’ Visium have also been used to analyze the global transcriptome in both the mouse and human brain.26,37
The application of spatially resolved transcriptomic methods is highly likely to improve our understanding of neurological diseases. For example, SRT has already revealed the spatiotemporal dynamics of molecular pathology in amyotrophic lateral sclerosis.78 Furthermore, in a murine model of Alzheimer disease, ST and in situ sequencing were used to better characterize the dysregulated cellular network in the early and late phases of disease. The results of this study provided profound insights into the pathogenesis of Alzheimer disease, and also provided clues concerning potential therapeutic targets.79
4.2. Cancer: profiling the spatial heterogeneity of tumors and the tumor microenvironment
Heterogeneity both within and between tumors is found in most types of cancer.80,81 The expansion of tumor subclones with varying genetic alterations and the interactions between tumor cells and the tumor microenvironment presents a complicated situation that influences disease development.82 Profiling the spatial heterogeneity of tumors and the tumor microenvironment will provide insights into tumor progression and responses to treatment, and will therefore assist in the diagnosis and treatment of this condition. ST technology has been used to identify intertumoral and intratumoral heterogeneity in multiple types of cancer. Furthermore, changes in gene expression during prostate cancer and cutaneous malignant melanoma progression have revealed a detailed landscape of tumor progression and metastases.83,84 Several features of the tumor microenvironment in primary liver cancer, including stromal and immune cell distribution, tumor cluster interaction, cancer stem cell-niche diversity and tertiary lymphoid structure composition have been characterized.85 Coupling SRT with other technologies also improves its capacity to characterize tumors: for example, when SRT was combined with scRNA-seq, tissue architecture and interactions between cell subpopulations were revealed in pancreatic ductal adenocarcinomas.86 A tumor-specific keratinocyte population, immune infiltrates, and heterogeneity at tumor leading edges have also been identified in human squamous cell carcinoma using multimodal analysis, which included scRNA-seq and ST.87 The heterogeneity of breast cancer has also been characterized, and several diagnostic and prognostic markers proposed.88,89 Additionally, by integrating spatial gene expression and HE-stained images of breast tumor, a deep learning algorithm was developed to predict the spatially resolved transcriptome of a tissue from the images, enabling image-based screening for molecular biomarkers.90 To guide clinical treatment, biomarkers associated with beneficial PD-1 checkpoint blockade in non-small-cell lung cancer were identified using digital spatial profiling.91 The ST analysis of human bladder cancer also identified an N-Cadherin 2-expressing epithelial cell subpopulation which could be used to predict therapeutic response, which could guide treatment decisions.92 In conclusion, these findings provide novel insights into the complex ecosystem of cancer. and have the potential to improve individualized cancer treatments and drug discovery.
4.3. Developmental biology: characterizing the temporal and spatial expression blueprint of embryonic development
Embryonic development is a complex process where dynamic changes rapidly occur at the biomolecular level.93,94 Lineage tracing and fate mapping are powerful tools in developmental biology, providing a picture of the lineage hierarchy and linking the initial cell position to future fate.95 With the emergence of SRT, analysis of both the temporal and spatial aspects of embryonic development can be performed, and a temporal and spatial expression blueprint is eagerly anticipated. In the early embryo, gastrulation is an essential developmental event, leading to the formation of gastruloids, which are aggregates of embryonic stem cells.96,97 Applying Geo-seq to the germ layers from pre- to late-stage gastrulation enabled the molecular genealogy of tissue lineages and the continuum of pluripotency states in time and space to be revealed. Furthermore, the molecular determinants that drive lineage specification and tissue patterning were identified.98 Tomo-seq and scRNA-seq revealed genome-wide gene expression patterns in mouse gastruloids and embryos, and identified various embryonic cell types that were not previously known to be present in gastruloids. Based on their results, the study’s authors proposed that somitogenesis occurs in gastruloids.99 In addition, ISH has determined the gene expression pattern in gastruloids generated from embryonic stem cells, to validate an in vitro model of early anteroposterior organization during human development.100 Stereo-seq has also reconstructed the spatially resolved developmental trajectories of cell-fate transitions and molecular changes during zebrafish embryogenesis.101
Some research has also focused on the formation of organs in embryonic development. The large field-of-view and cellular resolution of Stereo-seq has been used to map the spatiotemporal transcriptomic dynamics during mouse organogenesis.102 ST was also used to study cardiac morphogenesis, and the spatiotemporal organ-wide gene expression and cell atlas of the developing human heart was established during 3 developmental stages.103 Morphogenesis of the human intestine across time, location, and cellular compartments has been charted by both scRNA-seq and 10X Genomics’ Visum.104 The gene expression profile of fetal livers from 8 to 17 weeks post-conception in humans has also been illustrated using 10X Genomics’ Visium.105
4.4. Hematopoiesis: shedding light on the distinct waves of developmental hematopoiesis, and making breakthroughs in adult bone marrow
Hematopoiesis is the process by which blood cells are continually replenished throughout the lifetime of an organism, by virtue of the self-renewal and differentiation of hematopoietic stem cells (HSCs) and the regulation of the local microenvironment, or niche.106 Distinct waves of hematopoiesis have been defined during embryogenesis, including both primitive and definitive hematopoiesis.107,108 SRT technologies have been applied in developmental hematopoiesis to explore the generation and expansion of HSCs, as well as the regulation of their niche. Geo-seq, combined with bulk and single-cell RNA-seq, was used to examine the caudal hematopoietic tissue of zebrafish (the counterpart of the mammalian fetal liver) and a detailed spatiotemporal transcriptome of hematopoietic stem and progenitor cells (HSPCs) and niche during HSPCs expansion was characterized.109 Meanwhile, LCM-seq has been used in the human embryonic aorta-gonad-mesonephros (AGM) region to identify secretory factors that promote human HSC development.110 Another study combined with 10X Genomic’s Visium and Stereo-seq identified novel HSC/MPP pocket-like units (HSC PLUS) composed of niche cells and enriched with growth factors; providing an essential resource for understanding HSC/MPP development in the fetal liver.111 In addition, 10X Genomic’s Visium was used to visualize HSC emergence in human embryonic tissues, including the liver, AGM, gut, and vitelline and umbilical vessels, validating the generation of AGM-like definitive HSPCs from human pluripotent stem cells.112
scRNA-seq technologies have been widely performed in the adult hematopoietic system, and have confirmed the heterogeneity of HSPCs.113 However, SRT is yet to be applied in this field. Profiling the SRT of adult bone marrow would be of great significance because this will provide an opportunity to study the interplay between HSPCs and their niche under both normal conditions and stress. Most SRT methods require an intact slice of tissue, which may be challenging for adult bone marrow due to the presence of calcified bone. Therefore, tissue preparation and cryosectioning need to be optimized. Furthermore, bone marrow is composed of various hematopoietic cells and non-hematopoietic cells, including mesenchymal stromal cells, pericytes, fibroblasts, and endothelial cells.114,115 Importantly, bone marrow hematopoietic cells are usually smaller than other cell types, and the cells are usually distributed evenly in the bone marrow, lacking any apparent structure in terms of anatomical or functional division.116,117 These intrinsic features of bone marrow impede the experimental applications of SRT, as well as associated data processing. LCM-seq has been applied to study adult bone marrow to reveal bone marrow niche organization, setting an excellent example for future research.118 RNA-barcoding methods with improved resolution, such as Slide-seq, HDST and stereo-seq, are recommended for bone marrow, since the cellular or subcellular resolution is the key to compensating for smaller cells and diverse cell types. Furthermore, different types of HSPCs share high similarities in transcriptomics. Thus, we propose that SRT data can be properly analyzed only when combined with scRNA-seq data. Of note, in the past decade, single-cell transcriptomic data of bone marrow hematopoietic cells and niche cells has been collected, serving as a supporting resource for SRT.114,115
5. CONCLUSIONS AND FUTURE PERSPECTIVES
SRT technologies provide an opportunity to uncover the molecular architecture of tissues. Methods that are more scalable, have a higher resolution and easier workflows are expected to be developed, and more results produced by SRT should be published soon. For the remainder of this review, we anticipate future developments that will unlock the full potential of SRT in both biological research and the clinic.
5.1. Realizing single-cell SRT
The correct association of detected mRNAs and single cells is a great challenge for both imaging- and sequencing-based methods. For most sequencing-based methods, such as Slide-seq and ST,26,35,37 the spatial information is retained in the form of differentially barcoded spots. Although some methods, including HDST, Stereo-seq and Seq-Scope, have achieved cellular or subcellular resolution,39–41 the spot boundaries do not correspond to the cell boundaries, meaning that transcripts captured in the same spot do not always come from a single cell. Furthermore, it cannot be guaranteed that the 10 μm-thin slices usually obtained during tissue sampling will contain a single layer of cells. This means that the transcripts of all the cells on the Z-axis or the vertical axis of the slide will be captured on the same spot. In imaging-based methods, the lack of cell boundaries makes it hard to associate individual mRNA molecules with the correct cells. In such methods, DAPI-stained nuclei may be the only method to determine cell outlines, especially when different cell types are very close together. In Sci-Space, which leverages sci-Plex to transfer DNA oligos into permeabilized nuclei, DAPI-stained nuclei and spatial coordinates of oligo transfer are co-registered. The nuclei from the tissue on the slide are also extracted to enable sci-RNA-seq.29 Although Sci-Space retains a single-cell resolution, it is hard to avoid contamination from adjacent cells. Moreover, the nuclei-assembling strategy may collapse when applied to tissues that contain irregularly shaped cells, such as neurons and glial cells in the brain. Regardless, we expect that more robust approaches will ultimately permit the characterization of the whole transcriptome of single cells. The elaboration of some cytomembrane-outlining strategies, such as immunofluorescent staining or auto-fluorescent reporter animals to visualize the plasma membrane, may help define cell boundaries and realize the potential of single-cell SRT. Existing rapid membrane imaging methods, such as fluorescent probe-based methods, have already been used to facilitate in situ plasma membrane imaging of neurons and erythrocytes in the complex environment of the brain.119 However, optimization of in situ plasma membrane imaging in combination with SRT still requires further work.
5.2. Integrating spatially resolved multi-omics to advance our understanding of biological systems
SRT is far more than just transcriptomics; rather, the technique is empowered by the retained positional information provided by in situ methods. The locations of cells and the microenvironment play an important role in uncovering cell-to-cell and cell-to-extracellular matrix interactions, and the integration of functional and spatial information is of great benefit in biological studies. Spatially resolved transcriptomic methods provide several invaluable methods that can also be adapted to suit certain other in situ omics. We expect that SRT will inspire the development of other spatially resolved omics beyond transcriptomics, such as in situ epigenomics, epitranscriptomics, proteomics, and metabolomics, which will help us better understand different aspects of cell function. Epigenomics is a reliable tool to map chromatin accessibility dynamics and higher-order chromatin structure, and adding spatial information to this would surely enable new levels of understanding of cell-fate decisions, identity, and function.120,121 Epitranscriptomics concerns a wide range of post-transcriptional RNA modifications, and spatial epitranscriptomics enables studies under specific microenvironments, broadening our understanding of the spatial regulation of gene expression.122 In parallel, metabolomics transcends genomics and proteomics, representing the most downstream metabolic stage. The metabolome is widely accepted to constitute a dynamic and sensitive measure of phenotype, placing metabolomics at the forefront of biomarker discovery.123 Several spatially resolved metabolomic technologies have already been applied in several fields, such as tumor metabolism, allowing metabolites and metabolic enzymes to be directly discovered and studied in their native state.124–127
On the other hand, cell states that are defined based on data from one type of omics may not be accurate and comprehensive, and cell states defined by different omics may not always be in accordance. We therefore expect the emergence of more spatially resolved multimodal technologies, both at the experimental level and in computational analysis, to be able to simultaneously profile multiple data types in the same cell. For example, DBiT-seq can capture both proteins and transcripts on the same sequencing slide, so that proteomics and transcriptomics can be simultaneously performed on the same cells or tissues.38 We are also appreciative of and enthusiastic about methods that validate cell types, not just to predict transcriptomic states, but more importantly to provide more possibilities for multi-omic research. For example, some sequencing-based methods with high spatial resolution, such as Stereo-seq, involve the rapid staining and imaging of just the tissue nuclei to avoid mRNA degradation.40 Simultaneously retaining the mRNA, protein or metabolites in the same cells may be very challenging. Proper tissue preparation and fixation are essential to preserve the tissue, and so existing workflows containing these steps will need to be optimized to suit multimodal research.
This work was supported by grants from the Ministry of Science and Technology of China (2021YFA1100900 and 2020YFE0203000), the National Natural Science Foundation of China (81730006, 81922002, 81861148029, and 81870086), CAMS Innovation Fund for Medical Sciences (2021-I2M-1-040 and 2021-I2M-1-019), Haihe Laboratory of Cell Ecosystem Innovation Fund (HH22KYZX0016) and Distinguished Young Scholars of Tianjin (19JCJQJC63400).
. Hasin Y, Seldin M, Lusis A. Multi-omics approaches to disease. Genome Biol 2017;18:83.
. Tang F, Barbacioru C, Wang Y, et al. mRNA-Seq whole-transcriptome analysis of a single cell. Nat Methods 2009;6:377–382.
. van Galen P, Hovestadt V, Wadsworth Ii MH, et al. Single-cell RNA-Seq reveals AML hierarchies relevant to disease progression and immunity. Cell 2019;176:1265–1281.e24.
. Han X, Zhou Z, Fei L, et al. Construction of a human cell landscape at single-cell level. Nature 2020;581:303–309.
. Zheng GX, Terry JM, Belgrader P, et al. Massively parallel digital transcriptional profiling of single cells. Nat Commun 2017;8:14049.
. Macosko EZ, Basu A, Satija R, et al. Highly parallel genome-wide expression profiling of individual cells using nanoliter droplets. Cell 2015;161:1202–1214.
. Moor AE, Itzkovitz S. Spatial transcriptomics: paving the way for tissue-level systems biology. Curr Opin Biotechnol 2017;46:126–133.
. Marx V, Method of the Year 2020: spatially resolved transcriptomics
. Nat Methods 2021;18:1.
. Femino AM, Fay FS, Fogarty K, Singer RH. Visualization of single RNA transcripts in situ. Science 1998;280:585–590.
. Levsky JM, Singer RH. Fluorescence in situ hybridization: past, present and future. J Cell Sci 2003;116:2833–2838.
. Raj A, van den Bogaard P, Rifkin SA, van Oudenaarden A, Tyagi S. Imaging individual mRNA molecules using multiple singly labeled probes. Nat Methods 2008;5:877–879.
. Buxbaum AR, Haimovich G, Singer RH. In the right place at the right time: visualizing and understanding mRNA localization. Nat Rev Mol Cell Biol 2015;16:95–109.
. Kejiou NS, Palazzo AF. mRNA localization as a rheostat to regulate subcellular gene expression. Wiley Interdiscip Rev RNA 2017;8:10.1002/wrna.1416.
. Holt CE, Bullock SL. Subcellular mRNA localization in animal cells and why it matters. Science 2009;326:1212–1216.
. Dermit M, Dodel M, Lee FCY, et al. Subcellular mRNA localization regulates ribosome biogenesis in migrating cells. Dev Cell 2020;55:298–313.e10.
. Lubeck E, Cai L. Single-cell systems biology by super-resolution imaging and combinatorial labeling. Nat Methods 2012;9:743–748.
. Lubeck E, Coskun AF, Zhiyentayev T, Ahmad M, Cai L. Single-cell in situ RNA profiling by sequential hybridization. Nat Methods 2014;11:360–361.
. Shah S, Lubeck E, Zhou W, Cai L. In situ transcription profiling of single cells reveals spatial organization of cells in the mouse hippocampus. Neuron 2016;92:342–357.
. Chen KH, Boettiger AN, Moffitt JR, Wang S, Zhuang X. Spatially resolved, highly multiplexed RNA profiling in single cells. Science 2015;348:aaa6090.
. Moffitt JR, Zhuang X. RNA imaging with multiplexed error-robust fluorescence in situ hybridization (MERFISH). Methods Enzymol 2016;572:1–49.
. Chen F, Tillberg PW, Boyden ES. Optical imaging. Expansion microscopy. Science 2015;347:543–548.
. Coskun AF, Cai L. Dense transcript profiling in single cells by image correlation decoding. Nat Methods 2016;13:657–660.
. Eng CL, Lawson M, Zhu Q, et al. Transcriptome-scale super-resolved imaging in tissues by RNA seqFISH. Nature 2019;568:235–239.
. Trapnell C. Defining cell types and states with single-cell genomics. Genome Res 2015;25:1491–1498.
. Kruse F, Junker JP, van Oudenaarden A, Bakkers J. Tomo-seq: a method to obtain genome-wide expression data with spatial resolution. Methods Cell Biol 2016;135:299–307.
. Rodriques SG, Stickels RR, Goeva A, et al. Slide-seq: a scalable technology for measuring genome-wide expression at high spatial resolution. Science 2019;363:1463–1467.
. Holler K, Junker JP. RNA tomography for spatially resolved transcriptomics
(Tomo-Seq). Methods Mol Biol 2019;1920:129–141.
. Stickels RR, Murray E, Kumar P, et al. Highly sensitive spatial transcriptomics at near-cellular resolution with Slide-seqV2. Nat Biotechnol 2021;39:313–319.
. Srivatsan SR, McFaline-Figueroa JM, Ramani V, Saunders L. Massively multiplex chemical transcriptomics at single-cell resolution. Science 2020;367:45–51.
. Espina V, Wulfkuhle JD, Calvert VS, et al. Laser-capture microdissection. Nat Protoc 2006;1:586–603.
. Nichterwitz S, Benitez JA, Hoogstraaten R, Deng Q, Hedlund E. LCM-Seq: a method for spatial transcriptomic profiling using laser capture microdissection coupled with PolyA-Based RNA sequencing. Methods Mol Biol 2018;1649:95–110.
. Reuter JA, Spacek DV, Pai RK, Snyder MP. Simul-seq: combined DNA and RNA sequencing for whole-genome and transcriptome profiling. Nat Methods 2016;13:953–958.
. Chen J, Suo S, Tam PP, Han JJ, Peng G, Jing N. Spatial transcriptomic analysis of cryosectioned tissue samples with Geo-seq. Nat Protoc 2017;12:566–580.
. Picelli S, Faridani OR, Björklund AK, Winberg G, Sagasser S, Sandberg R. Full-length RNA-seq from single cells using Smart-seq2. Nat Protoc 2014;9:171–181.
. Picelli S, Björklund K, Faridani OR, Sagasser S, Winberg G, Sandberg R. Smart-seq2 for sensitive full-length transcriptome profiling in single cells. Nat Methods 2013;10:1096–1098.
. Medaglia C, Giladi A, Stoler-Barak L, et al. Spatial reconstruction of immune niches by combining photoactivatable reporters and scRNA-seq. Science 2017;358:1622–1626.
. Salmen F, Ståhl PL, Mollbrink A, et al. Barcoded solid-phase RNA capture for Spatial Transcriptomics profiling in mammalian tissue sections. Nat Protoc 2018;13:2501–2534.
. Liu Y, Yang M, Deng Y, et al. High-spatial-resolution multi-omics sequencing via deterministic barcoding in tissue. Cell 2020;183:1665–1681.e18.
. Vickovic S, Eraslan G, Salmén F, et al. High-definition spatial transcriptomics for in situ tissue profiling. Nat Methods 2019;16:987–990.
. Chen A, Liao S, Cheng M, et al. Spatiotemporal transcriptomic atlas of mouse organogenesis using DNA nanoball patterned arrays. bioRxiv 2021. doi: https://doi.org/10.1101/2021.01.17.427004
. Cho CS, Xi J, Si Y, et al. Microscopic examination of spatial transcriptome using Seq-Scope. Cell 2021;184:3559–3572.e22.
. Larsson C, Koch J, Nygren A, et al. In situ genotyping individual DNA molecules by target-primed rolling-circle amplification of padlock probes. Nat Methods 2004;1:227–232.
. Gyllborg D, Langseth CM, Qian X, et al. Hybridization-based in situ sequencing (HybISS) for spatially resolved transcriptomics
in human and mouse brain tissue. Nucleic Acids Res 2020;48:e112.
. Lee JH, Daugharthy ER, Scheiman J, et al. Highly multiplexed subcellular RNA sequencing in situ. Science 2014;343:1360–1363.
. Lee JH, Daugharthy ER, Scheiman J, et al. Fluorescent in situ sequencing (FISSEQ) of RNA for gene expression profiling in intact cells and tissues. Nat Protoc 2015;10:442–458.
. Wang X, Allen WE, Wright MA, et al. Three-dimensional intact-tissue sequencing of single-cell transcriptional states. Science 2018;361:eaat5691.
. Foley JW, Zhu C, Jolivet P, et al. Gene expression profiling of single cells from archival tissue with laser-capture microdissection and Smart-3SEQ. Genome Res 2019;29:1816–1825.
. Meylan M, Petitprez F, Becht E, et al. Tertiary lymphoid structures generate and propagate anti-tumor antibody-producing plasma cells in renal cell cancer. Immunity 2022;55:527–541.e5.
. Zeng Z, Li Y, Li Y, Luo Y. Statistical and machine learning methods for spatially resolved transcriptomics
data analysis. Genome Biol 2022;23:83.
. Edsgard D, Johnsson P, Sandberg R. Identification of spatial expression trends in single-cell gene expression data. Nat Methods 2018;15:339–342.
. Svensson V, Teichmann SA, Stegle O. SpatialDE: identification of spatially variable genes. Nat Methods 2018;15:343–346.
. Sun S, Zhu J, Zhou X. Statistical analysis of spatial expression patterns for spatially resolved transcriptomic studies. Nat Methods 2020;17:193–200.
. Zhu J, Sun S, Zhou X. SPARK-X: non-parametric modeling enables scalable and robust detection of spatial expression patterns for large spatial transcriptomic studies. Genome Biol 2021;22:184.
. Anderson A, Lundeberg J. sepal: Identifying transcript profiles with spatial patterns by diffusion-based modeling. Bioinformatics 2021;37:2644–2650.
. Hu J, Li X, Coleman K, et al. SpaGCN: integrating gene expression, spatial location and histology to identify spatial domains and spatially variable genes by graph convolutional network. Nat Methods 2021;18:1342–1351.
. Moehlin J, Mollet B, Colombo BM, Mendoza-Parra MA. Inferring biologically relevant molecular tissue substructures by agglomerative clustering of digitized spatial transcriptomes with multilayer. Cell Syst 2021;12:694–705.e3.
. Zhao E, Stone MR, Ren X, et al. Spatial transcriptomics at subspot resolution with BayesSpace. Nat Biotechnol 2021;39:1375–1384.
. Yang Y, Shi X, Liu W, et al. SC-MEB: spatial clustering with hidden Markov random field using empirical Bayes. Brief Bioinform 2022;23:bbab466.
. Dong R, Yuan GC. SpatialDWLS: accurate deconvolution of spatial transcriptomic data. Genome Biol 2021;22:145.
. Elosua-Bayes M, Nieto P, Mereu E, Gut I, Heyn H. SPOTlight: seeded NMF regression to deconvolute spatial transcriptomics spots with single-cell transcriptomes. Nucleic Acids Res 2021;49:e50.
. Cable DM, Murray E, Zou LS, et al. Robust decomposition of cell type mixtures in spatial transcriptomics. Nat Biotechnol 2022;40:517–526.
. Biancalani T, Scalia G, Buffoni L, et al. Deep learning and alignment of spatially resolved single-cell transcriptomes with Tangram. Nat Methods 2021;18:1352–1362.
. Korsunsky I, Millard N, Fan J, et al. Fast, sensitive and accurate integration of single-cell data with Harmony. Nat Methods 2019;16:1289–1296.
. Welch JD, Kozareva V, Ferreira A, Vanderburg C, Martin C, Macosko EZ. Single-cell multi-omic integration compares and contrasts features of brain cell identity. Cell 2019;177:1873–1887.e17.
. Satija R, Farrell JA, Gennert D, Schier AF, Regev A. Spatial reconstruction of single-cell gene expression data. Nat Biotechnol 2015;33:495–502.
. Abdelaal T, Mourragui S, Mahfouz A, Reinders MJ. T. SpaGE: spatial gene enhancement using scRNA-seq. Nucleic Acids Res 2020;48:e107.
. Shengquan C, Boheng Z, Xiaoyang C, Xuegong Z, Rui J. stPlus: a reference-based method for the accurate enhancement of spatial transcriptomics. Bioinformatics 2021;37:i299–i307.
. Maseda F, Cang Z, Nie Q. DEEPsc: a deep learning-based map connecting single-cell transcriptomics and spatial imaging data. Front Genet 2021;12:636743.
. Cang Z, Nie Q. Inferring spatial and signaling relationships between cells from single cell transcriptomic data. Nat Commun 2020;11:2084.
. Nitzan M, Karaiskos N, Friedman N, Rajewsky N. Gene expression cartography. Nature 2019;576:132–137.
. Ren X, Zhong G, Zhang Q, Zhang L, Sun Y, Zhang Z. Reconstruction of cell spatial organization from single-cell RNA sequencing data based on ligand-receptor mediated self-assembly. Cell Res 2020;30:763–778.
. Arnol D, Schapiro D, Bodenmiller B, Saez-Rodriguez J, Stegle O. Modeling cell-cell interactions from spatial molecular data with spatial variance component analysis. Cell Rep 2019;29:202–211.e6.
. Kleshchevnikov V, Shmatko A, Dann E, et al. Cell2location maps fine-grained cell types in spatial transcriptomics. Nat Biotechnol 2022;40:661–671.
. Li B, Zhang W, Guo C, et al. Benchmarking spatial and single-cell transcriptomics integration methods for transcript distribution prediction and cell type deconvolution. Nat Methods 2022;19:662–670.
. Zeng H, Sanes JR. Neuronal cell-type classification: challenges, opportunities and the path forward. Nat Rev Neurosci 2017;18:530–546.
. Codeluppi S, Borm LE, Zeisel A, et al. Spatial organization of the somatosensory cortex revealed by osmFISH. Nat Methods 2018;15:932–935.
. Moffitt JR, Bambah-Mukku D, Eichhorn SW, et al. Molecular, spatial, and functional single-cell profiling of the hypothalamic preoptic region. Science 2018;362:eaau5324.
. Maniatis S, ijö T, Vickovic S., Braine C. Spatiotemporal dynamics of molecular pathology in amyotrophic lateral sclerosis. Science 2019;364:89–93.
. Chen WT, Lu A, Craessaerts K, et al. Spatial transcriptomics and in situ sequencing to study Alzheimer’s disease. Cell 2020;182:976–991.e19.
. Janiszewska M. The microcosmos of intratumor heterogeneity: the space-time of cancer evolution. Oncogene 2020;39:2031–2039.
. Gerlinger M, Rowan AJ, Horswell S, et al. Intratumor heterogeneity and branched evolution revealed by multiregion sequencing. N Engl J Med 2012;366:883–892.
. Hinohara K, Polyak K. Intratumoral heterogeneity: more than just mutations. Trends Cell Biol 2019;29:569–579.
. Berglund E, Maaskola J, Schultz N, et al. Spatial maps of prostate cancer transcriptomes reveal an unexplored landscape of heterogeneity. Nat Commun 2018;9:2419.
. Thrane K, Eriksson H, Maaskola J, Hansson J, Lundeberg J. Spatially resolved transcriptomics
enables dissection of genetic heterogeneity in stage III cutaneous malignant melanoma. Cancer Res 2018;78:5970–5979.
. Wu R, Guo W, Qiu X, et al. Comprehensive analysis of spatial architecture in primary liver cancer. Sci Adv 2021;7:eabg3750.
. Moncada R, Barkley D, Wagner F, et al. Integrating microarray-based spatial transcriptomics and single-cell RNA-seq reveals tissue architecture in pancreatic ductal adenocarcinomas. Nat Biotechnol 2020;38:333–342.
. Ji AL, Rubin AJ, Thrane K, et al. Multimodal analysis of composition and spatial architecture in human squamous cell carcinoma. Cell 2020;182:497–514.e22.
. Andersson A, Larsson L, Stenbeck L, et al. Spatial deconvolution of HER2-positive breast tumors reveals novel intercellular relationships. bioRxiv 2020. doi:10.1101/2020.07.14.200600.
. Svedlund J, Strell C, Qian X, et al. Generation of in situ sequencing based OncoMaps to spatially resolve gene expression profiles of diagnostic and prognostic markers in breast cancer. EBioMedicine 2019;48:212–223.
. He B, Bergenstråhle L, Stenbeck L, et al. Integrating spatial gene expression and breast tumour morphology via deep learning. Nat Biomed Eng 2020;4:827–834.
. Zugazagoitia J, Gupta S, Liu Y, et al. Biomarkers associated with beneficial PD-1 checkpoint blockade in non-small cell lung cancer (NSCLC) identified using high-plex digital spatial profiling. Clin Cancer Res 2020;26:4360–4368.
. Gouin KH 3rd, Ing N, Plummer JT, et al. An N-Cadherin 2 expressing epithelial cell subpopulation predicts response to surgery, chemotherapy and immunotherapy in bladder cancer. Nat Commun 2021;12:4906.
. Peng G, Cui G, Ke J, Jing N. Using single-cell and spatial transcriptomes to understand stem cell lineage specification during early embryo development. Annu Rev Genomics Hum Genet 2020;21:163–181.
. Stapel LC, Lombardot B, Broaddus C, et al. Automated detection and quantification of single RNAs at cellular resolution in zebrafish embryos. Development 2016;143:540–546.
. VanHorn S, Morris SA. Next-generation lineage tracing and fate mapping to interrogate development. Dev Cell 2021;56:7–21.
. Beccari L, Moris N, Girgin M, et al. Multi-axial self-organization properties of mouse embryonic stem cells into gastruloids. Nature 2018;562:272–276.
. van den Brink SC, Baillie-Johnson P, Balayo T, et al. Symmetry breaking, germ layer specification and axial organisation in aggregates of mouse embryonic stem cells. Development 2014;141:4231–4242.
. Peng G, Suo S, Cui G, et al. Molecular architecture of lineage allocation and tissue organization in early mouse embryo. Nature 2019;572:528–532.
. van den Brink SC, Alemany A, van Batenburg V, et al. Single-cell and spatial transcriptomics reveal somitogenesis in gastruloids. Nature 2020;582:405–409.
. Moris N, Anlas K, van den Brink SC, et al. An in vitro model of early anteroposterior organization during human development. Nature 2020;582:410–415.
. Liu C, Li R, Li Y, et al. Spatiotemporal mapping of gene expression landscapes and developmental trajectories during zebrafish embryogenesis. Dev Cell 2022;57:1284–1298.e5.
. Chen A, Liao S, Cheng M, et al. Spatiotemporal transcriptomic atlas of mouse organogenesis using DNA nanoball-patterned arrays. Cell 2022;185:1777–1792.e21.
. Asp M, Giacomello S, Larsson L, et al. A spatiotemporal organ-wide gene expression and cell atlas of the developing human heart. Cell 2019;179:1647–1660.e19.
. Fawkner-Corbett D, Antanaviciute A, Parikh K, et al. Spatiotemporal analysis of human intestinal development at single-cell resolution. Cell 2021;184:810–826.e23.
. Hou X, Yang Y, Li P, et al. Integrating spatial transcriptomics and single-cell RNA-seq reveals the gene expression profling of the human embryonic liver. Front Cell Dev Biol 2021;9:652408.
. Crane GM, Jeffery E, Morrison SJ. Adult haematopoietic stem cell niches. Nat Rev Immunol 2017;17:573–590.
. Golub R, Cumano A. Embryonic hematopoiesis. Blood Cells Mol Dis 2013;51:226–231.
. Sturgeon CM, Ditadi A, Awong G, Kennedy M, Keller G. Wnt signaling controls the specification of definitive and primitive hematopoiesis from human pluripotent stem cells. Nat Biotechnol 2014;32:554–561.
. Xue Y, Liu D, Cui G, et al. A 3D atlas of hematopoietic stem and progenitor cell expansion by multi-dimensional RNA-Seq analysis. Cell Rep 2019;27:1567–1578.e5.
. Crosse EI, Gordon-Keylock S, Rybtsov S, et al. Multi-layered spatial transcriptomics identify secretory factors promoting human hematopoietic stem cell development. Cell Stem Cell 2020;27:822–839.e8.
. Gao S, Shi Q, Zhang Y, et al. Identification of HSC/MPP expansion units in fetal liver by single-cell spatiotemporal transcriptomics. Cell Res 2022;32:38–53.
. Calvanese V, Capellera-Garcia S, Ma F, et al. Mapping human haematopoietic stem cells from haemogenic endothelium to birth. Nature 2022;604:534–540.
. Zhang P, Li X, Pan C, et al. Single-cell RNA sequencing to track novel perspectives in HSC heterogeneity. Stem Cell Res Ther 2022;13:39.
. Tikhonova AN, Dolgalev I, Hu H, et al. The bone marrow microenvironment at single-cell resolution. Nature 2019;569:222–228.
. Baryawno N, Przybylski D, Kowalczyk MS, et al. A cellular taxonomy of the bone marrow stroma in homeostasis and leukemia. Cell 2019;177:1915–1932.e16.
. Bigildeev AE, Petinati NA, Drize NJ. How methods of molecular biology shape our understanding of the hematopoietic system. Mol Biol 2019;53:626–637.
. Travlos GS. Normal structure, function, and histology of the bone marrow. Toxicol Pathol 2006;34:548–565.
. Baccin C, Al-Sabah J, Velten L, et al. Combined single-cell and spatial transcriptomics reveal the molecular, cellular and spatial bone marrow niche organization. Nat Cell Biol 2020;22:38–48.
. Shi L, Liu Y-H, Li K, et al. An AIE-based probe for rapid and ultrasensitive imaging of plasma membranes in biosystems. Angew Chem Int Ed Engl 2020;59:9962–9966.
. Wang KC, Chang HY. Epigenomics: technologies and applications. Circ Res 2018;122:1191–1199.
. Deng Y, Bartosovic M, Kukanja P, et al. Spatial-CUT&Tag: spatially resolved chromatin modification profiling at the cellular level. Science 2022;375:681–686.
. Lee AC, Lee Y, Choi A, et al. Spatial epitranscriptomics reveals A-to-I editome specific to cancer stem cell microniches. Nat Commun 2022;13:2540.
. Johnson CH, Ivanisevic J, Siuzdak G. Metabolomics: beyond biomarkers and towards mechanisms. Nat Rev Mol Cell Biol 2016;17:451–459.
. He J, Sun C, Li T, et al. A sensitive and wide coverage ambient mass spectrometry imaging method for functional metabolites based molecular histology. Adv Sci (Weinh) 2018;5:1800250.
. Kompauer M, Heiles S, Spengler B. Atmospheric pressure MALDI mass spectrometry imaging of tissues and cells at 1.4-mum lateral resolution. Nat Methods 2017;14:90–96.
. Yuan Z, Zhou Q, Cai L, et al. SEAM is a spatial single nuclear metabolomics method for dissecting tissue microenvironment. Nat Methods 2021;18:1223–1232.
. Sun C, Li T, Song X, et al. Spatially resolved metabolomics to discover tumor-associated metabolic alterations. Proc Natl Acad Sci USA 2019;116:52–57.