Combining multiple spatial statistics enhances the description of immune cell localisation within tumours

Dissecting the Spatially Restricted Effects of Microenvironment-Mediated Resistance on Targeted Therapy Responses

The response of tumors to anti-cancer therapies is defined not only by cell-intrinsic therapy sensitivities but also by local interactions with the tumor microenvironment. Fibroblasts that make tumor stroma have been shown to produce paracrine factors that can strongly reduce the sensitivity of tumor cells to many types of targeted therapies. Moreover, a high stroma/tumor ratio is generally associated with poor survival and reduced therapy responses. However, in contrast to advanced knowledge of the molecular mechanisms responsible for stroma-mediated resistance, its effect on the ability of tumors to escape therapeutic eradication remains poorly understood. To a large extent, this gap of knowledge reflects the challenge of accounting for the spatial aspects of microenvironmental resistance, especially over longer time frames. To address this problem, we integrated spatial inferences of proliferation-death dynamics from an experimental animal model of targeted therapy responses with spatial mathematical modeling. With this approach, we dissected the impact of tumor/stroma distribution, magnitude and distance of stromal effects. While all of the tested parameters affected the ability of tumor cells to resist elimination, spatial patterns of stroma distribution within tumor tissue had a particularly strong impact.

Using random forests to uncover the predictive power of distance-varying cell interactions in tumor microenvironments

Tumor microenvironments (TMEs) contain vast amounts of information on patient's cancer through their cellular composition and the spatial distribution of tumor cells and immune cell populations. Exploring variations in TMEs between patient groups, as well as determining the extent to which this information can predict outcomes such as patient survival or treatment success with emerging immunotherapies, is of great interest. Moreover, in the face of a large number of cell interactions to consider, we often wish to identify specific interactions that are useful in making such predictions. We present an approach to achieve these goals based on summarizing spatial relationships in the TME using spatial K functions, and then applying functional data analysis and random forest models to both predict outcomes of interest and identify important spatial relationships. This approach is shown to be effective in simulation experiments at both identifying important spatial interactions while also controlling the false discovery rate. We further used the proposed approach to interrogate two real data sets of Multiplexed Ion Beam Images of TMEs in triple negative breast cancer and lung cancer patients. The methods proposed are publicly available in a companion R package funkycells.

Hypoxia coordinates the spatial landscape of myeloid cells within glioblastoma to affect survival

Myeloid cells are highly prevalent in glioblastoma (GBM), existing in a spectrum of phenotypic and activation states. We now have limited knowledge of the tumor microenvironment (TME) determinants that influence the localization and the functions of the diverse myeloid cell populations in GBM. Here, we have utilized orthogonal imaging mass cytometry with single-cell and spatial transcriptomic approaches to identify and map the various myeloid populations in the human GBM tumor microenvironment (TME). Our results show that different myeloid populations have distinct and reproducible compartmentalization patterns in the GBM TME that is driven by tissue hypoxia, regional chemokine signaling, and varied homotypic and heterotypic cellular interactions. We subsequently identified specific tumor subregions in GBM, based on composition of identified myeloid cell populations, that were linked to patient survival. Our results provide insight into the spatial organization of myeloid cell subpopulations in GBM, and how this is predictive of clinical outcome.

Semmaphorin 3 A causes immune suppression by inducing cytoskeletal paralysis in tumour-specific CD8+ T cells

Semaphorin-3A (SEMA3A) functions as a chemorepulsive signal during development and can affect T cells by altering their filamentous actin (F-actin) cytoskeleton. The exact extent of these effects on tumour-specific T cells are not completely understood. Here we demonstrate that Neuropilin-1 (NRP1) and Plexin-A1 and Plexin-A4 are upregulated on stimulated CD8+ T cells, allowing tumour-derived SEMA3A to inhibit T cell migration and assembly of the immunological synapse. Deletion of NRP1 in both CD4+ and CD8+ T cells enhance CD8+ T-cell infiltration into tumours and restricted tumour growth in animal models. Conversely, over-expression of SEMA3A inhibit CD8+ T-cell infiltration. We further show that SEMA3A affects CD8+ T cell F-actin, leading to inhibition of immune synapse formation and motility. Examining a clear cell renal cell carcinoma patient cohort, we find that SEMA3A expression is associated with reduced survival, and that T-cells appear trapped in SEMA3A rich regions. Our study establishes SEMA3A as an inhibitor of effector CD8+ T cell tumour infiltration, suggesting that blocking NRP1 could improve T cell function in tumours.

SpaceANOVA: Spatial Co-occurrence Analysis of Cell Types in Multiplex Imaging Data Using Point Process and Functional ANOVA

Multiplex imaging platforms have enabled the identification of the spatial organization of different types of cells in complex tissue or the tumor microenvironment. Exploring the potential variations in the spatial co-occurrence or colocalization of different cell types across distinct tissue or disease classes can provide significant pathological insights, paving the way for intervention strategies. However, the existing methods in this context either rely on stringent statistical assumptions or suffer from a lack of generalizability. We present a highly powerful method to study differential spatial co-occurrence of cell types across multiple tissue or disease groups, based on the theories of the Poisson point process and functional analysis of variance. Notably, the method accommodates multiple images per subject and addresses the problem of missing tissue regions, commonly encountered due to data-collection complexities. We demonstrate the superior statistical power and robustness of the method in comparison with existing approaches through realistic simulation studies. Furthermore, we apply the method to three real data sets on different diseases collected using different imaging platforms. In particular, one of these data sets reveals novel insights into the spatial characteristics of various types of colorectal adenoma.

Extended correlation functions for spatial analysis of multiplex imaging data

Imaging platforms for generating highly multiplexed histological images are being continually developed and improved. Significant improvements have also been made in the accuracy of methods for automated cell segmentation and classification. However, less attention has focused on the quantification and analysis of the resulting point clouds, which describe the spatial coordinates of individual cells. We focus here on a particular spatial statistical method, the cross-pair correlation function (cross-PCF), which can identify positive and negative spatial correlation between cells across a range of length scales. However, limitations of the cross-PCF hinder its widespread application to multiplexed histology. For example, it can only consider relations between pairs of cells, and cells must be classified using discrete categorical labels (rather than labeling continuous labels such as stain intensity). In this paper, we present three extensions to the cross-PCF which address these limitations and permit more detailed analysis of multiplex images: topographical correlation maps can visualize local clustering and exclusion between cells; neighbourhood correlation functions can identify colocalization of two or more cell types; and weighted-PCFs describe spatial correlation between points with continuous (rather than discrete) labels. We apply the extended PCFs to synthetic and biological datasets in order to demonstrate the insight that they can generate.

Stroma-Mediated Breast Cancer Cell Proliferation Indirectly Drives Chemoresistance by Accelerating Tumor Recovery between Chemotherapy Cycles

The ability of tumors to survive therapy reflects both cell-intrinsic and microenvironmental mechanisms. Across many cancers, including triple-negative breast cancer (TNBC), a high stroma/tumor ratio correlates with poor survival. In many contexts, this correlation can be explained by the direct reduction of therapy sensitivity induced by stroma-produced paracrine factors. We sought to explore whether this direct effect contributes to the link between stroma and poor responses to chemotherapies. In vitro studies with panels of TNBC cell line models and stromal isolates failed to detect a direct modulation of chemoresistance. At the same time, consistent with prior studies, fibroblast-produced secreted factors stimulated treatment-independent enhancement of tumor cell proliferation. Spatial analyses indicated that proximity to stroma is often associated with enhanced tumor cell proliferation in vivo. These observations suggested an indirect link between stroma and chemoresistance, where stroma-augmented proliferation potentiates the recovery of residual tumors between chemotherapy cycles. To evaluate this hypothesis, a spatial agent-based model of stroma impact on proliferation/death dynamics was developed that was quantitatively parameterized using inferences from histologic analyses and experimental studies. The model demonstrated that the observed enhancement of tumor cell proliferation within stroma-proximal niches could enable tumors to avoid elimination over multiple chemotherapy cycles. Therefore, this study supports the existence of an indirect mechanism of environment-mediated chemoresistance that might contribute to the negative correlation between stromal content and poor therapy outcomes.

SIGNIFICANCE

Integration of experimental research with mathematical modeling reveals an indirect microenvironmental chemoresistance mechanism by which stromal cells stimulate breast cancer cell proliferation and highlights the importance of consideration of proliferation/death dynamics. See related commentary by Wall and Echeverria, p. 3667.

SpaceANOVA: Spatial co-occurrence analysis of cell types in multiplex imaging data using point process and functional ANOVA

Motivation

Multiplex imaging platforms have enabled the identification of the spatial organization of different types of cells in complex tissue or tumor microenvironment (TME). Exploring the potential variations in the spatial co-occurrence or co-localization of different cell types across distinct tissue or disease classes can provide significant pathological insights, paving the way for intervention strategies. However, the existing methods in this context either rely on stringent statistical assumptions or suffer from a lack of generalizability.

Results

We present a highly powerful method to study differential spatial co-occurrence of cell types across multiple tissue or disease groups, based on the theories of the Poisson point process (PPP) and functional analysis of variance (FANOVA). Notably, the method accommodates multiple images per subject and addresses the problem of missing tissue regions, commonly encountered in such a context due to the complex nature of the data-collection procedure. We demonstrate the superior statistical power and robustness of the method in comparison to existing approaches through realistic simulation studies. Furthermore, we apply the method to three real datasets on different diseases collected using different imaging platforms. In particular, one of these datasets reveals novel insights into the spatial characteristics of various types of precursor lesions associated with colorectal cancer.

Availability

The associated R package can be found here, https://github.com/sealx017/SpaceANOVA.

The tumour ecology of quiescence: Niches across scales of complexity

Quiescence is a state of cell cycle arrest, allowing cancer cells to evade anti-proliferative cancer therapies. Quiescent cancer stem cells are thought to be responsible for treatment resistance in glioblastoma, an aggressive brain cancer with poor patient outcomes. However, the regulation of quiescence in glioblastoma cells involves a myriad of intrinsic and extrinsic mechanisms that are not fully understood. In this review, we synthesise the literature on quiescence regulatory mechanisms in the context of glioblastoma and propose an ecological perspective to stemness-like phenotypes anchored to the contemporary concepts of niche theory. From this perspective, the cell cycle regulation is multiscale and multidimensional, where the niche dimensions extend to extrinsic variables in the tumour microenvironment that shape cell fate. Within this conceptual framework and powered by ecological niche modelling, the discovery of microenvironmental variables related to hypoxia and mechanosignalling that modulate proliferative plasticity and intratumor immune activity may open new avenues for therapeutic targeting of emerging biological vulnerabilities in glioblastoma.

Immune cellular patterns of distribution affect outcomes of patients with non-small cell lung cancer

Studying the cellular geographic distribution in non-small cell lung cancer is essential to understand the roles of cell populations in this type of tumor. In this study, we characterize the spatial cellular distribution of immune cell populations using 23 makers placed in five multiplex immunofluorescence panels and their associations with clinicopathologic variables and outcomes. Our results demonstrate two cellular distribution patterns-an unmixed pattern mostly related to immunoprotective cells and a mixed pattern mostly related to immunosuppressive cells. Distance analysis shows that T-cells expressing immune checkpoints are closer to malignant cells than other cells. Combining the cellular distribution patterns with cellular distances, we can identify four groups related to inflamed and not-inflamed tumors. Cellular distribution patterns and distance are associated with survival in univariate and multivariable analyses. Spatial distribution is a tool to better understand the tumor microenvironment, predict outcomes, and may can help select therapeutic interventions.

Quantification of spatial and phenotypic heterogeneity in an agent-based model of tumour-macrophage interactions

We introduce a new spatial statistic, the weighted pair correlation function (wPCF). The wPCF extends the existing pair correlation function (PCF) and cross-PCF to describe spatial relationships between points marked with combinations of discrete and continuous labels. We validate its use through application to a new agent-based model (ABM) which simulates interactions between macrophages and tumour cells. These interactions are influenced by the spatial positions of the cells and by macrophage phenotype, a continuous variable that ranges from anti-tumour to pro-tumour. By varying model parameters that regulate macrophage phenotype, we show that the ABM exhibits behaviours which resemble the 'three Es of cancer immunoediting': Equilibrium, Escape, and Elimination. We use the wPCF to analyse synthetic images generated by the ABM. We show that the wPCF generates a 'human readable' statistical summary of where macrophages with different phenotypes are located relative to both blood vessels and tumour cells. We also define a distinct 'PCF signature' that characterises each of the three Es of immunoediting, by combining wPCF measurements with the cross-PCF describing interactions between vessels and tumour cells. By applying dimension reduction techniques to this signature, we identify its key features and train a support vector machine classifier to distinguish between simulation outputs based on their PCF signature. This proof-of-concept study shows how multiple spatial statistics can be combined to analyse the complex spatial features that the ABM generates, and to partition them into interpretable groups. The intricate spatial features produced by the ABM are similar to those generated by state-of-the-art multiplex imaging techniques which distinguish the spatial distribution and intensity of multiple biomarkers in biological tissue regions. Applying methods such as the wPCF to multiplex imaging data would exploit the continuous variation in biomarker intensities and generate more detailed characterisation of the spatial and phenotypic heterogeneity in tissue samples.

Paracrine enhancement of tumor cell proliferation provides indirect stroma-mediated chemoresistance via acceleration of tumor recovery between chemotherapy cycles

The ability of tumors to survive therapy is mediated not only by cell-intrinsic but also cell-extrinsic, microenvironmental mechanisms. Across many cancers, including triple-negative breast cancer (TNBC), a high stroma/tumor ratio correlates with poor survival. In many contexts, this correlation can be explained by the direct reduction of therapy sensitivity by stroma-produced paracrine factors through activating pro-survival signaling and stemness. We sought to explore whether this direct effect contributes to the link between stroma and poor responses to chemotherapies in TNBC. Our in vitro studies with panels of TNBC cell line models and stromal isolates failed to detect a direct modulation of chemoresistance. However, we found that fibroblasts often enhance baseline tumor cell proliferation. Consistent with this in vitro observation, we found evidence of stroma-enhanced TNBC cell proliferation in vivo , in xenograft models and patient samples. Based on these observations, we hypothesized an indirect link between stroma and chemoresistance, where stroma-augmented proliferation potentiates the recovery of residual tumors between chemotherapy cycles. To test this hypothesis, we developed a spatial agent-based model of tumor response to repeated dosing of chemotherapy. The model was quantitatively parameterized from histological analyses and experimental studies. We found that even a slight enhancement of tumor cell proliferation within stroma-proximal niches can strongly enhance the ability of tumors to survive multiple cycles of chemotherapy under biologically and clinically feasible parameters. In summary, our study uncovered a novel, indirect mechanism of chemoresistance. Further, our study highlights the limitations of short-term cytotoxicity assays in understanding chemotherapy responses and supports the integration of experimental and in silico modeling.

Detection and Spatiotemporal Analysis of In-vitro 3D Migratory Triple-Negative Breast Cancer Cells

The invasion of cancer cells into the surrounding tissues is one of the hallmarks of cancer. However, a precise quantitative understanding of the spatiotemporal patterns of cancer cell migration and invasion still remains elusive. A promising approach to investigate these patterns are 3D cell cultures, which provide more realistic models of cancer growth compared to conventional 2D monolayers. Quantifying the spatial distribution of cells in these 3D cultures yields great promise for understanding the spatiotemporal progression of cancer. In the present study, we present an image processing and segmentation pipeline for the detection of 3D GFP-fluorescent triple-negative breast cancer cell nuclei, and we perform quantitative analysis of the formed spatial patterns and their temporal evolution. The performance of the proposed pipeline was evaluated using experimental 3D cell culture data, and was found to be comparable to manual segmentation, outperforming four alternative automated methods. The spatiotemporal statistical analysis of the detected distributions of nuclei revealed transient, non-random spatial distributions that consisted of clustered patterns across a wide range of neighbourhood distances, as well as dispersion for larger distances. Overall, the implementation of the proposed framework revealed the spatial organization of cellular nuclei with improved accuracy, providing insights into the 3 dimensional inter-cellular organization and its progression through time.

Statistical Analysis of Multiplex Immunofluorescence and Immunohistochemistry Imaging Data

Advances in multiplexed single-cell immunofluorescence (mIF) and multiplex immunohistochemistry (mIHC) imaging technologies have enabled the analysis of cell-to-cell spatial relationships that promise to revolutionize our understanding of tissue-based diseases and autoimmune disorders. Multiplex images are collected as multichannel TIFF files; then denoised, segmented to identify cells and nuclei, normalized across slides with protein markers to correct for batch effects, and phenotyped; and then tissue composition and spatial context at the cellular level are analyzed. This chapter discusses methods and software infrastructure for image processing and statistical analysis of mIF/mIHC data.

Advantages of manual and automatic computer-aided compared to traditional histopathological diagnosis of melanoma: A pilot study

BACKGROUND

Cutaneous malignant melanoma (CMM) accounts for the highest mortality rate among all skin cancers. Traditional histopathologic diagnosis may be limited by the pathologists' subjectivity. Second-opinion strategies and multidisciplinary consultations are usually performed to overcome this issue. An available solution in the future could be the use of automated solutions based on a computational algorithm that could help the pathologist in everyday practice. The aim of this pilot study was to investigate the potential diagnostic aid of a machine-based algorithm in the histopathologic diagnosis of CMM.

METHODS

We retrospectively examined excisional biopsies of 50 CMM and 20 benign congenital compound nevi. Hematoxylin and eosin (H&E) stained WSI were reviewed independently by two expert dermatopathologists. A fully automated pipeline for WSI processing to support the estimation and prioritization of the melanoma areas was developed.

RESULTS

The spatial distribution of the nuclei in the sample provided a multi-scale overview of the tumor. A global overview of the lesion's silhouette was achieved and, by increasing the magnification, the topological distribution of the nuclei and the most informative areas of interest for the CMM diagnosis were identified and highlighted. These silhouettes allow the histopathologist to discriminate between nevus and CMM with an accuracy of 96% without any extra information.

CONCLUSION

In this study we proposed an easy-to-use model that produces segmentations of CMM silhouettes at fine detail level.

Different spatial distribution of inflammatory cells in the tumor microenvironment of ABC and GBC subgroups of diffuse large B cell lymphoma

Diffuse Large B-Cell Lymphoma (DLBCL) presents a high clinical and biological heterogeneity, and the tumor microenvironment chracteristics are important in its  progression. The aim of this study was to evaluate tumor T, B cells, macrophages and mast cells distribution in GBC and ABC DLBCL subgroups through a set of morphometric parameters allowing to provide a quantitative evaluation of the morphological features of the spatial patterns generated by these inflammatory cells.   Histological ABC and GCB samples were immunostained for CD4, CD8, CD68, CD 163, and tryptase in order to determine both percentage and position of positive cells in the tissue characterizing their spatial distribution. The results evidenced that cell patterns generated by CD4-, CD8-, CD68-, CD163- and tryptase-positive cell profiles exhibited a significantly higher uniformity index in ABC than in GCB subgroup. The positive-cell distributions appeared clustered in tissues from GCB, while in tissues from ABC such a feature was lower or absent. The combinations of spatial statistics-derived parameters can lead to better predictions of tumor cell infiltration than any classical morphometric method providing a more accurate description of the functional status of the tumor, useful for patient prognosis.

Multiparameter persistent homology landscapes identify immune cell spatial patterns in tumors

Highly resolved spatial data of complex systems encode rich and nonlinear information. Quantification of heterogeneous and noisy data-often with outliers, artifacts, and mislabeled points-such as those from tissues, remains a challenge. The mathematical field that extracts information from the shape of data, topological data analysis (TDA), has expanded its capability for analyzing real-world datasets in recent years by extending theory, statistics, and computation. An extension to the standard theory to handle heterogeneous data is multiparameter persistent homology (MPH). Here we provide an application of MPH landscapes, a statistical tool with theoretical underpinnings. MPH landscapes, computed for (noisy) data from agent-based model simulations of immune cells infiltrating into a spheroid, are shown to surpass existing spatial statistics and one-parameter persistent homology. We then apply MPH landscapes to study immune cell location in digital histology images from head and neck cancer. We quantify intratumoral immune cells and find that infiltrating regulatory T cells have more prominent voids in their spatial patterns than macrophages. Finally, we consider how TDA can integrate and interrogate data of different types and scales, e.g., immune cell locations and regions with differing levels of oxygenation. This work highlights the power of MPH landscapes for quantifying, characterizing, and comparing features within the tumor microenvironment in synthetic and real datasets.

Multi-Class Cell Detection Using Spatial Context Representation

In digital pathology, both detection and classification of cells are important for automatic diagnostic and prognostic tasks. Classifying cells into subtypes, such as tumor cells, lymphocytes or stromal cells is particularly challenging. Existing methods focus on morphological appearance of individual cells, whereas in practice pathologists often infer cell classes through their spatial context. In this paper, we propose a novel method for both detection and classification that explicitly incorporates spatial contextual information. We use the spatial statistical function to describe local density in both a multi-class and a multi-scale manner. Through representation learning and deep clustering techniques, we learn advanced cell representation with both appearance and spatial context. On various benchmarks, our method achieves better performance than state-of-the-arts, especially on the classification task. We also create a new dataset for multi-class cell detection and classification in breast cancer and we make both our code and data publicly available.

Challenges and Opportunities in the Statistical Analysis of Multiplex Immunofluorescence Data

Simple Summary

Immune modulation is considered a hallmark of cancer initiation and progression, and has offered promising opportunities for therapeutic manipulation. Multiplex immunofluorescence (mIF) technology has enabled the tumor immune microenvironment (TIME) to be studied at an increased scale, in terms of both the number of markers and the number of samples. Another benefit of mIF technology is the ability to measure not only the abundance but also the spatial location of multiple cells types within a tissue sample simultaneously, allowing for assessment of the co-localization of different types of immune markers. Thus, the use of mIF technologies have enable researchers to characterize patient, clinical, and tumor characteristics in the hope of identifying patients whom might benefit from immunotherapy treatments. In this review we outline some of the challenges and opportunities in the statistical analyses of mIF data to study the TIME.

Abstract

Immune modulation is considered a hallmark of cancer initiation and progression. The recent development of immunotherapies has ushered in a new era of cancer treatment. These therapeutics have led to revolutionary breakthroughs; however, the efficacy of immunotherapy has been modest and is often restricted to a subset of patients. Hence, identification of which cancer patients will benefit from immunotherapy is essential. Multiplex immunofluorescence (mIF) microscopy allows for the assessment and visualization of the tumor immune microenvironment (TIME). The data output following image and machine learning analyses for cell segmenting and phenotyping consists of the following information for each tumor sample: the number of positive cells for each marker and phenotype(s) of interest, number of total cells, percent of positive cells for each marker, and spatial locations for all measured cells. There are many challenges in the analysis of mIF data, including many tissue samples with zero positive cells or “zero-inflated” data, repeated measurements from multiple TMA cores or tissue slides per subject, and spatial analyses to determine the level of clustering and co-localization between the cell types in the TIME. In this review paper, we will discuss the challenges in the statistical analysis of mIF data and opportunities for further research.