- 1Department of Integrative Oncology, British Columbia Cancer Research Centre, Vancouver, BC, Canada
- 2Department of Biochemistry and Molecular Pharmacology, New York University (NYU) Langone Medical Center, New York, NY, United States
The growth and metastasis of solid tumours is known to be facilitated by the tumour microenvironment (TME), which is composed of a highly diverse collection of cell types that interact and communicate with one another extensively. Many of these interactions involve the immune cell population within the TME, referred to as the tumour immune microenvironment (TIME). These non-cell autonomous interactions exert substantial influence over cell behaviour and contribute to the reprogramming of immune and stromal cells into numerous pro-tumourigenic phenotypes. The study of some of these interactions, such as the PD-1/PD-L1 axis that induces CD8+ T cell exhaustion, has led to the development of breakthrough therapeutic advances. Yet many common analyses of the TME either do not retain the spatial data necessary to assess cell-cell interactions, or interrogate few (<10) markers, limiting the capacity for cell phenotyping. Recently developed digital pathology technologies, together with sophisticated bioimage analysis programs, now enable the high-resolution, highly-multiplexed analysis of diverse immune and stromal cell markers within the TME of clinical specimens. In this article, we review the tumour-promoting non-cell autonomous interactions in the TME and their impact on tumour behaviour. We additionally survey commonly used image analysis programs and highly-multiplexed spatial imaging technologies, and we discuss their relative advantages and limitations. The spatial organization of the TME varies enormously between patients, and so leveraging these technologies in future studies to further characterize how non-cell autonomous interactions impact tumour behaviour may inform the personalization of cancer treatment.
1 Introduction
The solid tumour microenvironment (TME) is a tumour-supporting niche that encompasses a diverse population of cells, including malignant, stromal, endothelial, bacterial, and immune cells. Within the larger TME, the population of immune cells is referred to as the tumour immune microenvironment (TIME). The TME is highly heterogeneous and can impact the progression of tumours, as well as their response to various therapeutics. The various cells of the TME influence tumour growth and invasiveness in numerous ways, including secreting cytokine factors that can activate cytotoxic effector cells, polarize immune cells, and drive both local and distal inflammation (1). The precise impacts of these cytokines on tumour biology vary between different cancer types and stages of tumour progression (1). Tumour microenvironments often have distinct physical and chemical traits, such as chronic hypoxia, which develops as a result of cancer cell growth outpacing the circulatory system’s capacity for expansion (2), and acidity, which arises largely as a consequence of high rates of glycolysis and the export of lactate and H+ ions by cancer cells (3). Many solid tumour microenvironments also contain distinct bacterial microbiomes, with taxonomic compositions that differ from those seen in healthy organs (4, 5).
Exposure to these unique characteristics of the TME, as well as to the cytokines and other signaling molecules secreted by malignant cells, frequently reprograms tumour-infiltrating and surrounding cells toward phenotypes and polarizations that further advantage the cancer cells. Reprogrammed stromal cells, such as cancer-associated fibroblasts (CAFs), create a highly cross-linked, stiff extracellular matrix (ECM) that promotes cell migration and invasion (6, 7). Reprogrammed cells of the TME, such as tumour-associated macrophages (TAMs), engage immunosuppressive programs that limit the anti-tumour activities of T cells (8). In this way there is extensive crosstalk through metabolites, cell surface receptor interactions, and secreted signaling molecules between cancer cells and other tumour-associated cells.
Specific features of an individual patient’s TME are now recognized to significantly impact their prognosis. Survival outcomes have been linked to a variety of factors, including diversity of species within the tumour microbiome (4), abundance or subtypes of CAFs (9, 10), or extent and distribution of tumour hypoxia (11), but the most well-described associations with patient outcomes are those pertaining to abundance, phenotypes, and distribution of immune cells both in the TME and within tumours themselves. DNA-, RNA-, and protein-level alterations acquired during malignant transformation can cause cancer cells to produce antigens with unique peptide sequences (“neoantigens”) (12, 13). Tumours that have greater genomic instability and a larger number of mutations (i.e. a higher mutational burden) are expected to display more neoantigens. Since neoantigens are not protected by central tolerance, their display on the surface of cancer cells can lead to T cell recognition and precipitate an anti-tumour immune response that improves patient outcomes (12, 14). Consequently, high quantities of tumour-infiltrating CD8+ T cells, and of cells that stimulate their function (e.g. CD4+ T cells, conventional type 1 dendritic cells, and B cells), are generally associated with relatively good prognoses (15–18). Conversely, large quantities of regulatory T cells (Tregs), myeloid-derived suppressor cells (MDSCs), and TAMs, all of which broadly impede anti-tumour immune responses, have been linked to poor prognoses in multiple cancer types (19–21). Furthermore, low levels of CD8+ T cells increase the risk of recurrence or progression for ductal carcinomas in situ and oral leukoplakias, suggesting that the composition of the immune cell infiltrate may also be predictive of the behaviour of precancerous lesions (22, 23).
Therapeutic strategies that target the TME, including immune checkpoint blockade (ICB) therapies, have profoundly impacted cancer patient survival rates, due in large part to their remarkable efficacy and even curative capability in some of the deadliest and most prevalent cancer types, such as melanoma and lung cancer. However, ICB induces responses in only a minority of patients, and is broadly ineffective against tumours that evade immune surveillance by preventing the generation or infiltration of tumour-specific T cells (24, 25). Patient responses to more widely-used treatments such as chemotherapy, radiation therapy, and various targeted therapies can also be impaired by TME-related factors, including low quantities of immune cell infiltration, dense stroma and poor circulation, and large areas of chronic hypoxia (2, 26, 27). The tremendous variance in the spatial organization of the TME between patients, even those with tumours driven by the same oncogene(s), remains a major confounding factor to treatment decisions and clinical outcomes. As an example, the confinement of immune cells to the tumour stroma and their exclusion from the tumour core has been associated with a decreased likelihood of response to anti-PD-L1 therapy (28), as well as poorer overall survival (29).
The spatial organization of TME-resident cells also influences the frequencies of non-cell autonomous interactions. These interactions can be driven by cell-cell contact, secreted proteins, or metabolites, and they contribute to cancer phenotypes. However, many methods used to study the TME, such as flow cytometry, only provide information about the abundance of cell subtypes and do not capture their spatial distribution. Imaging-based methods such as multiplex immunofluorescence (mIF) do record these spatial data, but the capacity for these technologies to deeply characterize interactions between specific cell subtypes is limited by their inability to visualize more than 6-8 markers in a given sample (30). Consequently, the recent creation of imaging technologies that enable highly multiplexed, high-resolution spatial profiling of solid TMEs is impactful, and has the potential to elevate the study of the malignant immune microenvironment from measuring cell abundance to investigating “cell sociology” – the myriad relationships, interactions, and communications between cells (31). In this article, we review the study of cell sociology in the TME through (a) summarizing the well-described non-cell autonomous interactions within the TME and their impact on tumour behaviour and therapeutic response, (b) discussing the image analysis programs and computational methods currently used to decipher cell sociology in the TME, and (c) reviewing recently developed highly-multiplexed spatial imaging technologies.
2 Non-cell autonomous interactions within tumors
2.1 Non-cell autonomous interactions mediated by cell-cell contact
The development of multiplex immunohistochemistry (mIHC) and mIF has made it possible to quantify how often cells of different types exist close to each other in the tumour microenvironment and consequently to estimate the frequency of specific cell-cell contacts. As we will discuss, physical contacts between cancer, immune, and bacterial cells within the TME play significant roles in the modulation of the anti-tumour immune response and cancer cell proliferation (Supplementary Table 1).
CD8+ T cells are among the most critical immune components of the TME, due to their ability to recognize, bind, and kill neoantigen-expressing tumour cells. However, tumour cells are only susceptible to this destruction if they come into direct contact with CD8+ cells. In metastatic colon cancer, having a high percentage of tumour epithelial cells located in proximity to CD8+ T cells has been linked to increased survival, despite overall CD8+ T cell infiltration being lower in these patients’ tumours (32). Similarly, the likelihood of lung adenocarcinoma recurrence has been observed to correlate more strongly with the frequency with which CD8+ cells neighbour tumour epithelial cells than with the overall CD8+ cell density (31). Underscoring the importance of immune cell subtyping, relapse-free survival in triple-negative breast cancer (TNBC) has been linked specifically to high levels of CD8+CD103+ T cells in immediate proximity to cancer cells (33). CD8+ T cell function is also mediated by other immune cells; for example, an increased proximity of CD8+ T cells to CD3+CD8-FOXP3- helper T cells (32) and B cells (34) has been associated with improved patient outcomes.
Direct interactions between T cells, B cells, and dendritic cells occur in tertiary lymphoid structures (TLSs), which are dense clusters of immune cells commonly present in chronically inflamed areas, including the tumour microenvironments of multiple cancers (35). TLSs contain both T and B cell regions and are transient sites of the germinal centre reactions, which lead to B cell differentiation (35). T follicular helper (Tfh) cells are significant components of TLSs, being involved in both TLS formation and the germinal centre reactions (36, 37). While the production of IL-21 and CXCL13 by Tfh cells can stimulate adaptive anti-tumour immunity at a distance (38), the direct engagement of Tfh cells with B cells through ICOS/ICOSL and CD40L/CD40 binding is also critical to the anti-tumour immune response (39).
The presence of TLSs has been linked to improved prognoses in a number of cancers, including melanoma (40), head and neck squamous cell carcinoma (41), pancreatic ductal adenocarcinoma (PDAC) (42), and muscle-invasive bladder cancer (43), and can have greater prognostic value than the bulk count of infiltrating CD8+ T cells (44–46). Tumours with TLSs tend to have distinct features, including increased infiltration by CD20+ B cells and both CD8+ and CD4+ T cells (44, 47), decreased infiltration of CD163+ M2 macrophages (48), decreased expression of the Treg markers FOXP3 and CD25 by CD4+ T cells (48, 49), and higher levels of Bcl-2 and lower levels of TIM3 expressed by T cells (40). These associations may in part reflect an increased likelihood of TLS development in the presence of chronic inflammation caused by a pre-existing anti-tumour immune response, as suggested by Cabrita et al. (40), but are also due to the local amplification of that response by TLSs through recruiting lymphocytes and facilitating antigen presentation and lymphocyte maturation (35, 50, 51) (Figure 1A). Supporting the idea that TLSs have a regional impact on the immune response, shorter distances between TLSs and the invasive front of bladder tumours have been linked to increased disease-specific survival (43). There is also some evidence that intra-tumoural TLSs are more prognostically favourable than peri-tumoural TLSs (48, 53).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1275890/fimmu-14-1275890-HTML/image_m/fimmu-14-1275890-g001.jpg)
Figure 1 Direct contacts between cells in the TME modulate tumour phenotypes. (A) Cell-cell interactions within TLSs, particularly the germinal centre, influence TME composition and immune phenotypes. TLS image is reproduced from (52) with permission of the authors; (B) Interactions between F nucleatum and TME-resident cells have pro-tumour effects; (C) T cell anti-tumour activity is attenuated by direct contact with Tregs and PD-L1+ cells.
Tumour-resident bacteria within the TME can influence tumour phenotypes in numerous ways, including through cell-cell contacts with cancer and immune cells. One of the best-characterized examples is Fusobacterium nucleatum, a Gram-negative bacterium that is normally present in the oral cavity, but that has been found within the TME of a number of cancer types, most prominently colorectal cancer (CRC) (54) (Figure 1B). The Fap2 protein of F. nucleatum can bind the immune checkpoint protein TIGIT, commonly expressed by NK and T cells, limiting their cytotoxic activity (55). F. nucleatum also expresses FadA adhesin, which can bind the E-cadherin expressed by CRC cells. This induces the nuclear translocation of β-catenin and downstream expression of oncogenic proteins, including Wnt, Myc, and cyclin D1, which stimulate tumour growth (56). F. nucleatum additionally exerts a number of pro-tumour effects through interactions between lipopolysaccharide (LPS) and host-expressed TLR4, including upregulation of hsa-miR-21-5p, which stimulates cancer cell growth and invasion (57), and downregulation of hsa-miR-18a-3p and hsa-miR-4802, which induces autophagy and confers resistance to oxaliplatin and 5-fluorouracil (58).
2.2 Influence of cell-cell contacts on immunotherapy response
As mentioned, the targeting of the PD-1/PD-L1 axis by ICB therapy represents a breakthrough in cancer therapy for thousands of patients, but the majority of treated patients do not experience responses (Figure 1C) (59). It has been frequently proposed that response to ICB is more common in the case where PD-L1 expression is not simply constitutive, but is instead induced by cytokines produced during an adaptive immune response (25, 60). Proxies for the presence or likelihood of an adaptive immune response, including tumour mutational burden, PD-L1 expression, various signatures of inflammatory gene expression, and biomarkers based on mIHC or mIF have all been correlated with response to anti-PD-1 or PD-L1 therapy, with a recent meta-analysis by Lu et al. finding mIHC/mIF biomarkers to be the most accurate predictors of therapeutic response (61). These spatial biomarkers offer unique insight into whether PD-L1+ cells are enriched in the vicinity of functional targets (in which case inhibition of the PD-1/PD-L1 axis is thought to be more likely to amplify anti-tumour immunity) (60) or randomly distributed throughout the TME, and can only be measured through multiplexed imaging. For instance, high levels of proximity between PD-L1+ cells and either PD-1+ cells (25, 62, 63), exhausted CD8+ cells (64), or cancer cells (65) have been associated with improved outcomes after ICB therapy. Similarly, high expression of PD-L1 by M1 macrophages that were located near both CD8+ T cells and the tumour-stroma boundary has been correlated with increased likelihood of response of metastatic melanoma to ICB therapy (66).
While adaptive expression of PD-L1 may represent a therapeutic opportunity, the PD-1/PD-L1 interaction is not uniformly associated with improved patient outcomes. Frequent interactions between PD-L1+ cells and CD8+ cells or PD-1+ cells have been correlated with poor outcomes in HPV- oral and oropharyngeal squamous cell carcinomas, while overall PD-L1+ abundance has not (67, 68). Interestingly, a study of diffuse large B cell lymphoma found that PD-1+/PD-L1+ interactions were associated with poor outcomes in patients with high infiltration of CD3+ cells, but improved outcomes in patients with low CD3+ infiltration (69). These improved outcomes may arise because PD-1+/PD-L1+ interactions indicate that the few CD3+ cells present are concentrated near tumour cells, which, despite the consequent immune checkpoint activation, is more favourable than them being isolated (69). This highlights how the prognostic interpretation of cell-cell interaction scores should be informed by the overall immune context of the tumour.
2.3 Non-cell autonomous interactions mediated by intermediaries
CD4+ and CD8+ regulatory T cells generally suppress immune responses, and consequently, their infiltration into tumours has been linked to worsened prognoses in many forms of cancer (19). Tregs attenuate anti-tumour immune activity through a variety of mechanisms, including depleting IL-2 from the TME (70), pushing antigen-presenting cells (APCs) towards tolerogenic phenotypes characterized by the downregulation of CD80/86 and the upregulation of IDO1 (71), inducing apoptosis in CD8+ cells via expression of the Fas ligand (72) or granzyme and perforin (73), and secreting immunosuppressive cytokines, including IL-10 and TGF-β (74). Links to poor outcomes have been observed for Tregs located near CD8+ T cells in HPV- oral squamous cell carcinoma (OSCC) (68) and in CRCs with microsatellite instability (75), and for both CD4+ and CD8+ Tregs located near non-small cell lung cancer (NSCLC) cells (76, 77). However, increased proximity between Tregs and CD8+ T cells has been associated with improved prognoses in gastric tumours and NSCLC (76, 78). Taken together, these results suggest not only that Treg/CD8+ interactions have distinct functional repercussions in different cancer types, but also that the dominant functions of tumour-infiltrating Tregs depend on their immediate immune context. As a further example of this, mouse models suggest that lung TLSs commonly contain Tregs, which limit the activation and proliferation of the TLS-localized T cells (79). A study of TLS-localized Tregs in breast cancer linked their presence to poorer patient outcomes, despite finding that the presence of Tregs in the tumour bed had no significant prognostic value (80).
High densities of TAMs are associated with poor outcomes in many types of cancer (21), owing to a range of functions that includes the secretion of growth factors, survival factors, and immunosuppressive cytokines (81). The impact of these growth and survival factors has been observed in NSCLC, where cancer cells undergoing apoptosis were on average located closer to pro-inflammatory M1-polarized macrophages than to anti-inflammatory M2-polarized macrophages, while the reverse was true for proliferative Ki67+ cancer cells (82). Accordingly, higher numbers of cancer cells located in proximity to M2 macrophages were found to be associated with poor patient outcomes, while the opposite was true for M1 macrophages (82). High M2 macrophage/cancer cell proximity has been correlated with poor outcomes in PDAC (83), although in gastric cancer it has been linked to improved outcomes (84). M2 macrophages can also engage CD8+ T cells in lengthy, antigen-specific interactions that do not instigate T cell proliferation, but instead induce an exhausted, PD-1+ phenotype (85). Short distances between CD8+ T cells and HLA-DR- (predominantly M2) macrophages have been associated with decreased survival in melanoma, potentially due to the immunosuppressive nature of the CD8+/M2 interaction (86).
2.4 Generation of pro-tumour metabolite profiles in the TME
Cells within the TME can induce pro-tumour metabolite profiles in a number of ways, including depleting metabolites required for the cytotoxic activity of immune cells, secreting immunosuppressive metabolites, and secreting metabolites that stimulate cancer cell growth and division.
Within the TME, the catabolic action of the IDO1 and TDO enzymes can convert the amino acid tryptophan into kynurenine (87, 88). Kynurenine engenders an immunosuppressive environment by promoting Treg development and decreasing the viability and IFNγ expression of CD8+ T cells, while also directly increasing cancer cell proliferation (Figure 2A) (89–91). IDO1 and TDO can be expressed constitutively by cancer cells, in which case they inhibit lymphocyte infiltration, or expressed by cancer and stromal cells in response to IFNγ, which limits the activity of tumour-infiltrating lymphocytes (92). A number of IDO1 inhibitors have been developed, but initial phase III trials have not shown any clinical benefit from adding IDO1 inhibitors to anti-PD-1 therapy regimens (93). It has been suggested that the relative benefit of IDO1 inhibition may be greater in immune-cold tumours that constitutively express IDO1 (94), as it could enable increased lymphocyte infiltration and in turn amplify the utility of PD-1 or PD-L1 inhibition.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1275890/fimmu-14-1275890-HTML/image_m/fimmu-14-1275890-g002.jpg)
Figure 2 TME metabolite profiles alter tumour behaviour. (A) Production of kynurenine by IDO1+ or TDO+ cells enhances cancer cell proliferation and impairs T cell functionality. (B) Conversion of dietary fibre into butyrate by intestinal bacteria suppresses colorectal cancer cell proliferation. (C) Lactate shuttling is induced by oxygen gradients, promoting resistance to angiogenesis-inhibiting treatments. (D) Pancreatic cancer cells induce autophagy in nearby pancreatic stellate cells, liberating alanine, which the cancer cells use to fuel the TCA cycle.
The metabolic activity of commensal bacteria can also influence the tumour microenvironment. Intratumoural bacteria that express the long isoform of cytidine deaminase, such as Enterobacteriaceae, are found commonly in human PDACs and have been shown to convert the chemotherapeutic agent gemcitabine into an inactive form, which confers drug resistance in mouse models (95). Conversely, various species of gut bacteria participate in the conversion of primary bile acids to secondary bile acids and the fermentation of dietary fibre into short-chain fatty acids (SCFAs). The SCFA butyrate, which is exclusively synthesized by the gut microbiome, has histone deacetylase inhibitor activity at high concentrations and is known to inhibit CRC cell proliferation and promote apoptosis through a range of mechanisms, including upregulating miR-203 (96), limiting ERK phosphorylation (97, 98), and upregulating p21 (Figure 2B) (97).
The metabolic composition of the TME is also altered by hypoxia, as it prevents cells from obtaining energy through oxidative phosphorylation and induces, among other transcriptional programs and adaptive phenotypes (99), compensatory upregulation of the glucose importer GLUT1 and the lactate exporter MCT4 (100). Expression of GLUT1 and MCT4 is thus highest in tumour regions that are far from blood vessels (101). This metabolic shift results in an excess of lactate within the TME, which can then be imported by MCT1-expressing cells within normoxic tumour regions and used as an energy source or anabolic building block (102–104) (Figure 2C). This lactate shuttle facilitates the survival of cancer cells in hypoxic regions through decreasing the glucose requirements of cells in nearby normoxic neighbourhoods (102). Establishment of this shuttle is an observed mechanism of resistance to angiogenesis-inhibiting therapy in a number of tumour types (105–107). The secretion of lactate by cancer cells also contributes to the reprogramming of immune cells, including through promoting the M2 polarization in macrophages (108, 109) and augmenting Treg induction (110).
Interactions between cancer cells and CAFs can induce HIF-1 and redox-mediated expression of MCT1 in cancer cells and both GLUT1 and MCT4 in CAFs, thereby creating a CAF-cancer cell lactate shuttle (103, 111, 112). Cancer cells involved in this shuttle generally exhibit MCT1-dependent increases in proliferative (103, 113, 114), invasive (115, 116), and migratory (116) capacity. Ippolito et al. have shown that the exposure of prostate cancer cells to CAF-conditioned media increases their mitochondrial mass and oxygen consumption, which were further augmented by the uptake of mitochondria from neighbouring CAFs along cytoplasmic bridges (115). Prostate cancer patients with both MCT1high cancer cells and MCT4high stromal cells tend to have later stage (pT3) tumours (117), and experience poor, stage-independent, biochemical failure-free survival (118). Neither of these associations were seen in patients with high expression of only one MCT (117, 118), which highlights the clinical significance of lactate shuttling.
Similarly, an alanine shuttle is induced by interactions between pancreatic stellate cells (PSCs) and PDAC cells. Exposure to PDAC-conditioned media causes PSCs to undergo autophagy and secrete alanine, which can be imported by PDAC cells and used as a carbon source for the synthesis of TCA cycle metabolites (119) (Figure 2D). Knockdown of the alanine importer SLC38A2 abrogates the increase in tumour growth seen when PSCs are injected into mouse xenograft models in addition to PDAC cells, which suggests that alanine shuttling is responsible for this growth advantage (120).
3 Emerging approaches for quantifying cell-cell interactions
Most of the multiplexed spatial data collected from tumour samples is in the form of multi-channel images derived from mIHC or mIF, imaged by either absorption or fluorescence microscopy. These images must undergo processing, typically including spectral unmixing, selection of regions of interest (ROIs), nuclear and cell segmentation, and cell phenotyping (Figure 3A) (30).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1275890/fimmu-14-1275890-HTML/image_m/fimmu-14-1275890-g003.jpg)
Figure 3 Distance and neighbour-based methods quantify cell-cell interactions differently. (A) Extracting rates of cell-cell interactions from imaging data requires image preprocessing, which includes cell segmentation and phenotyping. Raw multi-channel image reproduced from (31) with permission of the authors. (B) Distance-based metrics calculate interaction frequencies by considering a pair of cells to interact if they are in physical proximity (e.g. within a given radius of one another). (C) In contrast, neighbour-based metrics consider a pair of cells to interact if they are neighbours (e.g. if they share a Voronoi edge), regardless of the distance between them. In the examples illustrated here, frequencies of interaction between CD3+CD8- cells and three different cell types are calculated.
Cell segmentation and phenotyping performed during image processing enables subsequent extraction and analysis of cell sociology features, most notably the extent of interaction between different cell types (e.g. A+ and B+ cells). Cell-cell interactions are most commonly quantified using metrics that depend on intercellular distances, including the mean or median distance from an A+ cell to the nearest B+ cell (32, 65, 77, 86), the absolute number of B+ cells located within a specified distance of one or more A+ cells (63, 67, 68, 82, 121), the percentage of A+ cells with a B+ cell located within a specified distance (32, 84), and the area under the curve of the A+-B+ G-cross function (76, 122) (Figure 3B). The choice of metric, which can also require choosing a specific threshold distance (commonly between 10 and 40 µm) (32, 63, 67, 68, 82, 84, 121), and the decision of whether to apply cell density normalizations should be guided by the specific biological hypothesis being tested. Distance-based metrics are generally useful for determining the frequency of particular short-range cell-cell interactions in a sample, but they are influenced by cell size and do not identify direct contacts between cells.
The second class of metrics used to quantify the frequency of interaction between cell types can be referred to as ‘neighbour-based’ (Figure 3C). Starting from the positions of cell nuclei, Voronoi diagrams, or equivalently Delaunay triangulations, can be constructed to identify pairs of cells that are immediate neighbours and are likely in physical contact (31, 123). Alternatively, a cell’s neighbours can be identified through dilating its binary mask (62, 69). Neighbour-based metrics include the mean fraction of the neighbours of A+ cells that are B+ (31); the number of B+ cells that neighbour at least one A+ cell, normalized to the total number of immune cells (62, 69); and the ratio of the number of observed A+-B+ neighbour pairs to the expected number of such pairs, based on the proportions of all imaged cells that are A+ or B+ (123). These metrics are less dependent on cell size, but they are not suitable for the analysis of cell-cell interactions that do not require direct contact, and careful specification of ROIs is necessary to avert analysis artifacts, such as cells in isolated areas of a slide being assigned unreasonably large volumes (31). While they aim to measure the same quantity, distance-based and neighbour-based metrics can yield significantly different results (Figures 3B, C). Both classes of metrics have advantages in particular cases, and so the best metric for a given experiment depends on the hypothesis being tested.
With the growth of digital pathology, many bioimage analysis programs, both open and closed-source, have been developed to analyze these metrics (Table 1). In the open-source domain, in addition to ImageJ, researchers now have access to a range of alternatives, including QuPath, CellProfiler, Icy, MCMICRO, Cytomapper, and others (Table 1) (124–129). Open-source tools have become increasingly important, as they allow users to develop, customize, and share their own analysis solutions, commonly in the form of plug-ins, which makes these programs highly flexible and stimulate innovation (130). Closed-source software generally does not permit users to modify, augment, or share the software, but it frequently possesses superior ease of use, available online support, and vertically integrated support for common analysis tasks. HALO (Indica Labs, Albuquerque, NM, USA) and Visiopharm (Visiopharm A/S, Hørsholm, Denmark) are examples of closed-source tools, but there are a variety of others, frequently integrated with image acquisition hardware such as Aperio image analysis (Leica Biosystems, Buffalo Grove, IL, USA) and TissueGnostics Quest suites (TissueGnostics, Vienna, Austria) (Table 1). Each image analysis program offers its own advantages, and identifying the most suitable program for a given analysis project requires consideration of a number of factors, including project objectives, the scale of required analysis, and user experience and training. More recently, programs that focus on the statistical analysis of spatial cell distributions within tissues (cell neighbourhood analysis) have been developed. Many of these (SPIAT, Giotto, histoCAT/NeighbouRhood) (131–134) take only the spatial location of the cells and their quantified characteristics (such as protein marker expression and RNA sequence counts) as input and are predominantly written in R, allowing for user-specific optimization and modifications. Generally, these programs allow the user to define the size (in the number of cell layers, the distance from the cell, or both) of the neighborhood that they want to consider. Giotto in particular allows users to interactively refine cell marker positivity based upon neighbourhood characteristics. As these newer programs demonstrate, the capability of programs available for digital pathology analysis continues to expand, giving researchers access to a constantly evolving toolkit to meet the computational demands of highly-multiplexed imaging technologies.
4 Advantages and applications of highly-multiplexed imaging technologies
As discussed above, the advent of mIHC and mIF has greatly enhanced our understanding of non-cell autonomous mechanisms within the tumour microenvironment. Imaging with a fluorescence microscope can achieve sub-micrometer resolutions, and the use of secondary antibodies and tyramide signal amplification allows for the detection of even lowly-expressed molecular species (135). Yet even with the application of increasingly sophisticated computational analysis, the width and consequent overlapping of the spectra of common fluorophores makes the simultaneous fluorescence imaging of more than 6-8 markers infeasible (30, 136). The TME commonly contains a wide variety of cell types, which may be reprogrammed into various phenotypes or polarizations, and so many relevant questions cannot be answered with only these 6-8 markers. Recently developed imaging technologies can achieve considerably higher levels of multiplexing, which greatly expands their potential applications, but each technology has its own advantages and disadvantages in areas such as imaging resolution, sensitivity, sample throughput, and sample integrity (Figure 4, Tables 2, 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1275890/fimmu-14-1275890-HTML/image_m/fimmu-14-1275890-g004.jpg)
Figure 4 Technologies that preserve spatial information have broad applications. Tumours that carry identical densities of immune cells, and that would appear similar if evaluated by bulk methods, can vary significantly in their spatial organization. Imaging technologies provide a bevy of information on samples, both within multiple-cell ROIs (useful for e.g. assessing intra-tumour heterogeneity) and at the single-cell scale (ideal for e.g. subtyping individual cells or quantifying cell-cell interactions), greatly expanding the range of possible analyses.
4.1 Antibody-based methods
Cyclic immunofluorescence, which broadly involves repeated stain-image-bleach/wash cycles being performed on the same tissue sample, is the simplest extension of mIF. One of the furthest developed forms of this technology is t-CyCIF, in which each imaging cycle visualizes three fluorophore-conjugated primary antibodies and a DNA dye (137). Sample integrity declines with successive imaging cycles, but 60-fold multiplexing has been successfully demonstrated (137). t-CyCIF achieves this multiplexing at high, diffraction-limited resolutions and does not require a specialized microscope, but each imaging cycle is 6-8 hours long (137). Other variants of cyclic IF elute antibodies instead of photobleaching fluorophores (149). This enables the use of secondary antibodies, but harsh elution conditions can alter the structure of some epitopes, and so careful optimization of staining order and antibody choice is necessary (149).
Co-detection by indexing (CODEX) imaging is another cyclic technique, and involves staining with antibodies that are conjugated to oligonucleotide tags with 5’ overhangs of different lengths, followed by multiple cycles of imaging (123). In each cycle, a mix of unlabelled (‘indexing’) and fluorescently labelled dNTPs is added, causing each type of labelled dNTP to be incorporated into a unique tag. CODEX was later re-engineered so that each imaging cycle instead involved the addition of three fluorophore-conjugated oligonucleotides, each complementary to a region on a single antibody’s tag (150). This newer version of CODEX has demonstrated 56-fold multiplexing, and can be performed using a standard fluorescence microscope (150). However, species limitations make the use of secondary antibodies impossible at these high levels of multiplexing, and so imaging sensitivity is limited.
Similar to CODEX, DNA-Exchange-Imaging (DEI) is a cyclic technique in which a fluorophore-conjugated oligonucleotide is introduced during each imaging cycle, and binds to a single antibody-conjugated oligonucleotide (138). DEI is compatible with super-resolution platforms (138), but it shares CODEX’s drawback that the expansion of multiplexing beyond levels achievable through mIF makes the use of secondary antibodies impossible.
Exchange-SABER combines DEI with immunostaining with signal amplification by exchange reaction (Immuno-SABER). In Immuno-SABER, each antibody-conjugated oligonucleotide is bound by a DNA concatemer that has been synthesized via primer exchange reactions. This hybridization creates a lengthy overhang that contains numerous binding sites for short, fluorophore-conjugated ‘imager’ oligonucleotides (139). This approach allows for significant (5 to 180-fold) signal amplification, with the high-end being achieved when secondary and tertiary concatemers are used (Iterative-SABER) (139). Exchange-SABER holds the potential for high levels of multiplexing, but the imaging of more than 10 markers simultaneously has yet to be demonstrated (139).
Multiplexed ion beam imaging (MIBI) involves the preliminary staining of tissues with lanthanide-conjugated antibodies (151). During the imaging process a narrow primary ion beam is scanned across the sample pixel by pixel, which liberates secondary ions from the antibodies, allowing the lanthanide isotopes to be quantified by a mass spectrometer (151). MIBI-TOF, an improved version of MIBI that uses time-of-flight mass spectrometry, has demonstrated simultaneous imaging of up to 36 markers at sub-micrometer resolution, and can approach single-antibody sensitivity (140). However, despite improvements in MIBI-TOF’s throughput, achieving resolutions comparable to mIHC substantially increases the image acquisition time (140).
When used to assess protein expression, digital spatial profiling (DSP) employs antibodies that are conjugated via a UV-cleavable linker to unique oligonucleotides (141). After ROIs are determined, each is sequentially exposed to UV light, and the oligonucleotides released after each exposure are quantified through sequencing or the nCounter system (141). High levels of multiplexing (44 proteins) have been demonstrated, and the ability to establish ROIs of any shape allows for the profiling of highly specific areas of tissue (141). However, using very small (~single-cell) ROIs significantly impairs the detection of lowly-expressed proteins.
Collectively, these antibody-based techniques excel at providing information needed to answer questions about the prevalence, identity, and location of specific immune and stromal cell subtypes within the TME, and about how the proximity or interaction of different cell types impacts tumour behaviour. Demonstrated applications of these technologies include the identification of cellular phenotypes that are consistently located in close spatial proximity in the TNBC TME (121) and the screening of many protein markers simultaneously in distinct compartments of the non-small-cell lung cancer TME, in order to identify protein/compartment pairs that are associated with improved survival outcomes in patients treated with anti-PD-1 checkpoint blockade therapy (152). These technologies can also be combined with computational approaches that operate on a larger spatial scale through treating samples not only as mixtures of single cells but as aggregations of regions (cellular neighbourhoods) with distinct compositions and functions. For instance, the use of clustering and tensor decomposition techniques on CODEX data has enabled the characterization of a granulocyte-enriched neighbourhood within the CRC TME whose functional state, as defined by the expression of PD-1 and CD4 on T cells, was associated with patient outcomes (150). Another example is the use of a recently reviewed technique, imaging mass cytometry (153), to interrogate the TME of lung adenocarcinoma and identify spatial features from 5 µm sections that were predictive of recurrence (154). This group used the same technique to reveal cellular neighbourhoods that were associated with survival in glioblastoma, which were then used to identify a specific population of macrophages that were associated with long-term survival (155). Antibody-based techniques are a classical method to interrogate the spatial organization of the TME, and technological advances are facilitating highly-multiplexed applications to improve our understanding of this organization.
4.2 Oligonucleotide probe-based methods
While antibody-based protein analysis is suitable for many applications, antibodies against some proteins are challenging to create, and achieving whole-proteome multiplexing is currently unrealistic. Interrogating RNA expression instead, such as through the use of oligonucleotide probes, can circumvent these difficulties and streamline the investigation of transcript-level variations (141).
Conjugating the ‘indexing’ oligonucleotides used in DSP to single-stranded RNA probes instead of antibodies allows DSP to be used for the analysis of mRNAs (141). This method has achieved over 2000-fold mRNA multiplexing, but the detection of above-background signal for a given transcript requires around 600 copies of that transcript, and consequently most mRNAs can only be detected in relatively large (≥ 50 µm diameter) ROIs (140, 156).
seqFISH+ employs primary oligonucleotide probes that each contain a unique four-region barcode, with each region being complementary to a single ‘readout probe’ sequence (142). Imaging consists of four rounds, each with 20 hybridize-image-strip cycles in which three distinct readout probes, each conjugated to a different fluorescent dye, are hybridized to regions of the primary probes (142). Each primary probe’s barcode can be uniquely identified by the set of three readout probes that hybridize with it, one during each imaging round, with the final round used for correcting errors. In a given imaging cycle only 1 in 60 mRNAs is visualized, which minimizes optical crowding and enables the sub-diffraction limit localization of each mRNA molecule (142). seqFISH+ has demonstrated 10,000-fold multiplexing, with a sensitivity greater than that of single-cell RNA sequencing (scRNA-seq), but like most cyclic or sequential imaging methods, it requires lengthy workflow times (142).
In spatially-resolved transcript amplicon readout mapping (STARmap), DNA amplicons are constructed in hydrogel-embedded tissue through reverse transcription (RT), cDNA circularization, and rolling circle amplification (RCA), with primer and padlock probes used to prevent non-specific amplification (143). Each padlock probe is specific to a single gene and contains a unique five-base barcode, which can later be decoded by sequencing with error-reduction by dynamic annealing and ligation (SEDAL) (143). The degree of multiplexing depends on barcode length, and thus far 1020-fold multiplexing has been demonstrated at single-cell resolution, but high degrees of multiplexing limit STARmap’s sensitivity, due to the difficulty of optically resolving amplicons that are located in physical proximity (143). Under ideal conditions, STARmap’s sensitivity is comparable to that of scRNA-seq.
Oligonucleotide probes can also be used for the identification of bacterial taxa, through methods such as HiPR-FISH. In HiPR-FISH, several probes per interrogated taxon, each specific to the same unique sequence within that taxon’s rRNA, are first introduced to the sample (144). Each probe carries two flanking sequences, each complementary to one of ten fluorescent readout probes, chosen in such a way that each taxon’s probes collectively contain sequences complementary to a unique subset of the readout probes. After hybridization of the readout probes and fluorescence imaging, cells are segmented and then a machine learning classifier is used to assign each cell to a taxon based on its emission spectra (144). HiPR-FISH can potentially differentiate between 1023 distinct taxa in only one round of imaging, but notably it cannot detect unexpected or rare taxa for which probes were not designed (144).
4.3 Sequencing-based methods
The second major approach to spatially quantifying mRNA expression is to employ sequencing, either in situ or after spatially barcoding transcripts, which enables whole-transcriptome multiplexing.
With fluorescent in situ sequencing (FISSEQ), sequencing libraries are constructed in situ through a process of RT, cDNA circularization, and RCA (145). Cross-linking of aminoallyl dUTP residues introduced during RT prevents cDNA diffusion. Each amplicon is sequenced by oligonucleotide ligation and detection (SOLiD): sequencing primers are annealed, and fluorescently-tagged oligonucleotides are used to identify every fifth nucleotide (145). The primer is then stripped and the procedure is repeated with primers of incrementally shorter lengths, so that all nucleotides in the sequence are eventually ascertained. FISSEQ operates at single-cell resolutions, but with low sensitivity, as just ~200 mRNA reads are acquired from each cell (145). FISSEQ is also relatively slow, as the SOLiD sequencing step requires around 10 days on the microscope.
In spatial transcriptomics (ST), mRNAs are captured from tissue by oligo(dT)-containing probes that have been affixed to a slide in an array of 55-100 µm diameter spots (146). The probes in each spot carry a unique barcode, so that after all captured mRNAs are pooled and sequenced, individual reads can be assigned to their spot of origin. This strategy achieves whole-transcriptome multiplexing at levels of sensitivity comparable to standard next generation RNA sequencing techniques, but at ~10-cell resolutions (146). High-definition spatial transcriptomics (HD-ST) is a variant of ST that achieves significantly higher spatial resolution by affixing barcoded probes to 2 µm diameter beads, which reside in individual wells on a slide (147). HD-ST maintains the whole-transcriptome multiplexing capacity of ST, and operates at resolutions much closer to those of hybridization-based imaging techniques, but the decreased area of tissue that corresponds to each well greatly reduces its RNA capture efficiency and consequently sensitivity (147).
Slide-seq is a similar sequencing-based method that also relies upon barcoded oligo(dT)-containing probes (157). Barcoded probes are affixed to 10 µm diameter beads, which are deposited on a glass coverslip, and SOLiD is used to map each barcode to the corresponding bead (157). The tissue being studied is then transferred onto the coverslip for mRNA hybridization, after which the captured mRNA can be pooled and sequenced (157). The strengths and weaknesses of Slide-seq are similar to those of HD-ST: it achieves whole-transcriptome multiplexing at a moderate resolution, but with relatively low sensitivity (roughly 5% of that of scRNA-seq) (157). Slide-seq has been followed by Slide-seqV2, which replaces SOLiD encoding with a barcoding scheme that is more robust to errors and includes an additional second-strand synthesis step between the reverse transcription of captured mRNA and PCR, making Slide-seqV2 an order of magnitude more sensitive (148).
Successful applications of these RNA-interrogating techniques include investigating intra-patient heterogeneity in the expression of genes associated with neuroendocrine and androgen receptor activity in metastatic prostate cancer (156), characterizing the specific cell subpopulations and phenotypes present at the leading edge of squamous cell carcinomas (158), and identifying three recurring types of spatial cellular communities within PDAC tumours, each with a unique composition of malignant, stromal, and immune cell subtypes (159). Another recent application was the use of Slide-seqV2 to describe transcriptomic alterations induced by neighbouring cells in the immune-suppressive TME that promote tumourigenesis of prostate cancer (160). The authors further combined the spatial data obtained via Slide-seqV2 with computational approaches to ligand-receptor pair identification, which enabled them to identify candidate ligand-receptor pairs that were specifically expressed in neighbouring cells, and hence uncover a specific axis that contributes to prostate TME immunosuppression (160).
5 Conclusions
An increasing appreciation for the tremendous inter- and intra-tumor functional and phenotypic heterogeneity that exists in nearly every cancer type, both at the tumour and TME levels, underscores the importance of actualizing the goals of precision medicine for individual cancer patients. Achieving these goals will require expanded access to molecular profiling for more patient tumours and a more complete understanding of what these data mean in the context of spatially heterogeneous TMEs, including how they impact therapeutic response. Historically, the widespread use of mIF and mIHC to interrogate the TME has revealed how the abundance and degree of infiltration of immune and stromal cells impact cancer phenotypes. The application of increasingly sophisticated computational methods to mIF and mIHC data has shown that non-cell autonomous interactions within the TME, driven either through direct cell-cell contacts or indirect mechanisms such as metabolite secretion, greatly influence cellular functions and ultimately tumour behaviour. These advances have already led to the development of promising novel treatment modalities, including immune checkpoint blockade therapy. Recently developed profiling technologies offer the multiplexing capacity and the spatial resolution that are needed to further expand our understanding of tumour biology and therapeutic response. Of these technologies, those that are operable without prohibitive amounts of expertise and that can be made available to many patients by virtue of low cost or high sample throughput hold the most potential for effective translation into the clinic. A combination of increased use of these novel technologies and computational advances that improve the interpretation of multiplexed imaging data will give researchers and clinicians the opportunity to develop and apply treatment protocols tailored to the unique dynamics of every patient’s TME.
Author contributions
DC: Conceptualization, Visualization, Writing – original draft, Writing – review & editing. AF: Writing – original draft, Writing – review & editing. EM: Conceptualization, Visualization, Writing – original draft, Writing – review & editing. EV: Conceptualization, Writing – original draft, Writing – review & editing. GS: Writing – original draft, Writing – review & editing. KN: Writing – original draft, Writing – review & editing. WWL: Writing – review & editing. CM: Writing – original draft, Writing – review & editing. MG: Writing – original draft, Writing – review & editing. WLL: Conceptualization, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by funds from the Canadian Institutes of Health Research (CIHR FRN-143345 and 183775), the Terry Fox Foundation, the Lotte and John Hecht Memorial Foundation, and the BC Cancer Foundation. DC, AF, EM, EV and KN were supported by the CIHR Frederick Banting and Charles Best Canada Graduate Scholarships and University of British Columbia fellowships. EM was a Vanier Canada Graduate Scholar.
Acknowledgments
The authors would like to thank Dr. Kim Lonergan for her expert critique of this manuscript. All figures were created with BioRender.com. Figure 3 raw multi-channel image reproduced from (31) with permission of the authors under the terms of the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0/).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2023.1275890/full#supplementary-material
Glossary
References
1. Kartikasari AER, Huertas CS, Mitchell A, Plebanski M. Tumor-induced inflammatory cytokines and the emerging diagnostic devices for cancer detection and prognosis. Front Oncol (2021) 11:692142. doi: 10.3389/fonc.2021.692142
2. Jing X, Yang F, Shao C, Wei K, Xie M, Shen H, et al. Role of hypoxia in cancer therapy by regulating the tumor microenvironment. Mol cancer (2019) 18(1):157. doi: 10.1186/s12943-019-1089-9
3. Corbet C, Feron O. Tumour acidosis: from the passenger to the driver’s seat. Nat Rev Cancer (2017) 17(10):577–93. doi: 10.1038/nrc.2017.77
4. Riquelme E, Zhang Y, Zhang L, Montiel M, Zoltan M, Dong W, et al. Tumor microbiome diversity and composition influence pancreatic cancer outcomes. Cell (2019) 178(4):795–806 e12. doi: 10.1016/j.cell.2019.07.008
5. Nejman D, Livyatan I, Fuks G, Gavert N, Zwang Y, Geller LT, et al. The human tumor microbiome is composed of tumor type-specific intracellular bacteria. Science (2020) 368(6494):973–80. doi: 10.1126/science.aay9189
6. Winkler J, Abisoye-Ogunniyan A, Metcalf KJ, Werb Z. Concepts of extracellular matrix remodelling in tumour progression and metastasis. Nat Commun (2020) 11(1):5120. doi: 10.1038/s41467-020-18794-x
7. Kai F, Drain AP, Weaver VM. The extracellular matrix modulates the metastatic journey. Dev Cell (2019) 49(3):332–46. doi: 10.1016/j.devcel.2019.03.026
8. DeNardo DG, Ruffell B. Macrophages as regulators of tumour immunity and immunotherapy. Nat Rev Immunol (2019) 19(6):369–82. doi: 10.1038/s41577-019-0127-6
9. Herrera M, Herrera A, Dominguez G, Silva J, Garcia V, Garcia JM, et al. Cancer-associated fibroblast and M2 macrophage markers together predict outcome in colorectal cancer patients. Cancer science (2013) 104(4):437–44. doi: 10.1111/cas.12096
10. Brechbuhl HM, Finlay-Schultz J, Yamamoto TM, Gillen AE, Cittelly DM, Tan AC, et al. Fibroblast subtypes regulate responsiveness of luminal breast cancer to estrogen. Clin Cancer Res (2017) 23(7):1710–21. doi: 10.1158/1078-0432.CCR-15-2851
11. Yang L, Roberts D, Takhar M, Erho N, Bibby BAS, Thiruthaneeswaran N, et al. Development and validation of a 28-gene hypoxia-related prognostic signature for localized prostate cancer. EBioMedicine (2018) 31:182–9. doi: 10.1016/j.ebiom.2018.04.019
12. Stevanović S, Pasetto A, Helman SR, Gartner JJ, Prickett TD, Howie B, et al. Landscape of immunogenic tumor antigens in successful immunotherapy of virally induced epithelial cancer. Science (2017) 356(6334):200–5. doi: 10.1126/science.aak9510
13. Zhao Q, Laverdure JP, Lanoix J, Durette C, Côté C, Bonneil É, et al. Proteogenomics uncovers a vast repertoire of shared tumor-specific antigens in ovarian cancer. Cancer Immunol Res (2020) 8(4):544–55. doi: 10.1158/2326-6066.CIR-19-0541
14. Balachandran VP, Łuksza M, Zhao JN, Makarov V, Moral JA, Remark R, et al. Identification of unique neoantigen qualities in long-term survivors of pancreatic cancer. Nature (2017) 551(7681):512–6. doi: 10.1038/nature24462
15. Pages F, Mlecnik B, Marliot F, Bindea G, Ou FS, Bifulco C, et al. International validation of the consensus Immunoscore for the classification of colon cancer: a prognostic and accuracy study. Lancet (London England) (2018) 391(10135):2128–39. doi: 10.1016/S0140-6736(18)30789-X
16. Matsumoto H, Thike AA, Li H, Yeong J, Koo SL, Dent RA, et al. Increased CD4 and CD8-positive T cell infiltrate signifies good prognosis in a subset of triple-negative breast cancer. Breast Cancer Res Treat (2016) 156(2):237–47. doi: 10.1007/s10549-016-3743-x
17. Bottcher JP, Bonavita E, Chakravarty P, Blees H, Cabeza-Cabrerizo M, Sammicheli S, et al. NK Cells Stimulate Recruitment of cDC1 into the Tumor Microenvironment Promoting Cancer Immune Control. Cell (2018) 172(5):1022–37 e14. doi: 10.1016/j.cell.2018.01.004
18. Berntsson J, Nodin B, Eberhard J, Micke P, Jirstrom K. Prognostic impact of tumour-infiltrating B cells and plasma cells in colorectal cancer. Int J cancer (2016) 139(5):1129–39. doi: 10.1002/ijc.30138
19. Shang B, Liu Y, Jiang SJ, Liu Y. Prognostic value of tumor-infiltrating FoxP3+ regulatory T cells in cancers: a systematic review and meta-analysis. Sci Rep (2015) 5:15179. doi: 10.1038/srep15179
20. Ai L, Mu S, Wang Y, Wang H, Cai L, Li W, et al. Prognostic role of myeloid-derived suppressor cells in cancers: a systematic review and meta-analysis. BMC cancer (2018) 18(1):1220. doi: 10.1186/s12885-018-5086-y
21. Zhang QW, Liu L, Gong CY, Shi HS, Zeng YH, Wang XZ, et al. Prognostic significance of tumor-associated macrophages in solid tumor: a meta-analysis of the literature. PloS One (2012) 7(12):e50946. doi: 10.1371/journal.pone.0050946
22. Semeraro M, Adam J, Stoll G, Louvet E, Chaba K, Poirier-Colame V, et al. The ratio of CD8(+)/FOXP3 T lymphocytes infiltrating breast tissues predicts the relapse of ductal carcinoma in situ. Oncoimmunology (2016) 5(10):e1218106. doi: 10.1080/2162402X.2016.1218106
23. Chaves ALF, Silva AG, Maia FM, Lopes GFM, de Paulo LFB, Muniz LV, et al. Reduced CD8(+) T cells infiltration can be associated to a Malignant transformation in potentially Malignant oral epithelial lesions. Clin Oral investigations (2019) 23(4):1913–9. doi: 10.1007/s00784-018-2622-8
24. Chen DS, Mellman I. Elements of cancer immunity and the cancer-immune set point. Nature (2017) 541(7637):321–30. doi: 10.1038/nature21349
25. Tumeh PC, Harview CL, Yearley JH, Shintaku IP, Taylor EJ, Robert L, et al. PD-1 blockade induces responses by inhibiting adaptive immune resistance. Nature (2014) 515(7528):568–71. doi: 10.1038/nature13954
26. Denkert C, von Minckwitz G, Darb-Esfahani S, Lederer B, Heppner BI, Weber KE, et al. Tumour-infiltrating lymphocytes and prognosis in different subtypes of breast cancer: a pooled analysis of 3771 patients treated with neoadjuvant therapy. Lancet Oncol (2018) 19(1):40–50. doi: 10.1016/S1470-2045(17)30904-X
27. Olive KP, Jacobetz MA, Davidson CJ, Gopinathan A, McIntyre D, Honess D, et al. Inhibition of Hedgehog signaling enhances delivery of chemotherapy in a mouse model of pancreatic cancer. Science (2009) 324(5933):1457–61. doi: 10.1126/science.1171362
28. Herbst RS, Soria JC, Kowanetz M, Fine GD, Hamid O, Gordon MS, et al. Predictive correlates of response to the anti-PD-L1 antibody MPDL3280A in cancer patients. Nature (2014) 515(7528):563–7. doi: 10.1038/nature14011
29. Bocchialini G, Lagrasta C, Madeddu D, Mazzaschi G, Marturano D, Sogni F, et al. Spatial architecture of tumour-infiltrating lymphocytes as a prognostic parameter in resected non-small-cell lung cancer. Eur J cardio-thoracic Surg (2020) 58(3):619–28. doi: 10.1093/ejcts/ezaa098
30. Taube JM, Roman K, Engle EL, Wang C, Ballesteros-Merino C, Jensen SM, et al. Multi-institutional TSA-amplified multiplexed immunofluorescence reproducibility evaluation (MITRE) study. J Immunother Cancer (2021) 9(7):e002197. doi: 10.1136/jitc-2020-002197
31. Enfield KSS, Martin SD, Marshall EA, Kung SHY, Gallagher P, Milne K, et al. Hyperspectral cell sociology reveals spatial tumor-immune cell interactions associated with lung cancer recurrence. J Immunother Cancer (2019) 7(1):13. doi: 10.1186/s40425-018-0488-6
32. Lazarus J, Maj T, Smith JJ, Perusina Lanfranca M, Rao A, D’Angelica MI, et al. Spatial and phenotypic immune profiling of metastatic colon cancer. JCI Insight (2018) 3(22):e121932. doi: 10.1172/jci.insight.121932
33. Egelston CA, Avalos C, Tu TY, Rosario A, Wang R, Solomon S, et al. Resident memory CD8+ T cells within cancer islands mediate survival in breast cancer patients. JCI Insight (2019) 4(19):e130000. doi: 10.1172/jci.insight.130000
34. Hladikova K, Koucky V, Boucek J, Laco J, Grega M, Hodek M, et al. Tumor-infiltrating B cells affect the progression of oropharyngeal squamous cell carcinoma via cell-to-cell interactions with CD8(+) T cells. J Immunother Cancer (2019) 7(1):261. doi: 10.1186/s40425-019-0726-6
35. Sautes-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. doi: 10.1038/s41568-019-0144-6
36. Chaurio RA, Anadon CM, Lee Costich T, Payne KK, Biswas S, Harro CM, et al. TGF-β-mediated silencing of genomic organizer SATB1 promotes Tfh cell differentiation and formation of intra-tumoral tertiary lymphoid structures. Immunity (2022) 55(1):115–28.e9. doi: 10.1016/j.immuni.2021.12.007
37. Merkenschlager J, Finkin S, Ramos V, Kraft J, Cipolla M, Nowosad CR, et al. Dynamic regulation of T(FH) selection during the germinal centre reaction. Nature (2021) 591(7850):458–63. doi: 10.1038/s41586-021-03187-x
38. Lin X, Ye L, Wang X, Liao Z, Dong J, Yang Y, et al. Follicular helper T cells remodel the immune microenvironment of pancreatic cancer via secreting CXCL13 and IL-21. Cancers (2021) 13(15):3678. doi: 10.3390/cancers13153678
39. Cui C, Wang J, Fagerberg E, Chen PM, Connolly KA, Damo M, et al. Neoantigen-driven B cell and CD4 T follicular helper cell collaboration promotes anti-tumor CD8 T cell responses. Cell (2021) 184(25):6101–18.e13. doi: 10.1016/j.cell.2021.11.007
40. Cabrita R, Lauss M, Sanna A, Donia M, Skaarup Larsen M, Mitra S, et al. Tertiary lymphoid structures improve immunotherapy and survival in melanoma. Nature (2020) 577(7791):561–5. doi: 10.1038/s41586-019-1914-8
41. Ruffin AT, Cillo AR, Tabib T, Liu A, Onkar S, Kunning SR, et al. B cell signatures and tertiary lymphoid structures contribute to outcome in head and neck squamous cell carcinoma. Nat Commun (2021) 12(1):3349. doi: 10.1038/s41467-021-23355-x
42. Castino GF, Cortese N, Capretti G, Serio S, Di Caro G, Mineri R, et al. Spatial distribution of B cells predicts prognosis in human pancreatic adenocarcinoma. Oncoimmunology (2016) 5(4):e1085147. doi: 10.1080/2162402X.2015.1085147
43. Pfannstiel C, Strissel PL, Chiappinelli KB, Sikic D, Wach S, Wirtz RM, et al. The tumor immune microenvironment drives a prognostic relevance that correlates with bladder cancer subtypes. Cancer Immunol Res (2019) 7(6):923–38. doi: 10.1158/2326-6066.CIR-18-0758
44. Goc J, Germain C, Vo-Bourgais TK, Lupo A, Klein C, Knockaert S, et al. Dendritic cells in tumor-associated tertiary lymphoid structures signal a Th1 cytotoxic immune contexture and license the positive prognostic value of infiltrating CD8+ T cells. Cancer Res (2014) 74(3):705–15. doi: 10.1158/0008-5472.CAN-13-1342
45. Behr DS, Peitsch WK, Hametner C, Lasitschka F, Houben R, Schonhaar K, et al. Prognostic value of immune cell infiltration, tertiary lymphoid structures and PD-L1 expression in Merkel cell carcinomas. Int J Clin Exp pathol (2014) 7(11):7610–21.
46. Horeweg N, Workel HH, Loiero D, Church DN, Vermij L, Léon-Castillo A, et al. Tertiary lymphoid structures critical for prognosis in endometrial cancer patients. Nat Commun (2022) 13(1):1373. doi: 10.1038/s41467-022-29040-x
47. Petitprez F, de Reynies A, Keung EZ, Chen TW, Sun CM, Calderaro J, et al. B cells are associated with survival and immunotherapy response in sarcoma. Nature (2020) 577(7791):556–60. doi: 10.1038/s41586-019-1906-8
48. Hiraoka N, Ino Y, Yamazaki-Itoh R, Kanai Y, Kosuge T, Shimada K. Intratumoral tertiary lymphoid organ is a favourable prognosticator in patients with pancreatic cancer. Br J Cancer (2015) 112(11):1782–90. doi: 10.1038/bjc.2015.145
49. Germain C, Devi-Marulkar P, Knockaert S, Biton J, Kaplon H, Letaief L, et al. Tertiary lymphoid structure-B cells narrow regulatory T cells impact in lung cancer patients. Front Immunol (2021) 12:626776. doi: 10.3389/fimmu.2021.626776
50. Di Caro G, Bergomas F, Grizzi F, Doni A, Bianchi P, Malesci A, et al. Occurrence of tertiary lymphoid tissue is associated with T-cell infiltration and predicts better prognosis in early-stage colorectal cancers. Clin Cancer Res (2014) 20(8):2147–58. doi: 10.1158/1078-0432.CCR-13-2590
51. Meylan M, Petitprez F, Becht E, Bougouin A, Pupier G, Calvez A, et al. Tertiary lymphoid structures generate and propagate anti-tumor antibody-producing plasma cells in renal cell cancer. Immunity (2022) 55(3):527–41 e5. doi: 10.1016/j.immuni.2022.02.001
52. Ng KW, Marshall EA, Enfield KS, Martin SD, Milne K, Pewarchuk ME, et al. Somatic mutation-associated T follicular helper cell elevation in lung adenocarcinoma. Oncoimmunology (2018) 7(12):e1504728. doi: 10.1080/2162402X.2018.1504728
53. Shang T, Jiang T, Lu T, Wang H, Cui X, Pan Y, et al. Tertiary lymphoid structures predict the prognosis and immunotherapy response of cholangiocarcinoma. Front Immunol (2023) 14:1166497. doi: 10.3389/fimmu.2023.1166497
54. Brennan CA, Garrett WS. Fusobacterium nucleatum - symbiont, opportunist and oncobacterium. Nat Rev Microbiol (2019) 17(3):156–66. doi: 10.1038/s41579-018-0129-6
55. Gur C, Ibrahim Y, Isaacson B, Yamin R, Abed J, Gamliel M, et al. Binding of the Fap2 protein of Fusobacterium nucleatum to human inhibitory receptor TIGIT protects tumors from immune cell attack. Immunity (2015) 42(2):344–55. doi: 10.1016/j.immuni.2015.01.010
56. Rubinstein MR, Wang X, Liu W, Hao Y, Cai G, Han YW. Fusobacterium nucleatum promotes colorectal carcinogenesis by modulating E-cadherin/beta-catenin signaling via its FadA adhesin. Cell Host Microbe (2013) 14(2):195–206. doi: 10.1016/j.chom.2013.07.012
57. Yang Y, Weng W, Peng J, Hong L, Yang L, Toiyama Y, et al. Fusobacterium nucleatum increases proliferation of colorectal cancer cells and tumor development in mice by activating toll-like receptor 4 signaling to nuclear factor-kappaB, and up-regulating expression of microRNA-21. Gastroenterology (2017) 152(4):851–66 e24. doi: 10.1053/j.gastro.2016.11.018
58. Yu T, Guo F, Yu Y, Sun T, Ma D, Han J, et al. Fusobacterium nucleatum promotes chemoresistance to colorectal cancer by modulating autophagy. Cell (2017) 170(3):548–63 e16. doi: 10.1016/j.cell.2017.07.008
59. Gong J, Chehrazi-Raffle A, Reddi S, Salgia R. Development of PD-1 and PD-L1 inhibitors as a form of cancer immunotherapy: a comprehensive review of registration trials and future considerations. J Immunother Cancer (2018) 6(1):8. doi: 10.1186/s40425-018-0316-z
60. Ribas A, Hu-Lieskovan S. What does PD-L1 positive or negative mean? J Exp Med (2016) 213(13):2835–40. doi: 10.1084/jem.20161462
61. Lu S, Stein JE, Rimm DL, Wang DW, Bell JM, Johnson DB, et al. Comparison of biomarker modalities for predicting response to PD-1/PD-L1 checkpoint blockade: A systematic review and meta-analysis. JAMA Oncol (2019) 5(8):1195–204. doi: 10.1001/jamaoncol.2019.1549
62. Johnson DB, Bordeaux J, Kim JY, Vaupel C, Rimm DL, Ho TH, et al. Quantitative spatial profiling of PD-1/PD-L1 interaction and HLA-DR/IDO-1 predicts improved outcomes of anti-PD-1 therapies in metastatic melanoma. Clin Cancer Res (2018) 24(21):5250–60. doi: 10.1158/1078-0432.CCR-18-0309
63. Giraldo NA, Nguyen P, Engle EL, Kaunitz GJ, Cottrell TR, Berry S, et al. Multidimensional, quantitative assessment of PD-1/PD-L1 expression in patients with Merkel cell carcinoma and association with response to pembrolizumab. J Immunother Cancer (2018) 6(1):99. doi: 10.1186/s40425-018-0404-0
64. Farkkila A, Gulhan DC, Casado J, Jacobson CA, Nguyen H, Kochupurakkal B, et al. Immunogenomic profiling determines responses to combined PARP and PD-1 inhibition in ovarian cancer. Nat Commun (2020) 11(1):1459. doi: 10.1038/s41467-020-15315-8
65. Gide TN, Silva IP, Quek C, Ahmed T, Menzies AM, Carlino MS, et al. Close proximity of immune and tumor cells underlies response to anti-PD-1 based therapies in metastatic melanoma patients. Oncoimmunology (2020) 9(1):1659093. doi: 10.1080/2162402X.2019.1659093
66. Antoranz A, Van Herck Y, Bolognesi MM, Lynch SM, Rahman A, Gallagher WM, et al. Mapping the immune landscape in metastatic melanoma reveals localized cell-cell interactions that predict immunotherapy response. Cancer Res (2022) 82(18):3275–90. doi: 10.1158/0008-5472.CAN-22-0363
67. Tsakiroglou AM, Fergie M, Oguejiofor K, Linton K, Thomson D, Stern PL, et al. Spatial proximity between T and PD-L1 expressing cells as a prognostic biomarker for oropharyngeal squamous cell carcinoma. Br J Cancer (2020) 122(4):539–44. doi: 10.1038/s41416-019-0634-z
68. Feng Z, Bethmann D, Kappler M, Ballesteros-Merino C, Eckert A, Bell RB, et al. Multiparametric immune profiling in HPV- oral squamous cell cancer. JCI Insight (2017) 2(14):e93652. doi: 10.1172/jci.insight.93652
69. Li L, Sun R, Miao Y, Tran T, Adams L, Roscoe N, et al. PD-1/PD-L1 expression and interaction by automated quantitative immunofluorescent analysis show adverse prognostic impact in patients with diffuse large B-cell lymphoma having T-cell infiltration: a study from the International DLBCL Consortium Program. Modern Pathol (2019) 32(6):741–54. doi: 10.1038/s41379-018-0193-5
70. McNally A, Hill GR, Sparwasser T, Thomas R, Steptoe RJ. CD4+CD25+ regulatory T cells control CD8+ T-cell effector differentiation by modulating IL-2 homeostasis. Proc Natl Acad Sci United States America (2011) 108(18):7529–34. doi: 10.1073/pnas.1103782108
71. Jang JE, Hajdu CH, Liot C, Miller G, Dustin ML, Bar-Sagi D. Crosstalk between regulatory T cells and tumor-associated dendritic cells negates anti-tumor immunity in pancreatic cancer. Cell Rep (2017) 20(3):558–71. doi: 10.1016/j.celrep.2017.06.062
72. Strauss L, Bergmann C, Whiteside TL. Human circulating CD4+CD25highFoxp3+ regulatory T cells kill autologous CD8+ but not CD4+ responder cells by Fas-mediated apoptosis. J Immunol (Baltimore Md: 1950) (2009) 182(3):1469–80. doi: 10.4049/jimmunol.182.3.1469
73. Cao X, Cai SF, Fehniger TA, Song J, Collins LI, Piwnica-Worms DR, et al. Granzyme B and perforin are important for regulatory T cell-mediated suppression of tumor clearance. Immunity (2007) 27(4):635–46. doi: 10.1016/j.immuni.2007.08.014
74. Yamaguchi T, Wing JB, Sakaguchi S. Two modes of immune suppression by Foxp3(+) regulatory T cells under inflammatory or non-inflammatory conditions. Semin Immunol (2011) 23(6):424–30. doi: 10.1016/j.smim.2011.10.002
75. Bergsland CH, Jeanmougin M, Moosavi SH, Svindland A, Bruun J, Nesbakken A, et al. Spatial analysis and CD25-expression identify regulatory T cells as predictors of a poor prognosis in colorectal cancer. Modern Pathol (2022) 35(9):1236–46. doi: 10.1038/s41379-022-01086-8
76. Barua S, Fang P, Sharma A, Fujimoto J, Wistuba I, Rao AUK, et al. Spatial interaction of tumor cells and regulatory T cells correlates with survival in non-small cell lung cancer. Lung cancer (2018) 117:73–9. doi: 10.1016/j.lungcan.2018.01.022
77. Mezheyeuski A, Bergsland CH, Backman M, Djureinovic D, Sjoblom T, Bruun J, et al. Multispectral imaging for quantitative and compartment-specific immune infiltrates reveals distinct immune profiles that classify lung cancer patients. J pathol (2018) 244(4):421–31. doi: 10.1002/path.5026
78. Feichtenbeiner A, Haas M, Buttner M, Grabenbauer GG, Fietkau R, Distel LV. Critical role of spatial interaction between CD8(+) and Foxp3(+) cells in human gastric cancer: the distance matters. Cancer immunol immunother: CII (2014) 63(2):111–9. doi: 10.1007/s00262-013-1491-x
79. Joshi NS, Akama-Garren EH, Lu Y, Lee DY, Chang GP, Li A, et al. Regulatory T cells in tumor-associated tertiary lymphoid structures suppress anti-tumor T cell responses. Immunity (2015) 43(3):579–90. doi: 10.1016/j.immuni.2015.08.006
80. Gobert M, Treilleux I, Bendriss-Vermare N, Bachelot T, Goddard-Leon S, Arfi V, et al. Regulatory T cells recruited through CCL22/CCR4 are selectively activated in lymphoid infiltrates surrounding primary breast tumors and lead to an adverse clinical outcome. Cancer Res (2009) 69(5):2000–9. doi: 10.1158/0008-5472.CAN-08-2360
81. Chen Y, Song Y, Du W, Gong L, Chang H, Zou Z. Tumor-associated macrophages: an accomplice in solid tumor progression. J Biomed science (2019) 26(1):78. doi: 10.1186/s12929-019-0568-z
82. Zheng X, Weigert A, Reu S, Guenther S, Mansouri S, Bassaly B, et al. Spatial density and distribution of tumor-associated macrophages predict survival in non-small cell lung carcinoma. Cancer Res (2020) 80(20):4414–25. doi: 10.1158/0008-5472.CAN-20-0069
83. Vayrynen SA, Zhang J, Yuan C, Vayrynen JP, Dias Costa A, Williams H, et al. Composition, spatial characteristics, and prognostic significance of myeloid cell infiltration in pancreatic cancer. Clin Cancer Res (2021) 27(4):1069–81. doi: 10.1158/1078-0432.CCR-20-3141
84. Huang YK, Wang M, Sun Y, Di Costanzo N, Mitchell C, Achuthan A, et al. Macrophage spatial heterogeneity in gastric cancer defined by multiplex immunohistochemistry. Nat Commun (2019) 10(1):3928. doi: 10.1038/s41467-019-11788-4
85. Kersten K, Hu KH, Combes AJ, Samad B, Harwin T, Ray A, et al. Spatiotemporal co-dependency between macrophages and exhausted CD8(+) T cells in cancer. Cancer Cell (2022) 40(6):624–38 e9. doi: 10.1016/j.ccell.2022.05.004
86. Gartrell RD, Marks DK, Hart TD, Li G, Davari DR, Wu A, et al. Quantitative analysis of immune infiltrates in primary melanoma. Cancer Immunol Res (2018) 6(4):481–93. doi: 10.1158/2326-6066.CIR-17-0360
87. Uyttenhove C, Pilotte L, Theate I, Stroobant V, Colau D, Parmentier N, et al. Evidence for a tumoral immune resistance mechanism based on tryptophan degradation by indoleamine 2,3-dioxygenase. Nat Med (2003) 9(10):1269–74. doi: 10.1038/nm934
88. Opitz CA, Litzenburger UM, Sahm F, Ott M, Tritschler I, Trump S, et al. An endogenous tumour-promoting ligand of the human aryl hydrocarbon receptor. Nature (2011) 478(7368):197–203. doi: 10.1038/nature10491
89. Mezrich JD, Fechner JH, Zhang X, Johnson BP, Burlingham WJ, Bradfield CA. An interaction between kynurenine and the aryl hydrocarbon receptor can generate regulatory T cells. J Immunol (Baltimore Md: 1950) (2010) 185(6):3190–8. doi: 10.4049/jimmunol.0903670
90. Greene LI, Bruno TC, Christenson JL, D’Alessandro A, Culp-Hill R, Torkko K, et al. A role for tryptophan-2,3-dioxygenase in CD8 T-cell suppression and evidence of tryptophan catabolism in breast cancer patient plasma. Mol Cancer research: MCR (2019) 17(1):131–9. doi: 10.1158/1541-7786.MCR-18-0362
91. Bishnupuri KS, Alvarado DM, Khouri AN, Shabsovich M, Chen B, Dieckgraefe BK, et al. IDO1 and kynurenine pathway metabolites activate PI3K-akt signaling in the neoplastic colon epithelium to promote cancer cell proliferation and inhibit apoptosis. Cancer Res (2019) 79(6):1138–50. doi: 10.1158/0008-5472.CAN-18-0668
92. van Baren N, Van den Eynde BJ. Tryptophan-degrading enzymes in tumoral immune resistance. Front Immunol (2015) 6:34. doi: 10.3389/fimmu.2015.00034
93. Long GV, Dummer R, Hamid O, Gajewski TF, Caglevic C, Dalle S, et al. Epacadostat plus pembrolizumab versus placebo plus pembrolizumab in patients with unresectable or metastatic melanoma (ECHO-301/KEYNOTE-252): a phase 3, randomised, double-blind study. Lancet Oncol (2019) 20(8):1083–97. doi: 10.1016/S1470-2045(19)30274-8
94. Van den Eynde BJ, van Baren N, Baurain J-F. Is there a clinical future for IDO1 inhibitors after the failure of epacadostat in melanoma? Annu Rev Cancer Biol (2020) 4(1):241–56. doi: 10.1146/annurev-cancerbio-030419-033635
95. Geller LT, Barzily-Rokni M, Danino T, Jonas OH, Shental N, Nejman D, et al. Potential role of intratumor bacteria in mediating tumor resistance to the chemotherapeutic drug gemcitabine. Science (2017) 357(6356):1156–60. doi: 10.1126/science.aah5043
96. Han R, Sun Q, Wu J, Zheng P, Zhao G. Sodium butyrate upregulates miR-203 expression to exert anti-proliferation effect on colorectal cancer cells. Cell Physiol biochem: Int J Exp Cell physiol biochem Pharmacol (2016) 39(5):1919–29. doi: 10.1159/000447889
97. Zeng H, Taussig DP, Cheng WH, Johnson LK, Hakkak R. Butyrate inhibits cancerous HCT116 colon cell proliferation but to a lesser extent in noncancerous NCM460 colon cells. Nutrients (2017) 9(1):25. doi: 10.3390/nu9010025
98. Li Q, Ding C, Meng T, Lu W, Liu W, Hao H, et al. Butyrate suppresses motility of colorectal cancer cells via deactivating Akt/ERK signaling in histone deacetylase dependent manner. J Pharmacol Sci (2017) 135(4):148–55. doi: 10.1016/j.jphs.2017.11.004
99. Al Tameemi W, Dale TP, Al-Jumaily RMK, Forsyth NR. Hypoxia-modified cancer cell metabolism. Front Cell Dev Biol (2019) 7:4. doi: 10.3389/fcell.2019.00004
100. Ullah MS, Davies AJ, Halestrap AP. The plasma membrane lactate transporter MCT4, but not MCT1, is up-regulated by hypoxia through a HIF-1alpha-dependent mechanism. J Biol Chem (2006) 281(14):9030–7. doi: 10.1074/jbc.M511397200
101. Meijer TW, Schuurbiers OC, Kaanders JH, Looijen-Salamon MG, de Geus-Oei LF, Verhagen AF, et al. Differences in metabolism between adeno- and squamous cell non-small cell lung carcinomas: spatial distribution and prognostic value of GLUT1 and MCT4. Lung cancer (2012) 76(3):316–23. doi: 10.1016/j.lungcan.2011.11.006
102. Sonveaux P, Vegran F, Schroeder T, Wergin MC, Verrax J, Rabbani ZN, et al. Targeting lactate-fueled respiration selectively kills hypoxic tumor cells in mice. J Clin Invest (2008) 118(12):3930–42. doi: 10.1172/JCI36843
103. Fiaschi T, Marini A, Giannoni E, Taddei ML, Gandellini P, De Donatis A, et al. Reciprocal metabolic reprogramming through lactate shuttle coordinately influences tumor-stroma interplay. Cancer Res (2012) 72(19):5130–40. doi: 10.1158/0008-5472.CAN-12-1949
104. Guillaumond F, Leca J, Olivares O, Lavaut MN, Vidal N, Berthezene P, et al. Strengthened glycolysis under hypoxia supports tumor symbiosis and hexosamine biosynthesis in pancreatic adenocarcinoma. Proc Natl Acad Sci United States America (2013) 110(10):3919–24. doi: 10.1073/pnas.1219555110
105. Pisarsky L, Bill R, Fagiani E, Dimeloe S, Goosen RW, Hagmann J, et al. Targeting metabolic symbiosis to overcome resistance to anti-angiogenic therapy. Cell Rep (2016) 15(6):1161–74. doi: 10.1016/j.celrep.2016.04.028
106. Jimenez-Valerio G, Martinez-Lozano M, Bassani N, Vidal A, Ochoa-de-Olza M, Suarez C, et al. Resistance to antiangiogenic therapies by metabolic symbiosis in renal cell carcinoma PDX models and patients. Cell Rep (2016) 15(6):1134–43. doi: 10.1016/j.celrep.2016.04.015
107. Allen E, Mieville P, Warren CM, Saghafinia S, Li L, Peng MW, et al. Metabolic symbiosis enables adaptive resistance to anti-angiogenic therapy that is dependent on mTOR signaling. Cell Rep (2016) 15(6):1144–60. doi: 10.1016/j.celrep.2016.04.029
108. Zhao Y, Wang D, Xu T, Liu P, Cao Y, Wang Y, et al. Bladder cancer cells re-educate TAMs through lactate shuttling in the microfluidic cancer microenvironment. Oncotarget (2015) 6(36):39196–210. doi: 10.18632/oncotarget.5538
109. Mu X, Shi W, Xu Y, Xu C, Zhao T, Geng B, et al. Tumor-derived lactate induces M2 macrophage polarization via the activation of the ERK/STAT3 signaling pathway in breast cancer. Cell Cycle (2018) 17(4):428–38. doi: 10.1080/15384101.2018.1444305
110. Angelin A, Gil-de-Gomez L, Dahiya S, Jiao J, Guo L, Levine MH, et al. Foxp3 reprograms T cell metabolism to function in low-glucose, high-lactate environments. Cell Metab (2017) 25(6):1282–93 e7. doi: 10.1016/j.cmet.2016.12.018
111. Wu J, Hong Y, Wu T, Wang J, Chen X, Wang Z, et al. Stromal-epithelial lactate shuttle induced by tumorderived interleukin1beta promotes cell proliferation in oral squamous cell carcinoma. Int J Mol Med (2018) 41(2):687–96. doi: 10.3892/ijmm.2017.3267
112. Whitaker-Menezes D, Martinez-Outschoorn UE, Lin Z, Ertel A, Flomenberg N, Witkiewicz AK, et al. Evidence for a stromal-epithelial “lactate shuttle” in human tumors: MCT4 is a marker of oxidative stress in cancer-associated fibroblasts. Cell Cycle (2011) 10(11):1772–83. doi: 10.4161/cc.10.11.15659
113. Sanita P, Capulli M, Teti A, Galatioto GP, Vicentini C, Chiarugi P, et al. Tumor-stroma metabolic relationship based on lactate shuttle can sustain prostate cancer progression. BMC cancer (2014) 14:154. doi: 10.1186/1471-2407-14-154
114. Becker LM, O’Connell JT, Vo AP, Cain MP, Tampe D, Bizarro L, et al. Epigenetic reprogramming of cancer-associated fibroblasts deregulates glucose metabolism and facilitates progression of breast cancer. Cell Rep (2020) 31(9):107701. doi: 10.1016/j.celrep.2020.107701
115. Ippolito L, Morandi A, Taddei ML, Parri M, Comito G, Iscaro A, et al. Cancer-associated fibroblasts promote prostate cancer Malignancy via metabolic rewiring and mitochondrial transfer. Oncogene (2019) 38(27):5339–55. doi: 10.1038/s41388-019-0805-7
116. Jiang E, Xu Z, Wang M, Yan T, Huang C, Zhou X, et al. Tumoral microvesicle-activated glycometabolic reprogramming in fibroblasts promotes the progression of oral squamous cell carcinoma. FASEB J (2019) 33(4):5690–703. doi: 10.1096/fj.201802226R
117. Pertega-Gomes N, Vizcaino JR, Attig J, Jurmeister S, Lopes C, Baltazar F. A lactate shuttle system between tumour and stromal cells is associated with poor prognosis in prostate cancer. BMC cancer (2014) 14:352. doi: 10.1186/1471-2407-14-352
118. Andersen S, Solstad O, Moi L, Donnem T, Eilertsen M, Nordby Y, et al. Organized metabolic crime in prostate cancer: The coexpression of MCT1 in tumor and MCT4 in stroma is an independent prognosticator for biochemical failure. Urologic Oncol (2015) 33(8):338 e9–17. doi: 10.1016/j.urolonc.2015.05.013
119. Sousa CM, Biancur DE, Wang X, Halbrook CJ, Sherman MH, Zhang L, et al. Pancreatic stellate cells support tumour metabolism through autophagic alanine secretion. Nature (2016) 536(7617):479–83. doi: 10.1038/nature19084
120. Parker SJ, Amendola CR, Hollinshead KER, Yu Q, Yamamoto K, Encarnacion-Rosado J, et al. Selective alanine transporter utilization creates a targetable metabolic niche in pancreatic cancer. Cancer discovery (2020) 10(7):1018–37. doi: 10.1158/2159-8290.CD-19-0959
121. Keren L, Bosse M, Marquez D, Angoshtari R, Jain S, Varma S, et al. A structured tumor-immune microenvironment in triple negative breast cancer revealed by multiplexed ion beam imaging. Cell (2018) 174(6):1373–87 e19. doi: 10.1016/j.cell.2018.08.039
122. Vayrynen JP, Haruki K, Lau MC, Vayrynen SA, Ugai T, Akimoto N, et al. Spatial organization and prognostic significance of NK and NKT-like cells via multimarker analysis of the colorectal cancer microenvironment. Cancer Immunol Res (2022) 10(2):215–27. doi: 10.1158/2326-6066.CIR-21-0772
123. Goltsev Y, Samusik N, Kennedy-Darling J, Bhate S, Hale M, Vazquez G, et al. Deep profiling of mouse splenic architecture with CODEX multiplexed imaging. Cell (2018) 174(4):968–81 e15. doi: 10.1016/j.cell.2018.07.010
124. Rueden CT, Schindelin J, Hiner MC, DeZonia BE, Walter AE, Arena ET, et al. ImageJ2: ImageJ for the next generation of scientific image data. BMC Bioinf (2017) 18(1):529. doi: 10.1186/s12859-017-1934-z
125. Bankhead P, Loughrey MB, Fernandez JA, Dombrowski Y, McArt DG, Dunne PD, et al. QuPath: Open source software for digital pathology image analysis. Sci Rep (2017) 7(1):16878. doi: 10.1038/s41598-017-17204-5
126. McQuin C, Goodman A, Chernyshev V, Kamentsky L, Cimini BA, Karhohs KW, et al. CellProfiler 3.0: Next-generation image processing for biology. PloS Biol (2018) 16(7):e2005970. doi: 10.1371/journal.pbio.2005970
127. de Chaumont F, Dallongeville S, Chenouard N, Herve N, Pop S, Provoost T, et al. Icy: an open bioimage informatics platform for extended reproducible research. Nat Methods (2012) 9(7):690–6. doi: 10.1038/nmeth.2075
128. Schapiro D, Sokolov A, Yapp C, Chen YA, Muhlich JL, Hess J, et al. MCMICRO: a scalable, modular image-processing pipeline for multiplexed tissue imaging. Nat Methods (2022) 19(3):311–5. doi: 10.1038/s41592-021-01308-y
129. Eling N, Damond N, Hoch T, Bodenmiller B. cytomapper: an R/Bioconductor package for visualization of highly multiplexed imaging data. Bioinformatics (2020) 36(24):5706–8. doi: 10.1101/2020.09.08.287516
130. Eliceiri KW, Berthold MR, Goldberg IG, Ibanez L, Manjunath BS, Martone ME, et al. Biological imaging software tools. Nat Methods (2012) 9(7):697–710. doi: 10.1038/nmeth.2084
131. Feng Y, Yang T, Zhu J, Li M, Doyle M, Ozcoban V, et al. Spatial analysis with SPIAT and spaSim to characterize and simulate tissue microenvironments. Nat Commun (2023) 14(1):2697. doi: 10.1038/s41467-023-37822-0
132. Dries R, Zhu Q, Dong R, Eng CL, Li H, Liu K, et al. Giotto: a toolbox for integrative analysis and visualization of spatial expression data. Genome Biol (2021) 22(1):78. doi: 10.1186/s13059-021-02286-2
133. Schapiro D, Jackson HW, Raghuraman S, Fischer JR, Zanotelli VRT, Schulz D, et al. histoCAT: analysis of cell phenotypes and interactions in multiplex image cytometry data. Nat Methods (2017) 14(9):873–6. doi: 10.1038/nmeth.4391
134. BodenmillerGroup. neighbouRhood. (2019). Available at: https://github.com/BodenmillerGroup/neighbouRhood [Accessed July 4, 2023].
135. Gorris MAJ, Halilovic A, Rabold K, van Duffelen A, Wickramasinghe IN, Verweij D, et al. Eight-color multiplex immunohistochemistry for simultaneous detection of multiple immune checkpoint molecules within the tumor microenvironment. J Immunol (Baltimore Md: 1950) (2018) 200(1):347–54. doi: 10.4049/jimmunol.1701262
136. McNamara G, Gupta A, Reynaert J, Coates TD, Boswell C. Spectral imaging microscopy web sites and data. Cytometry Part A: J Int Soc Analytical Cytol (2006) 69(8):863–71. doi: 10.1002/cyto.a.20304
137. Lin JR, Izar B, Wang S, Yapp C, Mei S, Shah PM, et al. Highly multiplexed immunofluorescence imaging of human tissues and tumors using t-CyCIF and conventional optical microscopes. eLife (2018) 7:e31657. doi: 10.7554/eLife.31657
138. Wang Y, Woehrstein JB, Donoghue N, Dai M, Avendano MS, Schackmann RCJ, et al. Rapid sequential in situ multiplexing with DNA exchange imaging in neuronal cells and tissues. Nano letters (2017) 17(10):6131–9. doi: 10.1021/acs.nanolett.7b02716
139. Saka SK, Wang Y, Kishi JY, Zhu A, Zeng Y, Xie W, et al. Immuno-SABER enables highly multiplexed and amplified protein imaging in tissues. Nat Biotechnol (2019) 37(9):1080–90. doi: 10.1038/s41587-019-0207-y
140. Keren L, Bosse M, Thompson S, Risom T, Vijayaragavan K, McCaffrey E, et al. MIBI-TOF: A multiplexed imaging platform relates cellular phenotypes and tissue structure. Sci Adv (2019) 5(10):eaax5851. doi: 10.1126/sciadv.aax5851
141. Merritt CR, Ong GT, Church SE, Barker K, Danaher P, Geiss G, et al. Multiplex digital spatial profiling of proteins and RNA in fixed tissue. Nat Biotechnol (2020) 38(5):586–99. doi: 10.1038/s41587-020-0472-9
142. Eng CL, Lawson M, Zhu Q, Dries R, Koulena N, Takei Y, et al. Transcriptome-scale super-resolved imaging in tissues by RNA seqFISH. Nature (2019) 568(7751):235–9. doi: 10.1038/s41586-019-1049-y
143. Wang X, Allen WE, Wright MA, Sylwestrak EL, Samusik N, Vesuna S, et al. Three-dimensional intact-tissue sequencing of single-cell transcriptional states. Science (2018) 361(6400):eaat5691. doi: 10.1126/science.aat5691
144. Shi H, Shi Q, Grodner B, Lenz JS, Zipfel WR, Brito IL, et al. Highly multiplexed spatial mapping of microbial communities. Nature (2020) 588(7839):676–81. doi: 10.1038/s41586-020-2983-4
145. Lee JH, Daugharthy ER, Scheiman J, Kalhor R, Ferrante TC, Terry R, et al. Fluorescent in situ sequencing (FISSEQ) of RNA for gene expression profiling in intact cells and tissues. Nat Protoc (2015) 10(3):442–58. doi: 10.1038/nprot.2014.191
146. Stahl PL, Salmen F, Vickovic S, Lundmark A, Navarro JF, Magnusson J, et al. Visualization and analysis of gene expression in tissue sections by spatial transcriptomics. Science (2016) 353(6294):78–82. doi: 10.1126/science.aaf2403
147. Vickovic S, Eraslan G, Salmen F, Klughammer J, Stenbeck L, Schapiro D, et al. High-definition spatial transcriptomics for in situ tissue profiling. Nat Methods (2019) 16(10):987–90. doi: 10.1038/s41592-019-0548-y
148. Stickels RR, Murray E, Kumar P, Li J, Marshall JL, Di Bella DJ, et al. Highly sensitive spatial transcriptomics at near-cellular resolution with Slide-seqV2. Nat Biotechnol (2021) 39(3):313–9. doi: 10.1038/s41587-020-0739-1
149. Migliozzi D, Pelz B, Dupouy DG, Leblond AL, Soltermann A, Gijs MAM. Microfluidics-assisted multiplexed biomarker detection for in situ mapping of immune cells in tumor sections. Microsystems nanoengineering (2019) 5:59. doi: 10.1038/s41378-019-0104-z
150. Schurch CM, Bhate SS, Barlow GL, Phillips DJ, Noti L, Zlobec I, et al. Coordinated cellular neighborhoods orchestrate antitumoral immunity at the colorectal cancer invasive front. Cell (2020) 182(5):1341–59 e19. doi: 10.1016/j.cell.2020.07.005
151. Angelo M, Bendall SC, Finck R, Hale MB, Hitzman C, Borowsky AD, et al. Multiplexed ion beam imaging of human breast tumors. Nat Med (2014) 20(4):436–42. doi: 10.1038/nm.3488
152. Zugazagoitia J, Gupta S, Liu Y, Fuhrman K, Gettinger S, Herbst RS, et al. Biomarkers associated with beneficial PD-1 checkpoint blockade in non-small cell lung cancer (NSCLC) identified using high-plex digital spatial profiling. Clin Cancer Res (2020) 26(16):4360–8. doi: 10.1158/1078-0432.CCR-20-0175
153. Milosevic V. Different approaches to Imaging Mass Cytometry data analysis. Bioinf Adv (2023) 3(1):vbad046. doi: 10.1093/bioadv/vbad046
154. Sorin M, Rezanejad M, Karimi E, Fiset B, Desharnais L, Perus LJM, et al. Single-cell spatial landscapes of the lung tumour immune microenvironment. Nature (2023) 614(7948):548–54. doi: 10.1038/s41586-022-05672-3
155. Karimi E, Yu MW, Maritan SM, Perus LJM, Rezanejad M, Sorin M, et al. Single-cell spatial immune landscapes of primary and metastatic brain tumours. Nature (2023) 614(7948):555–63. doi: 10.1038/s41586-022-05680-3
156. Brady L, Kriner M, Coleman I, Morrissey C, Roudier M, True LD, et al. Inter- and intra-tumor heterogeneity of metastatic prostate cancer determined by digital spatial gene expression profiling. Nat Commun (2021) 12(1):1426. doi: 10.1038/s41467-021-21615-4
157. Rodriques SG, Stickels RR, Goeva A, Martin CA, Murray E, Vanderburg CR, et al. Slide-seq: A scalable technology for measuring genome-wide expression at high spatial resolution. Science (2019) 363(6434):1463–7. doi: 10.1126/science.aaw1219
158. Ji AL, Rubin AJ, Thrane K, Jiang S, Reynolds DL, Meyers RM, et al. Multimodal analysis of composition and spatial architecture in human squamous cell carcinoma. Cell (2020) 182(2):497–514 e22. doi: 10.1016/j.cell.2020.05.039
159. Hwang WL, Jagadeesh KA, Guo JA, Hoffman HI, Yadollahpour P, Reeves JW, et al. Single-nucleus and spatial transcriptome profiling of pancreatic cancer identifies multicellular dynamics associated with neoadjuvant treatment. Nat Genet (2022) 54(8):1178–91. doi: 10.1038/s41588-022-01134-8
Keywords: tumor immune microenvironment (TIME), tumor microenvironment (TME), tumor-infiltrating lymphocytes, immune checkpoint blockade, non-cell autonomous interactions, multiplexed imaging
Citation: Cohn DE, Forder A, Marshall EA, Vucic EA, Stewart GL, Noureddine K, Lockwood WW, MacAulay CE, Guillaud M and Lam WL (2023) Delineating spatial cell-cell interactions in the solid tumour microenvironment through the lens of highly multiplexed imaging. Front. Immunol. 14:1275890. doi: 10.3389/fimmu.2023.1275890
Received: 10 August 2023; Accepted: 11 October 2023;
Published: 23 October 2023.
Edited by:
Felix Marsh-Wakefield, Royal Prince Alfred Hospital, AustraliaReviewed by:
Aurobind Vidyarthi, Yale University, United StatesQingfei Pan, St. Jude Children’s Research Hospital, United States
Copyright © 2023 Cohn, Forder, Marshall, Vucic, Stewart, Noureddine, Lockwood, MacAulay, Guillaud and Lam. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Aisling Forder, aforder@bccrc.ca; Wan L. Lam, wanlam@bccrc.ca