{'fraction_non_alphanumeric': 0.026813511830635118, 'fraction_numerical': 0.023894769613947694, 'mean_word_length': 4.9830034924330615, 'pattern_counts': {'":': 0, '<': 7, '<?xml version=': 0, '>': 2, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Detailed resolution analysis reveals spatial T cell heterogeneity in the invasive margin of colorectal cancer liver metastases associated with improved survival

ABSTRACT On a broader scale, T cell density and localization in colorectal cancer liver metastases have prognostic and predictive implications. As T cell distribution at higher resolutions has not been fully investigated, a detailed resolution analysis of T cell distribution was performed. Patient tissues were divided into 10 µm distance classes between the tumor border and adjacent normal liver. Thereby, distinct density patterns of T cell localization in relation to the malignant tissue could be detected. At a distance of 20 to 30 µm to the tumor, a decrease of CD3 T cells is common. Within this area, cytotoxic Granzyme B and CD8+ T cells were found to be significantly reduced as well as CD163 macrophages were increased and identified to be in close contact with T cells. Our data suggests a physical or functional border within this region. Survival analysis revealed improved overall survival in patients with high T cells numbers at the direct tumor border. Interestingly, the decreased T cells in the 20 to 30 µm region were also found to be significantly associated with improved survival. Consequently, the detailed localization of T cells, despite blockade, could be associated with improved clinical outcome. The high-resolution analysis represents new insights into relevant heterogenous T cell distributions especially related to clinical responses. As the paradoxical observation of localization-dependent prognostic relevance of T cell densities is only detectable by detailed spatial analyses, this investigation of spatial profiles at higher resolutions is suggested as a new biomarker for survival and response to therapies.


Introduction
Metastatic colorectal cancer (CRC) accounts for more than 55.000 deaths per year in the United States alone. 1 Around 15% of the patients with liver metastases (LM) can be cured by surgical resection of the metastases. 2 Immune system interactions with the primary tumor and the metastases have an important influence on the clinical course of patients. Immune cells infiltrating or surrounding colorectal primary tumors 3 and especially tumor-infiltrating lymphocytes (TILs) represent a major prognostic factor. 4,5 On a broader scale, high densities of TILs were shown to be correlated with an improved survival in patients with primary CRC as well as LM. [6][7][8][9][10][11][12][13][14] Nevertheless, tumor heterogeneity in many cancers represents an obstacle in detailed analyses of T cell distributions, whereas colorectal cancer liver metastases (CRC-LM) are rather homogenous with a clear border between liver and tumor tissue allowing more precise analyses of spatial distributions. 15,16 Especially, T cells and their relation to tumor tissue are gaining increased attention, as the majority of the T cells can be found in the stromal part of primary CRC without direct contact to tumor epithelium. 17 Evaluation of LM revealed heterogeneous immune cell infiltration into the metastasis center dependent of histopathological tumor conditions like necrosis or fibrosis. Previous analyses, however, did not reveal a relationship of these lymphocytic infiltrates and the clinical course. 17 In contrast, in the invasive margin, immune cell infiltrations showed remarkably strong variations between different patients and allowed prediction of responses to chemotherapy. 13,17 While the only few T cells in contact with the tumor or within the epithelial compartment are generally regarded as being engaged in direct antitumor activity, more distant T cells might be blocked by immunosuppressive factors and exploited by the tumor. [18][19][20] It is therefore of high interest to analyze differences in this distribution. The spatial distribution of T cells at the invasive margin of colorectal cancer LM has not been investigated in detail. The aim of this study was to explore the pattern of detailed T cell distribution in the invasive margin of CRC-LM patients using automated high-throughput image quantification analyses on whole slide sections samples. This quantitative approach allows the identification of distribution patterns in 10 mm steps within 100 mm distance from the tumor border and represents a high-resolution map of tumor infiltrating cells. By using serial sections, T cell distance analysis not only to the tumor tissue but also to the existing barriers like immunosuppressive cells is possible. Surface tissue areas of 36 patients were analyzed, allowing robust statistical evaluation. Clustering analyses revealed distinct distribution patterns that could be associated with improved overall survival (OS).

Results
The classification of the 10 mm distance classes within the 100 mm area between the tumor and adjacent liver as outlined in Fig. 1 and the following cell quantification analyses could be performed successfully for all 36 patients. Thereby, bar plots for all distance classes were generated. This revealed distinct patterns that were detectable within each patient. For CD3, T cell quantification representative results of different patients are shown in Fig. S1. Apparently, CD3 T cells are located in all distance classes with maximum densities at different distances.
As shown in Fig. 2, the distribution of CD3 cell densities obviously varied between the different distance classes (Kruskal-Wallis test, p D 0.0001). The statistically significant differences between adjacent distance classes are marked by asterisks (post-hoc Mann-Whitney U tests). It was observed that in the 20-30 mm distance class median T cell densities significantly decreased below 100 cells/mm 2 , whereas in the distance classes closer than 20 mm and further away than 30 mm from the tumor epithelium median densities higher than 150 cells/mm 2 were found. To identify patterns of T cell densities, unsupervised clustering was performed. To avoid possible bias due to different lengths of the invasive margin between different patients, we randomly selected one single field from each patient for further pattern analysis. The clustering analysis revealed distinct groups as shown in Fig. 3 (non-informative distance classes were removed stepwise, informative distance classes shown). The adjacent distance classes 10-20 mm and 20-30 mm clustered, whereas the most proximal distance classes (< 10 mm) and the classes with the highest distance from the tumor epithelium (90-100 mm) were found to be least related to the 20-30 mm distance classes. Furthermore, clustering of patients according to their T cell densities within single distance classes was observed. Survival analysis comparing patients with high and low CD3 T cell densities within each distance class resulted in a significant prolongation of OS for patients with high CD3 cell numbers in the most proximal distance class (< 10 mm), while in the 20-30 mm region a significant prolongation of survival for patients with low T cell numbers was observed (Fig. 4). No significant differences between high and low T cell infiltrates were found in all other distance classes.
To get a more detailed picture of the CD3 cell distribution pattern, we did a similar distribution and survival analysis (only CD8 C ) in representative patient samples for CD8 C (n D 10), Granzyme B (n D 6), CD163 (n D 4) and CD20 (n D 4) as shown in Fig. 5 (and Fig. S2). Especially, GranB was found to significantly "decrease" at the distance of 20-30 mm (Wilcoxon signed-rank test, p D 0.03) following a significant increase further away from the tumor (p D 0.03). In line with this, also CD8 C cells slightly but significantly decreased in the 20-30 mm region (p D 0.049). CD8 C cells, however, were rather evenly distributed among the whole area analyzed with a median cell number of »87 cells/mm 2 . Furthermore, survival analysis showed that high CD8 C infiltrates in the >10 mm region significantly elongated OS similar to CD3 (Fig. 6). Nevertheless, no significant differences between CD8 C high and low infiltration were observed in the other distance classes. In contrast to the T cells, CD163    macrophages showed a trend to increase in the 20-30 mm distance class (Fig. 5). While in the other distance classes CD163 cells were rather evenly distributed with a median of »200 cells/mm 2 , the distribution of CD20 cells was found to be more heterogenous with least cells at the tumor border, and slightly increasing densities further away ( Fig. 5 and Fig. S2). From all cell types analyzed, GranB and CD20 cells appeared to be the smallest populations present in the tissues with median cell counts of 9 and 16 cells/mm 2 , respectively. In Figs. S3 and 4, representative images of CD3, CD8 C , GranB, CD163 and CD20 immunostainings are shown.
In the course of this analysis, we further wanted to understand whether CD163 cells achieve immunosuppressive effects against T cells. We hypothesized that both cell types might be in close contact with each other and therefore analyzed the localization of T cells in relation to CD163 macrophages using serial section analysis (Fig. 7). We identified that 52-81% (median 68%) of CD3 Figure 4. Kaplan-Meier plots for estimated overall survival probabilities of CD3 T cell high and CD3 T cell low patient groups within all distance classes. Graphs indicate cumulative survival (y-axis) and survival in months (x-axis). p values for statistically significant differences between groups are shown (Breslow test). Cutoff values are taken from previous works. 6,13 lymphocytes within the tumor microenvironment were in close contact (< 10 mm) with CD163 macrophages. In addition, immunofluorescent double staining of CD163 and programmed death-ligand 1 (PD-L1) showed that the majority of CD163 macrophages were PD-L1 positive (Fig. 8). Figure 6. Kaplan-Meier plots for estimated overall survival probabilities of CD8 C high and low patient groups within the first three distance classes. Graphs indicate cumulative survival (y-axis) and survival in months (x-axis). p values for statistically significant differences between groups are shown (Breslow test). Cutoff values are taken from previous works. 6,13 Figure 5. Boxplots of identified CD8 C , GranB, CD163 and CD20 cell densities in different distance classes. Statistically significant differences between adjacent distance classes are marked with asterisks (Wilcoxon signed-rank test; blue p values given for paired t-tests; patient numbers for CD8 C n D 10, GranB n D 6, CD163 n D 4 and CD20 n D 4).

Discussion
Image analysis is a powerful tool to analyze spatial distributions as well as patterns of cells in the tumor microenvironment 25,26 and allows in combination with virtual microscopy the automated analysis of large tissue areas. Among the inflammatory cells in the composition of the tumor stroma and epithelium, T cells play a crucial role in regulating tumor growth and progression. It was shown that the invasive margin and its immune cell composition are predictive for the clinical course of cancer, 6,13,17,27 and that immunophenotypes and distribution of those can change under therapy. 28 Whereas the sheer quantity of T cells is associated with clinical outcome, 6-9,13,29 it has not yet been established whether there are any significant differences in the spatial distribution of T cells at the invasive margin. For this reason, we have investigated the detailed localization of CD3eC T cells in the invasive margin of colorectal cancer LM. CD3e is a pan-T cell marker and shows great heterogeneity in distribution intratumorally and at the invasive margin. 17 The interface between metastatic liver lesions and adjacent liver allows the precise evaluation of T cell localization, while in the primary tumor the heterogeneity is much more pronounced and is therefore difficult for this analysis. 15 The analyzed regions within 100 mm distance from the tumor show a significant "drop" of T cell density in the region of 20-30 mm from tumor epithelium, whereas higher T cell densities are found in close proximity and more far away from the tumor as schematically shown in Fig. 9. Similar examination of CD8 C cells revealed that this CD3 T cells subset is likewise affected by showing a slight but significant "drop" at 20-30 mm from the tumor. The ratio between CD8 C and CD4 C cells varied across distance classes. In general, the CD8 C cell numbers observed among the complete analyzed regions from 0 to 100 mm were »45% (87 cells/mm 2 ) of the CD3 cell numbers (195 cells/mm 2 ) leading to the assumption that the other major CD3 cell subtype (CD4 C ) is evenly present in this area. This is, however, not the case at the location of the T cell "drop" in the 20-30 mm region. Here, a median of »90 CD3 cells/mm 2 and »67 CD8 C cells/mm 2 was counted, meaning that CD4 C cells are more affected to decrease in this region. Interestingly, when examining spatial cytotoxic activity by GranB immunostaining it was likewise found to be significantly   decreased at the 20-30 mm distance. As natural killer cells have been shown to be very rare in CRC-LM, cytotoxic CD8 C T cells are supposed to be the main source of GranB. 30 These findings strikingly confirm not only a reduction of cell quantities but also a reduction of cytotoxic T cell function at the 20-30 mm region.
Those new observations prompted us to further investigate the possible factors involved in this heterogenous T cell distribution patterns. As regulatory T cell numbers are generally low in CRC-LM (3.5 cells/mm 2 ), 13 we decided to repeat our spatial distribution analysis to examine the localization of CD20 B cells and CD163 macrophages. This latter tumor-associated macrophage (TAM) population, which is characterized by a tumorpromoting and immunosuppressive phenotype and associated with poor clinical outcome in many different tumors, 31-33 we found highly enriched in CRC-LM tissues (»200 cells/mm 2 ). In our exploratory analyses, we noticed an increase (p D 0.068) at the 20-30 mm region, whereas no relevant differences between the other examined distances classes were observed. This opposite distribution to the T cells leads to the assumption that CD163 cells fill the gap between the tumor distant and tumor proximal T cells. Our data suggests that CD163 macrophages are forming a physical barrier which only a specific T cell subset can overcome. Even more important, this barrier could have functional impact, because it also marks the border, where cytotoxic GranB significantly decreased at the tumor margin. Indeed, TAM have been shown to directly suppress T cell responses through immune inhibitory checkpoint molecules like PD-L1. 34,35 Interestingly, these PD-L1 positive TAM were also found to be in close contact with T cells. This data goes in line with our own results showing a high percentage (68%) of T cells being in direct contact with CD163 macrophages, and furthermore, the majority of CD163 cells in the invasive margin were found to express PD-L1. Our data, therefore, suggests itself that TAM directly suppress T cells and consequently GranB production via immune checkpoints. The fact that in the invasive margin of CRC-LM 98% of all CD3 T cells express the checkpoint receptor PD-1 further encourages this assumption. 19 In contrast to CD163, the distribution of CD20 B cells is very heterogenous among the distance classes with lower cell densities at the tumor border and increased densities further away (p D 0.062). The role of tumor-infiltrating CD20 B cells is controversial. Some studies report about a protective role against tumors, [36][37][38] whereas others connect B cells with enhanced tumor growth and suppression of antitumor immunity. 36,39,40 Notably, B cells can polarize macrophages toward immunosuppressive phenotypes. 41 Nevertheless, B cell numbers appeared to be very low with »16 cells/mm 2 among the complete tissues analyzed and especially compared with CD163 (»200 cells/mm 2 ) a minor role for CD20 B cells is presumed. Besides immunosuppressive cells, the distribution of the T cells in the tissues could also be accounted for other factors. Chemotactic fields, which might be generated by tumor, endothelial cells or even macrophages, could induce the observed T cell distribution patterns. 19,42 Another explanation could be cellmatrix interactions. 43,44 Physical properties of the extracellular matrix can influence directional cell movement and might govern the positioning of T cells in the direct vicinity of tumor cells. In addition, CD31, which is a marker for angiogenesis and expressed on endothelial cells, also has a functional role in leukocyte migration. 45,46 We analyzed two CRC-LM patient tissues for detailed distribution of CD31 (Fig. S5). Nevertheless, CD31 area was generally very high and showed strong heterogeneity among patients (»700-8,000 mm 2 ). In addition to that, CD31 can be expressed by leukocytes including T cells and thereby even having inhibitory effects on T cell activation 46,47 underscoring its highly diverse role in tumors. In this study, we performed a high-resolution analysis of different immune cell subsets. Consequently, the detailed functional mechanisms and cell interactions involved in the detected spatial T cell distribution patterns need to be investigated in future studies. The clear limitation of our study is the low sample number for some markers that limited our ability to identify statistically significant differences between the distance classes.
Whereas we could not see distinct patterns limited to single patients (data not shown), the clustering analysis revealed that specific groups of patients show patterns of CD3 T cell aggregates, which were limited to a given region (e.g., 0-10 mm or 20-30 mm, Fig. 3). Indeed, OS analysis showed that patients with high T cell densities (CD3 as well as CD8 C ) within the >10 mm distance class had a survival benefit compared with patients with low T cells in the aforementioned distance (Figs. 4-6). Interestingly, the opposite effect was found at the 20-30 mm distance from the tumor. Here, low CD3 T cell numbers had positive effects on patient survival. Consequently, despite the physical or functional T cells barrier, the observed pattern of T cell distribution at 20-30 mm distance obviously is beneficial for patients (under palliative therapy). When considering this observation from the other side, high T cell numbers unexpectedly might have a tumor promoting effect. Our own previous work illustrated the exploitation of CCL5-producing T cells in the invasive margin of CRC-LM patients by supporting CCL5 receptor (CCR5)-expressing tumor cells. 19 Furthermore, blockade of CCR5 led to repolarization of TAM resulting in beneficial clinical responses. These findings might explain, why especially in patients with low T cells in the 20-30 mm region a better survival is observed than in patients with high T cell numbers, and furthermore, it shows that reeducated macrophages, which were found to be increased at the 20-30 mm distance and which are in fact generally high in the invasive margin, can render antitumor effects under therapy. In addition to that, our results indicate that T cell induced antitumor effects can be conducted by T cells that have engaged the tumor margin at the 0-10 mm distance class, which is indeed also true for the cytotoxic CD8 C subtype, while T cells further away remain blocked. The contrary outcome of our detailed distribution analysis that both localization-dependent high as well as low T cell numbers have clinical benefit might explain startling observations for similar tumors, where in one study high 48 and in another study low T cell densities 49 indicate prognostic relevance. Our finding might further elucidate related controversial observations, where studies identify high densities of immunosuppressive cell populations being associated with favorable prognosis. 31,50,51 The reasons for the controversial observation of local T cell decrease being associated with favorable prognosis and if the observed patterns change under therapy are currently investigated. The paradoxical observation of localization-dependent prognostic relevance of T cell densities is only decipherable by detailed spatial analyses. Therefore, the analysis of spatial distribution patterns at higher resolution, and especially the distance classes <10 and 20-30 mm alone, could serve as more precise prognostic markers for responses to chemotherapy and OS compared with cell densities measured on a broader scale.
In summary, we present here for the first time detailed data on the local spatial distribution of immune cells in colorectal cancer LM. Thereby, distinct patterns of immune cells reveal a physical or functional T cell barrier at 20-30 mm distance from the tumor epithelium with decreased T cell and cytotoxic Gran B densities as well as increased CD163 macrophages. These specific distribution patterns with either high T cells at the direct tumor border or low T cells at the 20-30 mm distance are associated with improved OS. This detailed analysis of spatial profiles within the microenvironment, therefore, represents new insights into T cell distribution, function as well as influence of T cell localization on clinical responses and introduces possible spatial immunosuppressive hurdles that need to be overcome by successful immunotherapies. We suggest the analysis of spatial profiles within the < 10 and 20-30 mm distance from the tumor epithelium as new predictive biomarker for survival and response to therapies.

Tissue samples
The study was approved by the Medical Ethics Committee of the University of Heidelberg. Tumor samples were used after receiving a signed informed consent from all 36 patients. Patients were included in this study when diagnosed with stage UICC IV CRC after metastatic disease in a surgically resected liver metastasis was histologically proven. None of the patients received neoadjuvant treatment. Samples had to contain (a) metastasis, (b) invasive margin and (c) adjacent liver (Fig. 1A). Patients received standard of care palliative therapies and follow-up was obtained as described previously. 13 All samples were screened for microsatellite instability (MSI) using established protocols. 21 No MSI sample was included in this analysis. Clinical data as well as detailed pathological reports showed no evidence for the presence of ethyltoxic liver damage, chronic inflammation due to ongoing infections (e.g., hepatitis) or other chronic inflammatory conditions.

Evaluation of immunostainings and tissue classification
Using the NDP Nanozoomer system (Hamamatsu Photonics, Japan) complete microscopic images of full tissue sections were automatically obtained. Cell quantification and distance analysis were performed with the VisioMorph software (Visiopharm, Denmark). Image segmentation and subsequent classification according to color and morphology was established on a separate training set with manually counted cell numbers previously. 22 This approach in combination with serial cuts of tissue blocks allows robust and reproducible large scale histological evaluations with high precision across complete (adjacent) sections.
1 mm 2 regions of the invasive margin (Fig. 1A) were generated in an automated process, visually controlled and adapted correspondingly. Then, the positions of the T cells in relation to the tumor epithelium were identified in consecutive 10 mm steps (Fig. 1B). Distances to tumor epithelium were grouped into distinct classes (distance classes), e.g., "below 10 mm," "10-20 mm," "20-30 mm" and so on. So for each tissue section, T cell quantities in 10 distance classes within 100 mm from the tumor epithelium into the normal adjacent liver were established. Immune cells and other markers were quantified in this process by spectral and morphological criteria as described before 15,22 in combination with aggregation of data from the spatial registry.
To evaluate the distance between CD163 macrophages and T cells within the tumor microenvironment, serial sections were stained with the respective markers and after scanning of the sections, these were computationally analyzed using Visio-Morph software (Fig. 7).

Statistical analyses
Statistical analyses were performed using the SPSS 18.0 and 22 packages (IBM) as well as GraphPad Prism version 5. To discriminate CD3 cell counts between distance classes Kruskal-Wallis tests with post-hoc Mann-Whitney U tests were performed. Results with two-tailed p values <0.01 were judged to be statistically significant. Boxplots depict the minimum, first quartile, median, third quartile and maximum. Measurements with values between 1.5 and 3 lengths of the box away from the upper and lower border of the box were classified as outliers (white circles). The length of the box equals the interquartile range. To discriminate cell counts between single distance classes for CD8 C and GranB Wilcoxon signed-rank and for CD20 and CD163 paired t-tests were performed. Results with two-tailed p values <0.05 were judged to be statistically significant. Unsupervised hierarchical clustering analysis using Euclidean distance was performed with the R software package. 23 OS curves were computed using the Kaplan-Meier estimator (Breslow test) and significance was assumed for p values of 0.05 or less. Cutoff values were adapted (quantile-quantile class) from previous work, 13 which is in line with findings from primary CRC tumors. 24

Disclosure of potential conflicts of interest
No potential conflicts of interest were disclosed.