Skip to main content

An ecological measure of immune-cancer colocalization as a prognostic factor for breast cancer

Abstract

Introduction

Abundance of immune cells has been shown to have prognostic and predictive significance in many tumor types. Beyond abundance, the spatial organization of immune cells in relation to cancer cells may also have significant functional and clinical implications. However there is a lack of systematic methods to quantify spatial associations between immune and cancer cells.

Methods

We applied ecological measures of species interactions to digital pathology images for investigating the spatial associations of immune and cancer cells in breast cancer. We used the Morisita-Horn similarity index, an ecological measure of community structure and predator–prey interactions, to quantify the extent to which cancer cells and immune cells colocalize in whole-tumor histology sections. We related this index to disease-specific survival of 486 women with breast cancer and validated our findings in a set of 516 patients from different hospitals.

Results

Colocalization of immune cells with cancer cells was significantly associated with a disease-specific survival benefit for all breast cancers combined. In HER2-positive subtypes, the prognostic value of immune-cancer cell colocalization was highly significant and exceeded those of known clinical variables. Furthermore, colocalization was a significant predictive factor for long-term outcome following chemotherapy and radiotherapy in HER2 and Luminal A subtypes, independent of and stronger than all known clinical variables.

Conclusions

Our study demonstrates how ecological methods applied to the tumor microenvironment using routine histology can provide reproducible, quantitative biomarkers for identifying high-risk breast cancer patients. We found that the clinical value of immune-cancer interaction patterns is highly subtype-specific but substantial and independent to known clinicopathologic variables that mostly focused on cancer itself. Our approach can be developed into computer-assisted prediction based on histology samples that are already routinely collected.

Introduction

In recent years, the recognition that cancer is an evolutionary process has penetrated much of cancer biology and evolutionary biology [1, 2], and a variety of evolutionary approaches have been adapted for use in cancer biology (e.g., diversity measures for predicting progression) [37]. However, cancer is more than just an evolutionary process; it is also an ecological process [8]. Cancer cells utilize resources and construct habitats within the tissues of the body just as organisms do in the natural world. The ecology of cancer is therefore critical for our understanding of the natural forces that shape cancer development, yet has rarely been systematically investigated in tumors. This parallel between organismal ecology and the tumor microenvironment means that there are unrealized opportunities for adopting measures from ecology for understanding the dynamics and selective pressures on a tumor, which may lead to improved cancer prognosis, prediction, risk stratification and therapeutics.

Immune cell infiltration is one of the most important aspects of the ecology of neoplastic cells [911]. An array of studies has shown that the spatial location of immune cells relative to cancer cells is clinically important in many different cancer types [9, 1215]. However, the spatial patterns of associations between cancer cells and immune cells are rarely quantified, and the evolutionary and ecological processes underlying the role of immune infiltration in human tumors are poorly understood. Recently, we discovered that dense concentrations (“hotspots”) formed by both immune and cancer cells, rather than those formed by one cell type alone, are associated with good prognosis in estrogen receptor (ER)-negative breast cancer [16]. This highlights the importance of investigating how cancer and immune cells are spatially related and raises the question: can we characterize the spatial association between cancer and immune cells and thereby elucidate the ecological dynamics that could ultimately influence tumor progression and response to treatment?

Evolution and ecology can provide a framework for understanding these complex dynamics and predicting clinical outcomes [2]. The Morisita-Horn index [17], a measure that has been used to study community structure in ecology [18, 19], can be adapted for quantifying spatial colocalization of immune and cancer cells. The relationship between colocalization (as measured by the Morisita-Horn index) and patient survival can further reveal whether the immune cells are having a pro- or anti-tumor effect and so provide information about the types of ecological interaction occurring in the tumor. Thus, the aims of this study were 1) to use ecological measures for quantifying spatial associations between cancer cells and immune cells, 2) perform a comprehensive evaluation of the spatial parameters involved in that quantification and 3) to investigate the implications of immune cell colocalization with cancer cells on cancer progression across all subtypes of breast cancer in a large patient cohort.

Methods

Clinical samples and histology analysis

Hematoxylin and eosin (H&E)-stained section images representing 1002 consenting patients in the METABRIC study, under ethical approval by relevant review boards as reported in our previous publication [20], were analyzed using our tool CRImage [21]. More details on the sample set can be found in Additional file 1. The principle of our image analysis tool is to identify immune cells on H&E images based on their typical morphology of small, round and homogeneously basophilic nuclei, which differentiates them from other leukocytes such as neutrophils with more pleomorphic nuclei [2, 3, 12]. Cancer cells typically have nuclei of large size and greater variability in texture and shape. They can be differentiated from the generally more elongated nuclei of fibroblasts and endothelial cells. H&E images have previously been used for evaluating immune infiltration in breast cancer [12, 14, 22, 23]. We have performed three experiments to test the accuracy of our image analysis tool: 1) comparison with visual scoring of tumor cellularity and immune infiltrate (Additional file 2: Figure S1, Jonckheere-Terpstra (JT) trend test p <0.0001); 2) 10-fold cross-validation within the training set (90.1 % accuracy) [21]; 3) 10,000 single-cell annotation by an expert pathologist in random samples (correlation R 2 = 0.98; specifically for lymphocyte R 2 = 0.99) [21]. On average three sections from different tumor locations were obtained for each tumor to increase our ability to capture intra-tumor spatial heterogeneity. The patients were split into a discovery cohort (486 samples, 475 with survival data) and a validation cohort (516 samples, 514 with survival data), thus obtaining two cohorts of similar size. Information on chemotherapy (CT), radiotherapy (RT) and hormone therapy (HT) was available for all patients.

Immune infiltrate scores

Visual scores for immune abundance is available for 675 samples provided by three pathologists in the METABRIC consortium in three categories: absent, mild and severe. The automated score for immune abundance is defined as the percentage of immune cells identified by image analysis, and was highly correlated with the visual score (p <0.0001, JT trend test). We used a cutoff of 8 % as reported in our previous study for estrogen receptor-negative (ER−) cancer [21]. An optimal cutoff search did not yield a more significant result for survival analysis in human epidermal growth factor receptor 2-positive (Her2+) tumors. The visual score for immune-cancer cell colocalization was obtained semiquantitatively for 40 randomly selected cases on the same digital images, which were subjected to automated image analysis by a pathologist blinded to the results of the automated scoring, using a four-tiered scale comprising the following intratumoral distribution patterns (modified after the previous publication [24]): focal (disperse circumscribed immune cell aggregates), low (mild immune cell infiltrate irregularly distributed and spatially unassociated with the invasive tumor epithelia), moderate (immune cell infiltrate displaying some spatial association with the invasive tumor epithelia) and marked (dense immune cell infiltrate closely associated with the invasive tumor epithelia).

Other statistical methods

Monotone trend between a continuous variable and a categorical variable was tested using the JT trend test [25]. Association with clinicopathologic variables was tested using the Kruskal-Wallis test [26] or Fisher’s test in the discovery and validation cohorts. Survival analysis was performed with breast cancer-specific 10-year survival data. The Kaplan-Meier estimator was used for patient stratification and the log-rank test was used for testing for differences among groups. The Cox proportional hazards regression model was fitted to the survival data, and hazard ratios (HR) and 95 % confidence intervals were computed to determine the correlation with disease-specific survival, where the log-rank test with p <0.05 was considered significant. Optimal cutoff searching for dichotomizing the proposed indices was carried out by searching stepwise in the discovery cohort from 20 to 80 percentiles at an interval of 1.5. The cutoffs that displayed the highest prognostic significance with the log-rank test were selected and used for the validation cohort. Tests using a decreasing amount of tissue were performed by dividing slides into areas containing 75 %, 50 % and 25 % of square polygons that were used to compute the Morisita index (Additional file 1). R code and data for reproducing all our results are provided (Additional file 3: Sweave file).

Results

Immune-cancer cell colocalization was independent of known parameters of breast cancer

H&E-stained tumor section images representing 1002 primary breast tumors were analyzed using our image analysis tool to identify cancer cells and immune cells based on their morphology (Fig. 1a-b, Table 1, Additional file 1). To study the spatial distribution of cancer cells and immune cells, each H&E image was virtually divided into non-overlapping squares of 250 μm × 250 μm based on the effective cell − cell communication distance [27]. The number of cancer cells and immune cells within each square was counted (Fig. 1c). To calculate colocalization using cell counts, we then applied two statistics to the data: Morisita-Horn’s similarity index and Pearson correlation (Additional file 1). The value of the Morisita-Horn similarity index ranges from 0, indicating no similarity between two community structures, to 1, when the two structures are the same (an equal number of immune cells and cancer cells in each tessellation square), indicating that the two species are highly colocalized (Fig. 1d and e). We observed a good correlation between visual scoring by an expert pathologist and image analysis scores for immune-cancer colocalization in a subset of samples that were randomly selected (p = 0.008, Fig. 1f, “Methods”). Furthermore, the Morisita index was not correlated with clinical and molecular variables including grade, node, size, ER and Her2 status or PAM50 gene set expression subtypes [28] (p >0.05 in at least one cohort).

Fig. 1
figure 1

Measuring spatial colocalization of immune and cancer cells through image analysis and spatial statistics. a Example H&E image of a breast tumor. Three sections obtained from different locations of the tumor were stained with H&E. b Automated image analysis was used to identify cell types (cancer, immune and stromal cells including fibroblasts and endothelial cells) in this image. c Density of cancer cells and immune cells per square after applying a square tessellation to this image; squares with less than a predefined amount of tissue would be excluded from analysis. d and e Schematic over an arbitrary spatial plane demonstrating how colocalization statistics can discriminate a highly segregated cell pattern from a highly colocalized cell pattern. f Significant correlation between the Morisita index and visual scoring for immune-cancer colocalization in 40 randomly selected samples (JT-test p = 0.0084); focal (disperse immune cell aggregates), low (mild infiltrate unassociated with cancer), moderate (some spatial association with cancer) and marked (dense infiltrate closely associated with cancer)

Table 1 Distribution of the clinicopathological characteristics of all samples or Her2-amplified samples split into discovery and validation cohorts according to the Morisita index

Immune-cancer cell colocalization measured by the Morisita-Horn ecological index is associated with good prognosis in unselected breast cancer

We then examined these measures and their association with 10-year breast cancer disease-specific survival. The optimal cutoff for dichotomizing each measure was selected using the discovery cohort (n = 475) and then tested in the independent validation cohort (n = 514, “Methods”). For both measures, a higher score was associated with significantly better disease-specific survival in the discovery cohort (Morisita-Horn index p = 0.0052, HR = 0.6, 95 % CI = 0.42 − 0.86; Pearson correlation p = 0.00096, HR = 0.55, 95 % CI = 0.38 − 0.79; Fig. 2a), suggesting that a high degree of immune-cancer colocalization is associated with good prognosis. While this association for the Morisita-Horn index was confirmed in the validation cohort (p = 0.00067, HR = 0.41, 95 % CI = 0.25 − 0.7), Pearson correlation was not statistically significant in the validation cohort (p = 0.099; Fig. 2b). We further investigated different spatial parameters used for the calculation of the two measures, including a Voronoi vs square tilling/tessellation and different polygon sizes (Additional file 1, Additional file 2: Figure S2-S5). We found that the combination of the Morisita-Horn index with square tessellation is robust to polygon size and the most prognostic, although Voronoi tessellation yielded scores with a higher correlation with visual scoring of colocalization (Pearson correlation p = 5.6 × 10−5), indicating that the visual estimation of colocalization by the pathological criteria used in this study integrates an adjustment for cell density somewhat similar to the Voronoi tessellation. In summary, our data suggest that colocalization of cancer cells and immune cells, as indicated by high scores of the Morisita-Horn index (henceforth referred to as the Morisita index), predicts favorable prognosis in breast cancer.

Fig. 2
figure 2

Association between immune-cancer cell colocalization and disease-specific survival in breast cancer. Kaplan-Meier curves illustrate disease-specific survival probabilities of patient groups in two subsets stratified by the Morisita index and Pearson correlation in the discovery (a) and validation (b) cohorts. The thresholds for dichotomizing two indices were optimized in the discovery cohort and then used without modification in the validation cohort (0.6940734 for the Morisita index and 0.4884236 for Pearson correlation). Numbers in the legend show the number of patients in each group and numbers in brackets show the number of disease-specific deaths

Immune-cancer cell colocalization is associated with a good prognosis in Her2+ and luminal A tumors

We investigated the association between colocalization with the intrinsic molecular subtypes defined by the PAM50 gene set. Using an optimal threshold derived from the same subtype in the discovery cohort (Additional file 2: Figure S6), we found a significant association with favorable prognosis in the Her2 subtype for both cohorts (discovery p = 0.00017, HR = 0.29, 95 % CI = 0.15 − 0.57; validation p = 0.0022, HR = 0.12, 95 % CI = 0.02 − 0.61; Table 2; Fig. 3) and a significant association with favorable prognosis in the discovery and a borderline effect in the validation cohort for the luminal A subtype (discovery p = 0.0014, HR = 0.13, 95 % CI: 0.04–0.45; Validation p = 0.059, HR = 0.29, CI: 0.07–1.14; Table 2). In basal or luminal B subtypes, there was no association with prognosis in the discovery or validation cohort (Additional file 2: Table S2).

Table 2 Prognostic value of immune-cancer cell colocalization measures using univariate (white background) and multivariate (gray background) Cox regression in the discovery and validation cohorts for all breast cancers and Her2+ cancer defined by PAM50 and by SNP6
Fig. 3
figure 3

The association between immune-cancer cell colocalization and breast cancer prognosis is highly subtype-specific. Kaplan-Meier curves for the validation cohort alone are shown for human epidermal growth factor receptor 2-positive (Her2+) (a), basal (b), luminal A (c), and luminal B (d) Pam50 subtypes. The thresholds for dichotomizing two indices were optimized in the discovery cohort of the subtype and then used for both the discovery and validation cohorts (0.639985 for luminal A, 0.7439517 for luminal B, 0.699701 for basal and 0.7106531 for Her2)

Immune-cancer cell colocalization outperforms other scores of immune infiltrate in Her2+ cancer

To further examine the prognostic value of immune-cancer colocalization for Her2+ tumors, we compared the Morisita index with visual and automated scores of immune abundance in the Her2+ subtype defined using PAM50 gene expression-based subtyping or using Her2 amplification status from SNP6 copy number data (both cohorts combined, Fig. 4a-b). The Morisita index appeared to be superior to visual and automated scoring of immune abundance for patient stratification, highlighting the importance of investigating spatial patterns beyond cell abundance. In addition, two automated scores of immune response, namely intra-tumor lymphocyte ratio [29] and hotspots [16] were not found to be prognostic in these Her2+ subtypes (p >0.05 in the validation cohort). For both definitions of Her2+ cancer, the subtype-specific cutoffs for dichotomizing the Morisita index were the same (0.71). High Morisita index was again significantly correlated with good prognosis in the Her2-amplified samples (univariate analysis discovery: p = 0.039, HR = 0.49, 95 % CI = 0.25 − 0.98; validation: p = 9 × 10−5, HR = 0.05, 95 % CI = 0.01 − 0.4; Table 2).

Fig. 4
figure 4

Immune-cancer cell colocalization is a strong prognostic factor in the Her2+ subtype, while immune cell abundance is not. Kaplan-Meier curves to show differences in disease-specific survival stratified by the Morisita index, visual and automated scores of immune abundance in human epidermal growth factor receptor 2 (Her2) subtype defined by PAM50 (a) and Her2-amplified (b) samples. Automated immune abundance was estimated as the percentage of cells that are lymphocytes in H&E images using a cutoff of 8 % (“Methods”)

Immune-cancer colocalization is an independent prognostic factor in Her2+ cancer

Next, we examined the Morisita index with multivariate analyses for each cohort separately in Her2+ tumors. For both definitions of the Her2+ subtype, multivariate analysis showed prognostic value of the Morisita index independent of and stronger than those of node, size and grade in both cohorts (Her2-amplified discovery: p = 0.027, HR = 0.44, 95 % CI: 0.21–0.91; validation: p = 0.0031, HR = 0.03, 95 % CI = 0–0.3; Table 2). Multivariate analysis revealed additional prognostic value of the Morisita index given any of the clinicopathologic variables tested including grade, node size, ER status, immune abundance and TP53 mutation (Additional file 2: Figure S7). Moreover, this association was consistently observed when using bootstrap analysis to test robustness due to small sample sizes in some of the comparisons (Additional file 1, Additional file 2: Figure S8). In another test using decreasing amounts of tissue, the Morisita index remained prognostic in Her2+ cancer with as little as 50 % of tissue (Additional file 1, Additional file 2: Figure S9-10).

Immune-cancer cell colocalization predicts long-term outcome after chemotherapy and radiotherapy in Her2+ cancer

We investigated the influence of treatment options including CT, RT and HT on the Her2-amplified subtype and whether outcomes from these treatments were different according to the Morisita index. Patients did not receive anti-Her2 therapies. We found that effect of patient stratification by Morisita was independent of treatments given to these patients (multivariate analysis with Morisita, CT, RT and HT: Morisita p = 1.3 × 10−4, HR = 0.28, 95 % CI = 0.15 − 0.54; Additional file 2: Table S3). Subsequently, we examined outcome differences separately for each type of treatment given the Morisita index. Patients with late-stage disease were more likely to be given aggressive treatments such as CT and RT, and this was shown in their long-term prognosis (Fig. 5a-c). Nevertheless, CT-treated patients with a high Morisita index had a 65 % survival probability 10 years after diagnosis, significantly better than the 19 % survival probability for patients with a low Morisita index (Fig. 5d, solid red versus dashed red curve, p = 0.00065, HR = 0.26, 95 % CI = 0.11 − 0.59). Furthermore, the Morisita index further stratified patients who had not received CT (chemo-naïve) (solid blue versus dashed blue curve; p = 0.033, HR = 0.34, 95 % CI = 0.12 − 0.96) and had a 10-year survival probability of 80 % (solid black curve, Fig. 5d), an outlook rare for this subtype before the availability of Her2-inhibitors. Similarly, for other treatments, a low Morisita index also further pinpointed aggressive cancers in the RT-treated (90 patients) and HT-naïve groups (60 patients), after controlling for all these clinicopathologic variables (Fig. 5e and f, Additional file 2: Table S3).

Fig. 5
figure 5

Immune-cancer cell colocalization predicts long-term outcome following chemotherapy (CT), radiotherapy (RT) and hormone therapy (HT) in human epidermal growth factor receptor 2-positive (Her2+) breast cancer. Kaplan-Meier curves show differences in disease-specific survival between Her2+ patients with or without CT (a), RT (b) and HT (c). The Morisita index significantly stratifies disease-specific survival for both treated and untreated Her2+ patients, shown by CT treatment and Morisita index (d), RT treatment and Morisita index (e) and HT treatment and Morisita index (f)

Immune-cancer cell colocalization is clinical relevant for long-term treatment outcome in the luminal A subtype

In a similar vein, Morisita was predictive of outcome in patients with luminal A breast cancer who were chemo-naïve (273 patients, p = 0.0006, HR = 0.14, 95 % CI = 0.05 − 0.44), RT-treated (209 patients, p = 0.03, HR = 0.26, 95 % CI = 0.08-0.89), RT- naïve (119 patients, p = 0.005, HR = 0.13, 95 % CI = 0.03 − 0.54), and HT-treated (273 patients, p = 0.0035, HR = 0.25, 95 % CI = 0.1 − 0.64; Fig. 6a-c). In addition, when all treatments and clinical parameters including grade, node, size were considered together, Morisita was the only variable that was significant in a multivariate model (p = 0.0008, HR = 0.22, 95 % CI = 0.19 − 0.53) and was further validated with bootstrap analysis (>95 % of the times Morisita was significant with more than 75 % (n=225) of luminal A patients being randomly taken), (Fig. 6d). We found a clear advantage of the Morisita index compared with image analysis-based (Morisita: p=4.2x10-5, HR=0.19, 95% CI=0.08-0.47; IA: p=0.0057, HR=0.32, 95% CI=0.14-0.75; Fig. 6e and f) or visual scoring for immune abundance (p >0.1).

Fig. 6
figure 6

Immune-cancer cell colocalization predicts long-term outcome following radiotherapy (RT) and hormone therapy (HT) in luminal A breast cancer. Kaplan-Meier curves show differences in disease-specific survival stratified by CT treatment and Morisita index (a), RT treatment and Morisita index (b) and HT treatment and Morisita index (c). (d) Percentage of times where Morisita was statistically significant in univariate and multivariate analysis with different amounts of patient samples that were selected randomly 1000 times. Comparison of Morisita (e) and automated score of immune abundance (f) based on association with disease-specific survival in luminal A cancer

Discussion

There is a wide diversity of outcomes for breast cancer, and differential responses to treatments. Reproducible biomarkers that can identify which patients are likely to follow a benign course and which would benefit from chemotherapy or radiotherapy would dramatically improve clinical care. We combined digital pathology with an ecological measure, the Morisita index, to study the spatial colocalization of cancer cells and immune cells. This is analogous to analyzing predator–prey interactions in an ecosystem. Our study, in 1002 breast cancer patients using routine histology sections stained with H&E, reveals that colocalization of cancer cells and immune cells is a significant indicator of favorable survival, independent of and stronger than standard clinical variables, particularly in Her2+ breast cancer. It is a fully reproducible and highly prognostic measure of immune response. This biomarker could be developed into computer-assisted prediction tools for routine pathology and clinical use.

The association between immune-cancer cell colocalization and favorable prognosis is indicative of continued effectiveness of the immune system to identify and target cancer cells to limit their survival, proliferation or invasion. Just as prey can evolve complex predator avoidance adaptations, cancer cells can likewise evolve complex adaptations to evade immune predation [30]. The Morisita index, which has been used to study predator–prey relationships in ecology, thus provides a quantitative measure of immune-cancer colocalization and potential immune predation in specific breast cancers.

Our ecological measure further revealed clinically relevant information about long-term benefit of specific treatments in Her2+ breast cancer. Patients with late-stage disease are more likely to be treated with chemotherapy, yet these patients with a high Morisita index had significantly better outcomes than patients with a low Morisita index, revealing a strong link between immune predation and long-term benefit from chemotherapy even in aggressive Her2+ cancer. This is consistent with evidence that immune-infiltrated Her2+ cancer is highly sensitive to immune-mediated cytotoxic treatments such as chemotherapy [31, 32]. Although our automated scoring of immune-cancer colocalization is conceptually different from the stromal lymphocytic infiltration reported to correlate with benefit from chemotherapy in the neoadjuvant setting [14], our findings support an important role for immune cells in the success of therapy and provide a tool for identifying high-risk patients following chemotherapy for continuous monitoring, and the two types of scores remain to be compared.

Since these data were collected, the standard of care for Her2+ patients has changed with the introduction of Her2-targeted therapies. However, it has been argued that an active immune response is also critical for Her2-inibition treatment [33, 34]. Certain Her2 inhibitors such as trastuzumab are known to mediate their tumor killing effects through the immune system [35, 36]. Furthermore, promising new treatment paradigms for these cancers may involve the addition of anti-PD1 or anti-CD137 immunotherapies to stimulate IFNγ-producing CD8+ T cells and maximize benefits from Her2 inhibitors such as trastuzumab [34] or other therapies [37, 38]. Therefore, our proposed ecological measure of immune predation may be also effective for identifying patients who could benefit from such therapies to reactivate immune programs and for monitoring the success of those therapies in stimulating immune predation of the cancer cells.

Our data show a significant difference in prognosis or long-term benefit of therapies according to immune-cancer colocalization in Her2+ and luminal A tumors, but not triple-negative or basal tumors. However, abundance of lymphocytes has been found to be a prognostic biomarker in triple-negative breast cancer in multiple studies [39, 40], including our previous report on the triple-negative subtype in our cohort [29]. These suggest that there are different ecological dynamics with immune predation in basal versus Her2+ subtypes, and that our ecological measure reveals a distinction in a specific pattern of immune response that was not previously recognized. Our results also highlight the subtype specificity of immune infiltrate patterns and warrant a comprehensive future study on this, including the spatial clustering pattern we recently reported [16], across all breast cancer subtypes.

A limitation in our study is the lack of immunohistochemistry experiments to identify the types of immune cells that are infiltrating (or being excluded from) the tumor. Studying spatial colocalization in the context of immune subsets including T-regulatory and T-effector cells aided by multi-color staining will be a powerful approach for defining the spatial and molecular heterogeneity of the tumor microenvironment, and we are currently pursuing this. In addition, further validation is needed for stratification based on treatment options. We observed increasing variability in the Morisita scores with decreasing amounts of tissue in our test of intra-slide variability with decreasing amounts of tissue (Additional file 1, Additional file 2: Figure S9-10). This is not unexpected, because the Morisita index is a global statistic of colocalization pattern and so loses statistical power as the amount of tumor material decreases. We recommend using tumor sections and not needle biopsy or tissue microarrays for these measures. We found that Morisita remains prognostic with 50 % of tissue (mean 33,826 cells), suggesting that the Morisita index can be applied as a prognostic marker when tissue area is sufficiently large (>46,230 cells, third quantile of cell numbers in 50 % tissue).

Conclusions

Taken together, we postulate that immune predation exhibits different spatial patterns under specific contexts, which should be accounted for during biomarker development for a specific cancer subtype. The addition of a spatial analysis to the evaluation of the types of immune cells in the microenvironment with immunohistochemistry experiments can reveal further details in the ecological interactions between the neoplastic cells and the immune system. In summary, our study demonstrates the power of applying quantitative image analysis, spatial statistics and ecological theory to the study of tumor microenvironments and the routine pathological assessment of breast cancer. These approaches can be used both for clinical benefit and for revealing the ecological dynamics that are driving any cancer.

Data and software availability

All data have been deposited at the European Genome-Phenome Archive hosted by the European Bioinformatics Institute [EGAS00000000083]. R code and the dataset for performing the proposed methodology and reproducing reported results are provided as Additional file 3: Sweave file for full reproducibility.

Abbreviations

CT:

chemotherapy

ER:

estrogen receptor

H&E:

hematoxylin and eosin

Her2:

human epidermal growth factor receptor 2

HR:

hazard ratio

HT:

hormone therapy

IFN:

interferon

JT:

Jonckheere-Terpstra

RT:

radiotherapy

References

  1. Greaves M, Maley CC. Clonal evolution in cancer. Nature. 2012;481:306–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Aktipis CA, Boddy AM, Gatenby R, Brown JS, Maley CC. Life history trade-offs in cancer evolution. Nat Rev Cancer. 2013;13:883–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Maley CC, Galipeau PC, Finley JC, Wongsurawat VJ, Li X, Sanchez CA, et al. Genetic clonal diversity predicts progression to esophageal adenocarcinoma. Nat Genet. 2006;38:468–73.

    Article  CAS  PubMed  Google Scholar 

  4. Park SY, Gonen M, Kim HJ, Michor F, Polyak K. Cellular and genetic diversity in the progression of in situ human breast carcinomas to an invasive phenotype. J Clin Invest. 2010;120:636–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Merlo LMF, Shah NA, Li X, Blount PL, Vaughan TL, Reid BJ, et al. A comprehensive survey of clonal diversity measures in Barrett’s esophagus as biomarkers of progression to esophageal adenocarcinoma. Cancer Prev Res. 2010;3:1388–97.

    Article  Google Scholar 

  6. Anderson AR, Hassanein M, Branch KM, Lu J, Lobdell NA, Maier J, et al. Microenvironmental independence associated with tumor progression. Cancer Res. 2009;69:8797–806.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Gerlinger M, Horswell S, Larkin J, Rowan AJ, Salm MP, Varela I, et al. Genomic architecture and evolution of clear cell renal cell carcinomas defined by multiregion sequencing. Nat Genet. 2014;46:225–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Merlo LMF, Pepper JW, Reid BJ, Maley CC. Cancer as an evolutionary and ecological process. Nat Rev Cancer. 2006;6:924–35.

    Article  CAS  PubMed  Google Scholar 

  9. Galon J, Costes A, Sanchez-Cabo F, Kirilovsky A, Mlecnik B, Lagorce-Pagès C, et al. Type, density, and location of immune cells within human colorectal tumors predict clinical outcome. Science. 2006;313:1960–4.

    Article  CAS  PubMed  Google Scholar 

  10. Fridman WH, Pages F, Sautes-Fridman C, Galon J. The immune contexture in human tumours: impact on clinical outcome. Nat Rev Cancer. 2012;12:298–306.

    Article  CAS  PubMed  Google Scholar 

  11. Anitei MG, Zeitoun G, Mlecnik B, Marliot F, Haicheur N, Todosi AM, et al. Prognostic and predictive values of the immunoscore in patients with rectal cancer. Clin Cancer Res. 2014;20:1891–9.

    Article  PubMed  Google Scholar 

  12. Issa-Nummer Y, Darb-Esfahani S, Loibl S, Kunz G, Nekljudova V, Schrader I, et al. Prospective validation of immunological infiltrate for prediction of response to neoadjuvant chemotherapy in HER2-negative breast cancer--a substudy of the neoadjuvant GeparQuinto trial. PLoS One. 2013;8, e79775.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Mahmoud SM, Paish EC, Powe DG, Macmillan RD, Grainge MJ, Lee AHS, et al. Tumor-infiltrating CD8+ lymphocytes predict clinical outcome in breast cancer. J Clin Oncol. 2011;29:1949–55.

    Article  PubMed  Google Scholar 

  14. Loi S, Sirtaine N, Piette F, Salgado R, Viale G, Van Eenoo F, et al. Prognostic and predictive value of tumor-infiltrating lymphocytes in a phase III randomized adjuvant breast cancer trial in node-positive breast cancer comparing the addition of docetaxel to doxorubicin with doxorubicin-based chemotherapy: BIG 02–98. J Clin Oncol. 2013;31:860–7.

    Article  CAS  PubMed  Google Scholar 

  15. Denkert C, Loibl S, Noske A, Roller M, Muller BM, Komor M, et al. Tumor-associated lymphocytes as an independent predictor of response to neoadjuvant chemotherapy in breast cancer. J Clin Oncol. 2010;28:105–13.

    Article  CAS  PubMed  Google Scholar 

  16. Nawaz S, Heindl A, Koelble K, Yuan Y. Beyond immune density: critical role of spatial heterogeneity in estrogen receptor-negative breast cancer. Mod Pathol. 2015. doi:10.1038/modpathol.2015.37. United States & Canadian Academy of Pathology.

  17. Horn H. Measurement of “overlap” in comparative ecological studies. Am Nat. 1966;100:419–24.

    Article  Google Scholar 

  18. McIntosh AR, Peckarsky BL, Taylor BW. Predator-induced resource heterogeneity in a stream food web. Ecology. 2004;85:2279–90.

    Article  Google Scholar 

  19. Scalon JD, Avelar MBL, de Alves GF, Zacarias MS. Spatial and temporal dynamics of coffee-leaf-miner and predatory wasps in organic coffee field in formation. Ciência Rural. 2011;41:646–52.

    Article  Google Scholar 

  20. Curtis C, Shah SP, Chin S-F, Turashvili G, Rueda OM, Dunning MJ, et al. The genomic and transcriptomic architecture of 2,000 breast tumours reveals novel subgroups. Nature. 2012;486:346–52.

    CAS  PubMed  PubMed Central  Google Scholar 

  21. Yuan Y, Failmezger H, Rueda OM, Ali HR, Graf S, Chin S-FF, et al. Quantitative image analysis of cellular heterogeneity in breast tumors complements genomic profiling. Sci Transl Med. 2012;4:157ra143.

  22. Denkert C, von Minckwitz G, Brase JC, Sinn BV, Gade S, Kronenwett R, et al. Tumor-infiltrating lymphocytes and response to neoadjuvant chemotherapy with or without carboplatin in human epidermal growth factor receptor 2–positive and triple-negative primary breast cancers. J Clin Oncol. 2015;33:983–91.

    Article  CAS  PubMed  Google Scholar 

  23. Salgado R, Denkert C, Demaria S, Sirtaine N, Klauschen F, Pruneri G, et al. The evaluation of tumor-infiltrating lymphocytes (TILs) in breast cancer: recommendations by an International TILs Working Group 2014. Ann Oncol. 2014. doi:10.1093/annonc/mdu450. Epub 2014 Sep 11.

  24. Mohammed ZMA, Going JJ, Edwards J, Elsberger B, Doughty JC, McMillan DC. The relationship between components of tumour inflammatory cell infiltrate and clinicopathological factors and survival in patients with primary operable invasive ductal breast cancer. Br J Cancer. 2012;107:864–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Jonckheere AR. A test of significance for the relation between m rankings and k ranked categories. Br J Stat Psychol. 1954;7:93–100.

    Article  Google Scholar 

  26. Kruskal WH, Wallis WA. Use of ranks in one-criterion variance analysis. J Am Stat Assoc. 1952;47:583–621.

    Article  Google Scholar 

  27. Francis K, Palsson BO. Effective intercellular communication distances are determined by the relative time constants for cyto/chemokine secretion and diffusion. Proc Natl Acad Sci U S A. 1997;94:12258–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Sorlie T, Tibshirani R, Parker J, Hastie T, Marron JS, Nobel A, et al. Repeated observation of breast tumor subtypes in independent gene expression data sets. Proc Natl Acad Sci U S A. 2003;100:8418–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Yuan Y. Modelling the spatial heterogeneity and molecular correlates of lymphocytic infiltration in triple-negative breast cancer. J R Soc Interface. In press.

  30. Seliger B. Strategies of tumor immune evasion. BioDrugs. 2005;19:347–54.

    Article  CAS  PubMed  Google Scholar 

  31. Perez EA, Thompson EA, Ballman KV, Anderson SK, Asmann YW, Kalari KR, et al. Genomic analysis reveals that immune function genes are strongly linked to clinical outcome in the North Central Cancer Treatment Group N9831 Adjuvant Trastuzumab Trial. J Clin Oncol. 2015;33:701–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Andre F, Dieci MV, Dubsky P, Sotiriou C, Curigliano G, Denkert C, et al. Molecular pathways: involvement of immune pathways in the therapeutic response and outcome in breast cancer. Clin Cancer Res. 2013;19:28–33.

    Article  CAS  PubMed  Google Scholar 

  33. Park S, Jiang Z, Mortenson ED, Deng L, Radkevich-Brown O, Yang X, et al. The therapeutic effect of anti-HER2/neu antibody depends on both innate and adaptive immunity. Cancer Cell. 2010;18:160–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Stagg J, Loi S, Divisekera U, Ngiow SF, Duret H, Yagita H, et al. Anti-ErbB-2 mAb therapy requires type I and II interferons and synergizes with anti-PD-1 or anti-CD137 mAb therapy. Proc Natl Acad Sci U S A. 2011;108:7142–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Gennari R, Menard S, Fagnoni F, Ponchio L, Scelsi M, Tagliabue E, et al. Pilot study of the mechanism of action of preoperative trastuzumab in patients with primary operable breast tumors overexpressing HER2. Clin Cancer Res. 2004;10:5650–5.

    Article  CAS  PubMed  Google Scholar 

  36. Clynes RA, Towers TL, Presta LG, Ravetch JV. Inhibitory Fc receptors modulate in vivo cytotoxicity against tumor targets. Nat Med. 2000;6:443–6.

    Article  CAS  PubMed  Google Scholar 

  37. Kershaw MH, Devaud C, John LB, Westwood JA, Darcy PK. Enhancing immunotherapy using chemotherapy and radiation to modify the tumor microenvironment. Oncoimmunology. 2013;2, e25962.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Formenti SC, Demaria S. Combining radiotherapy and cancer immunotherapy: a paradigm shift. J Natl Cancer Inst. 2013;105:256–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Adams S, Gray RJ, Demaria S, Goldstein L, Perez EA, Shulman LN, et al. Prognostic value of tumor-infiltrating lymphocytes in triple-negative breast cancers from two phase III randomized adjuvant breast cancer trials: ECOG 2197 and ECOG 1199. J Clin Oncol. 2014;32(27):2959-66. JCO.2013.55.0491.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Loi S, Michiels S, Salgado R, Sirtaine N, Jose V, Fumagalli D, et al. Tumor infiltrating lymphocytes are prognostic in triple negative breast cancer and predictive for trastuzumab benefit in early breast cancer: results from the FinHER trial. Ann Oncol. 2014;25:1544–50.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This work was supported in part by Wellcome Trust 105104/Z/14/Z, NIH grants P01 CA91955, R01 CA149566, R01 CA170595, R01 CA185138 and R01 CA140657, F32 CA144331 and CDMRP Breast Cancer Research Program Award BC132057. RN is funded by a career development fellowship from Breast Cancer Now May 2011, SF01. We acknowledge National Health Service (NHS) funding from the National Institute for Health Research (NIHR) Biomedical Research Centre. We thank the Molecular Taxonomy of Breast Cancer International Consortium.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Yinyin Yuan.

Additional information

Competing interests

The authors declared that they have no competing interests.

Authors’ contributions

YY conceived the study; CCM, CAA and YY designed the experiments; KK provided histopathology expertise; CCM, CAA, KK, RN and YY analyzed the data, and CCM, CAA, KK, RN and YY wrote the manuscript. All authors have read and approve of the final version of the manuscript.

Additional files

Additional file 1:

Supplementary methods. Supplementary information on sample materials and methods. (DOCX 116 kb)

Additional file 2: Figure S1.

Supplementary Figures and Tables. This file contains the following figures and tables: Validation of image analysis with pathological scores. Figure S2. Comparison of Voronoi and Square tessellation. Figure S3. Comparing the Morisita-Horn index with Pearson correlation. Figure S4. Heatmap to show correlation among Pearson correlation and the Morisita-Horn index computed over different scales with square (S) or Voronoi (V) tessellation. Figure S5. Association of immune-cancer correlation with survival according to different spatial scales. Figure S6. Optimizing the Morisita index cutoff in the discovery cohort. Figure S7. Immune-cancer colocalization sub-stratifies clinical parameters and immune abundance, estrogen receptor (ER) status and TP53 mutation status in human epidermal growth factor receptor 2 (HER2) + subtype. Figure S8. Robustness of prognostic value of Morisita index in HER2-amplified subtype. Figure S9. Intra-slide score variability. Figure S10. Kaplan-Meier curves of the Morisita index estimated using decreasing amount of tissue. Table S1. Prognostic value of immune-cancer colocalization measures with Voronoi tessellation using univariate and multivariate Cox regression results in two independent cohorts. Table S2. Univariate and multivariate Cox regression analysis result with the Morisita index, node, size and grade in four PAM50 subtypes and HER2-amplified subtype defined by SNP6 microarray with two cohorts. Table S3. Groups of multivariate Cox proportional hazard analysis of the Morisita index and treatment, known clinical or other variables in HER2-amplified, HER2-amplified and chemotherapy (CT)-treated, and HER2-amplified and radiotherapy (RT)-treated patients. (DOCX 6725 kb)

Additional file 3:

Sweave file. R code used for reproducing all results presented in the paper. (PDF 600 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Maley, C.C., Koelble, K., Natrajan, R. et al. An ecological measure of immune-cancer colocalization as a prognostic factor for breast cancer. Breast Cancer Res 17, 131 (2015). https://doi.org/10.1186/s13058-015-0638-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13058-015-0638-4

Keywords