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

Digital pathology enables computational analysis algorithms to be applied at scale to histological images. An example is the identification of immune cells within solid tumours. Image analysis algorithms can extract precise cell locations from immunohistochemistry slides, but the resulting spatial coordinates, or point patterns, can be difficult to interpret. Since localisation of immune cells within tumours may reflect their functional status and correlates with patient prognosis, novel descriptors of their spatial distributions are of biological and clinical interest. A range of spatial statistics have been used to analyse such point patterns but, individually, these approaches only partially describe complex immune cell distributions. In this study, we apply three spatial statistics to locations of CD68+ macrophages within human head and neck tumours, and show that images grouped semi-quantitatively by a pathologist share similar statistics. We generate a synthetic dataset which emulates human samples and use it to demonstrate that combining multiple spatial statistics with a maximum likelihood approach better predicts human classifications than any single statistic. We can also estimate the error associated with our classifications. Importantly, this methodology is adaptable and can be extended to other histological investigations or applied to point patterns outside of histology.

[1]  A. Madabhushi,et al.  Spatial Architecture and Arrangement of Tumor-Infiltrating Lymphocytes for Predicting Likelihood of Recurrence in Early-Stage Non–Small Cell Lung Cancer , 2018, Clinical Cancer Research.

[2]  Bilwaj Gaonkar,et al.  Cell Counting and Segmentation of Immunohistochemical Images in the Spinal Cord: Comparing Deep Learning and Traditional Approaches , 2018, 2018 40th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC).

[3]  Ursula Klingmüller,et al.  Structural and practical identifiability analysis of partially observed dynamical models by exploiting the profile likelihood , 2009, Bioinform..

[4]  Franz Schweiggert,et al.  On the Classification of Prostate Carcinoma With Methods from Spatial Statistics , 2007, IEEE Transactions on Information Technology in Biomedicine.

[5]  A. Getis The Analysis of Spatial Association by Use of Distance Statistics , 2010 .

[6]  A. J. Baddeley,et al.  A non-parametric measure of spatial interaction in point patterns , 1996, Advances in Applied Probability.

[7]  Souptik Barua,et al.  Spatial interaction of tumor cells and regulatory T cells correlates with survival in non-small cell lung cancer. , 2018, Lung cancer.

[8]  Anant Madabhushi,et al.  A Deep Convolutional Neural Network for segmenting and classifying epithelial and stromal regions in histopathological images , 2016, Neurocomputing.

[9]  Crispin J. Miller,et al.  Prospective technical validation and assessment of intra-tumour heterogeneity of a low density array hypoxia gene profile in head and neck squamous cell carcinoma. , 2013, European journal of cancer.

[10]  Bernd Bodenmiller,et al.  Automatic single cell segmentation on highly multiplexed tissue images , 2015, Cytometry. Part A : the journal of the International Society for Analytical Cytology.

[11]  Sidra Nawaz,et al.  Beyond immune density: critical role of spatial heterogeneity in estrogen receptor-negative breast cancer , 2015, Modern Pathology.

[12]  Helen Pitman,et al.  Artificial intelligence in digital pathology: a roadmap to routine use in clinical practice , 2019, The Journal of pathology.

[13]  D. W. D. Cruz,et al.  Enhanced larval supply and recruitment can replenish reef corals on degraded reefs , 2017, Scientific Reports.

[14]  Anant Madabhushi,et al.  Accurate and reproducible invasive breast cancer detection in whole-slide images: A Deep Learning approach for quantifying tumor extent , 2017, Scientific Reports.

[15]  Marc Lartaud,et al.  Analyzing huge pathology images with open source software , 2013, Diagnostic Pathology.

[16]  Luis A. Escobar,et al.  Teaching about Approximate Confidence Regions Based on Maximum Likelihood Estimation , 1995 .

[17]  Yinyin Yuan,et al.  Modelling the spatial heterogeneity and molecular correlates of lymphocytic infiltration in triple-negative breast cancer , 2015, Journal of The Royal Society Interface.

[18]  Matti Pietikäinen,et al.  Identification of tumor epithelium and stroma in tissue microarrays using texture analysis , 2012, Diagnostic Pathology.

[19]  Andrew H. Beck,et al.  Systematic Analysis of Breast Cancer Morphology Uncovers Stromal Features Associated with Survival , 2011, Science Translational Medicine.

[20]  Jens Rittscher,et al.  The use of digital pathology and image analysis in clinical trials , 2019, The journal of pathology. Clinical research.

[21]  Dietrich Stoyan,et al.  Basic Ideas ofSpatial Statistics , 2000 .

[22]  Lassi Paavolainen,et al.  Systems pathology by multiplexed immunohistochemistry and whole-slide digital image analysis , 2017, Scientific Reports.

[23]  Aki Vehtari,et al.  Understanding predictive information criteria for Bayesian models , 2013, Statistics and Computing.

[24]  Z. Trajanoski,et al.  Type, Density, and Location of Immune Cells Within Human Colorectal Tumors Predict Clinical Outcome , 2006, Science.

[25]  E. Manias,et al.  Improving medication adherence in adult kidney transplantation (IMAKT): A pilot randomised controlled trial , 2019, Scientific reports.

[26]  S. Astley,et al.  Breast screening interval and the characteristics of screen-detected cancers , 2015, Breast Cancer Research.

[27]  Yinyin Yuan Spatial Heterogeneity in the Tumor Microenvironment. , 2016, Cold Spring Harbor perspectives in medicine.

[28]  Iman H Hewedi,et al.  Diagnostic value of progesterone receptor and p53 expression in uterine smooth muscle tumors , 2012, Diagnostic Pathology.

[29]  J. Møller,et al.  Statistical Inference and Simulation for Spatial Point Processes , 2003 .

[30]  Arvydas Laurinavicius,et al.  A methodology to ensure and improve accuracy of Ki67 labelling index estimation by automated digital image analysis in breast cancer tissue , 2014, Breast Cancer Research.

[31]  D. Stoyan,et al.  Fractals, random shapes and point fields : methods of geometrical statistics , 1996 .

[32]  F. Marincola,et al.  Towards the introduction of the ‘Immunoscore’ in the classification of malignant tumours , 2013, The Journal of pathology.

[33]  F. Marincola,et al.  International validation of the consensus Immunoscore for the classification of colon cancer: a prognostic and accuracy study , 2018, The Lancet.

[34]  Jack Cuzick,et al.  Relevance of Spatial Heterogeneity of Immune Infiltration for Predicting Risk of Recurrence After Endocrine Therapy of ER+ Breast Cancer , 2018, Journal of the National Cancer Institute.

[35]  L. Santamaría,et al.  Estimate of the pair correlation functions of nuclei in normal prostate, prostatic intraepithelial neoplasia, and cancer , 2015 .

[36]  Peter Bankhead,et al.  QuPath: Open source software for digital pathology image analysis , 2017, Scientific Reports.

[37]  Carlo C. Maley,et al.  An ecological measure of immune-cancer colocalization as a prognostic factor for breast cancer , 2015, Breast Cancer Research.

[38]  Johannes E. Schindelin,et al.  The ImageJ ecosystem: An open platform for biomedical image analysis , 2015, Molecular reproduction and development.

[39]  Helen M. Byrne,et al.  Techniques for analysing pattern formation in populations of stem cells and their progeny , 2011, BMC Bioinformatics.

[40]  Victoria Sanz-Moreno,et al.  Faculty Opinions recommendation of QuPath: Open source software for digital pathology image analysis. , 2019, Faculty Opinions – Post-Publication Peer Review of the Biomedical Literature.

[41]  A. Baddeley,et al.  A non-parametric measure of spatial interaction in point patterns , 1996, Advances in Applied Probability.

[42]  Torsten Mattfeldt,et al.  Statistical analysis of labelling patterns of mammary carcinoma cell nuclei on histological sections , 2009, Journal of microscopy.

[43]  Zhongyu Yuan,et al.  High infiltration of tumor-associated macrophages in triple-negative breast cancer is associated with a higher risk of distant metastasis , 2014, OncoTargets and therapy.

[44]  R. Rabadán,et al.  Quantitative Analysis of Immune Infiltrates in Primary Melanoma , 2018, Cancer Immunology Research.

[45]  Geert J. S. Litjens,et al.  Computer aided quantification of intratumoral stroma yields an independent prognosticator in rectal cancer , 2019, Cellular Oncology.

[46]  Rajarsi R. Gupta,et al.  Spatial Organization and Molecular Correlation of Tumor-Infiltrating Lymphocytes Using Deep Learning on Pathology Images. , 2018, Cell reports.

[47]  Bram van Ginneken,et al.  Epithelium segmentation using deep learning in H&E-stained prostate specimens with immunohistochemistry as reference standard , 2018, Scientific Reports.

[48]  Gang Wang,et al.  Automatic grading of prostate cancer in digitized histopathology images: Learning from multiple experts , 2018, Medical Image Anal..

[49]  Ching-Wei Wang Robust Automated Tumour Segmentation on Histological and Immunohistochemical Tissue Images , 2011, PloS one.

[50]  Fazly Salleh Abas,et al.  Cell Classification in ER-Stained Whole Slide Breast Cancer Images Using Convolutional Neural Network , 2018, 2018 40th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC).

[51]  Aki Vehtari,et al.  Practical Bayesian model evaluation using leave-one-out cross-validation and WAIC , 2015, Statistics and Computing.

[52]  Laurence Zitvogel,et al.  The immune contexture in cancer prognosis and treatment , 2017, Nature Reviews Clinical Oncology.

[53]  Darren Treanor,et al.  Digital pathology access and usage in the UK: results from a national survey on behalf of the National Cancer Research Institute’s CM-Path initiative , 2018, Journal of Clinical Pathology.

[54]  David B. A. Epstein,et al.  Fast and Accurate Tumor Segmentation of Histology Images using Persistent Homology and Deep Convolutional Features , 2018, Medical Image Anal..