Spatial transcriptomics explained

Spatial transcriptomics is a method for assigning cell types (identified by the mRNA readouts) to their locations in the histological sections. Recent work demonstrated that the subcellular localization of mRNA molecules, for example, in the nucleus can also be studied.[1] The first widely-adopted method, described by Ståhl et al., in a landmark 2016 paper in Science[2] implies positioning individual tissue samples on the glass slides or arrays including "spots" of spatially barcoded tails to capture the poly-adenylated mRNAs. Besides tail and spatial barcode, which indicates the x and y position on the arrayed slide, the probe contains an amplification and sequencing handle, and unique molecular identifier. Commonly, histological samples are cut using cryotome, then fixed, stained, and put on the microarrays. After that, it undergoes enzymatic permeabilization, so that RNA molecules can diffuse down to the slide, with further hybridization of polyadenilated mRNA molecules to the oligo (dT) probes. Reverse transcription is then carried out in situ. As a result, spatially marked complementary DNA (cDNA) is synthesized, providing information about gene expression in the exact location of the sample. From the cDNA, libraries are synthesised for short-read sequencing. In summary, the first widely-adopted spatial transcriptomics protocol combines paralleled sequencing and staining of the same sample. In the downstream analysis, bioinformatic tools allow to overlapp the tissue image with the gene expression. The output is a map of the transcriptome capture the gene expression of individual cells within a tissue section. It is important to mention that the first generation of the arrayed slides comprised about 1,000 spots of the 100-μm diameter, limiting resolution to ~10-40 cells per spot. Spatial transcriptomics includes methods that can be divided into two modalities, those based in next-generation sequencing for gene detection, and those based in imaging.[3] Some common approached to resolve spatial distribution of transcripts are microdisection techniques, Fluorescent in situ hybridization methods, in situ sequencing, in situ capture protocols and in silico approaches.[4]

Applications

Defining the spatial distribution of mRNA molecules allows for the detection of cellular heterogeneity in tissues, tumours, immune cells as well as determine the subcellular distribution of transcripts in various conditions. This information provides a unique opportunity to decipher both the cellular and subcellular architecture in both tissues and individual cells. These methodologies provide crucial insights in the fields of embryology, oncology, immunology and histology. The functioning of the individual cells in multicellular organisms can only be completely explained in the context of identifying their exact location in the body. Spatial transcriptomics techniques sought to elucidate cells’ properties this way. Below, we look into the methods that connect gene expression to the spatial organization of cells.

Microdissection

Laser capture microdissection

Laser capture microdissection enables capturing single cells without causing morphologic alterations.[5] [6] It exploits transparent ethylene vinyl acetate film apposed to the histological section and a low-power infrared laser beam. Once such beam is directed at the cells of interest, film directly above the targeted area temporarily melts so that its long-chain polymers cover and tightly capture the cells. Then, the section is removed and cells of interest remain embedded in the film. This method allows further RNA transcript profiling and cDNA library generation of the retrieved cells.

RNA sequencing of individual cryosections

RNA sequencing of the selected regions in individual cryosections is another method that can produce location-based genome-wide expression data.[7] This method is carried out without laser capture microdissection. It was first used to determine genome-wide spatial patterns of gene expression in cryo-sliced Drosophila embryos. Essentially, it implies simple preparation of the library from the selected regions of the sample. This method had difficulties in obtaining high-quality RNA-seq libraries from every section due to the material loss as a result of the small amount of total RNA in each slice. This problem was resolved by adding RNA of a distantly related Drosophila species to each tube after initial RNA extraction.

TIVA

Transcriptome in vivo analysis (TIVA) is a technique that enables capturing mRNA in live single cells in intact live tissue sections.[8] It uses a photoactivatable tag. The TIVA tag has several functional groups and a trapped poly(U) oligonucleotide coupled to biotin. A disulfide-linked peptide, which is adjacent to the tag, allows it to penetrate the cell membrane. Once inside, laser photoactivation is used to unblock poly(U) oligonucleotide in the cells of interest, so that TIVA tag hybridizes to mRNAs within the cell. Then, streptavidin capture of the biotin group is used to extract poly(A)-tailed mRNA molecules bound to unblocked tags, after which these mRNAs are analyzed by RNA sequencing. This method is limited by low throughput, as only a few single cells can be processed at a time.

tomo-seq

An advanced alternative for RNA Sequencing of Individual Cryosections described above, RNA tomography (tomo-seq) features better RNA quantification and spatial resolution.[9] It is also based on tissue cryosectioning with further RNA sequencing of individual sections, yielding genome-wide expression data and preserving spatial information. In this protocol, usage of carrier RNA is omitted due to linear amplification of cDNA in individual histological sections. The identical sample is sectioned in different directions followed by 3D transcriptional construction using overlapping data. Overall, this method implies using identical samples for each section and thus cannot be applied for processing clinical material.

LCM-seq

LCM-seq utilizes laser capture microdissection (LCM) coupled with Smart-Seq2 RNA sequencing and is applicable down to the single cell level and can even be used on partially degraded tissues. The workflow includes cryosectioning of tissues followed by laser capture microdissection, where cells are collected directly into lysis buffer and cDNA is generated without the need for RNA isolation, which both simplifies the experimental procedures as well as lowers technical noise.[10] [11] As the positional identity of each cell is recorded during the LCM procedure, the transcriptome of each cell after RNA sequencing of the corresponding cDNA library can be inferred to the position where it was isolated from.[10] LCM-seq has been applied to multiple cell types to understand their intrinsic properties, including oculomotor neurons, facial motor neurons, hypoglossal motor neurons, spinal motor neurons, red nucleus neurons,[12] interneurons,[13] dopamine neurons,[14] and chondrocytes.[15]

Geo-seq

Geo-seq is a method that utilizes both laser capture microdissection and single-cell RNA sequencing procedures to determine the spatial distribution of the transcriptome in tissue areas approximately ten cells in size.[16] The workflow involves removal and cryosectioning of tissue followed by laser capture microdissection.[17] The extracted tissue is then lysed, and the RNA is purified and reverse transcribed into a cDNA library. Library is sequenced, and the transcriptomic profile can be mapped to the original location of the extracted tissue. This technique allows the user to define regions of interest in a tissue, extract said tissue and map the transcriptome in a targeted approach.

NICHE-seq

The NICHE-seq method uses photoactivatable fluorescent markers and two-photon laser scanning microscopy to provide spatial data to the transcriptome generated.[18] The cells bound by the fluorescent marker are photoactivated, dissociated and sorted via fluorescence-activated cell sorting. This provides sorting specificity to only labeled, photoactivated cells. Following sorting, single-cell RNA sequencing generates the transcriptome of the visualized cells. This method can process thousands of cells within a defined niche at the cost of losing spatial data between cells in the niche.

ProximID

ProximID is a methodology based on iterative micro digestion of extracted tissue to single cells.[19] Initial mild digestion steps preserve small interacting structures that are recorded prior to continued digestion. The single cells are then separated from each structure and undergo sc-RNAseq and clustered using t-distributed stochastic neighbour embedding.[20] The clustered cells can be mapped to physical interactions based on the interacting structures prior to the micro digestions. While the throughput of this technique is relatively low it provides information on physical interaction between cells to the dataset.

Fluorescent in situ hybridization

NanoString CosMx

CosMx Spatial Molecular Imager is the first high-plex in situ analysis platform to provide spatial multiomics with formalin-fixed paraffin-embedded (FFPE) and fresh frozen (FF) tissue samples at cellular and subcellular resolution. It enables rapid quantification and visualization of up to 1,000 RNA and 64 validated protein analytes and is the flexible, spatial single-cell imaging platform for cell atlasing, tissue phenotyping, cell-cell interactions, cellular processes, and biomarker discovery.

smFISH

One of the first techniques able to achieve spatially resolved RNA profiling of individual cells was single-molecule fluorescent in situ hybridization (smFISH).[21] It implemented short (50 base pairs) oligonucleotide probes conjugated with 5 fluorophores which could bind to a specific transcript yielding bright spots in the sample. Detection of these spots provides quantitative information about expression of certain genes in the cell. However, usage of probes labeled with multiple fluorophores was challenged by self-quenching, altered hybridization characteristics, their synthesis and purification.

Later, this method was changed[22] by substituting the above described probes with those of 20 bp length, coupled to only one fluorophore and complementary in tandem to an mRNA sequence of interest, meaning that those would collectively bind to the targeted mRNA. One such probe itself wouldn't produce a strong signal, but the cumulative fluorescence of the congregated probes would show a bright spot. Since single misbound probes are unlikely to co-localize, false-positive signals in this method are limited.

Thus, this in situ hybridization (ISH) technique spots spatial localization of RNA expression via direct imaging of individual RNA molecules in single cells.

RNAscope

Another in situ hybridization technique termed RNAscope employs probes of the specific Z-shaped design to simultaneously amplify hybridization signals and suppress background noise.[23] It allows for the visualization of single RNA in a variety of cellular types.[24] Most steps of RNAScope are similar to the classic ISH protocol. The tissue sample is fixed onto slides and then treated with RNAscope reagents that permeate the cells. Z-probes are designed in a way that they are only effective when bound in pairs to the target sequence. This allows another element of this method (preamplifier) to connect to the top tails of Z-probes. Once affixed, preamplifier serves as a binding site for other elements: amplifiers which in turn bind to another type of probes: label probes. As a result, a bulky structure is formed on the target sequence. Most importantly, the preamplifier fails to bind to a singular Z-probe, thus, nonspecific binding wouldn't entail signal emission, thus, eliminating background noise mentioned in the beginning.

seqFISH

Sequential fluorescence in situ hybridization (seqFISH) is another method that provides identification of mRNA directly in single cells with preservation of their spatial context.[25] [26] This method is carried out in multiple rounds; each of them includes fluorescent probe hybridization, imaging, and consecutive probe stripping. Various genes are assigned different colors in every round, generating a unique temporal barcode. Thus, seqFISH distinguishes mRNAs by a sequential color code, such as red-red-green. Nevertheless, this technique has its flaws featuring autofluorescent background and high costs due to the number of probes used in each round.

MERFISH

Conventional FISH methods are limited by the small number of genes that can be simultaneously analyzed due to the small number of distinct color channels, so multiplexed error-robust FISH was designed to overcome this problem.[27] Multiplexed Error-Robust FISH (MERFISH) greatly increases the number of RNA species that can be simultaneously imaged in single cells employing binary code gene labeling in multiple rounds of hybridization.[28] This approach can measure 140 RNA species at a time using an encoding scheme that both detects and corrects errors. The core principle lies in identification of genes by combining signals from several consecutive hybridization rounds and assigning N-bit binary barcodes to genes of interest. The Code depends on specific probes and comprises “1” or “0” values and their combination is set differently for each gene. Errors are avoided by using six-bit or longer codes with any two of them differing by at least 3 bits. A specific probe is created for each RNA species. Each probe is a target-specific oligonucleotide that consists of 20-30 base pairs and complementary binds to mRNA sequence after permeating the cell. Then, multiple rounds of hybridization are conducted as follows: for each round, only a probe that includes “1” in the corresponding binary code position is added. At the end of each round, fluorescent microscopy is used to locate each probe. Expectedly, only those mRNAs which had “1” in the assigned position would be captured. Photos are then photobleached and a new subset is added. Thus, we retrieve combination of binary values which makes it possible to distinguish between numerous RNA species.

smHCR

Single-molecule RNA detection at depth by hybridization chain reaction (smHCR) is an advanced seqFISH technique that can overcome typical complication of autofluorescent background in thick and opaque tissue samples.[29] In this method, multiple readout probes are bound with the target region of mRNA. Target is detected by a set of short DNA probes which attach to it in defined subsequence. Each DNA probe carries an initiator for the same HCR amplifier. Then, fluorophore-labeled DNA HCR hairpins penetrate the sample and assemble into fluorescent amplification polymers attaching to initiating probes. In multiplexed studies, the same two-stage protocol described above is used: all probe sets are introduced simultaneously, just as all HCR amplifiers are; spectrally distinct fluorophores are used for further imaging.

osmFISH

Cyclic-ouroboros smFISH (osmFISH) is an adaptation of smFISH which aims to overcome the challenge of optical crowding.[30] In osmFISH, transcripts are visualized, and an image is acquired before the probe is stripped and a new transcript is visualized with a different fluorescent probe. After successive rounds the images are compiled to view the spatial distribution of the RNA. Due to transcripts being sequentially visualized it eliminates the issue of signals interfering with each other. This method allows the user to generate high resolution images of larger tissue sections than other related techniques.

ExFISH

Expansion FISH (ExFISH) leverages expansion microscopy to allow for super-resolution imaging of RNA location, even in thick specimens such as brain tissue.[31] It supports both single-molecule and multiplexed readouts.

EASI-FISH

Expansion-Assisted Iterative Fluorescence In Situ Hybridization (EASI-FISH) optimizes and builds on ExFISH with improved detection accuracy and robust multi-round processing across samples thicker (300 μm) than what was previously possible.[32] It also includes a turn-key computational analysis pipeline.[33]

seqFISH+

SeqFISH+ resolved optical issues related to spatial crowding by subsequent rounds of fluorescence.[34] First, a primary probe anneals to targeted mRNA and then subsequent probes bind to flanking regions of the primary probe resulting in a unique barcode. Each readout probe is captured as an image and collapsed into a super resolved image. This method allows the user to target up to ten thousand genes at a time.

DNA microscopy

DNA microscopy is a distinct imaging method for optics-free mapping of molecules’ positions with simultaneous preservation of sequencing data carried out in several consecutive in situ reactions.[35] First, cells are fixed and cDNA is synthesized. Randomized nucleotides then tag target cDNAs in situ, providing unique labels for each molecule. Tagged transcripts are amplified in the second in situ reaction, retrieved copies are concatenated, and new randomized nucleotides are added. Each consecutive concatenation event is labeled, yielding unique event identifiers. Algorithm then generates images of the original transcripts based on decoded molecular proximities from the obtained concatenated sequences, while target's single nucleotide information is being recorded as well.

in situ sequencing

ISS using padlock probes

The ISS padlock method[36] is based on padlock probing,[37] rolling-circle amplification (RCA),[38] and sequencing by ligation chemistry.[39] Within intact tissue sections, mRNA is reversely transcribed to cDNA, which is followed by mRNA degradation by RNase H. Then, there are two ways of how this method can be carried out. The first way, gap-targeted sequencing, involves padlock probe binding to cDNA with a gap between the ends of the probe which are targeted for sequencing by ligation. DNA polymerization then fills this gap and a DNA circle is created by DNA ligation. Another way, barcode-targeted sequencing, DNA circularization of a padlock probe with a barcode sequence is conducted by ligation only. In both versions of the method, the ends are ligated forming a circle of DNA. Target amplification is then performed by RCA, yielding micrometer-sized RCA products (RCPs). RCAs consist of repeats of the padlock probe sequence. These DNA molecules are then subjected to sequencing by ligation, decoding either a gap-filled sequence or an up to four-base-long barcode within the probe with adjacent ends, depending on the version. No-gap variant claims higher sensitivity, while gap-filled one implies reading out the actual RNA sequence of the transcript. Later, this method was improved by automatization on a microfluidic platform and substitution of sequencing by ligation with sequencing by hybridization technology.

FISSEQ

Fluorescent in situ sequencing (FISSEQ),[40] like ISS padlock, is a  method that uses reverse transcription, rolling-circle amplification, and sequencing by ligation techniques. It allows spatial transcriptome analysis in fixed cells. RNA is first reverse transcribed into cDNA with regular and modified amine-bases and tagged random hexamer RT primers. Amine-bases mediate the cross-linkage of cDNA to its cellular surrounding. Then cDNA is circulated by ligation and amplified by RCA. Single-stranded DNA nanoballs of 200–400 nm in diameter are obtained as a result. Thus, these nanoballs comprise numerous tandem repeats of the cDNA sequence. Then sequencing is performed via SOLiD sequencing by ligation. Positions of both product of reverse transcription and clonally amplified RCPs are maintained via cross-linkage to cellular matrix components mentioned previously, creating a 3D in situ RNA-seq library within the cell. Once bound with fluorescent probes featuring different colors, amplicons become highly fluorescent which allows visual detection of the signal; however, the image-processing algorithm relies on read alignment to reference sequences rather than signal intensity.

Barista-seq

Barcode in situ targeted sequencing (Barista-seq) is an improvement on the gap padlock probe methodology boasting a fivefold increase in efficiency, an increased read length of fifteen bases and is compatible with illumina sequencing platforms. The method also uses padlock probes and rolling circle amplification, however this approach uses sequencing-by-synthesis and crosslinking unlike the gap padlock method. The crosslinking to the cellular matrix in the same procedure is the same as FISSEQ.

STARmap

Spatially-resolved transcript amplicon readout mapping (STARmap) utilizes a padlock probe with an additional primer which allows for direct amplification of mRNA, forgoing the need for reverse transcription.[41] Similar to other padlock probe based methods amplification occurs via rolling circle amplification. The DNA amplicons are chemically modified and embedded into a polymerized hydrogel within the cell. Captured RNA can then be sequenced in situ providing three dimensional locations of the mRNA within each cell.

in situ capture

Stereo-seq(STOmics)

STOmics is a pioneer in advancing spatially-resolved transcriptomic analysis through its proprietary SpaTial Enhanced REsolution Omics-Sequencing (Stereo-seq) technology.[42] It combines in situ capture with DNB-seq, DNB sequencing is based on lithographically etched chips (patterned arrays) for in situ sequencing. Unlike other um-level in situ capture technologies, standard DNB chips have spots with approximately 220 nm diameter and a center-to-center distance of 500 nm, providing up to 20000 spots for tissue RNA capture per 10mm linear distance, or 4x108 spots per 1cm2. Therefore, STOmics can show higher resolution and wider field of view than other in situ capture technologies.[43]

Spatial transcriptomics

Spatial Transcriptomics is the ability to capture the positional context of transcriptional activity within intact tissue, either for regions or single cells.

When a tissue cryosection is attached to a spatial transcriptomic slide the barcoded primers bind and capture adjacent mRNAs from the tissue. While the tissue section is attached to the slide, reverse transcription of captured mRNA is initiated and the resulting cDNA incorporates the spatial barcode of the primer. Following mRNA capture and reverse transcription, sequencing libraries are prepared and analyzed with Illumina dye sequencing. The spatial barcode present within each generated sequence allows the data for each individual mRNA transcript to be mapped back to its point of origin within the tissue section.

NanoString GeoMx

NanoString's GeoMx Digital Spatial Profiler (DSP) allows the user to define a microscopic region of interest on an FFPE or frozen tissue slide due to a UV-photocleavable barcode engineered into the in-situ hybridization probes.[44] The region of interest is specifically exposed to UV light, the barcodes are cleaved and used to identify the RNA or protein present in the tissue. The defined regions of interest can vary in size between ten and six hundred micrometers allowing targeting of a wide variety of structures and cells in the histological sample.

Slide-seq

Slide-seq relies on the attachment of RNA binding, DNA-barcoded micro beads to a rubber coated glass coverslip. The microbeads are mapped to their spatial location via SOLiD sequencing. Tissue sections are transferred to this coverslip to capture extracted RNA. Captured RNA is amplified and sequenced. Transcript localization is determined by the barcode oligonucleotide sequence from the bead that captured it.

APEX-seq

APEX-seq allows the for assessment of the spatial transcriptome in different regions of a cell.[45] The method utilizes the APEX2 gene, expressed in live cells which are incubated with biotin-phenol and hydrogen peroxide. In these conditions the APEX2 enzymes catalyse the transfer of biotin groups to the RNA molecules and these can then be purified via streptavidin bead purification. The purified transcripts are then sequenced to determine which molecules were in close proximity to the biotin tagging enzyme.

HDST

High-Definition Spatial Transcriptomics (HDST) begins with decoding the location of mRNA capture beads in wells on a glass slide.[46] This is accomplished by sequential hybridization to the barcode oligonucleotide sequence of each bead. Once the location of each bead is decoded, a tissue sample can be placed on the slide and permeabilized. The captured transcripts are then sequenced. HDST uses smaller beads than Slide-seq and thus can resolve at a spatial resolution of two micrometers compared to ten micrometers of Slide-seq.

10X Genomics Visium

The 10X Genomics Visium assay is a newer and improved version of the Spatial Transcriptomics assay. It also utilizes spotted arrays of mRNA-capturing probes on the surface of glass slides but with increased spot number, minimized spot size and increased amount of capture probes per spot. Within each of the four capture areas of the Visium Spatial Gene Expression slides, there are approximately 5000 barcoded spots, which in turn contain millions of spatially barcoded capture oligonucleotides. Tissue mRNA is released upon permeabilization and binds to the barcoded oligos, enabling capture of gene expression information. Each barcoded spot is 55 μm in diameter, and the distance from the center of one spot to the center of another is approximately 100 μm. The spots are staggered to minimize the distance between them. On average, mRNA from anywhere between 1 and 10 cells are captured per spot which provides near single-cell resolution.[47]

in silico construction

Reconstruction using ISH

in silico Spatial Reconstruction with ISH implies computational spatial reconstruction of cells’ locations according to their expression profiles. Several similar methods of this principle exist.[48] [49] They co-analyze single-cell transcriptomics and available ISH-based gene expression atlases of the same cell type. Based on these data, cells are then assigned to their positions in the tissue. Obviously, this method is limited by the factor of availability of ISH references. Additionally, it becomes more complicated when assigning cells in complex tissues. This approach is not applicable for clinical samples due to the lack of paired references. Reported success rate for the exact allocation of cells in brain tissue was 81%.

DistMap

Mapping the transcriptome using the Distmap algorithm requires high-throughput single cell sequencing and an existing in situ hybridization atlas for the tissue of interest.[50] The Distmap algorithm generates a virtual 3D model of the tissue of interest using the transcriptomes of sequenced cells and said reference atlas. The transcriptomes can be clustered into cell types using t-distributed stochastic neighbour embedding and mapped to the 3D model using virtual in situ hybridization.[51] Essentially, this algorithm takes data generated from single cells in a dissociated tissue and is able to map individual transcripts to where the cell type exists in the tissue using virtual in situ hybridization.

History

in situ hybridization was developed in the late 60's and saw major developments in the 80's (smFISH) and 10's (RNAscope, seqFISH, MERFISH and osmFISH) with the most recent additions being seqFISH+ and DNA microscopy.[16] Microdisecction techniques were first developed in the late 90's (Laser Capture Microdissection) and the most recent addition was in 2016 with the development of ProximID. Spatial genomics/transcriptomics as a technique was invented and developed in 2000 by Michael Doyle (of Eolas), Maurice Pescitelli (of the University of Illinois at Chicago), Betsey Williams (of Harvard), and George Michaels (of George Mason University),[52] [53] [54] as part of the Visible Embryo Project[55] [56] It was later expanded upon in 2016 by Jonas Frisén, Joakim Lundeberg, Patrik Ståhl and their colleagues in Stockholm, Sweden. In 2019 at the Broad Institute, the labs of Fei Chen and Evan Macosko developed Slide-seq, which used barcoded oligos on beads.[57]

See also

Visible Embryo Project

Single cell sequencing

Single-cell transcriptomics

RNA-Seq

Fluorescence in situ hybridization

Notes and References

  1. Rodriques SG, Stickels RR, Goeva A, Martin CA, Murray E, Vanderburg CR, Welch J, Chen LM, Chen F, Macosko EZ . Slide-seq: A scalable technology for measuring genome-wide expression at high spatial resolution . Science . 363 . 6434 . 1463–1467 . March 2019 . 30923225 . 6927209 . 10.1126/science.aaw1219 .
  2. Ståhl PL, Salmén F, Vickovic S, Lundmark A, Navarro JF, Magnusson J, Giacomello S, Asp M, Westholm JO, Huss M, Mollbrink A, Linnarsson S, Codeluppi S, Borg Å, Pontén F, Costea PI, Sahlén P, Mulder J, Bergmann O, Lundeberg J, Frisén J . Visualization and analysis of gene expression in tissue sections by spatial transcriptomics . Science . 353 . 6294 . 78–82 . July 2016 . 27365449 . 10.1126/science.aaf2403 . 30942685 . 2016Sci...353...78S .
  3. Tian L, Chen F, Macosko EZ . The expanding vistas of spatial transcriptomics . Nature Biotechnology . 41 . 6 . 773–782 . June 2023 . 36192637 . 10091579 . 10.1038/s41587-022-01448-2 .
  4. Asp M, Bergenstråhle J, Lundeberg J . Spatially Resolved Transcriptomes-Next Generation Tools for Tissue Exploration . BioEssays . 42 . 10 . e1900221 . October 2020 . 32363691 . 10.1002/bies.201900221 . 218492475 . free .
  5. Emmert-Buck MR, Bonner RF, Smith PD, Chuaqui RF, Zhuang Z, Goldstein SR, Weiss RA, Liotta LA . Laser capture microdissection . Science . 274 . 5289 . 998–1001 . November 1996 . 8875945 . 10.1126/science.274.5289.998 . 220105962 . 1996Sci...274..998E .
  6. Simone NL, Bonner RF, Gillespie JW, Emmert-Buck MR, Liotta LA . Laser-capture microdissection: opening the microscopic frontier to molecular analysis . Trends in Genetics . 14 . 7 . 272–6 . July 1998 . 9676529 . 10.1016/S0168-9525(98)01489-9 .
  7. Combs PA, Eisen MB . Sequencing mRNA from cryo-sliced Drosophila embryos to determine genome-wide spatial patterns of gene expression . PLOS ONE . 8 . 8 . e71820 . 2013-08-12 . 23951250 . 3741199 . 10.1371/journal.pone.0071820 . 1302.4693 . 2013PLoSO...871820C . free .
  8. Lovatt D, Ruble BK, Lee J, Dueck H, Kim TK, Fisher S, Francis C, Spaethling JM, Wolf JA, Grady MS, Ulyanova AV, Yeldell SB, Griepenburg JC, Buckley PT, Kim J, Sul JY, Dmochowski IJ, Eberwine J . Transcriptome in vivo analysis (TIVA) of spatially defined single cells in live tissue . Nature Methods . 11 . 2 . 190–6 . February 2014 . 24412976 . 3964595 . 10.1038/nmeth.2804 .
  9. Junker JP, Noël ES, Guryev V, Peterson KA, Shah G, Huisken J, McMahon AP, Berezikov E, Bakkers J, van Oudenaarden A . Genome-wide RNA Tomography in the zebrafish embryo . Cell . 159 . 3 . 662–75 . October 2014 . 25417113 . 10.1016/j.cell.2014.09.038 . 2713635 . free .
  10. Nichterwitz S, Chen G, Aguila Benitez J, Yilmaz M, Storvall H, Cao M, Sandberg R, Deng Q, Hedlund E . Laser capture microscopy coupled with Smart-seq2 for precise spatial transcriptomic profiling . Nature Communications . 7 . 12139 . July 2016 . 27387371 . 4941116 . 10.1038/ncomms12139 . 2016NatCo...712139N .
  11. Book: Nichterwitz S, Aguila Benitez J, Hoogstraaten R, Deng Q, Hedlund E . LCM-Seq: A Method for Spatial Transcriptomic Profiling Using Laser Capture Microdissection Coupled with PolyA-Based RNA Sequencing . RNA Detection . Methods in Molecular Biology . 2018 . 1649 . 95–110 . 10.1007/978-1-4939-7213-5_6 . 29130192 . 978-1-4939-7212-8 .
  12. Nichterwitz S, Nijssen J, Storvall H, Schweingruber C, Comley LH, Allodi I, van der Lee M, Deng Q, Sandberg R, Hedlund E . LCM-seq reveals unique transcriptional adaptation mechanisms and protective pathways of resistant neurons in spinal muscular atrophy . Genome Research . 30 . 1083–1096 . August 2020 . 8 . 32820007 . 10.1101/gr.265017.120 . 7462070 .
  13. Pereira M, Birtele M, Shrigley S, Benitez JA, Hedlund E, Parmar M, Ottosson DR . Direct Reprogramming of Resident NG2 Glia into Neurons with Properties of Fast-Spiking Parvalbumin-Containing Interneurons . Stem Cell Reports . 9 . 3 . 742–751 . September 2017 . 28844658 . 5599255 . 10.1016/j.stemcr.2017.07.023 .
  14. Aguila J, Cheng S, Kee N, Cao M, Wang M, Deng Q, Hedlund E . Spatial RNA Sequencing Identifies Robust Markers of Vulnerable and Resistant Human Midbrain Dopamine Neurons and Their Expression in Parkinson's Disease . Frontiers in Molecular Neuroscience . 14 . 699562 . 2021 . 34305528 . 8297217 . 10.3389/fnmol.2021.699562 . free .
  15. Newton PT, Li L, Zhou B, Schweingruber C, Hovorakova M, Xie M, Sun X, Sandhow L, Artemov AV, Ivashkin E, Suter S, Dyachuk V, El Shahawy M, Gritli-Linde A, Bouderlique T, Petersen J, Mollbrink A, Lundeberg J, Enikolopov G, Qian H, Fried K, Kasper M, Hedlund E, Adameyko I, Sävendahl L, Chagin AS . A radical switch in clonality reveals a stem cell niche in the epiphyseal growth plate . Nature . 567 . 7747 . 234–238 . March 2019 . 30814736 . 10.1038/s41586-019-0989-6 . 2019Natur.567..234N. 71143703.
  16. Chen J, Suo S, Tam PP, Han JJ, Peng G, Jing N . Spatial transcriptomic analysis of cryosectioned tissue samples with Geo-seq . Nature Protocols . 12 . 3 . 566–580 . March 2017 . 28207000 . 10.1038/nprot.2017.003 . 3879096 .
  17. Peng G, Suo S, Chen J, Chen W, Liu C, Yu F, Wang R, Chen S, Sun N, Cui G, Song L, Tam PP, Han JD, Jing N . Spatial Transcriptome for the Molecular Annotation of Lineage Fates and Cell Identity in Mid-gastrula Mouse Embryo . Developmental Cell . 36 . 6 . 681–97 . March 2016 . 27003939 . 10.1016/j.devcel.2016.02.020 . free .
  18. Medaglia C, Giladi A, Stoler-Barak L, De Giovanni M, Salame TM, Biram A, David E, Li H, Iannacone M, Shulman Z, Amit I . Spatial reconstruction of immune niches by combining photoactivatable reporters and scRNA-seq . Science . 358 . 6370 . 1622–1626 . December 2017 . 29217582 . 7234837 . 10.1126/science.aao4277 . 2017Sci...358.1622M .
  19. Boisset JC, Vivié J, Grün D, Muraro MJ, Lyubimova A, van Oudenaarden A . Mapping the physical network of cellular interactions . Nature Methods . 15 . 7 . 547–553 . July 2018 . 29786092 . 10.1038/s41592-018-0009-z . 29166537 .
  20. Chen X, Sun YC, Church GM, Lee JH, Zador AM . Efficient in situ barcode sequencing using padlock probe-based BaristaSeq . Nucleic Acids Research . 46 . 4 . e22 . February 2018 . 29190363 . 5829746 . 10.1093/nar/gkx1206 .
  21. Femino AM, Fay FS, Fogarty K, Singer RH . Visualization of single RNA transcripts in situ . Science . 280 . 5363 . 585–90 . April 1998 . 9554849 . 10.1126/science.280.5363.585 . 1998Sci...280..585F .
  22. Raj A, van den Bogaard P, Rifkin SA, van Oudenaarden A, Tyagi S . Imaging individual mRNA molecules using multiple singly labeled probes . Nature Methods . 5 . 10 . 877–9 . October 2008 . 18806792 . 3126653 . 10.1038/nmeth.1253 .
  23. Wang F, Flanagan J, Su N, Wang LC, Bui S, Nielson A, Wu X, Vo HT, Ma XJ, Luo Y . RNAscope: a novel in situ RNA analysis platform for formalin-fixed, paraffin-embedded tissues . The Journal of Molecular Diagnostics . 14 . 1 . 22–29 . January 2012 . 22166544 . 3338343 . 10.1016/j.jmoldx.2011.08.002 .
  24. De Biase D, Prisco F, Piegari G, Ilsami A, d'Aquino I, Baldassarre V, Zito Marino F, Franco R, Papparella S, Paciello O . RNAScope in situ Hybridization as a Novel Technique for the Assessment of c-KIT mRNA Expression in Canine Mast Cell Tumor . Frontiers in Veterinary Science . 8 . 591961 . 16 February 2021 . 33665215 . 7921150 . 10.3389/fvets.2021.591961 . free .
  25. Lubeck E, Coskun AF, Zhiyentayev T, Ahmad M, Cai L . Single-cell in situ RNA profiling by sequential hybridization . Nature Methods . 11 . 4 . 360–1 . April 2014 . 24681720 . 4085791 . 10.1038/nmeth.2892 .
  26. 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 . 92 . 2 . 342–357 . October 2016 . 27764670 . 5087994 . 10.1016/j.neuron.2016.10.001 .
  27. Chen KH, Boettiger AN, Moffitt JR, Wang S, Zhuang X . RNA imaging. Spatially resolved, highly multiplexed RNA profiling in single cells . Science . 348 . 6233 . aaa6090 . April 2015 . 25858977 . 4662681 . 10.1126/science.aaa6090 .
  28. Zhuang X . Spatially resolved single-cell genomics and transcriptomics by imaging . Nature Methods . 18 . 1 . 18–22 . January 2021 . 33408406 . 10.1038/s41592-020-01037-8 . 9805800 . 230796841 .
  29. Shah S, Lubeck E, Schwarzkopf M, He TF, Greenbaum A, Sohn CH, Lignell A, Choi HM, Gradinaru V, Pierce NA, Cai L . Single-molecule RNA detection at depth by hybridization chain reaction and tissue hydrogel embedding and clearing . Development . 143 . 15 . 2862–7 . August 2016 . 27342713 . 5004914 . 10.1242/dev.138560 .
  30. Codeluppi S, Borm LE, Zeisel A, La Manno G, van Lunteren JA, Svensson CI, Linnarsson S . Spatial organization of the somatosensory cortex revealed by osmFISH . Nature Methods . 15 . 11 . 932–935 . November 2018 . 30377364 . 10.1038/s41592-018-0175-z . 53114385 .
  31. Chen F, Wassie AT, Cote AJ, Sinha A, Alon S, Asano S, Daugharthy ER, Chang JB, Marblestone A, Church GM, Raj A, Boyden ES . Nanoscale imaging of RNA with expansion microscopy . Nature Methods . 13 . 8 . 679–684 . August 2016 . 27376770 . 4965288 . 10.1038/nmeth.3899 .
  32. Wang Y, Eddison M, Fleishman G, Weigert M, Xu S, Wang T, Rokicki K, Goina C, Henry FE, Lemire AL, Schmidt U, Yang H, Svoboda K, Myers EW, Saalfeld S, Korff W, Sternson SM, Tillberg PW . EASI-FISH for thick tissue defines lateral hypothalamus spatio-molecular organization . English . Cell . 184 . 26 . 6361–6377.e24 . December 2021 . 34875226 . 10.1016/j.cell.2021.11.024 . 244906189 . free .
  33. Wang Y, Eddison M, Fleishman G, Weigert M, Xu S, Wang T, Rokicki K, Goina C, Henry FE, Lemire AL, Schmidt U, Yang H, Svoboda K, Myers EW, Saalfeld S, Korff W, Sternson SM, Tillberg PW . EASI-FISH for thick tissue defines lateral hypothalamus spatio-molecular organization . Cell . 184 . 26 . 6361–6377.e24 . December 2021 . 34875226 . 10.1016/j.cell.2021.11.024 . 244906189 . free .
  34. Eng CL, Lawson M, Zhu Q, Dries R, Koulena N, Takei Y, Yun J, Cronin C, Karp C, Yuan GC, Cai L . Transcriptome-scale super-resolved imaging in tissues by RNA seqFISH . Nature . 568 . 7751 . 235–239 . April 2019 . 30911168 . 6544023 . 10.1038/s41586-019-1049-y . 2019Natur.568..235E .
  35. Weinstein JA, Regev A, Zhang F . DNA Microscopy: Optics-free Spatio-genetic Imaging by a Stand-Alone Chemical Reaction . Cell . 178 . 1 . 229–241.e16 . June 2019 . 31230717 . 6697087 . 10.1016/j.cell.2019.05.019 .
  36. Ke R, Mignardi M, Pacureanu A, Svedlund J, Botling J, Wählby C, Nilsson M . In situ sequencing for RNA analysis in preserved tissue and cells . Nature Methods . 10 . 9 . 857–60 . September 2013 . 23852452 . 10.1038/nmeth.2563 . 205421785 .
  37. Szemes M, Bonants P, de Weerdt M, Baner J, Landegren U, Schoen CD . Diagnostic application of padlock probes--multiplex detection of plant pathogens using universal microarrays . Nucleic Acids Research . 33 . 8 . e70 . April 2005 . 15860767 . 1087788 . 10.1093/nar/gni069 .
  38. Ali MM, Li F, Zhang Z, Zhang K, Kang DK, Ankrum JA, Le XC, Zhao W . Rolling circle amplification: a versatile tool for chemical biology, materials science and medicine . Chemical Society Reviews . 43 . 10 . 3324–41 . May 2014 . 24643375 . 10.1039/c3cs60439j . 24839827 .
  39. Churko JM, Mantalas GL, Snyder MP, Wu JC . Overview of high throughput sequencing technologies to elucidate molecular pathways in cardiovascular diseases . Circulation Research . 112 . 12 . 1613–23 . June 2013 . 23743227 . 3831009 . 10.1161/CIRCRESAHA.113.300939 .
  40. Lee JH, Daugharthy ER, Scheiman J, Kalhor R, Yang JL, Ferrante TC, Terry R, Jeanty SS, Li C, Amamoto R, Peters DT, Turczyk BM, Marblestone AH, Inverso SA, Bernard A, Mali P, Rios X, Aach J, Church GM . Highly multiplexed subcellular RNA sequencing in situ . Science . 343 . 6177 . 1360–3 . March 2014 . 24578530 . 4140943 . 10.1126/science.1250212 . 2014Sci...343.1360L .
  41. Wang X, Allen WE, Wright MA, Sylwestrak EL, Samusik N, Vesuna S, Evans K, Liu C, Ramakrishnan C, Liu J, Nolan GP, Bava FA, Deisseroth K . Three-dimensional intact-tissue sequencing of single-cell transcriptional states . Science . 361 . 6400 . eaat5691 . July 2018 . 29930089 . 6339868 . 10.1126/science.aat5691 .
  42. Web site: STOmics Stereo-seq Spatial Transcriptomics Technology .
  43. Chen A, Liao S, Cheng M, Ma K, Wu L, Lai Y, Qiu X, Yang J, Xu J, Hao S, Wang X, Lu H, Chen X, Liu X, Huang X, Li Z, Hong Y, Jiang Y, Peng J, Liu S, Shen M, Liu C, Li Q, Yuan Y, Wei X, Zheng H, Feng W, Wang Z, Liu Y, Wang Z, Yang Y, Xiang H, Han L, Qin B, Guo P, Lai G, Muñoz-Cánoves P, Maxwell PH, Thiery JP, Wu QF, Zhao F, Chen B, Li M, Dai X, Wang S, Kuang H, Hui J, Wang L, Fei JF, Wang O, Wei X, Lu H, Wang B, Liu S, Gu Y, Ni M, Zhang W, Mu F, Yin Y, Yang H, Lisby M, Cornall RJ, Mulder J, Uhlén M, Esteban MA, Li Y, Liu L, Xu X, Wang J . Spatiotemporal transcriptomic atlas of mouse organogenesis using DNA nanoball-patterned arrays . Cell . 185 . 10 . 1777–1792.e21 . May 2022 . 35512705 . 10.1016/j.cell.2022.04.003 . 10230/53512 . free .
  44. Book: Zollinger DR, Lingle SE, Sorg K, Beechem JM, Merritt CR . In Situ Hybridization Protocols . Nielsen BS, Jones J . GeoMx™ RNA Assay: High Multiplex, Digital, Spatial Analysis of RNA in FFPE Tissue . Methods in Molecular Biology . Clifton, N.J. . Humana . 2148 . 331–345 . 2020 . 32394392 . 10.1007/978-1-0716-0623-0_21 . 978-1-0716-0623-0 . 218599146 .
  45. Fazal FM, Han S, Parker KR, Kaewsapsak P, Xu J, Boettiger AN, Chang HY, Ting AY . Atlas of Subcellular RNA Localization Revealed by APEX-Seq . Cell . 178 . 2 . 473–490.e26 . July 2019 . 31230715 . 6786773 . 10.1016/j.cell.2019.05.027 .
  46. Vickovic S, Eraslan G, Salmén F, Klughammer J, Stenbeck L, Schapiro D, Äijö T, Bonneau R, Bergenstråhle L, Navarro JF, Gould J, Griffin GK, Borg Å, Ronaghi M, Frisén J, Lundeberg J, Regev A, Ståhl PL . High-definition spatial transcriptomics for in situ tissue profiling . Nature Methods . 16 . 10 . 987–990 . October 2019 . 31501547 . 6765407 . 10.1038/s41592-019-0548-y . 1721.1/126032 .
  47. Web site: Spatial Gene Expression . 10x Genomics .
  48. Achim K, Pettit JB, Saraiva LR, Gavriouchkina D, Larsson T, Arendt D, Marioni JC . High-throughput spatial mapping of single-cell RNA-seq data to tissue of origin . Nature Biotechnology . 33 . 5 . 503–9 . May 2015 . 25867922 . 10.1038/nbt.3209 . 19535470 . free .
  49. Satija R, Farrell JA, Gennert D, Schier AF, Regev A . Spatial reconstruction of single-cell gene expression data . Nature Biotechnology . 33 . 5 . 495–502 . May 2015 . 25867923 . 4430369 . 10.1038/nbt.3192 .
  50. Karaiskos N, Wahle P, Alles J, Boltengagen A, Ayoub S, Kipar C, Kocks C, Rajewsky N, Zinzen RP . The Drosophila embryo at single-cell transcriptome resolution . Science . 358 . 6360 . 194–199 . October 2017 . 28860209 . 10.1126/science.aan3235 . 206659563 . 2017Sci...358..194K . free .
  51. Kobak D, Berens P . The art of using t-SNE for single-cell transcriptomics . Nature Communications . 10 . 1 . 5416 . November 2019 . 31780648 . 6882829 . 10.1038/s41467-019-13056-x . 2019NatCo..10.5416K .
  52. US. 7613571. Method and system for the multidimensional morphological reconstruction of genome expression activity. 2009-11-03. Doyle MD, Pescitelli Jr MJ, Williams BS, Michaels GS.
  53. US. 7894997. 2011-02-22. Multidimensional morphological reconstruction of genome expression activity. Doyle MD, Pescitelli Jr MJ, Williams BS, Michaels GS .
  54. US. 10011864. 2018-07-03. Eolas Technologies Inc. . Multidimensional microdissection and morphological reconstruction of genomic or proteomic expression activity. Doyle MD, Pescitelli Jr MJ, Williams BS, Michaels GS .
  55. Doyle MD, Noe A, Michaels GS . Oliver WR . The Visible Embryo Project: A Platform for Spatial Genomics . Proceedings of SPIE . 28th AIPR Workshop: 3D Visualization for Data Exploration and Decision Making . 3905 . 2000 . 248 . 10.1117/12.384880 . 2000SPIE.3905..248D . 85027703 .
  56. Al-Janabi A . Decoding the mechanisms of embryo development with spatial biology . BioTechniques . 75 . 5 . 179–182 . November 2023 . 37855245 . 10.2144/btn-2023-0093 . free .
  57. Rodriques SG, Stickels RR, Goeva A, Martin CA, Murray E, Vanderburg CR, Welch J, Chen LM, Chen F, Macosko EZ . Slide-seq: A scalable technology for measuring genome-wide expression at high spatial resolution . Science . 363 . 6434 . 1463–1467 . March 2019 . 30923225 . 6927209 . 10.1126/science.aaw1219 . 2019Sci...363.1463R .