Clinical Research Papers:

Tumors with high-density tumor infiltrating lymphocytes constitute a favorable entity in breast cancer: a pooled analysis of four prospective adjuvant trials

PDF |  HTML  |  Supplementary Files  |  How to cite  |  Order a Reprint

Oncotarget. 2016; 7:5074-5087. https://doi.org/10.18632/oncotarget.6231

Metrics: PDF 1839 views  |   HTML 1942 views  |   ?  

Vassiliki Kotoula _, Kyriakos Chatzopoulos, Sotiris Lakis, Zoi Alexopoulou, Eleni Timotheadou, Flora Zagouri, George Pentheroudakis, Helen Gogas, Eleni Galani, Ioannis Efstratiou, Thomas Zaramboukas, Angelos Koutras, Gerasimos Aravantinos, Epaminontas Samantas, Amanda Psyrri, Helen Kourea, Mattheos Bobos, Pavlos Papakostas, Paris Kosmidis, Dimitrios Pectasides and George Fountzilas


Vassiliki Kotoula1,2, Kyriakos Chatzopoulos2, Sotiris Lakis2, Zoi Alexopoulou3, Eleni Timotheadou4, Flora Zagouri5, George Pentheroudakis6, Helen Gogas7, Eleni Galani8, Ioannis Efstratiou9, Thomas Zaramboukas9, Angelos Koutras10, Gerasimos Aravantinos11, Epaminontas Samantas12, Amanda Psyrri13, Helen Kourea14, Mattheos Bobos2, Pavlos Papakostas15, Paris Kosmidis16, Dimitrios Pectasides17 and George Fountzilas2

1 Department of Pathology, Aristotle University of Thessaloniki School of Medicine, Thessaloniki, Greece

2 Laboratory of Molecular Oncology, Hellenic Foundation for Cancer Research, Thessaloniki, Greece

3 Health Data Specialists Ltd, Department of Biostatistics, Athens, Greece

4 Department of Medical Oncology, “Papageorgiou” Hospital , Thessaloniki, Greece

5 Department of Clinical Therapeutics, “Alexandra” Hospital, Athens, Greece

6 Department of Medical Oncology, Ioannina University Hospital, Ioannina, Greece

7 First Department of Medicine, “Laiko” General Hospital, Athens, Greece

8 Second Department of Medical Oncology, “Metropolitan” Hospital, Piraeus, Greece

9 Department of Pathology, “Papageorgiou” Hospital, Thessaloniki, Greece

10 Division of Oncology, Department of Medicine, University Hospital, Patras, Greece

11 Second Department of Medical Oncology, “Agii Anargiri” Cancer Hospital, Athens, Greece

12 Third Department of Medical Oncology, “Agii Anargiri” Cancer Hospital, Athens, Greece

13 Division of Oncology, Second Department of Internal Medicine, Attikon University Hospital, Athens, Greece

14 Department of Pathology, University Hospital of Patras, Rion, Greece

15 Oncology Unit, “Hippokration” Hospital, Athens, Greece

16 Second Department of Medical Oncology, Hygeia Hospital, Athens, Greece

17 Oncology Section, Second Department of Internal Medicine, “Hippokration” Hospital, Athens, Greece

Correspondence to:

Vassiliki Kotoula, email:

Keywords: breast cancer, tumor infiltrating lymphocytes, clinical breast cancer subtypes, trastuzumab, prognostic

Received: April 15, 2015 Accepted: October 09, 2015 Published: October 25, 2015


Background: Tumor infiltrating lymphocytes (TILs) are considered in the prognosis of breast cancer (BC) patients. Here, we investigated the prognostic/predictive effect of TILs in patients treated in the frame of four prospective trials with adjuvant anthracycline-based chemotherapy in the pre- and post-trastuzumab era.

Methods: TILs density was histologically assessed as percentage of stromal area on whole routine sections of 2613 BC (1563 Luminal A/B; 477 Luminal HER2; 246 HER2-enriched; 327 triple negative [TNBC]) and were evaluated as high/low at three cut-offs (c/o; 50% [lymphocytic predominance, LP], 35% and 25%), in separate training and validation sets.

Results: High TILs were present in 3.5%, 6.5% and 11.5% of all tumors, using the 50%, 35% and 25% c/o, respectively. TILs status did not interact with BC subtypes or trastuzumab treatment. LPBC patient outcome was not affected by nodal status, while high TILs were favorable in TNBC with unfavorable nodal status. When adjusted for standard clinicopathological parameters and treatment, high TILs independently predicted for favorable outcome, e.g., disease-free survival with the 35% c/o in the entire cohort (HR = 0.44, 95% CI 0.28-0.69, p < 0.001) and in specific subtypes.

Conclusions: High TILs tumors, especially LPBC seem worthy validating as a separate entity of favorable prognosis in breast cancer.


Tumor infiltrating lymphocytes (TILs) within the tumor stroma or within tumor nests reflect the host immune response against the tumor, described as cancer “immunoediting” [1]. Measurable tumors are in the “escape” phase during which immune cells are not able to eliminate tumor growth but their presence at least denotes stand-by immunocompetency [1-3] that can be reactivated by treatment. Indeed, in a per se non-immunogenic environment, such as the breast, the presence of stromal and/or intratumoral TILs seems to fully support that breast cancer (BC) is immunogenic, especially concerning estrogen receptor (ER) negative disease [3, 4].

Recent studies have shown that TILs presence and density are favorable prognosticators in breast cancer either with [4-11] or without [12] subtype specificity. Similar TILs effects were noticed in the neoadjuvant setting [4, 6, 11, 12] and in the adjuvant setting as well [5, 7, 8, 11], while the presence of TILs in residual disease after initial chemotherapy may also indicate favorable outcome [13]. The presence and density of TILs [7] and the expression signatures of immune function genes [14] may predict benefit from trastuzumab in HER2-positive patients, while TILs presence may also be predictive for benefit from cytotoxic drugs like docetaxel [8] and carboplatin [6].

Methods for TILs assessment in breast cancer greatly vary in different studies, ranging from simple mononuclear cell infiltrate counting on routine hematoxylin & eosin (H&E) stained slides [7, 8, 15]; immune cell typing with immunohistochemistry [11, 12, 15]; and, immune-related gene expression [4, 9, 14] to digital immune cell weighing [16]. In an effort to bring TILs closer to clinical application, recommendations for the morphological assessment of TILs on H&E sections have also been recently published [17].

In the present pooled analysis, we investigated the effect of TILs on the outcome of more than 2500 patients with operable breast cancer who were treated in the setting of prospective trials conducted by the Hellenic Cooperative Oncology Group (HeCOG). The study population allowed for the evaluation of TILs as prognostic parameter upon anthracycline - taxane regimens and as predictive for trastuzumab benefit in HER2 positive patients. TILs were morphologically assessed as a continuous variable but were finally analyzed at three different cut-offs according to breast cancer clinical subtypes and to nodal status.


The distribution of TILs rates is shown in Figure 1A. The mean stromal TILs density in all 2618 tumors was 12.7 (±SD 15.5) and the median was 7. In 1290 tumors (49.3% of total), stromal TILs rate was 0 - 5%, while additional 587 tumors (22.4%) had 10 - 20% stromal TILs.

Among all tumors, 91 (3.5%) were considered as LPBC since they had ≥50% stromal TILs; 170 (6.5%) had ≥35%, and 301 (11.5%) had ≥25% TILs. In LPBC and in tumors with TILs close to 0 (Figure 1B-1D), these infiltrates were homogeneously dense or absent in the same tumor, respectively. In comparison, TILs distribution was heterogeneous in cases with 5 - 35%, i.e., TILs density differed by >10% in different areas in the same section or among sections from the same tumor (Figure 1E & 1F).

Stromal TILs in breast carcinomas.

Figure 1: Stromal TILs in breast carcinomas. A. Distribution of TILs rates in the study population. B. < 5% (Bar: 10um); C. 50% (Bar: 100um); D. 75% (Bar: 100um) TILs homogeneously dispersed in these tumors. E and F: Same case, overall 35% TILs, heterogeneously populating the stroma of this tumor. E at 10%, F at 40%. Bars at 10um.

TILs and clinicopathological characteristics

Higher stromal TILs rates, either as continuous or binary variables as shown in Figure 2A, were significantly more frequent in the absence of ER/PgR, in HER2 positive, especially in HER2-enriched, in basal-like and in TNBC tumors (all p’s < 0.0001). For example, with the 35% cut-off, TNBC (13.6%), HER2-enriched (11.9%) and Luminal HER2 (9.3%) had significantly more often high TILs as compared to Luminal A/B (3.3%) tumors. Importantly though, Luminal A/B, HER2-positive and TNBC were represented at similar rates within the group of tumors with high TILs (Figure 2B).

Tumors with high TILs also had statistically significantly higher Ki67 (Mann-Whitney p < 0.001) and were more often of grade III (Pearson’s chi-square p < 0.001). However, Ki67 labeling values largely overlapped between high and low TILs categories, rendering the biological significance of this finding questionable. TILs density was not associated with nodal status and tumor size in our series; it also did not differ in comparison to age and menopausal status (Table S5).

With respect to specific histological breast cancer types that were adequately represented in the studied cohort for statistical comparisons (>15 cases per subtype), high TILs were present in the majority of the 24 medullary carcinomas (67%, 75%, and 96% with the 50%, 35% and 25% TILs cut-offs, respectively), while these were absent or very rare among the 221 lobular, 30 apocrine, and 17 mucinous carcinomas (Figure 2C).

TILs association with IHC breast cancer phenotypes and with specific histological types.

Figure 2: TILs association with IHC breast cancer phenotypes and with specific histological types. A. Numbers in bars indicate the actual population per category. High TILs bars (blue, green, violet) are shown as percentage per phenotype with all examined cut-offs, as indicated. A prevalence of high TILs in the non-ER/PgR context is evident, where the incidence of tumors with ≥50% (lymphocyte predominant, LPBC) was minimal. B. By contrast, the incidence of major subtypes within the group of high TILs tumors was not substantially different (chi square p = 0.87). C. Numbers in the Y-axis represent the actual population per category. The majority of high TILs tumors were of the non-specific type (NST). However, the rate of high TILs tumors among the 24 medullary carcinomas was the highest among all subtypes. Lobular, apocrine and mucinous carcinomas seldom exhibited high-TILs.

Favorable prognostic high TILs in breast cancer

Among all patients, even without distinguishing for subtype and trastuzumab treatment, high TILs at all examined cut-offs were associated with longer DFS as compared to low TILs (Figure 3). Among the 91 patients with LPBC (TILs ≥50%), only 6 (6.6%) relapsed during a period of more than 10yrs follow-up and out of them only 2 (2.2% of all LPBC) during the first 3 years (Table 1). None or rare events were noticed for these patients for overall survival (OS) as well (Table S6). Risk for relapse and death was remarkably constant in the training and validation sets; however, due to the small number of patients and events (0-2 patients with high TILs tumors in subtype categories), results in the validation set appeared as statistically non-significant. Similar, albeit less pronounced results were obtained for favorable high TILs at the 35% and 25% cut-offs in the entire patient cohort (Table 1).

Table 1: Univariable analyses (Log-rank and Cox) showing the effect of high TILs at all 3 cut-offs in the training and validation sets with respect to patient disease-free survival

TILs association with patient disease-free survival (DFS) in 2613 breast cancer patients.

Figure 3: TILs association with patient disease-free survival (DFS) in 2613 breast cancer patients. High TILs at all cut-offs examined conferred longer DFS.

As can be retrieved from Table 1, the rate of relapses for high TILs at any cut-off was similarly low among subtypes. However, when comparing the impact of TILs density with respect to each subtype, no statistically significant effect was observed for high TILs on the outcome of patients with Luminal A/B tumors. Still, such patients with high TILs, for example at the 35% cut-off, who did not relapse within the first 5 yrs (46 out of 52) remained relapse-free during the entire follow-up period (Figure 4A, Table 1). High TILs at all cut-offs conferred decreased risk for relapse and death in patients with HER2 positive tumors irrespectively of trastuzumab treatment (Figure 4B). This effect was particularly strong in LPBC where only 3/37 patients relapsed and none died within the first 3 years from diagnosis, as compared to 101/686 relapses and 30 deaths that were observed for non-LPBC HER2-positive patients during the same period (Table 1). Similarly, favorable high TILs effects at all cut-offs, again more pronounced for LPBC were observed for TNBC patients (Figure 4C). Only 1 out of 27 TNBC patients with LPBC relapsed and none died during the entire follow-up period, as compared to 85 and 71 out of the 300 non-LPBC TNBC patients who relapsed and died, respectively (Table 1).

In the entire cohort and with respect to clinical subtypes, high nodal burden (≥4 infiltrated lymph nodes) strongly predicted for shorter DFS and OS as compared to low nodal burden (0-3 nodes) (Table S7). High nodal burden was universally unfavorable for patients with low TILs tumors at all cut-offs (Table S8a for DFS; Table S8b for OS). In all settings, nodal status did not affect the outcome of LPBC patients, for whom 0 to 1 relapses were observed per subtype group. For patients with >35% and >25% TILs, high nodal burden remained unfavorable in the Luminal A/B and HER2-positive groups. Of note, the rare events in the favorable high TILs groups always concerned the same Luminal A/B and HER2-positive patients, which probably biased the statistical significance of the interaction for the 25% cut-off. In TNBC patients, high TILs conferred significantly better outcome for unfavorable nodal status (Figure 4D), other than in the Luminal A/B and HER2-positive groups. Again, however, the numbers of events were very small and no statistical significance was reached for this interaction.

Histologic subtypes were not separately evaluated for TILs effect on outcome, due to small numbers in each category.

Effect of high TILs on patient DFS with respect to clinical breast cancer subtypes.

Figure 4: Effect of high TILs on patient DFS with respect to clinical breast cancer subtypes. Results are shown for high TILs at the 35% cut-off. In all subtype categories A. - C. high TILs were associated with longer DFS. In A. although high TILs are not considered prognostic in Luminal A/B disease, patients with such tumors who did not relapse within the first 5 years remained relapse free for more than 15 years. In B. HER2-positive patients were treated with and without trastuzumab. D. Patients with unfavourable nodal status and high TILs (yellow star) fared equally well as patients with favourable nodal status, where TILs levels were not associated with outcome.

High TILs did not significantly predict for trastuzumab benefit

The effect of TILs in HER2 positive disease was also examined as an interaction with trastuzumab, by comparing the outcome of patients treated in the pre- and post-trastuzumab era. Among all HER2-positive patients, a significant trastuzumab benefit was noticed for those with low TILs. Thus, trastuzumab significantly benefited patients with < 35% TILs (HR: 0.37, 95%CI: 0.25-0.53), with < 25% TILs (HR: 0.39, 95%CI: 0.27-0.57), and, with < 50% TILs (HR: 0.36, 95%CI 0.25-0.52). In comparison, trastuzumab did not offer statistically significant outcome advantage in patients with high TILs tumors; with or without the drug, patients with high TILs fared better than those with low TILs (Table 2 and Figure 5 for the 35% cut-off; Table S9 [DFS] and Table S10 [OS] for LPBC and for the 25% cut-off). No interaction between trastuzumab and TILs was found among all HER2-positive patients (Figure 5A). Similar results were obtained for all TILs cut-offs, in the entire group of HER2-positive patients and separately in the test and validation sets for DFS (Table S9) and OS (Table S10). However, at all comparisons, HER2-positive patients treated with trastuzumab for tumors with high TILs density fared best. For example, none of 18 trastuzumab-treated patients with LPBC and only 2/61 patients with >25% TILs relapsed over a period of 5 years, as compared to 2/19 and 18/66 non-treated patients with respective TILs status (Table S9). In addition, no interaction was observed between TILs density and HER2-positive subtypes (Figure 5B and 5C), despite that trastuzumab treated patients with high TILs fared best with few relapses during the entire follow-up period. Relapse rate was proportionally lower for patients with high TILs Luminal-HER2 as compared to those with HER2-Enriched, but for both subtypes high TILs were favorable as compared to low TILs. Again, statistical results were probably biased by the small number of patients and events for each HER2-positive subtype with high TILs (Table 2).

Finally, no interaction was observed between TILs density, disease stage I vs. II, and nodal status with respect to trastuzumab treatment.

Table 2: Interaction tests for TILs, T treatment and HER2 positive subtypes on patient DFS

Effect of TILs with respect to trastuzumab (T) in HER2-positive patients.

Figure 5: Effect of TILs with respect to trastuzumab (T) in HER2-positive patients. Results are shown for high TILs at the 35% cut-off. Whether statistically significant or not, high TILs were favourable in every context examined (all HER2-positive and HER2-positive subtypes). Trastuzumab significantly benefitted patients with low TILs but this effect, although present, was insignificant for patients with high TILs.

High TILs as an independent favorable prognostic marker in breast cancer

Multivariable models including standard clinicopathological parameters affecting patient outcome, as described in the Methods section and in Table S7, revealed the independent favorable prognostic effect of high TILs in breast cancer patients treated with adjuvant anthracyclines and taxanes, as shown for DFS in Figure 6 and for OS in Table S11. In all models examined for each subtype separately, high TILs retained their favorable independent prognostic effect. When examined in the entire patient population, high TILs were the strongest favorable prognosticator in the examined clinical context.

Forest plot showing the strongly significant independent favourable effect of high TILs on patient DFS.

Figure 6: Forest plot showing the strongly significant independent favourable effect of high TILs on patient DFS. Multivariable models for the entire cohort and for breast cancer subtypes are shown. High TILs, here shown according to the 35% cut-off, had a constantly favourable effect on patient outcome in all settings (green circles).


This study confirmed the presence of high TILs as a robust favorable prognosticator in high-risk operable breast cancer, an effect that was more pronounced in HER2 positive and TNBC patients, in line with previous reports in the adjuvant [5, 7-9, 11, 24] and neoadjuvant setting [6, 9, 11, 13]. Herein, instead of using continuous TILs counts and reporting on added risk discounts, TILs were evaluated in a binary mode as high/low according to predefined cut-offs, all of which were validated in a second study set and in the entire population. These data and the fact that the patients studied represented a pooled cohort from different prospective trials that were conducted over a period spanning more than 10 years of clinical practice, support the robustness of morphologically assessed high TILs as a favorable prognostic marker in operable high-risk breast cancer, definitely in tumors with unfavorable subtypes.

In the better prognosis Luminal A/B tumors, high TILs may also offer DFS advantage; however, even with the present long follow-up, statistics were underpowered for the evaluation of high TILs in this group of patients. High TILs tumors, even at the lowest cut-off we applied, are an infrequent occasion in breast cancer, with lowest rates in Luminal A/B; the latter constitutes an additional reason for low statistical power and for inconclusive results in this group of tumors in most studies, as also recently noticed [24]. However, if we consider high TILs tumors as a separate entity the following apply: (a) major subtypes, i.e., Luminal A/B, HER2 positive irrespective of ER/PgR, and, TNBC are evenly distributed within this entity; (b) the rate of events per subtype within the entity is very similar; (c) the number of events is extremely limited for cases with LPBC. By setting the cut-off for calling LPBC at 50% TILs density, corresponding patients of any subtype had an excellent course of disease, as previously described for HER2 positive and TNBC [7, 9]. In our series, about 10% of patients with non-Luminal A/B tumors and less than 2% of patients with Luminal A/B had LPBC. Due to the high worldwide incidence of breast cancer, however, the corresponding absolute number of patients is not negligible. Our data support addressing LPBC as a distinct breast cancer subtype to be considered for alternative therapeutic interventions.

In the presented series, the cohort specific cut-off for distinguishing high/low TILs was 35%, while the favorable TILs effect statistically disappeared below 25% average density. This gradual declining of prognostic TILs significance supports the notion to evaluate TILs as a continuous parameter [17], which, however, may be impractical for the application of this marker in diagnostics. Although the need to consistently report TILs still remains [25], based on our findings, a three-scale grading system for reporting TILs status seems applicable, i.e., LPBC (high and homogeneous), non-LPBC with high TILs called at 25% as the lower cut-off, and non-LPBC with TILs below that level.

Of note, a parameter that may also need to be assessed along with TILs density is the heterogeneous distribution of TILs in the stroma of breast carcinomas. As shown here, apart from LPBC and from tumors with close to null immune infiltrates, TILs density is largely heterogeneous within the same tumor, which is challenging to report in routine diagnostics. Assessing the degree or patterns of spatial intratumoral TILs heterogeneity might also be an option for a more detailed description of the immune response status of the tumor, as recently suggested for ER/PgR negative breast cancer [16]. In order to approach the heterogeneous TILs distribution, we may consider anti-tumor immune response being triggered by neo-antigens produced by the tumor, which is the result of the underlying genetic alterations [26]. Following this concept, the observed spatial heterogeneity in TILs density may reflect the well established genetic heterogeneity in breast cancer. Although the issue is tremendously important, such data from clinical material are still missing. Immune response to neoantigens is likely a universal process, not confined to specific breast cancer subtypes or to clinically actionable targets, at least not in the way these parameters are currently considered. The fact that the majority of high TILs tumors are found among TNBC and HER2-positive cases may just reflect the extensive genomic alterations described for these tumors as compared to Luminal A/B [27]. Since, as described here, heterogeneous TILs density is observed in non-LPBC, it seems worthy including TILs heterogeneity as a binary parameter (presence vs. absence) for the evaluation of non-LPBC in future studies.

With respect to HER2-positive tumors, we did not observe any statistically significant interaction between TILs and trastuzumab treatment. Trastuzumab significantly benefitted patients with low TILs as determined here; patients with high TILs fared better with or without this drug. This appears in line with the findings by Perez et al for the N9831 trial [28] but in contrast to the previously reported predictive role of TILs for trastuzumab benefit in the FinHER trial [7]. In this trial, HER2-positive LPBC without trastuzumab fared worse than non-LPBC, which was not the case in our series from the pre-trastuzumab era involving a larger number of patients; this may explain the present discrepant statistical result with respect to the predictive role of TILs for trastuzumab benefit. Further, the same prognostic effect of high TILs was observed with all cut-offs used herein. However, although not statistically significant, patients with high TILs who were treated with trastuzumab fared best over a period of 5 years. Overall, these results support the prognostic role of high TILs in HER2-positive disease independently of trastuzumab treatment and underline the necessity for larger patient series and follow-up longer than 5 years for the evaluation of the predictive value of TILs with respect to trastuzumab.

A novel piece of data in the present study concerns the impact of nodal status in association with TILs density on patient outcome. Universal findings for all subtypes at all TILs cut-offs were (a) that patients with LPBC fared best irrespectively of nodal status, which further indicates LPBC as a distinct entity; and (b) that the aggravating prognostic impact of low vs. high nodal burden was effective in patients with low TILs tumors, which was constant in all comparisons. Inefficiency of the host to block tumor expansion may only partially explain this condition, since, for example, patients with high TILs TNBC and unfavorable nodes treated with adjuvant chemotherapy fared as well as patients with favorable nodes. As reported, immune changes in tumor draining lymph nodes may not be solely driven by local tumor invasion [29] but the biological basis of the role of intratumoral TILs in the behavior of locally spread disease has not yet been investigated. Nevertheless, it remains challenging and clinically important to dissect patient prognosis beyond the classic parameter of nodal status. If the present findings are validated in independent studies TILs density may be helpful in this context.

In conclusion, the presented pooled analysis supports the need for the histological assessment of TILs density as a major marker of outcome along with nodal status in breast cancer patients treated in the adjuvant setting. No clear predictive role was revealed for trastuzumab benefit with this marker but it is useful to know that trastuzumab-treated patients with high TILs fare best. Distinguishing LPBC as a distinct entity with excellent prognosis in the adjuvant setting appears clinically important for treatment decisions, while using high TILs cut-offs may identify patients who would benefit from agents recharging their anti-tumor immune response. The present findings merit validation in independent large cohorts by taking into account the low rate of high TILs tumors in patients with operable high-risk breast cancer.


Patients and tumors

Tumor tissue material from 2618 patients who had been diagnosed between 1997 and 2010 with operable breast cancer and had been treated with adjuvant chemotherapy (anthracyclines - taxanes) in the setting of four prospective clinical trials by the Hellenic Cooperative Oncology Group (HeCOG) was examined (HE 10/97 [18]; HE 10/00 [19]; HE 10/05 [20]; HE 10/08 [manuscript in preparation]). The basic trial characteristics are shown in Table S1. In HE10/05 and HE10/08 trastuzumab was administered sequentially for one year after the completion of chemotherapy. Patients had provided written consent for the use of their biologic material for research purposes and the study was approved by the Bioethics Committee of the Aristotle University of Thessaloniki School of Medicine (#77/10June2014) and by the Institutional Review Board of Papageorgiou Hospital of Thessaloniki (#725/10May2013). Paraffin blocks were collected retrospectively for HE10/97 and prospectively for the other 3 trials. The distribution of patients and tumors per clinical study and basic demographic and clinicopathological characteristics are shown in Table S2.

Tumors had been routinely diagnosed in local pathology labs, where they were also evaluated for ER/PgR/HER2 protein expression with immunohistochemistry (IHC) for patient stratification to receive hormone treatment and trastuzumab. HER2 FISH or CISH had also been applied locally in cases with ambiguous HER2 IHC for patients in trials HE 10/05 and HE 10/08. HER2-positive patients were treated with trastuzumab in these two trials (post-trastuzumab era) but not in HE 10/97 and HE 10/00 (pre-trastuzumab era). For the purposes of the present study, tumors were evaluated for combined ER/PgR positivity in the absence of HER2 protein overexpression and/or gene amplification as Luminal A/B tumors; Luminal-HER2 if ER/PgR positive and HER2 over-expressing or amplified; HER2-enriched if ER/PgR negative but HER2 pathology positive; and, as triple negative (TNBC) if ER/PgR/HER2 negative.

All tumors were also retrospectively centrally subtyped in the Laboratory of Molecular Oncology (HeFCR/HeCOG/AUTH, Thessaloniki, GR) with ER/PgR/HER2/Ki67 IHC and HER2 FISH, as previously described [21]. In addition, cytokeratin-5 (CK5) and epidermal growth factor receptor (EGFR) IHC were applied for typing basal-like carcinomas with 1% cut-off for positivity, as suggested [22].

TILs evaluation

Mononuclear infiltrates corresponding to TILs were assessed on whole routine hematoxylin and eosin (H&E) sections of 2613 breast carcinomas, by one pathologist (K.Ch.) as % of stromal area according to the recently published recommendation [17]. In 5 cases it was impossible to evaluate the stromal component of the tumor. Areas to be evaluated were identified under low power (X100), TILs morphology was validated under higher power, and TILs density was assessed as % of covered stromal area under low power if >50% and under higher power (X200) if very low and up to 50%. Efforts were paid to exclude intra-tumoral TILs but this distinction was not always possible due to tumor architecture patterns. Sections from all available blocks per case were evaluated (2 blocks in 43 and 3 blocks in 16 cases; all others, 1 block per case). The entire sections were scanned and average TILs content was recorded per case. Tertiary lymphoid structures, areas with in situ carcinoma or lobular cancerization that often exhibited regional high TILs, and necrotic areas were not considered. Semi-continuous values (5% increments) were recorded per case. Representative examples of the range of TILs densities encountered in the present study are given in Figure 1.

Statistical analysis

Categorical variables were presented as frequencies and percentages while various measures (mean, standard deviation, median, range) were used for continuous variables. Hormonotherapy and trastuzumab were administered based on local breast cancer typing; therefore, local tumor classification based on ER/PgR and HER2 status (IHC/FISH or CISH) was used for analysis in the present study. Concordance for calling HER2 positive tumors and TNBC between local and central pathology was 90.6% (Cohen’s Kappa = 53.1, 95%CI = 47.7-58.5) and 86.9% (Cohen’s Kappa = 64.6, 95%CI = 61.0-68.2), respectively.

Due to the exploratory nature of the study, no correction for multiple testing was applied. In order to assess the reproducibility/validity of the outcome analysis findings, the whole dataset was split into training and validation sets in a controlled fashion. For the latter purpose, nodal status, tumor size, menopausal status, hormone therapy and treatment with trastuzumab were equally assigned into both sets. Patient characteristics did not significantly differ between the two sets (Table S3).

Survival status was updated in June 2014. Disease-free survival (DFS) was set as the primary endpoint of the study and overall survival (OS) as the secondary. DFS was measured from the date of diagnosis until verified disease progression, death or last contact, whichever occurred first, while OS from diagnosis until death from any cause or date of last contact. The annual hazard function for the risk of recurrence (DFS) was estimated by subtypes using penalized B-splines. The number of knots used for the splines was 100. The estimated hazards are presented along with the 95% CIs (Table S4).

There is currently no standard cut-off for TILs classification in BC, while it is suggested that TILs be examined in a continuous mode [17]. Unlike previous reports using 10% increments of recorded TILs values for outcome comparisons (e.g. [7]), here we used three fixed cut-offs: 25%, 50%, and 35%. The 25% cut-off is arbitrary; the 50% cut-off is used for describing lymphocyte-predominant breast cancer (LPBC) but is also considered arbitrary [17]. By using TILs as a continuous variable, the 35% cut-off was obtained by ROC curve analysis in the training set, with DFS at three years as the outcome variable.

The selection of the 3-year DFS as the outcome variable for cut-off assessment was based on the following: a) the follow up period of the last study (HE 10/08) is still short; b) the number of patients at risk at 5 years compared to 3 years were excessively reduced for both HE 10/05 and HE 10/08 studies, leading to overestimation of the risk and probably to biased estimates; c) the annual hazards analysis for the first 5 years showed that for HER2 positive in the pre-trastuzumab era and for TNBC patients the risk of recurrence after 3 years was greatly reduced (Table S4).

Associations among demographic, clinical and treatment characteristics, as well as among TILs cut-offs, were examined. Chi-square tests were used in order to examine possible associations among categorical variables. For testing categorical with continuous variables, the Mann-Whitney or the Kruskal-Wallis test was used, where appropriate. Time-to-event distributions were estimated using the product limit method. Kaplan-Meier curves and log-rank tests were used for comparing time to event distributions and evaluating DFS and OS differences, while univariate Cox analysis was used for and reporting hazard ratios. Univariate Cox with interactions was used in order to identify factors that differentiated TILs’ effect, while Firth correction for monotone likelihood was used for reporting hazard ratios in the case of subgroups with no events.

Survival data including median follow-up for each trial and for the entire population are given in Table S2. Given the large difference in the follow-up periods between the pre- (HE 10/97, HE 10/00) and the post-trastuzumab trials (HE 10/05, HE 10/08), two approaches were used regarding the predictive analysis: a) use of the original data, and b) use of the normalized data. The latter were obtained by truncating the follow-up time of the pre-trastuzumab trials down to the maximum follow-up time of post-trastuzumab trials; this modification was adapted as a further validation step for the results that were obtained with the original data.

Univariate analysis was conducted in the whole dataset, in the training set and in the validation set. Interactions between TILs with the above described BC subgroups and the nodal status were also examined. All univariate tests were two-sided, while significance level was set at α = 0.05. In multivariate analysis the clinicopathological parameters were chosen by backward elimination among the following: age (>50 vs. ≤50), tumor size (>2 cm vs. ≤2 cm), histological grade (I vs. II vs. III-Undifferentiated), positive nodes (≥4 vs. 0-3), treatment group (E-T-CMF vs. ET-CMF vs. E-CMF-Doc vs. E-CMF-T vs. E-CMF), trastuzumab administration (Yes vs. No), adjuvant hormone therapy (Yes vs. No) and subtypes (Luminal A/B vs. Luminal-HER2 vs. HER2-enriched vs. TNBC). Significance threshold for keeping a variable in the final model was set at α = 0.15. Multivariate analysis was conducted in the whole dataset, as well as in each BC subgroup.

The analysis was fully compliant with the reporting recommendations for tumor marker prognostic studies [23]. The SAS software was used for statistical analysis (SAS for Windows, version 9.3, SAS Institute Inc., Cary, NC, USA).


The authors wish to thank Dr. Sofia Chrisafi for tumor block management and Mrs. Maria Moschoni for technical assistance with manuscript processing.


This study was supported by an internal Hellenic Cooperative Oncology Group (HeCOG) translational research grant (HE TRANS_BR). The funders played no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

This study was also partly supported by the Greek General Secretary for Research and Technology (GSRT) Program, Research in Excellence II, funded by 75 % from the European Union and the Operational Program ‘‘Education & Lifelong Learning’’ ESPA-THALIS#266 of the Ministry of Education, Lifelong Learning & Religious Affairs.

conflicts of interest



1. Schreiber RD, Old LJ and Smyth MJ. Cancer immunoediting: integrating immunity’s roles in cancer suppression and promotion. Science. 2011; 331:1565-1570.

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

3. Bianchini G and Gianni L. The immune system and response to HER2-targeted treatment in breast cancer. Lancet Oncol. 2014; 15:e58-68.

4. Denkert C, Loibl S, Noske A, Roller M, Muller BM, Komor M, Budczies J, Darb-Esfahani S, Kronenwett R, Hanusch C, von Torne C, Weichert W, Engels K, Solbach C, Schrader I, Dietel M, et al. Tumor-associated lymphocytes as an independent predictor of response to neoadjuvant chemotherapy in breast cancer. J Clin Oncol. 2010; 28:105-113.

5. Adams S, Gray RJ, Demaria S, Goldstein L, Perez EA, Shulman LN, Martino S, Wang M, Jones VE, Saphner TJ, Wolff AC, Wood WC, Davidson NE, Sledge GW, Sparano JA and Badve SS. 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:2959-66.

6. Denkert C, von Minckwitz G, Brase JC, Sinn BV, Gade S, Kronenwett R, Pfitzner BM, Salat C, Loi S, Schmitt WD, Schem C, Fisch K, Darb-Esfahani S, Mehta K, Sotiriou C, Wienert S, 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-991.

7. Loi S, Michiels S, Salgado R, Sirtaine N, Jose V, Fumagalli D, Kellokumpu-Lehtinen PL, Bono P, Kataja V, Desmedt C, Piccart MJ, Loibl S, Denkert C, Smyth MJ, Joensuu H and Sotiriou C. Tumor infiltrating lymphocytes are prognostic in triple negative breast cancer and predictive for trastuzumab benefit in early breast cancer: results from the FinHER trial. Annals of oncology. 2014; 25:1544-1550.

8. Loi S, Sirtaine N, Piette F, Salgado R, Viale G, Van Eenoo F, Rouas G, Francis P, Crown JP, Hitre E, de Azambuja E, Quinaux E, Di Leo A, Michiels S, Piccart MJ and Sotiriou C. 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-867.

9. Mao Y, Qu Q, Zhang Y, Liu J, Chen X and Shen K. The value of tumor infiltrating lymphocytes (TILs) for predicting response to neoadjuvant chemotherapy in breast cancer: a systematic review and meta-analysis. PLoS One. 2014; 9:e115103.

10. Seo AN, Lee HJ, Kim EJ, Kim HJ, Jang MH, Lee HE, Kim YJ, Kim JH and Park SY. Tumour-infiltrating CD8+ lymphocytes as an independent predictive factor for pathological complete response to primary systemic therapy in breast cancer. Br J Cancer. 2013; 109:2705-2713.

11. West NR, Milne K, Truong PT, Macpherson N, Nelson BH and Watson PH. Tumor-infiltrating lymphocytes predict response to anthracycline-based chemotherapy in estrogen receptor-negative breast cancer. Breast Cancer Res. 2011; 13:R126.

12. Rathore AS, Kumar S, Konwar R, Makker A, Negi MP and Goel MM. CD3+, CD4+ & CD8+ tumour infiltrating lymphocytes (TILs) are predictors of favourable survival outcome in infiltrating ductal carcinoma of breast. Indian J Med Res. 2014; 140:361-369.

13. Dieci MV, Criscitiello C, Goubar A, Viale G, Conte P, Guarneri V, Ficarra G, Mathieu MC, Delaloge S, Curigliano G and Andre F. Prognostic value of tumor-infiltrating lymphocytes on residual disease after primary chemotherapy for triple-negative breast cancer: a retrospective multicenter study. Annals of oncology. 2014; 25:611-618.

14. Perez EA, Thompson EA, Ballman KV, Anderson SK, Asmann YW, Kalari KR, Eckel-Passow JE, Dueck AC, Tenner KS, Jen J, Fan JB, Geiger XJ, McCullough AE, Chen B, Jenkins RB, Sledge GW, 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-708.

15. Salgado R, Denkert C, Demaria S, Sirtaine N, Klauschen F, Pruneri G, Wienert S, Van den Eynden G, Baehne FL, Penault-Llorca F, Perez EA, Thompson EA, Symmans WF, Richardson AL, Brock J, Criscitiello C, et al. Harmonization of the evaluation of tumor infiltrating lymphocytes (TILs) in breast cancer: recommendations by an international TILs-working group 2014. Annals of oncology. 2015; 26:259-71.

16. Nawaz S, Heindl A, Koelble K and Yuan Y. Beyond immune density: critical role of spatial heterogeneity in estrogen receptor-negative breast cancer. Mod Pathol. 2015; 28:766-77.

17. Salgado R, Denkert C, Demaria S, Sirtaine N, Klauschen F, Pruneri G, Wienert S, Van den Eynden G, Baehner FL, Penault-Llorca F, Perez EA, Thompson EA, Symmans WF, Richardson AL, Brock J, Criscitiello C, et al. The evaluation of tumor-infiltrating lymphocytes (TILs) in breast cancer: recommendations by an International TILs Working Group 2014. Annals of oncology. 2015; 26:259-271.

18. Fountzilas G, Skarlos D, Dafni U, Gogas H, Briasoulis E, Pectasides D, Papadimitriou C, Markopoulos C, Polychronis A, Kalofonos HP, Siafaka V, Kosmidis P, Timotheadou E, Tsavdaridis D, Bafaloukos D, Papakostas P, et al. Postoperative dose-dense sequential chemotherapy with epirubicin, followed by CMF with or without paclitaxel, in patients with high-risk operable breast cancer: a randomized phase III study conducted by the Hellenic Cooperative Oncology Group. Annals of oncology. 2005; 16:1762-1771.

19. Gogas H, Dafni U, Karina M, Papadimitriou C, Batistatou A, Bobos M, Kalofonos HP, Eleftheraki AG, Timotheadou E, Bafaloukos D, Christodoulou C, Markopoulos C, Briasoulis E, Papakostas P, Samantas E, Kosmidis P, et al. Postoperative dose-dense sequential versus concomitant administration of epirubicin and paclitaxel in patients with node-positive breast cancer: 5-year results of the Hellenic Cooperative Oncology Group HE 10/00 phase III Trial. Breast Cancer Res Treat. 2012; 132:609-619.

20. Fountzilas G, Dafni U, Papadimitriou C, Timotheadou E, Gogas H, Eleftheraki AG, Xanthakis I, Christodoulou C, Koutras A, Papandreou CN, Papakostas P, Miliaras S, Markopoulos C, Dimitrakakis C, Korantzopoulos P, Karanikiotis C, et al. Dose-dense sequential adjuvant chemotherapy followed, as indicated, by trastuzumab for one year in patients with early breast cancer: first report at 5-year median follow-up of a Hellenic Cooperative Oncology Group randomized phase III trial. BMC Cancer. 2014; 14:515.

21. Fountzilas G, Dafni U, Bobos M, Batistatou A, Kotoula V, Trihia H, Malamou-Mitsi V, Miliaras S, Chrisafi S, Papadopoulos S, Sotiropoulou M, Filippidis T, Gogas H, Koletsa T, Bafaloukos D, Televantou D, et al. Differential response of immunohistochemically defined breast cancer subtypes to anthracycline-based adjuvant chemotherapy with or without paclitaxel. PLoS One. 2012; 7:e37946.

22. Badve S, Dabbs DJ, Schnitt SJ, Baehner FL, Decker T, Eusebi V, Fox SB, Ichihara S, Jacquemier J, Lakhani SR, Palacios J, Rakha EA, Richardson AL, Schmitt FC, Tan PH, Tse GM, et al. Basal-like and triple-negative breast cancers: a critical review with an emphasis on the implications for pathologists and oncologists. Mod Pathol. 2011; 24:157-167.

23. McShane LM and Hayes DF. Publication of tumor marker research results: the necessity for complete and transparent reporting. J Clin Oncol. 2012; 30:4223-4232.

24. Dieci MV, Mathieu MC, Guarneri V, Conte P, Delaloge S, Andre F and Goubar A. Prognostic and predictive value of tumor-infiltrating lymphocytes in two phase III randomized adjuvant breast cancer trials. Annals of oncology. 2015; 26:1698-1704.

25. Loi S. Host antitumor immunity plays a role in the survival of patients with newly diagnosed triple-negative breast cancer. J Clin Oncol. 2014; 32:2935-2937.

26. Schumacher TN and Schreiber RD. Neoantigens in cancer immunotherapy. Science. 2015; 348:69-74.

27. Cancer Genome Atlas N. Comprehensive molecular portraits of human breast tumours. Nature. 2012; 490:61-70.

28. Perez EA, Ballman KV and Anderson SK. Stromal tumor-infiltrating lymphocytes(S-TILs): in the alliance N9831 trial S-TILs are associated with chemotherapy benefit but not associated with trastuzumab benefit. Presented at 2014 San Antonio Breast Cancer Symposium; December 9-13, 2014; San Antonio, Texas. Abstract S1-06. 2014.

29. Zuckerman NS, Yu H, Simons DL, Bhattacharya N, Carcamo-Cavazos V, Yan N, Dirbas FM, Johnson DL, Schwartz EJ and Lee PP. Altered local and systemic immune profiles underlie lymph node metastasis in breast cancer patients. International journal of cancer. 2013; 132:2537-2547.

Creative Commons License All site content, except where otherwise noted, is licensed under a Creative Commons Attribution 3.0 License.
PII: 6231