Topography of cancer-associated immune cells in human solid tumors

  1. Jakob Nikolas Kather  Is a corresponding author
  2. Meggy Suarez-Carmona
  3. Pornpimol Charoentong
  4. Cleo-Aron Weis
  5. Daniela Hirsch
  6. Peter Bankhead
  7. Marcel Horning
  8. Dyke Ferber
  9. Ivan Kel
  10. Esther Herpel
  11. Sarah Schott
  12. Inka Zörnig
  13. Jochen Utikal
  14. Alexander Marx
  15. Timo Gaiser
  16. Herrmann Brenner
  17. Jenny Chang-Claude
  18. Michael Hoffmeister
  19. Dirk Jäger
  20. Niels Halama  Is a corresponding author
  1. National Center for Tumor Diseases, University Hospital Heidelberg, Germany
  2. German Cancer Consortium, Germany
  3. German Cancer Research Center, Germany
  4. University Hospital RWTH Aachen, Germany
  5. University Medical Center Mannheim, Heidelberg University, Germany
  6. Centre for Cancer Research and Cell Biology, Queen’s University Belfast, United Kingdom
  7. University Hospital Heidelberg, Germany
  8. Tissue Bank of the National Center for Tumor Diseases, Germany
  9. German Cancer Research Center and National Center for Tumor Diseases, Germany
  10. German Cancer Research Centre, Germany
  11. University Cancer Center Hamburg, University Medical Center Hamburg-Eppendorf, Germany
7 figures, 1 table and 7 additional files

Figures

Figure 1 with 1 supplement
Semiautomatic image analysis defines immune cell topography.

(A) Manual delineation of three compartments: outer 500 µm invasive margin, inner 500 µm invasive margin, tumor core. (B) Example of automatic cell detection in a CD3-stained gastric carcinoma …

https://doi.org/10.7554/eLife.36967.002
Figure 1—figure supplement 1
Example images for cell count in all six immunostains.

For all slides with sufficient sample availability, we performed immunostaining for Foxp3, CD68, PD1, CD8, CD163 and CD3. QuPath was used for automatic cell detection and classification. Here, …

https://doi.org/10.7554/eLife.36967.003
Figure 2 with 5 supplements
Cell densities in the tumor core and in the outer invasive margin in the pan-cancer cohort.

Raw cell densities are plotted for each cell type and both major compartments. Gray lines indicate the median density for this cell type. Split at the median, tumors can be classified as cold, hot …

https://doi.org/10.7554/eLife.36967.004
Figure 2—figure supplement 1
Replication experiment for the full tissue analysis pipeline.

60 tissue specimens were randomly selected and stained for CD3 or CD163, scanned and analyzed independently and by blinded observers. The resulting cell counts were similar to the original analysis. …

https://doi.org/10.7554/eLife.36967.005
Figure 2—figure supplement 2
Normalized immune cell counts for cytotoxic T lymphocytes (CD8) and pro-tumor macrophages (CD163).

For all three compartments outer invasive margin [out], inner invasive margin [in] and tumor core [core], normalized cell density is shown. In accordance with previous reports, lymphocytes and …

https://doi.org/10.7554/eLife.36967.006
Figure 2—figure supplement 3
Average cell density percentile score for all compartments in ten tumor entities and six immunostains.

For all tumor types and all immunostains, the average target plot is shown based on the mean value of all tumor samples in each group.

https://doi.org/10.7554/eLife.36967.007
Figure 2—figure supplement 4
Optimal number of lymphocyte topography clusters arising in repeated optimization runs with different methods.

Using lymphocyte densities in three regions in tumors from 168 cancer patients, we looked for the optimal number of clusters (1 to 12) to group the topography phenotypes. We used CD3+ and CD8+ cell …

https://doi.org/10.7554/eLife.36967.008
Figure 2—figure supplement 5
Correlations between cell densities in different spatial compartments.

Pearson’s correlation coefficient is shown for all pairwise comparisons between cell densities in three spatial compartments in tumors from N = 144 cancer patients. Cell densities in the ‘inner …

https://doi.org/10.7554/eLife.36967.009
Distribution of immune topography phenotypes among different tumor types in the pan-cancer cohort.

Analysis for all six immune cell types (A–F) and for all analyzed tumor types (MEL = melanoma, LUAD = lung adeno, LUSC = lung squamous, BLCA = bladder, HNSC = head and neck squamous, STAD = stomach …

https://doi.org/10.7554/eLife.36967.010
Pairwise analysis of immune phenotypes for all immune cell types in the pan-cancer cohort.

For all tissue samples in all tumor types, a pairwise classification into cold-excluded-hot was done for all immune cell types. This analysis was based on the median cutoff for high and low cell …

https://doi.org/10.7554/eLife.36967.011
Bivariate immune phenotypes for each tumor type in the pan-cancer cohort.

We analyzed the concordance between hot-cold-excluded topographies for CD8+ lymphocytes and all other cell types for each tumor type separately. Absolute numbers of patients assigned to each of nine …

https://doi.org/10.7554/eLife.36967.012
Overall similarity between tumor entities based on full immune topography.

Hierarchical clustering based on all normalized cell densities of (A) PD1+exhausted lymphocytes; (B) Foxp3+regulatory T cells; (C) CD8+cytotoxic T lymphocytes; (D) CD68+monocytes/macrophages; (E) …

https://doi.org/10.7554/eLife.36967.013
Prognostic value of the myeloid-lymphoid topography in primary colorectal cancer (CRC) in the DACHS cohort.

In a validation cohort of N = 287 colorectal cancer patients (N = 286 with follow-up data) from the DACHS study, CD8 and CD163 staining of the primary surgical sample was correlated to clinical …

https://doi.org/10.7554/eLife.36967.014

Tables

Key resources table
Reagent type (species) or resourceDesignationSource/ReferenceIdentifierAdditional information
Software,
algorithm
QuPath v0.1.2Bankheadet al.DOI: 10.1038/s41598-017-17204-5-
AntibodyAnti-human
CD3
Leica
Novocastra
RRID:AB_563544Dilution 1:100
AntibodyAnti-human
CD8
Leica
Novocastra
RRID:AB_442068Dilution 1:50
AntibodyAnti-human
Foxp3
eBioscienceRRID:AB_467555Dilution 1:100
AntibodyAnti-human
CD163
BioRadRRID:AB_2074540Dilution 1:500
AntibodyAnti-human
CD68
Thermo Fisher
Scientific
RRID:AB_720547Dilution 1:2000
AntibodyAnti-human
PD1
AbcamRRID:AB_881954Dilution 1:50

Additional files

Supplementary file 1

Clinical characterization of the DACHS cohort.

This table lists summary statistics of all relevant clinico-pathological features of the DACHS cohort.

https://doi.org/10.7554/eLife.36967.015
Supplementary file 2

List of all image analysis parameters.

In this table, all parameters for cell detection and classification using the open source software QuPath are listed. Two sets of parameters are distinguished: ‘DAB’ (diaminobenzidine), used for blue-brown staining, and ‘Red’, used for blue-red staining in melanoma. OD = optical density. All parameters were used in the pan-cancer cohort unless labeled as ‘DACHS’, in which case they were used in the DACHS cohort.

https://doi.org/10.7554/eLife.36967.016
Supplementary file 3

List of all samples and all measurement values of the pan-cancer cohort.

In this table, we report all raw measurements for all samples that were used in this study. Column names are: ‘class’ (tumor type as listed above), ‘patient’ (patient pseudonym), ‘antigen’ (antigen for immunostain), ‘TU_CORE_cells_mm2’ (number of positively stained cells per square millimeter in the tumor core), ‘MARG_500_IN_cells_mm2’ (number of positively stained cells per square millimeter in the inner invasive margin, defined as ranging 0–500 µm to the inside from the tumor edge), ‘MARG_500_OUT_cells_mm2’ (number of positively stained cells per square millimeter in the inner invasive margin, defined as ranging 0–500 µm to the outside from the tumor edge).

https://doi.org/10.7554/eLife.36967.017
Supplementary file 4

List of all cutoff values for all cell types.

On the full data set of N = 965 tissue slides from N = 177 patients in 10 tumor types, we calculated the median cell density for each antigen, taking the compartments ‘outer invasive margin’ and ‘tumor core’ into account. These median values were subsequently used as cutoff values for low and high cell densities which were then used to define hot, cold and excluded phenotypes.

https://doi.org/10.7554/eLife.36967.018
Supplementary file 5

Continuous cell densities of CD8+ and CD163+ cells are not significantly associated with overall survival in colorectal cancer.

A multivariable Cox proportional hazard model was fitted to all variables listed in this table. N = 286 CRC patients in the DACHS cohort, number of events = 108, significance codes (sig): *<0.05, **<0.01, ***<0.001. HR = hazard ratio, UICC = Union internationale contre le cancer.

https://doi.org/10.7554/eLife.36967.019
Supplementary file 6

Bivariate immune phenotype predicts risk of death of any cause.

A multivariable Cox proportional hazard model was fitted to all variables listed in this table. N = 286 CRC patients in the DACHS cohort, number of events = 108, significance codes (sig): *<0.05, **<0.01, ***<0.001. HR = hazard ratio, UICC = Union internationale contre le cancer.

https://doi.org/10.7554/eLife.36967.020
Transparent reporting form
https://doi.org/10.7554/eLife.36967.021

Download links