Spatial multi-omics analyses of the tumor immune microenvironment
Journal of Biomedical Science volume 29, Article number: 96 (2022)
In the past decade, single-cell technologies have revealed the heterogeneity of the tumor-immune microenvironment at the genomic, transcriptomic, and proteomic levels and have furthered our understanding of the mechanisms of tumor development. Single-cell technologies have also been used to identify potential biomarkers. However, spatial information about the tumor-immune microenvironment such as cell locations and cell–cell interactomes is lost in these approaches. Recently, spatial multi-omics technologies have been used to study transcriptomes, proteomes, and metabolomes of tumor-immune microenvironments in several types of cancer, and the data obtained from these methods has been combined with immunohistochemistry and multiparameter analysis to yield markers of cancer progression. Here, we review numerous cutting-edge spatial ‘omics techniques, their application to study of the tumor-immune microenvironment, and remaining technical challenges.
The tumor-immune microenvironment and the therapeutic challenge
The inter-related, co-existing, and competitive nature of interactions between tumor cells, the surrounding tissue, and infiltrating innate and adaptive immune cells results in a unique environment that varies by tumor type and is highly adapted to the tumor behavior. This complex ecosystem is composed of tumor cells, immune cells, stromal cells, fibroblasts, extracellular matrix, and blood vessels, and is referred to as the tumor-immune microenvironment (TIME). In the TIME, dynamic and bidirectional interactions occur between cells of various types through communication signals such as secreted molecules, proteins, and vesicles . Tumor cells, as well as immune and stromal cells, utilize specific metabolic pathways to survive in the oxygen- and nutrient-limiting environment of a tumor .
Tumor-infiltrating T cells, which recognize tumor-specific antigens and kill tumor cells, are modulated by multiple signals emitted by tumor cells and by myeloid cells . Monocytes and macrophages, specifically tumor-associated macrophages, can suppress or modulate the tumor killing by T cells . Other non-immune cells, such as endothelial cells and fibroblasts, are critical regulators of tumor progression in the TIME [5, 6]. Endothelial cells build blood vessels necessary for transfer of nutrients to the tumor. Cancer-associated fibroblasts not only lay down the matrix but also actively participate in immune modulation.
Recently developed therapeutic strategies that leverage the immune system to inhibit tumor growth have proven effective in some solid tumors. Although immune checkpoint blockers (ICBs) such as anti-CTLA-4 antibody ipilimumab, anti-PD-1 antibodies nivolumab and pembrolizumab, and anti-PD-L1 antibody atezolizumab have been evaluated in the treatment of many different cancers , ICB therapy works only on certain types of cancers. Loss of tumor-expressing MHC class I molecules, which initiate the T cell response, insufficient numbers of T cells in the tumor mass, and the presence of immunosuppressive cells or factors can all limit efficacy of ICBs. Thus, the TIME strongly influences the effectiveness of ICB therapies . To understand how to harness the immune system to inhibit tumor progression, a detailed understanding of immune cell distribution and function is imperative.
Single-cell technologies reveal the heterogeneity of the tumor-immune microenvironment
The development of single-cell detection platforms has enabled deep profiling of the heterogeneity of tumors and the immune system both within individual tumors and between patients [9, 10]. Genomic single-cell sequencing methods, which detect mutations at the DNA level, have been used to confirmed co-existence of dangerous mutations in individual cells, a resolution that cannot be achieved by bulk sequencing [11, 12]. Single-cell RNA sequencing (scRNA-seq) enables profiling of transcriptomes in cells of the TIME . These technologies have yielded high-resolution and unbiased profiling of cancerous cells, T cells, myeloid cells, and stromal cells and have revealed a vast heterogeneity of immune profiles across tumor types [14, 15].
Single-cell proteome detection platforms with multiplex capacity have also been developed. For example, mass cytometry (also known as cytometry by time-of-flight or CyTOF) methods can detect over 40 cellular markers at one time [16, 17]. CD4+ and CD8+ T cells, tumor-associated macrophages, and cells that express immune checkpoint markers, the hierarchy of hematopoietic stem cell differentiation, functional heterogeneity and signaling in T and natural killer (NK) cells have been revealed using the CyTOF platform [18,19,20,21,22,23,24].
Experiments using single-cell technologies have convincingly demonstrated that tumor masses usually contain multiple genetically defined subclones of cells with distinct sets of gene mutations and different transcript profiles. Immune cells, with vastly different transcript profiles characteristic of progenitor, active, exhausted, and suppressed cells, coexist with tumor. These data have raised critical questions: How do these immune players with vastly different functions coexist in the same environment? Are these cells equally distributed in the tumor mass, or do they segregate into clusters with distinct spatial and biological features? To answer these questions, a single-cell dataset embedding histological structure information is needed.
A new era of single-cell-level histological research
Long before the development of the spatial ‘omics, it was known that certain histological patterns within tumors are highly linked to patient prognosis. For example, peritumoral T cell and B cell infiltration in colorectal cancer patients is correlated with positive prognosis, whereas the depletion of lymphocytes in the tumor core is an indicator of poor prognosis . Further, stromal infiltration of T cells is linked to good prognosis for specific types of breast cancer [26, 27]. To better understand spatial correlates with prognosis researchers have used multiplexed immunohistochemistry (IHC), immunofluorescence (IF) and laser capture microdissection-based IHC/IF with auxiliary tools to select tumor regions for analyses of expression patterns of protein and RNA transcripts . These types of histological analyses are limited by resolution and problems related to bias in sample selection. Recent technological advancements in solid-phase sequencing and multiplex imaging have enabled multiplexed detection of transcripts, proteins, and metabolites in high-resolution images (Fig. 1). In the following sections, we will review recent developments in spatial multi-omics, research applications in tumor tissue exploration, and current challenges.
Spatial transcriptomics methods
The most commonly used methods for spatially resolved transcriptomic analyses are listed in Table 1, and a schematic of the spatial transcriptomics workflow is shown in Fig. 1A (upper panel). Depending on the technique, these experiments use either fresh-frozen (FF) or formalin-fixed, paraffin-embedded (FFPE) tissues. The main strategies for spatial transcriptomics are based on next-generation sequencing (NGS) and on fluorescence in situ hybridization (FISH) [29,30,31]. With NGS techniques, transcripts are encoded with position information prior to sequencing. For example, with the 10 × Genomics Visium platform, chips containing spatially barcoded oligo(dT) are used to capture mRNA from the tissue overlaid on the chip and then processed for sequencing yielding unbiased spatial transcriptomic data . However, due to the low numbers of transcripts captured, transcript abundances from multiple neighboring cells are aggregated for downstream analyses, so the effective spatial resolution is not at the single-cell level.
Other spatial transcriptomic techniques rely on spatially barcoded probes conjugated to beads to capture RNA from the tissue samples. For example, 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 each RNA transcript to be inferred by sequencing . Slide-seq V2 has higher RNA capture efficiency than Slide-seq; the efficiency is about 50% that of scRNA-seq . XYZseq encodes spatial information by physically partitioning tissue sections into an array of microwells containing barcoded primers for reverse transcription primers . The High-Definition Spatial Transcriptomics (HDST) method uses microwell-based fluorescence spatial indexing beads to capture transcripts, and fluorescence signals reveal transcript positions in tissue samples. The barcoded RNA is sequenced, and sequencing data is mapped to fluorescence data. Cyclic sequencing of individual transcripts can be used to increase resolution. For both Slide-seq and HDST, the size of beads and the capture efficiency of RNA transcripts determines the resolution of spatial information: The diameter of beads used in HDST is around 2 μm, and the beads used for Slide-seq are 10 μm in diameter. HDST captures more targets and provides higher resolution than Slide-seq .
Recently, spatial transcriptomic approaches with higher transcriptome capture efficiencies have been developed to enable capture of sub-µm-resolution images. In Seq-Scope, tissue is attached to an RNA-capturing array with a dense arrangement of barcoded clusters . mRNAs captured from the tissue are used as the template to generate cDNAs for NGS analysis. In Stereo-seq, circular amplified DNA nanoballs containing barcode sequences are generated and dispersed onto the etched chips with patterned array . An mRNA capture sequence is linked to the spatial barcode and used to capture the mRNA released from the tissue section overlaid on the chip. Each dot is 220 nm in diameter, and there is 500 or 715 nm between the dots; however, the lateral diffusion of mRNA during the capture step is more than 5 μm. Optimized Stereo-seq can read-out 1450 unique molecular identifiers per 10 μm diameter (a bin of 14 × 14 DNA nanoballs) providing spatial transcriptomic data for a tissue area as large as several cm2. Although detection of low-abundance copy transcripts is not currently possible, this technique has been applied to analyze mouse embryonic development and the transcriptomes of solid tumors at near single-cell resolution [38, 39].
Spatially encoded barcodes can also be allocated to individual cells through light-based printing of DNA barcodes onto the surface of cells. In ZipSeq, patterned illumination and photocaged oligonucleotides are used to serially print barcodes (the so-called zipcodes) onto live cells within tissues . This annealed zipcode terminates in a polyA sequence that streamlines subsequent cDNA library construction. For sequencing, the tissue sample is dissociated into single cells and scRNA-seq is performed. The final analysis combines mapping transcripts onto zipcoded regions.
In Deterministic Barcoding in Tissue sequencing (DBiT-seq), two perpendicular microfluidic chips with parallel channels are sequentially placed against the tissue section to introduce oligo-dT-tagged with combinations of spatial barcodes for each tissue pixel . The tissue is then digested to recover spatially barcoded cDNAs, and libraries are prepared and sequenced. Spatial transcriptome mapping of the developing eye in an E10 mouse embryo was achieved using DBiT-seq with 10-μm resolution. By combining DBiT-seq with immunofluorescence staining or scRNA-seq, a better understanding of specific cell types could be achieved.
Methods to profile the transcriptome within a region of interest (ROIs) have also been developed. With the Digital Spatial Profiler (DSP), the abundances of proteins or RNAs within an ROI are measured by counting unique barcoded oligonucleotides assigned to each target of interest. Antibodies or mRNA hybridization probes are conjugated with barcoded photocleavage oligonucleotides and used to stain the tissues. Next, photocleavage is induced by UV light, and the released oligonucleotides, which encode spatial information, are collected for NGS analysis . An optimized DSP method was recently used to detect 44 proteins and 96 transcripts simultaneously . Similarly, Pick-seq can be used to sequence RNAs in an ROI of a tissue sample . The ROI, selected based on immunofluorescence, is isolated by aspiration into a liquid-filled 40-μm bore needle. Subsequently, cells are lysed, and cDNA is prepared and sequenced.
Imaging-based strategies, which combine interactive in situ hybridization or in situ sequencing with high-resolution microscopy, can achieve subcellular spatial resolution and can potentially provide genome-wide transcriptomic information. However, these methods are usually technically demanding, require iterative workflows, sophisticated image analysis processes, and large collections of probes. In methods based on in situ sequencing, probes are used to hybridize pre-selected RNA targets to allow reverse transcription. Rolling circle amplification, sequencing-by-ligation, sequencing-by synthesis, or sequencing by hybridization can be used as read-out . Fluorescent in situ sequencing (FISSEQ) is an untargeted gene profiling method in which transcripts are reverse transcribed, and rolling circle amplification is used to generate cDNA amplicon nanoballs of 200–400 nm in diameter within the cell. The amplicons are then sequenced using Supported Oligonucleotide Ligation and Detection (SOLiD) technology to yield in situ transcriptomic data . Targeted in situ transcriptomics images have been obtained using Spatially resolved Transcript Amplicon Readout mapping (STARmap), which is a combination of hydrogel tissue chemistry methods with targeted signal amplification- And a targeted method, Sequencing with Error-reduction by Dynamic Annealing and Ligation (SEDAL) method, was used to simultaneously map 1020 genes in mouse V1 neocortex .
The method known as MERFISH, for multiplexed error-robust FISH , is an advanced method derived from single-molecule FISH (smFISH) [49, 50]. Multiplexing results from rounds of hybridization, imaging, and stripping. In smFISH due to numerous single image merges, the error of each image also accumulates. The scheme used in MERFISH allows detection of errors and replacement with valid sequence [48, 51]. MERFISH has been used to spatially resolve 10,050 genes simultaneously at subcellular resolution and to identify specific gene subsets that are enriched in subcellular compartments. seqFISH + applies a complex in situ hybridization technique to obtain higher dimensional results than MERFISH . The primary oligonucleotide probes have sequence complementary to RNAs of interest and a readout region for secondary probe binding. The secondary probes are conjugated to combinatorial fluorophore labels used in imaging to potentially allow detection of 24,000 genes in single cells.
Spatial proteomics methods
The idea that arose from traditional IHC that multiple proteins could be detected in a single step staining procedure prompted development of highly multiplexed spatial proteomic detection (Table 2 and Fig. 1A, middle). Most proteomics platforms either employ fluorophores or metal tags. The method are similar in terms of spatial resolution, cell throughput, number of the molecular targets, and the temporal dynamics. For fluorescence-based methods, a procedure called iterative image acquisition is applied to the FF or FFPE tissue of interest through cyclic addition and removal of fluorescently labeled primary antibodies. Examples include MACSima Imaging Cyclic Staining (MICS) , tissue-based cyclic immunofluorescence (t-CycIF) [54, 55], Co-Detection by IndEXing (CODEX) , Signal Amplification By Exchange Reaction (Immuno-SABER) , and InSituPlex . The choice of fluorophores or design of oligonucleotides is important for the success of spatial detection using these iterative methods. The antibodies are conjugated to fluorophores that have minimal spectral overlap (maximal Stokes shift), to indexing oligonucleotides (CODEX), to orthogonal DNA concatemers (Immuno-SABER), or amplified oligonucleotides (InSituPlex) to provide specificity and sensitivity during the iterative detection process. To bleach or inactivate the label, t-CycIF and MICS use gentle conditions such as photobleaching with specialized antibodies (REAfinity or REAdye_lease) or a specialized removal reagent (REAlease). In Multi Omic Single-scan Assay with Integrated Combinatorial Analysis (MOSAICA), secondary probes are conjugated to the combinatorial fluorophore labels used in imaging . Spatial results are captured by fluorescence lifetime imaging and microscopy (FLIM). Using FLIM, fluorescence spectral detection and lifetime measurements are obtained for each pixel. Fluorescence spectral and lifetime data are processed by a machine learning-based decoding method, and phasor analysis is utilized for mapping spectral and temporal information to original images yielding three-dimensional tissue images. CODEX has the advantage of not requiring an amplification enzyme or specialized buffers, which makes it less costly and less time-consuming .
The metal tag-based spatial proteomics methods are imaging mass cytometry (IMC)  and Multiplex Ion Beam Imaging (MIBI) . In these strategies, metal-conjugated antibodies are used to stain tissue samples, and target proteins are detected using mass spectrometry based on the abundances of isotopic reporter masses released from the tissue upon ablation with a laser beam or ion beams. The output is as a non-overlapping mass signal integration for each measured cell. A single laser is used in IMC to ablate tissue. The two source ion beams employed in MIBI result in higher spatial resolution than is obtained with IMC [61,62,63]. Both IMC and MIBI have been used to study tumor samples from pre-clinical and clinical studies and have significantly contributed to our understanding of the highly complex architecture of TIME at a cellular level and its role in tumorigenesis [62, 64,65,66,67,68,69,70]. Multiplexed fluorescence- and metal-based tagging methods can be employed to detect protein and mRNA within the same tissue. Co-detection strategies aim to reveal both genotypic and phenotypic information about a cell simultaneously .
Spatial metabolic methods
Mass spectrometry is a robust technique for multiplexed analysis of proteins, natural products, and metabolic derivatives [72, 73]. Since classical mass spectroscopy methods do not provide spatial information, mass spectrometry-based imaging strategies have been developed that employ different ionization methods . These include matrix-assisted laser desorption/ionization mass spectrometry (MALDI) [75, 76], desorption electrospray ionization (DESI) [77, 78], and secondary ion mass spectrometry (SIMS) [79, 80]. These platforms allow label-free detection, quantification, and mapping of multiple metabolites including small molecules, lipids, peptides, organic compounds, and elemental ions in cells and tissues (Table 3, Fig. 1A, lower panel).
MALDI requires applying and crystallizing the sample onto a matrix. The resolution of current commercial MALDI instrumentation is about 10 μm, and atmospheric pressure MALDI provides resolution below 2 μm . To improve ionization efficiency, MALDI-2 combines a laser post-ionization with MALDI . Transmission-mode MALDI-2 (t-MALDI-2) provides resolution of 1–2 μm for detection of phospholipids and some biomolecules , and 0.6-μm spatial resolution was obtained by adapting the strategy to the Orbitrap mass analyzer .
DESI has been used to analyze drug, biological and metabolic molecules under ambient conditions [85,86,87,88]. For DESI, a matrix is not necessary; charged droplets and ions of solvent are sprayed directly onto the surface of analyte. The analyte on the surface is taken up by a stream of charged solvent droplets to form the multiply charged ions analyzed by mass spectroscopy. The spatial resolution of DESI is about 50–200 µm. Signal levels decrease and sensitivity is reduced when large areas of tissue are analyzed, such as whole-body sections of mice. To address this and expand the coverage of metabolites, air flow-assisted ionization was incorporated into the DESI workflow in a pipeline called air flow-assisted DESI (AFADESI) . The resolution of AFADESI is 300–500 µm. Another DESI-based protocol, nanoDESI, which incorporates a solvent bridge between primary and nanospray capillaries, has spatial resolution of 10–15 µm .
SIMS facilitates soft ionization of analytes via a primary ion beam. As SIMS is a high vacuum technique, sample preparation for SIMS usually requires chemical or cryogenic fixation to maintain tissue integrity . Several methods have been developed based on SIMS. TOF–SIMS combines the time-of-flight (TOF) and SIMS to obtain information on the molecular layers of a solid surface to increase resolution . The spatial resolution of TOF–SIMS can reach 1 μm. 3D OrbiSIMS is label-free and has subcellular lateral resolution (0.3 μm) and a high mass resolving power [92, 93]. Despite limitations in terms of sample type and complicated sample preparation, 3D OrbiSIMS enables visualization of exogenous and endogenous metabolites in tissues in three-dimensions. The spatial single nuclear metabolomics (SEAM) method was developed to solve problems in segmentation and representation in SIMS data . SEAM preserves the native state of samples with fast and minimal sample processing, providing in situ metabolic fingerprints and single nuclei clustering.
Studies of tumor progression using spatial ‘omics
Researchers have used spatial ‘omics to zoom in on cancer, uncovering various features of the TIME. Studies of the genomes, transcriptomes, proteomes, and metabolomes of normal and cancerous tissues have revealed differences in cellular compositions in different compartments, tumor-immune cell interactomes, and correlates of tumor progression (Fig. 2). The following sections describe how spatial ‘omics tools have been used to image normal and tumor tissues in preclinical and clinical samples.
Spatial imaging reveals compartmentalization of tumor and immune cells
Solid tumors are composed of tumor cells, immune cells, stromal cells, and a vascular system. These components are not evenly distributed. Tumor cell-rich and stroma-rich areas are evident under traditional histological examination. For example, breast cancer has well characterized histopathological features and serves a good model for spatial cellular profiling. This cancer has been studied using spatial transcriptomics [95,96,97] and high-dimensional antibody-based tissue profiling [67, 70]. Using scRNA-seq, Wu et al. profiled breast cancer tissue cell types and the spatial localization of tumor-associated immune and stromal cells . Based on the spatial annotation, the authors found that T cells were mostly located in an area rich in lymphocytes and in an area composed of stromal cells and lymphocytes. The locations of CD4+ and CD8+ T cells were positively correlated with APOE+ macrophages which also express PD-L1 and PD-L2, suggesting an active immunoregulation. At the invasive tumor front, both cxcl10-expressing macrophages and CD8+ T cells were detected.
In another study, Andersson et al. used spatial transcriptomics to profile HER2+ breast cancer tissue . They discovered that tumor and immune cells are compartmentalized. Immune cells are either intermingled with connective tissue or clustered as a lymphoid organ structure, some surrounding the ductal carcinoma in situ (DCIS) area. Similar findings were reported by Salmen et al., who also used a spatial transcriptomics approach to study HER2+ breast cancer specimens . They also found enrichment of T helper 2 and class-switched memory B cells in the DCIS area. In agreement with findings reported by Wu et al. , Salmen et al. detected strong associations between macrophages that express CXCL9 and CXCL10 and T cells and NK cells that express CXCR3.
Antibody-based, single-cell protein expression profiling using MIBI was used by Keren et al. to analyze locations of 36 proteins in specimens from patients with triple-negative breast cancer (TNBC) . The proteins included those expressed by tumor cells, immune cells, and stromal cells. Functional markers for proliferation and epigenetic and immune regulation were also included. In agreement with spatial transcriptomics approaches, B cell presence was highly correlated with CD4+ and CD8+ T cells. A cell neighborhood analysis revealed a clear tumor-to tumor neighbor and immune-to-immune neighbor patterns. These findings reflect the tumor tissue compartmentalization seen in spatial transcriptomics data.
Spatial proteomics studies have also revealed this compartmentalization of tumor and immune cells in several other cancers. Topological distribution of immune cells within the hepatocellular carcinoma (HCC) specimens was studied by Sheng et al. using IMC . Based on 36 antibodies identifying tumor, stromal, immune, and endothelial cells, they discovered several tumor-stromal-immune cell clusters that reflected coordinated cellular neighborhoods and compartmentalization of tumor and immune cells. Most immune cells (infiltrated macrophages, T and B cells) were absent from the cancerous region or limited to the perivascular region. Instead, these cells clustered within an immune-fibroblastic area. Sheng et al. also discovered a spatial relationship between macrophages and T cells. The liver-resident macrophages, the Kupffer cells, were enriched in the peritumoral area and were mainly immunosuppressive. T cells surrounding the Kupffer cells expressed higher level of PD-1 than those surrounding the infiltrating macrophages.
Chan et al. studied small cell lung cancer (SCLC) samples using the MIBI platform . SCLC has long been considered an “immune-cold” cancer due to little infiltration of leukocytes into the tumor. Quantitatively, the immune-mixing score determined using MIBI was indeed lower within the SCLC tumors than in other types of cancer. The SCLC tumors expressed NEUROD, which was associated with poor prognosis, and Chan et al. also found that CD8+ T cells from SCLC tumor samples were more exhausted and had lower effector-like gene expression than did tumors of other subtypes.
Immunological features are used to classify patients with colorectal cancer (CRC) into either the Crohn’s-like reaction (CLR) subtype or the diffuse inflammatory infiltration (DII) subtype. The immune cell topologies in CRC samples of both subtypes were studied by Schurch et al. using CODEX . Tertiary lymphoid structures (TLSs) are present at the tumor invasive front in CLR tumors but are absent in tumors from DII patients; the former patients have much longer overall survival. Cell neighborhood coordination was evaluated by analysis of the CODEX metadata. As suggested by histological analyses, CODEX analysis showed that immune-tumor cellular neighborhoods were more prevalent in the DII type than in the more compartmentalized CLR type. The CODEX images also showed that CD8+ T cell proliferation was more frequently observed in the CLR type tumors, whereas immunosuppressive regulatory T cell proliferation was more frequently found in the DII type. In CLR, direct communication between T cells and macrophages that express activity markers such as ICOS were frequently found in tumor boundary areas but not the bulk tumor areas. In contrast, immunosuppressive granulocytes were more active in the DII type than the CLR type tumor boundary areas. Schurch et al. also found higher frequencies of T cells in the bulk CLR-type tumors than in DII bulk tumors and enrichment of immunosuppressive granulocytes and macrophages in the DII subtype tumors. These data suggest that increased number of immunosuppressive cells and lack of compartmentalized interactions between immune subsets might lead to worse outcomes in DII patients compared to patients with CLR-type tumors.
Immune cells show more exhaustion markers when distributed in the tumor parenchyma
In a spatially resolved transcriptomic analysis of tumor from patients with cholangiocarcinoma, Wu et al. observed that T cells in the tumor core expressed higher levels of exhaustion markers than those in the tumor boundary or peritumoral normal tissue . The same pattern was detected in our recent unpublished CODEX analysis of a mouse model of B cell lymphoma (Fig. 3). We detected enrichment of exhausted CD8+ tumor-infiltrating lymphocytes (TILs) in the tumor compartment relative to the surrounding tissue (Fig. 3A and B). A cell neighborhood coordination analysis also showed that the cluster of exhausted CD8+ TILs (cluster 33) strongly interacts with two tumor cell clusters (clusters 28 and 32) as well as a cluster of CD4+ TILs (cluster 21) (Fig. 3C). In contrast, our data indicate that cells of the non-exhausted CD8+ TIL cluster (cluster 29) interact with macrophages (cluster 0) and non-tumor cells (cluster 24) (Fig. 3C).
Spatial multi-omics approaches reveal active tumor-immune interactions at the invasive front
Invasive fronts, where tumor cells border immune cells, are usually the areas of active immune reactions. A spatial transcriptomics stratification reported by Anderssen et al. revealed more interaction among various cell types at the front than in the tumor core in HER2+ breast cancer . Further, Anderssen et al. found that tumors in some patients are enriched for natural killer T cells, CD8+ effector memory T cells, T helper 1 cells, CD4+ naïve T cells, and memory B cells at the invasive front and that these patients have better survival outcomes than those without prominent infiltration. A spatial proteomics analysis of TNBC by Keren et al. revealed a gradient of histone methylation activity along the tumor-immune border, indicating an active chromatin status when tumor cells are close to the border . They also detected CD11c+/CD11b+ myeloid-derived suppressor cells expressing PD-1, PD-L1, and IDO at the invasive front, suggesting that the front is an immunosuppressive site in this type of breast cancer. In agreement with Anderssen et al., Keren et al. showed that immune cell infiltration at the tumor-immune border was correlated with better prognosis.
In a study of HCC and intrahepatic cholangiocarcinoma (ICC) tumors, Wu et al. carefully examined the cellular composition of the tumor border using Stereo-seq . They found that NK cells, T cells, and macrophages accumulate close to the boundary between the tumor and normal tissue. NK cells and T cells were detected within 250–500 μm from the border and express immune checkpoint genes. Macrophages, mainly anti-inflammatory (M2-like), were located within 250 μm of the border on the tumor side. These results indicate that there is a unique microenvironment at the invasive tumor front. Detailed subgrouping of the hepatocytes at the borders of these tumors revealed a group of hepatocytes that express cellular damaging markers including SAA1 accumulate in the border region. Spatially, SAA1-expressing hepatocytes were detected along the border in close vicinity to macrophages. Since the recruited macrophages were mostly the M2-like phenotype, this recruitment might facilitate invasion of cancerous cells into normal tissue and, in the long run, worsen patient prognosis.
Behaviors of the tumor-associated stromal cells and fibroblasts are spatially correlated
Cancer-associated stromal and fibroblastic cell compositions are heterogenous and have correlated spatial distributions in the tumor tissue. Combining scRNA-seq and spatial transcriptomics, Moncada et al. studied the cellular distribution and functional status of these cell types in pancreatic ductal adenocarcinoma tissue sections . They identified a cancer cell cluster highly associated with fibroblasts in which both cancer cells and fibroblasts are enriched for inflammatory gene signatures. The spatial transcriptomics data also revealed that M1 inflammatory macrophages are associated with a cancer cell subpopulation that is enriched for stress-response genes and associated with inflammatory fibroblasts. Similarly in ICC tumors, fibroblasts that express matrix-forming genes and inflammatory, antigen-presenting features were present in the peritumoral area .
In breast cancer, cancer-associated fibroblasts (CAFs) have been mapped using a spatial transcriptomics approach and scRNA-seq . The so-called inflammatory CAFs (iCAFs) that express ALDH1A1, KLF4, LEPR, and CXCL12 were distributed across the areas that contain both tumor and immune-stromal cells, whereas the myofibroblast CAFs were enriched only in the tumor cell-enriched area. The iCAFs were also anatomically correlated with areas enriched in memory and naïve B cells and CD4+ and CD8+ T cells. CAFs have long been considered immune regulatory cells that mobilize against tumor antigens. Analyses of spatial transcriptomic data revealed a strong chemokine, complement pathway and TGFβ-mediated cross-talk between the iCAFs and the neighboring T cells.
Risom et al. utilized the MIBI to profile the stromal component evolution from the early DCIS to invasive breast cancer . They found that DCIS without progression was characterized by fewer fibroblasts and discontinuous myoepithelium surrounding the DCIS. In samples of tumors that progressed into invasive cancer, the DCIS was associated with active fibroblasts, fiber formation, and absence of myoepithelium. Further, E-Cadherin expression by the myoepithelium and its continuity around the tumor were highly associated with tumor malignancy.
Spatial cellular profiling identifies tumor-associated tertiary lymphoid structures
TLSs are locations where tumor antigens are presented to T cells . Due to the heterogeneous composition of the lymphoid structures, the identification of the tumor-associated TLSs can be difficult without the power of multiplexed immunostaining and transcriptomics information. TLSs harbor multiple types of immune cells with distinct functions. Spatial analyses can identify aggregates of B cells, T cells, and dendritic cells that are characteristic of these structures. Studies of breast tumors revealed a B cell-enriched TLS in TNBC patients and a TLS in HER+ breast cancer patients with a high degree of co-localization between B and T cells [70, 96]. A TLS structure was also found in the CLR type of CRC that is not present in the DII type; in the latter the PD-1+/ICOS+/CD4+ T cells were enriched in TLSs . The presence of TLSs is an indicator of better prognosis for patients with CRC.
Tumor tissue topological profiling identifies novel tumor biomarkers
Like other multi-omics approaches, spatial transcriptomics and scRNA-seq methods have the potential to lead to discovery of new tumor biomarkers. As an example of this potential, Gouin et al. used scRNA-seq, spatial transcriptomics methods, and CODEX to study specimens collected from 25 muscle-invasive bladder cancer patients . From the single-cell clustering data, they discovered a novel biomarker, N-Cadherin 2 (encoded by CDH12), which is indicative of a cancerous epithelium cluster associated with stem-cell-like feature and poor prognosis. Cellular neighborhood analyses from spatial transcriptomics and CODEX both revealed that the T cells in the CDH12-expressing cancer cell vicinity expressed more exhaustion markers. These findings, supported by data from other types of cancer , show that tumor cells with more de-differentiated, more proliferative, and more aggressive features are usually associated in situ with exhausted T cells.
The immune landscape of the metastatic tumor
The spatial interactions of metastatic tumor cells and the immune cells have also been explored. Liver metastases of CRC were investigated by Wu et al. using spatial transcriptomics and multiplexed IHC . These authors discovered a specific macrophage population that expresses MRC1+ and CCL18+ that resembles the resident Kupffer cells. This finding coincides with the role of the resident Kupffer cells in the primary HCC tumors . These macrophages, which had anti-inflammatory M2-like expression features, as well as neutrophils were enriched in the liver metastases compared to the primary lesion. Spatially, these macrophages were scattered in the metastatic tumors and peritumoral boundaries. Ligand-receptor interaction analyses indicated extensive cross-talk between the metastatic tumor cells and the MRC1+/CCL18+ macrophages; an example is the “don’t-eat-me” signaling between CD47 on tumor cells and SIRPA on the macrophages.
The T cell landscape of solid tumor brain metastases was studied by Sudmeier et al. . Using a spatial transcriptomics approach, they found a consistent percentage of PD-1-expressing CD8+ T cells among the metastatic tumor cells. They also detected variable expression of exhaustion markers and lesser T cell receptor diversity in metastatic lesions than in the primary tumor. Spatially, exhausted T cells were distributed on the periphery of the metastatic tumor rather than in the tumor core. Progenitor-like T cells were located adjacent to the inflamed peritumoral area, suggesting that differential cytokine or chemokine signaling within the tumor could serve as a differential niche for CD8+ T cells.
Tumor and immune cell metabolism profiling
Tumor-specific metabolic pathways have been identified by enrichment analysis of spatially resolved transcriptomics data from ICC samples  and of spatially resolved metabolic analyses of esophageal squamous cell carcinomas . The former study identified activation of hypoxia-related pathways and metabolism through the tricarboxylic acid cycle as well as upregulation of fatty acid metabolism components (e.g., fatty acyl CoA synthesis and fatty acid beta-oxidation of tumor cells), high proliferative capacity, and high levels of apoptosis that likely reflect higher proliferative capacity, damaged states, and energy requirements of tumor cells. The spatial metabolic platform revealed upregulation of arginine and proline metabolism, fatty acid biosynthesis, and alanine, aspartate, glutamate, pyrimidine, and histidine metabolism in tumors. Six critical metabolic enzymes within four pathways were present at higher levels in tumor compartments than in normal tissue: the pyrroline-5-carboxylate reductase PYCR2, glutaminase, uridine phosphorylase UPase1, fatty acid synthase, and ornithine decarboxylase. These enzymes catalyze proline biosynthesis, glutamine catabolism, phosphorolytic cleavage of uridine to uracil, and decarboxylation of histidine to form histamine, respectively. PYCR2 is an essential enzyme in proline biosynthesis and promotes cancer proliferation and progression ; however, histamine is derived from the decarboxylation of histidine, which is catalyzed by an enzyme dramatically down-regulated in cancer , and histamine-based therapies can lead to cancer cell apoptosis and senescence and prolong survival in tumor-bearing animals .
scRNA-seq and spatial transcriptomics have also been applied to a recently developed mouse model of human neuroblastoma, revealing that the spatial relationship of CD4+ and CCR2+ macrophages play a pro-tumor role via the arginine metabolic pathway . Single-cell metabolic profiling also identified a tumor-specific metabolic phenotype characterized by high levels of amino acid transporter CD98 expression and showed that a tumor-associated metabolic T cell state is characterized by expression of exhaustion markers PD1 and CD39  as well as downregulated levels of TCF1, which is indicative of terminal exhaustion .
Intratumor cellular topologic patterns are correlated with patient prognosis
To understand the topological cellular interaction patterns in tissues, high-dimensional multi-omics data are analyzed computationally for correlations between locations of cells of particular types. The terms “cellular neighborhood” and “ecotype”, among others, have been used by independent research groups to describe the coordinated presence or absence of particular cell types [95, 98, 112]. Certain patterns of cellular coordination are linked to patient survival or prognosis in various cancers. Gouin et al. found that the CD8+ T cells within CDH12-enriched epithelial cellular neighborhoods in bladder tumors expressed high levels of CD49a, PD-1, or LAG3 . These data may explain why the patients who have higher CDH12 levels have better response to atezolizumab treatment and longer overall survival.
The spatial topology of the TIME is also correlated with overall survival, and compartmentalization of the tumor-immune architecture is usually positively correlated with better survival. In diffuse large B cell lymphoma, Colombo et al. showed that the more structured germinal center B cell subtype was correlated with better overall survival than the dispersed subtype . In CRC, tumors from the CLR subtype had more separated compartments than the tumors from patients with the DII subtype, and the presence of the compartments were statistically correlated with patient prognosis . The phenomenon was also observed in the breast cancer tumors where compartmentalized score is highly correlated with overall survival . In another example, Andersson et al. defined a TLS score in breast cancer specimens based on the degree of co-localization between B and T cells and found that a higher TLS score was also associated with better overall survival .
Changes in compartmentalization occur during tumor progression
The development of spatial ‘omics technologies has made it possible to evaluate the spatial distribution of cell subsets, cell characteristics, and metabolism among different compartments of tumors during the course of cancer progression. For example, in HCC patients, Sheng et al. found that levels and locations of dedifferentiation, proliferation, and immune checkpoint markers evolved during the process of tumorigenesis . In CRC tissue, Hartmann et al. observed metabolic polarization of immune cells toward the tumor-immune boundary with increased expression of CD98 and ASCT2 . CD98 and ASCT2 are transporters with prognostic value in human cancers [114,115,116].
Furthermore, key genes in the progression of tumors can be identified by multi-dimensional deconstruction of spatial architecture datasets. Risom et al. obtained MIBI data on DCIS and invasive breast cancer (IBC) samples as well as normal breast tissue and extracted the top 20 most important features during the progression from DCIS to IBC [67, 117]. Hypoxia, glycolysis, stromal immune density, and desmoplasia/remodeling of the extracellular matrix were enriched in DCIS. Myoepithelial and immunoregulatory markers (e.g., PDL1, IDO1, COX2, PD1) on tumor and immune cells were enriched in IBC. Interestingly, the expression patterns of ECAD, SMA, CK5, and myoepithelial markers were similar in normal tissue and in IBC tissues, and the highest myoepithelial expression of ECAD was observed in normal breast tissue. These findings suggest that the loss of “normal” features may serve a protective function in non-progressors. Further, myoepithelial loss in the stroma surrounding cancer cells appears to induce fibroblast and immune cell activation, playing a critical role in determining clinical outcome. In summary, spatial analyses suggest that invasive progression depends on an evolving spatial distribution of multiple cell types rather than alteration in levels or distribution of a single cell subset.
Studies using spatial multi-omics tools have revealed the complexity of the TIME and have shown that relative positions and interactions of cell types in the microenvironment of tumors, in addition to the cellular composition, strongly influence tumor development. A better understanding of spatial interactions is driving redefinition of tumor subtypes and shifting the focus of research to tumor-immune interaction units, the discovery of additional cell types, and the examination of the changes in the TIME compartment as cancer progresses.
Spatial ‘omics technologies are in the “boom” period of development. Scientists in this field are addressing technological obstacles that limit resolution of the platforms, multiplexing, sensitivity, and accuracy. In terms of resolution, currently commercialized spatial transcriptomic and proteomics technologies do not provide true single-cell level resolution. In several studies, parallel scRNA-seq analyses of the same specimen have been performed and used to re-assign and adjust the spatial transcriptomics data . This may be a way to solve the resolution problem until higher resolution methods become widely available.
Multiplex antibody staining was a great leap forward from traditional IHC, which only allows simultaneous imaging of three or four targets. Multiplexing results in loss of sensitivity, however. For example, traditional IHC has a significantly higher detection sensitivity than cytometry-based multiplexed imaging . Current multiplexed imaging protocols rely heavily on computational methods for data processing. The segmentation process, which translates the image data into cellular data, is critical and can suffer from inaccuracies. Computational segmentation requires prior selection of markers of cellular identity (e.g., nucleus and membrane border). Training data are also required to demarcate irregularly shaped or polymorphous cells (e.g., endothelial cells, fibroblasts). In addition, interpretation may be confounded by partially overlapping cells or non-specific background staining, which may affect the results of subsequent analyses and should be examined carefully before the data processing. At the time of preparation of this manuscript, new methods had recently been reported that yield three-dimensional images [119, 120]. In the future, we expect that a powerful and easy-to-handle protocol will be developed to investigate cellular heterogeneity from multilayered tumor specimens.
In the last decade, spatial ‘omics developments have led to progress toward personalized, precision medicine. Precision medicine, in which the treatment strategy is customized for each individual patient, relies heavily on detailed patient data. Thanks to technical advances, scientists are now able to analyze tumor tissue at the single-cell level with high data dimensions. The spatial multi-omics tools described in this review have revealed the heterogeneous composition of tumor and immune cells in the tumor microenvironment. These approaches have enabled comprehensive explorations of cancer and have furthered our understanding of the underlying mechanisms of tumor progression, which are critical for treatment and response. In the future, with the implementation of robust automated pipelines, clinicians and pathologists should be able to evaluate disease progression and tailor the therapeutic regimens to each patient, bringing us closer to the goal of precision medicine.
Availability of data and materials
The authors declare that all data supporting the findings of this study are available from the corresponding author on reasonable request.
Tumor immune microenvironment
Immune checkpoint blockers
Single-cell RNA sequencing
Cytometry by time-of-flight
Fluorescence in situ hybridization
High-Definition Spatial Transcriptomics
Deterministic Barcoding in Tissue sequencing
Regions of interest
Digital Spatial Profiler
Fluorescent in situ sequencing
Sequencing by Oligonucleotide Ligation and Detection
Spatially resolved Transcript Amplicon Readout mapping
Sequencing with Error-reduction by Dynamic Annealing and Ligation
Multiplexed error-robust FISH
Single-molecule fluorescence in situ hybridization
MACSima Imaging Cyclic Staining
Tissue-based cyclic immunofluorescence
Co-Detection by IndEXing
Immuno-signal amplification by exchange reaction
Multi Omic Single-scan Assay with Integrated Combinatorial Analysis
Fluorescence lifetime imaging and microscopy
Imaging mass cytometry
Multiplex Ion Beam Imaging
Matrix-assisted laser desorption/ionization–mass spectrometry
Desorption electrospray ionization
Secondary ion mass spectrometry
Transmission mode MALDI
Air flow-assisted DESI
Spatial single nuclear metabolomics
Myeloid derived suppressor cell
Ductal carcinoma in situ
Triple-negative breast cancer
Small cell lung cancer
Crohn’s-like reaction type
Diffuse inflammatory infiltration type
Tertiary lymphoid structure
Tumor infiltrating lymphocyte
Invasive breast cancer
Xia L, Oyang L, Lin J, Tan S, Han Y, Wu N, Yi P, Tang L, Pan Q, Rao S, et al. The cancer metabolic reprogramming and immune response. Mol Cancer. 2021;20(1):28.
Kaymak I, Williams KS, Cantor JR, Jones RG. Immunometabolic interplay in the tumor microenvironment. Cancer Cell. 2021;39(1):28–37.
O’Donnell JS, Teng MWL, Smyth MJ. Cancer immunoediting and resistance to T cell-based immunotherapy. Nat Rev Clin Oncol. 2019;16(3):151–67.
Ugel S, Canè S, De Sanctis F, Bronte V. Monocytes in the tumor microenvironment. Annu Rev Pathol. 2021;16:93–122.
Biffi G, Tuveson DA. Diversity and biology of cancer-associated fibroblasts. Physiol Rev. 2021;101(1):147–76.
Sobierajska K, Ciszewski WM, Sacewicz-Hofman I, Niewiarowska J. Endothelial cells in the tumor microenvironment. Adv Exp Med Biol. 2020;1234:71–86.
Sharma P, Siddiqui BA, Anandhan S, Yadav SS, Subudhi SK, Gao J, Goswami S, Allison JP. The next decade of immune checkpoint therapy. Cancer Discov. 2021;11(4):838–57.
Kubli SP, Berger T, Araujo DV, Siu LL, Mak TW. Beyond immune checkpoint blockade: emerging immunological strategies. Nat Rev Drug Discov. 2021;20(12):899–919.
Pan D, Jia D. Application of single-cell multi-omics in dissecting cancer cell plasticity and tumor heterogeneity. Front Mol Biosci. 2021;8: 757024.
Lei Y, Tang R, Xu J, Wang W, Zhang B, Liu J, Yu X, Shi S. Applications of single-cell sequencing in cancer research: progress and perspectives. J Hematol Oncol. 2021;14(1):91.
Trapnell C. Defining cell types and states with single-cell genomics. Genome Res. 2015;25(10):1491–8.
Bagger FO, Probst V. Single cell sequencing in cancer diagnostics. Adv Exp Med Biol. 2020;1255:175–93.
Zhang Y, Wang D, Peng M, Tang L, Ouyang J, Xiong F, Guo C, Tang Y, Zhou Y, Liao Q, et al. Single-cell RNA sequencing in cancer research. J Exp Clin Cancer Res. 2021;40(1):81.
Ren X, Zhang L, Zhang Y, Li Z, Siemers N, Zhang Z. Insights gained from single-cell analysis of immune cells in the tumor microenvironment. Annu Rev Immunol. 2021;39:583–609.
Bassler K, Schulte-Schrepping J, Warnat-Herresthal S, Aschenbrenner AC, Schultze JL. The myeloid cell compartment-cell by cell. Annu Rev Immunol. 2019;37:269–93.
Bendall SC, Simonds EF, Qiu P, el Amir AD, Krutzik PO, Finck R, Bruggner RV, Melamed R, Trejo A, Ornatsky OI, et al. Single-cell mass cytometry of differential immune and drug responses across a human hematopoietic continuum. Science. 2011;332(6030):687–96.
Bodenmiller B, Zunder ER, Finck R, Chen TJ, Savig ES, Bruggner RV, Simonds EF, Bendall SC, Sachs K, Krutzik PO, et al. Multiplexed mass cytometry profiling of cellular states perturbed by small-molecule regulators. Nat Biotechnol. 2012;30(9):858–67.
Spitzer MH, Nolan GP. Mass cytometry: single cells many features. Cell. 2016;165(4):780–91.
Hartmann FJ, Bendall SC. Immune monitoring using mass cytometry and related high-dimensional imaging approaches. Nat Rev Rheumatol. 2020;16(2):87–99.
Hsieh WC, Lai EY, Liu YT, Wang YF, Tzeng YS, Cui L, Lai YJ, Huang HC, Huang JH, Ni HC et al. NK cell receptor and ligand composition influences the clearance of SARS-CoV-2. J Clin Invest. 2021; 131(21).
Wei SC, Levine JH, Cogdill AP, Zhao Y, Anang NAS, Andrews MC, Sharma P, Wang J, Wargo JA, Pe’er D, et al. Distinct cellular mechanisms underlie anti-CTLA-4 and anti-PD-1 checkpoint blockade. Cell. 2017;170(6):1120-1133.e1117.
Newell EW, Sigal N, Bendall SC, Nolan GP, Davis MM. Cytometry by time-of-flight shows combinatorial cytokine expression and virus-specific cell niches within a continuum of CD8+ T cell phenotypes. Immunity. 2012;36(1):142–52.
Becher B, Schlitzer A, Chen J, Mair F, Sumatoh HR, Teng KW, Low D, Ruedl C, Riccardi-Castagnoli P, Poidinger M, et al. High-dimensional analysis of the murine myeloid cell system. Nat Immunol. 2014;15(12):1181–9.
Wagner J, Rapsomaniki MA, Chevrier S, Anzeneder T, Langwieder C, Dykgers A, Rees M, Ramaswamy A, Muenst S, Soysal SD, et al. A single-cell atlas of the tumor and immune ecosystem of human breast cancer. Cell. 2019;177(5):1330-1345.e1318.
Jass JR. Classification of colorectal cancer based on correlation of clinical, morphological and molecular features. Histopathology. 2007;50(1):113–30.
Loi S, Drubay D, Adams S, Pruneri G, Francis PA, Lacroix-Triki M, Joensuu H, Dieci MV, Badve S, Demaria S, et al. Tumor-infiltrating lymphocytes and prognosis: a pooled individual patient analysis of early-stage triple-negative breast cancers. J Clin Oncol. 2019;37(7):559–69.
Adams S, Goldstein LJ, Sparano JA, Demaria S, Badve SS. Tumor infiltrating lymphocytes (TILs) improve prognosis in patients with triple negative breast cancer (TNBC). Oncoimmunology. 2015;4(9): e985930.
Gruosso T, Gigoux M, Manem VSK, Bertos N, Zuo D, Perlitch I, Saleh SMI, Zhao H, Souleimanova M, Johnson RM, et al. Spatially distinct tumor immune microenvironments stratify triple-negative breast cancers. J Clin Invest. 2019;129(4):1785–800.
Lewis SM, Asselin-Labat ML, Nguyen Q, Berthelet J, Tan X, Wimmer VC, Merino D, Rogers KL, Naik SH. Spatial omics and multiplexed imaging to explore cancer biology. Nat Methods. 2021;18(9):997–1012.
Wang N, Li X, Wang R, Ding Z. Spatial transcriptomics and proteomics technologies for deconvoluting the tumor microenvironment. Biotechnol J. 2021;16(9): e2100041.
Rao A, Barkley D, França GS, Yanai I. Exploring tissue architecture using spatial transcriptomics. Nature. 2021;596(7871):211–20.
Ståhl PL, Salmén F, Vickovic S, Lundmark A, Navarro JF, Magnusson J, Giacomello S, Asp M, Westholm JO, Huss M, et al. Visualization and analysis of gene expression in tissue sections by spatial transcriptomics. Science. 2016;353(6294):78–82.
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. 2019;363(6434):1463–7.
Stickels RR, Murray E, Kumar P, Li J, Marshall JL, Di Bella DJ, Arlotta P, Macosko EZ, Chen F. Highly sensitive spatial transcriptomics at near-cellular resolution with Slide-seqV2. Nat Biotechnol. 2021;39(3):313–9.
Lee Y, Bogdanoff D, Wang Y, Hartoularos GC, Woo JM, Mowery CT, Nisonoff HM, Lee DS, Sun Y, Lee J et al. XYZeq: Spatially resolved single-cell RNA sequencing reveals expression heterogeneity in the tumor microenvironment. Sci Adv. 2021; 7(17).
Vickovic S, Eraslan G, Salmén F, Klughammer J, Stenbeck L, Schapiro D, Äijö T, Bonneau R, Bergenstråhle L, Navarro JF, et al. High-definition spatial transcriptomics for in situ tissue profiling. Nat Methods. 2019;16(10):987–90.
Cho CS, Xi J, Si Y, Park SR, Hsu JE, Kim M, Jun G, Kang HM, Lee JH. Microscopic examination of spatial transcriptome using Seq-Scope. Cell. 2021;184(13):3559-3572.e3522.
Chen A, Liao S, Cheng M, Ma K, Wu L, Lai Y, Qiu X, Yang J, Xu J, Hao S, et al. Spatiotemporal transcriptomic atlas of mouse organogenesis using DNA nanoball-patterned arrays. Cell. 2022;185(10):1777-1792.e1721.
Wu L, Yan J, Bai Y, Chen F, Xu J, Zou X, Huang A, Hou L, Zhong Y, Jing Z et al. Spatially-resolved transcriptomics analyses of invasive fronts in solid tumors. bioRxiv 2021:2021.2010.2021.465135.
Hu KH, Eichorst JP, McGinnis CS, Patterson DM, Chow ED, Kersten K, Jameson SC, Gartner ZJ, Rao AA, Krummel MF. ZipSeq: barcoding for real-time mapping of single cell transcriptomes. Nat Methods. 2020;17(8):833–43.
Liu Y, Yang M, Deng Y, Su G, Enninful A, Guo CC, Tebaldi T, Zhang D, Kim D, Bai Z, et al. High-spatial-resolution multi-omics sequencing via deterministic barcoding in tissue. Cell. 2020;183(6):1665-1681.e1618.
Geiss GK, Bumgarner RE, Birditt B, Dahl T, Dowidar N, Dunaway DL, Fell HP, Ferree S, George RD, Grogan T, et al. Direct multiplexed measurement of gene expression with color-coded probe pairs. Nat Biotechnol. 2008;26(3):317–25.
Merritt CR, Ong GT, Church SE, Barker K, Danaher P, Geiss G, Hoang M, Jung J, Liang Y, McKay-Fleisch J, et al. Multiplex digital spatial profiling of proteins and RNA in fixed tissue. Nat Biotechnol. 2020;38(5):586–99.
Maliga Z, Nirmal AJ, Ericson NG, Boswell SA, U’Ren L, Podyminogin R, Chow J, Chen Y-A, Chen AA, Weinstock DM et al. Micro-region transcriptomics of fixed human tissue using Pick-Seq. bioRxiv 2021:2021.2003.2018.431004.
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. Nat Methods. 2013;10(9):857–60.
Lee JH, Daugharthy ER, Scheiman J, Kalhor R, Yang JL, Ferrante TC, Terry R, Jeanty SS, Li C, Amamoto R, et al. Highly multiplexed subcellular RNA sequencing in situ. Science. 2014;343(6177):1360–3.
Wang X, Allen WE, Wright MA, Sylwestrak EL, Samusik N, Vesuna S, Evans K, Liu C, Ramakrishnan C, Liu J et al. Three-dimensional intact-tissue sequencing of single-cell transcriptional states. Science. 2018; 361(6400).
Chen KH, Boettiger AN, Moffitt JR, Wang S, Zhuang X. RNA imaging. Spatially resolved, highly multiplexed RNA profiling in single cells. Science. 2015;348(6233):6090.
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(10):877–9.
Femino AM, Fay FS, Fogarty K, Singer RH. Visualization of single RNA transcripts in situ. Science. 1998;280(5363):585–90.
Xia C, Fan J, Emanuel G, Hao J, Zhuang X. Spatial transcriptome profiling by MERFISH reveals subcellular RNA compartmentalization and cell cycle-dependent gene expression. Proc Natl Acad Sci U S A. 2019;116(39):19490–9.
Eng CL, Lawson M, Zhu Q, Dries R, Koulena N, Takei Y, Yun J, Cronin C, Karp C, Yuan GC, et al. Transcriptome-scale super-resolved imaging in tissues by RNA seqFISH. Nature. 2019;568(7751):235–9.
Kinkhabwala A, Herbel C, Pankratz J, Yushchenko DA, Rüberg S, Praveen P, Reiß S, Rodriguez FC, Schäfer D, Kollet J, et al. MACSima imaging cyclic staining (MICS) technology reveals combinatorial target pairs for CAR T cell treatment of solid tumors. Sci Rep. 2022;12(1):1911.
Lin JR, Fallahi-Sichani M, Sorger PK. Highly multiplexed imaging of single cells using a high-throughput cyclic immunofluorescence method. Nat Commun. 2015;6:8390.
Lin J-R, Izar B, Wang S, Yapp C, Mei S, Shah PM, Santagata S, Sorger PK. Highly multiplexed immunofluorescence imaging of human tissues and tumors using t-CyCIF and conventional optical microscopes. Elife. 2018;7: e31657.
Black S, Phillips D, Hickey JW, Kennedy-Darling J, Venkataraaman VG, Samusik N, Goltsev Y, Schürch CM, Nolan GP. CODEX multiplexed tissue imaging with DNA-conjugated antibodies. Nat Protoc. 2021;16(8):3802–35.
Saka SK, Wang Y, Kishi JY, Zhu A, Zeng Y, Xie W, Kirli K, Yapp C, Cicconet M, Beliveau BJ, et al. Immuno-SABER enables highly multiplexed and amplified protein imaging in tissues. Nat Biotechnol. 2019;37(9):1080–90.
Manesse M, Patel KK, Bobrow M, Downing SR. The InSituPlex(®) staining method for multiplexed immunofluorescence cell phenotyping and spatial profiling of tumor FFPE samples. Methods Mol Biol. 2020;2055:585–92.
Vu T, Vallmitjana A, Gu J, La K, Xu Q, Flores J, Zimak J, Shiu J, Hosohama L, Wu J, et al. Spatial transcriptomics using combinatorial fluorescence spectral and lifetime encoding, imaging and analysis. Nat Commun. 2022;13(1):169.
Giesen C, Wang HA, Schapiro D, Zivanovic N, Jacobs A, Hattendorf B, Schüffler PJ, Grolimund D, Buhmann JM, Brandt S, et al. Highly multiplexed imaging of tumor tissues with subcellular resolution by mass cytometry. Nat Methods. 2014;11(4):417–22.
Angelo M, Bendall SC, Finck R, Hale MB, Hitzman C, Borowsky AD, Levenson RM, Lowe JB, Liu SD, Zhao S, et al. Multiplexed ion beam imaging of human breast tumors. Nat Med. 2014;20(4):436–42.
Keren L, Bosse M, Thompson S, Risom T, Vijayaragavan K, McCaffrey E, Marquez D, Angoshtari R, Greenwald NF, Fienberg H, et al. MIBI-TOF: a multiplexed imaging platform relates cellular phenotypes and tissue structure. Sci Adv. 2019;5(10):eaax5851.
Rovira-Clavé X, Jiang S, Bai Y, Zhu B, Barlow G, Bhate S, Coskun AF, Han G, Ho CK, Hitzman C, et al. Subcellular localization of biomolecules and drug distribution by high-definition ion beam imaging. Nat Commun. 2021;12(1):4628.
Kuett L, Catena R, Özcan A, Plüss A, Schraml P, Moch H, de Souza N, Bodenmiller B. Three-dimensional imaging mass cytometry for highly multiplexed molecular and cellular mapping of tissues and the tumor microenvironment. Nat Cancer. 2022;3(1):122–33.
Hoch T, Schulz D, Eling N, Gómez JM, Levesque MP, Bodenmiller B. Multiplexed imaging mass cytometry of the chemokine milieus in melanoma characterizes features of the response to immunotherapy. Sci Immunol. 2022;7(70):eabk1692.
Jackson HW, Fischer JR, Zanotelli VRT, Ali HR, Mechera R, Soysal SD, Moch H, Muenst S, Varga Z, Weber WP, et al. The single-cell pathology landscape of breast cancer. Nature. 2020;578(7796):615–20.
Risom T, Glass DR, Averbukh I, Liu CC, Baranski A, Kagel A, McCaffrey EF, Greenwald NF, Rivero-Gutiérrez B, Strand SH, et al. Transition to invasive breast cancer is associated with progressive changes in the structure and composition of tumor stroma. Cell. 2022;185(2):299-310.e218.
Allam M, Cai S, Coskun AF. Multiplex bioimaging of single-cell spatial profiles for precision cancer diagnostics and therapeutics. NPJ Precis Oncol. 2020;4:11.
Devine RD, Behbehani GK. Mass cytometry, imaging mass cytometry, and multiplexed ion beam imaging use in a clinical setting. Clin Lab Med. 2021;41(2):297–308.
Keren L, Bosse M, Marquez D, Angoshtari R, Jain S, Varma S, Yang SR, Kurian A, Van Valen D, West R, et al. A structured tumor-immune microenvironment in triple negative breast cancer revealed by multiplexed ion beam imaging. Cell. 2018;174(6):1373-1387.e1319.
Jiang S, Chan CN, Rovira-Clavé X, Chen H, Bai Y, Zhu B, McCaffrey E, Greenwald NF, Liu C, Barlow GL, et al. Combined protein and nucleic acid imaging reveals virus-dependent B cell and macrophage immunosuppression of tissue microenvironments. Immunity. 2022;55(6):1118-1134.e1118.
Kowalczyk T, Ciborowski M, Kisluk J, Kretowski A, Barbas C. Mass spectrometry based proteomics and metabolomics in personalized oncology. Biochim Biophys Acta Mol Basis Dis. 2020;1866(5): 165690.
Zhou W, Liotta LA, Petricoin EF. Cancer metabolism and mass spectrometry-based proteomics. Cancer Lett. 2015;356(2 Pt A):176–83.
Murayama C, Kimura Y, Setou M. Imaging mass spectrometry: principle and application. Biophys Rev. 2009;1(3):131.
Stoeckli M, Chaurand P, Hallahan DE, Caprioli RM. Imaging mass spectrometry: a new technology for the analysis of protein expression in mammalian tissues. Nat Med. 2001;7(4):493–6.
Kaufmann R. Matrix-assisted laser desorption ionization (MALDI) mass spectrometry: a novel analytical tool in molecular biology and biotechnology. J Biotechnol. 1995;41(2–3):155–75.
Takáts Z, Wiseman JM, Gologan B, Cooks RG. Mass spectrometry sampling under ambient conditions with desorption electrospray ionization. Science. 2004;306(5695):471–3.
Takáts Z, Wiseman JM, Cooks RG. Ambient mass spectrometry using desorption electrospray ionization (DESI): instrumentation, mechanisms and applications in forensics, chemistry, and biology. J Mass Spectrom. 2005;40(10):1261–75.
Pacholski ML, Winograd N. Imaging with mass spectrometry. Chem Rev. 1999;99(10):2977–3006.
Sjövall P, Lausmaa J, Johansson B. Mass spectrometric imaging of lipids in brain tissue. Anal Chem. 2004;76(15):4271–8.
Kompauer M, Heiles S, Spengler B. Atmospheric pressure MALDI mass spectrometry imaging of tissues and cells at 14-μm lateral resolution. Nat Methods. 2017;14(1):90–6.
Soltwisch J, Kettling H, Vens-Cappell S, Wiegelmann M, Müthing J, Dreisewerd K. Mass spectrometry imaging with laser-induced postionization. Science. 2015;348(6231):211–5.
Zavalin A, Todd EM, Rawhouser PD, Yang J, Norris JL, Caprioli RM. Direct imaging of single cells and tissue at sub-cellular spatial resolution using transmission geometry MALDI MS. J Mass Spectrom. 2012;47(11):1473–81.
Niehaus M, Soltwisch J, Belov ME, Dreisewerd K. Transmission-mode MALDI-2 mass spectrometry imaging of cells and tissues at subcellular resolution. Nat Methods. 2019;16(9):925–31.
Kauppila TJ, Talaty N, Kuuranne T, Kotiaho T, Kostiainen R, Cooks RG. Rapid analysis of metabolites and drugs of abuse from urine samples by desorption electrospray ionization-mass spectrometry. Analyst. 2007;132(9):868–75.
Kauppila TJ, Wiseman JM, Ketola RA, Kotiaho T, Cooks RG, Kostiainen R. Desorption electrospray ionization mass spectrometry for the analysis of pharmaceuticals and metabolites. Rapid Commun Mass Spectrom. 2006;20(3):387–92.
Reynolds JC, Blackburn GJ, Guallar-Hoyas C, Moll VH, Bocos-Bintintan V, Kaur-Atwal G, Howdle MD, Harry EL, Brown LJ, Creaser CS, et al. Detection of volatile organic compounds in breath using thermal desorption electrospray ionization-ion mobility-mass spectrometry. Anal Chem. 2010;82(5):2139–44.
Wiseman JM, Ifa DR, Venter A, Cooks RG. Ambient molecular imaging by desorption electrospray ionization mass spectrometry. Nat Protoc. 2008;3(3):517–24.
Luo Z, He J, Chen Y, He J, Gong T, Tang F, Wang X, Zhang R, Huang L, Zhang L, et al. Air flow-assisted ionization imaging mass spectrometry method for easy whole-body molecular imaging under ambient conditions. Anal Chem. 2013;85(5):2977–82.
Roach PJ, Laskin J, Laskin A. Nanospray desorption electrospray ionization: an ambient method for liquid-extraction surface sampling in mass spectrometry. Analyst. 2010;135(9):2233–6.
Spraker JE, Luu GT, Sanchez LM. Imaging mass spectrometry for natural products discovery: a review of ionization methods. Nat Prod Rep. 2020;37(2):150–62.
Makarov A. Electrostatic axially harmonic orbital trapping: a high-performance technique of mass analysis. Anal Chem. 2000;72(6):1156–62.
Passarelli MK, Pirkl A, Moellers R, Grinfeld D, Kollmer F, Havelund R, Newman CF, Marshall PS, Arlinghaus H, Alexander MR, et al. The 3D OrbiSIMS-label-free metabolic imaging with subcellular lateral resolution and high mass-resolving power. Nat Methods. 2017;14(12):1175–83.
Yuan Z, Zhou Q, Cai L, Pan L, Sun W, Qumu S, Yu S, Feng J, Zhao H, Zheng Y, et al. SEAM is a spatial single nuclear metabolomics method for dissecting tissue microenvironment. Nat Methods. 2021;18(10):1223–32.
Wu SZ, Al-Eryani G, Roden DL, Junankar S, Harvey K, Andersson A, Thennavan A, Wang C, Torpy JR, Bartonicek N, et al. A single-cell and spatially resolved atlas of human breast cancers. Nat Genet. 2021;53(9):1334–47.
Andersson A, Larsson L, Stenbeck L, Salmén F, Ehinger A, Wu SZ, Al-Eryani G, Roden D, Swarbrick A, Borg Å, et al. Spatial deconvolution of HER2-positive breast cancer delineates tumor-associated cell type interactions. Nat Commun. 2021;12(1):6012.
Salmén F, Vickovic S, Larsson L, Stenbeck L, Vallon-Christersson J, Ehinger A, Häkkinen J, Borg Å, Frisén J, Ståhl PL et al. Multidimensional transcriptomics provides detailed information about immune cell distribution and identity in HER2+ breast tumors. bioRxiv 2018:358937.
Sheng J, Zhang J, Wang L, Tano V, Tang J, Wang X, Wu J, Song J, Zhao Y, Rong J, et al. Topological analysis of hepatocellular carcinoma tumour microenvironment based on imaging mass cytometry reveals cellular neighbourhood regulated reversely by macrophages with different ontogeny. Gut. 2022;71(6):1176–91.
Chan JM, Quintanal-Villalonga Á, Gao VR, Xie Y, Allaj V, Chaudhary O, Masilionis I, Egger J, Chow A, Walle T, et al. Signatures of plasticity, metastasis, and immunosuppression in an atlas of human small cell lung cancer. Cancer Cell. 2021;39(11):1479-1496.e1418.
Schürch CM, Bhate SS, Barlow GL, Phillips DJ, Noti L, Zlobec I, Chu P, Black S, Demeter J, McIlwain DR, et al. Coordinated cellular neighborhoods orchestrate antitumoral immunity at the colorectal cancer invasive front. Cell. 2020;182(5):1341-1359.e1319.
Moncada R, Barkley D, Wagner F, Chiodin M, Devlin JC, Baron M, Hajdu CH, Simeone DM, Yanai I. Integrating microarray-based spatial transcriptomics and single-cell RNA-seq reveals tissue architecture in pancreatic ductal adenocarcinomas. Nat Biotechnol. 2020;38(3):333–42.
Sautès-Fridman C, Petitprez F, Calderaro J, Fridman WH. Tertiary lymphoid structures in the era of cancer immunotherapy. Nat Rev Cancer. 2019;19(6):307–25.
Wu Y, Yang S, Ma J, Chen Z, Song G, Rao D, Cheng Y, Huang S, Liu Y, Jiang S, et al. Spatiotemporal immune landscape of colorectal cancer liver metastasis at single-cell level. Cancer Discov. 2022;12(1):134–53.
Sudmeier LJ, Hoang KB, Nduom EK, Wieland A, Neill SG, Schniederjan MJ, Ramalingam SS, Olson JJ, Ahmed R, Hudson WH. Distinct phenotypic states and spatial distribution of CD8(+) T cell clonotypes in human brain metastases. Cell Rep Med. 2022;3(5): 100620.
Sun C, Li T, Song X, Huang L, Zang Q, Xu J, Bi N, Jiao G, Hao Y, Chen Y, et al. Spatially resolved metabolomics to discover tumor-associated metabolic alterations. Proc Natl Acad Sci U S A. 2019;116(1):52–7.
Liu W, Hancock CN, Fischer JW, Harman M, Phang JM. Proline biosynthesis augments tumor cell growth and aerobic glycolysis: involvement of pyridine nucleotides. Sci Rep. 2015;5:17206.
Ruffell B, Coussens LM. Histamine restricts cancer: nothing to sneeze at. Nat Med. 2011;17(1):43–4.
Martinel Lamas DJ, Cortina JE, Ventura C, Sterle HA, Valli E, Balestrasse KB, Blanco H, Cremaschi GA, Rivera ES, Medina VA. Enhancement of ionizing radiation response by histamine in vitro and in vivo in human breast cancer. Cancer Biol Ther. 2015;16(1):137–48.
Van de Velde LA, Allen EK, Crawford JC, Wilson TL, Guy CS, Russier M, Zeitler L, Bahrami A, Finkelstein D, Pelletier S, et al. Neuroblastoma formation requires unconventional CD4 T cells and arginase-1-dependent myeloid cells. Cancer Res. 2021;81(19):5047–59.
Hartmann FJ, Mrdjen D, McCaffrey E, Glass DR, Greenwald NF, Bharadwaj A, Khair Z, Verberk SGS, Baranski A, Baskar R, et al. Single-cell metabolic profiling of human cytotoxic T cells. Nat Biotechnol. 2021;39(2):186–97.
Blank CU, Haining WN, Held W, Hogan PG, Kallies A, Lugli E, Lynn RC, Philip M, Rao A, Restifo NP, et al. Defining ‘T cell exhaustion.’ Nat Rev Immunol. 2019;19(11):665–74.
Gouin KH 3rd, Ing N, Plummer JT, Rosser CJ, Ben Cheikh B, Oh C, Chen SS, Chan KS, Furuya H, Tourtellotte WG, et al. An N-Cadherin 2 expressing epithelial cell subpopulation predicts response to surgery, chemotherapy and immunotherapy in bladder cancer. Nat Commun. 2021;12(1):4906.
Colombo AR, Hav M, Singh M, Xu A, Gamboa A, Lemos T, Gerdtsson E, Chen D, Houldsworth J, Shaknovich R, et al. Single-cell spatial analysis of tumor immune architecture in diffuse large B-cell lymphoma. Blood Adv. 2022;6(16):4675–90.
Kaira K, Oriuchi N, Imai H, Shimizu K, Yanagitani N, Sunaga N, Hisada T, Ishizuka T, Kanai Y, Nakajima T, et al. Prognostic significance of L-type amino acid transporter 1 (LAT1) and 4F2 heavy chain (CD98) expression in stage I pulmonary adenocarcinoma. Lung Cancer. 2009;66(1):120–6.
Furuya M, Horiguchi J, Nakajima H, Kanai Y, Oyama T. Correlation of L-type amino acid transporter 1 and CD98 expression with triple negative breast cancer prognosis. Cancer Sci. 2012;103(2):382–9.
Shimizu K, Kaira K, Tomizawa Y, Sunaga N, Kawashima O, Oriuchi N, Tominaga H, Nagamori S, Kanai Y, Yamada M, et al. ASC amino-acid transporter 2 (ASCT2) as a novel prognostic marker in non-small cell lung cancer. Br J Cancer. 2014;110(8):2030–9.
Risom T, Glass DR, Averbukh I, Liu CC, Baranski A, Kagel A, McCaffrey EF, Greenwald NF, Rivero-Gutierrez B, Strand SH, et al. Transition to invasive breast cancer is associated with progressive changes in the structure and composition of tumor stroma. Cell. 2022;185(2):299–310.
Halse H, Colebatch AJ, Petrone P, Henderson MA, Mills JK, Snow H, Westwood JA, Sandhu S, Raleigh JM, Behren A, et al. Multiplex immunohistochemistry accurately defines the immune context of metastatic melanoma. Sci Rep. 2018;8(1):11158.
Lin J-R, Wang S, Coy S, Tyler M, Yapp C, Chen Y-A, Heiser CN, Lau KS, Santagata S, Sorger PK. Multiplexed 3D atlas of state transitions and immune interactions in colorectal cancer. bioRxiv 2021:2021.2003.2031.437984.
Zhang Z, Tang R, Chen X, Waller L, Kau A, Fung AA, Gutierrez B, An C, Cho SH, Shi L, et al. A high-throughput technique to map cell images to cell positions using a 3D imaging flow cytometer. Proc Natl Acad Sci U S A. 2022;119(8):e2118068119.
Work in our laboratory has been supported by a Ministry of Science and Technology, Taiwan Grant (MOST109-2320-B-001–017-MY3) and Academia Sinica, Taiwan Grants (AS-CDA-110-L09, AS-GC-110-05, AS-KPQ-110-EIMD, AS-GCS-111-L03, VTA111-V3-1-2), and NTUH 110-T20.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Hsieh, WC., Budiarto, B.R., Wang, YF. et al. Spatial multi-omics analyses of the tumor immune microenvironment. J Biomed Sci 29, 96 (2022). https://doi.org/10.1186/s12929-022-00879-y